破裂影响在钻孔地应力测量中数值模拟.pdfVIP

·300· 中国地球物理2011 (8)地球内部结构及其动力学 破裂影响在钻孔地应力测量中的数值模拟 numericalsimulationoftheinfluenceoffracture The instressmeasurementofborehole 程成石耀霖 ShiYaolin ChengCheng 中国科学院研究生院计算地球动力学重点实验室北京100049 1．引言 地应力的钻孔测量是认识地球内部结构力学性质的重要手段，一直受到科学家的关注和研究的热点问 题。许多国家也开展地应力测量的科研钻探项目，如俄罗斯科拉半岛大陆超深钻、德国的KTB等。尤其 是今年来活跃的KTB钻孔提供的大量地应力数据推动了对引力随深度变化的认识。KTB应力测量结果表 明应力随深度增加，不超过摩擦系数为0．㈨．7的包络线，M．Brudy等人在讨论KTB的实际观测应力的解 释时，指出岩石破裂引起应力释放是引起观测现象的主要原因。 随着计算机运算能力和数值计算技术的发展，本文拟结合数值方法和半解析解方法计算存在的隐伏 断裂的应力释放对应力测量的影响，探讨KTB一类钻孔测到的应力随深度变化特征形成的地质历史演化 成因。 2．研究方法 本研究采用有限元模拟和半解析解相结合的方法，构建二维和三维分层模型，以有限元计算库仑应力 变化过程，以半解析解计算破裂过程引起的应力变化。 在具体讨论过程中，先建岩石圈粘弹性分层模型，上地壳为线弹性体，下地壳和上地幔盖层为Maxwell 体，模型每层的力学性质根据地球物理资料确定，但在模型的每层中随机加入材料参数小的扰动，体现地 层的不均匀性。仅在一个边界给定均匀速度边界条件，考察岩石圈在板块长期运动下的变形和应力演变。 其他侧边界边界水平位移受固定，垂直可以自由滑动；底边界垂向固定，水平可以自由滑动。计算的应力 在累积到一定时间后会超过岩石强度，因此每一定时间步长后要判断在层中一点尤其是材料不均一的位置 是否达到脆性破裂条件。如果达到破裂条件，运用半解析解程序计算破裂引起的应力释放量和引起的位移 变化，与原应力相加减得到破裂后的应力，并考察该应力变化是否会引起其他单元破裂。在不会再其他单 元破裂后，继续进行下一个时间步长应力增长量得计算。如此反复计算，观察地壳应力增长过程和形成的 应力一深度特征，与KTB等钻孔观测进行定性的比较。 以上计算开始在二维模型中进行，程序调通工作后进而推广到三维模型。单位模型中将模型分层迸一 步细化，模型四周边界与二维模型相当，除底面外垂向自由。按照二维模型给出位移速度条件，或使速度 随深度变化。在三轴引力作用下在隐伏断层处发生破裂，得到包含重力作用下应力变化。在二维模型中只 存在逆冲断层存在，在三维模型可以包含走滑分量，得到计算更精确的结果。 3．结果分析 利用上述模型计算了在速度边界条件下破裂在钻孔地应力测量中的影响，计算结果表明，在脆性地壳 中，最大最小主应力差与观测结果大体吻合，说明在有隐伏断层的情况下在达到断裂条件时就会发生破裂 产生应力释放，因此形成目前观测到的地壳应力深度变化特征。 参考文献 【1】M ofthe snesstensorto8b intheKTBscientificdrill forcrustal Brudy．Estimationcomplete depth holes：Implication strength[J]．JournalofGeophysicalResearch，1997，102：18453-18475． s仃eSseS 【2】BezalelHaimson．TrueTriaxialandtheBrittleFractureofRock[J]．Pure