共查询到16条相似文献，搜索用时 125 毫秒
1.
本文发展了积分方程法用于层状介质中三维不均匀体的大地电磁模拟算法（简称MT）.分为二个步骤:第一步,异常体用等效的散射电流代替,通过层状介质中的格林函数,建立以散射电流为未知参数的积分方程;第二步,把求得的散射电流乘上相应的格林函数,即得地面上的二次电磁场,由此而计算出各种MT响应. 文中采用了数值滤波与插值、群变换以及格林矩阵带状化三个方面的数值处理方法,提高了计算效率.通过与已发表的三维MT计算结果对比及格林函数互易性检验,表明了该算法的正确性.在此基础上,进行了数值模拟,初步讨论了三维MT曲线的畸变特点. 相似文献
2.
3.
4.
利用积分方程方法以及半空间并矢格林函数的快速算法对埋地目标体矢量电磁散射进行正演计算。首先,利用半空间电磁并矢格林函数建立起埋地目标体的体积分方程。然后通过将空间偏导转移至格林函数谱域积分的积分号之外,并采用离散复镜像方法来近似余下的零阶索末菲积分,进一步得到并矢格林函数各个分量的闭合形式。由于避免了对索末菲积分的繁琐数值计算,使得生成反应矩阵和计算散射场时由半空间并矢格林函数计算带来的瓶颈问题得到较好的克服,因而极大地提高了埋地目标体电磁响应正演计算效率,同时也能保证足够的精度。 相似文献
5.
系统地论述了大地电磁三维交错采样有限差分数值模拟算法实现过程中交错网格剖分、积分公式离散化、边界条件、方程组求解、三维张量阻抗的计算等内容. 由于提出了简洁的边界条件，采用了解大型系数矩阵方程组的双共轭梯度稳定解法，所实现的三维交错采样有限差分数值模拟算法具有迭代收敛稳定、计算精度高、速度快等特点. 通过两个理论模型的计算结果检验了算法的正确性和计算精度. 所实现的三维交错采样有限差分数值模拟算法为研究三维反演问题奠定了基础. 相似文献
6.
为了实现大地电磁测深三维高效快速正演算法,本文研究了大地电磁测深积分方程三维正演模拟的理论与方法.首先详细地推导了Fredholm型积分方程;然后给出积分方程的离散方法、并矢格林函数积分的求解方法以及大型矩阵方程的存储与求解方法;最后,通过理论地电模型的计算检验了方法和程序编制的正确性,并展示了几种理论地电模型的数值计算结果,分析其异常特征.算例结果表明,该正演算法在CPU主频为1.81GHz、内存为1GB的计算机上能正确地模拟多个异常体的异常特征;地表电性不均匀体会使地表视电阻率曲线产生十分严重的畸变;TE模式相对TM模式纵向分辨率比较好,而TM模式的横向分辨率较好;随着异常体埋深的增加,反映异常的分辨能力进一步降低. 相似文献
7.
8.
当地表存在三维非均匀电导率分布时,区域大地电磁响应发生畸变. 以往对这种畸变研究多假设近地表为三维,区域构造为一维或二维. 对于更一般的三维/三维构造,为了分析并消除这种畸变影响,真实反映地下三维区域构造信息,本文实现了三维大地电磁相位张量积分方程数值算法,并研究在不同地质模型下相位张量响应. 结果表明,相位张量不仅可以反映一般三维构造信息,亦可有效反映复杂近地表构造下三维区域构造信息,而无须假设区域构造为一维或二维,证明相位张量具有较强抗近地表局部非均匀构造干扰能力,能够保持更为一般的三维区域构造信息. 为了加快正演计算,同时保持一定精度,算法采用了积分方程多网格法. 相似文献
9.
10.
瞬变电磁法正演计算进展 总被引：4，自引：1，他引：3
详细介绍了瞬变电磁法正演计算的方法、现状和发展趋势.瞬变电磁法一维正演计算需要将电磁场从频率域转换至时间域,转换方法有三种,分别是Gaver-Stehfest算法、余弦变换和Guptasarma算法.在这三种方法中,使用较多的是Gaver-Stehfest算法和余弦变换,Gaver-Stehfest算法速度较快,但精度不及余弦变换.瞬变电磁法的数值模拟主要集中于2.5维和三维,使用的数值计算方法有积分方程法、有限差分法、有限单元法和SLDM法.积分方程法主要在三维数值模拟中使用,现已很少使用;有限差分法和有限单元法是目前瞬变电磁法2.5维和三维数值模拟的主要方法;SLDM法主要应用于三维数值模拟.我国瞬变电磁法正演计算成果主要集中在回线源激发的瞬变电磁场一维数值计算和利用有限单元法进行2.5维和三维数值模拟.瞬变电磁法正演计算的发展趋势有:数值算法的改进、提高计算效率和研究地形对瞬变电磁场的影响规律. 相似文献
11.
12.
本文给出一种用积分方程法配合有限差分法计算位于导电介质中三维导体的电磁响应的新途径,在二层大地条件下给出了具体计算方法和结果。详细讨论了围岩介质的导电性对导体异常的影响。这对于开展低阻覆盖层地区的电磁法工作具有一定的指导意义。 相似文献
13.
Electromagnetic induction of a three-dimensional conductivity inhomogeneity in the two-layered earth
Summary The calculation of the elements of Green's tensor function is presented for solving the problem of the electromagnetic induction by means of a vector integral equation. A two-layered Earth is considered as the medium, the surface layer including a three-dimensional conductivity inhomogeneity. Use is made of the boundary condition requiring the vertical component of the electric current to be zero at the Earth's surface which partly simplifies the theoretical computation. Long-period asymptotics of the individual complicated functions, occurring in Green's tensor function as well as in the tensor function required to calculate the components of the anomalous magnetic field at the surface of the halfspace, were effected. With the aid of these asymptotics one can obtain estimates of the functions occurring in the theoretical analysis of the problem. 相似文献
14.
格林函数在解决电磁场和直流电场问题的积分方程法和边界单元法等方法中有着广泛的应用.在应用格林函数时会遇到格林函数奇异性的问题,即其在源区域的内部,当场点和源点重合时,会碰到对格林函数为无穷大,积分为奇异积分.对于不同问题中遇到的不同形式的格林函数奇异性问题的各种处理方法进行了分析和评述.对于标量格林函数奇异性问题处理方法有：挖去法、级数展开法、绕开法和解析法等.对于并矢格林函数奇异性问题处理方法有：分量处理法、源并矢法和拟源并矢法等处理方法.通过对实际工作中所遇到的新方法新问题中的不同形式的格林函数的研究,提供了奇异性问题的处理的方法和途径. 相似文献
15.
16.
《地震科学（英文版）》2019,(Z1)
A two and a half dimensional(2.5 D) multidomain indirect boundary element method(IBEM) is developed to study the wave scattering of obliquely incident P-, SVand SH-waves by a hill-valley staggered topography in a multi-layered half-space. The IBEM algorithm includes using2.5 D full-space and half-space Green's functions to construct scattered fields in decomposed closed and opened half-space regions, respectively, and using the dynamic stiffness method to solve the free fields. All regions are finally integrated by introducing the compatibility conditions to obtain the total wave fields. The proposed 2.5 D IBEM has the flexibility in dealing with complex boundaries by directly applying the fictitious loads on the regions' boundaries, with a less storage requirement compared to the full 3 D models. Besides, by combining the specific advantages of the two kinds of Green's functions, the method is well suitable for handling coupled topographies with high accuracy. The method is validated by comparison with published results for a single valley as well as a single hill topography. The effects of height-to-width ratio of hill and layering on dynamic responses are further parametrically investigated by numerical implementations in frequency domain. Results show that the interaction between valley and hills can lead to a more significant amplification within the valley region, and dynamic responses are deeply influenced by the height-towidth of hill, simultaneously depending on incident angle and frequency. Besides, the site effects become more complex when the stratification feature is taken into account. 相似文献