基于岩石物理模型的页岩气储层地应力测井预测方法
【专利摘要】本发明涉及基于岩石物理模型的页岩气储层地应力测井预测方法，建立考虑干酪根颗粒的页岩气储层岩石物理模型，预测测井纵、横波速度；在此基础上计算储层的最大、最小水平主应力，及破裂压力，以在没有横波测井的情况下也能够对页岩气储层进行精确的应力评价；本发明的有益效果是：通过建立符合页岩气储层特征的岩石物理模型，提高了预测速度的精度；在岩石物理模型基础上求取最大、最小水平主应力、破裂压力，使得在没有实测横波测井数据的情况下，依据常规测井曲线就能够对地下应力进行预测，且预测结果精度较高。
【专利说明】基于岩石物理模型的页岩气储层地应力测井预测方法
【技术领域】
[0001]本发明属于页岩气储层的地应力测井预测领域，尤其涉及基于岩石物理模型的页岩气储层地应力测井预测方法。
【背景技术】
[0002]对于页岩气储层，由于压裂改造直接影响着气产量，而对地下应力的预测决定着压裂改造的好坏。储层地应力的测井评价模型有多种，比较成熟的有Mohr Coulomb模型，黄荣樽模型(参见参考文献I )，单轴应变模型等。通过这些模型求取最大、最小水平主应力，进而求取破裂压力，应力梯度等，可以对储层进行评价；但是这些模型的应用必须以较为精确的杨氏模量、泊松比等岩石力学参数作为输入参数，而精确的岩石力学参数，必须要有精确的纵、横波速度才能计算得到，通常，纵、横波速度可以由偶极子测井获得。但是，大部分已钻井中没有横波测井资料，这给地应力的评价造成了难度。Mullen (参见参考文献9)等使用一个综合的经验公式求算横波，并使用单轴应力公式计算最小闭合应力以对地下的应力场进行评价。经验公式求算横波受地区具体地质情况影响，很难求取准确。通过建立符合具体岩石物理特征的岩石物理模型，可以计算得到比较准确的横波。页岩气储层岩性复杂，微观孔隙类型多样，富含干酪根，针对这种非常规储层，常规的岩石物理模型不能够求算较为准确的速度，经典的Xu-White模型(参见参考文献2)，将砂岩和粘土混合，使用二维的Kuster- ToksBz模型(参见参考文献6)求取带孔隙的干岩石的弹性模量，使用Gassmann(参见参考文献10)方程求取饱和流体的岩石的弹性模量，但是这个模型只考虑了粘土和砂岩，岩性不复杂，且孔隙类型也较为简单；在此基础上发展的Xu-Payne岩石物理模型(参见参考文献7)及DEM-Gassmann岩石物理模型(参见参考文献3)都使用三维孔隙谱计算含不同孔隙类型岩石的体积模量和剪切模量，得到较好的效果，但是这两个模型都只适合于比较简单的岩性，且都没有考虑干酪根，不适合于页岩气储层。Bandyopadhyay等(参见参考文献8)使用Brown-Korriga (参见参考文献4)模型,采用固体替换的方式,计算含干酪根岩石的弹性模量，但是根据岩心观察，干酪根是以一定颗粒形态分散存在于泥页岩中的，这种方式模拟的不是有机质在岩石中的真实状态。
【发明内容】
[0003]本发明的目的是提出基于岩石物理模型的页岩气储层地应力测井预测方法的技术方案，建立考虑干酪根颗粒的页岩气储层岩石物理模型，真实的反映在岩石中分散存在的干酪根的的弹性特征；预测测井纵、横波速度；在此基础上计算储层的最大、最小水平主应力，及破裂压力，以在没有横波测井的情况下也能够对页岩气储层进行精确的应力评价；
[0004]为了实现上述目的，本发明的技术方案是:基于岩石物理模型的页岩气储层地应力测井预测方法，建立包括干酪根颗粒的页岩气储层岩石物理模型，在常规测井曲线的基础上预测测井纵、横波速度；在此基础上计算储层的最大、最小水平主应力，及破裂压力；所建立的岩石物理模型，包括矿物，干酪根，孔隙及其中所含的流体，所述孔隙中的流体包括油、气、水，其中将干酪根等效为硬币状裂缝形态的颗粒，使用颗粒纵横比表示其大小和形状，以微分方式计算含干酪根的岩石的弹性参数，得到在岩石中分散存在的干酪根的弹性特征；
[0005]所述方法包括以下步骤:
[0006]步骤1:通过测井复杂岩性分析程序进行测井解释，获得矿物体积、岩石孔隙度、流体饱和度、干酪根体积，进而求得矿物在岩石中的体积分数，所述矿物包括粘土、石英、方解石；使用交会分析获得每个矿物的弹性模量和密度，使用已知文献中的经验值作为孔隙中流体的弹性模量和密度，所述弹性模量指体积模量、剪切模量；
[0007]步骤2:根据步骤I中所得的矿物的体积模量和剪切模量，矿物体积分数，使用等应力和等应变平均的算术平均估算岩石基质的体积模量和剪切模量，所述岩石基质指不含孔隙和流体的固体岩石部分；在此基础上，将岩石基质与干酪根颗粒混合，即计算含干酪根颗粒的岩石的体积模量和剪切模量，将干酪根等效为硬币状裂缝形态的颗粒，使用颗粒纵横比表示干酪根颗粒的形状，使用微分等效介质方法计算含所述干酪根颗粒的岩石的体积模量和剪切模量；在所述的微分等效介质方法中，形状因子使用Berryman的三维包含物理论中的硬币状裂缝公式进行计算，所述颗粒纵横比指硬币状裂缝形状颗粒的短轴和长轴之t匕，在岩石物理模型中颗粒纵横比的取值范围为O?1，作为迭代计算的变量参数；
[0008]步骤3:将步骤I所得的孔隙度划分为砂岩矿物孔隙度、灰岩矿物孔隙度、粘土矿物孔隙度三种；将灰岩矿物孔隙度、砂岩矿物孔隙度等效为刚性孔隙占多数，柔性孔隙占少数的孔隙谱；将粘土矿物孔隙度等效为柔性孔隙占多数，刚性孔隙占少数的孔隙谱；将不同的孔隙度与步骤2所得到的含干酪根的岩石基质相混合，这个过程使用微分等效介质方法计算，得到的结果为包括孔隙、干酪根、岩石基质的干岩石的体积模量、剪切模量；所述干岩石为不含流体的岩石固体部分。
[0009]步骤4:将流体与步骤3所得到的干岩石相混合，该过程使用流体替换方程计算，所得到的结果为饱和流体岩石的体积模量和剪切模量；所述的流体指混合流体，是在步骤I得到的饱和度及流体的弹性模量的基础上，使用等应力平均公式计算不同饱和度的油、气、水混合时的体积模量，混合流体的剪切模量为O ;根据步骤I所得的矿物体积、孔隙度、饱和度、干酪根体积，使用算术平均计算饱含流体岩石的密度；所述流体替换方程为Gassmann方程,所述算术平均计算饱含流体岩石的密度的方法采用Wood公式；
[0010]步骤5:由步骤4所得的饱和流体岩石的体积模量、剪切模量、密度计算饱和流体岩石的纵、横波速度，即预测纵波速度和预测横波速度，对比实测纵波速度与预测纵波速度，相减求误差，当误差大于限定值时，返回步骤2，调整所述颗粒纵横比的取值范围，直到误差小于限定值或误差达到极小值，当误差小于限定值或达到极小值时，输出预测的纵、横波速度，这是一个迭代的过程；当测井有实测横波速度时，也需要对预测横波速度与实测横波速度进行对比求误差并用于判别迭代是否结束。
[0011]步骤6:使用步骤5所得到的预测纵、横波速度计算动态杨氏模量和泊松比，使用地区经验公式将动态杨氏模量和泊松比转换为静态杨氏模量和泊松比；所述的地区经验公式是张晋言和孙建孟关系式；使用测井数据求取上覆地层压力和孔隙压力，所述测井数据特指密度。[0012]步骤7:将步骤6所得的静态杨氏模量和泊松比，上覆地层压力、孔隙压力代入黄荣樽应力模型计算最大水平主应力、最小水平主应力，在这个过程中，需要实际工区岩心的应力测试数据对黄荣樽应力模型中的参数进行标定；根据应力参数之间的关系，计算破裂压力；综合最大水平主应力、最小水平主应力、破裂压力，进行地层应力评价。
[0013]更进一步，步骤I中所述复杂岩性分析程序，是Atlas公司的复杂岩性程序。
[0014]更进一步，步骤I中，所述的交会分析是指，在测量有纵波、横波、密度测井数据的井中，根据纵横波、密度计算体积模量、剪切模量，根据测井、录井、钻井及地质信息综合判断出井筒中不同矿物较纯的岩性段，分别将这些层段的体积模量、剪切模量、密度数据使用交会图的方法进行显示，选择数据点较为集中的位置作为纯的矿物的体积模量、剪切模量和密度数据。
[0015]本发明的有益效果是:所建立岩石物理模型可更真实的反映页岩气储层的复杂矿物以及复杂矿物造成的多孔隙类型，即以一定形态分散存在的干酪根的岩石物理特征，从而在使用该岩石物理模型进行储层速度预测的过程中提高预测速度的精度；在此基础上的应力预测使得在没有测井横波速度的情况下，根据常规测井曲线可以计算得到与测量值吻合的地下的最大、最小水平主应力，以及破裂压力，是一种可靠的页岩气储层地应力测井预测方法。
[0016]下面结合附图和实施例对本发明作一详细描述。
【专利附图】
【附图说明】
[0017]图1是基于岩石物理模型的页岩气储层地应力测井预测技术流程图；
[0018]图2是所建立包含干酪根的岩石物理模型的流程示意图；
[0019]图3是建立的新岩石物理模型与经典的Xu-White模型速度预测结果对比图，从左往右，图道依次为地层，深度，录井岩性，气测甲烷与烃类，新岩石物理模型预测纵波时差与实测纵波时差，Xu-White岩石物理模型预测纵波时差与实测纵波时差，新岩石物理模型预测横波时差与实测横波时差，Xu-White岩石物理模型预测横波时差与实测横波时差；岩心图片为箭头所指位置的储层特征；图道中，CH4指甲烧，Hydrocarbon指烃类，Acpre_new指本发明预测的纵波速度，Dspre_new指本发明预测的横波速度，Acpre_Xff指Xu-White模型预测的纵波速度，Dspre_Xff指Xu-White模型预测的横波速度，AC指测量纵波时差，DTS指测量横波时差；所述时差为速度的导数；
[0020]图4为建立的新岩石物理模型与经典的Xu-White岩石物理模型预测速度与实测速度之间的相对误差的对比图，a)，b)图为新模型预测速度的纵、横波误差；c)，d)为Xu-White模型预测速度的纵、横波误差；
[0021]图5是页岩气储层应力计算结果，从左往右，图道依次为地层，伽马与井径曲线，深度，录井岩性，气测甲烷与烃类，静态泊松比，静态杨氏模量，实测与预测上覆地层压力、孔隙压力，实测与预测最小水平主应力，实测与预测最大水平主应力，破裂压力；其中，实测的最大、最小水平主应力，上覆地层压力为杆状离散点，预测最大、最小水平主应力，上覆地层压力为曲线。
【具体实施方式】[0022]基于岩石物理模型的页岩气储层地应力测井预测方法，建立考虑干酪根颗粒的页岩气储层岩石物理模型，在常规测井曲线的基础上预测测井纵、横波速度；在此基础上计算储层的最大、最小水平主应力，及破裂压力，以在没有横波测井的情况下根据常规测井曲线对页岩气储层进行精确的应力评价；所建立的岩石物理模型，包括矿物，干酪根，干酪根，孔隙及其中所含的流体，所述孔隙中的流体包括油、气、水，其中将干酪根等效为硬币状裂缝形态的颗粒，使用颗粒纵横比表示其大小和形状，以微分方式计算含干酪根的岩石的弹性参数，得到在岩石中分散存在的干酪根的弹性特征。
[0023]以下实施例结合中国南方某页岩气盆地侏罗系下统页岩气储层为例来阐述本发明的【具体实施方式】:
[0024]如图1所示，为基于岩石物理模型的页岩气储层地应力测井预测方法的流程图，首先通过测井解释得到地层中的矿物体积、孔隙度、饱和度、干酪根体积；通过交会分析得到主要矿物的弹性模量和密度，通过文献调研得到流体的弹性模量和密度，此处弹性参数特指体积模量、剪切模量；将测井解释结果引入新建立的页岩气储层岩石物理模型，通过调整岩石物理模型中的孔隙类型分布谱参数，迭代计算纵波速度与横波速度。使用预测纵、横波速度计算动态杨氏模量和泊松比，并根据地区实测数据总结的经验关系式将动态的杨氏模量和泊松比转换为静态杨氏模量和泊松比；由测井数据计算上覆地层压力和孔隙压力，此处测井数据特指密度和声波时差；将静态的杨氏模量、泊松比，上覆地层压力、孔隙压力代入黄荣樽应力模型中计算最大水平主应力和最小水平主应力，在这个过程中需要根据实际工区岩心应力测试数据标定黄荣樽应力模型中的参数；在此基础上计算破裂压力；在最大水平主应力，最小水平主应力、破裂压力的基础上就可以对地层应力进行评价。
[0025]对图1所述流程共分7个步骤进行详细说明，包括两部分主要内容，第一部分为岩石物理模型的构建和基于岩石物理模型的速度预测，包括下述步骤1、2、3、4、5，图2为岩石物理模型的流程图；第二部分为应力的预测，包括下述步骤6、7。以下结合公式即实例进行详细说明。
[0026]步骤1:通过测井复杂岩性分析程序进行测井解释，获得矿物体积、岩石孔隙度、流体饱和度、干酪根体积，进而求得矿物在岩石中的体积分数，所述矿物包括粘土、石英、方解石；使用交会分析获得每个矿物的弹性模量和密度，使用已知文献中的经验值作为孔隙中流体的弹性模量和密度，所述弹性模量指体积模量、剪切模量。
[0027]所述复杂岩性分析程序，是Atlas公司的复杂岩性程序；通过该程序，可以根据常规测井曲线计算较为精细的泥质含量，矿物体积分数，孔隙度，流体饱和度等参数。其具体公式记载在雍世和、张超谟主编的于2002年由中国石油大学出版社出版的《测井数据处理与综合解释》书的313-319页，参见参考文献12。
[0028]所述的交会分析是指，在测量有纵波、横波、密度测井数据的井中，根据纵横波、密度计算体积模量、剪切模量，根据测井、录井、钻井及地质信息综合判断出井筒中不同矿物较纯的岩性段，分别将这些层段的体积模量、剪切模量、密度数据使用交会图的方法进行显示，选择数据点较为集中的位置作为纯的矿物的体积模量、剪切模量和密度数据。
[0029]本实施例中，得到的矿物、流体的弹性模量和密度如表I所示。
[0030]表I本发明中矿物及流体的性质参数
[0031]
【权利要求】
1.基于岩石物理模型的页岩气储层地应力测井预测方法，其特征在于，建立包括干酪根颗粒的页岩气储层岩石物理模型，在常规测井曲线的基础上预测测井纵、横波速度；在此基础上计算储层的最大、最小水平主应力，及破裂压力；所建立的岩石物理模型，包括矿物，干酪根，孔隙及其中所含的流体，所述孔隙中的流体包括油、气、水，其中将干酪根等效为硬币状裂缝形态的颗粒，使用颗粒纵横比表示其大小和形状，以微分方式计算含干酪根的岩石的弹性参数，得到在岩石中分散存在的干酪根的弹性特征； 所述方法包括以下步骤: 步骤1:通过测井复杂岩性分析程序进行测井解释，获得矿物体积、岩石孔隙度、流体饱和度、干酪根体积，进而求得矿物在岩石中的体积分数，所述矿物包括粘土、石英、方解石；使用交会分析获得每个矿物的弹性模量和密度，使用已知文献中的经验值作为孔隙中流体的弹性模量和密度，所述弹性模量指体积模量、剪切模量； 步骤2:根据步骤I中所得的矿物的体积模量和剪切模量，矿物体积分数，使用等应力和等应变平均的算术平均估算岩石基质的体积模量和剪切模量，所述岩石基质指不含孔隙和流体的固体岩石部分；在此基础上，将岩石基质与干酪根颗粒混合，即计算含干酪根颗粒的岩石的体积模量和剪切模量，将干酪根等效为硬币状裂缝形态的颗粒，使用颗粒纵横比表示干酪根颗粒的形状，使用微分等效介质方法计算含所述干酪根颗粒的岩石的体积模量和剪切模量；在所述的微分等效介质方法中，形状因子使用Berryman的三维包含物理论中的硬币状裂缝公式进行计算，所述颗粒纵横比指硬币状裂缝形状颗粒的短轴和长轴之比，在岩石物理模型中颗粒纵横比的取值范围为O~1，作为迭代计算的变量参数； 步骤3:将步骤I所得的孔隙度划分为砂岩矿物孔隙度、灰岩矿物孔隙度、粘土矿物孔隙度三种；将灰岩矿物孔隙度、砂岩矿物孔隙度等效为刚性孔隙占多数，柔性孔隙占少数的孔隙谱；将粘土矿物孔隙度等效为柔性孔隙占多数，刚性孔隙占少数的孔隙谱；将不同的孔隙度与步骤2所得到的含干酪根的岩石基质相混合，这个过程使用微分等效介质方法计算，得到的结果为包括孔隙、干酪根、岩石基质的干岩石的体积模量、剪切模量；所述干岩石为不含流体的岩石固体部分； 步骤4:将流体与步骤3所得到的干岩石相混合，该过程使用流体替换方程计算，所得到的结果为饱和流体岩石的体积模量和剪切模量；所述的流体指混合流体，是在步骤I得到的饱和度及流体的弹性模量的基础上，使用等应力平均公式计算不同饱和度的油、气、水混合时的体积模量，混合流体的剪切模量为O ;根据步骤I所得的矿物体积、孔隙度、饱和度、干酪根体积，使用算术平均计算饱含流体岩石的密度；所述流体替换方程为Gassmann方程，所述算术平均计算饱含流体岩石的密度的方法采用Wood公式； 步骤5:由步骤4所得的饱和流体岩石的体积模量、剪切模量、密度计算饱和流体岩石的纵、横波速度，即预测纵波速度和预测横波速度，对比实测纵波速度与预测纵波速度，相减求误差，当误差大于限定值时，返回步骤2，调整所述颗粒纵横比的取值范围，直到误差小于限定值或误差达到极小值，当误差小于限定值或达到极小值时，输出预测的纵、横波速度，这是一个迭代的过程；当测井有实测横波速度时，也需要对预测横波速度与实测横波速度进行对比求误差并用于判别迭代是否结束； 步骤6:使用步骤5所得到的预测纵、横波速度计算动态杨氏模量和泊松比，使用地区经验公式将动态杨氏模量和泊松比转换为静态杨氏模量和泊松比；所述的地区经验公式是张晋言和孙建孟关系式；使用测井数据求取上覆地层压力和孔隙压力，所述测井数据特指密度； 步骤7:将步骤6所得的静态杨氏模量和泊松比，上覆地层压力、孔隙压力代入黄荣樽应力模型计算最大水平主应力、最小水平主应力，在这个过程中，需要实际工区岩心的应力测试数据对黄荣樽应力模型中的参数进行标定；根据应力参数之间的关系，计算破裂压力；综合最大水平主应力、最小水平主应力、破裂压力，进行地层应力评价。
2.根据权利要求1所述的基于岩石物理模型的页岩气储层地应力测井预测方法，其特征在于，步骤I中所述复杂岩性分析程序，是Atlas公司的复杂岩性程序。
3.根据权利要求1所述的基于岩石物理模型的页岩气储层地应力测井预测方法，其特征在于，步骤I中，所述的交会分析是指，在测量有纵波、横波、密度测井数据的井中，根据纵横波、密度计算体积模量、剪切模量，根据测井、录井、钻井及地质信息综合判断出井筒中不同矿物较纯的岩性段，分别将这些层段的体积模量、剪切模量、密度数据使用交会图的方法进行显示，选择数据点较为集中的位置作为纯的矿物的体积模量、剪切模量和密度数据。
【文档编号】G01V1/40GK103792581SQ201310751429
【公开日】2014年5月14日 申请日期:2013年12月31日 优先权日:2013年12月31日
【发明者】刘致水, 孙赞东 申请人:刘致水, 孙赞东, 王招明, 韩剑发