一種黃土塬地震資料井控中長波長靜校正方法
【技術(shù)領(lǐng)域】
[0001] 本發(fā)明屬于石油地球物理勘探領(lǐng)域，具體涉及一種黃土塬地震資料井控中長波長 靜校正方法。
【背景技術(shù)】
[0002] 鄂爾多斯盆地南部黃土塬地區(qū)地表條件溝壑縱橫，地表起伏劇烈，地表巨厚的黃 土產(chǎn)生嚴(yán)重的靜校正問題，雖然在地震資料處理過程中采用一次靜校正（層析靜校正）和 剩余靜校正多次疊代，較好的解決了中波長和短波長靜校正問題，但仍然存在較為嚴(yán)重的 長波長靜校正問題，處理成果存在與地表高程相關(guān)的假構(gòu)造，影響了地震資料的構(gòu)造解釋。 如何消除剩余中長波長靜校正問題是黃土塬地震資料處理的一大難點(diǎn)，目前還沒有有效的 方法解決這個問題。

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

[0003] 本發(fā)明的目的在于解決上述現(xiàn)有技術(shù)中存在的難題，提供一種黃土塬地震資料井 控中長波長靜校正方法，利用地表高程和關(guān)鍵控制井鉆井分層數(shù)據(jù)，通過鉆井分層和地震 剖面同相軸時差，通過對地表高程進(jìn)行平滑，然后通過控制井位置標(biāo)志層時差，擬合中長波 長靜校正量函數(shù)，計算全區(qū)中長波長靜校正量。消除與地表高程相關(guān)的假構(gòu)造。
[0004] 本發(fā)明是通過以下技術(shù)方案實現(xiàn)的：
[0005] -種黃土塬地震資料井控中長波長靜校正方法，包括：
[0006] (1)輸入資料處理范圍的地表高程文件，包括X、Y坐標(biāo)及高程信息；
[0007] (2)計算高程平滑參數(shù)k、1，然后對輸入地表高程文件進(jìn)行平滑；
[0008] 所述步驟（2)是這樣實現(xiàn)的：
[0009] ③計算高程平滑參數(shù)k、1 :
[0010] k = Max_offset/Dcmp_inline,
[0011] I = Max_offset/Dcmp_xline,
[0012] 其中：k為沿測線方向的平滑點(diǎn)數(shù)，
[0013] 1為垂直測線方向平滑點(diǎn)數(shù)，
[0014] Max_offset為地震采集最大排列長度，
[0015] Dcmpjnline為沿測線方向道間距離
[0016] Dcmp_Xline為垂直測線方向道間距離
[0017] ④對輸入的地表高程文件進(jìn)行平滑：
[0018] 利用下式對地表高程文件進(jìn)行平滑：
[0020] 其中：XX(i，j)為平滑后的點(diǎn)（i，j)的地表高程；
[0021] x (m, η)為點(diǎn)（m, η)的地表高程；
[0022] (3)井控中長波長靜校正函數(shù)擬合；
[0023] 所述步驟（3)是這樣實現(xiàn)的：
[0024] 通過選擇擬合函數(shù) s (i, j) = a (XX (i, j) -X。) 3+b (XX (i, j) -X。) +c (2)
[0025] 其中：s (i, j)為點(diǎn)（i, j)中長波長靜校正量，
[0026] XX (i，j)為平滑后點(diǎn)（i，j)的地表高程，
[0027] X。為平滑后的地表高程的算術(shù)平均值，
[0028] a、b、c分別為相關(guān)系數(shù)，
[0029] 通過3 口井以上的數(shù)據(jù)，迭代求解得到系數(shù)a、b、c，然后剔除誤差較多的井，重新 計算求解系數(shù)a、b、c，最后分別求取a、b、c平均值。
[0030] (4)計算中長波長靜校正量。
[0031] 根據(jù)擬合函數(shù)s(i, j) = a(xx(i, j)-xQ)3+b(xx(i, j)_xQ)+c計算得到中長波長靜 校正量。
[0032] 與現(xiàn)有技術(shù)相比，本發(fā)明的有益效果是：本發(fā)明解決了與地表地形相關(guān)的假構(gòu)造 問題，效果非常明顯。本發(fā)明能夠通過少量井的控制有效計算黃土塬地震資料中長波長靜 校正問題。
【附圖說明】
[0033] 圖1是存在剩余中長波長靜校正量的等TO圖。
[0034] 圖2是井控深度圖。
[0035] 圖3是井控剩余中長波長靜校正處理后等TO圖。
[0036] 圖4是井控剩余中長波長靜校正處理前后疊加效果對比。
[0037] 圖5是本發(fā)明方法的步驟框圖。
【具體實施方式】
[0038] 下面結(jié)合附圖對本發(fā)明作進(jìn)一步詳細(xì)描述：
[0039] 剩余中長波長靜校正量與地表1?程具有相關(guān)性，本發(fā)明研制了 一種井控中長波長 靜校正技術(shù)來解決黃土塬地區(qū)剩余中長波長靜校正問題。
[0040] 由于三維地震資料采集，采用多條排列中間激發(fā)觀測系統(tǒng)，受到排列長度的影響， 存在長波長靜校正問題，采用層析靜校正處理技術(shù)，反演近地表模型計算一次靜校正量，反 演近地表模型時一般采用2500m以內(nèi)的排列初至?xí)r間進(jìn)行計算，以提高疊加成像效果為評 判依據(jù)，中長波長靜校正問題難以解決。
[0041] (1)計算高程平滑參數(shù)k、1
[0042] k = Max_offset/Dcmp_inline,
[0043] I = Max_offset/Dcmp_xline,
[0044] 其中：k為沿測線方向的平滑點(diǎn)數(shù)，
[0045] 1為垂直測線方向平滑點(diǎn)數(shù)，
[0046] Max_offset為地震采集最大排列長度，
[0047] Dcmpjnline為沿測線方向道間距離
[0048] Dcmp_xline為垂直測線方向道間距離
[0049] (2)地表高程平滑
[0050] 通過研究，分析了剩余中長波長靜校正問題和地表高程趨勢有相關(guān)性，因此采用 平滑地表高程，得到地表趨勢，計算分析中長波長靜校正量的合理性。
[0052] 其中：XX(i，j)為平滑后點(diǎn)（i，j)的地表高程
[0053] X (m, η)為點(diǎn)（m, η)地表高程
[0054] k為沿測線方向的平滑點(diǎn)數(shù)
[0055] 1為垂直測線方向平滑點(diǎn)數(shù)
[0056] (3)井控中長波長靜校正函數(shù)擬合
[0057] 通過選擇擬合函數(shù) s (i, j) = a (XX (i, j) -X。) 3+b (XX (i, j) -X。) +c (2)
[0058] 其中：s (i, j)為點(diǎn)（i, j)中長波長靜校正量
[0059] XX (i，j)為平滑后點(diǎn)（i，j)的地表高程
[0060] X。為平滑后高程的算術(shù)的平均值
[0061] a、b、c分別為相關(guān)系數(shù)
[0062] 通過3 口井以上的數(shù)據(jù)，就可以迭代求解系數(shù)a、b、c
[0063] 3 口井以上的數(shù)據(jù)，可以分別迭代求解系數(shù)a、b、c后，剔除誤差較多的井，重新計 算求解系數(shù)a、b、c，最后分別求取a、b、c平均值。
[0064] (4)計算中長波長靜校正量
[0065] 根據(jù)擬合函數(shù)s(i, j) = a(xx(i, j)-xQ)3+b(xx(i, j)_xQ)+c計算中長波長靜校正 量。
[0066] (5)技術(shù)實現(xiàn)的步驟
[0067] 如圖5所示，本發(fā)明完整的技術(shù)實現(xiàn)步驟如下：
[0068] ①輸入需要處理的地表高程文件，包括X、Y坐標(biāo)及高程信息。
[0069] ②按照以下公式計算平滑參數(shù)：
[0070] k = Max_offset/Dcmp_inline
[0071] I = Max_offset/Dcmp_xline
[0072] 其中：Max_〇ffSet為地震采集最大排列長度 [0073] k為沿測線方向的平滑點(diǎn)數(shù)
[0074] 1為垂直測線方向平滑點(diǎn)數(shù)
[0075] Dcmpjnline為沿測線方向道間距離
[0076] Dcmp_Xline為垂直測線方向道間距離
[0077] ③按照上文
【發(fā)明內(nèi)容】
（2)中的公式（1)方法，對輸入地表高程文件進(jìn)行平滑；
[0078] ④在地表高程平滑的基礎(chǔ)上，按照上文
【發(fā)明內(nèi)容】
（3)以及公式（2)中的方法計算 出相關(guān)系數(shù)a、b、c。
[0079] ⑤根據(jù)
【發(fā)明內(nèi)容】
（4)計算得到中長波長靜校正量。
[0080] 本發(fā)明是為了解決剩余中長波長靜校正問題，這些問題反映在與地表高程具有相 關(guān)性。本發(fā)明只要輸入原始的地表高程文件，指定對應(yīng)的平滑半徑，程序就能自動輸出高程 平滑文件。再輸入控制井位置及時差文件，就可以自動輸出中長波長靜校正量文件。
[0081] 圖1是剩余中長波長靜校正前的處理成果，得到的目的層等TO圖，圖2是利用鉆 井分層數(shù)據(jù)控制的深度圖，圖3是利用井控約束中長波長靜校正處理后得到的目的層等TO 圖。對比圖1與圖2,構(gòu)造差異較大，存在與地表具有相關(guān)性的假鼻隆構(gòu)造，經(jīng)過處理解釋多 次結(jié)合分析，認(rèn)為還存在剩余中長波長靜校正問題。對比圖2和圖3,等TO圖與鉆井分層數(shù) 據(jù)控制的深度圖趨勢、形態(tài)較為一致，說明該方法較好的解決了黃土塬覆蓋區(qū)嚴(yán)重的靜校 正問題。
[0082] 圖4上部所示是地表高程，中部、下部所示分別是剩余中長波長靜校正處理前后 的疊加剖面。通過對比分析，可以看到井控剩余中長波長靜校正處理后，剖面上對應(yīng)高程較 低位置的表現(xiàn)為"下凹"的假構(gòu)造得到消除。
[0083] 本發(fā)明針對黃土塬地區(qū)存在地較為嚴(yán)重的中長波長靜校正問題，通過輸入地表高 程
[0084] 和控制井位置的時間誤差差，計算工區(qū)內(nèi)的剩余中長波長靜校正量，實際資料表 明本發(fā)明較好的解決了與地表地形相關(guān)的假構(gòu)造問題，效果非常明顯。
[0085] 上述技術(shù)方案只是本發(fā)明的一種實施方式，對于本領(lǐng)域內(nèi)的技術(shù)人員而言，在本 發(fā)明公開了應(yīng)用方法和原理的基礎(chǔ)上，很容易做出各種類型的改進(jìn)或變形，而不僅限于本 發(fā)明上述【具體實施方式】所描述的方法，因此前面描述的方式只是優(yōu)選的，而并不具有限制 性的意義。
【主權(quán)項】
1. 一種黃±堀地震資料井控中長波長靜校正方法，其特征在于：所述方法包括： (1) 輸入資料處理范圍的地表高程文件，包括Χ、γ坐標(biāo)及高程信息； (2) 計算高程平滑參數(shù)k、l，然后對輸入的地表高程文件進(jìn)行平滑； (3) 井控中長波長靜校正函數(shù)擬合； (4) 計算中長波長靜校正量。2. 根據(jù)權(quán)利要求1所述的黃±堀地震資料井控中長波長靜校正方法，其特征在于：所 述步驟（2)是送樣實現(xiàn)的： ① 計算高程平滑參數(shù)k、l; k=Max_offset/Dcmp_inline, 1 =Max_offset/Dcmp_xline, 其中；k為沿測線方向的平滑點(diǎn)數(shù)， 1為垂直測線方向平滑點(diǎn)數(shù)， Max_offset為地震采集最大排列長度， 化mp_inline為沿測線方向道間距離， 化mp_xline為垂直測線方向道間距離； ② 對輸入的地表高程文件進(jìn)行平滑： 利用下式對地表高程文件進(jìn)行平滑：(1) 其中：xx(i，j)為平滑后的點(diǎn)（i，j)的地表高程；x(m,η)為點(diǎn)（m,η)的地表高程。3. 根據(jù)權(quán)利要求2所述的黃±堀地震資料井控中長波長靜校正方法，其特征在于：所 述步驟（3)是送樣實現(xiàn)的： 通過選擇擬合函數(shù)s(i,j) =a(xx(i,j)-x〇)3+b(xx(i,j)-x〇)+c(2) 其中：s(i，j)為點(diǎn)（i，j)中長波長靜校正量，xx(i，j)為平滑后點(diǎn)（i，j)的地表高程， X。為平滑后的地表高程的算術(shù)平均值，a、b、c分別為相關(guān)系數(shù)， 通過3 口井W上的數(shù)據(jù)，迭代求解得到系數(shù)a、b、C，然后剔除誤差較多的井，重新計算 求解系數(shù)a、b、C，最后分別求取a、b、C平均值。4. 根據(jù)權(quán)利要求3所述的黃±堀地震資料井控中長波長靜校正方法，其特征在于：所 述步驟（4)是送樣實現(xiàn)的： 根據(jù)擬合函數(shù)s(i,j) =a(xx(i,j)-X(j)3+b(xx(i，j)-X(j)+c計算得到中長波長靜校正 量。
【專利摘要】本發(fā)明提供了一種黃土塬地震資料井控中長波長靜校正方法，屬于石油地球物理勘探領(lǐng)域。本方法包括：(1)輸入資料處理范圍的地表高程文件，包括X、Y坐標(biāo)及高程信息；(2)計算高程平滑參數(shù)k、l，然后對輸入的地表高程文件進(jìn)行平滑；(3)井控中長波長靜校正函數(shù)擬合；(4)計算中長波長靜校正量。本發(fā)明能夠通過少量井的控制有效計算黃土塬地震資料中長波長靜校正問題。
【IPC分類】G01V1/36
【公開號】CN105319595
【申請?zhí)枴緾N201410374331
【發(fā)明人】居興國, 肖盈, 鄒少峰, 劉燕峰
【申請人】中國石油化工股份有限公司, 中國石油化工股份有限公司石油物探技術(shù)研究院
【公開日】2016年2月10日
【申請日】2014年7月31日