劉　杰,張忠濤,劉道理,王秀玲,徐　徽,屈　亮,袁　才

(中海石油(中國(guó))有限公司深圳分公司,廣東深圳518000)

?

強(qiáng)反射背景下沉積體邊界檢測(cè)及流體識(shí)別方法

劉杰,張忠濤,劉道理,王秀玲,徐徽,屈亮,袁才

(中海石油(中國(guó))有限公司深圳分公司,廣東深圳518000)

摘要:針對(duì)地震資料中存在強(qiáng)反射背景屏蔽有效儲(chǔ)層弱反射的情況,利用匹配追蹤算法的自適應(yīng)稀疏表示能力,提出了一種強(qiáng)反射背景下沉積體邊界檢測(cè)和流體識(shí)別的方法。該方法綜合利用測(cè)井資料、井旁道地震數(shù)據(jù)以及快速匹配追蹤算法實(shí)現(xiàn)了對(duì)強(qiáng)反射背景的最優(yōu)匹配和有效分離,削弱其對(duì)下伏地層造成的振幅“假亮點(diǎn)”以及砂體連續(xù)分布的假象;在此基礎(chǔ)上利用疊前彈性阻抗反演和流體因子識(shí)別技術(shù),精細(xì)落實(shí)砂體的沉積期次和展布規(guī)律,提高了巖性圈閉邊界和流體識(shí)別的精度。模型數(shù)據(jù)試算和實(shí)際地震資料處理的結(jié)果驗(yàn)證了該方法的可靠性和有效性。

關(guān)鍵詞：匹配追蹤;強(qiáng)反射背景;尖滅線;彈性阻抗反演;流體因子

強(qiáng)反射背景下有效信號(hào)的恢復(fù)以及優(yōu)質(zhì)儲(chǔ)層評(píng)價(jià)和含流體性質(zhì)檢測(cè),已成為國(guó)內(nèi)外不少學(xué)者研究的熱點(diǎn)和難點(diǎn)[1-4],如：煤系地層沉積時(shí),往往會(huì)造成煤層和目標(biāo)儲(chǔ)層存在上下疊置關(guān)系,從而導(dǎo)致有效砂體反射信息淹沒于煤層反射;地質(zhì)沉積間斷面會(huì)引起不整合破碎帶強(qiáng)反射，從而屏蔽不整合面附近的有效砂體反射;海侵體系域時(shí)期砂體頂部含灰質(zhì)成分或因泥巖脫水形成的高速強(qiáng)反射層對(duì)下伏地層反射的干擾等。針對(duì)上述特殊的地質(zhì)條件,如何剔除反射背景干擾,去偽存真,發(fā)展一種行之有效的地震數(shù)據(jù)目標(biāo)處理方法,進(jìn)而為優(yōu)質(zhì)儲(chǔ)層評(píng)價(jià)和流體識(shí)別奠定數(shù)據(jù)基礎(chǔ)是迫切需要解決的問題。

針對(duì)地震信號(hào)的稀疏表示方法——匹配追蹤(matching pursuit,MP)算法,國(guó)內(nèi)外學(xué)者進(jìn)行了大量研究,主要體現(xiàn)在非平穩(wěn)地震信號(hào)的時(shí)頻分析[5-9]、稀疏地震反演[10-13]以及地震稀疏去噪處理[14-16]等。MP算法可實(shí)現(xiàn)地震信號(hào)的自適應(yīng)分解,建立在超完備原子字典構(gòu)建的基礎(chǔ)上,將地震信號(hào)表示為所匹配原子的線性組合,從而利用有限個(gè)非零值組成的稀疏信號(hào)對(duì)連續(xù)信號(hào)進(jìn)行有效表征。因此利用MP地震信號(hào)稀疏表示的思想,結(jié)合測(cè)井資料進(jìn)行井旁道地震數(shù)據(jù)強(qiáng)反射干擾的時(shí)頻特征分析,在此基礎(chǔ)上通過強(qiáng)反射先驗(yàn)信息構(gòu)建超完備小波字典,最后實(shí)現(xiàn)強(qiáng)反射背景的匹配追蹤與識(shí)別。

本文針對(duì)地震資料存在上覆強(qiáng)反射背景屏蔽有效砂體反射的情況,提出了一種強(qiáng)反射背景下沉積體邊界檢測(cè)和流體識(shí)別方法,利用快速匹配追蹤算法,將強(qiáng)反射地震記錄從原始地震記錄中剝離,再通過建立不同的地質(zhì)、地球物理模型模擬不同地質(zhì)條件下反射特征的差異,對(duì)該方法的穩(wěn)定性和適用性進(jìn)行測(cè)試;最終結(jié)合已鉆井信息和實(shí)際地震資料的處理驗(yàn)證了該方法的有效性,為此類強(qiáng)反射背景下的儲(chǔ)層預(yù)測(cè)及流體識(shí)別提供了一種可靠的技術(shù)途徑。

1強(qiáng)反射背景的匹配追蹤與識(shí)別

1.1匹配追蹤算法原理

匹配追蹤算法可實(shí)現(xiàn)信號(hào)的自適應(yīng)分解,通過將信號(hào)在超完備時(shí)頻原子庫(字典D)中進(jìn)行分解,從而將地震信號(hào)表示為所匹配時(shí)頻原子的線性組合[16]。字典是一系列時(shí)頻原子的集合,通過對(duì)母小波進(jìn)行時(shí)間延遲、頻率調(diào)制、相位調(diào)制得到。經(jīng)典匹配追蹤算法的具體數(shù)學(xué)原理可以表示為：

(1)

(2)

式中：mγ(t)為子波字典中的母小波;mγ=(u,ω,φ)(t)為經(jīng)過時(shí)間u,頻率ω以及相位φ調(diào)制后的時(shí)頻原子;γ=(u,ω,φ)為小波字典的控制參數(shù)集合。

非平穩(wěn)地震信號(hào)時(shí)頻分析方法通常選擇高斯類時(shí)頻原子,但在實(shí)際地震強(qiáng)反射匹配過程中,創(chuàng)建超完備子波字典D使用的母小波mγ(t)與地震子波相關(guān)系數(shù)越大,則地震信號(hào)s的稀疏表示能力越強(qiáng),進(jìn)而可以更好地匹配出強(qiáng)反射地震背景。因此，通過井旁道強(qiáng)反射特征分析確定母小波mγ(t)的類型以及構(gòu)建小波字典控制集合Γ在地震信號(hào)稀疏重構(gòu)中起到了關(guān)鍵作用。原始地震數(shù)據(jù)s可表示為匹配原子的線性表達(dá)：

(3)

式中：a為稀疏表示系數(shù);N為匹配原子數(shù)目;dn為地震背景噪聲。

貪婪匹配追蹤算法是通過設(shè)置合理的迭代次數(shù)K或終止閾值σ,進(jìn)而估計(jì)信號(hào)s在字典D中的投影,從而實(shí)現(xiàn)對(duì)原信號(hào)的重構(gòu)。信號(hào)s可表達(dá)為沿mγ1方向的投影與垂直方向的分量之和：

(4)

式中：〈·〉代表向量?jī)?nèi)積運(yùn)算符;R1s表示初次自適應(yīng)分解后的殘差信號(hào)。若假設(shè)迭代算法已進(jìn)行N-1次,得到殘差信號(hào)為RNs,則第N次迭代算法為：

(5)

通過求解〈RN-1s,mγN〉的最大值，依次匹配出符合信號(hào)s結(jié)構(gòu)特征的時(shí)頻原子集合。

由于匹配追蹤算法每次匹配振幅最強(qiáng)、相關(guān)系數(shù)最優(yōu)的位置,考慮對(duì)時(shí)頻字典中時(shí)間參數(shù)u的約束,利用Hilbert變換求取復(fù)地震道信號(hào)并通過匹配振幅包絡(luò)最強(qiáng)點(diǎn)確定逐次迭代的時(shí)間中心點(diǎn)位置;此外,由于在匹配追蹤之前的原信號(hào)任意時(shí)刻的瞬時(shí)頻率和瞬時(shí)相位等先驗(yàn)信息皆給出了此時(shí)間點(diǎn)比較準(zhǔn)確的頻率和相位特征,從而進(jìn)一步縮小了頻率和相位的搜索半徑;最后動(dòng)態(tài)匹配追蹤算法的參數(shù)搜索方式可以表示為：

(6)

式中：A(t),φ(t)和ω(t)代表信號(hào)的瞬時(shí)振幅、瞬時(shí)相位以及瞬時(shí)頻率;u0為最大瞬時(shí)振幅包絡(luò)A(t)對(duì)應(yīng)的時(shí)間中心;U[ω(u0),δω]和U[φ(u0),δφ]表示以瞬時(shí)屬性為中心的頻率搜索鄰域和相位搜索鄰域;δω和δφ均為根據(jù)信號(hào)時(shí)頻特征定義的參數(shù)搜索半徑。此種動(dòng)態(tài)局部搜索相比全局搜索大大提高了計(jì)算效率,從而使匹配追蹤算法應(yīng)用于實(shí)際地震信號(hào)處理成為可能。

本文設(shè)計(jì)了由8個(gè)Ricker小波組成的理論信號(hào),如圖1a中黑線所示,圖中紅線表示無噪聲信號(hào)匹配重構(gòu)結(jié)果。分析可知,無噪聲情況下理論信號(hào)的匹配程度較高,基本無殘差信號(hào)。為了進(jìn)一步驗(yàn)證匹配追蹤算法的抗噪性,圖1b為含10%隨機(jī)噪聲情況下信號(hào)的匹配重構(gòu)結(jié)果,可以看出,即使在信號(hào)中含有噪聲的情況下,匹配追蹤算法仍然可以較好地重構(gòu)出理論信號(hào)。圖2為無噪聲信號(hào)匹配追蹤逐次迭代原子的結(jié)果,直觀地體現(xiàn)了8個(gè)Ricker子波的匹配過程。

圖1　理論信號(hào)匹配追蹤與重構(gòu)a 無噪聲信號(hào)匹配重構(gòu); b 含10%噪聲信號(hào)匹配重構(gòu)

圖2　無噪聲匹配追蹤逐次迭代結(jié)果

1.2強(qiáng)反射背景的匹配識(shí)別與分離

由于匹配追蹤算法具有自適應(yīng)的稀疏表示能力,并且具有較強(qiáng)的抗噪能力。因此可利用稀疏序列較好地匹配識(shí)別地震信號(hào),本文綜合利用測(cè)井資料、井旁道地震數(shù)據(jù)及快速匹配追蹤算法實(shí)現(xiàn)對(duì)強(qiáng)反射背景的最優(yōu)匹配和分離,從而使得淹沒在其中的弱小反射信息凸顯出來。基于快速匹配追蹤算法的強(qiáng)反射識(shí)別與分離算法的詳細(xì)流程如下：

1) 結(jié)合測(cè)井資料進(jìn)行井旁道地震數(shù)據(jù)強(qiáng)反射時(shí)頻特征分析,并提取相應(yīng)強(qiáng)反射層段地震子波;

2) 根據(jù)強(qiáng)反射干擾的時(shí)頻特征和地震子波提取結(jié)果,從子波字典的稀疏表示能力及地震子波擬合度兩方面來構(gòu)建超完備稀疏表示字典D={mγ(t)}γ∈Γ;

3) 輸入地震數(shù)據(jù)S,計(jì)算地震數(shù)據(jù)總道數(shù)Tn及采樣點(diǎn)數(shù)Sn,并求取其對(duì)應(yīng)的復(fù)數(shù)道;

4) 求取地震數(shù)據(jù)瞬時(shí)屬性,計(jì)算振幅包絡(luò)最大位置處的瞬時(shí)頻率、瞬時(shí)相位;

5) 檢索局部子波字典D={mγ(t)}γ∈Γ中的每一個(gè)原子,并計(jì)算子波相關(guān)系數(shù),選擇相關(guān)系數(shù)最大的原子為強(qiáng)反射識(shí)別原子,然后存儲(chǔ)匹配原子參數(shù)mγ=(u,ω,φ);

6) 利用阻尼最小二乘方法確定匹配原子對(duì)應(yīng)的振幅響應(yīng)信息,并獲取最終匹配的強(qiáng)反射干擾S0;

7) 強(qiáng)反射干擾匹配分離Snew=S-λ*S0,其中λ的選擇主要依據(jù)3個(gè)準(zhǔn)則：①同相軸極性不反轉(zhuǎn);②同相軸橫向連續(xù)性保持較好;③綜合不同地質(zhì)沉積模式確定相應(yīng)參數(shù)大小(注意λ參數(shù)的大小一般不超過1)。

由于地質(zhì)沉積模式不同引起的強(qiáng)反射干擾存在巨大差異,為了驗(yàn)證該方法的有效性和應(yīng)用性，本文從實(shí)際工區(qū)中總結(jié)出3種強(qiáng)反射背景模型并對(duì)其進(jìn)行了處理，理論模型及處理結(jié)果如圖3，圖4及圖5所示。

圖3a為含單一強(qiáng)反射界面理論模型,其主要對(duì)應(yīng)地質(zhì)沉積間斷或者因巖性突變引起的強(qiáng)反射界面,圖3b為理論反射系數(shù)模型,圖3c為30Hz零相位Ricker子波,圖3d為合成地震記錄,由圖可見，其中弱小反射信息幾乎被完全壓制,圖3e為利用強(qiáng)反射匹配追蹤方法識(shí)別出的強(qiáng)反射地震信息,可以看出匹配識(shí)別的結(jié)果與理論強(qiáng)反射干擾基本一致,圖3f為強(qiáng)反射分離后的地震記錄,可以觀察到被強(qiáng)反射掩蓋的弱小反射信息被凸現(xiàn)出來,并且同相軸與弱小反射界面吻合度較高。

圖4為含強(qiáng)反射層(即頂?shù)纂p強(qiáng)反射界面)理論模型及匹配追蹤結(jié)果,此種模型對(duì)應(yīng)著一套穩(wěn)定發(fā)育的高阻抗沉積層,由圖4e分析可知,高阻抗地層引起的強(qiáng)反射干擾可以被有效識(shí)別出來,雖然識(shí)別效果遜于單一強(qiáng)反射界面,但吻合度仍然較高。

圖3　含單一強(qiáng)反射界面理論模型及匹配追蹤結(jié)果a 理論縱波阻抗模型; b 理論反射系數(shù)模型; c 正演30Hz零相位Ricker子波; d 合成地震記錄; e 強(qiáng)反射匹配追蹤與識(shí)別結(jié)果,其中紅線為強(qiáng)反射識(shí)別結(jié)果,黑線為對(duì)應(yīng)地質(zhì)強(qiáng)反射層地震正演記錄; f 強(qiáng)反射分離后地震記錄

圖4　含強(qiáng)反射層(雙界面)理論模型及匹配追蹤結(jié)果a 理論縱波阻抗模型; b 理論反射系數(shù)模型; c 正演30Hz零相位Ricker子波; d 合成地震記錄; e 強(qiáng)反射匹配追蹤與識(shí)別結(jié)果,其中紅線為強(qiáng)反射識(shí)別結(jié)果,黑線為對(duì)應(yīng)地質(zhì)強(qiáng)反射層地震正演記錄; f 強(qiáng)反射分離后地震記錄

圖5為含漸變強(qiáng)反射層理論模型,此種模型與地質(zhì)上連續(xù)沉積背景且由于靠近砂巖的泥巖脫水形成的遞變強(qiáng)反射層對(duì)應(yīng),從圖5e和圖5f可見,對(duì)由漸變地層引起的強(qiáng)反射干擾的匹配識(shí)別效果仍然較好,弱小反射層的有效反射信息也得到了較好的恢復(fù)?？偨Y(jié)上述3種強(qiáng)反射背景模型處理結(jié)果可知,基于快速匹配追蹤算法的強(qiáng)背景反射識(shí)別與剝離技術(shù)具有較好的應(yīng)用效果和適用性。

為了進(jìn)一步驗(yàn)證強(qiáng)反射背景下匹配追蹤與識(shí)別方法的可行性和有效性,我們根據(jù)中國(guó)南部某探區(qū)的地質(zhì)沉積模式以及測(cè)井資料建立了二維縱波阻抗模型并用本文方法對(duì)其進(jìn)行分析。圖6為根據(jù)各期沉積體的含不同厚度泥巖隔層地震尖滅特征模型去強(qiáng)反射前后對(duì)比分析結(jié)果,圖7為薄砂互層不同巖性充填時(shí)尖滅點(diǎn)特征在去強(qiáng)反射前后的對(duì)比分析結(jié)果。其中圖6a,圖6d及圖6g和圖7a,圖7d及圖7g表示理論地質(zhì)模型,模型大小為橫向上CDP450,縱向上采樣時(shí)長(zhǎng)300ms,模型頂部為大套泥巖,泥巖下方發(fā)育因含灰或泥巖脫水形成的高速層,下覆兩期前積薄砂互層(上覆砂體尖滅點(diǎn)在CDP100位置處),模型中灰色表示上覆泥巖,黑色表示高速層,黃色表示前積砂巖,綠色表示下覆泥巖。

圖5　含漸變強(qiáng)反射層理論模型及匹配追蹤結(jié)果a 理論縱波阻抗模型; b 理論反射系數(shù)模型; c 正演30Hz零相位Ricker子波; d 合成地震記錄; e 強(qiáng)反射匹配追蹤與識(shí)別結(jié)果,其中紅線為強(qiáng)反射識(shí)別結(jié)果,黑線為對(duì)應(yīng)地質(zhì)強(qiáng)反射層地震正演記錄; f 強(qiáng)反射分離后地震記錄

圖6a,圖6d及圖6g分別為泥巖隔層厚度為3～13ms,2～8ms及1～3ms時(shí)的地質(zhì)模型(模型Ⅰ,模型Ⅱ,模型Ⅲ),砂體厚度變化范圍為1～12ms,高速層厚度為10～13ms。利用25Hz的零相位Ricker子波合成地震記錄如圖6b,圖6e及圖6h所示。分析可知,隨著薄砂互層的泥巖隔層厚度的逐漸減小,地震尖滅特征逐漸模糊(相位突變位置CDP140—CDP150—CDP165),尤其當(dāng)泥巖隔層厚度為1～3ms時(shí),地震剖面中無明顯尖滅現(xiàn)象。圖6c,圖6f及圖6i為強(qiáng)反射背景匹配分離后的地震記錄,圖中可見,本文方法消除了強(qiáng)反射干擾對(duì)下覆砂體有效反射信息的屏蔽作用,進(jìn)而地震尖滅特征更加明顯(相位突變位置分別為CDP120—CDP125—CDP135)。

圖6　含不同厚度泥巖隔層尖滅特征去強(qiáng)反射前后對(duì)比分析a理論地質(zhì)模型(模型Ⅰ); b 模型Ⅰ合成地震記錄; c 模型Ⅰ去強(qiáng)反射后地震記錄; d 理論地質(zhì)模型(模型Ⅱ); e 模型Ⅱ合成地震記錄; f 模型Ⅱ去強(qiáng)反射后地震記錄; g 理論地質(zhì)模型(模型Ⅲ);h 模型Ⅲ合成地震記錄; i 模型Ⅲ去強(qiáng)反射后地震記錄

圖7a,圖7d及圖7g為不同巖性充填時(shí)尖滅特征正演模型(模型Ⅳ,模型Ⅴ,模型Ⅵ),考慮兩套薄砂互層不同巖性充填時(shí)的地震尖滅特征識(shí)別精度,合成地震記錄如圖7b,圖7e及圖7h所示。隨著薄砂互層泥質(zhì)含量的不斷增加,砂體尖滅位置由清晰變得欠清晰(相位突變位置為CDP140—CDP140—CDP145)。圖7c,圖7f及圖7i為高速層強(qiáng)反射去除后的地震記錄,明顯可見,去強(qiáng)反射地震剖面中尖滅點(diǎn)位置向前移動(dòng)(相位突變位置為CDP115—CDP120—CDP125),相比強(qiáng)反射屏蔽有效反射時(shí)的尖滅點(diǎn)識(shí)別效果有所提高,進(jìn)而有助于落實(shí)砂體的平面展布,提高對(duì)巖性圈閉的識(shí)別精度。

圖7　不同巖性充填尖滅特征正演模型去強(qiáng)反射前后對(duì)比分析a 理論地質(zhì)模型(模型Ⅳ); b 模型Ⅳ合成地震記錄; c 模型Ⅳ去強(qiáng)反射后地震記錄; d 理論地質(zhì)模型(模型Ⅴ); e 模型Ⅴ合成地震記錄; f 模型Ⅴ去強(qiáng)反射后地震記錄; g 理論地質(zhì)模型(模型Ⅵ); h 模型Ⅵ合成地震記錄; i 模型Ⅵ去強(qiáng)反射后地震記錄

3實(shí)際資料處理

國(guó)內(nèi)某海上探區(qū)主要研究層段為新近系珠江組下段,由于受海侵期沉積環(huán)境變化的影響,造成了砂體頂部鈣質(zhì)含量的迅速增高,此外,由于下伏砂巖儲(chǔ)層孔滲條件好(屬中孔、中-高滲儲(chǔ)層),在越接近儲(chǔ)層的地方泥巖脫水越快,形成了一套穩(wěn)定的高速層,疊后地震剖面上主要表現(xiàn)為較強(qiáng)的反射同相軸淹沒了高速層下方有效儲(chǔ)層和流體響應(yīng)信息。圖8為去強(qiáng)反射前后地震剖面中砂體尖滅位置對(duì)比分析結(jié)果,其中圖8a和圖8b分別為L(zhǎng)ine1和Line2的去強(qiáng)反射前地震剖面,地震剖面中綠線為高速層頂界面(綠色箭頭指示紅色同相軸位置,高速層底界面即下伏砂巖頂界面)。強(qiáng)反射干擾導(dǎo)致了砂體連續(xù)分布的假象,黑色同相軸無明顯的相位突變,然而根據(jù)地質(zhì)沉積背景可知,此處含有兩期前積砂體,但砂體尖滅線位置難以確定。圖8c和圖8d分別為L(zhǎng)ine1和Line2去強(qiáng)反射后地震剖面,明顯看出,砂體頂?shù)捉缑婢霈F(xiàn)明顯的間斷點(diǎn),如圖中綠色箭頭指示位置,因此有助于劃分各期沉積砂體的平面展布范圍。

圖9為去除強(qiáng)反射前后儲(chǔ)層含氣性檢測(cè)差異性分析結(jié)果,其中圖9a和圖9b分別為去強(qiáng)反射前后小角度部分疊加地震數(shù)據(jù),分析可見,高速層強(qiáng)反射可以被有效的匹配分離出來,強(qiáng)反射附近的弱小反射信息顯現(xiàn)出來,并且去強(qiáng)反射后的地震數(shù)據(jù)同樣表現(xiàn)出較好的橫向連續(xù)性。

由于目標(biāo)儲(chǔ)層地震響應(yīng)受到強(qiáng)反射背景的壓制作用,地震剖面中出現(xiàn)的大量并非由于含氣引起的“假亮點(diǎn)”現(xiàn)象,圖9c和圖9d分別為去強(qiáng)反射前后疊前流體因子(Gassmann流體項(xiàng)f)直接反演結(jié)果,其中A井和B井對(duì)應(yīng)含氣井,C井對(duì)應(yīng)含水井,圖中測(cè)井解釋結(jié)果紅色表示含氣儲(chǔ)層,藍(lán)色表示含水儲(chǔ)層。明顯可見,原始地震流體因子反演剖面中C井出現(xiàn)含水“假亮點(diǎn)”,且流體因子出現(xiàn)含氣低值異常,與實(shí)鉆不符;而圖9d去強(qiáng)反射后的地震流體因子反演剖面可以看出,不僅A井和B井在含氣儲(chǔ)層位置吻合較好,而且在含水井C儲(chǔ)層發(fā)育位置反演結(jié)果與測(cè)井解釋結(jié)果也保持較高的一致性,可以有效識(shí)別儲(chǔ)層中的流體性質(zhì)。

圖8　地震去強(qiáng)反射前后砂體尖滅位置分析a Line1去強(qiáng)反射前地震剖面; b Line2去強(qiáng)反射前地震剖面; c Line1去強(qiáng)反射后地震剖面; d Line2去強(qiáng)反射后地震剖面

圖9　去強(qiáng)反射前后儲(chǔ)層含氣性檢測(cè)差異性對(duì)比a 小角度部分疊加地震數(shù)據(jù); b 強(qiáng)反射匹配分離后的小角度部分疊加地震數(shù)據(jù); c 原始流體因子反演結(jié)果; d 去強(qiáng)反射后流體因子反演結(jié)果

4結(jié)論

針對(duì)地震資料中存在上覆強(qiáng)反射背景屏蔽有效砂體反射的情況,本文發(fā)展了一種強(qiáng)反射背景下沉積體邊界檢測(cè)和流體識(shí)別方法。作為地震資料目標(biāo)處理環(huán)節(jié)之一,強(qiáng)反射背景匹配追蹤與分離方法將強(qiáng)反射地震記錄從原始地震記錄中剝離,進(jìn)而為沉積體邊界檢測(cè)和流體識(shí)別提供數(shù)據(jù)基礎(chǔ)。模型試算和實(shí)際資料處理表明：①利用匹配追蹤算法對(duì)強(qiáng)反射信息進(jìn)行識(shí)別和分離適用性強(qiáng),可廣泛應(yīng)用于多種地質(zhì)沉積模式;②該方法在強(qiáng)地震反射背景下有助于精確落實(shí)砂體的發(fā)育期次與疊置關(guān)系,能區(qū)分出由強(qiáng)反射干擾引起的目標(biāo)區(qū)儲(chǔ)層砂體連續(xù)分布的假象;③通過疊前流體因子反演能夠有效消除強(qiáng)反射對(duì)相鄰地層造成的非含氣振幅“假亮點(diǎn)”異常,進(jìn)一步降低了儲(chǔ)層流體檢測(cè)的多解性。

參考文獻(xiàn)

[1]秦成崗,施和生,張忠濤,等.珠江口盆地番禺低隆起-白云凹陷北坡SQ21.0層序陸架坡折帶沉積特征及油氣勘探潛力[J].中國(guó)海上油氣,2011,23(1):14-18

QIN C G,SHI H S,ZHANG Z T,et al.Sedimentary characteristics and hydrocarbon exploration potential along the SQ21.0 sequence shelf-break zone on Panyu low-uplift and the north slope of Baiyun sag,Pearl River Mouth Basin[J].China Offshore Oil and Gas,2011,23(1):14-18

[2]劉杰,秦成崗,張忠濤,等.低位楔三角洲砂體尖滅特征描述與識(shí)別技術(shù)——以珠江口盆地番禺地區(qū)為例[J].天然氣地球科學(xué),2013,24(6):1268-1273

LIU J,QIN C G,ZHANG Z T,et al.Sandbody pinchout description and identification technique of lowstand wedge-shaped delta：taking Panyu field in Pearl River Mouth Basin as an example[J].Natural Gas Geoscience,2013,24(6):1268-1273

[3]劉保國(guó).提高大牛地氣田三維地震儲(chǔ)層預(yù)測(cè)精度的思考[J].石油物探,2008,47(1):95-102

LIU B G.Considerations for improving accuracy of reservoir prediction from 3-D seismic data in Daniudi gas filed[J].Geophysical Prospecting for Petroleum,2008,47(1):95-102

[4]馬水平,王志宏,田永樂,等.地震反射同向軸的地質(zhì)含義[J].石油地球物理勘探,2010,45(3):454-457

MA S P,WANG Z H,TIAN Y L,et al.Studies on geological meaning of seismic reflection events[J].Oil Geophysical Prospecting,2010,45(3):454-457

[5]MALLAT S G,ZHANG Z F.Matching pursuit with time-frequency dictionaries[J].IEEE Transactions on Signal Processing,1993,41(12):3397-3415

[6]張顯文,韓立國(guó),王寧,等.地震信號(hào)譜分解匹配追蹤快速算法及其應(yīng)用[J].石油物探,2010,49(1):1-6

ZHANG X W,HAN L G,WANG N,et al.Seismic spectral decomposition fast matching pursuit algorithm and its application[J].Geophysical Prospecting for Petroleum,2010,49(1):1-6

[7]韓海英,王志章,王宗俊,等.基于Ricker類地震子波的匹配追蹤[J].石油物探,2014,53(1):17-25

HAN H Y,WANG Z Z,WANG Z J,et al.Matching pursuit based on Ricker-like seismic wavelet[J].Geophysical Prospecting for Petroleum,2014,53(1):17-25

[8]CASTAGNA J P,SUN S J,SIEGFRIED R W.Instantaneous spectral analysis:Detection of low-frequency shadows associated with hydrocarbons[J].The Leading Edge,2003,22(2):120-127

[9]李傳輝,張繁昌.地震信號(hào)可變分辨率匹配追蹤頻譜成像方法[J].石油物探,2012,51(3):213-218

LI C H,ZHANG F C.Variable resolution matching pursuit spectrum imaging of seismic signals[J].Geophysical Prospecting for Petroleum,2012,51(3):213-218

[10]張繁昌,李傳輝,印興耀.三角洲砂巖尖滅線的地震匹配追蹤瞬時(shí)譜識(shí)別方法[J].石油地球物理勘探,2012,47(1):82-88

ZHANG F C,LI C H,YIN X Y.Delta fringe line recognition based on seismic matching pursuit instantaneous spectral characteristics[J].Oil Geophysical Prospecting,2012,47(1):82-88

[11]張繁昌,李傳輝.地震信號(hào)復(fù)數(shù)域高效匹配追蹤分解[J].石油地球物理勘探,2013,48(2):171-175

ZHANG F C,LI C H.Complex domain efficient matching pursuit decomposition of seismic signals[J].Oil Geophysical Prospecting,2013,48(2):171-175

[12]NGUYEN T H.High resolution seismic reflectivity inversion[D].Houston:University of Houston,2008

[13]LIU X J,YIN X Y,WU G C.Pre-stack seismic inversion based on orthogonal matching pursuit algorithm[J].Expanded Abstracts of 84thAnnual Internat SEG Mtg,2014:576-579

[14]趙天姿,宋煒,王尚旭.基于匹配追蹤算法的時(shí)頻濾波去噪方法[J].石油物探,2008,47(4):367-371

ZHAO T Z,SONG W,WANG S X.Time-frequency filtering de-noise method based on matching pursuit algorithm[J].Geophysical Prospecting for Petroleum,2008,47(04):367-371

[15]張汛汛,張繁昌,劉漢卿.基于快速匹配追蹤算法的地震道集剩余時(shí)差校正[J].石油物探,2015,54(4):420-426

ZHANG X X,ZHANG F C,LIU H Q.Seismic gathers residual moveout correction based on fast matching pursuit algorithm[J].Geophysical Prospecting for Petroleum,2015,54(4):420-426

[16]李傳輝.地震信號(hào)快速匹配追蹤分解及瞬時(shí)譜分析[D].山東青島:中國(guó)石油大學(xué)(華東),2012

LI C H.Fast matching pursuit decomposition and instantaneous spectral analysis of seismic signals[D].Qingdao:China University of Petroleum(Hua dong),2012

(編輯：朱文杰)

Sediment boundary identification and fluid detection for the seismic data with strong background reflections

LIU Jie,ZHANG Zhongtao,LIU Daoli,WANG Xiuling,XU Hui,QU Liang,YUAN Cai

(CNOOCChinaLIMITED-Shenzhen,Guangzhou510240,China)

Abstract:Considering the ability of self-adaptive sparse representation for matching pursuit algorithm,we proposed a method for sedimentary boundary identification and fluid detection on the seismic data with strong background reflections shielding the effective information in weak reflections.Integrated of log data,uphole seismic trace and matching pursuit algorithm,the strong background reflections could be recognized and peeled away from the original seismic data;meanwhile,some interference of “false bright spots” and artificial continuous distribution of sand bodies could be weakened easily.Then,the elastic impedance inversion and fluid factor detection could be implemented to finely characterize the distribution of sand bodies and improve the accuracy of lithological trap detection and fluid identification.Applications on theoretical model and field seismic data demonstrate the reliability and effectiveness of our method.

Keywords:matching pursuit,strong background reflection,pinch-out boundary,elastic impedance inversion,fluid factor

文章編號(hào)：1000-1441(2016)01-0142-08

DOI:10.3969/j.issn.1000-1441.2016.01.018

中圖分類號(hào)：P631

文獻(xiàn)標(biāo)識(shí)碼：A

基金項(xiàng)目：國(guó)家科技重大專項(xiàng)(2011ZX05023-002-007)項(xiàng)目資助。

作者簡(jiǎn)介：劉杰(1982—),男,工程師,碩士,主要從事地震綜合解釋與儲(chǔ)層預(yù)測(cè)等研究工作。

收稿日期：2015-07-03;改回日期：2015-08-25。

劉杰,張忠濤,劉道理,等.強(qiáng)反射背景下沉積體邊界檢測(cè)及流體識(shí)別方法[J].石油物探,2016,55(1):-149

LIU Jie,ZHANG Zhongtao,LIU Daoli,et al.Sediment boundary identification and fluid detection for the seismic data with strong background reflections[J].Geophysical Prospecting for Petroleum,2016,55(1):-149

This research is financially supported by the National Science and Technology Major Project of China (Grant No.2011ZX05023-002-007).