CN106054248A - 一种基于大面积致密储层地震岩石物理反演方法 - Google Patents
一种基于大面积致密储层地震岩石物理反演方法 Download PDFInfo
- Publication number
- CN106054248A CN106054248A CN201610563338.7A CN201610563338A CN106054248A CN 106054248 A CN106054248 A CN 106054248A CN 201610563338 A CN201610563338 A CN 201610563338A CN 106054248 A CN106054248 A CN 106054248A
- Authority
- CN
- China
- Prior art keywords
- rock
- template
- reservoir
- well
- work area
- 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
Links
- 239000011435 rock Substances 0.000 title claims abstract description 173
- 238000000034 method Methods 0.000 title claims abstract description 50
- 239000012530 fluid Substances 0.000 claims abstract description 71
- 238000012360 testing method Methods 0.000 claims abstract description 24
- 230000008859 change Effects 0.000 claims abstract description 15
- 238000010206 sensitivity analysis Methods 0.000 claims abstract description 12
- 239000007789 gas Substances 0.000 claims description 51
- 239000011148 porous material Substances 0.000 claims description 30
- 230000035945 sensitivity Effects 0.000 claims description 18
- 238000010008 shearing Methods 0.000 claims description 18
- 238000002604 ultrasonography Methods 0.000 claims description 16
- 238000009826 distribution Methods 0.000 claims description 14
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 14
- 238000004458 analytical method Methods 0.000 claims description 9
- 230000004044 response Effects 0.000 claims description 9
- 238000012937 correction Methods 0.000 claims description 8
- 229910052500 inorganic mineral Inorganic materials 0.000 claims description 8
- 239000011707 mineral Substances 0.000 claims description 8
- 238000001615 p wave Methods 0.000 claims description 8
- 239000004575 stone Substances 0.000 claims description 8
- 238000004364 calculation method Methods 0.000 claims description 7
- 239000006185 dispersion Substances 0.000 claims description 7
- 230000008569 process Effects 0.000 claims description 7
- 229920006395 saturated elastomer Polymers 0.000 claims description 7
- 239000011159 matrix material Substances 0.000 claims description 6
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 claims description 6
- 239000000758 substrate Substances 0.000 claims description 6
- 239000000470 constituent Substances 0.000 claims description 5
- 238000009499 grossing Methods 0.000 claims description 5
- 238000004519 manufacturing process Methods 0.000 claims description 5
- 239000000203 mixture Substances 0.000 claims description 5
- 238000012545 processing Methods 0.000 claims description 4
- 238000005520 cutting process Methods 0.000 claims description 3
- 230000004927 fusion Effects 0.000 claims description 3
- 239000003345 natural gas Substances 0.000 claims description 3
- 230000035699 permeability Effects 0.000 claims description 3
- 241000208340 Araliaceae Species 0.000 claims description 2
- 238000012935 Averaging Methods 0.000 claims description 2
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 claims description 2
- 235000003140 Panax quinquefolius Nutrition 0.000 claims description 2
- 238000004422 calculation algorithm Methods 0.000 claims description 2
- 238000001514 detection method Methods 0.000 claims description 2
- 238000006073 displacement reaction Methods 0.000 claims description 2
- 238000005553 drilling Methods 0.000 claims description 2
- 238000000605 extraction Methods 0.000 claims description 2
- 235000008434 ginseng Nutrition 0.000 claims description 2
- 239000000463 material Substances 0.000 claims description 2
- 239000004568 cement Substances 0.000 claims 1
- 230000007613 environmental effect Effects 0.000 claims 1
- 238000012216 screening Methods 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 13
- 229910000514 dolomite Inorganic materials 0.000 description 10
- 239000010459 dolomite Substances 0.000 description 10
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 7
- 238000005516 engineering process Methods 0.000 description 7
- 229930195733 hydrocarbon Natural products 0.000 description 7
- 150000002430 hydrocarbons Chemical class 0.000 description 7
- 230000007246 mechanism Effects 0.000 description 6
- 239000004215 Carbon black (E152) Substances 0.000 description 5
- 238000002474 experimental method Methods 0.000 description 5
- 230000015572 biosynthetic process Effects 0.000 description 4
- 230000007812 deficiency Effects 0.000 description 3
- 235000013399 edible fruits Nutrition 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 230000003247 decreasing effect Effects 0.000 description 2
- 239000004744 fabric Substances 0.000 description 2
- 238000011002 quantification Methods 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 238000003860 storage Methods 0.000 description 2
- 208000035126 Facies Diseases 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000013499 data model Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 238000002637 fluid replacement therapy Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000008595 infiltration Effects 0.000 description 1
- 238000001764 infiltration Methods 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 238000004445 quantitative analysis Methods 0.000 description 1
- 238000010200 validation analysis Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/282—Application of seismic models, synthetic seismograms
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/61—Analysis by combining or comparing a seismic data set with other data
- G01V2210/616—Data from specific type of measurement
- G01V2210/6169—Data from specific type of measurement using well-logging
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种基于大面积致密储层地震岩石物理反演方法,该方法包括:建立多尺度岩石物理模型;开展声波测井与超声波两种尺度下岩石物理参数的流体敏感性分析,筛选出多个观测尺度下对孔隙度、含气饱和度最为敏感的岩石物理参数;构建单井岩石物理模板,优选其标准模板;考虑储层地质特征的横向变化与非均质性,在大工区内根据所有井试气情况对岩石物理模板的输入参数进行微调,优化整个工区并建立三维工区岩石物理模板数据体,将其与地震叠前反演结合计算目的层孔隙度、饱和度;对结果进行平滑处理最终输出储层参数反演数据体。通过与实际地层测井解释、试气结果的对比,证实了所预测结果的正确性及本发明的有效性。
Description
技术领域
本发明属于地震岩石物理反演技术领域,特别是涉及一种基于大面积致密储层地震岩石物理反演方法。
背景技术
目前,利用地震资料直接进行储层预测与油气特征的识别,已发展成为勘探地球物理学领域的研究热点,但其研究难点在于地震数据是地下储层多方面性质的综合反映,反问题的欠定性导致了地震解释、地震预测存在多解性。地震岩石物理作为沟通地面地震数据与地下储层参数的桥梁,是实现含油气储层定量解释的理论与实验基础。岩石物理驱动下的地震反演,有助于准确认识地下储层的油气特征与分布规律,从而应对地震反演问题的多解性与局限性。岩石物理学研究可建立岩石弹性参数与储层物性、孔隙流体参数之间的定量关系,将传统岩石物理实验研究的岩心尺度拓展到油田尺度,建立不同观测尺度(岩心、测井、地震)的岩石物理模型,实现不同尺度、多学科数据的融合(Tang,2008)。针对地震、测井以及超声波实验三个观测尺度,目前多数地球物理学者采用孔隙流体流动机制所诱导的波速频散机理解释不同尺度间观测结果的差异,三种类型的孔隙流体流动机制,分别为宏观尺度孔隙流体流动机制(Gassmann,1951;Biot,1956,1962)、颗粒尺度孔隙流体流动机制(Mavko,1975;Dvorkin,1993;Gurevich等,2010;Vinci等,2014;Papageorgiou等,2015)、中观尺度(大于孔隙尺寸,但远小于波长尺寸)流体流动机制(White,1975;Pride,2003;Ba等,2011,2012;Sun等,2015)。鉴于碳酸岩储层孔隙结构复杂,横向非均质性强,现有技术中的单一孔隙结构、单一尺度数据建立的常规岩石物理模板,无法适用于大型工区。
中国专利201210335739.9公开了“一种用于储层烃类检测的多尺度岩石物理图版方法及装置”,该方法包括:获得储层岩石基质模型,获得储层岩石骨架模型,根据储层岩石基质模型、储层岩石骨架模型和流体参数获得含流体的岩石模型,并根据含流体的岩石模型获得多尺度的初始储层岩石物理图版;根据岩石物理实验数据,在多个频段下校正初始储层岩石物理图版:根据测井解释结果,在声波测井频段校正初始储层岩石物理图版;根据地震解释结果,对初始储层岩石物理图版进行校正;在通过校正后,输出最终的储层岩石物理图版,以进行储层岩石物理参数的反演。该方法可提高岩石物理图版的精度并实现对岩石参数及流体饱和度的定量反演。但该方法主要考虑了单一类型的孔隙结构联合多尺度波资料建模,所存在的明显不足是,不能适用于地层孔隙结构类型发生变化的非均质性储层。
中国专利201310752436.1公开了“一种基于孔隙结构特征的碳酸盐岩储层的气藏预测方法及系统”,该方法包括:采集碳酸盐岩储层目标储层段的岩石样本;对岩石样本进行地质薄片鉴定,得到岩石基本参数,岩石基本参数包括岩石成分、孔隙形状、面孔率以及沉积相带;对岩石样本进行孔渗测量,得到孔渗基本参数,孔渗基本参数包括孔隙度、渗透率以及密度;根据岩石基本参数、孔渗基本参数以及微分等效介质模型构建岩石干骨架模型;对岩石干骨架模型进行流体替换,生成岩石物理图版;获取碳酸盐岩储层的叠前地震反演数据;将叠前地震反演数据与岩石物理图版进行交会,得到碳酸盐岩储层的孔隙度与含气饱和度的预测结果。实现了精确的气藏定量化预测。该方法相对于专利201210335739.9,适度考虑了孔隙结构的横向非均质性变化特征,分别分析了孔隙型、裂缝型储层两类结构特征,可适用于小面积工区的反演应用,但不能适用于大面积工区地质特征横向变化剧烈及多井联合反演储层参数的技术需求,如:一是没有联合多井资料建立工区标准模板以描述地层一般性特征,并在工区各个坐标点进行微调;二是不能实现三维工区模板数据体的制作,仅根据少量参考井对结构进行定性划分,因此对非均质性描述的精度较低;三是在储层烃类预测之前,没有开展多尺度资料的流体敏感性分析。
综上所述,如何克服现有技术所存在的不足,已成为当今地震岩石物理反演技术领域中亟待解决的重点难题之一。
发明内容
本发明的目的是为克服现有技术的不足而提供一种基于大面积致密储层地震岩石物理反演方法,本发明能够有效的利用多井资料、多尺度波形响应数据,切实提高烃类地震预测的精度;将岩石物理模型与叠前地震数据反演相结合,可实现储层的定量预测。
根据本发明提出的一种基于大面积致密储层地震岩石物理反演方法,其特征在于,包括如下具体步骤:
步骤101,基于孔隙弹性理论预测波响应频散,建立多尺度岩石物理模型的参量因素以联系多尺度数据;所述建立多尺度岩石物理模型的参量因素,是以岩石的矿物成分、孔隙结构及所处的地层环境会对岩石的波响应特征产生影响为依据,确定包括温度和压力为储层环境因素、包括矿物组分、孔隙形状、泥质含量和孔隙结构为储层岩性因素、包括流体黏性和气水斑状饱和为储层流体因素;
步骤102,基于模型及部分井试气情况,对测井解释结果进行分析和校正,开展声波测井与超声波两种尺度下岩石物理参数的流体敏感性分析,筛选出多个观测尺度下对孔隙度、含气饱和度最为敏感的岩石物理参数;所述超声波和声波测井两种尺度下岩石物理参数是指弹性参数及弹性参数的组合;所述弹性参数至少包括以下物理量:纵波速度Vp、横波速度Vs、纵波阻抗Zp、横波阻抗Zs、纵横波速度比Vp/Vs、拉梅常数λ、剪切模量μ、拉梅常数和密度乘积λρ、拉梅常数和剪切模量乘积λμ、准压力PR、剪切模量和密度乘积μρ;所述流体敏感性分析包括:超声波尺度下测量纵波速度Vp和横波速度Vs、以及气和水的变饱和度的情况下的波速,即Vp/Vs与波阻抗交会图;
步骤103,基于各参考井资料构建的单井岩石物理模板,优选各单井模板制作工区标准模板;所述工区标准模板是指以优选的敏感参数λρ为纵坐标,以纵波阻抗为横坐标;
步骤104,基于储层地质特征的横向变化与非均质性,在工区内根据所有井试气情况在各井坐标处对岩石物理模板的输入参数进行微调,并在整个工区内进行优化,建立三维工区岩石物理模板数据体,并将三维工区岩石物理模板数据体与地震叠前反演结合计算目的层孔隙度、饱和度;在整个工区开展大面积的三维岩石物理模板参数反演,并对反演结果进行平滑处理,最终输出储层参数反演数据体,进而实现储层孔隙度及饱和度的定量解释;所述建立三维工区岩石物理模板数据体是指将待反演、解释的三维地震数据体根据工程需要进行切割、整理,开展叠前三维地震反演,并进行储层孔隙度、饱和度的反演计算。
本发明提出的一种基于大面积致密储层地震岩石物理反演方法的进一步优选方案是:
步骤101所述建立多尺度岩石物理模型的建模,包括计算岩石基质弹性模量和岩石骨架弹性模量,采用Voigt-Reuss-Hill平均方程求取基质的等效弹性模量:
其中MVRH是矿物基质的弹性模量,fi和Mi分别是第i种组分的体积含量和弹性模量,N是矿物组分的总数量;采用微分等效介质(DEM)理论(Mavko,1998),求取白云岩的干岩石骨架的体积模量与剪切模量
(1-y)d/dy[K*(y)]=(K2-K*(y))P(*2)(y) (2a)
(1-y)d/dy[μ*(y)]=(μ2-μ*(y))Q(*2)(y) (2b)
初始条件是K*(0)=K1和μ*(0)=μ1,其中K1,μ1是初始主相材料的体积模量和剪切模量(相1),K2,μ2是逐渐加入的包含物的体积模量和剪切模量(相2),y为相2的含量,P(*2)和Q(*2)与嵌入体形状有关。
步骤101所述建立多尺度岩石物理模型的建模还包括考虑储层环境下的岩石,采用vanderWaals方程,估算储层流体在不同温度、压力条件下的天然气的密度和体积模量。
步骤101所述建立多尺度岩石物理模型的建模还包括考虑孔隙流体的非均质分布,忽略孔隙结构的非均质性,利用Biot-Rayleigh方程来实现饱和流体岩石纵横波速度的预测,同时还考虑弹性波速度频散的影响,以实现不同尺度多学科数据间的融合,Biot-Rayleigh方程具体形式如下:
其中u=[u1,u2,u3],分别代表三种组分的空间矢量位移(骨架、流体1、流体2),下标1、2、3表示矢量空间的三个方向;表示地震波激励过程中产生的局域流体变形增量,eij、和分别为:
其中:x1,x2与x3分别表示三个方向的坐标。φ1和φ2表示两类孔隙的绝对孔隙度,岩石的总孔隙度φ=φ1+φ2;φ10与φ20分别表示两个区域内的局部孔隙度,如岩石内部仅含有一种骨架,但却饱和有两种流体,则φ10=φ20=φ;假设φ1代表水孔(背景相流体),φ2代表气孔(包含物相流体),那么φ1/φ就是含水饱和度,φ2/φ就是含气饱和度。ρf1和η1表示背景相流体的密度与黏度,ρf2和η2表示包含物相流体的密度与黏度;R0表示气泡半径,κ10表示岩石渗透率;弹性参数A、N、Q1、R1、Q2与R2,密度参数ρ11、ρ12、ρ13、ρ22与ρ33,耗散参数b1与b2的确定方法。
步骤102还包括对比分析超声波尺度和测井尺度下岩石物理参数敏感性规律、次序基本一致,对孔隙流体最为敏感性参数为λ和λρ。
步骤103对所述工区标准岩石物理模板格点的横、纵坐标进行校正,是指保证工区标准岩石物理模板的描述结果基本涵盖标准井所有数据,则校正后的工区标准岩石物理模板表示为:
假设各单井岩石物理模板在二维交会图中的任意一个模板格点位置可根据其横、纵坐标用(Mk(i,j),Nk(i,j))表示,其中i、j分别对应于孔隙度和饱和度的渐变(i=1,2…11对应于孔隙度0.02渐变到0.12;j=1,2…11对应于饱和度0渐变到100%),k表示第k个参考井,则标准模板在该格点位置的值(Ms(i,j),Ns(i,j))可表示为:
其中:A(k)表示第k口井的权重,B(i,j)、C(i,j)分别表示对应于i、j的模板格点上基于综合各井参考数据对工区标准岩石物理模板格点的横、纵坐标进行的校正。
步骤104所述储层参数岩石物理反演,是指基于工区标准岩石物理模板,在各井位置开展地震反演与流体检测测试,具体方法为:
从三维地震数据体中抽取二维过井线,对井坐标附件的目的层进行储层、流体参数估算,对比估算结果与已知钻井、试气结果,适度调整岩石物理模板,试反演结果、解释结论与实际试气情况达到最佳吻合;在第k口井附近,工区标准岩石物理模板经二维地震反演调试后的地震反演模板可用(Mk′(i,j),Nk′(i,j))表示;
在工区内,基于各井的地震反演模板(Mk′(i,j),Nk′(i,j)),开展全区优化,生成工区岩石物理模型的三维数据体,具体对应于工区内每一个坐标位置(x,y),都存在一个独立的地震岩石物理模型,其模板(M3D(x,y,i,j),N3D(x,y,i,j))的确定依赖于各井位置的地震反演模板:
M3D(x,y,i,j)=Mk′(i,j),若x=xk,y=yk (5a)
N3D(x,y,i,j)=Nk′(i,j),若x=xk,y=yk (5b)
其中:(xk,yk)为第k口井坐标,L为总井数,Q(x,y,k)为计算工区内坐标(x,y)处模板时所采用的第k口井的权重系数,可由下式确定:
基于以上方法,即可兼顾工区内各井资料,基于各地理位置的测井观测控制储层参数的地震反演及解释,空间位置越靠近某口参考井,则受该口井的控制因素影响越为显著;(M3D(x,y,i,j),N3D(x,y,i,j))即工区岩石物理模型三维数据体,制作过程中按地震数据.sgy标准格式首先沿xline方向,继而沿inline方向逐条xline线进行处理。
步骤104所述采用加权平均法对反演的孔隙度、饱和度进行平滑处理,以减弱反演数据的野值、奇异值影响;
假定离目标点越近反演结果对目标点造成的影响越大,根据二维正态分布定义了三种加权模板,如式(7)所示;T1与T2是通过二维高斯离散模板演变而成,而T3就是3×3邻域的二维高斯模板,模板中最大权值为目标点所在位置,可根据最小均匀度所在邻域位置,将模板进行适当调整,以形成与该邻域相关的模板:
以T1为例,在邻域内对二维反演剖面data(i,j)进行平滑处理,具体算法可表示为:
其中dataS即为平滑处理后目标点的数值。
本发明的实现原理是:本发明基于孔隙弹性理论预测波响应频散,建立多尺度岩石物理模型以联系多尺度数据;基于模型及部分井试气情况,对测井解释结果进行分析和校正,开展声波测井与超声波尺度下岩石参数的流体敏感性分析,筛选出多个观测尺度下对孔隙度、含气饱和度最为敏感的岩石物理参数;基于各参考井资料构建的单井岩石物理模板,优选各单井模板制作工区的标准模板;考虑储层地质特征的横向变化与非均质性,在工区内根据所有井试气情况在各井坐标处对岩石物理模板的输入参数进行微调,并在整个工区内进行优化,建立三维工区岩石物理模板数据体,将其与地震叠前反演结合实现大面积工区的岩石物理参数定量计算。
本发明与现有技术相比其显著优点在于:
一是本发明在本领域首创了基于大面积致密储层地震岩石物理反演方法,具有大面积推广应用的重要现实意义。
二是本发明能够有效的利用多井资料、多尺度波形响应数据,切实提高烃类地震预测的精度;将岩石物理模型与叠前地震数据反演相结合,可实现储层的定量预测。
三是本发明开展岩石物理参数多尺度流体敏感性分析,基于多尺度观测选取最敏感的岩石物理参数用于储层烃类预测。
四是本发明结合多井资料建立了工区目的层的标准模板,可描述目的层的一般性地质及岩石物理特征,再将模板于各已知井进行地震反演测试,建立了三维岩石物理模板数据体,在一般性地质特征的基础上考虑了横向的非均质变化,从而可适用于大面积工区的地震参数的岩石物理反演。
附图说明
图1为本发明提出的一种基于大面积致密储层地震岩石物理反演方法的流程示意图。
图2(包括图2-1、图2-2)为白云岩骨架体积模量(a)和剪切模量(b)随孔隙度的变化关系计算结果示意图;其中:图2-1为白云岩骨架体积模量(a)随孔隙度的变化关系计算结果示意图;图2-2为白云岩骨架剪切模量(b)随孔隙度的变化关系计算结果示意图。
图3为MX区204井岩石物理参数的交汇示意图。
图4为MX9井岩石物理解释模板示意图。
图5为MX12井岩石物理解释模板示意图。
图6为白云岩样品Vp/Vs与波阻抗交汇示意图。
图7为不同岩样岩石物理参数的流体敏感性分析示意图。
图8(包括图8-1、图8-2)为不同观测尺度下岩石参数的流体敏感性分析示意图;其中:图8-1为超声波尺度下各岩石物理参数的流体敏感性分析;图8-2为声波测井尺度各岩石物理参数的流体敏感性分析。
图9为工区标准岩石物理模板及纵波阻抗、λρ交汇图(基于MX203、204、26、27井)示意图。
图10为整个工区三维岩石物理模型数据体制作示意图。
图11为平滑前、后过MX8井2维线孔隙度反演剖面对比示意图。
图12为大面积工区储层平均孔隙度反演结果示意图。
图13为大面积工区储层平均含气饱和度反演结果示意图。
图14(包括图14-1、图14-2)为过MX13号井储层孔隙度和含气饱和度反演结果示意图;其中:图14-1为储层孔隙度反演结果;图14-2为储层含气饱和度反演结果。
图15(包括图15-1、图15-2)为过MX17号井储层孔隙度和含气饱和度反演结果示意图;其中:图15-1为储层孔隙度反演结果;图15-2为储层含气饱和度反演结果。
具体实施方式
下面结合附图和实施例对本发明的具体实施方式作进一步的详细说明。
如图1所示,本发明提出的一种基于大面积致密储层地震岩石物理反演方法包括如下具体步骤:
步骤101,基于孔隙弹性理论预测波响应频散,建立多尺度岩石物理模型的参量因素以联系多尺度数据;所述建立多尺度岩石物理模型的参量因素,是以岩石的矿物成分、孔隙结构及所处的地层环境会对岩石的波响应特征产生影响为依据,确定包括温度和压力为储层环境因素、包括矿物组分、孔隙形状、泥质含量和孔隙结构为储层岩性因素、包括流体黏性和气水斑状饱和为储层流体因素;
步骤102,基于模型及部分井试气情况,对测井解释结果进行分析和校正,开展声波测井与超声波两种尺度下岩石物理参数的流体敏感性分析,筛选出多个观测尺度下对孔隙度、含气饱和度最为敏感的岩石物理参数;所述超声波和声波测井两种尺度下岩石物理参数是指弹性参数及弹性参数的组合;所述弹性参数至少包括以下物理量:纵波速度Vp、横波速度Vs、纵波阻抗Zp、横波阻抗Zs、纵横波速度比Vp/Vs、拉梅常数λ、剪切模量μ、拉梅常数和密度乘积λρ、拉梅常数和剪切模量乘积λμ、准压力PR、剪切模量和密度乘积μρ;所述流体敏感性分析包括:超声波尺度下测量纵波速度Vp和横波速度Vs、以及气和水的变饱和度的情况下的波速,即Vp/Vs与波阻抗交会图;
步骤103,基于各参考井资料构建的单井岩石物理模板,优选各单井模板制作工区标准模板;所述工区标准模板是指以优选的敏感参数λρ为纵坐标,以纵波阻抗为横坐标;
步骤104,基于储层地质特征的横向变化与非均质性,在工区内根据所有井试气情况在各井坐标处对岩石物理模板的输入参数进行微调,并在整个工区内进行优化,建立三维工区岩石物理模板数据体,并将三维工区岩石物理模板数据体与地震叠前反演结合计算目的层孔隙度、饱和度;在整个工区开展大面积的三维岩石物理模板参数反演,并对反演结果进行平滑处理,最终输出储层参数反演数据体,进而实现储层孔隙度及饱和度的定量解释;所述建立三维工区岩石物理模板数据体是指将待反演、解释的三维地震数据体根据工程需要进行切割、整理,开展叠前三维地震反演,并进行储层孔隙度、饱和度的反演计算。
本发明实施例基于大面积致密储层地震岩石物理反演方法,针对碳酸盐岩储层的强非均质性特征,分析储层的矿物组成,建立储层的干岩石骨架模型,根据储层环境建立储层流体模型,并构建了多尺度岩石物理模型;基于多尺度岩石物理模型,提出了大面积工区储层、流体参数岩石物理反演的工业化技术流程。
本发明实施例基于大面积致密储层地震岩石物理反演方法,首先开展测井资料分析及精细解释,制作单井岩石物理模型及模板;其次开展声波测井、超声波尺度下的岩石物理参数流体敏感性分析,结果表明:随着孔隙度的降低,各参数对流体的敏感性明显下降,但各参数之间的敏感性次序基本不变;两种尺度下岩石物理参数的敏感性次序基本一致,但同一参数在不同观测尺度下对孔隙流体的敏感性有所不同。基于多尺度敏感性分析优选出对含气饱和度最为敏感的参数λρ和λ。
本发明实施例基于大面积致密储层地震岩石物理反演方法,为描述目的层地质特征的一般性规律,优选各单井模拟结果制作工区的标准岩石物理模板,结合各过井二维线的地震反演测试,在大面积工区制作了三维岩石物理模型数据体,将其与地震叠前反演结合,估算储层孔隙度与含气饱和度,并采用了平滑矩阵对反演结果数据体进行了平滑处理;将反演结果与目的层的测井解释、试油结论对比显示,所预测结果与实际地层孔隙度、天然气产能状况吻合很好。
本发明实施例基于大面积致密储层地震岩石物理反演方法,结合岩石物理分析,研发了大面积工区储层、流体参数地震定量预测的工业化技术流程,并成功应用于中国西部某致密白云岩储层。考虑到各个工区的地质特征、储层类型、流体分布情况存在差异,为保证岩石物理模型的适用性与准确性,必须从储层实际情况出发,针对不同岩性、不同孔构以及不同流体开展分析,并结合实验数据、测井数据以及地震反演数据进行校正,从而实现地震岩石物理反演方法的有效应用。
以下以四川盆地致密白云岩气藏为例,说明基于大面积致密储层岩石地震岩石物理反演方法的实现过程。
如图2所示,为本发明实施例提供的白云骨架体积模量(a)和剪切模量(b)随孔隙度的变化关系计算结果。采用DEM方法计算的含不同纵横比(α)硬币状裂缝的白云岩骨架模量随孔隙度的变化关系曲线,结果中同孔隙度情况下随孔隙纵横比增大,骨架模量增大,同纵横比情况下随孔隙度的增大,骨架模量减小。
如图3所示,为本发明实施例给出了工区某井的测井岩石物理参数交汇图。可根据测井数据的分布规律、分散程度编辑测井曲线,去除野值,并将测井解释结论对比实际试气结果,开展测井数据精细解释,而精细解释所得的孔隙度、含气性(饱和度)则可用于标定岩石物理模型输入参数。
如图4和图5所示,为本发明实施例提供的MX9井、MX12井的岩石物理解释模板与实际储层孔隙度解释结果的对比。对比显示,由于碳酸岩储层类型的多样性,孔隙流体分布的非均质性,造成基于单一孔隙结构、分布模式建立的岩石物理模板,难以覆盖较分散的实际数据,不能有效辅助识别岩性与流体,并且由于井与井之间存在横向非均质性,造成各单井模板之间存在差异,因此基于单井建立的岩石物理模板,难以表征整个工区的有利储层,这个问题在大面积工区实际应用中尤其显著。
如图6所示,为本发明实施例提供的5块岩心白云岩样品Vp/Vs与波阻抗交会图。
如图7所示,为本发明实施例,为调查岩样孔隙度对岩石参数的影响,根据5块孔隙度分别为0.051、0.0534、0.0547、0.1208与0.1228的白云岩样品实测数据,对比分析超声波尺度岩样弹性参数及弹性参数组合对储层含气性的敏感性,统计结果表明,如图8所示,随着孔隙度下降,各参数敏感性明显下降,各参数之间的敏感性次序有微弱变化,但基本趋势不变;λ和λρ对含气饱和度的变化最为敏感,μ和Vs对含气饱和度的敏感性最差。
如图8所示,为本发明实施例,为调查不同观测尺度下岩石物理参数对孔隙流体的敏感性,从上述5块白云岩样品及7口井资料出发,对比分析超声波尺度与测井尺度的敏感性计算结果,如图9所示。
对比显示,两种尺度下的岩石物理参数敏感性规律、次序基本一致,λ和λρ对孔隙流体的变化最为敏感,少数参数在不同观测尺度下对流体的敏感性有微弱变化,如μ在超声波尺度对孔隙流体最不敏感,但在声波测井尺度对流体敏感性有所增强,这些尺度之间的变化增加了储层流体识别难度。通过上述对比分析,优选出两个观测尺度下对流体最为敏感的参数λ和λρ,这是后续储层预测与流体识别的基础。
如图9所示,为本发明实施例提供的工区标准岩石物理模板及纵波阻抗、λρ交会图。将测井解释结果投影到校正后的岩石物理模板上,结果显示,工区内测井解释的高含水储层,在模板上分布在饱水线附近,其孔隙度分布范围为2.8%~12%;解释的高含气储层,在模板上分布在饱气线附近,其孔隙度分布范围为3%~9.2%;解释的非储层,在模板上分布在低孔隙度区域,其孔隙度基本在4%以下。采用此模板,基本可描述工区所有重点井重点地层的试气情况,并反映了目的层孔隙度分布的一般性特征。
如图10所示,为本发明实施例(M3D(x,y,i,j),N3D(x,y,i,j))即工区岩石物理模型三维数据体,制作过程中按地震数据.sgy标准格式首先沿xline方向,继而沿inline方向逐条xline线进行处理。
如图11所示,为本发明实施例,平滑前、后过MX8井的孔隙度反演剖面,对比平滑前、后可知,孔隙度反演剖面在空间连续性上有明显改善。
如图12所示,为本发明实施例,为千余平方公里大面积工区储层平均孔隙度反演结果,反演结果显示:目的层储层具有低孔低渗特征,储层平均孔隙度分布范围为0~0.045,优质储层集中分布于工区中部的构造高部位。
如图13所示,为本发明实施例,为储层平均含气饱和度反演结果,预测结果显示含气富集区集中分布于上、下两个区,这一结果与前期针对该地区的石油地质学研究认识一致,此外反演结果也吻合了各井的试气产能情况,如上气下水的MX203井、高产气井MX8井,反演结果均符合实际情况。
如图14,图15所示,为本发明实施例,如图14、15所示,在三维工区的岩石物理反演数据体中选取2条过井地震反演剖面,将反演结果与测井解释结果、试气报告进行了对比。
图14为过MX13井储层孔隙度与含气饱和度预测结果,对比发现,地震反演的孔隙度与测井资料解释的岩心孔隙度基本一致,且地震反演有效识别了目的层中部的优质储层段,含气饱和度结果显示,MX13井在目的层有较强的含气性,与试气结果吻合(该井目的层日产气达到128.84万方/天),而在剖面的构造低部位有含水性显示,反演结果初步展示了气、水界面直接地震识别的可能性。图15为过MX17井储层孔隙度与含气饱和度反演结果,结果显示:地震反演结果精细识别了测井解释结论中给出的三个含气薄层,且吻合试气结果(该井日产气53.2万方/天)。
本发明的具体实施方式中未涉及的说明属于本领域公知的技术,可参考公知技术加以实施。
本发明经反复试验验证,取得了满意的试用效果。
以上具体实施方式及实施例是对本发明提出的一种基于大面积致密储层地震岩石物理反演方法技术思想的具体支持,不能以此限定本发明的保护范围,凡是按照本发明提出的技术思想,在本技术方案基础上所做的任何等同变化或等效的改动,均仍属于本发明技术方案保护的范围。
Claims (8)
1.一种基于大面积致密储层地震岩石物理反演方法,其特征在于,包括如下具体步骤:
步骤101,基于孔隙弹性理论预测波响应频散,建立多尺度岩石物理模型的参量因素以联系多尺度数据;所述建立多尺度岩石物理模型的参量因素,是以岩石的矿物成分、孔隙结构及所处的地层环境会对岩石的波响应特征产生影响为依据,确定包括温度和压力为储层环境因素、包括矿物组分、孔隙形状、泥质含量和孔隙结构为储层岩性因素、包括流体黏性和气水斑状饱和为储层流体因素;
步骤102,基于模型及部分井试气情况,对测井解释结果进行分析和校正,开展声波测井与超声波两种尺度下岩石物理参数的流体敏感性分析,筛选出多个观测尺度下对孔隙度、含气饱和度最为敏感的岩石物理参数;所述超声波和声波测井两种尺度下岩石物理参数是指弹性参数及弹性参数的组合;所述弹性参数至少包括以下物理量:纵波速度Vp、横波速度Vs、纵波阻抗Zp、横波阻抗Zs、纵横波速度比Vp/Vs、拉梅常数λ、剪切模量μ、拉梅常数和密度乘积λρ、拉梅常数和剪切模量乘积λμ、准压力PR、剪切模量和密度乘积μρ;所述流体敏感性分析包括:超声波尺度下测量纵波速度Vp和横波速度Vs、以及气和水的变饱和度的情况下的波速,即Vp/Vs与波阻抗交会图;
步骤103,基于各参考井资料构建的单井岩石物理模板,优选各单井模板制作工区标准模板;所述工区标准模板是指以优选的敏感参数λρ为纵坐标,以纵波阻抗为横坐标;
步骤104,基于储层地质特征的横向变化与非均质性,在大区内根据所有井试气情况在各井坐标处对岩石物理模板的输入参数进行微调,并在整个工区内进行优化,建立三维工区岩石物理模板数据体,并将三维工区岩石物理模板数据体与地震叠前反演结合计算目的层孔隙度、饱和度;在整个工区开展大面积的三维岩石物理模板参数反演,并对反演结果进行平滑处理,最终输出储层参数反演数据体,进而实现储层孔隙度及饱和度的定量解释;所述建立三维工区岩石物理模板数据体是指将待反演、解释的三维地震数据体根据工程需要进行切割、整理,开展叠前三维地震反演,并进行储层孔隙度、饱和度的反演计算。
2.根据权利要求1所述的一种基于大面积致密储层地震岩石物理反演方法,其特征在于,步骤101所述建立多尺度岩石物理模型的建模,包括计算岩石基质弹性模量和岩石骨架弹性模量,采用Voigt-Reuss-Hill平均方程求取基质的等效弹性模量:
其中MVRH是矿物基质的弹性模量,fi和Mi分别是第i种组分的体积含量和弹性模量,N是矿物组分的总数量;采用微分等效介质(DEM)理论(Mavko,1998),求取白云岩的干岩石骨架的体积模量与剪切模量
(1-y)d/dy[K*(y)]=(K2-K*(y))P(*2)(y) (2a)
(1-y)d/dy[μ*(y)]=(μ2-μ*(y))Q(*2)(y) (2b)
初始条件是K*(0)=K1和μ*(0)=μ1,其中K1,μ1是初始主相材料的体积模量和剪切模量(相1),K2,μ2是逐渐加入的包含物的体积模量和剪切模量(相2),y为相2的含量,P(*2)和Q(*2)与嵌入体形状有关。
3.根据权利要求1所述的一种基于大面积致密储层地震岩石物理反演方法,其特征在于,步骤101所述建立多尺度岩石物理模型的建模还包括考虑储层环境下的岩石,采用vander Waals方程,估算储层流体在不同温度、压力条件下的天然气的密度和体积模量。
4.根据权利要求1所述的一种基于大面积致密储层地震岩石物理反演方法,其特征在于,步骤101所述建立多尺度岩石物理模型的建模还包括考虑孔隙流体的非均质分布,忽略孔隙结构的非均质性,利用Biot-Rayleigh方程来实现饱和流体岩石纵横波速度的预测,同时还考虑弹性波速度频散的影响,以实现不同尺度多学科数据间的融合,Biot-Rayleigh方程具体形式如下:
其中u=[u1,u2,u3],分别代表三种组分的空间矢量位移(骨架、流体1、流体2),下标1、2、3表示矢量空间的三个方向;表示地震波激励过程中产生的局域流体变形增量,eij、和分别为:
其中:x1,x2与x3分别表示三个方向的坐标。φ1和φ2表示两类孔隙的绝对孔隙度,岩石的总孔隙度φ=φ1+φ2;φ10与φ20分别表示两个区域内的局部孔隙度,如岩石内部仅含有一种骨架,但却饱和有两种流体,则φ10=φ20=φ;假设φ1代表水孔(背景相流体),φ2代表气孔(包含物相流体),那么φ1/φ就是含水饱和度,φ2/φ就是含气饱和度;ρf1和η1表示背景相流体的密度与黏度,ρf2和η2表示包含物相流体的密度与黏度;R0表示气泡半径,κ10表示岩石渗透率;弹性参数A、N、Q1、R1、Q2与R2,密度参数ρ11、ρ12、ρ13、ρ22与ρ33,耗散参数b1与b2的确定方法。
5.根据权利要求1所述的一种基于大面积致密储层地震岩石物理反演方法,其特征在于,步骤102还包括对比分析超声波尺度和测井尺度下岩石物理参数敏感性规律、次序基本一致,对孔隙流体最为敏感性参数为λ和λρ。
6.根据权利要求1所述的一种基于大面积致密储层地震岩石物理反演方法,其特征在于,步骤103对所述工区标准岩石物理模板格点的横、纵坐标进行校正,是指保证工区标准岩石物理模板的描述结果基本涵盖标准井所有数据,则校正后的工区标准岩石物理模板表示为:
假设各单井岩石物理模板在二维交会图中的任意一个模板格点位置可根据其横、纵坐标用(Mk(i,j),Nk(i,j))表示,其中i、j分别对应于孔隙度和饱和度的渐变(i=1,2…11对应于孔隙度0.02渐变到0.12;j=1,2…11对应于饱和度0渐变到100%),k表示第k个参考井,则标准模板在该格点位置的值(Ms(i,j),Ns(i,j))可表示为:
其中:A(k)表示第k口井的权重,B(i,j)、C(i,j)分别表示对应于i、j的模板格点上基于综合各井参考数据对工区标准岩石物理模板格点的横、纵坐标进行的校正。
7.根据权利要求1所述的一种基于大面积致密储层地震岩石物理反演方法,其特征在于,步骤104所述储层参数岩石物理反演,是指基于工区标准岩石物理模板,在各井位置开展地震反演与流体检测测试,具体方法为:
从三维地震数据体中抽取二维过井线,对井坐标附件的目的层进行储层、流体参数估算,对比估算结果与已知钻井、试气结果,适度调整岩石物理模板,试反演结果、解释结论与实际试气情况达到最佳吻合;在第k口井附近,工区标准岩石物理模板经二维地震反演调试后的地震反演模板可用(M′k(i,j),N′k(i,j))表示;
在工区内,基于各井的地震反演模板(M′k(i,j),N′k(i,j)),开展全区优化,生成工区岩石物理模型的三维数据体,具体对应于工区内每一个坐标位置(x,y),都存在一个独立的地震岩石物理模型,其模板(M3D(x,y,i,j),N3D(x,y,i,j))的确定依赖于各井位置的地震反演模板:
M3D(x,y,i,j)=M′k(i,j),若x=xk,y=yk (5a)
N3D(x,y,i,j)=N′k(i,j),若x=xk,y=yk (5b)
其中:(xk,yk)为第k口井坐标,L为总井数,Q(x,y,k)为计算工区内坐标(x,y)处模板时所采用的第k口井的权重系数,可由下式确定:
基于以上方法,即可兼顾工区内各井资料,基于各地理位置的测井观测控制储层参数的地震反演及解释,空间位置越靠近某口参考井,则受该口井的控制因素影响越为显著;(M3D(x,y,i,j),N3D(x,y,i,j))即工区岩石物理模型三维数据体,制作过程中按地震数据.sgy标准格式首先沿xline方向,继而沿inline方向逐条xline线进行处理。
8.根据权利要求1所述的一种基于大面积致密储层地震岩石物理反演方法,其特征在于,步骤104所述采用加权平均法对反演的孔隙度、饱和度进行平滑处理,以减弱反演数据的野值、奇异值影响;
假定离目标点越近反演结果对目标点造成的影响越大,根据二维正态分布定义了三种加权模板,如式(7)所示;T1与T2是通过二维高斯离散模板演变而成,而T3就是3×3邻域的二维高斯模板,模板中最大权值为目标点所在位置,可根据最小均匀度所在邻域位置,将模板进行适当调整,以形成与该邻域相关的模板:
以T1为例,在邻域内对二维反演剖面data(i,j)进行平滑处理,具体算法可表示为:
其中dataS即为平滑处理后目标点的数值。
Priority Applications (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610563338.7A CN106054248B (zh) | 2016-07-15 | 2016-07-15 | 一种基于大面积致密储层地震岩石物理反演方法 |
EP17826964.3A EP3474045A4 (en) | 2016-07-15 | 2017-07-11 | SEISMIC STONE-PHYSICAL INVERSION PROCESS BASED ON A LARGE DENSITY RESERVOIR |
US16/317,822 US10983232B2 (en) | 2016-07-15 | 2017-07-11 | Seismic rock physics inversion method based on large area tight reservoir |
PCT/CN2017/092445 WO2018010628A1 (zh) | 2016-07-15 | 2017-07-11 | 一种基于大面积致密储层地震岩石物理反演方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610563338.7A CN106054248B (zh) | 2016-07-15 | 2016-07-15 | 一种基于大面积致密储层地震岩石物理反演方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106054248A true CN106054248A (zh) | 2016-10-26 |
CN106054248B CN106054248B (zh) | 2017-07-18 |
Family
ID=57187118
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610563338.7A Active CN106054248B (zh) | 2016-07-15 | 2016-07-15 | 一种基于大面积致密储层地震岩石物理反演方法 |
Country Status (4)
Country | Link |
---|---|
US (1) | US10983232B2 (zh) |
EP (1) | EP3474045A4 (zh) |
CN (1) | CN106054248B (zh) |
WO (1) | WO2018010628A1 (zh) |
Cited By (29)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106646668A (zh) * | 2016-11-30 | 2017-05-10 | 中石化石油工程技术服务有限公司 | 一种雷达测井标准井模型的建立方法 |
CN107367762A (zh) * | 2017-06-27 | 2017-11-21 | 中国石油天然气股份有限公司 | 一种确定储层参数的方法及装置 |
WO2018010628A1 (zh) * | 2016-07-15 | 2018-01-18 | 河海大学 | 一种基于大面积致密储层地震岩石物理反演方法 |
CN108254780A (zh) * | 2018-01-22 | 2018-07-06 | 河海大学 | 一种微地震定位及各向异性速度结构层析成像方法 |
CN108829717A (zh) * | 2018-05-07 | 2018-11-16 | 西南石油大学 | 一种基于地震数据进行深水水道构型量化分析和形态模拟的数据库系统及方法 |
CN108957534A (zh) * | 2018-06-12 | 2018-12-07 | 中国石油天然气股份有限公司 | 含气饱和度预测方法及装置 |
CN109100796A (zh) * | 2018-06-13 | 2018-12-28 | 中国石油天然气集团有限公司 | 一种天然气水合物地震数据处理方法及装置 |
CN109115987A (zh) * | 2018-07-20 | 2019-01-01 | 中国石油天然气股份有限公司 | 一种基于岩石物理模型的流体因子的评价方法及装置 |
CN109374497A (zh) * | 2018-10-29 | 2019-02-22 | 河海大学 | 一种岩石微观孔隙结构测试方法 |
CN109580679A (zh) * | 2017-09-29 | 2019-04-05 | 中国石油化工股份有限公司 | 基于岩石物理分析的粘土矿物类型及含量识别方法和系统 |
CN110275206A (zh) * | 2019-08-12 | 2019-09-24 | 河海大学 | 一种裂隙-孔隙型岩石物理弹性模板 |
CN110703330A (zh) * | 2018-07-10 | 2020-01-17 | 中国石油化工股份有限公司 | 基于岩石物理的碳酸盐岩测井校正和横波速度预测方法 |
CN110887772A (zh) * | 2018-09-11 | 2020-03-17 | 中国石油化工股份有限公司 | 一种碳酸盐岩储层渗透率识别方法、系统及装置 |
CN110954960A (zh) * | 2019-11-20 | 2020-04-03 | 中国石油大学(华东) | 一种海底冷泉地震响应数值模拟方法及系统 |
CN111158052A (zh) * | 2020-01-07 | 2020-05-15 | 吉林大学 | 基于三维储层反演井中含铀砂体井旁分布范围的圈定方法 |
CN111352172A (zh) * | 2020-02-21 | 2020-06-30 | 吉林大学 | 一种用井震联合法获取铀异常在砂体中空间分布位置的方法 |
CN111751900A (zh) * | 2019-03-29 | 2020-10-09 | 中国石油天然气股份有限公司 | 一种流体识别方法和装置 |
CN111852467A (zh) * | 2020-07-28 | 2020-10-30 | 核工业北京地质研究院 | 一种砂岩铀矿矿体延伸范围的圈定方法及系统 |
CN111897011A (zh) * | 2020-07-27 | 2020-11-06 | 中国石油大学(华东) | 储层孔隙特征确定方法、装置及设备 |
CN112285801A (zh) * | 2019-07-24 | 2021-01-29 | 中国石油天然气股份有限公司 | 碳酸盐岩储层渗透率的预测方法、设备及系统 |
CN112649855A (zh) * | 2019-10-11 | 2021-04-13 | 中国石油化工股份有限公司 | 三维含气饱和度预测方法及系统 |
CN112684503A (zh) * | 2020-11-16 | 2021-04-20 | 中国石油天然气股份有限公司 | 一种基于叠后地震属性的储层流体识别方法及装置 |
CN112946739A (zh) * | 2021-01-27 | 2021-06-11 | 中国石油天然气股份有限公司 | 裂隙-溶蚀孔双孔系统中深层碳酸盐岩储层地震岩石物理模板构建方法及储层参数预测方法 |
CN113049471A (zh) * | 2021-03-23 | 2021-06-29 | 中国石油大学(北京) | 一种碳酸盐岩层序地层的孔隙度演化过程的恢复方法 |
CN113311483A (zh) * | 2020-02-27 | 2021-08-27 | 中国石油天然气集团有限公司 | 基于页岩油的叠前弹性参数组合反演方法及系统 |
CN113589385A (zh) * | 2021-08-11 | 2021-11-02 | 成都理工大学 | 一种基于地震散射波场分析的储层特征反演方法 |
CN114492235A (zh) * | 2022-01-14 | 2022-05-13 | 成都理工大学 | 一种适用于碳酸盐岩地层的地层孔隙流体压力预测方法 |
CN115201910A (zh) * | 2022-07-07 | 2022-10-18 | 中国石油化工股份有限公司 | 一种岩石物理动态模板建立方法、系统、存储介质及设备 |
CN115586573A (zh) * | 2022-09-15 | 2023-01-10 | 河海大学 | 一种致密砂岩储层的动态约束物性参数地震反演方法 |
Families Citing this family (67)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109975876A (zh) * | 2019-03-20 | 2019-07-05 | 中国石油化工股份有限公司 | 一种基于构造层位的井震融合速度模型的建模方法 |
CN110276091B (zh) * | 2019-04-26 | 2022-10-14 | 中国石油化工股份有限公司 | 基于岩石多重孔隙结构的弹性波响应模型建模方法 |
CN112147696B (zh) * | 2019-06-27 | 2023-05-23 | 中国石油化工股份有限公司 | 基于储集性能指数的页岩储集性能地震评价方法及系统 |
CN112415588A (zh) * | 2019-08-22 | 2021-02-26 | 中国石油化工股份有限公司 | 基于三维地震网格的储层参数油气储量计算方法及系统 |
US20230367031A1 (en) * | 2019-08-26 | 2023-11-16 | Landmark Graphics Corporation | Building scalable geological property models using machine learning algorithms |
CN112485827A (zh) * | 2019-09-12 | 2021-03-12 | 中国石油天然气股份有限公司 | 裂缝孔隙型储层的岩石物理建模方法及装置 |
CN112649870B (zh) * | 2019-10-12 | 2024-08-16 | 中国石油化工股份有限公司 | 岩石物理建模中矿物弹性参数的确定方法及系统 |
CN112780253B (zh) * | 2020-01-20 | 2022-05-10 | 中国石油天然气集团有限公司 | 裂缝型储层的预测及评价方法 |
CN113267809B (zh) * | 2020-02-17 | 2024-05-28 | 中国石油天然气集团有限公司 | I类页岩储层预测方法及装置 |
CN111736233A (zh) * | 2020-05-11 | 2020-10-02 | 中国石油天然气集团有限公司 | 蒸汽腔分布的预测方法及装置 |
CN111751882B (zh) * | 2020-05-12 | 2022-12-02 | 中国石油天然气股份有限公司 | 生成敏感岩石物理参数的方法和装置 |
CN113655523B (zh) * | 2020-05-12 | 2024-02-06 | 中国石油化工股份有限公司 | 一种定量预测及评价深层碳酸盐岩局部盖层的地球物理方法与系统 |
CN113687412B (zh) * | 2020-05-18 | 2024-03-26 | 中国石油化工股份有限公司 | 盐间地层压力预测方法、装置、电子设备及介质 |
CN111880234B (zh) * | 2020-06-24 | 2022-09-30 | 西南石油大学 | 基于常规测井对致密储层裂缝的精细识别方法 |
CN113900142B (zh) * | 2020-07-06 | 2024-06-07 | 中国石油化工股份有限公司 | 横波速度预测方法、存储介质和计算机设备 |
CN113917531B (zh) * | 2020-07-08 | 2024-02-23 | 中国石油化工股份有限公司 | 基于扰动体积模量的压力系数预测方法及电子设备 |
CN113933910B (zh) * | 2020-07-13 | 2024-04-09 | 中国石油化工股份有限公司 | 物性参数同时计算的方法、系统、存储介质与电子设备 |
CN113945970B (zh) * | 2020-07-15 | 2024-06-07 | 中国石油化工股份有限公司 | 一种致密砂岩储层预测方法 |
CN114114406B (zh) * | 2020-08-26 | 2023-09-01 | 中国石油化工股份有限公司 | 一种储层渗透率估计方法及装置 |
CN112133376A (zh) * | 2020-08-27 | 2020-12-25 | 中国石油天然气集团有限公司 | 一种饱含水地层体积模量的流体识别方法 |
CN114114409B (zh) * | 2020-08-27 | 2024-08-30 | 中国石油化工股份有限公司 | 海域天然气水合物的岩石物理建模方法、电子设备及介质 |
CN112213778B (zh) * | 2020-08-27 | 2024-01-30 | 中国石油天然气股份有限公司 | 碳酸盐岩储层孔隙流体饱和度识别方法及量版建立方法 |
CN112014881B (zh) * | 2020-08-27 | 2023-10-27 | 中海石油(中国)有限公司 | 基于时移地震的水驱速度预测方法 |
CN114114454B (zh) * | 2020-08-31 | 2023-06-23 | 中国石油化工股份有限公司 | 考虑压实作用的岩石物理建模及纵波速度估算方法 |
CN114152982A (zh) * | 2020-09-08 | 2022-03-08 | 中国石油天然气股份有限公司 | 叠前储层分析方法及装置 |
CN114428372B (zh) * | 2020-09-09 | 2023-06-27 | 中国石油化工股份有限公司 | 一种自适应的岩石物理建模方法 |
CN112394416B (zh) * | 2020-10-20 | 2024-06-25 | 中国石油天然气集团有限公司 | 非均质断裂控制储层预测方法及装置 |
CN112347670B (zh) * | 2020-10-26 | 2022-09-27 | 青海大学 | 一种基于神经网络响应面的堆石料蠕变参数预测方法 |
CN112433248B (zh) * | 2020-11-10 | 2022-09-30 | 成都理工大学 | 一种碳酸盐岩沉积环境下隐蔽储层的探测方法 |
CN112392477B (zh) * | 2020-12-04 | 2023-08-22 | 中国石油天然气股份有限公司 | 一种单井潜力快速预测方法 |
CN112505772B (zh) * | 2020-12-10 | 2022-05-31 | 中国石油大学(华东) | 一种利用孔隙、裂隙介质弹性波理论反演岩石孔隙分布特征的方法 |
CN112651171B (zh) * | 2020-12-15 | 2024-05-14 | 清华大学 | 一种数据驱动岩石物理建模方法及系统 |
CN112835098B (zh) * | 2021-01-05 | 2024-03-01 | 中国石油天然气股份有限公司 | 风化壳岩溶储层的储能系数预测方法及装置 |
CN112882100B (zh) * | 2021-02-25 | 2024-01-12 | 中海石油深海开发有限公司 | 一种储层参数确定方法、装置、电子设备和存储介质 |
CN113075728B (zh) * | 2021-02-26 | 2021-11-09 | 河海大学 | 一种建立致密砂岩多尺度三维岩石物理图板的方法 |
CN112949065B (zh) * | 2021-03-04 | 2022-05-27 | 长江水利委员会长江科学院 | 模拟层状岩体力学行为的双尺度方法、装置、存储介质及设备 |
CN112987126A (zh) * | 2021-03-12 | 2021-06-18 | 中国石油大学(华东) | 一种岩石物理相的地球物理预测方法 |
CN113009565A (zh) * | 2021-03-24 | 2021-06-22 | 中国石油大学(北京) | 基于sca模型的地震波速度参数确定方法、装置及设备 |
CN113156502B (zh) * | 2021-03-31 | 2022-08-23 | 华能煤炭技术研究有限公司 | 一种基于地震资料的煤层流体压力预测方法 |
CN113219527A (zh) * | 2021-04-01 | 2021-08-06 | 中国石油化工股份有限公司 | 一种基于导航金字塔分解的油气储层反演方法及装置 |
CN113514891B (zh) * | 2021-04-27 | 2022-10-14 | 东北石油大学 | 一种基于地震反演的油水过渡带范围确定方法 |
CN113376699B (zh) * | 2021-06-08 | 2022-04-29 | 西南石油大学 | 一种基于单声波老测井系列的多矿物组分计算方法 |
CN113552627B (zh) * | 2021-07-01 | 2023-01-17 | 中国石油化工股份有限公司 | 一种自适应地震横波速度预测方法及装置 |
CN113552626B (zh) * | 2021-07-01 | 2023-01-17 | 中国石油化工股份有限公司 | 一种致密砂岩储层的地震横波速度预测方法及装置 |
CN113687419B (zh) * | 2021-08-24 | 2023-07-07 | 太原理工大学 | 基于主控要素地震-地质一体化分析构造煤分布预测方法 |
CN113740911B (zh) * | 2021-09-06 | 2023-09-26 | 北京海润联创石油科技有限公司 | 一种基于坐标旋转波阻抗反演提高储层预测精度的方法 |
CN113640883B (zh) * | 2021-09-16 | 2022-09-23 | 中国石油大学(北京) | 测井速度预测、基于井震融合的孔隙度确定方法及装置 |
CN113671565B (zh) * | 2021-09-17 | 2022-11-04 | 中国海洋石油集团有限公司 | 一种针对巨厚储层气藏开发的地震多尺度储层预测方法 |
CN113887046B (zh) * | 2021-09-30 | 2023-03-28 | 重庆科华安全设备有限责任公司 | 一种基于三维地质体的煤矿巷道建模方法 |
CN113960659B (zh) * | 2021-10-14 | 2023-02-03 | 中国矿业大学 | 一种地震岩石物理驱动的煤层气储层含气量预测方法 |
CN114325845B (zh) * | 2021-10-26 | 2024-03-15 | 重庆科技学院 | 一种基于数字岩心技术的非常规储层多尺度融合方法 |
US11746644B2 (en) * | 2021-12-02 | 2023-09-05 | Halliburton Energy Services, Inc. | Measuring low-frequency casing guided waves to evaluate cement bond condition behind casing in the presence of a tubing |
CN113989454B (zh) * | 2021-12-30 | 2022-03-18 | 山东省地质测绘院 | 适用于地质数据和地理信息数据的融合方法、装置及系统 |
CN114674934B (zh) * | 2022-02-18 | 2023-06-13 | 河海大学 | 一种建立饱和重油岩石波速随温度变化的理论模型的方法 |
CN115079261B (zh) * | 2022-06-06 | 2023-07-11 | 吉林大学 | 一种基于多参数定量解释模板的致密砂岩气储层评价方法 |
CN115061202B (zh) * | 2022-06-10 | 2024-03-15 | 吉林大学 | 一种页岩含气性地震储层直接检测方法 |
CN115308796B (zh) * | 2022-06-23 | 2023-07-07 | 吉林大学 | 基于各向异性频散属性的页岩含气性与层理结构预测方法 |
CN116840912B (zh) * | 2022-09-13 | 2024-07-19 | 上海电子信息职业技术学院 | 不同裂隙纵横比部分饱和孔隙介质中纵波速度预测方法 |
CN115392090B (zh) * | 2022-09-13 | 2023-05-23 | 中国矿业大学 | 基于三维数字岩心的地震波频散和衰减特征的预测方法 |
CN115508890B (zh) * | 2022-09-28 | 2023-05-12 | 北京中恒利华石油技术研究所 | 一种裂缝孔隙型储层叠前叠后反演方法 |
CN115857047B (zh) * | 2022-09-28 | 2023-07-21 | 北京中恒利华石油技术研究所 | 一种地震储层综合预测方法 |
CN115961950B (zh) * | 2023-02-03 | 2024-06-25 | 西南石油大学 | 基于岩石声学信息的水岩作用下钻井液密度下限预测方法 |
CN116738517B (zh) * | 2023-07-04 | 2024-07-26 | 中海石油(中国)有限公司深圳分公司 | 一种储层预测方法、装置、电子设备及介质 |
CN116522695B (zh) * | 2023-07-05 | 2023-09-05 | 矿冶科技集团有限公司 | 岩体裂纹生长轨迹构建方法、装置、设备及存储介质 |
CN116757557B (zh) * | 2023-08-15 | 2023-11-07 | 山东新巨龙能源有限责任公司 | 一种基于数据分析的原矸充填开采质量评估方法 |
CN117388919B (zh) * | 2023-08-29 | 2024-05-28 | 河海大学 | 一种基于dnn预测致密油储层横波速度的方法 |
CN117250658B (zh) * | 2023-11-17 | 2024-02-09 | 核工业北京地质研究院 | 建立研究区的地震数据集的方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5059907A (en) * | 1988-09-06 | 1991-10-22 | Amoco Corporation | Method for evaluating the water saturation of subterranean formations from dielectric permittivity measurements |
CN102854531A (zh) * | 2012-09-11 | 2013-01-02 | 中国石油天然气股份有限公司 | 用于储层烃类检测的多尺度岩石物理图版方法及装置 |
CN103163553A (zh) * | 2013-02-01 | 2013-06-19 | 中国石油天然气股份有限公司 | 基于多重孔隙介质模型的地震烃类检测方法及装置 |
CN103645516A (zh) * | 2013-11-20 | 2014-03-19 | 中国石油大学(北京) | 基于岩石物理相控油气作用确定油气产能的方法 |
CN103760081A (zh) * | 2013-12-31 | 2014-04-30 | 中国石油天然气股份有限公司 | 基于孔隙结构特征的碳酸盐岩储层的气藏预测方法及系统 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB2400438B (en) * | 2003-04-11 | 2005-06-01 | Westerngeco Ltd | Determination of waveguide parameters |
CN102749643B (zh) | 2011-04-22 | 2015-06-03 | 中国石油天然气股份有限公司 | 一种面波地震记录的频散响应获取方法及其装置 |
US10422922B2 (en) * | 2012-05-24 | 2019-09-24 | Exxonmobil Upstream Research Company | Method for predicting rock strength by inverting petrophysical properties |
CN103984010B (zh) | 2014-04-16 | 2016-02-17 | 岳航羽 | 一种基于三项频率依赖avo反演的流体识别方法 |
CN106054248B (zh) * | 2016-07-15 | 2017-07-18 | 河海大学 | 一种基于大面积致密储层地震岩石物理反演方法 |
-
2016
- 2016-07-15 CN CN201610563338.7A patent/CN106054248B/zh active Active
-
2017
- 2017-07-11 EP EP17826964.3A patent/EP3474045A4/en not_active Withdrawn
- 2017-07-11 WO PCT/CN2017/092445 patent/WO2018010628A1/zh unknown
- 2017-07-11 US US16/317,822 patent/US10983232B2/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5059907A (en) * | 1988-09-06 | 1991-10-22 | Amoco Corporation | Method for evaluating the water saturation of subterranean formations from dielectric permittivity measurements |
CN102854531A (zh) * | 2012-09-11 | 2013-01-02 | 中国石油天然气股份有限公司 | 用于储层烃类检测的多尺度岩石物理图版方法及装置 |
CN103163553A (zh) * | 2013-02-01 | 2013-06-19 | 中国石油天然气股份有限公司 | 基于多重孔隙介质模型的地震烃类检测方法及装置 |
CN103645516A (zh) * | 2013-11-20 | 2014-03-19 | 中国石油大学(北京) | 基于岩石物理相控油气作用确定油气产能的方法 |
CN103760081A (zh) * | 2013-12-31 | 2014-04-30 | 中国石油天然气股份有限公司 | 基于孔隙结构特征的碳酸盐岩储层的气藏预测方法及系统 |
Non-Patent Citations (1)
Title |
---|
巴晶: "非均质天然气藏的岩石物理模型及含气饱和度反演", 《地球物理学报》 * |
Cited By (46)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2018010628A1 (zh) * | 2016-07-15 | 2018-01-18 | 河海大学 | 一种基于大面积致密储层地震岩石物理反演方法 |
US10983232B2 (en) | 2016-07-15 | 2021-04-20 | Hohai University | Seismic rock physics inversion method based on large area tight reservoir |
CN106646668A (zh) * | 2016-11-30 | 2017-05-10 | 中石化石油工程技术服务有限公司 | 一种雷达测井标准井模型的建立方法 |
CN107367762A (zh) * | 2017-06-27 | 2017-11-21 | 中国石油天然气股份有限公司 | 一种确定储层参数的方法及装置 |
CN107367762B (zh) * | 2017-06-27 | 2019-06-11 | 中国石油天然气股份有限公司 | 一种确定储层参数的方法及装置 |
CN109580679A (zh) * | 2017-09-29 | 2019-04-05 | 中国石油化工股份有限公司 | 基于岩石物理分析的粘土矿物类型及含量识别方法和系统 |
CN108254780A (zh) * | 2018-01-22 | 2018-07-06 | 河海大学 | 一种微地震定位及各向异性速度结构层析成像方法 |
CN108829717A (zh) * | 2018-05-07 | 2018-11-16 | 西南石油大学 | 一种基于地震数据进行深水水道构型量化分析和形态模拟的数据库系统及方法 |
CN108829717B (zh) * | 2018-05-07 | 2021-10-08 | 西南石油大学 | 一种基于地震数据进行深水水道构型量化分析和形态模拟的数据库系统及方法 |
CN108957534A (zh) * | 2018-06-12 | 2018-12-07 | 中国石油天然气股份有限公司 | 含气饱和度预测方法及装置 |
CN108957534B (zh) * | 2018-06-12 | 2020-02-14 | 中国石油天然气股份有限公司 | 含气饱和度预测方法及装置 |
CN109100796A (zh) * | 2018-06-13 | 2018-12-28 | 中国石油天然气集团有限公司 | 一种天然气水合物地震数据处理方法及装置 |
CN110703330A (zh) * | 2018-07-10 | 2020-01-17 | 中国石油化工股份有限公司 | 基于岩石物理的碳酸盐岩测井校正和横波速度预测方法 |
CN110703330B (zh) * | 2018-07-10 | 2021-12-24 | 中国石油化工股份有限公司 | 基于岩石物理的碳酸盐岩测井校正和横波速度预测方法 |
CN109115987A (zh) * | 2018-07-20 | 2019-01-01 | 中国石油天然气股份有限公司 | 一种基于岩石物理模型的流体因子的评价方法及装置 |
CN110887772A (zh) * | 2018-09-11 | 2020-03-17 | 中国石油化工股份有限公司 | 一种碳酸盐岩储层渗透率识别方法、系统及装置 |
CN109374497A (zh) * | 2018-10-29 | 2019-02-22 | 河海大学 | 一种岩石微观孔隙结构测试方法 |
CN111751900A (zh) * | 2019-03-29 | 2020-10-09 | 中国石油天然气股份有限公司 | 一种流体识别方法和装置 |
CN112285801B (zh) * | 2019-07-24 | 2023-02-28 | 中国石油天然气股份有限公司 | 碳酸盐岩储层渗透率的预测方法、设备及系统 |
CN112285801A (zh) * | 2019-07-24 | 2021-01-29 | 中国石油天然气股份有限公司 | 碳酸盐岩储层渗透率的预测方法、设备及系统 |
CN110275206A (zh) * | 2019-08-12 | 2019-09-24 | 河海大学 | 一种裂隙-孔隙型岩石物理弹性模板 |
CN110275206B (zh) * | 2019-08-12 | 2021-09-28 | 河海大学 | 一种裂隙-孔隙型岩石物理弹性模板 |
CN112649855A (zh) * | 2019-10-11 | 2021-04-13 | 中国石油化工股份有限公司 | 三维含气饱和度预测方法及系统 |
CN112649855B (zh) * | 2019-10-11 | 2024-04-09 | 中国石油化工股份有限公司 | 三维含气饱和度预测方法及系统 |
CN110954960A (zh) * | 2019-11-20 | 2020-04-03 | 中国石油大学(华东) | 一种海底冷泉地震响应数值模拟方法及系统 |
CN111158052A (zh) * | 2020-01-07 | 2020-05-15 | 吉林大学 | 基于三维储层反演井中含铀砂体井旁分布范围的圈定方法 |
CN111158052B (zh) * | 2020-01-07 | 2021-07-13 | 吉林大学 | 基于三维储层反演井中含铀砂体井旁分布范围的圈定方法 |
CN111352172A (zh) * | 2020-02-21 | 2020-06-30 | 吉林大学 | 一种用井震联合法获取铀异常在砂体中空间分布位置的方法 |
CN111352172B (zh) * | 2020-02-21 | 2021-12-03 | 吉林大学 | 一种用井震联合法获取铀异常在砂体中空间分布位置的方法 |
CN113311483A (zh) * | 2020-02-27 | 2021-08-27 | 中国石油天然气集团有限公司 | 基于页岩油的叠前弹性参数组合反演方法及系统 |
CN113311483B (zh) * | 2020-02-27 | 2024-04-30 | 中国石油天然气集团有限公司 | 基于页岩油的叠前弹性参数组合反演方法及系统 |
CN111897011B (zh) * | 2020-07-27 | 2024-03-22 | 中国石油大学(华东) | 储层孔隙特征确定方法、装置及设备 |
CN111897011A (zh) * | 2020-07-27 | 2020-11-06 | 中国石油大学(华东) | 储层孔隙特征确定方法、装置及设备 |
CN111852467B (zh) * | 2020-07-28 | 2023-04-07 | 核工业北京地质研究院 | 一种砂岩铀矿矿体延伸范围的圈定方法及系统 |
CN111852467A (zh) * | 2020-07-28 | 2020-10-30 | 核工业北京地质研究院 | 一种砂岩铀矿矿体延伸范围的圈定方法及系统 |
CN112684503A (zh) * | 2020-11-16 | 2021-04-20 | 中国石油天然气股份有限公司 | 一种基于叠后地震属性的储层流体识别方法及装置 |
CN112946739A (zh) * | 2021-01-27 | 2021-06-11 | 中国石油天然气股份有限公司 | 裂隙-溶蚀孔双孔系统中深层碳酸盐岩储层地震岩石物理模板构建方法及储层参数预测方法 |
US11487045B2 (en) | 2021-03-23 | 2022-11-01 | China University Of Petroleum-Beijing | Method for recovering porosity evolution process of sequence stratigraphy of carbonate rocks |
CN113049471A (zh) * | 2021-03-23 | 2021-06-29 | 中国石油大学(北京) | 一种碳酸盐岩层序地层的孔隙度演化过程的恢复方法 |
CN113589385A (zh) * | 2021-08-11 | 2021-11-02 | 成都理工大学 | 一种基于地震散射波场分析的储层特征反演方法 |
CN113589385B (zh) * | 2021-08-11 | 2023-08-04 | 成都理工大学 | 一种基于地震散射波场分析的储层特征反演方法 |
CN114492235B (zh) * | 2022-01-14 | 2023-04-07 | 成都理工大学 | 一种适用于碳酸盐岩地层的地层孔隙流体压力预测方法 |
CN114492235A (zh) * | 2022-01-14 | 2022-05-13 | 成都理工大学 | 一种适用于碳酸盐岩地层的地层孔隙流体压力预测方法 |
CN115201910B (zh) * | 2022-07-07 | 2023-05-26 | 中国石油化工股份有限公司 | 一种岩石物理动态模板建立方法、系统、存储介质及设备 |
CN115201910A (zh) * | 2022-07-07 | 2022-10-18 | 中国石油化工股份有限公司 | 一种岩石物理动态模板建立方法、系统、存储介质及设备 |
CN115586573A (zh) * | 2022-09-15 | 2023-01-10 | 河海大学 | 一种致密砂岩储层的动态约束物性参数地震反演方法 |
Also Published As
Publication number | Publication date |
---|---|
US20200132869A1 (en) | 2020-04-30 |
EP3474045A1 (en) | 2019-04-24 |
EP3474045A4 (en) | 2020-01-15 |
US10983232B2 (en) | 2021-04-20 |
WO2018010628A1 (zh) | 2018-01-18 |
CN106054248B (zh) | 2017-07-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106054248B (zh) | 一种基于大面积致密储层地震岩石物理反演方法 | |
Paillet et al. | Acoustic waves in boreholes | |
CN103256046B (zh) | 非常规油气藏水平井全缝长压裂参数模拟的方法及装置 | |
US8688425B2 (en) | Transport property data calculated from derivative seismic rock property data for transport modeling | |
CN105510993B (zh) | 前陆盆地深埋挤压型复杂膏盐岩层识别和分布预测方法 | |
CN103258091B (zh) | 非常规油气藏水平井段三维岩体力学模型建立的方法及装置 | |
CN102445708B (zh) | 三维等效富泥质砂岩速度预测模型 | |
CN109613612A (zh) | 一种碳酸盐岩颗粒滩精细刻画与预测方法 | |
CN105182424B (zh) | 一种基于斑块饱和模型定量预测储层孔隙度的方法和装置 | |
CN106855636A (zh) | 基于碳酸盐岩储层露头的原型地质模型地震正演方法 | |
GB2421822A (en) | Method for predicting quantitative values of a rock or fluid property in a reservoir using seismic data | |
CN103914620B (zh) | 一种计算断层破裂带裂缝张开空间分布的方法 | |
CN106368691A (zh) | 基于岩石物理地震信息三维异常孔隙压力预测方法 | |
CN105629308A (zh) | 相控非均质力学参数地应力方法 | |
Kianoush et al. | Compressional and shear interval velocity modeling to determine formation pressures in an oilfield of SW Iran | |
Montone et al. | P-wave velocity, density, and vertical stress magnitude along the crustal Po Plain (Northern Italy) from sonic log drilling data | |
CN104975851B (zh) | 用于振幅随炮检距变化道集分析的油藏模型优化方法 | |
CN107728205B (zh) | 一种地层压力预测方法 | |
CN103643949A (zh) | 一种储层含油气的定量预测方法及装置 | |
Bigi et al. | Discrete fracture network of the Latemar carbonate platform | |
Shahin et al. | Optimal dynamic rock-fluid physics template validated by petroelastic reservoir modeling | |
Rolfs | Integrated geomechanical, geophysical, and geochemical analysis of the Bakken Formation, Elm Coulee field, Williston Basin, Montana | |
Dou et al. | Reservoir Modeling of Braided River Reservoirs Based on Geological Knowledge Database: A Case Study of P1x Formation of the Daniudi Gas Field, Ordos Basin, China | |
Soleimani et al. | Petroleum reservoir simulation, Ramin Oil Field, Zagros, Iran | |
Han et al. | The method and application of numerical simulation of high-precision stress field and quantitative prediction of multiperiod fracture in carbonate reservoir |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |