CN105093341B - The computational methods and system of the full period apparent resistivity of time-domain transition level electric field - Google Patents
The computational methods and system of the full period apparent resistivity of time-domain transition level electric field Download PDFInfo
- Publication number
- CN105093341B CN105093341B CN201510433717.XA CN201510433717A CN105093341B CN 105093341 B CN105093341 B CN 105093341B CN 201510433717 A CN201510433717 A CN 201510433717A CN 105093341 B CN105093341 B CN 105093341B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- electric field
- record
- time
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention provides the computational methods and system of a kind of full period apparent resistivity of time-domain transition level electric field, wherein method includes：（1）, obtain field acquisition obtained by time-domain record material：Including the electric field in two mutually perpendicular directionsWith,（2）, to step（1）In time-domain record material carry out temporal filtering, superposition, reject the processing such as flying spot；（3）, calculate the complete waveform of upper and lower step；（4）, calculateThe electric field of the upper and lower step in direction；（5）, pass through electric Field Calculation " full period " apparent resistivity of upper step。
Description
Technical field
The invention belongs to the electromagnetic survey neck that the geological prospectings such as oil, natural gas, untraditional reservoir, underground mine are explored
Domain, and in particular to the computational methods of the full period apparent resistivity of DEEP TRANSIENT ELECTROMAGNETIC SOUNDING.
Background technology
DEEP TRANSIENT ELECTROMAGNETIC SOUNDING is to launch pumping signal using artificial source, the electromagnetic field in research area observation time domain, with
A kind of electromagnetic exploration method of Study of The Underground electrical structure.Fig. 1 is showing for one of the field laying mode of DEEP TRANSIENT ELECTROMAGNETIC SOUNDING
It is intended to.Electric field in the physical quantity that observation station measures for two horizontal directions, i.e.,：exAnd e (t)y(t) (unit be volt/
Rice), wherein t is the time delay of observation (unit is the second).
In Fig. 1, AB is transmitting electric dipole, and by emitter to underground feed-in square wave current, current strength is I (peaces
Training), the dipole moment for launching electric dipole is PE=IAB (ampere per meter)),To receive and dispatch away from (AB midpoints O to receiving point P
Distance, unit rice), φ is the angle between transmitting electric dipole central point O and measuring point P line and transmitting electric dipole AB.
exAnd e (t)y(t) it is respectively and the horizontal electric field component of vertical direction parallel with transmitting electric dipole.
Usually, dutycycle is used as 1:1 square-wave waveform, as shown in Figure 2：T is the time delay (second) of a cycle,
In the positive power supply of the period of preceding a quarter, tight secondly power-off a quarter cycle, and tight secondly reverse power supply a quarter
Cycle (electric current is-I), last a quarter cycle power-off.
The conventional method of DEEP TRANSIENT ELECTROMAGNETIC SOUNDING data processing is to calculate early stage apparent resistance by the electromagnetic field component observed
Rate and late period apparent resistivity.But early stage apparent resistivity mainly reflects the resistivity on shallower stratum, late period apparent resistivity then only reflects
The situation of deep resistivity, both reflect unclear to transitional resistivity.
For example, according to the electric field observed in the power supply phase(in this symbol, x represents the direction of electric field, before+expression
The positive power-on time section of a quarter, " upper step " represent to carry out data collection from starting power supply) early stage for calculating regard it is electric
Resistance rateWith late period apparent resistivityExpression formula is respectively：
In order to obtain to " full period " apparent resistivity in late period from early stage to transitional period, intactly to reflect stratum
From shallow-layer to mid-deep strata again to deep layer change in resistance information, it is necessary to provide definition and the meter of " full period " apparent resistivity
Calculation method.Necessary to this " full period " apparent resistivity is also inverting.
The content of the invention
The present invention is for the insufficient of the early and late apparent resistivity of DEEP TRANSIENT ELECTROMAGNETIC SOUNDING and follow-up further inverting
Needs, provide a kind of method of the definition and calculating of " full period " apparent resistivity.
The technical solution adopted for the present invention to solve the technical problems is：
The invention provides a kind of computational methods of the full period apparent resistivity of time-domain transition level electric field, including it is following
Step：
The time-domain record material of nagneto-telluric field obtained by S1, acquisition field acquisition, obtains the electricity of two horizontal directions
Field component exAnd e (t)y(t)；
S2, to electric field component exAnd e (t)y(t) it is overlapped and filtering process, to eliminate interference signal；
S3, to electric field component exAnd e (t)y(t) upper and lower step electric field value carries out editing and processing, rejects what is deviated considerably from
Flying spot；
S4, calculatingThe horizontal component of electric field in direction
For in transmitting electric dipole
Angle between the line and transmitting electric dipole line of heart point and measuring point；
S5, further according to formula
Calculate full-time
Section apparent resistivity ρw(t), ρw(t) imply in the formula；
Wherein μ0=4 π × 10-7Henry/rice, t are time delay,It is to receive and dispatch away from erf (x) is that normal probability accumulates
Point or error function：
Wherein xpFor transition time parameter：
Wherein f (xp) it is normalized function, f (xp) it is on xpMonotonic function：
Present invention also offers a kind of computing system of the full period apparent resistivity of time-domain transition level electric field, including：
Data acquisition module, for obtaining the time-domain record material of the nagneto-telluric field obtained by field acquisition, obtain two
The electric field component e of individual horizontal directionxAnd e (t)y(t)；
Signal processing module, for electric field component exAnd e (t)y(t) it is overlapped and filtering process, to eliminate interference letter
Number；And for electric field component exAnd e (t)y(t) upper and lower step electric field value carries out editing and processing, rejects flying of deviating considerably from
Point；
Computing module, for calculatingThe horizontal component of electric field in direction
For in transmitting electric dipole
Angle between the line and transmitting electric dipole line of heart point and measuring point；
Further according to formula
Calculate full-time
Section apparent resistivity ρw(t), ρw(t) imply in the formula；
Wherein μ0=4 π × 10-7Henry/rice, t are time delay,It is to receive and dispatch away from erf (x) is that normal probability accumulates
Point or error function：
Wherein xpFor transition time parameter：
Wherein f (xp) it is normalized function, f (xp) it is on xpMonotonic function：
The beneficial effect comprise that：The present invention is when carrying out DEEP TRANSIENT ELECTROMAGNETIC SOUNDING exploration, by measuring transition
Electric field, the complete apparent resistivity of full period can be calculated, and then rational correctly resistivity-depth is obtained by inverting
Spend profile.
Brief description of the drawings
Below in conjunction with drawings and Examples, the invention will be further described, in accompanying drawing：
Fig. 1 is the schematic diagram of one of the field laying mode of DEEP TRANSIENT ELECTROMAGNETIC SOUNDING；
Fig. 2 be DEEP TRANSIENT ELECTROMAGNETIC SOUNDING to the dutycycle of underground feed-in be 1：1 square wave current；
Fig. 3 is that waveform is completely received in DEEP TRANSIENT ELECTROMAGNETIC SOUNDING a cycle；
Fig. 4 is the electric field wavy curve that theoretical model calculates；
Fig. 5 is the early stage that conventional method calculatesAnd late periodApparent resistivity curve
(step field is worth in use)；
Fig. 6 is the full period apparent resistivity curve (ρ calculated using the present inventionw(t), it is worth using upper step field)；
Fig. 7 is the comparison of apparent resistivity curve.
Embodiment
In order to make the purpose , technical scheme and advantage of the present invention be clearer, it is right below in conjunction with drawings and Examples
The present invention is further elaborated.It should be appreciated that specific embodiment described herein is only to explain the present invention, not
For limiting the present invention.
In embodiments of the invention, DEEP TRANSIENT ELECTROMAGNETIC SOUNDING is used, as shown in figure 1, the field for DEEP TRANSIENT ELECTROMAGNETIC SOUNDING
The schematic diagram of one of laying mode.Electric field in the physical quantity that observation station measures for two horizontal directions, i.e.,：exAnd e (t)y(t)
(unit is volts/meter), wherein t are the time delay of observation (unit is the second).
In Fig. 1, AB is transmitting electric dipole, and by emitter to underground feed-in square wave current, current strength is I (peaces
Training), the dipole moment for launching electric dipole is PE=IAB (ampere per meter)),To receive and dispatch away from (AB midpoints O to receiving point P
Distance, unit rice),To launch the angle between electric dipole central point O and measuring point P line and transmitting electric dipole AB.
exAnd e (t)y(t) it is respectively and the horizontal electric field component of vertical direction parallel with transmitting electric dipole.Dutycycle is used as 1:1
Square-wave waveform, as shown in Figure 2：T is the time delay (second) of a cycle, is powered in the period forward direction of preceding a quarter, tightly
Secondly power-off a quarter cycle, and tight secondly reverse power supply a quarter cycle (electric current is-I), last a quarter cycle
Power-off.
The present invention gives " full period " apparent resistivity (ρw(t), Ω M) definition, it is lain in below equation：
WhereinWhen being powered for upper stepThe horizontal component of electric field in direction, it is turned by the electric field of x and y both directions
Change, i.e.,：
Wherein μ0=4 π × 10-7Henry/rice；
Wherein t is time delay (unit：Second), the starting point time delay in each a quarter cycle is disposed as zero；
Wherein erf (x) is that normal probability integrates (or error function)：
Wherein xpFor transition time parameter：
Wherein f (xp) it is normalized function, f (xp) it is on xpMonotonic function：
The present invention can use " intermediate value search approximatioss " to calculate " full period " apparent resistivity ρ according to definition (2)w(t)。
The computational methods of the full period apparent resistivity of time-domain transition level electric field of the present invention comprise the following steps：
(1) the time-domain record material e obtained by field acquisition, is obtainedxAnd e (t)y(t)。
Because physical record is continuously record data, all records according to transmitted waveform, to be divided into the note of signal period
Record (N number of record altogether), and the record in each cycle is divided into the record in four a quarter cycles, each a quarter cycle again
The initial delay time of record is all set to zero.As shown in figure 3, completely to be received in DEEP TRANSIENT ELECTROMAGNETIC SOUNDING a cycle
Waveform (abscissa is the time, and ordinate is electric field), (first a quarter cycle records for step forward, second four
/ mono- cycle recorded for step downward, and the 3rd a quarter cycle records for reversely upper step, the 4th a quarter
Cycle is reversely lower step record.
(2), it is superimposed and filters.
Filtering refers to the influence (using numerical filtering method) for eliminating fixed frequency interference signal (50 or 60Hz civilian electricity).
Superposition refer to the process of to carry out the signal in multiple cycles of collection it is average after linear, additive, it is non-the purpose is to eliminate
The influence of the interference signal of fixed frequency, signal to noise ratio is improved, to obtain the electric field data of high-quality.
The specific additive process of the present invention is separately to be superimposed upper step and lower step, namely：It is upper when forward direction is powered
(the 3rd four points of the transient signal of upper step when the transient signal (in first a quarter cycle) and reverse power supply of step
One of in the cycle) superposition；By the transient signal of the lower step during being powered off after positive power (in second a quarter cycle)
With transient signal (in the 4th a quarter cycle) superposition of the lower step during being powered off after reverse power supply.Specific formula for calculation
For：
(3) editing and processing, is carried out to upper and lower step electric field value.This is the processing procedure of an auxiliary, only records superposition
In " flying spot " using manual type reject.Due to the only poor constant term of upper and lower step electric field value, therefore in editing and processing mistake
Can mutually it be referred in journey.
(4) calculateThe horizontal component of electric field in directionIt can be completed using formula (3).
(5), using " intermediate value search approximatioss " calculating " full period " apparent resistivity ρw(t)。
Specific calculating process is an iterative process.
The first step：Set ρw(t) two extreme value ρw(t)min(such as being taken as 0.0001 Ω Μ) and ρw(t)max(such as
It is taken as 100000 Ω Μ)；
Second step：Calculate both intermediate values：
3rd step：Calculated according to following formulaTheoretical value：
4th step：Calculate relative error error：
If error, which is less than preassigned error, is subjected to limit value (such as 10-20), then ρ noww(t)avgIt can recognize
It is set to the ρ finally requiredw(t).Otherwise, the iterative operation of (the 5th step) in next step is carried out.
5th step：CompareActual value and theoretical value, if：
Then make following displacement, even ρw(t)max=ρw(t)avg, otherwise, make ρw(t)min=ρw(t)avg.Return to second
Step.Formed and solve ρw(t) iterative cycles.
Based on the above method, the invention also provides a kind of meter of the full period apparent resistivity of time-domain transition level electric field
Calculation system, including：
Data acquisition module, for obtaining the time-domain record material of the nagneto-telluric field obtained by field acquisition, obtain two
The electric field component e of individual horizontal directionxAnd e (t)y(t)；
Signal processing module, for electric field component exAnd e (t)y(t) it is overlapped and filtering process, to eliminate interference letter
Number；And for electric field component exAnd e (t)y(t) upper and lower step electric field value carries out editing and processing, rejects flying of deviating considerably from
Point；
Computing module, for calculatingThe horizontal component of electric field in direction
For in transmitting electric dipole
Angle between the line and transmitting electric dipole line of heart point and measuring point；
Further according to formula
Calculate full-time
Section apparent resistivity ρw(t), ρw(t) imply in the formula；
Wherein μ0=4 π × 10-7Henry/rice, t are time delay,It is to receive and dispatch away from erf (x) is that normal probability accumulates
Point or error function：
Wherein xpFor transition time parameter：
Wherein f (xp) it is normalized function, f (xp) it is on xpMonotonic function：
Computing module is calculating full period apparent resistivity ρw(t) when, full-time is calculated using " intermediate value search approximatioss "
Section apparent resistivity ρw(t), specific calculating process is an iterative process.
The signal processing module eliminates the influence of fixed frequency interference signal especially by filtering process；Pass through superposition
Processing, will be average after the signal progress linear, additive in multiple cycles of collection, to eliminate the shadow of the interference signal of on-fixed frequency
Ring.
The time-domain record material of nagneto-telluric field obtained by data collecting module collected is continuously record data, according to hair
Ejected wave shape, all records are divided into the record of signal period, altogether N number of record, and the record in each cycle is divided into four four again
The record in/mono- cycle, the initial delay time of each a quarter periodic recording are all set to zero.
In order to illustrate the superiority of the above method and system in the present invention, to same theory model, using conventional method and
The method of the present invention calculates apparent resistivity respectively.The parameter of wherein theoretical model see the table below 1：
Table 1：The parameter of theoretical model
If Fig. 4 is the electric field wavy curve (abscissa that theoretical model calculates：Time delay；Ordinate：Electric field value).Fig. 5
The early stage calculated for conventional methodAnd late periodApparent resistivity curve (step field in use
Value, early and late apparent resistivity curve is two parallel curves in log-domain.Early stage apparent resistivity curve delay when compared with
Can reflect the essential characteristic of stratum Electrical distribution when early, but it is later in delay when, reflection is artefact, it is impossible to accurately reflects ground
Layer electrical property feature.Late period apparent resistivity curve can reflect the essential characteristic of stratum Electrical distribution when later in delay, but prolong
When more early when late, reflection is artefact, it is impossible to accurately reflects stratum electrical property feature.Both of which is defective)；Fig. 6 is using this hair
The full period apparent resistivity curve (ρ of bright calculatingw(t), it is worth using upper step field)；Fig. 7 is four kinds of apparent resistivity curve (triangles
Symbol is early stage apparent resistivity curve, and closed square is late period apparent resistivity curve, and sphere of movements for the elephants symbol is full period apparent resistance
Rate curve, round dot are mt apparent resistivity curve) comparison, full period apparent resistivity curve delay when it is more early when with
Early stage apparent resistivity curve overlaps, and is overlapped when later in delay with late period apparent resistivity curve, by early stage apparent resistivity curve
The information filling that the transit time section of stratum characteristic can not be accurately reflected with late period apparent resistivity curve is perfect.Therefore full-time
Section apparent resistivity curve either delay when it is more early when, or delay when it is later when, also or in transit time section, equal energy
Correct reflection stratum electrical property feature.Full period apparent resistivity curve is the same as the apparent resistivity curve difference of magnetotelluric sounding method
Be in：The former will tend to geometric sounding apparent resistivity in long delay, and the latter then tends to the resistivity of the bottom.They are all
Stratum electrical property feature correctly and is objectively reflected.It can be seen that the complete apparent resistivity that the present invention calculates the full period solves
The apparent resistivity definition of transit time section and it is continuous the problem of.
It should be appreciated that for those of ordinary skills, can according to the above description be improved or converted,
And all these modifications and variations should all belong to the protection domain of appended claims of the present invention.
Claims (7)
1. a kind of computational methods of the full period apparent resistivity of time-domain transition level electric field, it is characterised in that including following step
Suddenly：
The time-domain record material of nagneto-telluric field obtained by S1, acquisition field acquisition, obtains the electric field point of two horizontal directions
Measure exAnd e (t)y(t)；
S2, to electric field component exAnd e (t)y(t) it is overlapped and filtering process, to eliminate interference signal；
S3, to electric field component exAnd e (t)y(t) upper and lower step electric field value carries out editing and processing, rejects the flying spot deviated considerably from；
S4, calculatingThe horizontal component of electric field in direction
For transmitting electric dipole central point and
Angle between the line and transmitting electric dipole line of measuring point；
S5, further according to formula
Calculating the full period regards
Electricalresistivityρw(t), ρw(t) imply in the formula；Calculating full period apparent resistivity ρw(t) when, using " intermediate value search approaches
Method " calculates full period apparent resistivity ρw(t), specific calculating process is an iterative process；
Wherein μ0=4 π × 10-7Henry/rice, PETo launch the dipole moment of electric dipole, t is time delay,For transmitting-receiving
Away from erf (x) is normal probability integration or error function：
<mrow>
<mi>e</mi>
<mi>r</mi>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>2</mn>
<msqrt>
<mi>&pi;</mi>
</msqrt>
</mfrac>
<msubsup>
<mo>&Integral;</mo>
<mn>0</mn>
<mi>x</mi>
</msubsup>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<msup>
<mi>u</mi>
<mn>2</mn>
</msup>
</mrow>
</msup>
<mi>d</mi>
<mi>u</mi>
<mo>;</mo>
</mrow>
Wherein xpFor transition time parameter：
<mrow>
<msub>
<mi>x</mi>
<mi>p</mi>
</msub>
<mo>=</mo>
<mfrac>
<mi>r</mi>
<msqrt>
<mrow>
<mn>4</mn>
<msub>
<mi>&rho;</mi>
<mi>w</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>/</mo>
<msub>
<mi>&mu;</mi>
<mn>0</mn>
</msub>
</mrow>
</msqrt>
</mfrac>
<mo>;</mo>
</mrow>
Wherein f (xp) it is normalized function, f (xp) it is on xpMonotonic function：
<mrow>
<mi>f</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>p</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<msubsup>
<mi>x</mi>
<mi>p</mi>
<mn>2</mn>
</msubsup>
</mfrac>
<mo>&lsqb;</mo>
<mn>1</mn>
<mo>+</mo>
<mi>e</mi>
<mi>r</mi>
<mi>f</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>p</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mfrac>
<mn>2</mn>
<msqrt>
<mi>&pi;</mi>
</msqrt>
</mfrac>
<msub>
<mi>x</mi>
<mi>p</mi>
</msub>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<msubsup>
<mi>x</mi>
<mi>p</mi>
<mn>2</mn>
</msubsup>
</mrow>
</msup>
<mo>&rsqb;</mo>
<mo>></mo>
<mn>0.</mn>
</mrow>
2. according to the method for claim 1, it is characterised in that be specially in step S5：
The first step：Set ρw(t) two extreme value ρw(t)minAnd ρw(t)max；
Second step：Calculate both intermediate values：
<mrow>
<msub>
<mi>&rho;</mi>
<mi>w</mi>
</msub>
<msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mrow>
<mi>a</mi>
<mi>v</mi>
<mi>g</mi>
</mrow>
</msub>
<mo>=</mo>
<msup>
<mn>10</mn>
<mrow>
<mo>(</mo>
<msub>
<mi>log&rho;</mi>
<mi>w</mi>
</msub>
<msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mi>min</mi>
</msub>
<mo>+</mo>
<msub>
<mi>log&rho;</mi>
<mi>w</mi>
</msub>
<msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>)</mo>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msup>
</mrow>
3rd step：Calculated according to following formulaTheoretical value：
xpmin、xpmax、xpavgThe intermediate quantity calculated for auxiliary, PETo launch the dipole moment of electric dipole；
4th step：Calculate relative error error：
If error, which is less than preassigned error, is subjected to limit value, ρ noww(t)avgIt can regard as what is finally required
ρw(t)；Otherwise, the iterative operation of the 5th step is carried out；
5th step：CompareActual value and theoretical value, if：
Then make following displacement, make ρw(t)max=ρw(t)avg, otherwise, make ρw(t)min=ρw(t)avg, second step is returned to, formation is asked
Solve ρw(t) iterative cycles.
3. according to the method for claim 1, it is characterised in that filtering refers to eliminate fixed frequency interference signal in step S2
Influence；
Superposition refers to the process of to carry out the signal in multiple cycles of collection average after linear, additive, and the purpose is to eliminate on-fixed
The influence of the interference signal of frequency.
4. according to the method for claim 1, it is characterised in that the time-domain of the nagneto-telluric field of collection gained in step S1
Record material is continuously record data, foundation transmitted waveform, all records is divided into the record of signal period, altogether N number of record, and
The record in each cycle is divided into the record in four a quarter cycles again, during the initial delay of each a quarter periodic recording
Between be all set to zero.
A kind of 5. computing system of the full period apparent resistivity of time-domain transition level electric field, it is characterised in that including：
Data acquisition module, for obtaining the time-domain record material of the nagneto-telluric field obtained by field acquisition, obtain two water
Square to electric field component exAnd e (t)y(t)；
Signal processing module, for electric field component exAnd e (t)y(t) it is overlapped and filtering process, to eliminate interference signal；
And for electric field component exAnd e (t)y(t) upper and lower step electric field value carries out editing and processing, rejects the flying spot deviated considerably from；
Computing module, for calculatingThe horizontal component of electric field in direction
For transmitting electric dipole central point and
Angle between the line and transmitting electric dipole line of measuring point；
Further according to formula
Calculating the full period regards
Electricalresistivityρw(t), ρw(t) imply in the formula；
Wherein μ0=4 π × 10-7Henry/rice, PETo launch the dipole moment of electric dipole, t is time delay,For transmitting-receiving
Away from erf (x) is normal probability integration or error function：
<mrow>
<mi>e</mi>
<mi>r</mi>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>2</mn>
<msqrt>
<mi>&pi;</mi>
</msqrt>
</mfrac>
<msubsup>
<mo>&Integral;</mo>
<mn>0</mn>
<mi>x</mi>
</msubsup>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<msup>
<mi>u</mi>
<mn>2</mn>
</msup>
</mrow>
</msup>
<mi>d</mi>
<mi>u</mi>
<mo>;</mo>
</mrow>
Wherein xpFor transition time parameter：
<mrow>
<msub>
<mi>x</mi>
<mi>p</mi>
</msub>
<mo>=</mo>
<mfrac>
<mi>r</mi>
<msqrt>
<mrow>
<mn>4</mn>
<msub>
<mi>&rho;</mi>
<mi>w</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>/</mo>
<msub>
<mi>&mu;</mi>
<mn>0</mn>
</msub>
</mrow>
</msqrt>
</mfrac>
<mo>;</mo>
</mrow>
Wherein f (xp) it is normalized function, f (xp) it is on xpMonotonic function：
The computing module regards in the calculating full period
Electricalresistivityρw(t) when, full period apparent resistivity ρ is calculated using " intermediate value search approximatioss "w(t), specific calculating process changes for one
For process.
6. system according to claim 5, it is characterised in that the signal processing module disappears especially by filtering process
Except the influence of fixed frequency interference signal；By overlap-add procedure, put down after the signal in multiple cycles of collection is carried out into linear, additive
, to eliminate the influence of the interference signal of on-fixed frequency.
7. system according to claim 5, it is characterised in that the nagneto-telluric field obtained by the data collecting module collected
Time-domain record material be continuously record data, according to transmitted waveform, all record is divided into the record of signal period, common N
Individual record, and the record in each cycle is divided into the record in four a quarter cycles again, each a quarter periodic recording
Initial delay time is all set to zero.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510433717.XA CN105093341B (en) | 2015-07-22 | 2015-07-22 | The computational methods and system of the full period apparent resistivity of time-domain transition level electric field |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510433717.XA CN105093341B (en) | 2015-07-22 | 2015-07-22 | The computational methods and system of the full period apparent resistivity of time-domain transition level electric field |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105093341A CN105093341A (en) | 2015-11-25 |
| CN105093341B true CN105093341B (en) | 2017-12-12 |
Family
ID=54574195
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510433717.XA Active CN105093341B (en) | 2015-07-22 | 2015-07-22 | The computational methods and system of the full period apparent resistivity of time-domain transition level electric field |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105093341B (en) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110412656B (en) * | 2019-07-18 | 2021-05-04 | 长江大学 | Magnetotelluric sounding data time domain noise suppression method and system |
| CN112415606B (en) * | 2020-10-16 | 2022-12-13 | 山东大学 | Resistivity detection method and system for urban environment non-contact arbitrary pole arrangement |
| CN113671585B (en) * | 2021-08-18 | 2022-04-29 | 中国矿业大学 | Intelligent transient electromagnetic detection and real-time early warning method for excavation roadway |
| CN116400415B (en) * | 2023-04-12 | 2023-11-24 | 中国地质科学院地球物理地球化学勘查研究所 | Transient electromagnetic global apparent resistivity defining method based on attenuation process |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7538555B2 (en) * | 2003-11-05 | 2009-05-26 | Shell Oil Company | System and method for locating an anomaly ahead of a drill bit |
| EP1889200A4 (en) * | 2005-06-09 | 2017-08-16 | Exxonmobil Upstream Research Company | Method for determining earth vertical electrical anisotropy in marine electromagnetic surveys |
| CN101650443B (en) * | 2009-09-15 | 2012-03-21 | 重庆大学 | Back-propagation network calculating method of apparent resistivity |
| CN102062876A (en) * | 2010-11-17 | 2011-05-18 | 汤井田 | Electrical sounding method for whole-region couple source frequency domain |
| CN104360401B (en) * | 2014-12-10 | 2016-10-12 | 中国科学院地质与地球物理研究所 | A kind of transient electromagnetic B field descends objective body geological information method definitely |
| CN104777516A (en) * | 2015-04-15 | 2015-07-15 | 国网重庆市电力公司电力科学研究院 | Apparent resistivity calculating method on basis of non-linear equation solution modular form |
-
2015
- 2015-07-22 CN CN201510433717.XA patent/CN105093341B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN105093341A (en) | 2015-11-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105093341B (en) | The computational methods and system of the full period apparent resistivity of time-domain transition level electric field | |
| CN108873083B (en) | A kind of artificial field source frequency domain electromagnetism apparent resistivity measurement method | |
| CN102373923B (en) | Reservoir stratum identification method | |
| CN103064120B (en) | Underground coal mine magnetoelectricity composite geophysical methods | |
| CN101813725A (en) | Method for measuring phase difference of low-frequency signals | |
| CN103733089B (en) | For including the system and method that the underground of uncertainty estimation characterizes | |
| CN105093319B (en) | Ground micro-seismic static correcting method based on 3D seismic data | |
| CN106944719B (en) | Posture of welding torch recognition methods based on multi-pole capacitance sensor | |
| CN203149134U (en) | Ground penetrating radar for geological exploration | |
| CN104483705A (en) | Three-dimensional residual static correction method | |
| CN105445789A (en) | Three-dimensional Fresnel volume travel-time tomographic method based on multiple reflected refraction wave constraint | |
| CN104133245A (en) | Seismic data static correction method and system | |
| CN107345484B (en) | A kind of high seam lower leaf " mixing top plate " integrality detection and its control method | |
| CN104122592B (en) | Time shift magnetotelluric signal acquisition and inversion method | |
| CN104375172B (en) | Volcanic underlayer structural configuration correct imaging method | |
| CN106249297A (en) | Fracturing microseism seismic source location method and system based on Signal estimation | |
| CN107884825B (en) | Uncertainty modeling method based on seismic multi-attribute | |
| CN104155694B (en) | A kind of residual static corrections reflecting converted shear wave common geophone stack section | |
| CN109343134A (en) | A kind of Transient Electromagnetic Method in Mine data analysis interpretation method and system | |
| CN206177406U (en) | Real -time tracking measurement device of ground water level | |
| CN103969685B (en) | A kind of processing method of thin interbed seismic signal | |
| CN108121001A (en) | The quasi-static scene positioning accuracy optimization method differentiated based on fixed solution continuity | |
| CN202033478U (en) | Very low frequency electromagnetic instrument with orthogonal antenna | |
| CN104375181A (en) | Method for rapidly generating three-dimensional common-detection wave point first-break gather | |
| CN105425305B (en) | The method that electric layers resolution capability is improved using step response under time-domain transition level electric field |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |