汉江黄金峡站河系洪水预报研究.pdf

第10 期2023年10 月文章编号：16 7 3-9 0 0 0（2 0 2 3)10-0 0 7 3-0 4陕西水利Shaanxi Water ResourcesNo.10October,2023汉江黄金峡站河系洪水预报研究顾钊，宋淑红？，席小康？（1.西安水文水资源勘测中心，陕西西安7 10 0 6 8；2.陕西省水文水资源勘测中心，陕西西安7 10 0 6 8）摘要河系洪水预报对于流域防洪调度具有重要作用，在延长预报预见期上有显著优势。以汉江黄金峡断面以上流域为研究对象，采用新安江模型与垂向混合产流模型构建河系洪水预报方案，为检验模型选取与方案构建的合理性及预报精度，选取2 0 2 1年汛期的2 场典型洪水进行预报模拟。结果表明：以新安江模型与垂向混合产流模型构建的河系洪水预报精度良好，在洪峰预报与洪水过程预报上新安江模型略优于垂向混合产流模型，并且全流域暴雨洪水预报精度优于上游型暴雨洪水，预报预见期可延长10 h20h，可以为汉江汉中段及黄金峡、石泉、安康水库防洪调度决策提供重要依据。【关键词河系洪水预报；新安江模型；垂向混合产流模型；汉江；黄金峡中图分类号P338【文献标识码B流计算包括坡面汇流和河网汇流2-3。新安江模型对湿润地区1类数据来源与研究方法1.1研究流域概况汉江是长江中游的重要支流，发源于汉中市宁强县潘家山，流经汉中市宁强、勉县、南郑、汉台、城固、洋县、西乡7 县（区），于西乡县茶镇进入安康市，汉中市境内集水面积19 6 9 2 km，其中黄金峡水文站以上流域集水面积17 0 7 0 km多年平均降水量7 0 0 mm1700mm,属于典型的湿润地区，流域内暴雨洪水受季风影响较大，7 月 9 月易形成大洪水。汉江黄金峡以上流域水系发育、河网密度大、水利工程较多，集水面积大于10 0 km的一级支流有2 2 条，干流设有武侯镇、汉中、洋县、黄金峡4 处水文站，主要支流设有11处基本水文站，流域内建有石门、黄金峡2 座大型水库，另有多座中小型水库。1.2数据来源本研究选择的资料为汉江黄金峡断面以上流域19 9 5年2 0 2 1年汛期（4 月 10 月）整编的降雨、流量资料，从19 9 5年 2 0 2 0 年中筛选出3 9 场洪水，按时间前后顺序选取前2 9场为率定场次、其余10 场作为检验场次，分别进行场次洪水率定及模型验证；从2 0 2 1年汛期洪水场次中选取2 场典型洪水进行河系洪水预报应用分析。1.3研究方法1.3.1新安江模型新安江模型是由赵人俊教授等人提出的一种具有物理基础的概念性流域水文模型。该模型主要包括四部分：蒸散发计算、产流计算、水源划分以及汇流计算。模型采用三层蒸散发模式计算蒸发量；以蓄满产流方式计算产流量；将径流分为地面径流、地下径流和壤中流三类进行水源划分；汇与半湿润地区的洪水预报及径流模拟具有很好的模拟效果4-7。1.3.2垂向混合产流模型垂向混合产流模型是包为民教授8 提出的一种在垂向上将超渗产流与蓄满产流进行组合的流域水文模型。当流域内降雨到达地面，首先通过空间分布的下渗能力分布曲线，划分为地面径流和下渗水流；地面径流取决于雨强和下渗，为超渗产流模式；地面以下的径流，取决于前期土壤缺水量和下渗的水量，包含壤中流和地下径流，为蓄满产流模式。模型仍然采用三层蒸散发模式计算蒸发量；水源划分则按开式自由水箱按比例进行划分计算，坡面汇流采用线性水库计算，河道汇流采用马斯京根法。2河系洪水预报方案构建根据汉江汉中段干支流的洪水传播演进规律，采用新安江模型和垂向混合产流模型分别构建武侯镇、汉中、洋县、黄金峡4 处干流水文站，以及沮水河茶店、胥水河升仙村2 处支流水文站与褒河石门水库入库的单站预报方案。在单站预报的基础上，采用上断面的预报结果代替实测值进行下断面的连续预报，进一步构建了汉江源头至黄金峡的河系预报方案。2.1干流4 站预报方案构建（1）汉江源头武侯镇将汉江源头武侯镇区间划分为沮水河茶店站以上流域、玉带河铁锁关站以上流域、汉江源头至武侯镇站未控流域共3 个单元。沮水河茶店站以上流域、汉江源头至武侯镇未控流域2 个子区间的产流部分分别采用两种模型的产流模块计算；玉带河铁锁关站以上流域以及上述2 个子区间的汇流部分采用分段马斯京根连续演算，根据各子区间汇人河道收稿日期 2 0 2 3-0 8-0 3基金项目国家重点专项（2 0 18 YFC1508100）【作者简介顾钊（19 8 9-），男，江苏泰州人，工程师，主要从事洪水预警预报方面的工作。.73第10 期2023年10 月位置，进行先合后演，逐河段演算至武侯镇站断面。(2）武侯镇汉中将武侯镇汉中区间划分为干流武侯镇站以上流域、漾家河元墩站以上流域、濂水河江西营站以上流域、褒河石门水库以上流域、武侯镇至汉中未控流域共5个子区间。褒河石门水库以上流域、武侯镇至汉中未控流域2 个子区间的产流部分分别采用两种模型的产流模块计算；干流武侯镇站以上流域、漾家河元墩站以上流域、濂水河江西营站以上流域3个子区间及上述2 个子区间的汇流部分采用分段马斯京根连续演算，根据各子区间汇人河道位置，进行先合后演，逐河段演算至汉中站断面。(3)汉中 洋县将汉中洋县区间划分为干流汉中站以上流域、冷水河三华石站以上流域、溢水河长滩村站以上流域、胥水河升仙村站以上流域、汉中至黄金峡未控流域共5个子区间。胥水河升仙村站以上流域、汉中至黄金峡未控流域2 个子区间的产流部分分别采用两种模型的产流模块计算；干流汉中站以上流域、冷水河三华石站以上流域、溢水河长滩村站以上流域3 个子区间以及上述2 个子区间的汇流部分采用分段马斯京根连续演算，根据各子区间汇人河道位置，进行先合后演，逐河段演算至洋县站断面。（4）洋县 黄金峡洋县黄金峡区间划分为干流洋县站以上流域、酉水河酉水街站以上流域、金水河金水站以上流域和洋县至黄金峡未控流域共4 个子区间。金水河金水站以上流域子区间以及洋县至黄金峡未控流域子区间的产流部分分别采用两种模型的产流模块计算；干流洋县站以上流域、酉水河酉水街站以上流域2 个子区间及上述3 个子区间的汇流部分采用分段马斯京根连续演算，根据各子区间汇人河道位置，进行先合后演，逐河段演算至黄金峡水库坝址断面。玉带河铁锁关武侯镇漾家河元墩赛河石门水库江西营谦水河汉中冷水河三华石背水河溢水河升仙村长滩村洋县酉水河西水街金水河金水黄金峡图1汉江汉中段单站洪水预报拓扑结构表1新安江模型次洪方案参数预报断面KB茶店0.90武侯镇0.90石门水库0.90汉中0.90升仙村0.90洋县0.90黄金峡0.90.74陕西水利Shaanxi WaterResources2.2河系预报方案构建以汉江干流黄金峡水文站为预报对象，根据汉江汉中段流域内子流域的上下游拓扑关系，串联组合单站预报方案，并依据历史场次洪水进一步优化单站预报中的新安江模型与垂向混合产流模型汇流参数，分别构建汉江源头、武侯镇站、汉中站至黄金峡站的河系预报方案，组合结构见图2。武侯镇汉中汉江源头至黄金峡站河系预报图2 汉江源头、武侯镇站、汉中站至黄金峡站河系洪水预报方案结构2.3参数率定与结果分析对单站次洪水分别进行新安江模型与垂向混合产流模型参数率定，茶店、石门水库、升仙村、武侯镇、汉中、洋县站选择19 9 0 年 2 0 14 年资料用于参数率定，2 0 15年 2 0 2 0 年资料用于参数检验；黄金峡站选择2 0 10 年 2 0 18 年资料用于参数率定，2 0 19 年 2 0 2 0 年资料用于参数检验。两种模型的计算步长均为1h。模型参数率定采用人工优选法和SCE-UA自动优选法相结合的方式。对于土壤蓄水容量分布曲线指数B、深层蒸散发扩散系数C、流域不透水面积比例IM、土壤自由水容量分布曲线指数Ex、地下水出流系数Kc、地下水汇流系数Cc、流域下渗率分布曲线指数Br、土壤对下渗率影响系数K,等相对不敏感参数，根据流域特征与预报经验使用人工优选法沮水河对其进行率定；对于流域蒸散发折算系统K、土壤自由水蓄茶店水容量SM、地面水汇流系数Cs和流域实际下渗率Fc等敏感参数，采用SCE-UA自动优选法率定；次洪模型中张力水蓄构建单站预报实测输人CWMWuMWMIM0.200.130.300.130.200.130.300.130.200.130.300.130.300.13No.10October,2023武侯镇至黄金峡站河系预报洋县汉中至黄金峡站河系预报水容量WM、上层张力水蓄水容量WuM和下层张力水蓄水容量WLuw数值通过日模型计算0两种参数率定方法相结合有效避免了人工优选参数慢、参数精度低的缺点，同时可以避免自动优选参数取值不唯一、易陷入局部最优解的缺点。表1为新安江模型次洪模拟参数率定结果，表2 为垂向混合产流模型次洪模拟参数率定结果。表3 统计了两种模型的单站洪水模拟精度评定指标，表4 统计了两种模型的河系洪水模拟精度评定指标。SMEx15020150201502015020150201502015020黄金峡KcKCc800.0077.00800.007800.0077.00800.00722.0800.0077.00800.007800.009Cs1.300.308.001.301.301.301.3018.01.308.001.30Ke0.300.980.300.400.300.300.300.400.300.300.300.400.300.40X0.920.750.970.900.980.920.990.930.980.920.990.930.990.931.100.831.300.751.100.871.000.751.100.871.500.871.100.380.380.380.380.380.380.38第10 期2023年10 月预报断面K茶店0.90武侯镇0.90石门水库0.90汉中0.90升仙村0.90洋县0.90黄金峡0.90表3单站洪水模拟精度评定指标统计总评价指标率定期检验期单站预模型洪水Ro4H报合格合格DC场次场次率%率1%XAJ2085850.84茶店MIXXAJ3989.741000.88111001000.92武侯镇MIX石门XAJ1687.581.250.851090800.8水库MIXXAJ2993.11000.82131001000.79汉中MIXXAJ1610081.250.81升仙村MIXXAJ5298.081000.9125100 1000.94洋县MIXXAJ2910096.550.8610901000.85黄金峡MIXXAJ20994.1893.01平均值MIX20987.8794.290.81表4河系洪水模拟精度评定指标统计总评价指标率定期检验期单站模型洪水R4H预报合格合格DC场次场次率/%率/%源头 黄XAJ2996.5596.550.8310801000.83金峡MIX武侯镇XAJ2996.551000.8410801000.84黄金峡MIX汉中 黄XAJ2910096.550.8410901000.84金峡MIXXAJ17492.9992.250.819088.1594.820.85平均值MIX由两种模型的次洪模拟参数结果可以看出，7 站同处于汉江汉中段流域，年降雨量比较接近，下垫面条件大致相同，因此7 站的大部分模型参数差别不大，但是在新安江模型次洪模拟中汉中、洋县两站的土壤自由水蓄水容量SM取值较大，这主要是由于武侯镇以下至洋县小峡口为汉江平川段，土层陕西水利Shaanxi Water Resources表2垂向混合产流模型次洪方案参数BCWMWuM0.350.130.390.130.390.130.390.130.390.130.300.130.390.1380950.7187.1897.440.877581.250.7789.661000.827587.50.896.151000.8910093.10.840.848893.4491.7785.7393.330.8洪水RQ4H合格合格DC率/%率1%93.196.550.8293.196.550.821001000.8290.690.390.8No.10October,2023MSM15525150201653515020155251502015020洪水RQ4H合格合格DC率/%率1%887.587.50.887.51000.731001000.9370800.8492.311000.76610066.670.715066.670.65961000.92901000.770.83801000.81801000.83901000.8384.8298.150.83Ex80158020802080201.20801580268020Kc1.200.401.200.401.200.400.401.200.401.200.501.200.40较厚，土壤的蓄水量较山地明显增大。由次洪模拟精度评定结果可以看出，在率定期与检验期，单站预报与河系预报的洪峰误差、峰现时间误差的合格率以及确定性系数均达到乙等方案精度要求，表明了两种模型在汉江汉中段流域具有较好的适用性，也说明了单站预报与河系预报方案构造合理，同时新安江模型的洪水模拟精度高于垂向混合产流模型，说明了汉江汉中段流域产流方式以蓄满产流为主，超渗产流占比不高。在单站次洪模拟中，支流模拟精度小于干流，说明了支流洪水受中小型水利工程蓄泄影响较大，在蓄泄作用下，支流单站的场次洪水量级与洪水过程较天然状态下有差异，导致洪峰模拟合格率与确定性系数下降，同时还有雨量站数量不足与分布不均的影响。在干流次洪模拟中，洋县、黄金峡站模拟效果优于武侯镇及汉中站，说明了汉中洋县区间、洋县黄金峡区间来水对洋县、黄金峡站洪水的贡献小于上断面的干流洪水演进。黄金峡河系次洪模拟与单站次洪模拟对比，在洪峰模拟中，汉中黄金峡河系方案模拟精度与单站模拟精度基本一致，汉江源头黄金峡、武侯镇黄金峡两种河系方案精度较单站下降了7%，在峰现时间与过程拟合中，河系方案模拟精度与单站方案接近。主要原因在于河系连续预报中模型计算误差传递对精度影响较大，单站预报只进行了1次计算，汉江源头、武侯镇、汉中黄金峡站河系预报分别进行了4次、3 次、2 次连续计算，同时石门水库集水面积占汉中站控制流域面积的4 1.4%，汉中洪水主要来自干流武侯镇站与区间支流褒河，汉中洪水受石门水库蓄泄影响较大，部分场次洪水模拟精度较低，因此以汉中为预报依据站的河系方案较优，也说明了黄金峡站洪水主要来自洋县以上流域，洋县黄金峡区间来水贡献较低。3应用结果与分析2021年8 月10 月汉江汉中段出现连续7 场暴雨洪水，从中选取了“2 0 2 10 8 2 1”“2 0 2 10 9 2 5”两场典型洪水进行单站及河系预报计算，计算步长均为1h，源头黄金峡河系预报开始时间选择过程降雨开始时刻，武侯镇黄金峡、汉中黄金峡河系预报以及洋县黄金峡单站预报开始时间选择上边界的峰现时刻。表5为两场洪水河系预报结果评定指标统计，图4、图5分别是“2 0 2 10 8 2 1”“2 0 2 10 9 2 5”号洪水河系预报过程线。.75KCc0.300.390.310.390.970.310.390.310.390.310.390.310.390.310.39CCs0.920.850.920.970.920.990.930.970.920.800.920.970.92Ke0.701.300.801.100.801.100.871.000.801.100.971.500.801.10X0.380.480.480.480.480.480.48参考文献第10 期2023年10 月表5河系预报结果评定指标统计洪峰相对峰现时间确定性系数误差/%误差/小h河系预报洪号XAJMIXXA.JMIXXAJMIX源头2021082112.917.1黄金峡20210925-0.21.7武侯镇黄2021082112.114.7金峡20210925-0.4-1.12021082110.8汉中黄金峡20210925-0.6-0.1洋县 黄金20210821峡(单站)20210925-0.11.8一面平均雨量一实测流量一面平均雨量一实测流量XAJ计算流量-MIX计算流量XAJ计算流量-MIX计算流量8000070001060005000=400030005200010008/218/228/238/248/258/268/278/28（a）源头黄金峡河系预报过程线面平均雨量，一实测流量XAJ计算流量-MIX计算流量800007000060005000=40003000200010008/218/228/238/248/258/268/278/28（c）汉中黄金峡河系预报过程线图4“20210821”号洪水河系预报过程线一面平均雨量一实测流量一面平均雨量，一实测流量XAJ计算流量-MIX计算流量XAJ计算流量-MIX计算流量12000120007.100005“600080004000200009/24.9/259/269/279/289/299/30 10/110/2（a）源头黄金峡河系预报过程线一面平均雨量一实测流量-XAJ计算流量-MIX计算流量12000.10008000“6000面4 0 0 0200009/24.9/259/269/279/289/299/3010/110/2（c）汉中黄金峡河系预报过程线图5“20210925”号洪水河系预报过程线应用结果表明：在黄金峡站单站和河系预报中，两场洪水预报误差均在合格范围内，新安江模型预报结果的洪峰相对误差和确定性系数均好于垂向混合产流模型，并且黄金峡单站方案在洪峰预报中精度优于河系方案，汉中黄金峡河.76.陕西水利Shaanxi WaterResources系方案预报精度优于汉江源头黄金峡、武侯镇黄金峡方案。从黄金峡站两场洪水的河系预报结果分析，“2 0 2 10 9 2 5”洪水较“2 0 2 10 8 2 1”预报误差更小，研究认为：“2 0 2 10 8 2 1”场次中黄金峡站洪水均主要来自汉中站，汉中站至黄金峡站-1-10.9330.905-1-10.9680.906220.9330.925010.9180.96613009.912.1800070006000205000=4000300040200050100008/218/228/238/248/258/268/278/28时间（b）武侯镇黄金峡河系预报过程线一面平均雨量一实测流量XAJ计算流量-MIX计算流量8000700010?6000205000=400030200030004010001508/218/228/23 8/24 8/25 8/26 8/278/280时间（d）洋县黄金峡单站预报过程线,100002510=60008000400020002509/24.9/259/269/279/289/299/3010/110/2时间（b）武侯镇黄金峡河系预报过程线一面平均雨量，一实测流量XAJ计算流量-MIX计算流量120000uu/本510000800010=60001540002020002509/249/259/269/279/28 9/29 9/3010/110/2时间（d）洋县黄金峡单站预报过程线No.10October,2023区间流域产流不多，属于典型的上游型洪水，而“2 0 2 10 9 2 5”场次属于全流域暴雨洪水，在“2 0 2 10 8 2 1”场次中模型计算的区间产流量较实际偏大，同时对于上游型洪水来讲汉中站预报误差对下游洪水预报精度影响较全流域暴雨洪水更明00.9330.92500.9380.916110.9330.924110.9680.966050时间时间时间时间显。因此，与汉江汉中段上游型洪水相比，河系预报在全流域暴雨洪水应用中优势更明显。4结论本文结合历史洪水数据，采用新安江模型与垂向混合产流模型建立的河系预报方法在汉江黄金峡站洪水的洪峰预报、洪量预报、峰现时间误差上模拟精度良好，同时新安江模型总体优于垂向混合产流模型，河系预报在全流域暴雨洪水模拟精度优于上游型洪水。较单站预报方案，河系方案洪峰模拟精度下降了7%，峰现时间和洪水过程拟合中河系方案模拟精度与单站接近，河系方案预见期较单站延长了10 h20h。因此，汉江黄金峡河系预报在满足了流域内防洪需求之上，更大程度地延长了洪水预见期，可为下游石泉、安康水库的防洪调度提供更精确、更提前的洪水预报信息。但是，值得注意的是，石门水库调度以及流域暴雨洪水类型是汉中及以下断面河系洪水预报误差产生的重要因素，这正是下一步的研究内容。1范火生.基于栅格的分布式新安江模型在三岔河流域的洪水预报研究J.中国农村水利水电,2 0 19（2)：113-118.2ZHAO R J.The Xinanjiang model applied in ChinaJ.Journal ofHydrology,1992,135(1-4):371-381.3芮孝芳，凌哲，刘宁宁，等.新安江模型的起源及对其进一步发展的建议J.水利水电科技进展，2 0 12,（4）：1-5.4刘佩瑶，郝振纯，王国庆，等.新安江模型和改进BP神经网络模型在闽江水文预报中的应用J.水资源与水工程学报,2 0 17,(1):40-44.5姚成，李致家，张珂，等.基于栅格型新安江模型的中小河流精细化洪水预报J.河海大学学报（自然科学版),2 0 2 1,4 9(1):19-2 5.56宋晓猛，孔凡哲.新安江模型和人工神经网络的耦合应用J.水10152025土保持通报,2 0 10,3 0(6):13 5-13 8,14 4.7姚成，章玉霞，李致家，等.无资料地区水文模拟及相似性分析J.河海大学学报（自然科学版）,2 0 13,4 1（2）：10 8-113.8包为民，王从良垂向混合产流模型及应用J.水文,19 9 7,(3):18-21.9瞿思敏，包为民，张明，等。新安江模型与垂向混合产流模型的比较J.河海大学学报（自然科学版）,2 0 0 3,3 1（4）：3 7 4-3 7 7.10】霍文博.新安江模型和支持向量机模型实时洪水预报应用比较J.河海大学学报（自然科学版）,2 0 18，（4）：2 8 3-2 8 9.