技术领域
[0001] 本
发明属于
地震数据分析领域，具体地说，涉及一种三维三分量矢量波场分离的方法和系统。
背景技术
[0002] 在三维三分量矢量（3 DIMENSION 3 COMPONENT，3D3C）
地震勘探来中，通常需要
检波器来接受来自地下介质的
地震波，包括纵波P和横波S。在其定义的O-Z-R-T三维
正交坐标系中，其Z分量作为垂直分量的纵波P，R分量和T分量作为
水平分量接收横波S，即沿X或Y水平方向的横波S中的一横波旋转变换到指向
震源-检波器方位的R分量，将另外一横波变换到与震源-检波器方位正交的T分量，该过程可称为RT旋转变换。
[0003] 当地下介质为各向同性且水平时，由于低速带的存在，可假设纵波P完全被Z分量接收，即纵波P具备真振幅，而RT旋转变换后的横波S
能量可以集中到R分量上，与此同时，剩余的噪音分量则集中到T分量上。这种假设使早期的三维三分量地震勘探在岩性油藏、缝洞油气藏、
煤炭以及
煤层气等复杂储层的发挥了重要作用。
[0004] 但是，在勘探条件以及勘探目标不断复杂的情况下，三维三分量（3DIMENSION 3 COMPONENT，3D3C）地震勘探的上述这种假设无法满足实际的要求，比如检波器接收的地震波发生“波型
泄漏”现象，水平分量R、T与垂直分量Z都接收到纵波P与横波S的部分偏振投影。产生“波型泄漏”现象通常有如下两个原因：
[0005] 1.纵波P不完全被Z分量接收
[0006] 当低速带较浅、偏移距较大或者地维存在三维倾斜面时，纵波P会斜交地面出射，从而导致Z、R、T分量上均能接收到纵波P的能量。
[0007] 2.横波S不完全被R分量接收
[0008] 当低速带较浅、偏移距较大或者地维存在三维倾斜面时，横波S会斜交地面出射，导致Z、R、T分量上均能接收到横波S的能量，其中，R分量和T分量上接收到的横波S分别称为快横波SV、慢横波SH。当地下介质为
各向异性介质如水平对称轴的横观各向同性介质（Horizontal transverse isotropic，HTI），或者倾斜对称轴的横观各向同性介质（Tilted Transverse Isotropic，TTI）时，会有横波分裂现象的产生，即不同含裂隙
地层的快、慢横波复合在一起形成复合波，此时快横波SV、慢横波SH分别为横波S复合波在R、T分量上的投影。
[0009] 可见，压制“波型泄漏”的效应是多分量地震
数据处理的一个重要方面，也是难点所在，即如何有效分离纵波P、快横波SV、慢横波SH。
[0010] 目前发展的
算法主要有两大类：一类是基于波场“域”转换的方法，典型的方法有在f-k域滤波、tao-p域滤波等；另一类是基于地震波振动特性（极化特性）不同的滤波方法，例如偏振方向滤波、协方差矩阵特征值分析或解析
信号分析极化滤波等。
[0011] 例如，发明
专利（公开号CN1404582A）是基于
波动方程上下行波的分解，在Dankbaar.J.W.M.发表于1985年的文章（Geophys.Prosp.33:970-986）方法的
基础上，通过波动方程的求散度和旋度的变换，分解涨缩P波和剪切S波；只是该专利技术克服了
Dankbarr方法对于浅层低、降速带不适应问题，通过浅层速度结构的求取，校正了浅层各向异性或不均一性对深层纵、横波场分离所产生的误差。
[0012] 专利（公开号CN101630017A）公开保护的是针对VSP数据的波场分离技术，通过仿射坐标变换实现。针对地面三分量地震勘探数据，该专利提出了一种仿射坐标变换和波场分离的实现方法与技术路线，但是该专利技术只能实现波场的分离，不能实现之后的再合成，无法真正有效地恢复波场的全振幅属性。
[0013] 专利（公开号CN101251603A）只能对二分量进行处理，忽略了实际采集数据的三分量特征；在不断速度分析和抽道集过程中，速度比模型会破坏该方法的效果；该方法损伤了三分量地震信号所具有的矢量特征，保幅性差，不利于后续的解释与反演。此外，在波型泄漏严重的地区，或地表复杂
信噪比较差的地区，尤其是浅层低、降速带厚度较薄的地区，该技术无法实现完全的波场分离。
[0014] 以上方法都能在一定程度上实现波场分离，获得压制噪声后的地震
波数据，但是由于其无法实现完全的波场分离，造成了地震数据分析时精确性不足。
发明内容
[0015] 本发明所要解决的技术问题是提供一种三维三分量矢量波场分离的方法和系统，基于两步仿射坐标系的旋转变换压制“波型泄漏”，对3D3C纵波P、快横波SV、慢横波SH的矢量波场进行分离，提高矢量波场的信噪比。
[0016] 为了解决上述技术问题，本发明提供了一种三维三分量矢量波场分离的方法，该方法包括：
[0017] 建立正交的O-Z-R-T三维坐标系，其中，R-O-Z平面与R-O-T平面正交，Z轴方向与地面垂直向上，所述R轴指向震源方向，O点为地震波地表接收点，T轴方向与所述R轴方向正切；检测纵波的波矢量方向、Z-O-T平面下慢横波的波矢量初始方向、纵波在R-O-Z平面的投影波矢量的波矢量方向和R-O-Z平面下快横波的波矢量方向；
[0018] 根据所述纵波在R-O-Z平面的投影矢量的波矢量方向、所述快横波的波矢量方向建立第一仿射坐标系，将纵波在R-O-Z平面的投影的偏振投影矢量、R-O-Z平面内快横波的偏振投影矢量旋转变换到所述第一仿射坐标系中，从而获得纵波在R-O-Z平面的投影矢量、快横波；
[0019] 根据所述R-O-Z平面向所述第一仿射坐标系的旋转变换，获取所述Z轴旋转到所述纵波在R-O-Z平面的投影矢量的波矢量方向的
角度，依据所述角度将所述慢横波的波矢量初始方向旋转得到所述慢横波的波矢量参考方向；
[0020] 根据所述慢横波的波矢量参考方向、所述纵波的波矢量方向建立第二仿射坐标系，将所述纵波在R-O-Z平面的投影矢量和所述T轴上的偏振投影矢量旋转变换到所述第二仿射坐标系中，从而获得纵波和慢横波。
[0021] 进一步地，所述方法包括：根据纵波单位波矢量和横波单位波矢量，对纵波在R-O-Z平面的投影矢量和快横波进行分解，获取所述纵波在R-O-Z平面的投影矢量与纵波单位波矢量平行的部分、快横波与横波单位波矢量平行的部分，分别输出为纵波在R-O-Z平面的投影有效信号分量和快横波的有效信号分量。
[0022] 进一步地，所述方法包括：根据所述纵波、快横波、慢横波进行振幅分离。
[0023] 进一步地，输出为纵波在R-O-Z平面的投影有效信号分量和快横波的有效信号分量之前包括：
[0024] 对所述纵波在R-O-Z平面的投影矢量和所述快横波所在象限的所有波矢量进行
叠加，获得纵波单位矢量和快横波单位矢量。
[0025] 进一步地，输出为纵波在R-O-Z平面的投影有效信号分量和快横波的有效信号分量包括：
[0026] 将纵波在R-O-Z平面的投影矢量按照平行纵波单位矢量进行分解、快横波按照平行快横波单位矢量进行分解，分别获得纵波在R-O-Z平面的投影有效信号分量和快横波的有效信号分量；
[0027] 将纵波在R-O-Z平面的投影矢量按照垂直于纵波单位矢量进行分解、快横波按照垂直于快横波单位矢量进行分解，分别获得纵波在R-O-Z平面的投影的噪音分量和快横波的噪音分量。
[0028] 进一步地，获得纵波单位矢量和快横波单位矢量的步骤包括：
[0029] 确定纵波在R-O-Z平面的投影矢量和快横波所在象限的角平分线矢量；
[0030] 对纵波在R-O-Z平面的投影矢量和快横波所在象限的所有波矢量进行叠加分别得到叠加后的纵波在R-O-Z平面的投影矢量和叠加后的快横波；
[0031] 根据叠加后的纵波在R-O-Z平面的投影矢量、叠加后的快横波以及两者各自所在象限内的角平分线矢量，分别获得纵波单位矢量和快横波单位矢量。
[0032] 进一步地，根据叠加后的纵波在R-O-Z平面的投影矢量、快横波以及两者各自所在象限内的角平分线矢量分别获得纵波单位矢量和快横波单位矢量的步骤中，包括：
[0033] 将纵波在R-O-Z平面的投影矢量和快横波所在象限的角平分线矢量分别作为初始纵波单位矢量和初始快横波单位矢量；
[0034] 分别在不同象限内对叠加后的纵波在R-O-Z平面的投影矢量和叠加后的快横波进行矢量差值和单位化处理，以对初始纵波单位矢量和初始快横波单位矢量进行
迭代处理，分别获得纵波单位矢量和快横波单位矢量。
[0035] 进一步地，根据所述纵波在R-O-Z平面的投影矢量的波矢量方向、所述快横波的波矢量方向建立第一仿射坐标系，将纵波在R-O-Z平面的投影的偏振投影矢量、R-O-Z平面内快横波的偏振投影矢量旋转变换到所述第一仿射坐标系中，从而获得纵波在R-O-Z平面的投影矢量、快横波的步骤包括：
[0036] 根据所述纵波在R-O-Z平面的投影矢量的波矢量方向、所述快横波的波矢量方向确定非正交的第一仿射坐标系；
[0037] 根据纵波在R-O-Z平面的投影矢量、快横波叠加形成的复合矢量与纵波在R-O-Z平面的投影矢量、快横波在R轴、Z轴上形成的偏振投影矢量叠加产生同一复合矢量，得到R-O-Z平面到所述第一仿射坐标系的第一仿射坐标张量矩阵；
[0038] 根据第一仿射张量矩阵，将纵波在R-O-Z平面的投影的偏振投影矢量、R-O-Z平面内快横波的偏振投影矢量从R-O-Z平面旋转变换到到第一仿射坐标系中，从而获得纵波在R-O-Z平面的投影矢量和快横波。
[0039] 进一步地，根据所述慢横波的波矢量参考方向、所述纵波的波矢量建立第二仿射坐标系，将所述纵波在R-O-Z平面的投影矢量和T轴上的偏振投影矢量旋转变换到所述第二仿射坐标系中从而获得纵波和慢横波的步骤包括：
[0040] 根据所述慢横波的波矢量参考方向、所述纵波的波矢量方向建立第二仿射坐标系；
[0041] 根据在第一仿射坐标系下、第二仿射坐标系下慢横波的振幅值相等，得到纵波在R-O-Z平面的投影矢量与T轴上的偏振投影矢量到第二仿射坐标系下的第二仿射坐标张量矩阵；
[0042] 根据第二仿射坐标张量矩阵，将所述纵波在R-O-Z平面的投影矢量和T轴上的偏振投影矢量旋转变换到所述第二仿射坐标系中，从而获得纵波和慢横波。
[0043] 为了解决上述技术问题，本发明提供了一种三维三分量矢量波场分离的系统，包括：处理设备和存储设备，其中，所述处理设备包括：
[0044] 配置
采样单元，用于建立正交的O-Z-R-T三维坐标系，其中，R-O-Z平面与R-O-T平面正交，Z轴方向与地面垂直向上，所述R轴指向震源方向，O点为地震波地表接收点，T轴方向与所述R轴方向正切；检测纵波的波矢量方向、Z-O-T平面下慢横波的波矢量初始方向、R-O-Z平面下纵波在R-O-Z平面的投影矢量的波矢量方向和R-O-Z平面下快横波的波矢量方向，发送至所述存储设备保存；
[0045] 第一波场分离单元，用于根据所述纵波在R-O-Z平面的投影矢量的波矢量方向、所述快横波的波矢量方向建立第一仿射坐标系，将纵波在R-O-Z平面的投影的偏振投影矢量、R-O-Z平面内快横波的偏振投影矢量旋转变换到所述第一仿射坐标系中，从而获得纵波在R-O-Z平面的投影矢量、快横波，并发送至所述存储设备保存；
[0046] 参考辅助单元，用于根据所述R-O-Z平面向所述第一仿射坐标系的旋转变换，获取所述Z轴旋转到所述纵波在R-O-Z平面的投影矢量的波矢量方向的角度，依据所述角度将所述慢横波的波矢量初始方向旋转得到所述慢横波的波矢量参考方向，并发送至所述存储设备保存；
[0047] 第二波场分离单元，用于根据所述慢横波的波矢量参考方向、所述纵波的波矢量方向建立第二仿射坐标系，将所述纵波在R-O-Z平面的投影矢量和所述T轴上的偏振投影矢量旋转变换到所述第二仿射坐标系中，从而获得纵波和慢横波。
[0048] 进一步地，所述系统包括：第三波场分离单元，用于根据纵波单位波矢量和横波单位波矢量，对纵波在R-O-Z平面的投影矢量和快横波进行分解，获取所述纵波在R-O-Z平面的投影矢量与纵波单位波矢量平行的部分、快横波与横波单位波矢量平行的部分，分别输出为纵波在R-O-Z平面的投影矢量的有效信号分量和快横波的有效信号分量。
[0049] 进一步地，所述系统包括：单位矢量获取单元，对纵波在R-O-Z平面的投影矢量和快横波所在象限的所有波矢量进行叠加，获得纵波单位矢量和快横波单位矢量并发送至所述存储设备保存。
[0050] 进一步地，所述系统包括：振幅分离单元，用于根据所述纵波、快横波、慢横波进行振幅分离。
[0051] 与现有的方案相比，本发明所获得的技术效果：
[0052] (1)基于两步仿射坐标系的旋转变换能够较好实现3D3C地震数据的纵波P、快横波SV、慢横波SH波的矢量波场分离。
[0053] (2)利用角平分线上的基矢量与逐点滑动时窗方法计算每个样点的单位波矢量，减少了计算时间，也提高了波矢方向求解的
精度。
[0054] (3)该方法在仿射坐标系的旋转变压制了噪音分量，提高了矢量波场的信噪比，所以较为适用于较低信噪比的数据处理。
附图说明
[0055] 图1为本发明
实施例的三维三分量矢量波场分离的方法
流程图；
[0056] 图2为本发明实施例的二维坐标系内矢量波场的示意图；
[0057] 图3为本发明实施例获得纵波单位矢量和快横波单位矢量的流程示意图；
[0058] 图4为本发明实施例第一次仿射坐标旋转过程示意图；
[0059] 图5为本发明实施例的三维坐标系内矢量波场的示意图；
[0060] 图6为本发明实施例第二次仿射坐标旋转处理流程图；
[0061] 图7(a)为某气田3D3C地震数据中的Z分量原始数据；
[0062] 图7(b)为某气田3D3C地震数据中的R分量原始数据；
[0063] 图7(c)为某气田3D3C地震数据中的T分量原始数据；
[0064] 图7(d)为波场分离后得到的纵波在R-O-Z平面的投影矢量PV；
[0065] 图7(e)为波场分离后分别得到的纵波P；
[0066] 图7(f)为
波长分离后的快横波SV；
[0067] 图7(g)为波长分离后的慢横波SH；
[0068] 图7(h)为O-Z-R平面内从Z分量上分离出的横波波型泄漏；
[0069] 图7(i)为O-PV-T平面内从PV分量上分离出的横波波型泄漏；
[0070] 图7(j)为O-Z-R平面内从R分量上分离出的纵波波型泄漏；
[0071] 图7(k)为O-PV-T平面内从T分量上分离出的纵波波型泄漏
[0072] 图8为本发明实施例三维三分量矢量波场分离系统结构示意图。
具体实施方式
[0073] 以下将配合图式及实施例来详细说明本发明的实施方式，藉此对本发明如何应用技术手段来解决技术问题并达成技术功效的实现过程能充分理解并据以实施。
[0074] 本发明的核心构思在于：基于两步仿射坐标系的旋转变换压制“波型泄漏”，对3D3C纵波P、快横波SV、慢横波SH的矢量波场进行分离，提高了矢量波场的信噪比。本发明中，首先建立了正交的O-Z-R-T三维坐标系，其中，R-O-Z（ROZ）平面与R-O-T（ROT）平面正交，R-O-Z平面为纵波P的入射面，Z轴方向与地面垂直向上，R轴指向震源方向，O点为地震波地表接收点，T轴方向与R轴方向正切，快横波SV、慢横波SH分别为横波S在ROZ平面与ROT平面的投影矢量。
[0075] 如图1所示，为本发明实施例的三维三分量矢量波场分离的方法流程图，包括：
[0076] 步骤101、确定ROZ平面中纵波单位矢量和快横波单位矢量。
[0077] 如图2所示，为本发明实施例的二维坐标系内矢量波场的示意图。在ROZ平面内，Z轴方向与地面垂直向上，R轴指向震源方向，因此R轴方向也可称为径向分量方向，O点为地震波地表接收点。因此，在确定纵波单位矢量Dp和快横波单位矢量Ds时，可以对纵波在ROZ平面的投影矢量PV和快横波SV所在象限的所有波矢量进行叠加。本领域技术人员应当可以理解，建立ROZ平面并不局限于这种特定的方式。
[0078] 进一步地，如图3所示，为本发明实施例获得纵波单位矢量和快横波单位矢量的流程示意图，通过如下具体步骤来实现：
[0079] 步骤111、确定纵波在ROZ平面的投影矢量和快横波所在象限的角平分线矢量；
[0080] 如果所述纵波在ROZ平面的投影矢量的所有波矢量位于一、三象限，取第一象限的角平分线矢量 eR、eZ分别为R与Z的基矢量。
[0081] 如果所述快横波的所有波矢量位于一、三象限，取第一象限的角平分线矢量
[0082] 如果所述快横波的所有波矢量位于二、四象限，取第四象限的角平分线矢量
[0083] 如果所述纵波在ROZ平面的投影矢量的所有波矢量位于二、四象限，取第四象限的角平分线矢量
[0084] 步骤121、对纵波在ROZ平面的投影矢量和快横波所在象限的所有波矢量进行叠加分别得到叠加后的纵波在ROZ平面的投影矢量和叠加后的快横波；
[0085] 如果所述纵波在ROZ平面的投影矢量的所有波矢量位于一、三象限，取第一象限的角平分线矢量 把第一象限内与D夹角小于90°的全部波矢量叠加得到第一象限内叠加后的纵波在ROZ平面的投影矢量d1；把第三象限内与D夹角大于90°的全部波矢量叠加得到第三象限内叠加后的纵波在ROZ平面的投影矢量d2；
[0086] 所述快横波的所有波矢量位于二、四象限，取第四象限的角平分线矢量把第二象限内与D夹角小于90°的全部波矢量叠加得到第二象限内
叠加后的快横波d1；把第四象限内与D夹角大于90°的全部波矢量叠加得到第四象限内叠加后的快横波d2；
[0087] 如果所述快横波的所有波矢量位于一、三象限，取第一象限的角平分线矢量把第一象限内与D夹角小于90°的全部波矢量叠加得到第一象限内
叠加后的快横波d1；把第三象限内与D夹角大于90°的全部波矢量叠加得到第三象限内叠加后的快横波d2；
[0088] 如果所述纵波在ROZ平面的投影矢量的所有波矢量位于二、四象限，取第四象限的角平分线矢量 把第二象限内与D夹角小于90°的全部波矢量叠加得到第二象限内叠加后的纵波在ROZ平面的投影矢量d1；把第四象限内与D夹角大于90°的全部波矢量叠加得到第四象限内叠加后的纵波在ROZ平面的投影矢量d2。
[0089] 步骤131、根据叠加后的纵波在ROZ平面的投影矢量、叠加后的快横波，以及其两者各自所在象限内的角平分线矢量，分别获得纵波单位矢量和快横波单位矢量。
[0090] 进一步地，步骤131可以具体为：
[0091] 将纵波在ROZ平面的投影矢量和快横波所在象限的角平分线矢量分别作为初始纵波单位矢量和快横波单位矢量；
[0092] 即：根据纵波在ROZ平面的投影矢量、快横波所在的象限，选择第一象限的角平分线矢量 或第四象限的角平分线矢量 作为初始纵波单位矢量、初始快横波单位矢量。
[0093] 分别在不同象限内对叠加后的纵波在ROZ平面的投影矢量和快横波进行矢量差值和单位化处理，并对初始纵波单位矢量和快横波单位矢量进行迭代处理，分别获得纵波单位矢量和快横波单位矢量。
[0094] 即可通过差值和单位化公式：
[0095] DSV or PV=(d1-d2)/|d1-d2| （1）
[0096] 确定迭代纵波单位矢量或快横波单位矢量DPV，替换初始纵波单位矢量、初始快横波单位矢量，以在上述对应的象限内，与迭代纵波单位矢量或快横波单位矢量DPV夹角小于90°全部波矢量叠加得到叠加后的纵波在ROZ平面的投影矢量或快横波d1；以在上述对应的象限内，与迭代纵波单位矢量或快横波单位矢量DPV夹角大于90°的全部波矢量叠加得到叠加后的纵波在ROZ平面的投影矢量或快横波d2，再次重新根据公式（1）确定一个新的迭代纵波单位矢量或快横波单位矢量DPV，重复上述迭代过程2-3次即可确定最终的纵波单位矢量Dp和快横波单位矢量Ds。
[0097] 下面以纵波在ROZ平面的投影矢量的所有波矢量位于一、三象限，快横波的所有波矢量位于二、四象限为例详细如何确定精确的纵波单位矢量Dp和快横波单位矢量Ds。
[0098] 如果所述纵波在ROZ平面的投影矢量的所有波矢量位于一、三象限，取第一象限的角平分线矢量 把第一象限内与D夹角小于90°的全部波矢量叠加得到第一象限内叠加后的纵波在ROZ平面的投影矢量d1；把第三象限内与D夹角大于90°的全部波矢量叠加得到第三象限内叠加后的纵波在ROZ平面的投影矢量d2，然后对第一象限内叠加后的纵波在ROZ平面的投影矢量d1和第三象限内叠加后的纵波在ROZ平面的投
影矢量d2进行矢量差值和单位化处理公式（1），并对初始纵波单位矢量进行迭代处理，获得精确的纵波单位矢量Dp；
[0099] 如果所述快横波的所有波矢量位于二、四象限，取第四象限的角平分线矢量把第二象限内与D夹角小于90°的全部波矢量叠加得到第二象限内
叠加后的快横波d1；把第四象限内与D夹角大于90°的全部波矢量叠加得到第四象限内叠加后的快横波d2，然后对第二象限内叠加后的快横波d1和第四象限内叠加后的快横波d2进行矢量差值和单位化处理公式（1），并对初始快横波单位矢量进行迭代处理，获得精确的快横波单位矢量Ds。
[0100] 步骤102、在ROZ平面内进行第一次仿射坐标旋转。
[0101] 进一步地，在实现第一次仿射坐标变换时，根据纵波在ROZ平面的投影矢量的波矢量方向、快横波的波矢量方向建立第一仿射坐标系，将纵波在ROZ平面的投影的偏振投影矢量、ROZ平面内快横波的偏振投影矢量旋转变换到第一仿射坐标系中，从而获得纵波在ROZ平面的投影矢量、快横波。
[0102] 进一步地，如图4所示，为本发明实施例第一次仿射坐标旋转过程示意图，其可以具体包括：
[0103] 步骤112、根据纵波在ROZ平面的投影矢量的波矢量方向、快横波的波矢量方向确定非正交的第一仿射坐标系；
[0104] 纵波在ROZ平面的投影矢量PV的波矢量方向、快横波SV的波矢量方向如图2中PV和SV箭头方向所示，根据这两个方向组成非正交的第一仿射坐标系PV-O-SV。
[0105] 步骤122、根据纵波在ROZ平面的投影矢量、快横波叠加形成的复合矢量与纵波在ROZ平面的投影矢量、快横波在R轴、Z轴上形成的偏振投影矢量叠加产生同一复合矢量，得到ROZ平面到所述第一仿射坐标系的第一仿射坐标张量矩阵；
[0106] 具体地，如图2所示，可以根据R轴、Z轴上形成的偏振投影矢量AR、AZ叠加形成复合矢量A：
[0107] A=AReR+AzeZ （2）
[0108] R轴、Z轴上形成的偏振投影矢量AR、AZ由纵波在ROZ平面的投影矢量PV、快横波SV各自在R轴、Z轴上的偏振投影矢量PVR、PVZ、SVR、SVZ叠加而成：
[0109] AR+AZ＝(PVR+SVR)eR+(PVZ+SVZ)eZ （3）；
[0110] 同时，由于该复合矢量A也可以纵波在ROZ平面的投影矢量和快横波叠加形成：
[0111] A=PV+SV （4）
[0112] 因此，根据公式（2）（3）（4）可得到第一仿射坐标张量矩阵：
[0113]
[0114] 其中，eR、eZ、ePV、eSV分别为R轴的基矢量、Z轴的基矢量、纵波在ROZ平面的投影矢量的波矢量的基矢量、快横波量的波矢量的基矢量，（ePV，eZ)、（eSV，eZ)、（ePV，eR）、（eSV，eR）分别表示纵波在ROZ平面的投影矢量的波矢量方向与Z轴的夹角、快横波的波矢量方向与Z轴的夹角、纵波在ROZ平面的投影矢量的波矢量方向与R轴的夹角、纵波在ROZ平面的投影矢量的波矢量方向与R轴的夹角。
[0115] 步骤132、根据第一仿射坐标张量矩阵，将纵波在ROZ平面的投影矢量、ROZ平面内快横波的偏振投影矢量从ROZ平面旋转变换到第一仿射坐标系中，获得纵波在ROZ平面的投影矢量和快横波。
[0116] 进一步地，本发明实施例的方法中，可以在第一仿射坐标旋转变换得到纵波在ROZ平面的投影矢量PV和快横波SV后，或者在完成第二仿射坐标变换得到纵波和慢横波后，将纵波在ROZ平面的投影矢量PV按照平行纵波单位矢量Dp、快横波按照平行于横波单位矢量Ds进行分解，获得纵波在ROZ平面的投影矢量的有效信号分量PVt和有效快横波的有效信号分量SVt。如下所示：
[0117]
[0118] 进一步地，本发明实施例的方法中，可以在第一仿射坐标旋转变换得到纵波在ROZ平面的投影矢量PV和快横波SV后，或者在完成第二仿射坐标变换得到纵波和慢横波后，将纵波在ROZ平面的投影矢量PV按照垂直于纵波单位矢量Dp、快横波SV分别按照垂直于快横波单位矢量Ds进行分解，获得纵波在ROZ平面的投影矢量的噪音分量PVn和快横波的噪音分量SVn。
[0119]
[0120] 步骤103、O-Z-R-T三维坐标系下进行第二次仿射坐标旋转。
[0121] 进一步地，步骤103具体可以根据慢横波的波矢量参考方向、纵波的波矢量方向建立第二仿射坐标系，将纵波在ROZ平面的投影矢量和T轴上的偏振投影矢量旋转变换到第二仿射坐标系中，从而获得纵波和慢横波。
[0122] 本实施例中，如图5所示，为本发明实施例的三维坐标系内矢量波场的示意图，三维坐标系下O-Z-R-T三维坐标系中T轴方向与R轴方向正切，因此，T轴方向也可称为切向分量方向。O-Z-R-T三维坐标系中，ROT平面与ROZ平面正交。
[0123] 进一步地，如图6所示，为本发明实施例第二次仿射坐标旋转处理流程图。其可以具体包括：
[0124] 步骤113、根据慢横波的波矢量参考方向、纵波的波矢量方向建立第二仿射坐标系；
[0125] 对于三维空间内的横波来说，地下介质的TTI或者HTI特性会导致横波发生分裂形成横波S分别在O-Z-R面内的快横波SV与O-Z-T面内的慢横波SH，复合波S=SV+SH。
[0126] SV=SVZeZ+SVReR,SH=SHZeZ+SHTeT （8）
[0127] SVZeZ与SHZeZ是快横波SV和慢横波SH分别在Z分量上的偏振投影即为其波型泄漏。快横波SV在O-Z-R平面内第一步仿射坐标系旋转变换过程中已经得到，所以在第二步仿射坐标系旋转变化过程中需要分离慢横波SH。
[0128] 如图5所示，因为慢横波SH的振幅SHZ是其在O-Z-R平面内的唯一投影振幅，所以在O-Z-R平面内第一步仿射坐标系旋转变换过程中，SHZ已经转移到纵波在ROZ平面的投影矢量PV的波矢量方向上，而慢横波SH也发生了旋转变化，其旋转变换后的矢量为慢横波投影参考矢量SH’，该慢横波投影参考矢量SH’的方向即为慢横波SH的波矢量参考方向。
[0129] 在三维坐标系O-Z-R-T内，纵波P矢量为
[0130] P=PZeZ+PReR+PTeT (9)
[0131] 参见图2所示，纵波在ROZ平面的投影矢量PV在ROZ平面下可表示为:
[0132] PV=PZeZ+PReR (10)
[0133] 则三维空间内纵波P矢量又可写为：
[0134] P=PV+PTeT (11)
[0135] 其中，eT为T方向上的基矢量。
[0136] 因此，本发明实施例中，以三维空间内纵波P的波矢量方向与慢横波投影参考矢量SH’的波矢量方向即慢横波的波矢量参考方向，建立第二仿射坐标系P-O-SH’。
[0137] 步骤123、根据在第一仿射坐标系下、第二仿射坐标系下慢横波的振幅值相等，得到纵波在ROZ平面的投影矢量与T轴上的偏振投影矢量到所述第二仿射坐标系下的第二仿射坐标张量矩阵；
[0138] 慢横波SH在T轴上的偏振投影矢量SHT为：
[0139] SHT=SHTeT (12)
[0140] SHT始终正交平面ROZ，慢横波投影参考矢量SH’和慢横波SH两者振幅值相同，只是矢量的波矢量方向变换到了纵波在ROZ平面的投影矢量P的波矢量方向上，因此，可以根据纵波在ROZ平面的投影矢量PV的波矢量方向与T轴构成PV-O-T坐标系，得到PV-O-T坐标系到第二仿射坐标系P-O-SH’的第二仿射坐标张量矩阵。
[0141] 进一步地，本实施例中第二仿射坐标张量矩阵为：
[0142]
[0143] 其中，ST表示T轴上的偏振投影矢量，（eP，ePV）、（eP，eT）分别表示纵波与纵波在ROZ平面的投影矢量的波矢量方向、T轴的夹角，（e，SH，ePV）、（e，SH，eT）分别表示慢横波投影的波矢量参考方向与纵波在ROZ平面的投影矢量的波矢量方向、T轴的夹角，e’SH表示慢横波投影参考矢量的基矢量。
[0144] 步骤133、根据第二仿射坐标张量矩阵，将纵波在ROZ平面的投影矢量和T轴上的偏振投影矢量旋转变换到第二仿射坐标系中，从而获得纵波和慢横波。
[0145] 由上述公式（10）可得出纵波的波矢量方向与纵波在ROZ平面的投影矢量的波矢量方向、T轴的夹角分别为（eP，ePV）、（eP，eT）；慢横波的波矢量参考方向与纵波在ROZ平面的投影矢量的波矢量方向、T轴的夹角分别为（e，SH，ePV）、（e，SH，eT）。
[0146] 进一步地，在本发明的上述实施例中，采用逐点滑动时窗法，即在所述采样点上下各开半个时窗，对整个时窗长度内计算的波矢量应用到时窗中心点上。
[0147] 进一步地，在本发明的上述实施例中，在获得纵波P、快横波SV、慢横波SH后，可以对其振幅进行分离。
[0148] 图7(a)、(b)、(c)为某气田3D3C地震数据的Z、R、T分量原始数据，虽然原始数据信噪比很低导致波型泄漏现象难以用肉眼识别，但是它是客观存在的。通过本方法的应用可以帮助我们压制波型泄漏现象。图7(d)、(e)是波场分离后分别得到的纵波在ROZ平面的投影矢量PV与纵波P，可以看出快横波SV与慢横波SH的波型泄漏得到很大程度的压制，特别是在椭圆框标识的
位置波形泄漏压制作用更明显。从图7(f)、(g)可以看出在箭头标识的位置，纵波P的初至得到很大程度的压制。图7(h)、(i)分别为Z分量上压制的快横波SV与慢横波SH的波型泄漏。从图7(f)与(h)以及(g)与(i)之间可以看出泄漏的横波波场与纯横波波场之间在同相轴的时距特征上有很强的相似性，同样在7(j)与(k)显示的R、T分量上压制的纵波P波形泄漏与图7(e)的纯纵波波场之间也能够发现同相轴时距特征的相似性。这也说明，本发明的技术方案能够有效分离纵波P、快横波SV、慢横波SH。
[0149] 如图8所示，为本发明实施例三维三分量矢量波场分离的系统结构示意图。该系统包括：处理设备80和存储设备81，其中，所述处理设备80包括：配置采样单元801、第一波场分离单元802、参考辅助单元803、第二波场分离单元804。
[0150] 配置采样单元801用于建立正交的O-Z-R-T三维坐标系，其中，ROZ平面与ROT平面正交，Z轴方向与地面垂直向上，所述R轴指向震源方向，O点为地震波地表接收点，T轴方向与所述R轴方向正切；检测纵波的波矢量方向、Z-O-T平面下慢横波的波矢量初始方向、ROZ平面下纵波在ROZ平面的投影矢量的波矢量方向和ROZ平面下快横波的波矢量方向，发送至所述存储设备86保存；
[0151] 第一波场分离单元802用于从所述存储设备86获取纵波在ROZ平面的投影矢量的波矢量、所述快横波的波矢量，并根据纵波在ROZ平面的投影矢量的波矢量方向、所述快横波的波矢量方向建立第一仿射坐标系，将ROZ平面内纵波在ROZ平面的投影的偏振投影矢量、ROZ平面内快横波的偏振投影矢量旋转变换到所述第一仿射坐标系中，从而获得纵波在ROZ平面的投影矢量、快横波，并发送至所述存储设备86保存；
[0152] 参考辅助单元803用于根据所述ROZ平面向第一仿射坐标系的旋转变换，获取Z轴旋转到纵波在ROZ平面的投影矢量的波矢量方向的角度，依据所述角度将由所述存储设备86获取的慢横波的波矢量初始方向旋转得到慢横波的波矢量参考方向，并发送至所述存储设备86保存；
[0153] 第二波场分离单元804用于由所述存储设备86获取所述慢横波的波矢量参考方向、所述纵波的波矢量方向并据此建立第二仿射坐标系，将纵波在ROZ平面的投影矢量和所述T轴上的偏振投影矢量旋转变换到第二仿射坐标系中，从而获得纵波和慢横波并发送至存储设备86保存。
[0154] 进一步地，本发明实施例系统还可以包括：第三波场分离单元（图中未示出），该第三波场分离单元用于根据纵波单位矢量、快横波单位矢量对从所述存储设备86获取的纵波在ROZ平面的投影矢量、快横波进行分解，获取纵波在ROZ平面的投影矢量、快横波与纵波单位波矢量和横波单位波矢量平行的部分，输出为有效的纵波在ROZ平面的投影矢量的有效信号分量和快横波的有效信号分量。
[0155] 进一步地，本发明实施例系统还可以包括：单位矢量获取单元（图中未示出），该单位矢量获取单元用于对由所述存储设备86获取的纵波在ROZ平面的投影矢量和快横波所在象限的所有波矢量进行叠加，获得纵波单位矢量和快横波单位矢量并发送至所述存储设备86保存。
[0156] 进一步地，本发明实施例系统还可以包括：振幅分离单元（图中未示出），用于根据由所述存储设备86获取的纵波、快横波、慢横波进行振幅分离。
[0157] 上述说明示出并描述了本发明的若干优选实施例，但如前所述，应当理解本发明并非局限于本文所披露的形式，不应看作是对其他实施例的排除，而可用于各种其他组合、
修改和环境，并能够在本文所述发明构想范围内，通过上述教导或相关领域的技术或知识进行改动。而本领域人员所进行的改动和变化不脱离本发明的精神和范围，则都应在本发明所附
权利要求的保护范围内。
