[发明专利]一种地震子波的提取方法及处理终端

摘要

本发明涉及一种地震子波的提取方法及处理终端，所述提取方法包括依次进行的以下步骤：步骤S1：获取地震道中直达波s_j(t)，并经过傅里叶变换获得s_j(ω)；步骤S2：计算得出反射系数R谱；步骤S3：获取检波器到震源和虚震源的距离，计算出中间变量B；步骤S4：将步骤S3中的中间变量B、步骤S1中的s_j(ω)和步骤S2中的反射系数R代入公式计算出地震子波w(ω)。本发明不依赖于地震褶积理论模型，不需要测井资料，不需要对地震数据和地层反射系数的分布做任何假设，也不需要对地震子波的相位特性做任何假设；同时，本发明适用于海上拖缆地震资料的地震子波提取，计算结果精度高，可以支持超长序列地震子波的提取，满足地震资料批量处理的要求。

主权项

1.一种地震子波的提取方法，其特征在于：其包括依次进行的以下步骤：步骤S1：获取地震道中直达波s_j(t)，并经过傅里叶变换获得s_j(ω)；步骤S2：按公式①计算得出反射系数R：式中，R₀表示水平海面反射系数，为常数，射线参数θ为入射波的入射角，c表示水中的声波速度，为常数，ω表示入射波的圆周频率，σ表示海面波浪谱；步骤S3：获取检波器到震源和虚震源的距离，并代入公式②，计算出中间变量B：式中，l表示每个地震道检波器的组合个数，n表示地震道记录的个数，i为虚数单位，ω表示入射波的圆周频率，r_jk表示震源到第j个地震道中第k个检波器的直达波传播距离，r′_jk表示虚震源到第j个地震道中第k个检波器的直达波传播距离；步骤S4：将步骤S3中的中间变量B和步骤S1中的s_j(ω)代入公式③，计算出地震子波w(ω)：经过公式③即计算出地震子波w(ω)。