第３１卷第２期２ ０ １
３年２月水 电 能 源 科 学
Ｗａｔｅｒ Ｒｅｓｏｕｒｃｅｓ ａｎｄ ＰｏｗｅｒＶｏｌ．３１Ｎｏ．２
Ｆｅｂ．２ ０ １
３文章编号：１０００－７７０９（２０１３）０２－００２３－
华舒愉，顾圣平，贺 军，张 楠，郑贵桥
（河海大学水利水电学院，江苏南京２１００９８
）摘要：三水源新安江模型含有比其他模型更多的参数，参数间的相互作用使其成为影响模型模拟结果的重要因素。为使人工模拟结果与实际水文测量数据及水文现象较为一致，基于三水源新安江模型结构及参数敏感性建立参数优化率定的目标函数，并应用于竹竿河流域，取得了较好的率定及校验结果，证明该模型及优化方法有效可行。
关键词：水文模型；新安江模型；参数优化；应用中图分类号：ＴＶ１２４
文献标志码：Ａ
收稿日期：２０１２－０９－２７，修回日期：２０１２－１１－
１４基金项目：水利部公益行业科研专项经费基金资助项目（２０１１０１０１１
）作者简介：华舒愉（１９８７－），女，硕士研究生，研究方向为水利水电工程，Ｅ－ｍａｉｌ：４７６４７０２４３＠ｑｑ．ｃｏｍ通讯作者：顾圣平（１９５７－），男，教授，研究方向为水资源系统规划，Ｅ－ｍａｉｌ：ｓｐｇ
ｕ＠ｈｈｕ．ｅｄｕ．ｃｎ ２
０世纪７０～８０年代，赵人俊［１］
借鉴山坡水文学和产汇流理论提出了三水源新安江模型是概念性分散性模型，多适用于我国南方湿润与半湿润地区。使用该模型时需根据具体情况将全流域分为多个单元流域，按模型结构对每个单元流域进行蒸散发、产流、水源划分和汇流计算，得到单元流域出口的流量过程，再对其出口以下的河道进行汇流计算，将每个单元流域在全流域出口的流量过程线性叠加，即得全流域出口的总出流过程。新安江模型在湿润半湿润地区使用至今，其自身结构已基本得到肯定。但在不同模拟流域的产汇流参数不尽相同时，参数的影响作用都集中反映在模拟结果中，所以必须通过参数优化率定为模拟及产水预报计算提供适宜的参数，本文就此研究了新安江模型的参数优化问题，取得了较
满意的结果。
! 模型结构及参数率定
!$!新安江模型的分层结构
图１为三水源新安江水文模型结构图，包含有蒸散发、产流、分水源、汇流计算四个层次。四个层次均可通过模型中的公式进行计算，上一层次的计算结果可作为下一层次的数据输入，而下一层次的参数变化对上一层次的结果不产生影响。这种各层次间的相对独立性为模型中的参数
图１ 三水源新安江模型结构Ｆｉｇ．１ Ｓｔｒｕｃｔｕｒｅ ｄｉａｇ
ｒａｍ ｏｆ Ｔｈｒｅｅ－ｗａｔｅｒＳｏｕｒｃｅｓ Ｘｉｎ’ａｎｊｉａｎｇ
Ｍｏｄｅｌ确定提供了便利［
１］
。!$" 新安江模型中的参数及其优化率定
新安江模型四个层次的计算中一共涉及到
１６个参数［２，３］
，
第一层次蒸散发计算中的参数有流域蒸散发折算系数Ｋ、上层张力水容量ＷＵＭ、
下层张力水容量ＷＬＭ、
深层蒸散发折算系数Ｃ；第二层次产流计算中的参数有流域平均张力水容量ＷＭ、
张力水蓄水容量曲线方次Ｂ、不透水面积占全流域面积的比例ＩＭ；第三层次水源划分中的参数有自由水蓄水容量ＳＭ、
自由水蓄水容量曲线方次ＥＸ、壤中流出流系数ＫＳＳ、地下径流出流系数ＫＧ；第四层次汇流计算中的参数有壤中流消退系数ＫＫＳＳ、地下径流消退系数ＫＫＧ、地面径流单位线ＵＨ、马斯京根法中的参数ＫＥ、ＸＥ。由于模型中的参数系统反映了水文过程的物理规律，在取值范围内的增大或减小会对模拟结果造成一定的影响。经参数敏感性分析可知，模型对其中部分
参数的响应结果较不敏感，参数变化对输出影响不大，如ＷＭ、ＷＵＭ、ＷＬＭ、Ｂ、Ｃ、ＩＭ、ＥＸ、ＫＥ、ＸＥ敏
感度较低，
按经验值即可确定，仅在必要时做单独调整［４］，单位线是根据已知雨量及流量资料推求
出的适用于该流域的时段单位线，即从中选取几段时空分布较均匀，
历时较短的降雨形成的单峰洪水来推求适用于该单元流域及类似流域的无因
次单位线ＵＨ［５］
。为提高模拟的精确度，可按径流
量级划分确定单位线；其余参数均为敏感性较大的参数，沈圣哲
包括影响产流结果的Ｋ及影响径流变化过程的ＳＭ、ＫＧ、ＫＳＳ、ＫＫＳＳ、ＫＫＧ。针对以上敏感性较大的参数建立目标函数以便进行优化率定。（１
）在蒸散发和产流计算层次中较为敏感的参数是蒸散发折算系数Ｋ，输入值为各时段降雨量Ｐ，输出值为各时段产流量Ｒ，Ｒ的计算结果受
Ｋ的影响。由于各时段实际与计算值的误差可
正可负，若利用相对误差极小化确定参数Ｋ，简单求和可能将很大的误差抵消掉，所以基于最小二乘法原理（ＯＬＳ）建立目标函数能更好地反映在总体上的接近程度，表达式为：
Ｆ１＝ｍｉｎ∑ｎ
ｉ＝１
ＳＣＪＬ（ｉ）／Ｕ－Ｒ（ｉ［
］）２
ｉ＝１，２，…，ｎ（１
）
其中
Ｕ＝Ｆ／（３．６Δ
ｔ）式中，ＳＣＪＬ（ｉ）为实测径流量；Ｕ为单位换算系数；Ｆ为面积；Δｔ为计算时段长；Ｒ（ｉ）为产流量；ｉ为计算时段；ｎ为总时段数。
（２
）在进行水源划分和汇流计算时，这两部分的参数是独立的。由于地面径流与流量过程的高水部分相关联，参数ＳＭ、ＫＳＳ、ＫＧ影响地面径流值；地下径流与低水部分相关联，参数ＫＫＳＳ、ＫＫＧ影响地下径流值，二者共同影响ＴＲ（ｉ）及径流过程。通常情况下将高低水部分分开处理，可得到适宜的参数。
当采用以日为模拟时段时，为削减高水部分引起的误差，常选用相关的对数函数作为目标函数。当自变量值较大时，对数函数可将误差放大化，
利用这一特性使实测流量与计算流量之比的对数绝对值与实测流量对数值的比值极小化，即可优选出参数ＳＭ、ＫＳＳ、
ＫＧ，目标函数为：Ｆ２＝
ｍｉｎ∑ｎ
ｉ＝１
ａ
ｂｓ Ｌｏｇ ＳＣＪＬ（ｉ）／ＴＲ（ｉ（）［］）∑ｎｉ＝１
Ｌ
ｏｇ ＳＣＪＬ（ｉ（））（２
）式中，ＴＲ（ｉ
）为计算径流量。考虑到参数ＫＫＳＳ、
ＫＫＧ对流量过程线低水部分的影响，类似目标函数式（２），将实测流量与计算流量的误差绝对值与实测流量之比极小化作为约束，建立的目标函数为：
Ｆ３＝ｍｉｎ∑ｎ
ｉ＝１
ａ
ｂｓ ＳＣＪＬ（ｉ）－ＴＲ（ｉ［］）∑ｎ
ｉ＝１
ＳＣＪＬ（ｉ）（３
）（３
）为保证模拟过程与实测过程之间总体的吻合程度，参考《水文情报预报规范》［６］
精度评定中的确定性系数ＤＣ建立目标函数。ＤＣ在０～１之间变化，其值越接近１，说明两者的吻合程度越高，模型拟合效果越好。据此建立的目标函数为：
Ｆ４＝ｍｉｎ（１－ＤＣ）
（４
）其中
ＤＣ＝１－
∑ｎ
ｉ＝１
ＴＲ
（ｉ）－ＳＣＪＬ（ｉ［
］
）２
∑ｎ
ｉ＝１
ＴＲ
（ｉ）－Ｐ［
］
ＪＪＬ２
式中，ＰＪＪＬ为实测平均径流量。
当洪水预报精度达到甲（ＤＣ≥０．９０）、乙（０．９０＞ＤＣ≥０．７０
）两个等级才可用于正式预报，所以针对以上目标函数的约束条件为０．７０≤ＤＣ≤１。!$# 参数优化率定技术路线
以上目标函数建立起来的多目标极小化模型
（ＶＭＰ）
，可根据“理想点法”进行求解，即分别求出各目标函数的极小值Ｆｍｉｎ１、Ｆｍｉｎ２、Ｆｍｉｎ３、Ｆｍｉｎ
４，通常情况下每个目标函数对应的最优解不完全相同，
各个最小值Ｆｍｉｎ
ｍ（ｍ＝１，２，３，４）分别是对应目标函数最理想的值，点Ｇｍｉｎ＝（Ｆｍｉｎ１，Ｆｍｉｎ２，
Ｆｍｉｎ３，Ｆｍｉｎ４）Ｔ
即称为ＶＭＰ模型的理想点［
５］
。针对该问题采用“最短距离法”进行求解，令Ｘ＝（Ｋ，ＳＭ，ＫＧ，ＫＳＳ，
ＫＫＳＳ，
ＫＫＧ）Ｔ
，原多目标优化问题可写成Ｇ（Ｘ）＝Ｆ１（Ｘ），Ｆ２（Ｘ），Ｆ３（Ｘ），Ｆ４（Ｘ［］）Ｔ
的形式。分别求解Ｆｍ（Ｘ）（ｍ＝１，２，３，４）得到极小值Ｆｍｉｎ
ｍ（其中Ｘｍ为极小值Ｆｍｉｎ
ｍ
情况下的最优解，继而建立距离模评价函数为：
φ（）［
］Ｇ Ｘ＝‖（）Ｇ Ｘ－Ｇｍｉｎ
‖＝∑４
ｉ＝１
Ｆｉ
（）Ｘ－Ｆｍｉｎ［
］
ｉ
槡
２
（５
）
Ｘ∈Ｒ，Ｒ＝｛Ｘ｜Ｋｍｉｎ＜Ｋ＜Ｋｍａｘ，ＳＭｍｉｎ
＜ＳＭ＜ＳＭｍａｘ，ＫＳＳｍｉｎ＜ＫＳＳ＜ＫＳＳｍａｘ，ＫＫＳＳｍｉｎ
＜ＫＫＳＳ＜ＫＫＳＳｍａｘ，ＫＧｍｉｎ＜ＫＧ＜ＫＧｍａｘ，ＫＫＧｍｉｎ＜ＫＫＧ＜ＫＫＧｍａｘ
，０．７０≤ＤＣ≤１｝式（５）是取目标Ｇ与理想点Ｇｍｉｎ
之间的距离作为评价函数，即将原问题转换为求解φ（）［］
Ｇ Ｘ的极小值ｍｉｎ‖（）Ｇ Ｘ－Ｇｍｉｎ
‖问题。由于变量Ｘ的值不固定，故利用策略空间迭代法给变量赋予初值，通常情况下将各变量的下界值（即Ｋｍｉｎ、
·
４２·水 电 能 源 科 学 ２０１３年
第３１卷第２期华舒愉等：三水源新安江模型参数优化及其应用
ＳＭｍｉｎ、ＫＳＳｍｉｎ、ＫＧｍｉｎ、ＫＫＳＳｍｉｎ、ＫＫＧｍｉｎ
）作为初始值，设置合理的迭代步长即可在约束范围内得到不同的变量值。在约束范围内求各目标函数的极小值，再利用式（５）求得最优解必为原始多目标问题的有效解。
依据具体的参数优化目标函数及模型结构，使用ＶＢ语言编写模型模块及参数优化模块，最终可得到适合不同流域的最优解，使模拟结果与实际情况较为吻合，
为之后的水资源优化配置等综合利用规划提供有效的数据。
" 实例
竹竿河位于淮河南罗山县，为一条低山丘陵
区的河流，流域面积约１
６３９ｋｍ２
，属于半湿润地区。根据１９９７～２００４年的日雨量及流量资料，考虑资料的完整性及其他因素，从中选取２００２年为参数率定期，２００３年为校验期。
由于该流域面积较小且仅有一个控制站点，不需要再分成单元流域进行计算，可忽略河道演算部分，不考虑参数ＫＥ、ＸＥ。另外，在已知流量过程的情况下可利用多段完整的洪水流量过程图，推求该流域的退水曲线，继而利用经验公式ＫＫＧ＝－１／ｌｎＣｇ（其中Ｃｇ＝Ｑ
ｔ＋１／Ｑｔ）计算出该流域的地下径流消退系数ＫＫＧ＝０．９８５，该参数便可不参与率定，
减少了参数优化的数量和难度。根据该地区水文规律，结合新安江模型参数取值经验，
给各个参数赋予初始值，见表１。表１ 参数初始值
Ｔａｂ．１ Ｖａｌｕｅ ｏｆ ｉｎｉｔｉａｌ ｐ
ａｒａｍｅｔｅｒｓ
对蒸散发折算系数Ｋ，
以０．０１的步长进行迭代，
将部分Ｋ值相应得到的Ｆ１值列入表２中。当Ｋ取值为０．５５时，Ｆ１为最
小值，此时计算与实测流量的相对误差为５％。
表２ 不同Ｋ值情况下的Ｆ１值Ｔａｂ．２ Ｆ１ｖ
ａｌｕｅ ｕｎｄｅｒ ｄｉｆｆｅｒｅｎｔ Ｋｖａｌｕｅ ＫＧ和ＫＳＳ是
代表自由水出流快慢的两个参
数，一般情况下１ ０００ｋｍ２
左右的流域从雨止到壤中流止的时间为３～４ｄ，可令［１－（ＫＳＳ＋
ＫＧ）］ｔ＝０．０１，其中的ｔ为４ｄ，则有ＫＳＳ＋ＫＧ＝
０．
７。因此，令参数ＳＭ和ＫＧ在不同的取值范围内变化迭代，ＫＳＳ也随着ＫＧ而变化，得其对应目标函数Ｆ２值及其变化趋势图，
见图２。图２ 不同ＳＭ、ＫＧ值对应的Ｆ２变化趋势图
Ｆｉｇ．２ Ｃｈａｒｔ ｏｆ Ｆ２ｃｈａｎｇ
ｅ ｔｒｅｎｄ ｗｉｔｈ ｄｉｆｆｅｒｅｎｔ ＳＭａｎｄ ＫＧ令ＳＭ在２０～５０ｍｍ范围内以５ｍｍ的步长
迭代，ＫＧ在０．２～０．６范围内以０．０５的步长迭代，由图
２观察可得Ｆ２的变化规律，当ＳＭ＝２５ｍｍ、ＫＧ＝０．２０时，ＫＫＳ＝０．７０－ＫＧ＝０．
５０，Ｆ２＝
０．１７５，为Ｆ２的极小值。
壤中流消退系数ＫＫＳＳ通过不断迭代得到目标函数Ｆ３值，由图３可直观地发现其基本呈抛物线形的变化规律，当ＫＫＳＳ＝０．７０时，Ｆ３＝０．５１３，该值为Ｆ３的极小值。
图３ 不同ＫＫＳＳ值对应的Ｆ３变化趋势图Ｆｉｇ．３ Ｃｈａｒｔ ｏｆ Ｆ３ｃｈａｎｇ
ｅ ｔｒｅｎｄ ｗｉｔｈ ｄｉｆｆｅｒｅｎｔ ＫＫＳＳ在以上参数不断迭代的同时，通过模型计算亦可得到相应的确定性系数ＤＣ，当Ｋ、ＫＧ、ＳＭ、ＫＫＳＳ取值分别为０
．６５、０．４０、３５．００、０．８０时，确定性系数为０．９８，此时目标函数Ｆ４达到最小值０
．０２。每个目标函数相应的极小值一经确定，便可利用最短距离法建立的距离模评价函数，求得目标函数与理想点之间的最小距离ｍｉｎ‖（
）Ｇ Ｘ－Ｇｍｉｎ
‖及该最小距离情况下的各个
变量值，即为优化得到的参数值。以上过程均通过程序自动优化计算，最后确定优化率定得到的参数解见表３。
表３ 参数优化率定表
Ｔａｂ．３ Ｐａｒａｍｅｔｅｒ ｌｉｓｔ ａｆｔｅｒ ｏｐ
ｔｉｍｉｚａｔｉｏｎ 通过优化率定及计算的２
００２年计算径流与·
５２·
实测径流对比结果见图４（ａ）及表４。分析表明，２００２年计算与实测流量过程线较为吻合，年径
流相对误差为１．０８％，总误差较小，模拟得到确定性系数ＤＣ＝０．９０，精度为甲级，二者总体趋势较为一致，满足正式预报的要求。此外参考《水文情报预报规范》［６］，以实测值的２０％作为许可精度，月径流合格率（即月径流量相对误差小于２０％的月份占总月份的百分比）为７５％，精度达到乙级，可用于正式预报。根据参数优化率定结果，对校验期２００３年进行模拟，检验其合理性，得到２００３年的对比结果见表４及图４（ｂ）。模拟得到确定性系数ＤＣ＝０．９８，二者总体趋势较为吻合，满足
正
图４ 计算流量与实测流量过程对比
Ｆｉｇ．４ Ｃｏｍｐａｒｉｓｏｎ ｏｆ ｃａｌｃｕｌａｔｅｄ ａｎｄ ｍｅａｓｕｒｅｄ ｆｌｏｗｒａｔｅ式预报的要求；计算与实测年径流系数分别为０．５７７、０．６０７，相对误差为４．９３％，月径流合格率为７５％，精度达到乙级，可用于正式预报。
同时，针对该区域资料齐全的２０００、２００４、２００５年分别进行模拟校验，可得计算与实测总年径流量相对误差分别为８．１％、１２．８％、１１．５％，确定性系数分别为０．７７、０．７１、０．８６，均满足《水文情报预报规范》［６］中的精度要求，因此以上参数可用于正式预报。
# 结语
ａ．通过对模型结构及其参数敏感度分析，模拟结果表明三水源新安江模型在竹竿河流域具备较高的适用性。
ｂ．以多目标函数优选敏感度较高的参数方法可行，预报精度基本满足要求，利用计算机技术优选率定参数简单、易操作、计算效率高。
ｃ．原始资料的准确性、模拟过程中的计算方法及计算数值的误差、暴雨中心的不同、产流量级的分类都是导致计算与实测流量误差的原因，实际应用中需尽量减小误差，使模拟结果更接近实际径流。
ｄ．在确定参数最优解时会出现很多组不同的参数值使模型输出具有相同的拟合精度，这种“异参同效”的现象是水文过程的复杂性、模型参数之间的相互补偿作用及随机性所引起的。为避免该现象，可根据需求采用人为的主观设置目标函数的权重、迭代循环次数等方法，同时必须加强对水文过程物理规律及参数率定技术的研究。
表４ 计算与实测月（累加）径流量对比
Ｔａｂ．４ Ｃｏｍｐａｒｉｓｏｎ ｏｆ ｃａｌｃｕｌａｔｅｄ ａｎｄ ｍｅａｓｕｒｅｄ ｍｏｎｔｈｌｙ（ｃｕｍｕｌａｔｉｖｅ）ｆｌｏｗｒａｔｅ
月份
２００２年
累加降
雨量／ｍｍ
总径流量／（ｍ３·ｓ－１）
计算实测
相对误
差／％
２００３年
累加降
雨量／ｍｍ
总径流量／（ｍ３·ｓ－１）
计算实测
相对误
差／％
１ １３．７ １２２．９２８ １２８．４７－４．３１ １９．０ １６７．１７ ２ ０５．０９－１８．４９
２ ３５．５ １４０．０１０ １１７．９３ １８．７３ ８２．０ ４７９．９３ ４ ５８．９１ ４．５８
３ ７９．４ ５９５．３５９ ７３８．４７－１９．３８ １２０．０ １ ３４０．１８ １ ６４９．８０－１８．７７
４ ６８．０ ３２３．７０４ ２４５．１８ ３２．０３ ８８．８ ６９０．５６ ８２６．２０－１６．４２
５ １７６．３ １ ４４０．９９４ １ ７２４．００－１６．４２ ９０．９ １ １３３．９９ １ ５５４．０８－２７．０３
６ ３６４．６ ４ ４８１．５９５ ３ ９７３．９４ １２．７７ ２６２．７ １ １４０．６２ １ ７２９．１５－３４．０４
７ ２５２．５ ３ ３４６．９１０ ３ ８２８．３２－１２．５８ ２４７．５ ４ ６１８．３５ ４ ６２０．１０－０．０４
８ ５２．７ ８７１．８０８ ７３３．７０ １８．８２ １１９．４ １ ２４８．３６ ８９６．７３ ３９．２１
９ ２５．６ ４５８．３３３ ２９７．４６ ５４．０８ ６９．２ ９６５．３７ ９１１．３０ ５．９３
１０ ３６．６ ３７４．０７９ ２１２．０２ ７６．４４ １０２．７ １ ２７６．２７ １ ０７０．５９ １９．２１
１１ ９．７ １４２．２６２ １２６．３９ １２．５６ ６８．３ ７１０．００ ６０２．６０ １７．９３
１２ ７７．６ ３８３．７１５ ４１９．７５－８．５８ ２３．５ ３９１．２５ ３７２．２６ ５．１０
累计１ １９２．２ １２ ６８１．７１０ １２ ５４５．６３ １．０８ １ ２９４．０ １４ １６２．０７ １４ ８９６．８０－４．９３
（下转第２４２页）
（上接第２６页）
参考文献：
［１］ 赵人俊．流域水文模型———新安江模型与陕北模型［Ｍ］．北京：水利电力出版社，１９８４．
［２］ 王静，黄强，王义民．基于启发式算法新安江模型参数优选研究进展［Ｊ］．水电能源科学，２０１１，２９（４）：２８－３０．［３］ 赵人俊，王佩兰．新安江模型参数的分析［Ｊ］．水文，１９８８（６）：２－９．［４］ 孔凡哲．新安江模型（三水源）参数的调试［Ｊ］．水资源研究，２００２，２３（１）：２９－３１．
［５］ 周凤歧，崔敬波．单位线法推求流域出口洪水过程综述［Ｊ］．吉林水利，２００５（８）：５－７．
［６］ 水利部水利信息中心，水利部长江水利委员会水文局，水利部黄河水利委员会水文局，等．水文情报预
报规范（ＳＬ２５０－２０００）［Ｓ］．北京：中国水利水电出版
社，２０００．
Ｐａｒａｍｅｔｅｒｓ Ｏｐｔｉｍｉｚａｔｉｏｎ ａｎｄ Ａｐｐｌｉｃａｔｉｏｎ ｏｆ Ｔｈｒｅｅ－ｗａｔｅｒ
Ｓｏｕｒｃｅｓ Ｘｉｎ’ａｎｊｉａｎｇ Ｍｏｄｅｌ
ＨＵＡ Ｓｈｕｙｕ，ＧＵ Ｓｈｅｎｇｐｉｎｇ，ＨＥ Ｊｕｎ，ＺＨＡＮＧ Ｎａｎ，ＺＨＥＮＧ Ｇｕｉｑｉａｏ
（Ｃｏｌｌｅｇｅ ｏｆ Ｗａｔｅｒ Ｃｏｎｓｅｒｖａｎｃｙ ａｎｄ Ｈｙｄｒｏｐｏｗｅｒ，Ｈｏｈａｉ Ｕｎｉｖｅｒｓｉｔｙ，Ｎａｎｊｉｎｇ ２１００９８，Ｃｈｉｎａ）Ａｂｓｔｒａｃｔ：Ｔｈｒｅｅ ｗａｔｅｒ Ｓｏｕｒｃｅｓ Ｘｉｎ＇ａｎ ｊｉａｎｇ Ｍｏｄｅｌ ｈａｓ ｍｏｒｅ ｐａｒａｍｅｔｅｒｓ ｔｈａｎ ｏｔｈｅｒ ｍｏｄｅｌｓ．Ｔｈｅ ｉｎｔｅｒａｃｔｉｏｎ ｏｆ ｐａｒａｍ－ｅｔｅｒｓ ｍａｋｅｓ ｔｈｅｍ ｂｅｃｏｍｅ ｉｍｐｏｒｔａｎｔ ｆａｃｔｏｒｓ ｔｏ ａｆｆｅｃｔ ｔｈｅ ｓｉｍｕｌａｔｉｏｎ ｒｅｓｕｌｔｓ ｉｎ ｔｈｅ ｍｏｄｅｌ．Ｉｎ ｏｒｄｅｒ ｔｏ ｍａｋｅ ａｒｔｉｆｉｃｉａｌ ｓｉｍｕｌａ－ｔｉｏｎ ｒｅｓｕｌｔｓ ｃｏｎｓｉｓｔｅｎｔ ｗｉｔｈ ｐｒａｃｔｉｃａｌ ｈｙｄｒｏｇｒａｐｈｉｃ ｓｕｒｖｅｙ ｄａｔａ ａｎｄ ｈｙｄｒｏｌｏｇｉｃａｌ ｐｈｅｎｏｍｅｎａ，ｔｈｉｓ ｐａｐｅｒ ｅｓｔａｂｌｉｓｈｅｓ ｏｂｊｅｃ－ｔｉｖｅ ｆｕｎｃｔｉｏｎ
ａｂｏｕｔ ｔｈｅ ｐａｒａｍｅｔｅｒ ｏｐｔｉｍｉｚａｔｉｏｎ ｂａｓｅｄ ｏｎ ｔｈｒｅｅ－ｗａｔｅｒ ｓｏｕｒｃｅｓ Ｘｉｎ＇ａｎｊｉａｎｇ ｍｏｄｅｌ ｓｔｒｕｃｔｕｒｅ ａｎｄ ｐａｒａｍｅｔｅｒｓｓｅｎｓｉｔｉｖｉｔｙ．Ａｎｄ ｔｈｅｎ ｔｈｅ ｍｏｄｅｌ ｉｓ ａｐｐｌｉｅｄ ｔｏ Ｚｈｕｇａｎｈｅ Ｂａｓｉｎ．Ｉｔ ｇｅｔｓ ｔｈｅ ｂｅｔｔｅｒ ｃａｌｉｂｒａｔｉｏｎ ａｎｄ ｖｅｒｉｆｉｃａｔｉｏｎ ｒｅｓｕｌｔｓ，ｗｈｉｃｈｐｒｏｖｅｓ ｔｈｅ ｆｅａｓｉｂｉｌｉｔｙ ｏｆ ｔｈｅ ｍｏｄｅｌ ａｎｄ ｏｐｔｉｍｉｚａｔｉｏｎ ｍｅｔｈｏｄ．
Ｋｅｙ ｗｏｒｄｓ：ｈｙｄｒｏｌｏｇｉｃａｌ ｍｏｄｅｌ；Ｘｉｎ＇ａｎｊｉａｎｇ ｍｏｄｅｌ；ｐａｒａｍｅｔｅｒ ｏｐｔｉｍｉｚａｔｉｏｎ；ａｐｐｌｉｃａｔｉｏｎ
发布评论