CN113484210B - 一种强风化层弥散度现场尺度试验测定方法 - Google Patents
一种强风化层弥散度现场尺度试验测定方法 Download PDFInfo
- Publication number
- CN113484210B CN113484210B CN202110593679.XA CN202110593679A CN113484210B CN 113484210 B CN113484210 B CN 113484210B CN 202110593679 A CN202110593679 A CN 202110593679A CN 113484210 B CN113484210 B CN 113484210B
- Authority
- CN
- China
- Prior art keywords
- tracer
- underground water
- conductivity
- test
- boundary
- 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
- 238000012360 testing method Methods 0.000 title claims abstract description 67
- 238000000034 method Methods 0.000 title claims abstract description 39
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 118
- 239000000700 radioactive tracer Substances 0.000 claims abstract description 71
- 238000012544 monitoring process Methods 0.000 claims abstract description 46
- 239000006185 dispersion Substances 0.000 claims abstract description 25
- FAPWRFPIFSIZLT-UHFFFAOYSA-M sodium chloride Inorganic materials [Na+].[Cl-] FAPWRFPIFSIZLT-UHFFFAOYSA-M 0.000 claims description 44
- VEXZGXHMUGYJMC-UHFFFAOYSA-M Chloride anion Chemical compound [Cl-] VEXZGXHMUGYJMC-UHFFFAOYSA-M 0.000 claims description 26
- 239000003673 groundwater Substances 0.000 claims description 25
- 239000000523 sample Substances 0.000 claims description 24
- 239000011780 sodium chloride Substances 0.000 claims description 23
- SQGYOTSLMSWVJD-UHFFFAOYSA-N silver(1+) nitrate Chemical compound [Ag+].[O-]N(=O)=O SQGYOTSLMSWVJD-UHFFFAOYSA-N 0.000 claims description 8
- 238000013508 migration Methods 0.000 claims description 7
- 230000005012 migration Effects 0.000 claims description 7
- 238000005553 drilling Methods 0.000 claims description 6
- 230000035699 permeability Effects 0.000 claims description 5
- 230000000694 effects Effects 0.000 claims description 4
- 238000002474 experimental method Methods 0.000 claims description 4
- 229910001961 silver nitrate Inorganic materials 0.000 claims description 4
- 238000004448 titration Methods 0.000 claims description 4
- 238000005086 pumping Methods 0.000 claims description 3
- 239000007787 solid Substances 0.000 claims description 3
- 239000003795 chemical substances by application Substances 0.000 claims description 2
- 238000011160 research Methods 0.000 abstract description 3
- 238000013178 mathematical model Methods 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 238000002347 injection Methods 0.000 description 3
- 239000007924 injection Substances 0.000 description 3
- 239000003153 chemical reaction reagent Substances 0.000 description 2
- 238000012806 monitoring device Methods 0.000 description 2
- 238000012805 post-processing Methods 0.000 description 2
- 239000004576 sand Substances 0.000 description 2
- 239000002689 soil Substances 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 238000010200 validation analysis Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/08—Investigating permeability, pore-volume, or surface area of porous materials
-
- E—FIXED CONSTRUCTIONS
- E02—HYDRAULIC ENGINEERING; FOUNDATIONS; SOIL SHIFTING
- E02D—FOUNDATIONS; EXCAVATIONS; EMBANKMENTS; UNDERGROUND OR UNDERWATER STRUCTURES
- E02D1/00—Investigation of foundation soil in situ
- E02D1/02—Investigation of foundation soil in situ before construction work
- E02D1/027—Investigation of foundation soil in situ before construction work by investigating properties relating to fluids in the soil, e.g. pore-water pressure, permeability
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N13/00—Investigating surface or boundary effects, e.g. wetting power; Investigating diffusion effects; Analysing materials by determining surface, boundary, or diffusion effects
- G01N13/04—Investigating osmotic effects
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N27/00—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means
- G01N27/02—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance
- G01N27/04—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance by investigating resistance
- G01N27/06—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance by investigating resistance of a liquid
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N31/00—Investigating or analysing non-biological materials by the use of the chemical methods specified in the subgroup; Apparatus specially adapted for such methods
- G01N31/16—Investigating or analysing non-biological materials by the use of the chemical methods specified in the subgroup; Apparatus specially adapted for such methods using titration
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Health & Medical Sciences (AREA)
- Analytical Chemistry (AREA)
- Pathology (AREA)
- Physics & Mathematics (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Engineering & Computer Science (AREA)
- Dispersion Chemistry (AREA)
- Electrochemistry (AREA)
- Molecular Biology (AREA)
- Hydrology & Water Resources (AREA)
- Soil Sciences (AREA)
- Chemical Kinetics & Catalysis (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Mining & Mineral Resources (AREA)
- Paleontology (AREA)
- Civil Engineering (AREA)
- General Engineering & Computer Science (AREA)
- Structural Engineering (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种强风化层弥散度现场尺度试验测定方法，属于水文地质技术领域，在野外现场构建含有示踪剂投放源以及若干观测孔的试验场地；通过将电导率和水压力传感器布设在观测孔中，利用标定的研究区地下水体电导率与示踪剂浓度的相关曲线将监测所得电导率数据转换为示踪剂浓度，利用GMS软件搭建与试验平台对应的数值模型，通过拟合示踪剂浓度‑时间序列数据获取现场尺度下强风化层弥散度。本发明综合搭建了强风化层弥散度现场尺度试验平台和相应数值模型，实时观测孔地下水中示踪剂浓度变化情况，并通过数值模型拟合获取强风化层弥散度，精度较高。
Description
技术领域
本发明属于水文地质技术领域，一种强风化层弥散度现场尺度试验测定方法。
背景技术
弥散度是表征地下水中溶质迁移能力的重要指标，是预测地下水污染情况的重要参数。目前传统测定方法为土柱法或砂槽法，土柱法只能测定一维流的弥散度，砂槽法只能测定一维流下二维弥散度，与野外实际情况相差很大，且也不能反映溶质运移过程中弥散度在各个方向的关系。
目前，强风化层弥散度的确定方法尚无统一标准，且现场尺度的确定方法尚未见报导。因此，强风化层现场尺度弥散度试验测定方法具有重要的理论和实践意义。
发明内容
发明目的：本发明的目的在于提供一种强风化层弥散度现场尺度试验测定方法，可以实时监测现场所研究强风化层中示踪剂运移情况，并通过数值模型拟合获取强风化层中的弥散度，精度较高且较为方便。
技术方案：为实现上述目的，本发明提供如下技术方案：一种强风化层弥散度现场尺度试验测定方法，包括以下步骤：
1)利用现场钻孔和地下水位观测确定现场地下水主流向；
2)选择地表平整的区域作为试验场地，在强风化层中沿地下水主流向布置示踪剂投放孔和观测孔，另在地下水主流方向的侧向布置观测孔；测试钻孔和观测孔地下水水位；
3)采取现场地下水体，用其配置不同浓度的示踪剂-氯化钠样品，利用电导率仪测试氯化钠样品的电导率值，并利用硝酸银滴定法测定氯化钠样品中氯离子浓度，建立电导率值和氯离子浓度标准曲线；
4)将监测探头布设在示踪剂投放孔和监测孔中，确保监测孔中有一个监测探头位于初始地下水位以下，将监测探头与数据采集器连接；
5)配置固定浓度的示踪剂置于容器中，实验开始后将示踪剂泵入示踪剂投放孔中；利用所述监测探头连续监测水体的电导率和水压力；利用电导率值和氯离子浓度标准曲线将地下水体中的电导率监测数据转换为氯离子浓度数据；
6)示踪剂投放孔中的示踪剂沿着地下水流向通过对流弥散向前运移；描述地下水流动过程的为非均值各向异性等效连续介质地下水非稳定流方程，描述地下水中示踪剂运移的方程为对流弥散方程，具体如下：非均值各向异性等效连续介质地下水非稳定流方程为：
式中：Kx，Ky，Kz，为x，y，z方向渗透系数，m/d；H为地下水水头，m；W为单位时间从单位体积含水层流入或流出的水量，1/d；μ为给水度，无量纲；t-时间，d；Ω为渗流区域；B1为水头已知边界，第一类边界；H(x,y,z)为坐标为(x,y,z)处的地下水水头；H1(x,y,z,t)，任意大于零时刻，坐标(x,y,z)处的含水层边界地下水水头为H1，m；B2为流量已知边界，第二类边界；K为渗透系数，m/d；n是法线方向向量；q(x,y,z,t)为任意大于零时刻，坐标(x,y,z)处的含水层第二类边界单位面积过水断面补给流量为q，m/d；
示踪剂在地下水流中的对流弥散方程如下：
式中：R-阻滞系数；θ-介质孔隙度；C-地下水中示踪剂质量浓度，g/L；t-时间，d；D-水动力弥散系数，m2/d；α-弥散度，单位：m；ν-地下水渗流速度，m/d。
该数学模型难以用解析方法求解，可用数值法求解，采用有限差分数值法求解该数学模型，预测示踪剂从投放孔运移至观测孔的过程。建立与现场试验对应数值模型并拟合试验参数：利用GMS(Groundwater Modeling System)软件搭建与试验平台对应的数值模型，按照试验的物理边界设置数值模型边界，包括试验场地四周的水文地质边界条件、地下水初始流场和初始氯离子浓度。考虑示踪剂在地下水中的对流与弥散作用，通过拟合试验数据获取强风化层现场尺度弥散度。
进一步的：所述的步骤1)中，在现场利用至少三口井确定地下水主流方向，井的深度需要穿透目标层-强风化层，位于强风化层的井壁上设置孔眼，作为水井的花管过滤器。
进一步的：所述的步骤2)中，投放孔和观测孔的深度穿透强风化层到达弱风化层，侧面观测孔与主流向的夹角应不大于15°，距离投放孔直线距离2m处。
进一步的：所述的步骤3)中，用于标定标准曲线的氯化钠样品数应不少于五组，电导率值和氯离子浓度标准曲线的决定系数不低于0.99。
进一步的：所述的步骤4)中，孔内监测探头有一个位于初始地下水位10cm以下，另外一个位于强风化层中，自动记录电导率/水压力值的频率不低于10次/小时。
进一步的：所述的步骤5)中，应先往投放孔中投入固体氯化钠以使得孔中地下水中的电导率迅速上升至配置的示踪剂电导率，试验期间确保投放孔中的水位稳定。
进一步的：所述的步骤6)中，示踪剂浓度的单位要与步骤5)中的一致，拟合的相关系数不低于0.9。
发明原理：本发明通过在野外现场构建含有示踪剂投放源以及若干观测孔的试验场地；通过将电导率和水压力传感器布设在观测孔中，利用标定的研究区地下水体电导率与示踪剂浓度的相关曲线将监测所得电导率数据转换为示踪剂浓度，利用GMSGroundwater Modeling System)软件搭建与试验平台对应的数值模型，通过拟合示踪剂浓度-时间序列数据获取现场尺度下强风化层弥散度。本发明综合搭建了强风化层弥散度现场尺度试验平台和相应数值模型，实时观测孔地下水中示踪剂浓度变化情况，并通过数值模型拟合获取强风化层弥散度，精度较高。
有益效果：与现有技术相比，本发明的一种强风化层弥散度现场尺度试验测定方法，属于现场尺度方法，较之试验室尺度更贴近实际，所得系数更为精确；试验的示踪剂一般选择氯化钠，为环境友好型试剂，且制备简单，成本低廉；试验采用的监测装置可实现连续无扰动长时间监测，数据的后处理简单精确；利用GMS(Groundwater Modeling System)软件搭建与试验平台对应的数值模型，可以契合试验的各物理边界条件，同时拟合研究强风化地层中的溶质运移弥散度，拟合结果精确。
附图说明
图1利用三点法获得地下水主流向示意图；
图2试验区钻孔及监测孔布平面位置图；
图3示踪试验工作布置剖面图；
图4电导率与氯化钠浓度标准曲线；
图5实例数值模型网格及边界条件设置示意图；
图6实例区地下水渗流模型验证结果；
图7示踪试验拟合及弥散度获取。
具体实施方式
下面结合附图和具体实施例对本发明作更进一步的说明。
本发明提出一种强风化层弥散度现场尺度试验测定方法，包括以下步骤：
1)利用现场钻孔和地下水位观测确定现场地下水主流向；
2)选择地表平整的区域作为试验场地，在强风化层中沿地下水主流向布置示踪剂投放孔和观测孔，另在地下水主流方向的侧向布置观测孔；示踪剂投放孔和监测孔设置花管；
3)测试地下水背景电导率和氯离子浓度，测试钻孔和观测孔地下水水位；配置固定浓度的示踪剂置于容器中，实验开始后将示踪剂泵入示踪剂投放孔中；
4)搭建现场监测系统：将监测探头(连续监测水体的电导率和水压力)布设在示踪剂投放孔和监测孔中，确保监测孔中有一个监测探头位于初始地下水位以下，记录地下水中的电导率和水压力变化情况；将监测探头与数据采集器连接；
5)标定地下水体中的电导率与氯离子浓度标准曲线：采取现场地下水体，配置不同浓度的氯化钠样品，利用电导率仪测试其电导率值，利用硝酸银滴定法测定其氯离子浓度，建立电导率值和氯离子浓度标准曲线；利用上述标准曲线将地下水体中的电导率监测数据转换为氯离子浓度数据；
6)示踪剂投放孔中的示踪剂会沿着地下水流向通过对流弥散向前运移；描述地下水流动过程的为非均值各向异性等效连续介质地下水非稳定流方程，描述地下水中示踪剂运移的方程为对流弥散方程，具体如下：
非均值各向异性等效连续介质地下水非稳定流方程为：
式中：Kx，Ky，Kz，为x，y，z方向渗透系数，m/d；H为地下水水头，m；W为单位时间从单位体积含水层流入或流出的水量，1/d；μ为给水度；t-时间；Ω为渗流区域；B1为水头已知边界，第一类边界；H1为各层边界水位，m；B2为流量已知边界，第二类边界；q为含水层第二类边界单位面积过水断面补给流量，m/d。
示踪剂在三维地下水流中的对流弥散方程如下：
式中：R-阻滞系数；θ-介质孔隙度；C-地下水中组分质量浓度，g/L；t-时间，d；Di-水动力弥散系数，m2/d；α-弥散度，单位：m；νi-地下水渗流速度，m/d。
该数学模型难以用解析方法求解，可用数值法求解，采用有限差分数值法求解该数学模型，预测示踪剂从投放孔运移至观测孔的过程。本例建立对应数值模型并拟合试验参数：利用GMS(Groundwater Modeling System)软件搭建与试验平台对应的数值模型，按照试验的物理边界设置数值模型边界，包括试验场地四周的水文地质边界条件、地下水初始流场和初始氯离子浓度。考虑示踪剂在地下水中的对流与弥散作用，通过拟合试验数据获取强风化层现场尺度弥散度。
步骤1)，在现场利用至少三口井确定地下水主流方向，井的深度要穿透目标层-强风化层。步骤2)，投放孔和观测孔的深度穿透强风化层到达弱风化层，孔的花管应至少包含强风化层；侧面观测孔与主流向的夹角应不大于15°。步骤3)中，应测试枯水期和丰水期钻孔和观测孔中的地下水电导率、氯离子浓度和地下水水位，分别作为枯水期及丰水期的背景；应先往投放孔中投入固体氯化钠以使得孔中地下水中的电导率迅速上升至配置的示踪剂电导率；试验期间确保投放孔中的水位稳定。步骤4)中，孔内监测探头至少有一个位于初始地下水位10cm以下，另外一个位于强风化层中，自动记录电导率/水压力值的频率不低于10次/小时。步骤5)，地下水示踪剂样品数应不少于五组，电导率值和氯离子浓度标准曲线的决定系数不低于0.99。步骤6)中，示踪剂浓度的单位要与步骤5)中的一致，拟合的相关系数不低于0.9。
实施例
一种强风化层弥散度现场尺度试验测定方法，具体操作如下：
(1)利用现场钻孔ZK1、ZK2和ZK3，使用三点法确定现场地下水主流向，首先，测出三口钻孔的地下水位高程；然后三点连线作一个三角形，在ZK1和ZK2连线上内插出ZK3地下水位高程(7m)所在点ZK3’，连接ZK3-ZK3’即是等水位线，垂直等水位线方向由地下水位高程从高到低即是地下水主流向，可见图1。
(2)在现场选择地表平整的12m×18m的区域作为试验场地，在强风化层中沿地下水主流向布置示踪剂投放孔ZK1和观测孔JC 1、JC 3和JC 4，另在地下水主流方向的侧向夹角为15°方向上布置观测孔JC 2；JC 1、JC 3和JC 4距离ZK1的距离分别为1.5m、3m和6m，JC2距离ZK1的距离为2.0m；示踪剂投放孔和监测孔的深度到达微风化层，其中ZK1孔深30米，花管深度埋深5～35m，观测孔JC 1、JC 2、JC 3和JC 4孔深25米，花管深度埋深5～25m，可见图2和图3。
(3)测试地下水背景电导率和氯离子浓度，测试钻孔和观测孔地下水水位；配置固定浓度(20g/L)的示踪剂置于容器中，实验开始后先投氯化钠示踪剂11千克，使得孔中地下水中的氯化钠浓度迅速达到20g/L，接着持续泵入示踪剂至投放孔ZK 1中，利用溢流槽保持ZK 1的水位保持在孔口高程。
(4)搭建现场监测系统
将监测探头(连续监测水体的电导率和水压力)布设在示踪剂投放孔ZK 1和监测孔JC 1、JC 2、JC 3和JC 4中，确保监测孔中有一个监测探头位于初始地下水位以下，记录地下水中的电导率和水压力变化情况；将监测探头与数据采集器连接，可见图3；
(5)标定地下水体中的电导率与氯离子浓度标准曲线：采取现场地下水体，配置不同浓度的氯化钠样品，利用电导率仪测试其电导率值，利用硝酸银滴定法测定其氯离子浓度，建立电导率值和氯离子浓度标准曲线，试验所得曲线关系为：氯化钠浓度＝0.5848×电导率值，单位分别为g/L和ms/cm，决定系数为0.9995，利用上述标准曲线将地下水体中的电导率监测数据转换为氯离子浓度数据，可见图4；
(6)建立对应数值模型并拟合试验参数：
利用GMS(Groundwater Modeling System)软件搭建与试验平台对应的数值模型，按照试验的物理边界设置数值模型边界，包括试验场地四周的水文地质边界条件、地下水初始流场和初始氯离子浓度。考虑示踪剂在地下水中的对流与弥散作用，通过拟合试验数据获取强风化层现场尺度弥散度。利用GMS建立以试验区为中心的网格模型，模型的尺寸为180m×100m×70m，场地东西两侧设置为的定水头边界，模型分为6层，1-2层为覆盖层与强风化层，共厚18m，2-4层为中等风化层，共厚6m，5-6层为微风化岩体层，共厚46m。模型的网格按照钻孔所在位置进行加密，共有有效网格121632个，可见图5；
利用GMS建立了示踪剂持续注入35天的地下水渗流非稳定流模型，由模型的计算结果与JC01水位观测值可知，模型计算的水位较为准确，误差小于0.5m，可见图6；
接着建立ZK1持续注入浓度为20g/L示踪剂，连续注入35天后停止注入，继续监测至第60天的示踪剂运移数值模型。模型中的孔隙度分别设置为0.1(1-2层)、0.05(3-4层)和0.01(5-6层)。分别输入不同的弥散度α(单位：m)来模拟示踪剂的运移过程，，并以JC01井中监测的示踪剂浓度变化情况为例进行拟合分析。本次拟合分别计算了弥散度为0.01m～6m范围共20个工况，选取了其中决定系数较高的四种工况，即弥散度分别为0.3m、0.5m、1.0m和2.0m的情况，其中弥散度α＝0.5的决定系数最高为0.94，因此，实例强风化层弥散度建议值为0.5m，结果见图7。
本发明的一种强风化层弥散度现场尺度试验测定方法，属于现场尺度方法，较之试验室尺度更贴近实际，所得系数更为精确；试验的示踪剂一般选择氯化钠，为环境友好型试剂，且制备简单，成本低廉；试验采用的监测装置可实现连续无扰动长时间监测，数据的后处理精确；利用GMS(Groundwater Modeling System)软件搭建与试验平台对应的数值模型，可以契合试验的各物理边界条件，同时拟合研究强风化地层中的溶质运移弥散度，拟合结果精确。
Claims (4)
1.一种强风化层弥散度现场尺度试验测定方法，其特征在于，包括以下步骤：
1)利用现场钻孔和地下水位观测确定现场地下水主流向；
2)选择地表平整的区域作为试验场地，在强风化层中沿地下水主流向布置示踪剂投放孔和观测孔，另在地下水主流方向的侧向布置观测孔；测试钻孔和观测孔地下水水位；
3)采取现场地下水体，用其配置不同浓度的示踪剂-氯化钠样品，利用电导率仪测试氯化钠样品的电导率值，并利用硝酸银滴定法测定氯化钠样品中氯离子浓度，建立电导率值和氯离子浓度标准曲线；
4)将监测探头布设在示踪剂投放孔和监测孔中，确保监测孔中有一个监测探头位于初始地下水位以下，将监测探头与数据采集器连接；
5)配置固定浓度的示踪剂置于容器中，实验开始后将示踪剂泵入示踪剂投放孔中；利用所述监测探头连续监测水体的电导率和水压力；利用电导率值和氯离子浓度标准曲线将地下水体中的电导率监测数据转换为氯离子浓度数据；
6)示踪剂投放孔中的示踪剂沿着地下水流向通过对流弥散向前运移；描述地下水流动过程的为非均值各向异性等效连续介质地下水非稳定流方程，描述地下水中示踪剂运移的方程为对流弥散方程，具体如下：非均值各向异性等效连续介质地下水非稳定流方程为：
式中：Kx，Ky，Kz，为x，y，z方向渗透系数，m/d；H为地下水水头，m；W为单位时间从单位体积含水层流入或流出的水量，1/d；μ为给水度，无量纲；t-时间，d；Ω为渗流区域；B1为水头已知边界，第一类边界；H(x,y,z)为坐标为(x,y,z)处的地下水水头；H1(x,y,z,t)，任意大于零时刻，坐标(x,y,z)处的含水层边界地下水水头为H1，m；B2为流量已知边界，第二类边界；K为渗透系数，m/d；n是法线方向向量；q(x,y,z,t)为任意大于零时刻，坐标(x,y,z)处的含水层第二类边界单位面积过水断面补给流量为q，m/d；
示踪剂在地下水流中的对流弥散方程如下：
式中：R-阻滞系数；θ-介质孔隙度；C-地下水中示踪剂质量浓度，g/L；t-时间，d；D-水动力弥散系数，m2/d；α-弥散度，单位：m；ν-地下水渗流速度，m/d；
建立与现场试验对应数值模型并拟合试验参数：利用GMS(Groundwater ModelingSystem)软件搭建与试验平台对应的数值模型，按照试验的物理边界设置数值模型边界，包括试验场地四周的水文地质边界条件、地下水初始流场和初始氯离子浓度，考虑示踪剂在地下水中的对流与弥散作用，通过拟合试验数据获取强风化层现场尺度弥散度；
所述的步骤2)中，投放孔和观测孔的深度穿透强风化层到达弱风化层，侧面观测孔与主流向的夹角应不大于15°，距离投放孔直线距离2m处；
所述的步骤3)中，用于标定标准曲线的氯化钠样品数应不少于五组，电导率值和氯离子浓度标准曲线的决定系数不低于0.99；
所述的步骤4)中，孔内监测探头有一个位于初始地下水位10cm以下，另外一个位于强风化层中，自动记录电导率/水压力值的频率不低于10次/小时。
2.根据权利要求1所述的一种强风化层弥散度现场尺度试验测定方法，其特征在于：所述的步骤1)中，在现场利用至少三口井确定地下水主流方向，井的深度穿透目标层-强风化层，位于强风化层的井壁上设置孔眼，作为水井的花管过滤器。
3.根据权利要求1所述的一种强风化层弥散度现场尺度试验测定方法，其特征在于：所述的步骤5)中，应先往投放孔中投入固体氯化钠以使得孔中地下水中的电导率迅速上升至配置的示踪剂电导率，试验期间确保投放孔中的水位稳定。
4.根据权利要求1所述的一种强风化层弥散度现场尺度试验测定方法，其特征在于：所述的步骤6)中，示踪剂浓度的单位要与步骤5)中的一致，拟合的相关系数不低于0.9。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110593679.XA CN113484210B (zh) | 2021-05-28 | 2021-05-28 | 一种强风化层弥散度现场尺度试验测定方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110593679.XA CN113484210B (zh) | 2021-05-28 | 2021-05-28 | 一种强风化层弥散度现场尺度试验测定方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113484210A CN113484210A (zh) | 2021-10-08 |
| CN113484210B true CN113484210B (zh) | 2022-11-18 |
Family
ID=77933520
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110593679.XA Active CN113484210B (zh) | 2021-05-28 | 2021-05-28 | 一种强风化层弥散度现场尺度试验测定方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113484210B (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114859074B (zh) * | 2022-04-25 | 2024-04-16 | 无锡水文工程地质勘察院有限责任公司 | 一种地下水流向测量方法 |
| CN117610319B (zh) * | 2024-01-23 | 2024-04-30 | 南京大学 | 融合浓度和电导率的渗透系数与弥散度识别方法、系统 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5905036A (en) * | 1995-01-23 | 1999-05-18 | Board Of Regents, The University Of Texas System | Characterization of organic contaminants and assessment of remediation performance in subsurface formations |
| CN104794298A (zh) * | 2015-04-29 | 2015-07-22 | 南华大学 | 一种铀尾矿库核素对库区地下水污染风险评价的方法 |
| CN106600035A (zh) * | 2016-11-08 | 2017-04-26 | 北京师范大学 | 基于污染物迁移模拟的水源地水质安全预警方法 |
| CN108918388A (zh) * | 2018-07-18 | 2018-11-30 | 武汉大学 | 地下水含水层溶质弥散系数及孔隙速率测定方法 |
| CN110851991A (zh) * | 2019-11-18 | 2020-02-28 | 核工业二〇八大队 | 一种地下水流数值模拟方法 |
| AU2020102747A4 (en) * | 2020-10-16 | 2020-12-03 | North China Electric Power University | A decision-making method for in-situ remediation of petroleum-contaminated groundwater |
Family Cites Families (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103197039A (zh) * | 2013-03-15 | 2013-07-10 | 陈秋阳 | 一种人工流场下的野外地下水弥散试验方法 |
| CN103226732B (zh) * | 2013-03-21 | 2016-06-22 | 中南大学 | 一种基于gms的矿区不同开采中段的地下水渗流场预测方法 |
| CN103514489A (zh) * | 2013-10-09 | 2014-01-15 | 天津科技大学 | 正交设计和gms耦合的地下水污染物水力截获优化方法 |
| CN108614942B (zh) * | 2018-05-10 | 2020-05-05 | 河海大学 | 一种表征多孔介质中溶质运移时空尺度关联的方法 |
| CN110346247B (zh) * | 2019-07-29 | 2021-04-06 | 河海大学 | 一种吹填区底泥盐分扩散系数中尺度试验测定方法 |
| CN110608978B (zh) * | 2019-09-27 | 2024-04-05 | 中国地质大学（武汉） | 一种模拟含水层抽注水过程中细颗粒运移示踪试验装置 |
| CN111914447B (zh) * | 2020-07-13 | 2022-09-20 | 河海大学 | 模拟地下水溶质运移的新型有限体积多尺度有限元方法 |
| CN112507279A (zh) * | 2020-12-17 | 2021-03-16 | 天津科技大学 | 一种探究浅层咸水越流补给深层淡水的方法 |
-
2021
- 2021-05-28 CN CN202110593679.XA patent/CN113484210B/zh active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5905036A (en) * | 1995-01-23 | 1999-05-18 | Board Of Regents, The University Of Texas System | Characterization of organic contaminants and assessment of remediation performance in subsurface formations |
| CN104794298A (zh) * | 2015-04-29 | 2015-07-22 | 南华大学 | 一种铀尾矿库核素对库区地下水污染风险评价的方法 |
| CN106600035A (zh) * | 2016-11-08 | 2017-04-26 | 北京师范大学 | 基于污染物迁移模拟的水源地水质安全预警方法 |
| CN108918388A (zh) * | 2018-07-18 | 2018-11-30 | 武汉大学 | 地下水含水层溶质弥散系数及孔隙速率测定方法 |
| CN110851991A (zh) * | 2019-11-18 | 2020-02-28 | 核工业二〇八大队 | 一种地下水流数值模拟方法 |
| AU2020102747A4 (en) * | 2020-10-16 | 2020-12-03 | North China Electric Power University | A decision-making method for in-situ remediation of petroleum-contaminated groundwater |
Non-Patent Citations (5)
| Title |
|---|
| Differences in the scale dependence of dispersivity and retardation factors estimated from forced‐gradient and uniform flow tracer tests in three‐dimensional physically and chemically heterogeneous porous media;Fernàndez;《Water resources research》;20050316;第41卷(第3期);1-13 * |
| Influence of Thick Karst Vadose Zone on Aquifer Recharge in Karst Formations;Ruitong Liu;《Journal of Hydrology》;20210131;第592卷;3-8 * |
| 宏观弥散度和阻滞系数对地下水中核素迁移模拟的影响;谢水波等;《湖南大学学报(自然科学版)》;20070525(第05期);78-83 * |
| 岩溶区某磷石膏堆场地下水渗流与污染物运移数值模拟;张科正等;《中国煤炭地质》;20180525(第05期);46-48 * |
| 某岩溶区固废渣场地下水污染模拟;彭传路等;《中国煤炭地质》;20190930;第31卷(第9期);81-87 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113484210A (zh) | 2021-10-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Molz et al. | The impeller meter for measuring aquifer permeability variations: evaluation and comparison with other tests | |
| CN106437844B (zh) | 一种超前预报隧道涌水位置的方法 | |
| AU2009227997A1 (en) | System, apparatus and method for measuring soil moisture content | |
| CN113484210B (zh) | 一种强风化层弥散度现场尺度试验测定方法 | |
| CN106338457A (zh) | 一种模拟不同填充物裂隙的动水注浆试验装置 | |
| CN110793964A (zh) | 一种降雨诱发土质滑坡的野外模拟试验系统 | |
| Boaga et al. | Monitoring soil-plant interactions in an apple orchard using 3D electrical resistivity tomography | |
| CN110346533B (zh) | 用于指导盐碱地暗管布置组合的方法 | |
| CN110865008B (zh) | 基于有限尺度下圆形定水头边界非稳定流抽水的水文地质参数测定方法 | |
| CN109061105A (zh) | 一种土壤盐渍化临界地下水埋深的计算方法 | |
| CN107328909A (zh) | 结构性差异非饱和土壤水动力弥散系数现场测定方法 | |
| Lu et al. | Imaging and characterization of the preferential flow process in agricultural land by using electrical resistivity tomography and dual-porosity model | |
| CN108303512B (zh) | 一种原位测试土水特征曲线的方法 | |
| CN111766189B (zh) | 一种基于水力刺激的堤防隐伏渗漏通道三维层析扫描方法 | |
| CN106680454A (zh) | 一种具拦沙坝已治理崩岗土壤侵蚀模数测算方法 | |
| CN115979900A (zh) | 一种基于北方全排型岩溶泉的地下水循环要素监测方法 | |
| CN112650970B (zh) | 饱和-非饱和水分及溶质运移双重迭代耦合方法及装置 | |
| CN117408178A (zh) | 一种场地尺度温度-水动力耦合影响下污染物迁移关键参数获取方法 | |
| CN113901546B (zh) | 隧道涌水量计算方法及装置 | |
| CN110346247B (zh) | 一种吹填区底泥盐分扩散系数中尺度试验测定方法 | |
| CN115656468A (zh) | 一种基于示踪的原位非破坏地下水净补给量的测量方法 | |
| CN108982317B (zh) | 一种高填方土体大尺度渗透系数的确定方法 | |
| Aparicio et al. | TRANSPORT OF BROMIDE UNDER FIVE STEADY STATE WATER FLOW RATES AND THREE SOIL DEPTH: FIELD EXPERIMENT IN MOLISOL SOIL | |
| JP2002285535A (ja) | 水平透水係数の鉛直分布測定方法 | |
| JPH0355310A (ja) | 割れ目系岩盤の透水性試験方法 |
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 |