首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到16条相似文献,搜索用时 125 毫秒
1.
现有海洋可控源电磁三维数值模拟方法大多基于电导率各向同性介质理论,不能模拟海底地层电导率各向异性的实际情况.本文给出了电导率各向异性三维介质中电性源海洋可控源电磁二次电场的边值问题以及相应的变分问题,采用长方体单元对研究区域剖分,将场分量定义在剖分单元的边上,利用矢量有限单元法求解变分问题,实现了电导率任意各向异性海洋可控源电磁三维矢量有限元数值模拟.这个新的正演方法可以计算电导率任意各向异性三维地电模型的海洋可控源电磁响应,基于二次场矢量有限元法直接求解电磁场,避免了传统有限元方法可能遇到的伪解问题和难于处理电场法向分量不连续的问题,提高了数值模拟计算精度.一维电导率各向异性模型电磁场数值解与解析解吻合得相当好,无论在源附近还是远离源处相对误差均不超过1%.电导率各向异性二维模型的计算结果与已有文献采用的非结构有限元模拟结果十分吻合.三维地电模型数值模拟结果显示,电导率各向异性张量电导率主轴分量和欧拉角对不同装置海洋可控源电磁响应均有着明显的影响.  相似文献   

2.
目前,对于可控源电磁法各向异性介质2.5维问题,主要采用一次场、二次场分离的方法消除场源奇异性并降低截断边界对计算区域的影响.该方法数值计算精度高,但是很难适用于复杂地形条件下的数值模拟.针对复杂地形问题,基于总场的有限元方法表现出一定的优越性,然而,这种方法存在场源奇异性问题和截断边界问题.本文采用基于总场计算的方法对带地形的可控源电磁法2.5维各向异性介质进行模拟研究,推导了考虑电导率和介电常数各向异性的2.5维控制方程;引入网格加密-收缩算法降低场源奇异性的影响范围,提升数值计算效率;引入行波分解边界条件降低截断边界的影响;提出任意采样反傅里叶变换方法,快速、高精度地计算出空间域电磁场分量.理论模型数值算例中:首先,验证了本文算法的有效性;其次,对任意各向异性倾角产生的可控源电磁响应规律进行研究;最后,采用山丘模型对各向异性介质电磁场的响应规律进行了模拟和分析.  相似文献   

3.
在前人工作的基础上,本文推导了电导率任意各向异性介质的海洋可控源电磁三维谱元法正演方程.采用一次场/二次场分离算法结合混合阶矢量基函数,可以有效避免源点的奇异性的影响,从而提高数值解的精度.采用任意六面体单元离散研究区域,有利于模拟复杂地形和地电结构.利用不完全LU分解的Induced Dimension Reduction(IDR(s))迭代算法求解线性方程组,有效地提高了求解的效率.设计典型的地电模型进行正演计算,并将计算结果与有限元解进行对比,对比结果表明本文提出的基于混合阶矢量基函数的海洋可控源电磁三维谱元数值模拟算法是正确的、有效的.本文算法具有良好的通用性,可推广用于电导率呈任意各向异性的陆地电磁、井中电磁等数值模拟研究.  相似文献   

4.
李刚  李予国  韩波  段双敏 《地球物理学报》2017,60(12):4887-4900
在海洋可控源电磁法勘探中,接收站常置于海底.在进行海洋电磁场模拟时,由于海水和海底介质存在显著电性差异,这给海底接收点处场值的求取带来困难.本文提出一种新的接收点插值算法,该算法考虑到海底电场法向分量不连续性问题,用法向电流分量进行插值以准确求取海底任意接收点处电磁场值.本文利用交错网格有限差分法实现了二维介质中频率域海洋可控源法(CSEM)正演.对构造走向做傅里叶变换,将三维电磁模拟问题转换为波数域2.5维问题,即三维场源激励下针对二维地电模型的电磁模拟问题.使用交错网格有限差分法,基于一次场/二次场分离方法导出波数域二次电场离散形式,并进一步求得波数域电磁场.采用本文提出的改进的插值算法可求得海底任意接收点处波数域电磁场,采用傅里叶逆变换对波数域电磁场进行积分可得到接收点处空间域电磁场.模型算例表明,与常规的线性插值和严格插值算法相比,本文提出的改进的插值算法具有更高的精度.  相似文献   

5.
本文以非结构化网格剖分为基础,推导了电导率呈任意各向异性条件下的海洋可控源电磁法二次场磁矢量位、标量位所满足的有限元方程.将不完全LU分解预条件因子(ILU)与Induced dimension reduction(IDR(s))迭代算法相结合对有限元方程进行求解.之后,采用加权移动最小二乘法对二次场矢量位、标量位进行求导得到相应的电磁场各分量.最后,开展了两个地电模型的电磁场计算与分析.结果表明:本文所提出的电导率呈任意各向异性条件下的海洋可控源电磁法三维正演算法正确;不完全LU分解预处理的IDR(s)算法计算效率明显优于常用的ILU-QMR、ILU-BICGSTAB算法;算法具有良好的通用性,可用于陆地电磁、航空电磁、井中电磁等地球物理勘探方法在电导率任意各向异性方面的研究.  相似文献   

6.
罗鸣  李予国 《地球物理学报》2015,58(8):2851-2861
地下介质的电阻率常常表现为各向异性,海底裂隙地层和层状沉积序列可能形成宏观电阻率各向异性.在解释海洋电磁资料时,电阻率各向异性的影响不应该被忽略,否则可能会得到错误的海底地电模型.作者编写了电阻率任意各向异性一维层状介质海洋可控源电磁场计算程序,计算了电阻率各向异性层状模型的海洋可控源电磁响应,讨论了覆盖层和高阻储层分别具有电阻率各向异性时的电磁场响应特征.  相似文献   

7.
李勇  林品荣  刘祖鉴 《地球物理学报》2019,62(10):3923-3933
考虑地球介质电导率任意各向异性且随空间位置连续变化的情况,本文实现了直接求解电磁场的可控源音频大地电磁测深(CSAMT)三维有限元数值模拟.首先给出了电导率任意各向异性介质中CSAMT二次电场满足的控制方程及其相应变分问题,然后采用任意六面体单元对研究区域进行剖分,在网格单元中对任意各向异性电导率进行线性插值,解决了实际工作中岩矿石电导率各向异性且连续变化的情况,将变分问题转化为线性代数方程组的求解.电导率各向异性且连续变化一维模型三维有限元数值模拟结果与电导率各向异性且分层均匀渐进模型解析解结果对比验证了方法的有效性;三维地电模型电导率随位置线性变化且各向同性、主轴各向异性、方位各向异性和倾斜各向异性的数值模拟结果表明,电导率各向异性且连续变化对CSAMT视电阻率和相位数据均有明显的影响.  相似文献   

8.
三维任意各向异性介质中海洋可控源电磁法正演研究   总被引:20,自引:15,他引:5       下载免费PDF全文
殷长春  贲放  刘云鹤  黄威  蔡晶 《地球物理学报》2014,57(12):4110-4122
由于海底介质受沉积环境的影响,层理发育呈现明显各向异性特征.对于海洋可控源电磁法各向异性的研究以往主要局限于一维和二维模型,为更深入了解复杂情况下海底各向异性对海洋可控源电磁响应的影响规律,本文开展三维任意各向异性介质中海洋可控源电磁法正演研究.采用交错网格有限差分技术,通过对任意各向异性介质电导率张量实行体积和空间电流密度平均,完成海洋可控源电磁二次散射电场的离散化,成功实现任意各向异性介质中海洋可控源电磁正演模拟.通过对几种典型各向异性电性模型条件下海洋电磁电场多分量响应及分布特征和各向同性情况的对比分析,总结电各向异性对海洋电磁响应的影响规律和识别方法.本文算法研究及算例可为海洋可控源电磁数据精细化处理解释提供技术支撑.  相似文献   

9.
油气藏自身会产生激电效应,明确海洋可控源电磁场激电效应特征非常重要。本文基于自适应有限元方法实现了带激电效应的2.5D海洋可控源频率域电磁场正演。采用非结构化三角网格单元剖分地电模型,可以模拟任意起伏地形和复杂地电构造。利用Cole-Cole模型表述油气藏的激电效应,实现了考虑激电效应的海洋电磁复电阻率正演算法。将本文的计算结果与已发表的一维模型结果对比,检验了本文算法的正确性。最后,设计了不同的极化模型,针对不同的场源方位角、激电参数和地形,计算了的电磁场响应结果,分析了激电效应对海洋可控源电磁场的影响规律,结果表明油气藏的激电效应和海底地形对海洋可控源电磁场有一定的影响。  相似文献   

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

11.
电导率各向异性的海洋电磁三维有限单元法正演   总被引:10,自引:8,他引:2       下载免费PDF全文
本文提出了一种基于非结构化网格的海洋电磁有限单元正演算法.为了回避场源奇异性,文中选用二次场算法,将背景电阻率设置为水平层状且各向异性,场源在水平层状各向异性介质中所激发的一次场通过汉克尔积分得到.基于Coulomb规范得到二次矢量位和标量位所满足的Maxwell方程组,通过Galerkin加权余量法形成大型稀疏有限元方程,采用不完全LU分解(ILU)预条件因子的quasi-minimum residual(QMR)迭代解法对有限元方程进行求解得到二次矢量位和标量位;进而,利用滑动平均方法得到二次矢量位和标量位在空间的导数,由此得到二次电磁场;通过一维模型对算法的可靠性进行验证,与此同时,针对实际复杂海洋电磁模型,比较有限元模拟结果与积分方程模拟结果,进一步验证算法精度.若干计算结果均表明,文中算法具有良好的通用性,适用于井中电磁、航空电磁,环境地球物理等非均匀且各向异性介质中的电磁感应基础研究.  相似文献   

12.
频率域地空电磁探测方法是指在地面布设人工场源,在空中测量电磁场的一种高效的地球物理勘探技术.该方法具有大范围、高通过性、快速测量的优势,尤其适合崎岖山地、沙漠、沼泽、海陆交互带等复杂地貌区域的资源勘查.但是这些地区的地形起伏通常较大,因此分析地形对地空电磁响应的影响具有重要意义.本文利用有限元法对频率域地空电磁响应进行了正演计算,分析了起伏地表条件下的频率域地空电磁响应特征.首先利用傅里叶变换将2.5维问题转化成二维问题,利用伽辽金加权余量法推导了相应的离散有限元方程组.采用任意四边形单元对区域进行不均匀网格剖分,源和异常体附近网格加密处理,保证计算精度,远离目标区域网格逐渐稀疏,模拟无穷远边界,降低对计算资源的要求.在单元内进行插值,将有限元方程组变换为线性方程组,采用总场算法,利用具有一定面积的伪δ函数表达源电流分布,源项近似为分布在以电偶极源为中心的25个节点上.通过求解线性方程组得到波数域电磁响应,再对波数域电磁场响应进行反傅里叶变换从而获得空间域2.5维频率域电磁场值.通过对比2.5维正演结果与均匀半空间解析解,验证了本文算法的精度,同时本文还对地空电磁场与地面电磁场的响应特性进行了对比.  相似文献   

13.
为了分析海洋中含金属矿的储油气砂层的成分、结构对于海洋可控源电磁场的影响特征,引入多相微增等效介质模型,结合非结构化网格实现了海洋可控源电磁三维矢量有限元正演.首先对多相微增模型进行了介绍,并分析不同参数影响下的等效电导率变化特征.而后,详细推导了频率域可控源电磁法非结构化矢量有限元方程.接着,利用预处理的IDR (s)迭代算法求解线性方程组.最后设计典型的海洋地电模型,通过正演计算,验证了本文算法的正确性和有效性,同时研究了海洋储油砂层中的金属矿的含量、孔隙度、含水饱和度等参数对于可控源电磁场响应特征的影响.  相似文献   

14.
罗鸣  李予国  李刚 《地球物理学报》2016,59(11):4349-4359
本文提出了一维垂直各向异性(VTI)介质倾斜偶极源频率域海洋可控源电磁(CSEM)资料高斯-牛顿反演方法.在电阻率各向异性介质水平偶极源和垂直偶极源海洋CSEM正演算法的基础上,利用欧拉旋转方法,实现了各向异性介质倾斜偶极源海洋CSEM正演算法.海洋可控源电磁场关于地下介质横向电阻率(ρ_h)和垂向电阻率(ρ_v)的偏导数(即灵敏度矩阵)是解析计算的,结合垂直各向异性介质横向电阻率与垂向电阻率的关系,将各向异性率融入到正则化因子选择中,实现了正则化因子的自适应选择.理论模型合成数据和实测资料反演算例表明,我们提出的反演方法能够较准确的重构海底围岩和基岩的各向异性电阻率以及高阻薄层的埋藏深度、厚度和垂向电阻率.  相似文献   

15.
海洋可控源电磁三维非结构矢量有限元数值模拟   总被引:14,自引:12,他引:2       下载免费PDF全文
杨军  刘颖  吴小平 《地球物理学报》2015,58(8):2827-2838
本文实现了海洋可控源电磁三维矢量有限元数值模拟.由于采用非结构四面体单元进行三维网格剖分,该方法可以模拟复杂电性异常体和海底地形.一维模型的数值模拟结果表明,电场实、虚部均与解析解吻合得相当好,计算误差基本小于1%.二维模型的计算结果与已有的二维自适应非结构有限元模拟结果吻合很好.带地形模型的数值模拟结果显示,海底地形对电场影响很大,有可能掩盖海底油气藏产生的异常.  相似文献   

16.
频率域激电有限元数值模拟   总被引:7,自引:4,他引:3       下载免费PDF全文
在三维地电条件下,应用有限元方法模拟了不同频率谐变电流激发下的地表激电响应,首先根据麦克斯韦方程推导出谐变电磁场复电位的波动方程,利用广义变分原理得到复电位波动方程及其边界条件下边值问题等价的变分问题.尔后将区域剖分离散化,在单元中进行插值,得到一组线性方程组,通过合理存储刚度矩阵、解方程组,得到各个结点上的复电位值,最后得到表征频率域激电响应的幅频率等参数.模拟结果表明该方法合理正确.  相似文献   

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

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