用边界元法模拟同震形变场-中国地震信息网.pdfVIP

维普资讯 Vo1 11．No4 簟ll卷第4期 地 壳 形 变 与 地 震 ． ． 1991年 U 月 CRUSTAL DEFORMATION AND EARTHQuAKE Nov． ． 1991 用边界元法模拟同震形变场 孙 建 中 (国家地震局地震研究所 ) 摘 要 断层非均匀强度分布对地震 的孕育、发生和迁移具有十分重要 的作用。本文简要地阐述了边 界元方法的基本原理，利用该方法摸拟计 算了断层非均匀强度分布的同震形变场 ，并就利用形 变场观测资料探索断层带的强度分布展开 了讨论。 一 、 引 言 大量的物理实验和观测事实证明，地震断层面的强度非均匀分布是地震孕育、发生和迁 移的主要机制之一。这种非均匀强度分布可抽象为若干个线性障碍体将断层分为若干段的情 形，实际上这是对断层面几何非平滑性的物理牲能的描述。由于微观粗糙表面的相互嵌套， 引起宏观强度分布的不均匀。MYATAKEE”认为低强度区的破裂能绕过高强度区继续扩展， 在高强度块体中发生应力集中，从而孕育下次地震。因而断层面强度分布的确定有助于强余 震和复发地震的预测。 安艺敬一 (Aki)试图根据余震的分布 以及体波和面波地震矩的差异 来探求高强度区的存在 。 斯图阿特 (stuart)利用板块运动和历史地震资料，通过对一系 列大地震的拟台给出了圣安德烈斯南闭锁段的强度分布及其可能发震的时间表，这些工作都 是有益的探索。 地表形变场在一定程度上携带了地壳深部 的信息，反映了地壳深部物质流和能量流的交 换以及地质构造体的相互作用，因此有可能根据形变场的形态寻求高强度 区的存在。本文利 用边界元法模拟计算深埋高强度区对地表同震形变场的影响，并就利用实际形变场观测资料 确定断层的强度分布进行了探讨 。 二、边界元基本理论 边界元方法是近年来得 到迅速发展和应用的一种数疽计算方法，它是将边界广义位移和 广义力作为独立变量，以满足场方程的奇异源函数为加权函数，因而是一种特殊格式的加权 余量法。 奇异源函数的引入使边 界元具有独特的优点。 首先， 问题的解只涉及研究体的边 本文 1990年 n 月P4日骏 到 枣本文系国家地震局地霹研究所，硕士论文柏一部分 维普资讯 4期 孙建中：用边界元法摸拟 同震形变场 3l 界，这种降维特性减少了方程数 目和数据准备最，从而节省了计算机时间，其次，由于边界元 法是建立在 问题的物理场方程和边界条件基础之上的，无需事先寻找任何泛函，因而可以求 解经典数值方法所无法解决的有关无限域一类的问题J再次，由于引入算子基本解，边界元 具有解析与离散相结合的特点，因而具有较高的精度 。这些特点使其适合于研究地学问题。 对复杂的物理场方程的数值解 的研究，产生 了许 多数值方法，但其解不外乎满足两个条 件，一是物理场方程，另一个是边界条件 。边界元方法就是基于这两点提出来的。设一体积 为 V，边界为r的弹性体受到两组外力的作用，对其运用敦度定理和贝蒂互易定理，并考虑 物理场方程，可 以导出边界积分方程为； c·U·l(P)=IU t(P，0)。(口)a8(Q) J ， 一 IT l(P，口)地(口)d (口) J +IU·l(P，g)，(‘g)d (g)