≯震源区和接收区结构不同情况下体波.PDFVIP

维普资讯 第 35卷 第 2期 地 球 物 理 学 报 V0I．35 N0．2 l992年 3月 ACTA GEOPHYSI~A SINICA M ar．， l992 ≯ 震源区和接收区结构不同情况下体波 合成地震图的矩阵一射线方法 缒 ／ 2 ． (国隶地震局地质研究所，北京 nooo29) 摘 要 本文进一步扩展了合成地震图的广义射线方祛 ． 由于算法中将地球模型划分成震源区、 接收区和地幔区，因此可用于计算震源区与接收区速度结构不同时的体波合成地震图． 新算 法把广义射线理论与传播矩阵理论结合起来，并采用了分层9值模型，从而可以自动包括指定 介质层内所有多次波的棍响．为简便数据准备，设计了射线编码自动生成算法． 本文算法与 Kind扩展的反射率法用相同的地球模型进行了数值结果对比， 结臬表明， 矩阵一射线方法用于计算体波地震图时，不但有足够好的精度，而且还具有速窿快并适于单个 震相研究的优点．在计算 由l0十接收点组成的地震剖面时，其计算时间倪为反射率法的7宁矗 左右． 关罐词 龟堕 垦 ，亡 越些， · 一 ， 引 言 利用合成地震图技术拟合实际的地震波形已被广泛地应用于推断震源参数、地球 内 部结构和震相特征的研究 “ ．虽然地球 内部 (特别是岩石层 内)普遍存在着横向变化及速 度各向异性 ，但在远震体波研究中，人们仍广泛采用分层均匀、各向同性的地球模型． 广义射线方法是计算体波合成地震图的常用方法 “．它具有计算速度快和便于研究 单个震相的优点，但 由于它不能计算包含所有多次反射被的合成地震图，因此在实际应 用中受到一定的限制 ． 特别是在利用波形信息研究地球 内部的细结构时，考虑层内的 多次反射被是十分必l要的．这是因为地球内的间断面常常并不是简单的一级间断面，而 可能是由一个翦数个薄层组成的复杂结构 ，I】．在这种情况下，特别是层阃存在较大的速 度对 比度时，在合成地震图中包括层 内的多次波则变得十分必要． 反射率法 “，”用于计算分层均匀模型合成地震图．它能 自动包括各种层内及层间的混 响，数据准备也比较简单，但计算量较大，这在计算包含高频成分的远震地震图时尤为突 出．用这一方法所得的理论地震图几乎与实际的地震记录一样复杂 ，难 以对图中的震相 作 出解释． 对于远震事件，震源区与接收区的速度结构一般是不同的，经典的广义射线和反射率 奉文 1990年 j月3日收到，1991年 6月 日收到修改稿 维普资讯 地 球 物 理 学 报 。 ～ … ～ “ 一 的篓情况可能用于这种情况．Kind扩展了经典的反射率法，使之适用于简单横向非均句 ． 堂 ：苎，：。荸nfJ结合曾起将来地球模型分为薄层和厚层，并利用稳相近似积分方法将传播 ， 以便 自动包括薄层 内的混晌效应 求被 数在稳相点 ． 由于稳相 … 一萎 附近不能有任何奇点，这一方法不适用于临界区内震相的