大地电磁局部畸变异常特征三维数值模拟研究

局部畸变问题曾经困扰大地电磁资料反演解释几十年,大地电磁三维数值模拟技术的发展为剖析局部畸变特点和得到可靠的反演成像结果提供了技术基础。本文采用三维数值模拟成像方法对典型三维局部畸变模型进行模拟分析。三维数值模拟结果显示:电场分量垂直电性分界面的极化模式视电阻率曲线(对应二维情况下TM模式)在穿越低阻异常体界面时,曲线会先上移后下移,而在穿越高阻异常体界面时,曲线会先下移后上移,这与电性分界面处积累面电荷产生的二次电场有关。三维模型中XY模式、YX模式视电阻率和相位在三维异常体附近的水平变化是呈现近似垂向对称的,该现象与电场垂直跨越电性界面时视电阻率的变化规律是吻合的,当测线分别沿X方向和Y方向展布时,三维情况下的XY和YX模式分别对应二维情况下的TM模式。低阻小异常体对区域构造响应的畸变影响比高阻小异常体要严重。低阻小异常体对二维区域响应的两种极化模式视电阻率和相位都有非常明显的畸变影响,相比较而言对TE模式的畸变要大于TM模式,因此我们在做二维反演解释时,可优先考虑拟合TM模式数据。位于小异常体中心上方测点的三维畸变响应虽然与对应真实二维区域响应的差异比较大,但可以等效于某种二维模型响应,这种由局部畸变造成的假二维响应在实际野外数据的解释中是需要注意的。      

三维/三维构造下大地电磁相位张量数值模拟

当地表存在三维非均匀电导率分布时,区域大地电磁响应发生畸变. 以往对这种畸变研究多假设近地表为三维,区域构造为一维或二维. 对于更一般的三维/三维构造,为了分析并消除这种畸变影响,真实反映地下三维区域构造信息,本文实现了三维大地电磁相位张量积分方程数值算法,并研究在不同地质模型下相位张量响应. 结果表明,相位张量不仅可以反映一般三维构造信息,亦可有效反映复杂近地表构造下三维区域构造信息,而无须假设区域构造为一维或二维,证明相位张量具有较强抗近地表局部非均匀构造干扰能力,能够保持更为一般的三维区域构造信息. 为了加快正演计算,同时保持一定精度,算法采用了积分方程多网格法.     

一种简单的大地电磁阻抗张量畸变分解方法

提出了一种简单的大地电磁阻抗张量畸变分解方法，引入3d/2d模型，使得全畸变的大地电磁阻抗张量被生意人经，通过对比Swift旋转和Bahr分解后的资料，用联系二的畸变矩阵表示浅层三维电性异常体对区域二维构造的局部畸变影响，通过研究畸变张量矩阵，并对其进行分解，可将局部三维电性不均匀体产生的畸变与区域构造分离开来，从而了解区域构造的维数特征和局部畸变类型及其对区域构造的影响，最后，应用3个测点的实际资料，对这种简单的畸变分解方法进行了检验，结果表明，该方法简单且应用方便，解释效果较好。     

大地电磁三维畸变特征及校正新技术的应用研究

大地电磁研究和应用中,面临的实际地质结构多数是三维的,但是当前三维反演技术还不成熟,仍需要利用二维模型进行近似反演解释。在进行二维反演解释时,如何分析认识三维局部畸变的性质,最大限度地克服三维局部畸变的影响,获得尽可能可信的地质解释结果,如何认识三维模型的二维近似程度,     

基于神经网络的水力压裂监测数据校正方法研究

油气田井况地表复杂,为准确推断水力压裂裂缝方位,提出一种基于BP神经网络的压裂裂缝监测数据校正方法.本文将复杂地表条件下测点布置不均匀时的电位梯度响应,分解为正常场响应和畸变场响应.利用神经网络估算因测点布置不均匀引起的畸变场响应,通过剔除电位梯度响应中的畸变场量,获得测点布置均匀条件下的电位梯度响应,即正常场响应.采用井地电阻率法三维正演模拟裂缝推断过程,结果表明神经网络可有效估计畸变场响应,提高了利用电位梯度响应推断裂缝方位的准确度,具有良好的应用前景.     

大地电磁资料精细处理和二维反演解释技术研究(二)——反演数据极化模式选择

针对大地电磁二维反演中TE、TM极化模式的选择问题,设计了多个二维、三维模型进行了正演计算,从数据对比和反演结果对比两方面入手,分析研究了三维模型和二维模型响应数据的差异、反演中极化模式的选择、以及三维模型数据的二维反演近似等问题.研究结果表明:在三维模型条件下,利用二维模型进行反演时,TE模式对模型的二维的近似程度要求远高于TM模式;当三维结构影响较明显时,利用TM模式数据进行二维反演比利用TE模式或利用TE+TM模式联合反演都更合理,反演结果中的虚假结构明显减少;对于TM模式,相位受三维畸变影响较小,视电阻率较大,所以二维反演中可适当加大相位的权;对于实测数据的二维反演,应优先考虑采用TM模式数据进行二维反演,其次是TM+TE模式,一般不要单独采用TE模式.     

大地电磁Hx型波二维地形改正的方法与效果

用边界无法计算大地电磁H_x型波二维地形影响，并用有限元法计算了起伏地形下水平层状介质的大地电磁响应，用边界元法的结果对畸变后的大地电磁响应进行地形改正，基本恢复了水平地形下层状介质的大地电磁响应的特点．本文给出了物理意义明确、形式简单的视电阻车和阻抗相位地形改正公式．     

三维不规则地形河谷场地地震响应分析方法研究

本文提出了不规则地形地震响应分析的自由场计算方法,并以小湾拱坝场址作为三维场地模型,比较了三维不规则场地侧边界自由场的精确解和近似解;同时,采用大型有限元软件Ansys/Ls-Dyna计算了自由场对场地内部场点地震响应的影响.计算结果表明,三维场地模型的二维自由场近似解与精确解有很大差异,其结果对三维场地内部场点地震响应影响显著,因此,要准确求解三维不规则地形的地震响应必须采用真实的二维侧边界模型.     

瞬变电磁传感器阻尼特性的标定研究

传统瞬变电磁法（TEM）认为临界阻尼匹配时传感器响应快、无振荡,为最佳匹配方式.本文通过建立基于标定环的TEM系统标定模型,理论证明当略欠阻尼匹配时传感器响应最快、时域畸变最小,为最佳匹配方式.本文采用基于实测数据的背景场剔除方法将一次场和大地响应剔除,通过对不同阻尼状态下标定环早期响应畸变的定量评价,实验证明传感器在略欠阻尼匹配状态下二次场响应畸变最小.本文理论和标定实验均证明,在纯二次场观测方式下,采用略欠阻尼匹配方式能够显著降低早期响应畸变、提高系统浅层探测能力,是浅层探测时传感器最佳匹配方式.     

大地电磁测深二维反演方法求解复杂电性结构问题的适应性研究

为探讨二维反演方法在三维电性结构中的适应性问题,本研究中设计了一系列的二维/三维正演模型进行实验计算,分析了三维高阻/低阻异常体对模型响应的畸变作用,并从反演模式选择和数据旋转方向两个方面进行模型二维反演的对比分析,与三维反演的结果进行了比较,最后采用了实测数据进一步进行了二维和三维反演的比较实验.实验和研究结果表明,在剖面选择方面,在剖面方向与垂直主构造方向相差不大的情况下,截取剖面方向,将电性主轴旋转到垂直剖面方向的二维反演结果与垂直主构造方向的反演结果都可以较好地还原正演模型,在大的构造的反映上并无太大差异.在地下为二维或近三维条件时,正演模型的主要结构都可以较好地被二维和三维反演解析出来.二维的反演结果可能甚至会比三维的反演结果的边界更清晰,更精确.然而,对于具有较强的三维结构的模型而言,其二维反演结果与原始模型可能仍然存在较大差异,其中TM+TP或TM模式的二维反演结果相对更接近原始模型,而TE模式的结果往往会有较大误差,需要在解释时特别注意以免得出错误结论.     

海岸效应对近海地区大地电磁测深数据畸变作用研究

在近海地区采集的大地电磁测深数据通常受到海岸效应的影响,使得大地电磁测深数据发生畸变,因而很难利用大地电磁测深资料较为可靠地获得地下深部的电性结构.本文通过正演模拟方法,分析和总结海水深度变化和海底地形变化对近海地区大地电磁测深数据的畸变影响.当测区与海岸线的距离小于目标频率的大地电磁场趋肤深度时,高导海洋的存在会严重影响测区内电磁场的分布.由于海岸效应的影响,大地电磁测深视电阻率曲线和相位曲线均会发生不同程度的畸变,在低频部分,这种畸变作用尤为明显.大地电磁测深一维Occam反演方法和二维非线性共轭梯度反演方法,对近海地区浅部地层具有较好的反演效果.随着海水深度的增加和海底地形的复杂变化,两种反演方法均会出现不同程度的假异常,为地质解释工作造成了影响.近渤海地区的实测大地电磁测深数据在低频部分可能受到海岸效应的影响而导致视电阻率曲线的严重畸变.     

面向目标自适应有限元法的带地形三维大地电磁各向异性正演模拟

传统三维大地电磁各向异性模拟均是基于规则六面体网格,计算精度有限且较难拟合复杂地质条件.本文采用面向目标自适应非结构矢量有限元法,对三维大地电磁各向异性介质进行模拟.首先从电场双旋度方程出发,利用伽辽金方法建立变分方程;然后利用电流密度连续性条件构建适合大地电磁各向异性问题的加权后验误差估计方法,实现面向目标的网格自适应正演;最后通过典型算例分析各向异性对网格自适应和大地电磁响应的影响特征以及各向异性的识别方法.本文算法能够高精度地拟合起伏地表和任意各向异性介质,适用于分析复杂地电条件大地电磁响应特征,为提高大地电磁资料解释水平提供了理论基础.     

APPLYING 3D INVERSION OF SINGLE-PROFILE MAGNETOTELLURIC DATA TO IDENTIFY THE SHADE AND YUNONGXI FAULTS

Many synthetic model studies suggested that the best way to obtain good 3D interpretation results is to distribute the MT sites at a 2D grid array with regular site spacing over the target area. However, MT 3D inversion was very difficult about 10 years ago. A lot of MT data were collected along one profile and then interpreted with 2D inversion. How to apply the state-of-the-art 3D inversion technique to interpret the accumulated mass MT profiles data is an important topic. Some studies on 3D inversion of measured MT profile data suggested that 2D inversions usually had higher resolution for the subsurface than 3D inversions. Meanwhile, they often made their interpretation based on 2D inversion results, and 3D inversion results were only used to evaluate whether the overall resistivity structures were correct. Some researchers thought that 3D inversions could not resolute the local structure well, while 2D inversion results could agree with the surface geologic features much well and interpret the geologic structures easily. But in the present paper, we find that the result of 3D inversion is better than that of 2D inversion in identifying the location of the two local faults, the Shade Fault(SDF)and the Yunongxi Fault(YNXF), and the deep structures. In this paper, we first studied the electrical structure of SDF and YNXF based on a measured magnetotelluric(MT) profile data. Besides, from the point of identifying active faults, we compared the capacity of identifying deep existing faults between 2D inversion models and 3D models with different inversion parameters. The results show that both 2D and 3D inversion of the single-profile data could obtain reasonable and reliable electrical structures on a regional scale. Combining 2D and 3D models, and according to our present data, we find that both SDF and YNXF probably have cut completely the high resistivity layer in the upper crust and extended to the high conductivity layer in the middle crust. In terms of the deep geometry of the faults, at the profile's location, the SDF dips nearly vertically or dips southeast with high dip angle, and the YNXF dips southeast at depth. In addition, according to the results from our measured MT profile, we find that the 3D inversion of single-profile MT data has the capacity of identifying the location and deep geometry of local faults under present computing ability. Finally, this research suggests that appropriate cell size and reasonable smoothing parameters are important factors for the 3D inversion of single-profile MT data, more specifically, too coarse meshes or too large smoothing parameters on horizontal direction of 3D inversion may result in low resolution of 3D inversions that cannot identify the structure of faults. While, for vertical mesh size and data error thresholds, they have limited effect on identifying shallow tectonics as long as their changes are within a reasonable range. 3D inversion results also indicate that, to some extent, adding tippers to the 3D inversion of a MT profile can improve the model's constraint on the deep geometry of the outcropped faults.     

3D magnetotelluric modelling including surface topography

An edge finite‐element method has been applied to compute magnetotelluric (MT) responses to three‐dimensional (3D) earth topography. The finite‐element algorithm uses a single edge shape function at each edge of hexahedral elements, guaranteeing the continuity of the tangential electric field while conserving the continuity of magnetic flux at boundaries. We solve the resulting system of equations using the biconjugate gradient method with a Jacobian preconditioner. The solution gives electric fields parallel to the slope of a surface relief that is often encountered in MT surveys. The algorithm is successfully verified by comparison with other numerical solutions for a 3D‐2 model for comparison of modelling methods for EM induction and a ridge model. We use a 3D trapezoidal‐hill model to investigate 3D topographic effects, which are caused mainly by galvanic effects, not only in the Zxy mode but also in the Zyx mode. If a 3D topography were approximated by a two‐dimensional topography therefore errors occurring in the transverse electric mode would be more serious than those in the transverse magnetic mode.     

利用大地电磁三维反演方法获得二维剖面附近三维电阻率结构的可行性

大地电磁野外实测数据月前大多为二维剖面数据.如何反演这些二维剖面数据获得较为接近实际地电情况的结果,是多数大地电磁工作者关心的问题.我们通过对理论模型的三维响应进行分析和对合成数据及实测资料的反演结果进行对比研究,讨论了利用三维反演的方法来获得大地电磁二维剖面附近三维电阻率结构的可行性.结果表明:可用三维反演的方法来解...     

稀疏测线大地电磁资料三维反演研究:合成算例(英文)

受勘探成本和工区环境等因素的影响,当前大多数大地电磁实际工作采取布置稀疏测线采集数据和使用二维反演方法解释这些稀疏测线数据的方式。然而,二维反演方法在解释三维地电构造数据时存在局限性,有时甚至做出错误的地质解释。本文尝试了使用三维反演方法对大地电磁稀疏测线数据进行反演解释。使用大地电磁全信息资料三维共轭梯度反演程序对理论模型合成稀疏测线数据进行了三维反演。结果表明:这种反演方案是可行与有效的。同时,我们发现在不同数据的三维反演结果中,四个张量阻抗元素和两个倾子数据同时反演的结果相对更为准确,更接近理论模型。     

Mapping deep targets based on integrated 3D MT-gravity interpretation: a case study

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.     

Evaluating MT3DMS for Heat Transport Simulation of Closed Geothermal Systems

Owing to the mathematical similarities between heat and mass transport, the multi-species transport model MT3DMS should be able to simulate heat transport if the effects of buoyancy and changes in viscosity are small. Although in several studies solute models have been successfully applied to simulate heat transport, these studies failed to provide any rigorous test of this approach. In the current study, we carefully evaluate simulations of a single borehole ground source heat pump (GSHP) system in three scenarios: a pure conduction situation, an intermediate case, and a convection-dominated case. Two evaluation approaches are employed: first, MT3DMS heat transport results are compared with analytical solutions. Second, simulations by MT3DMS, which is finite difference, are compared with those by the finite element code FEFLOW and the finite difference code SEAWAT. Both FEFLOW and SEAWAT are designed to simulate heat flow. For each comparison, the computed results are examined based on residual errors. MT3DMS and the analytical solutions compare satisfactorily. MT3DMS and SEAWAT results show very good agreement for all cases. MT3DMS and FEFLOW two-dimensional (2D) and three-dimensional (3D) results show good to very good agreement, except that in 3D there is somewhat deteriorated agreement close to the heat source where the difference in numerical methods is thought to influence the solution. The results suggest that MT3DMS can be successfully applied to simulate GSHP systems, and likely other systems with similar temperature ranges and gradients in saturated porous media.     

Deep structure of the northeastern margin of the Parnaiba Basin, Brazil, from magnetotelluric imaging

The magnetotelluric (MT) method has been applied to the determination of the deep resistivity structure of the northeastern margin of the Parnaiba Basin. Transient electromagnetic (TEM) and MT data were collected in early 1999 along a 95 km long N–S line, extending from the coast across the projected subcrop position of a discontinuous fault found to the west of the study area that is believed to be a possible basin‐bounding fault. The MT data were processed to yield the TE‐ and TM‐mode responses and then corrected for static shift using central‐loop and single‐loop TEM data, respectively. Regularized 2D MT inversion was subsequently undertaken using a structured initial model with the near‐surface constrained by TEM inversion results. As a consistency check, we performed another set of 2D inversions using different smooth initial models. The various optimal 2D inversion models show clearly the presence of a major basement trough, over 2 km deep, located about 70 km from the coast. We interpret it as possibly marking the main basin margin and suggest that it may have implications for groundwater resource development in the area.