CN108983312A - 一种裂缝型储层各向异性岩石模量估算方法 - Google Patents
一种裂缝型储层各向异性岩石模量估算方法 Download PDFInfo
- Publication number
- CN108983312A CN108983312A CN201810517043.5A CN201810517043A CN108983312A CN 108983312 A CN108983312 A CN 108983312A CN 201810517043 A CN201810517043 A CN 201810517043A CN 108983312 A CN108983312 A CN 108983312A
- Authority
- CN
- China
- Prior art keywords
- modulus
- crack
- formula
- rock
- isotropism
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V99/00—Subject matter not provided for in other groups of this subclass
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明涉及一种裂缝型储层各向异性岩石模量估算方法,其特征在于,该方法由纵、横波速度和密度求得背景岩石动态模量,根据岩心实验数据提供的动、静态模量转换关系,可得到裂缝型储层背景各向同性模量;利用基于线性走滑模型的裂缝等效理论,建立速度各向异性参数与裂缝参数之间量化关系,求出裂缝引起的各向异性扰动量,进而可以得到裂缝型储层的完整弹性柔度张量;在给定任意观测角度条件下,能够估算出裂缝型储层各向异性岩石模量。本发明可预测裂缝型储层沿不同观测方位上的各向异性岩石模量,指导储层人工压裂。
Description
技术领域
本发明涉及一种裂缝型储层各向异性岩石模量估算方法,属于勘探物理地球学领域。
背景技术
天然裂缝或人工压裂缝是引起各向异性的主要原因之一,地层受构造应力作用,在脆性区域易于形成高角度裂缝。在油田开发前预知地下岩石模量具有极大的指导价值,大量的岩石物理理论与实验证实地层岩石往往具有各向异性和压力敏感性。
现有计算各向异性岩石模量的方法主要为室内测量方法和等效裂缝介质模型预测方法。室内实验方法采用超声纵横波综合测试仪测量沿0度、45度和90度方向取芯样本的轴向测得的纵、横波速度。需要严格保证测试过程是在高温高压同等边界条件下进行,使岩样在同等状态下完成纵波和横波的测试。然而,对多个方向操作中,取样和测量中都会引入较大不等条件误差。
等效裂缝介质模型是建立含裂缝地层岩石力学模型的基础。作为一种经典的裂缝等效模型,线性滑动理论描述了弹性波在含非理想联结边界介质中的传播。非理想联结边界上仅满足应力连续,而位移不连续,故界面上、下介质存在相对滑动,地下层理面、节理面、共面裂缝等厚度较薄的力学弱面都视作这类界面。Schoenberg提出了单界面的线性滑动模型和长波长近似下周期层状弹性介质的线性滑动模型,此后线性滑动模型得到了诸多改进与推广,可用于描述具有任意对称性的裂缝等效介质,模拟具有旋转对称性的多套裂缝等效性质,建立含单套裂缝VTI介质的等效各向异性模型等,反映裂缝含流体性质。线性滑动模型精度要高于Hudson一阶和二阶近似,但使用弱度参数不能直观反映裂缝物性,Bakulin证明了线性滑动模型与Hudson模型和Thomsen模型之间物理等价性,建立了裂缝弱度参数和裂缝物性参数量化关系。基于岩石物理模型的预测方法,需要预知孔隙纵横比、矿物含量、孔隙度等微观物性参数,实际很难获得完备的输入数据。目前,缺乏利用速度信息和裂缝弱度参数计算任意方向上岩石模量和泊松比等各向异性岩石模量的有效方法。
发明内容
针对上述问题,本发明的目的是提供一种裂缝型储层各向异性岩石模量估算方法。
为实现上述目的,本发明采取以下技术方案:一种裂缝型储层各向异性岩石模量估算方法,其包括以下步骤:
步骤1:向式(1)和式(2)中输入裂缝储层介质的纵波速度VP、横波速度VS和密度ρ数据,根据式(1)和式(2)求出裂缝储层介质中背景岩石的各向同性模量参数:
式中,E0为背景岩石的各向同性杨氏模量;υ0为背景岩石的各向同性泊松比;
步骤2:根据室内岩石力学测试和超声波速度测量,拟合目标裂缝型储层的动态模量和静态模量之间的量化关系,进而根据步骤1得出的各向同性模量参数间接计算出背景岩石的静态各向同性杨氏模量;
步骤3:向式(4)中输入快横波速度VS1和慢横波速度VS2,根据式(4)计算出横波各向异性参数γ(V):
步骤4:利用线性走滑理论,建立完整的裂缝型储层弹性柔度张量:
S=S0+△S (5)
式中,S0为背景岩石柔度张量;△S为裂缝扰动引起的附加柔度张量;
其中,KT和KN分别为裂缝切向柔度和法向柔度,二者与横波各向异性参数具有以下关系:
利用得到背景岩石的静态各向同性杨氏模量、各向同性泊松比和横波各向异性参数,根据式(8)求出两个裂缝柔度参数,然后根据式(5)~式(7)求出完整的裂缝型储层弹性柔度张量;
步骤5:在获得完整的裂缝型储层弹性柔度张量的基础上,利用Bond变换求出沿任意方位上的裂缝型储层各向异性岩石模量:
经Bond变换后的弹性刚度张量表示为:
其中,方位转换矩阵M为:
式中,为观测方位角;
裂缝型储层的各向异性岩石模量和泊松比:
式中,为裂缝型储层的各向异性岩石模量;和为裂缝型储层的泊松比。
本发明由于采取以上技术方案,其具有以下优点:1、本发明可估算出沿任意观测角度的各向异性岩石模量参数,研究裂缝对储层岩石力学性质的影响,有助于优化储层压裂等工程应用。2、本发明仅需要0度和90度两个方向的横波速度提供速度各向异性,极大减少了速度测量次数。3、本发明对输入数据要求不高,既可以用于计算岩石物理实验数据,也适用于测井数据,叠前地震反演方法能够得到的速度参数和裂缝参数,因此可以基于本发明可得到点、线、三维体的各向异性岩石模量。4、本发明可有效预测裂缝型储层的杨氏模量、泊松比等各向异性岩石模量,是一种简单、实用、高效的方法。
附图说明
图1是裂缝型储层沿不同方位上的单轴应力-应变示意图;
图2是本发明的流程图;
图3(a)和(b)是根据本发明所估算的各向异性杨氏模量和各向同性岩石模量的结果图;
图4(a)和(b)是根据本发明所估算的沿x1-x2方向的各向异性泊松比和各向同性泊松比的结果图;
图5(a)和(b)是根据本发明所估算的沿x1-x3方向的各向异性泊松比和各向同性泊松比的结果图。
具体实施方式
下面结合实施例对本发明进行详细的描述。
一种裂缝型储层各向异性岩石模量估算方法,该方法可预测裂缝型储层沿不同观测方位上的各向异性岩石模量,指导储层人工压裂。本发明的基本思路是由纵、横波速度和密度求得背景岩石动态模量,根据岩心实验数据提供的动、静态模量转换关系,可得到裂缝型储层背景各向同性模量;利用基于线性走滑模型的裂缝等效理论,建立速度各向异性参数与裂缝参数之间量化关系,求出裂缝引起的各向异性扰动量,进而可以得到裂缝型储层的完整弹性柔度张量;在给定任意观测角度条件下,能够估算出裂缝型储层各向异性岩石模量(如图1所示)。
如图2所示,本发明提出的裂缝型储层各向异性岩石模量估算方法,其包括以下步骤:
步骤1:向式(1)和式(2)中输入裂缝储层介质的纵波速度VP、横波速度VS和密度ρ数据,根据式(1)和式(2)求出裂缝储层介质中背景岩石的各向同性模量参数:
式中,E0为背景岩石的各向同性杨氏模量;υ0为背景岩石的各向同性泊松比。
步骤2:根据室内岩石力学测试和超声波速度测量,拟合目标裂缝型储层的动态模量和静态模量之间的量化关系,进而根据步骤1得出的各向同性模量参数间接计算出背景岩石的静态各向同性杨氏模量。在不具备室内试验的条件下,也可以直接利用式(3)计算背景岩石的静态各向同性杨氏模量:
Estatic=(Edynamic/3.3674)2.042 (3)
式中,Estatic为静态各向同性杨氏模量;Edynamic为动态各向同性杨氏模量。
步骤3:向式(4)中输入快横波速度VS1和慢横波速度VS2,根据式(4)计算出横波各向异性参数γ(V):
步骤4:利用线性走滑理论,建立完整的裂缝型储层弹性柔度张量:
S=S0+△S (5)
式中,S0为背景岩石柔度张量;△S为裂缝扰动引起的附加柔度张量。
其中,KT和KN分别为裂缝切向柔度和法向柔度,KT能够反映裂缝密度,而KN可以反映裂缝所含流体成分,二者与横波各向异性参数具有以下关系:
利用得到背景岩石的静态各向同性杨氏模量、各向同性泊松比和横波各向异性参数,根据式(8)可求出两个裂缝柔度参数,然后根据式(5)~式(7)可求出完整的裂缝型储层弹性柔度张量。
步骤5:在获得完整的裂缝型储层弹性柔度张量的基础上,利用Bond变换可求出沿任意方位上的裂缝型储层各向异性岩石模量:
经Bond变换后的弹性刚度张量可表示为:
其中,方位转换矩阵M为:
式中,为观测方位角。
岩石在单向受压(或受拉)时,杨氏模量是轴向应力和轴向应变的比值,泊松比是垂直轴向应变与轴向应变绝对值之比。根据上述定义,可以得到裂缝型储层的各向异性岩石模量和泊松比:
式中,为裂缝型储层的各向异性岩石模量;和为裂缝型储层的泊松比。
式(11)~式(13)给出了沿任意方位上观测得到的杨氏模量和泊松比,当裂缝柔度参数KN和KT为零时,各向异性岩石模量退化为各向同性值。
图3(a)和图3(b)分别为横波各向异性参数等于-0.05和-0.1时根据本发明所估算的各向异性杨氏模量(实线)和各向同性岩石模量(虚线)。图4(a)和图4(b)分别为横波各向异性参数等于-0.05和-0.1时根据本发明所估算的沿x1-x2方向的各向异性泊松比(实线)和各向同性泊松比(虚线)。图5(a)和图5(b)分别为横波各向异性参数等于-0.05和-0.1时根据本发明所估算的沿x1-x3方向的各向异性泊松比(实线)和各向同性泊松比(虚线)。可以看到,各向异性杨氏模量和泊松比沿方位展现出了显著的椭圆各向异性特征,随着横波各向异性参数的增大,椭圆率增大,即沿裂缝法向方向的岩石强度降低。
上述各实施例仅用于对本发明的目的、技术方案和有益效果进行示例性描述,并不局限于上述具体实施方式,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (1)
1.一种裂缝型储层各向异性岩石模量估算方法,其包括以下步骤:
步骤1:向式(1)和式(2)中输入裂缝储层介质的纵波速度VP、横波速度VS和密度ρ数据,根据式(1)和式(2)求出裂缝储层介质中背景岩石的各向同性模量参数:
式中,E0为背景岩石的各向同性杨氏模量;υ0为背景岩石的各向同性泊松比;
步骤2:根据室内岩石力学测试和超声波速度测量,拟合目标裂缝型储层的动态模量和静态模量之间的量化关系,进而根据步骤1得出的各向同性模量参数间接计算出背景岩石的静态各向同性杨氏模量;
步骤3:向式(4)中输入快横波速度VS1和慢横波速度VS2,根据式(4)计算出横波各向异性参数γ(V):
步骤4:利用线性走滑理论,建立完整的裂缝型储层弹性柔度张量:
S=S0+△S (5)
式中,S0为背景岩石柔度张量;△S为裂缝扰动引起的附加柔度张量;
其中,KT和KN分别为裂缝切向柔度和法向柔度,二者与横波各向异性参数具有以下关系:
利用得到背景岩石的静态各向同性杨氏模量、各向同性泊松比和横波各向异性参数,根据式(8)求出两个裂缝柔度参数,然后根据式(5)~式(7)求出完整的裂缝型储层弹性柔度张量;
步骤5:在获得完整的裂缝型储层弹性柔度张量的基础上,利用Bond变换求出沿任意方位上的裂缝型储层各向异性岩石模量:
经Bond变换后的弹性刚度张量表示为:
其中,方位转换矩阵M为:
式中,为观测方位角;
裂缝型储层的各向异性岩石模量和泊松比:
式中,为裂缝型储层的各向异性岩石模量;和为裂缝型储层的泊松比。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810517043.5A CN108983312B (zh) | 2018-05-25 | 2018-05-25 | 一种裂缝型储层各向异性岩石模量估算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810517043.5A CN108983312B (zh) | 2018-05-25 | 2018-05-25 | 一种裂缝型储层各向异性岩石模量估算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108983312A true CN108983312A (zh) | 2018-12-11 |
CN108983312B CN108983312B (zh) | 2019-09-24 |
Family
ID=64542121
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810517043.5A Active CN108983312B (zh) | 2018-05-25 | 2018-05-25 | 一种裂缝型储层各向异性岩石模量估算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108983312B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109596713A (zh) * | 2018-12-29 | 2019-04-09 | 西南石油大学 | 一种基于分形分维的煤岩超声波衰减系数评价方法 |
CN110909487A (zh) * | 2019-12-06 | 2020-03-24 | 西南石油大学 | 一种精细各向异性页岩岩石物理模型的建立方法 |
CN110954954A (zh) * | 2019-12-23 | 2020-04-03 | 西南石油大学 | 一种估计裂缝间距对地震各向异性影响的方法 |
CN111911146A (zh) * | 2020-08-17 | 2020-11-10 | 中国石油大学(华东) | 一种基于裂缝交汇图版的储层压裂效果评价方法 |
CN113791457A (zh) * | 2021-09-08 | 2021-12-14 | 中国海洋石油集团有限公司 | 一种天然气水合物储层岩石骨架模量的计算方法及装置 |
CN115758851A (zh) * | 2022-11-28 | 2023-03-07 | 中国海洋石油集团有限公司 | 一种含天然裂缝地层裂缝多尺度支撑剂的选择方法 |
CN116165054A (zh) * | 2023-04-21 | 2023-05-26 | 煤炭科学研究总院有限公司 | 岩体力学参数获取方法、装置和电子设备 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101553742A (zh) * | 2006-09-12 | 2009-10-07 | 普拉德研究及开发股份有限公司 | 使用图像和声波测井图的组合区分天然裂缝导致的声波各向异性和应力导致的声波各向异性 |
CN104406849A (zh) * | 2014-11-21 | 2015-03-11 | 中国石油天然气股份有限公司 | 一种储层岩石脆性的预测方法及装置 |
CN104573150A (zh) * | 2013-10-18 | 2015-04-29 | 中国石油化工股份有限公司 | 一种针对富含有机质页岩的各向异性岩石物理建模方法 |
US20160370499A1 (en) * | 2013-07-08 | 2016-12-22 | Schlumberger Technology Corporation | Determining Geomechanics Completion Quality |
CN107203005A (zh) * | 2016-03-18 | 2017-09-26 | 中国石油化工股份有限公司 | 一种定量化计算裂缝描述参数的方法 |
-
2018
- 2018-05-25 CN CN201810517043.5A patent/CN108983312B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101553742A (zh) * | 2006-09-12 | 2009-10-07 | 普拉德研究及开发股份有限公司 | 使用图像和声波测井图的组合区分天然裂缝导致的声波各向异性和应力导致的声波各向异性 |
US20160370499A1 (en) * | 2013-07-08 | 2016-12-22 | Schlumberger Technology Corporation | Determining Geomechanics Completion Quality |
CN104573150A (zh) * | 2013-10-18 | 2015-04-29 | 中国石油化工股份有限公司 | 一种针对富含有机质页岩的各向异性岩石物理建模方法 |
CN104406849A (zh) * | 2014-11-21 | 2015-03-11 | 中国石油天然气股份有限公司 | 一种储层岩石脆性的预测方法及装置 |
CN107203005A (zh) * | 2016-03-18 | 2017-09-26 | 中国石油化工股份有限公司 | 一种定量化计算裂缝描述参数的方法 |
Non-Patent Citations (1)
Title |
---|
吴国忱 等: "横观各向同性(TI)岩石弹性模量一阶扰动近似", 《中国科学:地球科学》 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109596713A (zh) * | 2018-12-29 | 2019-04-09 | 西南石油大学 | 一种基于分形分维的煤岩超声波衰减系数评价方法 |
CN110909487A (zh) * | 2019-12-06 | 2020-03-24 | 西南石油大学 | 一种精细各向异性页岩岩石物理模型的建立方法 |
CN110954954A (zh) * | 2019-12-23 | 2020-04-03 | 西南石油大学 | 一种估计裂缝间距对地震各向异性影响的方法 |
CN111911146A (zh) * | 2020-08-17 | 2020-11-10 | 中国石油大学(华东) | 一种基于裂缝交汇图版的储层压裂效果评价方法 |
CN113791457A (zh) * | 2021-09-08 | 2021-12-14 | 中国海洋石油集团有限公司 | 一种天然气水合物储层岩石骨架模量的计算方法及装置 |
CN113791457B (zh) * | 2021-09-08 | 2022-11-15 | 中国海洋石油集团有限公司 | 一种天然气水合物储层岩石骨架模量的计算方法及装置 |
CN115758851A (zh) * | 2022-11-28 | 2023-03-07 | 中国海洋石油集团有限公司 | 一种含天然裂缝地层裂缝多尺度支撑剂的选择方法 |
CN115758851B (zh) * | 2022-11-28 | 2024-01-05 | 中国海洋石油集团有限公司 | 一种含天然裂缝地层裂缝多尺度支撑剂的选择方法 |
CN116165054A (zh) * | 2023-04-21 | 2023-05-26 | 煤炭科学研究总院有限公司 | 岩体力学参数获取方法、装置和电子设备 |
Also Published As
Publication number | Publication date |
---|---|
CN108983312B (zh) | 2019-09-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108983312B (zh) | 一种裂缝型储层各向异性岩石模量估算方法 | |
Herwanger | Seismic geomechanics: how to build and calibrate geomechanical models using 3D and 4D seismic data | |
CN105068117B (zh) | 用于裂缝性介质的avo反演方法、装置和设备 | |
Qiao et al. | Determination of Biot's effective-stress coefficient for permeability of Nikanassin sandstone | |
CA2580570A1 (en) | Integrated anisotropic rock physics model | |
CN104181585A (zh) | 一种地球物理勘探横波估算方法和系统 | |
CN113075728B (zh) | 一种建立致密砂岩多尺度三维岩石物理图板的方法 | |
CN104570065B (zh) | 一种利用地震波阻抗定量反演孔隙度的方法 | |
CN104484573B (zh) | 一种确定地层刚性系数的方法 | |
Sun et al. | Fractured reservoir modeling by discrete fracture network and seismic modeling in the Tarim Basin, China | |
Deng et al. | Cross-bedding related anisotropy and its interplay with various boundary conditions in the formation and orientation of joints in an aeolian sandstone | |
CN110471129A (zh) | 一种深层页岩高温高压下的各向异性岩石物理建模方法 | |
Far et al. | Geomechanics of orthorhombic media | |
CN113589385B (zh) | 一种基于地震散射波场分析的储层特征反演方法 | |
Korneev et al. | Elastic properties of fluid-saturated fractured rock | |
CN102809761A (zh) | 一种获取煤岩的速度特征的控制方法 | |
WANG et al. | Multi‐Azimuth Three‐Component Surface Seismic Modeling in Cracked Monoclinic Media | |
Jiang et al. | The constructing of pore structure factor in carbonate rocks and the inversion of reservoir parameters | |
Nourifard et al. | The effects of stress and fluid on the anisotropy of reservoir rock: case study of a sandstone from the harvey 3 CCS site, Western Australia | |
Arora et al. | The Effect of Composite Material on Rayleigh Wave at Free Surface of Composite Matrix Saturated by Fluids | |
Gibson Jr et al. | An effective medium model for the stress-dependence of anisotropic seismic velocities in fractured rock | |
Bedayat et al. | The equivalent inclusion method for poroelasticity problems | |
Wang et al. | Multi-azimuth three-component surface seismic modeling for viscoelastic cracked monoclinic media | |
CN118131311A (zh) | 一种正交各向异性建模方法及装置 | |
Li et al. | Effects of rock anisotropy on integrated modelling |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |