排序方式: 共有21条查询结果,搜索用时 250 毫秒
1.
全局弱式无单元法是在有限单元法基础上发展起来的一种数值模拟方法,它采用局部支持域内的节点信息来构造形函数实现局部精确逼近,摆脱了单元,仅依赖于节点信息,具有预处理简单、模拟精度高、灵活性强的特点,适用于复杂地电条件下直流电阻率正演模拟.本文采用RPIM构造直流电阻率全局弱式无单元法形函数,利用RPIM形函数推导了直流电阻率全局弱式无单元法方程.然后,编制了直流电阻率全局弱式无单元法正演模拟Fortran程序,利用该程序对典型的地电模型进行了正演模拟,并将正演结果与基于线性插值的FEM正演结果及解析解进行对比,结果表明采用RPIM形函数的全局弱式无单元法用于直流电阻率正演模拟的正确性及有效性,且在同等条件下,全局弱式无单元法模拟精度高于矩形剖分的FEM,更有利于指导电法勘探的数据解译;利用该程序对复杂地电模型进行了正演模拟,结果表明全局弱式无单元法对复杂地电模型模拟效果良好,适应性强,灵活性高,可任意加密节点提高模拟精度.
相似文献2.
3.
复杂地形下高密度激电法2.5维有限单元法数值模拟 总被引:1,自引:0,他引:1
目前三维有限元模拟计算量大,计算效率低,对计算机要求较高,2.5维模拟是三维问题的简化,较好地克服了上述问题。在高密度激电法模拟中,模拟的点数一般较大,要求程序计算效率高,因此对高密度激电法2.5维模拟研究是必要的。首先给出了2.5维稳定电流场的边值问题及对应的变分问题;为了能更好模拟复杂地形对极化率异常的影响,采用三角形网格,有利于模拟复杂地形又有利于将其应用于反演计算。推导了基于连续电性介质的2.5维稳定电流场的有限元法,利用等效电阻率法,编制了2.5维高密度激电法有限元模拟程序,计算了水平地形及起伏地形下的极化率异常模型,并验证了方法的可行性,分析得出了山谷地形比山脊地形对极化率观测的影响大的结果。 相似文献
4.
Parker公式在起伏地层重力场正演方面以其简洁、快速等优点被广泛应用于地形改正以及界面反演等,但传统的Parker公式正演方法存在一定的缺陷,在起伏地层模型重力场的正演计算中,由于正演式中e的指数项数值稳定性较差,导致正演结果精度不高。在前人研究的基础上,对Parker公式进行了改进,通过增加上、下界面平均值的方法,提高正演的数值稳定性,改进后的正演精度得到大幅提升。在此基础上,提出了一种曲面观测算法,将观测面由传统的水平面推广到任意起伏面,实现了上下界面起伏地层在起伏观测面上的快速高精度重力正演计算和带地形的密度界面起伏模型的快速高精度重力正演模拟。为了保障数值精度,采用了Gauss-FFT算法,模型实验均取得很好效果。 相似文献
5.
基于有限单元法的二维/三维大地电磁正演模拟策略 总被引:1,自引:0,他引:1
对于二维和三维大地电磁正演问题,有限单元法最后形成了一个线性方程组KX=p。方程组中的K是大型稀疏的带状对称复系数矩阵,其条件数远大于1,为严重病态矩阵,求解其对应方程组会遇到很多困难。不完全LU分解处理的BICGSTAB算法,可用于该线性方程组的求解,并且具有速度快,精度高,稳定性好等优点。为了模拟无穷远边界及满足计算机的内存需求,在保证计算精度的情况下,设计了非均匀网格剖分。在程序编制中,因只存储有限元系数矩阵的非零元素,大大减少了正演计算的时间。通过对二维模型和三维模型电磁响应的计算,验证了该算法的正确性。 相似文献
6.
复杂二维/三维大地电磁的有限单元法正演模拟策略 总被引:1,自引:0,他引:1
复杂二维和三维大地电磁模型的正演数值模拟具有一定的挑战性。对于复杂的二维和三维大地电磁正演问题,我们采用有限单元法进行求解。有限单元法最后形成一个线性方程组,系数矩阵是大型稀疏的带状对称复系数矩阵,并且其条件数远大于1,为严重病态矩阵,求解其对应方程组会遇到很多困难。不完全LU分解处理的Bi-CGSTAB迭代方法可用于该线性方程组的求解,并且具有速度快、精度高和稳定性好等优点;为了模拟无穷远边界及满足计算机的内存需求,在保证计算精度的情况下设计了非均匀网格剖分;在程序编制中,只存储有限元系数矩阵的非零元素,大大减少了正演计算的时间。通过对二维和三维模型电磁响应的计算,验证了算法的正确性。 相似文献
7.
根据青海多龙恰柔勘查区已有的地质信息,以含辉钼矿化、黄铜矿化的印支期、燕山期花岗岩为主要研究对象,用地球化学土壤(岩屑)测量、地球物理双频激电扫面综合圈定了3个有意义的成矿远景区,用激电测深与EH4大地电磁测深初步固定了2个大的异常带的空间大致分布:均是走向以NNW向为主、倾向以SWW向为主、倾角较陡的大型异常带,总结了此矿区地质一地球物理一地球化学成矿模型,最后在推断出的异常部位布设钻孔并见矿。这次工作表明:因地制宜地结合各种物探方法,利用地质填图资料,综合全局考虑化探异常,扬长避短,发挥各自优势,在西部特殊景观区寻找金属矿科研与实践中具有重要意义。 相似文献
8.
作为一种新型开发的绿色能源,地热资源被越来越多的人所重视.开采埋藏较深的地热资源风险大,因此开发前的地质和地球物理勘查是十分必要的.采用单一的地球物理方法勘探具有很大的风险,因此采用多种方法进行综合调查可以降低单一方法的风险,取得较好的效果.本文以狮子湖温泉为例,研究了地球物理方法在温泉勘探的应用.本次勘探采用可控源音频大地电磁法(CSAMT-Controlled Source Audio-frequency Magnetotellurics)和微重力测量,CSAMT是针对大地电磁测深法场源的随机性和信号微弱,提出的一种采用可以控制的人工场源改进方案.仪器采用美国Zonge公司生产的GDP-32Ⅱ.该方法由人工向地下供入音频谐变电流建立电磁场,通过仪器在地面接收从地下反馈来的信息,根据不同时代、岩性地层电性特征达到勘查目的.为此我们在测区做了两条剖面,从CSAMT反演图推断,自上而下可分为3个电阻率层,该剖面视电阻率具有很好的层理特征,反映了厚大的中新生代地层(Edn、K)覆盖.剖面西部有明显的泥盆系地层(D)存在.根据视电阻率的变化特点,可以推断这条剖面的5条断层.在剖面中部距地表400~800米深处存在一明显低阻区,推测应为含水破碎带或低阻泥岩.对数据进行二维反演,可以清楚的看到利于储水的盆地构造.CSAMT方法受静态效应影响很大.静态效应位移可能是由地形和电阻率的浅部的横向变化引起的,既是不可避免的,也是不可预测的.因此对数据进行预处理是十分必要的.减少静态效应的影响办法有以下三种:(1)、对效应进行理论计算;(2)、采用空间滤波和相位积分等处理方法;(3)、使用独立的、无静态效应的测量方法.计算静态效应理论值在理论上是简捷的,但在实际的野外条件下,由于无法预测引起静态效应的物体的几何尺寸和电性参数,因此这种方法无法得到可靠的校正值.空间滤波处理是目前广泛采用的一类方法,Bostick(1986)提出了消除MT数据中静态效应的电磁列阵剖面法(EMAP).EMAP法由于采用连续的剖面测量,可采用窗口可变的自适应空间滤波器-汉宁窗(Hanning window)或叫余弦钟形滤波器消除静态效应.但是这种方法提供了静态效应的要求的数据密度,这就增大了大量的额外测量,提高了获得数据的代价.因此,我们做重力勘探与CSAMT相对比.微重力测量使用美国产LCR-D型重力仪.实测的微重力异常是地下由浅到深各类地质体的物性差异在地面综合叠加的效应,其中包括界面起伏、岩性不均匀、地壳与壳下物质的厚度变化等诸多地质因素在内.实测的重力异常值经过固体潮改正、零点漂移改正、布格改正、正常场改正之后,得到改正后的重力异常值.从微重力反演结果看来,自西往东重力异常逐渐减小,程台阶下降趋势,并趋于平缓,验证了CSAMT的异常结果.从而弥补了CSAMT法的不足.通过这两种方法,我们大致查明了新生代红层盆地的产出形态.该红层盆地西侧边缘位于青山小学一带,自西往东变深;并反演计算出新生代红层与晚古生代泥盆系地层的分界面.深部地球物理勘探的方法有很多,各有各自的优缺点,我们不能从单一的一种方法得到的结论来判断地下地质构造.CSAMT法勘探深度大,但由于本身的物理特性,导致静态效应、近场效应等影响甚大.严重影响我们对地下地质目标体深度的判断,重力和CSAMT法的相互验证,很好的说明了多种地球物理方法综合测定的优势. 相似文献
9.
大地电磁法正演中多重网格法求解的广义傅里叶谱分析 总被引:1,自引:0,他引:1
为了准确预测和分析多重网格法用于大地电磁法正演计算的收敛效果,对多重网格法的收敛性进行了广义傅里叶谱分析。通常情况下,系数矩阵的傅里叶谱是复数,为了直观地判断、提取收敛信息,将谱转换到实数域。在实数域内,特征向量谱定量解释了二重网格法收敛慢的原因(最粗网格用Gauss-Seidel法求解)。传统的局部傅里叶谱分析没有考虑边界条件和模型参数的变化,在分析二重网格法求解收敛性时得出的渐进收敛估计与数值解偏差大。针对这一问题提出广义傅里叶谱分析,其结果与数值解接近(比如2-V(0,1)的广义傅里叶谱分析为0.706,真实值为0.710)。对五重网格法求解的高重广义傅里叶谱分析结果表明,随广义傅里叶谱分析分析重数的增加, 所求渐近收敛估计趋于收敛的数值结果。基于此得出高重广义傅里叶谱分析的经验公式。求得低重广义傅里叶谱分析,通过矫正近似得到高重广义傅里叶谱分析,其结果在本文2个例子有效。 相似文献
10.
从电磁正演模拟结果出发,引入了地震勘探中“照明度”的概念,提出了电磁勘探的照明度定义及特征函数等。对照明度的地表平面特征、随深度变化特征和随频率与大地电阻率的变化特征进行了模拟分析,总结了其规律。通过分析,对广域电磁法的观测系统提出了优化建议,对野外测量的偏移角度和范围提供了借鉴。 相似文献