CN103389514A - MT (magnetotelluric) denoising device and method - Google Patents
MT (magnetotelluric) denoising device and method Download PDFInfo
- Publication number
- CN103389514A CN103389514A CN201310303585XA CN201310303585A CN103389514A CN 103389514 A CN103389514 A CN 103389514A CN 201310303585X A CN201310303585X A CN 201310303585XA CN 201310303585 A CN201310303585 A CN 201310303585A CN 103389514 A CN103389514 A CN 103389514A
- Authority
- CN
- China
- Prior art keywords
- electric field
- field amplitude
- pole span
- denoising
- amplitude 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.)
- Granted
Links
Images
Abstract
The invention relates to an MT denoising device which comprises an electromagnetic receiver, a magnetic induction sensor and a non-polarizable electrode group, wherein the electromagnetic receiver is placed on the ground; the magnetic induction sensor is placed on the ground in a direction parallel to an x axis or a y axis and connected with the electromagnetic receiver through a cable; the non-polarizable electrode group takes a coordinate system origin as a central point, is placed symmetrically along the y axis or the x axis and grounded, and comprises a first long pole pitch electric dipole and a first short pole pitch electric dipole; the first long pole pitch electric dipole is used for measuring amplitude data of a long pole pitch electric field, and the first short pole pitch electric dipole is used for measuring amplitude data of a short pole pitch electric field; and all non-polarizable electrodes are connected with the electromagnetic receiver through wires. The device is used for calculating the correlation coefficient and the standard deviation between the amplitude data of the long pole pitch electric field and the short pole pitch electric field, when the correlation coefficient is smaller than a set correlation threshold, electric field amplitude data with largest standard deviation are removed, the correlation coefficient is calculated repeatedly until the correlation coefficient is larger than or equal to the set correlation threshold, and near field interference data are effectively removed.
Description
Technical field
The present invention relates to the resistivity prospecting field, particularly relate to a kind of data correlation analysis method that utilizes and judge field source character, and remove mt denoising device and the mt denoising method of interference noise.
Background technology
Magnetotelluric method is as a kind of important exploration geophysics method, and application is extensive, and branch is numerous.After proposing first, audio-frequency magnetotelluric magnetic method (AMT) has been proposed again the sixties in 20th century.1971, D.W.Strangbway professor and his postgraduate Myron Goldstein of University of Toronto proposed again controllable source audio-frequency magnetotelluric magnetic method (CSAMT).These methods (MT, AMT, CSAMT) are all the natural plane electromagnetic wave Electric and magnetic fields signals that is caused by far field source by observation, with the resistivity value under definitely.Wherein Observation of The Geoelectric Field research is the earth electric field part of projection at the earth's surface;on the face of the globe, electric field component is to lay pair of electrodes to receive electric field signal on specific orientation, with this to the potential difference (PD) that records on electrode and the business of the distance between spread of electrodes, as the tolerance of electric field at component on this orientation. when the field source temporal evolution, earth electric field also will change thereupon.
In recent years, along with the progress of science and technology, electromagnetic sounding method (MT, AMT, CSAMT) is no matter be the field acquisition instrument, or data are processed interpretation software, all obtained develop rapidly, the exploration amount increases year by year, and application constantly enlarges, the development of still along with national economy and industrial or agricultural, building, the electromagnetism ground unrest is day by day serious, so how effectively to avoid, suppress or eliminate interference, improve the image data quality, be the problem that at first will consider in present electromagnetic prospecting work.
The noise of mt data is mainly near field noise, can be divided into natural noise and human noise.Natural noise comes from nagneto-telluric field noise hole and thunderstorm movable (as lightning closely); Human noise refers near the electromagnetic field noise of the unsettled nonplanar wave of observation station, as power transmission line, wireless station, electric locomotive, various industrial travelling powder currents etc., causes.Particularly in industry prosperity, densely populated areas, this near-field interference ubiquity, affect more serious.
At present conventional processing Noise Method has: narrow-band filtering method, coherent method, least square method, far with reference to method, Robust method, wavelet analysis method and time domain iterative technique etc., all impacts of " flying spot " in Removing Random No and apparent resistivity curve by all means, effectively suppress noise, improved the impedance estimation quality.But, above-mentioned not fully effective, sometimes at the inhibition interference noise, also suppressed simultaneously some useful signals simultaneously.Therefore, artificial denoising remains requisite a kind of means in the electromagnetic data processing, but the human factor of the method is higher, so operator's experience has played decisive role, method can not complied with accurately.
In the telluric electromagnetic sounding method, measurement parameter is mainly magnetic field and the electric field of a pair of time dependent quadrature, its metering system can be scalar measurement (measuring Hy, Ex or Hx, Ey) or vector measurement (measuring simultaneously Hy, Hx, Ex, Ey), because electricity, magnetic signal encourage each other, and the annoyance level of interference source depends primarily on the impact size on the electromagnetic field signal to noise ratio (S/N ratio).Because described noise generally has nearly source feature, be far-field signal and mt is mainly observed, for the signal source with same intensity, far away apart from observation station, the electric field intensity component difference that produces on two different interpolar distances on observation station is less; And near apart from observation station, the electric field intensity component difference that produces on two different interpolar distances on observation station is larger; To survey district inner when signal source is positioned at, and the earth electric field component value that observes on two pole spans even can opposite in sign.Therefore shown: the frequency electromagnetic sounding data quality is impacted, and is mainly in-plant interference, and at this moment the signal to noise ratio (S/N ratio) in magnetic field always is better than electric field, and this is also that people are with the reason of road, magnetic field as local library track or library track far away.Therefore from the angle of research disturbed sources of noise, the research electric field can be more responsive than research magnetic field.
On the basis of existing correlation technique, the present invention, from a kind of mt denoising measurement mechanism,, by correlation analysis method, judges the character (near field or far field) of field source, and then the undesired signal of artificial removal's different frequent points, thereby obtain reliable far field electromagnetic measurement data.
Summary of the invention
Technical matters to be solved by this invention is to provide a kind of mt denoising device and mt denoising method, is used for solving the limitation problem that prior art is measured in far field electromagnetic.
The technical scheme that the present invention solves the problems of the technologies described above is as follows: a kind of mt denoising measurement mechanism, this mt denoising measurement mechanism are arranged in the xy coordinate system that point is set up centered by measurement point, comprising:
An electromagnetism receiver that is placed in ground, be used for reading magnetic field amplitude data and electric field amplitude data;
A magnetic induction sensor that is positioned over ground along parallel x direction of principal axis or y direction of principal axis, it is connected with described electromagnetism receiver by cable, is used for measuring the magnetic field amplitude data;
One group of point centered by the xy coordinate origin, and the non-polarizing electrode group along the y axle vertical with described magnetic induction sensor placement direction or the placement of x rotational symmetry ground connection, this non-polarizing electrode group comprises for the first long pole span eelctric dipole of measuring long pole span electric field amplitude data with for the first short pole span eelctric dipole of measuring short pole span electric field amplitude data, the pole span length of the described first long pole span eelctric dipole is greater than the pole span length of the described first short pole span eelctric dipole, and each non-polarizing electrode is connected with described electromagnetism receiver with wire respectively.
On the basis of technique scheme, the present invention can also do following improvement.
Further, described mt denoising measurement mechanism is applied to the denoising measurement of natural source mt or controllable source mt denoising measurement.
Further, described mt denoising measurement mechanism is applied in traditional Electromagnetic Array Profiling (EMAP) denoising measurement.Wherein the EMAP harvester can be to measure a magnetic field amplitude data Hy or pair of orthogonal magnetic field amplitude data Hy and Hx, arrange a plurality of not polarization measurement electrodes along x axle or y axle simultaneously, make and measure electric dipole along the end to end continuous layout of survey line, can measure simultaneously several or tens electric dipoles electric field amplitude component along survey line.
Technical scheme of the present invention also comprises a kind of mt denoising measuring method, and it adopts above-mentioned mt denoising measurement mechanism, comprising:
Step 2, the long pole span electric field amplitude data under the calculating different frequency and correlation coefficient and the standard deviation between short pole span electric field amplitude data;
Step 3,, if the correlation coefficient that calculates is removed the electric field amplitude data of standard deviation maximum, and performed step 4 less than the degree of correlation threshold value of setting, otherwise keep corresponding electric field amplitude data;
Step 4, after removing the electric field amplitude data of standard deviation maximum, repeating step 2 carries out correlation coefficient and calculates, until the correlation coefficient of calculating is more than or equal to the degree of correlation threshold value of setting.
The span of the degree of correlation threshold value of further, setting in described step 1 is 0.8 to 1.
Based on same principle, technical scheme of the present invention also comprises another mt denoising measurement mechanism, and it also is arranged in the xy coordinate system that point is set up centered by measurement point, comprising:
An electromagnetism receiver that is placed in ground, be used for reading magnetic field amplitude data and electric field amplitude data;
The a pair of magnetic induction sensor that is placed in parallel in respectively ground along horizontal x, y coordinate axis, and described magnetic induction sensor is connected with the electromagnetism receiver with cable, is respectively used to measure the magnetic field amplitude data along x direction of principal axis and y direction of principal axis;
two groups of points centered by coordinate origin, the first non-polarizing electrode group and the second non-polarizing electrode group of placing along x axle and y rotational symmetry ground connection respectively, the first non-polarizing electrode group comprises for measuring the first long pole span eelctric dipole of long pole span electric field amplitude data along the x direction of principal axis and being used for along the first short pole span eelctric dipole of the short pole span electric field amplitude data of x direction of principal axis measurement, the pole span length of the described first long pole span eelctric dipole is greater than the pole span length of the described first short pole span eelctric dipole, described the second non-polarizing electrode group comprises for measuring the second long pole span eelctric dipole of long pole span electric field amplitude data along the y direction of principal axis and being used for along the second short pole span eelctric dipole of the short pole span electric field amplitude data of y direction of principal axis measurement, and the pole span length of the described second long pole span eelctric dipole is greater than the pole span length of the described second short pole span eelctric dipole, and each non-polarizing electrode is connected with described electromagnetism receiver with wire respectively.
Further, described mt denoising measurement mechanism is applied to the denoising measurement of natural source mt or controllable source mt denoising measurement.
Further, described mt denoising measurement mechanism is applied in traditional Electromagnetic Array Profiling denoising measurement.
Technical scheme of the present invention also comprises a kind of mt denoising measuring method, and it adopts above-mentioned another kind of mt denoising measurement mechanism, comprising:
Steps A, measure magnetic field amplitude data and electric field amplitude data under different frequency, and set a degree of correlation threshold value that is used for the judgement electric field amplitude quality of data;
Step B, under different frequency, calculate respectively along x, the long pole span electric field amplitude data of y direction of principal axis measurement and correlation coefficient and the standard deviation between short pole span electric field amplitude data;
Step C, if the correlation coefficient that calculates, less than the degree of correlation threshold value of setting, is removed the electric field amplitude data of standard deviation maximum, and execution step D, otherwise keep corresponding electric field amplitude data;
Step D, after removing the electric field amplitude data of standard deviation maximum, repeating step B carries out correlation coefficient and calculates, until the correlation coefficient of calculating is more than or equal to the degree of correlation threshold value of setting.
The span of the degree of correlation threshold value of further, setting in described steps A is 0.8 to 1.
The invention has the beneficial effects as follows: mt denoising measuring system of the present invention and mt denoising measuring method can effectively judge the field source feature (near field or far field) with the electric field amplitude data of each frequency in the magnetotelluric sounding curve of frequency change, effectively remove the near-field interference data; Simultaneously, only removed the data of near field noise due to the present invention, kept merely the far field source observation data, other data de-noising disposal route, fully kept useful signal relatively, improved the confidence level of observation data.
Description of drawings
Fig. 1 is a kind of scalar crustal magnetotelluric measurement device schematic diagram that the embodiment of the present invention one proposes, and Rx is the electromagnetism receiver, and Hx is magnetic induction sensor, M1, M2, M3, M4 are non-polarizing electrode, M1-M4 measures eelctric dipole for long pole span, and M2-M3 is that short pole span is measured eelctric dipole, and O is the measurement point position;
Fig. 2 is the implementation step schematic diagram of the device of employing Fig. 1;
A kind of vector crustal magnetotelluric measurement device schematic diagram that Fig. 3 embodiment of the present invention two proposes, in figure, symbol is identical with Fig. 1, Rx is the electromagnetism receiver, and Hx, Hy are magnetic induction sensor, and M1, M2, M3, M4, N1, N2, N3, N4 are non-polarizing electrode, wherein M1, M4 and N1, N4 are that two pairs of long pole spans are measured eelctric dipole, hereinafter all with M1-M4 and N1-N4, represent,, M2, M3 and N2, N3 are that two pairs of short pole spans are measured eelctric dipole, hereinafter all with M2-M3 and N2-N3, represent, O is the measurement point position;
Fig. 4 is the implementation step schematic diagram of the device of employing Fig. 3;
Magnetic field amplitude-the frequency curve chart of Fig. 5 for measuring;
Fig. 6 for long pole span and the short utmost point according to electric field amplitude-frequency curve chart;
Fig. 7 is the short utmost point according to removing electric field amplitude-frequency curve chart after the noise frequency;
Fig. 8 is the apparent resistivity-frequency curve after denoising;
Fig. 9 is scalar controllable source crustal magnetotelluric measurement device schematic diagram, O1, O2 ... O7 is the measuring center point, M1-0, M1-1, M1-2 ... ..M1-8, M2-1 is potential electrode, D1, D2 ... D7 is that the short utmost point is according to measuring eelctric dipole, C1, C2 ... C7 measures eelctric dipole for long pole span, and S1 is the far field source signal, and S2 is the near field sources signal, rs1 be far field source to distance between potential electrode, rs2 is that near field sources is to distance between potential electrode;
Figure 10 is a plurality of arrangement magnetic field amplitude-frequency curve charts based on the described device actual measurement of Fig. 9;
Figure 11 is the short pole span electric field amplitude-frequency curve based on the described device actual measurement of Fig. 9;
Figure 12 is the long pole span electric field amplitude-frequency curve based on the described device actual measurement of Fig. 9;
Figure 13 is for carrying out electric field amplitude-frequency curve after denoising based on the described device of Fig. 9;
Figure 14 is for carrying out apparent resistivity-frequency curve after denoising based on the described device of Fig. 9.
Embodiment
Below in conjunction with accompanying drawing, principle of the present invention and feature are described, example only is used for explaining the present invention, not is used for limiting scope of the present invention.
Embodiment one
As shown in Figure 1, embodiment one has provided a kind of mt denoising measurement mechanism, it is based on the mt denoising measurement mechanism of scalar measurement pattern, be arranged in the xy coordinate system that point is set up centered by measurement point, can be applicable to the denoising measurement of natural source mt or the denoising of controllable source mt and measure.Its structure comprises:
An electromagnetism receiver Rx who is placed in ground, be used for reading magnetic field amplitude data and electric field amplitude data;
A magnetic induction sensor Hy who is positioned over ground along parallel y change in coordinate axis direction, it is connected with described electromagnetism receiver Rx by cable, is used for measuring the magnetic field amplitude data;
One group of point centered by coordinate origin, and non-polarizing electrode group M1, M2, M3, M4 along the placement of x rotational symmetry ground connection, this non-polarizing electrode group comprises for the first long pole span eelctric dipole M1-M4 that measures long pole span electric field amplitude data with for the first short pole span eelctric dipole M2-M3 that measures short pole span electric field amplitude data, the pole span length of the described first long pole span eelctric dipole M1-M4 is greater than the pole span length of the described first short pole span eelctric dipole M2-M3, and each non-polarizing electrode is connected with described electromagnetism receiver Rx with wire respectively.
For this mt denoising measurement mechanism, also magnetic induction sensor can be positioned over ground along parallel x change in coordinate axis direction, accordingly the non-polarizing electrode group being placed along y rotational symmetry ground connection, is namely that above-mentioned mt denoising measurement mechanism has been rotated 90 degree along central point.
As shown in Figure 2, the present embodiment gives a kind of method that adopts above-mentioned mt denoising measurement mechanism to carry out mt denoising measurement, comprising:
Step 2, the long pole span electric field amplitude data under the calculating different frequency and correlation coefficient and the standard deviation between short pole span electric field amplitude data;
Step 3,, if the correlation coefficient that calculates is removed the electric field amplitude data of standard deviation maximum, and performed step 4 less than the degree of correlation threshold value of setting, otherwise keep corresponding electric field amplitude data;
Step 4, after removing the electric field amplitude data of standard deviation maximum, repeating step 2 carries out correlation coefficient and calculates, until the correlation coefficient of calculating is more than or equal to the degree of correlation threshold value of setting.
Except above-mentioned measuring method, the mt denoising measurement mechanism in the present embodiment is equally applicable to traditional Electromagnetic Array Profiling denoising and measures.
Embodiment two
As shown in Figure 3, embodiment two has provided a kind of mt denoising measurement mechanism of vector measurement pattern, be arranged in equally in the xy coordinate system that point is set up centered by measurement point, can be applicable to the denoising measurement of natural source mt or the denoising of controllable source mt and measure.The device of embodiment two is based on that mt denoising measurement mechanism that embodiment one is used for the scalar measurement pattern builds, and it needs to increase the measurement mechanism under a cover scalar measurement pattern, and makes it take measuring center o as the axle 90-degree rotation.Concrete structure is:
An electromagnetism receiver Rx who is placed in ground, be used for reading magnetic field amplitude data and electric field amplitude data;
A pair of magnetic induction sensor Hx, the Hy that is placed in parallel in respectively ground along horizontal x, y coordinate axis, and described magnetic induction sensor Hx, Hy be connected with electromagnetism receiver Rx with cable, is respectively used to measure the magnetic field amplitude data along x direction of principal axis and y direction of principal axis;
two groups of points centered by coordinate origin, the the first non-polarizing electrode group M1 that places along x axle and y rotational symmetry ground connection respectively, M2, M3, M4 and the second non-polarizing electrode group N1, N2, N3, N4, the first non-polarizing electrode group comprises for measuring the first long pole span eelctric dipole M1-M4 of long pole span electric field amplitude data along the x direction of principal axis and being used for along the first short pole span eelctric dipole M2-M3 of the short pole span electric field amplitude data of x direction of principal axis measurement, the pole span length of the described first long pole span eelctric dipole M1-M4 is greater than the pole span length of the described first short pole span eelctric dipole M2-M3, described the second non-polarizing electrode group N1, N2, N3, N4 comprises for measure the second long pole span eelctric dipole N1 of long pole span electric field amplitude data along the y direction of principal axis, N4 and be used for measuring along the y direction of principal axis the second short pole span eelctric dipole N2-N3 of short pole span electric field amplitude data, and the pole span length of the described second long pole span eelctric dipole N1-N4 is greater than the pole span length of the described second short pole span eelctric dipole N2-N3, and each non-polarizing electrode is connected with described electromagnetism receiver Rx with wire respectively.
As shown in Figure 4, adopt above-mentioned vector measurement pattern to carry out the method that the mt denoising is measured, comprising:
Steps A, measure magnetic field amplitude data and electric field amplitude data under different frequency, and set a degree of correlation threshold value that is used for the judgement electric field amplitude quality of data;
Step B, under different frequency, calculate respectively along x, the long pole span electric field amplitude data of y direction of principal axis measurement and correlation coefficient and the standard deviation between short pole span electric field amplitude data;
Step C, if the correlation coefficient that calculates, less than the degree of correlation threshold value of setting, is removed the electric field amplitude data of standard deviation maximum, and execution step D, otherwise keep corresponding electric field amplitude data;
Step D, after removing the electric field amplitude data of standard deviation maximum, repeating step B carries out correlation coefficient and calculates, until the correlation coefficient of calculating is more than or equal to the degree of correlation threshold value of setting.
Except above-mentioned measuring method, the mt denoising measurement mechanism in the present embodiment is equally applicable to traditional array measuring method.
For above-mentioned two embodiment, the span of the degree of correlation threshold value of wherein setting is 0.8 to 1, is taken as 0.99 in embodiment one and embodiment two.After getting and determining degree of correlation threshold value, can calculate correlation coefficient by following correlation analysis computing formula:
In formula, R is correlation coefficient, and x is short pole span electric field amplitude Exsi (or Eysi), and y is long pole span electric field amplitude Exli (or Eyli), and i=1,2,3,4n, represent different frequency, and N is for participating in the frequency number of calculating.
More than or equal to 0.99 the time, show that the signal source of observation data is mainly far-field signal as the correlation coefficient R that calculates, the data observation quality is good; As the correlation coefficient R that calculates less than 0.99 the time, the signal source that shows observation data has near-field signals, there is near field noise in observation data, utilize the computing method of standard deviation, calculate successively the standard deviation of a plurality of pole span measurement data of respective frequencies, remove the electromagnetic data of a frequency measurement of standard deviation maximum., until observation data all met R more than or equal to 0.99 o'clock, just reached the purpose of artificial elimination near field noise.
Have particular application as example with one, on known ground, carry out crustal magnetotelluric measurement, obtain simultaneously under the scalar observation mode with the magnetic field amplitude data Hy of frequency change, the short utmost point according to electric field amplitude data Ex
M1-N1, long pole span electric field amplitude data Ex
M2-N2
Due to the angle from the research disturbed sources of noise, magnetic field amplitude is eager to excel in whatever one does with respect to the electric field amplitude antijamming capability, and magnetic field amplitude is relatively stable in certain area, and actual measurement obtains 2
13, 2/3 * 2
13, 2
12, 2/3 * 2
12, 2
11, 2/3 * 2
11, 2
102
1, 2/3 * 2
1, 2
0Totally 27 measure frequencies magnetic field amplitude-frequency data curve, as shown in Figure 5, line smoothing, stable, show that the be interfered impact of noise of magnetic field amplitude Hy is less.
Fig. 6 is the long and short pole span actual measurement of actual measurement electric field amplitude-frequency data curve, and curve is 2
9, 2
5, 2/3 * 2
1Hertz (Hz) three frequencies exist zigzag to change, and the relatively short utmost point of long pole span electric field amplitude changes obviously and strengthens according to electric field amplitude.The short utmost point is carried out degree of correlation analytical calculation according to the electric field amplitude measurement data of 27 frequencies with the electric field amplitude measurement data of 27 frequencies of long pole span, and result of calculation correlation coefficient R is 0.968, shows and has the near field sources noise.
Calculate the standard deviation of the short utmost point according to electric field amplitude measurement data with the electric field amplitude measurement data of 27 frequencies of long pole span of 27 frequencies, wherein 2
9, 2
5, 2/3 * 2
1The standard deviation at hertz (Hz) frequency place is maximum, is respectively-26.54 ,-152.12,132.34.Remove 2
5The electric field amplitude observation data of hertz (Hz) frequency, adopt formula (1) to recalculate the short utmost point according to the electric field amplitude measurement data of 26 frequencies and the correlation coefficient of the electric field amplitude measurement data of 26 frequencies of long pole span, result of calculation correlation coefficient R is 0.978, show that noise obviously reduces, but still have the near field sources interference noise.
Remove again 2/3 * 2
1The electric field amplitude observation data of hertz (Hz) frequency, recalculate the short utmost point according to the electric field amplitude measurement data of 25 frequencies and the correlation coefficient of the electric field amplitude measurement data of 25 frequencies of long pole span, result of calculation correlation coefficient R is 0.999, the signal source that shows observation data is mainly far-field signal, the data observation quality is good, and wherein 2
9Hertz (Hz) electric field amplitude observation data is that far-field signal causes, should be the objective electric field response of geologic body.
Fig. 7 is for removing 2/3 * 2
1Hertz (Hz) and 2
5Hertz (Hz) short utmost point is according to electric field amplitude-frequency observation data curve, and curve is substantially round and smooth, is mainly far-field signal electric field measurement result.
Fig. 8 is according to the magnetic field amplitude data of the short utmost point after denoising according to electric field amplitude data and actual observation, according to apparent resistivity ρ
sComputing formula:
Wherein ω, μ, E
x, H
yBe respectively magnetic permeability, electric field amplitude, magnetic field amplitude in angular frequency, vacuum, the apparent resistivity that calculating obtains-frequency electromagnetic sounding curve map, curve shows along with frequency from 2
13Hertz is to 2
0The variation of hertz, investigation depth increases gradually, and underground apparent resistivity has the trend that raises gradually, can carry out the resistivity inverting on this basis, thereby obtain underground true resistance rate value.
Mt denoising measurement mechanism and the measuring method of above-mentioned two embodiment, be easy to expand to controllable source crustal magnetotelluric measurement field.As shown in Figure 9, the controllable source crustal magnetotelluric measurement of carrying out at Inner Mongol gold mine areal coverage for practical application, the scalar measurement device that relates to an arrangement, measure and adopt a magnetic-field measurement passage Hy and 7 electric field measurement passage D1, D2 ... D7, totally 8 potential electrode (are M1-0, M1-1, M1-2 ... ..M1-8), S1 is the artificial field Electromagnetic Sources that is placed in far field, r
s1To the distance of measuring eelctric dipole, S2 is the near field noise source that possible exist, r for far field source
s2For the distance of near field sources to the measurement eelctric dipole.
Measure in order to complete denoising, definition D1, D2 ... D7 is for respectively with O1, O2 ... O7 is that the short utmost point of measuring center point is according to electric field measurement, use last electric field measurement passage D0 of previous arrangement and first electric field measurement passage D8 of a rear arrangement, make the arithmetic mean of D0, D1, three short pole span electric field measurements of D2, the long pole span electric field measurement C1 of formation take O1 as measuring center point, by that analogy, C2, C3, C4, C5, C6, C7 are respectively the long pole span electric field measurement of corresponding measuring center point O2, O3, O4, O5, O6, O7.
As can be known, reality is completed 7 arrangements controllable source crustal magnetotelluric measurement of totally 49 measurement points altogether, and Figure 10 is that 7 magnetic field amplitude-frequency data are measured curve, and curve is round and smooth, and consistance is stronger, shows that the suffered noise jamming of magnetic field amplitude is less.
Figure 11 is the short pole span electric field amplitude-frequency curve of actual measurement, Figure 12 is the long pole span electric field amplitude-frequency curve of actual measurement that obtains after the arithmetic mean value is calculated, short pole span electric field amplitude-frequency curve and long pole span electric field amplitude-frequency curve are compared, curvilinear motion is not quite identical, may be subject to the interference near field noise.
Then according to the short pole span of aforementioned definitions and long pole span electric field amplitude measurement data, respectively the long and short pole span electric field amplitude measurement data of each depth measurement point is carried out related coefficient analysis and standard variance calculating, remove and to affect correlation coefficient less than 0.99 frequency electric field measurement data, the short pole span electric field amplitude-frequency curve (or data) after the noise of marching into the arena is removed in acquisition finally, as shown in figure 13.
, according to march into the arena short pole span electric field amplitude-frequency curve (or data) after noise of removal, obtain as calculated the apparent resistivity-frequency curve (data) of 49 depth measurement points, as shown in figure 14.
Finally 49 apparent resistivity-frequency curves (data) are carried out one-dimensional inversion, subsurface resistivity one-dimensional inversion sectional drawing can be obtained, spatial position or the form of changing condition with the Hidden Granite Body that may exist of latent mining area underground structure can be clearly reflected.
The foregoing is only preferred embodiment of the present invention, in order to limit the present invention, within the spirit and principles in the present invention not all, any modification of doing, be equal to replacement, improvement etc., within all should being included in protection scope of the present invention.
Claims (10)
1. a mt denoising measurement mechanism, is characterized in that, this mt denoising measurement mechanism is arranged in the xy coordinate system that point is set up centered by measurement point, comprising:
An electromagnetism receiver that is placed in ground, be used for reading magnetic field amplitude data and electric field amplitude data;
A magnetic induction sensor that is positioned over ground along parallel x direction of principal axis or y direction of principal axis, it is connected with described electromagnetism receiver by cable, is used for measuring the magnetic field amplitude data;
One group of point centered by the xy coordinate origin, and the non-polarizing electrode group along the y axle vertical with described magnetic induction sensor placement direction or the placement of x rotational symmetry ground connection, this non-polarizing electrode group comprises for the first long pole span eelctric dipole of measuring long pole span electric field amplitude data with for the first short pole span eelctric dipole of measuring short pole span electric field amplitude data, the pole span length of the described first long pole span eelctric dipole is greater than the pole span length of the described first short pole span eelctric dipole, and each non-polarizing electrode is connected with described electromagnetism receiver with wire respectively.
2. mt denoising measurement mechanism according to claim 1, is characterized in that, described mt denoising measurement mechanism is applied in the denoising measurement of natural source mt or controllable source mt denoising measurement.
3. mt denoising measurement mechanism according to claim 1, is characterized in that, described mt denoising measurement mechanism is applied in traditional Electromagnetic Array Profiling denoising measurement.
4. mt denoising measuring method, it adopts arbitrary described mt denoising measurement mechanism in claims 1 to 3, it is characterized in that, comprising:
Step 1, measure magnetic field amplitude data and electric field amplitude data under different frequency, and set a degree of correlation threshold value that is used for the judgement electric field amplitude quality of data;
Step 2, the long pole span electric field amplitude data under the calculating different frequency and correlation coefficient and the standard deviation between short pole span electric field amplitude data;
Step 3,, if the correlation coefficient that calculates is removed the electric field amplitude data of standard deviation maximum, and performed step 4 less than the degree of correlation threshold value of setting, otherwise keep corresponding electric field amplitude data;
Step 4, after removing the electric field amplitude data of standard deviation maximum, repeating step 2 carries out correlation coefficient and calculates, until the correlation coefficient of calculating is more than or equal to the degree of correlation threshold value of setting.
5. mt denoising measuring method according to claim 4, is characterized in that, the span of the degree of correlation threshold value of setting in described step 1 is 0.8 to 1.
6. a mt denoising measurement mechanism, is characterized in that, this mt denoising measurement mechanism is arranged in the xy coordinate system that point is set up centered by measurement point, comprising:
An electromagnetism receiver (Rx) that is placed in ground, be used for reading magnetic field amplitude data and electric field amplitude data;
The a pair of magnetic induction sensor (Hx, Hy) that is placed in parallel in respectively ground along horizontal x, y coordinate axis, and described magnetic induction sensor (Hx, Hy) is connected with electromagnetism receiver (Rx) with cable, is respectively used to measure the magnetic field amplitude data along x direction of principal axis and y direction of principal axis;
two groups of points centered by coordinate origin, the first non-polarizing electrode group (M1 that places along x axle and y rotational symmetry ground connection respectively, M2, M3, M4) and the second non-polarizing electrode group (N1, N2, N3, N4), the first non-polarizing electrode group (M1, M2, M3, M4) comprise for measure the first long pole span eelctric dipole (M1 of long pole span electric field amplitude data along the x direction of principal axis, M4) and be used for measuring along the x direction of principal axis the first short pole span eelctric dipole (M2 of short pole span electric field amplitude data, M3), the described first long pole span eelctric dipole (M1, M4) pole span length is greater than the described first short pole span eelctric dipole (M2, M3) pole span length, described the second non-polarizing electrode group (N1, N2, N3, N4) comprise for measure the second long pole span eelctric dipole (N1 of long pole span electric field amplitude data along the y direction of principal axis, N4) and be used for measuring along the y direction of principal axis the second short pole span eelctric dipole (N2 of short pole span electric field amplitude data, and the described second long pole span eelctric dipole (N1 N3), N4) pole span length is greater than the described second short pole span eelctric dipole (N2, N3) pole span length, and each non-polarizing electrode is connected with described electromagnetism receiver (Rx) with wire respectively.
7. mt denoising measurement mechanism according to claim 6, is characterized in that, described mt denoising measurement mechanism is applied to the denoising measurement of natural source mt or the denoising of controllable source mt is measured.
8. mt denoising measurement mechanism according to claim 6, is characterized in that, described mt denoising measurement mechanism is applied in traditional Electromagnetic Array Profiling denoising measurement.
9. mt denoising measuring method, it adopts arbitrary described mt denoising measurement mechanism in claim 6 to 8, it is characterized in that, comprising:
Steps A, measure magnetic field amplitude data and electric field amplitude data under different frequency, and set a degree of correlation threshold value that is used for the judgement electric field amplitude quality of data;
Step B, under different frequency, calculate respectively along x, the long pole span electric field amplitude data of y direction of principal axis measurement and correlation coefficient and the standard deviation between short pole span electric field amplitude data;
Step C, if the correlation coefficient that calculates, less than the degree of correlation threshold value of setting, is removed the electric field amplitude data of standard deviation maximum, and execution step D, otherwise keep corresponding electric field amplitude data;
Step D, after removing the electric field amplitude data of standard deviation maximum, repeating step B carries out correlation coefficient and calculates, until the correlation coefficient of calculating is more than or equal to the degree of correlation threshold value of setting.
10. mt denoising measuring method according to claim 9, is characterized in that, the span of the degree of correlation threshold value of setting in described steps A is 0.8 to 1.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310303585.XA CN103389514B (en) | 2013-07-18 | 2013-07-18 | A kind of mt denoising device and mt denoising method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310303585.XA CN103389514B (en) | 2013-07-18 | 2013-07-18 | A kind of mt denoising device and mt denoising method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN103389514A true CN103389514A (en) | 2013-11-13 |
| CN103389514B CN103389514B (en) | 2016-04-20 |
Family
ID=49533837
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201310303585.XA Active CN103389514B (en) | 2013-07-18 | 2013-07-18 | A kind of mt denoising device and mt denoising method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN103389514B (en) |
Cited By (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105629317A (en) * | 2016-04-08 | 2016-06-01 | 中国矿业大学(北京) | Magnetotelluric noise suppressing method based on intersite transfer function |
| RU2594641C1 (en) * | 2015-05-25 | 2016-08-20 | Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Омский государственный технический университет" | System of sensors for magnetotelluric sounding of the earth |
| CN106324433A (en) * | 2016-08-31 | 2017-01-11 | 国网山东省电力公司菏泽供电公司 | Quick maintenance method for power distribution line |
| CN107092033A (en) * | 2017-04-19 | 2017-08-25 | 中南大学 | It is a kind of to be used to regulate and control the method for measurement signal intensity |
| CN107329183A (en) * | 2017-07-14 | 2017-11-07 | 中国地质科学院地球物理地球化学勘查研究所 | A kind of controlled-source audiomagnetotellurics sounding collecting method and device |
| CN107657242A (en) * | 2017-10-10 | 2018-02-02 | 湖南师范大学 | A kind of mt noise identification and separation method |
| CN110488365A (en) * | 2019-09-11 | 2019-11-22 | 湖南五维地质科技有限公司 | A kind of multipolarization telluric electromagnetic sounding method |
| CN110596763A (en) * | 2019-08-23 | 2019-12-20 | 南方科技大学 | Three-dimensional acquisition method and device for magnetotelluric data and terminal equipment |
| CN111190234A (en) * | 2020-01-13 | 2020-05-22 | 山东大学 | Noise observation method and device for artificial electrical source frequency domain electromagnetic method |
| CN111983701A (en) * | 2020-09-01 | 2020-11-24 | 湖南五维地质科技有限公司 | Transverse geodetic electromagnetic wave depth measuring method and depth measuring device thereof |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4392109A (en) * | 1980-11-14 | 1983-07-05 | Mobil Oil Corporation | Static correction for magnetotelluric data |
| CN86103551A (en) * | 1986-04-25 | 1987-11-04 | 得克萨斯系统大学董事会 | Electromagnetic array profiling survey method |
| WO1997041457A1 (en) * | 1996-04-26 | 1997-11-06 | Anthony Charles Leonid Fox | Satellite synchronized 3-d magnetotelluric system |
| CN101071179A (en) * | 2006-05-11 | 2007-11-14 | 中国石油集团东方地球物理勘探有限责任公司 | Magnetotelluric impedance measuring method |
| CN102879828A (en) * | 2012-10-11 | 2013-01-16 | 甘肃铁道综合工程勘察院有限公司 | Audio magnetotelluric optimized collection method |
-
2013
- 2013-07-18 CN CN201310303585.XA patent/CN103389514B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4392109A (en) * | 1980-11-14 | 1983-07-05 | Mobil Oil Corporation | Static correction for magnetotelluric data |
| CN86103551A (en) * | 1986-04-25 | 1987-11-04 | 得克萨斯系统大学董事会 | Electromagnetic array profiling survey method |
| WO1997041457A1 (en) * | 1996-04-26 | 1997-11-06 | Anthony Charles Leonid Fox | Satellite synchronized 3-d magnetotelluric system |
| CN101071179A (en) * | 2006-05-11 | 2007-11-14 | 中国石油集团东方地球物理勘探有限责任公司 | Magnetotelluric impedance measuring method |
| CN102879828A (en) * | 2012-10-11 | 2013-01-16 | 甘肃铁道综合工程勘察院有限公司 | Audio magnetotelluric optimized collection method |
Non-Patent Citations (5)
| Title |
|---|
| T.D.GAMBLE: "Magnetotellurics with a remote magnetic reference", 《GEOPHYSICS》, vol. 44, no. 1, 31 January 1979 (1979-01-31), pages 53 - 68 * |
| 叶青: "地电场变化的基本要素研究及物理解释", 《中国优秀博硕士论文全文数据库（硕士） 基础科学辑》, no. 12, 15 December 2006 (2006-12-15) * |
| 杨生等: "远参考法大地电磁测深法应用研究", 《物探与化探》, vol. 26, no. 1, 28 February 2002 (2002-02-28) * |
| 汤井田: "庐枞矿集区大地电磁测深强噪声的影响规律", 《地球物理学报》, vol. 55, no. 12, 31 December 2012 (2012-12-31), pages 4147 - 4159 * |
| 郭志强等: "可控源音频大地电磁在铁路工程勘察中的应用", 《铁道建筑技术》, no. 1, 31 May 2008 (2008-05-31) * |
Cited By (16)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU2594641C1 (en) * | 2015-05-25 | 2016-08-20 | Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Омский государственный технический университет" | System of sensors for magnetotelluric sounding of the earth |
| CN105629317A (en) * | 2016-04-08 | 2016-06-01 | 中国矿业大学(北京) | Magnetotelluric noise suppressing method based on intersite transfer function |
| CN106324433B (en) * | 2016-08-31 | 2019-01-01 | 国网山东省电力公司菏泽供电公司 | A kind of distribution line Quick overhaul method |
| CN106324433A (en) * | 2016-08-31 | 2017-01-11 | 国网山东省电力公司菏泽供电公司 | Quick maintenance method for power distribution line |
| CN107092033A (en) * | 2017-04-19 | 2017-08-25 | 中南大学 | It is a kind of to be used to regulate and control the method for measurement signal intensity |
| CN107092033B (en) * | 2017-04-19 | 2019-06-11 | 中南大学 | A method of for regulating and controlling measuring signal intensity |
| CN107329183B (en) * | 2017-07-14 | 2019-03-22 | 中国地质科学院地球物理地球化学勘查研究所 | A kind of controlled-source audiomagnetotellurics sounding collecting method and device |
| CN107329183A (en) * | 2017-07-14 | 2017-11-07 | 中国地质科学院地球物理地球化学勘查研究所 | A kind of controlled-source audiomagnetotellurics sounding collecting method and device |
| CN107657242B (en) * | 2017-10-10 | 2018-08-21 | 湖南师范大学 | A kind of identification of magnetotelluric noise and separation method |
| CN107657242A (en) * | 2017-10-10 | 2018-02-02 | 湖南师范大学 | A kind of mt noise identification and separation method |
| CN110596763A (en) * | 2019-08-23 | 2019-12-20 | 南方科技大学 | Three-dimensional acquisition method and device for magnetotelluric data and terminal equipment |
| 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 |
| CN111190234A (en) * | 2020-01-13 | 2020-05-22 | 山东大学 | Noise observation method and device for artificial electrical source frequency domain electromagnetic method |
| CN111983701A (en) * | 2020-09-01 | 2020-11-24 | 湖南五维地质科技有限公司 | Transverse geodetic electromagnetic wave depth measuring method and depth measuring device thereof |
| CN111983701B (en) * | 2020-09-01 | 2023-10-10 | 湖南五维地质科技有限公司 | Transverse magnetotelluric wave sounding method and sounding device thereof |
Also Published As
| Publication number | Publication date |
|---|---|
| CN103389514B (en) | 2016-04-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103389514B (en) | A kind of mt denoising device and mt denoising method | |
| Jiracek | Near-surface and topographic distortions in electromagnetic induction | |
| CN103869371B (en) | Manual field source frequency domain full-gradient electromagnetic measuring method | |
| CN103913776B (en) | The detection method of spherical boulder in shield-tunneling construction | |
| US20060186887A1 (en) | Method for identifying subsurface features from marine transient controlled source electromagnetic surveys | |
| Di et al. | Geophysical survey over molybdenum mines using the newly developed M-TEM system | |
| CN104656156A (en) | Magnetic reference processing method for acquiring data in three dimensions audio-frequency magnetotelluric sounding | |
| US20110001482A1 (en) | Electromagnetic Survey Using Naturally Occurring Electromagnetic Fields as a Source | |
| CN102707323A (en) | Controllable source audio-frequency magnetic field sounding method for geological exploration | |
| CN104375195A (en) | Time-frequency electromagnetic multi-source multi-component three-dimensional joint inversion method | |
| CN103837898A (en) | High-density electrical method near-end dipole electromagnetic sounding method | |
| CN103995301A (en) | Method and device for evaluating total organic carbon content in shale gas reservoir | |
| CN105911595A (en) | Method and apparatus for obtaining controllable source audio-frequency magnetotelluric (CSAMT) apparent phase information | |
| US8447521B2 (en) | Subsurface electrical conductivity contrast reconstruction using narrow-band electromagnetic data and fullwave imaging | |
| Kuras et al. | Capacitive resistivity imaging with towed arrays | |
| EP2230534A1 (en) | Method for determining resistivity anistropy from earth electromagnetic responses | |
| CN105445805A (en) | Time-space array differential electromagnetic prospecting method | |
| CN105204073A (en) | Tensor apparent conductivity measurement method | |
| CA2834518A1 (en) | Method and system for surveying a distribution of chargeability in a volume of earth | |
| CN104749643A (en) | CSAMT (controlled source audio frequency magneto telluric) measuring method and device used for electromagnetic interference areas | |
| WO2021098667A1 (en) | Full-coverage anti-interference controlled-source electromagnetic detection method | |
| CN103135140A (en) | Computing method of center loop line transverse electric and magnetic field (TEM) whole period true resistivity without fringe effect | |
| Guillemoteau et al. | Evaluation of the 3-D Multichannel Deconvolution Method for the Case of Low S/N Inphase Data Collected With Loop–Loop FD-EMI Sensors | |
| CN113406707A (en) | Magnetotelluric multi-scale and multi-time-period detection method | |
| CN107329183B (en) | A kind of controlled-source audiomagnetotellurics sounding collecting method and device |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| CB03 | Change of inventor or designer information |
Inventor after: Xing Baoshan Inventor after: Yu Changming Inventor before: Yu Changming Inventor before: Xing Baoshan |
|
| COR | Change of bibliographic data | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20200522 Address after: 012000 science and technology research and development base of Chahar Industrial Park, Wulanchabu City, Inner Mongolia Autonomous Region Patentee after: Ulanqab Industrial Technology Research Institute Co., Ltd Address before: 100029 Beijing city Chaoyang District Beitucheng West Road No. 19 Patentee before: INSTITUTE OF GEOLOGY AND GEOPHYSICS, CHINESE ACADEMY OF SCIENCES |
|
| TR01 | Transfer of patent right |