双侧向测井有限元正演算法的GPU并行实现.pdfVIP

中国地球科学联合学术年会 2014 ·1327 · 双侧向测井有限元正演算法的 GPU 并行实现 苏 俊 尚景涛 邹长春* 1.中国地质大学（北京）地球物理与信息技术学院 北京 100083 ； 2.地下信息探测技术与仪器教育部重点实验室和地质过程与矿产资源国家重点实验室 北京 100083 双侧向测井是油气勘探中的重要的电阻率测井方法之一，它利用屏蔽电极对发射电流进行聚焦，所 以具有较大的径向探测深度和较好的垂向分辨率。由于该仪器结构复杂，所以即使对于简单的地层模型 一般也没有解析解。为了更好的认识不同地层模型的响应特征，人们利用许多数值算法进行双侧向的正 演模拟，例如：有限差分法、有限元法、边界元法等等。 有限元法由于具有易处理复杂边界条件，容易提高计算精度，对求解区域剖分较为灵活等优点，在 计算地球物理领域得到了广泛的应用。李大潜 （1984）对有限元在电法测井中的运用进行了系统的阐述。 运用有限元算法会形成大型对称正定稀疏线性方程组，针对这种稀疏线性方程组，迭代求解算法，例如 共轭梯度算法（CG ）、预处理共轭梯度算法（PCG ）收敛速度快，精度高。 传统上一般利用串行算法实现迭代求解算法，当数据量很大时，计算过程将会很耗时。由于显卡 （GPU ）技术的进步，使得计算并行化得到了长足的发展。GPU 相对于 CPU 的优势在于每一个 GPU 都有成百上千个计算单元，它本身就具有高度的并行性和高密度数据的计算优势。CUDA （Compute Unified Device Architecture ）是在GPU 上编程的一种统一框架，利用 CUDA 编程时，是将程序的逻辑 控制部分交由 CPU 执行，而将计算部分交由 GPU 执行。由于 GPU 的特点，一般是将它运用于计算密 集型的任务中，例如：刘国峰等（2009 ）利用GPU 实现了地震资料处理的并行化。 GPU 编程模型中最重要的就是计算资源的调用，而显卡中每一个计算单元均被赋予了唯一的编号。 这样我们就可以将相同的重复计算过程分别在每一个计算单元上实现，而这些计算单元是同时进行计算 的，这样我们就可以实现GPU 并行化。 本文根据非常规油气中的页岩气储层的地质与物性特点，建立了页岩储层的地层模型；对双侧向有 限元正演算法进行了深入研究，发现其中求解线性方程组的共轭梯度算法比较耗时，因为它涉及大型矩 阵与矩阵、矩阵与向量相乘，并且也发现大型稀疏矩阵的存储是很耗费资源的。我们利用压缩稀疏行 （CSR ）存储技术对大型稀疏矩阵进行存储；根据矩阵相乘的特点，对其进行 GPU 并行优化，实现了 求解线性方程组的并行共轭梯度算法；与串行算法进行对比，发现随着迭代次数的增多，串行算法所需 的时间急剧增加，而并行算法所需时间增长的并不快，迭代十次以上，并行算法至少比串行算法节约 50%的时间。 资金资助：中央高校基本科研业务费专项资金（编号：2-9-2014-088） 参考文献 [1] 李大潜. 有限元素法在电法测井中的应用. 北京:石油工业出版社,1980. [2] 刘国峰,刘钦,李博等. 油气勘探地震资料处理 GPU/CPU 协同并行计算. 地球物理学进展, 2009, 24(5): 1671-1678.