来自海底高速层径向波的理论地震图研究

本文利用各向异性反射率技术计算理论地震图,提出海底高速薄层会产生沿高速层水平传播的波（简称径向波）,这种波在水层中作为P波,在固液界面激发下行横波,该均匀横波以临界角入射高速薄层,在层内作为超临界角的非均匀横波水平传播,再以临界角转换为上行传播的均匀横波,最终在固液界面上行透射转换为水层中P波.高速薄层传播的径向波不同于界面折射波,也不同于具有频散的面波和通道波.理论地震图的研究表明,径向波具有线性时距,能与海底强反射具有同等振幅水平;径向波有其振幅、时距位置和斜率这些观测记录参数,分别对应高速层的厚度、深度和近似的横波速度;径向波可以克服折射波解释中遇到的振幅强弱和高速层速度等困难.径向波可作为探测海底高速薄层的有力工具,对于研究高速层屏蔽、海底反射类型的多样性和相应的资料处理解释有重要意义.     

地震反射法中薄煤层分辨能力的研究

本文以非线性地震层状构造理论为基础,对薄层分辨率进行了研究,研究中采用理论计算和合成地震记录,并用超声地震模型试验加以验证。结果表明,薄层可检测性准则,应是它的复合反射波的振幅不小于围岩中较强反射波的振幅。由Farr提出的可检测薄层的极限厚度为λ/12（λ为薄层中地震波波长）,是在薄层顶、底界面反射系数很小时的近似值。指出,当薄层顶、底界面反射系数较大时,如煤层,其可检测性高于Farr定义,且随界面反射系数的增大或频率的提高而提高;适当地提高地震反射波的频率,就能得到厚2-3m薄煤层的大振幅、高信噪比的反射波。上述结论均用煤田地震勘探实例证实。     

薄层时频特征的正演模拟

基于正演模拟的薄层时频特征响应分析是薄层定性识别与厚度定量预测的基础与关键.本文基于波场延拓理论,利用深度域相移法对不同厚度薄层进行了正演模拟.采用广义S变换对零偏移距地震道的反射复合波进行瞬时频谱分析,发现当薄层顶底反射系数极性相反时,复合波的干涉作用具有升频降幅作用,而极性相同时具有升幅降频作用;理论推导和实验分析均证明主峰值频率与薄层双程旅行时间厚度存在定量解析关系,这为薄层厚度定量预测提供了重要的技术手段;由于峰值频率对反射系数大小与极性变化不敏感,因此,峰值频率方法与时域振幅方法相比更具稳定性.     

煤层中流体地震可探测性的模拟分析

作为一种典型的强阻抗差低阻抗薄层,煤层中孔隙含流体时是否会引起地震反射产生明显的异常是回答地震检测流体是否可行的根本.为此,本文针对强阻抗差薄层模型,基于Biot双相介质理论,通过弹性波有限差分法数值模拟,与各向同性单相介质假设的煤层反射对比,探讨了反射复合波受煤层孔隙度及流体性质变化的影响程度.模拟分析发现：由于薄层孔隙度和孔隙流体属性的变化在Biot理论中表现为纵波速度的变化,PP波反射AVO（Amplitude Versus Offset,振幅随偏移距变化）特征对薄层是否含流体相对敏感;综合使用PP与PS波对比有利于薄层中流体的预测;孔隙度一定时,PP波反射振幅随着含气饱和度的增加而增大;受薄层调谐作用的影响,孔隙和流体变化对煤层反射的频谱特征影响不大,近似于单相介质时的情况.     

各向异性衰减薄层地震响应特征研究

薄储层的叠前地震响应特征研究,特别是针对具有速度各向异性的含流体薄层,对储层描述具有十分重要的意义.文中基于波动方程数值模拟方法,正演得到各向同性弹性、各向同性衰减、速度各向异性、各向异性衰减模型的地震波场,并对比分析了四种模型的纵波(PP)和转换横波(PS)地震反射特征.研究结果表明:在衰减介质背景下,引入各向异性,PP和PS波的反射波振幅较弹性介质均减弱,且衰减因素对薄层振幅的影响强于各向异性.同时,VTI各向异性衰减在单频PS振幅曲线上表现出强差异性,而HTI各向异性衰减则会影响PP和PS波单频振幅曲线的极值点幅值和位置,通过分析单频振幅曲线的极值点振幅和极值点位置对各向异性衰减薄层预测有指导作用,尤其对平时较难分辨的VTI各向异性衰减薄层,单频分析方法的优势更明显.     

天然地震PS转换波动力学特征的初步研究

本文根据合成理论地震图的广义射线理论,并应用自行编写的FORTRAN计算程序,研究了单一界面、薄层和迭层界面上形成的PS转换波的动力学特征,讨论了沉积盖层及上部地壳介质吸收对转换波的影响。     

三维倾斜界面P波共反射点道集NMO速度

为了实现三维倾斜界面真正的共反射点叠加,获得更加精确的P波速度模型,进行可靠的AVO/AVA分析,形成与三维多分量PS转换波共转换点道集配套的技术,在常速介质共反射点精确走时的基础上,推导并给出三维倾斜界面P波共反射点道集近似双曲时距关系及NMO速度;认识到近似时距在偏移距与深度比小于2时具有较高的拟合精度;指出P波共反射点道集NMO速度是P波速度、界面倾角、界面倾向和测线方位的函数,且随测线方位变化而呈现椭圆特征.与P波共中心点道集NMO速度不同.P波共反射点道集NMO速度小于介质速度.且随界面倾角增大而减小;椭圆长轴方向为倾斜界面走向方位,短轴方向为倾斜界面倾向方位.     

薄互层弹性波反演面临的困境

对于薄互层地震传播特征的研究我们已有了充分的认识,这是一种完全不同于现有的基于厚层单阻抗差界面地震反射理论的、具有频率依赖性的调谐反射.因而也使得对具有N层结构的弹性介质薄互层反演时,现有有限带宽的基于空间-频率域采样的地震观测难以规避多解性问题.长波长假设下的薄互层VTI近似似乎可以为我们提供一种待反演弹性参数只有5个的可行性方案,使得我们可以将薄互层背景下任一目标薄层的反演看作VTI厚层介质背景下单薄层的反演.但是,我们首先面临薄互层VTI近似的误差和适用性条件等问题;其次,既使对于单薄层厚度的反演,目前地震领域也未能提供一种通用有效的解决方案,实际地震勘探中薄层厚度的预测依然是一个悬而未决的难题.因此在此基础上再进行单薄层的物性参数、裂缝隙和流体的预测则更叠加了难以预知的多解性.显然单纯P波已难以应对,弹性波的介入成为必然;而且单薄层含裂缝与流体的预测更需要横波分裂和双相介质理论的支撑.可以预见的是,除了纵、横波速度,薄互层反射系数的频率依赖性是我们试图给出相对可靠的、单薄层物性精细刻画必须使用的一种重要地震属性.     

粘弹各向异性反射体模型方位地震AVO模拟及分析（英文）

基于岩石物理模型和广义各向异性Zoeppritz方法在频率域计算裂缝型反射体模型反射波方位地震AVO响应。反射体模型为粘弹各向异性、有限厚度的地质体,其地震反射波形序列包含如下动力学信息,即分界面处介质的波阻抗和非弹性差异、反射体内部波的各向异性传播、在传播路径上的频散与衰减,以及来自顶底界面的反射波的调谐与干涉等。计算表明,速度频散和衰减增顶界面反射波大入射角反射时的振幅,而减弱底界面反射振幅。对于固定入射角的方位地震方位地震响应,PP波反射特征表现为随方位角的增加反射波形序列延续时间变长,而PSV和PSH转换类型反射波的方位各向异性变化特征稳定且受储层厚度影响较小,表现为PSV波反射振幅随方位角增加而增加,PSH波在0°和90°方位无反射能量,在45°方位反射振幅最强。     

等厚薄互层时频特征的正演模拟

本文基于波动理论,利用深度域相移法对不同互层数、不同单层厚度的等厚薄互层进行正演模拟.采用广义S变换分析零偏移距地震道反射复合波的瞬时频谱.研究发现当单层厚度大于3/16波长时,等厚薄互层时域特征为中、高频等幅振荡,瞬时频谱为与单层厚度成反比的单一峰值梳状函数;当单层厚度小于3/16波长时,震源子波的有限带宽对瞬时频谱的高频频率峰有压制作用,导致等厚薄互层中部的时域振幅和频域最大幅度均为低值,薄互层的结构特征不易分辨,说明宽高频的震源子波是提高薄互层分辨能力的关键.此外,理论推导和实验分析均证明:当地层厚度大于1/8波长时,其峰值频率与薄互层单层厚度存在定量解析关系,这为等厚薄互层单层厚度定量预测提供了重要的技术手段.     

平均入射角道集PP波与PS波联合反演

在界面两侧地层的弹性参数弱反差的假设难以成立的情况下,本文提出用平均入射角道集进行PP波与PS波的联合反演.首先,在PP波与PS波AVA(amplitude versus angle,振幅随入射角变化)道集的基础上,分别选择小入射角范围与大入射角范围的AVA道集进行局部加权叠加,以获得由两个角度组成的平均入射角道集,并作为后续反演的输入数据.然后,再通过最小二乘原理建立了PP波与PS波联合反演目标函数,推导了模型修改量的向量公式,建立了平均入射角道集联合反演的流程.模型数据与实际数据的测试结果表明:在信噪比较低、地层弹性参数反差较大、层厚较薄的情况下,该反演方法的精度在很大程度上超过了基于近似反射系数的反演方法,为复杂油气藏勘探提供了新的思路.     

Sensitivity of reflection seismic data to oil-column height in high-porosity sandstones

Tuning is the effect of interference between the reflections from the top and bottom of a thin layer on the amplitude of the composite reflection. For a homogeneous sandstone reservoir containing an oil column overlying brine, interference between the reflection from the top reservoir and the oil/water contact is a function of the height of the oil column. If the properties of the sandstone do not vary across the oil/water contact, the SS, PS and SP reflection coefficients from the oil/water contact are small in comparison to the PP reflection coefficient. This allows analytic expressions for the effective PP and PS reflection coefficients from the reservoir to be derived that include all P‐wave multiples within the oil column. For a given source/receiver offset, the component of the wavevector inside the oil column normal to the interface is larger for the PPPP reflection than for the PPPS reflection, due to the asymmetry in the raypath for the PPPS reflection. The PPPS reflection is therefore useful for determining oil‐column heights larger than that discriminated by the PPPP reflection, especially when used at wider offsets. A convenient classification of the AVO response of hydrocarbon‐bearing sandstone reservoirs overlain by shale is the scheme of Rutherford and Williams. Class 1 sands have higher acoustic impedance for normal incidence than the overlying shale, Class 2 sands have nearly the same acoustic impedance as the shale and Class 3 sands have lower acoustic impedance. Synthetic shot gathers calculated for these three classes as a function of oil‐column height show that a combination of the PPPP and the PPPS amplitudes can be plotted as a tuning trajectory, which can be used to determine the oil‐column height. This method is most sensitive for reservoirs that belong to AVO classes 1 and 2, and therefore may be useful in AVO analysis of Class 1 and 2 reservoirs where the traditional AVO indicators (developed for Class 3 reservoirs) do not work very well. These results demonstrate the usefulness of shear waves recorded in the marine environment at wide offsets.     

REFLECTION-WINDOW MAPPING OF SHALLOW BEDROCK*

A field study was undertaken to evaluate the effectiveness of the high frequency seismic reflection technique for mapping of shallow and irregular bedrock. Bedrock reflections were obtained using a hammer source with both in-line and common offset field layouts. The recording equipment included 12-channel enhancement seismographs, 28 Hz vertical geo-phones and a microcomputer. The latter increased the overall versatility of the seismic system. Field sites for this study are typical of the geological settings of the tin mining areas of Malaysia. The topographical ‘lows’ of the irregular bedrock control the localization of tin ore. The subsurface geology consists of a thin low velocity layer (± 300 m/s) overlying the compact overburden (± 1700 m/s) which in turn lies on bedrock. This paper discusses various criteria for designing an optimum window for obtaining usable reflections between the first arrival and the leading edge of the ground roll cone. Detailed mapping of the overburden and the bedrock interface by the reflection method can be useful in delineating areas for exploratory drilling and for optimum planning of mining operations.     

层状各向异性介质转换波克希霍夫叠前时间偏移

在克希霍夫叠前时间偏移处理中,地震波走时的计算方法是决定大偏移距地震资料成像品质的重要因素.在常规的三维转换波各向异性叠前时间偏移公式中,走时的计算是基于等效单层各向异性介质的非双曲线方法.用这种方法处理的成像道集,在偏移/深度比超过一定阈值后,成像道集中的反射同相轴将出现过偏现象,这种偏移不平的同相轴将影响偏移叠加的最佳响应,使得偏移成像波组呈低频化特征,最终降低三维转换波偏移成像质量.我们采用层状介质的走时计算方法代替常规算法,并且利用了常规方法的转换波各向异性偏移速度模型.基于层状介质的算法能够提高大偏移距转换波走时计算精度,克服中浅地层大偏移距远道成像道集中反射同相轴逐渐上翘的问题.两个地区的三维转换波资料处理结果证实,基于层状各向异性介质的转换波克希霍夫叠前时间偏移方法,明显改善了反射成像剖面的连续性和分辨率,提高成像剖面构造的可解释性.     

南海东北部东沙海域海水层反射特征分析

本文基于对南海东北部东沙海域近期采集的多道反射地震资料进行重新处理获得新的地震海洋学数据,分析了该海域内孤立波/内孤立波包、沙丘上方和陡坎附近特殊反射结构特征,从而提供了新的海水层与海底相互作用依据.研究结果表明,除之前已发表文章中地震海洋学资料显示存在的第一模态内孤立波/波包和沙丘上方常见的反射样式-披毛状发射外,地震海洋学资料上还发现了第二模态内孤立波、陡坎上方的上抬型波动反射结构样式.在新的地震海洋学数据中,第一模态内孤立波振幅均小于50 m,宽度上都小于5 km,单个内孤立波的最大振幅约为45 m.内孤立波包的内孤立波振幅都相对较小,均小于40 m,并且与之前不同的是,彼此之间振幅相差不大,没有明显的排列规律.此次地震海洋学数据记录到的第二模态内孤立波,形态较为完整,上层和下层反射的振幅相差不大,在30 m左右;中间层大约在水深130 m位置处,垂向结构的整体大小大于200 m.沙丘上方反射结构普遍存在弱反射层,可能是湍流边界层,并且存在特殊反射样式-披毛状反射.但并不是沙丘上方都存在披毛状反射样式,本文分析它出现在地震海洋学资料上可能是受测线与沙丘走向之间夹角的影响.陡坎区域的水层反射结构则表现为上抬型波动,并常常伴随着同相轴连续性的变化.该波动的大小及反射同相轴的连续性可能取决于陡坎的高度/坡度及水层动力的强度,新数据中出现的一个上抬型波动,高度达20~30 m,它的附近水层还存在一个形态不完整的内孤立波.陡坎附近的水层反射也常常出现弱反射带和小的波动.     

Compositional variation and the origin of deep crustal reflections

Identifying the sources of crustal reflections is essential for deriving more geological information from deep crustal reflection profiles. Theoretical and model studies help place constraints on the role of compositional variation in producing deep crustal reflections. Analysis of laboratory-derived velocities and densities of rock types typical for the continental crust indicates that reflection coefficients are generally small, but significantly, 17% of the possible reflection coefficients have magnitudes between 0.1 and 0.2. Comparison between reflection coefficients derived from well logs and reflections observed in associated seismic profiles indicates that constructive interference associated with geological layering is at least as important as the magnitude of the reflection coefficients for producing detectable reflections. Constructive interference can increase reflection amplitude by two to three times but is limited to a relatively narrow range of layer thickness. For a typical 10–40 Hz seismic wavelet and typical crustal velocities of about 6 km/s, constructive interference occurs for layer thickness ranging between about 35 and 80 m. Layers thinner than 35 m interfere destructively. If reflections result from compositional variation, seismic models of hypothetical and observed geologic relations provide analogs for interpreting complex reflection patterns observed in deep crustal reflection profiles. Such models show reflection patterns similar to those observed in the reflection profiles. The models indicate that the reflections could originate in the complexly deformed and intruded terranes that are common in the crystalline crust and it may not be necessary to appeal to unobserved phenomena such as special lamellae or fluid-filled fractures to explain the reflections.     

Some SH-wave seismic reflections from depths of less than three metres

Shallow SH-wave reflections are far from routine, although their study can provide insights into important properties of near-surface materials that cannot be inferred from P-wave data alone. Difficulties in separating SH-wave reflections from Love waves are generally considered the major obstacle to progress in shallow SH-wave seismic reflection. This may be the case in surveys undertaken at great depths, but it is not necessarily true for reflection data gathered at shallow and ultra-shallow depths. This paper shows that when SH-wave data possess wavelengths greater than the thickness of the superficial layer, Love waves are not greatly dispersed. In this case, misinterpretation between parts of reflection hyperbolae and waveguide arrivals is sufficiently limited. In a one-layer model earth, which well approximates typical situations of the near-surface underground, the most energetic modes (the lowermost modes) of the dispersed surface waves have a dominant frequency band that falls below the wavelet spectrum of the shallow reflections; therefore, they can be filtered out in the frequency domain. Higher modes, although their spectral content overlaps that of the reflections, exhibit small amplitudes on seismograms and leave strong reflections unaffected.We present field examples from three different sites where we were able to obtain ultra-shallow reflections (< 3 m) in unconsolidated sediments. The high level of resolution (vertical resolution up to 15 cm) suggests that SH-wave reflection imaging has the potential to complement other high-resolution techniques, such as P-wave reflection and ground-penetrating radar (GPR) imaging, allowing a better and more complete characterization of the near-surface environments.     

Seismic Rheological Model and Reflection Coefficients of the Brittle–Ductile Transition

It is well established that the upper—cooler—part of the crust is brittle, while deeper zones present ductile behaviour. In some cases, this brittle–ductile transition is a single seismic reflector with an associated reflection coefficient. We first develop a stress–strain relation including the effects of crust anisotropy, seismic attenuation and ductility in which deformation takes place by shear plastic flow. Viscoelastic anisotropy is based on the eigenstrain model and the Zener and Burgers mechanical models are used to model the effects of seismic attenuation, velocity dispersion, and steady-state creep flow, respectively. The stiffness components of the brittle and ductile media depend on stress and temperature through the shear viscosity, which is obtained by the Arrhenius equation and the octahedral stress criterion. The P- and S-wave velocities decrease as depth and temperature increase due to the geothermal gradient, an effect which is more pronounced for shear waves. We then obtain the reflection and transmission coefficients of a single brittle–ductile interface and of a ductile thin layer. The PP scattering coefficient has a Brewster angle (a sign change) in both cases, and there is substantial PS conversion at intermediate angles. The PP coefficient is sensitive to the layer thickness, unlike the SS coefficient. Thick layers have a well-defined Brewster angle and show higher reflection amplitudes. Finally, we compute synthetic seismograms in a homogeneous medium as a function of temperature.     

Velocity/interface model building in a thrust belt by tomographic inversion of global offset seismic data

Between September and November 1999, two test seismic lines were recorded in the southern Apennine region of southern Italy using the global offset technique, which involves the acquisition of a wide offset range using two simultaneously active seismic spreads. One consisted of a symmetrical spread moving along the line, with geophone arrays every 30 m and a maximum offset of 3.6 km. The other one consisted of fixed geophone arrays every 90 m with a maximum offset of 18 km. This experimental acquisition project was carried out as part of the enhanced seismic in thrust belt (ESIT) research project, funded by the European Union, Enterprise Oil and Eni-Agip. An iterative and interactive tomographic inversion of refraction/reflection arrivals was carried out on the data from line ESIT700 to produce a velocity/interface model in depth, which used all the available offsets. The tomographic models allowed the reconstruction of layer interface geometries and interval velocities for the target carbonate platform (Apula) and the overburden sequence. The value of this technique is highlighted by the fact that the standard approach, based on near-vertical reflection seismic and a conventional processing flow, produced poor seismic images in both stack and migrated sections.