一种页岩各向异性岩石物理建模方法技术

本发明专利技术提供了一种页岩各向异性岩石物理建模方法，属于岩石物理研究领域。本方法包括：(1)计算混合矿物的弹性模量；(2)将各向异性自洽模型和各向异性微分有效介质模型结合起来计算粘土和流体复合物的有效弹性张量；(3)将粘土和流体复合物、层间缝与矿物组分进行混合，计算饱和岩石的有效弹性参数，求取等效各向异性介质的纵波速度和横波速度，即完成建模。本发明专利技术可以被应用于预测页岩的纵横波速度，计算精度优于各向同性有效介质模型，更加准确的描述了岩石的弹性特征。

【技术实现步骤摘要】

本专利技术属于岩石物理研究领域，具体涉及一种页岩各向异性岩石物理建模方法。
技术介绍
地球物理方法在页岩气勘探开发过程中起到了重要的技术支持作用，它在研究储层的复杂特征方面具有独特的优势。岩石物理技术是地球物理研究的基础。岩石物理建模和岩石物理分析是岩石物理研究的两个重要方面，岩石物理建模是为了模拟介质的弹性性质所提出的一种等效模型，而岩石物理分析则是基于合适的岩石物理模型，分析岩石的微观结构对宏观响应特征的影响。在地球物理领域，地震岩石物理的发展越来越受到重视。研究岩石物理性质通常将岩石简化为双相介质。岩石骨架由一种矿物组成，为一相；孔隙流体由液体或气体组成，为另一相。基于这种认识，在岩石物理模型建立方面，许多学者做了大量的工作，1928年至今已建立了近20种理论方法。最早的理论模型是70多年前Voigt的等应变模型与Reuss的等应力模型。从岩石的结构方面出发，各种理论模型大致可分为三类：对矿物的性质进行体积平均的有效介质理论，如Wood方程和Wyllie时间平均方程；研究岩石内部球形孔隙、椭球形孔隙及裂缝对岩石性质影响的自适应理论，如Gassman方程、Biot理论及模型；研究岩石颗粒的有效弹性的接触理论，如Hertz理论、Mindlin理论及Digby模型。有效介质理论首先应用在材料科学领域，并用以研究复合材料的弹性性质。其中，有些理论给出了准确的解析解，另一些理论给出了复合材料弹性参数的上下边界。目前应用较为广泛的有效介质理论仅适用于各向同性岩石，没有考虑各向异性。页岩的微观结构相对于常规储层更加复杂，它的孔隙和颗粒的尺度更小，矿物组分更加多样，而且由于岩石骨架结构中粘土的定向排列和层间缝的存在，页岩储层具有强各向异性特征，可见常规的岩石物理模型很难适用。
技术实现思路
本专利技术的目的在于解决上述现有技术中存在的难题，提供一种页岩各向异性岩石物理建模方法，构建合适的页岩岩石物理模型，为页岩气储层参数与地震响应之间建立联系，指导地震反演和解释。本专利技术是通过以下技术方案实现的：一种页岩各向异性岩石物理建模方法，包括：(1)计算混合矿物的弹性模量；(2)将各向异性自洽模型和各向异性微分有效介质模型结合起来计算粘土和流体复合物的有效弹性张量；(3)将粘土和流体复合物、层间缝与矿物组分进行混合，计算饱和岩石的有效弹性参数，求取等效各向异性介质的纵波速度和横波速度，即完成建模。所述步骤(1)是这样实现的：利用Voigt-Reuss-Hill平均计算混合矿物的弹性模量：Km=KV+KR2]]>Um=UV+UR2]]>ρm＝Vcal*ρcal+Vqua*ρqua+Vker*ρkerKV＝Vcal*Kcal+Vqua*Kqua+Vker*Kker1KR=VcalKcal+VquaKqua+VkerKker]]>UV＝Vcal*Ucal+Vqua*Uqua+Vker*Uker1UR=UcalUcal+UquaUqua+UkerUker]]>(1)其中，Kcal和Ucal分别为方解石的体积模量和剪切模量，Kqua和Uqua分别为石英的体积模量和剪切模量，Kker和Uker分别为干酪根的体积模量和剪切模量，Vcal、Vqua和Vker分别为基质体积百分比归一化后的方解石、石英和干酪根的含量，Km和Um分别为混合矿物的体积模量和剪切模量，ρ表示密度。所述步骤(2)是这样实现的：应用各向异性自洽模型，岩石的刚度张量表示为：其中，I是单位张量，是与孔隙几何形状有关的张量。Cn是第n种成分刚度张量，vn是第n相的体积含量，Cp是第p种成分刚度张量，vp是第p相的体积含量，CSCA是由各向异性SCA模型估算的弹性张量；采用各向异性微分有效介质模型计算岩石有效弹性张量的表达式如下：式中各个符号所代表的物理量与(2)式相同；将公式(2)的解作为求解微分方程(3)的初始值，不断迭代孔隙含量，来求解最终的有效弹性张量，即先用公式(2)计算孔隙度为50％时岩石的有效弹性模量，将该值作为初始值，然后运用公式(3)不断迭代孔隙含量，一直调整到岩石的真实孔隙度，求得最终的有效弹性张量。所述步骤(3)是利用Backus平均实现的，具体如下：在长波极限下，一个由多层横向各向同性材料组成的层状介质，是等效各向异性的，其等效刚度是：ABF000BAF000FFC000000D000000D000000M,M=12(A-B)---(4)]]>其中A＝〈a-f2c-1〉+〈c-1〉-1〈fc-1〉2B＝〈b-f2c-1〉+〈c-1〉-1〈fc-1〉2C＝〈c-1〉-1F＝〈c-1〉-1〈fc-1〉D＝〈d-1〉-1M＝〈m〉式中a，b，c，d和f是五个独立的弹性常数，当每一个单层是各向同性时，a＝c＝λ+2μ，b＝f＝λ，d＝m＝μ，〈·〉表示对括号内属性按体积比的加权平均；利用步骤(2)中获得的CDEM，求得每一个单层的λ和μ，然后再带入公式(4)，获得最终的等效参数，求取等效各向异性介质的纵波速度和横波速度，即完成建模。所述步骤(3)中的所述利用步骤(2)中获得的CDEM，求得每一个单层的λ和μ，然后再带入公式(4)，获得最终的等效参数是这样实现的：张量CDEM中包含了分量C33和C55，利用这个两个参数，获得粘土和流体混合后的纵波速度Vp=C33/ρ]]>和横波速度Vs=C55/ρVs=D/ρ,]]>其中Vclay为粘土体积，Vf为流体体积；利用下式求取每一个单层的λ和μ：λ=ρVP2-2ρVS2]]>和μ=ρVS2;]]>然后将该结果和将步骤(1)获得的Km和Um代入和μ＝Um得到的结果，一起带入公式(4)，获得最终的等效参数。所述步骤(3)中的所述求取等效各向异性介质的纵波速度和横波速度是这样实现的：纵波速度横波速度其中ρ＝Vmρm+Vcalyρclay+Vfρf，Vm为基质体积，Vclay为粘土体积，Vf为流体体积。与现有技术相比，本专利技术的有益效果是：本方法构建的页岩各向异性岩石物理模型综合了各向异性自洽模型(SCA)、各向异性微分等效介质模型(DEM)和Backus平均，考虑了粘土矿物的定向排列和层间缝，通过输入岩石基质的含量、体积模量、剪切模量和孔隙纵横比，可以被应用于预测页岩的纵横波速度，计算精度优于各向同性有效介质模型，更加准确的描述了岩石的弹性特征。附图说明图1是本专利技术方法的步骤框图。本文档来自技高网...

【技术保护点】
一种页岩各向异性岩石物理建模方法，其特征在于：所述方法包括：(1)计算混合矿物的弹性模量；(2)将各向异性自洽模型和各向异性微分有效介质模型结合起来计算粘土和流体复合物的有效弹性张量；(3)将粘土和流体复合物、层间缝与矿物组分进行混合，计算饱和岩石的有效弹性参数，求取等效各向异性介质的纵波速度和横波速度，即完成建模。

【技术特征摘要】
1.一种页岩各向异性岩石物理建模方法，其特征在于：所述方法包括：
(1)计算混合矿物的弹性模量；
(2)将各向异性自洽模型和各向异性微分有效介质模型结合起来计算粘土
和流体复合物的有效弹性张量；
(3)将粘土和流体复合物、层间缝与矿物组分进行混合，计算饱和岩石的
有效弹性参数，求取等效各向异性介质的纵波速度和横波速度，即完成建模。
2.根据权利要求1所述的页岩各向异性岩石物理建模方法，其特征在于：所述
步骤(1)是这样实现的：
利用Voigt-Reuss-Hill平均计算混合矿物的弹性模量：
Km=KV+KR2]]>Um=UV+UR2]]>ρm＝Vcal*ρcal+Vqua*ρqua+Vker*ρkerKV＝Vcal*Kcal+Vqua*Kqua+Vker*Kker1KR=VcalKcal+VquaKqua+VkerKker]]>UV＝Vcal*Ucal+Vqua*Uqua+Vker*Uker1UR=UcalUcal+UquaUqua+UkerUker---(1)]]>其中，Kcal和Ucal分别为方解石的体积模量和剪切模量，Kqua和Uqua分别为
石英的体积模量和剪切模量，Kker和Uker分别为干酪根的体积模量和剪切模量，
Vcal、Vqua和Vker分别为基质体积百分比归一化后的方解石、石英和干酪根的含量，

\tKm和Um分别为混合矿物的体积模量和剪切模量，ρ表示密度。
3.根据权利要求2所述的页岩各向异性岩石物理建模方法，其特征在于：所述
步骤(2)是这样实现的：
应用各向异性自洽模型，岩石的刚度张量表示为：
其中，I是单位张量，是与孔隙几何形状有关的张量。Cn是第n种成分
刚度张量，vn是第n相的体积含量，Cp是第p种成分刚度张量，vp是第p
相的体积含量，CSCA是由各向异性SCA模型估算的弹性张量；
采用各向异性微分有效介质模型计算岩石有效弹性张量的表达式如下：
式中各个符号所代表的物理量与(2)式相同；
将公式(2)的解作为求解微分方程(3)的初始值，不断迭代孔隙含量，
来求解最终的有效弹性张量，即先用公式(2)计算孔隙度为50％时岩石的有效
弹性模量，将该值作为初始值，然后运用公式(3)不断迭代孔隙含量，...

【专利技术属性】
技术研发人员：李呈呈，马中高，刘卫华，周枫，张克非，周中彪，
申请(专利权)人：中国石油化工股份有限公司，中国石油化工股份有限公司石油物探技术研究院，
类型：发明
国别省市：北京;11