Page 86 - 2022年第53卷第4期

式中：ET can 为渠系蒸发损失量，mm；Irr can 为进入渠道的灌溉水量，mm；Ls can 为渠系渗漏损失量，
mm；Surp can 为渠系退水量，mm；φ为渠系水有效利用系数；α为渠系蒸发系数；β为渠系渗漏系数。

（3）水源模块的衔接。多源供水模块首先识别计算单元的土地利用类型。若是居工地，则程序进
入城乡供水计算流程；若是农田用地，则程序进入灌溉供水计算流程；若是其他土地利用类型，则
该计算单元不进行供用水模拟。然后，根据设定的取水优先序规则识别计算单元的水源识别码和取
水优先序代码，依次调用各水源模块（图 4）。河道、水库、浅层地下水、深层地下水、外调水、坑塘
水和再生水的水源识别码分别是 1、2、3、4、5、6、7。

农田用地？ 计算单元 居工地？
是 是

灌溉水源代码及供水序识别 城乡水源代码及供水序识别

否 否 否 否 否 否
wsc=1 wsc=2 wsc=3 wsc=4 wsc=5 wsc=6 wsc=7
是 是 是 是 是 是 是
河道水 水库水 浅层地下水 深层地下水 外调水 坑塘水 再生水
调用 调用
调用
rchuse 模块 res 模块 watuse 模块

否 否 否 否 否
irrsc=1 irrsc=2 irrsc=3 irrsc=4 irrsc=5 irrsc=6
是 是 是 是 是 是
河道水 水库水 浅层地下水 深层地下水 外调水 坑塘水
调用 调用
调用
irr_rch 模块 irr_res 模块 irrsub 模块

图 4 多源供水模块中的水源模块调用流程

（4）地下水模块修改。对 gwmod 模块相关代码进行修改，描述渠系渗漏水量、城镇供水管网漏失
水量补给潜水这一过程。修改公式如下：
)
)
é
ù
rh t = 1 - exp(-1 delay · ( prc + Ls can + Ls pip + rh t - 1 ·exp(-1 delay ) （4）
û
ë
Ls pip = WSP· (1 - pip ) Area （5）
式中：rh t 为第 t 天的潜水补给量，mm；rh 为第 t-1 天的潜水补给量，mm；prc 为补给地下水的土壤
t-1
渗漏水量，mm；Ls can 为渠系漏损量，mm；Ls 为城镇供水管网漏失量，mm；delay 为地下水补给延
pip
2
迟系数，mm；WSP 为水源取水量，m ；pip 为城镇供水管网漏损率；Area 为计算单元的面积，m 。
3
（5）地表产流模块修改。对 surface 模块进行修改，在地表产流计算公式中添加渠道退水参数。这
样，渠道退水随地表径流进入河道，参与后续的河道汇流过程模拟。
（6）点源、面源模块衔接。为实现水资源配置模型的排水模块与 SWAT 模型的点源模块（recday 模
块和 recmon 模块）的衔接，将点源模块中的排水参数置换为排水模块计算出的排水参数，采用错日传
递法 ［14］ 将排水数据传递给点源模块，将产生的污水安排在下一日模拟时排放，进而实现点源的逐日
衔接模拟。

— 460 —