CN108229096B - 一种湿润地区土壤分层产流计算方法 - Google Patents
一种湿润地区土壤分层产流计算方法 Download PDFInfo
- Publication number
- CN108229096B CN108229096B CN201810204046.3A CN201810204046A CN108229096B CN 108229096 B CN108229096 B CN 108229096B CN 201810204046 A CN201810204046 A CN 201810204046A CN 108229096 B CN108229096 B CN 108229096B
- Authority
- CN
- China
- Prior art keywords
- soil
- runoff
- calculating
- water storage
- layer
- 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.)
- Active
Links
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16Z—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
- G16Z99/00—Subject matter not provided for in other main groups of this subclass
Abstract
本发明公开了一种湿润地区土壤分层产流计算方法,通过在流域不同地点及土壤不同分层上设置蓄水箱,提出一种垂向“水箱并联”、横向“水箱串联”的蓄满产流分水源计算方法,能够反映流域土壤实际的分层产流过程,解决了湿润地区蓄满产流模式中分水源计算缺乏明确物理意义的难题,具有较强的工程意义。
Description
技术领域
本发明涉及一种湿润地区土壤分层产流计算方法。
背景技术
新安江模型是我国赵人俊教授研制,在国内外湿润地区水文预报中应用广泛的概念性流域水文模型。新安江模型中产流计算采用蓄满产流理论,采用自由水蓄水库结构将总径流划分为地表径流、壤中流和地下径流三种。但是,模型中的自由水蓄水库是“虚拟的”,在实际土壤中难以测定其水分状态,也不能描述土壤不同深度壤中流运动的差异。为此,需要构建一种能够描述土壤不同深度水分运动特点、且在实际土壤中可测定的自由水蓄水库产流结构。
目前,概念性水文模型主要通过率定出流系数进行壤中流计算,因而不能反映土壤垂向及流域空间上的壤中流产流不均匀性。日本的水箱模型采用串联水箱模拟各个土层的产流,但其物理概念不是直接的,没有明确的土壤含水量概念,它的水箱结构是单纯依赖模型率定的,不能反映流域土壤实际的分层产流过程。
发明内容
本发明的目的在于克服现有技术中的不足,提供一种湿润地区土壤分层产流计算方法,解决现有技术中壤中流计算过度概化、不能描述壤中流产流在土壤不同深度及流域空间差异性的技术问题。
为解决上述技术问题,本发明所采用的技术方案是:一种湿润地区土壤分层产流计算方法,包括如下步骤:
将研究流域划分为若干个正方形网格;
将土壤垂向分为N层,N≥3,其中:1~N-1土壤分层用于模拟不同出流速度的壤中流径流成分,第N层土壤分层用于模拟地下径流;
在各个网格土壤分层上设置蓄水箱,形成同一土壤分层上不同网格间“水箱串联”及同一网格上不同土壤分层间“水箱并联”的非线性产流结构;
假定蒸散发只发生在上3层土壤中,计算各网格上的净雨量PE;
根据各网格不同土壤分层蓄水箱侧向出流系数KIi、底孔下渗系数KGi以及各土壤分层的初始土壤含水量Wi,计算地下径流RG、1~N-1土壤分层的壤中流RIi以及地面径流RS,计算公式分别如下:
RG=KIN·WN (1-1),
RIi=KIi·Wi (1-2),
式中:i表示土壤分层序号,计算地下径流时i=N,计算壤中流时i∈[1,N-1];WM表示研究流域平均蓄水容量;
根据水量平衡公式,计算各土壤分层的各网格上产流总量R,公式如下:
进一步的,所述非线性产流结构形成的具体方法是:
在每一层蓄水箱旁侧设置出流孔,底部设置下渗孔,最下层蓄水箱设置出流孔但不设置下渗孔。
进一步的,土壤分层数目的确定方法如下:
在研究流域内选择若干代表性实验点,记录分析土壤纵剖面分层的实际情况,在室内实验室测定不同分层土壤的机械组成,进而确定土壤分层数目N。
进一步的,所述净雨量PE的计算方法如下:
根据新安江模型三层蒸发模式计算上3层土壤的实际蒸发量,分别为EU、EL和ED,则实际蒸发总量为E=EU+EL+ED;
各网格上降雨P扣除实际蒸发总量即为净雨量,净雨量PE=P-E。
进一步的,研究流域平均蓄水容量WM的计算方法如下:
进一步的,各网格不同土壤分层蓄水箱侧向出流系数KIi、底孔下渗系数KGi根据流域土壤质地分布和各土层的水力传导度资料确定,具体计算公式如下:
KIi=α·κi,i=1,2,…,N (1-5)
KGi=β·(κi+κi+1)/2,i=1,2,…,N-1 (1-6)
式中:κi为第i层土壤的饱和水力传导度,α和β分别为反映土壤侧向和纵向与孔隙率相关的修正系数。
与现有技术相比,本发明所达到的有益效果是:
通过在流域不同地点及土壤不同分层上设置蓄水箱,提出一种垂向“水箱并联”、横向“水箱串联”的蓄满产流分水源计算方法,能够反映流域土壤实际的分层产流过程,解决了湿润地区蓄满产流模式中分水源计算缺乏明确物理意义的难题,具有较强的工程意义。
具体实施方式
本发明通过在流域不同地点及土壤不同分层上设置蓄水箱,提出一种垂向“水箱并联”、横向“水箱串联”的蓄满产流分水源计算方法,能够反映流域土壤实际的分层产流过程,解决了湿润地区蓄满产流模式中分水源计算缺乏明确物理意义的难题,具有较强的工程意义。
下面结合实施例对本发明作进一步描述。以下实施例仅用于更加清楚地说明本发明的技术方案,而不能以此来限制本发明的保护范围。
现有某湿润地区一小流域的降雨、蒸发皿蒸发等气象水文资料,依据本发明方法,该流域研究期内产流的计算过程为:
(1)将研究流域划分为若干个1km×1km的正方形网格,将降水、蒸发资料空间插值至每个网格,以网格为基本产流计算单元;
(2)在研究流域内选择若干代表性实验点,记录分析土壤纵剖面分层的实际情况,在室内实验室测定不同分层土壤的机械组成,进而确定土壤分层数目N(一般地,N≥3),根据各土层厚度和土壤孔隙率资料,确定各土壤分层的蓄水容量Wsmi,i=1,2,…,N,则流域平均蓄水容量为
(3)在各网格土壤分层上设置蓄水箱,每一层蓄水箱旁侧设有出流孔,底部设有下渗孔,最下层蓄水箱仅有出流孔无下渗孔,因此在流域上形成同一土壤分层上不同网格间“水箱串联”及同一网格上不同土壤分层间“水箱并联”的非线性产流结构;
(4)根据网格蒸发皿蒸发EM和折算系数k,计算得到流域蒸散发能力Ep=k·EM,假定蒸散发只发生在上3层土壤中,根据新安江模型三层蒸发模式计算上3层土壤的实际蒸发量,分别为EU、EL和ED,则实际蒸发总量为E=EU+EL+ED;
(5)网格上降雨P扣除实际蒸发总量,即得到净雨量PE=P-E;
(6)根据流域土壤质地分布和各土层的水力传导度资料,确定各网格不同土壤分层蓄水箱侧向出流系数KIi=α·κi(i=1,2,…,N)和底孔下渗系数KGi=β·(κi+κi+1)/2(i=1,2,…,N-1),其中,κi为第i层土壤的饱和水力传导度,α和β为参数;
(7)根据各土层的初始土壤含水量Wi(i=1,2,…,N),计算地下径流为RG=KIN·WN,1~N-1土壤分层的壤中流为RIi=KIi·Wi,地面径流为:
(8)根据水量平衡公式,采用新安江模型蓄满产流模式计算各网格上产流总量R,公式如下:
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明技术原理的前提下,还可以做出若干改进和变形,这些改进和变形也应视为本发明的保护范围。
Claims (4)
1.一种湿润地区土壤分层产流计算方法,其特征在于,包括如下步骤:
将研究流域划分为若干个正方形网格;
将土壤垂向分为N层,N≥3,其中:1~N-1土壤分层用于模拟不同出流速度的壤中流径流成分,第N层土壤分层用于模拟地下径流;
在各个网格土壤分层上设置蓄水箱,形成同一土壤分层上不同网格间“水箱串联”及同一网格上不同土壤分层间“水箱并联”的非线性产流结构;
假定蒸散发只发生在上3层土壤中,计算各网格上的净雨量PE;
根据各网格不同土壤分层蓄水箱侧向出流系数KIi、底孔下渗系数KGi以及各土壤分层的初始土壤含水量Wi,计算地下径流RG、1~N-1土壤分层的壤中流RIi以及地面径流RS,计算公式分别如下:
RG=KIN·WN (1-1),
RIi=KIi·Wi (1-2),
式中:i表示土壤分层序号,计算地下径流时i=N,计算壤中流时i∈[1,N-1];WM表示研究流域平均蓄水容量;
根据水量平衡公式,计算各土壤分层的各网格上产流总量R,公式如下:
研究流域平均蓄水容量WM的计算方法如下:
各网格不同土壤分层蓄水箱侧向出流系数KIi、底孔下渗系数KGi根据流域土壤质地分布和各土层的水力传导度资料确定,具体计算公式如下:
KIi=α·κi,i=1,2,…,N (1-5)
KGi=β·(κi+κi+1)/2,i=1,2,…,N-1 (1-6)
式中:κi为第i层土壤的饱和水力传导度,α和β分别为反映土壤侧向和纵向与孔隙率相关的修正系数。
2.根据权利要求1所述的湿润地区土壤分层产流计算方法,其特征在于,所述非线性产流结构形成的具体方法是:
在每一层蓄水箱旁侧设置出流孔,底部设置下渗孔,最下层蓄水箱设置出流孔但不设置下渗孔。
3.根据权利要求1所述的湿润地区土壤分层产流计算方法,其特征在于,土壤分层数目的确定方法如下:
在研究流域内选择若干代表性实验点,记录分析土壤纵剖面分层的实际情况,在室内实验室测定不同分层土壤的机械组成,进而确定土壤分层数目N。
4.根据权利要求1所述的湿润地区土壤分层产流计算方法,其特征在于,所述净雨量PE的计算方法如下:
根据新安江模型三层蒸发模式计算上3层土壤的实际蒸发量,分别为EU、EL和ED,则实际蒸发总量为E=EU+EL+ED;
各网格上降雨P扣除实际蒸发总量即为净雨量,净雨量PE=P-E。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810204046.3A CN108229096B (zh) | 2018-03-13 | 2018-03-13 | 一种湿润地区土壤分层产流计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810204046.3A CN108229096B (zh) | 2018-03-13 | 2018-03-13 | 一种湿润地区土壤分层产流计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108229096A CN108229096A (zh) | 2018-06-29 |
CN108229096B true CN108229096B (zh) | 2020-11-10 |
Family
ID=62659382
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810204046.3A Active CN108229096B (zh) | 2018-03-13 | 2018-03-13 | 一种湿润地区土壤分层产流计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108229096B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114491768B (zh) * | 2022-02-16 | 2022-11-18 | 河海大学 | 一种变动产流层统一模式产流计算方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20140031754A (ko) * | 2012-09-05 | 2014-03-13 | 경희대학교 산학협력단 | 취약성 평가 지표를 이용한 하수처리장 취약성 경감 시스템 |
CN104392111A (zh) * | 2014-11-10 | 2015-03-04 | 连承锴 | 基于水位样本的洪水预报方法 |
CN104846926A (zh) * | 2015-05-19 | 2015-08-19 | 浙江大学 | 一种体积法测量管网外源入侵的系统和方法 |
CN105912836A (zh) * | 2016-03-25 | 2016-08-31 | 环境保护部卫星环境应用中心 | 一种纯遥感数据驱动的流域水循环模拟方法 |
CN106548012A (zh) * | 2016-10-17 | 2017-03-29 | 河海大学 | 一种通过流量过程线估算流域平均自由水蓄水容量的方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050071139A1 (en) * | 2003-09-29 | 2005-03-31 | Patwardhan Avinash S. | Method and system for water flow analysis |
CN106951612B (zh) * | 2017-03-06 | 2020-07-14 | 河海大学 | 土壤冻融过程中动态蓄水容量产流计算方法 |
-
2018
- 2018-03-13 CN CN201810204046.3A patent/CN108229096B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20140031754A (ko) * | 2012-09-05 | 2014-03-13 | 경희대학교 산학협력단 | 취약성 평가 지표를 이용한 하수처리장 취약성 경감 시스템 |
CN104392111A (zh) * | 2014-11-10 | 2015-03-04 | 连承锴 | 基于水位样本的洪水预报方法 |
CN104846926A (zh) * | 2015-05-19 | 2015-08-19 | 浙江大学 | 一种体积法测量管网外源入侵的系统和方法 |
CN105912836A (zh) * | 2016-03-25 | 2016-08-31 | 环境保护部卫星环境应用中心 | 一种纯遥感数据驱动的流域水循环模拟方法 |
CN106548012A (zh) * | 2016-10-17 | 2017-03-29 | 河海大学 | 一种通过流量过程线估算流域平均自由水蓄水容量的方法 |
Non-Patent Citations (2)
Title |
---|
A statistically-based runoff-yield model;Zhong Minliang et al;《Proceedings of Chinese PUB International Symposium》;20081231;全文 * |
土壤不同含碳量情况下水分特征曲线的拟合及修正;王晓蕾;《广东水利水电》;20130731(第7期);23-27 * |
Also Published As
Publication number | Publication date |
---|---|
CN108229096A (zh) | 2018-06-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109492299B (zh) | 基于swmm与modflow耦合的水资源模拟方法 | |
CN111914432B (zh) | 一种基于大数据的水文预报方法 | |
CN113610264B (zh) | 一种精细化电网台风洪涝灾害预测系统 | |
Loiola et al. | Hydrological performance of modular-tray green roof systems for increasing the resilience of mega-cities to climate change | |
CN106202790B (zh) | 一种分布式河北模型构建方法及其应用 | |
Ahmad et al. | Sustainable use of groundwater for irrigation: a numerical analysis of the subsoil water fluxes | |
Tuong et al. | Quantifying flow processes during land soaking of cracked rice soils | |
CN112084671B (zh) | 城市时变增益降雨-径流过程模拟计算方法 | |
CN109614655B (zh) | 一种河流径流量的研究方法 | |
CN110188323B (zh) | 一种适宜于城市水生态保护的水面率值域测算方法 | |
CN113011685A (zh) | 一种无径流资料地区内陆湖泊水位变化模拟预测方法 | |
Zhao et al. | Application of a simple raster-based hydrological model for streamflow prediction in a humid catchment with polder systems | |
CN108182539A (zh) | 一种基于生态需水计算的河道治理需水量确定方法 | |
CN108229096B (zh) | 一种湿润地区土壤分层产流计算方法 | |
Zhao et al. | Rainwater harvesting capacity of soils subjected to reservoir tillage during rainfall on the Loess Plateau of China | |
Selbig et al. | Evaluation of turf-grass and prairie-vegetated rain gardens in a clay and sand soil, Madison, Wisconsin, Water years 2004-08 | |
CN110287595B (zh) | 一种城市不同下垫面减灾效应分析方法 | |
van den Elsen et al. | Discharge and sediment measurements at the outlet of a watershed on the Loess plateau of China | |
CN103116001B (zh) | 干旱半干旱高山区优先流观测方法 | |
CN113888710A (zh) | 一种基于dem的动力学流域产沙量模拟方法 | |
CN112663760A (zh) | 基于lid措施的城市雨水径流调控方法 | |
CN102590060B (zh) | 不同坡度土壤稳定入渗仪 | |
CN111982780B (zh) | 一种城市新区规划下垫面与旧有下垫面的关系测算方法 | |
CN112561205B (zh) | 一种考虑淤地坝和梯田对径流影响的定量测量方法 | |
Daniels et al. | Soil variability and productivity: future developments |
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 |