Application of the method of dynamic calibration for IMS stations in regions with lower seismic activity

The Method of Dynamic Calibration (MDC) of stations of the International Monitoring System (IMS) was developed for calibrating regions where no underground nuclear explosions were carried out, with the purpose of providing conditions for implementation of the Comprehensive Nuclear-Test-Ban Treaty (CTBT) in nontrivial cases. Initially, the MDC had been presented in [Kedrov, 2001; Kedrov et al., 2001; Kedrov O.K. and Kedrov E.O., 2003] and then considered in detail in [Kedrov et al., 2008]. The core of MDC relates to adapting diagnostic parameters for the identification of underground nuclear explosions (UNE) and earthquakes elaborated for the region of Eurasia, taken as a basic region (BR), for other researched regions that differ from BR in the character of the attenuation of seismic waves. The unique characteristic of this method lies in the fact that calibration of diagnostic parameters with the help of attenuation coefficients b _Δ at varied source-station traces is implemented using only natural seismicity data within the limits of an explored region and does not require special underground chemical explosions. The MDC algorithm is implemented in the research program ”Kalibr”, which was tested by using the experimental data from Eurasia region. It is shown in this work that MDC can be used for calibration of regions where a very low level of natural seismicity is observed. According to the results of the calibration of diagnostic parameters at IMS stations in several regions of North America, Africa, and Asia, the approximate classification of propagation conditions for seismic signals at source-station traces in platform and tectonically active regions is made. The results for the development of two research programs, “Spektr” and “Signal”, are presented; this software is intended for automation of calculation procedures for spectral diagnostic parameters of UNEs’ and earthquakes’ identification by amplitude spectra of P waves and by the maximal amplitudes of P, S, and LR signals. The application of these programs allowed us to accelerate the whole calibration procedure for a particular source-station trace using the ”Kalibr” program.     

形变台站数据处理系统在C语言下的实现方法

介绍了地震前兆（变形）台站日常数据处理系统的编程方法、新思路及一些具体实现过程。     

Application of array data processing techniques to the swedish seismograph stations

Summary An efficient tool in seismic signal analysis is offered by the array data processing techniques which are based upon similarity of pulse shapes between stations. Records at the six Swedish stations in general exhibit the signal resemblance required by the array methods. By applying a simple processing technique, the signal-to-noise ratio is increased by a factor of about 2, and the accuracy ofP-wave readings is improved. High similarity is found not only between the same phase (P, pP) at different stations, but also between different phases (P, pP) at the same or different stations. Altogether, records of eleven earthquakes were investigated in detail. The Swedish station network can be considered as a super-large aperture seismic array (SLASA), despite some differences in instrumental equipment. A general adaptation of array techniques to SLASA systems is discussed. Such methods will greatly improve the seismological information which can be extracted from ordinary station records, provided the stations are located on geologically homogeneous ground.     

网上前兆数据动态时间变化无人值守警视系统

在Windows环境下开发了应用软件“网上前兆数据动态时间变化无人值守警视系统”。该软件基于安徽省前兆数据传输网上的14个台站、各学科(117个测项)前兆数据格式等特点,随时间的递增自动循环完成检索、汇集、处理、初级分析、绘制和显示时间变化曲线、报告缺测缺记时效等任务。能快捷、可视化地监知各台站、各学科前兆数据的时间变化情况。     

微震监测系统标定信号的应用

徐州矿务集团有限公司三河尖煤矿建立了覆盖全矿井的KZ-1型微震监测系统,通过地面监测主机向井下各测站发时间同步校时信号,用GPS秒号起始去同步各个测站下位机时钟,保证了各测站之间相对时间精确度.KZ-1型微震监测系统的拾震器设置了双换能器,且各有一个独立的磁回路,通过地面监测主机向拾震器标定换能器发标定指令信号,信号检测换能器再将输出信号上传至监测主机,实现了即时定量标定.标定信号的应用使KZ-1型微震监测系统具备整个系统检测功能,包括灵敏度、频响特征以及各测站的时间同步误差.     

地震活动因子预测方法在华北地区的应用研究

为探索地震活动性前兆异常图像在1年尺度上的时空演化特征,应用地震活动因子预测方法,对华北地区进行了全时空扫描分析,筛选出对应效果较好的4个因子进行了全程统计。结果表明,震前1年左右扫描区域内异常小区由多到少到发震的现象,以及发震构造单元的异常小区所占比例较高的现象。归纳出根据异常小区数从多-少-发震的时间预测方法和根据构造单元的异常小区高比值的地点预测方法,它是利用地震活动因子进行预测的更易于识别的时空异常特征指标,经R值评分检验在一定程度上提高了预报效能。     

用多通道AD分级采集扩展地震数据采集器的动态范围

地震数据采集是地震信号数字化必不可少的环节,动态范围是其一个重要的性能指标.实际地震信号的动态范围在160 dB以上,而目前普遍使用的24位地震数据采集器动态范围相对较小且在50 Hz采样率时最大只达到135 dB,致使24位地震数据采集器在实际使用中对小信号分辨率不够,不能有效提取地震信息;在大地震时又容易使数据采集器出现饱和限幅失真的现象而失去地震监测记录功能.本文针对在地震监测和地震研究中需要具有高分辨率和高动态范围的地震数据采集器这个亟待解决的问题,提出一种采用多通道AD转换器并行分级采集的方法,讨论了通道间失配及其标定.对研制实验样机的测试表明,其动态范围在50 Hz采样时可以达到157 dB以上,线性度优于0.005%.     

永安无人值守地震台远程网络监控系统

针对永安无人值守地震台设备维护及安全巡视等问题,利用地震台现有网络搭建方式,使用网络摄像机及相应监控系统软件,通过合理配置,实现无人值守地震台远程网络视频监控。该系统的实现,提高了永安地震台管控无人值守台站的能力,对其他地震台有一定借鉴意义。     

Locations and magnitudes of earthquakes in Central Asia from seismic intensity data

We apply the Bakun and Wentworth (Bull Seism Soc Am 87:1502–1521, 1997) method to determine the location and magnitude of earthquakes occurred in Central Asia using MSK-64 intensity assignments. The attenuation model previously derived and validated by Bindi et al. (Geophys J Int, 2013) is used to analyse 21 earthquakes that occurred over the period 1885–1964, and the estimated locations and magnitudes are compared to values available in literature. Bootstrap analyses are performed to estimate the confidence intervals of the intensity magnitudes, as well as to quantify the location uncertainty. The analyses of seven significant earthquakes for the hazard assessment are presented in detail, including three large historical earthquakes that struck the northern Tien-Shan between the end of the nineteenth and the beginning of the twentieth centuries: the 1887, M 7.3 Verny, the 1889, M 8.3 Chilik and the 1911, M 8.2 Kemin earthquakes. Regarding the 1911, Kemin earthquake the magnitude values estimated from intensity data are lower (i.e. MI_LH?=?7.8 and MI_W?=?7.6 considering surface wave and moment magnitude, respectively) than the value M?=?8.2 listed in the considered catalog. These values are more in agreement with the value M _S?=?7.8 revised by Abe and Noguchi (Phys Earth Planet In, 33:1–11, 1983b) for the surface wave magnitude. For the Kemin earthquake, the distribution of the bootstrap solutions for the intensity centre reveal two minima, indicating that the distribution of intensity assignments do not constrain a unique solution. This is in agreement with the complex source rupture history of the Kemin earthquake, which involved several fault segments with different strike orientations, dipping angles and focal mechanisms (e.g. Delvaux et al. in Russ Geol Geophys 42:1167–1177, 2001; Arrowsmith et al. in Eos Trans Am Geophys Union 86(52), 2005). Two possible locations for the intensity centre are obtained. The first is located on the easternmost sub-faults (i.e. the Aksu and Chon-Aksu segments), where most of the seismic moment was released (Arrowsmith et al. in Eos Trans Am Geophys Union 86(52), 2005). The second location is located on the westernmost sub-faults (i.e. the Dzhil'-Aryk segment), close to the intensity centre location obtained for the 1938, M 6.9 Chu-Kemin earthquake (MI_LH?=?6.9 and MI_W?=?6.8).     

应变能积累分析方法在冀鲁豫交界区地震活动趋势分析中的应用

利用应变能积累分析方法对冀鲁豫交界区地震活动特征进行分析研究,得到该区地震活动的一些规律,并对该区未来地震活动趋势进行分析判断,认为该区当前处于应变能积累一释放阶段后期,剩余能量较少,因此未来几年内,发生大地震的可能性很小,发生的最大地震应为中等或稍强地震。     

Application of array data processing techniques to a network of ordinary seismograph stations

Summary Array data processing techniques, which are based upon signal similarity between the stations, offer important tools in seismic signal analysis.P waves recorded in the Fennoscandian network in general exhibit the signal resemblance required for applying array processing methods. The largest station separation in this network is about 1800 km, and the site structures range from sedimentary rocks in Denmark to the homogeneous crust of the Baltic shield. Inother words, we can consider the Fennoscandian station network as a super-large aperture seismic array (SLASA). An analysis of the background noise shows that application of even very simple processing techniques results in very effective noise suppression.A general adaptation of array processing techniques to ordinary seismograph stations is considered possible. Such methods will greatly improve our ability to extract reliable information from ordinary seismograms. Although theC-D (compression-dilatation) readings are usually correct, they often do not refer to the very first signal onset. Records of nuclear explosions from stations all over the world exhibit a high degree of signal similarity.     

辽宁遥测数字地震台网和有人值守台站数据联合定位的方法研究

1　问题的提出辽宁省地震台网由辽宁遥测数字地震台网 (以下简称遥测台网 )和有人值守测震台 (下称台站 )组成。遥测台网由分布在全省各地的 1 5个遥测地震台站组成 ,但这些台站分布不够均匀 ,使得一些地区监测能力较弱。台站有的和遥测台在一起 ,有的则正好分布在没有遥测台的地区。为能在使用交互分析处理系统进行地震定位时 ,充分利用这些有人值守台的数据 ,使台网分布不够均匀的问题在一定程度上得到弥补 ,使地震定位精度得到进一步提高 ,本项目编制了面向台站速报的 LN-LHDW-1检测管理软件 ,以实现遥测台网与台站地震数据的联合定位…     

Application of mechanical and statistical models to the study of seismicity of synthetic earthquakes and the prediction of natural ones

Based on a cellular automata model, which obeys rules produced by a spring-damper-slider model of fault network, a synthetic earthquake catalog of deterministic chaos can be produced. It is a challenge to inverse the model parameters and predict the synthetic earthquakes from the synthetic catalog. We apply an improved statistical model of coupled stress release to the synthetic catalog, it does successfully inverse the stress evolution of the system, but effectively prediction of the earthquakes from the stresses remains a problem. This research is partially supported by the Chinese Joint Seismological Science Foundation and Asian 2000 Foundation of New Zealand.     

利用地脉动信息约束沉积层区域台站下方速度结构反演

利用接收函数和地震面波频散联合反演台站下方速度结构,已成为一种常用技术.联合反演结果能同时匹配两个不同数据集合,使得解的非唯一性在一定程度上得到了有效抑制.然而对于现有地震台站分布,短周期面波由于受多种因素的影响,可以获得的有效频散资料较少,层析成像的横向分辨率较低,难以反映台站附近局部介质真实的频散特性,而且这些面波...     

Temporal variations of seismicity parameters in the Central Alborz,Iran

Temporal changes of b-value, fractal (correlation) dimensions of epicenters (D _e ) and occurrence time of earthquakes (D _t ) and relations between these parameters were calculated to investigate precursory changes before 28 May 2004, Baladeh-Kojour earthquake (M _w = 6.3) of Central Alborz, Iran. 2086 events with M _N ≥ 1.7 were selected for our analyses. A wide range of variation was seen in these parameters: b-value ～ 0.6–1.11, D _e ～ 0.97–1.64, and D _t ～ 0.13–0.93. The results showed decreases in all fractal parameters several months before the main shock. This decrease, which might have arisen due to clusters of events occurred between 2002–2003, was followed by a systematic increase, corresponding to the increased level of low-magnitude seismicity. It seems that changes in fractal parameters may be precursors of Baladeh-Kojour earthquake which was caused by seismic activation and quiescence. Furthermore, a positive correlation between b-value and D _e was detected before the main shock (D _e = 0.87 + 0.7b) and during aftershock sequences (D _e = 2b ± 0.09), which was further on changed to a negative one (D _e = 2.56–1.32b).     

双极化Sentinel-1数据在昆明市沉降监测中的应用

针对多数学者用VV极化的Sentinel-1数据对城市区域进行监测而忽略VH极化数据的问题,本文利用SBAS-InSAR方法,结合28景双极化(VV极化和VH极化)Sentinel-1数据,得到了昆明市2018年8月—2019年7月间的地面沉降情况,并且对两种极化的沉降监测结果进行对比.结果表明昆明市主城区存在着较为明显的沉降现象,部分区域沉降严重形成明显的沉降漏斗.最大年平均沉降速率达到-53.63 mm/a(VV极化)和-53.67 mm/a (VH极化).对两种极化的数据集进行对比,发现在SDFP点的筛选和获取方面,VV极化数据集要优于VH极化数据集;其次,VV极化数据和VH极化数据的沉降速率、沉降分布及趋势等结果呈现出一致性,但是二者在同一时间段的各个监测节点的沉降时序变化存在着差异.     

Constrained tomography of realistic velocity models in microseismic monitoring using calibration shots

The knowledge of the velocity model in microseismic jobs is critical to achieving statistically reliable microseismic event locations. The design of microseismic networks and the limited sources for calibration do not allow for a full tomographic inversion. We propose optimizing a priori velocity models using a few active shots and a non‐linear inversion, suitable to poorly constrained systems. The considered models can be described by several layers with different P‐ and S‐wave velocities. The velocities may be constant or have 3D gradients; the layer interfaces may be simple dipping planes or more complex 3D surfaces. In this process the P‐ and S‐ wave arrival times and polarizations measured on the seismograms constitute the observed data set. They are used to estimate two misfit functions: i) one based on the measurement residuals and ii) one based on the inaccuracy of the source relocation. These two functions are minimized thanks to a simulated annealing scheme, which decreases the risk of converging to a local solution within the velocity model. The case study used to illustrate this methodology highlights the ability of this technique to constrain a velocity model with dipping layers. This was performed by jointly using sixteen perforation shots recorded during a multi‐stage fracturing operation from a single string of 3C‐receivers. This decreased the location inaccuracies and the residuals by a factor of six. In addition, the retrieved layer dip was consistent with the pseudo‐horizontal trajectories of the wells and the background information provided by the customer. Finally, the theoretical position of each calibration shot was contained in the uncertainty domain of the relocation of each shot. In contrast, single‐stage inversions provided different velocity models that were neither consistent between each other nor with the well trajectories. This example showed that it is essential to perform a multi‐stage inversion to derive a better updated velocity model.     

Induced seismicity in large-scale mining in the kola peninsula and monitoring to reveal informative precursors

Large volumes of rock mass, mined-out and moved within these deposits, resulted in irreversible changes in the geodynamic regime in the upper earth's crust of the adjacent territory. These changes manifest themselves in a more frequent occurrence of such intensive dynamic phenomena as tectonic rock bursts due to fault movement adjacent to the area which is mined-out and man-made earthquakes which sharply decrease mining safety and result in great material losses.To develop the prediction techniques of such phenomena, a monitoring system is created, based on the program of the Kola Complex of geodynamic measuring stations. Most of this system is realized in the region of the Khibiny apatite mines. The system provides regional seismological monitoring, local prediction of seismicity in separate areas of a rock mass and, determination of stress and strain in rock masses, local geophysical monitoring over the state of rocks in a rock mass as well as physical and mathematical modelling of geodynamic processes in the upper earth's crust.The investigations have resulted in the distinguishing of some regularities in manifestations of induced seismicity and tectonic rock bursts and in the determination of strain precursors of intensive seismic events in the Khibiny mines.The mechanism is provided by the induced seismicity which resulted from the anthropogenic impact on the geological medium. A geodynamic monitoring complex is described, which is used to reveal the precursors of powerful seismic eventsin situ, and monitoring results are shown, obtained in the Kola Complex of geodynamic stations. Methods of preventing tectonic rock bursts and induced earthquakes are presented.     

上海测震、地球物理台站隔离电源系统应用与改进

上海测震、地球物理台站运行中,因电源系统和雷电而导致的故障占了较大比例。为提升台站运行质量而研发的隔离电源系统,可有效地将市电与台站设备间进行物理隔离,杜绝了雷击通过交流电供电途径的串入。该系统研制成功后,在上海测震台站试运行,其间发现了系统的一些缺陷,技术人员进行多次改进,提升了系统的整体性能。该系统已在多个测震、地球物理台站应用,取得了显著的防雷、防浪涌效果。