本文涉及地震數(shù)據(jù)處理，尤指一種基于三維預(yù)測反褶積的地震數(shù)據(jù)處理方法和裝置。

背景技術(shù)：

1、多次波是地震資料處理，特別是海上資料處理中一個(gè)很難、很棘手的問題。因此，海洋地震數(shù)據(jù)處理迫切需要一項(xiàng)技術(shù)進(jìn)行淺水多次波壓制，提高數(shù)據(jù)分辨率。預(yù)測反褶積不同于拋物線拉東變換、srme等多次波去噪方法,壓制多次波不需要速度等輔助信息，而預(yù)測步長又很容易通過自相關(guān)得到，加上預(yù)測反褶積具有參數(shù)設(shè)置簡單的優(yōu)點(diǎn)，預(yù)測反褶積已成為壓制多次波的一種常用方法，并在消除虛反射、交混回響或其它簡單形式的多次波中得到了較好地應(yīng)用。一些地震資料處理方法中，采用一維預(yù)測反褶積進(jìn)行地震資料處理；因此，實(shí)現(xiàn)三維預(yù)測反褶積的地震數(shù)據(jù)處理方法是亟待解決的問題。

技術(shù)實(shí)現(xiàn)思路

1、本技術(shù)提供了一種基于三維預(yù)測反褶積的地震數(shù)據(jù)處理方法和裝置，該方法能夠考慮三維地震數(shù)據(jù)的多個(gè)方向，預(yù)測出三維預(yù)測反褶積的算子，進(jìn)而進(jìn)行地震數(shù)據(jù)處理，從而獲得在空間上更具有銜接性的地震數(shù)據(jù)；同時(shí)，相較于一維預(yù)測反褶積的地震數(shù)據(jù)處理方法，提高地震數(shù)據(jù)的分辨率，改進(jìn)地震數(shù)據(jù)處理和偏移成像的質(zhì)量。

2、第一方面，本技術(shù)提供了一種基于三維預(yù)測反褶積的地震數(shù)據(jù)處理方法，所述方法包括：

3、基于研究區(qū)目標(biāo)地震數(shù)據(jù)，采用預(yù)先建立的三維預(yù)測反褶積模型確定三維預(yù)測反褶積的算子；

4、根據(jù)所述三維預(yù)測反褶積的算子對所述目標(biāo)地震數(shù)據(jù)進(jìn)行三維預(yù)測反褶積處理，得到處理后的地震數(shù)據(jù)。

5、一種示例性的實(shí)施例中，所述目標(biāo)地震數(shù)據(jù)是三維共檢波點(diǎn)地震道集數(shù)據(jù)。

6、一種示例性的實(shí)施例中，所述預(yù)先建立的三維預(yù)測反褶積模型為：

7、

8、上述模型中，db,m,j表示炮線序號為b，炮點(diǎn)序號為m，時(shí)間方向樣點(diǎn)序號為j、且含多次波的共檢波點(diǎn)地震道集數(shù)據(jù)；pb,m,j表示炮線序號為b，炮點(diǎn)序號為m，時(shí)間方向樣點(diǎn)序號為j、且無多次波的共檢波點(diǎn)地震道集數(shù)據(jù)；fc,n,k為三維預(yù)測反褶積的算子；c表示與b相對應(yīng)的三維預(yù)測算子樣點(diǎn)序號，n表示與m相對應(yīng)的三維預(yù)測算子樣點(diǎn)序號，k表示與j相對應(yīng)的三維預(yù)測算子樣點(diǎn)序號；w為crossline方向半寬點(diǎn)數(shù)，u為inline方向半寬點(diǎn)數(shù)；g是表示三維預(yù)測反褶積計(jì)算的步長；db-c,m-n,j-k表示炮線序號為b-c、炮點(diǎn)序號為m-n、時(shí)間方向樣點(diǎn)序號為j-k、且含多次波的共檢波點(diǎn)地震道集數(shù)據(jù)。

9、一種示例性的實(shí)施例中，所述三維預(yù)測反褶積模型的建立過程，包括：

10、根據(jù)反射系數(shù)和地震子波，確定一維地震數(shù)據(jù)；其中，所述反射系數(shù)為白噪序列，所述地震子波為最小相位的子波；

11、根據(jù)地震數(shù)據(jù)的統(tǒng)計(jì)特性，確定一維預(yù)測反褶積模型；

12、根據(jù)所述一維預(yù)測反褶積模型和三維地震數(shù)據(jù)的三維屬性確定三維預(yù)測反褶積模型。

13、一種示例性的實(shí)施例中，所述根據(jù)反射系數(shù)和地震子波，確定一維地震數(shù)據(jù)，包括：

14、根據(jù)反射系數(shù)和地震子波，利用地震褶積公式計(jì)算一維地震數(shù)據(jù)，其中，所述地震褶積公式為：

15、d(t)＝b(t)*f(t)

16、上述公式中，d(t)為時(shí)間域的地震數(shù)據(jù)，b(t)為時(shí)間域的地震子波，f(t)為時(shí)間域的反射系數(shù)，t為地震數(shù)據(jù)的時(shí)間。

17、一種示例性的實(shí)施例中，所述根據(jù)地震數(shù)據(jù)的統(tǒng)計(jì)特性，確定一維預(yù)測反褶積模型，包括：

18、根據(jù)地震數(shù)據(jù)的統(tǒng)計(jì)特性確定地震數(shù)據(jù)與時(shí)間的關(guān)系；

19、根據(jù)所確定的地震數(shù)據(jù)與時(shí)間的關(guān)系建立一維預(yù)測反褶積模型；

20、其中，所述一維預(yù)測反褶積模型為：

21、

22、上述一維預(yù)測反褶積模型中，dj表示時(shí)間方向樣點(diǎn)序號為j、且含多次波的共檢波點(diǎn)地震道集數(shù)據(jù)；pj表示時(shí)間方向樣點(diǎn)序號為j、且無多次波的共檢波點(diǎn)地震道集數(shù)據(jù)；dj-k表示時(shí)間方向樣點(diǎn)序號為j-k、且含多次波的共檢波點(diǎn)地震道集數(shù)據(jù)；fk表示一維預(yù)測反褶積的算子。

23、第二方面，本發(fā)明實(shí)施例提供一種基于三維預(yù)測反褶積的地震數(shù)據(jù)處理裝置，所述裝置包括：確定三維預(yù)測反褶積的算子模塊和地震數(shù)據(jù)處理模塊；

24、所述確定三維預(yù)測反褶積的算子模塊，用于基于研究區(qū)目標(biāo)地震數(shù)據(jù)，采用預(yù)先建立的三維預(yù)測反褶積模型確定三維預(yù)測反褶積的算子；

25、所述地震數(shù)據(jù)處理模塊，用于根據(jù)所述三維預(yù)測反褶積的算子對所述目標(biāo)地震數(shù)據(jù)進(jìn)行三維預(yù)測反褶積處理，得到處理后的地震數(shù)據(jù)。

26、一種示例性的實(shí)施例中，所述預(yù)先建立的三維預(yù)測反褶積模型為：

27、

28、上述模型中，db,m,j表示炮線序號為b，炮點(diǎn)序號為m，時(shí)間方向樣點(diǎn)序號為j、且含多次波的共檢波點(diǎn)地震道集數(shù)據(jù)；pb,m,j表示炮線序號為b，炮點(diǎn)序號為m，時(shí)間方向樣點(diǎn)序號為j、且無多次波的共檢波點(diǎn)地震道集數(shù)據(jù)；fc,n,k為三維預(yù)測反褶積的算子；c表示與b相對應(yīng)的三維預(yù)測算子樣點(diǎn)序號，n表示與m相對應(yīng)的三維預(yù)測算子樣點(diǎn)序號，k表示與j相對應(yīng)的三維預(yù)測算子樣點(diǎn)序號；w為crossline方向半寬點(diǎn)數(shù)，u為inline方向半寬點(diǎn)數(shù)；g是表示三維預(yù)測反褶積計(jì)算的步長；db-c,m-n,j-k表示炮線序號為b-c、炮點(diǎn)序號為m-n、時(shí)間方向樣點(diǎn)序號為j-k、且含多次波的共檢波點(diǎn)地震道集數(shù)據(jù)。

29、一種示例性的實(shí)施例中，所述三維預(yù)測反褶積模型的建立過程，包括：

30、根據(jù)反射系數(shù)和地震子波，確定一維地震數(shù)據(jù)；其中，所述反射系數(shù)為白噪序列，所述地震子波為最小相位的子波；

31、根據(jù)地震數(shù)據(jù)的統(tǒng)計(jì)特性，確定一維預(yù)測反褶積模型；

32、根據(jù)所述一維預(yù)測反褶積模型和三維地震數(shù)據(jù)的三維屬性確定三維預(yù)測反褶積模型。

33、一種示例性的實(shí)施例中，所述根據(jù)反射系數(shù)和地震子波，確定一維地震數(shù)據(jù)，包括：

34、根據(jù)反射系數(shù)和地震子波，利用地震褶積公式計(jì)算一維地震數(shù)據(jù)，其中，所述地震褶積公式為：

35、d(t)＝b(t)*f(t)

36、上述公式中，d(t)為時(shí)間域的地震數(shù)據(jù)，b(t)為時(shí)間域的地震子波，f(t)為時(shí)間域的反射系數(shù)，t為地震數(shù)據(jù)的時(shí)間。

37、一種示例性的實(shí)施例中，所述根據(jù)地震數(shù)據(jù)的統(tǒng)計(jì)特性，確定一維預(yù)測反褶積模型，包括：

38、根據(jù)地震數(shù)據(jù)的統(tǒng)計(jì)特性確定地震數(shù)據(jù)與時(shí)間的關(guān)系；

39、根據(jù)所確定的地震數(shù)據(jù)與時(shí)間的關(guān)系建立一維預(yù)測反褶積模型；

40、其中，所述一維預(yù)測反褶積模型為：

41、

42、上述一維預(yù)測反褶積模型中，dj表示時(shí)間方向樣點(diǎn)序號為j、且含多次波的共檢波點(diǎn)地震道集數(shù)據(jù)；pj表示時(shí)間方向樣點(diǎn)序號為j、且無多次波的共檢波點(diǎn)地震道集數(shù)據(jù)；dj-k表示時(shí)間方向樣點(diǎn)序號為j-k、且含多次波的共檢波點(diǎn)地震道集數(shù)據(jù)；fk表示一維預(yù)測反褶積的算子。

43、第三方面，本發(fā)明實(shí)施例提供一種具備基于三維預(yù)測反褶積的地震數(shù)據(jù)處理功能的計(jì)算機(jī)程序產(chǎn)品，包括計(jì)算機(jī)程序/指令，其特征在于，該計(jì)算機(jī)程序/指令被處理器執(zhí)行時(shí)實(shí)現(xiàn)上述實(shí)施例中任一項(xiàng)所述的基于三維預(yù)測反褶積的地震數(shù)據(jù)處理方法。

44、與相關(guān)技術(shù)相比，本技術(shù)提供一種基于三維預(yù)測反褶積的地震數(shù)據(jù)處理方法和裝置，所述方法包括：基于研究區(qū)目標(biāo)地震數(shù)據(jù)，采用預(yù)先建立的三維預(yù)測反褶積模型確定三維預(yù)測反褶積的算子；根據(jù)所述三維預(yù)測反褶積的算子對所述目標(biāo)地震數(shù)據(jù)進(jìn)行三維預(yù)測反褶積處理，得到處理后的地震數(shù)據(jù)。本技術(shù)考慮三維地震數(shù)據(jù)的多個(gè)方向，預(yù)測出三維預(yù)測反褶積的算子，進(jìn)而進(jìn)行地震數(shù)據(jù)處理，從而獲得在空間上更具有銜接性的地震數(shù)據(jù)；同時(shí)，相較于一維預(yù)測反褶積的地震數(shù)據(jù)處理方法，提高地震數(shù)據(jù)的分辨率，改進(jìn)地震數(shù)據(jù)處理和偏移成像的質(zhì)量。

45、本技術(shù)的其它特征和優(yōu)點(diǎn)將在隨后的說明書中闡述，并且，部分地從說明書中變得顯而易見，或者通過實(shí)施本技術(shù)而了解。本技術(shù)的其他優(yōu)點(diǎn)可通過在說明書以及附圖中所描述的方案來實(shí)現(xiàn)和獲得。