本發(fā)明涉及地震數(shù)據(jù)分析領(lǐng)域，更具體地，涉及地震數(shù)據(jù)等效主頻計算的方法和地震數(shù)據(jù)等效主頻計算的裝置。

背景技術(shù)：

人工地震法是石油天然氣勘探的主要技術(shù)之一。人工地震法石油天然氣勘探數(shù)據(jù)(以下簡稱：地震數(shù)據(jù))包含了大量的地下石油天然氣藏的地質(zhì)信息，其品質(zhì)高低，直接影響油氣勘探的成功率。地震數(shù)據(jù)包含的信息量多少及其真實性，是衡量地震數(shù)據(jù)品質(zhì)的客觀標準。地震數(shù)據(jù)的信噪比高低、地層構(gòu)造形態(tài)的準確程度以及縱向厚度分辨能力的高低是衡量地震數(shù)據(jù)品質(zhì)好壞的三個主要標準。

縱向厚度分辨能力一般稱為分辨力Tr(Resolution)，其倒數(shù)稱為分辨率Pa(Resolving Power)，兩個定義本質(zhì)上是一致的。

隨著油氣田開發(fā)進入中晚期，需要尋找的油氣藏規(guī)模更小，也更加隱蔽，對地震數(shù)據(jù)的分辨率要求就越來越高。分辨率成為地震數(shù)據(jù)品質(zhì)高低的越來越重要的判別標準。如何計算地震數(shù)據(jù)的分辨率，其結(jié)果是否真實，將直接影響地震數(shù)據(jù)品質(zhì)的判斷。

地震子波是構(gòu)成地震數(shù)據(jù)的基本信號，地震子波的分辨率就代表了地震數(shù)據(jù)的分辨率。

地震子波是一個穩(wěn)定并有一定時間延遲(即有一定的時間長度)的信號。地震子波長度越大，相鄰地層反射的相互干涉就會越嚴重，從地震數(shù)據(jù)中得到地層界面反射系數(shù)及地層厚度信息的難度也就越大。只有相鄰的反射子波完全分離，才不會相互干涉，才可能完全準確得到地層界面反射系數(shù)及地層厚度信息。所以，嚴格的分辨率必須以地震子波的時間延續(xù)度來定義。

圖1(a)為一個薄砂體儲層地質(zhì)模型的反射系數(shù)和相應(yīng)的振幅譜；圖1(b)為采用短子波(主頻35Hz)得到的人工合成地震記錄和相應(yīng)的振幅譜；圖1(c)為采用長子波(主頻35Hz)得到的人工合成地震記錄和相應(yīng)的振幅譜，即圖1(b)和圖1(c)的主頻相同但子波波長不同。如圖1(b)所示，短子波合成記錄中，由于子波長度足夠短，所有的砂體界面反射子波沒有相互干涉，砂體的頂、底反射可以明確界定，即使是反射非常弱的“砂體3”頂界面反射也清晰可見。如圖1(c)所示，長子波合成記錄中，由于子波長度太長，“砂體2”和“砂體3”界面反射子波相互干涉，砂體的頂、底反射的位置和極性已無法界定。

圖2(a)為一個薄砂體儲層地質(zhì)模型的反射系數(shù)和相應(yīng)的振幅譜；圖2(b)為采用短子波(主頻30Hz)得到的人工合成地震記錄和相應(yīng)的振幅譜；圖2(c)為采用長子波(主頻15Hz)得到的人工合成地震記錄和相應(yīng)的振幅譜，即圖2(a)和圖2(b)的子波波長和主頻都不相同，圖2(b)為短子波高主頻，圖2(c)為長子波低主頻。如圖2(b)所示，主頻高、長度短的子波合成記錄中，由于子波長度足夠短，所有砂體界面反射都能清晰界定。如圖2(c)所示，主頻低、長度長的子波合成記錄中，由于子波長度太長，無法判定“砂體2”處反射特征是子波相位變化的反應(yīng)，還是多于一個以上界面反射疊加的反應(yīng)；也無法判定“砂體3”處存在的是一個反射界面還是兩個反射界面。

嚴格通過地震子波長度來定義地震數(shù)據(jù)分辨率是最準確的定義方法，但在實際的油氣勘探中難以得到有效應(yīng)用。由于地下反射界面非常復(fù)雜，加上各種噪聲的干擾，不論是在時間域還是在頻率域，地震數(shù)據(jù)中子波的長度都無法較為準確的求得，也就無法用子波長度來精確估算地震數(shù)據(jù)的分辨率。因此，通常很少采用嚴格地震子波長度來定義地震數(shù)據(jù)分辨率。

現(xiàn)有技術(shù)中，常用的地震數(shù)據(jù)分辨率定義方法較多。這些方法的計算結(jié)果不完全一致，但差異并不大。實際生產(chǎn)中，最直觀、最方便使用的分辨率定義方法應(yīng)該有以下兩種，其均具有直觀、方便使用的特點，成為了人們最習慣并經(jīng)常使用的方法。

一、基于時間域波形特征的主頻定義法

正弦波通常采用波長、周期和頻率來描述其特性，但地震子波不同于正弦波，沒有單一不變的波長、周期和頻率。油氣地震勘探中就仿照正弦波的定義方式，采用主波長(也稱為視波長)、主周期(也稱為視周期)和主頻率(也稱為視頻率)來描述地震子波的特征。一般定義子波波形中相鄰波峰(或波谷)的時間間隔為主周期T_b，其倒數(shù)稱為主頻率f_b＝1/T_b。如果已知地震波的速度v，得到主波長為

λ_b＝vT_b＝v/f_b， (1)

主波長、主周期和主頻率實際上是地震子波中起主要作用的頻率成份的反映。

實際生產(chǎn)中，人們習慣先求得地震子波(也就是地震數(shù)據(jù))的主頻率，再計算出其主周期或者主波長，估算出地震數(shù)據(jù)的分辨力。分辨力可表示為厚度或者時間長度。一般認為，相鄰反射子波時間間隔大于或等于半個主周期，就認為是可以分辨的。由于地震數(shù)據(jù)記錄的是反射子波的雙程旅行時間，所以時間分辨力為T_r：

T_r＝(T_b/2)/2＝T_b/4＝1/(4f_b)。 (2)

二、基于頻率域頻譜特征的平均中心頻率定義法

首先定義平均中心頻率f_m為振幅譜曲線所包含面積的1/2分界處的頻率。即設(shè)振幅譜曲線振幅譜為a(f)，包含的面積為A：

$A={\∫}_{f_l}^{f_h}a\left(f\right)df,---\left(3\right)$

其中，f_l、f_h分別為振幅譜的低、高截止頻率。

則平均中心頻率f_m為滿足以下方程的頻率：

$A/2={\∫}_{f_l}^{f_m}a\left(f\right)df.---\left(4\right)$

如本領(lǐng)域技術(shù)人員所公知的，式(3)和(4)中的積分上、下限可根據(jù)地震數(shù)據(jù)的實際情況(比如信噪比等因素)進行調(diào)整，以減少其他因素的干擾，提高估算精度。

與主頻定義法計算分辨力的方式相同，可得到f_m對應(yīng)的主周期T_b和時間分辨力T_r。

T_b＝1/f_m，T_r＝1/(4f_m)。 (5)

平均中心頻率定義法以所有頻率成份的振幅值為基礎(chǔ)，受單一頻率振幅值的影響小，是一個相對穩(wěn)健的計算方法。該方法的計算非常容易，它與主頻定義法的計算結(jié)果非常接近，使用也很方便。

圖3(a)是主頻法的估算方法示意圖，其主周期T_b＝20ms，主頻率f_b＝50Hz，時間分辨力T_r＝10ms；圖3(b)是平均中心頻率法的估算方法示意圖，其平均中心頻率f_m＝50Hz，主周期T_b＝20ms，時間分辨力T_r＝10ms。如圖所示，當子波振幅譜曲線關(guān)于平均中心頻率完全對稱時，兩種方法的結(jié)果完全一致。

但進一步分析可以發(fā)現(xiàn)，子波振幅譜曲線的形狀(如：子波有效頻帶寬度)與子波時間長度直接相關(guān)，可以準確反映子波的分辨力；而主頻和平均中心頻率主要反映了子波波峰或者波谷的寬窄程度，與子波長度沒有直接的關(guān)系，不能真實反映子波的分辨力。因此，在一些情況下，采用諸如主頻定義法和平均中心頻率定義法的現(xiàn)有技術(shù)來評估地震數(shù)據(jù)的分辨率并不可靠。

圖4(a)是主頻為50Hz的子波的波形及其振幅譜；圖4(b)為與圖4(a)時間長度相同但主頻為100Hz的子波的波形及其振幅譜，主頻為100Hz子波的振幅譜相當于將主頻為50Hz子波的振幅譜平移50Hz，所以兩個子波有相同的振幅譜曲線形狀和有效頻帶寬度，也就有相同的時間長度。如本領(lǐng)域技術(shù)人員所公知的，時間長度相同，則兩個子波嚴格意義上的分辨力就應(yīng)該相同的。而無論采用主頻定義法或平均中心頻率定義法，其計算得到的兩個子波的分辨力卻相差一倍，存在很大的誤差，不能反映地震數(shù)據(jù)的真實情況，給地震數(shù)據(jù)品質(zhì)的評估帶來很大的不確定性。

技術(shù)實現(xiàn)要素：

為了克服現(xiàn)有技術(shù)中存在的缺陷，本發(fā)明提出了能夠可靠評估地震數(shù)據(jù)的方法，本發(fā)明還提出了相應(yīng)的裝置。

根據(jù)本發(fā)明的一方面，提出了一種地震數(shù)據(jù)等效主頻計算的方法，該方法包括：得到地震子波的振幅譜a(f)；基于所述振幅譜a(f)確定所述地震子波的等效主頻f_eb：其中，f_l為振幅譜a(f)的低截止頻率，f_h為振幅譜a(f)的高截止頻率，f_p為振幅譜a(f)的極大值對應(yīng)的峰值頻率，f_es是根據(jù)確定的等效起始頻率；基于等效主頻f_eb越高則品質(zhì)越高的原則評估地震數(shù)據(jù)的品質(zhì)。

根據(jù)本發(fā)明的另一方面，在所述振幅譜a(f)關(guān)于平均中心頻率完全對稱或基本對稱的情況下，可以通過下式確定等效主頻f_eb：其中，f_vh是振幅譜a(f)的終止有效頻率，f_vl是振幅譜a(f)的起始有效頻率。

根據(jù)本發(fā)明的另一方面，可以通過下式確定等效主頻f_eb：f_eb＝f_b-f_es，其中，f_b是振幅譜a(f)的主頻率，f_es是根據(jù)確定的等效起始頻率。

根據(jù)本發(fā)明的另一方面，可以通過下式確定等效主頻f_eb：f_eb＝f_b-f_vl，其中，f_b是振幅譜a(f)的主頻率，f_vl是振幅譜a(f)的起始有效頻率。

根據(jù)本發(fā)明的一方面，提出了一種地震數(shù)據(jù)等效主頻計算的裝置，該裝置包括：振幅譜獲取單元，用于得到地震子波的振幅譜a(f)；等效主頻獲取單元，用于基于所述振幅譜a(f)確定所述地震子波的等效主頻f_eb：其中，f_l為振幅譜a(f)的低截止頻率，f_h為振幅譜a(f)的高截止頻率，f_p為振幅譜a(f)的極大值對應(yīng)的峰值頻率，f_es是根據(jù)確定的等效起始頻率；地震數(shù)據(jù)品質(zhì)評估單元，用于基于等效主頻f_eb越高則品質(zhì)越高的原則評估地震數(shù)據(jù)的品質(zhì)。

根據(jù)本發(fā)明的另一方面，在所述振幅譜a(f)關(guān)于平均中心頻率完全對稱或基本對稱的情況下，等效主頻獲取單元可以通過下式確定等效主頻f_eb：其中，f_vh是振幅譜a(f)的終止有效頻率，f_vl是振幅譜a(f)的起始有效頻率。

根據(jù)本發(fā)明的一方面，等效主頻獲取單元可以通過下式確定等效主頻f_eb：f_eb＝f_b-f_es，其中，f_b是振幅譜a(f)的主頻率，f_es是根據(jù)確定的等效起始頻率。

根據(jù)本發(fā)明的另一方面，等效主頻獲取單元可以通過下式確定等效主頻f_eb：f_eb＝f_b-f_vl，其中，f_b是振幅譜a(f)的主頻率，f_vl是振幅譜a(f)的起始有效頻率。

本發(fā)明的各方面通過獲取與地震子波時間長度直接相關(guān)的等效主頻、并基于該等效主頻來評估地震數(shù)據(jù)的品質(zhì)，能夠得到可靠的、確定性的評估結(jié)果。理論模型和實際地震數(shù)據(jù)的計算分析結(jié)果也表明，等效主頻是一個能準確反映地震數(shù)據(jù)真實分辨率的參數(shù)，完全可以表征地震數(shù)據(jù)對地下地質(zhì)體的分辨能力。進一步地，等效主頻的計算和使用簡潔、直觀，是一個方便實用的品質(zhì)衡量參數(shù)。

附圖說明

通過結(jié)合附圖對本發(fā)明示例性實施方式進行更詳細的描述，本發(fā)明的上述以及其它目的、特征和優(yōu)勢將變得更加明顯，其中，在本發(fā)明示例性實施方式中，相同的參考標號通常代表相同部件。

圖1(a)示出了一個薄砂體儲層地質(zhì)模型的反射系數(shù)和相應(yīng)的振幅譜；圖1(b)示出了采用短子波(主頻35Hz)得到的人工合成地震記錄和相應(yīng)的振幅譜；圖1(c)示出了采用長子波(主頻35Hz)得到的人工合成地震記錄和相應(yīng)的振幅譜。

圖2(a)示出了一個薄砂體儲層地質(zhì)模型的反射系數(shù)和相應(yīng)的振幅譜；圖2(b)示出了采用短子波(主頻30Hz)得到的人工合成地震記錄和相應(yīng)的振幅譜；圖2(c)示出了采用長子波(主頻15Hz)得到的人工合成地震記錄和相應(yīng)的振幅譜。

圖3(a)是主頻法的估算方法示意圖；圖3(b)是平均中心頻率法的估算方法示意圖。

圖4(a)是主頻為50Hz的子波的波形及其振幅譜；圖4(b)為與圖4(a)時間長度相同但主頻為100Hz的子波的波形及其振幅譜。

圖5(a)有旁瓣的子波波形及振幅譜曲線；圖5(b)是具有相同時間長度的無旁瓣的子波的波形及振幅譜曲線。

圖6(a)～(c)示出了獲取調(diào)整后的振幅譜曲線a'(f)的示意圖。

圖7(a)為主頻25Hz的帶通理論子波的時域波形圖；圖7(b)為圖7(a)所示子波的振幅譜曲線a(f)；圖7(c)為振幅譜曲線a'(f)；圖7(d)為振幅譜曲線a'(f)對應(yīng)的時域子波波形，等效主頻為13Hz；圖7(e)為圖7(a)所示的子波生成的鍥形體地質(zhì)模型合成地震記錄；圖7(f)為圖7(d)所示的子波生成的鍥形體地質(zhì)模型合成地震記錄。

圖8(a)為主頻25Hz的帶通理論子波的時域波形圖；圖8(b)為圖8(a)所示子波的振幅譜曲線a(f)；圖8(c)為振幅譜曲線a'(f)；圖8(d)為振幅譜曲線a'(f)對應(yīng)的時域子波波形，等效主頻為13Hz；圖8(e)為圖8(a)所示的子波生成的鍥形體地質(zhì)模型合成地震記錄；圖8(f)為圖8(d)所示的子波生成的鍥形體地質(zhì)模型合成地震記錄。

圖9(a)、(b)和(c)分別示出了不同地震子波下的某一河道砂體油氣儲層的地震數(shù)據(jù)。

圖10(a)、(b)和(c)分別示出了不同地震子波下的某一地址異常體地層的地震數(shù)據(jù)。

具體實施方式

下面將參照附圖更詳細地描述本發(fā)明的優(yōu)選實施方式。雖然附圖中顯示了本發(fā)明的優(yōu)選實施方式，然而應(yīng)該理解，可以以各種形式實現(xiàn)本發(fā)明而不應(yīng)被這里闡述的實施方式所限制。相反，提供這些實施方式是為了使本發(fā)明更加透徹和完整，并且能夠?qū)⒈景l(fā)明的范圍完整地傳達給本領(lǐng)域的技術(shù)人員。

此處，先對本發(fā)明涉及的基本原理進行介紹。

發(fā)明人研究發(fā)現(xiàn)，影響子波長度有兩個因素：一是波峰、波谷的寬度；二是波峰、波谷的個數(shù)。一般習慣把主波峰和主波谷之外的波峰和波谷稱為旁瓣，低頻成份的缺失是子波產(chǎn)生旁瓣的原因?，F(xiàn)有技術(shù)中用于評價地震數(shù)據(jù)品質(zhì)的主頻率或平均中心頻率主要反映了子波波峰、波谷的寬度，而不能反映波峰、波谷的個數(shù)。這就是其計算結(jié)果誤差大的根本原因。

如果可以使子波波形只有一個波峰，就可以消除了波峰、波谷個數(shù)這一影響因素，子波時間長度和主頻率(或平均中心頻率)之間就可以建立一種直接的關(guān)系。這樣的主頻率(或平均中心頻率)就可以準確描述地震數(shù)據(jù)的分辨力了。

圖5(a)有旁瓣的子波波形及振幅譜曲線；圖5(b)是具有相同時間長度的無旁瓣的子波的波形及振幅譜曲線。有旁瓣的子波一個特征就是其低頻成份缺失，子波主頻率(或平均中心頻率)與子波時間長度沒有直接關(guān)系，不能反映嚴格的分辨力。無旁瓣子波的特征就是低頻成份不缺失，只有單一主瓣，波峰的寬度就是子波的長度，其主頻率(或平均中心頻率)直接決定了波峰的寬度，也就決定了子波的長度。所以，無旁瓣子波的主頻率(或平均中心頻率)和子波長度直接相關(guān)，完全可以反映出子波嚴格的分辨力。

根據(jù)有旁瓣和無旁瓣子波各自的振幅譜曲線特征，可考慮通過下面的思路來得到地震子波的等效主頻f_eb。

首先可定義等效起始頻率f_es。子波振幅譜曲線經(jīng)過適當?shù)钠交幚碇?，相對平滑，其極大值(通常為1)對應(yīng)的峰值頻率為f_p?？啥x滿足以下條件的頻率f_es為等效起始頻率：如圖6(a)所示，如果頻率f_es的垂直線分割的區(qū)域(1)和區(qū)域(2)面積相等，則f_es可稱為有效起始頻率，即：

${\∫}_{f_l}^{f_{es}}a\left(f\right)df={\∫}_{f_{es}}^{f_p}(1-a(f\left)\right)df.---\left(6\right)$

得到f_es后，可將振幅譜曲線a(f)向左移動f_es，得到圖6(b)所示的新振幅譜曲線。圖6(b)所示曲線中，從-f_es到0Hz的頻率成份是物理不可實現(xiàn)的，沒有物理意義。由f_es的定義可知，區(qū)域(1)和區(qū)域(2)的面積相等。因此，可考慮將區(qū)域(2)面積移動并填充到區(qū)域(1)面積，就可以得到圖6(c)所示的調(diào)整后的振幅譜曲線a'(f)，根據(jù)上述過程可知

與公式(4)類似，可通過下式得到等效平均中心主頻f_em：

$\frac{A}{2}={\∫}_{f_l}^{f_{em}}{a}^{\′}\left(f\right)df.---\left(7\right)$

考慮到難以求得a'(f)在時間域的主頻率，另外考慮到實際生產(chǎn)中人們的使用習慣，可以直接將頻率f_em作為等效主頻f_eb，即可通過下式得到等效主頻f_eb：

$\frac{A}{2}={\∫}_{f_l}^{f_{eb}}{a}^{\′}\left(f\right)df.---\left(8\right)$

由等效主頻f_eb的計算方法可知，等效主頻f_eb與子波時間長度直接相關(guān)。理論上，等效主頻f_eb相同的子波有著相同的子波長度和相同的時間分辨力。等效主頻f_eb符合嚴格的分辨率定義要求，可以真實反映地震數(shù)據(jù)的分辨率。

以上定義了嚴格的等效主頻f_eb。在實際生產(chǎn)中，可對其計算方法進行適當調(diào)整，以犧牲很小的計算精度為代價獲取計算量的顯著減小。

在上述發(fā)明構(gòu)思的基礎(chǔ)上可以得到，當振幅譜曲線關(guān)于平均中心頻率(參見公式(4)的定義)完全對稱或者基本對稱時，可以得到如下的非常簡潔的等效主頻f_eb的近似計算公式：

$f_{eb}=\frac{f_{vh}-f_{vl}}{2},---\left(9\right)$

其中，f_vh是振幅譜a(f)的終止有效頻率，f_vl是振幅譜a(f)的起始有效頻率。有效頻帶寬度是實際應(yīng)用中經(jīng)常使用的概念，如本領(lǐng)域技術(shù)人員所公知的，低、高頻端半幅值點的頻率分別被定義為起始有效頻率f_vl和終止有效頻率f_vh。

發(fā)明人經(jīng)過大量驗證，發(fā)現(xiàn)在一般情況下，等效主頻f_eb可以近似表示為主頻率f_b(參見主頻定義法中定義的主頻率f_b)與等效起始頻率f_es的差，即：

f_eb＝f_b-f_es。 (10)

另外，在實際中，起始有效頻率f_vl和上文中定義的等效起始頻率f_es通常是很接近的，用f_vl代替f_es，其引入的誤差很小，可以滿足大多數(shù)實際資料的精度要求，因此，也可以通過下式得到等效主頻f_eb：

f_eb＝f_b-f_vl。 (11)

實施例1

根據(jù)本發(fā)明的一個實施例，公開了一種地震數(shù)據(jù)等效主頻計算的方法，該方法包括：

步驟101，得到地震子波的振幅譜a(f)；

步驟102，基于所述振幅譜a(f)確定所述地震子波的等效主頻f_eb：

其中，f_l為振幅譜a(f)的低截止頻率，f_h為振幅譜a(f)的高截止頻率，f_p為振幅譜a(f)的極大值對應(yīng)的峰值頻率，f_es是根據(jù)確定的等效起始頻率；

步驟103，基于等效主頻f_eb越高則品質(zhì)越高的原則評估地震數(shù)據(jù)的品質(zhì)。

上述實施例中，通過獲取與地震子波時間長度直接相關(guān)的等效主頻、并基于該等效主頻來評估地震數(shù)據(jù)的品質(zhì)，能夠得到可靠的、確定性的評估結(jié)果。上文中已經(jīng)從理論的角度對本發(fā)明的有益技術(shù)效果進行了詳細說明，以下兩個示例將從實際數(shù)據(jù)實驗的角度進一步證明本發(fā)明的有益技術(shù)效果。

圖7(a)為主頻25Hz的帶通理論子波的時域波形圖；圖7(b)為圖7(a)所示子波的振幅譜曲線a(f)；圖7(c)為振幅譜曲線a'(f)；圖7(d)為振幅譜曲線a'(f)對應(yīng)的時域子波波形，等效主頻為13Hz；圖7(e)為圖7(a)所示的子波生成的鍥形體地質(zhì)模型合成地震記錄；圖7(f)為圖7(d)所示的子波生成的鍥形體地質(zhì)模型合成地震記錄。對比圖7(e)和(f)，可以看出，兩個子波的合成記錄對鍥形體地質(zhì)體的分辨情況非常接近，等效主頻相同的不同子波具有非常接近的分辨率。

圖8(a)為主頻25Hz的帶通理論子波的時域波形圖；圖8(b)為圖8(a)所示子波的振幅譜曲線a(f)；圖8(c)為振幅譜曲線a'(f)；圖8(d)為振幅譜曲線a'(f)對應(yīng)的時域子波波形，等效主頻為13Hz；圖8(e)為圖8(a)所示的子波生成的鍥形體地質(zhì)模型合成地震記錄；圖8(f)為圖8(d)所示的子波生成的鍥形體地質(zhì)模型合成地震記錄。對比圖8(e)和(f)，可以看出，兩個子波的合成記錄對鍥形體地質(zhì)體的分辨情況非常接近，等效主頻相同的不同子波具有非常接近的分辨率。

實施例2

根據(jù)本發(fā)明的一個實施例，公開了一種地震數(shù)據(jù)等效主頻計算的方法，該方法包括：

步驟201，得到地震子波的振幅譜a(f)，其中所述振幅譜a(f)關(guān)于平均中心頻率完全對稱或基本對稱；

步驟202，基于所述振幅譜a(f)確定所述地震子波的等效主頻f_eb：

其中，f_vh是振幅譜a(f)的終止有效頻率，f_vl是振幅譜a(f)的起始有效頻率；

步驟203，基于等效主頻f_eb越高則品質(zhì)越高的原則評估地震數(shù)據(jù)的品質(zhì)。

實施例3

根據(jù)本發(fā)明的一個實施例，公開了一種地震數(shù)據(jù)等效主頻計算的方法，該方法包括：

步驟301，得到地震子波的振幅譜a(f)；

步驟302，基于所述振幅譜a(f)確定所述地震子波的等效主頻f_eb：

f_eb＝f_b-f_es，其中，f_b是振幅譜a(f)的主頻率，f_es是根據(jù)確定的等效起始頻率；

步驟303，基于等效主頻f_eb越高則品質(zhì)越高的原則評估地震數(shù)據(jù)的品質(zhì)。

實施例4

根據(jù)本發(fā)明的一個實施例，公開了一種地震數(shù)據(jù)等效主頻計算的方法，該方法包括：

步驟401，得到地震子波的振幅譜a(f)；

步驟402，基于所述振幅譜a(f)確定所述地震子波的等效主頻f_eb：

f_eb＝f_b-f_vl，其中，f_b是振幅譜a(f)的主頻率，f_vl是振幅譜a(f)的起始有效頻率；

步驟403，基于等效主頻f_eb越高則品質(zhì)越高的原則評估地震數(shù)據(jù)的品質(zhì)。

實施例5

根據(jù)本發(fā)明的一個實施例，公開了一種地震數(shù)據(jù)等效主頻計算的裝置，該裝置包括：

振幅譜獲取單元，用于得到地震子波的振幅譜a(f)；

等效主頻獲取單元，用于基于所述振幅譜a(f)確定所述地震子波的等效主頻f_eb：

其中，f_l為振幅譜a(f)的低截止頻率，f_h為振幅譜a(f)的高截止頻率，f_p為振幅譜a(f)的極大值對應(yīng)的峰值頻率，f_es是根據(jù)確定的等效起始頻率；

地震數(shù)據(jù)品質(zhì)評估單元，用于基于等效主頻f_eb越高則品質(zhì)越高的原則評估地震數(shù)據(jù)的品質(zhì)。

實施例6

根據(jù)本發(fā)明的一個實施例，公開了一種地震數(shù)據(jù)等效主頻計算的裝置，該裝置包括：

振幅譜獲取單元，用于得到地震子波的振幅譜a(f)，其中所述振幅譜a(f)關(guān)于平均中心頻率完全對稱或基本對稱；

等效主頻獲取單元，用于基于所述振幅譜a(f)確定所述地震子波的等效主頻f_eb：

其中，f_vh是振幅譜a(f)的終止有效頻率，f_vl是振幅譜a(f)的起始有效頻率；

地震數(shù)據(jù)品質(zhì)評估單元，用于基于等效主頻f_eb越高則品質(zhì)越高的原則評估地震數(shù)據(jù)的品質(zhì)。

實施例7

根據(jù)本發(fā)明的一個實施例，公開了一種地震數(shù)據(jù)等效主頻計算的裝置，該裝置包括：

振幅譜獲取單元，用于得到地震子波的振幅譜a(f)；

等效主頻獲取單元，用于基于所述振幅譜a(f)確定所述地震子波的等效主頻f_eb：

f_eb＝f_b-f_es，其中，f_b是振幅譜a(f)的主頻率，f_es是根據(jù)確定的等效起始頻率；

地震數(shù)據(jù)品質(zhì)評估單元，用于基于等效主頻f_eb越高則品質(zhì)越高的原則評估地震數(shù)據(jù)的品質(zhì)。

實施例8

根據(jù)本發(fā)明的一個實施例，公開了一種地震數(shù)據(jù)等效主頻計算的裝置，該裝置包括：

振幅譜獲取單元，用于得到地震子波的振幅譜a(f)；

等效主頻獲取單元，用于基于所述振幅譜a(f)確定所述地震子波的等效主頻f_eb：

f_eb＝f_b-f_vl，其中，f_b是振幅譜a(f)的主頻率，f_vl是振幅譜a(f)的起始有效頻率；

地震數(shù)據(jù)品質(zhì)評估單元，用于基于等效主頻f_eb越高則品質(zhì)越高的原則評估地震數(shù)據(jù)的品質(zhì)。

應(yīng)用示例

為便于理解本發(fā)明實施例的方案及其效果，以下給出一個具體應(yīng)用示例。本領(lǐng)域技術(shù)人員應(yīng)理解，該示例僅為了便于理解本發(fā)明，其任何具體細節(jié)并非意在以任何方式限制本發(fā)明。

圖9(a)、(b)和(c)分別示出了不同地震子波下的某一河道砂體油氣儲層的地震數(shù)據(jù)。表1為圖9(a)、(b)和(c)各自的地震數(shù)據(jù)主頻和等效主頻的對比列表。

表1地震數(shù)據(jù)主頻與等效主頻對比表

圖9(a)數(shù)據(jù)的主頻不是最高，但等效主頻最高，其分辨力最高。該數(shù)據(jù)中，河道砂體的形態(tài)和邊界最清楚(橢圓標示處)，箭頭指示的地質(zhì)界面反射同樣清晰，易于追蹤。

圖9(b)數(shù)據(jù)的主頻、等效主頻均較低，其分辨力相對較低。該數(shù)據(jù)中，河道砂體的形態(tài)和邊界已無法識別(橢圓標示處)，但紅色箭頭指示的地質(zhì)界面反射還是清晰的，易于追蹤。

圖9(c)數(shù)據(jù)主頻最高，等效主頻最低，其分辨力最低。該數(shù)據(jù)中，河道砂體的形態(tài)和邊界已很難識別(橢圓標示處)，箭頭指示的地質(zhì)界面反射不清晰，難于追蹤。

圖10(a)、(b)和(c)分別示出了不同地震子波下的某一地址異常體地層的地震數(shù)據(jù)。表2為圖10(a)、(b)和(c)各自的地震數(shù)據(jù)主頻和等效主頻的對比列表。

表2地震數(shù)據(jù)主頻與等效主頻對比表

圖10(a)數(shù)據(jù)的主頻不是最高，但等效主頻最高，其分辨力最高。該數(shù)據(jù)中，地質(zhì)異常體的形態(tài)最清楚(橢圓標示處)，箭頭指示的地質(zhì)界面反射同樣清晰，易于追蹤。

圖10(b)數(shù)據(jù)的主頻、等效主頻均較低，其分辨力相對較低。該數(shù)據(jù)中，地質(zhì)異常體無法識別(橢圓標示處)，箭頭指示的地質(zhì)界面反射可以追蹤，但界面位置不能落實。

圖10(c)數(shù)據(jù)主頻最高，等效主頻最低，其分辨力最低。該數(shù)據(jù)中，地質(zhì)異常體無法識別(橢圓標示處)，箭頭指示的地質(zhì)界面反射難于追蹤。

本發(fā)明可以是系統(tǒng)、方法和/或計算機程序產(chǎn)品。計算機程序產(chǎn)品可以包括計算機可讀存儲介質(zhì)，其上載有用于使處理器實現(xiàn)本發(fā)明的各個方面的計算機可讀程序指令。

計算機可讀存儲介質(zhì)可以是可以保持和存儲由指令執(zhí)行設(shè)備使用的指令的有形設(shè)備。計算機可讀存儲介質(zhì)例如可以是――但不限于――電存儲設(shè)備、磁存儲設(shè)備、光存儲設(shè)備、電磁存儲設(shè)備、半導體存儲設(shè)備或者上述的任意合適的組合。計算機可讀存儲介質(zhì)的更具體的例子(非窮舉的列表)包括：便攜式計算機盤、硬盤、隨機存取存儲器(RAM)、只讀存儲器(ROM)、可擦式可編程只讀存儲器(EPROM或閃存)、靜態(tài)隨機存取存儲器(SRAM)、便攜式壓縮盤只讀存儲器(CD-ROM)、數(shù)字多功能盤(DVD)、記憶棒、軟盤、機械編碼設(shè)備、例如其上存儲有指令的打孔卡或凹槽內(nèi)凸起結(jié)構(gòu)、以及上述的任意合適的組合。這里所使用的計算機可讀存儲介質(zhì)不被解釋為瞬時信號本身，諸如無線電波或者其他自由傳播的電磁波、通過波導或其他傳輸媒介傳播的電磁波(例如，通過光纖電纜的光脈沖)、或者通過電線傳輸?shù)碾娦盘枴?/p>

這里所描述的計算機可讀程序指令可以從計算機可讀存儲介質(zhì)下載到各個計算/處理設(shè)備，或者通過網(wǎng)絡(luò)、例如因特網(wǎng)、局域網(wǎng)、廣域網(wǎng)和/或無線網(wǎng)下載到外部計算機或外部存儲設(shè)備。網(wǎng)絡(luò)可以包括銅傳輸電纜、光纖傳輸、無線傳輸、路由器、防火墻、交換機、網(wǎng)關(guān)計算機和/或邊緣服務(wù)器。每個計算/處理設(shè)備中的網(wǎng)絡(luò)適配卡或者網(wǎng)絡(luò)接口從網(wǎng)絡(luò)接收計算機可讀程序指令，并轉(zhuǎn)發(fā)該計算機可讀程序指令，以供存儲在各個計算/處理設(shè)備中的計算機可讀存儲介質(zhì)中。

用于執(zhí)行本發(fā)明操作的計算機程序指令可以是匯編指令、指令集架構(gòu)(ISA)指令、機器指令、機器相關(guān)指令、微代碼、固件指令、狀態(tài)設(shè)置數(shù)據(jù)、或者以一種或多種編程語言的任意組合編寫的源代碼或目標代碼，所述編程語言包括面向?qū)ο蟮木幊陶Z言—諸如Smalltalk、C++等，以及常規(guī)的過程式編程語言—諸如“C”語言或類似的編程語言。計算機可讀程序指令可以完全地在用戶計算機上執(zhí)行、部分地在用戶計算機上執(zhí)行、作為一個獨立的軟件包執(zhí)行、部分在用戶計算機上部分在遠程計算機上執(zhí)行、或者完全在遠程計算機或服務(wù)器上執(zhí)行。在涉及遠程計算機的情形中，遠程計算機可以通過任意種類的網(wǎng)絡(luò)—包括局域網(wǎng)(LAN)或廣域網(wǎng)(WAN)—連接到用戶計算機，或者，可以連接到外部計算機(例如利用因特網(wǎng)服務(wù)提供商來通過因特網(wǎng)連接)。在一些實施例中，通過利用計算機可讀程序指令的狀態(tài)信息來個性化定制電子電路，例如可編程邏輯電路、現(xiàn)場可編程門陣列(FPGA)或可編程邏輯陣列(PLA)，該電子電路可以執(zhí)行計算機可讀程序指令，從而實現(xiàn)本發(fā)明的各個方面。

這里參照根據(jù)本發(fā)明實施例的方法、裝置(系統(tǒng))和計算機程序產(chǎn)品的流程圖和/或框圖描述了本發(fā)明的各個方面。應(yīng)當理解，流程圖和/或框圖的每個方框以及流程圖和/或框圖中各方框的組合，都可以由計算機可讀程序指令實現(xiàn)。

這些計算機可讀程序指令可以提供給通用計算機、專用計算機或其它可編程數(shù)據(jù)處理裝置的處理器，從而生產(chǎn)出一種機器，使得這些指令在通過計算機或其它可編程數(shù)據(jù)處理裝置的處理器執(zhí)行時，產(chǎn)生了實現(xiàn)流程圖和/或框圖中的一個或多個方框中規(guī)定的功能/動作的裝置。也可以把這些計算機可讀程序指令存儲在計算機可讀存儲介質(zhì)中，這些指令使得計算機、可編程數(shù)據(jù)處理裝置和/或其他設(shè)備以特定方式工作，從而，存儲有指令的計算機可讀介質(zhì)則包括一個制造品，其包括實現(xiàn)流程圖和/或框圖中的一個或多個方框中規(guī)定的功能/動作的各個方面的指令。

也可以把計算機可讀程序指令加載到計算機、其它可編程數(shù)據(jù)處理裝置、或其它設(shè)備上，使得在計算機、其它可編程數(shù)據(jù)處理裝置或其它設(shè)備上執(zhí)行一系列操作步驟，以產(chǎn)生計算機實現(xiàn)的過程，從而使得在計算機、其它可編程數(shù)據(jù)處理裝置、或其它設(shè)備上執(zhí)行的指令實現(xiàn)流程圖和/或框圖中的一個或多個方框中規(guī)定的功能/動作。

附圖中的流程圖和框圖顯示了根據(jù)本發(fā)明的多個實施例的系統(tǒng)、方法和計算機程序產(chǎn)品的可能實現(xiàn)的體系架構(gòu)、功能和操作。在這點上，流程圖或框圖中的每個方框可以代表一個模塊、程序段或指令的一部分，所述模塊、程序段或指令的一部分包含一個或多個用于實現(xiàn)規(guī)定的邏輯功能的可執(zhí)行指令。在有些作為替換的實現(xiàn)中，方框中所標注的功能也可以以不同于附圖中所標注的順序發(fā)生。例如，兩個連續(xù)的方框?qū)嶋H上可以基本并行地執(zhí)行，它們有時也可以按相反的順序執(zhí)行，這依所涉及的功能而定。也要注意的是，框圖和/或流程圖中的每個方框、以及框圖和/或流程圖中的方框的組合，可以用執(zhí)行規(guī)定的功能或動作的專用的基于硬件的系統(tǒng)來實現(xiàn)，或者可以用專用硬件與計算機指令的組合來實現(xiàn)。

以上已經(jīng)描述了本發(fā)明的各實施例，上述說明是示例性的，并非窮盡性的，并且也不限于所披露的各實施例。在不偏離所說明的各實施例的范圍和精神的情況下，對于本技術(shù)領(lǐng)域的普通技術(shù)人員來說許多修改和變更都是顯而易見的。本文中所用術(shù)語的選擇，旨在最好地解釋各實施例的原理、實際應(yīng)用或?qū)κ袌鲋械募夹g(shù)的改進，或者使本技術(shù)領(lǐng)域的其它普通技術(shù)人員能理解本文披露的各實施例。