共查询到20条相似文献，搜索用时 15 毫秒
1.
2.
3.
4.
5.
利用有限元方法进行大地电磁正演数值模拟时,由于是在有限网格区域上的数值计算,模拟计算时的网格边界为截断边界,而有限元数值模拟时的大地电磁场边界条件需要在足够远处才能够满足,所以截断边界的存在可能会使大地电磁正演模拟的边界条件无法满足,致使对计算结果和计算精度产生影响。利用有限元二维正演程序,在网格边界处加载一维情况下的大地电磁场,然后固定研究区域的网格剖分,并对一维地电模型和二维地电模型在改变有限元网格边界大小的情况下进行计算。在对一维模型进行模拟计算时,截断边界对边界条件没有影响,边界条件自然满足。而对二维模型进行模拟计算时,截断边界的存在对计算结果有较大影响。利用趋肤深度作为有限元网格边界变化的量度,通过改变网格边界大小,对不同的二维地电模型进行计算比较,总结出适合大地电磁有限元正演模拟的参考网格边界。 相似文献
6.
利用快速标量传递算法进行了一维层状模型的地震面波频散曲线正演计算,得到了基阶频散曲线;利用传统方法对一维层状模型进行了大地电磁正演计算,得到大地电磁的视电阻率曲线数据和相位数据。利用遗传算法对合成数据分别进行了单独反演和联合反演,得到拟合差分布图来分析两种数据的兼容性并研究联合反演的可行性。通过对比单独反演和联合反演的结果,可以看出地震面波频散曲线和大地电磁数据的联合反演结果较之单方法的反演结果要更加精确,收敛速度更快,两种数据的耦合情况较好,同时联合反演能够有效地限制多解性。 相似文献
7.
音频大地电磁法是深埋金属矿勘探中的一种常用物探方法,研究音频大地电磁对深埋、薄层锰矿的勘探效果具有重要的理论意义及实践价值。论文首先通过建立典型地电模型,完成了大地电磁正演模拟、反演计算及大地电磁的磁反射系数成像,从理论上提出磁反射系数可以提高大地电磁法的深部分辨率;然后以贵州松桃杨立掌锰矿床毗邻区的次级成锰槽盆为勘探对象,完成了横跨成锰槽盆的音频大地电磁勘探及资料处理,获得了沿测线方向的二维电阻率模型和磁反射系数成像等值线图;最后参考研究区物性资料及钻井柱状图,进行了电阻率模型与磁反射系数成像结果的地质地球物理综合解释。论文的研究成果说明,磁反射系数能够有效提高大地电磁法对深部电性层的划分能力,合理运用磁反射系数开展大地电磁多参数解释有利于提高大地电磁法的解释分辨率。 相似文献
8.
《物探化探计算技术》2019,(5)
有限元求解大地电磁正演问题时,对研究区域的剖分常规做法是采用规则化的网格。但规则化网格在剖分地形、断层、褶皱等复杂模型会产生较大的几何离散误差。针对上述情况,这里采用非结构化的四边形网格对二维地电模型进行网格剖分,并与自适应有限元相结合,由剖分的粗网格出发,利用每个频点下网格单元的后验误差估计值指导网格的局部加密,优化网格质量和数量,从而提高正演模拟的精度。通过一维K型地电模型利用本文方法算出数值解与解析解进行对比分析,验证了自适应有限元法在求解大地电磁二维正演模拟中的有效性。并通过对断层、褶皱模型的正演模拟,分析了其大地电磁正演响应特征。 相似文献
9.
薄层识别一直是大地电磁测深法(MT)数据反演解释的重点和难点.约束反演能有效提高大地电磁测深对薄层的识别能力,这里采用在井资料约束上覆地层与围岩的电性条件下,进行MT约束反演.通过约束浅部地层,提高对深部地层的反演识别能力,达到对目标层的准确识别与精确探测.通过模型设计与计算,讨论了地电参数约束类型、约束数量、以及测量... 相似文献
10.
利用一维正则化反演进行大地电磁测深数据拟二维反演解释 总被引：1，自引：0，他引：1
大地电磁测深的反演问题是不适定的,其反演结果不稳定,且具有非唯一性。通过在目标函数中采用正则化方法,可以使得不适定反演问题具有稳定的反演结果,并改善解的稳定性和非唯一性问题。为了提高野外大地电磁测深数据的处理效率和初步解释的精度,提出了大地电磁测深数据的一维正则化反演进行拟二维反演解释方法。这里所述的大地电磁测深一维反演解释,与以往的解释方法不同,其思路首先用Bostick反演的深度来控制层参数,使反演计算的模型参数仅存在电阻率;最后采用阻尼高斯-牛顿算法进行反演计算,并将Bostick反演结果作为反演计算的初始模型。通过模型试算,结果表明其处理速度快、解释直观,对野外大地电磁测深数据进行初步反演解释是可行的。 相似文献
11.
电各向异性在自然界中普遍存在,特别是沉积盆地中的部分岩层经过压实变质作用,表现出很强的各向异性导电性,采用基于各向同性模型对实测资料进行正反演解释必然会造成困难甚至结果的错误,只有基于各向异性理论的正反演解释才更为合理准确。本文通过自适应有限元对大地电磁各向异性进行正演计算,分别模拟水平、垂直和倾斜各向异性介质在不同偏转角和主轴电阻率下的响应结果。结果表明：自适应有限元能够在后验误差的控制下得到合理的网格,使计算结果更加接近解析解;在各向异性介质中,大地电磁TE极化模式的视电阻率和阻抗相位与垂直于层面的电阻率无关;二维电各向异性结构中,大地电磁TM极化模式响应结果总是由主轴上的电阻率在y轴方向上的分量所决定。 相似文献
12.
13.
《物探化探计算技术》2017,(3)
在山区进行音频大地电磁测深(AMT)时,电磁场易受地形起伏影响发生畸变。为研究起伏地形对AMT数据的影响,采用WinGlink软件基于有限差分算法的正演模块,建立二维起伏地形和平地地电模型,并对比分析了两个理论模型的正演响应。结果表明:TE模式的视电阻率曲线在高频段受干扰程度相对较大,在低频段影响减小甚至消失;TM模式的视电阻曲线变化正好相反;二者相位受地形影响则较小。对含异常体的带地形地电模型进行正演计算,在正演响应中加入2.5%的随机误差模拟实测数据,利用WinGlink软件基于非线性共轭梯度(NLCG)算法的反演模块,分别进行带地形和不带地形的TE、TM和TE+TM模式的二维反演。各反演结果的对比表明,带地形的二维反演结果明显优于不带地形的反演结果,能够较好地反映出异常体的位置、形态及电阻率值。依据得到的认识,对青海祁连隧道工程的AMT实测数据进行带地形二维反演,取得了理想的应用效果。 相似文献
14.
贵州镇远马坪“东方一号”岩体是我国首次发现的含金刚石原生矿母岩,马坪地区发现的岩体属于金伯利岩浆体系的浅部相,其深部可能存在规模较大的隐伏岩管或岩筒。为揭示马坪地区深部含金刚石隐伏岩管或岩筒的空间展布特征,在区内开展了80 m×40 m高密集网度的音频大地电磁勘探工作;利用三维正演技术模拟研究区纯地形响应并在实测数据中去除,得到的定性解释结果在一定程度上恢复了被静态效应扭曲的AMT阻抗相位不变量分布形态;使用AR-QN拟牛顿反演方法对数据进行三维反演,根据研究区岩性统计结果设定地下单元的电阻率变化区间,获得了可靠的三维电性结构;最后依据地表发现的岩筒、钻孔揭露的多条岩脉等地质资料对该地电模型进行精细解释,勾画出了隐伏岩管(岩筒)的形态,为区内下一步金刚石原生矿找矿方向及预测提供了地球物理依据。 相似文献
15.
《物探与化探》2015,(6)
在地-井瞬变电磁正演模拟中,前人均以总场为研究对象,分析响应特征。此次选择地-井瞬变电磁异常场为研究对象,对单异常体和多异常体以及在"定源异井"和"动源定井"两种观测方式下典型地电模型进行地-井瞬变电磁法三维正演模拟。首先推导了二次磁场的扩散方程,并用有限差分法进行数值计算。通过分析以地-井瞬变电磁异常场为研究对象,认为异常场能反映目标体空间几何位置参数。研究结果表明:在不同钻孔中接收的异常响应随时间推移而衰减减弱,但异常变化特征保持不变;在不同方位布置场源得到的异常响应差异明显,根据其测量结果可判定异常体空间位置信息。最后总结影响地-井瞬变电磁异常的主要因素有:目标体的埋藏深度、目标体与钻孔之间的耦合关系以及场源的位置关系。研究结果为相关资料解释提供理论参考。 相似文献
16.
积分方程法三维模拟井地电法并行算法研究 总被引：1，自引：1，他引：1
在积分方程法模拟三维电磁响应原理的基础上,对井地电法三维正演模拟并行算法做了研究。同时,还利用Fortran 6.5和消息传递接口(MPI)开发了井地电法三维正演模拟并行计算程序。在求解大型线性方程组时,采用的是稳定的双共轭梯度(BICGSTAB)算法,单机正演程序与并行正演程序的计算结果有着很好的一致性,这说明了并行算法的正确性,为进一步开发井地电法三维反演并行程序奠定了基础,也为井地电法三维正、反演程序实用化提供了技术支持。 相似文献
17.
本文以实例结合数值模拟说明CSAMT中非平面波效应对资料解释的重要影响。实际数据测量于河南卢氏县,在同一条测线测量,发射源长度1.3 km,收发距离分别为4.5 km和7 km。通过对数据进行必要的校正后进行了光滑模型反演,为了对比反演结果,整个处理过程的参数都保持一致。反演的电阻率断面显示,两种收发距离反演的电阻率断面其异常形态基本一致;但在电阻率的数值上,4.5 km的收发距离反演的电阻率较7 km的大。分析发现上述现象主要是源的非平面波效应造成的,并通过三维正演数值模拟得到了验证。由于在有限收发距离时,计算卡尼亚视电阻率靠近源受非平面波效应影响偏大,远离源后,该视电阻率趋于平面波视电阻率。因此,为了采用目前的可控源音频大地电磁测深解释理论,要求收发距离尽可能大。否则,必须采用带源效应的真正可控源电磁反演技术完成数据分析。 相似文献
18.
19.
给出各向异性介质二维地电断面大地电磁场的边值问题以及等价的变分问题。对计算区域采用矩形网格中进一步三角细化的剖分方式并在三角单元内进行线性插值,解出有限单元法数值解。通过典型模型的正演模拟,得到大地电磁测深曲线,并与前人的研究工作对比,验证了该方法的有效性。 相似文献