物探、化探软件交流

有限单元法、有限差分法及积分方程法程序这些大型实用电法勘探正演程序适用于计算复杂地形条件下,任意剖面、各种装置的视电阻率和视极化率值,这些程序分别在DUAL—68000、IBM—PC/XT、UNIVERSE—68等微机上调试成功,备有程序软盘供用户使用。深部地震数字回放系统软件深部地震数字回放系统软件在DUAL—68000微机上调试成功,并用于生产。该系统能将野外的调频模拟磁带记录回放成数字磁盘记录,经带速数字校正后,自动绘出地震剖面图。本系统备程序软盘。     

电法程序交流

本刊编辑部备有点源二维有限单元法、有限差分法和点源三维积分方程法的Fortran语言程序及BD—200语言程序各一套。这些程序适用于复杂地电断面的正演计算,也可用于反演计算。有限单元法程序还适用于计算复杂地形条件下、任意剖面(包括主剖面和旁侧剖面)的各种装置的视电阻率和视极化率值。该程序具有功能强,操作简便等特点。其计算精度和速度都超过了国外同类型程序的水平,并于1982年由地矿部物化探局评审通过。凡需要该程序以及咨询的用户,请来函向成都地质学院内《物化探计算技术》编辑部联系。     

海洋大地电磁二维正演及结果分析

研究了海洋大地电磁法的二维正演并编制了相应程序 ,程序采用有限单元法求出海底各节点处的电磁场值 ,进一步计算得到视电阻率。为了验证程序的正确性与有效性 ,根据卡尼亚理论推导出计算一维海洋模型的阻抗公式 ,并得到视电阻率的解析解 ,将其与有限元法得出的数值解进行了对比 ,表明了此方法的正确性。通过对正演结果的分析 ,表明海洋大地电磁测深能够成功地反映海底的电性结构     

点源二维电法正演的有限单元法

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

点源场电阻率法二维地形改正的边界元法

边界元法是一种新的数值计算方法。本文讨论了该法的基本原理。由于边界元法可降低数值模拟的空间维数,因而在地形改正中,边界元法与区域型解法(例如。有限单元法和有限差分法)比较,该法具有基本方程组阶数低、输入数据少、精度高、速度快等优点。本文用解析法验证了角域地形上视电阻率数值解的正确性。实践表明:用边界元法获得的地形改正的视电阻率曲线突出了矿异常;点源场电阻率法2—D地形改正的边界元解法是一种经济有效的方法。     

地下水向水平井三维流动的有限层分析

基于地下水向水平井三维流动数学模型,利用半解析数值分析原理推导了层状各向异性承压含水层有限层方程,给出了水平井源汇项计算公式,编制了相应的计算程序,实现了有限层方程的求解。通过算例的解析解和有限差分解验证了有限层法及程序的正确性,对比分析了有限层法与有限差法的计算效率。结果表明,有限层法适于层状各向异性承压含水层中地下水向水平井三维流动问题的模拟分析,且具有较高的计算效率。     

PC—1500微型计算机电测深曲线数字解释程序

利用电阻率转换函数便于在电子计算机上实现视电阻率测深曲线的数字解释。许多野外物探工作者常因野外条件所限,极需有一种采用简单计算工具即可进行视电阻率测深曲线数字解释的方法。本文介绍改进的电阻率转换函数消层法及其PC——1500微型计算机程序。本方法在理论模型和实际地电断面上的使用表明,其计算程序简单可靠,适于野外队使用。一、数学原理简述如图1所示,水平层状介质表面任一点0上的视电阻率函数,在施伦贝尔热装置情况下,其积分表达式为:     

非线性问题的有限单元法在工程地质中的应用

使用有限单元法对岩石基础、边坡和巷道等问题进行稳定分析在国内外做了大量的工作。本文结合我们自己的工作,以平面问题的位移法为例,介绍物理上非线性物体的有限单元分析的基本原理,并对某些具有非线性应力应变关系的岩体介质以及诸如节理、断层等间断面提出一些简单的力学模型。此外,还介绍几个和这些力学模型有关的通用程序以及使用它们得到的某些工程的计算成果。     

一种基于有限元应力场的边坡滑动面搜索方法

如何将有限单元法计算得到的结果是应力场、应变场和位移场等计算成果应用于边坡稳定安全系数及滑裂面位置的确定,是有限单元法边坡稳定分析中的一个重要研究课题.提出一种基于有限元应力场的最危险滑裂面的搜索方法,编制了搜索程序,采用澳大利亚边坡稳定分析程序的考核题进行对比验算,验证了该方法的正确性及程序的可行性.     

适用于袖珍计算机频率测深数字滤波方法

详细推导了电偶源频率测深单分量电场定义的视电阻率数字滤波计算公式,并且应用最小二乘法设计了只有30个系数的滤波器,编制出适用于袖珍计算机的视电阻率计算程序,使计算速度大为提高.与用均匀半空间的精确公式计算的视电阻率对比,最大相对误差<2%.     

双排桩支护结构的计算方法研究及工程应用

总结了目前基于弹性土抗力法的双排桩支护结构常用计算模型,分别编制了平面杆系有限差分法计算程序,介绍了双排桩支护结构在某实际工程中的应用,比较了这些模型之间计算结果的差异。通过与实测数据的对比,提出了较为合理的模型,并对双排桩支护结构的诸多影响因素进行了计算分析,所得结论对双排桩支护结构的优化设计具有参考价值。     

电阻率随深度变化的电偶源频测视电阻率曲线

本文给出了当第二层电阻率随深度变化时，二层及三层大地垂直磁场Ｈｚ表达的视电阻率公式。对第二层电阻率随深度呈线性变化和指数变化的情况，利用线性数字滤波法编写了ＦＯＲＴＡＮ－７７计算程序。计算了一些视电阻率曲线，对部分视电阻率曲线进行了讨论，得到了一些有意义的结果。     

磁性源瞬变电磁法视电阻率计算方法

瞬变电磁法视电阻率的定义和计算主要有早期、晚期和全区视电阻率定义和计算方法,然而根据瞬变电磁法激发场源、工作装置和工作环境的不同,其视电阻率的定义和计算也不相同.笔者总结和评述了各种磁性源不同工作装置的早期、晚期以及全区视电阻率计算方法,同时对隧道、矿井瞬变电磁法中的全空间视电阻率的计算以及考虑关断时间效应的视电阻率的计算等进行总结.在分析各种定义和计算方法的基础上,结合瞬变电磁法的发展应用,对视电阻率的定义和计算提出了几点建议.     

MT数据处理中静校正方法对比

在大地电磁测深数据处理中,静态效应会影响解释结果.这里主要采用有限单元法进行正演计算,采用阻抗相位校正法、空间滤波法以及视电阻率导数解释法,对存在静态效应的理论模型进行数据处理.对比结果表明,前二种方法对静态效应起到了一定的压制和改善作用,且阻抗相位校正法效果相对较好,而通过视电阻率导数解释法,可以准确地反映地下真实异常形态.     

CSAMT一维全频视电阻率反演

本文采用全频率视电阻率反演,使用2~(-3)~2~(13)Hz频率,并且不做近区和过渡区校正,直接用卡尼亚视电阻率作为反演参数。正演采用虚界面法计算有限长导线在水平层状介质的电磁场分量,反演采用有限内存拟牛顿法。数值模拟采用三层模型和实测数据进行反演,得到了符合地下电性变化的反演模型结果,证实了CSAMT一维全频视电阻率反演具有可实施性。     

有限差分法(第一讲)

电法数据处理,目前正处在研究阶段。这些方法包括有限差分法、有限单元法和积分方程法。由于电法计算工作的复杂性,计算工作量大,通常需要用大中型计算机进行处理。这些方法除了解决电法问题外,还能用于其他物探方法及地质领域内,例如地震、重磁、测井、水文工程地质、岩石力学等问题,而且可以预料,它们将在地质领域内得到更加广泛的应用。由于这些计算方法具有较普遍的意义,为此,本刊特邀请有关同志介绍电法的正演计算方法,以达到触类旁通的目的。从本期开始,以连载讲座形式向读者介绍这些方法。     

井地电位2.5维有限元数值模拟异常分析CSCD

从井地电位法满足的边值问题出发,利用有限元推导了井地电位2.5维的详细解法,并编制了相应的程序。然后利用该程序对一水平层状的地电模型进行了正演计算,并与解析解做了对比分析,验证了程序的正确性和有效性。在此基础上对均匀介质中分别存在高阻板、低阻板以及高阻板与低阻板同时存在的模型,进行了正演计算,得到了视电阻率拟断面图及点电流源在不同深度的视电阻率曲线分布图。结果表明,该方法能够对地下介质进行有效的探测,能够精确地确定异常体上边界、下边界的埋深位置。     

Zohdy方法在三维电阻率测深数据反演中的应用

三维电阻率测深数据采用Zohdy方法近似反演,正演计算采用有限单元法,反演初始模型由测量视电阻率数据给出,通过比较实测视电阻率值和预测模型计算的视电阻率值对数差来修改模型网格电阻率.为了解决任意电极距测深数据的反演,采用大、小双网格剖分,大网格反映地下电性分布情况,小网格用于实际有限元正演计算.文章对加5%随机噪声的模型理论视电阻率测深断面数据和实测视电阻率测深断面数据进行反演.结果表明该反演方法需要的内存小、计算量少,并能够取得很好的反演效果.     

电测深法视电阻率对介质电阻率偏导数的计算

在前人工作基础上,对电测深法一维正演解析解递推电阻率转换函数求偏导,给出了电测深法视电阻率对介质电阻率偏导数的解析计算方法;根据Seigel体激发极化理论,采用视电阻率对介质电阻率偏导数矩阵的线性组合,计算出了电测深法的视极化率,其结果与"等效电阻率法"计算视极化率最大相对误差小于1.1%。在一次正演计算过程中就可以同时得到视电阻率、视极化率以及视电阻率对介质电阻率的偏导数矩阵,提高了计算效率,也为后续的多层地电模型电测深法数据优化反演计算提供了关键性技术。     

井地电位2.5维有限元数值模拟异常分析

从井地电位法满足的边值问题出发,利用有限元推导了井地电位2.5维的详细解法,并编制了相应的程序.然后利用该程序对一水平层状的地电模型进行了正演计算,并与解析解做了对比分析,验证了程序的正确性和有效性.在此基础上对均匀介质中分别存在高阻板、低阻板以及高阻板与低阻板同时存在的模型,进行了正演计算,得到了视电阻率拟断面图及点电流源在不同深度的视电阻率曲线分布图.结果表明,该方法能够对地下介质进行有效的探测,能够精确地确定异常体上边界、下边界的埋深位置.