[发明专利]一种油气储层中的可控层位分阶层间多次波模拟方法

摘要

本发明公布了一种油气储层中的可控层位分阶层间多次波模拟方法，属于地震波场正演模拟领域，包括自适应延拓步长计算、二维反周期延拓边界反射处理和波场模拟中的双重层位约束。该方法基于自适应变步长波场延拓，以递归循环的方式实现分阶层间多次波的模拟；包括输入数据准备过程、递归循环波场模拟过程和压制边界反射过程。通过对模型添加双重层位约束，可以模拟指定地层产生的各阶层间多次波；利用二维反周期延拓方法压制波场延拓过程中产生的人工边界反射，得到了比传统衰减边界更佳的边界反射压制效果；提出自适应变步长波场延拓技术，大大提升了波场模拟的效率。

主权项

1.一种油气储层中的可控层位分阶层间多次波模拟方法，通过构建层位约束矩阵、自适应延拓步长计算、震源设置、递归循环波场模拟和边界反射压制，实现油气储层中指定地层之间的各阶层间多次波模拟；包括以下步骤：A.数据准备，执行如下操作：A1.构建层位约束矩阵：在地质模型中拾取发生下行反射的层位和发生上行反射的层位，用于构建下行反射层位约束矩阵和上行反射层位约束矩阵n为层位位置；d和u分别表示下行和上行；A2.自适应延拓步长计算，即计算地质模型随深度变化的延拓步长；A3.震源设置：根据时间采样率、时间记录长度、检波点数、震源位置及震源子波定义原始震源矩阵s，对s进行二维反周期延拓，形成反周期延拓后的反周期延拓震源矩阵B.递归循环波场模拟，执行如下操作：B1.初始化：将原始震源向量赋值给地表的下行波场向量并将地下各个层位下方的上行波场向量赋值为零向量；B2.将地表的下行波场向量根据步骤A2计算得到的自适应延拓步长逐步向下延拓至地质模型底部，得到各个层位上方的下行波场向量和层位下方的下行波场向量及模型底部层位上方的下行波场层位e表示模型最深部的层位；B3.将B2中得到的模型底部上方的下行波场乘以模型底部层位的下行反射系数矩阵得到上行波场B4.将模型底部层位上方的上行波场逐步向上延拓至地表，得到各个层位上方的上行波场向量和层位下方的上行波场向量及地表层位下方的上行波场为即模拟的采集数据；B5.递归执行B2‑B4，每循环一次得到更高一阶的层间多次波波场，直到得到所需阶次的层间多次波波场；B6.反周期延拓波场模拟，用反周期延拓震源替换原始震源执行B2‑B5，模拟得到二维反周期延拓后的层间多次波波场；C.压制边界反射：将B5和B6中对应阶次的层间多次波波场求平均，得到压制边界反射后的层间多次波波场；通过上述步骤，实现油气储层中的可控层位分阶层间多次波模拟。