CN107167576A - 非饱和土壤紊散动力学系数测定方法及装置 - Google Patents
非饱和土壤紊散动力学系数测定方法及装置 Download PDFInfo
- Publication number
- CN107167576A CN107167576A CN201710322421.XA CN201710322421A CN107167576A CN 107167576 A CN107167576 A CN 107167576A CN 201710322421 A CN201710322421 A CN 201710322421A CN 107167576 A CN107167576 A CN 107167576A
- Authority
- CN
- China
- Prior art keywords
- soil
- mrow
- solute
- section
- spike
- 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
- 239000002689 soil Substances 0.000 title claims abstract description 204
- 238000003556 assay Methods 0.000 title claims abstract description 8
- 238000012544 monitoring process Methods 0.000 claims abstract description 62
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 45
- 238000000034 method Methods 0.000 claims abstract description 38
- 230000004907 flux Effects 0.000 claims abstract description 36
- 238000012360 testing method Methods 0.000 claims abstract description 21
- 239000011159 matrix material Substances 0.000 claims abstract description 20
- 230000008859 change Effects 0.000 claims abstract description 19
- 230000005012 migration Effects 0.000 claims abstract description 18
- 238000013508 migration Methods 0.000 claims abstract description 18
- 230000008878 coupling Effects 0.000 claims abstract description 3
- 238000010168 coupling process Methods 0.000 claims abstract description 3
- 238000005859 coupling reaction Methods 0.000 claims abstract description 3
- 238000002347 injection Methods 0.000 claims description 14
- 239000007924 injection Substances 0.000 claims description 14
- 230000008569 process Effects 0.000 claims description 14
- 238000002474 experimental method Methods 0.000 claims description 13
- 230000035699 permeability Effects 0.000 claims description 12
- 238000006073 displacement reaction Methods 0.000 claims description 9
- 239000000700 radioactive tracer Substances 0.000 claims description 9
- 239000004927 clay Substances 0.000 claims description 7
- 229920006395 saturated elastomer Polymers 0.000 claims description 6
- 238000004260 weight control Methods 0.000 claims description 6
- 230000001914 calming effect Effects 0.000 claims description 3
- 230000006835 compression Effects 0.000 claims description 3
- 238000007906 compression Methods 0.000 claims description 3
- 238000013178 mathematical model Methods 0.000 abstract 1
- 238000005259 measurement Methods 0.000 abstract 1
- 230000000007 visual effect Effects 0.000 abstract 1
- 239000006185 dispersion Substances 0.000 description 22
- 150000003839 salts Chemical class 0.000 description 10
- 238000009826 distribution Methods 0.000 description 7
- 238000004458 analytical method Methods 0.000 description 6
- 238000009792 diffusion process Methods 0.000 description 5
- 239000011148 porous material Substances 0.000 description 5
- 239000003337 fertilizer Substances 0.000 description 4
- 230000008595 infiltration Effects 0.000 description 4
- 238000001764 infiltration Methods 0.000 description 4
- 230000002123 temporal effect Effects 0.000 description 4
- 239000003905 agrochemical Substances 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 239000003673 groundwater Substances 0.000 description 3
- 230000007246 mechanism Effects 0.000 description 3
- 239000000203 mixture Substances 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 239000000126 substance Substances 0.000 description 3
- 238000012546 transfer Methods 0.000 description 3
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000010924 continuous production Methods 0.000 description 2
- 238000009795 derivation Methods 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 239000007788 liquid Substances 0.000 description 2
- 239000000243 solution Substances 0.000 description 2
- 229920005479 Lucite® Polymers 0.000 description 1
- 241001139441 Metopion Species 0.000 description 1
- 238000012271 agricultural production Methods 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000003344 environmental pollutant Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000003973 irrigation Methods 0.000 description 1
- 230000002262 irrigation Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 231100000719 pollutant Toxicity 0.000 description 1
- 239000004926 polymethyl methacrylate Substances 0.000 description 1
- 230000001737 promoting effect Effects 0.000 description 1
- -1 salinity Substances 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
- 239000002351 wastewater Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/24—Earth materials
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/24—Earth materials
- G01N33/245—Earth materials for agricultural purposes
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/24—Earth materials
- G01N33/246—Earth materials for water content
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Food Science & Technology (AREA)
- Analytical Chemistry (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- Medicinal Chemistry (AREA)
- Physics & Mathematics (AREA)
- Remote Sensing (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Sampling And Sample Adjustment (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明一种非饱和土壤紊散动力学系数测定方法及装置，该方法包括以下步骤：进行土样装填；按照一定间隔布设监测位点，各监测位点均布设土壤基质势监测系统和土壤示踪溶质浓度监测系统，用于对各监测位点的土壤含水率变化以及土壤中示踪溶质浓度变化进行连续监测；划分断面，测定各断面上、下边界土壤基质势确定水势差和各断面的水分运动通量；根据各断面水分运动通量，计算任一时刻t的平均断面水分通量；将土壤负压转化为土壤含水率，进一步耦合测定的断面平均水分通量，最终通过公式求得各测试位点的溶质迁移的紊散动力学系数。本发明建立数学模型，易于操作、测量精度高，实验模型容易搭建，条件易于实现，样本丰富多样，结果直观。
Description
技术领域
本发明属于土壤质量检测技术领域，具体的说是一种非饱和土壤紊散动力学系数测定方法及装置。
背景技术
随着农业生产的发展和科学技术水平的提高，要求农田合理施肥，经济灌溉和有效预测、控制土壤盐碱化，防治农药，化肥及废水灌溉引起的土壤和地下水污染，以及土壤溶质随地表径流迁移，引起水环境污染和土壤生产力降低，已经成为当今世界普遍关注的问题。对此研究涉及到的对土壤溶质的研究。对土壤溶质研究最广泛的是利用数学模型，模拟土壤溶质运动，其中，水动力弥散系数是不可或缺的参数之一。
溶解于土壤水的物质(化肥、农药、盐分、污染物等)在土壤空隙中运移，称之为土壤溶质运移。促使可溶性物质在土壤中运移的机制有：质流(也称对流)，水动力弥散。水动力弥散系数是表征一定流速下，多孔介质对某种污染物质弥散能力的参数。在宏观上反映了多孔介质中地下水流动过程和空隙结构特征对溶质运移过程的影响。包括机械弥散系数与分子扩散系数。是土壤水盐运移的一个重要参数。
非饱和土壤紊散动力学系数是化肥、农药在农田的移动规律、盐碱地水盐运动检测、地下水资源保护中不可缺少的参数。目前，虽然有关饱和条件下土壤弥散系数溶质运移各类参数的测定已较为成熟，但非饱和条件下由于土壤中实际孔隙流速随含水率的变化表现出数量级的变化，以及非饱和条件下溶质迁移非线性驱动机制，由于各类问题本身的复杂特性，使得对水动力紊散系数参数的测定相当困难，需要发展测试理论和测试技术。
国内外测定非饱和土壤弥散动力学系数的方法有：公式法、水平土柱吸渗法、瞬时剖面法。这些方法的不同在于测试的边界条件和控制过程的差异，以及与边界条件和控制过程对应的解析方法的差异，然而就测试理念来说，并没有差异，都是测定能够表述溶质迁移所表现出的宏观特性的等效参数，而不是正真意义的参数。而对等效参数，仅适合在与测试条件下相似的情况下，使用参数描述土壤中的溶质迁移，而当实际土壤条件与测试条件有较大差异的情况下，参数的适用性则显著降低。而本项研究提出的非饱和土壤水动力紊散参数测定方法则能够有效的解决这一问题。
发明内容
为了克服上述现有技术的不足，本发明提供了一种非饱和土壤紊散动力学系数测定的方法，具有操作简便、费用低、测定周期短、误差小的特点。
本发明的方法所采用的技术方案是：
一种非饱和土壤紊散动力学系数测定的方法，包括以下步骤：
步骤1：进行土样装填，在装填过程中，通过土壤初始含水率控制，土壤容重控制和土壤非饱和渗透性控制消除土壤分层和各断面不均匀现象；
步骤2：根据土样的质地，以及水流在土壤中的流动速度，按照一定间隔(5-10cm)布设监测位点，各监测位点均布设土壤基质势监测系统和土壤示踪溶质浓度监测系统，用于对各监测位点的土壤含水率变化以及土壤中示踪溶质浓度变化进行连续监测；
步骤3：基于各监测位点进行断面划分，各断面位于相邻两个监测位点中间位置，将各断面相邻的上、下监测位点所在断面分别作为各断面的上、下边界，通过各断面上、下边界的土壤基质势监测系统测定的土壤基质势确定水势差(负压坡降)，然后采用达西定律确定各断面的水分运动通量；
步骤4：在任一时刻t，根据各断面水分运动通量，计算该时刻平均断面水分通量；
步骤5：基于步骤2各监测位点测定的土壤负压，根据土壤基质势与土壤含水率之间土壤水分特征曲线关系，将土壤负压转化为土壤含水率，进一步耦合步骤4测定的断面平均水分通量，采用如下公式求得各测试位点的溶质迁移的紊散动力学系数：
其中，D为溶质迁移的紊散动力系数，J为溶质的质量通量，θ为土壤体积含水率，c为土壤示踪溶质浓度，t为时间，z为垂直坐标，q为水分运动通量，相邻监测位置节点深度差值Δz＝zi+1-zi，相邻监测时间节点差值Δt＝tj+1-tj，i表示监测位置节点，i＝1,2,3…,j表示监测时间节点，j＝1，2，3…，zi和zi+1分别为i和i+1监测位置节点的深度，1/2表示监测位置或时间节点的中点位置，tj和tj+1分别表示j和j+1时刻。
上述步骤1中，土壤初始含水率控制：土壤在105℃条件下完全烘干后，用孔径为2mm的震动筛去除影响土壤结构的大孔径土粒后，根据土壤重量，确定添加水量，采用表面喷洒拌合的方式实现土壤初始含水率，要求黏土达到10～15％的饱和含水率，要求壤土达到20～30％的饱和含水率。
上述步骤1中，土壤容重控制：根据所设定的土壤容重，确定分层装填的土壤重量，采用气动压缩机根据设定容重所确定的断面位置进行土样装填。
上述步骤1中，为控制装填土壤的质量，在106Pa的压强下进行压气实验，装填土壤控制在12小时，压缩位移和非饱和渗透性指标需要达到以下控制要求：黏土：位移小于装填土壤长度的1％，非饱和渗透性小于10-5m/s；壤土：位移小于装填土壤长度的5％，非饱和渗透性小于10-4m/s。
本发明还提供基于上述方法的一种非饱和土壤紊散动力学系数测定装置，所述装置包括边界控制系统和监测系统，边界控制系统用于控制实验边界条件以及实现土壤的示踪溶质的控制注入，监测系统包括土壤基质势监测系统和土壤示踪溶质浓度监测系统；所述边界控制系统能够实现多种条件下的多参数测试控制，包括：①定水头、定浓度示踪溶质连续注入，②定水头、脉冲示踪溶质注入，③定流量、定浓度示踪溶质连续注入，以及④定流量、脉冲示踪溶质注入；所述监测系统对示踪溶质迁移过程中的土壤含水率变化以及土壤中示踪溶质浓度变化进行全过程连续监测。
作为优选，所述土壤基质势监测系统选用负压计，所述土壤示踪溶质浓度监测系统选用盐桥传感器。
定流量边界，根据达西定律通过动态控制外压水头边界实现，土壤中的流动符合达西定律，通量与压力正比，根据土壤负压计测定的土壤基质，控制外压水头，实现定流量过程入渗。
定水头直接通过供平水装置(马氏瓶)控制外压水头即可。
根据测定的土壤含水率，以及示踪溶质(盐分)的分布，能够基于式(13)采用有限差分法确定土壤紊散动力学系数。通过对不同的工况条件下的溶质迁移动力学过程进行测试，各监测位置的土壤含水率以及示踪溶质不断发生变化，而这种变化，是由于对流和弥散作用所造成的，而紊散动力学系数的测定需要通过对土柱的所有检测点的连续过程监测，确定弥散作用的贡献，进而确定参数。
作为优选，本方法为了更加准确的得到剖面含水率，计算每一时段土柱含水量的增值ΔW：
式中，θ为计算时段第i点的含水率，xi为第i点与第i+1点的间距，t和t+1分别表示前后两个时刻。
土柱含水量的增值ΔW应等于边界处在这一时段的供水量。如两者相差悬殊，则必须检查是否计算有误，在相对误差大于5％的情况下则应修正含水率剖面，直至计算值和测定值相对误差在5％以内。
本发明方法包括测试方法，以及数据解析方法两个相互联系的部分，根据测定的土壤含水率，以及示踪溶质(盐分)的分布，基于式(13)采用有限差分法确定土壤紊散动力学系数，与现有技术相比，本发明具有以下有益效果：
(1)边界条件和控制过程对应的解析方法的差异，然而就测试理念来说，并没有差异，都是测定能够表述溶质迁移所表现出的宏观特性的等效参数，而不是正真意义的参数。等效参数所表现出的宏观特性在很大程度上依赖于实验条件，因此，当求解问题与参数的测定条件有显著差异的情况，等效参数的有效性显著降低。而本发明的方法则在程度上克服了这一问题。
(2)非饱和流动条件下，溶质迁移通量中对流通量和弥散通量的解析是非常困难的，本发明方法能够同步实现对流通量的解析，所提出的紊散水动力参数则在机理上更为准确。
附图说明
图1为非饱和弥散参数测试试验装置示意图；
图2为定流量入渗条件下土壤水分和盐分分布；
图3为定流量入渗条件下各点弥散系数和对应点的平均孔隙流速关系曲线；
图4为定水头入渗条件下土壤水分和盐分分布；
图5为定水头入渗条件下各点弥散系数和对应点的平均孔隙流速关系曲线。
具体实施方式
下面结合附图及实施例对本发明进行详细说明。以下实施例将有助于本领域技术人员进一步理解和实施本发明，但不以任何形式限制本发明。应当指出的是，对本领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干变形和改进。这些都属于本发明的保护范围。
步骤1：进行土样装填，实验装置如附图1所示，实验装填时，通过土壤初始含水率控制、容重控制和非饱和渗透性控制，消除由于土壤结构具有较强的空间变异性和各项异性而造成的分层和各断面不均匀等现象：
1.1，土壤壤初始含水率控制：土壤在105℃条件下完全烘干后，用孔径为2mm的震动筛去除影响土壤结构的大孔径土粒后，根据土壤重量，以及控制含水率(黏土：10～15％的饱和含水率，壤土20～30％的饱和含水率)，确定添加水量，采用表面喷洒拌合的方式实现土壤初始含水率。
1.2，容重控制根据设定土壤容重，确定分层装填的土壤重量，采用气动压缩机根据设定容重所确定的断面位置进行土样装填。
1.3，为控制装填土壤的质量，在106Pa的压强下进行压气实验，装填土壤在12hr时间，压缩位移和非饱和渗透性指标需要达到控制要求(黏土：位移小于装填土壤长度的1％，非饱和渗透性小于10-5m/s，壤土位移小于装填土壤长度的5％，非饱和渗透性小于10-4m/s)
实验装置(附图1)包括(1)边界控制系统，控制实验边界条件，以及实现土壤的示踪溶质的控制注入，(2)监测系统，包括土壤基质势监测系统(如负压计)和土壤中示踪溶质的监测系统(如盐桥传感器)。非饱和条件下，土壤中流速随土壤含水率变化而显著变化，为了使得测定的水动力紊散系数能够覆盖整个非饱和条件下流速变化范围，至少需要对4种边界控制条件下的2种条件进行测试。边界控制系统能够实现：①定水头、定浓度示踪溶质连续注入，②定水头、脉冲示踪溶质注入，③定流量、定浓度示踪溶质连续注入，以及④定流量、脉冲示踪溶质注入等多种条件下的多参数测试控制。监测系统对示踪溶质迁移过程中，土壤中含水率以及示踪溶质浓度变化进行全过程连续监测。
定流量边界，通过动态控制外压水头边界实现，土壤中的流动符合达西定律，通量与压力正比，根据土壤负压计测定的土壤基质，控制外压水头，实现定流量过程入渗。
定水头直接通过供平水装置(马氏瓶)控制外压水头即可。
根据测定的土壤含水率，以及示踪溶质(盐分)的分布，能够基于式(13)采用有限差分法确定土壤紊散动力学系数。
通过对不同的工况条件下的溶质迁移动力学过程进行测试，各监测点的土壤含水率以及示踪溶质不断发生变化，而这种变化，是由于对流和弥散作用所造成的，而紊散动力学系数的测定需要通过对土柱的所有检测点的连续过程监测，确定弥散作用的贡献，进而确定参数。
填装土样的容器为直径为12cm，长100cm的有机玻璃，土壤按照预定的初始含水率和浓度拌匀，按照容重进行装填，在试验土柱中按一定距离埋设负压计和盐分传感器，观测水盐运动过程。
步骤2：根据测试土样的质地以及水流在土壤中的流动速度，确定负压计和盐分传感器的安装间隔；对于黏土，按照5cm间隔要装负压计和盐分传感器，对于壤土，按照10cm间隔要装负压计和盐分传感器。在测试阶段，盐分传感器对各监测点的土壤示踪溶质浓度进行连续监测，负压计对各监测点的土壤水势进行连续监测，土壤基质势根据相应溶质条件下的土壤水势、示踪溶质浓度得到。
步骤3：基于传感器布设进行断面划分，断面位于相邻两个传感器中间位置，根据断面上、下边界位置所布设的负压计的所测定的土壤基质势确定水势差(负压坡降)后，采用达西定律确定断面的水分运动通量。
步骤4：在任一时刻t，根据各断面水分运动通量，计算该时刻断面平均水分运动通量。
步骤5：基于步骤2测定的各监测点的土壤水势，根据土壤基质势与土壤体积含水率之间土壤水分特征曲线关系，将土壤水势转化为土壤体积含水率，如附图2所示，确定所有传感器布设节点位置的土壤含水率，进一步考核土壤测定电导率关系确定土壤溶液浓度。在不同的控制边界条件下，土壤含水率以及浓度的变化范围应覆盖土壤非标和区域(附图2和附图4)。
根据所确定的土壤含水率和孔隙率，进一步耦合断面平均水分运动通量，确定断面平均流速v：
其中，A为断面面积，n为土壤孔隙率。其他符号意义同上。
采用式(13)，采用有限差分法，求得各监测点紊散动力学系数；
其中，D为溶质迁移的紊散动力系数，J为溶质的质量通量，θ为土壤体积含水率，c为土壤示踪溶质浓度，t为时间，z为垂直坐标，q为水分运动通量，相邻监测位置节点深度差值Δz＝zi+1-zi，相邻监测时间节点差值Δt＝tj+1-tj，i表示监测位置节点，i＝1,2,3…,j表示监测时间节点，j＝1，2，3…，zi和zi+1分别为i和i+1监测位置节点的深度，1/2表示监测位置或时间节点的中点位置，tj和tj+1分别表示j和j+1时刻。在此基础上，根据各个节点的土壤物理状态以及D的测定(附图3和附图5)关系，采用解析法或瞬时剖面法确立紊散动力学函数关系。
本实验采用了解析法和瞬时剖面法，分别对应实验条件为连续溶质注入和脉冲溶质注入的条件下确定紊散动力学系数。
(一)解析法
解析法是根据影响非饱和紊散动力系数的主要因素，提出的确定紊散动力学系数的直接方法。水动力弥散系数是综合反映溶质、土壤特性的参数,既和多孔介质的状况及溶质的性质有关,又受含水量和孔隙水流速度的影响。水动力弥散系数的一般表达式为:
Dsh(θ,V)＝Ds(θ)+Dh(V) 式(1)
式中：Ds(θ)为分子扩散系数；Dh(V)为机械弥散系数。
情况1：土壤介质中，当流速对应的雷诺数超过3.0后，当对流速度相当大时，机械弥散作用会大大超过分子扩散作用，以致水动力弥散中只考虑机械弥散作用，此时紊散动力学系数的表达式为：
式中αL、αT为非饱和土壤的纵向、横向弥散度，；Vi、Vj，i、j两个方向的孔隙水流速度；δij为delta函数，当两个下标相同时，取1，相异时取值为零。说明：两个方向为纵向和横向(纵向指流动方向，横向为垂直于流动方向，两个下标表示方向，相同即表示方向相同，相异即表示方向不同)
在进行数值计算，为减少数值弥散的影响，常用D′sh作为弥散系数：
D′sh＝Dsh-VΔx/2 式(3)
式中，V＝q/θ，(见式2)，Δx为渗流长度。
当土壤溶液静止时，则机械弥散完全不起作用，只剩下分子扩散作用，此时，Dsh可以用以下表达较好：
Dsh＝D0aebθ
式中，D0为溶质在自由溶液中的扩散系数，取值与溶质分子类型有关；θ为土壤体积的含水量；a，b为经验系数，a＝0.001～0.005，b＝10。
(二)瞬时剖面法
用瞬时剖面法求解紊散动力学系数，根据水分及溶质质量守恒原理,便可建立起剖面上各点的水分、溶质均衡方程,从而求出剖面上各点的水动力弥散系数。
如果已知土壤水、盐运动过程中任意两个时刻的剖面水盐分布，则根据水分以及盐分质量守恒原理，即可建立起剖面上各点的水分及盐分均衡方程，从而求出剖面各点的弥散系数。
根据质量守恒原理，一维非饱和溶质迁移方程为：
式中，J为溶质的质量通量，表示为：
其中，θ为土壤体积含水率，c为溶液浓度，t为时间，D为溶质迁移的紊散动力系数，z为垂直坐标，q为水分运动通量。由水流连续方程得：
式(4)两边积分，并在测点i位置采用有限差分法进行展开，得：
其中，i和j分别表示测点的位置和时间步长节点，zi和zi+1分别为i和i+1位置节点的深度，1/2表示位置或时间节点的重点位置，tj和tj+1分别表示j和j+1时刻。
式(7)左边可表示为：
表示zi处tj时刻的溶质的质量通量。式(7)右边采用一阶积分近似，得：
式(9)为式(7)的进一步推导，为中间步骤。
对(9)时间导数采用中心差分近似，得：
Δt＝tj+1-tj，Δz＝zi+1-zi
式(10)为式(9)的进一步推导，为中间步骤。
综合式(7)～式(10)得：
由式(5)得：
由式(11)和式(12)得：
式(13)即为求解紊散动力系数的计算公式，其中为：
式(14)为式(13)中的的表述，对应实验中所确定的各个节点的水分通量。
根据步骤1-5得到的土壤示踪溶质浓度土壤体积含水率i＝1,2,3,…m，j＝1,2,3,…n，计算：
土壤体积含水率
土壤示踪溶质浓度
土壤体积含水率
土壤示踪溶质浓度
土壤含水率和示踪剂浓度额中点值可直接采用几何平均方法计算。
因此，在已知上边界条件(zi＝0，tj，i＝1，j＝1,2,3,…n)的溶质的质量通量i＝1，j＝1,2,3,…n的条件下，
将上述数据代入式(11)中，
求得溶质的质量通量i＝1，j＝1,2,3,…n，再将相关数据代入式(11)，可求得溶质的质量通量i＝1，j＝1,2,3,…n，根据该方法可得到i＝1,2,3,…m，j＝1,2,3,…n。
根据步骤1-5得到的土壤体积含水率i＝1,2,3,…m，j＝1,2,3,…n，
因此，在已知上边界条件(zi＝0，tj，i＝1，j＝1,2,3,…n)的水分运动通量i＝1，j＝1,2,3,…n的条件下，
将上述数据代入式(14)中，
求得水分运动通量i＝1，j＝1,2,3,…n，再将数据代入式(14)，可求得水分运动通量i＝1，j＝1,2,3,…n，根据该方法可得到i＝1,2,3,…m，j＝1,2,3,…n。
当保持上边界(zi＝0，i＝1)为定流量、定浓度入渗时，
J(0,t)＝q(0,t)，c(0,t)＝q0c0 式(15)
式(15)中的J(0,t)、q(0,t)、c(0,t)、q0、c0，
J(0,t)为上边界位置(zi＝0，i＝1)、测试时间(tj，j＝1,2,3,…n)的溶质的质量通量，
q(0,t)为上边界位置(zi＝0，i＝1)、测试时间(tj，j＝1,2,3,…n)的水分运动通量，
c(0,t)为上边界位置(zi＝0，i＝1)、测试时间(tj，j＝1,2,3,…n)的土壤示踪溶质浓度，
q0和c0分别为边界的初始通量和溶度。
最简单的边界条件为q(0,t)＝0，J(0,t)＝0，即土壤水盐的再分布期间的任意两个时刻的盐分剖面分布均可对紊散动力系数进行求解。土壤水盐的再分布，也就是实验开始后，土壤中各点的水分和溶质(盐分)会不断发生变化。
通量边界和定水头边界条件下，连续注入与脉冲注入的求解方法相同。
①定水头、定浓度示踪溶质连续注入，
②定水头、脉冲示踪溶质注入，
③定流量、定浓度示踪溶质连续注入，以及
④定流量、脉冲示踪溶质注入
应当理解的是，本说明书未详细阐述的部分均属于现有技术。
应当理解的是，上述针对较佳实施例的描述较为详细，并不能因此而认为是对本发明专利保护范围的限制，本领域的普通技术人员在本发明的启示下，在不脱离本发明权利要求所保护的范围情况下，还可以做出替换或变形，均落入本发明的保护范围之内，本发明的请求保护范围应以所附权利要求为准。
Claims (6)
1.一种非饱和土壤紊散动力学系数测定方法，其特征在于，包括以下步骤：
步骤1：进行土样装填，在装填过程中，通过土壤初始含水率控制，土壤容重控制和土壤非饱和渗透性控制消除土壤分层和各断面不均匀现象；
步骤2：根据土样的质地，以及水流在土壤中的流动速度，按照一定间隔布设监测位点，各监测位点均布设土壤基质势监测系统和土壤示踪溶质浓度监测系统，用于对各监测位点的土壤含水率变化以及土壤中示踪溶质浓度变化进行连续监测；
步骤3：基于各监测位点进行断面划分，各断面位于相邻两个监测位点中间位置，将各断面相邻的上、下监测位点所在断面分别作为各断面的上、下边界，通过各断面上、下边界的土壤基质势监测系统测定的土壤基质势确定水势差，然后采用达西定律确定各断面的水分运动通量；
步骤4：在任一时刻t，根据各断面水分运动通量，计算该时刻平均断面水分通量；
步骤5：基于步骤2各监测位点测定的土壤负压，根据土壤基质势与土壤含水率之间土壤水分特征曲线关系，将土壤负压转化为土壤含水率，进一步耦合步骤4测定的断面平均水分通量，采用如下公式求得各测试位点的溶质迁移的紊散动力学系数：
<mrow>
<mi>D</mi>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<mrow>
<msubsup>
<mi>J</mi>
<mi>i</mi>
<mrow>
<mi>j</mi>
<mo>+</mo>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msubsup>
<mo>+</mo>
<mfrac>
<mrow>
<mi>&Delta;</mi>
<mi>z</mi>
</mrow>
<mrow>
<mi>&Delta;</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>&lsqb;</mo>
<msubsup>
<mrow>
<mo>(</mo>
<mi>&theta;</mi>
<mi>c</mi>
<mo>)</mo>
</mrow>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
</mrow>
<mi>j</mi>
</msubsup>
<mo>-</mo>
<msubsup>
<mrow>
<mo>(</mo>
<mi>&theta;</mi>
<mi>c</mi>
<mo>)</mo>
</mrow>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
</mrow>
<mrow>
<mi>j</mi>
<mo>+</mo>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msubsup>
<mo>&rsqb;</mo>
<mo>-</mo>
<msubsup>
<mrow>
<mo>(</mo>
<mi>q</mi>
<mi>c</mi>
<mo>)</mo>
</mrow>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>j</mi>
<mo>+</mo>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msubsup>
</mrow>
<mrow>
<msubsup>
<mi>Q</mi>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>j</mi>
<mo>+</mo>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msubsup>
<msubsup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mi>&Delta;</mi>
<mi>c</mi>
</mrow>
<mrow>
<mi>&Delta;</mi>
<mi>z</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>j</mi>
<mo>+</mo>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msubsup>
</mrow>
</mfrac>
<mo>,</mo>
</mrow>
其中，D为溶质迁移的紊散动力系数，J为溶质的质量通量，θ为土壤体积含水率，c为土壤示踪溶质浓度，t为时间，z为垂直坐标，q为水分运动通量，相邻监测位置节点深度差值Δz＝zi+1-zi，相邻监测时间节点差值Δt＝tj+1-tj，i表示监测位置节点，i＝1,2,3…,j表示监测时间节点，j＝1，2，3…，zi和zi+1分别为i和i+1监测位置节点的深度，1/2表示监测位置或时间节点的中点位置，tj和tj+1分别表示j和j+1时刻。
2.根据权利要求1所述的一种非饱和土壤紊散动力学系数测定方法，其特征在于：
上述步骤1中，土壤初始含水率控制：土壤在105℃条件下完全烘干后，用孔径为2mm的震动筛去除影响土壤结构的大孔径土粒后，根据土壤重量，确定添加水量，采用表面喷洒拌合的方式实现土壤初始含水率，要求黏土达到10～15％的饱和含水率，要求壤土达到20～30％的饱和含水率。
3.根据权利要求1所述的一种非饱和土壤紊散动力学系数测定方法，其特征在于：
上述步骤1中，土壤容重控制：根据所设定的土壤容重，确定分层装填的土壤重量，采用气动压缩机根据设定容重所确定的断面位置进行土样装填。
4.根据权利要求1所述的一种非饱和土壤紊散动力学系数测定方法，其特征在于：
上述步骤1中，为控制装填土壤的质量，在106Pa的压强下进行压气实验，装填土壤控制在12小时，压缩位移和非饱和渗透性指标需要达到以下控制要求：黏土：位移小于装填土壤长度的1％，非饱和渗透性小于10-5m/s；壤土：位移小于装填土壤长度的5％，非饱和渗透性小于10-4m/s。
5.基于权利要求1-4任一项所述方法的一种非饱和土壤紊散动力学系数测定装置，其特征在于：
所述装置包括边界控制系统和监测系统，边界控制系统用于控制实验边界条件以及实现土壤的示踪溶质的控制注入，监测系统包括土壤基质势监测系统和土壤示踪溶质浓度监测系统；所述边界控制系统能够实现多种条件下的多参数测试控制，包括：①定水头、定浓度示踪溶质连续注入，②定水头、脉冲示踪溶质注入，③定流量、定浓度示踪溶质连续注入，以及④定流量、脉冲示踪溶质注入；所述监测系统对示踪溶质迁移过程中的土壤含水率变化以及土壤中示踪溶质浓度变化进行全过程连续监测。
6.根据权利要求5所述的一种非饱和土壤紊散动力学系数测定装置，其特征在于：
所述土壤基质势监测系统选用负压计，所述土壤示踪溶质浓度监测系统选用盐桥传感器。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710322421.XA CN107167576B (zh) | 2017-05-09 | 2017-05-09 | 非饱和土壤紊散动力学系数测定方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710322421.XA CN107167576B (zh) | 2017-05-09 | 2017-05-09 | 非饱和土壤紊散动力学系数测定方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107167576A true CN107167576A (zh) | 2017-09-15 |
| CN107167576B CN107167576B (zh) | 2020-02-14 |
Family
ID=59813727
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710322421.XA Expired - Fee Related CN107167576B (zh) | 2017-05-09 | 2017-05-09 | 非饱和土壤紊散动力学系数测定方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107167576B (zh) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108759949A (zh) * | 2018-05-30 | 2018-11-06 | 河海大学 | 测量浅水区地下水出流量的装置及其操作方法 |
| CN110850064A (zh) * | 2019-12-03 | 2020-02-28 | 安徽省（水利部淮河水利委员会）水利科学研究院（安徽省水利工程质量检测中心站） | 一种检测坡面流和壤中流新老水的方法 |
| CN112161897A (zh) * | 2020-09-27 | 2021-01-01 | 武汉大学 | 大型土壤水运动实验系统负压边界控制方法 |
| CN117019863A (zh) * | 2023-08-29 | 2023-11-10 | 大庆永铸石油技术开发有限公司 | 一种基于微生物处理的盐碱地改良系统 |
| CN117091997A (zh) * | 2023-10-13 | 2023-11-21 | 四川省生态环境科学研究院 | 一种河道胶体纵向动力学分散系数直接测定方法及装置 |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1793881A (zh) * | 2005-12-20 | 2006-06-28 | 中国科学院水利部水土保持研究所 | 土壤水分、溶质运移参数自动测定装置 |
| CN202041444U (zh) * | 2011-04-14 | 2011-11-16 | 中国地质大学（武汉） | 一种测定三维水动力弥散系数的地下水污染模拟装置 |
| CN103197039A (zh) * | 2013-03-15 | 2013-07-10 | 陈秋阳 | 一种人工流场下的野外地下水弥散试验方法 |
| JP2014081233A (ja) * | 2012-10-15 | 2014-05-08 | Taisei Corp | トレーサー試験方法 |
| CN104359802A (zh) * | 2014-11-20 | 2015-02-18 | 中国地质大学（武汉） | 基于一维水平流的多段式包气带水动力弥散系数测定仪 |
| CN204594980U (zh) * | 2015-03-10 | 2015-08-26 | 湖南农业大学 | 一种测定土壤中污染物迁移通量的装置 |
| CN104865165A (zh) * | 2015-06-05 | 2015-08-26 | 武汉大学 | 全过程通量守恒的现场土壤溶质运移弥散系数测定的方法 |
-
2017
- 2017-05-09 CN CN201710322421.XA patent/CN107167576B/zh not_active Expired - Fee Related
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1793881A (zh) * | 2005-12-20 | 2006-06-28 | 中国科学院水利部水土保持研究所 | 土壤水分、溶质运移参数自动测定装置 |
| CN202041444U (zh) * | 2011-04-14 | 2011-11-16 | 中国地质大学（武汉） | 一种测定三维水动力弥散系数的地下水污染模拟装置 |
| JP2014081233A (ja) * | 2012-10-15 | 2014-05-08 | Taisei Corp | トレーサー試験方法 |
| CN103197039A (zh) * | 2013-03-15 | 2013-07-10 | 陈秋阳 | 一种人工流场下的野外地下水弥散试验方法 |
| CN104359802A (zh) * | 2014-11-20 | 2015-02-18 | 中国地质大学（武汉） | 基于一维水平流的多段式包气带水动力弥散系数测定仪 |
| CN204594980U (zh) * | 2015-03-10 | 2015-08-26 | 湖南农业大学 | 一种测定土壤中污染物迁移通量的装置 |
| CN104865165A (zh) * | 2015-06-05 | 2015-08-26 | 武汉大学 | 全过程通量守恒的现场土壤溶质运移弥散系数测定的方法 |
Non-Patent Citations (1)
| Title |
|---|
| 黄康乐: "求解非饱和纵向弥散系数的一种简便方法", 《水利学报》 * |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108759949A (zh) * | 2018-05-30 | 2018-11-06 | 河海大学 | 测量浅水区地下水出流量的装置及其操作方法 |
| CN110850064A (zh) * | 2019-12-03 | 2020-02-28 | 安徽省（水利部淮河水利委员会）水利科学研究院（安徽省水利工程质量检测中心站） | 一种检测坡面流和壤中流新老水的方法 |
| CN110850064B (zh) * | 2019-12-03 | 2021-12-03 | 安徽省（水利部淮河水利委员会）水利科学研究院（安徽省水利工程质量检测中心站） | 一种检测坡面流和壤中流新老水的方法 |
| CN112161897A (zh) * | 2020-09-27 | 2021-01-01 | 武汉大学 | 大型土壤水运动实验系统负压边界控制方法 |
| CN112161897B (zh) * | 2020-09-27 | 2021-10-22 | 武汉大学 | 大型土壤水运动实验系统负压边界控制方法 |
| CN117019863A (zh) * | 2023-08-29 | 2023-11-10 | 大庆永铸石油技术开发有限公司 | 一种基于微生物处理的盐碱地改良系统 |
| CN117019863B (zh) * | 2023-08-29 | 2024-03-15 | 大庆永铸石油技术开发有限公司 | 一种基于微生物处理的盐碱地改良系统 |
| CN117091997A (zh) * | 2023-10-13 | 2023-11-21 | 四川省生态环境科学研究院 | 一种河道胶体纵向动力学分散系数直接测定方法及装置 |
| CN117091997B (zh) * | 2023-10-13 | 2023-12-22 | 四川省生态环境科学研究院 | 一种河道胶体纵向动力学分散系数直接测定方法及装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107167576B (zh) | 2020-02-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107167576A (zh) | 非饱和土壤紊散动力学系数测定方法及装置 | |
| CN108760601B (zh) | 一种模拟冻融循环快速测试土体强度及渗透系数的试验装置 | |
| CN201615869U (zh) | 一种原位测量田间土壤饱和导水率的装置 | |
| CN108195723B (zh) | 一种加固松散砾石土的渗透注浆试验系统及方法 | |
| CN204228377U (zh) | 一种水文地质参数综合测定的多功能实验装置 | |
| CN109752306A (zh) | 动荷载扰动过程岩石渗透率测试方法及其测试系统 | |
| CN105628896B (zh) | 模拟沟床堆积物在上游泥石流扰动下的起动试验系统 | |
| Durner et al. | 75: Determining Soil Hydraulic Properties | |
| CN105986582B (zh) | 强夯加固不同地下水位地基室内模型装置及测试方法 | |
| CN108035301B (zh) | 一种渗流应力耦合下悬挂式防渗墙防渗效果模拟的试验装置和试验方法 | |
| CN105715233B (zh) | 碳酸岩缝洞单元体注采模拟实验评价装置 | |
| CN107328909B (zh) | 结构性差异非饱和土壤水动力弥散系数现场测定方法 | |
| CN105810075A (zh) | 抽水触发岩溶塌陷过程实验装置 | |
| CN101825646B (zh) | 电解质示踪测量薄层水流流速的方法 | |
| CN108593502A (zh) | 一种地下水渗流模拟装置及监测方法 | |
| CN203881646U (zh) | 一种土壤物质垂直迁移的模拟装置 | |
| CN104297124A (zh) | 一种基于电势的自发渗吸测量方法 | |
| CN104535476A (zh) | 一种野外土壤入渗实验实时监测装置 | |
| CN105544494A (zh) | 一种可控温洗-补盐土壤试验装置及其测试土样的方法 | |
| CN205538499U (zh) | 用于测定化学污染土剪切强度的全自动直剪仪 | |
| CN106644892A (zh) | 一种可加载的渗透试验装置 | |
| CN104614151A (zh) | 海岸带咸淡水突变界面的砂槽渗流模拟装置及方法 | |
| CN105862933A (zh) | 动态承压水作用的地基模型试验装置 | |
| CN106979922A (zh) | 土壤淋溶柱及土壤淋溶模拟系统 | |
| CN105866024B (zh) | 一种渗压作用下岩体超低摩擦效应的动摩擦系数测定装置 |
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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200214 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |