首页 理论教育储层裂缝建模数值模拟成功!

储层裂缝建模数值模拟成功!

【摘要】:本次研究的构造应力场模拟流程如图31。图34 模型边界条件的设置3.3.4.5 数值模拟结果根据以上模拟方法和步骤,对研究区的构造应力场进行了模拟分析。

3.3.4.1 构造应力场数值模拟流程

由于角点网格在储层地质三维建模中应用最为广泛,且处理断层和层面的能力强,对应的建模过程也非常灵活,存储方式也非常节约,所以角点网格模型在储层建模和油藏模拟中被广泛使用。本次构造应力场数值模拟研究,也开发了基于角点网格的构造应力场模拟流程。基于角点网格的构造应力场数值模拟,不仅可以直接利用现有的大量角点网格地质模型,其模拟结果更可以非常方便地导入后期油藏模拟程序。本次研究的构造应力场模拟流程如图31。首先根据井位、层位和平面构造图,建立储层层面模型和断层模型,进而构建储层构造模型,然后根据岩石力学参数,建立角点网格属性模型,给每一个角点网格赋岩石力学参数,之后开始对角点网格进行转化,使角点网格转化为有限元网格形式,进而获取储层角点有限元模型。紧接着,可以基于静水有限元模型进行数值模拟。之后根据FMI分析的地应力分布结果,给模拟边界条件一个试探解。对比数值模拟结果和井点结果,分析构造应力场模拟的效果。如果模拟结果和井点结果差别较大,则返回修改模型边界条件和模型岩石力学参数,否则输出模拟结果。

图31 角点有限元构造应力模拟流程图

3.3.4.2 研究区块角点有限元模型

基于前面FRESS建立的研究区构造模型,通过角点网格向有限元网格转化程序,将研究区的角点网格模型转化为适用于有限元模拟的角点有限元模型。其核心步骤是角点网格模型中断层的处理。让角点网格模型中的断层,在角点有限元模型中成为接触边界。

由于角点网格模型,直接从角点网格模型转化而来,所以不需要对模型进行网格划分。并且模型边界可以通过直接拾取角点网格的壁面获得,也不需要增加额外的边框单元。本次模拟所获得的角点有限元网格模型共划分了14000个角点有限单元(图32)。

图32 角点有限元模型

角点有限元模型不仅可以进行常规的有限元模型建模,还有很强的断层处理能力。断层处理效果和角点网格完全相同(图33)。

图33 角点有限元断层处理效果

3.3.4.3 岩石力学参数设置

根据研究区的数字阵列声波测井资料和岩心测试资料,将模拟角点有限元模型的岩石力学参数设置如表7。

表7 岩石力学参数设置

3.3.4.4 边界条件

模型Z方向为深度方向,X轴正方向指向东,Y轴正方向指向北。根据FMI井壁崩落和压裂缝产状分析,研究区应力主方向为东西方向。约束应变边界条件加载在东西边界角点有限元网格壁面上。约束位移方向和主应力方向相同。模型加载重力体载荷;模型底部设定X-Y两个自由度,Z向位移为0;模型顶部加载上层重力载荷,如图34。

图34 模型边界条件的设置

3.3.4.5 数值模拟结果

根据以上模拟方法和步骤,对研究区的构造应力场进行了模拟分析。模拟结果包括应力张量和应变张量分布,以及米塞斯应力分布图。部分模拟结果见图35-图43。

图35 应力场 xσ分量平面图

图36 应力场 yσ分量平面图

图37 应力场 zσ分量平面图

图38 应力场 xyσ分量平面图

图39 应力场 xzσ分量平面图

图40 应力场 yzσ分量平面图

图41 米塞斯应力分布图

图42 压力分布图

图43 第三不变量分布图

3.3.4.6 构造应力场分布规律

研究区的构造应力场数值模拟结果表明,研究区的应力场分布特征和规律为:

(1)研究区应力分布以中间断裂带为中心分为东西两个部分,断裂带的应力场幅度最高,东西两个区的应力分布较低。

(2)断层尖端应力场奇异性。研究区 xσ和 yσ在断层尖端出现局部高值,而zσ值在断层尖端没有奇异性。

(3)剪应力瓣状分布。研究区 xyσ和 xzσ剪应力,在断层尖端呈现花瓣状分布。 一般每个裂缝尖端,以裂缝延伸线为中轴,左右各一个应力瓣。

(4)应力分布不对称性。断层上下两盘的应力不变量分布呈现不对称性,上盘的压力分布大于下盘,下盘的第三不变量分布大于上盘。在断层两盘总有一侧的应力较附近应力场幅值低。

(5)在两条断层交错处,米塞斯应力的局部极大值由断层尖端移动到断层交错区域。

(6)断层尖端应力场奇异可形成局部高压,断层附近应力场幅值较大,易形成次生裂缝。