[发明专利]一种具有差异特征的地震数据融合的方法及处理终端

摘要

本发明涉及一种具有差异特征的地震数据融合的方法，包括依次进行的以下步骤：步骤S1：获取OBS数据，并对OBS数据进行预处理；步骤S2：根据预先从OBS数据的反射成像叠加剖面中选取的数据进行计算，获得平均地震子波；步骤S3：进行波阻抗反演，获得OBS波阻抗；步骤S4：获取拖缆地震数据，并对拖缆地震数据进行预处理；步骤S5：与步骤S2一样，经过计算后，获得拖缆地震数据的平均地震子波；步骤S6：进行波阻抗反演，获得地下界面的反射系数；步骤S7：进行一致性处理；骤S8：进行加权叠加生成新的振幅谱，并获得新的地震子波；步骤S9：进行褶积运算，获得融合后的结果。本发明实现OBS和拖缆数据的融合，丰富了低频信息，拓宽了频带，提高了成像精度。

主权项

1.一种具有差异特征的地震数据融合的方法，其特征在于：包括依次进行的以下步骤：步骤S1：获取OBS数据，并对OBS数据进行预处理，获得均方根速度和OBS数据的反射成像叠加剖面；步骤S2：根据预先从OBS数据的反射成像叠加剖面中选取的数据按公式①至④进行计算，获得OBS数据的平均地震子波和平均地震子波的振幅谱：式中，x(t)λ表示选取的数据中的第λ个地震道，ω为傅里叶变换变量，i为虚数单位，N表示地震道的个数，t为时间，X(ω)λ为x(t)λ的傅里叶变换，|X(ω)λ|为第λ个地震道的振幅谱，w(t)λ为第λ个单道地震子波，w(t)ave为OBS数据的平均地震子波，|X(ω)obs|为OBS数据的平均地震子波的振幅谱；步骤S3：以步骤S1获得的均方根速度作为初始模型、步骤S2中的平均地震子波和步骤S1预处理后的反射成像叠加剖面作为波阻抗反演需要的三个输入参数进行波阻抗反演，获得OBS波阻抗；步骤S4：获取拖缆地震数据，并对拖缆地震数据进行预处理，获得拖缆地震数据的反射成像叠加剖面；步骤S5：根据预先从拖缆地震数据的反射成像叠加剖面中选取的数据按步骤S2的公式①至④进行计算，获得拖缆地震数据的平均地震子波和拖缆地震数据的平均地震子波的振幅谱；步骤S6：以步骤S3中的OBS波阻抗作为初始模型、步骤S5中的拖缆地震数据的平均地震子波、步骤S4预处理后的反射成像叠加剖面作为波阻抗反演需要的三个输入参数进行波阻抗反演，获得地下界面的反射系数；步骤S7：对OBS数据的平均地震子波的振幅谱与拖缆地震数据的平均地震子波的振幅谱进行一致性处理，一致性处理依次按公式⑤和⑥进行：式中，max表示取最大值，coe为OBS振幅一致性处理系数，|X(ω)|_mcs为步骤S5中的拖缆地震数据的平均地震子波的振幅谱，为一致性处理后的OBS数据的平均地震子波的振幅谱；步骤S8：将步骤S7中的和|X(ω)|_mcs按公式⑦进行加权叠加并生成新的振幅谱：式中，α和β分别表示OBS振幅谱权重和拖缆振幅谱权重，|X(ω)|‑mcs表示新构建地震子波的振幅谱；经过公式⑦计算后，将|X(ω)|‑mcs代入公式⑧进行计算，并获得新的地震子波w′(t)：步骤S9：将步骤S8中获得的新的地震子波w′(t)与步骤S6中获得的反射系数进行褶积运算，获得OBS数据与拖缆地震数据融合后的结果，褶积运算如公式⑨：x′(t)＝w′(t)*r(t)‑‑‑‑‑‑⑨式中，x′(t)为融合后的地震数据，r(t)为步骤S6获得的反射系数。