基于目標(biāo)層可視性分析的激發(fā)位置確定方法
【專利摘要】本發(fā)明涉及一種通過地震波針對目標(biāo)反射層的可視性分析來確定地震波最優(yōu)激發(fā)位置的方法。本發(fā)明的基于目標(biāo)層可視性分析的激發(fā)位置確定方法包括：（a）從地質(zhì)模型中選取目標(biāo)反射層；（b）在地質(zhì)模型的整個(gè)表面上均勻地布置多個(gè)勘探炮和多個(gè)檢波器；（c）并且當(dāng)多個(gè)勘探炮中的一個(gè)勘探炮產(chǎn)生地震波時(shí)，記錄每個(gè)檢波器接收到的能量值，從而得到每個(gè)檢波器分別針對當(dāng)前產(chǎn)生地震波的勘探炮的炮檢對的能量值；（d）在所有勘探炮分別產(chǎn)生相同次數(shù)的地震波之后，基于多個(gè)檢波器分別針對每個(gè)勘探炮的炮檢對的能量值，繪制基于地質(zhì)模型的炮檢對可視性分析圖；（e）基于炮檢對可視性分析圖，確定在實(shí)際的地震勘探中應(yīng)該增加的勘探炮的最佳位置。
【專利說明】基于目標(biāo)層可視性分析的激發(fā)位置確定方法
【技術(shù)領(lǐng)域】
[0001]本發(fā)明涉及一種從定量的角度去指導(dǎo)地質(zhì)觀測系統(tǒng)的方法，，尤其涉及一種通過地震波針對目標(biāo)反射層的可視性分析來確定地震波激發(fā)位置的方法。
【背景技術(shù)】
[0002]隨著我國油氣勘探向西部和西南部的轉(zhuǎn)移，復(fù)雜地表及復(fù)雜構(gòu)造問題變得非常突出，地震資料信噪比低、地震成像困難等問題嚴(yán)重影響著這些地區(qū)的勘探效果。通過研究地震波在復(fù)雜地質(zhì)條件下的傳播規(guī)律與能量分布特征，進(jìn)而能夠?qū)崿F(xiàn)觀測系統(tǒng)的優(yōu)化設(shè)計(jì)具有重要的現(xiàn)實(shí)意義。
[0003]針對地震觀測系統(tǒng)的設(shè)計(jì)和采集參數(shù)的優(yōu)化是整個(gè)石油勘探中非常重要的環(huán)節(jié)，是地震資料品質(zhì)好壞的先決條件，也是地震資料處理和解釋的基礎(chǔ)。常規(guī)地震觀測系統(tǒng)的設(shè)計(jì)的覆蓋次數(shù)只依賴于炮點(diǎn)和檢波點(diǎn)的排列范圍，而與地下結(jié)構(gòu)、目的反射層深度無關(guān)。然而，在構(gòu)造復(fù)雜的地區(qū)，覆蓋次數(shù)除了與炮點(diǎn)和檢波點(diǎn)排列有關(guān)外，還與目的反射層的深度以及上覆地質(zhì)結(jié)構(gòu)有關(guān)，此時(shí)，常規(guī)的基于水平層狀假設(shè)的針對CMP點(diǎn)(共中心點(diǎn))進(jìn)行設(shè)計(jì)的觀測系統(tǒng)精度降低。由于目的發(fā)射層各CRP點(diǎn)(共反射點(diǎn))上的覆蓋次數(shù)和地震波照射能量共同決定了該點(diǎn)的成像質(zhì)量，但是復(fù)雜的上覆地質(zhì)結(jié)構(gòu)會造成下伏地層地震照明陰影區(qū)，從而引起下伏勘探目的層照明強(qiáng)度顯著下降，最終使得這些目的層界面成像的效果不理想?；诓▌臃匠痰挠^測系統(tǒng)分析方法目前主要是照明度分析，而地震波照明又可以分為基于射線追蹤的地震波照明方法和基于波動方程的地震波照明方法兩類。射線追蹤地震波照明方法目前應(yīng)用最為普遍，該方法簡單直觀、計(jì)算效率高，但由于射線理論本身的缺陷(如高頻近似、射線盲區(qū)等)和復(fù)雜介質(zhì)中的奇異性問題，基于射線的照明分析方法在復(fù)雜構(gòu)造區(qū)會出現(xiàn)很大的誤差，照明精度較低。相反，基于波動方程的地震波照明方法能有效地克服射線照明分析方法在復(fù)雜構(gòu)造區(qū)的不足，尤其是利用波動方程中的單程波方程的地下照明分析是近幾年的一個(gè)研究熱點(diǎn)，它克服了射線照明方法的缺陷，能適應(yīng)強(qiáng)橫向變速介質(zhì)，使得目標(biāo)照明分析更加合理、準(zhǔn)確。但是，基于波動方程的照明度計(jì)算是分析能量在地下介質(zhì)分布的情況，并且不會考慮從目標(biāo)反射層和檢波器接收能量針對觀測系統(tǒng)的影響，使得照明度分析在指導(dǎo)觀測系統(tǒng)設(shè)計(jì)更多的是定性的角度，無法從定量的角度指導(dǎo)觀測系統(tǒng)的優(yōu)化設(shè)計(jì)。
[0004]因此，需要一種能夠適應(yīng)復(fù)雜地表結(jié)構(gòu)且可以從定量的角度去指導(dǎo)地質(zhì)觀測系統(tǒng)的地震波激發(fā)位置的方法。

【發(fā)明內(nèi)容】

[0005]本發(fā)明的目的在于提出一種基于目標(biāo)層可視性分析的激發(fā)位置確定方法。
[0006]根據(jù)本發(fā)明的一實(shí)施例的基于目標(biāo)層可視性分析的激發(fā)位置確定方法，包括以下步驟:(a)從地質(zhì)模型中選取目標(biāo)反射層；(b)在地質(zhì)模型的整個(gè)表面上均勻地布置多個(gè)勘探炮和多個(gè)檢波器；(C)使所有勘探炮分別產(chǎn)生地震波，并且當(dāng)多個(gè)勘探炮中的一個(gè)勘探炮產(chǎn)生地震波時(shí)，多個(gè)檢波器分別接收該地震波被目標(biāo)反射層反射回來的反射波，同時(shí)記錄每個(gè)檢波器接收到的能量值，從而得到每個(gè)檢波器分別針對當(dāng)前產(chǎn)生地震波的勘探炮的炮檢對的能量值；(d)在所有勘探炮分別產(chǎn)生相同次數(shù)的地震波之后，基于多個(gè)檢波器分別針對每個(gè)勘探炮的炮檢對的能量值，繪制基于地質(zhì)模型的炮檢對可視性分析圖；(e)基于炮檢對可視性分析圖，確定在實(shí)際的地震勘探中應(yīng)該增加勘探炮數(shù)量的范圍，其中，炮檢對的能量值是指一個(gè)檢波器接收到的一個(gè)地震波的能量值 。
[0007]優(yōu)選地，在步驟(C)中，利用有限差分算子，對當(dāng)前產(chǎn)生的地震波進(jìn)行正向延拓到目標(biāo)反射層，然后再從目標(biāo)反射層反向延拓到檢波器，來模擬地震波被目標(biāo)反射層反射回來的反射波，從而記錄每個(gè)檢波器接收到的反射波的能量值。
[0008]優(yōu)選地，在步驟(C)中，當(dāng)產(chǎn)生地震波時(shí)，地震波在頻率域中以下面的公式(I)的形式傳播，
[0009](k2+ Δ ) U (rs, r, ω) =-f (ω) δ (r_rs)公式(I)
[0010]其中，k= ω Zc0 (r)為波數(shù)，ω表示地震波的頻率，C0(r)為位于目標(biāo)反射層的點(diǎn)r=(x, y, z)處的背景場速度值，Λ為拉普拉斯算子，f(?)為勘探炮函數(shù)，δ (r-rs)是單位脈沖函數(shù)，rs表示勘探炮的三維坐標(biāo)。
[0011]優(yōu)選地，在步驟(C)中，利用有限差分算子，對當(dāng)前產(chǎn)生的地震波進(jìn)行正向延拓到目標(biāo)反射層，然后再從目標(biāo)反射層反向延拓到檢波器，
[0012]其中，檢波器接收到的反向延拓的反射波的地震波場通過公式(2)來表示，
[0013]P (rs, rg, ω)=
[0014]2k2 / v(r)G(r，rs; co)m(r)G(r，rg; co)d3r 公式(2)
[0015]其中，G(r，rs; ω)和G(r，rg; ω)分別是頻率域中勘探炮和檢波器位置的格林函數(shù)，&表示檢波器位置的三維坐標(biāo)，m(r)是位于目標(biāo)反射層的點(diǎn)r=(x，y，z)處速度模型擾動，V(r)為目標(biāo)反射層的所有點(diǎn)的集合。數(shù)學(xué)上還可有如下的表示m(r)= δ c/c(r), c (r)是地質(zhì)模型的速度值，δ c是對速度值的求導(dǎo)，但是本發(fā)明中沒有使用該式來求出m(r)。
[0016]優(yōu)選地，多個(gè)檢波器中的一個(gè)檢波器針對當(dāng)前產(chǎn)生地震波的勘探炮的炮檢對的能量值通過公式(3)來表示，
[0017]Psr (rs, rg, ω) = |P(rs, rg, ω) 公式(3)。
[0018]本發(fā)明從基于目標(biāo)反射層可視性分析的角度對勘探炮位置(即，激發(fā)位置)進(jìn)行了定量評估，對在后續(xù)的地質(zhì)勘探中提高山地地震采集設(shè)計(jì)效率及優(yōu)化提供了保障，具有重要應(yīng)用價(jià)值。
[0019]本發(fā)明基于可視性分析計(jì)算角度，確定了對目標(biāo)反射層的影響較大的觀測系統(tǒng)的激發(fā)位置，對進(jìn)行山地復(fù)雜構(gòu)造的地震資料采集觀測系統(tǒng)設(shè)計(jì)的激發(fā)位置的優(yōu)選提供了保障，進(jìn)而為山地復(fù)雜條件下的偏移成像、油氣儲層預(yù)測、油藏描述的可靠性奠定了基礎(chǔ)，具有重要應(yīng)用價(jià)值。
[0020]根據(jù)本發(fā)明，可以確定能夠更好地獲取目標(biāo)反射層的反射信息的地表位置，即在這些優(yōu)選的地表上布置勘探炮時(shí)可以得到針對目標(biāo)反射層的更好的成像效果。
【專利附圖】

【附圖說明】
[0021]通過下面結(jié)合附圖進(jìn)行的對實(shí)施例的描述，本發(fā)明的上述和/或其他目的和優(yōu)點(diǎn)將會變得更加清楚，其中:
[0022]圖1是示出本發(fā)明的基于目標(biāo)層可視性分析的激發(fā)位置確定方法的流程圖。
[0023]圖2是示出二維地質(zhì)模型的示意圖。
[0024]圖3是示出地震波傳播和反射過程的示意圖。
[0025]圖4是示出一個(gè)勘探炮所產(chǎn)生的地震波被目標(biāo)反射層反射之后由檢波器所接收到的反射波的可視性能量分布圖。
[0026]圖5是示出在顯示單元上顯示基于地質(zhì)模型的炮檢對可視性分析圖。
[0027]圖6是示出根據(jù)本發(fā)明的方法對圖2的二維地質(zhì)模型的目標(biāo)反射層進(jìn)行的實(shí)際應(yīng)用的圖，圖6的(a)是示出在優(yōu)選的激發(fā)位置增加勘探炮后得到的偏移剖面的圖，圖6的(b)是示出在優(yōu)選的激發(fā)位置之外增加勘探炮后得到的偏移剖面的圖。
【具體實(shí)施方式】
[0028]以下，參照附圖來詳細(xì)說明本發(fā)明的實(shí)施例。
[0029]本發(fā)明利用有限差分算子，對勘探炮所產(chǎn)生的進(jìn)行地震波場正向延拓到目標(biāo)反射層，然后再反向延拓到檢波器，記錄每個(gè)檢波器接收到的反射波的能量值進(jìn)而繪制當(dāng)前地質(zhì)模型的炮檢對可視性分析圖，從而基于炮檢對可視性分析圖來確定在檢波器處的能量貢獻(xiàn)較大的地震波激發(fā)位置。
[0030]以下，詳細(xì)說明本發(fā)明的各個(gè)流程。
[0031]圖1是示出本發(fā)明的基于目標(biāo)層可視性分析的激發(fā)位置確定方法的流程圖。圖2是示出二維地質(zhì)模型的示意圖。
[0032]在圖1的步驟S101，從地質(zhì)模型中選擇目標(biāo)反射層。本發(fā)明中選擇如圖2所示的二維地質(zhì)模型中的某一目標(biāo)反射層，并且該地質(zhì)模型的任何地層都可被選擇為目標(biāo)反射層。并且，根據(jù)已知地質(zhì)信息建立二維地質(zhì)模型進(jìn)而分析各種參數(shù)的選取是地球物理領(lǐng)域通用的方法。
[0033]在步驟S102，在地質(zhì)模型的整個(gè)表面上均勻地布置多個(gè)勘探炮和多個(gè)檢波器。在此，只有在地質(zhì)模型的整個(gè)表面上均勻地布置勘探炮和檢波器，才能使每個(gè)勘探炮產(chǎn)生的地震波對炮檢對可視性分析圖的影響是相同的。并且，在整個(gè)表面上僅布置一個(gè)勘探炮和一個(gè)檢波器是沒有太大價(jià)值的，因此在此步驟中將不會考慮這種情況。
[0034]在步驟S103，使所有勘探炮分別產(chǎn)生地震波，并且當(dāng)多個(gè)勘探炮中的一個(gè)勘探炮產(chǎn)生地震波時(shí)，多個(gè)檢波器分別接收該地震波被目標(biāo)反射層反射回來的反射波，同時(shí)記錄每個(gè)檢波器接收到的能量值，從而得到每個(gè)檢波器分別針對當(dāng)前產(chǎn)生地震波的勘探炮的炮檢對的能量值。其中，炮檢對的能量值是指一個(gè)檢波器接收到的一個(gè)地震波的能量值。并且，在地質(zhì)模型的整個(gè)表面上均勻地布置多個(gè)勘探炮和多個(gè)檢波器的前提下，為使每個(gè)勘探炮產(chǎn)生的地震波對炮檢對可視性分析圖起到相同權(quán)值的影響，使所有勘探炮分別產(chǎn)生相同次數(shù)的地震波。
[0035]在本發(fā)明中，是基于勘探炮處產(chǎn)生的地震波對目標(biāo)反射層的貢獻(xiàn)來確定貢獻(xiàn)度較大的勘探炮，因此為了使多個(gè)勘探炮對目標(biāo)反射層有相同權(quán)值的影響，需要使每個(gè)勘探炮產(chǎn)生的地震波的數(shù)量相同。在此，雖然在本發(fā)明中使所有勘探炮各產(chǎn)生一次地震波，但是根據(jù)需要，可使所有勘探炮分別產(chǎn)生多次(例如，2次，3次。。。η次等)地震波。[0036]另外，在產(chǎn)生地震波時(shí)，勘探炮位置以及地震波的激發(fā)位置(B卩，地震炮和地震波的產(chǎn)生位置)實(shí)際上是不同的，但是相比其他諸如數(shù)Km的距離，位置差別很小，因此在本領(lǐng)域中一般認(rèn)為是相同的。
[0037]并且，當(dāng)產(chǎn)生地震波時(shí)，在頻率域中且密度固定的地震波能夠以下面的公式(I)形式傳播，
[0038](k2+ Δ ) U (rs, r, ω) =-f (ω) δ (r_rs)公式(I)
[0039]其中，k= ω /c0 (r)為波數(shù)，ω表示地震波的頻率，C0 (r)為點(diǎn)r=(x, y, z)處的背景場速度值，Λ為拉普拉斯算子，f(co)為勘探炮函數(shù)，δ (r-rs)是單位脈沖函數(shù)，rs表示勘探炮的三維坐標(biāo)。而且，公式(I)為本領(lǐng)域公知的地震波的傳播數(shù)學(xué)式，因此在此不對此進(jìn)行詳細(xì)的說明。
[0040]圖3是示出地震波傳播和反射過程的示意圖。因此，在步驟S103中，利用有限差分算子，對當(dāng)前產(chǎn)生的地震波進(jìn)行正向延拓到目標(biāo)反射層，然后再從目標(biāo)反射層反向延拓到檢波器，此時(shí)記錄每個(gè)檢波器接收到的針對當(dāng)前產(chǎn)生的地震波的能量值。其中，有限差分算子作為地震波場的傳播算子，其實(shí)際數(shù)學(xué)表示為一個(gè)矩陣，表示地震波的實(shí)際正向傳播過程。并且反向延拓可以看成，將從勘探炮正向延拓到目標(biāo)反射層的地震波當(dāng)作二次勘探炮，進(jìn)而從目標(biāo)反射層到檢波器的正向延拓。
[0041]此時(shí)，檢波器接收到的地震波場(被目標(biāo)反射層反射回來的反射波場)通過公式
(2)來表示，
[0042]P (rs, rg, ω)=
[0043]2k2 / v(r)G(r，rs; co)m(r)G(r，rg; co)d3r 公式(2)
[0044]其中，G(r, rs; ω)和G(r，rg; ω)分別是頻率域中勘探炮和檢波器位置的格林函數(shù)，&表示檢波器位置的三維坐標(biāo)，m(r)是位于目標(biāo)反射層的點(diǎn)r=(x，y，z)處速度模型擾動，V (r)為目標(biāo)反射層的所有點(diǎn)的集合。數(shù)學(xué)上還可有如下的表示m(r)= δ c/c(r), c (r)是地質(zhì)模型的速度值，δ c是對速度值的求導(dǎo)，但是本發(fā)明中沒有使用該式來求出m(r)。
[0045]并且，多個(gè)檢波器中的一個(gè)檢波器針對當(dāng)前產(chǎn)生地震波的勘探炮的炮檢對的能量值通過公式(3)來表示，
[0046]Psr (rs, rg, ω) = | P (rs, rg, ω) 公式(3)。
[0047]圖4是示出一個(gè)勘探炮所產(chǎn)生的地震波被目標(biāo)反射層反射之后由檢波器所接收到的反射波的能見度(可視性能量)分布圖。圖5是示出基于地質(zhì)模型的炮檢對可視性分析圖。
[0048]在步驟S104，在所有勘探炮分別產(chǎn)生相同次數(shù)的地震波之后，基于多檢波器分別針對每個(gè)勘探炮的炮檢對的能量值，繪制根據(jù)地質(zhì)模型的炮檢對可視性分析圖。
[0049]具體地，一個(gè)勘探炮所產(chǎn)生的地震波被目標(biāo)反射層反射之后由檢波器所接收到的反射波的可視性能量分布圖如圖4所示，并且圖5以圖2的二維地質(zhì)模型的左上角為原點(diǎn)，且檢波點(diǎn)離原點(diǎn)的距離為橫坐標(biāo)，勘探炮離原點(diǎn)的距離為縱坐標(biāo)。因此，在繪制圖5所示的炮檢對可視性分析圖時(shí)，將同一個(gè)檢波器在所有勘探炮產(chǎn)生地震波時(shí)所接收的炮檢對的能量值加起來(簡稱，能量值總和)，然后將該檢波器離原點(diǎn)的距離作為橫坐標(biāo)，將當(dāng)前產(chǎn)生地震波的勘探炮離原點(diǎn)的距離作為縱坐標(biāo)，并且，針對每個(gè)檢波器的能量值總和的大小在圖5上能夠用不同的顏色來表示(在顯示單元上用顏色來表示能量大小是可行的)，例如可用紅色代表能量最強(qiáng)，藍(lán)色代表能量最弱，此時(shí)為了提高用戶體檢，可使用顯示器來顯示彩色的炮檢對可視性分析圖。當(dāng)然，也可以用諸如直接標(biāo)注能量值、點(diǎn)的密度等來表示能量值的大小。由此，圖5的每個(gè)點(diǎn)表示同一個(gè)檢波器針對所有地震波(從目標(biāo)反射層反射回來的反射波)記錄的能量值的總和。
[0050]在步驟S105，基于炮檢對可視性分析圖，確定在實(shí)際的地震勘探中應(yīng)該增加勘探炮數(shù)量的范圍。在圖5中，用半圓形實(shí)線包圍的位置即為應(yīng)該增加勘探炮數(shù)量的范圍。
[0051]圖6是示出根據(jù)本發(fā)明的方法對圖2的二維地質(zhì)模型的目標(biāo)反射層進(jìn)行的實(shí)際應(yīng)用的圖，圖6的(a)是示出在優(yōu)選的激發(fā)位置增加勘探炮后得到的偏移剖面的圖，圖6的(b)是示出在優(yōu)選的激發(fā)位置之外增加勘探炮后得到的偏移剖面的圖。
[0052]通過實(shí)際應(yīng)用可知，通過本發(fā)明可以確定更好地獲取目標(biāo)反射層的反射信息的地表位置(即，優(yōu)選的激發(fā)位置)，因此在后續(xù)實(shí)際的地質(zhì)勘探中向這些優(yōu)選的激發(fā)位置增加勘探炮能夠獲取更準(zhǔn)確的目標(biāo)反射層的信息。
[0053]本發(fā)明不限于上述實(shí)施例，在不脫離本發(fā)明范圍的情況下，可以進(jìn)行各種變形和修改。
[0054]產(chǎn)業(yè)上的可利用性
[0055]本發(fā)明可用在對于地質(zhì)觀測系統(tǒng)的設(shè)計(jì)以及評價(jià)。
【權(quán)利要求】
1.一種基于目標(biāo)層可視性分析的激發(fā)位置確定方法，其特征在于，包括以下步驟: Ca)從地質(zhì)模型中選取目標(biāo)反射層； (b)在地質(zhì)模型的整個(gè)表面上均勻地布置多個(gè)勘探炮和多個(gè)檢波器； (C)使所有勘探炮分別產(chǎn)生地震波，并且當(dāng)多個(gè)勘探炮中的一個(gè)勘探炮產(chǎn)生地震波時(shí)，多個(gè)檢波器分別接收該地震波被目標(biāo)反射層反射回來的反射波，同時(shí)記錄每個(gè)檢波器接收到的能量值，從而得到每個(gè)檢波器分別針對當(dāng)前產(chǎn)生地震波的勘探炮的炮檢對的能量值； (d)在所有勘探炮分別產(chǎn)生相同次數(shù)的地震波之后，基于多個(gè)檢波器分別針對每個(gè)勘探炮的炮檢對的能量值，繪制基于地質(zhì)模型的炮檢對可視性分析圖； (e)基于炮檢對可視性分析圖，確定在實(shí)際的地震勘探中應(yīng)該增加勘探炮數(shù)量的范圍， 其中，炮檢對的能量值是指一個(gè)檢波器接收到的一個(gè)地震波的能量值。
2.根據(jù)權(quán)利要求1所述的基于目標(biāo)層可視性分析的激發(fā)位置確定方法，其特征在于，在步驟(C)中，利用有限差分算子，對當(dāng)前產(chǎn)生的地震波進(jìn)行正向延拓到目標(biāo)反射層，然后再從目標(biāo)反射層反向延拓到檢波器，來模擬地震波被目標(biāo)反射層反射回來的反射波，從而記錄每個(gè)檢波器接收到的反射波的能量值。
3.根據(jù)權(quán)利要求1所述的基于目標(biāo)層可視性分析的激發(fā)位置確定方法，其特征在于，在步驟(c)中，當(dāng)產(chǎn) 生地震波時(shí)，地震波在頻率域中以下面的公式(I)的形式傳播，
(k2+ Δ ) U (rs, r, ω) =-f (ω) δ (r-rs)公式(I) 其中，k= ω /C(l (r)為波數(shù)，ω表示地震波的頻率，C0 (r)為位于目標(biāo)反射層的點(diǎn)r=(x, y, z)處的背景場速度 值，Λ為拉普拉斯算子，f(?)為勘探炮函數(shù)，δ (r-rs)是單位脈沖函數(shù)，rs表示勘探炮的三維坐標(biāo)。
4.根據(jù)權(quán)利要求3所述的基于目標(biāo)層可視性分析的激發(fā)位置確定方法，其特征在于，在步驟(C)中，利用有限差分算子，對當(dāng)前產(chǎn)生的地震波進(jìn)行正向延拓到目標(biāo)反射層，然后再從目標(biāo)反射層反向延拓到檢波器， 其中，檢波器接收到的反向延拓的反射波的地震波場通過公式(2)來表示，
P (rs，rg, ω)=
2k2 f v(r)G(r，rs; co)m(r)G(r，rg; co)d3r 公式(2) 其中，G(r, rs; ω)和G(r，rg; ω)分別是頻率域中勘探炮和檢波器位置的格林函數(shù)，&表示檢波器位置的三維坐標(biāo)，m(r)是位于目標(biāo)反射層的點(diǎn)r=(X，y，z)處速度模型擾動，V(r)為目標(biāo)反射層的所有點(diǎn)的集合。
5.根據(jù)權(quán)利要求4所述的基于目標(biāo)層可視性分析的激發(fā)位置確定方法，其特征在于，多個(gè)檢波器中的一個(gè)檢波器針對當(dāng)前產(chǎn)生地震波的勘探炮的炮檢對的能量值通過公式(3)來表示，
Psr (rs, rg, ω) = IP (rs, rg, ω) 公式(3 )。
【文檔編號】G01S5/20GK103645460SQ201310596834
【公開日】2014年3月19日 申請日期:2013年11月21日 優(yōu)先權(quán)日:2013年11月21日
【發(fā)明者】胡善政, 李亞林, 張曉斌, 何光明, 敬龍江, 陳燕雄, 童濤, 袁光銀, 蔡力, 管敏, 王曉陽 申請人:中國石油集團(tuán)川慶鉆探工程有限公司地球物理勘探公司