[发明专利]提高地震分辨率的保构造井控空间反褶积方法和装置

权利要求书

1.一种提高地震分辨率的保构造井控空间反褶积方法，该方法包括：

步骤101，基于测井曲线得到第一反射系数矩阵；

步骤102，给出初始的空间正则化参数序列λ_i,i＝1,2,...,N，N为原始的多道地震记录中包含的总道数，进入步骤103；

步骤103，采用沿空间方向的正则化，基于当前的空间正则化参数序列对所述原始的多道地震记录进行空间反褶积，得到第二反射系数矩阵，进入步骤104；

步骤104，如果得到的第一反射系数矩阵和第二反射系数矩阵间的误差符合预定期望，则进入步骤106，否则，进入步骤105；

步骤105，调整当前的空间正则化参数序列，然后回到步骤103；

步骤106，将当前得到的第二反射系数矩阵作为最终的反射系数矩阵，将其中每一道对应的反射系数序列与该道的子波进行褶积，并叠加多道的褶积结果，以得到所需的叠后地震数据体。

2.根据权利要求1所述的方法，其中，在步骤103中，得到第二反射系数矩阵包括：

设定目标函数E(m)：

其中，d是对原始的多道地震记录进行叠加后得到的原始叠后地震数据体，G为子波褶积矩阵；m表示待求的第二反射系数矩阵，λ_i为第i道对应的空间正则化参数，D_i为第i道对应的表示空间结构的频率-波数域算子，||*||₁表示1范数，||*||₂表示2范数；

迭代求解上式得到使得目标函数E(m)小于第一预设阈值的矩阵m，该矩阵m即为在当前的空间正则化参数序列下得到的第二反射系数矩阵。

3.根据权利要求1所述的方法，其中，在步骤104中，通过下列方法判断第一反射系数矩阵和第二反射系数矩阵间的误差是否符合预定期望：

设定目标函数min(m)：

其中，m表示第二反射系数矩阵，S表示第一反射系数矩阵，C表示协方差矩阵，λ_i为第i道对应的空间正则化参数，D_i为第i道对应的表示空间结构的频率-波数域算子，||*||₁表示1范数，||*||_p表示p范数；

如果min(m)小于第二预设阈值，则判断第一反射系数矩阵和第二反射系数矩阵间的误差符合预定期望，反之，则判断第一反射系数矩阵和第二反射系数矩阵间的误差不符合预定期望。

4.一种提高地震分辨率的保构造井控空间反褶积装置，该装置包括第一反射系数矩阵确定单元、初始化单元、第二反射系数矩阵确定单元、误差判断单元、空间正则化参数调整单元、褶积叠加单元，其中：

第一反射系数矩阵确定单元，用于基于测井曲线得到第一反射系数矩阵；

初始化单元，用于给出初始的空间正则化参数序列λ_i,i＝1,2,...,N，其中N为原始的多道地震记录中包含的总道数；

第二反射系数矩阵确定单元，采用沿空间方向的正则化，用于基于当前的空间正则化参数序列对所述原始的多道地震记录进行空间反褶积以得到第二反射系数矩阵，开始时从初始化单元接收空间正则化参数序列，在误差判断单元判断误差不符合预定期望的情况下，从空间正则化参数调整单元接收空间正则化参数序列；

误差判断单元，分别从第一反射系数矩阵确定单元和第二反射系数矩阵确定单元接收第一反射系数矩阵和第二反射系数矩阵，用于判断第一反射系数矩阵和第二反射系数矩阵间的误差是否符合预定期望；

空间正则化参数调整单元，在误差判断单元判断所述误差不符合预定期望的情况下，用于调整当前的空间正则化参数序列；

褶积叠加单元，在误差判断单元判断所述误差符合预定期望的情况下，用于将第二反射系数矩阵确定单元当前得到的第二反射系数矩阵作为最终的反射系数矩阵，并将其中每一道对应的反射系数序列与该道的子波进行褶积，以及叠加多道的褶积结果，以得到所需的叠后地震数据体。