ISSN 1004-4140
CN 11-3017/P
| Citation: |
LIU W G, TU W M, MOU Q S, et al. Seismic Numerical Simulation Based on Half-precision Floating-point Number Optimization and OpenMP[J]. CT Theory and Applications, 2024, 33(3): 289-297. DOI: 10.15953/j.ctta.2023.177. (in Chinese).
|
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.
| [1] |
BOHLEN T, WITTKAMP F. Higher order FDTD seismic modelling using the staggered adams-bashforth time integrator[C]//77 EAGE Conference and Exhibition 2015, Madrid, Spain. European Association of Geoscientists & Engineers, 2015: 1-5.
|
| [2] |
李斌, 温明明, 牟泽霖. 用于声波正演模拟的时空域高精度交错网格有限差分方法[J]. CT理论与应用研究, 2017, 26(3): 259−266. DOI: 10.15953/j.1004-4140.2017.26.03.01.
LI B, WEN M M, MU Z L. Finite difference method of high-precision interleaved grid in spatiotemporal domain for acoustic wave forward simulation[J]. CT Theory and Applications, 2017, 26(3): 259−266. DOI: 10.15953/j.1004-4140.2017.26.03.01. (in Chinese).
|
| [3] |
ZHAO L, GU L L, HAO L, et al. Study on characteristics of underground horizontal electric field caused by lightning return stroke based on FDTD method[J]. Meteorological and Environmental Research, 2019, 10(1): 39−42.
|
| [4] |
傅旦丹, 刘洋, 陈遵德, 等. 粘弹介质中声波的伪谱法模拟[J]. 江汉石油学院学报, 1993, 10(4): 32−39.
FU D D, LIU Y, CHEN Z D, et al. Pseudospectral simulation of sound waves in viscoelastic medium[J]. Journal of Jianghan Petroleum Institute, 1993, 10(4): 32−39. (in Chinese).
|
| [5] |
董良国, 马在田, 曹景忠. 一阶弹性波方程交错网格高阶差分解法稳定性研究[J]. 地球物理学报, 2000, 56(6): 856−864.
DONG L G, MA Z T, CAO J Z. Study on stability of high-order difference decomposition method of staggered grid of first-order elastic wave equation[J]. Chinese Journal of Geophysics, 2000, 56(6): 856−864. (in Chinese).
|
| [6] |
裴正林, 牟永光. 地震波传播数值模拟[J]. 地球物理学进展, 2004, 18(4): 933−941.
PEI Z L, MOU Y G. Numerical simulation of seismic wave propagation[J]. Progress in Geophysics, 2004, 18(4): 933−941. (in Chinese).
|
| [7] |
井涌泉, 高红伟, 王维红. 二维各向同性介质P波和S波分离方法研究[J]. 地球物理学进展, 2008, 22(5): 1412−1416.
JING Y Q, GAO H W, WANG W H. Study on separation method of P wave and S wave in two-dimensional isotropic medium[J]. Progress in Geophysics, 2008, 22(5): 1412−1416. (in Chinese).
|
| [8] |
龚明平, 张军华, 刘文革, 等. 各向异性介质地震波场数值模拟及特征分析[J]. CT理论与应用研究, 2017, 26(6): 661−668. DOI: 10.15953/j.1004-4140.2017.26.06.01.
GONG M P, ZHANG J H, LIU W G, et al. Numerical simulation and characteristic analysis of seismic wave field in anisotropic media[J]. CT Theory and Applications, 2017, 26(6): 661−668. DOI: 10.15953/j.1004-4140.2017.26.06.01. (in Chinese).
|
| [9] |
宋聂平, 韩俊刚, 李涛. 半精度浮点数处理器的设计与验证[J]. 信息技术, 2018, 41(1): 136−140.
SONG N P, HAN J G, LI T. Design and verification of half-precision floating-point processor[J]. Information Technology, 2018, 41(1): 136−140. (in Chinese).
|
| [10] |
曹克乾, 赵荣彩, 王磊. 基于申威1621的半精度超越函数设计[J]. 计算机应用, 2022, 42(S1): 221−225.
CAO K Q, ZHAO R C, WANG L. Design of half-precision transcendental function based on Shenwei 1621[J]. Journal of Computer Applications, 2022, 42(S1): 221−225. (in Chinese).
|
| [11] |
GABRIEL F O. Seismic modeling and inversion using half-precision floating-point numbers[J]. Geophysics, 2020, 85(3): 65−76. DOI: 10.1190/geo2018-0760.1.
|
| [12] |
LIU W G, WANG F, ZHOU H W. Parallel seismic modeling based on OpenMP+AVX and optimization strategy[J]. Journal of Earth Science, 2019, 30(4): 843−848. DOI: 10.1007/s12583-018-0831-3.
|
| [13] |
谭嘉言, 刘国峰, 高敬语. 基于GPU计算的地震波场高阶有限差分正演研究[J]. CT理论与应用研究, 2016, 25(1): 1−12. DOI: 10.15953/j.1004-4140.2016.25.01.01.
TAN J Y, LIU G F, GAO J Y. Research on high-order finite difference forward modeling of seismic wave field based on GPU calculation[J]. CT Theory and Applications, 2016, 25(1): 1−12. DOI: 10.15953/j.1004-4140.2016.25.01.01. (in Chinese).
|
| [14] |
肖调杰, 王赟, 刘云, 等. 大地电磁二维正演的MPI及OpenMP的并行计算对比[C]//中国地质学会2015学术年会论文摘要汇编(下册). 2015, 1(1): 219-220.
XIAO T J, WANG Y, LIU Y, et al. Parallel computing comparison of MPI and OpenMP with two-dimensional forward evolution of geodetic electromagnetism[C]//Abstracts of the 2015 Annual Conference of the Geological Society of China (Volume II). 2015, 1(1): 219-220. (in Chinese).
|
| [15] |
唐佳, 王凡, 刘福烈. 三维波动方程正演的三级并行加速[J]. 石油地球物理勘探, 2016, 51(5): 1049−1054, 1840.
TANG J, WANG F, LIU F L. Three-level parallel acceleration of three-dimensional wave equation forward[J]. Oil Geophysical Prospecting, 2016, 51(5): 1049−1054,1840. (in Chinese).
|
| [16] |
王凡. 时间域FWI的高性能计算方法研究[D]. 成都: 西南石油大学, 2018.
WANG F. Research on high-performance computing method of time domain FWI[D]. Chengdu: Southwest Petroleum University, 2018. (in Chinese).
|
| [17] |
FRANKE A, BORNER R U, SPITZER K. Adaptive unstructured grid finite element simulation of two-dimensional magneto telluric fields for arbitrary surface and seafloor topography[J]. Geophysical Journal International, 2007, 171(1): 71−86. DOI: 10.1111/j.1365-246X.2007.03481.x.
|
| [18] |
张伟, 陈晓非. 并行计算在计算理论地震图中的应用: 应用MPI/OpenMP在SMP-cluster机群采用三维有限差分计算地震图[C]//中国地球物理学会第二十届年会论文集, 西安: 西安地图出版社, 2004: 521.
ZHANG W, CHEN X F. Application of parallel computing in computational theoretical seismogram: Application of MPI/OpenMP to calculate seismogram with 3D finite difference in SMP-cluster fleet[C]//Proceedings of the 20th Annual Meeting of the Chinese Geophysical Society, Xi'an: Xi'an Map Press, 2004: 521. (in Chinese).
|
| [19] |
董良国. 二维声波传播波动方程数值模拟[M]. 上海: 同济大学出版社, 1999.
DONG L G. Numerical simulation of two-dimensional acoustic wave propagation wave equation[M]. Shanghai: Tongji University Press, 1999. (in Chinese).
|
| [20] |
裴正林. 三维各向同性介质弹性波方程交错网格高阶有限差分法模拟[J]. 石油物探, 2005, (4): 308−315, 15.
PEI Z L. Three-dimensional isotropic medium elastic wave equation staggered grid high-order finite difference method simulation[J]. Geophysical Prospecting for Petroleum, 2005, (4): 308−315, 15. (in Chinese).
|
| [1] | FENG Weiyu, SHEN Xuzhang, TAN Jiawei, WANG Wentian. Influence of Ambient Noise Source Distributions on Cross-correlation Derived Seismic Waveforms and Dispersion Curve Retrieval[J]. CT Theory and Applications. DOI: 10.15953/j.ctta.2025.181 |
| [2] | PAN Jishun, SUN Kaixuan, ZHANG Niu, LI Penghui. Non-SIRT Engineering Seismic Wave CT Technology[J]. CT Theory and Applications, 2020, 29(2): 119-127. DOI: 10.15953/j.1004-4140.2020.29.02.01 |
| [3] | SI Youqiang, GUO Runhua, LI Mengru. Application of Optimal Noise Reduction Smooth Model Based on EMD and EEMD in Seismic Waves[J]. CT Theory and Applications, 2020, 29(1): 11-21. DOI: 10.15953/j.1004-4140.2020.29.01.02 |
| [4] | GONG Ming-ping, ZHANG Jun-hua, LIU Wen-ge, LIU Yang, WU Xi-yuan. The Seismic Wave Numerical Simulation and Characteristics Analysis in Anisotropic Medium[J]. CT Theory and Applications, 2017, 26(6): 661-668. DOI: 10.15953/j.1004-4140.2017.26.06.01 |
| [5] | PAN Ji-shun, SONG Chao-yang, LENG Yuan-bao, LI Yan-zhuo, WANG Xin-jian, ZHAO Xiang, GAO Dong-pan, LI Di. The Application of Seismic Wave CT to Qualities Detection of Concrete Cut-off Walls[J]. CT Theory and Applications, 2016, 25(3): 311-317. DOI: 10.15953/j.1004-4140.2016.25.03.07 |
| [6] | TAN Jia-yan, LIU Guo-feng, GAO Jing-yu. High Order Finite Difference Forward Simulation of 2D and 3D Seismic Wave Field Based on GPU[J]. CT Theory and Applications, 2016, 25(1): 1-12. DOI: 10.15953/j.1004-4140.2016.25.01.01 |
| [7] | LIU Rui-he, ZHOU Jian-ke, ZHOU Xue-feng, YIN Xing-yao, CAO Dan-ping. Finite Element Method with Arbitrary Quadrilateral Meshes for Numerical Modeling of Seismic Wave Based on Compact Storage[J]. CT Theory and Applications, 2015, 24(5): 667-680. DOI: 10.15953/j.1004-4140.2015.24.05.04 |
| [8] | ZHENG Hai-liang, LI Xing-dong, WANG Zhe, WEI Cun-feng, CHANG Tong. A Practice on Parallel Reconstruction Algorithm of High Resolution Cone Beam Micro-CT Based on NVDIA GPU Graphic Card[J]. CT Theory and Applications, 2014, 23(5): 805-814. |
| [9] | FAN Hong-jun, FAN Ting-en, WANG Hui, GAO Yu-fei. Application of Seismic Waveform Classification Technique in the Study of Fluvial Reservoir[J]. CT Theory and Applications, 2014, 23(1): 71-80. |
| [10] | DAN Na-lin, ZHU De-bing, HUANG Xi-liu. The Test Method of Concrete Elements Integrity Using Seismic Image Profiling[J]. CT Theory and Applications, 2002, 11(2): 1-4. |
Supported by: Beijing Renhe Information Technology Co. Ltd 
DownLoad: