CN108169800B - Controlled-source audiomagnetotellurics method apparent resistivity near-field calibrating method - Google Patents
Controlled-source audiomagnetotellurics method apparent resistivity near-field calibrating method Download PDFInfo
- Publication number
- CN108169800B CN108169800B CN201711443327.6A CN201711443327A CN108169800B CN 108169800 B CN108169800 B CN 108169800B CN 201711443327 A CN201711443327 A CN 201711443327A CN 108169800 B CN108169800 B CN 108169800B
- Authority
- CN
- China
- Prior art keywords
- data
- resistivity
- electric field
- asian
- field data
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- 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/38—Processing data, e.g. for analysis, for interpretation, for correction
-
- 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/083—Controlled source electromagnetic [CSEM] surveying
-
- 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/083—Controlled source electromagnetic [CSEM] surveying
- G01V2003/086—Processing
Abstract
The invention discloses a kind of controlled-source audiomagnetotellurics method apparent resistivity near-field calibrating methods, belong to the artificial field source frequency domain electromagnetic method data processing field of geophysical exploration, the method includes: step A, using Fourier transformation the current data of time-domain, electric field data, magnetic field data is transformed to frequency domain data；Step B, according to the electric field data and magnetic field data in frequency domain data, Ka Ni Asian TV Station resistivity and phase is calculated；Step C, near-field effect is judged whether there is according to Ka Ni Asian TV Station resistivity and phase property, such as there is near-field effect；Correction parameter is obtained by normalized theoretical the ratio between electric field data and the normalized electric field data of electric current, to correct the form partition of apparent resistivity curve, restore low-frequency data form, the data that can visit source audio magnetotelluric method are enable greatly to be applied in data inversion, to obtain more underground electrical structure information.
Description
Technical field
The present invention relates to the artificial field source frequency domain electromagnetic method data processing fields of geophysical exploration, more specifically,
The invention mainly relates to a kind of controlled-source audiomagnetotellurics method apparent resistivity near-field calibrating methods.
Background technique
Controlled-source audiomagnetotellurics method used at present cannot work near field region, limit controllable source audio the earth
The exploration construction and data processing of electromagnetic method.Currently, the work of controlled-source audiomagnetotellurics method is carried out, it is ensured that transmitting-receiving is away from big
In 3 to 5 times of skin depth, to achieve the purpose that far field is observed.But since underground apparent resistivity is our exploration targets, because
This, skin depth is a unknown quantity, can only generally be estimated according to existing data.This has resulted in practical transmitting-receiving away from inclined
It is small, it can not meet far zone condition, finally obtained apparent resistivity curve has apparent near-field effect, and (apparent resistivity curve is in 45 degree
Rise, 0) phase levels off to.At this moment the apparent resistivity curve obtained, cannot truly corresponsively under electrical structure.Therefore, it needs
A kind of near-field calibrating method is researched and developed, the true apparent resistivity information in apparent resistivity curve is restored.
Summary of the invention
One of the aims of the present invention is to provide...in response to the above drawbacks a kind of controlled-source audiomagnetotellurics method apparent resistivities
Near-field calibrating method is likely to occur near-field effect generation it is expected to solve controlled-source audiomagnetotellurics method apparent resistivity curve
Form partition causes the technical problems such as 45 degree of low-frequency data risings.
In order to solve the above technical problems, the invention adopts the following technical scheme:
A kind of controlled-source audiomagnetotellurics method apparent resistivity near-field calibrating method provided by the present invention, the method
Include the following steps:
Step A, the current data of time-domain, electric field data, magnetic field data are transformed to frequency domain using Fourier transformation
Data；
Step B, according to the electric field data and magnetic field data in frequency domain data, Ka Ni Asian TV Station resistivity and phase is calculated
Position；
Step C, near-field effect is judged whether there is according to Ka Ni Asian TV Station resistivity and phase property, such as there is near field effect
It answers, then continues next step, on the contrary then step terminates；
Step D, using the current data of frequency domain, electric current normalization is carried out to the electric field data of frequency domain, is obtained
The normalized electric field data of the electric current of frequency domain；
Step E, according to the Ka Ni Asian TV Station resistivity for having near-field effect, normalized theoretical electric field data is calculated；
Step F, according to normalized theoretical the ratio between electric field data and the normalized electric field data of electric current, correction ginseng is obtained
Number, is corrected the Ka Ni Asian TV Station resistivity data for having near-field effect using correction parameter.
Preferably, further technical solution is: the current data of time-domain, electric field data, magnetic field in the step A
Data are to survey to obtain by instrument.
Further technical solution is: judging whether to deposit to Ka Ni Asian TV Station resistivity and phase property in the step C
In near-field effect, to judge according to whether apparent resistivity rises at 45 degree with phase saturation.
Further technical solution is: the method further includes step G, according to the Ka Ni Asian TV Station resistivity after correction
Data corrected after Ka Ni Asian TV Station resistivity curve.
Further technical solution is: the method further includes step H, according to the Ka Ni Asian TV Station resistivity after correction
Curve calculates phase, thus the phase after being corrected.
Further technical solution is: normalized theoretical electric field data is calculate by the following formula in step E:
In formula,For normalized theoretical electric field data, Idl is unit electric dipole moment, and ρ is Ka Ni Asian TV Station resistivity, and e is
Natural constant, i are imaginary unit, and r is transmitting-receiving away from k is wave number；Eelctric dipole is enabled to be located at coordinate origin, measuring point is located at (x, y), formula
Middle y is the ordinate of measuring point.
Further technical solution is: the step D carries out electric current normalizing by electric field data of the following formula to frequency domain
Change correction:
In formula,For the normalized electric field data of electric current, ExFor the electric field data of frequency domain,It is normalized for electric current
Magnetic field data, I are the current data of frequency domain, HyFor the magnetic field data of frequency domain.
Further technical solution is: the step B is calculate by the following formula Ka Ni Asian TV Station resistivity:
In formula, ρ is Ka Ni Asian TV Station resistivity, and μ is magnetic conductivity, and ω is angular frequency, ExFor the electric field data of frequency domain, HyFor
The magnetic field data of frequency domain.
Compared with prior art, beneficial effects of the present invention first is that: pass through normalized theoretical electric field data and electric current
The ratio between normalized electric field data obtains correction parameter, to correct the form partition of apparent resistivity curve, restores low-frequency data
Form enables the data that can visit source audio magnetotelluric method to be greatly applied in data inversion, to obtain more
More underground electrical structure information, while a kind of controlled-source audiomagnetotellurics method apparent resistivity near field school provided by the present invention
Correction method is easily achieved, and the resistivity that can be surveyed in all kinds of places is corrected, and application range is wide.
Detailed description of the invention
Fig. 1 is for illustrating the electric field of model 1 and apparent resistivity calibration coordinate figure in one embodiment of the invention；
Fig. 2 is the phasing coordinate diagram for illustrating model 1 in one embodiment of the invention；
Fig. 3 is for illustrating the electric field of model 2 and apparent resistivity calibration coordinate figure in one embodiment of the invention；
Fig. 4 is the phasing coordinate diagram for illustrating model 2 in one embodiment of the invention；
Specific embodiment
The present invention is further elaborated with reference to the accompanying drawing.
One embodiment of the present of invention is a kind of controlled-source audiomagnetotellurics method apparent resistivity near-field calibrating method, the party
Method includes the following steps:
Step S1, the current data of time-domain, electric field data, magnetic field data are transformed to frequency domain using Fourier transformation
Data；
In this step, current data above-mentioned, electric field data, magnetic field data are both needed to survey to obtain by instrument；
Step S2, according to the electric field data and magnetic field data in frequency domain data, be calculated Ka Ni Asian TV Station resistivity and
Phase；
In this step, Ka Ni Asian TV Station resistivity can be calculate by the following formula:
In formula, ρ is Ka Ni Asian TV Station resistivity, and μ is magnetic conductivity, and ω is angular frequency, ExFor the electric field data of frequency domain, HyFor
The magnetic field data of frequency domain.
Step S3, near-field effect is judged whether there is according to Ka Ni Asian TV Station resistivity and phase property, such as there is near field effect
It answers, then continues next step, on the contrary then step terminates；
In this step, near-field effect is judged whether there is to Ka Ni Asian TV Station resistivity and phase property, basis can be used
Whether apparent resistivity rises at 45 degree and phase is saturated to judge；
Step S4, using the current data of frequency domain, electric current normalization is carried out to the electric field data of frequency domain, is obtained
The normalized electric field data of the electric current of frequency domain；
In this step, electric current normalization can be carried out by electric field data of the following formula to frequency domain:
In formula,For the normalized electric field data of electric current, ExFor the electric field data of frequency domain,It is normalized for electric current
Magnetic field data, I are the current data of frequency domain, HyFor the magnetic field data of frequency domain.
Step S5, according to the Ka Ni Asian TV Station resistivity for having near-field effect, normalized theoretical electric field data is calculated；
In this step, normalized theoretical electric field data can be calculate by the following formula:
In formula,For normalized theoretical electric field data, Idl is unit electric dipole moment, and ρ is Ka Ni Asian TV Station resistivity, e
For natural constant, i is imaginary unit, and r is transmitting-receiving away from k is wave number；Eelctric dipole is enabled to be located at coordinate origin, measuring point is located at (x, y),
Y is the ordinate of measuring point in formula.
Step S6, according to normalized theoretical the ratio between electric field data and the normalized electric field data of electric current, correction ginseng is obtained
Number, is corrected the Ka Ni Asian TV Station resistivity data for having near-field effect using correction parameter.
According to another embodiment of the invention, it can also continue to carry out following steps:
Step S7, the Ka Ni Asian TV Station resistivity curve after being corrected according to the Ka Ni Asian TV Station resistivity data after correction；
Step S8, phase is calculated according to the Ka Ni Asian TV Station resistivity curve after correction, thus the phase after being corrected.
Due to above-described embodiment it is found that method provided by the present invention substantially step are as follows:
The first step, when controlled-source audiomagnetotellurics method field data collection, electric current transmitted by record current transmitter
Signal i (t), the geographical position coordinates (x, y) of transmitter, the geographical position coordinates (m, n) of receiver, orthogonal electric field Ex(t)
With magnetic field Hy(t) signal, thus to obtain field initial data.
Second step, by the current signal i (t) of time-domain, electric field Ex(t) and magnetic field Hy(t) signal is turned by Fourier Tranform
It is changed to the current signal I of frequency domain, electric field ExWith magnetic field Hy。
Third step, to electric field ExWith magnetic field HyCarry out current correction:
4th step calculates Ka Ni Asian TV Station resistivity
5th step, according to Ka Niya computation of apparent resistivity theory electric field data
6th step, according to theoretical electric field dataWith the electric field data of correctionThe ratio between obtain correction parameter.
7th step, according to correction parameter correction card Buddhist nun Asian TV Station resistivity.
Based on the above-mentioned method of the present invention, carry out proving program using two-layer model:
Verify model 1: middle layer is resistive formation, and resistivity is respectively: 10 Ω m, 100 Ω m, 10 Ω m, thickness point
Other 1000m, 2000m are received and dispatched away from r=10000m.Wherein ExcFor theoretical electric field, ExFor survey electric field, verifying coordinate diagram such as Fig. 1 with
Shown in Fig. 2.
Verify model 2: middle layer is low resistivity layer.Resistivity is respectively: 100 Ω m, 10 Ω m, 100 Ω m, thickness
It is respectively as follows: 1000m, 2000m, it is as shown in Figure 3 and Figure 4 to verify coordinate diagram away from r=10000m for transmitting-receiving.
Than that described above, it is also necessary to which explanation is " one embodiment " spoken of in the present specification, " another implementation
Example ", " embodiment " etc. refer to that specific features, structure or the feature of embodiment description is combined to be included in the application summary
Property description at least one embodiment in.It is not centainly to refer to the same reality that statement of the same race, which occur, in multiple places in the description
Apply example.Furthermore, it is understood that is advocated is knot when describing a specific features, structure or feature in conjunction with any embodiment
Other embodiments are closed to realize that this feature, structure or feature are also fallen within the scope of the present invention.
Although reference be made herein to invention has been described for multiple explanatory embodiments of the invention, however, it is to be understood that
Those skilled in the art can be designed that a lot of other modification and implementations, these modifications and implementations will fall in this Shen
It please be within disclosed scope and spirit.More specifically, disclose in the application, drawings and claims in the range of, can
With the building block and/or a variety of variations and modifications of layout progress to theme combination layout.In addition to building block and/or layout
Outside the variations and modifications of progress, to those skilled in the art, other purposes also be will be apparent.
Claims (7)
1. a kind of controlled-source audiomagnetotellurics method apparent resistivity near-field calibrating method, it is characterised in that the method includes such as
Lower step:
Step A, the current data of time-domain, electric field data, magnetic field data are transformed to frequency domain data using Fourier transformation；
Step B, according to the electric field data and magnetic field data in frequency domain data, Ka Ni Asian TV Station resistivity and phase is calculated；
Step C, near-field effect is judged whether there is according to Ka Ni Asian TV Station resistivity and phase property, such as there is near-field effect, then
Continue next step, on the contrary then step terminates；
Step D, using the current data of frequency domain, electric current normalization is carried out to the electric field data of frequency domain, obtains frequency
The normalized electric field data of the electric current in domain；
Step E, it according to the Ka Ni Asian TV Station resistivity for having near-field effect, is calculate by the following formula and calculates normalized theoretical electric field number
According to；
In formula,For normalized theoretical electric field data, Idl is unit electric dipole moment, and ρ is Ka Ni Asian TV Station resistivity, and e is nature
Constant, i are imaginary unit, and r is transmitting-receiving away from k is wave number；Eelctric dipole is enabled to be located at coordinate origin, measuring point is located at (x, y), wherein y
For the ordinate of measuring point；
Step F, according to normalized theoretical the ratio between electric field data and the normalized electric field data of electric current, correction parameter is obtained, benefit
The Ka Ni Asian TV Station resistivity data for having near-field effect is corrected with correction parameter.
2. controlled-source audiomagnetotellurics method apparent resistivity near-field calibrating method according to claim 1, it is characterised in that:
The current data of time-domain, electric field data, magnetic field data are to survey to obtain by instrument in the step A.
3. controlled-source audiomagnetotellurics method apparent resistivity near-field calibrating method according to claim 1, it is characterised in that:
Near-field effect is judged whether there is to Ka Ni Asian TV Station resistivity and phase property in the step C, for according to apparent resistivity whether
Rise at 45 degree and phase is saturated to judge.
4. controlled-source audiomagnetotellurics method apparent resistivity near-field calibrating method according to claim 1, it is characterised in that:
The method further include step G, corrected according to the Ka Ni Asian TV Station resistivity data after correction after Ka Ni Asian TV Station resistance
Rate curve.
5. controlled-source audiomagnetotellurics method apparent resistivity near-field calibrating method according to claim 4, it is characterised in that:
The method further includes step H, according to the Ka Ni Asian TV Station resistivity curve calculating phase after correction, thus after being corrected
Phase.
6. controlled-source audiomagnetotellurics method apparent resistivity near-field calibrating method according to claim 1, it is characterised in that
The step D carries out electric current normalization by electric field data of the following formula to frequency domain:
In formula,For the normalized electric field data of electric current, ExFor the electric field data of frequency domain,For the normalized magnetic field of electric current
Data, I are the current data of frequency domain, HyFor the magnetic field data of frequency domain.
7. controlled-source audiomagnetotellurics method apparent resistivity near-field calibrating method according to claim 6, it is characterised in that
The step B is calculate by the following formula Ka Ni Asian TV Station resistivity:
In formula, ρ is Ka Ni Asian TV Station resistivity, and μ is magnetic conductivity, and ω is angular frequency, ExFor the electric field data of frequency domain, HyFor frequency
The magnetic field data in domain.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711443327.6A CN108169800B (en) | 2017-12-27 | 2017-12-27 | Controlled-source audiomagnetotellurics method apparent resistivity near-field calibrating method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711443327.6A CN108169800B (en) | 2017-12-27 | 2017-12-27 | Controlled-source audiomagnetotellurics method apparent resistivity near-field calibrating method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108169800A CN108169800A (en) | 2018-06-15 |
| CN108169800B true CN108169800B (en) | 2019-10-11 |
Family
ID=62521967
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711443327.6A Active CN108169800B (en) | 2017-12-27 | 2017-12-27 | Controlled-source audiomagnetotellurics method apparent resistivity near-field calibrating method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108169800B (en) |
Families Citing this family (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109470135A (en) * | 2018-11-12 | 2019-03-15 | 吉林大学 | CSAMT data inactivity bearing calibration |
| CN109884714B (en) * | 2019-03-05 | 2021-11-16 | 中国地质科学院地球物理地球化学勘查研究所 | Controllable source electromagnetic measurement method, device and storage medium thereof |
| CN110348568B (en) * | 2019-07-16 | 2021-02-23 | 山东科技大学 | Deep goaf detection method suitable for strong electromagnetic interference area |
| CN110989006B (en) * | 2019-12-25 | 2022-07-22 | 中国地质科学院地球物理地球化学勘查研究所 | Static displacement correction method and device in controllable source electromagnetic method measurement and intelligent terminal |
| CN111273356B (en) * | 2020-03-29 | 2021-04-16 | 中南大学 | Active power conduction method based on monitoring potential correction |
| CN111693881B (en) * | 2020-05-28 | 2021-09-10 | 江苏大学 | Method for estimating state of health of battery on line under wide temperature range based on standardized temperature |
| CN111965712B (en) * | 2020-10-21 | 2021-03-02 | 国网江西省电力有限公司电力科学研究院 | Method for correcting static effect of controllable source audio magnetotelluric method |
| CN112379449B (en) * | 2020-10-30 | 2023-05-26 | 中国石油天然气集团有限公司 | Processing method and device for electromagnetic data of controllable source |
| CN114047551A (en) * | 2021-10-13 | 2022-02-15 | 中国舰船研究院 | Real-time monitoring implementation method based on WEM data |
| CN115437024B (en) * | 2022-11-07 | 2022-12-30 | 四川省冶勘设计集团有限公司 | Artificial source frequency domain electromagnetic field signal of double-parameter processing method of plane wave data |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CA2611006C (en) * | 2005-06-09 | 2014-09-30 | Exxonmobil Upstream Research Company | Method for determining earth vertical electrical anisotropy in marine electromagnetic surveys |
| CN104656156A (en) * | 2013-11-22 | 2015-05-27 | 中国石油天然气集团公司 | Magnetic reference processing method for acquiring data in three dimensions audio-frequency magnetotelluric sounding |
| US10458230B2 (en) * | 2015-02-27 | 2019-10-29 | Halliburton Energy Services, Inc. | Formation resistivity measurement apparatus, systems, and methods |
| CN105911595A (en) * | 2016-02-02 | 2016-08-31 | 中国科学院地质与地球物理研究所 | Method and apparatus for obtaining controllable source audio-frequency magnetotelluric (CSAMT) apparent phase information |
-
2017
- 2017-12-27 CN CN201711443327.6A patent/CN108169800B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN108169800A (en) | 2018-06-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108169800B (en) | Controlled-source audiomagnetotellurics method apparent resistivity near-field calibrating method | |
| CN101573635B (en) | Systems and methods for measuring sea-bed resistivity | |
| Yu-Guo et al. | Transient electromagnetic in shallow water: insights from 1D modeling | |
| CN104375195B (en) | Many source multi-component three-dimensional joint inversion methods of time-frequency electromagnetism | |
| CN110058315B (en) | Three-dimensional anisotropic radio frequency magnetotelluric adaptive finite element forward modeling method | |
| CN106873041B (en) | A method of apparent resistivity is obtained by any level electric field component | |
| CN108873083A (en) | A kind of artificial field source frequency domain electromagnetism apparent resistivity measurement method | |
| NO342941B1 (en) | Method and computer product for identifying a boundary of a geological feature, as well as a drill string management system | |
| CN102707323A (en) | Controllable source audio-frequency magnetic field sounding method for geological exploration | |
| CN113051779B (en) | Numerical simulation method of three-dimensional direct-current resistivity method | |
| CN106291719A (en) | A kind of array artificial source's field frequency depth detecting method | |
| Jang et al. | Step-off, vertical electromagnetic responses of a deep resistivity layer buried in marine sediments | |
| CN109407160A (en) | A method of D integral pin-fin tube is carried out based on WEM vector finite element method | |
| CN104088630B (en) | With boring the method that azimuthal resistivity logger obtains stratigraphic boundary distance | |
| Hou et al. | Comparison between different apparent resistivity definitions of CSAMT | |
| US9020205B2 (en) | Methods of multinary inversion for imaging objects with discrete physical properties | |
| WO2021096877A1 (en) | Removing the effect of near-surface inhomogeneities in surface-to-borehole measurements | |
| CN111965712A (en) | Method for correcting static effect of controllable source audio magnetotelluric method | |
| CN109416411A (en) | Anisotropy and stratigraphic dip are determined using electromagnetism (EM) pulse from inclined antenna | |
| Li et al. | Calibration method of azimuth electromagnetic boundary detection logging while drilling based on air–seawater double-layer dielectric as reflection interface | |
| ZHANG et al. | Calibration study on the damping characteristics of TEM sensor | |
| Peacock et al. | Repeat magnetotelluric measurements to monitor The Geysers steam field in northern California | |
| CAI | Calculation and application of B-field for time-domain airborne EM | |
| Watts et al. | Limitations of MT static shift corrections using time-domain EM data. | |
| CN113566686A (en) | Method and device for verifying buried depth position based on ultra-large buried depth pipeline |
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 |