CN111931440B - 一种考虑地表积水深变化的坡面多流向漫流分析方法 - Google Patents
一种考虑地表积水深变化的坡面多流向漫流分析方法 Download PDFInfo
- Publication number
- CN111931440B CN111931440B CN202010885571.3A CN202010885571A CN111931440B CN 111931440 B CN111931440 B CN 111931440B CN 202010885571 A CN202010885571 A CN 202010885571A CN 111931440 B CN111931440 B CN 111931440B
- Authority
- CN
- China
- Prior art keywords
- water
- flow
- unit
- grid unit
- grid
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 39
- 239000002352 surface water Substances 0.000 title claims abstract description 29
- 238000005206 flow analysis Methods 0.000 title claims abstract description 6
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 145
- 238000004088 simulation Methods 0.000 claims abstract description 59
- 238000009826 distribution Methods 0.000 claims abstract description 20
- 238000009825 accumulation Methods 0.000 claims abstract description 14
- 239000002689 soil Substances 0.000 claims description 45
- 238000004364 calculation method Methods 0.000 claims description 36
- 230000008595 infiltration Effects 0.000 claims description 34
- 238000001764 infiltration Methods 0.000 claims description 34
- 238000009736 wetting Methods 0.000 claims description 24
- 230000001186 cumulative effect Effects 0.000 claims description 9
- 238000004458 analytical method Methods 0.000 claims description 6
- 229920006395 saturated elastomer Polymers 0.000 claims description 6
- 238000004422 calculation algorithm Methods 0.000 claims description 5
- 238000001704 evaporation Methods 0.000 claims description 4
- 230000008020 evaporation Effects 0.000 claims description 4
- 238000001556 precipitation Methods 0.000 claims description 3
- 238000009834 vaporization Methods 0.000 claims description 3
- 230000008016 vaporization Effects 0.000 claims description 3
- 239000013256 coordination polymer Substances 0.000 claims description 2
- 230000005855 radiation Effects 0.000 claims description 2
- 238000004519 manufacturing process Methods 0.000 abstract description 7
- 238000006243 chemical reaction Methods 0.000 abstract description 3
- 230000002093 peripheral effect Effects 0.000 abstract description 3
- 238000010276 construction Methods 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000009792 diffusion process Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 229920002472 Starch Polymers 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 235000019698 starch Nutrition 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Fluid Mechanics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computing Systems (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及一种考虑地表积水深变化的坡面多流向漫流分析方法,包括:地表积水深的确定、坡面动态汇流路径确定、坡面多流向水量分配权重确定和坡面多流向漫流水量确定。本发明针对平原洼地区地表起伏小,积水面积大、变化快的特点,实现了对地表积水深动态变化及其与周边单元的水量流入流出关系转化的准确刻画,解决了常规数字高程模型应用失真的问题,精细模拟了坡面多流向漫流过程,有效改善了平原低洼区地表产汇流过程模拟。该方法面向平原洼地区水文模拟,相关结果可为平原洼地区水文模型构建与应用提供基础支撑。
Description
技术领域
本发明涉及水利工程技术领域,尤其涉及平原低洼区地表产汇流过程模拟,具体为一种考虑地表积水深变化的坡面多流向漫流分析方法。
背景技术
平原低洼区具有地势低平、地表水流流向随机的特点,流向多不单一且不固定,易受地表积水变化影响而发生改变,长期以来对平原低洼区地表产汇流过程的准确模拟是个难题。为此,坡面多流向漫流理论被提出并不断完善,其基本思想是认为水流应向邻域中所有高程较低的方向进行分配。目前,相关理论方法均是依靠常规数字地形数据(一般为数字高程模型“DEM”),研究坡度、地表形态、粗糙度等对地表水流多流向流量分配的影响。然而,从平原区特别是平原湿地沼泽区来说,当地洼地星罗棋布,地表起伏小,积水面积大、变化快,采用传统不变的数字高程模型(DEM)难以准确反映地形地貌变化,会对流向和流量确定产生较大的干扰,影响地表产汇流过程模拟效果。以洼地蓄滞水过程为例,随着洼地积水,水流流出单元与汇入单元间的水力梯度逐渐得到消弭,直接影响到单元出流水量。
发明内容
本发明所要解决的技术问题,在于克服现有平原低洼区坡面多流向漫流过程模拟技术存在的缺陷,构建考虑地表积水深度变化的动态数字地形,并基于此创造性提出包含坡面动态汇流路径确定和多流向漫流水量计算的坡面汇流计算方法。该发明可有效刻画平原洼地积水变化及其对地表产汇流过程的影响,提高平原低洼区地表产汇流过程模拟精度。
本发明的目的是这样实现的:
一种考虑地表积水深变化的坡面多流向漫流分析方法,所述方法包括四个部分:地表积水深确定、坡面动态汇流路径确定、坡面多流向水量分配权重确定和坡面多流向漫流水量确定,各部分的具体步骤如下:
地表积水深确定:
步骤一:获取研究区内气象站点在模拟期内的气象数据,所述气象数据包括逐日降水、气温、风速、日照、湿度,选择参证气象站,利用反距离平方加权插值算法将参证气象站点的气象数据空间展布到研究区模拟栅格单元上;
步骤二:采用Green-Ampt入渗模型模拟计算栅格单元的降雨日累计入渗量;
步骤三:利用Penman公式模拟计算栅格单元地表当日最大可蒸发量;
步骤四:利用水量平衡方程计算模拟单元日积水深,计算公式如下:
式中,H为日积水深/mm;H0为当日初始积水深/mm;P为当日降水深/mm;E为当日最大可蒸发量/mm;F为降雨日累计入渗量/mm;
坡面动态汇流路径确定的具体步骤:
步骤1:以日为汇流模拟时段,赋予研究区内每个模拟栅格单元5个基本属性:坐标x、坐标y、地形高程z、日积水深H、是否已参与汇流计算;已参与汇流计算标记为ture,否则标记为False;
步骤2:逐汇流模拟时段内,引入当前模拟栅格单元日积水深H,将其与模拟栅格单元的地形高程z相加,以此作为坡面漫流、河道汇流和漫溢计算所需的“虚拟高程”,记作DH;
步骤3:遍历模拟栅格单元,寻找DH最大值所在的模拟栅格单元,确定流向和出流量,然后更新模拟栅格单元的DH,并将结束汇流的模拟栅格单元标记为True,即表示本时段内该栅格不再参与汇流计算;
步骤4:逐次遍历,直至所有模拟栅格单元的汇流路径和出流量全部计算完成;
坡面多流向水量分配权重确定的具体步骤:
步骤1):比较模拟栅格单元和该单元相邻的8个栅格单元的DH大小,若模拟栅格单元DH大于相邻栅格单元,则认为水流流向该临近栅格单元,研究认为模拟栅格单元的水流应向邻域中所有DH较低的方向进行分配,由此对流向进行编码;
步骤2):结合地表Manning糙率,确定模拟栅格单元在临近多个DH较低栅格单元上的水量分配权重;
坡面多流向漫流水量确定的具体步骤:
步骤Ⅰ:使用运动波模型计算水流流出栅格单元在本汇流时段内的流量,记为可出流量;计算公式如下:
Sf=S0 (20)
式中,A为过流断面面积,单位为m2;Q为断面流量,单位为m3/s,即为可出流量;qL为单宽流入量,单位为m3/s/m;Sf为摩擦比降;S0为地面坡降或河道坡降;R为过流断面水力半径,单位为m;n为坡面曼宁糙率系数;t为时间;x为河段长度;
步骤Ⅱ:计算水流流出栅格单元与流入栅格单元“虚拟高程”出现持平时的出流量dc,以此作为该计算单元在本汇流时段的最大出流量;
步骤Ⅲ:取Q与dc两者之间的较小值作为单元出流量最终模拟值,并按水量分配权重对各方向单元流入量进行修正。
进一步的优化,步骤一中利用反距离平方加权插值算法将参证气象站点的气象数据空间展布到研究区模拟栅格单元上,计算公式如下:
式中,z和M分别表示气象站的高程和气象数据值,dM/dz为气象数据随高程增加的线性变化率;m表示参证气象站点的总数;D和Di分别为模拟栅格单元和第i个参证气象站点的气象数据值;z0和zi分别为模拟栅格单元和第i个参证气象站点的高程;λi为第i个参证气象站点所占权重;x0和xi分别为模拟栅格单元和第i个参证气象站点的横坐标;y0和yi分别为模拟栅格单元和第i个参证气象站点的纵坐标;di为模拟栅格单元和第i个参证气象站点的平面距离;
进一步的优化,步骤二中采用Green-Ampt入渗模型模拟计算栅格单元的降雨日累计入渗量,计算方法如下:
当入渗湿润锋到达第m层土壤时,入渗能力f为:
若湿润锋进入第m-1层土壤时地表持续积水,降雨日累计入渗量F由式(5)计算;若前一时段tn-1地表无积水,当前时段tn开始积水,降雨日累计入渗量F则用式(6)计算;
其中,
Fp=Am-1(Ip/km-1)-Bm-1 (10)
tp=tn-1+(Fp-Fn-1)/Ip (11)
Δθi=θs,i-θi,0 (12)
式中,m为湿润锋最终到达的土层数;i表示湿润锋当前所在土层数;km为湿润锋所在m层时的土壤导水系数,单位为mm/day;ki为第i层土壤导水系数,单位为mm/day;Δθi为第i层土壤含水量的变化量,cm3/cm3;Δθm为湿润锋所在m层的土壤含水量变化量,cm3/cm3;Am-1和Bm-1为参数,无量纲;Fm-1为m-1层以上的土壤累积入渗量,mm;FP为当前时段开始积水时的累积入渗量,单位为mm;t为时刻;tn-1为第n-1时段,min;Fn-1为第n-1时段的土壤累积入渗量,mm;tm-1为湿润锋进入m-1和m层土壤界面时的时间;tp为开始积水时刻;Li为第i层土厚,单位为mm;SWm为第m层湿润锋处的毛管吸力,单位为mm;θ表示土壤饱和含水量,cm3/cm3;Ip为积水开始时的雨强,mm/min;θs,i为第i层土壤饱和含水量,cm3/cm3;θi,0为土壤初始含水量,cm3/cm3。
进一步的优化,步骤三:利用Penman公式模拟计算栅格单元地表当日最大可蒸发量,计算公式如下:
式中:E为当日最大可蒸发量,单位为mm;Δ为饱和水汽压曲线斜率,单位为kPa/℃;RN为净放射量,单位为MJ/m;G为地中热传导,单位为MJ/m2;ρa为空气密度,单位为kg/m3;CP为空气定压比热,单位为J/kg/℃;δe为实际与饱和水蒸气压之差,单位为kPa;ra为空气动力学阻抗,单位为s/m;λ为水体的气化潜热,MJ/kg;γ为空气湿度常数,kPa/℃;Pa为大气压,kPa。
进一步的,步骤2):确定模拟栅格单元在临近多个DH较低栅格单元上的水量分配权重,计算公式如下:
式中,i和j表示两个不同栅格单元;DHi和DHj分别为栅格单元i和栅格单元j的虚拟高程;wji表示从DH较高栅格单元j到相邻DH较低单元i的分配水量权重;sji表示栅格单元j到栅格单元i间的水力梯度;nji为栅格单元j到栅格单元i间的坡面曼宁糙率系数;ni和nj分别为栅格单元i和栅格单元j的坡面曼宁糙率系数;xi和xj分别为栅格单元i和栅格单元j的横坐标;yi和yj分别为栅格单元i和栅格单元j的的纵坐标。
本发明产生的有益效果是:
本发明一种考虑地表积水深变化的坡面多流向漫流分析方法,在模拟坡面多流向漫流过程中考虑地表积水深的变化及其影响,可精细刻画平原洼地单元的积水变化及其与周边单元的水量流入流出关系转化,提高平原低洼区地表产汇流过程模拟精度。
本发明方法创造性提出了坡面动态汇流路径确定和坡面多流向漫流水量确定的具体流程,准确描述了平原洼地区域地表起伏小,积水面积大、变化快的特点,深入刻画了洼地对坡面径流的调蓄作用以及平原地表径流多流向漫散过程,实现了对地表积水深动态变化及其与周边单元的水量流入流出关系转化的准确刻画,解决了常规数字高程模型应用失真的问题,有效改善了平原低洼区地表产汇流过程模拟。该方法面向平原洼地区水文模拟,相关结果可为平原洼地区水文模型构建与应用提供基础支撑。
附图说明
下面结合附图和具体实施例,对本发明作进一步详细说明。
图1是本发明的实施例1所述方法的流程图;
图2是本发明考虑地表积水深的动态汇流路径确定示意图;
图3是本发明不同流向对应编码示意图;
图4是本发明地表漫流多流向水量分配计算示意图。
具体实施方式
本实施例是一种考虑地表积水深变化的坡面多流向漫流分析方法,计算流程如图1所示。本实施例所述方法的计算过程如下:
第一步:计算确定研究区模拟栅格单元的地表积水深
1、获取研究区内气象站点在模拟期内的逐日降水、气温、风速、日照、湿度等气象数据,然后利用反距离平方加权插值算法将气象站点数据空间展布到研究区模拟栅格单元上。计算公式如下:
式中:z和M分别表示气象站的高程和气象数据值,dM/dz为气象数据随高程增加的线性变化率,若采用平面插值则dM/dz=0;下标0表示待插值栅格点,i表示参证气象站点;m表示参证气象站点的总数;D和Di分别为模拟栅格单元和第i个参证气象站点的数据值;z0和zi分别为模拟栅格单元和第i个参证气象站点的高程;λi为第i个参证气象站点所占权重;x0和xi分别为模拟栅格单元和第i个参证气象站点的横坐标;y0和yi分别为模拟栅格单元和第i个参证气象站点的纵坐标;di为模拟栅格单元和第i个参证气象站点的平面距离。
dM/dz为气象数据随高程增加的线性变化率,z和D分别为参证气象站点的高程和数据值;λi为第i个参证气象站点所占权重;di为第i个参证气象站点与待插值栅格点的平面投影距离。
2、采用Green-Ampt入渗模型模拟计算栅格单元的降雨日累计入渗量。当入渗湿润锋到达第m层土壤时,入渗能力f为:
若湿润锋进入第m-1层土壤时地表持续积水,降雨日累计入渗量F由式(5)计算;若前一时段tn-1地表无积水,当前时段tn开始积水,降雨日累计入渗量F则用式(6)计算。
其中,
Fp=Am-1(Ip/km-1)-Bm-1 (10)
tp=tn-1+(Fp-Fn-1)/Ip (11)
Δθi=θs,i-θi,0 (12)
式中,m为湿润锋最终到达的土层数;i表示湿润锋当前所在土层数;km为湿润锋所在m层时的土壤导水系数,单位为mm/day;ki为第i层土壤导水系数,单位为mm/day;Δθi为第i层土壤含水量的变化量,cm3/cm3;Δθm为湿润锋所在m层的土壤含水量变化量,cm3/cm3;Am-1和Bm-1为参数,无量纲;Fm-1为m-1层以上的土壤累积入渗量,mm;FP为当前时段开始积水时的累积入渗量,单位为mm;t为时刻;tn-1为第n-1时段,min;Fn-1为第n-1时段的土壤累积入渗量,mm;tm-1为湿润锋进入m-1和m层土壤界面时的时间;tp为开始积水时刻;Li为第i层土厚,单位为mm;SWm为第m层湿润锋处的毛管吸力,单位为mm;θ表示土壤饱和含水量,cm3/cm3;Ip为积水开始时的雨强,mm/min;θs,i为第i层土壤饱和含水量,cm3/cm3;θi,0为土壤初始含水量,cm3/cm3。
3、利用Penman公式模拟计算栅格单元地表当日最大可蒸发量。计算公式如下:
式中,Δ为饱和水汽压曲线斜率(kPa/℃);RN为净放射量(MJ/m2);G为地中热传导(MJ/m2);ρa为空气密度(kg/m3);CP为空气定压比热(J/kg/℃);δe为实际与饱和水蒸气压之差(kPa);ra为空气动力学阻抗(s/m);λ为水体的气化潜热(MJ/kg);γ为空气湿度常数(kPa/℃);Pa为大气压(kPa)。
4、利用水量平衡方程计算模拟单元日积水深H。计算公式如下:
式中,H为日积水深/mm;H0为当日初始积水深/mm;P为当日降水深/mm;E为当日最大可蒸发量/mm;F为降雨日累计入渗量/mm。
第二步:确定坡面动态汇流路径
1、以日为汇流模拟时段,赋予研究区内每个模拟栅格单元5个基本属性:坐标x、坐标y、地形高程z、地表积水深H、是否已参与汇流计算(True或False,时段开始时初始默认为False);
2、逐汇流模拟时段内,引入当前模拟栅格单元的地表积水深H,将其与模拟栅格单元的地形高程z相加,以此作为坡面漫流、河道汇流和漫溢计算所需的“虚拟高程”,记作DH(图2);
3、遍历模拟栅格单元,寻找DH最大值所在的模拟栅格单元,确定流向和出流量(见“坡面多流向水量分配权重确定和坡面多流向漫流水量确定”),然后更新模拟栅格单元的DH,并将结束汇流的模拟栅格单元标记为True,即表示本时段内该栅格不再参与汇流计算。
4、逐次遍历,直至所有模拟栅格单元的汇流路径和出流量全部计算完成。
第三步:确定坡面多流向水量分配权重确定
1、比较模拟栅格单元和该单元相邻的8个栅格单元的DH大小,若模拟栅格单元DH大于相邻栅格单元,则认为水流流向该临近栅格单元,研究认为模拟栅格单元的水流应向邻域中所有DH较低的方向进行分配,由此对流向进行编码(图3)。
2、确定模拟栅格单元在临近多个DH较低栅格单元上的水量分配权重。计算公式如下:
式中,i和j表示两个不同栅格单元;DHi和DHj分别为栅格单元i和栅格单元j的虚拟高程;wji表示从DH较高栅格单元j到相邻DH较低单元i的分配水量权重;sji表示栅格单元j到栅格单元i间的水力梯度;nji为栅格单元j到栅格单元i间的坡面曼宁糙率系数;ni和nj分别为栅格单元i和栅格单元j的坡面曼宁糙率系数;xi和xj分别为栅格单元i和栅格单元j的横坐标;yi和yj分别为栅格单元i和栅格单元j的的纵坐标。
第四步:确定坡面多流向漫流水量
1、使用运动波模型计算水流流出栅格单元在本汇流时段内的可出流量。计算公式如下:
Sf=S0 (20)
式中,A为过流断面面积,单位为m2;Q为断面流量,单位为m3/s;qL为单宽流入量,单位为m3/s/m;Sf为摩擦比降;S0为地面坡降或河道坡降;R为过流断面水力半径,单位为m;n为坡面曼宁糙率系数;t为时间;x为河段长度;
2、计算水流流出栅格单元与流入栅格单元“虚拟高程”出现持平时的出流量dc,以此作为该计算单元在本汇流时段的最大出流量。计算时已知时段初流出栅格单元和流入栅格单元的积水深,以及流向不同栅格单元的水量分配权重,通过建立水量平衡方程即可得到不同流向上的流入量以及栅格的流出量。以双向流举例如下(图4):
式中,Dc0、D10分别表示汇流时段初流出单元与方向1流入单元的积水深,由于单元面积相等,以积水深作为单元水量;w1、w2分别表示两方向上的水量分配权重;dc表示时段末水流流出栅格单元与流入栅格单元“虚拟高程”出现持平时的出流量;d1、d2分别表示时段末出流方向1和2上的栅格单元入流量。由于汇流时段初积水和水量分配权重已知,可求得dc、d1、d2。
3、取Q与dc两者之间的较小值作为单元出流量最终模拟值,并按水量分配权重w1、w2对各方向单元流入量进行修正,即各方向单元流入量等于单元出流量最终模拟值乘以各方向上水量分配权重。
最后应说明的是,以上仅用以说明本发明的技术方案而非限制,尽管参照较佳布置方案对本发明进行了详细说明,本领域的普通技术人员应当理解,可以对本发明的技术方案(比如各种公式的运用、步骤的先后顺序等)进行修改或者等同替换,而不脱离本发明技术方案的精神和范围。
Claims (5)
1.一种考虑地表积水深变化的坡面多流向漫流分析方法,其特征在于:所述方法包括四个部分:地表积水深确定、坡面动态汇流路径确定、坡面多流向水量分配权重确定和坡面多流向漫流水量确定,各部分的具体步骤如下:
地表积水深确定:
步骤一:获取研究区内气象站点在模拟期内的气象数据,所述气象数据包括逐日降水、气温、风速、日照、湿度,选择参证气象站,利用反距离平方加权插值算法将参证气象站点的气象数据空间展布到研究区模拟栅格单元上;
步骤二:采用Green-Ampt入渗模型模拟计算栅格单元的降雨日累计入渗量;
步骤三:利用Penman公式模拟计算栅格单元地表当日最大可蒸发量;
步骤四:利用水量平衡方程计算模拟单元日积水深,计算公式如下:
式中,H为日积水深/mm;H0为当日初始积水深/mm;P为当日降水深/mm;E为当日最大可蒸发量/mm;F为降雨日累计入渗量/mm;
坡面动态汇流路径确定的具体步骤:
步骤1:以日为汇流模拟时段,赋予研究区内每个模拟栅格单元5个基本属性:坐标x、坐标y、地形高程z、日积水深H、是否已参与汇流计算;已参与汇流计算标记为ture,否则标记为False;
步骤2:逐汇流模拟时段内,引入当前模拟栅格单元日积水深H,将其与模拟栅格单元的地形高程z相加,以此作为坡面漫流、河道汇流和漫溢计算所需的“虚拟高程”,记作DH;
步骤3:遍历模拟栅格单元,寻找DH最大值所在的模拟栅格单元,确定流向和出流量,然后更新模拟栅格单元的DH,并将结束汇流的模拟栅格单元标记为True,即表示本时段内该栅格不再参与汇流计算;
步骤4:逐次遍历,直至所有模拟栅格单元的汇流路径和出流量全部计算完成;
坡面多流向水量分配权重确定的具体步骤:
步骤1):比较模拟栅格单元和该单元相邻的8个栅格单元的DH大小,若模拟栅格单元DH大于相邻栅格单元,则认为水流流向相邻栅格单元,模拟栅格单元的水流应向邻域中所有DH较低的方向进行分配,由此对流向进行编码;
步骤2):结合地表曼宁糙率,确定模拟栅格单元在临近多个DH较低栅格单元上的水量分配权重;
坡面多流向漫流水量确定的具体步骤:
步骤Ⅰ:使用运动波模型计算水流流出栅格单元在本汇流时段内的流量,记为可出流量;计算公式如下:
Sf=S0 (20)
式中,A为过流断面面积,单位为m2;Q为断面流量,单位为m3/s,即为可出流量;qL为单宽流入量,单位为m3/s/m;Sf为摩擦比降;S0为地面坡降或河道坡降;R为过流断面水力半径,单位为m;n为坡面曼宁糙率系数;t为时间;x为河段长度;
步骤Ⅱ:计算水流流出栅格单元与流入栅格单元“虚拟高程”出现持平时的出流量dc,以此作为所述水流流出栅格单元在本汇流时段的最大出流量;
步骤Ⅲ:取Q与dc两者之间的较小值作为单元出流量最终模拟值,并按水量分配权重对各方向单元流入量进行修正。
2.根据权利要求1所述的考虑地表积水深变化的坡面多流向漫流分析方法,其特征在于:步骤一中利用反距离平方加权插值算法将参证气象站点的气象数据空间展布到研究区模拟栅格单元上,计算公式如下:
式中,z和M分别表示气象站的高程和气象数据值,dM/dz为气象数据随高程增加的线性变化率;m表示参证气象站点的总数;D和Di分别为模拟栅格单元和第i个参证气象站点的气象数据值;z0和zi分别为模拟栅格单元和第i个参证气象站点的高程;λi为第i个参证气象站点所占权重;x0和xi分别为模拟栅格单元和第i个参证气象站点的横坐标;y0和yi分别为模拟栅格单元和第i个参证气象站点的纵坐标;di为模拟栅格单元和第i个参证气象站点的平面距离。
3.根据权利要求1所述的考虑地表积水深变化的坡面多流向漫流分析方法,其特征在于:步骤二中采用Green-Ampt入渗模型模拟计算栅格单元的降雨日累计入渗量,计算方法如下:
当入渗湿润锋到达第m层土壤时,入渗能力f为:
若湿润锋进入第m-1层土壤时地表持续积水,降雨日累计入渗量F由式(5)计算;若前一时段tn-1地表无积水,当前时段tn开始积水,降雨日累计入渗量F则用式(6)计算;
其中,
Fp=Am-1(Ip/km-1)-Bm-1 (10)
tp=tn-1+(Fp-Fn-1)/Ip (11)
Δθi=θs,i-θi,0 (12)
式中,m为湿润锋最终到达的土层数;i表示湿润锋当前所在土层数;km为湿润锋所在m层时的土壤导水系数,单位为mm/day;ki为第i层土壤导水系数,单位为mm/day;Δθi为第i层土壤含水量的变化量,cm3/cm3;Δθm为湿润锋所在m层的土壤含水量变化量,cm3/cm3;Am-1和Bm-1为参数,无量纲;Fm-1为m-1层以上的土壤累积入渗量,mm;FP为当前时段开始积水时的累积入渗量,单位为mm;t为时刻;tn-1为第n-1时段,min;Fn-1为第n-1时段的土壤累积入渗量,mm;tm-1为湿润锋进入m-1和m层土壤界面时的时间;tp为开始积水时刻;Li为第i层土厚,单位为mm;SWm为第m层湿润锋处的毛管吸力,单位为mm;θ表示土壤饱和含水量,cm3/cm3;Ip为积水开始时的雨强,mm/min;θs,i为第i层土壤饱和含水量,cm3/cm3;θi,0为土壤初始含水量,cm3/cm3。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010885571.3A CN111931440B (zh) | 2020-08-28 | 2020-08-28 | 一种考虑地表积水深变化的坡面多流向漫流分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010885571.3A CN111931440B (zh) | 2020-08-28 | 2020-08-28 | 一种考虑地表积水深变化的坡面多流向漫流分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111931440A CN111931440A (zh) | 2020-11-13 |
CN111931440B true CN111931440B (zh) | 2021-03-05 |
Family
ID=73308556
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010885571.3A Expired - Fee Related CN111931440B (zh) | 2020-08-28 | 2020-08-28 | 一种考虑地表积水深变化的坡面多流向漫流分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111931440B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116821564B (zh) * | 2023-08-28 | 2024-01-09 | 中建八局第三建设有限公司 | 面向积水区沙化处理的添加剂质量计算系统、方法及存储介质 |
Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105760666A (zh) * | 2016-02-04 | 2016-07-13 | 南京信息工程大学 | 一种基于dem的汇水区临界雨量计算方法 |
CN106599605A (zh) * | 2017-02-22 | 2017-04-26 | 中国水利水电科学研究院 | 一种石灰岩土石山区山坡尺度水文过程模拟方法 |
CN106599473A (zh) * | 2016-12-15 | 2017-04-26 | 中国水利水电科学研究院 | 一种坡面水沙过程耦合模拟方法 |
CN109492259A (zh) * | 2018-10-15 | 2019-03-19 | 华北水利水电大学 | 一种城市水文模拟系统 |
CN109871621A (zh) * | 2019-02-25 | 2019-06-11 | 中国水利水电科学研究院 | 城市暴雨内涝汇水区分析方法 |
CN110188476A (zh) * | 2019-05-31 | 2019-08-30 | 青海大学 | 一种基于分层土壤的水沙过程计算方法 |
CN110298076A (zh) * | 2019-05-27 | 2019-10-01 | 广州奥格智能科技有限公司 | 一种基于gis和swmm的城市内涝智能建模及分析方法 |
CN110442988A (zh) * | 2019-08-12 | 2019-11-12 | 水利部交通运输部国家能源局南京水利科学研究院 | 一种基于元胞自动机的城市地表径流流向计算方法及装置 |
CN110717232A (zh) * | 2019-09-11 | 2020-01-21 | 中国水利水电科学研究院 | 一种考虑梯田淤地坝影响的子流域汇流模拟方法 |
CN110955924A (zh) * | 2019-09-11 | 2020-04-03 | 中国水利水电科学研究院 | 一种考虑淤地坝影响的子流域汇流模拟方法 |
CN111090902A (zh) * | 2019-12-12 | 2020-05-01 | 中国水利水电科学研究院 | 一种基于地下水模型的坎儿井数值模拟方法 |
CN111364379A (zh) * | 2020-03-25 | 2020-07-03 | 宗秀伟 | 一种隧道路面警示器 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106777688B (zh) * | 2016-12-15 | 2018-05-04 | 中国水利水电科学研究院 | 一种水土保持工程减水减沙效益定量评价方法 |
CN106599471B (zh) * | 2016-12-15 | 2018-07-20 | 中国水利水电科学研究院 | 一种分布式流域水沙过程耦合模拟方法 |
WO2019064037A1 (en) * | 2017-09-26 | 2019-04-04 | Total Sa | METHOD FOR DEFINING LOCATIONS OF A PLURALITY OF WELLS IN A FIELD, ASSOCIATED SYSTEM AND COMPUTER PROGRAM PRODUCT |
CN108108556B (zh) * | 2017-12-20 | 2019-03-01 | 中国水利水电科学研究院 | 一种基于耗散-汇合结构的灌区水循环模型构建方法 |
-
2020
- 2020-08-28 CN CN202010885571.3A patent/CN111931440B/zh not_active Expired - Fee Related
Patent Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105760666A (zh) * | 2016-02-04 | 2016-07-13 | 南京信息工程大学 | 一种基于dem的汇水区临界雨量计算方法 |
CN106599473A (zh) * | 2016-12-15 | 2017-04-26 | 中国水利水电科学研究院 | 一种坡面水沙过程耦合模拟方法 |
CN106599605A (zh) * | 2017-02-22 | 2017-04-26 | 中国水利水电科学研究院 | 一种石灰岩土石山区山坡尺度水文过程模拟方法 |
CN109492259A (zh) * | 2018-10-15 | 2019-03-19 | 华北水利水电大学 | 一种城市水文模拟系统 |
CN109871621A (zh) * | 2019-02-25 | 2019-06-11 | 中国水利水电科学研究院 | 城市暴雨内涝汇水区分析方法 |
CN110298076A (zh) * | 2019-05-27 | 2019-10-01 | 广州奥格智能科技有限公司 | 一种基于gis和swmm的城市内涝智能建模及分析方法 |
CN110188476A (zh) * | 2019-05-31 | 2019-08-30 | 青海大学 | 一种基于分层土壤的水沙过程计算方法 |
CN110442988A (zh) * | 2019-08-12 | 2019-11-12 | 水利部交通运输部国家能源局南京水利科学研究院 | 一种基于元胞自动机的城市地表径流流向计算方法及装置 |
CN110717232A (zh) * | 2019-09-11 | 2020-01-21 | 中国水利水电科学研究院 | 一种考虑梯田淤地坝影响的子流域汇流模拟方法 |
CN110955924A (zh) * | 2019-09-11 | 2020-04-03 | 中国水利水电科学研究院 | 一种考虑淤地坝影响的子流域汇流模拟方法 |
CN111090902A (zh) * | 2019-12-12 | 2020-05-01 | 中国水利水电科学研究院 | 一种基于地下水模型的坎儿井数值模拟方法 |
CN111364379A (zh) * | 2020-03-25 | 2020-07-03 | 宗秀伟 | 一种隧道路面警示器 |
Non-Patent Citations (5)
Title |
---|
"Development and validation of a physically-based, national-scale hydrological model in China";Huan Liu et.al;《Journal of Hydrology》;20200825;论文第1-14页 * |
"The treatment of flat areas and depressions in automated drainage analysis of raster digital elevation models:";Lawrence W. Martz and Jurgen Garbrecht;《Hydrological Processes》;19981231;第843-855页 * |
"基于多流向流量分配的水流阻力算法试验研究";张升堂等;《中国农村水利水电》;20191231(第10期);第153-157页 * |
"无填洼格网DEM数字流域分析算法研究";徐昕;《中国优秀硕士学位论文全文数据库 基础科学辑》;20180315(第03期);论文第2-4章 * |
"考虑产流机制空间差异的全国尺度分布式水文模型构建及应用";刘欢;《中国博士学位论文全文数据库 基础科学辑》;20190815(第08期);论文第1-6章 * |
Also Published As
Publication number | Publication date |
---|---|
CN111931440A (zh) | 2020-11-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109492299B (zh) | 基于swmm与modflow耦合的水资源模拟方法 | |
CN112036093A (zh) | 一种陆面水文耦合模型 | |
Jia et al. | Development of the WEP-L distributed hydrological model and dynamic assessment of water resources in the Yellow River basin | |
CN105678067B (zh) | 城市河涌防洪排涝动态预报控制方法及系统 | |
CN107045568B (zh) | 基于动态规划逐次逼近法的河道糙率反演方法 | |
CN113011685A (zh) | 一种无径流资料地区内陆湖泊水位变化模拟预测方法 | |
CN115859676B (zh) | 一种顾及气候要素的多层次城市内涝耦合模拟方法 | |
CN113610264A (zh) | 精细化电网台风洪涝灾害预测模型 | |
CN111723505B (zh) | 一种流域水质水量监测系统 | |
CN112270115B (zh) | 一种基于元胞自动机的复杂地形洪水淹没进程模拟方法 | |
CN115115262A (zh) | 一种洪水风险灾害评估的方法 | |
CN111931440B (zh) | 一种考虑地表积水深变化的坡面多流向漫流分析方法 | |
CN114399103A (zh) | 一种基于cnn的陆水一体化河流水质时空连续预测方法 | |
Gusev et al. | Northern Dvina runoff simulation using land-surface model SWAP and global databases | |
CN117494616B (zh) | 一种闸泵群调度模拟方法及系统 | |
CN112257286B (zh) | 一种多年冻土区温度主导的变源产流模式模拟方法 | |
CN111680886A (zh) | 一种内涝风险预测方法及其系统 | |
CN111783323A (zh) | 一种陆面水文耦合模拟系统及方法 | |
Jhorar | Estimation of effective soil hydraulic parameters for water management studies in semi-arid zones: integral use of modelling, remote sensing and parameter estimation | |
Guo et al. | A study of two mass-consistent models: problems and possible solutions | |
Gusev et al. | Physically based modeling of the long-term dynamics of water balance and snow water storage components in the Ob–Irtysh river basin | |
Al-Aboodi et al. | Estimation of groundwater recharge in Safwan-Zubair area, South of Iraq, using water balance and inverse modeling methods | |
Melas et al. | The impact of sea breeze on air quality in Athens area | |
CN118036893B (zh) | 一种三层结构嵌套的月尺度径流预测方法 | |
CN115062559B (zh) | 可嵌套水文模型的岩溶含水系统蓄水库的定量表征方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20210305 Termination date: 20210828 |
|
CF01 | Termination of patent right due to non-payment of annual fee |