本發(fā)明涉及油田勘探開發(fā)技術(shù)領(lǐng)域，特別是涉及到一種基于深變子波的合成記錄制作方法。

背景技術(shù)：
：

疊前深度偏移技術(shù)的研究是近10多年來全球油氣地球物理勘探領(lǐng)域的熱點，是目前地震資料處理技術(shù)中歸位相對準(zhǔn)確的一種偏移方法，特別是在速度復(fù)雜或者構(gòu)造復(fù)雜地區(qū)具有較明顯的優(yōu)勢，主要用來解決鹽丘、復(fù)雜基底、復(fù)雜斷塊、生物礁、泥底辟、流體底辟等復(fù)雜地質(zhì)體的波場成像問題。隨著偏移算法的成熟、計算機(jī)能力的進(jìn)步，規(guī)模化的疊前深度偏移成為可能。深度域地震剖面與時間域地震剖面相比，具有地下構(gòu)造真實、直觀、便于解釋等特點，并且無需進(jìn)行時深轉(zhuǎn)換就能得到初始深度構(gòu)造。目前，疊前深度偏移已經(jīng)是工業(yè)界的標(biāo)準(zhǔn)處理流程。因此如何對疊前深度偏移后的地震數(shù)據(jù)進(jìn)行后續(xù)的井震標(biāo)定、構(gòu)造解釋及其它的應(yīng)用就成為必須面對的問題，并且，制作深度域合成記錄是對疊前深度偏移后的地震數(shù)據(jù)進(jìn)行一切后續(xù)工作的基礎(chǔ)和關(guān)鍵。

張雪建(2000)等認(rèn)為在時間域合成地震記錄是地震子波與反射系數(shù)的褶積，深度域的合成地震記錄也與時間域類似，只不過先要將時間域的子波轉(zhuǎn)換成深度域的子波；然后利用聲波時差曲線和密度曲線計算出深度域的反射系數(shù),二者褶積便可得到深度域的合成地震記錄。該方法的關(guān)鍵在于提取子波和計算反射系數(shù)，要在深度偏移剖面上選取品質(zhì)較好的井旁地震道計算子波的主頻,要選用零相位子波，在計算反射系數(shù)之前，要對聲波時差曲線和密度曲線進(jìn)行校正、去野值等處理。

何惺華(2005)對深度域地震記錄的應(yīng)用開展了討論，從時間域和深度域的地震波場描述出發(fā)，對深度域的子波、褶積和Fourier變換等基本問題進(jìn)行分析探討。證明了時間域與深度域中的地球物理概念、原理和方法是相互聯(lián)系的，其數(shù)學(xué)表達(dá)方式是相似的。對深時轉(zhuǎn)換中采樣間隔的選擇、頻率譜與波數(shù)譜數(shù)據(jù)的關(guān)系和深時轉(zhuǎn)換剖面上視頻率變化等實際問題，給出了合理的解釋。

王永剛(2005)對深度域的地震資料合成記錄制作方法進(jìn)行了探討。總結(jié)歸納深度域合成記錄制作過程應(yīng)該包括：(1)利用聲波、密度等測井曲線進(jìn)行層位劃分建立深時關(guān)系，根據(jù)多層介質(zhì)反射系數(shù)的波動方程推導(dǎo)公式，求取頻率域的反射系數(shù)；(2)對時間域的井旁道記錄利用同態(tài)反褶積方法提取地震子波，得到頻率域的子波譜；(3)利用頻率域褶積模型計算合成記錄，作反傅里葉變換得到時間域合成地震記錄，再根據(jù)時深關(guān)系求出每個時域采樣點對應(yīng)的深度點，按照所要求的深度采樣點間隔進(jìn)行重采樣，得到深度域合成地震記錄。

林伯香等(2006)認(rèn)為，諧波在時間域中經(jīng)過反射系數(shù)“系統(tǒng)”后不改變頻率,符合線性時不變系統(tǒng)的條件,因此在時間域中可以用子波與反射系數(shù)序列褶積的方法計算時間域合成地震記錄。由于在深度域中子波的波數(shù)隨介質(zhì)速度變化,因此地震子波在深度域中傳播并不滿足線性時不變系統(tǒng)的條件,不能直接用褶積的方法計算深度域合成地震記錄。林伯香等利用簡單的理論模型,通過對深度域中的速度、深度參數(shù)進(jìn)行適當(dāng)?shù)淖儞Q,使其在新深度域中符合線性時不變系統(tǒng)的條件,在保證層間旅行時不變的情況下,定義深度域地震子波的介質(zhì)速度為標(biāo)準(zhǔn)速度，將其他各層的速度都以標(biāo)準(zhǔn)速度進(jìn)行調(diào)整,并對層間厚度也做相應(yīng)的調(diào)整。通過上述變換,可把時間域中的計算方法引用到深度域中,再通過反變換處理,將新深度域中的褶積結(jié)果反變換到原深度域中。如果不考慮深度域與時間域之間的重采樣誤差,常規(guī)時間域合成地震記錄經(jīng)時深轉(zhuǎn)換到深度域后的結(jié)果應(yīng)該與深度域合成地震記錄結(jié)果一致。

以上方法都沒有考慮到當(dāng)密度曲線欠缺時，如何根據(jù)實際情況由速度曲線轉(zhuǎn)換得出密度曲線，并且當(dāng)欠缺密度曲線時，一味地使用Gardner公式來轉(zhuǎn)換，準(zhǔn)確率太低。另外，在深度域，子波在某一時刻的傳播長度受地層速度的影響，子波的長度隨著速度的變化而變化，是地層速度的函數(shù)，不是固定值，這就給深度域合成記錄的制作帶來困難。

為此我們發(fā)明了一種新的基于深變子波的合成記錄制作方法，解決了以上技術(shù)問題。

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

本發(fā)明的目的是提供一種可以根據(jù)具體的沉積相態(tài)得到更合理的密度曲線，從而大大提高了深度域合成記錄制作的準(zhǔn)確率的基于深變子波的合成記錄制作方法。

本發(fā)明的目的可通過如下技術(shù)措施來實現(xiàn)：基于深變子波的合成記錄制作方法，該基于深變子波的合成記錄制作方法包括：步驟1，輸入時間域子波、速度曲線和密度曲線，在沒有密度曲線時，不輸入密度曲線；步驟2，建立時間域子波與速度的函數(shù)關(guān)系，將時間域子波轉(zhuǎn)換成深度域子波，即生成深變子波；步驟3，在步驟1輸入了密度曲線時，使用速度和密度曲線，計算反射系數(shù)；步驟4，將深變子波與反射系數(shù)進(jìn)行褶積，制作深度域合成記錄。

本發(fā)明的目的還可通過如下技術(shù)措施來實現(xiàn)：

在步驟2中，設(shè)定w(t)是時間域雷克子波，w(h)是深度域子波，v(h)是子波在深度h時的傳播速度，計算深度為h時，深度域子波w(h)的形態(tài)，

假設(shè)，t時刻的時間域子波w(t)與深度h的深度域子波w(h)相對應(yīng)，t^/時刻的時間域子波w(t^/)相對應(yīng)的深度域子波是w(h^/)，

若將t與t^/之間的時間平均分成n份，那么時間間隔Δt＝(t^/-t)/n，

假設(shè)，時刻t₁＝t-Δt，t₂＝t-2Δt，與t₁時刻對應(yīng)的深度是h₁，與t₂時刻對應(yīng)的深度是h₂，那么

如公式1和公式2所述，對于時間域子波w(t)周期T內(nèi)的任意時刻t^/都可計算出所對應(yīng)的深度值h^/，從而找準(zhǔn)任意時刻的w(t^/)所對應(yīng)的深度域子波w(h^/)，并且由公式3計算出由于速度變化而時刻改變的深度域子波周期H(v)，

將時間域子波轉(zhuǎn)換到深度域后，由于速度的影響，波形隨深度發(fā)生壓縮或拉伸兩種變化。

在步驟3中，在步驟1沒有輸入密度曲線時，根據(jù)陸相沉積環(huán)境當(dāng)中不同沉積相態(tài)，給定針對性的擬合公式，計算密度曲線。

在步驟3中，在步驟1沒有輸入密度曲線，只有速度曲線時，首先根據(jù)實際工區(qū)資料，擬合出比較符合實際的密度曲線，采用的密度與速度的經(jīng)驗公式是

D＝αv_p^β (公式4)

其中D是密度，v_p是縱波速度，α、β為參數(shù)值，根據(jù)實際工區(qū)，針對α、β做出相應(yīng)的調(diào)整。

在步驟3中，當(dāng)有真實的速度曲線和密度曲線時，由如下公式得到每一深度的的反射系數(shù)：

其中R_i是第i層與第i+1層界面的深度域反射系數(shù)，ρ_i是第i層的密度，v_i是第i層的速度。

在步驟4中，在時間域，合成記錄是子波與反射系數(shù)的褶積：

$f\left(t\right)=w\left(t\right)*R\left(t\right)=\underset{0}{\overset{T}{\∫}}w\left(\mathrm{\τ}\right)R(t-\mathrm{\τ})d\mathrm{\τ}$

其中，w(t)和R(t)分別是時間域子波和反射系數(shù)，T是時間域子波周期，現(xiàn)在，將時間域子波w(t)轉(zhuǎn)換成深度域子波w(h)，深度域子波周期H(v)由步驟2計算得到，由速度和密度曲線得出深度域反射系數(shù)R(h)，用τ代表時間延遲，基于深變子波的合成記錄公式如下：

本發(fā)明以子波、速度、密度為輸入，建立子波與速度的準(zhǔn)確的函數(shù)關(guān)系，即深變子波，進(jìn)而正確計算子波在每一時刻的長度；當(dāng)欠缺密度曲線時，可以根據(jù)陸相沉積環(huán)境當(dāng)中不同沉積相態(tài)，給定針對性的擬合公式，計算出更真實的密度曲線，有利于反射系數(shù)的準(zhǔn)確合理計算，然后，將深度子波與反射系數(shù)進(jìn)行褶積，得到深度域的合成記錄結(jié)果。該裝置解決了深度域子波長度隨速度變化的難題，并且可以根據(jù)具體的沉積相態(tài)得到更合理的密度曲線，從而大大提高了深度域合成記錄制作的準(zhǔn)確率。

附圖說明：

圖1為基于深變子波的合成記錄制作裝置示意圖；

圖2為時間域子波轉(zhuǎn)換成深度域子波示意圖；

圖3為不同參數(shù)的Gardner公式擬合密度曲線對比圖；

圖4為時間域和深度域合成記錄結(jié)果對比圖。

具體實施方式：

為使本發(fā)明的上述和其他目的、特征和優(yōu)點能更明顯易懂，下文特舉出較佳實施例，并配合附圖所示，作詳細(xì)說明如下。

如圖1所示，圖1為本發(fā)明的基于深變子波的合成記錄制作方法的流程圖。

在步驟101中，收集時間域子波、速度曲線和密度曲線，并且進(jìn)行輸入。如果沒有密度曲線，可以在后面的步驟104中生成。

在步驟102中，通過建立時間域子波與速度的函數(shù)關(guān)系，將時間域子波轉(zhuǎn)換成深度域子波，即生成深變子波。

設(shè)定w(t)是時間域雷克子波，w(h)是深度域子波，v是速度曲線，計算深度為h時，深度域子波w(h)的形態(tài)。

假設(shè)，t時刻的時間域子波w(t)與深度h的深度域子波w(h)相對應(yīng)，t^/時刻的時間域子波w(t^/)相對應(yīng)的深度域子波是w(h^/)。

若將t與t^/之間的時間平均分成n份，那么時間間隔Δt＝(t^/-t)/n。

假設(shè)，時刻t₁＝t-Δt，t₂＝t-2Δt，與t₁時刻對應(yīng)的深度是h₁，與t₂時刻對應(yīng)的深度是h₂，那么

如公式1和公式2所述，對于時間域子波w(t)周期T內(nèi)的任意時刻t^/都可以很方便地計算出所對應(yīng)的深度值h^/，從而找準(zhǔn)任意時刻的w(t^/)所對應(yīng)的深度域子波w(h^/)，并且由公式3可以計算出由于速度變化而時刻改變的深度域子波周期H(v)。

將時間域子波轉(zhuǎn)換到深度域后，由于速度的影響，波形隨深度發(fā)生壓縮或拉伸兩種變化，如附圖2所示。圖2中，時間域子波是零相位雷克子波，通過本發(fā)明將其轉(zhuǎn)換成圖2中的深度域子波，可以看出，雷克子波在時間域是“時不變”子波，轉(zhuǎn)換到深度域后為“深變”子波，由于速度的影響，波形隨深度發(fā)生壓縮或拉伸兩種變化。

在步驟103中，判斷步驟101是否輸入密度曲線，如果沒有密度曲線，則進(jìn)入到步驟104；如果有密度曲線，則進(jìn)入到步驟105。

在步驟104中，根據(jù)陸相沉積環(huán)境當(dāng)中不同沉積相態(tài)，給定針對性的擬合公式，計算密度曲線。

當(dāng)只有速度曲線，欠缺密度曲線時，必須首先根據(jù)實際工區(qū)資料，擬合出比較符合實際的密度曲線。Gardner在1974年給出的密度與速度的經(jīng)驗公式是

D＝αv_p^β (公式4)

其中D是密度，v_p是縱波速度，α＝0.31，β＝0.25。該公式直到今天仍被廣泛使用，但是它并不是放之四海皆準(zhǔn)的公式，本專利會根據(jù)實際工區(qū)，針對α、β做出相應(yīng)的調(diào)整。當(dāng)α、β取不同的參數(shù)值時，計算出的密度曲線存在差異，如附圖3所示。該圖中，(a)是速度曲線；(b)是α＝0.31，β＝0.25時，計算的密度曲線；(c)是α＝0.25，β＝0.27時，計算的密度曲線。可以看出，當(dāng)α、β取不同值時，計算出來的密度曲線存在差異，需要根據(jù)實際工區(qū)的情況，對α、β這兩個參數(shù)進(jìn)行調(diào)整。

在步驟105中，根據(jù)實際工區(qū)資料，使用速度和密度曲線，得到準(zhǔn)確的反射系數(shù)。

當(dāng)有真實的速度曲線和密度曲線時，由如下公式可以得到每一深度的的反射系數(shù)：

其中R_i是第i層與第i+1層界面的深度域反射系數(shù)，ρ_i是第i層的密度，v_i是第i層的速度。

在步驟106中，將深變子波與反射系數(shù)進(jìn)行褶積，制作深度域合成記錄。

在時間域，合成記錄是子波與反射系數(shù)的褶積：

$f\left(t\right)=w\left(t\right)*R\left(t\right)=\underset{0}{\overset{T}{\∫}}w\left(\mathrm{\τ}\right)R(t-\mathrm{\τ})d\mathrm{\τ}$

其中，w(t)和R(t)分別是時間域子波和反射系數(shù)，T是時間域子波周期?，F(xiàn)在，將時間域子波w(t)轉(zhuǎn)換成深度域子波w(h)，深度域子波周期H(v)由步驟2計算得到，由速度和密度曲線得出深度域反射系數(shù)R(h)，用τ代表時間延遲?；谏钭冏硬ǖ暮铣捎涗浌饺缦拢?/p>

因為基于深變子波的合成記錄考慮了地層速度的影響，所以能準(zhǔn)確刻劃出地層界面的真實深度，并且由于受到地層速度的影響，深度域合成記錄反射軸的胖瘦與時間域不同，如附圖4所示。該圖中，上圖是時間域合成記錄結(jié)果，下圖是深度域合成記錄結(jié)果。從兩種合成記錄結(jié)果的對比來看，深度域合成記錄結(jié)果能刻劃出地層界面的真實深度，并且由于受到地層速度的影響，深度域合成記錄反射軸的胖瘦與時間域不同。

該基于深變子波的合成記錄制作方法解決了深度域子波長度隨速度變化的難題，并且可以根據(jù)具體的沉積相態(tài)得到更合理的密度曲線，從而大大提高了深度域合成記錄制作的準(zhǔn)確率。