基于正则参数后验策略的PCS迭代正则化反演

采用改进的正则参数后验策略的迭代正则化方法对单峰分布和双峰分布颗粒系的电场自相关函数进行反演,并且与正则化方法反演结果进行了比较.实验结果表明,在噪声小于0.005的情况下,正则化方法和迭代正则化方法都能较好地反演出单峰和双峰分布的颗粒,当噪声水平为0.05时,正则化反演算法已无法得出粒径分布,迭代正则化方法仍能比较准确地反演出粒径分布.     

二维直流电阻率与音频大地电磁自适应渐进联合反演

开展了非结构化三角网格条件下的二维直流电阻率与音频大地电磁自适应渐进正则化联合反演研究。通过研发以模型灵敏度信息为依据的反演网格自适应优化技术,构建由粗网格到细网格逐步反演的自适应渐进反演策略,减少了反问题对稳定因子的依赖,进而降低了正则化因子搜索的计算量;通过最小二乘算法求解非结构化三角网格的模型粗糙度,构建了非结构化网格条件下的最小结构稳定因子;采用高斯-牛顿法优化求解正则化反演目标函数,通过双共轭梯度稳定算法求解高斯-牛顿方程,确保反演稳定性,同时有效减少了反问题对内存的需求。理论模型与实测数据的联合反演试算表明,直流电阻率与音频大地电磁数据联合反演能够有效减少反演多解性,提高反演效果,实现音频大地电磁静态效应异常体的直接反演。     

二维MT电阻率与磁导率同时反演研究

开展了同时考虑介质电阻率与磁导率变化的大地电磁二维正则化反演问题研究。不考虑位移电流影响的大地电磁勘探,电阻率与磁导率是影响观测电磁场的主要物理参数;考虑岩石的电阻率变化范围相比磁导率大的多,所以大地电磁反演解释常常只考虑电阻率;然而,二维正演模拟表明,当地下介质含有较多铁磁性矿物时,磁导率的局部异常体对大地电磁观测的影响不能忽略。因此,开展了同时反演电阻率与磁导率的正则化反演,建立了同时反演两个参数的正则化反演目标函数,并通过高斯-牛顿法对正则化目标函数进行最优化求解;采用模型空间到反演空间的函数映射技术,将电阻率与磁导率限定在具有实际物理意义的变化区间,并通过双共轭梯度稳定算法求解高斯-牛顿系统,以确保反演稳定性。对合成数据与实测数据进行了反演,表明从大地电磁数据可以恢复地下介质磁导率,开展磁导率与电阻率的同时反演提高了反演效果,丰富了数据解释依据。     

椭球体上双频激电法的正演与反演算法

基于柯尔-柯尔频散理论,采用中间梯度上椭球体的复电阻率计算相应的幅频率,讨论椭球体的充电率、中心埋深和长轴倾角对幅频率异常和曲线形态的影响;研究幅频率资料的反演算法。研究结果表明:由于中间梯度剖面数据量非常有限,选择最小二乘迭代拟合算法是合适的;偏导数矩阵采用差分的方法求得,对模型参数进行无量纲处理,提高了反演效率;该反演算法迭代速度快,能够稳定收敛,反演效果较好。     

一种适合磁异常反演的正则化因子选取方法

为更精确地确定地下磁性物质的分布,结合工程勘察中的磁法探勘和位场反演特点,借鉴Zhdanov等提出的自适应正则化因子选取方式,提出一种适合磁异常反演的完全自适应正则化方法。在进行一次常规自适应算法的基础上,得到估计的正则化因子初始值,根据数据拟合的下降效率确定合适的正则化因子的计算值,并采用较大的衰减因子,通过自动衰减正则化因子在考虑计算效率的条件下得到较好的数据拟合和稳定的反演结果。模型实验表明,此方法的数据拟合效果比L曲线法更好,同时在解稳定的前提下减少了迭代次数。     

基于自适应粒子群优化和MPI的大地电磁并行反演算法

针对大地电磁粒子群反演算法存在的计算时间过长的问题,基于自适应粒子群优化和消息传递接口提出一种新的大地电磁并行反演算法.在曙光4000L大型机平台上,利用该并行反演算法进行一维大地电磁层状介质反演实验,实验结果表明,新的并行反演算法能有效解决大地电磁粒子群反演计算时间过长的问题.     

大地电磁主轴各向异性二维理论模型反演分析

地下介质普遍存在各向异性,为反应地下真实地质情况,需建立各向异性模型进行研究。本文主要进行主轴各向异性大地电磁二维反演研究,采用有限单元法作为正演模拟方法,将正演响应结果与前人计算结果进行对比分析,验证了算法的正确性,并对不同方向上电阻率响应灵敏度进行分析。采用非线性共轭梯度法进行反演研究,该方法不需要直接计算雅克比矩阵,相对于其他反演方法节省了计算时间,计算效率高;目标函数的建立对于不同方向上的电阻率采用不同的正则化参数来进行约束。通过复杂的各向异性体模型进行反演,结果显示对于不同方向上电阻率的恢复以及异常体位置的圈定都取得了较好的效果,说明了反演算法的有效性。     

大地电磁反演中改进的自适应正则化因子选取

通过建立大地电磁(MT)层状地电模型,利用共轭梯度法求解反问题,在给定不同初始模型的条件下对多种正则化因子选取方法进行了计算比较,分析了各种方法的特点和使用条件.结果表明,自适应正则化算法的效果与传统的定值方法如L曲线法相近,但反演过程远比传统方法便捷.为了解决反演依赖于初始模型的局限并增强解的稳定性,基于多种自适应正则化方案的对比分析,提出了改进的自适应正则化方案,选取数据拟合泛函与模型稳定泛函较大的比值为正则化因子的初始值,并提出相应的调整方案自动控制正则化因子衰减.模型试验表明,该方法对初始模型的依赖性低于其他几种自适应的算法,反演结果的稳定性较强,可以进一步提高正则化反演的效率.     

大地电磁反演中改进的自适应正则化因子选取研究

通过建立大地电磁(MT)层状地电模型,利用共轭梯度法求解反问题,在给定不同初始模型的条件下对多种正则化因子选取方法进行了计算比较,分析了各种方法的特点和使用条件.结果表明,自适应正则化算法的效果与传统的定值方法如L曲线法相近,但反演过程远比传统方法便捷.为了解决反演依赖于初始模型的局限并增强解的稳定性,基于多种自适应正则化方案的对比分析,提出了改进的自适应正则化方案,选取数据拟合泛函与模型稳定泛函较大的比值为正则化因子的初始值,并提出相应的调整方案自动控制正则化因子衰减.模型试验表明,该方法对初始模型的依赖性低于其他几种自适应的算法,反演结果的稳定性较强,可以进一步提高正则化反演的效率.     

地下水对瞬变电磁反演解释的影响

反演的多解性是物探工作的致命弱点。在这一问题短时间内无法解决的前提下,综合其他学科的优势来进行反演解释,是一非常行之有效的方法。本文结合地下渗流力学、地下地质学等对瞬变电磁所采取的数据进行分析解释,阐明了地下水对岩石电阻率的影响规律,并着重分析了泥灰岩与灰岩在瞬变电磁反演断面上的可分性和可分条件。     

井中瞬变电磁全区视电阻率解释方法研究

针对井中瞬变电磁全区视电阻率定义求解和解释方法开展研究。依据均质全空间瞬变电磁场解析式计算不同采样时间或装置条件下感应电压响应,认识到响应值在电阻率取值范围的全区间或部分区间内保持单调性。根据"牛顿-下降"最优化原理,给出一种井中瞬变电磁全区视电阻率求解方法,通过均质全空间瞬变电磁场响应试算结果与介质真电阻率对比验证该方法的精度和适用性。建立典型径向及垂向分层地电介质模型,通过时域有限差分法模拟计算瞬态场感应电压响应,并进行全区视电阻率换算分析。结果表明:视电阻率是不同时刻电磁场作用范围内所有电性体的综合反映,由全时段视电阻率结果识别全空间介质的径向及垂向电性情况具有一定可行性,各时道的视电阻率值可作为推断地电结构等信息的依据。     

基于二维高密度电阻率勘探数据的三维反演及应用

针对目前高密度电阻率法所采用的数据处理方法主要是把地质结构体视为二度体进行二维处理,从而二维数据资料处理结果只是一种近似解释,其计算精度与反演效果达不到精确反演的要求,设计一种典型的地质体模型,利用有限单元法进行正演计算,对含有高斯随机误差的各测线二维正演数据分别进行高密度电阻率法二维、三维反演,对比二维和三维的反演结果.最后在广东某场地采用高密度电阻率法对岩层划分,利用二维勘探的数据进行二、三维反演.研究结果表明:三维反演受高斯随机误差的影响更小,反演结果与实际地质模型更接近;三维反演在岩层划分中更能凸显岩层分界面区域,岩层的空间分布及岩层、覆盖层的整体连续性能更真实地体现,与钻探结果吻合程度更好.     

电法测井的快速整体反演及其应用

为了提高电法测井反演的速度，提出了一种快速计算雅可比矩阵的方法，即只考虑每个待反演参数的微小改变对其所在层及其邻层视电阻率有影响，而对其他层的视电阻率没有影响，从而节约了大量计算雅可比矩阵的时间，其反演的时间与待反演的参数个数基本上呈线性关系。在计算精度基本不受影响的情况下，使得应用阻尼最小二乘方法反演侵入半径和原状地层真电阻率的速度提高了至少一个数量级。人工模型和实际资料的处理结果表明，该方法能有效地消除环境因素的影响，恢复地层的实际电阻率。同时也证实了阵列感应测井数据具有探测深度大、分辨率高的特点，在层不太薄（层厚超过1m）、侵入不太深（侵入半径不超过0．25m）时，深探测数据能较准确地反映地层真电阻率。     

利用高密度电法勘探盐边县格萨拉巨型滑坡堆积体

探讨盐边县格萨拉巨型滑坡堆积体的滑动面的位置和形态并分析滑坡稳定性进而为拟建高速公路路线设计提供有效的建议。采用高密度电法作为初步勘察阶段的滑坡勘探手段,建立不同采集装置的滑坡正演模型进行数值分析,初步了解拟勘探滑坡的视电阻率拟断面图并确定采集装置类型,结合地质调查结果,确定主要的反演参数,并根据反演结果修正正演模型进行对比判别。结果表明研究所用的反演参数效果较好,正演模型修正之后的视电阻率剖面与反演结果相匹配,使反演解译成果能够较好地反应出真实滑动面的位置和形态,且通过钻探得到了验证,计算结果表明该滑坡处于不稳定状态。最终,拟建高速公路的选线设计采纳了研究成果,对该巨型滑坡堆积体采取绕避措施,说明高速公路初步勘察阶段采用高密度电法作为巨型滑坡堆积体滑动面勘察的地球物理勘探手段是经济且高效的。     

M地区水淹层含水饱和度的计算

以M地区测井、实验、试油等资料为研究基础,采用电阻率相对值技术,建立含水饱和度解释模型,避开了地层混合液电阻率的求取.探讨了研究区两种含水饱和度计算方法的应用过程,计算结果表明电阻率相对值法计算的含水饱和度与岩心分析饱和度匹配较好.实际应用表明,该模型效果良好,为M地区饱和度计算提供了一种行之有效的方法.     

基于有限元数值模拟的直流电阻率激电找矿探测技术

为提高直流电阻率激电探测成果地质解释准确性,促进直流电阻率激电在找矿勘查中的应用,采用有限元数值模拟对直流电阻率激电探测中低阻高极化体、直立板状体、低阻盖层、高阻盖层、地表不均匀体、地形起伏及实测数据反演解释几个主要方面进行了对比分析。研究结果表明,对于实际勘查中的反演解释结果,根据异常特征构建地质体模型并进行有限元数值模拟,能快速判别探测目标体异常属性、区分虚假异常,提高定量解释水平及找矿勘查效率。研究结果为直流电阻率激电找矿探测定量解释提供了参考。     

基于非结构化网格的2.5维时间域激发极化反演

针对我国西部山区起伏地形,提出通过Delaunay三角剖分生成非结构化三角网格,可以根据实际地形情况,较好地控制节点分布,从而降低初始模型的几何离散误差。在反演中对模型参数施加最小光滑和背景约束信息,从而提高反演的分辨率。通过在正演求解稀疏矩阵时采用稳定型双共轭梯度法和计算系数矩阵时利用电位互换等手段,可有效地提高反演计算速度。最后,为验证反演方法的有效性,编制了跨平台反演软件;并对2例地电模型和实测数据进行反演试算。计算结果表明:反演耗费时间少,仅迭代3次拟合差便趋于稳定,反演结果能较好地刻画异常体形态。     

The ground and excitation levels of the trimer

The ground and excited levels of the trimer (⁴He)₃ are obtained theoretically by solving the 2-D Faddeev integral equation in the momentum space. In solving the equation, the direct matrix inversion method (I. M.) is used instead of the iteration and Pade approximate method. The new HFD-B Aziz potential ν is used in the calculation. You Pu: born in 1938, Associate professor     

基于CUDA架构并行算法的带地形AMT二维反演实现与应用

并行计算是提高音频大地电磁（audio-frequency magnetotelluric method,AMT）数据反演效率的有效途径。本文在统一计算设备架构(compute unified device architecture,CUDA) 下开展带地形的AMT数据二维反演并行算法研究,旨在利用GPU强大的计算能力及并行计算技术实现高精度、快速度的AMT数据二维反演。首先利用有限元和自适应正则化反演算法实现AMT数据二维反演的串行化计算;然后在PGI Visual Fortran+ CUDA5.5环境下编写基于CPU+GPU的CUDA并行代码,将正演中的频率循环、反演中的模型灵敏度矩阵计算和反演方程正则化求解部分进行并行化处理;通过不同复杂程度的理论模型正反演模拟验证了该并行算法的有效性和准确性。不同模型和不同模式下的数值模拟结果对比表明,基于CPU+GPU的CUDA并行算法相较于传统的CPU串行算法,在灵敏度矩阵计算和反演方程正则化方面耗时更少,加速比最高可达10倍以上。最后将该并行算法应用于某矿区实测AMT数据的二维反演中,取得了较好的应用效果。     

基于多种位场数据及3D反演结果的综合解释方法

近期位场梯度数据在勘探地球物理观测数据解释中的作用日益增加,随着计算机等硬件设施的进步,处理大型数据能力的增强,重力数据3D反演方法也越来越多的应用到重力数据解释工作中。本文提出结合重力数据、梯度数据及重力数据3D反演结果的综合解释方法,并将此解释方法应用到虎林盆地观测数据、梯度数据及3D反演结果的综合解释工作中。本文首先讨论单个棱柱体的重力数据、梯度数据的计算方法以及重力数据3D反演方法;然后,编制了matlab程序计算重力数据、梯度张量数据并处理重力数据的3D反演;最后,根据计算的数据详细讨论了梯度九分量之间的关系及重力数据反演结果的数值特征,提出结合梯度九分量数据及重力数据3D反演方法联合解释方法。模型试算及解释结果表明重力梯度数据对进一步提高勘探地球物理观测数据的解释精度及确定埋深较浅的异常源位置及水平尺寸等信息起到重要的作用,将其与重力数据3D反演方法联合解释可更准确的刻画异常源的分布特征。最后经过此综合解释方法将虎林盆地划分为8个构造单元。