大慶油田二類油層聚合物驅(qū)產(chǎn)油量模型應(yīng)用

張雪玲

(中國(guó)石油大慶油田有限責(zé)任公司，黑龍江 大慶 163712)

?

大慶油田二類油層聚合物驅(qū)產(chǎn)油量模型應(yīng)用

張雪玲

(中國(guó)石油大慶油田有限責(zé)任公司，黑龍江 大慶 163712)

針對(duì)大慶油田不同區(qū)塊的二類油層產(chǎn)油量差異大，影響因素復(fù)雜，產(chǎn)量預(yù)測(cè)難等問題，根據(jù)聚合物驅(qū)油墻推進(jìn)理論和產(chǎn)油量動(dòng)態(tài)變化特征的對(duì)應(yīng)關(guān)系，建立分階段產(chǎn)油量預(yù)測(cè)模型。在模型求解過程中，通過引入支持向量機(jī)方法，選取已開發(fā)區(qū)塊25組模型參數(shù)和影響因素進(jìn)行學(xué)習(xí)訓(xùn)練，5組數(shù)據(jù)樣本進(jìn)行預(yù)測(cè)檢驗(yàn)，最終確定了模型參數(shù)與影響因素之間的非線性擬合關(guān)系式，并將整個(gè)算法過程編制成了簡(jiǎn)便操作的應(yīng)用程序。現(xiàn)場(chǎng)應(yīng)用表明，新的模型預(yù)測(cè)精度在95%以上，便于跟蹤和調(diào)試，完全可以滿足二類油層聚合物驅(qū)開發(fā)方案和油田中長(zhǎng)期規(guī)劃方案的編制要求。

二類油層；聚合物驅(qū)；開發(fā)指標(biāo)；支持向量機(jī)；大慶油田

0 引 言

聚合物驅(qū)是大慶油田三次采油的重要技術(shù)之一，隨著區(qū)塊一類主力油層聚合物驅(qū)的逐步結(jié)束，目前二類油層已成為聚合物驅(qū)的主要開采對(duì)象。與一類主力油層開發(fā)情況相比，二類油層的砂體發(fā)育規(guī)模、連通程度和驅(qū)替方式上都發(fā)生了較大變化[1-4]。

現(xiàn)有的聚合物驅(qū)產(chǎn)油量預(yù)測(cè)方法，如模式圖法、累計(jì)液、累計(jì)水法和驅(qū)替特征曲線法等[5-9]，主要適用于一類主力油層的發(fā)育條件和開發(fā)方式，已不能滿足二類油層的產(chǎn)量預(yù)測(cè)精度要求。為此，根據(jù)產(chǎn)油量動(dòng)態(tài)變化特征與聚合物驅(qū)地下驅(qū)替過程的對(duì)應(yīng)關(guān)系，建立了分階段的產(chǎn)量預(yù)測(cè)模型，探索支持向量機(jī)方法解決模型參數(shù)與主控影響因素之間的非線性問題，現(xiàn)場(chǎng)應(yīng)用驗(yàn)證了預(yù)測(cè)模型及算法的準(zhǔn)確性。

1 產(chǎn)油量模型的建立

1.1 產(chǎn)油量曲線特征

統(tǒng)計(jì)大慶油田二類油層進(jìn)入聚合物驅(qū)中、后期開采的區(qū)塊有24個(gè)，利用歸類統(tǒng)計(jì)和單井跟蹤模擬得到二類油層產(chǎn)油量動(dòng)態(tài)曲線具有以下4個(gè)特征。

(1) 二類油層產(chǎn)油量同樣符合生長(zhǎng)曲線的規(guī)律。產(chǎn)油量曲線可分為4個(gè)階段，分別為低產(chǎn)量低速上升階段、高產(chǎn)量高速上升階段、高產(chǎn)量高速遞減階段和低產(chǎn)量低速遞減階段。由于二類油層的砂體規(guī)模小、油層厚度薄、產(chǎn)量接替能力差，導(dǎo)致穩(wěn)產(chǎn)階段時(shí)間短，產(chǎn)量達(dá)到峰值后迅速進(jìn)入遞減階段(圖1)。

(2) 產(chǎn)油量曲線在高產(chǎn)量高速上升階段和高產(chǎn)量高速遞減階段具有對(duì)稱性。根據(jù)單井跟蹤模擬，聚合物溶液在地下的推進(jìn)過程中逐漸形成油墻[10-11]。隨著地下油墻逐步地推進(jìn)至油井井底，直至全部被采出，一般不實(shí)施增產(chǎn)措施，因此，產(chǎn)油量曲線保持了很好的對(duì)稱性。

(3) 產(chǎn)油量曲線在高產(chǎn)量高速遞減階段和低產(chǎn)量低速遞減階段具有明確的分界點(diǎn)。統(tǒng)計(jì)進(jìn)入低速遞減的14個(gè)區(qū)塊的二類油層，雖然各油區(qū)到達(dá)分界點(diǎn)的時(shí)間存在差異，但是到達(dá)分界點(diǎn)處的開采時(shí)間一般是到達(dá)產(chǎn)油量高峰值時(shí)開采時(shí)間的1.33倍左右(表1)。

圖1 二類油層聚合物驅(qū)產(chǎn)油量動(dòng)態(tài)曲線

油區(qū)區(qū)塊數(shù)量/個(gè)高峰值開采時(shí)間/月分界點(diǎn)開采時(shí)間/月時(shí)間比值采油一礦515～2020～241.20～1.33采油二礦317～1822～241.29～1.33采油三礦217～2223～301.35～1.36采油四礦427～3035～421.29～1.40

(4) 產(chǎn)油量曲線在低速遞減階段符合遞減指數(shù)為0.5的雙曲遞減規(guī)律。對(duì)于二類油層，由于各區(qū)塊開發(fā)后期實(shí)施增產(chǎn)措施的規(guī)模、數(shù)量不同，初始遞減率存在一定的差別。

1.2 產(chǎn)油量預(yù)測(cè)模型

根據(jù)聚合物地下油墻推進(jìn)理論和產(chǎn)油量動(dòng)態(tài)變化特征的對(duì)應(yīng)關(guān)系，對(duì)聚合物驅(qū)見效到高產(chǎn)量高速遞減階段和低產(chǎn)量低速遞減階段，采用不同的函數(shù)進(jìn)行描述。

由于產(chǎn)油量曲線在高產(chǎn)量高速上升階段和高產(chǎn)量高速遞減階段具有對(duì)稱性，因此，對(duì)低產(chǎn)量低速上升階段、高產(chǎn)量高速上升階段和高產(chǎn)量高速遞減階段的產(chǎn)量曲線，采用三角正弦函數(shù)描述：

(1)

對(duì)低速遞減階段的產(chǎn)量曲線采用雙曲遞減函數(shù)描述：

(2)

其中，

ti=tQmax+1.33(tQmax-tQo)

(3)

(4)

式中：q為聚合物驅(qū)月產(chǎn)油量，104t/月；qo為空白水驅(qū)末月產(chǎn)油量，104t/月；t為聚合物驅(qū)開采時(shí)間，月；ΔQmax為聚合物驅(qū)最大增油量，104t；tQo為聚合物驅(qū)開始見效時(shí)的開采時(shí)間，月；tQmax為達(dá)到最大產(chǎn)油量時(shí)的開采時(shí)間，月；qi為開采時(shí)間為ti時(shí)月產(chǎn)油量，104t/月；Di為初始遞減率，月-1；ti為由高速遞減進(jìn)入低速遞減時(shí)的開采時(shí)間，月。

模型中最大增油量ΔQmax定義為最大月產(chǎn)油量與空白水驅(qū)末月產(chǎn)油量的差值。初始遞減率Di由增油階段的采出程度確定，不同的增油階段采出程度對(duì)應(yīng)不同的初始遞減率(圖2)。采用試湊法，先對(duì)實(shí)際區(qū)塊初始遞減率賦一個(gè)值，計(jì)算產(chǎn)油量變化，當(dāng)月產(chǎn)油量遞減到空白水驅(qū)末月產(chǎn)油量時(shí)停止，得到增油階段采出程度，然后與給定的增油階段采出程度對(duì)比。依次循環(huán)，當(dāng)2個(gè)階段采出程度相當(dāng)時(shí)確定出對(duì)應(yīng)的初始遞減率。

2 模型參數(shù)計(jì)算方法

產(chǎn)油量預(yù)測(cè)模型中存在4個(gè)重要的參數(shù)：聚合物驅(qū)見效時(shí)開采時(shí)間、產(chǎn)油量最大時(shí)開采時(shí)間、最大增油量、增油階段采出程度。由于每個(gè)區(qū)塊的油層發(fā)育狀況、剩余油分布以及注入?yún)?shù)不同，這4個(gè)模型參數(shù)也存在很大的不確定性。根據(jù)影響因素相關(guān)性大和易于獲取的原則，確定了油層有效厚度、滲透率、滲透率變異系數(shù)、聚合物驅(qū)控制程度、河道砂控制程度、非河道砂控制程度、連通厚度和注聚合物前初期含水等8項(xiàng)指標(biāo)作為模型參數(shù)的主要影響因素。

圖2 不同初始遞減率對(duì)應(yīng)的產(chǎn)油量曲線

為了求解模型參數(shù)影響因素的函數(shù)關(guān)系式，引入支持向量機(jī)方法。支持向量機(jī)具有嚴(yán)格的統(tǒng)計(jì)學(xué)理論和數(shù)學(xué)基礎(chǔ)，主要思想為應(yīng)用核函數(shù)技術(shù)，將因變量與多個(gè)自變量樣本，通過非線性函數(shù)映射到高維特征空間中，在高維空間中構(gòu)造線性判別函數(shù)，巧妙地解決算法復(fù)雜度與多影響因素的非線性問題。該方法以訓(xùn)練誤差作為優(yōu)化問題的約束條件，以置信范圍值最小化作為優(yōu)化目標(biāo)，將算法轉(zhuǎn)化為一個(gè)求凸二次優(yōu)化問題的全局最優(yōu)解，從總體樣本中挑選出少數(shù)具有代表性的樣本，即所謂支持向量，構(gòu)建擬合函數(shù)。

擬合函數(shù)為：

(5)

式中：j為輸入樣本的個(gè)數(shù)，K(x，xj)為核函數(shù)，常用的核函數(shù)有多項(xiàng)式核函數(shù)、RBF核函數(shù)和sigmoid核函數(shù)。系數(shù)b0，β1，…，βn變成下述二次規(guī)劃問題的解：

(6)

(7)

(8)

其中，式(6)～(8)體現(xiàn)了支持向量機(jī)預(yù)測(cè)誤差最小化和推廣能力最強(qiáng)的準(zhǔn)則。將求得的系數(shù)和確定的核函數(shù)類型帶入式(5)可得到擬合函數(shù)。

此次以8項(xiàng)影響因素作為自變量、每個(gè)模型參數(shù)作為因變量。根據(jù)支持向量機(jī)方法原理，在Libsvm-2.85支持向量機(jī)軟件的基礎(chǔ)上，改進(jìn)了數(shù)據(jù)準(zhǔn)備方式，編制了方便操作的模型參數(shù)預(yù)測(cè)程序。選取已開發(fā)區(qū)塊的25組模型參數(shù)和影響因素作為學(xué)習(xí)樣本，5組數(shù)據(jù)作為檢驗(yàn)樣本，利用支持向量機(jī)軟件計(jì)算得到模型參數(shù)的擬合值和預(yù)測(cè)值。通過做預(yù)測(cè)點(diǎn)與實(shí)際點(diǎn)的關(guān)系圖，可以看出2類點(diǎn)基本上在45°直線附近，說明由學(xué)習(xí)樣本得到的模型參數(shù)預(yù)測(cè)精度較高(圖3、4)。

圖3 最大增油量擬合值與預(yù)測(cè)值

圖4 階段采出程度擬合值與預(yù)測(cè)值

3 產(chǎn)油量預(yù)測(cè)模型的應(yīng)用

以L油田NZKX區(qū)為例說明產(chǎn)油量預(yù)測(cè)模型的應(yīng)用。NZKX區(qū)開采層段為薩Ⅲ1-7，含油面積為7.72 km2，石油地質(zhì)儲(chǔ)量為1 034×104t，2013年2月開始注聚合物，注聚合物時(shí)初期含水為95.3%，目前累計(jì)產(chǎn)油量為62.16×104t，數(shù)值模擬計(jì)算增油階段采出程度為9.58%。

收集整理已開發(fā)區(qū)塊的4個(gè)模型參數(shù)和8項(xiàng)影響因素指標(biāo)，形成學(xué)習(xí)樣本，通過支持向量機(jī)計(jì)算得到模型參數(shù)的擬合函數(shù)。將NZKX區(qū)的油層厚度、滲透率、控制程度、初期含水等8項(xiàng)靜、動(dòng)態(tài)數(shù)據(jù)輸入到支持向量機(jī)軟件，計(jì)算得到2類核函數(shù)的預(yù)測(cè)結(jié)果(表2)。選擇精度較高的徑向級(jí)核函數(shù)預(yù)測(cè)結(jié)果，帶入產(chǎn)油量預(yù)測(cè)模型，得到NZKX區(qū)的產(chǎn)油量預(yù)測(cè)曲線(圖5)。在NZKX區(qū)產(chǎn)油量預(yù)測(cè)模型中，綜合遞減率取0.05，遞減階段分界點(diǎn)在第30個(gè)月，產(chǎn)油量總體預(yù)測(cè)精度達(dá)到了97.94%，預(yù)測(cè)精度完全能夠滿足聚合物驅(qū)開發(fā)方案和中長(zhǎng)期規(guī)劃方案的編制。

表2 L油田NZKX區(qū)模型參數(shù)預(yù)測(cè)結(jié)果

圖5 L油田NZKX區(qū)產(chǎn)油量動(dòng)態(tài)曲線與預(yù)測(cè)結(jié)果

4 結(jié) 論

(1)根據(jù)二類油層地質(zhì)特征和產(chǎn)油量動(dòng)態(tài)曲線變化規(guī)律，建立了分階段的產(chǎn)油量預(yù)測(cè)模型，模型中各項(xiàng)參數(shù)物理意義明確，模型簡(jiǎn)捷、易于操作，預(yù)測(cè)精度達(dá)到95%以上，完全可以滿足規(guī)劃的要求。

(2)通過引入支持向量機(jī)方法，應(yīng)用大慶油田實(shí)際數(shù)據(jù)進(jìn)行了學(xué)習(xí)和預(yù)測(cè)，解決了產(chǎn)油量模型參數(shù)與多項(xiàng)影響因素之間的非線性問題，并編制了擬合函數(shù)算法軟件。結(jié)果表明，支持向量機(jī)方法是一項(xiàng)構(gòu)建開發(fā)指標(biāo)與多項(xiàng)影響因素函數(shù)關(guān)系的新方法。

[1] 李潔，武力軍，邵振波.大慶油田二類油層聚合物驅(qū)油技術(shù)要點(diǎn)[J].石油天然氣學(xué)報(bào)(江漢石油學(xué)院學(xué)報(bào))，2007，27(2)：394-396.

[2] 付京，宋考平，王志華，等.高分子質(zhì)量聚合物溶液與二類油層匹配性研究[J].特種油氣藏，2015，22(2)：129-132.

[3] 李宜強(qiáng)，賈敬鋒.薩南二類油層深度調(diào)剖配方篩選[J].特種油氣藏，2010，17(4)：88-90.

[4] 張曉芹.大慶油田二類油層聚合物驅(qū)注入?yún)?shù)的優(yōu)選[J].大慶石油學(xué)院學(xué)報(bào)，2005，29(4)：40-43.

[5] 孫強(qiáng)，鄧兵，馬麗梅.廣義翁氏與瑞利模型在聚合物驅(qū)產(chǎn)量預(yù)測(cè)中的應(yīng)用[J].大慶石油地質(zhì)與開發(fā)，2003，22(5)：58-59.

[6] 張平，王磊，唐放，等.喇嘛甸油田北北塊薩Ⅲ4-10油層高濃度聚合物驅(qū)數(shù)值模擬研究[J].特種油氣藏，2010，17(6)：69-82.

[7] 陳福明，盧金鳳，陳鵬.聚合物驅(qū)開采指標(biāo)測(cè)算方法研究[J].大慶石油地質(zhì)與開發(fā)，1999，18(2)：33-37.

[8] 孔祥亭，唐莉，周學(xué)民.聚合物驅(qū)開發(fā)規(guī)劃指標(biāo)預(yù)測(cè)方法研究[J].大慶石油地質(zhì)與開發(fā)，2001，20(5)：46-49.

[9] 宋考平，陳銳，鄧軍，等.聚合物驅(qū)產(chǎn)量和含水率變化規(guī)律的預(yù)測(cè)[J].大慶石油學(xué)院學(xué)報(bào)，2002，26(1)：97-99.

[10] 計(jì)秉玉，戰(zhàn)劍飛，蘇致新.油井見效時(shí)間和見水時(shí)間計(jì)算公式[J].大慶石油地質(zhì)與開發(fā)，2000，19(5)：24-26.

[11] 王錦梅，陳國(guó)，歷燁，等.聚合物驅(qū)油過程中形成油墻的動(dòng)力學(xué)機(jī)理研究[J].大慶石油地質(zhì)與開發(fā)，2007，26(6)：64-66.

編輯 朱雅楠

20151102；改回日期：20151223

國(guó)家科技重大專項(xiàng)“大型油氣田及煤層氣開發(fā)”子項(xiàng)目“大慶長(zhǎng)垣特高含水油田提高采收率示范工程”(2011ZX05052)

張雪玲(1977-)，女，工程師，2004年畢業(yè)于中國(guó)石油大學(xué)(華東)石油工程專業(yè)，2008年畢業(yè)于中國(guó)石油大學(xué)(北京)油氣田開發(fā)工程專業(yè)，獲碩士學(xué)位，現(xiàn)從事大慶油田開發(fā)規(guī)劃及戰(zhàn)略研究工作。

10.3969/j.issn.1006-6535.2016.02.031

TE33

A

1006-6535(2016)02-0128-04