面波有限频率法3D灵敏度核函数的计算一.pdfVIP

第34卷 第3期 成都理工大学学报 (自然科学版) VbI．34 No．3 2007年 6月 JOURNALOFCHENGDUUNIVERSITYOFTECHNOLOGY(Science＆TechnologyEdition) Jun．2007 [文章编号]1671—9727(2007)03—0340—08 面波有限频率法 3一D灵敏度核函数的计算 刘 魁 程先琼 朱介寿 孔春玉 (成都理工大学信息工程学院，成都 610059) [摘要]高分辨率面波成像应克服射线理论的高频近似假设，考虑有限频率效应。该文给 出了 一 种计算面波有限频率3-D灵敏度核函数的方法。该方法可直接计算球壳层状地球模型中面 波振型特征值和特征函数，而不再需要采用球化平变换来近似计算球壳层状地球模型中面波 频散参数。作者计算了相位延迟、振幅扰动和群到时延迟的3一D灵敏度核函数，并对面波有限 频率效应进行了讨论。分析结果显示：面波相位延迟的3-D灵敏度核函数同2一D体波走时的 核函数形状类似，说明面波具有二维传播特性。对相速度数据而言，Fresnel带方法适用于低 频面波的反演，射线方法适用于高频面波的反演，而对于中等频率的面波反演应采用有限频率 方法。频率的高低是相对于所要取得的分辨率而言的。在反演地壳上地幔细结构方面，群频 散数据更具有优越性 ；然而，群到时延迟核函数 的极强的旁瓣可能使得反演变得不稳定。 [关键词]有限频率；3-D灵敏度核函数；Frcsncl带方法；射线理论 [分类号]P631．443 [文献标识码]A 传统的面波分析方法如面波频散反演基于射 像必须克服高频近似的假设 ，考虑有限频率效应 ， 线理论 ，射线理论采用高频近似 ，从而可采用大圆 把由于介质非均匀性而产生的散射波考虑在 内。 定律 。该理论认为 ：当面波穿过足够平滑的介质 Nolet等 认为在 弱散射 的限制条件下 ，引入 时，实测数据仅受连接震源和接收点的大圆下方 Born近似，采用射线理论可以超越射线理论本身 的构造的影响。足够平滑意味着非均匀性的尺度 的局限性 ，即超越无限频率的假设 ，而只需要采用 要大于面波的波长和 Fresne／带 的宽度 ¨】_』．L引。这 相对弱得多的单散射近似假设 ]，从而利用射线 个假设限制 了射线理论所能取得的分辨率，同时 理论既可以考虑有 限频率效应 ，又具有较高的计 也在实际中造成 了悖论 ：基于射线理论 的高分辨 算效率 。这个观点 由Hung等 经过数值模拟 率面波层析成像首先假设介质非均匀性 的尺度要 证实。事实上，关于有限频率效应的讨论许多都 大于面波的波长和 Fresnel带的宽度，而成像结 隐含地利用 了这一观点。Dahlen等 u基于这个 果却显示出非均匀性的大小要远小于假设的尺 观点，利用体波的WKBJ理论给出了一个完整的 度Ea3。同射线理论类似，面波的WKBJ理论 】．_ 体波有限频率理论，并引入旁轴近似极大地提高 同样采用高频近似 ，要求介质足够平滑。射线理 了计算效率。Meier等 。采用面波的WKBJ理 论和 WKBJ理论在焦散处将不再适用，Tromp和 论来考虑有限频率效应对面波 的影响。Yoshiza— DahlenE6基于 Maslov理论 7¨给出了考虑焦散效 wa¨l。基于面波的WKBJ理论给出了一个二维的 应的面波 Green函数 。 灵敏度核函数的讨论 ，并为提高计算效率给 出了 摆在人们面前 的问题是 ，高分辨率的面波成 旁轴近似条件下核 函数的计算方法。Zhao等 ] [收稿 日