大地电磁二维反演方法对比研究

由反演方程出发, 从构造目标函数和方程求解方面, 比较了国际上通用的OCCAM (奥可姆法)、RRI(快速松弛法)和REBOCC(简化基奥可姆法), 并由两个合成模型分析比较了OCCAM、RRI和REBOCC的反演结果。试验表明, REBOCC是继RRI之后又一种较快的反演方法, 可为电法勘探资料处理提供新的途径。而在RRI和REBOCC反演时, 网格的划分对结果很有影响, 这是值得注意的地方。     

大地电磁三维模型二维反演计算数值模拟分析

设计了多个1D/3D、2D/3D模型进行数值模拟分析,从模型三维响应和二维响应的差异大小和反演结果对比两个方面,分析了利用二维反演解释方法来逼近解释三维地质模型的适用性问题.分析结果表明:在异常体沿走向长度逐渐缩短的过程中,模型的三维特性逐渐增强,采用TM极化模式的二维反演结果基本保持不变,TE极化模式的二维反演结果逐渐产生畸变,表明TM极化模式对模型二维性的要求远远低于TE极化模式,TE/TM联合模式二维反演效果通常介于前两者之间,有时优于TM极化模式的二维反演结果.综合得出:一般条件下采用TM极化模式进行二维反演通常可以较好地重建原始模型,TE/TM联合二维反演结果在与TM极化模式反演结果整体电性特征相似的情况下,其对模型电阻率值的反演优于TM极化模式.     

二维大地电磁主轴各向异性理论模型反演

系统阐述了广泛认可的二维有限差分正演算法理论,再考虑一种特殊情况,将任意各向异性情况的电导率简化为只有对角线的主轴电导率,得到电场E x与磁场H x解耦的方程组,并利用TE模式只与σxx有关、TM模式只与σyy和σzz有关的特点,使用最小二乘法将TM模式单独反演计算.在y轴方向上,反演结果较为准确地还原了各向异性异常体...     

起伏地形大地电磁二维反演

为了适应地形起伏的实际地质情况,开展带地形的最小二乘二维反演研究。鉴于大地电磁（ MT）反演的不适定问题,引入Tikhonov的正则化方法,从而获得关于总目标函数的方程,利用光滑约束最小二乘法求解总目标函数方程。由于正则化因子值与反演精度以及稳定性相关,采用主动约束平衡方法获取最优化的正则化因子,以确保反演精度和稳定性都达到最佳。与此同时,利用电磁场互易定理以节省反演迭代过程求解雅可比矩阵的计算时间。构建了若干地质构造模型进行试算,分别讨论TE、TM模式以及二者联合模式的反演结果,并与前人研究工作对比以说明本文方法的反演效果。     

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

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

井地大地电磁非线性共轭梯度二维反演研究

本文提出并验证了地下大地电磁测深数据的正、反演算法和程序,为井地大地电磁测深和井地可控源音频大地电磁测深这一前沿研究课题提供了地球物理解释方法.使其更加完善和有效.在Rodi(2001)发表的大地电磁场二维NLCG反演理论和算法的基础上,推导了井地大地电磁场二维NLCG反演算法,改编反演的代码,开发出二维NLCG反演程...     

电性可控源音频大地电磁法二维反演方法技术研究

在多年来电性可控源音频大地电磁法方法实际应用的基础上,对电性可控源音频大地电磁法二维反演方法技术进行研究。用美国ZONGE公司的SCS2D软件对一个应用实例资料的处理,给出不同的计算结果以及一般的反演计算工作流程。     

在大地电磁二维Occam反演中求取拉格朗日乘子方法改进

Occam反演具有可靠性和加上各种先验信息可行性的优点,但是因为求取拉格朗日乘子需要很长时间,所以前人提出了几个直接选取方法。笔者改进了二维大地电磁测深Occam反演中求取拉格朗日乘子的方法。在搜索过程中加上观测数据拟合差与其期望值的比较。这样,在Occam反演的拟合差下降阶段的最后迭代和模型光滑阶段,能够排除多余的正演计算。这种改进对反演结果没有影响,但却使得反演速度加快并且更符合于Occam思想。根据模型实验和野外数据处理,可以减少正演次数20%~50%。     

3维电磁跨孔反演

有效利用电磁法圈定和监测油气储，要求在其他方法的数据，臂如象地震数据和测井记录约束下的大的跨孔数据组的快速计算机化反演。理想条件下，一种反演算法可以对３维模型数据反演而不需要对模型的几何或电阻率分布有理论上限制。     

3—D结构上大地电磁数据的2—D反演

大地电磁法及其发展－音频大地电磁法与可控源音频大地电磁法,已广泛应用于矿产,石油及地热资源勘探中。尽管所测量的结构是三维的,但通常所应用的仅是二维或一维解释。在本项工作中,我们探讨在三维结构以大地电磁合成数据的三维反演适用性。     

通用河网二维水流模拟模式研究

通过对流域河网二维水流模拟问题的分析,提出了一个通用化的河网二维模拟概化模式,该模式将河网二维概化为"树状"河道计算单元、"环状"河道计算单元、"十字型"河道计算单元等河网二维基本单元,对这些河道计算单元的求解,构建流域型河网一二维耦合模式的求解模型,并由此建立了通用化的河网二维水流数学模型软件系统。通过"树状"河道二维计算单元求解模型的推导与实现,验证了概化模式的可行性,表明该思路可以应用于大型流域河网一、二维水流模拟,具有很好的应用推广价值。     

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

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

CSAMT一维全区反演

由汉克尔变换和连分式求解水平层状大地电偶源的水平电场和垂直磁场,然后由电场和磁场求得卡尼亚视电阻率,不进行近区和中间区校正,直接用阻尼最小法对近区和中间区发生畸变的卡尼亚视电阻率反演,得到地下电性参数。模型试算结果表明:虽然近区和中间区的视电阻率不能正确反映地下电性垂向变化,但由其反演仍能得到正确的结果。     

电导率分块连续变化的二维MT有限元模拟（Ⅰ）

用有限元法进行电导率分块连续变化的二维大地电磁的模拟。本文第一部分给了有限元法系数矩阵的计算公式，可用于编制有限元法的计算程度。模型上的计算结果和计算精度的分析交城本文的第二部分介绍。     

多沙河流平面二维泥沙数学模型研究

针对冲积河流河床演变速度快、边界调整幅度大的特点,采用积分方法数值求解水流泥沙基本方程,建立了适合于多沙河流泥沙及河床演变规律的二维数学模型,在数值方法方面考虑了物质守恒性及对多变边界的适应性,采用了比较简单的动边界处理模式。通过工程实例的计算验证,证实了模型所采用的基本方法是合理的,模型有较强的河床边界适应能力,可用以解决有关多沙河流河床演变预测及浑水洪水预报的有关问题。     

平面二维水流泥沙数值模拟

基于三角形网格划分,采用有限元方法,建立了河道平面二维水流泥沙数学模型。采用质量集中的简化处理和预估校正的时间推进算法,较好地解决了有限元计算存储量和计算速度问题。以空腔流和突扩段两种情况为例对模型进行了检验计算,结果表明,模型能较好地模拟河道水流泥沙运动及河床变形情况,且计算稳定性好、速度快、精度较高。     

基于球体层状介质模型的大地电磁正演

目前大地电磁测深法的一维正演理论,是基于平面波垂直入射水平层状介质的假设模型。但由于地球是一个球体,因此有必要研究基于球状介质模型的大地电磁正演理论。这里详细推导了基于球体层状介质模型的大地电磁正演公式,计算了若干理论模型。通过同基于水平层状介质模型的大地电磁正演结果对比,验证了正演公式的正确性。同时,指出当探测周期增加到上万秒时,阻抗相位会增大;而当探测周期增加到数十万秒时,视电阻率会减小。