CN104422970A - Amplitude-preserving static correcting method for electromagnetic data - Google Patents
Amplitude-preserving static correcting method for electromagnetic data Download PDFInfo
- Publication number
- CN104422970A CN104422970A CN201310397685.3A CN201310397685A CN104422970A CN 104422970 A CN104422970 A CN 104422970A CN 201310397685 A CN201310397685 A CN 201310397685A CN 104422970 A CN104422970 A CN 104422970A
- Authority
- CN
- China
- Prior art keywords
- rho
- data
- apparent resistivity
- log
- curves
- 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.)
- Granted
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
- Measurement Of Resistance Or Impedance (AREA)
Abstract
The invention belongs to the technology of electromagnetic exploitation data processing. The technology comprises the following steps of collecting electromagnetic data, determining an electric layer with stable depth and resistivity, averaging two branches of apparent resistivity curves of all detecting points on an electromagnetic data detecting line on a corresponding frequency range, and calculating average apparent resistivity data; obtaining the primary correction factors of two branches of curves of each detecting point, correcting the primary correction factors to obtain data of the two branches of apparent resistivity curve, selecting a frequency range again, respectively calculating the average values of apparent resistivity of the two branches of curves of each detecting point within the frequency range before filtration and after filtration, obtaining amplitude-preserving correction factors of the two branches of curves of each detecting point, and calculating to obtain the apparent resistivity data of the two branches of curves after final amplitude-preserving correction. The defect that the configuration of the apparent resistivity curves is changed after the traditional filtration method is used can be overcome, particularly at the places, such as a mountain front and an earth surface lithologic character mutation region, the effect is obviously improved, and the aim of keeping amplitude is fulfilled.
Description
Technical field
The present invention relates to geophysical exploration method, is a kind of electromagnetic survey data, comprises the data processing technique of mt (MT), audio magnetotelluric method (AMT) and continuous profiling mt (CEMP).
Background technology
Static relays in electromagnetic prospecting refer to due to closely show to exist local electrically heterogeneous body time, on resistivity interphase, the heap of polarization charge causes the distortion of electric field, produces an additional electric field be directly proportional to external electric field thus, and irrelevant with frequency.For magnetotelluric data (MT), show on single-point curve, be exactly that log-domain resistivity curve produces translation along the longitudinal axis, phase curve is substantially unaffected.Show on apparent resistivity section diagram, be exactly that upright abrupt change band appears in resistivity value, be commonly called as " vermicelli bar " phenomenon, and corresponding phase place pseudosection map does not just have this phenomenon.Showing on the apparent resistivity planimetric map of single-frequency point, is exactly there is a lot " distortional point " phenomenon.The intensity of static relays can reach two or three orders of magnitude, can cause very big error, and making structure elucidation complicated, therefore needing to carry out static correction before inversion interpretation when inferring the degree of depth.
At present, the method of static correction is a lot, in order to suppress static relays, improves constantly the effect of MT, lot of domestic and international expert has delivered correlative theses, proposes filter method, theoretical calculation, transient electromagnetic correction method, subsurface resistivities correction method, high-frequency electromagnetic function upward continuation method etc.These method major parts are all based on two-dimentional MT static correction, and the principle of various bearing calibration is different with to the requirement of data, are all subject to restriction to a certain extent in actual applications.Because filter method is easy to realize, processing speed is fast, is most widely used in the static corrections processing of MT/CEMP data.But after filtering after method process, the apparent resistivity curve form of a lot of measuring point can be changed, especially in local performances such as mountain front, Lithology sudden change bands obviously, in practice, static correction effect is bad.
Summary of the invention
The present invention is directed to and adopt traditional filter method to carry out the bad defect of effect after electromagnetic data static correction, provide a kind of can overcome apparent resistivity curve form can reformed defect, reach guarantor's width static correcting method of the object electromagnetic data of guarantor's width.
The present invention adopts following steps to realize:
1) according to the magnetotelluric data that work area gathers, determine a degree of depth and the stable electric layers of resistivity, obtain frequency range corresponding on electromagnetic data curve, do on average to two apparent resistivity curves of measuring points all on electromagnetic data survey line in frequency range, calculate Average apparent resistivity data
Described electromagnetic data comprises magnetotelluric data (MT), audio-frequency magnetotelluric magnetic data (AMT), continuous profiling magnetotelluric data (CEMP).
Described on average comprises arithmetic mean and geometric mean.
Described calculating Average apparent resistivity data
be:
Wherein,
be respectively i-th measuring point jth frequency f
jxY and YX pattern measured visual resistivity value, n and m is selected frequency point number.
2) following formula is utilized to ask for the preliminary corrections factor of each measuring point two curves;
Wherein,
with
be respectively the preliminary corrections factor of i-th measuring point XY and YX mode data,
the calculating mean value that (4) formula of serving as reasons obtains, the implication of other parameter is identical with (4) formula.
3) carry out preliminary corrections according to following formula, obtain two apparent resistivity curve data after correcting
with
Right again
carry out filtering respectively, obtain filtered apparent resistivity curve data
with
the implication of other parameter and (4), (5) formula are identical.
Described filtering comprises two dimension median filter, two-dimensional space filtering, three-dimensional medium filtering and three dimensions filtering.
4) frequency band is reselected, according to step 1) to 3) calculate before filtering respectively and the arithmetic of each measuring point two curves apparent resistivity in this band limits or geometrical mean after filtering
with
Wherein, nf1 and nf2 is respectively selected frequency point number,
be respectively a si measuring point jth frequency f
jxY and YX pattern measured visual resistivity value,
be respectively apparent resistivity arithmetic or the geometrical mean of XY and YX pattern before the filtering of calculating,
with
be respectively apparent resistivity arithmetic or the geometrical mean of XY and YX pattern after the filtering of calculating.
5) guarantor's width correction factor of each measuring point two curves is asked for according to following formula:
Wherein,
with
be respectively guarantor's width correction factor of XY and the YX curve of i-th measuring point.
6) obtain final guarantor's width according to following formulae discovery and correct rear two curve apparent resistivity datas
with
complete static correction.
Wherein,
with
the guarantor's width being exactly final XY and YX curve corrects rear apparent resistivity data.
Compared with background technology, the present invention can to overcome after the traditional filtering method of application apparent resistivity curve form can reformed defect, especially in places such as mountain front, Lithology sudden change bands, obviously effect improved, reaches the object of guarantor's width.
Accompanying drawing explanation
Fig. 1 is then for electromagnetic data protects width static correcting method basic flow sheet.
Embodiment
Be below the specific embodiment of the invention, Fig. 1 is the process flow diagram protecting the enforcement of width static correcting method.
Concrete grammar is:
1) according to the magnetotelluric data that work area gathers, determine a degree of depth and the stable electric layers of resistivity, obtain frequency range corresponding on electromagnetic data curve, as chosen No. 5 to No. 9 frequency band, then n=5, m=9 in formula (4).Do on average to two apparent resistivity curves of measuring points all on electromagnetic data survey line in frequency range, calculate Average apparent resistivity data
Described electromagnetic data comprises magnetotelluric data (MT), audio-frequency magnetotelluric magnetic data (AMT), continuous profiling magnetotelluric data (CEMP).
Described on average comprises arithmetic mean and geometric mean.
Described calculating Average apparent resistivity data
be:
Wherein,
be respectively i-th measuring point jth frequency f
jxY and YX pattern measured visual resistivity value, n and m is selected frequency point number, is respectively 5 and 9 in the present example.
2) following formula is utilized to ask for the preliminary corrections factor of each measuring point two curves;
Wherein,
with
be respectively the preliminary corrections factor of i-th measuring point XY and YX mode data,
the calculating mean value that (4) formula of serving as reasons obtains, the implication of other parameter is identical with (4) formula.
3) carry out preliminary corrections according to following formula, obtain two apparent resistivity curve data after correcting
with
Right again
carry out filtering respectively, as in the present example, adopt two dimension median filter, medium filtering window selection is (5,3), obtains filtered apparent resistivity curve data
with
the implication of other parameter and (4), (5) formula are identical.
Described filtering comprises two dimension median filter, two-dimensional space filtering, three-dimensional medium filtering and three dimensions filtering.
4) frequency band is reselected, if No. 7 frequency is to No. 13 frequency band, namely the nf1=7 in formula (7) and formula (8), nf2=13, according to step 1) to 3) calculate before filtering respectively and the arithmetic of each measuring point two curves apparent resistivity in this band limits or geometrical mean after filtering
with
Wherein, nf1 and nf2 is respectively selected frequency point number, is respectively 7 and 13 in the present example,
be respectively a si measuring point jth frequency f
jxY and YX pattern measured visual resistivity value,
be respectively apparent resistivity arithmetic or the geometrical mean of XY and YX pattern before the filtering of calculating,
with
be respectively apparent resistivity arithmetic or the geometrical mean of XY and YX pattern after the filtering of calculating.
5) guarantor's width correction factor of each measuring point two curves is asked for according to (9) with the formula:
Wherein,
with
be respectively guarantor's width correction factor of XY and the YX curve of i-th measuring point.
6) calculate the final guarantor's width of acquisition according to formula (10) and correct rear two curve apparent resistivity datas
with
complete and protect width static correction.
Wherein,
with
the guarantor's width being exactly final XY and YX curve corrects rear apparent resistivity data.
Claims (5)
1. guarantor's width static correcting method of electromagnetic data, feature adopts following steps to realize:
1) according to the magnetotelluric data that work area gathers, determine a degree of depth and the stable electric layers of resistivity, obtain frequency range corresponding on electromagnetic data curve, do on average to two apparent resistivity curves of measuring points all on electromagnetic data survey line in frequency range, calculate Average apparent resistivity data
2) following formula is utilized to ask for the preliminary corrections factor of each measuring point two curves;
Wherein,
with
be respectively the preliminary corrections factor of i-th measuring point XY and YX mode data,
the calculating mean value that (4) formula of serving as reasons obtains, the implication of other parameter is identical with (4) formula;
3) carry out preliminary corrections according to following formula, obtain two apparent resistivity curve data after correcting
with
Right again
carry out filtering respectively, obtain filtered apparent resistivity curve data
with
the implication of other parameter and (4), (5) formula are identical;
4) frequency band is reselected, according to step 1) to 3) calculate before filtering respectively and the arithmetic of each measuring point two curves apparent resistivity in this band limits or geometrical mean after filtering
with
Wherein, nf1 and nf2 is respectively selected frequency point number,
be respectively a si measuring point jth frequency f
jxY and YX pattern measured visual resistivity value,
be respectively apparent resistivity arithmetic or the geometrical mean of XY and YX pattern before the filtering of calculating,
with
be respectively apparent resistivity arithmetic or the geometrical mean of XY and YX pattern after the filtering of calculating;
5) guarantor's width correction factor of each measuring point two curves is asked for according to following formula:
Wherein,
with
be respectively guarantor's width correction factor of XY and the YX curve of i-th measuring point;
6) obtain final guarantor's width according to following formulae discovery and correct rear two curve apparent resistivity datas
with
complete static correction;
Wherein,
with
the guarantor's width being exactly final XY and YX curve corrects rear apparent resistivity data.
2. method according to claim 1, feature is that the electromagnetic data described in step 1) comprises magnetotelluric data, audio-frequency magnetotelluric magnetic data, continuous profiling magnetotelluric data.
3. method according to claim 1, feature on average comprises arithmetic mean and geometric mean described in step 1).
4. method according to claim 1, feature is the calculating Average apparent resistivity data described in step 1)
be:
Wherein,
be respectively i-th measuring point jth frequency f
jxY and YX pattern measured visual resistivity value, n and m is selected frequency point number.
5. method according to claim 1, feature is that the filtering described in step 3) comprises two dimension median filter, two-dimensional space filtering, three-dimensional medium filtering and three dimensions filtering.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310397685.3A CN104422970B (en) | 2013-09-04 | 2013-09-04 | A kind of guarantor width static correcting method of electromagnetic data |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310397685.3A CN104422970B (en) | 2013-09-04 | 2013-09-04 | A kind of guarantor width static correcting method of electromagnetic data |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104422970A true CN104422970A (en) | 2015-03-18 |
CN104422970B CN104422970B (en) | 2017-03-15 |
Family
ID=52972515
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310397685.3A Active CN104422970B (en) | 2013-09-04 | 2013-09-04 | A kind of guarantor width static correcting method of electromagnetic data |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104422970B (en) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104903750A (en) * | 2014-07-23 | 2015-09-09 | 王雅苹 | Amplitude preservation statics method of electromagnetic data |
CN110488365A (en) * | 2019-09-11 | 2019-11-22 | 湖南五维地质科技有限公司 | A kind of multipolarization telluric electromagnetic sounding method |
CN111965712A (en) * | 2020-10-21 | 2020-11-20 | 国网江西省电力有限公司电力科学研究院 | Method for correcting static effect of controllable source audio magnetotelluric method |
CN113484920A (en) * | 2021-08-17 | 2021-10-08 | 成都理工大学 | Two-dimensional structured inversion method for frequency domain electromagnetic sounding data |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4392109A (en) * | 1980-11-14 | 1983-07-05 | Mobil Oil Corporation | Static correction for magnetotelluric data |
CN102495431A (en) * | 2011-11-16 | 2012-06-13 | 中煤科工集团西安研究院 | Data processing method for carrying out static correction on transient electromagnetic data |
-
2013
- 2013-09-04 CN CN201310397685.3A patent/CN104422970B/en active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4392109A (en) * | 1980-11-14 | 1983-07-05 | Mobil Oil Corporation | Static correction for magnetotelluric data |
CN102495431A (en) * | 2011-11-16 | 2012-06-13 | 中煤科工集团西安研究院 | Data processing method for carrying out static correction on transient electromagnetic data |
Non-Patent Citations (6)
Title |
---|
ALAN G. JONES: "Static shift of magnetotelluric data and its removal in a sedimentary basin environment", 《GEOPHYSICS》 * |
WILLIAM CUMMING ET AL.: "Resistivity Imaging of Geothermal Resources Using 1D, 2D and 3D MT Inversion and TDEM Static Shift Correction Illustrated by a Glass Mountain Case History", 《PROCEEDINGS WORLD GEOTHERMAL CONGRESS 2010 》 * |
杨妮妮,等: "两种新方法实现CSAMT静态校正", 《华南地震》 * |
汤井田,等: "静效应校正的波数域滤波方法", 《物探与化探》 * |
罗延钟,等: "可控源音频大地电磁法的静态效应校正", 《物探与化探》 * |
罗志琼: "用电磁阵列剖面法压制MT静态效应影响的研究", 《地球科学 中国地质大学学报》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104903750A (en) * | 2014-07-23 | 2015-09-09 | 王雅苹 | Amplitude preservation statics method of electromagnetic data |
WO2016011627A1 (en) * | 2014-07-23 | 2016-01-28 | 王雅苹 | Amplitude-preserved static correction method for magnetic data |
CN110488365A (en) * | 2019-09-11 | 2019-11-22 | 湖南五维地质科技有限公司 | A kind of multipolarization telluric electromagnetic sounding method |
CN110488365B (en) * | 2019-09-11 | 2021-05-07 | 湖南五维地质科技有限公司 | Multi-polarization magnetotelluric depth measurement method |
CN111965712A (en) * | 2020-10-21 | 2020-11-20 | 国网江西省电力有限公司电力科学研究院 | Method for correcting static effect of controllable source audio magnetotelluric method |
CN111965712B (en) * | 2020-10-21 | 2021-03-02 | 国网江西省电力有限公司电力科学研究院 | Method for correcting static effect of controllable source audio magnetotelluric method |
CN113484920A (en) * | 2021-08-17 | 2021-10-08 | 成都理工大学 | Two-dimensional structured inversion method for frequency domain electromagnetic sounding data |
CN113484920B (en) * | 2021-08-17 | 2023-05-19 | 成都理工大学 | Two-dimensional structured inversion method for frequency domain electromagnetic sounding data |
Also Published As
Publication number | Publication date |
---|---|
CN104422970B (en) | 2017-03-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102495431B (en) | Data processing method for carrying out static correction on transient electromagnetic data | |
Abd Allah et al. | Three-dimensional resistivity characterization of a coastal area: Application of Grounded Electrical-Source Airborne Transient Electromagnetic (GREATEM) survey data from Kujukuri Beach, Japan | |
CN104903750A (en) | Amplitude preservation statics method of electromagnetic data | |
Fedi et al. | Multiscale analysis of potential fields by a ridge consistency criterion: the reconstruction of the Bishop basement | |
CN102998704A (en) | Geophysical exploration seismic data processing method | |
Bastani et al. | CSRMT measurements in the frequency range of 1–250 kHz to map a normal fault in the Volvi basin, Greece | |
CN104422970A (en) | Amplitude-preserving static correcting method for electromagnetic data | |
Rajagopalan | Analytic signal vs. reduction to pole: solutions for low magnetic latitudes | |
CN104656156A (en) | Magnetic reference processing method for acquiring data in three dimensions audio-frequency magnetotelluric sounding | |
CN103984011A (en) | Dynamic Q compensation shifting method | |
Gnaneshwar et al. | Very-low-frequency electromagnetic (VLF-EM) measurements in the Schirmacheroasen area, East Antarctica | |
Ramazi et al. | Contribution of geophysical inversion theory and geostatistical simulation to determine geoelectrical anomalies | |
CN104330826A (en) | A method for removing various noises under the condition of complex surface | |
Dennis et al. | Transient electromagnetic surveys for the measurement of near-surface electrical anisotropy | |
Siemon | Levelling of helicopter-borne frequency-domain electromagnetic data | |
CN103913770A (en) | Method for processing seismic data based on VSP data | |
An et al. | Auto-pick first breaks with complex raypaths for undulate surface conditions | |
CN108008456B (en) | A method of delineation mesothermal gold deposits deep three-dimensional emphasis U metallogeny Favourable Target Areas | |
CN105093318A (en) | Self-adaptive wave equation wave field continuation static correction method | |
CN102608665A (en) | Fine treatment technology of data information measured by geophysical prospecting time-domain transient electromagnetic system | |
Domenzain et al. | Joint full-waveform ground-penetrating radar and electrical resistivity inversion applied to field data acquired on the surface | |
Aizebeokhai et al. | Field evaluation of 3D geo-electrical resistivity imaging for environmental and engineering studies using parallel 2D profiles | |
CN105093324B (en) | A kind of residual static corrections | |
CN102937725A (en) | Potential field anomaly edge enhancement method based on transition zone and phase congruency | |
Emujakporue et al. | Spectral depth estimate of subsurface structures over parts of offshore Niger Delta, Nigeria |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |