均衡重力异常和地壳表、浅层地质结构

从理论和试验角度论证了利用均衡重力异常研究地壳表、浅层地质结构的可行性,从而指出:所有合理的地壳均衡模型产生的均衡效应几乎是一样的;地壳均衡模型类型和模型参数的变化,仅使均衡改正值在长波长上发生很小变化,而对均衡重力异常的局部异常影响不大;均衡重力异常与地壳表、浅层地质结构具有明显的对应关系。因此,利用均衡重力异常研究地表和地壳浅层地质结构,以选择简单的地壳均衡模型为好。本方法的优点是不受地形格架影响,便于进行广泛对比和定量研究     

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.     

Seismicity in the platform regions of Ukraine in the zones of anomalous electrical conductivity

It is established for the first time that there are several regions in Ukraine, in which the earthquakes occurring within platform territory are correlated to the anomalous conductive structures in the Earth’s crust and upper mantle. These regions are identified as (1) Donbass and the eastern part of the Dnieper-Donetsk Depression (DDD); (2) eastern margin of the Ingulets-Krivoi Rog suture zone in the area of the Krivoi Rog-Kremenchug fault zone; (3) the western part of the Cis-Azov megablock; (4) the western boundary of the Ukrainian Shield and its slope; (5) North Dobruja and Pre-Dobrujan Depression. The reconstructed tree-dimensional (3D) geoelectrical models of the Earth’s crust and upper mantle feature anomalously low values of electric resistivity. The earthquake sources in the platform areas of Ukraine are localized above the top and in the upper parts of the crustal anomalies of electrical conductivity.     

Comparison of gravimetric and seismic constraints on the structure of the Aegean lithosphere in the forearc of the Hellenic subduction zone in the area of Crete

The island of Crete is located in the forearc of the Hellenic subduction zone, where the African lithospheric plate is subducting beneath the Eurasian one. The depth of the plate contact as well as the internal structure of the Aegean plate in the area of Crete have been a matter of debate. In this study, seismic constrains obtained by wide-angle seismic, receiver function and surface wave studies are discussed and compared to a 3D density model of the region.The interface between the Aegean continental lithosphere and the African one is located at a depth of about 50 km below Crete. According to seismic studies, the Aegean lithosphere in the area of Crete is characterised by strong lateral, arc–parallel heterogeneity. An about 30 km thick Aegean crust is found in central Crete with a density of about 2850 kg/m³ for the lower Aegean continental crust and a density of about 3300 kg/m³ for the mantle wedge between the Aegean crust and the African lithosphere. For the deeper crust in the area of western Crete two alternative models have been proposed by seismic studies. One with an about 35 km thick crust and another one with crustal velocities down to the plate contact. A grid search is performed to test the consistency of these models with gravimetric constraints. For western Crete a model with a thick lower Aegean crust and a density of about 2950 kg/m³ is favoured. The inferred density of the lower Aegean crust in the area of Crete correlates well with S-wave velocities obtained by surface wave studies.Based on the 3D density model, the weight of the Aegean lithosphere is estimated along an E–W oriented profile in the area of Crete. Low weights are found for the region of western Crete.     

首都圈上地壳高精度三维P波速度模型——基于石油地震叠加速度和人工地震测深剖面

本文收集了首都圈地区40个测点的石油地震叠加速度资料,经常规处理后得到各测点下方速度随深度变化的曲线;对9条人工地震测深剖面的解释结果进行数字化处理获得各剖面下方离散的速度数据;应用上述资料和专业地质建模软件构建了首都圈地区(115.50°E-117.60°E,38.40°N-40.75°N)范围内上地壳高精度三维P波速度模型.结果表明:华北盆地为隆坳相间区,从东至西依次是黄骅坳陷、沧县隆起和冀中坳陷,上地壳速度结构十分复杂;结晶基底的埋深变化剧烈,冀中坳陷下最深处可达10 km,沿构造走向整体呈西南深、东北浅的趋势,沧县隆起下埋深约2～4 km,黄骅坳陷下最深处则达9 km,剧烈的基底起伏反映出盆地内部不同次级构造单元的差异沉降和中、新生代以来强烈的拉张构造运动.太行山、燕山隆起下的基底埋深较盆地区浅,体现出隆起区新生代以来的抬升构造运动.本文首次将石油地震叠加速度资料用于首都圈地壳速度模型的构建,与以往用人工地震测深资料得到的模型相比,本文结果对华北盆地复杂的上地壳结构刻画得更为细致.     

Comparison of gravimetric and seismic constraints on the structure of the Aegean lithosphere in the forearc of the Hellenic subduction zone in the area of Crete

The island of Crete is located in the forearc of the Hellenic subduction zone, where the African lithospheric plate is subducting beneath the Eurasian one. The depth of the plate contact as well as the internal structure of the Aegean plate in the area of Crete have been a matter of debate. In this study, seismic constrains obtained by wide-angle seismic, receiver function and surface wave studies are discussed and compared to a 3D density model of the region.The interface between the Aegean continental lithosphere and the African one is located at a depth of about 50 km below Crete. According to seismic studies, the Aegean lithosphere in the area of Crete is characterised by strong lateral, arc–parallel heterogeneity. An about 30 km thick Aegean crust is found in central Crete with a density of about 2850 kg/m³ for the lower Aegean continental crust and a density of about 3300 kg/m³ for the mantle wedge between the Aegean crust and the African lithosphere. For the deeper crust in the area of western Crete two alternative models have been proposed by seismic studies. One with an about 35 km thick crust and another one with crustal velocities down to the plate contact. A grid search is performed to test the consistency of these models with gravimetric constraints. For western Crete a model with a thick lower Aegean crust and a density of about 2950 kg/m³ is favoured. The inferred density of the lower Aegean crust in the area of Crete correlates well with S-wave velocities obtained by surface wave studies.Based on the 3D density model, the weight of the Aegean lithosphere is estimated along an E–W oriented profile in the area of Crete. Low weights are found for the region of western Crete.     

Crustal and uppermost mantle structure of SE Tibetan plateau from Rayleigh-wave group-velocity measurements

A shear-wave velocity model of the crust and uppermost mantle beneath the SE Tibetan plateau was derived by inverting Rayleigh-wave group-velocity measurements of periods between 10 and 70 s. Rayleigh-wave group-velocity dispersions along more than 3,000 interstation paths were measured based on analysis of teleseismic waveform data recorded by temporary seismic stations. These observations were then utilized to construct 2D group-velocity maps in the period range of 10–70 s. The new group-velocity maps have an enhanced resolution compared with previous global and regional group-velocity models in this region because of the denser and more uniform data coverage. The lateral resolution across the region is about 0.5° for the periods used in this study. Local dispersion curves were then inverted for a 3D shear-wave velocity model of the region by applying a linear inversion scheme. Our 3D shear-wave model confirms the presence of low-velocity zones (LVZs) in the crust beneath the northern part of this region. Our imaging shows that the upper-middle crustal LVZ beneath the Tengchong region is isolated from these LVZs beneath the eastern and northern part of this region. The upper–middle crustal LVZ may be regarded as evidence of a magma chamber in the crust beneath the Tengchong Volcanoes. Our model also reveals a slow lithospheric structure beneath Tengchong and a fast shield-like mantle beneath the stable Yangtze block.     

Three Dimensional Crustal Density Structure of Central Asia and its Geological Implications

This paper introduces the scale-depth law of multi-scale wavelet analysis for regional gravity data processing, and presents the results of its application to Central Asia for computation of the 3D crustal density structures. The wavelet analysis method is applied for characterizing 3D crustal density structure, producing five maps of density disturbance corresponding to different depths of equivalent layers in the crust. The results provide important evidence for the study of crustal structures and mass movement in Central Asia: (i) the small-scale and intensive linear density disturbances in the upper crust indicate Phanerozoic orogenic belts; (ii) there exists a horseshoe-shaped low-density belt in the middle crust coinciding with the Kazakhstan orocline; (iii) there is a very low density zone in the lower crust, extending from western Kunlun to Tianshan, probably indicating a lower-crust flow; (iv) there are a few lowdensity spots in the middle crust, which might be caused by low-density mass squeezing upward from the lower crust flows.     

3D upper crustal density structure of the Deccan Syneclise,Central India

A constrained 3D density model of the upper crust along a part of the Deccan Syneclise is carried out based on the complete Bouguer anomaly data. Spectral analysis of the complete Bouguer gravity anomaly map of the study region suggests two major sources: short wavelength anomalies (<100 km) caused primarily due to the density inhomogeneities at shallow crustal level and long wavelength anomalies (>100 km) produced due to the sources deeper than the upper crust. A residual map of the short wavelength anomalies is prepared from the complete Bouguer anomaly using Butterworth high‐pass filter (100 km cut‐off wavelength). Utilizing the constraints from deep resistivity sounding, magnetotellurics and deep seismic sounding studies, 2.5D density models have been generated along 39 profiles of this region. The mismatch between the calculated response of the a priori 2.5D model with the residual (short wavelength) gravity anomalies is minimized by introducing high‐density intrusive bodies (≥2.81 g/cm³) in the basement. With these 2.5D density models, the initial geometry of our 3D density model, which includes alluvium, Deccan trap, Mesozoic sediment and high‐density intrusive bodies in the basement up to a depth of 7 km (upper crust), is generated. In the final 3D model, Deccan trap extends from 200 m to nearly 1700 m below the 90–150 m thick Quaternary sediment. Further down, the sub‐trappean Mesozoic sediment is present at a depth range of 600–3000 m followed by the basement. The derived 3D density model also indicates six intrusive bodies of density 2.83 g/cm³ in the basement at an average depth of about 4–7 km that best fits the residual gravity anomaly of the study area.     

Three-dimensional velocity model of the crust of the Bohemian Massif and its effects on seismic tomography of the upper mantle

We have compiled a representative three-dimensional P-velocity model of the crust of the Bohemian Massif (BM) to provide a basis for removing effects of the crustal structure in teleseismic tomography of the upper mantle. The model is primarily based on recently published 2D velocity models from findings of wide-angle refraction and near-vertical reflection seismic profiles of CELEBRATION 2000, ALP 2002, and SUDETES 2003 experiments. The best fitting 3D model of the BM crust (NearNeighbour model) is complemented by velocities according to the reference Earth model at sites where data are sparse, which precludes creating artificial heterogeneities that are products of interpolation method. To test the model, we have performed tomographic inversions of the P-wave travel times measured during the BOHEMA II experiment and compared the results obtained with and without crustal corrections. The tests showed that the presented crustal model decreases magnitudes of velocity perturbations leaking from the crust to the mantle in the western part of the BM. The tomographic images also indicated a highvelocity anomaly in the lower crust or just beneath the crust in the Brunovistullian unit. Such anomaly is not described by our model of the crust since no seismic profile intersects this part of the unit. The tests also indicated that crustal corrections are of the great importance especially for interpretations of the uppermost mantle down to depths of about 100 km.     

One-Dimensional qP-Wave Velocity Model of the Upper Crust for the West Bohemia/Vogtland Earthquake Swarm Region

The western part of the Bohemian Massif (West Bohemia/Vogtland region) is characteristic in the relatively frequent recurrence of intraplate earthquake swarms and in other manifestations of past-to-recent geodynamic activity. In this study we derived 1D anisotropic qP-wave model of the upper crust in the seismogenic West Bohemia/Vogtland region by means of joint inversion of two independent data sets - travel times from controlled shots and arrival times from local earthquakes extracted from the WEBNET seismograms. We derived also simple 1-D P-wave and S-wave isotropic models. Reasons for deriving these models were: (a) only simplified crustal velocity models, homogeneous half-space or 1D isotropic layered models of this region, have been derived up to now and (b) a significant effective anisotropy of the upper crust in the region which was indicated recently by S-wave splitting. Both our anisotropic qP-wave and isotropic P-and S-wave velocity models are constrained by four layers with the constant velocity gradient. Weak anisotropy for P-waves is assumed. The isotropic model is represented by 9 parameters and the anisotropic one is represented by 24 parameters. A new robust and effective optimization algorithm - isometric algorithm - was used for the joint inversion. A two-step inversion algorithm was used. During the first step the isotropic P- and S-wave velocity model was derived. In the second step, it was used as a background model and the parameters of anisotropy were sought. Our 1D models are adequate for the upper crust in the West Bohemia/Vogtland swarm region up to a depth of 15 km. The qP-wave velocity model shows 5% anisotropy, the minimum velocity in the horizontal direction corresponds to an azimuth of 170°. The isotropic model indicates the V_P/V_S ratio variation with depth. The difference between the hypocentre locations based on the derived isotropic and anisotropic models was found to be several hundreds of meters.     

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.     

Density and magnetic models of the lithosphere along CELEBRATION 2000 profile CEL01

The paper presents 2D density and magnetic models of the crust and upper mantle along the DSS line profile of the CELEBRATION 2000 project that crosses the most important geological units in Central Europe. These are the Alps-Carpathians-Pannonian (ALCAPA) region, the SE part of the Paleozoic Platform (PLZ), the Trans-European Suture Zone (TESZ) and a fragment of the SW portion of the East European Craton (EEC). The density and magnetic models were constructed on the basis of a 2D model of P-wave velocity converted into density model, geological data as well as geothermal data and the results of integrated geophysical modelling for the lithosphere-asthenosphere boundary.     

Application of a Global 3D Model to Improve Regional Event Locations

Accurate location of weak seismic events is crucial for monitoring clandestine nuclear tests, for studying local seismic structures, and for assessing possible seismic hazards. Outside of a few regions with dense seismic networks, weak seismic events (with magnitude less than 4) are usually sparsely recorded at epicentral distances less than 20°. Because of lateral variations in crustal and upper mantle structures, observed travel times of seismic phases deviate significantly from predictions based on 1-dimensional (1D) seismic models. Accurately locating weak seismic events remains a difficult task for modern seismology. Perhaps the most promising solution to this problem is the use of a 3-dimensional (3D) model of the Earth. Here we present the results of a validation test in which, using the 3D model SR2002 of the crust and upper mantle and regional phase data alone, we relocate 200 earthquakes and nuclear explosions in Eurasia. The 3D model is constructed using surface wave dispersion data. The event locations using the 3D model are compared with so-called Ground Truth data, either known by non-seismic means or validated by cluster analysis, with location accuracy mostly 5 km or better. Typically, the 3D model reduces the location errors to about half the values attained with the 1D model; i.e., 18 km location errors are reduced to about 9 km. This test indicates that the location of regional events can be significantly improved by using a global 3D model.     

基于沉积过程的数字岩石建模方法研究

本文提出模拟地层沉积及成岩过程的矿物沉积算法,建立数字岩石模型,并通过对比Micro-CT扫描图像和数值模型的局部孔隙度及平均渗流概率函数分布特征,评价建模的准确性.结果表明,由二维扫描提取的粒径信息作为输入参数,模拟矿物沉积过程建模得到的三维数字岩石模型,能够准确重构原始岩心的非均质性及渗流特性,成功应用于泥质砂岩、碳酸盐岩、页岩等存在多矿物或多尺度孔隙的数字岩石建模中.数字岩石物理是正在兴起的重要技术.数字岩石采用超高分辨率先进成像装备,采集和表征微纳尺度岩石结构,在岩石弹性、电性、核磁、渗流特性等数值计算中发挥重要作用.但是,由于三维直接成像在有限视域内难以表征足够的岩石非均质性,提取二维结构统计特征,利用统计或地质过程法重构具有代表性的三维岩石结构成为十分有价值的研究课题,而且,对业界大量存在的岩石薄片及电镜高清二维图像的深度开发应用也具有重要的现实意义.本文发展的新方法,复原沉积过程,较好地解决了孔隙尺度岩石物理定量研究中数值建模与理论计算的技术瓶颈.     

上月壳中的散射引起月震尾波的数值模拟研究

在阿波罗月震记录中普遍存在着强烈持久的尾波信号,这样的波形特征无法用均匀分层月球模型解释.一个普遍被接受的解释是月震尾波由月球浅层结构对月震波的散射引起.我们采用基于交错网格的伪谱和有限差分混合方法模拟研究非均匀上月壳对月震波的散射效应,在此基础上解释月震尾波的形成机制,并估计出上月壳速度扰动的强度.我们发现,在均匀分层模型基础上,进一步考虑上月壳中的非均匀结构对月震波的散射效应,能有效地解释月震信号中强烈持久的尾波.我们认为月震尾波可能是由上月壳中的低波速、低衰减和散射这三个因素的共同作用所引起.采用不同的扰动标准差模拟上月壳的非均匀性,并比较模拟波形与真实月震图的相似程度,我们发现上月壳中速度扰动的标准差应该在3%到5%之间,很可能接近于3%.     

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

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

Joint inversion of teleseismic P waveforms and surface-wave group velocities from ambient seismic noise in the Bohemian Massif

Joint inversion of teleseismic P-waveforms and local group velocities of surface waves retrieved from ambient seismic noise has been performed to model velocity structure of the crust and uppermost mantle of the Bohemian Massif. We analysed P-waveforms of 381 teleseismic earthquakes recorded at 54 broadband seismic stations located on the territory of the Czech Republic and in its close surroundings. Group velocities of Rayleigh and Love surface waves were obtained by cross-correlating long-term recordings of seismic noise. The basis for waveform inversion is the well-known methodology of P-to-S receiver functions constructed from converted phases. Due to instabilities in direct inversion of receiver functions caused by the necessity of applying deconvolution, we propose an alternative formulation to fit observed and calculated radial components of P waveforms. The joint inversion is transformed into a search for the minimum of the cost function defined as a weighted sum of waveform and group velocity misfits. With the use of the robust stochastic optimizer (Differential Evolution Algorithm), neither derivatives nor a starting model are needed. The task was solved for 1D layered isotropic models of the crust and the uppermost mantle. We have performed a sequence of inversions with models containing one, two, three and four layers above a half-space. By using statistical criteria (F-test) we were able to select the simplest velocity models satisfying data and representing local geological structures. Complex crustal models are typical for stations located close to boundaries of major tectonic units. The relatively low average P to S wave-velocity ratio is in agreement with the generally accepted view that the BM crust is predominantly felsic.     

中国深部构造研究的进展

本文扼要地总结了我国对深部构造研究的进展情况。曾经在几个不同地区进行过地震测深的工作,有些地区发现地壳中存在高速或是低速的夹层。曾经利用数字计算求得不同地区的地震面波频散,并用阻尼最小二乘方法反演地壳的参数。根据地震面波频散的特征,可以研究我国大地构造单元的区分。     

基于多震相联合反演江苏地区一维速度模型研究

利用地震走时数据,采用联合反演方法获取了江苏地区的一维P波速度模型。与仅采用初至波走时的传统天然地震走时获取方法相比,该方法充分利用了大量存在的续至波参与反演,能有效改进中下地壳的反演能力。针对地震震相目录中常存在震相标识错误的问题,采用的自动判别筛选震相方法能最大限度提高数据走时的精度,可以对不同震相进行有效区分。与其他常用一维速度模型相比,本文反演的模型对Pg、Pn震相走时拟合效果最佳,残差最小。当所用走时数据拥有较高定位精度时,该反演方法能为研究区三维速度结构成像和地震定位提供较可靠的一维速度模型。