[发明专利]一种从深度域地震数据中提取深度域地震子波的方法

权利要求书

1.一种从深度域地震数据中提取深度域地震子波的方法，其包括以下步骤：

1)从某口井测井数据中获得深度、速度和密度信息，由此计算出对应的反射系数r；

2)从叠前深度偏移地震数据中挑选这口井的井旁地震道集s，同时选取一个常速度v_c，并以该常速度为标准速度，对选取的井旁地震道集s进行速度、深度参数的变换，使其符合线性时不变的条件，得到变换后的井旁地震道集s'；

3)将反射系数r映射到变换后的深度位置上成为r'；

4)利用s'和r'提取深度域地震子波，按如下过程实现：

首先，将经过变换后的深度域地震褶积模型用向量的形式表达为：

式中，是由变换后的井旁地震道集s'构建的列向量，R是由深度变换后的反射系数r'构建的Toeplitz矩阵，是待求取的深度域地震子波向量；

为了得到上式中的唯一解，对下式无约束最小化问题进行优化：

式中，α为岭参数，α＞0；

在实际计算中，根据实际数据的特点确定合适的α值，利用下式来求解上式最小化问题的解析解

式中，R^T是矩阵R的转置矩阵，I是单位矩阵，(·)^-1表示对矩阵的求逆运算。