共查询到18条相似文献，搜索用时 187 毫秒
1.
电磁场数值模拟的背景场/异常场算法是三维正演的有效策略之一,优点为采用解析法计算电磁场背景场代替场源项、克服了场源奇异性,缺点为不适用于发射源布置于起伏地表或背景模型复杂的情形.总场算法是直接对电磁场总场开展数值模拟,其难点是有效加载场源、保证近区与过渡区数值解精度.本文以水平电偶源形式分段加载接地长导线源,并以电场总场Helmholtz方程为矢量有限元法控制方程,实现了基于非结构化四面体网格剖分的接地长导线源频率域电磁法三维正演.通过与均匀全空间中水平电偶源产生的电场解析解对比,验证了本文算法的正确性,并分析了四面体外接圆半径与其最短棱边的最大比值和四面体二面角最小值对数值解精度的影响规律.通过与块状高导体地电模型的积分方程法、有限体积法和基于磁矢量势Helmholtz方程的有限元法数值解对比,进一步验证了本文算法正确性,同时说明了非结构化四面体网格能够更加精细地剖分电性异常体,利于获得精确数值解. 相似文献
2.
与可控源音频大地电磁（CSAMT）相比,广域电磁法通过采用全区视电阻率定义,突破了卡尼亚视电阻率所需的远区条件限制,极大拓展了可控源电磁观测区域和探测深度.考虑到电偶源激发场的三维特征以及地下复杂三维结构,为提高广域电磁数据解释精度,本文实现了基于二次耦合势的广域电磁法三维正演计算.该算法利用Helmholtz定理将麦克斯韦方程转化为库伦规范下的磁矢势和电标势耦合方程,有效改善了离散所得大型线性方程组的谱性质,并通过强加散度条件来消除电场伪解的影响.此外,采用散射场方法,其中一次场使用准解析法求解,二次场使用有限体积法求解,克服了局部激发场源奇异性问题.通过与一维层状模型下电偶源产生的电磁场准解析解对比,验证了本文算法的正确性.在此基础上,利用本文的正演算法对比分析了广域电磁法与CSAMT对典型三维目标体的探测能力,结果表明在相同的观测条件下,广域电磁法能够更准确地反映地下目标体信息,拥有更优的分辨能力. 相似文献
3.
《地球物理学进展》2015,(6)
为了应用时间域电磁法进行深部勘探,需要计算大尺度模型的晚时时域响应.本文将待求解的电磁响应分解为一次场和二次场之和,实现了波形为拟高斯脉冲的大功率脉冲源激励的一次场的计算.对于三维异常体产生的二次场,采用基于非均匀步长交错网格的时域有限差分(FDTD)和非等时步长的迭代算法求解,在层状介质模型的上阶跃响应计算结果对比以及与积分方程计算的瞬断响应结果对比验证基础上,实现了大功率脉冲源激励下晚时时域电磁响应的计算.通过设计一个简单异常体的三维模型,采用电偶极源脉冲电流激发,计算的时间域响应很好的揭示了大功率脉冲源激励的场在地中随时间扩散,以及三维异常体产生的异常场并二次扩散的过程. 相似文献
4.
从电偶源三维地电断面可控源电磁法的二次电场边值问题及其变分问题出发,采用任意六面体单元对研究区域进行剖分,并且在单元分析中同时对电导率及二次电场进行三线性插值,实现电导率分块连续变化情况下,基于二次场的可控源电磁三维有限元数值模拟.这个新的可控源电磁三维正演方法可以模拟实际勘探中地下任意形状及电性参数连续变化的复杂模型.理论模型的计算结果表明,均匀大地计算的视电阻率误差和相位误差分别为0.002%和0.0005°.分层连续变化模型的有限元计算结果表明,其与对应的分层均匀模型解析结果有明显差异.三维异常体组合模型以及倾斜异常体等复杂模型的有限元计算结果也有效地反映了异常形态. 相似文献
5.
在一些地层层理发育的地区,地下介质存在显著的电各向异性,此时基于各向同性模型解释含各向异性效应的可控源音频大地电磁（CSAMT）测深观测数据会导致错误的结果.本文通过引入3×3的对称正定张量表征电导率各向异性,采用非结构四面体网格和矢量有限元方法离散电场满足的矢量Helmholtz方程,并将电磁场源等效为系列电偶极子,实现任意各向异性介质中CSAMT高效数值模拟.本文首先通过层状各向异性模型检验三维有限元算法的精度和有效性,进一步建立三维地电模型研究异常体各向异性和围岩各向异性对CSAMT响应的影响,最后使用视电阻率极性图来识别各向异性电导率主轴方向.数值模拟结果表明,各向异性电导率对CSAMT视电阻率幅值及分布规律都有很大影响,视电阻率极性图能够很好地识别各向异性主轴方向. 相似文献
6.
地空电磁法已经成为深部资源勘探的重要地球物理方法,但对频率域地空系统的三维多源电磁响应特征研究较少.本文设计了多种激励源组合方式,采用非结构化有限元数值模拟方法,对三维地电模型的空中垂直磁场的响应特征进行了研究.首先推导了基于电场的双旋度公式及其变分形式,加入罚项以减少伪解的影响.接着把有限元稀疏矩阵方程转换为频率的函数,采用Krylov子空间投影方法,通过模型降阶算法降低稀疏矩阵的阶数,实现多频点的快速计算.建立了三维低阻体模型、高阻体模型以及两个相邻低阻体模型,分别采用单源、双源、三源和四源激励模式,从垂直磁场的总场、二次场响应和全域视电阻率等方面进行分析比较.结果表明:多源地空电磁法不仅可以增加总场的强度,而且可以改变异常体的二次电磁响应分布规律.各电偶源延长线呈正三角形分布的三源和矩形分布的四源激励模式在增强信号强度以及削弱异常体的边界效应方面具有一定的优势,是一种优化的多源激励方式. 相似文献
7.
本文以非结构化网格剖分为基础,推导了电导率呈任意各向异性条件下的海洋可控源电磁法二次场磁矢量位、标量位所满足的有限元方程.将不完全LU分解预条件因子(ILU)与Induced dimension reduction(IDR(s))迭代算法相结合对有限元方程进行求解.之后,采用加权移动最小二乘法对二次场矢量位、标量位进行求导得到相应的电磁场各分量.最后,开展了两个地电模型的电磁场计算与分析.结果表明:本文所提出的电导率呈任意各向异性条件下的海洋可控源电磁法三维正演算法正确;不完全LU分解预处理的IDR(s)算法计算效率明显优于常用的ILU-QMR、ILU-BICGSTAB算法;算法具有良好的通用性,可用于陆地电磁、航空电磁、井中电磁等地球物理勘探方法在电导率任意各向异性方面的研究. 相似文献
8.
可控源电磁法具有分辨率高及抗干扰能力强等特点,是一种重要的地电磁勘探方法.目前,可控源电磁法的高精度正演计算一直是其核心研究问题之一.传统积分方程法一般采用近似积分公式、简单矩形网格和近似的奇异性体积分计算技术,制约了体积分方程法处理复杂地下异常体的能力,降低了计算精度.针对上述问题,本文基于完全积分公式、四面体非结构化网格和奇异体积分的精确解析解来高精度求解复杂可控源电磁模型的正演响应.首先,从电场积分公式出发,推导了可控源电磁问题满足的积分方程;其次,借助于非结构化四面体网格离散技术,实现了地下复杂异常体的有效模拟.最后,利用散度定理把强奇异值体积分转换为一系列弱奇异性的面积分公式,并通过推导获得了这些弱奇异性的面积分公式的解析解,从而最终实现三维可控源电磁问题的高精度积分求解.以块状低阻体地电模型为测试模型,采用本文提出的积分方程方法获得的数值解与其他公开数值算法解进行对比分析,其对比结果具有高度的吻合性,验证了算法的正确性;同时,设计了球状及复杂地电模型进行算法收敛性测试,进一步验证算法的正确性以及能够处理地下复杂模型的能力. 相似文献
9.
本文基于非结构网格实现了海洋可控源电磁法三维有限元正演模拟.该算法采用完全非结构网格剖分,可以模拟任意起伏地形和复杂地电模型.为了避免场源的奇异性,采用一次场/二次场分解算法,一次场由基于Schelkunoff势函数的一维解析公式得到.为了提高算法的精度和效率,采用对测点附近单元和异常体区域进行体积约束加密的方法,实现了非结构网格的局部加密.一、二维模型计算和分析表明,本文采用的局部加密方法能够明显地改善算法的精度,最大相对误差基本在1%以内.对三维模型计算及对比分析,说明了该算法对三维可控源电磁正演的实用性.复杂海底地形模型的正演模拟表明,海底地形对电磁场的影响很大,在进行海洋可控源电磁资料解释时,地形的影响有必要考虑在内. 相似文献
10.
11.
Munirdin Tohti Yibo Wang Evert Slob Yikang Zheng Xu Chang Yi Yao 《Geophysical Prospecting》2020,68(6):1927-1943
Seismoelectric coupling in an electric isotropic and elastic anisotropic medium is developed using a primary–secondary formulation. The anisotropy is of vertical transverse isotropic type and concerns only the poroelastic parameters. Based on our finite difference time domain algorithm, we solve the seismoelectric response to an explosive source. The seismic wavefields are computed as the primary field. The electric field is then obtained as a secondary field by solving the Poisson equation for the electric potential. To test our numerical algorithm, we compared our seismoelectric numerical results with analytical results obtained from Pride's equation. The comparison shows that the numerical solution gives a good approximation to the analytical solution. We then simulate the seismoelectric wavefields in different models. Simulated results show that four types of seismic waves are generated in anisotropic poroelastic medium. These are the fast and slow longitudinal waves and two separable transverse waves. All of these seismic waves generate coseismic electric fields in a homogenous anisotropic poroelastic medium. The tortuosity has an effect on the propagation of the slow longitudinal wave. The snapshot of the slow longitudinal wave has an oval shape when the tortuosity is anisotropic, whereas it has a circular shape when the tortuosity is isotropic. In terms of the Thomsen parameters, the radiation anisotropy of the fast longitudinal wave is more sensitive to the value of ε, while the radiation anisotropy of the transverse wave is more sensitive to the value of δ. 相似文献
12.
常规的三维时间域航空电磁模拟通常采用隐式步长方法进行时间离散,需要几次矩阵分解和上百次右端源项回带,计算效率较低.为了提高正演计算效率,本文提出使用有理Krylov方法求解时间域电场扩散方程.首先使用非结构四面体网格进行空间离散,采用Nédélec矢量基函数近似四面体单元内的电场;然后基于有限元离散给出矩阵指数和矢量乘积表示的电场显式解;最后采用有理Arnoldi算法构造Krylov子空间内的正交基函数并进一步求解矩阵指数与矢量的乘积,直接得到任意时刻的电场解向量,避免步长离散过程.此外,本文还提出一种指数加权偏移参数优化方法,使得有理Arnoldi近似在瞬变衰减晚期具备更高的精度,从而降低Krylov子空间阶数并提高计算效率.通过和层状模型解析解的对比验证了有理Krylov方法的精度.针对三维异常体模型使用全局网格和局部网格剖分并和其他数值方法比较,进一步说明了有理Krylov方法的有效性. 相似文献
13.
现有海洋可控源电磁三维数值模拟方法大多基于电导率各向同性介质理论,不能模拟海底地层电导率各向异性的实际情况.本文给出了电导率各向异性三维介质中电性源海洋可控源电磁二次电场的边值问题以及相应的变分问题,采用长方体单元对研究区域剖分,将场分量定义在剖分单元的边上,利用矢量有限单元法求解变分问题,实现了电导率任意各向异性海洋可控源电磁三维矢量有限元数值模拟.这个新的正演方法可以计算电导率任意各向异性三维地电模型的海洋可控源电磁响应,基于二次场矢量有限元法直接求解电磁场,避免了传统有限元方法可能遇到的伪解问题和难于处理电场法向分量不连续的问题,提高了数值模拟计算精度.一维电导率各向异性模型电磁场数值解与解析解吻合得相当好,无论在源附近还是远离源处相对误差均不超过1%.电导率各向异性二维模型的计算结果与已有文献采用的非结构有限元模拟结果十分吻合.三维地电模型数值模拟结果显示,电导率各向异性张量电导率主轴分量和欧拉角对不同装置海洋可控源电磁响应均有着明显的影响. 相似文献
14.
根据库伦规范势的定义,推导出关于磁矢量势和电标量势的偏微分方程,为了克服由电流源引起的奇异性和数值模拟计算困难,将电磁总场分解为一次场和二次场,一次场由基于Schelkunoff势函数的一维正演算法得到,二次场由有限元法计算得到,实现了海洋可控源电磁法三维有限元正演算法。通过一维数值模拟实例,验证该算法的计算精度。然后,利用该算法对带海底地形的三层储层模型进行正演,分析了海底地形对海洋控源电磁场各分量产生的影响。 相似文献
15.
本文提出了一种基于非结构化网格的海洋电磁有限单元正演算法.为了回避场源奇异性,文中选用二次场算法,将背景电阻率设置为水平层状且各向异性,场源在水平层状各向异性介质中所激发的一次场通过汉克尔积分得到.基于Coulomb规范得到二次矢量位和标量位所满足的Maxwell方程组,通过Galerkin加权余量法形成大型稀疏有限元方程,采用不完全LU分解(ILU)预条件因子的quasi-minimum residual(QMR)迭代解法对有限元方程进行求解得到二次矢量位和标量位;进而,利用滑动平均方法得到二次矢量位和标量位在空间的导数,由此得到二次电磁场;通过一维模型对算法的可靠性进行验证,与此同时,针对实际复杂海洋电磁模型,比较有限元模拟结果与积分方程模拟结果,进一步验证算法精度.若干计算结果均表明,文中算法具有良好的通用性,适用于井中电磁、航空电磁,环境地球物理等非均匀且各向异性介质中的电磁感应基础研究. 相似文献
16.
本文将大地电磁场分解为一次场和二次场,应用交错网格有限差分法模拟计算大地电磁二次场,并引入各向异性最佳匹配层(PML)吸收边界条件作为二次场边界条件,实现了耦合PML吸收边界条件的三维大地电磁二次场有限差分正演模拟.为了确保正演的稳定性和效率,QMR求解器和磁感应矢量散度校正技术被用于PML吸收边界条件下系数矩阵的快速求解.三维模型正演响应表明,基于二次场的三维大地电磁有限差分算法具有较高的计算精度和可靠性.通过计算分析不同PML吸收因子条件的大地电磁正演结果,显示在适当的吸收因子下,PML吸收边界条件可较大幅度的减小外边界距离,从而有效的压缩模型求解空间,最终提高三维大地电磁正演模拟的效率. 相似文献
17.
本文针对多辐射场源地空瞬变电磁法理论,建立了相应的多分量全域视电阻率定义方法.分析了利用磁场强度进行多辐射场源地空系统全域视电阻率定义的优点,针对磁场强度的各个分量,提出了各自的全域视电阻率算法,实现了多分量、全时域、全空域视电阻率计算,并分析了偏移距对全域视电阻率的影响.通过调整源的相对位置及电流方向等参数,多辐射源瞬变电磁地空系统不仅可以加强不同分量信号强度,削弱随机干扰,还可更好地分辨地下异常体的位置.通过对所设计模型的处理,证实了多辐射场源地空系统多分量全域视电阻率算法的有效性,也验证了多辐射场源地空系统的优势. 相似文献
18.
传统三维大地电磁各向异性模拟均是基于规则六面体网格,计算精度有限且较难拟合复杂地质条件.本文采用面向目标自适应非结构矢量有限元法,对三维大地电磁各向异性介质进行模拟.首先从电场双旋度方程出发,利用伽辽金方法建立变分方程;然后利用电流密度连续性条件构建适合大地电磁各向异性问题的加权后验误差估计方法,实现面向目标的网格自适应正演;最后通过典型算例分析各向异性对网格自适应和大地电磁响应的影响特征以及各向异性的识别方法.本文算法能够高精度地拟合起伏地表和任意各向异性介质,适用于分析复杂地电条件大地电磁响应特征,为提高大地电磁资料解释水平提供了理论基础. 相似文献