基阶及高阶面波频散和接收函数的联合反演.pdfVIP

基阶及高阶面波频散和接收函数的联合反演 沈玉松徐果明倪四道朱良保 中国科学技术大学地球与空间科学学院安徽合肥230026 面波频散曲线含有丰富的地下层状结构信息，对地下横波速度结构的绝对值和连续变化特别敏感，而且在同… 反演深度上，高阶面波比基阶面波有更高的分辨率，能够较好地分辨出地壳和上地幔一些局部性的较薄的速度变 化层。而接收函数含有接收器附近地下层状结构的信息，对地下横波速度结构的局部突变特别敏感。所以联合基 阶及高阶面波频散曲线和接收函数的反演有助于双方的优势互补。 基阶及高阶面波频散曲线和接收函数联合反演的具体过程。由于接收函数反应的是接收器附近的地下结构信 息，而我们提取的频散曲线反应的是射线经过的大圆路径下的地下结构信息，因此我们首先必须使用提取的频散 曲线通过层析成像的方法得到局部地区地下结构的频散曲线。我们把面波频散曲线和接收函数的正演表示为： 【儿，…，虬，y】=Fix]，其中yl，．一，y，，是分别是1到n阶面波频散数据向量，Y是接收函数数据向量，x是模型 空间向量，FI1是把模型空间向量映射到数据空间向量的非线性算子，为了求解模型空间向量z，我们使用阻尼 最小二乘法把上述非线性问题的反演化为线性的迭代过程：05,=vF．况。，x。=x。+国¨，因此可以通过迭代 ’ ” ” lJ¨ 的方法由初始模型得到反演结果。与单一反演相比，由于联合反演中面波频散曲线数据和接收函数数据有不同的 物理单位和数据点数，因此在反演中必须对数据进行归一化处理，否则，反演结果很有可能被单方蕊的数据所支 配。为此我们定义一个联合反演预测误差： ， +号封率 ￡：酬雩M 其中P．和M分别为i阶频散的权重因子和数据点数，P和JV为接收函数的权重因子和 数据点数，匕、％和yf、只分别是频散曲线和接收函数的偏导数矩阵和数据残差，盯，．是对应数据的协方差， P。+…+P。+P=l。反演时我们可根据各类观测数据的误差大小，在肛l之间调节各个权重因子的大小，以便 合理地反应各类观测数据对反演结果的影响。当其中一个因子权重为l，其它为O时，就退化为单一的反演。 为了验证基阶及高阶面波频散和接收函数联合反演的有效性，我们在理论上进行了一些反演模拟。即先给定 一个模型，用正演计算出3阶频散曲线和接收函数，并假定其为实际观测的频散曲线和接收函数，然后用上面介 绍的方法进行反演拟合，观察反演结果和给定模型之间的偏差。考虑到反演问题的非线性性，我们在一个给定范 围的模型空间中随机产生了50个初始模型，分别反演后取误差E为最小的反演结果为最终结果。通过调节权重 因子，发现反演结果和已知模型很吻合，并且有效地改善了在单一接收函数反演中的速度深度的耦合问题和在面 波频散反演中的界面模糊问题。 (国家自然科学基金资助项 参考文献 inversionofreceiverfunctionandsurface [1]J．Juliaet．a1．Geophys．J．Int．(2000)143，99一112，Joint wave observations dispersion · · 521