复杂地表条件下共反射面元(CRS)叠加方法研究

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

共反射面元叠加的应用理论--从共反射点到共反射面元

共反射面元（CRS）叠加是目前认为最好的生成零炮检距剖面的方式. 共反射面元 意指地下某一反射点邻近的一个反射弧段,该弧段在时空域内的走时响应称为CRS叠加面,该 叠加面可视为反射弧段上各共反射点（CRP）的时空域内走时响应的组合. 在一般的共反射 点走时关系基础上,引入两种特征波——Normal波和Normal Incidence Point波,就可以在 傍轴近似假设下,将CRP走时关系推广到反射点邻近的各反射点,将这些反射点的CRP走时关系 加以组合就得到了关于该反射点的共反射面元的走时关系. 考察从共反射点（CRP）到共反 射面元(CRS)的过渡,这一过程提供了CRS叠加的应用理论基础.     

胜利探区基于最优孔径的共反射面元(CRS)叠加的成功应用(英文)

由于CRS叠加考虑了反射层的局部特征和第一菲涅耳带内的全部反射,从而更充分地利用了多次覆盖反射数据的信息。就目前的地震资料处理技术而言,它是最佳的零偏移距成像方式。本论文利用改进型的参数优化技术,得到高质量的CRS运动学参数剖面,并利用参数剖面计算出叠加孔径,实现了基于最优孔径的CRS叠加,使CRS参数的用途得到了充分利用。模型数据和实际资料的试算表明,基于最优孔径的CRS叠加的成像剖面与传统CRS叠加剖面相比,有着较高的信噪比和同相轴的连续性。     

部分共反射面元叠加技术在海上数据中的应用

本文对共反射面元(CRS)叠加方法做改进,利用得到的波场参数来提高叠前地震资料的质量.利用CRS波场参数做部分CRS叠加,对菲涅尔带内的多个相邻CMP道集做倾角、曲率等校正后合并为一个道集即CRS超道集,可以补齐缺失地震道,实现叠前数据规则化,并提高信噪比.从而使得叠前道集中的同相轴尤其是来自深层的反射有更好的连续性,有利于识别和追踪.提高质量后的叠前道集可用于后续的速度分析、叠加、偏移等常规处理中,效果好于原始CMP道集.模型和实际数据的计算结果验证了该方法的正确性和有效性.该方法在低信噪比资料的处理中将会有广阔的应用前景.     

共反射面道集偏移速度建模

共反射面（CRS）叠加是一种与宏观速度模型无关，仅依赖于近地表速度的地震成像方法．其通过地震三参数的优化实现地震成像．本文推导了基于CRS叠加得出的优化三参数与偏移速度之间的解析关系，提出了在CRS道集通过优化三参数实现速度估计的CRS道集偏移速度建模方法．模型试算表明，这种速度建模方法效率较高，速度分析精度取决于优化三参数的精度，适于较复杂地质体的速度建模．      

基于CRP不规则二维地震采集方法研究（英文）

中国西部复杂山地是地震勘探的难点区域,地表和地下都复杂,从而导致地表激发困难,地震资料干扰严重、信噪比低。地下构造形态复杂,导致反射点分散严重,前人尝试过弯线采集,但很少考虑地下实际构造特征。本文章提出基于CRP的不规则采集设计方法,使地表避开难以激发的区域,地下考虑构造特征及CRP点分布位置,优化地表接收线位置,通过地表有目的布置接收点位置,最大程度的保证地下CRP点不分散,从而提高地震数据的信噪比.并以四川盆地某工地的实际资料通过模型数据验证了方法的合理性与可行性。表明该方法能够解决复杂地区地震采集存在的部分难题,为复杂地区的地震勘探采集提供一个新的思路。     

叠前三维共反射面叠加在实际地震资料处理中的应用

针对复杂地表地质条件下地震资料覆盖次数不均、部分数据缺失而导致的地震资料品质下降、陡倾角地区成像质量变差等问题,文中采用了叠前三维共反射面叠加技术.共反射面叠加是一种与宏观速度无关且考虑了反射点曲率的地震成像方法.它利用共反射点道集一个邻域内(菲涅尔带)道之间的相关性,并将相干区域内道集的能量相加来增强地震反射信号和绕射波能量,并借助于相邻CMP道集数据形成CRS超道集,使得面元内的覆盖次数更加均匀;同时利用超道集的高覆盖次数来压制噪声,最终得到高信噪比道集数据.实际地震资料处理表明,该方法可大幅度地提高地震资料信噪比和分辨率,增强地震同相轴的连续性,此外还可为叠前偏移成像处理、AVO属性分析、叠前地震反演等提供高质量的输入数据.     

二维菲涅耳带共反射面元叠加方法研究

大量研究证明CRS叠加能提高地震勘探的信噪比,但是能否提高分辨率的关键在于如何确定CRS叠加孔径.本文详细探讨了地震波反射过程中菲涅耳带的特征,认为起伏地形下菲涅耳带可以采用椭圆予以近似,在此基础上提出了一种通过菲涅耳带来确定CRS叠加孔径的方法,并应用于泌阳凹陷陡坡带的地震剖面.结果表明,由于菲涅耳带确定的叠加范围使地震信号的能量达到最佳,相对于CMP叠加,菲涅耳带CRS叠加同时提高了地震资料的信噪比和分辨率,特别显示了中深部较弱的地震信息,而常规的CRS叠加则只在于提高资料的信噪比和改善浅部的地震信息.     

共反射面元叠加的应用实践

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

起伏地表煤田地震资料静校正

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

Multifocusing revisited ‐ inhomogeneous media and curved interfaces*

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.     

Handling dip discrimination phenomenon in common‐reflection‐surface stack via combination of output‐imaging‐scheme and migration/demigration

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.     

Directional statics in common reflection surface stack for rugged surface topography

Seismic data acquired along rugged topographic surfaces present well‐known problems in seismic imaging. In conventional seismic data processing, datum statics are approximated by the surface consistence assumption, which states that all seismic rays travel vertically in the top layer. Hence, the datum static for each single trace is constant. In case this assumption does not apply, non‐constant statics are required. The common reflection surface (CRS) stack for rugged surface topography provides the capability to deal with this non‐vertical static issue. It handles the surface elevation as a coordinate component and treats the elevation variation in the sense of directional datuming. In this paper I apply the CRS stack method to a synthetic data set that simulates the acquisition along an irregular surface topography. After the CRS stack, by means of the wavefield attributes, a simple algorithm for redatuming the CRS stack section to an arbitrarily chosen planar surface is performed. The redatumed section simulates a stack section whose acquisition surface is the chosen planar surface.     

共反射面元走时曲面计算方法

共反射面元走时曲面计算是共反射面元叠加的关键.常规共反射面元叠加必须通过相干搜索和优化确定共反射面元叠加公式中的三个属性参数（二维），从而确定共反射面元走时曲面，该类算法具有三点不足：①相干搜索及优化法计算量大；②共反射面元叠加公式仅适用小炮检距；③波前曲率半径取负号且较小时，共反射面元叠加公式基本不适用.为此，本文提出了利用共反射点射线追踪拟合共反射面元走时曲面的计算方法.模型计算证明该方法比传统共反射面元叠加走时曲面计算精度高,适用性强.     

Common reflection surface stack using dip decomposition for rugged surface topography

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.     

The common focal point and common reflection surface methodologies: subsets or complements?

The common focal point (CFP) method and the common reflection surface (CRS) stack method are compared. The CRS method is a fast, highly automated procedure that provides high S/N ratio simulation of zero‐offset (ZO) images by combining, per image point, the reflection energy of an arc segment that is tangential to the reflector. It uses smooth parametrized two‐way stacking operators, based on a data‐driven triplet of attributes in 2D (eight parameters in 3D). As a spin‐off, the attributes can be used for several applications, such as the determination of the geometrical spreading factor, multiple prediction, and tomographic inversion into a smooth background velocity model. The CFP method aims at decomposing two‐way seismic reflection data into two full‐aperture one‐way propagation operators. By applying an iterative updating procedure in a half‐migrated domain, it provides non‐smooth focusing operators for prestack imaging using only the energy from one focal point at the reflector. The data‐driven operators inhibit all propagation effects of the overburden. The CFP method provides several spin‐offs, amongst which is the CFP matrix related to one focal point, which displays the reflection amplitudes as measured at the surface for each source–receiver pair. The CFP matrix can be used to determine the specular reflection source–receiver pairs and the Fresnel zone at the surface for reflection in one single focal point. Other spin‐offs are the prediction of internal multiples, the determination of reflectivity effects, velocity‐independent redatuming and tomographic inversion to obtain a velocity–depth model. The CFP method is less fast and less automated than the CRS method. From a pointwise comparison of features it is concluded that one method is not a subset of the other, but that both methods can be regarded as being to some extent complementary.