技术特征：
1.一种基于改进新安江模型的实际径流量预报方法，包括如下步骤：步骤1：基于新安江模型研发半分布式天然径流量预报模块，对天然径流量进行预报；步骤2：研发用水总量统计模块，对农业灌溉用水量和非农业灌溉用水量分别进行统计；步骤3：研发基于改进新安江模型的流域水循环模型，计算实际径流量；步骤4：收集基础资料，对水循环模型参数进行率定。2.根据权利要求1所述的基于改进新安江模型的实际径流量预报方法，其特征在于：所述研发半分布式天然径流量预报模块，具体通过数字流域构建、产流计算、汇流计算进行：1)所述数字流域构建过程是：(1)流域要素数字化：利用流域地形图等高线和实际水系数据，采用arcgis软件平台生成1km网格上的标量矩阵，以便在生成数字水系、子流域边界、数字坡度、数字河长要素时使用；(2)计算单元划分：运用矩形网格划分法或自然子流域-水文响应单元划分法，将流域细分为若干个计算单元，以考虑流域下垫面和气候因素的时空变异性；(3)计算单元雨量插值：采用降水空间分布点雨量插值方法，将雨量代表站实测雨量插值到各计算单元；2)所述产流计算是：采用新安江模型分别对上述数字流域的各计算单元产流过程进行半分布式模拟计算；所述汇流计算过程是：(1)计算单元汇流计算：对于覆盖子流域的每一个网格单元，计算网格单元的汇流路径、汇流时间；将汇流时间相的网格单元合并作为一个子流域，采用流时线法汇流模型进行子流域产汇流计算；流量计算公式如下：式中：q
i
为单元流域出流断面在第i时段出流量；h
i
为第i时段地面净雨；ω
j
为单元流域第j块流时面积；δt为时段长；(2)河系汇流计算：采用马斯京根分段连续演算法将每个计算单元的流量进行汇流演算到子流域出口并进行算术叠加，再将该子流域入流进行马斯京根河道汇流演算至子流域出口，将这两部分合并作为该子流域的出流，即天然径流。3.根据权利要求2所述的基于改进新安江模型的实际径流量预报方法，其特征在于：所述研发用水总量统计模块，对农业灌溉用水量和非农业灌溉用水量分别进行统计；具体按照如下步骤进行：(1)首先统计农业灌溉用水量，包括水稻灌溉用水量和旱作物灌溉用水量：
①
水稻灌溉用水量统计，包括灌排期用水量和插秧泡田期用水量；灌排期用水量计算公式为：
h
i+1
＝h
i
+p
i
+d
i-e
i-rcf
i
e
i
＝k
r
×
k1×
ep
i
插秧泡田期用水量计算公式为：δv＝(wm-w0+h
t
)/d+(r
f
+r
e
)w
i+1
＝w
i
+h
i
+p
i-et
i-rf
i
式中：h
i
为水稻田i时段初水深；h
min
、h
max
、h
p
为水稻生长的适宜水深下限、适宜水深上限和最大水深，按参考灌溉试验站数据取值；d
i
为水稻田i时段灌溉量；p
i
为i时段降雨量；e
i
为水稻i时段生长需水量；rcf
i
为水稻田i时段下渗量；ep
i
为蒸发皿i时段蒸发量；k1为水面蒸发量折算系数；k
r
为水稻生长需水系数；w
i
为水稻田i时段初土壤湿润层储水量；et
i
为i时段水面蒸发量；rf
i
为泡田期i时段渗漏量；δv为泡田期日最大灌水量；wm为田间持水量；w0为泡田期初土壤含水量；h
t
为泡田期末田间蓄水深；d为泡田期天数；r
f
为泡田期平均每天渗漏量；r
e
为泡田期平均每天水面蒸发能力量；
②
旱作物灌溉用水量统计；对于蓄满产流为主的南方湿润地区，旱作物灌溉用水量计算公式如下：w
i
'
+1
＝w
i
'+p
i-k
×
ep
i
式中：dh
i
为旱作物i时段灌溉用水量；w
max
为土壤湿润层旱作物生长适宜含水量；w
min
为土壤湿润层旱作物生长最小含水量，大于凋萎含水量，土壤湿润层含水量低于该值就应灌溉；w
i
'为旱田i时段初土壤湿润层储水量；k为流域综合蒸散发折算系数；β
i
为旱田i时段土壤湿润层相对于田间持水量的相对含水量；β
min
为农作物生长土壤湿润层最小相对含水量；β
max
为农作物生长土壤湿润层适宜相对含水量；(2)然后统计非农业灌溉用水量，包括生活、工业、林牧渔畜用水量：
①
生活用水量统计，包括城镇生活和农村生活用水量；城镇生活用水量根据城镇水厂取水计量监控数据，结合城镇生活供水比例进行统计；计算公式如下：农村生活用水量根据农村生活历史用水量数据；计算公式如下：式中：为城镇生活第i日用水量；为城镇水厂第i日取水量；k
sc
为城镇生活年供水比例；为农村生活第i年度第j天用水量，n为历史序列年份数，f(x)为用水量预测模型；
②
工业用水量统计，包括管网供水企业用水量和自备取水企业用水量；管网供水企业用水量从城镇水厂取水量中扣除城镇生活用水量得到；计算公式如下：自备取水企业用水量中在线监控的用水量根据监控数据直接统计；而离线计量的用水量根据历史用水量数据；计算公式如下：根据历史用水量数据；计算公式如下：式中：d
gg
为管网工业用第i日水量；在线监测的自备企业第i日用水量；第j家在线监测的自备企业第i日取水量；t为在线监测自备企业数量；为未在线监控的自备工业第i年度第j月用水量；
③
林牧渔畜用水量统计，包括林地灌溉、园地灌溉、草地灌溉、渔场养殖、大牲畜养殖、小牲畜养殖用水量；首先根据各类用水定额，结合经济统计指标对其年用水量进行统计，然后根据典型用水过程分配得到逐日用水量；计算公式如下：式中：为林牧渔畜第i日用水量，e
lmyc
为林牧渔畜养殖定额，p
lmyc
为林牧渔畜养殖数量，k
i
为林牧渔畜第i日用水量占全年用水量的比例；(3)最后对用水总量统计结果进行校正，具体过程如下：设置生活、工业、农业用水量校正参数分别为k
s
、k
g
、k
n
，以此来校准本次统计结果与水资源公报数据的偏差；采用ar自回归模型方法构建各行业用水量校正参数的预测模型；以参数k
g
为例，设其历史序列为则未来时段的校正参数计算公式为：式中：为t年度工业用水量校正参数；{c1,c2,c3,
…
,c
p
}为自回归参数系列；p为自回归阶数，值小于校正参数系列长度；ξ
t+1
为t+1年度校正计算后的系统残差。4.根据权利要求3所述的基于改进新安江模型的实际径流量预报方法，其特征在于：所述基于改进新安江模型研发流域水循环模型，计算实际径流量，具体按照如下步骤进行：1)搭建流域二元水循环结构；基于天然径流量预报模型对流域降水-产流-入渗-蒸发-汇流的自然水循环过程的模拟，再根据水源工程与各用水户之间的取供耗排关系，搭建流域自然-社会二元水循环结构；2)基于改进新安江模型研发流域水循环模型；水循环模型在天然径流量预报模块基础上，添加了小水库塘坝调蓄模拟模块、控制性水库调蓄模拟模块和取供耗排模拟模块；(1)小水库塘坝调蓄模拟；山丘区子流域内部建有较多小水库、小山塘及小池塘；为模拟上述蓄水工程对子流域产汇流过程的影响，将多个小型工程概化为一个虚拟水库，通过模拟虚拟水库的调蓄过程，将小水库塘坝的调蓄影响加入到新安江模型的产汇流过程，计算公式如下：
wr
i+1
＝wr
i
+r
fg
×
qr
i-sp
i-wx
i
式12式中：r
fg
为概化水库集水面积占子流域面积的比例；wr
i
为概化水库i时段初蓄水量；qr
i
为子流域i时段天然径流量；sp
i
为概化水库i时段供水量；wx
i
为概化水库i时段弃水量；w
max
为概化水库总库容；w
xl
为概化水库兴利库容；α为概化水库溢流系数；(2)控制性水库调蓄模拟；子流域出口控制性水库会显著改变子流域出流过程，进而影响河系汇流演算；为模拟控制性水库的影响，设置两种调蓄模拟方式，可择一选用：其一为水库调度规则模拟，即水库按照既定供水调度规则进行调蓄模拟；其二为水库自由蓄水模拟，即当水库蓄满才会向下游溢流；(3)取供耗排模拟；对水源供水过程、用户取水过程和用户退水过程进行模拟，具体过程如下：
①
水源供水过程模拟；基于水源天然径流量预报结果和其用水户用水总量统计结果，对水源向不同用水户的供水过程进行模拟，为水源实际径流量预报奠定基础；模型根据计算单元内某特定水源向城乡生活、工业、服务业、农业灌溉用水行业的实际供水次序进行调节计算，直至当日的可供水量依次全部分配给各用水部门；水源供水过程模拟计算公式如下：式中：d
i
和u
i
分别为计算单元内用水优先序为i的用水部门当日的用水量和供水量；n为用水部门个数；sp为计算单元当日的取水量；pip为城镇生活、工业供水管网漏损率；或农业灌溉水利用率；
②
用户取水过程模拟；基于用水户用水总量统计结果及其供水水源天然径流量预报结果，对用水户从不同水源的取水过程进行模拟；与水源供水过程模拟结合，实现区域水资源供用过程的完整模拟；模型根据计算单元内某特定用水户从水库水、塘坝水、河道水、外调水多种水源的实际取水顺序进行调节计算，直至满足该计算单元的当日需水量或者到最后一个供水源供水完毕为止；用户取水过程模拟计算公式如下：式中：sp
i
为计算单元中供水优先序为i的水源当日的取水量；m为供水源个数；d为计算单元当日的用水量；f
m
为取水工程的供水能力；sc
i
为水源当日的可利用水量；mx为区域取水限制总量；
③
用户退水过程模拟；计算单元内生活、工业、农业用水户取用水之后，经内部消耗后退水至河道中；用户退水计算公式如下：
式中：分别为i计算单元生活、工业、农业用水户当日退水量；分别为i计算单元生活、工业、农业用水户当日退水量；分别为i计算单元生活、工业、农业取水户当日取水总量、c
s
、c
g
、c
n
为生活、工业、农业用水户退水系数；3)利用流域水循环模型计算实际径流量；基于流域水循环模拟结果，特定水源工程集水区域实际径流量是在其天然径流量的基础上，扣除水源供水量，再加上区域用水户退水量，计算公式如下：qs
i
＝q
i-u
i
+b
i 式17式中：qs
i
为水源实际径流量；q
i
为水源天然径流量；u
i
为水源向外供水总量；b
i
为水源工程集水区域或上游区域内用水户退水量。5.根据权利要求4所述的基于改进新安江模型的实际径流量预报方法，其特征在于：(1)所述基础资料，包括水雨情资料、水利工程特征参数、取用水计量统计资料、作物种植结构、灌溉制度；水雨情资料包括日降水量、日均流量、日均水位、日均蒸发量；水利工程特征参数包括水库、塘坝的特征水位及库容和引调水工程的设计规模；取用水计量统计资料包括公共取水户、自备取水户、农饮水工程、农业灌区渠首取用水计量监测资料和公共水厂分行业供水量、林牧渔畜养殖数量统计资料；作物种植结构包括水稻、旱作物、经济作物各类农作物播种面积；灌溉制度为作物不同生育阶段的起止日期；(2)所述水循环模型参数，包括蒸发折算系数、张力水容量、自由水容量、壤中流出流系数、地下水库消退系数、稳定下渗率、渗透系数、河网蓄水消退系数、滞时共21个参数；水面日蒸发量折算系数、水稻生长需水系数、流域综合蒸散发折算系数、生活用水量校正参数、工业用水量校正参数、农业用水量校正参数共6个参数；生活用水户退水系数、工业用水户退水系数和农业用水户退水系数共3个参数。6.根据权利要求6所述的基于改进新安江模型的实际径流量预报方法，其特征在于：所述模型参数率定用人工经验率定或自动优化率定两种方法；自动优化率定方法以模型精度最优为目标，以各参数的取值范围为约束条件，构建模型参数优化模型；并采用群体复合形进化算法对模型进行求解，得到参数最优取值；模型参数优化模型描述如下：设模型的参数θ，将降雨过程p(t)、蒸发过程e(t)、用水户取水过程d(t)输入模型，使模型计算得到的径流过程与实测径流过程qs(t)误差平方和最小，表达式如下：
技术总结
本发明涉及一种水资源及节水管理方法。目的是提供一种基于改进新安江模型的实际径流量预报方法。该方法具有适用于中小流域、预报精度高的特点，可以为水资源及节水管理提供精确的实际径流量成果。技术方案是一种基于改进新安江模型的实际径流量预报方法，包括如下步骤：步骤1：基于新安江模型研发半分布式天然径流量预报模块，对天然径流量进行预报：所述研发半分布式天然径流量预报模块，包括数字流域构建、产流计算、汇流计算；计算流程参见图2。计算流程参见图2。计算流程参见图2。
技术研发人员：温进化 王贺龙 蔡晨凯 李进兴 戚核帅 杨辉斌
受保护的技术使用者：浙江省水利河口研究院（浙江省海洋规划设计研究院）
技术研发日：2022.12.16
技术公布日：2023/3/14