[发明专利]地震数据拓频方法及系统

权利要求书

1.一种地震数据拓频方法，其特征在于，包括：

根据原始地震数据获取一级降频地震波形、二级降频地震波形、一级升频地震波形和二级升频地震波形；

根据所述一级降频地震波形、所述二级降频地震波形、所述一级升频地震波形和所述二级升频地震波形计算四分频合成地震数据；

根据预设的四分频合成地震数据配比权重、预设的原始地震数据配比权重、所述原始地震数据和所述四分频合成地震数据计算多频融合地震数据；

根据原始地震数据获取一级降频地震波形、二级降频地震波形、一级升频地震波形和二级升频地震波形，包括：

对所述原始地震数据的初始极值特征点序列进行极值特征点筛选，得到一级降频极值特征点序列；

对所述一级降频极值特征点序列的每对相邻极值特征点之间进行余弦函数插值，获取一级降频地震波形；

对所述一级降频极值特征点序列进行极值特征点筛选，得到二级降频极值特征点序列；

对所述二级降频极值特征点序列的每对相邻极值特征点之间进行余弦函数插值，获取二级降频地震波形；

提取所述原始地震数据的地震中频成分，获取一级升频地震波形；

对所述原始地震数据的地震高频成分进行高频重构，获取二级升频地震波形；

根据所述一级降频地震波形、所述二级降频地震波形、所述一级升频地震波形和所述二级升频地震波形计算四分频合成地震数据，包括：

根据一级降频地震波形的高截频、一级降频地震波形的低截频、二级升频地震波形的高截频和二级降频地震波形的低截频计算一级降频地震波形配比权重；

根据二级降频地震波形的高截频、二级降频地震波形的低截频和二级升频地震波形的高截频计算二级降频地震波形配比权重；

根据一级升频地震波形的高截频、一级升频地震波形的低截频、二级升频地震波形的高截频和二级降频地震波形的低截频计算一级升频地震波形配比权重；

根据二级升频地震波形的高截频、二级升频地震波形的低截频和二级降频地震波形的低截频计算二级升频地震波形配比权重；

根据所述一级降频地震波形、所述二级降频地震波形、所述一级升频地震波形、所述二级升频地震波形、所述一级降频地震波形配比权重、所述二级降频地震波形配比权重、所述一级升频地震波形配比权重和所述二级升频地震波形配比权重计算所述四分频合成地震数据。

2.根据权利要求1所述的地震数据拓频方法，其特征在于，

通过如下公式计算一级降频地震波形配比权重：

通过如下公式计算二级降频地震波形配比权重：

通过如下公式计算一级升频地震波形配比权重：

通过如下公式计算二级升频地震波形配比权重：

其中，w_d1为一级降频地震波形配比权重，w_d2二级降频地震波形配比权重，w_u1为一级升频地震波形配比权重，w_u2为二级升频地震波形配比权重，为一级降频地震波形的高截频，为一级降频地震波形的低截频，为二级降频地震波形的高截频，为二级降频地震波形的低截频，为一级升频地震波形的高截频，为一级升频地震波形的低截频，为二级升频地震波形的高截频，为二级升频地震波形的低截频。

3.根据权利要求2所述的地震数据拓频方法，其特征在于，还包括：

通过如下公式计算四分频合成地震数据：

xs4_k＝w_d2*xd2_k+w_d1*xd1_k+w_u1*xu1_k+w_u2*xu2_k；k＝1,2,3,…nx；

其中，xs4_k为第k个四分频合成地震数据；nx为四分频合成地震数据的样点个数，xd1_k为第k个一级降频地震波形，xd2_k为第k个二级降频地震波形，xu1_k为第k个一级升频地震波形，xu2_k为第k个二级升频地震波形。