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

上接119页 9 反演方法推导中的有关公式 为了对复场展开式作简化表示,将会多次用到2个递推公式,这里预先导出.     

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

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

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

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

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

上接211页13 n个无限薄板叠加场反演13.1 复场展开以Mj代替2Mmj和2fMgj,以F(S)代替G′(S)和T(S),得复场F(S)=nj=1(Mj)/(σj-S)=-nj=1(Mj)/(S-σj).(13.1)展开上式,即     

二度体的重力张量有限元正演模拟

介绍了重力梯度张量,并将有限单元法应用于二维重力梯度张量的正演计算.为了验证有限元正演方法的精度,对截面为矩形的两个二度体组合模型进行有限元正演模拟,结果表明正演曲线与理论曲线形态一致,拟合情况好.通过对截面形状不规则、密度分块均匀的二度体进行正演模拟,说明有限元法可通过网格剖分来逼近不规则目标体的边界,并对剖分单元赋予不同的密度值来实现对复杂二度体的重力张量的正演模拟.     

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

上接 38页5　复场直接积分法和板状体复磁场公式5 .1　薄板状体复磁场公式5 .1.1　水平薄板设 :水平薄板厚为Ｄ、宽为 2ｂ ,左端点为σ1,右端点为σ2 =σ1+2ｂ ,板内变点σ =σ1+ｌ,0≤ｌ≤2ｂ。该板复场Ｔ(Ｓ)相当小水平圆柱体的复场沿水图 6　水平薄板断面示意平方向积分 (图 6)。在 ( 1 4 7)式中 ,取Ｍｍ =ＤｄｌＭｓ,　σ0 =σ1+ｌ ,( 5 1)则水平薄板状体复磁场为Ｔ(Ｓ) =∫2ｂ0μ04π2ＤＭｓｄｌ{Ｓ- (σ1+ｌ) } 2 =μ04π2Ｍｍ1σ1-Ｓ- 1σ2 -Ｓ ,( 5 2 )Ｍｍ=ＤＭｓ。 ( 5 3)5 .1.2　有限延深倾斜薄板设板宽为 2ｂ( Ｉｍσ1)、…     

线源频率电磁测探二维正演（二）

随着频率电磁测深（包括可控源声频大地电磁测深ＣＳＡＭＴ）的深入发展,必须要进行三维情况下的研究。线源频率电磁一维天演模拟研究为最终实现三维正演打下了基础。     

基于FlexPDE考虑退磁作用的有限元法磁场正演

FlexPDE是PDE Solutions Inc.开发的用有限元法求解偏微分方程的应用软件,是—个脚本有限要素模型建立工具和数字处理器.基于FlexPDE平台,这里用有限元法求解静磁场拉普拉斯方程,并模拟二度圆柱体、板状体模型的磁场.模拟结果表明:强磁性体的退磁作用使异常幅值变小、形态畸变,磁性体内部磁化强度大小、表面磁荷密度分布不均匀,磁化强度方向往磁性体长轴方向偏转.     

板状体瞬变电磁正演计算系统

为了便于研究地质体中常见层状地质结构的瞬变电磁响应特点,开发了针对板状体的正演计算软件。软件设计采用了面向对象与面向过程的方法,通过动态链接库将两者结合起来。系统界面友好、稳定,可以实现各种板状体的正演计算。二维计算采用有限差分的方法,实现了不同厚度、倾角、电导率的板状体及其组合的计算,计算结果以直观的图形曲线形式展现出来,并可利用其他成图软件进行深加工与研究。     

TI—59计算器多边形截面二度体ΔT、ΔZ正演程序

TI—59可编程序计算器在我国已有相当大的数量,为了充分发挥它在野外生产、科研、教学中的作用,本刊将择优选登该计算器的常用计算方法,物探、化探、物探测量方面的优化程序。     

基于泊松方程的空间波数混合域二度体磁异常数值模拟

正演数值模拟是反演成像的基础。为了实现磁法勘探精细化反演成像与定量解释,本文利用把一个大问题分解为多个小问题的思路,提出一种基于泊松方程的高效、高精度空间波数混合域二度体磁异常数值模拟方法。该方法利用傅里叶变换把二度体磁位偏微分方程转换为一维常微分方程,并采用基于二次插值的一维有限单元法求解该方程,进而通过反傅里叶变换得到空间域磁异常。在模型算例中,分别设计截面为矩形的常磁化率和变磁化率二度体模型,针对本文算法的计算精度和计算效率进行了验证。模型算例结果表明：该算法计算精度高,相对误差绝对值均小于1%;计算速度快,网格节点剖分2 501×2 501的模型模拟时间为4.18 s;适用于任意复杂地形模型。     

线源频率电磁测深二维正演(一)

随着频率电磁测深(包括可控源声频大地电磁测深 CSAMT)的深入发展,必然要进行三维情况下的研究。线源频率电磁测深二维正演模拟研究为最终实现三维正演打下了基础。     

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

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

线源频率电磁测深二维正演(二)

随着频率电磁测深（包括可控源声频大地电磁测深ＣＳＡＭＴ）的深入发展,必然要进行三维情况下的研究。线源频率电磁测深二维正演模拟研究为最终实现三维正演打下了基础。     

二维重磁资料联合正演模拟

重、磁异常的人机联作模拟反演(也称人机联作选择法),是根据选定地质模型及重、磁异常的正演计算公式,由给定的初始模型参数计算出理论场值,并与实际观测场值进行对比,根据对比中确定的不一致性或残差,依靠解释人员的经验来修改模型形状及物性参数,重新计算理论场值进行下一次对比,反复多次直至达到满意为止,人机联作选择法整体为一套计算机软件,并依据支撑条件的不同可分为在人机联作工作站上使用或在一般的微机上使用。人机联作选择法既可应用于二度体重、磁异常,也可以应用于三度体重、磁异常的模拟计算解释。     

矿井电穿透法二维正演数值模拟

通过研究巷道间电穿透二维正演数值模拟理论及复杂模型正演实例分析, 结果表明工作面内顺层电穿透技术, 仍属体积勘探, 探测的有效深度是一个深度段。不同建场方法有不同的建场特性且能分别突出各目的层 (煤层内或顶底板围岩内) 的异常, 为解决采煤工作面附近的具体问题使用相应的方法提供了理论依据。     

地表大回线源在任意层状介质中产生磁场的计算

在地面核磁共振数值模拟中,地下介质中任意一点处磁场快速计算非常重要。作者然对水以状介质,推导出利用递推方法求得任意一层中电磁波的透射和反射系数的简单方法,从而为精确、确定地计算位于地表的大回线源在地下任意一点处产生的磁场创造了有利条件。     

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

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