共查询到20条相似文献，搜索用时 156 毫秒
1.
在吸收并修正部分前人工作的基础上，给出了电导率分块均匀的瞬变电磁25维有限单元法正演算法. 在该算法中，（1）采用在矩形网格中进一步三角细化的剖分方式，同时考虑到矩形单元中心节点的性质，求解线性方程组前用高斯消元法消除矩形网格中心节点上的待求未知量，这样一方面较准确地模拟出任意二维复杂地电断面，另一方面又大大地节省了计算量；（2） 采用计算二次场的算法，事先仅需获得电场和磁场沿地电构造走向方向两个分量的一次场，故而不但计算精度较高，而且不明显增加计算量. 最后，对几例典型模型进行正演，并与解析解对比，验证了本文方法的有效性. 相似文献
2.
瞬变电磁法的野外资料解释基本上停留在一维水平，本文从生产实践中常用的中心回线法入手，基于电磁场本身的叠加原理，从麦克斯韦方程组出发，导出了中心回线瞬变电磁2.5维二次场（纯异常）的有限单元计算公式. 该算法采用三角形有限元网格，在尽可能拟合地下电性断面的情况下减少有限元网格的节点数和单元数；用选主元的LU分解法求解线性方程组，做到了在移动场源时只需改变右端项；在反傅氏变换中，使用新的波数选取方案，让波数m随时间t滑动. 最后给出的算例表明，该算法直接计算异常场，计算速度快、精度高. 相似文献
3.
假定垂直轴为二维对称各向异性介质主轴之一，构造走向与另一主轴方向成任一夹角，用伽勒金（Ｇａｌｅｒｋｉｎ）有限元法和矩形网格，优化地合成总体刚度矩阵，使二维各向异性介质的基本方程形成有限元代数方程组，求出各节点场值，并利用ＭＯＭ法求出辅助场，进而求出张量阻抗等响应函数．同时还对有关文献的计算模型进行了数值模拟检验 相似文献
4.
为了易于模拟野外复杂地形和地下任意形状地电体模型,将有限元单元网格设计为三角单元;并考虑到野外实际勘探中,地球介质的电性参数均是连续变化的情况,单元内的场值和电性参数被设计为双线性变化;推导出二维起伏地形条件下大地电磁法有限元数值模拟算法;根据单元节点主场值和线性插值形函数间的关系,计算出单元节点的辅助场值;在二维起伏地形情况下,定义TE、TM模式视电阻率和阻抗相位.4个模型的计算的结果与解析法的均方根误差小于1%,地形模拟与前人的计算结果相符,模拟倾斜界面异常体,能有效的反映出其异常形态. 相似文献
5.
《地球物理学进展》2017,(5)
二维波动方程谱元法通常采用四边形网格,网格质量和属性建模方式对波场模拟精度和效率有较大影响.常规属性建模需要先建立几何模型再对其进行贴体网格剖分.但存在已知稠密属性控制点却难以建立复杂几何模型的问题.针对此问题,提出对模型进行矩形网格剖分并计算GLL点,利用属性控制点对所有GLL点进行插值以提高属性建模的效率.将属性建模方式归纳为单元属性建模方式和节点属性建模方式.在节点属性建模方式中分析了双线性插值、快速Gauss径向基函数插值法的计算效率,将两种插值算法集成于SPECFEM2D的属性建模程序中.实例表明两种属性建模方式对谱元法波场模拟都是有效的. 相似文献
6.
本文基于非结构网格实现了海洋可控源电磁法三维有限元正演模拟.该算法采用完全非结构网格剖分,可以模拟任意起伏地形和复杂地电模型.为了避免场源的奇异性,采用一次场/二次场分解算法,一次场由基于Schelkunoff势函数的一维解析公式得到.为了提高算法的精度和效率,采用对测点附近单元和异常体区域进行体积约束加密的方法,实现了非结构网格的局部加密.一、二维模型计算和分析表明,本文采用的局部加密方法能够明显地改善算法的精度,最大相对误差基本在1%以内.对三维模型计算及对比分析,说明了该算法对三维可控源电磁正演的实用性.复杂海底地形模型的正演模拟表明,海底地形对电磁场的影响很大,在进行海洋可控源电磁资料解释时,地形的影响有必要考虑在内. 相似文献
7.
直流电阻率数值模拟对于实测电阻率数据反演、解释地下真实电阻率结构具有重要意义.为解决复杂地电模型剖分以及处理截断边界时面临的困难，本文提出一种新的三维自然单元-无限元耦合法.该算法在自然单元区域灵活布设自然节点以填充复杂模型，并引入无限单元替代截断边界条件以促使刚度矩阵与电源位置无关.首先推导了稳定电流场基本方程；然后介绍了二、三维自然单元法及其形函数的构建过程；再引入三维多向映射无限单元并实现自然单元-无限元的有效耦合；最后通过4个数值模型，结合传统有限单元法、有限元-无限元耦合法以及COMSOL Multiphysics软件的计算结果，验证了算法的正确性及有效性.数值结果表明，本文所提算法满足精度要求，能有效减小计算区域，无需二次更新刚度矩阵，对复杂模型适应性强. 相似文献
8.
9.
为了计算带任意地形的各向异性介质中二维大地电磁响应，本文在非结构化网格的基础上，采用有限体积法，开发了二维大地电磁各向异性正演模拟的新算法.首先，从Maxwell方程出发，推导二维各向异性介质中大地电磁场的边值问题；然后，采用三角网格自动生成技术对求解区域进行非结构化网格剖分，进而构建节点中心控制体积单元，利用有限体积方法，得到求解边值问题的大型稀疏线性方程组；最后，利用Pardiso精确地计算了大地电磁响应值.三个各向异性模型的计算结果表明，本文开发的有限体积算法，不仅能够高精度求解带任意地形的大地电磁电导率各向异性问题，而且对于同一模型，该方法的计算消耗和精度都与有限单元法相当.因此，有限体积法是处理电磁法各向异性问题的一种有效方法. 相似文献
10.
为进一步提高大地电磁非线性反演的稳定性、运算效率及准确度，将遗传神经网络算法引入大地电磁反演.首先针对大地电磁二维地电模型建立BP（Back Propagation）神经网络基本框架进行学习训练，网络输入为已知地电模型的视电阻率参数，输出为该地电模型参数；再利用遗传算法对神经网络学习训练过程进行优化，计算出多种地电模型网络连接权值和阈值的最优解；最后将最优连接权值和阈值对未知模型进行反演测试，网络输入为未知地电模型的视电阻率参数，输出为该地电模型参数.模型实验表明：遗传神经网络算法充分结合了遗传算法的全局寻优性和神经网络的局部寻优性，相比单一神经网络算法，在网络学习训练中提高了解的收敛成功率和计算速度，在反演测试中能更准确地逼近真实模型.将遗传神经网络算法与最小二乘正则化反演进行对比，理论模型和实测数据都验证了遗传神经网络算法在大地电磁反演中的可行性和有效性. 相似文献
11.
大地电性结构的横向变化会对磁暴时的感应地电流和地面电磁场产生影响.本文假设扰动地磁场变化的源为地面以上一定高度的面电流,以某一典型层状大地电导率结构为基础,构造含有电导率横向突变的地电模型.针对感应电流的方向与横向分界面平行的情形,采用伽辽金有限元法对电导率横向突变处的感应地电场进行了分析,揭示了电导率横向差异产生的趋肤效应和邻近效应的机理,针对与电性结构分界面平行的输电线路,从评估地磁感应电流的角度讨论了影响的严重程度和范围. 相似文献
12.
通过甘肃省嘉峪关台地磁场观测资料,研究嘉峪关台、瓜州台磁静日地电场日变化的时频特征波;由地电场分钟值观测数据的时序叠加残差方法,研究嘉峪关、瓜州山的地电暴变化。结果表明:(1)两台地电场静日变化以两次起伏变化为主,无相位差,但两台之间日变幅差异较大;(2)地电场分量变化与地磁场正交分量变化显著相关;地电场与地磁场日变波形不同,极值时间有差异。2个台存在很明显的高频成分,在去除了高频变化后,其优势周期也相同,从大到小依次为12 h、8 h、24 h。地磁场H分量因存在磁暴影响,故高频变化较多,在去除了磁暴影响后,其优势周期从大到小依次为24 h、12 h、8 h;(3)当电磁暴扰动剧烈时,两台可以较清晰地记录到地电暴的完整变化。在发生电磁暴时,地电场与地磁场的相关性明显降低,且不同台、不同测向之间的变化幅度也不尽相同。两台东分量E_Y暴日的日变幅较静日明显增大,磁暴期间Y分量变化率与地电场东分量E_Y观测数据显著相关,由此说明:两台日变幅的不同与台站台址电导率有关,太阳风引起的电离层活动是引起了地电场日变化主因。引起电暴的原因可能不同于引起日变化的原因,主要是两台之间及不同测向之间的浅、深层电阻率和地质构造等诸多因素的结果。 相似文献
13.
14.
本文应用中国大陆地区11个地电场台站5.5年（2012-01-01—2017-05-31）中磁静日和磁扰日的观测数据，并联系地电场变化敏感响应地下介质非均匀性和空间电磁活动的特点，发展了地电扰动指数，即消除静日地电变化的幅值比地电指数（GEI）.对比各台站的地电指数与地磁指数，87.43%的单台地电指数D与单台地磁指数K在数值上完全相同或相差1，93.75%的单台平均地电指数Ds与单台平均地磁指数[K+0.5]（四舍五入）相同或仅相差1，说明地电指数能用于量化界定地电扰动水平，如同地磁指数量化界定磁扰水平.地电指数将服务于地球电磁环境监测、科学研究等领域，同时也可以用来评估地电场观测数据的质量. 相似文献
15.
16.
认识地震电信号特征对地震短临预报具有重要意义.通过陕甘宁晋地区、川滇地区、山东及其周边地区和上海地区的换流站接地极向地下注入2100—4600 A大电流，研究各地区地电场的变化特征，以深入地认识地电场信号的区域性变化实质.研究显示:① 在陕甘宁晋地区，高沙窝大电流源的附加地电场可被分布于鄂尔多斯地块边缘弧形地带的490 km范围内的10个台站观测到，而距该信号源以西150 km的台站均未观测到该信号；② 在川滇地区，丽江信号源的南部地区观测到的信号强度和距离远大于北部地区.对于西昌大电流源，只有距其10 km处的小庙台和540 km处的仙女山台可观测到信号，而其它十几个台站均未观测到；③ 在山东和上海地区，可被观测到的人工源信号的最远距离及其特征各不相同，亦与陕甘宁晋和川滇地区大为不同；④ 对于本文所研究的地区，信号的波形变化存在着随台站与信号源距离的增加畸变程度增大的现象，一些台站只能记录到大电流开始和结束时的尖峰脉冲信息.针对上述变化特征，本文从机理和地质条件等方面进行了一些有益的探讨，以期加深对地震电信号变化特征的认识. 相似文献
17.
18.
使用地电扰动幅度与静日变化幅值比的地电活动指数(GEI)计算方法,连续计算江苏、安徽两省6个地电场台站2018年第一季度的单台地电指数D、平均地电指数Ds。同时将地电指数D、Ds与单台地磁指数K、平均地磁指数K(即[K+0.5])对比得到:86.48%的指数D与指数K完全相同或仅相差1,89.44%的指数D s与指数完全相同或仅相差1。两类电、磁指数的高度吻合说明:(1)连续产出的地电指数是可靠的。(2)地电指数D、Ds可以量化界定地电场活动水平。极个别的地电指数D、Ds与地磁指数K、仍存在一定差异,这与地电场台站观测有关,提升地电场观测资料的质量能更好地确定地电指数。(3)地电指数与地磁指数相互弥补能更准确地确定地球电磁活动水平,其在相关的地球科学技术领域有广泛的应用前景。 相似文献
19.
地空电磁法已经成为深部资源勘探的重要地球物理方法,但对频率域地空系统的三维多源电磁响应特征研究较少.本文设计了多种激励源组合方式,采用非结构化有限元数值模拟方法,对三维地电模型的空中垂直磁场的响应特征进行了研究.首先推导了基于电场的双旋度公式及其变分形式,加入罚项以减少伪解的影响.接着把有限元稀疏矩阵方程转换为频率的函数,采用Krylov子空间投影方法,通过模型降阶算法降低稀疏矩阵的阶数,实现多频点的快速计算.建立了三维低阻体模型、高阻体模型以及两个相邻低阻体模型,分别采用单源、双源、三源和四源激励模式,从垂直磁场的总场、二次场响应和全域视电阻率等方面进行分析比较.结果表明:多源地空电磁法不仅可以增加总场的强度,而且可以改变异常体的二次电磁响应分布规律.各电偶源延长线呈正三角形分布的三源和矩形分布的四源激励模式在增强信号强度以及削弱异常体的边界效应方面具有一定的优势,是一种优化的多源激励方式. 相似文献
20.
在地震电磁学研究中,电流在地下流动特征与构造关系是一个非常重要的问题.多个强电流源的存在和大范围地电场观测站的分布格局为我们研究这个问题提供了更好的机会.本文通过利用上海及其周边地区4个地点不同的高压换流站接地极向地下注入的大电流(1200~4780A)信号,研究了华东地区8个台站组成的地电场观测台网接收到的附加地电场信号,推进了对地电流在地下传播特征的认识.研究显示:1从幅度特征而言,对于源自不同的大电流信号,各台站观测到的附加地电场信号幅度特征有很大不同,而近处的浦东台只能观测到南桥大电流源的信号而观测不到其他3个大电流源发出的信号,存在着所谓的"敏感点"效应;2对于上海地区4个接地极大电流而言,能够观测到源自它们的附加地电场的台站最远为350km,在更远的台站则观测不到.3从均匀度方面而言,在每个台站对于不同的信号源,记录到的附加地电场信号在各方向上的长极距和短极距信号幅度之比可以有很大变化;4在方位特征上,只有崇明台站接收到的同里接地极和华新接地极大电流信号的计算方位度数与实际方位度数误差最小,分别为0.2°和0.8°.而在距离各个信号源最近的青浦台,这种误差要比崇明台站的大;5极化方向特征显示,在青浦台和崇明台所观测到的附加地电场信号中,存在着同一方向上长、短极距信号极化方向是反向的现象,这与常理相悖,怎样解释该现象有待于今后进一步研究.本文从点电流源与台站测线布局构成的装置系统、大区域介质非均匀性、台站测区以及附近地区的介质细结构的影响三个方面,对上述特征进行了较深入的分析和解释. 相似文献