矢量凸集投影多分量三维地震数据重建方法及装置制造方法及图纸

本发明专利技术提供了一种矢量凸集投影多分量三维地震数据重建方法及装置，其中，该方法包括：将三或四分量三维地震数据的各分量数据变换至频率域；在频率域提取各分量正频率部分的单频率切片数据；将提取的正频率部分的单频率切片数据矢量联合表示为复四元数；基于复四元数傅里叶变换，利用矢量凸集投影迭代方法对单频率切片对应的复四元数矢量联合重建，得到正频率部分单频率切片的重建结果；对正频率部分单频率切片的重建结果中的各分量数据取共轭，得到负频率部分单频率切片的重建结果；根据正频率部分和相应负频率部分的重建结果，得到频率域所有频率切片的重建结果，并变换至时间域。通过上述方案能实现三或四分量数据联合重建而且能提高重建效率。而且能提高重建效率。而且能提高重建效率。

【技术实现步骤摘要】
矢量凸集投影多分量三维地震数据重建方法及装置


[0001]本专利技术涉及地质勘探
，尤其涉及一种矢量凸集投影多分量三维地震数据 重建方法及装置。

技术介绍

[0002]在陆上地震数据的采集阶段，受障碍物、采集禁区及采集成本因素的影响，采集 到的地震数据通常存在稀疏或者不规则缺失问题。在海上地震数据的采集阶段，采用 海洋拖缆(Tower Streamer，简称TC)采集地震数据时会受到羽状漂移影响，采用海 底电缆(Ocean Bottom Cable，简称OBC)采集数据时会受海流和潮汐的影响，使用 海底节点地震仪(Ocean Bottom Node，简称OBN)采集数据时会受到海底崎岖复杂 地形和仪器成本因素的影响，这些影响使得采集到的地震数据也存在稀疏采样或者不 规则采样问题。另外，在室内地震资料处理阶段，由于剔除废炮和废道，以及数据处 理域转换也会引起地震道的不规则缺失。不规则缺失地震数据会对常规多道处理技术 的正确运行产生不良影响，降低地震资料品质，因此，需要对不规则地震数据进行规 则化重建处理。
[0003]此外，随着油气本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种矢量凸集投影多分量三维地震数据重建方法，其特征在于，包括：对时间域的多分量三维地震数据的各分量数据做关于时间变量的普通傅里叶变换，得到频率域的多分量三维地震数据的各分量数据；其中，所述多分量为三分量或四分量；从频率域的多分量三维地震数据的各分量数据中分别提取正频率部分的频率切片数据；利用复四元数框架将从各分量数据中提取的正频率部分的单频率切片数据进行矢量联合表示，得到相应正频率部分的单频率切片对应的复四元数；基于复四元数傅里叶变换，利用矢量凸集投影迭代方法对单频率切片对应的复四元数进行矢量联合重建，得到相应正频率部分的单频率切片的多分量三维地震数据的重建结果；对正频率部分的单频率切片的多分量三维地震数据的重建结果中的各分量数据取共轭，得到相应负频率部分的单频率切片的多分量三维地震数据重建结果；根据正频率部分的单频率切片的多分量三维地震数据的重建结果和相应负频率部分的单频率切片的多分量三维地震数据的重建结果，得到整个频率域的多分量三维地震数据的重建结果；对整个频率域的多分量三维地震数据的重建结果中的各分量数据分别做普通傅里叶反变换，得到时间域的多分量三维地震数据的重建结果。2.如权利要求1所述的矢量凸集投影多分量三维地震数据重建方法，其特征在于，利用复四元数框架将从各分量数据提取的正频率部分的单频率切片数据进行矢量联合表示，得到相应正频率部分的单频率切片对应的复四元数，包括：在多分量为三分量的情况下，将从各分量数据提取的正频率部分的单频率切片数据中的每个分量作为复四元数框架中的一个虚部，得到相应正频率部分的单频率切片对应的复四元数。3.如权利要求1所述的矢量凸集投影多分量三维地震数据重建方法，其特征在于，利用复四元数框架将从各分量数据提取的正频率部分的单频率切片数据进行矢量联合表示，得到相应正频率部分的单频率切片对应的复四元数，包括：在多分量为四分量的情况下，将从各分量数据提取的正频率部分的单频率切片数据中任选其中一个分量作为复四元数框架中的实部，其余三个分量中的每个分量作为复四元数框架中的一个虚部，得到相应正频率部分的单频率切片对应的复四元数。4.如权利要求1所述的矢量凸集投影多分量三维地震数据重建方法，其特征在于，基于复四元数傅里叶变换，利用矢量凸集投影迭代方法对单频率切片对应的复四元数进行矢量联合重建，得到相应正频率部分的单频率切片的多分量三维地震数据的重建结果，包括：在多分量三维地震数据的各分量数据的采样算子相同的情况下，基于包含复四元数傅里叶正变换和复四元数傅里叶反变换的第一重建表达式，利用矢量凸集投影迭代方法对单频率切片对应的复四元数进行矢量联合重建，得到相应正频率部分的单频率切片的多分量三维地震数据的重建结果的矢量联合表示，从重建结果的矢量联合表示中提取得到重建结果的各分量数据；第一重建表达式为：
其中，Q
l
(f
n
,x1,x2)表示单频率切片对应的复四元数的第l次重建结果，l＝1,2,...,N
iter
，N
iter
表示最大迭代次数，f
n
表示第n个频率切片，x1和x2表示检波器在地面坐标的两个空间变量，α表示加权因子，Q
obs
(f
n
,x1,x2)表示单频率切片对应的初始复四元数，I表示全1矩阵，S表示采样算子，表示复四元数傅里叶反变换，T
l
表示第l次重建的阈值算子，F
Q
表示复四元数傅里叶变换，Q
l
‑1(f
n
,x1,x2)表示单频率切片对应的复四元数的第l
‑
1次重建结果，Q0(f
n
,x1,x2)＝Q
obs
(f
n
,x1,x2)。5.如权利要求1所述的矢量凸集投影多分量三维地震数据重建方法，其特征在于，基于复四元数傅里叶变换，利用矢量凸集投影迭代方法对单频率切片对应的复四元数进行矢量联合重建，得到相应正频率部分的单频率切片的多分量三维地震数据的重建结果，包括：在多分量三维地震数据的各分量数据的采样算子存在不同的情况下，基于包含复四元数傅里叶正变换和复四元数傅里叶反变换的第二重建表达式，利用矢量凸集投影迭代方法对单频率切片对应的复四元数进行矢量联合重建，得到相应正频率部分的单频率切片的多分量三维...

【专利技术属性】
技术研发人员：高建军，李帆，王赟，芦俊，
申请(专利权)人：中国地质大学北京，
类型：发明
国别省市：