杜斌山， 賀振華， 王緒本 ， 雍學(xué)善, 劉應(yīng)如

(1.成都理工大學(xué) 地球物理學(xué)院，成都 610059；2.中國石油勘探開發(fā)研究院 西北分院，蘭州 730020)

基于地震偶極子波多重積分的初始阻抗模型建立方法

杜斌山1，2， 賀振華1， 王緒本1， 雍學(xué)善2, 劉應(yīng)如2

(1.成都理工大學(xué) 地球物理學(xué)院，成都 610059；2.中國石油勘探開發(fā)研究院 西北分院，蘭州 730020)

地震阻抗反演需要適當(dāng)?shù)牡皖l儲(chǔ)層模型，常規(guī)反演方法是基于地震、測井資料和地震速度預(yù)測來建立儲(chǔ)層低頻模型。由于儲(chǔ)層縱橫向變化快，常規(guī)方法易出現(xiàn)低頻模型的多解性和精度低等問題。綜合考慮儲(chǔ)層頂?shù)椎卣鸱瓷浣缑孀杩购蛢?chǔ)層厚度特征，提出了整體反映這些特征的地震偶極子波概念和相應(yīng)的地震響應(yīng)特征分析方法。以此為基礎(chǔ)，形成了偶極子波多重積分約束的地震反演初始模型建立新方法。該方法充分利用地震多重積分對儲(chǔ)層位置和厚度識別與表示功能，再結(jié)合測井聲波曲線的分頻譜等數(shù)據(jù)構(gòu)建地震初始阻抗低頻模型。該建模方法在柴達(dá)木盆地某工區(qū)實(shí)際應(yīng)用結(jié)果表明，所建模型與實(shí)際復(fù)雜的風(fēng)化殼裂縫儲(chǔ)層地震剖面吻合好，較完整地反映了地震剖面上的地質(zhì)層位與構(gòu)造特征，該模型對提高地震反演的穩(wěn)定性和精度十分有利。

地震偶極子波； 多重積分； 頻譜特征； 儲(chǔ)層厚度； 初始阻抗模型

0 引言

在油氣儲(chǔ)層預(yù)測中，地震反演技術(shù)應(yīng)用較為廣泛。其類型有遞推反演、測井約束地震反演、線性和非線性反演等多種，地震反演的小波變換與遺傳算法和神經(jīng)網(wǎng)絡(luò)算法都是非線性的。雍學(xué)善[1]等提出了逐道遺傳外推的反演方法；孟憲軍等[2]提出地震反演的精度與地震資料的品質(zhì)、層位和斷層解釋的精度以及地層沉積模式與地層接觸關(guān)系等都有緊密相關(guān)；崔巖等[3]指出初始模型建立對提高儲(chǔ)層預(yù)測的效果非常重要?，F(xiàn)階段儲(chǔ)層地震波阻抗反演存在的主要問題是常規(guī)波阻抗反演預(yù)測儲(chǔ)層時(shí)只簡單運(yùn)用地震速度和測井聲波資料建立低頻初始模型，多解性較強(qiáng)。

實(shí)際上，對砂泥巖薄互層、碳酸鹽巖、火山巖等特殊儲(chǔ)層來說，地震預(yù)測存在分辨率低、多解性強(qiáng)兩大關(guān)鍵問題[4-5]。對非均質(zhì)儲(chǔ)層預(yù)測精度和薄儲(chǔ)層的分辨率方面要求更高，但隨著預(yù)測分辨率的提高，各種預(yù)測方法以及迭代算法等都有較強(qiáng)的要求，會(huì)加入更多的人為信息，測井資料處理、井震標(biāo)定和層位解釋的精度都會(huì)影響反演結(jié)果，造成反演的多解性。為了克服地震預(yù)測的不確定性及多解性，世界上各大軟件研發(fā)公司和學(xué)者相繼對頻譜分解、地震屬性分析、小波分析、疊前分頻AVO等方法技術(shù)進(jìn)行了認(rèn)真研究[6-12]，但效果并不太理想，還是具有較強(qiáng)的多解性。

筆者提出了一種基于偶極地震子波多重積分約束的儲(chǔ)層初始阻抗模型建立的方法(DWMIC)。該方法把儲(chǔ)層頂、底雙反射系數(shù)，儲(chǔ)層時(shí)間厚度及其偶合波形作為一個(gè)基本完整研究單元，再結(jié)合測井和地震資料的頻譜特征，構(gòu)建儲(chǔ)層初始阻抗模型與地震波形、能量和多重積分變化的模型，在地震波形相關(guān)分析與多重積分約束下實(shí)現(xiàn)地震數(shù)據(jù)初始阻抗模型建立，最終形成一種適用于非均質(zhì)儲(chǔ)層的偶極子波多重積分約束的地震反演初始阻抗模型建立新方法。

1 模型建立方法原理

1.1 偶極子與地震偶極子波

關(guān)于偶極子的概念[13]，應(yīng)用地球物理百科詞典的作者Robert. E. Sheriff定義偶極子(Dipole)為“距離無限接近、電量相等的一對電荷或者是距離無限接近、極性相反的一對電極”。按此定義，在數(shù)學(xué)上，可將偶極子d(x) 表示為式(1)。

d(x)=q(x)+a·q(x-b)

a=±1,b→0

(1)

式中：q(x)為電荷或者極性場；a=1對應(yīng)關(guān)于電荷的偶極子，a=-1對應(yīng)關(guān)于電極的偶極子；b為兩極之間的距離，對于應(yīng)用地球物理勘探的實(shí)際問題，b→0是不可能存在的，在地震記錄中，由于較薄儲(chǔ)層頂、底之間的時(shí)間距離(一般為數(shù)十毫秒)遠(yuǎn)小于地震記錄長度L(一般為數(shù)千毫秒)，可以近似假設(shè)距離b→0而不會(huì)帶來明顯誤差，從而認(rèn)為它們是偶極子。我們把對應(yīng)于薄層(地層)頂、底的一對反射系數(shù)視為偶極子，定義其地震響應(yīng)為地震偶極子波，其數(shù)學(xué)表示式如式(2)所示。

W(t,b)=W(t)+f·W(t-b)

f∈R,b→0

(2)

式(2)為地震偶極子波與反射系數(shù)褶積合成地震記錄的計(jì)算模型。其中，S(t)為地震記錄(合成地震記錄)，t為時(shí)間，b為地層頂?shù)追瓷湎禂?shù)之間的時(shí)間差，W(t,b)表示為地震偶極子波函數(shù)，R(t,b)為儲(chǔ)層反射系數(shù)，R為實(shí)數(shù)。f=1，等價(jià)于前面公式a=1的情況，f=-1，等價(jià)于R(t,b)=Rt(t)的情況。但是f不再局限于正負(fù)1的特定值，而是任意實(shí)數(shù)。因此，可稱偶極子波為擴(kuò)展偶極子。擴(kuò)展的目的是為了使薄層(或地層)的上下接觸層的巖性任意變化，以適應(yīng)各種復(fù)雜的接觸關(guān)系，并簡稱擴(kuò)展偶極子波為偶極子波，R(t,b)為儲(chǔ)層反射系數(shù)。當(dāng)儲(chǔ)層厚度b=λ/4時(shí)，達(dá)到諧振。儲(chǔ)層反射系數(shù)為式(3)。

R(t,b)=Rt(t)

(3)

式(4)為單個(gè)地震子波表達(dá)式[14]。

W(t)=Ae-β(f*(t-τ))2sin(2πfi)

(4)

式中：w表示為地震子波；i為時(shí)間域函數(shù)；A為子波最大振幅；f為地震子波的主頻；β為衰減系數(shù)；τ為延遲時(shí)。

針對單個(gè)儲(chǔ)層頂、底反射系數(shù)的極性特征，可分為4類見圖1。

圖1 儲(chǔ)層頂?shù)鬃杩沟乃姆N組合模式及其對應(yīng)地震響應(yīng)特征(偶極子波)Fig.1 Four kinds of reservoir top and bottom interface impedance models and the corresponding seismic responses(seismic dipole wavelet )(a)儲(chǔ)層頂?shù)鬃杩沟乃姆N組合模式； (b)對應(yīng)地震響應(yīng)特征

圖1(a)為四個(gè)儲(chǔ)層模型結(jié)構(gòu)，其中(1)、(2)為頂?shù)追瓷湎禂?shù)極性相反的組合形式，(3)、(4)為頂?shù)追瓷湎禂?shù)極性相同組合形式。圖1(b)為模型對應(yīng)的地震偶極子波。地震偶極子波的波形有時(shí)比較簡單(凸、 凹型)，有時(shí)相對復(fù)雜(上、下臺(tái)階型)。

基于偶極地震子波的儲(chǔ)層初始阻抗模型建立方法的特點(diǎn):①把儲(chǔ)層頂、底界面雙反射系數(shù)特征、儲(chǔ)層厚度等信息整體融入地震響應(yīng)特征(偶極子波)中。而偶極地震子波的精確構(gòu)建首先需要通過測井、地震資料和井-震標(biāo)定等技術(shù)獲得儲(chǔ)層頂、底界面的反射系數(shù)或者波阻抗信息以及地層厚度信息;②偶極地震子波包含的時(shí)間厚度信息以及層位信息，可通過對含有偶極子波地震道的積分或者多重積分來預(yù)測或者表示;③由于采用了偶極子波和多重積分方法，DWMIC有利于低頻建模和消除隨機(jī)噪音，從而提高反演的多解性和穩(wěn)定性。

1.2 基于地震多重積分的阻抗數(shù)據(jù)轉(zhuǎn)換

通常地震道積分是通過地震反射系數(shù)的積分求和，實(shí)現(xiàn)地震數(shù)據(jù)轉(zhuǎn)換為地震相對阻抗剖面，是波阻抗反演的重要方法之一，且應(yīng)用廣泛[15]。地震多重積分是指對(合成)地震記錄進(jìn)行多輪次的地震道振幅累加運(yùn)算，即從起始t0時(shí)間開始不斷累加地震振幅的數(shù)值，在時(shí)間t進(jìn)行記錄，得到一個(gè)新的數(shù)據(jù)體，在t0～t區(qū)間累加地震振幅的次數(shù)即為積分的重?cái)?shù)。多重積分的表達(dá)式如下：

sm(t)=∫m∫s(t)dt

(5)

k=(N-1)，(N-2),…2,1,0

(6)

式(5)為偶極地震子波多重積分表達(dá)式。其中，S(t)為地震記錄，m為積分重?cái)?shù),sm(t)為積分結(jié)果。

在多重積分的離散式(6)中，N為數(shù)據(jù)總樣點(diǎn)數(shù)，n為當(dāng)前位置計(jì)算時(shí)間樣點(diǎn)數(shù)，Δt為采樣間隔，k為正整數(shù)。在式(6)中，規(guī)定s0(nΔt)為原始地震記錄或者子波。

圖2為一頂、底反射系數(shù)相等但極性相反的地震偶極子波的五重積分，圖2(a)為偶極子波，圖2(b)為其五重積分圖。五重積分的主瓣能量突出，能較好地反映地層的實(shí)際位置和地層厚度(兩紅線之間)，體現(xiàn)了多重積分的基本功能。需要特別注意的是多重積分重?cái)?shù)的選擇和保證獲得良好積分效果的數(shù)據(jù)處理方法：①在利用式(6)進(jìn)行的多重積分計(jì)算時(shí),其偶次積分結(jié)果中包含隨時(shí)間增加而增大的異常低的低頻分量，必須通過平滑濾波，加以消除，然后進(jìn)行相應(yīng)的奇數(shù)次積分；②建模過程中，我們常采用奇數(shù)次積分，這是因?yàn)橐恢胤e分(即通常的道積分)是對地震振幅數(shù)據(jù)進(jìn)行累加運(yùn)算，反映地震相對阻抗特征，且存在90°相位轉(zhuǎn)換， 使其主瓣能量與儲(chǔ)層所在位置對應(yīng)，便于地震地質(zhì)解釋[16]。二重積分，是在一重積分的基礎(chǔ)上對地震響應(yīng)數(shù)據(jù)進(jìn)行積分運(yùn)算，存在兩個(gè)(或者偶數(shù)次)90°相位轉(zhuǎn)換。因此，偶數(shù)重積分反映的是反射界面特征，而經(jīng)過三個(gè)(奇數(shù)個(gè))90°的相位轉(zhuǎn)換后，其結(jié)果又具有相對阻抗特征，但符號相反。所以奇數(shù)重積分反映儲(chǔ)層相對阻抗特性，這正是阻抗反演和阻抗建模所需要的。

圖2 偶極子波地震記錄和偶極子波5重積分Fig.2 Seismic dipole wavelet (left) and its 5th multiple integral (right)(a)偶極子波；(b)五重積分圖

1.3 地震多重積分與地層位置和厚度的對應(yīng)關(guān)系

為了明確多重積分與地層位置和厚度的對應(yīng)關(guān)系，我們制作一個(gè)儲(chǔ)層厚度不同，隔層厚度變化的水平分層地質(zhì)模型。從上到下儲(chǔ)層時(shí)間厚度(ms)分別為：90ms、80ms、70ms、60ms、50ms、40ms、30ms、20ms、10ms、5ms、2.5ms,見圖3(d)。對該模型用50HzRicker子波制作合成地震記錄，1/4波長相當(dāng)于10ms時(shí)間長度。對合成地震記錄分別進(jìn)行一、三、五重次積分的結(jié)果示于圖3(a)、圖3(b)和圖3(c)。若從能否識別頂、底界面的分辨能力來看，一重積分結(jié)果能辨識10ms薄層，相當(dāng)于1/4波長；三重積分結(jié)果能辨識30ms儲(chǔ)層，相當(dāng)于3/4波長儲(chǔ)層厚度；五重積分也能辨識70ms儲(chǔ)層，相當(dāng)于7/4波長儲(chǔ)層，但對儲(chǔ)層位置和厚度的表示和識別而言，則積分重?cái)?shù)多的效果更好，因此可用多重積分結(jié)果直接進(jìn)行相對地震阻抗建模。

圖3 多重積分與儲(chǔ)層頂?shù)追瓷湮恢煤皖A(yù)測儲(chǔ)層厚度圖Fig.3 The seismic trace multiple integral and resrvoiv's top and bottom reflection interface impedances and predict thickness(a)合成記錄的一重積分; (b)三重積分; (c)五重積分剖面;(d)儲(chǔ)、隔層厚度模型

1.4 地震多重積分的頻譜特征

為認(rèn)識地震數(shù)據(jù)多重積分相對阻抗的頻譜特征。我們需要對積分結(jié)果進(jìn)行頻譜分析。圖4是同一合成地震記錄的多重積分頻譜。由圖4可見，隨著積分重?cái)?shù)的增加，出現(xiàn)兩個(gè)明顯變化，頻帶愈來愈窄，主頻愈來愈低，最低頻率達(dá)到3Hz，說明高積分重?cái)?shù)的結(jié)果對構(gòu)建低頻模型有利。

2 初始阻抗模型建立的基本步驟

基于偶極子波和多重積分進(jìn)行初始阻抗模型構(gòu)建包括以下關(guān)鍵步驟：

1)測井基礎(chǔ)數(shù)據(jù)分析與整理。首先從單井測井資料解釋中獲取準(zhǔn)確反映儲(chǔ)層地震響應(yīng)特征的有關(guān)參數(shù)，然后在充分利用地震、測井和錄井資料的基礎(chǔ)上，結(jié)合鉆井、測井資料提取儲(chǔ)層厚度信息以及儲(chǔ)層上下介質(zhì)的聲波、密度和波阻抗等儲(chǔ)層信息。采用大層位標(biāo)定、分頻標(biāo)定、合成記錄精細(xì)標(biāo)定等技術(shù)手段，提取最優(yōu)地震子波，同時(shí)利用地震地質(zhì)分層、地震解釋層位以及波形特征等多種信息獲得最佳的時(shí)深關(guān)系。

2)基于地震多重積分的相對阻抗預(yù)測。獲取井點(diǎn)處儲(chǔ)層的偶極地震子波響應(yīng)特征函數(shù)，對其進(jìn)行地震多重積分運(yùn)算。對多個(gè)井點(diǎn)處地震多重積分結(jié)果進(jìn)行比較，最終確定適合該區(qū)的最佳滑動(dòng)時(shí)窗長度和積分重?cái)?shù)等計(jì)算參數(shù)。一般來說，積分計(jì)算時(shí)窗越大，低頻成分越多；計(jì)算時(shí)窗越小，高頻成分較豐富。 利用所確定的最佳計(jì)算參數(shù)和滑動(dòng)時(shí)窗長度對實(shí)際工區(qū)的地震數(shù)據(jù)體進(jìn)行地震多重積分運(yùn)算，以生成相對阻抗預(yù)測數(shù)據(jù)體。

圖4 同一地震數(shù)據(jù)的多重?cái)?shù)積分頻譜特征Fig.4 Seismic multiple integral frequency spectrum characteristics of the same seismic data(a)一重積分頻譜；(b)三重積分頻譜；(c)五重積分頻譜； (d)七重積分頻譜；(e)九重積分頻譜；(f)十一重積分頻譜

3)測井阻抗數(shù)據(jù)的分頻段濾波。對測井阻抗曲線進(jìn)行不同帶寬的頻率濾波，獲得分頻測井阻抗曲線。由于該曲線與地震道不同次數(shù)的重積分結(jié)果有相似性，可為井-震資料的有效結(jié)合提供依據(jù)。

4)多重積分與地震速度初始阻抗模型建立。地震數(shù)據(jù)多重積分轉(zhuǎn)換的相對阻抗頻譜特征具有明顯分段性(見圖4和圖5)，通過結(jié)合地震連續(xù)速度譜預(yù)測的頻譜，將其疊加(合并)，可逐步逼近測井資料的寬帶阻抗頻譜，并保持原始地震頻譜中低頻信息的完整性，為低頻建模做好準(zhǔn)備。在圖5中，d1井井旁地震道速度譜的頻譜帶寬為0Hz～3Hz，井旁地震道五重積分譜的帶寬為3Hz～7Hz，地震道三重積分具有7Hz～14Hz的頻譜特性，地震道一重積分具有14Hz～48Hz的頻譜特性，一重積分與地震道頻譜(12Hz～4 8Hz)有一定的相似性，于是形成了一個(gè)具有完整的連續(xù)的低頻特征的頻譜系列，有效凸顯3Hz～40Hz的隱蔽信息。

3 應(yīng)用實(shí)例分析

為檢驗(yàn)上述方法的效果，下面結(jié)合柴達(dá)木盆地某工區(qū)地震偶極子波多重積分阻抗建模的例子，說明建模的具體方法與效果。根據(jù)該區(qū)d1井等鉆井情況看，有多口井在儲(chǔ)層段獲得工業(yè)氣流，儲(chǔ)層類型為基巖風(fēng)化殼裂縫型。測井曲線表現(xiàn)出高聲波時(shí)差儲(chǔ)層特征[17-18]，d1井有多段儲(chǔ)層，其中目的層段壓裂試氣后，獲高產(chǎn)天然氣。在該區(qū)建模有以下步驟：

1)需要進(jìn)行測井-地震資料特征分析(如對d1關(guān)鍵井進(jìn)行精細(xì)井震標(biāo)定，利用井震資料聯(lián)合提取子波，并得到最佳的時(shí)-深關(guān)系和最優(yōu)子波，以獲得高質(zhì)量合成記錄)。對井旁地震剖面進(jìn)行多重積分，并進(jìn)行譜分析，獲得地震頻譜基本特征。圖6顯示了上述資料的綜合特征。為下一步研究提供基礎(chǔ)。

2)構(gòu)建不同頻譜特征的儲(chǔ)層模型。如圖7所示，測井曲線儲(chǔ)層段的聲波頻譜特性及地震速度譜的頻譜特征之間有一定的相似性，但是地震頻譜缺乏低頻段信息，測井反射系數(shù)的頻譜也缺乏低頻分量，而地震速度頻譜的低頻成分豐富，它們之間又有一定的互補(bǔ)性。將這些特征相結(jié)合，對建立具有連續(xù)頻譜特性的模型十分重要。

3)進(jìn)行合成地震記錄多重積分與測井分頻阻抗的計(jì)算，確定它們之間的相關(guān)性(圖8)。除個(gè)別道外，總體上多重積分相對阻抗結(jié)果與不同頻段的分頻濾波結(jié)果具有相似性。考慮到多重積分結(jié)果在頻率域具有分段性，可以將多個(gè)多重積分的頻譜合并，如圖9(b)所示。同理，也可將測井阻抗分頻譜合并，得到圖9(c)所示的結(jié)果。圖9(b)、圖9(c)在有效頻段較為相似，地震多重積分合并數(shù)據(jù)頻譜與測井分頻合并的相對阻抗數(shù)據(jù)較為相似，說明地震多重積分相對阻抗預(yù)測有較好的效果。

圖5 偶極子波地震多重分頻譜、地震速度頻譜特征對比圖Fig.5 The comparisons among seismic data spectrum, multiple integral and seismic velocity spectrum(a)測井?dāng)?shù)據(jù)頻譜；(b)地震速度譜的頻譜；(c)五重積分； (d)三重積分；(e)一重積分；(f)地震數(shù)據(jù)頻譜

圖6 d1井合成記錄標(biāo)定圖及多重積分剖面圖Fig.6 The well d1 synthetic seismogram calibration and multiple integral sections

圖7 測井阻抗-地震速度數(shù)據(jù)頻譜特征分析Fig.7 The comparisons among well logging impedance spectrum, seismic velocity frequency spectrum and seismic data spectrum(a)測井波阻抗頻譜(0.1ms)；(b)地震數(shù)據(jù)頻譜；(c)測井反射系數(shù)頻譜(0.1ms)； (d)地震速度頻譜；(e)地震速度轉(zhuǎn)換層速度

圖8 d1井合成地震記錄多重積分與測井分頻濾波結(jié)果對比圖Fig.8 The comparisons of seismic multiple integral traces and well logging fractional filter traces

4)建立相對阻抗初始模型。由于偶極地震子波和儲(chǔ)層地震響應(yīng)的多重積分均包括地層厚度和層位信息。因此基于地震偶極子波多重積分的初始阻抗模型可以直接在經(jīng)過前述處理的地震剖面上建立，而無需輸入常規(guī)地震層位和地層厚度。圖10為一重積分相對阻抗剖面圖。圖11為多重積分合并的相對阻抗剖面。多重積分剖面合并后，低頻特征較為清晰，剖面背景更清楚，能夠很好地反映目的層基巖風(fēng)化殼地質(zhì)特征。

圖9 測井阻抗頻譜與合并相對阻抗數(shù)據(jù)頻譜對比分析圖Fig.9 The comparisons among well logging impedance spectrum, seismic multiple integral spectrum and well logging fractional filter trace spectrum(a)原始測井阻抗數(shù)據(jù)頻譜圖；(b)多重積分?jǐn)?shù)據(jù)合并相對阻抗頻譜圖； (c)測井分頻阻抗數(shù)據(jù)合并相對阻抗頻譜圖

圖10 地震一重積分相對阻抗數(shù)據(jù)剖面模型Fig.10 The seismic 1st integral relative impedance section model

圖11 地震一、三、五、七重積分合并相對阻抗剖面模型Fig.11 The seismic relative impedance section model by stacking the 1st,3rd ,5th and 7th multiple integral profiles

圖12為地震一、三、五、七重積分合并相對阻抗數(shù)據(jù)與地震速度合并剖面圖，加入地震速度頻譜的低頻成分后，初始阻抗剖面模型能較全面地反映基巖頂面特征。由于該初始模型與實(shí)際地震剖面十分接近，完全不是常規(guī)線條型模型樣式，因此，有助于降低最終的地震參數(shù)反演的多解性，提高反演的穩(wěn)定性。

圖12 地震一、三、五、七重積分合并相對阻抗數(shù)據(jù)與地震速度合并的初始阻抗模型Fig.12 The stacking relative impedance section model with seismic 1st, 3rd ,5th and 7thmultiple integral profiles and seismic velocity analysis

4 結(jié)論

1)筆者提出的基于地震偶極子波多重積分的初始阻抗模型建立方法，充分挖掘了地震記錄蘊(yùn)含地震地質(zhì)信息，綜合考慮儲(chǔ)層偶極子波的響應(yīng)特征、偶極子波雙反射系數(shù)整體波形特征，提高了儲(chǔ)層預(yù)測可靠性，有利于降低反演的多解性。

2)明確了地震多重積分儲(chǔ)層的含義，構(gòu)建了地震偶極子波多重積分?jǐn)?shù)學(xué)表達(dá)式，明確了儲(chǔ)層頂?shù)捉缑娣瓷湎禂?shù)估算以及多重積分相對阻抗與地層界面的奇偶特性，證實(shí)地震奇數(shù)重積分具有相對阻抗預(yù)測的效果。

3)地震數(shù)據(jù)多重積分相對阻抗的頻譜特性具有明顯分段性，結(jié)合地震資料速度預(yù)測，采用多重積分合并預(yù)測相對阻抗，逐步逼近儲(chǔ)層阻抗頻譜，提高低頻模型預(yù)測效果，提高了初始模型建立的精度。

[1] 雍學(xué)善, 余建平, 石蘭亭.一種三維高精度儲(chǔ)層參數(shù)反演方法[J].石油地球物理勘探,1997，32(6):852-856.YONGXS,YUJP，SHILT.Anaeeuratemethodfor3-DreservoirParameterinversion[J].OGP,1997,32(6):852-856.(InChinese)

[2] 孟憲軍，金翔龍，鈕學(xué)民,等.地震反演中的三維復(fù)雜約束模型[J].石油大學(xué)學(xué)報(bào)(自然科學(xué)版),2004, 28(6):21-26.MENGXJ,JINXL,NIUXM.Researchof3-Dcomplexconstrainedmodelinseismicinversion[J].JournaloftheUniversityofPetroleum，2014, 28(6):21-26.(InChinese)

[3] 崔巖, 王彥飛, 楊長春. 帶先驗(yàn)知識的波阻抗反演正則化方法研究[J].地球物理學(xué)報(bào),2009，52(8)：2135-2141.CUIY,WANGYF,YANGCC.Regu1arizingmethodwithaprioriknow1edgeforseismicimpedanceinversion[J]．ChineseJ.Geophys，2009,52(8): 2135-2141.(InChinese)

[4] 雍學(xué)善. 砂泥巖薄互層儲(chǔ)層預(yù)測的難點(diǎn)與對策[J]. 天然氣工業(yè)，2005，25(增B)：96-100.YONGXS.Difficultiesandcountermeasurestopredictreservoirswiththinsand-shalealternatinglayers[J].NaturalGasIndustry，2005,25(SB):96-100.(InChinese)

[5]WIDESSM.B.Howthinisathinbed[J].Geophysics, 38(6):1176-1180.

[6]PIERRETHORE,OLIVIERDUPLANTIER，YUNQUIXU. 2Seismicinversionuncertainty:whatdoesreallymatter[C].2013SEGHoustonAnnualMeeting.

[7]CHARLESI.PURYEAR,JOHNP.CASTAGNA.Layer-thicknessdeterminationandstratigraphicinterpretationusingspectralinversion:Theoryandapplication[J].Geophysics，2008，73(2):37-48.

[8]WANGX,ZHANGYQ.Pre-stackinversioncombinedwithgeostatisticalsimulationtopredictthinreservoirs[C]. 2010SEGDenverAnnualMeeting，2010.

[9]WHITCOMBEDN.Elasticimpedancenormalization[J].Geophysics,2002,67(1):63-67.

[10]劉曉晶, 印興耀, 吳國忱,等.基于正交匹配追蹤算法的疊前地震反演方法[J].石油地球物理勘探，2015,50(5):925-935.LIUXJ，YINXY，WUGC，etal.ZongZhaoyun;Prestackseismicinversionbasedonorthogonalmatchingpursuitalgorithm[J].OilGeophysicalProspecting,2015,50(5):925-935.(InChinese)

[11]馬勁風(fēng). 地震勘探中廣義彈性阻抗的正反演[J]. 地球物理學(xué)報(bào),2003, 46(1)：118-124.MAJF.Forwardmodelingandinversionmethodofgeneralizedelasticimpedanceinseismicexploration[J].ChineseJournalofGeophysic. 2003, 46(1)：118-124.(InChinese)

[12]CONNOLLY.Elasticimpedance[J].TheLeadingEdge,1999, 4:438-452.

[13]ROBERT.E.SHERIFF.Encyclopedicdictionaryofappliedgeophysics,FourthEdition[C].SEG,TulsaOklahoma,USA,2006.

[14]雍學(xué)善, 吳勝和.儲(chǔ)集層厚度譜的建立及其意義[J].新疆石油地質(zhì),2005,26(6):647-649.YONGXS,WUSH.EstablishmentandSignificanceofReservoirThicknessSpectrum[J].XinjangPetroleumGeology,2005,26(6):647-649.(InChinese)

[15]陸基孟. 地震勘探原理[M]. 北京:石油大學(xué)出版社，1996.LUJM.Principleofseismicexploration[M].Beijing:petroleumuniversitypress,1996.(InChinese)

[16]HONGLIUZENG.Seismicanalysisofverythinbeds:whichattributetouse[C]. 2012SEGAnnualMeeting. (phaseandfrequencyslices)

[17]賀振華, 鄧英爾, 劉樹根，等.巖石彈性參數(shù)對滲流測試分析的影響[J]. 天然氣工業(yè)，2006，26(6):44-46.HEZH,DENGYE,LIUSG,etal.Influenceofrockelasticparametersontestanalysesofflowinrock[J].NaturalGasIndustry,2006,26(6):44-46.(InChinese)

[18]賀振華, 黃德濟(jì), 文曉濤.裂縫油氣藏地球物理預(yù)測[M]. 成都:四川科學(xué)技術(shù)出版社，2007.HEZH,HUANGDJ,WENXT.Fracturedreservoirgeophysicalprediction[M].Chengdu:Sichuanscienceandtechnologypress,2007. (InChinese)

Initial impedance model establishment based on the seismic dipole wavelet and multiple integral methods

DU Binshan1,2, HE Zhenhua1, WANG Xuben1, YONG Xueshan2, LIU Yingru2

(1.Chengdu university of technology Geophysical Institute，Chengdu 610059， China；2.Chinese petroleum exploration and development research institute ，northwest branch,Nanzhou 730020， China)

The seismic impedance inversion requries a proper low-frequency reservoir model. The conventional inversion method is based on logging data and seismic velocity prediction to establish reservoir low-frequency model. Due to the reservoir parameters horizontally and vertically quickly change, the uncertainty and low precision of the seismic inversion by using conventional low frequency model are not avoidable. In order to overcome these problems, a dipole wavelet multiple integral constraint method (DWMIC) for seismic inversion initial model establishment is proposed in this paper. The seismic dipole wavelet of a reservoir is a comprehensive representation of reservoir's top and bottom reflection interface impedances and thickness and the seismic trace multiple integral combined well logging acoustic fractional spectrum data can obtain reservoir positions and seismic impedances. The DWMIC modelling method has been applied to a work area in Qaidam basin. The result shows that the model constructed by the DWMIC fits the real seismic profile quite well for a complex weathering crust fracture reservoir. It means the model can be fully reflected the geological horizon and structure characteristics on the seismic section, suggesting that it is possible to improve the stability and precision of seismic inversion.

seismic dipole wavelet; multiple integral; spectrum characteristics; reservoir thickness; initial impedance model

2016-11-02 改回日期：2016-12-21

國家自然科學(xué)基金(41374111);國家重大專項(xiàng)(2016zx05007-006);柴達(dá)木重大專項(xiàng)(2016E-03)

杜斌山(1970-)，男，博士，主要從事地震反演、儲(chǔ)層預(yù)測、油藏描述等方面研究,E-mail：dubs@petrochina.com.cn。

1001-1749(2017)01-0071-10

P 631.4

A

10.3969/j.issn.1001-1749.2017.01.11