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

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

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

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