2.5-D稳定电流场h-自适应有限元的后验误差分析

本文针对2.5-D稳定电流场问题,用Galerkin方法推导了其变分问题;基于残值推导了h-自适应有限元方法求解2.5-D稳定电流场边值问题的后验误差公式,阐述了2.5-D稳定电流场的基于超收敛恢复的后验误差估计方法.为2.5-D稳定电流场h-自适应有限元方法的后验误差的实现奠定了理论基础.     

基于非结构化网格的2.5-D直流电阻率自适应有限元数值模拟

2.5-D直流电阻率有限元数值模拟中,模型的剖分及加密主要通过手动实现.另外,采用的单元类型比较规则如矩形单元等,不易实现复杂模型的模拟.为解决上述问题,文中提出了一种自适应有限元算法.算法中采用稳健的后验误差估计来自动预测下一次网格的单元尺寸,直到设定的迭代条件满足为止.另外,采用非结构化三角形单元实现了任意复杂模型的灵活剖分.基于此,利用垂直接触面模型分析和对比了不同自适应策略的效率.通过对比发现,点源附近的单元得到了加密以消除源的奇异性.另外,对于任意一种策略,有限元结果均能最终收敛到精确解.最后,模拟了两个模型：2-D单个异常体模型和2-D地形模型.     

可控源电磁场三维自适应矢量有限元正演模拟

本文实现了可控源电磁（CSEM）场三维自适应矢量有限元正演算法,该算法采用非结构四面体单元进行三维网格剖分,能够真实模拟地形起伏和复杂电性异常体.采用一次场和二次场分离的方式计算电磁场响应,能够有效解决有限元模拟中的源点奇异性,提高场源附近电磁场数值精度,其中一次场利用CSEM一维正演算法解析求得,二次场采用矢量有限元方法求得.并利用基于后验误差估计的自适应网格细化算法指导网格细化,以减少人为设计网格导致的误差.通过一维和三维模型的数值模拟,验证了本文算法的有效性：一维模型有限元解与解析解吻合得很好,电磁场振幅相对误差在1%左右,相位差整体小于1°;三维模型有限元解与有限体积解吻合得也很好.模拟了一个含三维倾斜板状异常体的可控源电磁响应,表明了本文算法模拟复杂地电结构电磁场的能力和有效性.

     

基于自适应有限元的复电阻率法2.5维复杂构造电磁场模拟

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

     

基于Coulomb规范势的电导率呈任意各向异性海洋可控源电磁三维非结构化有限元数值模拟

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

     

基于电流密度连续性条件的直流电阻率各向异性问题自适应有限元模拟

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

     

基于自适应有限元的起伏地形MT二维正演模拟

本文利用自适应有限元方法模拟起伏地形对大地电磁数据的影响.其中非结构化网格的自动剖分采用基于后处理技术的后验误差方法(又称Z-Z方法)进行控制、加密;通过地堑、山脊、地垒、山谷四种模型讨论了起伏地形对大地电磁TE和TM两种极化模式的视电阻率和阻抗相位曲线的影响,其中仅取10 Hz(高频)和0.01 Hz(低频)两个频点进行分析.证实了复杂的地形条件对TM模式的影响要比对TE模式的影响大得多,在TM极化模式中,低频部分的视电阻率比高频部分的更容易受起伏地形影响,而在TE极化模式中,则刚好相反等重要结论,为实际工作提供了理论指导作用.     

一种激发极化法2.5维正演的自适应有限元方法（英文）

传统的基于结构化网格有限元法采用的单元比较规则如矩形等,且网格剖分和加密要靠手动实现,所以传统的基于结构化网格有限元法不能准确和灵活地模拟复杂介质。本文采用易于模拟复杂介质模型的非结构化三角形网格进行剖分,且利用对偶加权后验误差估计指导网格自动细化过程,然后在电位模拟的基础上计算雅可比偏导矩阵,并依据Seigel（1959）理论实现激发极化法2.5维自适应有限元正演模拟算法。通过对垂直接触面模型进行正演分析,接收点附近网格得到了明显加密,电位数值解平均相对误差收敛到0.4%,视极化率平均相对误差收敛到1.2%,表明经自适应网格细化后,该算法数值解最终能收敛到精确解附近。最后对两个较复杂模型进行了正演计算与分析,进一步验证了该算法的准确性和灵活性。     

面向目标自适应有限元法的带地形三维大地电磁各向异性正演模拟

传统三维大地电磁各向异性模拟均是基于规则六面体网格,计算精度有限且较难拟合复杂地质条件.本文采用面向目标自适应非结构矢量有限元法,对三维大地电磁各向异性介质进行模拟.首先从电场双旋度方程出发,利用伽辽金方法建立变分方程;然后利用电流密度连续性条件构建适合大地电磁各向异性问题的加权后验误差估计方法,实现面向目标的网格自适应正演;最后通过典型算例分析各向异性对网格自适应和大地电磁响应的影响特征以及各向异性的识别方法.本文算法能够高精度地拟合起伏地表和任意各向异性介质,适用于分析复杂地电条件大地电磁响应特征,为提高大地电磁资料解释水平提供了理论基础.

     

数值模拟静态应力扰动下的断层失稳:结果分析兼与Dieterich模型和Coulomb模型的对比

基于大量岩石力学实验,Dieterich和其他研究者(Dieterich, 1978; Ruina, 1983)首先提出了描述岩石摩擦过程的速率-状态摩擦定律(R-S摩擦定律).如今R-S摩擦定律已成为研究地震成核等地震演化机制的有效手段.Dieterich (1992, 1994)最早提出了描述受静态剪应力扰动后断层失稳时间提前或推后的余震触发机制的解析模型.现在Dieterich模型已经成为解释余震随时间衰减规律的Omori定律等地震观测现象的有力工具.与之相对应,广泛使用的Coulomb应力失稳模型也可以给出断层受到静态剪应力扰动后,断层失稳时间的提前和推后量.Dieterich模型和Coulomb应力失稳模型基于不同的物理方法,所以在进行地震危险性评估时,二者均有各自的局限性.本文利用R-S摩擦定律控制的1-D弹簧-滑块模型,模拟计算了理论地震循环以及在不同静态剪切应力扰动下,失稳时间的提前和推后量的变化情况,然后将计算得到的时间提前和推后量分别与Dieterich模型和Coulomb模型的相应计算结果进行了定量化对比和差异性分析,并给出了相应的解释.数值模拟结果显示,对于R-S摩擦定律在参数不变的条件下,断层模型失稳时间的提前和推后量的大小强烈依赖于静态剪应力扰动的大小和作用时间,而且绝对值相同的正、负向静态剪应力扰动造成的失稳时间的提前和推后量的变化情况并不完全一致.在震后松弛/滑移阶段和闭锁阶段,时间提前和推后量是常数,且随静态剪应力扰动绝对值的增大而增大,两者的比值接近于1.0,这与Dieterich模型和Coulomb模型的结果是一致的,相应的差值小于两模型结果的10%.而在自加速阶段,模拟计算结果则存在与Dieterich模型和Coulomb模型结果不同的特征.首先,在自加速阶段模拟计算结果均偏离Coulomb模型,而且时间提前和推后量的比值小于1.0,相异于Coulomb模型的论断.不过当受到正向静态剪应力扰动后,Dieterich模型的结果和模拟计算结果是一致的,最大相差量不超过Dieterich模型结果的7%,可接近0.对于负向静态剪应力扰动,当其绝对值较小时,Dieterich模型的结果很接近模拟计算结果,相差量不超过该结果的14%.但对于绝对值较大的情况,模拟计算结果远大于Dieterich模型的结果,最大可达Dieterich模型结果的35倍,这是由于负向静态剪应力扰动后使得Ω=δθ/dc>>1的条件不再成立,进而使得Dieterich模型不再成立.总的来说,与模拟计算相比Dieterich模型可以很好地描述1-D断层受扰后失稳时间提前和推后量的变化情况,并且可以体现出正、负静态剪切应力扰动后失稳时间提前量和推后量变化的差异性,而Coulomb模型则不能完整地给出受到静态剪应力扰动后断层失稳时间提前或推后的估计值.     

基于自适应有限元正演的大地电磁法三维反演算法研究

目前大地电磁法三维反演大多基于有限差分正演.反演使用规则六面体网格无法有效模拟复杂地形, 同时正反演网格同套存在严重影响反演可靠性的问题.针对上述两个问题, 本文利用自适应有限元算法, 开发了基于非规则六面体的高精度大地电磁三维正演方法; 反演中, 使用独立的反演网格和正演网格来提高正演计算精度和反演可靠性.这一思路既保证了正演响应和灵敏度矩阵计算的精确性, 又降低了因反演参数过多造成的不唯一性.最后, 通过算例验证了正演算法的精确性和反演算法对地形处理的有效性.

     

有限元模拟弹性位错的等效体力方法

计算地震位错造成的位移场和应力场,对于估计大地震引起的后续地震活动发展趋势十分重要.本文提出在弹性位错问题的有限元模拟中,用等效体力代替位错源,从而在构建几何模型时不用包含断层,却可以处理包含任意复杂断层的问题,极大降低建模的难度.使用此方法,本文计算并讨论了在球形地球模型下2011年日本Tohoku-Oki特大地震对华北地区断层的影响,结果表明此次地震使华北主要断层趋于稳定.     

复电阻率法二维有限元数值模拟

伴随着复电阻率法的广泛应用,发展精确和快速的正演和反演算法成为复电阻率法研究的重点.本文采用基于三角单元剖分的有限单元法进行了复电阻率二维数值模拟研究.为了提高计算速度,对无穷远边界进行了近似处理.整个正演计算过程分为两步,首先采用有限单元法计算四个不同频率的视复电阻率数据,然后对前一步得到的视复电阻率数据采用递推算法计算视Cole-Cole参数.采用这种正演算法与一维正演的结果进行了对比,验证了本文方法的正确性.设计了两个二维极化模型,数值模拟结果表明视复电阻率和Cole-Cole视参数等值线断面图对于异常目标体都有比较明显的反映.     

自适应非结构有限元MT二维起伏地形正反演研究

在山区进行MT勘探时,用规则网格有限元方法模拟起伏地形会受到限制.本文采用非结构三角网格可以有效地模拟任意二维地质结构,如起伏地形、倾斜岩层和多尺度构造等.正演引入自适应有限元方法,其在网格剖分过程中能根据单元误差自动细化网格,保证了正演结果的精度.将自适应有限元与Occam算法结合,且引用并行处理技术提高正反演计算速度.通过对比两个理论模型,讨论了地形对MT正演响应的影响;其次进行了不同地电模型带地形反演展示了本文算法的正确性和适用性;最后将该方法应用于实测MT数据处理,证明了自适应非结构有限元方法是复杂地形下处理MT数据的有力工具.     

可控源电磁三维频率域有限元模拟

本文采用电磁场的磁矢量位和电标量势,将Maxwell方程组化为位势的类似于Helmholtz型方程,并引入罚项及稳定化方法克服了电磁三维有限元算法中的伪解及数值不稳定性;采用人工边界把计算区域局域化,将均匀半空间中水平电偶极子源产生的位势值作为人工边界上的第一类边界条件以表示源的作用,减少了实际的计算区域.理论模型和复杂模型的计算结果均表明,可控源电磁三维有限元数值模拟给出了稳定、可靠的电磁场分布.     

三维三分量CSAMT法有限元正演模拟研究初探

首先从麦克斯韦方程出发,用伽里金方法推导了三维三分量CSAMT法的有限元方程.在研究过程中,认识到加入散度条件的必要性,在公式中强加了散度条件,提高了解的完备性.其次将成功应用于二维线源频率域电磁法有限元模拟中的两种技术推广到三维中,一是边界条件统一采用一阶吸收边界,使线源产生的电磁波在边界上按波的传播规律被吸收,以降低平面波假设造成的影响;二是总体系数矩阵的存储,用两个二维数组分别记录总体系数矩阵的非零元素及其在总体结点编号中所处的位置,使总体系数矩阵的存储量达到最小的同时,物理意义明确,迭代求解时迅速简便.最后用均匀半空间模型进行了验证.