识别地震与爆炸的一种时间域新方法

利用地震观测波形直达纵横波能量的差异,在时间域创建标志地震事件震源性质的新参量a,计算182个台站记录的24个地震事件的a值,并用以分析2017年9月3日敏感地区发生的地震事件性质。结果表明,判据a具有良好的震源性质依附效应,从而为地震与爆炸的识别引进了一种新方法。     

震源性质的倒谱分析ht

介绍了倒谱的概念．从天然地震与人工爆炸的震源差异入手，推断震源差异在各域的反映形式，探索表征这种震源差异的有效途径．将震源识别的研究方法从时间域、频率域扩展到倒谱域，提出了识别震源性质的倒谱分析方法，建立了一个震源类型的识别判据．对近年来发生的一些地震及爆炸事件进行了倒谱分析，结果表明该方法具有良好的效能．      

用ISODATA方法对兰州——古浪地区潜在震源的识别

潜在震源区的判定是地震危险性分析工作的基础,潜在震源的准确性直接影响着地震危险性分析的结果的合理与否;本文采用多维实测数据类似性和差异性相对识别的ISODATA算法,对兰州——古浪地区的潜在震源进行了识别,给出了该地区的潜在震源区及其震级上限,为潜在震源区的综合判定提供了可靠的科学依据,使得潜在震源区的综合判定更为合理。     

近期国外震源研究综述

本文简要介绍了近年来国外在地震震源、台阵（网）资料在震源参数研究中的应用,全球震源谱对比、地震破裂传播方向、遗传算法在震源参数研究中的应用,地震矩、拐角频率、震源尺度、应力降与震源模型的探讨,以及确定初动震源机制的新方法等方面的研究情况和动态,并对其中的研究结果给予了一定的评述。     

利用全波形匹配方法确定水力压裂诱发地震震源机制

准确确定水力压裂诱发地震的震源机制对于描述裂缝破裂类型以及工区地应力状态十分重要.本文采用全波形匹配的方法确定诱发地震震源机制解,在拟合波形的同时对实际数据P波初动极性和纵横波振幅比也进行匹配,并在此基础上发展了一种基于邻域算法分级优化确定震源机制解的新方法.新方法的优势在于它对解具有更强的约束,并且能够解决现有方法中存在的权值系数选择问题.合成数据测试结果表明新方法能够得到可靠的震源机制解.本文将该方法应用于国内某页岩气田水力压裂诱发地震监测数据,反演结果表明诱发地震的震源类型以走滑破裂为主,表明该地区的最大水平主应力大于垂直主应力.     

云计算环境下震前震源异常次声波自动识别方法

云计算下采用三点阵次声源定位方法,在自动识别震前震源次声波过程中不能自动筛选识别大量的异常次声波数据,导致震前监测准确度不高且效率低下。因此提出一种云计算环境下震前震源异常次声波自动识别方法,构建JNS异常次声波数据采集筛查模组,全天候实时扫描访问端口,快速反馈异常次声波数据,采用NDS异常次声波数据序列异常检测算法快速识别错误序阵,准确回查、定位和锁定异常次声波数据;利用震前震源异常次声波自动识别方法分类识别异常次声波信号,判断该信号是否是地震可疑信号。实验结果表明,所提方法可有效自动识别震前震源异常次声波信号类型,信号分类准确率最大值达到99.99%;多次识别耗时最大均值仅为1.3min,具有准确率高和效率快的优势。     

一种矩张量反演新方法及其对2008年汶川Ms 8.0地震序列的应用

本文利用全球地震台网（GSN）的宽频带与长周期地震波形资料,采用矩张量反演的新方法,反演得到了2008年5月12日四川省汶川县MS8.0地震及其7个较大余震（MS5.0-6.0）的矩张量解与震源时间函数等震源参数.文章首先简要叙述矩张量反演新方法的理论背景和技术途径,并以汶川大地震的一个余震为例阐述了具体的实现过程;然后给出包括主震在内的8次地震的矩张量解和震源时间函数;最后分析探讨这些结果的构造意义.本文提出的矩张量反演新方法,不但与全球矩心矩张量（GCMT）一样,可以给出点源矩张量解,而且还可以给出点源的震源时间函数. 反演得到的汶川大地震的7个较大余震的震源时间函数表明,即使是中等强度的地震也可能有复杂的震源过程;汶川大地震的多数余震发生在以逆冲为主、兼具小量走滑分量的龙门山断裂带的主断裂上,但很可能有些余震则发生在主断裂附近的次级走滑断裂上.      

多种方法测定广东阳西3.7级地震震源深度

应用CAP和深度震相方法,对2018年3月20日发生在广东阳西的M3.7级地震震源深度进行了测定。首先通过CAP方法反演获得震源机制解,拟合最佳震源深度为12km。然后在震中距100km以内的近台识别出清晰的sPL和PmP、sPmP震相,利用频率-波数(F-K)方法,计算出深度震相在不同深度下的理论地震图,与实际观测波形对比测定震源深度为12km。再利用250~400km震中距范围内台站上识别出的sPn与Pn震相的走时差,测得震源深度12.6km。多种方法的研究结果一致,表明该地震震源深度为12km比较可靠。     

区域地震事件的纵横波差异与震源类型

简要介绍地震与爆炸的震源物理性质及地震波激发机制特征。利用地震观测波形直达纵横波震相能量释放的差异,在时间域引进标志事件震源性质的参量a。132个台站记录的22个不同震源类型事件的计算结果表明,a值具有震源性质依附效应,展现出在事件识别领域的应用前景。     

四川攀枝花2008年8月31日Ms 5.6地震震源深度的确定

利用四川区域数字地震台网记录,采用滑动时窗相关法识别sPn震相,测定了四川攀枝花2008年8月31日Ms 5.6地震的震源深度.研究结果显示,震源深度为17 km,表明震源位于中上地壳.元谋-绿汁江断裂带上中强地震分布与断裂走向吻合,震源机制解以走滑为主,这些表明此断裂近乎直立.结合此次地震的震源深度,我们推测此断裂可能达到20 km深.     

地震波双向照明的概念及计算方法

根据WRW模型理论,从地震波传播的物理过程给出了双向照明的具体概念,基于该概念提出了计算地震波双向照明强度的基本思路(称为DUC双向照明方法).为了有效提高DUC双向照明方法的计算效率,又提出了DUC-DC双向照明方法,并将提出的上述方法与传统的双向照明计算方法从计算精度和计算效率等方面进行了对比.试验结果表明,本文提...     

基于共姿态道集的静校正方法

静校正问题是地震勘探的关键问题,直接影响地震勘探精度和准确性.实际地震采集过程中,当在相同接收点位置上不同时间内插拔布设了不同的检波器时,对于目前基于地表一致性理论假设的基准面静校正和剩余静校正,以及非地表一致性剩余静校正都不具备适用条件.为解决这一问题,本文提出了基于共姿态道集的静校正方法,将相同接收点位置上不同时间布设的检波点所接收的地震数据抽成不同的共姿态道集,在共姿态道集内实施地表一致性静校正;当某接收点位置上具有若干个共姿态道集时,该接收点位置上可能会存在多个检波点静校量;炮点静校问题仍然采用地表一致性静校正方法解决.该方法解决了同一接收点位置上不同共姿态道集之间的非地表一致性静校正问题,同时也解决了全区的检波点和炮点的地表一致性静校正问题,在实际数据应用效果明显.     

炸药震源的地震安全性野外实验研究

选定一块场地,针对相同井深、不同药量的情况,进行炸药震源的地震安全性野外实验,用Etna数字强震动加速度仪记录每次爆破时房屋振动的加速度.在选定房屋结构的地基、窗台、屋顶等3个位置分别布设了仪器,记录结构响应的加速度波形.为了与《爆破安全规程》的参考标准对应,将加速度值转换成速度值,并对位于屋顶的结构响应速度峰值进行分析.通过对结构响应、安全距离、频率的研究,分析建筑物结构对每次爆破的响应,并对其安全性进行讨论.     

Detection of seismic reflections from seismic attributes through fractal analysis

The concept of fractals is used here for the identification of seismic reflectors with special emphasis on thin‐bed delineation, which is generally overlooked during standard data processing. A new fractal analysis scheme is applied to both synthetic and real field seismic data. The fractal dimensions of the three seismic attributes – amplitude, phase, and instantaneous frequency – have been analysed and evaluated. A change in fractal dimension is found to occur whenever there is a reflection. However, the resolution in the delineation of reflectors varies, depending on the attribute under consideration and the method of fractal dimension estimation used. Fractal analysis is performed on both noise‐free and noisy synthetic data to establish the noise tolerance limit for both the ‘divider method’ and the ‘Hurst method’. It is then tested with different peak frequencies of the source wavelet to establish the criteria for using the divider method and the Hurst method. The divider method is found to be suitable for high peak frequency source wavelets (> 25 Hz), while the Hurst method is best suited for low peak frequency source wavelets (< 25 Hz). Finally, when applied to the digitally processed and migrated field seismic data, it could even delineate reflectors which otherwise went undetected on the migrated time section.     

Separation method for multi-source blended seismic data

Multi-source seismic technology is an efficient seismic acquisition method that requires a group of blended seismic data to be separated into single-source seismic data for subsequent processing. The separation of blended seismic data is a linear inverse problem. According to the relationship between the shooting number and the simultaneous source number of the acquisition system, this separation of blended seismic data is divided into an easily determined or overdetermined linear inverse problem and an underdetermined linear inverse problem that is difficult to solve. For the latter, this paper presents an optimization method that imposes the sparsity constraint on wavefields to construct the object function of inversion, and the problem is solved by using the iterative thresholding method. For the most extremely underdetermined separation problem with single-shooting and multiple sources, this paper presents a method of pseudo-deblending with random noise filtering. In this method, approximate common shot gathers are received through the pseudo-deblending process, and the random noises that appear when the approximate common shot gathers are sorted into common receiver gathers are eliminated through filtering methods. The separation methods proposed in this paper are applied to three types of numerical simulation data, including pure data without noise, data with random noise, and data with linear regular noise to obtain satisfactory results. The noise suppression effects of these methods are sufficient, particularly with single-shooting blended seismic data, which verifies the effectiveness of the proposed methods.     

Moment magnitude determination of local seismic events recorded at selected Polish seismic stations

The paper presents the method of local magnitude determination used at Polish seismic stations to report events originating in one of the four regions of induced seismicity in Poland or its immediate vicinity. The method is based on recalculation of the seismic moment into magnitude, whereas the seismic moment is obtained from spectral analysis. The method has been introduced at Polish seismic stations in the late 1990s but as of yet had not been described in full because magnitude discrepancies have been found between the results of the individual stations. The authors have performed statistics of these differences, provide their explanation and calculate station corrections for each station and each event source region. The limitations of the method are also discussed. The method is found to be a good and reliable method of local magnitude determination provided the limitations are observed and station correction applied.     

一种地震P波和S波初至时间自动拾取的新方法

地震P波、S波初至时间的拾取是地震波分析的一项基础性工作.本文提出了一种新的地震波初至时间自动拾取的方法:首先,把地震波的三分量时程曲线变换为一组空间向的能量变化率时程曲线;然后对能量变化率时程曲线进行STA/LTA(Short Time Average/Long Time Average,短时间的均值/长时间的均值)处理,拾取地震P波和S波的大致初至时间;最后提出采用一种二次方自回归模型对初至附近的能量变化率曲线进行二次方自回归处理,精确拾取出P波和S波的初至时间.本文采用了10组芦山地震的记录数据和150组汶川地震的记录数据对此方法的可靠性进行了检验.以人工拾取结果为参考,此方法具有很高的准确率和稳定性,同时,相比于常用的STA/LTA方法和AIC(Akaike Information Criterion,Akaike信息准则)方法,此方法在计算时间效率方面稍微逊色,但是对S波初至时间的拾取精度和可靠性更高.此方法丰富了地震P波、S波初至时间的自动拾取方法.     

龙门山和成都地震构造区的划分

新编制的地震动参数区划图采用了潜在震源区三级划分方案,以体现背景地震活动空间分布的不均匀性,并在地震构造区内归纳出统一的地震构造模型.本文根据西南地区潜在震源区三级划分的成果,分析了龙门山地震统计区内的龙门山和成都地震构造区的基本特征,历史地震活动强度及频度,主要活动构造的构造变形样式,建立了地震构造区的发震构造模型,确定了构造区的本底地震及划分构造源的地震构造标志.同时,提出了确定背景源空间分布函数的简单方法.     

Making the transition from discrete shot records to continuous seismic records and source wavefields,and its potential impact on survey efficiency and environmental footprint

A marine seismic method based on continuous source and receiver wavefields has been developed. The method requires continuous recording of the seismic data. The source that may consist of multiple source elements can emit signals continuously while moving. The ideal source wavefield to be used with this method should be as white as possible both in a temporal and a spatial sense to avoid deep notches in the spectrum enabling a stable multi-dimensional deconvolution. White noise has such properties. However, equipment that can generate white noise does not exist. In order to generate a continuous source wavefield that is approaching the properties of white noise using existing equipment onboard marine seismic vessels, individual air-guns can be triggered with short randomized time intervals in a near-continuous fashion. The main potential benefits with the method are to reduce the environmental impact of marine seismic surveys and to improve acquisition efficiency. The peak sound pressure levels are significantly reduced by triggering one air-gun at a time compared to conventional marine seismic sources. Sound exposure levels are also reduced in most directions. Since the method is based on continuous recording of seismic data and the air-guns are triggered based on time and not based on position, there are less vessel speed limitations compared to conventional marine seismic data acquisition. Also, because the source wavefield is spread out in time, the wavefields emitted from source elements in different cross-line positions can be designed such that the emitted wavefield is spatially white in this direction. This means that source elements in multiple cross-line positions can be operated simultaneously, potentially improving the cross-line sampling and/or the acquisition efficiency.     

判定潜在震源的模式识别方法的改进

本文给出了在非均匀网格划分的情况下,利用改进的连续Hamming方法判定潜在震源的方法和计算过程。实际计算结果表明,用此方法计算,避免了由于潜在震源形状和空间参数分配所带来的一些问题,更适用于各种形状的构造和各种形状的网格。