CN113657047B - 一种基于双自由水库的产流计算方法、装置及存储介质 - Google Patents

一种基于双自由水库的产流计算方法、装置及存储介质 Download PDF

Info

Publication number
CN113657047B
CN113657047B CN202110930898.2A CN202110930898A CN113657047B CN 113657047 B CN113657047 B CN 113657047B CN 202110930898 A CN202110930898 A CN 202110930898A CN 113657047 B CN113657047 B CN 113657047B
Authority
CN
China
Prior art keywords
runoff
aeration zone
reservoir
free reservoir
period
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
Application number
CN202110930898.2A
Other languages
English (en)
Other versions
CN113657047A (zh
Inventor
陈新宇
张珂
张企诺
黄鹏年
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hohai University HHU
Original Assignee
Hohai University HHU
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Hohai University HHU filed Critical Hohai University HHU
Priority to CN202110930898.2A priority Critical patent/CN113657047B/zh
Publication of CN113657047A publication Critical patent/CN113657047A/zh
Application granted granted Critical
Publication of CN113657047B publication Critical patent/CN113657047B/zh
Priority to US17/877,164 priority patent/US20230059267A1/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V20/00—Geomodelling in general
    • 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
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • E—FIXED CONSTRUCTIONS
    • E03—WATER SUPPLY; SEWERAGE
    • E03F—SEWERS; CESSPOOLS
    • E03F1/00—Methods, systems, or installations for draining-off sewage or storm water
    • E03F1/006—Pneumatic sewage disposal systems; accessories specially adapted therefore
    • E03F1/007—Pneumatic sewage disposal systems; accessories specially adapted therefore for public or main systems
    • 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
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
    • Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Algebra (AREA)
  • General Engineering & Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Computational Mathematics (AREA)
  • Geometry (AREA)
  • Evolutionary Computation (AREA)
  • Computer Hardware Design (AREA)
  • Fluid Mechanics (AREA)
  • Computing Systems (AREA)
  • Operations Research (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geophysics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了一种基于双自由水库的产流计算方法、装置及存储介质，计算方法包括：在新安江模型的三层包气带结构的基础上设置张力水蓄满层，形成包括上层包气带、下层包气带、深层包气带和张力水蓄满层的四层包气带结构；建立双自由水库产流结构：将包气带内自由水能够占据的空间划分为上层自由水库和下层自由水库，其中上层包气带、下层包气带和深层包气带占据上层自由水库，张力水蓄满层占据下层自由水库；使用蓄满产流算法，计算时段产流量；根据建立的双自由水库产流结构将时段产流量划分为地表径流、壤中流和地下径流三种径流成分。本发明提出了基于双自由水库的产流结构，可以很好地适用于浅层地下水埋深较深的半干旱半湿润流域。

Description

一种基于双自由水库的产流计算方法、装置及存储介质
技术领域
本发明属于水文预报领域，具体涉及一种基于双自由水库的产流计算方法。
背景技术
水文模型是模拟流域水文过程的重要方法，一般主要由产流计算和汇流计算组成。产流计算目前主要有两种算法：一种是蓄满产流算法；另一种是超渗产流算法。使用蓄满产流算法的水文模型典型代表是新安江模型。新安江模型中的蓄满产流算法包括产流量计算和分水源计算两部分。目前新安江模型在国内的湿润区流域得到了广泛的应用，并取得了良好的模拟效果。在半干旱半湿润区，新安江模型也有较多的应用，但模拟效果不如在湿润区的模拟效果好。出现这种现象的一个主要原因便是，由于气候条件及人类对地下水的开采，半干旱半湿润区的浅层地下水埋深较深，包气带较厚，这导致了半干旱半湿润区的产流规律与湿润区不同。基于湿润区的产流规律构建的产流计算方法不适用于半干旱半湿润区的实际情况，这使得新安江模型在半干旱半湿润区的模拟效果较差。因此，需要提出一种新的产流计算方法来提高半干旱半湿润区的水文模拟效果。
发明内容
本发明的目的是提出一种更好地模拟半干旱半湿润区流域的产流过程的基于双自由水库的产流计算方法。
为了实现上述目的，本发明具体采用以下技术方案：
一种基于双自由水库的产流计算方法，其特征在于，包括以下步骤：
建立双自由水库产流结构：在新安江模型的三层包气带结构的基础上，在深层包气带之下设置张力水蓄满层，形成包括上层包气带、下层包气带、深层包气带和张力水蓄满层的四层包气带结构；将包气带内自由水能够占据的空间划分为上层自由水库和下层自由水库，其中上层包气带、下层包气带和深层包气带占据上层自由水库，张力水蓄满层占据下层自由水库；
使用蓄满产流算法，计算时段产流量；
根据建立的双自由水库产流结构将时段产流量划分为地表径流、壤中流和地下径流三种径流成分。
使用蓄满产流算法，计算一次降雨过程各时段产流量：
Figure BDA0003210652480000021
式中，R(t)为第t时段的时段产流量，W(t)为第t时段的初始时刻的张力水蓄水量，WM为张力水蓄水容量，Pε(t)为扣除蒸散发损耗、植被冠层截留损耗后的第t时段的时段净降雨量。
根据建立的双自由水库产流结构分别计算地表径流、壤中流和地下径流三种径流的方法是：
Figure BDA0003210652480000022
式中，Rs(t)为第t时段的地表径流，R(t)为第t时段的时段产流量，S(t)为第t时段初的上层自由水库蓄水量，SM为上层自由水蓄水容量；Ri(t)和S(t)的计算方法是：
Ri(t)＝Ki*(R(t)+S(t)-Rs(t)-Fd(t)) (3)
Figure BDA0003210652480000023
式中，Ri(t)为第t时段的壤中流，Ki为壤中流出流系数；Fd(t)为第t时段的下层自由水库入流；S(t-1)为第t-1时段的上层自由水库蓄水量；R(t-1)为第t-1时段的时段产流量；Ri(t-1)为第t-1时段的壤中流；Rs(t-1)为第t-1时段的地表径流；Fd(t-1)为第t-1时段的下层自由水库入流；S(0)为初始时刻上层自由水库蓄水量，根据流域的初始状态观测值或估计值来设定；
Figure BDA0003210652480000024
式中，K为土壤饱和水力传导度，Ψ为湿润锋处土壤吸力，△θ为土壤饱和含水率与田间持水量之差，F(t)为第t时段初的累计渗漏量；
Figure BDA0003210652480000025
式中，Fd(i)为第i时段的下层自由水库入流，F0为第1时段初的累计渗漏量，可设置为某一极小值，如0.001。
Figure BDA0003210652480000031
式中，Rg(t)为第t时段的地下径流，Sl(t)为第t时段初的下层自由水库蓄水量，SLM为下层自由水蓄水容量，Kg为地下径流出流系数。
第t时段初的下层自由水库蓄水量Sl(t)、下层自由水蓄水容量SLM的计算方法为：
Figure BDA0003210652480000032
Figure BDA0003210652480000033
式中，Sl(0)为初始时刻的下层自由水库蓄水量，Zr为初始时刻的河道水位高程，Zg为初始时刻的地下水位高程，Zi为河道下伏含水层底边界高程，μ为地下水位变动带给水度；以上初始时刻状态变量值根据流域的初始状态观测值或估计值来设定。
本发明的有益效果：本发明提供一种基于双自由水库的产流计算方法，在新安江模型的三层包气带结构的基础上，在深层包气带之下设置张力水蓄满层，形成包括上层包气带、下层包气带、深层包气带和张力水蓄满层的四层包气带结构；将包气带内自由水能够占据的空间划分为上、下两层自由水库，其中上层、下层、深层包气带占据上层自由水库，张力水蓄满层占据下层自由水库，构成双自由水库的产流结构；计算方法为，首先，获得扣除蒸散发损耗、植被冠层截留损耗后的时段降雨量，基于蓄满产流理论计算时段产流量；接着，基于双自由水库产流结构，将时段产流量划分为地表径流、壤中流、地下径流三种径流成分。本发明提出了基于双自由水库的产流结构，可以很好地适用于浅层地下水埋深较深的半干旱半湿润区的产流计算，提高了半干旱半湿润区的水文模拟效果。
附图说明
图1是本发明提供的基于双自由水库的产流计算方法流程示意图；
图2是本发明提供的基于双自由水库的产流结构示意图；
图3是本发明提供的基于双自由水库的产流计算方法与原始的基于单自由水库的产流计算方法的洪水模拟效果对比图。
具体实施方式
下面结合附图和具体实施例对本发明作进一步描述。
应当理解的是，此处所描述的具体实施方式仅仅用以解释本发明，并不用于限定本发明。
实施例1
如图1所示，本发明实施例提供了一种基于双自由水库的产流计算方法，包括以下步骤：
步骤1，建立双自由水库产流结构：
基于双自由水库的产流结构的示意图如图2所示。是在新安江模型的三层包气带结构的基础上，在深层包气带之下设置张力水蓄满层，形成包括上层包气带、下层包气带、深层包气带和张力水蓄满层的四层包气带结构；将包气带内自由水能够占据的空间划分为上层自由水库和下层自由水库，其中上层包气带、下层包气带和深层包气带占据上层自由水库，张力水蓄满层占据下层自由水库；
步骤2，使用蓄满产流算法，计算时段产流量，获得扣除蒸散发损耗、植被冠层截留损耗后的时段降雨量：
Figure BDA0003210652480000041
式中，R(t)为第t时段的时段产流量，W(t)为第t时段初始时刻的张力水蓄水量，WM为张力水蓄水容量，Pε(t)为扣除蒸散发损耗、植被冠层截留损耗后的第t时段的时段净降雨量。
步骤3，根据建立的双自由水库产流结构将时段产流量划分为地表径流、壤中流和地下径流三种径流成分：
Figure BDA0003210652480000042
式中，Rs(t)为第t时段的地表径流，R(t)为第t时段的时段产流量，S(t)为第t时段的自由水量，SM为上层自由水蓄水容量；Ri(t)和S(t)的计算方法是：
Ri(t)＝Ki*(R(t)+S(t)-Rs(t)-Fd(t))(3)
Figure BDA0003210652480000043
式中，Ri(t)为第t时段的壤中流，Ki为壤中流出流系数；Fd(t)为第t时段的下层自由水库入流；S(t-1)为第t-1时段的上层自由水库蓄水量；R(t-1)为第t-1时段的时段产流量；Ri(t-1)为第t-1时段的壤中流；Rs(t-1)为第t-1时段的地表径流；Fd(t-1)为第t-1时段的下层自由水库入流；S(0)为初始时刻上层自由水库蓄水量，根据流域的初始状态观测值或估计值来设定；
Figure BDA0003210652480000051
式中，K为土壤饱和水力传导度，Ψ为湿润锋处土壤吸力，△θ为土壤饱和含水率与田间持水量之差，F(t)为第t时段初的累计渗漏量；
Figure BDA0003210652480000052
式中，Fd(i)为第i时段的下层自由水库入流，F0为降雨初始时刻渗漏量，可设置为某一极小值，如0.001。
Figure BDA0003210652480000053
式中，Rg(t)为第t时段的地下径流，Sl(t)为第t时段初的下层自由水库蓄水量，SLM为下层自由水蓄水容量，Kg为地下径流出流系数。
步骤3中，计算下层自由水库蓄水量、下层自由水蓄水容量的计算方法为：
Figure BDA0003210652480000054
Figure BDA0003210652480000055
式中，Sl(0)为初始时刻的下层自由水库蓄水量，Zr为初始时刻的河道水位高程，Zg为初始时刻的地下水位高程，Zi为河道下伏含水层底边界高程，μ为地下水位变动带给水度；以上初始时刻状态变量值根据流域的初始状态观测值或估计值来设。
选择典型半干旱半湿润区海河流域内的清水河流域作为实施对象，采用本发明提出的基于双自由水库的产流算法和原始的基于单自由水库的产流算法分别计算两种不同产流算法下对应的3种径流成分，并采用统一的汇流算法(坡面及河网汇流均采用马斯京根算法)进行汇流计算。通过参数优选及调整，获得两种算法的最优洪水模拟效果。参考《水文情报预报规范GB/T 22482-2008》，选择洪峰流量相对误差、峰现时间误差、Nash-Sutcliffe效率系数和径流深误差作为评价指标，可以发现基于双自由水库的产流算法的模拟效果要优于原始的基于单自由水库的产流算法，如图3所示。
实施例2
本发明还提供了一种装置，用于计算产流。装置包括处理器和存储器；存储器中存储有程序或指令，程序或指令由处理器加载并执行以实现实施例1的基于双自由水库的产流计算方法。
实施例3
本发明还提供一种计算机可读存储介质，该计算机可读存储介质可以为非易失性计算机可读存储介质，该计算机可读存储介质也可以为易失性计算机可读存储介质，计算机可读存储介质中存储有指令，当指令在计算机上运行时，使得计算机执行实施例1的基于双自由水库的产流计算方法。
所属领域的技术人员可以清楚地了解到，本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的全部或部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器(read-only memory，ROM)、随机存取存储器(random access memory，RAM)、磁碟或者光盘等各种可以存储程序代码的介质。
以上所述仅是本发明的优选实施方式，并非用以限定本发明。本领域技术人员在不脱离本发明构思的前提下做出若干改进和优化，都应当视为本发明保护的范围。

Claims (9)

1.一种基于双自由水库的产流计算方法，其特征在于，包括以下步骤：
建立双自由水库产流结构：在新安江模型的三层包气带结构的基础上，在深层包气带之下设置张力水蓄满层，形成包括上层包气带、下层包气带、深层包气带和张力水蓄满层的四层包气带结构；将包气带内自由水能够占据的空间划分为上层自由水库和下层自由水库，其中上层包气带、下层包气带和深层包气带占据上层自由水库，张力水蓄满层占据下层自由水库；
使用蓄满产流算法，计算时段产流量；
根据建立的双自由水库产流结构，将时段产流量划分为地表径流、壤中流和地下径流三种径流成分并计算。
2.根据权利要求1所述的产流计算方法，其特征在于，根据建立的双自由水库产流结构，计算地表径流、壤中流和地下径流三种径流的方法是：
Figure FDA0003500805640000011
式中，Rs(t)为第t时段的地表径流，SM为上层自由水蓄水容量，R(t)为第t时段的时段产流量，S(t)为第t时段初的上层自由水库蓄水量；
Ri(t)＝Ki*(R(t)+S(t)-Rs(t)-Fd(t))
式中，Ri(t)为第t时段的壤中流，Ki为壤中流出流系数；Fd(t)为第t时段的下层自由水库入流；
Figure FDA0003500805640000012
式中，Rg(t)为第t时段的地下径流，Sl(t)为第t时段初的下层自由水库蓄水量，SLM为下层自由水蓄水容量，Kg为地下径流出流系数。
3.根据权利要求2所述的产流计算方法，其特征在于，第t时段的时段产流量R(t)的计算方法是：
Figure FDA0003500805640000013
式中，W(t)为第t时段的初始时刻的张力水蓄水量，WM为张力水蓄水容量，Pε(t)为扣除蒸散发损耗、植被冠层截留损耗后的第t时段的时段净降雨量。
4.根据权利要求2所述的产流计算方法，其特征在于，第t时段初的上层自由水库蓄水量S(t)的计算方法是：
Figure FDA0003500805640000014
式中，S(t-1)为第t-1时段的上层自由水库蓄水量；R(t-1)为第t-1时段的时段产流量；Ri(t-1)为第t-1时段的壤中流；Rs(t-1)为第t-1时段的地表径流；Fd(t-1)为第t-1时段的下层自由水库入流；S(0)为初始时刻上层自由水库蓄水量，根据流域的初始状态观测值或估计值来设定。
5.根据权利要求2所述的产流计算方法，其特征在于，第t时段的下层自由水库入流Fd(t)的计算方法为：
Figure FDA0003500805640000021
式中，K为土壤饱和水力传导度，Ψ为湿润锋处土壤吸力，△θ为土壤饱和含水率与田间持水量之差，F(t)为第t时段初的累计渗漏量。
6.根据权利要求5所述的产流计算方法，其特征在于，第t时段初的累计渗漏量F(t)的计算方法为：
Figure FDA0003500805640000022
式中，Fd(i)为第i时段的下层自由水库入流，F0为初始时刻的渗漏量。
7.根据权利要求2所述的产流计算方法，其特征在于，第t时段初的下层自由水库蓄水量Sl(t)、下层自由水蓄水容量SLM的计算方法为：
Figure FDA0003500805640000023
Figure FDA0003500805640000024
式中，Sl(0)为初始时刻的下层自由水库蓄水量，Zr为初始时刻的河道水位高程，Zg为初始时刻的地下水位高程，Zi为河道下伏含水层底边界高程，μ为地下水位变动带给水度；以上初始时刻状态变量值根据流域的初始状态观测值或估计值来设定。
8.一种产流计算装置，其特征在于，包括处理器和存储器；所述存储器中存储有程序或指令，所述程序或指令由所述处理器加载并执行以实现如权利要求1至7任一所述产流计算方法的步骤。
9.一种计算机可读存储介质，所述可读存储介质上存储程序或指令，所述程序或指令被处理器执行时实现如权利要求1至7任一所述产流计算方法的步骤。
CN202110930898.2A 2021-08-13 2021-08-13 一种基于双自由水库的产流计算方法、装置及存储介质 Active CN113657047B (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN202110930898.2A CN113657047B (zh) 2021-08-13 2021-08-13 一种基于双自由水库的产流计算方法、装置及存储介质
US17/877,164 US20230059267A1 (en) 2021-08-13 2022-07-29 Runoff yield calculation method and device based on double free reservoirs, and storage medium

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110930898.2A CN113657047B (zh) 2021-08-13 2021-08-13 一种基于双自由水库的产流计算方法、装置及存储介质

Publications (2)

Publication Number Publication Date
CN113657047A CN113657047A (zh) 2021-11-16
CN113657047B true CN113657047B (zh) 2022-03-22

Family

ID=78479813

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110930898.2A Active CN113657047B (zh) 2021-08-13 2021-08-13 一种基于双自由水库的产流计算方法、装置及存储介质

Country Status (2)

Country Link
US (1) US20230059267A1 (zh)
CN (1) CN113657047B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117494402A (zh) * 2023-10-24 2024-02-02 中国水利水电科学研究院 一种基于土壤蓄水量外溅的产流模式

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106548012A (zh) * 2016-10-17 2017-03-29 河海大学 一种通过流量过程线估算流域平均自由水蓄水容量的方法

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050230317A1 (en) * 2004-04-20 2005-10-20 Environmental Equipment Service Corp. Watershed runoff drainage device & method
WO2005116906A2 (en) * 2004-05-24 2005-12-08 Ch2M Hill, Inc. Method and system for water flow analysis
CN106777724B (zh) * 2016-12-23 2020-03-31 吉林大学 一种针对半湿润半干旱地区及湿润地区的基流分割方法
CN108874936B (zh) * 2018-06-01 2020-06-16 河海大学 一种基于改进新安江模型的适用于山丘区的水文预报方法
CN110928965B (zh) * 2018-12-13 2022-11-04 河海大学 一种基于流域精细分类的多模型灵活架构的模拟方法
CN109960861B (zh) * 2019-03-14 2023-05-09 南京大学 基于蓄满产流模式的地表地下水源分割方法
CN112905949B (zh) * 2021-03-05 2023-10-24 河海大学 一种基于流域下垫面特征的分布式产流参数估算方法
CN113011685A (zh) * 2021-04-27 2021-06-22 黄河勘测规划设计研究院有限公司 一种无径流资料地区内陆湖泊水位变化模拟预测方法

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106548012A (zh) * 2016-10-17 2017-03-29 河海大学 一种通过流量过程线估算流域平均自由水蓄水容量的方法

Also Published As

Publication number Publication date
CN113657047A (zh) 2021-11-16
US20230059267A1 (en) 2023-02-23

Similar Documents

Publication Publication Date Title
CN106884405B (zh) 一种无资料地区溃堤型山洪灾害分析评价方法
CN106202790B (zh) 一种分布式河北模型构建方法及其应用
Sun et al. Development and application of long-term root zone salt balance model for predicting soil salinity in arid shallow water table area
CN106971013A (zh) 山洪灾害小流域设计暴雨洪水计算方法
CN113657047B (zh) 一种基于双自由水库的产流计算方法、装置及存储介质
Kijne Water productivity under saline conditions.
CN114139258A (zh) 城市内涝建模评估方法、系统及计算机可读存储介质
Tao et al. Theoretical analysis and experimental verification of the improved subsurface drainage discharge with ponded water
Zhu et al. Estimation of groundwater residence time and recession rate in watershed of the Loess Plateau
CN109948220B (zh) 闸坝多目标泄流估算方法及系统
Koch et al. Continuous hydrologic modeling with HMS in the Aggtelek Karst region
CN105976103B (zh) 一种基于动态蓄水容量的洪水预报方法
Sana et al. Feasibility study of using treated wastewater to mitigate seawater intrusion along northern coast of Oman
Smith Rainfall and irrigation controls on groundwater rise and salinity risk in the Ord River Irrigation Area, northern Australia
Kupper et al. Using the SIMGRO regional hydrological model to evaluate salinity control measures in an irrigation area
Banerjee et al. Hydrogeological component assessment for water resources management of semi-arid region: a case study of Gwalior, MP, India
Kovacs et al. Conjunctive water management to sustain agricultural economic returns and a shallow aquifer at the landscape level
Kumar Recent studies on impact of climate change on groundwater resources
Luna Bharati et al. Analysis of the Inter-basin water transfer scheme in India: a case study of the Godavari–Krishna Link
Tayyab et al. Optimal application of conceptual rainfall–runoff hydrological models in the Jinshajiang River basin, China
Hublart et al. Modelling current and future trends in water availability for agriculture on a semi-arid and mountainous Chilean catchment
Temesgen et al. The impact of climate change on the availability of irrigation water at the rift valley lakes basin in Southern Ethiopia: A review
Dassargues Sustainable groundwater resources: how to quantify them, why is it so ‘location-dependent’?
Balasubramanian et al. Hydrology and its branches
I Naif A et al. Minimizing evaporation losses from small reservoirs (Wadi Houran-case study)

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