本申請(qǐng)涉及油氣勘探技術(shù)領(lǐng)域，特別涉及一種深度域地震子波的確定方法和裝置。

背景技術(shù)：

在油氣勘探中，地震子波，尤其是深度域地震子波是一種重要的地震數(shù)據(jù)。通常，地震子波指的是一段具有確定的起始時(shí)間，能量有限，且具有一定延續(xù)長(zhǎng)度的信號(hào)。它是地震記錄中的基本單元。一般認(rèn)為，地震震源激發(fā)時(shí)所產(chǎn)生的地震波僅是一個(gè)延續(xù)時(shí)間極短的尖脈沖，隨著尖脈沖在粘彈性介質(zhì)中傳播，尖脈沖的高頻成分會(huì)很快衰減，波形會(huì)隨之增長(zhǎng)，進(jìn)而形成了地震子波，然后以地震子波的形式在地下傳播。由于地震子波的上述特點(diǎn)，在具體實(shí)施時(shí)，經(jīng)常需要使用深度域地震子波。例如，在正演問(wèn)題的分析中，需要通過(guò)波動(dòng)方程或褶積模型結(jié)合深度域地震子波產(chǎn)生模擬數(shù)據(jù)，進(jìn)而可以得到合成記錄，以合成記錄作為指導(dǎo)，進(jìn)行具體的油氣勘探或者儲(chǔ)層預(yù)測(cè)。因此，如何準(zhǔn)確地確定深度域的地震子波一直是人們關(guān)注的一個(gè)問(wèn)題。

現(xiàn)有的深度域地震子波的確定方法通常不能直接地提取得到深度域地震子波，常常是先獲取時(shí)間域地震子波，再根據(jù)褶積定理將時(shí)間域地震子波轉(zhuǎn)化為深度域地震子波。但是，具體實(shí)施時(shí)，由于深度域地震子波與反射系數(shù)的褶積并不滿(mǎn)足褶積定理，即不能夠滿(mǎn)足子波隨深度遞變的特性。因此，現(xiàn)有的深度域地震子波的確定方法具體實(shí)施時(shí)，往往存在無(wú)法直接確定深度域地震子波、所確定的深度域地震子波準(zhǔn)確度差、確定過(guò)程復(fù)雜的技術(shù)問(wèn)題。

針對(duì)上述問(wèn)題，目前尚未提出有效的解決方案。

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

本申請(qǐng)實(shí)施方式提供了一種深度域地震子波的確定方法和裝置，以解決現(xiàn)有的深度域地震子波確定方法中存在的無(wú)法直接確定深度域地震子波、所確定的深度域地震子波準(zhǔn)確度差、確定過(guò)程復(fù)雜的技術(shù)問(wèn)題。

本申請(qǐng)實(shí)施方式提供了一種深度域地震子波的確定方法，包括：

獲取待測(cè)區(qū)域的聲波測(cè)井?dāng)?shù)據(jù)；

根據(jù)所述聲波測(cè)井?dāng)?shù)據(jù)，確定深度域反射系數(shù)序列和井位置的深度域成像數(shù)據(jù)；

根據(jù)所述深度域反射系數(shù)序列和所述井位置的深度域成像數(shù)據(jù)，確定深度域地震子波。

在一個(gè)實(shí)施方式中，根據(jù)聲波測(cè)井?dāng)?shù)據(jù)，確定深度域反射系數(shù)序列，包括：

根據(jù)所述聲波測(cè)井?dāng)?shù)據(jù)，獲取速度場(chǎng)數(shù)據(jù)和密度場(chǎng)數(shù)據(jù)；

根據(jù)所述速度場(chǎng)數(shù)據(jù)和所述密度場(chǎng)數(shù)據(jù)，確定所述深度域反射系數(shù)序列。

在一個(gè)實(shí)施方式中，根據(jù)所述速度場(chǎng)數(shù)據(jù)和所述密度場(chǎng)數(shù)據(jù)，確定所述深度域反射系數(shù)序列，包括：按照以下公式確定所述深度域反射系數(shù)序列：

其中，r^d(zi)為深度域反射系數(shù)序列中的第i個(gè)采樣點(diǎn)的深度域反射系數(shù)，υ(zi)為速度場(chǎng)數(shù)據(jù)中的第i個(gè)采樣點(diǎn)的速度場(chǎng)數(shù)據(jù)，ρ(zi)為密度場(chǎng)數(shù)據(jù)中的第i個(gè)采樣點(diǎn)的密度場(chǎng)數(shù)據(jù)，zi為測(cè)井?dāng)?shù)據(jù)中的第i個(gè)采樣點(diǎn)的深度值，i為測(cè)井?dāng)?shù)據(jù)中深度方向上的采樣點(diǎn)的編號(hào)。

在一個(gè)實(shí)施方式中，根據(jù)所述深度域反射系數(shù)序列和所述井位置的深度域成像數(shù)據(jù)，確定深度域地震子波，包括：

根據(jù)所述深度域反射系數(shù)序列，確定深度域反射系數(shù)序列的自相關(guān)參數(shù)；

根據(jù)所述深度域反射系數(shù)序列和所述井位置的深度域成像數(shù)據(jù)，確定深度域反射系數(shù)序列和深度域成像數(shù)據(jù)的互相關(guān)參數(shù)；

根據(jù)所述自相關(guān)參數(shù)和所述互相關(guān)參數(shù)，確定匹配算子，并將所述匹配算子作為所述深度域地震子波。

在一個(gè)實(shí)施方式中，根據(jù)所述深度域反射系數(shù)序列，確定深度域反射系數(shù)序列的自相關(guān)參數(shù)，包括：按照以下公式確定所述自相關(guān)參數(shù)：

其中，為第i個(gè)采樣點(diǎn)的自相關(guān)參數(shù)，r^d(n)為深度為n時(shí)的深度域反射系數(shù)，r^d(n+zi)為深度為第i個(gè)采樣點(diǎn)的深度和n之和時(shí)的深度域反射系數(shù)，n為深度變量，取值為1至l范圍內(nèi)的自然數(shù)，zi為測(cè)井?dāng)?shù)據(jù)中的第i個(gè)采樣點(diǎn)的深度值，l為算子的長(zhǎng)度。

在一個(gè)實(shí)施方式中，根據(jù)所述深度域反射系數(shù)序列和所述井位置的深度域成像數(shù)據(jù)，確定深度域反射系數(shù)序列和深度域成像數(shù)據(jù)的互相關(guān)參數(shù)，包括：按照以下公式確定所述互相關(guān)參數(shù)：

其中，μrx(zi)為第i個(gè)采樣點(diǎn)的互相關(guān)參數(shù)，x^d(n)為深度為n時(shí)的成像數(shù)據(jù)，r^d(n+z)為深度為第i個(gè)采樣點(diǎn)的深度和n之和時(shí)的深度域反射系數(shù)，n為深度變量，取值為1至l范圍內(nèi)的自然數(shù)，zi為測(cè)井?dāng)?shù)據(jù)中的第i個(gè)采樣點(diǎn)的深度值，l為算子的長(zhǎng)度。

在一個(gè)實(shí)施方式中，根據(jù)所述自相關(guān)參數(shù)和所述互相關(guān)參數(shù)，確定匹配算子，包括：

根據(jù)所述自相關(guān)參數(shù)和所述互相關(guān)參數(shù)，構(gòu)建求解矩陣方程；

通過(guò)求解所述求解矩陣方程，確定所述匹配算子。

在一個(gè)實(shí)施方式中，根據(jù)所述自相關(guān)參數(shù)和所述互相關(guān)參數(shù)，構(gòu)建求解矩陣方程，包括：

根據(jù)所述自相關(guān)參數(shù)和所述互相關(guān)參數(shù)，構(gòu)建以下方程作為所述求解矩陣方程：

其中，為深度為z時(shí)的自相關(guān)參數(shù)，μrx(z)為深度為z時(shí)的互相關(guān)參數(shù)，w^d(z)為深度為z時(shí)的匹配算子，l為匹配算子的長(zhǎng)度，z為深度，z的取值為0至l范圍內(nèi)的自然數(shù)。

在一個(gè)實(shí)施方式中，在確定所述深度域地震子波后，所述方法還包括：

根據(jù)所述深度域地震子波，對(duì)疊前深度偏移數(shù)據(jù)進(jìn)行反演處理，得到處理結(jié)果；

根據(jù)所述處理結(jié)果，進(jìn)行油氣勘探。

基于相同的發(fā)明構(gòu)思，本申請(qǐng)實(shí)施方式還提供了一種深度域地震子波的確定裝置，包括：

獲取模塊，用于獲取待測(cè)區(qū)域的聲波測(cè)井?dāng)?shù)據(jù)；

第一確定模塊，用于根據(jù)所述聲波測(cè)井?dāng)?shù)據(jù)，確定深度域反射系數(shù)序列和井位置的深度域成像數(shù)據(jù)；

第二確定模塊，用于根據(jù)所述深度域反射系數(shù)序列和所述井位置的深度域成像數(shù)據(jù)，確定深度域地震子波。

在一個(gè)實(shí)施方式中，所述第一確定模塊包括：

獲取單元，用于根據(jù)所述聲波測(cè)井?dāng)?shù)據(jù)，獲取速度場(chǎng)數(shù)據(jù)和密度場(chǎng)數(shù)據(jù)；

第一確定單元，用于根據(jù)所述速度場(chǎng)數(shù)據(jù)和所述密度場(chǎng)數(shù)據(jù)，確定所述深度域反射系數(shù)序列和所述井位置的深度域成像數(shù)據(jù)。

在一個(gè)實(shí)施方式中，所述第二確定模塊包括：

第二確定單元，用于根據(jù)所述深度域反射系數(shù)序列，確定深度域反射系數(shù)序列的自相關(guān)參數(shù)；

第三確定單元，用于根據(jù)所述深度域反射系數(shù)序列和所述井位置的深度域成像數(shù)據(jù)，確定深度域反射系數(shù)序列和深度域成像數(shù)據(jù)的互相關(guān)參數(shù)；

第四確定單元，用于根據(jù)所述自相關(guān)參數(shù)和所述互相關(guān)參數(shù)，確定匹配算子，并將所述匹配算子作為所述深度域地震子波。

在本申請(qǐng)實(shí)施方式中，通過(guò)聲波測(cè)井?dāng)?shù)據(jù)，分別確定深度域反射系數(shù)序列和井位置的深度域成像數(shù)據(jù)；進(jìn)而可以根據(jù)深度域反射系數(shù)序列和井位置的深度域成像數(shù)據(jù)，確定深度域地震子波。從而解決了現(xiàn)有的深度域地震子波的確定方法中存在的無(wú)法直接求取深度域地震子波、確定過(guò)程復(fù)雜、所確定的深度域地震子波準(zhǔn)確度差的技術(shù)問(wèn)題，達(dá)到了準(zhǔn)確、簡(jiǎn)便地確定深度域地震子波的技術(shù)效果。

附圖說(shuō)明

為了更清楚地說(shuō)明本申請(qǐng)實(shí)施例或現(xiàn)有技術(shù)中的技術(shù)方案，下面將對(duì)實(shí)施例或現(xiàn)有技術(shù)描述中所需要使用的附圖作簡(jiǎn)單地介紹，顯而易見(jiàn)地，下面描述中的附圖僅僅是本申請(qǐng)中記載的一些實(shí)施例，對(duì)于本領(lǐng)域普通技術(shù)人員來(lái)講，在不付出創(chuàng)造性勞動(dòng)性的前提下，還可以根據(jù)這些附圖獲得其他的附圖。

圖1是根據(jù)本申請(qǐng)實(shí)施方式的深度域地震子波的確定方法的處理流程圖；

圖2是根據(jù)本申請(qǐng)實(shí)施方式的深度域地震子波的確定裝置的組成結(jié)構(gòu)圖；

圖3是在一個(gè)場(chǎng)景示例中獲取的井?dāng)?shù)據(jù)得到的深度域反射系數(shù)示意圖；

圖4是在一個(gè)場(chǎng)景示例中獲取的時(shí)間域成像示意圖；

圖5是在一個(gè)場(chǎng)景示例中獲取的深度域子波示意圖；

圖6是在一個(gè)場(chǎng)景示例中獲取的深度域等效子波合成記錄示意圖；

圖7是在一個(gè)場(chǎng)景示例中獲取的深度域等效子波合成記錄與時(shí)深轉(zhuǎn)換記錄的對(duì)比示意圖。

具體實(shí)施方式

為了使本技術(shù)領(lǐng)域的人員更好地理解本申請(qǐng)中的技術(shù)方案，下面將結(jié)合本申請(qǐng)實(shí)施例中的附圖，對(duì)本申請(qǐng)實(shí)施例中的技術(shù)方案進(jìn)行清楚、完整地描述，顯然，所描述的實(shí)施例僅僅是本申請(qǐng)一部分實(shí)施例，而不是全部的實(shí)施例?；诒旧暾?qǐng)中的實(shí)施例，本領(lǐng)域普通技術(shù)人員在沒(méi)有做出創(chuàng)造性勞動(dòng)前提下所獲得的所有其他實(shí)施例，都應(yīng)當(dāng)屬于本申請(qǐng)保護(hù)的范圍。

考慮到現(xiàn)有的深度域地震子波的確定方法通常是先獲取時(shí)間域的地震子波，再將時(shí)間域的地震子波轉(zhuǎn)化為深度域的地震子波。但是，由于子波與反射系數(shù)的褶積并不滿(mǎn)足褶積定理。因此，這種將時(shí)間域地震子波轉(zhuǎn)換成深度域地震子波的方法具體實(shí)施時(shí)，往往會(huì)存在無(wú)法直接求取深度域地震子波、確定過(guò)程復(fù)雜、所確定的深度域地震子波準(zhǔn)確度差的技術(shù)問(wèn)題。針對(duì)產(chǎn)生上述技術(shù)問(wèn)題的根本原因，本申請(qǐng)?zhí)岢隹梢愿鶕?jù)測(cè)井?dāng)?shù)據(jù)先獲取深度域反射系數(shù)序列和深度域成像數(shù)據(jù)，再通過(guò)根據(jù)深度域反射系數(shù)序列和深度域成像數(shù)據(jù)，求取匹配算子作為深度域地震子波。即可以直接根據(jù)深度域地震數(shù)據(jù)確定深度域地震子波。從而解決了現(xiàn)有的深度域地震子波的確定方法中存在的無(wú)法直接求取深度域地震子波、確定過(guò)程復(fù)雜、所確定的深度域地震子波準(zhǔn)確度差的技術(shù)問(wèn)題，達(dá)到準(zhǔn)確、簡(jiǎn)便地確定深度域地震子波的技術(shù)效果。

基于上述思考思路，本申請(qǐng)實(shí)施方式提供了一種深度域地震子波的確定方法。請(qǐng)參閱圖1的根據(jù)本申請(qǐng)實(shí)施方式的深度域地震子波的確定方法的處理流程圖。本申請(qǐng)實(shí)施方式提供的深度域地震子波的確定方法，具體可以包括以下步驟。

步驟101：獲取待測(cè)區(qū)域的聲波測(cè)井?dāng)?shù)據(jù)。

在本實(shí)施方式中，聲波測(cè)井?dāng)?shù)據(jù)具體可以是測(cè)井?dāng)?shù)據(jù)中的一種。具體的，聲波測(cè)井?dāng)?shù)據(jù)又可以包括：速度場(chǎng)數(shù)據(jù)和密度場(chǎng)數(shù)據(jù)等。

步驟102：根據(jù)所述聲波測(cè)井?dāng)?shù)據(jù)，確定深度域反射系數(shù)序列和井位置的深度域成像數(shù)據(jù)。

在一個(gè)實(shí)施方式中，為了能夠分別確定深度域反射系數(shù)序列和井位置的深度域成像數(shù)據(jù)，具體實(shí)施時(shí)，可以按照以下步驟執(zhí)行。

s102-1：根據(jù)聲波測(cè)井?dāng)?shù)據(jù)，確定深度域反射系數(shù)序列。

在一個(gè)實(shí)施方式中，為了獲取上述深度域反射系數(shù)序列，具體可以按照以下方式通過(guò)分別確定對(duì)個(gè)深度域反射系數(shù)，進(jìn)而確定上述深度域反射系數(shù)序列。

s102-1-1：根據(jù)所述聲波測(cè)井?dāng)?shù)據(jù)，獲取速度場(chǎng)數(shù)據(jù)和密度場(chǎng)數(shù)據(jù)；

s102-1-2：根據(jù)所述速度場(chǎng)數(shù)據(jù)和所述密度場(chǎng)數(shù)據(jù)，確定所述深度域反射系數(shù)序列。

在一個(gè)實(shí)施方式中，具體實(shí)施時(shí)，可以按照以下公式根據(jù)速度場(chǎng)數(shù)據(jù)和密度場(chǎng)數(shù)據(jù)分別確定各個(gè)深度域反射系數(shù)，進(jìn)而確定上述深度域反射系數(shù)序列。

其中，r^d(zi)為深度域反射系數(shù)序列中的第i個(gè)采樣點(diǎn)的深度域反射系數(shù)，υ(zi)為速度場(chǎng)數(shù)據(jù)中的第i個(gè)采樣點(diǎn)的速度場(chǎng)數(shù)據(jù)，ρ(zi)為密度場(chǎng)數(shù)據(jù)中的第i個(gè)采樣點(diǎn)的密度場(chǎng)數(shù)據(jù)，zi為測(cè)井?dāng)?shù)據(jù)中的第i個(gè)采樣點(diǎn)的深度值，i為測(cè)井?dāng)?shù)據(jù)中深度方向上的采樣點(diǎn)的編號(hào)。

在本實(shí)施方式中，按照上述公式計(jì)算得到各個(gè)深度域反射系數(shù)后，將得到的多個(gè)深度域反射系數(shù)組合成深度域反射系數(shù)序列。其中，上述深度域反射系數(shù)序列具體可以用以下方式表示：r^d(z)，z＝z1,z2...zi...zn，其中，zi為對(duì)應(yīng)的深度。

在本實(shí)施方式中，采樣點(diǎn)的個(gè)數(shù)可以根據(jù)具體情況和施工要求，靈活設(shè)定。例如，測(cè)井的深度方向?yàn)?000米，可以每4米作為一個(gè)采樣間隔。如此，一共有750個(gè)采樣點(diǎn)，也就是i具體可以為1到750中的自然數(shù)。

s102-2：根據(jù)聲波測(cè)井?dāng)?shù)據(jù)，確定井位置的深度域成像數(shù)據(jù)。

在一個(gè)實(shí)施方式中，具體實(shí)施時(shí)，可以根據(jù)測(cè)井?dāng)?shù)據(jù)確定井位置的深度域成像數(shù)據(jù)。其中，上述井指的可以是具體施工時(shí)現(xiàn)場(chǎng)設(shè)置好井。相應(yīng)的，井位置是指具體施工時(shí)現(xiàn)場(chǎng)設(shè)置好井后就確定的井的位置。上述井位置的深度域成像數(shù)據(jù)具體可以用以下方式表示：x^d(z)，z＝z1,z2...zi...zn，其中，zi為對(duì)應(yīng)的深度。

步驟103：根據(jù)所述深度域反射系數(shù)序列和所述井位置的深度域成像數(shù)據(jù)，確定深度域地震子波。

在一個(gè)實(shí)施方式中，為了避免先確定時(shí)間域地震子波，再將時(shí)間域地震子波轉(zhuǎn)化為深度域地震子波所導(dǎo)致的誤差，本申請(qǐng)實(shí)施方式提出可以根據(jù)所述深度域反射系數(shù)序列和所述井位置的深度域成像數(shù)據(jù)，直接確定深度域地震子波。具體實(shí)施時(shí)，可以按照以下步驟執(zhí)行。

s103-1：根據(jù)所述深度域反射系數(shù)序列，確定深度域反射系數(shù)序列的自相關(guān)參數(shù)。

在一個(gè)實(shí)施方式中，為了確定深度域反射系數(shù)序列的自相關(guān)參數(shù)，具體實(shí)施時(shí)，可以按照以下公式計(jì)算得到上述的反射系數(shù)序列的自相參數(shù)：

其中，為第i個(gè)采樣點(diǎn)的自相關(guān)參數(shù)，r^d(n)為深度為n時(shí)的深度域反射系數(shù)，r^d(n+zi)為深度為第i個(gè)采樣點(diǎn)的深度和n之和時(shí)的深度域反射系數(shù)，n為深度變量，取值為1至l范圍內(nèi)的自然數(shù)(包括1和l)，zi為測(cè)井?dāng)?shù)據(jù)中的第i個(gè)采樣點(diǎn)的深度值，l為算子的長(zhǎng)度。

s103-2：根據(jù)所述深度域反射系數(shù)序列和所述井位置的深度域成像數(shù)據(jù)，確定深度域反射系數(shù)序列和深度域成像數(shù)據(jù)的互相關(guān)參數(shù)。

在一個(gè)實(shí)施方式中，為了確定深度域反射系數(shù)序列和深度域成像數(shù)據(jù)的互相關(guān)參數(shù)，具體實(shí)施時(shí)，可以按照以下公式計(jì)算得到上述深度域反射系數(shù)序列和深度域成像數(shù)據(jù)的互相關(guān)參數(shù)：

其中，μrx(zi)為第i個(gè)采樣點(diǎn)的互相關(guān)參數(shù)，x^d(n)為深度為n時(shí)的成像數(shù)據(jù)，r^d(n+z)為深度為第i個(gè)采樣點(diǎn)的深度和n之和時(shí)的深度域反射系數(shù)，n為深度變量，取值為1至l范圍內(nèi)的自然數(shù)(包括1和l)，zi為測(cè)井?dāng)?shù)據(jù)中的第i個(gè)采樣點(diǎn)的深度值，l為算子的長(zhǎng)度。

s103-3：根據(jù)所述自相關(guān)參數(shù)和所述互相關(guān)參數(shù)，確定匹配算子，并將所述匹配算子作為所述深度域地震子波。

在一個(gè)實(shí)施方式中，為了可以直接確定深度域地震子波，具體實(shí)施時(shí)可以根據(jù)計(jì)算出的深度域反射系數(shù)序列的自相關(guān)參數(shù)、計(jì)算出的深度域反射系數(shù)序列和深度域成像數(shù)據(jù)的互相關(guān)參數(shù)，求取出以深度域成像數(shù)據(jù)為期望輸出，且基于井?dāng)?shù)據(jù)的深度域反射系數(shù)序列的最佳匹配算子。即直接根據(jù)深度域數(shù)據(jù)確定深度域地震子波。其中，該算子可以認(rèn)為是最小二乘意義下的深度域等效子波。因此，可以將該算子作為上述深度域地震子波。需要說(shuō)明的是，本申請(qǐng)實(shí)施方式提出的確定方法是基于最佳匹配的算法，使子波與反射系數(shù)褶積得到的合成記錄與深度域地震道在最佳匹配情況下得到的子波作為最佳深度域地震子波。如此，使用所得到的地震域子波可以更好的適用于后期合成記錄制作合反演，即使用上述深度域地震子波可以得到更為準(zhǔn)確的處理結(jié)果。

在一個(gè)實(shí)施方式中，為了確定上述匹配算子，具體實(shí)施時(shí)，可以按照以下步驟執(zhí)行：

s103-3-1：根據(jù)所述自相關(guān)參數(shù)和所述互相關(guān)參數(shù)，構(gòu)建求解矩陣方程；

s103-3-2：通過(guò)求解所述求解矩陣方程，確定所述匹配算子。

在一個(gè)實(shí)施方式中，根據(jù)所述自相關(guān)參數(shù)和所述互相關(guān)參數(shù)，構(gòu)建求解矩陣方程，具體可以包括以下內(nèi)容：

根據(jù)所述自相關(guān)參數(shù)和所述互相關(guān)參數(shù)，構(gòu)建以下方程作為所述求解矩陣方程：

其中，為深度為z時(shí)的自相關(guān)參數(shù)，μrx(z)為深度為z時(shí)的互相關(guān)參數(shù)，w^d(z)為深度為z時(shí)的匹配算子，l為匹配算子的長(zhǎng)度(即深度域地震子波波長(zhǎng))，z為深度，z的取值為0至l范圍內(nèi)的自然數(shù)。

在一個(gè)實(shí)施方式中，在構(gòu)建出了求解矩陣方程后，為了提高求解速率，達(dá)到快速求解上述方程以確定對(duì)應(yīng)匹配算子的目的。具體實(shí)施時(shí)，可以通過(guò)遞推算法等計(jì)算方法對(duì)上述求解矩陣方程進(jìn)行求解。具體的，例如，可以通過(guò)萊文森(levinson)快速遞推算法求解上述求解矩陣方程，進(jìn)而獲得上述匹配算子作為上述深度域地震子波。當(dāng)然，需要說(shuō)明的是，具體實(shí)施時(shí)，還可以根據(jù)具體情況和實(shí)施條件，選擇除上述所列舉的萊文森算法以外的算法對(duì)上述方程進(jìn)行求解。對(duì)此，本申請(qǐng)不作限定。

在一個(gè)實(shí)施方式中，在確定出上述深度域地震子波后，為了可以進(jìn)一步地進(jìn)行油氣勘探。具體實(shí)施時(shí)，上述方法還可以包括以下步驟。

s104-1：根據(jù)所述深度域地震子波，對(duì)疊前深度偏移數(shù)據(jù)進(jìn)行反演處理，得到處理結(jié)果；

s104-2：根據(jù)所述處理結(jié)果，進(jìn)行油氣勘探。

在本實(shí)施方式中，在根據(jù)深度域地震子波，對(duì)疊前深度偏移數(shù)據(jù)進(jìn)行反演得到處理結(jié)果后，除了可以利用上述處理結(jié)果作為指導(dǎo)依據(jù)進(jìn)行油氣勘探，還可以以上述處理結(jié)果作為參考依據(jù)，進(jìn)行儲(chǔ)層預(yù)測(cè)等具體施工。

在本申請(qǐng)實(shí)施例中，相較于現(xiàn)有方法，通過(guò)根據(jù)聲波測(cè)井?dāng)?shù)據(jù)，分別確定深度域反射系數(shù)序列和井位置的深度域成像數(shù)據(jù)；進(jìn)而可以根據(jù)深度域反射系數(shù)序列和井位置的深度域成像數(shù)據(jù)，確定深度域地震子波。即直接根據(jù)深度域地震數(shù)據(jù)，確定深度域地震子波。從而解決了現(xiàn)有的深度域地震子波的確定方法中存在的無(wú)法直接求取深度域地震子波、確定過(guò)程復(fù)雜、所確定的深度域地震子波準(zhǔn)確度差的技術(shù)問(wèn)題，達(dá)到了準(zhǔn)確、簡(jiǎn)便地確定深度域地震子波的技術(shù)效果。

基于同一發(fā)明構(gòu)思，本發(fā)明實(shí)施方式中還提供了一種深度域地震子波的確定裝置，如下面的實(shí)施方式所述。由于裝置解決問(wèn)題的原理與深度域地震子波的確定方法相似，因此裝置的實(shí)施可以參見(jiàn)深度域地震子波的確定方法的實(shí)施，重復(fù)之處不再贅述。以下所使用的，術(shù)語(yǔ)“單元”或者“模塊”可以實(shí)現(xiàn)預(yù)定功能的軟件和/或硬件的組合。盡管以下實(shí)施例所描述的裝置較佳地以軟件來(lái)實(shí)現(xiàn)，但是硬件，或者軟件和硬件的組合的實(shí)現(xiàn)也是可能并被構(gòu)想的。請(qǐng)參閱圖2，是本申請(qǐng)實(shí)施方式的深度域地震子波的確定裝置的一種組成結(jié)構(gòu)圖，該裝置具體可以包括：獲取模塊201、第一確定模塊202、第二確定模塊203，下面對(duì)該結(jié)構(gòu)進(jìn)行具體說(shuō)明。

獲取模塊201，具體可以用于獲取待測(cè)區(qū)域的聲波測(cè)井?dāng)?shù)據(jù)；

第一確定模塊202，具體可以用于根據(jù)所述聲波測(cè)井?dāng)?shù)據(jù)，確定深度域反射系數(shù)序列和井位置的深度域成像數(shù)據(jù)；

第二確定模塊203，具體可以用于根據(jù)所述深度域反射系數(shù)序列和所述井位置的深度域成像數(shù)據(jù)，確定深度域地震子波。

在一個(gè)實(shí)施方式中，為了確定上述深度域反射系數(shù)序列，具體實(shí)施時(shí)，所述第一確定模塊202具體可以包括以下幾個(gè)單元。

獲取單元，具體可以用于根據(jù)所述聲波測(cè)井?dāng)?shù)據(jù)，獲取速度場(chǎng)數(shù)據(jù)和密度場(chǎng)數(shù)據(jù)；

第一確定單元，具體可以用于根據(jù)所述速度場(chǎng)數(shù)據(jù)和所述密度場(chǎng)數(shù)據(jù)，確定所述深度域反射系數(shù)序列和所述井位置的深度域成像數(shù)據(jù)。

在一個(gè)實(shí)施方式中，為了能夠根據(jù)深度域地震數(shù)據(jù)確定深度域地震子波，具體實(shí)施時(shí)，所述第二確定模塊203具體可以包括以下幾個(gè)單元。

第二確定單元，具體可以用于根據(jù)所述深度域反射系數(shù)序列，確定深度域反射系數(shù)序列的自相關(guān)參數(shù)；

第三確定單元，具體可以用于根據(jù)所述深度域反射系數(shù)序列和所述井位置的深度域成像數(shù)據(jù)，確定深度域反射系數(shù)序列和深度域成像數(shù)據(jù)的互相關(guān)參數(shù)；

第四確定單元，具體可以用于根據(jù)所述自相關(guān)參數(shù)和所述互相關(guān)參數(shù)，確定匹配算子，并將所述匹配算子作為所述深度域地震子波。

本說(shuō)明書(shū)中的各個(gè)實(shí)施例均采用遞進(jìn)的方式描述，各個(gè)實(shí)施例之間相同相似的部分互相參見(jiàn)即可，每個(gè)實(shí)施例重點(diǎn)說(shuō)明的都是與其他實(shí)施例的不同之處。尤其，對(duì)于系統(tǒng)實(shí)施例而言，由于其基本相似于方法實(shí)施例，所以描述的比較簡(jiǎn)單，相關(guān)之處參見(jiàn)方法實(shí)施例的部分說(shuō)明即可。

需要說(shuō)明的是，上述實(shí)施方式闡明的系統(tǒng)、裝置、模塊或單元，具體可以由計(jì)算機(jī)芯片或?qū)嶓w實(shí)現(xiàn)，或者由具有某種功能的產(chǎn)品來(lái)實(shí)現(xiàn)。為了描述的方便，在本說(shuō)明書(shū)中，描述以上裝置時(shí)以功能分為各種單元分別描述。當(dāng)然，在實(shí)施本申請(qǐng)時(shí)可以把各單元的功能在同一個(gè)或多個(gè)軟件和/或硬件中實(shí)現(xiàn)。

此外，在本說(shuō)明書(shū)中，諸如第一和第二這樣的形容詞僅可以用于將一個(gè)元素或動(dòng)作與另一元素或動(dòng)作進(jìn)行區(qū)分，而不必要求或暗示任何實(shí)際的這種關(guān)系或順序。在環(huán)境允許的情況下，參照元素或部件或步驟(等)不應(yīng)解釋為局限于僅元素、部件、或步驟中的一個(gè)，而可以是元素、部件、或步驟中的一個(gè)或多個(gè)等。

從以上的描述中，可以看出，本申請(qǐng)實(shí)施方式提供的深度域地震子波的確定方法和裝置通過(guò)根據(jù)聲波測(cè)井?dāng)?shù)據(jù)，分別確定深度域反射系數(shù)序列和井位置的深度域成像數(shù)據(jù)；進(jìn)而可以根據(jù)深度域反射系數(shù)序列和井位置的深度域成像數(shù)據(jù)，確定深度域地震子波。從而解決了現(xiàn)有的深度域地震子波的確定方法中存在的無(wú)法直接求取深度域地震子波、確定過(guò)程復(fù)雜、所確定的深度域地震子波準(zhǔn)確度差的技術(shù)問(wèn)題，達(dá)到了準(zhǔn)確、簡(jiǎn)便地確定深度域地震子波的技術(shù)效果。又通過(guò)利用萊文森遞推算法求解匹配算子，提高了求解速率，達(dá)到快速、簡(jiǎn)便地求解深度域地震子波的技術(shù)效果。

在一個(gè)具體實(shí)施場(chǎng)景，應(yīng)用本申請(qǐng)?zhí)峁┥疃扔虻卣鹱硬ǖ拇_定方法/裝置對(duì)某區(qū)域的深度域地震子波進(jìn)行提取。具體實(shí)施可以參閱下面內(nèi)容執(zhí)行。

步驟(1)：利用聲波測(cè)井?dāng)?shù)據(jù)的速度場(chǎng)數(shù)據(jù)v(z)和密度場(chǎng)數(shù)據(jù)ρ(z)計(jì)算得到基于井?dāng)?shù)據(jù)的深度域反射系數(shù)序列：r^d(z),z＝z1,z2,…,zn，其中，z代表深度。具體可以參閱圖3的在一個(gè)場(chǎng)景示例中獲取的井?dāng)?shù)據(jù)得到的深度域反射系數(shù)示意圖。

步驟(2)：選取井位置處的深度域成像數(shù)據(jù)：x^d(z),z＝z1,z2,…,zn。

步驟(3)：計(jì)算深度域反射系數(shù)序列r^d(z)的自相關(guān)，(即深度域反射系數(shù)序列的自相關(guān)參數(shù))。

步驟(4)：計(jì)算深度域反射系數(shù)序列r^d(z)與深度域成像數(shù)據(jù)x^d(z)的互相關(guān)μrx(0),μrx(1),…,μrx(l)(即深度域反射系數(shù)序列與深度域成像數(shù)據(jù)的互相關(guān)參數(shù))。

步驟(5)：具體實(shí)施時(shí)，可以再根據(jù)步驟(3)計(jì)算出的深度域反射系數(shù)序列的自相關(guān)與步驟(4)計(jì)算出的深度域反射系數(shù)序列和深度域成像數(shù)據(jù)的互相關(guān)μrx(0),μrx(1),…,μrx(l)，求取出以深度域成像數(shù)據(jù)x^d(z)為期望輸出，且基于井?dāng)?shù)據(jù)的深度域反射系數(shù)序列r^d(z)的最佳匹配算子w^d(z)，該算子就是最小二乘意義下的深度域等效子波，即為深度域地震子波。

具體施工時(shí)，

步驟(1)的計(jì)算方法如下：

深度zi處反射系數(shù)r^d(zi)可以由速度場(chǎng)v(z)和密度場(chǎng)ρ(z)計(jì)算得到，表示為：

步驟(3)的計(jì)算方法如下：

深度域反射系數(shù)序列r^d(z)的自相關(guān)計(jì)算可以表示為：

步驟(4)的計(jì)算方法如下：

深度域反射系數(shù)序列r^d(z)與深度域成像數(shù)據(jù)x^d(z)的互相關(guān)計(jì)算可以表示為：

步驟(5)的計(jì)算方法如下：

為求取最小二乘意義下的深度域等效子波w^d(z)，根據(jù)步驟(3)和步驟(4)求取的和μrx(z)構(gòu)建矩陣方程表達(dá)式如下：

需要說(shuō)明的是，上式矩陣方程具體可以通過(guò)levinson快速遞推算法求解出w^d(z)即為深度地震子波，其中，深度域地震子波的波長(zhǎng)度為l。且通過(guò)現(xiàn)有方法所提取的時(shí)間域地震子波，參閱圖4的在一個(gè)場(chǎng)景示例中獲取的時(shí)間域成像示意圖，與本實(shí)施方式通過(guò)深度域地震子波確定方法/裝置確定的深度域地震子波，參閱圖5的在一個(gè)場(chǎng)景示例中獲取的深度域子波示意圖，比較發(fā)現(xiàn)兩者波形并不相同。進(jìn)而，可以根據(jù)深度域地震子波，確定深度域地震子波合成記錄，參閱圖6的在一個(gè)場(chǎng)景示例中獲取的深度域等效子波合成記錄。與通過(guò)時(shí)間域地震子波的合成記錄進(jìn)行比較，參閱圖7的在一個(gè)場(chǎng)景示例中獲取的深度域等效子波合成記錄與時(shí)深轉(zhuǎn)換記錄的對(duì)比示意圖。根據(jù)深度域地震子波所確定的合成記錄要更加的準(zhǔn)確。

因此，通過(guò)上述場(chǎng)景示例，驗(yàn)證了應(yīng)用本申請(qǐng)實(shí)施方式提供的深度域地震子波的確定方法和裝置確實(shí)可以解決現(xiàn)有的深度域地震子波的確定方法中存在的無(wú)法直接求取深度域地震子波、確定過(guò)程復(fù)雜、所確定的深度域地震子波準(zhǔn)確度差的技術(shù)問(wèn)題，達(dá)到準(zhǔn)確、簡(jiǎn)便地確定深度域地震子波的技術(shù)效果。實(shí)現(xiàn)了可以直接在深度域提取地震子波，解決了深度域地震資料無(wú)法提取子波的難題，為疊前深度域資料在反演方面的直接應(yīng)用奠定了基礎(chǔ)。

盡管本申請(qǐng)內(nèi)容中提到不同的深度域地震地震波的確定方法或裝置，但是，本申請(qǐng)并不局限于必須是行業(yè)標(biāo)準(zhǔn)或?qū)嵤├枋龅那闆r等，某些行業(yè)標(biāo)準(zhǔn)或者使用自定義方式或?qū)嵤├枋龅膶?shí)施基礎(chǔ)上略加修改后的實(shí)施方案也可以實(shí)現(xiàn)上述實(shí)施例相同、等同或相近、或變形后可預(yù)料的實(shí)施效果。應(yīng)用這些修改或變形后的數(shù)據(jù)獲取、處理、輸出、判斷方式等的實(shí)施例，仍然可以屬于本申請(qǐng)的可選實(shí)施方案范圍之內(nèi)。

雖然本申請(qǐng)?zhí)峁┝巳鐚?shí)施例或流程圖所述的方法操作步驟，但基于常規(guī)或者無(wú)創(chuàng)造性的手段可以包括更多或者更少的操作步驟。實(shí)施例中列舉的步驟順序僅僅為眾多步驟執(zhí)行順序中的一種方式，不代表唯一的執(zhí)行順序。在實(shí)際中的裝置或客戶(hù)端產(chǎn)品執(zhí)行時(shí)，可以按照實(shí)施例或者附圖所示的方法順序執(zhí)行或者并行執(zhí)行(例如并行處理器或者多線程處理的環(huán)境，甚至為分布式數(shù)據(jù)處理環(huán)境)。術(shù)語(yǔ)“包括”、“包含”或者其任何其他變體意在涵蓋非排他性的包含，從而使得包括一系列要素的過(guò)程、方法、產(chǎn)品或者設(shè)備不僅包括那些要素，而且還包括沒(méi)有明確列出的其他要素，或者是還包括為這種過(guò)程、方法、產(chǎn)品或者設(shè)備所固有的要素。在沒(méi)有更多限制的情況下，并不排除在包括所述要素的過(guò)程、方法、產(chǎn)品或者設(shè)備中還存在另外的相同或等同要素。

上述實(shí)施例闡明的裝置或模塊等，具體可以由計(jì)算機(jī)芯片或?qū)嶓w實(shí)現(xiàn)，或者由具有某種功能的產(chǎn)品來(lái)實(shí)現(xiàn)。為了描述的方便，描述以上裝置時(shí)以功能分為各種模塊分別描述。當(dāng)然，在實(shí)施本申請(qǐng)時(shí)可以把各模塊的功能在同一個(gè)或多個(gè)軟件和/或硬件中實(shí)現(xiàn)，也可以將實(shí)現(xiàn)同一功能的模塊由多個(gè)子模塊的組合實(shí)現(xiàn)等。以上所描述的裝置實(shí)施例僅僅是示意性的，例如，所述模塊的劃分，僅僅為一種邏輯功能劃分，實(shí)際實(shí)現(xiàn)時(shí)可以有另外的劃分方式，例如多個(gè)模塊或組件可以結(jié)合或者可以集成到另一個(gè)系統(tǒng)，或一些特征可以忽略，或不執(zhí)行。

本領(lǐng)域技術(shù)人員也知道，除了以純計(jì)算機(jī)可讀程序代碼方式實(shí)現(xiàn)控制器以外，完全可以通過(guò)將方法步驟進(jìn)行邏輯編程來(lái)使得控制器以邏輯門(mén)、開(kāi)關(guān)、專(zhuān)用集成電路、可編程邏輯控制器和嵌入微控制器等的形式來(lái)實(shí)現(xiàn)相同功能。因此這種控制器可以被認(rèn)為是一種硬件部件，而對(duì)其內(nèi)部包括的用于實(shí)現(xiàn)各種功能的裝置也可以視為硬件部件內(nèi)的結(jié)構(gòu)。或者甚至，可以將用于實(shí)現(xiàn)各種功能的裝置視為既可以是實(shí)現(xiàn)方法的軟件模塊又可以是硬件部件內(nèi)的結(jié)構(gòu)。

本申請(qǐng)可以在由計(jì)算機(jī)執(zhí)行的計(jì)算機(jī)可執(zhí)行指令的一般上下文中描述，例如程序模塊。一般地，程序模塊包括執(zhí)行特定任務(wù)或?qū)崿F(xiàn)特定抽象數(shù)據(jù)類(lèi)型的例程、程序、對(duì)象、組件、數(shù)據(jù)結(jié)構(gòu)、類(lèi)等等。也可以在分布式計(jì)算環(huán)境中實(shí)踐本申請(qǐng)，在這些分布式計(jì)算環(huán)境中，由通過(guò)通信網(wǎng)絡(luò)而被連接的遠(yuǎn)程處理設(shè)備來(lái)執(zhí)行任務(wù)。在分布式計(jì)算環(huán)境中，程序模塊可以位于包括存儲(chǔ)設(shè)備在內(nèi)的本地和遠(yuǎn)程計(jì)算機(jī)存儲(chǔ)介質(zhì)中。

通過(guò)以上的實(shí)施方式的描述可知，本領(lǐng)域的技術(shù)人員可以清楚地了解到本申請(qǐng)可借助軟件加必需的通用硬件平臺(tái)的方式來(lái)實(shí)現(xiàn)?；谶@樣的理解，本申請(qǐng)的技術(shù)方案本質(zhì)上或者說(shuō)對(duì)現(xiàn)有技術(shù)做出貢獻(xiàn)的部分可以以軟件產(chǎn)品的形式體現(xiàn)出來(lái)，該計(jì)算機(jī)軟件產(chǎn)品可以存儲(chǔ)在存儲(chǔ)介質(zhì)中，如rom/ram、磁碟、光盤(pán)等，包括若干指令用以使得一臺(tái)計(jì)算機(jī)設(shè)備(可以是個(gè)人計(jì)算機(jī)，移動(dòng)終端，服務(wù)器，或者網(wǎng)絡(luò)設(shè)備等)執(zhí)行本申請(qǐng)各個(gè)實(shí)施例或者實(shí)施例的某些部分所述的方法。

本說(shuō)明書(shū)中的各個(gè)實(shí)施例采用遞進(jìn)的方式描述，各個(gè)實(shí)施例之間相同或相似的部分互相參見(jiàn)即可，每個(gè)實(shí)施例重點(diǎn)說(shuō)明的都是與其他實(shí)施例的不同之處。本申請(qǐng)可用于眾多通用或?qū)Ｓ玫挠?jì)算機(jī)系統(tǒng)環(huán)境或配置中。例如：個(gè)人計(jì)算機(jī)、服務(wù)器計(jì)算機(jī)、手持設(shè)備或便攜式設(shè)備、平板型設(shè)備、多處理器系統(tǒng)、基于微處理器的系統(tǒng)、置頂盒、可編程的電子設(shè)備、網(wǎng)絡(luò)pc、小型計(jì)算機(jī)、大型計(jì)算機(jī)、包括以上任何系統(tǒng)或設(shè)備的分布式計(jì)算環(huán)境等等。

雖然通過(guò)實(shí)施例描繪了本申請(qǐng)，本領(lǐng)域普通技術(shù)人員知道，本申請(qǐng)有許多變形和變化而不脫離本申請(qǐng)的精神，希望所附的權(quán)利要求包括這些變形和變化而不脫離本申請(qǐng)。