某气场储气砂层组地应力场模拟

吴 屹 张国强

（成都理工大学地质灾害防治与地质环境保护国家重点实验室，四川 成都610059）

0 引言

地应力是引起采矿、水利水电、土木建筑、铁道、公路、军事和其他各种地下或露天岩土开挖工程变形和破坏的根本作用力，是确定工程岩土力学属性，进行围岩稳定性分析，实现岩土工程开挖设计和决策科学化的必要前提。地应力状态对地震预报、区域地壳稳定性评价、油田油井的稳定性、核废料的储存、岩爆、煤和瓦斯突出的研究以及地球动力学的研究等也具有重要意义。

本文以某地区储气层为例，拟利用三维有限元方法，在充分考虑研究区地质、构造等背景，同时充分考虑模拟区几何模型的形状及大小、区域构造格架、边界条件下对储气砂层地应力场进行数值模拟，并对地应力分布规律进行分析。

研究区域内JP气藏层中砂层组是储气层，其勘探与开采过程中，砂层组及其周围岩层的应力分布是确定勘探与开采方法与确保开采效果的重要参数，因此，地应力场的研究显得尤为重要。

根据圣维南原理，分布于弹性体上一小块面积（或体积）内的荷载所引起的物体中的应力，在离荷载作用区稍远的地方，基本上只同荷载的合力和合力矩有关，荷载的具体分布只影响荷载作用区附近的应力分布。作为厚度仅有几米到几十米的砂层，相对于研究区域，其引起的荷载变化整体而言可视作微小，其对整个地区现今地应力场的影响在远离储气砂层组的地方可以忽略。因此，在研究储气砂层的应力分布时，可以将现今地应力场反演得到的结果作为边界条件（远离砂层组的地方），加在考虑了砂层的模型之上，利用数值模拟的方法，得到较为准确的砂层组及其附近岩体的应力分布。遵循这一思路，本文将利用已知现今地应力场计算的结果，得到JP37+8砂层组地应力模拟的边界条件，建立包含砂层组的有限元模型，进而对其地应力场进行数值模拟。

1 JP37+8砂层组地应力场模拟

1.1 模型建立

根据之前对场区的地应力的反演恢复作为依据对储气砂层组的地应力分布特征进行模拟，但因砂体分布不规律，且相对于整个区域，砂体厚度太薄（场区长50000米，宽40000，厚约3000米，砂体厚度5～20米），限于计算机硬件及软件的约束，模型在划分网格过程中无法成功。因此本文对场区某局部含储气砂层区域进行建模分析。模型如图1，模型长16500米，宽9000米，厚度为2800米，砂体厚度为5米，在模型中的展布如图2。

1.2 地质情况和材料参数

JP37+8层的砂体分布区域及厚度范围如图1、图2所示，砂体主要岩性为砂岩，砂体周围为泥岩。

砂体材料计算参数取值如表1所示：

表1 JP37+8层计算物理力学参数表

1.3 单元划分和网格类型

建模分网时采用Ansys中的Solid45单元，建模后采用该Solid45单元进行分网，利用软件接口将分网后的模型方便地转换为FLAC3d可以识别的模型。整个计算域内共有112877个节点，623450个实体单元，研究区域的单元模型见图3与图4所示。

1.4 数值模拟

根据已得到的现今地应力场模拟的结果进行处理，得到模型在X轴及Y轴平面上应力边界条件施加如图5、图6，以此为条件计算得到JP37+8砂层组的应力场分布。

2 JP37+8计算结果分析

根据该地区已知地应力场的反演计算结果，为目标层JP37+8砂体采用小模型进行应力分布分析提供了边界条件，通过计算得到JP37+8砂层组的应力场分布，下面对计算结果进行分析，为满足加砂压裂储层改造时优化设计方案的要求，对JP37+8计算结果中的砂岩及其与周围泥岩的应力分布情况进行分析讨论。

2.1 垂直向应力

图7、图8分别为JP37+8砂层的垂直向应力的内部剖面云图与垂直向应力的表面云图。

2.2 最大水平应力

图9 、图10分别为JP37+8最大水平应力的内部垂直向与水平向剖面云图及表面云图。

2.3 最小水平应力

图11、图12分别为JP37+8最小水平应力的内部垂直向与水平向剖面云图及表面的云图。

通过计算得到JP37+8层目标区域地应力场垂直应力分布的特征，从云图上均可见该区域的垂直应力整体上呈现由浅向深逐渐增大的特征，砂体的存在对垂直应力的分布有一定的影响，从以上各种云图都可以看出砂体上的垂直应力都较周围泥岩的垂直应力小。在局部存在应力集中现象。

该区域的最大与最小水平应力整体上呈现由浅向深逐渐增大的特征，但是同时从各云图可以看出 砂体以及周围的最大与最小水平应力明显要比相同深度没有砂体分布的区域要小。

3 JP37+8砂层对地应力分布的影响

为研究砂体对目标地层地应力分布的影响，分析不含砂体地层与含砂体地层的的垂直、最大水平、最小水平地应力的变化，对JP37+8层如图13所示进行数据取样分析，得到砂体对地应力分布的影响如图14所示。

由图14可以看出JP37+8应力场的泥岩到砂岩中部的应力分布变化特征，由泥岩向砂体变化，垂直应力与水平应力在边界地段均发生急剧的衰减，水平应力（最大与最小）总体上呈现减小，然后在砂体中趋于平衡，垂直应力有一减小-增大-减小-平衡的过程，总体上砂体内的应力小于砂体外围泥岩上的应力。水平应力减小的量值较大，而垂直应力在砂体中部比砂体外围泥岩小约5-7MPa。

4 结论

本文综合运用Ansys与FLAC3D软件，通过反演计算得到研究区域合理可靠的现今地应力场，利用该研究成果确定了该研究地区JP37+8模型的边界条件，根据地质模型建立数值模型，对含气砂层中的砂体及其周围地层的地应力场进行数值模拟，得到了砂层及其周围岩层的垂直应力与最大、最小水平应力的方向与大小等分布特征，发现了岩性及地形对地应力分布的影响。为该地区的气藏勘探与开采提供数据支持。

［1］李志明,张金珠.地应力与油气勘探开发[M].北京：石油工业出版社,1997.

［2］周文.裂缝性油气储集层评价方法[M].成都：四川科学技术出版社,1998.

［3］李志明.地应力与油层改造方案[J].石油钻采工艺,1998,20(6)：47-52.

［4］徐向荣,马利成,唐汝众.地应力及其在致密砂岩气藏压裂开发中的应用[J].钻采工艺,2000,23(6)：17-21.

［5］余琪祥,吴爱军,李善良,等.地应力剖面在储层开发中的应用[J].石油与天然气地质,2001,22(1)：42-47.

［6］周文,闫长辉,王世泽,等.油气藏现今地应力场评价方法及应用[M].北京：地质出版社,2007.