CN105093331A - 获取岩石基质体积模量的方法 - Google Patents

获取岩石基质体积模量的方法 Download PDF

Info

Publication number
CN105093331A
CN105093331A CN201510410200.9A CN201510410200A CN105093331A CN 105093331 A CN105093331 A CN 105093331A CN 201510410200 A CN201510410200 A CN 201510410200A CN 105093331 A CN105093331 A CN 105093331A
Authority
CN
China
Prior art keywords
rock
depth point
bulk modulus
target
logging depth
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
CN201510410200.9A
Other languages
English (en)
Other versions
CN105093331B (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 National Petroleum Corp
BGP Inc
Original Assignee
Geophysical Prospecting Co of CNPC Chuanqing Drilling Engineering 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 Geophysical Prospecting Co of CNPC Chuanqing Drilling Engineering Co Ltd filed Critical Geophysical Prospecting Co of CNPC Chuanqing Drilling Engineering Co Ltd
Priority to CN201510410200.9A priority Critical patent/CN105093331B/zh
Publication of CN105093331A publication Critical patent/CN105093331A/zh
Application granted granted Critical
Publication of CN105093331B publication Critical patent/CN105093331B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

提供一种获取岩石基质体积模量的方法,包括:(A)利用目标测井深度点及其上下邻近的预定数量个测井深度点的测井资料确定目标测井深度点的等效岩石基质体积模量;(B)确定目标测井深度点的干岩石泊松比的变化区间和变化增量;(C)确定目标测井深度点的岩石基质体积模量的变化区间和变化增量;(D)在双收敛条件下利用自适应基质矿物模量反演方法反演得到目标测井深度点的最优自适应岩石基质体积模量;(E)将目标测井深度点的等效岩石基质体积模量和最优自适应岩石基质体积模量进行平均处理来获取目标测井深度点的岩石基质体积模量。根据所述方法,能够获取具有较高精度和稳定性的岩石基质体积模量,为后续预测油气藏奠定了有利基础。

Description

获取岩石基质体积模量的方法
技术领域
本发明总体说来涉及岩石物理学领域,更具体地讲,涉及一种获取岩石基质体积模量的方法。
背景技术
随着油气勘探程度的提高,我国陆上已进入了地层岩性油气藏勘探阶段,其中,岩石物理参数模型在油气藏勘探中起着重要的桥梁作用:岩石物理参数模型能够有效地将地震解释和油气藏特性参数联系起来,例如,通过岩石物理参数模型分析建立的基本准则,可定量或半定量地从地震数据中导出储层岩石、流体特性及油气藏等特征参数(例如,岩石总孔隙度、孔隙结构特征等),并可通过上述特征参数来预测油气藏。因此,岩石物理学研究是地震资料预测油气藏的物理基础。
通常,岩石物理参数模型建立的正确与否直接关系到相关后续数据的准确性。通常,岩石物理参数模型的应用可有下面几种情况:(1)用V-R-H模型(Hill,1952)和H-S模型(Hashin和Shtrikman,1963)计算岩石等效模量的上下界限;(2)用SC模型(Budiansky,1965)、模型(1974)和DEM模型(Berryman,1992)计算不同孔隙度及对应孔隙形状饱和岩石或干燥状态岩石等效模量;(3)用Hudson模型(Hudson,1980)和Eshelby-Cheng模型(1978)计算各向异性岩石等效弹性模量;(4)用Biot理论(1941)和Gassmann理论(1951)计算高低频极限时的纵波速度。其中,这些岩石物理参数模型的应用都需要对背景岩石进行描述,也就是说,需要确定背景岩石基质体积模量,如果基质体积模量取值不准确,将会极大地影响岩石物理参数模型的应用效果。
具体说来,岩石基质体积模量反映岩石矿物组分和胶结物等固体成分的性质,但是,由于基质类型复杂多变且受成岩作用、成岩后改造作用、地层压力和温度等综合影响容易造成同一基质模量值在不同地区差异较大,因此,直接求取岩石基质模量十分困难。
目前,求取岩石基质体积模量的方法主要有:
(1)基于实验室测试值的方法:该方法主要针对矿物类型单一的情况,具体说来,该方法主要参考国内外文献发表的岩石物理实验室测试值进行选取,也可通过在由所述测试值组成的取值范围内进行取值来获取。
(2)基于等效介质理论的方法:该方法主要针对矿物类型非单一的情况,具体说来,该方法需要已知组成岩石各成分的百分比含量和基质体积模量来获取岩石基质体积模量,但是,由于矿物各成分的百分比含量取值不准确,且受地区差异的影响,计算结果往往误差较大。
(3)测井资料统计法:该方法可根据测井资料通过选取岩石物性较均匀的地层进行参数统计分析来获取均匀岩性所对应的岩石基质体积模量,但是,对于物性非均匀、岩石结构复杂的碳酸盐岩地层,往往误差较大。
(4)等效换算法:该方法可利用岩石物理实验室测试得到干岩石骨架岩石物理参数,然后,以孔隙度、固结系数等为桥梁,通过换算得到等效的岩石基质体积模量。但是,由于近似公式受地域性影响较大,如果应用于其它的工区,估算结果的精度会受到较大影响,并且,岩石基质体积模量需要经过地区校正后才能使用。
综上所述,现有的获取岩石基质体积模量的方法无法准确而方便地获取精度和稳定性高的岩石基质体积模量。
发明内容
本发明示例性实施例的目的在于提供一种获取岩石基质体积模量的方法,以克服现有技术中难以准确而方便地获取岩石基质体积模量,并且获取到的岩石基质体积模量精度低和稳定性差的问题。
根据本发明示例性实施例的一方面,提供一种获取岩石基质体积模量的方法,包括:(A)利用目标测井深度点及其上下邻近的预定数量个测井深度点的测井资料确定目标测井深度点的等效岩石基质体积模量;(B)基于目标测井深度点的测井资料,确定目标测井深度点的干岩石泊松比的变化区间和变化增量;(C)基于目标测井深度点的测井资料,利用岩石骨架统一模型和岩石基质体积模量的相对关系确定目标测井深度点的岩石基质体积模量的变化区间和变化增量;(D)基于目标测井深度点的测井资料,在双收敛条件下利用自适应基质矿物模量反演方法反演得到目标测井深度点的最优自适应岩石基质体积模量,其中,反演所用到的自适应岩石基质体积模量在所述岩石基质体积模量的变化区间内依照岩石基质体积模量的变化增量进行取值,并且,反演所用到的干岩石泊松比在所述干岩石泊松比的变化区间内依照干岩石泊松比的变化增量进行取值;(E)将所述目标测井深度点的等效岩石基质体积模量和最优自适应岩石基质体积模量进行平均处理来获取目标测井深度点的岩石基质体积模量。
可选地,所述测井资料可包括以下项中至少一项:纵波速度、横波速度、密度、含水饱和度、孔隙流体体积模量以及岩石孔隙度。
可选地,步骤(A)可包括:(A1)利用目标测井深度点及其上下邻近的预定数量个测井深度点的测井资料分别确定目标测井深度点及其上下邻近的预定数量个测井深度点的饱和岩石体积模量;(A2)基于目标测井深度点及其上下邻近的预定数量个测井深度点各自对应的饱和岩石体积模量和岩石孔隙度,利用线性拟合法拟合得到饱和岩石体积模量与岩石孔隙度的线性关系式,并确定所述线性关系式的斜率和截距;(A3)利用所述线性关系式的斜率和截距确定目标测井深度点的等效岩石基质体积模量。
可选地,步骤(D)可包括:(D1)基于目标测井深度点的测井资料以及Gassmann方程和Gassmann-Boit-Geertsman方程确定Gassmann流体因子;(D2)基于目标测井深度点的测井资料以及Russell流体因子公式确定Russell流体因子;(D3)将Gassmann流体因子与Russell流体因子之差的绝对值作为第一反演目标函数,采用寻找全局最优解的方法反演得到满足第一收敛条件的包括自适应岩石基质体积模量的集合,其中,当第一反演目标函数的值小于第一预定阈值时,确定为满足第一收敛条件;(D4)将所述目标测井深度点的等效岩石基质体积模量与最优自适应岩石基质体积模量之差的绝对值作为第二反演目标函数,在所述包括自适应岩石基质体积模量的集合中采用寻找全局最优解的方法反演得到满足第二收敛条件的最优自适应岩石基质体积模量,其中,当第二反演目标函数的值小于第二预定阈值时,确定为满足第二收敛条件。
在根据本发明示例性实施例的获取岩石基质体积模量的方法中,能够准确且方便地获取具有较高精度和稳定性的岩石基质体积模量,为后续预测油气藏奠定了有利基础。
附图说明
通过下面结合附图进行的详细描述,本发明示例性实施例的上述和其它目的、特点和优点将会变得更加清楚,其中:
图1示出根据本发明示例性实施例的获取岩石基质体积模量的方法的流程图;
图2示出根据本发明示例性实施例的确定目标测井深度点的等效岩石基质体积模量的步骤的流程图;
图3示出根据本发明示例性实施例的利用自适应基质矿物模量反演方法反演得到目标测井深度点的最优自适应岩石基质体积模量的步骤的流程图;
图4示出根据本发明示例性实施例获取的岩石基质体积模量曲线与建模时输入的岩石基质体积模量的曲线的对比图;
图5示出根据本发明示例性实施例获取的岩石基质体积模量与建模输入模量的交会图。
具体实施方式
现将详细参照本发明的示例性实施例,所述实施例的示例在附图中示出,其中,相同的标号始终指的是相同的部件。以下将通过参照附图来说明所述实施例,以便解释本发明。
图1示出根据本发明示例性实施例的获取岩石基质体积模量的方法的流程图。这里,作为示例,所述方法可由用于获取岩石基质体积模量的设备来实现,也可以完全通过计算机程序来实现。
如图1所示,在步骤S100,利用目标测井深度点及其上下邻近的预定数量个测井深度点的测井资料确定目标测井深度点的等效岩石基质体积模量。这里,作为示例,所述测井资料可包括以下项中至少一项:纵波速度、横波速度、密度、含水饱和度、孔隙流体体积模量以及岩石孔隙度。下面,将结合图2来具体描述根据本发明示例性实施例的确定目标测井深度点的等效岩石基质体积模量的步骤的流程图。
图2示出根据本发明示例性实施例的确定目标测井深度点的等效岩石基质体积模量的步骤的流程图。如图2所示,在步骤S110,利用目标测井深度点及其上下邻近的预定数量个测井深度点的测井资料分别确定目标测井深度点及其上下邻近的预定数量个测井深度点的饱和岩石体积模量。这里,作为示例,可选取与目标测井深度点上下邻近的一个测井深度点,并分别确定目标测井深度点及其上下邻近的一个测井深度点所对应的饱和岩石体积模量。
具体说来,在步骤S110,可利用下面的等式(1-1)确定目标测井深度点的或目标测井深度点上下邻近的一个测井深度点的饱和岩石体积模量:
K s a t , i = ρ i ( V p , i 2 - 4 3 V s , i 2 ) , 等式(1-1)
其中,Ksat,i表示测井深度点i的饱和岩石体积模量,Vp,i表示测井深度点i的纵波速度,Vs,i表示测井深度点i的横波速度,ρi表示测井深度点i的密度。
这里,可通过等式(1-1)分别获取目标测井深度点及其上下邻近的一个测井深度点的饱和岩石体积模量。
接下来,在步骤S120,基于目标测井深度点及其上下邻近的预定数量个测井深度点各自对应的饱和岩石体积模量和岩石孔隙度,利用线性拟合法拟合得到饱和岩石体积模量与岩石孔隙度的线性关系式,并确定所述线性关系式的斜率和截距。
具体说来,可采用线性拟合法通过对Gassmann(加斯曼)理论进行合理简化,并引入Eshelly-Walsh干燥岩石椭球包体近似公式,利用目标测井深度点及其上下邻近的一个测井深度点的饱和岩石体积模量与岩石孔隙度,拟合推导出饱和岩石体积模量与岩石孔隙度的线性关系式,这里,所述线性关系式可用下面的等式(1-2)来表示:
φ i = A 1 K s a t , i - B , 等式(1-2)
其中,φi表示测井深度点i的岩石孔隙度,Ksat,i表示测井深度点i的饱和岩石体积模量,A表示线性关系式(1-2)所表示的直线的斜率,B表示线性关系式(1-2)所表示的直线的截距。
在步骤S130,利用所述线性关系式的斜率和截距确定目标测井深度点的等效岩石基质体积模量。
具体说来,由于线性关系式(等式(1-2))所表示的直线的斜率A和截距B可分别被表示为下面的等式:
A = K f , i + α m K o , i , 等式(1-3)
B = ( K f , i + α m K o , i ) K o , i , 等式(1-4)
其中,Kf,i表示测井深度点i的孔隙流体体积模量,Ko,i表示测井深度点i的岩石基质体积模量,m表示干燥岩石椭球体的结构参数,α表示岩石孔隙扁率。
因此,由等式(1-3)和等式(1-4)可以得出:
K o , i = A B , 等式(1-5)
其中,假设测井深度点i为目标测井深度点,则Ko,i可表示目标测井深度点i的等效岩石基质体积模量,A表示线性关系式(1-2)所表示的直线的斜率,B表示线性关系式(1-2)所表示的直线的截距。
通过等式(1-5)可确定通过线性拟合法获取的目标测井深度点的等效岩石基质体积模量。
再次返回图1,在步骤S200,基于目标测井深度点的测井资料,确定目标测井深度点的干岩石泊松比的变化区间和变化增量。
具体说来,基于目标测井深度点的测井资料,可确定目标测井深度点的干岩石泊松比的变化区间和变化增量,例如,基于目标测井深度点的测井资料,可确定目标测井深度点的干岩石泊松比的变化区间为[0.00.4],初始值可取变化区间的下限值0.0,变化为递增变化,向上的变化增量为0.03。
在步骤S300,基于目标测井深度点的测井资料,利用岩石骨架统一模型和岩石基质体积模量的相对关系确定目标测井深度点的岩石基质体积模量的变化区间和变化增量。
这里,可通过如下的方式来确定目标测井深度点的岩石基质体积模量的变化区间和变化增量:
具体说来,在岩石物理学中,岩石体积模量存在这样的相对关系:
Ko>Ksat>Kdry,不等式(3-1)
其中,Ko表示岩石基质体积模量,Ksat表示饱和岩石体积模量,Kdry表示干岩石骨架体积模量。这里,其中,φ表示岩石孔隙度,p、q表示调节函数,其中,p、q可由岩性条件和地区经验所确定。
由上面的分析可得:岩石体积模量的变化区间可表示为:
K s a t < K o = K d r y ( ( 1 + q &phi; ) ) ( 1 - p &phi; ) < K s a t ( 1 + q &phi; ) ( 1 - p &phi; ) , 不等式(3-2)
由不等式(3-2),可确定目标测井深度点的岩石基质体积模量的变化区间和变化增量,其中,初始值可选变化区间的上限值,变化为递减变化,变化增量为0.1Gpa。
通过上述处理,可获得目标测井深度点的岩石基质模量的变化区间和变化增量以及干岩石泊松比的变化区间和变化增量,为后续继续处理进一步限定了目标测井深度点的岩石基质体积模量和干岩石泊松比的取值范围。
在步骤S400,基于目标测井深度点的测井资料,在双收敛条件下利用自适应基质矿物模量反演方法反演得到目标测井深度点的最优自适应岩石基质体积模量,其中,反演所用到的自适应岩石基质体积模量在所述岩石基质体积模量的变化区间内依照岩石基质体积模量的变化增量进行取值,并且,反演所用到的干岩石泊松比在所述干岩石泊松比的变化区间内依照干岩石泊松比的变化增量进行取值。
下面,将结合图3来具体描述根据本发明示例性实施例的在双收敛条件下利用自适应基质矿物模量反演方法反演得到目标测井深度点的最优自适应岩石基质体积模量的步骤的流程图。
如图3所示,在步骤401,基于目标测井深度点的测井资料以及Gassmann方程和Gassmann-Boit-Geertsman方程确定Gassmann流体因子。
具体说来,Gassmann方程可表示为:
f G = ( 1 - K d r y K o ) &phi; K f + 1 - &phi; K o + K d r y K o , 公式(4-1)
其中,fG表示Gassmann流体因子,φ表示岩石孔隙度,Ko表示岩石基质体积模量,Kf表示孔隙流体体积模量,Kdry表示干岩石骨架体积模量,其中,Kdry=(1+β)×Ko,这里,β为中间变量,这里,β可基于目标测井深度点的测井资料应用Gassmann-Boit-Geertsman方程来获取:
这里,Gassmann-Boit-Geertsman方程可表示为:
( Y - 1 ) &beta; 2 + &lsqb; Y &phi; ( K o K f - 1 ) - Y + M K o &rsqb; &beta; - &phi; ( Y - M K o ) ( K o K f - 1 ) = 0 , 公式(4-2)
其中,其中,ρsat表示密度,Vp表示纵波速度,σdry表示干岩石泊松比,Kf表示孔隙流体体积模量,Ko表示岩石基质体积模量。
具体说来,可将目标测井深度点的相关测井资料(例如,目标测井深度点所对应的密度、纵波速度、孔隙流体体积模量等)代入相关公式,其中,与公式(4-2)相关的岩石基质体积模量或干岩石泊松比可在上述提及的目标测井深度点的岩石基质体积模量或干岩石泊松比的变化区间内,依照各自的变化增量,选取合适的岩石基质体积模量或干岩石泊松比作为已知参数来求取中间变量β,从而进一步获取目标测井深度点的干岩石骨架体积模量Kdry,然后,在公式(4-1)中以岩石基质体积模量Ko为未知参数,使得公式(4-1)可被表示成Gassmann流体因子fG关于岩石基质体积模量Ko的关系式。
另一方面,可采用另一种方法获取流体因子的不同表达式,这里,在步骤S402,可基于目标测井深度点的测井资料以及Russell流体因子公式确定Russell流体因子。
具体说来,Russell流体因子公式可表示为:
f R = ( Z p 2 - cZ s 2 ) &rho; s a t , 公式(4-3)
其中,fR表示Russell流体因子,Zp表示岩石纵波阻抗,Zs表示岩石横波阻抗、ρsat表示密度,c表示中间变量c系数,其中,Zp=ρsatVp,Zs=ρsatVsVp表示纵波速度,σdry表示干岩石泊松比。
具体说来,可利用目标测井深度点的相关测井资料(例如,目标测井深度点所对应的密度、纵波速度、横波速度等)代入相关公式,其中,干岩石泊松比可在上述提及的干岩石泊松比的变化区间内,依照变化增量,选取合适的干岩石泊松比作为已知参数来求取c系数,从而进一步代入公式(4-3),然后,在公式(4-3)中以干岩石泊松比σdry为未知参数,使得公式(4-3)可被表示成Russell流体因子fR关于干岩石泊松比σdry的关系式。
在步骤S403,将Gassmann流体因子与Russell流体因子之差的绝对值作为第一反演目标函数。
在步骤S404,分别在岩石基质体积模量的变化区间或干岩石泊松比的变化区间内按照各自的变化增量进行取值,并将选取到的第一组自适应岩石基质体积模量和干岩石泊松比代入第一反演目标函数。
在步骤S405,确定第一反演目标函数的值是否小于第一预定阈值。具体说来,可设置第一收敛条件,来缩小获取岩石基质体积模量的选取范围。这里,第一收敛条件指示第一反演目标函数的值小于第一预定阈值。
当第一反演目标函数的值小于第一预定阈值时,在步骤S406,获取满足第一收敛条件的自适应岩石基质体积模量。
接下来,在步骤S407,确定是否是分别在岩石基质体积模量的变化区间或干岩石泊松比的变化区间内按照各自的变化增量进行取值后得到的最后一组自适应岩石基质体积模量和干岩石泊松比。
当不是最后一组自适应岩石基质体积模量和干岩石泊松比时,在步骤S408,选取下一组自适应岩石基质体积模量和干岩石泊松比,并将该组自适应岩石基质体积模量和干岩石泊松比代入第一反演目标函数,并返回执行步骤S405,继续判断第一反演目标函数是否小于第一预定阈值。
当是最后一组自适应岩石基质体积模量和干岩石泊松比时,在步骤S409,获取满足第一收敛条件的包括自适应岩石基质体积模量的集合,这里,作为示例,在理想情况下,所述包括自适应岩石基质体积模量的集合可包括满足第一收敛条件的多组自适应岩石基质体积模量和干岩石泊松比。
在步骤S410,将所述目标测井深度点的等效岩石基质体积模量与最优自适应岩石基质体积模量之差的绝对值作为第二反演目标函数。
在步骤S411,在满足第一收敛条件的包括自适应岩石基质体积模量的集合中选取第一个自适应岩石基质体积模量并代入第二反演目标函数。
在步骤S412,确定第二反演目标函数的值是否小于第二预定阈值。具体说来,可设置第二收敛条件,来进一步缩小获取岩石基质体积模量的选取范围。
当第二反演目标函数的值小于第二预定阈值时,在步骤S413,获取满足第二收敛条件的最优自适应岩石基质体积模量。
在步骤S414,确定当前的自适应岩石基质体积模量是否是所述满足第一收敛条件的包括自适应岩石基质体积模量的集合中的最后一个自适应岩石基质体积模量。
当是最后一个自适应岩石基质体积模量时,结束所述方法。
当不是最后一个自适应岩石基质体积模量时,将下一个满足第一收敛条件的自适应岩石基质体积模量代入第二反演目标函数,并返回执行步骤S412,继续确定第二反演目标函数是否小于第二预定阈值。
再次返回参照图1,在步骤S500,将所述目标测井深度点的等效岩石基质体积模量和最优自适应岩石基质体积模量进行平均处理来获取目标测井深度点的岩石基质体积模量。
具体说来,可求取目标测井深度点的等效岩石基质体积模量和最优自适应岩石基质体积模量平均值作为目标测井深度点的岩石基质体积模量。
这里,应理解,针对预定工区目标测井深度段内的其他的目标测井深度点,同样可利用本发明所述的方法来获取预定工区目标测井深度段内所有目标测井深度点的岩石基质体积模量。
图4示出根据本发明示例性实施例获取的岩石基质体积模量曲线与建模时输入的岩石基质体积模量的曲线的对比图。
如图4所示,标记A指示的曲线为根据本发明示例性实施例所述的方法获取的岩石基质体积模量曲线,标记B指示的曲线为建模时输入的岩石基质体积模量曲线,由图4可以看出,根据本发明示例性实施例所述方法获取的岩石基质体积模量的曲线与建模时输入的岩石基质体积模量的总体特征完全一致,两者之差的绝对值介于0.0-2.5Gpa之间,因此,根据本发明示例性实施例所述方法获取的岩石基质体积模量具有较高的精度。
图5示出根据本发明示例性实施例获取的岩石基质体积模量与建模输入模量的交会图。这里,通过图5可以有效地验证根据本发明示例性实施例所述方法获取的岩石基质体积模量的正确性。
综上所述,在根据本发明示例性实施例的获取岩石基质体积模量的方法中,能够准确且方便地获取具有较高精度和稳定性的岩石基质体积模量,为后续预测油气藏奠定了有利基础。
尽管已经参照其示例性实施例具体显示和描述了本发明,但是本领域的技术人员应该理解,在不脱离权利要求所限定的本发明的精神和范围的情况下,可以对其进行形式和细节上的各种改变。

Claims (4)

1.一种获取岩石基质体积模量的方法,包括:
(A)利用目标测井深度点及其上下邻近的预定数量个测井深度点的测井资料确定目标测井深度点的等效岩石基质体积模量;
(B)基于目标测井深度点的测井资料,确定目标测井深度点的干岩石泊松比的变化区间和变化增量;
(C)基于目标测井深度点的测井资料,利用岩石骨架统一模型和岩石基质体积模量的相对关系确定目标测井深度点的岩石基质体积模量的变化区间和变化增量;
(D)基于目标测井深度点的测井资料,在双收敛条件下利用自适应基质矿物模量反演方法反演得到目标测井深度点的最优自适应岩石基质体积模量,其中,反演所用到的自适应岩石基质体积模量在所述岩石基质体积模量的变化区间内依照岩石基质体积模量的变化增量进行取值,并且,反演所用到的干岩石泊松比在所述干岩石泊松比的变化区间内依照干岩石泊松比的变化增量进行取值;
(E)将所述目标测井深度点的等效岩石基质体积模量和最优自适应岩石基质体积模量进行平均处理来获取目标测井深度点的岩石基质体积模量。
2.如权利要求1所述的方法,其中,所述测井资料包括以下项中至少一项:纵波速度、横波速度、密度、含水饱和度、孔隙流体体积模量以及岩石孔隙度。
3.如权利要求2所示的方法,其中,步骤(A)包括:
(A1)利用目标测井深度点及其上下邻近的预定数量个测井深度点的测井资料分别确定目标测井深度点及其上下邻近的预定数量个测井深度点的饱和岩石体积模量;
(A2)基于目标测井深度点及其上下邻近的预定数量个测井深度点各自对应的饱和岩石体积模量和岩石孔隙度,利用线性拟合法拟合得到饱和岩石体积模量与岩石孔隙度的线性关系式,并确定所述线性关系式的斜率和截距;
(A3)利用所述线性关系式的斜率和截距确定目标测井深度点的等效岩石基质体积模量。
4.如权利要求1所述的方法,其中,步骤(D)包括:
(D1)基于目标测井深度点的测井资料以及Gassmann方程和Gassmann-Boit-Geertsman方程确定Gassmann流体因子;
(D2)基于目标测井深度点的测井资料以及Russell流体因子公式确定Russell流体因子;
(D3)将Gassmann流体因子与Russell流体因子之差的绝对值作为第一反演目标函数,采用寻找全局最优解的方法反演得到满足第一收敛条件的包括自适应岩石基质体积模量的集合,其中,当第一反演目标函数的值小于第一预定阈值时,确定为满足第一收敛条件;
(D4)将所述目标测井深度点的等效岩石基质体积模量与最优自适应岩石基质体积模量之差的绝对值作为第二反演目标函数,在所述包括自适应岩石基质体积模量的集合中采用寻找全局最优解的方法反演得到满足第二收敛条件的最优自适应岩石基质体积模量,其中,当第二反演目标函数的值小于第二预定阈值时,确定为满足第二收敛条件。
CN201510410200.9A 2015-07-13 2015-07-13 获取岩石基质体积模量的方法 Active CN105093331B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510410200.9A CN105093331B (zh) 2015-07-13 2015-07-13 获取岩石基质体积模量的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510410200.9A CN105093331B (zh) 2015-07-13 2015-07-13 获取岩石基质体积模量的方法

Publications (2)

Publication Number Publication Date
CN105093331A true CN105093331A (zh) 2015-11-25
CN105093331B CN105093331B (zh) 2018-07-10

Family

ID=54574185

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510410200.9A Active CN105093331B (zh) 2015-07-13 2015-07-13 获取岩石基质体积模量的方法

Country Status (1)

Country Link
CN (1) CN105093331B (zh)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105370274A (zh) * 2015-12-14 2016-03-02 长江大学 井下地层孔隙度确定方法
CN105589110A (zh) * 2015-12-10 2016-05-18 中国石油集团川庆钻探工程有限公司 孔隙流体压缩系数判别致密砂岩气层的方法
CN106596908A (zh) * 2016-12-16 2017-04-26 中国水利水电科学研究院 混凝土等孔隙类材料的基质体积模量测定方法及装置
CN111323823A (zh) * 2019-12-27 2020-06-23 中国石油天然气股份有限公司 测井孔隙度曲线确定方法及系统
CN111368245A (zh) * 2020-03-12 2020-07-03 成都理工大学 一种考虑多影响因素的定量计算岩石Biot系数方法
CN112578456A (zh) * 2019-09-27 2021-03-30 中国石油化工股份有限公司 井中流体识别方法及系统

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102096107A (zh) * 2009-12-09 2011-06-15 中国石油天然气股份有限公司 一种根据声波时差和密度反演孔隙扁度进行储层渗透性评价的方法
CN102353989A (zh) * 2011-08-24 2012-02-15 成都理工大学 基于自适应基质矿物等效弹性模量反演的横波速度估算方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102096107A (zh) * 2009-12-09 2011-06-15 中国石油天然气股份有限公司 一种根据声波时差和密度反演孔隙扁度进行储层渗透性评价的方法
CN102353989A (zh) * 2011-08-24 2012-02-15 成都理工大学 基于自适应基质矿物等效弹性模量反演的横波速度估算方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
HE XI-LEI 等: "Calculations of rock matrix modulus based on a linear regression relation", 《APPLIED GEOPHYSICS》 *
LIN KAI 等: "Self-adapting extraction of matrix mineral bulk modulus and verification of fluid substitution", 《APPLIED GEOPHYSICS》 *
熊晓军 等: "基于等效弹性模量反演的横波速度预测方法", 《石油地球物理勘探》 *

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105589110A (zh) * 2015-12-10 2016-05-18 中国石油集团川庆钻探工程有限公司 孔隙流体压缩系数判别致密砂岩气层的方法
CN105370274A (zh) * 2015-12-14 2016-03-02 长江大学 井下地层孔隙度确定方法
CN106596908A (zh) * 2016-12-16 2017-04-26 中国水利水电科学研究院 混凝土等孔隙类材料的基质体积模量测定方法及装置
CN106596908B (zh) * 2016-12-16 2018-11-20 中国水利水电科学研究院 孔隙类材料的基质体积模量测定方法及装置
CN112578456A (zh) * 2019-09-27 2021-03-30 中国石油化工股份有限公司 井中流体识别方法及系统
CN111323823A (zh) * 2019-12-27 2020-06-23 中国石油天然气股份有限公司 测井孔隙度曲线确定方法及系统
CN111323823B (zh) * 2019-12-27 2022-05-10 中国石油天然气股份有限公司 测井孔隙度曲线确定方法及系统
CN111368245A (zh) * 2020-03-12 2020-07-03 成都理工大学 一种考虑多影响因素的定量计算岩石Biot系数方法
CN111368245B (zh) * 2020-03-12 2023-04-14 成都理工大学 一种考虑多影响因素的定量计算岩石Biot系数方法

Also Published As

Publication number Publication date
CN105093331B (zh) 2018-07-10

Similar Documents

Publication Publication Date Title
CN105093331B (zh) 获取岩石基质体积模量的方法
CN109655903B (zh) 页岩层横波速度预测方法及系统
US8359184B2 (en) Method, program and computer system for scaling hydrocarbon reservoir model data
US8670288B2 (en) Velocity model for well time-depth conversion
US20130229892A1 (en) Method of predicting the pressure sensitivity of seismic velocity within reservoir rocks
US20030028325A1 (en) Method of constraining by dynamic production data a fine model representative of the distribution in the reservoir of a physical quantity characteristic of the subsoil structure
AU2013399120B2 (en) Static earth model calibration methods and systems
CN103982179A (zh) 一种油藏储层的古压力定量反演探测方法
RU2709047C1 (ru) Способ адаптации гидродинамической модели продуктивного пласта нефтегазоконденсатного месторождения с учетом неопределенности геологического строения
CN104272140A (zh) 用于校准在储藏层建模中使用的渗透性的系统和方法
CN103913774A (zh) 基于微地震事件的储层地质力学参数反演方法
US8599644B2 (en) Velocity models for a single well and for a set of wells
CN107797139A (zh) 页岩储层游离气含气量地震预测方法及系统
CN113552621B (zh) 页岩气地应力确定方法和装置
CN105223616A (zh) 一种页岩储层的孔隙纵横比反演方法
CN104632206B (zh) 一种石油测井标准层确定方法及装置
CN111967677B (zh) 一种非常规资源甜点分布的预测方法及装置
CN109839676A (zh) 一种基于记忆模拟退火的基质模量估算方法及电子设备
NO20200978A1 (en) Optimized methodology for automatic history matching of a petroleum reservoir model with ensemble kalman filter
CN114563823B (zh) 一种储层岩性和物性参数的联合贝叶斯反演方法和系统
CN110954949A (zh) 一种致密砂岩软孔隙度分布反演方法
CN110927787A (zh) 基于虚拟岩性的横波速度计算方法及系统
CN105893674B (zh) 采用全局协方差进行地质属性预测的方法
CN114076987B (zh) 基于致密砂岩油气藏有效储层的测井识别方法及装置
CN111381279B (zh) 储层孔隙度定量预测方法及装置

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
TA01 Transfer of patent application right

Effective date of registration: 20180418

Address after: No. 189, fan Yangxi Road, Zhuozhou, Hebei Province, Hebei

Applicant after: BGP INC., CHINA NATIONAL PETROLEUM Corp.

Address before: Shuangliu County Huayang Huayang Road in Chengdu city of Sichuan Province in 610213 section of No. 216, Igawa geophysical exploration company of the Ministry of science and technology

Applicant before: GEOPHYSICAL EXPLORATION COMPANY OF CNPC CHUANQING DRILLING ENGINEERING Co.,Ltd.

TA01 Transfer of patent application right
GR01 Patent grant
GR01 Patent grant
TR01 Transfer of patent right

Effective date of registration: 20201111

Address after: 100007 Beijing, Dongzhimen, North Street, No. 9, No.

Patentee after: CHINA NATIONAL PETROLEUM Corp.

Patentee after: BGP Inc., China National Petroleum Corp.

Address before: No. 189, fan Yangxi Road, Zhuozhou, Hebei Province, Hebei

Patentee before: BGP Inc., China National Petroleum Corp.

TR01 Transfer of patent right