刘中春;吕心瑞;李玉坤;张辉
【摘 要】断层附近区域地应力方向相对于远离断层区域存在明显变化.前人利用实验及有限元模拟等手段研究了断层对其附近区域地应力方向的影响规律,但是缺少从力学机理角度解释这些规律的研究.研究总结了单一断层附近最大水平主应力方向变化规律,利用孔口问题弹性力学解析解揭示了这种规律的力学机理,并以某单元为例,研究了多断层区域最大主应力方向沿断层走向上的变化.结果表明,断层端部区域最大水平主应力朝着断层走向偏转,应用孔口问题弹性力学求解可以获得理论证明,塔河某单元的地应力计算结果也验证了这种规律. 【期刊名称】《石油与天然气地质》
【年(卷),期】2016(037)003
【总页数】7页(P387-393)
【关键词】椭圆孔口问题;地应力;断层;塔河油田
【作 者】刘中春;吕心瑞;李玉坤;张辉
【作者单位】中国石化石油勘探开发研究院,北京100083;中国石化海相油气藏开发重点实验室,北京100083;中国石化石油勘探开发研究院,北京100083;中国石化海相油气藏开发重点实验室,北京100083;中国石油大学(华东)储运与建筑工程学院,山东青岛266580;中国石油大学(华东)储运与建筑工程学院,山东青岛266580
【正文语种】中 文
【中图分类】TE122.1
地壳中较大范围内地应力场方向通常趋于一致，但在断层存在的周边区域水平最大主应力方向会发生较大变化。前人针对断层对周边地应力场方向影响规律做了大量工作。苏生瑞通过对国内外大量实测地应力资料的分析，分别研究了活动断裂和非活动断裂对地应力场的影响[1-2]。孙礼健等用有限元方法模拟研究了断层规模、岩体和断层弹性模量、断层走向与区域最大水平主应力的夹角、边界应力比、断层几何形态等因素影响下断层附近应力场分布规律[3]。彭学军等进行了断层区域应力场对围岩稳定性影响研究[4]。以往的研究多
集中在断层对地应力场方向影响规律的描述上，未从力学角度对这种影响机理做深入研究。
本文在前人研究基础上，用有限元软件建立单个断层平面计算模型，进一步总结了断层对水平最大主应力方向的影响规律。运用椭圆孔口问题的弹性力学复变函数解，分析了断层附近水平最大主应力方向变化机理。最后针对塔河某区块地应力场反演结果，分析了多断层区域最大主应力方向沿断层走向上的变化规律。
由于实施实际地应力方向测量代价高昂，针对断层附近地应力方向的测量数据较少，无法对其规律进行系统总结。而且实际地层中影响地应力方向因素众多，单靠实测值无法排除其他因素的干扰[5]。因此，在对已有资料分析的基础上，利用数值模拟手段研究断层附近地应力方向的变化规律十分必要。为了逐步揭示断层对地应力场方向影响规律，本文先从单断层情况入手进行讨论。
由于地层受上覆岩层压力，可以认为在水平挤压力下其垂直方向应变为零，因此为了更直观地研究断层附近地应力方向变化规律，采用平面应变模型建模。建立100 m×100 m平面方型地层区块模型，在区块中间建立延伸长度为40 m，最大断距为2 m的椭圆型断层，断
层走向与区块边界呈45°夹角(图1a)。
为了使计算结果更有针对性，本文在塔河油田缝洞型油藏条件下进行研究[6-8]。根据塔河油田某区块地应力测试结果，该区块平均最大水平主应力为120 MPa，平均最小水平主应力为90 MPa。因此，在模型左右边界施加最大水平主应力120 MPa，上边界施加最小水平主应力90 MPa。为了使模型能够求解，将底部边界上下方向位移进行约束，建立的ANSYS几何模型及施加的边界条件如图1b所示。该区块油藏岩心实测周围岩体弹性模量为45 GPa，泊松比为0.22。断层物质多为填充物，物性较软，所以断层弹性模量设为5 GPa，泊松比设为0.35。选择平面应变单元进行网格划分，其中断层端部采用较精细的网格(图1c)。该模型假设断层只产生微小位移，断层活动产生的次生应力场可以忽略。
垂直主应力主要受上覆岩层压力控制，一般为竖直方向，断层的存在对垂直主应力方向几乎没有影响。最小水平主应力与最大水平主应力垂直，确定了最大水平主应力方向也就确定了最小水平主应力方向。因此，本文主要研究断层周围最大水平主应力的变化规律。
对上述模型进行求解，得到最大水平主应力大小和方向云图(图2)。
从与最大水平主应力夹角为45°的断层模型的计算结果可以看出，断层内部以及断层边缘最大水平主应力较远离断层区域偏小，但是断层端部区域最大水平主应力出现明显应力集中现象。靠近断层端部区域最大水平主应力方向出现偏转，偏转趋势为沿着断层向断层端部偏转。断层内部最大水平主应力方向总是垂直于断层走向。
为了研究断层走向与最大水平主应力夹角是否影响断层附近区域最大水平主应力偏转程度，按照上述建模方法，分别建立断层与最大水平主应力方向夹角为30°,60°以及75°的计算模型。求解后提取模型应力场方向，图3为断层走向与最大水平主应力方向不同夹角计算模型应力场方向结果。
从对比图3a至图3d结果可知，断层走向与最大水平主应力方向夹角α对断层附近应力场有显著影响：α越大，断层对附近地应力场方向影响范围越大，其附近最大水平主应力方向朝断层端部偏转越明显。
本文总结的上述规律已被很多学者的研究所证实[9-12]，如何解释这些规律是本文的研究重点。如果不考虑断层的填充情况，可以将断层在地层中的受力情况用无限大受压平板椭圆孔口平面应变模型来描述。对于椭圆孔口问题，弹性力学给出了复变函数解答。
如图4a所示的椭圆孔口平板受力模型，对其采用保角变换：
u魅
(1)
式中：z为z平面复数；ω为变换函数；ζ为ζ平面复数。
把弹性体在z平面上所占的区域变换为ζ平面上的区域。这样，在ζ平面上令
(2)
式中：ρ为椭圆孔边缘一点和椭圆中心连线的长度，m；θ为椭圆孔边缘一点和椭圆中心连线与椭圆长轴的夹角，(°)。ρ和θ构成极坐标系。
则ρ，θ与x，y转换关系为：zeidan
(3)
(4)
式中：x和y分别为以椭圆长、短轴为直角坐标系中的横、纵坐标；，其中a为椭圆长半轴长
度，b为椭圆短半轴长度，m；。
最终求得应力分量复变函数表达式为：
(5)
式中：σρ为极坐标中径向应力，Pa；σθ为极坐标中环向应力，Pa；τρθ极坐标系中切应力，Pa；q为板两端拉应力，Pa；Re(z)为复数z的实部；α为椭圆长轴与拉应力夹角,(°)。
隧道定位
图4 椭圆孔口问题理论解分析Fig.4 Theoretical solution analysis of the elliptic orifice problema.椭圆孔口平板受力模型;b.无限大平面椭圆孔口四周受压模型;c.孔边一点应力状态
将公式(2)代入公式(3)—公式(5)，分开实部和虚部，可得σρ，σθ和τρθ的表达式，但是该过程复杂不易求解。根据保角变换可知，在椭圆孔边，ρ方向为垂直于椭圆曲线切线方向，θ方向为椭圆曲线切线方向。则在孔边界上有σρ=0，ζ=cosθ+isinθ=eiθ，因此由(5)式得：
(7)
磨具制造
式中：σ代指eiθ。
传送侦测
公式(7)为单向受拉或受压无限大平板椭圆孔口环向应力计算公式，根据叠加原理可求得如图4b四周受压无限大平板椭圆孔口环向应力计算公式：
(8)
式中：q1,q2分别为两个方向拉应力，Pa。
从图4c孔边一点应力状态可以看出，在椭圆孔口边缘τρθ=0，所以σρ和σθ为两个主应力，由于σρ=0，所以σθ是椭圆孔口最主要的应力。上文模拟得知，断层附近应力场方向沿着断层边缘向端部偏转，而椭圆孔口σθ方向正是椭圆边界切线方向。下面通过考察椭圆孔口σθ变化规律来研究σθ对断层附近应力场方向的影响程度。
2.2 椭圆孔口环向应力分析
对于图4所示力学模型，如果不考虑椭圆断层填充物的影响，根据公式(8)可计算出断层边缘环向应力分布。将q1=120 MPa，q2=90 MPa，m=0.9，α=π/2代入公式(7)，计算不同θ下σθ的值，并绘制出断层边缘一周环向应力曲线(图5)。
古马隆树脂从图5可以看出，如果不考虑椭圆断层填充物的影响，断层边缘σθ变化非常剧烈。当θ=0°(360°)或180°时，σθ达到边界应力的近40倍；当θ在40°～160°以及220°～340°范围内时，σθ不到边界应力值的一半甚至接近于0。由此可知，断层端部区域σθ非常大，导致最大水平主应力方向趋于与σθ方向一致，即沿着断层走向方向偏转。而远离端部边界处σθ很小，不足以改变最大水平主应力方向，因此远离断层端部边缘最大水平主应力几乎没有偏转。
豫公网安备41160202000603
站长QQ:729038198
关于我们
投诉建议