CN106446408A - Rapid forward and inversion processing method for compensating-while-drilling electromagnetic wave instrument - Google Patents
Rapid forward and inversion processing method for compensating-while-drilling electromagnetic wave instrument Download PDFInfo
- Publication number
- CN106446408A CN106446408A CN201610847357.2A CN201610847357A CN106446408A CN 106446408 A CN106446408 A CN 106446408A CN 201610847357 A CN201610847357 A CN 201610847357A CN 106446408 A CN106446408 A CN 106446408A
- Authority
- CN
- China
- Prior art keywords
- electromagnetic wave
- resistivity
- processing method
- wave instrument
- different
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N5/00—Computing arrangements using knowledge-based models
- G06N5/04—Inference or reasoning models
- G06N5/042—Backward inferencing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N5/00—Computing arrangements using knowledge-based models
- G06N5/04—Inference or reasoning models
- G06N5/046—Forward inferencing; Production systems
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Computing Systems (AREA)
- Data Mining & Analysis (AREA)
- Computational Linguistics (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Artificial Intelligence (AREA)
- Computer Hardware Design (AREA)
- Geometry (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention discloses a rapid forward and inversion processing method for a compensating-while-drilling electromagnetic wave instrument, comprising the steps of S1 acquiring a compensating electromagnetic wave signal; S2 selecting a multi-interface multi-parameter geological initial model according to structural parameters of the compensating-while-drilling electromagnetic wave instrument, well inclination data, and geological information; S4 acquiring corresponding amplitude ratio and phase difference signal according to the principle of measurement, and converting them to obtain different types of resistivity; S5 iterating continuously and comparing with measured results to obtain resistivity parameter values of different layers in a forward model; S6 performing human-machine interaction multi-parameter joint inversion according to stratigraphic model parameter obtained by longitudinal inversion in conjunction with electromagnetic wave data to obtain resistive values of the various layers at different probing depths. The method enables resistivity curves for different transverse strata at different longitudinal depths to be displayed in real time, and geosteering and reservoir assessment are guaranteed.
Description
Technical field
The present invention relates to petroleum drilling log well Exploration Domain and in particular to a kind of with bore compensation electromagnetic wave instrument quickly positive and negative
Drill processing method.
Background technology
At present, carrying out in a deep going way with world's oil-gas exploration and development, emphasis has turned to fine and close oil gas, shale gas, fracture hole
The unconventionaloil pool energy such as type oil-gas reservoir, coal bed gas and deep-sea oil gas, with the evaluating reservoir skill boring geosteering and well logging
Art becomes extremely important, and for boring compensation electromagnetic wave instrument, measurement obtains Amplitude Ratio and the phase place of diverse location receiving coil
After the recovery, needs quickly to be calculated the formation parameter of drill bit and wellbore.So promoting supporting with boring compensation electromagnetic wave instrument
The Fast Solution of Forward And Inverse processing method of data just seems more necessary.
Therefore, present situation in accordance with the above, be highly desirable to propose a kind of for boring the quick of compensation electromagnetic wave instrument
FORWARD AND INVERSE PROBLEMS processing method.
Content of the invention
It is an object of the invention to provide a kind of Fast Solution of Forward And Inverse processing method with brill compensation electromagnetic wave instrument, Neng Goushi
When display transversely Different Strata, the resistivity curve of longitudinal direction upper difference investigation depths, be geosteering and evaluating reservoir offer
Ensure.
In order to achieve the above object, the present invention is achieved through the following technical solutions：A kind of with boring compensation electromagnetic wave instrument
Fast Solution of Forward And Inverse processing method, is characterized in, comprises the steps of：
S1, will with bore compensation electromagnetic wave instrument obtain compensation electromagnetic wave signal be changed, obtain apparent resistivity, according to
Apparent resistivity carries out the mud well plate correction of each Self Matching；
S2, according to boring the structural parameters of compensation electromagnetic wave instrument, hole deviation data and geological information, from log database
Choose the multiple solutions multiparameter geology initial model of coupling, that is, obtain one-dimensional stratiform single shaft isotropic medium inclined shaft just drills mould
Type；
S3, Maxwell equation is applied in the forward model of one-dimensional stratiform single shaft isotropic medium inclined shaft and is imitated
True calculating, obtains many stratigraphic models and produces magnetic field along any direction in each stratum with the coil boring compensation electromagnetic wave instrument
Analytic solutions；
S4, basis, with the measuring principle boring compensation electromagnetic wave instrument, obtain corresponding Amplitude Ratio and phase signal, then by
Different types of resistivity is obtained after the conversion of resistivity chained list；
S5, continuous iteration are simultaneously compared with measured result, until meeting convergence threshold values condition or exceeding the iterationses upper limit,
Obtain the resistance parameter value of different layers position in forward model；
S6, according to the stratigraphic model parameter that longitudinally upper inverting obtains, combine containing length not homology away from, phase contrast and amplitude
The ratio electromagnetic wave data of different qualities, height different frequency, carries out man-machine interaction multi-parameter joint inversion, obtains each layer position different
The resistivity value of investigation depth, thus extrapolating radially invaded zone, virgin zone resistivity, has completed at Fast Solution of Forward And Inverse
Reason.
Step S7 is also comprised after described step S6；
S7, drafting resistivity curve, with display radially invaded zone, virgin zone resistivity directly perceived.
Described comprises 4 emitter coils and 2 receiver coils with brill compensation electromagnetic wave instrument.
Described is 2MHz or 400kHz with the frequency boring the transmission of compensation electromagnetic wave instrument or the electromagnetic wave signal receiving.
Described geological information comprises the drilling well well logging information of geological structure information, seismic profiling information and offset well.
Comprise in described step S5：
When simulation calculation is carried out to the forward model of one-dimensional stratiform single shaft isotropic medium inclined shaft, Magnetic Dipole Source is divided
Solve the vertical magnetic dipole for relative level stratum and horizontal magnetic dipole, add the condition of continuity on layered medium border, knot
Close the Helmholtz equation that electromagnetic field meets, to derive vertical magnetic dipole and horizontal magnetic dipole solution in layers respectively
Analysis solution.
The present invention a kind of with bore compensation electromagnetic wave instrument Fast Solution of Forward And Inverse processing method compared with prior art have with
Lower advantage：It is iterated convergence with reference to quick forward model in longitudinal resistivity inverting to calculate, calculate time-consuming very short, completely
Drilling well can be met calculate in real time, can comprehensively use with bore in compensation electromagnetic wave instrument length not homology away from, phase contrast and width
Degree, than the electromagnetic wave data of different qualities, height different frequency, more fully to verify the |input paramete of model iteration；Using benefit
Repay electromagnetic wave data, the stratigraphic model parameter obtaining with longitudinal inverting carries out man-machine interaction joint inversion, obtain each layer position different
The resistivity value of investigation depth, thus extrapolating the information of invaded zone, undisturbed formation, and is intuitively shown by a plurality of resistivity curve
Illustrate it becomes possible to real-time helping to geosteering and evaluating reservoir in drilling process.
Brief description
Fig. 1 is with the signal transmitting and receiving curve synoptic diagram boring compensation electromagnetic wave instrument；
Fig. 2 is a kind of flow chart with the Fast Solution of Forward And Inverse processing method boring compensation electromagnetic wave instrument of the present invention；
Fig. 3 is conventional resistive rate borehole correction charts schematic diagram；
Fig. 4 is the Magnetic Dipole Source decomposing schematic representation of one-dimensional stratiform single shaft isotropic medium inclined shaft forward model.
Specific embodiment
Below in conjunction with accompanying drawing, by describing a preferably specific embodiment in detail, the present invention is further elaborated.
As shown in figure 1, being that specifically, this instrument comprises 4 with the signal transmitting and receiving curve synoptic diagram boring compensation electromagnetic wave instrument
Individual emitter coil T and 2 receiver coil R, 4 emitter coils are symmetrically arranged amongst the both sides of 2 receiver coils, permissible
Measure the resistivity curve of the phase shift 4 different investigation depths each with amplitude fading of 2MHz with 400kHz height different frequency.
The a kind of of the present invention comprises the steps of with the Fast Solution of Forward And Inverse processing method boring compensation electromagnetic wave instrument, such as Fig. 2
Shown：
S1, will with bore compensation electromagnetic wave instrument obtain compensation electromagnetic wave signal be changed, obtain apparent resistivity, according to
Apparent resistivity carries out the mud well plate correction of each Self Matching.
Because different well sizes and mud resistivity value have different impact to electromagnetic response, therefore for
Weaken or eliminate the effects of the act, need to carry out well plate correction, remove and drawn because of the difference of well size and mud resistivity
The data deviation rising, so that closer to the resistivity on true stratum, as shown in figure 3, well 8in, mud resistivity 0.1ohm-m,
Long (L) short (S) not homology conventional resistive rate borehole correction charts away from, phase contrast (Rps) and Amplitude Ratio (Rad) different qualities.
S2, according to boring the structural parameters of compensation electromagnetic wave instrument, hole deviation data and geological information, from log database
Choose the multiple solutions multiparameter geology initial model of coupling, that is, obtain one-dimensional stratiform single shaft isotropic medium inclined shaft just drills mould
Type.
Geological information comprises the drilling well well logging information of geological structure information, seismic profiling information and offset well.
S3, Maxwell equation is applied in the forward model of one-dimensional stratiform single shaft isotropic medium inclined shaft and is imitated
True calculating, obtains many stratigraphic models and produces magnetic field along any direction in each stratum with the coil boring compensation electromagnetic wave instrument
Analytic solutions.
As shown in figure 4, being the Magnetic Dipole Source exploded view of one-dimensional stratiform single shaft isotropic medium inclined shaft forward model, magnetic
Dipole transmitting-receiving is to calculate the basic model that compensation electromagnetic wave is propagated and received.Because with the line boring compensation electromagnetic wave instrument
Enclose radius with respect to field domain very little to be analyzed, when it passes to the operating current converting in time it is possible to be equivalent to
One electric little ring.In view of the equivalence of electric little ring and magnetic dipole, go approximately to survey with brill using the magnetic dipole being placed on axle center
Transmitting coil in well, receiving coil is the coil of range transmission coil certain intervals, and the induction electromotive force that it produces is according to method
Draw electromagnetic induction principle, the magnetic field in coil array can be obtained, be calculated through the magnetic flux in coil, then to time derivation
Draw induced electromotive force.
Maxwell equation group is expressed as：
AndIn formulaRepresent signature partners
Calculate,For electric field intensity, unit V/m；Represent magnetic field intensity, unit A/m；Represent electric current density, unit A/m2；ω is frequency
Rate, μ represents pcrmeability, thus can obtain：
In formula：
σhRepresent horizontal direction electrical conductivity, σvRepresent vertical direction electrical conductivity, introduce hertz bit function：
Then have：
Lorentz gauge：In formula,Represent Hertzian vector potential, φ represents scalar potential,Represent gradiometer
Calculate,. represent that divergence calculates.
In boundary value problem, hertz position only exists the component in a direction, and direction keeps constant in spatial variations.
Obtain the solution of this component, then given boundary value problem can be obtained according to the relation of hertz position and electric field intensity and magnetic field intensity
Solution.
Vector hertz position and the relation of electric field intensity and magnetic field intensity：
Hertz bit function wave equation is accordingly：
The magnetic dipole of any direction always can be decomposed into vertical perpendicular and parallel two components managing direction in layer, thus appointing
The response problem of meaning direction tilted stratum can be decomposed into vertical dipole (VMD) and two kinds of situations of horizontal dipole (HMD) again
Electromagnetic Calculation problem, respectively Mv, Mh, wherein,
Mv=MTSin θ=TmSTITsinθ
Mh=MTCos θ=TmSTITcosθ
Formula Instrumental axle and stratum relative inclination are θ, and transmitting coil dipole moment is MT, TmFor the transmitting coil number of turn, STFor
Transmitting coil cross-sectional area, ITFor transmitting coil electric current.According to Faraday law of electromagnetic induction, the electromotive force V at receiving coilR
For：
VR=i ω μ SR((Hhx+Hvx)cosθsinα+(Hhz+Hvz)sinθcosα)
In formula, each layer dielectric constant, pcrmeability, electrical conductivity are respectively εi、μi、σi, receiving coil normal direction and instrument press from both sides
Angle is α, and transmitting coil frequency is ω, and receiving coil area is SR, and Hhx, Hvx, Hhz, HvzIt is respectively vertical magnetic dipole, level
The x direction magnetic-field component that magnetic dipole produces and z direction magnetic-field component.
Solve the wave equation that layered medium magnetic dipole can strictly solve hertz magnetic vector potential, obtain in each layer medium
The general solution form of magnetic-dipole field, its solution is the improper integral containing Bessel function.Vector equation is converted into by mark according to model
Amount equation, and then the analytic solutions according to electromagnetic field in boundary condition derivation layered medium.
S4, basis, with the measuring principle boring compensation electromagnetic wave instrument, obtain corresponding Amplitude Ratio and phase signal, then by
Different types of resistivity is obtained after the conversion of resistivity chained list.
S5, continuous iteration are simultaneously compared with measured result, until meeting convergence threshold values condition or exceeding the iterationses upper limit,
Obtain the resistance parameter value of different layers position in forward model.
S6, according to the stratigraphic model parameter that longitudinally upper inverting obtains, combine containing length not homology away from, phase contrast and amplitude
The ratio electromagnetic wave data of different qualities, height different frequency, carries out man-machine interaction multi-parameter joint inversion, obtains each layer position different
The resistivity value of investigation depth, thus extrapolating radially invaded zone, virgin zone resistivity, has completed at Fast Solution of Forward And Inverse
Reason.
With bore in compensation electromagnetic wave data containing length not homology away from, phase contrast and Amplitude Ratio different qualities, just different
The electromagnetic wave data of frequency, the formation characteristics that these data reflect is different, including different investigative ranges, different shape
Feature, the measurement data of different response intensity, can obtain each layer with the stratigraphic model parametric joint inverting that longitudinal FORWARD AND INVERSE PROBLEMS obtains
The resistivity value of the different investigation depth in position, thus extrapolate the information of radially invaded zone, undisturbed formation.
S7, drafting resistivity curve, with display radially invaded zone, virgin zone resistivity directly perceived.
Although present disclosure has been made to be discussed in detail by above preferred embodiment, but it should be appreciated that above-mentioned
Description is not considered as limitation of the present invention.After those skilled in the art have read the above, for the present invention's
Multiple modifications and substitutions all will be apparent from.Therefore, protection scope of the present invention should be limited to the appended claims.
Claims (6)
1. a kind of Fast Solution of Forward And Inverse processing method with brill compensation electromagnetic wave instrument is it is characterised in that comprise the steps of：
S1, will with bore compensation electromagnetic wave instrument obtain compensation electromagnetic wave signal be changed, obtain apparent resistivity, according to regarding electricity
Resistance rate carries out the mud well plate correction of each Self Matching；
S2, basis, with the structural parameters, hole deviation data and the geological information that bore compensation electromagnetic wave instrument, are chosen from log database
The multiple solutions multiparameter geology initial model of coupling, that is, obtain the forward model of one-dimensional stratiform single shaft isotropic medium inclined shaft；
S3, Maxwell equation is applied in the forward model of one-dimensional stratiform single shaft isotropic medium inclined shaft carry out emulation meter
Calculate, obtain the parsing that many stratigraphic models produce magnetic field in each stratum along any direction with the coil boring compensation electromagnetic wave instrument
Solution；
S4, basis, with the measuring principle boring compensation electromagnetic wave instrument, obtain corresponding Amplitude Ratio and phase signal, then by resistance
Different types of resistivity is obtained after the conversion of rate chained list；
S5, continuous iteration are simultaneously compared with measured result, until meeting convergence threshold values condition or exceeding the iterationses upper limit, obtain
The resistance parameter value of different layers position in forward model；
S6, according to the stratigraphic model parameter that longitudinally upper inverting obtains, combine that homology is not away from, phase contrast and Amplitude Ratio containing length
Same characteristic, the electromagnetic wave data of height different frequency, carry out man-machine interaction multi-parameter joint inversion, obtain the different detection in each layer position
The resistivity value of depth, thus extrapolating radially invaded zone, virgin zone resistivity, having completed Fast Solution of Forward And Inverse and having processed.
2. Fast Solution of Forward And Inverse processing method as claimed in claim 1 is it is characterised in that also comprise one after described step S6
Step S7；
S7, drafting resistivity curve, with display radially invaded zone, virgin zone resistivity directly perceived.
3. Fast Solution of Forward And Inverse processing method as claimed in claim 1 it is characterised in that described with boring compensation electromagnetic wave instrument
Comprise 4 emitter coils and 2 receiver coils.
4. Fast Solution of Forward And Inverse processing method as claimed in claim 1 it is characterised in that described with boring compensation electromagnetic wave instrument
The frequency of the electromagnetic wave signal sending or receiving is 2MHz or 400kHz.
5. Fast Solution of Forward And Inverse processing method as claimed in claim 1 is it is characterised in that described geological information inclusively quality structure
Make the drilling well well logging information of information, seismic profiling information and offset well.
6. Fast Solution of Forward And Inverse processing method as claimed in claim 1 is it is characterised in that comprise in described step S5：
When simulation calculation is carried out to the forward model of one-dimensional stratiform single shaft isotropic medium inclined shaft, Magnetic Dipole Source is decomposed into
The vertical magnetic dipole on relative level stratum and horizontal magnetic dipole, add the condition of continuity on layered medium border, in conjunction with electricity
The Helmholtz equation that magnetic field meets, to derive the parsing in layers of vertical magnetic dipole and horizontal magnetic dipole respectively
Solution.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610847357.2A CN106446408A (en) | 2016-09-23 | 2016-09-23 | Rapid forward and inversion processing method for compensating-while-drilling electromagnetic wave instrument |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610847357.2A CN106446408A (en) | 2016-09-23 | 2016-09-23 | Rapid forward and inversion processing method for compensating-while-drilling electromagnetic wave instrument |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN106446408A true CN106446408A (en) | 2017-02-22 |
Family
ID=58167259
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610847357.2A Pending CN106446408A (en) | 2016-09-23 | 2016-09-23 | Rapid forward and inversion processing method for compensating-while-drilling electromagnetic wave instrument |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106446408A (en) |
Cited By (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107621629A (en) * | 2017-09-12 | 2018-01-23 | 重庆梅安森科技股份有限公司 | A kind of underground Precise Position System and underground location method |
| CN108019206A (en) * | 2017-11-16 | 2018-05-11 | 中国石油集团长城钻探工程有限公司 | With brill electromagnetic wave resistivity instrument Range Extension method under a kind of high-k |
| CN108019207A (en) * | 2017-12-15 | 2018-05-11 | 北京六合伟业科技股份有限公司 | A kind of measuring method of SYMMETRIC ELECTROMAGNETIC ripple resistivity |
| CN108547611A (en) * | 2018-03-01 | 2018-09-18 | 杭州迅美科技有限公司 | Horizontal well complex environment is with brill electromagnetic wave resistivity logging rapid simulation method |
| CN109838229A (en) * | 2017-11-27 | 2019-06-04 | 中石化石油工程技术服务有限公司 | A kind of electromagnetic resistivity data processing method |
| CN110005398A (en) * | 2019-04-04 | 2019-07-12 | 中国石油大学(北京) | With brill electromagnetic wave resistivity logging Instrument Design parameter determination method and device |
| CN111075432A (en) * | 2019-12-21 | 2020-04-28 | 中国电波传播研究所（中国电子科技集团公司第二十二研究所） | Mud cake correction method based on electromagnetic wave propagation measuring instrument |
| CN107630697B (en) * | 2017-09-26 | 2020-07-10 | 长江大学 | Formation resistivity joint inversion method based on electromagnetic wave resistivity logging while drilling |
| CN112034257A (en) * | 2020-09-11 | 2020-12-04 | 斯伦贝谢油田技术（山东）有限公司 | Calculation method of downhole resistivity |
| CN113138425A (en) * | 2020-01-19 | 2021-07-20 | 中国石油化工股份有限公司 | Logging-while-drilling electromagnetic wave data acquisition method and device |
| CN113803060A (en) * | 2020-06-01 | 2021-12-17 | 中石化石油工程技术服务有限公司 | Correction method and device for electromagnetic wave resistivity conversion curve scale while drilling |
| CN113803061A (en) * | 2020-06-01 | 2021-12-17 | 中石化石油工程技术服务有限公司 | Coding method and device for orientation-while-drilling electromagnetic wave boundary detection signals |
| CN113868919A (en) * | 2021-09-30 | 2021-12-31 | 西南石油大学 | Simplified method for electromagnetic wave logging while drilling 3D simulation |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20030051914A1 (en) * | 1999-01-28 | 2003-03-20 | Bittar Michael S. | Electromagnetic wave resistivity tool having a tilted antenna for geosteering within a desired payzone |
| US20110074428A1 (en) * | 2009-09-29 | 2011-03-31 | Smith International, Inc. | Apparatus and Method for Downhole Electromagnetic Measurement While Drilling |
| CN102678106A (en) * | 2012-05-02 | 2012-09-19 | 中国电子科技集团公司第二十二研究所 | Data processing method for LWD (Logging While Drilling) electromagnetic wave resistivity logging instrument |
| CN103774988A (en) * | 2013-06-05 | 2014-05-07 | 中国石油大学(华东) | Real-time while-drilling gamma forward modeling geo-steering drilling method |
-
2016
- 2016-09-23 CN CN201610847357.2A patent/CN106446408A/en active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20030051914A1 (en) * | 1999-01-28 | 2003-03-20 | Bittar Michael S. | Electromagnetic wave resistivity tool having a tilted antenna for geosteering within a desired payzone |
| US20110074428A1 (en) * | 2009-09-29 | 2011-03-31 | Smith International, Inc. | Apparatus and Method for Downhole Electromagnetic Measurement While Drilling |
| CN102678106A (en) * | 2012-05-02 | 2012-09-19 | 中国电子科技集团公司第二十二研究所 | Data processing method for LWD (Logging While Drilling) electromagnetic wave resistivity logging instrument |
| CN103774988A (en) * | 2013-06-05 | 2014-05-07 | 中国石油大学(华东) | Real-time while-drilling gamma forward modeling geo-steering drilling method |
Non-Patent Citations (3)
| Title |
|---|
| 杨震，刘庆成，岳步江，马慧斌: "《随钻电磁波电阻率测井仪器响应影响因素数值模拟》", 《测井技术》 * |
| 赵琳琳等: "《补偿型随钻电磁波电阻率测井仪的仿真研究》", 《系统仿真学报》 * |
| 高杰，岳云雷，冯启宁: "《随钻电阻率测井的正演和实用反演方法研究》", 《中国石油学会第十四届测井年会论文集》 * |
Cited By (20)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107621629B (en) * | 2017-09-12 | 2020-07-07 | 重庆梅安森科技股份有限公司 | Underground accurate positioning system and underground positioning method |
| CN107621629A (en) * | 2017-09-12 | 2018-01-23 | 重庆梅安森科技股份有限公司 | A kind of underground Precise Position System and underground location method |
| CN107630697B (en) * | 2017-09-26 | 2020-07-10 | 长江大学 | Formation resistivity joint inversion method based on electromagnetic wave resistivity logging while drilling |
| CN108019206A (en) * | 2017-11-16 | 2018-05-11 | 中国石油集团长城钻探工程有限公司 | With brill electromagnetic wave resistivity instrument Range Extension method under a kind of high-k |
| CN108019206B (en) * | 2017-11-16 | 2018-10-12 | 中国石油集团长城钻探工程有限公司 | With boring electromagnetic wave resistivity instrument Range Extension method under a kind of high-k |
| CN109838229A (en) * | 2017-11-27 | 2019-06-04 | 中石化石油工程技术服务有限公司 | A kind of electromagnetic resistivity data processing method |
| CN108019207A (en) * | 2017-12-15 | 2018-05-11 | 北京六合伟业科技股份有限公司 | A kind of measuring method of SYMMETRIC ELECTROMAGNETIC ripple resistivity |
| CN108547611B (en) * | 2018-03-01 | 2022-02-11 | 杭州迅美科技有限公司 | Method for rapidly simulating logging of electromagnetic wave resistivity while drilling in complex environment of horizontal well |
| CN108547611A (en) * | 2018-03-01 | 2018-09-18 | 杭州迅美科技有限公司 | Horizontal well complex environment is with brill electromagnetic wave resistivity logging rapid simulation method |
| CN110005398A (en) * | 2019-04-04 | 2019-07-12 | 中国石油大学(北京) | With brill electromagnetic wave resistivity logging Instrument Design parameter determination method and device |
| CN110005398B (en) * | 2019-04-04 | 2024-03-22 | 中国石油大学(北京) | Method and device for determining design parameters of electromagnetic wave resistivity logging instrument while drilling |
| CN111075432A (en) * | 2019-12-21 | 2020-04-28 | 中国电波传播研究所（中国电子科技集团公司第二十二研究所） | Mud cake correction method based on electromagnetic wave propagation measuring instrument |
| CN111075432B (en) * | 2019-12-21 | 2022-11-04 | 中国电波传播研究所（中国电子科技集团公司第二十二研究所） | Mud cake correction method based on electromagnetic wave propagation measuring instrument |
| CN113138425A (en) * | 2020-01-19 | 2021-07-20 | 中国石油化工股份有限公司 | Logging-while-drilling electromagnetic wave data acquisition method and device |
| CN113138425B (en) * | 2020-01-19 | 2024-02-27 | 中国石油化工股份有限公司 | Logging-while-drilling electromagnetic wave data acquisition method and device |
| CN113803060A (en) * | 2020-06-01 | 2021-12-17 | 中石化石油工程技术服务有限公司 | Correction method and device for electromagnetic wave resistivity conversion curve scale while drilling |
| CN113803061A (en) * | 2020-06-01 | 2021-12-17 | 中石化石油工程技术服务有限公司 | Coding method and device for orientation-while-drilling electromagnetic wave boundary detection signals |
| CN112034257A (en) * | 2020-09-11 | 2020-12-04 | 斯伦贝谢油田技术（山东）有限公司 | Calculation method of downhole resistivity |
| CN112034257B (en) * | 2020-09-11 | 2023-09-29 | 斯伦贝谢油田技术（山东）有限公司 | Underground resistivity calculation method |
| CN113868919A (en) * | 2021-09-30 | 2021-12-31 | 西南石油大学 | Simplified method for electromagnetic wave logging while drilling 3D simulation |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106446408A (en) | Rapid forward and inversion processing method for compensating-while-drilling electromagnetic wave instrument | |
| US8478530B2 (en) | Using multicomponent induction data to identify drilling induced fractures while drilling | |
| EP2697669B1 (en) | Method for real-time downhole processing and detection of bed boundary for geosteering application | |
| US8364404B2 (en) | System and method for displaying data associated with subsurface reservoirs | |
| US7202670B2 (en) | Method for characterizing a subsurface formation with a logging instrument disposed in a borehole penetrating the formation | |
| US8278931B2 (en) | Anisotropy orientation image from resistivity measurements for geosteering and formation evaluation | |
| US6950749B2 (en) | Method for resistivity anisotropy determination in near vertical wells | |
| US8407005B2 (en) | Method and apparatus for gradient electromagnetic induction well logging | |
| US10295697B2 (en) | Determination of true formation resistivity | |
| US6925384B2 (en) | Method for resistivity anisotropy determination in conductive borehole environments | |
| CN106407574A (en) | A rapid forward modeling and inversion processing method for a multi-component while-drilling azimuth electromagnetic wave instrument | |
| US10317560B2 (en) | Systems and methods of robust determination of boundaries | |
| US8931335B2 (en) | Processing of multi-component induction data in the presence of borehole abnormalities | |
| WO2012037340A1 (en) | Real-time fracture detection and fracture orientation estimation | |
| WO2001067135A1 (en) | Method of interpreting induction logs in horizontal wells | |
| US20180347353A1 (en) | Improved mci logging for processing downhole measurements | |
| Zhang et al. | Directional LWD resistivity tools and their business impacts | |
| US11294092B2 (en) | Low frequency complex resistivity measurement in a formation | |
| CN107939385B (en) | Method for quantitatively calculating polarization value and application | |
| US10295696B2 (en) | Multi-component induction logging data processing in non-circular boreholes | |
| Diren et al. | Numerical analysis of effects of low-resistivity-annulus due to invasion on dual-induction logging responses in deviated hole | |
| GB2417783A (en) | Method for characterising a subsurface formation | |
| Epov et al. | High-frequency EM soundings in horizontal wells: basic theory and applications |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170222 |
|
| RJ01 | Rejection of invention patent application after publication |