擬合初至波獲取短波長靜校正量的方法
【專利摘要】本發(fā)明提供了一種擬合初至波獲取短波長靜校正量的方法。所述方法包括：A)選取測線上的接收點，獲取和所述接收點對應的炮點集合；B)拾取對應于所述接收點和炮點集合的共接收點道集的第一初至時間集；C)對所拾取的共接收點道集的第一初至時間集進行樣條約束擬合以得到樣條約束擬合后的第二初至時間集；D)基于所述共接收點道集的第一初至時間集和樣條約束擬合后的第二初至時間集獲取所述接收點的剩余靜校正量。通過本方法可使地震剖面成像效果好并且計算效率高。
【專利說明】擬合初至波獲取短波長靜校正量的方法
【技術領域】
[0001]本發(fā)明涉及地震勘探資料處理領域。更具體地講，涉及一種獲取復雜地區(qū)二維資料的靜校正量的方法。
【背景技術】
[0002]靜校正問題是制約我國復雜地區(qū)油氣勘探發(fā)展的一個瓶頸問題。靜校正問題分為長波長、短波長校正問題，而短波長靜校正問題解決的好壞直接影響地震剖面成像效果，所以大批的地震勘探工作者對剩余靜校正問題進行了大量研究，形成了初至波、折射波、反射波剩余靜校正方法。
[0003]基于反射波的剩余靜校正方法與速度拾取相結合多次迭代計算，在資料品質較好的地區(qū)可得到較好效果，但在諸如山地、沙漠、厚黃土塬等地區(qū)，地形起伏大，近地表速度、厚度變化劇烈，激發(fā)、接收條件差，資料信噪比低，導致反射波剩余靜校正方法因構建的模型道難以見到有效反射而失效，因此在做反射波剩余靜校正之前有必要做初至波、折射波剩余靜校正。
[0004]折射波剩余靜校正方法在一些地區(qū)取得了一定效果，但是這種方法也需要對時窗內的地震數據做互相關，計算量較大，計算耗時長。另外，折射波盡管能量較強、易識別，但有的地震道因為接收條件差，折射波出現(xiàn)奇異值，影響計算的精準性。而初至波剩余靜校正方法，可人機聯(lián)合拾取初至時，去除不可靠的奇異點，保證計算結果準確、可靠。
[0005]基于初至波的剩余靜校正方法拾取延遲時的前提是將初至時間擬合成一條趨勢線，有學者應用最小二乘法等數學方法將初至時間擬合成一條直線，或分步長擬合成一條折線。前者需要根據偏移距對初至時間進行限制，那么可用的初至時間數量會減少，統(tǒng)計信息的可靠性會降低；后者的不足之處在于計算機無法自動識別直達波、折射波、反射波，無法判斷哪些初至時間來自同一層，導致計算結果不是很精確。

【發(fā)明內容】

[0006]本發(fā)明的目的在于提供一種地震剖面成像效果好并且計算效率高的獲取初至波剩余靜校正量的方法。
[0007]為了實現(xiàn)上述目的，根據本發(fā)明的一方面，提供了一種擬合初至波獲取短波長靜校正量的方法，包括:A)選取測線上的接收點，獲取和所述接收點對應的炮點集合；B)拾取對應于所述接收點和炮點集合的共接收點道集的第一初至時間集；C)對所拾取的共接收點道集的第一初至時間集進行樣條約束擬合以得到樣條約束擬合后的第二初至時間集；D)基于所述共接收點道集的第一初至時間集和樣條約束擬合后的第二初至時間集獲取所述接收點的剩余靜校正量。
[0008]可選地，步驟D)可包括:D1)將所述共接收點道集的第一初至時間集包括的每個初至時間和擬合后的第二初至時間集中相應的初至時間做差，并將每個做差后的差值累加；D2)基于與所述接收點對應的炮點數對Dl)中所累加的差值做平均，獲取所述接收點的剩余靜校正量。
[0009]可選地，所述的方法還可包括:E)檢測所述接收點是否是測線上最后的點，當所述接收點不是測線上最后的點時，改變接收點，返回執(zhí)行步驟A)以對測線上的接收點進行遍歷，并獲取測線上所有接收點的剩余靜校正量。
[0010]根據本發(fā)明的另一方面，提供了一種擬合初至波獲取短波長靜校正量的方法，包括:A)選取測線上的炮點，獲取和所述炮點對應的接收點集合；B)拾取對應于所述炮點和接收點集合的共炮點道集的第一初至時間集；C)對所拾取的共炮點道集的第一初至時間集進行樣條約束擬合以得到樣條約束擬合后的第二初至時間集；D)基于所述共炮點道集的第一初至時間集和樣條約束擬合后的第二初至時間集獲取所述炮點的剩余靜校正量。
[0011]可選地，步驟D)可包括:D1)將所述共炮點道集的第一初至時間集包括的每個初至時間和擬合后的第二初至時間集中相應的初至時間做差，并將每個做差后的差值累加；D2)基于與所述炮點對應的接收點數量對Dl)中所累加的差值做平均，獲取所述炮點的剩余靜校正量。
[0012]可選地，所述的方法還可包括:E)檢測所述炮點是否是測線上最后的點，當所述炮點不是測線上最后的點時，改變炮點，返回執(zhí)行步驟A)以對測線上的炮點進行遍歷，并獲取測線上所有炮點的剩余靜校正量。
[0013]將在接下來的描述中部分闡述本發(fā)明另外的方面和/或優(yōu)點，還有一部分通過描述將是清楚的，或者可以經過本發(fā)明的實施而得知。
【專利附圖】

【附圖說明】
[0014]通過下面結合附圖進行的詳細描述，本發(fā)明的上述和其它目的、特點和優(yōu)點將會變得更加清楚，其中:
[0015]圖1示出了根據本發(fā)明示例性實施例的擬合初至波獲取接收點短波長靜校正量的方法的流程圖；
[0016]圖2示出了根據本發(fā)明示例性實施例的擬合初至波獲取炮點短波長靜校正量的方法的流程圖。
【具體實施方式】
[0017]以下，將參照附圖更充分地描述本發(fā)明的示例性實施例，示例性實施例在附圖中示出。然而，可以以許多不同的形式實施示例性實施例，并且不應被解釋為局限于在此闡述的示例性實施例。相反，提供這些實施例從而本公開將會徹底和完整，并將完全地將示例性實施例的范圍傳達給本領域的技術人員。
[0018]圖1示出了擬合初至波獲取接收點短波長靜校正量的方法的流程圖。
[0019]參照圖1，在步驟S101，選取測線上的接收點，獲取和所述接收點對應的炮點集
八
口 ο
[0020]例如，可通過在地震資料的炮點集合信息中檢索具有所述接收點號的炮點的方式來獲取與所述接收點對應的炮點集合。
[0021]在步驟S102，拾取對應于所述接收點和炮點集合的共接收點道集的第一初至時間集。[0022]在初至時間信息中檢索同時具有所述接收點號和在步驟SlOl中獲取的炮點集合包括的炮點號的初至時間，以獲取第一初至時間集。
[0023]在步驟S103，對所拾取的共接收點道集的第一初至時間集進行樣條約束擬合以得到樣條約束擬合后的第二初至時間集。
[0024]例如，可以使用炮點坐標作為變量建立初至時間和炮點坐標的函數，對初至時間和炮點坐標的函數進行樣條約束擬合。
[0025]在步驟S104，基于所述共接收點道集的第一初至時間集和樣條約束擬合后的第二初至時間集獲取所述接收點的剩余靜校正量。
[0026]例如，可將所述共接收點道集的第一初至時間集包括的每個初至時間和擬合后的第二初至時間集中相應的初至時間做差，并將每個做差后的差值累加；可基于與所述接收點對應的炮點數對所累加的差值做平均，獲取所述接收點的剩余靜校正量。[0027]作為一個示例，第一初至時間集和第二初至時間集包括一一對應的元素，可將第二初至時間集中的每個元素減第一初至時間集中相應的元素得到差值，將每個元素相減的差值累加后除以炮點數獲得平均值，所獲得的平均值是所述接收點的剩余靜校正量。
[0028]此外，所述方法還可包括以下步驟(未示出):可檢測所述接收點是否是測線上最后的點，當所述接收點不是測線上最后的點時，改變接收點，返回執(zhí)行步驟S102至步驟S104以對測線上的接收點進行遍歷，并獲取測線上所有接收點的剩余靜校正量。
[0029]下面將示出擬合初至波獲取短波長靜校正量的方法的詳細實施例以對步驟SlOl至步驟S104進行說明。
[0030]作為一個示例，以一個接收點&為例，首先獲取與該接收點對應的所有炮點S”i=l, 2，…，N。
[0031]拾取接收點&與各炮點所對應的共接收點道集的初至時間Tijt5
[0032]提取所有所對應的共炮點道集，對道集中的初值時間進行擬合處理，得到接收點Rj與各炮點所對應的共接收點道集新的初至時間T’ ij。
[0033]擬合的具體實施方法如下:
[0034]將坐標點及初至時間按如下表示
[0035]f (Xi)=Yi, i=0, I,..., η
[0036]其中，Xi表示炮點坐標，Yi表示對應的初至時間，且令
[0037]x0<x1<...<χη
[0038]要擬合函數f (Xi)=yi得到擬合后的函數g(x)=y須滿足
[0039]min Jg"(x)2dx(I)
[0040]且?guī)в屑s束
[0041]
【權利要求】
1.一種擬合初至波獲取短波長靜校正量的方法，包括: A)選取測線上的接收點，獲取和所述接收點對應的炮點集合； B)拾取對應于所述接收點和炮點集合的共接收點道集的第一初至時間集； C)對所拾取的共接收點道集的第一初至時間集進行樣條約束擬合以得到樣條約束擬合后的第二初至時間集； D)基于所述共接收點道集的第一初至時間集和樣條約束擬合后的第二初至時間集獲取所述接收點的剩余靜校正量。
2.根據權利要求1所述的方法，其中，步驟D)包括: Dl)將所述共接收點道集的第一初至時間集包括的每個初至時間和擬合后的第二初至時間集中相應的初至時間做差，并將每個做差后的差值累加； D2)基于與所述接收點對應的炮點數對Dl)中所累加的差值做平均，獲取所述接收點的剩余靜校正量。
3.根據權利要求1所述的方法，還包括: E)檢測所述接收點是否是測線上最后的點，當所述接收點不是測線上最后的點時，改變接收點，返回執(zhí)行步驟A)以對測線上的接收點進行遍歷，并獲取測線上所有接收點的剩余靜校正量。
4.一種擬合初至波獲取短波長靜校正量的方法，包括: A)選取測線上的炮點，獲取和所述炮點對應的接收點集合； B)拾取對應于所述炮點和接收點集合的共炮點道集的第一初至時間集； C)對所拾取的共炮點道集的第一初至時間集進行樣條約束擬合以得到樣條約束擬合后的第二初至時間集； D)基于所述共炮點道集的第一初至時間集和樣條約束擬合后的第二初至時間集獲取所述炮點的剩余靜校正量。
5.根據權利要求4所述的方法，其中，步驟D)包括: Dl)將所述共炮點道集的第一初至時間集包括的每個初至時間和擬合后的第二初至時間集中相應的初至時間做差，并將每個做差后的差值累加； D2)基于與所述炮點對應的接收點數量對Dl)中所累加的差值做平均，獲取所述炮點的剩余靜校正量。
6.根據權利要求4所述的方法，還包括: E)檢測所述炮點是否是測線上最后的點，當所述炮點不是測線上最后的點時，改變炮點，返回執(zhí)行步驟A)以對測線上的炮點進行遍歷，并獲取測線上所有炮點的剩余靜校正量。
【文檔編號】G01V1/36GK103837895SQ201410086448
【公開日】2014年6月4日 申請日期:2014年3月10日 優(yōu)先權日:2014年3月10日
【發(fā)明者】周強, 陳愛萍, 何光明, 曹琳昱, 楊海濤 申請人:中國石油集團川慶鉆探工程有限公司地球物理勘探公司