CN110990753A - 一种暗管系统淋洗效果的计算方法 - Google Patents

一种暗管系统淋洗效果的计算方法 Download PDF

Info

Publication number
CN110990753A
CN110990753A CN201911125274.2A CN201911125274A CN110990753A CN 110990753 A CN110990753 A CN 110990753A CN 201911125274 A CN201911125274 A CN 201911125274A CN 110990753 A CN110990753 A CN 110990753A
Authority
CN
China
Prior art keywords
salt
calculating
mass transfer
formula
concealed
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.)
Granted
Application number
CN201911125274.2A
Other languages
English (en)
Other versions
CN110990753B (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 CN201911125274.2A priority Critical patent/CN110990753B/zh
Publication of CN110990753A publication Critical patent/CN110990753A/zh
Application granted granted Critical
Publication of CN110990753B publication Critical patent/CN110990753B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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
    • 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
    • Y02A20/00—Water conservation; Efficient water supply; Efficient water use
    • Y02A20/152—Water filtration

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Data Mining & Analysis (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Computational Mathematics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Algebra (AREA)
  • Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)

Abstract

本发明公开了一种暗管系统淋洗效果的计算方法，包括以下步骤：(1)建立二维暗管排盐系统地下水渗流模型；(2)推导地下水渗流复势函数；(3)利用复势函数计算地下水渗流量；(4)引入双重孔隙介质模型，计算流动区内盐分淋洗时间；(5)建立双重孔隙介质传质模型；(6)根据双重孔隙介质传质模型，推导非流动区脱盐过程；(7)推导盐分在暗管系统内受动力学传质过程影响的淋洗时间，总需水量。本发明可以计算受动力学传质过程影响下暗管系统淋洗的时间和总水量，能够评估脱盐地块是否适合使用淋洗方式进行改良，为盐碱地开发及规划提供指导。

Description

一种暗管系统淋洗效果的计算方法
技术领域
本发明涉及一种暗管系统的计算方法，尤其涉及一种暗管系统淋洗效果的计算方法。
背景技术
土壤盐碱化是一种在全世界范围内普遍产生的土壤退化现象。一方面由于岩石、土壤受到风化、降雨侵蚀作用，产生大量无机盐离子从而富集在表层土壤，造成原生盐碱化。另一方面，由于地下咸水水位上升和灌溉水将盐分带入表层土壤，经过蒸发过程，大量盐分来不及向深层土壤运动便富集在表层土壤，造成次生盐碱化。土壤盐碱化严重威胁土壤生物圈生态平衡，减少作物产量和减少土壤微生物多样性，造成土壤进一步退化。据统计，全球盐碱地总面积约为955万 km2，约占据全球陆地总面积的10％。值得指出的是，我国耕地盐碱化问题尤为严重，据不完全统计我国盐碱地总面积约有3.6亿亩，占据我国可耕地面积的1/5。
目前，暗管排水系统被广泛应用于淋洗土壤中盐分，改善土壤质地。然而，水力淋洗的效率往往受到非流动域的动力学传质过程的影响。众所周知，土壤渗透性具有空间非均质性，包括高渗透性的大孔隙和低渗透性的聚合小孔隙土体。当淋洗效果发生在有结构的土壤时，聚合小孔隙土体通常作为非流动区，而大孔隙可以作为优先流通道，快速传递水量和溶质。研究表明，70％-85％的水流通过大孔隙从表层土壤传导到深层土壤。这意味着聚合小孔隙土壤内水流流速较慢且流量较小，从而造成对此区域盐分淋洗的严重滞后，影响脱盐地块的淋洗效率。
尽管不能完全刻画表层土壤的非均质性，然而双重孔隙介质模型能够合理地描述大孔隙优先流作用和流动域与非流动域之间的传质作用，能够显示盐分淋洗效果中受到动力学传质过程影响。描述受到传质限制过程中对盐分运动的影响能够合理地反应淋洗策略的效果，能够为盐碱地开发规划提供指导。然而，目前尚未有耦合暗管系统地下水流动和双重孔隙介质传质模型的解析解，而水文地质参数变化(如土壤饱和渗透系数，非流动区孔隙率，传质速率等)对暗管排盐系统脱盐效果的影响更不得而知。
发明内容
发明目的：上述现有技术的不足，本发明旨在提供一种暗管系统盐分淋洗效果受动力学传质过程影响的解析计算方法。
技术方案：本发明的一种暗管系统淋洗效果的计算方法，包括如下步骤：
(1)建立二维暗管排盐系统地下水渗流模型；假设饱和稳态的地下水流动，各向同性均质的渗透介质，土壤表面入渗过程忽略横向流动；渗流模型的下边界和左右边界设置为无流动边界，暗管外壁设置为等势面，上边界设置为定水头边界，水头为0。
(2)利用保角变换将渗流平面进行保角映射，同时根据渗流平面绘制势能平面，建立保角映射后平面与势能平面之间的关系，推导出地下水渗流复势函数为：
Figure BDA0002276631510000021
式中，
Figure BDA0002276631510000023
为势函数；ψ为流函数；q单位时间单位长度的暗管的渗流量，L2/T；sn(u)和dn(u)为雅克比椭圆函数；z为复函数，即z＝x+iy，x为直角坐标系x轴方向，y为直角坐标系y轴方向，i2＝-1；S为暗管间距的一半长度，L；d为暗管埋设深度，L；r为圆柱形排水暗管直径，L；H为不透水层距离土壤表面的距离， L；K为模数为m的第一类完全椭圆积分；K'为模数为m'＝1-m的第一类完全椭圆积分。
(3)利用复势函数计算地下水渗流量；根据地下水渗流复势函数w推导入渗量q的计算公式为
Figure BDA0002276631510000022
式中，k为土壤饱和渗透系数，L/T。
特别地，在步骤2和3中，m可通过如下公式计算
K/K'＝2H/S
将求解的m带入方程q中，即可求得模型渗流总量；将求解的m带入方程 w中，求出入渗量q；进一步求得x＝L和x＝L+△S之间的入渗量
f＝(qL-qL+ΔS)/q
其中，qL为距离暗管距离L处土壤表面入渗量，qL+△S为距离暗管距离L+ △S处土壤表面入渗量，ΔS为土壤表面任意两条流线之间的宽度。
(4)考虑脱盐系统内总盐度(包括流动区与非流动区)，计算流动区内盐分淋洗时间；在双重孔隙介质中，流动区内盐分淋洗时间tl计算公式为
Figure BDA0002276631510000031
式中，θm为流动区孔隙度，θim为非流动区孔隙度，h为盐分目标淋洗深度，L。
(5)建立双重孔隙介质的传质模型，流动区与非流动区浓度差为非流动区内的盐分扩散到流动区提供动力，控制方程为
Figure BDA0002276631510000032
中，cm为流动区内盐分浓度，M/L3；cim为非流动区内盐分浓度，M/L3；α为传质速率。
(6)根据双重孔隙介质的传质模型，推导盐分从非流动区通过传质方式扩散到流动区的时间。特别地，假设流动区被水流冲刷过后，其盐分会迅速随着水流向深层土层运动，则可将该区域盐分浓度简化为0，即cm＝0
将cm＝0代入方程θim中，求解可推导盐分从非流动区通过传质方式扩散到流动区的时间计算公式为
Figure BDA0002276631510000033
特别地，式中无量纲参数β和τ的计算公式为
Figure BDA0002276631510000034
Figure BDA0002276631510000035
Figure BDA0002276631510000036
式中，c0 crit为非流动区内初始盐分浓度，M/L3；cim 0为非流动区内初始溶质浓度；ccrit为适合作物生长的盐分浓度，M/L3；tre为某一常数时间，β为双重孔隙介质的容量系数，τim为盐分传质的时间尺度，τ无量纲时间参数。
(7)推导盐分在暗管系统内受动力学传质过程影响的淋洗时间为，
Figure BDA0002276631510000041
总需水量为，
Figure BDA0002276631510000042
有益效果：与现有技术相比，本发明具有如下显著优点：
(1)本发明首次推导了暗管淋洗效果受传质限制下淋洗时间和总需水量解析；可以为盐碱地开发与改良提供指导，可以准确判断盐碱地是否适合使用淋洗法改良。
(2)本发明的解析解能刻画水文地质参数(如含水层介质饱和渗透系数、传质速率、非流动区孔隙度等)变化导致的盐分淋洗效率的变化，增加了该解析解的适应性；
附图说明
图1为暗管系统地下水渗流示意图；
图2为双重孔隙介质传质模型示意图；
图3为解析解推导与计算流程图；
图4为解析解与数值解的拟合结果。
具体实施方式
下面结合附图对本发明的技术方案作进一步说明。
如图1所示，为淋洗效果中地下水渗流规律。将淡水或再生水通过漫灌方式输送到土壤表面11，使土壤表面一次性完全被淹没。淋洗水源穿过土壤表面11 进入土壤层12，溶解并冲刷盐分离子，混合盐分随着流线13的轨迹，排入到暗管10中，并最终排出脱盐系统。
根据根据图1概念模型完成如下步骤：
包括如下步骤：
(1)建立二维暗管排盐系统地下水渗流模型；假设饱和稳态的地下水流动，各向同性均质的渗透介质，土壤表面入渗过程忽略横向流动；渗流模型的下边界和左右边界设置为无流动边界，暗管外壁设置为等势面，上边界设置为定水头边界，水头为0。
(2)利用保角变换将渗流平面进行保角映射，同时根据渗流平面绘制势能平面，建立保角映射后平面与势能平面之间的关系，推导出地下水渗流复势函数为：
Figure BDA0002276631510000051
式中，
Figure BDA0002276631510000053
为势函数；ψ为流函数；q单位时间单位长度的暗管的渗流量，L2/T；sn(u)和dn(u)为雅克比椭圆函数；z为复函数，即z＝x+iy，x为直角坐标系x轴方向，y为直角坐标系y轴方向，i2＝-1；S为暗管间距的一半长度，L；d为暗管埋设深度，L；r为圆柱形排水暗管直径，L；H为不透水层距离土壤表面的距离， L；K为模数为m的第一类完全椭圆积分；K'为模数为m'＝1-m的第一类完全椭圆积分。
(3)利用复势函数计算地下水渗流量；根据地下水渗流复势函数w推导入渗量q的计算公式为
Figure BDA0002276631510000052
式中，k为土壤饱和渗透系数，L/T。
特别地，在步骤2和3中，m可通过如下公式计算
K/K'＝2H/S (3)
将求解的m带入方程2中，即可求得模型渗流总量；将求解的m带入方程 (1)中，求出入渗量q；进一步求得x＝L和x＝L+△S之间的入渗量
f＝(qL-qL+ΔS)/q (4)
其中，qL为距离暗管距离L处土壤表面入渗量，qL+△S为距离暗管距离L+ △S处土壤表面入渗量，ΔS为土壤表面任意两条流线之间的宽度。
如图2所示，假设土壤中每一单元体内孔隙均有两部分构成，一部分为流动区20，水流流动仅仅发生在流动区20，在该区域内溶质的运动受到水动力弥散作用驱使；另一部分为非流动区22，水流在非流动区22内不发生流动，溶质的迁移只通过流动区20与非流动区22之间的动力学传质过程21。考虑系统双重孔隙性，完成如下步骤：
(4)考虑脱盐系统内总盐度(包括流动区与非流动区)，计算流动区内盐分淋洗时间；在双重孔隙介质中，流动区内盐分淋洗时间tl计算公式为
Figure BDA0002276631510000061
式中，θm为流动区孔隙度，θim为非流动区孔隙度，h为盐分目标淋洗深度， L；△S为土壤表面任意两条流线之间的宽度，L。
(5)建立双重孔隙介质的传质模型，流动区与非流动区浓度差为非流动区内的盐分扩散到流动区提供动力，控制方程为
Figure BDA0002276631510000062
中，cm为流动区内盐分浓度，M/L3；cim为非流动区内盐分浓度，M/L3；α为传质速率。
(6)根据双重孔隙介质的传质模型，推导盐分从非流动区通过传质方式扩散到流动区的时间。特别地，假设流动区被水流冲刷过后，其盐分会迅速随着水流向深层土层运动，则可将该区域盐分浓度简化为0，即cm＝0
将方程6代入方程5中，求解方程5可推导盐分从非流动区通过传质方式扩散到流动区的时间计算公式为
Figure BDA0002276631510000063
特别地，式中无量纲参数β和τ的计算公式为
Figure BDA0002276631510000064
Figure BDA0002276631510000065
Figure BDA0002276631510000066
式中，c0 crit为非流动区内初始盐分浓度，M/L3；cim 0为非流动区内初始溶质浓度；ccri为适合作物生长的盐分浓度，M/L3；tre为某一常数时间，β为双重孔隙介质的容量系数，τim为盐分传质的时间尺度，τ无量纲时间参数。
(7)推导盐分在暗管系统内受动力学传质过程影响的淋洗时间为，
Figure BDA0002276631510000071
总需水量为，
Figure BDA0002276631510000072
为了验证解析解的准确性，将解析计算结果与数值模拟计算结果进行比对。根据崇明岛东部区域的现场勘测数据，脱盐地块渗透系数为1m/d，流动区孔隙度为0.3，非流动区孔隙度为0.15，纵向弥散系数DL＝0.1m，横向弥散系数DT为0.01m,含水层地下水初始浓度(流动域与非流动域)为10kg/m3。暗管埋深为1m，暗管直径为8cm，暗管间距为20m，不透水层位置为-5m，土壤表面定水头高5cm。选取脱盐地块二维横切剖面的一半进行数值模拟，模型的尺寸为 10m×5m。
首先计算解析结果。根据关系式K/K'＝2H/S计算m的值，将m带入到公式 1与公式2中，计算渗流量与土壤表面流量分配关系。特别地，选取两条暗管中间区域间隔为0.04m的两条流线之间的流量为f，计算流动域内盐分淋洗时间。值得注意的是，由于自然界中传质速率很难确定，根据方程11和方程12，本案例将计算传质速率为4d--0.00004d-范围内所有淋洗时间和淋洗总需水量。
然后，我们运行了12组数值模拟去拟合解析结果。其中，12组数值模型的传质速率分布为4d-,0.4d-,0.1269d-,0.04d-,0.0225d-,0.0127d-,0.0071d-,0.004 d-,0.00225d-,0.00127d-,0.00071d-,0.0004d-,0.000225d-,0.000127d-,0.000071 d-,0.00004d-，其余参数与解析模型相同。如图3所示，解析计算结果与数值模拟结果拟合程度相当好，显示了解析解的准确性。此外，结果还展示了双重孔隙介质中动力学传质过程会增加淋洗周期和总需水量，尤其是传质速率比较小时，会大大减少淋洗脱盐的效率，造成水资源和资金的大量浪费。
本发明首次推导了暗管淋洗效果受传质限制下淋洗时间和总需水量解析；本发明的解析解能刻画水文地质参数(如含水层介质饱和渗透系数、传质速率、非流动区孔隙度等)变化导致的盐分淋洗效率的变化，增加了该解析解的适应性；本发明可以为盐碱地开发与改良提供指导，可以准确判断盐碱地是否适合使用淋洗法改良。

Claims (10)

1.一种暗管系统淋洗效果的计算方法，其特征在于，包括以下步骤：
(1)建立二维暗管排盐系统地下水渗流模型；
(2)推导地下水渗流复势函数；
(3)利用复势函数计算地下水渗流量；
(4)引入双重孔隙介质模型，计算流动区内盐分淋洗时间；
(5)建立双重孔隙介质传质模型；
(6)根据双重孔隙介质传质模型，推导非流动区脱盐过程；
(7)推导盐分在暗管系统内受动力学传质过程影响的淋洗时间，总需水量。
2.根据权利要求1所述的暗管系统淋洗效果的计算方法，其特征在于，在步骤(1)中，地下水为饱和稳态流动，渗透介质为各向同性均质，忽略土壤表面入渗过程横向流动；渗流模型的下边界和左右边界设置为无流动边界，暗管外壁设置为等势面，上边界设置为定水头边界，水头为0。
3.根据权利要求1所述的暗管系统淋洗效果的计算方法，其特征在于：在步骤(2)中，利用保角变换将渗流平面进行保角映射，同时根据渗流平面绘制势能平面，建立保角映射后平面与势能平面之间的关系，推导出地下水渗流复势函数为：
Figure FDA0002276631500000011
式中，
Figure FDA0002276631500000012
为势函数；ψ为流函数；q为单位时间单位长度的暗管的渗流量；sn(u)和dn(u)为雅克比椭圆函数；z为复函数，即z＝x+iy；x为直角坐标系x轴方向，y为直角坐标系y轴方向，i2＝-1；S为暗管间距的一半长度；d为暗管埋设深度；r为圆柱形排水暗管直径；H为不透水层距离土壤表面的距离；K为模数为m的第一类完全椭圆积分；K'为模数为m'的第一类完全椭圆积分，其中m和m'为椭圆积分参数。
4.根据权利要求1所述的暗管系统淋洗效果的计算方法，其特征在于：在步骤(3)中，根据地下水渗流复势函数w推导入渗量q的计算公式为
Figure FDA0002276631500000021
式中，k为土壤饱和渗透系数。
5.根据权利要求1所述的暗管系统淋洗效果的计算方法，其特征在于：在步骤(2)和(3)中，m通过如下公式计算
K/K'＝2H/S
其中，m'通过如下公式计算
m'＝1-m
将求解的m和m'带入q方程与w方程中，求出入渗量q；进一步求得x＝L和x＝L+△S之间的入渗量
f＝(qL-qL+ΔS)/q
其中，qL为距离暗管距离L处土壤表面入渗量，qL+△S为距离暗管距离L+△S处土壤表面入渗量，△S为土壤表面任意两条流线之间的宽度。
6.根据权利要求1所述的暗管系统淋洗效果的计算方法，其特征在于：在步骤(4)中，双重孔隙介质中，流动区内盐分淋洗时间tl计算公式为
Figure FDA0002276631500000022
式中，θm为流动区孔隙度，θim为非流动区孔隙度，h为盐分目标淋洗深度。
7.根据权利要求1所述的暗管系统淋洗效果的计算方法，其特征在于：在步骤(5)中，建立双重孔隙介质的传质模型，流动区与非流动区浓度差为非流动区内的盐分扩散到流动区提供动力，控制方程为
Figure FDA0002276631500000023
式中，cm为流动区内盐分浓度，cim为非流动区内盐分浓度，α为传质速率。
8.根据权利要求1所述的暗管系统淋洗效果的计算方法，其特征在于：流动区内盐分浓度为0，即
cm＝0
将cm＝0代入θim方程中，求得盐分从非流动区通过传质方式扩散到流动区的时间tm计算公式为
Figure FDA0002276631500000031
式中无量纲参数β和τ的计算公式为
Figure FDA0002276631500000032
Figure FDA0002276631500000033
Figure FDA0002276631500000034
式中，c0 crit为非流动区内初始盐分浓度，cim 0为非流动区内初始溶质浓度，ccrit为适合作物生长的盐分浓度，tre为某一常数时间，β为双重孔隙介质的容量系数，τim为盐分传质的时间尺度，τ为无量纲时间参数。
9.根据权利要求1所述的暗管系统淋洗效果的计算方法，其特征在于：在步骤(6)中，盐分在暗管系统内受动力学传质过程影响的淋洗总时间计算公式为
Figure FDA0002276631500000035
10.根据权利要求1所述的暗管系统淋洗效果的计算方法，其特征在于：总需水量Q计算公式为
Figure FDA0002276631500000036
CN201911125274.2A 2019-11-18 2019-11-18 一种暗管系统淋洗效果的计算方法 Active CN110990753B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911125274.2A CN110990753B (zh) 2019-11-18 2019-11-18 一种暗管系统淋洗效果的计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911125274.2A CN110990753B (zh) 2019-11-18 2019-11-18 一种暗管系统淋洗效果的计算方法

Publications (2)

Publication Number Publication Date
CN110990753A true CN110990753A (zh) 2020-04-10
CN110990753B CN110990753B (zh) 2022-03-29

Family

ID=70084698

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911125274.2A Active CN110990753B (zh) 2019-11-18 2019-11-18 一种暗管系统淋洗效果的计算方法

Country Status (1)

Country Link
CN (1) CN110990753B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111898257A (zh) * 2020-07-16 2020-11-06 武汉大学 区域暗管布局及排水排盐数值模拟方法和装置
CN112462032A (zh) * 2020-11-12 2021-03-09 武汉大学 适用于盐渍土地区评价暗管排水排盐效果的方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107155431A (zh) * 2017-05-18 2017-09-15 深圳市微润灌溉技术有限公司 微润驱盐管及其局部淡化盐碱地方法
CN108112294A (zh) * 2017-12-23 2018-06-05 广东粤壮园林建设有限公司 一种盐碱地改良方法
CN108990452A (zh) * 2018-05-25 2018-12-14 河海大学 一种耦合潜水层咸淡水替换及淋洗脱盐装置
CN110230342A (zh) * 2019-04-26 2019-09-13 天津生态城市政景观有限公司 盐碱区域地上地下联合水资源调蓄利用系统

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107155431A (zh) * 2017-05-18 2017-09-15 深圳市微润灌溉技术有限公司 微润驱盐管及其局部淡化盐碱地方法
CN108112294A (zh) * 2017-12-23 2018-06-05 广东粤壮园林建设有限公司 一种盐碱地改良方法
CN108990452A (zh) * 2018-05-25 2018-12-14 河海大学 一种耦合潜水层咸淡水替换及淋洗脱盐装置
CN110230342A (zh) * 2019-04-26 2019-09-13 天津生态城市政景观有限公司 盐碱区域地上地下联合水资源调蓄利用系统

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
MGUIDICH BELHAJ AMEL ET.AL: "Assessing HYDRUS-2D model to simulate water content and salt accumulation in a soil irrigated with a subsurface drip system: Application in a semiarid area of central tunisia", 《IRRIGATION AND DRAINAGE》 *
瞿兴业: "渠道和水田渗漏下的暗管排水计算", 《水利学报》 *
石培君等: "基于HYDRUS模型的暗管排水水盐运移模拟", 《干旱地区农业研究》 *
韩翀等: "暴雨条件下小流域出口氮污染物流失特征及分析", 《科技视界》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111898257A (zh) * 2020-07-16 2020-11-06 武汉大学 区域暗管布局及排水排盐数值模拟方法和装置
CN111898257B (zh) * 2020-07-16 2022-03-04 武汉大学 区域暗管布局及排水排盐数值模拟方法和装置
CN112462032A (zh) * 2020-11-12 2021-03-09 武汉大学 适用于盐渍土地区评价暗管排水排盐效果的方法
CN112462032B (zh) * 2020-11-12 2022-03-04 武汉大学 适用于盐渍土地区评价暗管排水排盐效果的方法

Also Published As

Publication number Publication date
CN110990753B (zh) 2022-03-29

Similar Documents

Publication Publication Date Title
Shamrukh et al. Modeling the effect of chemical fertilizers on ground water quality in the Nile Valley Aquifer, Egypt
CN106600035A (zh) 基于污染物迁移模拟的水源地水质安全预警方法
CN105426568B (zh) 一种估算喀斯特地区土壤流失量的方法
CN110990753B (zh) 一种暗管系统淋洗效果的计算方法
JP2013535597A (ja) 液体注入によるその位置での溶脱方法
CN109034656A (zh) 一种地下水脆弱性评价方法
CN113239598B (zh) 一种应用数值模拟的地下水污染源空间综合识别方法
Bouwer Ground water recharge design for renovating waste water
Saatsaz et al. The application of groundwater modelling to simulate the behaviour of groundwater resources in the Ramhormooz Aquifer, Iran
Hu et al. Development of a groundwater flow and reactive solute transport model in the Yongding River alluvial fan, China
Merdun et al. Effect of the initial soil moisture content on the spatial distribution of the water retention
CN110771291B (zh) 一种提高淋洗排盐效率的覆盖层分区方法
Yu et al. Method for improving leaching efficiency of coastal subsurface drainage systems
Blanco et al. Water erosion
CN106242029A (zh) 一种坎儿井式地下水污染修复系统
Weiler et al. Connectivity due to preferential flow controls water flow and solute transport at the hillslope scale
CN116467958A (zh) 盐湖卤水数值模型构建与补水溶矿效率计算方法
CN109147047A (zh) 一种渗渠型傍河水源地数值模型构建方法
Turgunbaeva et al. The problem of construction on salted soils due to insufficient use of Underground Water of the Bukhara Region
CN201802363U (zh) 原地浸取注液系统
CN209368013U (zh) 地下水硝酸盐串联式原位抽出耦合处理系统
CN110616772B (zh) 一种用于雨水存储和监测的回灌系统及该系统的建设方法
Abo Shaeshaa et al. Control of Contaminant transport through Different Horizontal Layered Soil Using Sheet Pile.
Osmanaj et al. Determination of Groundwater Protection Zones of the Pozharan Wellfield Using Hydrogeological Modflow Model
Bekbayev et al. Factors that intensify soil degradation in the Kazakhstan part of the Golodnostepsky irrigation massif

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
CB03 Change of inventor or designer information

Inventor after: Zhang Jiaxu

Inventor after: Lu Chunhui

Inventor after: Chen Qian

Inventor after: Chen Qinyu

Inventor before: Lu Chunhui

Inventor before: Zhang Jiaxu

Inventor before: Chen Qian

Inventor before: Chen Qinyu

CB03 Change of inventor or designer information
GR01 Patent grant
GR01 Patent grant