CN107329183A - A kind of controlled-source audiomagnetotellurics sounding collecting method and device - Google Patents
A kind of controlled-source audiomagnetotellurics sounding collecting method and device Download PDFInfo
- Publication number
- CN107329183A CN107329183A CN201710574423.8A CN201710574423A CN107329183A CN 107329183 A CN107329183 A CN 107329183A CN 201710574423 A CN201710574423 A CN 201710574423A CN 107329183 A CN107329183 A CN 107329183A
- Authority
- CN
- China
- Prior art keywords
- field signal
- magnetic field
- electric field
- frequency
- value
- 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
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
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention provides a kind of controlled-source audiomagnetotellurics sounding collecting method and device, this method includes：Receiver/transmitter transmitting propagated through earth medium after electric field signal and magnetic field signal；The frequency of frequency or magnetic field signal based on electric field signal, and fundamental frequency frequency, obtain sample rate and sampling length；The electric field signal is sampled using the sample rate and using length, electric field information is obtained, the magnetic field signal is sampled using the sample rate and using length, Magnetic Field is obtained；Based on the electric field information and Magnetic Field, the electrical structure of the earth medium is calculated, wherein, the fundamental wave is the artificial electric power interference signal of influence electric field signal and magnetic field signal.Invention additionally discloses corresponding controlled-source audiomagnetotellurics sounding data acquisition device.
Description
Technical field
The present invention relates to geology detecting technical field, in particular to a kind of controlled-source audiomagnetotellurics sounding data
Acquisition method and device.
Background technology
Controlled-source audiomagnetotellurics sounding (Controlled Source Audio-frequency
Magnetotellurics, abbreviation CSAMT) it is a kind of artificial source frequency domain electromagnetic sounding method in geophysical exploration, lead to
Development CSAMT measurement work is crossed, the electrical structure of underground medium can be obtained, is mainly used in concealed mineral resource, underground water money
The detection in source and rift structure, earth formation etc..But in CSAMT measurements, the not only transmission signal of someone's industry control, simultaneously
Also include various noises, for example, based on the random signal of Natural electromagnetic field, receiver background noise, 50Hz and its harmonic wave
Power frequency electromagnet is disturbed.Particularly, when the power frequency electromagnet interference based on the 50Hz and its harmonic wave is very big, the accuracy meeting of data is obtained
It is greatly lowered, it is impossible to objectively respond the electrical property feature of underground medium.
The content of the invention
In view of this, it is an object of the invention to provide a kind of controlled-source audiomagnetotellurics sounding collecting method and
Device, the inaccurate problem of data for solving to cause in fundamental wave and its larger harmonic wave interference.
In a first aspect, embodiments of the invention provide a kind of controlled-source audiomagnetotellurics sounding collecting method, bag
Include：
Receiver/transmitter transmitting propagated through earth medium after electric field signal and magnetic field signal；
The frequency of frequency or magnetic field signal based on electric field signal, and fundamental frequency frequency, obtain sample rate and adopting
Sample length；
The electric field signal is sampled using the sample rate and using length, electric field information is obtained, utilizes the sampling
Rate and use length are sampled to the magnetic field signal, obtain Magnetic Field；
Based on the electric field information and Magnetic Field, the electrical structure of the earth medium is calculated, wherein, the fundamental wave is
Influence the artificial electric power interference signal of electric field signal and magnetic field signal.
Alternatively, the calculating obtains sample rate and sampling length, including：
The frequency of electric field signal or the frequency and the product of predetermined threshold value of magnetic field signal are calculated, sample rate is obtained, it is described pre-
If threshold value is the value more than or equal to the first numerical value；
The ratio and the product of second value of sample rate and the frequency of fundamental wave are calculated, sampling length, second number is obtained
Value is by calculating the ratio of the frequency of fundamental wave and the frequency of electric field signal or the frequency of magnetic field signal, being obtained with the product of third value
Arrive, the second value and the third value are positive integer.
Alternatively, it is described that the electric field signal is sampled using the sample rate and using length, electric field information is obtained,
The magnetic field signal is sampled using the sample rate and using length, Magnetic Field is obtained, including：
Conversion process is carried out to the electric field signal after sampling using discrete Fourier transform；
Conversion process is carried out to the magnetic field signal after sampling using discrete Fourier transform.
Alternatively, the electrical structure at least includes resistivity, described to be based on the electric field information and Magnetic Field, calculates
The electrical structure of the earth medium, including：
Calculate electric field information and Magnetic Field ratio square, the cycle of the cycle of electric field signal or magnetic field signal and the
The product of four numerical value, obtains the resistivity of earth medium.
Second aspect, a kind of controlled-source audiomagnetotellurics sounding data acquisition device, including：
Receiving unit, for receiver/transmitter launch propagated through earth medium after electric field signal and magnetic field signal；
First processing units, the frequency for frequency or the frequency of magnetic field signal based on electric field signal, and fundamental frequency
Rate, obtains sample rate and sampling length；
Second processing unit, for being sampled using the sample rate and using length to the electric field signal, obtains electricity
Field information, is sampled to the magnetic field signal using the sample rate and using length, obtains Magnetic Field；
3rd processing unit, for based on the electric field information and Magnetic Field, calculating the electrical knot of the earth medium
Structure, wherein, the fundamental wave is the artificial electric power interference signal of influence electric field signal and magnetic field signal.
Alternatively, the first processing units specifically for：
The frequency of electric field signal or the frequency and the product of predetermined threshold value of magnetic field signal are calculated, sample rate is obtained, it is described pre-
If threshold value is the value more than or equal to the first numerical value；
The ratio and the product of second value of sample rate and the frequency of fundamental wave are calculated, sampling length, second number is obtained
Value is by calculating the ratio of the frequency of fundamental wave and the frequency of electric field signal or the frequency of magnetic field signal, being obtained with the product of third value
Arrive, the second value and the third value are positive integer.
Alternatively, the second processing unit specifically for：
Conversion process is carried out to the electric field signal after sampling using discrete Fourier transform；
Conversion process is carried out to the magnetic field signal after sampling using discrete Fourier transform.
Alternatively, the electrical structure at least include resistivity, the 3rd processing unit specifically for：
Calculate electric field information and Magnetic Field ratio square, the cycle of the cycle of electric field signal or magnetic field signal and the
The product of four numerical value, obtains the resistivity of earth medium.
Technique according to the invention scheme, in the frequency and the sample rate of receiver of optimization design emitter, is effectively eliminated
The influence of fundamental wave and its harmonic wave to gathered data, improves acquisition precision, reduces work degree.
To enable the above objects, features and advantages of the present invention to become apparent, preferred embodiment cited below particularly, and coordinate
Appended accompanying drawing, is described in detail below.
Brief description of the drawings
In order to illustrate the technical solution of the embodiments of the present invention more clearly, below will be attached to what is used required in embodiment
Figure is briefly described, it will be appreciated that the following drawings illustrate only certain embodiments of the present invention, therefore is not construed as pair
The restriction of scope, for those of ordinary skill in the art, on the premise of not paying creative work, can also be according to this
A little accompanying drawings obtain other related accompanying drawings.
Fig. 1 is a kind of flow of controlled-source audiomagnetotellurics sounding collecting method provided in an embodiment of the present invention
Figure；
Fig. 2 is a kind of schematic diagram handled electric field signal and magnetic field signal provided in an embodiment of the present invention；
Fig. 3 is that a kind of structure of controlled-source audiomagnetotellurics sounding data acquisition device provided in an embodiment of the present invention is shown
It is intended to.
Embodiment
To make the purpose, technical scheme and advantage of the embodiment of the present invention clearer, below in conjunction with the embodiment of the present invention
Middle accompanying drawing, the technical scheme in the embodiment of the present invention is clearly and completely described, it is clear that described embodiment is only
It is a part of embodiment of the invention, rather than whole embodiments.The present invention being generally described and illustrated herein in the accompanying drawings is real
Applying the component of example can be arranged and be designed with a variety of configurations.Therefore, it is of the invention to what is provided in the accompanying drawings below
The detailed description of embodiment is not intended to limit the scope of claimed invention, but is merely representative of the selected reality of the present invention
Apply example.Based on embodiments of the invention, the institute that those skilled in the art are obtained on the premise of creative work is not made
There is other embodiment, belong to the scope of protection of the invention.
Fig. 1 is according to a kind of knot of controlled-source audiomagnetotellurics sounding collecting method provided in an embodiment of the present invention
Structure schematic diagram.As shown in figure 1, this method starts from step S110.
In step s 110, receiver/transmitter transmitting propagated through earth medium after electric field signal and magnetic field signal.
In the step s 120, based on electric field signal frequency or the frequency of magnetic field signal, and fundamental frequency frequency, obtain
To sample rate and sampling length.
In step s 130, the electric field signal is sampled using the sample rate and using length, obtains electric field letter
Breath, is sampled to the magnetic field signal using the sample rate and using length, obtains Magnetic Field.
In step S140, based on the electric field information and Magnetic Field, the electrical structure of the earth medium is calculated, its
In, the fundamental wave is the artificial electric power interference signal of influence electric field signal and magnetic field signal.
The present invention carries out data acquisition based on geophysics electromagnetic method, enters preferably by controlled-source audiomagnetotellurics method
Row data acquisition, the present invention refuses this limitation.
Controlled-source audiomagnetotellurics method (Controlled Source Audio-frequency
Magnetotellurics, abbreviation CSAMT), also known as controlled-source audiomagnetotellurics sounding utilizes artificial field source excitation underground to be situated between
Matter, and electric field and the magnetic field that different frequency of supplies are produced are received away from lower in appropriate transmitting-receiving, because the field of different frequency is on stratum
In propagation depth it is different, reflect the depth also just mathematical relationship certain with frequency composition, the rock of different electrical conductivity is in electricity
Stream excites lower produced electric field and magnetic field to be different, and CSAMT methods are exactly by observing electric field and magnetic field intensity on measuring point
Change, and then the apparent resistivity and impedance phase of different frequency are calculated, surveyed with a kind of electromagnetism for obtaining underground medium Electrical distribution
Spy method.
Data collecting system based on CSAMT is constituted by launching (Tx) portion and receiving (Rx) portion two parts, and artificial field source has
Two kinds of grounded source and magnetic source, overcome the weak deficiency of natural field source signal, lateral resolution larger with depth of exploration scope
Height, can delicately find tomography, further, since acceptance division can receive magnetic field while electric field is received, therefore shield layer with high resistivity is made
With small, operating efficiency is substantially increased, labor intensity is alleviated.
Specifically, emission part such as emitter is using artificial dipole source to earth medium transmission signal, and the signal is generally electricity
Magnetic field signal, by changing the frequency of electromagnetic field, reaches the electrical structure of Underground medium.Wherein, electromagnetic field include sky wave,
Surface wave, formation wave etc., the frequency of electromagnetic field signal can be configured according to actual conditions, not be limited herein.It should be appreciated that
The emission part that electromagnetic field signal can be launched is all within the scope of the present invention.
Acceptance division receives the electric field signal and magnetic field signal of earth medium transmitting, using all electric fields (E) measured and
The amplitude and phase of magnetic field (H) component, calculate resistivity and phase difference etc..
With reference to Fig. 2, the electric field signal of each launching frequency, power frequency and its harmonic wave interference are believed using the time-domain method of sampling
Number sampled, sample rate is generally the integral multiple of the electromagnetic field signal of correspondence each launching frequency, utilizes time-domain sampling side
Method is sampled to the magnetic field signal of launching frequency, power frequency and its harmonic interference signals., can be to adopting when carrying out time-domain sampling
Sample frequency, sample rate, sampling length optimize design.Fourier transformation is carried out to the data after sampling, electric field signal is obtained
Real part imaginary part, the real part imaginary part of magnetic field signal.The electricity that the frequency of the electromagnetic field signal of emitter transmitting is received with receiver
The frequency of field signal or magnetic field signal is identical.It is described in detail below.
Alternatively, it is described that the electric field signal is carried out using the sample rate and using length when performing step S120
Sampling, before obtaining electric field information, in addition to：
The electric field signal is filtered；
The electric field signal after after filtering is amplified respectively；
The electric field signal after amplification is digitized respectively.
Alternatively, it is described that the magnetic field signal is carried out using the sample rate and using length when performing step S120
Sampling, before obtaining Magnetic Field, in addition to：
The magnetic field signal is filtered；
The magnetic field signal after after filtering is amplified respectively；
The magnetic field signal after amplification is digitized respectively.
Specifically, can be using high-pass filtering or LPF etc. when being filtered to electric field signal and magnetic field signal
Deng of the invention to refuse this limitation.
The explanation by taking electric field signal as an example.
The electric field signal received includes useful signal (EEffectively), random interfering signal (EAt random) and 50Hz fundamental frequencies and
Interference signal (E based on its harmonic waveInterference).Electric field expression formula is as follows：
Wherein, f1For the frequency of electric field signal；f0For the frequency of 50Hz fundamental frequencies and its harmonic wave, (frequency of harmonic wave is base
The integral multiple of wave frequency rate), t is time, A1It is generally constant for the amplitude of electric field signal；A0For shaking for 50Hz fundamental frequencies
Width, generally constant；ψ1For the phase of electric field signal, ψ0For the phase of 50Hz fundamental frequencies, generally constant.
Random interfering signal is mainly the random signal of Natural electromagnetic field and the noise of receiver background, random disturbances letter
It is number relatively low to the influence degree of gathered data, it can be neglected herein, main compacting 50Hz fundamental frequencies and its harmonic wave are to collection
The influence of data.
Influence for compacting fundamental frequency and its harmonic wave to gathered data, when carrying out discrete Fourier transform so that Fu
In leaf transformation real and imaginary parts and be zero, that is, it is measurement frequency f to ensure that n discrete data point will be met simultaneously1Complete cycle
Phase, f0+f1Complete cycle and f0-f1Complete cycle.The real part A of discrete Fourier transformG real partsWith imaginary part AG imaginary partsRespectively such as the He of formula 2
Shown in formula 3：
Wherein, f1For the frequency of electric field signal, f0For the frequency of 50Hz fundamental frequencies and its harmonic wave, (frequency of harmonic wave is base
The integral multiple of wave frequency rate), t is time, A0For the amplitude of 50Hz fundamental frequencies, generally constant, ψ0For the phase of 50Hz fundamental frequencies
Position, generally constant.
Alternatively, when performing step S120, the calculating obtains sample rate and sampling length, including：
The frequency of electric field signal or the frequency and the product of predetermined threshold value of magnetic field signal are calculated, sample rate is obtained, it is described pre-
If threshold value is the value more than or equal to the first numerical value；
The ratio and the product of second value of sample rate and the frequency of fundamental wave are calculated, sampling length, second number is obtained
Value is by calculating the ratio of the frequency of fundamental wave and the frequency of electric field signal or the frequency of magnetic field signal, being obtained with the product of third value
Arrive, the second value and the third value are positive integer.
Upper two formula is respectively the real and imaginary parts of the discrete Fourier transform to fundamental frequency, wants to obtain measurement frequency point f1
Amplitude, phase is by PERIODIC INTERFERENCE frequency f0Influence, it is summed and should be 0, that is, needs to ensure n discrete data point not
Only it is measurement frequency f1Complete cycle, while it is f that need to meet0+f1And f0-f1The complete cycle of frequency, such PERIODIC INTERFERENCE frequency
f0Influence just can obtain very big compacting.
For example, the sample rate of receiver is samp, sampling length is n, then the sampling time is n/samp, sample frequency
Samp/n need to meet below equation：
Wherein, n1、n2、n3For positive integer.
Formula (4) is set up, and premise isFor integer.
Sample rate is generally the integral multiple of electric field signal frequency, it is preferred that sample rate is more than or equal to electric field signal frequency
4 times of rate, i.e. the first numerical value are set to 4.Such as, when electric field signal frequency is 10000Hz, sample rate can be 40000Hz.
Formula 4 is simplified, obtaining sampling length n is：
Wherein, f0The generally integral multiple of 50Hz fundamental frequencies, herein f0Take 50Hz.
A is to causeObtained value is positive integer and n3For the value of positive integer, it is preferred that a value should be all desirable
Positive integer in volume minimum positive integer, for example, f0=50Hz, f1=4000Hz, then n3=8, a=1.
Alternatively, it is described that the electric field signal is carried out using the sample rate and using length when performing step S130
Sampling, obtains electric field information, and the magnetic field signal is sampled using the sample rate and using length, Magnetic Field is obtained,
Including：
Conversion process is carried out to the electric field signal after sampling using discrete Fourier transform；
Conversion process is carried out to the magnetic field signal after sampling using discrete Fourier transform.
Alternatively, when performing step S140, the electrical structure at least includes resistivity, described based on electric field letter
Breath and Magnetic Field, calculate the electrical structure of the earth medium, including：
Calculate electric field information and Magnetic Field ratio square, the cycle of the cycle of electric field signal or magnetic field signal and the
The product of four numerical value, obtains the resistivity of earth medium.
For example, to the data of collection carry out discrete Fourier transform after, respectively obtain electric field information such as electric-field intensity E with
Magnetic Field such as magnetic field intensity H, resistivity is calculated by below equation：
Wherein, E is electric-field intensity, and H is magnetic field intensity, and T is the cycle of electric field signal or the cycle of magnetic field signal.
Technique according to the invention scheme, is grown by the frequency of optimization design emitter and the sample rate of receiver and sampling
Degree, can effectively suppress the influence of fundamental frequency and its harmonic wave to gathered data, improve acquisition precision, it is ensured that gathered data it is accurate
Property.
Fig. 3 is that a kind of structure of controlled-source audiomagnetotellurics sounding data acquisition device provided in an embodiment of the present invention is shown
It is intended to.With reference to Fig. 3, the device includes：Receiving unit 31, first processing units 32, the processing of second processing unit 33 and the 3rd are single
Member 34.
Receiving unit 31, for receiver/transmitter launch propagated through earth medium after electric field signal and magnetic field signal；
First processing units 32, for frequency or the frequency of magnetic field signal based on electric field signal, and fundamental frequency
Frequency, obtains sample rate and sampling length；
Second processing unit 33, for being sampled using the sample rate and using length to the electric field signal, is obtained
Electric field information, is sampled to the magnetic field signal using the sample rate and using length, obtains Magnetic Field；
3rd processing unit 34, for based on the electric field information and Magnetic Field, calculating the electrical of the earth medium
Structure, wherein, the fundamental wave is the artificial electric power interference signal of influence electric field signal and magnetic field signal.
Alternatively, the first processing units 32 specifically for：
The frequency of electric field signal or the frequency and the product of predetermined threshold value of magnetic field signal are calculated, sample rate is obtained, it is described pre-
If threshold value is the value more than or equal to the first numerical value；
The ratio and the product of second value of sample rate and the frequency of fundamental wave are calculated, sampling length, second number is obtained
Value is by calculating the ratio of the frequency of fundamental wave and the frequency of electric field signal or the frequency of magnetic field signal, being obtained with the product of third value
Arrive, the second value and the third value are positive integer.
Alternatively, the second processing unit 33 specifically for：
Conversion process is carried out to the electric field signal after sampling using discrete Fourier transform；
Conversion process is carried out to the magnetic field signal after sampling using discrete Fourier transform.
Alternatively, the electrical structure at least include resistivity, the 3rd processing unit 34 specifically for：
Calculate electric field information and Magnetic Field ratio square, the cycle of the cycle of electric field signal or magnetic field signal and the
The product of four numerical value, obtains the resistivity of earth medium.
Alternatively, the second processing unit 33 is additionally operable to：
The electric field signal is filtered；
The electric field signal after after filtering is amplified respectively；
The electric field signal after amplification is digitized respectively.
The technique effect and preceding method embodiment phase of the device that the embodiment of the present invention is provided, its realization principle and generation
Together, to briefly describe, device embodiment part does not refer to part, refers to corresponding contents in preceding method embodiment.Affiliated neck
The technical staff in domain can be understood that, for convenience and simplicity of description, the tool of system described above, device and unit
Body running process, may be referred to the corresponding process in above method embodiment, will not be repeated here.
, can be by others side in embodiment provided by the present invention, it should be understood that disclosed apparatus and method
Formula is realized.Device embodiment described above is only schematical, for example, the division of the unit, only one kind are patrolled
Collect function to divide, there can be other dividing mode when actually realizing, in another example, multiple units or component can combine or can
To be integrated into another system, or some features can be ignored, or not perform.It is another, it is shown or discussed each other
Coupling or direct-coupling or communication connection can be the INDIRECT COUPLING or communication link of device or unit by some communication interfaces
Connect, can be electrical, machinery or other forms.
The unit illustrated as separating component can be or may not be it is physically separate, it is aobvious as unit
The part shown can be or may not be physical location, you can with positioned at a place, or can also be distributed to multiple
On NE.Some or all of unit therein can be selected to realize the mesh of this embodiment scheme according to the actual needs
's.
In addition, each functional unit in the embodiment that the present invention is provided can be integrated in a processing unit, also may be used
To be that unit is individually physically present, can also two or more units it is integrated in a unit.
If the function is realized using in the form of SFU software functional unit and is used as independent production marketing or in use, can be with
It is stored in a computer read/write memory medium.Understood based on such, technical scheme is substantially in other words
The part contributed to prior art or the part of the technical scheme can be embodied in the form of software product, the meter
Calculation machine software product is stored in a storage medium, including some instructions are to cause a computer equipment (can be individual
People's computer, server, or network equipment etc.) perform all or part of step of each of the invention embodiment methods described.
And foregoing storage medium includes：USB flash disk, mobile hard disk, read-only storage (ROM, Read-Only Memory), arbitrary access are deposited
Reservoir (RAM, Random Access Memory), magnetic disc or CD etc. are various can be with the medium of store program codes.
It should be noted that：Similar label and letter represents similar terms in following accompanying drawing, therefore, once a certain Xiang Yi
It is defined in individual accompanying drawing, then it further need not be defined and explained in subsequent accompanying drawing, in addition, term " the
One ", " second ", " the 3rd " etc. are only used for distinguishing description, and it is not intended that indicating or implying relative importance.
Finally it should be noted that：Embodiment described above, is only the embodiment of the present invention, to illustrate the present invention
Technical scheme, rather than its limitations, protection scope of the present invention is not limited thereto, although with reference to the foregoing embodiments to this hair
It is bright to be described in detail, it will be understood by those within the art that：Any one skilled in the art
The invention discloses technical scope in, it can still modify to the technical scheme described in previous embodiment or can be light
Change is readily conceivable that, or equivalent substitution is carried out to which part technical characteristic；And these modifications, change or replacement, do not make
The essence of appropriate technical solution departs from the spirit and scope of technical scheme of the embodiment of the present invention.The protection in the present invention should all be covered
Within the scope of.Therefore, protection scope of the present invention described should be defined by scope of the claims.
Claims (8)
1. a kind of controlled-source audiomagnetotellurics sounding collecting method, it is characterised in that including：
Receiver/transmitter transmitting propagated through earth medium after electric field signal and magnetic field signal；
The frequency of frequency or magnetic field signal based on electric field signal, and fundamental frequency frequency, obtain sample rate and sampling grow
Degree；
The electric field signal is sampled using the sample rate and using length, electric field information is obtained, using the sample rate with
The magnetic field signal is sampled using length, Magnetic Field is obtained；
Based on the electric field information and Magnetic Field, the electrical structure of the earth medium is calculated, wherein, the fundamental wave is influence
The artificial electric power interference signal of electric field signal and magnetic field signal.
2. the method as described in claim 1, it is characterised in that the calculating obtains sample rate and sampling length, including：
The frequency of electric field signal or the frequency and the product of predetermined threshold value of magnetic field signal are calculated, sample rate, the default threshold is obtained
It is worth for the value more than or equal to the first numerical value；
The ratio and the product of second value of sample rate and the frequency of fundamental wave are calculated, sampling length is obtained, the second value leads to
Cross and calculate the ratio of the frequency of fundamental wave and the frequency of electric field signal or the frequency of magnetic field signal, obtained with the product of third value,
The second value and the third value are positive integer.
3. the method as described in claim 1, it is characterised in that described to believe using the sample rate and using length the electric field
Number sampled, obtain electric field information, the magnetic field signal is sampled using the sample rate and using length, magnetic field is obtained
Information, including：
Conversion process is carried out to the electric field signal after sampling using discrete Fourier transform；
Conversion process is carried out to the magnetic field signal after sampling using discrete Fourier transform.
4. the method as described in claim 1, it is characterised in that the electrical structure at least includes resistivity, described to be based on institute
Electric field information and Magnetic Field are stated, the electrical structure of the earth medium is calculated, including：
Calculate electric field information and Magnetic Field ratio square, the cycle of electric field signal or the cycle of magnetic field signal and the 4th number
The product of value, obtains the resistivity of earth medium.
5. a kind of controlled-source audiomagnetotellurics sounding data acquisition device, it is characterised in that including：
Receiving unit, for receiver/transmitter launch propagated through earth medium after electric field signal and magnetic field signal；
First processing units, for the frequency of frequency or the frequency of magnetic field signal based on electric field signal, and fundamental frequency, are obtained
To sample rate and sampling length；
Second processing unit, for being sampled using the sample rate and using length to the electric field signal, obtains electric field letter
Breath, is sampled to the magnetic field signal using the sample rate and using length, obtains Magnetic Field；
3rd processing unit, for based on the electric field information and Magnetic Field, calculating the electrical structure of the earth medium, its
In, the fundamental wave is the artificial electric power interference signal of influence electric field signal and magnetic field signal.
6. device as claimed in claim 5, it is characterised in that the first processing units specifically for：
The frequency of electric field signal or the frequency and the product of predetermined threshold value of magnetic field signal are calculated, sample rate, the default threshold is obtained
It is worth for the value more than or equal to the first numerical value；
The ratio and the product of second value of sample rate and the frequency of fundamental wave are calculated, sampling length is obtained, the second value leads to
Cross and calculate the ratio of the frequency of fundamental wave and the frequency of electric field signal or the frequency of magnetic field signal, obtained with the product of third value,
The second value and the third value are positive integer.
7. device as claimed in claim 5, it is characterised in that the second processing unit specifically for：
Conversion process is carried out to the electric field signal after sampling using discrete Fourier transform；
Conversion process is carried out to the magnetic field signal after sampling using discrete Fourier transform.
8. device as claimed in claim 5, it is characterised in that the electrical structure is at least included at resistivity, the described 3rd
Manage unit specifically for：
Calculate electric field information and Magnetic Field ratio square, the cycle of electric field signal or the cycle of magnetic field signal and the 4th number
The product of value, obtains the resistivity of earth medium.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710574423.8A CN107329183B (en) | 2017-07-14 | 2017-07-14 | A kind of controlled-source audiomagnetotellurics sounding collecting method and device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710574423.8A CN107329183B (en) | 2017-07-14 | 2017-07-14 | A kind of controlled-source audiomagnetotellurics sounding collecting method and device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107329183A true CN107329183A (en) | 2017-11-07 |
| CN107329183B CN107329183B (en) | 2019-03-22 |
Family
ID=60226982
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710574423.8A Active CN107329183B (en) | 2017-07-14 | 2017-07-14 | A kind of controlled-source audiomagnetotellurics sounding collecting method and device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107329183B (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109470135A (en) * | 2018-11-12 | 2019-03-15 | 吉林大学 | CSAMT data inactivity bearing calibration |
| CN114325851A (en) * | 2021-12-31 | 2022-04-12 | 华中科技大学 | Method and system for detecting underwater ferromagnetic target based on multi-frequency radiation source |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5229765A (en) * | 1991-05-08 | 1993-07-20 | Halliburton Logging Services, Inc. | SP noise cancellation technique |
| CN101251606A (en) * | 2008-04-09 | 2008-08-27 | 吉林大学 | Interference suppression circuit for industrial frequency harmonic of desired signal in frequency band of weak signal detection apparatus |
| CN102426393A (en) * | 2011-11-16 | 2012-04-25 | 中国地质大学（北京） | Electric prospecting method and device |
| CN103389514A (en) * | 2013-07-18 | 2013-11-13 | 中国科学院地质与地球物理研究所 | MT (magnetotelluric) denoising device and method |
| CN103795473A (en) * | 2012-11-02 | 2014-05-14 | 华为技术有限公司 | Method and system for eliminating power-frequency interference |
-
2017
- 2017-07-14 CN CN201710574423.8A patent/CN107329183B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5229765A (en) * | 1991-05-08 | 1993-07-20 | Halliburton Logging Services, Inc. | SP noise cancellation technique |
| CN101251606A (en) * | 2008-04-09 | 2008-08-27 | 吉林大学 | Interference suppression circuit for industrial frequency harmonic of desired signal in frequency band of weak signal detection apparatus |
| CN102426393A (en) * | 2011-11-16 | 2012-04-25 | 中国地质大学（北京） | Electric prospecting method and device |
| CN103795473A (en) * | 2012-11-02 | 2014-05-14 | 华为技术有限公司 | Method and system for eliminating power-frequency interference |
| CN103389514A (en) * | 2013-07-18 | 2013-11-13 | 中国科学院地质与地球物理研究所 | MT (magnetotelluric) denoising device and method |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109470135A (en) * | 2018-11-12 | 2019-03-15 | 吉林大学 | CSAMT data inactivity bearing calibration |
| CN114325851A (en) * | 2021-12-31 | 2022-04-12 | 华中科技大学 | Method and system for detecting underwater ferromagnetic target based on multi-frequency radiation source |
| CN114325851B (en) * | 2021-12-31 | 2024-05-10 | 华中科技大学 | Method and system for detecting underwater ferromagnetic target based on multi-frequency radiation source |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107329183B (en) | 2019-03-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105911603B (en) | Four-dimensional geophysical prospecting method based on natural electric field | |
| RU2349936C2 (en) | Electromagnetic study of hydrocarbon collectors | |
| CN103869371B (en) | Manual field source frequency domain full-gradient electromagnetic measuring method | |
| CN104656157B (en) | A kind of method and device in identification shale gas dessert area | |
| CN105116453B (en) | The transient electromagnetic exploitation method and device of a kind of gas hydrate in permafrost | |
| CN103389514B (en) | A kind of mt denoising device and mt denoising method | |
| CN105607124B (en) | Seismic wave near surface interval quality factors compensation method and device | |
| WO2009037419A2 (en) | Electromagnetic surveying | |
| CN109782360A (en) | A kind of low-resistance area of coverage deep sand body detection method | |
| CN109085653A (en) | A kind of detection method of geology of deep part, sulfide ore body resource | |
| CN104714254B (en) | A kind of method and apparatus that colliery ponding goaf is detected based on transient electromagnetic method | |
| MX2011009449A (en) | Electromagnetic survey using naturally occurring electromagnetic fields as a source. | |
| CN102707323A (en) | Controllable source audio-frequency magnetic field sounding method for geological exploration | |
| Wu et al. | Ground-source Airborne Time-domain ElectroMagnetic (GATEM) modelling and interpretation method for a rough medium based on fractional diffusion | |
| CN107329183A (en) | A kind of controlled-source audiomagnetotellurics sounding collecting method and device | |
| Gloaguen et al. | Stochastic borehole radar velocity and attenuation tomographies using cokriging and cosimulation | |
| CN108169802A (en) | A kind of time domain electromagnetic data slow diffusion imaging method of harsh media model | |
| CN108897050A (en) | A kind of controlled-source audiomagnetotellurics sounding method is most preferably received and dispatched away from determining method and system | |
| CN108594301B (en) | A kind of method and processing terminal of the seismic data fusion with difference characteristic | |
| CN112327377B (en) | Method for calculating resistivity of transition region and near region by controllable source electromagnetic method | |
| WO2021098667A1 (en) | Full-coverage anti-interference controlled-source electromagnetic detection method | |
| Xu et al. | Groundwater resources survey of Tongchuan city using the audio magnetotelluric method | |
| Campaña et al. | Inversion of TEM data and analysis of the 2D induced magnetic field applied to the aquifers characterization in the Paraná basin, Brazil | |
| CN106610509B (en) | A kind of time-domain processing method of transient electromagnetic field data | |
| CN115128679A (en) | Frequency domain electromagnetic sounding method and system and electronic equipment |
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 |