网格剖分及边界对MT Occam反演的影响研究

在大地电磁有限元模拟计算中网格剖分及边界放置是否得当,影响着最终的反演结果。基于Occam反演法从网格剖分及边界方面,设计不同模型,讨论不同网格对反演精度的影响。研究表明,对于二维介质体,只要将网格边界置于适当倍数趋肤深度的距离处,便可忽略截断边界的影响,无需放置在无穷远。在同一频段及不均匀剖分的前提下,粗细网格均可获得良好的反演结果。细网格的反演精度高于粗网格,但未能很好地圈定异常体的边界。粗网格剖分下,MT Occam法对低阻异常体反演效果较好,而对高阻异常体效果较差。     

基于非结构三角网格的海洋CSEM和MT二维联合反演研究

采用非结构三角形网格结合快速Occam算法对海洋可控源电磁(CSEM)与大地电磁(MT)数据进行二维联合反演试验研究。非结构三角形网格能够准确地模拟起伏地形和复杂地质构造,对反演目标区域采用精细网格剖分,其他区域采用粗网格剖分,在满足精度的前提下减少了不必要的计算量。将CSEM与MT数据加入到同一反演数据集中,通过对联合反演数据权重公式进行推导,构建CSEM及MT数据相关权重因子,控制不同数据的拟合权重来实现联合反演。最后对不同模型进行了反演计算,结果表明联合反演比单一反演对海底构造与异常体的还原度更高,验证了联合反演算法的可靠性。     

基于非结构网格的带地形MT二维Occam反演及应用

本文对基于自适应非结构三角形网格的带地形MT数据二维Occam反演进行应用研究.自适应非结构三角形网格能够准确地模拟起伏地形和复杂地质构造,正演网格由有限元解的后验误差估计指导自动细化,保证了模型响应的准确性,反演网格在反演目标区域采用精细网格剖分,在模型边界区域采用粗网格剖分,在满足反演精度的前提下减少了不必要的反演...     

电阻率层析成像及演中若干问题的讨论

以双异性点电流源电流场二维地电断面为例，着重讨论了电阻率展析成像反演中有限单元法数值模拟边界条件的选择、网格剖分以及模拟精度等问题，同时也讨论了成像反演中最优化及演方法及拟合参数的评价等。文中提出的新的网格剖分方法和选择的边界条件可有效地提高计算速度和数值模拟的精度。     

网格尺寸及边界对大地电磁有限元正演精度的影响

大地电磁正演的精度对后续反演意义重大。决定大地电磁有限元正演精度的因素,主要有插值函数、方程组的求解误差,以及网格的剖分程度。而网格剖分的合适与否,为有限元的求解提供了先决条件。因此,从边界条件和横向网格以及纵向网格三方面,讨论了网格剖分对大地电磁正演精度的影响。经研究表明,对于一维和左右对称的二维介质,左右和下边界都无需放置在无穷远处,并且横向网格对精度影响不大。而由于纵向网格与近地表的最小趋肤深度密切相关,所以地表网格不得大于1/4个趋肤深度。     

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

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

基于矢量有限元法的三维CSAMT模拟对近场效应影响因素的研究

可控源音频大地电磁测深法(CSAMT)的人工可控场源会带来近场效应的问题,如果在实际数据反演时利用了近场数据,将会导致对地质构造的错误推断,从而给实际工作带来麻烦和经济损失。本文采用六面体网格剖分的矢量有限元方法进行三维CSAMT正演数值模拟,通过建立不同的地质模型,研究近场效应的影响因素及规律。层状模型的正演结果与理论解对比验证了算法以及程序的正确性。设计了不同电阻率的均匀半空间、场源下存在异常体、场源与接收点之间存在地形和地表异常薄层等地质模型,分析了视电阻率和阻抗相位的响应特征。得到如下结论:(1)均匀大地电阻率大小、场源与接收点之间存在低阻异常体,场源与接收点之间存在山峰地形以及地表异常薄层均会影响近场效应;(2)场源和接收点之间存在高阻异常体,场源与接收点之间存在山谷地形,场源下方存在异常体虽然影响近区的数据,但是对近区频点基本无影响。研究结果可为实际工作中收发距的选取提供参考,具有一定的理论价值。     

起伏地形条件下的磁异常二维聚焦反演及应用

地形对磁异常的反演和解释有很大的影响,特别是磁异常体贴近地表分布时更是如此,为了降低这种影响,笔者推导出了基于有限延深的二度厚板状体的带地形的磁异常正演公式,并提出了一种基于起伏地形的地下网格模型剖分的新方法。同时,聚焦反演对物性体有很强的聚焦效果,如果能够预先判断磁异常主要由地下浅层的强磁性体所引起,则可以考虑采用聚焦反演的方法得到磁性体的位置,显然这种方法非常适合铁矿体的反演。笔者将带地形的二维磁异常正演与聚焦反演结合起来,进行了相应的模型试验,并与水平地形下的聚焦反演及起伏地形下的光滑反演进行对比,验证了这种结合的有效性。最后,将该方法应用于新疆特克斯的磁异常剖面,对该剖面的磁性体分布做出了合理的解释。     

覆盖层下三维板状体地-井瞬变电磁响应

地-井瞬变电磁法是地质找矿的有效方法之一,它是利用地表发射、井中接收瞬变场响应来实现查找深部矿产资源的一种方法。在时间域有限差分算法的基础上,以回线源为激发源,采用非均匀网格剖分技术,对均匀半空间和倾斜板状体、含有低阻覆盖层的倾斜板状体模拟。通过对比分析表明:均匀半空间中,响应极大值随观测时间延迟而减小;场源位于接收钻孔正上方时,异常体响应特征明显;当场源位于旁侧位置,需要较长的观测时间获得异常体的响应;低阻覆盖层对异常响应的影响在整个观测时间范围内都存在,减弱了异常体的响应;覆盖层的影响随发射场源与接收钻孔位置距离的增加而减弱。研究工作为定性分析地-井TEM响应特征和资料解释时提供参考。     

基于非结构化网格的大地电磁2D自适应有限元正演模拟

有限元求解大地电磁正演问题时,对研究区域的剖分常规做法是采用规则化的网格。但规则化网格在剖分地形、断层、褶皱等复杂模型会产生较大的几何离散误差。针对上述情况,这里采用非结构化的四边形网格对二维地电模型进行网格剖分,并与自适应有限元相结合,由剖分的粗网格出发,利用每个频点下网格单元的后验误差估计值指导网格的局部加密,优化网格质量和数量,从而提高正演模拟的精度。通过一维K型地电模型利用本文方法算出数值解与解析解进行对比分析,验证了自适应有限元法在求解大地电磁二维正演模拟中的有效性。并通过对断层、褶皱模型的正演模拟,分析了其大地电磁正演响应特征。     

任意形状三度体磁场可视化人机交互反演

采用任意形状三度体人机交互反演方法解释复杂形态磁性体。其作法是采用直立线元法近似计算地质体产生的磁异常,利用可视化技术实现人机交互反演。通过理论模型分析了影响精度的因素:当剖分剖面数由17减为3时,近似计算的均方误差由0.099nT增大到3.160nT,说明剖分剖面数越多,近似计算的误差越小;当剖分剖面数固定时,地质体中心深度由100m增加到500m,近似计算的最大相对误差从39.19%下降到7.84%,说明中心深度越大,近似计算的误差越小。因此,合理确定剖分剖面数是提高精度的关键。提出了提高近似计算精度的2种措施:模型剖分准则和位场曲线约束模型准则。将该方法应用于澳大利亚艾尔半岛某矿区条带状含铁建造(BIF)的评价,控制钻孔证实了反演解释的结果。     

基于区域Delaunay自动剖分的含水层参数反演

针对地下水有限元数值模拟中区域三角网格剖分复杂难以处理的情况,提出适合其特点的Delaunay三角网格自动剖分方法,并对含有多个参数分区的含水层进行网格剖分,最后利用遗传算法和有限元程序相耦合来反演含水层水文地质参数。结果表明此方法可大大简化地下水数值模拟的前处理工作,并能提高有限元网格剖分的有效性和准确性,从而得到令人满意的数值模拟结果。     

基于OpenMP的重力张量并行正演CSCD

重力张量是重力位的二阶空间导数,对密度体的变化和细节部份反映更为灵敏。对于复杂的重力密度模型来说,张量的解析公式是很难推导的。为了模拟这种情况,将复杂模型进行有限元剖分,计算每个单元的重力异常对测点张量的影响,最后叠加得到整个复杂模型在测点处的重力张量。采用Delaunay四面体非结构化网格对密度体进行剖分,经分析表明,这种剖分方式具有较高的精度。有限元剖分的方式拟合复杂形体缺点之一是计算量比较大,而并行化能有效解决这个问题。这里基于OpenMP并行模型编写了重力张量并行正演程序,分析了不同情况下并行的执行性能,为大规模复杂模型的重力张量正演提供了一种并行策略和思路。     

起伏地表下的复杂三维地质模型建立与重力异常计算

由于地质体和矿体的形态非常复杂，使用长方体网格离散建立正演模型时可能和真实情况有很大差别，因此计算结果可靠性差。本文提出一种基于约束Delaunay网格剖分的方法对地质体进行离散并进行重力建模，在模型边界等复杂区域使用网格自适应加密技术，将三维地质体离散为有限个四面体；并详细推导出针对四面体网格的重力正演公式，实现了基于约束Delaunay网格剖分技术的三维重力数值模拟；最后，针对一个合成数据模型，将计算解与解析解对比。结果表明，细化网格的模拟结果比粗糙网格更好，满足数值模拟的精度要求。将该方法应用到金川矿区实际地质体建模中，根据局部需要，建立各处网格密度不均匀的三维模型，并计算该模型的地表重力场，而后对比模拟数据与实测数据，结果表明Delaunay网格建模方法具有很强的适用性，能够模拟复杂的地质体重力异常。     

深岩溶区浅部结构CSAMT探测解译方法探索 ——以桂林兴安县谢必湾村—上界富坪村段为例

利用CSAMT(可控源音频大地电磁法)对桂林市兴安西部地区的岩溶地质结构进行了探测研究.为了获得更多的岩溶地质信息,通过设定不同深度的反演参数,对比分析了网格剖分大小对反演结果分辨率的影响,以此获得较高分辨率的数据反演方法.利用小网格剖分反演获得的浅部地质结构,比利用大网格反演获得的浅部地质结构的分辨率高,通过小网格可...     

基于OpenMP的重力张量并行正演

重力张量是重力位的二阶空间导数,对密度体的变化和细节部份反映更为灵敏.对于复杂的重力密度模型来说,张量的解析公式是很难推导的.为了模拟这种情况,将复杂模型进行有限元剖分,计算每个单元的重力异常对测点张量的影响,最后叠加得到整个复杂模型在测点处的重力张量.采用Delaunay四面体非结构化网格对密度体进行剖分,经分析表明,这种剖分方式具有较高的精度.有限元剖分的方式拟合复杂形体缺点之一是计算量比较大,而并行化能有效解决这个问题.这里基于OpenMP并行模型编写了重力张量并行正演程序,分析了不同情况下并行的执行性能,为大规模复杂模型的重力张量正演提供了一种并行策略和思路.     

基于积分方程法的大地电磁三维反演

采用积分方程法计算均匀半空间中三维异常体的大地电磁响应,然后利用阻尼最小二乘法反演得到地下三维异常体电阻率分布。在计算灵敏度矩阵时,将灵敏度矩阵分为线性项和体现反演过程非线性的偏微分项进行求解,提高了计算精度。理论模型反演结果表明,反演较准确,对初值依赖性小,且只针对异常体剖分单元进行反演,反演效率较高。     

地—井瞬变电磁响应三维时域有限差分模拟

地—井瞬变电磁法是利用地表发射、井中接收瞬变场响应来实现查找深部矿产资源的一种方法。在时间域有限差分算法的基础上,以回线源为激发源,采用非均匀网格剖分技术,对均匀半空间和倾斜板状体、含有低阻覆盖层的倾斜板状体进行了数值模拟。通过对比分析表明:均匀半空间中,响应极大值随观测时间延迟而减小;早期主要反映围岩响应特征,异常体响应特征在晚期表现出来;场源位于中心位置激发时,异常体响应可以很好地体现出来,当场源位于右侧位置,则要选取较长的观测时间降低周围围岩的影响,使异常体的响应特征显示出来;低阻覆盖层对异常响应的影响在整个观测时间范围内都存在,削弱了异常体的响应特征。研究工作为定性分析地—井TEM响应特征和资料解释提供了参考。     

基于泊松方程的空间波数混合域二度体磁异常数值模拟

正演数值模拟是反演成像的基础。为了实现磁法勘探精细化反演成像与定量解释,本文利用把一个大问题分解为多个小问题的思路,提出一种基于泊松方程的高效、高精度空间波数混合域二度体磁异常数值模拟方法。该方法利用傅里叶变换把二度体磁位偏微分方程转换为一维常微分方程,并采用基于二次插值的一维有限单元法求解该方程,进而通过反傅里叶变换得到空间域磁异常。在模型算例中,分别设计截面为矩形的常磁化率和变磁化率二度体模型,针对本文算法的计算精度和计算效率进行了验证。模型算例结果表明：该算法计算精度高,相对误差绝对值均小于1%;计算速度快,网格节点剖分2 501×2 501的模型模拟时间为4.18 s;适用于任意复杂地形模型。     

二度体的重力张量有限元正演模拟

介绍了重力梯度张量,并将有限单元法应用于二维重力梯度张量的正演计算。为了验证有限元正演方法的精度,对截面为矩形的两个二度体组合模型进行有限元正演模拟,结果表明正演曲线与理论曲线形态一致,拟合情况好。通过对截面形状不规则、密度分块均匀的二度体进行正演模拟,说明有限元法可通过网格剖分来逼近不规则目标体的边界,并对剖分单元赋予不同的密度值来实现对复杂二度体的重力张量的正演模拟。