2013年四川芦山7.0级地震烈度遥感评估

2013年4月20日四川芦山MS7.0级地震发生后,在灾区应急获取了多种高分辨率航空和无人机遥感影像,并快速解译提取了灾区建筑物震害信息.采用地震烈度遥感定量评估方法,利用2008年汶川8.0级地震等震后震害遥感解译和现场调查研究确定的经验震害遥感定量评估模型,获得了芦山地震灾区126个主要居民点的地震烈度遥感评估结果,并据此圈画了地震烈度分布遥感评估图.结果显示,本次地震Ⅸ度区面积约150km2,Ⅷ度区面积约900km2.该结果在第一时间(4月21日晚)提供给了中国地震局地震现场应急指挥部.对比分析显示,地震烈度遥感快速评估结果与中国地震局4月25日公布的地震烈度图,以及与笔者在现场实地进行的建筑物震害详细调查结果基础上评定的地震烈度具有较高的一致性.表明强烈地震发生后,借助于快速获取的灾区高分辨率遥感影像,可以快速估计地震烈度分布,对地震灾区灾情估计和抗震救灾工作具有十分重要的参考意义.     

利用三维高斯射线束成像进行地震定位

常规的地震定位方法通常需要拾取地震记录的初至,当初至不明显或被较高水平的噪声淹没时精度较低.本文采用基于三维高斯射线束的偏移成像方法对震源进行定位,较好地解决了该问题.通过三维高斯射线束对台站记录进行偏移归位,并将各台站成像结果的交点作为地震能量释放的中心位置;当各台站成像结果不能交于一点时,采用三维空间高斯滤波方法可实现震源位置的自动获取.提出的变网格计算方案极大地减少了计算量,显著地提高了成像精度和计算效率.利用首都圈地震台网数据,对涿鹿、滦县以及房山三个地震事件进行试算,结果表明:基于变网格三维高斯束偏移成像的地震定位方法自动化程度很高,而且具有较好的抗噪能力,特别适合处理低信噪比资料的地震定位问题.     

礼县—罗家堡断裂晚第四纪活动特征:兼论1654年礼县8级地震孕震机制

1654年7月21日甘肃省礼县8级地震发生在南北地震带的中北段,该地区的构造变形和构造活动与青藏高原向北东方向的扩展密切相关,复杂的构造几何特征主要受控于东昆仑断裂、西秦岭北缘断裂和一系列北东向断裂.礼县—罗家堡断裂为一条北东东向的左旋走滑活动断裂,错断了含有仰韶文化红色陶瓷片的一级阶地堆积物,阶地面上断层陡坎高约1.5m.沿断裂带发现冲沟的左旋位错量为3~10m,晚更新世黄土中残留的断层陡坎高4.5~8m.其中两条冲沟中发现裂点,高3.5~3.9m,距断层陡坎的距离为16m.礼县—盐关—罗家堡—天水镇一带发育大量的滑坡,长轴走向与礼县—罗家堡断裂一致,滑坡体后缘、侧壁陡峭,出露晚第四纪黄土,鲜有植被覆盖.礼县—罗家堡断裂为1654年8级地震的发震断裂.综合分析认为,受青藏高原向北东方向的扩展,被西秦岭北缘断裂、礼县—罗家堡断裂和岷县—宕昌断裂围限的礼县次级地块向南东滑动可能导致了1654年礼县8级地震的发生.断裂北、南两侧地壳电性结构的差异为强震的孕育提供了深部构造条件.     

GNMF小波谱分离在地震勘探噪声压制中的应用

地震勘探资料噪声压制及信噪比提高是整个地震勘探信号处理过程中的重要任务,随着地震勘探深度的增加及其复杂性,人们对地震数据质量的要求越来越高.勘探环境的复杂化使得采集到的地震资料中有效信号被大量噪声淹没,无法清晰辨识,严重影响后续的数据处理与解释.小波去噪是地震勘探中常用且发展较成熟的一种方法,但是其涉及到的阈值函数选取问题一直令人困扰,虽然已有多种阈值函数被提出,但仍存在各自的缺陷.本文利用小波分解在时域及频域良好的信号细节体现特性,引入模式识别中的非负矩阵分解(NMF)谱分离思想,针对小波系数阈值优化问题,提出了一种小波域图非负矩阵分解(GNMF)消噪算法.该方法首先在小波分解基础上,利用GNMF算法实现小波分解系数谱中信号分量与噪声分量的谱分离,然后通过反变换重构各分离子谱对应的子信号,最后利用K均值聚类算法将得到的多个子信号划分为信号类及噪声类,最终得到重构信号及分离噪声.合成记录和实际地震资料的消噪结果验证了新方法在提高信号与噪声分离准确性和精度方面的有效性,同时新方法避免了阈值选取造成的噪声压制不理想或有效成分损失问题.与小波消噪结果的对比及数值分析也说明了新方法在噪声压制及有效成分保持方面的优势.     

Sensitivity of the elastic anisotropy and seismic reflection amplitude of the Eagle Ford Shale to the presence of kerogen

The Eagle Ford Shale of Central and South Texas is currently of great interest for oil and gas exploration and production. Laboratory studies show that the Eagle Ford Shale is anisotropic, with a correlation between anisotropy and total organic carbon. Organic materials are usually more compliant than other minerals present in organic‐rich shales, and their shapes and distribution are usually anisotropic. This makes organic materials an important source of anisotropy in organic‐rich shales. Neglecting shale anisotropy may lead to incorrect estimates of rock and fluid properties derived from inversion of amplitude versus offset seismic data. Organic materials have a significant effect on the PP and PS reflection amplitudes from the Austin Chalk/Upper Eagle Ford interface, the Upper Eagle Ford/Lower Eagle Ford interface, and the Lower Eagle Ford/Buda Limestone interface. The higher kerogen content of the Lower Eagle Ford compared with that of the Upper Eagle Ford leads to a negative PP reflection amplitude that dims with offset, whereas the PS reflection coefficient increases in magnitude with increasing offset. The PP and PS reflection coefficients at the Austin Chalk/Upper Eagle Ford interface, the Upper Eagle Ford/Lower Eagle Ford interface, and the Lower Eagle Ford/Buda Limestone interface all increase in magnitude with increasing volume fraction of kerogen.     

Retrieval of super‐virtual refraction by cross‐correlation

Topography and severe variations of near‐surface layers lead to travel‐time perturbations for the events in seismic exploration. Usually, these perturbations could be estimated and eliminated by refraction technology. The virtual refraction method is a relatively new technique for retrieval of refraction information from seismic records contaminated by noise. Based on the virtual refraction, this paper proposes super‐virtual refraction interferometry by cross‐correlation to retrieve refraction wavefields by summing the cross‐correlation of raw refraction wavefields and virtual refraction wavefields over all receivers located outside the retrieved source and receiver pair. This method can enhance refraction signal gradually as the source–receiver offset decreases. For further enhancement of refracted waves, a scheme of hybrid virtual refraction wavefields is applied by stacking of correlation‐type and convolution‐type super‐virtual refractions. Our new method does not need any information about the near‐surface velocity model, which can solve the problem of directly unmeasured virtual refraction energy from the virtual source at the surface, and extend the acquisition aperture to its maximum extent in raw seismic records. It can also reduce random noise influence in raw seismic records effectively and improve refracted waves’ signal‐to‐noise ratio by a factor proportional to the square root of the number of receivers positioned at stationary‐phase points, based on the improvement of virtual refraction's signal‐to‐noise ratio. Using results from synthetic and field data, we show that our new method is effective to retrieve refraction information from raw seismic records and improve the accuracy of first‐arrival picks.     

Point‐spread functions for interferometric imaging

Interferometric redatuming is a data‐driven method to transform seismic responses with sources at one level and receivers at a deeper level into virtual reflection data with both sources and receivers at the deeper level. Although this method has traditionally been applied by cross‐correlation, accurate redatuming through a heterogeneous overburden requires solving a multidimensional deconvolution problem. Input data can be obtained either by direct observation (for instance in a horizontal borehole), by modelling or by a novel iterative scheme that is currently being developed. The output of interferometric redatuming can be used for imaging below the redatuming level, resulting in a so‐called interferometric image. Internal multiples from above the redatuming level are eliminated during this process. In the past, we introduced point‐spread functions for interferometric redatuming by cross‐correlation. These point‐spread functions quantify distortions in the redatumed data, caused by internal multiple reflections in the overburden. In this paper, we define point‐spread functions for interferometric imaging to quantify these distortions in the image domain. These point‐spread functions are similar to conventional resolution functions for seismic migration but they contain additional information on the internal multiples in the overburden and they are partly data‐driven. We show how these point‐spread functions can be visualized to diagnose image defocusing and artefacts. Finally, we illustrate how point‐spread functions can also be defined for interferometric imaging with passive noise sources in the subsurface or with simultaneous‐source acquisition at the surface.     

A framework for the evaluation of ground motion selection and modification procedures

This study develops a framework to evaluate ground motion selection and modification (GMSM) procedures. The context is probabilistic seismic demand analysis, where response history analyses of a given structure, using ground motions determined by a GMSM procedure, are performed in order to estimate the seismic demand hazard curve (SDHC) for the structure at a given site. Currently, a GMSM procedure is evaluated in this context by comparing several resulting estimates of the SDHC, each derived from a different definition of the conditioning intensity measure (IM). Using a simple case study, we demonstrate that conclusions from such an approach are not always definitive; therefore, an alternative approach is desirable. In the alternative proposed herein, all estimates of the SDHC from GMSM procedures are compared against a benchmark SDHC, under a common set of ground motion information. This benchmark SDHC is determined by incorporating a prediction model for the seismic demand into the probabilistic seismic hazard analysis calculations. To develop an understanding of why one GMSM procedure may provide more accurate estimates of the SDHC than another procedure, we identify the role of ‘IM sufficiency’ in the relationship between (i) bias in the SDHC estimate and (ii) ‘hazard consistency’ of the corresponding ground motions obtained from a GMSM procedure. Finally, we provide examples of how misleading conclusions may potentially be obtained from erroneous implementations of the proposed framework. Copyright © 2014 John Wiley & Sons, Ltd.     

黏弹性介质中瑞利波频散曲线和衰减系数曲线的反演

瑞利波具有能量大、信噪比高等特点,可以用来反演介质内部的力学信息,近年来在浅层地球物理勘探、深层地震学研究以及超声波无损检测等多个领域都有较广泛的应用。目前大多数瑞利波的应用中都假设介质是弹性的,然而实际中岩石、土壤和金属等介质都在一定程度上体现出了黏弹性。当介质的黏弹性较强时仍然采用弹性假设研究其中瑞利波的反演将增大误差,因此有必要考虑黏弹性介质中的瑞利波反演,但是目前这方面的研究仍不够深入。本文研究黏弹性介质中瑞利波频散曲线和衰减系数曲线的反演问题,给出其在半空间中联合反演的方法,并对该方法的误差进行分析。     

2014年于田M_S7.3地震序列的频谱特征分析及其前震识别

利用新疆区域数字地震台网的波形资料,对2014年2月12日于田MS7.3地震前震和余震序列的频谱特征进行了分析。基于Brune模型,对观测记录进行仪器响应、传播路径和场地响应的校正,使用遗传算法对于田地震序列103次ML≥3.0地震的震源谱进行了反演计算。结果表明:1前震序列的拐角频率相对偏低,扣除地震矩影响后,余震的拐角频率有明显增高的趋势;2距离震中最近的YUT台资料分析显示,于田地震几次前震活动的台站观测谱相对较为独立,与主震以及彼此之间的相关系数均不高;3序列事件和主震的台站观测谱相关系数的变化反映了序列发展的动态破裂过程,相关系数较低的事件基本上偏离了主震破裂区而发生在分支断裂上;4前震的震源谱中存在异常的高频成份。