电和磁信号中有噪声时估计地球阻抗函数问题

一、引言我们介绍一种当输入和输出信号含有相加性噪声时求电路的传递函数或地球的阻抗函数的渐近无偏估计值的方法。基本的输入-输出关系是 e(t)=Z_(ei)(t)*i(t) (1.1)其中*是褶积运算,Z_(ei)(t)是脉冲响应,图1a中表示的是系统。当测量i(t)和e(t)的器件如图1b所示含有噪声时,方程(1.1)中的Z_(ei)(t)应按下式估算:     

有限元法在大地电磁场正演计算中的应用及改进

本文讨论了有限元法在大地电磁测深二维正演问题中的应用.研究了现有文献中的方法,提出了以矩形单元为主,结合采用三角形单元的方案,作为对现有方法的改进,并以潜伏倾斜界面为例,给出了具体的计算公式.     

华北平原GB大地电磁剖面的二维解释

大地电磁测深资料的二维解释主要有两种方法,即二维反演拟合和二维数值模拟。作者参考二维三次样条快速反演方法的思路,在一维样条插值的基础上作形式推广,建立了三维三次样条插值的一般表达式及其计算过程,给出了三维三次样条插值的反演方法。二维数值模拟是在一维形式化解释的基础上构制二维初始模型,然后计算各测点上的MT响应,与实测值拟合,以便求得最佳拟合结果。所谓一维化解释是把剖面上各测点测得的纵向视电阻率ρ_⊥和横向视电阻率ρ₁₁分别按一维模型曲线作解释,求得地电剖面的二维模型和参数。用所述方法对华北地区的GB剖面的大地电磁测深资料进行试算,获得了地下导电率分布的基本情况。     

用数字滤波方法计算电、磁偶极子的场

在计算水平层状地层上由电、磁偶极子产生的电磁场分量的公式中,涉及一积分限为无穷大的积分。其被积函数是两个函数的乘积。一个称为核函数,它和地电断面的参数、频率及积分常数λ有关,在数字滤波计算中被视为输入函数。另一个是不同阶数的第一类贝塞耳函数。由于它们具有振荡性质且衰减很慢,故用数值方法求这一积分是很费时间的。在用数字滤波方法时,这一部分构成数字滤波器系数。由于此滤波器系数只与装置的类型有关,不随地电断面的性质变化,故对各种装置只需计算一次,因而,使用起来比较方便。与数字积分方法相比,数字滤波方法大约可以提高功效四倍。本文包括三部分内容:第一部分总结了文献中用数字滤波器计算垂直磁偶极子场的方法;第二部分给出了利用已有的数字滤波器计算水平电偶极子的电场;第三部分介绍具体实现步骤及计算结果的对比。     

处理大地电磁资料的瞬时谱方法

大地电磁资料处理方法大多数都是对一段记录做傅氏变换,构成功率谱后再求张量阻抗元素。本文通过对记录上的实信号做窄带复数滤波或窄带递归滤波和希尔伯特变换,转换成窄带复信号,再由复信号的瞬时谱求出张量阻抗元素。利用这种方法对实际资料处理的结果与常规方法处理的结果基本一致。由于采用了快速褶积实现希尔伯特变换,因而处理时间也与常规方法大体相当。     

有限元法在大地电磁场二维正演中的改进

本文讨论了有关有限元法在大地电磁测深二维正演中的应用问题,论述了现有文献中采用的网格剖分方案。矩形网格具有易于使公式规则化和便于进行计算的优点,但失去对界面几何形状较好的适应能力;再有,由于在矩形单元上采用双线性插值公式,因而场相对z的偏导数仅是y的(线性)函数,而沿z方向是常数,这就影响确定辅助场值的准确程度。基于这种考虑,作者采用三角形-矩形单元综合剖分方案,即在地面和界面(特别是倾斜界面)处,以对角线将矩形划分为两个直角三角形,然后采用拉格朗日二次插值公式,计算泛函,再按矩形单元构成组合的刚度矩阵,最后合成总刚度矩阵。这样,既可较好地适应各种几何形状的界面,又可提高计算界面附近场值及地面处辅助场值的准确程度。这时,辅助场值的表达式是y、z的线性函数,并可表示成插值点处函数值的简单的线性组合。这在反演计算中会带来很大的方便。文中给出了相应的计算公式。     

油气田地区地球化学作用形成的地电异常的勘探技术在匈牙利的应用

文中编译了有关烃类聚集后出现的地球化学过程及其周围引起岩石变化的理论假说和应用电磁方法检测相应电异常的可能性,给出了在匈牙利进行这项工作的经验及实例,着重介绍了为此专门设计的WEGA—D自动控制电磁勘探系统的性能、工作原理及实际效果。     

二维地电构造大地电磁测深资料的解析方法

对于包含垂直界面的二维大地电磁测深资料的广义逆反演问题,采用解析法求偏导数矩阵的方法,在一定程度上克服了用差商法求雅可比偏导数矩阵带来的精度低、反演参数少、用机时间长等缺点。使用Rodi的模型及参数按本方法进行了试算,证明程序是正确的。应用理论模型数据对二维大地电磁测深广义逆反演方法进行验证,表明该方法是可行的。利用所编制的MT二维广义逆解析法反演程序,对华北某地MT试验剖面的E_x极化和H_x极化资料进行解释,并做了联合反演,取得了良好的效果。     

大地电磁场二维正演计算中若干问题的讨论

在用有限元法研究二维构造上大地电磁测深曲线时,采用Crout分解法取得了满意的结果。使用该法不仅可利用一维紧缩存储,而且可按列块分别求解,从而充分地利用了计算机的存储单元,而无须顾忌网格过多所造成的计算困难。本文讨论了网格设计的原则,并就地垒模型进行了计算。试算结果和文献[5]、[6]中提供的结果是一致的。文章最后提出了解释实际大地电磁测深曲线应该遵循的程序。