GPU-based simulation of seismic wave propagation with hybrid PSM/FDM method
-
摘要: 伪谱和高阶有限差分混合方法， 在垂直方向采用交错网格有限差分算子， 利用其并行程度高的特点， 在水平方向采用伪谱算子， 保留其高精度的优势， 是计算地震波场的有效方法. 图形处理器(graphic processing unit， 简写为GPU) 由于其高度并行性， 在计算此类问题中有显著的优势. 由英伟达(NVIDIA)公司推出的统一计算设备架构(compute unified device architecture， 简写为CUDA)平台极大地简化了GPU编程的难度. 为提高计算效率， 本文实现了基于CUDA 平台的混合方法二维地震波场模拟. 然后基于二维均匀介质模型将CPU与GPU版本的运行时间进行对比. 实际测试结果表明， 基于CUDA 的并行模拟方法在保证计算精度的同时显著地提高了计算速度， 为开展大规模非均匀地球介质地震波传播数值模拟提供了一种可选的方法.Abstract: Hybrid PSM/FDM method, approximating the horizontal spatial derivative with pseudospectral method (PSM) while using finite difference method (FDM) to approximate the vertical spatial derivative, is an efficient technique in seismic wave simulation because it makes usage of both the high level of parallelism in FDM and the great accuracy in PSM. With natural parallelism, graphic processing unit (GPU) is very suitable for calculating such kind of problem, and compute unified device architecture (CUDA) drastically simplifies GPU programming. In order to improve the computing efficiency, in this paper, we implemented the 2D hybrid method for seismic wave simulation based on CUDA. The GPU-based hybrid method brings about significant speedup according to numerical tests performed in this paper and proves to be a reliable way for simulating large scale seismic wave propagation.
-
-
![]()
图 5 兰州盆地模型地表水平分量ux和垂直分量uz的合成地震图
Figure 5. Horizontal and vertical component of synthetic waveforms in Lanzhou basin model
![]()
图 6 不同时刻时地震波传播的波场快照. 图中红色和绿色分别表示P和SV波
Figure 6. Wavefield snapshots at t=4，8，10，14，22，30 s，where red and green represent P and SV waves，respectively
表 1 测试环境
Table 1 Test environment
CPU 内存/GB GPU 显存/GB GCC版本 NVCC版本 Intel Core i5-4570 4 NVIDIA GeForce GTX 750 1 4.6.3 6 
下载: 导出CSV
表 2 计算耗时表
Table 2 Elapsed time
网格数 B/T GPU耗时/s CPU耗时/s 加速比 128×128 2 0.84 16.03 190.8 256×256 4 1.96 67.10 34.23 512×512 8 6.42 266.91 41.57 1024×1024 16 22.79 1186.99 52.08 2048×2048 32 87.98 5223.65 59.37 注： B为线程块数，T为每个线程块中的线程数.
下载: 导出CSV
-
郑亮, 张怀, 尹凤玲, 石耀霖. 2013. 基于GPU的大规模三维地震波传播有限元程序实现[C]//中国地球物理2013: 第七分会场论文集. 昆明, 2013年10月13日. Zheng L, Zhang H, Yin F L, Shi Y L. 2013. Finite element program implementation of large scale seismic wave propagation modeling based on GPU[C]//Chinese Geophysics 2013: Proceedings of the 7th Branch. Kunming, October 13, 2013.
Aoi S, Nishizawa N, Aoki T. 2012. Large scale simulation of seismic wave propagation using GPGPU[C]//15th World Conference on Earthquake Engineering. Lisbon, Portugal, September 24 to 28, 2012.
Cai C, Chen H Q, Deng Z, Chen D, Khan S U, Zeng K, Wu M X. 2012. GPGPU-aided 3D staggered-grid finite-difference seismic wave modeling[C/OL]//12th International Conference on Scalable Computing and Communications (ScalCom)/. [2014-05-12]. http://newtheme.hgpu.org/?p=8610.
NVIDIA. 2014. CUFFT Library User’s Guide[M/OL]. [2014-04-26]. http://docs.nvidia.com/cuda/cufft/.
Okamoto T, Takenaka H, Nakamura T, Aoki T. 2012. Large-scale simulation of seismic wave propagation of the 2011 Tohoku-Oki M9 earthquake[C]//Proceedings of the International Symposium on Engineering Lessons Learned from the 2011 Great East Japan Earthquake. Tokyo: Japan Association for Earthquake Engineering: 349-360.
-
期刊类型引用(17)
1. 吴丽慧，廖武林，黎金玲，包萨日娜. 黄石市地震台背景噪声特征分析. 地震工程学报. 2025(01): 152-159 . 
百度学术
2. 林丽萍，王宇航，吴朋，谌亮. 2022-06-01芦山M_S6.1地震流动地震台网部署及余震序列监测. 大地测量与地球动力学. 2024(02): 195-201 . 百度学术
3. 包文超，胡玮，申影，郭晔，甄齐，王西. 内蒙古基准站噪声特征与监测能力. 华北地震科学. 2024(02): 94-101 . 百度学术
4. 安全，韩晓明，包文超，翟浩，赵铁锁，赵星. 内蒙古地区强震动台站背景噪声与数据质量分析. 地震学报. 2024(03): 490-501 . 本站查看
5. 包文超，赵铁锁，郭雷，郭晔，王西. 额肯呼都格地震台地震观测波形干扰分析. 地震地磁观测与研究. 2024(03): 94-99 . 百度学术
6. 李文超，殷锴，刘炜，陈永新，王天琦，张子俊. 大同、朔州地震预警基准、基本站运行质量分析. 山西地震. 2024(03): 28-34 . 百度学术
7. 张小艳，安全，熊峰，苏亚梅，赵星. 内蒙古VP型垂直摆倾斜仪背景噪声特征分析. 大地测量与地球动力学. 2024(11): 1193-1198 . 百度学术
8. 唐浩，徐金银，何秋菊，马瑞，李青梅. 宁夏区域地震观测背景噪声特征. 大地测量与地球动力学. 2024(11): 1199-1206 . 百度学术
9. 白玛桑布，单增曲珠，次旦卓玛，洛桑罗布，普穷，格桑旦珍. 观测环境变化对日喀则地震台观测质量的影响. 地震地磁观测与研究. 2024(05): 67-73 . 百度学术
10. 方伟，韦进，黄雅，辜有启，胡敏章，刘高川，崔庆谷. 重力建站测试的影响因素研究——以云南燕子岩站为例. 地球物理学进展. 2024(06): 2126-2136 . 百度学术
11. 丁大业，董春丽，宫卓宏，陈永新，吕飞亚. 山西地震预警新建基准站数据质量评估. 山西地震. 2023(02): 15-22 . 百度学术
12. 文金龙，包文超，王西，范东海，甄齐，李男. 锡林浩特地震台台基背景噪声分析. 地震地磁观测与研究. 2023(01): 83-91 . 百度学术
13. 邓明文，许鑫，李金. 乌鲁木齐及周边地区新冠疫情前后背景噪声分析. 华南地震. 2023(03): 7-18 . 百度学术
14. 尹康达，李小军，张晓刚，丁成. 台基噪声功率谱估计中的Welch参数选择. 地震地磁观测与研究. 2023(04): 63-67 . 百度学术
15. 文金龙，贾彦杰，包文超，白少奇，郭晔. 锡林郭勒盟三个测震台站台基背景噪声对比分析. 河南科技. 2023(21): 99-103 . 百度学术
16. 曾均，赵寅甫，胡耀. 不同型号地震烈度仪振动测试对比分析. 四川地震. 2023(04): 13-18 . 百度学术
17. 刘艺璇，周元泽. 2013年芦山M_S7.0地震震源区背景噪声源特征分析. 中国地震. 2022(04): 721-736 . 百度学术
其他类型引用(0)
计量
- 文章访问数: 643
- HTML全文浏览量: 360
- PDF下载量: 10
- 被引次数: 17
