一种垂向和横向耦合混合产流计算方法技术

本发明专利技术属于径流预报相关技术领域，其公开了一种垂向和横向耦合混合产流计算方法，方法包括：根据下渗能力分布曲线获得超渗地表径流量和当前时段下渗量；构建流域蓄水容量分布曲线和相对不透水层下渗能力分布曲线；将流域蓄水容量分布曲线和相对不透水层下渗能力分布曲线耦合于同一坐标系，以两条曲线的交点和所述当前时段下渗量为基准获得总地下径流量；采用流域平均稳定下渗能力和流域稳渗分布不均匀系数将总地下径流量分割为壤中流径流量和地下径流量；将超渗地表径流量、壤中流径流量和地下径流量汇总获得流域出口总径流量。本申请解决了现有的混合产流方法开展洪水预报时存在径流产生及分割机制不清晰、洪水预报误差较大的技术问题。较大的技术问题。较大的技术问题。

【技术实现步骤摘要】
一种垂向和横向耦合混合产流计算方法


[0001]本专利技术属于径流预报相关
，更具体地，涉及一种垂向和横向耦合混合产流计算方法。

技术介绍

[0002]高精度水文预报可为水库调度、防洪减灾及水资源优化配置等提供决策依据，在水资源管理、水资源开发利用和国民经济建设等领域发挥着重要作用。产流计算是降雨径流建模的关键过程之一，受自然条件胁迫、陆面水文过程驱动，流域产流呈现异常复杂的时空变异特性，单一蓄满、超渗产流模式已经无法精确刻画产流过程的时空变异规律，开展基于混合产流的流域洪水预报研究是未来的重要发展方向。
[0003]目前常用的混合产流模式有垂向混合额产流、横向兼容混合产流、VIC
‑
3L混合产流等，地面径流、壤中流和地下水径流构成了不同产流机制下的全部径流成分，与单一产流模式相比，采用混合产流模式的预报精度有明显提升，但现有混合产流理论仍不足以精确刻画各种径流成分的组成变异规律，主要原因如下：现有混合产流理论假定土壤含水量达到田间持水量后产生壤中流，进而将采用蓄水容量分布曲线生产的蓄满径流量划分为壤中流和地下径流，忽略了未饱和带中也可能产生壤中流；其次，现有产流理论根据流域稳定下渗率进行地下径流分割时，其认为包气带岩土结构为水平方向空间分布均匀，该假定与实际情况不符，理论上，流域稳定下渗率是空间分布不均的，与流域植被覆盖、地形、土壤等下垫面特性相关。因此，现有的混合产流计算方法应用时存在径流产生机制不准确以及径流分割机制不合理，造成洪水预报精度偏低的技术问题。

技术实现思路
<br/>[0004]针对现有技术的以上缺陷或改进需求，本专利技术提供了一种垂向和横向耦合混合产流计算方法，其目的在于解决现有的混合产流方法开展洪水预报时存在径流产生及分割机制不清晰、洪水预报误差较大的技术问题。
[0005]为实现上述目的，按照本专利技术的一个方面，提供了一种垂向和横向耦合混合产流计算方法，所述方法包括：S1：根据下渗能力分布曲线获得超渗地表径流量，进而获得当前时段下渗量；S2：以流域平均蓄水容量以及流域蓄水容量分布不均系数为参数构建流域蓄水容量分布曲线，并以相对不透水层平均下渗率和相对不透水层下渗能力分布不均系数为参数构建相对不透水层下渗能力分布曲线；S3：将所述流域蓄水容量分布曲线和相对不透水层下渗能力分布曲线耦合于同一坐标系，以两条曲线的交点和所述当前时段下渗量为基准获得总地下径流量；S4：采用流域平均稳定下渗能力和流域稳渗分布不均匀系数将所述总地下径流量分割为壤中流径流量和地下径流量；S5：将所述超渗地表径流量、壤中流径流量和地下径流量汇总获得流域出口总径流量。
[0006]优选地，所述超渗地表径流量RS和当前时段下渗量FA的计算公式为：
[0007][0008]FA＝P
‑
RS
[0009]其中，P为降雨量，为下渗能力分布曲线围成的面积，表示流域平均下渗率；BF为下渗能力分布不均系数；Δt代表计算时段。
[0010]优选地，步骤S2还包括根据所述流域蓄水容量分布情况和相对不透水层下渗能力分布曲线分别获得流域蓄水容量小于某固定值WM
′
的面积比例α，以及根据相对不透水层下渗率小于某固定值f
s
的面积比例γ：
[0011][0012][0013]其中，WM
′
为流域点蓄水容量，WM为流域水容量分布曲线围成的面积，表示流域平均蓄水容量，B为流域蓄水容量分布不均系数，f
s
为相对不透水层点下渗能力，为相对不透水层下渗能力分布曲线围成的面积，表示相对不透水层平均下渗率，BF
s
为相对不透水层下渗能力分布不均系数。
[0014]优选地，步骤S3具体为：将所述流域蓄水容量分布曲线和相对不透水层下渗能力分布曲线耦合于同一坐标系，以下渗水量为输入，根据两条曲线的相交情况，计算总地下径流量RR：
[0015]①
当所述流域蓄水容量分布曲线和相对不透水层下渗能力分布曲线存在交点且则：
[0016][0017]②
当所述流域蓄水容量分布曲线和相对不透水层下渗能力分布曲线存在交点且则：
[0018][0019]③
当a+FA≤y，则：
[0020][0021]④
当所述流域蓄水容量分布曲线和相对不透水层下渗能力分布曲线不存在交点且a+FA＞WM(1+B)，则：
[0022]RR＝FA
‑
WM+W
[0023]其中，FA为当前时段下渗量，y为相对不透水层下渗能力分布曲线与流域蓄水容量分布曲线交点的纵坐标，a为当前时段流域平均土壤含水两的最大值，WM
′
为流域点蓄水容量，WM为流域水容量分布曲线围成的面积，表示流域平均蓄水容量，W为上一时段土壤含水量，α为流域蓄水容量小于某固定值WM
′
的面积比例，γ为相对不透水层下渗率小于某固定值f
s
的面积比例，WM
′
为流域点蓄水容量，f
s
为相对不透水层点下渗能力，B为流域蓄水容量分布不均系数，为相对不透水层下渗能力分布曲线围成的面积，表示相对不透水层平均下渗率。
[0024]优选地，步骤S4中采用流域平均稳定下渗能力和流域稳渗分布不均匀系数将所述总地下径流量分割为壤中流径流量RI和地下径流量RG具体公式为：
[0025][0026]RG＝RR
‑
RI
[0027]其中，RR为总地下径流量，为流域平均稳定下渗能力；BF
c
为流域稳渗分布不均系数，Δt为计算时段。
[0028]优选地，步骤S5中采用线性水库及马斯京根方法将所述超渗地表径流量、壤中流径流量和地下径流量汇总获得流域出口总径流量。
[0029]优选地，采用线性水库及马斯京根方法将所述超渗地表径流量、壤中流径流量和地下径流量汇总获得流域出口总径流量具体表达式为：
[0030][0031][0032][0033]QT(t)＝QS(t)+QI(t)+QI(t)
[0034]Q(t)＝C0·
QT(t)+C1·
QT(t
‑
1)+C2·
Q(t
‑
1)
[0035]式中，RS(t)、RI(t)、RG(t)分别代表t时段超渗地表径流量、壤中流径流量和地下径流量；QS(t)、QI(t)、QG(t)分别代表t时段超渗地表径流、壤中流和地下径流；QS(t
‑
1)、QI(t
‑
1)、QG(t
‑
1)分别代表t
‑
1时段超渗地表径流、壤中流和地下径流；CS、CI、CG分别代表超渗地表径流消退系数、壤中流消退系数和地下径流消退系数；F为流域面积；QT(t
‑
1)、QT(t)分别代表t
‑
1时段和t时段单元面积河网总入流；Q(t
‑
1)、Q(t)分别本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种垂向和横向耦合混合产流计算方法，其特征在于，所述方法包括：S1：根据下渗能力分布曲线获得超渗地表径流量，进而获得当前时段下渗量；S2：以流域平均蓄水容量以及流域蓄水容量分布不均系数为参数构建流域蓄水容量分布曲线，并以相对不透水层平均下渗率和相对不透水层下渗能力分布不均系数为参数构建相对不透水层下渗能力分布曲线；S3：将所述流域蓄水容量分布曲线和相对不透水层下渗能力分布曲线耦合于同一坐标系，以两条曲线的交点和所述当前时段下渗量为基准获得总地下径流量；S4：采用流域平均稳定下渗能力和流域稳渗分布不均匀系数将所述总地下径流量分割为壤中流径流量和地下径流量；S5：将所述超渗地表径流量、壤中流径流量和地下径流量汇总获得流域出口总径流量。2.根据权利要求1所述的方法，其特征在于，所述超渗地表径流量RS和当前时段下渗量FA的计算公式为：FA＝P
‑
RS其中，P为降雨量，为下渗能力分布曲线围成的面积，表示流域平均下渗率；BF为下渗能力分布不均系数；Δt代表计算时段。3.根据权利要求1所述的方法，其特征在于，步骤S2还包括根据所述流域蓄水容量分布情况和相对不透水层下渗能力分布曲线分别获得流域蓄水容量小于某固定值WM
′
的面积比例α，以及根据相对不透水层下渗率小于某固定值f
s
的面积比例γ：的面积比例γ：其中，WM
′
为流域点蓄水容量，WM为流域水容量分布曲线围成的面积，表示流域平均蓄水容量，B为流域蓄水容量分布不均系数，f
s
为相对不透水层点下渗能力，为相对不透水层下渗能力分布曲线围成的面积，表示相对不透水层平均下渗率，BF
s
为相对不透水层下渗能力分布不均系数。4.根据权利要求1或3所述的方法，其特征在于，步骤S3具体为：将所述流域蓄水容量分布曲线和相对不透水层下渗能力分布曲线耦合于同一坐标系，以下渗水量为输入，根据两条曲线的相交情况，计算总地下径流量RR：
①
当所述流域蓄水容量分布曲线和相对不透水层下渗能力分布曲线存在交点且则：
②
当所述流域蓄水容量分布曲线和相对不透水层下渗能力分布曲线存在交点且则：
③
当a+FA≤y，则：
④
当所述流域蓄水容量分布曲线和相对不透水层下渗能力分布曲线不存在交点且a+FA＞WM(1+B)，则：RR＝FA
‑
WM+W其中，FA为当前时段下渗量，y为相对不透水层下渗能力分布曲线与流域蓄水容量分布曲线交点的纵坐标，a为当前时段流域平均土壤含水两的最大值...

【专利技术属性】
技术研发人员：陈璐，易彬，陈桂亚，冯宝飞，许银山，张俊宏，郭鹤翔，刘一卓，冷致远，覃叶红萍，甘晓雪，李思明，杨彬林，
申请(专利权)人：长江水利委员会水文局，
类型：发明
国别省市：