基于构造应力场模拟的页岩气藏裂缝表征与建模方法.pdf

oetroleum Engineering石油工程基于构造应力场模拟的页岩气藏裂缝表征与建模方法刘礼军1周昕任启强李洁辛陈克勇1高俊华陈雪1.成都理工大学能源学院四川成都6 10 0 59；2.中石油川庆钻探工程有限公司页岩气勘探开发项目经理部四川成都6 10 0 51；3.成都理工大学沉积地质研究院四川成都6 10 0 59摘要在构造模型、测井解释力学参数及岩相分布基础上，结合古今构造应力场特征，构建包含地层和断层的研究区三维地质模型，通过弹塑性有限元数值模拟法，进行了构造应力场模拟。以岩石破裂准则为基础，结合岩石应变能理论、能量守恒定律，建立了构造应力作用下裂缝开度、线密度、体密度等参数的定量预测模型，实现裂缝参数空间分布定量预测。通过对某页岩气田构造应力模拟及裂缝定量预测，获取了目的层段古今应力场分布特征，在此基础上，明确了研究区的裂缝参数展布特征，为页岩气田高效开发提供了基础裂缝地质模型。关键词页岩气藏裂缝建模构造应力场中图分类号TE355文献标志码B杨数值模拟文章编号16 7 2-9 3 2 3(2 0 2 3)0 7-0 16 6-0 3随着天然气需求日益增加，页岩气成为天然气的换代资源之一 1-2 。我国具有丰富的页岩气资源，页岩气富集区的研究对我国能源后续供给具有战略意义。天然裂缝的发育程度直接影响了页岩气藏产量 3-4。裂缝为游离态的页岩气提供了储存空间和流动通道，有效提高了页岩的渗透率和孔隙度。因此，页岩气藏裂缝表征对页岩气高效开发具有重要意义。目前，页岩气藏的裂缝表征已成为世界热门课题。但是页岩具有复杂特征，利用常规的油气藏表征和建模手段不能有效描述页岩的裂缝参数和展布。国外对页岩裂缝研究较早，取得了一定的进展。相比于国外，我国页岩气藏条件更加复杂，页岩储层大多经历了多旋回的构造演化，导致页岩气藏断裂系统极为发育。目前页岩气藏开发地质模型过于简化，影响模型准确性，尚无适用于我国1661石油化工建设2 0 2 3.0 7地质特征的成型的页岩气藏裂缝建模方法和开发理论。针对以上问题，基于构造应力场模拟建立页岩气藏裂缝表征与建模方法。此方法依据研究区所处的地质及构造条件，建立地质力学模型，开展构造应力场数值模拟，结合构造应力与裂缝参数间的关系，对研究区裂缝发育情况进行表征，并采用离散裂缝建模技术建立天然裂缝模型，为研究区页岩气高效开发提供支撑。1构造应力场模拟以四川盆地W区块页岩气藏为研究对象进行构造应力场模拟。W区块位于四川盆地西南低褶构造带，主力目的层段为龙马溪组龙一1亚段，整体为深水陆棚相沉积，埋深大多小于3 50 0 m。构造应力场包括古构造应力场和现今构造应力场，其中古构造应力场需结合构造演化、oetroleum Engineering石油工程构造应力场演化和构造裂缝发育史进行确定。构造演化加5 6 MPa和1 2 0 MPa的挤压力。现今模型边界条件：底部史和构造应力场演化史表明，区块构造裂缝可分为4 期：设置约束，南北方向施加4 5 MPa挤压力，东西方向施加第1 期为燕山早期，由于NW挤压形成NE构造体系；第18MPa挤压力，垂向设置重力。2期为燕山晚期，由于NS挤压形成EW构造体系；第3期形成于喜山期，即库EW挤压形成NS构造体系；第4期为印支期，由于龙山北段关闭，由北向南挤压，距离远导致构造变形不大可以忽略。因此，古构造应力场的数值模拟主要是依据燕山期和喜山期的地质构造来构造应力场进行。1.1地质模型依据W区块页岩气田的顶面构造图，并结合区块地震解释成果提取断层及大尺度裂缝几何，在有限元分析软件中建立页岩气田三维地质模型。模型采用四面体网格进行剖分,形成构造应力模拟模型。1.2岩石力学参数岩石力学参数是进行准确有限元数值模拟的前提和保障，其参数的准确性将决定构造应力场模拟的精度和可靠性。以下利用研究区内页岩气井钙质、硅质和泥质含量测试数据对页岩岩相进行了划分，获取目标层段岩相分布模型。而后，依据研究区页岩气井测井解释数据，明确各井点的岩石力学参数空间分布，并将岩石力学参数赋予到对应的地质模型中。采用线弹性八节点单元进行网格划分，共划分4 7 0 5 9 个节点和1 4 0 8 4 5 个单元，每个单元赋予对应的杨氏模量、泊松比和密度。1.3边界条件应力加载方式及边界约束是开展构造应力场研究的基础。基于W区块古今构造应力场特征，确定构造应力模拟的边界条件。根据燕山期和喜山期古应力地图对模型施加力学边界条件，即在模型西南方和东南方分别施0,=45MPa03=18MPa03=56MPa0,=120MPa(a）燕山期和喜山期古构造应力场图1 古今应力场施加边界条件1.4数值模拟结果根据燕山期古构造应力场数值模拟结果（图2），水平最小主应力方向为近NE-SW向，水平最大主应力方向为近NW-SE向，沿断层出现转弯现象，垂向主应力最大,整体表现为Ia型应力状态。平面上最大主应力由南至北方向挤压作用逐渐变小，且应力高值区呈集中在工区中南部位置，受断层和地层影响，最小主应力高值区位于断层附近区域和工区中部；应力分布受成岩微相的影响，主要分布在工区中部及中西部区域。现今构造应力场分布特征与古构造应力场类似，但水平最小主应力方向为近EW向，水平最大主应力方向为近NS向，平面上最大主应力由东至西方向挤压作用逐渐变小。2页岩储层裂缝定量预测2.1裂缝参数计算模型针对泥质页岩，结合库伦一莫尔、格里菲斯等屈服一破坏复合破裂准则以及裂缝力学模型，根据岩石特点，即不论实砂岩是脆性或者是塑性，应变能达到表面能的情况下，就会产生相应的规模破裂。由此推导出挤压应力状态(b）现今构造应力场（a）水平最小主应力方向(b）水平最大主应力方向（c）垂直主应力方向-485E+00.425E00（d）水平最小主应力-.918E+01356E+97.206E+-.3SE+00(e）水平最大主应力(f）垂向主应变图2 燕山期构造应力场模拟结果石油化工建设 2 0 2 3.7 /1 6 7-706557559E119E+0oetroleum Engineering石油工程下,构造裂缝的定量参数计算模型，见式(1)。Du=0i+0+0-2(0,+0,+0)-0,+2(0,+02+0.)0,/(2EJ)2Dv,L,L,sin cos -L,sin -L,cos eD80D.Lsin+L;cose(1)式中：Dv构造裂缝体密度,m/m；Dir一构造裂缝线密度，条/m；b一构造裂缝有效开度,m;岩石破裂应力，MPa；1、2、3 最大、中间、最小有效主应力，MPa;81、8 2、3LI、L、L 沿1、2、3 的单元长度,m;E一杨氏模量,MPa;M一一泊松比；一裂缝破裂角；80岩石承受最大弹性张应变；J一裂缝表面能,J/m。2.2裂缝定量预测结果基于古今构造应力场数值模拟结果，将应力、应变模拟结果代入裂缝参数计算模型，对W区块页岩气田裂缝进行定量预测,预测结果如图3 所示。由图3(a)和(b)可知，裂缝线密度受断层和大裂缝影响较大，在断层边缘地带裂缝线密度较大，但整体裂缝呈现零散发育，集中在工区中部和北东部位。由图3(c)中裂缝开度分布可知,裂缝开度较大的地区为西北方向较集中，由于西南方向喜山期挤压，造成裂缝发育开度较大。由图3(d)可知,裂缝呈北东方向发育，且集中在北东部位。依据构造应力场数值模拟获取的裂缝密度、走向等定量预测结果，利用三维离散裂缝网络建模技术，建立页岩气田裂缝网络模型，为研究区页岩气藏压裂施工和高效开发提供基础模型。3结论(1）基于储层地质和岩石力学参数模型，以及古今TOE-.125B-05-500E-04(a）裂缝线密度一最大、中间、最小主应变；.05构造应力场特征，建立了构造应力场数值模拟方法，并结合岩石应变能理论和能量守恒定律，形成了构造应力作用下裂缝开度、线密度、体密度等参数的定量预测模型，可实现裂缝参数空间分布定量预测。(2）通过页岩气田构造应力数值模拟，获取研究区古今应力场分布特征，结合裂缝定量预测模型对裂缝参数进行计算，明确了研究区裂缝参数分布特征，为页岩气田高效开发提供支撑。参考文献1邹才能，赵群,董大忠，等.页岩气基本特征、主要挑战与未来前景 J.天然气地球科学,2 0 1 7,2 8(1 2)：1 7 8 1-1 7 9 6.2邹才能，张国生，杨智，等.非常规油气概念、特征、潜力及技术一兼论非常规油气地质学 J.石油勘探与开发，2 0 1 3,4 0(0 4)：385-399,454.3徐政语，梁兴，鲁慧丽，等.四川盆地南缘昭通页岩气示范区构造变形特征及页岩气保存条件 J.天然气工业，2 0 1 9,3 9(1 0)：2 2-3 1.4郭旭升，胡东风，魏祥峰，等.四川盆地焦石坝地区页岩裂缝发育主控因素及对产能的影响 J石油与天然气地质,2 0 1 6,3 7(0 6)：799-808.(收稿日期：2 0 2 3-0 5-1 7)005.200E03.4008-0325303550(c）裂缝开度图3 裂缝参数定量预测结果1.04.05(b）裂缝体密度65164.8(d）裂缝走向18516b2145.7.952:05165.91702265822.26592100360168石油化工建设2 0 2 3.0 7