gPhone重力仪的面波频段响应实测研究

以中国大陆构造环境监测网络昆明台和恩施台gPhone相对重力仪连续重力潮汐观测数据为基础,研究了gPhone重力仪在1 mHz以上频段的高频响应。从瑞利面波角度获得gPhone重力仪的高频响应,并且通过与同址观测的STS-1地震仪LHZ分量数据进行对比,验证了gPhone重力仪高频观测结果的可靠性。对比从gPhone重力仪和STS-1地震仪观测数据中提取到的面波波形和群速度频散曲线,发现昆明台两类仪器观测到面波信号的振幅和相位都较为一致,而恩施台仅振幅较为一致,相位上存在较明显差异,gPhone重力仪记录的面波信号在各频段存在不同的时间延迟。用两类仪器观测数据获得了大地震激发的自由振荡,结果表明两类仪器观测到的基频球型模态自由振荡的频率和振幅都吻合较好,进一步验证了gPhone重力仪对高频频段信号振幅响应的可靠性。以上研究结果表明:利用gPhone重力仪能够准确地观测到大地震激发的面波和自由振荡等高频信号的振幅,但在记录信号的相位信息时,有些仪器会有相位偏移产生,如果研究中需要考虑信号的相位,则必须获得仪器相位偏移量,再进行仪器相位校正。     

Physically compliant,conditionally simulated spatially variable seismic ground motions for performance‐based design

The performance‐based design of lifeline systems requires spatially variable seismic excitations at the structures' supports that are consistent with prescribed seismic ground motion characteristics and an appropriate spatial variability model—such motions can be obtained through conditional simulation. This work revisits the concept of conditional simulation and critically examines the conformity of the generated motions with the characteristics of the target random field and observations from data recorded at dense instrument arrays. Baseline adjustment processing techniques for recorded earthquake accelerograms are extended to fit the requirements of simulated and conditionally simulated spatially variable ground motions. Emphasis is placed on the use of causal vs acausal filtering in the data processing. Acceleration, velocity and displacement time histories are evaluated in two example applications of the approach. The first application deals with a prescribed synthetic time history that incorporates nonstationarity in the amplitude and frequency content of the motions and depends on earthquake magnitude, source–site distance and local soil conditions; this example results in zero residual displacements. The second application considers as prescribed time history a recording in the vicinity of a fault and yields nonzero residual displacements. It is shown that the conditionally simulated time histories preserve the characteristics of the prescribed ones and are consistent with the target random field. The results of this analysis suggest that the presented methodology provides a useful tool for the generation of spatially variable ground motions to be used in the performance‐based design of lifeline systems. Copyright © 2006 John Wiley & Sons, Ltd.     

柴达木盆地东部地震地面运动放大效应

柴达木盆地是青藏高原东北部大型断陷山间盆地,该地区的流动观测记录了2008年11月10日发生于大柴旦附近的M_W6.3地震。和附近的基岩上的记录相比,盆地内部的记录显示出非常显著的地面运动放大效应,表现为峰值速度的增大、持续时间的延长,其呈现出长持续时间的后续震相。傅里叶频谱分析表明盆地内部显著的后续震相的频率和直达波相比较低,地面质点运动轨迹图显示后续震相为面波运动特征。为了解释地面运动的差异,构建二维模型,通过交错网格高阶有限差分方法计算了地震波在盆地内部的传播过程,结果显示盆地内部低速层的存在造成直达波的放大以及多次反射与转换,盆地边缘结构造成的波的相干叠加产生了强烈的次生面波,其低频、大振幅、长持续时间的特征是盆地内部地面运动放大的主要原因。     

利用区域地震波形振幅包络约束朝鲜地下核试验的埋深和当量

核爆当量和埋藏深度是地下核试验的两个重要参数.根据中国东北地区区域范围内地震台站观测记录,利用Pn,Pg,Sn和Lg波波形,计算水平分量的尾波振幅包络,调查2006年10月9日,2009年5月25日,2013年2月13日和2016年1月6日四次朝鲜地下核试验的爆炸当量和埋藏深度.以牡丹江(MDJ)台站记录为例,对两个水平分量波形进行带通滤波,计算平均波形的振幅包络.最终得到区域地震台站水平分量振幅包络,振幅稳定,包络振幅的变化清晰地显示区域震相的位置.区域震相的时间域振幅包络由震源谱函数、传播效应、台基响应和传递函数及尾波形状函数构成.通过网格搜索的方法,拟合水平分量记录的波形包络,可以获得核爆当量和埋藏深度的估计.结果表明,朝鲜四次地下核试验爆炸当量以时间为序从0.6±0.2kt到3.0±1.5kt,增加至10.0±2.0kt,再降到8.0±2.0kt.2006年核爆的埋藏深度较浅,为150±100m,2009年朝鲜核试验的埋藏深度为350±100m,2013年和2016年朝鲜核试验的深度非常一致,均为500±200m.这些结果与前人的调查结果一致性较高,说明使用单一地震台站时间域水平分量尾波振幅包络是可能同时约束地下核试验爆炸当量和埋藏深度的.     

相关多点地震动合成HOP方法的分析与改进

解释了相关多点地震动合成的HOP方法具有局部场地收敛性,分析了其生成各点地震动幅值相差较大的原因。除第一点外,HOP方法生成的各点地震动自功率谱围绕真实值波动,导致某些频率的能量过大而某些频率变小。由于结构频谱的离散性,如果恰好与能量大的频率吻合,则会激发结构强烈的振动响应,反之结构的响应值偏小,这与实际不符;进行非线性响应计算时或因为某些频率能量过大而使结构进入非线性阶段,则响应的期望值也是错误的。通过对实际地震波传播机制的模拟,提出了改进的HOP的方法。算例表明,此方法具有局部场地收敛性,各点自功率谱吻合给定的目标值,进而保证了各点幅值的一致性,可以用于结构抗震非线性响应分析。     

Lateral variations of the local magnitude at Ebre station,Northeastern Iberian Peninsula

Since 1984, a seismic surveillance at the Baix Ebre Region (NE of the Iberian Peninsula) has been carried out from a seismic network consisting of two digital and one analog stations. Large discrepancies between the recorded amplitudes and the magnitude estimates given by different agencies have been observed. To explain these discrepancies, that consist of a large reduction of amplitude and a strong dependence of the azimuth, a local magnitude formula for the area has been computed, using the analog records obtained at EBR short-period station. Global results show an average station correction factor between the computed local magnitude and the reported magnitude of the order of 1, presenting large variations around the mean value with a clear dependence on the azimuth. A detailed analysis of this variation around the mean value reveals that earthquakes can be classified according to their location in different geological units, thus obtaining for each zone a distinct station correction factor that can fluctuate as much as 0.6 around the average. These large lateral variations of the magnitude are in good agreement with previous studies of coda-Q attenuation.     

首都圈数字地震台网对微弱爆破信号的检测能力

利用首都圈数字地震台网接收人工地震信号,进行地下结构研究具有重要意义.但人工震源释放的能量小,激发的地震波以短周期为主,因此本文较全面地研究了地震台网对短周期微弱信号（1～20 Hz）的检测能力：(1) 分析了台网的背景噪声,结果表明基岩台址的地震台噪声比沉积盖层台址的地震台噪声低约13 dB,这相当于近1个震级的检测阈值;夜间的噪声比白天低约5 dB;噪声有逐年增高的趋势,2006年比2001年噪声提高约4 dB.(2 )分析了在台网内进行的药量为25 kg的陆地井下爆破实验,一次爆破相当于0.69级（M_L）的天然地震,有18个地震台可辨认爆炸产生的Pg、Pm或Pc波;离爆破点218 km的基岩台,仍可以接收到振幅只有1.6 nm 的Pm波,这个结果可为地震勘探实际工作提供参考.(3) 研究了台网外核爆试验的信号特征,2006年发生在朝鲜的地下核试验是一次检验台网检测微弱信号能力的好机会.波形记录经1～5Hz滤波后,台网中噪声小的18个基岩台可以清晰辨认核爆破产生的P波或Lg波,P波平均振幅为16 nm,计算的平均震级为m_b4.3,和NEIC给出的震级相同;分析还表明背景噪声是影响台站信号检测能力的主要因素之一.     

Variations of the lunar geomagnetic tide in the Indian region

The seasonal variations of the lunar geomagnetic tide are analysed at a group of stations around the longitude of India. The position of the focus of the lunar ionospheric current system is found to move with season and to disappear in northern winter. A year-to-year variability of the phases is also demonstrated. The phase relationships indicate the presence of antisymmetric atmospheric tidal modes of various strengths in different seasons. The different lunar Hough mode amplitudes, derived earlier, may also explain the observed September maximum in the lunar horizontal magnetic component at the equator.     

Increment of Macroseismic Intensity at Seismic Stations of the Kamchatka Region Relative to the Petropavlovsk Seismic Station

The paper presents the method and results of calculating the increment of macroseismic intensity at seismic stations of Kamchatka. Calculation is based on measurement of the relative level of maximum accelerations of intense earth vibrations in the phase of S-waves of comparatively strong regional earthquakes and the root-mean-square deviation of acceleration in the phase of P-waves of a strong distant earthquake. In the latter case, records of an earthquake with a magnitude of 9.1, which occurred in Japan on March 11, 2011, were used. The Petropavlovsk seismic station was used as the reference station. At the foundation of this station rests on rocky soil composed of siliceous shales. An estimate of the increment for the majority of digital stations is presented. Anomalously high intensity values were noted at a number of stations. The data obtained are used to assess the properties of soils in the investigated area. At several stations, the intensity of the horizontal component of soil vibrations above the intensity of the vertical component is much greater than the corresponding design value, which is probably due to the presence of resonant soil layers under these stations. The discrepancy in the incremental intensity estimates from records of intense oscillations from regional earthquakes and from records of a very strong remote earthquake obtained from sensors located in basements of heavy-frame concrete structures is revealed. To avoid distortion in recording ground vibrations, it is desirable to place seismic instruments far from such structures. The results obtained in the study can be used for seismic microzoning of construction sites in the investigated territory.     

Effects of lateral velocity heterogeneity under the Nevada Test Site on short-periodP wave amplitudes and travel times

Short-period teleseismicP waves from the Nevada Test Site (NTS) show systematic variations in amplitudes and travel times, with low amplitudes corresponding to fast travel times, suggesting elastic focussing-defocussing effects. Also, the azimuthal amplitude and travel time patterns for events at the Pahute Mesa subsite are systematically different from those at the Yucca Flat subsite, indicating the presence of a near-source component in both the amplitude and travel-time variations. This component is isolated by removing the mean station pattern for all of NTS from the observations. A very-near-source component in the Pahute Mesa observations is also isolated by removing subsite station means from the measurements, whereas the Yucca Flat observations exhibited no coherent very-near-source component. These anomalies are back-projected through laterally homogeneous structure to form thin lens models at various depths. Travel-time delays are predicted from the amplitude variations using the equation for wavefront curvature. The long-wavelength components of the predicted and observed time delays correlate well, at depths of 25 km for the very-near-source component under Pahute Mesa and 160 km for the regional component under NTS. The time delay surfaces predicted by the amplitudes at these depths are mapped into warped velocity discontinuities suitable for the calculation of synthetic seismograms using the Kirchhoff-Helmholtz integral formulation. Both the intersite (near-source) and intrasite (very-near-source) differences in amplitudes are qualitatively predicted very well, although the range of variation is somewhat underpredicted. This deficiency is likely due to the destructive interference of anomalies inherent in back-projection to a single layer.     

一种震相拾取误差的自动判别和计算方法

本文给出了一个主要用于深地震测深数据的震相识别误差（不确定性）的判别和计算方法.该方法集中讨论从记录截面拾取震相这一过程所引起的判别误差.以震相前后一定时窗内的地震记录振幅的均方根之比为判别依据,找出误差分布范围并给出走时误差与振幅比的分级相关函数.由此,当震相确定后,计算程序将根据记录数据自动算出识别误差.实践证明该方法不仅更加客观真实、方便快捷,而且为今后震相提取工作的进一步科学规范打下了基础.     

甘肃省测震台网地震台站地震计方位角检验与校正

利用2007年8月—2013年9月甘肃省测震台网记录的549个MS≥6远震事件的P波资料,对所属44个地震台站分量方位等问题进行检核计算。考虑到地震计维修、更换等因素,提供了甘肃省测震台网地震台站按时间序列计算的方位角偏差及方位校正变化情况,以保障测震台网数据的连续性与可靠性。研究结果表明,反演的方位偏差与文献[1]的结果基本一致,表明中国"十五"数字地震网络确实存在部分台站方位偏差较大等方面问题。因此,在进行现代地震学研究中应充分考虑台站地震计方位误差较大和研究时间段内方位角变动等因素的影响。     

Experimental analysis of seismic torsional ground motion recorded by the LSST‐Lotung array

The purpose of this study was to perform an experimental analysis of the amplitude of full‐scale spatial variability of seismic motions with regard to earthquake engineering. The LSST‐Lotung array in Taiwan provides a good set of records for this type of study. Of interest are the free‐field torsional seismic components induced by the spatial variability of seismic motions. In this study, three events have been considered: LSST‐06, LSST‐07 and LSST‐16. In time domain analysis, the experimental results obtained show that when the separation distance increases, the induced torsion decreases, and the normalized differential motion parameter increases. Also, the results show that the stronger the event, in terms of maximum PGA, the larger the induced torsional amplitude and the smaller the induced normalized differential motion parameter. Copyright © 2002 John Wiley & Sons, Ltd.     

Modelling of depth phases and source processes of some central Aleutian earthquakes

Short-period synthetic seismograms are computed to determine the relative amplitudes and arrival times of P, pP, pwP (water surface reflection), and sP phases. Except along nodal planes of upgoing p, pwP is of greater amplitude than sP. For central Aleutian earthquakes, pwP dominates over sP and pP at North American stations for thrust mechanisms and modelled crustal structures. The pwP phase is clearly identified in three Aleutian events and can be used to constrain focal depths. Complex fracturing processes are identified in two of the events. In simple events, smaller phases which are consistent from station to station can be identified as sub-surface reflections and used for modelling the structure of the forearc.     

Preliminary study of crust-upper mantle structure of the Tibetan Plateau by using broadband teleseismic body waveforms

As part of a joint Sino-U.S. research project to study the deep structure of the Tibetan Plateau, 11 broadband digital seismic recorders were deployed on the Plateau for one year of passive seismic recording. In this report we use teleseimic P waveforms to study the seismic velocity structure of crust and upper mantle under three stations by receiver function inversion. The receiver function is obtained by first rotating two horizontal components of seismic records into radial and tangential components and then deconvolving the vertical component from them. The receiver function depends only on the structure near the station because the source and path effects have been removed by the deconvolution. To suppress noise, receiver functions calculated from events clustered in a small range of back-azimuths and epicentral distances are stacked. Using a matrix formalism describing the propagation of elastic waves in laterally homogeneous stratified medium, a synthetic receiver function and differential receiver functions for the parameters in each layer can be calculated to establish a linearized inversion for one-dimensional velocity structure. Preliminary results of three stations, Wen-quan, Golmud and Xigatze (Coded as WNDO, TUNL and XIGA), located in central, northern and southern Plateau are given in this paper. The receiver functions of all three stations show clear P-S converted phases. The time delays of these converted phases relative to direct P arrivals are: WNDO 7.9s (for NE direction) and 8.3s (for SE direction), TUNL 8.2s, XIGA 9.0s. Such long time delays indicate the great thickness of crust under the Plateau. The differences between receiver function of these three station shows the tectonic difference between southern and north-central Plateau. The waveforms of the receiver functions for WNDO and TUNL are very simple, while the receiver function of XIGA has an additional midcrustal converted phase. The S wave velocity structures at these three stations are estimated from inversions of the receiver function. The crustal shear wave velocities at WNDO and TUNL are vertically homogeneous, with value between 3.5–3.6 km/s down to Moho. This value in the lower crust is lower than the normal value for the lower crust of continents, which is consistent with the observed strong Sn attenuation in this region. The velocity structure at XIGA shows a velocity discontinuity at depth of 20 km and high velocity value of 4.0 km/s in the midcrust between 20–30 km depth. Similar results are obtained from a DSS profile in southern Tibet. The velocity under XIGA decreases below a depth of 30 km, reaching the lowest value of 3.2 km/s between 50–55 km. depth. This may imply that the Indian crust underthrusts the low part of Tibetan crust in the southern Plateau, forming a “double crust”. The crustal thickness at each of these sites is: WNDO, 68 km; TUNL, 70 km; XI-GA, 80 km. The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,14, Supp., 581–592, 1992.     

Can tidal perturbations associated with sea level variations in the western Pacific Ocean be used to understand future effects of tidal evolution?

This study examines connections between mean sea level (MSL) variability and diurnal and semidiurnal tidal constituent variations at 17 open-ocean and 9 continental shelf tide gauges in the western tropical Pacific Ocean, a region showing anomalous rise in MSL over the last 20 years and strong interannual variability. Detrended MSL fluctuations are correlated with detrended tidal amplitude and phase fluctuations, defined as tidal anomaly trends (TATs), to quantify the response of tidal properties to MSL variation. About 20 significant amplitude and phase TATs are found for each of the two strongest tidal constituents, K₁ (diurnal) and M₂ (semidiurnal). Lesser constituents (O₁ and S₂) show trends at nearly half of all gauges. Fluctuations in MSL shift amplitudes and phases; both positive and negative responses occur. Changing overtides suggest that TATs are influenced by changing shallow water friction over the equatorial Western Pacific and the eastern coast of Australia (especially near the Great Barrier Reef). There is a strong connection between semidiurnal TATs at stations around the Solomon Islands and changes in thermocline depth, overtide generation, and the El Niño Southern Oscillation (ENSO). TATs for O₁, K₁, and M₂ are related to each other in a manner that suggests transfer of energy from M₂ to the two diurnals via resonant triad interactions; these cause major tidal variability on sub-decadal time scales, especially for M₂. The response of tides to MSL variability is not only spatially complex, it is frequency dependent; therefore, short-term responses may not predict long-term behavior.     

Seismic Wavefields in the Deep Seafloor Area from a Submarine Landslide Source

We use the finite difference method to simulate seismic wavefields at broadband land and seafloor stations for a given terrestrial landslide source, where the seafloor stations are located at water depths of 1,900–4,300 m. Our simulation results for the landslide source explain observations well at the seafloor stations for a frequency range of 0.05–0.1 Hz. Assuming the epicenter to be located in the vicinity of a large submarine slump, we also model wavefields at the stations for a submarine landslide source. We detect propagation of the Airy phase with an apparent velocity of 0.7 km/s in association with the seawater layer and an accretionary prism for the vertical component of waveforms at the seafloor stations. This later phase is not detected when the structural model does not consider seawater. For the model incorporating the seawater, the amplitude of the vertical component at seafloor stations can be up to four times that for the model that excludes seawater; we attribute this to the effects of the seawater layer on the wavefields. We also find that the amplification of the waveform depends not only on the presence of the seawater layer but also on the thickness of the accretionary prism, indicating low amplitudes at the land stations and at seafloor stations located near the trough but high amplitudes at other stations, particularly those located above the thick prism off the trough. Ignoring these characteristic structures in the oceanic area and simply calculating the wavefields using the same structural model used for land areas would result in erroneous estimates of the size of the submarine landslide and the mechanisms underlying its generation. Our results highlight the importance of adopting a structural model that incorporates the 3D accretionary prism and seawater layer into the simulation in order to precisely evaluate seismic wavefields in seafloor areas.     

地震自动定位的综合解决方案

针对已有两类地震自动定位方法各自的优点和限制，提出了综合利用这两类方法以形成一个统一系统的方案，希望以此提高地震自动定位能力。该方案的基本流程是：利用波形相关方法在连续地震波形上检测出地震事件并初定其震源位置和发震时刻；利用初定的震源位置和发震时刻，预测各台站各震相的到时；以预测的各震相到时为参考点，利用单台波形(单道或三分向)在该参考点附近精确测定震相到时；利用修定后的震相到时测定震源位置和发震时刻；根据震中距判定地震类型(地方震、近震、远震)，在预期的波形段自动测定波形最大振幅和相应周期，进而测定相应类型的震级(ML、Ms或mb)。     

Microtremors and volcanic explosions at Stromboli (Aeolian Islands)

Seismic observations were carried out at Stromboli from October 16^th to 24^th, 1972. One three-component seismic station was set up at Semaforo Nuovo, about 2 km on the E of the crater plain. Seismic observations showed the occurrence of microtremors, with variable amplitude but constant frequency content, and different types of explosion earthquakes. Spectral and vibration orbit analyses of microtremors and quakes were made in order to infer their nature and the depth of origin. Changes in microtremor amplitudes corresponding to different phases of volcanic activity indicate that their study may be useful in forecasting strongly explosive volcanic phases at Stromboli.