CN110471122B - Three-dimensional resistivity depth sounding application method based on high-density electrical method Wener device - Google Patents

Three-dimensional resistivity depth sounding application method based on high-density electrical method Wener device Download PDF

Info

Publication number
CN110471122B
CN110471122B CN201910784248.4A CN201910784248A CN110471122B CN 110471122 B CN110471122 B CN 110471122B CN 201910784248 A CN201910784248 A CN 201910784248A CN 110471122 B CN110471122 B CN 110471122B
Authority
CN
China
Prior art keywords
dimensional
grid
resistivity
density
section
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
Application number
CN201910784248.4A
Other languages
Chinese (zh)
Other versions
CN110471122A (en
Inventor
李忠平
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Individual
Original Assignee
Individual
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Individual filed Critical Individual
Priority to CN201910784248.4A priority Critical patent/CN110471122B/en
Publication of CN110471122A publication Critical patent/CN110471122A/en
Application granted granted Critical
Publication of CN110471122B publication Critical patent/CN110471122B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V3/00Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
    • G01V3/38Processing data, e.g. for analysis, for interpretation, for correction

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

The invention discloses a three-dimensional resistivity sounding application method based on a high-density electrical method wenna device, which is characterized in that a three-dimensional resistivity data body is constructed by using a high-density 2D resistivity instrument wenna device and a large-network-degree large-polar-distance electrical sounding profile, an electrode arrangement is converted into a symmetrical quadrupole device, the symmetrical quadrupole device is reconstructed into an uneven measuring network, and a least square method is used for carrying out three-dimensional inversion.

Description

Three-dimensional resistivity depth sounding application method based on high-density electrical method Wener device
Technical Field
The invention relates to the technical field of geological exploration, in particular to a three-dimensional resistivity sounding application method based on a high-density electrical method Wener device.
Background
The high-resolution geological mineral three-dimensional exploration is a target for electrical method development, and the collection of high-resolution electrical method three-dimensional observation data and the development of three-dimensional inversion are the keys of the detection of a three-dimensional fine structure of resistivity in a measurement area. At present, a resistivity sounding three-dimensional observation mode is an important aspect of research, and the resistivity sounding three-dimensional observation mode not only relates to the efficiency and cost of field observation, but also relates to the resolution and inversion effect of data.
The existing data show that the common resistivity three-dimensional inversion can only meet the requirements of shallow three-dimensional exploration and is mainly used for shallow hydrological, engineering and environmental exploration, Loke provides a simple three-dimensional observation mode, under the condition that the resolution is not influenced, the data points are reduced by one third, so that the field workload is reduced, the calculated amount is also reduced, the three-dimensional inversion interpretation can be implemented only by an 80486DX2/66 microcomputer, and the application prospect is good. From the perspective of satisfying the mineral resource deep exploration, the common three-dimensional exploration depth is difficult to satisfy the requirements, and the three-dimensional inversion by using the large polar distance electrical sounding data is more feasible. Because the electrical depth data measured on the plane is three-dimensional, the field observation is easy to implement, and the existing large amount of electrical depth data can be reused.
3DRES was used to process three-dimensional resistivity imaging measurements (Li and Oldenburg 1992, White at al.2001) data, which enabled automatic formation of three-dimensional resistivity models based on the measured data. In this type of measurement, the electrodes are arranged in a rectangular grid. It is emphasized that three-dimensional resistivity imaging measurements are not merely superimposed from a series of two-dimensional data, but rather are mature three-dimensional inversion methods, which have their own application features. Three-dimensional electrode arrangements such as pole-pole, pole-dipole, and dipole-dipole, etc., are often used in practice. Other arrangements are rarely used because of less effective data coverage. When the computer has 1.5GB RAM, the program can support 77 x 77 (or 5929) electrode points (Loke 2002).
Experiments prove that when a simple 2D method (or 1D electrical depth) is not good in effect, the high-density electrical method instrument supporting multiple parallel cables and a multiple covering technology are used for converting 2D measurement into 3D measurement, and an obvious shallow three-dimensional inversion model can be obtained.
In order to reduce the cost, increase the detection depth, effectively reflect the spatial characteristics of a geologic body and provide detailed information for the arrangement of drill holes, the invention provides a three-dimensional resistivity sounding application method based on a high-density electrical method Wener device.
Disclosure of Invention
The invention mainly solves the technical problem of how to provide a three-dimensional resistivity sounding application method based on a high-density electrical method Wener device, a three-dimensional resistivity data body is constructed by using a high-density 2D resistivity instrument Wener device and a large-network-degree large-polar-distance electrical sounding profile, electrode arrangement is converted into a symmetrical quadrupole device, the symmetrical quadrupole device is reconstructed into an uneven survey network, and three-dimensional inversion is carried out by using a least square method.
In order to solve the technical problems, the invention adopts a technical scheme that: the three-dimensional resistivity depth measurement application method based on the high-density electrical method Wener device is provided, and the specific method is as follows:
1) intercepting a data volume along a rectangular section at equal intervals by adopting a 2D high-density resistivity instrument;
2) constructing a three-dimensional uneven grid electrode arrangement method and a grid on the basis of the existing 2D high-density resistivity profile data;
3) and (4) converting two-dimensional data formats by combining the requirements of a three-dimensional inversion program, and compiling a 2D high-density electrical method non-uniform network measurement three-dimensional inversion data format.
Preferably, the specific interception method comprises: 6 high-density electrical method sections of 2D are arranged in the selected area, the sections are numbered as 11, 12, 13, 14, 15 and 16, the section spacing is 50m, the section point spacing is 10m, the section orientation is 295 degrees, the observation data points of the inverted trapezoid of the 2D high-density electrical method section wenner device are intercepted and adjusted to be rectangular data points, the starting point and the ending point are 38 and 62 respectively, the section length is 240m, the minimum AB distance is 60m, and the maximum AB distance is 660 m.
Preferably, the method for constructing the three-dimensional uneven grid electrode arrangement and the specific method for constructing the grid are as follows:
converting the high-density electrical method Wener arrangement data format into a symmetric quadrupole electrical depth AB and MN polar distance relationship as follows:
AB/2:30m,60m,90m,120m,150m,180m,210m,240m,270m,300m,330m
MN/2:10m,20m,30m,40m,50m,60m,70m,80m,90m,100m,110m
the three-dimensional electrical sounding measures the uneven grid spacing, each section electrical sounding point corresponds to A, B, M, N electrode positions, and the arrangement is as follows: taking the number 38 of the initial sounding point of the No. 16 section as the origin of coordinates 0 and 0, taking the position of an A, B, M, N electrode corresponding to each sounding point as the abscissa, taking the equal intervals of 0m, 50m, 100m, 150m, 200m and 250m of the No. 11, 12, 13, 14, 15 and 16 sections as the ordinate to establish the corresponding relation between all electrode points and resistivity and elevation, wherein the grid interval of the abscissa is as follows:
Figure GDA0002978120390000031
Figure GDA0002978120390000041
the grid spacing is symmetrically arranged as follows: 0m 20m 10m 10m … … 10m 20m
The vertical coordinate grid distance is: 50 m.
Preferably, the writing of the 2D high-density electrical method non-uniform measurement network three-dimensional inversion data format is as follows: the non-uniform grid electrode arrangement and grid spacing are combined with the three-dimensional inversion program requirement, non-uniform grid three-dimensional inversion data format conversion software is compiled, and the high-density Wener data format is converted into the non-uniform grid three-dimensional inversion data format as follows:
1 Block
89 6
Nonuniform grid
x-location of grid-lines
0 20 30 40……870 880 900
y-location of grid-lines
0 50 100 150 200 250
7
858
300 0 360 0 320 0 340 0 251
270 0 390 0 310 0 350 0 327
240 0 420 0 300 0 360 0 392
……………………………………………………
Topography
2
244 244 243………………………………221
0
0
0
0
0。
the invention has the following beneficial effects: the high-density 2D resistivity instrument temperature-sodium device and the large-network-degree large-polar-distance electrical depth profile are used for building a three-dimensional resistivity data body, the electrodes are arranged and converted into a symmetrical quadrupole device, the symmetrical quadrupole device is built into an uneven measurement network, three-dimensional inversion is carried out by using a least square method, the detection depth is increased by the observation method, the cost is reduced, the spatial characteristics of a geologic body are effectively reflected, detailed data are provided for borehole arrangement, the uneven measurement network high-density electrical method observation system is suitable for a rapid rectangular grid high-density electrical method three-dimensional inversion technology, the application effect is limited due to instrument precision and construction conditions, the data body can meet the three-dimensional inversion requirement through conversion and arrangement, and the inversion result is more accurate under the condition of encrypted profile spacing and polar distance.
Drawings
In order to more clearly illustrate the technical solutions in the embodiments of the present invention, the drawings needed to be used in the description of the embodiments are briefly introduced below, it is obvious that the drawings in the following description are only some embodiments of the present invention, and other drawings can be obtained by those skilled in the art without inventive efforts, wherein:
FIG. 1 is a 2D high density electrical profile of the present invention;
FIG. 2 is a three-dimensional electrical depth measurement non-uniform grid diagram of the present invention;
FIG. 3 is a cross-sectional view of a plot of apparent resistivity values in the x and y directions according to the present invention;
FIG. 4 is a first three-dimensional inversion result diagram of the present invention;
FIG. 5 is a second diagram of the three-dimensional inversion result of the present invention;
FIG. 6 is a comprehensive three-dimensional view of a geological geophysical survey of the present invention;
FIG. 7 is a three-dimensional inversion slice diagram of the invention at a depth of 40m, 80m, 120 m.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
The three-dimensional resistivity sounding data format of the high-density electrical method Wener device is similar to an uneven network measuring mode, the requirements of high-density electrical method Wener or equidistant quadrupole arrangement are not required to be followed, the uneven network measuring data format can be formed for three-dimensional inversion as long as the relative position or the coordinate of the electrode is determined, and electrical depth profiles can be arranged at equal distance or unequal distance. The electrical depth profile can not be too sparse, so that the boundary error of the three-dimensional inversion geologic body can be increased.
Referring to fig. 1-7, in an embodiment of the present invention, a method for applying a three-dimensional resistivity depth measurement based on a high-density electrical method wenna device includes the following steps:
1) intercepting a data volume along a rectangular section at equal intervals by adopting a 2D high-density resistivity instrument;
2) constructing a three-dimensional uneven grid electrode arrangement method and a grid on the basis of the existing 2D high-density resistivity profile data;
3) and (4) converting two-dimensional data formats by combining the requirements of a three-dimensional inversion program, and compiling a 2D high-density electrical method non-uniform network measurement three-dimensional inversion data format.
The specific implementation process comprises the following steps:
1) in the gory large willow gold mine area, 6 high-density electrical method sections are arranged, the sections are numbered as 11, 12, 13, 14, 15 and 16, the section spacing is 50m, the section point distance is 10m, the section orientation is 295 degrees, observation data points of the inverted trapezoid of the 2D high-density electrical method section wenna device are intercepted and adjusted to be rectangular data points, as shown in a 2D high-density electrical method section layout diagram shown in fig. 1, the starting and ending point numbers are 38 and 62 respectively, the section length is 240m, the minimum AB distance is 60m, and the maximum AB distance is 660 m.
2) In the course of measurement related to resistivity imaging, the three-dimensional electrode arrangement grid comprises monopole-monopole measurement in different directions, 2D parallel line measurement, broken line grid measurement and different direction line measurement; the test is defined as an uneven measurement network electrode arrangement grid, as shown in a three-dimensional electrical depth measurement uneven grid diagram of fig. 2:
the uneven measuring net grid is arranged based on symmetric quadrupole depth measurement equal ratio arrangement AB and MN; taking the goaf detection in the great salix matsudana area as an example, the data format of the high-density electrical method wenna arrangement is converted into the relationship between the polar distances of the symmetric quadrupole electrical depths AB and MN as follows:
converting the high-density electrical method Wener arrangement data format into a symmetric quadrupole electrical depth AB and MN polar distance relationship as follows:
AB/2:30m,60m,90m,120m,150m,180m,210m,240m,270m,300m,330m
MN/2:10m,20m,30m,40m,50m,60m,70m,80m,90m,100m,110m
the three-dimensional electrical sounding measures the uneven grid spacing, each section electrical sounding point corresponds to A, B, M, N electrode positions, and the arrangement is as follows: taking the number 38 of the initial sounding point of the No. 16 section as the origin of coordinates 0 and 0, taking the position of an A, B, M, N electrode corresponding to each sounding point as the abscissa, taking the equal intervals of 0m, 50m, 100m, 150m, 200m and 250m of the No. 11, 12, 13, 14, 15 and 16 sections as the ordinate to establish the corresponding relation between all electrode points and resistivity and elevation, wherein the grid interval of the abscissa is as follows:
Figure GDA0002978120390000081
the grid spacing is symmetrically arranged as follows: 0m 20m 10m 10m … … 10m 20m
The vertical coordinate grid distance is: 50 m.
3) The three-dimensional inversion data format for compiling the 2D high-density electrical method non-uniform measurement network is as follows: the non-uniform grid electrode arrangement and grid spacing are combined with the three-dimensional inversion program requirement, non-uniform grid three-dimensional inversion data format conversion software is compiled, and the high-density Wener data format is converted into the non-uniform grid three-dimensional inversion data format as follows:
1 Block
89 6
Nonuniform grid
x-location of grid-lines
0 20 30 40……870 880 900
y-location of grid-lines
0 50 100 150 200 250
7
858
300 0 360 0 320 0 340 0 251
270 0 390 0 310 0 350 0 327
240 0 420 0 300 0 360 0 392
……………………………………………………
Topography
2
244 244 243………………………………221
0
0
0
0
0。
(4)2D high-density electrical method non-uniform net-measuring three-dimensional inversion scheme column
The medium generation quartz vein distributed in the large willow gold mine area has high resistance characteristic, the resistivity is more than 9000 ohm m, the granite is in medium resistance characteristic, the resistivity is between 3000-4000 ohm m, the fractured rock is in low resistance characteristic, the resistivity is less than 1000 ohm m, theoretically, the resistivity of dry rock and air is relatively large, the resistivity of the rock is rapidly reduced due to the moisture of pores and cracks and along with the increase of the humidity or the saturation, and meanwhile, the resistivity of different rocks has certain difference due to the difference of the mineralization degrees of moisture under the condition that the moisture content is the same. Therefore, the resistivity of a fault depends mainly on the fracture degree of the fault and the saturation degree of water; the resistivity of the formation depends primarily on the saturation and water mineralization of the rock. In general, when faults, formation fractures and goafs are hydrated, their resistivity is much less than that of water-free surrounding rock; when faults, formation fractures and goafs are free of water, their resistivity is generally characterized by high resistivity. In addition, the lithology of the rock stratum can change along with the change of the depth, the corresponding resistivity is greatly changed along the longitudinal direction (depth), the overall resistivity trend is increased along with the increase of the depth, and the change of the resistivity of the bedrock top interface weathering zone is reflected.
Through on-site investigation, 6 high-density electrical method sections of No. 11, No. 12, No. 13, No. 14, No. 15 and No. 16 are distributed in a section where a goaf possibly exists, the section spacing is 50m, the point spacing is 20m, after format conversion is carried out on data, Swedish RES3D software is used for carrying out three-dimensional inversion, low-resistance abnormity exists below the middle ground surface of the 6 sections, and the data are deduced to be caused by the goaf, namely: the three-dimensional inversion result graph of the X-direction and the Y-direction is shown in a cross section drawing of an apparent resistivity value in the X-direction and the Y-direction in fig. 3, the three-dimensional inversion result graph I in fig. 4 and the three-dimensional inversion result graph II in fig. 5, the geological and geophysical exploration comprehensive three-dimensional stereogram in fig. 6 and the three-dimensional inversion slice graphs of the depths of 40, 80 and 120m in fig. 7. The transverse scales of the above graphs are labeled differently, the x direction in the x and y direction apparent resistivity value simulated section graph is drawn according to the grid distance of 0 … … 900m, the three-dimensional inversion result graph is drawn according to the actual length of the section of 0-240m, and the inversion abnormal relative positions are correspondingly consistent.
WZK14-1 verification drill holes are arranged in the middle of the No. 14 section (see figure 6), faults are seen at a depth of 39.5-43.6m, lithology is weak limonite mineralized and fragmented granite, and the core sampling rate is 89%; a goaf is formed at 120.3-122.3m, lithology mainly comprises chlorite petrochemical property and kaolin property, a rock core is broken and has strong alteration, and the average rock core sampling rate is 25%; a goaf at 140-141.8m, a lithology of granite, broken rock core, and a rock core sampling rate of 100%; and the area is a goaf at 217.6-219.7m, the lithology is granite, the core is broken and angular, and the core sampling rate is 19%.
(5) Precision error analysis
The No. 11, 12, 13, 14, 15 and 16 high-density electrical method sections in the area are subjected to resistivity two-dimensional inversion, the iteration times are 2 times as low as possible and 4 times as high as possible, the mean square error is 19.7% as low as possible and 60.4% as high as possible, the mean square error is 40.3%, the three-dimensional inversion iteration times are 6 times and 49.8% of the mean square error, the difference between the two is 9.5%, and the resistivity two-dimensional inversion method sections are shown in Table 1
TABLE 1 precision error statistics table
Figure GDA0002978120390000111
Due to the fact that the iteration times are different, inversion mean square errors are different, inversion instability can be caused by excessively high iteration times, and 2D high-density electrical method uneven network measurement three-dimensional inversion can be converged after 6 times of iteration. And drilling verification shows that the result of the three-dimensional inversion after 6 iterations reflects the spreading condition of the goaf.
The above description is only an embodiment of the present invention, and not intended to limit the scope of the present invention, and all modifications of equivalent structures and equivalent processes, which are made by the present specification, or directly or indirectly applied to other related technical fields, are included in the scope of the present invention.
It will be evident to those skilled in the art that the invention is not limited to the details of the foregoing illustrative embodiments, and that the present invention may be embodied in other specific forms without departing from the spirit or essential attributes thereof. The present embodiments are therefore to be considered in all respects as illustrative and not restrictive, the scope of the invention being indicated by the appended claims rather than by the foregoing description, and all changes which come within the meaning and range of equivalency of the claims are therefore intended to be embraced therein. Any reference sign in a claim should not be construed as limiting the claim concerned.
Furthermore, it should be understood that although the present description refers to embodiments, not every embodiment may contain only a single embodiment, and such description is for clarity only, and those skilled in the art should integrate the description, and the embodiments may be combined as appropriate to form other embodiments understood by those skilled in the art.

Claims (3)

1. A three-dimensional resistivity depth measurement application method based on a high-density electrical method Wener device is characterized by comprising the following steps: the specific method comprises the following steps:
1) intercepting a data volume along a rectangular section at equal intervals by adopting a 2D high-density resistivity instrument;
2) constructing a three-dimensional uneven grid electrode arrangement method and a grid on the basis of the existing 2D high-density resistivity profile data;
3) combining the requirements of a three-dimensional inversion program, performing two-dimensional data format conversion, and compiling a 2D high-density electrical method non-uniform network measurement three-dimensional inversion data format; the three-dimensional inversion data format for compiling the 2D high-density electrical method non-uniform measurement network is as follows: the non-uniform grid electrode arrangement and grid spacing are combined with the three-dimensional inversion program requirement, non-uniform grid three-dimensional inversion data format conversion software is compiled, and the high-density Wener data format is converted into the non-uniform grid three-dimensional inversion data format as follows:
Figure FDA0002978120380000011
Figure FDA0002978120380000021
2. the three-dimensional resistivity depth measurement application method based on the high-density electrical method Wener device according to claim 1, characterized in that: the specific interception method comprises the following steps: 6 high-density electrical method sections of 2D are arranged in the selected area, the sections are numbered as 11, 12, 13, 14, 15 and 16, the section spacing is 50m, the section point spacing is 10m, the section orientation is 295 degrees, the observation data points of the inverted trapezoid of the 2D high-density electrical method section wenner device are intercepted and adjusted to be rectangular data points, the starting point and the ending point are 38 and 62 respectively, the section length is 240m, the minimum AB distance is 60m, and the maximum AB distance is 660 m.
3. The three-dimensional resistivity depth measurement application method based on the high-density electrical method Wener device according to claim 1, characterized in that: the method for constructing the three-dimensional uneven grid electrode arrangement and the specific method of the grid are as follows:
converting the high-density electrical method Wener arrangement data format into a symmetric quadrupole electrical depth AB and MN polar distance relationship as follows:
AB/2:30m,60m,90m,120m,150m,180m,210m,240m,270m,300m,330m
MN/2:10m,20m,30m,40m,50m,60m,70m,80m,90m,100m,110m
the three-dimensional electrical sounding measures the uneven grid spacing, each section electrical sounding point corresponds to A, B, M, N electrode positions, and the arrangement is as follows: taking a No. 16 section No. 38 sounding point as a coordinate origin (0,0) of a three-dimensional electrical sounding grid, taking A, B, M, N electrode positions corresponding to each sounding point as abscissa, taking equal intervals of 0m, 50m, 100m, 150m, 200m and 250m of No. 11, 12, 13, 14, 15 and 16 sections as ordinate to establish corresponding relations between all electrode points and resistivity and elevation, wherein the intervals of the abscissa grid are as follows:
Figure FDA0002978120380000031
the grid spacing is symmetrically arranged as follows: 0m 20m 10m 10m … … 10m 20m
The vertical coordinate grid distance is: 50 m.
CN201910784248.4A 2019-08-23 2019-08-23 Three-dimensional resistivity depth sounding application method based on high-density electrical method Wener device Active CN110471122B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910784248.4A CN110471122B (en) 2019-08-23 2019-08-23 Three-dimensional resistivity depth sounding application method based on high-density electrical method Wener device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910784248.4A CN110471122B (en) 2019-08-23 2019-08-23 Three-dimensional resistivity depth sounding application method based on high-density electrical method Wener device

Publications (2)

Publication Number Publication Date
CN110471122A CN110471122A (en) 2019-11-19
CN110471122B true CN110471122B (en) 2021-07-30

Family

ID=68512074

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910784248.4A Active CN110471122B (en) 2019-08-23 2019-08-23 Three-dimensional resistivity depth sounding application method based on high-density electrical method Wener device

Country Status (1)

Country Link
CN (1) CN110471122B (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111427092B (en) * 2020-05-13 2022-12-27 天津华北地质勘查总院 Method for improving geophysical prospecting electrical sounding deep field construction efficiency
CN111856589B (en) * 2020-06-18 2021-07-06 山东大学 Three-dimensional imaging method and system for ground surface comprehensive geophysical prospecting
CN113125517B (en) * 2021-04-09 2023-11-28 浙江省水利河口研究院(浙江省海洋规划设计研究院) Device and method for monitoring root infiltrating irrigation by utilizing three-dimensional resistivity

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008128994A (en) * 2006-11-27 2008-06-05 Shinku Jikkenshitsu:Kk Ionization vacuum gage, and mass spectrometer using the same
CN102587984A (en) * 2012-02-20 2012-07-18 太原理工大学 Distribution type water information monitoring system for underground working face of coal mine
CN106405250A (en) * 2016-08-31 2017-02-15 山东电力工程咨询院有限公司 High-density earth resistivity measurement system suitable for complex terrain conditions and high-density earth resistivity measurement method thereof
CN109375271A (en) * 2018-09-29 2019-02-22 山东大学 A kind of the multi -components resistivity detection method and system of any cloth pole
CN110045436A (en) * 2019-03-31 2019-07-23 云南省环境科学研究院(中国昆明高原湖泊国际研究中心) A kind of accurate investigation method in solid waste stockpiling place

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9759829B2 (en) * 2015-02-28 2017-09-12 Glenn Dale Chubak System and method for geophysical data collection

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008128994A (en) * 2006-11-27 2008-06-05 Shinku Jikkenshitsu:Kk Ionization vacuum gage, and mass spectrometer using the same
CN102587984A (en) * 2012-02-20 2012-07-18 太原理工大学 Distribution type water information monitoring system for underground working face of coal mine
CN106405250A (en) * 2016-08-31 2017-02-15 山东电力工程咨询院有限公司 High-density earth resistivity measurement system suitable for complex terrain conditions and high-density earth resistivity measurement method thereof
CN109375271A (en) * 2018-09-29 2019-02-22 山东大学 A kind of the multi -components resistivity detection method and system of any cloth pole
CN110045436A (en) * 2019-03-31 2019-07-23 云南省环境科学研究院(中国昆明高原湖泊国际研究中心) A kind of accurate investigation method in solid waste stockpiling place

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
RES2DINV在江西清江盆地常规电测深数据反演中的应用;刘迁 等;《地球物理学进展》;20181231;第33卷(第6期);第2416-2427页 *

Also Published As

Publication number Publication date
CN110471122A (en) 2019-11-19

Similar Documents

Publication Publication Date Title
CN110471122B (en) Three-dimensional resistivity depth sounding application method based on high-density electrical method Wener device
CN101930083B (en) United inversion method for multiple inter-well electromagnetic wave tomography hole pairs
CN107742015A (en) DC lasering electric method Three-dimensional Numerical Simulation Method based on any dipole dipole device
CN102866417A (en) Device and method for seismic cross hole computed tomography (CT) detection and tomography of underground cave
CN109271738A (en) It is a kind of for obtaining the Method for Numerical Inversion of roadway surrounding rock Weibull distribution parameter
CN102508303A (en) Advanced detection method for focusing chromatography induced polarization of underground engineering
CN108693572A (en) A kind of ground stress field inversion method based on three-dimensional modeling
CN107144380A (en) The dynamic testing method of country rock disturbance stress field during constructing tunnel
CN105604066A (en) Application of resistivity profiling method to detection of building foundation pit support structure water leakage
CN106706715A (en) Polluted soil detection method based on three-dimensional high-density electrical resistivity method
CN111830578A (en) Electromagnetic wave CT method for accurately detecting underground pipeline
Wang et al. Determination of discontinuity persistent ratio by Monte-Carlo simulation and dynamic programming
Jin et al. Elliptical fracture network modeling with validation in Datong Mine, China
Lin et al. Rock mass characterization for shallow granite by integrating rock core indices and seismic velocity
CN108614307B (en) Inverted arch quality detection method and system by using comprehensive geophysical prospecting technology
Wang et al. Research on multi-frequency ultrasonic scanning detecting technology of cavity in the test borehole
CN113779811A (en) Quantitative analysis method for structural deformation to current ground stress orientation disturbance
Lu et al. Tunnel concealed karst cave joint detection by tunnel seismic and transient electromagnetic
CN103335169A (en) Geotechnical engineering investigation method
Bandpey et al. Validation of 3D discrete fracture network model focusing on areal sampling methods-a case study on the powerhouse cavern of Rudbar Lorestan pumped storage power plant, Iran
Zhifa et al. Three-dimensional back-analysis of displacements in exploration adits—principles and application
Ni et al. Simulation of borehole radar responses to rough fractures based on 3-D conformal FDTD
CN110488366B (en) Three-dimensional resistivity sounding application method based on non-uniform measuring network
CN212364603U (en) Tunnel unfavorable geologic body advanced geological forecast detection simulation device
Esmaeilzadeh et al. Shape effect of fractures on intensity and density of discreet fracture networks

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