地下核爆炸的主要非爆炸源机制

通常将地下封闭式爆炸近似为球对称爆炸震源来研究,然而利用这一模型来解释由地下核爆炸激发的近场和远场地震波的全部特征是远远不够的,通过理论和利用地下核爆炸产生的地震波数据的研究,学者们总结出了一些重要的辅助震源模式,作为地下核爆炸的二次源机制,如地震波转换震源模式、层裂震源模式、构造应变能释放震源模式等.其中的一些模型成功地解释了诸如地下核爆炸Lg波的激发机制等重要问题.本文介绍了这些重要模式的研究背景及理论依据,并且详细地推导了层裂源的数学描述过程.     

不同模型下爆炸源和层裂源的 Lg波特性研究

地下核爆炸区域及远区地震波信号中,有二次源即层裂源的贡献,其对区域震相Lg波具有调制作用.本文利用理论地震图数值模拟方法,分析爆炸源和层裂源在3种典型介质速度模型下的区域震相特征,得到了爆炸源和层裂源所激发的Lg波与介质速度模型的关系,并结合CLVD源所激发的S波随震中距的演化过程,对Lg波的形成机制进行了阐述.结果表明,CLVD源是激发地下核爆炸低频Lg波的主要因素.      

基于Lg波识别判据的特性分析

Lg波是检测、识别地下核爆炸的重要震相.通常情况下,地下核爆炸的低频段Lg波能量相对较强,基于Lg波的识别判据在高频段表现较好.然而在低频范围内却表现出难以用球对称爆炸源模式解释的现象,如Lg/Pg判据在1 Hz附近失效等现象.本文利用理论地震图方法,分析了伴随层裂过程产生的重要辅助源——CLVD源对Lg波低频成分的调...     

地下核爆炸Lg波的激发机制

地下核爆炸Lg波的激发机制,目前尚未十分清楚.其中主要观点之一是源区附近由补偿线性矢量偶极源（CLVD）激发的Rg波的散射形成的S波是Lg波能量的主要来源.本文利用理论地震图模拟方法,基于东哈萨克斯坦地区速度模型,分析比较了东哈萨克斯坦地壳速度模型下的爆炸源、张裂源及CLVD源对区域震相Lg波的影响.结果表明,从能量的角度来看,CLVD源是激发Lg波的主要因素.模拟计算结果也进一步证实了CLVD源激发的Lg波振幅谱具有低谷点的特征源自于该震源所激发的Rg波;在此基础之上,检验了Patton提出的估算震源埋藏深度的经验公式.结果发现,该公式仅适用于震源埋藏深度较浅的情况（<500m）.这些结果对于进一步理解及更好地利用Lg波具有重要理论指导意义.     

地下爆炸Rg波低谷点激发机理

目前普遍认为Rg_S散射是地下核爆炸短周期Lg波的主要能量来源,其主要依据为由线性偶极补偿源（CLVD）激发的Rg波谱的低谷点随深度变化的特征与区域震相Lg波谱的低谷点随深度变化的特征具有一致性. 本文利用简正振型理论,分析了Rg波谱中低谷点的形成机理,给出了各种速度模型下CLVD源深度与低谷点频率之间的关系,并分析了利用公式hCLVD=V/(16fNull)估计震源深度的局限性. 本文的研究结果进一步支持了CLVD源是激发Lg波的主要因素的观点,并得到重要结论：CLVD源形式下的本征位移函数及其导数的叠加所形成的极小值,即源的空间分布特征,是形成Rg波低谷点的原因,而不能仅用水平向基频本征位移函数过零点来解释.     

地下核爆炸地震波二次源特征

层裂是地下核爆炸震源物理过程中最重要的非球对称源过程.对其机制的研究,多基于面波资料,本文基于远场体波资料,利用经验格林函数方法,消除路径及接收场地效应,利用地下核爆炸远场地震波P波数据对震源时间函数进行了反演,结果显示明显的双脉冲特征,该二次源特征不随方位角及震中距变化,说明它源自震源本身,其双脉冲特征与层裂机制相符.     

核爆、化爆、地震识别研究综述

地震学方法是监测地下核试验最为重要的方法.为准确鉴别地下核试验,不仅包括识别核爆炸与天然地震,还包括对化学爆炸震源类型的识别.论文从核爆炸与天然地震的识别、化学爆炸与地震的识别、核爆炸与化学爆炸的识别三个方面对国内外的研究状况分别进行了综述,较为详尽地阐述了从统计分析和信号处理的角度进行研究的进展状况,并对目前的研究进展进行了分析总结,最后对未来的核爆监测识别研究的发展趋势进行了展望,提出了一些建议.     

基于理论地震图方法的波导结构对 地下核爆炸Lg波传播特性研究

基于频率-波数域算法的理论地震波形图方法, 可以数值模拟频率达到10 Hz、 震中距达1000 km的区域理论地震波形图. 该算法适用于计算大量分层地壳结构中激发的导波Lg波. 本文在前苏联东哈萨克斯坦地下核试验场至我国乌鲁木齐台站间的地球介质速度模型中, 引入速度梯度结构、 速度扰动分布的薄叠加层结构、 降低Q值结构以及速度扰动与Q值变化的综合结构来模拟实际地壳波导结构, 较好地模拟出东哈萨克斯坦地下核爆炸地震在乌鲁木齐台站记录的宽频带地震波形图, 模拟出完整的Lg波序列, 该序列符合Lg波能量分布特征, 且能够解释Lg波波尾的特征. 结果表明, Lg波的形状和峰值结构均依赖于地壳的不同波导结构.      

震源性质的时频分析与事件识别

应用时频分析技术研究了近年来发生在朝鲜半岛的核爆炸、化学爆炸与天然地震事件的识别问题,计算了上述三者的瞬时频率、瞬时振幅谱与群延迟等多项参数,提出地震与爆炸的识别指标.初步研究结果表明,对较大事件,上述三者来自震源的激发频率信息有明显差异性,核爆炸源激发的频率明显高于地震,化学爆炸最低.特别是,这些信息指标对于识别小当...     

层裂对区域震相Lg波的影响

地下核爆炸Lg波的激发机制, 目前尚未十分清楚. 普遍接受的观点是: 伴随层裂源的补偿线性矢量偶极源激发的Rg波的散射形成的S波, 是Lg波的最主要贡献因素. 本文利用频谱比对方法, 结合合成理论地震图方法, 基于东哈萨克斯坦地区速度模型, 分析了东哈萨克斯坦地壳速度模型下的层裂对区域震相Lg波的影响, 并进一步分析了该地区的实际观测资料. 结果表明, 层裂时间函数对Lg波具有扇贝形调制作用; Lg波振幅谱的低谷点源自于补偿线性矢量偶极（CLVD）源激发的Rg波. 这一结果也表明, 东哈萨克斯坦速度模型下的Lg波受层裂的调制作用, 与内华达试验区的Lg波是一样的, 进一步支持了Patton和Taylor的观点, 对于进一步理解及更好地利用Lg波具有重要理论指导意义.      

Source Parameters of the 8 October, 2005 M<Subscript>w</Subscript>7.6 Kashmir Earthquake

During the last six years, the National Geophysical Research Institute, Hyderabad has established a semi-permanent seismological network of 5 broadband seismographs and 10 accelerographs in the Kachchh seismic zone, Gujarat, with the prime objective to monitor the continued aftershock activity of the 2001 M_w7.7 Bhuj mainshock. The reliable and accurate broadband data for the M_w 7.6 (8 Oct., 2005) Kashmir earthquake and its aftershocks from this network, as well as from the Hyderabad Geoscope station, enabled us to estimate the group velocity dispersion characteristics and the one-dimensional regional shear-velocity structure of peninsular India. Firstly, we measure Rayleigh- and Love-wave group velocity dispersion curves in the range of 8 to 35 sec and invert these curves to estimate the crustal and upper mantle structure below the western part of peninsular India. Our best model suggests a two-layered crust: The upper crust is 13.8-km thick with a shear velocity (Vs) of 3.2 km/s; the corresponding values for the lower crust are 24.9 km and 3.7 km/sec. The shear velocity for the upper mantle is found to be 4.65 km/sec. Based on this structure, we perform a moment tensor (MT) inversion of the bandpass (0.05–0.02 Hz) filtered seismograms of the Kashmir earthquake. The best fit is obtained for a source located at a depth of 30 km, with a seismic moment, Mo, of 1.6 × 10²⁷ dyne-cm, and a focal mechanism with strike 19.5°, dip 42°, and rake 167°. The long-period magnitude (M_A ~ M_w) of this earthquake is estimated to be 7.31. An analysis of well-developed sPn and sSn regional crustal phases from the bandpassed (0.02–0.25 Hz) seismograms of this earthquake at four stations in Kachchh suggests a focal depth of 30.8 km.     

Crustal structure of Tushka Region, Abu-Simbel, Egypt, inferred from spectral ratios of P waves of local earthquakes

The crustal structure of North Abu-Simbel area was studied using spectral ratios of short-period P waves. Three-component short period seismograms from the Masmas seismic station of the Egyptian National Seismic Network Stations were used. The Thomson-Haskell matrix formulation was applied for linearly elastic, homogeneous crustal layers. The obtained model suggests that the crust under the study region consists of a thin (0.8 km) superficial top layer with a P-wave velocity of 3.8±0.7 km/s and three distinct layers with a mean P-wave velocity of 6.6 km/s, overlaying the upper mantle with a P-wave velocity of 8.3 km/s (fixed). The results were obtained for 14 different earthquakes. The P-wave velocities of the three layers are: 5.8±0.6 km/s, 6.5±0.4 km/s and 7.2±0.3 km/s. The total depth to the Moho interface is 32±2 km. The crustal velocity model estimated using observations is relatively simple, being characterized by smooth velocity variations through the middle and lower crust and normal crustal thickness. The resultant crustal model is consistent with the model obtained from previous deep seismic soundings along the northern part of Aswan lake zone.     

Source Parameters of the Deadly M<Subscript>w</Subscript> 7.6 Kashmir Earthquake of 8 October, 2005

During the last six years, National Geophysical Research Institute, Hyderabad has established a semi-permanent seismological network of 5–8 broadband seismographs and 10–20 accelerographs in the Kachchh seismic zone, Gujarat with a prime objective to monitor the continued aftershock activity of the 2001 M_w 7.7 Bhuj mainshock. The reliable and accurate broadband data for the 8 October M_w 7.6 2005 Kashmir earthquake and its aftershocks from this network as well as Hyderabad Geoscope station enabled us to estimate the group velocity dispersion characteristics and one-dimensional regional shear velocity structure of the Peninsular India. Firstly, we measure Rayleigh-and Love-wave group velocity dispersion curves in the period range of 8 to 35 sec and invert these curves to estimate the crustal and upper mantle structure below the western part of Peninsular India. Our best model suggests a two-layered crust: The upper crust is 13.8 km thick with a shear velocity (Vs) of 3.2 km/s; the corresponding values for the lower crust are 24.9 km and 3.7 km/sec. The shear velocity for the upper mantle is found to be 4.65 km/sec. Based on this structure, we perform a moment tensor (MT) inversion of the bandpass (0.05–0.02 Hz) filtered seismograms of the Kashmir earthquake. The best fit is obtained for a source located at a depth of 30 km, with a seismic moment, Mo, of 1.6 × 10²⁷ dyne-cm, and a focal mechanism with strike 19.5°, dip 42°, and rake 167°. The long-period magnitude (M_A ~ M_w) of this earthquake is estimated to be 7.31. An analysis of well-developed sPn and sSn regional crustal phases from the bandpassed (0.02–0.25 Hz) seismograms of this earthquake at four stations in Kachchh suggests a focal depth of 30.8 km.     

An approach of refraction seismology for processing and modeling of local earthquake seismogram sections of virtual sources at multiple depths in seismogenic regions—application to Koyna-Warna region,India, for upper crustal P and S velocity structure

We present an approach based on controlled source seismology (CSS) methods, especially developed for processing and modeling of the local earthquake seismograms. Record sections of the local earthquake seismograms generated for multiple source depths illuminate the upper crustal velocity structure in the region. Extensive travel times and synthetic seismograms modeling of the observed record sections reveal the P and S velocity structure in the region. The strength of this approach essentially lies with the possibility of validating the upper crustal velocity models inferred in various subregions of the seismogenic region. A redundant and significantly large number of virtual source local earthquake seismogram sections, gathered for multiple source depths and varying source mechanisms in each of the subregions, validate the same set of P and S velocity models in that region. Further, those models are found to generate the synthetic seismograms consistent with the observed sections. The proposed approach effectively utilizes a reliable dataset from a great volume of well-located local earthquake recordings of a state-of-the-art digital seismograph network. Such a dataset of local earthquake seismograms in the Koyna-Warna active earthquake zone is used here to demonstrate this approach and obtained subregion-specific models of upper crustal P and S velocity structure in the epicentral region. The results indicate that the technique presented here is efficient for processing and modeling the local earthquake seismograms and deriving upper crustal velocity models in the seismogenic regions.     

Response spectral attenuation relationships for Singapore and the Malay Peninsula due to distant Sumatran‐fault earthquakes

As part of the effort to assess the seismic hazards of Singapore and the Malay Peninsula, representative ground motion prediction models have to be established. Seven existing attenuation relationships developed for shallow crustal earthquakes in stable continent and active tectonic regions are examined, and they are found to consistently over‐predict the ground motions of Sumatran‐fault earthquakes recently recorded in Singapore. This may be attributed to the differences in the regional crustal structures and distance ranges considered. Since the number of recorded ground motions in the region is very limited, a new set of attenuation relationships is derived based on synthetic seismograms. The uncertainties in rupture parameters, such as stress drop, focal depth, dip and rake angles, are defined according to the regional geological and tectonic settings as well as the ruptures of previous earthquakes. Ground motions are simulated for earthquakes with Mw ranging from 4.0 to 8.0, within a distance range from 174 to 1379km. Besides magnitude and distance, source‐to‐station azimuth is found to influence the amplitudes of the ground motions simulated. Thus, the azimuth is taken as an independent variable in the derived ground motion attenuation relationships. The Sumatran‐fault segments that have the potential to generate a specified level of response spectral accelerations in Singapore and Kuala Lumpur are identified based on the newly derived ground motion models. Copyright © 2003 John Wiley & Sons, Ltd.     

Spectral element analysis on the characteristics of seismic wave propagation triggered by Wenchuan <Emphasis Type="Italic">M</Emphasis><Subscript>s</Subscript>8.0 earthquake

The 2008 Wenchuan earthquake occurred in an active earthquake zone, i.e., Longmenshan tectonic zone. Seismic waves triggered by this earthquake can be used to explore the characteristics of the fault rupture process and the hierarchical structure of the Earth’s interior. We employ spectral element method incorporated with large-scale parallel computing technology, to investigate the characteristics of seismic wave propagation excited by Wenchuan earthquake. We calculate synthetic seismograms with one-point source model and three-point source model respectively. The AK135 model is employed as a prototype of our numerical global Earth model. The Earth’s ellipticity, Earth’s medium attenuation, and topography data are taken into consideration. These wave propagation processes are simulated by solving three-dimensional elastic wave governing equations. Three-dimensional visualization of our numerical results displays the profile of the seismic wave propagation. The three-point source, which is proposed from the latest investigations through field observation and reverse estimation, can better demonstrate the spatial and temporal characteristics of the source rupture process than the one-point source. We take comparison of synthetic seismograms with observational data recorded at 16 observatory stations. Primary results show that the synthetic seismograms calculated from three-point source agree well with the observations. This can further reveal that the source rupture process of Wenchuan earthquake is a multi-rupture process, which is composed by at least three or more stages of rupture processes. Supported by National Basic Research Program of China (Grant No. 2004CB418406), National Natural Science Foundation of China (Grant Nos. 40774049 and 40474038), and Computer Network Information Center, Chinese Academy of Sciences (Grant No. INF105-SCE-02-12)     

Structure of the Lower Crust Revealed by One- and Two-dimensional Modeling of Wide-angle Reflections — West Bengal Basin, India

--The crustal structure in the West Bengal basin, India has been investigated by means of wide-angle reflection data recorded along (i) Bishnupur-Palashi-Kandi, 227-km long profile in the north-south direction and (ii) Taki-Arambagh, 120-km long profile in the east-west direction. The data were acquired using multichannel digital seismic instruments with close station spacing. The crustal model, initially derived by 1-D forward modeling of the wide-angle reflection data, has been iteratively refined by 2-D ray tracing and modeling of travel-time observations and the corresponding synthetic seismograms computation. The structural contour map of the Moho prepared from the present data set, indicates the crustal thickness of about 37 km in the western margin of the basin, thinning to about 28 km in the east with an upwarp in the Moho boundary. The upwarp in the Moho and the inferred structural features may be indicative of crustal rifting. The well-known gravity anomaly in the West Bengal basin, 'Calcutta gravity high,' appears to have resulted from the Moho upwarp in combination with the huge thickness of sediments deposited east of the steep flexure of the crystalline basement representing the 'Hinze zone.'     

SEISMOGRAM SYNTHESIS AND RECOMPRESSION OF DISPERSIVE IN-SEAM SEISMIC MULTIMODE DATA USING A NORMAL-MODE SUPERPOSITION APPROACH

In spite of a geometrical rotation into radial and transverse parts, two- or three-component in-seam seismic data used for underground fault detection often suffer from the problem of overmoding ‘noise’. Special recompression filters are required to remove this multimode dispersion so that conventional reflection seismic data processing methods, e.g. CMP stacking techniques, can be applied afterwards. A normal-mode superposition approach is used to design such multimode recompression filters. Based on the determination of the Green's function in the far-field, the normal-mode superposition approach is usually used for the computation of synthetic single- and multi-mode (transmission) seismograms for vertically layered media. From the filter theory's point of view these Green's functions can be considered as dispersion filters which are convolved with a source wavelet to produce the synthetic seismograms. Thus, the design of multimode recompression filters can be reduced to a determination of the inverse of the Green's function. Two methods are introduced to derive these inverse filters. The first operates in the frequency domain and is based on the amplitude and phase spectrum of the Green's function. The second starts with the Green's function in the time domain and calculates two-sided recursive filters. To test the performance of the normal-mode superposition approach for in-seam seismic problems, it is first compared and applied to synthetic finite-difference seismograms of the Love-type which include a complete solution of the wave equation. It becomes obvious that in the case of one and two superposing normal modes, the synthetic Love seam-wave seismograms based on the normal-mode superposition approach agree exactly with the finite-difference data if the travel distance exceeds two dominant wavelengths. Similarly, the application of the one- and two-mode recompression filters to the finite-difference data results in an almost perfect reconstruction of the source wavelet already two dominant wavelengths away from the source. Subsequently, based on the dispersion analysis of an in-seam seismic transmission survey, the normal-mode superposition approach is used both to compute one- and multi-mode synthetic seismograms and to apply one- and multimode recompression filters to the field data. The comparison of the one- and two-mode synthetic seismograms with the in-seam seismic transmission data reveals that arrival times, duration and shape of the wavegroups and their relative excitation strengths could well be modelled by the normal-mode superposition approach. The one-mode recompressions of the transmission seismograms result in non-dispersive wavelets whose temporal resolution and signal-to-noise ratio could clearly be improved. The simultaneous two-mode recompressions of the underground transmission data show that, probably due to band-limitation, the dispersion characteristics of the single modes could not be evaluated sufficiently accurately from the field data in the high-frequency range. Additional techniques which overcome the problem of band-limitation by modelling all of the enclosed single-mode dispersion characteristics up to the Nyquist frequency will be mandatory for future multimode applications.