CN110187394A - Double field source electromagnetic depth methods obtain the anisotropic method and device of formation resistivity - Google Patents
Double field source electromagnetic depth methods obtain the anisotropic method and device of formation resistivity Download PDFInfo
- Publication number
- CN110187394A CN110187394A CN201910536766.4A CN201910536766A CN110187394A CN 110187394 A CN110187394 A CN 110187394A CN 201910536766 A CN201910536766 A CN 201910536766A CN 110187394 A CN110187394 A CN 110187394A
- Authority
- CN
- China
- Prior art keywords
- source
- signal
- field
- frequency
- stratum
- 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/08—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with magnetic or electric fields produced or modified by objects or geological structures or by detecting devices
- G01V3/083—Controlled source electromagnetic [CSEM] surveying
-
- 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
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Electromagnetism (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
- Measurement Of Resistance Or Impedance (AREA)
Abstract
The present invention relates to a kind of double field source electromagnetic depth methods to obtain the anisotropic method and device of formation resistivity, which comprises is surveying outside area, is laying two field sources along orthogonal both direction；Two field source substep alternate emissions, two irrelevant signals；Record the current time sequence of two field sources transmitting；Horizontal component of electric field response signal orthogonal on each measuring point is acquired surveying in area, and/or the vertical magnetic field response signal on each measuring point of acquisition, and is recorded as time series；According to the current time sequence and equipment parameter and electromagnetic field response signal time series, the APPARENT RESISTIVITY for surveying stratum below area's measuring point is obtained；It is calculated further according to the APPARENT RESISTIVITY and surveys area's formation resistivity anisotropy coefficient.The method of the invention observational components are few, and treatment process is easy；Strong antijamming capability, measurement accuracy is high, and the frequency point of the APPARENT RESISTIVITY of acquisition is abundant, has higher resolution ratio to stratum, and can improve working efficiency.
Description
Technical field
The present invention relates to technical field of geophysical exploration, and in particular to a kind of double field source electromagnetic depth methods acquisitions stratum electricity
The anisotropic method and device of resistance rate.
Background technique
Due to the shallow earth's surface media components uneven distribution of the earth, all kinds of tomographies and construction, ground are formd under geologic process
Lower complex geologic body all has three-dimensional feature, and has anisotropy.In electromagnetic prospecting field, can be obtained by surface exploration
Underground geologic bodies anisotropic character is taken, technical parameter can be provided for the fine inverting of electromagnetic sounding data.But it can accurately obtain at present
Take the anisotropic observation method of earth resistivity few.Five component magnetotelluric sounding methods (MT, AMT) are by observing orthogonal water
Ordinary telegram field (Ex, Ey) and horizontal magnetic field (Hx, Hy) and perpendicular magnetic component (Hz), estimate two orthogonal electrical main shafts
Apparent resistivity can be used to analytically resistivity anisotropy feature.But since natural field signal is faint, estimate on some frequency points
It is inadequate to count precision.
In addition, existing artificial source's electromagnetic depth method exploration engineering is needed using multiple field sources or tensor controllable source electromagnetic sounding
Method frequency variation method obtains earth resistivity anisotropic character, and complex disposal process, working efficiency is lower, and anti-interference ability is not strong, special
It is not in mining area operation, it is very difficult to collect quality datas.
Summary of the invention
In view of this, it is an object of the invention to overcome the deficiencies of the prior art and provide a kind of double field source electromagnetic depth methods
Obtain the anisotropic method and device of formation resistivity.
In order to achieve the above object, the present invention adopts the following technical scheme: a kind of double field source electromagnetic depth methods obtain stratum electricity
The anisotropic method of resistance rate, comprising:
It is surveying outside area, is laying two field sources: source x and source y along orthogonal x-axis and y-axis both direction；
Described two field source substep alternate emissions go out two irrelevant signals；
The current time sequence Ix (t), Iy (t) of record source x and source y transmitting；
It is surveying in area, is acquiring horizontal component of electric field response signal Ex and Ey orthogonal on each measuring point, and/or, it acquires on each measuring point
Vertical magnetic field response signal Hz, and collected response signal is recorded as time series Ex (t) and Ey (t) and/or Hz
(t)；
According to the current time sequence and equipment parameter and electromagnetic field response signal time series, by following
Ring cross-correlation method obtains the APPARENT RESISTIVITY for surveying area stratum；
Resistivity anisotropy coefficient is calculated according to the APPARENT RESISTIVITY for surveying area stratum.
Optionally, described two field source substep alternate emissions go out two irrelevant signals, specifically include:
Source x and source y needs substep alternately to excite twice；When source x sends the first signal, source y sends second signal；Source x hair
When sending second signal, source y sends the first signal；Source x and source y substep alternate emission the first signal and the second signal, and described
One signal and the second signal are irrelevant.
Optionally, first signal is single-frequency square-wave signal, and the second signal is that the inverse M sequence pseudorandom that repeats is believed
Number, the frequency of the single-frequency square-wave signal is lower than the inverse low-limit frequency for repeating M sequence pseudo-random signal.
Optionally, described according to the current time sequence and equipment parameter and the electromagnetic field response signal time
Sequence obtains the APPARENT RESISTIVITY for surveying area stratum by recycling cross-correlation method, comprising:
The current time sequence Ix (t) and Ex (t) or Hz (t) time series that source x is emitted are matched, the electricity that source y is emitted
It flows time series Iy (t) and Ey (t) or Hz (t) time series is matched, then the survey area stratum is obtained by circulation cross-correlation method
Frequency response；
According to the frequency response for surveying area stratum and equipment parameter, excitation current and receive electromagnetic field component
Time series, iterative method obtain the APPARENT RESISTIVITY for surveying area stratum.
Optionally, anisotropy coefficient is calculated according to the APPARENT RESISTIVITY on the stratum to be measured by formula (1),
λ=ρx(wi)/ρy(wi) formula (1)
Wherein, λ is anisotropy coefficient, ρx(wi) and ρy(wi) it is wiApparent resistivity under frequency in x and y both direction.
Optionally, source x is approximately equal to source y at a distance from survey district center and surveys at a distance from district center, and each field source and survey area
The value range of the distance at center is 4km~20km.
Optionally, the stratum includes: land and ocean.
The present invention also provides a kind of double field source electromagnetic depth methods to obtain the anisotropic device of formation resistivity, comprising:
First field source and the second field source go out two irrelevant signals for substep alternate emission；Wherein, described first
Field source and the second field source are surveying the orthogonal x-axis in area outer and the laying of y-axis both direction；
Electric current timing logging modle, for recording the current time sequence Ix of first field source and the transmitting of the second field source
(t), Iy (t)；
Acquisition module is used to survey in area, acquires horizontal component of electric field response signal Ex and Ey orthogonal on each measuring point, and/or,
The vertical magnetic field response signal Hz on each measuring point is acquired, and collected response signal is recorded as time series Ex (t) and Ey
(t) and/or Hz (t)；
Module is obtained, for according to the current time sequence and equipment parameter, excitation current and reception electromagnetism
Field component time series obtains the APPARENT RESISTIVITY for surveying area stratum by recycling cross-correlation method；
Computing module surveys area's formation resistivity respectively to different for calculating according to the APPARENT RESISTIVITY for surveying area stratum
Property.
Optionally, the acquisition module includes:
Unpolarizable electrode and electric field recorder, for acquiring horizontal component of electric field response signal Ex and Ey orthogonal on each measuring point；
And/or magnetic field acquisition module, for acquiring the vertical magnetic field response signal Hz on each measuring point；
Wherein, the magnetic field acquisition module includes magnetograph and broadband induction type magnetic sensor or fluxgate
Sensor or high-temperature superconducting magnetometer.
Optionally, first field source and the second field source substep alternately excite twice: first field source sends the first letter
Number when, second field source sends second signal；When first field source sends second signal, second field source sends first
Signal；
Wherein, first signal is single-frequency square-wave signal, and the second signal is inverse repetition M sequence pseudo-random signal,
The frequency of the single-frequency square-wave signal is lower than the inverse low-limit frequency for repeating M sequence pseudo-random signal.
The invention adopts the above technical scheme, and double field source electromagnetic depth methods obtain the anisotropic side of formation resistivity
Method includes: to survey outside area, lays two field sources, source x and source y along orthogonal x-axis and y-axis both direction；Described two field sources point
Step alternate emission goes out two irrelevant signals；The current time sequence Ix (t), Iy (t) of record source x and source y transmitting；It is surveying
In area, horizontal component of electric field response signal Ex and Ey orthogonal on each measuring point are acquired, and/or, the vertical magnetic field acquired on each measuring point is rung
Induction signal Hz, and collected response signal is recorded as time series Ex (t) and Ey (t) and/or Hz (t)；According to the electricity
It flows time series and equipment parameter and receives electromagnetic field component time series, by described in circulation cross-correlation method acquisition
Survey the APPARENT RESISTIVITY on area stratum；It is calculated according to the APPARENT RESISTIVITY for surveying area stratum and surveys area's formation resistivity respectively to different
Property coefficient.The method of the invention feature irrelevant using single-frequency square-wave signal and inverse repetition M sequence pseudo-random signal, and
Electromagnetic prospecting is carried out using circulation cross-correlation method, only acquire the horizontal component of electric field (Ex, Ey) of both direction on measuring point or is only measured
Vertical magnetic field (Hz) can obtain the APPARENT RESISTIVITY for surveying area stratum, survey area's formation resistivity anisotropy system to calculate
Number.Method of the present invention reduces observational components, and treatment process is easier；Due to inverse repetition M sequence pseudo-random signal
Power spectrum is line spectrum, bandwidth and uncorrelated to Hz noise, therefore this method has strong anti-interference ability, measurement essence
Degree is high, and the frequency point of the APPARENT RESISTIVITY of acquisition is abundant, has higher resolution ratio to stratum, and can improve working efficiency.In addition,
This method carries out field source excitation simultaneously from two orthogonal directions, is more advantageous to the electromagnetic response for obtaining deep formation.
Detailed description of the invention
In order to more clearly explain the embodiment of the invention or the technical proposal in the existing technology, to embodiment or will show below
There is attached drawing needed in technical description to be briefly described, it should be apparent that, the accompanying drawings in the following description is only this
Some embodiments of invention for those of ordinary skill in the art without creative efforts, can be with
It obtains other drawings based on these drawings.
Fig. 1 is that double field source electromagnetic depth methods of the invention obtain the anisotropic offer of embodiment of the method one of formation resistivity
Flow diagram；
Fig. 2 is that schematic layout pattern is detected involved in the embodiment of the present invention one；
Fig. 3 (a) is the inverse waveform diagram for repeating M sequence pseudo-random signal described in the embodiment of the present invention one；
Fig. 3 (b) is the inverse power spectrum for repeating M sequence pseudo-random signal described in the embodiment of the present invention one；
Fig. 3 (c) is the inverse auto-correlation function for repeating M sequence pseudo-random signal described in the embodiment of the present invention one；
Fig. 4 is that single-frequency square-wave signal described in the embodiment of the present invention one is covered with the inverse M sequence pseudo-random signal that repeats
Frequency spectrum；
Fig. 5 is that double field source electromagnetic depth methods of the invention obtain the anisotropic offer of Installation practice one of formation resistivity
Structural schematic diagram.
In figure: 1, the first field source；2, the second field source；3, electric current timing logging modle；4, acquisition module；5, module is obtained；
6, computing module.
Specific embodiment
To make the object, technical solutions and advantages of the present invention clearer, technical solution of the present invention will be carried out below
Detailed description.Obviously, described embodiments are only a part of the embodiments of the present invention, instead of all the embodiments.Base
Embodiment in the present invention, those of ordinary skill in the art are obtained all without making creative work
Other embodiment belongs to the range that the present invention is protected.
Fig. 1 is that double field source electromagnetic depth methods of the invention obtain the anisotropic offer of embodiment of the method one of formation resistivity
Flow diagram.
As depicted in figs. 1 and 2, method described in the present embodiment includes:
S11: surveying outside area, lays two field sources, source x and source y along orthogonal x-axis and y-axis both direction；
Further, source x and source y is two orthogonal electric dipole sources；Source x is approximately equal to source y and is surveyed at a distance from survey district center
The distance of district center, and value range of each field source at a distance from survey district center is 4km~20km.
S12: described two field source substep alternate emissions go out two irrelevant signals；
Further, described two field source substep alternate emissions go out two irrelevant signals, specifically include:
Source x and source y needs substep alternately to excite twice: when source x sends the first signal, source y sends second signal；Source x hair
When sending second signal, source y sends the first signal；Source x and source y substep alternate emission the first signal and the second signal, and described
One signal and the second signal are irrelevant.
Specifically, as shown in Figure 3 and Figure 4, first signal is single-frequency square-wave signal, the second signal is inverse repetition
M sequence pseudo-random signal, the frequency of the single-frequency square-wave signal is lower than the inverse lowest frequency for repeating M sequence pseudo-random signal
Rate.
The current time sequence Ix (t), Iy (t) of S13: record source x and source y transmitting；
S14: surveying in area, acquiring horizontal component of electric field response signal Ex and Ey orthogonal on each measuring point, and/or, acquire each survey
Vertical magnetic field response signal Hz on point, and collected response signal is recorded as time series Ex (t) and Ey (t), and/or
Hz(t)；
S15: according to the current time sequence and equipment parameter and electromagnetic field response signal time series, lead to
It crosses circulation cross-correlation method and obtains the APPARENT RESISTIVITY for surveying area stratum；
Further, it is described according to the current time sequence and equipment parameter and electromagnetic field response signal when
Between sequence, by recycle cross-correlation method obtain it is described survey area stratum APPARENT RESISTIVITY, comprising:
The current time sequence Ix (t) and Ex (t) or Hz (t) time series that source x is emitted are matched, the electricity that source y is emitted
It flows time series Iy (t) and Ey (t) or Hz (t) time series is matched, then the survey area stratum is obtained by circulation cross-correlation method
Frequency response；
According to the frequency response for surveying area stratum and equipment parameter and electromagnetic field component is received, described in acquisition
Survey the APPARENT RESISTIVITY ρ on area stratumx(wi) and ρy(wi)。
It is described to obtain the frequency response for surveying area stratum by recycling cross-correlation method, further according to the frequency for surveying area stratum
Rate response and equipment parameter and reception electromagnetic field component obtain the APPARENT RESISTIVITY for surveying area stratum, the processing
Specific calculating process is referring to calculating process disclosed in existing literature (patent No.: CN201610013962.X).
Wherein, the equipment parameter includes: the two poles of the earth spacing of horizontal electric dipole, the electrode spacing of receiving electrode, connects
The midpoint and the distance between the two poles of the earth midpoint of horizontal electric dipole for receiving electrode, with the angle of horizontal electric dipole the two poles of the earth perpendicular bisector, with
And the parameters such as frequency of transmitting signal.
S16: resistivity anisotropy coefficient is calculated according to the APPARENT RESISTIVITY for surveying area stratum.
Further, be in step S16 by formula (1) according to the APPARENT RESISTIVITY on the stratum to be measured calculate it is each to
Anisotropic coefficient,
λ=ρx(wi)/ρy(wi) formula (1)
Wherein, λ is anisotropy coefficient, ρx(wi) and ρy(wi) it is wiApparent resistivity under frequency in x and y both direction.
Method described in the present embodiment is being surveyed outside area, is being laid along orthogonal x-axis and y-axis both direction when actually executing
Two field sources: source x and source y, described two field source substep alternate emissions go out two irrelevant signals.Surveying each measuring point in area
Receiving electrode is laid along x and y both direction in position, for acquiring two orthogonal horizontal electric field components of Ex and Ey, or can be with
With broadband induction type magnetic sensor measurement vertical magnetic field Hz, or horizontal component of electric field Ex, Ey and vertical magnetic field Hz are measured simultaneously, with
Upper three kinds of modes can select measurement electric field or magnetic field signal according to actual place grounding requirement.Acquisition instrument is with identical sampling
Ex and Ey horizontal component of electric field or vertical magnetic field Hz component recording are time series Ex (t) and Ey (t) or Hz (t) by rate, and use phase
The source x and source y current signal sent is recorded as current time sequence Ix (t) respectively with identical sample rate by same acquisition instrument
With Iy (t).
Source x and source y are worked at the same time, and send single-frequency square-wave signal and inverse repetition M sequence pseudo-random signal, two field sources respectively
Need substep alternating excitation signal twice: the first step is that source x sends single-frequency square wave, and source y, which sends the inverse M sequence pseudorandom that repeats, to be believed
Number；Then second step is carried out, when source x sends inverse repetition M pseudo-random signal, source y sends single-frequency square wave.Two field source substeps are alternately
Emit two kinds of signal, field source transmitting signal and acquisition instrument acquisition signal are to carry out both ends synchronized sampling by satellite.
Horizontal component of electric field and magnetic field at measuring point are the stack responses on the stratum under exciting jointly to source x and source y.Respectively by source
The current time sequence Ix (t) and Ex (t) or Hz (t) time series of x transmitting are matched, the current time sequence Iy that source y is emitted
(t) matched with Ey (t) or Hz (t) time series, then by circulation cross-correlation method, using existing literature (patent No.:
CN201610013962.X the calculation method disclosed in) calculates measuring point both direction APPARENT RESISTIVITY ρx(wi) and ρy(wi),
Recycle identical frequency wiLower x and y both direction apparent resistivity ρx(wi) and ρy(wi) calculate anisotropy coefficient.
The working foundation of the present embodiment the method is: when Ix is single-frequency square-wave signal electric current, Iy is that inverse repetition M sequence is pseudo-
Random signal electric current；When Ix is inverse repetition M sequence pseudo-random signal electric current, Iy is single-frequency square-wave signal electric current；As shown in figure 3,
The inverse power spectrum for repeating M sequence is line spectrum, single-frequency square-wave signal and it is inverse repeat M sequence pseudo-random signal be it is irrelevant,
That is, Iy is mono- when making circulation cross-correlation with the inverse repetition M sequence pseudo-random signal electric current of Ix and electric field Ex or vertical magnetic field Hz
The response of frequency square wave current is pressed；Equally, make circulation cross-correlation with Ix single-frequency square wave current and electric field Ex or vertical magnetic field Hz
When, the inverse response for repeating M sequence pseudo-random signal electric current of Iy is pressed；Pass through the excitation twice of two field sources, and acquisition mould
Block receives corresponding response signal in measuring point, so that it may obtain the APPARENT RESISTIVITY ρ that x and y both direction is calculated by horizontal component of electric fieldx
(wi) and ρy(wi), or the APPARENT RESISTIVITY ρ calculated by measurement vertical magnetic fieldx(wi) and ρy(wi), wherein frequency range packet
Include single-frequency square wave frequency and the inverse frequency for repeating M sequence pseudo-random signal and being covered.
In actual operation, the inverse repetition M sequence is arranged according to the difference of place grounding requirement and medium in engineering staff
The sequence order and corresponding frequencies of pseudo-random signal, when the inverse sequence order for repeating M sequence pseudo-random signal and frequency are true
After fixed, then the frequency of single-frequency square-wave signal is set, the frequency of the single-frequency square-wave signal repeats M sequence pseudo-random signal lower than inverse
Low-limit frequency, and the frequency of single-frequency square-wave signal is adjustable.
It should be noted that stratum described in the present embodiment includes land and ocean stratum；It is announced in the present embodiment
It about the measurement of vertical magnetic field, is not limited to that aircraft also can be used and measure in the sky in ground observation, or uses submariner device
It is measured in marine environment, or electromagnetic field receivers is launched and measure vertical magnetic field in seabed.Equally, field source excitation can also be with
It is laid in ocean, carries out marine electromagnetic detection.
In addition, in the technical solution that the present embodiment is announced about signal behavior be not limited to it is inverse repeat M sequence puppet with
The incoherent single-frequency square-wave signal of machine signal, other also can be used with the inverse incoherent signal of M sequence pseudo-random signal that repeats.
The present embodiment the method spy irrelevant using single-frequency square-wave signal and inverse repetition M sequence pseudo-random signal
Point, and electromagnetic prospecting is carried out using circulation cross-correlation method, only acquire the horizontal component of electric field (Ex, Ey) or only of both direction on measuring point
Measurement vertical magnetic field (Hz) can obtain the APPARENT RESISTIVITY for surveying area stratum, survey area's formation resistivity respectively to different to calculate
Property coefficient.Method described in the present embodiment reduces observational components, and treatment process is easier；M sequence pseudorandom is repeated due to inverse
The power spectrum of signal is line spectrum, bandwidth and uncorrelated to Hz noise, therefore this method has strong anti-interference ability,
Measurement accuracy is high, and the frequency point of the APPARENT RESISTIVITY of acquisition is abundant, has higher resolution ratio to stratum, and can improve work effect
Rate.In addition, this method carries out field source excitation simultaneously from two orthogonal directions, it is more advantageous to the electromagnetic response for obtaining deep formation.
Fig. 5 is that double field source electromagnetic depth methods of the invention obtain the anisotropic offer of Installation practice one of formation resistivity
Structural schematic diagram.
As shown in figure 5, device described in the present embodiment, comprising:
First field source 1 and the second field source 2 go out two irrelevant signals for substep alternate emission；Wherein, described
One field source and the second field source are surveying the orthogonal x-axis in area outer and the laying of y-axis both direction；
Electric current timing logging modle 3, for recording the current time sequence Ix of first field source and the transmitting of the second field source
(t), Iy (t)；
Acquisition module 4 acquires horizontal component of electric field response signal Ex and Ey orthogonal on each measuring point for surveying in area, and/
Or, acquire the vertical magnetic field response signal Hz on each measuring point, and by collected response signal be recorded as time series Ex (t) and
Ey (t) and/or Hz (t)；
Module 5 is obtained, when for according to the current time sequence and equipment parameter and electromagnetic field response signal
Between sequence, by recycle cross-correlation method obtain it is described survey area stratum APPARENT RESISTIVITY；
Computing module 6, for calculating resistivity anisotropy according to the APPARENT RESISTIVITY for surveying area stratum.
Further, the acquisition module 4 includes:
Unpolarizable electrode and electric field recorder, for acquiring horizontal component of electric field response signal Ex and Ey orthogonal on each measuring point；
And/or magnetic field acquisition module, for acquiring the vertical magnetic field response signal Hz on each measuring point；
Wherein, the magnetic field acquisition module includes magnetograph and broadband induction type magnetic sensor or fluxgate
Sensor or high-temperature superconducting magnetometer.
Further, first field source and the second field source substep alternately excite twice: first field source sends first
When signal, second field source sends second signal；When first field source sends second signal, second field source sends the
One signal；
Wherein, first signal is single-frequency square-wave signal, and the second signal is inverse repetition M sequence pseudo-random signal,
The frequency of the single-frequency square-wave signal is lower than the inverse low-limit frequency for repeating M sequence pseudo-random signal.
Device described in the present embodiment in actual use, as shown in Fig. 2, first survey area outside, first field source (source
X) it lays along the x-axis direction, (source y) is laid second field source along the y-axis direction, and source x and source y are worked at the same time, and sends single-frequency respectively
Square-wave signal and inverse repetition M sequence pseudo-random signal, two field sources need substep alternating excitation signal twice: the first step is source x
Single-frequency square wave is sent, source y sends inverse repetition M sequence pseudo-random signal；Then second step is carried out, source x sends inverse repetition M pseudorandom
When signal, source y sends single-frequency square wave.Two two kinds of signals of field source substep alternate emission, the electric current timing record mould
Block 3 and the acquisition module 4 send and receive by satellite time signals the synchronized sampling at both ends.
The acquisition module respectively matches current time sequence Ix (t) and Ex (t) or Hz (t) time series that source x emits
Right, the current time sequence Iy (t) and Ey (t) or Hz (t) time series that source y is emitted are matched, then by recycling cross-correlation method,
Measuring point APPARENT RESISTIVITY ρ is calculated using calculation method disclosed in existing literature (patent No.: CN201610013962.X)x
(wi) and ρy(wi)；Identical frequency w is finally utilized by the computing moduleiLower x and y both direction apparent resistivity ρx(wi) and ρy
(wi) calculate anisotropy coefficient.
The working principle of the present embodiment described device and the working principle of method as described above are identical, and this will not be detailed here.
It should be noted that stratum described in the present embodiment includes but is not limited to land and ocean stratum；In the present embodiment
The device can also be laid in ocean in addition to that can detect to land, carry out marine electromagnetic detection.
In addition, in the technical solution that the present embodiment is announced about signal behavior be not limited to it is inverse repeat M sequence puppet with
The incoherent single-frequency square-wave signal of machine signal, other also can be used with the inverse incoherent signal of M sequence pseudo-random signal that repeats.
Device of the present invention reduces observational components, and device is lighter；M sequence pseudo-random signal is repeated due to inverse
Power spectrum be line spectrum, bandwidth and uncorrelated to Hz noise, thus the collected response signal of the device have it is relatively strong
Anti-interference ability, measurement accuracy is high, and the frequency point of the APPARENT RESISTIVITY of acquisition is abundant, have higher resolution ratio to stratum, and
Working efficiency can be improved.Swash in addition, the present apparatus carries out field source from two orthogonal directions by the first field source and the second field source simultaneously
Hair is more advantageous to the electromagnetic response for obtaining deep formation.
It is understood that same or similar part can mutually refer in the various embodiments described above, in some embodiments
Unspecified content may refer to the same or similar content in other embodiments.
It should be noted that in the description of the present invention, term " first ", " second " etc. are used for description purposes only, without
It can be interpreted as indication or suggestion relative importance.In addition, in the description of the present invention, unless otherwise indicated, the meaning of " multiple "
Refer at least two.
Any process described otherwise above or method description are construed as in flow chart or herein, and expression includes
It is one or more for realizing specific logical function or process the step of executable instruction code module, segment or portion
Point, and the range of the preferred embodiment of the present invention includes other realization, wherein can not press shown or discussed suitable
Sequence, including according to related function by it is basic simultaneously in the way of or in the opposite order, Lai Zhihang function, this should be of the invention
Embodiment person of ordinary skill in the field understood.
It should be appreciated that each section of the invention can be realized with hardware, software, firmware or their combination.Above-mentioned
In embodiment, software that multiple steps or method can be executed in memory and by suitable instruction execution system with storage
Or firmware is realized.It, and in another embodiment, can be under well known in the art for example, if realized with hardware
Any one of column technology or their combination are realized: having a logic gates for realizing logic function to data-signal
Discrete logic, with suitable combinational logic gate circuit specific integrated circuit, programmable gate array (PGA), scene
Programmable gate array (FPGA) etc..
Those skilled in the art are understood that realize all or part of step that above-described embodiment method carries
It suddenly is that relevant hardware can be instructed to complete by program, the program can store in a kind of computer-readable storage medium
In matter, which when being executed, includes the steps that one or a combination set of embodiment of the method.
It, can also be in addition, each functional unit in each embodiment of the present invention can integrate in a processing module
It is that each unit physically exists alone, can also be integrated in two or more units in a module.Above-mentioned integrated mould
Block both can take the form of hardware realization, can also be realized in the form of software function module.The integrated module is such as
Fruit is realized and when sold or used as an independent product in the form of software function module, also can store in a computer
In read/write memory medium.
Storage medium mentioned above can be read-only memory, disk or CD etc..
In the description of this specification, reference term " one embodiment ", " some embodiments ", " example ", " specifically show
The description of example " or " some examples " etc. means specific features, structure, material or spy described in conjunction with this embodiment or example
Point is included at least one embodiment or example of the invention.In the present specification, schematic expression of the above terms are not
Centainly refer to identical embodiment or example.Moreover, particular features, structures, materials, or characteristics described can be any
One or more embodiment or examples in can be combined in any suitable manner.
Although the embodiments of the present invention has been shown and described above, it is to be understood that above-described embodiment is example
Property, it is not considered as limiting the invention, those skilled in the art within the scope of the invention can be to above-mentioned
Embodiment is changed, modifies, replacement and variant.
Claims (10)
1. a kind of double field source electromagnetic depth methods obtain the anisotropic method of formation resistivity characterized by comprising
It is surveying outside area, is laying two field sources: source x and source y along orthogonal x-axis and y-axis both direction；
Described two field source substep alternate emissions go out two irrelevant signals；
The current time sequence Ix (t), Iy (t) of record source x and source y transmitting；
It is surveying in area, is acquiring horizontal component of electric field response signal Ex and Ey orthogonal on each measuring point, and/or, acquire hanging down on each measuring point
Straight magnetic responsiveness signal Hz, and collected response signal is recorded as time series Ex (t) and Ey (t) and/or Hz (t)；
It is mutual by recycling according to the current time sequence and equipment parameter and electromagnetic field response signal time series
Correlation method obtains the APPARENT RESISTIVITY for surveying area stratum；
It is calculated according to the APPARENT RESISTIVITY for surveying area stratum and surveys area's formation resistivity anisotropy coefficient.
2. the method according to claim 1, wherein described two field source substep alternate emissions go out two mutual not phases
The signal of pass, specifically includes:
Source x and source y needs substep alternately to excite twice；When source x sends the first signal, source y sends second signal；Source x sends the
When binary signal, source y sends the first signal；Source x and source y substep alternate emission the first signal and the second signal, and first letter
It is number irrelevant with the second signal.
3. according to the method described in claim 2, it is characterized in that, first signal be single-frequency square-wave signal, described second
Signal is inverse repetition M sequence pseudo-random signal, and the frequency of the single-frequency square-wave signal is lower than the inverse repetition M sequence pseudorandom letter
Number low-limit frequency.
4. the method according to claim 1, wherein described join according to the current time sequence and equipment
Several and electromagnetic field response signal time series obtains the APPARENT RESISTIVITY for surveying area stratum by recycling cross-correlation method,
Include:
Current time sequence Ix (t) that source x is emitted and Ex (t) or Hz (t) time series are matched, when the electric current that source y is emitted
Between sequence Iy (t) and Ey (t) or Hz (t) time series match, then the frequency on the survey area stratum is obtained by recycling cross-correlation method
Rate response；
According to the frequency response for surveying area stratum and equipment parameter, excitation current and receive electromagnetic field component, iteration
Method obtains the APPARENT RESISTIVITY for surveying area stratum.
5. the method according to claim 1, wherein by formula (1) according to whole district's apparent resistance on the stratum
Rate calculates anisotropy coefficient,
λ=ρx(wi)/ρy(wi) formula (1)
Wherein, λ is anisotropy coefficient, ρx(wi) and ρy(wi) it is wiApparent resistivity under frequency in x and y both direction.
6. method according to any one of claims 1 to 5, which is characterized in that source x is equal to source y at a distance from survey district center
At a distance from survey district center, and value range of each field source at a distance from survey district center is 4km~20km.
7. method according to any one of claims 1 to 5, which is characterized in that the stratum includes: land and ocean.
8. a kind of double field source electromagnetic depth methods obtain the anisotropic device of formation resistivity characterized by comprising
First field source and the second field source go out two irrelevant signals for substep alternate emission；Wherein, first field source
The x-axis and y-axis both direction orthogonal in survey area outer with the second field source is laid；
Electric current timing logging modle, for recording the current time sequence Ix (t), Iy of first field source and the transmitting of the second field source
(t)；
Acquisition module is used to survey in area, acquires horizontal component of electric field response signal Ex and Ey orthogonal on each measuring point, and/or, acquisition
Vertical magnetic field response signal Hz on each measuring point, and collected response signal is recorded as time series Ex (t) and Ey (t),
And/or Hz (t)；
Module is obtained, for according to the current time sequence and equipment parameter and reception electromagnetic field component time sequence
Column obtain the APPARENT RESISTIVITY for surveying area stratum by recycling cross-correlation method；
Computing module, for calculating resistivity anisotropy according to the APPARENT RESISTIVITY for surveying area stratum.
9. device according to claim 8, which is characterized in that the acquisition module includes:
Unpolarizable electrode and electric field recorder, for acquiring horizontal component of electric field response signal Ex and Ey orthogonal on each measuring point；With/
Or, magnetic field acquisition module, for acquiring the vertical magnetic field response signal Hz on each measuring point；
Wherein, the magnetic field acquisition module includes magnetograph and broadband induction type magnetic sensor or fluxgate sensing
Device or high-temperature superconducting magnetometer.
10. device according to claim 8 or claim 9, which is characterized in that first field source and the second field source substep alternating two
Secondary excitation: when first field source sends the first signal, second field source sends second signal；First field source sends the
When binary signal, second field source sends the first signal；
Wherein, first signal is single-frequency square-wave signal, and the second signal is inverse repetition M sequence pseudo-random signal, described
The frequency of single-frequency square-wave signal is lower than the inverse low-limit frequency for repeating M sequence pseudo-random signal.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910536766.4A CN110187394B (en) | 2019-06-20 | 2019-06-20 | Method and device for acquiring formation resistivity anisotropy by double-field source electromagnetic sounding method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910536766.4A CN110187394B (en) | 2019-06-20 | 2019-06-20 | Method and device for acquiring formation resistivity anisotropy by double-field source electromagnetic sounding method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110187394A true CN110187394A (en) | 2019-08-30 |
| CN110187394B CN110187394B (en) | 2020-11-27 |
Family
ID=67722692
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910536766.4A Active CN110187394B (en) | 2019-06-20 | 2019-06-20 | Method and device for acquiring formation resistivity anisotropy by double-field source electromagnetic sounding method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110187394B (en) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111983704A (en) * | 2020-09-28 | 2020-11-24 | 西安石油大学 | Method and system for three-dimensional electromagnetic detection between wells |
| CN112696195A (en) * | 2019-10-23 | 2021-04-23 | 中国石油天然气股份有限公司 | Method and device for determining azimuthal anisotropy of formation resistivity |
| CN113917541A (en) * | 2021-09-30 | 2022-01-11 | 湖南科技大学 | Method and device for acquiring electrical anisotropy of underground medium based on direct current electrical prospecting |
| CN116088060A (en) * | 2023-04-12 | 2023-05-09 | 山东大学 | Artificial source electromagnetic exploration system and method based on orthogonal field source-error frequency excitation |
| CN116088059A (en) * | 2023-04-10 | 2023-05-09 | 山东省煤田地质规划勘察研究院 | Artificial source frequency domain electromagnetic exploration method and system based on double-source error frequency emission |
Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2008085063A1 (en) * | 2007-01-11 | 2008-07-17 | Rocksource Geotech As | Method for data acquisition |
| CN101303415A (en) * | 2008-06-24 | 2008-11-12 | 吉林大学 | Combined type field source transmitting apparatus with field source mixing electromagnetical method |
| WO2012106483A2 (en) * | 2011-02-02 | 2012-08-09 | Geco Technology B.V. | Electromagnetic source to produce multiple electromagnetic components |
| CN102721982A (en) * | 2012-06-26 | 2012-10-10 | 中国科学院电工研究所 | Ground electromagnetic prospecting method based on SPSP (Spread Spectrum) coding technology and detection system thereof |
| CN104360403A (en) * | 2014-11-08 | 2015-02-18 | 长安大学 | Multi-component universe apparent resistivity defining method based on multiple-radiation-field-source transient electromagnetic method |
| CN104391332A (en) * | 2014-11-14 | 2015-03-04 | 吉林大学 | Shallow sea double-frequency controllable source electromagnetic prospecting method |
| CN104597506A (en) * | 2015-01-26 | 2015-05-06 | 吉林大学 | Frequency domain ground-to-air electromagnetic prospecting method |
| CN104749640A (en) * | 2015-03-26 | 2015-07-01 | 吉林大学 | Multi-source multi-frequency ground-air electromagnetic detection emission method |
| CN104794313A (en) * | 2015-05-18 | 2015-07-22 | 甘肃省有色地质调查院 | Method and device for obtaining frequency response function of system to be identified |
| CN105607131A (en) * | 2016-01-11 | 2016-05-25 | 甘肃省有色地质调查院 | Method and apparatus for obtaining full-region apparent resistivity spectrum of ground by using coding source electromagnetic sounding method |
| CN105974487A (en) * | 2016-06-12 | 2016-09-28 | 何继善 | Device and method for performing full-area measurement on electric field component Ephi of current source frequency domain |
-
2019
- 2019-06-20 CN CN201910536766.4A patent/CN110187394B/en active Active
Patent Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2008085063A1 (en) * | 2007-01-11 | 2008-07-17 | Rocksource Geotech As | Method for data acquisition |
| CN101303415A (en) * | 2008-06-24 | 2008-11-12 | 吉林大学 | Combined type field source transmitting apparatus with field source mixing electromagnetical method |
| WO2012106483A2 (en) * | 2011-02-02 | 2012-08-09 | Geco Technology B.V. | Electromagnetic source to produce multiple electromagnetic components |
| CN102721982A (en) * | 2012-06-26 | 2012-10-10 | 中国科学院电工研究所 | Ground electromagnetic prospecting method based on SPSP (Spread Spectrum) coding technology and detection system thereof |
| CN104360403A (en) * | 2014-11-08 | 2015-02-18 | 长安大学 | Multi-component universe apparent resistivity defining method based on multiple-radiation-field-source transient electromagnetic method |
| CN104391332A (en) * | 2014-11-14 | 2015-03-04 | 吉林大学 | Shallow sea double-frequency controllable source electromagnetic prospecting method |
| CN104597506A (en) * | 2015-01-26 | 2015-05-06 | 吉林大学 | Frequency domain ground-to-air electromagnetic prospecting method |
| CN104749640A (en) * | 2015-03-26 | 2015-07-01 | 吉林大学 | Multi-source multi-frequency ground-air electromagnetic detection emission method |
| CN104794313A (en) * | 2015-05-18 | 2015-07-22 | 甘肃省有色地质调查院 | Method and device for obtaining frequency response function of system to be identified |
| CN105607131A (en) * | 2016-01-11 | 2016-05-25 | 甘肃省有色地质调查院 | Method and apparatus for obtaining full-region apparent resistivity spectrum of ground by using coding source electromagnetic sounding method |
| CN105974487A (en) * | 2016-06-12 | 2016-09-28 | 何继善 | Device and method for performing full-area measurement on electric field component Ephi of current source frequency domain |
Non-Patent Citations (2)
| Title |
|---|
| 袁博，等: "正交激励双电极源的电场分量表达式和视电阻率提取优势", 《中国地球物理2011》 * |
| 高曙德，等: "编码源地电阻率观测试验", 《地震学报》 * |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112696195A (en) * | 2019-10-23 | 2021-04-23 | 中国石油天然气股份有限公司 | Method and device for determining azimuthal anisotropy of formation resistivity |
| CN112696195B (en) * | 2019-10-23 | 2023-11-28 | 中国石油天然气股份有限公司 | Stratum resistivity azimuth anisotropy determining method and device |
| CN111983704A (en) * | 2020-09-28 | 2020-11-24 | 西安石油大学 | Method and system for three-dimensional electromagnetic detection between wells |
| CN111983704B (en) * | 2020-09-28 | 2023-09-12 | 西安石油大学 | Method and system for three-dimensional electromagnetic detection between wells |
| CN113917541A (en) * | 2021-09-30 | 2022-01-11 | 湖南科技大学 | Method and device for acquiring electrical anisotropy of underground medium based on direct current electrical prospecting |
| CN113917541B (en) * | 2021-09-30 | 2024-02-23 | 湖南科技大学 | Method and device for acquiring electric anisotropy of underground medium based on direct current electric exploration |
| CN116088059A (en) * | 2023-04-10 | 2023-05-09 | 山东省煤田地质规划勘察研究院 | Artificial source frequency domain electromagnetic exploration method and system based on double-source error frequency emission |
| CN116088060A (en) * | 2023-04-12 | 2023-05-09 | 山东大学 | Artificial source electromagnetic exploration system and method based on orthogonal field source-error frequency excitation |
| US12123998B1 (en) | 2023-04-12 | 2024-10-22 | Shandong University | System and method for controlled-source electromagnetic exploration based on staggered frequency excitation of crossed sources |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110187394B (en) | 2020-11-27 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110187394A (en) | Double field source electromagnetic depth methods obtain the anisotropic method and device of formation resistivity | |
| CA2654442C (en) | Method for acquiring and interpreting seismoelectric and electroseismic data | |
| CA2649370C (en) | Integrated earth formation evaluation method using controlled source electromagnetic survey data and seismic data | |
| Anderson et al. | An integrated approach to marine electromagnetic surveying using a towed streamer and source | |
| CA2383931C (en) | Remote reservoir resistivity mapping | |
| RU2358293C2 (en) | System and method of imaging subsurface conductivity | |
| CA2741011C (en) | Method for determining electromagnetic survey sensor orientation | |
| US20060186887A1 (en) | Method for identifying subsurface features from marine transient controlled source electromagnetic surveys | |
| CN104656157B (en) | A kind of method and device in identification shale gas dessert area | |
| WO2007089486A2 (en) | Method for acquiring transient electromagnetic survey data | |
| CN109085653A (en) | A kind of detection method of geology of deep part, sulfide ore body resource | |
| Bechtel et al. | Geophysical methods | |
| MX2011009449A (en) | Electromagnetic survey using naturally occurring electromagnetic fields as a source. | |
| WO2017155737A1 (en) | Magneto-seismic exploration method and system | |
| CN109061750A (en) | A kind of method and device of detection hot dry rock probing target area | |
| CN109407144A (en) | A kind of single hole boulder three-dimensional probe method based on more waves | |
| CN114791633B (en) | Method, system and medium for monitoring shale gas fracturing | |
| Rabbel et al. | Seismoelectric Field Measurements in Unconsolidated Sediments in Comparison With Other Methods of Near‐Surface Prospecting | |
| CN112162323A (en) | Underground unfavorable geological survey method and system based on artificial field source frequency domain electrical method | |
| CN110927632A (en) | Frequency domain horizontal X-direction magnetic field component observation and data processing method | |
| CN117310824B (en) | Thermocouple source excitation controllable source electromagnetic field time sequence forward modeling method and system | |
| CN115032714B (en) | Karst space detection method and system based on comprehensive geophysical prospecting | |
| Rapids et al. | Pilot experiment for the acquisition of marine sediment properties via small scale tomography system | |
| Manataki et al. | Monitoring of water leakages from pipes through geophysical imaging methods | |
| Huang et al. | Limitations of small EM sensors in resistive terrain |
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 |