Full‐wavefield migration: using surface and internal multiples in imaging

Multiple scattering is usually ignored in migration algorithms, although it is a genuine part of the physical reflection response. When properly included, multiples can add to the illumination of the subsurface, although their crosstalk effects are removed. Therefore, we introduce full‐wavefield migration. It includes all multiples and transmission effects in deriving an image via an inversion approach. Since it tries to minimize the misfit between modeled and observed data, it may be considered a full waveform inversion process. However, full‐wavefield migration involves a forward modelling process that uses the estimated seismic image (i.e., the reflectivities) to generate the modelled full wavefield response, whereas a smooth migration velocity model can be used to describe the propagation effects. This separation of modelling in terms of scattering and propagation is not easily achievable when finite‐difference or finite‐element modelling is used. By this separation, a more linear inversion problem is obtained. Moreover, during the forward modelling, the wavefields are computed separately in the incident and scattered directions, which allows the implementation of various imaging conditions, such as imaging reflectors from below, and avoids low‐frequency image artefacts, such as typically observed during reverse‐time migration. The full wavefield modelling process also has the flexibility to image directly the total data (i.e., primaries and multiples together) or the primaries and the multiples separately. Based on various numerical data examples for the 2D and 3D cases, the advantages of this methodology are demonstrated.     

Crustal structure of northeastern Tibetan plateau and Ordos block: Waveform interpretation of the Maqen-Jingbian seismic refraction profile

The Maqen-Jingbian wide-angle seismic reflection and refraction experiment was carried out in 1998, which aims at determining detailed structure in the crust and top of the upper mantle and understanding structural relation between the northeastern Tibetan plateau and the Ordos block. The 1-D crustal models inferred by waveform inversion show strong variations in crustal structure, which can be classified into four different types: ① an Ordos platform with the Proterozoic crust and two high-velocity layers in the northeast section, ② a transitional crust between the northeastern Tibetan plateau and the Ordos block across the Haiyuan earthquake zone, ③ the Qilian orogenic zone in the central part, and 4 the Qinling orogenic zone in the southwestern section. The Moho depth increases from ~42 km to ~62 km from the NE part to the SW part of the profile. The crystalline crust consists of the upper crust and lower crust in northeastern Tibetan plateau. There is an obviously low P-wave velocity layer dipping northeastward, which is 12–13 km thick, at the bottom of the upper crust in Qinling orogenic zone and Haiyuan earthquake zone. The lower crust is characterized by alternating high and low P-wave velocity layers. Beneath Ordos block, i.e., the NE part of the profile, the crust shows quite a smooth increase in P-wave velocity down to the Moho at a depth of about 42 km.     

有限频率地震层析成像方法及研究进展

本文从基于射线理论的走时地震层析成像的发展引出有限频率地震层析成像, 说明两者最为本质的差别是有限频率灵敏度核函数。 概述了有限频率层析成像方法和基本原理, 详细阐述了国内外有关有限频率层析成像方法的研究进展和应用成果, 强调有限频率层析成像对全球地幔柱研究的重要意义, 最后指出有限频率层析成像方法在理论上的诸多优点目前还没有在实际应用中得到充分体现, 还存在一些问题需要解决。     

Source mechanisms of the June 2004 Tabuk earthquake sequence, Eastern Red Sea margin, Kingdom of Saudi Arabia

A sequence of earthquakes took place in June 2004 approximately 60 km southeast of Tabuk, Saudi Arabia. The first felt event (M _W = 3.9) occurred on June 9 and caused minor damage in the epicentral area according to the National Earthquake Information Center and the local reports. Another moderate size event occurred on June 22 (M _W = 5.1) and was followed by a few felt aftershocks without any reported damage. This earthquake sequence caused considerable alarm at Tabuk and highlights the fact that damaging earthquakes can occur in this region away from the major plate boundary in the Red Sea. Being the largest well-recorded event in the area for which the digital and broadband records from Saudi Arabia, Egypt, Israel, Jordan, Turkey, Cyprus, and Kuwait are available, it provides an excellent opportunity to study the tectonic process and present day stress field acting on this area. The digital records from these regional networks were used to relocate the largest three events of this sequence. Focal mechanisms were obtained from full waveform inversion and indicate normal faulting mechanisms with two nodal planes oriented NW–SE in parallel to the faults bounding the Tabuk graben and the Red Sea rift axis. These events originated at shallow focal depths of 4–5 km, possibly contributing to the widely felt ground motions. These events offer a unique opportunity to study the active tectonics of the region as well as inform future studies of seismic hazard in northwestern Saudi Arabia, the Gulf of Aqaba, and eastern Egypt.     

使用双线性变换构造实时地震波形仿真滤波器的方法研究

介绍了直接使用双线性变换把模拟系统的传递函数转换为离散系统的传递函数,构造出由宽带数字地震记录,仿真出传统模拟地震仪器记录波形的递归滤波器的方法,从而使在线实时波形仿真成为可能。并通过比较利用实际宽带数字地震记录,分别采用频域仿真和本文算法构造的仿真滤波器仿真得到的波形,验证了本文算法的正确性。     

北京遥测地震台网数字记录与模拟记录在波形对比和震级对比方面的研究

利用“九五”期间北京遥测地震台网部分地震台站模拟信号和数字信号并行的条件,通过开展模拟记录、数字记录与仿真处理记录在波形对比和震级对比方面的研究工作,我们认为可以实现区域地震台网从模拟记录系统到数字记录系统在常规数据处理方面的平稳过渡,建议在近震的数字记录中采用DD-1仿真震级。     

国家测震台网中心强震动数据处理系统的设计与实现

为满足国家测震台网海量强震动观测数据的处理时效性、格式标准化、产品丰富度等需求,开发了兼容多类强震动观测站点且具备数据快速汇集、处理及归档等功能的强震动数据处理系统。该系统提供地震波形人机交互数据预处理界面,对加速度记录进行预处理,进一步分析预处理后的加速度事件波形数据,计算得到地震动各项参数,包括峰值加速度(PGA)、峰值速度(PGV)、峰值位移(PGD)、仪器烈度、持时、傅氏谱、反应谱和三联谱等,可以导出地震元数据、地震记录波形,对各类数据进行归档存储。该系统具有平台统一性、功能集成性、数据完备性等特点,有效提升了日常数据处理和管理能力,能在地震应急、震害评估和科学研究中发挥实效。     

The October 11, 1999 and November 08, 2006 Beni Suef Earthquakes

Two felt moderate-sized earthquakes with local magnitudes of 4.9 on October 11, 1999 and 4.3 on November 08, 2006 occurred southeast of Beni Suef and Cairo cities. Being well recorded by the digital Egyptian National Seismic Network (ENSN) and some regional broadband stations, they provided us with a unique opportunity to study the tectonic process and present-day stress field acting on the northern part of the Eastern Desert of Egypt. In this study, we analyze the main shocks of these earthquakes and present 15 well recorded aftershocks (0.9 ≤ ML ≤ 3.3) which have small errors on both horizontal and vertical axes. The relocation analysis using the double difference algorithm clearly reveals a NW trending fault for the 1999 earthquake. The spatial distribution of its aftershocks indicates a propagation of rupture from the SW towards the NW along a fault length ~5 km dipping nearly ~40°SW. We also determined the focal mechanisms of the two main shocks by two methods (polarities and amplitudes ratios of P, S_V and S_H and regional waveform inversion). Our results indicate a normal faulting mechanism with a slight shear component for the first event, while pure normal faulting for the second one. The spatial distribution of the 1999 aftershocks sequence along with the retrieved focal mechanism confirmed the NW plane as the true fault plane. While for the 2006 event, the few aftershocks do not reveal any fault geometry; its focal mechanism indicated a pure normal fault nearly trending WNW-ESE that corresponds more likely to the extension of the 1999 earthquake fault. The seismicity distribution between the two earthquake sequences reveals a noticeable gap that may be a site of a future event. The NNE-SSW extensional stress indicated by the mechanisms of these events is in agreement with the regional stress field and the rifting of the northern Red Sea in its northern branches (Gulf of Suez and Gulf of Aqaba). The source parameters (seismic moment, moment magnitude, fault radius, stress drop and displacement across the fault) were also estimated and compared based on both the regional waveform inversion and the displacement spectra and interpreted in the context of the tectonic setting. The obtained results imply a reactivation of the pre-exiting NW-SE faults as a result of extensional deformation from the northern Red Sea-Gulf of Suez rifts.     

Synthetic study on SS waveform splitting and complication

Some influential factors on the complication of SS waveform (with epicentral distance within 40°～180°) are ana- lyzed quantitatively by calculating the full-wave synthetic seismogram using propagation matrix method. Our re- sults show that the transmission-conversion and reflection-conversion phases of S wave at the interface of Moho and free surface beneath bounce points are mainly responsible for the complication of SS waveform, the velocity contrast between the two sides of Moho boundary under SS bounce point also has great effects on the amplitudes of all the conversion and the reverberation phases; the properties of the crust at seismic station also play a role in the complication of SS waveform while the crustal thickness beneath bounce point is thinner relatively. At the same time, two sets of real SS waveform data at the two positions in eastern and western China are analyzed, and the splitting time between SS transverse component and radial component is measured by cross-correlation. Our analysis demonstrates that there is a positive correlation between crustal thickness and the splitting time because of the influences of adjacent conversion and reverberation phases, the splitting time in west with thick crust is ob- viously greater than that in the east with thin crust. Moreover, It is promising that one new method of measurement of crustal thickness will be developed by using the observed SS splitting time.     

Region-related features of crustal and upper-mantle velocity structure of the Chinese mainland detected by surface waveform modeling

Surface waveform modeling has played an important role on many continental-scale studies of upper mantle ve-locity structure, but it was seldom applied to the Chinese mainland study. The present study firstly analyzed sur-face waveform fittings for eight wave paths crossing through four different regions of the Chinese mainland (east-ern, central, northern and western China), and then inverted for 1D path-averaged S-velocities for these paths. The inverted crustal and upper-mantle S-velocities showed obvious region-related features, which are well consistent with known geotectonic units and previous research results. These results indicate that surface waveform modeling is a reliable method to get crustal and upper-mantle velocity structure. Furthermore, this method has a prominent advantage in detecting upper-mantle structure compared with fundamental-mode surface-wave dispersion method.