CN102288986A - 一种地震尺度下碳酸盐岩储层弹性模量获取方法 - Google Patents
一种地震尺度下碳酸盐岩储层弹性模量获取方法 Download PDFInfo
- Publication number
- CN102288986A CN102288986A CN2011101263994A CN201110126399A CN102288986A CN 102288986 A CN102288986 A CN 102288986A CN 2011101263994 A CN2011101263994 A CN 2011101263994A CN 201110126399 A CN201110126399 A CN 201110126399A CN 102288986 A CN102288986 A CN 102288986A
- Authority
- CN
- China
- Prior art keywords
- rock
- elastic modulus
- pore
- sub
- modulus
- 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
Images
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明实施例提供了一种地震尺度下碳酸盐岩储层弹性模量获取方法,方法包括:将非均匀碳酸盐岩储层岩石网格化以划分成有限的岩石子块,获取每一所述岩石子块的独立物理参数;根据差分有效介质模型获取含粒内孔隙的所述岩石子块的骨架弹性模量;根据硬岩模型获取含粒内孔隙与粒间孔隙的所述岩石子块的干岩石弹性模量;根据改进的差分有效介质模型获取含粒内孔隙、粒间孔隙与裂缝的所述岩石子块的干岩石弹性模量;根据Gassmann方程获取含粒内孔隙、粒间孔隙与裂缝的所述岩石子块的饱和流体岩石弹性模量;根据边界方法获取地震尺度下碳酸盐岩储层弹性模量边界与弹性模量。本发明实施例能较高精度地实现地震尺度下碳酸盐岩储层弹性参数的获取。
Description
技术领域
本发明涉及碳酸盐岩油气田地震勘探与油藏预测领域,特别是地震尺度下碳酸盐岩储层弹性模量获取方法。
背景技术
碳酸盐岩储层是世界油气勘探的最重要目标之一。从碳酸盐岩储层中发现的油气储量已接近世界油气储量的一半,产量则已达总产量的60%以上。碳酸盐岩储层的主要岩石类型包括石灰岩、白云岩、粒屑灰岩、礁灰岩等,其储集空间通常包括孔隙、溶洞和裂缝三类。孔隙和溶洞是主要的储集空间,裂缝是主要的渗滤通道,也是储集空间。碳酸盐岩储集空间的形成过程是一个复杂而长期的过程,它贯穿在整个沉积过程及其以后的各个地质历史时期。它除了受沉积环境的控制外,地下热动力场、地下或地表水化学场、构造应力场等因素均对它们的形成和发展有巨大的影响。由于碳酸盐岩的特殊性(易溶性和不稳定性),使碳酸盐岩储集空间的演化相当复杂,孔隙类型多、变化快,在同一储集层内存在多种类型的孔隙,各种孔隙又往往经受几种因素的作用和改造。因此碳酸盐岩储层具有很强的非均质型,这使得其地震尺度下弹性模量的获取十分困难,而如果没有合理的碳酸盐岩储层弹性模量获取方法就难以提高碳酸盐岩储层预测与描述精度。
目前,碳酸盐岩储层弹性模量获取方法研究存在的主要问题是:
1、碳酸盐岩储层弹性模量获取方法都基于岩芯和测井数据,不能考虑地震尺度下碳酸盐岩储层的非均质性,难以应用于地震尺度下的储层测量结果。
2、目前碳酸盐岩储层弹性模量获取方法不能对碳酸盐岩储层孔隙类型进行成因、特征详细分类分析,通常根据主导孔隙类型简化为裂缝介质或孔隙介质,而考虑孔隙类型影响只限于高频信号,如喷射流等,往往在地震频段难以发生或影响很小。
3、对岩石模量采用平均模量描述,如平均孔隙度、平均孔隙纵横比,不能很好表征岩石物理模量变化特征。
发明内容
本发明实施例提供一种地震尺度下碳酸盐岩储层弹性模量获取方法用于弥补现有技术中的上述缺陷。
本发明实施例提供了一种地震尺度下碳酸盐岩储层弹性模量获取方法,包括:将非均匀碳酸盐岩储层岩石网格化以划分成有限的岩石子块,获取每一所述岩石子块的独立物理参数;利用所述独立物理参数并根据差分有效介质模型获取含粒内孔隙的所述岩石子块的骨架弹性模量;利用所述独立物理参数以及所述含粒内孔隙岩石子块的骨架弹性模量并根据硬岩模型获取含粒内孔隙与粒间孔隙的所述岩石子块的干岩石弹性模量;利用所述独立物理参数以及所述含粒内、粒间孔隙岩石子块的骨架弹性模量并根据改进的差分有效介质模型获取含粒内孔隙、粒间孔隙与裂缝的所述岩石子块的干岩石弹性模量;利用所述独立物理参数以及所述含粒内、粒间孔隙与裂缝岩石子块的骨架弹性模量并根据Gassmann方程获取含粒内孔隙、粒间孔隙与裂缝的所述岩石子块的饱和流体岩石弹性模量;利用所述含粒内孔隙、粒间孔隙与裂缝的所述岩石子块的饱和流体岩石弹性模量并根据边界方法获取地震尺度下碳酸盐岩储层弹性模量边界与弹性模量。
优选的,本发明实施例中根据边界方法获取地震尺度下碳酸盐岩储层弹性模量边界与弹性模量包括:利用改进的Hashin-Shtrikman边界方法获取地震尺度下饱和碳酸盐岩弹性模量上、下边界;采用平均边界方法计算地震尺度下饱和碳酸盐岩弹性模量。
优选的,本发明实施例中获取每一所述岩石子块的独立物理参数包括:利用有限测井数据和岩心数据为输入通过Monte-Carlo方法模拟得到每一所述岩石子块的独立物理参数。
优选的,本发明实施例中根据差分有效介质模型获取含粒内孔隙的所述岩石子块的骨架弹性模量包括:将粒内孔隙的孔隙纵横比设置为1,含水饱和度设置为100%,根据差分有效介质模型获取含粒内孔隙的所述岩石子块的骨架弹性模量。
优选的,本发明实施例中根据硬岩模型获取含粒内孔隙与粒间孔隙的所述岩石子块的干岩石弹性模量包括:利用如下公式获得胶结砂岩干岩石有效体积模量Keff和剪切模型Geff:
其中φ和φ0为岩石孔隙度和临界孔隙度,K和G为岩石基质体积模量和剪切模量;以及,利用如下公式获得干岩石在临界孔隙度φ0时的体积模量KHM和剪切模量GHM:
其中n为颗粒接触数,G为岩石颗粒的剪切模量,v为岩石颗粒的泊松比K为岩石颗粒的体积模量,P为有效压力。
优选的,本发明实施例中根据Gassmann方程获取含粒内孔隙、粒间孔隙与裂缝的所述岩石子块的饱和流体岩石弹性模量包括:利用下述公式获取饱和流体岩石体积模量KWet和饱和流体岩石剪切模量GWet:
GWet=GDry;
其中KS、KDry和Kf分别为岩石基质、干岩石和孔隙流体体积模量,GDry为干岩石剪切模量。
本发明实施例由于采取以上技术方案,因此具有以下优点:1、采用岩石网格化方法和储层岩石子块的独立参数表征,从而实现地震尺度下非均匀油藏的有效表征,并为地震尺度下非均匀油藏流体替换奠定基础;2、根据岩石孔隙成因、结构特征和孔隙联通性以及孔隙之间相互作用,采用不同的合适模型并进行有效组合进行岩石物理弹性参数计算,提高岩石物理弹性参数的计算精度。3、采用边界计算方法将岩石子块组合有效弹性模量上、下边界,并采用平均边界方法实现地震尺度下碳酸盐岩储层弹性参数的获取。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的一种地震尺度下碳酸盐岩储层弹性模量获取方法流程示意图;
图2a是孔隙度网格化处理前的示意图;
图2b是孔隙度网格化处理后的示意图;
图3a-3b为在岩石粒内孔隙均为0.05,粒间孔隙平均值为0.15(图3a)和0.05两种情况下,岩石子块体积模量和剪切模量随孔隙度变化情况示意图;
图4a-4b为在岩石粒内孔隙均为0.05,粒间孔隙平均值为0.15、裂缝孔隙平均值为0.05和粒间孔隙平均值为0.05、裂缝孔隙平均值为0.15两种情况下,岩石子块体积模量和剪切模量随孔隙度变化情况示意图;
图5a-5b为在岩石粒内孔隙均为0.05,粒间孔隙平均值为0.15、裂缝孔隙平均值为0.05和粒间孔隙平均值为0.05、裂缝孔隙平均值为0.15两种情况下,岩石子块流体替换前体积模量和流体替换后体积模量随孔隙度变化情况示意图;
图6a-6b为在岩石粒内孔隙均为0.05,粒间孔隙平均值为0.15、裂缝孔隙平均值为0.05和粒间孔隙平均值为0.05、裂缝孔隙平均值为0.15两种情况下,岩石子块体积模量与孔隙度(粒间空隙+裂缝孔隙)、裂缝孔隙纵横比以及泥质含量的关系示意图;
图7为在岩石粒间孔隙平均值为0.05、裂缝孔隙平均值为0.15情况下,岩石子块流体替换前后体积模量与孔隙度(粒间空隙+裂缝孔隙)、裂缝孔隙纵横比的关系示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
由于碳酸盐岩的特殊性(易溶性和不稳定性),使碳酸盐岩储集空间的演化相当复杂,孔隙类型多、变化快,在同一储集层内存在多种类型的孔隙,各种孔隙又往往经受几种因素的作用和改造。本发明将碳酸盐岩石孔隙根据其岩石物理特点分为三个大类,即粒内孔隙、粒间孔隙和裂缝孔隙,其的目是根据孔隙结构、孔隙连通性和孔隙间相互作用采用不同类型岩石物理模型进行描述。并采用岩石网格化方法,将地震尺度下非均质碳酸盐岩储层岩石划分为有限块具有独立岩石物理参数的均质岩石子块,根据岩石孔隙特征采用不同岩石物理模型分步获得岩石子块干岩石与饱和流体岩石弹性参数,并采用弹性边界计算方法实现地震尺度下碳酸盐岩储层弹性参数的获取。
如图1所示为本发明实施例提供的一种地震尺度下碳酸盐岩储层弹性模量获取方法流程示意图,该方法包括如下步骤:
S101:将非均匀碳酸盐岩储层岩石网格化以划分成有限的岩石子块,获取每一所述岩石子块的独立物理参数。
地震尺度下碳酸盐岩孔隙度、孔隙结构和孔隙纵横比等参数都是非均匀的,难以统一描述。为此将碳酸盐岩储层岩石通过网格化划分成有限的岩石子块,而每一岩石子块等效为均匀体,且具有独立的岩石物理参数进行表征。每一岩石子块的独立物理参数包括粒内孔隙度、粒间孔隙度、裂缝孔隙度、裂缝孔隙纵横比以及泥质含量等,而参数的来源可以利用有限测井数据和岩心数据为输入通过Monte-Carlo方法模拟得到,模拟过程中充分考虑不同油藏参数之间的相关性,从而得到合理的岩石子块参数组合。
在本实施例中,岩石网格化过程如图2a和图2b所示,其中图2a是孔隙度网格化处理前的示意图,而图2b是孔隙度网格化处理后的示意图,其中图中明暗度代表岩石孔隙。在本实施例中,可以将孔隙度非均匀碳酸岩石网格化为10*10*10个具有不同孔隙度的均匀岩石子块组合,而每一岩石子块可以看做是足够大均匀岩石的一部分,从而进行弹性模量计算。其他参数的网格化可以参考孔隙度网格化进行。关于描述岩石子块的物理参数,根据该参数平均值,如平均孔隙度,和岩石该物理参数分布特征,如孔隙纵横比分布特征,采用Monte-Carlo方法产生,产生过程中充分考虑不同岩石物理参数间的相关性,如泥质含量与裂缝的负相关性。Monte-Carlo模拟是一种通过设定随机过程,反复生成时间序列,计算参数估计量和统计量,进而研究其分布特征的方法。具体的,当***中各个单元的可靠性特征量已知,但***的可靠性过于复杂,难以建立可靠性预计的精确数学模型或模型太复杂而不便应用时,可用随机模拟法近似计算出***可靠性的预计值;随着模拟次数的增多,其预计精度也逐渐增高。
S102:利用所述独立物理参数并根据差分有效介质模型获取含粒内孔隙的所述岩石子块的骨架弹性模量。
对于相对独立的粒内孔隙,考虑其形成原因、孔隙结够特征和孔隙联通性以及孔隙之间相互作用,将其作为岩石颗粒的一部分,并且将其孔隙纵横比规定为1,含水饱和度100%,利用差分有效介质模型计算含粒内孔隙岩石骨架弹性参数,由于岩石粒内孔隙为外界极少联通,因此在进行流体替换计算时,不对粒内孔隙进行流体替换。基于岩石子块的独立岩石物理参数,获得所用岩石子块对应含粒内孔隙岩石骨架弹性模量。
差分有效介质模型假定复杂介质可以通过对现有复杂介质通过无限小改变获得,如果初始介质体积模量和剪切模量分别为K*(y)和G*(y),y为被包含物质的体积百分含量,则粒内孔隙包含地下盐水介质的体积模量K2和剪切模量G2满足下面的方程(1)和(2):
式中P和Q是表征新包含物质(岩石孔隙)几何特征的因子。
S103:利用所述独立物理参数以及所述含粒内孔隙岩石子块的骨架弹性模量并根据硬岩模型获取含粒内孔隙与粒间孔隙的所述岩石子块的干岩石弹性模量。
碳酸盐岩粒间孔隙形成过程、特征与砂岩油藏相似,因此采用硬岩模型计算含粒间孔隙碳酸盐岩子块干岩石模量,而组成岩石子块骨架颗粒的弹性模量是由上一步得到的。
粒间孔隙形成过程及结构特征与碎屑砂岩储层岩石相似,因此可以利用硬岩模型计算含粒间孔隙碳酸盐岩弹性模型,硬岩模型基于Hashin-Shtrikman的上边界推导得到,可以获得胶结砂岩干岩石有效体积模量(Keff)和剪切模型(Geff)。岩石物理模型方程如方程(3)和(4)所示:
方程中φ和φ0为岩石孔隙度和临界孔隙度,K和G为岩石基质体积模和剪切模量,KHM和GHM为一定有效压力下随机排列同规格球堆体积模量与剪切模量,可以利用如式(5)和式(6)所示Hertz-Mindlin接触理论模型计算。
其中,KHM,GHM是干岩石在临界孔隙度φ0时的体积模量和剪切模量,n为颗粒接触数,G为岩石颗粒的剪切模量,v为岩石颗粒的泊松比K为岩石颗粒的体积模量,P为有效压力。图3显示了岩石粒内孔隙均为0.05,而粒间孔隙平均值为0.15(图3a)和0.05(图3b)两种情况下,岩石子块体积模量(图中浅色圆点)和剪切模量(图中深色圆点)随孔隙度变化情况。
S104:利用所述独立物理参数以及所述含粒内、粒间孔隙岩石子块的骨架弹性模量并根据改进的差分有效介质模型获取含粒内孔隙、粒间孔隙与裂缝的所述岩石子块的干岩石弹性模量。
对于变化特征复杂的碳酸盐岩储层裂缝,可等效为均匀储层岩石子块嵌入不同纵横比裂缝,嵌入裂缝后岩石弹性可以利用改进的有效介质模型进行计算,而裂缝嵌入前的母岩弹性参数是由上一步基于硬岩模型计算的含粒内与粒间孔隙碳酸盐岩子块干岩石弹性模量。由于差分有效介质模型没有考虑裂缝孔隙之间的相互作用,因此在母岩中切入的裂缝为干裂缝,及体积与剪切模量都为0,从而得到干岩石子块含粒内、粒间和裂缝时的弹性模量。
改进的差分有效介质理论是将包含物质由原来的流体改变为具有临界孔隙度ΦC的临界相物质,其弹性模量为KC和GC,而总孔隙度为Φ=yΦC。因此有效介质弹性模量计算过程为:
式中P和Q是依然为表征新包含物质(岩石孔隙)几何特征的因子,但在计算中将K2和G2相应替换为KC和GC。
图4显示了岩石粒内孔隙均为0.05,粒间孔隙平均值为0.15、裂缝孔隙平均值为0.05(图4a)和粒间孔隙平均值为0.05、裂缝孔隙平均值为0.15(图4b)两种情况下,岩石子块体积模量(图中浅色圆点)和剪切模量(图中深色圆点)随孔隙度变化情况。左图和右图虽然具有相同的孔隙度,但由于左图所示岩石孔隙以粒间孔隙为主而右图所示岩石孔隙以裂缝孔隙为主,因此岩石子块干岩石弹性模量分布不同。这与实验室岩石物理测量特征是相一致的。
S105:利用所述独立物理参数以及所述含粒内、粒间孔隙与裂缝岩石子块的骨架弹性模量并根据Gassmann方程获取含粒内孔隙、粒间孔隙与裂缝的所述岩石子块的饱和流体岩石弹性模量。
基于上面得到的含粒内孔隙、粒间孔隙与裂缝的碳酸盐岩储层子块干岩石弹性模量,利用Gassmann方程计算储层饱和流体时子块岩石弹性模量,值得注意的是基于Gassmann方程进行流体替换计算时,由于粒内孔隙相对独立,难以与其它孔隙联通,且通常饱和水,因此进行流体替换时,该部分孔隙不进行流体替换。对于低孔低渗岩石子块也可以根据实际油藏开发认识确定是否进行流体替换。
岩石饱和不同流体时体积模量与剪切模量采用Gassmann方程进行计算,如方程(9)和(10)所示。
GWet=GDry (10)
方程中,KWet、KS、KDry和Kf分别为饱和流体岩石体积模量、岩石基质、干岩石和孔隙流体体积模量,GWet和GDry分别为饱和流体岩石和干岩石剪切模量。实际流体替换过程中,由于粒内孔隙相对独立,难以与其它孔隙联通,且通常饱和水,因此进行流体替换时,该部分孔隙不进行流体替换。对于低孔低渗岩石子块,由于流体难以驱动,因此也可以根据实际油藏开发认识确定是否对孔隙度或裂缝孔隙度低的岩石子块进行流体替换,以保证岩石流体替换情况与实际油藏储存开发情况一致。
S106:利用所述含粒内孔隙、粒间孔隙与裂缝的所述岩石子块的饱和流体岩石弹性模量并根据边界方法获取地震尺度下碳酸盐岩储层弹性模量边界与弹性模量。
在确定碳酸盐岩储层岩石子块的弹性模量后,利用改进的Hashin-Shtrikman边界方法计算地震尺度下岩石整体弹性模量上下边界。由于实验表明上下边界差异较小,因此利用平均边界方法确定岩石弹性模量,从而实现了地震尺度下碳酸盐岩储层弹性参数计算。
岩石弹性参数与岩石孔隙度并不是简单的线性关系,它依靠于岩石颗粒、孔隙的几何形状、粘土含量等很多因素。为了估计岩石弹性模量的可能变化范围,Hasin-Shtrikman边界方法可以估计岩石弹性模量的上下边界,其中对于计算岩石体积模量的方程如下:
Km和KP分别为岩石矿物和岩石岩石孔隙填充物的体积模量,Φ为岩石的孔隙度,K(HS)为岩石体积模量的边界值,计算岩石体积模量上边界时有:
其中,Gm为岩石矿物的剪切模量。计算岩石体积模量下边界时有:
其中,Gp为岩石孔隙填充物的剪切模量。对于计算岩石的剪切模量公式有:
其中G(HS)为岩石体积模量的边界值。计算岩石剪切模量上边界时有:
计算岩石剪切模量下边界时有:
同时根据以上计算方法,获得饱和水岩石弹性模量上、下边界。
对于具有多种弹性体组成的复杂岩石,可以采用广义的Hashin-Shtrikman-Walpole边界条件的进行计算,其计算公式如下所示:
KHS+=Λ(Gmax)KHS-=Λ(Gmin)
GHS+=Γ(ζ(Kmax,Gmax))GHS-=Γ(ζ(Kmin,Gmin)) (17)
上式中KHS+、KHS-、GHS+和GHS-分别表示复杂岩石体积模量上、下边界和剪切模量上下边界,K和G分别表示岩石子块体积模量和剪切模量,「表示取平均值。
对于给定岩石,其孔隙度一定,弹性模量的上、下界含有岩石孔隙强度的信息,即弹性模量越接近上边界,说明岩石孔隙强度大;反之,说明岩石的胶结程度低,岩石孔隙强度小。因此,可以用平均边界描述岩石的弹性模量,其计算关系式可写为:
K=K-+ω(K+-K-)(18)
G=G-+λ(G+-G-)(19)
式中,K+,K-和G+,G-分别表示岩石体积模量和剪切模量的上、下边界;K,G为所求岩石的体积模量与剪切模量;ω,λ为权因子,它们的取值在0和1之间,是表示岩石孔隙强度的参数。对于同一岩石,弹性模量的上、下边界取决于流体的弹性模量。但是流体种类对岩石孔隙结构的强度不会产生影响。因此,对于同一岩石,权因子ω,λ是不变的值。这样,可以根据实验室岩芯测量和测井数据得到的弹性模量来计算权因子ω,λ,并利用其计算其它流体饱和岩石的弹性模量。
图5显示了岩石粒内孔隙均为0.05,粒间孔隙平均值为0.15、裂缝孔隙平均值为0.05(图5a)和粒间孔隙平均值为0.05、裂缝孔隙平均值为0.15(图5b)两种情况下,岩石子块流体替换前(图中深色圆点)体积模量和流体替换后(图中浅色圆点)体积模量随孔隙度变化。图中浅色正三角代表流体替换前岩石体积模量下界;深色正三角代表流体替换前岩石体积模量上界;浅色倒三角代表流体替换后岩石体积模量下界;深色倒三角代表流体替换后岩石体积模量上界。从图中可以看出岩石上下边界差异不大,因此利用平均边界方法计算地震尺度下碳酸盐岩岩石弹性模量误差不大。对比图5a和图5b可以看出,虽然岩石总体平均孔隙度相同,但由于5a所示岩石孔隙以粒间孔隙为主,而图5b所示岩石孔隙以裂缝孔隙为主,因此两者的弹性不同,而岩石子块弹性参数分布特征也不相同,裂缝为主岩石弹性参数分布更为分散,这与实际油藏状况是相符合的。
基于本发明实施例的方法,可以对不同岩石物理参数对碳酸盐岩岩石弹性参数影响进行多维分析,比如,图6显示了岩石粒内孔隙均为0.05,粒间孔隙平均值为0.15、裂缝孔隙平均值为0.05(图6a)和粒间孔隙平均值为0.05、裂缝孔隙平均值为0.15(图6b)两种情况下,岩石子块体积模量与孔隙度(粒间空隙+裂缝孔隙)、裂缝孔隙纵横比以及泥质含量的关系。图7显示了岩石粒间孔隙平均值为0.05、裂缝孔隙平均值为0.15情况下,岩石子块流体替换前后体积模量与孔隙度(粒间空隙+裂缝孔隙)、裂缝孔隙纵横比的关系。
本发明实施例由于采取以上技术方案,因此具有以下优点:1、采用岩石网格化方法和储层岩石子块的独立参数表征,从而实现地震尺度下非均匀油藏的有效表征,并为地震尺度下非均匀油藏流体替换奠定基础;2、根据岩石孔隙成因、结构特征和孔隙联通性以及孔隙之间相互作用,采用不同的合适模型并进行有效组合进行岩石物理弹性参数计算,提高岩石物理弹性参数的计算精度。3、采用边界计算方法将岩石子块组合有效弹性模量上、下边界,并采用平均边界方法实现地震尺度下碳酸盐岩储层弹性参数的获取。
以上所述的具体实施例,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施例而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (6)
1.一种地震尺度下碳酸盐岩储层弹性模量获取方法,其特征在于,所述方法包括:
将非均匀碳酸盐岩储层岩石网格化以划分成有限的岩石子块,获取每一所述岩石子块的独立物理参数;
利用所述独立物理参数并根据差分有效介质模型获取含粒内孔隙的所述岩石子块的骨架弹性模量;
利用所述独立物理参数以及所述含粒内孔隙岩石子块的骨架弹性模量并根据硬岩模型获取含粒内孔隙与粒间孔隙的所述岩石子块的干岩石弹性模量;
利用所述独立物理参数以及所述含粒内、粒间孔隙岩石子块的骨架弹性模量并根据改进的差分有效介质模型获取含粒内孔隙、粒间孔隙与裂缝的所述岩石子块的干岩石弹性模量;
利用所述独立物理参数以及所述含粒内、粒间孔隙与裂缝岩石子块的骨架弹性模量并根据Gassmann方程获取含粒内孔隙、粒间孔隙与裂缝的所述岩石子块的饱和流体岩石弹性模量;
利用所述含粒内孔隙、粒间孔隙与裂缝的所述岩石子块的饱和流体岩石弹性模量并根据边界方法获取地震尺度下碳酸盐岩储层弹性模量边界与弹性模量。
2.如权利要求1所述的方法,其特征在于,所述根据边界方法获取地震尺度下碳酸盐岩储层弹性模量边界与弹性模量包括:
利用改进的Hashin-Shtrikman边界方法获取地震尺度下饱和碳酸盐岩弹性模量上、下边界;
采用平均边界方法计算地震尺度下饱和碳酸盐岩弹性模量。
3.如权利要求2所述的方法,其特征在于,所述获取每一所述岩石子块的独立物理参数包括:
利用有限测井数据和岩心数据为输入通过Monte-Carlo方法模拟得到每一所述岩石子块的独立物理参数。
4.如权利要求2所述的方法,其特征在于,所述根据差分有效介质模型获取含粒内孔隙的所述岩石子块的骨架弹性模量包括:
将粒内孔隙的孔隙纵横比设置为1,含水饱和度设置为100%,根据差分有效介质模型获取含粒内孔隙的所述岩石子块的骨架弹性模量。
5.如权利要求2所述的方法,其特征在于,所述根据硬岩模型获取含粒内孔隙与粒间孔隙的所述岩石子块的干岩石弹性模量包括:
利用如下公式获得胶结砂岩干岩石有效体积模量Keff和剪切模型Geff:
其中φ和φ0为岩石孔隙度和临界孔隙度,K和G为岩石基质体积模量和剪切模量;以及,
利用如下公式获得干岩石在临界孔隙度φ0时的体积模量KHM和剪切模量GHM:
其中n为颗粒接触数,G为岩石颗粒的剪切模量,v为岩石颗粒的泊松比K为岩石颗粒的体积模量,P为有效压力。
6.如权利要求2所述的方法,其特征在于,所述根据Gassmann方程获取含粒内孔隙、粒间孔隙与裂缝的所述岩石子块的饱和流体岩石弹性模量包括:
利用下述公式获取饱和流体岩石体积模量KWet和饱和流体岩石剪切模量GWet:
GWet=GDry;
其中KS、KDry和Kf分别为岩石基质、干岩石和孔隙流体体积模量,GDry为干岩石剪切模量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110126399 CN102288986B (zh) | 2011-05-16 | 2011-05-16 | 一种地震尺度下碳酸盐岩储层弹性模量获取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110126399 CN102288986B (zh) | 2011-05-16 | 2011-05-16 | 一种地震尺度下碳酸盐岩储层弹性模量获取方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102288986A true CN102288986A (zh) | 2011-12-21 |
CN102288986B CN102288986B (zh) | 2013-01-30 |
Family
ID=45335548
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201110126399 Expired - Fee Related CN102288986B (zh) | 2011-05-16 | 2011-05-16 | 一种地震尺度下碳酸盐岩储层弹性模量获取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102288986B (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105468886A (zh) * | 2014-09-04 | 2016-04-06 | 中国石油化工股份有限公司 | 基于岩石物理学参数计算地层压力的方法 |
CN106251321A (zh) * | 2015-06-12 | 2016-12-21 | 中国石油化工股份有限公司 | 用于获取岩石裂缝信息的方法和装置 |
CN106324672A (zh) * | 2015-07-03 | 2017-01-11 | 中国石油化工股份有限公司 | 一种岩石物理建模方法及*** |
CN108181192A (zh) * | 2017-12-25 | 2018-06-19 | 万向钱潮(上海)汽车***有限公司 | 一种制动盘的弹性模量获取方法 |
CN109100801A (zh) * | 2018-07-04 | 2018-12-28 | 中国科学院地质与地球物理研究所 | 一种高压介质地震波传播模拟方法 |
WO2019119676A1 (zh) * | 2017-12-21 | 2019-06-27 | 西南石油大学 | 一种实验确定酸液对致密碳酸盐岩杨氏模量影响的方法 |
CN111090125A (zh) * | 2019-12-30 | 2020-05-01 | 中国石油大学(北京) | 确定致密砂岩储层的弹性模量的方法和装置 |
CN111366452A (zh) * | 2020-03-26 | 2020-07-03 | 北京科技大学 | 一种自蓄能岩体蓄能水平的测量方法 |
CN111426616A (zh) * | 2020-02-20 | 2020-07-17 | 中国石油大学(北京) | 碳酸盐岩弹性性质与孔隙结构获取方法、装置及存储介质 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1151021A (zh) * | 1995-11-22 | 1997-06-04 | 中国海洋石油总公司海洋石油勘探开发研究中心 | 确定、展示和使用地下岩石弹性模量和密度相对变化的方法 |
WO2003065013A1 (en) * | 2002-01-29 | 2003-08-07 | Baker Hughes Incorporated | Reservoir monitoring in a laminated reservoir using 4-d time lapse data and multicomponent induction data |
WO2007100827A2 (en) * | 2006-02-24 | 2007-09-07 | Saudi Arabian Oil Company | Monte carlo simulation of well logging data |
-
2011
- 2011-05-16 CN CN 201110126399 patent/CN102288986B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1151021A (zh) * | 1995-11-22 | 1997-06-04 | 中国海洋石油总公司海洋石油勘探开发研究中心 | 确定、展示和使用地下岩石弹性模量和密度相对变化的方法 |
WO2003065013A1 (en) * | 2002-01-29 | 2003-08-07 | Baker Hughes Incorporated | Reservoir monitoring in a laminated reservoir using 4-d time lapse data and multicomponent induction data |
WO2007100827A2 (en) * | 2006-02-24 | 2007-09-07 | Saudi Arabian Oil Company | Monte carlo simulation of well logging data |
Non-Patent Citations (3)
Title |
---|
刘春园,魏修成等: "地球物理方法在碳酸盐岩储层预测中的应用综述", 《地球物理学进展》 * |
张金强,曲寿利等: "一种碳酸盐岩储层中流体替换的实现方法", 《石油地球物理勘探》 * |
李景叶,陈小宏等: "高孔低胶结砂岩储层饱和度和压力变化的时移地震AVO响应", 《中国石油大学学报(自然科学版)》 * |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105468886B (zh) * | 2014-09-04 | 2018-07-27 | 中国石油化工股份有限公司 | 基于岩石物理学参数计算地层压力的方法 |
CN105468886A (zh) * | 2014-09-04 | 2016-04-06 | 中国石油化工股份有限公司 | 基于岩石物理学参数计算地层压力的方法 |
CN106251321A (zh) * | 2015-06-12 | 2016-12-21 | 中国石油化工股份有限公司 | 用于获取岩石裂缝信息的方法和装置 |
CN106324672A (zh) * | 2015-07-03 | 2017-01-11 | 中国石油化工股份有限公司 | 一种岩石物理建模方法及*** |
CN106324672B (zh) * | 2015-07-03 | 2018-02-13 | 中国石油化工股份有限公司 | 一种岩石物理建模方法及*** |
US11199483B2 (en) | 2017-12-21 | 2021-12-14 | Southwest Petroleum University | Method for experimentally determining influence of acid liquor on Young's modulus of compact carbonate rock |
WO2019119676A1 (zh) * | 2017-12-21 | 2019-06-27 | 西南石油大学 | 一种实验确定酸液对致密碳酸盐岩杨氏模量影响的方法 |
CN108181192A (zh) * | 2017-12-25 | 2018-06-19 | 万向钱潮(上海)汽车***有限公司 | 一种制动盘的弹性模量获取方法 |
CN109100801A (zh) * | 2018-07-04 | 2018-12-28 | 中国科学院地质与地球物理研究所 | 一种高压介质地震波传播模拟方法 |
CN111090125A (zh) * | 2019-12-30 | 2020-05-01 | 中国石油大学(北京) | 确定致密砂岩储层的弹性模量的方法和装置 |
CN111426616A (zh) * | 2020-02-20 | 2020-07-17 | 中国石油大学(北京) | 碳酸盐岩弹性性质与孔隙结构获取方法、装置及存储介质 |
CN111426616B (zh) * | 2020-02-20 | 2022-05-24 | 中国石油大学(北京) | 碳酸盐岩弹性性质与孔隙结构获取方法、装置及存储介质 |
CN111366452A (zh) * | 2020-03-26 | 2020-07-03 | 北京科技大学 | 一种自蓄能岩体蓄能水平的测量方法 |
CN111366452B (zh) * | 2020-03-26 | 2021-01-29 | 北京科技大学 | 一种自蓄能岩体蓄能水平的测量方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102288986B (zh) | 2013-01-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102288986B (zh) | 一种地震尺度下碳酸盐岩储层弹性模量获取方法 | |
Wei et al. | Geometrical, fractal and hydraulic properties of fractured reservoirs: A mini-review | |
CN110824556B (zh) | 一种非常规致密砂岩储层的岩石物理模型建立方法及应用 | |
CN102426390B (zh) | 一种非均质泥砂岩储层储量确定方法 | |
CN104847340B (zh) | 一种水淹层测井定量评价方法 | |
CN102175832B (zh) | 一种确定典型储层最佳饱和度计算模型的方法 | |
CN102565858B (zh) | 一种多孔介质含水饱和度的计算方法 | |
CN105095631A (zh) | 一种页岩各向异性岩石物理建模方法 | |
Li et al. | Estimation of relative permeability by assisted history matching using the ensemble Kalman filter method | |
CN105464654A (zh) | 确定致密砂岩饱和度指数的方法及装置 | |
CN102156297B (zh) | 基于砂岩油藏叠后地震数据的流体替换方法 | |
CN104730089B (zh) | 并联岩心中流动泡沫层析成像检测装置及其工作方法 | |
WO2012035036A1 (en) | Method of predicting the pressure sensitivity of seismic velocity within reservoir rocks | |
CN106289964B (zh) | 一种泥页岩古抗压强度确定方法 | |
CN104712330A (zh) | 一种测井渗透率解释方法 | |
CN103713327B (zh) | 基于测井及地震资料的页岩气储层最小闭合压力评价方法 | |
CN106285664A (zh) | 基于逾渗网络模拟的双重介质储层岩石含水饱和度计算法 | |
Li et al. | Ensemble-based relative permeability estimation using B-spline model | |
CN105425280A (zh) | 一种矿物模量和孔隙结构预测方法 | |
CN106285652B (zh) | 确定页岩游离气体饱和度的方法 | |
CN102536223B (zh) | 煤层裂隙孔隙度快速计算方法 | |
CN111706318A (zh) | 一种确定低渗储层剩余油分布状况的方法 | |
CN102288993B (zh) | 基于砂岩油藏叠前地震数据的流体替换方法 | |
CN105275459A (zh) | 一种确定页岩地层可动水体积含量的方法 | |
CN105116466A (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 | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20130130 Termination date: 20200516 |
|
CF01 | Termination of patent right due to non-payment of annual fee |