CN116359277A - Method for in-situ monitoring of volume weight of non-rigid soil - Google Patents
Method for in-situ monitoring of volume weight of non-rigid soil Download PDFInfo
- Publication number
- CN116359277A CN116359277A CN202310086020.4A CN202310086020A CN116359277A CN 116359277 A CN116359277 A CN 116359277A CN 202310086020 A CN202310086020 A CN 202310086020A CN 116359277 A CN116359277 A CN 116359277A
- Authority
- CN
- China
- Prior art keywords
- soil
- probe
- rigid
- heat pulse
- target non
- 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.)
- Pending
Links
- 239000002689 soil Substances 0.000 title claims abstract description 278
- 238000000034 method Methods 0.000 title claims abstract description 35
- 238000012544 monitoring process Methods 0.000 title claims abstract description 34
- 238000011065 in-situ storage Methods 0.000 title claims abstract description 23
- 239000000523 sample Substances 0.000 claims abstract description 155
- 239000004519 grease Substances 0.000 claims abstract description 12
- 229920001296 polysiloxane Polymers 0.000 claims abstract description 12
- 238000013213 extrapolation Methods 0.000 claims abstract description 5
- 238000010438 heat treatment Methods 0.000 claims description 41
- 238000001514 detection method Methods 0.000 claims description 26
- 239000004576 sand Substances 0.000 claims description 12
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 12
- 239000004927 clay Substances 0.000 claims description 9
- 239000011888 foil Substances 0.000 claims description 8
- 229920001817 Agar Polymers 0.000 claims description 6
- 230000035945 sensitivity Effects 0.000 claims description 6
- 238000001035 drying Methods 0.000 claims description 5
- 238000004364 calculation method Methods 0.000 claims description 4
- 238000012937 correction Methods 0.000 claims description 2
- 238000013461 design Methods 0.000 abstract description 17
- 238000005516 engineering process Methods 0.000 abstract description 6
- 230000008602 contraction Effects 0.000 abstract description 5
- 230000007774 longterm Effects 0.000 abstract description 4
- 238000005336 cracking Methods 0.000 abstract description 3
- 230000007547 defect Effects 0.000 abstract 1
- 238000005259 measurement Methods 0.000 description 7
- 238000002310 reflectometry Methods 0.000 description 3
- 230000005284 excitation Effects 0.000 description 2
- 230000006698 induction Effects 0.000 description 2
- 238000009434 installation Methods 0.000 description 2
- 239000002734 clay mineral Substances 0.000 description 1
- 230000018044 dehydration Effects 0.000 description 1
- 238000006297 dehydration reaction Methods 0.000 description 1
- 238000001739 density measurement Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000005429 filling process Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N25/00—Investigating or analyzing materials by the use of thermal means
- G01N25/20—Investigating or analyzing materials by the use of thermal means by investigating the development of heat, i.e. calorimetry, e.g. by measuring specific heat, by measuring thermal conductivity
-
- 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/22—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance by investigating capacitance
- G01N27/221—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance by investigating capacitance by investigating the dielectric properties
-
- 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/22—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance by investigating capacitance
- G01N27/223—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance by investigating capacitance for determining moisture content, e.g. humidity
Landscapes
- Chemical & Material Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Electrochemistry (AREA)
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
Abstract
Description
技术领域technical field
本发明涉及一种原位监测非刚性土壤容重的方法,属于土壤物理性质测定技术领域。The invention relates to a method for in-situ monitoring of non-rigid soil bulk density, and belongs to the technical field of soil physical property measurement.
背景技术Background technique
土壤容重是表征土壤结构最常用的指标之一,在调控土壤水气传输和作物生长过程中发挥着重要作用。目前,获取土壤容重的动态变化主要通过田间间隔采样的方法进行,费时费力,并且无法得到土壤结构的连续动态过程。Soil bulk density is one of the most commonly used indicators to characterize soil structure, and plays an important role in regulating soil moisture transport and crop growth. At present, the dynamic change of soil bulk density is mainly obtained through the method of sampling at intervals in the field, which is time-consuming and laborious, and the continuous dynamic process of soil structure cannot be obtained.
热脉冲-时域反射技术(Thermo-TDR)能够利用同步获取的土壤热特性及含水量,反推出土壤容重,实现土壤容重的连续监测。但是目前,热脉冲-时域反射技术还未实现商业化应用,而自制探头本身制作工艺复杂,且探头之间存在一些不确定性差异,限制了该项技术的广泛应用。有研究结合市场化的土壤水分和热脉冲探头,在不同刚性土壤类型中实现了准确的容重测定,进一步提高了该技术的应用前景。Thermal pulse-time domain reflectometry (Thermo-TDR) can use the synchronously acquired soil thermal characteristics and water content to reverse the soil bulk density and realize continuous monitoring of soil bulk density. However, at present, the thermal pulse-time domain reflectometry technology has not yet been applied commercially, and the manufacturing process of self-made probes is complex, and there are some uncertainties between the probes, which limit the wide application of this technology. Studies have combined commercial soil moisture and heat pulse probes to achieve accurate bulk density measurements in different rigid soil types, further improving the application prospects of this technology.
然而,对于富含膨胀性粘土矿物的非刚性土壤而言,土壤的剧烈收缩不仅造成探针本身明显的变形,并且失水过程中土壤裂隙逐渐发育。探针和土壤间形成的空气间隙显著增大了土壤的接触热阻,进而对容重预测产生较大误差。而在田间监测过程中,容重预测具有更多的不确定性,比如除脉冲加热过程外,野外土壤热特性的测定更易受到太阳辐射和环境温度漂移的影响,造成昼夜之间或不同天气条件下容重预测的较大偏差。因此,亟需研发准确便捷且易于推广的非刚性土壤容重原位监测方法,填补多探头联用监测土壤容重动态的技术空白。However, for the non-rigid soil rich in expansive clay minerals, the severe contraction of the soil not only caused the obvious deformation of the probe itself, but also the gradual development of soil cracks during the dehydration process. The air gap formed between the probe and the soil significantly increases the contact thermal resistance of the soil, which leads to a large error in the prediction of the bulk density. In the process of field monitoring, the prediction of bulk density has more uncertainties. For example, in addition to the pulse heating process, the determination of soil thermal properties in the field is more susceptible to the influence of solar radiation and ambient temperature drift, resulting in bulk density between day and night or under different weather conditions. large deviation from the forecast. Therefore, it is urgent to develop an accurate, convenient and easy-to-promote non-rigid soil bulk density in-situ monitoring method to fill the technical gap in monitoring soil bulk density dynamics with multiple probes.
发明内容Contents of the invention
本发明所要解决的技术问题是提供一种原位监测非刚性土壤容重的方法,针对非刚性土胀缩过程产生的接触热阻增大、田间温度漂移等问题,采用全新设计,能够高效准确测定非刚性土壤容重的动态变化。The technical problem to be solved by the present invention is to provide a method for in-situ monitoring of non-rigid soil bulk density. Aiming at problems such as increased contact thermal resistance and field temperature drift caused by non-rigid soil expansion and contraction, a new design can be used to efficiently and accurately measure Dynamic changes in bulk density of non-rigid soils.
本发明为了解决上述技术问题采用以下技术方案:本发明设计了一种原位监测非刚性土壤容重的方法,基于土壤热脉冲探针与土壤水分探针对目标非刚性土壤热导率和体积含水量的同时检测,执行如下步骤A至步骤C,实现对目标非刚性土壤的土壤容重值ρb的检测;In order to solve the above technical problems, the present invention adopts the following technical solutions: the present invention designs a method for in-situ monitoring of non-rigid soil bulk density, based on the measurement of the thermal conductivity and volume content of the target non-rigid soil by the soil heat pulse probe and the soil moisture probe. For the simultaneous detection of the amount of water, the following steps A to C are performed to realize the detection of the soil bulk density value ρ b of the target non-rigid soil;
步骤A.基于土壤热脉冲探针的探针灵敏度S、土壤热脉冲探针中加热丝长度L、土壤热脉冲探针中加热丝电阻R、土壤热脉冲探针中加热丝两端电压Uheater、以及土壤热脉冲探针对应预设探针加热时长的加热前后探针输出电压U0与U1,获得目标非刚性土壤的土壤热导率值λ,然后进入步骤B;Step A. Based on the probe sensitivity S of the soil heat pulse probe, the length L of the heating wire in the soil heat pulse probe, the resistance R of the heating wire in the soil heat pulse probe, and the voltage U heater at both ends of the heating wire in the soil heat pulse probe , and the probe output voltage U 0 and U 1 before and after heating corresponding to the preset probe heating time of the soil heat pulse probe to obtain the soil thermal conductivity value λ of the target non-rigid soil, and then enter step B;
步骤B.根据土壤水分探针所获目标非刚性土壤的土壤介电常数Eb,获得目标非刚性土壤的土壤体积含水量值θv,然后进入步骤C;Step B. According to the soil dielectric constant Eb of the target non-rigid soil obtained by the soil moisture probe, the soil volumetric water content value θ v of the target non-rigid soil is obtained, and then enter step C;
步骤C.根据目标非刚性土壤的土壤热导率值λ、土壤体积含水量值θv,针对如下模型:Step C. According to the soil thermal conductivity value λ and soil volumetric water content value θ v of the target non-rigid soil, aim at the following model:
迭代计算获得目标非刚性土壤的土壤总孔隙度φ、以及目标非刚性土壤的土壤容重值ρb,其中,α=0.67fclay+0.24,β=1.97fsand+1.87ρb-1.36fsandρb-0.95,fsand和fclay分别表示目标非刚性土壤的砂粒和粘粒含量,λdry表示干土热导率。Iterative calculation to obtain the total soil porosity φ of the target non-rigid soil and the soil bulk density value ρ b of the target non-rigid soil, where α=0.67f clay +0.24, β=1.97f sand +1.87ρ b -1.36f sand ρ b -0.95, f sand and f clay represent the sand and clay content of the target non-rigid soil, respectively, and λ dry represents the thermal conductivity of dry soil.
作为本发明的一种优选技术方案:所述步骤A中,基于土壤热脉冲探针的探针灵敏度S、土壤热脉冲探针中加热丝长度L、土壤热脉冲探针中加热丝电阻R、土壤热脉冲探针中加热丝两端电压Uheater、以及土壤热脉冲探针对应预设探针加热时长的加热前后探针输出电压U0与U1,按如下公式:As a preferred technical solution of the present invention: in the step A, based on the probe sensitivity S of the soil heat pulse probe, the heating wire length L in the soil heat pulse probe, the heating wire resistance R in the soil heat pulse probe, The voltage U heater at both ends of the heating wire in the soil heat pulse probe, and the output voltages U 0 and U 1 of the probe before and after heating corresponding to the preset probe heating time of the soil heat pulse probe are as follows:
获得目标非刚性土壤的土壤热导率值λ,其中,S表示土壤热脉冲探针的探针灵敏度,Q表示土壤热脉冲探针中每米加热丝的电功率,ΔU表示U1与U0之间的差值,然后进入步骤B。Obtain the soil thermal conductivity value λ of the target non-rigid soil, where S represents the probe sensitivity of the soil thermal pulse probe, Q represents the electric power of the heating wire per meter in the soil thermal pulse probe, and ΔU represents the difference between U 1 and U 0 The difference between, and then go to step B.
作为本发明的一种优选技术方案:所述步骤B中,根据土壤水分探针所获目标非刚性土壤的土壤介电常数Eb,按如下公式:As a preferred technical solution of the present invention: in the step B, according to the soil dielectric constant Eb of the target non-rigid soil obtained by the soil moisture probe, according to the following formula:
θv=4.3×10-6×Eb3-5.5×10-4×Eb2+2.92×10-2×Eb-5.3×10-2 θ v =4.3×10 -6 ×Eb 3 -5.5×10 -4 ×Eb 2 +2.92×10 -2 ×Eb-5.3×10 -2
获得目标非刚性土壤的土壤体积含水量值θv,然后进入步骤C。Obtain the soil volume moisture content value θ v of the target non-rigid soil, and then enter step C.
作为本发明的一种优选技术方案:所述步骤A中关于土壤热脉冲探针对目标非刚性土壤的检测,剔除预设白天时间段的检测数据,进而根据预设夜间时间段的检测数据,执行所述步骤A,获得目标非刚性土壤的土壤热导率值λ。As a preferred technical solution of the present invention: In the step A, regarding the detection of the target non-rigid soil by the soil heat pulse probe, the detection data of the preset daytime time period is eliminated, and then according to the detection data of the preset nighttime time period, Execute the step A to obtain the soil thermal conductivity value λ of the target non-rigid soil.
作为本发明的一种优选技术方案:所述步骤A中关于土壤热脉冲探针对目标非刚性土壤的检测,针对预设夜间时间段检测数据,根据如下线性外推法公式:As a preferred technical solution of the present invention: in the step A, regarding the detection of the target non-rigid soil by the soil heat pulse probe, for the preset night time period detection data, according to the following linear extrapolation formula:
U0(t)=at+bU 0 (t)=at+b
拟合获得土壤热脉冲探针对应加热前探针输出电压随时间(t)变化的趋势方程U0(t),即校正预设夜间时间段监测过程中环境背景温度对检测数据的影响,进而根据其中at,应用U1-at的结果更新U1,更新获得土壤热脉冲探针对应预设探针加热时长加热后的校正探针输出电压U1,其中,a、b表示方程拟合参数。Fitting to obtain the trend equation U 0 (t) of the output voltage of the soil heat pulse probe corresponding to the change of the probe output voltage with time (t) before heating, that is, to correct the influence of the ambient background temperature on the detection data during the monitoring process of the preset night time period, and then According to where at, update U 1 with the result of U 1 -at, update and obtain the corrected probe output voltage U 1 after the heating time of the soil heat pulse probe corresponding to the preset probe heating time, where a and b represent the equation fitting parameters .
作为本发明的一种优选技术方案:所述土壤水分探针对目标非刚性土壤检测前,分别应用土壤水分探针与烘干法获得土壤环刀样品的测定值,并获得两者之间的线性关系,用于后续校正土壤水分探针对目标非刚性土壤检测的误差。As a preferred technical solution of the present invention: before the soil moisture probe detects the target non-rigid soil, the measured value of the soil ring knife sample is obtained by applying the soil moisture probe and the drying method respectively, and the difference between the two is obtained. The linear relationship is used to subsequently correct the error of the soil moisture probe for the detection of the target non-rigid soil.
作为本发明的一种优选技术方案:所述步骤A中,基于目标非刚性土壤的土壤热导率值λ的获得,进一步根据K=Kreference×Treference/T,获得目标非刚性土壤的土壤热扩散率K,进而根据C=λ/K,获得目标非刚性土壤的土壤体积热容量C;其中,Kreference表示标准参考物琼脂凝胶20℃条件下的热扩散率,Treference表示20℃条件下标准参考物琼脂凝胶加热后输出电压降低至ΔU的37%时所需的时长,ΔU表示U1与U0之间的差值,T表示土壤热脉冲探针对应预设探针加热时长加热后、探针输出电压降低至ΔU的37%的响应时长。As a preferred technical solution of the present invention: in the step A, based on the acquisition of the soil thermal conductivity value λ of the target non-rigid soil, further according to K=K reference ×T reference /T, the soil of the target non-rigid soil is obtained The thermal diffusivity K, and then according to C=λ/K, the soil volumetric heat capacity C of the target non-rigid soil is obtained; where, K reference represents the thermal diffusivity of the standard reference agar gel at 20°C, and T reference represents the condition at 20°C The time required for the output voltage to drop to 37% of ΔU after heating the standard reference agar gel, ΔU represents the difference between U 1 and U 0 , T represents the heating time of the soil heat pulse probe corresponding to the preset probe After heating, the response time for the probe output voltage to decrease to 37% of ΔU.
作为本发明的一种优选技术方案:应用土壤热脉冲探针对目标非刚性土壤检测之前,采用高导热硅脂涂抹于土壤热脉冲探针箔片表面,并且因热脉冲探针准确测量土壤热特性的空间范围的影响(半径4mm),高导热硅脂的涂抹厚度应不超过1mm。As a preferred technical solution of the present invention: before using the soil heat pulse probe to detect the target non-rigid soil, use high thermal conductivity silicone grease to smear the surface of the soil heat pulse probe foil, and because the heat pulse probe can accurately measure the soil heat Due to the influence of the spatial range of the characteristics (radius 4mm), the thickness of the high thermal conductivity silicone grease should not exceed 1mm.
本发明所述一种原位监测非刚性土壤容重的方法,采用以上技术方案与现有技术相比,具有以下技术效果:A method for in-situ monitoring of non-rigid soil bulk density described in the present invention, compared with the prior art by adopting the above technical scheme, has the following technical effects:
(1)本发明所设计一种原位监测非刚性土壤容重的方法,针对非刚性土胀缩过程产生的接触热阻增大、田间温度漂移等问题,采用全新设计,在热脉冲-时域反射技术原理上,利用土壤热脉冲探针与土壤水分探针,填补多探头联用长期原位监测土壤容重动态的技术空白,基于土壤热导率值、土壤体积含水量值的获得,按设计模型测定非刚性土壤的土壤容重值,并针对土壤热脉冲探针,采用高导热硅脂改善土壤热脉冲探针与土壤之间的热接触,克服了原方法因非刚性土壤开裂造成容重预测误差增大的缺点;以及设计中通过线性外推法,校正了野外监测时因环境温度漂移造成的不同天气状况下的容重预测偏差;整个设计方法能够实现室内及野外条件下非刚性土壤容重的准确预测。(1) The present invention designs a method for in-situ monitoring of non-rigid soil bulk density. Aiming at problems such as increased contact thermal resistance and field temperature drift caused by non-rigid soil expansion and contraction, a new design is adopted. In the thermal pulse-time domain In principle, the reflection technology uses soil heat pulse probes and soil moisture probes to fill in the technical gaps in the long-term in-situ monitoring of soil bulk density dynamics with multiple probes. The model determines the soil bulk density value of non-rigid soil, and for the soil thermal pulse probe, uses high thermal conductivity silicone grease to improve the thermal contact between the soil thermal pulse probe and the soil, which overcomes the bulk density prediction error caused by the cracking of non-rigid soil in the original method The shortcomings of the increase; and the linear extrapolation method in the design corrected the deviation of the bulk density prediction under different weather conditions caused by the ambient temperature drift during field monitoring; the entire design method can realize the accuracy of the non-rigid soil bulk density under indoor and field conditions predict.
附图说明Description of drawings
图1a、1b、1c分别是TP01热脉冲探针和5TE水分探针连接CR3000数据采集器结构示意图;Figures 1a, 1b, and 1c are schematic diagrams of the structure of TP01 thermal pulse probe and 5TE moisture probe connected to CR3000 data collector;
图2是5TE土壤水分探针测定θv值与烘干θv之间的关系;Figure 2 is the relationship between the θ v value measured by the 5TE soil moisture probe and the drying θ v ;
图3是校正前后日间和夜间监测数据与环境温度间关系;Figure 3 is the relationship between daytime and nighttime monitoring data and ambient temperature before and after correction;
图4是容重预测值与实测值间的关系。Figure 4 is the relationship between the predicted value of bulk density and the measured value.
具体实施方式Detailed ways
下面结合说明书附图对本发明的具体实施方式作进一步详细的说明。The specific implementation manners of the present invention will be further described in detail below in conjunction with the accompanying drawings.
本发明所设计一种原位监测非刚性土壤容重的方法,实际应用当中,首先在土壤水分探针对目标非刚性土壤检测前,分别应用土壤水分探针与烘干法获得土壤环刀样品的测定值,并获得两者之间的线性关系,用于后续校正土壤水分探针对目标非刚性土壤检测的误差,具体将土壤水分探针插入利用一系列容重梯度和含水量梯度的PVC环刀样品中进行测定,环刀样品尺寸要求大于水分探针空间测定范围,且环刀材质不应对土壤介电常数测量产生影响;测定完成后PVC样品105℃烘干至恒重,通过烘干值与测定值间的线性关系,校正水分探针误差,如图2所示。The present invention designs a method for in-situ monitoring of non-rigid soil bulk density. In practical applications, firstly, before the soil moisture probe detects the target non-rigid soil, the soil moisture probe and the drying method are used to obtain the soil ring knife sample respectively. Measure the value and obtain the linear relationship between the two, which is used to correct the error of the soil moisture probe for the target non-rigid soil detection. Specifically, the soil moisture probe is inserted into the PVC ring knife using a series of bulk density gradients and water content gradients The measurement is carried out in the sample. The size of the ring knife sample is required to be larger than the spatial measurement range of the moisture probe, and the material of the ring knife should not affect the measurement of the soil dielectric constant. After the measurement is completed, the PVC sample is dried at 105°C to constant weight. The linear relationship between the measured values is corrected for the error of the moisture probe, as shown in Figure 2.
并且应用土壤热脉冲探针对目标非刚性土壤检测之前,采用高导热硅脂涂抹于土壤热脉冲探针箔片表面,并且在具体实施中,因热脉冲探针准确测量土壤热特性的空间范围的影响(半径4mm),高导热硅脂的涂抹厚度应不超过1mm;由于高导热硅脂选用导热系数14.8WmK-1、热阻<0.0028℃W-1、密度大于2.6gcm-3、粘度大于2500poise的绝缘硅脂,将其均匀地薄涂于土壤热脉冲探针箔片表面,能够降低非刚性土壤失水开裂后土壤与探针间产生的接触热阻。And before using the soil heat pulse probe to detect the target non-rigid soil, use high thermal conductivity silicone grease to smear the surface of the soil heat pulse probe foil, and in the specific implementation, because the heat pulse probe accurately measures the spatial range of the soil thermal characteristics (radius 4mm), the thickness of the high thermal conductivity silicone grease should not exceed 1mm; because the high thermal conductivity silicone grease is selected with a thermal conductivity of 14.8WmK -1 , a thermal resistance <0.0028°CW -1 , a density greater than 2.6gcm -3 , and a viscosity greater than 2500poise insulating silicone grease, evenly and thinly applied to the surface of the soil heat pulse probe foil, can reduce the contact thermal resistance between the soil and the probe after the non-rigid soil loses water and cracks.
按如下,执行探针的安装。Probe installation is performed as follows.
(1)对室内土柱实验,先填装土柱至一半高度,将TP01探针箔片放置于土面上并保证良好接触后,再填装另一半土柱至目的高度。填装过程中分层按压土柱,保证探针箔片与上下土层良好接触。(1) For the indoor soil column experiment, first fill the soil column to half the height, place the TP01 probe foil on the soil surface and ensure good contact, and then fill the other half of the soil column to the target height. During the filling process, the soil column is pressed layer by layer to ensure that the probe foil is in good contact with the upper and lower soil layers.
(2)对野外不同土层的原位监测,在安装前先用刀片在土层中划出一条间隙,将涂抹过高导热硅脂的探针箔片小心插入,用湿土将探针底座固定好,防止填土时脱落。安装填埋好后,在地表浇1-2次水使土层落实,达到土壤与探针充分接触的目的。(2) For in-situ monitoring of different soil layers in the field, use a blade to draw a gap in the soil layer before installation, carefully insert the probe foil coated with high thermal conductivity silicone grease, and use wet soil to seal the probe base Fix it well to prevent it from falling off when filling the soil. After the landfill is installed, pour water on the surface 1-2 times to make the soil layer solid, so as to achieve the purpose of full contact between the soil and the probe.
然后基于土壤热脉冲探针与土壤水分探针对目标非刚性土壤热导率和体积含水量的同时检测,执行如下步骤A至步骤C,实现对目标非刚性土壤的土壤容重值ρb的检测。Then, based on the simultaneous detection of the thermal conductivity and volumetric water content of the target non-rigid soil by the soil heat pulse probe and the soil moisture probe, the following steps A to C are performed to realize the detection of the soil bulk density value ρ b of the target non-rigid soil .
步骤A.基于土壤热脉冲探针的探针灵敏度S、土壤热脉冲探针中加热丝长度L、土壤热脉冲探针中加热丝电阻R、土壤热脉冲探针中加热丝两端电压Uheater、以及土壤热脉冲探针对应预设探针加热时长的加热前后探针输出电压U0与U1,按如下公式:Step A. Based on the probe sensitivity S of the soil heat pulse probe, the length L of the heating wire in the soil heat pulse probe, the resistance R of the heating wire in the soil heat pulse probe, and the voltage U heater at both ends of the heating wire in the soil heat pulse probe , and the probe output voltage U 0 and U 1 of the soil heat pulse probe corresponding to the preset probe heating time before and after heating, according to the following formula:
获得目标非刚性土壤的土壤热导率值λ,其中,S表示土壤热脉冲探针的探针灵敏度,Q表示土壤热脉冲探针中每米加热丝的电功率,ΔU表示U1与U0之间的差值,然后进入步骤B。Obtain the soil thermal conductivity value λ of the target non-rigid soil, where S represents the probe sensitivity of the soil thermal pulse probe, Q represents the electric power of the heating wire per meter in the soil thermal pulse probe, and ΔU represents the difference between U 1 and U 0 The difference between, and then go to step B.
实际应用当中,步骤A中关于土壤热脉冲探针对目标非刚性土壤的检测过程,进一步设计首先剔除预设白天时间段的检测数据,然后针对预设夜间时间段的检测数据,根据如下线性外推法公式:In practical applications, in step A, regarding the detection process of the target non-rigid soil by the soil heat pulse probe, the further design first eliminates the detection data of the preset daytime time period, and then for the detection data of the preset nighttime time period, according to the following linear formula Inference formula:
U0(t)=at+bU 0 (t)=at+b
拟合获得土壤热脉冲探针对应加热前探针输出电压随时间(t)变化的趋势方程U0(t),即校正预设夜间时间段监测过程中环境背景温度对检测数据的影响,如图3所示,进而根据其中at,应用U1-at的结果更新U1,更新获得土壤热脉冲探针对应预设探针加热时长加热后的校正探针输出电压U1,最后执行所述步骤A,获得目标非刚性土壤的土壤热导率值λ,其中,a、b表示方程拟合参数。Fitting to obtain the trend equation U 0 (t) of the output voltage of the soil heat pulse probe corresponding to the change of the probe output voltage with time (t) before heating, that is, to correct the influence of the ambient background temperature on the detection data during the monitoring process of the preset night time period, such as As shown in Fig. 3,
上述关于土壤热脉冲探针对目标非刚性土壤的检测校验,具体根据加热前土壤背景温度(T0)随时间(t)的变化的趋势方程:The above-mentioned detection and verification of the target non-rigid soil by the soil heat pulse probe is specifically based on the trend equation of the change of the soil background temperature (T 0 ) with time (t) before heating:
T0=at+bT 0 =at+b
扣除加热过程中热脉冲感应针温度(Tv)中的土壤背景温度的变化,得到脉冲加热单一作用下的感应针温度变化(Tc),进而拟合得到土壤热特性:Subtracting the change of soil background temperature from the temperature of the thermal pulse induction needle (T v ) during the heating process, the temperature change of the induction needle under the single action of pulse heating (T c ) is obtained, and then the thermal characteristics of the soil are obtained by fitting:
Tc=Tv-atT c =T v -at
由于本专利中TP01热脉冲探针将土壤温度变化(ΔT)转化为微小的电压输出值(ΔU),故将Tv(t)的变化趋势方程转化为:Since the TP01 heat pulse probe in this patent converts the soil temperature change (ΔT) into a small voltage output value (ΔU), the change trend equation of T v (t) is transformed into:
U0(t)=at+bU 0 (t)=at+b
将Tc(t)方程转化为U180(t)方程:Transform the T c (t) equation into the U 180 (t) equation:
U1(t)=U1-atU 1 (t) = U 1 -at
在实际应用实施中,基于目标非刚性土壤的土壤热导率值λ的获得,进一步设计根据K=Kreference×Treference/T,获得目标非刚性土壤的土壤热扩散率K,进而根据C=λ/K,获得目标非刚性土壤的土壤体积热容量C;其中,Kreference表示标准参考物琼脂凝胶20℃条件下的热扩散率,通常为0.14×10-6m2s-1,Treference表示20℃条件下标准参考物琼脂凝胶加热后输出电压降低至ΔU的37%时所需的时长,通常为19s,ΔU表示U1与U0之间的差值,T表示土壤热脉冲探针对应预设探针加热时长加热后、探针输出电压降低至ΔU的37%的响应时长。In the actual application implementation, based on the acquisition of the soil thermal conductivity value λ of the target non-rigid soil, the further design is based on K=K reference ×T reference /T to obtain the soil thermal diffusivity K of the target non-rigid soil, and then according to C= λ/K, to obtain the soil volumetric heat capacity C of the target non-rigid soil; where, K reference represents the thermal diffusivity of the standard reference agar gel at 20°C, usually 0.14×10 -6 m 2 s -1 , T reference Indicates the time required for the output voltage to drop to 37% of ΔU after the standard reference agar gel is heated at 20°C, usually 19s, ΔU represents the difference between U 1 and U 0 , T represents the soil heat pulse probe Aiming at the preset heating time of the probe, it is the response time when the output voltage of the probe drops to 37% of ΔU after heating.
步骤B.根据土壤水分探针所获目标非刚性土壤的土壤介电常数Eb,按如下公式:Step B. According to the soil dielectric constant Eb of the target non-rigid soil obtained by the soil moisture probe, according to the following formula:
θv=4.3×10-6×Eb3-5.5×10-4×Eb2+2.92×10-2×Eb-5.3×10-2 θ v =4.3×10 -6 ×Eb 3 -5.5×10 -4 ×Eb 2 +2.92×10 -2 ×Eb-5.3×10 -2
获得目标非刚性土壤的土壤体积含水量值θv,然后进入步骤C。Obtain the soil volume moisture content value θ v of the target non-rigid soil, and then enter step C.
实际应用中,这里应用最初土壤水分探针与烘干法分别所获土壤环刀样品测定值之间的线性关系,校正土壤水分探针对目标非刚性土壤检测的误差。In practical applications, the linear relationship between the initial soil moisture probe and the measured values of the soil ring knife sample obtained by the drying method is used here to correct the error of the soil moisture probe in the detection of the target non-rigid soil.
步骤C.根据目标非刚性土壤的土壤热导率值λ、土壤体积含水量值θv,针对如下模型:Step C. According to the soil thermal conductivity value λ and soil volumetric water content value θ v of the target non-rigid soil, aim at the following model:
迭代计算获得目标非刚性土壤的土壤总孔隙度φ、以及目标非刚性土壤的土壤容重值ρb,其中,α=0.67fclay+0.24,β=1.97fsand+1.87ρb-1.36fsandρb-0.95,fsand和fclay分别表示目标非刚性土壤的砂粒和粘粒含量,λdry表示干土热导率。Iterative calculation to obtain the total soil porosity φ of the target non-rigid soil and the soil bulk density value ρ b of the target non-rigid soil, where α=0.67f clay +0.24, β=1.97f sand +1.87ρ b -1.36f sand ρ b -0.95, f sand and f clay represent the sand and clay content of the target non-rigid soil, respectively, and λ dry represents the thermal conductivity of dry soil.
将上述设计应用于实际当中,土壤热脉冲探针采用Hukseflux公司TP01热脉冲探针,如图1a所示,用于长期监测土壤热特性值,该探针采用箔片式设计(宽度20mm,长度60mm,厚度0.15mm)。加热丝位于箔片中心轴上,其两侧4mm处分别有一个热电堆探针,两热电堆探针相互串联。Applying the above design to practice, the soil thermal pulse probe adopts the TP01 thermal pulse probe of Hukseflux Company, as shown in Fig. 1a, which is used for long-term monitoring of soil thermal characteristic values. 60mm, thickness 0.15mm). The heating wire is located on the central axis of the foil, and there is a thermopile probe at 4 mm on both sides, and the two thermopile probes are connected in series.
土壤水分探针用于测定土壤体积含水量,可采用常用的水分探针如TDR100、5TE、EC-5、315H等,具体选用METER公司5TE土壤水分探针,如图1b所示,该探针采用频域反射技术(FDR),根据电磁波在土壤中传播频率测得土壤的表观介电常数,进而得到土壤体积含水量(θv,cm3cm-3)。The soil moisture probe is used to measure the volumetric water content of the soil. Commonly used moisture probes such as TDR100, 5TE, EC-5, 315H, etc. can be used. Specifically, the 5TE soil moisture probe of METER company is used, as shown in Figure 1b. The probe Using frequency domain reflectometry (FDR), the apparent dielectric constant of the soil was measured according to the propagation frequency of electromagnetic waves in the soil, and then the volumetric water content of the soil (θ v , cm 3 cm -3 ) was obtained.
实际应用中关于上述步骤A至步骤C的执行中,采用数据采集器分别与土壤热脉冲探针和土壤水分探针连接,用于同步获取土壤热特性及土壤含水量值,具体应用中,使用Campbell公司CR系列数据采集器,数据采集器在室内使用时,数据采集器通过12v稳压器供电;在野外时外使用时,安装12v电池及太阳能板,这里设计具体选用Campbell公司CR3000数据采集器,如图1c所示。所述数据采集器分别与TP01土壤热脉冲探针和5TE水分探针连接,用于同步获取土壤热特性及土壤含水量值。In actual application, in the implementation of the above steps A to C, the data collector is used to connect with the soil heat pulse probe and the soil moisture probe respectively to obtain the soil thermal characteristics and soil water content values synchronously. In specific applications, use Campbell company CR series data collector, when the data collector is used indoors, the data collector is powered by a 12v voltage regulator; when it is used in the field, install a 12v battery and solar panel, the design here specifically uses the Campbell company CR3000 data collector , as shown in Figure 1c. The data collector is respectively connected with the TP01 soil thermal pulse probe and the 5TE moisture probe for synchronous acquisition of soil thermal characteristics and soil moisture content.
进一步采用LoggerNet软件进行编写控制程序,用于控制所述热脉冲探针的加热时间、测定间隔、土壤热特性值以及土壤体积含水量值的同步计算与输出。Further use LoggerNet software to write a control program, which is used to control the simultaneous calculation and output of the heating time, measurement interval, soil thermal characteristic value and soil volumetric water content value of the thermal pulse probe.
实际应用中,TP01探针加热丝电源线正极与150Ω电阻串联后与开关激发通道相连(SW12V),负极接地;热电堆信号线正负极与一对差分端口相连;加热器电压输出线正极与差分端口相连,负极接地;地线接地。每个TP01探针占据1.5对差分通道;5TE土壤水分探针电源线正负极与开关激发通道相连;数据信号线与单端相连。由于水分探针接头为3.5mm圆柱形耳机接口,为不破坏探针本身构造,外加相应尺寸的接线器连接数据采集器,如图1c所示。In practical applications, the positive pole of the TP01 probe heating wire power line is connected in series with the 150Ω resistor and then connected to the switch excitation channel (SW12V), and the negative pole is grounded; the positive and negative poles of the thermopile signal line are connected to a pair of differential ports; the positive pole of the heater voltage output line is connected to The differential ports are connected, the negative pole is grounded; the ground wire is grounded. Each TP01 probe occupies 1.5 pairs of differential channels; the positive and negative poles of the 5TE soil moisture probe power line are connected to the switch excitation channel; the data signal line is connected to the single end. Since the connector of the moisture probe is a 3.5mm cylindrical earphone interface, in order not to damage the structure of the probe itself, a corresponding size connector is added to connect the data collector, as shown in Figure 1c.
实际操作中,关于上述设计步骤A至步骤C,设定的时间间隔,每3小时启动程序。程序启动时,令加热丝加热180s,使其在周围产生稳定的径向温差,进而基于土壤热脉冲探针与土壤水分探针对目标非刚性土壤的同时检测,执行如下步骤A至步骤C,容重预测值与实测值间的关系如图4所示。In actual operation, regarding the above-mentioned design steps A to C, the set time interval starts the program every 3 hours. When the program starts, heat the heating wire for 180s to generate a stable radial temperature difference around it, and then perform the following steps A to C based on the simultaneous detection of the target non-rigid soil by the soil heat pulse probe and the soil moisture probe, The relationship between the predicted value of bulk density and the measured value is shown in Figure 4.
上述技术方案所设计原位监测非刚性土壤容重的方法,针对非刚性土胀缩过程产生的接触热阻增大、田间温度漂移等问题,采用全新设计,在热脉冲-时域反射技术原理上,利用土壤热脉冲探针与土壤水分探针,填补多探头联用长期原位监测土壤容重动态的技术空白,基于土壤热导率值、土壤体积含水量值的获得,按设计模型测定非刚性土壤的土壤容重值,并针对土壤热脉冲探针,采用高导热硅脂改善土壤热脉冲探针与土壤之间的热接触,克服了原方法因非刚性土壤开裂造成容重预测误差增大的缺点;以及设计中通过线性外推法,校正了野外监测时因环境温度漂移造成的不同天气状况下的容重预测偏差;整个设计方法能够实现室内及野外条件下非刚性土壤容重的准确预测。The method of in-situ monitoring of non-rigid soil bulk density designed by the above technical scheme adopts a new design for the problems of increased contact thermal resistance and field temperature drift caused by the expansion and contraction of non-rigid soil. Based on the principle of thermal pulse-time domain reflection technology , use the soil heat pulse probe and soil moisture probe to fill the technical gap of long-term in-situ monitoring of soil bulk density dynamics with multiple probes, based on the acquisition of soil thermal conductivity and soil volumetric water content, the non-rigid The soil bulk density value of the soil, and for the soil thermal pulse probe, high thermal conductivity silicone grease is used to improve the thermal contact between the soil thermal pulse probe and the soil, which overcomes the shortcomings of the original method of increasing the bulk density prediction error due to non-rigid soil cracking ; and the linear extrapolation method in the design corrected the bulk density prediction deviation under different weather conditions caused by ambient temperature drift during field monitoring; the entire design method can realize accurate prediction of non-rigid soil bulk density under indoor and field conditions.
上面结合附图对本发明的实施方式作了详细说明,但是本发明并不限于上述实施方式,在本领域普通技术人员所具备的知识范围内,还可以在不脱离本发明宗旨的前提下做出各种变化。The embodiments of the present invention have been described in detail above in conjunction with the accompanying drawings, but the present invention is not limited to the above embodiments, and can also be made without departing from the gist of the present invention within the scope of knowledge possessed by those of ordinary skill in the art. Variations.
Claims (8)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310086020.4A CN116359277A (en) | 2023-02-09 | 2023-02-09 | Method for in-situ monitoring of volume weight of non-rigid soil |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310086020.4A CN116359277A (en) | 2023-02-09 | 2023-02-09 | Method for in-situ monitoring of volume weight of non-rigid soil |
Publications (1)
Publication Number | Publication Date |
---|---|
CN116359277A true CN116359277A (en) | 2023-06-30 |
Family
ID=86926435
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310086020.4A Pending CN116359277A (en) | 2023-02-09 | 2023-02-09 | Method for in-situ monitoring of volume weight of non-rigid soil |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116359277A (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117191894A (en) * | 2023-11-06 | 2023-12-08 | 湖南省交通规划勘察设计院有限公司 | Soil physical property testing system and method based on thermal response and storage medium |
CN118090823A (en) * | 2024-04-26 | 2024-05-28 | 江西省水利科学院(江西省大坝安全管理中心、江西省水资源管理中心) | Method for measuring volume weight of mesoscale soil in field by using heat pulse distributed optical fiber |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113008935A (en) * | 2021-03-10 | 2021-06-22 | 中国农业大学 | In-situ measuring device and method for water content, thermal characteristics and electric conductivity of soil |
-
2023
- 2023-02-09 CN CN202310086020.4A patent/CN116359277A/en active Pending
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113008935A (en) * | 2021-03-10 | 2021-06-22 | 中国农业大学 | In-situ measuring device and method for water content, thermal characteristics and electric conductivity of soil |
Non-Patent Citations (6)
Title |
---|
YUEKAI WANG等: "Determination of soil bulk density dynamic in a Vertisol during wetting and drying cycles using combined soil water content and thermal property sensors", 《GEODERMA》, vol. 428, 15 December 2022 (2022-12-15) * |
ZHENGCHAO TIAN等: "Estimating soil bulk density with combined commercial soil water content and thermal property", 《SOIL AND TILLAGE RESEARCH》, vol. 196, 28 February 2022 (2022-02-28) * |
卢奕丽: "基于土壤热导率定位监测容重的Thermo-TDR技术", 《中国博士学位论文全文数据库 农业科技辑》, 15 August 2017 (2017-08-15), pages 11 - 12 * |
张晓: "感热平衡原理测定土壤水分蒸发的改进及分解农田蒸散的实验研究", 《中国博士学位论文全文数据库 农业科技辑》, 15 August 2014 (2014-08-15), pages 28 - 34 * |
陈照峰: "《真空绝热板300问》", 31 January 2022, 中国建材工业出版社, pages: 28 - 29 * |
高福平等: "《土石坝物理场感知与数值模拟》", 31 December 2019, 河海大学出版社, pages: 29 - 30 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117191894A (en) * | 2023-11-06 | 2023-12-08 | 湖南省交通规划勘察设计院有限公司 | Soil physical property testing system and method based on thermal response and storage medium |
CN117191894B (en) * | 2023-11-06 | 2024-02-23 | 湖南省交通规划勘察设计院有限公司 | Soil physical property testing system and method based on thermal response and storage medium |
CN118090823A (en) * | 2024-04-26 | 2024-05-28 | 江西省水利科学院(江西省大坝安全管理中心、江西省水资源管理中心) | Method for measuring volume weight of mesoscale soil in field by using heat pulse distributed optical fiber |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN116359277A (en) | Method for in-situ monitoring of volume weight of non-rigid soil | |
CN100594378C (en) | Equipment for checking-out effect of heat insulation for building material, checking-out method and evaluation method | |
CN204832095U (en) | Alternation of wetting and drying crack measuring device of indoor soil property side slope model | |
CN109668922B (en) | Monitoring device for frozen soil model test and using method thereof | |
CN108333481A (en) | A kind of detection method of insulator | |
CN106872529A (en) | The method that content of unfrozen water in frozen earth is determined by measured resistivity | |
CN110297017A (en) | Soil body electrical parameter spatial and temporal distributions indoor test system and working method under Frozen-thawed cycled | |
CN201876432U (en) | Energy environmental static sounding probe capable of detecting soil heat conduction coefficient | |
Liu et al. | Effect of soil moisture content on thermal performance of ground source heat exchangers: An electromagnetism topology-based analysis | |
CN102778473A (en) | Field detection method for thermal resistance of building envelope | |
CN108287185A (en) | A kind of sounding humidity sensor, preparation method, sounding humidity measurement system and measurement method | |
CN202204566U (en) | Novel liquid level sensor | |
CN106018492A (en) | Method, device and system for detecting insulator surface dirt and moisture | |
CN107402243A (en) | A kind of method of capacitance type sensor measurement plant water content | |
CN203672806U (en) | Heat insulation effect evaluation device for mineral powder material | |
CN103837582B (en) | A kind of micro-damage rammed earth salinity measuring nose structure for earthen ruins protection | |
CN110568262A (en) | A device and method for steady-state detection of thermal insulation coatings | |
CN113252724B (en) | Method for detecting heat preservation performance of external wall | |
CN109425551B (en) | Method for calibrating undisturbed soil of FDR frequency domain reflection type soil moisture sensor | |
CN103245833B (en) | Conductive asphalt concrete resistivity measurement method based on two-electrode method | |
JPS6324267B2 (en) | ||
CN210604785U (en) | Equipment for detecting heat-insulating coating in stable state | |
CN113970518B (en) | An in-situ high-throughput atmospheric corrosion monitoring system and monitoring method thereof | |
Zhitao et al. | Design and performance analysis of composite sensor for multilayer soil profile | |
CN203100776U (en) | Multi-parameter data sensor for soil |
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 |