針對海上非均質(zhì)儲層提液措施的歷史擬合策略研究

姜瑞忠，喬杰

(中國石油大學(xué)(華東) 石油工程學(xué)院，山東 青島 266555)

孫輝

(中國石油大學(xué)(華東) 計算機與通信工程學(xué)院，山東 青島 266555)

彭國強

(中國石油大學(xué)(華東) 石油工程學(xué)院，山東 青島 266555)

陳冠中，林春陽，李強

(海洋石油高效開發(fā)國家重點實驗室；中海油研究總院，北京 100028)

由于在實際生產(chǎn)中，低滲透層存在啟動壓力梯度，但是，目前大多數(shù)較為成熟的商業(yè)軟件在針對嚴重非均質(zhì)地層的模擬過程中，都是基于達西定律研制開發(fā)的[1]；加之在多小層的產(chǎn)量劈分上，現(xiàn)場在沒有產(chǎn)液剖面測試資料的情況下，常用的方法按照地層系數(shù)(即滲透率與厚度之積)進行。而海上油田，往往會通過產(chǎn)液量大幅度改變來達到一定的開發(fā)目的。因此，在商業(yè)化軟件中，必定會導(dǎo)致歷史擬合工作很難進行下去。為此，筆者將討論如何在Eclispe軟件中解決非均質(zhì)儲層產(chǎn)液量大幅度變化引起的歷史擬合困難問題。

1 非均質(zhì)儲層中的非線性滲流模型研究現(xiàn)狀

圖1 低滲透儲層壓力梯度與流體速度示意圖

隨著海上區(qū)塊實際生產(chǎn)過程中的非均質(zhì)問題的日益突出，越來越多的研究學(xué)者[2～8]開始注意非均質(zhì)儲層開發(fā)過程中的非線性滲流問題。在非均質(zhì)的多孔介質(zhì)中，低滲透層發(fā)生的流動往往為非線性流動。壓力梯度與流體速度的關(guān)系如圖1所示。

由于商業(yè)軟件Eclispe中的基本流動理論都是基于達西定律研制開發(fā)的，所以針對這種流動，往往是很難實現(xiàn)的。即使能實現(xiàn)，也是采用了一種近似處理的方法，設(shè)置“門限壓力”。其原理就是將儲層設(shè)置n個平衡分區(qū)，然后運用關(guān)鍵字“THPRES”設(shè)置每一個平衡分區(qū)與相鄰平衡分區(qū)之間的門限壓力，如果將這種方法應(yīng)用于現(xiàn)場實際模型，不僅繁瑣而且難度較大。

2 提液/降液生產(chǎn)提高采收率機理

在實際的非均質(zhì)儲層的開發(fā)生產(chǎn)過程中，現(xiàn)場經(jīng)常會采用提液或者換小泵降低產(chǎn)液量生產(chǎn)，以求達到提高采收率的目的[9～14]。

提液生產(chǎn)能提高采收率的機理是：在縱向上非均質(zhì)嚴重的儲層中，開發(fā)初期的低滲透層由于啟動壓力的存在，并不會真正的投入生產(chǎn)。當進行提液操作之后，產(chǎn)液量的大幅度上升必會導(dǎo)致井底壓力的迅速下降，當?shù)蜐B層位內(nèi)部的地層壓力p0與井底流壓的差值大于啟動壓力的時候，低滲層位便開始有流體產(chǎn)出，達到了動用低滲儲層、提高采收率的目的(圖2)。

換小泵生產(chǎn)，其主要機理是在實際的儲層生產(chǎn)中，存在低滲透層開發(fā)效果變差、含水率升高、高滲透層并沒有完全動用的現(xiàn)象。當該現(xiàn)象出現(xiàn)時，現(xiàn)場就會采用換小泵生產(chǎn)，降低產(chǎn)液速率。其機理在于，當降液生產(chǎn)的時候，產(chǎn)油井的井底流壓迅速提高，低滲透層的生產(chǎn)壓差變小，當生產(chǎn)壓差低于該低滲透層的啟動壓力梯度時，低滲透層便會停止生產(chǎn)，達到降低含水率的目的(圖3)。

圖2 提液生產(chǎn)提高采收率機理

圖3 降液生產(chǎn)改善開采效果機理

3 等效模擬方法

在進行相關(guān)操作之前，首先需要對實際的地質(zhì)模型和生產(chǎn)動態(tài)進行仔細的分析與研究，在該基礎(chǔ)上，運用Eclispe中的WPIMULT關(guān)鍵字[15,16]，對每一個射孔層位進行設(shè)置。

3.1 WPIMULT描述

WPIMULT的主要作用是在模型內(nèi)部可以分別針對某一層的傳導(dǎo)率乘積因子進行設(shè)置。當設(shè)置的乘積因子大于1的時候，表示該小層的滲流效果變好；反之，當乘積因子小于1時，則表示該小層的滲流效果變差。運用WPIMULT的這一特殊功能，可以在一定程度上模擬出因為提液生產(chǎn)或者降液生產(chǎn)帶來的低滲透層、高滲透層的滲流效果變化情況。

3.2 產(chǎn)液量變化設(shè)置描述

在實際生產(chǎn)過程中進行提液操作時，會導(dǎo)致井底流壓迅速下降，生產(chǎn)壓差變大，當生產(chǎn)壓差大于啟動壓力時，低滲透儲層開始有流體產(chǎn)出，有效地動用了低滲透儲層中的油層。

在實際生產(chǎn)過程中進行降液生產(chǎn)時，井底流壓瞬間增大，生產(chǎn)壓差變小，當生產(chǎn)壓差小于啟動壓力時,低滲透層便會停止生產(chǎn)，在現(xiàn)場這往往是為了解除由于低滲層的不良生產(chǎn)情況導(dǎo)致的含水率過高的問題。

3.3 WPIMULT的優(yōu)勢對比

在一般的概念模型中，根據(jù)生產(chǎn)機理，經(jīng)常還可以用WELOPEN關(guān)鍵字設(shè)置某一小層和某幾小層的開、關(guān)情況，這在機理上完全可以[17，18]。但是，在現(xiàn)場實際模型當中，有兩方面因素限制了該關(guān)鍵字的應(yīng)用：其一，由于地質(zhì)模型內(nèi)部的網(wǎng)格劃分和隨機建模的原因，某一油層內(nèi)部可能既包含高滲透層，又包含低滲透層，因此，并不能僅僅用WELOPEN將這一層關(guān)閉，而當采用這個WPIMULT時，就完全可以設(shè)置一個介于0～1之間的乘積因子，來達到完全模擬的效果；其二，在進行降液生產(chǎn)時，可能存在高滲透層的產(chǎn)能提高現(xiàn)象，而僅僅采用WELOPEN，并不能將該效果體現(xiàn)出來。

4 應(yīng)用實例分析

以渤海SZ36-1油田中的J02井為例。該油田分布廣，埋藏淺，層系多，油層厚，縱向上呈現(xiàn)嚴重的非均質(zhì)特性。J02井為一口后期轉(zhuǎn)注井，射開層位為23、24、26、27、29、30、34、36小層。

圖4為該油井的生產(chǎn)歷史曲線。在2005～2006年之間，現(xiàn)場進行了換小泵大幅度降液生產(chǎn)操作。顯然，含水率開始下降明顯，而2006年以后，又開始了大幅度提液生產(chǎn)操作。

圖4 目標井J02井含水率、產(chǎn)液量變化曲線

以現(xiàn)場實際地質(zhì)模型為基礎(chǔ)，采用Eclispe軟件對該區(qū)塊進行模擬，針對J02井的含水率變化擬合曲線見圖5。

很明顯，當不進行任何設(shè)置的時候，J02井在2005～2006年的時候，含水率擬合出現(xiàn)很大偏差，整體呈偏高趨勢。并且在2005年之后，當現(xiàn)場進行提液生產(chǎn)的時候，含水率與實際又出現(xiàn)了整體偏低的現(xiàn)象。這充分說明，在商業(yè)軟件Eclispe中，針對嚴重非均質(zhì)儲層，提液/降液生產(chǎn)時，模型內(nèi)部并不會根據(jù)提液/降液機理進行分析計算，忽略了低滲透層的啟動壓力問題，導(dǎo)致含水率擬合出現(xiàn)很大偏差。

根據(jù)儲層的沉積構(gòu)造和儲層特性，得到該井附近的儲層29、30小層滲透率特別高,其余小層滲透率很低的情況，滲透率級差大于40。根據(jù)先前的機理研究，當進行降液生產(chǎn)時，則在降液生產(chǎn)的時間段，針對除高滲透層每一個射孔層位，采用WPIMULT設(shè)置，乘積因子設(shè)置為0.7，此處的0.7為經(jīng)驗值，可以根據(jù)實際的含水率變化情況和小層產(chǎn)液結(jié)構(gòu)進行適當調(diào)整；針對29、30小層，采用WPIMULT設(shè)置乘積因子為4.0，此處4.0依舊為經(jīng)驗值，可以根據(jù)實際的含水率變化情況和小層產(chǎn)液結(jié)構(gòu)進行適當調(diào)整。

調(diào)整之后的含水率變化曲線見圖6。從圖6中可以看出，修改之后的歷史擬合效果變得很好，這也進一步驗證了上述機理的準確性和Eclispe中的WPIMULT關(guān)鍵字在處理現(xiàn)場提液和降液生產(chǎn)方面的優(yōu)越性。

圖5 目標井J02井初期歷史擬合含水率變化曲線

5 結(jié)論

給出了如何運用商業(yè)軟件Eclispe中的WPIMULT關(guān)鍵字實現(xiàn)對現(xiàn)場提液、降液生產(chǎn)過程中帶來的含水率擬合困難的解決方法，并從機理分析入手，采用現(xiàn)場的實際生產(chǎn)數(shù)據(jù)，進行了該方法的優(yōu)越性闡述。結(jié)果表明，采用該方法可以有效地模擬出非均質(zhì)儲層大幅度提液、降液生產(chǎn)帶來的擬合困難問題，對研究指導(dǎo)具有重要意義。

[1]王公昌, 李林凱, 王平, 等. 啟動壓力在數(shù)值模擬軟件中等效模擬的新方法[J]. 科學(xué)技術(shù)與工程, 2012, 12(1):165～167.

[2]黃延章. 低滲透油層非線性滲流特征[J]. 特種油氣藏, 1997, 4(1):9～14.

[3]李洋, 雷群, 劉先貴,等. 微尺度下的非線性滲流特征[J]. 石油勘探與開發(fā), 2011, 38(3):336～340.

[4]王恩志, 韓小妹, 黃遠智. 低滲巖石非線性滲流機理討論[J]. 巖土力學(xué), 2003,(S2):120～124.

[5]劉建軍, 劉先貴, 胡雅衽. 低滲透巖石非線性滲流規(guī)律研究[J]. 巖石力學(xué)與工程學(xué)報, 2003, 22(04):556～561.

[6]時宇, 楊正明, 黃延章. 低滲透儲層非線性滲流模型研究[J]. 石油學(xué)報, 2009, 30(5):731～734.

[7]李宗翔, 孟憲臣, 紀傳仁. 冒落非均質(zhì)平面流場非達西滲流新迭代算法[J]. 力學(xué)與實踐, 2007, 29(2):28～30.

[8]于春生, 楊曉培, 蘇娜,等. 網(wǎng)絡(luò)模型在非線性滲流模擬中的應(yīng)用[J]. 重慶科技學(xué)院學(xué)報(自然科學(xué)版), 2012, 14(6):88～92.

[9]馮其紅, 王守磊, 白軍偉, 等. 層間非均質(zhì)油藏提液效果數(shù)值模擬[J]. 油氣地質(zhì)與采收率, 2013, 20(3):49～52.

[10]張超. 水驅(qū)油藏提液提高采收率理論分析與實例[J]. 科技信息, 2012,(32):395～396.

[11]王濤. 響應(yīng)曲面法在油藏提液效果影響因素研究中的應(yīng)用[J]. 重慶科技學(xué)院學(xué)報(自然科學(xué)版), 2015, 17(3):67～70.

[12]張超, 鄭川江, 肖武, 等. 特高含水期提液效果影響因素及提高采收率機理——以勝坨二區(qū)沙二段74-81單元為例[J]. 油氣地質(zhì)與采收率, 2013, 20(5):88～91.

[13]楊浩, 楊進, 周長所, 等. 生產(chǎn)壓差對稠油儲層產(chǎn)量的影響[J]. 石油鉆采工藝, 2008, 30(3):63～65.

[14]趙建勇. 勝坨油田產(chǎn)液規(guī)模調(diào)整技術(shù)經(jīng)濟界限研究[D].東營：中國石油大學(xué)(華東), 2012.

[15]朱焱, 謝進莊, 楊為華, 等. 提高油藏數(shù)值模擬歷史擬合精度的方法[J]. 石油勘探與開發(fā), 2008, 35(2):225～229.

[16]陳兆芳, 張建榮, 陳月明, 等. 油藏數(shù)值模擬自動歷史擬合方法研究及應(yīng)用[J]. 石油勘探與開發(fā), 2003, 30(4):82～84.

[17]Burgoyne M W, Little A J H. From high perm oil to tight gas-a practical approach to model hydraulically fractured well performance in coarse grid reservoir simulators[J].SPE156610-MS,2012.

[18]Kuzevanov M A, Buchinskiy S V, Glumov D N,et al. Building of a multiple layer oil and gas condensate field integrated simulation model (Russian)[J].SPE171255-RU,2014.