起伏地形条件下的磁异常二维聚焦反演及应用

地形对磁异常的反演和解释有很大的影响,特别是磁异常体贴近地表分布时更是如此,为了降低这种影响,笔者推导出了基于有限延深的二度厚板状体的带地形的磁异常正演公式,并提出了一种基于起伏地形的地下网格模型剖分的新方法。同时,聚焦反演对物性体有很强的聚焦效果,如果能够预先判断磁异常主要由地下浅层的强磁性体所引起,则可以考虑采用聚焦反演的方法得到磁性体的位置,显然这种方法非常适合铁矿体的反演。笔者将带地形的二维磁异常正演与聚焦反演结合起来,进行了相应的模型试验,并与水平地形下的聚焦反演及起伏地形下的光滑反演进行对比,验证了这种结合的有效性。最后,将该方法应用于新疆特克斯的磁异常剖面,对该剖面的磁性体分布做出了合理的解释。     

复杂条件二维重力场及重力张量场空间波数域正演方法

随着传感器技术的发展,重力场与重力张量场测量技术发展迅速,为实现地下密度分布精细反演提供了数据保障。正演是反演的基础,解决任意密度分布复杂地质体重力场与重力张量正演高效、高精度计算问题,是实现重力高效、精细反演、人机交互反演解释的关键。针对起伏地形和任意密度分布这种复杂条件下二维重力场及重力张量场高效高精度正演问题,这里提出了一种空间波数混合域正演算法,其关键环节包括:①结合新的矩形二度体组合模型波数域表达式和一维Gauss-FFT算法,提出了一种任意密度分布和起伏地形下重力场及重力张量高效、高精度正演算法;②采用新的二维正演算法,计算观测最高点和最低点之间多个不同高度水平网格重力场及重力张量,结合三次样条插值方法,实现了起伏地形上重力场及重力张量场高效、高精度正演。模型算例结果表明,新方法具有高效、高精度的显著特点。     

基于非结构化网格的大规模磁测数据正则化反演研究

磁测数据解析正演的内存需求与模型规模、观测点数成正比,以其为依托开展大规模磁测数据反演时,计算成本较大。文章利用非结构化网格有限单元法进行磁测正演,分析了地表观测网格、边界节点数和局部网格密度等因素对正演计算精度的影响;在此基础上,基于灵敏度矩阵隐式存储的高斯-牛顿法求解正则化反演目标函数,实现了大规模磁测数据三维反演。理论模型试算和铀矿区域实测数据案例表明,文章所述方法可保证正演计算的精确度,提高了计算效率、减少了反演的内存消耗,可用于复杂条件下大规模磁测数据的三维反演。     

利用小波变换快速反演大型磁数据

我们应用小波变换反演大型磁数据,恢复三维磁化率分布问题,快速小波变换通过给出小波系数的门槛值来产生灵敏度矩阵表达式,在波域进行正演模拟计算,通过小波域内稀疏矩阵矢量的乘积运算,和压缩后的矩阵来实施快速正演模拟,所用ＣＰＵ的时间的减少与矩阵压缩在正比,而且,我们采用对数阻尼法,用共轭梯度解算器来获得反演结果,这种求解手段的组合,直接减少了正演模拟的时间,进而加快了反演速度,产生了更的三维磁反演算法。     

起伏地形条件下FPS不同解释方法对比

固定点源激电测深(FPS)是一种有效的激发极化测深勘探方法,但在数据解释中,常规的定量解释法只适应较为简单的地电条件,在起伏地形条件下其解释效果会受到较为明显的影响。针对常规的定量解释方法所存在的不足与缺陷,这里从有限单元法正演模拟入手,采用最小二乘法进行反演计算,实现了在起伏地形条件下的对地电模型进行模拟计算,并将圆弧交汇-相对强度计算结果与最小二乘反演计算结果进行对比。结果表明:最小二乘反演方法能适应起伏地形对模型进行反演计算,其反演效果优于圆弧交汇-相对强度法计算结果。     

起伏地形下偶极-偶极激电测深二维反演软件开发及应用

在有限元点源二维正演模拟中,采用第二类齐次边界条件结合相应的网格剖分技术,并将二维正演模拟程序引入到偶极-偶极激电测深二维反演中,编制相应反演程序并进行集成,形成集数据输入、反演计算、数据成图和结果输出为一体的二维反演软件。通过对起伏地形下的两个理论模型合成数据及一例实测数据的反演,反演结果表明反演软件有效。     

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

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

重磁异常数据三维人机联作模拟

在研究三角形多面体模型重、磁异常三维正演和反演技术的基础上,吸取人机交互正演模拟的优点,实现了三角形多面体模型重、磁异常数据三维人机联作反演。通过研究三角形多面体模型节点偏导数的计算方法,对目标函数进行线性化处理,形成了计算机自动迭代修改模型体的技术。利用计算机图形技术,在三维空间显示重、磁场和模型体,开发了模型的交互修改技术,使数据解释过程中,可以结合已知信息及人的推断和经验,完成重、磁异常数据的三维模拟,减少了数据解释结果的不确定性。     

起伏地形上规则二度体复重磁场正演和直接反演

"起伏地形上规则二度体复重磁场正演和直接反演"由4篇连续文章组成。第一篇和第二篇阐述复重磁场正演问题,主要采用前苏联学者的方法,在复数域内,严格而简明地推导出了规则二度体及其组合体复重磁场表达式。这些复场表达式与直角坐标系内对应形体重磁场表达式相比,具有非常简单的形式。第三篇和第四篇阐述复重磁场直接反演问题,是本系列文章的核心,从规则二度体及其组合体复重磁场表达式出发,完全在复坐标系内,采用独创的线性化方法,建立起求解中间变量的线性方程组,再结合复高阶方程解法,成功地解决了起伏地形上(或者观测线处于二维场源任意方位甚至封闭时)几乎所有规则二度体重磁场的直接反演问题,推动了重磁场解释理论的发展。     

起伏地形上规则二度体复重磁场正演和直接反演

"起伏地形上规则二度体复重磁场正演和直接反演"由4篇连续文章组成.第一篇和第二篇阐述复重磁场正演问题,主要采用前苏联学者的方法,在复数域内,严格而简明地推导出了规则二度体及其组合体复重磁场表达式.这些复场表达式与直角坐标系内对应形体重磁场表达式相比,具有非常简单的形式.第三篇和第四篇阐述复重磁场直接反演问题,是本系列文章的核心,从规则二度体及其组合体复重磁场表达式出发,完全在复坐标系内,采用独创的线性化方法,建立起求解中间变量的线性方程组,再结合复高阶方程解法,成功地解决了起伏地形上(或者观测线处于二维场源任意方位甚至封闭时)几乎所有规则二度体重磁场的直接反演问题,推动了重磁场解释理论的发展.