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

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

用改进的佐迪反演方法进行二维电阻率图像重建

针对佐迪反演方法在二维拟断面应用中的缺陷进行了三个方面的改进：固定佐迪反演方法中用于转换电极距和深度关系的深度转换因子,可使绝大多数数字模型的反演过程稳定收敛;在调整视电阻率的公式中加入一个迭代系数,可使计算值快速趋向于观测值,减少了迭代次数;对于含有噪声的数据,加入滤波因子可以在一定程度上抑制噪声。通过数字模拟和野外实测数据对改进后的佐迪反演方法在二维电阻率反演应用中的可行性进行了验证。模拟研究包含了简单模型,也包含一些复合模型,二者均取得较理想的应用效果。     

三维地质体对AMT二维反演的影响研究

为探究三维地质体对二维反演结果地影响以及二维反演中极化模式的选择问题,设计了一系列二维、三维地电模型进行正演计算,并对正演结果进行二维反演计算。研究结果表明:相位受到三维畸变的影响较小,视电阻率较大;当三维异常体模型走向延伸较小时,TE模式的反演结果误差较大,TM模式的反演结果相对较好,但是反演得到了虚假的下覆构造,TE+TM联合反演结合了TE和TM模式反演的优点,可以较好地反映异常体的赋存位置及电阻率值;当三维异常体模型走向延伸长度增大时,三维数据的二维反演结果是可靠的,用TM和TE+TM模式反演都比较合理。     

中条山铜矿某工区激发极化法三极测深及其三维反演效果

文章简单描述了时域激电法原理、技术和测区地质概况及岩(矿)石电性特征.单极-偶极装置测深电阻率和充电率二维反演模型能够有效地圈定电性体在纵向和横向的分布特点,依据二维测深数据进行三维反演可获得电阻率和充电率异常曲线图,推测了不同深度范围内断裂和主要岩性分界,极化体空间分布位置和走向,钻探验证见铜矿化于划分的激电异常带,激电二维测深及其数据三维反演取得了较好的地质效果.     

直流电阻率测深二维自适应正则化反演

直流电阻率测深二维反演中,正则化参数的选取影响反演结果分辨率及反演过程稳定性。利用主动约束平衡正则化因子,进行直流电阻率光滑约束最小二乘二维自适应反演,改善直流电阻率测深二维反演的分辨率与稳定性。在反演迭代过程中,正则化因子根据模型参数的空间展布函数进行自适应计算、正则化参数的自适应计算。模拟数据反演结果验证了该方法的有效性与可行性,反演结果能准确地反映地下模型的真实电性结构。     

基于积分方程法的大地电磁三维反演

采用积分方程法计算均匀半空间中三维异常体的大地电磁响应,然后利用阻尼最小二乘法反演得到地下三维异常体电阻率分布。在计算灵敏度矩阵时,将灵敏度矩阵分为线性项和体现反演过程非线性的偏微分项进行求解,提高了计算精度。理论模型反演结果表明,反演较准确,对初值依赖性小,且只针对异常体剖分单元进行反演,反演效率较高。     

航空瞬变电磁数据一维Occam反演

分析了奥克姆(Occam)反演法的原理和计算方法,并将其应用于时间域航空电磁一维反演中,给出了一种简便的正则化参数搜索方法,对迭代过程进行了改进,大大提高了反演速度,5次迭代约20s。在反演初始模型中,电阻率数据由原始感应电动势数据通过自动迭代法转换而来,深度数据由初始电阻率计算得出。经算例表明,收敛过程中满足拟合差最小的拉格朗日乘子μ值随迭代次数增加而递减;反演过程非常稳定,一般5次迭代至10次迭代就能达到收敛,对低阻层的反演深度和电阻率值都接近真实值,但对高阻层的反演深度和电阻率值误差较大。     

起伏地形三维激电连续介质模型快速反演

为提高起伏地形三维激电数据的反演解释精度和处理效率,笔者开展了连续介质模型快速反演方法研究。针对起伏地形连续电性介质模型,提出了四面体单元电性参数分块线性连续变化的参数化方法。在反演中对模型参数施加光滑和背景约束信息,以提高反演的稳定性和分辨率。通过采用降低反问题的维数、压缩存储线性反演方程各矩阵的元素以及将互换原理和拟牛顿法相结合的方式计算偏导数矩阵等手段,可有效地加快反演的计算速度。最后,为验证反演方法的有效性,对2例地电模型进行反演试算,计算结果表明：反演耗费时间较少,仅迭代6次拟合差便趋于稳定,反演结果能较好地刻画异常体形态,编制的反演解释软件可用于实际生产。     

水平层状介质电磁法数据时频联合反演

通过一维层状介质电磁法数据的时域和频域正演结果,对不同电阻率大小的电阻层进行了反演计算,说明了在这种情况下电磁法数据时频联合反演的合理性.正演计算了一维层状模型下水平电偶极频率域水平电场响应及利用G-S变换计算时间域垂直磁场导数响应;利用最光滑模型约束反演方法建立了时频联合反演迭代公式.通过理论模型反演试算表明,该反演效果较好.     

改进的遗传算法在电测深反演中的应用

遗传算法是地球物理数据非线性反演中常用的方法,但在使用中所需内存较多,容易"早熟"或陷入局部极小。为了在电阻率测深反演中更好地使用遗传算法,提出了利用统计学建模获得精细初始模型后,再利用改进的遗传算法反演。改进了传统遗传算法评价解的方法、增加了局部目标函数并限制遗传基因的繁殖方向。对理论模型的反演计算证明改进的遗传算法可以减少计算量,改善电阻率测深的反演效果。     

电阻率三维反演稳定性和可靠性研究

电阻率三维反演参数太多,将导致反演不稳定性和高度非唯一性等诸多方面的问题。比较传统的最小二乘反演方法,利用共轭梯度的电阻率三维最小构造反演显示了更好的效果。讨论了初始模型选择及数据误差对三维反演结果的影响,说明了反演方法的稳定性和可靠性。      

水下直流电阻率法数值模拟

利用有限单元法对水下电阻率测深进行数值模拟。与地面测量不同的是，水下电极所处的深度、水体本身和水底地形都会给电阻率数据造成较大影响，给数据的解释带来困难。若先通过计算，分析不同电阻率装置在水下测量时的视电阻率与水深的关系，再用有限单元法对水下的三维地电断面进行数值模拟，讨论水底地形对视电阻率的影响．即可用比较法消除地形影响。     

直流电阻率测深中二维反演程序对三维数据的近似解释

直流电阻率测深剖面数据解释最好的方法是二维反演，这是由于它反演速度快，而三维反演则需要很多的实测数据、很长的计算时间。而且通过对几种三维模型测深数据的二维反演结果的研究，发现尽管对无走向的良导矿体，二维反演方法所圈定的异常中心与真实矿体中心有一定的误差，但基本上还是能反映矿体的真实形态和中心位置。当矿体具有一定的走向长度时，则无论对高阻矿体还是良导矿体二维反演的效果都很好。     

三维地电断面激发极化法有限元数值模拟

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

三维井地电阻率有限元数值模拟及反演

油气资源和矿产资源勘探中,常需要较大的勘探深度和高的深部分辨能力,充分利用现有钻孔的三维井地电阻率成像技术能够适应于这种要求。通过设计一个高阻水平板状体模型,利用Ansys有限元软件对模型进行三维井地电阻率法有限元正演计算,并将模拟的响应视电阻率作为反演的输入,结合平滑约束最小二乘反演,分析了异常特征和分布规律,其结论和认识具有一定实用价值。     

大地电磁与地震反射波走时数据二维联合反演

由于地球物理反演存在固有的多解性问题,如何减小多解性,提高地球物理方法反演结果的可靠性和准确性是当今地球物理领域研究的热点。为了减小大地电磁与地震反射波走时资料反演的多解性,在Gallardo和Meju研究的交叉梯度理论的基础上,实现大地电磁与地震反射波走时数据的二维联合反演算法;设计理论地电模型和速度模型进行合成数据的单独反演和联合反演试算,并验证大地电磁和地震反射波走时数据二维联合反演算法的正确性和有效性。反演结果表明,利用交叉梯度方法耦合大地电磁电阻率参数和地震波速度参数能使联合反演过程中的这两个模型在更新的时候相互约束,联合反演结果比单一反演结果更接近真实模型,从而使反演结果更可靠,减少了反演解的多解性。     

井地2.5维DC反演研究

开展了不同观测方式的井地2.5D直流电阻率反演研究。①从2.5维直流电阻率满足的边值问题出发,采用变分原理结合节点线性基函数推导了2.5维井地DC满足的积分弱解形式;②构建了二阶最大平滑稳定泛函的2.5维井地直流电阻率正则化目标函数,采用共轭梯度算法对正则化目标函数进行最优化求解,并采用逐步衰减正则化因子的求解策略来提高反演的稳定性;③设计了均匀半空间模型得到的数值解与解析解对比,电位的相对误差在2%以内,阐述正演算法的正确性和高精度。另外,分析了不同观测装置的2.5维井地直流电阻率异常体特征,并对不同观测方式对2.5维井地DC理论数据进行反演研究。研究结果表明,井中数据的引入提高了2.5DC对纵向探测的分辨率能力,同时提高了2.5DC反演有效性以及准确性。     

电阻率三维复杂结构的快速反演

运用共轭梯度迭代算法解三维反问题中的线性方程组,并结合求Jacobi矩阵 G 的Rodi算法,则每次反演迭代仅需一次正演计算,大大加快了计算速度,实现了直流电阻率三维快速反演。另外,由于避免了存储 G 和 G T G 所需的庞大存储量,以及在三维反演中加入光滑约束,有利于精细网格下的复杂模型反演。      

高密度电法勘探数据三维处理技术及城市勘察中的应用

在以往的高密度数据反演解释过程中,将地质结构简化成二维,通过地电断面图、切片图进行分析.实际上,地质体具有三维地电结构,二维反演难免会受到旁侧效应的影响,二维图件也不能清晰直观地反映地质体的空间分布情况.本文以某陶瓷厂空洞探测为例,利用实测资料构建三维数据场,在反演参数优化的基础上,选择基于最小二乘法的三维反演技术对数据场进行处理,将反演结果利用三维可视化技术显示出来,对采空区位置及空间形态给出结论.研究结果显示,三维处理技术能够从空间上,多手段、多视角分析物探数据,有效地消除了旁侧效应的影响; 物探结果直观清晰地呈现出空洞的分布位置、范围以及空间特征,与钻探结果吻合较好.     

重磁三维自动反演软件系统的分析与设计

由于三维反演计算非常耗时,导致目前的重磁三维自动反演系统的工作方式效率低。为了改善工作效率和提高反演结果的质量,提出实时可视化自动反演方式,设计并开发了一个重磁三维实时可视化反演系统。该系统在三维反演计算的过程中,将中间反演结果同时进行成图可视化,使得工作人员能及时分析和评价中间反演结果并能修改参数,以此来介入和引导复杂的反演过程,快速有效地得到比较理想的反演结果。该系统具有实时性、交互性和引导性特点。对该系统进行了分析和设计,并给出了其中的关键技术与设计方案。