一種精細剩余靜校正方法
【專利摘要】本發(fā)明公開了一種精細剩余靜校正方法，該方法用于處理原始地震數(shù)據(jù)，該原始地震數(shù)據(jù)包括至少一個CMP道集，每個CMP道集包含多個地震道，該方法包括：A）對每個CMP道集中的地震道進行子波變形道預(yù)處理，以消除子波變形嚴重的地震道；B）對每個CMP道集中的地震道進行自適應(yīng)分權(quán)疊加，以生成新的模型道G2(t)；C）將每個CMP道集中的地震道與對應(yīng)的新的CMP疊加模型道G2(t)進行互相關(guān)計算，以得到該CMP道集中的地震道的延遲時；以及D）在對所述至少一個CMP道集都執(zhí)行步驟A）至步驟C）后，采用統(tǒng)計互相關(guān)算法分離炮點和檢波點靜校正量。該方案能極大壓制地震子波變形道對模型道的影響，有效地提高了剩余靜校正方法的計算精度，使地震資料的疊加剖面構(gòu)造成像更加清晰。
【專利說明】一種精細剩余靜校正方法
【技術(shù)領(lǐng)域】
[0001]本發(fā)明涉及地震資料數(shù)字信號處理領(lǐng)域，具體地，涉及一種壓制子波變形道的剩余靜校正方法。
【背景技術(shù)】
[0002]在地震資料的剩余靜校正處理領(lǐng)域中，通常CMP (共中心點)道集中采用統(tǒng)計互相關(guān)算法求解炮點和檢波點靜校正量。該方法認為，在共反射點道集中，來自同一反射點的反射波同相軸應(yīng)為一條直線，通過各道記錄與該道集的疊加道的互相關(guān)即可求得各地震道的剩余靜校正量。
[0003]統(tǒng)計相關(guān)法的基本實現(xiàn)步驟主要有三步:首先生成共反射點道集疊加模型道?’然后利用道集內(nèi)地震道和模型道進行互相關(guān)，通過相關(guān)曲線確定當前地震道靜校正量；最后對整個工區(qū)數(shù)據(jù)進行靜校正量分解，得到每個炮點和每個檢波點的靜校正量。
[0004]在這三步中，模型道的生成占有重要的地位。生成的模型道質(zhì)量直接影響著靜校正量的計算效果。對于建立模型道的最簡單的方法是，把CMP道集內(nèi)的所有各道加起來形成模型道，然后利用每個單道去跟模型道進行互相關(guān)。
[0005]但是，在經(jīng)過多種地震資料的實踐觀察中，發(fā)現(xiàn)地震信號在中深層區(qū)域，由于高頻分量嚴重缺失，導(dǎo)致部分地震道的子波衰減程度很大，發(fā)生嚴重的畸變，由此生成的模型道不再代表準確的構(gòu)造時間，嚴重影響了剩余靜校正的求解精度。
[0006]Chun和Jacewitz (19 77)提出，在求解剩余靜校正量的迭代過程中，還可以在某次迭代時沿CMP增加方向建立模型道，而在下一次迭代時沿CMP減小方向建立模型道，這樣做可以改善低信噪比區(qū)域的模型道，也可以改善測線兩端的模型道。Kirchheimer (1987)提出利用各道與模型道的互相關(guān)系數(shù)作為該道的加權(quán)因子，采用加權(quán)疊加形成模型道，這可確保模型道中質(zhì)量差的道比質(zhì)量高的道的貢獻少。Tjan等(1994)，Baixas (1995)，Larner和Tjan (1995)提出一種較新的建立模型道的方法，這一方法使用傾角時差校正(DMO)和疊前深度偏移形成疊加數(shù)據(jù)，在復(fù)雜構(gòu)造地區(qū)這些方法特別適用，疊加數(shù)據(jù)經(jīng)過正演后可以形成相應(yīng)的參考道或模型道，疊加前的每一道和自己相應(yīng)的參考道進行互相關(guān)可以求取該道的時移，像常規(guī)方法那樣對時移進行分解便可求取剩余靜校正量，該方法利用真實的反射點形成模型道，因而模型道更加標準。
[0007]但是上述提出的優(yōu)化模型道的方法，并沒有對地震子波變形道對模型道影響作出分析，也沒有提出針對性的有效措施，因此不能解決當?shù)卣鹱硬ㄗ冃螄乐貢r，剩余靜校正精度無法提高的問題。

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

[0008]為解決上述問題，本發(fā)明的目的是提供一種能夠壓制地震子波變形道對疊加生成模型道影響的改進的剩余靜校正方法。
[0009]為了實現(xiàn)上述目的，本發(fā)明提供一種精細剩余靜校正方法，該方法用于處理原始地震數(shù)據(jù)，該原始地震數(shù)據(jù)包括至少一個CMP道集，每個CMP道集包含多個地震道，該方法包括:A)對每個CMP道集中的地震道進行子波變形道預(yù)處理，以消除子波變形嚴重的地震道。B)對每個CMP道集中的地震道進行自適應(yīng)分權(quán)疊加，以生成新的模型道G2(t) ;C)將每個CMP道集中的地震道與對應(yīng)的新的CMP疊加模型道G2 (t)進行互相關(guān)計算，以得到該CMP道集中的地震道的延遲時；以及D)在對所述至少一個CMP道集都執(zhí)行步驟A)至步驟
C)后，采用統(tǒng)計互相關(guān)算法分離炮點和檢波點靜校正量。
[0010]通過上述技術(shù)方案，能極大壓制地震子波變形道對模型道的影響，有效地提高了剩余靜校正方法的計算精度，使地震資料的疊加剖面構(gòu)造成像更加清晰。
[0011]本發(fā)明的其他特征和優(yōu)點將在隨后的【具體實施方式】部分予以詳細說明。
【專利附圖】

【附圖說明】
[0012]附圖是用來提供對本發(fā)明的進一步理解，并且構(gòu)成說明書的一部分，與下面的【具體實施方式】一起用于解釋本發(fā)明，但并不構(gòu)成對本發(fā)明的限制。在附圖中:
[0013]圖1是本發(fā)明的實施方式提供的剩余靜校正方法的流程圖；
[0014]圖2a_e是CMP道集中壓制地震子波變形道效果圖；
[0015]圖3a_c是剩余靜校正處理后的疊加效果圖；
[0016]圖4a_c是剩余靜校正處理后的CMP道集圖；
[0017]圖5a和5b是本發(fā)明與常用商業(yè)軟件執(zhí)行的剩余靜校正后疊加效果對比圖；以及
[0018]圖6a_c是本發(fā)明與常用商業(yè)軟件執(zhí)行的剩余靜校正后CMP道集效果對比圖。
【具體實施方式】
[0019]以下結(jié)合附圖對本發(fā)明的【具體實施方式】進行詳細說明。應(yīng)當理解的是，此處所描述的【具體實施方式】僅用于說明和解釋本發(fā)明，并不用于限制本發(fā)明。
[0020]本文中提到的炮點是指激發(fā)震源的位置，檢波點是指接收地震波信號的儀器位置。
[0021]如圖1所示，根據(jù)本發(fā)明的實施方式，提供了一種精細剩余靜校正方法，該方法用于處理原始地震數(shù)據(jù)，該原始地震數(shù)據(jù)包括至少一個CMP道集，每個CMP道集包含多個地震道，例如CMP道集的個數(shù)可以為M個，每個CMP道集中包含的地震道個數(shù)可以為N個，可以記為 gi (t), g2(t)...gn(t), I ≤ η ≤ N。
[0022]該方法可以包括步驟:
[0023]A.對每個CMP道集中的地震道&(0，g2(t)...gn(t)進行子波變形道預(yù)處理，以去除子波變形嚴重的地震道，其目的是不讓這些地震道參與之后的疊加模型道和計算。
[0024]具體來說，步驟A可以包括以下步驟:
[0025]Al:對于每個CMP道集，根據(jù)等式(I)將CMP道集中的地震道例如gi (t)，g2 (t)...gn (t)疊加生成模型道 Gtl (t):
[0026]G0 (O =(O( 1 )
[0027]其中，gj (i)表示CMP道集中第j道地震道，η表示CMP道集中的地震道數(shù)量，(t1； t2)表示地震道的時間范圍；[0028]A2:給定相關(guān)時窗(4。，0，按照等式(2)將地震道81(0，&(0..|(0分別與G0(t)在相關(guān)時窗內(nèi)進行互相關(guān)計算，將每一地震道的最大互相關(guān)系數(shù)記錄成Cl，C2...cn,:
【權(quán)利要求】
1.一種精細剩余靜校正方法，該方法用于處理原始地震數(shù)據(jù)，該原始地震數(shù)據(jù)包括至少一個共中心點(CMP)道集，每個CMP道集包含多個地震道，該方法包括: A)對每個CMP道集中的地震道進行子波變形道預(yù)處理，以消除子波變形嚴重的地震道； B)對每個CMP道集中的地震道進行自適應(yīng)分權(quán)疊加，以生成新的模型道G2(t)； C)將每個CMP道集中的地震道與對應(yīng)的新的CMP疊加模型道G2(t)進行互相關(guān)計算，以得到該CMP道集中的地震道的延遲時；以及 D)在對所述至少一個CMP道集都執(zhí)行步驟A)至步驟C)后，采用統(tǒng)計互相關(guān)算法分離炮點和檢波點靜校正量。
2.根據(jù)權(quán)利要求1所述的方法，其中，所述步驟A)包括: Al)根據(jù)等式(I)將所述CMP道集中的地震道疊加生成模型道GtlU):
3.根據(jù)權(quán)利要求2所述的方法，其中，所述步驟B)包括: BI)將所述CMP道集中沒有在步驟A4)被去除的地震道重新疊加生成新的模型道G1 (t)； B2)在零時刻相關(guān)位置上，根據(jù)等式(3)求取模型道G1 (t)的自相關(guān)系數(shù)M:
4.根據(jù)權(quán)利要求1-3中任意一項所述的方法，其中所述步驟D)包括: Dl)對所述至少一個CMP道集中與每個檢波點對應(yīng)的CMP道集的地震道的時間延遲求和后進行平均，以得到該檢波點對應(yīng)靜校正量；以及 D2)對 所述至少一個CMP道集中與每個炮點對應(yīng)的CMP道集的地震道的時間延遲求和后進行平均，以得到該炮點對應(yīng)靜校正量。
【文檔編號】G01V1/36GK103728665SQ201210387423
【公開日】2014年4月16日 申請日期:2012年10月12日 優(yōu)先權(quán)日:2012年10月12日
【發(fā)明者】吳波, 唐建明, 楊克明, 徐天吉, 甘其剛, 王薦, 馬昭軍, 胡治權(quán), 張虹, 李顯貴, 曹波, 鄭公營, 江瑩瑩, 斯興炎, 丁蔚楠 申請人:中國石油化工股份有限公司, 中國石油化工股份有限公司西南油氣分公司