频率测深正演问题滤波算法及人机联作反演

频率测深一维问题正演常用的计算方法是数值积分方法,但是这种方法由于利用求和代替积分使得计算速度很慢,精度也较低。本文提出了频率测深一维问题正演计算的一种有效方法——快速汉克尔变换数字滤波法,较系统地讨论了快速汉克尔变换滤波系数的解析计算及其在频率测深正演计算中的应用。最后应用上述正演计算方法进行了频率测深的人机联作反演的研究。     

磁偶源频率测深视电阻率曲线的正演计算

磁偶源频率测深法是一种有效的地球物理勘探方法。多年来,人们对这种方法的正演理论和反演解释进行了研究。文中利用快速汉克尔数值滤波算法计算了层状介质上直立磁偶极子的电磁场响应以及频率测深视电阻率的振幅响应,这种计算方法与常用的数值积分算法相比,具有计算速度快、精度高等特点。同时详细地讨论了磁偶源频率测深的特点及视电阻率理论曲线的特征。     

直流电阻率三维正演的预条件共轭梯度法研究

三维电阻率法对反演的精度和速度的要求越来越高,而正演是反演的基础,因此直流电阻率三维正演计算的速度和精度是三维电阻率反演实用化的关键。这里利用对称超松弛预条件共轭梯度法(SSOR-PCG),求解有限差分法离散生成的大型稀疏线性方程组,预条件矩阵的选择大大降低了系数矩阵的条件数,结合矩阵的一维非零元素压缩存储模式,使得正演计算速度得以提高,而内存占用量明显减小。在直流电阻率三维正演中采用异常场法,提高了电源点附近的解的精度。利用编制的有限差分正演程序,对两层模型、垂直接触带模型和低阻异常体模型进行了数值模拟,计算结果表明该算法是可行的,且可以明显提高正演计算的速度和精度。     

频率测深资料数字解释

针对我国目前存在的频率测深资料解释问题,本文提出了数字解释技术.首先实现了频率测深曲线的数字滤波正演计算,较数值积分大大缩短了计算时间.在此基础上,进而完成了广义逆矩阵方法的自动反演解释,取得了良好的效果.     

起伏地形重磁三维快速正演计算

正演是反演技术的基础,正演速度和求解反演问题的系数矩阵存放一直是起伏地形下重、磁三维反演的关键技术问题。这里提出了一种起伏地形下重磁快速正演计算方法,其计算原理是根据反演在垂向的剖分层数,利用水平地形正演计算形成二个不同大小刻度标尺矩阵,然后在模型空间,使用分段线性插值的方式,直接计算出起伏地形观测点的正演值。该方法的主要特点是在保持很高的计算精度下计算速度可提高二倍,且节省计算内存,适合起伏地形下重磁三维反演技术研究。     

基于全区视电阻率的瞬变电磁一维Occam反演中雅克比矩阵的解析算法

在目前的工程应用中,瞬变电磁测深一维反演结果依然是地质解释的重要依据。为提高一维反演的计算效率和计算精度,解析推导出一维正演的雅克比矩阵就非常必要。本文详细推导了中心回线装置下瞬变电磁一维正演全区视电阻率的雅克比矩阵的解析计算公式,并通过Occam反演测试理论模型和实测数据,证明该推导正确,值得采用。     

基于LSQR算法的二维声波方程频率域正演模拟与数值实现

频率域正演是频率域全波形反演的基础,高效准确地完成频率域正演计算是目前该领域的一个热门问题。这里利用优化17点差分格式和PML完全匹配层构建波场,并使用LSQR算法求解波场,通过在简单层状模型和Overthrust模型上计算验证可知,LSQR算法可大幅提高频率域正演的计算速度和模拟精度,从而为此方向的研究提供一个可行的技术手段。     

电阻率测深自动直接解释法的若干技术问题

本文对应用Koefoed法自动直接解释电测深曲线中的数据采集、归他采样值的计算、反演计算、正演计算、图象处理以及计算精度等问题作了讨论.研制了一套完整的电测深自动解释软件.     

电阻率测深正反演交互解释系统

本文对应用Ｋｏｅｆｏｅｄ法自动直接解释电测深曲线中的数据采集、归化采样值的计算、反演计算、正演计算、图象处理以及计算精度等问题作了讨论。研制了一套完整的电测深自动解释软件。     

电阻率测深自动直接解释法的若干技术问题

讨论了应用Ｋｏｅｆｏｅｄ法自动直接解释电测深曲线中的数据采集、归化采样值的计算、反演计算、正演计算、图形显示以及计算精度等问题，并介绍了由作者开发的一套完整的电测深自动解释软件。     

Zohdy方法在三维电阻率测深数据反演中的应用

三维电阻率测深数据采用Zohdy方法近似反演,正演计算采用有限单元法,反演初始模型由测量视电阻率数据给出,通过比较实测视电阻率值和预测模型计算的视电阻率值对数差来修改模型网格电阻率.为了解决任意电极距测深数据的反演,采用大、小双网格剖分,大网格反映地下电性分布情况,小网格用于实际有限元正演计算.文章对加5%随机噪声的模型理论视电阻率测深断面数据和实测视电阻率测深断面数据进行反演.结果表明该反演方法需要的内存小、计算量少,并能够取得很好的反演效果.     

基于两步变异差分进化算法的激电测深一维反演

激电测深一维反演是一个多参数非线性优化问题。本文基于改进的两步变异差分进化全局最优化算法实现了激电测深的一维反演。传统的差分进化算法包含变异、交叉和选择操作,变异为单步变异。本文提出的两步变异法分步将最优个体与随机选取的两个个体经变异后产生新个体。加强了最优个体的影响度,提高了全局寻优能力。通过模型试算结果表明,两步变异法比传统方法寻优能力更强。利用等效电阻率法加载极化率参数,通过数字滤波算法可快速正演计算层状模型表面激电测深电阻率曲线,在此基础上应用两步变异差分进化算法不断变异产生新个体,正演计算电阻率与观测值进行拟合,选择适应度值趋近于最大适应度值的个体作为反演结果。本文反演方法操作简便,计算速度快。通过对H型和KH型地电模型进行计算,得出本反演方法有较高的拟合精度。     

2.5维复电阻率反演并行算法设计

近年来2.5维复电阻率反演取得了一些理论和应用方面的成果,但计算过程复杂、计算量巨大,特别是正演计算过程中计算量大导致了反演速度较慢。这里在2.5维复电阻率反演计算基础上,将反演算法中的正演过程的排列和频率循环分割成多个规模较小的程序段,并发执行正演过程中的每个排列和频率循环的求解计算过程,并根据并行计算的节点数目或者计算核心数目,动态调整每个节点或者计算核心的计算量,从而提高了算法的计算效率和适应性。     

用电子计算机自动解释电测深曲线的尝试

电测深曲线的定量解释工作历来是使用传统的量板法。这种人工对比曲线的方法带有人的主观因素。同一条实测曲线由于解释人员的不同,解释结果常有较大差异,精度较差,效率也低。近十年来,国內外发展了电测深曲线的自动解释技术。使电测深资料的解释实现了数字化。在这方面已有不少资料。作者根据现有资料,编制成解电测深曲线正演问题和反演问题的程序。正演计算时根据核函数的递推公式求得任意层数的核函数,再用数字滤波法求得理     

磁异常和梯度的频率域三维成像方法

三维反演在磁数据定量解释中具有重要作用。常用的空间域三维反演方法通常需要大量的正演和反演计算,因此对大规模数据的反演效率较低。三维成像是另一种定性和定量解释磁数据的重要方法。文中给出了一种磁异常与梯度三维成像的频率域迭代方法,该方法可以提高成像效率,适用于大规模数据的快速成像。笔者推导了磁总场异常和异常梯度频率域正演公式和成像公式,并将一种深度尺度因子引入成像公式中,提高了深度精度;笔者采用了迭代优化算法,减小了拟合误差,进一步提高了成像精度。通过理论模型数据试验和中国新疆某金属矿床实测数据,验证了本文方法的有效性、可行性。     

在大地电磁二维Occam反演中求取拉格朗日乘子方法改进

Occam反演具有可靠性和加上各种先验信息可行性的优点,但是因为求取拉格朗日乘子需要很长时间,所以前人提出了几个直接选取方法。笔者改进了二维大地电磁测深Occam反演中求取拉格朗日乘子的方法。在搜索过程中加上观测数据拟合差与其期望值的比较。这样,在Occam反演的拟合差下降阶段的最后迭代和模型光滑阶段,能够排除多余的正演计算。这种改进对反演结果没有影响,但却使得反演速度加快并且更符合于Occam思想。根据模型实验和野外数据处理,可以减少正演次数20%~50%。     

基于磁导率的频率域航空电磁法双频反演方法

以往的频率域航空电磁法数据反演方法中,通常在模型中以真空磁导率代替介质磁导率,忽略了磁性层影响。但在高磁性区域,这样做将会导致计算电阻率值偏高。针对这一问题,首先根据层状介质模型研究了考虑磁导率和电阻率因素条件下的电磁场正演计算方法,并且分析了不同磁性介质对频率域电磁响应的影响规律。在此基础上,提出了磁导率和电阻率的双频反演方法,电阻率是通过2个频率的虚分量之比计算得到的,计算结果更加可靠,提高了磁性区域视电阻率填图的精度。通过对模型计算和实测数据进行反演计算,表明这种方法具有较高的反演精度,为强磁性区域航空电磁法视电阻率填图提了供一种行之有效的反演方法。     

任意角度瞬变电磁测深正演问题计算方法

针对目前比较普遍的瞬变测深晚期场计算精度较低的情况,本文依据频率域电磁场实虚分量的奇偶性及付氏变换的导数性质,利用余弦变换实现了由频率域到时间域电磁场的转换。计算表明,本算法计算晚期瞬变电磁及视电阻率响应具有较高的精度。文中重点计算了单一偶极供电,多点测量的任意角度瞬变测深的正演问题,并对其规律性作了总结。     

可控源音频磁大地电流法一维正演及精度评价

研究发现，两种目前比较常用的可控源音频磁大地电流法（ＣＳＡＭＴ）的正演计算方法，即简化积分法和二次插值法之间计算结果差异较大。为了正确地评价这两种方法的计算精度，本文提出了利用直接积分进行，ＣＳＡＭＴ一维正演的算法，即通过将任意角度频率测深￣［１］正演结果沿发射源方向直接进行数值积分，得到有限长度发射源的ＣＳＡＭＴ正演结果。大量计算和精度评价发现，在这些算法中，简化积分法和本文提出的直接积分法计算结果吻合很好，是精度较高的算法，而二次插值法与前两种算法结果差别较大，是精度较低的算法。     

改进的重力模拟退火反演研究

通过对模拟退火方法中常用的VFSA算法的分析，针对其在模型扰动及退火计划上存在的缺陷，针对性地作了更合理的改进VFSA算法，提高了VFSA算法的计算速度与稳健性，并结合重力反演过程中对模型的修正，克服了常规重力反演方法所存在的等值效应。研究中还把重力正演计算的单元体改为垂直侧边梯形单元，并合单同一层位的垂直侧边梯形单元成为一大单元再作计算，加快了正反演计算速度并提高了正反演精度。通过模型试验及实际资料的处理，验证了作者所提出的改进的模拟退火重力反演方法的有效性。