礫巖油藏采收率理論公式推導與應用

楊新平，陳光喜，姜瑞忠，王玉多，張利鋒

(1.中國石油新疆油田分公司，新疆 克拉瑪依 834000；2.中國石油大學(華東)，山東 青島 266580；3.克拉瑪依地質工程有限責任公司，新疆 克拉瑪依 834000)

?

礫巖油藏采收率理論公式推導與應用

楊新平1，陳光喜1，姜瑞忠2，王玉多3，張利鋒1

(1.中國石油新疆油田分公司，新疆 克拉瑪依 834000；2.中國石油大學(華東)，山東 青島 266580；3.克拉瑪依地質工程有限責任公司，新疆 克拉瑪依 834000)

為從理論角度分析礫巖油藏各影響因素與采收率之間的函數關系，根據油田實際生產情況，結合非活塞式水驅油基礎滲流理論，綜合考慮滲透率、孔隙度、儲層非均質性、油水黏度比、相滲參數、井網參數等因素對波及系數的影響，推導出變截面B-L方程，在此基礎上得到了礫巖油藏采收率的理論公式。計算克拉瑪依油田13個礫巖油藏的采收率可知，中、低滲礫巖油藏理論公式結果與動態分析法的平均相對誤差僅為3.5%，表明方法準確、可靠。該方法可以從理論角度計算不同儲層條件不同井距下的理論采收率，從而為礫巖油藏合理井網、井距的確定提供決策依據。

礫巖油藏；采收率；滲流理論；變截面B-L方程；克拉瑪依油田

0 引 言

礫巖油藏是克拉瑪依油田最早發現和投入開發的油藏類型，目前年產油占全油田的23.4%，所占比重較高。可采程度為71.1%，綜合含水為73.5%，已進入高采出程度、高含水的“雙高”階段，因此，對礫巖油藏采收率開展深入研究具有非常重要的意義。

1 變截面B-L方程理論推導

在兩相滲流區中，油水界面為非等面積驅動[10]，任取一微小單元體，在dt時間內流入流出的水相體積差值為：

(1)

式中：Vb為水相體積變化量，m3；QL為每個剖面上的流體流量，m3/s；fw為含水率，%；A(x)為橫截面積，m2；t為流入流出時間，s；dx為單元體厚度，m。

由于dt時間內存在流入流出的水相體積差值，因此，微小單元體的含水飽和度將發生變化。其水相體積變化值為：

(2)

式中：Sw為地層含水飽和度，%；φ為地層孔隙度，%。

聯立式(1)、(2)可得：

(3)

NASA-TLX采用6個維度分析任務負荷，每個維度以1到21的刻度表示任務負荷。6個維度下的分數被重新量化為0(低任務負荷)到100(高任務負荷)的值。NASA-TLX 6個維度的平均值和標準偏差見表2。

注水量基本穩定時，QL為常數，由物質平衡可得到等飽和度面移動規律。

(4)

式(4)即為考慮變截面的B-L方程，實現了橫截面從常量向變量轉變的理論模型。

2 變橫截面礫巖油藏采收率理論模型推導

根據變截面B-L方程，對式(4)分離變量并積分可得：

(5)

(6)

(7)

式中：H為歸整后高度，m。

非均質性越強，平面波及系數越小，以滲透率變異系數表征平面非均質性，則H可表示為[13]：

(8)

式中：d為當前井網下的井排距，m；K為空氣滲透率，10-3μm2；VK為平面滲透率變異系數；Cs為砂體形態系數(砂體的周長與砂體的面積之比)，m-1；Kro(Swi)為原始含水飽和度下油相相對滲透率；Krw(Sor)為殘余油飽和度下水相相對滲透率；μo、μw分別為地下原油、地層水黏度，mPa·s；Sc為井網密度，口/km2；a1、a2、a3、a4為井網系數；b為與油層非均質性有關的系數；nw、no為井網的注水、采油井數，口；n為空間維數；m為相對滲透率系數。

層間非均質性對縱向波及程度影響嚴重，主要表現為高滲透層動用程度高，低滲透層動用程度低甚至未動用[14]。根據縱向非均質概念模型分析縱向波及系數與層間非均質性的合理函數關系，得到縱向波及系數與層間滲透率突進系數近似呈冪律關系，則見水時間為：

(9)

式中：he為縱向波及系數；TK為層間滲透率突進系數。

在定液量生產的條件下，累計產油量為：

(10)

式中：Np為累計產油量，104t；A、B、C、D為原油物性相關常數；Swc為平衡飽和度，%；Sor為殘余油飽和度，%；Δp為生產壓差，MPa；Bo為原油體積系數；rw為井半徑，m。

采出程度為：

(11)

式中：N為石油地質儲量，104t；R為采出程度，%。

式(11)即為考慮了礫巖油藏滲透率、平面滲透率變異系數、層間突進系數、孔隙度、相滲參數、油水黏度比、井網參數等多個影響因素的采出程度公式。

3 應用實例

克拉瑪依油田共有230個礫巖油藏，統計百21井區百口泉組、八區上克拉瑪依組等多個礫巖油藏的滲透率、變異系數、相滲參數、油水黏度比等參數(表1)，利用上述推導的理論公式與動態法分別計算不同類型礫巖油藏的采收率，得到礫巖中、低滲油藏平均絕對誤差為1.4%，平均相對誤差為3.5%，誤差較小，礫巖特低滲油藏平均絕對誤差為8.7%，平均相對誤差為29.2%(表2)，表明基于變截面B-L方程推導的礫巖油藏采收率理論公式對克拉瑪依油田中、低滲礫巖油藏具有很好的適用性，但對特低滲礫巖油藏適用性較差。

表1 理論采收率計算參數

表2 理論計算與動態法預測采收率值對比

圖1、2分別為百21百口泉組、三疊系上克拉瑪依組油藏井距為100～500 m和100～400 m的理論采收率曲線，極限含水率時的采收率即為對應井距下的理論采收率。初期分別采用500、350 m井距開發，在采出程度分別達到21%、17%時進行了井網加密，井距分別調整為320、175 m，生產情況得到改善，并逐漸分別與300、200 m井距的理論曲線相吻合，對應理論最終采收率分別為36%、30%，與動態分析法計算的最終采收率基本一致，證明理論公式符合油藏實際生產情況。同時，也可以根據采收率反算合理井網井距，從而為合理井距的確定提供理論基礎和參考方法。

圖1 百21百口泉組油藏不同井距下采出程度與含水率關系

圖2 三疊系上克拉瑪依組油藏不同井距下采出程度與含水率關系

4 結 論

(1) 以基礎滲流理論為切入點，通過推導出變截面B-L方程，得到了考慮滲透率、平面滲透率變異系數、層間突進系數、孔隙度、相滲參數、油水黏度比、井網井距等參數的礫巖油藏采出程度理論公式。

(2) 根據推導出的理論公式，通過計算克拉瑪依油田不同類型礫巖油藏的理論采收率，表明基于變截面B-L方程推導的礫巖油藏采收率理論公式對克拉瑪依油田中低滲礫巖油藏具有很好的準確性和適用性，對特低滲礫巖油藏適用性較差。

(3) 通過對克拉瑪依油田多個實際油藏的采收率進行理論公式計算和動態分析法計算，結合井網井距圖版，驗證了模型的適用性，為礫巖油藏合理井網井距評價及預測提供了參考方法。

[1] 溫靜.遼河油田特高含水期油藏二元復合驅井網井距優化研究[J].特種油氣藏，2012，19(2)：73-75.

[2] 巴拉基洛夫 Ю A.關于原油采收率的計算[J].斷塊油氣田，1996，3(2)：48-49.

[3] 陳元千，鄒存友.對謝爾卡切夫(Щелкачев)公式的推導及拓展[J]. 斷塊油氣田，2010，17(6)：729-732.

[4] 郭蘭磊.聚驅后油藏化學驅提高采收率技術及先導試驗[J].大慶石油地質與開發，2014，33(1)：122-126.

[5] 高淑玲.聚驅后殘留聚合物對提高采收率的影響[J].大慶石油地質與開發，2014，33(6)：112-116.

[6] 張曉亮，楊仁鋒，李娜.考慮啟動壓力梯度的彈性采收率計算方法及影響因素[J].油氣地質與采收率，2015，22(2)：72-77，87.

[7] 侯吉瑞，張麗，李海波，等.碳酸鹽巖縫洞型油藏氮氣驅提高采收率的影響因素[J].油氣地質與采收率，2015，22(5)：64-68.

[8] JERRY E Ongerth. The concentration of cryptosporidium and giardia in water-the role and importance of recovery efficiency[J]. Water Research，2013， 47(7)：2478-2488.

[9] 姚約東，雍潔，朱黎明，等. 砂礫巖油藏采收率的影響因素與預測[J]. 石油天然氣學報，2010，32(4)：108-113.

[10] 王增林，王敬，劉慧卿，等.非均質油藏開發規律研究[J]. 油氣地質與采收率，2011，18(5)：63-66.

[11] 齊與峰.砂巖油田注水開發合理井網研究中的幾個理論問題[J].石油學報，1990，11(4)：51-60.

[12] 范江，張子香，等.非均質油層波及系數計算模型[J].石油學報，1993，14(1)：92-98.

[13] 唐洪明，文鑫，張旭陽，等. 層間非均質礫巖油藏水驅油模擬實驗[J]. 西南石油大學學報(自然科學版)，2014，36(5)：129-135.

[14] 程曉倩，劉華勛，熊偉，等. 新疆低滲透砂礫巖油藏自發滲吸實驗研究[J]. 科學技術與工程，2013，13(26)：7793-7797.

編輯 劉 巍

10.3969/j.issn.1006-6535.2016.05.023

20151228；改回日期：20160720

中國石油天然氣股份有限公司“新疆大慶”科技重大專項“砂礫巖油藏提高采收率技術研究與應用”(2012E-34-07)

楊新平(1971-)，男，高級工程師，1995年畢業于石油大學(華東)油藏工程專業，現從事油氣田開發方面的工作。

TE349

A

1006-6535(2016)05-0096-04