轉(zhuǎn)換波共檢波點(diǎn)疊加靜校正方法
【專利摘要】本發(fā)明公開一種轉(zhuǎn)換波共檢波點(diǎn)疊加靜校正方法，包括以下步驟：（a）在縱波和轉(zhuǎn)換波共檢波點(diǎn)疊加剖面上分段拾取滿足預(yù)定條件的層位；（b）執(zhí)行自動(dòng)層拉平，并在設(shè)置的時(shí)窗內(nèi)自動(dòng)精細(xì)拾取相應(yīng)的層位段；（c）計(jì)算拾取的層位段的速度比值；（d）以所述速度比值將縱波剖面上拾取的層位段投影到轉(zhuǎn)換波剖面，并將轉(zhuǎn)換波剖面上拾取的層位段與投影到轉(zhuǎn)換波剖面上的層位段之間的時(shí)間差確定為檢波點(diǎn)靜校正量；（e）將確定的檢波點(diǎn)靜校正量應(yīng)用到轉(zhuǎn)換波剖面，以便對(duì)檢波點(diǎn)進(jìn)行靜校正。
【專利說(shuō)明】轉(zhuǎn)換波共檢波點(diǎn)疊加靜校正方法
【技術(shù)領(lǐng)域】
[0001]本發(fā)明涉及石油天然氣地震勘探領(lǐng)域，更具體地講，涉及一種能夠精確地獲取轉(zhuǎn)換波的檢波點(diǎn)靜校正量的基于迭代的最優(yōu)轉(zhuǎn)換波共檢波點(diǎn)疊加靜校正方法。
【背景技術(shù)】
[0002]近年來(lái)，多波勘探得到了越來(lái)越多的重視。特別是多波勘探在巖性分析、流體識(shí)另IJ、裂縫檢測(cè)、氣云區(qū)成像等方面顯示出其獨(dú)特優(yōu)勢(shì)，在實(shí)際工業(yè)生產(chǎn)在得到了越來(lái)越多的重視。轉(zhuǎn)換波處理與傳統(tǒng)的縱波處理有很多不同，其關(guān)鍵步驟通常是靜校正量的求取。轉(zhuǎn)換波靜校正量是由炮點(diǎn)的縱波靜校正量和檢波點(diǎn)的橫波靜校正量組成，轉(zhuǎn)換波炮點(diǎn)靜校正可以用傳統(tǒng)方法求得，而檢波點(diǎn)橫波靜校正則需要專門的方法進(jìn)行求取，這是因?yàn)闄M波與縱波在近地表傳播速度差異極大，相同的位置上橫波靜校正量通常比縱波靜校正量大2?10倍。與縱波相比，橫波實(shí)際上不受近地表潛水面的影響，因此縱波和橫波靜校正量不存在簡(jiǎn)單的比例關(guān)系，使用縱波靜校正量乘以某一系數(shù)來(lái)對(duì)橫波靜校正量進(jìn)行近似通常是不能用于精細(xì)成像的。
[0003]檢波點(diǎn)橫波靜校正方法可以分為兩大類:一類是基于建立速度模型再求取靜校正量的方法，另一類是基于數(shù)據(jù)直接得到靜校正量的方法。模型法主要依賴于近地表速度模型?；诘卣饠?shù)據(jù)本身的靜校正方法，即，共檢波點(diǎn)疊加相干法，雖然在求取靜校正量的方法上得到了最好的處理效果，但需要很多的手工操作，不利于規(guī)模應(yīng)用，而且僅適用于構(gòu)造較為簡(jiǎn)單的情況。
【發(fā)明內(nèi)容】
[0004]因此，本發(fā)明的目的在于提供一種基于迭代的最優(yōu)轉(zhuǎn)換波共檢波點(diǎn)疊加靜校正方法，從而為轉(zhuǎn)換波處理提供更精細(xì)和更高效的靜校正處理方法。
[0005]根據(jù)本發(fā)明的一方面，提供一種轉(zhuǎn)換波共檢波點(diǎn)疊加靜校正方法，包括以下步驟:(a)在縱波和轉(zhuǎn)換波共檢波點(diǎn)疊加剖面上分段拾取滿足預(yù)定條件的層位；(b)執(zhí)行自動(dòng)層拉平，并在設(shè)置的時(shí)窗內(nèi)自動(dòng)精細(xì)拾取相應(yīng)的層位段；(C)計(jì)算拾取的層位段的速度比值；Cd)以所述速度比值將縱波剖面上拾取的層位段投影到轉(zhuǎn)換波剖面，并將轉(zhuǎn)換波剖面上拾取的層位段與投影到轉(zhuǎn)換波剖面上的層位段之間的時(shí)間差確定為檢波點(diǎn)靜校正量；(e)將確定的檢波點(diǎn)靜校正量應(yīng)用到轉(zhuǎn)換波剖面，以便對(duì)檢波點(diǎn)進(jìn)行靜校正。
[0006]優(yōu)選地，所述轉(zhuǎn)換波共檢波點(diǎn)疊加靜校正方法還包括以下步驟:在步驟(a)之前執(zhí)行縱波和轉(zhuǎn)換波共檢波點(diǎn)疊加。
[0007]優(yōu)選地，步驟(e)包括:(el)將確定的檢波點(diǎn)靜校正量從低頻至高頻劃分為多個(gè)頻率成分；(e2)按照從低頻至高頻的順序?qū)澐值臋z波點(diǎn)靜校正量的多個(gè)頻率成分依次應(yīng)用到轉(zhuǎn)換波剖面，其中，在步驟(e2)中，每當(dāng)應(yīng)用劃分的檢波點(diǎn)靜校正量的一個(gè)頻率成分之后，判斷檢波點(diǎn)靜校正量是否穩(wěn)定，如果檢波點(diǎn)靜校正量不穩(wěn)定，則應(yīng)用劃分的檢波點(diǎn)靜校正量的下一個(gè)頻率成分，如果檢波點(diǎn)靜校正量穩(wěn)定，則終止所述轉(zhuǎn)換波共檢波點(diǎn)疊加靜校正方法，其中，如果在應(yīng)用劃分的檢波點(diǎn)靜校正量的最高頻率成分之后，檢波點(diǎn)靜校正量仍然不穩(wěn)定，則返回步驟(a)，以進(jìn)行下一輪的檢波點(diǎn)靜校正。
[0008]優(yōu)選地,如果檢波點(diǎn)靜校正量的當(dāng)前頻率成分與檢波點(diǎn)靜校正量的前一個(gè)頻率成分之間的差值大于預(yù)定閾值，則確定檢波點(diǎn)靜校正量不穩(wěn)定。
[0009]優(yōu)選地，所述預(yù)定條件表示通過(guò)拾取層位段獲得的二維曲線段具有良好的連續(xù)性。
[0010]優(yōu)選地，在步驟(C)中，根據(jù)如下等式來(lái)計(jì)算層位段的速度比值:
[0011]r (CRP) =2T0ps/T0pp-l,
[0012]其中，Ttlps表示與轉(zhuǎn)換波剖面上拾取的層位段對(duì)應(yīng)的時(shí)間，T0pp表示與縱波剖面上拾取的層位段對(duì)應(yīng)的時(shí)間，r (CRP)表示拾取的層位段的速度比值。
[0013]優(yōu)選地，對(duì)層位段的速度比值r (CRP)進(jìn)行空間上的平滑處理，來(lái)消除空間非構(gòu)造因素的干擾。
【專利附圖】
【附圖說(shuō)明】
[0014]通過(guò)下面結(jié)合附圖對(duì)實(shí)施例進(jìn)行的描述，本發(fā)明的這些和/或其他方面和優(yōu)點(diǎn)將會(huì)變得清楚和更易于理解，其中:
[0015]圖1是示出根據(jù)本發(fā)明的示例性實(shí)施例的轉(zhuǎn)換波共檢波點(diǎn)疊加靜校正方法的流程圖。
【具體實(shí)施方式】
[0016]現(xiàn)在對(duì)本發(fā)明實(shí)施例進(jìn)行詳細(xì)的描述，其示例表示在附圖中，其中，相同的標(biāo)號(hào)始終表示相同部件。下面通過(guò)參照附圖對(duì)實(shí)施例進(jìn)行描述以解釋本發(fā)明。
[0017]圖1是示出根據(jù)本發(fā)明的示例性實(shí)施例的轉(zhuǎn)換波共檢波點(diǎn)疊加靜校正方法的流程圖。
[0018]參照?qǐng)D1，在步驟SlOl中，執(zhí)行縱波和轉(zhuǎn)換波共檢波點(diǎn)疊加。這里，可根據(jù)現(xiàn)有技術(shù)容易地執(zhí)行縱波和轉(zhuǎn)換波共檢波點(diǎn)疊加。
[0019]在步驟S102中，在縱波和轉(zhuǎn)換波共檢波點(diǎn)疊加剖面上分段拾取滿足預(yù)定條件的層位，即，針對(duì)共檢波點(diǎn)，同時(shí)在縱波剖面和轉(zhuǎn)換波剖面上分段拾取相應(yīng)的層位。這里，對(duì)于構(gòu)造簡(jiǎn)單的地震剖面，一個(gè)層位就是地下某一時(shí)期沉積而形成的反射界面，在視覺(jué)上是從左到右的一個(gè)波峰或波谷連續(xù)拾取而獲得的曲線，而由于沉積的復(fù)雜性、后期構(gòu)造的破壞以及轉(zhuǎn)換波采集技術(shù)的原因，通常一些層位連續(xù)性不佳。因此，本發(fā)明提出“分段”的概念，也就是拾取不同層位(時(shí)間)在空間上連續(xù)性較好的段。在這種情況下，層位段就是一個(gè)時(shí)間-空間的二維曲線段。所述預(yù)定條件表示通過(guò)拾取層位段獲得的二維曲線段具有良好的連續(xù)性。因此，在步驟S102中，可以在縱波和轉(zhuǎn)換波共檢波點(diǎn)疊加剖面上拾取具有良好連續(xù)性的層位段。
[0020]此外，由于一個(gè)剖面在同一空間區(qū)域從淺到深會(huì)有很多連續(xù)的反射同相軸，因此每一個(gè)反射同相軸對(duì)應(yīng)一個(gè)層位。在這種情況下，拾取的層位段可以不是同一層位。
[0021]接下來(lái)，在步驟S103中，執(zhí)行自動(dòng)層拉平，并在設(shè)置的時(shí)窗內(nèi)自動(dòng)精細(xì)拾取相應(yīng)的層位段。由于在步驟S102中拾取的層位段一般都是扭曲的，特別是對(duì)于轉(zhuǎn)換波來(lái)說(shuō)，上下跳動(dòng)比較嚴(yán)重。因此，可通過(guò)根據(jù)現(xiàn)有技術(shù)執(zhí)行自動(dòng)層拉平，通過(guò)整體上下移動(dòng)拾取的地震道，使得拾取的層位點(diǎn)都在同一水平面上。通過(guò)層拉平以后，轉(zhuǎn)換波本身上下跳動(dòng)比較嚴(yán)重的情況可以得到改善。因此，在此基礎(chǔ)上給定一個(gè)較小的時(shí)窗范圍進(jìn)行精細(xì)拾取，可以獲得更好的成像效果，以便更好地進(jìn)行后續(xù)處理。
[0022]其后，在步驟S104中，計(jì)算拾取的層位段的速度比值。具體地講，可根據(jù)如下等式(I)來(lái)計(jì)算層位段的速度比值:
[0023]r (CRP) =2T0ps/T0pp-l(I),
[0024]其中，Ttlps表示與轉(zhuǎn)換波剖面上拾取的層位段對(duì)應(yīng)的時(shí)間，T0pp表示與縱波剖面上拾取的層位段對(duì)應(yīng)的時(shí)間，r (CRP)表示空間-時(shí)間位置的速度比值，即，拾取的層位段的速度比值。這里，r (CRP)需要通過(guò)空間上的平滑處理，來(lái)消除空間非構(gòu)造因素(主要是靜校正)的干擾，以便后面計(jì)算的靜校正更準(zhǔn)確。
[0025]接著，在步驟S105中，以所述速度比值將縱波剖面上拾取的層位段投影到轉(zhuǎn)換波剖面，并將轉(zhuǎn)換波剖面上拾取的層位段與投影到轉(zhuǎn)換波剖面上的層位段之間的時(shí)間差確定為檢波點(diǎn)靜校正量。
[0026]然后，在步驟S106中，將確定的檢波點(diǎn)靜校正量應(yīng)用到轉(zhuǎn)換波剖面，以便對(duì)檢波點(diǎn)進(jìn)行靜校正。具體地講，在步驟S106中，首先將確定的檢波點(diǎn)靜校正量從低頻至高頻劃分為多個(gè)頻率成分，然后按照從低頻至高頻的順序?qū)澐值臋z波點(diǎn)靜校正量的多個(gè)頻率成分依次應(yīng)用到轉(zhuǎn)換波剖面，以便對(duì)檢波點(diǎn)進(jìn)行靜校正。這里，可通過(guò)橫向平滑方法來(lái)提取檢波點(diǎn)靜校正量的低頻成分， 并通過(guò)從原始值減去經(jīng)橫向平滑得到的低頻值來(lái)提取檢波點(diǎn)靜校正量的中高頻成分。更具體地，每當(dāng)應(yīng)用劃分的檢波點(diǎn)靜校正量的一個(gè)頻率成分之后，需要判斷檢波點(diǎn)靜校正量是否穩(wěn)定。如果檢波點(diǎn)靜校正量不穩(wěn)定，則應(yīng)用劃分的檢波點(diǎn)靜校正量的下一個(gè)頻率成分，如果檢波點(diǎn)靜校正量穩(wěn)定，則結(jié)束所述轉(zhuǎn)換波共檢波點(diǎn)疊加靜校正方法。此外，如果在應(yīng)用劃分的檢波點(diǎn)靜校正量的最高頻率成分之后，檢波點(diǎn)靜校正量仍然不穩(wěn)定，則返回步驟S102，以進(jìn)行下一輪的檢波點(diǎn)靜校正，從而進(jìn)行更多的迭代處理。這里，如果檢波點(diǎn)靜校正量的當(dāng)前頻率成分與檢波點(diǎn)靜校正量的前一個(gè)頻率成分之間的差值小于預(yù)定閾值，則可確定檢波點(diǎn)靜校正量穩(wěn)定，如果所述差值大于預(yù)定閾值，則可確定檢波點(diǎn)靜校正量不穩(wěn)定。
[0027]根據(jù)本發(fā)明的示例性實(shí)施例的轉(zhuǎn)換波共檢波點(diǎn)疊加靜校正方法能夠應(yīng)用于二維和三維轉(zhuǎn)換波高精度檢波點(diǎn)靜校正處理，對(duì)油氣、礦產(chǎn)資源勘探有重要應(yīng)用價(jià)值。此外，該方法特別適用于復(fù)雜山地的構(gòu)造轉(zhuǎn)換波靜校正處理，具有計(jì)算步驟簡(jiǎn)單、計(jì)算效率高和處理效果好等特點(diǎn)，在目前西部的復(fù)雜山地地震資料處理中具有廣闊的應(yīng)用前景。
[0028]將理解:可按照硬件、軟件或硬件和軟件的組合的形式來(lái)實(shí)現(xiàn)如上所述的本發(fā)明的示例性實(shí)施例?？稍诜菚簳r(shí)性計(jì)算機(jī)可讀存儲(chǔ)介質(zhì)中存儲(chǔ)任何所述軟件。非暫時(shí)性計(jì)算機(jī)可讀存儲(chǔ)介質(zhì)存儲(chǔ)一個(gè)或多個(gè)程序(軟件模塊)，所述一個(gè)或多個(gè)程序包括指令，其中，當(dāng)由在電子裝置中的一個(gè)或多個(gè)處理器執(zhí)行所述指令時(shí)，所述指令促使電子裝置執(zhí)行本發(fā)明的方法。任何所述軟件可按照易失性或非易失性存儲(chǔ)器(諸如，不論是否是可擦除的或可重寫的存儲(chǔ)裝置，比如只讀存儲(chǔ)器(ROM))的形式被存儲(chǔ)，或者按照存儲(chǔ)器(諸如，隨機(jī)存取存儲(chǔ)器(RAM)、存儲(chǔ)芯片、裝置或集成電路)的形式被存儲(chǔ)，或者在光學(xué)或磁可讀介質(zhì)(諸如，⑶、DVD、磁盤或磁帶等)上被存儲(chǔ)。將理解:存儲(chǔ)裝置和存儲(chǔ)介質(zhì)是適于存儲(chǔ)包括指令的程序的機(jī)器可讀存儲(chǔ)器的實(shí)施例，其中，當(dāng)所述指令被執(zhí)行時(shí)，所述指令實(shí)施本發(fā)明的實(shí)施例。因此，本發(fā)明的實(shí)施例提供一種包括代碼的程序和存儲(chǔ)所述程序的機(jī)器可讀存儲(chǔ)器，其中，所述代碼用于實(shí)施如上所述的根據(jù)本發(fā)明的示例性實(shí)施例的轉(zhuǎn)換波共檢波點(diǎn)疊加靜校正方法。進(jìn)一步地說(shuō)，所述程序可經(jīng)由任何介質(zhì)(諸如通過(guò)有線或無(wú)線連接攜帶的通信信號(hào))被電子地傳輸，并且本發(fā)明的實(shí)施例適當(dāng)?shù)匕龀绦颉?br>
[0029]雖然已經(jīng)顯示和描述了一些實(shí)施例，但是本領(lǐng)域技術(shù)人員應(yīng)該理解，在不脫離本發(fā)明的原理和精神的情況下，可以對(duì)這些實(shí)施例進(jìn)行修改，本發(fā)明的范圍由權(quán)利要求及其等同物限定。
【權(quán)利要求】
1.一種轉(zhuǎn)換波共檢波點(diǎn)疊加靜校正方法，包括以下步驟: Ca)在縱波和轉(zhuǎn)換波共檢波點(diǎn)疊加剖面上分段拾取滿足預(yù)定條件的層位； (b)執(zhí)行自動(dòng)層拉平，并在設(shè)置的時(shí)窗內(nèi)自動(dòng)精細(xì)拾取相應(yīng)的層位段； (C)計(jì)算拾取的層位段的速度比值； Cd)以所述速度比值將縱波剖面上拾取的層位段投影到轉(zhuǎn)換波剖面，并將轉(zhuǎn)換波剖面上拾取的層位段與投影到轉(zhuǎn)換波剖面上的層位段之間的時(shí)間差確定為檢波點(diǎn)靜校正量；(e)將確定的檢波點(diǎn)靜校正量應(yīng)用到轉(zhuǎn)換波剖面，以便對(duì)檢波點(diǎn)進(jìn)行靜校正。
2.根據(jù)權(quán)利要求1所述的轉(zhuǎn)換波共檢波點(diǎn)疊加靜校正方法，還包括以下步驟:在步驟(a)之前執(zhí)行縱波和轉(zhuǎn)換波共檢波點(diǎn)疊加。
3.根據(jù)權(quán)利要求1所述的轉(zhuǎn)換波共檢波點(diǎn)疊加靜校正方法，其中，步驟(e)包括: (el)將確定的檢波點(diǎn)靜校正量從低頻至高頻劃分為多個(gè)頻率成分； (e2)按照從低頻至高頻的順序?qū)澐值臋z波點(diǎn)靜校正量的多個(gè)頻率成分依次應(yīng)用到轉(zhuǎn)換波剖面， 其中，在步驟(e2)中，每當(dāng)應(yīng)用劃分的檢波點(diǎn)靜校正量的一個(gè)頻率成分之后，判斷檢波點(diǎn)靜校正量是否穩(wěn)定，如果檢波點(diǎn)靜校正量不穩(wěn)定，則應(yīng)用劃分的檢波點(diǎn)靜校正量的下一個(gè)頻率成分，如果檢波點(diǎn)靜校正量穩(wěn)定，則終止所述轉(zhuǎn)換波共檢波點(diǎn)疊加靜校正方法， 其中，如果在應(yīng)用劃分的檢波點(diǎn)靜校正量的最高頻率成分之后，檢波點(diǎn)靜校正量仍然不穩(wěn)定，則返回步驟(a)，以進(jìn)行下一輪的檢波點(diǎn)靜校正。
4.根據(jù)權(quán)利要求3所述的轉(zhuǎn)換波共檢波點(diǎn)疊加靜校正方法，其中，如果檢波點(diǎn)靜校正量的當(dāng)前頻率成分與檢波點(diǎn)靜校正量的前一個(gè)頻率成分之間的差值大于預(yù)定閾值，則確定檢波點(diǎn)靜校正量不穩(wěn)定。
5.根據(jù)權(quán)利要求1所述的轉(zhuǎn)換波共檢波點(diǎn)疊加靜校正方法，其中，所述預(yù)定條件表示通過(guò)拾取層位段獲得的二維曲線段具有良好的連續(xù)性。
6.根據(jù)權(quán)利要求1所述的轉(zhuǎn)換波共檢波點(diǎn)疊加靜校正方法，其中，在步驟(c)中，根據(jù)如下等式來(lái)計(jì)算層位段的速度比值:
r (CRP) =2T0ps/T0pp-l， 其中，Ttlps表示與轉(zhuǎn)換波剖面上拾取的層位段對(duì)應(yīng)的時(shí)間，T0pp表示與縱波剖面上拾取的層位段對(duì)應(yīng)的時(shí)間，r (CRP)表示拾取的層位段的速度比值。
7.根據(jù)權(quán)利要求6所述的轉(zhuǎn)換波共檢波點(diǎn)疊加靜校正方法，其中，對(duì)層位段的速度比值r (CRP)進(jìn)行空間上的平滑處理，來(lái)消除空間非構(gòu)造因素的干擾。
【文檔編號(hào)】G01V1/36GK103487834SQ201310408061
【公開日】2014年1月1日 申請(qǐng)日期:2013年9月9日 優(yōu)先權(quán)日:2013年9月9日
【發(fā)明者】羅紅明, 何光明, 李忠, 陳愛(ài)萍, 陳三平, 王棟, 張華
申請(qǐng)人:中國(guó)石油集團(tuán)川慶鉆探工程有限公司地球物理勘探公司