Born近似快速三维反演井地电法数据

本篇研究了井中电偶极激发地面接收的井地电法的快速反演成像问题.我们采用了Born近似方法和重加权正则化共轭梯度法(RRCG)算法.数值计算的结果表明Born近似是一种有效的井地电法三维快速反演方法,同时也说明井地电法监测油水前驱和储层边界预测的观测数据可以用该方法进行快速三维反演成像.     

可控源音频大地电磁三维共轭梯度反演研究

可控源音频大地电磁法在资源勘探等领域中发挥着重要的作用.我们把有限差分数值模拟方法用于可控源音频大地电磁三维正演,结合正则化反演方案和共轭梯度反演的思路,将反演中的雅可比矩阵计算问题转为求解两次"拟正演"问题,得到模型参数的更新步长,形成反演迭代,实现了可控源音频大地电磁三维共轭梯度反演算法.该反演算法可用于对有限长度电偶源激发下采集到的可控源音频大地电磁全区(近区、过渡区和远区)视电阻率和相位资料进行三维反演定量解释,获得地下三维模型的电阻率结构.理论模型合成数据的反演算例验证了所实现的可控源音频大地电磁三维共轭梯度反演算法的有效性和稳定性.     

大地电磁非线性共轭梯度拟三维反演

提出了非线性共轭梯度法大地电磁拟三维反演.该方法选取共轭梯度反演算法为拟三维反演的核心.在计算灵敏度（Jacobian）矩阵时,吸取近似灵敏度矩阵思想,采用一维灵敏度矩阵来代替三维灵敏度矩阵,并对非测点的灵敏度元素提出一种近似方法.在第一次反演之后,采用拟牛顿法更新灵敏度矩阵.拟三维反演法在很大程度上节省了计算时间,并且理论模型和实际资料的反演试算结果表明大地电磁拟三维反演法具有一定的实用价值.     

基于数据空间和稀疏约束的三维重力和重力梯度数据联合反演

随着重力和重力梯度测量技术的日趋成熟,基于重力和重力梯度数据的反演技术得到了广泛关注.针对反演多解性严重、计算效率低和内存消耗大等难点问题,本文开展了三维重力和重力梯度数据的联合反演研究,该方法结合重力和重力梯度两种数据,将L0范数正则化项加入到目标函数中,并在数据空间下采用改进的共轭梯度算法求解反演最优化问题.同时,本文摒弃了依赖先验信息的深度加权函数,引入了自适应模型积分灵敏度矩阵,用来克服因重力和重力梯度数据核函数随深度增加而衰减引起的趋肤效应问题.为了提高反演计算效率,本文又推导出基于规则网格化的重力和重力梯度快速正演计算方法.模拟试算表明,改进的共轭梯度法可以降低反演的迭代次数,提高反演的收敛速度;自适应模型积分灵敏度矩阵,可以有效解决趋肤效应,提高反演纵向分辨能力;数据空间和改进的共轭梯度算法结合,可以更好地降低反演求解方程的维度,避免存储灵敏度矩阵,有效地降低反演计算时间和内存消耗量.野外实例表明,该算法可以在普通计算机下快速地获得地下密度分布模型,表现出较强的稳定性和适用性.

     

基于非结构网格剖分的海洋可控源电磁三维正则化反演研究

为更好地处理与解释复杂海底地形条件下测得的海洋可控源电磁数据,本文提出了一种基于非结构网格剖分的频率域海洋可控源电磁数据三维正则化反演方法.该方法首先对海洋地电模型以非结构四面体单元进行离散,然后基于矢量有限元方法获得海洋可控源电磁响应和灵敏度信息,最后采用共轭梯度法求解高斯-牛顿反演方程计算模型修正量.为提高反演的稳定性,通过在反演过程中采用对数转换方法实现反演模型参数的上下限约束.本文分别测试了单测线水平海底地形反演算例和面积性测量的起伏海底地形反演算例.反演结果表明,本文提出的频率域海洋可控源电磁三维反演能够准确地恢复高阻储油层的位置和电阻率信息,且计算效率较高,可用于实测海洋电磁资料的处理与解释.

     

2.5维海洋可控源电磁反演算法及影响参数研究

海洋可控源电磁法作为一种较新的地球物理勘探方法,在国外已被成功应用于海洋油气资源与天然气水合物的探测中,而我国在这方面的研究开始较晚.本文编写了基于非线性共轭梯度算法的海洋可控源电磁2.5维反演程序,结合理论数据对反演代码的正确性进行验证并讨论了反演参数对反演运算速度及效果的影响.研究结果表明,1)本文编写的2.5维反演程序正确可靠,计算结果与理论模型一致,2)在非线性共轭梯度反演计算中,正则化因子、预条件矩阵、线性搜索及共轭梯度更新因子等参数对实际的反演速度及精度都存在一定影响,可根据不同的勘探需求,调整反演参数以达到较好的反演结果.     

基于矢量有限元的带地形大地电磁三维反演研究

研究了基于矢量有限元方法的大地电磁带地形三维反演算法并开发了三维反演计算程序代码.在大地电磁场正演数值模拟方面,采用并行直接稀疏求解器PARDISO且无需进行散度校正的快速正演方案,对典型地形模型,在中等规模计算条件下,与双共轭梯度法（BICG）计算结果比较,发现PARDISO比BICG快10倍以上;通过理论模型试算,并与前人的有限元法计算结果对比,验证了带地形三维正演计算程序的正确性.在反演方面,本研究基于共轭梯度方法编写了大地电磁带地形三维反演代码,为了避免直接求取雅可比矩阵,将反演中的雅可比矩阵计算问题转为求解两次"拟正演"问题,进而将PARDISO的快速正演方案应用于"拟正演"问题的求解,以提高反演计算效率.利用开发的反演算法对多个带地形地电模型的合成数据进行了三维反演,反演结果能很好地重现理论模型的电性结构,验证了本文开发的三维反演算法的正确性和可靠性.最后,利用该算法反演了某矿区大地电磁实测数据,反演得到的三维电性结构清晰地反映了研究区的地电特征,将反演结果与该区已有地质资料结合进行解释,应用效果明显,进一步验证了本文算法的有效性.     

井地电磁法三维正演模拟和Born近似反演

本文提出一种数值模拟井地电磁法的方法。用体积分方程法对层状介质中的垂直长导线源三维电磁响应做了三维模拟。模拟的结果与解析解对比误差很小,说明算法是正确的。开发了井地电磁法Born近似反演程序,理论模型合成的数据反演结果非常好。     

起伏地形下大地电磁L-BFGS三维反演方法

为推进大地电磁三维反演的实用化,本文实现了基于L-BFGS算法的带地形大地电磁三维反演.首先推导了大地电磁法三维反演的Tikhonov正则化目标函数以及Hessian矩阵逆矩阵近似表达式和计算方法,然后设计了一种既能保证空气电阻率固定不变又能保证模型平滑约束的协方差矩阵统一表达式,解决带地形反演问题.在反演算法中采用正则化因子冷却法以及基于Wolf条件的步长搜索策略,提升了反演的稳定性.利用开发的算法对多个带地形地电模型（山峰地形下的单个异常模型、峰-谷地形下的棋盘模型）的合成数据进行了三维反演,并与已有大地电磁三维反演程序（ModEM）进行对比,验证了本文开发的三维反演算法的正确性和可靠性.最后,利用该算法反演了华南某山区大地电磁实测数据,得到该区三维电性结构,揭示了研究区以高阻介质为基底,中间以低阻不整合面和相对低阻介质连续分布,浅部覆盖高阻介质的电性结构特征,进一步验证了本文算法的实用性.

     

双参数混合正则化方法及在电导率反演成像中的应用

通过引入带有二阶正则算子的正则化项,建立了一种双参数混合正则化方法.为确定最佳正则化参数,这里主要应用L-曲线法、偏差原理和广义交叉校验准则的优化组合来确定.首先对理论模型进行了数值模拟,通过与截断奇异值分解法、共轭梯度法及标准Tikhonov正则化法的结果比较,表明该方法不仅精度高,而且对于数据的随机扰动具有稳定性.然后将此方法应用于对电导率反演成像的计算中,得到的电导率反演成像精确细致可靠,也符合实际情况.这表明该方法求解地球物理反问题时,在增强反演的稳定性及减少多解性的同时,还能进一步提高反演计算速度.     

3D Modeling and Born approximation inversion for the borehole surface electromagnetic method

We present a 3D approach to numerical modeling of the borehole-surface electromagnetic （BSEM） method. The 3D electromagnetic response created by a vertical line current source in a layered medium is modeled using the 3D integral equation method. The modeling results are consistent with analytical solutions. 3D Born approximation inversion of BSEM data is also conducted for reservoir delineation. The inversion method is verified by a synthetic reservoir model.     

基于MPI并行算法的电阻率法多种装置数据的三维联合反演

由于地表电阻率法受到浅地表局部异常体的干扰,反演精度受到影响,井中装置数据资料参与反演虽然可以减小浅地表局部异常体的干扰,但是由于钻井位置的局限性,数据量得不到保障,也会导致反演精度降低.为此,本文开发了一套结合地表、地-井、井-地和井-井装置数据的三维联合反演算法.首先,利用有限差分法实现正演模拟,采用非线性共轭梯度法（NLCG）恢复电阻率结构;其次,调用Message Passing Interface（MPI）函数库加速正演模拟和灵敏度矩阵运算,当开辟12个进程时,反演程序获得最大加速比4.51;最后,通过合成数据和实测数据算例证明该反演算法的有效性,也证实了多种装置组合数据体反演结果明显优于单一地表装置数据体反演结果.

     

基于GPU并行的重力、重力梯度三维正演快速计算及反演策略

利用NVIDIA CUDA编程平台,实现了基于GPU并行的重力、重力梯度三维快速正演计算方法.采用当前在重力数据约束反演或联合反演中流行的物性模型(密度大小不同、规则排列的长方体单元)作为地下剖分单元,对任意三维复杂模型体均可用很多物性模型进行组合近似,利用解析方法计算出所有物性模型在计算点的异常值并累加求和,得到整个模型体在某一计算点引起的重力(或重力梯度)值.针对精细的复杂模型体产生的问题,采用GPU并行计算技术,主要包括线程有效索引与优化的并行归约技术进行高效计算.在显卡型号为NVIDIA Quadro 2000相对于单线程CPU程序,重力和重力梯度Uxx、Uxy正演计算可以分别达到60与50倍的加速.本文还讨论了GPU并行计算在两种反演方法中的策略,为快速三维反演技术提供了借鉴.     

可控源音频大地电磁数据的反演方法

从反演方程、构造目标函数和求解三方面对用于可控源音频大地电磁法(CSAMT)的实用反演方法中的四种进行了描述．水平层状地层CSAMT法资料的直接反演法首次尝试了一维空间的全资料CSAMT反演，效果较好，但该方法尚难应用于2D、3D复杂介质中；奥克姆反演方法既考虑了横向的光滑函数，又考虑了纵向的光滑函数，得到比较光滑的横向、纵向变化的背景电性结果，但有可能把一些小构造光滑掉．快速松驰反演算法和共轭梯度算法由于计算速度快，占内存少而被用于三维反演中，二者相比，快速松驰算法在求解雅可比矩阵时只做一次正演计算，在更新模型时解小型方程组，所以在速度上更胜一筹．在后三种算法中，由于复杂电性结构无解析解，正演计算都采用数值计算．数值计算的可靠性、速度影响着反演算法的有效性，这方面的研究也将是2D、3D复杂电性结构反演的研究方向之一.     

An application of three-dimensional electrical resistivity imaging for the detection of an underground waste-water system

A 3D electrical resistivity imaging survey is presented in this paper. The objective was to investigate an underground wastewater system at the University of Malaya, Malaysia. Apparent resistivity data were collected along ten parallel lines using a Wenner-Schlumberger configuration; electrode cables were oriented in the x-direction with 3 m spacing. Roll-along measurements using a line spacing of 3 m were carried out covering a grid of 20 × 10 electrodes. All data sets were merged into a single data file in order to perform a 3D inversion. Two different 3D least squares algorithms, based on the robust inversion method and the smoothness-constrained technique, were used for the inversion of the apparent resistivity data. Both the horizontal and vertical extents of the anomalous zones found by inversion are displayed. The results indicate the superiority of the robust inversion method over the smoothness-constrained technique at this site. The results are in sufficient accordance with previously known information about the investigation area. The results show that 3D electrical resistivity imaging surveys, in combination with an appropriate 3D inversion method, can be highly useful for engineering and archaeological investigations as well as for environmental applications.     

井-地与井间电位技术联合数值模拟研究剩余油分布(英文)

电法测井有很高的分辨率,但是它的探测半径仅限于井孔周围;井-地电位技术虽然可以探测到足够大的范围,但是它的分辨率却受到很大的限制,特别是对于油水分布或者结构复杂的储层。本文试图通过井-地电位技术和井间电位技术的联合来研究地下储层油水分布范围。具体方法是,采用井-地电位技术研究油水分布在在横向上的展布情况,利用井间电位技术研究油水储层在垂向上的分布,然后采用井间电位结果标定井地电位结果,两者结合提高纵向分辨率,从而确定剩余油的三维空间分布。研究中通过研究注水初期数值模拟结果与水淹期数值模拟结果之差,求取剩余油分布范围。有限差分方法数值模拟表明:井-地电位技术与井间电位技术联合方法可以有效地确定剩余油分布。     

基于不等式约束的最小二乘法三维电阻率反演及其算法优化

基于光滑约束的最小二乘法是三维电阻率反演的主要方法,但该方法在某些情况下存在着多解性较强的问题,且普遍耗时较长,严重制约了三维反演方法的推广与发展.为改善上述问题,将表征模型参数变化范围的不等式约束作为先验信息引入最小二乘线性反演方法中,有效地改善了反演结果的精度,降低了反演的多解性问题.为了解决耗时较长的问题,基于预条件共轭梯度(PCG)算法和Cholesky分解法的特点提出了一套优化三维电阻率反演计算效率的计算方案.在该方案中,Cholesky分解法被用来求解敏感度矩阵计算中的多个点源场的正演问题,Cholesky分解法只需对总体系数矩阵进行一次分解,然后对不同的右端向量进行回代即可.将预条件共轭梯度法引入到三维电阻率反演方程的求解中,将雅可比迭代中的对角阵作为预处理矩阵,其具有求逆方便、无需内存空间的特点,有效地加快了收敛速度.对合成数据以及实测数据的反演算例表明,借助不等式约束和反演效率优化方案,最小二乘反演方法可得到较为精确的反演结果,有效地提高了反演计算效率,具有良好的推广前景.