潮汕坳陷多次波组合压制技术研究

潮汕坳陷中生界发育,是我国海上重要的油气勘探后备区。但该区多次波非常发育,多次波识别和压制一直是该区地震资料处理的关键和难点。合理地选取有效的多次波压制方法是获得高品质地震剖面的关键所在。针对原始资料中多次波异常发育的情况,文章分析了该区多种类型多次波的特点,尝试了零偏移距截距时间-射线参数(τ-p)域预测反褶积、自由界面多次波衰减法(Surface-Related Multiple Elimination,SRME)和双曲线Radon变换,对多次波进行组合压制,克服了仅依靠一种手段难以压制多种类型多次波的局限。从资料处理效果看,τ-p域预测反褶积方法对短周期的鸣震和微曲多次波的压制效果较好,SRME有效地压制自由表面多次波,且在近炮检距范围内压制多次波效果很好,而双曲线Radon变换是衰减长周期的层间多次波和未被充分压制的自由表面多次波很好的选择。组合压制多次波后的剖面质量明显得到改善。     

炮集域压制海水鸣震的τ-p变换法

利用Radon变换压制多次波的传统做法一般都是先进行叠加速度分析,然后在CMP道集上做相应滤波处理.实际资料处理中,带有多次波的速度分析难以准确地提取速度,进而造成多次波压制的效果不够理想.有别于传统的做法,本文提出的方法通过对炮集记录进行同相轴重排,运用变换域自适应切除的τ-p变换(线性Radon变换)滤波,实现了直接在炮集记录中压制海水鸣震.实际资料处理表明,该方法对于海水鸣震多次波压制是一种行之有效的方法.     

三步组合法衰减深水多次波

深水地震资料多次波具有能量强、周期长的特点,降低了资料的信噪比,是处理人员面临的首要问题。多次波衰减效果的好坏,直接影响着最终成像效果及后期储层反演结果的真实性。分析深水多次波的类型、成因和特征,并研究各种多次波衰减技术的优缺点和适用性,最终合理组合多方法、分步衰减多次波是关键。以南海A区深水地震资料为例,通过组合优选基于模型的SRME技术和高分辨率Radon技术进行了实际应用,并就出现的剩余多次波给予分析和改进,提出了“三步法”衰减深水多次波,应用后取得了显著的多次波衰减效果。     

A case study: imaging OBS multiples of South China Sea

The subseafloor structure offshore South China Sea was imaged using first-order water-layer multiples from ocean-bottom seismometer data and the results were compared to conventional imaging using primary reflections. The mirror-imaging method employs a primaries-only reverse time pre-stack depth migration algorithm to image the receiver ghosts. The additional travel path of the multiples through the water layer is accounted for by a simple manipulation of the velocity model and processing datum: the receivers lie not on the sea floor but on a sea surface twice as high as the true water column. Migration results show that the multiple-migrated image provides a much broader illumination of the subsurface than the conventional image using the primaries, especially for the very shallow reflections. The resulting image from mirror imaging has illumination comparable to the vertical incidence surface streamer (single-channel) reflection data.     

SRME技术在澳大利亚Timer Sea地区的应用

表面多次波压制（SRME）技术是目前应用效果较好的多次波衰减技术。在介绍多次波类型和常见压制方法的基础上,引出SRME方法。基于平面波传播模型讨论SRME技术的基本原理,以及针对实际地震数据的迭代预测和减去方法,结合该技术的特点给出实际地震数据应用流程。澳大利亚Timer Sea地区海洋地震资料应用实例表明SRME这种数据驱动的表面多次波压制方法非常实用,尤其对崎岖海底多次波有显著效果。     

Estimating saturation of gas hydrates using conventional 3D seismic data,Gulf of Mexico Joint Industry Project Leg II

We present a methodology for generating pre-drill estimations of hydrate saturations using conventional 3D seismic data. These seismic-based estimates will be compared with well log derived saturations from the subsequently drilled wells of the Gulf of Mexico Gas Hydrate Joint Industry Project Leg II (JIP Leg II) expedition.Predicting saturation of gas hydrates (S_h-seismic) combines pre-stack seismic inversion, rock physics modeling and stratigraphic interpretation. Before the wells were drilled, no nearby sonic and density logs were available to define and calibrate the elastic property trends for the shallow target interval containing the gas hydrate stability zone. Therefore, rock property trends were established by applying principles of rock physics and shallow sediment compaction, constrained by known regional geological parameters. S_h-seismic volumes were generated by inverting pre-stack data to acoustic impedance (PI) and shear impedance (SI) volumes, and then analyzing deviations from modeled impedance trends. To enhance the quality of the inversion, the signal-to-noise ratio of the offset data was maximized by conditioning the seismic prior to inversion. Seismic stratigraphic interpretation plays an important role by identifying the more promising strata and structures for the presence of gas hydrates.The pre-drill S_h-seismic results are compared with saturations calculated from log data, S_h-log, of the wells drilled in the JIP Leg II campaign. Due to weaker seismic reflections, predictions may be less accurate for low concentrations, such as saturations less than 40%, and for thin intervals below the vertical resolution of the seismic data (about 15 m). However, the integrated geophysical workflow is very effective for identifying and quantifying significant hydrate concentrations, making the method a promising prospecting technique.     

Optimum parameter selection in offset-dependent predictive deconvolution: testing on multichannel marine seismic data

Marine Geophysical Research - Predictive deconvolution is an effective way to suppress multiple reflections, especially short path multiples, in seismic data. However, the effectiveness of the...     

浅水多次波衰减技术在多道地震数据处理中的应用

拖缆采集的海上多道地震数据受海水间传播的多次波影响,往往在共炮点道集上发育周期性的强振幅干扰波,这些噪音会掩盖海底以下及强反射层以下地层的面貌,严重影响多道地震数据的成像效果。为了消除海上多道地震数据中存在的海水间传播的多次波对地震数据的干扰,首先分析浅水多次波在海水间多次震荡的产生机制,采用了τ-P域静校正延迟技术对其进行压制,处理的核心是在τ-P域对周期性多次波模型进行预测,再通过自适应相减以去除地震数据中的多次波。实际资料处理结果表明,方法对海水间震荡的浅水多次波具有很好的压制效果,经过浅水多次波去除后的叠加剖面信噪比得到有效提高,剖面的品质得以提升。     

Problems associated with seismic facies analysis of diamicton-dominated,shelf glacigenic sequences

Diamicton-dominated shelf glacigenic sequences are characterized on high-resolution seismic records by a structureless to chaotic acoustic texture, often with abundant hyperbolic (point-source) reflections. Existing depositional models, based primarily on seismic data, have been constructed on the basis that this acoustic signature is indicative of subglacial till. However, borehole data from the UK continental shelf indicate that glaciomarine facies are also present. Thus, while this acoustic character may be indicative of diamicton lithofacies, it should not be taken as a predictive genetic indicator of a specific depositional process.     

深海地震资料叠前去噪方法研究

以琼东南深海三维地震资料为例,全面分析深海地震资料特点,深入分析深海干扰波的类型和特点。针对不同类型的干扰波,采用不同的处理方法,开展针对深海地震资料叠前去噪处理技术研究。通过研究深海地震多次波压制、水鸟干扰压制、涌浪干扰和异常振幅压制等处理方法,总结出针对深海地震资料的叠前去噪技术,为今后的深海地震资料处理工作奠定了技术基础。     

Predicting saturation of gas hydrates using pre-stack seismic data, Gulf of Mexico

A promising method for gas hydrates exploration incorporates pre-stack seismic inversion data, elastic properties modeling, and seismic interpretation to predict saturation of gas hydrates (Sgh). The technology can be modified slightly and used for predicting hydrate concentrations in shallow arctic locations as well. Examples from Gulf of Mexico Walker Ridge (WR) and Green Canyon (GC) protraction areas illustrate how Sgh was derived and used to support the selection of well locations to be drilled for gas hydrates in sand reservoirs by the Chevron-led Joint Industry Project (JIP) Leg II cruise in 2009. Concentrations of hydrates were estimated through the integration of seismic inversion of carefully conditioned pre-stack data, seismic stratigraphic interpretation, and shallow rock property modeling. Rock property trends were established by applying principles of rock physics and shallow sediment compaction, constrained by regional geological knowledge. No nearby sonic or density logs were available to define the elastic property trends in the zone of interest. Sgh volumes were generated by inverting pre-stack data to acoustic and shear impedance (PI and SI) volumes, and then analyzing deviations from modeled impedance trends. In order to enhance the quality of the inversion, we stress the importance of maximizing the signal to noise ratio of the offset data by conditioning seismic angle gathers prior to inversion. Seismic interpretation further plays an important role by identifying false anomalies such as hard, compact strata, which can produce apparent high Sgh values, and by identifying the more promising strata and structures for containing the hydrates. This integrated workflow presents a highly promising methodology, appropriate for the exploration of gas hydrates.     

地震成像道集域的变周期预测反褶积

针对残余多次波的识别与压制问题,提出了1种应用于成像道集域的变周期预测反褶积方法.与传统的预测反褶积方法不同,该方法首先根据所建立的强波阻抗界面的深度-速度模型,在地震偏移剖面上追踪识别出残余多次波,并根据其传播时间确定出相应变化的多次波周期或预测步长,然后在成像道集上对含有残余多次波的成像道进行变周期预测反褶积压制处理,最后再将剔除残余多次波后的成像道集合成为新的地震偏移剖面.实验处理的结果表明:该方法不仅对微屈多次波有效.而且对残余的长周期鸣震多次波也可达到令人满意的压制效果.     

Characteristics of Gas Hydrate and Free Gas Offshore Southwestern Taiwan from a Combined MCS/OBS Data Analysis

In this study, we present the results of the combined analyses of ocean bottom seismometer and multi-channel seismic reflection data collection offshore southwestern Taiwan, with respect to the presence of gas hydrates and free gas within the accretionary wedge sediments. Estimates of the compressional velocities along EW9509-33 seismic reflection profile are obtained by a series of pre-stack depth migrations in a layer stripping streamlined Deregowski loop. Strong BSR is imaged over most of the reflection profile while low velocity zones are imaged below BSR at several locations. Amplitude versus angle analysis that are performed within the pre-stack depth migration processes reveal strong negative P-impedance near the bottom of the hydrate stability zone, commonly underlain by sharp positive P impedance layers associated with negative pseudo-Poisson attribute areas, indicating the presence of free gas below the BSR. Ray tracing of the acoustic arrivals with a model derived from the migration velocities generally fits the vertical and hydrophone records of the four ocean-bottom seismographs (OBS). In order to estimate the Poisson’s ratios in the shallow sediments at the vicinity of the OBSs, we analyze the mode-converted arrivals in the wide-angle horizontal component. P-S mode converted reflections are dominant, while upward P-S transmissions are observed at large offsets. We observe significant compressional velocity and Poisson’s ratio pull-down in the sediment below the BSR likely to bear free gas. When compared to Poisson’s ratio predicted by mechanical models, the values proposed for the OBSs yield rough estimates of gas hydrate saturation in the range of 0–10% in the layers above the BSR and of free gas saturation in the range of 0–2% just below the BSR.     

Evidence of gas hydrate accumulation and its resource estimation in Andaman deep water basin from seismic and well log data

2D and 3D seismic reflection and well log data from Andaman deep water basin are analyzed to investigate geophysical evidence related to gas hydrate accumulation and saturation. Analysis of seismic data reveals the presence of a bottom simulating reflector (BSR) in the area showing all the characteristics of a classical BSR associated with gas hydrate accumulation. Double BSRs are also observed on some seismic sections of area (Area B) that suggest substantial changes in pressure–temperature (P–T) conditions in the past. The manifestation of changes in P–T conditions can also be marked by the varying gas hydrate stability zone thickness (200–650 m) in the area. The 3D seismic data of Area B located in the ponded fill, west of Alcock Rise has been pre-stack depth migrated. A significant velocity inversion across the BSR (1,950–1,650 m/s) has been observed on the velocity model obtained from pre-stack depth migration. The areas with low velocity of the order of 1,450 m/s below the BSR and high amplitudes indicate presence of dissociated or free gas beneath the hydrate layer. The amplitude variation with offset analysis of BSR depicts increase in amplitude with offset, a similar trend as observed for the BSR associated with the gas hydrate accumulations. The presence of gas hydrate shown by logging results from a drilled well for hydrocarbon exploration in Area B, where gas hydrate deposit was predicted from seismic evidence, validate our findings. The base of the hydrate layer derived from the resistivity and acoustic transit-time logs is in agreement with the depth of hydrate layer interpreted from the pre-stack depth migrated seismic section. The resistivity and acoustic transit-time logs indicate 30-m-thick hydrate layer at the depth interval of 1,865–1,895 m with 30 % hydrate saturation. The total hydrate bound gas in Area B is estimated to be 1.8 × 10¹⁰ m³, which is comparable (by volume) to the reserves in major conventional gas fields.     

东沙海域复杂海况下单源单缆长排列地震资料处理*

大量研究表明南海北部东沙海域的中生界地层具有很好的油气成藏条件, 是油气勘探的有利区域。但中生界地层构造成像不清晰限制了该区域油气资源的准确评价。本文利用在该海域采用单源单缆长排列采集的三维地震数据, 提出了对该类数据的优化处理技术方案, 以期获得中生界地层的清晰成像。处理过程主要包括各种噪音干扰和不同类型多次波的压制, 尤其是比较严重的海洋湍流干扰压制, 提高中深层反射信号信噪比; 通过压制气泡和压缩子波提高地震反射的分辨率; 同时利用高密度速度分析进行各向异性双谱非双曲线动校正, 提高成像质量, 使波组特征、振幅特征更为明显。针对密集二维采集数据特点, 通过合理的面元规则化处理, 使得覆盖次数均匀和横向分辨率提高, 实现高精度三维成像。本研究针对噪音干扰较严重的单源单缆长排列加密地震资料, 按照三维地震资料进行数据处理, 处理后明显突出了浅、中、深各层的有效信息, 特别是针对中生界的成像质量改善更为明显, 形成了一套行之有效的处理流程, 为今后邻近区块的勘探和资料处理提供了技术基础, 也为东沙海域的油气勘探与评价提供了科学依据。     

Formation and evolution of the tertiary carbonate reefs in the Madura Strait Basin of Indonesia

Analysis of 2 D seismic data over 4 500 km in length from the Madura Strait Basin in the East Java Sea reveals seismic re?ection characteristics of reefs and associated sedimentary bodies, including asymmetrical or symmetrical dome re?ections, slope progradational re?ections, chaotic re?ections and discontinuous strong re?ections inside the reef, which onlap the ?ank of the reef. It is concluded that the developmental paleo-environment of most reefs is mainly conducive to shallow marine carbonate platform facies and platform margin facies, based on well core data, variations in seismic facies and strata thickness.The formation and evolution of all reefs are primarily in?uenced by the tectonic framework of the Madura Strait Basin. Platform margin reefs are principally controlled by two types of structures: one is a series of E-W trending Paleogene normal faults, and the other is an E-W trending Neogene inversion structures. In addition, wave actions, tidal currents and other ocean currents play an accelerated role in sorting, rounding and redeposition for the accumulation and evolution of reefs. Tertiary reefs in the MSB can be divided into four types: 1) an open platform coral reef of Late Oligocene to Early Miocene, 2) a platform margin coral reef controlled by normal faults in Late Oligocene to Early Miocene, 3) a platform margin Globigerina moundreef controlled by a "hidden" inversion structure in Early Pliocene, and 4) a platform margin Globigerina mound-reef controlled by thrust faults in the early Pliocene. Patterns of the formation and evolution of reefs are also suggested.     

南黄海地震资料邻船干扰特征及压制方法

受过往船只的影响,南黄海地震原始资料普遍存在邻船干扰,使地震数据的品质受到严重影响。在系统总结前人研究成果的基础上,指出邻船干扰在炮集上具有直线型、双曲线型和不规则近双曲线形态,同时具有传播路径相对简单、相对高振幅、宽频带特征。利用将邻船干扰拉平后去除及共偏移距域内f—x预测滤波方法对南黄海二维地震资料进行了测试,结果表明,2种方法均能较好地去除邻船干扰,其中后者更易于实现。     

Mapping gas hydrate and fluid flow indicators and modeling gas hydrate stability zone (GHSZ) in the Ulleung Basin,East (Japan) Sea: Potential linkage between the occurrence of mass failures and gas hydrate dissociation

The Ulleung Basin, East (Japan) Sea, is well-known for the occurrence of submarine slope failures along its entire margins and associated mass-transport deposits (MTDs). Previous studies postulated that gas hydrates which broadly exist in the basin could be related with the failure process. In this study, we identified various features of slope failures on the margins, such as landslide scars, slide/slump bodies, glide planes and MTDs, from a regional multi-channel seismic dataset. Seismic indicators of gas hydrates and associated gas/fluid flow, such as the bottom-simulating reflector (BSR), seismic chimneys, pockmarks, and reflection anomalies, were re-compiled. The gas hydrate occurrence zone (GHOZ) within the slope sediments was defined from the BSR distribution. The BSR is more pronounced along the southwestern slope. Its minimal depth is about 100 m below seafloor (mbsf) at about 300 m below sea-level (mbsl). Gas/fluid flow and seepage structures were present on the seismic data as columnar acoustic-blanking zones varying in width and height from tens to hundreds of meters. They were classified into: (a) buried seismic chimneys (BSC), (b) chimneys with a mound (SCM), and (c) chimneys with a depression/pockmark (SCD) on the seafloor. Reflection anomalies, i.e., enhanced reflections below the BSR and hyperbolic reflections which could indicate the presence of gas, together with pockmarks which are not associated with seismic chimneys, and SCDs are predominant in the western-southwestern margin, while the BSR, BSCs and SCMs are widely distributed in the southern and southwestern margins. Calculation of the present-day gas-hydrate stability zone (GHSZ) shows that the base of the GHSZ (BGHSZ) pinches out at water depths ranging between 180 and 260 mbsl. The occurrence of the uppermost landslide scars which is below about 190 mbsl is close to the range of the GHSZ pinch-out. The depths of the BSR are typically greater than the depths of the BGHSZ on the basin margins which may imply that the GHOZ is not stable. Close correlation between the spatial distribution of landslides, seismic features of free gas, gas/fluid flow and expulsion and the GHSZ may suggest that excess pore-pressure caused by gas hydrate dissociation could have had a role in slope failures.     

Fast static correction methods for high-frequency multichannel marine seismic reflection data: A high-resolution seismic study of channel-levee systems on the Bengal Fan

Very high-frequency marine multichannel seismic reflection data generated by small-volume air- or waterguns allow detailed, high-resolution studies of sedimentary structures of the order of one to few metres wavelength. The high-frequency content, however, requires (1) a very exact knowledge of the source and receiver positions, and (2) the development of data processing methods which take this exact geometry into account. Static corrections are crucial for the quality of very high-frequency stacked data because static shifts caused by variations of the source and streamer depths are of the order of half to one dominant wavelength, so that they can lead to destructive interference during stacking of CDP sorted traces. As common surface-consistent residual static correction methods developed for land seismic data require fixed shot and receiver locations two simple and fast techniques have been developed for marine seismic data with moving sources and receivers to correct such static shifts. The first method – called CDP static correction method – is based on a simultaneous recording of Parasound sediment echosounder and multichannel seismic reflection data. It compares the depth information derived from the first arrivals of both data sets to calculate static correction time shifts for each seismic channel relative to the Parasound water depths. The second method – called average static correction method – utilises the fact that the streamer depth is mainly controlled by bird units, which keep the streamer in a predefined depth at certain increments but do not prevent the streamer from being slightly buoyant in-between. In case of calm weather conditions these streamer bendings mainly contribute to the overall static time shifts, whereas depth variations of the source are negligible. Hence, mean static correction time shifts are calculated for each channel by averaging the depth values determined at each geophone group position for several subsequent shots. Application of both methods to data of a high-resolution seismic survey of channel-levee systems on the Bengal Fan shows that the quality of the stacked section can be improved significantly compared to stacking results achieved without preceding static corrections. The optimised records show sedimentary features in great detail, that are not visible without static corrections. Limitations only result from the sea floor topography. The CDP static correction method generally provides more coherent reflections than the average static correction method but can only be applied in areas with rather flat sea floor, where no diffraction hyperbolae occur. In contrast, the average static correction method can also be used in regions with rough morphology, but the coherency of reflections is slightly reduced compared to the results of the CDP static correction method.     

南海北部揭阳凹陷天然气水合物的地震异常特征分析

大量研究表明南海北部珠江口盆地是天然气水合物发育区,但是该盆地东部揭阳凹陷水合物研究较少。本文利用揭阳凹陷新采集三维地震资料,对该三维地震资料进行成像道集优化和叠前时间偏移处理,得到针对水合物的新处理地震数据体,并通过高精度网格层析反演得到层速度数据体。利用该数据开展叠后约束稀疏脉冲反演,获得含天然气水合物地层波阻抗异常,综合分析反演与地震属性识别水合物。从新处理地震资料看,该区域似海底反射(bottom simulation reflection,BSR)反射呈连续、不连续与地层斜交等特征,BSR发育在一个继承性小型水道上,且下部断裂和气烟囱发育。通过分析BSR特征及BSR上下地层的速度、波阻抗、振幅、频率、相干等属性异常,结合水合物成藏条件,发现了南海北部新的天然气水合物有利富集区,为该区域水合物勘探提供基础。