劉晶晶，張小莉

(西北大學(xué)地質(zhì)學(xué)系，陜西 西安 710069)

子長(zhǎng)油田長(zhǎng)6儲(chǔ)層物性測(cè)井解釋模型

劉晶晶，張小莉

(西北大學(xué)地質(zhì)學(xué)系，陜西 西安 710069)

建立物性測(cè)井解釋模型是為了提高子長(zhǎng)油田測(cè)井解釋的精度，對(duì)于儲(chǔ)層描述和油層評(píng)價(jià)具有重要意義。利用測(cè)井?dāng)?shù)據(jù)，結(jié)合物性分析化驗(yàn)資料和取芯資料，先對(duì)儲(chǔ)層的四性關(guān)系進(jìn)行簡(jiǎn)要的介紹，然后在其基礎(chǔ)上利用層點(diǎn)對(duì)應(yīng)的讀值方法建立適合子長(zhǎng)油田長(zhǎng)6儲(chǔ)層的泥質(zhì)含量和孔隙度的測(cè)井解釋模型。滲透率的解釋利用泥質(zhì)含量進(jìn)行分類(lèi)和多參數(shù)的回歸，其解釋精度得到明顯的提高。經(jīng)過(guò)驗(yàn)證，所做關(guān)于子長(zhǎng)油田長(zhǎng)6儲(chǔ)層物性參數(shù)的解釋模型準(zhǔn)確率高，效果好。

子長(zhǎng)油田;長(zhǎng)6儲(chǔ)集層;測(cè)井解釋模型;孔隙度;滲透率

子長(zhǎng)油田地處鄂爾多斯盆地東部伊陜斜坡中部，其主力產(chǎn)層為長(zhǎng)6油層組，一般厚度為125～135 m。長(zhǎng)6期為三角洲建設(shè)的高潮期，研究區(qū)長(zhǎng)6儲(chǔ)層早期為三角洲前緣沉積，晚期為三角洲平原沉積，沉積物粒度較細(xì)［2］，孔隙度主要分布在6% ～12%之間，滲透率主要分布于 0.1～0.5×10－3μm2以及1～5×10－3μm2之間。屬于復(fù)雜的低孔低滲儲(chǔ)層。利用單一測(cè)井?dāng)?shù)據(jù)對(duì)儲(chǔ)層滲透性進(jìn)行解釋亦或是利用孔隙度與滲透率的相關(guān)關(guān)系計(jì)算滲透率值，在研究區(qū)范圍內(nèi)的精度不夠。而泥質(zhì)含量對(duì)滲透率的影響在該區(qū)表現(xiàn)相當(dāng)明顯。所以將多參數(shù)回歸和泥質(zhì)含量的控制因素結(jié)合才是提高物性解釋精度的關(guān)鍵。

1 儲(chǔ)層“四性”關(guān)系研究

“四性”關(guān)系指的是“儲(chǔ)集層巖性、物性、含油性與電性”之間的關(guān)系。儲(chǔ)層巖性、物性、含油性之間既存在內(nèi)在聯(lián)系又相互制約［3］。圖1是研究區(qū) Z203井的四性關(guān)系圖，由圖可見(jiàn)研究區(qū)的四性關(guān)系如下:

測(cè)井曲線能夠反映不同的巖性。研究區(qū)長(zhǎng)6儲(chǔ)集層砂巖電性特征表現(xiàn)為:自然電位曲線為負(fù)異常，自然伽馬低值，微電極兩條曲線分開(kāi)，聲波時(shí)差曲線相對(duì)較低，而且比較穩(wěn)定。泥巖部位表現(xiàn)為自然電位為基線，自然伽馬高值，微電極兩條曲線重合，聲波時(shí)差曲線相對(duì)較高，且有波動(dòng)，電阻率曲線表現(xiàn)為中 － 高阻［4］。

子長(zhǎng)油田長(zhǎng)6儲(chǔ)層巖性致密。聲波時(shí)差曲線能較好地反映儲(chǔ)層的物性，物性較好的儲(chǔ)層聲波時(shí)差增大。儲(chǔ)層聲波時(shí)差曲線大致為218～254 μs/m。聲波時(shí)差對(duì)儲(chǔ)層孔隙度的反映較為靈敏，但是對(duì)于儲(chǔ)層的滲透性，因?yàn)槠溆绊懸蛩剌^為復(fù)雜，單就聲波時(shí)差難以判斷儲(chǔ)層滲透性的好壞。

圖1 子長(zhǎng)油田Z203井長(zhǎng)6油層組四性關(guān)系圖

油水層的識(shí)別一般用深感應(yīng)曲線。感應(yīng)曲線在油層的電阻率值高于水層，而且油層深、淺感應(yīng)電阻率幅度差小于水層［4］。油層的深感應(yīng)值一般大于13.5Ω·m。

2 定量測(cè)井解釋模型

2.1 泥質(zhì)含量解釋模型

泥質(zhì)含量不僅影響儲(chǔ)層的巖性，而且也影響儲(chǔ)層其他參數(shù)的求取?，F(xiàn)在所使用的測(cè)井系列都不同程度的受泥質(zhì)含量的影響［3，6］。因此，準(zhǔn)確的求取泥質(zhì)含量對(duì)利用測(cè)井?dāng)?shù)據(jù)求取儲(chǔ)層的其他參數(shù)至關(guān)重要［7－10］。在研究區(qū)，自然伽馬曲線能夠靈敏的反映儲(chǔ)層的泥質(zhì)含量，故選用自然伽馬資料來(lái)進(jìn)行定量分析。泥質(zhì)含量計(jì)算公式如下:

式中:c為希爾奇指數(shù)，通常新地層為3.7，老地層為2;GRmin為純砂巖處的自然伽馬值;GRmax

為純泥巖處的自然伽馬值。

根據(jù)上圖可以看出用自然伽馬值所計(jì)算出的泥質(zhì)含量能夠反映出巖心的巖性。通過(guò)對(duì)本區(qū)19口取芯井的巖性與泥質(zhì)含量的對(duì)比研究，認(rèn)為計(jì)算泥質(zhì)含量大于30%的部位在巖心上表現(xiàn)為泥質(zhì)含量較高的儲(chǔ)層，而計(jì)算泥質(zhì)含量小于30%的部位巖心巖性顯示為泥質(zhì)含量較低的砂巖儲(chǔ)層。

2.2 孔隙度解釋模型

此次研究利用了聲波時(shí)差曲線與實(shí)驗(yàn)室物性分析資料來(lái)建立孔隙度模型。為了使解釋模型的精度更高，首先對(duì)研究區(qū)內(nèi)61口井的曲線進(jìn)行了分層，將長(zhǎng)6油層組分為3段，每段的物性有所差異，所以分別建立模型;然后對(duì)所有曲線進(jìn)行了巖心歸位，以消除取芯所造成的巖—電差［10－13］;最后對(duì)聲波時(shí)差采用層點(diǎn)取值的方法，即對(duì)巖性物性較相似的同一層，視為一個(gè)層點(diǎn)進(jìn)行處理［12－15］，這樣既能剔除異常數(shù)據(jù)點(diǎn)消除干擾，又能反映巖性較一致處的整體特征，而不是某一特殊點(diǎn)，該方法更具實(shí)踐性［12］。

圖2 Z111井長(zhǎng)6油層組計(jì)算泥質(zhì)含量與取芯巖性對(duì)比圖

利用聲波時(shí)差與實(shí)測(cè)孔隙度回歸后所得結(jié)果見(jiàn)圖3。

利用上述孔隙度解釋公式，經(jīng)驗(yàn)證計(jì)算孔隙度與分析孔隙度的絕對(duì)誤差的平均值分別為:長(zhǎng)61誤差為0.87%，長(zhǎng)62誤差為0.91%，長(zhǎng) 63誤差為 0.96%，均小于 1%，符合精度要求。

2.3 滲透率解釋模型

考慮到研究區(qū)為低孔低滲儲(chǔ)層，微觀孔隙結(jié)構(gòu)復(fù)雜，僅用孔隙度來(lái)預(yù)測(cè)滲透率的相關(guān)性差［15，18］。說(shuō)明研究區(qū)滲透率的影響因素復(fù)雜，所以在巖心分析、歸位處理、測(cè)井參數(shù)和滲透率影響因素分析的基礎(chǔ)上，選用反映泥質(zhì)含量的自然伽馬相對(duì)值、反映孔隙度相對(duì)大小的聲波時(shí)差值、反映儲(chǔ)層的滲流能力的自然電位幅度差值和滲透率值進(jìn)行回歸分析［4］，期望達(dá)到更高的解釋精度。滲透率與其它參數(shù)的關(guān)系如下:

式中:K－儲(chǔ)集層滲透率，Δt－聲波時(shí)差讀值，ΔGR，ΔSP分別為測(cè)井自然伽馬相對(duì)值與自然電位計(jì)算幅度差值，A、B、C、D 均為常數(shù)。

其中，ΔSP即自然電位幅度差值與自然伽馬相對(duì)值的計(jì)算方法相類(lèi)似，能夠較好的反映儲(chǔ)層的滲流能力。其計(jì)算如下:

圖3 研究區(qū)各亞油層組聲波時(shí)差－孔隙度交會(huì)圖

圖4 研究區(qū)各亞油層組實(shí)測(cè)滲透率－初次計(jì)算滲透率交會(huì)圖

同樣采用層點(diǎn)取值法，對(duì)研究區(qū)38口取心的長(zhǎng)61、62、63亞油層組分別進(jìn)行讀值，并計(jì)算 ΔGR，ΔSP等參數(shù)，最終利用多元線性回歸分析得出各亞油層組的滲透率計(jì)算公式:

根據(jù)上圖可以看出計(jì)算值與分析值其相關(guān)性不好，不能作為滲透率解釋模型直接使用。利用初次解釋模型計(jì)算出的滲透率與分析滲透率做對(duì)比時(shí)發(fā)現(xiàn)，在泥質(zhì)含量較高的部位，計(jì)算滲透率與分析滲透率相差較大，一般情況下計(jì)算滲透率會(huì)大于分析滲透率，而且在這些部位，反映儲(chǔ)層孔隙度的聲波時(shí)差就與滲透率的相關(guān)關(guān)系很差，因此這些部位如果利用初次回歸的公式來(lái)計(jì)算滲透率的話，必然會(huì)造成很不準(zhǔn)確的結(jié)果。

對(duì)于低孔低滲儲(chǔ)層，泥質(zhì)含量的多少影響著儲(chǔ)層的孔隙結(jié)構(gòu)，而孔隙結(jié)構(gòu)的復(fù)雜程度又影響著儲(chǔ)層的滲透性［16，17］，儲(chǔ)層的泥質(zhì)含量如若在某些部位稍有升高就會(huì)造成該部位的細(xì)小喉道的堵塞和使孔隙間的連通性降低，在這些部位的滲透率會(huì)大大減?。?9，20］，因此說(shuō)明泥質(zhì)含量對(duì)滲透率的影響雖是間接但是影響較大，而回歸出的公式雖然添加了能夠反映泥質(zhì)含量的參數(shù)ΔGR，但是要想提高滲透率解釋模型的精度，必須在計(jì)算滲透率時(shí)加強(qiáng)泥質(zhì)含量對(duì)其的控制［14，15，20］。

為了解決這一問(wèn)題并以此來(lái)提高解釋精度的話，采取分類(lèi)的方法來(lái)重新回歸滲透率的計(jì)算公式，運(yùn)用自然伽馬值所計(jì)算的泥質(zhì)含量的大小作為分類(lèi)的標(biāo)準(zhǔn)，將兩種方法計(jì)算的泥質(zhì)含量均小于30%的部位認(rèn)為是砂巖部位，而泥質(zhì)含量大于30%的部位認(rèn)為是泥巖部位，分別進(jìn)行回歸。

圖5 研究區(qū)長(zhǎng)6油層組實(shí)測(cè)滲透率分類(lèi)計(jì)算滲透率交會(huì)圖

分類(lèi)解釋的結(jié)果如下:

長(zhǎng)61亞油層組:泥質(zhì)含量小于30%的部位，二次分類(lèi)回歸公式如下:

泥質(zhì)含量大于30%的部位，二次分類(lèi)回歸公式如下:

LgK=0.032× Δt－0.075 × ΔSP－0.216× ΔGR －7.687 R=0.697

二次分類(lèi)之后，再計(jì)算滲透率與實(shí)測(cè)滲透率的相關(guān)關(guān)系圖(圖5A):

長(zhǎng)62亞油層組:泥質(zhì)含量小于30%的部位，二次分類(lèi)回歸公式如下:

泥質(zhì)含量大于30%的部位，二次分類(lèi)回歸公式如下:

二次分類(lèi)之后，再計(jì)算滲透率與實(shí)測(cè)滲透率的相關(guān)關(guān)系圖(圖5B):

長(zhǎng)63亞油層組:泥質(zhì)含量小于30%的部位，二次分類(lèi)回歸公式如下:

泥質(zhì)含量大于30%的部位，二次分類(lèi)回歸公式如下:

二次分類(lèi)之后，再計(jì)算滲透率與實(shí)測(cè)滲透率的相關(guān)關(guān)系圖(圖5C):

根據(jù)圖5可以看出，二次分類(lèi)之后，利用分類(lèi)后的數(shù)據(jù)回歸的公式算出的滲透率與實(shí)測(cè)滲透率的相關(guān)性很好，滿足研究區(qū)解釋精度的要求，可以作為研究區(qū)的滲透率測(cè)井解釋模型。

3 結(jié)論

1)利用自然伽馬值和自然電位幅度值計(jì)算出的泥質(zhì)含量能夠?qū)r性有較好地反映，并通過(guò)對(duì)比得出儲(chǔ)層砂巖的泥質(zhì)含量下限值。

2)利用聲波時(shí)差值回歸出了符合精度要求的孔隙度的計(jì)算模型。

3)先根據(jù)泥質(zhì)含量大小對(duì)儲(chǔ)層的滲透率進(jìn)行分類(lèi)，分類(lèi)后運(yùn)用聲波時(shí)差值，自然伽馬差值及自然電位的幅度差值進(jìn)行回歸，分類(lèi)后的計(jì)算模型明顯比未分類(lèi)前更加接近實(shí)際滲透率值，且達(dá)到了解釋的精度。

［1］趙軍，陳福宣，燕軍.巖心刻度法在測(cè)井解釋模型中的應(yīng)用［J］.新疆石油地質(zhì).1998，19(05):421－423.

［2］郭蘭，丁超，閆繼福.鄂爾多斯盆地子長(zhǎng)油田安定區(qū)延長(zhǎng)組長(zhǎng)6段沉積微相［J］.油氣地質(zhì)與采收率.2009，16(03):11－14.

［3］譚成仟.鄂爾多斯盆地白豹地區(qū)長(zhǎng)6油藏地質(zhì)模型研究［D］.西安:西北大學(xué)，2008:63－71.

［4］張小莉.延長(zhǎng)油區(qū)中生界測(cè)井系列優(yōu)選及解釋技術(shù)研究［R］.西安:陜西延長(zhǎng)石油(集團(tuán))有限責(zé)任公司.2010:115－146.

［5］杜奉屏.油礦地球物理測(cè)井［M］.北京:地質(zhì)出版社.1984:81－110.

［6］歐陽(yáng)健，王貴文，吳繼余等.測(cè)井地質(zhì)分析與油氣層定量評(píng)價(jià)［M］.北京:石油工業(yè)出版社.1999:142－170.

［7］朱紅濤，黃眾，劉浩冉等.利用測(cè)井資料識(shí)別層序地層單元技術(shù)與方法進(jìn)展及趨勢(shì)［J］.地質(zhì)科技情報(bào).2011，30(04):29－36.

［8］任培罡，夏存銀，李媛等.自組織神經(jīng)網(wǎng)絡(luò)在測(cè)井儲(chǔ)層評(píng)價(jià)中的應(yīng)用［J］.地質(zhì)科技情報(bào).2010，29(03):114－117.

［9］劉秀娟，陳超，曾沖等.利用測(cè)井?dāng)?shù)據(jù)進(jìn)行巖性識(shí)別的多元統(tǒng)計(jì)方法［J］.地質(zhì)科技情報(bào).2007，26(03):109－112.

［10］趙彥超，吳春萍，吳東平.致密砂巖氣層的測(cè)井評(píng)價(jià)—以鄂爾多斯盆地大牛地山西組一段氣田為例［J］.地質(zhì)科技情報(bào).2003，22(04):65－70.

［11］趙俊英，羅靜蘭，雷曉嵐等.鄂爾多斯盆地子長(zhǎng)油田楊家園則探區(qū)長(zhǎng) 6物源分析［J］.中國(guó)地質(zhì).2007，34(03):422－427.

［12］高鵬，高勝利，汶鋒剛等.鄂爾多斯盆地子長(zhǎng)油田長(zhǎng)6油層組非均質(zhì)性研究［J］.天然氣勘探與開(kāi)發(fā).2008，31(03):14－18.

［13］王文明.測(cè)井參數(shù)定量化技術(shù)在榆樹(shù)林油田葡萄花儲(chǔ)層物性解釋中的應(yīng)用［J］.石油與天然氣地質(zhì).2009，30(02):230－235.

［14］高潤(rùn)峰.鄂爾多斯盆地南部鎮(zhèn)涇油田電測(cè)井資料在油層物性解釋中的應(yīng)用［J］.科技情報(bào)開(kāi)發(fā)與經(jīng)濟(jì).2007，17(21):176－177.

［15］王志章，熊琦華，張一偉等.非均質(zhì)復(fù)雜斷塊油田儲(chǔ)層參數(shù)測(cè)井解釋模型［J］.石油學(xué)報(bào).1994，15(S1):89－97.

［16］Chunming Xu，Thomas P. Cronin，Thalbert E. McGinness，et al.Middle Atokan sediment gravity flows in the Red Oak field，Arkoma Basin，Oklahoma:A sedimentary analysis using electrical borehole images and wireline logs［J］. AAPG Bulletin，2009，93:1 － 29.

［17］Flavio S. Anselmetti，Gregor P. Eberli，The velocity － deviation log;a tool to predict pore type and permeability trends in carbonate drill holes from sonic and porosity or density logs［J］. AAPG Bulletin，1999，83:450－466.

［18］Jordan G. Mimoun，Carlos Torres － Verdín，William E. Preeg.Quantitative interpretation of pulsed neutron capture logs［J］. Geophysics，2011，76:E81－E103.

［19］Romain Prioul，Adam Donald，Randy Koepsell，et al. Forward modeling of fracture－induced sonic anisotropy using a combination of borehole image and sonic logs［J］. Geophysics，2007，72:E135－E147.

［20］L.de Witte，M.R.J.Wyllie，Log interpretation in sandstone reservoirs;discussion and reply［J］. Geophysics，1961，26:101 －102.

Log Interpretation Model of Physical Property of Chang 6 Reservoir，Zichang Oilfield

LIU Jing－jing，ZHANG Xiao－li
(Department of Geology，Northwest University，Xi’an 710069，Shaanxi)

The interpretation models of physical parameters are established to enhance the accuracy of log interpretation of Zichang oilfield，which has important implication for the formation assessment and reservoir description. Firstly，it briefly introduces the four property relations by using the logging data，core physical property analysis and core description. Then，the interpretation models of shale content and porosity of reservoir are established by using the layer－point reading data method based on the first step. The permeability interpretation is classified by shale content and then regression analyzed by multiparameters，so that the accuracy of permeability interpretation is enhanced. The model validation test has been shown with high accuracy and satisfactory result.

Zichang oilfield;Chang 6 reservoir;log interpretation model;porosity and permeability

P631.8+19

A

1004－1184(2012)05－0081－03

2012－05－21

“十一五”國(guó)家科技支撐計(jì)劃資助項(xiàng)目(2007BAB17B01)

劉晶晶(1987－)，女，江西蓮花人，在讀碩士研究生，主攻方向:測(cè)井資料處理與解釋方面的研究。