共查询到20条相似文献，搜索用时 421 毫秒
1.
徐世浙 《物探化探计算技术》1986,(1)
本文给出水平地形下二维地电断面的电阻率及面激发极化异常的数值计算方法。首先用镜象法将半空间问题变为全空间问题,然后用付氏变换将点源二维地电断面的三维电场边值问题变换成二维边值问题。再用格林公式将二维边值问题转变成积分方程。用边界单元法解积分方程,并由付氏反变换,获得三维空间中的电位。用本方法计算了几种模型的电阻率异常及激发极化异常,并与模型实验对比,效果是良好的。与起伏地形二维地电断面的边界单元解法比较,本文所需的计算机内存减少很多,计算时间节省几倍。 相似文献
2.
边界元法(BEM)是一种新的数值方法。由于该法可以降低所研究场问题的维数,因而边界元法较之域型法(FEM和FDM)具有应用简便、数据量少、计算快、精度高等优点。从而,用该法解决了域型法难以实现的三维地电模型视电阻率异常的计算问题。本文论述了用边界无法求解点源场地表水平和起伏下三维地电体位场问题的方法原理和数值处理技术,并给出了若干算例:导电球状矿体上视电阻率数值解与解析解结果;导电球状矿体上不同测线上视电阻率平剖曲线;三维山脊地形下导电椭球状矿体上视电阻率联剖曲线及其地形改正结果等。由本文内视电阻率的边界元法数值解与解析解结果对比的一致性和三维地形校正的显著效果,表明了边界元法是求解任意三维地电模型上位场问题的有效方法。由于用该法实现了对任意三线地质体上空间位场计算,必将推动三维电法勘探工作的开展与深入研究,无疑会对提高电法勘探的地质效果发挥重要作用。 相似文献
3.
三维地电断面激发极化法有限元数值模拟 总被引：5，自引：0，他引：5
用有限单元法求解三维地电断面激发极化法正演模拟算法.首先给出了三维构造中点源电场异常电位的边值问题与变分问题, 简化无穷远边界上的边界条件以提高计算速度及计算精度.以此为基础计算视电阻率对模型电阻率的偏导数矩阵, 并进行三维地电断面激发极化法正演模拟, 与等效电阻率法相比节约了计算时间.对几例较典型的地电模型进行计算, 结果表明本方法是正确可行的. 相似文献
4.
本文以二层介质中存在球体的模型为例，讨论了用边界元方法解决层状介质中存在三维不均匀体情况下点源场视电阻率的计算问题。文中通过对基本解的选取，避免了对水平界面的积分计算，并且在积分方程中消去了位的法向导数，使得积分方程的形式更为简单，也大大地减少了所需的边界剖分节点数和未知量的个数，节省了计算机内存，提高了计算精度。计算的结果表明，边界元法对本问题的研究是成功的。 相似文献
5.
变步长激发极化法有限元数值模拟 总被引：1，自引：0，他引：1
用改变剖分步长的方法对三维地电模型进行激发极化法正演模拟计算.首先,在给出三维地电模型异常电位的边值问题与变分问题的基础上,通过对三维激电模型进行变步长的六面体剖分,将异常体部分用较短步长的单元进行加密剖分,并在各个单元上进行插值,于是得到其总刚度矩阵;然后用变带宽方式对其进行编号存储,并选用成熟的LDLT方法求解大型稀疏矩阵,运用等效视电阻率方法求取地表视极化率值.通过用中间梯度装置对几例典型的激电模型进行了正演计算,结果表明,该方法能够达到满意的精度和效果,从而为激发极化正演解释提供有益的帮助. 相似文献
6.
基于非结构化网格的2．5D直流电阻率模拟 总被引：2，自引：1，他引：1
结合大型通用有限元软件ANSYS所提供的强大网格剖分功能,实现了2D地质模型的非结构化三角形网格剖分,并编写了2.5D直流电阻率有限元法FORTRAN语言计算程序.经计算对比表明,在节点数基本相同时,非结构化三角形网格比传统的矩形网格、矩形～三角形网格具有更高的计算精度,可以很好地拟合复杂地形和地质体边界.在此基础上,采用非结构化三角形网格,计算了起伏地形下2.5D地电模型的视电阻率异常,并利用比较法进行了地形改正.与水平地形时的结果对比表明,比较法可以较好地消除地形影响,突出局部地质体的异常. 相似文献
7.
点源二维起伏界面二层介质边界元解 总被引：3，自引：0，他引：3
本文以二层介质下界面上隆模型为例,讨论了用边界元法在解决层状介质中存在二维异常体情况下点源场视电阻率计算问题。文中对基本解的选取避免了对水平层状界面的积分计算,从而大大减小了所需的边界剖分点和未知量的数目,大量地节省了计算机内存,改善了计算精度,使其从解决二层介质中横向不均匀性问题推广到解决三层以上介质横向不均匀性的问题。结果表明,边界元法对本问题的研究是成功的。 相似文献
8.
复杂地形下高密度激电法2.5维有限单元法数值模拟 总被引：1，自引：0，他引：1
目前三维有限元模拟计算量大,计算效率低,对计算机要求较高,2.5维模拟是三维问题的简化,较好地克服了上述问题。在高密度激电法模拟中,模拟的点数一般较大,要求程序计算效率高,因此对高密度激电法2.5维模拟研究是必要的。首先给出了2.5维稳定电流场的边值问题及对应的变分问题;为了能更好模拟复杂地形对极化率异常的影响,采用三角形网格,有利于模拟复杂地形又有利于将其应用于反演计算。推导了基于连续电性介质的2.5维稳定电流场的有限元法,利用等效电阻率法,编制了2.5维高密度激电法有限元模拟程序,计算了水平地形及起伏地形下的极化率异常模型,并验证了方法的可行性,分析得出了山谷地形比山脊地形对极化率观测的影响大的结果。 相似文献
9.
本文重点阐述了点源二维电阻率测深边界单元方法原理、方法技术。针对电阻率测深装置特点。提出的区间截断剖分法,有利于节省计算机内存,提高计算速度,降低成本。文中用边界元法计算了若干地电形体上对称四极测深和等比装置的ρs测深曲线,获得了一些有益的结果。 相似文献
10.
COMSOL是一个多物理场仿真软件,笔者尝试将其应用于大地电磁法的数值模拟,研究起伏地形对大地电磁响应的影响及解决方案。本文由麦克斯韦方程组出发,施加第一类边界条件,将地形起伏数字化为参数化曲面,导入COMSOL进行计算。在对COMSOL的精度进行了验证的基础上,分析不同频点所需内存以及自由度,然后选择合适的频点进行研究。将三维起伏地形的模拟结果与无地形起伏时的结果进行对比,总结出山谷、山峰地形对三维大地电磁响应的影响,发现在XY模式下,地形对视电阻率的影响大于相位,在视电阻率等值线图中能够清晰地看出起伏地形的轮廓,而在YX模式下,地形对相位的影响大于对视电阻率的影响。最后,利用“比值法”进行地形校正,校正前由于起伏地形影响,在异常体附近会出现部分视电阻率数值跳跃现象,通过校正有效抑制了起伏地形对异常体附近视电阻率的影响,校正后的视电阻率曲线整体上移,并且异常体周围的视电阻率值也趋于背景场值,获得了较好的效果。 相似文献
11.
三维地形大地电磁场的边界元模拟方法 总被引：1，自引：0，他引：1
提出了一种用边界元法计算大地电磁场三维地形影响的数值模拟方法.首先用矢量积分理论和电磁场边界条件, 将上半空间(空气)和下半空间(地下介质)两个区域电磁场边值问题变为仅对地形界面的两个矢量面积分方程, 其中一个计算磁场, 称磁场方程; 另一个计算电场, 称电场方程.然后将对地形界面的积分剖分为一系列的三角单元积分.在三角单元积分中, 假设单元中电磁场为水平均匀大地空间电磁场与地形影响的迭加, 并假设地形影响为常项, 这样既保证了计算精度又使得计算方法简便.通过分解和计算, 每一个矢量面积分方程分解为对应3个坐标方向的3个常量线性方程, 这些线性方程组成了对角占优的线性方程组, 可用SSOR方法求解.文中给出了2个三维地形上大地电磁视电阻率曲线的计算结果. 相似文献
12.
复杂地形条件下激发极化有限单元法三维数值模拟 总被引：2，自引：1，他引：1
野外激电勘查过程中,据激电异常进行钻探验证时,常发生验证目标矿体擦边而过现象.究其原因是在实际勘查中,地形的起伏变化改变了地下电场的分布规律,导致激电异常发生位移及幅值改变.为了改善、克服这种现象,必须从理论上进行研究.运用四面体单元三维有限单元数值正演计算方法对此进行了理论研究,选用典型山谷地形水平层状极化体模型进行数值模拟计算.研究结果表明:①起伏地形不产生假的激电异常,仅对激电异常的形态和幅值产生影响;②采用四面体剖分单元有限单元法研究复杂地形条件对激电异常的影响是行之有效的. 相似文献
13.
讨论了小型脉状金矿的地质及激发极化异常特征,通过已知脉状金矿体的实例,说明寻找隐伏金矿时,要注意研究低缓的,幅度较小但稳定连续的异常,这些激发极化异常有可能是小型脉状金矿体的反映。由于成矿条件的限制,在金属硫化物含量不同的情况下,矿体可以表现为低阻高极化体,也可能表现为高阻高极化体。电阻率法受地形影响大,不易分辨小型脉状体的异常,可以作为评价和解释的参考信息,而激发极化法受地形影响小,是寻找含硫化物的小型脉状金矿体的有效方法。 相似文献
14.
15.
采用易于模拟复杂地形起伏和倾斜界面的非结构三角单元剖分网格,并利用对偶加权后验误差估计指导网格自动细化过程,实现了2.5维直流电阻率法自适应有限元数值模拟。在实例模型分析中,分别计算了层状模型和垂直岩脉模型的直流电阻率响应,并与其解析解进行了比较。对比结果表明,该算法所得数值解精度很高,解的相对误差小于0.5%。最后,计算了起伏地形2.5维地电模型视电阻率异常,并利用比较法进行了地形改正。地形改正结果与水平地形时的结果对比表明,比较法可以较好地消除地形影响,突出局部地质体异常。 相似文献
16.
一般情况下,边界元法所建立的线性方程组系数矩阵为一满置矩阵。针对地形剖面的计算特点,本文提出了满置矩阵转换为带宽矩阵,用一维压缩存贮该带宽矩阵的方法,编制了点源二维地形视电阻率响应的计算程序,极大地减少了计算机内存和时间。 相似文献
17.
Zohdy方法在三维电阻率测深数据反演中的应用 总被引：2，自引：0，他引：2
三维电阻率测深数据采用Zohdy方法近似反演,正演计算采用有限单元法,反演初始模型由测量视电阻率数据给出,通过比较实测视电阻率值和预测模型计算的视电阻率值对数差来修改模型网格电阻率.为了解决任意电极距测深数据的反演,采用大、小双网格剖分,大网格反映地下电性分布情况,小网格用于实际有限元正演计算.文章对加5%随机噪声的模型理论视电阻率测深断面数据和实测视电阻率测深断面数据进行反演.结果表明该反演方法需要的内存小、计算量少,并能够取得很好的反演效果. 相似文献
18.
19.
三维地层的三棱柱剖分与土方计算 总被引：5，自引：0，他引：5
三维地层剖分算法是实现地层三维可视化与土方计算的重要基础。采用三棱柱体建模方法构造三维地层,分析任意平面切割三棱柱的不同情形,得出了针对三棱柱地层的剖分算法;提出了以切割面为分界面,将单个三棱柱剖分成多个小三棱柱的三角投影分割法--将三棱柱顶面、底面及切割截面按一定规则进行合理三角化划分,使划分后三角形在水平面上的投影相重合;将该3层三角形组合成以切割面为分界面的2层小三棱柱,从而完成剖分。土方计算中首先给出开挖体的边界控制方程,根据该控制方程所限定的内边界得到开挖区域,结合上述剖分算法,对该区域内所有小三棱柱体积求和,得到开挖土方量。将以上算法应用到广东省某重大水利工程中,实现了基坑模拟开挖及其土方计算。 相似文献
20.
起伏地形对可控源音频大地电磁(CSAMT)响应具有强烈的影响,因此在CSAMT数据处理解释时需要考虑地形。同时,实际的地下地质情况和地表的地形情况通常比较复杂,地质结构和地形大部分情况下都是三维的。在水平地表三维有限差分CSAMT数值模拟算法的基础上,推导了利用地下交错网格采样点处的总磁场计算起伏地形下空气-地下介质分界面处的总电场和总磁场的表达式,从而实现了起伏地形下三维CSAMT数值模拟算法。在算法实现过程中,采用伪δ函数代替麦克斯韦方程中的场源项和直接计算总场的策略,避免了原有的将总场分离成背景场和二次场的策略在复杂地质条件下难以选择合适背景电阻率的问题。为了直接模拟总场,起伏地形下三维CSAMT数值模拟算法给出了新的三维正演方程的边界条件。将模拟水平地表三维异常体和三维山峰地形两个理论模型得到的响应结果与前人算法的计算结果进行对比,验证了所实现算法的正确性和有效性。 相似文献