起伏地形条件下二维声波频率域全波形反演.pdf

2015 年 9 月 APPLIED GEOPHYSICS （应用地球物理） 第 12 卷，第 3 期 起伏地形条件下二维声波频率域全波形反演//Two-dimensional frequency-domain acoustic full-waveform inversion with rugged topography, APPLIED GEOPHYSICS, 2015, 12(3), P. 378-388. DOI: 10.1007/s11770-015-0510-4 起伏地形条件下二维声波频率域全波形反演 1 ♦1 1 1 1 2 张钱江 ，戴世坤 ，陈龙伟 ，李昆 ，赵东东 ，黄兴兴 （1. 中南大学 地球科学与信息物理学院，长沙 410083 ；2. 中国石油大学(北京) 地球物理与信息工程学院，北京 102249 ） 摘要：本文以有限单元法为基础，对起伏地形条件下二维声波频率域全波形反演 进行了研究。在正演算法中，针对截断边界问题，并考虑多频率联合反演中计算 区域采用同一套剖分网格的需求，提出了一种适用于起伏地形的衰减边界条件算 法。该算法的核心思想是在控制方程波数项中引入衰减因子，通过一定方式调节 衰减因子使得声波在衰减层中充分衰减，达到压制截断边界影响的目的。根据指 数衰减规律，文中推导出了一种新的衰减因子计算公式，并给出了不同频率条件 下衰减层厚度计算公式；在反演算法中，采用共轭梯度法求解高斯牛顿反演迭代 方程组，避免直接求解雅克比矩阵和 Hessian 矩阵带来的巨额计算量，并采用相同 的反演模型，对比分析了不同初始模型和频率组合对全波形反演结果的影响。起 伏地形模型数值模拟和全波形反演数值试验表明，本文提出的指数衰减边界条件 算法和基于该算法的全波形反演算法具有很好的应用效果。 关键词：全波形反演，起伏地形，衰减边界条件，有限单元法 引言 2011; Hu et al., 2013) 。 波形反演从计算方式上分为三类:时间-空间 域波形反演(Lailly, 1983; Tarantola, 1984)，频率 随着资源勘探的深入，勘探区域的地形条件 域波形反演(Prattet al., 1990)和拉普拉斯域波形 和地质构造越来越复杂，地震勘探中常规的速度 反演(Shin et al., 2008, 2009) 。相比其它两种波形 分析不能够满足高精度成像的要求。而诞生于 反演，频率域波形反演具有如下优势:克服了时 20 世纪 80 年代的全波形反演能够提供精确的速 间方向上的局部极小值问题；只需要有限个频率 度 模 型 (Lailly, 1983; Tarantola, 1984, 1986; 就能反演得到精度很高的速度模型；由于边界条 Gauthier,1986)，因此近些年全波形反演得到了大 件与炮点位置无关，多炮条件下线性方程组可以 量的研究，特别是挪威北海油田(Valhall 地区) 的 回代求解，大大节省了计算量和计算时