CN101887132A - Method for quantificationally predicting sandstone reservoir fluid saturation by combining well and seism - Google Patents

Method for quantificationally predicting sandstone reservoir fluid saturation by combining well and seism Download PDF

Info

Publication number
CN101887132A
CN101887132A CN200910084537XA CN200910084537A CN101887132A CN 101887132 A CN101887132 A CN 101887132A CN 200910084537X A CN200910084537X A CN 200910084537XA CN 200910084537 A CN200910084537 A CN 200910084537A CN 101887132 A CN101887132 A CN 101887132A
Authority
CN
China
Prior art keywords
modulus
curve
volume
elastic parameter
saturation
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
Application number
CN200910084537XA
Other languages
Chinese (zh)
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 Petroleum and Natural Gas Co Ltd
Original Assignee
China Petroleum and Natural Gas Co Ltd
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 Petroleum and Natural Gas Co Ltd filed Critical China Petroleum and Natural Gas Co Ltd
Priority to CN200910084537XA priority Critical patent/CN101887132A/en
Publication of CN101887132A publication Critical patent/CN101887132A/en
Pending legal-status Critical Current

Links

Images

Abstract

The invention discloses a method for quantificationally predicting sandstone reservoir fluid saturation by combining a well and seism in geophysical prospecting for petroleum. The method comprises the following steps of: acquiring and processing seismic data and log data to obtain each parameter curve; respectively determining the relationships among the shale content, the porosity and an elastic parameter of the well; establishing prediction models of the shale content and the porosity; calculating a pore volume modulus and a dry rock volume modulus to establish a modulus prediction model; establishing an effective reservoir water saturation prediction model; performing prestack seismic inversion by utilizing prestack seismic trace gathers and the log data to obtain an inversion data volume of each elastic parameter; calculating the dry rock volume modulus to solve the pore volume modulus by the obtained elastic inversion volume; and solving the water saturation by using the pore volume modulus. In the method, prestack seismic data and the log data are fully utilized for quantificationally predicting reservoir saturation; and good geological effects on predicting oil and gas saturation are achieved.

Description

A kind of method of quantificationally predicting sandstone reservoir fluid saturation by combining well and seism
Technical field
The present invention relates to the geophysical prospecting for oil technology, specifically is a kind of method of the quantificationally predicting sandstone reservoir fluid saturation by combining well and seism based on the volume of voids modulus.
Background technology
Utilize seismic data the pore fluid of underground rock is predicted it is the target of seismic work.One of means that AVO (amplitude changes with geophone offset) attributive analysis directly detects as oil gas, once in bright spot and the detection of flat spot type gas reservoir, brought into play vital role, yet the AVO attributive analysis is subjected to tuning influence, can only carry out qualitative description to reservoir fluid, and has a lot of traps.Seismic inversion, the especially appearance of pre-stack seismic inversion technique make that utilizing seismic data that the reservoir pore space fluid is carried out quantitative forecast becomes possibility.In recent years, utilize geological data to carry out petroleum-gas prediction and become a focus, but main achievement in research still concentrates on the qualitative description and the prediction aspect of convection cell, as the AVO intersection that Ruthorford and Castanga etc. propose, the wavelet energy analysis method of propositions such as Peters etc. all can only carry out qualitative forecasting to the gas-bearing property of reservoir.Li Hongbing has carried out the gas-bearing property prediction with the wavelet scale spectrometry to Caidamu Basin gas reservoir; Sweet sharp lamps etc. have been discussed the potentiality of elastic impedance in lithology and fluid prediction in detail; Li Zongjie etc. utilize various seismic properties in system in Tahe Oilfield the oiliness of Ordovician system oil-containing reservoir to be predicted, and have obtained certain effect; Comprehensive utilization infiltration rate frequency dispersions (AVD) such as Li dignitary, dynamic power spectrum technology such as (DR) have been carried out the gas-bearing property detection to the non-homogeneous tight sand of Xu Jia river, gas field, new field gas reservoir; Li Aishan etc. have showed the application of the synchronous prestack inversion method of AVA in Shengli Oil Field Jiyang depression middle level gas reservoir is effectively discerned; Gao Jianhu etc. have proposed the application of poststack oil and gas testing technique in oil and gas detection such as comprehensive utilization effective absorption coefficient, seismic event kinetic parameter and spectral amplitude identification.
Also there are many scholars to explore domestic aspect the fluid saturation quantitative forecast, as, Zhang Eyong multiple regression procedure and more than 20 poststack seismic properties have carried out predicting (petroleum prospecting and exploitation to the sandstone thickness and the zone of interest of study area, 2000,27 rolled up for 1 phases), and obtained effect preferably.But because this method is based on poststack seismic properties rather than seismic inversion parameter, make this method have resolution low, the quality of seismic data is relied on shortcomings such as serious.The thinking that Li Lailin has proposed to utilize earthquake and log data to ask for static remaining oil saturation based on time-average equation (grand celebration oil geology and exploitation,, 21 3 phases of volume in 2000) has the certain experiences meaning to the remaining oil prediction in oil and gas development district.Liu Yang has proposed a kind of method and thinking (petroleum journal that utilizes geological data estimation rock porosity and fluid saturation, 2005,26 2 phases of volume), its method is at first to utilize the Changing Pattern of elastic parameters such as Gassmann Model Calculation p-and s-wave velocity, density, speed compare with factor of porosity, saturation degree, and set up the respective explanations plate, obtain elastic parameter by various inversion methods then, and with the interpretation chart comparative analysis, and then factor of porosity and saturation degree carried out the sxemiquantitative estimation, its method and thinking have certain promotion value.Jin Long etc. have proposed a kind of new method (petroleum journal based on petrophysical model and hybrid optimization algorithm while inverting factor of porosity and saturation degree, 2006,27 4 phases of volume), this method adds refutation process as constraint with factor of porosity and saturation degree, its stability and validity have been passed through the check of The model calculation, but are not applied in the actual seismic data.As seen, the method that does not also have a kind of practicality, effective quantitative forecast saturation degree at present both at home and abroad.
Summary of the invention
The object of the present invention is to provide a kind of method of the quantificationally predicting sandstone reservoir fluid saturation by combining well and seism based on the volume of voids modulus.
The present invention is achieved through the following technical solutions:
1) adopt arrangement greatly, small distance between receivers, high-density sampling is gathered high-quality geological data, seismic data is protected the width of cloth handle, and forms the common midpoint gather of fidelity;
It is to keep the AVO feature in processing procedure that the described guarantor's width of cloth of step 1) is handled.
2) gather compressional wave time difference, density, neutron, spontaneous potential or GR, resistivity conventional logging data, obtain logging trace, and processing obtains shale index, factor of porosity, saturation curves.
Step 2) described collection data is gathered the shear wave slowness curve to the well more than a bite at least.
3) obtain longitudinal and transverse wave velocity curve from compressional wave time difference, shear wave slowness curve, calculate p-wave impedance according to p-and s-wave velocity, densimetric curve, the shear wave impedance, p-and s-wave velocity is than elastic parameter curve;
The described shear wave velocity curve of step 3) utilizes Xu-White Model Calculation shear wave velocity curve when not gathering the well of shear wave slowness curve.
4) utilize multivariate statistical method to determine relation between shale index and factor of porosity and the elastic parameter respectively at well point place, set up the forecast model of following shale index and factor of porosity:
V Sh=f sh(L 1,L 2,L 3,...) (1)
φ=f φ(L 1,L 2,L 3,...) (2)
In the formula: V ShThe expression shale index, φ represents factor of porosity, f Sh, f φAll represent the multivariate function, independent variable is that elastic parameter, dependent variable are respectively shale index and factor of porosity, L 1, L 2, L 3... L nBe the elastic parameter curve;
5) adopt following formula to calculate volume of voids modulus and dry rock bulk modulus:
Figure B200910084537XD0000031
K d=K sat-K p (4)
In the formula: K SatThe bulk modulus of representing complete saturated rock, K dThe bulk modulus of expression dry rock, K pExpression hole bulk modulus, φ represents factor of porosity, K sThe bulk modulus of expression rock matrix, K fThe bulk modulus of expression pore fluid;
6) utilize dry rock bulk modulus and elastic parameter curve, set up the forecast model of dry rock bulk modulus:
K d=f kd(L 1,L 2,L 3,...) (5)
In the formula: L 1, L 2, L 3... L nBe elastic parameter curve, f KdExpression is that independent variable, dry rock bulk modulus are the multivariate function of dependent variable with the elastic parameter;
7) utilize multiple regression procedure to set up the forecast model of water saturation in effective reservoir:
S w=f sw(K p,L 4,L 5,...) (6)
In the formula: L 1, L 2, L 3Expression elastic parameter curve, S wBe water saturation, f SwExpression is that independent variable, water saturation are the multivariate function of dependent variable with the elastic parameter.
8) utilize pre-stack seismic road collection and well-log information to carry out the pre-stack seismic inverting, obtain the inverting data volume of various elastic parameters, and utilize formula (1) and formula (2) to calculate shale index and factor of porosity data volume from the elastic parameter inversion data volume;
9) forecast model that utilizes step 5) to set up calculates the dry rock bulk modulus from step 8) gained elastic parameter inversion data volume;
10) utilize inverting elastic parameter data volume to ask the bulk modulus of saturated rock, and then ask volume of voids modulus K by formula (4) p
11) utilize formula (6) to ask water saturation S from the volume of voids modulus w
The present invention has made full use of pre-stack seismic data and well-log information, has realized the quantitative forecast to the reservoir pore space fluid saturation, and on-the-spot test has good geological effect to oil-containing, gas saturation prediction.
Description of drawings
Fig. 1 is the technology of the present invention process flow diagram;
Fig. 2 for the present invention in shale index, factor of porosity and the saturation curves of somewhere well logging interpretation with by the comparison diagram of the curve of elastic parameter prediction; (solid line is a measured curve, and dotted line is a prediction curve)
Fig. 2 (a) is shale index prediction effect figure;
Fig. 2 (b) is the porosity prediction design sketch;
Fig. 2 (c) is water saturation prediction effect figure;
Fig. 3 (a) protects the width of cloth for certain survey line and handles the CDP road collection that obtains;
Fig. 3 (b) is the stack result of CDP road collection shown in Fig. 3 (a);
Fig. 4 (a) is certain survey line shale index sectional view;
Fig. 4 (b) is for laterally following the trail of effective reservoir of gained according to the shale index section;
Fig. 4 (c) is the factor of porosity section that calculates at effective reservoir;
Fig. 5 is the volume of voids modulus sectional view of certain survey line;
Fig. 6 is certain survey line water saturation sectional view;
Fig. 7 is certain survey line gas saturation and factor of porosity product sectional view.
Embodiment
The present invention utilizes well logging and pre-stack seismic data quantitative forecast sandstone reservoir fluid saturation, at first by the water saturation forecast model of well-log information foundation based on elastic parameter, utilize the pre-stack seismic inverting to obtain required elastic parameter data volume then, calculate water saturation in conjunction with forecast model, realize the quantitative forecast of reservoir fluid.Fig. 1 is the technology of the present invention process flow diagram.
The present invention will be further described below in conjunction with example:
1) adopt arrangement greatly, small distance between receivers, high-density sampling is gathered high-quality geological data, seismic data is protected the width of cloth handle, and forms the common midpoint gather of fidelity, as shown in Figure 3.Fig. 3 (a) and (b) be respectively certain and distinguish certain survey line and protect CDP road collection and corresponding stacked section after width of cloth is handled.Zone of interest is 8 sections in a box, and its corresponding volume lineups in end, top mark (among the patent figure color can not be arranged, interchangeable saying) in Fig. 3 (b).
2) gather this district's compressional wave time difference, density, neutron, spontaneous potential or GR, resistivity conventional logging data, obtain logging trace, and explained shale index, factor of porosity, saturation curves; This area has Liang Koujing that the shear wave slowness curve of actual measurement is arranged.
3) surveying shear wave slowness with two causes for gossip is constraint, the shear wave slowness curve of all the other wells that utilized the Xu-White Model Calculation.According to compressional wave time difference, the longitudinal and transverse wave velocity curve of shear wave slowness curve calculation, and then calculate p-wave impedance, shear wave impedance, p-and s-wave velocity according to p-and s-wave velocity, densimetric curve and the elastic parameter curve such as compare.
4) utilize multivariate statistical method to determine relation between shale index and factor of porosity and the elastic parameter respectively at well point place, set up the forecast model of following shale index and factor of porosity.
Shale index (V Sh) available following nonlinear function calculating:
V Sh=f(V p,V s,ρ),R=0.836 (7)
In the formula: V ShThe expression shale index, V pThe expression velocity of longitudinal wave, V SThe expression shear wave velocity, ρ represents density, R represents related coefficient.
In effective reservoir, factor of porosity can directly carry out linear transformation from density and obtain:
φ=85.642-30.83ρ,R=0.79 (8)
In the formula: φ represents factor of porosity, and ρ represents density, and R represents related coefficient.
Fig. 2 (a) and Fig. 2 (b) are respectively the comparison with result of log interpretation of predicting the outcome of shale index and factor of porosity, and both coincide fine, illustrates predicting the outcome of well point shale index and factor of porosity reliably, and this accurately predicts for water saturation and lays a good foundation.
5) calculate volume of voids modulus (K according to formula (3) p), and then according to formula (4) calculating dry rock bulk modulus (K d).
6) utilize dry rock bulk modulus and elastic parameter curve, set up the forecast model of dry rock bulk modulus:
K d=30.26+0.1721V sh-101.37φ,R=0.80 (9)
In the formula: K dBe dry rock bulk modulus, V ShBe shale index, φ is a hole, and R is a related coefficient;
7) by statistics, the relation of water saturation and volume of voids modulus has difference slightly in local area I class sandstone and the II class sandstone,
In I class sandstone:
S w=20.42K p+18.26,R=0.93 (10a)
In II class sandstone:
S w=14.11K p+32.95,R=0.79 (10b)
According to formula (10a) and (10b) prediction saturation curves and well logging interpretation saturation curves relatively see Fig. 2 (c), the two is the variation tendency unanimity not only, and numerical value is close, has illustrated that the saturation degree forecast model of setting up is effective.
8) utilize pre-stack seismic road collection and well-log information to carry out the pre-stack seismic inverting, obtained the inversion result of velocity of longitudinal wave, shear wave velocity and density, and calculated p-wave impedance, shear wave impedance, p-and s-wave velocity and elastic parameter such as compared.
By the pre-stack seismic inverting gained elastic parameter calculating of the forecast model shown in (7) by formula shale index data volume, the result is shown in Fig. 4 (a).According to the shale index section reservoir is laterally followed the trail of and to be obtained effective reservoir distribution figure, as Fig. 4 (b).At effective reservoir by elastic parameter by formula (8) calculate factor of porosity data volume, result such as Fig. 4 (c).
Fig. 4 (a) has shown the shale index curve among 4 (b), and the curve that shows among Fig. 4 (c) is a porosity curve.As seen from the figure, shale index and factor of porosity predict the outcome and logging trace coincide fine.
9) utilize forecast model shown in the formula (9) to calculate dry rock bulk modulus data volume from step 8) gained elastic parameter inversion data volume;
10) utilize inverting elastic parameter data volume to ask the bulk modulus of saturated rock, and then ask the volume of voids modulus by formula (4).Fig. 5 is for calculating volume of voids modulus sectional view, and the curve that shows among the figure is the water saturation curve.From scheming as seen, it is good with water saturation curve corresponding relation to predict the outcome, and illustrates that the volume of voids modulus section of asking can be used for the fluid prediction.
11) utilize formula (10a) and (10b) calculated the water saturation data volume, as shown in Figure 6 from the volume of voids modulus.The curve that shows among the figure is a well logging interpretation water saturation curve.As seen from the figure, earthquake prediction water saturation and well logging interpretation saturation degree are coincide good.
Press S g=100-S wWith water saturation (S w) be converted into gas saturation (S g) section, and and factor of porosity
Figure B200910084537XD0000081
Section multiplies each other and obtains factor of porosity saturation degree product
Figure B200910084537XD0000082
Section, this section is described more directly perceived than water saturation section, more effective to reservoir gas-bearing property.
Fig. 7 is gas saturation and the factor of porosity product sectional view that is calculated, and the curve that shows among the figure is the water saturation curve of well logging interpretation.According to log well in advance A and well B place zone of interest of this section is the high yield gas-bearing formation, after 150,000 sides are produced in drilling well and gas testing confirmation well A gas testing daily, though not gas testing of well B, well logging interpretation is the high-quality payzone, and predicting the outcome has obtained the check of drilling well, illustrates that prediction effect is good.

Claims (4)

1. method based on the quantificationally predicting sandstone reservoir fluid saturation by combining well and seism of volume of voids modulus is characterized in that being achieved through the following technical solutions:
1) adopt arrangement greatly, small distance between receivers, high-density sampling is gathered high-quality geological data, seismic data is protected the width of cloth handle, and forms the common midpoint gather of fidelity;
2) gather compressional wave time difference, density, neutron, spontaneous potential or GR, resistivity conventional logging data, obtain logging trace, and processing obtains shale index, factor of porosity, saturation curves.
3) obtain longitudinal and transverse wave velocity curve from compressional wave time difference, shear wave slowness curve, calculate p-wave impedance according to p-and s-wave velocity, densimetric curve, the shear wave impedance, p-and s-wave velocity is than elastic parameter curve;
4) utilize multivariate statistical method to determine relation between shale index and factor of porosity and the elastic parameter respectively at well point place, set up the forecast model of following shale index and factor of porosity:
V Sh=f sh(L 1,L 2,L 3,...) (1)
φ=f φ(L 1,L 2,L 3,...) (2)
In the formula: V ShThe expression shale index, φ represents factor of porosity, f Sh, f φAll represent the multivariate function, independent variable is that elastic parameter, dependent variable are respectively shale index and factor of porosity, L 1, L 2, L 3... L nBe the elastic parameter curve;
5) adopt following formula to calculate volume of voids modulus and dry rock bulk modulus:
Figure F200910084537XC0000011
K d=K sat-K p (4)
In the formula: K SatThe bulk modulus of representing complete saturated rock, K dThe bulk modulus of expression dry rock, K pExpression hole bulk modulus, φ represents factor of porosity, K sThe bulk modulus of expression rock matrix, K fThe bulk modulus of expression pore fluid;
6) utilize dry rock bulk modulus and elastic parameter curve, set up the forecast model of dry rock bulk modulus:
K d=f kd(L 1,L 2,L 3,...) (5)
In the formula: L 1, L 2, L 3... L nBe elastic parameter curve, f KdExpression is that independent variable, dry rock bulk modulus are the multivariate function of dependent variable with the elastic parameter;
7) utilize multiple regression procedure to set up the forecast model of water saturation in effective reservoir:
S w=f sw(K p,L 4L 5,...) (6)
In the formula: L 1, L 2, L 3Expression elastic parameter curve, S wBe water saturation, f SwExpression is that independent variable, water saturation are the multivariate function of dependent variable with the elastic parameter;
8) utilize pre-stack seismic road collection and well-log information to carry out the pre-stack seismic inverting, obtain the inverting data volume of various elastic parameters, and utilize formula (1) and formula (2) to calculate shale index and factor of porosity data volume from the elastic parameter inversion data volume;
9) forecast model that utilizes step 5) to set up calculates the dry rock bulk modulus from step 8) gained elastic parameter inversion data volume;
10) utilize inverting elastic parameter data volume to ask the bulk modulus of saturated rock, and then ask volume of voids modulus K by formula (4) p
11) utilize formula (6) to ask water saturation S from the volume of voids modulus w
2. the method for the quantificationally predicting sandstone reservoir fluid saturation by combining well and seism based on the volume of voids modulus according to claim 1 is characterized in that it is to keep the AVO feature in processing procedure that the described guarantor's width of cloth of step 1) is handled.
3. the method for the quantificationally predicting sandstone reservoir fluid saturation by combining well and seism based on the volume of voids modulus according to claim 1, step 2) described collection data is gathered the shear wave slowness curve to the well more than a bite at least.
4. the method for the quantificationally predicting sandstone reservoir fluid saturation by combining well and seism based on the volume of voids modulus according to claim 1, the described shear wave velocity curve of step 3) utilizes Xu-White Model Calculation shear wave velocity curve when not gathering the well of shear wave slowness curve.
CN200910084537XA 2009-05-15 2009-05-15 Method for quantificationally predicting sandstone reservoir fluid saturation by combining well and seism Pending CN101887132A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN200910084537XA CN101887132A (en) 2009-05-15 2009-05-15 Method for quantificationally predicting sandstone reservoir fluid saturation by combining well and seism

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN200910084537XA CN101887132A (en) 2009-05-15 2009-05-15 Method for quantificationally predicting sandstone reservoir fluid saturation by combining well and seism

Publications (1)

Publication Number Publication Date
CN101887132A true CN101887132A (en) 2010-11-17

Family

ID=43073123

Family Applications (1)

Application Number Title Priority Date Filing Date
CN200910084537XA Pending CN101887132A (en) 2009-05-15 2009-05-15 Method for quantificationally predicting sandstone reservoir fluid saturation by combining well and seism

Country Status (1)

Country Link
CN (1) CN101887132A (en)

Cited By (28)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102156297A (en) * 2011-05-16 2011-08-17 中国石油大学(北京) Fluid substitution method based on sandstone reservoir post-stack seismic data
CN102426390A (en) * 2011-10-21 2012-04-25 中国石油大学(北京) Method for determining reserve volume of nonhomogeneous sandstone reservoir
CN102759757A (en) * 2012-07-06 2012-10-31 中国石油天然气股份有限公司 Diabase alteration zone reservoir prediction method based on combination of well and seism
CN103278850A (en) * 2013-05-24 2013-09-04 西安石油大学 Transverse wave time difference curve construction method based on coal rock industrial component physical volume model
CN103412335A (en) * 2013-08-20 2013-11-27 成都晶石石油科技有限公司 Method for predicting reservoir stratums through seismic phase bodies
CN103913429A (en) * 2014-03-21 2014-07-09 南京大学 Method for quantifying fluid saturation of two-phase flow in pore medium by light transmission method
CN104007467A (en) * 2014-04-16 2014-08-27 孙赞东 Pre-stack three-parameter inversion implementation reservoir stratum and fluid prediction method based on mixed norm regularization
CN104330822A (en) * 2014-10-23 2015-02-04 中国石油天然气股份有限公司 Method and device for determining remaining oil and gas distribution through coupling four-dimensional seismic inversion
CN104343445A (en) * 2014-10-29 2015-02-11 中国石油集团川庆钻探工程有限公司 Method for identifying gas reservoir, water-bearing stratum, reservoir with high gas-oil ratio and gas saturation
CN104570127A (en) * 2013-10-09 2015-04-29 中国石油化工股份有限公司 Method for simultaneously inverting porosity and shaliness by using seismic wave velocity
CN105093333A (en) * 2014-05-14 2015-11-25 中国石油天然气股份有限公司 Fracture-cavity type carbonate reservoir evaluation method based on earthquake-logging calibration
CN105116466A (en) * 2015-07-30 2015-12-02 中国石油天然气股份有限公司 Method and apparatus for determining oil field reservoir physical characteristics
CN106204302A (en) * 2016-07-05 2016-12-07 中国石油天然气股份有限公司 A kind of method calculating original water saturation and application
CN106353813A (en) * 2016-10-25 2017-01-25 中国海洋石油总公司 Method for identifying fluid properties based on array acoustic logging
CN107367762A (en) * 2017-06-27 2017-11-21 中国石油天然气股份有限公司 A kind of method and device for determining reservoir parameter
CN107422371A (en) * 2016-05-23 2017-12-01 中国石油化工股份有限公司 A kind of pre stack data AVO (A) attributes extraction method along line direction
CN108957534A (en) * 2018-06-12 2018-12-07 中国石油天然气股份有限公司 Gas saturation prediction technique and device
CN109001803A (en) * 2018-09-26 2018-12-14 中国石油化工股份有限公司 The method that seismic resolution is improved using fitting micronormal
CN110414129A (en) * 2011-10-18 2019-11-05 沙特阿拉伯石油公司 The modeling of 4D saturation degree
CN111323814A (en) * 2019-12-31 2020-06-23 中国石油天然气股份有限公司 Method for quantitatively determining water saturation of sandstone reservoir based on rock physical template
CN112213770A (en) * 2019-07-09 2021-01-12 中国石油天然气股份有限公司 Method and device for identifying reservoir hydrocarbon-containing sandstone based on log domain difference distribution characteristics
CN112213780A (en) * 2020-09-10 2021-01-12 同济大学 Reservoir parameter inversion method considering elastic impedance second-order gradient
CN112415588A (en) * 2019-08-22 2021-02-26 中国石油化工股份有限公司 Reservoir parameter oil-gas reserve calculation method and system based on three-dimensional seismic grid
CN112415616A (en) * 2019-08-22 2021-02-26 中国石油天然气集团有限公司 Deep-buried reservoir porosity inversion method and device
CN112649855A (en) * 2019-10-11 2021-04-13 中国石油化工股份有限公司 Three-dimensional gas saturation prediction method and system
CN112698399A (en) * 2020-12-02 2021-04-23 中国石油天然气股份有限公司 Gravel well seismic-logging linkage constraint efficient reservoir quantitative prediction method and system
CN113514891A (en) * 2021-04-27 2021-10-19 东北石油大学 Oil-water transition zone range determination method based on seismic inversion
CN113742655A (en) * 2021-09-18 2021-12-03 中国石油大学(北京) Method and device for determining oil saturation

Cited By (44)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102156297A (en) * 2011-05-16 2011-08-17 中国石油大学(北京) Fluid substitution method based on sandstone reservoir post-stack seismic data
CN110414129A (en) * 2011-10-18 2019-11-05 沙特阿拉伯石油公司 The modeling of 4D saturation degree
CN102426390A (en) * 2011-10-21 2012-04-25 中国石油大学(北京) Method for determining reserve volume of nonhomogeneous sandstone reservoir
CN102426390B (en) * 2011-10-21 2013-07-03 中国石油大学(北京) Method for determining reserve volume of nonhomogeneous sandstone reservoir
CN102759757A (en) * 2012-07-06 2012-10-31 中国石油天然气股份有限公司 Diabase alteration zone reservoir prediction method based on combination of well and seism
CN103278850A (en) * 2013-05-24 2013-09-04 西安石油大学 Transverse wave time difference curve construction method based on coal rock industrial component physical volume model
CN103278850B (en) * 2013-05-24 2014-08-20 西安石油大学 Transverse wave time difference curve construction method based on coal rock industrial component physical volume model
CN103412335A (en) * 2013-08-20 2013-11-27 成都晶石石油科技有限公司 Method for predicting reservoir stratums through seismic phase bodies
CN103412335B (en) * 2013-08-20 2015-12-09 成都晶石石油科技有限公司 A kind of method utilizing earthquake thing phase body predicting reservoir
CN104570127A (en) * 2013-10-09 2015-04-29 中国石油化工股份有限公司 Method for simultaneously inverting porosity and shaliness by using seismic wave velocity
CN103913429A (en) * 2014-03-21 2014-07-09 南京大学 Method for quantifying fluid saturation of two-phase flow in pore medium by light transmission method
CN104007467A (en) * 2014-04-16 2014-08-27 孙赞东 Pre-stack three-parameter inversion implementation reservoir stratum and fluid prediction method based on mixed norm regularization
CN104007467B (en) * 2014-04-16 2016-04-13 张远银 A kind of reservoir of the prestack three parametric inversion realization based on mixing norm regularization and fluid prediction method
CN105093333A (en) * 2014-05-14 2015-11-25 中国石油天然气股份有限公司 Fracture-cavity type carbonate reservoir evaluation method based on earthquake-logging calibration
CN104330822A (en) * 2014-10-23 2015-02-04 中国石油天然气股份有限公司 Method and device for determining remaining oil and gas distribution through coupling four-dimensional seismic inversion
CN104330822B (en) * 2014-10-23 2017-01-18 中国石油天然气股份有限公司 Method and device for determining remaining oil and gas distribution through coupling four-dimensional seismic inversion
CN104343445A (en) * 2014-10-29 2015-02-11 中国石油集团川庆钻探工程有限公司 Method for identifying gas reservoir, water-bearing stratum, reservoir with high gas-oil ratio and gas saturation
CN104343445B (en) * 2014-10-29 2017-07-14 中国石油集团川庆钻探工程有限公司 The method for recognizing gas-bearing formation and water layer, high gas-oil ratio (HGOR) reservoir and gas saturation
CN105116466B (en) * 2015-07-30 2017-12-19 中国石油天然气股份有限公司 A kind of method and apparatus for determining oilfield reservoir physical property characteristic
CN105116466A (en) * 2015-07-30 2015-12-02 中国石油天然气股份有限公司 Method and apparatus for determining oil field reservoir physical characteristics
CN107422371A (en) * 2016-05-23 2017-12-01 中国石油化工股份有限公司 A kind of pre stack data AVO (A) attributes extraction method along line direction
CN106204302A (en) * 2016-07-05 2016-12-07 中国石油天然气股份有限公司 A kind of method calculating original water saturation and application
CN106204302B (en) * 2016-07-05 2019-09-10 中国石油天然气股份有限公司 A kind of method and application calculating original water saturation
CN106353813A (en) * 2016-10-25 2017-01-25 中国海洋石油总公司 Method for identifying fluid properties based on array acoustic logging
CN106353813B (en) * 2016-10-25 2018-12-25 中国海洋石油总公司 Flow net model method based on array sonic log
CN107367762A (en) * 2017-06-27 2017-11-21 中国石油天然气股份有限公司 A kind of method and device for determining reservoir parameter
CN107367762B (en) * 2017-06-27 2019-06-11 中国石油天然气股份有限公司 A kind of method and device of determining reservoir parameter
CN108957534A (en) * 2018-06-12 2018-12-07 中国石油天然气股份有限公司 Gas saturation prediction technique and device
CN108957534B (en) * 2018-06-12 2020-02-14 中国石油天然气股份有限公司 Method and device for predicting gas saturation
CN109001803A (en) * 2018-09-26 2018-12-14 中国石油化工股份有限公司 The method that seismic resolution is improved using fitting micronormal
CN112213770A (en) * 2019-07-09 2021-01-12 中国石油天然气股份有限公司 Method and device for identifying reservoir hydrocarbon-containing sandstone based on log domain difference distribution characteristics
CN112213770B (en) * 2019-07-09 2023-11-28 中国石油天然气股份有限公司 Method and device for identifying hydrocarbon-containing sandstone in reservoir based on logarithmic domain difference distribution characteristics
CN112415588A (en) * 2019-08-22 2021-02-26 中国石油化工股份有限公司 Reservoir parameter oil-gas reserve calculation method and system based on three-dimensional seismic grid
CN112415616A (en) * 2019-08-22 2021-02-26 中国石油天然气集团有限公司 Deep-buried reservoir porosity inversion method and device
CN112649855A (en) * 2019-10-11 2021-04-13 中国石油化工股份有限公司 Three-dimensional gas saturation prediction method and system
CN112649855B (en) * 2019-10-11 2024-04-09 中国石油化工股份有限公司 Three-dimensional gas saturation prediction method and system
CN111323814A (en) * 2019-12-31 2020-06-23 中国石油天然气股份有限公司 Method for quantitatively determining water saturation of sandstone reservoir based on rock physical template
CN111323814B (en) * 2019-12-31 2022-07-05 中国石油天然气股份有限公司 Method for quantitatively determining water saturation of sandstone reservoir based on rock physical template
CN112213780A (en) * 2020-09-10 2021-01-12 同济大学 Reservoir parameter inversion method considering elastic impedance second-order gradient
CN112698399B (en) * 2020-12-02 2023-08-22 中国石油天然气股份有限公司 Method and system for quantitatively predicting efficient reservoir based on vibration measurement linkage constraint of conglomerate well
CN112698399A (en) * 2020-12-02 2021-04-23 中国石油天然气股份有限公司 Gravel well seismic-logging linkage constraint efficient reservoir quantitative prediction method and system
CN113514891A (en) * 2021-04-27 2021-10-19 东北石油大学 Oil-water transition zone range determination method based on seismic inversion
CN113742655B (en) * 2021-09-18 2023-06-06 中国石油大学(北京) Method and device for determining oil saturation
CN113742655A (en) * 2021-09-18 2021-12-03 中国石油大学(北京) Method and device for determining oil saturation

Similar Documents

Publication Publication Date Title
CN101887132A (en) Method for quantificationally predicting sandstone reservoir fluid saturation by combining well and seism
CN105134195A (en) Shale gas reservoir quality evaluation method based on logging information
CN106368691B (en) Three-dimensional abnormal pore pressure prediction method based on rock physics seismic information
CN102736107B (en) Energy constraint heterogeneous reservoir thickness identification system
KR101415199B1 (en) Method for estimating slowness, Young's modulus, Poisson's ratio and brittleness of horizontal drilling zone in sweet spot range using resistivity and density logging data in shale play
CN108363100A (en) Based on the coal bed gas dessert area Seismology and Geology recognition methods for sealing unit and rough set up for safekeeping
Sil et al. Dynamic site characterization and correlation of shear wave velocity with standard penetration test ‘N’values for the city of Agartala, Tripura state, India
Noorlandt et al. Characterisation of ground motion recording stations in the Groningen gas field
CA2867583C (en) Fracking method for fracking intervals of a horizontal drilling zone in a sweet spot range based on measurements of resistivity and neutron logging data in the horizontal drillingzone
KR101415196B1 (en) Method for selecting fracking intervals of horizontal drilling zone in sweet spot range using resistivity and density logging data in shale play
CN107436452A (en) Hydrocarbon source rock Forecasting Methodology and device based on probabilistic neural network algorithm
CN104632202A (en) Method and device for determining dry clay three-porosity logging parameter values
Yilmaz et al. A case study of seismic zonation in municipal areas
CN107843927A (en) Shale formation pressure prediction method and device based on well shake joint speed
CN110424956A (en) Evaluation unit saves coefficient weights quantization assignment method in shale oil Resources calculation
CN103777245A (en) Method of quantitative evaluation on oil and gas accumulation conditions based on seismic data
CN105064986A (en) Method for building reservoir four-property relationship spectrum by using conventional well detection and logging information
CN101852864B (en) Method for processing mass seismic data by using surface consistent statistical spectrum analysis technology
CN109061752A (en) A kind of resistivity curve bearing calibration on the stratum containing grey matter
CN104749618A (en) Shale low-angle crack post-stack probability quantitative characterization method
Zou et al. Seismic studies of coal bed methane content in the west coal mining area of Qinshui Basin
Efetobore et al. Reservoir Characterization and Fluid Analysis of the Majosa Field in the Niger Delta
Nurhandoko Shale Gas Reservoir Characterization Strategy
Babak et al. Facies trend modeling based on geobody-based binning of seismic impedance
CN110297280A (en) A kind of prediction technique of carbonate rock superpressure spatial distribution characteristic

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C02 Deemed withdrawal of patent application after publication (patent law 2001)
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20101117