背景技術(shù)：

通常，多波地震勘探利用了地震勘探過程中的多種矢量波場信息，在巖性分析、流體識別、裂縫檢測、氣云區(qū)成像等方面顯示出其獨(dú)特優(yōu)勢，在實(shí)際工業(yè)生產(chǎn)中越來越受到重視。例如，近年來，四川盆地已采集超過5000km²三維多波地震資料，勘探目標(biāo)也從中淺層平緩構(gòu)造拓展到深層具有陡傾地層的復(fù)雜構(gòu)造。目前，陸上多波地震主要是通過炸藥激發(fā)縱波，三分量檢波器同時接收縱波(可簡稱為PP波)和轉(zhuǎn)換橫波(可簡稱為轉(zhuǎn)換波，或PS波)，其中，近地表橫波靜校正處理是影響陸上多波地震勘探數(shù)據(jù)成像的最關(guān)鍵的因素之一。

一般來講，在生產(chǎn)實(shí)踐中，橫波低降速帶橫向變化異常劇烈，與縱波低降速帶在厚度和速度兩方面都具有很大差異，因此，建立近地表橫波速度模型難度較大，接收點(diǎn)處的近地表橫波靜校正量求取誤差也較大?？傮w來講，近地表橫波靜校正量求取大致有三類方法。第一類是通過表層結(jié)構(gòu)調(diào)查建立橫波速度模型而獲得，其主要問題是，野外地表?xiàng)l件復(fù)雜多樣，施工難度大，獲得的橫波資料差，調(diào)查點(diǎn)密度低，導(dǎo)致橫波速度模型精度低；第二類是將縱波靜校正量乘以一固定的比例系數(shù)作為橫波靜校正量，該方法是預(yù)設(shè)一系列比例系數(shù)，分別與已知的縱波靜校正量相乘作為橫波靜校正量，并分別應(yīng)用到轉(zhuǎn)換波疊加處理中，根據(jù)疊加剖面成像的好壞來選擇最佳的比例系數(shù)，該方法不適用于橫向變化的縱、橫波速度模型結(jié)構(gòu)；第三類延遲時靜校正方法需要在單炮上拾取縱波、橫波初至，適用于水平層狀地層，拾取工作量大且橫波初至不易拾取。

上述三種近靜校正量已達(dá)到縱波的2-8倍，嚴(yán)重影響PS轉(zhuǎn)換波地震資料的聚焦成像。目前，針對橫波剩余靜校正量大的問題，常規(guī)方法是采用共接收波點(diǎn)疊加反射層位匹配法求取接收點(diǎn)處大的橫波剩余靜校正量。該方法是通過分別拾取縱波、轉(zhuǎn)換波共接收點(diǎn)全偏移距疊加剖面上同一反射層位并匹配的方法求取接收點(diǎn)處縱橫波速度比，再根據(jù)已知的縱波靜校正量和縱橫波速度比求取大的橫波剩余靜校正量。但是，該方法在遇到陡傾地層或明顯的地層各向異性時，由于轉(zhuǎn)換波的共轉(zhuǎn)換點(diǎn)在傾斜地層界面上非常分散，導(dǎo)致全偏移距道集中數(shù)據(jù)間的反射時差較大，使傾斜地層在疊加剖面上不能成像，因而在疊加剖面上無法拾取相應(yīng)的反射層，導(dǎo)致該方法無法應(yīng)用。

技術(shù)實(shí)現(xiàn)要素：

本發(fā)明的目的在于解決現(xiàn)有技術(shù)存在的上述不足中的至少一項(xiàng)。

本發(fā)明的目的之一在于提供一種適合于求取傾斜地層或各向異性地層的橫波剩余靜校正量的方法。

本發(fā)明的一方面提供了一種求取傾斜地層橫波剩余靜校正量的方法。所述方法采用共接收波點(diǎn)疊加反射層位匹配法求取橫波剩余靜校正量，并包括以下步驟：收集目標(biāo)傾斜地層的轉(zhuǎn)換波地震道集數(shù)據(jù)、轉(zhuǎn)換波疊加速度和切除參數(shù)；保留轉(zhuǎn)換波地震道集中的正偏移距數(shù)據(jù)或負(fù)偏移距數(shù)據(jù)，并利用轉(zhuǎn)換波疊加速度和切除參數(shù)對所保留的正偏移距數(shù)據(jù)或負(fù)偏移距數(shù)據(jù)進(jìn)行疊加處理，以獲得半偏移距疊加剖面；采用所述共接收波點(diǎn)疊加反射層位匹配法處理所述半偏移距疊加剖面，從而求得目標(biāo)傾斜地層橫波剩余靜校正量。

本發(fā)明的另一方面提供了一種求取各向異性地層橫波剩余靜校正量的方法。所述方法采用共接收波點(diǎn)疊加反射層位匹配法求取橫波剩余靜校正量，并包括以下步驟：收集目標(biāo)各向異性地層的轉(zhuǎn)換波地震道集數(shù)據(jù)、轉(zhuǎn)換波疊加速度和切除參數(shù)；保留轉(zhuǎn)換波地震道集中的正偏移距數(shù)據(jù)或負(fù)偏移距數(shù)據(jù)，并利用轉(zhuǎn)換波疊加速度和切除參數(shù)對所保留的正偏移距數(shù)據(jù)或負(fù)偏移距數(shù)據(jù)進(jìn)行疊加處理，以獲得半偏移距疊加剖面；采用所述共接收波點(diǎn)疊加反射層位匹配法處理所述半偏移距疊加剖面，從而求得目標(biāo)各向異性地層橫波剩余靜校正量。

與現(xiàn)有技術(shù)相比，本發(fā)明的有益效果包括：能夠適合于求取傾斜地層或各向異性地層的橫波剩余靜校正量；能夠提高傾斜地層或各向異性地層的橫波剩余靜校正量求取的準(zhǔn)確度。

具體實(shí)施方式

在下文中，將結(jié)合示例性實(shí)施例來詳細(xì)說明本發(fā)明的求取傾斜地層或各向異性地層橫波剩余靜校正量的方法。

總地來講，為了能夠在陡傾地層轉(zhuǎn)換波資料靜校正處理中有效地求取橫波剩余靜校正量，發(fā)明人采用半偏移距疊加方法，以削弱共轉(zhuǎn)換點(diǎn)在傾斜地層界面上分散的效應(yīng)，從而保證陡傾地層(或稱傾斜地層)或各向異性明顯地層能夠在轉(zhuǎn)換波疊加剖面上成像，進(jìn)而可通過常規(guī)的共接收波點(diǎn)疊加反射層位匹配法求得傾斜地層或各向異性明顯地層的橫波剩余靜校正量。

常規(guī)的共接收波點(diǎn)疊加反射層位匹配法求取接收點(diǎn)處大的橫波剩余靜校正量方法要求在共接收波點(diǎn)疊加剖面上拾取對應(yīng)的地震反射層位，但是地下陡傾地層在轉(zhuǎn)換波全偏移距疊加剖面上無法成像，導(dǎo)致該方法無法應(yīng)用。本發(fā)明有效地解決了該問題。

在常用的中間激發(fā)、兩邊接收、多次覆蓋的多波地震觀測系統(tǒng)中，轉(zhuǎn)換波共轉(zhuǎn)換點(diǎn)分布在共接收點(diǎn)兩側(cè)，在地震資料處理中可分為正偏移距道集和負(fù)偏移距道集，兩者的集合統(tǒng)稱為全偏移距道集。常規(guī)的共檢收點(diǎn)疊加剖面就是將所有共檢收點(diǎn)處全偏移距道集疊加，在橫向上可形成疊加剖面。對于各向同性、水平層狀地層，動校正處理后的共檢收點(diǎn)道集中各轉(zhuǎn)換點(diǎn)具有相同的反射時差，具有同相疊加的特點(diǎn)，在疊加剖面上可以看到較為清晰的地層反射，并能夠方便地拾取相應(yīng)的地震反射層位。

然而，發(fā)明人發(fā)現(xiàn)：對于地下地層，普遍存在各向異性且常常發(fā)育有復(fù)雜的傾斜地層，因此，在傾斜地層界面上，共接收點(diǎn)兩側(cè)正、負(fù)偏移距道集的轉(zhuǎn)換波共轉(zhuǎn)換點(diǎn)呈分散狀分布，地層傾斜度越大，各向異性越明顯，共轉(zhuǎn)換點(diǎn)的分散度越高，正、負(fù)偏移距道集數(shù)據(jù)之間的反射時差也越大。如果正偏移距道集的共轉(zhuǎn)換點(diǎn)分布在傾斜地層的上傾方向，則正偏移距道集的共轉(zhuǎn)換點(diǎn)分布在傾斜地層的下傾方向；如果正偏移距道集的共轉(zhuǎn)換點(diǎn)分布在傾斜地層的下傾方向，則正偏移距道集的共轉(zhuǎn)換點(diǎn)分布在傾斜地層的上傾方向。因此，動校正處理后的共接收點(diǎn)道集中各轉(zhuǎn)換點(diǎn)具有不同的反射時差，尤其是正、負(fù)偏移距道集數(shù)據(jù)之間的反射時差相對最大，導(dǎo)致包含正、負(fù)偏移距道集的全偏移距道集不能同相疊加或疊加效果差，在疊加剖面上無法看到較為清晰的傾斜地層反射，因而無法拾取相應(yīng)的地震反射層位。

針對地層傾角和各向異性因素造成的轉(zhuǎn)換波共接收點(diǎn)全偏移距道集中正、負(fù)偏移距數(shù)據(jù)之間存在反射時差，全偏移距疊加不能成像的問題，發(fā)明人采取僅選擇全偏移距數(shù)據(jù)中正或負(fù)偏移距數(shù)據(jù)進(jìn)行疊加，以消除或削弱數(shù)據(jù)間的反射時差，使傾斜地層在轉(zhuǎn)換波疊加剖面上可以成像，為常規(guī)的共接收波點(diǎn)疊加反射層位匹配法求取接收點(diǎn)處大的橫波剩余靜校正量方法能夠在陡傾地層和各向異性明顯的探區(qū)應(yīng)用奠定了基礎(chǔ)。

在本發(fā)明的一個示例性實(shí)施例中，求取傾斜地層橫波剩余靜校正量的方法采用共接收波點(diǎn)疊加反射層位匹配法求取橫波剩余靜校正量。具體來講，求取傾斜地層橫波剩余靜校正量的方法通過以下步驟來實(shí)現(xiàn)：

收集目標(biāo)傾斜地層的轉(zhuǎn)換波地震道集數(shù)據(jù)、轉(zhuǎn)換波疊加速度和切除參數(shù)。這里，轉(zhuǎn)換波地震道集數(shù)據(jù)可以為三維數(shù)據(jù)或二維數(shù)據(jù)。另外，為了方便計算，可利用常規(guī)計算軟件平臺來實(shí)現(xiàn)本步驟。例如，輸入轉(zhuǎn)換波共接收點(diǎn)地震道集、轉(zhuǎn)換波疊加速度、切除參數(shù)文件到常規(guī)的地震數(shù)據(jù)處理軟件中。常規(guī)的地震數(shù)據(jù)處理軟件中都具有切除功能，數(shù)據(jù)處理人員根據(jù)業(yè)內(nèi)專業(yè)知識和經(jīng)驗(yàn)從地震道集處理過程中獲得切除文件。

可僅保留轉(zhuǎn)換波地震道集中的正偏移距數(shù)據(jù)，并利用轉(zhuǎn)換波疊加速度和切除參數(shù)對所保留的正偏移距數(shù)據(jù)進(jìn)行疊加處理，以獲得半偏移距疊加剖面。另外，也可以僅保留轉(zhuǎn)換波地震道集中的負(fù)偏移距數(shù)據(jù)，并利用轉(zhuǎn)換波疊加速度和切除參數(shù)對所保留的負(fù)偏移距數(shù)據(jù)進(jìn)行疊加處理，以獲得半偏移距疊加剖面。常規(guī)的疊加處理是同時保留地震道集中的正偏移距數(shù)據(jù)和負(fù)偏移距數(shù)據(jù)進(jìn)行疊加。

采用常規(guī)的共接收波點(diǎn)疊加反射層位匹配法處理前述半偏移距疊加剖面，從而求得目標(biāo)傾斜地層橫波剩余靜校正量。例如，(1)拾取縱波半偏移距疊加剖面上的橫向連續(xù)性好、地震反射特征明顯的某一地震地質(zhì)反射層位，獲得反射時間T_p；(2)在轉(zhuǎn)換橫波共接收點(diǎn)半偏移距疊加剖面上拾取同一地震地質(zhì)反射層位，獲得反射時間T_ps(由于橫波速度小于縱波速度，導(dǎo)致同一地質(zhì)層位在縱波和轉(zhuǎn)換橫波共接收點(diǎn)疊加剖面上具有不同的地震反射時間，因此，反射層位匹配是指必須對比解釋縱波和轉(zhuǎn)換橫波共接收點(diǎn)半偏移距疊加剖面，才能準(zhǔn)確拾取轉(zhuǎn)換橫波共接收點(diǎn)疊加剖面上對應(yīng)的層位)；(3)根據(jù)以下公式計算縱橫波速度比r₀：

r₀＝(2T_ps-1)/T_p

其中，r₀、T_p、T_ps分別是縱橫波速度比、拾取層位的縱波反射時間、拾取層位的轉(zhuǎn)換橫波反射時間；(4)根據(jù)已知的縱波靜校正量和縱橫波速度比求取目標(biāo)傾斜地層橫波剩余靜校正量，例如，求取接收點(diǎn)處大的橫波剩余靜校正量。

在本發(fā)明的另一個示例性實(shí)施例中，求取各向異性地層橫波剩余靜校正量的方法采用共接收波點(diǎn)疊加反射層位匹配法求取橫波剩余靜校正量。

具體來講，求取各向異性地層橫波剩余靜校正量的方法通過以下步驟來實(shí)現(xiàn)：

收集目標(biāo)各向異性地層的轉(zhuǎn)換波地震道集數(shù)據(jù)、轉(zhuǎn)換波疊加速度和切除參數(shù)。這里，轉(zhuǎn)換波地震道集數(shù)據(jù)可以為三維數(shù)據(jù)或二維數(shù)據(jù)。另外，為了方便計算，可利用常規(guī)計算軟件平臺來實(shí)現(xiàn)本步驟。例如，輸入轉(zhuǎn)換波共接收點(diǎn)地震道集、轉(zhuǎn)換波疊加速度、切除參數(shù)文件到常規(guī)的地震數(shù)據(jù)處理軟件中。常規(guī)的地震數(shù)據(jù)處理軟件中都具有切除功能，數(shù)據(jù)處理人員根據(jù)業(yè)內(nèi)專業(yè)知識和經(jīng)驗(yàn)從地震道集處理過程中獲得切除文件。

僅保留轉(zhuǎn)換波地震道集中的正偏移距數(shù)據(jù)，并利用轉(zhuǎn)換波疊加速度和切除參數(shù)對所保留的正偏移距數(shù)據(jù)進(jìn)行疊加處理，以獲得半偏移距疊加剖面。常規(guī)的疊加處理是同時保留地震道集中的正偏移距數(shù)據(jù)和負(fù)偏移距數(shù)據(jù)進(jìn)行疊加。另外，也可以僅保留轉(zhuǎn)換波地震道集中的負(fù)偏移距數(shù)據(jù)，并利用轉(zhuǎn)換波疊加速度和切除參數(shù)對所保留的負(fù)偏移距數(shù)據(jù)進(jìn)行疊加處理，以獲得半偏移距疊加剖面。

采用常規(guī)的共接收波點(diǎn)疊加反射層位匹配法處理前述半偏移距疊加剖面，從而求得目標(biāo)各向異性地層橫波剩余靜校正量。例如，(1)拾取縱波半偏移距疊加剖面上的橫向連續(xù)性好、地震反射特征明顯的某一地震地質(zhì)反射層位，獲得反射時間T_p；(2)在轉(zhuǎn)換橫波共接收點(diǎn)半偏移距疊加剖面上拾取同一地震地質(zhì)反射層位，獲得反射時間T_ps(由于橫波速度小于縱波速度，導(dǎo)致同一地質(zhì)層位在縱波和轉(zhuǎn)換橫波共接收點(diǎn)疊加剖面上具有不同的地震反射時間，因此，反射層位匹配是指必須對比解釋縱波和轉(zhuǎn)換橫波共接收點(diǎn)半偏移距疊加剖面，才能準(zhǔn)確拾取轉(zhuǎn)換橫波共接收點(diǎn)疊加剖面上對應(yīng)的層位)；(3)根據(jù)以下公式計算縱橫波速度比r₀：

r₀＝(2T_ps-1)/T_p

其中，r₀、T_p、T_ps分別是縱橫波速度比、拾取層位的縱波反射時間、拾取層位的轉(zhuǎn)換橫波反射時間；(4)根據(jù)已知的縱波靜校正量和縱橫波速度比求取目標(biāo)各向異性地層橫波剩余靜校正量，例如，求取接收點(diǎn)處大的橫波剩余靜校正量。

綜上所述，本發(fā)明的方法針對因地層傾角或各向異性因素造成的轉(zhuǎn)換波共接收點(diǎn)全偏移距道集中正、負(fù)偏移距數(shù)據(jù)之間存在反射時差，全偏移距疊加不能成像的問題，通過采取僅選擇全偏移距數(shù)據(jù)中的正偏移距數(shù)據(jù)或負(fù)偏移距數(shù)據(jù)進(jìn)行疊加，從而消除或削弱了數(shù)據(jù)間的反射時差，使傾斜地層在轉(zhuǎn)換波疊加剖面上可以成像，進(jìn)而通過常規(guī)的共接收波點(diǎn)疊加反射層位匹配法求取接收點(diǎn)處大的橫波剩余靜校正量。本發(fā)明的方法尤其適合于在陡傾地層和各向異性明顯的探區(qū)應(yīng)用。本發(fā)明的方法能應(yīng)用于二維和三維轉(zhuǎn)換波多波地震資料處理，對油氣、頁巖氣、礦產(chǎn)資源勘探有重要的應(yīng)用價值。

盡管上面已經(jīng)結(jié)合示例性實(shí)施例描述了本發(fā)明，但是本領(lǐng)域普通技術(shù)人員應(yīng)該清楚，在不脫離權(quán)利要求的精神和范圍的情況下，可以對上述實(shí)施例進(jìn)行各種修改。