时延地震互约束频率一致性处理方法

频率是地震波的一个重要属性,地震频率一致性处理是时延地震处理过程中关键的一步.本文提出的时延地震互约束频率一致性处理方法,能够消除两期资料非油藏部分带有的不应该有的频率差异,保证不同时期地震资料频率的差异仅来自油气储层的流体变化,为后续的处理研究奠定了基础.     

时延地震处理中的共反射点道集抽取方法

非重复采集地震资料进行时延地震一致性处理时,首先需要消除两期资料观测系统的差异.在分析同一区块不同时期地震资料观测系统差异的基础上,本文提出共反射点道集抽取方法,大幅度提高了两期资料观测系统的一致性,为后续的时延地震一致性处理打下基础.     

非刚性匹配技术在地震资料一致性处理中的应用

受地震勘探资料的采集施工方法、地表激发接收条件、复杂的地下传播介质等因素的影响会导致同一区块同一地层反射波组的频率、波形不同,静校正不彻底、速度不准会导致剩余时差,动校正拉伸会导致波形畸变,通常在常规一致性处理后的CMP/CRP道集数据仍然存在一些不一致性问题,会影响资料潜力的发挥,资料的一致性处理对处理成果品质的提高起到至关重要的作用,尤其是在多块三维连片处理、时延地震处理、过渡带探区资料处理、叠前反演处理、高分辨率处理、高精度成像处理等方面均对资料的一致性提出了更高的要求.针对该现实问题,介绍了依赖于图像的相似和有局部光滑变形(或位移)场的非刚性匹配处理技术.该技术在常规一致性处理的基础上,首先将叠前数据分选成炮检距数据,优选部分炮检距数据形成模型,然后对每一炮检距数据逐一与模型进行非刚性匹配,实现数据的一致性处理.给出了实际资料非刚性匹配一致性处理的效果,证实了该技术在提高数据的一致性、改善数据处理成果品质、保持分辨率不降低和振幅保真方面具有一定的优势.因此,该技术可适用于时延地震匹配处理、叠前角道集拉平、叠前CMP/CRP道集拉平、振幅随角度变化(AVA)处理、PS波与PP波同相轴匹配处理等.     

非重复采集时移地震正演模拟及可行性分析

本文主要研究前期为常规三维地震勘探和后期为高精度三维地震勘探的这种典型非重复采集时移地震的可行性,主要使用正演模拟和匹配处理相结合的方法进行.本文推导了基于PML边界的三维声波波动方程时空域高阶有限差分正演方程,根据岩石物理流体替代分析方法计算得到时移前后储层的弹性参数,建立了时移地震地球物理模型,设计了常规三维和高精度三维观测系统,进行地震正演模拟,获得了非重复采集时移地震理论模拟数据.利用匹配处理法对非重复采集时移地震正演模拟数据和实际资料进行可行性分析,匹配处理过程主要包括初始数据匹配、一致性处理和互均化处理等.理论模拟数据和实际资料的匹配处理结果表明,通过针对性的分析和处理,利用非重复采集的地震资料进行时移地震研究是可行的,能够得到较好的时移地震响应,在一定程度上反映了油藏的变化情况.     

致密砂岩油气藏地震地质研究关键技术

致密砂岩油气藏在当前油气勘探开中占有重要地位.地震地质综合研究是此类油气藏勘探开发工作的重要组成部分.多数致密砂岩油气藏地表条件复杂,导致地震资料品质较差.油气藏岩性致密、低孔、低渗的地质特征导致地球物理异常微弱.这两个不利条件给地震地质综合研究造成了较大的困难.结合致密砂岩油气藏地震地质研究工作实践,本文对地震资料特殊处理、储层和特殊地质体地震预测与描述的关键技术进展及应用效果进行了阐述.面对黄土塬、山地等复杂地表采集条件和致密砂岩储层等特殊地质体地震地质研究的需要,介绍了高精度静校正、叠前多域去噪、地表一致性振幅补偿、小波变换提高分辨率处理技术及在鄂尔多斯盆地一些探区的应用情况.围绕致密储层等特殊地质体的地震预测与描述需求,着重介绍了致密砂泥岩地层岩石物理分析、煤系等特殊地层的地震属性分析与优化、致密碎屑岩储层反演、叠前地震资料流体检测、致密储层裂缝发育带地震预测和储层综合建模等关键技术在部尔多斯及川西探区实际应用中面临的具体问题、改进对策和实施效果.     

广西地震地下流体观测现状分析

简要回顾了广西地震地下流体观测发展历程,着重分析了广西地震地下流体观测资料质量.包括稳定性与连续性、动态特征、影响因素和信息反映能力,并对广西地震地下流体观测网存在的问题进行了分析,并提出改进措施及建议.     

湖南保靖X区页岩气储层“甜点”区地震资料非线性反演预测方法

针对湖南保靖地区地表覆盖层复杂、地下高陡构造多变、地震勘探资料信噪比低,以及页岩气储层与周围泥岩层波阻抗对比度小等特点,研究敏感测井曲线重构与优化处理、地震资料相控非线性反演和叠后地震属性流体检测的综合储层评价方法。通过对多条测线的地震资料反演和流体检测并与实钻井结果对比,获得了良好的地质效果,初步验证了所研究方法的有效性和实用性。      

新疆沙湾5.2级地震地下水流体前兆异常特征的分析

列举了１９９６年１月９日新疆沙湾５．２级震地下流体前兆观测的最新资料，总结了地震孕育由稳态进行非稳态过程中流体异常的特征。结果表明，地下流体前兆异常以中短期为主，订集中分布在距震中２００ｋｍ范围内，异常形态和幅度复杂多样，并对这次地震的跟踪预测，前兆的可靠性问题进行了讨论。     

利用振幅谱差识别水淹层的变化

储层在长期水驱前后,其中的流体及储层本身会发生变化,导致速度和密度的变化,这一变化为应用四维地震监测油藏的变化提供了依据.本文基于平面波入射水平层状介质,推导出了用振幅谱差（四维振幅谱差）预测油藏变化的公式,并给出数值模拟的结果.理论计算表明,不同程度水驱和不同层位发生水淹其振幅谱差呈现明显的不同,从而在理论上证实了用振幅谱差预测油藏油水变化的可能.这为实际工作中利用振幅谱差预测地下油水变化奠定了理论基础.     

1990年乌苏5.2级地震前地下流体异常

１９９０年乌苏５．２级地震是北天山地区地下流体观测网改造以来，发生在监视区内的第一次中强地震。地震前出现了大面积、多井孔的地下流体异常，显示出了异常在空间上的一些变化。（１）震前异常形态复杂多样。由于不同的观测井其水文地质条件、水文地球化学环境以及观测参数本身的物理化学性质不同，表现出的异常形态也各不相同。另外，异常井孔的受力方式和受力性质的不同也是形成不同异常形态的重要原因。（２）位于震中区的观测井出现异常的比例显著高于外围地区；震中区异常出现时间较晚，外围地区异常出现较早，表现出异常从外围向震中迁移的变化图像。（３）与地震震中处于同一构造带的观测井（泉）容易受到震源应力的牵动和影响，出现异常的可能性大。     

河道砂体宽/窄方位三维观测系统地震物理模型分析

The effect of the wide and narrow azimuth 3D observation systems on seismic imaging precision is becoming a hot area for studies of high precision 3D seismic acquisition methods in recent years. In this paper we utilize 3D physical seismic modeling experiments. A 3D channel sand body physical seismic model is constructed and two acquisition systems are designed with wide azimuth （16 lines） and narrow azimuth （8 lines） to model 3D seismic data acquisition and processing seismic work flows. From analysis of migrated time slice data with high quality and small size, we conclude that when the overlying layers are smooth and lateral velocities have little change, both wide and narrow azimuth observation systems in 3D acquisition can be used for obtaining high precision imaging and equivalent resolution of the channel sand body.     

海上时移地震中多道匹配的观测系统重复性研究

观测系统的重复性是时移地震的关键问题.而对于全面评价,多道重复性度量又是很有意义的.基于归一化均方根差异计算方法,地震多道数据的重复性一般为单道数据重复性的加权均方根.针对海上拖缆监测数据与基线数据多道匹配不确定的问题,采取了定向最佳匹配原则和失配虚拟技术方法.基于地震数据重复性与观测系统重复性的线性相关模型,研究了时移地震观测系统的多道加权均方根重复性特征.结合实际数据研究表明,基于面元的观测系统最佳定向匹配的多道加权均方根重复性基本反映全部观测系统的重复性和覆盖程度.     

Model-based optimization of source locations for 3D acoustic seismic full-waveform inversion

Besides classical imaging techniques, full-waveform inversion is an increasingly popular method to derive elastic subsurface properties from seismic data. High-resolution velocity models can be obtained, and spatial sampling criteria are less strict than for imaging methods, because the entire information content of the seismic waveforms is used. As high operational costs arise from seismic surveys, the acquirable data volume is often limited by economic criteria. By selecting optimal locations for seismic sources, the information content of the data can be maximized, and the number of sources and thus the acquisition costs can be reduced compared with standard acquisition designs. The computation of such optimized designs for large-size 3D inverse problems at affordable computational cost is challenging. By using a sequential receiver-wise optimization strategy, we substantially reduce the computational requirements of the optimization process. We prove the applicability of this method by means of numerical 3D acoustic examples. Optimized source designs for different receiver patterns are computed for a realistic subsurface model, and the value of the designs is evaluated by comparing checkerboard inversion tests with different acquisition designs. Our examples show that inversion results with higher accuracy can be obtained with the optimized designs, regardless of the number of sources, the number of receivers, or the receiver distribution. Larger benefits of the optimized designs are visible when a sparse receiver geometry is used.     

A hexagonal sampling grid for 3D recording and processing of 3D seismic data

3D seismic data are usually recorded and processed on rectangular grids, for which sampling requirements are generally derived from the usual 1D viewpoint. For a 3D data set, the band region (the region of the Fourier space in which the amplitude spectrum is not zero) can be approximated by a domain bounded by two cones. Considering the particular shape of this band region we can use the 3D sampling viewpoint, which leads to weaker sampling requirements than does the 1D viewpoint; i.e. fewer sample points are needed to represent data with the same degree of accuracy. The 3D sampling viewpoint considers regular nonrectangular sampling grids. The recording and processing of 3D seismic data on a hexagonal sampling grid is explored. The acquisition of 3D seismic data on a hexagonal sampling grid is an advantageous economic alternative because it requires 13.4% fewer sample points than a rectangular sampling grid. The hexagonal sampling offers savings in data storage and processing of 3D seismic data. A fast algorithm for 3D discrete spectrum evaluation and trace interpolation in the case of a 3D seismic data set sampled on a hexagonal grid is presented and illustrated by synthetic examples. It is shown that by using this algorithm the hexagonal sampling offers, approximately, the same advantage of saving 13.4% in data storage and computational time for 3D phase-shift migration.     

Towards quantitative evaluation of gas injection using time‐lapse seismic data

Of particular concern in the monitoring of gas injection for the purposes of storage, disposal or improved oil recovery is the exact spatial distribution of the gas volumes in the subsurface. In principle this requirement is addressed by the use of 4D seismic data, although it is recognized that the seismic response still largely provides a qualitative estimate of moved subsurface fluids. Exact quantitative evaluation of fluid distributions and associated saturations remains a challenge to be solved. Here, an attempt has been made to produce mapped quantitative estimates of the gas volume injected into a clastic reservoir. Despite good results using three accurately repeated seismic surveys, time‐delay and amplitude attributes reveal fine‐scale differences though large‐scale agreement in the estimated fluid movement. These differences indicate disparities in the nature of the two attributes themselves, which can be explained by several possible causes. Of most impact are the effects of processing and migration, wave interference effects and noise from non‐repeatability of the seismic surveys. This subject highlights the need for a more careful consideration in 4D acquisition, amplitude processing and use of true amplitude preserving attributes in quantitative interpretation.     

Application of 3D vertical seismic profile multi‐component data to tight gas sands‡

Due to the complicated geophysical character of tight gas sands in the Sulige gasfield of China, conventional surface seismic has faced great challenges in reservoir delineation. In order to improve this situation, a large‐scale 3D‐3C vertical seismic profiling (VSP) survey (more than 15 000 shots) was conducted simultaneously with 3D‐3C surface seismic data acquisition in this area in 2005. This paper presents a case study on the delineation of tight gas sands by use of multi‐component 3D VSP technology. Two imaging volumes (PP compressional wave; PSv converted wave) were generated with 3D‐3C VSP data processing. By comparison, the dominant frequencies of the 3D VSP images were 10–15 Hz higher than that of surface seismic images. Delineation of the tight gas sands is achieved by using the multi‐component information in the VSP data leading to reduce uncertainties in data interpretation. We performed a routine data interpretation on these images and developed a new attribute titled ‘Centroid Frequency Ratio of PSv and PP Waves’ for indication of the tight gas sands. The results demonstrated that the new attribute was sensitive to this type of reservoir. By combining geologic, drilling and log data, a comprehensive evaluation based on the 3D VSP data was conducted and a new well location for drilling was proposed. The major results in this paper tell us that successful application of 3D‐3C VSP technologies are only accomplished through a synthesis of many disciplines. We need detailed analysis to evaluate each step in planning, acquisition, processing and interpretation to achieve our objectives. High resolution, successful processing of multi‐component information, combination of PP and PSv volumes to extract useful attributes, receiver depth information and offset/ azimuth‐dependent anisotropy in the 3D VSP data are the major accomplishments derived from our attention to detail in the above steps.     

Crooked line, rough topography: advancing towards the correct seismic image

Seismic exploration in mountainous areas imposes serious compromises on both acquisition and processing. Access restrictions usually result in profiles that are not straight and are not recorded along the true dip direction (if there is a true dip direction!). Processing constraints often result in very poor approximate corrections for elevations and for deviations from a straight line. Most fundamentally, 2D acquisition and processing assumes that the earth is 2D; this assumption is often seriously violated in mountainous areas. While we cannot efficiently correct 2D seismic data for the effects of a fully 3D subsurface, we can improve the data quality in thrust areas where the assumption of 2D subsurface variation is reasonable. We do this in a series of small steps, which improves the accuracy of several approximations made in processing 2D land data.     

利用深反射地震数据构建的多阶面波频散曲线反演近地表横波速度结构——以跨班公湖—怒江缝合带深反射地震资料为例

深反射地震剖面法为了获取深部结构特征常常采取大的偏移距采集数据.目前公开发表的相关资料中,鲜有利用深反射地震炮集数据获取近地表的结构特征.为此,本文通过正演测试了相关数据处理流程,即利用有限差分正演了起伏地表模型的大偏移距地震单炮弹性波场特征,通过共检波点域面波信号F-K频谱叠加构建新方法,从深反射地震数据集中提取了高品质的多阶面波频散曲线,再利用多阶面波联合反演获得了近地表的结构特征.在前述正演流程基础上,利用跨越班公湖—怒江缝合带的SinoProbe深反射地震剖面中的实际炮集数据,求取了基阶和一阶瑞利波频散曲线,联合反演后得到近地表横波速度结构.该结果与初至波走时反演获取的纵波速度结构具有较好的一致性,且在近地表的浅层分辨率较纵波速度结构特征更高,而更与已有地质认识相吻合.本文提供的相关数据处理流程表明利用深反射地震炮集数据,也能够获取近地表浅层的横波速度结构.     

Seeing below the basalt – offshore Faroes

Exploration in the basalt covered areas of the Faroes offshore has always suffered from poor seismic imaging below the basalt. Long offset 2D and 3D seismic data were acquired and a significant improvement in the seismic image below top basalt has been achieved. Deep towing of the source and receiver cables helped by extending the seismic bandwidth towards lower frequencies. Bubble‐tuned rather than conventional peak‐tuned source arrays gave little, if any, incremental benefit. The improvement in the imaging comes primarily from the approach to processing the data. High frequencies (dominantly noise) are filtered out of the data early in the processing to concentrate on the low frequency data. Careful multiple removal is important with several passes of demultiple being applied to the data using both Surface‐Related Multiple Elimination (SRME) and Radon techniques. Velocity analysis is performed as an iterative process taking into account the geological model. Reprocessing legacy 2D surveys, acquired with wide‐ranging parameters, using these processing techniques improved these datasets significantly, indicating that sub‐basalt imaging seems to be more sensitive to processing than to the choice of acquisition parameters.     

Intra-survey reservoir fluctuations – implications for quantitative 4D seismic analysis

During the time taken for seismic data to be acquired, reservoir pressure may fluctuate as a consequence of field production and operational procedures and fluid fronts may move significantly. These variations prevent accurate quantitative measurement of the reservoir change using 4D seismic data. Modelling studies on the Norne field simulation model using acquisition data from ocean-bottom seismometer and towed streamer systems indicate that the pre-stack intra-survey reservoir fluctuations are important and cannot be neglected. Similarly, the time-lapse seismic image in the post-stack domain does not represent a difference between two states of the reservoir at a unique base and monitor time, but is a mixed version of reality that depends on the sequence and timing of seismic shooting. The outcome is a lack of accuracy in the measurement of reservoir changes using the resulting processed and stacked 4D seismic data. Even for perfect spatial repeatability between surveys, a spatially variant noise floor is still anticipated to remain. For our particular North Sea acquisition data, we find that towed streamer data are more affected than the ocean-bottom seismometer data. We think that this may be typical for towed streamers due to their restricted aperture compared to ocean-bottom seismometer acquisitions, even for a favourable time sequence of shooting and spatial repeatability. Importantly, the pressure signals on the near and far offset stacks commonly used in quantitative 4D seismic inversion are found to be inconsistent due to the acquisition timestamp. Saturation changes at the boundaries of fluid fronts appear to show a similar inconsistency across sub-stacks. We recommend that 4D data are shot in a consistent manner to optimize aerial time coverage, and that additionally, the timestamp of the acquisition should be used to optimize pre-stack quantitative reservoir analysis.