共查询到20条相似文献，搜索用时 156 毫秒
1.
《地球物理学进展》2016,(6)
受计算机硬件水平限制,大地电磁(MT)三维反演难以在实际中推广应用,MT数据解释仍以二维反演为主.地质构造具有一定的二维性特征,但不同的地质构造走向往往不同、MT测线无法与每个构造的走向垂直,这使得实际MT数据并不满足严格的二维反演条件.因此,有必要开展大地电磁测深二维反演对"准二维"地质构造的适应性研究.本文设计三个理论地电模型,通过三维正演计算获得各测点大地电磁响应,以此模拟实际观测数据.对MT数据进行相位张量分析,结果表明设计的地电模型主要表现为二维性特征.利用非线性共轭梯度(NLCG)反演算法对理论MT数据进行二维反演研究,重点讨论了测线方向、电性主轴旋转策略及反演模式选择,对反演结果的影响.对比分析反演结果,得到如下认识:1.测线方向对反演结果影响较小;2.电性主轴旋转角度对反演结果影响较大;3.TE+TM联合模式及单独TM模式的反演效果较好.研究结果表明:当一条测线下方在横向与纵向上存在走向不同的多个异常体时,对整条剖面分测点、分频段进行电性主轴旋转,反演所得结果最可靠. 相似文献
2.
传统的一维反演技术已经被广泛应用于航空电磁数据解释中.然而,利用单点水平层状介质模型模拟地下复杂地电结构有时会遇到困难.突出表现在反演参数的横向不连续性,即使相邻测点的反演结果也会出现突变.本文针对航空电磁直升机吊舱系统可进行密集采样,相邻测点地下电性结构应具有某种程度连续性的特点,研究航空电磁数据横向约束反演理论,并提出参数加权约束方法.首先阐述频率域航空电磁正演和加权横向约束反演理论,着重介绍这种拟二维反演方法的基本原理和实施步骤,以及将该方法成功应用于频率域航空电磁数据反演处理的方法技术.最后,通过对理论和实测数据反演处理,并与传统的一维反演结果进行对比,验证加权横向约束反演方法的有效性. 相似文献
3.
4.
航空瞬变电磁法适用于浅表地层的大面积探测,由于实测数据量巨大,一维反演作为一种快速高效的方法被广泛应用.在层状地电条件下,横向约束反演能有效的改善传统单点阻尼最小二乘反演结果中界面横向不连续的情况,但其存在反演过程依赖初始模型,且在复杂层状情况下迭代不稳定的问题.本文提出一种基于OCCAM反演并在正则化项中引入横向光滑约束的光滑拟二维反演方法,在保持OCCAM方法对初始模型依赖性小优势的同时,也使得相邻测点电性横向连续.在反演过程中采用令满足迭代误差的测点不再参与后续迭代的策略,以此减少不必要的正演和雅可比计算,提高反演效率.通过对三个不同地电情况的理论模拟数据进行光滑拟二维反演,并与其他传统反演方法对比,验证了该方法能有效的提高反演稳定性,得到与理论模型高度吻合的反演结果. 相似文献
5.
稀疏测线大地电磁资料三维反演研究:合成算例(英文) 总被引：2，自引：1，他引：1
受勘探成本和工区环境等因素的影响,当前大多数大地电磁实际工作采取布置稀疏测线采集数据和使用二维反演方法解释这些稀疏测线数据的方式。然而,二维反演方法在解释三维地电构造数据时存在局限性,有时甚至做出错误的地质解释。本文尝试了使用三维反演方法对大地电磁稀疏测线数据进行反演解释。使用大地电磁全信息资料三维共轭梯度反演程序对理论模型合成稀疏测线数据进行了三维反演。结果表明:这种反演方案是可行与有效的。同时,我们发现在不同数据的三维反演结果中,四个张量阻抗元素和两个倾子数据同时反演的结果相对更为准确,更接近理论模型。 相似文献
6.
大地电阻率分布信息是影响磁共振地下水探测反演结果准确性的重要因素.在众多电磁法勘探技术中,瞬变电磁法具有高分辨率、高效率和大探测深度等优势,能准确探测地下几百米范围内的电阻率分布信息.因此磁共振与瞬变电磁联合解释方法具有重要意义.然而,利用单一测点拼接的磁共振与瞬变电磁联合解释方法进行模拟二维反演时存在解释结果不唯一,容易出现错误异常体等问题,尤其在复杂地质情况下,同一测线上相邻测点探测结果连续性差,解释结果偏离实际.基于此,本文提出磁共振与瞬变电磁横向约束联合反演方法（Laterally Constrained Inversion,简称LCI）,重点引入外推积分法（quadrature with extrapolation,简称QWE）,解决了传统正演过程中基于直接数值积分方法引起的求解效率低的问题,保证了联合反演方法的顺利实施,进而以相邻测点地下结构应具备连续性为依据,引入横向约束反演思想,通过在联合反演目标函数中加入相邻测点间各模型参数约束矩阵,提高磁共振解释结果准确性,加强探测剖面地质结构和含水模型连续性.经过理论模型证实,本文提出的LCI方法能有效提高传统一维反演结果的稳定性和唯一性.最后,对安徽黄山野外实际探测数据进行横向约束联合反演,验证了磁共振与瞬变电磁LCI联合反演方法的实用性.本文的研究成果将为磁共振与瞬变电磁空间约束联合反演奠定基础. 相似文献
7.
在大地电磁二维反演中,如果网格不均匀或者测点在网格中的位置不恰当,则将出现测点偏离现象,影响反演的精度.此外,由于正演精度和稳定性的需要,反演中网格间距不能差别太大,需要具有一定的光滑性.本文提出了几种自动构建测点中心网格的方法来消除测点偏离现象,并以正则化约束条件来保证网格的光滑程度.本文还提出了一套网格质量检验指标.算例对比研究表明,对于不同的测点集,最佳网格构建的方法可能不相同,但其中“最平模型约束动态配点法”、“最光模型约束动态配点法”拥有最强的适应能力,它们在各种情况下所构建的网格一般同时具备较光滑的网格间距变化和较小的测点偏离程度.此外,采用高质量（光滑性好、测点偏离程度低）网格在一定程度上不仅可提高反演结果的可靠性,还能提高反演的计算效率.本文的方法为获得高质量的大地电磁二维反演网格提供了简单有效的途径. 相似文献
8.
《地球物理学报》2017,(2)
大地电阻率分布信息是影响磁共振地下水探测反演结果准确性的重要因素.在众多电磁法勘探技术中,瞬变电磁法具有高分辨率、高效率和大探测深度等优势,能准确探测地下几百米范围内的电阻率分布信息.因此磁共振与瞬变电磁联合解释方法具有重要意义.然而,利用单一测点拼接的磁共振与瞬变电磁联合解释方法进行模拟二维反演时存在解释结果不唯一,容易出现错误异常体等问题,尤其在复杂地质情况下,同一测线上相邻测点探测结果连续性差,解释结果偏离实际.基于此,本文提出磁共振与瞬变电磁横向约束联合反演方法(Laterally Constrained Inversion,简称LCI),重点引入外推积分法(quadrature with extrapolation,简称QWE),解决了传统正演过程中基于直接数值积分方法引起的求解效率低的问题,保证了联合反演方法的顺利实施,进而以相邻测点地下结构应具备连续性为依据,引入横向约束反演思想,通过在联合反演目标函数中加入相邻测点间各模型参数约束矩阵,提高磁共振解释结果准确性,加强探测剖面地质结构和含水模型连续性.经过理论模型证实,本文提出的LCI方法能有效提高传统一维反演结果的稳定性和唯一性.最后,对安徽黄山野外实际探测数据进行横向约束联合反演,验证了磁共振与瞬变电磁LCI联合反演方法的实用性.本文的研究成果将为磁共振与瞬变电磁空间约束联合反演奠定基础. 相似文献
9.
《地球物理学进展》2016,(3)
航空电磁探测数据量大,二维、三维反演算法复杂、计算速度慢,通常采用一维反演,利用层状模型拼接描述地下复杂结构,但航空电磁数据信噪比低,容易引起一维反演结果横向连续性差等问题.本文针对上述问题,基于一维反演算法,通过整合测线观测数据,建立了测线数据整体的目标函数,并根据Tikhonov正则化反演理论,引入包含空间粗糙度和先验信息的模型参数约束项,确定了拟二维整体反演的目标函数,推导了反演迭代方程组,利用超松弛共轭梯度算法,求得由于整条测线整体反演所致的大型稀疏矩阵的极小化解,实现了对整条测线数据同时反演的固定翼航空电磁数据的拟二维整体反演算法.在反演迭代过程中,正则化因子采用线性搜索自适应迭代的方法自动选取,提高了反演结果的稳定性.对比分析了仿真数据的一维反演与拟二维整体反演结果,得出拟二维整体反演算法横向连续性较好,对高导覆盖层下的导体分辨率优于一维反演,同时受高斯噪声的影响较小.最后,将直升机飞行实测噪声加入仿真数据中,拟二维整体反演结果平均相对误差较一维反演结果降低了31.6%,进一步验证了拟二维整体反演算法的有效性. 相似文献
10.
针对大地电磁二维反演中TE、TM极化模式的选择问题,设计了多个二维、三维模型进行了正演计算,从数据对比和反演结果对比两方面入手,分析研究了三维模型和二维模型响应数据的差异、反演中极化模式的选择、以及三维模型数据的二维反演近似等问题.研究结果表明:在三维模型条件下,利用二维模型进行反演时,TE模式对模型的二维的近似程度要求远高于TM模式;当三维结构影响较明显时,利用TM模式数据进行二维反演比利用TE模式或利用TE+TM模式联合反演都更合理,反演结果中的虚假结构明显减少;对于TM模式,相位受三维畸变影响较小,视电阻率较大,所以二维反演中可适当加大相位的权;对于实测数据的二维反演,应优先考虑采用TM模式数据进行二维反演,其次是TM+TE模式,一般不要单独采用TE模式. 相似文献
11.
对于三维可控源电磁,反演计算效率、张量测量、旁侧效应以及阴影效应是目前研究的热点.本文正演采用基于库仑规范条件的耦合势有限体积算法,反演采用有限内存BFGS(L-BFGS)算法.合成数据反演结果表明:(1)有限内存BFGS法比非线性共轭梯度法,在反演计算效率上具有一定的优势,更适合求解大规模三维可控源电磁反演问题.(2)张量可控源电磁法相对于标量可控源电磁法,前者在模型分辨率上优于后者.(3)在某个区域无法布置测网的情况下,我们可利用旁侧效应在异常体周围布置测网进行三维反演,从而获得真实异常体的信息.同时,为避免阴影效应,我们应在测网外增加可控源电磁控制点,使得三维反演的数据更加完备. 相似文献
12.
系统地论述了大地电磁三维交错采样有限差分数值模拟算法实现过程中交错网格剖分、积分公式离散化、边界条件、方程组求解、三维张量阻抗的计算等内容. 由于提出了简洁的边界条件，采用了解大型系数矩阵方程组的双共轭梯度稳定解法，所实现的三维交错采样有限差分数值模拟算法具有迭代收敛稳定、计算精度高、速度快等特点. 通过两个理论模型的计算结果检验了算法的正确性和计算精度. 所实现的三维交错采样有限差分数值模拟算法为研究三维反演问题奠定了基础. 相似文献
13.
Zhi-Yong Zhang Han-Dong Tan Kun-Peng Wang Chang-Hong Lin Bin Zhang Mao-Bi Xie 《应用地球物理》2016,13(1):13-24
Traditional two-dimensional (2D) complex resistivity forward modeling is based on Poisson’s equation but spectral induced polarization (SIP) data are the coproducts of the induced polarization (IP) and the electromagnetic induction (EMI) effects. This is especially true under high frequencies, where the EMI effect can exceed the IP effect. 2D inversion that only considers the IP effect reduces the reliability of the inversion data. In this paper, we derive differential equations using Maxwell’s equations. With the introduction of the Cole–Cole model, we use the finite-element method to conduct 2D SIP forward modeling that considers the EMI and IP effects simultaneously. The data-space Occam method, in which different constraints to the model smoothness and parametric boundaries are introduced, is then used to simultaneously obtain the four parameters of the Cole—Cole model using multi-array electric field data. This approach not only improves the stability of the inversion but also significantly reduces the solution ambiguity. To improve the computational efficiency, message passing interface programming was used to accelerate the 2D SIP forward modeling and inversion. Synthetic datasets were tested using both serial and parallel algorithms, and the tests suggest that the proposed parallel algorithm is robust and efficient. 相似文献
14.
We present a new, fast and versatile method, the lateral parameter correlation method, of invoking lateral smoothness in model sections of one-dimensional (1D) models. Modern, continuous electrical and electromagnetic methods are capable of recording very large data sets and except for a few cases, standard inversion methodology still relies on 1D models. In environments where the lateral rate of change of resistivity is small, 1D inversion can be justified but model sections of concatenated 1D models do not necessarily display the expected lateral smoothness.
The lateral parameter correlation method has three steps. First, all sounding data are inverted individually. Next, a laterally smooth version of each model parameter, one at a time, is found by solving a simple constrained inversion problem. Identity is postulated between the uncorrelated and correlated parameters and the equations are solved including a model covariance matrix. As a last step, all sounding data are inverted again to produce models that better fit the data, now subject to constraints by including the correlated parameter values as a priori values. Because the method separates the inversion from the correlation it is much faster than methods where the inversion and correlation are solved simultaneously, typically with a factor of 200–500.
Theoretical examples show that the method produces laterally smooth model sections where the main influence comes from the well-determined parameters in such a way that problems with equivalence and poor resolution are alleviated. A field example is presented, demonstrating the improved resolution obtained with the lateral parameter correlation method. The method is very flexible and is capable of coupling models from inversion of different data types and information from boreholes. 相似文献
The lateral parameter correlation method has three steps. First, all sounding data are inverted individually. Next, a laterally smooth version of each model parameter, one at a time, is found by solving a simple constrained inversion problem. Identity is postulated between the uncorrelated and correlated parameters and the equations are solved including a model covariance matrix. As a last step, all sounding data are inverted again to produce models that better fit the data, now subject to constraints by including the correlated parameter values as a priori values. Because the method separates the inversion from the correlation it is much faster than methods where the inversion and correlation are solved simultaneously, typically with a factor of 200–500.
Theoretical examples show that the method produces laterally smooth model sections where the main influence comes from the well-determined parameters in such a way that problems with equivalence and poor resolution are alleviated. A field example is presented, demonstrating the improved resolution obtained with the lateral parameter correlation method. The method is very flexible and is capable of coupling models from inversion of different data types and information from boreholes. 相似文献
15.
Yanling Shi Zhanxiang He Zuzhi Hu Qiang Wei Dechun Li Cuixian Meng Liansheng Ji Sheng Zhang 《Studia Geophysica et Geodaetica》2018,62(1):167-185
Mapping deep geological hydrocarbon targets is of significant importance in basin exploration. In areas lacking reliable seismic data, magnetotelluric (MT) and gravity explorations are helpful to delineate the distribution of potential deep geological hydrocarbon targets. Here we investigate the effectiveness of the integrated 3D MT and gravity explorations for mapping the potential deep hydrocarbon source rocks. The result based on the data from the W Basin (part of the Ordes Basin) of China demonstrates that the method is efficient and economical for basin exploration. The method is particularly useful in target areas which are of great interest for oil and gas exploration but lack high quality seismic data. In our method, we first use the high-precision 3D small-bin MT data acquisition to improve the data accuracy. Then we perform datum static correction method and apply 3D inversion to obtain the3D resistivity distribution. We also develop a layered resistivity model based on resistivity logging to assist the interpretation of the inverted 3D resistivity data so as to derive an initial 3D geological model. Starting from the initial model, we use 2D gravity data to update the model via 2D inversion line by line, and then pass the updated model for the next round of the 3D MT inversion. The integrated inversion is implemented iteratively so the model converges to satisfy the need of final geological analysis. The application to the W Basin shows that we could successfully delineate the geological distribution of the potential deep hydrocarbon source rocks within the basin and map the thickness of the upper Paleozoic. 相似文献
16.
瞬变电磁法的严格三维反演计算复杂、占用资源多,在普通计算机上难以实现.本文引入瞬变电磁矩变换的概念,提出一种快速三维反演方法.该方法基于阻性限制(resistive limit)特性,建立包含异常体的三维大地的一阶矩响应正演算法,根据不同约束条件,选择优化的最速下降法实现瞬变电磁快速三维反演.文中通过含异常体的三维大地正演一阶矩与仿真数据一阶矩的对比,验证了快速三维正演算法的有效性,之后在不同约束条件下,利用优化的最速下降法实现了对含噪声的仿真瞬变电磁数据的快速三维反演.结果表明,该方法能够在普通计算机上短时间内较为准确地反演出地下异常体的体积和位置,在瞬变电磁数据的实时解释工作中具有良好的应用前景. 相似文献
17.
在主流的线性最优化大地电磁二维反演中,如何合理构建初始模型是一个亟待解决的问题.常用的是采用均匀半空间或一维反演结果构建初始模型,不易获得稳定可靠的反演效果.实践表明,尽管基于不同初始模型的大地电磁二维反演结果差别较大,但均较初始模型更为接近真实模型.基于这样一种认识,经过反复的理论和实践探索,我们提出构建大地电磁二维反演初始模型的印模法.印模法的基本思想是依据已有反演结果和均匀半空间模型之间的加权来确定下一步二维反演的初始模型,它一方面保留了已有反演结果中关于真实模型的宏观轮廓信息,另一方面,保证了深部电性结构的均匀性,从而满足大地电磁二维正演所要求的底边界条件.基于印模法,本文进一步提出了迭代重构的反演思想.通过多个理论模型和实测数据的反演计算,验证了上述方法可在很大程度上压制初始模型对反演结果的影响. 相似文献
18.
从二维麦克斯韦方程组出发推导出反演介电常数和电导率等二维介质物性参数的反演公式.反演的步骤是: 建立初始猜测模型，利用电磁波时间域有限差分法模拟正演数据，用正演数据与观测数据之间的数据残差建立目标函数，通过引入一个由麦克斯韦方程计算的伴随场，将目标函数对介质参数的导数表示成显式形式，应用最优化理论得出对初始猜测模型的修改，用共轭梯度法迭代，最终得到反演结果.用合成数据反演具有粗糙地表的非导电介质的介电常数，用实验数据同时反演介电常数和电导率，并比较了麦克斯韦方程反演结果与声波方程反演结果、波动方程偏移剖面的差异. 相似文献
19.
重力数据包含较多的低频信息,重力梯度数据包含较多的高频信息,将重力数据和重力梯度数据进行联合反演得到的结果更加可信.本文基于聚焦反演方法,实现了这一过程.因为联合反演中分量种类增加,所以计算灵敏度矩阵所需要的时间增加,为此,本文提出了一种快速计算灵敏度矩阵的方法.因为联合反演对内存的要求增大,本文选择有限内存BFGS拟牛顿法求解反演问题.本文通过再加权的方法实现深度加权.文中利用单一分量的反演结果来预测异常体的埋深信息,随后将埋深信息结合到深度加权函数中,将其用于多分量组合反演计算.给出了模型试验,发现预测得到的异常体的埋深信息与其实际埋深存在偏差,但是将这一信息应用到反演计算,能够得到与真实模型一致的结果.之后,本文通过模型试验来探究重力和重力梯度联合反演的优势,发现将重力和重力梯度数据联合,能够识别出额外的噪声,反演得到的模型更加合理.但是,对于不同分量组合得到的反演结果是相近的,反演模型的提高很小.最后,将联合反演方法应用到美国路易斯安那州Vinton岩丘的实际数据中,结果显示,将重力和重力梯度数据联合反演,反演模型得到了提高,反演得到的结果与地质资料吻合. 相似文献
20.
Many joint inversion schemes use 1D forward modelling in the integrated interpretation of various geophysical data. In extending the joint inversion approach to the investigation of 2D structures, the discretization of the model parameters and the appropriate choice of the forward‐modelling procedure play a very important role. In this paper, a hybrid seismic–geoelectric joint inversion method is proposed for the investigation of 2D near‐surface geological structures. The electric and seismic models are coupled together through the use of common boundaries between the adjacent layers. Assuming a 2D model composed of homogeneous layers with curved boundaries, a fast ray‐tracing algorithm is used for the calculation of refraction seismic traveltime data. In the geoelectric forward modelling, a locally 1D approximation is used. The boundary surfaces are written in the form of series expansion; the inversion algorithms are formulated for the expansion coefficients and the petrophysical parameters as unknowns. Two versions of the inversion method are proposed: in versions A and B, interval‐wise constant functions and Chebyshev polynomials are, respectively, used as basis functions of the series expansion. The versions are tested by means of synthetic and in situ measured data. The tests show that both methods are stable and accurate. 相似文献