CN116858749B - Determination method for relationship between barrier curtain complex resistivity signal and permeability - Google Patents
Determination method for relationship between barrier curtain complex resistivity signal and permeability Download PDFInfo
- Publication number
- CN116858749B CN116858749B CN202310924264.5A CN202310924264A CN116858749B CN 116858749 B CN116858749 B CN 116858749B CN 202310924264 A CN202310924264 A CN 202310924264A CN 116858749 B CN116858749 B CN 116858749B
- Authority
- CN
- China
- Prior art keywords
- barrier curtain
- permeability
- complex resistivity
- barrier
- representing
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 230000004888 barrier function Effects 0.000 title claims abstract description 81
- 230000035699 permeability Effects 0.000 title claims abstract description 53
- 238000000034 method Methods 0.000 title claims abstract description 36
- 239000000463 material Substances 0.000 claims abstract description 65
- 238000007596 consolidation process Methods 0.000 claims abstract description 50
- 238000001228 spectrum Methods 0.000 claims abstract description 23
- 238000012360 testing method Methods 0.000 claims abstract description 6
- 230000004044 response Effects 0.000 claims abstract description 5
- 230000000903 blocking effect Effects 0.000 claims description 18
- 238000005070 sampling Methods 0.000 claims description 18
- 238000012544 monitoring process Methods 0.000 claims description 11
- 238000005553 drilling Methods 0.000 claims description 9
- 239000002245 particle Substances 0.000 claims description 8
- 238000005266 casting Methods 0.000 claims description 4
- 238000002247 constant time method Methods 0.000 claims description 4
- 230000010287 polarization Effects 0.000 claims description 4
- 238000012937 correction Methods 0.000 claims description 3
- 210000003041 ligament Anatomy 0.000 claims description 3
- 238000007789 sealing Methods 0.000 claims description 3
- 238000005259 measurement Methods 0.000 abstract description 5
- 238000006243 chemical reaction Methods 0.000 abstract description 3
- 238000004445 quantitative analysis Methods 0.000 abstract description 2
- 238000003860 storage Methods 0.000 description 7
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 5
- 229910000278 bentonite Inorganic materials 0.000 description 4
- 239000000440 bentonite Substances 0.000 description 4
- 238000001514 detection method Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- 238000002347 injection Methods 0.000 description 3
- 239000007924 injection Substances 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 239000000243 solution Substances 0.000 description 3
- 230000003595 spectral effect Effects 0.000 description 3
- 238000003756 stirring Methods 0.000 description 3
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 description 2
- 238000002441 X-ray diffraction Methods 0.000 description 2
- SVPXDRXYRYOSEX-UHFFFAOYSA-N bentoquatam Chemical compound O.O=[Si]=O.O=[Al]O[Al]=O SVPXDRXYRYOSEX-UHFFFAOYSA-N 0.000 description 2
- 238000009933 burial Methods 0.000 description 2
- 238000002591 computed tomography Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000005336 cracking Methods 0.000 description 1
- 239000003344 environmental pollutant Substances 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- QSHDDOUJBYECFT-UHFFFAOYSA-N mercury Chemical compound [Hg] QSHDDOUJBYECFT-UHFFFAOYSA-N 0.000 description 1
- 229910052753 mercury Inorganic materials 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 229910052757 nitrogen Inorganic materials 0.000 description 1
- 238000001225 nuclear magnetic resonance method Methods 0.000 description 1
- 231100000719 pollutant Toxicity 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 238000005204 segregation Methods 0.000 description 1
- 239000002002 slurry Substances 0.000 description 1
- 239000002689 soil Substances 0.000 description 1
- 238000007711 solidification Methods 0.000 description 1
- 230000008023 solidification Effects 0.000 description 1
- 238000001179 sorption measurement Methods 0.000 description 1
- 238000005507 spraying Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000001502 supplementing effect Effects 0.000 description 1
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
-
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Landscapes
- Chemical & Material Sciences (AREA)
- Biochemistry (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Analytical Chemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Dispersion Chemistry (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Electrochemistry (AREA)
- Investigating Or Analyzing Materials By The Use Of Electric Means (AREA)
Abstract
The invention discloses a method for determining the relation between a barrier curtain complex resistivity signal and permeability, which comprises the following steps: obtaining a barrier curtain material sample for microstructure test, and analyzing microstructure characteristics of the barrier curtain material; acquiring complex resistivity response data of the barrier curtain under different consolidation states; based on the structural characteristic parameters of the barrier curtain material, establishing a relation model of the barrier curtain complex resistivity spectrum parameter and the permeability coefficient; and taking the microstructure parameters as the bands, and evaluating the permeability characteristics of the barrier curtain material in different consolidation states according to the complex resistivity spectrum parameters. According to the invention, the difference of permeability coefficients of the barrier curtain in different consolidation states is considered, time-shifting complex resistivity measurement is carried out on the barrier material, a quantitative analysis method of the permeability of the barrier curtain material based on relaxation time of complex resistivity spectrum characteristic parameters in different consolidation states is provided, conversion from the relaxation time to the permeability coefficient is completed, and a novel method is provided for evaluating the permeability of the barrier curtain.
Description
Technical Field
The invention relates to the technical field of pollution site investigation, in particular to a method for determining the relation between a barrier curtain complex resistivity signal and permeability.
Background
The blocking technology is an important means for repairing the underground water of the polluted land. The conventional blocking mode mainly controls the migration of pollutants by paving blocking curtains, and different quality problems often occur in the construction operation process of different types of blocking curtains, such as wall cracking, hollowness and the like of plastic concrete blocking curtains, wall continuity poor of deep stirring concrete curtains and wall segregation and overhead problems of high-pressure spraying curtains. Aiming at various quality problems possibly occurring, quality detection needs to be carried out on the barrier curtain, so that the integrity of the wall body and the exertion of the seepage-proofing function are ensured.
The traditional method mainly takes drilling sampling and water injection experiment research on the permeability characteristics, has the problems of long time consumption, high cost, easiness in damaging a wall body, incapability of realizing long-time monitoring and the like, utilizes geophysical exploration technology to detect the barrier curtain, is lossless and rapid, can realize continuous detection, and can perform time-shifting complex resistivity measurement on the consolidation characteristics of the barrier curtain to acquire the permeability characteristics of the barrier curtain in different consolidation states.
Disclosure of Invention
The invention aims to overcome the defects in the prior art, and provides a method for determining the relation between a barrier curtain complex resistivity signal and permeability.
The first aspect of the invention provides a method for determining the relation between a barrier curtain complex resistivity signal and permeability, which comprises the following steps:
Obtaining a barrier curtain material sample, performing microstructure test on the barrier curtain material sample, and analyzing microstructure characteristics of the barrier curtain material;
Acquiring complex resistivity response data of the barrier curtain under different consolidation states;
Based on the structural characteristic parameters of the barrier curtain material, establishing a relation model of the barrier curtain complex resistivity spectrum parameter and the permeability coefficient;
and taking the microstructure parameters as a ligament, and evaluating the permeability characteristics of the barrier curtain material in different consolidation states by using the complex resistivity spectrum parameters.
In the scheme, the blocking curtain material sample is grouting filling material fully stirred before pouring or is obtained after pouring through sampling and drilling, drilling and sampling are performed immediately after pouring is completed, the sampling position is right above the blocking curtain, the sampling depth is 2-4 m, obvious damage to the curtain cannot be caused, hole sealing is performed in time after sampling, and the influence on the anti-seepage effect of the blocking curtain material sample is avoided.
In this scheme, the parameters that characterize the microstructure characteristics of the material include: particle diameter, porosity.
In this scheme, the different consolidation states are determined according to the following formula:
；
Wherein: Representing a monitoring period (h); /(I) Representing different consolidation phases,/>Take the value of/>N is the consolidation completion time; /(I)Representing the temperature equivalent coefficient,/>，/>Representing the temperature of the casting site,/>The standard temperature is represented, and the value is 20 ℃; /(I)The correction coefficient is generally 1.02 to 1.04.
In the scheme, the complex resistivity signal is acquired by a cross-hole resistivity CT method, and resistivity and polarization rate signals are acquired at the same time; the acquisition method is a time-shifting complex resistivity method, namely, the same data acquisition system is adopted at the same position to acquire complex resistivity signals at different times.
In this embodiment, the porosity of the material in the different consolidation states is determined by the following formula:
；
Wherein: Representing the initial porosity,/> Representing different consolidation phases,/>Take the value of/>N is the consolidation completion time.
In this scheme, the step of establishing the relationship model is:
Step S1: the complex resistivity parameter is obtained using the following formula:
；
Wherein: Representing complex resistivity; /(I) The resistivity of the porous medium when the power supply time is infinitely long is represented; /(I)Representing the polarizability of the porous medium; /(I)Representing a relaxation time; /(I)Represents angular frequency; /(I)Representing the frequency correlation coefficient/>；/>Representing imaginary units;
extracting complex resistivity spectrum characteristic parameters such as relaxation time and the like;
step S2: based on microstructure parameters of the early-stage blocking porous medium material, the following formula is utilized:
；
Wherein: Represents the permeability coefficient; /(I) Represents particle diameter; /(I)Indicating a cementation index; /(I)Representing porosity;
acquiring the relation between the microstructure parameter of the barrier porous medium material and the permeability coefficient k;
And combining the early-stage complex resistivity data acquisition result, and based on the formula:
；
acquiring the relation between the complex resistivity spectrum parameter and the microstructure of the material;
the relationship between microstructure parameters of the barrier porous medium material and relaxation time of corresponding parameters of complex resistivity spectrum characteristics is summarized as follows:
；
Step S3: based on the steps S1-S2, taking the structural characteristic parameters of the barrier porous medium material as the bands, solving the permeability coefficients of the barrier curtain material in different consolidation states through the established relation model of the complex resistivity spectrum parameters and the permeability coefficients, and evaluating the permeability characteristics of the barrier curtain material.
The invention solves the defects existing in the background technology, and has the following beneficial effects:
According to the invention, a polluted site grouting blocking porous medium material is taken as a research object, a microstructure test is carried out on a blocking material sample based on the microstructure and complex resistivity theory of the blocking porous medium material, time-shifting complex resistivity measurement is carried out aiming at different permeability of blocking curtains in different consolidation states, a complex resistivity signal response research of blocking curtain materials in different consolidation states is carried out by fitting complex resistivity measurement data and a frequency domain Cole-Cole model through a least square method, characteristic spectral parameter relaxation time is extracted, microstructure characteristics of the material are taken as a ligament, a quantitative relation between the complex resistivity spectral characteristic parameter relaxation time and the permeability coefficient is established, a quantitative analysis method of the blocking material based on the permeability of the complex resistivity spectral characteristic parameter relaxation time in different consolidation states is provided, conversion from the relaxation time to the permeability coefficient is completed, and a novel method is provided for evaluating the blocking curtain permeability. Compared with the method for analyzing the curtain permeability characteristics by utilizing the traditional drilling sampling and water injection experiments, the method has the advantages of nondestructive detection, high working efficiency, convenience in operation and capability of realizing continuous detection.
Drawings
In order to more clearly illustrate the technical solutions of embodiments or examples of the present invention, the drawings that are required to be used in the embodiments or examples of the present invention will be briefly described below, and it is apparent that the drawings in the following description are only some embodiments of the present invention, and other drawings may be obtained according to the drawings without inventive efforts for those skilled in the art.
FIG. 1 is a flow chart of a method of determining the relationship between barrier curtain complex resistivity signal and permeability according to the present invention;
FIG. 2 is a schematic diagram of a method of determining the relationship between barrier curtain complex resistivity signal and permeability according to the present invention;
FIG. 3 shows a plot of monitoring period versus different consolidation phases in an embodiment of the present invention;
The achievement of the objects, functional features and advantages of the present invention will be further described with reference to the accompanying drawings, in conjunction with the embodiments.
Detailed Description
In order that the above-recited objects, features and advantages of the present application will be more clearly understood, a more particular description of the application will be rendered by reference to the appended drawings and appended detailed description. It should be noted that, without conflict, the embodiments of the present application and features in the embodiments may be combined with each other.
In the following description, numerous specific details are set forth in order to provide a thorough understanding of the present invention, but the present invention may be practiced in other ways than those described herein, and therefore the scope of the present invention is not limited to the specific embodiments disclosed below.
FIGS. 1 and 2 are a flow chart and a schematic diagram of a method for determining the relationship between the barrier curtain complex resistivity signal and permeability according to the present invention.
As shown in fig. 1, a first aspect of the present invention provides a method for determining a relationship between a barrier curtain complex resistivity signal and permeability, including:
S102, obtaining a barrier curtain material sample, performing microstructure test on the barrier curtain material sample, and analyzing microstructure characteristics of the barrier curtain material;
It should be noted that, the sample of the barrier curtain material is a grouting filling material which is fully stirred before pouring, or is obtained after pouring is completed through sampling and drilling, only the characteristics of the barrier curtain material need to be reflected, drilling and sampling should be immediately carried out after pouring is completed, the sampling position is right above the barrier curtain, the sampling depth is preferably 2-4 m, obvious damage to the curtain cannot be caused, and the sealing of holes should be carried out in time after sampling, so that the anti-seepage effect of the barrier curtain is prevented from being influenced.
In the embodiment of the invention, typical grouting barrier curtain types of the polluted site are soil-bentonite wall and concrete-bentonite wall. When the bentonite is hydrated, the proportion of the bentonite to the water is 1: and 7, standing for 24 hours after stirring is completed. And then, respectively configuring the barrier curtain materials by using undisturbed soil and concrete of the field, and completing pouring by using the stirring piles and the jet piles.
The microstructure characteristics of the sample can be obtained by adopting a CT scanning technology, a nuclear magnetic resonance method, a low-temperature nitrogen adsorption method, a mercury injection method and the like. The parameters characterizing the microstructure features of the material include: particle diameter, porosity.
In the embodiment of the invention, a CT scanning method is adopted to carry out X-ray diffraction test on the barrier curtain material, the structure scanning is carried out on the 1/3 height of the middle part of the sample, the binarization processing is carried out on the image by using threshold segmentation after the scanning is finished, and the denoising is carried out on the obtained CT image.
S104, collecting response data of complex resistivity of the barrier curtain under different consolidation states;
it should be noted that the different consolidation states are determined according to the following formula:
；
Wherein: Representing a monitoring period (h); /(I) Representing different consolidation phases,/>Take the value of/>N is the consolidation complete time (N is the consolidation complete time, where 1 represents the consolidation time of 1 day, which is the first consolidation stage, 2 represents the next consolidation stage after completion of monitoring the first consolidation stage, and so on); /(I)Representing the temperature equivalent coefficient,/>，/>Representing the temperature of the casting site,/>The standard temperature is represented, and the value is 20 ℃; /(I)The correction coefficient is generally 1.02 to 1.04.
In the embodiment of the invention, the value of a is 1.02, the casting site temperature is 25 ℃, and the monitoring period is carried outDifferent consolidation phases/>The relationship between these is shown in fig. 3, and the monitoring period is gradually longer as the consolidation phase progresses, that is, the monitoring frequency is gradually reduced, and the monitoring period is monitored every 15 days at the 15 th consolidation phase.
Specifically, the complex resistivity signals are acquired by a cross-hole resistivity CT method, the acquisition method is a time-shifting complex resistivity method, namely, the same data acquisition system is adopted at the same position to acquire the complex resistivity signals with different time, the acquired signals are resistivity and polarization rate signals, and the acquisition time of each time is the monitoring period T.
In the embodiment of the invention, the inter-hole resistivity CT measurement electrode spacing is 0.5m, 18 single measuring line electrodes are provided, the total measuring line length is 8.5m, the first electrode burial depth is determined according to the water level burial depth and is 2.5m, namely the drilling hole depth is 11m. The measuring array adopts an AM-BN type array, and data inversion is carried out through Res2dinv software, and the maximum iteration number is 7.
S106, establishing a relation model of the complex resistivity spectrum parameter of the barrier curtain and the permeability coefficient based on the structural characteristic parameter of the barrier curtain material;
in an embodiment of the invention, both the particle diameter and the initial porosity have been obtained by X-ray diffraction. As the particle diameter is unchanged and the porosity is gradually reduced when the curtain is solidified, the porosity of the barrier curtain in different solidification states needs to be calculated.
It should be noted that the porosity of the material in the different consolidation states is determined by the following formula:
；
Wherein: Representing the initial porosity,/> Representing different consolidation phases,/>Take the value of/>N is the consolidation completion time.
From the initial porosityAnd the porosity/>, under different consolidation states, is calculated by the formula。
S108, taking the microstructure parameters as the bands, and evaluating the permeability characteristics of the barrier curtain material in different consolidation states by using the complex resistivity spectrum parameters.
It should be noted that, the step of establishing the relationship model is as follows:
Step S1: the complex resistivity parameter is obtained using the following formula:
；
Wherein: representing complex resistivity/> ；/>Representing the resistivity of the porous medium at an infinitely long power supply time；/>Representing the polarizability of the porous medium; /(I)Expressed as relaxation time/>；
Represent angular frequency/>；/>Representing the frequency correlation coefficient/>；/>Representing imaginary units;
extracting complex resistivity spectrum characteristic parameters such as relaxation time and the like;
step S2: based on microstructure parameters of the early-stage blocking porous medium material, the following formula is utilized:
；
Wherein: Represents the permeability coefficient; /(I) Represents particle diameter; /(I)Indicating a cementation index; /(I)Representing porosity;
acquiring the relation between the microstructure parameter of the barrier porous medium material and the permeability coefficient k;
And combining the early-stage complex resistivity data acquisition result, and based on the formula:
；
acquiring the relation between the complex resistivity spectrum parameter and the microstructure of the material;
Based on the formula, the relationship between the microstructure parameters of the barrier porous medium material and the relaxation time of the corresponding parameters of the complex resistivity spectrum characteristics is summarized as follows:
；
Step S3: based on the steps S1-S2, taking the structural characteristic parameters of the barrier porous medium material as the bands, solving the permeability coefficients of the barrier curtain material in different consolidation states through the established relation model of the complex resistivity spectrum parameters and the permeability coefficients, and evaluating the permeability characteristics of the barrier curtain material.
In the embodiment of the invention, complex resistivity parameters are fitted by using a Cole-Cole model, and characteristic parameters such as complex resistivity, zero frequency resistivity, polarization rate and relaxation time are obtained. And (3) according to the obtained cementation index, wherein q and p are cementation constants and are porosities. The consolidation states are different, the acquired complex resistivity spectrum parameters and the porosity of the barrier curtain are different, the permeability coefficient of the barrier curtain in different consolidation states is obtained according to the acquired relaxation time in different consolidation states and the porosity in different consolidation states calculated by using a formula, and the conversion from the relaxation time of the complex resistivity spectrum characteristic parameters to the permeability coefficient is completed. If the permeability coefficient is smaller than cm/s, the barrier curtain is considered to still play a role, otherwise, the barrier curtain is required to be repaired and treated through plugging, slurry supplementing and other operations.
In the several embodiments provided by the present application, it should be understood that the disclosed apparatus and method may be implemented in other ways. The above described device embodiments are only illustrative, e.g. the division of the units is only one logical function division, and there may be other divisions in practice, such as: multiple units or components may be combined or may be integrated into another system, or some features may be omitted, or not performed. In addition, the various components shown or discussed may be coupled or directly coupled or communicatively coupled to each other via some interface, whether indirectly coupled or communicatively coupled to devices or units, whether electrically, mechanically, or otherwise.
The units described above as separate components may or may not be physically separate, and components shown as units may or may not be physical units; can be located in one place or distributed to a plurality of network units; some or all of the units may be selected according to actual needs to achieve the purpose of the solution of this embodiment.
In addition, each functional unit in each embodiment of the present invention may be integrated in one processing unit, or each unit may be separately used as one unit, or two or more units may be integrated in one unit; the integrated units may be implemented in hardware or in hardware plus software functional units.
Those of ordinary skill in the art will appreciate that: all or part of the steps for implementing the above method embodiments may be implemented by hardware related to program instructions, and the foregoing program may be stored in a computer readable storage medium, where the program, when executed, performs steps including the above method embodiments; and the aforementioned storage medium includes: a mobile storage device, a Read-Only Memory (ROM), a random access Memory (RAM, random Access Memory), a magnetic disk or optical disk, or the like, which can store program codes.
Or the above-described integrated units of the invention may be stored in a computer-readable storage medium if implemented in the form of software functional modules and sold or used as separate products. Based on such understanding, the technical solutions of the embodiments of the present invention may be embodied in essence or a part contributing to the prior art in the form of a software product stored in a storage medium, including several instructions for causing a computer device (which may be a personal computer, a server, or a network device, etc.) to execute all or part of the methods described in the embodiments of the present invention. And the aforementioned storage medium includes: a removable storage device, ROM, RAM, magnetic or optical disk, or other medium capable of storing program code.
The foregoing is merely illustrative of the present invention, and the present invention is not limited thereto, and any person skilled in the art will readily recognize that variations or substitutions are within the scope of the present invention. Therefore, the protection scope of the present invention shall be subject to the protection scope of the claims.
Claims (5)
1. A method for determining the relationship between a barrier curtain complex resistivity signal and permeability, comprising the steps of:
Obtaining a barrier curtain material sample, performing microstructure test on the barrier curtain material sample, and analyzing microstructure characteristics of the barrier curtain material;
Acquiring complex resistivity response data of the barrier curtain under different consolidation states;
Based on the structural characteristic parameters of the barrier curtain material, establishing a relation model of the barrier curtain complex resistivity spectrum parameter and the permeability coefficient;
taking microstructure parameters as a ligament, and evaluating permeability characteristics of the barrier curtain material in different consolidation states by using complex resistivity spectrum parameters;
Parameters characterizing the microstructure characteristics of the barrier curtain material sample include: particle diameter, porosity;
porosity of material under different consolidation conditions Is determined by the following formula:
；
Wherein, Representing the initial porosity,/>Representing different consolidation phases,/>Take the value of/>N is the consolidation completion time;
The relation model building step is as follows:
Step S1: the complex resistivity parameter is obtained using the following formula:
；
Wherein: Representing complex resistivity; /(I) The resistivity of the porous medium when the power supply time is infinitely long is represented; /(I)Representing the polarizability of the porous medium; /(I)Representing a relaxation time; /(I)Represents angular frequency; /(I)Representing the frequency correlation coefficient/>；/>Representing imaginary units;
extracting complex resistivity spectrum characteristic parameters;
step S2: based on microstructure parameters of the early-stage blocking porous medium material, the following formula is utilized:
；
Wherein: Represents the permeability coefficient; /(I) Represents particle diameter; /(I)Indicating a cementation index; /(I)Representing porosity;
acquiring the relation between the microstructure parameter of the barrier porous medium material and the permeability coefficient k;
And combining the early-stage complex resistivity data acquisition result, and based on the formula:
；
acquiring the relation between the complex resistivity spectrum parameter and the microstructure of the material;
the relationship between microstructure parameters of the barrier porous medium material and relaxation time of corresponding parameters of complex resistivity spectrum characteristics is summarized as follows:
；
Step S3: based on the steps S1-S2, taking the structural characteristic parameters of the barrier porous medium material as the bands, solving the permeability coefficients of the barrier curtain material in different consolidation states through the established relation model of the complex resistivity spectrum parameters and the permeability coefficients, and evaluating the permeability characteristics of the barrier curtain material.
2. The method for determining the relation between the complex resistivity signal and the permeability of the barrier curtain according to claim 1, wherein the barrier curtain material sample is a grouting filling material which is fully stirred before pouring or is obtained after pouring through a sampling drilling hole, the characteristic of the barrier curtain material only needs to be reflected, sampling by the sampling drilling hole is performed immediately after pouring, the sampling position is right above the barrier curtain, the sampling depth is 2-4 m, and hole sealing is performed in time after sampling.
3. The method of determining the relationship between barrier curtain complex resistivity signal and permeability of claim 1, wherein the different consolidation states are determined according to the following equation:
；
Wherein: representing a monitoring period; /(I) Representing different consolidation phases,/>Take the value of/>N is the consolidation completion time; /(I)Representing the temperature equivalent coefficient,/>，/>Representing the temperature of the casting site,/>The standard temperature is represented, and the value is 20 ℃; /(I)Representing a correction coefficient, and taking 1.02-1.04; /(I)Representing porosity.
4. The method for determining the relationship between a barrier curtain complex resistivity signal and permeability of claim 1, wherein the complex resistivity signal is acquired by a cross-hole resistivity CT method while the resistivity and polarization rate signals are acquired.
5. The method for determining the relation between the barrier curtain complex resistivity signal and the permeability according to claim 4, wherein the complex resistivity signal is acquired by a cross-hole resistivity CT method, and the complex resistivity signals with different times are acquired by adopting the same data acquisition system at the same position for a time-shifting complex resistivity method.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310924264.5A CN116858749B (en) | 2023-07-26 | 2023-07-26 | Determination method for relationship between barrier curtain complex resistivity signal and permeability |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310924264.5A CN116858749B (en) | 2023-07-26 | 2023-07-26 | Determination method for relationship between barrier curtain complex resistivity signal and permeability |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN116858749A CN116858749A (en) | 2023-10-10 |
| CN116858749B true CN116858749B (en) | 2024-04-23 |
Family
ID=88221497
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310924264.5A Active CN116858749B (en) | 2023-07-26 | 2023-07-26 | Determination method for relationship between barrier curtain complex resistivity signal and permeability |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116858749B (en) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110988999A (en) * | 2019-12-02 | 2020-04-10 | 山东大学 | Detection method and system for analyzing pile foundation based on cross-hole resistivity method CT inversion imaging |
| CN111535373A (en) * | 2020-05-12 | 2020-08-14 | 中建六局水利水电建设集团有限公司 | Deep foundation pit water-proof curtain leakage detection method |
| CN112985718A (en) * | 2021-02-05 | 2021-06-18 | 中建一局集团建设发展有限公司 | Waterproof curtain evaluation system and method based on high-density resistivity method |
| CN114993918A (en) * | 2022-08-03 | 2022-09-02 | 石家庄铁道大学 | Triaxial pressure chamber soil sample permeability coefficient measuring method |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| FR2884923B1 (en) * | 2005-04-26 | 2007-10-19 | Inst Francais Du Petrole | METHOD AND DEVICE FOR EVALUATING FLOW PARAMETERS AND ELECTRIC PARAMETERS OF A POROUS MEDIUM. |
| US7825659B2 (en) * | 2005-06-03 | 2010-11-02 | Baker Hughes Incorporated | Pore-scale geometric models for interpretation of downhole formation evaluation data |
| US7324899B2 (en) * | 2005-07-22 | 2008-01-29 | University Of Utah | Geophysical technique for mineral exploration and discrimination based on electromagnetic methods and associated systems |
-
2023
- 2023-07-26 CN CN202310924264.5A patent/CN116858749B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110988999A (en) * | 2019-12-02 | 2020-04-10 | 山东大学 | Detection method and system for analyzing pile foundation based on cross-hole resistivity method CT inversion imaging |
| CN111535373A (en) * | 2020-05-12 | 2020-08-14 | 中建六局水利水电建设集团有限公司 | Deep foundation pit water-proof curtain leakage detection method |
| CN112985718A (en) * | 2021-02-05 | 2021-06-18 | 中建一局集团建设发展有限公司 | Waterproof curtain evaluation system and method based on high-density resistivity method |
| CN114993918A (en) * | 2022-08-03 | 2022-09-02 | 石家庄铁道大学 | Triaxial pressure chamber soil sample permeability coefficient measuring method |
Non-Patent Citations (2)
| Title |
|---|
| 关联帷幕防渗性能变化的重力坝结构稳定评价方法;彭琦;高大水;高江林;张宇驰;;水利与建筑工程学报(第03期);68-71 * |
| 含甲烷水合物多孔介质的复电阻率模型对比;王彩程;陈强;邢兰昌;刘昌岭;郑金吾;;海洋地质前沿(第05期);67-73 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN116858749A (en) | 2023-10-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Kuras et al. | Monitoring hydraulic processes with automated time-lapse electrical resistivity tomography (ALERT) | |
| Barker et al. | The application of time-lapse electrical tomography in groundwater studies | |
| CN110735635B (en) | Method for determining content of lost gas in shale gas content test | |
| EP2543769B1 (en) | Method of consolidating a foundation soil | |
| Lu et al. | Using electrical resistivity tomography and surface nuclear magnetic resonance to investigate cultural relic preservation in Leitai, China | |
| CN110208487A (en) | A kind of Shale Hydration damage measure method based on CT scan | |
| Muñoz-Carpena et al. | TDR estimation of electrical conductivity and saline solute concentration in a volcanic soil | |
| CN107300521A (en) | The assay method and measure device of swelled ground cracking water content | |
| CN105158108A (en) | Method for calibrating capacitive soil moisture sensor on basis of undisturbed soil columns | |
| RU2467315C1 (en) | Method to detect spatial distribution and concentration of clay in core sample | |
| CN102183447A (en) | Test system and test method for permeability coefficients of aquifer | |
| CN116858749B (en) | Determination method for relationship between barrier curtain complex resistivity signal and permeability | |
| Shi et al. | History matching of CO2 core flooding CT scan saturation profiles with porosity dependent capillary pressure | |
| Pineda et al. | Performance of an innovative direct-push piston sampler in soft clay | |
| CN207066937U (en) | The measure device of swelled ground cracking water content | |
| US10260316B2 (en) | Facility for monitoring geological gas storage sites | |
| CN111505041A (en) | Shale oil saturation well site rapid evaluation method | |
| CN117033864A (en) | Method for calculating critical depth of instability of expansive soil pit wall under excavation unloading and predicting deformation | |
| CN109613041B (en) | Frozen soil permeability coefficient determination system and method based on nuclear magnetic resonance | |
| Pineda et al. | Testing block samples from silty deposits | |
| CN107121336A (en) | The method for building up of Sandy Silt uniaxial compressive strength curve | |
| Dumont et al. | The effect of initial water distribution and spatial resolution on the interpretation of ERT monitoring of water infiltration in a landfill cover | |
| CN113916193A (en) | Method for calculating hydrogeological parameters of aquifer by inversion | |
| CN116609242B (en) | Barrier material permeability characteristic evaluation method and system based on resistivity signal | |
| Quinnan et al. | Hydrostratigraphic and permeability profiling for groundwater remediation projects |
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 |