CN106407574A - A rapid forward modeling and inversion processing method for a multi-component while-drilling azimuth electromagnetic wave instrument - Google Patents
A rapid forward modeling and inversion processing method for a multi-component while-drilling azimuth electromagnetic wave instrument Download PDFInfo
- Publication number
- CN106407574A CN106407574A CN201610847347.9A CN201610847347A CN106407574A CN 106407574 A CN106407574 A CN 106407574A CN 201610847347 A CN201610847347 A CN 201610847347A CN 106407574 A CN106407574 A CN 106407574A
- Authority
- CN
- China
- Prior art keywords
- electromagnetic wave
- information
- processing method
- model
- dimensional
- 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
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention provides a rapid forward modeling and inversion processing method for a multi-component while-drilling azimuth electromagnetic wave instrument. The method comprises the steps of: S1, acquiring compensation electromagnetic wave signals; S2, according to the structure parameters, multi-component information, well deflection data and geological information of a while-drilling azimuth electromagnetic wave instrument, selecting a multi-interface and multi-parameter geological initial model; S3, applying the Maxwell equation to a one-dimension stratiform double-axis anisotropic medium inclined shaft forward modeling model for simulating calculation to obtain conductivity tensor values in three directions; S4, according to the formation model parameters obtained through longitudinal inversion, the electromagnetic wave data and difference signals pointing to different azimuths, performing man-machine interaction combined inversion to obtain specific resistivity values of different investigation depths of each horizon and thus to calculate radial invaded zone, undisturbed formation and formation boundary distance information. The method can realize real-time imagining display of inversed formation interfaces and anisotropic characteristics and provide guarantee for geosteering and reservoir assessment.
Description
Technical field
The Exploration Domain the present invention relates to petroleum drilling is logged well is and in particular to a kind of multi -components are fast with brill orientation electromagnetic wave instrument
Fast FORWARD AND INVERSE PROBLEMS 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, not only inherits traditional resistivity logging while drilling instrument for the multi -components compared with forward position with boring orientation electromagnetic wave instrument
The mentality of designing away from, multi-frequency and symmetric coil structure for the device multi-coil and use advantage, increased crossed coil simultaneously or inclines
Skew coil, obtains bearing signal using coil with the wheel measuring of drill collar, by the voltage ratio of diverse location receiving coil or
Phase contrast and Amplitude Ratio, reflect formation azimuth information.Multi -components are with brill orientation electromagnetic wave logging instrument to resistivity anisotropy
The very high information of orientation unsymmetry such as stratum, fractured reservoir, tilted stratum border has the sensitivity of height, and is wanted
Seek real-time display log data in drilling process to realize geosteering and evaluating reservoir function, so promoting many
Component just seems more necessary with the immediate processing method boring the supporting data of orientation electromagnetic wave instrument.
Therefore, present situation in accordance with the above, be highly desirable to propose a kind of for multi -components with boring orientation electromagnetic wave instrument
Fast Solution of Forward And Inverse processing method.
Content of the invention
It is an object of the invention to provide a kind of multi -components are with the Fast Solution of Forward And Inverse processing method boring orientation electromagnetic wave instrument,
It is capable of realtime imaging and shows the bed boundary of inverting and anisotropic character, be geosteering and evaluating reservoir provides safeguard.
In order to achieve the above object, the present invention is achieved through the following technical solutions：A kind of multi -components are with brill orientation electromagnetic wave
The Fast Solution of Forward And Inverse processing method of instrument, is characterized in, comprises the steps of：
S1, will be changed with boring the compensation electromagnetic wave signal that orientation electromagnetic wave instrument obtains, be obtained apparent resistivity, according to
Apparent resistivity carries out the mud well plate correction of each Self Matching；
S2, according to structural parameters, multi -components information, hole deviation data and the geological information boring orientation electromagnetic wave instrument, from
The multiple solutions multiparameter geology initial model of coupling is chosen, i.e. one-dimensional stratiform biaxial anisotropy's medium inclined shaft in log database
Forward model；
S3, Maxwell equation is applied in the forward model of one-dimensional stratiform biaxial anisotropy's medium inclined shaft and is imitated
True calculating, continuous iteration is simultaneously compared with measured result, until meeting convergence threshold values condition or exceeding the iterationses upper limit
Till, obtain the conductivity tensor value of three different directions；
S4, 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, in conjunction with the difference signal pointing to different azimuth, carries out man-machine interaction connection
Close inverting, obtain the resistivity value of each layer position difference investigation depth, thus extrapolating radially invaded zone, undisturbed formation and ground
Layer frontier distance information, has completed Fast Solution of Forward And Inverse and has processed.
Step S5 is also comprised after described step S4；
S5, combination point to the difference signal of different azimuth, by being imaged bed boundary and the anisotropy spy of real time inversion
Levy and intuitively show.
In described step S3, Maxwell equation is just being applied to one-dimensional stratiform biaxial anisotropy's medium inclined shaft
Drill and also comprise before carrying out simulation calculation in model：Take most basic multilayer parallel sedimentary formation model, by multi -components information
The initial conductivity tensor of setting anisotropic medium different directions.
Adopt separate variables in described step S3, use Fourier transform in horizontal cross-section direction, in Vertical Square
To using One Dimensional Finite Difference Calculation, fast inversion calculates three-dimensional conductivity tensor value.
Described geological information comprises the drilling well well logging information of geological structure information, seismic profiling information and offset well.
The described Maxwell equation containing three-dimensional conductivity tensor is expressed as：
Wherein, three-dimensional conductivity tensor can be embodied as：
In formula,Represent that curl calculates,Represent electric field intensity,Represent magnetic field intensity,Represent dipole moment, ω table
Show frequency, μ represents pcrmeability, ε represents dielectric constant, σx, σy, σzRepresent x, the electrical conductivity on y, z different directions respectively.
A kind of multi -components of the present invention are with the Fast Solution of Forward And Inverse processing method boring orientation electromagnetic wave instrument compared with prior art
There is advantages below：Multi -components are with boring the bearing signal that orientation electromagnetic wave instrument provides to the bearing sense of stratum form and each
Anisotropy feature has natural advantage, and its data processing method is also unique, and there is gun parallax in the geological model storehouse of employing
Different parameter, the complexity of whole model library and comprehensive have very big lifting；For multi -components with brill orientation electromagnetic wave instrument
The three-dimensional magnetic field strength component providing, while comprising more information, simple parsing resolving Algorithm cannot meet complexity
Design conditions, and Finite Element Method (FEM), finite difference calculus (FDM), integral equation method (IE), pattern matching method (NMM) or
Many or limited by calculating time, calculating scale or computation model less, be not suitable for the complicated real-time FORWARD AND INVERSE PROBLEMS in scene
Calculate, the forward model of one-dimensional stratiform biaxial anisotropy's medium inclined shaft that the present invention uses, using separate variables, in water
Plane section direction uses Fourier transform, adopts One Dimensional Finite Difference Calculation in vertical direction, can calculate three by fast inversion
Dimension conductivity tensor value, can reflect the anisotropic character on stratum, the such as fractured reservoir in straight well and inclined shaft well
In tilted stratum；Using with bore in orientation electromagnetic wave data length not homology away from, phase contrast and Amplitude Ratio different qualities, just
The electromagnetic wave data of different frequency, the stratigraphic model parameter obtaining with longitudinal inverting carries out man-machine interaction joint inversion, obtains each
The resistivity value of layer position difference investigation depth, thus extrapolate the information such as invaded zone, undisturbed formation and stratigraphic boundary distance, knot
Close the difference signal pointing to different azimuth, eventually through imaging, the bed boundary of real time inversion and anisotropic character are intuitively shown
Illustrate it becomes possible to real-time in drilling process realize geosteering and evaluating reservoir function.
Brief description
Fig. 1 is a kind of flow chart with the Fast Solution of Forward And Inverse processing method boring orientation electromagnetic wave instrument for the multi -components of the present invention；
Fig. 2 is conventional resistive rate borehole correction charts schematic diagram；
Fig. 3 is simple schematic diagram in the inclined shaft of stratiform anisotropic stratum for the three axle vertical coil；
Fig. 4 is the forward model schematic diagram of one-dimensional stratiform biaxial anisotropy's medium inclined shaft；
Fig. 5 is the disk schematic diagram being close to orientation delineation during bed boundary on the downside of horizontal well Instrumental.
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, a kind of multi -components, with the Fast Solution of Forward And Inverse processing method boring orientation electromagnetic wave instrument, comprise following step
Suddenly：
S1, will be changed with boring the compensation electromagnetic wave signal that orientation electromagnetic wave instrument obtains, be obtained 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 Fig. 2 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.
Fig. 3 is simple schematic diagram in the inclined shaft of stratiform anisotropic stratum for the three axle vertical coil, and sending and receiving coil all contains
Three axle vertical coil will produce multi -components with boring orientation electromagnetic wave signal, specific three-dimensional magnetic field strength component be represented by as
Lower matrix form：
Wherein, H represents magnetic field intensity, and subscript first letter represents transmitting coil, and second letter represents receiving coil.
S2, according to structural parameters, multi -components information, hole deviation data and the geological information boring orientation electromagnetic wave instrument, from
The multiple solutions multiparameter geology initial model of coupling is chosen, i.e. one-dimensional stratiform biaxial anisotropy's medium inclined shaft in log database
Forward model.
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 biaxial anisotropy's medium inclined shaft and is imitated
True calculating, continuous iteration is simultaneously compared with measured result, until meeting convergence threshold values condition or exceeding the iterationses upper limit
Till, obtain the conductivity tensor value of three different directions.
Also comprise to take most basic multilayer parallel sedimentary formation model in step S3, anisotropic medium different directions are set
Initial conductivity tensor.As shown in figure 4, the forward model schematic diagram for one-dimensional stratiform biaxial anisotropy's medium inclined shaft.?
While comprising more information, simple parsing resolving Algorithm cannot meet the design conditions of complexity, and Finite Element Method
(FEM), finite difference calculus (FDM), integral equation method (IE), pattern matching method (NMM) are more or less subject to calculating time, calculating
Scale or the restriction of computation model, the real-time FORWARD AND INVERSE PROBLEMS being not suitable for scene complexity calculates, the one-dimensional layer that the present invention uses
The forward model of shape biaxial anisotropy's medium inclined shaft, using separate variables, is become using Fourier in horizontal cross-section direction
Change, adopt One Dimensional Finite Difference Calculation in vertical direction, three-dimensional conductivity tensor value can be calculated by fast inversion, can be well
Reflect the anisotropic character on stratum, such as the tilted stratum in fractured reservoir and inclined shaft in straight well.
Maxwell (Maxwell) equation containing three-dimensional conductivity tensor is shown in following formula：
Three-dimensional conductivity tensor can be embodied as：
In formulaRepresent that curl calculates,Represent electric field intensity,Represent magnetic field intensity,Represent dipole moment, ω table
Show frequency, μ represents pcrmeability, ε represents dielectric constant, σx, σy, σzRepresent x, the electrical conductivity on y, z different directions respectively.
S4, 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, in conjunction with the difference signal pointing to different azimuth, carries out man-machine interaction connection
Close inverting, obtain the resistivity value of each layer position difference investigation depth, thus extrapolating radially invaded zone, undisturbed formation and ground
Layer frontier distance information, has completed Fast Solution of Forward And Inverse and has processed.
With bore in orientation 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, the stratigraphic model parametric joint inverting obtaining with longitudinal inverting can get each layer position
The resistivity value of different investigation depths, thus extrapolate the information such as invaded zone, undisturbed formation and stratigraphic boundary distance.
S5, combination point to the difference signal of different azimuth, by being imaged bed boundary and the anisotropy spy of real time inversion
Levy and intuitively show it becomes possible to real-time in drilling process realize geosteering and evaluating reservoir function.
As shown in figure 5, being the disk schematic diagram being close to orientation delineation during bed boundary on the downside of horizontal well Instrumental, indication wire
501 position of orientation that the center of circle is located relatively, represent the direction that the stratigraphic boundary closed on is relative to well track, center of circle position is arrived
The distance of indication wire 501, represents instrument track to the distance of stratigraphic boundary, is just close on the downside of marginal data horizontal well track
Bed boundary.
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 multi -components are with boring the Fast Solution of Forward And Inverse processing method of orientation electromagnetic wave instrument it is characterised in that comprising following step
Suddenly：
S1, will be changed with boring the compensation electromagnetic wave signal that orientation electromagnetic wave instrument obtains, be obtained apparent resistivity, according to regarding electricity
Resistance rate carries out the mud well plate correction of each Self Matching；
S2, according to structural parameters, multi -components information, hole deviation data and the geological information boring orientation electromagnetic wave instrument, from well logging
The multiple solutions multiparameter geology initial model of coupling is chosen, that is, one-dimensional stratiform biaxial anisotropy's medium inclined shaft is being just in data base
Drill model；
S3, Maxwell equation is applied in the forward model of one-dimensional stratiform biaxial anisotropy's medium inclined shaft carry out emulation meter
Calculate, continuous iteration is simultaneously compared with measured result, till meeting convergence threshold values condition or exceeding the iterationses upper limit,
Obtain the conductivity tensor value of three different directions；
S4, 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, in conjunction with the difference signal pointing to different azimuth, carry out man-machine interaction joint anti-
Drilling, obtaining the resistivity value of each layer position difference investigation depth, thus extrapolating radially invaded zone, undisturbed formation and stratum side
Boundary's range information, has completed Fast Solution of Forward And Inverse and has 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 S4
Step S5；
S5, combination point to the difference signal of different azimuth, straight by the bed boundary and anisotropic character of imaging real time inversion
Sight shows.
3. Fast Solution of Forward And Inverse processing method as claimed in claim 1 is it is characterised in that by Max in described step S3
Wei Fangcheng is applied to before carrying out simulation calculation in the forward model of one-dimensional stratiform biaxial anisotropy's medium inclined shaft and also comprises：Take
Most basic multilayer parallel sedimentary formation model, arranges the initial conductance of anisotropic medium different directions by multi -components information
Rate tensor.
4. Fast Solution of Forward And Inverse processing method as claimed in claim 1 is it is characterised in that become using separating in described step S3
Amount method, uses Fourier transform in horizontal cross-section direction, adopts One Dimensional Finite Difference Calculation, fast inversion meter in vertical direction
Calculate three-dimensional conductivity tensor value.
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 it is characterised in that described containing three-dimensional conductivity tensor
Maxwell equation is expressed as：
Wherein, three-dimensional conductivity tensor can be embodied as：
In formula,Represent that curl calculates,Represent electric field intensity,Represent magnetic field intensity,Represent dipole moment, ω represents frequency
Rate, μ represents pcrmeability, and ε represents dielectric constant, σx, σy, σzRepresent x, the electrical conductivity on y, z different directions respectively.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610847347.9A CN106407574A (en) | 2016-09-23 | 2016-09-23 | A rapid forward modeling and inversion processing method for a multi-component while-drilling azimuth electromagnetic wave instrument |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610847347.9A CN106407574A (en) | 2016-09-23 | 2016-09-23 | A rapid forward modeling and inversion processing method for a multi-component while-drilling azimuth electromagnetic wave instrument |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN106407574A true CN106407574A (en) | 2017-02-15 |
Family
ID=57997200
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610847347.9A Pending CN106407574A (en) | 2016-09-23 | 2016-09-23 | A rapid forward modeling and inversion processing method for a multi-component while-drilling azimuth electromagnetic wave instrument |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106407574A (en) |
Cited By (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107944159A (en) * | 2017-11-29 | 2018-04-20 | 中国石油集团长城钻探工程有限公司 | One kind is with brill electromagnetic wave forward modeling simulation data base high accuracy compression method |
| CN108019206A (en) * | 2017-11-16 | 2018-05-11 | 中国石油集团长城钻探工程有限公司 | With brill electromagnetic wave resistivity instrument Range Extension method under a kind of high-k |
| CN108073765A (en) * | 2017-12-11 | 2018-05-25 | 中海石油(中国)有限公司湛江分公司 | A kind of horizontal well routine well logging bed boundary identification is with border away from inversion method |
| CN108547611A (en) * | 2018-03-01 | 2018-09-18 | 杭州迅美科技有限公司 | Horizontal well complex environment is with brill electromagnetic wave resistivity logging rapid simulation method |
| CN108643895A (en) * | 2018-04-03 | 2018-10-12 | 中国石油天然气股份有限公司 | Resistivity logging response calculation method and device for gradient stratum |
| CN108875122A (en) * | 2018-04-25 | 2018-11-23 | 杭州迅美科技有限公司 | The artificial intelligence approach and system of geologic parameter are calculated using well logging data |
| CN110552689A (en) * | 2018-05-15 | 2019-12-10 | 中国石油化工股份有限公司 | method for determining distance from while-drilling instrument to stratum boundary |
| CN107630697B (en) * | 2017-09-26 | 2020-07-10 | 长江大学 | Formation resistivity joint inversion method based on electromagnetic wave resistivity logging while drilling |
| CN113570726A (en) * | 2021-08-10 | 2021-10-29 | 中海油田服务股份有限公司 | Multi-buckle while-drilling electrical imaging image generation method and device and computing equipment |
| CN113803061A (en) * | 2020-06-01 | 2021-12-17 | 中石化石油工程技术服务有限公司 | Coding method and device for orientation-while-drilling electromagnetic wave boundary detection signals |
| CN116122794A (en) * | 2023-01-10 | 2023-05-16 | 中国石油大学(华东) | Complex stratum environment array resistivity logging decoupling dimension-reducing forward modeling method |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102678106A (en) * | 2012-05-02 | 2012-09-19 | 中国电子科技集团公司第二十二研究所 | Data processing method for LWD (Logging While Drilling) electromagnetic wave resistivity logging instrument |
| US20130046474A1 (en) * | 2010-07-16 | 2013-02-21 | Haliburton Energy Services, Inc. | Efficient Inversion Systems and Methods for Directionally-Sensitive Resistivity Logging Tools |
| CN104088630A (en) * | 2014-07-09 | 2014-10-08 | 中国石油集团长城钻探工程有限公司 | Method for obtaining stratum boundary distance through azimuth resistivity logging instrument while drilling |
| CN104895561A (en) * | 2015-04-28 | 2015-09-09 | 中国海洋石油总公司 | Electromagnetic propagation resistivity logging joint inversion method based on parallel computing technology |
| CN205207515U (en) * | 2015-12-23 | 2016-05-04 | 上海神开石油化工装备股份有限公司 | Along with boring inclinometer bumper shock absorber |
-
2016
- 2016-09-23 CN CN201610847347.9A patent/CN106407574A/en active Pending
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20130046474A1 (en) * | 2010-07-16 | 2013-02-21 | Haliburton Energy Services, Inc. | Efficient Inversion Systems and Methods for Directionally-Sensitive Resistivity Logging Tools |
| CN102678106A (en) * | 2012-05-02 | 2012-09-19 | 中国电子科技集团公司第二十二研究所 | Data processing method for LWD (Logging While Drilling) electromagnetic wave resistivity logging instrument |
| CN104088630A (en) * | 2014-07-09 | 2014-10-08 | 中国石油集团长城钻探工程有限公司 | Method for obtaining stratum boundary distance through azimuth resistivity logging instrument while drilling |
| CN104895561A (en) * | 2015-04-28 | 2015-09-09 | 中国海洋石油总公司 | Electromagnetic propagation resistivity logging joint inversion method based on parallel computing technology |
| CN205207515U (en) * | 2015-12-23 | 2016-05-04 | 上海神开石油化工装备股份有限公司 | Along with boring inclinometer bumper shock absorber |
Non-Patent Citations (2)
| Title |
|---|
| 杨震等: ""随钻方位电磁波仪器探测电阻率各向异性新方法"", 《石油钻探技术》 * |
| 杨震等: ""随钻电磁波电阻率测井仪器响应影响因素数值模拟"", 《测井技术》 * |
Cited By (16)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| 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 |
| CN107944159A (en) * | 2017-11-29 | 2018-04-20 | 中国石油集团长城钻探工程有限公司 | One kind is with brill electromagnetic wave forward modeling simulation data base high accuracy compression method |
| CN107944159B (en) * | 2017-11-29 | 2018-10-02 | 中国石油集团长城钻探工程有限公司 | One kind is with brill electromagnetic wave forward modeling simulation data base high-precision compression method |
| CN108073765A (en) * | 2017-12-11 | 2018-05-25 | 中海石油(中国)有限公司湛江分公司 | A kind of horizontal well routine well logging bed boundary identification is with border away from inversion method |
| 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 |
| CN108643895A (en) * | 2018-04-03 | 2018-10-12 | 中国石油天然气股份有限公司 | Resistivity logging response calculation method and device for gradient stratum |
| CN108643895B (en) * | 2018-04-03 | 2022-03-01 | 中国石油天然气股份有限公司 | Resistivity logging response calculation method and device for gradient stratum |
| CN108875122A (en) * | 2018-04-25 | 2018-11-23 | 杭州迅美科技有限公司 | The artificial intelligence approach and system of geologic parameter are calculated using well logging data |
| CN110552689A (en) * | 2018-05-15 | 2019-12-10 | 中国石油化工股份有限公司 | method for determining distance from while-drilling instrument to stratum boundary |
| CN113803061A (en) * | 2020-06-01 | 2021-12-17 | 中石化石油工程技术服务有限公司 | Coding method and device for orientation-while-drilling electromagnetic wave boundary detection signals |
| CN113570726A (en) * | 2021-08-10 | 2021-10-29 | 中海油田服务股份有限公司 | Multi-buckle while-drilling electrical imaging image generation method and device and computing equipment |
| CN116122794A (en) * | 2023-01-10 | 2023-05-16 | 中国石油大学(华东) | Complex stratum environment array resistivity logging decoupling dimension-reducing forward modeling method |
| CN116122794B (en) * | 2023-01-10 | 2024-06-11 | 中国石油大学(华东) | Complex stratum environment array resistivity logging decoupling dimension-reducing forward modeling method |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106407574A (en) | A rapid forward modeling and inversion processing method for a multi-component while-drilling azimuth electromagnetic wave instrument | |
| Liu | Principles and applications of well logging | |
| US8478530B2 (en) | Using multicomponent induction data to identify drilling induced fractures while drilling | |
| CN102678106B (en) | Data processing method for LWD (Logging While Drilling) electromagnetic wave resistivity logging instrument | |
| US8407005B2 (en) | Method and apparatus for gradient electromagnetic induction well logging | |
| US8624969B2 (en) | Methods of electromagnetic migration imaging of geologic formation | |
| CN106446408A (en) | Rapid forward and inversion processing method for compensating-while-drilling electromagnetic wave instrument | |
| CN106170605A (en) | Proppant in induced breakage is positioned and the system and method for imaging | |
| US20110254552A1 (en) | Method and apparatus for determining geological structural dip using multiaxial induction measurements | |
| WO2004029664A1 (en) | A method for resistivity anisotropy determination in near vertical wells | |
| CN105074505A (en) | Determination of true formation resistivity | |
| Khoshbakht et al. | Investigating induction log response in the presence of natural fractures | |
| Epov et al. | Electromagnetic tool for high-resolution logging: theoretical and experimental studies | |
| Mikhaylov et al. | Joint inversion of induction and galvanic logging data in axisymmetric geological models | |
| CN107939385B (en) | Method for quantitatively calculating polarization value and application | |
| AU2013356645B2 (en) | Method and system of using a data weighted electromagnetic source for time-lapse monitoring of a reservoir production or hydraulic fracturing | |
| Jahani et al. | Limits of 3D detectability and resolution of LWD deep-sensing borehole electromagnetic measurements acquired in the Norwegian Continental Shelf | |
| Epov et al. | Electromagnetic sounding in deviated and horizontal wells: mathematical modeling and real data interpretation | |
| Nan et al. | Study on the effect of fractures on conventional induction log based on physical and numerical simulation | |
| Lubis et al. | Elastic and electrical properties evaluation of low resistivity pays in Malay Basin clastics reservoirs | |
| Gao et al. | Borehole petrophysical imaging using induction and acoustic measurements | |
| Shetty et al. | New multiphysics, multiscale inversion for imaging petrophysical properties in anisotropic, laminated formations | |
| Bittar et al. | First LWD Co-Located Antenna Sensors for Real-Time Anisotropy and Dip Angle Determination, Yielding Better Look-Ahead Detection | |
| Dutta et al. | Modeling tools for drilling, reservoir navigation, and formation evaluation | |
| CN104612661A (en) | Electromagnetic wave well logging device and method along with drilling |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| WD01 | Invention patent application deemed withdrawn after publication | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20170215 |