高铁震源地震信号的挤压时频分析应用

我国每天有数千趟高铁列车驰骋在纵横交错的高铁线路上,构成了十分理想的均布震源,但寻找适合高铁震源地震信号的处理方法是充分挖掘信息的关键.传统的频谱分析结果表明高铁震源所产生的地震信号具有明显的窄带分立谱特征,但无法精确获得高铁震源地震信号的时频变化规律.本文首次将挤压时频分析这种分析工具引入到高铁震源地震信号处理中,对中国南方某高铁沿线采集到的高铁震源地震数据进行了分析.处理结果表明:利用挤压时频分析能够更加精确地刻画频率成分随时间的变化,能够利用单检波器精确刻画高铁列车的运行状态(匀速、加速等);同时利用挤压时频变换还可高精度地重构出所需频带的信号,为提取高铁震源地震信号的特征成分提供了一种有力工具.     

HVSR方法用于地震作用下场地效应分析的适用性研究

以美国Garner Valley Downhole Array(GVDA)竖向观测台阵的强震观测记录为基础,利用水平与竖向谱比(Horizontal-to-Vertical Spectral Ratio,HVSR)方法与传统谱比方法研究了多次地震作用下竖向台阵场地的波动传播规律,探讨了该场地HVSR曲线与传递函数(Transfer Function,TF)曲线的差异,综合利用频域求解一维层状场地P、SV波入射情况下场地响应的土层地震反应分析方法和合成理论地震图的波数积分法解释了二者差异可能的原因,在此基础上研究了HVSR方法应用于地震作用下场地效应分析的适用性.研究结论表明:竖向地震动场地效应是HVSR与传递函数产生差异的主要原因,竖向地震动场地效应主要是由竖向地震动中的P波成分引起的.在场地竖向放大可以忽略的频段范围内,HVSR可以作为传递函数研究地震作用下的场地效应.     

Illumination effects on the differenced Normalized Burn Ratio's optimality for assessing fire severity

The influence of illumination effects on the optimality of the dNBR (differenced Normalized Burn Ratio) was evaluated for the case of the 2007 Peloponnese (Greece) wildfires using a pre/post-fire Landsat TM (Thematic Mapper) image couple. Well-illuminated pixels (south and south-east facing slopes) exhibited more optimal displacements in the bi-spectral feature space than more shaded pixels (north and north-west exposed slopes). Moreover, pixels experiencing a small image-to-image difference in illumination obtained a higher optimality than pixels with a relatively large difference in illumination. To correct for illumination effects, the c-correction method and a modified c-correction technique were applied. The resulting median dNBR optimality of uncorrected, c-corrected and modified c-correction data was respectively 0.58, 0.60 and 0.71 (differences significant for p < 0.001). The original c-correction method improved the optimality of badly illuminated pixels while deteriorating the optimality of well-illuminated pixels. In contrast, the modified c-correction technique improved the optimality of all the pixels while retaining the prime characteristic of topographic correction techniques, i.e. detrending the illumination–reflectance relationship. For a minority of the data, for shaded pixels and/or pixels with a high image-to-image difference in illumination, the original c-correction outperformed the modified c-correction technique. In this study conducted in rugged terrain and with a bi-temporal image acquisition scheme that deviated up to two months from the ideal anniversary date scheme the modified c-correction technique resulted in a more reliable change detection.     

淮河流域洪涝灾害时间对称性研究

准确判断重大自然灾害发生趋势对于防灾减灾意义重大.基于时间对称性理论,本文构建“可公度降噪—对称性趋势判断—信息结构系稳定性检验”研究模式,判断了淮河流域(HRB)洪涝灾害发生趋势.结果表明:(1)通过可公度降噪处理去除了噪声元素,提高了预测的精确性;(2)运用三元可公度、蝴蝶结构图、信息结构系等方法判断,2019年HRB发生洪涝灾害的信号较强,该年发生洪涝灾害的随机性概率为66.7%,不漏报置信水平为57.1%;(3)通过洪涝灾害与太阳黑子活动的相关性分析可知,HRB在太阳黑子活动谷年附近易发生洪涝灾害,进一步提高了预测结果的可信度.“可公度降噪—对称性趋势判断—信息结构系稳定性检验”研究模式是对现有时间对称性方法的补充和完善,以期为中国灾害预测起到一定的推动作用.     

A high accurate automated first-break picking method for seismic records from high-density acquisition in areas with a complex surface

As the application of high-density high-efficiency acquisition technology becomes more and more wide, the areas with complex surface conditions gradually become target exploration areas, and the first-break picking work of massive low signal-to-noise ratio data is a big challenge. The traditional method spends a lot of manpower and time to interactively pick first breaks, a large amount of interactive work affects the accuracy and efficiency of picking. In order to overcome the shortcoming that traditional methods have weak anti-noise to low signal-to-noise ratio primary wave, this paper proposes a high accurate automated first-break picking method for low signal-to-noise ratio primary wave from high-density acquisition in areas with a complex surface. Firstly, this method determines first-break time window using multi-azimuth spatial interpolation technology; then it uses the improved clustering algorithm to initially pick first breaks and then perform multi-angle comprehensive quality evaluation to first breaks according to the following sequence: ‘single trace → spread → single shot → multiple shots’ to identify the abnormal first breaks; finally it determines the optimal path through the constructed evaluation function and using the ant colony algorithm to correct abnormal first breaks. Multi-azimuth time window spatial interpolation technology provides the base for accurately picking first-break time; the clustering algorithm can effectively improve the picking accuracy rate of low signal-to-noise ratio primary waves; the multi-angle comprehensive quality evaluation can accurately and effectively eliminate abnormal first breaks; the ant colony algorithm can effectively improve the correction quality of low signal-to-noise ratio abnormal first breaks. By example analysis and comparing with the commonly used Akaike Information Criterion method, the automated first-break picking theory and technology studied in this paper has high picking accuracy and the ability to stably process low signal-to-noise ratio seismic data, has a significant effect on seismic records from high-density acquisition in areas with a complex surface and can meet the requirements of accuracy and efficiency for massive data near-surface modelling and statics calculation.     

Prestack data enhancement with phase corrections in time-frequency domain guided by local multidimensional stacking

We present a new approach to enhancing weak prestack reflection signals without sacrificing higher frequencies. As a first step, we employ known multidimensional local stacking to obtain an approximate ‘model of the signal’. Guided by phase spectra from this model, we can detect very weak signals and make them visible and coherent by ‘repairing’ corrupted phase of original data. Both presented approaches – phase substitution and phase sign corrections – show good performance on complex synthetic and field data suffering from severe near-surface scattering where conventional processing methods are rendered ineffective. The methods are mathematically formulated as a special case of time-frequency masking (common in speech processing) combined with the signal model from local stacking. This powerful combination opens the avenue for a completely new family of approaches for multi-channel seismic processing that can address seismic processing of land data with nodes and single sensors in the desert environment.     

A high-precision time–frequency analysis for thin hydrocarbon reservoir identification based on synchroextracting generalized S-transform

Improving the seismic time–frequency resolution is a crucial step for identifying thin reservoirs. In this paper, we propose a new high-precision time–frequency analysis algorithm, synchroextracting generalized S-transform, which exhibits superior performance at characterizing reservoirs and detecting hydrocarbons. This method first calculates time–frequency spectra using generalized S-transform; then, it squeezes all but the most smeared time–frequency coefficients into the instantaneous frequency trajectory and finally obtains highly accurate and energy-concentrated time–frequency spectra. We precisely deduce the mathematical formula of the synchroextracting generalized S-transform. Synthetic signal examples testify that this method can correctly decompose a signal and provide a better time–frequency representation. The results of a synthetic seismic signal and real seismic data demonstrate that this method can identify some reservoirs with thincknesses smaller than a quarter wavelength and can be successfully applied for hydrocarbon detection. In addition, examples of synthetic signals with different levels of Gaussian white noise show that this method can achieve better results under noisy conditions. Hence, the synchroextracting generalized S-transform has great application prospects and merits in seismic signal processing and interpretation.     

2019年6月17日长宁MS6.0地震前地震活动和水位的加卸载响应比异常研究

加卸载响应比（LURR）方法是通过固体潮加卸载过程中的某一物理参数的响应差异探查区域应力场演化。本文将这一方法应用于2019年6月17日四川长宁M_S6.0地震,根据LURR时序演化探查震源区介质的应力状态变化,提取可能的震前异常信息。首先采用贝尼奥夫应变作为响应量计算LURR,通过固体潮在地震断层面的最优滑动方向上引起的库仑破坏应力变化来判断加载还是卸载,结果显示LURR值自2018年年初开始快速增加并在2018年年中达到峰值,此后异常持续至地震。在此基础上对该地区的大地形变和地下井水位资料进行了LURR分析,发现水富水平摆倾斜仪的东西和南北分量在LURR出现异常的同期发生明显偏转,而地震周边的地下井水位高值变化过程与长宁地震也有较好的对应关系。研究表明在长宁地震发生前,存在地震活动和地下水位的LURR异常,且异常时间与水富倾斜同步,暗示震源区介质存在明显的应力积累过程。     

Three-parameter Radon transform in layered transversely isotropic media

Radon transform is a powerful tool with many applications in different stages of seismic data processing, because of its capability to focus seismic events in the transform domain. Three-parameter Radon transform can optimally focus and separate different seismic events, if its basis functions accurately match the events. In anisotropic media, the conventional hyperbolic or shifted hyperbolic basis functions lose their accuracy and cannot preserve data fidelity, especially at large offsets. To address this issue, we propose an accurate traveltime approximation for transversely isotropic media with vertical symmetry axis, and derive two versions of Radon basis functions, time-variant and time-invariant. A time-variant basis function can be used in time domain Radon transform algorithms while a time-invariant version can be used in, generally more efficient, frequency domain algorithms. Comparing the time-variant and time-invariant Radon transform by the proposed basis functions, the time-invariant version can better focus different seismic events; it is also more accurate, especially in presence of vertical heterogeneity. However, the proposed time-invariant basis functions are suitable for a specific type of layered anisotropic media, known as factorized media. We test the proposed methods and illustrate successful applications of them for trace interpolation and coherent noise attenuation.     

Seismic receiver coupling to the seafloor

The presence of geophysical receivers on the seafloor changes the local wave field due to the receiver seafloor interaction. The resulting PP- and PS-wave distortion of the wave field is often referred to as receiver coupling to the seafloor and can make data processing challenging and sometimes impossible. This paper provides an overview of the mathematical approaches to describe receiver coupling, how to estimate receiver coupling and what the difficulties and possible solutions are. The first section shows how the mathematical approach developed from a simple model considering only the vertical receiver component to include all three receiver components and their interactions with the seafloor. In the second section, I show how receiver coupling can be measured and how it can be improved using mathematical and data-driven approaches.