基于数据组合的叠前深度偏移方法及装置制造方法及图纸

本发明专利技术提供的基于数据组合的叠前深度偏移方法及装置，通过数值方法产生新的震源和与其对应的新的地震记录(将每个地震道的所述频率域数据按照地震道所处的方位角和入射角划分，生成多组合成数据)，并使新的震源和新的记录满足以下条件：1、满足波动方程(对数据进行了傅里叶变换得到频率域数据和波数域数据)；2、与原数据相比，数据量显著降低，从而提高地震数据成像的计算效率；3、生成和使用方便。生成和使用方便。生成和使用方便。

【技术实现步骤摘要】
基于数据组合的叠前深度偏移方法及装置


[0001]本专利技术涉及石油勘探
，更具体的，涉及一种基于数据组合的叠前深度偏移方法及装置。

技术介绍

[0002]针对地震数据偏移成像的质量，基于波动方程的炮域叠前深度偏移优于基于Kirchhoff积分法叠前深度偏移，而炮域叠前深度偏移的计算量是随着炮数的增加而增大的，总的计算时间等于单炮的计算时间乘以总炮数。如果炮数较多则计算效率较低，特别是经过五维数据规则化处理的地震数据，炮数显著增加，此时应用炮域叠前深度偏移变得不现实。

技术实现思路

[0003]为了解决上述不足的至少一个，本专利技术第一方面提供一种基于数据组合的叠前深度偏移方法，包括：
[0004]对包括设定数量地震道的地震道集数据进行傅里叶变换，得到每个地震道的频率域数据和波数域数据；
[0005]将每个地震道的所述频率域数据按照地震道所处的方位角和入射角划分，生成多组合成数据；
[0006]结合每个所述波数场数据对所述合成数据进行叠前深度偏移，生成各合成数据的偏移结果；
[0007]将各合成数据的偏移结果叠加形成偏移成像。
[0008]在某些实施例中，所述对包括设定数量地震道的地震道集数据进行傅里叶变换，包括：
[0009]对所述地震道集数据中的时间域数据进行傅里叶变换，得到频率域数据；
[0010]对所述地震道集数据中的空间域数据进行傅里叶变换，得到波数域数据。
[0011]在某些实施例中，在进行结合每个所述波数场数据对所述合成数据进行叠前深度偏移之前，所述方法还包括：
[0012]将每个所述合成数据中的同一接收点数据合成在一个地震道中，并在每个接收点位置输入一预设地震子波。
[0013]在某些实施例中，所述将每个所述合成数据中的同一接收点数据合成在一个地震道中，包括：
[0014]使用第一公式将每个所述合成数据中的同一接收点数据合成在一个地震道中；其中所述第一公式为：
[0015][0016]其中，r是接收点，s是炮点，d
r
(ω)是共接收点数据w是角频率，i是复数，是相移因子。
[0017]在某些实施例中，所述波数场数据包括：高波数波场数据和低波数波场数据，所述结合每个所述波数场数据对所述合成数据进行叠前深度偏移，生成各合成数据的偏移结果，包括：
[0018]将高波数波场数据在频率波数域对导入的参考速度做波场相位移延拓；
[0019]将高波数波场数据的相位移延拓结果转换到空间域，并用当地速度作差值，生成高波数波场数据的偏移结果；
[0020]将低波数场数据用当地速度在空间域做波场相位移得到低波数场数据的偏移结果；其中，所述高波数场数据偏移结果和低波数场数据偏移结果共同组成对应合成数据的偏移结果。
[0021]本专利技术第二方面实施例提供一种基于数据组合的叠前深度偏移装置，包括：
[0022]傅里叶变换模块，对包括设定数量地震道的地震道集数据进行傅里叶变换，得到每个地震道的频率域数据和波数域数据；
[0023]合成数据生成模块，将每个地震道的所述频率域数据按照地震道所处的方位角和入射角划分，生成多组合成数据；
[0024]偏移结果生成模块，结合每个所述波数场数据对所述合成数据进行叠前深度偏移，生成各合成数据的偏移结果；
[0025]偏移成像模块，将各合成数据的偏移结果叠加形成偏移成像。
[0026]在某些实施例中，所述傅里叶变换模块，包括：
[0027]频率域数据变换单元，对所述地震道集数据中的时间域数据进行傅里叶变换，得到频率域数据；
[0028]波数域数据变换单元，对所述地震道集数据中的空间域数据进行傅里叶变换，得到波数域数据。
[0029]在某些实施例中，所述装置还包括：
[0030]归一地震道模块，将每个所述合成数据中的同一接收点数据合成在一个地震道中，并在每个接收点位置输入一预设地震子波。
[0031]在某些实施例中，所述归一地震道模块使用第一公式将每个所述合成数据中的同一接收点数据合成在一个地震道中；其中所述第一公式为：
[0032][0033]其中，r是接收点，s是炮点，d
r
(ω)是共接收点数据w是角频率，i是复数，是相移因子。
[0034]在某些实施例中，所述波数场数据包括：高波数波场数据和低波数波场数据，所述偏移结果生成模块，包括：
[0035]高波数波场数据相位移单元，将高波数波场数据在频率波数域对导入的参考速度做波场相位移延拓；
[0036]高波数波场数据偏移单元，将高波数波场数据的相位移延拓结果转换到空间域，并用当地速度作差值，生成高波数波场数据的偏移结果；
[0037]低波数波场数据偏移单元，将低波数场数据用当地速度在空间域做波场相位移得到低波数场数据的偏移结果；其中，所述高波数场数据偏移结果和低波数场数据偏移结果共同组成对应合成数据的偏移结果。
[0038]本专利技术第三方面实施例提供一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现如上所述的基于数据组合的叠前深度偏移方法的步骤。
[0039]本专利技术第四方面实施例提供一种计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现如上所述的基于数据组合的叠前深度偏移方法的步骤。
[0040]本专利技术的有益效果如下：
[0041]本专利技术提供的基于数据组合的叠前深度偏移方法及装置，通过数值方法产生新的震源和与其对应的新的地震记录(将每个地震道的所述频率域数据按照地震道所处的方位角和入射角划分，生成多组合成数据)，并使新的震源和新的记录满足以下条件：1、满足波动方程(对数据进行了傅里叶变换得到频率域数据和波数域数据)；2、与原数据相比，数据量显著降低，从而提高地震数据成像的计算效率；3、生成和使用方便。解决了由于炮数较大，或者五维数据规则化处理后炮数进一步增大，导致炮域偏移变得不现实的问题，可以使用本方法，在保证成像质量的前提下，显著降低计算效率。在实际试验中，理论模型和实际数据的测试结果均表明该项专利技术不仅非常有效，而且随着计算效率的显著提高会带来一定的经济效益。
附图说明
[0042]为了更清楚地说明本专利技术实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本专利技术的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
[0043]图1示出本专利技术实施例中基于数据组合的叠前深度偏移方法流程示意图。
[0044]图2示出本专利技术实施例中二维模型入射角为零时的合成数据示意图。
[0045]图3示出本专利技术实施例中二维模型入射角从-45度到45度，间隔为1度，91个合成数据的叠前深度偏移结果示意图。
[0046]图4示出本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于数据组合的叠前深度偏移方法，其特征在于，包括：对包括设定数量地震道的地震道集数据进行傅里叶变换，得到每个地震道的频率域数据和波数域数据；将每个地震道的所述频率域数据按照地震道所处的方位角和入射角划分，生成多组合成数据；结合每个所述波数场数据对所述合成数据进行叠前深度偏移，生成各合成数据的偏移结果；将各合成数据的偏移结果叠加形成偏移成像。2.根据权利要求1所述的叠前深度偏移方法，其特征在于，所述对包括设定数量地震道的地震道集数据进行傅里叶变换，包括：对所述地震道集数据中的时间域数据进行傅里叶变换，得到频率域数据；对所述地震道集数据中的空间域数据进行傅里叶变换，得到波数域数据。3.根据权利要求1所述的叠前深度偏移方法，其特征在于，在进行结合每个所述波数场数据对所述合成数据进行叠前深度偏移之前，所述方法还包括：将每个所述合成数据中的同一接收点数据合成在一个地震道中，并在每个接收点位置输入一预设地震子波。4.根据权利要求2所述的叠前深度偏移方法，其特征在于，所述将每个所述合成数据中的同一接收点数据合成在一个地震道中，包括：使用第一公式将每个所述合成数据中的同一接收点数据合成在一个地震道中；其中所述第一公式为：其中，r是接收点，s是炮点，d
r
(ω)是共接收点数据w是角频率，i是复数，是相移因子。5.根据权利要求2所述的叠前深度偏移方法，其特征在于，所述波数场数据包括：高波数波场数据和低波数波场数据，所述结合每个所述波数场数据对所述合成数据进行叠前深度偏移，生成各合成数据的偏移结果，包括：将高波数波场数据在频率波数域对导入的参考速度做波场相位移延拓；将高波数波场数据的相位移延拓结果转换到空间域，并用当地速度作差值，生成高波数波场数据的偏移结果；将低波数场数据用当地速度在空间域做波场相位移得到低波数场数据的偏移结果；其中，所述高波数场数据偏移结果和低波数场数据偏移结果共同组成对应合成数据的偏移结果。6.一种基于数据组合的叠前深度偏移装置，其特征在于，包括：傅里叶变换模块，对包括设定数量地震道的地震道集数据进行傅里叶变换，得到每个地震道的频率域数据...

【专利技术属性】
技术研发人员：吴晓丰，戴南浔，葛志文，李文卉，王东，武威，
申请(专利权)人：中国石油集团东方地球物理勘探有限责任公司，
类型：发明
国别省市：