首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 15 毫秒
1.
We present an extension of the Common Reflection Surface (CRS) stack that provides support for an arbitrary top surface topography. CRS stacking can be applied to the original prestack data without the need for any elevation statics. The CRS-stacked zero- offset section can be corrected (redatumed) to a given planar level by kinematic wave field attributes. The seismic processing results indicate that the CRS stacked section for rugged surface topography is better than the conventional stacked section for S/N ratio and better continuity of reflection events. Considering the multiple paths of zero-offset rays, the method deals with reflection information coming from different dips and performs the stack using the method of dip decomposition, which improves the kinematic and dynamic character of CRS stacked sections.  相似文献   

2.
复杂地表条件下共反射面元(CRS)叠加方法研究   总被引:17,自引:3,他引:17       下载免费PDF全文
在地表地形复杂的情况下,静校正不易做好,这是制约山地资料处理质量的一个很重要的因素.复杂地表共反射面元(CRS)叠加不需对叠前数据做静校正,而且在得到叠加剖面后可以利用叠加得到的波场参数剖面实现基准面重建.地震数据的试算表明,复杂地表CRS叠加得出的剖面与常规处理剖面相比有着较高的信噪比和同相轴连续性.与水平地表CRS叠加不同的是,在复杂地表CRS叠加的时距公式中,波场三参数耦合,难以通过简化CRS道集的方法将它们全部分离并逐个优化.引入模拟退火算法后,有效地解决了这一组合优化的难题.  相似文献   

3.
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.  相似文献   

4.
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.  相似文献   

5.
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.  相似文献   

6.
When topography and low velocity zone differences vary greatly, conventional vertical static time shifts will cause wavefield distortion and influence wave equation seismic imaging for seismic data acquired on a complex near surface. In this paper, we propose an approach to datum correction that combines a joint tomography inversion with wavefield continuation to solve the static problem for seismic data on rugged acquisition topography. First, the near surface model is obtained by refracted wave tomography inversion. Second, the wavefield of sources and receivers are continued downward and upward to accomplish datum correction starting from a flat surface and locating the datum above topography. Based on the reciprocal theorem, Huygens' and Fresnel principles, the location of sources and receivers, and regarding the recorded data on the surface as a secondary emission, the sources and receivers are upward-continued to the datum above topography respectively. Thus, the datum correction using joint tomography inversion and wavefield continuation with the condition of a complex near surface is accomplished.  相似文献   

7.
We apply a redatuming methodology, designed to handle rugged topography and the presence of high‐velocity layers near the acquisition surface, to a 2D land seismic data set acquired in Saudi Arabia. This methodology is based on a recently developed prestack operator, which we call the topographic datuming operator (TDO). The TDO, unlike static corrections, allows for the movement of reflections laterally with respect to their true locations, corresponding to the new datum level. Thus, it mitigates mispositioning of events and velocity bias introduced by the assumption of surface consistency and the time‐invariant time shifts brought about by static corrections. Using the shallow velocities estimated from refracted events, the TDO provides a superior continuity of reflections and better focusing than that obtained from conventional static corrections in most parts of the processed 2D line. The computational cost of applying the TDO is only slightly higher than static corrections. The marginal additional computational cost and the possibility of estimating, after TDO redatuming, stacking velocities that are not affected by a spurious positive bias, as in the case of static corrections, are further advantages of the proposed methodology. The likelihood of strong heterogeneities in the most complex part of the line limits the applicability of any approach based upon geometrical optics; however, the TDO produces results that are slightly better than those obtained from static corrections because of its ability to partially collapse diffractions generated in the near surface.  相似文献   

8.
We review the multifocusing method for traveltime moveout approximation of multicoverage seismic data. Multifocusing constructs the moveout based on two notional spherical waves at each source and receiver point, respectively. These two waves are mutually related by a focusing quantity. We clarify the role of this focusing quantity and emphasize that it is a function of the source and receiver location, rather than a fixed parameter for a given multicoverage gather. The focusing function can be designed to make the traveltime moveout exact in certain generic cases that have practical importance in seismic processing and interpretation. The case of a plane dipping reflector (planar multifocusing) has been the subject of all publications so far. We show that the focusing function can be generalized to other surfaces, most importantly to the spherical reflector (spherical multifocusing). At the same time, the generalization implies a simplification of the multifocusing method. The exact traveltime moveout on spherical surfaces is a very versatile and robust formula, which is valid for a wide range of offsets and locations of source and receiver, even on rugged topography. In two‐dimensional surveys, it depends on the same three parameters that are commonly used in planar multifocusing and the common‐reflection surface (CRS) stack method: the radii of curvature of the normal and normal‐incidence‐point waves and the emergence angle. In three dimensions the exact traveltime moveout on spherical surfaces depends on only one additional parameter, the inclination of the plane containing the source, receiver and reflection point. Comparison of the planar and spherical multifocusing with the CRS moveout expression for a range of reflectors with increasing curvature shows that the planar multifocusing can be remarkably accurate but the CRS becomes increasingly inaccurate. This can be attributed to the fact that the CRS formula is based on a Taylor expansion, whereas the multifocusing formulae are double‐square root formulae. As a result, planar and spherical multifocusing are better suited to model the moveout of diffracted waves.  相似文献   

9.
In the case of onshore data sets, the acquired reflection events can be strongly impaired due to rough top‐surface topography and inhomogeneities in the uppermost low‐velocity layer, the so‐called weathering layer. Without accounting for these influences, the poor data quality will make data processing very difficult. Usually, the correction for the top‐surface topography is not perfect. The residuals from this correction and the influence of the weathering layers lead to small distortions along the reflection events. We integrated a residual static correction method into our data‐driven common‐reflection‐surface‐stack‐based imaging workflow to further eliminate such distortions. The moveout‐corrected traces and the stacked pilot trace are cross‐correlated to determine a final estimate of the surface‐consistent residual statics in an iterative manner. As the handling of top‐surface topography within the common‐reflection‐surface stack is discussed in a separate paper in this special issue, the corresponding residual static correction will be explained in more detail. For this purpose, the results obtained with a data set from the Arabian Peninsula will be presented.  相似文献   

10.
共反射面元叠加的应用实践   总被引:19,自引:5,他引:14       下载免费PDF全文
共反射面元(Common Reflection Surface)叠加是一种不依赖于宏观速度模型的零炮检距剖面成像方法,实现共反射面元叠加依赖于3个波场属性参数的确定,它们分别是零偏移距射线的出射角α、Normal波和Normal Incident Point波出射到地表的波前曲率半径RN和RNIP. 在CRS叠加的理论基础上,本文阐述如何在实际数据上实现CRS叠加. 首先,通过简洁的一维相关性分析在常规叠加剖面上找到对应该共反射面元的一组初始波场属性参数(α,RN,RNIP),然后在对应的叠前数据上应用最优化算法对这组参数进行优化处理,相比初始属性参数,优化后的属性参数能够更好地聚集来自地下反射层的能量,最后应用优化后的属性参数实现最优CRS叠加.  相似文献   

11.
复杂地表条件下高斯波束叠前深度偏移(英文)   总被引:6,自引:0,他引:6  
在复杂地表条件的区域,地震数据的采集和处理是一项极大的挑战。虽然可以通过静校正来消除起伏地表的影响,然而当地表高程以及近地表速度剧烈变化时,简单的垂直时移对地震波场造成的畸变会严重降低偏移成像的质量。基于射线的偏移方法可以直接在起伏地表面进行波场的延拓成像,是解决上述问题的有效手段。本文针对复杂地表条件下的高斯波束叠前深度偏移进行研究,对倾斜叠加公式进行修改,使之包含地表高程以及速度的信息,通过直接在复杂地表面进行平面波的合成,得到了一种具有更高成像精度的改进方法。首先简单介绍常规高斯波束偏移的基本原理和计算流程,并以此为基础,给出复杂地表条件下高斯波束偏移原有的实现方法以及本文的改进方法,最后通过模型和实际资料的试算验证本文方法的有效性。  相似文献   

12.
起伏地表煤田地震资料静校正   总被引:1,自引:1,他引:0       下载免费PDF全文
由于地表起伏和近地表结构变化产生的静校正问题严重影响了煤田地震资料的成像质量.为此,首先利用低速带分片拟合的广义线性反演技术进行折射波静校正,解决长波长静校正问题和部分短波长静校正问题,然后,利用叠加能量最大静校正技术进一步解决剩余静校正问题,最后,利用非地表一致性剩余时差校正技术,解决速度和射线等误差引起的非地表一致性剩余时差问题.实验结果表明,在以串连的方式应用了三种校正方法之后,在共炮点道集上,折射渡同相轴的线性形态得到了恢复;在动校正后的共中心点道集上,煤层反射的双曲线同相轴被拉平;在叠加剖面上,煤层反射的信噪比得到了改善.  相似文献   

13.
Statics are an effective approach to correct for complex velocity variations in the near surface, but so far, to a large extent, a general and robust automatic static correction method is still lacking. In this paper, we propose a novel two‐phase automatic static correction method, which is capable of handling both primary wave statics (PP statics) and converted‐wave statics (S‐wave statics). Our method is purely data driven, and it aims at maximizing stacking power in the target zone of the stack image. Low‐frequency components of the data are analysed first using an advanced genetic algorithm to estimate seed statics and the time structure for an event of interest, and then the original full‐band data are further aligned via the back‐and‐forth coordinate descent method using the seed statics as initial values and the time structure for event alignment guidance. We apply our new method to two field datasets, i.e., one for 2D PP static correction and the other for 3D S‐wave static correction.  相似文献   

14.
Static correction is a common step in a seismic data proccessing flowchart for land data. Here we propose a new algorithm for automatic short‐period static correction. The algorithm is based on the assumption that seismic events after short‐period static correction should be locally plane nearly everywhere. No other assumptions are made. Therefore the proposed method does not require a preliminary velocity analysis. The algorithm consists in two main parts: evaluation of second spatial differences of trajectories and subsequent regularized integration of these differences. The proposed method proves its robustness and shows results comparable with conventional residual static correction based on improving common‐midpoint stacking. In contrast to the conventional residual static, the proposed algorithm can estimate short‐period statics in complex cases where common‐midpoint stacking fails because of non‐hyperbolic events.  相似文献   

15.
共反射面元(CRS)叠加考虑了反射层的局部特征和第一菲涅耳带内的全部反射,更充分挖掘了多次覆盖数据的潜力.但在地下介质复杂并存在倾斜层时存在反射点分散的情况,从而影响了CRS叠加效果.本文从射线理论出发,在考虑反射层局部特征的情况下,推导了水平地表和起伏地表情况下计算真实反射点分散程度的公式,最终将反射点分散程度定量表达出来.通过对反射点分散程度的控制,从CRS道集中抽取出共反射点(CRP)道集,在CRP道集中而不是在CRS道集中实现叠加,其效果应比传统的CRS叠加效果要好.利用水平地表和起伏地表的模型验证了本文所推导的公式的正确性和有效性.该公式在实际资料处理中的运用尚待进一步研究.  相似文献   

16.
厚风化层覆盖区转换波静校正方法   总被引:3,自引:1,他引:2       下载免费PDF全文
P-SV转换波处理与传统的P-P波处理有很大的不同,如S波静校正、CCP叠加、P-SV速度分析和偏移等,其中最大的难题就是S波静校正问题.S波速度基本不受潜水面的影响,与纵波静校正没有直接相关性,有时横波静校正量能达到纵波静校正量的十倍,用纵波静校正量乘以比例系数来解决横波静校正问题将导致较大误差.同一接收点X和Z分量存在一定的初至时差,该时差代表了P波和S波在低降速带的走时差,可以利用该时差和近地表纵横波速度比信息去除低降速带对横波的影响,得到准确的静校正量.本文利用多分量初至时差推导了较为精确的横波静校正公式,再结合共检波点叠加求取剩余静校正量的方法,形成了完整的转换波静校正配套方法.利用该方法对苏里格气田二维及三维多波地震资料进行了实际处理,数据处理结果证明了该方法的有效性,该方法尤其适用于其他方法难以奏效的风化层较厚地区的横波静校正量求解,该方法也同时考虑了长波长横波静校正问题.  相似文献   

17.
In the application of a conventional common‐reflection‐surface (CRS) stack, it is well‐known that only one optimum stacking operator is determined for each zero‐offset sample to be simulated. As a result, the conflicting dip situations are not taken into account and only the most prominent event contributes to any a particular stack sample. In this paper, we name this phenomenon caused by conflicting dip problems as ‘dip discrimination phenomenon’. This phenomenon is not welcome because it not only leads to the loss of weak reflections and tips of diffractions in the final zero‐offset‐CRS stacked section but also to a deteriorated quality in subsequent migration. The common‐reflection‐surface stack with the output imaging scheme (CRS‐OIS) is a novel technique to implement a CRS stack based on a unified Kirchhoff imaging approach. As far as dealing with conflicting dip problems is concerned, the CRS‐OIS is a better option than a conventional CRS stack. However, we think the CRS‐OIS can do more in this aspect. In this paper, we propose a workflow to handle the dip discrimination phenomenon based on a cascaded implementation of prestack time migration, CRS‐OIS and prestack time demigration. Firstly, a common offset prestack time migration is implemented. Then, a CRS‐OIS is applied to the time‐migrated common offset gather. Afterwards, a prestack time demigration is performed to reconstruct each unmigrated common offset gather with its reflections being greatly enhanced and diffractions being well preserved. Compared with existing techniques dealing with conflicting dip problems, the technique presented in this paper preserves most of the diffractions and accounts for reflections from all possible dips properly. More importantly, both the post‐stacked data set and prestacked data set can be of much better quality after the implementation of the presented scheme. It serves as a promising alternative to other techniques except that it cannot provide the typical CRS wavefield attributes. The numerical tests on a synthetic Marmousi data set and a real 2D marine data set demonstrated its effectiveness and robustness.  相似文献   

18.
常规的转换波静校正的基本思想都是从地震波的运动学特征出发,基于地表一致性假设.在地表条件复杂和地表高程相差较大的地区,它不仅无法解决严重的静校正问题,反而会带来新的畸变.本文基于频率波数域波动方程偏移原理,采用波场延拓方法实现转换波静校正,其关键点在于时间空间域和频率波数域的相对应.文中通过坐标变换将起伏地表转化为新坐标系下的水平地表,把炮点和检波点映射到同一水平面上,然后在新坐标系下推导频率域波动方程延拓公式,接着对下行波P和上行转换波SV分别利用近地表速度向上延拓到基准面,恢复起伏地表到基准面之间的真实波场,最后转换到原始坐标系取出基准面数据完成转换波静校正.通过对模拟和实际数据处理,证明该方法是正确和有效的.  相似文献   

19.
The rough sea surface causes perturbations in the seismic data that can be significant for time‐lapse studies. The perturbations arise because the reflection response of the non‐flat sea perturbs the seismic wavelet. In order to remove these perturbations from the received seismic data, special deconvolution methods can be used, but these methods require, as input, the time varying wave elevation above each hydrophone in the streamer. In addition, the vertical displacement of the streamer itself must also be known at the position of each hydrophone and at all times. This information is not available in conventional seismic acquisition. However, it can be obtained from the hydrophone measurements provided that the hydrophones are recorded individually (not grouped), that the recording bandwidth is extended down to 0.05 Hz and that data are recorded without gaps between the shot records. The sea surface elevation, and also the wave‐induced vertical displacement of the streamer, can be determined from the time‐varying pressure that the sea waves cause in the hydrophone measurements. When this was done experimentally, using a single sensor seismic streamer without a conventional low cut filter, the wave induced pressure variations were easily detected. The inversion of these experimental data gives results for the sea surface elevation that are consistent with the weather and sea state at the time of acquisition. A high tension approximation allows a simplified solution of the equations that does not demand a knowledge of the streamer tension. However, best results at the tail end of the streamer are obtained using the general equation.  相似文献   

20.
静校正问题是地震勘探的关键问题,直接影响地震勘探精度和准确性.实际地震采集过程中,当在相同接收点位置上不同时间内插拔布设了不同的检波器时,对于目前基于地表一致性理论假设的基准面静校正和剩余静校正,以及非地表一致性剩余静校正都不具备适用条件.为解决这一问题,本文提出了基于共姿态道集的静校正方法,将相同接收点位置上不同时间布设的检波点所接收的地震数据抽成不同的共姿态道集,在共姿态道集内实施地表一致性静校正;当某接收点位置上具有若干个共姿态道集时,该接收点位置上可能会存在多个检波点静校量;炮点静校问题仍然采用地表一致性静校正方法解决.该方法解决了同一接收点位置上不同共姿态道集之间的非地表一致性静校正问题,同时也解决了全区的检波点和炮点的地表一致性静校正问题,在实际数据应用效果明显.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号