有限单元法与有限差分法在MT一维正演模拟中的对比分析

首先从电磁场所满足的麦克斯韦方程组出发,介绍了大地电磁测深正演的基本理论,并针对一维大地电磁模型加以讨论.运用有限单元法及有限差分法分别推导了大地电磁测深一维正演算法,并运用Matlab 7.0软件编写了相应的程序.为了检验这两种一维正演算法的准确性,设计了均匀半空间模型和层状介质模型,并给出了由本文程序、解析解得到的相应结果和图件,从而对正演结果进行对比分析.结果表明:两种方法的正演结果均真实地反映了模型的地电参数.     

有限长接地导线源电磁场三维正演研究

长接地导线源的一维解析解可以通过对电偶极源解沿导线方向积分获得,而三维模型中则难以直接加载。这里将长接地导线源看作多个电偶极源的组合,利用交错采样有限差分法,进行了直接计算电场总场的有限长接地导线源电磁场三维正演研究,并验证了算法的正确性。正演模拟时为消除场源奇异性,用电流密度呈伪delta函数分布的等效源来代替点源。     

水平层状介质电磁法数据时频联合反演

通过一维层状介质电磁法数据的时域和频域正演结果,对不同电阻率大小的电阻层进行了反演计算,说明了在这种情况下电磁法数据时频联合反演的合理性.正演计算了一维层状模型下水平电偶极频率域水平电场响应及利用G-S变换计算时间域垂直磁场导数响应;利用最光滑模型约束反演方法建立了时频联合反演迭代公式.通过理论模型反演试算表明,该反演效果较好.     

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

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

有耗媒质中脉冲响应的研究

本文详细论述了一维有耗媒质中脉冲响应的正演理论及其模拟技术。深入研究了导电率、介电常数对脉冲响应的影响。以平面脉冲电磁波在一维介质中的传播理论为基础，编制脉冲响应的仿真模拟软件。最后，用实际模型验证了该模拟技术的正确性和可靠性。     

差分法在电阻率测深数据处理中的应用

地球物理的测深方法广泛使用于水电工程与固体矿产地质勘探工作中。应用差分法原理进行数据处理，信息提取是提高电阻率测深的解释推断，它具有去掉背景场，突出局部场的作用；cm测深断面正演概念形象鲜明，能提供各种地质体的电场信息．从而有利于从场的形态、场的结构表研究异常性质。文中探索了二阶、三阶差分法的基本原理．并应用于水电工程坝址勘奋．获得较好的效果。     

LANCZOS迭代算法及其在三维地电场数值模拟计算中的应用研究

在三维电阻率的正演计算中往往涉及到快速、准确求解大型线性方程纽Ax=b的问题。通过采用有限差分法来构造出求解点电源三维地电场的大型稀疏对称线性方程组。并引入Lanczos迭代技术，构造出三对角阵方程组，然后采用正交分解法进行求解，它是Krylov子空间方法中的一种。与传统迭代算法相比，它占用内存少，收敛速度快且稳定。针对大型稀疏矩阵及MATLAB语言的特点，采用简单记录矩阵的非零元素值及其所在行、列值的方法存储大型稀疏矩阵，可大大节省机器内存，提高运算速度。理论分析和计算实例显示，此算法是地电三维正演计算的有效方法，为下一步的反演计算打好基础。     

海洋电磁法勘探油气的三维正演模拟

针对海底可控源电磁法( MCSEM) 在油气资源勘探中的实用性,应用矢量有限元进行了多种频率的水平电偶源发射电磁场识别油气层的实例研究,并利用矢量有限元对海底的3D 模型进行了正演数值模拟。利用MCSEM 能较好地识别油气层; 归一化后的电场曲线特征表明,浅海下不易测得油气层,在实验频率内,电场源频率越高对油气层识别能力越强。同时证实了MCSEM 识别油气层是建立在波导理论基础上的。     

地球物理场正演计算的广义差分法

地球物理场正演计算的广义差分法姜佩仁，王新民（长春地质学院综合信息矿产预测研究所长春１３００２６）在地球物理场（或其他类似问题）的正演计算中，目前常用有限差分法或有限单元法。有限差分法计算简单，计算工作量小，但逼近精度比较低；而有限单元法虽然逼近精度...     

探地雷达时域有限差分法正演模拟

在以往雷达波的正演模拟中,借鉴地震波中成熟的正演模拟方法,均采用模拟声波方程的方法,精度不够高。因此从麦克思韦方程组出发,采用目前电磁场理论中最常用的时域有限差分法,对探地雷达进行了正演模拟。为了进一步说明其正确性,对正演模拟的结果又进行了偏移处理,从结果看出该方法的正确性及可行性。     

考虑IP效应的可控源电磁法电磁响应畸变分析

在野外观测的频率域电磁法的电磁场响应,不仅包含电磁感应响应也含有激发极化(IP)响应。为对IP效应引起的电磁场响应畸变进行分析,首先在纯电磁感应有限长导线源E_x、H_y、H_z电磁场响应表达式的基础上引入Cole-Cole模型,其次设置不同的地电模型,分别对含IP效应、纯电磁感应的有限长导线源频率域电磁法进行一维正演。通过理论模拟结果表明,加入IP效应后电磁场响应会产生畸变,随着极化层埋深的增加,响应畸变出现在更低的频率处,同时受IP效应的影响减弱。根据计算的最大相对异常值可知,水平电场E_x对极化层最为灵敏,而水平磁场H_y对深部极化层的反映能力最弱。     

时间谱电阻率法的一维快速正演算法研究

基于偶极－偶极电极装置的瞬变响应表达式,建立了双极源电场的等效偶极和三次插值的快速近似看法。通过在一些典型地电新面上对算法的试算,表明三次插值法是一种较稳定的快速算法,从而实现了时间谱电阻率法一维正演的快速计算。     

巷道层测深理论曲线数值模拟及资料解释方法

建立了层测深理论曲线正常规模拟的数学模型，结合三层和四层地电模型层测深视电阻及视电阻率理论曲线的正演计算，研究了层测深曲线特征和解释方法，同时介绍了煤矿井下层测深法的应用实例。     

海水电阻率分层对海洋CSEM电磁响应的影响研究

海洋频率域可控源电磁法作为海底石油勘探一项新的技术手段,在深海油气勘探中已经取得许多成功的实例。但是,目前在海底地层电磁正演、反演技术中,很少考虑上部海水的电阻率分层对接收数据的影响。这里利用海洋电磁一维正演程序,研究了海水分层情况对理论数据的影响。通过比较分析,可以发现考虑海水电阻率分层对海洋电磁响应有较强的影响,因此,在实际海洋正反演过程中,海水电阻率层应该与地下地层一起参与反演解释,才能提高整体反演解释精度。     

二维随机介质模型正演指导复杂地层地震解释——以薄层砂泥岩交互地层为例

大港油田砂泥岩交互地层存在着非均质复杂地质形态,为保证其精细解释效果,根据随机过程理论,提出了选用二维指数型椭圆自相关函数为复杂非均质地质体进行简化建模的方法。正演算法使用交错网格有限差分法模拟弹性波在二维随机介质模型中的传播,激发方式使用自激自收模式。正演结果表明,随机介质模型模拟大港油田薄层砂岩互层地层具有合理性,交错网格有限差分算法对于该模型的正演数值模拟,算法合适。     

矿井电法三维有限元正演模拟--直流电透视方法技术研究

通过工作面两巷道间直流电透视三维正演模拟、资料解释方法以及正演实例分析，表明不同装置形式的直流电透视技术仍属于体积勘探。但不同的测量方式对应于不同的电场分布特性，且能分别突出探测目的层(煤层内或顶底板岩层内)的异常特征。为解决采煤工作面附近的具体地质问题，使用相应的方法可提供理论依据。     

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

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

基于姿态参数的海洋可控源电磁数据一维反演

基于姿态信息的海洋可控源数据一维反演采用Occam方法进行,反演数据为不同收发距离观测到的电场振幅(MVO)。接收器和发射器姿态参数同时引入到反演中,正演过程引入发射源姿态参数,借助矢量分解将任意姿态的发射源分解为3个正交坐标轴分量源,分别计算接收点处电场响应,叠加得到总电场响应;根据接收器姿态信息,将计算坐标系中的电场矢量进行欧拉旋转获取任意姿态时测量到的电场响应。理论模型与实际数据反演结果表明,基于姿态信息进行海洋可控源勘探数据反演是可行的。     

海洋可控源电磁法数值滤波解及算例

讨论了海洋可控源电磁法水平电偶极子源发射时电磁场的正演计算问题,采用Guptasarma线性数值滤波公式求解。对比验证了两组Guptasarma滤波系数的计算精度,分析了电场表达式核函数的变化规律。模型算例表明,高阻薄层的存在能使电场响应增强,若要提高海洋可控源电磁法对浅部高阻薄层的分辨能力,应增大发射电流的频率并采用短偏移距测量。振幅随偏移距变化曲线同样表明该方法对高阻水平薄层有较高的分辨率,适用于高阻油气藏的识别,这些都说明开展一维正演问题的研究,可以了解地质体的电磁响应特征及变化规律,为野外资料的反演解释奠定基础。     

垂直接触带上电阻率异常的空间分布特征

针对垂直层状结构岩层的电阻率分布特征问题,采用直立多层介质地面电场递推公式,计算岩层为多个垂直接触面时的地面视电阻率。利用Hankel线性滤波对零至无穷大区间的积分进行数值计算,使用多种电阻率测量装置进行正演计算,得出垂直接触带上多种测量装置的理论视电阻率异常特征。