【技术实现步骤摘要】
本专利技术涉及石油地球物理勘探资料处理，特别涉及一种多采样率地震数据体重构方法和装置。
技术介绍
1、我国东部地区由于城镇的快速发展，大量的生活和工农业设施严重影响了地震资料采集的质量；西部地区复杂的地表和地质构造严重影响地震激发点和接收点的布置。因此，目前常规的规则地震采集方法已无法保证高质量的地震信号采集，势必需要利用非规则地震采集稀疏优化的特点，采用多采样率采样方式经济高效地采集地震数据，而非规则地震数据重构方法是制约非规则地震采集的关键因素，目前尚未有成熟有效的三维数据重构方法进行实际工区应用。
2、基于稀疏变换的重构方法计算效率高，数值效果稳定，是非规则地震数据重构方法重要的研究方向。该方法利用信号在变换域的稀疏性为约束建立反演模型，并且采用稀疏反演方法求解。稀疏变换的性质是决定这类方法质量的首要因素。地震信号在变换域中的表达越稀疏，则越有利于提高重构质量。在地震信号处理领域经典的稀疏变换有傅里叶变换和radon类变换。最近曲波变换等多尺度变换被用于地震数据重构且展示了很好的数值效果。多尺度变换比单一尺度的变换能够展示更多的信息，因此是信号分析的有力工具。曲波变换打破了传统傅里叶变换和小波变换只能表示单一尺度或者不能稀疏表示多尺度信号的缺点，可对曲线或曲面形状信号进行多方向和各项异性的表示，为边缘不连续的目标提供了几乎最优的表达。张华等提出了对时间切片利用曲波变换进行重构进而实现三维地震数据重构；王本锋等在频率域对有效频率切片进行二维重构，最终实现三维地震数据重构。曹静杰等提出一种三维低冗余曲波变换的地震数
技术实现思路
1、专利技术人发现，虽然曲波变换有多尺度性、多方向性以及无须对数据划分时间窗口等优点，但是曲波变换的冗余度(冗余度指的是变换域系数的维数与原始信号系数维数之比)很大。在二维情况下的原始曲波变换的冗余度约为8，三维数据变换的冗余度为24～32，需要很多的内存和计算时间，影响大规模数据计算的效率。鉴于上述问题，提出了本专利技术以便提供一种克服上述问题或者至少部分地解决上述问题的一种多采样率地震数据体重构方法和装置，能够高保真高效率的完成多采样率地震数据体的重构。
2、第一方面，本专利技术实施例提供一种多采样率地震数据体重构方法，包括：
3、根据多采样率地震数据体建立观测系统规则的空数据体；
4、为所述空数据体中的道赋值所述地震数据体中匹配道的数据或设定初始值，构成初始数据体；
5、将所述初始数据体划分为多个小数据体，对每个小数据体作时间变量的一维傅里叶变换，得到f-x-y域的数据体；
6、针对每个f-x-y域的数据体，通过三维曲波迭代变换得到变换后数据体；
7、对变换后数据体进行组合得到初始重构数据体；
8、根据所述地震数据体修正所述初始重构数据体，得到重构数据体。
9、第二方面，本专利技术实施例提供一种多采样率地震数据体重构装置，包括：
10、空数据体建立模块，用于根据多采样率地震数据体建立观测系统规则的空数据体；
11、初始数据体建立模块，用于为所述空数据体中的道赋值所述地震数据体中匹配道的数据或设定初始值，构成初始数据体；
12、f-x-y域的数据体建立模块，用于将所述初始数据体划分为多个小数据体，对每个小数据体作时间变量的一维傅里叶变换，得到f-x-y域的数据体；
13、变换后数据体建立模块，用于针对每个f-x-y域的数据体，通过三维曲波迭代变换得到变换后数据体；
14、初始重构数据体建立模块，用于对变换后数据体进行组合得到初始重构数据体；
15、重构数据体建立模块，用于根据所述地震数据体修正所述初始重构数据体，得到重构数据体。
16、第三方面，本专利技术实施例提供一种计算机程序产品，包括计算机程序/指令，其中，该计算机程序/指令被处理器执行时实现上述多采样率地震数据体重构方法。
17、第四方面，本公开实施例提供一种服务器，包括：存储器、处理器及存储于存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现上述多采样率地震数据体重构方法。
18、本专利技术实施例提供的上述技术方案的有益效果至少包括：
19、(1)本专利技术实施例提供的多采样率地震数据体重构方法，首先由多采样率的地震数据体建立初始数据体；将初始数据体划分为多个小数据体，对每个小数据体作时间变量的一维傅里叶变换，得到f-x-y域的数据体；针对每个f-x-y域的数据体，通过曲波迭代变换得到变换后数据体；对变换后数据体进行组合得到初始重构数据体；由原始的多采样率地震数据体对初始重构数据体进行修正，得到最终的重构数据体。分线束在f-x-y域三维数据体进行三维曲波变换，解决了大数据体三维曲波变换占内存大的问题，同时分线束也为并行计算提供基础，提高了计算效率。
20、(2)本专利技术实施例提供的多采样率地震数据体重构方法，针对初始重构数据体中的每道，若该道为非待重构道，将该道数据替换为地震数据体中匹配道的数据，若该道为待重构道，根据地震数据体中相关道的数据修正该道数据，得到重构数据体。保证了重构数据的保真度。
21、本专利技术的其它特征和优点将在随后的说明书中阐述，并且，部分地从说明书中变得显而易见，或者通过实施本专利技术而了解。本专利技术的目的和其他优点可通过在所写的说明书、权利要求书、以及附图中所特别指出的结构来实现和获得。
22、下面通过附图和实施例，对本专利技术的技术方案做进一步的详细描述。
本文档来自技高网...【技术保护点】
1.一种多采样率地震数据体重构方法，其特征在于，包括：
2.如权利要求1所述的方法，其特征在于，所述根据多采样率地震数据体建立观测系统规则的空数据体，具体包括：
3.如权利要求1所述的方法，其特征在于，所述为所述空数据体中的道赋值所述地震数据体中匹配道的数据或设定初始值，具体包括：
4.如权利要求3所述的方法，其特征在于，所述距离阈值通过下述方式预先确定：
5.如权利要求1所述的方法，其特征在于，所述将所述初始数据体划分为多个小数据体，具体包括：
6.如权利要求5所述的方法，其特征在于，所述对变换后数据体进行组合得到初始重构数据体，具体包括：
7.如权利要求1所述的方法，其特征在于，所述针对每个f-x-y域的数据体，通过三维曲波迭代变换得到变换后数据体，具体包括针对每个f-x-y域的数据体执行下述步骤：
8.如权利要求7所述的方法，其特征在于，所述根据所述曲波系数集中的最大曲波系数和平均曲波系数确定迭代阈值的最大值和最小值，具体包括：
9.如权利要求7所述的方法，其特征在于，所述根据当
10.如权利要求7所述的方法，其特征在于，所述将当前曲波系数集中小于迭代阈值参数的曲波系数替换为设定值，具体包括：
11.如权利要求1所述的方法，其特征在于，为所述空数据体中的道赋值设定初始值后，还包括：
12.如权利要求11所述的方法，其特征在于，所述根据所述地震数据体中相关道的数据修正该道数据，具体包括：
13.如权利要求12所述的方法，其特征在于，所述根据该待重构道的位置确定搜索范围，具体包括：
14.如权利要求12所述的方法，其特征在于，所述判断确定的最大振幅绝对值与该待重构道的最大振幅绝对值的比值是否满足设定条件，具体包括：
15.一种多采样率地震数据体重构装置，其特征在于，包括：
16.一种计算机程序产品，包括计算机程序/指令，其特征在于，该计算机程序/指令被处理器执行时实现权利要求1～15中任一项所述的多采样率地震数据体重构方法。
...【技术特征摘要】
1.一种多采样率地震数据体重构方法，其特征在于，包括：
2.如权利要求1所述的方法，其特征在于，所述根据多采样率地震数据体建立观测系统规则的空数据体，具体包括：
3.如权利要求1所述的方法，其特征在于，所述为所述空数据体中的道赋值所述地震数据体中匹配道的数据或设定初始值，具体包括：
4.如权利要求3所述的方法，其特征在于，所述距离阈值通过下述方式预先确定：
5.如权利要求1所述的方法，其特征在于，所述将所述初始数据体划分为多个小数据体，具体包括：
6.如权利要求5所述的方法，其特征在于，所述对变换后数据体进行组合得到初始重构数据体，具体包括：
7.如权利要求1所述的方法，其特征在于，所述针对每个f-x-y域的数据体，通过三维曲波迭代变换得到变换后数据体，具体包括针对每个f-x-y域的数据体执行下述步骤：
8.如权利要求7所述的方法，其特征在于，所述根据所述曲波系数集中的最大曲波系数和平均曲波系数确定迭代阈值的最大值和最小值，具体包括：
9.如...
【专利技术属性】
技术研发人员：许银坡，倪宇东，张纯，杨丽霞，隆波，柳兴刚，
申请(专利权)人：中国石油天然气集团有限公司，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。