CN104865165A - 全过程通量守恒的现场土壤溶质运移弥散系数测定的方法 - Google Patents
全过程通量守恒的现场土壤溶质运移弥散系数测定的方法 Download PDFInfo
- Publication number
- CN104865165A CN104865165A CN201510307246.8A CN201510307246A CN104865165A CN 104865165 A CN104865165 A CN 104865165A CN 201510307246 A CN201510307246 A CN 201510307246A CN 104865165 A CN104865165 A CN 104865165A
- Authority
- CN
- China
- Prior art keywords
- solute
- dispersion coefficient
- soil
- partiald
- formula
- 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
Links
Abstract
本发明公开了一种全过程通量守恒的现场土壤溶质运移弥散系数测定的方法。首先进行地下水流场识别,确定地下水中溶质在土壤中运动的主要和次要迁移方向;然后基于溶质的通量守恒,建立模拟土壤溶质运动的数学模型;确定数学模型的现场实验条件;在控制边界条件下,以注入源的形式注入非转化型溶质,采用径向流动观测方式,测定非转化型溶质在主要和次要迁移方向上浓度随时间变化的情况;绘制弥散系数与孔隙平均流速的关系图并分析土壤和孔隙平均流速对弥散系数的影响,基于非转化型溶质全过程的通量守恒确定弥散系数。本发明建立数学模型,模拟土壤溶质运动,易于操作、测量精度高,实验模型容易搭建,条件易于实现,样本丰富多样,结果直观。
Description
技术领域
本发明涉及水环境保护领域,特别是涉及一种全过程通量守恒的现场土壤溶质运移弥散系数测定的方法。
背景技术
目前,地下水质的污染,例如排放污水或高浓度回灌水对地下水水质影响的预测与评价的研究,对于控制或预防地下水污染的发生。水质数学模型对于地下水的水质预测与评价至关重要,而地下水弥散系数即是建立这个模型的一个重要而不可缺少的参数。地下水弥散系数是表征溶质态污染物在地下水中运移传播的一种参数。
弥散系数的测定离不开以下几个方面技术的进步与支撑:(1)地下水含水层信息,包括含水层厚度、性质、顶底板埋深,隔水层与含水层结构,含水层形状、延伸方向和稳定性以及有关参数难以确定;(2)地下水补给排泄条件,地下水流速流向,主要含水层的富水性与透水性及溶质态污染物的水动力弥散过程等也对测定水动力弥散系数起到了约束作用;(3)应用数学模型,根据土壤溶质运动进行参数反演是实现准确测定水动力弥散系数的重要前提。由于以上三个方面的技术需求,应用数学模型,模拟土壤溶质运动实现土壤水动力弥散系数的精确计算,在现场条件下,对于影响弥散系数的地下水的流域进行直接测定基本上不可能的。因此现有的方法,通常是根据地下水运动方向污染物浓度的变化,依据地下水运动方程和溶质迁移方程,进行反演。然而这种方法没有办法检验溶质的质量平衡,反演确定的弥散系数在很大程度上取决于测点的浓度。或者说是部分溶质在地下水中运动所形成的弥散度。此外,对于地下水,污染物在水势梯度较低的地方可能发生聚集,而依据部分溶质测定的弥散系数在这种情况下对于溶质迁移的真实性是值得怀疑的。目前多数方法还没有形成公认而成熟的规范,假设偏于理想化,求得的参数难以体现实际情况。
如上所述,一套完整的水动力弥散系数测定方法还没有规范起来,无法为制定合理有效的地下水污染防治措施,控制污染的发展,改善地下水水质环境,选择最佳治理方案提供科学依据。
发明内容
为了克服上述现有技术的不足,本发明提供了一种全过程通量守恒的现场土壤溶质运移弥散系数测定的方法,为制定合理有效地地下水污染防治措施,控制污染的发展,改善地下水水质环境提供支持参数支持。
本发明所提供的技术方案具体如下:
一种全过程通量守恒的现场土壤溶质运移弥散系数测定的方法,包括以下步骤:
(1)对地下水进行流场识别,确定地下水中的溶质在土壤中运动的主要和次要迁移方向;
(2)根据步骤(1)中确定的溶质的主要和次要迁移方向,基于溶质的通量守恒原理,建立模拟土壤溶质运动的数学模型;
(3)确定步骤(2)中的模拟土壤溶质运动的数学模型的现场实验条件;
(4)在控制边界条件下,以注入源的形式注入非转化型溶质,采用径向流动观测方式,测定非转化型溶质在主要和次要迁移方向上浓度随时间变化的情况;
(5)根据步骤(4)中的实验结果绘制弥散系数与孔隙平均流速的关系图并分析土壤和孔隙平均流速对弥散系数的影响,基于非转化型溶质全过程的通量守恒原理确定弥散系数。
所述的模拟土壤溶质运动的数学模型为:
u=Q/(2πrhn0) 式(2)
D=BVA 式(3)
式中:
c表示注入非转化型溶质浓度;t为时间坐标,u为水流的孔隙流速,x、y为方向坐标,αL、αT分别表示沿x方向和沿y方向的的动力弥散扩展度,u为水流的孔隙流速,r表示测点距径向观测孔距离;θ为径向观测孔测点与边界条件控制孔连线和水流轴线的夹角;Q为控制孔非转化型溶质注入强度,h为含水层厚度;n0为含水层有效孔隙率;D为弥散系数,V为平均孔隙流速,A、B为与介质特征有关的常数。
其中,式(1)是根据质量守恒推导出来的,式(2)是根据达西定律的原理推出来的,式(3)是自己根据实际情况设定的函数形式。
非转化型溶质在径向流动监测点的浓度为:
t′rmax=(a2+4)1/2-2 式(6)
tr'=t/tc 式(7)
其中,tr'为相对时刻,t为从非转化型溶质注入时间,tc为在纯对流情况下非转化型溶质由注入孔到达径向观测孔的时间,为tr'的前一个时刻,a为径向相对位置,cr(a,tr')为tr'时刻径向相对位置a的浓度,K为径向最大移动相对距离,t'rmax为达到径向最大相对移动距离的相对时刻。
式(4)~式(7)是为了根据实验数据确定参数而推导出来的公式。
步骤(5)所述的关系图为折线图。
与现有技术相比,本发明具有以下优点和有益效果:
(1)解决了现有方法以主要流动方向确定弥散系数的缺点,充分考虑了地下水的流动特点以及所确定的弥散系数在机理上具有完全的物理意义,有效地解决了弥散系数的确定很大程度上依赖测点位置选择的问题;
(2)现场条件下,建立易于操作、误差小以及便于测量的数学模型,模拟土壤溶质运动,实验条件针对研究现场,能够直接确定现场的弥散系数;
(3)采用径向流动观测方式,在主要和次要迁移方向监测注入的非转化型溶质浓度变化;
(4)弥散系数关系图采用折线图方式,直观反映土壤和孔隙流速对弥散系数的影响;
(5)实验结果误差小,分析方法准确、可靠。
具体实施方式
为了便于本领域普通技术人员理解和实施本发明,下面结合实施例对本发明作进一步的详细描述,应当理解,此处所描述的实施示例仅用于说明和解释本发明,并不用于限定本发明。
本发明全过程通量守恒的现场土壤溶质运移弥散系数测定的方法,包括以下步骤:
(1)对地下水进行流场识别,确定地下水中的溶质在土壤中运动的主要和次要迁移方向;
(2)根据步骤(1)中确定的溶质的主要和次要迁移方向,基于溶质的通量守恒原理,建立模拟土壤溶质运动的数学模型;
(3)确定步骤(2)中的模拟土壤溶质运动的数学模型的现场实验条件;
(4)在控制边界条件下,以注入源的形式注入非转化型溶质,采用径向流动观测方式,测定非转化型溶质在主要和次要迁移方向上浓度随时间变化的情况;
(5)根据步骤(4)中的实验结果绘制弥散系数与孔隙平均流速的关系图并分析土壤和孔隙平均流速对弥散系数的影响,基于非转化型溶质全过程的通量守恒原理确定弥散系数。
以下通过对饱和条件下弥散系数的测定来详细阐述本发明的技术方案,具体实施步骤如下:
1.在采用径向流动观测方式的情况下,饱和土壤水动力弥散方程为:
u=Q/(2πrhn0) 式(2)
式中:c为注入非转化型溶质浓度;t为时间坐标,r为测点距径向观测孔距离;θ为径向观测孔测点与边界条件控制孔连线和水流轴线的夹角;u为水流的孔隙流速,Q为控制孔非转化型溶质注入强度,h为含水层厚度;n0为含水层有效孔隙率。
推断弥散系数D与平均孔隙流速V的依赖关系式,即:
D=BVA 式(3)
其中:D为弥散系数,V为平均孔隙流速,A、B为与介质特征有关的常数,为表征多孔介质对于溶质分散能力的基本参数,要研究溶质在介质中的运移过程,须首先确定出这些参数。
2.在控制边界条件下,在现场注入非转化型溶质,采用径向流动观测方式,在主要和次要迁移方向监测注入的非转化型溶质浓度变化,测定非转化型溶质浓度在主要和次要移动方向随时间的变化。
3.根据控制条件,径向流动监测点非转化型溶质浓度为:
式中:
t′rmax=(a2+4)1/2-2 式(6)
tr'=t/tc 式(7)
其中,tr'为相对时刻,t为从非转化型溶质注入时间,tc为在纯对流情况下非转化型溶质由注入孔到达径向观测孔的时间,为tr'的前一个时刻,a为径向相对位置,cr(a,tr')为tr'时刻径向相对位置a的浓度,K为径向最大移动相对距离,t'rmax为达到径向最大相对移动距离的相对时刻。
4.根据(4)~(7)式,采用在主要和次要监测孔的非转化型浓度随时间的变化过程监测,推求满足全过程溶质质量均衡的弥散系数。
5.通过观察直观反映弥散系数与孔隙平均流速的关系图来分析土壤和孔隙平均流速对弥散系数的影响。
本实施例根据土壤溶质运移的水动力弥散方程推断弥散系数与平均孔隙流速的依赖关系式,建立数学模型,通过现场条件下的主要流动方向和次要流动方向,在控制边界条件下,采用径向观测孔对注入非转化型溶质浓度变化进行监测,推求满足全过程溶质均衡的弥散系数。绘制弥散系数与孔隙平均流速的关系图,分析弥散系数的影响因素。该方法易于操作、测量精度高,实验结果直观,且在该领域具有独创性。
应当理解的是,本说明书未详细阐述的部分均属于现有技术。
应当理解的是,上述针对较佳实施例的描述较为详细,并不能因此而认为是对本发明专利保护范围的限制,本领域的普通技术人员在本发明的启示下,在不脱离本发明权利要求所保护的范围情况下,还可以做出替换或变形,均落入本发明的保护范围之内,本发明的请求保护范围应以所附权利要求为准。
Claims (4)
1.一种全过程通量守恒的现场土壤溶质运移弥散系数测定的方法,其特征在于,包括以下步骤:
(1)对地下水进行流场识别,确定地下水中的溶质在土壤中运动的主要和次要迁移方向;
(2)根据步骤(1)中确定的溶质的主要和次要迁移方向,基于溶质的通量守恒原理,建立模拟土壤溶质运动的数学模型;
(3)确定步骤(2)中的模拟土壤溶质运动的数学模型的现场实验条件;
(4)在控制边界条件下,以注入源的形式注入非转化型溶质,采用径向流动观测方式,测定非转化型溶质在主要和次要迁移方向上浓度随时间变化的情况;
(5)根据步骤(4)中的实验结果绘制弥散系数与孔隙平均流速的关系图并分析土壤和孔隙平均流速对弥散系数的影响,基于非转化型溶质全过程的通量守恒原理确定弥散系数。
2.根据权利要求1所述的全过程通量守恒的现场土壤溶质运移弥散系数测定的方法,其特征在于:所述的模拟土壤溶质运动的数学模型为:
u=Q/(2πrhn0) 式(2)
D=BVA 式(3)
式中:
c表示注入非转化性溶质浓度,t为时间坐标,u为水流的孔隙流速,x、y为方向坐标,αL、αT分别表示沿x方向和沿y方向的的动力弥散扩展度,r表示测点距径向观测孔距离,Q为控制孔非转化性溶质注入强度,h为含水层厚度,n0为含水层有效孔隙率,D为弥散系数,V为平均孔隙流速,A、B为与介质特征有关的常数。
3.根据权利要求1所述的全过程通量守恒的现场土壤溶质运移弥散系数测定的方法,其特征在于:非转化型溶质在径向流动监测点的浓度为:
t′rmax=(a2+4)1/2-2 式(6)
tr'=t/tc 式(7)
其中,tr'为相对时刻,t为从非转化型溶质注入时间,tc为在纯对流情况下非转化型溶质由注入孔到达径向观测孔的时间,为tr'的前一个时刻,a为径向相对位置,cr(a,tr')为tr'时刻径向相对位置a的浓度,K为径向最大移动相对距离,t'rmax为达到径向最大相对移动距离的相对时刻。
4.根据权利要求1所述的全过程通量守恒的现场土壤溶质运移弥散系数测定的方法,其特征在于:步骤(5)所述的关系图为折线图。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510307246.8A CN104865165B (zh) | 2015-06-05 | 2015-06-05 | 全过程通量守恒的现场土壤溶质运移弥散系数测定的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510307246.8A CN104865165B (zh) | 2015-06-05 | 2015-06-05 | 全过程通量守恒的现场土壤溶质运移弥散系数测定的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104865165A true CN104865165A (zh) | 2015-08-26 |
CN104865165B CN104865165B (zh) | 2018-01-12 |
Family
ID=53911148
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510307246.8A Expired - Fee Related CN104865165B (zh) | 2015-06-05 | 2015-06-05 | 全过程通量守恒的现场土壤溶质运移弥散系数测定的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104865165B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105181531A (zh) * | 2015-08-29 | 2015-12-23 | 西安科技大学 | 黄土水分迁移规律室内模拟系统及特征参数测定方法 |
CN105181533A (zh) * | 2015-08-29 | 2015-12-23 | 西安科技大学 | 基于室内模拟试验的降雪入渗参数测定系统及方法 |
CN107167576A (zh) * | 2017-05-09 | 2017-09-15 | 武汉大学 | 非饱和土壤紊散动力学系数测定方法及装置 |
CN108680470A (zh) * | 2018-07-17 | 2018-10-19 | 武汉大学 | 微喷灌条件下土壤水入渗分布测试方法 |
CN108918388A (zh) * | 2018-07-18 | 2018-11-30 | 武汉大学 | 地下水含水层溶质弥散系数及孔隙速率测定方法 |
CN109297867A (zh) * | 2018-10-15 | 2019-02-01 | 成都理工大学 | 一种污染区现场快速检测管及检测方法 |
CN109668809A (zh) * | 2019-02-28 | 2019-04-23 | 水利部牧区水利科学研究所 | 河谷平原带溶质累积与迁移模拟实验装置及方法 |
CN109917103A (zh) * | 2019-03-13 | 2019-06-21 | 中国地质科学院水文地质环境地质研究所 | 一种原位土柱灌溉淋溶试验及溶质迁移定量描述方法 |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108614942B (zh) * | 2018-05-10 | 2020-05-05 | 河海大学 | 一种表征多孔介质中溶质运移时空尺度关联的方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101908100A (zh) * | 2010-07-26 | 2010-12-08 | 中国科学院生态环境研究中心 | 一种地下水环境的建模及数值模拟方法 |
CN102073796A (zh) * | 2011-02-21 | 2011-05-25 | 南京大学 | 一种模拟溶质三维运移过程的格子行走方法 |
CN104359802A (zh) * | 2014-11-20 | 2015-02-18 | 中国地质大学(武汉) | 基于一维水平流的多段式包气带水动力弥散系数测定仪 |
-
2015
- 2015-06-05 CN CN201510307246.8A patent/CN104865165B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101908100A (zh) * | 2010-07-26 | 2010-12-08 | 中国科学院生态环境研究中心 | 一种地下水环境的建模及数值模拟方法 |
CN102073796A (zh) * | 2011-02-21 | 2011-05-25 | 南京大学 | 一种模拟溶质三维运移过程的格子行走方法 |
CN104359802A (zh) * | 2014-11-20 | 2015-02-18 | 中国地质大学(武汉) | 基于一维水平流的多段式包气带水动力弥散系数测定仪 |
Non-Patent Citations (3)
Title |
---|
孙讷正: "地下水水质的数学模拟(三)——水动力弥散方程的解析解法及其应用", 《水文地质工程地质》 * |
杨武成 等: "非饱和带地下水溶质运移模拟理论研究", 《节水灌溉》 * |
杨金钟: "弥散系数的室内实验研究", 《工程勘察》 * |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105181531A (zh) * | 2015-08-29 | 2015-12-23 | 西安科技大学 | 黄土水分迁移规律室内模拟系统及特征参数测定方法 |
CN105181533A (zh) * | 2015-08-29 | 2015-12-23 | 西安科技大学 | 基于室内模拟试验的降雪入渗参数测定系统及方法 |
CN105181531B (zh) * | 2015-08-29 | 2016-05-11 | 西安科技大学 | 黄土水分迁移规律室内模拟系统及特征参数测定方法 |
CN107167576A (zh) * | 2017-05-09 | 2017-09-15 | 武汉大学 | 非饱和土壤紊散动力学系数测定方法及装置 |
CN107167576B (zh) * | 2017-05-09 | 2020-02-14 | 武汉大学 | 非饱和土壤紊散动力学系数测定方法及装置 |
CN108680470A (zh) * | 2018-07-17 | 2018-10-19 | 武汉大学 | 微喷灌条件下土壤水入渗分布测试方法 |
CN108918388A (zh) * | 2018-07-18 | 2018-11-30 | 武汉大学 | 地下水含水层溶质弥散系数及孔隙速率测定方法 |
CN109297867A (zh) * | 2018-10-15 | 2019-02-01 | 成都理工大学 | 一种污染区现场快速检测管及检测方法 |
CN109297867B (zh) * | 2018-10-15 | 2023-10-13 | 成都理工大学 | 一种污染区现场快速检测管及检测方法 |
CN109668809A (zh) * | 2019-02-28 | 2019-04-23 | 水利部牧区水利科学研究所 | 河谷平原带溶质累积与迁移模拟实验装置及方法 |
CN109917103A (zh) * | 2019-03-13 | 2019-06-21 | 中国地质科学院水文地质环境地质研究所 | 一种原位土柱灌溉淋溶试验及溶质迁移定量描述方法 |
CN109917103B (zh) * | 2019-03-13 | 2021-09-07 | 中国地质科学院水文地质环境地质研究所 | 一种原位土柱灌溉淋溶试验及溶质迁移定量描述方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104865165B (zh) | 2018-01-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104865165A (zh) | 全过程通量守恒的现场土壤溶质运移弥散系数测定的方法 | |
Bouma et al. | Water infiltration and redistribution in a silt loam subsoil with vertical worm channels | |
Heiss et al. | Coupled surface‐subsurface hydrologic measurements reveal infiltration, recharge, and discharge dynamics across the swash zone of a sandy beach | |
Lewandowski et al. | A heat pulse technique for the determination of small‐scale flow directions and flow velocities in the streambed of sand‐bed streams | |
Hamm et al. | Relationship between transmissivity and specific capacity in the volcanic aquifers of Jeju Island, Korea | |
Abraham et al. | Validation of bed-load transport measurements with time-sequenced bathymetric data | |
Xu et al. | Analysis of single-ring infiltrometer data for soil hydraulic properties estimation: Comparison of BEST and Wu methods | |
Ata Amini et al. | A local scour prediction method for pile caps in complex piers | |
Lassen et al. | Monitoring CO2 gas-phase migration in a shallow sand aquifer using cross-borehole ground penetrating radar | |
CN104652347A (zh) | 山区非静态水体水位与淹没影响人口关系评价方法 | |
Cheng et al. | Water infiltration underneath single-ring permeameters and hydraulic conductivity determination | |
Cho et al. | Estimating natural gas emissions from underground pipelines using surface concentration measurements☆ | |
Schulz et al. | Monitoring of a simulated CO2 leakage in a shallow aquifer using stable carbon isotopes | |
Su et al. | Using multiple environmental methods to estimate groundwater discharge into an arid lake (Dakebo Lake, Inner Mongolia, China) | |
Nicosia et al. | Testing a theoretical resistance law for overland flow on a stony hillslope | |
Moret-Fernández et al. | Sequential infiltration analysis of infiltration curves measured with disc infiltrometer in layered soils | |
Caputo et al. | Measurement of field-saturated hydraulic conductivity on fractured rock outcrops near Altamura (Southern Italy) with an adjustable large ring infiltrometer | |
Hird et al. | Migration of sodium chloride in dry porous materials | |
Baldock et al. | Numerical study of the flow structure at a swash tip propagating over a rough bed | |
Dong et al. | Modelling soil solute release and transport in run‐off on a loessial slope with and without surface stones | |
Fisher et al. | Heat flow, sediment and pore fluid chemistry, and hydrothermal circulation on the east flank of Alarcon Ridge, Gulf of California | |
Eggleston et al. | Can we predict subsurface mass transport? | |
Külls et al. | Channel infiltration study using dye tracers | |
Charbeneau et al. | Physical modeling of sheet flow on rough impervious surfaces | |
CN113295590A (zh) | 一种改进的确保内环水分垂直入渗的双环渗透试验方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
EXSB | Decision made by sipo to initiate 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: 20180112 Termination date: 20180605 |
|
CF01 | Termination of patent right due to non-payment of annual fee |