古河道河床侵蚀深度自动识别方法
【技术领域】
[0001] 本发明涉及地球物理地震储层预测技术领域，特别涉及一种古河道河床侵蚀深度 自动识别方法。
【背景技术】
[0002] 暗河岩溶系统是高效油藏发育的有利区。暗河岩溶系统的发育受控于地表古水 系。地表水系深切部位是暗河岩溶系统的泄水场所，为岩溶作用提供水动力条件，有利于形 成大规模的暗河岩溶储集体。
[0003] 虽然河床侵蚀作用对河床两侧及底部位置的碳酸盐岩岩溶储层发育有建设作用， 但是侵蚀深度过深，容易破坏河床底部位置的盖层的封盖能力。即使河床底部下方的岩溶 储层发育，也不利于油气成藏。因此，古河道及其侵蚀深度的识别能够指导识别有利于油气 成藏的暗河岩溶储层。
[0004]目前，三维地震数据的相干、振幅、分频等地震属性结合切片技术：水平切片、沿层 切片及地层切片（雷明，王建功，刘彩艳等：地震切片技术在安达工区的应用、石油地球物 理勘探，2010,45(3) :418~422 ;李清仁，陈守田，张财等：海拉尔盆地的古河道砂体识别方 法、石油地球物理勘探，2006,41 (4) :402~4046 ;韩喜，高兴友，车延信等：利用地震属性沿 层分析方法研宄河流相沉积环境、石油地球物理勘探，2007,42(1) :120~124)虽然能够刻 画河道平面或空间展布，但是不能估算古河道河床的侵蚀深度。碳酸盐岩古河道河床底部 与围岩地层岩性差异大，波阻抗呈连续的强波峰反射，河道边界处呈杂乱弱反射，整个河床 底部地震反射轴连续，可以沿河道连续追踪。但是河道分布广、曲折度高、分支小河道发育， 人工追踪工作量极大。
【发明内容】
[0005] 本发明实施例提供了一种古河道河床侵蚀深度自动识别方法，用以实现古河道河 床侵蚀深度的自动识别，该方法包括：
[0006] 在河床地震剖面反射波的最大波峰处选取一个地震响应点作为种子点；
[0007] 计算种子点所在地震反射层面的相干和倾角；
[0008] 根据种子点所在地震反射层面的相干和倾角，驱动种子点向位于同一主测线上、 联络测线增大和减少的两个方向，自动追踪当前主测线上、与种子点相位相同的外推地震 响应点；
[0009] 计算每个外推地震响应点所在地震反射层面的相干和倾角，计算每个外推地震响 应点和与该外推地震响应点相邻的8个地震响应点，共9个地震响应点的相干平均值和方 差；
[0010] 根据相干平均值和方差，判断每个外推地震响应点是否到了河道边界；当判断外 推地震响应点未到河道边界时，根据每个外推地震响应点所在地震反射层面的相干和倾 角，继续以每个外推地震响应点为新的种子点，向联络测线增大和减少的两个方向，自动追 踪与新的种子点相位相同的外推地震响应点；当判断外推地震响应点到了河道边界处时， 停止追踪，得到当前主测线上所有的种子点；
[0011] 根据当前主测线上所有种子点的相干和倾角，驱动当前主测线上的所有种子点， 向主测线增大和减少的方向，自动追踪与当前主测线相邻主测线上、与种子点相位相同的 外推地震响应点，判断每个外推地震响应点是否到了河道边界，当判断外推地震响应点到 了河道边界处时，停止追踪，得到相邻主测线上的所有种子点，直至得到河道工区范围内所 有主测线上的种子点；
[0012] 根据河道工区范围内所有主测线上的种子点，计算古河道河床的侵蚀深度。
[0013] 在一个实施例中，根据种子点所在地震反射层面的相干和倾角，驱动种子点向位 于同一主测线上、联络测线增大和减少的两个方向，自动追踪当前主测线上、与种子点相位 相同的外推地震响应点，包括：
[0014] 取联络测线增大方向的6个地震响应点中，最大的相干值及其对应的倾角，向位 于同一主测线上、联络测线增大的方向，自动追踪当前主测线上、与种子点相位相同的外推 地震响应点，取联络测线减少方向的6个地震响应点中，最大的相干值及其对应的倾角，向 位于同一主测线上、联络测线减少的方向，自动追踪当前主测线上、与种子点相位相同的外 推地震响应点。
[0015] 在一个实施例中，在根据当前主测线上所有种子点的相干和倾角，驱动当前主测 线上的所有种子点，向主测线增大和减少的方向，自动追踪与当前主测线相邻主测线上、与 种子点相位相同的外推地震响应点中，驱动一个种子点向主测线增大和减少的方向，自动 追踪与当前主测线相邻主测线上、与种子点相位相同的外推地震响应点包括：
[0016] 取主测线增大方向的6个地震响应点中，最大的相干值及其对应的倾角，向位于 同一联络测线上、主测线增大的方向，自动追踪与当前主测线相邻主测线上、与种子点相位 相同的外推地震响应点，取主测线减少方向的6个地震响应点中，最大的相干值及其对应 的倾角，向位于同一联络测线上、主测线减少的方向，自动追踪与当前主测线相邻主测线 上、与种子点相位相同的外推地震响应点。
[0017] 在一个实施例中，根据当前主测线上所有种子点的相干和倾角，驱动当前主测线 上的所有种子点，向主测线增大和减少的方向，自动追踪与当前主测线相邻主测线上、与种 子点相位相同的外推地震响应点，判断每个外推地震响应点是否到了河道边界，当判断外 推地震响应点到了河道边界处时，停止追踪，得到相邻主测线上的所有种子点，直至得到河 道工区范围内所有主测线上的种子点，包括：
[0018] 根据当前主测线上边界种子点的相干和倾角，驱动边界种子点向相邻主测线自动 追踪，得到相邻主测线上的、与边界种子点相位相同的种子点；
[0019] 根据当前主测线上非边界种子点的相干和倾角，驱动非边界种子点向相邻主测线 自动追踪，得到相邻主测线上的、与非边界种子点相位相同的种子点。
[0020] 在一个实施例中，根据当前主测线上边界种子点的相干和倾角，驱动边界种子点 向相邻主测线自动追踪，得到相邻主测线上的、与边界种子点相位相同的种子点，包括：
[0021] 判断边界种子点的外推地震响应点位于河道上的位置；
[0022] 当判断边界种子点的外推地震响应点位于河道内部时，以外推地震响应点为新的 种子点，向同一主测线上、联络测线增大或减少的方向自动追踪相位相同的外推地震响应 点，直到判断外推地震响应点到了河道边界处，停止追踪，得到相邻主测线上的种子点。
[0023] 在一个实施例中，根据当前主测线上边界种子点的相干和倾角，驱动边界种子点 向相邻主测线自动追踪，得到相邻主测线上的、与边界种子点相位相同的种子点，包括： [0024] 当判断边界种子点的外推地震响应点位于河道边界或外部时，确定边界种子点为 无效种子点，重新确定有效种子点，以有效种子点为边界种子点，向相邻主测线自动追踪相 位相同的外推地震响应点，得到相邻主测线上的种子点；有效种子点为：与最后确定的无 效种子点相邻、且其外推地震响应点位于河道内部。
[0025] 在一个实施例中，判断外推地震响应点到了河道边界处，具体为：
[0026] 当外推地震响应点和与外推地震响应点相邻的8个地震响应点，共9个地震响应 点的相干平均值小，方差值大时，证明到了河道边界。
[0027] 本发明实施例技术方案基于地震波相似性原理，在设置河床种子点的基础上，以 地层倾角和地震波形相似性驱动种子点外推，同时采用多窗口旋转扫描识别河床边界，河 床边界约束外推范围，最终实现多条多期河道的侵蚀深度预测。因此，通过本发明实施例提 供的古河道河床侵蚀深度自动识别方法，可以间接指导暗