CN105352870B - 一种土壤导气率瞬态测算方法 - Google Patents
一种土壤导气率瞬态测算方法 Download PDFInfo
- Publication number
- CN105352870B CN105352870B CN201510798333.8A CN201510798333A CN105352870B CN 105352870 B CN105352870 B CN 105352870B CN 201510798333 A CN201510798333 A CN 201510798333A CN 105352870 B CN105352870 B CN 105352870B
- Authority
- CN
- China
- Prior art keywords
- gas
- cylinder
- air
- soil
- pressure
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 34
- 239000002680 soil gas Substances 0.000 title claims abstract description 27
- 230000001052 transient effect Effects 0.000 title claims abstract description 19
- 239000007789 gas Substances 0.000 claims abstract description 58
- 239000002689 soil Substances 0.000 claims abstract description 47
- 239000013074 reference sample Substances 0.000 claims abstract description 14
- 230000008859 change Effects 0.000 claims description 6
- 229920001971 elastomer Polymers 0.000 claims description 6
- 230000001419 dependent effect Effects 0.000 claims description 3
- 230000014509 gene expression Effects 0.000 claims description 3
- 239000012188 paraffin wax Substances 0.000 claims description 3
- 229940099259 vaseline Drugs 0.000 claims description 3
- 239000010985 leather Substances 0.000 claims 1
- 238000005259 measurement Methods 0.000 abstract description 11
- 239000007787 solid Substances 0.000 abstract description 2
- 239000000523 sample Substances 0.000 description 6
- 230000005540 biological transmission Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 238000012625 in-situ measurement Methods 0.000 description 2
- 238000000691 measurement method Methods 0.000 description 2
- 230000035699 permeability Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 229920006395 saturated elastomer Polymers 0.000 description 1
- 238000007789 sealing Methods 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Landscapes
- Investigating Or Analyzing Materials By The Use Of Electric Means (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
本发明公开了一种土壤导气率瞬态测算方法，其通过引入已知导气率的低渗透性参照样品，在同样的供气压力下能够有效降低气体的流量，从而有效提高了野外土壤导气率的测量效率；根据土壤空气传导基本理论，依据体积守恒原理得到基于低渗透性参照样品导气率的土壤导气率计算模型，为本方法的应用奠定了坚实的理论基础。
Description
技术领域：
本发明涉及土壤物理技术领域，具体地讲是一种土壤导气率瞬态测算方法。
背景技术：
土壤时刻与外接进行能量与物质的交换，土壤内部空气也在不停的运动，并不断的与外界大气进行着交换。土壤导气率是指单位面积单位时间上土壤通过的气体数量, 是反映土壤特性对土壤空气更新速率的综合性影响指标。由于空气与水分共同存在于土壤孔隙中，土壤水分的变化必然导致土壤中空气含量的变化，因此土壤空气和水之间渗透性关系紧密，通过导气率的测定可以间接获得给定土壤饱和导水率和非饱和导水率相关信息。考虑到土壤导气率的测量比导水率的测量省时、省力且对土壤结构扰动极小，特别是在需要大规模数据测量的大田土壤水力特征研究中，基于土壤导气率的水力参数求解方法具有明显优势。
国内外学者提出不同的土壤导气率测定模型与方法。根据模型供气方式不同可分为稳态法和瞬态法两大类；根据测量过程中气体运动边界条件不同，分为一维和三维两种形式。一维测定方式通常为扰动土导气率测量，将待测土样处理后填装进土柱管测量；三维测量方式通常为原状土导气率测量，将测量仪器插入待测土样中测量。一维和三维测量方式通常又与稳态、瞬态法相结合，称一维或三维稳态法、一维或三维瞬态法。传统土壤导气率测量方法多为稳态法，该方法在研究土体的一端持续施加稳定气压，然后测量通过土体的空气数量，根据土壤空气对流方程得出土壤导气率。稳态法测量技术成熟，试验时需压缩机等相关设备提供稳定气压，并用流量计等仪器测量通过土壤的空气数量，试验设备复杂昂贵。瞬态测量方法是指通过记录被测土样密封端压力动态变化，得到土样密封端压力随时间的变化关系，根据相关模型计算得出土壤导气率。瞬态模型的优势是无需测量通过土样的空气流量，试验所需供气压力小，对土样扰动小，但计算稍显复杂。
土壤导气率原位测量中会受到一些实际的问题的困扰，比如利用于瞬态法时，当测量透气性良好土壤时，每次测量用气量很大，须用打气筒或者其他供气源对储气筒连续的供气，非常耗电或者消耗体力，单位时间内测量的样点数目有限，导致测量效率大打折扣。如何提高测量效率是土壤导气率原位测量中亟待解决的关键技术难题。
发明内容：
本发明的目的是克服已有技术的不足，而提供一种土壤导气率瞬态测算方法；主要解决现有的土壤导气率瞬态测算法测量效率低的问题。
本发明的技术方案是：一种土壤导气率瞬态测算方法，采用如下装置，装置包括供气源，供气源通过导气软管连通储气筒上的第一导气管，储气筒上的另一支第二导气管通过导气软管连通测筒上的第三导气管，测筒上另一支第四导气管通过软管连通气压测量计；所述的第一导气管和第二导气管穿过密封在储气筒上的第一橡皮塞与储气筒相通，第三导气管和第四导气管穿过密封在测筒上的第二橡皮塞与测筒相通；其特殊之处在于，测算方法包括以下步骤：
1）将风干的被测土壤按照预定的容重均匀填装于内径为D（m）的测筒，填装的深度记为L2(m)；
2）将直径为D(m)、高度为L1(m)、导气率为Ka1(m2)的低渗透性参照样品放置于测筒内被测土壤之上，在测筒内壁涂抹石蜡或者凡士林，防止测筒与低渗透性参照样品之间漏气；
3）利用供气源向测筒供气，当气压测量计测定的测筒内空气相对压强（测筒内空气压强减去外界大气压）达到某一预定值△P(pa)时停止供气；
4）停止供气同时开始计时，并观测气压测量计所测测筒压强的变化，按照先疏后密的原则记录测筒内空气压强P(t)变化过程对应的累积时间t，当测筒内空气的压强接近大气压时，停止观测；
5）以时间t为因变量，以测筒内空气的相对压强（测筒内空气压强减去外界大气压即P(t)-Patm）的自然对数值ln（P(t)-Patm)为自变量，确定二者的数量关系，并以二者确定的直线斜率作为被测土壤的空气传导特征参数s；
6）将得到的空气传导特征参数s代入公式 即可计算土壤的导气率Ka2(m2)；
公式中V是储气筒和测筒内空间所占的容积m3，A表示直径为D的测筒的横截面积m2；为空气粘滞系数，与所处环境温度T有关，代表周围的大气压力；D、L1、L2、Ka1的含义同步骤（1）、（2）、（3）。
进一步的，该方法计算土壤的导气率采用的土壤导气率计算模型为：
利用供气源向测筒供气，当气压测量计测定的测筒内空气相对压强（测筒内空气压强减去外界大气压）达到某一预定值△P(pa)时，停止供气；停止供气同时开始计时，并观测气压测量计所测测筒内空气压强的变化，按照先疏后密的原则记录测筒内空气压强P(t)变化过程对应的累积时间t，当测筒内空气压强接近外界大气压时，停止观测；以时间t为因变量，以测筒内空气的相对压强（测筒内空气压强减去外界大气压即P(t)-Patm）的自然对数值Ln（P(t)-Patm)为自变量，确定二者的数量关系，并以二者确定的直线斜率作为被测土壤的空气传导特征参数s；已知低渗透性参照样品的导气率为Ka1，设被测土壤的导气率为Ka2，设低渗透性参照样品所对应的供气空间体积为V1，被测土壤对应的供气空间体积为V2，根据瞬态导气率传输理论可得：
式（1）（2）中的各有关参数的含义同前；
由式（1）可得：
考虑到
则:
代式(5)入式(2)得到：。
本发明所述的一种土壤导气率瞬态测算方法与已有技术相比具有突出的实质性特点和显著进步：1、通过引入已知导气率的低渗透性参照样品，在同样的供气压力下能够有效降低气体的流量，从而大为减少在土壤导气率测量过程中的用气量，达到同样的用气量条件下显著增加测量样本数目的目的，从而有效提高了野外土壤导气率的测量效率；2、根据土壤空气传到基本理论，在对供应气体体积分解的基础上，对被测土壤和低渗透性参照样品分别利用一维瞬态和三维瞬态导气理论，依据体积守恒原理得到基于参照样品导气率的土壤导气率计算模型，为本方法的应用奠定了坚实的理论基础。
附图说明：
图1是本发明的测量装置连接示意图。
具体实施方式：
下面结合附图和实施例，对本发明的具体实施方式作进一步详细描述；以下实施例用于说明本发明，但不用来限制本发明的范围。
实施例1，参见图1，供气源1采用打气筒，打气筒的气门嘴2通过导气软管连通储气筒3上的第一导气管5，储气筒3上的另一支第二导气管4通过导气软管连通测筒7上的第三导气管11，测筒7上另一支第四导气管12通过软管连通气压测量计8；第一导气管5和第二导气管4穿过密封在储气筒3上的第一橡皮塞6与储气筒3相通，第三导气管11和第四导气管12穿过密封在测筒7上的第二橡皮塞13与测筒7相通，形成测量装置。
采用上述测量装置的土壤导气率计算模型为：
利用打气筒向测筒供气，当气压测量计测定的测筒内空气相对压强达到某一预定值△P(pa)时，停止供气；停止供气同时开始计时，并观测气压测量计所测测筒内空气压强的变化，按照先疏后密的原则记录测筒内空气压强P(t)变化过程对应的累积时间t，当P(t)接近外界大气压时，停止观测；以时间t为因变量，以测筒内空气的相对压强（测筒内空气压强减去大气压即P(t)-Patm）的自然对数值Ln（P(t)-Patm)为自变量，确定二者的数量关系，并以二者确定的直线斜率作为被测土壤的空气传导特征参数s；已知低渗透性参照样品的导气率为Ka1，设被测土壤的导气率为Ka2，设低渗透性参照样品所对应的供气空间体积为V1，被测土壤对应的供气空间体积为V2，根据瞬态导气率传输理论可得：
式（1）（2）中的各有关参数的含义同前；
由式（1）可得：
考虑到
则:
代式(5)入式(2)得到：。
利用上述装置及土壤导气率计算模型进行测算，测算方法如下：
将风干的被测土壤9按照预定的容重均匀填装于内径为D（m）的测筒，填装的深度记为L2(m)；将直径为D(m)、高度为L1(m)、导气率为Ka1(m2)的低渗透性参照样品10放置于测筒内被测土壤9之上，在测筒内壁涂抹石蜡或者凡士林，防止测筒与低渗透性参照样品之间漏气，按照图1所示连接测量装置；利用打气筒（或者其他供气源）向测筒供气，当气压测量计测定的测筒内空气相对压强达到某一预定值△P(pa)时停止供气；停止供气同时开始计时，并观测气压测量计所测测筒内空气压强的变化，按照先疏后密的原则记录测筒内空气压强P(t)变化过程对应的累积时间t，当P(t)接近外界大气压时，停止观测；以时间t为因变量，以测筒内空气的相对压强的自然对数值ln（P(t)-Patm)为自变量,确定二者的数量关系，并以二者确定的直线斜率作为被测土壤的空气传导特征参数s；将得到的空气传导特征参数s代入公式即可计算土壤的导气率Ka2(m2)，公式中V是储气筒和测筒内空间所占的容积m3，A表示直径为D的低渗透性参照样品的横截面积m2；为空气粘滞系数，与所处环境温度T有关，代表周围的大气压力，。
Claims (1)
1.一种土壤导气率瞬态测算方法，采用如下装置，装置包括供气源（1），供气源（1）通过导气软管连通储气筒（3）上的第一导气管（5），储气筒（3）上的另一支第二导气管（4）通过导气软管连通测筒（7）上的第三导气管（11），测筒（7）上另一支第四导气管（12）通过软管连通气压测量计（8）；所述的第一导气管（5）和第二导气管（4）穿过密封在储气筒（3）上的第一橡皮塞（6）与储气筒（3）相通，第三导气管（11）和第四导气管（12）穿过密封在测筒（7）上的第二橡皮塞（13）与测筒（7）相通；其特征在于，测算方法包括以下步骤：
1）将风干的被测土壤（9）按照预定的容重均匀填装于内径为D（m）的测筒（7），填装的深度记为L2(m)；
2）将直径为D(m)、高度为L1(m)、导气率为Ka1(m2)的低渗透性参照样品（10）放置于测筒（7）内被测土壤（9）之上，在测筒（7）内壁涂抹石蜡或者凡士林，防止测筒（7）与低渗透性参照样品（10）之间漏气；
3）利用供气源（1）向测筒（7）供气，当气压测量计（8）测定的测筒（7）内空气相对压强即测筒内空气压强减去外界大气压达到某一预定值△P(Pa)时停止供气；
4）停止供气同时开始计时，并观测气压测量计（8）所测测筒（7）压强的变化，按照先疏后密的原则记录测筒内空气压强P(t)变化过程对应的累积时间t，当测筒（7）内空气的压强接近大气压时，停止观测；
5）以时间t为因变量，以测筒（7）内空气的相对压强即测筒内空气压强减去外界大气压即P(t)-Patm的自然对数值ln（P(t)-Patm)为自变量，确定二者的数量关系，并以二者确定的直线斜率作为被测土壤（9）的空气传导特征参数s；
6）将得到的空气传导特征参数s代入公式即可计算土壤的导气率Ka2(m2)；
公式中V是储气筒和测筒内空间所占的容积（m3 ），A表示直径为D的测筒的横截面积（m2）；为空气粘滞系数，与所处环境温度T有关，（），代表周围的大气压（）；D、L1、L2、Ka1的含义同步骤1）、2）、3）。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510798333.8A CN105352870B (zh) | 2015-11-19 | 2015-11-19 | 一种土壤导气率瞬态测算方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510798333.8A CN105352870B (zh) | 2015-11-19 | 2015-11-19 | 一种土壤导气率瞬态测算方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105352870A CN105352870A (zh) | 2016-02-24 |
| CN105352870B true CN105352870B (zh) | 2018-07-03 |
Family
ID=55328876
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510798333.8A Expired - Fee Related CN105352870B (zh) | 2015-11-19 | 2015-11-19 | 一种土壤导气率瞬态测算方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105352870B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109709020B (zh) * | 2019-01-24 | 2023-11-07 | 中南大学 | 非饱和冻土的气透性测试系统及测试方法 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP0577300A3 (zh) * | 1992-06-24 | 1994-04-06 | Halliburton Co | |
| CN102288528A (zh) * | 2011-08-05 | 2011-12-21 | 张振华 | 一种土壤饱和导水率和导气率原位测量系统 |
| CN102323180A (zh) * | 2011-08-12 | 2012-01-18 | 张振华 | 一种瞬态三维土壤导气率测算方法及测试装置 |
| CN102680352A (zh) * | 2012-05-10 | 2012-09-19 | 张振华 | 一种冻土导气率瞬态测量系统及测量方法 |
| CN102680351A (zh) * | 2012-05-10 | 2012-09-19 | 张振华 | 一种冻土导气率稳态测量系统及测量方法 |
| CN102692359A (zh) * | 2012-05-10 | 2012-09-26 | 张振华 | 一种土壤水气传导率原位测量系统与测量方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6843106B2 (en) * | 2000-03-08 | 2005-01-18 | Rensselaer Polytechnic Institute | Differential permeometer |
-
2015
- 2015-11-19 CN CN201510798333.8A patent/CN105352870B/zh not_active Expired - Fee Related
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP0577300A3 (zh) * | 1992-06-24 | 1994-04-06 | Halliburton Co | |
| CN102288528A (zh) * | 2011-08-05 | 2011-12-21 | 张振华 | 一种土壤饱和导水率和导气率原位测量系统 |
| CN102323180A (zh) * | 2011-08-12 | 2012-01-18 | 张振华 | 一种瞬态三维土壤导气率测算方法及测试装置 |
| CN102680352A (zh) * | 2012-05-10 | 2012-09-19 | 张振华 | 一种冻土导气率瞬态测量系统及测量方法 |
| CN102680351A (zh) * | 2012-05-10 | 2012-09-19 | 张振华 | 一种冻土导气率稳态测量系统及测量方法 |
| CN102692359A (zh) * | 2012-05-10 | 2012-09-26 | 张振华 | 一种土壤水气传导率原位测量系统与测量方法 |
Non-Patent Citations (1)
| Title |
|---|
| 土壤导气率测量方法与计算模型研究;李陆生;《中国优秀硕士学位论文全文数据库信息科技辑》;20120915(第9期);第2.2.2节、第2.3.1节、第2.3.2.2节、第4.1.1.1节、图2.3、图4.2 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105352870A (zh) | 2016-02-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106970000B (zh) | 煤/页岩超高压气体吸附和渗流实验评价页岩气吸附方法 | |
| CN102323180B (zh) | 一种瞬态三维土壤导气率测算方法及测试装置 | |
| CN103674799B (zh) | 一种测定气体在多孔隙介质中轴向扩散系数的装置及方法 | |
| CN104359819A (zh) | 低渗致密岩心气水相对渗透率测定装置及测定方法 | |
| CN201747363U (zh) | 煤层气完井方式评价实验装置 | |
| CN104297126B (zh) | 低渗透储层气体渗流启动压力梯度测量装置及测量方法 | |
| CN208140705U (zh) | 一种可自动采集孔压的渗透固结仪 | |
| CN105806712B (zh) | 孔内水压力测试装置及测试系统 | |
| CN103927913B (zh) | 一种深部地层环境二氧化碳地质储存模拟实验系统 | |
| CN108072749A (zh) | 一种隧道盾构开挖高精度渗流示踪模拟试验装置 | |
| CN105954499B (zh) | 湿陷黄土场地进行劈裂注浆加固后湿陷场地评测方法及装置 | |
| WO2020048071A1 (zh) | 一种变水头渗透系数测量系统和测量方法 | |
| CN105352871B (zh) | 一种土壤导气率原位高效测量方法 | |
| CN104389592A (zh) | 底水带油环凝析气藏水淹层油损失评价实验测试方法 | |
| CN105352870B (zh) | 一种土壤导气率瞬态测算方法 | |
| CN102680352A (zh) | 一种冻土导气率瞬态测量系统及测量方法 | |
| CN102692359A (zh) | 一种土壤水气传导率原位测量系统与测量方法 | |
| CN105424543B (zh) | 一种土壤导气率稳态测算方法 | |
| CN102323017B (zh) | 一种基于瞬态通气原理的地下滴灌堵塞程度原位定量诊断方法及诊断系统 | |
| CN204944999U (zh) | 一种简易增湿土样孔隙水压力实时测试装置 | |
| CN207964563U (zh) | 尾波特征动态监测致密岩心流体侵入过程的测量装置 | |
| CN105352845B (zh) | 一种基于稳态原理的土壤导气率原位测量方法 | |
| CN107314952B (zh) | 一种测量极低毛细数下动态接触角的方法及系统 | |
| CN208076324U (zh) | 页岩现场取芯过程中损失气量模拟系统 | |
| CN205691428U (zh) | 孔内水压力测试装置及测试系统 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| CB03 | Change of inventor or designer information | ||
| CB03 | Change of inventor or designer information |
Inventor after: Zhang Zhenhua Inventor after: Zhu Jinjin Inventor after: Yu Junbao Inventor after: Yang Runya Inventor after: Sun Junna Inventor before: Zhang Zhenhua Inventor before: Yang Runya |
|
| TA01 | Transfer of patent application right |
Effective date of registration: 20171129 Address after: 264025 Hongqi Road, Zhifu District, Shandong, China, No. 186, No. Applicant after: Ludong University Address before: 264025 School of geography and planning, Ludong University, 186 Hongqi Road, Zhifu District, Shandong, Yantai Applicant before: Zhang Zhenhua |
|
| TA01 | Transfer of patent application right | ||
| CB03 | Change of inventor or designer information | ||
| CB03 | Change of inventor or designer information |
Inventor after: Li Xueyan Inventor after: Zhu Jinjin Inventor after: Zhang Zhenhua Inventor after: Yu Junbao Inventor after: Sun Junna Inventor before: Zhang Zhenhua Inventor before: Zhu Jinjin Inventor before: Yu Junbao Inventor before: Yang Runya Inventor before: Sun Junna |
|
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180703 Termination date: 20201119 |