共查询到20条相似文献，搜索用时 15 毫秒
1.
2.
3.
提出一种可对起伏地表采集的三维地震资料直接进行偏移成像的叠前时间偏移方法和流程.它用两个等效速度描述近地表和上覆层对地震波传播的影响,可对炮、检点不在同一水平面的三维地震资料直接进行叠前时间偏移处理.该方法不对近地表地震波传播做垂直出、入射假定,因此可适应高速层出露等不存在明显低、降速带情况.描述近地表和上覆层的两个等效速度参数可依据偏移道集的同相轴是否平直来确定,避免了确定近地表速度的困难;而对已知近地表速度的情况,则可进一步修正近地表速度,获得更好的成像效果.用三维起伏地表的理论数据和中国东部某工区实际数据验证了所发展方法和处理流程的有效性和实用性. 相似文献
4.
5.
Land seismic data quality can be severely affected by near‐surface anomalies. The imprint of a complex near‐surface can be removed by redatuming the data to a level below the surface, from where the subsurface structures are assumed to be relatively smooth. However, to derive a velocity‐depth model that explains the propagation effects of the near‐surface is a non‐trivial task. Therefore, an alternative approach has been proposed, where the redatuming operators are obtained in a data‐driven manner from the reflection event related to the datum. In the current implementation, the estimation of these redatuming operators is done in terms of traveltimes only, based on a high‐frequency approximation. The accompanying amplitudes are usually derived from a local homogeneous medium, which is obviously a simplification of reality. Such parametrization has produced encouraging results in the past but cannot completely remove the near‐surface complexities, leaving artefacts in the redatumed results. In this paper we propose a method that estimates the redatuming operators directly from the data, i.e., without using a velocity model, in a full waveform manner, such that detailed amplitude and phase variations are included. The method directly outputs the inverse propagation operators that are needed for true‐amplitude redatuming. Based on 2D synthetic data it is demonstrated that the resulting redatuming quality is improved and artefacts are reduced. 相似文献
6.
In this case study we consider the seismic processing of a challenging land data set from the Arabian Peninsula. It suffers from rough top‐surface topography, a strongly varying weathering layer, and complex near‐surface geology. We aim at establishing a new seismic imaging workflow, well‐suited to these specific problems of land data processing. This workflow is based on the common‐reflection‐surface stack for topography, a generalized high‐density velocity analysis and stacking process. It is applied in a non‐interactive manner and provides an entire set of physically interpretable stacking parameters that include and complement the conventional stacking velocity. The implementation introduced combines two different approaches to topography handling to minimize the computational effort: after initial values of the stacking parameters are determined for a smoothly curved floating datum using conventional elevation statics, the final stack and also the related residual static correction are applied to the original prestack data, considering the true source and receiver elevations without the assumption of nearly vertical rays. Finally, we extrapolate all results to a chosen planar reference level using the stacking parameters. This redatuming procedure removes the influence of the rough measurement surface and provides standardized input for interpretation, tomographic velocity model determination, and post‐stack depth migration. The methodology of the residual static correction employed and the details of its application to this data example are discussed in a separate paper in this issue. In view of the complex near‐surface conditions, the imaging workflow that is conducted, i.e. stack – residual static correction – redatuming – tomographic inversion – prestack and post‐stack depth migration, leads to a significant improvement in resolution, signal‐to‐noise ratio and reflector continuity. 相似文献
7.
Wave‐equation redatuming can be a very efficient method of overcoming the overburden imprint on the target area. Owing to the growing amount of 3D data, it is increasingly important to develop a feasible method for the redatuming of 3D prestack data. Common 3D acquisition designs produce relatively sparse data sets, which cannot be redatumed successfully by applying conventional wave‐equation redatuming. We propose a redatuming approach that can be used to perform wave‐equation redatuming of sparse 3D data. In this new approach, additional information about the medium velocity below the new datum is included, i.e. redatumed root‐mean‐square (RMS) velocities, which can be extracted from the input data set by conventional velocity analysis, are used. Inclusion of this additional information has the following implications: (i) it becomes possible to simplify the 4D redatuming integral into a 2D integral such that the number of traces needed to calculate one output time sample and the computational effort are both reduced; (ii) the information about the subsurface enables an infill of traces which are needed for the integral calculation but which are missing in the sparse input data set. Two tests applying this new approach to fully sampled 2D data show satisfactory results, implying that this method can certainly be used for the redatuming of sparse 3D data sets. 相似文献
8.
提出了一种新的叠前时间偏移方法和流程,可不必应用野外静校正,直接对起伏地表采集的地震数据进行叠前时间偏移.本文采用输入道成像方式,通过基于稳相点原理给出单道数据的走时和振幅计算方法,发展了一个表驱动的叠前时间偏移算法.偏移方法可依据同相轴是否被拉平确定叠加速度和修正近地表速度模型,也可依据拟成像的构造倾角,自适应地确定偏移孔径;后者既减少了偏移计算量,也压制了偏移噪声.文中用二维起伏地表的断陷盆地模型的理论数据验证了所发展方法的成像效果. 相似文献
9.
由于地面微地震监测台站布设在地表,会受到地表起伏、低降速带厚度和速度变化的影响,降低了微地震事件的识别准确度和定位精度,限制了地面微地震监测技术在复杂地表地区的应用.因此,将三维地震勘探技术的思路引入到地面微地震监测中,提出了三维地震与地面微地震联合校正方法,将油气勘探和开发技术更加紧密地结合在一起.根据三维地震数据和低降速带测量数据,通过约束层析反演方法建立精确的近地表速度模型,将地面微地震台站从起伏地表校正到高速层中的平滑基准面上,有效消除复杂近地表的影响.其次,根据射孔数据和声波测井速度信息,通过非线性反演方法建立最优速度模型,由于已经消除复杂近地表的影响,在进行速度模型优化时不需要考虑近地表的影响,因而建立的速度模型更加准确.最后,在精确速度模型的基础上,通过互相关方法求取剩余静校正量,进一步消除了复杂近地表和速度模型近似误差的影响.三维地震与地面微地震联合校正方法采用逐步校正的思路,能够有效消除复杂近地表的影响,提高微地震数据的品质和速度模型的精确度,保证了微地震事件的定位精度,具有良好的应用前景. 相似文献
10.
Static shifts from near‐surface inhomogeneities very often represent the key problem in the processing of seismic data from arid regions. In this case study, the deep bottom fill of a wadi strongly degrades the image quality of a 2D seismic data set. The resulting static and dynamic problems are solved by both conventional and common‐reflection‐surface (CRS) processing. A straightforward approach derives conventional refraction statics from picked first breaks and then goes through several iterations of manual velocity picking and residual statics calculation. The surface‐induced static and dynamic inhomogeneities, however, are not completely solved by these conventional methods. In CRS processing, the local adaptation of the CRS stacking parameters results in very detailed dynamic corrections. They resolve the local inhomogeneities that were not detected by manual picking of stacking velocities and largely compensate for the surface‐induced deterioration in the stack. The subsequent CRS residual statics calculations benefit greatly from the large CRS stacking fold which increases the numbers of estimates for single static shifts. This improves the surface‐consistent averaging of static shifts and the convergence of the static solution which removes the remaining static shifts in the 2D seismic data. The large CRS stacking fold also increases the signal‐to‐noise ratio in the final CRS stack. 相似文献
11.
Paul L. Stoffa A. Mauffret M. Truchan Peter Buhl 《Earth and Planetary Science Letters》1981,53(2):131-146
Multichannel seismic data in the Aruba Gap region near JOIDES/DSDP Site 153 verify the presence of a deep sub-B″ reflection. One multichannel seismic line trends NE-SW on and along the edge of Beata Ridge and passes within 1 km of Site 153, and another line runs N-S across the entire Aruba Gap with the drill site 4 km east from its northern end. Closely spaced velocity analyses indicate the presence of deep primary reflection events and enable calculation of interval velocities between the A″-B″ marker horizons. Deconvolved, true amplitude, normal incidence profiles sharply delineate the A″-B″ marker horizons and indicate the presence of the sub-B″ reflection event. On the NE-SW line, this deep reflector is best described as a “diffuse” discontinuous zone, relatively horizontal, lying about 0.8 seconds of two-way travel time below the B″ reflector, with an interval velocity of approximately 5.0 km/s between Horizon B″ and this reflection. The N-S line is more complex since the sub-B′' reflection event is masked by a strong internal multiple from the A″-B″ interval. In the central and western Venezuela Basin, deep primary reflections beneath Horizon B″ are also observed on the northern and western sides of what appears to be a major fault zone. This fault zone separates the smooth B″ and sub-B″ reflectors on the northern and western sides of this fault zone from what appears to be typical oceanic basement. The widespread presence of sub-B″ reflections yielding high interval velocities for the section between these events and Horizon B″ suggest that this material is probably igneous in origin. 相似文献
12.
A. J. HARDING 《Geophysical Prospecting》1985,33(5):674-695
Conventional velocity analysis, based on the ideas of rms velocity and hyperbolic reflection events in the x-t domain, is restricted in validity to near vertical incidence. Thus analysis of near-offset datasets usually requires the muting of wide-angle reflections from shallow interfaces before the rms velocities are determined. The ray-theoretical integral for the delay time τ, which depends on the slowness p and the velocity function, is valid for all angles. The wide-angle reflections can be used to improve the accuracy of the derived velocity function in the near surface region, if the recorded x-t data are mapped into the τ-p domain. By representing the velocity function between reflectors as a series of gradient zones, i.e. regions with a uniform increase in velocity with depth, the recovery of the velocities may be posed as a matrix linear inverse problem for the slopes of the gradient zones. In order to convert the problem to a linear one, the velocity discontinuities at the reflecting interfaces must be fixed in advance. Their positions are based on the behaviour of the τ-p map of the data. Finding a stable velocity model may require several iterations with the reflecting interfaces at different positions. An understanding of the workings of the inversion algorithm allied with an analysis of the causes of instability aids the search for a stable model. 相似文献
13.
14.
Peter Bergmann Can YangStefan Lüth Christopher JuhlinCalin Cosma 《Journal of Applied Geophysics》2011,75(1):124-139
The Ketzin project provides an experimental pilot test site for the geological storage of CO2. Seismic monitoring of the Ketzin site comprises 2D and 3D time-lapse experiments with baseline experiments in 2005. The first repeat 2D survey was acquired in 2009 after 22 kt of CO2 had been injected into the Stuttgart Formation at approximately 630 m depth. Main objectives of the 2D seismic surveys were the imaging of geological structures, detection of injected CO2, and comparison with the 3D surveys. Time-lapse processing highlighted the importance of detailed static corrections to account for travel time delays, which are attributed to different near-surface velocities during the survey periods. Compensation for these delays has been performed using both pre-stack static corrections and post-stack static corrections. The pre-stack method decomposes the travel time delays of baseline and repeat datasets in a surface consistent manner, while the latter cross-aligns baseline and repeat stacked sections along a reference horizon.Application of the static corrections improves the S/N ratio of the time-lapse sections significantly. Based on our results, it is recommended to apply a combination of both corrections when time-lapse processing faces considerable near-surface velocity changes. Processing of the datasets demonstrates that the decomposed solution of the pre-stack static corrections can be used for interpretation of changes in near-surface velocities. In particular, the long-wavelength part of the solution indicates an increase in soil moisture or a shallower groundwater table in the repeat survey.Comparison with the processing results of 2D and 3D surveys shows that both image the subsurface, but with local variations which are mainly associated to differences in the acquisition geometry and source types used. Interpretation of baseline and repeat stacks shows that no CO2 related time-lapse signature is observable where the 2D lines allow monitoring of the reservoir. This finding is consistent with the time-lapse results of the 3D surveys, which show an increase in reflection amplitude centered around the injection well. To further investigate any potential CO2 signature, an amplitude versus offset (AVO) analysis was performed. The time-lapse analysis of the AVO does not indicate the presence of CO2, as expected, but shows signs of a pressure response in the repeat data. 相似文献
15.
A single‐layer model of the near surface throughout the Kingdom of Saudi Arabia is available. While this simple model suffices for most areas and large subsurface structures, it fails in situations where the surface topography is complex, the base of weathering is below the datum, or where the time structural closure is less than the uncertainty in the static correction. In such cases, multiple‐layered models that incorporate velocities derived from analysis of first arrivals picked from seismic shot records have proved to be successful in defining the lateral heterogeneity of the near surface. The additional velocity information obtained from this first‐arrival analysis (direct as well as refracted arrivals) vastly improves the velocity–depth model of the near surface, regardless of the topography. Static corrections computed from these detailed near‐surface velocity models have significantly enhanced subsurface image focusing, thereby reducing the uncertainty in the closure of target structures. Other non‐seismic methods have been used either to confirm qualitatively or to enhance the layer models previously mentioned. Gravity data may be particularly useful in sandy areas to confirm general structure, while geostatistical modelling of vibrator base‐plate attributes has yielded information that enhances the velocity field. In the global context, exploration targets of the oil and gas industry are seeking smaller and lower relief‐time structures. Thus, near‐surface models will need to enhance and integrate these methods, particularly in areas where the assumption of flat‐lying near‐surface layers cannot be met. 相似文献
16.
In many land seismic situations, the complex seismic wave propagation effects in the near‐surface area, due to its unconsolidated character, deteriorate the image quality. Although several methods have been proposed to address this problem, the negative impact of 3D complex near‐surface structures is still unsolved to a large extent. This paper presents a complete 3D data‐driven solution for the near‐surface problem based on 3D one‐way traveltime operators, which extends our previous attempts that were limited to a 2D situation. Our solution is composed of four steps: 1) seismic wave propagation from the surface to a suitable datum reflector is described by parametrized one‐way propagation operators, with all the parameters estimated by a new genetic algorithm, the self‐adjustable input genetic algorithm, in an automatic and purely data‐driven way; 2) surface‐consistent residual static corrections are estimated to accommodate the fast variations in the near‐surface area; 3) a replacement velocity model based on the traveltime operators in the good data area (without the near‐surface problem) is estimated; 4) data interpolation and surface layer replacement based on the estimated traveltime operators and the replacement velocity model are carried out in an interweaved manner in order to both remove the near‐surface imprints in the original data and keep the valuable geological information above the datum. Our method is demonstrated on a subset of a 3D field data set from the Middle East yielding encouraging results. 相似文献
17.
A redatuming operation is used to simulate the acquisition of data in new levels, avoiding distortions produced by near-surface irregularities related to either geometric or material property heterogeneities. In this work, the application of the true-amplitude Kirchhoff redatuming (TAKR) operator on homogeneous media is compared with conventional Kirchhoff redatuming (KR) operator restricted to the zero-offset case. The TAKR and the KR operators are analytically and numerically compared in order to verify their impacts on the data at a new level. Analyses of amplitude and velocity sensitivity of the TAKR and KR were performed: one concerning the difference between the weight functions and the other related to the velocity variation. The comparisons between operators were performed using numerical examples. The feasibility of the KR and TAKR operators was demonstrated not only kinematically but also dynamically for their purposes. In other words, one preserves amplitude (KR), and the other corrects the amplitude (TAKR). In the end, we applied the operators to a GPR data set. 相似文献
18.
We describe an integrated method for solving the complex near‐surface problem in land seismic imaging. This solution is based on an imaging approach and is obtained without deriving a complex near‐surface velocity model. We start by obtaining from the data the kinematics of the one‐way focusing operators (i.e. time‐reversed Green's functions) that describe propagation between the acquisition surface and a chosen datum reflector using the common‐focus‐point technology. The conventional statics solutions obtained from prior information about the near surface are integrated in the initial estimates of the focusing operators. The focusing operators are updated iteratively until the imaging principle of equal traveltime is fulfilled for each subsurface gridpoint of the datum reflector. Therefore, the seismic data is left intact without any application of time shifts, which makes this method an uncommitted statics solution. The focusing operators can be used directly for wave‐equation redatuming to the respective reflector or for prestack imaging if determined for multiple reflecting boundaries. The underlying velocity model is determined by tomographic inversion of the focusing operators while also integrating any hard prior information (e.g. well information). This velocity model can be used to perform prestack depth imaging or to calculate the depth of the new datum level. We demonstrate this approach on 2D seismic data acquired in Saudi Arabia in an area characterized by rugged topography and complex near‐surface geology. 相似文献
19.
Geostatistical integration of near-surface geophysical data 总被引：1，自引：0，他引：1
Accurate statics calculation and near‐surface related noise removal require a detailed knowledge of the near‐surface velocity field. Conventional seismic surveys currently are not designed to provide this information, and 3D high‐resolution reflection/refraction acquisition is not feasible for large survey areas. Satellite images and vibrator plate attributes are dense low‐cost data, which can be used in spatially extrapolating velocities from sparse uphole data by geostatistics. We tested this approach in two different areas of Saudi Arabia and found that the optimal recipe depends on the local geology. 相似文献
20.
Tariq Alkhalifah 《Geophysical Prospecting》2005,53(5):643-653
Imaging pre‐salt reflections for data acquired from the coastal region of the Red Sea is a task that requires prestack migration velocity analysis. Conventional post‐stack time processing lacks the lateral inhomogeneity capability, necessary for such a problem. Prestack migration velocity analysis in the vertical time domain reduces the velocity–depth ambiguity that usually hampers the performance of prestack depth‐migration velocity analysis. In prestack τ‐migration velocity analysis, the interval velocity model and the output images are defined in τ (i.e. vertical time). As a result, we avoid placing reflectors at erroneous depths during the velocity analysis process and thus avoid inaccurately altering the shape of the velocity model, which in turn speeds up the convergence to the true model. Using a 1D velocity update scheme, the prestack τ‐migration velocity analysis produces good images of data from the Midyan region of the Red Sea. For the first seismic line from this region, only three prestack τ‐migration velocity analysis iterations were required to focus pre‐salt reflections in τ. However, the second line, which crosses the first line, is slightly more complicated and thus required five iterations to reach the final, reasonably focused, τ‐image. After mapping the images for the two crossing lines to depth, using the final velocity models, the placements of reflectors in the two 2D lines were consistent at their crossing point. Some errors occurred due to the influence of out‐of‐plane reflections on 2D imaging. However, such errors are identifiable and are generally small. 相似文献