P wave velocity structure of upper and middle crust beneath the Three Gorges reservoir dam and adjacent regions

The three-dimensional P wave velocity structures of upper and middle crust beneath the Three Gorges reservoir dam and adjacent regions were reconstructed by means of the seismic tomography, using the absolute P-arrivals from the events observed by the Three Gorges seismic network. Our preliminary study of the deep structure shows: (1) the lithology of the Huangling anticline, integral and highly solidified, exhibits a uniformity. Down to a depth of 20 km, the Huangling anticline still shows high-velocity anomalies. Note that the characteristics of the velocity images with the different depths are not the same; (2) the depth of the sedimentary rocks beneath the entire Zigui basin is not constant, which is less than 6 km; (3) the top boundary of the crystalline basement varies from 0 to 7 km, whereas the bottom one between 14 and 16 km; (4) the partially molten substances of the upper mantle are located in the Badong-Zigui and its adjacent regions as well as nearby Peishi area. The development and evolution of the current tectonics are associated with the cooling and solidifying process of the molten substances within the upper mantle; (5) the different velocity structures of the deep crust in both east and west side of the Jiuwanxi fault are to be investigated in order to understand the origin of gravitational gradient zone in the Three Gorge areas; and (6) the velocities within the upper crust would be affected by reservoir water permeability. Supported by National Natural Science Foundation of China (Grant No. 40574039), National Key Technology R& D Program (Grant No. 2008BAC38B02)     

巧家地震和鲁甸地震震源区周边P波速度结构

本文利用2010年1月至2020年6月巧家地震和鲁甸地震震源区周围发生地震事件的走时观测资料,应用双差层析成像方法获得了2014年MS6.5鲁甸地震和2020年5月18日巧家MS5.0地震周边区域中上地壳的P波速度结构.成像结果显示:整个研究区域的速度结构存在很强的非均匀性:在鲁甸地震震源区附近,浅部存在速度高达6.4...     

Searching for Regional Crustal Velocity-Density Relations with the Use of 2-D Gravity Modelling – Central Europe Case

In this paper we search for a reference relation between seismic P-wave velocity V and density ρ _ref for the continental crust. Based on the results of modern seismic experiments, we compiled 2-D seismic models into a network of four, each about 1100–1400 km long, continental-scale seismic transects cutting all main tectonic units in Central Europe. The Moho depth (about 52 km beneath the TESZ in SE Poland, to about 25 km beneath the Pannonian Basin) and the crustal structure of this area are characterised by a large variation. This structural variation provides an interesting basis for gravity studies and especially for analysing the difference of the density structure between two major tectonic provinces of distinctive age difference: Precambrian and Phanerozoic. The 2-D gravity modelling applied for crustal cross-sections representing the regional structure, based on a unified gravity anomaly map of the area, allows for a stable determination of some general features of the regional reference velocity-density relation for the continental crust. In general three major seismo-petrological types of rocks can be distinguished: sediments, crystalline crust and mantle. In compacted sediments the reference velocity-density relation is well described by the Gardner or Nafe-Drake model. Calculated gravity anomalies, using unified velocity-density relation for the whole crystalline crust, well describe observed anomalies, with an average difference of 14 mGal. However, calculated gravity anomalies, using separated velocity-density relations for the crystalline crust of Precambrian and Phanerozoic Europe, describe observed anomalies better than for the entire crust, with an average difference 12 mGal. The most important feature of these relations is the large differentiation of the derivative dρ _ref /dV in the crystalline crust, being about 0.3 g s/m⁴ for Precambrian, and about 0.1 g s/m⁴ for the Phanerozoic crystalline crust. The modelling suggests a very small density value in the uppermost mantle ρ = 3.11 g/cm³ below the younger area, while for the older area it is ρ = 3.3 g/cm³.     

华北北部地区地壳三维P波速度结构的双差地震层析成像

利用区域固定台站和华北科学探测台阵记录的10 461个近震事件的183 909个Pg波绝对走时和495 753个相对走时数据,采用双差地震层析成像获得华北北部(37.5°～41.5°N,111.5°～119.5°E)范围内的地壳三维P波速度结构模型。结果表明:研究区内各主要构造单元具有明显不同的速度结构特征,速度异常的走向与区域构造的走向一致,浅层速度图像很好地反映了地表地质和岩性的变化;重定位后的大部分地震集中在0～20km的深度上,主要位于低速区的内部或高速和低速区的交界部位;三河—平谷和唐山地震震源区中、下地壳的低速异常可能是流体的显示。结合前人成果和本文模型所揭示的深、浅结构,我们认为太平洋板块在中国东部之下的俯冲和滞留引起板块脱水、软流圈物质上涌等一系列过程,软流圈热物质到达上地幔顶部并沿超壳断裂上侵进入地壳,致使上地幔顶部和下地壳中的含水矿物发生脱水作用产生流体,流体继续上移造成中、上地壳发震层的弱化,从而导致大地震的发生。因此华北北部地区的强震活动,以及地壳结构的非均匀性应是与太平洋板块俯冲、滞留引起的深部过程密切相关的。     

利用区域双差层析成像方法研究川滇南部地壳结构特征

本文联合使用云南、四川和贵州地震台网的85个地震台站在2008年1月—2017年12月期间记录的49130个地震、317366个初至Pg震相绝对到时数据和2674110条高精度的相对到时数据,采用区域双差地震层析成像方法联合反演了川滇南部地壳三维P波速度结构和39621个地震的震源参数,探究了川滇南部中下地壳流和腾冲火山区岩浆囊的分布特征.研究结果表明：（1）川滇南部上地壳的速度异常特征与地表地形密切相关;（2）小江断裂带的中下地壳存在一条绵延近二百多公里的低速异常结构,最南端受到红河断裂带的阻挡而终止于断裂带南段北侧,这可能是川滇南部的一条中下地壳流,低速异常结构在红河断裂带南段转而向南东流动反映了红河断裂带可能为川滇菱形块体的西南边界;（3）红河断裂带各段速度异常存在明显的差异,重定位后的震源分布显示红河断裂带中段和南段虽然不如北段地震活动强烈,但地震震源深度分布较北段深;（4）腾冲火山区西侧和北侧下方10~20 km深度范围内存在的低速异常体推测为通过怒江断裂带形成的岩浆通道从中地壳涌入上地壳的岩浆囊,可能反映了自更新世延续至今的以橄榄玄武岩和安山岩为主要岩性的壳内岩浆活动,持续的岩浆活动为地表热活动提供了主要动力.     

四川盆地东西陆块中下地壳结构存在差异

四川盆地是中上扬子克拉通的主要组成部分.作为我国三大稳定克拉通之一,扬子克拉通经历了自太古代以来的长期演化,直到新元古代晚期与华夏板块发生碰撞拼合前,一直被认为是一个稳定的统一陆块.基底包括了新太古宙-新元古代岩层,其上广泛被新元古代晚期至显生宙地层覆盖,仅有~2.9—2.95Ga基底岩石零星出露于四川盆地的西缘、西南缘和三峡地区,使得对于沉积盖层之下的中下地壳的性质和分布规模的认识十分有限.重力异常则能够宏观揭示区域结构特征.本文通过刨除沉积盖层和莫霍面起伏引起的重力异常而获得了中下地壳的重力异常,反映了四川盆地东西陆块中下地壳存在结构差异,结合深地震反射资料、航磁异常和地球化学资料,证实了该分界线位于重庆—华蓥一线,故而推测中上扬子克拉通在太古宙-古元古代可能存在东西两个陆核.     

四川芦山MS7.0地震震源区及其周边区域P波三维速度结构研究

利用四川数字地震台网和流动地震台站在芦山M_S7.0地震震后(2013年4月20日—6月23日)记录到的2026次区域地震事件的28188条P波到时资料,采用地震层析成像方法反演得到了芦山地震震源区及其周边区域中上地壳P波三维速度结构. 结果表明,浅部地壳的P波速度异常分布特征与地表地质构造、 地形和岩性密切相关,即成都断陷盆地表现出与第四纪沉积有关的低速异常区;犍为、 乐山一带的川中微升区和川青块体龙门山以西的邻近地带均表现为与构造抬升有关的高速异常;宝兴、 康定附近分布的基性火山岩及火山碎屑岩均呈局部高速异常分布. 芦山地震震源位于高低速异常分界线附近且偏向高速体一侧,其下方存在明显的低速异常分布,可能与流体的存在有关. 流体的作用导致中上地壳内部发震层的弱化,使孕震断层易于破裂,可能对芦山地震起到了触发作用. 芦山地震与汶川地震两次地震的余震密集区相距50 km,这50 km地震空区震源体的深度范围附近目前正处于高速异常区内,加之龙门山断裂带西南段又具有比较典型的断错地貌发育,使得该段地震空区(大邑—邛崃活动断裂破裂空段)现在所处的深浅部构造环境变得复杂,其潜在的地震危险性仍值得进一步关注.      

基于背景噪声的晋冀鲁豫交界地区瑞利波群速度层析成像

晋冀鲁豫交界地区位于华北克拉通的东南部, 研究该区地壳结构有助于进一步了解河北南部及其邻区地质构造演化特征、 发震机制、 发震规律。 本文收集了晋冀鲁豫交界地区两年84个地震台连续波形观测资料, 通过台站对波形互相关, 得到各台站对间的经验格林函数, 再利用图像分析技术提取瑞利波频散曲线。 本文基于面波直接成像方法得到不同深度剪切波速度结构。 结果显示, 该区部分速度异常与该区地表地质构造单元存在相关性, 速度异常呈区域分布特征, 山区和盆地呈现非常显著的速度差异, 浅层上地壳结构速度变化特征与沉积层厚度、 结晶基底埋深等区域构造密切相关, 中下地壳及上地幔存在不规则的速度异常, 与该区的地震发生存在相关性。     

小江断裂带周边地区三维P波速度结构及其构造意义

作为青藏高原的东南边界,小江断裂带在高原物质的侧向逃逸中发挥着重要的作用.本文利用流动地震台阵及固定台站的走时观测资料,对小江断裂带及周边区域的壳幔三维P波速度结构进行了研究.结果表明,在中上地壳,小江断裂带内部主要为低速异常,其东侧主要为高速异常.在中下地壳,小江断裂带中部为低速异常,北部和南部主要为高速异常,其中北部的高速异常可延伸到地表附近,南部的高速异常可一直延伸到上地幔.我们推测,小江断裂带中部的低速异常与深部热作用有关;北部的高速异常可能是晚古生代地幔柱活动导致大量基性和超基性幔源物质侵入地壳引起的,它的存在对青藏高原物质向南逃逸起到了一定的阻挡作用,可能是导致川滇活动块体北部次级块体快速抬升的重要因素;南部顶界面向北倾斜的高速异常体对川滇活动块体向南滑移起到了进一步的阻挡作用,导致其上覆的中上地壳低速异常区发生较强的变形和强烈的地震活动,同时在上地幔深度范围起到了稳定的作用,使其南部区域的介质受青藏高原物质向南挤出的影响明显减小.     

Rock physics modeling of heterogeneous carbonatereservoirs： porosity estimation and hydrocarbon detection

In heterogeneous natural gas reservoirs, gas is generally present as small patch-like pockets embedded in the water-saturated host matrix. This type of heterogeneity, alsocalled ＂patchy saturation＂, causes significant seismic velocity dispersion and attenuation. Toestablish the relation between seismic response and type of fluids, we designed a rock physicsmodel for carbonates. First, we performed CT scanning and analysis of the fluid distributionin the partially saturated rocks. Then, we predicted the quantitative relation between the waveresponse at different frequency ranges and the basic lithological properties and pore fluids.A rock physics template was constructed based on thin section analysis of pore structuresand seismic inversion. This approach was applied to the limestone gas reservoirs of the rightbank block of the Amu Darya River. Based on poststack wave impedance and prestack elasticparameter inversions, the seismic data were used to estimate rock porosity and gas saturation.The model results were in ~ood a~reement with the production regime of the wells.     

The application of well logging and seismic modeling to assess the degree of gas saturation in Miocene strata (Carpathian Foredeep, Poland)

The presence of natural gas in the pore space of reservoir rocks results in a significant decrease in P-wave velocity. Even low gas saturation can generate seismic anomalies (DHI) and false image of gas accumulation of economic importance. This article presents an attempt to evaluate gas saturation from 2D seismic section in the Miocene sandstone strata in the south-eastern part of the Carpathian Foredeep. The ESTYMACJA program and the Biot-Gassmann model were used to study the dependence between elastic parameters and saturating fluids (water and gas) recorded in wells. Series of calculations was carried out using a method of fluid substitution for various gas saturation. The applicability of seismic data for evaluating gas saturation of reservoir beds was assessed with the use of 1D modelling (synthetic seismograms) and 2D modelling (theoretical seismic section) calculated for different gas saturation. The proposed methodology can be used to identify low and high gas-saturated zones and contour the reservoir.     

岩石孔隙结构的地震信号特征

岩石孔隙结构是控制砂岩和碳酸盐岩的地震波速度和渗透率的重要参数之一。如果两种类似的岩石其孔隙度一定,而渗透率不同,那么它们的声波速度相差2km/s,而渗透率两者则可能相差近6个数量级,即从0.01mD到20mO。在本文中我们总结了由一个广义孔隙弹性理论简化的一个双参数弹性速度模型，以描述孔隙结构对弹性波的影响。由于矿物和储层的流体是给定的,我们用孔隙度和骨架的柔性(挠性)因子来确定速度模型,这个模型可以用于地震反演和储层表征,已改善孔隙度和储量的计算骨架柔性因子可用于岩石结构(PST)类型的定量分类,并可以利用叠前、叠后的两种地震资料将其与孔隙的连通性和渗透率联系起来。本项研究同时也有助于说明为什么振幅与偏移距(AVO)分析用于流体检测在某些情况下失败。这是由于孔隙结构对地震波的影响能够掩盖所有流体效应，特别是在碳酸岩中。     

上海地区地壳精细结构的综合地球物理探测研究

通过在上海地区开展深、浅地震反射、地震宽角反射/折射、高分辨地震折射和大地电磁测深等联合剖面探测, 获得了该地区近地表至Moho面的精细速度结构、电性结构和深浅构造关系.结果表明, 该地区地壳可划分为上、中、下三个组成部分.其中,上地壳厚为12～14 km,波速为57～59 km/s;中地壳厚度约为10 km,波速为59～62 km/s;下地壳厚为10～11 km, 波速为62～63 km/s,Moho面深度约为32 km.剖面浅部地质构造复杂,共解释出12条特征明显的断裂.其中,除3条断裂错断结晶基底（G界面）并向下延伸至上地壳底界面外,其他断裂均在深度3～5 km以上终止或收敛于G界面之上.此外,仅在剖面西侧基底下部约13～15 km埋深处发现一厚度在2 km左右的壳内高导层.所以,在综合各方面资料后分析认为,在剖面经过地区不存在发生大地震的深部构造条件,近地表所存在的活动断层是未来产生对该区有影响地震的震源区.     

A conceptual model of the copper–porphyry ore formation based on joint analysis of deep 3D geophysical models: Sorskoe complex (Russia) case study

Joint analysis of deep three-dimensional models of the electrical resistivity, seismic velocity, and density of the complex hosting the Sorskoe Cu–Mo deposit (Russia) is carried out aimed at finding geophysical markers characterizing the areas of ore generation, transportation and deposition. The three-dimensional lithology model of the study area is built based on the empirical relationship between the silica content of the rocks and seismic velocities. It is in agreement with geological and geochemical studies provided in this area earlier and could be used as a basis for forecasting locations of the copper–molybdenum ore deposits at depth. A conceptual model of the copper–porphyry complex explaining the mechanisms of ore generation, transportation from the lower to the upper crust and deposition in the upper crust is suggested. In particular, it is supposed that post-magmatic supercritical gas–water ore-bearing fluids are upwelling through the plastic crust due to the sliding of the fluid films along the cleavage planes of the foliated rocks while at the depths of the brittle upper crust this mechanism could be changed by volumetric fluid transportation along the network of large pores and cracks.     

Inversion of single-station teleseismic P-wave polarization data for the velocity structure of Beijing area

Broadband three-component seismic data recorded by Beijing station (BJI) of CDSN were used to calculate P-wave polarization of teleseismic events. These polarization data were then used in the inversion for the underground structure around the Beijing station, especially for the details of velocity discontinuities. The result shows that a conspicuous low velocity zone exists in the crust on the west of the station, which is in good agreement with previous studies. It proves the theory that polarization data could be applied to inversion for velocity structures, especially for boundaries with large velocity gradient. It also demonstrates the feasibility of velocity structure inversion with polarization data from high-quality broadband data recorded by a single station. Therefore, travel-times and polarization data can be jointly used to study velocity structure. Polarization data are more suitable for delineating the boundary of velocity anomalies. Moreover, if the polarization method is combined with receiver function method to fully exploit their complementarity, it is possible to obtain the lateral velocity variation around the station as well as the detailed vertical variation below the station.     

Rayleigh wave phase velocities of South China block and its adjacent areas

Using records of continuous seismic waveforms from 609 broadband seismic stations in the South China Block and its adjacent areas in 2010–2012, empirical Green's functions of surface waves were obtained from cross-correlation functions of ambient noise data between these stations. High quality phase velocity dispersion curves of Rayleigh waves were obtained using time-frequency analysis. These interstation dispersion curves were then inverted to build Rayleigh wave phase velocity maps at periods of 6–50 s. The results of phase velocity maps indicate that phase velocities at 6–10 s periods are correlated with the geological features in the upper crust. Major basins and small-scale grabens and basins display slow velocity anomalies; while most of the orogenic belts and the fold belts display high velocity anomalies. With the gravity gradient zone along Taihang Mountain to Wuling Mountain as the boundary for the phase velocity maps at period of 20–30 s, the western area mainly displays low velocity anomalies, while the eastern side shows high velocity anomalies. Phase velocities in the eastern South China Block south to the Qinling-Dabie orogenic belt is higher than that in the eastern North China Block to the north, which is possibly due to the differences of tectonic mechanisms between the North China Craton and the South China Block. The phase velocities at periods of40–50 s are possibly related to the lateral variations of the velocity structure in the lower crust and upper mantle: The low-velocity anomalies in the eastern part of the Tibetan Plateau are caused by the thick crust; while the Sichuan Basin and the southern part of the Ordos Basin display distinct high-velocity anomalies, reflecting the stable features of the lithosphere in these blocks. The lateral variation pattern of phase velocities in the southern part of the South China Block is not consistent with the surface trace of the block boundary in the eastern Yunnan Province and its vicinities. The phase velocities in the Sichuan Basin are overall slow at short periods and gradually increase with period from the central part to the edge of the basin, indicating the features of shallower basement in the center and overall stable lithospheric mantle of the basin. The middle and upper crust of the southern Ordos Basin in the North China Block is heterogeneous, while in lower crust and the uppermost mantle the phase velocities mainly exhibit high anomalies. High-velocity anomalies are widespread at the middle of the Qinling-Dabie orogenic belt, as well as the areas in southeastern Guangxi with Caledonian granite explosion, but its detailed mechanism is still unclear.     

大别造山带上部地壳结构的有限差分层析成像

将有限差分层析反演方法用于大别造山带人工地震测深剖面上部地壳初至波的走时层析成像，得到这一剖面上部地壳横向不均匀结构图像表明，在南大别地区的超高压变质带发现榴辉岩相岩石，其下方3km深度以内的基底仍具有大陆地壳构造中正常的结晶基底速度（6．00km／s左右）；北大别地区在整个上地壳保持正常的速度，而在超高压带下方3km深度以下为6．20-630km／s的相对高速异常区，这一现象可能与超高压变质岩的含量增大有关．同时还表明，南大别和北大别之间至少在上地壳已有明显差异，它们之间的水吼一五河断裂可能是大别造山带内部的主要构造分界线．     

华北地壳岩石波速类型及其地质意义

根据结晶岩的矿物成分、成因类型与高压室温下波速绝对值的关系,以及高温高压下岩石波速与物相变化的关系,把岩石分为７种波速类型。在一定的构造条件下,根据岩石波速随地壳深度的变化关系,提出了岩石级数的定量概念,它反映了岩石的酸性（或基性）程度,并可根据地震测深值获得相应深度下综合岩石的平均级数。在此基础上,对华北地壳的组成,角闪岩在壳内的生成和存在条件,地壳内石英岩的相和壳内低速层的成因等问题进行了讨论     

Influence of fluid composition on the elastic properties of sandstone and quartzite at high temperature and pressure in application to the problem of crustal waveguides

Recent geophysical studies have revealed zones with anomalously low seismic velocities and increased electric conductivity in the middle part of the Earth’s crust. It is suggested that these zones can probably be related with the variations in the porosity and permeability of the rocks and the presence of fluids, and that they can explain the existence of crustal waveguides. In order to verify this model, we carried out an experimental study. Among the studied metamorphic processes, the acid metasomatose most completely explains the silicification of rocks. The present work considers the physico-chemical causes of the redeposition of quartz from the solutions, the influence of the pH of the solutions on the mass transfer and deposition of quartz during silicification, and alteration of the physical properties of the rocks. Silicification leads to the increase in density, change in microtexture, and cementation of rocks. It also causes the change in seismic velocities in the Earth’s crust.     

汶川地震破裂带断层岩纵波速度与孔隙度关系的实验研究

分布于地震破裂带上的断层岩具有高孔隙度的特征.该特点造成了其弹性波速度与结晶岩石和沉积岩存在明显的差异.确定断层岩的弹性波速度与孔隙度和矿物组成的关系对于利用地震资料探测深部断层和测井资料的解释至关重要.在10~600 MPa条件下,本文对地震断层岩的纵波波速（V_p）和总孔隙度（φ_t）进行了测量,并深入分析了V_p与孔隙度的关系.结果表明在10~600 MPa的压力范围内,V_p（p）随着压力的增高呈现对数增加,其增长率随着压力的上升而逐渐减小,遵从∂V_p（p）/∂p=a_v/p的变化规律.断层岩中的孔隙度随着压力的增高呈对数减小.与传统的认识不同,实验发现在压力高达600 MPa以上,大多数断层岩中仍然可以残留可观的孔隙量.分析显示V_p与总孔隙度及总粘土含量呈负线性相关.该发现有助于认识深部流体的活动通道特征,有助于理解断层带中存在大量粘土矿物、断层带内的物质可被大量带出、围陷波的形成等地质和地球物理现象.