本發(fā)明涉及石油天然氣地震勘探技術(shù)領(lǐng)域，更具體地講，涉及一種球物理勘探地震資料處理過程中的轉(zhuǎn)換波地震資料處理的剩余靜校正方法。

背景技術(shù)：
隨著多波勘探逐漸得到重視，轉(zhuǎn)換波的地震資料處理也變得越來越頻繁。在實(shí)際工業(yè)生產(chǎn)中轉(zhuǎn)換波處理與傳統(tǒng)的縱波處理有很多不同，其中靜校正量的求取是關(guān)鍵步驟之一。轉(zhuǎn)換波靜校正量是由炮點(diǎn)的縱波靜校正量和檢波點(diǎn)的橫波靜校正量組成，轉(zhuǎn)換波炮點(diǎn)靜校正可以用傳統(tǒng)方法求得，而檢波點(diǎn)橫波靜校正則需要專門的方法進(jìn)行求取。這是因?yàn)榕c縱波相比，橫波受近地表潛水面的影響非常小，因此縱波和橫波靜校正量不存在簡單的比例關(guān)系，使用縱波靜校正量乘以某一系數(shù)來對橫波靜校正量進(jìn)行近似通常無法精細(xì)成像。檢波點(diǎn)橫波靜校正方法可以分為兩大類，一類方法是通過獲得地表的速度結(jié)構(gòu)，將橫向速度不均勻體產(chǎn)生的異常時(shí)差予以校正，這種方法一般又統(tǒng)稱為一次靜校正或長波長靜校正；另一類方法是直接基于地震數(shù)據(jù)得到，利用相鄰地震道上信號的相關(guān)性，通過相關(guān)分析求取各道間的時(shí)差以達(dá)到消除這種時(shí)差的方法，這種方法在地震信號處理領(lǐng)域一般又稱為剩余靜校正。目前對于轉(zhuǎn)換波剩余靜校正，一般的做法是利用基于共接收點(diǎn)道集的相關(guān)統(tǒng)計(jì)剩余靜校正方法解決，即共檢波點(diǎn)疊加相干法。該方法雖然在求取靜校正量的方法上得到了最好的處理效果，但是在轉(zhuǎn)換波資料信噪比低的情況下，需要大量的手工操作和校正，在實(shí)際生產(chǎn)中不利于規(guī)模應(yīng)用，僅適用于構(gòu)造較為簡單，資料信噪比高的情況。在公開號為CN101907727B，公開日為2012年5月30日的中國專利文獻(xiàn)公開了一種面波多分量轉(zhuǎn)換波靜校正方法，屬于勘探石油領(lǐng)域中對地震資料處理方法，該方法基于炮記錄中的瑞雷面波，通過面波的提取、頻散曲線反演求取淺層的地層分層及橫波速度，進(jìn)而精確計(jì)算轉(zhuǎn)換波檢波點(diǎn)靜校正量。Cary和Eaton(CaryP.W.andEatonD.W.S.，1993，Asimplemethodforresolvinglargeconverted-wave(P-SV)statics.Geophysics，58(3)：p429-433)提出利用共接收點(diǎn)疊加道相干法來初步確定較大的靜校正量，該法適用于構(gòu)造平緩地區(qū)，直接求出接收點(diǎn)的橫波靜校正量。郭桂紅等(郭桂紅、王德利、何樵登等，轉(zhuǎn)換波的靜校正，吉林大學(xué)學(xué)報(bào)(地球科學(xué)版)，2003，33(4)：542-544頁)提出的利用共接收點(diǎn)疊加道確定總的接收點(diǎn)靜校正量，然后用圓滑的方法得到對構(gòu)造起伏的估計(jì)，再用總的接收點(diǎn)靜校正量減去構(gòu)造起伏便得到了接收點(diǎn)剩余靜校正量。由于這一方法的原理是基于共接收點(diǎn)道集的模型道相關(guān)統(tǒng)計(jì)分析，它也是轉(zhuǎn)換波靜校正目前常用的方法之一。但是由于橫波反射信號較弱，頻率較低，難以同噪音和縱波反射準(zhǔn)確區(qū)分開，而且在地表復(fù)雜地區(qū)，例如山區(qū)，由于高差的變化使得相鄰檢波點(diǎn)的走時(shí)差異很大，或者因?yàn)榈叵陆橘|(zhì)存在橫向速度不均勻而造成相鄰檢波點(diǎn)的時(shí)差較大難以直接利用相關(guān)分析法進(jìn)行轉(zhuǎn)換橫波的剩余靜校正或者求取出的剩余靜校正量精度難以保證。

技術(shù)實(shí)現(xiàn)要素：
針對現(xiàn)有技術(shù)中存在的不足，本發(fā)明的目的之一在于解決上述現(xiàn)有技術(shù)中存在的一個(gè)或多個(gè)問題。例如，本發(fā)明的目的之一在于提供一種在地表起伏較大或低、降速帶復(fù)雜的地區(qū)，提取轉(zhuǎn)換波的剩余靜校正量以對轉(zhuǎn)換橫波進(jìn)行剩余靜校正的方法。為了實(shí)現(xiàn)上述目的，本發(fā)明提供了一種反射轉(zhuǎn)換橫波共檢波點(diǎn)疊加剖面的剩余靜校正方法。所述方法包括以下步驟：在反射縱波共檢波點(diǎn)疊加剖面上選擇滿足預(yù)定條件的層位線，在所述層位線拾取多個(gè)種子點(diǎn)且所述種子點(diǎn)位于層位線走向發(fā)生改變之處，然后連接所述多個(gè)種子點(diǎn)生成一條連線，將所述連線沿著反射縱波共檢波點(diǎn)疊加剖面表示時(shí)間的縱軸上下移動預(yù)定時(shí)間后構(gòu)成第一時(shí)窗，然后在第一時(shí)窗內(nèi)進(jìn)行反射縱波層位的自動追蹤，得到反射縱波層位線；根據(jù)給定的反射轉(zhuǎn)換橫波層位與反射縱波層位的速度比值將所述反射縱波層位線投影到反射轉(zhuǎn)換橫波共檢波點(diǎn)疊加剖面上得到反射轉(zhuǎn)換橫波的比擬層位線，將所述反射轉(zhuǎn)換橫波的比擬層位線沿著反射轉(zhuǎn)換橫波共檢波點(diǎn)疊加剖面表示時(shí)間的縱軸上下移動預(yù)定時(shí)間后構(gòu)成第二時(shí)窗，然后在所述第二時(shí)窗內(nèi)進(jìn)行反射轉(zhuǎn)換橫波層位的自動追蹤，得到第一反射轉(zhuǎn)換橫波層位線；將所述第一反射轉(zhuǎn)換橫波層位線和所述反射轉(zhuǎn)換橫波的比擬層位線之間的時(shí)間差確定為反射轉(zhuǎn)換橫波的初次剩余靜校正量；利用所述反射轉(zhuǎn)換橫波的初次剩余靜校正量對反射轉(zhuǎn)換橫波共檢波點(diǎn)疊加剖面進(jìn)行初次剩余靜校正，得到初次剩余靜校正后的臨時(shí)剖面；在所述初次剩余靜校正后的臨時(shí)剖面上，將所述第一反射轉(zhuǎn)換橫波層位線沿著所述初次剩余靜校正后的臨時(shí)剖面表示時(shí)間的縱軸上下移動預(yù)定時(shí)間后構(gòu)成第三時(shí)窗，然后在所述第三時(shí)窗內(nèi)進(jìn)行反射轉(zhuǎn)換橫波層位的自動追蹤，得到第二反射轉(zhuǎn)換橫波層位線；將所述第二反射轉(zhuǎn)換橫波層位線和所述第一反射轉(zhuǎn)換橫波層位線之間的時(shí)間差確定為反射轉(zhuǎn)換橫波的第二次剩余靜校正量；將初次剩余靜校正量和第二次剩余靜校正量進(jìn)行疊加，得到反射轉(zhuǎn)換橫波共檢波點(diǎn)疊加剖面的最終剩余靜校正量；利用所述反射轉(zhuǎn)換橫波共檢波點(diǎn)疊加剖面的最終剩余靜校正量對反射轉(zhuǎn)換橫波共檢波點(diǎn)疊加剖面進(jìn)行第二次剩余靜校正。根據(jù)本發(fā)明反射轉(zhuǎn)換橫波共檢波點(diǎn)疊加剖面的剩余靜校正方法的一個(gè)實(shí)施例，所述反射縱波層位的自動追蹤通過以下步驟實(shí)現(xiàn)：拾取所述第一時(shí)窗內(nèi)每一道地震數(shù)據(jù)的最大振幅能量所對應(yīng)的時(shí)間并投影到反射縱波共檢波點(diǎn)疊加剖面上，得到所述反射縱波層位線；或者，將所述第一時(shí)窗內(nèi)的第i道地震數(shù)據(jù)與第i+1道地震數(shù)據(jù)做互相關(guān)運(yùn)算，i從1取到N-1且i為正整數(shù)，N為地震道總數(shù)，然后拾取互相關(guān)運(yùn)算后的每一道地震數(shù)據(jù)的最大相關(guān)量所對應(yīng)的時(shí)間并投影到反射縱波共檢波點(diǎn)疊加剖面上，得到所述反射縱波層位線；或者，將所述第一時(shí)窗內(nèi)的第i-n道至第i+n道地震數(shù)據(jù)進(jìn)行疊加得到種子道，然后將所述種子道與第i道地震數(shù)據(jù)做互相關(guān)運(yùn)算，i從n取到N-n且i為正整數(shù)，1≤n≤N且為正整數(shù)，N為地震道總數(shù)，然后拾取互相關(guān)運(yùn)算后的每一道地震數(shù)據(jù)的最大相關(guān)量所對應(yīng)的時(shí)間并投影到反射縱波共檢波點(diǎn)疊加剖面上，得到所述反射縱波層位線。根據(jù)本發(fā)明反射轉(zhuǎn)換橫波共檢波點(diǎn)疊加剖面的剩余靜校正方法的一個(gè)實(shí)施例，所述第二時(shí)窗內(nèi)進(jìn)行的反射轉(zhuǎn)換橫波層位的自動追蹤步驟包括：拾取所述第二時(shí)窗內(nèi)每一道地震數(shù)據(jù)的最大振幅能量所對應(yīng)的時(shí)間并投影到反射轉(zhuǎn)換橫波共檢波點(diǎn)疊加剖面上，得到所述第一反射轉(zhuǎn)換橫波層位線；或者，將所述第二時(shí)窗內(nèi)的第i道地震數(shù)據(jù)與第i+1道地震數(shù)據(jù)做互相關(guān)運(yùn)算，i從1取到N-1且i為正整數(shù)，N為地震道總數(shù)，然后拾取互相關(guān)運(yùn)算后的每一道地震數(shù)據(jù)的最大相關(guān)量所對應(yīng)的時(shí)間并投影到反射轉(zhuǎn)換橫波共檢波點(diǎn)疊加剖面上，得到所述第一反射轉(zhuǎn)換橫波層位線；或者，將所述第二時(shí)窗內(nèi)的第i-n道至第i+n道地震數(shù)據(jù)進(jìn)行疊加得到種子道，然后將所述種子道與第i道地震數(shù)據(jù)做互相關(guān)運(yùn)算，i從n取到N-n且i為正整數(shù)，1≤n≤N且為正整數(shù)，N為地震道總數(shù)，拾取互相關(guān)運(yùn)算后的每一道地震數(shù)據(jù)的最大相關(guān)量所對應(yīng)的時(shí)間并投影到反射轉(zhuǎn)換橫波共檢波點(diǎn)疊加剖面上，得到所述第一反射轉(zhuǎn)換橫波層位線。根據(jù)本發(fā)明反射轉(zhuǎn)換橫波共檢波點(diǎn)疊加剖面的剩余靜校正方法的一個(gè)實(shí)施例，所述第三時(shí)窗內(nèi)進(jìn)行的反射轉(zhuǎn)換橫波層位的自動追蹤步驟包括：拾取所述第三時(shí)窗內(nèi)每一道地震數(shù)據(jù)的最大振幅能量所對應(yīng)的時(shí)間并投影到所述初次剩余靜校正后的臨時(shí)剖面上，得到所述第二反射轉(zhuǎn)換橫波層位線；或者，將所述第三時(shí)窗內(nèi)的第i道地震數(shù)據(jù)與第i+1道地震數(shù)據(jù)做互相關(guān)運(yùn)算，i從1取到N-1且i為正整數(shù)，N為地震道總數(shù)，然后拾取互相關(guān)運(yùn)算后的每一道地震數(shù)據(jù)的最大相關(guān)量所對應(yīng)的時(shí)間并投影到所述初次剩余靜校正后的臨時(shí)剖面上，得到所述第二反射轉(zhuǎn)換橫波層位線；或者，將所述第三時(shí)窗內(nèi)的第i-n道至第i+n道地震數(shù)據(jù)進(jìn)行疊加得到種子道，然后將所述種子道與第i道地震數(shù)據(jù)做互相關(guān)運(yùn)算，i從n取到N-n且i為正整數(shù)，1≤n≤N且為正整數(shù)，N為地震道總數(shù)，然后拾取互相關(guān)運(yùn)算后的每一道地震數(shù)據(jù)的最大相關(guān)量所對應(yīng)的時(shí)間并投影到所述初次剩余靜校正后的臨時(shí)剖面上，得到所述第二反射轉(zhuǎn)換橫波層位線。根據(jù)本發(fā)明反射轉(zhuǎn)換橫波共檢波點(diǎn)疊加剖面的剩余靜校正方法的一個(gè)實(shí)施例，所述第三時(shí)窗內(nèi)進(jìn)行的反射轉(zhuǎn)換橫波層位的自動追蹤通過以下方式實(shí)現(xiàn)：將所述第三時(shí)窗內(nèi)的第i道地震數(shù)據(jù)與第i+1道地震數(shù)據(jù)做互相關(guān)運(yùn)算，i從1取到N-1且i為正整數(shù)，N為地震道總數(shù)，然后拾取互相關(guān)運(yùn)算后的每一道地震數(shù)據(jù)的最大相關(guān)量所對應(yīng)的時(shí)間并投影到所述初次剩余靜校正后的臨時(shí)剖面上，得到所述第二反射轉(zhuǎn)換橫波層位線；或者，將所述第三時(shí)窗內(nèi)的第i-n道至第i+n道地震數(shù)據(jù)進(jìn)行疊加得到種子道，然后將所述種子道與第i道地震數(shù)據(jù)做互相關(guān)運(yùn)算，i從n取到N-n且i為正整數(shù)，1≤n≤N且為正整數(shù)，N為地震道總數(shù)，然后拾取互相關(guān)運(yùn)算后的每一道地震數(shù)據(jù)的最大相關(guān)量所對應(yīng)的時(shí)間并投影到所述初次剩余靜校正后的臨時(shí)剖面上，得到所述第二反射轉(zhuǎn)換橫波層位線。根據(jù)本發(fā)明反射轉(zhuǎn)換橫波共檢波點(diǎn)疊加剖面的剩余靜校正方法的一個(gè)實(shí)施例，所述方法還包括對所述種子道的地震數(shù)據(jù)進(jìn)行求平均處理。根據(jù)本發(fā)明反射轉(zhuǎn)換橫波共檢波點(diǎn)疊加剖面的剩余靜校正方法的一個(gè)實(shí)施例，所述預(yù)定條件表示在所述反射縱波共檢波點(diǎn)疊加剖面上選擇的層位線所對應(yīng)的同相軸連續(xù)清晰。根據(jù)本發(fā)明反射轉(zhuǎn)換橫波共檢波點(diǎn)疊加剖面的剩余靜校正方法的一個(gè)實(shí)施例，所述方法包括在所述反射縱波共檢波點(diǎn)疊加剖面表示時(shí)間的縱軸上分層和/或在表示地震道的橫軸上分段選擇滿足預(yù)定條件的層位線，并分別求出所選各層位線對應(yīng)的反射轉(zhuǎn)換橫波共檢波點(diǎn)的剩余靜校正量。根據(jù)本發(fā)明反射轉(zhuǎn)換橫波共檢波點(diǎn)疊加剖面的剩余靜校正方法的一個(gè)實(shí)施例，所述速度比值在所述反射縱波共檢波點(diǎn)疊加剖面表示地震道的橫軸上分段取不同值，或者在所述反射縱波共檢波點(diǎn)疊加剖面所有地震道上取同一個(gè)值。根據(jù)本發(fā)明反射轉(zhuǎn)換橫波共檢波點(diǎn)疊加剖面的剩余靜校正方法的一個(gè)實(shí)施例，所述方法還包括對所述反射縱波共檢波點(diǎn)疊加剖面和所述反射轉(zhuǎn)換橫波共檢波點(diǎn)疊加剖面進(jìn)行去噪和振幅預(yù)處理，以突出同相軸。與現(xiàn)有技術(shù)相比，本發(fā)明的有益效果包括：充分提取了PS波共接收點(diǎn)疊加剖面上的短波長剩余靜校正量，使PS波同相軸更連續(xù)，更收斂，提高了PS波疊加剖面的質(zhì)量。附圖說明通過下面結(jié)合附圖進(jìn)行的描述，本發(fā)明的上述和其他目的和特點(diǎn)將會變得更加清楚，其中：圖1示出了本發(fā)明示例性實(shí)施例的反射轉(zhuǎn)換橫波共檢波點(diǎn)疊加剖面的剩余靜校正方法的在反射縱波共檢波點(diǎn)疊加剖面選取的層位線上拾取種子點(diǎn)步驟的示意圖。圖2示出了本發(fā)明示例性實(shí)施例的反射轉(zhuǎn)換橫波共檢波點(diǎn)疊加剖面的剩余靜校正方法的構(gòu)成第一時(shí)窗步驟的示意圖。附圖標(biāo)記說明：1-層位線，2、3、4、5、6-種子點(diǎn)，7-第一時(shí)窗的上邊界線、8-第一時(shí)窗的下邊界線以及9-種子點(diǎn)連線。具體實(shí)施方式在下文中，將結(jié)合附圖和示例性實(shí)施例詳細(xì)地描述根據(jù)本發(fā)明的反射轉(zhuǎn)換橫波共檢波點(diǎn)疊加剖面的剩余靜校正方法。根據(jù)本發(fā)明，在地表起伏較大或低、降速帶復(fù)雜的地區(qū)，通過轉(zhuǎn)換橫波和縱波的信息聯(lián)合求取轉(zhuǎn)換波PS共接收點(diǎn)疊加剖面剩余靜校正量，并采用層位自動追蹤的方法來求取橫波和縱波的層位信息，能夠更高效的提取轉(zhuǎn)換波的剩余靜校正量。根據(jù)本發(fā)明示例性實(shí)施例的反射轉(zhuǎn)換橫波共檢波點(diǎn)疊加剖面的剩余靜校正方法包括以下步驟：步驟一，在反射縱波(即PP波)共檢波點(diǎn)疊加剖面上選擇滿足預(yù)定條件的層位線，在層位線拾取多個(gè)種子點(diǎn)且所述種子點(diǎn)位于層位線走向發(fā)生改變之處，然后連接所述多個(gè)種子點(diǎn)生成一條連線(可以為折線或直線)，將所述連線沿著PP波共檢波點(diǎn)疊加剖面表示時(shí)間的縱軸上下移動預(yù)定時(shí)間后構(gòu)成第一時(shí)窗，然后在第一時(shí)窗內(nèi)進(jìn)行反射縱波層位的自動追蹤，得到反射縱波層位線。這里，預(yù)定條件表示在PP波共檢波點(diǎn)疊加剖面上按照處理項(xiàng)目的任務(wù)要求，目標(biāo)勘探層的時(shí)間附近，選擇的層位線所對應(yīng)的同相軸連續(xù)清晰。預(yù)定時(shí)間的選取原則：一個(gè)或者幾個(gè)子波長度，或者選取層位線所在構(gòu)造的頂部到構(gòu)造底部之差的二分之一，即層位線在縱向上位于最高處(取最大值時(shí))的時(shí)間和位于最低處(取最小值時(shí))的時(shí)間之差的1/2，并且，在下文中形成第二時(shí)窗和第三時(shí)窗時(shí)的預(yù)定時(shí)間的選擇原則與形成第一時(shí)窗時(shí)的預(yù)定時(shí)間選取原則一樣，在下文不再贅述。另外，第二時(shí)窗和第三時(shí)窗不一定與第一時(shí)窗取一樣大，可以小一些的，比如選擇一個(gè)子波長度。另外，一個(gè)層位就是地下某一時(shí)期沉積而形成的反射界面，在視覺上從左到右的相鄰地震道波峰或波谷連續(xù)拾取而獲得的曲線，而由于沉積的復(fù)雜性、后期構(gòu)造的破壞以及轉(zhuǎn)換波采集技術(shù)的原因，通常一些層位連續(xù)性不佳，因此全部人工的判斷和拾取是非常困難的。本發(fā)明對此進(jìn)行了改善僅僅需要拾取幾個(gè)種子點(diǎn)，將各種子點(diǎn)連接，兩種子點(diǎn)之間的點(diǎn)采用線性插值，得到一條折線(或者直線)。然后用以此折線為基準(zhǔn)線，沿著時(shí)間軸上下拓展預(yù)定時(shí)間形成時(shí)窗，并在該時(shí)窗內(nèi)進(jìn)行層位的自動追蹤。其中，PP波層位的自動追蹤通過以下三種方式之一實(shí)現(xiàn)：第一種方式：拾取第一時(shí)窗內(nèi)每一道地震數(shù)據(jù)的最大振幅能量所對應(yīng)的時(shí)間并投影到PP波共檢波點(diǎn)疊加剖面上，得到PP波層位線。第二種方式：采用互相關(guān)法求取。具體地，將第一時(shí)窗內(nèi)的第i道地震數(shù)據(jù)與第i+1道地震數(shù)據(jù)做互相關(guān)運(yùn)算，i從1取到N-1且i為正整數(shù)，N為地震道總數(shù)，然后拾取互相關(guān)運(yùn)算后的每一道地震數(shù)據(jù)的最大相關(guān)量所對應(yīng)的時(shí)間并投影到PP波共檢波點(diǎn)疊加剖面上，得到PP波層位線。第三種方式：采用種子道互相關(guān)法求取。具體地，將第一時(shí)窗內(nèi)的第i-n道至第i+n道地震數(shù)據(jù)進(jìn)行疊加得到種子道，然后將種子道與第i道地震數(shù)據(jù)做互相關(guān)運(yùn)算，i從n取到N-n且i為正整數(shù)，1≤n≤N且為正整數(shù)，N為地震道總數(shù)，然后拾取互相關(guān)運(yùn)算后的每一道地震數(shù)據(jù)的最大相關(guān)量所對應(yīng)的時(shí)間并投影到PP波共檢波點(diǎn)疊加剖面上，得到PP波層位線(或稱為PP波層位)。其中，假設(shè)PP波共檢波點(diǎn)疊加剖面的總道數(shù)為N，以各種子點(diǎn)的連線為基準(zhǔn)線，沿時(shí)間軸上下各移動預(yù)定時(shí)間T形成的第一時(shí)窗，第一時(shí)窗內(nèi)包括N*(2T+1)的地震數(shù)據(jù)。步驟二，根據(jù)給定的PS波層位與PP波層位的速度比值將步驟一得到的PP波層位線投影到反射轉(zhuǎn)換橫波(PS波)共檢波點(diǎn)疊加剖面上得到PS波的比擬層位線，即通過將PP波共檢波點(diǎn)疊加剖面上自動計(jì)算出的PP波層位線和人工給予的速度比值相乘得到PS波的比擬層位線。然后，將PS波的比擬層位線沿著PS波共檢波點(diǎn)疊加剖面表示時(shí)間的縱軸上下移動預(yù)定時(shí)間后構(gòu)成第二時(shí)窗，然后在第二時(shí)窗內(nèi)進(jìn)行PS波層位的自動追蹤，得到第一PS波層位線(或稱為PS波層位)。這里，速度比值可以在PP波共檢波點(diǎn)疊加剖面表示地震道的橫軸上分段取不同值，或者在PP波共檢波點(diǎn)疊加剖面所有地震道上取同一個(gè)值。即可以利用某個(gè)固定的或者分段變化的速度比值將PP波層位投影到PS波共檢波點(diǎn)剖面上。速度比值一般先采用經(jīng)驗(yàn)值，給一個(gè)常數(shù)，或者給一組數(shù)，一組數(shù)是指道數(shù)與速度比值匹配的數(shù)；當(dāng)然也可以采用粗略拾取出來的PS波層位除以PP波層位，得到每一道的商，對商的集合進(jìn)行大尺寸的平滑，平滑后的值就是每一道的速度比值。后者比前者精確但是復(fù)雜，在實(shí)際應(yīng)用中通常就用的第一種方法根據(jù)經(jīng)驗(yàn)值給定速度比值。另外，在第二時(shí)窗內(nèi)進(jìn)行的PS波層位的自動追蹤可以采用以上三種方式之一：第一種方式：拾取所述第二時(shí)窗內(nèi)每一道地震數(shù)據(jù)的最大振幅能量所對應(yīng)的時(shí)間并投影到PS波共檢波點(diǎn)疊加剖面上，得到所述第一PS波層位線。第二種方式：采用互相關(guān)法求取。將第二時(shí)窗內(nèi)的第i道地震數(shù)據(jù)與第i+1道地震數(shù)據(jù)做互相關(guān)運(yùn)算，i從1取到N-1且i為正整數(shù)，N為地震道總數(shù)，然后拾取互相關(guān)運(yùn)算后的每一道地震數(shù)據(jù)的最大相關(guān)量所對應(yīng)的時(shí)間并投影到PS波共檢波點(diǎn)疊加剖面上，得到第一PS波層位線；第三種方式：采用種子道互相關(guān)法求取。將第二時(shí)窗內(nèi)的第i-n道至第i+n道地震數(shù)據(jù)進(jìn)行疊加得到種子道，然后將所述種子道與第i道地震數(shù)據(jù)做互相關(guān)運(yùn)算，i從n取到N-n且i為正整數(shù)，1≤n≤N且為正整數(shù)，N為地震道總數(shù)，拾取互相關(guān)運(yùn)算后的每一道地震數(shù)據(jù)的最大相關(guān)量所對應(yīng)的時(shí)間并投影到PS波共檢波點(diǎn)疊加剖面上，得到第一PS波層位線。步驟三，將第一PS波層位線和PS波的比擬層位線之間的時(shí)間差確定為PS波的初次剩余靜校正量。步驟四，利用步驟三得到的PS波的初次剩余靜校正量對PS波共檢波點(diǎn)疊加剖面進(jìn)行初次剩余靜校正，得到初次剩余靜校正后的臨時(shí)剖面。步驟五，在初次剩余靜校正后的臨時(shí)剖面上，將第一PS波層位線沿著初次剩余靜校正后的臨時(shí)剖面表示時(shí)間的縱軸上下移動預(yù)定時(shí)間后構(gòu)成第三時(shí)窗，然后在第三時(shí)窗內(nèi)進(jìn)行PS波層位的自動追蹤，得到第二PS波層位線。這里，在第三時(shí)窗內(nèi)進(jìn)行的反射轉(zhuǎn)換橫波層位的自動追蹤步驟可以采用以下三種方式之一實(shí)現(xiàn)：方式1：拾取所述第三時(shí)窗內(nèi)每一道地震數(shù)據(jù)的最大振幅能量所對應(yīng)的時(shí)間并投影到所述初次剩余靜校正后的臨時(shí)剖面上，得到所述第二反射轉(zhuǎn)換橫波層位線。方式2：采用互相關(guān)法求取。將所述第三時(shí)窗內(nèi)的第i道地震數(shù)據(jù)與第i+1道地震數(shù)據(jù)做互相關(guān)運(yùn)算，i從1取到N-1且i為正整數(shù)，N為地震道總數(shù)，然后拾取互相關(guān)運(yùn)算后的每一道地震數(shù)據(jù)的最大相關(guān)量所對應(yīng)的時(shí)間并投影到所述初次剩余靜校正后的臨時(shí)剖面上，得到所述第二反射轉(zhuǎn)換橫波層位線；方式3：采用種子道互相關(guān)法求取。將所述第三時(shí)窗內(nèi)的第i-n道至第i+n道地震數(shù)據(jù)進(jìn)行疊加得到種子道，然后將所述種子道與第i道地震數(shù)據(jù)做互相關(guān)運(yùn)算，i從n取到N-n且i為正整數(shù)，1≤n≤N且為正整數(shù)，N為地震道總數(shù)，然后拾取互相關(guān)運(yùn)算后的每一道地震數(shù)據(jù)的最大相關(guān)量所對應(yīng)的時(shí)間并投影到所述初次剩余靜校正后的臨時(shí)剖面上，得到所述第二反射轉(zhuǎn)換橫波層位線。優(yōu)選地，采用方式2的互相關(guān)法或者是方式3的種子道互相關(guān)法，以使層位拾取更準(zhǔn)確，可以有效的規(guī)避當(dāng)某一道能量變化出現(xiàn)異常導(dǎo)致拾取不準(zhǔn)確的情況。當(dāng)然，當(dāng)層位平緩，并且資料信噪比較高時(shí)，方式1和方式2、方式3的效果基本一致。步驟六，將第二PS波層位線和第一PS波層位線之間的時(shí)間差確定為PS波的第二次剩余靜校正量(或稱為層位抖動量)。步驟七，將步驟三得到的初次剩余靜校正量和步驟六得到的第二次剩余靜校正量進(jìn)行疊加，得到PS波共檢波點(diǎn)疊加剖面的最終剩余靜校正量(短波長剩余靜校正量)。步驟八，利用步驟七得到的PS波共檢波點(diǎn)疊加剖面的最終剩余靜校正量對PS波共檢波點(diǎn)疊加剖面進(jìn)行第二次剩余靜校正。作為上述示例性實(shí)施例的進(jìn)一步優(yōu)選方案，在地表起伏較大或低、降速帶復(fù)雜的地區(qū)，在步驟一和步驟二的層位自動追蹤過程中，采用模型種子道的計(jì)算可以提高相關(guān)分析法的精度，因而優(yōu)選采用種子道相關(guān)法，即第三種方式。作為上述示例性實(shí)施例的進(jìn)一步優(yōu)選方案，在步驟一、二和五的層位自動追蹤過程中，若采用種子道相關(guān)法，還可以對種子道的地震數(shù)據(jù)進(jìn)行求平均處理，即將第i-n到i+n道的地震數(shù)據(jù)疊加后求平均(再除以總道數(shù))得到種子道。作為上述示例性實(shí)施例的進(jìn)一步優(yōu)選方案，還根據(jù)PS波共檢波點(diǎn)疊加剖面質(zhì)量，在PP波共檢波點(diǎn)疊加剖面上分層分段的拾取同相軸相對最清晰和連續(xù)性最好的層位，自動計(jì)算PP波和PS波層位。即，在PP波共檢波點(diǎn)疊加剖面表示時(shí)間的縱軸上分層和/或在表示地震道的橫軸上分段選擇滿足預(yù)定條件的層位線，并分別求出所選各層位線對應(yīng)的PS波共檢波點(diǎn)的剩余靜校正量，即分層分段選擇層位線，并分別進(jìn)行上述步驟一至步驟八。例如，可以在步驟一和步驟二對層位走向種子點(diǎn)的拾取，可以選取分層分段拾取，原則是選取同相軸清晰和連續(xù)的層位。作為上述示例性實(shí)施例的進(jìn)一步優(yōu)選方案，還可以在步驟一和步驟二之前，對PP波共檢波點(diǎn)疊加剖面和PS波共檢波點(diǎn)疊加剖面進(jìn)行去噪和振幅預(yù)處理，以突出同相軸。這里，去噪和振幅預(yù)處理是本領(lǐng)域的常規(guī)處理步驟，即是對疊加剖面進(jìn)行一種質(zhì)量提高的處理手段，提高疊加剖面的分辨率。根據(jù)本發(fā)明，①不需要在人工在縱波和轉(zhuǎn)換波共檢波點(diǎn)疊加剖面上分段拾取滿足預(yù)定條件的層位，而是采用上述步驟一和步驟二，與以往的剩余靜校正方法相比，不需要提供準(zhǔn)確的層位線，僅僅需要提供幾個(gè)構(gòu)造走向(層位走向)位置信息，就自動計(jì)算了PP波和PS波的層位。②不需要提供PP波和PS波的層位段或者層位的拾取來計(jì)算速度比。通過給定一個(gè)常量或者分段變化的速度比值，再自動進(jìn)行計(jì)算，計(jì)算結(jié)果更準(zhǔn)確，剩余靜校正量的計(jì)算迭代次數(shù)減少，提高了剩余靜校正量提取的效率。③本發(fā)明在計(jì)算得到初次剩余靜校正量之后，再次通過步驟七計(jì)算層位的抖動量，能夠消除由于橫波反射信號較弱，頻率較低，難以同噪音和縱波反射準(zhǔn)確區(qū)分開，層位拾取不精確，影響剩余靜校正量計(jì)算的準(zhǔn)確性問題。為了更好地理解本發(fā)明的上述示例性實(shí)施例，下面結(jié)合具體示例對其進(jìn)行進(jìn)一步說明。示例1圖1是本發(fā)明示例性實(shí)施例的反射轉(zhuǎn)換橫波共檢波點(diǎn)疊加剖面的剩余靜校正方法的在反射縱波共檢波點(diǎn)疊加剖面選取的層位線上拾取種子點(diǎn)步驟的示意圖。圖2是本發(fā)明示例性實(shí)施例的反射轉(zhuǎn)換橫波共檢波點(diǎn)疊加剖面的剩余靜校正方法的構(gòu)成第一時(shí)窗步驟的示意圖。參照圖1，本示例的反射轉(zhuǎn)換橫波共檢波點(diǎn)疊加剖面的剩余靜校正方法包括以下步驟：(1)首先，在PP波共檢波點(diǎn)疊加剖面Stack1上(該剖面有200道數(shù)據(jù))，選擇一條同相軸連續(xù)清晰的層位線1，根據(jù)所選層位線構(gòu)造走向(層位線走向)發(fā)生改變之處拾取了5個(gè)種子點(diǎn)2、3、4、5、6，即在圖1層位線1的構(gòu)造底部和頂部拾取五個(gè)點(diǎn)作為種子點(diǎn)。如圖2所示，將各種子點(diǎn)2、3、4、5、6連接，兩種子點(diǎn)之間的點(diǎn)采用線性插值，得到一條折線9。以此折線9為基準(zhǔn)線，沿著時(shí)間軸向上移動T(T＝50ms)形成第一時(shí)窗的上邊界線7，沿著時(shí)間軸向下移動T(T＝50ms)形成第一時(shí)窗的下邊界線8，也就是說，沿著時(shí)間軸上下各取T＝50ms的數(shù)據(jù)，形成了新的200道*101ms的新數(shù)據(jù)。求取新數(shù)據(jù)(即第一時(shí)窗內(nèi)的數(shù)據(jù))的每一道的最大絕對值(最大振幅能量)對應(yīng)的位置Time，將Time記錄下來并投影到PP波共檢波點(diǎn)疊加剖面上，就得到自動追蹤的PP波層位HorPP。(2)在所有地震道利用固定的速度比值1.9將步驟(1)拾取的HorPP層位投影到PS波共檢波點(diǎn)剖面上得到HorPS比擬。在PS波共檢波點(diǎn)疊加剖面上Stack2上(該剖面有200道數(shù)據(jù))，將HorPS比擬沿著PS波共檢波點(diǎn)疊加剖面表示時(shí)間的縱軸上下移動30ms時(shí)間后構(gòu)成第二時(shí)窗，然后在第二時(shí)窗內(nèi)進(jìn)行PS波層位的自動追蹤。即以HorPS比擬為中心，上下各取30ms的數(shù)據(jù)，形成了新的200道*61ms的新數(shù)據(jù)，求取新數(shù)據(jù)(即第二時(shí)窗內(nèi)的數(shù)據(jù))的每一道的最大絕對值(最大振幅能量)對應(yīng)位置Time，將Time記錄下來并投影到PS波共檢波點(diǎn)疊加剖面上，就得到自動追蹤的PS波層位HorPS。(3)將HorPS和HorPS比擬之差記做PS波初次剩余靜校正量Statics1。(4)利用初次的剩余靜校正量Statics1，對PS波共檢波點(diǎn)疊加剖面進(jìn)行初次剩余靜校正，得到初次靜校正后的臨時(shí)剖面Stack3。(5)在臨時(shí)剖面上Stack3通過PS波層位自動追蹤的互相關(guān)法求取出層位抖動Statics2。具體地，在初次剩余靜校正后的臨時(shí)剖面Stack3上，將PS波層位HorPS沿著臨時(shí)剖面Stack3表示時(shí)間的縱軸上下移動25ms時(shí)間后構(gòu)成第三時(shí)窗，然后在第三時(shí)窗內(nèi)進(jìn)行PS波層位的自動追蹤，即將第三時(shí)窗內(nèi)的第i道地震數(shù)據(jù)與第i+1道地震數(shù)據(jù)做互相關(guān)運(yùn)算，i從1取到N-1(即200-1＝199)且i為正整數(shù)，N為地震道總數(shù)，N＝200，然后拾取互相關(guān)運(yùn)算后的每一道地震數(shù)據(jù)的最大相關(guān)量所對應(yīng)的時(shí)間并投影到初次剩余靜校正后的臨時(shí)剖面Stack3上，得到PS波層位線HorPS′。將PS波層位線HorPS′和PS波層位線HorPS之間的時(shí)間差確定為PS波的第二次剩余靜校正量(或稱為層位抖動量)Statics2。(6)將步驟(5)中求得的層位抖動量Statics2和步驟(3)求得的初次剩余靜校正量Statics1進(jìn)行疊加，得到當(dāng)前PS波共檢波點(diǎn)疊加剖面的最終剩余靜校正量。(7)利用步驟(6)得到的PS波共檢波點(diǎn)疊加剖面的最終剩余靜校正量對PS波共檢波點(diǎn)疊加剖面進(jìn)行第二次剩余靜校正。示例2(1)首先，與示例1類似的，在PP波共檢波點(diǎn)疊加剖面Stack1上(該剖面有200道數(shù)據(jù))，選擇一條同相軸連續(xù)清晰的層位線1，根據(jù)所選層位線構(gòu)造走向發(fā)生改變之處拾取了5個(gè)種子點(diǎn)。將各種子點(diǎn)連接，兩種子點(diǎn)之間的點(diǎn)采用線性插值，得到一條折線。以此折線為基準(zhǔn)線，沿著時(shí)間軸向上移動T(T＝100ms)形成第一時(shí)窗的上邊界線，沿著時(shí)間軸向下移動T(T＝100ms)形成第一時(shí)窗的下邊界線，也就是說，沿著時(shí)間軸上下各取T＝100ms的數(shù)據(jù)，形成了新的200道*201ms的新數(shù)據(jù)。在第一時(shí)窗內(nèi)利用互相關(guān)法求取PP波層位線。具體地，將第一時(shí)窗內(nèi)的第i道地震數(shù)據(jù)與第i+1道地震數(shù)據(jù)做互相關(guān)運(yùn)算，i從1取到N-1且i為正整數(shù)，N為地震道總數(shù)，N＝200，然后拾取互相關(guān)運(yùn)算后的每一道地震數(shù)據(jù)的最大相關(guān)量所對應(yīng)的時(shí)間并投影到PP波共檢波點(diǎn)疊加剖面上，得到自動追蹤的PP波層位線HorPP。(2)根據(jù)下表1給出的PS波層位與PP波層位的速度比值將步驟(1)得到的PP波層位線HorPP投影到PS波共檢波點(diǎn)疊加剖面上得到PS波的比擬層位線HorPS比擬。表1分段速度比值道1～5051～100101～150151～200速度比值1.91.61.61.9以HorPS比擬為中心(即基準(zhǔn)線)，上下各取50ms的數(shù)據(jù)，形成了新的200道*101ms的新數(shù)據(jù)。對PS波層位的自動追蹤采用互相關(guān)法求取，在形成了新的200道*101ms的新數(shù)據(jù)中，將第i道數(shù)據(jù)與第i+1道數(shù)據(jù)做互相關(guān)(i從1取到N-1道)，得到互相關(guān)結(jié)果組成的新的數(shù)據(jù)，求取互相關(guān)后的新數(shù)據(jù)的每一道的最大絕對值(最大相關(guān)量)對應(yīng)是位置Time，將Time記錄下來并投影到PS波共檢波點(diǎn)疊加剖面上，就得到自動追蹤的HorPS。(3)將HorPS和HorPS比擬之差記做PS波初次剩余靜校正量Statics1。(4)利用初次的剩余靜校正量Statics1，對PS波共檢波點(diǎn)疊加剖面進(jìn)行初次剩余靜校正，得到初次靜校正后的臨時(shí)剖面Stack3。(5)在臨時(shí)剖面上Stack3通過PS波層位自動追蹤的互相關(guān)法求取出層位抖動Statics2。該步驟與示例1的操作相同，在此不再贅述。(6)將步驟(5)中求得的層位抖動量Statics2和步驟(3)求得的初次剩余靜校正量Statics1進(jìn)行疊加，得到當(dāng)前PS波共檢波點(diǎn)疊加剖面的最終剩余靜校正量。(7)利用步驟(6)得到的PS波共檢波點(diǎn)疊加剖面的最終剩余靜校正量對PS波共檢波點(diǎn)疊加剖面進(jìn)行第二次剩余靜校正。示例3(1)與示例1和示例2采用相同操作，即在PP波共檢波點(diǎn)疊加剖面Stack1上取5個(gè)種子點(diǎn)并連線，然后拓展形成第一時(shí)窗。在第一時(shí)窗內(nèi)利用種子道互相關(guān)法求取PP波層位線。具體地，將第i-5道到第i+5道數(shù)據(jù)進(jìn)行疊加得到種子道，再將第i道數(shù)據(jù)與種子道做互相關(guān)，i＝n到N-n道，即i從5取到195道，得到互相關(guān)結(jié)果組成的新的數(shù)據(jù)，求取新數(shù)據(jù)的每一道的最大絕對值(最大相關(guān)量)對應(yīng)是位置Time，將Time記錄將下來并投影到PP波共檢波點(diǎn)疊加剖面上，就得到自動追蹤的HorPP。(2)根據(jù)表1給出的PS波層位與PP波層位的速度比值將步驟(1)得到的PP波層位線HorPP投影到PS波共檢波點(diǎn)疊加剖面上得到PS波的比擬層位線HorPS比擬。以HorPS比擬為中心(即基準(zhǔn)線)，上下各取30ms的數(shù)據(jù)，形成了新的200道*61ms的新數(shù)據(jù)。對PS波層位的自動追蹤采用種子道互相關(guān)法求取，在形成了新的200道*61ms的新數(shù)據(jù)中，將第i-5道數(shù)據(jù)與第i+5道數(shù)據(jù)進(jìn)行疊加得到種子道，再將第i道數(shù)據(jù)與種子道做互相關(guān)，i＝n到N-n道，即i從5取到195道，得到互相關(guān)結(jié)果組成的新的數(shù)據(jù)，求取新數(shù)據(jù)的每一道的最大絕對值(最大相關(guān)量)對應(yīng)是位置Time，將Time記錄將下來并投影到PS波共檢波點(diǎn)疊加剖面上，就得到自動追蹤的HorPS。(3)將HorPS和HorPS比擬之差記做PS波初次剩余靜校正量Statics1。(4)利用初次的剩余靜校正量Statics1，對PS波共檢波點(diǎn)疊加剖面進(jìn)行初次剩余靜校正，得到初次靜校正后的臨時(shí)剖面Stack3。(5)在臨時(shí)剖面上Stack3通過PS波層位自動追蹤的互相關(guān)法求取出層位抖動Statics2。該步驟與示例1的操作相同，在此不再贅述。(6)將步驟(5)中求得的層位抖動量Statics2和步驟(3)求得的初次剩余靜校正量Statics1進(jìn)行疊加，得到當(dāng)前PS波共檢波點(diǎn)疊加剖面的最終剩余靜校正量。(7)利用步驟(6)得到的PS波共檢波點(diǎn)疊加剖面的最終剩余靜校正量對PS波共檢波點(diǎn)疊加剖面進(jìn)行第二次剩余靜校正。綜上所述，本發(fā)明提供了一種基于自動拾取層位的轉(zhuǎn)換波剩余靜校正方法，在地表起伏較大或低、降速帶復(fù)雜的地區(qū)，通過轉(zhuǎn)換橫波和縱波的信息聯(lián)合求取轉(zhuǎn)換波PS共接收點(diǎn)疊加剖面剩余靜校正量，并采用層位自動追蹤的方式來求取轉(zhuǎn)換橫波和縱波的層位信息，能夠更高效的提取轉(zhuǎn)換波的剩余靜校正量。盡管上面已經(jīng)通過結(jié)合示例性實(shí)施例描述了本發(fā)明，但是本領(lǐng)域技術(shù)人員應(yīng)該清楚，在不脫離權(quán)利要求所限定的精神和范圍的情況下，可對本發(fā)明的示例性實(shí)施例進(jìn)行各種修改和改變。