大地电磁三维矢量有限元正演

采用广义变分原理,基于矢量基函数详细推导了大地电磁三维矢量有限元方程。为了提高计算精度和效率,应用直接法强加边界条件改善总体系数矩阵的条件数,同时使用SSOR(symmetric successive over relaxation)预处理的双共轭稳定梯度法求解复对称大型稀疏线性方程组。并利用国际标准模型与相关参考文献的结果进行了对比,验证了算法的准确性。对一个典型的三维低阻体模型进行正演,得到了不同测线的视电阻率和相位断面图,并与二维正演结果进行对比分析。结果表明:在x方向测线上,ρ_(yx)变化幅度较ρ_(xy)小,中心测线上的ρ_(yx)和ρ_(xy)响应均与二维TM模式条件下的响应特征相似。     

五极纵轴激电测深三维有限元正演模拟

本文研究了五极纵轴激电测深三维有限元正演模拟方法。首先,从三个点电流源总电位的边值问题出发,导出了异常电位的边值问题,证明了与异常电位的边值问题对应的变分问题。然后,基于三维连续电性介质模型,推导了有限元法求解变分问题的计算过程,编制了五极纵轴激电测深的三维模拟程序。最后,通过模拟水平层状介质模型的五极纵轴激电测深曲线,发现最大相对误差小于0.25%,耗费时间为18 s;从精度和效率的角度验证了算法和程序是正确的,为后续模拟复杂地电模型五极纵轴激电测深曲线、分析曲线异常特征规律以及开展反演工作奠定了基础。     

矿井层测深三维地电体异常规律研究

依据电磁场理论，构建地电模型，利用有限差分对矿井三维地电异常体进行数值模拟，得到异常体在电阻率不同、至巷道距离不同、及测点不同情况下的层测深曲线，并对其分布规律进行分析研究。通过模拟矿井三维地质异常体地电异常响应特征，并与井下实测层测深曲线对比和分析，论证了矿井三维地电体数值模拟的正确性。     

水平电偶源频率电磁测深理论曲线的正演计算

水平电偶源频率电磁测深是采用水平接地电偶极子在地下激发电磁场,在地面上分别测量电场分量和磁场分量（也可测相应场分量的相位）,以提取反映地质构造的信息。     

三维地电断面电阻率测深有限元数值模拟

用有限单元法进行了电导率分块均匀的三维点源电场电阻率测深的数值模拟.首先给出了三维构造中点源电场的边值问题、变分问题; 然后, 用有限单元法求解变分问题, 将区域剖分成六面体单元, 在单元中进行三线性函数插值, 将变分方程化为线性代数方程组; 最后解方程组, 得各节点的电位值, 进而计算出地表的视电阻率.对几例较典型的地电模型进行试算, 结果表明本方法是行之有效的      

一种二维地电模型有限元网格自动形成方法

本文给出了一种二维地电模型有限元网格自动形成方法,并详细介绍了该方法的数学模型及程序框图。     

线源频率电磁测深二维正演(一)

随着频率电磁测深(包括可控源声频大地电磁测深 CSAMT)的深入发展,必然要进行三维情况下的研究。线源频率电磁测深二维正演模拟研究为最终实现三维正演打下了基础。     

线源频率电磁测深二维正演(一)

随着频率电磁测深(包括可控源声频大地电磁测深CSAMT)的深入发展,必然要进行三维情况下的研究。线源频率电磁测深二维正演模拟研究为最终实现三维正演打下了基础     

线源频率电磁测深二维正演(二)

随着频率电磁测深（包括可控源声频大地电磁测深ＣＳＡＭＴ）的深入发展,必然要进行三维情况下的研究。线源频率电磁测深二维正演模拟研究为最终实现三维正演打下了基础。     

基于矢量有限元的大地电磁快速三维正演研究

本文采用直接求解器PARDISO且无需散度校正的正演方案,求解矢量有限元法对应的大型线性方程组,获得不同地形(水平和起伏)条件下地电模型的大地电磁响应,提高了大地电磁三维正演计算的速度。在中等规模计算条件下,通过本文的计算方法与带散度校正的迭代求解法对比,计算速度可提高十倍以上。     

频率测深二维地形影响的边界元素法正演模拟

本文用边界元素法模拟二维地形对频率测深的影响，给出了典型地形赤道偶极装置ＥＸ分量的纯地形影响视电阻率曲线，列举了用此值法进行地形改正的实例。     

非均匀介质中微侧向测井仪器的三维有限元数值分析

在这篇文章中，分析了一个具有井眼的非均匀介质中似稳电流场的边值问题。采用了三维有限元的直接解法求解三维似稳电流场；给出了微侧向测井仪器的电极表面电流分布、假几何因子以及８英寸井眼的泥饼校正曲线；求出了在具有倾斜地层条件下，微侧向测井仪器对界面的响应。数值结果与发表的文献非常吻合，从而证明了这种方法的正确性。因此，它为新的仪器研制和微侧向测井解释提供了可靠的数据和理论依据。     

电法测井的三维有限元模拟

研究了不同测量方式(井-地,地-井,井-井)下点源场井中电法的三维有限元数值模拟.考虑到深度方向上大范围的网格剖分和井眼的影响及井-井测量等因素,采用放射状三棱柱单元的网格剖分方式,以提高网格质量,减少剖分单元数;给出了三棱柱单元的坐标变换公式,进行精确的单元积分,减少了单元积分时间;结合非结构化网格技术,实现了复杂模型的模拟;开发出相应的程序实现了复杂条件下(如考虑井眼影响、井井测量、倾斜井情形、地形起伏等)电法测井的三维有限元模拟,数值算例验证了方法的可靠性及计算效率,并对不同情形下的异常响应进行了分析,为进一步的反演工作奠定了基础.     

水平多层介质上水平电偶源频率电磁测深的阻抗实部等效电阻率

本文讨论了水平电偶源频率电磁测深中阻抗实部等效电阻率的具体定义和算法。用该方法计算的等效电阻率是电磁响应的单值函数,并极大地改善了视电阻率随频率变化的关系曲线,能直观地反映出地电断面随深度的客观变化,数值较接近地层真电阻率,且接近速度很快,假极值效应明显压低。数值计算和实际应用表明该方法是灵活、可行的,可很好地消除频率电磁测深中波区视电阻率曲线在非波区的畸变效应。     

频率域中有限元素法的探讨

本文指出了用拉普拉斯变换处理有限元方程和用富里叶变换所得到的结果是等价的。从物理方面讨论了频率域有限元的局限性。     

地下水模型的Neumann展开Monte-Carlo随机有限元法

讨论了求解地下水模型(水流模型和水质模型)的Neumann展开Monte-Carlo随机有限元法。从基本的随机变量入手,避免了过程中随机变量的增多问题,给出结点水头(浓度)的均值、方差和水头在某区间的概率计算方法;改进了矩阵求逆的效率,对输入随机变量较多、随机变量变异较大的非稳定地下水问题特别有效。同时选取二维承压地下水水流问题(有解析解)作为例子,进行了随机数值模拟实验。     

崇文门站施工允许地表沉降及三维有限元模拟分析

建设中的北京地铁五号线崇文门站下穿既有地铁一号线区间隧道,为保证既有地下铁道正常运营和地下结构的安全,需严格控制新建车站施工引起的地层位移。结合新建车站暗挖施工及既有线地铁的实际情况,分析和计算穿越既有隧道施工地表允许沉降值。为了控制地表沉降,确保施工安全,拟采用小刚度管棚或大刚度顶管作超前预支护,用3D-S igm a三维有限元软件进行施工效应的计算模拟,掌握大管棚和顶管预支护洞室的力学效应,预测车站施工引起既有隧道的沉降量。最后分析比较两种预支护的支护效果,为选择合理的预支护方案提供参考。     

用异常电位边界单元法做电测深资料地形改正

求解异常电位的边界单元法对电测深资料进行地形改正,应用效果明显。在异常电位法中将异常电位与正常电位相加,得总电位,用这种方法得到的总电位比直接求总电位的精度高得多。     

贮量集中有限元方法在渗流计算中的优点

普通有限元方法和有限差分方法是渗流计算中常用的数值方法,而贮量集中有限元方法则很少被实际应用。为了使渗流计算者全面了解贮量集中有限元方法,并能够较多地应用这种方法,本文将渗流计算中的贮量集中有限元方法、普通有限元方法和有限差分方法进行了对比,具体比较和分析了有限元方程组的系数矩阵特征、计算工作量,通过具体的算例分析了它们的计算精度,由此指出贮量集中有限元方法的计算精度与有限差分法基本相同,而远好于普通有限元法,它的计算工作量小于其他二种方法。