CN107024724A - The method that deconvolution rapid Optimum helicopter aviation electromagnetic detects early time data - Google Patents
The method that deconvolution rapid Optimum helicopter aviation electromagnetic detects early time data Download PDFInfo
- Publication number
- CN107024724A CN107024724A CN201710479110.4A CN201710479110A CN107024724A CN 107024724 A CN107024724 A CN 107024724A CN 201710479110 A CN201710479110 A CN 201710479110A CN 107024724 A CN107024724 A CN 107024724A
- Authority
- CN
- China
- Prior art keywords
- deconvolution
- aviation electromagnetic
- data
- response
- electromagnetic
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/38—Processing data, e.g. for analysis, for interpretation, for correction
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/15—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for use during transport, e.g. by a person, vehicle or boat
- G01V3/165—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for use during transport, e.g. by a person, vehicle or boat operating with magnetic or electric fields produced or modified by the object or by the detecting device
Abstract
The present invention relates to a kind of method that deconvolution rapid Optimum helicopter aviation electromagnetic detects early time data, the present invention is that the computational methods of deconvolution are incorporated into the time domain aviation electromagnetic Data processing that transmitted waveform is trapezoidal wave, and reasonably improved, remove the characteristic response of system in itself, SMS message data track effective time is chosen to shift to an earlier date, and then it is helpful to detection near the ground, aviation electromagnetic investigative range is added, to solving the problems, such as that environment, the energy and military aspect are of great importance.Remove system by deconvolution to respond, it is simple and quick, efficiency high, solve due to the problem of the presence of system then responds concussion in early days, optimizing the quality of data in SMS message road, switch off current section early stage data available is shifted to an earlier date into 4 to 5 points, pre-set time is about 0.13ms, the data application in early stage road is become possible to, through skin depth formula scales, the effective scope of detection of aviation electromagnetic detection is added.
Description
Technical field
The present invention relates to a kind of deconvolution rapid Optimum helicopter aviation electromagnetic detection system early time data method, when belonging to
Between domain aviation electromagnetic data processing method.
Background technology
Aeroelectromagnetic method (Airborne Electromagnetic, AEM) is for carrying with aircraft (being usually aircraft)
Platform, a kind of Geophysical Exploration Method based on electromagnetic coupling theory, with speed is fast, cost is low, road ability good, can
The advantages such as large area covering, are widely used in basic geological survey, mineral resources, oil-gas exploration, and the hydrology, engineering, ring
The field such as reconnoitre in border.
The primary frequency range of Time-domain AEM is 5 to 25kHz, during data acquisition inevitably
Can be by various interference, and these noises are usually related to the factor such as acquisition time, collecting location and acquisition system.Aviation
The basic goal of electromagnetic surveying be by measure come data analyze, and then seek the distribution feelings of underground conducting medium
Condition, thus estimates the resource distribution of underground.Therefore, aviation electromagnetic data processing method is the main research work of aviation electromagnetic detection
One of make.But when instrument or set of system are being observed and recording a physical change, resulting measurement result shows thing
Phenomenon and process are managed, the characteristic of instrument or system is also included.The non-ideal characteristic of instrument system, such as system bandwidth, inside are made an uproar
Sound, non ideal operation etc., all make obtained observation and record degrade, this to degrade so that it is not true to test obtained data
As a result.
In the measurement for Time-domain AEM data, more severe, the involved test system of this influence
Characteristic, such as coupling of transmitting coil, receiving coil and aircraft skin all can be influence data factor, these factors
It is referred to as system performance, these system performances have impact on the early time data or frequency domain medium-high frequency letter in time-domain after switch off current
Number data quality.Early time data in time-domain after switch off current is contained in the data of frequency domain high frequency signal
The dampening information of the secondary field of substantial amounts of aviation electromagnetic detection, the underground effective information that this part attenuation curve is included is most,
It is especially most to subaerial message reflection.Only because in actual conditions, the transmitting magnetic moment of emission current causes greatly to close very much
Disconnected delay is longer, also causes the factors such as the system performances such as circuit concussion to result in early signal using ignored.
So the data of measurement are because there is the influence of much noise, obtained test result is not the sense of preferable measurement point
Electromotive force is answered, it is due to greatly error responses caused by the characteristic of the systems such as system bandwidth in itself to also have.This mistake
Difference can make the early time data time invalid, make effective information omit so that aviation electromagnetic occurs in that blind on the problem of shallow-layer is detected
Area.It is all increasingly ripe in the research of hardware aspect in the world, such as how to reduce the drift of transmitting coil, reduce capacity coupled
Influence etc..
In the world, the data area that aviation electromagnetic detection in recent years can be corrected accurately gradually increases.Macnae and Baron-
Hay is in 100 μ s after the power-off by VTEM system early time data raisings in 2010, and in the same year, they solve transmitting electricity using deconvolution
Capacitance Coupled and bandwidth problem caused by stream, represent that model can improve 10-20 μ s after a loss of power.
It is domestic to study relatively fewer specifically designed for Time-domain AEM early time data, should not by Deconvolution Technique yet
Use in aviation electromagnetic detection.
The content of the invention
Present invention aims at the deficiency for above-mentioned limit purchase technology, there is provided a kind of simple, quick helicopter aviation electricity
Magnetic detection early time data optimization method.
This method pin is to be responded influence by system to the measurement data of helicopter airborne electromagnetic system to cause early time data
Disabled problem.First, response and transmitting during the high-altitude flight detected using aviation electromagnetic, system is obtained using deconvolution
Response, secondly using the response on survey line, deconvolution calculating obtains truly responding greatly to remove system response again, last profit
With design current, high-altitude and survey line Current calculation correction factor standardize it to correct the earth response.
The present invention is that the computational methods of deconvolution are incorporated into the time domain aviation electromagnetic data that transmitted waveform is trapezoidal wave
In processing, and reasonably improved, remove the characteristic response of system in itself, chosen SMS message data track effective time and carry
Before, and then it is helpful to detection near the ground, aviation electromagnetic investigative range is added, to solving environment, the energy and military affairs side
The problem of face, is of great importance.
The method that deconvolution rapid Optimum helicopter aviation electromagnetic detects early time data, comprises the following steps：
A, typing aviation electromagnetic detect the induced electromotive force in high-altitude and are superposed to a measuring point；
B, typing aviation electromagnetic detect the emission current in high-altitude and are superposed to a measuring point；
C, deconvolution calculate the system response of airborne electromagnetic system；
D, each survey line of typing aviation electromagnetic induced electromotive force；
E, deconvolution is utilized to calculate the response of removal system；
F, typing aviation electromagnetic flight survey line emission current；
G, calculating correction factor；
The response that h, correction are gone after system response, obtains really responding greatly；
I, export result, are depicted as profile.
Beneficial effect：The present invention is to remove system by deconvolution to respond, and simple and quick, efficiency high makes what is finally obtained
Response is more nearly reality greatly, and Fig. 2 gives the result for taking off the earth response after system response；Correction factor can quantify
The height of this data quality；The present invention reduction is true respond greatly on the basis of, solve due to system presence then
The problem of early stage response concussion, the quality of data in SMS message road is optimized, the earth that Fig. 3 gives after the response of removal system rings
Should.As a result switch off current section early stage data available is shifted to an earlier date into 4 to 5 points, pre-set time is about 0.13ms, makes the data in early stage road
Using becoming possible to, Fig. 4 is given for same SMS message, Profile Correlation figure before optimization and after optimization；Early time data quality
Raising to it is near the ground detection have an unexpected effect, through skin depth formula scales, in typical the earth (100 Ω/m),
The effective of aviation electromagnetic detection can be added by the investigation depth of aviation electromagnetic detection system to 50 meters of shallow-layer extension in theory
Investigative range.
Brief description of the drawings
Fig. 1 is the method flow diagram that deconvolution rapid Optimum helicopter aviation electromagnetic detects early time data
Fig. 2 is that the earth is responded, survey line emission current and side line response curve.
Fig. 3 attenuation curves switch off current section early stage induced electromotive force amplification comparison diagram
Fig. 4 is that attenuation curve takes out into profile
Fig. 4 a surveys line remove profile when starting point is 184 after system response
Fig. 4 b surveys line do not go system to respond profile when Hou Chou roads starting point is 184
Embodiment
The present invention is described in further detail with reference to the accompanying drawings and examples：
The method that deconvolution rapid Optimum helicopter aviation electromagnetic detects early time data, comprises the following steps：
A, typing aviation electromagnetic detect the induced electromotive force in high-altitude and are superposed to a measuring point；
B, typing aviation electromagnetic detect the emission current in high-altitude and are superposed to a measuring point；
C, deconvolution calculate the system response of airborne electromagnetic system；
D, each survey line of typing aviation electromagnetic induced electromotive force；
E, deconvolution is utilized to calculate the response of removal system；
F, typing aviation electromagnetic flight survey line emission current；
G, calculating correction factor；
The response that h, correction are gone after system response, obtains really responding greatly；
I, export result, are depicted as profile.
First, the typing aviation electromagnetic detection system induced electromotive force data that measurement is received in high-altitude test flight, if
For V0(t).Emission current of the typing aviation electromagnetic detection system in high-altitude test flight, is set to I again0(t)。
In high hollow testing, aviation electromagnetic detection system equally opens emitter emission current, and receiver receives big underground
Conductor is produced because of electromagnetic induction.High space-time regards remote enough apart from ground as, does not receive secondary field, so now in height
It is emission current through being responded obtained from system performance function that the response received in the air, which is considered as excitation, if aviation electromagnetic is detected
The unit impact response of the system of system is H0(t), then its frequency domain representation is exactly H0(ω)。
To overcome the problem of time domain input signal is with system convolution length, we are instant with convolution theorem here
Domain convolution is multiplied equal to frequency domain, therefore the present invention in frequency domain be related to the related operation of convolution.By high aerial transmitting
Electric current I0(t) Fast Fourier Transform (FFT) is done for I0(ω), the induced electromotive force V that high-altitude is received0(t) fast Flourier change is done
Get V in return0(ω).It can be seen from Signals ＆ Systems knowledge and linear convolution theorem, the transmission function of aviation electromagnetic detection system
It can be expressed as
H0(ω)=V0(ω)/I0(ω) (1)
Secondly, the induced electromotive force data that typing aviation electromagnetic detection system is measured in survey line flight are set to V (t), done
FFT obtains V (ω).The response that system during due to survey line flight is received is that a part comes from underground
Conductor senses the secondary field that curent change is produced, and another part is metal shell from detection system etc., so root again
Understood according to convolution theorem, be after the response of removal system
It is preliminary to have obtained the secondary field that only subsurface anomaly body is produced.
Finally, correction factor is calculated, correction the earth is responded.
Emission current data I (t) of the typing aviation electromagnetic detection system in survey line flight, doing FFT is
I(ω).The standard emission electric current being in-flight pre-designed is C (t) by we, (fast fourier transformed into C (ω)) to make
For reference.No matter in high-altitude or survey line, each actual electric current is that small drift and jump is undergone on reference current
Jump, so, the Fourier transform of actual current is not exclusively equal to C (ω).In order to correct this error, one is multiplied by frequency domain
Individual correction factor, this correction factor is defined as ratio (I (ω)/C (ω)) * (C (ω)/A (ω)) by we.This correction because
Son mainly has two effects, and one is to obtain the matter that size may determine that emission current on high-altitude and side line by this rate value
Amount, can also rough estimate in advance response quality, second, the response of the earth can be corrected to a certain degree, make the earth respond more
Close to actual value.So, the principle formula of deconvolution is just write as
R (ω)=H (ω) * (I (ω)/C (ω)) * (C (ω)/A (ω)) (3)
It can be seen that when ratio C (ω)/A (ω) corrects drift and saltus step, ratio I (ω)/C (ω) provides ideal
The deconvolution of waveform.Possible C (ω) seems to remove on the position of molecule and denominator, but needs in last warp
Once such processing is carried out before product.
Embodiment 1
By taking the good field measured data of a flight progress as an example：
First, the induced electromotive force number that entry time domain helicopter type aviation electromagnetic detection system is measured in high-altitude flight
According to, if one group of m measuring point, the induction electric received during the aviation electromagnetic detection system high-altitude flight of each n sampled point of measuring point
Gesture data V0(n × m), after be superimposed as a measuring point V0(n)。
The emission current of typing aviation electromagnetic detection system correspondence induced electromotive force in high-altitude flight, if one group of m survey
Point, the current data I launched during the aviation electromagnetic detection system high-altitude flight of each n sampled point of measuring point0(n×m)。
The induced electromotive force V that high-altitude is received0(n) do Fourier transform and obtain V0(ω), high aerial emission current I0
(n) do Fourier transform and obtain I0(ω).If the unit impact response of the system of aviation electromagnetic detection system is H0(t), then it
Frequency domain representation is exactly H0(ω), then system response now can be expressed as
H0(ω)=V0(ω)/I0(ω) (6)
Secondly, the induced electromotive force data that typing aviation electromagnetic detection system is measured in survey line flight, if one group of m survey
Point, the induced electromotive force data V (n × m) received during the aviation electromagnetic detection system survey line flight of each n sampled point of measuring point,
FFT is done into Vn×m(ω), the response that system during due to survey line flight is received is that a part comes from underground
Conductor sense curent change produce secondary field, another part is metal shell from detection system etc., so again
It can be seen from convolution theorem, it is after the response of removal system:
Finally, correction factor is calculated, correction the earth is responded.
Emission current data I (n) of the typing aviation electromagnetic detection system in survey line flight, it is fast fourier transformed into
I(ω).The standard emission electric current being in-flight pre-designed is C (n) by we, (fast fourier transformed into C (ω)) to make
For reference.No matter in high-altitude or survey line, each actual electric current is that small drift and jump is undergone on reference current
Jump, so, the Fourier transform of actual current is not exclusively equal to C (ω).In order to correct this error, one is multiplied by frequency domain
Individual correction factor, this correction factor is defined as ratio (I (ω)/C (ω)) * (C (ω)/A (ω)) by we.This correction because
Son mainly has two effects, and one is to obtain the matter that size may determine that emission current on high-altitude and side line by this rate value
Amount, can also rough estimate in advance response quality, second, the response of the earth can be corrected to a certain degree, make the earth respond more
Close to actual value.So, the principle formula of deconvolution is just write as
R (ω)=H (ω) * (I (ω)/C (ω)) * (C (ω)/A (ω)) (8)
It can be seen that when ratio C (ω)/A (ω) corrects drift and saltus step, ratio I (ω)/C (ω) provides ideal
The deconvolution of waveform.Possible C (ω) seems to remove on the position of molecule and denominator, but needs in last warp
Once such processing is carried out before product.
Claims (1)
1. a kind of deconvolution rapid Optimum helicopter aviation electromagnetic detection system early time data method, comprises the following steps：
A, typing aviation electromagnetic detect the induced electromotive force in high-altitude and are superposed to a measuring point；
B, typing aviation electromagnetic detect the emission current in high-altitude and are superposed to a measuring point；
C, deconvolution calculate the system response of airborne electromagnetic system；
D, each survey line of typing aviation electromagnetic induced electromotive force；
E, deconvolution is utilized to calculate the response of removal system；
F, typing aviation electromagnetic flight survey line emission current；
G, calculating correction factor；
The response that h, correction are gone after system response, obtains really responding greatly；
I, export result, are depicted as profile.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710479110.4A CN107024724A (en) | 2017-06-22 | 2017-06-22 | The method that deconvolution rapid Optimum helicopter aviation electromagnetic detects early time data |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710479110.4A CN107024724A (en) | 2017-06-22 | 2017-06-22 | The method that deconvolution rapid Optimum helicopter aviation electromagnetic detects early time data |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN107024724A true CN107024724A (en) | 2017-08-08 |
Family
ID=59531729
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710479110.4A Pending CN107024724A (en) | 2017-06-22 | 2017-06-22 | The method that deconvolution rapid Optimum helicopter aviation electromagnetic detects early time data |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107024724A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110794464A (en) * | 2019-11-07 | 2020-02-14 | 中国地质科学院地球物理地球化学勘查研究所 | Method for removing time domain aeronautical electromagnetic data background field |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20090309599A1 (en) * | 2008-06-11 | 2009-12-17 | Antoni Marjan Ziolkowski | Method for subsurface electromagnetic surveying using two or more simultaneously actuated electromagnetic sources |
| CN102621584A (en) * | 2011-12-14 | 2012-08-01 | 吉林大学 | Aeromagnetic vector data processing method |
| CN103675927A (en) * | 2013-12-20 | 2014-03-26 | 吉林大学 | Correction method for pendulum angle of receiving pod of airborne electromagnetic system in fixed wing aircraft |
| CN103760614A (en) * | 2014-02-24 | 2014-04-30 | 中国科学院电子学研究所 | Transient electromagnetic forward modeling method applicable to irregular transmitted waveforms |
| CN103995292A (en) * | 2014-06-09 | 2014-08-20 | 桂林电子科技大学 | Transient electromagnetic early signal reconstruction method |
| CN106226828A (en) * | 2016-07-18 | 2016-12-14 | 吉林大学 | Time domain aviation electromagnetic data constraint fitting of a polynomial leveling method |
-
2017
- 2017-06-22 CN CN201710479110.4A patent/CN107024724A/en active Pending
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20090309599A1 (en) * | 2008-06-11 | 2009-12-17 | Antoni Marjan Ziolkowski | Method for subsurface electromagnetic surveying using two or more simultaneously actuated electromagnetic sources |
| CN102621584A (en) * | 2011-12-14 | 2012-08-01 | 吉林大学 | Aeromagnetic vector data processing method |
| CN103675927A (en) * | 2013-12-20 | 2014-03-26 | 吉林大学 | Correction method for pendulum angle of receiving pod of airborne electromagnetic system in fixed wing aircraft |
| CN103760614A (en) * | 2014-02-24 | 2014-04-30 | 中国科学院电子学研究所 | Transient electromagnetic forward modeling method applicable to irregular transmitted waveforms |
| CN103995292A (en) * | 2014-06-09 | 2014-08-20 | 桂林电子科技大学 | Transient electromagnetic early signal reconstruction method |
| CN106226828A (en) * | 2016-07-18 | 2016-12-14 | 吉林大学 | Time domain aviation electromagnetic data constraint fitting of a polynomial leveling method |
Non-Patent Citations (1)
| Title |
|---|
| J.M.LEGAULT, ETC.: "Recent AEM Case Study Examples of a Full Waveform Time-Domain System for Near-Surface and Groundwater Applications", 《SECOND INTERNATIONAL CONFERENCE ON ENGINEERING GEOPHYSICS》 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110794464A (en) * | 2019-11-07 | 2020-02-14 | 中国地质科学院地球物理地球化学勘查研究所 | Method for removing time domain aeronautical electromagnetic data background field |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104597506B (en) | The air-ground electromagnetic exploration method of frequency domain | |
| CN102096113B (en) | Time-domain ground-air electromagnetic detection system and calibration method | |
| US20200348438A1 (en) | Method for Suppressing Airborne Transient Electromagnetic In-Band Vibration Noise | |
| US20110095763A1 (en) | Geologic mapping using lightning | |
| CN102832908A (en) | Wavelet transform and variable-step-size LMS (least mean square) adaptive filtering based signal denoising method | |
| CN103941254A (en) | Soil physical property classification recognition method and device based on geological radar | |
| CN104076404B (en) | Use multichannel to be concerned with and suppress the magnetic anomaly detection method of earth magnetism background noise | |
| CN108072910B (en) | A kind of distribution magnetic anomaly detection system environment magnetic compensation method | |
| CN103823244A (en) | Magnetic resonance three-component noise removing device and noise removing method | |
| CN104267440A (en) | Common middle point (CMP) detection method used for ground penetrating radar (GPR) | |
| CN103852785A (en) | Evaluation method for stratum anisotropy | |
| CN103743969A (en) | Measurement noise suppression method for near-field electromagnetic radiation interference measurement signals | |
| CN116484177A (en) | Motion-induced noise prediction elimination method for electromagnetic detection of flight platform | |
| CN103675927A (en) | Correction method for pendulum angle of receiving pod of airborne electromagnetic system in fixed wing aircraft | |
| CN109581481B (en) | Portable high-frequency controllable seismic source seismic signal harmonic interference elimination method | |
| CN101825722B (en) | Robust method for estimating instantaneous frequency of seismic signal | |
| CN107024724A (en) | The method that deconvolution rapid Optimum helicopter aviation electromagnetic detects early time data | |
| CN110794464B (en) | Method for removing time domain aeronautical electromagnetic data background field | |
| Fee et al. | Characterizing infrasound station frequency response using large earthquakes and colocated seismometers | |
| Legault | Ten years of passive airborne AFMAG EM development for mineral exploration | |
| CN103376443B (en) | Ground penetrating radar terrestrial interference detecting and fast eliminating method | |
| Keating | Improved use of the local wavenumber in potential-field interpretation | |
| CN105866552B (en) | The measurement method of the impedance of Aircraft Cable Shield layer | |
| EP2634604A2 (en) | Determination of particle displacement or velocity from particle acceleration measurements | |
| CN104863578A (en) | Electrical conductive rate well measuring method of radar in well |
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 | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20170808 |
|
| WD01 | Invention patent application deemed withdrawn after publication |