， ，， ，，

(1.山東科技大學(xué) 地球科學(xué)與工程學(xué)院，山東 青島 266590；2.中國石化江蘇油田分公司 第二采油廠，江蘇 淮安 223000)

洛倫茲系數(shù)在儲層非均質(zhì)性評價中的應(yīng)用

鄭晨晨1，謝俊1，王金凱1，趙璇1，段雅君1，孫燕2

(1.山東科技大學(xué) 地球科學(xué)與工程學(xué)院，山東 青島 266590；2.中國石化江蘇油田分公司 第二采油廠，江蘇 淮安 223000)

針對目前油氣儲層宏觀非均質(zhì)性評價參數(shù)眾多、評價指標(biāo)無邊界的問題，以江蘇金湖凹陷王龍莊油田為例，綜合運用巖心觀察、室內(nèi)分析試驗和測井解釋等資料，通過計算巖石樣品的滲透率洛倫茲系數(shù)對研究區(qū)阜寧組二段儲層層內(nèi)和平面非均質(zhì)性進行研究。洛倫茲系數(shù)為0～1的數(shù)值，可以作為評價儲層非均質(zhì)程度的指標(biāo)，從而定量評價儲層的宏觀非均質(zhì)性。結(jié)果表明，洛倫茲系數(shù)能夠彌補傳統(tǒng)儲層非均質(zhì)性評價指標(biāo)無邊界的問題，可在其他地區(qū)的儲層宏觀非均質(zhì)性評價中進行推廣。

非均質(zhì)性；洛倫茲系數(shù)；低滲透儲層；王龍莊油田

儲層非均質(zhì)性研究是儲層評價的核心內(nèi)容，也是影響油氣水分布及注水開發(fā)的重要因素，在低滲透油田中尤其重要[1]。目前，碎屑巖儲層非均質(zhì)性研究在層序地層學(xué)、地質(zhì)統(tǒng)計學(xué)、地質(zhì)建模等方面取得了一些進展，但主要還是通過計算滲透率變異系數(shù)Vk、突進系數(shù)Tk、級差Jk等參數(shù)進行表征[2]。采用上述參數(shù)從不同角度對儲層的層內(nèi)非均質(zhì)性進行評價，但這些數(shù)值的計算在理論上都是無界的，不能對儲層的非均質(zhì)程度進行定量評價。目前的解決辦法是將這些數(shù)值定性分成三個區(qū)間(弱非均質(zhì)性、中等非均質(zhì)性及嚴(yán)重非均質(zhì)性)，從而大致描述儲層的層內(nèi)非均質(zhì)性級別[3-5]。所以傳統(tǒng)的手段無法定量描述儲層的非均質(zhì)程度。

洛倫茲系數(shù)又稱基尼系數(shù)，最早由意大利經(jīng)濟學(xué)家基尼[5-6]提出，原是經(jīng)濟學(xué)判斷收入分配是否均衡的指標(biāo)，后被石油地質(zhì)學(xué)家用來定量描述儲層非均質(zhì)性。如岳崇旺等[2]利用洛倫茲系數(shù)識別和評價儲層非均質(zhì)性；Peter等[7]將其引入碳酸鹽巖儲層的非均質(zhì)性評價中；劉超等[8]對洛倫茲系數(shù)評價儲層非均質(zhì)性的方法進行改進；王慶等[9]利用洛倫茲系數(shù)研究油藏的產(chǎn)液、吸水剖面。該方法多利用測井曲線計算洛倫茲系數(shù)。本次研究結(jié)合取心井巖心樣品數(shù)據(jù)和常規(guī)測井資料，計算滲透率洛倫茲系數(shù)來評價儲層的非均質(zhì)程度，實現(xiàn)對儲層非均質(zhì)性的定量評價。

1 地質(zhì)概況

王龍莊油田位于蘇北盆地金湖凹陷汊澗次凹東側(cè)，物源來自西側(cè)的張八嶺隆起[10]，是一系列斷鼻構(gòu)造組成的斷塊群(圖1)。研究區(qū)王龍莊斷塊位于王龍莊油田西南方向，該區(qū)主要含油層位為阜寧組二段(以下簡稱阜二段)，地質(zhì)儲量116×104t。按研究區(qū)原有油礦地質(zhì)學(xué)劃分，阜二段自上而下共分為Ⅰ、Ⅱ、Ⅲ三個砂組，Ⅱ、Ⅲ砂組為主要研究層位，Ⅱ砂組發(fā)育2個含油小層，Ⅲ砂組發(fā)育7個含油小層。王龍莊斷塊阜二段儲層平均滲透率為20.13×10-3μm2，按照碎屑巖儲層分類標(biāo)準(zhǔn)[11]屬于低滲透儲層，本次主要應(yīng)用洛倫茲系數(shù)法研究低滲透儲層非均質(zhì)性。

圖1 金湖凹陷王龍莊斷塊區(qū)域位置圖Fig.1 Location map of Wanglongzhuang-fault-block in Jinhu Depression

2 洛倫茲系數(shù)計算原理

洛倫茲系數(shù)法是采用儲層巖石樣品的實際滲透率貢獻值曲線和理想儲層絕對均質(zhì)曲線之間的面積與儲層絕對均質(zhì)曲線和縱坐標(biāo)軸圍成面積的比值來表示儲層的非均質(zhì)程度[12]。因此，繪制洛倫茲曲線的前提是求取巖樣滲透率貢獻值累積百分?jǐn)?shù)，計算步驟如下：

1) 將研究區(qū)儲層的滲透率值由大到小按降序排列，保證計算滲透率貢獻值累積百分?jǐn)?shù)從大到小開始累積。

2) 利用公式(1)、(2)計算每一塊巖石樣品的滲透率貢獻值累積百分?jǐn)?shù)。

(1)

其中：Li為第i個巖石樣品所對應(yīng)的滲透率值對整體滲透率，%；Ki為第i個巖石樣品所對應(yīng)的滲透率值，10-3μm2；n為此次研究所用巖石樣品的總個數(shù)。

(2)

其中：Yi為第i個巖石樣品所對應(yīng)的滲透率貢獻值累積百分?jǐn)?shù)，%。

3) 同理，可以利用公式(3)和(4)求出每一塊巖石樣品的巖樣累積百分?jǐn)?shù)。

(3)

(4)

其中，Qi為第i個巖石樣品對整體的貢獻值，%；Pi為第i個巖石樣品；n為此次研究所用巖石樣品的總個數(shù)；Xi為第i個巖石樣品所對應(yīng)的巖樣累積百分?jǐn)?shù)，%。

圖2 滲透率洛倫茲曲線Fig. 2 Lorentz curve of permeability

將求出的巖樣滲透率貢獻值累積百分?jǐn)?shù)和巖樣累積百分?jǐn)?shù)(Xi，Yi)放在直角坐標(biāo)系中得到相應(yīng)儲層的洛倫茲曲線(圖2)。對于完全均質(zhì)的儲層，洛倫茲曲線應(yīng)為一條直線Y=X(圖2，AC)，表示儲層中每一個數(shù)據(jù)點的滲透率值相等；而對于極端非均質(zhì)的儲層，洛倫茲曲線為一條折線(圖2，ABC)，表示儲層滲透率差異非常大；通常情況下，儲層滲透率洛倫茲曲線分布于兩者之間，為一條呈上凸形態(tài)的曲線(圖2，L)。

設(shè)洛倫茲曲線L與儲層絕對均質(zhì)曲線AC所圍成的面積為S1，絕對均質(zhì)曲線AC與極端非均質(zhì)曲線ABC所圍成的面積為(S1+S2)，則洛倫茲系數(shù)即為S1/(S1+S2)。其中，完全均質(zhì)儲層的洛倫茲系數(shù)為0，極端非均質(zhì)儲層的洛倫茲系數(shù)為1，常規(guī)儲層滲透率洛倫茲系數(shù)在0～1之間。因此，洛倫茲系數(shù)可以作為儲層的非均質(zhì)程度系數(shù)，以此定量評價儲層的宏觀非均質(zhì)程度。

3 儲層非均質(zhì)性評價

根據(jù)研究目的、儲層非均質(zhì)性可有多種分類方法，其中裘亦楠(1992)[13]將碎屑巖的儲層非均質(zhì)性劃分為層內(nèi)、層間、平面和孔隙非均質(zhì)性四類，其中儲層宏觀非均質(zhì)性研究主要從層內(nèi)、層間以及平面三方面著手。沉積因素是影響儲層非均質(zhì)性的最根本因素，因不同沉積環(huán)境沉積時水體能量、沉積物的性質(zhì)(粒度、成分、分選性、磨圓度)等不同，導(dǎo)致儲層內(nèi)部沉積結(jié)構(gòu)和構(gòu)造的各向異性，從而造成儲層非均質(zhì)程度的差別[14]。本研究從沉積微相劃分入手，利用洛倫茲系數(shù)評價儲層的層內(nèi)和平面非均質(zhì)性。

3.1 層內(nèi)非均質(zhì)性評價

儲層的層內(nèi)非均質(zhì)性是指各小層內(nèi)的儲層特征變換，分別運用傳統(tǒng)的變異系數(shù)法和洛倫茲系數(shù)法對研究區(qū)儲層的層內(nèi)非均質(zhì)性進行評價。

3.1.1 變異系數(shù)法

根據(jù)王龍莊斷塊33口井的測井?dāng)?shù)據(jù)，結(jié)合6口取心井的巖心資料綜合分析，建立回歸公式，對研究區(qū)各井的孔滲數(shù)據(jù)進行校正，分別計算每口井阜二段Ⅱ、Ⅲ砂組中7個主要含油層系的變異系數(shù)、突進系數(shù)等數(shù)值(表1)。

表1 王龍莊斷塊阜二段層內(nèi)變異系數(shù)統(tǒng)計Tab.1 Statistics of variation coefficient in the E1f2 of Wanglongzhuang-fault-block

3.1.2 洛倫茲系數(shù)法

表2 王龍莊斷塊阜二段層內(nèi)洛倫茲系數(shù)統(tǒng)計Tab.2 Statistics of Lorentz coefficient in the E1f2 of Wanglongzhuang-fault-block

表3 王龍莊斷塊儲層非均質(zhì)性洛倫茲系數(shù)評價標(biāo)準(zhǔn)Tab. 3 Reservoir heterogeneity evaluation standard of Lorentz coefficient in Wanglongzhuang-fault-block

由表2、表3可以看出，研究區(qū)各小層洛倫茲系數(shù)0.27～0.51，與變異系數(shù)呈正比，屬嚴(yán)重非均質(zhì)性。洛倫茲系數(shù)取值在0(完全均質(zhì))到1(極端非均質(zhì))之間，可以作為儲層非均質(zhì)程度系數(shù)來描述儲層的非均質(zhì)程度。

圖3 王龍莊斷塊滲透率洛倫茲曲線Fig.3 Permeability Lorentz curve of Wanglongzhuang-fault-block

此外，不同沉積環(huán)境控制下儲層的韻律性不同，而不同的韻律性將導(dǎo)致儲層物性縱向上的差異[15]。研究區(qū)主要為濱淺湖灘壩沉積，由于沉積過程受湖浪的頻繁改造，具有反韻律、正韻律及復(fù)合韻律多種韻律模式。通過單井縱向測井曲線、滲透率分布，研究區(qū)主要以反韻律為主，細砂巖或粉砂巖與泥巖頻繁互層，層內(nèi)滲透率變化較大、導(dǎo)致層內(nèi)非均質(zhì)性嚴(yán)重(圖4)。

兩種方法對比與應(yīng)用可以發(fā)現(xiàn)，相比于變異系數(shù)法，洛倫茲系數(shù)法具有以下優(yōu)點：①數(shù)值計算有界，能夠作為儲層非均質(zhì)程度系數(shù)使用，可以定量評價儲層的非均質(zhì)性；②可以克服傳統(tǒng)方法各參數(shù)表征角度各異的問題，作為統(tǒng)一評價非均質(zhì)程度的標(biāo)準(zhǔn)；③結(jié)合洛倫茲曲線使得結(jié)論更加簡潔、直觀。

3.2 平面非均質(zhì)性評價

平面非均質(zhì)性是指各小層滲透率在平面上的變化情況，由于不同沉積微相在平面上的分布不同所表現(xiàn)出的不同物性特征，可以從沉積微相平面分布、滲透率平面分布等反映出來[16]。目前廣泛采用的滲透率等值線圖雖然可以直觀展示平面非均質(zhì)性變化的特點，但依然難以定量評價平面非均質(zhì)程度，結(jié)合洛倫茲系數(shù)可以解決這一難題。

圖4 王龍莊斷塊阜二段Ⅲ砂組層內(nèi)非均質(zhì)性定量表征(W10井)Fig. 4 Quantitative characterization of the inlayer heterogeneity of the 3rd sand sets in the second member of Funing formation in Wanglongzhuang-fault-block (well W10)

3.3 結(jié)果評價及措施

圖5 王龍莊斷塊阜二段小層平面非均質(zhì)性定量表征Fig.5 Quantitative characterization of the plane heterogeneity of E1 and E1 layers in Wanglongzhuang-fault-block

4 結(jié)論

1) 洛倫茲系數(shù)是在儲層滲透率累積百分?jǐn)?shù)基礎(chǔ)上計算出的一種新參數(shù)，與變異系數(shù)法相比具有定量化、結(jié)論直觀和計算簡便等優(yōu)點；洛倫茲系數(shù)的引入可以完善目前儲層非均質(zhì)性評價體系，為儲層非均質(zhì)程度提供了統(tǒng)一的標(biāo)準(zhǔn)，改變了以往評價方法中缺乏統(tǒng)一評價指標(biāo)的問題。

2) 以王龍莊斷塊阜二段儲層滲透率參數(shù)為基礎(chǔ)，結(jié)合儲層層內(nèi)韻律性變化，利用洛倫茲系數(shù)法對研究區(qū)儲層層內(nèi)非均質(zhì)程度進行評價，結(jié)果顯示研究區(qū)儲層洛倫茲系數(shù)在0.27～0.51，非均質(zhì)程度嚴(yán)重；研究區(qū)發(fā)育濱淺湖灘壩沉積，沙灘砂體物性較差但分布廣泛，沙壩砂體物性較好但零星分布，導(dǎo)致儲層平面非均質(zhì)程度嚴(yán)重。

3) 洛倫茲系數(shù)法在王龍莊斷塊儲層非均質(zhì)性評價中效果較好，評價結(jié)果與現(xiàn)場生產(chǎn)實踐相吻合，可以在其他類似油田推廣應(yīng)用。

[1]張金亮,謝俊.油田開發(fā)地質(zhì)學(xué)[M].北京：石油工業(yè)出版社,2011:142-153.

[2]岳崇旺,楊小明,鐘曉勤,等.利用測井曲線的洛倫茲系數(shù)評價地層的非均質(zhì)性[J].吉林大學(xué)學(xué)報(地球科學(xué)版),2015,45(5):1539-1546.

YUE Chongwang,YANG Xiaoming,ZHONG Xiaoqin,et al.Evaluation of formation heterogeneity using Lorentz coefficient of logging curves[J].Journal of Jilin University (Earth Science Edition),2015,45(5):1539-1546.

[3]嚴(yán)科,楊少春,任懷強.儲層宏觀非均質(zhì)性定量表征研究[J].石油學(xué)報,2008,29(6):870-874.

YAN Ke,YANG Shaochun,REN Huaiqiang.Research on quantitative characterization of macroscopic heterogeneity of reservoir[J].Acta Petrolei Sinica,2008,29(6):870-874.

[4]羅曉容,王忠楠,雷裕紅,等.特超低滲砂巖油藏儲層非均質(zhì)性特征與成藏模式:以鄂爾多斯盆地西部延長組下組合為例[J].石油學(xué)報,2016,37(增1):87-98.

LUO Xiaorong,WANG Zhongnan,LEI Yuhong,et al.Heterogeneity characteristics and accumulation model of ultra-low permeability sandstone reservoirs:A case study of the lower part of Yanchang formation in the western Ordos Basin,China[J].Acta Petrolei Sinica,2016,37(S1):87-98.

[5]游士兵,王原君.研究收入分配問題的一種新洛倫茲曲線模型:建構(gòu)與應(yīng)用[J].經(jīng)濟評論,2014(2):3-15.

YOU Shibing,WANG Yuanjun.A new Lorenz model for income distribution analysis:Model construction and application[J].Economic Review,2014,(2):3-15.

[6]何邦強,洪興健.基尼系數(shù)計算與分解方法研究綜述[J].統(tǒng)計與決策,2016,14:13-17.

HE Bangqiang,HONG Xingjian.A review of Gini coefficient calculation and decomposition[J].Statistics & Decision,2016,14:13-17.

[7]FITCH P,DAVIES S,LOVELL M,et al.Heterogeneity in carbonate petrophysical properties:Application to fluid flow units and sampling strategies[C]//SPWLA 51st Annual Logging Symposium.Xi’an:CNPC Logging,2010:19-23.

[8]劉超,馬奎前,李紅英,等.基于洛倫茲曲線法定量表征儲層非均質(zhì)性方法的改進與應(yīng)用[J].中國海上油氣,2012,24(2):36-38.

LIU Chao,MA Kuiqiang,LI Hongying，et al.An improvement of quantitative reservoir-heterogeneity characterization based on Lorenz curve method and its application[J].China Offshore Oil and Gas,2012,24(2):36-38.

[9]王慶,劉慧卿,殷方好.洛倫茲曲線在油藏產(chǎn)液、吸水剖面研究中的應(yīng)用[J].特種油氣藏,2010,17(1):71-74.

WANG Qing,LIU Huiqing,YIN Fanghao.Application of Lorenz curve to the study of production profile and injection profile[J].Special Oil and Gas Reservoirs,2010,17(1):71-74.

[10]丁玉盛,周巨標(biāo),石先達,等.王龍莊油田阜寧組二段碳酸鹽與碎屑巖混合沉積灘壩體系[J].山東科技大學(xué)學(xué)報(自然科學(xué)版),2015,34(2):58-63.

DING Yusheng,ZHOU Jubiao,SHI Xianda,et al.Beach bar system mixed carbonate and siliciclastic sediments in the second member of Funing formation in Wanglongzhuang oilfield[J].Journal of Shandong University of Science and Techno-logy(Natural Science),2015,34(2):58-63.

[11]中華人民共和國地質(zhì)礦產(chǎn)行業(yè)標(biāo)準(zhǔn).DZ/T0217-2005,石油天然氣儲量計算規(guī)范[S].

[12]馬玉玲.勞倫茲方法在儲層非均質(zhì)評價中的應(yīng)用[J].長江大學(xué)學(xué)報(自然科學(xué)版),2013,10(10):90-91.

MA Yuling.Application of Lorenz method in the evaluation of reservoir heterogeneity[J].Journal of Yangtze University(Nature Science Edition),2013,10(10):90-91.

[13]裘亦楠.中國陸相碎屑巖儲層沉積學(xué)的進展[J].沉積學(xué)報,1992,10(3):16-24.

QIU Yinan.Developments in reservoir sedimentology of continental clastic rocks in China[J].Acta Sedimentologica Sinica,1992,10(3):16-24.

[14]孟寧寧,劉懷山,張金亮,等.黑46斷塊低滲油藏非均質(zhì)性研究及其對剩余油分布的影響[J].海洋地質(zhì)與第四紀(jì)地質(zhì),2016,36(1):143-150.

MENG Ningning,LIU Huaishan,ZHANG Jinliang,et al.Heterogeneity and its impact on remaining oil distribution of low-permeability reservoirs in the HEI46-Block[J].Marine Geology & Quaternary Geology,2016,36(1):143-150.

[15]中華人民共和國石油天然氣行業(yè)標(biāo)準(zhǔn).SY/T5717-95,單井碎屑巖儲層評價[S].

[16]張世廣,盧雙舫,張雁,等.高分辨率層序地層學(xué)在儲層宏觀非均質(zhì)性研究中的應(yīng)用:以松遼盆地朝陽溝油田朝1-朝氣3區(qū)塊扶余油層為例[J].沉積學(xué)報,2009,27(3):458-469.

ZHANG Shiguang,LU Shuangfang,ZHANG Yan,et al.Application of high-resolution sequence stratigraphy in the study of the reservoir macro-heterogeneity:Taking the Fuyu oil layer of Chao1-Chaoqi3-Blocks in Chaoyanggou oil reservoir of Songliao basin as an example[J].Acta Sedimentologica Sinica,2009,27(3):458-469.

ApplicationofLorentzCoefficientinReservoirHeterogeneityEvaluation

ZHENG Chenchen1, XIE Jun1, WANG Jinkai1, ZHAO Xuan1, DUAN Yajun1, SUN Yan2

(1. College of Earth Science and Engineering, Shandong University of Science and Technology, Qingdao, Shandong 266590, China; 2. Second Oil Production Plant of Jiangsu Oilfield Company, SINOPEC, Huai′an, Jiangsu 223000, China)

There are some defects in the existing evaluation systems of macroscopic heterogeneity such as the great number and unboundedness of evaluation parameters. Taking Wanglongzhuang Oilfield in Jinhu Depression, Jiangsu province as an example, this paper studied the reservoir layers and plane heterogeneity of the second number of Funing formation through the calculation of Lorentz coefficient of rock samples’ permeability in the study area based on core observation, laboratory analysis and logging interpretation. Lorentz coefficient value 0-1 can be used as the evaluation index of reservoir heterogeneity degree so as to make quantitative evaluations of the reservoir’s macro-scopic heterogeneity. Application results show that Lorentz coefficient can make up for the unboundedness of traditional reservoir hetero-geneity evaluation indexes and can be extended to the evaluation of reservoir’s macroscopic heterogeneity in other places.

heterogeneity； Lorentz coefficient; low permeability reservoir； Wanglongzhuang Oilfield

鄭晨晨,謝俊,王金凱,等.洛倫茲系數(shù)在儲層非均質(zhì)性評價中的應(yīng)用[J].山東科技大學(xué)學(xué)報(自然科學(xué)版),2018,37(1):103-110.

2016-11-29

國家自然科學(xué)基金項目(51674156，51504143)

鄭晨晨(1992—)，男，山東濟寧人，碩士研究生，主要從事油氣田開發(fā)方面的研究.E-mail：441760455@qq.com

謝 俊(1968—)，男，湖北京山人，教授，博士生導(dǎo)師，主要從事油氣儲層評價與油藏描述方面的研究，本文通信作者.E-mail：xiejun0532@163.com

ZHENG Chenchen, XIE Jun, WANG Jinkai,et al. Application of Lorentz coefficient in reservoir heterogeneity evaluation[J].Journal of Shandong University of Science and Technology(Natural Science),2018,37(1):103-110.

TE348

A

1672-3767(2018)01-0103-08

10.16452/j.cnki.sdkjzk.2018.01.010

呂海亮)