渭北地区某区块二维构造应力场模拟.docxVIP

第三章 渭北地区某区块二维构造应力场模拟3.1渭北地区某区块概况渭北地区位于鄂尔多斯盆地南部，包括渭北隆起及其临区，大地构造位置上处于秦岭造山带与鄂尔多斯稳定地块过渡地带。其演化、改造与其南部秦岭造山带和渭北盆地形成演化息息相关。渭北隆起现今构造面貌为一近东西向展布，略向东呈弧形凸出，并具有相对隆起的地貌特征。渭北隆起以渭河盆地北缘断裂为南界，北部与陕北斜坡呈过渡关系。渭北地区地处多个性质不同的大地构造单元的复合部位。中新生代以来经历了多期次构造变动，遭受了强烈的剥蚀改造，致使渭北地区的构造特征变得复杂。研究区域全图如下：图3-1 区域长3断面全图渭北2井区断裂相对发育，高楼坬-旬邑县北二级断裂穿越三维区的东南部，区内延伸长度3km。其次，区内发育一些三级或更次一级的小断裂，其中西北部发育两条延伸长度大于10km的断层，断层走向北东东，倾向分别为北西和南东，断面倾角大于700，为拉张应力条件下产生的正断层。相干及倾角切片分析认为，本区裂缝不发育。图3-2 渭北2井区长3断裂要素表图3-3长3研究区域垂向剖面图3.2最小二乘法反演区域地应力3.2.1 区域地应力资料整理经过对已有的地应力数据结果整理和计算，目标井点的资料整理如下：井号层位深度/m水平最大主应力/(MPa)水平最小主应力/(MPa)最大水平主应力平均方向/MPa/MPa/MPa渭北46(C点)长355015.8314.2N25°E15.4214.60.71渭北2(B点)长354512.0510.58N45.75°E11.6810.950.636渭北39（E点）长349014.3513.18N5°E14.0613.470.51渭北20（D点）长3505N60.25°E渭北33（A点）长3480N50.5°E3.2.2 力学模型的建立首先确定地质模型，地质模型是数值模拟的基础。对平面二维模型，依据渭北地区断裂系统图以及渭北地质构造图（参见上图），模型东西向（X向）取20km、南北向（Y向）取20km的矩形区域。本模型主要针对渭北长3区块进行研究，其中为了消除边界效应而增加的隔离面。边界模型取东西向（X向）取30km、南北向（Y向）取30km的矩形区域。比例尺为1:10000，本文研究区域为为渭北地区黄陵县下面的长3地区一部分区块，通过对已有的区块地质资料整理，判断该区块的最大主应力的平均方向在N53.5°E，所以分析区域如下图所示的红色部分所示 图3-4 渭北长3研究区域 图3-5 区域地质模型断裂以逆断层性质为主。建立模型时，根据地质资料建立F1，F2，F3，F4，F5，F6六个断层进行分析，根据地质资料分析，对断层带的岩石力学参数进行调整，该区域的断层带的弹性模量为区域弹性模量的40%左右，泊松比比区域泊松比的108%左右。整理如下：岩石力学参数弹性模量/MPa泊松比区域与边界130000.23断层52000.2484对模型进行网格划分，研究区域采用四节点等参矩形单元，边界采用三角形单元。共划分6222个单元，5398个节点。渭北区域模型共划分4412个单元，断裂带共划分1810个单元，有限单元网格模型如下图：图3-6 力学模型网格3.3.3 最小二乘法反演过程在实际工作中对二维地应力反演具体操作步骤如下：（1）在同一地层选择不在同一直线上的3个点作为反演测试点，用试验的方法测得其最大最小主应力大小与方向。再根据公式算得这3个点的σx, σy,τxy大小。如下面给出3个不同方向节点的应力值。（2）我们设有m个边界载荷参数（该步骤中边界载荷参数越多，计算越精确），它们是i=1,2,3….m。我们用abaqus对m个基本载荷工况进行分析（也就是计算在假设的100MPa外部压力条件下工况内的应力场分布，假设的100MPa也可以为其他值）。从而可得m组基本压力场σ(i)x, σ(i)y,τ(i)xy。图3-7加载工况下表为m组不同的压力场下节点46、39、2的压力分布情况，此算例中给出了16组载荷工况的结果。（3）若用λi（i=0,1,2…m）来表示某个基本载荷工况的放大系数，即载荷参数。将这些基本工况应力乘以相应的放大系数，并相加，得到在m个载荷工况共同作用下的应力场的表达式引入矢量和矩阵上式可写成如下的矩阵形式其中σ(x)为实验所测的数据为已知数据，J矩阵为abaqus所计算得出数据。根据上述公式可计算得出载荷边界系数矩阵λ。将载荷边界系数乘以100MPa即得出应力场的边界条件。（4）根据反演计算出的边界条件计算地层整体应力场分布，与实验所测得的3个点的应力值进行核对分析，若精度不够，可增加边界工况个数，直至计算结果与实验值吻合为止。本文不断调试，最终选择边况个数为16个，提取的矩阵系数如下3.3.4 初始地应力模型的选取已知的应力参数选取实验测得的主应力值和方位角。