首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到18条相似文献,搜索用时 140 毫秒
1.
起伏海底地形时间域海洋电磁三维自适应正演模拟   总被引:1,自引:0,他引:1       下载免费PDF全文
本文基于自适应非结构有限元算法实现海洋电磁起伏海底地形三维正演模拟.通过采用隐式后推欧拉时间离散技术,保证在较大的时间步长条件下获得正确结果.为获得多时间道海洋电磁正演模拟的有效网格,我们采用基于法向电流连续的后验误差估计的自适应方法和网格融合技术;同时为了控制网格数量和保证正演模拟稳定性,我们还在网格融合过程中应用了随机网格挑选技术.对于方程组求解我们使用MUMPUS直接求解器.当时间步长不变时,只需对系数矩阵进行一次分解,大大提高计算速度.将本文计算结果与半空间模型解析解进行对比,验证了本文算法精度.针对海洋电磁半拖曳式和双船拖曳式工作方式,我们通过典型模型的模拟计算,研究海底地形影响及海底高阻层识别特征.  相似文献   

2.
徐凯军  李猛 《地球物理学报》2018,61(7):3102-3111
复电阻率法在矿产、油气勘探调查中发挥着重要作用.为了认识复杂构造的复电阻率法电磁场的变化规律,本文基于自适应有限元方法,采用非结构化网格,引入Cole-Cole模型,实现了电偶源2.5D复电阻率法电磁场正演,可以模拟复杂地形和地电结构,正演结果更符合野外实际地质情况.通过将本文的计算结果与半空间模型解析解、层状介质和起伏模型结果进行对比,验证了本文算法的正确性.最后,基于复杂地电模型,通过正演模拟,系统分析了地形、激电参数、复杂构造对复电阻率法电磁场的影响特征.  相似文献   

3.
面向目标自适应三维大地电磁正演模拟   总被引:3,自引:3,他引:0       下载免费PDF全文
本文将面向目标的自适应算法应用于三维大地电磁数值模拟.使用基于非结构网格的矢量有限单元法对起伏地表大地电磁正演模拟问题进行求解.使用利用垂向电流密度在物性界面上的连续性对后验误差进行估算的算法指导网格优化.由于全局自适应算法针对观测点优化网格的能力较差,本文通过求解正演问题的对偶问题计算后验误差的加权系数,并对相关加权系数进行改进,从而实现了面向目标的自适应算法.与传统基于结构化网格的电磁正演算法相比,采用非结构网格能够更好地拟合起伏地表和地下不规则异常体.由于使用了面向目标的自适应算法,本文能够使用更少的网格达到较高的计算精度.通过对比本文模拟结果与半空间响应和全局自适应算法计算结果,并通过对比使用改进前和改进后加权系数得到的网格剖分结果验证了本文算法的有效性.  相似文献   

4.
本文基于非结构网格实现了海洋可控源电磁法三维有限元正演模拟.该算法采用完全非结构网格剖分,可以模拟任意起伏地形和复杂地电模型.为了避免场源的奇异性,采用一次场/二次场分解算法,一次场由基于Schelkunoff势函数的一维解析公式得到.为了提高算法的精度和效率,采用对测点附近单元和异常体区域进行体积约束加密的方法,实现了非结构网格的局部加密.一、二维模型计算和分析表明,本文采用的局部加密方法能够明显地改善算法的精度,最大相对误差基本在1%以内.对三维模型计算及对比分析,说明了该算法对三维可控源电磁正演的实用性.复杂海底地形模型的正演模拟表明,海底地形对电磁场的影响很大,在进行海洋可控源电磁资料解释时,地形的影响有必要考虑在内.  相似文献   

5.
起伏地形和关断时间对地面瞬变电磁响应影响严重,这给传统基于水平地表模型和理论阶跃波形的瞬变电磁数据解释技术带来很大困难.为此,本文开展考虑起伏地形和关断时间的地面瞬变电磁三维反演算法研究.正演采用基于非结构网格和后退欧拉隐式时间离散格式的时间域有限元算法,快速模拟起伏地表模型瞬变电磁响应.反演采用L-BFGS算法,减少每次反演迭代的计算量.针对发射线圈随地表起伏变化的特点,利用基于偶极子离散的场源处理技术模拟发射源的实际形状,采用瞬时电流脉冲技术实现考虑关断时间的地面瞬变电磁三维正演模拟.我们首先将本文开发的三维反演算法应用于理论模型的反演计算中,检验本文算法的可靠性,并分析地形和关断时间对反演结果的影响特征.在此基础上,进一步将本文算法应用于实测数据反演,验证本文算法的实用性.  相似文献   

6.
可控源电磁场三维自适应矢量有限元正演模拟   总被引:1,自引:1,他引:0       下载免费PDF全文
刘颖  李予国  韩波 《地球物理学报》2017,60(12):4874-4886
本文实现了可控源电磁(CSEM)场三维自适应矢量有限元正演算法,该算法采用非结构四面体单元进行三维网格剖分,能够真实模拟地形起伏和复杂电性异常体.采用一次场和二次场分离的方式计算电磁场响应,能够有效解决有限元模拟中的源点奇异性,提高场源附近电磁场数值精度,其中一次场利用CSEM一维正演算法解析求得,二次场采用矢量有限元方法求得.并利用基于后验误差估计的自适应网格细化算法指导网格细化,以减少人为设计网格导致的误差.通过一维和三维模型的数值模拟,验证了本文算法的有效性:一维模型有限元解与解析解吻合得很好,电磁场振幅相对误差在1%左右,相位差整体小于1°;三维模型有限元解与有限体积解吻合得也很好.模拟了一个含三维倾斜板状异常体的可控源电磁响应,表明了本文算法模拟复杂地电结构电磁场的能力和有效性.  相似文献   

7.
直流电阻率法被广泛应用在工程和环境及水文地球物理、野外采矿、地热探测等领域.地下岩石常具有层理面和裂缝等具有方向依赖性的结构,岩石电导率常常呈现各向异性特征,因此研究复杂直流电阻率各向异性问题的高精度正演算法具有迫切的理论和学术需求.本文利用面向目标的自适应有限元算法和非结构化网格相结合的方式,解决了带地形任意复杂直流电阻率各向异性问题的高精度正演这一难题.有别于前人的研究成果,本文提出了一种特别的二次虚拟场算法实现带源的任意起伏地形问题模拟;另外,本文第一次基于电流密度连续性条件构建适合直流电阻率各向异性问题的后验误差估计算法,有效地驱动面向目标有限元网格自适应加密过程.最后,通过三组电阻率各向异性模型验证本文提出算法的正确性和适应性,测试结果表明:对于任意复杂直流电阻率各向异性问题,本文提出的算法具有精度高、适应性强等特点;另外,我们还发现电流密度连续性条件可用于设计直流电阻率问题的有效后验误差估计算法.  相似文献   

8.
基于非结构网格的电阻率三维带地形反演   总被引:6,自引:3,他引:3       下载免费PDF全文
吴小平  刘洋  王威 《地球物理学报》2015,58(8):2706-2717
地表起伏地形在野外矿产资源勘察中不可避免,其对直流电阻率法勘探影响巨大.近年来,电阻率三维正演取得诸多进展,特别是应用非结构网格我们能够进行任意复杂地形和几何模型的电阻率三维数值模拟,但面向实际应用的起伏地形下电阻率三维反演依然困难.本文基于非结构化四面体网格,并考虑到应用GPS/GNSS时,区域地球物理调查中可非规则布设测网的实际特点,实现了任意地形(平坦或起伏)条件下、任意布设的偶极-偶极视电阻率数据的不完全Gauss-Newton三维反演.合成数据的反演结果表明了方法的有效性,可应用于复杂野外环境下的三维电法勘探.  相似文献   

9.
各向异性介质模型电性结构复杂,如何进行合理的网格剖分成为获得高精度正演结果的关键,为此本文开展时间域航空电磁各向异性大地三维自适应有限元正演算法研究.通过结合非结构时间域有限元算法和自适应网格优化技术,实现各向异性介质条件下三维时间域航空电磁自适应正演.考虑到时间域航空电磁响应随时间的衰减特性,为了综合评价不同时刻的后验误差,本文将时间作为加权因子,调整各个时刻后验误差的相对权重,进而实现对浅部和深部网格的同步优化.通过与一维解析结果进行对比验证了本文算法的可靠性.数值实验结果显示电导率各向异性对自适应网格影响严重,其最大主轴电导率的数值及其分布特征直接决定了网格加密效果.此外,各向异性对时间域航空电磁三分量响应的分布形态和异常幅值也会产生严重影响,利用全域视电阻率极性图,可以很好地识别各向异性主轴方向.  相似文献   

10.
起伏地形下线源频率域电磁场二维正演模拟   总被引:1,自引:1,他引:0  
为了研究地形对线源频率域电磁测深测量结果的影响,本文介绍了起伏地形下的线源频率域电磁场的二维正演模拟方法.在正演模拟过程中采用有限单元法,首先为了模拟起伏地形,用三角形单元对研究区域进行网格剖分;然后再对每个单元进行线性插值、积分、总体合成;最后通过求解线性方程组得出地表的电磁响应.文章最后给出了山谷与山脊模型的数值模...  相似文献   

11.
研究了基于矢量有限元方法的大地电磁带地形三维反演算法并开发了三维反演计算程序代码.在大地电磁场正演数值模拟方面,采用并行直接稀疏求解器PARDISO且无需进行散度校正的快速正演方案,对典型地形模型,在中等规模计算条件下,与双共轭梯度法(BICG)计算结果比较,发现PARDISO比BICG快10倍以上;通过理论模型试算,并与前人的有限元法计算结果对比,验证了带地形三维正演计算程序的正确性.在反演方面,本研究基于共轭梯度方法编写了大地电磁带地形三维反演代码,为了避免直接求取雅可比矩阵,将反演中的雅可比矩阵计算问题转为求解两次“拟正演”问题,进而将PARDISO的快速正演方案应用于“拟正演”问题的求解,以提高反演计算效率.利用开发的反演算法对多个带地形地电模型的合成数据进行了三维反演,反演结果能很好地重现理论模型的电性结构,验证了本文开发的三维反演算法的正确性和可靠性.最后,利用该算法反演了某矿区大地电磁实测数据,反演得到的三维电性结构清晰地反映了研究区的地电特征,将反演结果与该区已有地质资料结合进行解释,应用效果明显,进一步验证了本文算法的有效性.  相似文献   

12.
We present a 2-D inversion code incorporating a damped least-squares and a minimum-model approach for plane wave electromagnetic (EM) methods using an adaptive unstructured grid finite element forward operator. Unstructured triangular grids permit efficient discretization of arbitrary 2-D model geometries and, hence, allow for modeling arbitrary topography. The inversion model is parameterized on a coarse parameter grid which constitutes a subset of the forward modeling grid. The mapping from parameter to forward modeling grid is obtained by adaptive mesh refinement. Sensitivities are determined by solving a modified sensitivity equation system arising from the derivative of the finite element equations with respect to the model parameters. Firstly, we demonstrate that surface topography may induce significant effects on the EM response and in the inversion result, and that it cannot be ignored when the scale length of topographic variations is in the order of magnitude of the skin depth. Secondly, the dependency of the inversion on the starting model is discussed for VLF and VLF-R data. Thirdly, we demonstrate the inversion of a synthetic data set obtained from a model with topography. Finally, the inversion approach is applied to field data collected in a region with undulating topography.  相似文献   

13.
电性参数分块连续变化二维MT有限元数值模拟   总被引:4,自引:2,他引:2       下载免费PDF全文
刘云  王绪本 《地球物理学报》2012,55(6):2079-2086
为了易于模拟野外复杂地形和地下任意形状地电体模型,将有限元单元网格设计为三角单元;并考虑到野外实际勘探中,地球介质的电性参数均是连续变化的情况,单元内的场值和电性参数被设计为双线性变化;推导出二维起伏地形条件下大地电磁法有限元数值模拟算法;根据单元节点主场值和线性插值形函数间的关系,计算出单元节点的辅助场值;在二维起伏地形情况下,定义TE、TM模式视电阻率和阻抗相位.4个模型的计算的结果与解析法的均方根误差小于1%,地形模拟与前人的计算结果相符,模拟倾斜界面异常体,能有效的反映出其异常形态.  相似文献   

14.
We describe an algorithm for inverting magnetotelluric (MT) data in the presence of strong bathymetry or topography. Instead of correcting distortions due to bathymetry or topography we incorporate them directly into the inversion. To achieve a high accuracy in computing MT responses we use finite‐difference approximations that permit sloping discrete boundaries inside model elements. The same approach is applicable to any seafloor electromagnetic calculation and can also be used to incorporate steep topography on land. We test our approach on various topographic features and compare our results to that of a finite element approach. Finally, we present inversion examples that illustrate the effectiveness of our inversion algorithm in recovering true subsurface structures in the presence of strong bathymetry and topography.  相似文献   

15.
射频大地电磁法(RMT)是以潜艇天线发射的射线源等作为场源的一种地球物理勘探方法,目前被广泛应用于近地表工程和环境地球物理勘探.RMT数据解释常采用基于静态假设的大地电磁法(MT)程序,往往会反演出不真实的浅层目标体.为解决这一长期困扰RMT资料解释的问题,本文实现了考虑位移电流的RMT有限元数值模拟.为了处理任意复杂模型,如起伏地形,非结构的三角形网格被用于离散RMT模型.首先通过算例对比,验证了程序的正确性和可靠性.通过Dike模型讨论了空气层厚度对RMT数值解的影响,结果表明当空气层厚度大于1/4波长即可满足精度要求.以山脊模型为例计算了位移电流对RMT响应的影响,表明位移电流的影响会随着山脊高程的增加而增大.最后通过舒家店实际模型进一步验证了位移电流在RMT中的重要性.  相似文献   

16.
基于有限差分正演的带地形三维大地电磁反演方法   总被引:4,自引:4,他引:0  
本研究实现了一套基于有限差分(FD)方法的大地电磁测深数据带地形三维反演算法及代码.其中,在大地电磁场正演数值模拟方面,开发了起伏地形条件下基于交错网格剖分、有限差分方法的大地电磁测深三维正演代码;在满足平面波场假设的前提下,使用长方体网格剖分模拟三维起伏地形,实现了带地形三维正演计算;并设计理论模型进行试算,经试算结果与前人的有限元法计算结果对比,验证了所研发的带地形三维正演计算的正确性与可靠性.在反演方面,本研究基于非线性共轭梯度方法编写了大地电磁测深带地形三维反演代码,试验了不同的共轭梯度搜索因子β,避免了目标函数对海森矩阵(参数二次导数矩阵)的显式计算和存储,初步实现了大地电磁资料的带地形三维反演.最后,对一系列理论模型进行正演计算,利用其生成的合成数据模拟实测数据进行反演,并与现有的不带地形大地电磁测深三维反演结果比较,检验了所研发的带地形三维反演计算的可靠性与稳定性.  相似文献   

17.
无网格局部Petrov-Galerkin法大地电磁场二维正演模拟   总被引:1,自引:0,他引:1       下载免费PDF全文
卢杰  李予国 《地球物理学报》2017,60(3):1189-1200
有限差分法和有限单元法在大地电磁场数值模拟中已经得到了广泛的应用,但其数值结果的精度在很大程度上依赖于网格的离散程度.当模拟起伏地形、弯曲界面等复杂地电模型大地电磁场响应时,常常需要花费大量的时间以便得到较合理的离散网格.无网格局部Petrov-Galerkin法(MLPG)不同于有限差分法和有限元法,其形函数和权函数脱离了网格的束缚.本文详细推导了二维大地电磁场边值问题的弱式形式,并将其离散为局部积分域内的表达形式.通过模拟二维海洋地电模型大地电磁场响应,并与结构网格有限元结果进行对比,验证了本文算法和程序的正确性及精度.设计了一个含有弯曲界面的二维地电模型,讨论了不同离散网格对MLPG无网格法模拟结果的影响,并与结构有限元法结果进行了比较,结果表明MLPG无网格法模拟结果受离散网格影响较小.最后利用MLPG无网格法计算了两个海洋起伏地形模型的大地电磁响应,讨论了海底起伏地形对大地电磁响应的影响.  相似文献   

18.
电阻率和速度随机分布的MT与地震联合反演   总被引:10,自引:5,他引:5       下载免费PDF全文
在已有研究成果的基础上,为了适应物性参数剧烈变化的复杂模型并满足联合反演的要求,开发了速度和电阻率随机分布共网格单元模型的建模技术.基于这种统一的物性随机分布的网格介质模型,利用有限元方法和改进的射线追踪法分别正演计算大地电磁场和地震走时,结合改进的模拟退火算法,研究实现了电阻率和速度随机分布条件下的大地电磁与地震资料的同步联合反演.对物性界面不完全一致和物性变化剧烈的带地形复杂模型的试验,表明了该方法在精细反演复杂电阻率和速度结构方面的效果,克服了以往研究局限于简单模型的不足.对地震资料品质差的地区开展的实际资料联合反演,表明了方法的适用性,先验信息约束下的联合反演提高了反演精度.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号