首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 93 毫秒
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种网格的线性插值基函数和单元系数矩阵;采用SSOR–PCG迭代算法求解异常电位满足的大型线性方程组,得到各观测点的响应。在保证计算精度的前提下,应用混合网格有效地离散钻孔地电模型,探讨钻孔因素对井–地、地–井观测方式视电阻率数据的影响特征。对于井–地观测方式:钻孔附近视电阻率受钻孔因素影响最大,严重影响视电阻率资料的合理解释,但随着测点收发距的增加,视电阻率响应逐渐趋于围岩电阻率;而且随着发射源深度的增加,钻孔对地表视电阻率响应的影响也逐渐减小。对于地–井观测方式:钻孔因素对浅部数据影响较大,对深部数据影响小,相比于井液电阻率,视电阻率响应更易受钻孔孔径的影响。直流电阻率法混合网格有限元三维正演算法的提出,对实际井–地、地–井电阻率法勘探具有理论指导意义,可结合钻孔信息进行正演模拟分析,选择合适的收发距,有效压制钻孔因素对实测资料的影响。   相似文献   

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

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

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

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

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

17.
电偶源频率电磁测深三维地电模型有限元正演   总被引:9,自引:1,他引:8  
在有源变频测深(包括可控源音频大地电磁测深CSAMT)的三维有限元正演中,应用了吸收边界条件和边元有限元算法,计算精度为3%左右,基本上可满足遇极源三维地电模型情况下频率电磁测深响应特征的分析。同时由于广义变分概念的引入,可使泛函的变分原理应用于有耗煤质电磁场问题。  相似文献   

18.
频率域电磁剖面有限差分法2.5维正演数值模拟   总被引:1,自引:0,他引:1  
对于实际生产中遇到的2.5维问题,从麦克斯韦方程组出发,把地电参数变化小的走向方向转化到波数域,用一系列波数模拟三维源的特征,并在波数域中,得到2组关于Hx(kx,y,z)和Ex(kx,y,z)的偏微分方程。选取适当的kx值,用有限差分法在y-z平面的网格中求解,再通过反傅里叶变换得到空间域中的电磁场。在验证了算法的正确性之后,对不同埋深的直立异常体、倾斜异常体及断陷模型进行了数值模拟,其结果直观地显示了异常体磁异常的形态,同时研究了程序对低阻异常体的横向分辨率。对层状大地模型,用井间收发方式进行了模拟研究,并取得了较好的结果,这对今后实际勘探应用奠定了基础。  相似文献   

19.
在胶莱盆地,其红层所属地层岩性为粉砂岩、泥岩,含、透水性较差,查找低阻断裂构造是其地下水勘探的唯一途径。而红层地区的压扭性断层,会因挤压密实的断层泥存在,导致单纯利用低电阻率参数判定地下水的赋存位置出现差错。根据该区地下水的自然电位异常值约为-100mV、含水破碎砂岩极化率大于2%、衰减度D大于50%的特点,提出了自然电位观测法与激发极化电测深法相结合的组合方法。即利用自然电场法进行剖面测量,判定含水异常位置,然后利用激发极化电测深法准确判定地下水所在层位及埋深。在胶州市某镇政府水源勘探中,利用该组合方法,成功的圈定了3个井位,并在预定的深度内打出了单井出水量为1000t/d的优质井。  相似文献   

20.
为了研究起伏地表条件下三维直流电场问题,基于全空间的拟解析近似理论,推导出起伏地表条件下三维直流电场的拟解析近似理论公式。将起伏地表条件下地下存在异常体时所产生的异常电场分为两部分,其中一部分为地下三维异常体所产生的体异常电场,通过对异常体的体积分求得;另一部分为由起伏地表所产生的面异常电场,可以利用对起伏地表地空界面的面积分求得。利用拟解析理论,对两部分异常分别推导体电反射函数和面电反射函数,得到起伏地表下电场积分方程的拟解析近似公式。该理论公式为解决起伏地表问题提供了一种思路,丰富了拟解析理论的应用,为直接进行正反演模拟、消除起伏地表对电法勘探的影响打下了基础。  相似文献   

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

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