CN110927632A - Frequency domain horizontal X-direction magnetic field component observation and data processing method - Google Patents
Frequency domain horizontal X-direction magnetic field component observation and data processing method Download PDFInfo
- Publication number
- CN110927632A CN110927632A CN201910905194.2A CN201910905194A CN110927632A CN 110927632 A CN110927632 A CN 110927632A CN 201910905194 A CN201910905194 A CN 201910905194A CN 110927632 A CN110927632 A CN 110927632A
- Authority
- CN
- China
- Prior art keywords
- horizontal
- magnetic field
- field component
- frequency domain
- normalization function
- 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.)
- Withdrawn
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/02—Measuring direction or magnitude of magnetic fields or magnetic flux
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/08—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with magnetic or electric fields produced or modified by objects or geological structures or by detecting devices
- G01V3/10—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with magnetic or electric fields produced or modified by objects or geological structures or by detecting devices using induction coils
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geophysics (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention provides a method for measuring horizontal X-direction magnetic field components of a frequency domain and simultaneously calculating the apparent resistivity of a whole area, which comprises the following steps: (1) at the approximation
Measuring frequency spectrum data of a horizontal X-direction magnetic field component in a frequency domain in an exploration area in the direction; (2) converting the frequency spectrum data of the horizontal X-direction magnetic field component into an actual measurement normalization function; (3) and finding out the point where the theoretical normalization function is equal to the actually measured normalization function by using the numerical value, and then acquiring the apparent resistivity of the whole area. And obtaining resistivity distribution data of the target stratum at the specific underground depth through frequency spectrum data processing and inversion. The invention only collects the magnetic field component in the horizontal X direction, and has convenient instrument arrangement, high construction efficiency and low engineering cost; the resolving power of the geological target in the horizontal direction is enhanced; solves the problem of the common phenomenon under the condition of complex landformTechnical difficulties of the conventional electromagnetic method; the method overcomes the defect that the long delay signal-to-noise ratio is lower than the short delay signal-to-noise ratio in the conventional time domain method.
Description
Technical Field
The invention belongs to the field of electromagnetic exploration of geological and underground mineral exploration such as petroleum, natural gas and unconventional oil and gas reservoirs, and particularly relates to a frequency domain horizontal X-direction magnetic field component observation and data processing method.
Background
The ground electromagnetic method is an electromagnetic exploration method which adopts an artificial source to emit an excitation signal and observes primary and induced electric field and magnetic field total fields in an exploration area so as to research the electrical structure characteristics of underground strata. The existing electromagnetic method has the technical defects of 'space observation blind area' and difficulty in finding apparent resistivity, and therefore, a novel field measurement and distribution mode of the magnetic field component in the horizontal X direction of the frequency domain ground is provided as shown in figure 1. The physical quantity measured directly or indirectly at the observation point (point P) is the magnetic field strength in the horizontal X-direction (i.e. in the direction parallel to the artificial emission source AB), denoted hx (f) (in nano-nt), where f is the observation frequency (in Hz). Compared with the prior related background art, the technical advantages of the invention are analyzed and described in detail in the section of specific implementation methods.
The field arrangement method adopted by the invention is shown in figure 1, AB is a grounded transmitting electric dipole, sinusoidal current is fed into the ground through a transmitter, the amplitude of the current intensity is I (unit ampere), and the dipole moment of the transmitting electric dipole is PEI.ab (in ampere meters)),
is the receiving-transmitting distance (distance from the midpoint O to the receiving point P in meters), azimuth angle
Connecting line for transmitting electric dipole central point O and measuring point P andthe angle between the transmitting electric dipoles AB. Hx (f) is the horizontal X-direction magnetic field component parallel to the emitting electric dipole. Typically the dot pitch of the dots (dots in figure 1) is between 2 and 100 metres. A plurality of measuring points are arranged along a straight line to form a two-dimensional section; a plurality of measuring points are arranged in a certain plane area to form a three-dimensional exploration surface. The measuring points are arranged on
Azimuthal region, in general
The transmitting-receiving distance r is 4000 to 10000 meters.
Generally, the frequency spectrum of the horizontal X-direction magnetic field component is obtained on all the measurement points frequency point by frequency point within a finite time range, namely: hx (Xi, Yi, ω j), i ═ 1.., M; j ═ 1.., N. M is the total number of measuring points actually observed on the ground, and N is the total number of frequency points observed on each measuring point. ω j is the circle frequency of the j-th frequency point: omegaj＝2πfj. (Xi, Yi) is the plane coordinates of the ith measurement point.
In order to obtain the resistivity distribution data of the target stratum with the specific underground depth by processing and inverting the frequency spectrum data, the invention also provides a method for defining and calculating the apparent resistivity of the whole region based on the actually measured normalization function of the horizontal X-direction magnetic field component.
Disclosure of Invention
The invention provides a method for observing frequency domain horizontal X direction magnetic field component and processing data, which aims to solve the problems in the prior art and comprises the following steps:
(1) the field data is collected by adopting an artificial source method and the arrangement mode shown in figure 1. The measuring points can be arranged along a straight line or in a certain range of plane areas. The selected measurement area approximately falls within
In the direction of (1), usually take
Due to the fact that
Thereby absolutely avoiding
The data processing can obtain stronger stability, and simultaneously, the maximization of both the amplitude of the magnetic field intensity and the signal-to-noise ratio of the data is ensured. The point distance is between 2 meters and 100 meters. The transmitting-receiving distance r is 4000 to 10000 meters.
(2) An independent magnetic probe for the magnetic field component in the horizontal X direction is arranged on each measuring point (for example, an induction coil (magnetic bar) which can indirectly convert the magnetic field component after measuring the induced electromotive force, a superconducting quantum interferometer for directly measuring the magnetic field component, and the like).
(3) The frequency spectrum of the horizontal X-direction magnetic field component is obtained on all measuring points in a frequency point by frequency point within a limited and short time range (for example, a time period of 1 day to 1 week required for oil field reservoir dynamic monitoring).
(4) The invention provides the 'whole-area' apparent resistivity based on the horizontal X-direction magnetic field component (
Ohm-meters), which is implied in the following equation:
wherein
I and K are each a molar amount of
The first and second class of modified bessel functions.
due to strict guarantee of
Therefore, it is
The calculation of (a) is stable and non-singular.
(5) The invention adopts a method (such as 'dichotomy') for solving a nonlinear equation in a numerical method to process the actually measured normalization function of the horizontal X-direction magnetic field component, namely, the 'whole-area' apparent resistivity is calculated in a mode of finding out the point where the actually measured normalization function and the theoretical normalization function of the horizontal X-direction magnetic field component are equal
The invention is realized by the following technical scheme:
(1) according to the field surveying mode shown in fig. 1, acquiring frequency domain horizontal X-direction magnetic field recording data, namely frequency spectrum Hx (Xi, Yi, ω j), i is 1. j ═ 1.., N;
(2) calculating the actually measured normalization function of the horizontal X-direction magnetic field component through the conversion processing of the formula (3)
(3) For actually measured
Calculating the "total zone" apparent resistivity ρ using a suitable numerical method (e.g., "dichotomy") of solving a nonlinear equationHAx(ω)。
The method provided by the invention only collects the magnetic field component in the horizontal X direction, so that the instrument layout is convenient, the construction efficiency is improved, the engineering cost is reduced, the signal-to-noise ratio is optimized, and the data processing is stable and reliable.
Drawings
FIG. 1 is a schematic diagram of a field layout mode of magnetic field components in the horizontal X direction of a frequency domain;
FIG. 2 is a "full field" apparent resistivity curve corresponding to the horizontal X-direction magnetic field component calculated by the present invention (
This example N = 43);
FIG. 3 is a plan view of the resistivity distribution of the oil layer in the fireflood burned area obtained by the horizontal X-direction magnetic field component measurement technique in the embodiment of the present invention.
Detailed Description
The technical scheme of the invention is further explained by the specific implementation mode in combination with the attached drawings. It is to be understood that the specific embodiments described herein are merely illustrative of the invention and are not limiting of the invention. It should be further noted that, for the convenience of description, only some but not all of the relevant aspects of the present invention are shown in the drawings.
The embodiment provides an embodiment of a method for observing a magnetic field component in a horizontal X direction of a frequency domain and processing data, which comprises the following steps:
firstly, acquiring horizontal X-direction magnetic field recording data of a frequency domain acquired in the field, namely a frequency spectrum Hx (Xi, Yi, and ω j), wherein i is 1. j ═ 1.., N.
By adopting the arrangement mode shown in FIG. 1, an independent magnetic probe of the horizontal X-direction magnetic field component is arranged on each measuring point. In a limited and short time frame (e.g. 1 day to 1 week), on all stationsThe frequency spectrum of the horizontal X-direction magnetic field component is obtained from frequency point to frequency point. The selected measurement area approximately falls within
In the direction of (1), in general
The point distance is between 2 meters and 100 meters. The transmitting-receiving distance r is 4000 to 10000 meters.
Secondly, converting the frequency spectrum obtained by actual measurement into an actually measured normalization function of the horizontal X-direction magnetic field component
The method can be completed by adopting the formula (3).
Thirdly, actually measuring a normalization function of the horizontal X-direction magnetic field component, and calculating the 'whole-area' apparent resistivity of the corresponding stratum by adopting a numerical method for solving a nonlinear equation
The results of theoretical calculations and analysis show that: theoretical normalization function
The relationship with the variable a ═ real (kr) (a is a real part of kr which is a frequency domain parameter) is monotonous, so that the "dichotomy" in the numerical method can be selected. The specific calculation procedure is as follows.
(1) Taking the initial two end points of the resistivity value as (i.e. defining the root separation interval of the solution)
Ommi sum
Ohm-rice.
(2) Finding out the actually measured normalization function by the dichotomy
And a theoretical normalization function
Points of equal value (kr)s＝nI.e. satisfy the equality condition
Point (2) of (c).
(3) The apparent resistivity of the whole area is calculated by the following formula
(4) The specific calculation process of the 3 steps is repeated for j being 1, …, and N is circulated, so that the apparent resistivity of the whole region corresponding to the horizontal X-direction magnetic field component at all frequency points can be calculated
(5) By implementing the above 4-step specific calculation process on each measuring point (i is 1.., M cycles), the "full-area" apparent resistivity corresponding to the horizontal X-direction magnetic field component at all the measuring points can be calculated.
In order to verify the effectiveness of the method for calculating the apparent resistivity of the whole region, a theoretical one-dimensional laminar electrical model is set, and the parameters are as follows: 8 electrical layers, 1 st to 8 th layer resistivity 50, 10, 1000, 5, 200, 10, 100 and 10 ohm meters respectively, thickness 100, 500, 200 and 250 meters. The transmit-receive distance r is 8100 meters. The "full zone" apparent resistivity calculated for the theoretical model is shown in FIG. 2. The method is a curve with continuity, and keeps an approximate positive corresponding relation with the resistivity characteristics of the underground electric stratum (similar to an apparent resistivity curve of a Magnetotelluric (MT) method and finally has nothing to do with the transceiving distance and the azimuth angle). It has the main functions as follows: firstly, an initial model for inversion is convenient to construct; and secondly, the method can be used for constructing an objective function suitable for inversion, so that the inversion can be ensured to normally start and continue to be carried out and converge to a true solution.
Compared with the prior art, the invention has the following beneficial effects:
(1) only the magnetic field component in the horizontal X direction is collected, instruments are conveniently arranged (especially in a survey area with complex terrain), the construction efficiency is high, and the construction difficulty and the construction cost of the project are reduced.
(2) The acquisition device is a magnetic probe. Magnetic probes (e.g., induction coils) are small in size and have better pinpoint ability to detect subsurface signals than methods that measure electric fields. The horizontal resolution of the engineered survey to the subsurface electrical target zone is improved with the smallest possible measurement point spacing (typically twice the length of the magnetic probe).
(3) Only the magnetic probes for collecting the magnetic field components in the horizontal X direction are arranged, so that the electromagnetic mutual inductance effect between signal transmission lines is reduced to the maximum extent, and the anti-interference capability and the signal-to-noise ratio of collected data are improved. If the magnetic probe is designed to be tightly integrated with the digitized signal acquisition board (type without signal transmission line), the signal-to-noise ratio can be further improved.
(4) The observation mode of the horizontal X-direction magnetic component is based on the characteristics of low cost and strong expansion capacity of a set of instruments and hardware (magnetic probes) and the like, so that the data acquisition of all measuring points of a measuring area in a synchronous or short period is feasible, and the time efficiency is high (the time efficiency is important for improving the accuracy of dynamic monitoring of an oil reservoir and can also meet other engineering requirements with short construction period). In general, by properly adapting the capacity of the instrument suite, all frequency point data at all survey points in the survey area can be obtained within a pre-specified limited short time frame (e.g., 1 day to 1 week).
(5) Compared with the electric field component in the conventional electromagnetic method, the accuracy of distinguishing the underground electric target layer is improved by the geometric depth measuring effect, the static displacement effect and the volume effect when the low frequency does not exist in the horizontal X-direction magnetic component.
(6) The formula of the horizontal X-direction magnetic component in a uniform half-space medium is simple, so that the definition and calculation of accurate but not approximate 'full-area' apparent resistivity are easy to realize and stable in calculation, and a foundation is laid for solving key technical difficulties in the subsequent data processing and inversion processes.
(7) The method solves the technical bottleneck of the conventional electromagnetic method adopted when geological structure investigation is carried out in complex terrain areas, cold freezing surfaces, dry desert areas and seabed. Compared with the conventional electromagnetic method in which the measurement of the electric field component is extremely difficult, the distribution of the measurement of the magnetic component in the horizontal X direction is simple and easy.
(8) The method can improve the precision of shallow geological structure investigation. A strong primary field is excited based on a high-power artificial source, so that the human electromagnetic interference can be suppressed to the maximum extent, and a high signal-to-noise ratio signal is obtained. The horizontal resolution capability of the underground electric target layer can be ensured due to the fixed point of the magnetic probe. The method can get rid of the dilemma of the conventional shallow geological survey pattern (due to the limitation of the method, the ground penetrating radar, the vertical magnetic field measurement method of the overlapped loop and the high-density electrical method have the defects of low signal-to-noise ratio and low resolution capability).
(9) And designing a magnetic probe: different lengths and different frequency (band frequency) characteristics are selected according to different depths of exploration targets. Generally, for shallow exploration targets: f 105～102Hz, the exploration depth is about several to 100 meters; for the middle layer exploration target: f 102～100Hz, the exploration depth is about 100-1000 m; for deep exploration targets: f 100～10-3Hz and the exploration depth is about 1000-5000 meters. Generally, the higher the frequency band, the shorter the magnetic probe length. On the premise of mature technology, a universal magnetic probe with a passband can be designed.
(10) The method used by the invention is a frequency domain method, adopts a frequency sweeping mode, and ensures the consistency of the signal-to-noise ratio of the data of each frequency point. The method overcomes the defect that the long delay signal-to-noise ratio is lower than the short delay signal-to-noise ratio in the conventional time domain method.
The specific embodiment is as follows:
the selected exploration area is a certain fire flooding oil production area in Xinjiang, the oil deposit burial depth is 310.5-328.5 m, and the block ignition has been carried out for two years. The borehole logging shows that the formation resistivity of the burnt area can be 10 times higher than that of the normal unburnt area, and the physical property foundation for performing dynamic monitoring of the fireflood by applying the method is formed.
Around the ignition well, net-shaped measuring points 811 are arranged in total in a mode shown in figure 1, the distance between the measuring points is 2m-5m, the magnetic component in the horizontal X direction is measured, the dipole distance of an AB pole is 2000m, the transmitting and receiving distance of the central point of the measuring area is 4000m, the current is 20A, the frequency range is 1000 Hz to 0.01 Hz, and the recording time is 60 minutes/point.
By applying the technology, the technical difficulties of low construction efficiency, long period, space observation blind area, incapability of meeting production requirements on the longitudinal and transverse resolution capability of a geological target layer, weak industrial electromagnetic interference resistance, incapability of acquiring the visual resistivity of the whole area, low data processing precision and the like in similar exploration areas in the conventional technology are overcome.
The 'whole-area' apparent resistivity is obtained through field measurement, normalized function conversion and dichotomy numerical fitting. And finally obtaining the three-dimensional space resistivity data volume of the measuring area by adopting a proper inversion method. And extracting resistivity data of the ignition layer section along the layer, drawing a high resistivity abnormal distribution plan (figure 3) around the ignition well, and determining the range of the burnt area according to the plan. The result also obtains the indirect verification of the production dynamic data and the digital-analog data, and then the data are used for production regulation and control in the subsequent fireflood development, thereby obtaining the expected oil recovery and yield increase effect.
The present invention is not limited to the above preferred embodiments, and any modifications, equivalent replacements, improvements, etc. within the spirit and principle of the present invention should be included in the protection scope of the present invention.
Claims (8)
1. A method for observing the magnetic field component in the horizontal X direction of a frequency domain and processing data is characterized in that: the method comprises the following steps:
step 1: determining the horizontal X direction in the exploration area according to the specified method, and approximating
Determining the azimuth range as a data measurement area;
step 2: measuring the frequency spectrum data of the horizontal X-direction magnetic field component of the frequency domain one by using a magnetic probe one by one and frequency point by one in a measurement region by using an artificial source electric dipole excitation frequency domain electromagnetic wave field;
and step 3: processing the actually measured horizontal X-direction magnetic field component frequency spectrum data into an actually measured normalization function through conversion;
and 4, step 4: and (3) searching out points with the same theoretical normalization function and the actually measured normalization function by adopting a numerical value dichotomy calculation method, and further acquiring the apparent resistivity of the whole measuring point in the exploration area.
2. The method of claim 1, wherein: in step 1, the horizontal X-direction of the data measurement is defined as the direction parallel to the artificial emission source (AB).
3. The method of claim 1, wherein: in the step 1, the data measurement area falls on the azimuth
The range defined.
4. The method of claim 1, wherein: in the step 1, the measuring points are arranged along a straight line or in a certain plane area, the point distance of the measuring points is between 2 meters and 100 meters, and the transceiving distance r is between 4000 meters and 10000 meters.
5. The method of claim 1, wherein: in step 2, the measured data are: and (3) exciting a frequency domain signal by using an artificial source electric dipole, and measuring the frequency spectrum of the horizontal X-direction magnetic field component of the frequency domain one by using a magnetic probe one by one and one by one frequency point in a measurement region.
6. The method of claim 1, wherein: in step 3, the measured normalization function is calculated according to the frequency spectrum data of the horizontal X-direction magnetic field component in the frequency domain, and the conversion formula for calculation is:
7. the method of claim 1, wherein: in step 4, the calculation formula of the theoretical normalization function is:
8. the method of claim 1, wherein: in the step 4, the numerical value "dichotomy" calculation method includes:
(1) taking the initial two end values of the resistivity as
And
(2) finding out the actually measured normalization function by the dichotomy
And a theoretical normalization function
Points of equal value (kr)s＝nI.e. satisfy the equality condition
A point of (a);
(3) the apparent resistivity of the whole area is calculated by the following formula
(4) The specific calculation processes of 3 steps in the above (1) - (3) are repeated for j being 1
(5) And (3) circulating the specific calculation processes of 4 steps in the steps (1) to (4) for i to 1.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910905194.2A CN110927632A (en) | 2019-09-24 | 2019-09-24 | Frequency domain horizontal X-direction magnetic field component observation and data processing method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910905194.2A CN110927632A (en) | 2019-09-24 | 2019-09-24 | Frequency domain horizontal X-direction magnetic field component observation and data processing method |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN110927632A true CN110927632A (en) | 2020-03-27 |
Family
ID=69848822
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910905194.2A Withdrawn CN110927632A (en) | 2019-09-24 | 2019-09-24 | Frequency domain horizontal X-direction magnetic field component observation and data processing method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110927632A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115292660A (en) * | 2022-08-04 | 2022-11-04 | 中国自然资源航空物探遥感中心 | Adaptive iterative bit field separation method, system, equipment and computer readable storage medium |
Citations (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS6449995A (en) * | 1987-08-20 | 1989-02-27 | Maakurando Denshi Kogyo Kk | Electromagnetic surveying method |
| CN101216566A (en) * | 2008-01-15 | 2008-07-09 | 湖南继善高科技有限公司 | Couple source frequency domain electric sounding method and apparatus suitable for quasi far zone |
| CN101261327A (en) * | 2008-04-15 | 2008-09-10 | 何继善 | All-zone couple source frequency domain magnetic depth measuring method and device |
| CN102062876A (en) * | 2010-11-17 | 2011-05-18 | 汤井田 | Electrical sounding method for whole-region couple source frequency domain |
| CN102707323A (en) * | 2012-07-03 | 2012-10-03 | 长沙五维地科勘察技术有限责任公司 | Controllable source audio-frequency magnetic field sounding method for geological exploration |
| US20130319659A1 (en) * | 2010-11-19 | 2013-12-05 | Schlumberger Technology Corporation | Open-Hole Logging Instrument And Method For Making Ultra-Deep Magnetic And Resistivity Measurements |
| CN104360401A (en) * | 2014-12-10 | 2015-02-18 | 中国科学院地质与地球物理研究所 | Method for determining geological information of underground target body by transient electromagnetic field B |
| CN104597506A (en) * | 2015-01-26 | 2015-05-06 | 吉林大学 | Frequency domain ground-to-air electromagnetic prospecting method |
| CN104614774A (en) * | 2015-02-12 | 2015-05-13 | 中国科学院地质与地球物理研究所 | Transient electromagnetic detection method, device and system |
| CN105891892A (en) * | 2016-06-12 | 2016-08-24 | 何继善 | Device and method for whole-region measurement of vertical components of current source frequency domain magnetic fields |
| CN106291719A (en) * | 2016-08-09 | 2017-01-04 | 中南大学 | A kind of array artificial source's field frequency depth detecting method |
| CN108873083A (en) * | 2018-05-06 | 2018-11-23 | 东华理工大学 | A kind of artificial field source frequency domain electromagnetism apparent resistivity measurement method |
| CN109917470A (en) * | 2019-04-25 | 2019-06-21 | 中南大学 | A kind of vertical magnetism source wave number apparent resistivity measurement method and device |
-
2019
- 2019-09-24 CN CN201910905194.2A patent/CN110927632A/en not_active Withdrawn
Patent Citations (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS6449995A (en) * | 1987-08-20 | 1989-02-27 | Maakurando Denshi Kogyo Kk | Electromagnetic surveying method |
| CN101216566A (en) * | 2008-01-15 | 2008-07-09 | 湖南继善高科技有限公司 | Couple source frequency domain electric sounding method and apparatus suitable for quasi far zone |
| CN101261327A (en) * | 2008-04-15 | 2008-09-10 | 何继善 | All-zone couple source frequency domain magnetic depth measuring method and device |
| CN102062876A (en) * | 2010-11-17 | 2011-05-18 | 汤井田 | Electrical sounding method for whole-region couple source frequency domain |
| US20130319659A1 (en) * | 2010-11-19 | 2013-12-05 | Schlumberger Technology Corporation | Open-Hole Logging Instrument And Method For Making Ultra-Deep Magnetic And Resistivity Measurements |
| CN102707323A (en) * | 2012-07-03 | 2012-10-03 | 长沙五维地科勘察技术有限责任公司 | Controllable source audio-frequency magnetic field sounding method for geological exploration |
| CN104360401A (en) * | 2014-12-10 | 2015-02-18 | 中国科学院地质与地球物理研究所 | Method for determining geological information of underground target body by transient electromagnetic field B |
| CN104597506A (en) * | 2015-01-26 | 2015-05-06 | 吉林大学 | Frequency domain ground-to-air electromagnetic prospecting method |
| CN104614774A (en) * | 2015-02-12 | 2015-05-13 | 中国科学院地质与地球物理研究所 | Transient electromagnetic detection method, device and system |
| CN105891892A (en) * | 2016-06-12 | 2016-08-24 | 何继善 | Device and method for whole-region measurement of vertical components of current source frequency domain magnetic fields |
| CN106291719A (en) * | 2016-08-09 | 2017-01-04 | 中南大学 | A kind of array artificial source's field frequency depth detecting method |
| CN108873083A (en) * | 2018-05-06 | 2018-11-23 | 东华理工大学 | A kind of artificial field source frequency domain electromagnetism apparent resistivity measurement method |
| CN109917470A (en) * | 2019-04-25 | 2019-06-21 | 中南大学 | A kind of vertical magnetism source wave number apparent resistivity measurement method and device |
Non-Patent Citations (5)
| Title |
|---|
| 何继善: "广域电磁测深法研究", 《中南大学学报(自然科学版)》 * |
| 崔江伟: "电性源短偏移距瞬变电磁法全程视电阻率计算研究", 《中国优秀硕士学位论文全文数据库 基础科学辑》 * |
| 汤井田 等: "水平电偶源频率测深中全区视电阻率定义的新方法", 《地球物理学报》 * |
| 蒋奇云: "广域电磁测深仪关键技术研究", 《中国博士学位论文全文数据库 基础科学辑》 * |
| 闫国翔 等: "电性源瞬变电磁全区视电阻率定义", 《物探与化探》 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115292660A (en) * | 2022-08-04 | 2022-11-04 | 中国自然资源航空物探遥感中心 | Adaptive iterative bit field separation method, system, equipment and computer readable storage medium |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108873083B (en) | A kind of artificial field source frequency domain electromagnetism apparent resistivity measurement method | |
| CN102121374B (en) | Borehole time domain pulse electromagnetic method for detecting electrical resistivity of stratum out of metal bushing | |
| CN101520517B (en) | Method for accurately evaluating targets containing oil gas in clastic rock basin | |
| US10955578B2 (en) | Device and method for ground source transient electromagnetic near-field detection and related device | |
| CN104656157B (en) | A kind of method and device in identification shale gas dessert area | |
| CN110208866B (en) | Ground well array type optical fiber time-frequency electromagnetic data acquisition device and data acquisition method thereof | |
| CN102866417A (en) | Device and method for seismic cross hole computed tomography (CT) detection and tomography of underground cave | |
| CN104280782B (en) | One-dimensional joint inversion method for time-frequency electromagnetic data and magnetotelluric data | |
| CN105044792A (en) | Ground-well time-frequency electromagnetic exploration data acquisition apparatus and method | |
| CN202837558U (en) | Underground karst cave earthquake cross-hole CT (computer tomography) detection and tomographic imaging device | |
| CN111580170A (en) | Time-frequency electromagnetic data acquisition device and method based on extremely-low-frequency electromagnetic source | |
| CN109782360A (en) | A kind of low-resistance area of coverage deep sand body detection method | |
| CN114236621A (en) | Multi-technology cooperation based nondestructive detection method for potential hazards of dikes | |
| CN110187394A (en) | Double field source electromagnetic depth methods obtain the anisotropic method and device of formation resistivity | |
| CA1305216C (en) | Areal method of magnetotelluric exploration | |
| CN102182437B (en) | Method for determining and eliminating hydraulic fracture stress boundary of coal mine underground drilling | |
| CN105204073A (en) | Tensor apparent conductivity measurement method | |
| CN103323883A (en) | Near-field magnetic dipole source high-density geometric resistivity sounding method | |
| CN114217354A (en) | Electromagnetic data acquisition system and method based on optical fiber electromagnetic sensor | |
| CN212255732U (en) | Time-frequency electromagnetic data acquisition device based on extremely-low-frequency electromagnetic source | |
| CN110927632A (en) | Frequency domain horizontal X-direction magnetic field component observation and data processing method | |
| CN108761540A (en) | A kind of frequency domain natural electric field three-dimensional exploitation method | |
| CN210038187U (en) | Ground well array type optical fiber time-frequency electromagnetic data acquisition device | |
| RU2733095C2 (en) | Method of three-dimensional objects searching by tm-polarization geoelectrics methods | |
| CN110967767A (en) | Method for detecting karst cave by cross-hole induced polarization |
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 | ||
| WW01 | Invention patent application withdrawn after publication |
Application publication date: 20200327 |
|
| WW01 | Invention patent application withdrawn after publication |