| 基于MPI并行算法的电阻率法多种装置数据的三维联合反演 |
| |
| 引用本文： | 马欢,郭越,吴萍萍,谭捍东.基于MPI并行算法的电阻率法多种装置数据的三维联合反演[J].地球物理学报,2018,61(12):5052-5065. |
| |
| 作者姓名： | 马欢 郭越 吴萍萍 谭捍东 |
| |
| 作者单位： | 1. 防灾科技学院地球科学学院, 河北三河 065201;2. 防灾科技学院电子科学与控制工程学院, 河北三河 065201;3. 中国地质大学(北京)地球物理与信息技术学院, 北京 100083 |
| |
| 基金项目： | 国家自然科学基金（41804119，41704055）和中央高校基本科研业务费资助计划项目（ZY20180205）联合资助. |
| |
| 摘 要： | 由于地表电阻率法受到浅地表局部异常体的干扰，反演精度受到影响，井中装置数据资料参与反演虽然可以减小浅地表局部异常体的干扰，但是由于钻井位置的局限性，数据量得不到保障，也会导致反演精度降低.为此，本文开发了一套结合地表、地-井、井-地和井-井装置数据的三维联合反演算法.首先，利用有限差分法实现正演模拟，采用非线性共轭梯度法（NLCG）恢复电阻率结构；其次，调用Message Passing Interface（MPI）函数库加速正演模拟和灵敏度矩阵运算，当开辟12个进程时，反演程序获得最大加速比4.51；最后，通过合成数据和实测数据算例证明该反演算法的有效性，也证实了多种装置组合数据体反演结果明显优于单一地表装置数据体反演结果.
|
| 关 键 词： | 多种装置数据体 电阻率法 三维 联合反演 MPI |
| 收稿时间： | 2018-01-08 |
3-D joint inversion of multi-array data set in the resistivity method based on MPI parallel algorithm |
| |
| MA Huan,GUO Yue,WU PingPing,TAN HanDong.3-D joint inversion of multi-array data set in the resistivity method based on MPI parallel algorithm[J].Chinese Journal of Geophysics,2018,61(12):5052-5065. |
| |
| Authors: | MA Huan GUO Yue WU PingPing TAN HanDong |
| |
| Institution: | 1. School of Earth Sciences, Institute Disaster of Prevention Science and Technology, Hebei Sanhe 065201, China;2. School of Electronic Science and Control Engineering, Institute Disaster of Prevention Science and Technology, Hebei Sanhe 065201, China;3. School of Geophysics and Information Technology, China University of Geosciences, Beijing 100083, China |
| |
| Abstract: | Because of interference from local anomaly bodies in the shallow subsurface, the inversion accuracy of ground-based resistivity method is limited. Although such interference can be reduced by virtue of borehole data in the inversion, due to the limitation of drilling position, the amount of data is not guaranteed, thus precision of inversion can be reduced. In this paper, a set of three-dimensional joint inversion algorithms that can simultaneously invert surface, surface-borehole, borehole-surface and borehole-borehole array data is developed. Firstly, forward modeling is conducted by the finite different method, and restoration of the resistivity structure is made by the Nonlinear Conjugate Gradient method (NLCG). Secondly, the Message Passing Interface (MPI) function library is used to accelerate forward modeling and sensitivity matrix calculation. The maximum speedup ratio of 4.51 is achieved by the inversion program when the 12 processes are opened up. Finally, the validity of the inversion algorithm is proved by numerical examples of synthetic data and measured data, and it is also proved that the inversion results of the combined data set of multiple arrays are obviously superior to those of a single surface array data set. |
| |
| Keywords: | Multi-array data set Resistivity method 3-D Joint inversion MPI |
|
| 点击此处可从《地球物理学报》浏览原始摘要信息 |
| 点击此处可从《地球物理学报》下载免费的PDF全文 |
|
