区域横向非均匀地球模型地震波传播的富里叶伪谱法模拟.pdfVIP

区域横向非均匀地球模型地震波传播的富里叶伪谱法模拟 王彦宾罗扬陈明雅刘启明 陈晓非 北京大学地球与空间科学学院地球物理学系，北京，100871 1．引言 不断丰富的全球地震观测资料为地球深部精细结构研究提供了条件。和仅仅利用不同震相的到时信息相比， 同对考虑振幅和波形信息将会提高结构反演的精度。目前的研究主要是利用单一震相的异常来反演特定区域(如 核幔边界)的非均匀结构，对于更多的其它震相所提供的信息没有考虑。本研究中，我们试图通过全波场模拟 研究地球深部非均匀结构对全波场的影响，从而探讨未来利用更多震相反演深部非均匀结构的可行性。我们利 用错格富里叶伪谱法(Kosloff订aL，1984；赵志新等，2003)计算二维柱坐标系下横向非均匀模型中地震波动 场。通过和直接解方法对比，检验了该方法的精度。提出实际发生地震计算时的矩张量投影和线性叠加方法， 提高实际地震计算的效率。通过模拟2003年6月20日发生在巴西的一个地震的波形，检验了该方法的计算精 度。本研究的正演结果说明错格富里叶伪谱法可以比较精确地计算区域地震的全波场，用于研究地球深部非均 匀结构对全波场的影响。 2．波动方程与数值模型 考虑弹性、各向同性地球介质，在柱坐标系(厂，秒，z)下，设所有空间变量沿z方向保持不变，N-维柱坐 标系p，护)下的P-SV波动方程可以写为： ∥，：三昙(％)+!等一旦盟+f， ∥一=71石0(F2GrO)+÷等+厶， 该点介质密度，厶=‘(r，0，f)(p=r，0)是体力，仃冈=仃M(，．，0，f)(p，q=r，0)是应力分量。 考虑的物理模型为通过某一大圆截取的三维球状地球模型的部分削面，模型中的横向非均匀设为沿方位角 方向对称。将该物理模型沿矢径(r)方向和横向(臼)进行空间离散化。在计算实际发生地震的波场时，使用 了哈佛大学CMT震源机制解给出的三维情况下的六个矩张量分量，通过坐标变换从六个分量中通过投影得到 所需的对应二维模型的三个分量。由于不同台站与震源所在的剖面互不相同，导致矩张量的投影分量也互不相 同。由于方程是线性的，所以采用叠加的方法来得到不同台站的记录，丽不必对每一个台站都进行计算，从而 大大减少实际地震的计算量。 3．计算结果与讨论 通过和直接解方法的对比，说明该方法的精度可靠。对2003年6月20日发生在巴西的一个地震的波形模 拟，说明对于实际发生地震，大多数震相的波形计算比较精确。通过PREM模型和已有层析成像模型模拟结果的 对比研究了横向非均匀结构对全波场的影响。本研究结果说明该方法应用于非均匀区域地球模型的地震波场计 算是可行的，对于进一步利用地震全波形信息研究地球深部精细结构提供了正演方法依据。 参考文献 赵志新等，地球物理学报，第46卷，第2期，2003。 Kosloffet a1．，砌肛Se括m．SocAm．，V01．74，875—891，1984． 本研究由国家自然科学基金项目和教育部留学回国启动基金资助． · · 525