基于快速汉克尔变换算法的CSAMT一维正演

水平层状介质上的电偶源产生的电磁场数值解是含有零阶和一阶贝塞尔函数的无穷积分,这类积分实质上是汉克尔变换式。通过3种方案的计算实例,说明快速汉克尔变换滤波系数的选取对正演结果的影响。计算结果显示:加长滤波系数能明显提高计算精度,使正演结果真实地接近实际;相反,滤波系数区间过短会出现场值的畸变,视电阻率会受到较大影响。因此,合理选取滤波系数可以提高计算精度和速度。     

Guptasarma线性滤波算法在CSAMT一维正演中的应用

CSAMT法一维正演问题,其核心运算是对汉克尔变换的求解,目前主要使用数值滤波方法进行求解,求解过程中使用给定长度的滤波系数进行累加运算.滤波系数的选取,影响着正演计算的精度和速度.这里使用Guptasarma线性滤波系数,并给出具体实现方法.经过正演模型验证,发现这套方案可以在较短的滤波系数长度基础上,与正演模型拟合...     

电场差分法一维正演研究

通过研究电场差分法的一维正演理论，建立了一维正演的数学模型，导出了相应的理论计算公式和算法，编制了计算机程序，对一些典型的均匀、二层和三层可极化地电断面的计算结果进行了分析，概括了Ｐ参数的异常性态。在此基础上，对电场差分法的工作方法技术、资料解释和实际应用提出了建议。     

时间域激发极化法非结构化三角网格有限元正演模拟

这里分析了传统网格剖分方式的不足,提出采用非结构化渐变三角网格对求解区域进行离散,采用固定带宽LDLT分解法求解线性方程组。设计了非结构化渐变三角网格自动剖分程序,并针对固定点源测深装置的激发极化2.5维有限元正演模拟程序。通过理论模型试算,验证了程序的正确性与有效性,并对正演模拟的速度和精度进行了初步评价。     

基于Matlab与Fortran混合编程的一维fCSEM正演可视化

这里介绍了Matlab与Fortran混合编程在可控源音频大地电磁法(fCSEM)一维正演中的应用,实现了一维fCSEM正演的可视化,建立了操作界面窗口,加入了可以调节H_y取值的功能,使曲线更接近真实值。并以双极源全区视电阻率计算为例,验证了该方法的可实用性。     

微机网络并行计算及2.5维CSAMT正演的并行实现

本简要介绍了如何建立一个简单、易行的微机网络并行计算的开发环境。在该并行环境上完成了２．５维ＣＳＡＭＴ有限单元正演的并行程序，并且获得较好的速比。     

CSAMT一维全区反演

由汉克尔变换和连分式求解水平层状大地电偶源的水平电场和垂直磁场,然后由电场和磁场求得卡尼亚视电阻率,不进行近区和中间区校正,直接用阻尼最小法对近区和中间区发生畸变的卡尼亚视电阻率反演,得到地下电性参数。模型试算结果表明:虽然近区和中间区的视电阻率不能正确反映地下电性垂向变化,但由其反演仍能得到正确的结果。     

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

通过研究工作面两巷道间电穿透建场理论和一维正演数值模拟,结果表明6种电穿透技术对工作面顶底板内的电性变化特别敏感,而对工作面煤层内部的电性变化反应不敏感。为该方法探测工作面顶底板内的含水、导水构造提供了理论依据。     

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

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

MT二维正演的混合剖分算法及其在线源CSAMT正演计算中的应用

在矩形有限元和三角形有限元方法大地电磁测深（ＭＴ）二维正演计算的基础上，设计了矩形单元粗剖分、程序自动进行三角形单元细剖分的混合剖分算法程序，使得ＭＴ二维模型正演的计算速度大为提高。该方法同时被应用于线电流源变频电磁测深（ＣＳＡＭＴ）ＴＥ极化方式的二维正演计算。文章给出了在不同模型上计算的ＭＴ和线源ＣＳＡＭＴ的响应，分析了计算ＣＳＡＭＴ响应时，计算点离源的距离对视电阻率和相位值计算的影响，比较了ＭＴ和ＣＳＡＭＴ网格剖分的特点，提出了计算ＣＳＡＭＴ响应的改进的边界条件设计     

基于DCT有限元二维重力张量的正演

介绍了国内、外张量测量技术及数据处理技术的发展情况。把离散余弦变换与有限单元方法结合起来,进行重力张量的正演计算。首先,利用有限单元法求解重力场的一阶导数,然后利用离散余弦变换求解重力张量分量。这样既保留了有限单元法对复杂模型的正演优势,又避免了利用有限元法直接求解重力张量时,因多次插值引起的数据量过大的问题。模型试验证明,用基于DCT的有限元法进行重力张量的正演,不但具有很高的计算精度,而且相对于利用有限单元法直接正演,其计算速度得到了显著的提高。     

电法二维有限元正演计算的若干新结果

文章揭示出最优化选取波数的方法给电法二维有限元正演计算带来了显著的效果，并经大量试算发现了误差随波数个数及网格剖分参数变化的规律——变化曲线上存在突变点，突变点之后误差趋于平稳变化。这一发现使得我们可以在保证精度的前提下大大节省计算时间，这将对正演更充分地应用于反演有着重要意义。     

二维弹性波动方程F—K正演与偏移

根据Ｈｏｏｋｅ定律，消除应力矢量后得到弹性波动方程位移矢量和波势矢量的非耦合形式递推解，在这个解答中各种类型的波相互分离，但在界面上满足能量守恒定律。正演可变步长延拓；偏移是将震源脉冲和地面记录波场同时下延，然后在每个深度作时间零延迟互相关来完成。由于延拓和成像在频率波数域（Ｆ—Ｋ）中进行，因而算法快速、偏移和正演可方便地扩展到三维空间，能处理各向同性任意非均匀介质中的弹性波。     

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

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

基于球体层状介质模型的大地电磁正演

目前大地电磁测深法的一维正演理论,是基于平面波垂直入射水平层状介质的假设模型。但由于地球是一个球体,因此有必要研究基于球状介质模型的大地电磁正演理论。这里详细推导了基于球体层状介质模型的大地电磁正演公式,计算了若干理论模型。通过同基于水平层状介质模型的大地电磁正演结果对比,验证了正演公式的正确性。同时,指出当探测周期增加到上万秒时,阻抗相位会增大;而当探测周期增加到数十万秒时,视电阻率会减小。     

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

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

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

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

基于八叉树结构的可视化三维位场正反演

本文采用八叉树数据结构提出了一种新的三维可视化人机交互正反演方法,应用该方法可以构造复杂的三维位场模型,计算出其正演场值,通过人机交互界面,使模型与数据在整个正反演过程中可视化、形象化。该方法不仅可以进行均匀物性分布正反演,也可以进行非均匀物性分布正反演,而且在多种地球物理参数可视化联合反演中有着很好的应用前景。     

基于OpenMP一维CSAMT正演并行计算研究

电偶极子场源在水平层状介质中传播的电磁场的数值解是通过求解含有零阶和一阶贝塞尔函数的积分(汉克尔变换)所得,可控源音频大地电磁法(CSAMT)一维正演计算就是计算各个频率的视电阻率,使用传统串行方式对频点逐一计算,速度相对较慢。本文旨在使用Intel Visual Fortran 2013和VS2012软件,编写基于OpenMP的CSAMT一维正演并行计算程序。与一维CSAMT串行程序的对比:相同的计算量下,并行程序运行时间少于串行运行时间,计算量越大并行加速比越明显;验证了并行程序的可行性,为CSAMT多维复杂模型正演计算提供了并行研究基础。