一种大倾角地区的垂直地震数据桥式标定方法
【专利摘要】本发明是油气地震勘探的垂直地震桥式标定方法，将NMO数据分离成两部分，其一为数据格式SEGY的卷头、道头等说明部分，其二为纯数值部分，对NMO数据纯数值部分进行NMO数据转换，卷头、道头等说明信息，与转换后矩阵数据拼合，得到SEGY格式的新NMO地震数据，对VSP桥式标定。本发明与目前的VSP桥式标定方法不同的是，标定结果中不再需要VSP走廊叠加剖面，这是因为通过上述系列转换，已将井口附近的地震反射轴转移至NMO剖面一侧，与地面地震剖面可直接进行对比，在回避大倾角地区走廊不能同相叠加问题的同时，得到更直观、更可靠的桥式标定结果。
【专利说明】一种大倾角地区的垂直地震数据桥式标定方法
【技术领域】
[0001] 本发明属于非常规油气地震勘探领域，是一种利用转换正常动校正（NMO)数据进 行大倾角地区的垂直地震桥式标定方法。
【背景技术】
[0002] 桥式标定是垂直地震（VSP)资料的重要用途之一，是利用VSP数据建立地面地震 反射与井中地层之间的关系，通常的做法是基于VSP时深关系将VSP的NMO剖面和走廊叠 加剖面、过井地震剖面、测井曲线、岩合柱状图组合起来，得到VSP桥式标定成果。
[0003] 随着近些年来VSP应用范围的不断扩大，在许多大倾角地区进行了 VSP采集，现有 的VSP走廊标定流程为：首先处理得到的NMO数据，然后在近井口位置沿某一时间（或深度） 进行窄走廊切除，进而利用数据叠加的方法得到走廊叠加剖面。由于地层倾角的存在，数据 NMO波场中上行波反射同相轴与时间轴(或深度轴）存在夹角，因此利用现有的方法不能同 相叠加得到可靠的VSP走廊剖面，并且现有的桥式标定方法主要是基于VSP走廊和过井地 震剖面的波阻特征进行对比标定，走廊不可靠就大大的影响了 VSP桥式标定的准确性。
【发明内容】
[0004] 本发明的目的是提供一种标定准确、可靠性高的利用转换NMO数据进行大倾角地 区的垂直地震数据桥式标定方法。
[0005] 本发明通过以下步骤实现：
[0006] 1)采集VSP数据，在VSP处理中，进行TAR补偿(真振幅恢复）和道均衡处理，得到 保真的转换前NMO数据；
[0007] 将NMO数据分离成两部分，其一为数据格式SEGY的卷头、道头等说明部分，其二为 纯数值部分；
[0008] 所述的TAR补偿和道均衡处理是保证剖面中不同位置的能量相对稳定，不因转换 处理出现能量不均衡现象。
[0009] 所述的处理在时间域和深度域中方法完全相同。
[0010] 所述的数据格式为SEGY格式。
[0011] 2)按照以下方式对NMO数据纯数值部分进行NMO数据转换：
[0012] (1)将NMO纯数值部分以矩阵的方式排列，矩阵的行数为最大采样点数，列数为 NMO数据所包含的总道数；
[0013] (2)利用地震处理中初至拾取的方法，拾取NMO纯数值部分数据的原始初至时间 序列，每一道对应一个初至时间；
[0014] (3)对初至时间序列进行整数采样计算，使各点时间值均为采样率的整数倍，得到 新初至时间序列T' :
[0015] 所述的计算依下式进行：
[0016] r, = int(会)* ddt U )
[0017] 其中T'为计算得到的新初至时间序列，T为步骤（2)拾取的原始初至时间序列， A t为数据的采样率，int为取整函数；
[0018] 所述的采样率为数据记录时间除以最大采样点数。
[0019] (4)依据新初至时间序列，将步骤（1)矩阵的每一道进行分段，得到包含系列子矩 阵的新矩阵；
[0020] (5)将得到新矩阵进行数据变换运算；
[0021] 所述的数据变换运算依下述方法进行：
[0022] 首先取出矩阵的右下三角阵，然后对三角阵进行左右翻转，并将翻转后三角阵中 各子矩阵顺序不变的重新排布回原三角阵的子矩阵位置，进而将新三角阵替换步骤（5)中 新矩阵的右下三角部分，最后展开新矩阵中各个子矩阵，得到与步骤（1)大小完全相同的转 换后矩阵数据。
[0023] (6)将步骤1)中分离得到的卷头、道头等说明信息，与转换后矩阵数据拼合，得到 SEGY格式的新NMO地震数据；
[0024] 3)按照以下方式进行VSP桥式标定：
[0025] 将过井地震剖面、测井数据、录井剖面、以及步骤2)得到的转换后矩阵数据以图像 方式输出，然后拼合，再以实钻分层数据进行VSP标定。
[0026] 本发明与目前的VSP桥式标定方法不同的是，标定结果中不再需要VSP走廊叠加 剖面，这是因为通过上述系列转换，已将井口附近的地震反射轴转移至NMO剖面一侧，与地 面地震剖面可直接进行对比，在回避大倾角地区走廊不能同相叠加问题的同时，得到更直 观、更可靠的桥式标定结果。
【专利附图】
【附图说明】
[0027] 图1原始NMO数据(左）和转换后NMO数据(右）示意图；
[0028] 图2实际VSP资料的原始的NMO波场(左）和转换后的NMO波场(右）；
[0029] 图3原始的NMO波场的近井口切除剖面(左)、VSP走廊叠加剖面（中）、VSP走廊镶 嵌到过井地震剖面(右）；
[0030] 图4转换后NMO波场与过井地震剖面拼接对比；
[0031] 图5实验地区某井转换NMO波场的VSP标定实例。
【具体实施方式】
[0032] 本发明【具体实施方式】通过以下步骤实现：
[0033] 1)利用采集得到的VSP数据，进行常规地震资料处理，其中TAR补偿和道均衡处理 两步流程影响到剖面中的能量分布，要通过对比观察的方式确定最佳参数，进而处理得到 保真的转换前NMO数据；
[0034] 将NMO数据分离成两部分，其一为数据格式SEGY的卷头、道头等说明部分，其二为 纯数值部分。
[0035] 2)按照以下方式对NMO数据纯数值部分进行数据转换：
[0036] (1)将上步得到的NMO纯数值部分记作X，它由mXn个元素x构成，其中x为数据 样点值，双下标指示元素在矩阵中的位置，如X 11在矩阵中的位置是第1行第1列，其物理意 义是NMO数据中第1道的第1个样点值，以此类推，n为NMO数据所包含的总道数，m为数据 最大采样点数，如下式所示：
[0037]
【权利要求】
1. 一种大倾角地区的垂直地震数据桥式标定方法，特点是通过以下步骤实现： 1) 采集VSP数据，在VSP处理中，进行TAR补偿和道均衡处理，得到保真的转换前NMO 数据； 将NMO数据分离成两部分，其一为数据格式SEGY的卷头、道头等说明部分，其二为纯数 值部分； 2) 按照以下方式对NMO数据纯数值部分进行NMO数据转换： (1) 将NMO纯数值部分以矩阵的方式排列，矩阵的行数为最大采样点数，列数为NMO数 据所包含的总道数； (2) 利用地震处理中初至拾取的方法，拾取NMO纯数值部分数据的原始初至时间序列， 每一道对应一个初至时间； (3) 对初至时间序列进行整数采样计算，使各点时间值均为采样率的整数倍，得到新初 至时间序列T': (4) 依据新初至时间序列T'，将步骤（1)矩阵的每一道进行分段，得到包含系列子矩 阵的新矩阵； (5) 将得到新矩阵进行数据变换运算； (6) 将步骤1)中分离得到的卷头、道头等说明信息，与转换后矩阵数据拼合，得到SEGY 格式的新NMO地震数据； 3) 按照以下方式进行VSP桥式标定： 将过井地震剖面、测井数据、录井剖面、以及步骤2)得到的转换后矩阵数据以图像方式 输出，然后拼合，再以实钻分层数据进行VSP标定。
2. 根据权利要求1的方法，特点是步骤1)所述的TAR补偿和道均衡处理是保证剖面中 不同位置的能量相对稳定，不因转换处理出现能量不均衡现象。
3. 根据权利要求1的方法，特点是步骤1)所述的道均衡处理在时间域和深度域中方法 完全相同。
4. 根据权利要求1的方法，特点是步骤1)所述的NMO数据格式为SEGY格式。
5. 根据权利要求1的方法，特点是步骤2)所述的整数采样计算依下式进行：
其中：T'为计算得到的新初至时间序列，T为步骤（2)拾取的原始初至时间序列，At 为数据的采样率，int为取整函数； 所述的采样率为数据记录时间除以最大采样点数。
6. 根据权利要求1的方法，特点是步骤2)所述的数据变换运算依下述方法进行： 首先取出矩阵的右下三角阵，然后对三角阵进行左右翻转，并将翻转后三角阵中各子 矩阵顺序不变的重新排布回原三角阵的子矩阵位置，进而将新三角阵替换步骤（5)中新矩 阵的右下三角部分，最后展开新矩阵中各个子矩阵，得到与步骤（1)大小完全相同的转换后 矩阵数据。
【文档编号】G01V1/32GK104375174SQ201310356189
【公开日】2015年2月25日 申请日期:2013年8月15日 优先权日:2013年8月15日
【发明者】蔡志东, 张庆红, 彭继新, 王冲 申请人:中国石油天然气集团公司, 中国石油集团东方地球物理勘探有限责任公司