Page 116 - 2025年第56卷第1期

表 １ 校正后模型参数
模型参数 第 １ 层 第 ２ 层 第 ３ 层河西 第 ３ 层河东
水平渗透系数∕（ｍ∕ｄ） ３５．０ ５．０ ３０．０ １．５
垂直各向异性 ４．０ ３．０ ３．０ ３．０
给水度 ０．０７０ ０．０５５ ０．０４０ ０．０４０
贮水率∕（１∕ｍ） ０．０００１ ０．０００２ ０．０００３ ０．０００３
孔隙度 ０．２５ ０．３２ ０．２５ ０．４０
纵向弥散度∕ｍ ６２．６ ３５．０ ６０．０ ３０．０
横向与纵向弥散度比值 ０．０１ ０．０１ ０．０１ ０．０１
降雨入渗系数 ０．１５
２
河流水力传导率∕（ｍ ∕ｄ） ２００ ２００ ２００ ２００
注： 河西表示研究区黄水河以西部分， 河东表示黄水河以东部分。

图 ５ ＳＥＡＷＡＴ 模型模拟结果与监测结果对比

模型量化不同抽咸注淡策略对地下坝内陆侧咸水回退的影响， 进而提出应用混合整数多目标优化模型
来优化井群抽注的时空分配。 为尽可能降低地下水抽注水量， 同时最大程度减缓内陆侧海水入侵程
度， 以抽（注）水井的抽（注）水量和抽（注）水启闭状态作为决策变量， 以地下坝内陆侧区域的盐分质
量作为状态变量。 通过计算地下坝内陆侧区域的盐分质量变化来衡量海水入侵的平均程度， 优化模型
具体形式如下。
目标函数 １： 使咸水的抽取量和淡水的注入量总和达到最小， 以减少能源消耗节约成本

ｍｉｎ Ｆ ＝ ｔ ｉ ｔ ｉ ｔ ｔ （９）
１ ∑∑ （ｚ Ｑ ＋ｚ Ｉ ）
ｊ ｊ
ｉ， ｊ ｔ
目标函数 ２： 使地下坝内陆侧区域盐分质量减少程度达到最大， 从而降低海水入侵平均程度
ｍａｘＦ ＝ ［（ｍ －ｍ ）∕ｍ ］×１００ （１０）
２ ｉｎｉ ｅｎｄ ｉｎｉ
约束条件 １： 模拟模型得到的地下坝内陆侧盐分质量的等式约束
ｔ ｔ ｔ ｔ
ｍ ＝ ｆ（Ｑ ， Ｉ ， ｚ ， ｚ ） （１１）
ｉ
ｊ
ｉ
ｊ
约束条件 ２： 表征区内所有注水量之和不能超过一个最大供水量
ｔ ｔ
∑ ｚ Ｉ ≤Ｉ ｍａｘ （１２）
ｊ ｊ
ｊ， ｔ
ｔ
ｔ
式中： Ｑ 为第 ｉ 口抽水井在 ｔ 时间内的抽水量； Ｉ 为第 ｊ 口注水井在 ｔ 时间内的注水量； ｚ 和 ｚ 为二进
ｔ
ｔ
ｉ
ｊ
ｊ
ｉ
ｔ ｔ
制变量， 取 ０ 或 １， 用来控制第 ｉ 口抽水井或第 ｊ 口注水井在时间 ｔ 内的启闭状态， 如果 ｚ ＝ １ 或 ｚ ＝ １，
ｉ ｊ
表示第 ｉ 口抽水井或第 ｊ 口注水井在时间 ｔ 内处于运行状态， 反之， 如果取值为 ０， 则代表抽（注）水井
在时间 ｔ 内处于关闭状态； ｍ 为修复管理期初始时刻地下坝内陆侧的盐分质量； ｍ 为修复管理期末
ｉｎｉ ｅｎｄ
— １ １ １ —