共查询到18条相似文献，搜索用时 125 毫秒
1.
基于非结构化网格的2．5D直流电阻率模拟 总被引：2，自引：1，他引：1
结合大型通用有限元软件ANSYS所提供的强大网格剖分功能,实现了2D地质模型的非结构化三角形网格剖分,并编写了2.5D直流电阻率有限元法FORTRAN语言计算程序.经计算对比表明,在节点数基本相同时,非结构化三角形网格比传统的矩形网格、矩形～三角形网格具有更高的计算精度,可以很好地拟合复杂地形和地质体边界.在此基础上,采用非结构化三角形网格,计算了起伏地形下2.5D地电模型的视电阻率异常,并利用比较法进行了地形改正.与水平地形时的结果对比表明,比较法可以较好地消除地形影响,突出局部地质体的异常. 相似文献
2.
角域是描绘起伏地形的基本单元,通过角域叠加可以获得复杂纯地形产生的视电阻率异常,进而可以进行地形校正。将简单角域视电阻率畸变进行组合叠加得到复杂角域地形的视电阻率畸变,分析了组合叠加误差,并与ANSYS软件数值模拟结果进行了对比。结果表明,在地形起伏不大的情况下,可以通过组合叠加获得复杂角域地形的视电阻率异常。最后给出了利用组合叠加法获得纯地形影响后应用比较法进行地形改正的模型算例;与水平地形时的结果对比表明,比较法可以削弱地形影响,突出局部地质体的异常。 相似文献
3.
为了解起伏地表对复电阻率法数据的影响,发展了起伏地表条件下复电阻率法2.5维有限元数值模拟方法。经数值模拟结果与解析解对比表明,该模拟方法在计算精度上得到了较满意的结果;几种常见纯地形异常的数值结果表明:起伏地形对复电阻率法数据的畸变特征与地形坡度大小以及地形分布范围有很强的相关性,坡度越大,地形影响越为剧烈;坡度角为45°的山脊地形、山谷地形和斜坡地形所引起的复电阻率法实分量异常分别达到背景值的21.2倍、3.9倍和3.4倍,虚分量异常分别达到背景值的2.3倍、2.2倍和2.1倍,振幅异常分别达到背景值的21.2倍、3.9倍和3.4倍;而相位受地形影响非常微弱。这些工作为复电阻率法的地形改正和反演奠定了基础。 相似文献
4.
复杂地形下高密度激电法2.5维有限单元法数值模拟 总被引：1，自引：0，他引：1
目前三维有限元模拟计算量大,计算效率低,对计算机要求较高,2.5维模拟是三维问题的简化,较好地克服了上述问题。在高密度激电法模拟中,模拟的点数一般较大,要求程序计算效率高,因此对高密度激电法2.5维模拟研究是必要的。首先给出了2.5维稳定电流场的边值问题及对应的变分问题;为了能更好模拟复杂地形对极化率异常的影响,采用三角形网格,有利于模拟复杂地形又有利于将其应用于反演计算。推导了基于连续电性介质的2.5维稳定电流场的有限元法,利用等效电阻率法,编制了2.5维高密度激电法有限元模拟程序,计算了水平地形及起伏地形下的极化率异常模型,并验证了方法的可行性,分析得出了山谷地形比山脊地形对极化率观测的影响大的结果。 相似文献
5.
《物探化探计算技术》2019,(6)
为了提高二维直流电阻率法的数值模拟精度,这里实现了2.5维各向异性介质点源直流电场异常场的数值模拟。首先推导出2.5维各向异性介质总场和异常场满足的变分问题;然后采用矩形内三角网格的方式进行单元剖分,三角单元内二次场值和异常电阻率均设计为线性变化。为了简化背景场的计算,将背景介质的电性主轴与观测坐标系的夹角设计为零,从而避免了直接计算任意各向异性介质电位的傅里叶变换步骤。同时对计算中点源波数的选择和适用情况进行了讨论。最后通过对简单层状模型的计算,与解析解进行对比,验证了本数值算法的正确性和可靠性。比较常规各向同性介质模拟方式,本数值模拟发现直流电阻率法对横向电阻率的变化更为明显。 相似文献
6.
为了模拟实际的起伏地形,在前人研究的基础上,单元网格设计为任意四边形网格,单元内的场值双线性变化,采用了高斯数值积分法计算单元系数矩阵,并给出二维大地电场的辅助场表达式。通过模型算例表明,模拟结果与解析解的均方误差在1%以内,地形模型与前人的模拟结果相吻合。分析对比了两个不同坡度模型对视电阻率和相位的影响。采用任意四边形网格剖分法,降低了编程难度,可以方便地适应野外地形的起伏变化。 相似文献
7.
边界元法(BEM)是一种新的数值方法。由于该法可以降低所研究场问题的维数,因而边界元法较之域型法(FEM和FDM)具有应用简便、数据量少、计算快、精度高等优点。从而,用该法解决了域型法难以实现的三维地电模型视电阻率异常的计算问题。本文论述了用边界无法求解点源场地表水平和起伏下三维地电体位场问题的方法原理和数值处理技术,并给出了若干算例:导电球状矿体上视电阻率数值解与解析解结果;导电球状矿体上不同测线上视电阻率平剖曲线;三维山脊地形下导电椭球状矿体上视电阻率联剖曲线及其地形改正结果等。由本文内视电阻率的边界元法数值解与解析解结果对比的一致性和三维地形校正的显著效果,表明了边界元法是求解任意三维地电模型上位场问题的有效方法。由于用该法实现了对任意三线地质体上空间位场计算,必将推动三维电法勘探工作的开展与深入研究,无疑会对提高电法勘探的地质效果发挥重要作用。 相似文献
8.
基于二次场二维起伏地形MT有限元数值模拟 总被引：4，自引：1，他引：3
通过计算二次场来进行二维大地电磁数值模拟;导出了二维大地电磁二次场的微分方程,利用有限单元法来解微分方程;对矩形网格进行对角线的二次剖分,更容易且真实地模拟起伏地形。对几个典型模型进行了试算,与前人总场法的计算结果做了比较,两者视电阻率曲线一致,证明本文算法是正确的;通过2个简单的算例说明复杂地表下2种极化模式的MT观测资料都有明显的异常,视电阻率在TM模式下比TE模式更易受地形影响,TE模式下视电阻率曲线形态与地形呈"正相关",TM模式下反之。 相似文献
9.
《物探化探计算技术》2017,(1)
利用交错网格有限差分技术实现了频率域海洋可控源电磁三维正演算法,并探讨了海底地形对海洋可控源电磁响应的影响和校正方法及效果。通过与海底二维山峰地形的二维自适应有限元解对比,验证了文中算法的正确性和模拟海底起伏地形的可行性,并分析了海底二维山峰地形对三维油气储层电场响应的影响。关于海底起伏地形对海洋可控源电场响应的影响,采用比较法对其进行了地形校正,结果表明,在深水区比较法可以很好地消除地形影响,恢复海底高阻油气储层的电场响应。 相似文献
10.
徐世浙 《物探化探计算技术》1986,(1)
本文给出水平地形下二维地电断面的电阻率及面激发极化异常的数值计算方法。首先用镜象法将半空间问题变为全空间问题,然后用付氏变换将点源二维地电断面的三维电场边值问题变换成二维边值问题。再用格林公式将二维边值问题转变成积分方程。用边界单元法解积分方程,并由付氏反变换,获得三维空间中的电位。用本方法计算了几种模型的电阻率异常及激发极化异常,并与模型实验对比,效果是良好的。与起伏地形二维地电断面的边界单元解法比较,本文所需的计算机内存减少很多,计算时间节省几倍。 相似文献
11.
付德俊 《物探化探计算技术》2019,(3):365-370
开展了不同观测方式的井地2.5D直流电阻率反演研究。①从2.5维直流电阻率满足的边值问题出发,采用变分原理结合节点线性基函数推导了2.5维井地DC满足的积分弱解形式;②构建了二阶最大平滑稳定泛函的2.5维井地直流电阻率正则化目标函数,采用共轭梯度算法对正则化目标函数进行最优化求解,并采用逐步衰减正则化因子的求解策略来提高反演的稳定性;③设计了均匀半空间模型得到的数值解与解析解对比,电位的相对误差在2%以内,阐述正演算法的正确性和高精度。另外,分析了不同观测装置的2.5维井地直流电阻率异常体特征,并对不同观测方式对2.5维井地DC理论数据进行反演研究。研究结果表明,井中数据的引入提高了2.5DC对纵向探测的分辨率能力,同时提高了2.5DC反演有效性以及准确性。 相似文献
12.
2.5维起伏地表条件下坐标变换法直流电场数值模拟 总被引：3，自引：1，他引：2
研究起伏地表对视电阻率分布的影响是进行地形校正的基础。由于很难处理不规则边界,计算简洁且效率高的有限差分方法很少用于解决起伏地表问题。为了解决该问题,引入曲化平思想,实现了一种基于坐标变换法的起伏地表条件下的直流电场数值模拟方法。方法从传统的2.5维基本方程及边界条件出发,通过坐标变换将起伏地表问题转化为水平地表问题,并利用有限差分法求解水平地表问题,最后再通过坐标映射得到起伏地表条件下的电位和视电阻率分布。精度分析及计算实例表明：本方法实现简洁,计算误差主要集中在震源附近,整个计算区域内的平均相对误差为1.39%,计算结果满足起伏地表条件下的电位和视电阻率的分布规律。 相似文献
13.
针对巷道直流电阻率法超前探测的三维数值模拟,应用代数多重网格快速算法对二次场的有限差分问题进行求解。将均匀全空间板状体数值结果与解析解对比,最大误差不超过0.28%。模拟结果表明,巷道空腔对测量结果存在一定影响,采用比值曲线消除该影响,修正后的视电阻率极值与理论值误差为2.4%。进一步研究受旁侧影响的前方异常识别方法,可通过在不同巷道面进行定点源测量或掘进过程中的定点测量来识别旁侧影响,达到准确探测前方异常的目的。 相似文献
14.
为了研究复杂地电模型的航空瞬变电磁法全波形响应特征,需要开发考虑发射波形的三维数值模拟算法。本研究基于非结构四面体网格和位移逆Krylov子空间(Shift-and-Invert Krylov,简称SAI Krylov)方法,采用基于电偶极子离散的场源处理方法模拟场源,在时间域进行计算实现了全波形航空瞬变电磁法矢量有限元三维数值模拟。使用均匀半空间模型在阶跃波、半正弦波、三角波和梯形波激发下的全波形解析解、VTEM实际激发波形的后推欧拉算法计算结果,检验了本研究开发的数值模拟算法的正确性。设计地表起伏异常体模型,计算和分析了航空瞬变电磁响应特征。开发的基于位移逆Krylov子空间的全波形航空瞬变电磁法三维数值模拟算法适合模拟复杂地电模型的响应,具有较高的计算精度。 相似文献
15.
高精度、高效率的深海勘探是当前国家海洋资源勘探开发的研究热点之一,设计观测效率高、勘探成本低的探测方式,开发可精确计算海底复杂电性环境的正演模拟算法有助于推进该研究的进展。本文基于前人的研究经验,引入了全拖曳式深海直流电阻率观测方式,同时考虑到海底沉积环境的电各向异性,开发了基于非结构有限元方法的三维任意各向异性深海直流电阻率正演算法,实现了对海底任意电各向异性情况的仿真模拟。本文算法与层状一维模型半解析解的对比验证了算法的精度。对典型各向异性电性模型的仿真模拟和分析,证明了全拖曳式深海探测方式对海底沉积层和矿产资源电各向异性的高分辨率。在此基础上,总结了电导率各向异性主轴沿笛卡尔坐标旋转时视电阻率分布也发生同方向旋转的视电阻率分布特征,并根据海底地形对各向异性探测的影响进行了简单探讨。 相似文献
16.
COMSOL是一个多物理场仿真软件,笔者尝试将其应用于大地电磁法的数值模拟,研究起伏地形对大地电磁响应的影响及解决方案。本文由麦克斯韦方程组出发,施加第一类边界条件,将地形起伏数字化为参数化曲面,导入COMSOL进行计算。在对COMSOL的精度进行了验证的基础上,分析不同频点所需内存以及自由度,然后选择合适的频点进行研究。将三维起伏地形的模拟结果与无地形起伏时的结果进行对比,总结出山谷、山峰地形对三维大地电磁响应的影响,发现在XY模式下,地形对视电阻率的影响大于相位,在视电阻率等值线图中能够清晰地看出起伏地形的轮廓,而在YX模式下,地形对相位的影响大于对视电阻率的影响。最后,利用“比值法”进行地形校正,校正前由于起伏地形影响,在异常体附近会出现部分视电阻率数值跳跃现象,通过校正有效抑制了起伏地形对异常体附近视电阻率的影响,校正后的视电阻率曲线整体上移,并且异常体周围的视电阻率值也趋于背景场值,获得了较好的效果。 相似文献
17.
三维地形大地电磁场的边界元模拟方法 总被引：1，自引：0，他引：1
提出了一种用边界元法计算大地电磁场三维地形影响的数值模拟方法.首先用矢量积分理论和电磁场边界条件, 将上半空间(空气)和下半空间(地下介质)两个区域电磁场边值问题变为仅对地形界面的两个矢量面积分方程, 其中一个计算磁场, 称磁场方程; 另一个计算电场, 称电场方程.然后将对地形界面的积分剖分为一系列的三角单元积分.在三角单元积分中, 假设单元中电磁场为水平均匀大地空间电磁场与地形影响的迭加, 并假设地形影响为常项, 这样既保证了计算精度又使得计算方法简便.通过分解和计算, 每一个矢量面积分方程分解为对应3个坐标方向的3个常量线性方程, 这些线性方程组成了对角占优的线性方程组, 可用SSOR方法求解.文中给出了2个三维地形上大地电磁视电阻率曲线的计算结果. 相似文献
18.
By using the numerical method to model the ter-rain effect on the magnetotelluric field,few resultshave been obtained. The finite element method(FEM) was used by Chouteau and Bouchard (1988)and Wannamaker et al .(1986) ,andthe boundary el-ement method (BEM) was used by Xu and Zhou(1997) and Xu (1995) to model 2Dtopographyinflu-ences on magnetotelluric surveys . The BEM methodwas also used to model the 3Dtopographic effect onmagnetotelluric deep sounding (Xu et al .,1997 ;Xu,1995) .In t… 相似文献