一種基于蟻群算法的地震數(shù)據(jù)剩余靜校正方法
【專利摘要】本發(fā)明提供了一種基于蟻群算法的地震數(shù)據(jù)剩余靜校正方法。所述方法充分吸收基于全局高效自適應(yīng)尋優(yōu)的蟻群算法的優(yōu)點，采用在炮-檢域進(jìn)行剩余靜校正量非線性反演的思路，不但保證了剩余靜校正估算對輸入低信噪比道集的高的成像分辨率，又突出了剩余靜校正最優(yōu)值收斂的效率，形成基于炮-檢域全局最優(yōu)蟻群算法的地震數(shù)據(jù)剩余靜校正方法，為地震數(shù)據(jù)處理提供更精細(xì)和高效的剩余靜校正計算方法。
【專利說明】一種基于蟻群算法的地震數(shù)據(jù)剩余靜校正方法

【技術(shù)領(lǐng)域】
[0001]本發(fā)明涉及石油天然氣地震勘探【技術(shù)領(lǐng)域】，更具體地講，涉及一種基于炮-檢域全局最優(yōu)蟻群算法的地震數(shù)據(jù)剩余靜校正方法，主要應(yīng)用于石油地震勘探的野外地震資料處理解釋。

【背景技術(shù)】
[0002]靜校正估算是地震數(shù)據(jù)處理中非常重要和及其關(guān)鍵的技術(shù)手段，而復(fù)雜地表條件下的剩余靜校正問題是一個非常復(fù)雜的非線性反演問題，難度大，研究工作少，長期未得到很好解決。隨著我國地震勘探向戈壁、黃土塬、山地等地表復(fù)雜地區(qū)發(fā)展，靜校正問題成了迫切需要解決的任務(wù)。原始地震資料經(jīng)過野外一次靜校正后，剩余的靜校正量可用傳統(tǒng)的方法求取。但當(dāng)?shù)卣鹳Y料的剩余靜校正量大、信噪比低或剩余靜校正量大于地震波波形的1/2周期時，采用傳統(tǒng)的剩余靜校正方法會出現(xiàn)周波跳躍現(xiàn)象，影響地震資料的成像效果，必須尋找新的解決思路。
[0003]從20世紀(jì)70年代開始剩余靜校正方法的研究受到廣大學(xué)者的關(guān)注。先后提出了分解旅行時計算炮點和檢波點靜校正量的廣義線性反演法，統(tǒng)計相關(guān)法，最大能量法等。針對剩余靜校正量本質(zhì)上是一個非線性反演問題，將全局最優(yōu)化方法的模擬退火算法，遺傳算法引入剩余靜校正量的計算中，并更好的結(jié)合模擬退火算法和遺傳算法的優(yōu)勢，提出了綜合并行遺傳算法。這些方法基本上都是在疊前共中心點道集中進(jìn)行計算。常規(guī)的剩余靜校正量的計算基于的假設(shè):①是剩余靜校正量隨機分布原則，在CMP道集(即共中心點道集)中疊加后剩余靜校正量的影響可以消除。②是地表一致性原則，即炮點靜校正量和檢波點靜校正量僅與地表位置有關(guān)。統(tǒng)計相關(guān)法、最大能量法和模擬退火法等剩余靜校正算法都是基于這兩點假設(shè)。但對低信噪比、大剩余靜校正量(例如，剩余靜校正量大于地震波波形的1/2周期)的地震資料，應(yīng)用這種常規(guī)的通過CMP道集獲得模型道的方法無法解決這類資料的靜校正問題。


【發(fā)明內(nèi)容】

[0004]針對現(xiàn)有技術(shù)中存在的不足，本發(fā)明的目的之一在于解決上述現(xiàn)有技術(shù)中存在的一個或多個問題。例如，本發(fā)明的目的之一在于提供一種能夠?qū)崿F(xiàn)對低信噪比、大剩余靜校正量(例如，剩余靜校正量大于地震波波形的1/2周期)的地震資料進(jìn)行剩余靜校正的方法。
[0005]為了實現(xiàn)上述目的，本發(fā)明提供了一種基于蟻群算法的地震數(shù)據(jù)剩余靜校正方法。所述方法包括第一操作或第二操作，其中，所述第一操作包括將單炮地震數(shù)據(jù)經(jīng)過野外靜校正和動校正之后的共中心點道集選排為共炮點道集，并對所述共炮點道集依次進(jìn)行以下處理:
[0006]步驟a:
[0007]Al、根據(jù)預(yù)設(shè)的剩余靜校正量的范圍為共炮點道集的每一道隨機生成一個初始剩余靜校正量，形成共炮點道集剩余靜校正量的初始解集；A2、根據(jù)所述共炮點道集剩余靜校正量的初始解集并利用蟻群算法進(jìn)行搜索，得到共炮點道集剩余靜校正量的最優(yōu)解集并應(yīng)用于對共炮點道集的每一道進(jìn)行剩余靜校正。
[0008]步驟b:
[0009]判斷是否滿足最終終止條件:若滿足最終終止條件，則停止運算；若不滿足最終終止條件，則將利用所述共炮點道集剩余靜校正量的最優(yōu)解集進(jìn)行剩余靜校正處理后的數(shù)據(jù)選排為共檢波點道集，并執(zhí)行步驟C。
[0010]步驟c:
[0011]Cl、根據(jù)預(yù)設(shè)的剩余靜校正量的范圍為共檢波點道集的每一道隨機生成一個初始剩余靜校正量，形成共檢波點道集剩余靜校正量的初始解集；C2、根據(jù)所述共檢波點道集剩余靜校正量的初始解集并利用蟻群算法進(jìn)行搜索，得到共檢波點道集剩余靜校正量的最優(yōu)解集并應(yīng)用于對共檢波點道集的每一道進(jìn)行剩余靜校正。
[0012]步驟d:
[0013]判斷是否滿足最終終止條件:若滿足最終終止條件，則停止運算；若不滿足最終終止條件，則返回執(zhí)行步驟a ;所述第二操作包括將單炮地震數(shù)據(jù)經(jīng)過野外靜校正和動校正后的共中心點道集選排為共檢波點道集，并對所述共檢波點道集依次進(jìn)行以下處理:
[0014]步驟a':采用與所述第一操作中的步驟c相同的操作。
[0015]步驟b':判斷是否滿足最終終止條件:若滿足最終終止條件，則停止運算；若不滿足最終終止條件，則將利用所述共檢波點道集剩余靜校正量的最優(yōu)解集進(jìn)行剩余靜校正處理后的數(shù)據(jù)選排為共炮點道集，并執(zhí)行所述步驟c'。
[0016]步驟C,:采用與所述第一操作中的步驟a相同的操作。
[0017]步驟d':判斷是否滿足最終終止條件:若滿足最終終止條件，則停止運算；若不滿足最終終止條件，則返回執(zhí)行步驟a'。
[0018]根據(jù)本發(fā)明基于蟻群算法的地震數(shù)據(jù)剩余靜校正方法的一個實施例，所述步驟A2包括:
[0019]A21、將m只螞蟻分別放置于共炮點道集剩余靜校正量的第N解集中對應(yīng)道剩余靜校正量所在的節(jié)點上，進(jìn)行炮域內(nèi)的第t次循環(huán)迭代(簡稱迭代)，其中，m等于共炮點道集的總道數(shù)，t和N的初始值均為1，并且當(dāng)N= I時，第N解集為所述共炮點道集剩余靜校正量的初始解集。
[0020]A22、利用第N解集對共炮點道集進(jìn)行第t次剩余靜校正，按下面的等式I計算并比較第t次剩余靜校正之后的共炮點道集的疊加道的振幅能量Et和第t-Ι次剩余靜校正之后的共炮點道集的疊加道的振幅能量Ew，其中，當(dāng)t = I時，Etl為第一次剩余靜校正步驟之前的共炮點道集的疊加道的振幅能量，
[0021]等式1:

【權(quán)利要求】
1.一種基于蟻群算法的地震數(shù)據(jù)剩余靜校正方法，其特征在于，所述方法包括第一操作或第二操作，其中， 所述第一操作包括將單炮地震數(shù)據(jù)經(jīng)過野外靜校正和動校正之后的共中心點道集選排為共炮點道集，并對所述共炮點道集依次進(jìn)行以下處理: 步驟a: Al、根據(jù)預(yù)設(shè)的剩余靜校正量的范圍為共炮點道集的每一道隨機生成一個初始剩余靜校正量，形成共炮點道集剩余靜校正量的初始解集； A2、根據(jù)所述共炮點道集剩余靜校正量的初始解集并利用蟻群算法進(jìn)行搜索，得到共炮點道集剩余靜校正量的最優(yōu)解集并應(yīng)用于對共炮點道集的每一道進(jìn)行剩余靜校正； 步驟b: 判斷是否滿足最終終止條件:若滿足最終終止條件，則停止運算；若不滿足最終終止條件，則將利用所述共炮點道集剩余靜校正量的最優(yōu)解集進(jìn)行剩余靜校正處理后的數(shù)據(jù)選排為共檢波點道集，并執(zhí)行步驟c ; 步驟c: Cl、根據(jù)預(yù)設(shè)的剩余靜校正量的范圍為共檢波點道集的每一道隨機生成一個初始剩余靜校正量，形成共檢波點道集剩余靜校正量的初始解集； C2、根據(jù)所述共檢波點道集剩余靜校正量的初始解集并利用蟻群算法進(jìn)行搜索，得到共檢波點道集剩余靜校正量的最優(yōu)解集并應(yīng)用于對共檢波點道集的每一道進(jìn)行剩余靜校正； 步驟d: 判斷是否滿足最終終止條件:若滿足最終終止條件，則停止運算；若不滿足最終終止條件，則返回執(zhí)行步驟a; 所述第二操作包括將單炮地震數(shù)據(jù)經(jīng)過野外靜校正和動校正后的共中心點道集選排為共檢波點道集，并對所述共檢波點道集依次進(jìn)行以下處理: 步驟a':采用與所述第一操作中的步驟c相同的操作； 步驟b':判斷是否滿足最終終止條件:若滿足最終終止條件，則停止運算；若不滿足最終終止條件，則將利用所述共檢波點道集剩余靜校正量的最優(yōu)解集進(jìn)行剩余靜校正處理后的數(shù)據(jù)選排為共炮點道集，并執(zhí)行所述步驟c,; 步驟c,:采用與所述第一操作中的步驟a相同的操作； 步驟d':判斷是否滿足最終終止條件:若滿足最終終止條件，則停止運算；若不滿足最終終止條件，則返回執(zhí)行步驟a'。
2.根據(jù)權(quán)利要求1所述的基于蟻群算法的地震數(shù)據(jù)剩余靜校正方法，其特征在于，所述步驟A2包括: A21、將m只螞蟻分別放置于共炮點道集剩余靜校正量的第N解集中對應(yīng)道剩余靜校正量所在的節(jié)點上，進(jìn)行炮域內(nèi)的第t次循環(huán)迭代，其中，m等于共炮點道集的總道數(shù)，t和N的初始值均為1，并且當(dāng)N = I時，第N解集為所述共炮點道集剩余靜校正量的初始解集； A22、利用第N解集對共炮點道集進(jìn)行第t次剩余靜校正，按下面的等式I計算并比較第t次剩余靜校正之后的共炮點道集的疊加道的振幅能量Et和第t-Ι次剩余靜校正之后的共炮點道集的疊加道的振幅能量Ew，其中，當(dāng)t = I時，Etl為第一次剩余靜校正步驟之前的共炮點道集的疊加道的振幅能量， 等式1: ε；=ΣΣ<,
k=l /-1 在等式I中，m為共炮點道集的總道數(shù)，T為地震道的總樣點數(shù)，d'為第t次剩余靜校正之后的共炮點道集的第j道的第k個樣點的振幅能量； 當(dāng)Et彡Et^1時，則對第N解集中每一道的剩余靜校正量進(jìn)行隨機擾動并根據(jù)蟻群算法的概率轉(zhuǎn)移方式來選擇共炮點道集中每一道擾動后的剩余靜校正量，形成共炮點道集剩余靜校正量的第N+1解集，其中，所述蟻群算法的概率轉(zhuǎn)移方式為: 首先，根據(jù)下面的等式2計算第t次循環(huán)迭代時位于節(jié)點i的第k只螞蟻轉(zhuǎn)移到下一節(jié)點j的轉(zhuǎn)移概率Ρ? (O，其中，節(jié)點i是第N解集中第k道的剩余靜校正量所對應(yīng)的節(jié)點，節(jié)點j是第N解集中第k道的剩余靜校正量經(jīng)隨機擾動后得到的剩余靜校正量所對應(yīng)的節(jié)點， 等式2:
Γηα Γ?β

(O [η^.(θ],,J----jr./ E allow；E allow,
辦)=ΣΜ)Τ[ιω]Ρ
0，其他 在等式2中，τ Jj (t)表示在第t次迭代時所有螞蟻遺留在i節(jié)點到j(luò)節(jié)點路徑上的信息素量，τ is(t)表示在第t次迭代時所有螞蟻遺留在i節(jié)點到s節(jié)點路徑上的信息素量，i節(jié)點到s節(jié)點路徑代表除i節(jié)點到j(luò)節(jié)點路徑之外其它所有可選路徑，nu(t)為啟發(fā)函數(shù)，是i節(jié)點到j(luò)節(jié)點歐幾里德距離的倒數(shù)，nis(t)是i節(jié)點到s節(jié)點歐幾里德距離的倒數(shù)，α表示轉(zhuǎn)移軌跡的重要性，β表示能見度的相對重要性；然后，比較乂(O與預(yù)設(shè)的接受概率，判斷是否接受螞蟻k從第i節(jié)點轉(zhuǎn)移到第j節(jié)占-^ \\\.當(dāng)> Pucp時，則接受將第N解集中第k道剩余靜校正量經(jīng)擾動后得到的剩余靜校正量作為第N+1解集中第k道的剩余靜校正量； 當(dāng)!)郵!、時，則拒絕將第N解集中第k道剩余靜校正量經(jīng)擾動后得到的剩余靜校正量作為第N+1解集中第k道的剩余靜校正量，并重新對第N解集中第k道的剩余靜校正量進(jìn)行擾動，然后根據(jù)所述蟻群算法的概率轉(zhuǎn)移方式判斷是否接受將第N解集中第k道剩余靜校正量經(jīng)重新擾動后的剩余靜校正量作為第N+1解集中第k道的剩余靜校正量，直到P^)〉Pacp ; 最后，當(dāng)所有螞蟻均完成下一節(jié)點的選擇并形成第N+1解集之后，更新信息素，令t =t+1，N = N+1并返回執(zhí)行步驟A21 ； 當(dāng)Et > Et^1時，則將第N解集作為共炮點道集剩余靜校正量的當(dāng)前最優(yōu)解集； A23、判定是否滿足炮域終止條件:若不滿足炮域終止條件則，則對第N解集中每一道的剩余靜校正量進(jìn)行隨機擾動并將擾動后的各道剩余靜校正量所構(gòu)成的解集直接作為共炮點道集剩余靜校正量的第N+1解集，令t = t+1, N = N+1并返回執(zhí)行步驟A21 ;若滿足炮域終止條件則將第N解集作為共炮點道集剩余靜校正量的最優(yōu)解集。
3.根據(jù)權(quán)利要求2所述的基于蟻群算法的地震數(shù)據(jù)剩余靜校正方法，其特征在于，所述步驟A23的炮域終止條件為: 確定t是否達(dá)到炮域最大迭代次數(shù)和確定第t次剩余靜校正之后的共炮點道集數(shù)據(jù)的疊加道的振幅能量Et和第t-Ι次剩余靜校正步驟之后的共炮點道集數(shù)據(jù)的疊加道的振幅能量Ew的差值是否達(dá)到最小， 當(dāng)Et-Ew沒有達(dá)到最小且t未達(dá)到炮域最大迭代次數(shù)時，判定未滿足炮域終止條件； 當(dāng)Et-Ew達(dá)到最小或者t達(dá)到炮域最大迭代次數(shù)時，判定滿足炮域終止條件； 其中，多次重復(fù)所述步驟A21和A22且無法獲得比Et-Ew更小的差值時判定Et-Ew達(dá)到最小。
4.根據(jù)權(quán)利要求2所述的基于蟻群算法的地震數(shù)據(jù)剩余靜校正方法，其特征在于，所述步驟A22按照下面的等式3進(jìn)行信息素更新， 等式3:
Tij (t+Ι) = (1- ζ ) τ Jj (t) + Δ τ (t)


τη 在等式3中，Δ?) = ΣΔ? —I是信息素衰減系數(shù)，ξ = Pt,并


k=l且P是O?I的常數(shù)，t是迭代次數(shù)。
5.根據(jù)權(quán)利要求2或4中任意一項所述的基于蟻群算法的地震數(shù)據(jù)剩余靜校正方法，其特征在于，所述步驟C2包括: C21、將m'只螞蟻分別放置于共檢波點道集的第N'解集中對應(yīng)道剩余靜校正量所在的節(jié)點上，m'等于共檢波點道集的總道數(shù)，進(jìn)行檢波點域內(nèi)的第t,次循環(huán)迭代，其中，t'和N,的初始值均為1，并且當(dāng)N, = I時，共檢波點道集的第N,解集為所述共檢波點道集剩余靜校正量的初始解集； C22、利用第N,解集對共檢波點道集進(jìn)行第t,次剩余靜校正，利用所述等式I計算并比較第t,次剩余靜校正之后的共檢波點道集的疊加道的振幅能量Et'和第t, -1次剩余靜校正之后的共檢波點道集的疊加道的振幅能量Et_/，其中，Ec/為第一次剩余靜校正步驟之前的共檢波點道集的疊加道的振幅能量， 當(dāng)E/彡Et_/時，則對第N'解集中每一道的剩余靜校正量進(jìn)行隨機擾動并根據(jù)所述蟻群算法的概率轉(zhuǎn)移方式來選擇共檢波點道集中每一道擾動后的剩余靜校正量，形成共檢波點道集剩余靜校正量的第N' +1解集，并更新信息素，令t' =t' +1，N' =N' +1并返回執(zhí)行步驟C21 ； 當(dāng)E/ >Et_/時，則將第N,解集作為共檢波點道集剩余靜校正量的當(dāng)前最優(yōu)解集； C23、判定是否滿足檢波點域終止條件:若不滿足檢波點域終止條件，則對第N'解集中每一道的剩余靜校正量進(jìn)行隨機擾動并將擾動后的各道剩余靜校正量所構(gòu)成的解集直接作為共檢波點道集剩余靜校正量的第N' +1解集，令t' =t' +1，N' =N' +1并返回執(zhí)行步驟C21 ;若滿足檢波點域終止條件則將第N,解集作為共檢波點道集的剩余靜校正量的最優(yōu)解集。
6.根據(jù)權(quán)利要求5所述的基于蟻群算法的地震數(shù)據(jù)剩余靜校正方法，其特征在于，所述步驟C23的檢波點域終止條件為:確定t'是否達(dá)到檢波點域最大迭代次數(shù)和確定第t'次剩余靜校正之后的共檢波點道集數(shù)據(jù)的疊加道的振幅能量Et'和第t, -1次剩余靜校正步驟之后的共檢波點道集數(shù)據(jù)的疊加道的振幅能量Et_/的差值是否達(dá)到最小， 當(dāng)Et' -Et_/沒有達(dá)到最小且t'未達(dá)到檢波點域最大迭代次數(shù)時，判定未滿足檢波點域終止條件； 當(dāng)Et' -Et_/達(dá)到最小或者t'達(dá)到檢波點域最大迭代次數(shù)時，判定滿足檢波點域終止條件； 其中，多次重復(fù)所述步驟C21和C22且無法獲得比E/ -Et_/更小的差值時判定E/ _Et_/達(dá)到最小。
7.根據(jù)權(quán)利要求1所述的基于蟻群算法的地震數(shù)據(jù)剩余靜校正方法，其特征在于，所述步驟b和d'包括:判定是否達(dá)到預(yù)設(shè)的運行時間和確定第t次剩余靜校正之后的共炮點道集數(shù)據(jù)的疊加道的振幅能量Et和第t-Ι次剩余靜校正步驟之后的共炮點道集數(shù)據(jù)的疊加道的振幅能量Ew的差值是否達(dá)到預(yù)定闕值，其中， 當(dāng)Et-Ew大于預(yù)定闕值且未達(dá)到預(yù)設(shè)的運行時間時，判定不滿足最終終止條件； 當(dāng)Et-Ew小于或等于預(yù)定闕值或者達(dá)到預(yù)設(shè)的運行時間時，判定滿足最終終止條件。
8.根據(jù)權(quán)利要求1所述的基于蟻群算法的地震數(shù)據(jù)剩余靜校正方法，其特征在于，所述步驟c^Pb'包括:判定是否達(dá)到預(yù)設(shè)的運行時間和確定第t'次剩余靜校正之后的共檢波點道集數(shù)據(jù)的疊加道的振幅能量Et'和第t, -1次剩余靜校正步驟之后的共檢波點道集數(shù)據(jù)的疊加道的振幅能量Et_/的差值是否達(dá)到預(yù)定闕值，其中， 當(dāng)E/ -Et_/大于預(yù)定闕值且未達(dá)到預(yù)設(shè)的運行時間時，判定不滿足最終終止條件； 當(dāng)E/ -Et_/小于或等于預(yù)定闕值或者達(dá)到預(yù)設(shè)的運行時間時，判定滿足最終終止條件。
【文檔編號】G01V1/36GK104199106SQ201410460617
【公開日】2014年12月10日 申請日期:2014年9月11日 優(yōu)先權(quán)日:2014年9月11日
【發(fā)明者】羅紅明, 王勤耕, 何光明, 李忠, 王棟, 楊海濤, 金德剛, 張華
申請人:中國石油集團川慶鉆探工程有限公司地球物理勘探公司