琼东南海域水合物OBS资料成像处理关键技术

海底地震仪(OBS)资料同时具备纵波和转换横波的优势,但其接收站点少、接收间距大的特点导致其处理方式相对常规海洋地震资料而言具有极大的挑战.对本文主要针对于OBS资料处理的关键技术:利用直达波进行OBS重定位处理,保证后续处理的正确性;利用水检资料和陆检资料进行波场分离,获得上行波场和下行波场;利用下行波进行P波镜像叠前时间偏移处理;对转换横波利用欧拉重构进行重定向处理及极化旋转处理、进行转换横波叠前时间偏移处理.最终获得海底成像清楚、海底以下成像分辨率高、偏移归位准确的OBS纵剖面及相应的转换横波剖面.同时具备的纵、横波剖面利于解释专家对水合物成矿区域进行综合分析.     

三维地震与OBS联合勘探揭示的神狐海域含水合物地层声波速度特征

以三维高分辨地震与海底高频地震仪（OBS）联合勘探数据为基础,获得海底之下沉积层的地震反射成像剖面及多波信息,并以此确定研究区含天然气水合物沉积层的纵、横波速度的变化特征.根据走时反演获得的横波速度与纵波速度对比分析发现,研究区海底之下500 m深度范围内的某些沉积层具有较高的纵横波速度,这一纵波速度升高区域与水合物稳定带对应,而纵波速度下降并且横波速度变化较小的区域,可能与游离气的存在相关.游离气的可能存在与基于这一区域2007年钻探测井结果的普遍认识不完全相符.     

台湾西南海域地震数据处理及天然气水合物地震属性

对台湾西南海域增生楔部位长排列多道地震数据进行地震成像、速度分析、AVO分析、AVO反演处理,获得了天然气水合物多属性地震特征.在偏移剖面上,BSR与海底近似平行,极性与海底相反,穿越沉积层.AVO分析显示,强BSR振幅部位,BSR振幅随偏移距增大而增大.精细速度分析表明强BSR振幅下方存在纵波低速层.对应于强BSR振幅部位,AVO反演的P波、G波为相对高负值区,位于P、G交会图的第三象限,该部位泊松比变化率为负值,横波反射系数接近于零.以上多属性地震特征均预示着该区域可能存在天然气水合物层,且天然气水合物层下方可能存在游离甲烷气层.     

海域天然气水合物的地震研究进展

海域天然气水合物的地震研究是在岩石物性分析成果的指导下,用地震技术进行真假BSR的识别和含水合物沉积层物性预测.海域天然气水合物地震研究的进展主要表现在：物性分析理论模型由简单模型发展到最大程度模拟实际情况的复杂模型;地震勘探方法已由常规的单道、多道地震发展到多频地震、高分辨率二维、三维地震和海底多分量地震;地震资料处理由常规处理发展到突出BSR特征的“三高”和叠前时间偏移处理;从利用速度、振幅结构研究识别天然气水合物发展到AVO、多属性判别、多弹性参数和多物性参数反演识别天然气水合物、预测其物性参数.这些新技术、新方法的应用,加快了海域天然气水合物调查进度,提高了天然气水合物地球物理识别的可靠性.     

低信噪比海底多分量地震数据波场分离方法

以多分量地震观测为基础,联合纵波和转换横波数据能更有效地估计地下介质的弹性和物性参数,提升地质构造成像与油气储层描述的精度.在海底多分量地震数据处理过程中,观测记录的上-下行波分解和P/S波分离可压制水层鸣震以及P与S波之间的串扰,对偏移成像和纵横波速度建模至关重要.但受海底环境、仪器与观测因素共同影响,许多海底多分量地震资料都无法基于现有的海底波场分离方法与流程取得合理的结果.本文以海底声波场与弹性波场分离基本原理为基础,通过对方法流程的修正,摆脱常规流程对中小偏移距直达波信号的依赖性.借助模拟数据实验讨论了波场分离对海底介质参数、噪声的敏感性.结合东海YQ探区海底多分量地震资料上-下行P/S波分离及其叠前深度偏移处理,验证了本文方法流程的可行性.     

基于海底电缆数据的声弹耦合介质多参数最小二乘逆时偏移（英文）

在海洋地震勘探中,海底电缆技术采集的多分量地震数据,涉及到多参数的反演与成像问题,本文针对海底电缆多分量地震数据提出了一种弹性波多参数最小二乘逆时偏移方法。该方法的波场延拓算子为混合方程,即在海水介质中采用声波方程进行波场计算,而海底固体介质的波场由纵横波形分离的矢量弹性波方程得到。在海底界面中采用声弹耦合控制方程将两种类型的方程结合起来。通过推导纵横波形分离的弹性波偏移算子、反偏移算子和梯度公式实现基于纵横波形分离的弹性波最小二乘逆时偏移方法。通过模型试算证明了该方法能够得到高质量的纵波速度和横波速度分量的成像剖面,相比与传统弹性波最小二乘逆时偏移方法,多参数耦合造成的成像串扰噪音得到了很好地压制。     

南海北部东沙海域天然气水合物的初步研究

利用地震、测井与地温资料综合分析了南海北部东沙海域可能存在的天然气水合物的分布特征.研究表明,在东沙海域地震剖面上出现似海底反射层、弱振幅带等天然气水合物分布标志,在声波测井曲线上呈现高速、速度倒转等天然气水合物存在特征.似海底反射层的深度与1144站位,及平均地温梯度资料得出的稳定带厚度较吻合.1144站位与1148站位似海底反射层距海底较深,分别为654m与475m.在1144站位附近,弱振幅带的顶界可能代表含天然气水合物沉积层的顶界,约在450m左右.     

海域天然气水合物勘测的地球物理方法

到目前为止,已经探明的天然气水合物储量要远小于预测的水合物储量.这种现状不同程度地反映出我们对天然气水合物赋存规律认识不足,和对天然气水合物勘探方法认识不足.本文根据作者多次参加水合物地球物理调查国际航次的认识及文献资料,综述了海域天然气水合物勘探方面一些有效的地球物理技术方法,以利于我国海域天然气水合物的勘探工作以及勘探方法的创新.本文指出,地震勘探是目前进行天然气水合物勘探最常用、也是最重要的方法.地震方法主要包括传统的单道、多道地震方法、高分辨地震方法、深拖多道地震探测方法、海底地震仪方法、多道-多分量海底地震电缆方法、海底地震检波器方法等.此外,根据水合物发育区特有的海底地形地貌特征和水体异常特征,根据水合物发育所需要的温度-压力场特征、电磁特征和含水合物地层的剪切模量特征发展的多波束方法、旁扫声纳方法、海底热流探测方法、海底电磁方法以及海底重力测量方法等都在海域天然气水合物勘探中有着很好的效果.     

基于稀疏反演的OBS数据多次波压制方法

自由表面多次波压制是海底地震仪（Ocean Bottom Seismometer,OBS）数据处理和成像中的难点,OBS数据多次波能量强,周期长,严重影响深层一次反射波的处理和成像.不同于常规拖缆观测系统,OBS数据站点一般相隔较远,仅仅利用检波点稀疏的波场信息难以压制OBS数据中的自由表面多次波.本文采用拖缆数据与OBS数据联合,利用稀疏反演估计（Estimation of Primaries and Multiples by Sparse Inversion,EPSI）方法,研究了OBS数据自由表面多次波压制理论,分析了OBS多次波产生的机理,详细推导了拖缆数据与OBS数据联合预测OBS多次波的EPSI方法基本原理.通过利用拖缆数据的信息,实现了OBS检波点稀疏数据多次波的压制问题.EPSI方法通过稀疏反演直接估计一次反射波,避免了SRME（Surface Related Multiple Elimination）方法中自适应相减对有效信号的损害,保真了一次反射有效信号,理论模拟OBS数据验证了方法的有效性.     

井间透射-反射波场分离方法及应用

近十几年来,井间地震在石油勘探领域发挥着重要作用.然而井间地震波极其复杂,给同相轴拾取、层析成像、偏移成像等工作带来了困难.因此,波场分离成为井间地震数据处理的一项重要内容.目前,不少学者研究了井间上下行波分离,纵横波分离等多种方法,但尚未出现相关文献研究透射波和反射波的分离方法.本文根据井间透射波和反射波在不同道集(共炮点集或共接收点集)中视速度极性的差异,首次提出了一种"透射.反射"波场分离方法.该方法首先在不同道集中利用τ-p变换进行视速度分离,得到反射上行波、反射下行波、透射上行波、透射下行波四个波场,然后选择对应的波场叠加,从而得到透射波和反射波.其中,透射波为井间层析成像提供了相对清晰的初至同相轴,尤其是横渡初至同相轴,以便进行到时提取;所得反射波为井中偏移成像提供了可靠的反射波数据,比人工切除透射波方法更为准确有效.本文通过理论模型试验和实际数据处理,证明了该方法的有效性和实用性.此外,在本文基础上,再进行纵横波分离处理可以得到更多的波场分量,为井间成像提供更有效的数据.     

多波高斯束叠前深度偏移速度分析(英文)

多分量地震资料叠前深度偏移技术可以对地下复杂地质构造进行更准确的成像,精确成像的前提是获取准确的纵横波偏移速度。本文采用高斯束偏移方法对多波地震数据进行偏移速度分析,首先分别给出纵波和转换波共偏移距域高斯束叠前深度偏移方法原理,在此基础上抽取纵波和转换波偏移距域共成像点道集;然后根据共成像点道集拉平准则,分别对纵波和横波速度进行更新;当两种波成像深度不一致时,对纵波和转换波成像剖面进行深度匹配,完成高精度的纵横波偏移速度分析。模型数据和实际资料试算表明,该方法是一种有效的多波偏移速度分析方法。     

椭圆展开共反射点叠加方法的应用研究

本文详细介绍了均匀介质条件下椭圆展开共反射点(CRP)叠加原理,并引入双参数(上行波与下行波的速度比和平均速度)来解决非均匀介质条件下的叠加成像,严密论证了所求得的速度是真正的共反射点叠加速度,并结合理论模型计算和地震资料处理证实,利用椭圆展开CRP方法可以对复杂地质剖面求取准确的共反射点叠加速度和正确的零偏移距剖面,得到的成像效果远优于传统共中心点(CMP)方法.     

海底地震仪对海上地震定位精度的提升作用——以南澳外海域M_L2.8级地震为例

由于海底地震台站的缺乏,如何提高海域地震定位的精度是目前地震学家面临的问题之一。南海北部陆缘作为板内地震的活跃带,其微震活动在华南沿海和近海海域频发,海上由于未布设相应的地震台网,因此对近海海域地震尤其是微震的定位精度极低。利用2012年南海东北部海陆地震联测期间两台海底地震仪(Ocean Bottom Seismographs,OBS)所记录到的南澳外海域发生的ML2.8级地震为例,对比分析OBS数据加入前后对地震定位结果的影响。这是首次在南海北部陆缘利用OBS联合陆上地震台网对天然地震进行定位的对比研究,结果说明了近海海域OBS的布设将会对我国沿海地震定位精度的提升、地震活动性的特征以及海域断裂构造的研究具有极其重要的意义。     

海底天然气水合物OBS多分量地震正演模拟（英文）

根据Ecker的水合物沉积物的三种微观模式,计算含水合物沉积层和含游离气沉积物的弹性模量,分析对比了水合物的不同微观模式、不同水合物饱和度以及不同游离气饱和度对沉积物弹性模量的影响;从纵横波分离的弹性波动方程出发,采用交错网格空间有限差分方法模拟地震波在海底天然气水合物沉积地层的传播,得到纵、横波的海底地震（OBS）共接收点道集。数值算例表明,当水合物作为流体的一部分或胶结颗粒骨架时,仅纵波记录上存在BSR;当水合物胶结颗粒接触,纵、横波记录上均存在BSR。并且,OBS会接收到上行纵波和上行横波在海底界面形成的转换波,干扰横波记录上BSR的识别。     

宽频带海底地震仪的地震重定位对马尼拉俯冲板片形态的约束

将宽频带OBS用于海底天然地震长期观测,在国内尚处于实验阶段.2015年在马尼拉俯冲带北段开展了为期6个月的宽频带海底天然地震观测试验.根据回收的1台海底地震仪(OBS04)与国际地震台网的718台陆地地震台站,共记录到7562个P波走时和5002个S波走时数据,利用Hyposat地震定位方法,对马尼拉俯冲带北部(119°E—123°E,19°N—22°N)在2015年8月至2016年2月期间的264个地震进行了重定位.地震重定位的结果及定位误差分析表明,在海域布设的OBS04台站让地震观测的空间分布更为合理,提高了地震定位精度;重定位后的震中分布更为集中,与地质构造吻合良好;浅部的地震活动较为活跃,分布密集,与浅部断层发育有关;重定位后的4条震源深度投影剖面,从不同角度较好地约束了俯冲板片上边界的板片形态,板片倾角在浅部0～30km区间约为10°～22°,随着深度的增加,俯冲板片逐渐变陡,在深度120～180km处倾角约为41°～58°.该项研究为马尼拉俯冲带北段的板片形态提供了重要约束,而且为今后长期天然地震观测提供了重要而宝贵的经验与借鉴.     

最小二乘反演迭代算法在压制海上变深度缆采集数据虚反射中的应用

海水面的虚反射(鬼波)引起海上拖缆采集数据陷波,导致地震记录频带变窄,而近年发展的变深度缆采集技术,具有多样的陷波特征,通过专门的去虚反射处理方法可获得宽频数据.本文基于已有研究成果,将最小二乘反演迭代压制虚反射算法应用于某海上变深度缆宽频处理.基于频率波数域镜像记录生成方法获得镜像炮集记录,并采用最小二乘解从变深度缆原始和镜像炮集记录中提取上行波.针对镜像炮集记录生成受初始速度模型精度的影响,使得某深度缆接收的上行波和下行波之间的实际延迟时间存在误差,采用最小二乘反演迭代算法最优化计算下行波与上行波之间的平均延迟时间和上行波记录,并采用时空数据窗口滑动克服延迟时间随炮检距和目的层深度变化问题.合成数据及某海上实际变深度缆数据处理测试结果表明,该方法能较好地压制变深度缆由海水面产生的虚反射,能达到拓宽地震记录频带目的.     

角度域弹性波Kirchhoff叠前深度偏移速度分析方法

为提高地震成像结果的准确性并真实反映实际地震波场在介质中的传播特性,应该充分利用多分量地震数据的矢量特征进行弹性波成像,其中,最为棘手的问题是纵横波偏移速度场的确定,为此,本文提出了直接利用多分量地震数据进行弹性波角度域偏移速度分析的方法.基于空移成像条件的弹性波Kirchhoff偏移方程提取了弹性波局部偏移距域共成像...     

多波多分量高斯束叠前深度偏移

本文对基于弹性波动理论的多波多分量高斯束偏移进行了完整且详细的分析和公式推导,实现了3D空间多分量(矢量)波场的直接成像.由于当前多数基于弹性波动方程的偏移方法只是假设应力边界条件为自由地表边界条件,这种假设不符合垂直地震剖面(VSP)和海底电缆(OBC)等地震数据.为此本文详细分析了实际应用中常见的三种弹性各向同性介质模型的应力边界条件:自由空间、海底和自由地表模型.在上行传播假设情况下,获得了应力边界条件与位移边界条件的关系式.在此基础上,准确推导了3D多波多分量高斯束波场延拓和偏移成像公式,并在偏移过程中实现了完整的多波型自动分离.由于常规的互相关成像条件不适用于矢量波场成像,本文引用了散度/旋度互相关成像条件.通过约定PS转换波的正向旋转方向解决了3D空间PS成像极性翻转问题.利用2D和3D模型数据偏移成像验证了我们所提出的多波多分量高斯束偏移方法的可行性.     

海域天然气水合物全波形反演研究进展

在海域天然气水合物勘探过程中,全波形反演可以在缺少井资料情况下提供水合物及游离气区域精细速度结构,与AVO分析相比,全波形反演可以有效降低薄层效应影响.概述了全波形反演的研究进展及其在海域天然气水合物参数结构分析中的应用现状;讨论了天然气水合物全波形反演策略,包括基于层状介质假设的全局优化与局部优化相结合反演策略和基于有限差分、有限元正演模拟的多参数结构直接反演策略等;探讨了天然气水合物全波形反演中一些具有代表性的技术特点,包括地震子波估计方式、海底多分量地震数据的应用、多参数同时反演;最后给出天然气水合物全波形反演发展趋势和可能改进建议.随着三维全波形反演算法的不断改进,海域天然气水合物勘探精度将进一步提高,为水合物商业开发提供依据.     

西南海海槽地震资料处理及其似海底反射层特征

对西南海海槽的NT2-2地震测线进行了重新处理,突出了天然气水合物似海底反射层(BSR)的特征.比较了近炮检距剖面、叠加剖面、偏移剖面、瞬时振幅剖面与瞬时相位剖面上BSR的显示特征.指出了BSR在炮道集与共中心点道集上的分布.虽然在近炮检距剖面上能分辨部分BSR,但通过叠加与偏移可以得到BSR较清晰的图像.BSR的高振幅特征在瞬时振幅剖面上最为明显,而在瞬时相位剖面上容易勾画BSR的横向展布.