共查询到20条相似文献，搜索用时 15 毫秒
1.
解重磁反演问题的人机结合选择法 总被引：1，自引：0，他引：1
林振民 《物探化探计算技术》1983,(1)
(一)选择法的简要评述选择法是解重磁反演问题的一种老方法。它的做法是根据某一剖面上所观测的异常的特点,假定地下地质体的形状和物性参数;对这个假设的地质体做正演计算,并将计算的异常与观测的异常进行比较;如果它们之间相差很大,就修改原来所假定的地质体模型, 相似文献
2.
《物探化探计算技术》2017,(2)
解决任意密度分布复杂地质体重力异常三维正演快速、高精度计算问题,是实现重力三维反演、人机交互解释建模的关键。针对该问题,从积分方程出发,提出一种波数域重力异常三维正演方法,其关键环节包括三个方面:(1)将研究区域剖分成许多规则小棱柱体,每个小棱柱体密度值可以任意给定,以此刻画任意密度分布和起伏地形条件下的复杂地质体;(2)给出一种新的高精度均匀棱柱体重力异常二维波数域的计算公式,用于计算组合棱柱体模型的重力异常;(3)采用Gauss-FFT法将重力异常从波数域转换到空间域,保证计算效率的同时,有效克服了传统FFT法引起的边界效应问题。模型算例检验结果表明,该算法计算速度快、精度高,对于剖分为百万个棱柱体的模型,耗时只需几秒。 相似文献
3.
位场相关成像是根据实测异常与地下不同位置地质体所产生异常之间的相关系数来快速获得地质体的空间位置.现有的相关成像方法是利用球体模型来模拟地下地质体的形状,当场源体的实际形状与球体相差较大时,计算结果势必出现较大误差.为了解决这一问题,对该方法进行改进,以不同模型来模拟地下地质体形状,计算其产生异常与实测数据的相关系数,理论上使相关系数取得最大值的模型与实际地质体情况一致.因此,改进后的方法不仅可以获得地质体的位置参数,还可以对地质体的类型(构造指数)进行估计.磁异常的相关成像计算采用异常的解析信号来完成,这样可有效地避免磁化方向的干扰,且计算公式相对简单.通过理论模型试验,证明此方法可以成功地完成位场数据的反演工作,且稳定性较高.最后将其应用于上海实测磁异常数据的解释中,获得了地下未爆炸物的分布情况. 相似文献
4.
为探究三维地质体对二维反演结果地影响以及二维反演中极化模式的选择问题,设计了一系列二维、三维地电模型进行正演计算,并对正演结果进行二维反演计算。研究结果表明:相位受到三维畸变的影响较小,视电阻率较大;当三维异常体模型走向延伸较小时,TE模式的反演结果误差较大,TM模式的反演结果相对较好,但是反演得到了虚假的下覆构造,TE+TM联合反演结合了TE和TM模式反演的优点,可以较好地反映异常体的赋存位置及电阻率值;当三维异常体模型走向延伸长度增大时,三维数据的二维反演结果是可靠的,用TM和TE+TM模式反演都比较合理。 相似文献
5.
重磁场二次导数法是对复杂重磁异常进行滤波的一种重要手段,研究基底构造,突出隐伏在区域背景场中的局部异常,都能取得良好效果.在已发表的文献中,对各类地质体的二次导数异常的特征,均未作过详细介绍.本文从理论上并结合一些实践经验,探讨几种形状规则的地质体的二次导数异常的特征,并对如何利用重磁场二次导数异常来推断地质体的度数、埋深、边界、构造线部位等问题提出一些看法.另外还讨论了在计算时选择最佳计算半径等问题. 相似文献
6.
根据理论模型对剩余异常进行直接推断解释的问题,在于选择模型的大小和密度差,以使它的重力效应尽可能地与剩余异常图形吻合.这基本上是一种实行曲线或曲面拟合的方法,然而,它绝不是一种简单的方法.由于即使用公式来描述最简单模型的位场都过于繁杂,因此,给任何一种直接拟合法都带来了困难.为了说明曲线直接拟合方法的困难,让我们来研究一个经常用来表示岩脉或矿脉的模型.一般认为,岩脉或矿脉是一种宽度远小于其长度的地质体,通常用薄板模型来表示这种地质体,其图形如图1所示.因为重力效应有可能比其他类型的位场观测结果对这种地质体的走向长度更为敏感,所 相似文献
7.
重力资料定位地质体边界问题的探讨 总被引：4，自引：0，他引：4
针对重力资料应用于地质体的边界定位问题进行了系统的研究,通过大量二维、三维模型的正演计算,统计分析了重力各异常特征值点边界定位的误差,最后给出了每种方法精确定位边界的应用前提。 相似文献
8.
对于地面实测的重磁异常剖面,在准确知晓其对应地质剖面中地质体的密度及磁化强度等物性的前提下,结合区域地质条件、区内地质背景及剖面上(周边)钻探(坑道开采)资料,建立初始地质剖面模型,利用重磁异常精细正反演,可较准确地确定剖面上的地层、岩矿等地质体界线。对于已开采的(磁)铁矿,可以通过正演,计算出磁异常并与实测异常进行比较,确定是否有剩余异常存在,从而判定是否有(磁)铁矿体存在。以梅山铁矿为例,对其重磁异常资料重新进行数据处理及精细正反演计算,结合地质及物性资料进行地质解释,利用其成果布设钻孔,取得了明显的地质找矿效果。 相似文献
9.
激电扫面实际测量时,由于测线布设方位和地下地质异常体的走向不尽一致,从而导致测量结果不同,给物探异常的解释带来了困难,甚至会引起误导,因此这里利用数值模拟方法,计算了激电中间梯度法测量条件下二度异常体在不同方位测量条件下的三维响应规律。首先实现了有限单元法三维正演的公式推导及其数值解求解;其次通过一维模型计算了数值解和解析解的误差,验证了本文中理论计算方法的有效性;最后设计了三维模型进行了理论计算,发现了激电中间梯度装置条件下二度异常体的走向和测量剖面方向成不同夹角时的面上响应规律。研究结果表明:随着两者夹角的增加,激电异常体的规模和形态逐渐与设计异常体一致,且两者夹角在60°~90°之间探测的实测数据较为可信。 相似文献
10.
利用层析成像的基本原理对地质异常体进行速度层析成像正、反演研究,计算结果可以确定出地质模型中速度异常体的大小、位置和慢度值等物性参数。首先通过建立正演模型,将含有异常体的模型区域剖分成规则的矩形网格,在区域两边的钻孔中分别安置发射和接受装置,采用多点发射、多点接收的方法获得已知地质条件下的每条速度射线的走时值,而后根据线性回归的原则对所得的时间数据加入20%的噪声,采用医学中成熟的ART算法进行模型的反演计算,计算表明反演结果与正演模型十分地逼近,计算精度和速度都能达到满意的效果。 相似文献
11.
中间梯度法是野外常用的一种电阻率勘探方法,当其测量范围在供电电极的中间地段(即长度为AB/3)时,电场即可近似地看成是均匀的。如果地质断面是二维的,则可采用本文的方法做视电阻率异常的正演计算。本方法与国外已有的方法相比,具有模拟功能强和计算量小等优点,可以对二维复杂地形条件下存在多个地质体的情况进行模拟,特别是略去了对偏导数项((?)u/(?)n)的计算,从而大幅度节省了计算量。 相似文献
12.
为探究激发极化中赛格尔公式的适用性及饱和效应的问题,笔者分析了赛格尔公式推导过程中所做的近似,通过正演模拟分析赛格尔公式的适用范围;并研究了饱和效应问题。结果表明:赛格尔公式不适合用于高极化地质体或复杂极化地电结构模型计算,特别是电性变化剧烈的地电结构情况;饱和效应在理论上存在,且其饱和值与极化体大小有关系。极化体越大,饱和时所对应的地质体电阻率相对围岩电阻率就越低,饱和效应的存在使得激电异常的极大值可能会偏离富集矿体中心。 相似文献
13.
《物探化探计算技术》2016,(2)
由于重力梯度张量测量相对于重力异常测量有许多优点,因此对重力梯度的研究十分必要。在重力梯度各张量的正演计算中,复杂地质体的计算式难以直接推导,而利用有限元技术在复杂形体体积积分的优势,可以较为快速和简便地进行复杂模型的重力梯度全张量正演计算。通过模型的建立和正演计算,分析了球体模型梯度张量异常的平面和剖面特征,以及重力梯度张量与球体模型位置的规律;最后进一步利用复杂模型的重力梯度正演模拟,说明了重力梯度识别地下地质体位置的优势和不足。 相似文献
14.
正演数值模拟是反演成像的基础。为了实现磁法勘探精细化反演成像与定量解释,本文利用把一个大问题分解为多个小问题的思路,提出一种基于泊松方程的高效、高精度空间波数混合域二度体磁异常数值模拟方法。该方法利用傅里叶变换把二度体磁位偏微分方程转换为一维常微分方程,并采用基于二次插值的一维有限单元法求解该方程,进而通过反傅里叶变换得到空间域磁异常。在模型算例中,分别设计截面为矩形的常磁化率和变磁化率二度体模型,针对本文算法的计算精度和计算效率进行了验证。模型算例结果表明：该算法计算精度高,相对误差绝对值均小于1%;计算速度快,网格节点剖分2 501×2 501的模型模拟时间为4.18 s;适用于任意复杂地形模型。 相似文献
15.
众所周知,米可夫量板是一种适用于走向长度无限的二度磁性体异常正演计算的简便工具.若将米可夫量板用于走向长度有限的三度磁性体,则由于二度模型与实际的三度体相差甚远,使计算结果产生较大的误差,甚至导致错误的解释.似二度量板克服了米可夫量板的局限性,可以计算任意截面、任意有限走向长度的似二度体异常,在实际工作中提高了异常解释精度.但是这种计算方法要用较多的量板,而且量板的构制数据与磁性体的走向长度和 相似文献
16.
二度体磁异常快速反演方法 总被引：5，自引：0，他引：5
本文提出一种二度磁异常快速反演方法。该方法不是利用磁异常直接反演,而是通过构造一个新函数,该函数与异常源地质体各欲求参数间存在比较简单的函数关系,从而将复杂的非线性问题化为简单的线性问题,提高了反演速度。基于这一原理,本文就二度岩脉、接触带、垂直断层磁异常的反问题导出其表达式。具体解法归结为:先求解一个线性方程组,再通过代数运算,求出欲求异常体的参数。通过理论模型与实际资料的反演,证实了该方法是正确而有效的,且适用于横向具有一定间隔的组合模型的反演。 相似文献
17.
《吉林地质》2016,(2)
无网格方法是伴随着工程计算领域的深入发展而应运而生的一种新兴数值计算方法,作为有限元等传统方法的补充和发展,在近十几年里得到了广泛的关注。本文将无单元Galerkin法(EFGM)应用到二维大地电磁正演模拟中,推导了对应无单元法的等价线性方程组,采用滑动最小二乘法构造形函数,本质边界条件采用罚因子法处理,编制了对应EFGM的程序,通过一维层状模型验证了算法的正确性,通过倾斜接触面大地电磁模型的计算和分析,验证了无网格法计算复杂模型便利的特点,对于地下倾斜低阻异常体,TM模式能较好的反映异常体的产状,而TE模式则不能反映。对于地下倾斜高阻异常体,TM、TE模式均不能反映出其产状。 相似文献
18.
地质体“基底构造”的计算是区域重力解释中的重要任务之一。本文以深部二维重力异常的正演计算为例,用矩量法(MOMENT METHODS)求解积分方程反演其“基底构造”。计算结果表明,用矩量法中的分段基和点选配法反演深部“基底构造”,特别是提供“基底”的平均深度是一种卓有成效的直接方法。 相似文献
19.
20.
对成矿过程进行数值模拟实验是当前开展成矿作用研究的新兴手段之一,其前处理阶段即计算用地质网格模型的建立往往成为制约后续实验过程的一个瓶颈,原因在于：一、地质体的形态不规则,甚至是复杂的,对其进行三维建模难度较大,所成模型也趋于简化处理,从而影响计算结果的精确性;二、专业的数值计算软件其建模模块稍弱,不足以应对复杂地质体模拟,而专业的建模软件往往又不具备模拟成矿过程的数值计算功能。因此需取长补短,对不同软件加以综合运用,实现对地质体的计算用网格模型构建。 相似文献