联合利用走时与波形反演技术研究地壳三维速度结构(Ⅱ)——应用

使用阻尼最小二乘法进行震源参数和地壳三维速度结构的走时联合反演．所用资料为Ｓ波和Ｐ波到时差,并用人工地震资料的二维解释结果作为三维速度模型的特定约束条件．为建立初始模型,又利用天然地震构成了准二维剖面．在走时反演基础上,利用遗传算法进行了几个地震事件的波形反演尝试,并对走时反演获得的地壳速度结构模型的局部进行了修正．以３４°～４２°Ｎ,９４°～１１２°Ｅ作为研究区域,在该区域中收集了１９８６年以来大量地震的Ｓ波和Ｐ波到时差资料,７条人工地震二维速度剖面资料和２个数字化地震台的几个地震的三分向记录资料．对这些资料进行了处理,最后得出了０～２５ｋｍ深度不同截面的速度分布,并对所得结果进行了分析．     

TIME-TO-DEPTH MIGRATION USING WAVEFRONT CURVATURE*

A well-known technique for the migration of normal-incidence two-way travel-time maps is extended to common-source-point travel-time data. The travel time and the travel-time gradient are used to compute the parameters defining the tangent plane of the reflecting interface. It is also shown how the curvature matrix of the received wavefront can be used to compute the curvature of the reflecting interface. The method is initially derived for common-source-point data and then extended to common-midpoint data. In a three-dimensional medium the wavefront curvature matrix is computed by solving a 2 × 2 symmetric matrix Riccati equation. In a two-dimensional medium and in a medium with constant velocity gradient, the wavefront curvature matrix is computed by solving a scalar Riccati equation and two linear equations. The migration procedures are also simplified. When the velocity function is unknown, the migration procedures cannot be used. An inverse modeling algorithm which simultaneously performs the migration and estimates the velocity function must then be applied. Two different inversion schemes are discussed briefly.     

Commentary on “Earthquake Depths in the Crimea–Black Sea Region” by V.E. Kulchitsky,B.G. Pustovitenko and V.A. Svidlova

It is well known that the results of determining earthquake parameters depend to a large extent on data processing algorithms and velocity models of the seismic wave propagation medium used in solving hypocenter problems. In 1992, V.Yu. Burmin developed a hypocentric algorithm that minimizes the functional of distances between the points corresponding to the theoretical and observed travel times of seismic waves from an earthquake source to recording stations. The determination of the coordinates of earthquake hypocenters in this case is much more stable than for the commonly used minimization of the functional of discrepancies in the seismic wave arrival times at a station. Using this algorithm and the refined velocity model of the medium, V.Yu. Burmin and L.A. Shumlyanskaya reinterpreted the earthquake parameters for the Crimea–Black Sea region. The most important result of this reinterpretation was the conclusion about the occurrence of deep earthquakes with a source depth of more than 60 km in the region. This result contradicts the conventional beliefs about the seismicity of the region and therefore aroused strong criticism from experts directly involved in compiling the existing catalogs of regional earthquakes. These comments and criticisms are presented by V.E. Kulchitsky with coauthors in a work published in this issue of the journal. In the present paper, we analyze the comments in detail and respond. In particular, we show that the previously used methods of seismic data processing made it highly unlikely by default that deep earthquakes would appear in the results. As an example, we refer to the use of travel-time curves for depths down to 35 km. It is clear that deep earthquakes could not have been found with this approach.     

关于网外地震的定位问题--利用福建区域数字地震台网资料测定台湾地震震源位置

针对地壳内部的网外地震震源的定位问题，由震相组合成不同的方程组，利用均值寻查定位法迭代求解。在选取地球物理参数模型时，考虑了台湾地震部门定位时所用的模型和华南地区地震走时表等。利用福建数字地震台网记录到的绕射波Pn、Sn的到时资料，通过定位程序测定了134个位于台湾区域的壳内地震的震源位置。测定结果与台湾地震部门在其网址上公布的结果进行比对．平均震中位詈相差9km．震源深度相差8km．     

地震定位研究及应用综述

地震定位是地震学中最经典、最基本的问题之一，提高定位精度也一直是地震学应用研究的重要课题之一。本文就目前各种地震定位方法进行了大体分类并概述了基本原理及其应用，介绍了目前在国内各实时地震观测系统中应用的台网定位方法，同时将部分台网方法应用于黑龙江省绥棱爆破的定位。     

Crustal structure of the Caucasus from the Stepnoe-Bakuriani and Volgograd-Nakhichevan DSS profiles (reinterpretation of the primary data)

This work presents the results of reinterpretation of the deep seismic sounding (DSS) data for the Stepnoe-Bakuriani profile and the southern part of the Volgograd-Nakhichevan profile, carried out using new processing methods. Both the profiles cut the trend of the Greater Caucasus. They were acquired in the 1960s by multichannel continuous profiling, which provided high-quality records; however, only the travel-time curves for the main wave types have survived till now. The waves recorded on these profiles have a rather complex origin and their processing by the methods existing at that time was a challenge. At present, the modern computer technologies allowed us to invert the preserved travel-time curves for the velocity models of the Earth??s crust and the very tops of the mantle down to 80 km. It is shown that the crustal thickness increases under the Greater Caucasus up to 50?C60 km and this increase is not gradual, as implied in the previous reconstructions, but occurs through a system of deep dislocations. Traced by the oblique reflections and sharply contrasting seismic velocities, these dislocations extend into the crust. An extended, north-dipping boundary is revealed at a depth of 50?C80 km. The velocity model of the Greater Caucasian crust exhibits slightly decreased velocities compared to the surrounding platform regions. At the same time, the velocities sharply increase in the middle and even upper parts of the crust in the Kura depression.     

华北地区上地幔速度结构

本文选用了北京台网记录的来自不同方位并通过华北地区的140个地震,使用平滑法和dT/d△法,得到P波慢度曲线;采取走时和慢度结合使用的方法,得到走时曲线;从Gerver-Markushevich公式出发,导出一个易于计算的简捷形式,并由此反演初始速度模型;通过模型理论走时曲线的计算,选出与实测数据符合得较好的模型。平均模型中,低速层出现的深度较浅,约60公里;在360-420公里和660-680公里深处,存在两个速度急剧增加的梯度层;500公里处可能还有一个较小的梯度层。华北的东南和西南方向上,相当于660公里梯度层的深度的差异还较为显著。     

Results of deep seismic soundings along international profile VII in Czechoslovakia and poland

Summary The paper deals with the results of DSS measurements along international profile VII, carried out by Czechoslovak and Polish geophysicists in 1970 – 71. The profile situation is shown in Fig. 1. By 1971 part of the profile in the region of the Bohemian Massif between points 1 and 3 and in Poland between points 5 and 7 had been surveyed (Fig. 2). The seismograms were used to construct the travel-time curves of the fundamental types of waves P^K, P^M, P_n (Fig. 4). The mean velocities were computed from the travel-time curves of the reflected waves (P^M and P^K) and the refracted waves (P_g). Isolines of the mean velocities could be constructed for the region of the Bohemian Massif (Fig. 6). The velocity data found were used for the depth interpretation of the travel-time curves of the principal types of waves and to construct a seismic section (Fig. 8). In the region of the Pre-Sudeten block the thickness of the crust was found to be 34–37 km, and in the Sudeten it increased to 40 km. Towards the south the thickness of the crust gradually reduces to 30 km in the system of the Luice faults. In the Bohemian Cretaceous the thickness of the crust is about 30 km. Further to the south, in the region of the Moldanubicum, the thickness of the Earth's crust increases rapidly, and at the southern border of the Central Bohemian pluton it reaches values of about 42 km.     

2012年6月24日宁蒗-盐源 MS5.7地震的精确定位

采用双差定位方法对2012年6月24日宁蒗-盐源M_S5.7地震主震及其余震序列进行了重新定位. 首先采用时间域互相关技术对波形数据进行处理, 获得高精度的事件对走时差; 然后分别使用目录数据、 波形互相关数据和目录数据+波形互相关数据等3种不同的数据组合, 对宁蒗-盐源M_S5.7地震主震及其余震序列进行重新定位. 3种数据组合的主震定位结果、 参与重定位的地震事件数、 最终被定位的地震事件数及重定位结果标准差等对比结果表明, 同时使用目录数据和波形互相关数据所获得的重定位结果最佳. 重定位结果显示, 宁蒗-盐源M_S5.7地震主震的震中位置为(27.790°N、 100.707°E), 震源深度为10.4 km, 发震时刻为北京时间2012年6月24日15时59分32.74秒. 本文结果表明, 震源区附近NW走向的永宁断裂为宁蒗-盐源M_S5.7地震的发震构造.      

Analysis of the IMS Location Accuracy in Northern Eurasia and North America Using Regional and Global Pn Travel-time Tables

—?Joint Research Program of Seismic Calibration of the International Monitoring System (IMS) in Northern Eurasia and North America has been signed by the Nuclear Treaty Programs Office (NTPO), Department of Defense USA, and the Special Monitoring Service (SMS) of the Ministry of Defense, Russian Federation (RF). Under the Program historical data from nuclear and large chemical explosions of known location and shot time, together with appropriate geological and geophysical data, has been used to derive regional Pn/P travel-time tables for seismic event location within the lower 48 States of the USA and the European part of the RF. These travel-time tables are up to 5?seconds faster in shields than the IASPEI91 tables, and up to 5?seconds slower in the Western USA. Relocation experiments using the regional Pn travel-time curves and surrogate networks for the IMS network generally improved locations for regional seismic events. The distance between true and estimated location (mislocation) was decreased from an average of 18.8?km for the IASPEI91 tables to 10.1?km for the regional Pn travel-time tables. However, the regional travel-time table approach has limitations caused by travel-time variations inside major tectonic provinces and paths crossing several tectonic provinces with substantially different crustal and upper mantle velocity structure.¶The RF members of the Calibration Working Group (WG): Colonel Vyacheslav Gordon (chairman); Dr. Prof. Marat Mamsurov, and Dr. Nikolai Vasiliev. The US members of the WG: Dr. Anton Dainty (chairman), Dr. Douglas Baumgardt, Mr. John Murphy, Dr. Robert North, and Dr. Vladislav Ryaboy.     

Earthquake-Generating Structure During the Period of the 1992 Sichuan Lima Earthquake Sequence

In this paper,a test or alternative scheme for studying large earthquake sequences through the study of small earthquake sequences is suggested,and a small earthquake sequence,the Lima earthquake sequence for which analogue records have been turned into digital data,is used here.In order to provide the deep construction background and the spatial distribution of structure for generating earthquakes,the P-wave and S-wave layered velocity models in this area are obtained by using mine explosion and earthquake observed records; then,the hypocenter locations and focal depths of the Lima earthquake sequence are determined adopting the velocity models given above and using a location method with numerical properties for a microseismic monitoring network(Zhao et al.,1994)and a new method for determining focal depth from data of a local seismographic network(Zhao,1992); finally,based on this,the variation of quality factor Q of the crustal medium during the period of the sequence is estimated.The obtained resul     

Fresnel zone inversion for lateral heterogeneities in the earth

We propose a different kind of seismic inversion from travel-time or waveform inversion for lateral heterogeneities in the earth: Fresnel zone inversion. Amplitude and phase delay of data in several frequency ranges are inverted for model space around ray paths with a width corresponding to the considered frequency so that primary effect of finiteness of wavelength be included. For vertically heterogeneous media, Fréchet derivatives for inversion are obtained very efficiently using the paraxial ray approximation, with nearly similar amounts of computation compared to travel-time inversion. As an example, Fréchet derivatives are computed for a teleseismic observation system for a three-dimensional structure in the lithosphere beneath an array of seismic stations. Even if the used frequency is around 2 Hz, the width of Fréchet derivatives cannot be neglected, particularly near the bottom of the lithosphere. Sensitivity of model parameters to observations is, moreover, different in our approach from conventional travel-time inversion: it is zero along ray paths but large slightly away from them. Some model calculations show that travel-time inversion, particularly with models divided into very fine meshes or blocks, might give misleading results. An example of inversion for a simple Camembert model, in the event that travel-time inversion gives no reliable results, shows how this technique works with much smaller data sets and computation than waveform inversions.     

Quantitative analysis of travel-time curves of seismic events in the Ostrava-Karviná coal mining district

Summary The paper is intended as a contribution to the quantitative analysis of travel-time curves of seismic events recorded in the Ostrava-Karviná District (OKD). The input data represent a set of 2621 seismic events, recorded by the local seismological network of 26 mine stations DSLA and a regional diagnostic polygon consisting of five surface Lennartz stations. All the events were processed automatically in the Operational Seismological Centre of the Czechoslovak Army Mine in Karviná and stored in the seismological data base. The results are presented in the form of graphs of arrival times versus distance for the whole OKD, for two mines and one tectonic block.Travel-time curves of direct P and S waves, as well as of reflected and refracted waves are given. The direct P and S waves propagate well practically throughout the whole region studied, but their apparent velocities of propagation are affected by the properties of the rock medium.As a result of the complicated geological conditions, the recorded wave image is quite complicated. Methods of mathematical modelling, using kinematic and dynamic parameters of seismic waves, will have to be applied to identify the separate wave groups uniquely.     

一种修定震源参数的方法

本文描述了空间域中的一种定位方法。此方法的主要特点是:1.震中位置与震源深度、发震时刻分离求解;2.震中位置的确定是采用Romney提出的曾融生改进的台偶到时差来建立最小平方条件方程,几乎不受结构的影响;3.震源深度和发震时刻的确定是采用慢度定深度、到时曲线时间截距定发震时刻的方法来实现的。本文用四川境内的5个工业爆破资料,作了方法试验并修定了近年发生的50个天然地震。修定后的震中位置更加靠近断裂带,走时残差减小了,震源深度呈现出沿龙门山断裂由北东向南西逐渐加深的规律性。     

The earthquake of February 3, 2015, near Sumy,Ukraine: Source parameters and focal mechanism

The parameters of the earthquake that took place February 3, 2015, near the city of Sumy, Ukraine, were calculated from an analysis of records obtained by both Russian and Ukrainian seismic stations (Poltava, Skvira, Nikolaev, Dneropetrovsk, and Desna). The calculated hypocenter depth of 54 km was verified by several approaches: isolation of deep PP, SP phases from the records of remote stations and solution of the kinematic problem for the Poltava station. The focal mechanism as shear with a complex fault component was determined by the first arrivals of P-waves. The data on the azimuthal travel-time curve confirm the focal mechanism. We have calculated the earthquake parameters; they are as follows: length gap L₁ = 8.08 km, L₂ = 6.68 km, a destruction rate of C = 2 km/s. We have obtained the dynamic parameters of the event. The calculated fault length (L = 5.46 km) within the accuracy limits of the method coincides with the early result obtained by the azimuthal travel-time curve. On the basis of these results, we suggest that elastic energy release and formation of the dislocations in the earthquake source occurred on a smooth, prefractured fault (σ_r > 0). Association of the hypocenter with the tectonic node of the northern marginal fault of the Dnieper–Donets graben and northern branch of Kryvyi Rih–Kremenchuk suture confirm this. Here, we observe a considerable Moho depth, structural alteration, and high gradients of the temperature and magnetic and electric rock properties in the lower Earth’s crust and upper mantle. These circumstances are favorable for the earthquake occurring here.     

On the Problem of Deep Earthquakes in Crimea–Black Sea Region

It is a common opinion that only crustal earthquakes can occur in the Crimea–Black Sea region. Since the existence of deep earthquakes in the Crimea–Black Sea region is extremely important for the construction of a geodynamic model for this region, an attempt is made to verify the validity of this widespread view. To do this, the coordinates of all earthquakes recorded by the stations of the Crimean seismological network are reinterpreted with an algorithm developed by one of the authors. The data published in the seismological catalogs and bulletins of the Crimea–Black Sea region for 1970–2012 are used for the analysis. To refine the coordinates of hypocenters of earthquakes in the Crimea–Black Sea region, in addition to the data from stations of the Crimean seismological network, information from seismic stations located around the Black Sea coast are used. In total, the data from 61 seismic stations were used to determine the hypocenter coordinates. The used earthquake catalogs for 1970–2012 contain information on ~2140 events with magnitudes from–1.5 to 5.5. The bulletins provide information on the arrival times of P- and S-waves at seismic stations for 1988 events recorded by three or more stations. The principal innovation of this study is the use of the original author’s hypocenter determination algorithm, which minimizes the functional of distances between the points (X, Y, H) and (x, y, h) corresponding to the theoretical and observed seismic wave travel times from the earthquake source to the recording stations. The determination of the coordinates of earthquake hypocenters is much more stable in this case than the usual minimization of the residual functional for the arrival time of an earthquake wave at a station (the difference between the theoretical and observed values). Since determination of the hypocenter coordinates can be influenced by the chosen velocity column beneath each station, special attention is focused on collecting information on velocity profiles. To evaluate the influence of the upper mantle on the results of calculating the velocity model, two different low-velocity and high-velocity models are used; the results are compared with each other. Both velocity models are set to a depth of 640 km, which is fundamentally important in determining hypocenters for deep earthquakes. Studies of the Crimea–Black Sea region have revealed more than 70 earthquakes with a source depth of more than 60 km. The adequacy of the obtained depth values is confirmed by the results of comparing the initial experimental data from the bulletins with the theoretical travel-time curves for earthquake sources with depths of 50 and 200 km. The sources of deep earthquakes found in the Crimea–Black Sea region significantly change our understanding of the structure and geotectonics of this region.     

Improving Regional Seismic Event Location in China

—?In an effort to improve our ability to locate seismic events in China using only regional data, we have developed empirical propagation path corrections and applied such corrections using traditional location routines. Thus far, we have concentrated on corrections to observed P arrival times for crustal events using travel-time observations available from the USGS Earthquake Data Reports, the International Seismic Centre Bulletin, the preliminary International Data Center Reviewed Event Bulletin, and our own travel-time picks from regional data. Location ground truth for events used in this study ranges from 25?km for well-located teleseimic events, down to 2?km for nuclear explosions located using satellite imagery. We also use eight events for which depth is constrained using several waveform methods. We relocate events using the EvLoc algorithm from a region encompassing much of China (latitude 20°–55°N; longitude 65°–115°E). We observe that travel-time residuals exhibit a distance-dependent bias using IASPEI91 as our base model. To remedy this bias, we have developed a new 1-D model for China, which removes a significant portion of the distance bias. For individual stations having sufficient P-wave residual data, we produce a map of the regional travel-time residuals from all well-located teleseismic events. Residuals are used only if they are smaller than 10?s in absolute value and if the seismic event is located with accuracy better than 25?km. From the residual data, correction surfaces are constructed using modified Bayesian kriging. Modified Bayesian kriging offers us the advantage of providing well-behaved interpolants and their errors, but requires that we have adequate error estimates associated with the travel-time residuals from which they are constructed. For our P-wave residual error estimate, we use the sum of measurement and modeling errors, where measurement error is based on signal-to-noise ratios when available, and on the published catalog estimate otherwise. Our modeling error originates from the variance of travel-time residuals for our 1-D China model. We calculate propagation path correction surfaces for 74 stations in and around China, including six stations from the International Monitoring System. The statistical significance of each correction surface is evaluated using a cross-validation technique. We show relocation results for nuclear tests from the Balapan and Lop Nor test sites, and for earthquakes located using interferometric synthetic aperture radar. These examples show that the use of propagation path correction surfaces in regional relocations eliminates distance bias in the residual curves and significantly improves the accuracy and precision of seismic event locations.     

The velocity structure of the Carpathian zone from the ambient noise surface wave tomography

The dispersion curves of the Rayleigh wave group velocities are constructed along 60 interstation seismic paths in Central Europe based on the cross-correlation function of seismic noise. Together with the previous data (Yanovskaya and Lyskova, 2013), this information was used for reconstructing the three-dimensional distribution of S-wave velocities in the upper mantle of the Carpathian region. In the present work, the previous results are refined by expanding the data set by the additional seismic paths that intersect the Carpathian region and by modifying the procedure for constructing the locally averaged dispersion curves so as to obtain a more compact resolution. The results of the study suggest the complex, multidirectional character of the plate motion in the region.     

1992年四川力马地震序列期间激发的孕震结构

对小震序列的深入研究或许可以作为对大震序列震后趋势研究的替代，本文将地震模拟量已转化为数字记录的马力小震序列来尝试这种替代。。     

Subcrustal structure of the black sea basin from seismological data

The P-wave travel time data from the earthquakes offshore and onshore around the Black Sea are used for the tomographic reconstruction of the three-dimensional (3D) velocity distribution in the lithosphere of the region. The preliminary refinement of the foci parameters (the coordinates and origin time) has reduced the random errors in the travel-time data. The earthquake data were supplemented by the previous deep seismic sounding (DSS) data on the profiles in Crimea and offshore off the Black Sea. The dataset included more than 4000 travel times overall. In order to eliminate the crustal effect, the travel times were reduced to a surface at a depth of 35 km corresponding to the mean Moho depth in the region. The improved crustal model was used for removing the contribution of the crust from the initial data. The new tomography method, which was recently developed by one of the authors and which relies on the assumption of smoothness of the lateral velocity variations, was applied for reconstructing the velocity structure of the upper mantle beneath the Black Sea up to a depth of 95 km. The lateral velocity variation maps at different depths and the vertical velocity distributions along the meridional and sublatitudinal cross sections across the Black Sea were constructed. High velocities were revealed in the subcrustal lithosphere, and the structural difference below two subbasins—the West Black Sea (WBS) and the East Black Sea (EBS) ones—was established. It shows that the high-velocity body below the WBS is located deeper than below the EBS and is distinguished by higher velocities. Based on these results, it is concluded that the lithosphere beneath the Black Sea has a continental origin.