Simultaneous Denoising and Interpolation of Seismic Data via the Deep Learning Method

Utilizing data from controlled seismic sources to image the subsurface structures and invert the physical properties of the subsurface media is a major effort in exploration geophysics. Dense seismic records with high signal-to-noise ratio (SNR) and high fidelity helps in producing high quality imaging results. Therefore, seismic data denoising and missing traces reconstruction are significant for seismic data processing. Traditional denoising and interpolation methods rarely occasioned rely on noise level estimations, thus requiring heavy manual work to deal with records and the selection of optimal parameters. We propose a simultaneous denoising and interpolation method based on deep learning. For noisy records with missing traces, we adopt an iterative alternating optimization strategy and separate the objective function of the data restoring problem into two sub-problems. The seismic records can be reconstructed by solving a least-square problem and applying a set of pre-trained denoising models alternatively and iteratively.We demonstrate this method with synthetic and field data.     

天山东北部地震的重新定位和一维地壳速度模型的改善

We apply three methods to relocate 599 earthquake events that occurred from August 2004 to August 2005 in the northeastern Tianshan Mountains area ( 85°30’ ～ 88°30’E,43°00’ ～ 44°40’ N ) by using travel times recorded by regional seismic network and 10 portable seismic stations deployed around the Urumqi city. By comparing the reliability of different results,we determined a suitable location method,and an improved 1-D crustal velocity model of the study area. The uncertainty of earthquake location is significantly reduced with combined data of seismic network and portable stations. The relocated events are clearly associated with regional tectonics of the northeastern Tianshan Mountains area, and are also in agreement with the existence of active faults imaged by deep seismic reflection profile. The relocated seismicity discovers some potential traces of buried active faults,which need to be validated further.     

Suppressing random noise to broaden the useful bandwidth of seismic reflection data

The signal-to-noise ratio （SNR） of seismic reflection data in many areas is rather poor and conventional two-dimensional filters designed to suppress noise with different moveout from the signal tend to generate artifacts. We have extended a method of multichannel filtering, based on the hypothesis that signals on adjacent channels are similar, for enhancing the SNR on stacked sections. Using only the mid-range frequencies where the SNR is highest, the event trend is found for overlapping windows on the section and the average signal vector is calculated. Then the data from the full bandwidth section are projected onto the spatially varying unit similarity vectors and the results are merged for the overlapping windows. Application of the method to synthetic data containing steeply dipping events and to a stacked section for a marine 2D line has produced good results. The modifications we have introduced carry a small overhead in computing time but they should enable the method to be used effectively even on sections containing steep dips.     

Seismic inversion for underground fractures detection based on effective anisotropy and fluid substitution

Underground fractures play an important role in the storage and movement of hydrocarbon fluid. Fracture rock physics has been the useful bridge between fracture parameters and seismic response. In this paper, we aim to use seismic data to predict subsurface fractures based on rock physics. We begin with the construction of fracture rock physics model. Using the model, we may estimate P-wave velocity, S-wave velocity and fracture rock physics parameters. Then we derive a new approximate formula for the analysis of the relationship between fracture rock physics parameters and seismic response, and we also propose the method which uses seismic data to invert the elastic and rock physics parameters of fractured rock. We end with the method verification, which includes using well-logging data to confirm the reliability of fracture rock physics effective model and utilizing real seismic data to validate the applicability of the inversion method. Tests show that the fracture rock physics effective model may be used to estimate velocities and fracture rock physics parameters reliably, and the inversion method is resultful even when the seismic data is added with random noise. Real data test also indicates the inversion method can be applied into the estimation of the elastic and fracture weaknesses parameters in the target area.     

Estimating instantaneous concentration of suspended sediment using acoustic backscatter from an ADV

The measurement of instantaneous sediment concentration remains a challenging task. In this study, a three-step procedure is proposed to estimate instantaneous sediment concentration using acoustic backscatter from Acoustic Doppler Velocimeter (ADV). The influences of acoustic noise and particle diameter on an ADV's performance was first tested in the laboratory, then the three-step procedure was verified based on field measurements in the Zhongxian and Fengjie reaches in the Three Gorges Reservoir. The first step involves reconstructing the backscatter signal time series. Due to contamination from both the noise floor and spurious spikes, the denoising-despiking method was applied instead of the traditional velocity-despiking methods, and this approach performed well based on spectrum ana-lysis. The second step involves calibrating the sediment concentration against the backscatter signal. A linear relation, whose slope and intercept were calibrated to be dependent on particle diameter, is proposed in double logarithmic coordinates. The third step involves calculating the instantaneous sediment concentration using reconstructed instantaneous backscatter based on the proposed relation. The accuracy of the proposed method was evaluated through consideration of the concentration spec-trum and sediment flux, indicating that the proposed three-step procedure is effective for the mea-surement of instantaneous sediment concentration.     

对台湾强震的预警时间研究

The paper collects the records by the Fujian Digital Seismic Network of 40 shallow earthquakes in Taiwan with M_S≥5.0 from 1999 to 2013,analyzes the seismic phase(Pn,Sn phase)characteristics and travel-time rules,determines travel-time models and develops a seismic phase travel-time equation based on the two-step fitting algorithm.With the deduction of processing time and network delay time,this method can provide an accurate estimation of early warning time of Taiwan earthquakes for the Fujian region,and has been officially employed in the earthquake early warning system of Fujian Province.     

Serviceability evaluation of water supply networks under seismic loads utilizing their operational physical mechanism

The serviceability of water supply networks(WSNs)under seismic loads has significant importance for estimating the probable losses and the impact of diminished functionality on affected communities.The innovation presented in this paper is suggesting a new strategy to evaluate the seismic serviceability of WSNs,utilizing their operational physical mechanism.On one hand,this method can obtain the seismic serviceability of each node as well as entire WSNs.On the other hand,this method can dynamically reflect the propagation of randomness from ground motions to WSNs.First,a finite element model is established to capture the seismic response of buried pipe networks,and a leakage model is suggested to obtain the leakage area of WSNs.Second,the transient flow analysis of WSNs with or without leakage is derived to obtain dynamic water flow and pressure.Third,the seismic serviceability of WSNs is analyzed based on the probability density evolution method(PDEM).Finally,the seismic serviceability of a real WSN in Mianzhu city is assessed to illustrate the method.The case study shows that randomness from the ground motions can obviously affect the leakage state and the probability density of the nodal head during earthquakes.     

面状发震构造在地震构造法中的应用：以大姚—姚安发震构造鉴定为例

In seismic hazard analysis of nuclear power plant of China there is a need to identify both seismogenic structures and seismotectonic zones. In past practice, the identification of the seismogenic structures was often based on the surface active faults and characterization of linear seismic source. In a situation which shows quite strong non-random seismic activity and lacks surface active faults, it is difficult to evaluate the seismic hazard reasonably. Taking seismogenic structures in the Dayao-Yao＇an area as a case study in this paper, we discuss the need and the possibility to apply the planar seismogenic structure to the seismotectonic method. We suggest that the planar seismogenic structure should be considered when applying the seismotectonic method to the seismic risk assessment of nuclear engineering in future.     

Comparative study on the deviation characteristics of post-earthquake intensity evaluation based on two intensity attenuation relationship models

We collect the isoseismal data of 488 earthquakes in western China and 182 earthquakes in eastern China after 1900. The least square method is used to establish two models of the seismic intensity attenuation relationship partitions in China, based on the major and minor axis ellipse models, and the area and the major axis radius,respectively. The two models are applied to the calculation of the intensity circle of the earthquake events with a magnitude M_S≥5.0 from 2008 to 2019, and the actual intensity is compared with the model intensity value as an indicator to verify the consistency between the actual intensity and the value from the empirical statistical model.Three real earthquake results are selected to calculate the major and minor axis and area of the ellipse using the two intensity attenuation relationship models. After comparison, we summarize the deviation characteristics of the intensity value, and put forward corresponding improvement suggestions.     

Weak Seismic Signal Extraction Based on the Curvelet Transform

Seismic signal denoising is a key step in seismic data processing. Airgun signals are easy to be interfered with by noise when it travels a long distance due to the weak energy of active source signal of the airgun. Aiming to solve this problem, and considering that the conventional Curvelet transform threshold processing method does not use the seismic spectrum information, we independently process the Curvelet scale layer corresponding to valid data based on the characteristics of the Curvelet transform of multi-scale, multi-direction and capable of expressing the sparse seismic signals in order to fully excavate the information features. Combined with the Curvelet adaptive threshold denoising the algorithm, we apply the Curvelet transform to denoising seismic signals while retaining the weak information in the signal as much as possible. The simulation experiments show that the improved threshold denoising method based on Curvelet transform is superior to the frequency domain filtering, wavelet denoising and traditional Curvelet denoising method in detailed information extraction and signal denoising of low SNR signals. The calculation accuracy of the relative wave velocity variation of underground medium is improved.     

自适应多道匹配追踪去强屏蔽方法研究与应用

阻抗差异较大的地层在地震剖面上呈现为强振幅同相轴,会掩盖附近储层有效信息,需要做去除强屏蔽的目标处理。针对常规匹配追踪方法在构造复杂地区匹配精度和空间连续性不佳的问题,本文提出一种基于自适应权值的多道匹配追踪去强屏蔽方法。首先利用层位构造信息对强反射层进行局部拉平处理,减弱地层空间构造对提取强反射层的影响,然后引入相邻地震道与中心道之间的相关系数作为多道平均的权值,以提高匹配结果的稳定性和横向连续性。同时使用解释得到的层位时间作为子波初始时间,有效提高运算效率。模型试算和实际地震资料应用分析表明,改进方法能有效剥离强反射并突出有效储层信息,剥离后的井震吻合度得到明显提升,相较于常规匹配追踪算法匹配精度更高,同时具有更好的空间连续性与强反射去除效果。      

Well tie,fluid substitution and AVO modelling: a North Sea example

Accurate well ties are essential to practical seismic lithological interpretation. As long as the geology in the vicinity of the reservoir is not unduly complex, the main factors controlling this accuracy are the processing of the seismic data and the construction of the seismic model from well logs. This case study illustrates how seismic data processing to a near-offset stack, quality control of logs and petrophysical modelling improved a well tie at an oil reservoir. We demonstrate the application of a predictive petrophysical model in the preparation and integration of the logs before building the seismic model and we quantify our improvements in well-tie accuracy. The data for the study consisted of seismic field data from a 3D sail line through a well in a North Sea oilfield and a suite of standard logs at the well. A swathe of fully processed 3D data through the well was available for comparison. The well tie in the shallow section from first-pass seismic data processing and a routinely edited sonic log was excellent. The tie in a deeper interval containing the reservoir was less satisfactory: the phase errors within the bandwidth of the seismic wavelet were of the order of 20°, which we consider too large for subsequent transformation of the data to seismic impedance. Reprocessing the seismic data and revision of the well-log model reduced these phase errors to less than 10° and improved the consistency of the deep and shallow well ties. The reprocessing included densely picked iterative velocity analysis, prestack migration, beam-forming multiple attenuation, stacking the near-offset traces and demigration and remigration of the near-offset data. The petrophysical model was used to monitor and, where necessary, replace the P-wave sonic log with predictions consistent with other logs and to correct the sonic log for mud-filtrate invasion in the hydrocarbon-bearing sand. This editing and correction of the P-wave transit times improved the normal-incidence well tie significantly. The recordings from a monopole source severely underestimated the S-wave transit times in soft shale formations, including the reservoir seal, where the S-wave velocity was lower than the P-wave velocity in the drilling mud. The petrophysical model predicted an S-wave log that matched the valid recordings and interpolated between them. The subsequent seismic modelling from the predicted S-wave log produced a class II AVO anomaly seen on the CDP gathers around the well.     

基于复数道地震记录的匹配追踪算法及其在储层预测中的应用

匹配追踪算法是一种高精度的非线性的时频分析方法.该方法自从提出以来在地球物理学领域得到了广泛的应用.在以前的应用中,作者们都是在实信号的基础上对算法原理作了适当的简化.本文在复数地震道的基础上,基于复Morlet子波,阐述了匹配追踪算法的原理及其算法实现过程,并给出了具体的算法流程图.通过与传统的时频分析的比较,验证了该算法具有较高的时频分析精度.对于实际资料的应用表明,该算法能够很好的显示储层的位置与边界.     

The attenuated Ricker wavelet basis for seismic trace decomposition and attenuation analysis

The widely used wavelets in the context of the matching pursuit are mostly focused on the time–frequency attributes of seismic traces. We propose a new type of wavelet basis based on the classic Ricker wavelet, where the quality factor Q is introduced. We develop a new scheme for seismic trace decomposition by applying the multi-channel orthogonal matching pursuit based on the proposed wavelet basis. Compared with the decomposition by the Ricker wavelets, the proposed method could use fewer wavelets to represent the seismic signal with fewer iterations. Besides, the quality factor of the subsurface media could be extracted from the decomposition results, and the seismic attenuation could be compensated expediently. We test the availability of the proposed methods on both synthetic seismic record and field post-stack data.     

东祁连山北缘断裂带基于深度学习的密集台阵地震事件快速检测与定位研究

为监测东祁连山北缘断裂带附近的地震活动性,布设包含240台短周期地震仪的面状密集台阵,进行约30 d的连续观测。首先使用基于深度学习的多台站地震事件检测算法（CNNDetector）进行地震事件检测,然后使用震相拾取网络（PhaseNet）对地震事件进行P波和S波到时拾取,其次使用震相关联算法（REAL）进行震相关联及初定位,最后使用双差定位（hypoDD）进行地震重定位,最终的精定位地震目录中共有517个地震。在密集台阵观测期间,中国地震台网正式地震目录中共有39个位于台阵内的地震事件,相比而言,密集台阵检测到大量小于0级的地震。因此通过布设密集台阵,可提高活动断裂微地震活动性的监测能力。与历史地震空间分布相比,密集台阵地震精定位分布具有较好的一致性,表现出更明显的线性分布特征。基于地震分布,发现研究区域存在与地表断层迹线走向不同的隐伏活跃断裂。     

小波变换的过零点特性与地震勘探信号的信噪比和分辨率

利用小波变换研究地震勘探信号小波变换的过零点特性，本文提出了用小波变换的过零点特性和地震勘探信号相邻道的横向相关性提高信号分辨率和信噪比的新方法．该方法包括两个主要步骤：①利用相邻地震道信号具有很好相关性，而噪音相关性差的特点以及小波变换的过零点特性得到有效反射波同相轴随空间坐标的变化信息．②利用奇异值分解和最小二乘（SVD-TLS）方法沿同相轴对振幅进行多项式拟合去噪并增加信号高频提高信号分辨率．     

High‐definition frequency decomposition

Spectral decomposition is a widely used technique in analysis and interpretation of seismic data. According to the uncertainty principle, there exists a lower bound for the joint time–frequency resolution of seismic signals. The highest temporal resolution is achieved by a matching pursuit approach which uses waveforms from a dictionary of functions (atoms). This method, in its pure mathematical form can result in atoms whose shape and phase have no relation to the seismic trace. The high‐definition frequency decomposition algorithm presented in this paper interleaves iterations of atom matching and optimization. It divides the seismic trace into independent sections delineated by envelope troughs, and simultaneously matches atoms to all peaks. Co‐optimization of overlapping atoms ensures that the effects of interference between them are minimized. Finally, a second atom matching and optimization phase is performed in order to minimize the difference between the original and the reconstructed trace. The fully reconstructed traces can be used as inputs for a frequency‐based reconstruction and red–green–blue colour blending. Comparison with the results of the original matching pursuit frequency decomposition illustrates that high‐definition frequency decomposition based colour blends provide a very high temporal resolution, even in the low‐energy parts of the seismic data, enabling a precise analysis of geometrical variations of geological features.     

Seismic interferometry by tangent‐phase correction

We present a modified interferometry method based on local tangent‐phase analysis, which corrects the cross‐correlated data before summation. The approach makes it possible to synthesize virtual signals usually vanishing in the conventional seismic interferometry summation. For a given pair of receivers and a set of different source positions, a plurality of virtual traces is obtained at new stationary projected points located along the signal wavefronts passing through the real reference receiver. The position of the projected points is estimated by minimizing travel times using wavefront constraint and correlation‐signal tangent information. The method uses mixed processing, which is partially based on velocity‐model knowledge and on data‐based blind interferometry. The approach can be used for selected events, including reflections with different stationary conditions and projected points with respect to those of the direct arrivals, to extend the interferometry representation in seismic exploration data where conventional illumination coverage is not sufficient to obtain the stationary‐phase condition. We discuss possible applications in crosswell geometry with a velocity anomaly and a time lapse.     

Restoration of missing offsets by parabolic Radon transform1

Restoration of missing offsets and trace interpolation is an interesting and important problem in seismic data processing. Based on the parabolic Radon transform, a method is presented for missing offset restoration, resampling and regularization of prestack individual CMP gathers. The method is also valid for resampling spatially aliased seismic data. The method is based on the parabolic assumption of the seismic events which is generally verified after a partial NMO correction in the CMP organization of the data. The essence of the method consists of a band-limited forward parabolic Radon transform of the data containing zero traces at the missing offset locations. The curvature range is chosen to map properly the coherent energy while the zero traces map beyond that range. After inverse transform the originally zero traces are partly filled with information. Several iterations of forward and inverse transform, every time replacing the zero traces in the original gather with the partially reconstructed ones, almost fully restore the zero traces. Efficient and fast algorithms can be built up to process data having a uniform geometry. Examples on synthetic as well as on field data demonstrate clearly the robustness of the method.     

利用叠前CMP资料估计介质品质因子

本文提出了一种利用叠前CMP道集资料估计介质品质因子的方法,同时提出利用地震道间有效信号的相干性估计层位信息的方法.首先采用具有4个待定系数的函数去逼近震源子波,利用黏弹介质中单程波传播理论推导出地震子波包络峰值处的瞬时频率(EPIF)与不同偏移距处走时的关系;其次利用该关系,外推出该同相轴零偏移距处的EPIF,用小波域包络峰值处瞬时频率法(WEPIF)结合层位信息估计Q值.合成数据的结果表明,EPIFVO法(子波包络峰值瞬时频率随偏移距变化)无边界效应,计算精度相对较高;将该方法用于实际资料算例,结果表明,衰减强弱与储层的吸收有较好的对应关系.