3D声波全波形反演计算策略与研究.pdfVIP

·1160 · 中国地球科学联合学术年会 2014 3D 声波全波形反演计算策略研究 1、2* 1 2 1 1、2 尹 力 方伍宝 李振春 胡光辉 王振宇 1.中国石化石油物探技术研究院 江苏南京 210014 ；2. 中国石油大学（华东） 山东青岛 266580 油气工业界采集的地震资料大部分为三维资料，频率域全波形反演的提出使得全波形反演在 3D 实 际资料处理上具有了可行性。全波形反演对计算资源的要求是制约其发展的瓶颈之一。近年来随着计算 机技术的发展，3D 全波形反演实际资料的应用开始陆续出现。全波形反演包含两个过程即，正演过程 与反演过程。三维情况下仅仅使用时间域算法进行全波形的正反演，其反演过程需要巨大的存储空间， 耗费较长的 CPU 时间；而仅仅使用频率域算法，直接求解法的正演过程对内存存储要求极高，炮集的 LU 分解处理需要的内存通常达到几百 GB 甚至 TB 数量级，并且难以实现并行，耗费时间长。为优化 计算问题，本文提出频率域反演联合时间域正演的算法实现 3D 声波的全波形反演。 正演是反演的重要组成部分，正演方法的选择决定着反演的精度和效率。因此，能否正确的选择正 演方法决定了全波形反演的可行性，尤其是 3D 多炮点勘探情况下。 在 2D 情况下，频率域正演有不可 替代的优越性，3D 情况下则完全不同，3D 频率域正演采用 LU 直接分解法对内存的要求极高，甚至难 以实现。Operto 在其文章中指出在使用 SEG/EAGE 3D 推覆体模型进行波场模拟时，使用空间步长为 50m 的网格，正演频率为 10Hz 波场时，需要 320Gb 内存存储单精度的LU 分解结果，一般的计算机群 不能满足要求。为避免这一问题，我们在处理 3D 情况下的全波形反演将采取时间域正演方法。在正演 得到时间域地震记录后灵活运用离散傅里叶变换，将时间域波场转换到频率域来实现频率域的全波形反 演，这样我们不仅得到了频率域的波场值，还可以在时间域对数据进行所需要的预处理，并且一次正演 可以得到多个所需频率的波场值。反演也可以在时间域或频率域实现。时间域反演，数据量庞大，尤其 对 3D 大工区来说数据存储量和计算量都尤为巨大甚至不可实现。频率域反演的一个优势在于其利用一 个复数代替整个时间域各采样点，同时仅使用少数的几个离散频率即可得到高分辨率速度模型，从而大 大减小了存储空间和计算量。本文使用频率域反演联合时间域正演的混合算法，取二者的长处实现 3D 的全波形反演。其实现步骤如下：1)时间域正演获得时间域地震记录；2)运用离散傅里叶变换从时间域 地震记录中提取出几个所需频率的波场信息；3)频率域全波形反演。 实际资料验证表明，频率域反演联合时间域正演的混合算法降低了内存需求，节省了计算时间。此 外，采用伴随状态法求取梯度，避免了 Frechet 矩阵导数的计算，降低了反演的计算量以及程序开发复 杂度。并且随着伪保守形式的一阶速度-压力场声波方程的引入，使得梯度计算更加简单有效。 参考文献 [1] Operto S,Virieux J,Amestoy P,et al., 3D finite-difference frequency-domain modeling of visco-acoustic wave propagation using a massively parallel direct solver: A feasibility study.Geophysics, 2007,72(5):SM195-SM211. [2] 胡光辉，贾春梅，夏洪瑞等.3D 声波全波形反演的实现及应用[J].石油物探，2013,52(4) ：417-425