The assessment and analysis of ecosystem service functions are of great significance for the rational allocation and optimization of regional resources. Taking Cangzhou City as the research object, the ecosystem service functions in 2005, 2015, and 2020 were evaluated based on the InVEST model. The Spearman correlation coefficient was used to analyze the trade-off and synergy between ecosystem services, and the self-organizing mapping method was used to identify ecosystem service bundles. The research results indicate that the water yield in Cangzhou City shows a trend of first increasing and then decreasing, while the N/P output first decreases and then increases, while carbon storage and habitat quality remain stable with some improvement. There is a significant synergistic relationship between water yield and N output, water yield and P output, N output and P output, carbon storage and habitat quality. There is a significant trade-off between water yield and carbon storage, water yield and habitat quality, N output and carbon storage, N output and habitat quality, P output and carbon storage, and P output and habitat quality. There are three ecosystem service bundles in Cangzhou City, namely ESB1 focused on water yield and N/P output, ESB2 focused on habitat quality and carbon storage, and ESB3 focused on cultural tourism.The research results provide decision-making basis for improving the overall efficiency of ecosystem services and differentiated management of ecological functional areas in Cangzhou City.