CN104360389B - 一种致密砂岩储层岩石弹性模量计算方法 - Google Patents
一种致密砂岩储层岩石弹性模量计算方法 Download PDFInfo
- Publication number
- CN104360389B CN104360389B CN201410635327.6A CN201410635327A CN104360389B CN 104360389 B CN104360389 B CN 104360389B CN 201410635327 A CN201410635327 A CN 201410635327A CN 104360389 B CN104360389 B CN 104360389B
- Authority
- CN
- China
- Prior art keywords
- rock
- modulus
- high frequency
- modelling quantity
- under
- 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.)
- Active
Links
Abstract
本发明首先从岩石物理参数出发分析了常规砂岩储层与致密砂岩储层各自特点,针对致密砂岩孔隙度低、孔隙连通性差、束缚水饱和度高等主要特征进行研究,基于对包含物形状敏感的岩石物理模型分析岩石中的孔隙形状和孔隙分布特征,考虑流体均匀混合与斑状饱和分布同时存在,结合衰减模型分析低渗透率、高流体粘滞性产生的地球物理数据在不同频带下不匹配的问题,最终将针对致密砂岩储层建立的岩石物理模型应用于储层岩石速度估算。
Description
技术领域
本申请涉及测井领域,更具体地说,涉及一种致密砂岩储层岩石弹性模量计算方法。
背景技术
随着地球物理勘探研究的深入,我们所面临的地下条件变的日益复杂,针对非常规储层的地震勘探研究可以在一定程度上提高勘探的成功率,提高生产效率以及减少勘探开发的成本。
致密储层与常规储层相比,其主要的区别在于致密砂岩储层的孔隙度和渗透率较小,在地震波穿过致密砂岩储层的过程中,岩石孔隙流体的压力不平衡,因此采用Gassmann方程对致密砂岩储层进行岩石物理分析时精度较低。
发明内容
有鉴于此,本申请提供一种致密砂岩储层岩石模量计算方法,以提高对致密砂岩储层的岩石物理分析精度。
为了实现上述目的,现提出的方案如下:
一种致密砂岩储层岩石弹性模量计算方法,包括:
获取测井参数和致密砂岩储层的岩石物性参数;
根据所述测井参数和所述岩石物性参数,结合预先建立的孔隙裂缝衰减模型计算高频下的岩石弹性模量;
根据所述高频下的岩石弹性模量,结合预先建立的自洽模型,计算岩石的高频基质弹性模量;
所述岩石的高频基质弹性模量,结合所述孔隙裂缝衰减模型,计算任意频率下的饱和岩石弹性模量。
优选的,所述岩石的高频基质弹性模量,结合所述孔隙裂缝衰减模型,计算任意频率下的饱和岩石弹性模量,之后还包括:
根据所述任意频率下的饱和岩石弹性模量,结合开始模拟退火寻优的方法,计算岩石的声波速度。
优选的,所述根据所述测井参数和所述岩石物性参数,结合预先建立的孔隙裂缝衰减模型计算高频下的岩石弹性模量包括:
根据所述测井参数与岩石物性参数,利用wood公式和patchy模型计算岩石流体体积模量;
根据所述岩石物性参数与测井参数,利用VRH模型计算岩石基质弹性模量;
根据所述岩石流体体积模量和所述岩石基质弹性模量,结合所述孔隙裂缝衰减模型计算高频下的岩石弹性模量;
其中,所述孔隙裂缝衰减模型的计算公式为:
和
其中,K、λ和μ分别代表岩石基质体积模量、拉梅常数和岩石基质剪切模量,ε为裂缝密度,φP为刚性孔隙度,ω为角频率,r为软孔隙纵横比,Kf为岩石流体体积模量,
τ为时间常数且η为流体粘滞性,α为柔性孔隙半径,υ为岩石基质泊松比,κ为渗透率,ζ为特征矿物颗粒大小,φ为总孔隙度,K1为高频下岩石弹性模量的体积模量,μ1为高频下岩石弹性模量的剪切模量。
优选的,所述高频下的岩石弹性模量,结合预先建立的自洽模型,计算岩石的高频基质弹性模量,包括:
利用自洽模型计算公式:
(1-φp-φs)(Km-Ksc)P*m+φp(Kf-Ksc)P*p+φs(Kf-Ksc)P*s=0
和
(1-φp-φs)(μm-μsc)Q*m+φp(μf-μsc)Q*p+φs(μf-μsc)Q*s=0
计算高频下的岩石基质弹性模量中的体积模量Km和μm;
其中,Ksc为高频下的岩石体积模量,μsc为高频下的岩石剪切模量,φp刚性孔隙度,φs为柔性孔隙度,P和Q为与岩石基质和孔隙结构有关的几何因子。
优选的,利用所述岩石的高频基质弹性模量,结合所述孔隙裂缝衰减模型,计算任意频率下的饱和岩石弹性模量,包括:
利用所述孔隙裂缝衰减模型的计算公式:
计算任意频率下的岩石弹性模量中的岩石体积模量K2和剪切模量μ2,式中:
优选的,所述根据所述任意频率下的岩石模量,结合开始模拟退火寻优的方法,计算岩石的声波速度,包括:
根据计算公式和计算岩石的纵波速度VP和横波速度Vs。
优选的,所述根据所述岩石物性参数与测井参数,利用VRH模型计算岩石基质弹性模量,包括:
利用VRH模型计算公式:
其中,其中MV为模型上限,MR为模型下限。且fi和Mi分别为构成岩石的第i中成分的体积分数和模量;M可以是任何模量,比较常用的有体积模量K、剪切模量μ和杨氏模量E等。
优选的,根据所述任意频率下的岩石模量,结合开始模拟退火寻优的方法,计算岩石的声波速度,之后包括:对所述声波速度进行迭代,寻求最优解。
优选的,所述对所述声波速度进行迭代,寻求最优解,包括:
利用目标函数
其中,为实际测量的纵波速度,为每一次迭代之后计算得到的纵波速度,计算精度满足F<ε时,认为结果收敛,得到最优解。
优选的,所述快速模拟退火法的初始退火温度为:
Tj+1=Tj/(1+a)
其中,a为衰减系数,j为退火过程中的迭代次数;
升温退火温度为:
其中,k0为在上一过程中的迭代次数,β为温度增幅因子与T0呈反比。
经由上述技术方案可知,本发明公开了一种致密砂岩储层岩石弹性模量计算方法。该方法针对致密砂岩储层的孔隙度的特点,利用建立的适用于致密砂岩储层的岩石物理模型,如:孔隙裂缝衰减模型和自洽模型对致密砂岩储层进行分析,得到致密砂岩储层的岩石弹性模量,提高了对致密砂岩储层岩石模量的计算精度。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据提供的附图获得其他的附图。
图1示出了本发明一种致密砂岩储层岩石弹性模量计算方法的一个实施例的流程示意图。
具体实施方式
针对于现有技术中采用Gassmann方程对致密砂岩储层进行岩石物理分析时精度较低的问题,本申请公开了一种致密砂岩储层岩石模量计算方法。
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
参见图1示出了本发明一种致密砂岩储层岩石弹性模量计算方法的一个实施例的流程示意图。
如图1可知,该方法包括:
A:获取测井参数和致密砂岩储层的岩石物性参数。
B:根据所述测井参数和所述岩石物性参数,结合预先建立的孔隙裂缝衰减模型计算高频下的岩石弹性模量。
其中,该步骤具体包括以下步骤:
B1:根据所述测井参数与岩石物性参数,利用wood公式和patchy模型计算岩石流体体积模量。
B2:根据所述岩石物性参数与测井参数,利用VRH模型计算岩石基质弹性模量。
其中,VRH模型的计算公式为:
其中,其中MV为模型上限,MR为模型下限。且fi和Mi分别为构成岩石的第i中成分的体积分数和模量;M可以是任何模量,比较常用的有体积模量K、剪切模量μ和杨氏模量E等。
B3:根据所述岩石流体体积模量和所述岩石基质弹性模量,结合所述孔隙裂缝衰减模型计算高频下的岩石弹性模量;
其中,所述孔隙裂缝衰减模型的计算公式为:
和
其中,K、λ和μ分别代表岩石基质体积模量、拉梅常数和岩石基质剪切模量,ε为裂缝密度,φP为刚性孔隙度,ω为角频率,r为软孔隙纵横比,Kf为岩石流体体积模量,
τ为时间常数且η为流体粘滞性,α为柔性孔隙半径,υ为岩石基质泊松比,κ为渗透率,ζ为特征矿物颗粒大小,φ为总孔隙度,K1为高频下岩石弹性模量的体积模量,μ1为高频下岩石弹性模量的剪切模量。
C:根据所述高频下的岩石弹性模量,结合预先建立的自洽模型,计算岩石的高频基质弹性模量。
该自洽模型的计算公式为:
(1-φp-φs)(Km-Ksc)P*m+φp(Kf-Ksc)P*p+φs(Kf-Ksc)P*s=0
和
(1-φp-φs)(μm-μsc)Q*m+φp(μf-μsc)Q*p+φs(μf-μsc)Q*s=0
计算高频下的岩石基质弹性模量中的体积模量Km和μm;
其中,Ksc为高频下的岩石体积模量,μsc为高频下的岩石剪切模量,φp刚性孔隙度,φs为柔性孔隙度,P和Q为与岩石基质和孔隙结构有关的几何因子。
D:利用所述高频下的饱和岩石基质弹性模量,结合所述孔隙裂缝衰减模型,计算任意频率下的饱和岩石弹性模量。
其中,在本步骤中该孔隙裂缝衰减模型的计算公式为:
计算任意频率下的岩石弹性模量中的岩石体积模量K2和剪切模量μ2,式中:
由以上实施例可知:本发明公开了一种致密砂岩储层岩石弹性模量计算方法。该方法针对致密砂岩储层的孔隙度的特点,利用建立的适用于致密砂岩储层的岩石物理模型,如:孔隙裂缝衰减模型和自洽模型对致密砂岩储层进行分析,得到致密砂岩储层的岩石弹性模量,提高了对致密砂岩储层岩石模量的计算精度。
另外,在实际的勘探中,工程人员需要根据致密砂岩储层的岩石弹性模型计算岩石的声波速度,包括:横波速度和纵波速度。
在本发明的其他实施例中,可采用模拟退火寻优的方法,提高计算声波的准确率。
具体的过程为:
根据计算公式和计算岩石的纵波速度VP和横波速度Vs。
进而,利用目标函数进行寻优。
其中,为实际测量的纵波速度,为每一次迭代之后计算得到的纵波速度,计算精度满足F<ε时,认为结果收敛,得到最优值。
可选的,在本实施例中所述快速模拟退火法的初始退火温度为:
Tj+1=Tj/(1+a)
其中,a为衰减系数,j为退火过程中的迭代次数;
升温退火温度为:
其中,k0为在上一过程中的迭代次数,β为温度增幅因子与T0呈反比。
最后,还需要说明的是,在本文中,诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来,而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
本说明书中各个实施例采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似部分互相参见即可。对所公开的实施例的上述说明,使本领域专业技术人员能够实现或使用本发明。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的,本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下,在其它实施例中实现。因此,本发明将不会被限制于本文所示的这些实施例,而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。
Claims (9)
1.一种致密砂岩储层岩石弹性模量计算方法,其特征在于,包括:
获取测井参数和致密砂岩储层的岩石物性参数;
根据所述测井参数和所述岩石物性参数,结合预先建立的孔隙裂缝衰减模型计算高频下的岩石弹性模量;
根据所述高频下的岩石弹性模量,结合预先建立的自洽模型,计算岩石的高频基质弹性模量;
利用所述岩石的高频基质弹性模量,结合所述孔隙裂缝衰减模型,计算任意频率下的饱和岩石弹性模量;
所述根据所述测井参数和所述岩石物性参数,结合预先建立的孔隙裂缝衰减模型计算高频下的岩石弹性模量包括:
根据所述测井参数与岩石物性参数,利用wood公式和patchy模型计算岩石流体体积模量;
根据所述岩石物性参数与测井参数,利用VRH模型计算岩石基质弹性模量;
根据所述岩石流体体积模量和所述岩石基质弹性模量,结合所述孔隙裂缝衰减模型计算高频下的岩石弹性模量;
其中,所述孔隙裂缝衰减模型的计算公式为:
和
其中,K、λ和μ分别代表岩石基质体积模量、拉梅常数和岩石基质剪切模量,ε为裂缝密度,φP为刚性孔隙度,ω为角频率,r为软孔隙纵横比, Kf为岩石流体体积模量,
τ为时间常数且η为流体粘滞性,α为柔性孔隙半径,υ为岩石基质泊松比,κ为渗透率,ζ为特征矿物颗粒大小,φ为总孔隙度,K1为高频下岩石弹性模量的体积模量,μ1为高频下岩石弹性模量的剪切模量。
2.根据权利要求1所述的方法,其特征在于,所述岩石的高频基质弹性模量,结合所述孔隙裂缝衰减模型,计算任意频率下的饱和岩石弹性模量,之后还包括:
根据所述任意频率下的饱和岩石弹性模量,结合快速模拟退火法,计算岩石的声波速度。
3.根据权利要求1所述的方法,其特征在于,所述高频下的岩石弹性模量,结合预先建立的自洽模型,计算岩石的高频基质弹性模量,包括:
利用自洽模型计算公式:
(1-φp-φs)(Km-Ksc)P*m+φp(Kf-Ksc)P*p+φs(Kf-Ksc)P*s=0
和
(1-φp-φs)(μm-μsc)Q*m+φp(μf-μsc)Q*p+φs(μf-μsc)Q*s=0
计算高频下的岩石基质弹性模量中的体积模量Km和μm;
其中,Ksc为高频下的岩石体积模量,μsc为高频下的岩石剪切模量,φp刚性孔隙度,φs为柔性孔隙度,P和Q为与岩石基质和孔隙结构有关的几何因子。
4.根据权利要求3所述的方法,其特征在于,所述利用所述岩石的高频基质弹性模量,结合所述孔隙裂缝衰减模型,计算任意频率下的饱和岩石弹性模量,包括:
利用所述孔隙裂缝衰减模型的计算公式:
计算任意频率下的岩石弹性模量中的岩石体积模量K2和剪切模量μ2,式中:
5.根据权利要求2所述的方法,其特征在于,所述根据所述任意频率下的饱和岩石弹性模量,结合快速模拟退火法,计算岩石的声波速度,包括:
根据计算公式和计算岩石的纵波速度VP和横波速度Vs;
其中,K2为所述岩石弹性模量中的岩石体积模量,μ2为所述岩石弹性模量中的剪切模量。
6.根据权利要求1所述的方法,其特征在于,所述根据所述岩石物性参数与测井参数,利用VRH模型计算岩石基质弹性模量,包括:
利用VRH模型计算公式:
其中,其中MV为模型上限,MR为模型下限,且 fi和Mi分别为构成岩石的第i中成分的体积分数和模量;M可以是任何模量,比较常用的有体积模量K、剪切模量μ和杨氏模量E等。
7.根据权利要求5所述的方法,其特征在于,根据所述任意频率下的岩石模量,结合快速模拟退火法,计算岩石的声波速度,之后包括:对所述声波速度进行迭代,寻求最优解。
8.根据权利要求7所述的方法,其特征在于,所述对所述声波速度进行迭代,寻求最优解,包括:
利用目标函数
其中,为实际测量的纵波速度,为每一次迭代之后计算得到的纵波速度,计算精度满足F<ε时,认为结果收敛,得到最优解。
9.根据权利要求2所述的方法,其特征在于,所述快速模拟退火法的初始退火温度为:
Tj+1=Tj/(1+a)
其中,a为衰减系数,j为退火过程中的迭代次数;
升温退火温度为:
其中,k0为在上一过程中的迭代次数,β为温度增幅因子与T0呈反比。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410635327.6A CN104360389B (zh) | 2014-11-12 | 2014-11-12 | 一种致密砂岩储层岩石弹性模量计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410635327.6A CN104360389B (zh) | 2014-11-12 | 2014-11-12 | 一种致密砂岩储层岩石弹性模量计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104360389A CN104360389A (zh) | 2015-02-18 |
CN104360389B true CN104360389B (zh) | 2017-04-26 |
Family
ID=52527669
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410635327.6A Active CN104360389B (zh) | 2014-11-12 | 2014-11-12 | 一种致密砂岩储层岩石弹性模量计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104360389B (zh) |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105425280A (zh) * | 2015-11-21 | 2016-03-23 | 西南石油大学 | 一种矿物模量和孔隙结构预测方法 |
CN106353813B (zh) * | 2016-10-25 | 2018-12-25 | 中国海洋石油总公司 | 基于阵列声波测井的流体性质识别方法 |
CN109490963B (zh) * | 2017-09-13 | 2021-04-23 | 中国石油化工股份有限公司 | 裂缝储层岩石物理建模方法及系统 |
CN108664754B (zh) * | 2018-06-22 | 2020-02-14 | 成都理工大学 | 一种由缓冲层密实度计算缓冲层弹性模量的方法 |
CN109001426B (zh) * | 2018-06-22 | 2021-04-27 | 福建省地质工程勘察院 | 缓冲层弹性模量计算方法 |
CN108845109B (zh) * | 2018-06-22 | 2021-01-12 | 成都理工大学 | 一种通过缓冲层均匀度计算缓冲层弹性模量的方法 |
CN109115987B (zh) * | 2018-07-20 | 2021-01-29 | 中国石油天然气股份有限公司 | 一种基于岩石物理模型的流体因子的评价方法及装置 |
CN112485827A (zh) * | 2019-09-12 | 2021-03-12 | 中国石油天然气股份有限公司 | 裂缝孔隙型储层的岩石物理建模方法及装置 |
CN110824556B (zh) * | 2019-10-22 | 2022-03-01 | 中国石油天然气股份有限公司 | 一种非常规致密砂岩储层的岩石物理模型建立方法及应用 |
CN111208565B (zh) * | 2020-03-04 | 2021-08-27 | 中国石油大学(北京) | 基于kt模型的孔缝参数反演方法、装置及存储介质 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8451683B2 (en) * | 2009-04-03 | 2013-05-28 | Exxonmobil Upstream Research Company | Method for determining the fluid/pressure distribution of hydrocarbon reservoirs from 4D seismic data |
US20130229892A1 (en) * | 2010-09-14 | 2013-09-05 | Statoil Asa | Method of predicting the pressure sensitivity of seismic velocity within reservoir rocks |
CN102445708B (zh) * | 2010-10-14 | 2013-10-23 | 中国石油大学(北京) | 三维等效富泥质砂岩速度预测模型 |
CN104101904B (zh) * | 2013-04-03 | 2017-08-18 | 中国石油化工股份有限公司 | 一种快速求取地层横波速度的方法 |
CN103713320B (zh) * | 2013-12-31 | 2015-04-01 | 孙赞东 | 一种富有机质泥页岩岩石物理模型的建立方法 |
-
2014
- 2014-11-12 CN CN201410635327.6A patent/CN104360389B/zh active Active
Non-Patent Citations (2)
Title |
---|
基于Biot相洽理论的致密砂岩弹性参数计算方法;印兴耀,等;《石油物探》;20130930;第52卷(第5期);第445-451页 * |
横波速度计算方法与应用;郭栋,等;《石油地球物理勘探》;20071031;第42卷(第5期);535-538 * |
Also Published As
Publication number | Publication date |
---|---|
CN104360389A (zh) | 2015-02-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104360389B (zh) | 一种致密砂岩储层岩石弹性模量计算方法 | |
Li et al. | A review on hydraulic fracturing of unconventional reservoir | |
Smith et al. | Layered modulus effects on fracture propagation, proppant placement, and fracture modeling | |
CN103258091B (zh) | 非常规油气藏水平井段三维岩体力学模型建立的方法及装置 | |
Wang et al. | Experimental investigations of hydraulic fracture growth complexity in slickwater fracturing treatments | |
CN103713327B (zh) | 基于测井及地震资料的页岩气储层最小闭合压力评价方法 | |
CN105089615A (zh) | 一种基于油藏模型的测井数据历史回归处理方法 | |
Song et al. | Hydraulic fracturing simulation of complex fractures growth in naturally fractured shale gas reservoir | |
CN104914470A (zh) | 碳酸盐岩缝洞储层储量校正方法 | |
CN105527652A (zh) | 一种岩石脆性的测井方法和装置 | |
Tang et al. | Geomechanics evolution integrated with hydraulic fractures, heterogeneity and anisotropy during shale gas depletion | |
Bui et al. | A Coupled Geomechanics-Reservoir Simulation Workflow to Estimate the Optimal Well-Spacing in the Wolfcamp Formation in Lea County | |
CN104316977B (zh) | 一种页岩气储层的岩石模量计算方法 | |
Iverson | Closure stress calculations in anisotropic formations | |
Jiang et al. | Estimation of rock mass deformation modulus using variations in transmissivity and RQD with depth | |
Sehbi et al. | Optimizing fracture stages and completions in horizontal wells in tight gas reservoirs using drainage volume calculations | |
Zhao et al. | Theoretical and numerical analyses of pore-fluid flow patterns around and within inclined large cracks and faults | |
Suarez et al. | Fracturing-to-Production Simulation Approach for Completion Optimization in the Vaca Muerta Shale | |
Iwere et al. | Numerical Simulation of thick, tight fluvial sands | |
Pandey et al. | Vertical growth of hydraulic fractures in layered formations | |
Wallace et al. | Understanding Completion Performance in Niobrara-Codell Reservoirs Through the Use of Innovative Software-Guided Workflows and Models | |
Cherian et al. | Understanding development drivers in horizontal wellbores in the Midland Basin | |
Sharma et al. | Optimizing polymerflood via rate control | |
Abdelmoneim et al. | Modeling hydraulic fractures in finite difference simulators using amalgam LGR (Local Grid Refinement) | |
Swami et al. | A novel approach to history matching and optimization of shale completions and EUR-A case study of Eagle Ford Well |
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 |