台湾地区的壳幔P波速度和Vp／Vs波速比结构研究

利用中国台湾地震台网1992-2004年间记录到的P波和S波到时数据,使用地震层析成像方法联合反演了台湾及邻近地区地壳上地幔的P波速度和Vp/Vs波速比结构。与已有结果相比,采用的球坐标系有限差分算法适合于地壳波速结构强烈不均匀的情况,能够提高走时计算和三维射线追踪的准确程度,而联合反演得到的P波速度、Vp/Vs波速比及其重新定位后的地震分布为揭示壳幔结构的横向变化和深部特征提供了更多约束。结果表明：在地壳浅表层,沉积盆地和造山带具有明显不同的波速特征,沉积盆地P波速度偏低、Vp/Vs波速比偏高,造山带P波速度偏高、Vp/Vs波速比偏低;台东纵谷作为欧亚大陆和菲律宾海板块的碰撞缝合带,Vp/Vs波速比明显偏高,并延伸至中、下地壳和上地幔,反映了岩石破裂、部分熔融和流体活动对壳幔深部物质的影响;台湾东北部琉球岛弧附近的P波速度、Vp/Vs波速比和重新定位后的地震分布勾画出了菲律宾海板块俯冲至欧亚大陆下方的形态,高速的海洋俯冲板片及板片上方P波低速和高Vp/Vs波速比的特征一直延伸至80km;沿着东西向的纵向波速剖面,地壳厚度在中央山脉东侧达到60km左右,并呈现向东倾的趋势,反映出欧亚大陆板块的向东俯冲和挤压使造山带中下地壳物质变形增厚。     

巷道前方煤体波速与瓦斯参数相关分析研究

为了获取掘进巷道前方煤体地震波速度和瓦斯参数间的对应关系,以阳泉五矿为实施矿井.采用矿井地震波速度超前探测方法对掘进煤巷进行探测,通过初至波拾取获得地震波原始记录中纵、横波速,并与巷道前方30m范围内平均瓦斯参数作相关分析.分析结果表明:一定条件下,煤体瓦斯含量越高,其地震波纵、横波速度越低.运用曲线拟合对纵、横波视速度与瓦斯绝对涌出量和瓦斯浓度进行拟合发现:纵、横波视速度与瓦斯绝对涌出量间存在ΔP=-4.6415Vp-2.7048Vs+16.1353的负线性关系;纵、横波视速度与瓦斯浓度间存在η=-0.9302Vp-0.5704Vs+3.3132的负线性关系,并对线性拟合结果进行验证表明:在误差允许范围内,计算值与实测值基本吻合.最后计算得到适用于太原组15#煤层及变质程度更高煤层的视速度瓦斯预警指标:Vp=2.00m/ms和Vs=0.73m/ms.     

基于邻域算法的瑞利面波垂直-水平振幅比及频散曲线联合反演及应用

瑞利面波垂直-水平振幅比(或ZH振幅比)是一个随频率变化的函数,对于台站下方浅层地壳结构非常敏感,且具有和频散资料不同的深度敏感核,是传统频散反演方法的一个很好的补充,从而可以将基阶瑞利面波的ZH振幅比和面波频散数据联合起来更好地反演获得观测台站下方的速度结构.本文提出了基于邻域算法的面波频散曲线与ZH振幅比联合反演方法,我们进行了基于理论模型的模拟测试,证明了联合反演是一种更为可靠的反演方法,且能更好地约束浅层地壳结构.相比于频散曲线单独反演,联合反演不仅可以精确反演获得地壳的Vs结构,对分层地壳的Vp/Vs也能很好地约束.然后我们将联合反演算法应用于实际测量数据,获得了中国西南昆明台(KMI)下方更为准确的地壳横波速度结构及Vp/Vs模型.     

岩石的波速比与静水压的关系

应用高压容器以油为传压介质对岩石样品施加静水压.在加压过程中应用超声脉冲法测量岩石的纵波速度 Vp 和横波速度 Vs 随压力的变化,从而得到岩石的波速度比 Vp/Vs 与静水压的关系.对15种地壳典型岩石的研究表明,当压力超过200MPa 后,由于 vp 和 Vs 随压力变化都很小,Vp/Vs 趋于稳定.当压力在200MPa 以下,虽然 vp 和 Vs 均随压力增加而单调地增加,VP/Vs与压力的关系分为四种类型:(1)Vp/Vs,随压力增加明显地下降;(2)Vp/Vs 随压力增加缓慢地上升:(3)Vp/Vs 随压力增加缓慢地下降;(4)Vp/Vs 随压力增加缓慢地波动.得到的几种岩石在常压下和围压下的波速比的数据,在地壳介质物性参数的研究中具有参考价值.      

长白山天池火山区上地壳三维速度层析成像

利用三维地震透射观测系统所获得的莫霍界面反射波对长白山天池火山区下方上地壳的三维速度细结构进行了反演计算，得到了研究区下方地壳内不同深度的P、S波速度层析成像和VP/VS扰动分布图像，结果表明，天池火山区及其周边速度结构所呈现出的非均匀性变化与该区的地质构造和岩浆系统的存在有着密切的关系，异常变化比较明显的低P波速度异常体分布和高VP/VS扰动体意味这些介质体处在“软化高温”的状态，这与该地区广泛分布的地热异常和低Q值结构的研究结果是一致的.     

上地壳纵横波速度结构相关反演成像方法

基于纵横波初至走时数据的层析成像方法越来越广泛地被应用于揭示不同构造域壳幔速度结构特征.我们从同一地质体的纵横波速度属性相关这一基本思想出发,提出一种相关反演成像的方法:纵横波速度反演交替进行,在迭代反演过程中每通过一次反演获得相应的纵波速度(或横波速度)结构后,更新相应的纵横波速度比模型以及相应的横波(或纵波)速度反演的初始模型,然后继续开展后续横波(或纵波)速度反演工作.在反演过程中依据纵横波速度的相关性信息和射线路径长度将走时残差以不同权重分配到射线路径经过的单元,依据网格节点周围平均的慢度扰动更新速度模型.正反演过程分别基于有限差分走时计算方法和反投影成像方法.两种典型模型试验表明,该技术应用于上地壳速度结构反演成像过程,可有效提高反演结果的可靠性,在很大程度上避免了常规单独反演纵波和横波速度过程容易带来的畸变和失真.该方法应用于重建青藏高原西部札达—泉水沟深地震测深(DSS)剖面下方的上地壳速度结构,揭示出与青藏高原西缘板块碰撞相关的上地壳速度结构特征.     

岩石临破坏前波速奇异变化的再研究

从实验上确证：岩石在主破裂前一瞬间，纵波速度升高甚至超过原最大值的两倍以上，纵、横波速比Vp/Vs升高甚至达2.5～3.0，之所以造成这种现象：对前一个原因是主破裂前有些张裂缝闭合，另一个原因是主破裂前微破裂急剧增多，在某一时刻同时发生几个小破裂，各检波器所记讯息可能不是同一事件。对后，P、S波的走时在破裂孕育过程中的下降--回升（相对变化量）并不同步而是P波的变化先于或大于S波也是重要原因。     

偏离大圆路径传播对四川西部面波相速度成像的影响

面波成像是研究地壳上地幔横波速度结构的一种重要方法．通常,面波相速度或群速度成像假设面波沿大圆路径传播．但是,在地下介质速度结构变化较大时,面波会偏离大圆路径传播,从而导致基于大圆路径假设下的面波成像结果存在一定的误差. 我们采用基于射线追踪的面波成像方法,研究了面波的偏离大圆路径传播对四川西部地区面波相速度成像结果的影响．使用快速行进法(fast marching method)进行面波传播路径的射线追踪,采用子空间反演法(subspace inversion)进行迭代反演,对理论模型合成数据和川西台阵的短周期背景噪声相速度频散数据进行成像分析,并与使用大圆路径传播的成像结果进行对比．对理论模型的测试结果表明,当速度结构变化较大时,基于偏离大圆路径传播的面波成像能够更好地恢复模型异常．对川西台阵的真实数据反演结果显示:在短周期为6 s时,基于偏离大圆路径传播的反演方法较基于大圆路径传播的反演方法所获得的相速度异常的幅度更大些,在四川盆地区域两者的差异接近0.2 km/s;在周期为10 s时,两种反演方法的差异显著减小,基本都在0.1 km/s以内．这主要是因为6 s周期的面波相速度对复杂的上地壳浅层结构更为敏感,从而使得面波的偏离大圆路径传播效应对反演结果的影响更为显著．本文结果表明,当某一周期不同路径的面波相速度测量值变化较大,例如相对于平均相速度的异常超过10％时,则需考虑采用基于偏离大圆路径传播的面波成像方法,否则速度异常较大区域的反演结果可能会造成较大的偏差．      

基于交叉梯度约束的地震初至纵波与瑞雷面波联合反演

本文提出了一种初至纵波(P波)与瑞雷面波的交叉梯度联合反演策略.通过对初至P波进行全波形反演可以获得近地表P波速度结构;通过对仅含瑞雷面波信息的地震数据转换到频率-波数域进行加窗振幅波形反演(Windowed-Amplitude Waveform Inversion,w-AWI)可获得近地表横波(S波)速度结构.在二者反演的目标函数中均加入P波速度和S波速度的交叉梯度作为正则化约束项,使得在反演过程中P波速度和S波速度相互制约,相互约束,从而实现对地震初至P波与瑞雷面波的联合反演.数值模拟结果表明交叉梯度联合反演可以提高S波速度反演分辨率,而P波速度反演结果并没有得到提高.实际资料的反演结果表明,交叉梯度联合反演能够获得更加可信的近地表速度结构.     

基于贝叶斯框架的各向异性页岩储层岩石物理反演技术

页岩岩石物理建模旨在建立页岩矿物组分、微观结构、流体填充与岩石弹性参数的关系.对四川盆地龙马溪组页岩进行岩石物理建模研究,针对页岩黏土含量高、层间微裂缝发育等特点,利用Backus平均理论描述页岩黏土矿物弹性参数,利用Chapman理论计算与水平微裂缝有关的VTI各向异性,并利用Bond变换考虑地层倾角的影响.提出以黏土矿物纵、横波速度和孔隙纵横比为拟合参数进行岩石物理反演的方法,并引入贝叶斯框架减小反演的多解性.由已知的黏土矿物纵、横波速度和孔隙纵横比作为先验信息,并以测井纵、横波速度作为约束条件建立反演的目标函数,同时利用粒子群算法进行最优化搜索.计算结果表明,基于先验约束和粒子群算法的反演方法能够较准确地反演黏土矿物的弹性参数、孔隙形态参数以及裂缝密度等参数.计算得到的黏土纵、横波速度较高,并且在一定范围内变化,这可能与龙马溪组页岩的黏土矿物组分中具有较高弹性模量的伊利石含量较高有关,同时也与黏土定向排列等微观物性特征有关.反演得到的裂缝密度与纵波各向异性参数ε呈明显的正相关,而与横波各向异性参数γ相关性较小.另外,页岩各向异性参数与黏土垂向的纵横波速度有较强的相关性.     

Seismic tomography of Central São Miguel, Azores

We determine the three-dimensional distribution of P- and S-wave velocities for Central São Miguel Island (Azores, Portugal) by tomographic inversion of local earthquake arrival times. We use P- and S-phases from 289 earthquakes recorded by a network of 20 seismometers. The model shows good resolution in the shallowest 5–6 km, as illustrated by different resolution tests. There are several velocity anomalies, interpreted as pyroclastic deposits, intrusive bodies, geothermal fields, and the effects of tectonics. A low Vp zone marks Furnas caldera, probably evidencing volcaniclastic sediments with development of intense geothermal activity. Another low Vp zone extends in correspondence of the highly fractured area between Fogo and the north coast. Conversely, strong positive anomalies are found south of Fogo and northwest of Furnas. They are interpreted in terms of high-density deposits and remnants of a plutonic intrusion. These interpretations are supported by the distribution of Vp/Vs, and are consistent with previous geological, geochemical, and geophysical data.     

The lithosphere and asthenosphere system in Italy as inferred from the Vp and Vs 3D velocity model and Moho map

We present an updated high-resolution tomographic P- and S-wave velocity model of the lithosphere and asthenosphere system in Italy, obtained by adding the observations from ∼25,000 earthquakes (recorded between 2003 and April 2012 at three-component stations) to the previously inverted dataset (165,000 P-wave phases from ∼8000 events recorded between 1988 and 2002 at single-component). The final dataset includes ∼586,000 P- and ∼290,000 S-wave arrival times. The main strength of this research is the use of a method able to model P- and S-seismic phases refracted at the Moho discontinuity. We use a new and original map of the 3D Moho geometry obtained by integrating selected high quality controlled source seismology (CSS) and teleseismic receiver function data. Resolution strongly benefits also from the fast increase in number and quality of INGV National Seismic Network since year 2003 and from its integration with several permanent regional seismic networks. This study confirms the main structural features in the best-resolved parts of the inverted volume and much better images details in the less resolved areas of previous Vp models, due to both the larger number of inverted phases and the more even distribution of seismic stations. A new 3D Vp/Vs model is presented and discussed. Shallow structures and the relationships between the Adriatic, Tyrrhenian, and Adria plates are better imaged.     

Tomographic Imaging of the Seismic Structure Beneath the East Anatolian Plateau, Eastern Turkey

The high level of seismic activity in eastern Turkey is thought to be mainly associated with the continuing collision of the Arabian and Eurasian tectonic plates. The determination of a detailed three-dimensional (3D) structure is crucial for a better understanding of this on-going collision or subduction process; therefore, a body wave tomographic inversion technique was performed on the region. The tomographic inversion used high quality arrival times from earthquakes occurring in the region from 1999 to 2001 recorded by a temporary 29 station broadband IRIS-PASSCAL array operated by research groups from the Universities of Bo?azi?i (Turkey) and Cornell (USA). The data was inverted and consisted of 3,114 P- and 2,298 S-wave arrival times from 252 local events with magnitudes (M _D) ranging from 2.5 to 4.8. The stability and resolution of the results were qualitatively assessed by two synthetic tests: a spike test and checkerboard resolution test and it was found that the models were well resolved for most parts of the imaged domain. The tomographic inversion results reveal significant lateral heterogeneities in the study area to a depth of ~20?km. The P- and S-wave velocity models are consistent with each other and provide evidence for marked heterogeneities in the upper crustal structure beneath eastern Turkey. One of the most important features in the acquired tomographic images is the high velocity anomalies, which are generally parallel to the main tectonic units in the region, existing at shallow depths. This may relate to the existence of ophiolitic units at shallow depths. The other feature is that low velocities are widely dispersed through the 3D structure beneath the region at deeper crustal depths. This feature can be an indicator of the mantle upwelling or support the hypothesis that the Anatolian Plateau is underlain by a partially molten uppermost mantle.     

Prestack density inversion using the Fatti equation constrained by the P- and S-wave impedance and density

Simultaneous prestack inversion is based on the modified Fatti equation and uses the ratio of the P- and S-wave velocity as constraints. We use the relation of P-wave impedance and density (PID) and S-wave impedance and density (SID) to replace the constant Vp/Vs constraint, and we propose the improved constrained Fatti equation to overcome the effect of P-wave impedance on density. We compare the sensitivity of both methods using numerical simulations and conclude that the density inversion sensitivity improves when using the proposed method. In addition, the random conjugate-gradient method is used in the inversion because it is fast and produces global solutions. The use of synthetic and field data suggests that the proposed inversion method is effective in conventional and nonconventional lithologies.     

基于Love波相速度反演南北地震带地壳上地幔结构

收集了南北地震带区域地震台网中292个地震台站2008年1月至2011年3月期间的地震波形数据,由频时分析方法提取了Love波相速度频散曲线,经过反演得到了研究区内的Love波相速度分布.根据Love波纯路径频散,采用线性反演方法对0.25°×0.25°的网格点进行了一维S波速度结构反演,利用线性插值获取了南北地震带地区的三维S波速度结构.结果显示了松潘—甘孜地体和川滇菱形块体地区的下地壳具有明显的S波低速层分布,该异常分布特征支持解释青藏高原隆升及其地壳物质运移的下地壳流模型.在100至120km深度上,川滇菱形块体西北部呈现较强的S波高速异常,这可能是印度岩石圈板块沿喜马拉雅东构造结下插至该区域所致,该区域下地壳的低速软弱物质与上地幔的高速强硬物质形成了鲜明对比,暗示了地壳和上地幔可能具有不同的构造运动和变形方式,这为该区域的壳幔动力学解耦提供了条件.     

Tomographic imaging of P- and S-wave velocity structure beneath Costa Rica

45287 P-wave and 26813 S-wave arrival times from the data base of the Costa Rica network have been tomographically inverted to image the structure beneath Costa Rica. A regularized recursive least squares inverse method was used to produce the high resolution and minimum variance model parameter estimates. The first arrival times are calculated using a finite difference technique, which allows for flexible parameterization of the velocity model and easy inclusion of topography and source-receiver geometry. The P wave velocity structure and hypocenters are determined simultaneously, while the S wave velocity structure is determined using the relocated seismicity and an initial model derived from the P wave model assuming an average P to S wave velocity ratio of 1.78. The most prominent features in the inverted model are a low velocity structure under the volcanic chain in the center of the country, which is related to the hot material connected with the active volcanoes; and a high velocity zone in the mantle, which is related to the Cocos plate subducted under Costa Rica.     

Seismic imaging of the entire arc of Tohoku and Hokkaido in Japan using P-wave, S-wave and sP depth-phase data

In order to better understand seismic structure and seismotectonics of the entire arc of Tohoku and Hokkaido in Japan, we combined arrival time data from earthquakes beneath Tohoku and Hokkaido land areas, and beneath the Pacific Ocean to determine the three-dimensional (3D) velocity structures (Vp and Vs) under the entire Northeast (NE) Japan-Kuril arc. We adopted 176,431 P-wave and 110,953 S-wave arrival times, from 5123 local earthquakes, and 2843 sP depth-phase data from 385 events that occurred beneath the Pacific Ocean. The 385 suboceanic events were accurately relocated by using P-wave, S-wave and sP depth-phase arrival time data jointly. The obtained results confirmed the major features delineated by previous studies and revealed some new features of the structural heterogeneity beneath NE Japan and the Kuril forearcs. High-velocity anomalies of the cold subducting Pacific slab and low-velocity anomalies in the hot mantle wedge were imaged clearly. Strong lateral heterogeneities were revealed on the upper boundary of the Pacific slab under the forearc region, which showed a good correlation with the spatial distribution of large interplate earthquakes. These results indicated that strong coupling sections (or asperities) and weak-coupled or decoupled patches might exist along the upper boundary of the Pacific slab. Widespread low-velocity anomalies were visible in the forearc mantle above the subducting Pacific slab, which might reflect serpentinization of the forearc mantle associated with the dehydration process of the subducting slab. Our results also showed a general tendency for seismic coupling in the asperities to be located around low-velocity areas on the slab boundary under the suboceanic region.     

邢台震源区波速比异常与地震的关系

通过对忻州─泰安人工地震测深剖面Ｐ波、Ｓ波的联合解释，得到沿剖面不同地质单元隆起区与裂陷区、震源区与非震区的速度和波速比结构．鲁西隆起和太行、山西隆起为较均匀的成层构造，地壳厚度分别为３２ｋｍ和４０－４３ｋｍ，波速比为１．７４．中段裂陷区构造变化较大，地壳厚度约３０－３３ｋｍ，波速比为１．７５－１．７７．邢台地震区上地壳下部和中地壳出现高波速比１．７７的异常，与裂陷区东的１．７３形成明显的差异．由此推测，地震的发生不仅与震源区的构造有关，更主要是与震源区岩石的性质有关。     

Merging active and passive data sets in traveltime tomography: The case study of Campi Flegrei caldera (Southern Italy)

We propose a strategy for merging both active and passive data sets in linearized tomographic inversion. We illustrate this in the reconstruction of 3D images of a complex volcanic structure, the Campi Flegrei caldera, located in the vicinity of the city of Naples, southern Italy. The caldera is occasionally the site of significant unrests characterized by large ground uplifts and seismicity. The P and S velocity models of the caldera structure are obtained by a tomographic inversion based on travel times recorded during two distinct experiments. The first data set is composed of 606 earthquakes recorded in 1984 and the second set is composed of recordings for 1528 shots produced during the SERAPIS experiment in 2001. The tomographic inversion is performed using an improved method based on an accurate finite‐difference traveltime computation and a simultaneous inversion of both velocity models and earthquake locations. In order to determine the adequate inversion parameters and relative data weighting factors, we perform massive synthetic simulations allowing one to merge the two types of data optimally. The proper merging provides high resolution velocity models, which allow one to reliably retrieve velocity anomalies over a large part of the tomography area. The obtained images confirm the presence of a high P velocity ring in the southern part of the bay of Pozzuoli and extends its trace inland as compared to previous results. This annular anomaly represents the buried trace of the rim of the Campi Flegrei caldera. Its shape at 1.5 km depth is in good agreement with the location of hydrothermalized lava inferred by gravimetric data modelling. The Vp/Vs model confirms the presence of two characteristic features. At about 1 km depth a very high Vp/Vs anomaly is observed below the town of Pozzuoli and is interpreted as due to the presence of rocks that contain fluids in the liquid phase. A low Vp/Vs body extending at about 3–4 km depth below a large part of the caldera is interpreted as the top of formations that are enriched in gas under supercritical conditions.