顧 雯，章 雄，徐 敏，梁 虹，張洞君，羅 晶，鄭 虹

（中國(guó)石油集團(tuán)川慶鉆探工程有限公司地球物理勘探公司，四川成都610213）

強(qiáng)屏蔽下薄儲(chǔ)層高精度預(yù)測(cè)研究
——以松遼盆地三肇凹陷為例

顧 雯，章 雄，徐 敏，梁 虹，張洞君，羅 晶，鄭 虹

（中國(guó)石油集團(tuán)川慶鉆探工程有限公司地球物理勘探公司，四川成都610213）

松遼盆地三肇凹陷扶余油層以淺水三角洲沉積體系為主，單層砂體厚度薄，橫向變化快，受T2界面強(qiáng)反射層屏蔽作用影響，目的層砂體無明顯地震反射特征。如何降低強(qiáng)屏蔽層的調(diào)諧影響，突出薄砂體并進(jìn)行精細(xì)刻畫，一直是困擾本區(qū)勘探開發(fā)的技術(shù)難題。研究發(fā)現(xiàn)強(qiáng)屏蔽現(xiàn)象主要與反射系數(shù)的干涉作用有關(guān)，地震波是反射系數(shù)干涉作用形成的，當(dāng)?shù)貙又写嬖谳^大反射系數(shù)時(shí)，由于地震波主頻較低，會(huì)突出強(qiáng)阻抗界面的影響，弱反射系數(shù)被淹沒，解決強(qiáng)屏蔽效應(yīng)問題主要應(yīng)減少地震波干涉效應(yīng)的影響。由于提高分辨率處理手段局限性較大，經(jīng)研究分析及試驗(yàn)，采用地震波形指示反演技術(shù)，即利用地震波形的橫向變化信息代替?zhèn)鹘y(tǒng)變差函數(shù)，實(shí)現(xiàn)高頻成分的井震聯(lián)合反演，降低地震強(qiáng)阻抗界面的影響，有效解決了研究區(qū)強(qiáng)屏蔽下的薄儲(chǔ)層預(yù)測(cè)問題。

扶余油層；強(qiáng)屏蔽；干涉效應(yīng)；變差函數(shù)；地震波形指示反演

近年來，隨著松遼盆地北部整體勘探程度的提高，油氣勘探主要目標(biāo)逐步從構(gòu)造油氣藏轉(zhuǎn)向巖性油氣藏或構(gòu)造 巖性復(fù)合油氣藏。薄儲(chǔ)層勘探面對(duì)的問題依然是地震分辨率問題［1－3］，即當(dāng)?shù)卣鸱直媛市∮?／4波長(zhǎng)時(shí)，在地震剖面上難以分辨出薄儲(chǔ)層，使得測(cè)井與地震數(shù)據(jù)標(biāo)定困難。正演分析表明，強(qiáng)屏蔽主要是地震分辨率較低，強(qiáng)反射系數(shù)界面在調(diào)諧過程中“淹沒”薄層弱反射系數(shù)形成的。為了克服地震數(shù)據(jù)成像分辨率不足的問題，LEVY等［4］提出了稀疏脈沖地震反演方法，但該反演方法受干涉作用影響，地震信息與薄儲(chǔ)層沒有直接對(duì)應(yīng)關(guān)系，無法對(duì)薄層橫向變化做出有效預(yù)測(cè)［5－8］。由于儲(chǔ)層井間變化大，井間薄儲(chǔ)層仍存在很大不確定性。已有的各種疊后反演、疊前反演以及振幅、頻率、相位、波形等地震屬性分析技術(shù)雖然在儲(chǔ)層預(yù)測(cè)和油氣識(shí)別中取得了一定效果，但是在薄儲(chǔ)層研究方面的應(yīng)用還存在各自的局限性。疊后波阻抗反演是一種比較成熟的技術(shù)，但其反演結(jié)果的分辨能力依賴于反射波主頻高低和頻帶寬度。同樣地，振幅、頻率、相位、波形等地震屬性的分辨率也受到反射波主頻和帶寬的限制，不能滿足薄儲(chǔ)層預(yù)測(cè)的要求［910］。HAAS等［11］于20世紀(jì)90年代引入了地質(zhì)統(tǒng)計(jì)學(xué)反演方法，該方法基于空間域樣點(diǎn)分布，對(duì)井均勻分布的要求較高。此外，變差函數(shù)的統(tǒng)計(jì)尤其是變程的確定往往不能精細(xì)反映儲(chǔ)層空間沉積相的變化，導(dǎo)致模擬結(jié)果平面地質(zhì)規(guī)律性差、隨機(jī)性強(qiáng)［12－15］。

松遼盆地三肇凹陷扶余油層以淺水三角洲沉積體系為主，單砂體厚度大多在2．0m左右，非均質(zhì)性強(qiáng)、橫向變化較大，受上覆T2界面強(qiáng)反射干涉影響，儲(chǔ)層無明顯地震反射特征。強(qiáng)屏蔽下薄互層砂體的精細(xì)刻畫成為制約該地區(qū)有利巖性圈閉識(shí)別以及水平井部署的主要瓶頸，需要建立一套針對(duì)薄互層儲(chǔ)層及河道砂體的儲(chǔ)層預(yù)測(cè)方法。本文采用井震聯(lián)合的地震波形指示反演方法有效拓寬地震頻帶，降低強(qiáng)阻抗界面的屏蔽影響，有效解決了研究區(qū)強(qiáng)屏蔽薄儲(chǔ)層預(yù)測(cè)難題。

1 強(qiáng)屏蔽薄砂體預(yù)測(cè)難點(diǎn)分析

三肇凹陷肇35井區(qū)F11油層組屬于致密薄砂巖儲(chǔ)層，由于青一段泥巖與泉四段存在較大的波阻抗差異，在上覆T2界面形成了強(qiáng)反射，受其干涉作用影響，F(xiàn)11油層表現(xiàn)為T2下伏橫向不連續(xù)的弱波峰反射，地層內(nèi)部為單一波谷反射，儲(chǔ)層無明顯地震反射特征（圖1）。

圖1 過井地震剖面

由測(cè)井聲波曲線分析可知，青山口組地震層速度為2 500m／s，泉四段地震層速度為3 280m／s，上下兩層速度差異較大，形成了一個(gè)較強(qiáng)的波阻抗界面。由于地震分辨率低，泉四段內(nèi)部地層形成的若干較小的反射與該界面發(fā)生干涉作用時(shí)，主要體現(xiàn)了較大反射系數(shù)的影響，小的反射系數(shù)作用不明顯。由此可見，強(qiáng)阻抗界面對(duì)下伏地層的屏蔽作用主要是由存在較大反射系數(shù)時(shí)地震主頻低、分辨率不夠造成的。如果能夠提高地震分辨率，就可以有效避免強(qiáng)阻抗界面蔽屏作用的影響。

進(jìn)一步通過正演模擬來分析強(qiáng)阻抗界面的屏蔽作用。首先制作Zhao43井合成地震記錄，子波采用標(biāo)準(zhǔn)零相位雷克子波，主頻為40Hz，T2反射界面表現(xiàn)為強(qiáng)波峰反射（圖2）；然后采用一系列不同頻率的子波進(jìn)行正演實(shí)驗(yàn)，子波主頻從40Hz逐步提高到100Hz（圖3），隨著分辨率的提高，強(qiáng)反射界面的屏蔽作用逐漸減小，下伏弱反射層逐漸顯現(xiàn)，說明強(qiáng)屏蔽現(xiàn)象主要是由層間干涉作用造成的。如果信號(hào)頻帶足夠?qū)?，各反射系?shù)不發(fā)生干涉，就可以有效解決強(qiáng)阻抗界面的屏蔽問題。因此，解決強(qiáng)屏蔽問題的關(guān)鍵在于拓寬頻帶，增加高頻成分。

圖2 Zhao43井合成地震記錄

圖3 強(qiáng)屏蔽作用與地震分辨率關(guān)系分析

提高信號(hào)分辨率一般通過地震資料高分辨率處理來實(shí)現(xiàn)，但受采集手段的制約，處理方法僅能在一定程度上提高分辨率，提高幅度一般局限在15%左右，很難獲取50Hz以上的高頻成分。本文采用井震聯(lián)合反演的方法，利用井資料補(bǔ)充高頻成分，去除強(qiáng)屏蔽作用的影響，實(shí)現(xiàn)強(qiáng)屏蔽下的薄儲(chǔ)層精細(xì)預(yù)測(cè)。

2 地震波形指示反演

近年來，井震聯(lián)合的反演技術(shù)取得了迅速發(fā)展，其大體可以分為兩類。一類是以地震為主的經(jīng)典反演方法，例如稀疏脈沖等確定性反演方法。這種反演方法從地震資料出發(fā)，利用相關(guān)公式計(jì)算阻抗，不能有效突破地震分辨率的限制，無法獲取薄層信息。另一類是以井模擬為主的廣義反演方法，例如地質(zhì)統(tǒng)計(jì)學(xué)隨機(jī)反演?；趨f(xié)克里金、序貫高斯、模擬退火等地質(zhì)統(tǒng)計(jì)學(xué)算法的隨機(jī)反演，其本質(zhì)是利用儲(chǔ)層參數(shù)的空間分布特征實(shí)現(xiàn)測(cè)井參數(shù)模擬，獲得一組等概率的儲(chǔ)層參數(shù)模型。

地震波形指示反演［16－18］與地質(zhì)統(tǒng)計(jì)學(xué)反演一樣，利用井的信息補(bǔ)充高頻成分。早期這種模擬方法一般是用井信息插值模型（克里金、自然鄰域等算法），后期人們采用地質(zhì)統(tǒng)計(jì)學(xué)進(jìn)行隨機(jī)模擬。由于統(tǒng)計(jì)學(xué)模擬的方法基于變差函數(shù)表征儲(chǔ)層空間結(jié)構(gòu)特征，模擬時(shí)樣本井的優(yōu)選參照變程控制，高頻信息完全來自于井，沒有考慮地震信息，地震信息只是起到對(duì)隨機(jī)模擬結(jié)果進(jìn)行優(yōu)選的作用?；隈薹e理論的正演實(shí)踐表明，地震波的干涉特征與反射系數(shù)結(jié)構(gòu)和分布具有密切的關(guān)系，反射系數(shù)的垂向分布，包括反射系數(shù)的間距、大小、個(gè)數(shù)，決定了地震波形干涉樣式。在相似的地質(zhì)條件下，反射波波形與反射系數(shù)結(jié)構(gòu)能形成良好的匹配關(guān)系，可以根據(jù)波形的干涉特征優(yōu)化反射系數(shù)的分布。因此，地震波形指示反演可以看作廣義的反演過程，是在傳統(tǒng)地質(zhì)統(tǒng)計(jì)學(xué)反演基礎(chǔ)上發(fā)展起來的一種高精度模擬表征的方法。

地震波形指示反演的主要思想是在等時(shí)地層格架約束下，將地震波形的薄層干涉特征作為判別、優(yōu)化反射系數(shù)結(jié)構(gòu)的控制條件代替變差函數(shù)，優(yōu)選有效樣本井，模擬砂體縱向分布結(jié)構(gòu)，將井約束地震反演與地震指示的井模擬相結(jié)合，實(shí)現(xiàn)井震聯(lián)合反演。變差函數(shù)能夠表征區(qū)域化變量的空間結(jié)構(gòu)性和隨機(jī)性，反映了區(qū)域化變量在某個(gè)方向上某一距離范圍內(nèi)的變化程度，但由于地質(zhì)結(jié)構(gòu)的非均質(zhì)性，變差函數(shù)所描述的空間變異性不能精確表征空間沉積環(huán)境的變化。而地震波形的橫向變化與沉積作用有關(guān)，并且地震具有空間密集分布的優(yōu)勢(shì)，能夠較好地體現(xiàn)沉積環(huán)境的變化及其對(duì)儲(chǔ)層組合結(jié)構(gòu)的控制作用。因此，利用地震波形橫向變化特征來表征儲(chǔ)層空間變化規(guī)律，能更好地體現(xiàn)沉積要素的影響，實(shí)現(xiàn)相控隨機(jī)反演。在進(jìn)行隨機(jī)模擬時(shí)，傳統(tǒng)變差函數(shù)是根據(jù)所有井統(tǒng)計(jì)出的變程優(yōu)選統(tǒng)計(jì)樣本進(jìn)行模擬，考慮的是空間變化程度，和距離相關(guān)（圖4a）。波形指示反演主要是根據(jù)波形相似性優(yōu)選統(tǒng)計(jì)樣本，通過將預(yù)測(cè)道地震波形與所有已知井旁道地震波形進(jìn)行對(duì)比，優(yōu)選出最相似的若干井樣本（圖4b），再對(duì)這些井進(jìn)行不同頻段下的曲線濾波比較，尋找共性結(jié)構(gòu)特征并建立初始模型，最后在貝葉斯框架下根據(jù)樣本井的分布特征進(jìn)行克里金概率模擬，得到每一個(gè)樣點(diǎn)值的概率分布。模擬過程充分利用了空間密集分布的地震數(shù)據(jù)并體現(xiàn)了相控模擬的思想，其頻率成分是一個(gè)由低到高逐步確定的過程，高頻成分的整體確定性相比傳統(tǒng)隨機(jī)模擬得到大幅提高。與傳統(tǒng)隨機(jī)模擬不同的是，地震波形指示模擬建立初始模型的過程不是采用序貫的方式進(jìn)行，即待模擬點(diǎn)的模擬數(shù)據(jù)并不參與下一個(gè)未知點(diǎn)的模擬過程。這樣做使得模擬過程更加符合地震波形相似的樣本優(yōu)選原則，使模擬結(jié)果的中頻（地震頻帶）符合地震反演結(jié)果，超出地震頻帶的高頻成分與樣本井結(jié)構(gòu)特征一致，得到寬頻帶阻抗輸出。

圖4 統(tǒng)計(jì)樣本優(yōu)選

地震波形指示反演是在空間結(jié)構(gòu)化數(shù)據(jù)指導(dǎo)下不斷尋優(yōu)的過程，即參照空間分布距離和地震波形相似性兩個(gè)因素對(duì)所有井按關(guān)聯(lián)度排序，優(yōu)選與預(yù)測(cè)點(diǎn)關(guān)聯(lián)度高的井作為初始模型，對(duì)高頻成分進(jìn)行無偏最優(yōu)估計(jì)，并保證最終反演的地震波形與原始地震波形一致。

地震波形指示反演的流程如圖5所示，具體步驟為：

1）按照地震波形特征對(duì)已知井進(jìn)行分析，優(yōu)選與待判別道波形關(guān)聯(lián)度高的井樣本建立初始模型，并統(tǒng)計(jì)其縱波阻抗作為先驗(yàn)概率分布。為了避免大距離范圍的樣本優(yōu)選誤差，優(yōu)選過程增加最大空間距離控制，在已知井中利用波形相似性和空間距離雙變量?jī)?yōu)選中低頻結(jié)構(gòu)相似的井作為空間估值樣本。傳統(tǒng)變差函數(shù)受井位分布的影響，難以精確表征儲(chǔ)層的非均質(zhì)性，而分布密集的地震波形則可以精確表征空間結(jié)構(gòu)的低頻變化。

2）將初始模型與地震頻帶阻抗進(jìn)行匹配濾波，計(jì)算得到似然函數(shù)。如果兩口井的地震波形相似，表明這兩口井大的沉積環(huán)境相似，雖然其高頻成分可能來自不同的沉積微相差異較大，但其低頻成分具有共性，且經(jīng)過測(cè)井曲線統(tǒng)計(jì)證明，其共性頻帶范圍大幅度超出了地震有效頻帶。根據(jù)似然函數(shù)的定義，此處通過匹配濾波計(jì)算得到的似然函數(shù)描述了某一個(gè)空間位置取值為某一特定值的概率。匹配濾波將確定性信息（源自有色反演）和逐步確定性信息（源自波形指示模擬）進(jìn)行融合之后產(chǎn)生了一個(gè)新的概率分布空間。利用這一特性可以增強(qiáng)反演結(jié)果低頻段的確定性，同時(shí)約束了高頻的取值范圍，使反演結(jié)果確定性更強(qiáng)。

圖5 地震波形指示反演流程

3）在貝葉斯框架下聯(lián)合似然函數(shù)分布和先驗(yàn)分布得到后驗(yàn)概率分布并將其作為目標(biāo)函數(shù)，不斷擾動(dòng)模型參數(shù)，將后驗(yàn)概率分布函數(shù)最大時(shí)的解作為有效的隨機(jī)實(shí)現(xiàn)，取多次有效實(shí)現(xiàn)的均值作為期望值輸出。

地震波形指示反演結(jié)果在空間上體現(xiàn)了沉積相帶的約束，平面上更符合沉積規(guī)律和特點(diǎn)：①在貝葉斯框架下將地震、地質(zhì)和測(cè)井信息有效結(jié)合，利用地震信息指導(dǎo)井參數(shù)高頻模擬，是一種全新的井震結(jié)合方式，較好地減少了地震噪聲對(duì)反演結(jié)果的影響；②利用地震波形特征代替變差函數(shù)分析儲(chǔ)層空間結(jié)構(gòu)變化，提高了橫向分辨率，且更符合平面地質(zhì)規(guī)律，具有相控意義；③采用全局優(yōu)化算法使反演確定性大大增強(qiáng)（從完全隨機(jī)到逐步確定）；④對(duì)井位分布沒有嚴(yán)格要求，適用性更廣。

3 應(yīng)用效果分析

3．1 研究區(qū)地質(zhì)背景

肇35井區(qū)位于松遼盆地三肇凹陷西部，宋芳屯油田以南，永樂油田以北，面積約60km2，共有探井12口。地震剖面分析可知，泉頭組頂面T2反射同相軸為強(qiáng)波峰，在泉四段內(nèi)部發(fā)育多套砂體，由于單層厚度薄，與圍巖之間阻抗差異小，因此，在地震記錄上沒有明顯的反射界面，受調(diào)諧作用影響主要對(duì)應(yīng)強(qiáng)波峰下的波谷反射，砂巖發(fā)育的厚度及期次對(duì)波形形狀有一定影響。圖6是Zhao43井和Zhao48井聲波曲線和井旁地震道，其中Zhao43井F11油層組共發(fā)育兩套砂巖，平均厚度為3．2m，Zhao48井發(fā)育3套砂巖，平均厚度為2．3m。兩口井泉四段內(nèi)部砂巖雖然均表現(xiàn)為波谷特征，但地震波形仍存在一定差異，表明波形的橫向變化指示了其高頻反射結(jié)構(gòu)的變化。分析T2目的層的地質(zhì)和地球物理響應(yīng)特征可知，扶余油層為砂泥巖互層，儲(chǔ)層特征為高伽馬、低聲波時(shí)差?；谏皫r的這種測(cè)井響應(yīng)特征，可以通過阻抗反演區(qū)分砂泥巖。

圖6 Zhao43井（a）和Zhao48井（b）波形差異分析

3．2 測(cè)井?dāng)?shù)據(jù)分析處理

分析、編輯與校正測(cè)井曲線形成高質(zhì)量的測(cè)井資料，是定量化油藏描述研究的基礎(chǔ)。測(cè)井?dāng)?shù)據(jù)用于地震數(shù)據(jù)體標(biāo)定，然而收集到的實(shí)測(cè)測(cè)井曲線通常會(huì)有部分?jǐn)?shù)據(jù)缺失、質(zhì)量較差、多井之間缺乏一致性等問題。因此，在與地震資料進(jìn)行聯(lián)合應(yīng)用之前，必須對(duì)測(cè)井資料進(jìn)行質(zhì)量控制，以提供一致性較好、相對(duì)完整的測(cè)井?dāng)?shù)據(jù)集。針對(duì)儲(chǔ)層敏感曲線開展多井一致性檢查和標(biāo)準(zhǔn)化處理，以獲得縱向質(zhì)量可靠、橫向具有一致性的測(cè)井資料。

本文將研究區(qū)Zhao43井選為測(cè)井評(píng)價(jià)的標(biāo)準(zhǔn)井，將T2作為標(biāo)志層（測(cè)井曲線形態(tài)特征穩(wěn)定并且全部鉆遇的目的層），采用直方圖法對(duì)全區(qū)12口井進(jìn)行標(biāo)準(zhǔn)化處理，剔除異常值，消除了測(cè)井資料對(duì)地震反演結(jié)果的影響。由于同一個(gè)區(qū)塊相同的地層沉積環(huán)境、沉積年代相近且地層性質(zhì)相似，在直方圖上的分布概型相似，因此可以對(duì)不同的測(cè)井曲線進(jìn)行相應(yīng)調(diào)整，進(jìn)而實(shí)現(xiàn)多井曲線標(biāo)準(zhǔn)化的目的。通過分析目的層段測(cè)井曲線直方圖分布，統(tǒng)計(jì)出各井的偏離值并加以校正（圖7）。

圖7 聲波曲線標(biāo)準(zhǔn)化前（a）后（b）直方圖對(duì)比

3．3 反演參數(shù)分析

地震波形指示反演是一種改進(jìn)的統(tǒng)計(jì)學(xué)反演方法，反演結(jié)果具有“低頻確定、高頻隨機(jī)”的特點(diǎn)。低頻主要是受地震頻帶及地震相的影響，高頻則主要受同沉積結(jié)構(gòu)樣本的控制，頻率越高，隨機(jī)性就越強(qiáng)。因此，地震波形指示反演的關(guān)鍵參數(shù)為有效樣本數(shù)和最佳截止頻率。

3．3．1 有效樣本數(shù)

有效樣本數(shù)是地震波形指示反演中的一個(gè)重要參數(shù)，主要表征地震波形空間變化對(duì)儲(chǔ)層的影響程度，一般參照已知井的統(tǒng)計(jì)結(jié)果，在質(zhì)控中利用樣本數(shù)和地震相關(guān)性統(tǒng)計(jì)分析得到。由相關(guān)性曲線可見，相關(guān)性隨樣本數(shù)的增加而逐漸增大，達(dá)到一定程度后不再隨樣本數(shù)增加而增加，表明更多的樣本無助于預(yù)測(cè)精度的提高，其相關(guān)性最大時(shí)的樣本數(shù)就是最佳有效樣本數(shù)。該參數(shù)也和總樣本數(shù)有關(guān)，通常較大的樣本數(shù)表明地層彈性參數(shù)變化小，非均質(zhì)性弱，在橫向變化快、非均質(zhì)性強(qiáng)的地區(qū)可適當(dāng)減小樣本數(shù)。本文根據(jù)研究區(qū)的樣本分析擬合結(jié)果，優(yōu)選有效樣本數(shù)為2（圖8）。

3．3．2 最佳截止頻率

圖8 地震波形指示反演有效樣本數(shù)優(yōu)選與質(zhì)控

圖9 地震波形指示反演最佳截止頻率優(yōu)選與質(zhì)控

在確定了有效樣本數(shù)后，需進(jìn)行最佳截止頻率的確定。與有效樣本數(shù)分析方法類似，計(jì)算各井不同頻率情況下與其它井的波形相關(guān)性，擬合所有相關(guān)性曲線。如果處理的目標(biāo)偏向于反演的確定性，該參數(shù)不宜設(shè)置太高，反之，如果處理的目標(biāo)偏向于反演分辨率，能夠接受隨機(jī)的結(jié)果，可以設(shè)置較高的截止頻率。本次為了提高反演高頻成分的確定性，通過參數(shù)測(cè)試確定最大截止頻率為200～250Hz（圖9）。

3．4 反演結(jié)果對(duì)比分析

對(duì)井約束稀疏脈沖、地質(zhì)統(tǒng)計(jì)學(xué)隨機(jī)反演和地震波形指示反演三種方法的反演結(jié)果進(jìn)行了對(duì)比分析（圖10～圖13）。圖10為過Zhao43，Zhao48井的井約束稀疏脈沖反演剖面，圖中紅、黃色代表高阻抗，綠、藍(lán)色等代表低阻抗，可以看出剖面分辨率低，薄層砂體無法有效識(shí)別，且強(qiáng)阻抗界面影響嚴(yán)重。圖11為疊后地質(zhì)統(tǒng)計(jì)學(xué)反演剖面，垂向分辨率較高，但橫向連續(xù)性較差，這是由于該地區(qū)井?dāng)?shù)較少，無明顯統(tǒng)計(jì)學(xué)規(guī)律，且塊金值較大，無法擬合準(zhǔn)確的變差函數(shù)模型（圖12）。因此，該地區(qū)不建議采用地質(zhì)統(tǒng)計(jì)學(xué)方法。圖13為利用地震波形指示方法進(jìn)行的波阻抗反演剖面，垂向分辨率較高，橫向變化自然，與沉積規(guī)律相符，反演剖面縱向上大套地層速度結(jié)構(gòu)合理，砂巖高速異常較為突出。該方法反演結(jié)果垂向上由于充分利用了測(cè)井信息，分辨率遠(yuǎn)高于地震剖面，在目的層段能夠有效反映2m以上的薄層；橫向上，高阻抗層的展布受地震波形橫向變化的影響，表明該方法能夠充分利用地震、測(cè)井信息，有效提高井間薄儲(chǔ)層預(yù)測(cè)能力，同時(shí)可有效避免上覆T2強(qiáng)阻抗層的干擾。3．5 反演結(jié)果驗(yàn)證

圖10 過Zhao43，Zhao48井的井約稀疏脈沖反演剖面

圖11 過Zhao43，Zhao48井地質(zhì)統(tǒng)計(jì)學(xué)反演剖面

圖12 地質(zhì)統(tǒng)計(jì)學(xué)變差函數(shù)擬合

圖13 過Zhao43，Zhao48井地震波形指示方法反演剖面

地質(zhì)統(tǒng)計(jì)學(xué)反演常用后驗(yàn)井的方法來驗(yàn)證反演結(jié)果的可靠性。工區(qū)內(nèi)共有12口井，選出2口井（Zhao43，F(xiàn)176－19）作為后驗(yàn)井不參與反演。圖14為所有井參與地震波形指示反演的剖面，圖15為部分井參與地震波形指示反演的剖面，可以看出，后驗(yàn)井Zhao43，F(xiàn)176－19與反演剖面吻合度較高，且兩次反演剖面面貌一致，說明地震波形指示方法反演結(jié)果可靠，兼具描述性與預(yù)測(cè)性。

圖14 所有井參與的地震波形指示反演剖面

圖15 部分井參與的地震波形指示反演剖面

4 結(jié)論與認(rèn)識(shí)

1）松遼盆地三肇凹陷普遍發(fā)育T2強(qiáng)反射界面，給扶余油層砂體預(yù)測(cè)帶來較大困難。分析表明，T2強(qiáng)反射界面主要是由于上下地層阻抗差異大，形成較大的反射系數(shù)，在地震波調(diào)諧過程中占據(jù)主導(dǎo)作用。

2）地震波形指示反演利用地震波形的橫向變化優(yōu)選反射結(jié)構(gòu)相似的樣本井建立初始模型，較好地克服了變差函數(shù)對(duì)井分布要求高、相控性差的缺點(diǎn)，有效提高了薄儲(chǔ)層垂向分辨率和橫向識(shí)別能力。

［1］ RICKER N．The form and nature of seismic waves and the structure of seismograms［J］．Geophysics，1940，5（4）：348－366

［2］ WIDESS M A．How thin is a thin bed？［J］．Geophysics，1973，38（6）：1176－1180

［3］ SHERIFF R E．Limitations on resolution of seismic reflections and geologic detail derivable from them：section1，fundamentals of stratigraphic interpretation of seismic data［J］．AAPG Memoir，1977，26：3－14

［4］ LEVY S，F(xiàn)ULLAGAR P K．Reconstruction of a sparse spike train from a portion of its spectrum and application to high－resolution deconvolution［J］．Geophysics，1981，46（9）：1235－1243

［5］ 郭朝斌，楊小波，陳紅岳．等．約束稀疏脈沖反演在海通集洼陷中的應(yīng)用［J］．石油物探，2006，45（4）：397－400

GUO C B，YANG X B，CHEN H Y，et al．Constrained sparse pulse inversion research in north of Haitongji depression［J］．Geophysical Prospecting for Petroleum，2006，45（4）：397－400

［6］ 王權(quán)鋒，郭科．約束稀疏脈沖反演在儲(chǔ)層預(yù)測(cè)中的應(yīng)用［J］．測(cè)井技術(shù)，2008，32（1）：33－36

WANG Q F，GUO K．Constrained sparse spike inversion applied in the reservoir prediction［J］．Well Logging Technology，2008，32（1）：33－36

［7］ 王權(quán)鋒，王元君．寬帶約束地震反演在蘭尕地區(qū)儲(chǔ)層預(yù)測(cè)中的應(yīng)用［J］．成都理工大學(xué)學(xué)報(bào)（自然科學(xué)版），2008，35（1）：48－51

WANG Q F，WANG Y J．Application of the broadband constrained seismic inversion to the prediction of reservoirs in the Langa region of the Tahe oilfield，Xinjiang，China［J］．Journal of Chengdu University of Technology（Science＆Technology Edition），2008，35（1）：48－51

［8］ 呂鐵良，王永剛，謝萬(wàn)學(xué)．等．稀疏脈沖反演技術(shù)在井間地震反演中的應(yīng)用［J］．石油物探，2007，46（1）：58－63

LV T L，WANG Y G，XIE W X，et al．Application of sparse pulse inversion in cross well seismic data［J］．Geophysical Prospecting for Petroleum，2007，46（1）：58－63

［9］ DUBRULE O，THIBAUT M，LAMY P，et al．Geostatistical reservoir characterization constrained by 3D seismic data［J］．Petroleum Science，1998，4（2）：121－128

［10］ 胡瑋，韓立國(guó)，尚帥．頻散屬性在致密薄儲(chǔ)層識(shí)別中的應(yīng)用［J］．石油物探，2014，53（3）：367－372

HU W，HAN L G，SHANG S．Application of frequency－dispersion attributes in identification of tight thin interbedded reservoirs［J］．Geophysical Prospecting for Petroleum，2014，53（3）：367－372

［11］ HAAS A，DUBRULE O．Geostatistical inversion：a sequential method for stochastic reservoir modeling constrained by seismic data［J］．First Break，1994，12（11）：561－569

［12］ 郭見樂，唐文榜，李宗杰，等．薄儲(chǔ)層（體）含油氣性識(shí)別的HFC技術(shù)及其應(yīng)用研究［J］．石油物探，2013，52（1）：79－88

GUO J L，TANG W B，LI Z J，et al．The harmonic frequencies characteristics（HFC）analysis for identifying and predicting hydrocarbon in thin reservoir strata（bodies）［J］．Geophysical Prospecting for Petroleum，2013，52（1）：79－88

［13］ TORRES－VERDíN C，VICTORIA M，MERLETTI G，et al．Trace－based and geostatistical inversion of 3－D seismic data for thin－sand delineation：an application in San Jorge Basin，Argentina［J］．The Leading Edge，18（9）：1070－1077

［14］ 孫思敏，彭仕宓．地質(zhì)統(tǒng)計(jì)學(xué)反演方法及其在薄層砂體預(yù)測(cè)中的應(yīng)用［J］．西安石油大學(xué)學(xué)報(bào)（自然科學(xué)版），2007，22（1）：41－48

SUN S M，PENG S M．Geostatistical inversion method and its application in the prediction of thin reservoirs［J］．Journal of Xi’an Shiyou University（Natural Science Edition），2007，22（1）：41－48

［15］ 孫思敏，彭仕宓．地質(zhì)統(tǒng)計(jì)學(xué)反演及其在吉林扶余油田儲(chǔ)層預(yù)測(cè)中的應(yīng)用［J］．物探與化探，2007，31（1）：51－56

SUN S M，PENG S M．Geostatistical inversion method and its application to reservoir prediction of Fuyu oilfield in Jilin province［J］．Geophysical＆Geochemical Exploration，2007，31（1）：51－56

［16］ 王西文，劉全新，周嘉璽，等．精細(xì)儲(chǔ)層預(yù)測(cè)技術(shù)在板南5—3井區(qū)的應(yīng)用［J］．石油物探，2003，42（3）：389－394

WANG X W，LIU Q X，ZHOU J X，et al．Application of fine reservoir prediction techniques in Bannan 5—3 oilfied［J］．Geophysical Prospecting for Petroleum，2003，42（3）：389－394

［17］ 盛述超，畢建軍，李維振，等．關(guān)于地震波形指示模擬反演（SMI）方法的研究［J］．內(nèi)蒙古石油化工，2015（17）：147－151

SHENG S C，BI J J，LI W Z，et al．The study of seismic motion inversion（SMI）method［J］．Inner Mongolia Petrochemical Engineering，2015（17）：147－151

［18］ 劉宗堡，呂延防，付曉飛，等．三肇凹陷扶余油層沉積特征及油氣成藏模式［J］．吉林大學(xué)學(xué)報(bào)（地球科學(xué)版），2009，39（6）：998－1005

LIU Z B，LV Y Y，F(xiàn)U X F，et al．Sedimentary chracteristics and hydrocarbon accumulation model of Fuyu reservoir in Sanzhao depression［J］．Journal of Jilin University（Earth Science Edition），2009，39（6）：998－1005

（編輯：戴春秋）

High precision prediction of thin reservoir under strong shielding effect and its application：a case study from Sanzhao Depression，Songliao Basin

GU Wen，ZHANG Xiong，XU Min，LIANG Hong，ZHANG Dongjun，LUO Jing，ZHENG Hong
（Geophysical Prospecting Company of Chuanqing Drilling Engineering Company Limited，CNPC，Chengdu610213，China）

Fuyu oil layer in Sanzhao Depression，Songliao Basin developed typical sedimental system of shallow water delta that had thin layer of sandbody and fast transverse change．There is no obvious seismic reflection characteristics which affected by the shield effect from strong T2reflection．How to reduce the strong shield tuning effect，highlight the responses from the fuyu oil layer thin sandbodies and has a fine description，it’s always been a problem in exploration and development in the target．Strong impedance shielding phenomenon is mainly associated with the tuning effect of strong reflection coefficient in the study．When the reflection coefficient is high in the formation，it causes the low seismic frequency．Thus the low reflection coefficient will be"swamped"，and the influence of the strong impedance interface will be highlighted．Therefore，the main consideration of weakening shielding effect should be given to reduce the influence of the seismic wave interference effect．Due to the limitations of increased resolution by seismic processing，through a large number of trials，we adopted a new high－resolution inversion technology that called seismic meme inversion．It makes full use of lateral seismic waveform space change information instead of conventional variogram，and realizes seismic and well joint characterization of high frequency components．It effectively reduces the influence of strong seismic impedance interface and achieves the high precision thin reservoir prediction under strong shielding effect in the target area．

Fuyu oil layer，strong reflection shielding，interference effect，variogram，seismic motion inversion

P631

A

1000－1441（2017）03－0439－10

10．3969／j．issn．1000－1441．2017．03．014

顧雯，章雄，徐敏，等．強(qiáng)屏蔽下薄儲(chǔ)層高精度預(yù)測(cè)研究［J］．石油物探，2017，56（3）：439－448

GU Wen，ZHANG Xiong，XU Min，et al．High precision prediction of thin reservoir under strong shielding effect and its application：a case study from Sanzhao Depression，Songliao Basin［J］．Geophysical Prospecting for Petroleum，2017，56（3）：439－448

2016－08－31；改回日期：2017－01－23。

顧雯（1986—），女，碩士，工程師，主要從事儲(chǔ)層預(yù)測(cè)方面的研究工作。

國(guó)家科技重大專項(xiàng)（2011ZX05013－001）與CNPC重大科技專項(xiàng)“陸上地震勘探關(guān)鍵裝備和軟件研發(fā)”（2013E－3808）聯(lián)合資助。

This research is financially supported by the National Science and Technology Major Project of China（Grant No．2011ZX05013－001）and the Major Science and Technology Project of CNPC"key equipment and software development in land seismic exploration”（Grant No．2013E－3808）．