首页 / 专利库 / 地球科学 / 水文学 / 一种可定量描述优先流现象的运动波下渗方法

一种可定量描述优先流现象的运动波下渗方法

阅读:211发布:2020-05-19

专利汇可以提供一种可定量描述优先流现象的运动波下渗方法专利检索,专利查询,专利分析的服务。并且本 发明 公开了一种可定量描述优先流现象的运动波下渗方法,属于 土壤 水 文学领域。本方法通过将土壤水分剖面划分为优先流区域与基质流区域,根据运动波下渗理论提出一种优先流区域与基质流区域通过自由水面双向耦合的下渗模型。本发明提供的下渗方法能够反映湿润锋变化以及优先流与基质流的土壤水分互馈,解决了当前下渗方法中难以定量描述优先流现象的难题,具有较强的工程意义。,下面是一种可定量描述优先流现象的运动波下渗方法专利的具体信息内容。

1.一种可定量描述优先流现象的运动波下渗方法,其特征在于,包括如下步骤:
(1)采用染色实验方法将研究区域土壤剖面纵向划分为优先流区域p与基质流区域m,两区域面积比例为a∶(1-a),其中a为优先流区域面积比例,在此基础上采用双环入渗实验方法确定区域p与m的饱和传导度Konp与Konm;
(2)确定区域p与m饱和传导度随深度衰减系数fp与fm;
(3)假定蒸发E发生在土壤表层,净雨量为R=P-E,其中P为水文站观测降雨量;
(4)将运动波下渗理论分别应用于区域p与m,根据水量平衡原理,求得p区域上湿润锋位置 下湿润锋位置 m区域上湿润锋位置 下湿润锋位置 以及初始地下
水面位置
2.根据权利要求1所述的一种可定量描述优先流现象的运动波下渗方法,其特征在于,所述步骤(1)具体包括以下步骤:
在研究区域放置双环入渗仪,在两个克杯中放入亮蓝色染色示踪剂,每隔一分钟读取马克杯染色剂变化,根据达西定律计算得到Konp;
待染色剂完全渗入后,平衡15分钟,开挖土壤剖面,拍照记录染色剂染色范围,求得平均染色比例,即得优先流面积比例a;
用环刀切取染色区域土壤,在实验室做浸润实验,获得饱和水力传导度Konm。
3.根据权利要求1所述的一种可定量描述优先流现象的运动波下渗方法,其特征在于,所述步骤(2)中所述区域p与m饱和传导度随深度衰减系数fp与fm之间的关系为fp=α·fm,其中,α为小于0.5的参数,fm是经验系数,取值范围为0.0001/m-0.01/m。
4.根据权利要求1所述的一种可定量描述优先流现象的运动波下渗方法,其特征在于,所述步骤(3)所述净雨量由以下方式确定:
使用单层蒸发模型计算实际蒸发量,E=k·Ep,其中k为蒸散发折算系数,Ep为蒸发皿观测值,因此净雨量为R=P-E,其中:E为实际蒸发量,P为水文站观测降雨量。
5.根据权利要求1所述的一种可定量描述优先流现象的运动波下渗方法,其特征在于,所述步骤(4)包含以下步骤:
当湿润锋下界 时,其中 为初始地下水深,通过钻探获得,计算区域p与m
的上、下湿润锋位置 均使用如下公式:
式中:T为当前时刻,Δt为计算时间步长, 为T时刻湿润锋上界位置, 为T-1时刻湿润锋上界位置, 湿润锋下界, 为T-1时刻湿润锋下界, 为 在T时刻的下移速度, 为 在T时刻的下移速度;
其中
式中:Nt湿润锋上界,Nf湿润锋下界,Vt为Nt的下移速度,Vf为Nf的下移速度,Keq()等价饱和水力传导度的函数,Kon为饱和水力传导度,f为饱和水力传导度随深度衰减系数,N*为湿润锋的临界值,α为表面坡度,ε为土壤空隙排列指数,取4,R为净雨量,Ri为初始土壤含水量,取0.05,θ()为土壤含水量的函数,θs为饱和土壤含水量,θs通过实验室测得,其他参数均与上述相同;
当 且 时,根据水量平衡原理,计算公式如下:
式中: 为T时刻的地下水面位置, 为p区域T时刻的等价雨强,Konp为p区域的饱和水力传导度,fp为p区域饱和水力传导度随深度衰减系数,f1与f2为过程量见式(1-8)与(1-9);
其中
其中:Mu()为累积入渗量,θr为凋萎系数, 为T-1时刻的地下水面位置, 为T时刻的累积入渗量, 为T时刻等价雨强,ε为土壤孔隙指数;
又有如下关系
式中: 为p区域T时刻湿润锋上界的位置, 为初始时刻的地下水水位, 为p区域T时刻的湿润锋下界的位置;
方程组(1-7)为非线性方程组,方程组中 与 为未知解,用迭代法求数值解,其中 与 的求解用式(1-1)与(1-2)求解;
当 且 时,根据水量平衡公式计算如下:
其中
其中:fm为m区域的饱和水力传导度随深度衰减系数, 为m区域T时刻的等价雨强,为m区域T时刻湿润锋上界的位置, 为m区域T时刻湿润锋下界的位置;
又有如下关系
式中:方程组(1-13)为非线性方程组,方程组中 为未知变量,用牛顿
迭代法求数值解,其余参数与前述相同。

说明书全文

一种可定量描述优先流现象的运动波下渗方法

技术领域

[0001] 本发明涉及一种可定量描述优先流现象的运动波下渗方法,属于土壤文学领域。

背景技术

[0002] 优先流是指通过空隙、虫洞植物根系等通道,绕过土壤基质在土壤中快速通过的现象。优先流在自然界中广泛存在,且对土壤水文过程以及生态环境产生重要影响。但因其运动规律复杂,不遵循达西定律描述下的土壤水运动规律,因此难以对优先流进行模拟。
[0003] 目前,对优先流现象的认识主要停留在概念描述上,如下渗速度快,受到的阻较小等,但缺乏对其下渗过程的模拟。一般下渗理论如运动波下渗理论将土壤水运动看成一个整体,难以区分土壤分中的“快速流”与“慢速流”,因此无法对优先流现象进行有效刻画,不能反映实际的土壤水运动规律。

发明内容

[0004] 为解决现有技术的不足,本发明提出了一种可定量描述优先流现象的运动波下渗方法,解决现有技术中土壤水模拟过于概化、不能有效描述优先流现象的技术问题。
[0005] 本发明为解决其技术问题采用如下的技术方案:
[0006] 一种可定量描述优先流现象的运动波下渗方法,包括如下步骤:
[0007] (1)采用染色实验方法将研究区域土壤剖面纵向划分为优先流区域p与基质流区域m,两区域面积比例为a∶(1-a),其中a为优先流区域面积比例,在此基础上采用双环入渗实验方法确定区域p与m的饱和水力传导度Konp与Konm;
[0008] (2)确定区域p与m饱和传导度随深度衰减系数fp与fm;
[0009] (3)假定蒸发E发生在土壤表层,净雨量为R=P-E,其中P为水文站观测降雨量;
[0010] (4)将运动波下渗理论分别应用于区域p与m,根据水量平衡原理,求得p区域上湿润锋位置 下湿润锋位置 m区域上湿润锋位置 下湿润锋位置 以及初始地下水面位置
[0011] 所述步骤(1)具体包括以下步骤:
[0012] 在研究区域放置双环入渗仪,在两个克杯中放入亮蓝色染色示踪剂,每隔一分钟读取马克杯染色剂变化,根据达西定律计算得到Konp;
[0013] 待染色剂完全渗入后,平衡15分钟,开挖土壤剖面,拍照记录染色剂染色范围,求得平均染色比例,即得优先流面积比例a;
[0014] 用环刀切取染色区域土壤,在实验室做浸润实验,获得饱和水力传导度Konm。
[0015] 所述步骤(2)中所述区域p与m饱和传导度随深度衰减系数fp与fm之间的关系为fp=α·fm,其中,α为小于0.5的参数,fm是经验系数,取值范围为0.0001/m–0.01/m。
[0016] 所述步骤(3)所述净雨量由以下方式确定:
[0017] 使用单层蒸发模型计算实际蒸发量,E=k·Ep,其中k为蒸散发折算系数,Ep为蒸发皿观测值,因此净雨量为R=P-E,其中:E为实际蒸发量,P为水文站观测降雨量。
[0018] 所述步骤(4)包含以下步骤:
[0019] 当湿润锋下界 时,其中 为初始地下水深,通过钻探获得,计算区域p与m的上、下湿润锋位置 均使用如下公式:
[0020]
[0021]
[0022] 式中:T为当前时刻,Δt为计算时间步长, 为T时刻湿润锋上界位置, 为T–1时刻湿润锋上界位置, 湿润锋下界, 为T–1时刻湿润锋下界, 为 在T时刻的下移速度, 为 在T时刻的下移速度;
[0023] 其中
[0024]
[0025]
[0026]
[0027]
[0028] 式中:Nt湿润锋上界,Nf湿润锋下界,Vt为Nt的下移速度,Vf为Nf的下移速度,Keq()等价饱和水力传导度的函数,Kon为饱和水力传导度,f为饱和水力传导度随深度衰减系数,N*为湿润锋的临界值,α为表面坡度,ε为土壤空隙排列指数,一般取4,R为净雨量,Ri为初始土壤含水量,取0.05,θ()为土壤含水量的函数,θs为饱和土壤含水量,θs通过实验室测得,其他参数均与上述相同;
[0029] 当 且 时,根据水量平衡原理,计算公式如下:
[0030]
[0031] 式中: 为T时刻的地下水面位置, 为p区域T时刻的等价雨强,Konp为p区域的饱和水力传导度,fp为p区域饱和水力传导度随深度衰减系数,f1与f2为过程量见式(1–8)与(1–9),其余参数均与上述相同
[0032] 其中
[0033]
[0034]
[0035]
[0036]
[0037] 其中:Mu()为累积入渗量,θr为凋萎系数, 为T–1时刻的地下水面位置,为T时刻的累积入渗量, 为T时刻等价雨强,ε为土壤孔隙指数;
[0038] 又有如下关系
[0039]
[0040] 式中: 为p区域T时刻湿润锋上界的位置, 为初始时刻的地下水水位,为p区域T时刻的湿润锋下界的位置,
[0041] 方程组(1–7)为非线性方程组,方程组中 与 为未知解,用迭代法求数值解,其中 与 的求解用式(1–1)与(1–2)求解;
[0042] 当 且 时,根据水量平衡公式计算如下:
[0043]
[0044] 其中
[0045]
[0046]
[0047]
[0048]
[0049] 其中:fm为m区域的饱和水力传导度随深度衰减系数, 为m区域T时刻的等价雨强, 为m区域T时刻湿润锋上界的位置, 为m区域T时刻湿润锋下界的位置;
[0050] 又有如下关系
[0051]
[0052]
[0053] 式中:方程组(1–13)为非线性方程组,方程组中 为未知变量,用牛顿迭代法求数值解,其余参数与前述相同。
[0054] 本发明的有益效果如下:
[0055] 本发明公开了一种可定量描述优先流现象的运动波下渗方法,通过将土壤水分剖面划分为优先流区域与基质流区域,根据运动波下渗理论提出一种优先流区域与基质流区域通过自由水面双向耦合的下渗模型。本发明提供的下渗方法能够反映湿润锋变化以及优先流与基质流的土壤水分互馈,解决了当前下渗方法中无法定量描述优先流现象的难题,具有较强的工程意义。附图说明
[0056] 图1为 时优先流区域与基质流区域下渗过程示意图。
[0057] 图2为 且 时优先流区域与基质流区域下渗过程示意图。
[0058] 图3为 且 时优先流区域与基质流区域下渗过程示意图。
[0059] 图4为根据南方某山坡区域一场降雨所模拟的优先流与基质流湿润锋随深度方向的变化图。

具体实施方式

[0060] 下面结合附图和具体实施方式对本发明作进一步描述。以下实施例仅用于更加清楚地说明本发明的技术方案,而不能以此来限制本发明的保护范围。
[0061] 下面结合实例对本发明作更进一步的说明。
[0062] 现有南方某一植被覆盖较好的山坡,根据本发明方法,该下渗方法的计算方法如下:
[0063] (1)在研究区域放置双环入渗仪,在两个马克杯中放入亮蓝色染色示踪剂,保证内外环水位稳定,每隔一分钟读取马克杯染色剂变化,根据达西定律计算得到Konp=230mm/h;
[0064] (2)待染色剂完全渗入后,平衡15分钟,开挖土壤剖面,拍照记录染色剂染色范围,求得平均染色比例,即得优先流面积比例a=32%;
[0065] (3)用环刀切取染色区域土壤,在实验室做浸润实验,获得饱和水力传导度Konm=70mm/h。
[0066] (4)根据经验取fm=0.01/m,α=0.1,则fp=0.001/m。
[0067] (5)取其中一场降雨过程的降雨资料P与蒸发皿观测资料Ep,
[0068] 根据R=P-k·Ep,蒸散发折算系数k可取1,可得净雨量过程。
[0069] (6)根据电钻获得土壤厚度,可将土壤底层深度近似认为地下水深度,即得[0070] (7)根据以下公式即可计算得到优先流湿润锋 以及基质流湿润锋随时间变化,最终模拟结果如图4所示,
[0071] 当下湿润锋 时,模拟结果如图1所示,其中 为初始地下水深,可通过钻探获得,计算区域p与m的上、下湿润锋位置 均可使用如下公式:
[0072]
[0073]
[0074] 式中:T为当前时刻,Δt为计算时间步长, 为T时刻湿润锋上界位置, 为T–1时刻湿润锋上界位置, 湿润锋下界, 为T–1时刻湿润锋下界, 为 在T时刻的下移速度, 为 在T时刻的下移速度。
[0075] 其中
[0076]
[0077]
[0078]
[0079]
[0080] 式中:Nt湿润锋上界,Nf湿润锋下界,Vt为Nt的下移速度,Vf为Nf的下移速度,Keq()等价饱和水力传导度的函数·,Kon为饱和水力传导度,f为饱和水力传导度随深度衰减系数,N*为湿润锋的临界值,α为表面坡度,ε为土壤空隙排列指数,一般可取4,R为净雨量,Ri为初始土壤含水量,取0.05,θ()为土壤含水量的函数,θs为饱和土壤含水量,θs可通过实验室测得,其他参数均与上述相同。
[0081] 当 且 时,模拟结果如图2所示,根据水量平衡原理,计算公式如下:
[0082]
[0083] 式中: 为T时刻的地下水面位置, 为p区域T时刻的等价雨强,Konp为p区域的饱和水力传导度,fp为p区域饱和水力传导度随深度衰减系数,f1与f2为过程量见式(1–8)与(1–9),其余参数均与上述相同。
[0084] 其中
[0085]
[0086]
[0087]
[0088]
[0089] 其中:Mu()为累积入渗量,θr为凋萎系数, 为T–1时刻的地下水面位置,为T时刻的累积入渗量, 为T时刻等价雨强,ε为土壤孔隙指数。
[0090] 又有如下关系
[0091]
[0092] 式中: 为p区域T时刻湿润锋上界的位置, 为初始时刻的地下水水位,为p区域T时刻的湿润锋下界的位置,其余参数与前述相同。
[0093] 方程组(1–7)为非线性方程组,方程组中 与 为未知解,可用牛顿迭代法求数值解,其中 与 的求解可用式(1–1)与(1–2)求解。
[0094] 当 且 时,模拟结果如图3所示,根据水量平衡公式计算如下:
[0095]
[0096] 其中
[0097]
[0098]
[0099]
[0100]
[0101] 其中:fm为m区域的饱和水力传导度随深度衰减系数, 为m区域T时刻的等价雨强, 为m区域T时刻湿润锋上界的位置, 为m区域T时刻湿润锋下界的位置。
[0102] 又有如下关系
[0103]
[0104]
[0105] 式中:方程组(1–13)为非线性方程组,方程组中 为未知变量,可用牛顿迭代法求数值解,其余参数与前述相同。
[0106] 以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明技术原理的前提下,还可以做出若干改进和变形,这些改进和变形也应视为本发明的保护范围。
高效检索全球专利

专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台,是提供专利分析,专利查询,专利检索等数据服务功能的知识产权数据服务商。

我们的产品包含105个国家的1.26亿组数据,免费查、免费专利分析。

申请试用

分析报告

专利汇分析报告产品可以对行业情报数据进行梳理分析,涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度,系统通过AI智能系统对图表进行解读,只需1分钟,一键生成行业专利分析报告。

申请试用

QQ群二维码
意见反馈