声波方程的最优近似解析离散化方法与其理论分析.pdfVIP

声波方程的最优近似解析离散化方法及其理论分析 杨顶辉 清华大学数学科学系，北京100084 数值仿真技术的研究一直是地震学和勘探地球物理学研究领域中的重要研究课题。大量研究表明数值模拟地 震波场已成为研究复杂介质中地震传播的重要手段。随着地震学和地震勘探理论与实用研究的深入，我们所丽临 的实际问题越来越复杂，现有的数值模拟技术有时很难解决我们所面对的复杂问题，为了更好地认识地球内部构 造和石油勘探的实际需要，必须研究求解声波和弹性波方程的更有效、更精确的数值模拟新方法，特别是需要发 展复杂介质中波传播模拟的快速、弱数值频散、低存储量的数值仿真新技术、新方法。 纵观国内外对数值仿真技术及其波场模拟的研究现状和发展趋势，不难发现，目前比较流行的数值模拟方法 主要包括；有限元法与边界元法、射线追踪法、反射率法、有限差分法、谱(或伪谱、或谱元)法，等等。每一 eta1．，2003)。 种方法都有其优点和不足之处，如：多数数值方法在粗网格条件下存在严重的数值频散(详见Yang 因此发展有别于上述方法的高效、高保真、弱网格频散、低存储费用，且适用于复杂介质的数值模拟新方法和新 技术具有特别重要的理论和实际意义。 近似解析离散化方法(NA蹦)是近年发展的一种新方法，此方法最初用于求解双曲型和抛物型方程(Konddoh et eta1．．1994)，最近被发展用于求解二维声波和弹性波方程(Yanga1．，2003)。此方法的基本思想是：在时间 上采用截断泰勒展式来解析地表示波位移及其一阶偏导数，在空间上利用截断的泰勒公式来构造高阶插值函数以 逼近高阶空间偏导数，这种处理使得因离散表示地震位移而导致的高阶地震信息得以补偿，有效地压制了因数值 离散波方程而产生的网格频散，从而与上述被普遍应用的方法相比，NADM具有粗网格条件下弱数值频散、且能有 et 效压制源噪声、并有较高的数值精度和较好的数值稳定性等优点(Yanga1．，2003)。然而，这种方法尚有某 些不足，如；需要计算速度，从而增加了存储空间，计算效率也有待提高：同时，应用NADM求解波动方程时遇到 的人为边界问题也尚未解决。为此，更最近的研究工作优化了这种方法，获得了二维情况下优化形式的NA阴，并 et 解决了人为边界问题，极大地提高了这种方法的计算精度、分辨率以及计算效率，并节省存储量约53％(Yang a1．。2004)，适用于大尺度波场数值模拟，具有巨大的应用前景。 然而，迄今为止，对复杂介质中的优化NADM尚无任何报道，其相关的理论问题(如：稳定性、数值频散理论 分析等)也未解决。为此本文给出了最优近似解析离散化方法应用于声波方程时的数值频散理论分析、数值方法 稳定性条件，以及应用于非均匀介质的波场模拟结果，并与目前被普遍使用的数值方法进行了比较，为优化NADM 的进一步研究和应用迈出了重要一步。 参考文献 。 Konddoh，Y．，Hosaka，Y．，andIshii，K．，1994，Kernel discretization 11l optimumnearly—analytical algorithm to and appliedparabolichyperbolicequations，ComputersMath．Applic．，27(3)。59-90． methodforacousticandelastic