CN104090303A - Seismic inversion method and device - Google Patents

Seismic inversion method and device Download PDF

Info

Publication number
CN104090303A
CN104090303A CN201410328284.7A CN201410328284A CN104090303A CN 104090303 A CN104090303 A CN 104090303A CN 201410328284 A CN201410328284 A CN 201410328284A CN 104090303 A CN104090303 A CN 104090303A
Authority
CN
China
Prior art keywords
wave impedance
value
inversion
curve
inverting
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
Application number
CN201410328284.7A
Other languages
Chinese (zh)
Other versions
CN104090303B (en
Inventor
朱卫星
杨玉卿
田召琛
张文博
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Oilfield Services Ltd
China National Offshore Oil Corp CNOOC
Original Assignee
China Oilfield Services Ltd
China National Offshore Oil Corp CNOOC
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China Oilfield Services Ltd, China National Offshore Oil Corp CNOOC filed Critical China Oilfield Services Ltd
Priority to CN201410328284.7A priority Critical patent/CN104090303B/en
Publication of CN104090303A publication Critical patent/CN104090303A/en
Application granted granted Critical
Publication of CN104090303B publication Critical patent/CN104090303B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Abstract

The invention discloses a seismic inversion method and device. The method comprises the steps that according to wave impedance curves of all drilled wells in a research block, it is judged that inversion is feasible, and a wave impedance value of a while-drilling point is obtained in the process of logging while drilling; according to the wave impedance curves and the wave impedance value of the while-drilling point, statistics is performed on wave impedance probability distribution curves of the research block; inversion is performed according to the wave impedance probability distribution curves to obtain an inversion result; a wave impedance value of a next while-drilling point is obtained, and statistics is continued to be performed on the wave impedance probability distribution curves of the research block until the drilling process is completed. According to the seismic inversion method and device, inversion accuracy can be improved. Furthermore, a hole drift angle is adjusted according to the inversion result, and the drill encounter rate of a reservoir of a horizontal well is increased.

Description

A kind of seismic inversion method and apparatus
Technical field
The present invention relates to geosteering field, espespecially a kind of seismic inversion method and apparatus.
Background technology
Seismic inversion refers to and utilizes earth's surface observation seismic data, taking known geologic rule and drilling well, well-log information as constraint, and the process that subterranean strata space structure and physical property are carried out to imaging.
Existing seismic inversion method comprises resistivity inverting spy limit and the density imaging technique etc. that poststack Sparse Pulse seismic inversion, conventional well logging while drilling (LWD, Logging While Drilling) measurement plots changes method of discrimination, connector for logging while drilling carry.
In existing inversion method, for example, poststack Sparse Pulse seismic inversion, because this inversion technique depends on earthquake pure wave data, and the lower earthquake pure wave data of signal to noise ratio (S/N ratio) can cause seismic inversion error larger, Horizon Prognosis is caused to erroneous judgement, and earthquake pure wave data are limited to seismic band, reservoir longitudinal frame is lower, and horizontal physical property that cannot the thin reservoir of Accurate Prediction changes.For example, shale layer only has 2 meters of left and right, and poststack Sparse Pulse seismic inversion is subject to this shale layer of accuracy limitations None-identified.And because well track can not touch shale layer, if this shale layer of inversion method None-identified adopting, well track is difficult to guarantee and does not touch shale layer.
Summary of the invention
In order to address the above problem, the present invention proposes a kind of seismic inversion method and apparatus, can improve inversion accuracy.
In order to achieve the above object, the present invention proposes a kind of seismic inversion method, comprising: judge inverting according to the wave impedance curve of all wells that bored in research block feasible, in well logging while drilling process, obtain with the wave impedance value of boring point;
According to wave impedance curve with the wave impedance probability distribution curve that bores the wave impedance Data-Statistics research block of putting;
According to wave impedance probability distribution curve, inverting obtains inversion result;
Obtain the next wave impedance value with boring point, continue the wave impedance probability distribution curve of statistical research block, until drilling process completes.
It is described that according to wave impedance probability distribution curve, inverting obtains inversion result and comprises:
Described research block is set up to stratigraphic model screen work;
According to described wave impedance curve and the described wave impedance value with boring point, the predictably impedance initial value value of the unknown point in layer model screen work, makes the impedance initial value value in described stratigraphic model screen work meet wave impedance probability distribution curve;
Calculate synthetic seismic record according to the impedance initial value value of the stratigraphic model screen work of prediction, calculate again the residual error between synthetic seismic record and earthquake pure wave data, if the residual error calculating is less than preset value, described inversion result is that the residual error calculating is less than the wave impedance value that preset value is corresponding.
In the time that described residual error is more than or equal to described preset value, the method also comprises:
Revise the wave impedance value in described stratigraphic model screen work, continue to calculate synthetic seismic record according to amended wave impedance value, until the residual error calculating is less than preset value.
Impedance initial value value in described stratigraphic model screen work also meets longitudinal range value and the horizontal range value of the variogram of described research block; Wherein,
Described longitudinal range value obtains according to described wave impedance curve;
Described horizontal range value is according to the horizontal scale information acquisition of reservoir of prediction.
Judge inverting when infeasible when the wave impedance curve of all wells that bored in block according to research, the method also comprises:
Described wave impedance curve is reconstructed, makes inverting feasible.
The method also comprises: adjust hole drift angle according to described inversion result.
Described according to inversion result adjust hole drift angle comprise:
If the reservoir that described inversion result obtains with last inversion result or thickness, inclination angle and the degree of depth of non-reservoir are identical, described hole drift angle is not adjusted; Otherwise calculate new hole drift angle.
The present invention also provides a kind of seismic inversion device, at least comprises:
Acquisition module, feasible for judging inverting according to the wave impedance curve of all wells of having bored of research block, in well logging while drilling process, obtain with the wave impedance value of boring point; Receive circulation indication information, obtain the next one with the wave impedance value of boring point, the next one of acquisition is sent to the wave impedance probability distribution curve of statistical module continuation statistical research block with the wave impedance value of boring point, until drilling process completes;
Statistical module, for according to wave impedance curve with the wave impedance probability distribution curve of wave impedance Data-Statistics research block that bores point;
Inverting module, for inverting obtains inversion result according to wave impedance probability distribution curve, sends circulation indication information to acquisition module.
Described inverting module, specifically for: described research block is set up to stratigraphic model screen work; According to described wave impedance curve and the described impedance initial value value with the unknown point in the wave impedance value prediction stratigraphic model screen work of brill point, make the impedance initial value value in described stratigraphic model screen work meet wave impedance probability distribution curve;
Calculate synthetic seismic record according to the impedance initial value value of the stratigraphic model screen work of prediction, calculate again the residual error between synthetic seismic record and earthquake pure wave data, judge the residual error calculating and be less than preset value, described inversion result is that the residual error calculating is less than the wave impedance value that preset value is corresponding.
Described inverting module, also for: judge described residual error and be more than or equal to described preset value, revise the wave impedance value in described stratigraphic model screen work, continue to calculate synthetic seismic record according to amended wave impedance value, until the residual error calculating is less than preset value.
Described acquisition module, also for: judge that to judge inverting according to the wave impedance curve of all wells of having bored of research block infeasible, described wave impedance curve be reconstructed, make inverting feasible.
Also comprise: adjusting module, for adjusting hole drift angle according to described inversion result, sends described circulation indication information to described acquisition module.
Described adjusting module, specifically for: thickness, inclination angle and the degree of depth of judging reservoir that described inversion result obtains with last inversion result or non-reservoir are identical, described hole drift angle are not adjusted; Judge the thickness of reservoir that described inversion result obtains with last inversion result or non-reservoir or inclination angle or the degree of depth not identical, calculate new hole drift angle.
Compared with prior art, the present invention includes: judge inverting according to the wave impedance curve of all wells that bored in research block feasible, in well logging while drilling process, obtain with the wave impedance value of boring point; According to wave impedance curve with the wave impedance probability distribution curve that bores the wave impedance Data-Statistics research block of putting; According to wave impedance probability distribution curve, inverting obtains inversion result; Obtain the next wave impedance value with boring point, continue the wave impedance probability distribution curve of statistical research block, until drilling process completes.By the solution of the present invention, utilize in well logging while drilling process, obtain with the wave impedance value of boring point, and the wave impedance curve of all wells of having bored in block of research carries out inverting to the wave impedance value of research block, can improve inversion accuracy.
Further, adjust hole drift angle according to inversion result, the reservoir that has improved horizontal well bores encountering rate.
Brief description of the drawings
Below the accompanying drawing in the embodiment of the present invention is described, the accompanying drawing in embodiment is for a further understanding of the present invention, is used from explanation the present invention with instructions one, does not form limiting the scope of the invention.
Fig. 1 is the process flow diagram of seismic inversion method of the present invention;
Fig. 2 is the structure composition schematic diagram of seismic inversion device of the present invention.
Embodiment
For the ease of those skilled in the art's understanding, below in conjunction with accompanying drawing, the invention will be further described, can not be used for limiting the scope of the invention.
Referring to Fig. 1, the present invention proposes a kind of seismic inversion method, comprising:
It is feasible that the wave impedance curve of step 100, all wells of having bored in block according to research is judged inverting, obtains with the wave impedance value of boring point in well logging while drilling process.
In this step, wave impedance curve is the product of densimetric curve and sound wave curve, and densimetric curve and sound wave curve are obtained by logging instrumentation collection.The concrete acquisition of wave impedance curve belongs to those skilled in the art's conventional techniques means, repeats no more here.
In this step, in the time that wave impedance curve can be distinguished reservoir and non-reservoir, think that research block can carry out seismic inversion.For example, the wave impedance value of reservoir is generally 1000~3500, and the wave impedance value of non-reservoir is generally 6000 left and right.
In this step, in the time judging inverting infeasible (being that wave impedance curve cannot be distinguished reservoir and non-reservoir) according to the wave impedance curve of all wells that bored in research block, wave impedance curve is reconstructed, makes inverting feasible (being that wave impedance curve can be distinguished reservoir and non-reservoir).Wherein, can adopt existing algorithm to be reconstructed wave impedance curve, the protection domain being not intended to limit the present invention.
Step 101, the wave impedance Data-Statistics of putting according to wave impedance curve with brill are studied the wave impedance probability distribution curve of block.
In this step, can adopt the wave impedance probability distribution curve of existing method statistic research block, the protection domain being not intended to limit the present invention.
Step 102, according to wave impedance probability distribution curve, inverting obtains inversion result.
In this step, research block is set up to stratigraphic model screen work; According in wave impedance curve and well logging while drilling process, obtain with the wave impedance value of boring point, the predictably impedance initial value value of the unknown point in layer model screen work, makes that primary wave resistance in stratigraphic model screen work is anti-meets wave impedance probability distribution curve; Calculate synthetic seismic record according to the impedance initial value value of the stratigraphic model screen work of prediction, calculate the residual error between synthetic seismic record and earthquake pure wave data, if the residual error calculating is less than preset value, inversion result is that the residual error calculating is less than the wave impedance value that preset value is corresponding.In the time that the residual error calculating is more than or equal to preset value, the wave impedance value in amendment stratigraphic model screen work, continues to calculate synthetic seismic record according to amended wave impedance value, until the residual error calculating is less than preset value.
Wherein, can adopt the wave impedance of existing method (as bayesian criterion) prediction unknown point, the protection domain being not intended to limit the present invention.
Wherein, the impedance initial value value in stratigraphic model screen work also meets longitudinal range value and the horizontal range value of the variogram of research block; Longitudinally range value obtains according to wave impedance curve; Laterally range value is according to the horizontal scale information acquisition of reservoir of prediction.
Wherein, obtain longitudinal range value of the variogram of reservoir and non-reservoir according to formula (1).
v ( h ) = 1 2 N Σ i = 1 N [ z ( x i ) - z ( x i + h ) ] 2 - - - ( 1 )
Wherein, v (h) is variogram value, and h is hysteresis distance, the number that N is location point, x ifor the height of i location point of wave impedance curve, z (x i) be the wave impedance value of i location point of wave impedance curve.
Wherein, when h is increased to infinitely-great process from 0, v (h) tends towards stability gradually, probably equals 1, and when v (h) just starts to stablize, corresponding h is range value.
Wherein, first adopt the wave impedance of existing inversion method (as taking convolution model as basic seismic inversion etc.) inverting research block to distribute, the transverse width that the wave impedance distributed acquisition that obtains according to inverting is each layer, this transverse width is the horizontal range value of the variogram of each layer.
Wherein, adopt the earthquake pure wave data of research block to extract wavelet, then convert the impedance initial value value of stratigraphic model screen work to stratum reflectance value, finally wavelet and stratum reflectance value convolution are obtained to synthetic seismic record.
Wherein, can adopt the wave impedance in existing method (as monte carlo method) amendment stratigraphic model screen work, the protection domain being not intended to limit the present invention.
In this step, can also adjust hole drift angle according to inversion result.
Wherein, if thickness, inclination angle and the degree of depth of the reservoir that inversion result obtains with last inversion result or non-reservoir are identical, hole drift angle is not adjusted; Otherwise calculate new hole drift angle, carry out drilling well according to new hole drift angle, specifically can adopt existing method (as minimum-curvature method) to calculate new hole drift angle, the protection domain being not intended to limit the present invention.
Step 103, obtain the next wave impedance value with boring point, continue the wave impedance probability distribution curve of statistical research block, until drilling process completes.
Referring to Fig. 2, the invention allows for a kind of seismic inversion device, at least comprise:
Acquisition module, feasible for judging inverting according to the wave impedance curve of all wells of having bored of research block, in well logging while drilling process, obtain with the wave impedance value of boring point; Receive circulation indication information, obtain the next one with the wave impedance value of boring point, the next one of acquisition is sent to the wave impedance probability distribution curve of statistical module continuation statistical research block with the wave impedance value of boring point, until drilling process completes;
Statistical module, for according to wave impedance curve with the wave impedance probability distribution curve of wave impedance Data-Statistics research block that bores point;
Inverting module, for inverting obtains inversion result according to wave impedance probability distribution curve, sends circulation indication information to acquisition module.
In seismic inversion device of the present invention, inverting module, specifically for:
Research block is set up to stratigraphic model screen work; According to the impedance initial value value with the unknown point in the wave impedance value prediction stratigraphic model screen work of brill point obtaining in wave impedance curve and well logging while drilling process, make the impedance initial value value in stratigraphic model screen work meet wave impedance probability distribution curve;
Calculate synthetic seismic record according to the impedance initial value value of the stratigraphic model screen work of prediction, calculate again the residual error between synthetic seismic record and earthquake pure wave data, judge the residual error calculating and be less than preset value, inversion result is that the residual error calculating is less than the wave impedance value that preset value is corresponding.
In seismic inversion device of the present invention, inverting module, also for:
Judge residual error and be more than or equal to preset value, the wave impedance value in amendment stratigraphic model screen work, continues to calculate synthetic seismic record according to amended wave impedance value, until the residual error calculating is less than preset value.
In seismic inversion device of the present invention, acquisition module, also for:
It is infeasible that the wave impedance curve of judging all wells that bored in block according to research is judged inverting, and wave impedance curve is reconstructed, and makes inverting feasible.
In seismic inversion device of the present invention, also comprise: adjusting module, for adjusting hole drift angle according to inversion result, sends circulation indication information to acquisition module.
In seismic inversion device of the present invention, adjusting module, specifically for:
Thickness, inclination angle and the degree of depth of judging reservoir that inversion result obtains with last inversion result or non-reservoir are identical, hole drift angle are not adjusted; Judge the thickness of reservoir that inversion result obtains with last inversion result or non-reservoir or inclination angle or the degree of depth not identical, calculate new hole drift angle.
It should be noted that; above-described embodiment understands for the ease of those skilled in the art; be not limited to protection scope of the present invention; do not departing under the prerequisite of inventive concept of the present invention, any apparent replacement and improvement etc. that those skilled in the art make the present invention are all within protection scope of the present invention.

Claims (13)

1. a seismic inversion method, is characterized in that, comprising:
Judge inverting according to the wave impedance curve of all wells that bored in research block feasible, in well logging while drilling process, obtain with the wave impedance value of boring point;
According to wave impedance curve with the wave impedance probability distribution curve that bores the wave impedance Data-Statistics research block of putting;
According to wave impedance probability distribution curve, inverting obtains inversion result;
Obtain the next wave impedance value with boring point, continue the wave impedance probability distribution curve of statistical research block, until drilling process completes.
2. seismic inversion method according to claim 1, is characterized in that, described according to wave impedance probability distribution curve, inverting obtains inversion result and comprises:
Described research block is set up to stratigraphic model screen work;
According to described wave impedance curve and the described wave impedance value with boring point, the predictably impedance initial value value of the unknown point in layer model screen work, makes the impedance initial value value in described stratigraphic model screen work meet wave impedance probability distribution curve;
Calculate synthetic seismic record according to the impedance initial value value of the stratigraphic model screen work of prediction, calculate again the residual error between synthetic seismic record and earthquake pure wave data, if the residual error calculating is less than preset value, described inversion result is that the residual error calculating is less than the wave impedance value that preset value is corresponding.
3. seismic inversion method according to claim 2, is characterized in that, in the time that described residual error is more than or equal to described preset value, the method also comprises:
Revise the wave impedance value in described stratigraphic model screen work, continue to calculate synthetic seismic record according to amended wave impedance value, until the residual error calculating is less than preset value.
4. seismic inversion method according to claim 2, is characterized in that, the impedance initial value value in described stratigraphic model screen work also meets longitudinal range value and the horizontal range value of the variogram of described research block; Wherein,
Described longitudinal range value obtains according to described wave impedance curve;
Described horizontal range value is according to the horizontal scale information acquisition of reservoir of prediction.
5. seismic inversion method according to claim 1, is characterized in that, judges inverting when infeasible when the wave impedance curve of all wells that bored in block according to research, and the method also comprises:
Described wave impedance curve is reconstructed, makes inverting feasible.
6. according to the seismic inversion method described in claim 1~5 any one, it is characterized in that, the method also comprises: adjust hole drift angle according to described inversion result.
7. seismic inversion method according to claim 6, is characterized in that, described according to inversion result adjust hole drift angle comprise:
If the reservoir that described inversion result obtains with last inversion result or thickness, inclination angle and the degree of depth of non-reservoir are identical, described hole drift angle is not adjusted; Otherwise calculate new hole drift angle.
8. a seismic inversion device, is characterized in that, at least comprises:
Acquisition module, feasible for judging inverting according to the wave impedance curve of all wells of having bored of research block, in well logging while drilling process, obtain with the wave impedance value of boring point; Receive circulation indication information, obtain the next one with the wave impedance value of boring point, the next one of acquisition is sent to the wave impedance probability distribution curve of statistical module continuation statistical research block with the wave impedance value of boring point, until drilling process completes;
Statistical module, for according to wave impedance curve with the wave impedance probability distribution curve of wave impedance Data-Statistics research block that bores point;
Inverting module, for inverting obtains inversion result according to wave impedance probability distribution curve, sends circulation indication information to acquisition module.
9. seismic inversion device according to claim 8, is characterized in that, described inverting module, specifically for:
Described research block is set up to stratigraphic model screen work; According to described wave impedance curve and the described impedance initial value value with the unknown point in the wave impedance value prediction stratigraphic model screen work of brill point, make the impedance initial value value in described stratigraphic model screen work meet wave impedance probability distribution curve;
Calculate synthetic seismic record according to the impedance initial value value of the stratigraphic model screen work of prediction, calculate again the residual error between synthetic seismic record and earthquake pure wave data, judge the residual error calculating and be less than preset value, described inversion result is that the residual error calculating is less than the wave impedance value that preset value is corresponding.
10. seismic inversion device according to claim 9, is characterized in that, described inverting module, also for:
Judge described residual error and be more than or equal to described preset value, revise the wave impedance value in described stratigraphic model screen work, continue to calculate synthetic seismic record according to amended wave impedance value, until the residual error calculating is less than preset value.
11. seismic inversion devices according to claim 8, is characterized in that, described acquisition module, also for:
It is infeasible that the wave impedance curve of judging all wells that bored in block according to research is judged inverting, and described wave impedance curve is reconstructed, and makes inverting feasible.
Seismic inversion device described in 12. according to Claim 8~11 any one, is characterized in that, also comprises:
Adjusting module, for adjusting hole drift angle according to described inversion result, sends described circulation indication information to described acquisition module.
13. seismic inversion devices according to claim 12, is characterized in that, described adjusting module, specifically for:
Thickness, inclination angle and the degree of depth of judging reservoir that described inversion result obtains with last inversion result or non-reservoir are identical, described hole drift angle are not adjusted; Judge the thickness of reservoir that described inversion result obtains with last inversion result or non-reservoir or inclination angle or the degree of depth not identical, calculate new hole drift angle.
CN201410328284.7A 2014-07-10 2014-07-10 A kind of seismic inversion method and device Active CN104090303B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410328284.7A CN104090303B (en) 2014-07-10 2014-07-10 A kind of seismic inversion method and device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410328284.7A CN104090303B (en) 2014-07-10 2014-07-10 A kind of seismic inversion method and device

Publications (2)

Publication Number Publication Date
CN104090303A true CN104090303A (en) 2014-10-08
CN104090303B CN104090303B (en) 2016-11-02

Family

ID=51638034

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410328284.7A Active CN104090303B (en) 2014-07-10 2014-07-10 A kind of seismic inversion method and device

Country Status (1)

Country Link
CN (1) CN104090303B (en)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104612660A (en) * 2015-01-23 2015-05-13 中南大学 MRF-based oil and gas reservoir lithofacies stochastic simulation and achievement method
CN104847264A (en) * 2015-05-05 2015-08-19 中国海洋石油总公司 Method and device for achieving geological steering
CN105629310A (en) * 2014-11-05 2016-06-01 中国石油天然气股份有限公司 Well-constraint-free geostatistics inversion method and device for carbonate reservoir
CN110967757A (en) * 2018-09-30 2020-04-07 中国石油化工股份有限公司 Real-time pre-stack seismic inversion method and system while drilling
CN111123359A (en) * 2019-12-24 2020-05-08 同济大学 Logging while drilling and stratum grid constrained well periphery seismic imaging detection method and device
CN112462428A (en) * 2020-11-13 2021-03-09 中国地质科学院 Multi-component seismic data migration imaging method and system

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8116979B2 (en) * 2005-03-09 2012-02-14 Baker Hughes Incorporated System and method for determining a more accurate resistivity model of a geological formation using time-lapse well logging data
WO2013116099A1 (en) * 2012-01-30 2013-08-08 Schlumberger Canada Limited Improving efficiency of pixel-based inversion algorithms
CN103454685A (en) * 2013-08-09 2013-12-18 中国石油天然气股份有限公司 Method and device for predicating sand body thicknesses through logging constraint wave impedance inversion
CN103487835A (en) * 2012-06-12 2014-01-01 中国石油化工股份有限公司 Multi-resolution wave impedance inversion method based on model constraints

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8116979B2 (en) * 2005-03-09 2012-02-14 Baker Hughes Incorporated System and method for determining a more accurate resistivity model of a geological formation using time-lapse well logging data
WO2013116099A1 (en) * 2012-01-30 2013-08-08 Schlumberger Canada Limited Improving efficiency of pixel-based inversion algorithms
CN103487835A (en) * 2012-06-12 2014-01-01 中国石油化工股份有限公司 Multi-resolution wave impedance inversion method based on model constraints
CN103454685A (en) * 2013-08-09 2013-12-18 中国石油天然气股份有限公司 Method and device for predicating sand body thicknesses through logging constraint wave impedance inversion

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
余峰等: "综合地质导向技术及其在海上油田的应用", 《测井技术》 *
徐明华: "苏里格气田苏5区块储层地震反演预测技术", 《天然气工业》 *
朱卫星等: "基于模型的地震波阻抗反演在水平井随钻跟踪中的应用", 《海洋石油》 *

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105629310A (en) * 2014-11-05 2016-06-01 中国石油天然气股份有限公司 Well-constraint-free geostatistics inversion method and device for carbonate reservoir
CN105629310B (en) * 2014-11-05 2018-07-13 中国石油天然气股份有限公司 The method and apparatus of the no restriction from borehole data geostatistical inversion of carbonate reservoir
CN104612660A (en) * 2015-01-23 2015-05-13 中南大学 MRF-based oil and gas reservoir lithofacies stochastic simulation and achievement method
CN104612660B (en) * 2015-01-23 2017-04-26 中南大学 MRF-based oil and gas reservoir lithofacies stochastic simulation and achievement method
CN104847264A (en) * 2015-05-05 2015-08-19 中国海洋石油总公司 Method and device for achieving geological steering
CN104847264B (en) * 2015-05-05 2017-03-08 中国海洋石油总公司 A kind of method and apparatus realizing geosteering
CN110967757A (en) * 2018-09-30 2020-04-07 中国石油化工股份有限公司 Real-time pre-stack seismic inversion method and system while drilling
CN111123359A (en) * 2019-12-24 2020-05-08 同济大学 Logging while drilling and stratum grid constrained well periphery seismic imaging detection method and device
CN112462428A (en) * 2020-11-13 2021-03-09 中国地质科学院 Multi-component seismic data migration imaging method and system
CN112462428B (en) * 2020-11-13 2024-02-20 中国地质科学院 Multi-component seismic data migration imaging method and system

Also Published As

Publication number Publication date
CN104090303B (en) 2016-11-02

Similar Documents

Publication Publication Date Title
Grana Joint facies and reservoir properties inversion
Jacobi et al. Integrated petrophysical evaluation of shale gas reservoirs
Jiménez‐Martínez et al. Temporal and spatial scaling of hydraulic response to recharge in fractured aquifers: Insights from a frequency domain analysis
US11802985B2 (en) Method and system for analyzing filling for karst reservoir based on spectrum decomposition and machine learning
Dutta Deepwater geohazard prediction using prestack inversion of large offset P-wave data and rock model
US10445669B2 (en) System and method for mapping reservoir properties away from the wellbore
US10386531B2 (en) Geological model analysis incorporating cross-well electromagnetic measurements
CN101414013B (en) Method for determining underground fluid by seismic data
CN104090303A (en) Seismic inversion method and device
US11092714B2 (en) Fluid substitution method for T2 distributions of reservoir rocks
WO2007019497A2 (en) Method and system for pre-drill pore pressure prediction
CN105277982B (en) A kind of mud shale total content of organic carbon earthquake prediction method
CN104213899B (en) A kind of Logging Identification Method of formation rock skeleton
CN103592680B (en) A kind of log data based on FORWARD AND INVERSE PROBLEMS and Depth Domain seismic profile synthetic method
CN106154322A (en) Log curve correction method and apparatus
Soleymani et al. Velocity based pore pressure prediction—A case study at one of the Iranian southwest oil fields
CN107861162B (en) Microelectrode logging data-based natural crack identification method and system
Sullivan et al. Permeability from production logs—method and application
CN111077578B (en) Rock stratum distribution prediction method and device
Chou et al. Prediction of hydraulically transmissive fractures using geological and geophysical attributes: a case history from the mid Jhuoshuei River basin, Taiwan
Rotimi et al. Reservoir characterization and modeling of lateral heterogeneity using multivariate analysis
Haustveit* et al. Far-Field Proppant Imaging Offsetting Depletion: A STACK Case History
Johansen et al. Use of Advanced Wireline Logs To Reduce Uncertainties in a Complex Reservoir: A Case Study From the Ivar Aasen Oilfield in the Norwegian Central North Sea
Al Shueili et al. Microfracs and Frac Diagnostic in Tight Gas Reservoirs; Continuous Surveillance to Update Frac Design
Torres-Verdin et al. A Comparison between Geostatistical Inversion and Conventional Geostatistical-simulation Practices for Reservoir Delineation

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CP01 Change in the name or title of a patent holder
CP01 Change in the name or title of a patent holder

Address after: 100010 Chaoyangmen North Street, Dongcheng District, Dongcheng District, Beijing

Co-patentee after: China Oilfield Services Limited

Patentee after: China Offshore Oil Group Co., Ltd.

Address before: 100010 Chaoyangmen North Street, Dongcheng District, Dongcheng District, Beijing

Co-patentee before: China Oilfield Services Limited

Patentee before: China National Offshore Oil Corporation