东辛油田地震勘探历程——一个复杂断块构造的三维地震解释实例

东辛油田是一个典型的断块复杂构造。在这个油田上,地震精查的方法在搞清断裂系统方面,取得了较好的成效。本文简述了在该油田早期所采用的一套三维地震解释（立体归位）方法,其中有些方法对当前搞三维地震勘探,及断层面的地质解释方面都还有着参考价值。近年来,又用数字地震方法,在此浅层的复杂构造的下面,发现了比较简单的深层构造,在几张附图中展示了深浅层构造之间的关系。     

汶川地震发震构造带过WSFD-2钻井深反射地震数据3-D处理技术研究

2008年5月12日汶川发生里氏8.0级地震的发震机制与龙门山断裂带的构造特征紧密相关,应用反射地震探测方法精细探测发震构造内部的结构特征,对研究汶川地震的发震机理有重要意义.5.12地震后在地表错断最剧烈的区域之一虹口乡白庙村穿过发震构造和深钻科研井WSFD-1和WSFD-2布置了兼顾浅中深层信息的反射地震探测剖面,该区域断裂带硬岩出露,地形地貌复杂多变,起伏剧烈,只能采取弯线地震测线布设方式采集二维反射地震数据,而且弯曲度大,造成CDP点分散严重.区域内构造复杂,褶皱逆推构造发育,地层和构造倾角大,采用常规二维地震数据处理方法进行叠加处理时,易将不同地层的反射信息叠加到同一反射层.弯线叠加剖面上侧面波,混波干扰严重,会显著地降低原本低信噪比数据的分辨率,或剖面可靠性低,容易在二维剖面中造成解释陷阱.本文利用弯线地震采集的三维特性,用拟三维地震叠加技术处理汶川地震科学钻井附近横跨北川-映秀断裂带的二维弯线地震数据,弥补常规二维弯线地震数据处理技术的不足,将不同地层的反射信息分离归位到不同的三维叠加剖面上,解决二维弯线地震数据处理时混波干扰严重的难题.通过理论分析和实验选取合适的共中心点面元,获得了高分辨率的三维叠加剖面.相比于二维弯线叠加剖面,三维叠加剖面切片成像更真实,剔除了不同地层反射信息混叠的影响,能得到更准确的断点信息,并可以获得沿断层走向横向的信息,显著提高了构造解释的可靠性和精度.应用拟三维地震数据处理方法处理龙门山断裂带的二维弯线地震数据,获取高分辨率的构造信息,有利于断裂带内汶川地震发震机制和龙门山隆起机制的解释.     

三维叠前深度偏移速度分析及蒙特卡洛自动层速度拾取

是否能够正确地建立深度域三维速度模型是三维叠前深度偏移成败的关键 .本文根据Deregowski循环 ,利用叠前深度域地震成像对速度模型变化的敏感性 ,采用偏移迭代逐次逼近最佳成像速度 ,研究开发了一套快捷有效的三维叠前深度偏移深度域速度模型建立技术 .借鉴时间域CDP(共深度点 )道集上常规叠加速度分析的策略 ,在深度域CRP(共反射点 )道集上 ,提出剩余慢度平方谱的概念并建立相应的实现技术 .导出深度域中均方根速度与层速度之间的关系 ;按照串级偏移原理确定偏移循环过程中初始速度、剩余速度及修改后速度之间的关系 ;采用蒙特卡洛非线性优化算法实现从剩余慢度平方谱中自动拾取层速度 ,讨论了其地质速度约束条件和蒙特卡洛非线性优化的收敛准则 ,使得所拾取的层速度模型具有合理的地质意义并获得最佳偏移成像效果 .SEG EAGE理论模型数值试算验证了方法的有效性 ,在海拉尔盆地霍多莫尔工区 ,5 8km2 三维资料的速度模型建立并获得满意的三维叠前深度偏移成像 .     

CDP RAY MODELING IN THE PRESENCE OF 3-D PLANE ISOVELOCITY LAYERS OF VARYING DIP AND STRIKE*

Various Dix-type formulae are derived, which are useful to approximate travel time functions that can be observed while modeling the common depth point (CDP) technique for 3-D isovelocity layers of varying dip and strike. All formulae can be used to compute interval velocities and recover the depth model from surface measurements. They are established by making use of the concept of wavefront curvature. Many similarities with known formulae valid for the 2-D plane isovelocity layer case exist.     

EFFICIENT 3-D MIGRATION IN TWO STEPS*

The addition of a third dimension to seismic data acquisition offers significant improvement of our reflection images of geologic structures, but potentially at a substantial increase in data processing cost. Here we study an approach to migration of three-dimensional (3-D) stacked data that can considerably reduce the burden of computation and data manipulation. Cascading two-dimensional (2-D) migrations in orthogonal lateral directions, we obtain an efficient 3-D scheme that is exact for homogeneous media. As applied to examples of synthetic data from inhomogeneous media, the scheme introduces errors well below those attributable, in practice, to uncertainties in migration velocity. Application of this two-step approach, using both Kirchhoff-summation and finite-difference algorithms, to a 3-D seismic land survey over an area of known hydrocarbon traps additionally proves the importance of 3-D migration in stratigraphic imaging.     

利用三维高斯射线束成像进行地震定位

常规的地震定位方法通常需要拾取地震记录的初至,当初至不明显或被较高水平的噪声淹没时精度较低.本文采用基于三维高斯射线束的偏移成像方法对震源进行定位,较好地解决了该问题.通过三维高斯射线束对台站记录进行偏移归位,并将各台站成像结果的交点作为地震能量释放的中心位置;当各台站成像结果不能交于一点时,采用三维空间高斯滤波方法可实现震源位置的自动获取.提出的变网格计算方案极大地减少了计算量,显著地提高了成像精度和计算效率.利用首都圈地震台网数据,对涿鹿、滦县以及房山三个地震事件进行试算,结果表明:基于变网格三维高斯束偏移成像的地震定位方法自动化程度很高,而且具有较好的抗噪能力,特别适合处理低信噪比资料的地震定位问题.     

Multidimensional GPR array processing using Kirchhoff migration

We compare the ability of several practical ground-penetrating radar (GPR) array processing methods to improve signal-to-noise ratio (SNR), increase depth of signal penetration, and suppress out-of-plane arrivals for data with SNR of roughly 1. The methods include two-dimensional (2-D) monostatic, three-dimensional (3-D) monostatic, and 3-D bistatic Kirchhoff migration. The migration algorithm is modified to include the radiation pattern for interfacial dipoles. Results are discussed for synthetic and field data. The synthetic data model includes spatially coherent noise sources that yield nonstationary signal statistics like those observed in high noise GPR settings. Array results from the model data clearly indicate that resolution and noise suppression performance increases as array dimensionality increases. Using 50-MHz array data collected on a temperate glacier (Gulkana Glacier, AK), we compare 2-D and 3-D monostatic migration results. The data have low SNR and contain reflections from a complex, steeply dipping bed. We demonstrate that the glacier bed can only be accurately localized with the 3-D array. In addition, we show that the 3-D array increases SNR (relative to a 2-D array) by a factor of three.     

三维地震勘探数据处理的问题及其解决方法

三维地震勘探在十年内取得了显著的地质效果和经济效益。三维地震数据处理技术也得到了迅速的发展。但是,与二维相比,三维地震数据处理技术有许多不完善之处。本文对其中几个主要问题进行分析并提出解决方法。这些问题是三维叠加速度分析、三维静校正和叠前三维偏移等处理方法。     

波动方程三维吸收边界分解与拟合深度偏移方法

在偏移问题中引入吸收边界条件，既可以消除由人工边界激发的虚假反射，从而提高剖面质量。又可以减少计算工作量．本文讨论了三维吸收边界条件方程，提出了求解具有吸收边界条件的三维波动方程偏移定解问题的分解与拟合方法。理论分析与合成记录及野外实际地震资料处理结果表明，本文方法为一有效的三维吸收边界深度偏移方法。     

表驱三维角度域Kirchhoff叠前时间偏移成像方法

针对传统Kirchhoff叠前时间偏移方法的一些不足,以及振幅随入射角、方位角变化(AVA/AVAZ)分析的需要,本文提出一种基于射线理论的三维叠前时间偏移角度域成像方法.它通过横向均匀介质中稳健的射线追踪建立单程波走时和传播角度的数值表,然后在此基础上估算反射波双程走时以及在界面处传播的方位角和入射角,最后基于脉冲响应叠加原理获得三维构造图像和方位\|角度域共成像点道集.与传统方法不同之处在于,上述过程均考虑了地震波在垂向变速介质中的射线弯曲效应和三维传播特征,有利于准确提取随入射角和方位角变化的振幅和时差信息.理论模型合成数据和实际地震资料测试结果展示了方法的优越性与实用性.     

基于信号子空间分解的三维地震资料高分辨率处理方法

提出了一种新的基于信号子空间分解技术的3 D地震资料高分辨率处理方法.利用信号子空间分解技术,不仅可以分离3 D地震资料中的信号和噪声,而且可以进一步根据地震同相轴的倾角不同,将混合信号进行分解得到单个同相轴信号.通过丢弃噪声子空间,只对不同信号子空间重构的信号利用谱白化技术进行高分辨率处理,然后累加所有处理结果,从而达到既提高地震资料的分辨率,又提高地震资料信噪比的目的.对合成资料和实际地震资料的处理结果表明,此法具有好的应用前景.     

3-D Prestack Kirchhoff Migration of the ISO89-3D Data Set

— This paper presents an overview of the results obtained from a 3-D prestack depth migration of the ISO89-3D data set. The algorithm is implemented as a Kirchhoff-type migration, in which the migrated image is generated by weighted summation along diffraction surfaces through the shot record section. The diffraction surfaces are computed by a 3-D finite difference solution of the eikonal equation. A 3-D macro-velocity model derived mainly from wide-angle tomographic inversion served as input for the travel-time calculations. The results of the migration are presented as slices through a volume covering an area of 21 km × 21 km in the horizontal and 15 km in the vertical direction, centered around the KTB drill hole. In these slices the continuation of the Franconian Lineament or SE1 reflector, respectively, can be identified over most of the survey area as a northeast dipping reflector plane. Its signature appears partly curved and discontinuous and with different strength of reflection down to a maximum depth of 9 km. About 5 km to the south-southeast of the KTB drill hole the uppermost top reflection of the Erbendorf body (EB) can be recognized at approximately the same depth. The slices clearly show its complicated internal structure consisting of several apparently separated reflective parts. Moreover, the geometry and the shape of a few other subsurface structures are described.     

2-D Versus 3-D Magnetotelluric Data Interpretation

In recent years, the number of publications dealing with the mathematical and physical 3-D aspects of the magnetotelluric method has increased drastically. However, field experiments on a grid are often impractical and surveys are frequently restricted to single or widely separated profiles. So, in many cases we find ourselves with the following question: is the applicability of the 2-D hypothesis valid to extract geoelectric and geological information from real 3-D environments? The aim of this paper is to explore a few instructive but general situations to understand the basics of a 2-D interpretation of 3-D magnetotelluric data and to determine which data subset (TE-mode or TM-mode) is best for obtaining the electrical conductivity distribution of the subsurface using 2-D techniques. A review of the mathematical and physical fundamentals of the electromagnetic fields generated by a simple 3-D structure allows us to prioritise the choice of modes in a 2-D interpretation of responses influenced by 3-D structures. This analysis is corroborated by numerical results from synthetic models and by real data acquired by other authors. One important result of this analysis is that the mode most unaffected by 3-D effects depends on the position of the 3-D structure with respect to the regional 2-D strike direction. When the 3-D body is normal to the regional strike, the TE-mode is affected mainly by galvanic effects, while the TM-mode is affected by galvanic and inductive effects. In this case, a 2-D interpretation of the TM-mode is prone to error. When the 3-D body is parallel to the regional 2-D strike the TE-mode is affected by galvanic and inductive effects and the TM-mode is affected mainly by galvanic effects, making it more suitable for 2-D interpretation. In general, a wise 2-D interpretation of 3-D magnetotelluric data can be a guide to a reasonable geological interpretation.     

2-D and 3-D processing and interpretation of multi-fold ground penetrating radar data: a case history from an archaeological site

A 2.5-D and 3-D multi-fold GPR survey was carried out in the Archaeological Park of Aquileia (northern Italy). The primary objective of the study was the identification of targets of potential archaeological interest in an area designated by local archaeological authorities. The second geophysical objective was to test 2-D and 3-D multi-fold methods and to study localised targets of unknown shape and dimensions in hostile soil conditions. Several portions of the acquisition grid were processed in common offset (CO), common shot (CSG) and common mid point (CMP) geometry. An 8×8 m area was studied with orthogonal CMPs thus achieving a 3-D subsurface coverage with azimuthal range limited to two normal components. Coherent noise components were identified in the pre-stack domain and removed by means of FK filtering of CMP records. Stack velocities were obtained from conventional velocity analysis and azimuthal velocity analysis of 3-D pre-stack gathers. Two major discontinuities were identified in the area of study. The deeper one most probably coincides with the paleosol at the base of the layer associated with activities of man in the area in the last 2500 years. This interpretation is in agreement with the results obtained from nearby cores and excavations. The shallow discontinuity is observed in a part of the investigated area and it shows local interruptions with a linear distribution on the grid. Such interruptions may correspond to buried targets of archaeological interest. The prominent enhancement of the subsurface images obtained by means of multi-fold techniques, compared with the relatively poor quality of the conventional single-fold georadar sections, indicates that multi-fold methods are well suited for the application to high resolution studies in archaeology.     

SEG/EAGE盐丘和推覆体模型的波动方程三维叠前深度偏移成像

三维波动方程叠前深度偏移是复杂介质中进行构造成像、弹性参数反演的重要环节．由于其技术实现不仅涉及波场延拓理论的创新,而且需要大规模计算,因而研究难度较大．本文以实验效果的取得为目的,完整地实现了SEG／EAEG盐丘和推覆体模型的三维波动方程辛几何算法的叠前深度偏移成像计算．文中详细考察了所研制的波动方程三维叠前深度偏移软件系统及其对复杂地质构造的成像能力,具体包括：1)对于盐丘模型,文中讨论了成像参数的选择、地震子波对成像精度的影响、完成二维及三维叠前深度偏移的比较;2)对推覆体模型,文中进行了脉冲响应测试;3)由两个模型的成像结果可见本文的波动方程三维叠前深度偏移软件系统已具有适应强速度横向变化、复杂构造的成像能力。     

IMPROVEMENT OF SIGNAL-TO-NOISE RATIO IN RESISTIVITY PROFILES1

The cross-correlation technique for increasing the anomaly-to-noise ratio is applied to the interpretation of resistivity profiles. To verify the method theoretically, resistivity profiles above a body having the shape of a parallelepiped were simulated with correlated and uncorrelated noise. Seven different electrode geometries are discussed. As a practical test, we considered profiles of a geoelectrical survey to locate tombs at a site of archaeological interest. When the shape and dimensions of the anomaly can be foreseen, the cross-correlation method can be applied and it gives a substantial improvement in the signal-to-noise ratio.     

拉东投影法三维叠前深度偏移

对地下地质构造进行正确成像是地震勘探的最终目的,由于三维地震资料采集不可能都沿垂直构造走向的方向进行,为地震资料的三维处理带来了许多困难. 本文将三维叠后拉东投影偏移思想应用于三维叠前处理中,提出了三维叠前投影偏移算法. 利用拉东投影变换的原理,将整个三维叠前数据体投影到一系列各方向的径向线上,各方位角的构造都包含在其中某条或多条径向剖面上. 投影完成后,形成一系列的独立的二维叠前测线,可采用二维叠前深度偏移成像方法来实现各径向线的叠前偏移,当各径向剖面偏移完成后,在时间切片上进行反投影,从而最终形成三维叠前深度偏移结果. 实际应用表明,用本方法进行三维叠前深度偏移时,深度偏移剖面对横向分辨率有所提高,对陡地层和小断层的成像效果有所改善.     

稀疏采样下陡角度构造的波动方程深度偏移成像

针对稀疏采样,特别是crossline方向稀疏采样的三维地震数据,提出了一个基于波动方程方法的陡角度构造深度偏移成像方法.文中针对一类混合域（空间和波数）的波动方程偏移方法展开研究.通过建立反假频的单程波算子和通过（形式上）填充空道重建理想采样的地震数据,有效地实现了陡倾角地层的准确成像.这一算法策略解决了波动方程偏移方法在实际应用中的一个关键问题,从而更好地发挥了波动方程偏移方法在复杂构造成像上的优势.二维理论数据（Marmousi模型）和三维实际地震资料成像结果表明本文方法是有效的.这一方法对现行的三维地震数据处理有重要的现实意义.     

三维高阶深度偏移方程及其数值求解方法

从三维非均匀介质中的波动方程出发，利用拟微分算子理论，Ｐａｄｅ逼近方法及因式分解技巧，获得了非均匀介质的三维高阶深度偏移方程，相应地提出了逐次低阶方程方法、低阶方程组方法及分裂方法等３种求解方法．与二维情形不同，以上每一种方法在数值求解时均存在由测线坐标ｙ的出现而带来的困难．为了克服这一困难，我们提出了差分算子分解方法，避免了近年来人们竞相研究的ｘ－ｙ方向微分算子分裂带来的分裂误差，保持了应有的相容性，解决了这一令人烦恼的问题．     

STUDY AND APPLICATIONS OF SPATIAL DIRECTIONAL FILTERING IN THREE-COMPONENT RECORDINGS1

Three-component recordings permit the construction of particle trajectories. These three-dimensional pictures of particle motion show successive predominant directions of polarization and allow wave modes with distinct polarization directions to be recognized. A polarization selection, called ‘spatial directional filtering’, can be accomplished by several methods; four techniques are described. The application of these polarization filters to a noise shot, offset VSPs and CDP stack are also presented. This type of filtering is shown to cancel waves with undue polarization and to enhance the signal-to-noise ratio.