CN106443669B - A kind of four antenna receiving array radar logging system imaging methods - Google Patents
A kind of four antenna receiving array radar logging system imaging methods Download PDFInfo
- Publication number
- CN106443669B CN106443669B CN201610596655.9A CN201610596655A CN106443669B CN 106443669 B CN106443669 B CN 106443669B CN 201610596655 A CN201610596655 A CN 201610596655A CN 106443669 B CN106443669 B CN 106443669B
- Authority
- CN
- China
- Prior art keywords
- antenna
- dipole
- angle
- logging system
- antennas
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Geophysics And Detection Of Objects (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The present invention relates to a kind of four antenna receiving array radar logging system imaging methods, beam radar logging system of the invention uses a transmitting antenna Tx and four receiving antenna Rx, pass through waveform minor differences received by four receiving antennas of analysis, so that it is determined that the specific azimuth of underground target geological body, and then three-dimensional imaging by well is carried out using azimuth.The radar logging data processing algorithm for the aerial array that the present invention is suitble to single-shot four to receive, the orientation of objective body can be calculated, and for being imaged, which is appropriate for various related geophysical explorations, can be used in petroleum, coalfield, metallic ore and engineering investigation.
Description
Technical field
The present invention relates to beam radar well logging field, especially a kind of four antenna receiving array radar logging system imaging sides
Method.
Background technique
The transmitting of traditional radar logging and receive all be it is comprehensive, therefore, target geological body depth can only be detected, with
And its distance to wellhole, its orientation can not be detected.Currently, carrying out three-dimensional imaging to underground medium is still
One very difficult thing.
Application No. is 200880119215.8 Chinese patent applications to disclose a kind of be used for stratum progress electromagnetic logging
Method, and the logging equipment being related to.The structure of antenna is to utilize the successful importance of electromagnetic wave logging.The skill of this document
Art scheme uses multiple transmittings and receiving antenna array, is combined to multiple antennas received signal.These arrays include handing over
Pitch dipole antenna and dual crossing dipole antenna.
With the development of technology, further types of antenna is developed.It, also can be with correspondingly, corresponding different antenna
The corresponding imaging algorithm of generation.This invention address that providing a kind of imaging effect more preferably logging system and method.
Summary of the invention
The object of the present invention is to provide the radar logging systems that a kind of single-shot four is received, and devise suitable for the well logging system
The imaging method of system, the logging system can be realized high-precision well logging three-dimensional imaging.
To achieve the above object, the solution of the present invention includes:
A kind of imaging method of four antennas receiving array radar logging system, four antenna receiving array radar logging system packets
A transmitting antenna and four dipole receiving antennas are included, four dipole receiving antennas are parallel, set at equal intervals by circumferencial direction
It sets, each dipole antenna is equal at a distance from two adjacent dipole antennas；Four dipole receiving antennas with
One transmitting antenna is arranged up and down；Steps are as follows for imaging method:
Quadrature coil magnetic-field component is calculated by the electric field component from four dipole antennas；
Target bearing is calculated by orthogonal coil magnetic field.
Further, four dipoles are denoted as East, West, South, North antenna respectively；Electric and magnetic fields are discrete for time sequence
Column, west are to electric fieldEast orientation electric field isThe surrounded region equivalent coil of thing antenna
Surveyed magnetic field isAntenna surrounded the surveyed magnetic field of region equivalent coil in north and south isRoot
According to It calculatesSimilarly, it finds out
Further, with the BewAnd BnsFor x, y；According to Si=xicosαi+yisinαi,It obtains
Its
In, αiIndicate the incident angle of signal, xi、yiRespectively indicate the component in signal level direction and vertical direction, SiIndicate synthesis
Signal, subscript i indicate the sampling number of time series, i1、i2Indicate the first sampled point serial number and most for the time window chosen
The latter sampled point serial number；Two tangent values of the corresponding angle of every bit are acquired by the formula, by seeking arc tangent, often
One sampled point obtains two angle values, the two angle values are brought into formulaTake that angle value so that Q (α i) big, the angle as the point
Angle value.
Beam radar logging system of the invention uses a transmitting antenna Tx and four receiving antenna Rx, passes through analysis four
Waveform minor differences received by a receiving antenna, so that it is determined that the specific azimuth of underground target geological body, and then utilization side
Parallactic angle carries out three-dimensional imaging by well.Radar imagery logger can be placed in well, emitted using transmitting antenna and utilization connects
It receives antenna and receives electromagnetic wave, and judge that the information such as the location fix of objective body, the present invention include being suitble to according to the signal received
In the synthesis of the algorithm and magnetic field data of quadrature receiving coil, the depth of objective body can be judged by well and away from wellbore distance, simultaneously
The orientation of objective body by well can also be calculated, and for being imaged, this method is appropriate for various related geophysical explorations, can
For in petroleum, coalfield, metallic ore and engineering investigation.Common radar logging system can only judge objective body depth and
Away from well spacing from, can not judge target bearing, the radar logging data processing algorithm for the aerial array that the present invention is suitble to single-shot four to receive,
The orientation of objective body can be calculated, and for being imaged, which is appropriate for various related geophysical explorations, can be used for stone
In oil, coalfield, metallic ore and engineering investigation.
Detailed description of the invention
Fig. 1 radar logging directional reception schematic diagram；
Fig. 2 receiving antenna array floor plan, ringlet representative antennas；
Two quadrature coil perspective views of Fig. 3 a, two quadrature coil plan views of Fig. 3 b；
Fig. 4 is with the radar signal schematic diagram of α incidence；
Fig. 5 coil field integral schematic diagram (by taking East and West direction coil as an example).
Specific embodiment
The present invention will be further described in detail with reference to the accompanying drawing.
Beam radar logging system as shown in Figure 1, a transmitting antenna Tx and four dipole receiving antenna Rx, transmitting
Antenna is below receiving antenna.Four electric dipole receiving antennas are distributed on a circumference in the plane perpendicular to wellbore,
And it is angularly distributed, such as Fig. 2.By waveform minor differences received by four receiving antennas of analysis, so that it is determined that underground target
The specific azimuth of geologic body, and then three-dimensional imaging by well is carried out using azimuth.The east, south, west, north provided in Fig. 2 four
Direction is not actual geographic orientation.North orientation need to be only measured in practice, is subject to the antenna for being directed at north orientation, so that it may be corrected it
His aerial position.
It further include the various equipment such as processor about logging system, but these equipment belong to the prior art, therefore no longer
Here it is described in detail.
The imaging method of logging system is specifically described below.
Imaging method mainly includes two aspects:
First, quadrature coil magnetic-field component is calculated by the electric field component from four antennas.Second, by orthogonal coil
Magnetic field calculates target bearing.
It is illustrated in turn below.
(1) it is suitable for the method for quadrature receiving coil
Two are considered first two coil of same position omnidirectional distribution the case where, received physical quantity is magnetic field signal,
See Fig. 3 a, Fig. 3 b.
Since the distance of target geological body in vertical direction can have been provided by ordinary radar well logging, so only discussing
In a plane, the orientation of target geological body.We use magnetic field data, and α is incident at an angle for signal, such as Fig. 4 institute
Show, in this way, incoming signal can be broken down into the component in two vertical direction using incidence point as origin, be respectively labeled as x and
y.In this way, the component in two vertical direction can synthesize incoming signal.If indicating the sampling number of time series with i,
The then available following expression of composite signal.
Si=xicosαi+yisinαi (1)
So the key of this research is to seek α i, we choose the length of a wavelet as time window, one
In a time window, it would be desirable to seek incidence angle using least square method, i.e., in a time window, meet following public affairs
Formula is maximum.
Entire window includes n point, calculates since the last one point of window enters first point together, moves every time
A sampled point is moved, when the last one point of window moves to the last one sampled point of the road, the calculating in the road stops.Entirely
The process of operation is divided into window and does not have completely passed into two parts completely into window, when window not completely into
When, calculate the part between first, road sampled point and window the last one sampled point, when window completely into when, calculate one
Sampled point in a window.
Formula (2) differentiates to α i and can obtain following formula
It can be in the hope of two tangent values of the corresponding angle of every bit, in this way by seeking arc tangent, often by the formula
Available two angle values of one sampled point, the two angle values are brought into formula (2), and we take so that big that of Q (α i)
A angle value, the angle value as the point.
(2) synthesis of magnetic field data
Consider that actual reception system is four dipole antennas, receives electric field signal, need to be converted to electric field signal
Magnetic field signal could use above-mentioned algorithm.
Since in practice, we can only obtain the electric field data in the vertical direction of four receiving antennas, below with
Maxwell equation group synthesizes the magnetic field data in both direction, and analyzes it.
The relationship of variation magnetic field and electric field is expressed as:
Wherein E is electric field intensity, and ds is small surface element vector, and B is magnetic flux density vector, and H is magnetic intensity vector,
Dl is small line element vector.
By taking East and West direction antenna as an example, as shown in figure 5, choose is positive direction clockwise.It enablesWherein
East orientation is that electric field is Ee, west to electric field be Ew, consider the polarization direction of transmitting antenna, the electric field on upper and lower side is Eh, clockwise
Direction is positive direction.Enabling dl is antenna length, lxFor antenna distance.By electric and magnetic fields it is discrete be time series, west is to electric field
ForEast orientation electric field isThe surveyed magnetic field B of the surrounded region equivalent coil of thing antennaewForThe surveyed magnetic field B of the surrounded region equivalent coil of north and south antennansFor
Discrete, a certain moment is carried out to equation (4),
Dt is time step in formula, and ds is area coil.
Two formulas are merged, following relational expression can be obtained,
It considersIt can calculate
Similarly, it can find out,BewAnd Bns, i.e. x and y in Fig. 4.
Specific embodiment of the present invention is presented above, but the present invention is not limited to described embodiment.
Under the thinking that the present invention provides, to the skill in above-described embodiment by the way of being readily apparent that those skilled in the art
Art means are converted, are replaced, are modified, and play the role of with the present invention in relevant art means it is essentially identical, realize
Goal of the invention it is also essentially identical, the technical solution formed in this way is to be finely adjusted to be formed to above-described embodiment, this technology
Scheme is still fallen in protection scope of the present invention.
Claims (1)
1. a kind of four antenna receiving array radar logging system imaging methods, four antenna receiving array radar logging systems include one
A transmitting antenna and four dipole receiving antennas, four dipole receiving antennas are parallel, are arranged at equal intervals by circumferencial direction, often
A dipole antenna is equal at a distance from two adjacent dipole antennas；Four dipole receiving antennas and a hair
Antenna is penetrated to arrange up and down；Steps are as follows for imaging method:
Quadrature coil magnetic-field component is calculated by the electric field component from four dipole antennas；
Target bearing is calculated by orthogonal coil magnetic field；
Four dipoles are denoted as East, West, South, North antenna respectively；By electric and magnetic fields it is discrete be time series, west to electric field isEast orientation electric field isSurrounded the surveyed magnetic field of region equivalent coil of thing antenna isAntenna surrounded the surveyed magnetic field of region equivalent coil in north and south isAccording to It calculatesSimilarly, it finds out
With the BewAnd BnsFor x, y；According to Si=xi cosαi+yi sinαi,It obtainsWherein, αiIt indicates
The incident angle of signal, xi、yiRespectively indicate the component in signal level direction and vertical direction, SiIndicate composite signal, subscript
I indicates the sampling number of time series, i1、i2Indicate that first sampled point serial number of the time window chosen is adopted with the last one
Sample point number；Two tangent values that the corresponding angle of every bit is acquired by the formula, by seeking arc tangent, each sampling
Point obtains two angle values, the two angle values are brought into formulaIt takes so that Q (αi) big that angle value, the angle as the point
Angle value.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610596655.9A CN106443669B (en) | 2016-07-26 | 2016-07-26 | A kind of four antenna receiving array radar logging system imaging methods |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610596655.9A CN106443669B (en) | 2016-07-26 | 2016-07-26 | A kind of four antenna receiving array radar logging system imaging methods |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106443669A CN106443669A (en) | 2017-02-22 |
| CN106443669B true CN106443669B (en) | 2019-02-26 |
Family
ID=58185154
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610596655.9A Active CN106443669B (en) | 2016-07-26 | 2016-07-26 | A kind of four antenna receiving array radar logging system imaging methods |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106443669B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115561743A (en) * | 2020-07-17 | 2023-01-03 | 深圳市安卫普科技有限公司 | Nonlinear node detection method and detector |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1769645A (en) * | 2004-11-03 | 2006-05-10 | 长沙中联重工科技发展股份有限公司 | Guide positioning method and positioning instrument for horizontal directional drilling |
| CN102239633A (en) * | 2008-09-27 | 2011-11-09 | 韦特里西提公司 | Wireless energy transfer systems |
| CN104727812A (en) * | 2013-12-24 | 2015-06-24 | 中国石油化工集团公司 | Measuring set and method of while-drilling azimuthal electromagnetic valve resistivity |
| CN104931924A (en) * | 2015-06-12 | 2015-09-23 | 华南理工大学 | Electromagnetic field angular momentum-based antenna positioning method |
-
2016
- 2016-07-26 CN CN201610596655.9A patent/CN106443669B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1769645A (en) * | 2004-11-03 | 2006-05-10 | 长沙中联重工科技发展股份有限公司 | Guide positioning method and positioning instrument for horizontal directional drilling |
| CN102239633A (en) * | 2008-09-27 | 2011-11-09 | 韦特里西提公司 | Wireless energy transfer systems |
| CN104727812A (en) * | 2013-12-24 | 2015-06-24 | 中国石油化工集团公司 | Measuring set and method of while-drilling azimuthal electromagnetic valve resistivity |
| CN104931924A (en) * | 2015-06-12 | 2015-09-23 | 华南理工大学 | Electromagnetic field angular momentum-based antenna positioning method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106443669A (en) | 2017-02-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Vozoff | The magnetotelluric method | |
| Stolz et al. | Superconducting sensors and methods in geophysical applications | |
| CA2895018C (en) | Deep azimuthal system with multi-pole sensors | |
| Ebihara et al. | Super-resolution of coherent targets by a directional borehole radar | |
| Takayama et al. | A novel direction-finding algorithm for directional borehole radar | |
| Sheinker et al. | Remote tracking of a magnetic receiver using low frequency beacons | |
| Gaballah et al. | Characterizing subsurface archaeological structures with full resolution 3D GPR at the early dynastic foundations of Saqqara Necropolis, Egypt | |
| CN102944905A (en) | Gravity and magnetic anomaly handling method based on directional wavelet analysis | |
| CN103439707A (en) | Method for migration imaging of windowing distance of ground penetrating radar | |
| CN106443669B (en) | A kind of four antenna receiving array radar logging system imaging methods | |
| Gürer et al. | A VLF survey using current gathering phenomena for tracing buried faults of Fethiye–Burdur Fault Zone, Turkey | |
| CN106593420B (en) | Interface detection method by a kind of multi-source combination sound reflecting logging instrument and well | |
| Hager et al. | Charge rearrangement by sprites over a north Texas mesoscale convective system | |
| CN106772677A (en) | A kind of method for finding area of coverage skarn type Eastern Pamir | |
| CN107884821B (en) | A method of estimation microseism direction of arrival and distance | |
| Thiesson et al. | Detection of resistive features using towed Slingram electromagnetic induction instruments | |
| CN104407386A (en) | Electromagnetic exploration data acquiring and processing method and system for executing method | |
| CN105006648B (en) | Double frequency for high frequency marine radar receives and dispatches spike antenna system altogether | |
| Sternberg et al. | Numerical studies and potential applications of the vertical array-differential target antenna coupling (DTAC) method for rapid sensing and imaging of subsurface targets | |
| Lymburner et al. | A procedure for collecting electromagnetic data using multiple transmitters and receivers capable of deep and focused exploration | |
| Chuang et al. | Prediction of antenna induced biases for precision GNSS receivers | |
| Buzin et al. | Search for Chelyabinsk Meteorite Fragments in Chebarkul Lake Bottom (GPR and Magnetic Data) | |
| GB2458959A (en) | Electromagnetic surveying using a crossed electric dipole transmitter with selectable orientations | |
| Ebihara et al. | Singlehole borehole rader measurement using dipole array antenna fed by coaxial cable | |
| Owen et al. | Borehole directional radar detection of subsurface cavities |
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 | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20220208 Address after: 100028 Chaoyang District, Beijing Hui Xin Street 6, Twelfth level. Patentee after: SINOPEC OILFIELD SERVICE Corp. Patentee after: Sinopec Jingwei Co.,Ltd. Patentee after: SINOPEC NORTH CHINA PETROLEUM ENGINEERING Co.,Ltd. Patentee after: North China Measurement and control company of Sinopec Jingwei Co.,Ltd. Address before: 100728 9 Jishi Kou Road, Chaoyang District, Beijing Patentee before: SINOPEC OILFIELD SERVICE Corp. Patentee before: Sinopec North China Petroleum Engineering Co., Ltd |
|
| TR01 | Transfer of patent right |