塔里木盆地结晶基底的反射地震调查

大型克拉通内部沉积盆地基底组构与盆地起源和油气聚集有密切关系.除了钻井岩心及周缘露头研究可提供盆地结晶基底类型证据外,反射地震信号记录了沉积盆地起源时期的有关大地构造作用信息,盆地基底组构可通过记录长度大的反射地震剖面研究.2007年中石化在塔里木盆地将1400 km的地震剖面接收记录从6 s加长到12 s,为研究克拉通盆地结晶基底的组构和类型提供了难得的第一手资料.这篇文章主要介绍这次调查的反射地震剖面,讨论深反射地震数据处理的关键技术,展示塔里木盆地巴楚—塔中地区的四条12 s反射地震剖面,并对反射剖面的散射模式作初步分析.本次调查表明,将地震剖面接收记录从6 s加长到12 s,采集处理的成本只增加了3％左右,但是可为研究克拉通盆地结晶基底和上地壳不均匀性提供很有价值的第一手资料.     

基于蒙特卡罗法的三源锥束CT散射分析

在医用X-CT中,散射现象对重建图像的质量有重要的影响。多源锥束CT是医用CT的一种发展趋势,可成倍提高时间分辨率。但与单源CT相比,多源CT额外的光源会导致散射现象的加重,因此在实际应用中需要着重考虑这一问题。本文在开源工具EGSnrc的基础上,开发了一套基于蒙特卡罗法的CT仿真软件SimProj,用该软件获得了贴近...     

Numerical inversion of deformation caused by pressure sources: application to Mount Etna (Italy)

The interpretation of geodetic data in volcanic areas is usually based on analytical deformation models. Although numerical finite element (FE) modelling allows realistic features such as topography and crustal heterogeneities to be included, the technique is not computationally convenient for solving inverse problems using classical methods. In this paper, we develop a general tool to perform inversions of geodetic data by means of 3-D FE models. The forward model is a library of numerical displacement solutions, where each entry of the library is the surface displacement due to a single stress component applied to an element of the grid. The final solution is a weighted combination of the six stress components applied to a single element-source. The pre-computed forward models are implemented in a global search algorithm, followed by an appraisal of the sampled solutions. After providing extended testing, we apply the method to model the 1993–1997 inflation phase at Mt Etna, documented by GPS and EDM measurements. We consider four different forward libraries, computed in models characterized by homogeneous/heterogeneous medium and flat/topographic free surface. Our results suggest that the elastic heterogeneities of the medium can significantly alter the position of the inferred source, while the topography has minor effect.     

Crustal structure of Iraq from receiver functions and surface wave dispersion: implications for understanding the deformation history of the Arabian–Eurasian collision

We report the crustal structure for two locations in Iraq estimated by joint inversion of P -wave receiver functions (RFs) and surface (Rayleigh) wave group velocity dispersion. RFs were computed from teleseismic recordings at two temporary broad-band seismic stations located in Mosul (MSL) in the Zagros Fold Belt and Baghdad (BHD) in the Mesopotamian Foredeep. Group velocity dispersion curves at the sites were derived from continental-scale tomography. The inversion results show that the crustal thicknesses are 39 km at MSL and 43 km at BHD. We observe a strong Ps _Moho at BHD consistent with a sharp Moho discontinuity. However, at MSL we observe a weak Ps _Moho suggesting a transitional Moho where crustal thickening is likely to be occurring in the deep crust. Both sites reveal low velocity surface layers consistent with sedimentary thickness of about 3 km at station MSL and 7 km at BHD and agreeing well with the previous reports. Ignoring the sediments, the crystalline crustal velocities and thicknesses are remarkably similar at both stations. The similarity of crustal structure suggests that the crust of the northeastern proto-Arabian Platform was uniform before subsidence and deposition of the sediments in the Cenozoic. If crystalline crustal structure is uniform across the northern Arabian Platform then crustal thickness variations in the Zagros Fold Belt and Thrust Zone should reveal the history of deformation and crustal shortening in the Arabian–Eurasian collision zone and not reflect pre-existing crustal thickness variations in the Arabian Plate.     

New seismogenic source and deep structures revealed by the 1999 Chia-yi earthquake sequence in southwestern Taiwan

In a tectonically active setting large earthquakes are always threats; however, they may also be useful in elucidating the subsurface geology. Instrumentally recorded seismicity is, therefore, widely utilized to extend our knowledge into the deeper crust, especially where basement is involved. It is because the earthquakes are triggered by underground stress changes that usually corresponding to the framework of geological structures. Hidden faults, therefore, can be recognized and their extension as well as orientation can be estimated. Both above are of relevance for assessment on seismic hazard of a region, since the active faults are supposed to be re-activated and cause large earthquakes. In this study, we analysed the 1999 October 22 earthquake sequence that occurred in southwestern Taiwan. Two major seismicity clusters were identified with spatial distribution between depths of 10 and 16 km. One cluster is nearly vertical and striking 032°, corresponding to the strike-slip Meishan fault (MSF) that generated the 1906 surface rupture. Another cluster strikes 190° and dips 64° to the west, which is interpreted as west-vergent reverse fault, in contrast to previous expectation of east vergence. Our analysis of the focal solutions of all the larger earthquakes in the 1999 sequence with the 3-D distribution of all the earthquakes over the period 1990–2004 allows us reinterpret the structural framework and suggest previously unreognized seismogenic sources in this area. We accordingly suggest: (1) multiple detachment faults are present in southwestern Taiwan coastal plain and (2) additional seismogenic sources consist of tear faults and backthrust faults in addition to sources associated with west-vergent fold-and-thrust belt.     

On crustal corrections in surface wave tomography

Mantle models from surface waves rely on good crustal corrections. We investigated how far ray theoretical and finite frequency approximations can predict crustal corrections for fundamental mode surface waves. Using a spectral element method, we calculated synthetic seismograms in transversely isotropic PREM and in the 3-D crustal model Crust2.0 on top of PREM, and measured the corresponding time-shifts as a function of period. We then applied phase corrections to the PREM seismograms using ray theory and finite frequency theory with exact local phase velocity perturbations from Crust2.0 and looked at the residual time-shifts. After crustal corrections, residuals fall within the uncertainty of measured phase velocities for periods longer than 60 and 80 s for Rayleigh and Love waves, respectively. Rayleigh and Love waves are affected in a highly non-linear way by the crustal type. Oceanic crust affects Love waves stronger, while Rayleigh waves change most in continental crust. As a consequence, we find that the imperfect crustal corrections could have a large impact on our inferences of radial anisotropy. If we want to map anisotropy correctly, we should invert simultaneously for mantle and crust. The latter can only be achieved by using perturbation theory from a good 3-D starting model, or implementing full non-linearity from a 1-D starting model.     

Surface wave array tomography in SE Tibet from ambient seismic noise and two-station analysis – II. Crustal and upper-mantle structure

We determine the 3-D shear wave speed variations in the crust and upper mantle in the southeastern borderland of the Tibetan Plateau, SW China, with data from 25 temporary broad-band stations and one permanent station. Interstation Rayleigh wave (phase velocity) dispersion curves were obtained at periods from 10 to 50 s from empirical Green's function (EGF) derived from (ambient noise) interferometry and from 20 to 150 s from traditional two-station (TS) analysis. Here, we use these measurements to construct phase velocity maps (from 10 to 150 s, using the average interstation dispersion from the EGF and TS methods between 20 and 50 s) and estimate from them (with the Neighbourhood Algorithm) the 3-D wave speed variations and their uncertainty. The crust structure, parametrized in three layers, can be well resolved with a horizontal resolution about of 100 km or less. Because of the possible effect of mechanically weak layers on regional deformation, of particular interest is the existence and geometry of low (shear) velocity layers (LVLs). In some regions prominent LVLs occur in the middle crust, in others they may appear in the lower crust. In some cases the lateral transition of shear wave speed coincides with major fault zones. The spatial variation in strength and depth of crustal LVLs suggests that the 3-D geometry of weak layers is complex and that unhindered crustal flow over large regions may not occur. Consideration of such complexity may be the key to a better understanding of relative block motion and patterns of seismicity.