CN103576195B - 一种随压力变化的裂隙介质横波速度预测方法 - Google Patents
一种随压力变化的裂隙介质横波速度预测方法 Download PDFInfo
- Publication number
- CN103576195B CN103576195B CN201310518774.9A CN201310518774A CN103576195B CN 103576195 B CN103576195 B CN 103576195B CN 201310518774 A CN201310518774 A CN 201310518774A CN 103576195 B CN103576195 B CN 103576195B
- Authority
- CN
- China
- Prior art keywords
- wave velocity
- rock
- medium
- pressure
- phi
- 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
- 238000000034 method Methods 0.000 title claims abstract description 46
- 230000001419 dependent effect Effects 0.000 title claims abstract description 7
- 239000011435 rock Substances 0.000 claims abstract description 148
- 239000011148 porous material Substances 0.000 claims abstract description 99
- 239000012530 fluid Substances 0.000 claims abstract description 27
- 230000008859 change Effects 0.000 claims description 21
- 238000005259 measurement Methods 0.000 claims description 20
- 239000002245 particle Substances 0.000 claims description 14
- 230000008569 process Effects 0.000 claims description 13
- 238000012360 testing method Methods 0.000 claims description 13
- 238000012937 correction Methods 0.000 claims description 6
- 230000009466 transformation Effects 0.000 claims description 5
- 230000014509 gene expression Effects 0.000 claims description 2
- 230000001131 transforming effect Effects 0.000 claims description 2
- 229920006395 saturated elastomer Polymers 0.000 description 8
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 7
- 238000002347 injection Methods 0.000 description 6
- 239000007924 injection Substances 0.000 description 6
- 238000004458 analytical method Methods 0.000 description 4
- 238000007596 consolidation process Methods 0.000 description 4
- 230000006870 function Effects 0.000 description 4
- 238000004519 manufacturing process Methods 0.000 description 3
- 238000013459 approach Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 1
- 235000019738 Limestone Nutrition 0.000 description 1
- 229910052925 anhydrite Inorganic materials 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000012267 brine Substances 0.000 description 1
- OSGAYBCDTDRGGQ-UHFFFAOYSA-L calcium sulfate Chemical compound [Ca+2].[O-]S([O-])(=O)=O OSGAYBCDTDRGGQ-UHFFFAOYSA-L 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 229910052799 carbon Inorganic materials 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 239000006028 limestone Substances 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000013031 physical testing Methods 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 230000009919 sequestration Effects 0.000 description 1
- HPALAKNZSZLMCH-UHFFFAOYSA-M sodium;chloride;hydrate Chemical compound O.[Na+].[Cl-] HPALAKNZSZLMCH-UHFFFAOYSA-M 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
Landscapes
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
- Geophysics And Detection Of Objects (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
一种随压力变化的裂隙介质横波速度预测方法,采集测井数据,建立随压力变化的孔隙介质横波速度预测模型;将校正过的各个参数代入随压力变化的孔隙介质横波速度预测模型得到孔隙介质的与压力相关的纵横波速度;根据孔隙介质下的与压力相关的纵波速度、横波速度得到裂隙介质对称方向的与压力相关的纵波速度、横波速度并代入Thomson公式中,得到平行于裂隙方向的与压力相关的横波速度和垂直于裂隙方向的与压力相关的横波速度。本发明既考虑了压力对于介质横波速度的影响,也考虑到了储层中含有裂隙,更加符合实际。根据本发明得到的裂隙介质横波速度,可以更好的评估油气藏岩石属性,包括孔隙度、密度、岩性与流体含量,可以建立AVO模。
Description
技术领域
本发明属于地震勘探中岩石物理领域,具体涉及一种随压力变化的裂隙介质横波速度预测方法。
背景技术
人们通常用纵横波速度比为常数来代替横波速度,但是对于不同的介质纵横波速度比应该是不同的,这样的假设不合理。因此很多横波速度的预测方法随之而出,有经验公式,也有基于岩石物理理论的。
比较著名的经验公式包括Castagna等人在1985年提出的著名的泥岩线为Vs=0.862Vp-1.172。Gardner于1974年给出了不同的岩性之间的速度与密度的关系,其中它的平均变换式为ρ=0.23V0.25,这个平均变换式是对所有岩性的速度与密度关系的最佳拟合,它适合于所有岩性,不仅仅适用用某种岩性。而Castagna于1993年又对Gardner的公式进行了扩充,得到了不同岩性的速度与密度之间的关系:对于砂岩有对于页岩有对于石灰岩有对于白云岩有对硬石膏则是Wyllie等人在1958年和1963年陆续提出了充满盐水的孔隙介质的孔隙度与速度之间的经验关系:1/V=(1-φ)/Vma+φ/Vf1,其中V为岩石的整体速度,Vma为岩石骨架的速度,Vf1为孔隙流体的速度,φ为孔隙度。当某些测井曲线缺失的时候,或者地震振幅异常都可以应用这些期限进行质量监控,但是这些公式对于岩性的依赖很强,而且依赖于局部条件,而Mavko等人多次在他们的书中提到:“这些关系公式都是经验公式,因而严格来说它们只能用在当时研究的岩石上”,因此这些经验公式不具有普遍性。
随着岩石物理理论的完善,基于岩石物理理论的横波速度预测方法逐渐成为研究的主流。如Greenberg和Castagna于1992年利用Biot-GassmannTheory(BGT)进行横波速度的预测,也就是假设在纵横波速度之间存在一个稳固的关系同时假设固体岩石成分之间的混合定律是线性的。Xu和White将Kuster和于1974年建立的理论与差分有效介质理论结合,进行岩石弹性模量的计算,具体表现为利用孔隙纵横比来表征砂泥成分之间的关系,然后利用Gassmann方程来预测横波速度,该方法后来又被引申到预测裂隙介质的横波速度预测。Nolen-Hoeksema与WangZhijing于1996年根据实验室测得的干岩石的横波速度,利用Gassmann方程计算出干岩石的体变模量与切变模量,进而用到流体饱和岩石的横波速度预测中。在2006年的时候Lee提出的用固结参数联系基质弹性模量与骨架弹性模量之间的关系,通过实测纵波速度与预测纵波速度的比较得出固结参数,然后利用固结参数计算横波速度。2008年孙福利等人利用实际的数据对Lee的方法进行了验证,并且提出了固结参数的取值范围。
但是现有的这些方法没有考虑裂隙介质中横波速度分裂问题,更没有考虑到压力对裂隙介质属性的影响。而在实际生产中,储层存在裂隙,或者随着油田开发过程中的注水压裂等而产生裂隙,这些裂隙显然是受储层压力而发生变化的。在油田注水开发及CCS(碳捕捉与封存)技术中,注入井点的压力很高,随着水或者CO2的高压注入,注入井附近可能产生新裂缝或者将原有裂缝撑大。而注CO2驱油时,生产井点则由于控制孔隙压力而可能导致裂隙的闭合。此时,注入井与生产井点储层内压力诱导下的快慢横波速度会发生变化。如何预测出不同孔隙压力下裂隙介质的快慢波速度,是利用四维地震监测CO2在地下封存安全性的基础。无论是利用四维地震解释裂隙介质还是四维方位AVO反演、各向异性反演,以及四维地震转换波解释,都需要利用随孔隙压力变化的快慢横波速度才能进行。而同时考虑了压力变化和裂隙的快慢横波速度的地震解释也才更加符合实际。
发明内容
本发明的目的在于提供一种随压力变化的孔隙加裂隙介质的横波速度预测方法,预测得到的横波速度更加符合实际。
为实现上述目的,本发明采用如下的技术方案:
本发明包括以下步骤:
1)采集测井数据:采集总孔隙度φz、岩石的体积密度ρ以及实际纵波速度Vp、流体饱和度;采集岩石骨架的体变模量Kma,岩石骨架的切变模量μma,差异压力p,测量配位数Cp',变形之前接触区域的半径a与颗粒的半径R;
分别对总孔隙度φz、岩石的体积密度ρ进行校正得到孔隙介质的孔隙孔隙度φp、孔隙介质的体积密度ρp,利用流体饱和度得到混合流体的体变模量Kf;计算岩石骨架的泊松比νma;
对测量配位数Cp'进行线性拟合并且加权,得到含有加权系数W的配位数Cp,然后利用Digby、Mindlin公式以及Gassmann方程建立随压力变化的孔隙介质横波速度预测模型;
2)对岩心进行岩石物理测试,得到岩心的纵波速度、岩心的横波速度、岩心的密度以及岩心的孔隙度,然后利用随压力变化的孔隙介质横波速度预测模型对岩心的纵波速度、岩心的横波速度、岩心的密度以及岩心的孔隙度进行计算得到岩心的加权系数Wcore,然后利用岩心的加权系数Wcore对加权系数W进行校正,得到校正后的各向同性介质的加权系数Wisotropy,如公式(14)所示:
其中,Wlog为岩心对应层位的加权系数;
3)将校正后的各向同性介质的加权系数Wisotropy、孔隙介质的孔隙孔隙度φp以及孔隙介质下的密度ρp,代入Gassmann方程的变形公式,得到孔隙介质下随压力变化的纵波速度α以及横波速度β;
4)根据孔隙介质下随压力变化的纵波速度α以及横波速度β,得到对称方向的弹性介质的纵波速度α0、横波速度β0;将对称方向的弹性介质的纵波速度α0、横波速度β0代入公式(15)-(17)中,得到平行于裂隙方向的横波速度Vs1和垂直于裂隙方向的横波速度Vs2;
其中θ为波前与对称轴的夹角,ε、γ和δ为各向异性参数,Vp(θ)为预测纵波速度;
5)根据步骤4)得到的平行于裂隙方向的横波速度Vs1和垂直于裂隙方向的横波速度Vs2,建立各向异性介质的AVO模型,预测油气藏岩石属性。
所述步骤1)中若取岩心的井与进行采集测井数据的井是两口井,则标定岩心在进行采集数据的井中的层位。
所述步骤1)中随压力变化的孔隙介质横波速度预测模型建立的具体过程为:
a)对测量配位数Cp'进行线性拟合,得到测量配位数Cp'与孔隙度φ的关系为:Cp'=11.759e1-φ-12.748(1)
对测量配位数Cp'进行进行加权,得到配位数Cp,如公式(2)所示,
Cp=W*Cp'(2)
b)在均匀的孔隙介质的条件下,利用Gassmann方程的变形进行纵横波速度的预测;公式(3)-(5)为Gassmann方程的变形为:
μsat=μdry(4)
其中分别为预测的纵波速度、横波速度,μsat为孔隙介质的切变模量,Kdry为干岩石的体变模量、μdry为干岩石的切变模量,Kma为岩石骨架的体变模量;Kf为混合流体的体变模量,ρ为岩石的体积密度;
将配位数Cp代入Digby公式中得到公式(7):
将配位数Cp代入Mindlin公式中得到公式(9):
其中,νma为岩石骨架的泊松比,μma为岩石骨架的切变模量;φ为孔隙度;Cp'为测量配位数;a为岩石颗粒变形之前接触区域的半径,b为岩石颗粒变形之后接触区域的半径,R为岩石颗粒的半径;
将干岩石的体变模量Kdry与干岩石的切变模量μdry代入Gassmann方程的变形公式中计算预测纵波速度利用预测的纵波速度值等于实际测量的纵波速度值,得到加权系数W;
c)将加权系数W代入公式(2),得到配位数Cp,再把配位数Cp值代入计算随压力变化的干岩石切变模量的Mindlin公式中得到干岩石的切变模量μdry,再将干岩石的切变模量μdry代入Gassmann方程计算得到预测横波速度即建立起了随压力变化的孔隙介质横波速度预测模型。
所述步骤b)中Digby公式为:
所述步骤b)和步骤c)中Mindlin公式为:
所述步骤b)中
其中p为差异压力。
所述步骤3)中将总孔隙度φz校正为孔隙介质下的孔隙孔隙度φp的具体过程为:将总孔隙度φz减去裂隙孔隙度φc得到孔隙孔隙度φp。
所述步骤1)中对岩石的体积密度ρ进行校正,得到孔隙介质情况下的孔隙介质的密度ρp的具体过程如下:
首先利用公式ρd=ρ*(1-φz)计算得到干岩石的体积密度,其中ρd为干岩石的体积密度,然后利用公式ρp=ρd/(1-φp)得到孔隙介质的体积密度ρp。
所述步骤4)中对称方向的弹性介质的纵波速度α0、横波速度β0是通过以下过程得到的:
根据公式(18)和(19),计算得到对称方向的弹性介质的纵波速度α0、横波速度β0;
其中α、β分别为孔隙介质下随压力变化的纵波速度、横波速度;ν为孔隙介质的泊松比,ε以及γ为各向异性参数。
所述步骤4)中各向异性参数ε、γ和δ的表达式分别为:
公式(20)和(21)中的
公式(20)中Dci为流体影响因子,在低频和高频声波测量时分别如下:
公式(23)中
其中,Kma为岩石骨架的体变模量,Kf为混合流体的体变模量,ν为孔隙介质的泊松比,E为杨氏模量,νd和Ed为相应的各向同性孔隙干岩石的泊松比和杨氏模量,νma则为岩石骨架的泊松比,φz为总孔隙度,Dci(lo)为低频声波测量时的流体影响因子,Dci(mh)为高频声波测量时的流体影响因子。
本发明的有益效果是:本发明针对现有技术中利用孔隙介质横波速度预测模型预测所有储层横波速度与实际情况不符的问题,本发明考虑了储层中发育裂隙的情况,以建立的随压力变化孔隙介质横波速度预测模型为基础,利用岩石物理测试资料对预测的孔隙介质纵横波速度等参数进行校正,然后利用Thomson裂隙介质理论进行横波速度的预测,这样更加符合存在裂隙的储层。由于横波速度对压力是敏感的,本发明既考虑了压力对于介质横波速度的影响,同时也考虑到了储层中含有裂隙,更加符合实际。根据本发明得到的裂隙介质横波速度,可以建立AVO模型,进行转换波分析,通过AVO分析,可以更好地评估油气藏岩石属性,包括孔隙度、密度、岩性与流体含量。
附图说明
图1为C语言编程流程图。
图2为岩心为干岩石进行校正,得到的快横波速度与实测快横波速度的对比,其中实线为实测值,虚线为预测值。
图3为岩心为干岩石进行校正,得到的慢横波速度与实测慢横波速度的对比,其中实线为实测值,虚线为预测值。
图4为岩心为油水饱和岩石进行校正,得到的快横波速度与实测快横波速度的对比,其中实线为实测值,虚线为预测值。
图5为岩心为油水饱和岩石进行校正后得到的慢横波速度与实测慢横波速度的对比,其中实线为实测值,虚线为预测值。
具体实施方式
下面结合附图对本发明做详细描述。本发明中的*表示乘号。本发明利用采集数据,得到测井资料,进行储层横波速度的预测,具体涉及一种利用Thomson裂隙介质理论结合随压力变化的孔隙介质模型以及岩石物理测试资料进行含裂隙的储层横波速度预测的方法。
本发明具体包括以下步骤:
1)本发明利用考虑压力变化的Digby、Mindlin公式以及Gassmann方程的变形建立的随压力变化孔隙介质横波速度预测模型来预测储层介质的配位数,然后利用岩石物理测试资料来校正该配位数,得到更为准确的孔隙介质的配位数,以及纵横波速度。以孔隙介质的纵横波速度为基础,利用Thomson裂隙介质理论进行横波速度的预测。
需要说明一下孔隙度的问题,在裂隙介质中总孔隙度为裂隙孔隙度和孔隙孔隙度的和;而下面建立的是孔隙介质模型,介质为孔隙介质,所以只有一个孔隙度;建立的随压力变化孔隙介质横波速度预测模型的具体过程为:
a)采集测井数据:采集总孔隙度φz、岩石的体积密度ρ以及实际纵波速度Vp、饱和度;采集岩石骨架的体变模量Kma,岩石骨架的切变模量μma,差异压力p,测量配位数Cp',岩石颗粒变形之前接触区域的半径a与岩石颗粒的半径R;根据公式计算得到岩石骨架的泊松比νma,其中Kma为岩石骨架的体变模量,μma为岩石骨架的切变模量。
根据测量配位数Cp'与e1-φ成正比的结论,进行线性拟合得到测量配位数与孔隙度的关系为:Cp'=11.759e1-φ-12.748(1)
其中Cp'为测量配位数,φ为孔隙度;
在实际情况中,对测量配位数Cp'进行加权,得到配位数Cp,如公式(2)所示,其中W为加权系数:
Cp=W*Cp'(2);
b)假设介质为均匀的孔隙介质,利用Gassmann方程的变形进行纵横波速度的预测;公式(3)-(5)为Gassmann方程的变形:
μsat=μdry(4)
其中与分别为预测的纵波速度、横波速度,μsat为孔隙介质的切变模量,Kdry与μdry分别为干岩石的体变模量和切变模量,Kma为岩石骨架的体变模量;φ为岩石的孔隙度;Kf为混合流体的体变模量,ρ为岩石的体积密度;其中岩石骨架的体变模量Kma、孔隙度φ以及岩石的体积密度ρ均从采集的测井数据资料中获得,混合流体的体变模量Kf由流体的饱和度得到;
计算随压力变化的压实紧密的干岩石的体变模量的Digby公式为:
将公式(2)中的配位数Cp公式代入计算随压力变化的压实紧密的干岩石的体变模量的Digby公式得到公式(7):
计算随压力变化的干岩石切变模量的Mindlin公式为:
将公式(2)中的的配位数Cp公式计算随压力变化的干岩石切变模量的Mindlin公式中得到公式(9):
其中,
φ为孔隙度;Cp'为测量配位数;p为差异压力;a为岩石颗粒变形之前接触区域的半径,b为岩石颗粒变形之后接触区域的半径,R为岩石颗粒的半径;差异压力p的变化会体现在中,即公式(11),无论是在公式(11)中求得差异压力p值还是在Digby、Mindlin公式中利用差异压力p值,在实际计算中是作为一个整体的;
将Kdry与μdry代入Gassmann方程的变形公式中计算预测纵波速度由于实际采集的数据中有纵波速度,预测的纵波速度值要接近实际测量的纵波速度值,因此利用公式(13)即实际测量纵波速度值Vpmeasured减去预测纵波速度值Vpsat(W)等于零(从理论上来说预测的横波速度要无限接近实际的横波速度,而在实际上,无限接近无法实际操作,因此假设实际测量纵波速度值Vpmeasured与预测纵波速度值Vpsat(W)完全相等),得到加权系数W;
Vpmeasured-Vppor(W)→0(13)
c)将加权系数W代入公式(2),得到配位数Cp,得到配位数Cp之后,把配位数Cp值代入计算随压力变化的干岩石切变模量的Mindlin公式中得到干岩石的切变模量μdry,进而将干岩石的切变模量μdry代入Gassmann方程即代入公式(5)计算得到预测横波速度即建立起了随压力变化的孔隙介质横波速度预测模型。
本发明利用已经建立好的随压力变化的孔隙介质横波速度预测模型,预测实际资料的加权系数W、配位数Cp以及预测随压力变化的孔隙介质横波速度,预测的加权系数W见表1。
2)由于在孔隙介质模型中,假设储层为孔隙介质,因此所用到的实际的采集的数据也假设为孔隙介质的测量结果,而当储层含有裂隙时这与实际情况是不符合的,因此本方法利用岩石物理测试资料对上述预测结果进行校正。
进行岩石物理测试的岩心相对于整个储层来说是各向同性的;如果该岩心所在的井与测井所在的井为同一口井,那么不需要标定层位,如果是两口井(即取岩心的井与进行数据采集的井不是同一口井),则需要根据地质资料标定岩心在采集测井数据所在井中的层位。对岩心进行岩石物理测试,得到岩心的纵波速度、横波速度、密度以及孔隙度参数,然后利用随压力变化的孔隙介质横波速度预测模型对这些数据进行计算得到岩心的加权系数Wcore,然后利用岩心的加权系数Wcore对加权系数W进行校正,得到校正后的各向同性介质的加权系数Wisotropy,同时,在利用岩石物理测试数据的时候,使岩心的差异压力与储层的差异压力要相同。标定好层位之后,利用公式(14)对加权系数W进行校正,得到校正后各向同性介质的加权系数Wisotropy:
其中Wisotropy为校正后的各向同性介质的加权系数,Wcore为计算得到的岩心的加权系数,Wlog为岩心对应层位的加权系数。校正之前的加权系数W利用随压力变化的孔隙介质模型预测出来。利用干岩石物理测试结果校正后得到的校正后各向同性介质的加权系数Wisotropy与利用饱和岩石物理测试结果校正后得到的Wisotropy见表1。表1为校正之前的加权系数W以及利用干岩石以及油水饱和岩石物理测试结果校正得到的加权系数W。校正之前的加权系数为利用孔隙介质模型预测的加权系数W值。该储层进行岩心岩石物理测试的时候,既对干岩石进行了测试,也对油水饱和岩石进行了测试。因此分别利用干岩石的预测结果和油水饱和岩石的预测结果进行校正。
表1预测的加权系数W
3)利用裂隙孔隙度将总孔隙度校正为孔隙孔隙度,同时对岩石的体积密度进行校正,得到孔隙介质的体积密度,具体过程为:将由测井数据得到的总孔隙度φz减去裂隙孔隙度φc得到孔隙孔隙度φp;然后对体积密度进行校正,方法如下:首先利用下述公式ρd=ρ*(1-φz)计算得到干岩石的体积密度,其中ρd为干岩石的体积密度,ρ为测井数据资料得到的体积密度,然后利用公式ρp=ρd/(1-φp)得到孔隙介质的体积密度,其中ρp为校正后得到的孔隙介质的体积密度,φp为孔隙孔隙度;然后将得到的校正后的各向同性介质的加权系数Wisotropy,以及校正过的孔隙度,体积密度代入孔隙介质横波速度预测模型,得到孔隙介质下随压力变化的纵波速度α以及横波速度β;
4)然后利用Thomson裂隙介质理论进行横波速度的预测;上述孔隙介质纵横波速度的预测以及Thomson裂隙介质预测横波速度的实现是通过C语言编程实现的,具体流程见图1。Thomsen于1986年提出了裂隙为垂向裂隙且有一个对称轴的弱各向异性的弹性介质的相速度公式如下:
其中α0、β0分别为对称方向的孔隙的纵波速度、横波速度,可以利用上述利用孔隙介质模型预测的纵波速度、横波速度以及公式(18)和公式(19)获得;ε、γ和δ为各向异性参数;θ为波前与对称轴的夹角。Vs1是平行于裂隙方向的横波速度,而Vs2则是垂直于裂隙方向的横波速度;这两个横波速度不同(即快横波速度与慢横波速度,在裂隙介质中横波会发生分裂,产生两个横波,按照速度划分,就可以成为快横波与慢横波,相应的速度就是快横波速度与慢横波速度),实际测井时,是先后接收到两个横波速度,因此称其为快横波和慢横波。
其中α,β为上述孔隙介质模型得到的随压力变化的孔隙介质的纵横波速度,α0、β0分别为储层对称方向的纵波速度、横波速度,ν为孔隙介质的泊松比,ε以及γ为各向异性参数。
其中各向异性参数分别为:
其中Dci为流体影响因子,在低频和高频声波测量时分别如下:
在上述公式中,Kma为岩石骨架的体变模量,Kf为裂隙中混合流体的体变模量,ν为孔隙介质的泊松比,E为杨氏模量,νd和Ed为相应的各向同性孔隙干岩石的泊松比和杨氏模量,νma则为岩石骨架的泊松比,φz为总孔隙度也就是孔隙孔隙度φp与裂隙孔隙度φc的和。
由于θ是未知的,公式(15)中Vp(θ)为预测纵波速度,由于预测纵波速度是趋近与实测纵波速度Vp,可以利用Vp与Vp(θ)相减为0,得到θ。因此首先利用实测纵波速度和公式(15)预测出θ值,然后将所得到的θ代入公式(16)和公式(17)得到两个横波速度。
上述孔隙介质纵横波速度的预测、角度θ的求取以及Thomson裂隙介质预测横波速度的实现是通过C语言编程实现的,具体流程见图1,具体过程为:输入差异压力p以及校正过的参数,包括孔隙介质的孔隙孔隙度φp、孔隙介质的体积密度ρp以及校正过的加权系数Wisotropy,利用利用考虑压力变化的Digby、Mindlin公式以及Gassmann方程的变形公式计算随压力变化的孔隙介质的纵波速度α、横波速度β;然后利用一个循环来求出波前与对称轴的夹角θ,首先设定一个初值即定义波前与对称轴的夹角θ=1(角度),将该θ值和孔隙介质的孔隙孔隙度φp、孔隙介质的密度ρp以及孔隙介质的纵波速度Vpmeasured代入预测纵波速度Vp(θ),得到预测纵波速度Vp(θ),然后将该波前与对称轴的夹角θ加上1得到新的θ值,然后判断新的θ值是否小于90度,如果小于90度,则继续将该新的θ值和孔隙介质的孔隙孔隙度φp、孔隙介质的密度ρp以及孔隙介质的纵波速度Vpmeasured代入预测纵波速度Vp(θ),然后在将θ值加上1判断新的θ值是否小于90度,若小于90度重复上述步骤,即将大于0度小于90度的每一个整数角度下的预测纵波速度Vp(θ)计算出来;若新的θ值不小于90度,则再继续下一步,即得到这89个预测纵波速度之后,将它们分别与实际纵波速度Vp相减,得到误差,找到误差最小的那个角度下的预测纵波速度Vp(θ),这个角度即为波前与对称轴的夹角θ;然后将得到的角度以及孔隙介质的孔隙孔隙度φp、孔隙介质的密度ρp、随压力变化的孔隙介质的纵波速度α、横波速度β代入Thomson公式,得到裂隙介质的横波速度并且输出,C语言程序结束。
图2是岩心为干岩石进行校正即利用表1中校正过的Wisotropy值,得到的快横波速度与实测快横波速度的对比,其中实线为实测值,虚线为预测值,其平均误差为7.742%。图3是岩心为干岩石进行校正即利用表1中校正过的加权系数W值,得到的慢横波速度与实测慢横波速度的对比,其中实线为实测值,虚线为预测值,其平均误差为7.633%。图4是岩心为饱和岩石进行校正即利用表1中校正过的W值,得到的快横波速度与实测快横波速度的对比,其中实线为实测值,虚线为预测值,其平均误差为7.478%。图5中岩心为干岩石进行校正即利用表1中校正过的W值,得到的慢横波速度与实测慢横波速度的对比,其中实线为实测值,虚线为预测值,其平均误差为7.274%。
5)根据步骤4)得到的平行于裂隙方向的横波速度Vs1和垂直于裂隙方向的横波速度Vs2,建立各向异性介质的AVO模型,预测油气藏岩石属性,包括孔隙度、密度、岩性与流体含量。
通过AVO分析,地球物理学家可以更好地评估油气藏岩石属性,包括孔隙度、密度、岩性与流体含量,而横波速度是建立AVO模型、转换波分析过程中的一个不可缺少的弹性参数,而在多数情况下研究工区是没有横波速度数据的,因此本发明给出了一个横波速度预测的方法。储层中还可能发育裂隙,那么如果利用孔隙介质横波速度预测模型与实际情况不符,而本发明以已经建立好的随压力变化孔隙介质横波速度预测模型为基础,利用岩石物理测试资料对预测的孔隙介质纵横波速度等参数进行校正,然后利用Thomson裂隙介质理论进行横波速度的预测,这样更加符合存在裂隙的储层。横波速度对压力是敏感的,本发明既考虑了压力对于介质横波速度的影响,同时也考虑到了储层中含有裂隙,更加符合实际。
本发明采集测井数据,建立起随压力变化的孔隙介质横波速度预测模型;利用由岩心的测试资料得到的加权系数对实际采集到的数据得到的加权系数进行校正,得到校正后的各向同性介质的加权系数;将校正过的各个参数代入随压力变化的孔隙介质横波速度预测模型得到孔隙介质的与压力相关的纵横波速度;根据孔隙介质下的与压力相关的纵波速度、横波速度得到裂隙介质对称方向的与压力相关的纵波速度、横波速度并带入Thomson提出的公式中,得到平行于裂隙方向的与压力相关的横波速度Vs1和垂直于裂隙方向的与压力相关的横波速度Vs2。
Claims (10)
1.一种随压力变化的裂隙介质横波速度预测方法,其特征在于,包括以下步骤:
1)采集测井数据:采集总孔隙度φz、岩石的体积密度ρ以及实际纵波速度Vp、流体饱和度;采集岩石骨架的体变模量Kma,岩石骨架的切变模量μma,差异压力p,测量配位数Cp',岩石颗粒变形之前接触区域的半径a与岩石颗粒的半径R;
分别对总孔隙度φz、岩石的体积密度ρ进行校正得到孔隙介质的孔隙孔隙度φp、孔隙介质的体积密度ρp,利用流体饱和度得到混合流体的体变模量Kf;计算岩石骨架的泊松比νma;
对测量配位数Cp'进行线性拟合并且加权,得到含有加权系数W的配位数Cp,然后利用Digby、Mindlin公式以及Gassmann方程的变形公式建立随压力变化的孔隙介质横波速度预测模型;
2)对岩心进行岩石物理测试,得到岩心的纵波速度、岩心的横波速度、岩心的密度以及岩心的孔隙度,然后利用随压力变化的孔隙介质横波速度预测模型对岩心的纵波速度、岩心的横波速度、岩心的密度以及岩心的孔隙度进行计算得到岩心的加权系数Wcore,然后利用岩心的加权系数Wcore对加权系数W进行校正,得到校正后的各向同性介质的加权系数Wisotropy,如公式(14)所示:
其中,Wlog为岩心对应层位的加权系数;
3)将校正后的各向同性介质的加权系数Wisotropy、孔隙介质的孔隙孔隙度φp以及孔隙介质下的密度ρp代入Gassmann方程的变形公式,得到孔隙介质下随压力变化的纵波速度α以及横波速度β;
4)根据孔隙介质下随压力变化的纵波速度α以及横波速度β,得到对称方向的弹性介质的纵波速度α0、横波速度β0;将对称方向的弹性介质的纵波速度α0、横波速度β0代入公式(15)-(17)中,得到平行于裂隙方向的横波速度Vs1和垂直于裂隙方向的横波速度Vs2;
其中θ为波前与对称轴的夹角,ε、γ和δ为各向异性参数,Vp(θ)为预测纵波速度;
5)根据步骤4)得到的平行于裂隙方向的横波速度Vs1和垂直于裂隙方向的横波速度Vs2,建立各向异性介质的AVO模型,预测油气藏岩石属性;
公式(3)-(5)为Gassmann方程的变形为:
μsat=μdry(4)
其中分别为预测的纵波速度、横波速度,μsat为孔隙介质的切变模量,Kdry为干岩石的体变模量、μdry为干岩石的切变模量,Kma为岩石骨架的体变模量;Kf为混合流体的体变模量,ρ为岩石的体积密度。
2.一种随压力变化的裂隙介质横波速度预测方法,其特征在于,所述步骤1)中若取岩心的井与进行采集测井数据的井是两口井,则标定岩心在进行采集数据的井中的层位。
3.根据权利要求1所述的一种随压力变化的裂隙介质横波速度预测方法,其特征在于,所述步骤1)中随压力变化的孔隙介质横波速度预测模型建立的具体过程为:
a)对测量配位数Cp'进行线性拟合,得到测量配位数Cp'与孔隙度φ的关系为:Cp'=11.759e1-φ-12.748(1)
对测量配位数Cp'进行加权,得到配位数Cp,如公式(2)所示,
Cp=W*Cp'(2)
其中,W为加权系数;
b)在均匀的孔隙介质的条件下,利用Gassmann方程的变形进行纵横波速度的预测;
将配位数Cp代入Digby公式中得到公式(7):
将配位数Cp代入Mindlin公式中得到公式(9):
其中,νma为岩石骨架的泊松比,μma为岩石骨架的切变模量;φ为孔隙度;Cp'为测量配位数;a为岩石颗粒变形之前接触区域的半径,b为岩石颗粒变形之后接触区域的半径,R为岩石颗粒的半径;
将干岩石的体变模量Kdry与干岩石的切变模量μdry代入Gassmann方程的变形公式中计算预测纵波速度利用预测的纵波速度值等于实际测量的纵波速度值,得到加权系数W;
c)将加权系数W代入公式(2),得到配位数Cp,再把配位数Cp值代入计算随压力变化的干岩石切变模量的Mindlin公式中得到干岩石的切变模量μdry,再将干岩石的切变模量μdry代入Gassmann方程的变形公式中计算得到预测横波速度即建立起了随压力变化的孔隙介质横波速度预测模型。
4.根据权利要求3所述的一种随压力变化的裂隙介质横波速度预测方法,其特征在于,所述步骤b)中Digby公式为:
5.根据权利要求3所述的一种随压力变化的裂隙介质横波速度预测方法,其特征在于,所述步骤b)和步骤c)中Mindlin公式为:
6.根据权利要求3所述的一种随压力变化的裂隙介质横波速度预测方法,其特征在于,所述步骤b)中
其中p为差异压力。
7.根据权利要求1所述的一种随压力变化的裂隙介质横波速度预测方法,其特征在于,所述步骤3)中将总孔隙度φz校正为孔隙介质下的孔隙孔隙度φp的具体过程为:将总孔隙度φz减去裂隙孔隙度φc得到孔隙孔隙度φp。
8.根据权利要求1所述的一种随压力变化的裂隙介质横波速度预测方法,其特征在于,所述步骤1)中对岩石的体积密度ρ进行校正,得到孔隙介质情况下的孔隙介质的密度ρp的具体过程如下:
首先利用公式ρd=ρ*(1-φz)计算得到干岩石的密度,其中ρd为干岩石的密度,然后利用公式ρp=ρd/(1-φp)得到孔隙介质的密度ρp。
9.根据权利要求1所述的一种随压力变化的裂隙介质横波速度预测方法,其特征在于,所述步骤4)中对称方向的弹性介质的纵波速度α0、横波速度β0是通过以下过程得到的:
根据公式(18)和(19),计算得到对称方向的弹性介质的纵波速度α0、横波速度β0;
其中α、β分别为孔隙介质下随压力变化的纵波速度、横波速度;ν为孔隙介质的泊松比,ε以及γ为各向异性参数。
10.根据权利要求1所述的一种随压力变化的裂隙介质横波速度预测方法,其特征在于,所述步骤4)中各向异性参数ε、γ和δ的表达式分别为:
公式(20)和(21)中的
公式(20)中Dci为流体影响因子,在低频和高频声波测量时分别如下:
公式(23)中
其中,Kma为岩石骨架的体变模量,Kf为混合流体的体变模量,ν为孔隙介质的泊松比,E为杨氏模量,νd和Ed为相应的各向同性孔隙干岩石的泊松比和杨氏模量,νma则为岩石骨架的泊松比,φz为总孔隙度,Dci(lo)为低频声波测量时的流体影响因子,Dci(mh)为高频声波测量时的流体影响因子。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310518774.9A CN103576195B (zh) | 2013-10-28 | 2013-10-28 | 一种随压力变化的裂隙介质横波速度预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310518774.9A CN103576195B (zh) | 2013-10-28 | 2013-10-28 | 一种随压力变化的裂隙介质横波速度预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103576195A CN103576195A (zh) | 2014-02-12 |
CN103576195B true CN103576195B (zh) | 2016-06-08 |
Family
ID=50048361
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310518774.9A Active CN103576195B (zh) | 2013-10-28 | 2013-10-28 | 一种随压力变化的裂隙介质横波速度预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103576195B (zh) |
Families Citing this family (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106368691B (zh) * | 2015-07-24 | 2021-05-28 | 中国石油化工股份有限公司 | 基于岩石物理地震信息三维异常孔隙压力预测方法 |
CN105607118B (zh) * | 2015-09-22 | 2018-01-09 | 西北大学 | Co2地质封存中四维多分量地震监测的转换波avo计算方法 |
CN106842326B (zh) * | 2015-12-04 | 2020-10-13 | 中国石油化工股份有限公司 | 无横波速度测井时砂泥互层co2地质封存时移地震正演模拟方法 |
CN106772608B (zh) * | 2017-02-16 | 2019-05-03 | 甘肃省地震局 | 等效孔隙裂缝介质的弹性阻抗及广义流体因子分析方法 |
CN107203006B (zh) * | 2017-06-14 | 2019-01-29 | 广东石油化工学院 | 基于纵横波速度的脆性矿物指数预测方法 |
CN108387931B (zh) * | 2018-01-22 | 2019-09-10 | 中国石油天然气股份有限公司 | 储层孔隙压力的确定方法和装置 |
CN109459498B (zh) * | 2018-09-14 | 2019-10-01 | 中国石油大学(华东) | 一种含倾斜裂缝横向各向同性岩石的声波速度计算方法 |
CN110954949A (zh) * | 2018-09-27 | 2020-04-03 | 中国石油化工股份有限公司 | 一种致密砂岩软孔隙度分布反演方法 |
CN109374497B (zh) * | 2018-10-29 | 2020-12-11 | 河海大学 | 一种岩石微观孔隙结构测试方法 |
CN109471168B (zh) * | 2018-11-06 | 2020-11-03 | 河海大学 | 一种孔裂隙介质中纵波速度与衰减的预测方法 |
CN110275206B (zh) * | 2019-08-12 | 2021-09-28 | 河海大学 | 一种裂隙-孔隙型岩石物理弹性模板 |
CN110596757B (zh) * | 2019-08-14 | 2021-01-05 | 西南石油大学 | 一种页岩地层纵波及横波速度的校正方法 |
CN114136556B (zh) * | 2021-11-30 | 2024-03-08 | 北京卫星环境工程研究所 | 一种基于波速修正的航天器复材结构泄漏定位方法 |
CN114674934B (zh) * | 2022-02-18 | 2023-06-13 | 河海大学 | 一种建立饱和重油岩石波速随温度变化的理论模型的方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4686657A (en) * | 1984-09-24 | 1987-08-11 | Mobil Oil Corporation | Matching P wave and shear wave sections in seismic exploration |
CN102353989A (zh) * | 2011-08-24 | 2012-02-15 | 成都理工大学 | 基于自适应基质矿物等效弹性模量反演的横波速度估算方法 |
-
2013
- 2013-10-28 CN CN201310518774.9A patent/CN103576195B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4686657A (en) * | 1984-09-24 | 1987-08-11 | Mobil Oil Corporation | Matching P wave and shear wave sections in seismic exploration |
CN102353989A (zh) * | 2011-08-24 | 2012-02-15 | 成都理工大学 | 基于自适应基质矿物等效弹性模量反演的横波速度估算方法 |
Non-Patent Citations (5)
Title |
---|
A Biot-Gassmann perspective;Russell H.B.等;《CREWES Research Report》;Fluid-property discrimination with AVO:;20011231;第13卷;第403-419页 * |
A Tutorial on Gassmann Fluid Substitution: Formulation, Algorithm and Matlab Code;Dhananjay Kumar;《Geohorizons》;20060430;第1-9页 * |
Gassmann"s equation and fluid-saturation effects on seismic velocities;Han De-hua等;《Geophysics》;20040430;第69卷(第2期);第398-405页 * |
有效压力对岩石纵横波速度的影响;马中高 等;《勘探地球物理进展》;20060630;第29卷(第3期);第0183-0186页 * |
横波速度预测方法;孙福利 等;《地球物理学进展》;20080430;第23卷(第2期);第470-474页 * |
Also Published As
Publication number | Publication date |
---|---|
CN103576195A (zh) | 2014-02-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103576195B (zh) | 一种随压力变化的裂隙介质横波速度预测方法 | |
Li et al. | Permeability evolution of shale under anisotropic true triaxial stress conditions | |
CN107203005B (zh) | 一种定量化计算裂缝描述参数的方法 | |
Da‐Xing | A study on the rock physics model of gas reservoir in tight sandstone | |
CN104155693A (zh) | 储层流体流度的角道集地震响应数值计算方法 | |
CN103713327B (zh) | 基于测井及地震资料的页岩气储层最小闭合压力评价方法 | |
CN104181585A (zh) | 一种地球物理勘探横波估算方法和系统 | |
CN103576196B (zh) | 一种随压力变化的孔隙介质横波速度预测方法 | |
CN110515126B (zh) | 一种含随机分布裂缝横向各向同性岩石的声速计算方法 | |
Jin et al. | Rock-physics-model-based pore type characterization and its implication for porosity and permeability qualification in a deeply-buried carbonate reservoir, Changxing formation, Lower Permian, Sichuan Bain, China | |
CN110275206B (zh) | 一种裂隙-孔隙型岩石物理弹性模板 | |
CN105588883A (zh) | 三维岩石力学参数获取方法和系统 | |
King | Recent developments in seismic rock physics | |
CN107728205A (zh) | 一种地层压力预测方法 | |
CN102230973A (zh) | 一种三维分步傅立叶粘滞声波深度偏移方法 | |
CN110456412B (zh) | 一种基于叠后地震数据识别碳酸盐岩储层流体饱和度的方法 | |
CN103364834B (zh) | 一种利用叠前地震频散分析预测储层渗透率的方法 | |
CN102914797B (zh) | 一种获得地层各向异性系数的方法及装置 | |
DONG et al. | Test on elastic anisotropic coefficients of gas coal | |
CN114994752A (zh) | 一种弹性波地震频率谐振勘探方法 | |
CN115061202B (zh) | 一种页岩含气性地震储层直接检测方法 | |
CN109709610B (zh) | 一种岩石裂缝探测方法及系统 | |
He et al. | Rock skeleton models and seismic porosity inversion | |
CN110471129A (zh) | 一种深层页岩高温高压下的各向异性岩石物理建模方法 | |
CN112505766B (zh) | 一种评价井外不同方位下裂缝发育程度的方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |