DMO技术在苏北二维地震资料处理中的应用效果

由于苏北盆地历经了多次地质运动,造成地层构造复杂,倾角多变,断层多等系列地质因素。以往这一地区的地震资料处理,通常采用的是常规水平叠加,作为地层界面水平的情况是较为适用的。当地层界面倾斜时,它不能实现真正的共反射点叠加,在构造复杂地区,会使叠加剖面分辨率降低。DMO技术（倾角时差校正技术）是目前资料处理中的一项非常重要的技术,它主要是针对叠前能识别出倾角的道集中（如共炮检距道集）,经特定的速度分析和倾角时差计算,来消除正常时差校正（NMO）无法消除的地层倾角的影响,实现叠前部分偏移,得到最佳叠加速度,以达到提高剖面叠加质量的目的。     

利用地震资料进行煤层厚度解释预测

利用C++语言基于Windows操作环境开发了适用于煤田地震资料解释的煤层厚度辅助解释系统。该系统采用人工神经网络非线性反演方法对煤层厚度变化进行解释,对煤层厚度和地震属性参数之间的非线性关系给出了定量描述,具有较高的解释精度。理论模型和实际资料试算结果表明:利用该系统进行煤层厚度解释,得到了较好的地质效果。     

Modeling of seismic hazard for Turkey using the recent neotectonic data

Recent developments in the neotectonic framework of Turkey introduced new tectonic elements necessitating the reconstruction of Turkey's seismic hazard map. In this regard, 14 seismic source zones were delineated. Maximum earthquake magnitudes for each seismic zones were determined using the fault rupture length approximation. Regression coefficients of the earthquake magnitude–frequency relationships for the seismic zones were compiled mostly from earlier works. Along with these data, a strong ground motion attenuation relationship developed by Joyner and Boore [Joyner, W.B., Boore, D.M., 1988. Measurement, characterization, and prediction of strong ground motion. Earthquake Engineering and Soil Dynamics, 2. Recent Advances Ground Motion Evaluation, pp. 43–102.] was utilized to model the seismic hazard for Turkey using the probabilistic approach. For the modeling, the “earthquake location uncertainty” concept was employed. A grid of 5106 points with 0.2° intervals was constituted for the area encompassed by the 25–46°E longitudes and 35–43°N latitudes. For the return periods of 100 and 475 years, the peak horizontal ground acceleration (pga) in bedrock was computed for each grid point. Isoacceleration maps for the return periods of 100 and 475 years were constructed by contouring the pga values at each node.     

History matching time-lapse seismic data using the ensemble Kalman filter with multiple data assimilations

The ensemble Kalman filter (EnKF) has become a popular method for history matching production and seismic data in petroleum reservoir models. However, it is known that EnKF may fail to give acceptable data matches especially for highly nonlinear problems. In this paper, we introduce a procedure to improve EnKF data matches based on assimilating the same data multiple times with the covariance matrix of the measurement errors multiplied by the number of data assimilations. We prove the equivalence between single and multiple data assimilations for the linear-Gaussian case and present computational evidence that multiple data assimilations can improve EnKF estimates for the nonlinear case. The proposed procedure was tested by assimilating time-lapse seismic data in two synthetic reservoir problems, and the results show significant improvements compared to the standard EnKF. In addition, we review the inversion schemes used in the EnKF analysis and present a rescaling procedure to avoid loss of information during the truncation of small singular values.     

Converted phases analysis of the Campi Flegrei caldera using active seismic data

We applied a type of depth migration for converted seismic phases to active and passive seismic data sets from the Campi Flegrei volcanic region in southern Italy. The migration method is based on the diffraction summation migration technique. Travel times to grid points in a volume were calculated in smooth P and S-velocity models and trace energy near the calculated converted phase time was stacked over multiple sources at one receiver. Weighting factors based on Snell's Law at an arbitrarily oriented local interface were applied to better focus trace energy. PP reflection images from the active data set provide increased detail to images of the caldera rim from other studies. The migrated images also show features near 2–3 km in depth beneath Pozzuoli City, which may be associated with an over-pressured gas volume, as suggested by other geophysical investigations. Possible deeper features near 4 km depth may be related to the presence of the carbonate basement or may image a previously undetected feature, such as a small body of strongly thermometamorphosed volcanic rock. The current passive earthquake data set from the 1984 ground uplift episode was not well suited to the converted phases analysis due to narrow P–S windows and high noise levels in the traces. However, two stations provide confirmation and extension of imaged features in the active data.     

Bayesian techniques for seismic hazard assessment using imprecise data

Some Bayesian methods of dealing with inaccurate or vague data are introduced in the framework of seismic hazard assessment. Inaccurate data affected by heterogeneous errors are modeled by a probability distribution instead of the usual value plus a random error representation; these data are generically called imprecise. The earthquake size and the number of events in a certain time are modeled as imprecise data. Imprecise data allow us to introduce into the estimation procedures the uncertainty inherent in the inaccuracy and heterogeneity of the measuring systems from which the data were obtained. The problem of estimating the parameter of a Poisson process is shown to be feasible by the use of Bayesian techniques and imprecise data. This background technique can be applied to a general problem of seismic hazard estimation. Initially, data in a regional earthquake catalog are assumed imprecise both in size and location (i.e errors in the epicenter or spreading over a given source). By means of scattered attenuation laws, the regional catalog can be translated into a so-called site catalog of imprecise events. The site catalog is then used to estimate return periods or occurrence probabilities, taking into account all sources of uncertainty. Special attention is paid to priors in the Bayesian estimation. They can be used to introduce additional information as well as scattered frequency-size laws for local events. A simple example is presented to illustrate the capabilities of this methodology.     

利用地震P波确定煤层瓦斯富集带的分布

查明煤层瓦斯富集区域,对可能的瓦斯突出点进行预报,是当前煤矿生产中亟待解决的重要课题。利用地震P波对裂缝性地层所表现出的方位各向异性特征,根据地震属性随方位角变化可以预测裂隙发育方向和密度的基本原理,应用多种地震P波方位属性预测裂隙发育带。通过对淮南张集煤矿西三采区三维地震P波资料的处理,获得6个方位地震数据体,从中提取多种与煤层和围岩裂隙相关的地震属性,并计算出裂隙的发育方向和密度,为确定瓦斯富集带的分布提供了一种新的途径。     

Random seismic noise attenuation data using the discrete and the continuous wavelet transforms

In this paper, we propose a technique of random noise attenuation from seismic data using discrete and continuous wavelet transforms. Firstly, the discrete wavelet transform (DWT) is applied to denoise seismic data using the threshold method. After, we calculate the continuous wavelet transform of the denoised seismic seismogram, the final denoised seismic seismogram is the continuous wavelet transform coefficients at the lower scale. Application to a synthetic seismogram shows the robustness of the proposed tool for random noise attenuation. Application to real vertical seismic profile recorded in Algeria clearly shows the efficiency of the proposed tool for random noise attenuation.     

Crustal studies of the Koyna region using explosion data from deep seismic soundings

This study is based on the seismic data collected as a result of explosions carried out during the 1976 and 1978 Deep Seismic Sounding (DSS) field operations in the Koyna region. These shots were exploded from twelve shot points by the National Geophysical Research Institute along the Guhagar-Chorochi and Kelsi-Loni profiles.Refraction studies of the records reveal a two-layered crust. The top layer consists 17 km of granite and the second layer 19 km of basalt, giving the average depth of the Moho as 36 km in the region. The velocities of the phases P_g, P* and P_n have been computed as 5.82 ± 0.01, 6.61 ± 0.05 and 8.23 ± 0.05 km/sec respectively and those of S_g, S* and S_n as 3.41 ± 0.00, 4.09 ± 0.07 and 4.60 ± 0.08 km/sec respectively. The shear wave velocity in the basement rock has been found to be lower in comparison with other regions of the peninsular India.In some cases reflections were recorded both from the Moho as well as from the intermediate layer. These reveal a crustal thickness of 39 km with 19 km of granitic and 20 km of basaltic layers.Coda signal durations from DSS explosions recorded by microearthquake seismographs indicate a lateral heterogeneity in the crust on either side of Karad in an east-west direction.     

Active tectonics of the Yizre'el valley, Israel, using high-resolution seismic reflection data

We present a series of high-resolution seismic reflection lines across the Yizre'el valley, which is the largest active depression in Israel, off the main trend of the Dead Sea rift. The new seismic reflection data is of excellent quality and shows that the valley is dissected into numerous small blocks, separated by active faults. The Yizre'el valley is found to consist of a series of half grabens, rather than a single half graben, or a symmetrical graben. The faults are generally vertical and appear to have a dominant strike-slip component, but some dip-slip is also evident. A marked zone of compression near Megido is associated with the intersection of the two largest faults in the valley, the Carmel fault and the Gideon fault. Variable trend of the faults reflects the complexity of the local geology along the boundary between the wide NW–SE trending Farah–Carmel fault zone and the E–W trending basins and ranges in the Lower Galilee. This tectonic complexity is likely to result from a highly variable stress pattern, modified by the structures inside it. Normal faulting in the valley occurred at an early stage of its development as a tectonic depression. However, strike-slip motion on the Carmel fault, and possibly also on some of the other faults, appears to have started together with the onset of normal faulting. Earthquake hazard in the area appears to be uniform as faults are distributed throughout the Yizre'el valley.     

基于小波变换的阈值降噪方法在地震资料处理中的应用

这里主要讨论了小波变换阈值降噪方法的基本原理,并且将广泛使用的通用阈值方法与Birge-Massart阈值方法进行了比较,展现了Birge-Massart阈值选取方法的优越性。通过对模拟数据和实际地震资料的处理结果表明,基于小波变换的Birge-Massart阈值降噪方法对地震信号降噪具有很好的效果。     

利用地震资料预测煤层古河道冲刷带

为预测沁水盆地西缘某矿区煤层古河道冲刷带,开展了地震波数值正演模拟,利用地震属性分析技术,结合波阻抗反演技术、三维可视化技术等手段,研究了煤层古河道冲刷带分布的特征。后期钻孔验证表明,利用地震资料预测煤层古河道冲刷带的方法是可行的。     

DMO方法综述

DMO（倾角时差校正）是一种将非零偏移距地震道转换为零偏移距地震道的数据处理方法，主要针对地层倾角不为零的情况。通过应用DMO，共中心点道集变成以自激自收（垂直入射）定义的共深度点道集。DMO是地震资料处理中一个非常重要的流程，对DMO的发展状况和一些主要算法进行总结，各向异性DMO将有望成为今后的发展方向。     

ARMA order model detection using minimum of Kurtosis: application on seismic data

The objective of this study is to find the order and coefficients of non-low-phase causal filters for ARMA (auto regressive moving average) filter model, using the Kurtosis minimization criterion. This method is based on the Kurtosis calculation of the treated sample at the input level and its identification at the output of the ARMA model. For this purpose, the order and coefficients of the AR (auto regressive) part are identified using the Yule-Walker algorithm at order two and then extended to order four. To obtain the MA (moving average) part, the AR components are calculated at first from the ARMA filter by deconvolution. Then, spectrally equivalent and minimum phase (SEMP) MA filter is identified using the Durbin algorithm at second and fourth order. Finally, the correct filter is found when the Kurtosis value of the output ARMA filter reconstituted is the closest to the Kurtosis of introduced signal. The proposed method is then tested on simulated processes and applied to real seismic data to perform blind deconvolution and obtain the reflectivity coefficients of subsoil studied.     

综合应用测井及地震资料建立复杂断块地质模型

针对滩海油田具有完钻井少,以地震资料为主的特点,采用井震结合的方法进行地质建模。复杂断块油藏具有许多内部断层,充分利用地震解释层位和断层成果建立合理的构造模型。采用“多步建模”的思路,选用序贯高斯同位协同模拟算法建立储层参数模型。以测井资料为硬数据,地震反演数据为软数据,先后建立自然伽玛模型、泥质含量模型和砂泥岩相模型,然后以砂泥岩相模型为约束,建立储层参数模型。经研究表明,综合使用地震解释成果和三维地震数据体,有助于建立准确、合理的构造模型;而采用“多步建模”的思路,利用序贯高斯同位协同模拟算法,并综合测井资料、地震资料建立储层参数模型,有助于从地质的角度对储层参数模拟进行约束,提高储层随机建模精度,保证储层参数模型的准确和可靠性。     

2-D velocity structure in Kerala-Konkan basin using traveltime inversion of seismic data

The existence of gas-hydrates in marine sediments increases the seismic velocity, whereas even a small amount of underlying free-gas reduces the velocity considerably. The change in velocities against the background (without gas-hydrates and free-gas) velocity can be used for identification and assessment of gas-hydrates. Traveltime inversion of identifiable reflections from large offset multi channel seismic (MCS) experiment is an effective method to derive the 2-D velocity structure in an area. We apply this method along a seismic line in the Kerala-Konkan (KK) offshore basin for delineating the gas-hydrates and free-gas bearing sediments across a bottom simulating reflector (BSR). The result reveals a four layer 2-D shallow velocity model with the topmost sedimentary layer having velocity of 1,680–1,740 m/s and thickness of 140–190 m. The velocity of the second layer of uniform thickness (110 m) varies from 1,890 to 1,950 m/s. The third layer, exhibiting higher velocity of 2,100–2,180 m/s, is interpreted as the gas-hydrates bearing sediment, the thickness of which is estimated as 100 to 150 m. The underlying sedimentary layer shows a reduction in seismic velocity between 1,620 to 1,720 m/s. This low-velocity layer with 160–200 m thickness may be due to the presence of free-gas below the gas-hydrates layer.     

Denoising and improving the quality of seismic data using combination of DBM filter and FX deconvolution

Seismic data denoising, random noise attenuation (RNA) and spike-like noise suppression, is a main consideration for improving the quality of records. RNA could increase signal to noise ratio (S/N) to avoid misinterpretation of seismic data. In this research, a novel method is created by using the combination of frequency-offset deconvolution (FXD) and decision-based median (DBM) filter for RNA from seismic data. The method is applied in two main phases; FXD is focused to remove the Gaussian noise and DBM filter is focused to attenuate the impulsive noise and spikes. To implement and verify the method, three types of data are used: two synthetic models (a model with linear events and a model with hyperbolic events) and an observed seismic section. The ability of the proposed method (FXD-DBM) in comparison of applying each in seismic RNA application is proven. The noise level is reduced obviously, and hence, the S/N of all examined seismic records is increased considerably after denoising by the combination of FX deconvolution and DBM filter. About the real seismic section, suppressing random noise and spikes show up improving the seismic reflector continuity and hence enhancing the interpretability of data. Moreover, some masked events by random noise are clarified in different parts of data after denoising using the planned method.     

Spatial analysis of two-dimensional orientation data

A frequently encountered problem in structural geology is the analysis and interpretation of directional data with a geographical distribution. Such data may arise, for example, as samples of fracture strikes collected at a number of locations (rocky outcrops, road cuttings, etc.), trends of geological structures mapped in underground mines where they are exposed in the tunnel roof, or as lineaments annotated on aerial or satellite photographs. From such data, it is of interest to determine the structural domains in the area, that is, to partition the area into subareas of homogeneous structural pattern. This paper presents a method for performing such an analysis, and applies it to a data set collected from an underground coal mine.     

Laboratory seismic measurements: an aid in the interpretation of seismic field data

This contribution reports data and evidence from laboratory measurements of compressional and shear-wave velocities and mineralogical, microstructural and chemical analyses in crustal rocks that may provide an aid to understanding and interpreting geophysical data. It is shown that the correlation of laboratory seismic data with stratigraphic field data can provide the necessary data base for constructing a seismic model of the Earth's crust in a particular geological area.     

Wavefield analysis of crosswell seismic data

Despite the numerous advantages of crosswell seismic data over surface seismic data, crosswell seismic geophysics is still underutilized and underdeveloped. The factors limiting the full utilization of crosswell data include the lack of standardized methods for processing and imaging the data. This is because crosswell data is not completely understood. To improve the understanding of crosswell data, we performed acoustic and elastic modeling of a west Texas carbonate oilfield data using finite difference methods and crosswell geometry. To account for the different wave modes in the field data, we decomposed the full data into its constituent wave modes. Results of the forward modeling show that elastic synthetic data is a better representation of crosswell field data than the popular acoustic synthetic data. Wavefield decomposition gave insight into the time-space kinematics behavior of the different wave modes that constitute the full data. Overall, the study improved our understanding of crosswell field data. The learning from this study has been utilized to perform data-driven reflection enhancement processing where the discerned characteristic of different seismic arrival is utilize to suppress unwanted and enhanced the desired wave modes. The processing reduced the complex data to only up-going P-P reflections that can be imaged to reveal the subtle geological structures of the oilfield.