近震源破裂过程反演研究水-中国地震信息网.PDFVIP

维普资讯 -．r 目 柳 堇 (D 辑) 第 33卷 第 5期 SCIENCEINCHINA (SeriesD) 2003年 5月 近震源破裂过程反演研究水 — — I．方法和数字试验 周仕勇 陈晓非料 刘金朝 盖增喜 (北京大学地球物理系，北京 100871) 摘要 利用非线性规划研 究的最新成果，设计 了一种全新的震源破裂过程 的反演方 法，同时反演震源破裂的时间图像和滑动在破裂面上 的分布．该方法与 目前震源过程 破裂反演计算所用 的两类主要方法 (线性矩阵方法和全局搜索寻优方法)相 比，在线性 矩 阵方法对初始模型依赖大、解不稳定和全局搜索寻优方法在高维解空间寻优效率不 高、易出现伪解等方面有较显著改善．数字试验 的测试表明，本方法计算效率高，反演 结果稳定、可靠，可有效应用于实际地震震源过程的反演研究中． 关键 词 震源破裂过程 反演 非线性规划 信赖域 方法 波形拟合 震源破裂过程 的反演 目前多采用波形拟合反演的思想．大多是 由先期研究设定震源位置 和断层方位，将断层面划分成若干子单元 (子断层)，把各子断层滑动产生的地震波场按一定时 间延迟叠加，得到整个震源破裂引起的地震波场，通过与观测到的地震波形拟合，反演计算断 层面上各子单元上 的破裂参数 (如破裂时间、滑动角及位错等)，勾画出震源破裂过程 ． 震源破裂过程反演所依据 的基本方程可表达为 A(r) ‘X =b， (1) 式 中， 为需反演的各子破裂单元破裂延迟时间向量，A()为包含 自变量因子 的格林函数矩阵，b 为由真实地震波形记录构造的观测数据 向量，X 为需反演的各子破裂单元上的位错矢量矩阵．需 指 出的是，格林函数 A()是关于 的非线性 函数并且一般情形下没法给出解析关系式． 基本方程(1)可 以进一步详细写为 A(r)11 A(r)12 … …A(r)l，，l A(r)21 A()22 … …A()2m X2 ● (2) ： ： A()1A()2 … …A() l ljj l 式中列号为子破裂单元编号，m为子单元数，行号为资料点编号，n为资料点数 ． 现在 国际上采用 的种种反演计算方法大体可分为线性和非线性两类 ．线性方法 是通过 假定断层面各子单元的破裂可能发生在预先设定 (通过假定破裂传播方式和传播速度)的多个 2002—06—25收稿 国家 自然科学基金 (批准号4叭34叭O)和国家重点基础研究发展规划项 目(批准号：95—13—05一O4)联合资助 联系人 E—mail：xfchen@pku．edu．ca 维普资讯 第 5期 周仕勇等：近震源破裂过程反演研究 483 子时间窗内，这样格林函数矩阵简化为了不包含破裂延迟时间、由若干子单元和子单元在多个 子时间窗发生的子破裂组成 的常系数矩 阵，基本公式为 AX=b， (3) 其 中A为 由各子破裂单元在预先设定的多个时间窗发生破裂构成的格林函数矩阵，因而与(1) 式不同的是格林函数矩阵 A 不包含破裂延迟时间因子，且由于假定子破裂单元在每一个预先 设定的时间窗中都可能发生了破裂，则A矩阵变为 rtxml维矩阵(其 中，m和 n的意义与(2)式 同；l为预先设定的时间窗个数)；式中x为各