Page 76 - 2022年第53卷第11期

水位变化进行匀速处理。模型采用循环迭代法计算，输入初始参数计算后得到新的参数，作为下一次
运算的初始参数，如此循环往复直至运算结束。模型借助 Ｐｙｔｈｏｎ软件编程运算。

３ 结果分析

３．１ 浅层地下水超采量 依据研究整理的 １９８０、１９８５、１９９０、１９９５、１９９８—２０１９年浅层地下水开采
量数据，绘制浅层地下水开采量变化图（见图 ５（ａ）），缺失年份参考降雨量及相邻年份开采量数据进
３
行插值。１９８０—２０１９年，海河平原区浅层地下水累计开采量为 ５２１５亿ｍ ，开采量整体呈现先增加后
３
３
减少趋势，２０１９年开采量为 １３２．９亿ｍ ，较年开采量最多的年份（２００１年）减少了 ３５．９亿ｍ 。将单位
面积开采量大小定义为开采强度，１９８０—２０１９年海河平原区各三级区累计开采量和开采强度分布如
图 ５（ｂ）所示（图中：①滦河及冀东沿海诸河，②北四河下游平原，③大清河淀西平原，④大清河淀东
平原，⑤子牙河平原，⑥漳卫河平原，⑦黑龙港及运东平原，⑧徒骇马颊河平原，下同），大清河淀
２
３
西平原累计开采量和累计开采强度均最大，累计开采强度为８．５ｍ ?ｍ ；漳卫河平原虽然累计开采量较
３
２
小，但开采强度仅低于大清河淀西平原，为 ６．６ｍ ?ｍ 。

图 ５ １９８０—２０１９年海河平原区浅层地下水开采量

１９５９—２０１９年的浅层地下水位变差图如图 ６（ａ）所示，山前平原洪冲积扇缘一带浅层地下水位下
降最为严重，形成 “高蠡清－ 肃宁” “宁柏隆漏斗” “雄县霸州” 等一连串浅层漏斗，地下水位降幅从
山前平原向东部滨海平原逐渐减小，整个平原区地下水位平均下降约 １１．５ｍ。

图 ６ 研究区 １９５９—２０１９年浅层地下水水位差及累计超采强度
４
— １ ３ ２ —