CN109164033A - A method of the acquisition unsaturated soil waterpower conductance based on resistivity method - Google Patents
A method of the acquisition unsaturated soil waterpower conductance based on resistivity method Download PDFInfo
- Publication number
- CN109164033A CN109164033A CN201811295465.9A CN201811295465A CN109164033A CN 109164033 A CN109164033 A CN 109164033A CN 201811295465 A CN201811295465 A CN 201811295465A CN 109164033 A CN109164033 A CN 109164033A
- Authority
- CN
- China
- Prior art keywords
- soil
- resistivity
- acquisition
- waterpower
- conductance
- 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
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
- G01N15/088—Investigating volume, surface area, size or distribution of pores; Porosimetry
-
- 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/041—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance by investigating resistance of a solid body
Abstract
The present invention relates to hydrogeologys, field of geophysical exploration, specifically disclose a kind of method of acquisition unsaturated soil waterpower conductance based on resistivity method, the following steps are included: the topography and geomorphology feature of (1) according to area to be evaluated, apparent resistivity acquisition and pedotheque collection scheme are formulated, area's apparent resistivity and corresponding pedotheque are studied in acquisition.(2) Inversion Calculation is carried out to collected apparent resistivity data, obtains the soil true resistance rate distribution in research area.(3) according to Archie equation, the quantitative relationship between soil resistivity and soil moisture is established.(4) porosity of pedotheque is determined by laboratory experiment, to obtain the distribution of research area's soil porosity.The present invention can greatly reduce evaluated error, parameter explicit physical meaning, measurement is simple, convenience of calculation, provides more accurate supplemental characteristic for hydrogeology, geophysical exploration, feasibility is high to provide convenience in the research of large area.
Description
Technical field
The present invention relates to hydrogeology, field of geophysical exploration, specifically a kind of acquisition based on resistivity method is non-full
With the method for soil waterpower conductance.
Background technique
The selection determination of soil hydraulic parameters is basis and the computer mould of soil moisture movement and solute migration prediction
The key that type calculates, and influence the main factor of prediction order of accuarcy.Soil unsaturation hydraulic conductivity is soil waterpower
One important parameter of characteristic, has closely with the continuity of soil particle partial size, the distribution of hole, pore shape and hole
Connection.Since these parameters have strong Spatial Variability, so that Direct Determination is past in the research of large area
It is past relatively difficult, and influenced by soil texture space structure variability, there are biggish errors for measurement.To overcome these to lack
Point, numerous researchers attempt to set up certain functional relation between soil texture and hydraulic parameters, to find a kind of letter
The method that soil hydraulic parameters easily are effectively predicted.But therefore the disadvantages of such methods are indefinite there are parameter physical significance utilizes
New geophysical exploration means, establishing a kind of simple unsaturated soil waterpower conductance acquisition methods still has stronger reason
Refer to realistic meaning.
Summary of the invention
The method for obtaining unsaturated soil waterpower conductance based on resistivity method that the purpose of the present invention is to provide a kind of,
To solve the problems mentioned in the above background technology.
To achieve the above object, the invention provides the following technical scheme:
A method of the acquisition unsaturated soil waterpower conductance based on resistivity method, comprising the following steps:
(1) according to the topography and geomorphology feature in area to be evaluated, apparent resistivity acquisition and pedotheque collection scheme, acquisition are formulated
Study area's apparent resistivity and corresponding pedotheque.
(2) Inversion Calculation is carried out to collected apparent resistivity data, obtains the soil true resistance rate distribution in research area.
(3) it according to Archie equation, by way of parameter calibration, determines Archie equation relevant parameter, establishes electric resistance of soil
Quantitative relationship between rate and soil moisture.
(4) porosity of pedotheque is determined by laboratory experiment, to obtain the distribution of research area's soil porosity.
(5) according to fractal theory, the quantitative model of unsaturated conductance and soil moisture, porosity is established.
Further: four parameters that Archie equation is related to use the mode of parameter calibration to obtain, and are built using matlab
Vertical four parameter fit-procedures.
It is further: to obtain the apparent resistivity distribution in research area using high-density resistivity instrument in step (1).
Compared with prior art, the beneficial effects of the present invention are: the present invention can for large area research in provide just
Benefit greatly reduces evaluated error, and parameter explicit physical meaning, measurement is simple, convenience of calculation, for hydrogeology, earth object
Reason exploration provides more accurate supplemental characteristic, and feasibility is high.
Detailed description of the invention
Fig. 1 is the flow chart of a kind of method of acquisition unsaturated soil waterpower conductance based on resistivity method.
Specific embodiment
Referring to Fig. 1, in the embodiment of the present invention, a kind of acquisition unsaturated soil waterpower conductance based on resistivity method
Method, the specific steps are as follows:
Step 1: according to the topography and geomorphology feature in area to be evaluated, establishing reasonable dimensional high-density resistance for region to be measured
Rate measures cloth pole mode, is distributed using the apparent resistivity that high-density resistivity instrument obtains research area, and according to the measurement of resistivity
The pedotheque of principle acquisition corresponding position.
Step 2: being established according to global newton-Gauss equation and minimize disturbing for the available resistivity inversion of target equation
Dynamic equation (1):
S is sensitivity matrix in formula (1), and d is observation data, and ρ is the resistivity of model, and G is forward operator, RddWith Rmm
Respectively indicate the covariance matrix of observation data and analogue data, ρ0For initial resistivity model, Δ ρkAfter this iteration of kth
Artifical resistance rate.It is distributed according to soil true resistance rate can be obtained after the iterative calculation of this model.
Step 3: according to Archie equation (2) (3), establishing the quantitative relationship between resistivity and soil moisture.
Formula (2): ρ=ρ0S-β
Formula (3):
ρ is unsaturated soil resistivity, ρ in formula (2) (3)0Soil resistivity when saturation, θ are soil moisture content, θrWith
θsRespectively saturated soil moisture content and soil residue moisture content, β is empirical parameter, related with the physical property of soil.
To avoid measuring bring error, four parameters that the Archie equation in the present invention is related to use parameter rate
Fixed mode obtains, and establishes four parameter fit-procedures using matlab, that is, can determine the parameter of Archie equation, so that it is determined that resistance
Quantitative relationship between rate and soil moisture.
Step 4: by laboratory experiment, determining the porosity ε of pedotheque.
Step 5: according to soil fractal theory formula (4) (5), determining the hydraulic conductivity power of unsaturated soil;
Formula (4)
Formula (5)
K is unsaturation waterpower conductance in formula (4) (5), and θ is soil moisture content, and D is fractal parameter, and ε is soil aperture
Degree.
The porosity ε of the soil moisture content θ and soil that obtain in step 3,4, which are brought into formula (4) (5), can be obtained
Unsaturated soil waterpower conductance.
Unsaturated soil waterpower conductance based on resistivity method obtains case verification:
Choosing Linan City of Zhejiang Province prosperousization Catchment, nearby one piece of test site carries out, and about 9 square metres of test site area,
It is measuring point number 100, each at a distance of 0.3m, apparent resistivity measurement is carried out using proud roc E60DN, measurement uses secondary device, can lead to
The distance of change measuring electrode is crossed to obtain the apparent resistivity value of different depth.Meanwhile corresponding depth is acquired according to test
Pedotheque.Inversion Calculation is iterated to two-dimentional apparent resistivity data is measured according to formula (1).According to resistivity inversion result
Available each measuring point soil true resistance rate value ρ.According to the soil moisture data of actual measurement and soil resistivity data, to Ah
Odd formula (2) (3) carries out parameter fitting, determines parameter ρ0、θr、θs, β value, establish soil resistivity and soil moisture content
Transforming relationship obtains soil moisture content distribution.Then at the porosity ε of indoor measurement soil, soil moisture content and hole will be obtained
Degree brings formula (4) (5) into, and the unsaturated soil moisture conductivity for obtaining research area is 2.321m/d.Pass through the bicyclic test of tradition
The soil moisture conductivity for measuring and being calculated research area is 2.305m/d, and the two differs 0.6%, as a result reliably, invented party
Case feasibility is high.
In addition, it should be understood that although this specification is described in terms of embodiments, but not each embodiment is only wrapped
Containing an independent technical solution, this description of the specification is merely for the sake of clarity, and those skilled in the art should
It considers the specification as a whole, the technical solutions in the various embodiments may also be suitably combined, forms those skilled in the art
The other embodiments being understood that.
Claims (6)
1. a kind of method of the acquisition unsaturated soil waterpower conductance based on resistivity method, which is characterized in that including following step
It is rapid:
(1) according to the topography and geomorphology feature in area to be evaluated, apparent resistivity acquisition and pedotheque collection scheme, acquisition research are formulated
Area's apparent resistivity and corresponding pedotheque;
(2) Inversion Calculation is carried out to collected apparent resistivity data, obtains the soil true resistance rate distribution in research area;
(3) according to Archie equation, by way of parameter calibration, determine Archie equation relevant parameter, establish soil resistivity with
Quantitative relationship between soil moisture;
(4) porosity of pedotheque is determined by laboratory experiment, to obtain the distribution of research area's soil porosity;
(5) according to fractal theory, the quantitative model of unsaturated conductance and soil moisture, porosity is established.
2. a kind of method of acquisition unsaturated soil waterpower conductance based on resistivity method according to claim 1,
It is characterized in that: obtaining the apparent resistivity distribution in research area using high-density resistivity instrument in step (1).
3. a kind of method of acquisition unsaturated soil waterpower conductance based on resistivity method according to claim 2,
It is characterized in that: the formula of resistivity data Inversion Calculation in step (2) are as follows:
Wherein S is sensitivity matrix, and d is observation data, and ρ is the resistivity of model, and G is forward operator, RddWith RmmRespectively indicate sight
The covariance matrix of measured data and analogue data, ρ0For initial resistivity model, Δ ρkFor the artifical resistance rate after this iteration of kth.
4. a kind of method of acquisition unsaturated soil waterpower conductance based on resistivity method according to claim 3,
It is characterized in that: establishing the quantitative relationship formula between resistivity and soil moisture in step (3) are as follows:
ρ=ρ0S-β,
ρ is unsaturated soil resistivity, ρ in formula0Soil resistivity when saturation, θ are soil moisture content, θrAnd θsIt is respectively full
With soil moisture content and soil residue moisture content, β is empirical parameter, related with the physical property of soil.
5. a kind of method of acquisition unsaturated soil waterpower conductance based on resistivity method according to claim 4,
Be characterized in that: four parameters that Archie equation is related to use the mode of parameter calibration to obtain, and establish four parameters using matlab
Fit procedure.
6. a kind of method of acquisition unsaturated soil waterpower conductance based on resistivity method according to claim 5,
It is characterized in that: determining the formula of the hydraulic conductivity power of unsaturated soil in step (5) are as follows:
K is unsaturation waterpower conductance in formula, and θ is soil moisture content, and D is fractal parameter, and ε is soil porosity;By step (3)
The soil moisture content θ of the middle acquisition and porosity ε of soil, which is brought into formula, can be obtained unsaturated soil waterpower conductance.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811295465.9A CN109164033B (en) | 2018-11-01 | 2018-11-01 | Method for acquiring hydraulic conductivity of unsaturated soil based on resistivity method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811295465.9A CN109164033B (en) | 2018-11-01 | 2018-11-01 | Method for acquiring hydraulic conductivity of unsaturated soil based on resistivity method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109164033A true CN109164033A (en) | 2019-01-08 |
CN109164033B CN109164033B (en) | 2020-12-22 |
Family
ID=64876595
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811295465.9A Expired - Fee Related CN109164033B (en) | 2018-11-01 | 2018-11-01 | Method for acquiring hydraulic conductivity of unsaturated soil based on resistivity method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109164033B (en) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109884131A (en) * | 2019-02-23 | 2019-06-14 | 机械工业勘察设计研究院有限公司 | One kind soil layer in submerging test of testing pits soaks lossless continuous probe method |
CN110836799A (en) * | 2019-11-15 | 2020-02-25 | 辽宁工程技术大学 | Suction control type rapid unsaturated soil sample preparation control system and method |
CN110988045A (en) * | 2019-12-24 | 2020-04-10 | 西安理工大学 | Soil body conductivity measuring method based on active heating optical fiber method |
CN111208042A (en) * | 2020-02-07 | 2020-05-29 | 中国科学院武汉岩土力学研究所 | Device and method for inverting hydraulic parameters of unsaturated waste soil |
CN112162014A (en) * | 2020-09-02 | 2021-01-01 | 塔里木大学 | Cotton field soil profile water data processing method based on electromagnetic induction data |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5861750A (en) * | 1995-01-09 | 1999-01-19 | Anderson; Dennis M. | Geophysical methods and apparatus for determining the hydraulic conductivity of porous materials |
CN101105489A (en) * | 2006-07-12 | 2008-01-16 | 中国农业大学 | Soil infiltration performance real-time automatic measuring system |
CN103439237A (en) * | 2013-08-30 | 2013-12-11 | 山东大学 | Blocking test device and method for water-permeable concrete pile in piping soil foundation |
CN105785476A (en) * | 2016-05-17 | 2016-07-20 | 浙江水利水电学院 | High density electrical resistivity imaging method |
CN105938204A (en) * | 2016-06-27 | 2016-09-14 | 浙江水利水电学院 | Resistivity exploration method-based groundwater vulnerability assessment method |
CN206208821U (en) * | 2016-11-01 | 2017-05-31 | 西安理工大学 | A kind of native water feature and resistivity joint test device |
CN108204939A (en) * | 2017-12-15 | 2018-06-26 | 东南大学 | A kind of water penetration analysis survey meter for evaluating pollutant diffusion and method of work |
CN108267394A (en) * | 2018-03-09 | 2018-07-10 | 浙江省水利河口研究院 | A kind of earth-rock dam seepage farm monitoring system and its method for early warning |
-
2018
- 2018-11-01 CN CN201811295465.9A patent/CN109164033B/en not_active Expired - Fee Related
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5861750A (en) * | 1995-01-09 | 1999-01-19 | Anderson; Dennis M. | Geophysical methods and apparatus for determining the hydraulic conductivity of porous materials |
CN101105489A (en) * | 2006-07-12 | 2008-01-16 | 中国农业大学 | Soil infiltration performance real-time automatic measuring system |
CN103439237A (en) * | 2013-08-30 | 2013-12-11 | 山东大学 | Blocking test device and method for water-permeable concrete pile in piping soil foundation |
CN105785476A (en) * | 2016-05-17 | 2016-07-20 | 浙江水利水电学院 | High density electrical resistivity imaging method |
CN105938204A (en) * | 2016-06-27 | 2016-09-14 | 浙江水利水电学院 | Resistivity exploration method-based groundwater vulnerability assessment method |
CN206208821U (en) * | 2016-11-01 | 2017-05-31 | 西安理工大学 | A kind of native water feature and resistivity joint test device |
CN108204939A (en) * | 2017-12-15 | 2018-06-26 | 东南大学 | A kind of water penetration analysis survey meter for evaluating pollutant diffusion and method of work |
CN108267394A (en) * | 2018-03-09 | 2018-07-10 | 浙江省水利河口研究院 | A kind of earth-rock dam seepage farm monitoring system and its method for early warning |
Non-Patent Citations (2)
Title |
---|
G. E. ARCHIE: "The Electrical Resistivity Log as an Aid in Determining Some Reservoir Characteristics", 《DALLAS MEETING》 * |
周明耀 等: "基于孔隙分形维数的土壤大孔隙流水力特征参数研究", 《水科学进展》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109884131A (en) * | 2019-02-23 | 2019-06-14 | 机械工业勘察设计研究院有限公司 | One kind soil layer in submerging test of testing pits soaks lossless continuous probe method |
CN110836799A (en) * | 2019-11-15 | 2020-02-25 | 辽宁工程技术大学 | Suction control type rapid unsaturated soil sample preparation control system and method |
CN110988045A (en) * | 2019-12-24 | 2020-04-10 | 西安理工大学 | Soil body conductivity measuring method based on active heating optical fiber method |
CN111208042A (en) * | 2020-02-07 | 2020-05-29 | 中国科学院武汉岩土力学研究所 | Device and method for inverting hydraulic parameters of unsaturated waste soil |
CN111208042B (en) * | 2020-02-07 | 2021-06-29 | 中国科学院武汉岩土力学研究所 | Device and method for inverting hydraulic parameters of unsaturated waste soil |
CN112162014A (en) * | 2020-09-02 | 2021-01-01 | 塔里木大学 | Cotton field soil profile water data processing method based on electromagnetic induction data |
Also Published As
Publication number | Publication date |
---|---|
CN109164033B (en) | 2020-12-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109164033A (en) | A method of the acquisition unsaturated soil waterpower conductance based on resistivity method | |
Min et al. | Estimating groundwater recharge using deep vadose zone data under typical irrigated cropland in the piedmont region of the North China Plain | |
CN104375195B (en) | Many source multi-component three-dimensional joint inversion methods of time-frequency electromagnetism | |
Cui et al. | Modeling tree root diameter and biomass by ground-penetrating radar | |
Calamita et al. | Electrical resistivity and TDR methods for soil moisture estimation in central Italy test-sites | |
CN204832095U (en) | Alternation of wetting and drying crack measuring device of indoor soil property side slope model | |
WO2009117784A1 (en) | System, apparatus and method for measuring soil moisture content | |
Tromp-van Meerveld et al. | Assessment of multi-frequency electromagnetic induction for determining soil moisture patterns at the hillslope scale | |
Wang et al. | A laboratory study of the correlation between the thermal conductivity and electrical resistivity of soil | |
Bagarello et al. | Physical and hydraulic characterization of a clay soil at the plot scale | |
Zhao et al. | Electrical imaging of plant root zone: A review | |
Vanderborght et al. | Geophysical methods for field‐scale imaging of root zone properties and processes | |
CN106018970B (en) | The measurement method of Resistance of Grounding Grids in a kind of double-layer soil | |
Giambastiani et al. | Indirect root distribution characterization using electrical resistivity tomography in different soil conditions | |
Pan et al. | Estimating field-scale soil water dynamics at a heterogeneous site using multi-channel GPR | |
Descloitres et al. | Study of water tension differences in heterogeneous sandy soils using surface ERT | |
Batlle-Aguilar et al. | Axisymetrical infiltration in soil imaged by noninvasive electrical resistivimetry | |
Gu et al. | Soil salinity simulation based on electromagnetic induction and deep learning | |
JP2005337746A (en) | Electric exploration method | |
CN106321087B (en) | A method of obtaining lithostratigraphy factor | |
Gnecchi et al. | Soil water infiltration measurements using electrical impedance tomography | |
Zhang et al. | Application of resistivity measurement to stability evaluation for loess slopes | |
Wu et al. | Response of three-dimensional spatial variability of soil salinity to change of season of Xinjiang based on electromagnetic induction | |
CN102087093B (en) | Method for detecting soil mass cracking depth and detection device | |
CN106226218B (en) | A kind of method of determining tight sand circumferential direction permeability principal direction |
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: 20201222 Termination date: 20211101 |
|
CF01 | Termination of patent right due to non-payment of annual fee |