首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 109 毫秒
1.
直流电法正演计算中的有限单元法(第三讲)   总被引:1,自引:1,他引:1  
为了解决在复杂地电条件下电法勘探中场的计算,近十年来,国内外发展了电法勘探正演计算的数值模拟方法,或者称为计算机模拟,其中一个比较活跃的领域就是有限单元法。电法勘探中场的正演计算就是给定地电条件和供电电源的位置,计算地下和地面的电位分布。例如对电阻率法就是给定地下的电阻率(或导电率)的分布和地形形状,计算电位的分布,进而对于具体的供电和测量电极装置计算视电阻率值,这种正演计算是反演计算的基础。但遗憾的是直流电法正演的解析计算往往很困难,对许多地电剖面情况甚至不可能求解  相似文献   

2.
利用传统有限单元法在有限空间范围内开展三维地球物理场正演模拟时,由于截断边界的影响,会引起局部异常的畸变,影响数值模拟的精度,对该问题通常采用扩边的办法加以解决,但需要的范围较大,从而大大增加运算成本,影响正演模拟效率。本文基于COMSOL Multiphysics软件,在求解域外部边界设置无限元以替代传统边界条件,达到减小计算区域目的。通过孤立球体和组合模型磁场正演模拟,考虑退磁、剩磁及地表起伏条件,与传统有限元方法相比,有限元—无限元耦合算法能够有效克服边界效应,提高计算精度,降低运算量,从而提高了有限单元法正演数值模拟效率。  相似文献   

3.
在弹性波传播问题的研究中,伪谱法是一种很好的模拟方法。该法的最大优点是节省存储空间,运算速度优于其它方法。但是伪谱法也存在一定的问题,由于折叠效应的影响,普通的吸收边界条件不能达到良好的效果。所以,这里采用一种"反周期延拓"与普通吸收边界相结合的新办法,用来消除折叠效应和边界反射的影响,并在TTI介质的伪谱法正演模拟中加以应用。模拟结果表明,这种新的方法能有效地达到消除伪谱法正演模拟过程中所产生的折叠效应,具有一定的推广价值。  相似文献   

4.
三维地电断面激发极化法有限元数值模拟   总被引:5,自引:0,他引:5  
用有限单元法求解三维地电断面激发极化法正演模拟算法.首先给出了三维构造中点源电场异常电位的边值问题与变分问题, 简化无穷远边界上的边界条件以提高计算速度及计算精度.以此为基础计算视电阻率对模型电阻率的偏导数矩阵, 并进行三维地电断面激发极化法正演模拟, 与等效电阻率法相比节约了计算时间.对几例较典型的地电模型进行计算, 结果表明本方法是正确可行的.   相似文献   

5.
在编写电磁法正演模拟代码这类高级且复杂的计算程序时,通常需要采用动态语言和静态语言相结合的双重编程语言来完成,这给研究工作者带来了诸多不便.为了提高编程效率,降低编程复杂性,这里采用全新的、能有效地简化编写工作量的Julia语言来编写大地电磁三维正演程序,正演代码实现由有限体积法(MFV)建立的基于各向同性介质的三维模...  相似文献   

6.
傅氏变换有限差分(即FFD)混合波场延拓法波动方程正演模拟精度高,但三维FFD因其计算复杂、数据存储巨大而难以推广。这里提出三维FFD一步法正演模拟并行算法,把一个串行的大规模频率循环分割成并行的多个小规模频率循环,并发执行多个波场延拓来提高运算速度,并通过由给定的资源设定并行数,以及适时重排三维数据,大大增强了算法的适应性和时效性。通过对比并行和非并行正演实验,验证了算法的有效性。  相似文献   

7.
F-K反偏移法正演模拟是在给定地质模型后通过构建人工偏移场,再对人工偏移场进行F-K反偏移来实现。人工偏移场的构建是在地质模型上直接放上经正确拉伸和幅度校正后的地震脉冲子波。F-K反偏移中插值映射采用两点sinc插值,其能避免假同相轴的出现且具有较快的速度。奇异点处理使用泰勒近似展开。最后对凹陷模型进行了F-K反偏移法正演模拟实现,得到了预期的效果。  相似文献   

8.
时间域有限元算法已被广泛应用于航空电磁三维正演模拟当中,然而由于航空电磁测区面积大,且采样密集,造成正演计算量巨大,传统的串行算法已经无法满足计算效率要求,为此,开展了并行加速算法研究以解决计算效率不足的问题。基于航空电磁系统的影响范围有限,采用局部网格技术将计算任务划分成多个子网格,即每个发射源一套网格,各网格的正演计算相互独立,不存在数据依赖性,具有很好的可并行性;利用MPI技术对多个子网格正演任务进行分配,在各个进程上进行并行计算;针对每个正演子网格,在进行时间域有限元算法正演模拟过程中,采用OpenMP技术对单元矩阵进行并行计算。典型地电模型的数值模拟结果表明本文开发的MPI+OpenMP并行正演算法可以有效提高正演速度,最高加速比可达10倍。  相似文献   

9.
点源二维电法正演的有限单元法   总被引:3,自引:3,他引:3  
点源二维任意地电断面电阻率法的正演计算是电法勘探中的一个重要问题。对于这个问题,我们已经实现了有限差分方法的计算,获得了较好的效果。有限单元法是解决电法正演计算又一重要途径,在1971年Coggon首先将有限单元法引进到电法正演的模拟计算中来,以后Rijo等人又发展了这一方法,有限单元法进行电阻率法正演计算的优点,是比较容易拟合实际复杂的地电断电和起伏不平的地形。在实现点源二维地电断面视电阻率异常正演计算的有限单元法中,我们将混合边界条件引入求解,并采朋了若干优化措施,使有限单元法计算视电阻率异常的精度和速度均有进一步提高,获得较好的计算效果。在文献中我们简单地介绍了这个方法,下面我们较详细  相似文献   

10.
矿井电法三维有限元正演模拟--直流电透视方法技术研究   总被引:1,自引:2,他引:1  
通过工作面两巷道间直流电透视三维正演模拟、资料解释方法以及正演实例分析,表明不同装置形式的直流电透视技术仍属于体积勘探。但不同的测量方式对应于不同的电场分布特性,且能分别突出探测目的层(煤层内或顶底板岩层内)的异常特征。为解决采煤工作面附近的具体地质问题,使用相应的方法可提供理论依据。  相似文献   

11.
2.5维起伏地表条件下坐标变换法直流电场数值模拟   总被引:3,自引:1,他引:2  
研究起伏地表对视电阻率分布的影响是进行地形校正的基础。由于很难处理不规则边界,计算简洁且效率高的有限差分方法很少用于解决起伏地表问题。为了解决该问题,引入曲化平思想,实现了一种基于坐标变换法的起伏地表条件下的直流电场数值模拟方法。方法从传统的2.5维基本方程及边界条件出发,通过坐标变换将起伏地表问题转化为水平地表问题,并利用有限差分法求解水平地表问题,最后再通过坐标映射得到起伏地表条件下的电位和视电阻率分布。精度分析及计算实例表明:本方法实现简洁,计算误差主要集中在震源附近,整个计算区域内的平均相对误差为1.39%,计算结果满足起伏地表条件下的电位和视电阻率的分布规律。  相似文献   

12.
影响井下直流电法超前探测效果的因素较多,巷道空腔对全空间电流场分布的影响就是其中之一。论述了有限单元法模拟井下稳定电流场正演模型,介绍了全空间条件下三维模型中的点源电场的边值及变分问题;模拟在巷道内加载定点电流源,通过巷道影响因子研究巷道空腔对全空间稳定电流场分布的影响。数值模拟结果表明,巷道影响与巷道几何尺寸大小、供电电极布置及巷道围岩导电性有关。   相似文献   

13.
三维地质建模具有直观、形象、便捷等优点,在地质、矿产、地下施工等领域有着突出贡献。COMSOL Multiphysics作为多物理场耦合模拟仿真软件,具备一定的三维建模能力。通过充分挖掘该软件自身潜力,结合Kriging插值方法,尝试建立了山西某井田的三维地质模型。研究结果显示,基于COMSOL Multiphysics的三维地质建模方法可实现地质体模型的缩放、旋转和平移等操作,生成地质剖切面图、透视图,展示地质体内部的各个细节。实例证明了COMSOL Multiphysics软件的三维地质建模能力及其在地质解译、矿产评估和多物理场求解等方面的应用价值。   相似文献   

14.
直流电阻率法因效率高、成本低等特点广泛应用于有色金属、煤田等矿产资源勘探和地质调查行业。但在井–地、地–井电阻率法勘探过程中,钻孔因素(井液电阻率和钻孔孔径)对不同探测模式下视电阻率响应的影响尚未明确,是否影响视电阻率资料的解释是值得探讨的一个问题。因此,提出采用混合网格有限元法实现直流电阻率三维正演。给出异常电位法满足的边值问题及有限元变分问题,应用三棱柱和四面体混合网格实现对计算区域的快速离散,并建立2种网格的线性插值基函数和单元系数矩阵;采用SSOR–PCG迭代算法求解异常电位满足的大型线性方程组,得到各观测点的响应。在保证计算精度的前提下,应用混合网格有效地离散钻孔地电模型,探讨钻孔因素对井–地、地–井观测方式视电阻率数据的影响特征。对于井–地观测方式:钻孔附近视电阻率受钻孔因素影响最大,严重影响视电阻率资料的合理解释,但随着测点收发距的增加,视电阻率响应逐渐趋于围岩电阻率;而且随着发射源深度的增加,钻孔对地表视电阻率响应的影响也逐渐减小。对于地–井观测方式:钻孔因素对浅部数据影响较大,对深部数据影响小,相比于井液电阻率,视电阻率响应更易受钻孔孔径的影响。直流电阻率法混合网格有限元三维正演算法的提出,对实际井–地、地–井电阻率法勘探具有理论指导意义,可结合钻孔信息进行正演模拟分析,选择合适的收发距,有效压制钻孔因素对实测资料的影响。   相似文献   

15.
二维起伏地表条件下坐标变换法直流电场数值模拟   总被引:6,自引:2,他引:4  
针对有限差分方法处理不规则边界难的问题,引入曲化平的思想,实现了一种新的基于坐标变换法的起伏地表条件下直流电场数值模拟方法.通过坐标变换,将物理空间的直流电场满足的基本方程和边界条件变换到计算空间;在计算空间内,利用有限差分方法对变换后的方程和边界条件进行数值分析;最后,将计算结果由计算空间映射回物理空间.精度分析及计算实例表明:本文提出的算法精度高,计算误差主要集中在源附近,且当源离地表较近时平均相对误差在1.0%以下;计算结果满足起伏地表条件下电位及视电阻率的分布规律.  相似文献   

16.
电阻率图像重建的佐迪(Zohdy)反演及其应用效果   总被引:19,自引:0,他引:19  
王兴泰  李晓芹 《物探与化探》1996,20(3):228-233,208
本文讨论了利用佐迪反演方法对地表高密度电阻率法的数据采集结果进行电阻率图像重建。通过对多种复杂地电模型的数值模拟和野外实例的反演分析,表明该方法简单易行,计算成本较低,成像的分辨率高,效果较理想,将其用于地表测量的电阻率层析成像完全可行。作者用Fortran语言编制了实现上述过程的程序,对所有数学模型的正演计算均采用2.5维有限元法。为突出成像效果,本文在结果图示时均采用了伪彩色图像处理技术。  相似文献   

17.
电偶源是研究电磁法探测原理中常被使用的一种主动源。在实际工作中,常采用供电导线向大地发射电流,当发射端与接收点距离远大于供电导线两极长度时,才可将供电导线看作是电偶极子。但在很多情况下并不满足此条件,使基于电偶极子的场计算公式不再适用,因此正反演中需要先将长导线源切分成若干个电偶极子,然后对各偶极子的场响应进行叠加。对此,这里研究并实现了将长导线源切分成多个电偶极子的算法,以二分法为原理的递归切分算法切分出的电偶极子数目相对较多;以穷举法为原理的切分算法在穷举步长较小(如0.1m)时能较准确地切分出电偶极子的位置。两种算法都可应用于长导线源正反演场的计算。  相似文献   

18.
山西某矿位于构造堆积盆地,属黄土匠陵地貌,区内被第四系黄土覆盖,为解决该矿地质构造及水文地质问题,采用三维地震勘探与瞬变电磁法、直流电测深法相结合的勘探方案。利用三维地震勘探解决目的层赋存形态及其地质构造问题,再利用瞬变电磁法进行平面控制,以使平面及深度解释误差达到勘探要求;应用直流电测深法及瞬变电磁进行勘探,在正反演解释的基础上对其资料综合分析、对比,并结合三维地震资料,确定地质构造的富水及导水性。  相似文献   

19.
三维坑道直流聚焦法超前探测技术   总被引:1,自引:1,他引:0  
对坑道超前探测方法中的三维坑道直流聚焦法的电极组合作了研究.计算、分析了2个同性点电源的电流分布规律,指出由于电流的相互排斥作用,电流密度在这两个点电源的中垂线下一定深度处产生最大值,这个深度与这两个点源之间的距离成正比.基于上述理由,设计出3种电极组合:四点电源、五点电源和九点电源布极方式.用三维有限单元数值模拟的方法计算了上述3种电极系的空间电位分布,发现九点电源方式对电流聚束效果较好,而且当主电极电流小于屏蔽电极电流时效果最好.  相似文献   

20.
高精度、高效率的深海勘探是当前国家海洋资源勘探开发的研究热点之一,设计观测效率高、勘探成本低的探测方式,开发可精确计算海底复杂电性环境的正演模拟算法有助于推进该研究的进展。本文基于前人的研究经验,引入了全拖曳式深海直流电阻率观测方式,同时考虑到海底沉积环境的电各向异性,开发了基于非结构有限元方法的三维任意各向异性深海直流电阻率正演算法,实现了对海底任意电各向异性情况的仿真模拟。本文算法与层状一维模型半解析解的对比验证了算法的精度。对典型各向异性电性模型的仿真模拟和分析,证明了全拖曳式深海探测方式对海底沉积层和矿产资源电各向异性的高分辨率。在此基础上,总结了电导率各向异性主轴沿笛卡尔坐标旋转时视电阻率分布也发生同方向旋转的视电阻率分布特征,并根据海底地形对各向异性探测的影响进行了简单探讨。  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号