施例提供的一種面元覆蓋次數(shù)的計(jì)算方法的流程示意圖。
[0048] 圖4是觀測(cè)系統(tǒng)1在點(diǎn)位（10000, 3000)處的偏移覆蓋次數(shù)圖。
[0049] 圖5是觀測(cè)系統(tǒng)1在點(diǎn)位（11000, 3000)處的偏移覆蓋次數(shù)圖。
[0050] 圖6是觀測(cè)系統(tǒng)2在點(diǎn)位（11000, 3000)處的偏移覆蓋次數(shù)圖。
[0051] 圖7是本申請(qǐng)實(shí)施例提供的一種面元覆蓋次數(shù)的計(jì)算裝置的模塊示意圖。
【具體實(shí)施方式】
[0052] 本申請(qǐng)實(shí)施例提供一種面元覆蓋次數(shù)的計(jì)算方法及裝置。
[0053] 為了使本技術(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ù) 的范圍。
[0054] 地震偏移成像是地震數(shù)據(jù)處理三大技術(shù)（地震反褶積、疊加和偏移成像）之一。地 震偏移就是指要采用特定的方法實(shí)現(xiàn)反射點(diǎn)和繞射點(diǎn)正確的空間歸位，即成像，故又可以 稱(chēng)為偏移歸位或偏移成像。
[0055] 地震偏移可以使傾斜反射歸位到它們真正的地下界面位置，并使繞射波收斂，以 此提高空間分辨率和得到地下界面的地震圖像。地震偏移的方法可以分為二維偏移和三維 偏移；又可以分為疊后偏移和疊前偏移；還可以分為深度偏移和時(shí)間偏移。
[0056] 對(duì)于地震偏移成像，偏移算法主要是指將反射能量或繞射能量（也可以統(tǒng)稱(chēng)為歸 位能量）進(jìn)行分配再疊加。
[0057] 本申請(qǐng)實(shí)施例中主要是利用疊前偏移方法來(lái)計(jì)算面元覆蓋次數(shù)。
[0058] 下面結(jié)合附圖對(duì)本申請(qǐng)實(shí)施例所提供的一種面元覆蓋次數(shù)的計(jì)算方法進(jìn)行詳細(xì) 說(shuō)明。雖然本申請(qǐng)?zhí)峁┝巳缦率鰧?shí)施例或流程圖所述的方法操作步驟，但基于常規(guī)或者無(wú) 需創(chuàng)造性的勞動(dòng)在所述方法中可以包括更多或者更少的操作步驟，操作步驟之間的執(zhí)行順 序并沒(méi)有限制。如圖3所示，本申請(qǐng)實(shí)施例所提供的一種面元覆蓋次數(shù)的計(jì)算方法包括如 下步驟：
[0059] S110:將所獲取的預(yù)設(shè)觀測(cè)系統(tǒng)中每個(gè)炮檢對(duì)所對(duì)應(yīng)的反射能量進(jìn)行疊前偏移計(jì) 算。
[0060] 在對(duì)目標(biāo)區(qū)域進(jìn)行勘探施工前期，可以根據(jù)目標(biāo)區(qū)域的野外勘探任務(wù)和對(duì)以往該 目標(biāo)區(qū)域的勘探資料的分析，針對(duì)該目標(biāo)區(qū)域預(yù)設(shè)一個(gè)或多個(gè)（例如N個(gè)，N為正整數(shù)）觀 測(cè)系統(tǒng)，然后可以對(duì)每個(gè)觀測(cè)系統(tǒng)中每個(gè)炮檢對(duì)所對(duì)應(yīng)的反射能量進(jìn)行疊前偏移計(jì)算。具 體的，
[0061] 可以首先依次建立第一觀測(cè)系統(tǒng)中每個(gè)炮檢對(duì)（即一個(gè)炮點(diǎn)與對(duì)應(yīng)的一個(gè)檢波 點(diǎn)所形成的組合）所對(duì)應(yīng)的反射能量進(jìn)行疊前偏移計(jì)算的偏移路徑；然后按照所建立的偏 移路徑，依次利用疊前偏移算法對(duì)該觀測(cè)系統(tǒng)中每個(gè)炮檢對(duì)所對(duì)應(yīng)的反射能量進(jìn)行疊前偏 移計(jì)算；再按照對(duì)該觀測(cè)系統(tǒng)中每個(gè)炮檢對(duì)所對(duì)應(yīng)的反射能量進(jìn)行疊前偏移計(jì)算的方法， 依次對(duì)其他N-I個(gè)觀測(cè)系統(tǒng)中每個(gè)炮檢對(duì)所對(duì)應(yīng)的反射能量進(jìn)行疊前偏移計(jì)算。
[0062] 也可以首先建立第一觀測(cè)系統(tǒng)中第一炮檢對(duì)所對(duì)應(yīng)的反射能量進(jìn)行偏移計(jì)算的 偏移路徑；然后按照所建立的偏移路徑，利用疊前偏移算法對(duì)該炮檢對(duì)所對(duì)應(yīng)的反射能量 進(jìn)行偏移計(jì)算；再按照對(duì)該炮檢對(duì)所對(duì)應(yīng)的反射能量進(jìn)行偏移計(jì)算的方法，依次建立該觀 測(cè)系統(tǒng)中其他炮檢對(duì)（例如剩余的M-I個(gè)炮檢對(duì)，M為正整數(shù)）所對(duì)應(yīng)的反射能量的偏移 路徑并進(jìn)行疊前偏移計(jì)算；最后按照對(duì)該觀測(cè)系統(tǒng)中每個(gè)炮檢對(duì)所對(duì)應(yīng)的反射能量進(jìn)行疊 前偏移計(jì)算的方法，依次對(duì)其他N-I個(gè)觀測(cè)系統(tǒng)中每個(gè)炮檢對(duì)所對(duì)應(yīng)的反射能量進(jìn)行疊前 偏移計(jì)算。
[0063] 所述觀測(cè)系統(tǒng)可以是根據(jù)目標(biāo)區(qū)域的野外勘探任務(wù)和對(duì)以往該目標(biāo)區(qū)域的勘探 資料的分析來(lái)進(jìn)行預(yù)設(shè)的。具體的，可以根據(jù)目標(biāo)區(qū)域的勘探任務(wù)，設(shè)計(jì)合適的觀測(cè)系統(tǒng)參 數(shù)（包括面元、炮線距、最大偏移距、覆蓋次數(shù)、接收線距等）；然后將目標(biāo)區(qū)域中的炮點(diǎn)和 檢波點(diǎn)按照所設(shè)計(jì)的觀測(cè)系統(tǒng)參數(shù)進(jìn)行布設(shè)。
[0064] 所述目標(biāo)區(qū)域可以是整個(gè)勘探施工區(qū)域，也可以是勘探施工區(qū)域中的部分區(qū)域。 所述勘探任務(wù)可以包括所述目標(biāo)區(qū)域的覆蓋范圍、采樣密度以及地震波的激發(fā)方式等。
[0065] 可以按照下述要求來(lái)設(shè)計(jì)所述觀測(cè)系統(tǒng)參數(shù)：
[0066] (1)面元
[0067] 可以根據(jù)所述目標(biāo)區(qū)域中目標(biāo)地質(zhì)體的大小來(lái)設(shè)計(jì)面元，例如所設(shè)計(jì)面元的邊長(zhǎng) 不得大于目標(biāo)地質(zhì)體邊長(zhǎng)的三分之一。
[0068] 所設(shè)計(jì)的面元可能除了滿足最高無(wú)混疊頻率法則、橫向分辨率的要求、考慮去噪 無(wú)空間假頻要求、陡傾角地層及斷層偏移成像要求外，還可能要確保波場(chǎng)空間采樣充分、有 利于地震波場(chǎng)歸位以及提高目的層成像精度等。
[0069] (2)最大偏移距
[0070] 最大偏移距也可以稱(chēng)為最大炮檢距?？梢栽跐M足目的層埋深、大于6 %的速度分析 精度、小于12. 5%動(dòng)校拉伸、反射系數(shù)穩(wěn)定性對(duì)排列長(zhǎng)度的限制以及偏移時(shí)可以歸位95% 的繞射波能量等需求的基礎(chǔ)上，結(jié)合波動(dòng)方程正演、波場(chǎng)照明、實(shí)際地震資料道集、速度分 析、以往二維地震資料的偏移距成像剖面等來(lái)設(shè)計(jì)最大偏移距。
[0071] (3)接收線距
[0072] 可以根據(jù)菲涅爾帶半徑和最大非縱距來(lái)設(shè)計(jì)接收線距。所設(shè)計(jì)的接收線距一般不 大于第一菲涅爾帶半徑。所述最大非縱距需要滿足三維地震資料同一面元內(nèi)不同方位角的 反射同相疊加。
[0073] (4)覆蓋次數(shù)
[0074] 可以考慮目標(biāo)區(qū)域中構(gòu)造的復(fù)雜程度、表層激發(fā)、接收條件、激發(fā)方式以及接收方 式來(lái)預(yù)設(shè)每個(gè)面元的覆蓋次數(shù)。所述覆蓋次數(shù)可以是根據(jù)呈水平層狀的介質(zhì)來(lái)預(yù)設(shè)的。
[0075] 所述將目標(biāo)區(qū)域中的炮點(diǎn)和檢波點(diǎn)按照所設(shè)計(jì)的觀測(cè)系統(tǒng)參數(shù)進(jìn)行布設(shè)可以是 將所述目標(biāo)區(qū)域劃分成多個(gè)子區(qū)域，然后在某一個(gè)子區(qū)域中按照所設(shè)計(jì)的觀測(cè)系統(tǒng)參數(shù)， 對(duì)炮點(diǎn)和檢波點(diǎn)進(jìn)行排布，并建立炮點(diǎn)和檢波點(diǎn)之間的對(duì)應(yīng)關(guān)系，最后按照在該子區(qū)域中 布設(shè)炮點(diǎn)和檢波點(diǎn)的方法，在其他子區(qū)域中對(duì)對(duì)炮點(diǎn)和檢波點(diǎn)進(jìn)行布設(shè)。所述將目標(biāo)區(qū)域 中的炮點(diǎn)和檢波點(diǎn)按照所設(shè)計(jì)的觀測(cè)系統(tǒng)參數(shù)進(jìn)行布設(shè)也可以是首先在整個(gè)目標(biāo)區(qū)域中 按照所設(shè)計(jì)的觀測(cè)系統(tǒng)參數(shù)對(duì)炮點(diǎn)和檢波點(diǎn)進(jìn)行排布，然后建立整個(gè)目標(biāo)區(qū)域中炮點(diǎn)和檢 波點(diǎn)之間的對(duì)應(yīng)關(guān)系。所述對(duì)應(yīng)關(guān)系可以是指每個(gè)炮點(diǎn)與接收該炮點(diǎn)激發(fā)的地震波的檢波 點(diǎn)之間的對(duì)應(yīng)關(guān)系。
[0076] 所述反射能量可以是指炮點(diǎn)處激發(fā)的地震波經(jīng)過(guò)反射后，檢波器所接收到的地震 波能量。在本申請(qǐng)實(shí)施例中，所述反射能量可以為預(yù)設(shè)的觀測(cè)系統(tǒng)參數(shù)中的覆蓋次數(shù)；也可 以為炮檢對(duì)的個(gè)數(shù)。一般對(duì)于某一個(gè)炮檢對(duì)，可以將該炮檢對(duì)所對(duì)應(yīng)的反射能量視為1。
[0077] 所述偏移路徑可以是指對(duì)所述反射能量進(jìn)行偏移計(jì)算時(shí)，將反射點(diǎn)和繞射點(diǎn)歸位 到真實(shí)位置所模擬的虛擬路徑。在所述偏移路徑上可以分布有多個(gè)面元。所述偏移路徑可 以是以預(yù)設(shè)點(diǎn)位為圓心，以偏移孔徑為半徑所構(gòu)建的橢圓路徑。所述預(yù)設(shè)點(diǎn)位可以為給定 的繞射點(diǎn)，也可以為其他給定坐標(biāo)位置的點(diǎn)。
[0078] 所述偏移孔徑可以定義為對(duì)成像點(diǎn)有貢獻(xiàn)的所有地震數(shù)據(jù)在地面投影點(diǎn)的集合。 所述偏移孔徑可以根據(jù)以下3個(gè)因素來(lái)進(jìn)行設(shè)定：①第一菲涅爾帶半徑；②傾斜層偏移歸 位橫向移動(dòng)距離（傾角），即Zxtan機(jī) m，其中Z為最深目的層的深度，胃為最深目的層的 最大傾角；③滿足30°繞射收斂，即大于ZXtan30°。所述偏移孔徑一般可以為滿足上述 3個(gè)因素的最大值，也可以在滿足上述3個(gè)因素的基礎(chǔ)上，根據(jù)實(shí)際情況來(lái)進(jìn)行設(shè)定。
[0079] 按照所述偏移路徑，對(duì)每一個(gè)炮檢對(duì)所對(duì)應(yīng)的反射能量進(jìn)行偏移計(jì)算可以是指利 用疊前偏移算法，將該炮檢對(duì)所對(duì)應(yīng)的反射能量分散到所述偏移路徑中的每個(gè)面元內(nèi)。
[0080] 所述疊前偏移算法可以是基于繞射掃描疊加原理的Kirchhoff (克?；舴颍┓e分 法；也可以是基于遞推波場(chǎng)延拓的波動(dòng)方程有限差分法。所述疊前偏移算法不限于上述兩 種方法，也可以是其他方法。
[0081] 所述利用疊前偏移算法，將該炮檢對(duì)所對(duì)應(yīng)的反射能量分散到所述偏移路徑中的 每個(gè)面元內(nèi)可以是將該炮檢對(duì)所對(duì)應(yīng)的反射能量按照疊前偏移公式的權(quán)系數(shù)分散到每個(gè) 面元內(nèi)，即將反射能量與所述偏移路徑中每個(gè)面元所對(duì)應(yīng)的權(quán)系數(shù)進(jìn)行乘積運(yùn)算，得到該 反射能量分散到每個(gè)面元內(nèi)的反射能量分量。
[0082] 所述權(quán)系數(shù)可以是地層傾角與地震波傳播速度和炮檢對(duì)的旅行路徑之間的比值， 其可以表示如下：
[0083]
(1)
[0084] 上式中，Wei