華曉琴 史樂 梁濤

摘 ?要：塔里木盆地塔中Ⅲ區(qū)奧陶系碳酸鹽巖儲(chǔ)層是近年來勘探重點(diǎn)區(qū)域，良里塔格組儲(chǔ)層內(nèi)沉積泥質(zhì)灰?guī)r條帶，以及一間房組洞穴頂部存在的沉積泥質(zhì)條帶，給儲(chǔ)層預(yù)測(cè)帶來了較大難度。為提高儲(chǔ)層預(yù)測(cè)精度，分析奧陶系良里塔格組和一間房組儲(chǔ)層內(nèi)泥質(zhì)展布規(guī)律，建立沉積泥質(zhì)灰?guī)r條帶背景縱波阻抗體，完成了無井約束疊后地質(zhì)統(tǒng)計(jì)學(xué)反演。研究結(jié)果表明，無井反演儲(chǔ)層預(yù)測(cè)縱向發(fā)育位置、儲(chǔ)層級(jí)別和已鉆井吻合較好。利用疊前同時(shí)反演得到縱橫波速比，通過縱波阻抗與縱橫波速比交會(huì)圖分析法，劃分了一間房組優(yōu)質(zhì)儲(chǔ)層、泥質(zhì)充填儲(chǔ)層、基質(zhì)灰?guī)r分布區(qū)間，對(duì)洞穴型儲(chǔ)層泥質(zhì)充填、疊前流體預(yù)測(cè)反演進(jìn)行探索，確保儲(chǔ)層預(yù)測(cè)的可靠性。

關(guān)鍵詞：碳酸鹽巖;測(cè)井;疊前反演;地質(zhì)統(tǒng)計(jì)學(xué)反演;儲(chǔ)層預(yù)測(cè)

塔中Ⅲ區(qū)位于塔里木盆地塔中I號(hào)坡折帶西端內(nèi)帶，構(gòu)造較平緩，有利于油氣大面積成藏。奧陶系良里塔格組、一間房組和鷹山組為油氣富集有利層位，橫向上多類型、多套儲(chǔ)集體連片發(fā)育，縱向上不同層段儲(chǔ)層疊置發(fā)育，具多層系含油特征。作為塔里木盆地油氣勘探的主戰(zhàn)場(chǎng)之一，塔中地區(qū)近年來加快了油氣鉆探進(jìn)程，為提高鉆井成功率，對(duì)儲(chǔ)層的預(yù)測(cè)和識(shí)別提出了更高要求。地震反演預(yù)測(cè)是儲(chǔ)層預(yù)測(cè)的重要技術(shù)手段，地震反演技術(shù)方法較多，針對(duì)碳酸鹽巖儲(chǔ)層目前應(yīng)用較廣的為疊后波阻抗反演[1]。

地震反演技術(shù)在國(guó)內(nèi)外類似地質(zhì)條件區(qū)塊已開展大量工作，取得了良好效果。朱超等針對(duì)碳酸鹽巖儲(chǔ)層預(yù)測(cè)的復(fù)雜性，基于古地貌控相、相控儲(chǔ)層、儲(chǔ)層控藏的地質(zhì)認(rèn)識(shí)，采用分頻非線性反演方法，實(shí)現(xiàn)對(duì)復(fù)雜地質(zhì)體的準(zhǔn)確預(yù)測(cè)和描述[2]。任麗丹等針對(duì)塔里木盆地順南地區(qū)儲(chǔ)層埋藏深、地震資料信噪比低的問題，利用AVO同步反演技術(shù)，提高碳酸鹽巖儲(chǔ)層預(yù)測(cè)精度[3]。王保才等采用疊前彈性參數(shù)反演技術(shù)在順南地區(qū)開展碳酸鹽巖儲(chǔ)層預(yù)測(cè)，提高縱向上儲(chǔ)層預(yù)測(cè)精度[4]。張磊等利用地質(zhì)統(tǒng)計(jì)學(xué)反演對(duì)輪古油田碳酸鹽巖儲(chǔ)層進(jìn)行預(yù)測(cè)，較好地解決了該區(qū)儲(chǔ)層預(yù)測(cè)難的問題[5]。

塔中Ⅲ區(qū)奧陶系良里塔格組沉積泥質(zhì)灰?guī)r條帶廣泛分布，一間房組洞穴頂部常見泥質(zhì)充填，在地震上儲(chǔ)層與泥質(zhì)響應(yīng)難以區(qū)分。鑒于奧陶系地質(zhì)條件復(fù)雜，受沉積泥質(zhì)灰?guī)r影響，在常規(guī)反演剖面上，良里塔格組儲(chǔ)層難以識(shí)別，一間房組儲(chǔ)層頂部難以刻畫，而水平井鉆探又對(duì)一類儲(chǔ)層縱向發(fā)育位置及二類儲(chǔ)層識(shí)別提出了迫切要求。因此，在借鑒相關(guān)研究經(jīng)驗(yàn)的基礎(chǔ)上，針對(duì)奧陶系儲(chǔ)層實(shí)際地質(zhì)特征，以巖石物理研究為基礎(chǔ)、充分利用地質(zhì)、地震、測(cè)井等資料，開展儲(chǔ)層預(yù)測(cè)與流體分布預(yù)測(cè)，解決了以往儲(chǔ)層泥質(zhì)條帶對(duì)儲(chǔ)層預(yù)測(cè)的干擾，利用地質(zhì)統(tǒng)計(jì)學(xué)反演和疊前同時(shí)反演，識(shí)別了洞穴型儲(chǔ)層充填特征，有效提高儲(chǔ)層預(yù)測(cè)精度。

1 ?泥質(zhì)特征分析

奧陶系良里塔格組底部發(fā)育一套沉積泥質(zhì)灰?guī)r，一間房-鷹山組泥質(zhì)洞穴儲(chǔ)層較發(fā)育。據(jù)泥質(zhì)充填程度的不同，分為未完全充填泥洞和完全充填泥洞。目前高自然伽馬的泥質(zhì)發(fā)育層段多在生產(chǎn)層段內(nèi)，有些高自然伽馬的充填程度低的泥質(zhì)洞穴儲(chǔ)層試油為見油氣。

良里塔格組底部泥質(zhì)灰?guī)r條帶、一間房-鷹山組充填泥洞以及地層水高含鈾等多種因素，使工區(qū)內(nèi)自然伽馬測(cè)井在良里塔格及一間房組呈高伽馬異常。因此，通過高伽馬成因分析，分析對(duì)應(yīng)的測(cè)井響應(yīng)特征，消除含鈾地層水的影響，識(shí)別沉積泥質(zhì)與充填泥質(zhì)的發(fā)育特征，探索泥質(zhì)完全充填洞穴與泥質(zhì)半充填洞穴在測(cè)井及地球物理上的差異，指導(dǎo)含泥質(zhì)灰?guī)r背景阻抗的建立，提高儲(chǔ)層反演精度。

1.1 ?泥質(zhì)層段測(cè)井響應(yīng)特征

高質(zhì)量的測(cè)井資料是定量化油藏描述的基本條件，在和地震資料進(jìn)行聯(lián)合應(yīng)用之前，必須對(duì)測(cè)井資料進(jìn)行質(zhì)量控制，以提供一致性高、相對(duì)完整的測(cè)井?dāng)?shù)據(jù)集[6]。在此基礎(chǔ)之上，以測(cè)井曲線建立的體積模型作為油藏描述的特征參數(shù)，結(jié)合巖石彈性參數(shù)，對(duì)油藏特征進(jìn)行解釋。

1.1.1 ?沉積成因泥質(zhì)灰?guī)r測(cè)井響應(yīng)

高GR測(cè)井響應(yīng)對(duì)應(yīng)的沉積泥質(zhì)灰?guī)r條帶主要分布于良里塔格組，尤其良三段底部沉積泥質(zhì)灰?guī)r，全區(qū)發(fā)育穩(wěn)定，最大厚度15 m。伽馬曲線呈鐘形，表現(xiàn)為反旋回沉積特點(diǎn)，聲波時(shí)差47～70 μs/ft，密度多大于2.55 g/cm3，波阻抗約為14 500～ ? ? ? ? ? 17 200 g/cm3?m/s，成像測(cè)井呈水平層狀特征。

1.1.2 ?充填泥質(zhì)測(cè)井響應(yīng)

一間房組頂高GR測(cè)井響應(yīng)對(duì)應(yīng)于洞穴型、孔洞型儲(chǔ)層的充填泥質(zhì)，厚度約為5～10 m。據(jù)成像、測(cè)井曲線特征等資料，分析泥質(zhì)的充填程度，將該套泥質(zhì)充填孔洞劃分為全充填和半充填。

一間房組泥質(zhì)充填程度高層段 ?伽馬值為15～100 API，曲線多呈指狀或刺刀狀，無旋回特征，聲波時(shí)差為50～70 μs/ft，密度多大于2.55 g/cm3，波阻抗為12 500～16 600 g/cm3?m/s，成像測(cè)井呈堆積特征、頂?shù)壮拭黠@不整合接觸。

一間房組泥質(zhì)充填程度低層段 ?伽馬值為40～130 API，聲波時(shí)差為60～100 μs/ft，密度一般小于2.55 g/cm3，波阻抗小于14 000 g/cm3?m/s，與成像測(cè)井相比呈顏色較深的暗斑狀。與充填程度高層段相比，低充填層段伽馬值略小，但其速度與密度均偏低，波阻抗具有明顯的差異，說明伽馬不能直接反映洞穴的泥質(zhì)充填程度。

對(duì)研究區(qū)測(cè)井資料統(tǒng)計(jì)分析，得到沉積泥、全充填泥及部分充填泥層段測(cè)井曲線特征（表1）。總之，沉積泥質(zhì)和充填程度不同泥洞，其電性、物性和彈性特征等均表現(xiàn)出較大差異。沉積泥質(zhì)和全充填泥洞，兩者常規(guī)測(cè)井響應(yīng)特征差異較小，但FMI全井眼電成像上，沉積泥質(zhì)呈水平或低角度層狀，具明顯的旋回沉積，而全充填泥洞表現(xiàn)為快速堆積特征，頂?shù)壮拭黠@不整合接觸，無明顯旋回。未完全充填泥洞在常規(guī)測(cè)井響應(yīng)特征上表現(xiàn)為異于前兩者的特征，密度明顯降低、聲波時(shí)差變大、中子變大，深淺電阻率具明顯分異，井徑普遍擴(kuò)徑，F(xiàn)MI電成像上表現(xiàn)為顏色較深的暗色斑塊狀。

1.2 ?泥質(zhì)層段地震響應(yīng)特征

良里塔格沉積泥質(zhì)條帶地球物理響應(yīng)特征，泥質(zhì)含量小于20%，厚度小于15 m，對(duì)應(yīng)的地震反射較弱，常湮沒于背景反射，總體呈較弱腹波反射特征。一間房組泥質(zhì)充填程度高層段，對(duì)應(yīng)的地震反射特征接近沉積泥質(zhì)地震反射特征，對(duì)應(yīng)的地震反射弱，常湮沒于背景反射。一間房組泥質(zhì)充填程度低層段，地震剖面通常表現(xiàn)為強(qiáng)反射特征[2]。

通過去除泥質(zhì)條帶和泥質(zhì)半充填阻抗，去除沉積泥質(zhì)條帶后，同相軸變?yōu)槿醴瓷?，表明半充填泥質(zhì)層段在地震上有明顯儲(chǔ)層響應(yīng)特征。當(dāng)一間房沉積泥質(zhì)條帶與良里塔格組部分充填層段都存在時(shí)，地震剖面上為強(qiáng)振幅反射，洞頂位于波谷到波峰零相位，當(dāng)去除沉積泥質(zhì)條帶后，洞頂對(duì)應(yīng)位置及地震反射特征無變化;當(dāng)去除部分充填洞穴后，呈弱反射特征。

2 ?“相控”疊后地質(zhì)統(tǒng)計(jì)學(xué)反演

為建立適合本區(qū)地質(zhì)特點(diǎn)的儲(chǔ)層預(yù)測(cè)方法，先采取疊后反演實(shí)現(xiàn)儲(chǔ)層的準(zhǔn)確預(yù)測(cè)，消除儲(chǔ)層內(nèi)泥質(zhì)灰?guī)r條帶的影響，再利用疊前反演判別儲(chǔ)層內(nèi)泥質(zhì)充填特征，用于預(yù)測(cè)儲(chǔ)層內(nèi)部泥質(zhì)充填特征，提高儲(chǔ)層預(yù)測(cè)可靠性。

2.1 ?沉積泥質(zhì)巖性背景阻抗體構(gòu)建

良里塔格組碳酸鹽巖臺(tái)地泥質(zhì)灰?guī)r廣泛發(fā)育，而一間房組儲(chǔ)層頂部往往被泥質(zhì)半充填或全充填，沉積泥質(zhì)灰?guī)r條帶與泥質(zhì)充填的儲(chǔ)層頂部的縱波阻抗相互疊置，不具明顯差異，因此從常規(guī)地質(zhì)統(tǒng)計(jì)學(xué)反演上無法區(qū)分沉積泥質(zhì)灰?guī)r條帶與儲(chǔ)層頂部。在水平井鉆進(jìn)過程中，當(dāng)鉆遇泥質(zhì)灰?guī)r層段，很難從隨鉆曲線或常規(guī)反演剖面上區(qū)分是否鉆遇儲(chǔ)層頂部;其次消除沉積泥質(zhì)地震響應(yīng)，預(yù)測(cè)真正儲(chǔ)層頂部，規(guī)避泥質(zhì)灰?guī)r條帶，也是水平井地質(zhì)及工程設(shè)計(jì)所急需解決的難題[3]。

由于地震響應(yīng)是泥質(zhì)灰?guī)r（含泥質(zhì)）與儲(chǔ)層的疊合響應(yīng)，可通過沉積相約束建立沉積泥質(zhì)灰?guī)r背景阻抗，然后利用地質(zhì)統(tǒng)計(jì)學(xué)反演將與背景阻抗相耦合的儲(chǔ)層阻抗求出，得到去除沉積泥質(zhì)灰?guī)r橫向變化影響的縱波阻抗體。

2.2 ?巖性劃分、巖性比例和信噪比

巖性劃分是地質(zhì)統(tǒng)計(jì)學(xué)的一個(gè)重要部分，因?yàn)榈刭|(zhì)統(tǒng)計(jì)學(xué)反演中幾乎所有的參數(shù)都需分巖性進(jìn)行分析。一間房和鷹山組為碳酸鹽巖地層，石灰?guī)r儲(chǔ)集空間發(fā)育大量裂縫并伴生溶蝕孔洞。通過對(duì)地質(zhì)條件和測(cè)井曲線特征的綜合分析，確定地質(zhì)統(tǒng)計(jì)學(xué)反演的巖性劃分標(biāo)準(zhǔn)如下：Ⅰ類儲(chǔ)層：PORT>0.05;Ⅱ類儲(chǔ)層：0.045>PORT>0.02;非儲(chǔ)層：PORT<0.02，PORT為總孔隙度。

地質(zhì)統(tǒng)計(jì)學(xué)反演的這幾個(gè)參數(shù)要考慮到實(shí)際的地質(zhì)情況和客觀資料評(píng)價(jià)，在認(rèn)識(shí)的基礎(chǔ)上進(jìn)行多重組合測(cè)試，選擇一組最合理的組合作為最終反演的參數(shù)（表2）[7]。

2.3 ?地質(zhì)統(tǒng)計(jì)學(xué)反演

泥質(zhì)背景曲線構(gòu)建之后，需要將充填泥質(zhì)層段阻抗替代為基質(zhì)巖性阻抗，完成了沉積泥質(zhì)背景曲線構(gòu)建。對(duì)巖石物理建模正演、沉積泥質(zhì)特征曲線構(gòu)建，沉積泥質(zhì)背景阻抗體的建立、確定性反演質(zhì)控、地質(zhì)統(tǒng)計(jì)學(xué)反演參數(shù)優(yōu)選，進(jìn)行去除泥質(zhì)的地質(zhì)統(tǒng)計(jì)學(xué)反演，并得到了高分辨率的反演結(jié)果。與常規(guī)的地質(zhì)統(tǒng)計(jì)學(xué)反演結(jié)果相比，去除泥質(zhì)的地質(zhì)統(tǒng)計(jì)學(xué)反演消除了沉積泥質(zhì)灰?guī)r條帶影響，預(yù)測(cè)儲(chǔ)層與真實(shí)儲(chǔ)層相吻合，精度更高，更加接近地下實(shí)際情況（圖1）。

2.4 ?反演吻合性分析

利用工程異常、試油試采資料、測(cè)井解釋成果分別統(tǒng)計(jì)41口鉆井不同儲(chǔ)層類型鉆遇長(zhǎng)度與對(duì)應(yīng)的反演預(yù)測(cè)儲(chǔ)層長(zhǎng)度，對(duì)反演儲(chǔ)層吻合率進(jìn)行分析。反演剖面上預(yù)測(cè)洞穴型儲(chǔ)層有28口井，裂縫-孔洞型儲(chǔ)層有27口，其中純Ⅱ類（不含Ⅰ類儲(chǔ)層）有9口井，發(fā)生漏失的井有13口;洞穴型儲(chǔ)層反演預(yù)測(cè)吻合率約為92%，裂縫-孔洞型儲(chǔ)層反演預(yù)測(cè)吻合率約為80%。

“相控”疊后統(tǒng)計(jì)學(xué)反演解決了良里塔格沉積泥質(zhì)對(duì)儲(chǔ)層的影響，反演的一間房?jī)?chǔ)層頂發(fā)育位置與實(shí)鉆水平井吻合（圖2）?？傮w來說，反演結(jié)果中所顯示出的儲(chǔ)層比井上的相應(yīng)厚度稍大些，主要是測(cè)井與地震的分辨率不同造成的，依據(jù)電阻率、自然伽馬等測(cè)井曲線，能將小于0.5 m的儲(chǔ)層解釋出來，而通過前面正演分析，波阻抗反演體對(duì)Ⅰ類儲(chǔ)層識(shí)別精度約為5 m，Ⅱ類儲(chǔ)層識(shí)別精度約為10 m，很難與測(cè)井儲(chǔ)層厚度對(duì)應(yīng)。

3 ?疊前同時(shí)反演

基于部分疊加數(shù)據(jù)的疊前同時(shí)反演是目前疊前反演的主要方法，通過輸入3個(gè)不同角度疊加地質(zhì)數(shù)據(jù)和對(duì)應(yīng)子波，給出不同數(shù)據(jù)（縱波阻抗、橫波阻抗和密度）的縱向變化趨勢(shì)及橫向約束范圍，最后利用反演得到的縱橫波阻抗速度比等來對(duì)儲(chǔ)層內(nèi)泥質(zhì)充填進(jìn)行判別[8]。

3.1 ?輸入數(shù)據(jù)分析

疊前時(shí)間偏移道集的目的層段最大偏移距達(dá) 5 700 m，增加部分疊加子道集的數(shù)量，疊前AVO同時(shí)反演的算法穩(wěn)定性可增強(qiáng)。但矛盾的是，增加疊加子道集的數(shù)量，降低了單個(gè)部分疊加數(shù)據(jù)體覆蓋次數(shù)，即降低了部分疊加數(shù)據(jù)體信噪比[9]。

根據(jù)實(shí)際情況，設(shè)計(jì)了5個(gè)子疊加方案。在原始CRP道集檢查中，發(fā)現(xiàn)道集中振幅水平方向不一致，存在斜行向下的噪音，這種非地質(zhì)響應(yīng)的數(shù)據(jù)信息將會(huì)對(duì)反演結(jié)果產(chǎn)生不確定影響，對(duì)原始CRP道集進(jìn)行了Ranna去噪處理，然后進(jìn)行疊加，各子疊加數(shù)據(jù)信噪比都得到一定程度的提高，信噪比會(huì)有一定程度的上升。最終選用按變偏移距方式用5個(gè)子疊加地震數(shù)據(jù)作為疊前同時(shí)反演的輸入，對(duì)本工區(qū)的疊前彈性屬性進(jìn)行研究（圖3）。

通過VSP 速度和地震最大偏移距計(jì)算出的目的層最大地震入射角接近34°，考慮到遠(yuǎn)道的信噪比較低，有效的地震入射角也在34°左右，比較適合于疊前反演工作的開展（表3）。

3.2 洞穴型儲(chǔ)層泥質(zhì)充填性分析

通過巖石物理分析，隨孔隙度變大和粘土含量增加，彈性阻抗值降低，縱橫波速度比也越來越低，泥質(zhì)充填洞和基質(zhì)阻抗與縱橫波速比有明顯界限，縱橫波速比約1.8，縱阻抗約16 000 g/cc?m/s。當(dāng)縱橫波速比小于1.7，縱波阻抗小于15 000 g/cc?m/s時(shí)，該區(qū)間主要為洞穴型儲(chǔ)層，非儲(chǔ)層在該區(qū)所占比例不超過5%;當(dāng)縱橫波速比大于1.7小于1.8，縱波阻抗大于15 000 g/cc?m/s小于15 000 g/cc?m/s區(qū)間，非儲(chǔ)層在該區(qū)所占比例不超過15%;當(dāng)縱橫波速比大于1.8，縱波阻抗小于16 000 g/cc?m/s時(shí)，該區(qū)間主要為裂縫-孔洞型儲(chǔ)層;當(dāng)縱橫波速比大于1.8，縱波阻抗大于16 000 g/cc?m/s時(shí)，基質(zhì)主要集中于該區(qū)間。

在縱波阻抗、縱橫波速比與充填程度交會(huì)圖上，當(dāng)縱橫波速比小于1.7，縱波阻抗小于15 000 g/cc?m/s時(shí)，由于該區(qū)主要為洞穴型儲(chǔ)層，完全充填層段的縱橫波速比都高于1.7，當(dāng)縱橫波速比大于1.8，縱波阻抗小于16 000 g/cc?m/s時(shí)，該區(qū)間主要為裂縫～孔洞型儲(chǔ)層，散點(diǎn)集中于縱橫波阻抗1.8～1.9，半充填相較于全充填縱波阻抗較低。當(dāng)縱波阻抗位于 ? 15 000 g/cc?m/s～16 000 g/cc?m/s時(shí)，充填儲(chǔ)層段集中發(fā)育該區(qū)間。在縱橫波速比、縱波阻抗上，充填與半充填儲(chǔ)層相互疊置。

3.3 ?洞穴型儲(chǔ)層泥質(zhì)充填性預(yù)測(cè)應(yīng)用效果

ZG17-H2是一口鉆遇泥巖的典型井，結(jié)合該井的實(shí)際生產(chǎn)情況進(jìn)行檢驗(yàn)，巖性最大概率剖面上顯示該井全井段為泥質(zhì)充填，全井段測(cè)井解釋無儲(chǔ)層，從而證實(shí)了巖性預(yù)測(cè)的可靠性（圖4）。

TZ45-H1高自然伽馬段解釋為泥質(zhì)半充填，在巖性最大概率剖面上的泥巖分布與井上實(shí)測(cè)高自然伽馬曲線比較一致，且儲(chǔ)層位置與井上出油層段較吻合，從而證實(shí)了巖性預(yù)測(cè)的可靠性（圖5）。統(tǒng)計(jì)全區(qū)實(shí)鉆井的塊狀泥巖預(yù)測(cè)精度，共有 19 口井洞穴型儲(chǔ)層解釋為泥質(zhì)充填，17口井吻合，吻合率為 89.5%（表4）。

4 ?結(jié)語

（1） 嚴(yán)格的地震資料質(zhì)控與精細(xì)的巖石物理分析是提高反演精度的保證;以“相控”地質(zhì)統(tǒng)計(jì)學(xué)反演為核心的疊前、疊后聯(lián)合反演流程適用于工區(qū)地質(zhì)情況，建議在類似地區(qū)推廣使用。

（2） 通過分層系、分類型采用不同反演流程可較好解決塔中Ⅲ區(qū)儲(chǔ)層內(nèi)泥質(zhì)識(shí)別與充填特征判別問題，實(shí)鉆表明，消除泥質(zhì)影響后的反演成果有效提高Ⅰ、Ⅱ 類儲(chǔ)層的識(shí)別精度。

參考文獻(xiàn)

[1] ? 雷芬麗，賀振華，文曉濤，等.ZH區(qū)碳酸鹽巖儲(chǔ)層孔隙度預(yù)測(cè)方法研究[J].石油天然氣學(xué)報(bào)，2010，32（3）：236-239.

[2] ? 朱超，劉占國(guó)，楊少勇，等.利用相控分頻反演預(yù)測(cè)英西湖相碳酸鹽巖儲(chǔ)層[J].石油地球物理勘探，2018，53（4）：832-841+656-657.

[3] ? 任麗丹，王鵬，劉成芳，等.疊前AVO反演技術(shù)在順南地區(qū)碳酸鹽巖儲(chǔ)層含油氣性預(yù)測(cè)中的應(yīng)用[J].工程地球物理學(xué)報(bào)，2018，15（3）：292-298.

[4] ? 王保才，劉軍，陳黎，等.疊前密度反演技術(shù)在順南地區(qū)碳酸鹽巖儲(chǔ)層預(yù)測(cè)中的應(yīng)用[J].工程地球物理學(xué)報(bào)，2016，13（4）：508-512.

[5] ? 張磊，安海亭，但光箭，等.地質(zhì)統(tǒng)計(jì)學(xué)反演在輪古油田碳酸鹽巖儲(chǔ)層預(yù)測(cè)中的應(yīng)用[J].石油地質(zhì)與工程，2016，30（2）：1-4+147.

[6] ? 趙銘海.常用疊后波阻抗反演技術(shù)評(píng)析[J].油氣地質(zhì)與采收率，2004，11（1）：36-38.

[7] ? 劉百紅，李建華.測(cè)井和地震資料寬帶約束反演的應(yīng)用[J].石油物探，2004，43（1）：76-78.

[8] ? 馬勁風(fēng)，王學(xué)軍，賈春環(huán)，等.波阻抗約束反演中的約束方法研究[J].石油物探，2000，39（2）：52-63.

[9] ? 張永剛.地震波阻抗反演技術(shù)的現(xiàn)狀和發(fā)展[J].石油物探，2002，41（4）：385-390.

Abstract：The study area of lianglitage formation reservoir sedimentary argillaceous limestone belt，and a room at the top of the cave group sedimentary shale strip，which brings great difficulty to reservoir prediction.The Lianglitag formation and the Yijianfang Formation reservoir in mud distribution analysis，establish sedimentary argillaceous limestone belt background P-wave impedance，the non well constrained poststack geostatistical inversion.The results show that the predicted vertical development location and reservoir grade are in good agreement with the drilling well.Through the analysis of drilling and rock physics research，using prestack simultaneous inversion obtained by shear wave velocity，shear wave velocity and wave impedance crossplot analysis method，a real group of high quality reservoir，mud filling reservoirs，matrix limestone Distribution Division of cave type reservoir mud filling，pre stack the fluid prediction and inversion are explored，also provide reference for such mud filling in carbonate reservoir prediction.

Key words：Carbonate rock;Logging;Pre stack inversion;Geostatistical inversion;Reservoir prediction