基于多线程并行计算的有限差分法弹性波数值模拟.pdfVIP

技术研发 TECHNOLOGY RD 基于多线程并行计算的 有限差分法弹性波数值模拟 周 洲 （国核电力规划设计研究院，北京 100095） 摘 要 地震波场的数值模拟技术是在已知地下介质结构和参数的情况下，利用理论计算的方法研究地震波在地下介 质中的传播规律，从而合成地震记录的一种技术方法。有限差分法作为最常用的正演模拟方法之一，其基本原理是将 波动方程中波场函数关于空间和时间的导数用相应差分来代替，也正是由于这种离散近似，不可避免地降低了数值模 拟结果分辨率。为了有效减小数值频散造成的影响，常用的方法即是提高模型的剖分精度，采用并行计算的方法能有 效提高求解效率。 关键词 并行计算；有限差分；弹性波数值模拟 中图分类号 ：P6 文献标识码 ：A 文章编号 ：1671-7597（2015）04-0076-01 近年来，计算机图像处理器 GPU 发展迅速，并参与到 CPU 的影响。 的大数据计算任务中，在特定应用如计算密集应用上有远超CPU 计算能力的优势。 本文以CUDA为应用开发平台，主要有以下几个研究重点： 1）在三维各向同性介质中，采用交错网格进行有限差分法 弹性波数值模拟，并采用 PML 完全匹配层对边界条件进行有效 处理。 2）研究基于 CUDA 的有限差分法弹性波数值模拟的具体实 现方法，并探讨多线程并行计算对数值模拟提速的实际效能以 进一步改进的潜能。 1 基于弹性波场数值理论进行正演模拟 图1 三维域中弹性波有限差分数值模拟效果图（未添加边界条件） 1.1 弹性理论 本文采用了完全匹配层（Perfectly Matched Layer， PML） 基于弹性理论基础及理想弹性体基本假设，论文从广义胡 吸收边界条件。其基本思想是在网格区域物理边界上设置一定 克定律出发，考虑连续介质中任意点应力与应变之间的关系， 厚度含有特殊属性的媒质，即吸收层，其波阻抗与相邻研究区 获得三维笛卡尔坐标系中，均匀各向同性线弹性介质的应力 - 内的媒质的波阻抗完全匹配，波从网格区域内部传到截断边界 应变关系写为分量形式： 处时，将无反射地穿过分界面。 σ (λ=+2µε) +λε +λε   xx xx yy zz σ λε =+(λ+2µε) +λε  yy xx yy zz  =+ + + σ λε λε (λ 2µε)  zz xx yy zz  σ σ 2µε 2µε  yz zy yz zy  σ σ 2µε 2µε  xz zx xz zx  σ σ 2µε 2µε  xy yx xy yx 1.2 波动方程的有限差分离散并添加边界条件 由于数值模拟计算的要求，需将弹性波支配方程在离散空 （a） （b） 间做有限差分离散。 图2 相同条件下 ，应用边界条件前后对比图 于是，对 f (x, y ,t ) 关于时间和空间的一阶偏导数取中心差 图 2 显示应用 PML