一种用于页岩气井储层脆性综合评价的方法技术

本发明专利技术公开了一种用于页岩气井储层脆性综合评价的方法，包括如下步骤：利用待测的页岩岩样进行三轴压缩测试实验，获取全应力‑应变曲线；计算脆性跌落系数R、应力降落系数P、软化模量M，并分别将三者对应的脆性指标进行归一化处理，得到综合脆性指数Bd；再根据Bd计算模型进行修正，以修正后的脆性指数Bd3为目标函数，以测井数据中页岩动态弹性模量E等参数为自变量，进行函数回归，得出Bd3的多元回归函数Y，并代入不同层段储层测井参数，得到不同井段储层的综合脆性指数。本方法考虑了页岩破坏的整体情况，可以综合反映页岩的脆性特征，同时可以结合页岩气井测井资料对不同深度、不同围压、不同井段储层的脆性进行分析，适应性较强。

【技术实现步骤摘要】

本专利技术涉及一种非常规油气储层评价
，尤其涉及一种用于页岩气井储层脆性综合评价的方法。
技术介绍
页岩气藏在国内分布广泛，可采储量为36.08×1012m3，居世界第一，但是页岩本身具有低孔低渗特征，一般都需经过大规模压裂改造才能获得商业产量，研究发现页岩的脆性能够显著影响井壁的稳定性以及压裂效果，是评价储层力学特性的关键指标，是遴选射孔改造层段和设计压裂规模的重要基础，因此对页岩脆性的研究具有重要意义。脆性是反映页岩可压性的关键参数，亦是材料的综合特性，实质是在非均匀力下产生的一种由局部破坏演变为多维破裂面的动态破坏过程。目前页岩脆性评价指标是采用矿物脆性指数还是力学脆性指数尚无统一认识，且研究成果大多出于学者对于各自研究目的提出的，缺乏统一的标准和方法，反映了页岩在矿物组成或岩石力学上的脆性特征比较单一，难以全面反映页岩的可压性特征。以岩石矿物组分、杨氏模量及泊松比等传统静态参数定义的脆性特征相对单一，一方面无法描述岩石破裂过程中的脆性变化和能量释放，另一方面未考虑不同围压条件对岩石脆塑性影响，主要从宏观上定性比较脆性的相对强弱，对相同储层地质条件的页岩较为适用，用于分析不同埋深、构造应力环境下的岩石脆性则存在较大的局限性。以力学脆性指数为代表的专利申请号CN104865124A，题名《基于岩石应力-应变曲线和超声波纵波速度的页岩脆性指数测定方法》，其特征在于基于三轴压缩应力-应变全曲线和超声波纵波速度进行页岩脆性指数测定，具体采用动静态结合的三轴压缩实验，获取实验过程中各个时间点的应力-应变全曲线和超声波纵波速度；利用纵波速度的变化曲线确定页岩微裂缝开始发生损伤破坏的时间点；根据纵波速度变化和应力-应变全曲线的形状，把应力-应变全曲线分成微裂隙闭合、压实后至微裂缝开始破坏、微裂隙扩展至破坏失稳和岩石破坏失稳后等4个阶段；利用应力-应变全曲线，计算相应阶段的页岩试件吸收的单位体积能量；利用弹性阶段吸收的单位体积能量与吸收的总单位体积能量的比值计算页岩的脆性指数。上述专利的目的是提供页岩各个阶段力学性能的综合计算方法，提高岩石脆性评价的准确性和合理性。此方法主要基于室内试验方法对某一深度点的岩样进行脆性指数测试评价，对于不同深度、不同井段页岩气储层脆性进行综合评价。其主要缺点如下：(1)主要基于室内试验测试数据进行分析，无法应用于现场工程参数设计。(2)用于分析不同深度、地质力学环境下的页岩岩石脆性则存在较大的局限性。
技术实现思路
本专利技术所要解决的技术问题在于针对上述现有技术存在的不足提供一种用于页岩气井储层脆性综合评价的方法，该方法可以方便的计算页岩气井不同位置储层的脆性指数，可用于指导页岩气井压裂施工射孔层段选择。本专利技术所采用的技术方案为：一种用于页岩气井储层脆性综合评价的方法，其特征在于：包括如下步骤：(1)、利用待测的页岩岩样进行三轴压缩测试实验，获取全应力-应变曲线；(2)、根据实验获得的应力-应变全曲线，在曲线中读取峰值应力σa、峰值应变εA、残余应力σr、残余应变εB，计算弹性模量E，计算脆性跌落系数R、应力降落系数P、软化模量M；(3)、定义脆性跌落系数R相应的脆性指标为B1、应力降落系数P相应的脆性指标为B2、软化模量M相应的脆性指标为B3，并分别进行归一化处理；(4)、通过脆性指标B1、B2、B3，计算得到综合脆性指数Bd；(5)、通过某区域代表性页岩气井地质力学参数对综合脆性指数Bd计算模型进行修正，得到适用于本区域的脆性指数计算公式Bd3；(6)以修正后的脆性指数计算公式Bd3为目标函数，以测井数据中页岩动态弹性模量E、动态泊松比μ、自然伽玛API、断裂韧性K II为自变量，进行函数回归，得出Bd3的多元回归函数Y，其中，式中，a、b、c、n为常数；(7)、应用该回归函数代入不同层段储层测井解释的地质、力学参数，得到不同井段储层的综合脆性指数。按上述技术方案，步骤3中：B2＝(σa-σr)/σa；B3＝1-exp(M/E)；式中，μ-泊松比，σ3-试验加载围压。按上述技术方案，步骤4中：当脆性指标B1、B2、B3中至少一项权重系数大时，将B1、B2、B3分别赋值权重系数α、β、γ，综合脆性指数Bd1＝αB1+βB2+γB3，其中，α+β+γ＝1。按上述技术方案，还可以定义综合脆性指数Bd2＝B1*B2*B3。按上述技术方案，Bd3＝αXB1+βYB2+γZB3，式中αX表示不同的储层中B1的权重系数，βY表示不同的储层中B2的权重系数，γZ表示在不同的储层中B3的权重系数。按上述技术方案，300＜a＜350，1＜b＜2，200＜c＜220，0.1＜n＜0.3，根据不同储层条件选取不同数值。本专利技术所取得的有益效果为：1、本方法能够反映材料在破坏前后抵抗非弹性变形的能力和丧失承载力的情况，同时可实现室内试验数据与水平井段测井解释数据的相互结合，使得该脆性指数具有综合性，适应性较强。2、本方法考虑了页岩破坏的整体情况，综合脆性指数Bd为多因素定量评价指标，可更全面地反映不同压力条件下岩石的脆性破坏过程和特征，可以根据不同的目的选取不同参数，同时可以结合页岩气井测井资料对不同深度、不同围压、不同井段储层的脆性进行分析，突破了传统脆性评价仅考虑矿物成分、静态参数等对脆性影响的局限，具有较好的实用性。附图说明图1为本专利技术提供的三轴压缩试验示意图。图1中，1-三轴压缩实验主体框架；2-管线；3-管线；4-轴向应变测量传感器；5-接收器；6-实验试件；7-径向应变测量传感器；8-数据线；9-实验装置总控制台；10-增压系统；11-计算机。图2为待测页岩岩样的加卸载应力-应变曲线。图3为应力应变示意图。具体实施方式下面结合附图对本专利技术作进一步说明。本实施例提供了一种用于页岩气井储层脆性综合评价的方法，包括如下步骤：(1)利用待测的页岩岩样进行三轴压缩测试实验，获取全应力-应变曲线，其中，根据图1所示的实验装置连接图，连接实验管线和数据传输线，其中，装置为实验的核心部分；连接后管线2、3和数据线8后，安装实验试件6(待测页岩岩样)，实验试件6的安装主要包括径向应变传感器和轴向应变传感器、上下压头的安装，其中，上下压头中装有接收器5，连接密封腔内部的管线，保证数据能够传入计算机11；安装好试件后，放下三轴压力试验装置的密封腔，利用增压系统10加载围压，围压加载至设定值后，稳定围压2分钟后，利用三轴压缩实验主体框架1上部的液压增压泵进行轴向应力加载，实验过程中，轴向应变测量装置4和径向应变测量传感器7记录实验过程中的应变值；利用数据传输线8将数据传输实验装置总控制台9及计算机11。实验过程中应力-应变曲线测至获得页岩的残余强度，记录页岩实验过程中的应力-应变全曲线；(2)如图2所示，根据实验获得的应力-应变全曲线，读取峰值应力σa、峰值应变εA、残余应力σr、残余应变εB，计算弹性模量E，E＝σr/εB (1)根据岩石加载及卸载过程定义了脆性跌落系数R，从峰值强度对应点A到残余强度对应点B变化过程中，计算脆性跌落系数R为：R＝-(εB-εA)/(εM-εA) (2)其中εA、εB可在图2中直接读取，根据广义的胡克定律σ＝Eε+本文档来自技高网...

【技术保护点】
一种用于页岩气井储层脆性综合评价的方法，其特征在于：包括如下步骤：(1)、利用待测的页岩岩样进行三轴压缩测试实验，获取全应力‑应变曲线；(2)、根据实验获得的应力‑应变全曲线，在曲线中读取峰值应力σa、峰值应变εA、残余应力σr、残余应变εB，计算弹性模量E，计算脆性跌落系数R、应力降落系数P、软化模量M；(3)、定义脆性跌落系数R相应的脆性指标为B1、应力降落系数P相应的脆性指标为B2、软化模量M相应的脆性指标为B3，并分别进行归一化处理；(4)、通过脆性指标B1、B2、B3，计算得到综合脆性指数Bd；(5)、通过某区域代表性页岩气井地质力学参数对综合脆性指数Bd计算模型进行修正，得到适用于本区域的脆性指数计算公式Bd3；(6)以修正后的脆性指数计算公式Bd3为目标函数，以测井数据中页岩动态弹性模量E、动态泊松比μ、自然伽玛API、断裂韧性KⅡ为自变量，进行函数回归，得出Bd3的多元回归函数Y，其中，式中，a、b、c、n为常数，；(7)、应用该回归函数代入不同层段储层测井解释的地质、力学参数，得到不同井段储层的综合脆性指数。

【技术特征摘要】
1.一种用于页岩气井储层脆性综合评价的方法，其特征在于：包括如下步骤：(1)、利用待测的页岩岩样进行三轴压缩测试实验，获取全应力-应变曲线；(2)、根据实验获得的应力-应变全曲线，在曲线中读取峰值应力σa、峰值应变εA、残余应力σr、残余应变εB，计算弹性模量E，计算脆性跌落系数R、应力降落系数P、软化模量M；(3)、定义脆性跌落系数R相应的脆性指标为B1、应力降落系数P相应的脆性指标为B2、软化模量M相应的脆性指标为B3，并分别进行归一化处理；(4)、通过脆性指标B1、B2、B3，计算得到综合脆性指数Bd；(5)、通过某区域代表性页岩气井地质力学参数对综合脆性指数Bd计算模型进行修正，得到适用于本区域的脆性指数计算公式Bd3；(6)以修正后的脆性指数计算公式Bd3为目标函数，以测井数据中页岩动态弹性模量E、动态泊松比μ、自然伽玛API、断裂韧性KⅡ为自变量，进行函数回归，得出Bd3的多元回归函数Y，其中，式中，a、b、c、n为常数，；(7)、应用该回归函数代入不同层段储层测井解释的地质、力学参数，得到不同井段储层的综合脆性指数。2.根据权利要求1所述的...

【专利技术属性】
技术研发人员：肖佳林，李远照，吴雷泽，李奎东，陈希茜，王小军，黄进，张辉，刘炜，刘龙，
申请(专利权)人：中国石油化工股份有限公司江汉油田分公司石油工程技术研究院，
类型：发明
国别省市：湖北;42