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

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

The upper and middle crustal velocity structure of the northern part of Hebei plain inferred from short period surface wave dispersion

IntroductionThenorthernpartofNorthChinaplainisoneofthemajorregionsforearthquakemonitoringandprediction.AlotofdeepstructureinvestigationhasbeencarriedolltinBeding,Tianjin,TangshanandZhangiiakouinthepastyears,andalotofimportantresultshavebeenachievedby...     

大陆活化克拉通与裂谷地堑地区的地震构造模型

经常观察到这样的现象：活化的克拉通在其演化的不同阶段经历过裂谷作用，在这些地区有分割的地堑构造，有些地堑已变成现代的活跃地震带，如中国山西临汾地堑，中国河北邢台地堑和美国的新马德里地堑。大量资料表明，这些地区的地震成因与许多地质、地球物理条件有关，而不仅仅是活动断层。以临汾地堑为例，地震发生过程涉及到下列因素：（１）地堑周围强度较大的变质岩，这些岩体在区域应力场作用下产生应力积累；（２）地堑内５－     

陕西省数字地震台下方壳幔速度结构研究

利用广义反、透射系数的传递矩阵计算理论地震图，用快速模拟退火法搜索最佳模型，通过与地震台站实际记录的P波波形进行拟合对比，反演了陕西省13个数字地震台下方的P波速度结构。台站分布区域涉及到秦岭造山带、渭河盆地及鄂尔多斯地台三个构造区域。结果显示，3个构造区的地下速度结构差异明显，呈现了各自的构造特点。其中秦岭造山带上地壳乃至地表速度高，中下地壳出现多个低速层，可能与岩石脱水、构造滑脱带等有关；渭河盆地内上中地壳高低速度层交替出现，反映了断层活动在浅部的影响；鄂尔多斯地台有部分台站保持稳定的匀速构造。反映了古老地台地层稳定发展的特性，部分台站受到后期改造，低速层位增多。因为变化剧烈的台站分布在盆地的南部和西部。而且震源机制解显示该区主应力场方向为近东西向，可以说秦岭北缘的断裂作用以及后来来自青藏高原的挤压应力是导致渭河盆地形成和发展的主要因素。最后我们讨论了速度结构与强震、地震分布以及不同构造演化的关系。     

Grounded-source TEM modelling of some deep-seated 3D resistivity structures

Long-offset transient electromagnetics (LOTEM) is now regarded as a suitable electrical method for deep exploration along with magnetotellurics (MT). In this method, the vertical magnetic-field impulse response and, occasionally, the horizontal electric-field step response of a grounded-wire source on the surface of the earth are measured. Here, these two responses are computed for 3D models of three deep resistivity structures of interest in hydrocarbon exploration: (i) a faulted graben in a resistive basement rock at a depth of 4 km beneath a conductive overburden; (ii) a facies change in a resistive layer buried at a depth of 2 km in the conductive overburden above a resistive basement; and (iii) an anticlinal uplift of a resistive layer at a depth of 1 km in the conductive overburden above a resistive basement. The results show that the sensitivity of the electric-field response to model perturbation is generally greater than that of the magnetic-voltage response. Further, the electric-field sensitivity is confined to early and intermediate times while that of the magnetic-voltage response is confined to intermediate and late times. Hence it is recommended that both electric and magnetic recordings are made in a LOTEM survey so that the final results can be presented as apparent-resistivity curves derived from the two responses jointly as well as separately.     

Comparative Magnetotelluric Modeling of Smooth 2D and 3D Conducting Bodies Using Rayleigh-Fourier Codes

—Recently, a method for 3D magnetotelluric modeling was developed, which is based on the application of the Rayleigh scattering theory. This method, RF-3D, is especially capable of modeling multilayered structures with smooth irregular boundaries. The formulation allows inclusion of layers with vertically anisotropic electrical conductivity.¶Using RF-3D, the response of smooth structures of practical interest is calculated and the importance of 3D effects is evaluated. Two models consisting of a 3D conductive body in the presence of a 2D shallow distortion are analyzed. In the first model, the direction of maximum elongation of the body is perpendicular to the strike direction of the 2D upper structure, and in the second one both directions coincide. In addition, the case of a small 3D shallow conductor over a regional 2D structure is also considered.¶3D effects are compared to those generated by 2D models with identical cross sections. In all the cases, the 3D responses differ from those of the 2D, especially directly over the bodies. A good agreement between the 2D transverse magnetic response and the corresponding components of the 3D response, along centrally located transverse profiles, is expected for elongate, prismatic conductors. Then, the differences obtained for the models considered in this study, particularly for the second and third models, are a consequence of the smooth geometry. They can be explained in terms of galvanic effects produced by boundary charges, which are greater near the vertical sides of a prism than on the sides of a body with smooth contours.¶Equivalent 2D models of the first and second structures are also obtained. In these models, the thickness of the conductor is underestimated, respectively, by about 30% and 24%.¶For the third model, when vertical anisotropy is analyzed, it is found that only the anisotropy of the first layer can be detected. This is because the effect of vertical anisotropy decreases strongly with depth and appears to be important only near the 3D anomaly.     

起伏地形下各向异性的2D大地电磁无网格法数值模拟

作为一种基于节点的计算方法,无网格法具有构造高阶导数方便,自适应分析便利的优点,特别适合复杂地质构造的数值模拟.本文针对实际地球物理勘探中存在的起伏地形和各向异性的地电结构,提出用无网格法来模拟大地电磁响应,采用复合二次径向基函数构造了形函数,推导了大地电磁无网格法等价线性方程组,研究了系数矩阵的压缩存储方法以及大型稀疏复线性方程组快速求解算法,实现了起伏地形下各向异性的2D大地电磁高精度数值模拟.基于层状模型验证了算法的正确性,计算结果表明:无论是TM模式还是TE模式,计算相对误差均小于1%;通过对地垒和地堑模型的模拟,得出了起伏地形对视电阻率和相位的影响规律;对起伏地形下含有不同各向异性系数异常体的模型进行了数值模拟,为开展复杂地质结构的电磁场特征研究以及地形校正奠定了理论基础.随着计算科学的快速发展,无网格法必将成为新的高精度电磁场数值模拟方法.     

Selecting the appropriate hydraulic model structure using low-resolution satellite imagery

This study aims at investigating the value added by an explicit representation of floodplain processes in hydraulic models, by comparing the results of a pure 1D model to the ones of an integrated 1D-2D model. The analysis is performed on a 96 km reach of the River Po (Italy). For this test site, two different model structures (1D and 1D-2D) are implemented and compared. The models are first calibrated using high water marks of a high magnitude event and then validated by means of a coarse resolution space-borne flood extent map of a low-magnitude event. The explicit inclusion of 2D features strongly improves the performance of the model, both in terms of accuracy and robustness. The possible causes of this behaviour are investigated and some statements are drawn based on the analysis of the features of the two flood events. Some general conclusions on the effect of the choice between 1D and 2D models arise from the analysis of the calibration-validation procedure. In this context a new measure of fit is proposed, to be used if binary (wet/dry) observations are available and is proven to be able to better distinguish among different sets of model parameters. Lastly, the study also shows the ability of timely low resolution satellite imagery to assist the selection of the most appropriate hydraulic model structure.     

ON THE POSSIBILITY OF THE TELLURIC METHOD: SOME RESULTS ON FAULTED STRUCTURES*

The telluric method of geophysical investigation is re-evaluated and several improvements are suggested. Besides the traditional Jacobian, several new parameters that bear information on the structure under investigation are proposed as sounding tools. These parameters are calculated for several two-dimensional models and plotted as pseudo-sections. The model results, calculated over a wide frequency band (10^-4-10^-2 Hz) suggest that these new parameters can give very detailed information about the subsurface structures. By way of example, we present experimental results obtained along two profiles across familiar structures. The interpretation of the results leads to two-dimensional models of the subsurface that are in good agreement with known geophysical and geological features of these two sections.     

带地形的大地电磁各向异性二维模拟及实例对比分析

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

论华北地区的均衡状态(一)——方法和局部补偿

根据经典均衡的原理,本文重点分析了五种局部均衡补偿模式.在计算方法上利用频域三维重力场的理论公式,快速求得均衡校正值,据此计算了66个模型的均衡重力异常场,注意到对这许多不同参量的模型而言,它们的均衡异常场在形态与分布特征上基本一致,其中计入了地表和上下地壳密度差异分布的 Airy 模型具有最佳的补偿效果,它的均衡面在莫霍界面以下的上地幔中,标准深度50km.从整体上看,华北地台处于亚均衡状态,均衡异常的均值为1810-5m/s2.均衡重力异常的分布表现出明显的块体特征,正均衡异常区主要分布在东部胶辽地块和冀中平原北缘,在汾渭裂谷区存在负异常.模型对比表明,以莫霍界面作为均衡补偿面的模型是不可取的;Airy 模型比 Pratt 模型的补偿效果略好,这同地壳构造以层状为主而侧向变化有限的特征相符.有关复合补偿、均衡重力异常的基本特征和深部构造的关系等结果,将在文章的第二部分发表.      

Quasi-2D hybrid joint inversion of seismic and geoelectric data

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.     

The structure of the crust in TESZ area by kriging interpolation

A precise 3D model of the crust is necessary to start any tectonic or geodynamic interpretation. It is also essential for seismic interpretations of structures lying below as well as for correct analysis of shallow structures using reflection seismics. During the last decades, a number of wide-angle refraction experiments were performed on the territory of central and eastern Europe (POLONAISE’97, CELEBRATION 2000, SUDETES 2003), resulting in many high quality 2D models. It is an interesting and complicated transition zone between Precambrian and Palaeozoic Platforms. This paper presents 3D model of the velocity distribution in the crust and upper mantle interpolated from 2D models of the structure along 33 profiles. The obtained model extends to a depth of 50 km and accurately describes the main features of the crustal structures of Poland and surrounding areas. Different interpolation techniques (Kriging, linear) are compared to assure maximum precision. The final model with estimated uncertainty is an interesting reference of the area for other studies.     

Can Contaminant Transport Models Predict Breakthrough?

A solute breakthrough curve measured during a two-well tracer test was successfully predicted in 1986 using specialized contaminant transport models. Water was injected into a confined, unconsolidated sand aquifer and pumped out 125 feet (38.3 m) away at the same steady rate. The injected water was spiked with bromide for over three days; the outflow concentration was monitored for a month. Based on previous tests, the horizontal hydraulic conductivity of the thick aquifer varied by a factor of seven among 12 layers. Assuming stratified flow with small dispersivities, two research groups accurately predicted breakthrough with three-dimensional (12-layer) models using curvilinear elements following the arc-shaped flowlines in this test.
Can contaminant transport models commonly used in industry, that use rectangular blocks, also reproduce this breakthrough curve? The two-well test was simulated with four MODFLOW-based models, MT3D (FD and HMOC options), MODFLOWT, MOC3D, and MODFLOW-SURFACT.
Using the same 12 layers and small dispersivity used in the successful 1986 simulations, these models fit almost as accurately as the models using curvilinear blocks. Subtle variations in the curves illustrate differences among the codes. Sensitivities of the results to number and size of grid blocks, number of layers, boundary conditions, and values of dispersivity and porosity are briefly presented. The fit between calculated and measured breakthrough curves degenerated as the number of layers and/or grid blocks decreased, reflecting a loss of model predictive power as the level of characterization lessened. Therefore, the breakthrough curve for most field sites can be predicted only qualitatively due to limited characterization of the hydrogeology and contaminant source strength.     

Crust and upper mantle structure of East Asia from ambient noise and earthquake surface wave tomography

The complex tectonic background of East Asia makes it an ideal region for investigating the evolution of the continental lithosphere, for which high-resolution lithospheric structural models are essential. In this study, we measured Rayleigh-wave phase-velocity dispersion curves at periods of 10–120 s and group velocity dispersion curves at periods of 10–140 s using event records from more than 1,000 seismic stations in and around China. By jointly inverting new and previously published dispersion data from ambient noise and earthquakes, we developed a high-resolution shear-wave velocity model down to a depth of ～300 km beneath East Asia. Our model revealed heterogeneous lithospheric structures beneath East Asia, and allowed us to investigate the velocity structure of the entire lithosphere. We also derived crustal and lithospheric thickness models from the three-dimensional (3D) shear-wave model, revealing strong spatial heterogeneity and a general thinning trend of lithospheric thickness from west to east across the study region. Overall, our models reveal important lithospheric features beneath East Asia and provide a valuable baseline dataset for understanding continental-scale dynamics and evolution.     

2-D Modeling of Site Effects Along the EURO-SEISTEST Array (Volvi Graben, Greece)

v--vAn estimation of local site effects in the Volvi basin as derived from observation and modeling is presented in this paper. The Volvi basin is located in the Mygdonian graben in northern Greece near the city of Thessaloniki. This test site has been studied and instrumented in the framework of the "EURO-SEISTEST" and "EURO-SEISMOD" projects funded by the European Union, aimed at improving knowledge of the influence of the local geology on the seismic response of a target area. In this context we calculate synthetic seismograms along a 2-D profile intersecting the graben, instrumented and accurately investigated with a geophysical survey and geotechnical tests. The seismic wavefield from the source to the target area has been computed with the modal summation method, while inside two of the investigated models representing the 2-D section, the wavefield has been numerically propagated with the finite-difference method. We compare the results of the two simulations, both in the time and frequency domain. We also compare the results with experimental data related to an event recorded by the Reftek network installed in the target area. This permits a better understanding of how the structural features of the 2-D models affect the seismic wavefield, especially in the frequency range between 2 Hz and 4 Hz, where one can observe differences between the simulations and the observations. This means that the general features of the models are able to reproduce observed amplification effects, apart from some discrepancies due to still unresolved structural features of the site.     

Isostatic status in North China (1) Method and local compensation

Following Airy and Pratt principles, five kinds of local-compensation models are analysed and a rapid 3-D gravity formula is utilized to calculate isostatic anomalies for 66 models with different parameters. It is noted that isostatic gravity maps appear nearly identical in their main patterns and features. The optimum compensation model in North China is one of modified Airy models in which the different density distribution in the surface, upper crust and lower crust is taken into account and the standard crustal thickness is about 50km. The position of the complete compensation interface is located in the upper mantle. The North China platform as a whole is under sub-isostatic equilibrium status with an isostatic anomaly of about 18·10⁻⁵ m/s² on an average. The distribution of isostatic gravity anomaly shows an obvious blockwise pattern. Most positive anomaly areas occur over the eastern part, the Jiao-Liao Block, Mt. Yan block and northern margin of the Hebei-Shandong block, whereas a negative area occurs in the Shanxi graben. The comparison of models indicates that the Moho discontinuity is not suitable to be taken as a compensation interface, and the compensation effects in Airy model are better than that in Pratt model, which is consistent with the feature of dominant layered structure and less lateral inhomogeneity in crust. Some results about composite compensation, the basic characteristics of isostatic anomaly and deep stucture will be published later in the second part of this paper. Wang Bowen took part in some work in this paper.     

Broadband ground-motion simulation of the 24 May 2014 Gokceada (North Aegean Sea) earthquake (Mw 6.9) in NW Turkey considering local soil effects

On 24 May 2014, a Mw 6.9 earthquake occurred in the west of Gokceada Island, northern Aegean Sea. The earthquake was close to Canakkale, Enez, Tekirdag cities, and damaged 300 buildings in the Marmara Region, NW Turkey. We simulated its broadband (0.1–10 Hz) ground motions including 1D deep and shallow structures soil amplification effects at the 12 strong ground motion stations in the western Marmara Region. The 1D deep velocity structures from the focal layer to the engineering bedrock with an S-wave velocity of 0.78 km/s in different azimuthal directions were tuned by comparing the observed group-velocity dispersion curves of Rayleigh and Love waves from the mainshock with theoretical ones. We also added the shallow parts from previous surveys into the 1D models. Synthetic seismograms on the engineering bedrock were generated using the discrete wave number method with a source model and the 1D deep velocity structures. Then the surface motion was generated considering shallow soil amplification. The synthetic seismograms are generally in good agreement with the observed low and high-frequency parts at most of the stations indicating an appropriateness of the source model and the 1D structural model.     

A STUDY OF THE SEISMIC SIGNATURES OF SEDIMENTATION MODELS USING SYNTHETIC SEISMOGRAMS*

In modern exploration for hydrocarbons there is a great emphasis on the location of stratigraphic traps and estimation of lithologic information like sand-shale ratios from seismic data. In order to investigate the possibilities of success in this endeavour we have studied the synthetic seismograms for wave form and spectral characteristic for four basic sedimentation models: (I) interbedded sand-shale model representing the sediments of generally fluviatile origin, (2) interbedded coal-shale model representing deltaic deposits, (3) sedimentary models representing transgression and regression of shore lines, and (4) a basal sand model. The results have shown that for the first two models a change in the sand-shale or coal-shale ratio results in a characteristically different seismogram. The nature of the seismogram, however, is also strongly dependent on how the sand-shale or coal shale layers are arranged to ultimately give the same number of total layers, thus implying the same coal-shale or sand-shale ratios. The transgression, regression, and basal sand models also produce characteristically different seismic signatures. The spectra of these seismograms show attendant characteristic changes. However, it seems that in the case of real data which are disturbed by noise and the effects of overlying layers these characteristic features may not always be distinguishable.