T实现，大大提高了计算效率，具有重要的实际应用价值。
[0256] 本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序 产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实 施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机 可用存储介质（包括但不限于磁盘存储器、CD-ROM、光学存储器等）上实施的计算机程序产 品的形式。
[0257] 本发明是参照根据本发明实施例的方法、设备（系统）、和计算机程序产品的流程 图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流 程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序 指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产 生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实 现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
[0258] 这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特 定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指 令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或 多个方框中指定的功能。
[0259]这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计 算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或 其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一 个方框或多个方框中指定的功能的步骤。
[0260]以上所述仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技 术人员来说，本发明实施例可以有各种更改和变化。凡在本发明的精神和原则之内，所作的 任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
【主权项】
1. 一种不规则地震数据的五维插值处理方法，其特征在于，包括： 将不规则地震数据变换到频率域和空间域； 对于每一个变换到频率域和空间域的地震数据的频率切片，按照快速傅里叶变换算 法，对地震数据的矩阵和向量进行乘积运算，求取五维插值后的频率波数域数据； 将五维插值后的频率波数域数据变换到频率域和空间域，再变换到时间域和空间域。2. 如权利要求1所述的不规则地震数据的五维插值处理方法，其特征在于，所述将不规 则地震数据变换到频率域和空间域之前，包括： 对不规则地震数据进行去噪、静校正和线性动校正处理； 从处理后的不规则地震数据中抽取待五维插值处理的不规则地震数据； 所述将不规则地震数据变换到频率域和空间域，包括:将抽取的待五维插值处理的不 规则地震数据变换到频率域和空间域。3. 如权利要求1所述的不规则地震数据的五维插值处理方法，其特征在于，对于每一个 变换到频率域和空间域的地震数据的频率切片，按照快速傅里叶变换算法，对地震数据的 矩阵和向量进行乘积运算，求取五维插值后的频率波数域数据，包括： 对每一个频率切片i ω = {〇，1，2，…，i ω max}，按照如下公式进行处理：其中，i ?max是根据实际需要选取的频率上限索引，选取原则是使得地震数据d(Xix，tit) 的有效信号能量刚好处于Q~巧。，"的频带内； W0为一个包含罐邮叫个元素的列向量，且W0的每一个元素都等于打.其中，W0 是Wiu中下角标i ω =0时的情况； F = F(xix，kik)，为一个Nx行、n日nm2n3列的矩阵；其中，F(xix，kik)为不等间隔离散逆傅里 叶变换矩阵，ix={〇，l，2，.··，Νχ-1}为矩阵的行索弓|，化={0，1，2, ...，110111112113-1}为矩阵的 列索引.，其中，kik=化o(i〇)，ki(ii)，k2(i2)，k3(i3))ik是四维空间中的 一个点，knQn) = in-Cn，Cn = floo;r((nn-l)/2)，in= {0,1，2,...，ηη-1}，ri = 〇，l，2,3,函数 floorO表示向下取整，kikXix表示四维向量kik与Xix的内积，ik与i〇、ii、i2、i3保持关系式ik = ((io · ni+ii) · η2+?2) · ns+is；为一个包含Νχ个元素的列向量； 函数conjO表示对输入向量的每一个元素求复共辆，并输出一个向量;运算符·*表示 两个向量的对应元素相乘，并得到一个向量；函数sqdO表示对输入向量的每一个元素开 平方，并输出一个向量;函数sumO表示对输入向量的所有元素求和，并输出一个标量;运算 符Η表示求矩阵的共辆转置； 函数y = CG(b，w)表示如下预条件共辆梯度算法，①设定常量，包括:迭代终止相对容差 ?ο1 = 1(Τ4,最大迭代次数maxit = 30，A = FHF;②计算迭代初始变量，包括：'ro = b，P-1 = 1;③计算迭代终止绝对容差化巧=化:/ ^/b'b，.对 于1 = {0，1，2，'''，111日义;[1：-1}，执行如下处理操作： 如果< /"化或者i =maxit，跳出循环，返回y = yi，否则继续循环： zi=w ·本ri; 化二 riHzi; qi 二 Api; Qi 二化/pAii; yi+i 二 yi+QiPi; η+ι 二 n-ctiqi; Pi+l 二 Zi+(化/化-l)PiD4.如权利要求3所述的不规则地震数据的五维插值处理方法，其特征在于，所述函数y = CG(b，w)中，每次迭代需要计算一次矩阵与向量乘积运算q = Ap，其中Α = ρΗρ具有多级分 块To巧litz矩阵结构，采用如下快速算法计算q=Ap: ① iSN日二 η?Π 2η3、Νι 二 112113、化二 Π 3、Ν3二 1, π?η二化 η-1, η 二 0,1,2,3,技如下方式生成长度 为町的向量i;和i;，其中11=0，1，2，3:② 设向量曰二{a(i) I i = ia(io，ii，i2，i3)，iii = 0,1，2，···，mn-l ,11 = 0,1,2,3}是一个四维 数组的向量化表示，其中，13(：[〇414243) = ((：[〇?1]11 + ；[1)?1]12 + 12)?1]13+134为向量索引， 10，ii，i2，i3为四维数组索引，按如下方式对a的元素赋值，初始化变量《=〇，-1， 对于io={〇, 1,2，···，m〇-l}，执行如下循环（1)至（12) (1) 替=咕+邮。） (2) if=i^+i；(4) (3) 对于ii={0,1,2，···，执行如下循环(4)至（12) (4) r：=/；'+i：-〇；) 巧） (6) 对于i2={0,1,2，···，m2-l}，执行如下循环(7)至（12) (7) Γ = (8) /; :=/;;+Γ;(Λ) (9) 对于i3= {0,1，2,…，执行如下循环（10)至（12) (!0)< =這+ ;.'(/') (11) ?;二 ?; +;班5) (12 )3(4(/,,，/|，/:，/;)) = AK'，則； 其中，i(i)表示向量i的第i个元素，a(i)表示向量a的第i个元素，A(ir，i。)表示矩阵A的 第ir行、i。列的元素； ③ 对向量a所表示的四维数组做四维快速傅里叶变换，并返回向量a ; 其中，向量a = {撕非=4咕击?.2石)，V=化U:，…，奶。-]拂。化U，3}是一个四维数组的向 量化表示，4 的，/'i，/':，/.;) = ?(4 . "!1 + /'I) · "!2 + + Z;，i 为向量索引，i0，i 1，i 2，i3为四维数组 索引，;1 =护F&(a)，函数FFT40表示对四维数组做四维快速傅里叶变换，并返回一个四维数 组； ④ 设向量p={p(i) I i = ip(i〇，ii，i2，i3)，in=0,1,2,--,1111-1,11 = 0,1,2,3}是一个四维 数组的向量化表示，其中山（：1〇414243) = ((：[〇?111 + ；[1)?]12 + 12)?113+13，：[为向量索引， 1〇4142 43为四维数组索弓|;设向量口'={口'（1)|1 = 1。'。〇4142 43)411=0，1，2，'''，1]111-1， n = 0,1,2,3}是一个四维数组的向量化表示，其中，ip' (i〇，ii，i2，i3) = ((io · mi+ii) · Π12+ ?2) .m3+i3，i为向量索引，1〇4142 43为四维数组索引；按照如下方式对向量1/的各个元素 赋值：⑤ 计算q'的)，其中，函数FFT4()和IFFT4()分别表示对四维数组做四 维快速傅里叶变换和逆变换，并返回一个四维数组;运算符·*表示两个四维数组的对应元 素相乘，并得到一个四维数组； 向量q' = {q' (i) I i = iq' (i〇，ii，i2，i3)，iii=0，l，2，···，mn-l，n = 〇，l ,2,3}是一个四维数 组的向量化表示，其中，iq' (io，ii，i2，i3) = ((io · mi+ii) · m2+i2) · m3+i3，i为向量索引，io, ii，i2，i3为四维数组索引； ⑥ 设向量q= {q(i) I i = iq(i〇，ii，i2，i3)，iii = 0,1，2，···，ηη-1，n = 〇, 1,2,3}是一个四维 数组的向量化表示，其中山（：1〇414243) = ((：[〇?111 + ：[1)?112 + 12)?113+13，：[为向量索引， io，il，i2，i3为四维数组索引，按照如下方式对向量q的各个元素赋值，q(iq(io，il，i2，i3)) = q' (iq' (io+n〇,ii+ni,i2+n2，i3+n3))，其中0 < iii<nn-l，n = 〇，l，2,3。5.如权利要求3所述的不规则地震数据的五维插值处理方法，其特征在于，还包括:对 71。执行函数永,,。=.乂尸(>-,,,)，其中，1" = {〇,1,2，.''4"。3山 函数失=义巧次)表示对四维数组y的元素位置进行调整，并返回一个四维数组旁，其中： y = {y(i) I i = iy(i〇，ii，i2，i3)，iii=0,1，2，..·，町-1,11=0,1,2,3}是一个四维数组的向 量化表示，iy(io，il，i2,i3) = ((io·nl+il)·n2+i2)·n3+i3，i为向量索引，io,il,i2,i3为四 维数组索引； y 二{>''('')|/ = ';、('"，'.1，'.2，'.,)，'.,,=(化2，...，"。-^二化1义^^ /、'(7,,，/'|，/2，/;：)二（(/,,，"|十/,)，"2+/?)。";十/;，功向量索引山4142 43为四维数组索引； y，/'2, /;)) = y，':2，/;)) ' 其中.'·。= {0,1,2 ,…，ηη-1 }，η = 〇，1，2,3,函数mod(vi, V2)表示VI对V2求Ο~V2-1之间的余数，函数 floor 0表示向下取整。6. 如权利要求5所述的不规则地震数据的五维插值处理方法，其特征在于，还包括:对 索畑执行函数二/厂巧巧".,），其中，?ω={〇，1，2,…，iCO max}; 函数:S = /FFr4巧)表示对四维数组t做四维快速傅里叶逆变换，并返回一个四维数组s， 其中， y = {抑)，?':0, I，2,…- i,/? = 0,1，2,3,[是一个四维数组的向量化表示， 十+/;，i为向量索引，i〇，il，i2，i3为四维数组索引； s = {s(i) I i = is(io，ii，i2，i3)，iii=0,1，2，···，町-1,11 = 0,1,2,3}是一个四维数组的向 量化表示，15(：[〇414243) = ((：[〇?111 + ：[1)?112+12)?113+13，：[为向量索引，：[〇414243为四 维数组索引。7. 如权利要求6所述的不规则地震数据的五维插值处理方法，其特征在于，将插值后的 频率波数域数据变换到频率和空间域，再变换到时间和空间域，完成对不规则地震数据的 五维插值处理，包括： 利用Si。生成插值后的频率域五维数据翁其中，为插值后的空间采样坐标， ix = {0,1,2,... ,η曲 imM-l}，iw = {〇,l,2,... ， 至z,=柄化)，)-1〇。，&〇'2)，乂3巧从、，其中，.皆,;化?）=、"。''巧"·ΗΑν?η= {〇，l，2，... ,ηη-υ，n = 0, 1，2,3, Μ。为插值后各个空间维的采样间隔，ix与i〇、ii、i2、i3保持关系式 .&: = (( /(，?巧+ /, ) ·+ Z'2 ) · "3 + Z'3 ; 战f，巧,。)生成方式如下，用向量S*。代表名战、，巧0 )的一部分 数据，其中，函数