基于半精度浮点数优化与OpenMP的地震数值模拟  

Seismic Numerical Simulation Based on Half-precision Floating-point Number Optimization and OpenMP

作　　者：刘文革[1] 涂文茂 牟其松 陈康 周觅路 LIU Wenge;TU Wenmao;MOU Qisong;CHEN Kang;ZHOU Milu(School of Geosciences,Southwest Petroleum University,Chengdu 610500,China;SINOPEC Geophysical Corporation Nanfang Branch,Chengdu 610000,China;Exploration and Development Research Institute of CNPC Southwest Oil and Gas Field Branch,Chengdu,Sichuan 610400,China;CNOOC(China)Limited Hainan Branch,Haikou,Hainan 570100,China)

机构地区：[1]西南石油大学地球科学与技术学院,成都610500 [2]中石化石油工程地球物理有限公司南方分公司,成都610000 [3]中国石油西南油气田分公司勘探开发研究院,成都610400 [4]中海石油(中国)有限公司海南分公司,海口570100

出　　处：《CT理论与应用研究（中英文）》2024年第3期289-297,共9页Computerized Tomography Theory and Applications

基　　金：中国石油-西南石油大学创新联合体科技合作(地震自聚焦成像、多信息约束波形反演与地质解释一体化关键技术研究(2020CX010203));国家自然科学基金(羌塘盆地格架及沉积结构地球物理成像理论与识别技术研究(42241206))。

摘　　要：研究地震波场的传播特征是地震勘探的基础,利用有限差分法求解波动方程进行地震正演,其优势在于占用内存低、计算速度快、易于实现。但在三维情况下,计算量和数据量会急剧增加,传统串行地震数值模拟将无法满足计算效率的需求。因此,本文提出一种基于半精度浮点数优化与OpenMP的三维波动方程地震数值模拟方法,该方法首先利用半精度浮点数对地震常用的浮点型数据进行优化;其次利用应用程序接口OpenMP在多核CPU下通过以分割波场计算区域的方式实现并行计算;在保证计算结果满足精度需求的同时,能有效提高三维地震数值模拟的计算效率,并减少近一半的内存需求。通过数值试验证明该方法的有效性和实用性。Studying the propagation law of seismic wavefields is the basis of seismic exploration.The finite difference method is commonly used to solve the wave equation for seismic forward modeling,which has advantages of low memory usage,fast computation speed,and easy implementation.However,in the case of three dimensions(3D),the computational and data requirements increase dramatically and traditional,serial,seismic,numerical simulations can no longer meet efficiency demands.This study proposes a 3D wave equation,seismic,numerical simulation method based on half-precision floating-point numbers and OpenMP.This method optimizes commonly-used floating-point data in seismic calculations using half-precision floating-point numbers.It also utilizes the OpenMP application interface to achieve parallel computation on multi-core central processing units by dividing the wavefield calculation region.This strategy ensures the accuracy of the computation results while effectively improving the computation efficiency of 3D seismic numerical simulations and reducing memory requirements by almost half.The effectiveness and practicality of this method are demonstrated through numerical experiments.

关 键 词：有限差分法 半精度浮点数 OPENMP 并行计算 地震波 

分 类 号：O242[理学—计算数学] P315[理学—数学] P631[天文地球—地震学]