Wavelet correlation filter for wide-angle seismic data

A new filtering technique for single‐fold wide‐angle reflection/refraction seismic data is presented. The technique is based on the wavelet decomposition of a set of adjacent traces followed by coherence analysis. The filtering procedure consists of three steps. In the first, a wavelet decomposition of traces into different detail levels is performed. In the second, the coherence attributes for each level are evaluated by calculating cross‐correlation functions of detail portions contained in a space–time moving window. Finally, the filtered traces are obtained as a weighted reconstruction of the trace details. Each weight is obtained from the coherence‐attributes distribution estimated in a proper interval. A sequence of tests is then conducted in order to select possible optimum or unsuitable wavelet bases. The efficiency of the filter proposed was assessed by calculating some properly designed parameters in order to compare it with other standard de‐noising techniques. The proposed method produced a clear signal enhancement in high‐density wide‐angle seismic data, thus proving that it is a useful processing tool for a reliable correlation of seismic phases.     

模板匹配滤波技术在地震数据处理中的应用

模板匹配滤波技术(Matched Filter Technique,MFT)是一种近几年发展迅速的可用于检测遗漏信号的算法,在图像识别、行星探测等众多领域均有广泛应用。在地震学中,该技术通过扫描连续波形找出模板事件附近发生的、与模板事件相关性较好的遗漏地震事件,可用于完善地震目录,实现对地震活动性以及区域应力状态的更好刻画。本文综述了模板匹配滤波技术在地震数据处理中的应用,首先概述基本原理,其次介绍该方法用于探测各类地震事件的实例,包括前震和余震、远程动态触发地震、诱发地震和非火山震颤等,最后展望了用于提高算法效率及探测精度的一些改进方法。     

Magnetotelluric data inversion with seismic data constraint

Introduction The research on the structure and physical property of ancient hidden hill, igneous rocks and basement is relatively difficult by using seismic data only. If we combine seismic data, magneto-telluric (MT) data and geophysical data together, better results can be obtained for the above problem. A number of geophysicists at home and abroad, such as CHEN and WANG (1990), Siri-punvarapor and Egbert (2000) have tried many methods to solve the problem by the inversion of seismic da…     

模糊嵌套多级中值滤波方法及其在地震数据处理中的应用

多级中值滤波器在地震数据处理应用中,其滤波长度越大,消噪效果越好,但同时也会破坏有效信息. 如果为了保护有效信息而减少滤波长度,又会造成大量的噪声不能消除. 本文提出一种新的模糊嵌套多级中值滤波器,设计一个阀值作为判断参数,使滤波器能够在消除随机噪声时采用长滤波器滤波,而在保留有效信息时采用短滤波器滤波,从而既能很好地消除随机噪声,又能最大限度地保护有效信息,保留有效信息的细节结构. 经过模型分析和实际资料处理都取得了很好的效果.     

Vertical seismic behaviour of suspension bridges

The vertical response of suspension bridges to multi-support seismic excitations is investigated. A frequency-domain random-vibration approach is utilized to take into account not only the differences in ground motion inputs, but also the correlation among the various input motions. The earthquake response of the Golden Gate Suspension Bridge, California, is analysed. The ground motion inputs were taken from an array of time histories recovered from the Imperial Valley, California, earthquake (M_L = 6.6) of 15 October 1979. Mean square displacement, stress and cable tension responses are obtained. An estimation of the expected peak response is made by comparing existing methods to the results of a time-domain analysis; it was found that a peak factor of 3.5 times the r.m.s. response is a good estimate of the expected peak response for this multiple-input seismic problem.     

A seismic interpolation and denoising method with curvelet transform matching filter

A new seismic interpolation and denoising method with a curvelet transform matching filter, employing the fast iterative shrinkage thresholding algorithm (FISTA), is proposed. The approach treats the matching filter, seismic interpolation, and denoising all as the same inverse problem using an inversion iteration algorithm. The curvelet transform has a high sparseness and is useful for separating signal from noise, meaning that it can accurately solve the matching problem using FISTA. When applying the new method to a synthetic noisy data sets and a data sets with missing traces, the optimum matching result is obtained, noise is greatly suppressed, missing seismic data are filled by interpolation, and the waveform is highly consistent. We then verified the method by applying it to real data, yielding satisfactory results. The results show that the method can reconstruct missing traces in the case of low SNR (signal-to-noise ratio). The above three problems can be simultaneously solved via FISTA algorithm, and it will not only increase the processing efficiency but also improve SNR of the seismic data.     

灌芯混凝土石膏墙板结构抗震性能试验研究

本文通过对2块灌芯混凝土纤维增强石膏板与8个节点构件在低周反复荷载下的抗震性能试验,研究了构件的抗剪承载能力、构件的延性、恢复力特性、耗能能力和破坏特征等特性;通过实测试点工程结构的动力反应,获得结构的自振频率、振型和阻尼比等模态参数。研究表明,该结构在抗震方面具有良好的性能。     

L1范数约束被动源数据稀疏反演一次波估计

对于被动源地震数据,运用常规的互相关算法得到的虚拟炮记录中,不仅含有一次波反射信息,还包括了表面相关多次波.然而,通过传统的被动源数据稀疏反演一次波估计(EPSI)方法,可以求得只含有一次波,不含表面相关多次波的虚拟炮记录.本文改进了传统的被动源数据稀疏反演一次波估计问题的求解方法,将被动源稀疏反演一次波估计求解问题转化为双凸L1范数约束的最优化求解问题,避免了在传统的稀疏反演一次波估计过程中用时窗防止反演陷入局部最优化的情况.在L1范数约束最优化的求解过程中,又结合了2D Curvelet变换和小波变换,在2D Curvelet-wavelet域中,数据变得更加稀疏,从而使求得的结果更加准确,成像质量得到了改善.通过简单模型和复杂模型,验证了本文提出方法的有效性.     

Computing and plot seismic data with electronic computer

Summary The paper gives a simple method for computing of seismic data with electronic computer. The method described is suitable not only in the case of linear ray path but also for curved ray path. There is also a seismic profile included, plotted by digital-computer according to the explained process.

---

Zusammenfassung Die Arbeit gibt eine einfache Methode für Berechnung der seismischen Daten mit elektronischer Rechenmaschine. Die gegebene Methode ist nicht nur für geradlinigen Strahlenweg, sondern auch für gekrümmten Strahlenweg anwendbar. Es ist auch ein mit elektronischer Maschine berechnetes und gedrucktes Profil beiliegend.

     

地震资料小波剖面制作系统

对不同的小波变换算法和小波基函数进行了研究，得到了适用于地震资料处理的最佳小波基函数，并用Visual C 语言开发了基于Windows操作平台下的地震资料小波剖面制作系统。利用小波剖面制作系统和最佳小波基函数对煤田实际地震资料进行了处理，取得了令人满意的地质效果。     

非稳态地震记录时变子波估计

地震子波估计是地震资料处理和解释中的一个关键问题,子波估计的可靠性会直接影响反褶积和反演的准确度.现有的子波估计方法分为确定型和统计型两种类型,本文通过结合这两类方法,利用确定型的谱分析法和统计型的偏度最大化方法,分别提取时变子波的振幅和相位信息,得到估计的时变子波.这种方法不需要对子波进行任何时不变或相位等的假设,具有对时变相位的估计能力.进而利用估计时变子波进行非稳态反褶积,提高地震记录的保真度,为精细储层预测和描述提供高质量的剖面.理论模型试算验证了方法的可行性,通过实际地震资料的处理应用,表明该方法能有效地提取出子波时变信息.

     

基于S变换的软阈值滤波在深地震反射数据处理中的应用

深地震反射原始单炮数据是非平稳的弱能量反射信号, 信噪比较低. 如何提高信噪比一直是深地震反射数据前处理中的一大难题. S变换是一种适用于分析非平稳信号的时频变换方法. 同其他分析时变信号的方法相比, S变换的基本小波不必满足小波在时间域均值为零的容许性条件, 它的时频分辨率与分析信号的频率有关, 且其在时间域的积分可以得到傅里叶频谱,其反变换也简单. 因此, S变换容易表示深地震反射信号复杂的时频特性. 本文在S变换的基础上, 利用软阈值滤波方法对深地震反射数据进行处理, 实验结果表明, 该方法有效地提高了信噪比, 压制了有效频带范围内的混频干扰, 突出了弱反射信号, 使得波组信息更加丰富, 有利于连续追踪有效反射波组和识别薄地层, 特别是提高了深部Moho界面反射层位的分辨率, 为深地震反射剖面后续处理和准确解释奠定了基础.     

Time-domain identification of dynamic properties of layered soil by using extended Kalman filter and recorded seismic data

A novel time-domain identification technique is developed for the seismic response analysis of soil-structure interaction. A two-degree-of-freedom (2DOF) model with eight lumped parameters is adopted to model the frequency-dependent behavior of soils. For layered soil, the equivalent eight parameters of the 2DOF model arc identified by the extended Kalman filter (EKF) method using recorded seismic data. The polynomial approximations for derivation of state estimators are applied in the EKF procedure. A realistic identification example is given for the layered-soil of a building site in Anchorage, Alaska in the United States. Results of the example demonstrate the feasibility and practicality of the proposed identification technique. The 2DOF soil model and the identification technique can be used for nonlinear response analysis of soil-structure interaction in the time-domain for layered of complex soil conditions. The identified parameters can be stored in a database for use in other similar soil conditions. If a universal database that covers information related to most soil conditions is developed in the future, engineers could conveniently perform time history analyses of soil-structural interaction.     

The effects of source and receiver motion on seismic data1

The effects of source and receiver motion on seismic data are considered using extensions of the standard convolutional model. In particular, receiver motion introduces a time-variant spatial shift into data, while source motion converts the effect of the source signature from a single-channel convolution in time to a multichannel convolution in time and space. These results are consistent with classical Doppler theory and suggest that Doppler shifting can introduce distortions into seismic data even at relatively slow acquisition speeds. It is shown that, while both source and receiver motion are known to be important for marine vibroseis acquisition, receiver motion alone can produce significant artifacts in marine 3D data. Fortunately, the convolutional nature of the distortions renders them amenable to correction using simple deconvolution techniques. Specifically, the effects of receiver motion can be neutralized by applying an appropriate reverse time-variant spatial shift, while those due to source motion can be addressed by introducing time-variant spatial shifts both before and after standard, deterministic, signature deconvolution or correlation.     

基于结构自适应中值滤波器的随机噪声衰减方法

本文提出一种保护断层、裂缝等地层边缘特征的结构自适应中值滤波器,用于衰减地震资料中的随机噪声.基于地震反射同相轴局部呈线型结构的假设,采用梯度结构张量估计地层倾向,分析地层结构的规则程度,在此基础上引入地震剖面中线型和横向不连续性两种结构特征的置信度量.结构自适应中值滤波器根据这两种置信度量调整滤波器窗函数的尺度和形状,根据地层倾角调整滤波器窗函数的方向,从而使得滤波操作窗能够最佳匹配信号的局部结构特征.将本文方法用于合成和实际数据的处理,并与两种常用中值滤波方法进行对比,结果表明,该方法能够更好地解决地震剖面的随机噪声衰减和有效信号保真的问题,在增强反射同相轴的横向一致性的同时有效保持了剖面内的地层边缘和细节特征,显著改善了地震资料的品质.     

用稳定高效的反Q滤波技术提高地震资料分辨率

地震波在地下传播时受到衰减影响,衰减会导致地震波场高频能量的损失和相位畸变.反Q滤波可补偿大地衰减效应.已有的反Q滤波方法存在下列不足:频率域的算法由于算子长度较长,所以计算效率较低;时间域的算法,或者对地震记录上到时较晚的同相轴进行了过度的补偿,或者为防止过度补偿后来的振幅而在最大增益处进行限制,导致振幅的多解性,而且还会影响滤波器的相位效应.本文给出一种通过直接求解时间域的Q模型方程来进行反Q滤波的算法.由于采用带状矩阵解算器,所以具有较高的计算效率,理论数据和实际地震资料的试算结果证明,本方法对地震波的吸收衰减进行了出色的补偿,提高了地震资料的分辨率.     

The seismic behaviour of cross-braced steel frames

Analytical studies on the inelastic behaviour of concentrically braced steel frames for low-rise buildings are described in this paper. The bracing members which provide energy dissipation were used to provide information on the ductility levels that are likely to occur under differing levels of earthquake excitation. An indication of the relative performance of cross bracing is provided in terms of suitable SM values for use in the seismic provisions of New Zealand loadings code NZS 4203.     

地震资料多重分形裂缝检测技术方法

在讨论了多重分形理论后,结合地震多重分形的存在性、可行性等问题的分析,针对灰岩裂缝带地震记录特征,研究了适于研究区灰岩裂缝地震多重分形检测方法.在所研究方法中,针对以往地震记录离差、求和、极值等分形测度存在的问题,设计了能够揭示地震变化的全信息分形测度刻划方法.同时对窗口函数的选取、配分函数计算中的幂次的选择及拟合计算问题进行了讨论.在理论模型验证方法正确合理性的基础上,研究了实际资料的应用.应用结果表明:1)多重分形计算结果是保幅的,并且能够揭示地震振幅异常信息;2)对较小断层甚至是裂缝带在相干体等一般技术方法不能识别这一问题上,得到了较好地解决,并通过多口井资料证实了所解决问题的有效性和效果.     

在SUSE Linux系统下通过ADSL网络实现地震波形浏览

针对SUSELinux系统下非固定IP接入的ADSL入网方式不能直接浏览地震波形的问题,介绍了解决这一问题的方法,该方法已在沈阳地震台和牡丹江地震台实践取得了成功。     

基于Dreamlet变换的地震数据压缩理论与方法

为达到更加有效地表示地震数据的目的,仅仅将地震数据当作普通的图像数据处理是远远不够的,地震数据中蕴含的地震波的运动学特性也应作为重要因素而被考虑到.本文讨论了利用Dreamlet变换方法实现地震数据压缩的方法,并针对地震数据本身所蕴含的频散关系特性进一步提出了多尺度Dreamlet变换压缩方法.Dreamlet变换由2个一维局部谐波变换的张量积构成,它在提供地震波场时间-空间局部化性质的同时可以保留波场的运动学特性.通过对二维SEG/EAGE叠前、叠后数据的算例说明了Dreamlet变换用于地震数据压缩的有效性.利用压缩后的数据进行成像的结果更表明,与Curvelet变换方法相比,Dreamlet与多尺度Dreamlet方法可以提供更高的压缩比;在相同压缩比的条件下,使用Dreamlet与多尺度Dreamlet方法压缩重建后的数据进行成像能更好地保留成像结果中的重要结构.