CN105549101B - A kind of transient electromagnetic data differential conductance means of interpretation - Google Patents

A kind of transient electromagnetic data differential conductance means of interpretation Download PDF

Info

Publication number
CN105549101B
CN105549101B CN201610061215.3A CN201610061215A CN105549101B CN 105549101 B CN105549101 B CN 105549101B CN 201610061215 A CN201610061215 A CN 201610061215A CN 105549101 B CN105549101 B CN 105549101B
Authority
CN
China
Prior art keywords
transient electromagnetic
differential conductance
data
curve
electromotive force
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
CN201610061215.3A
Other languages
Chinese (zh)
Other versions
CN105549101A (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.)
China University of Mining and Technology CUMT
Original Assignee
China University of Mining and Technology CUMT
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 China University of Mining and Technology CUMT filed Critical China University of Mining and Technology CUMT
Priority to CN201610061215.3A priority Critical patent/CN105549101B/en
Publication of CN105549101A publication Critical patent/CN105549101A/en
Application granted granted Critical
Publication of CN105549101B publication Critical patent/CN105549101B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/30Assessment of water resources

Abstract

A kind of transient electromagnetic data differential conductance means of interpretation, belongs to geophysical exploration data interpretation method.Differential conductance means of interpretation is the induced electromotive force based on the acquisition of transient electromagnetic prospecting method in the different means of interpretation of the different regional decay characteristics of resistivity;S1, Different Strata model designed according to electrical property difference, forward simulation is carried out to various stratigraphic models, and data acquired collected;S2, the transient electromagnetic induced electromotive force data for acquisition, are explained using differential conductance method to transient electromagnetic data, the fluctuating situation for the curve that observation transient electromagnetic differential conductance value changes over time.Used when judging suitable for ground and underground strata division, water zone.It is different that electromagnetic field decay rate is spread in stratum based on different resistivity, line translation is entered to the induced electromotive force V in Transient Electromagnetic Apparatus data measured and time t, divide electric layers by the analysis to transformation results and judge water zone, enrich transient electromagnetic data means of interpretation.

Description

A kind of transient electromagnetic data differential conductance means of interpretation
Technical field
The present invention relates to a kind of geophysical exploration data interpretation method, particularly a kind of transient electromagnetic data differential conductance Means of interpretation.
Background technology
China is in fast development, and resource is the basic of guarantee national development, in resource exploration and recovery process In, geophysical exploration serves significant role.Particularly in terms of mine water inrush preventing and treating, transient electromagnetic method has obtained widely Using reducing the generation of mine water inrush accident, reduce economic loss, ensure that the safety of the people.Transient electromagnetic method is one The method of the great development prospect of kind, it is sensitive to low-resistance abnormal response, can find out aqueous geology such as Cave and passage, colliery Goaf, the irregular water body in deep etc..Transient electromagnetic method has that investigation depth is big, sensitive to low-resistance geologic body, with detecting target Couple, the advantages that exception response is strong, and form is simple, and resolution capability is strong.
Transient electromagnetic method is also referred to as TDEM, and abbreviation TEM, it is using earth-free loop line or is grounded line source to underground Launch pulsatile once magnetic field, in pulsatile once magnetic field tempus intercalare, utilize the side of coil or grounding electrode observation secondary turbulence field Method.Briefly, the general principle of transient electromagnetic method is exactly the law of electromagnetic induction.Attenuation process is generally divided into early, neutralization late period. The electromagnetic field of early stage is equivalent to the radio-frequency component in frequency domain, and decay is fast, and skin depth is small;And late period composition is then equivalent to frequency Low-frequency component in domain, decay is slow, and skin depth is big.The secondary field of each period changes over time rule after being powered off by measurement Rule, can obtain the ground electrical feature of different depth.
During transient electromagnetic detecting, electromagnetic field is decayed very fast in high resistant object, and is decayed in low-resistance object It is relatively slow, difference between the two is more obvious.
But for transient electromagnetic data common version be usually using apparent resistivity section diagram, multiple instrument channels profile figure as Main, means of interpretation is single.
The content of the invention
The invention aims to provide a kind of transient electromagnetic data differential conductance means of interpretation, solves current transient electromagnetic Data interpretation method relative to other geophysical exploration methods it is relatively simple the problem of.
The object of the present invention is achieved like this:Differential conductance means of interpretation is obtained based on transient electromagnetic prospecting method Induced electromotive force is in the different means of interpretation of the different regional decay characteristics of resistivity;Comprise the following steps that:
S1, Different Strata model designed according to electrical property difference, forward simulation is carried out to various stratigraphic models, and to obtained number According to being collected.
In the S1 steps, to the different types of stratigraphic model of design in the case where keeping other factors constant, become Change its certain Graph One factor and carry out forward simulation, and resulting transient electromagnetic data is subjected to Classifying Sum, specific steps include:
(1) homogeneous half space transient electromagnetic model, resistivity ρ, depth h, and h=∞ are designed, changes the big of ρ values It is small, obtain multigroup homogeneous half space transient electromagnetic data, including period, time t and corresponding induced electromotive force value V;
(2) two layers of stratum transient electromagnetic model of design level, the first layer resistivity are ρ1, the first layer depth is h1;The second layer Resistivity is ρ2, the second layer depth is h2=∞, changes ρ respectively1、h1、ρ2Size, obtain two layers of transient electromagnetic data of level, Including period, time t and corresponding induced electromotive force value V;
(3) three layers of stratum transient electromagnetic model of design level, the first layer resistivity are ρ1, the first layer depth is h1;The second layer Resistivity is ρ2, the second layer depth is h2;Third layer resistivity is ρ3, third layer depth is h3=∞, changes ρ respectively1、h1、ρ2、 h2、ρ3Size, obtain three layers of transient electromagnetic data of level, including period, time t and corresponding induced electromotive force value V.
The fluctuating situation for the curve that S2, observation transient electromagnetic differential conductance value change over time, obtains transient electromagnetic sensing Electromotive force data, transient electromagnetic data is explained using differential conductance method.
In the S2 steps, differential conductance method, specific steps include:
S21, the cubic spline interpolation on time t is done to gained transient electromagnetic induced electromotive force value V, obtained more abundant Transient electromagnetic data, ensure by differential conductance processing after curve obtained it is round and smooth as far as possible, the solution more conducively to data Release;
S22, under log-log coordinate system, draw out the curve that transient electromagnetic induced electromotive force V changes with time t;
S23, the committed step of differential conductance processing is carried out to transient electromagnetic data to ask for transient electromagnetic differential conductance value K, expression formula are:
Its differential conductance value K is asked with the time t curves changed to the induced electromotive force V under log-log coordinate system;
S24, the curve that differential conductance value K changes with time t, i.e. transient electromagnetic differential are drawn out under single pair number coordinate system Conductance plots;
S25, same type of differential conductance value K is pooled to together with the curve that time t changes according to its changed factor In one coordinate system, it is fabricated to and collects curve map;
When to three layers of H type Stratigraphic Curves, only changing its second layer depth h2, keep other specification it is constant, then its change because Son is the second layer depth h2, the obtained curve map that collects is properly termed as three layers of H types stratum second layer change in depth differential conductance and converges Master curve figure;
S26, according to transient electromagnetic differential conductance curve fluctuating situation, observing its extreme point position accurately can divide ground Underground low-resistance anomalous body is found in bed boundary.
Traditional transient electromagnetic data is generally based on the size of apparent resistivity, is provided by apparent resistivity isoline figure Material explains that Comparatively speaking, advantages and advantages of the invention are:Explanation to transient electromagnetic data obtains in terms of layering It is progressive, the ability to bed boundary division is improved, can effectively judge the interface of anomalous body and surrounding country rock;Based on transition Electromagnetic induction electromotive force is not done any near to initial data the different regional decay characteristics of resistivity are widely different the characteristics of Seemingly, the authenticity and integrality of data are ensure that, enriches transient electromagnetic means of interpretation.
Brief description of the drawings:
Fig. 1 is the curve that three layers of H types formation induction electromotive force V of transient electromagnetic change with time t under log-log coordinate system.
Fig. 2 is the curve that three layers of H types formation induction electromotive force V of transient electromagnetic after coordinate transform change with time t.
Fig. 3 is the curve that three layers of H types stratum differential conductance value K of transient electromagnetic change with time t under single pair number coordinate system.
Fig. 4 a are three layers of H types stratum second layer change in depth induced electromotive force congruence of curves of transient electromagnetic under single pair number coordinate system Total figure.
Fig. 4 b are that three layers of H types stratum second layer change in depth differential conductance curve of transient electromagnetic collect under single pair number coordinate system Figure.
Fig. 5 a be single pair number coordinate system under transient electromagnetic with half space change in resistance induced electromotive force curve summary view.
Fig. 5 b be single pair number coordinate system under transient electromagnetic with half space change in resistance differential conductance curve summary view.
Fig. 6 a are two layers of D types stratum first layer change in depth induced electromotive force congruence of curves of transient electromagnetic under single pair number coordinate system Total figure.
Fig. 6 b are that two layers of D types stratum first layer change in depth differential conductance curve of transient electromagnetic collects under single pair number coordinate system Figure.
Fig. 7 a are three layers of Q types stratum first layer change in depth induced electromotive force congruence of curves of transient electromagnetic under single pair number coordinate system Total figure.
Fig. 7 b are that three layers of Q types stratum first layer change in depth differential conductance curve of transient electromagnetic collect under single pair number coordinate system Figure.
Fig. 8 a are three layers of Q types stratum second layer change in depth induced electromotive force congruence of curves of transient electromagnetic under single pair number coordinate system Total figure.
Fig. 8 b are that three layers of Q types stratum second layer change in depth differential conductance curve of transient electromagnetic collect under single pair number coordinate system Figure.
Fig. 9 is transient electromagnetic data differential conductance means of interpretation flow chart.
Embodiment
Embodiments of the invention are further described below in conjunction with the accompanying drawings:
Differential conductance means of interpretation is that the induced electromotive force based on the acquisition of transient electromagnetic prospecting method is different in resistivity The different means of interpretation of regional decay characteristics;Comprise the following steps that:
S1, Different Strata model designed according to electrical property difference, forward simulation is carried out to various stratigraphic models, and to obtained number According to being collected.
In the S1 steps, to the different types of stratigraphic model of design in the case where keeping other factors constant, become Change its certain Graph One factor and carry out forward simulation, and resulting transient electromagnetic data is subjected to Classifying Sum, specific steps include:
(1) homogeneous half space transient electromagnetic model, resistivity ρ, depth h, and h=∞ are designed, changes the big of ρ values It is small, obtain multigroup homogeneous half space transient electromagnetic data, including period, time t and corresponding induced electromotive force value V;
(2) two layers of stratum transient electromagnetic model of design level, the first layer resistivity are ρ1, the first layer depth is h1;The second layer Resistivity is ρ2, the second layer depth is h2=∞, changes ρ respectively1、h1、ρ2Size, obtain two layers of transient electromagnetic data of level, Including period, time t and corresponding induced electromotive force value V;
(3) three layers of stratum transient electromagnetic model of design level, the first layer resistivity are ρ1, the first layer depth is h1;The second layer Resistivity is ρ2, the second layer depth is h2;Third layer resistivity is ρ3, third layer depth is h3=∞, changes ρ respectively1、h1、ρ2、 h2、ρ3Size, obtain three layers of transient electromagnetic data of level, including period, time t and corresponding induced electromotive force value V.
The fluctuating situation for the curve that S2, observation transient electromagnetic differential conductance value change over time, obtains transient electromagnetic sensing Electromotive force data, transient electromagnetic data is explained using differential conductance method.
In the S2 steps, differential conductance method, specific steps include:
S21, the cubic spline interpolation on time t is done to gained transient electromagnetic induced electromotive force value V, obtained more abundant Transient electromagnetic data, ensure by differential conductance processing after curve obtained it is round and smooth as far as possible, the solution more conducively to data Release;
S22, under log-log coordinate system, draw out the curve that transient electromagnetic induced electromotive force V changes with time t;
S23, according to differential conductance expression formula
Its differential conductance value K is asked with the time t curves changed to the induced electromotive force V under log-log coordinate system;
S24, the curve that differential conductance value K changes with time t, i.e. transient electromagnetic differential are drawn out under single pair number coordinate system Conductance plots;
S25, same type of differential conductance value K is pooled to together with the curve that time t changes according to its changed factor In one coordinate system, it is fabricated to and collects curve map;
When to three layers of H type Stratigraphic Curves, only changing its second layer depth h2, keep other specification it is constant, then its change because Son is the second layer depth h2, the obtained curve map that collects is properly termed as three layers of H types stratum second layer change in depth differential conductance and converges Master curve figure;
S26, according to transient electromagnetic differential conductance curve fluctuating situation, observing its extreme point position accurately can divide ground Underground low-resistance anomalous body is found in bed boundary.
Embodiment 1:
The first step, numerical simulation, transient electromagnetic data needed for acquisition are carried out by transient electromagnetic forward simulation software;
Second step, the transient electromagnetic data obtained in the first step is arranged, for every group of data, the first row of data Time t is classified as period, second, the 3rd is classified as induced electromotive force value V;
3rd step, change curves of the induced electromotive force value V on time t is drawn in log-log coordinate system;
4th step, to the curve obtained in the 3rd step, handled in differential conductance method, specific method is as follows:
Keep the curve shapes that change over time of induced electromotive force value V in Fig. 1 constant, curve respective coordinates value is taken pair Number, then obtain the arithmetic coordinate of decile, such as referred to as coordinate transform, Fig. 2;
Its gradient is asked to the curve in Fig. 2 again, then its Grad is exactly required differential conductance value K.
5th step, in single pair number coordinate system (x-axis is logarithmic coordinates), draw out transient electromagnetic differential conductance value K at any time Between the curve that changes, such as Fig. 3.
6th step, the curve that same type of induced electromotive force value V, differential conductance value K are changed with time t is according to it Changed factor is pooled in the same coordinate system, is fabricated to and is collected curve map, such as Fig. 4 a, Fig. 4 b:To three layers of H type Stratigraphic Curves, only Change changes its first layer depth h2, holding other specification is constant, then its changed factor is the second layer depth h2, obtain collecting song Line chart, Fig. 4 a are properly termed as three layers of H types stratum second layer change in depth induced electromotive force and collect curve map, and Fig. 4 b are properly termed as three Layer H types stratum second layer change in depth differential conductance collects curve map.
7th step, qualitative division is carried out to stratum according to the extreme point of curve, as two extreme values be present in the curve in Fig. 3, Maximum is the line of demarcation of first layer and the second layer, and minimum is the line of demarcation of the second layer and third layer.
It is of the invention entirely different with existing transient electromagnetic data means of interpretation, with transient electromagnetic signal in resistivity differently Area's attenuation degree is not all theoretical foundation, it is proposed that transient electromagnetic data differential conductance means of interpretation.
It is described above, it is only the embodiment of the present invention more for simplicity, protection scope of the present invention is not limited to This, any one skilled in the art can become apparent within the technical scope of present disclosure The simple change or equivalence replacement of this technology are each fallen within protection scope of the present invention.

Claims (1)

1. a kind of transient electromagnetic data differential conductance means of interpretation, it is characterized in that:Differential conductance means of interpretation is to be based on transient electrical The induced electromotive force that magnetic prospecting method obtains is in the different means of interpretation of the different regional decay characteristics of resistivity;Specific steps are such as Under:
S1, Different Strata model designed according to electrical property difference, forward simulation is carried out to various stratigraphic models, and data acquired entered Row collects;
The fluctuating situation for the curve that S2, observation transient electromagnetic differential conductance value change over time, obtains transient electromagnetic induction electric Gesture data, transient electromagnetic data is explained using differential conductance method;
, need to be to the different types of stratigraphic model of design in the case where keeping other factors constant in the S1 steps, change Its certain Graph One factor carries out forward simulation, and resulting transient electromagnetic data is carried out into Classifying Sum, and specific steps include:
(1)Homogeneous half space transient electromagnetic model, resistivity ρ, depth h, and h=∞ are designed, changes the size of ρ values, is obtained Multigroup homogeneous half space transient electromagnetic data, including period, time t and corresponding induced electromotive force value V;
(2)Two layers of stratum transient electromagnetic model of design level, the first layer resistivity are ρ1, the first layer depth is h1;Second layer resistance Rate is ρ2, the second layer depth is h2=∞, changes ρ respectively1、h1、ρ2Size, obtain two layers of transient electromagnetic data of level, including point Number, time t and corresponding induced electromotive force value V;
(3)Three layers of stratum transient electromagnetic model of design level, the first layer resistivity are ρ1, the first layer depth is h1;Second layer resistance Rate is ρ2, the second layer depth is h2;Third layer resistivity is ρ3, third layer depth is h3=∞, changes ρ respectively1、h1、ρ2、h2、ρ3 Size, obtain three layers of transient electromagnetic data of level, including period, time t and corresponding induced electromotive force value V;
In the S2 steps, differential conductance method, specific steps include:
S21, the cubic spline interpolation on time t is done to gained transient electromagnetic induced electromotive force value V, obtain more abundant wink Become electromagnetic data, ensure that curve obtained is round and smooth as far as possible after differential conductance processing, the explanation more conducively to data;
S22, under log-log coordinate system, draw out the curve that transient electromagnetic induced electromotive force V changes with time t;
S23, the committed step of differential conductance processing is carried out to transient electromagnetic data to ask for transient electromagnetic differential conductance value K, table It is up to formula:
Its differential conductance value K is asked with the time t curves changed to the induced electromotive force V under log-log coordinate system;
S24, the curve that differential conductance value K changes with time t, i.e. transient electromagnetic differential conductance are drawn out under single pair number coordinate system Curve;
S25, same type of differential conductance value K is pooled to same seat with the curve that time t changes according to its changed factor In mark system, it is fabricated to and collects curve map;
When to three layers of H type Stratigraphic Curves, only changing its second layer depth h2, keep other specification constant, then its changed factor is the Two layer depth h2, the obtained curve map that collects is properly termed as three layers of H types stratum second layer change in depth differential conductance and collects curve Figure;
S26, according to transient electromagnetic differential conductance curve fluctuating situation, ground stratum boundary can accurately be divided by observing its extreme point position Underground low-resistance anomalous body is found in face.
CN201610061215.3A 2016-01-28 2016-01-28 A kind of transient electromagnetic data differential conductance means of interpretation Expired - Fee Related CN105549101B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610061215.3A CN105549101B (en) 2016-01-28 2016-01-28 A kind of transient electromagnetic data differential conductance means of interpretation

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610061215.3A CN105549101B (en) 2016-01-28 2016-01-28 A kind of transient electromagnetic data differential conductance means of interpretation

Publications (2)

Publication Number Publication Date
CN105549101A CN105549101A (en) 2016-05-04
CN105549101B true CN105549101B (en) 2018-01-05

Family

ID=55828400

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610061215.3A Expired - Fee Related CN105549101B (en) 2016-01-28 2016-01-28 A kind of transient electromagnetic data differential conductance means of interpretation

Country Status (1)

Country Link
CN (1) CN105549101B (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106405665B (en) * 2016-11-18 2018-09-28 厦门大学 Transient electromagnetic conductivity inversion method based on DBIM
CN106772644B (en) * 2016-12-14 2018-06-26 中国矿业大学 mine transient electromagnetic three-component detection method
CN115542408B (en) * 2022-12-05 2023-03-28 成都理工大学 Ocean transient electromagnetic data preprocessing and rapid regional imaging method

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2352041A1 (en) * 2008-10-23 2011-08-03 Obshchestvo S Ogranichennoy Otvetstvennostyu "Sibirskaya Geofizicheskaya Nauchno-Proizvodstvennaya Kompaniya" Method for quantitatively separating the effects of electromagnetic induction and induced polarization
CN104360401A (en) * 2014-12-10 2015-02-18 中国科学院地质与地球物理研究所 Method for determining geological information of underground target body by transient electromagnetic field B
CN104749647A (en) * 2015-04-09 2015-07-01 中国科学院地质与地球物理研究所 Geological interface structure detection method, device and system
CN105277994A (en) * 2015-07-06 2016-01-27 中煤科工集团西安研究院有限公司 Determination method and device of differential conductance data synthetic aperture

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2352041A1 (en) * 2008-10-23 2011-08-03 Obshchestvo S Ogranichennoy Otvetstvennostyu "Sibirskaya Geofizicheskaya Nauchno-Proizvodstvennaya Kompaniya" Method for quantitatively separating the effects of electromagnetic induction and induced polarization
CN104360401A (en) * 2014-12-10 2015-02-18 中国科学院地质与地球物理研究所 Method for determining geological information of underground target body by transient electromagnetic field B
CN104749647A (en) * 2015-04-09 2015-07-01 中国科学院地质与地球物理研究所 Geological interface structure detection method, device and system
CN105277994A (en) * 2015-07-06 2016-01-27 中煤科工集团西安研究院有限公司 Determination method and device of differential conductance data synthetic aperture

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
An automatic technique for presentation of coincident-loop,impulse-response,transient electromagnetic data;Richadrd s.smith 等;《GEOPHYSICS》;19941031;第59卷(第10期);第1542-1550页 *
瞬变电磁微分成像方法的关键技术研究;马索尼;《中国优秀硕士学位论文全文数据库(基础科学辑)》;20131215(第S2期);第7-8页第2.2.1节,第13页第5段,第15-20第2.3-2.4节,26-28页第2.4.2节 *
瞬变电磁测深的微分电导成像;李貅 等;《煤田地质与勘探》;20030630;第31卷(第3期);第59-61页 *
瞬变电磁视纵向电导微分成像应用研究;李超;《中国优秀硕士学位论文全文数据库(基础科学辑)》;20100215(第02期);第1-48页 *

Also Published As

Publication number Publication date
CN105549101A (en) 2016-05-04

Similar Documents

Publication Publication Date Title
CN104278991B (en) Saline Lake Facies hydrocarbon source rock organic carbon and the polynary well logging computational methods of hydrocarbon potential
CN103510947B (en) The method for building up of dam, beach sandstone microfacies recognition mode
CN102419455B (en) Interwell parallel resistivity CT (computed tomography) testing method
CN104865614B (en) Complicated reservoirs Fluid Identification Method based on variable matrix parameter
CN104863574B (en) A kind of Fluid Identification Method suitable for tight sandstone reservoir
CN105549101B (en) A kind of transient electromagnetic data differential conductance means of interpretation
CN103698811A (en) Carbonate rock structure ingredient well logging quantitative recognition method and purpose thereof
CN107065019A (en) Applied to road disaster and the 3 D electromagnetic imaging device and application method that collapse detection
CN105064986A (en) Method for building reservoir four-property relationship spectrum by using conventional well detection and logging information
CN105891896A (en) Feature information recognition and analysis method for underground mined area
West et al. The borehole transient electromagnetic response of a three-dimensional fracture zone in a conductive half-space
CN105842738B (en) A kind of buried target body quantitatively judges method and apparatus extremely
CN112861890A (en) Reservoir evaluation model construction method and reservoir identification method
CN105676299A (en) Grounded source short offset transient electromagnetic method apparent resistivity determination method and device
CN104122590A (en) Oil and gas detection method and system based on electromagnetic survey
CN104793268B (en) The blind depth measurement method and device of a kind of transient electromagnetic detecting
CN105929457B (en) A kind of frequency spectrum discerning method of interflow subsurface drainage multidate information
AU2017279838B1 (en) Method for classifying deep rock geofacies based on data mining
CN104749645B (en) A kind of method and apparatus of coal seam substrate interface undulation form detection in ash top difficult to understand
CN104991281B (en) A kind of detection method and device of coal seam buried depth
CN106383363B (en) A kind of geologic prospecting method based on two-dimentional mise-a-la-masse method
CN105700026A (en) An underground target body conductive index extraction method and an apparatus thereof
CN106483580B (en) A kind of quick, intelligent measuring method of coal and rock stress characteristics index
Gao et al. Synergetic system for water body detection in coal mine: a case study
Hadush Magneto Telluric Methods in Geothermal Exploration at Aluto-Langano, Ethiopia

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into 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: 20180105

Termination date: 20200128

CF01 Termination of patent right due to non-payment of annual fee