首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
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.  相似文献   

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

3.
4.
Hydrocarbon exploration in the Papuan fold belt is made extremely difficult by mountainous terrain, equatorial jungle and thick karstified Miocene limestones at the surface. The high-velocity karstified limestones at or near the surface often render the seismic technique useless for imaging the subsurface. In such areas magnetotellurics (MT) provides a valuable capability for mapping subsurface structure. The main structural interface which can be mapped with MT, due to the large electrical contrast, is the contact between the resistive Darai limestone and the underlying conductive sediments of the Ieru Formation. In some areas the base of the Darai can be mapped with reasonable accuracy by fitting 1D models to the observed MT data. However, in many cases where 2D and 3D effects are severe, 1D interpretations can yield dramatically incorrect results. Numerical and field data examples are presented which demonstrate the severity of the 1D errors and the improvements in accuracy which can be achieved using a 2D inverse solution. Two MT lines over adjacent anticlines, both with well control and seismic data, are used to demonstrate the application of 1D and 2D inversions for structural models. In both cases the seismic data provide no aid in the interpretations. The example over the Hides anticline illustrates a situation where 1D inversion of either TE or TM mode provides essentially the same depth to base of Darai as 2D inversion of both TE and TM. Both models provide base Darai depth estimates which are within 10% of that measured in the well. The example over the Angore anticline illustrates the inadequacy of 1D inversion in structurally complex geology complicated by electrical statics. The TE mode fits a 1D Darai thickness of 800 metres while the TM mode fits a 1D Darai thickness of 3500 metres, bracketing the thickness of 2450 metres observed in the well. The final 2D inversion model provides a depth estimate of 2250 metres. Four MT lines along the Angore anticline have been interpreted using 2D inversion. A high degree of correlation exists between lineaments observed on an airborne radar image and zones of low resistivity within the high-resistivity material interpreted as Darai limestone. These low-resistivity zones are interpreted as fault zones. Three-dimensional modelling has been used to simulate 3D statics in an otherwise 2D earth. These data were used to test the Groom-Bailey (GB) decomposition for possible benefits in reducing static effects and estimating geoelectric strike in the Papua New Guinea (PNG) field data. It has been found that the GB decomposition can provide improved regional 2D strike estimates in 3D contaminated data. However, in situations such as PNG, where the regional 2D strike is well established and hence can be fixed, the GB decomposition provides apparent resistivities identical to those simply rotated to strike.  相似文献   

5.
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.  相似文献   

6.
带地形的大地电磁二维正演数值模拟多数基于电性各向同性理论,由于地球内部电性各向异性现象的普遍存在,基于电性各向异性理论研究地形起伏情况下大地电磁二维正演数值模拟就显得非常迫切.本文首先由麦克斯韦方程出发,引入张量电导率,求得一组关于平行走向的电场分量Ex和磁场分量Hx的二阶偏微分方程,使用有限差分法求解出Ex和Hx的近似解,并以此求得其他场分量;其次,引入地形因素,改变变量在网格节点中的排列方式,选择交错排列方式从而给有限差分系数矩阵的最大带宽分配合理的存储空间;最后,使用Weaver的方法解决TM模式下,在地-空分界面垂直于构造走向的一些区域存在不同电导率的问题.通过对带地形的二维电性各向异性结构做正演模拟,研究地形因素对大地电磁响应的影响;以电性各向异性理论为基础,将地形因素引入对实测大地电磁资料的处理中,通过做二维正演拟合和未引入地形因素的结果做对比,说明电性各向异性现象的普遍存在,认识地形因素对观测大地电磁场的影响,为今后分析解释实测大地电磁资料包含地形因素和电性各向异性情况提供理论基础和技术指导.  相似文献   

7.
本文利用大地电磁测深数据,对穿过兰坪—思茅地块和川滇菱形地块以及进入扬子地块的云南兰坪—贵州贵阳大地电磁测深剖面展开了深部电性结构研究.采用大地电磁数据处理分析以及反演技术,对观测资料进行了由定性到定量全面地分析,通过二维非线性共轭梯度反演得到了沿剖面的较为详细的地壳上地幔电性结构,结合其他地质和地球物理资料的分析,对该剖面的二维电性结构进行解释,确定了主要断裂带和边界带的位置以及深部延伸情况,同时确定了壳内低阻层的分布位置,最后进行了区域动力学和孕震构造环境的探讨.研究表明:剖面壳幔电性结构分块性特征与区域地质构造分布特征基本一致,不同地块的电性结构存在显著差异,其中川滇菱形地块的结构相对复杂,上地壳的电性结构为高低阻相间分布特征,电阻率的突变带与地表断裂具有很好的对应关系;兰坪—思茅地块存在中上地壳低阻层,川滇菱形地块中西部存在下地壳低阻层,川滇菱形地块东部和华南地块西部存在中上地壳的低阻层;川滇菱形地块中部攀枝花附近的低阻层埋深最深,而华南地块西部会泽附近的低阻层埋深则最浅;兰坪—思茅地块和川滇菱形地块的中下地壳的低阻层可能与青藏高原物质的东南逃逸有关;华南块体的宣威以东的下地壳不存在低阻层,华南块体下地壳和上地幔的电阻率较高;攀枝花附近的高阻体可能是峨眉山玄武岩喷发导致底侵作用及幔源物质上侵的结果.  相似文献   

8.
The main goal of our study is to investigate 3D topography of the Moho boundary for the area of the northern Red Sea including Gulf of Suez and Gulf of Aqaba. For potential field data inversion we apply a new method of local corrections. The method is efficient and does not require trial-and-error forward modeling. To separate sources of gravity and magnetic field in depth, a method is suggested, based on upward and downward continuation. Both new methods are applied to isolate the contribution of the Moho interface to the total field and to find its 3D topography. At the first stage, we separate near-surface and deeper sources. According to the obtained field of shallow sources a model of the horizontal layer above the depth of 7 km is suggested, which includes a density interface between light sediments and crystalline basement. Its depressions and uplifts correspond to known geological structures. At the next stage, we isolate the effect of very deep sources (below 100 km) and sources outside the area of investigation. After subtracting this field from the total effect of deeper sources, we obtain the contribution of the Moho interface. We make inversion separately for the area of rifts (Red Sea, Gulf of Suez and Gulf of Aqaba) and for the rest of the area. In the rift area we look for the upper boundary of low-density, heated anomalous upper mantle. In the rest of the area the field is satisfied by means of topography for the interface between lower crust and normal upper mantle. Both algorithms are applied also to the magnetic field. The magnetic model of the Moho boundary is in agreement with the gravitational one.  相似文献   

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

10.
Since the true Earth is 3-D in nature, a three-dimensional (3-D) inversion has clear advantages over lower dimensional inversions. We utilized a 3-D magnetotelluric (MT) inversion code, WSINV3DMT, to obtain a realistic resistivity model using a long period MT data set collected in the Northwest Sabalan geothermal field in Ardabil, Iran. The apparent resistivity and phase curves, the magnetic induction vectors, the impedance polar diagrams and the rotational invariant of impedance tensor, indicate a complex 3-D conductivity structure. After setting up the model parameters and designing the appropriate block discretization, we performed the 3-D inversions for two sets of observed data; one set includes the full MT impedance tensor and another set contains only off-diagonal elements of the MT impedance. The final model was selected according to the relative magnitude of the data misfit and the model norm with respect to various Lagrangian multipliers. The results of this study illustrate the 3-D inversion of the off-diagonal elements of MT impedance tensor is precisely enough to explain the structures related to the geothermal source. The obtained results were compared with the results of available 2-D models and they are then interpreted using all of the geological and drilling data of the area. The main outcome of this study is the precise delineation of the geometry of geothermal source that is located at the center of the study area with a surface coverage of about 7 km2.  相似文献   

11.
孙传文  王光锷 《地震地质》1992,14(2):176-182
通过反演受近地表不均匀介质影响的大地电磁测深资料的新方法,可获得较准确的一维地电断面参数、不均匀体埋藏深度及其电性特征。本方法计算速度快,单测ρxy,ρyx可独立反演并相互验证结果的准确性  相似文献   

12.
二维大地电磁资料的偏移   总被引:3,自引:2,他引:3       下载免费PDF全文
本文采用反射映像(即U/D成像)原理对大地电磁资料进行成像,将地面观测的波场分解成上行波和下行波,并分别向下延拓,利用上、下行波的时间一致性原理确定地下电性界面的位置.对TE和TM模式的MT响应分别处理可以得到两个深度剖面,两剖面的一致性便为实际反射界面的最佳估计.理论和实际资料的计算表明,MT偏移技术是比较有效的,它可获得地下界面直观的映像.和常规MT反演方法相比较,该方法具有算法简单等优点,并且还能获得表示地下地质构造真实映像的MT深度剖面.  相似文献   

13.
The coal-bearing strata of the deep Upper Paleozoic in the GS Sag have high hydrocarbon potential. Because of the absence of seismic data, we use electromagnetic (MT) and gravity data jointly to delineate the distribution of deep targets based on well logging and geological data. First, a preliminary geological model is established by using three-dimensional (3D) MT inversion results. Second, using the formation density and gravity anomalies, the preliminary geological model is modified by interactive inversion of the gravity data. Then, we conduct MT-constrained inversion based on the modified model to obtain an optimal geological model until the deviations at all stations are minimized. Finally, the geological model and a seismic profile in the middle of the sag is analysed. We determine that the deep reflections of the seismic profile correspond to the Upper Paleozoic that reaches thickness up to 800 m. The processing of field data suggests that the joint MT–gravity modeling and constrained inversion can reduce the multiple solutions for single geophysical data and thus improve the recognition of deep formations. The MT-constrained inversion is consistent with the geological features in the seismic section. This suggests that the joint MT and gravity modeling and constrained inversion can be used to delineate deep targets in similar basins.  相似文献   

14.
This study deals with two-dimensional (2D) inversions of synthetic and observed radiomagnetotelluric (RMT) data on typical buried conductive waste sites in Europe, and with the practical aspects of different inversion algorithms. In the inversion calculations, we used smoothing and L2-norm stabilizers and compared the results. The resolution of the geometry of the highly conductive waste site, in particular, was investigated. In the inversion with the L2-norm stabilizer, we used the least-squares solution with singular value decomposition (LSSVD) and conjugate gradient (CG), whereas only the conjugate gradient solver was used in the 2D-inversion with the smoothing stabilizer. The inversion results of the synthetic data showed a better resolution of the geometry of the highly conductive waste site when a L2-norm stabilizer was applied in the inversion; in particular, a better detection of the bottom of the waste deposit was achieved. Additional model studies were carried out using synthetic RMT data in order to investigate the 2D inversion of RMT data observed on a 3D structure; these studies showed that the use of TM mode data yields a better resolution of the structure than joint inversion of TE and TM modes.2D inversions of RMT data on a waste site near Cologne showed that the inversion of the TM mode could resolve well the geometry, especially the bottom of the waste site, if information about the background conductivity structure was available. In this case study, inversion with the L2-norm stabilizer produced a sharper image of the waste site than inversion with the smoothing stabilizer, as indicated also by the inversion study that used synthetic data.  相似文献   

15.
青藏高原东缘川滇构造区深部电性结构特征   总被引:4,自引:2,他引:2       下载免费PDF全文
本文对位于青藏高原东缘川滇构造区的贡山一绥江大地电磁测深(MT)剖面数据进行反演,获得沿剖面的深部电性结构,为研究喜马拉雅东构造结、川滇菱形地块与华南地块的构造变形特征、壳幔耦合关系、地块间接触关系以及相互作用等问题,提供电性结构的依据.研究发现:(1)电性结构揭示澜沧江断裂带和小金河断裂带为深大断裂带,控制着研究区的深部结构特征和形变机制;(2)澜沧江断裂带和金沙江断裂带之间的高阻体,可能是扬子古地块的残留部分;小金河断裂带和安宁河断裂带之间的高阻体,则是峨眉山大火山省喷发形成的冕宁一越西杂岩带;(3)在滇西地块、川滇地块和大凉山地块均存在低阻层,它们的介质属性有所不同,滇西地块下的低阻层"疑似"高热状态的岩浆囊,主要由缅甸弧向东俯冲运动引起的,中上地壳的高热状态使地块的活动性增强;川滇地块内部的壳内低阻层的成因为:理塘断裂带和小金河断裂带之间的地表低阻层由破碎带充水所致,而金沙江断裂带和理塘断裂带之间的中地壳低阻层可能是由局部熔融物质或含盐流体导致的,其为壳内物质运移的通道.从而在地下物质发生大规模走滑运动的过程中起到引导作用;川滇地块东部和大凉山地块西部的壳内低阻层可能与地慢物质的上涌有关;马边断裂带附近的低阻体可能与破碎带变宽和破碎带内的流体有关.  相似文献   

16.
《Geofísica Internacional》2014,53(3):259-275
A three-dimensional crustal model for Eastern Cuba, obtained through a process of gravity data inversion is presented. The study area cover a rectangular area of 64 600 km2. The initial model for the inversion was constrained by surface geology, seismic and drilling data. The inversion algorithm uses gravity data to estimate 3-D topographies from the main geological units. The model provides quantitative information on the depths and thicknesses of the geological formations. The resulting model provides new information about the regional composition of the crust. Alien sequences are observed with different compositions and origin over the basement of Bahamas carbonate platform. Most of the maximum gravity anomalies are associated with presence of dense shallow ophiolite sheets. The most remarkable detail is the gravity “southwest” maximum, related to the presence of denser oceanic crust generated in the Cayman spreading center.  相似文献   

17.
重震反演中国东北地壳上地幔三维密度结构   总被引:5,自引:3,他引:2       下载免费PDF全文
本文利用重力和地震P波到时数据反演得到了中国东北地区地壳上地幔三维密度结构.与单一的重力或地震反演相比,重震反演一方面有效地克服了重力反演结果垂向分辨率低的问题,另一方面也提高了地震反演结果的可靠性.结果显示:中国东北地区的地壳及上地幔剩余密度异常分布与构造单元具有明显的相关性,造山带对应低密度异常,盆地对应高密度异常;区域内火山下方有明显的低密度体存在,可能是由于太平洋板块俯冲进入上地幔并部分滞留,在滞留板块深部脱水和软流圈热物质共同作用下产生了上涌岩浆,喷发后形成了火山.  相似文献   

18.
使用阻尼最小二乘法进行震源参数和地壳三维速度结构的走时联合反演.所用资料为S波和P波到时差,并用人工地震资料的二维解释结果作为三维速度模型的特定约束条件.为建立初始模型,又利用天然地震构成了准二维剖面.在走时反演基础上,利用遗传算法进行了几个地震事件的波形反演尝试,并对走时反演获得的地壳速度结构模型的局部进行了修正.以34°~42°N,94°~112°E作为研究区域,在该区域中收集了1986年以来大量地震的S波和P波到时差资料,7条人工地震二维速度剖面资料和2个数字化地震台的几个地震的三分向记录资料.对这些资料进行了处理,最后得出了0~25km深度不同截面的速度分布,并对所得结果进行了分析.  相似文献   

19.
Magnetotelluric (MT) data along two profiles were collected in Tulargen Cu-Ni mining area, Xinjiang, China. The MT profiles, crossing the known ore-bearing magmatic body, are interpreted to have further understanding of the ore-bearing magmatic bodies' stretch and shape in depth, and to evaluate the prospect in the mining area. The values of phase sensitivity skewness indicate a contaminated two-dimensional (2D) structure (there exist local three-dimensional (3D) structures at shallow level) over 3D structures. G-B decomposition of the impedance tensors indicates the dominant geoelectrical strike in survey area is in East-West direction. To obtain reliable electrical structures, four different 2D inversion codes as REBOCC, NLCG, OCCAM, and ABIC are adopted, inversions of a synthetic 3D model using these codes indicate TM data are more suitable to invert the electrical structures along the profiles in Tulargen Cu-Ni mining area, while TE and TE + TM inversions can be used to verify the transverse variation of electrical structures at shallow level. The 2D electrical structures show high-resistivity basement covered by relatively low-resistivity strata in the mining area, the conductive anomalies are imbedded in the upper covers. The known ore-bearing magmatic body is coincided with the positions of two conductive anomalies in both profiles at shallow level; its extension at depth is deduced according to the shapes of the anomalies. Other conductive anomalies, isolated and existing only at limited depth as the known ore-bearing magmatic body, are presumed of good prospect.  相似文献   

20.
穿过郯庐断裂带中段(沂沭断裂带,36°N)所做的大地电磁测深(MT)剖面长约150km.使用Robust技术和远参考道大地电磁方法处理观测数据.通过分析视电阻率、阻抗相位、Swift二维偏离度和区域走向,定性确定测区的电性结构.二维反演解释中选择非线性共轭梯度(NLCG)方法,使用TE、TM两种模式资料联合反演,沿剖面的二维电性结构显示:自西向东,鲁西隆起、郯庐断裂带、胶莱坳陷及鲁东隆起4个电性区块分别对应,鲁东和鲁西隆起区为高阻,郯庐断裂带电性结构复杂,高、低阻相间,胶莱坳陷为低阻(高导)区.沿MT剖面附近曾发生3个地震,其震源区处在电性变化剧烈部位,并在震源区附近存在高导体.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号