横波源三维地震资料矢量横波四分量旋转和快慢波分离技术
时间：2023-01-26 16:05:03 来源：柠檬阅读网 本文已影响 人 
公 亭 王兆磊 罗文山 王永生 岳媛媛 顾小弟
(①东方地球物理公司研究院，河北涿州 072750；
②青海油田勘探事业部，甘肃敦煌 736202；
③东方地球物理公司物探技术研发中心，河北涿州 072750)
弹性介质中能够传播两种不同类型的体波，一种是质点振动方向与传播方向一致的纵波(压缩波)，另一种是质点振动方向与传播方向垂直的横波(剪切波)[1-2]。油气储层的横波和纵波响应特征不同：一方面，横波传播只与岩石骨架有关，对流体不敏感，对存在气云区的特殊油藏，能更准确地对地下构造成像；
另一方面综合利用纵波和横波信息可以提高储层反演和预测的精度[3-4]。基于上述优势，从二十世纪三十年代开始，前苏联科学家就开展横波勘探方法研究，美国、法国、西德也相继开展了大量横波勘探方法的试验。二十世纪七、八十年代后，受震源的限制，横波勘探逐步被成本更低的转换波勘探取代，针对横波勘探的研究几乎停滞。但由于转换波勘探射线路径不对称，上、下行波类型不一致[5-6]，给资料处理、解释带来许多困难，2010年后也开始归于沉寂[7]。
近年来，东方地球物理公司(BGP)通过不断地研究和探索，成功研发了大吨位横波可控震源EV-56S，采集了横波源三维地震资料，横波勘探取得重大技术突破。横波勘探不同于纵波，多分量检波器记录的地震波场与横波震源激发方向、地下介质的裂隙方向直接相关，有着明显的矢量特性，而以往横波处理技术研究局限于二维横波处理或者建立在三维转换波地震资料处理基础上，对横波三维空间上复杂波场的矢量关系考虑不足，无法应用于三维横波地震资料，因此迫切需要形成一套完整的矢量处理技术。
在柴达木盆地三湖地区开展的可控震源纵、横波联合激发的三维多波矢量采集，即在同一位置，纵波震源、x方向横波震源(振动方向与接收排列平行)、y方向横波震源(振动方向与接收排列垂直)各激发一次，三分量数字检波器接收，获得了三维九分量(3D9C)地震数据。
本文仅对横波震源激发、水平分量接收到的四个分量数据开展针对横波勘探多方向性的矢量处理技术研究。
图1 四分量数据定义示意图
为方便后续阐述，在野外观测系统坐标系中，定义平行于检波线方向为x方向，垂直于检波线方向为y方向，如用SxRy表示震源沿x方向水平激发、三分量检波器沿y方向接收的地震道；
在室内处理R、T坐标系中，定义平行于炮检点连线方向为径向R，垂直于连线方向为切向T，如用SRRT表示震源径向投影能量、检波点切向接收能量的地震道。
对于三维横波源地震资料，在各向同性水平层状介质假设条件下，如果横波震源激发方向与穿过炮点和接收点的垂直面即不平行也不垂直，则检波器两个水平分量上均会接收到能量，即在SxRx、SxRy、SyRx和SyRy四个水平分量上都会接收到横波能量。另一方面，受地表条件影响，实际野外采集时，震源车行进方向无法严格遵照施工设计中规定的检波线方向，不同的震源车车头朝向导致接收到的横波相位存在差异，当震源车朝向相反时，子波相位差180°，极性相反。基于以上两方面原因，任意水平分量未经矢量处理直接叠加均无法获得较好的成像效果。
因此，三维横波地震资料处理应首先将四个初始分量数据旋转到与炮、检点连线平行或正交的R、T方向，以寻求R、T方向能量最大化。借鉴横波分裂分析与校正公式[8-9]，四分量旋转可表示为
(1)
式中θ为炮检点连线与x方向的夹角。
如图2所示，使用原始四个水平分量数据进行叠加，成像效果较差。通过四分量旋转后，SRRR和STRT分量的叠加能量得以恢复(图3a和图3b)，但由于实际介质的各向异性，SRRT和STRR分量的叠加剖面仍有有效能量残留(图3c和图3d)，将在后续快、慢波分离处理中得到进一步解决。
图2 四分量旋转前叠加剖面
图3 四分量旋转后叠加剖面
当横波在含有裂缝的各向异性介质中传播时，就会明显分裂为两个偏振互相正交而速度不同的剪切波[10-12]，一个是偏振方向平行于裂缝的快横波，另一个是偏振方向垂直于裂缝的慢横波[13]。
2.1 裂缝方位角求取
图4 旋转互相关法求取的裂缝方位角(a)、Alford旋转扫描法求取的裂缝方位角(b)及对应的SRRR分量叠加剖面(c)
2.2 快慢波分离技术
假设地下含裂缝地层和预处理得到横波地震记录满足：①选定分析时窗内的目标地层只发育单一方向的裂缝；
②地下裂缝地层可等效于HTI介质，即裂缝垂直于地层。
快、慢波分离就是根据横波在各向异性介质中传播特性，应用传播矩阵旋转法，将混合了快、慢横波能量的R、T分量转换为快、慢横波记录
(2)
式中：β为裂缝方向与炮检连线方向的夹角；
SS1RS1分量表示分离后得到的快横波；
SS2RS2分量表示分离后得到的慢横波；
SS1RS2和SS2RS1分量能量在快、慢波分离彻底的情况下为0。
为了验证上述结论，设计了一个三层水平层状模型，第一层、第三层是各向同性介质，第二层是裂缝方位角为60°的HTI介质，具体参数如表1所示。
表1 模型参数
应用反射系数法合成第2个界面的横波反射数据，再将R、T方向四个分量按10°增量划分扇区进行叠加，得到第2界面的横波反射方位角道集(图5a～图5d)。受方位各向异性的影响，SRRR和STRT分量同相轴呈现周期性波浪状，而SRRT和STRR分量每隔90°则会出现极性反转现象。根据SRRR或STRT分量出现最大时差位置可推知，裂缝方位角β=60°。将裂缝方位角代入式(2)，旋转得到快(SS1RS1分量)、慢(SS2RS2分量)波分离结果(图5e、图5f)。分离后的SS1RS2分量(图5g)和SS2RS1分量(图5h)能量为0，说明利用该方法能够比较彻底地分离快、慢横波。
图5 单层HTI介质快、慢波分离结果
图6为实际地震数据SRRR分量和分离的快横波分量的速度谱及道集，可以看出，分离后道集中不同炮检距同相轴错断的异常现象得到有效消除，连续性增强，道集质量明显提高，速度谱能量也更为聚焦。图7、图8分别是快、慢横波分离前、后的叠加剖面和叠前时间偏移剖面对比，从中可以看到，相较于SRRR和STRT分量，快横波和慢横波数据的波阻关系清楚，消除了同一套地层横向上能量强弱不一致的异常现象，同相轴更光滑，层位可解释性明显提高，且快、慢横波波组特征基本一致，同一套层位仅存在时延，随深度增加而增大。
图6 SRRR分量(a)和快横波分量(b)的速度谱(左)及道集(右)
图7 快、慢波分离前、后叠加剖面对比
柴达木盆地三湖地区生物天然气藏发育。但受气云影响，纵波地震剖面成像效果较差，构造变形，无法落实气云区构造形态特征，严重制约该区后续勘探开发进程。为此，BGP在研究区部署了纵、横波联合激发的三维多波地震勘探项目。针对横波资料明显的矢量特性，初步形成了以四分量旋转和快慢波分离技术为核心的横波源三维地震资料矢量横波配套处理技术流程(图9)。通过该套技术流程的应用，获得了高品质的横波地震剖面。
图9 矢量横波四分量处理技术流程
研究区东侧位于A号气田，从图10 的纵、横波时间域处理成果可以看出，纵波受气云的影响，存在明显的低频下拉现象，成像质量极差，而横波在含气异常区成像优势明显，较好地恢复了A号气田区的构造形态。图11是根据纵、横波联合解释方案目的层K9所在位置平均时间提取的时间切片，图12是提取的目的层K13的沿层相干切片，可以看出，相较于纵波，横波除在含气异常区构造成像方面的优势外，对于微小断裂及地质异常体的刻画能力明显增强。
图10 纵波(a)、快横波(b)叠前时间偏移剖面对比
图11 纵波(a)、快横波(b)叠前时间偏移切片对比
图12 纵波(a)、慢横波(b)叠前时间偏移相干切片对比
从图13利用纵、横波成果获得的叠前波阻抗反演结果可见，在相同域(横波标定到纵波域)中，纵波仅能识别研究区厚度约为8.6m的薄储层，而横波能够识别的厚度约为4.6m，相较于纵波，横波分辨率提高了近一倍，这为落实该区域构造背景、助力开发阶段薄储层预测奠定了资料基础。
图13 纵波(a)、慢横波(b)叠前反演剖面对比黑色曲线为伽马测井曲线
(1)在柴达木盆地三湖地区，成功应用了横波源三维地震资料矢量横波配套处理技术，获得了可靠的快、慢横波成果资料。
(2)在研究区，相比于纵波，横波资料不仅能够有效刻画气云区的构造形态，而且分辨率更高。
(3)四分量旋转和快慢波分离技术的联合应用，基本解决了研究区横波勘探的处理难题。但目前仅能应用于单组裂缝发育地区，当实际地质条件更复杂、在多层系发育多组裂缝时，如何逐层开展矢量处理尚需进一步研究。
猜你喜欢 纵波横波震源 花岗岩物理参数与纵波波速的关系分析矿产与地质(2022年2期)2022-09-06基于横波分裂方法的海南地幔柱研究北京大学学报（自然科学版）(2022年2期)2022-04-08横波技术在工程物探中的应用分析建材发展导向(2021年19期)2021-12-06Pusher端震源管理系统在超高效混叠采集模式下的应用*石油管材与仪器(2020年2期)2020-05-11《光的偏振》探究指导书的设计与实现物理教学探讨(2018年11期)2018-01-15扬眉一顾，妖娆横波处传奇故事(破茧成蝶)(2015年6期)2015-02-28横波一顾，傲杀人间万户侯火花(2015年1期)2015-02-27给纵波演示器的弹簧加保护装置物理教学探讨(2014年10期)2015-01-091988年澜沧—耿马地震前震源区应力状态分析地震研究(2014年3期)2014-02-27自制纵波演示仪的方法及其作用探讨中小学实验与装备(2009年5期)2009-12-02
