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

摘要

本发明实现了一种利用岭回归法，从深度域地震数据中提取深度域地震子波的方法，它包括以下步骤：1)从某口井的测井数据中获得深度、速度和密度信息，由此计算出对应的反射系数r；2)从叠前深度偏移地震数据中抽取该井的井旁地震道集s，并选取一个常速度v_c作为标准速度，对选取的井旁地震道集s进行速度、深度参数的变换，使其符合线性时不变条件，得到变换后的井旁地震道集s'；3)将反射系数r映射到变换后的深度位置上成为r'；4)选取合适的岭参数α，用岭回归方法提取地震子波。本发明的方法是直接利用深度域地震数据提取深度域地震子波，且提取的深度域地震子波精度高，而不需要对深度域地震数据进行时间至深度域、深度至时间域的相互转换。

主权项

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表示对矩阵的求逆运算。