二阶弹性波动方程高精度交错网格波场分离数值模拟-物探与化探.pdf

　第３３卷第６期 物　探　与　化　探 Ｖｏｌ．３３，Ｎｏ．６　 　２００９年１２月 ＧＥＯＰＨＹＳＩＣＡＬ＆ＧＥＯＣＨＥＭＩＣＡＬＥＸＰＬＯＲＡＴＩＯＮ Ｄｅｃ．，２００９　 二阶弹性波动方程高精度交错网格波场分离数值模拟 １ １ １ ２ 陈可洋 ，杨微 ，刘洪林 ，吴清岭 （１．大庆石油学院地球科学学院，黑龙江 大庆　１６３３１８；２．大庆油田有限责任公司勘探开发研究 院，黑龙江 大庆　１６３７１２） 摘 要：给出了一种等价的二阶弹性波动方程，以解决弹性波场中完全弹性波动方程不能完全分离耦合的纵、横波 波场问题。应用高阶交错网格有限差分法求解该波动方程，并使用通量校正技术（ＦＣＴ）进一步压制频散，采用均 匀介质模型和层状介质模型进行波场分离数值试验，精确得到了混合波场、完全分离的纯纵波及纯横波波场。数 值结果分析表明，本文方法在均匀介质情况下准确可靠，在分离后的纯纵、横波波场中可观察到较为丰富的能量转 换信息，这对认识复杂弹性波的传播规律及弹性波理论具有重要意义。 关键词：地震波场分离；高阶交错网格；等价二阶弹性波动方程；数值模拟；通量校正技术 中图分类号：Ｐ６３１．４　　　文献标识码：Ａ　　　文章编号：１０００－８９１８（２００９）０６－０７００－０４ 　　地震波场数值模拟技术一直是勘探地球物理领 高阶交错网格有限差分法求解，仅需求解６个方程， 域内最为活跃的研究内容之一。利用精确的波动方 模型纵横向网格节点数不受方法限制，可精确得到 程数值解模拟复杂地下波场，为研究地震波传播机 混合波场及完全分离的纯Ｐ波和纯Ｓ波波场，且完 理、地震资料的特殊处理方法以及复杂地层的解释 全保留纵横波能量转换信息，这对认识弹性波的传 等许多方面提供更为科学的数学物理依据［１－３］。 播规律，指导实际资料处理与解释具有重要意义。 在实际勘探中仍然存在许多复杂问题值得深入 １　算法原理 探讨，这对研究波动理论，指导地震资料的特殊处理 与复杂构造解释具有重要意义。例如，多波多分量 １．１　等价的弹性波动方程 地震记录的每一个分量均包含不同的波型，简单地 在三维各向同性完全弹性介质中的弹性波动方 把垂直分量看成Ｐ波，水平分量看成Ｓ波是很不合 程由如下的运动方程描述［９］ 理的。 ２Ｓ  ２ ２ ＝ｖ （ ·Ｓ）－ｖ ×（ ×Ｓ） （１）     ２ ｐ ｓ ｔ 常规地震资料处理总是希望处理单一分量的标  Ｔ 量波场，因此从混合波场中解析出纯Ｐ波和纯Ｓ波 其中，Ｓ＝（ｕ，ｖ，ｗ）为位移场，ｖ为纵波速度，ｖ为 ｐ ｓ 场是进行速度分析、偏移成像等地震资料数据处理 横波速度。 的前提。