批量自动井震标定的方法

技术领域

本发明涉及油气田开发和油藏描述技术领域，特别是涉及到一种批量自动井震标定的方法。

背景技术

井震标定是油藏描述与开发中的一项关键技术，在地震资料解释、储层预测等过程中起着关键作用。地震资料的特征是具有较好的横向连续性但纵向分辨率较低，测井资料虽然有较高的纵向精度，但不具备横向展布特征。随着勘探的深入，依靠单一资料解决勘探问题十分困难，因此通过井震标定建立时间域地震资料与深度域测井资料的桥梁。

常规的井震标定流程是通过声波测井资料与密度测井资料求得反射系数序列并与子波褶积制作合成记录，通过手动调整与井旁地震道完成对齐，该方法在井较少时仍具有一定可行性。以胜利油田为代表的一批老油田特点是井数量多且密，井震标定耗时较长且标定的精度受人为因素影响较大。

在申请号：CN202210015020.0的中国专利申请中，涉及到一种基于多道地震叠加的井震标定方法，所述方法包括：获得待标定井的时间域合成记录；利用自动翘曲算法获得一般井旁地震道与参考井旁地震道的时移量，并使用该时移量对各井旁地震道进行校正；将校正后的井旁地震道与所述参考井旁地震道进行相关系数加权叠加，形成初始复合井旁地震道；对所述初始复合井旁道进行有效信号增强处理和/或标准化处理，得到处理后复合井旁道；将所述时间域合成记录与所述处理后复合井旁道进行对比处理，实现井震标定。本发明可在地震资料信噪比较低时，显著提升井震标定结果的准确性及合理性。

在申请号：CN202010115378.1的中国专利申请中，涉及到一种基于多子波的井震标定方法及装置，该方法包括：获取目标井位置处采集的测井数据和井旁地震道数据；根据地震数据纵向分辨率的不同，将井旁地震道数据划分为多个时窗；在不同的时窗内，采用不同频率的子波对测井数据进行井震标定，生成测井数据对应的合成地震记录；对不同时窗内的子波进行相位调整，使得测井数据对应的合成地震记录与井旁地震道数据满足预设相关性条件。本发明采用不同频率不同相位的地震子波进行井震标定，能够使得生成的合成地震记录与井旁地震道数据具有良好的匹配，从而实现更准确的层位标定。

在申请号：CN201910970415.4的中国专利申请中，涉及到一种井震标定方法，它包括:一、地震数据与测井数据准备；二、基于地震解释方案选取顶底标准层；三、顶部标准层标定与质控；四、底部标准层标定与质控；五、基于顶底标准层标定结果建立速度场；六、生成合成地震记录，开展井震层间小层标定。本发明以顶底参考层为约束，层间采用声波速度代替声波转换合成地震记录速度，利用平均速度拟合和趋势面去低频效应等手段消除空间误差，获得顶底及标准层间精确速度，本发明能够弥补井信息与地震数据存在时间和波形上不匹配问题，减小横向非均质性对井标定速度外推的影响，对油气资源勘探具有重要应用价值。

以上现有技术均与本发明有较大区别，未能解决我们想要解决的技术问题，为此我们发明了一种新的批量自动井震标定的方法。

发明内容

本发明的目的是提供一种能够在保证标定精度的情况下，有效提高井震标定的效率的批量自动井震标定的方法。

本发明的目的可通过如下技术措施来实现：批量自动井震标定的方法，该批量自动井震标定的方法包括：

步骤1，选择资料品质较好的单井作为参考井，制作完成单井时间域合成记录；

步骤2，抽取井旁地震道，利用地震层位与测井分层的绑定关系或初始时间关系对合成记录进行移动对齐，然后用动态时间规整算法计算合成记录与井旁地震道的新匹配关系，完成拉伸压缩，实现单井的标定；

步骤3，利用参考井参数，经过同样的步骤1和2完成其余井的合成记录制作及标定。

本发明的目的还可通过如下技术措施来实现：

在步骤1中，选择资料品质好的单井作为参考井，利用经过预处理的声波时差测井曲线与密度测井曲线相乘得到反射系数序列，将反射系数序列与子波做褶积，制作完成单井时间域合成记录。

步骤1包括：

S11：对待标定的井进行测井资料预处理，利用单井的声波时差测井曲线与密度测井曲线计算得到反射系数序列；

S12：将深度域反射系数序列转换到时间域，选择合适的地震子波与时间序列褶积。

在步骤S11中，计算反射系数序列的公式为：

r(i)＝vel(i)*den(i)

其中，r(i)表示深度域反射系数序列；vel(i)表示由声波时差测井曲线计算所得的速度序列；den(i)表示测井密度曲线；i表示深度域采样点。

在步骤S12中，地震子波与时间序列褶积的公式为：

W(t)＝r(t)*ω(t)

其中，W(t)表示合成记录序列；ω(t)表示地震子波序列；r

步骤2包括：

S21：根据井点位置从地震数据中抽取井旁单道地震道S(t)；

S22：利用地震层位与井分层绑定或者已有的初始时深关系对制作的合成记录进行移动对齐；

S23：通过井旁地震道与合成记录建立误差矩阵；

S24：利用改进的动态路径规划实现对误差矩阵的最小累计，建立累计最小误差矩阵；

S25：利用反向的动态路径规划进行路径回溯，寻找序列u；

S26：将序列u中时移量施加至原始合成记录中，得到新的合成记录

在步骤S23中，建立的误差矩阵为：

d(t,l)＝|S(t)-W(t+l)|

其中，d(t,l)表示由两序列建立的误差矩阵；l表示移动量，取值范围为-L至L；L为时窗，人为给定；W(t)表示合成记录序列。

在步骤S24中，建立的累计最小误差矩阵为：

D(0,l)＝d(0,l)

其中，D(t,l)表示累计最小误差矩阵；t表示时间域的采样点。

在步骤S25中，序列u为：

其中，u(t)表示合成记录每个采样点标定时的时移量，

在步骤S26中，新的合成记录

其中，

步骤3包括：

S31：采用与参考井相同的子波参数批量制作合成记录；

S32：利用地震层位与井分层绑定或者已有的初始时深关系对制作的合成记录进行批量移动对齐；

S33：利用动态时间规整完成对剩余的所有合成记录的拉伸压缩。

本发明中的批量自动井震标定的方法，满足了老油田高密度井网条件下标定的挑战，能够在保证标定精度的情况下，有效提高井震标定的效率。

有益效果如下：

本发明提出的改进的动态时间规整算法，可以自动实现批量井震标定，能够提高标定的精度、效率，特别是针对密井网的情况下，大大提高井震标定效率，从而提高后续的反演精度。

附图说明

图1为本发明的批量自动井震标定的方法的一具体实施例的流程图；

图2为本发明的一具体实施例中参考井合成记录制作及通过层位绑定移动对齐图；

图3为本发明的一具体实施例中参考井合成记录与井旁道建立的误差矩阵的示意图；

图4为本发明的一具体实施例中参考井合成记录与井旁道建立的累计最小误差矩阵及路径回溯的示意图；

图5为本发明的一具体实施例中参考井通过该方法自动标定的结果的示意图；

图6为本发明的一具体实施例中剩余井批量自动标定的结果的示意图。

具体实施方式

应该指出，以下详细说明都是示例性的，旨在对本发明提供进一步的说明。除非另有指明，本文使用的所有技术和科学术语具有与本发明所属技术领域的普通技术人员通常理解的相同含义。

需要注意的是，这里所使用的术语仅是为了描述具体实施方式，而非意图限制根据本发明的示例性实施方式。如在这里所使用的，除非上下文另外明确指出，否则单数形式也意图包括复数形式，此外，还应当理解的是，当在本说明书中使用术语“包含”和/或“包括”时，其指明存在特征、步骤、操作和/或它们的组合。

本发明的批量自动井震标定的方法包括：选择资料品质好的单井作为参考井，利用经过预处理的声波时差测井曲线与密度测井曲线相乘得到反射系数序列，将反射系数序列与子波做褶积，制作完成单井时间域合成记录；抽取井旁地震道，利用地震层位与测井分层的绑定关系或初始时间关系对合成记录进行移动对齐，然后用动态时间规整算法计算合成记录与井旁地震道的新匹配关系，完成拉伸压缩，实现单井的标定；利用参考井参数，完成其余井的合成记录制作及标定。本发明提出的改进的动态时间规整算法，自动实现批量井震标定，能够提高标定的精度、效率。

以下为应用本发明的几个具体实施例

实施例1

在应用本发明的一具体实施例1中，该批量自动井震标定的方法包括了以下步骤：

S1：选择资料品质好的单井作为参考井，利用经过预处理的声波时差测井曲线与密度测井曲线相乘得到反射系数序列，将反射系数序列与子波做褶积，制作完成单井时间域合成记录；

S1包括如下子步骤：

S11：对待标定的井进行测井资料预处理，利用单井的声波时差测井曲线与密度测井曲线计算得到反射系数序列：

r(i)＝vel(i)*den(i)

其中，r(i)表示深度域反射系数序列；vel(i)表示由声波时差测井曲线计算所得的速度序列；den(i)表示测井密度曲线；i表示深度域采样点。

S12：将深度域反射系数序列转换到时间域，选择合适的地震子波与时间序列褶积,

W(t)＝r(t)*ω(t)

其中，W(t)表示合成记录序列；ω(t)表示地震子波序列；r

S2：抽取井旁地震道，利用地震层位与测井分层的绑定关系或初始时间关系对合成记录进行移动对齐，然后用动态时间规整算法计算合成记录与井旁地震道的新匹配关系，完成拉伸压缩，实现单井的标定；

S2包括如下子步骤：

S21：根据井点位置从地震数据中抽取井旁单道地震道S(t)；

S22：利用地震层位与井分层绑定或者已有的初始时深关系对制作的合成记录进行移动对齐；

S23：通过井旁地震道与合成记录建立误差矩阵：

d(t,l)＝|S(t)-W(t+l)|

其中，d(t,l)表示由两序列建立的误差矩阵；l表示移动量，取值范围为-L至L；L为时窗，人为给定。

S24：利用改进的动态路径规划实现对误差矩阵的最小累计，建立累计最小误差矩阵：

D(0,l)＝d(0,l)

其中，D(t,l)表示累计最小误差矩阵；

S25：利用反向的动态路径规划进行路径回溯，寻找序列u：

其中，u(t)表示合成记录每个采样点标定时的时移量；

S26：将序列u中时移量施加至原始合成记录中，得到新的合成记录

其中，

S3：利用参考井参数，经过步骤S1和S2完成其余井的合成记录制作及标定。

S31：采用与参考井相同的子波参数批量制作合成记录；

S32：利用地震层位与井分层绑定或者已有的初始时深关系对制作的合成记录进行批量移动对齐；

S33：利用动态时间规整完成对剩余的所有合成记录的拉伸压缩。

实施例2

在应用本发明的一具体实施例中，以胜利油田某工区为例子，该区块共有井分层的直井700余口，如图1所示，图1为本发明的批量自动井震标定的方法的流程图。研究批量标定自动标定方法分为四步：

第一步，选取井资料较好的参考井A，利用将该井声波时差测井曲线和密度测井曲线进行运算得到反射系数序列，再将反射系数序列与子波进行褶积运算，得到该井时间域合成记录，通过地震层位与测井层位绑定，将其移动对齐标定，如图2所示。

第二步，提取该参考井A的井旁道道，通过专家经验给出时窗L的值，L＝5，建立A井的合成记录与井旁道的误差矩阵，如图3所示。利用改进的动态路径规划实现对误差矩阵的最小累计，并建立累计最小误差矩阵，并通过反向的路径规划寻找井A的序列u，如图4所示。

第三步，将A井的序列u中时移量施加至原始合成记录中，得到该井新的合成记录，完成参考井A的井震标定，如图5所示。

第四步，采用第一步相同的方法及与参考井A相同的子波参数批量制作剩余井的合成记录，并利用第二、三步同样的步骤及子波参数、时窗L值完成工区其余井的快速批量自动标定，如图6所示。

实施例3

在应用本发明的一具体实施例中，以孤岛工区某区块为例子，该区块共有的直井300余口，研究批量标定自动标定方法分为四步：

第一步，选取井资料较好的参考井B，利用参考井B的声波测井曲线和密度测井曲线进行运算得到反射系数序列，再将反射系数序列与子波进行褶积运算，得到该井时间域合成记录，通过地震层位与测井层位绑定，将其移动对齐。

第二步，提取该参考井B的井旁地震道，给出时窗L的值,L＝10，建立该参考井的合成记录与井旁道的误差矩阵。利用改进的动态路径规划实现对误差矩阵的最小累计，并建立累计最小误差矩阵，并通过反向的路径规划寻找参考井的序列u。

第三步，将A井的序列u中时移量施加至原始合成记录中，得到该井新的合成记录，完成参考井b的井震标定。

第四步，采用第一步相同的方法及与参考井B相同的子波参数批量制作剩余井的合成记录，并利用第二、三步同样的步骤及子波参数、时窗L值完成工区其余井的快速批量自动标定。

最后应说明的是：以上所述仅为本发明的优选实施例而已，并不用于限制本发明，尽管参照前述实施例对本发明进行了详细的说明，对于本领域技术人员来说，其依然可以对前述实施例记载的技术方案进行修改，或者对其中部分技术特征进行等同替换。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

除说明书所述的技术特征外，均为本专业技术人员的已知技术。