[发明专利]DDW高精度深度域井震匹配方法

权利要求书

1.DDW高精度深度域井震匹配方法，其特征在于，该DDW高精度深度域井震匹配方法包括：

步骤1，输入地震数据与地震地质层位，并对数据进行预处理；

步骤2，选择目标井的井旁地震道X，进行深度域合成记录制作得到Y，参照地质分层进行合成记录样点选取；

步骤3，在地质分层约束下扫描井旁道；

步骤4，逐点进行距离计算，寻找与合成记录样点距离最短的点；

步骤5，求取累计距离，规划最短路径；

步骤6，建立逐点对应的井震匹配关系，计算井震误差校正量；

步骤7，建立偏差约束体校正地震数据；

在步骤5中，采用基于递推的动态规划的方法求出该累计距离，从后向前递推，若(N,M)是一个匹配点，即地震深度N对应测井深度M，则该路径点前一个格点可能是(N-1,M)，(N-1,M-1)，(N,M-1)，找出这三个点中距离最小者即为其上一最优路径点，其累计距离公式为

TD(i,j)＝D(i,j)+min{D(i-1,j)，D(i-1,j-1),D(i,j-1)} (式2)

其中，TD(i,j)为某一路径上的累计距离，D(i,j)表示两记录中相应采样点i与j的距离；从最后一点或某个确定点出发搜索，反复递推，直到到达地震深度的起始点，构建出一条累计距离最小的路径，该路径即为两信号匹配的最优路径；在距离计算上定义为

D(i，j)＝A₁*|R(i)-S(j)|+A₂*|F(R(i))-F(S(j))|+A₃*|T(R(i))-T(S(j))| (式3)

其中，A₁、A₂、A₃为加权因子，F(i)为主频，T(i)为时频，R(i)与S(j)分别为两列信号各个元素的振幅值，D(i,j)为R(i)与S(j)的振幅值之差，即两帧信号采样点之间的距离。

2.根据权利要求1所述的DDW高精度深度域井震匹配方法，其特征在于，在步骤1中，输入的数据包括：深度域地震数据、测井数据、地震解释层位、地质分层数据。

3.根据权利要求2所述的DDW高精度深度域井震匹配方法，其特征在于，在步骤1中，进行深度域地震数据预处理，对深度域地震数据开展高阶累计量时频分析和频谱分析，获取相应的频谱类参数数据体，提取地震子波。

4.根据权利要求2所述的DDW高精度深度域井震匹配方法，其特征在于，在步骤1中，进行深度域测井数据预处理，提取声波时差、密度这些测井曲线，完成拼接、测井环境校正这些曲线预处理。

5.根据权利要求1所述的DDW高精度深度域井震匹配方法，其特征在于，在步骤2中，根据地震数据和测井数据的位置信息自动将测井位置定位在相应的地震位置上，将测井的分层信息直接标在深度域地震剖面上，统计井点处地质分层与地震解释的深度偏差表，建立各个解释层位的深度偏差量版。

6.根据权利要求1所述的DDW高精度深度域井震匹配方法，其特征在于，在步骤2中，对地震子波按照声波时差及其对应深度关系开展深变处理，利用褶积原理完成深度域合成记录制作。

7.根据权利要求1所述的DDW高精度深度域井震匹配方法，其特征在于，在步骤4中，首先计算地震信号与合成记录信号这两个序列中每个采样点之间的距离，即相似度，振幅值差越小，距离越小，二者越相似；定义D(i,j)为R(i)与S(j)的振幅值之差，即两帧信号采样点之间的距离，表示为D(i,j)＝|R(i)-S(j)| (式1)R(i)与S(j)分别为两列信号各个元素的振幅值，这样D(i,j)构成一个ixj的矩阵。

8.根据权利要求1所述的DDW高精度深度域井震匹配方法，其特征在于，在步骤4中，进行的逐点距离计算包括波形幅度、频率、时频分布。