水平井雙管注聚合理配注量的優化方法

孫 哲，孫 巍，曾 晟，萬慶生

(1.東北石油大學提高油氣采收率教育部重點實驗室，黑龍江大慶 163318；2. 中國石油集團工程設計有限責任公司北京分公司，北京 100000；3. 中國石油大慶油田有限責任公司儲運銷售分公司，黑龍江大慶 163000)

?

水平井雙管注聚合理配注量的優化方法

孫 哲1，孫 巍1，曾 晟2，萬慶生3

(1.東北石油大學提高油氣采收率教育部重點實驗室，黑龍江大慶 163318；2. 中國石油集團工程設計有限責任公司北京分公司，北京 100000；3. 中國石油大慶油田有限責任公司儲運銷售分公司，黑龍江大慶 163000)

為研究水平井配注量問題，以聚合物波及效率最大為目標函數，以注入井、生產井的井底流壓為約束條件，建立了水平井注聚-直井采油混合井網水平井雙管注聚合理配注量優化的數學模型，給出了遼河油田水驅高246塊高2-蓮H16水平井的合理配注量。結果表明，根據上述方法計算出的水平井的合理配注量與油田實際配注量基本吻合，驗證了此方法的正確性。

滲流力學；水平井與直井；水平井雙管注聚；合理配注量；優化方法

近年來，人們在水平井滲流理論分析方面取得了許多研究成果。美國的JOSHI[1]利用電場流理論，假定水平井泄油體是以水平井兩端點為焦點的橢圓體，將三維滲流問題簡化為垂直及水平漸內的二維間題，利用勢能理論詳細推導了均質油藏水平井產能公式；郎兆新等[2]用一種形式簡單的茹可夫斯基函數作變換，獲得了壓力分布、產量計算、見水時間及面積掃油系數公式；曲德斌等[3]應用復位勢理論及迭加原理研究典型的開發井網——一般的水平井五點法面積井網，其中包括產量、見水時間、波及系數的計算以及流函數的導出等；胡軍暉等[4]綜述了水平井產能的4種計算方法，并給出了當儲層存在滲透率各向異性、井偏心或地層損害時修正的產能計算公式；尹洪軍等[5]建立了非均質氣藏水平井三維滲流產能預測的數學模型，采用有限元方法對其進行求解；于金萍[6]應用復位勢理論及疊加原理對水平井注水-直井采油井組、直井注水-水平井和直井采油井組進行了研究，給出了這兩個典型井組水井合理配注量的確定方法。

目前，研究水平井合理配注量的優化方法多是基于水驅的，對于水平井注聚合理配注量優化的問題考慮的較少。水平井可以增大油藏泄油面積，提高注劑的注入速度和波及效率，水平井結合聚合物驅油技術可以提高原油采收率，改善油田開發效果。因此，對于水平井注聚合理配注量優化的研究具有十分重要的實際意義。

1 數學模型

1.1 油層壓力場分布

假設無限大地層一個混合井網中1口水平井雙管注聚，周圍有n口直井采油，供給邊界處的壓力為pe，注聚水平井水平段長度為2l，地層平均滲透率為K，平均有效厚度為h。

1)基本關系式

在Zo平面建立以注聚水平井外管的水平段軸線為xo軸，其水平段中點為坐標原點o的直角坐標系(xo,yo)。在Zi平面建立以注聚水平井內管的水平段軸線為xi軸，其水平段中點為坐標原點o的直角坐標系(xi,yi)。第i口生產井在Zo平面和Zi平面的坐標轉換關系為

xii=xoi-l ，

(1)

yii=yoi。

(2)

應用儒可夫斯基保角變換將真實平面上的復雜流動問題變換為像平面上的單向流動問題。第i口生產井在真實平面和像平面的坐標轉換關系為

(3)

(4)

式中：(xii,yii)為第i口生產井在Zi平面下的坐標；(xoi,yoi)為第i口生產井在Zo平面下的坐標；l為注聚水平井水平段長度的1/2，m；(ξi,ηi)為第i口生產井在像平面下的坐標；(xi,yi)為第i口生產井在真實平面下的坐標。

2)油層壓力場分布函數

根據鏡像理論及勢的疊加原理，推導出水平井與直井混合井網油層壓力場分布函數為

(5)

由式(5)可得第i口生產井井底流壓pwfi為

(6)

其中

(7)

(8)

(9)

(10)

(11)

式中：μs為聚合物溶液在地層中滲流時的剪切黏度，Pa·s；Qoi為第i口生產井受注聚水平井外管注聚驅替的產量，m3·s-1；(ξoi,ηoi)為第i口生產井經保角變換后在外管像平面下的坐標；K為地層滲透率，m2；h為油層的有效厚度，m；Qii為第i口生產井受注聚水平井內管注聚驅替的產量，m3·s-1；(ξii,ηii)為第i口生產井經保角變換后在內管像平面下的坐標；C為常數，可通過供給邊界處的勢確定，m2·s-1；ρowi為第i口生產井在外管像平面上的半徑，m；ρiwi為第i口生產井在內管像平面上的半徑，m；k為稠度系數，Pa·sn；n為流性指數，無因次；υ為滲流速度，m/s；c′為與毛細管迂曲度有關的系數，其變化范圍是25/12～2.5；φ為地層孔隙度，無因次；rwi為第i口生產井在W平面下的半徑，m；Qho為注聚水平井外管注入量，m3·s-1；Qhi為注聚水平井內管注入量，m3·s-1。

1.2 目標函數

目標函數是聚合物波及效率最大，即

(12)

其中

(13)

(14)

(15)

(16)

式中：toi為注入聚合物由注聚水平井外管的井壁沿主流線到達第i口生產井井底的時間，s；tii為注入聚合物由注聚水平井內管的井壁沿主流線到達第i口生產井井底的時間，s；A為井網面積，m2；vyo為任意一點從注聚水平井的外管到生產井的主流線上沿yo方向的滲流速度，m·s-1；vyi為任意一點從注聚水平井的內管到生產井的主流線上沿yi方向的滲流速度，m·s-1。

1.3 約束條件

約束條件為

pR≤piwf≤piwfmax，

(17)

其中

piwfmax=pf(1-X),

(18)

pf=HmY，

(19)

piwfmin=pp+0.01γm(Hm-Hp) ，

(20)

式中：pR為目前地層壓力，MPa；piwfmax為最大合理注聚井井底流壓，MPa；pf為地層破裂壓力，MPa；X為破裂概率，無因次；Hm為油層中部深度，m；Y為破裂壓力梯度，MPa·m-1；piwfmin為最小合理生產井井底流壓，MPa；pp為合理泵口壓力，MPa；γm為生產井井筒內混合液相對密度，無因次；Hp為泵掛深度，m。

2 優化方法

1)輸入水平井雙管注聚-直井采油井組的基礎數據：N，A，h，K，φ，n，k，Hm，pR，Y，2l，piwfmax，Zo，(xoi,yoi)，rwi，γm，Hp，pwf；

2)由式(7)計算聚合物溶液在地層中滲流時的剪切黏度μs；

3)由式(1)—式(4)計算直井生產井在外管、內管像平面下對應的直角坐標(ξoi,ηoi)、(ξii,ηii)；

4)由式(18)計算注聚水平井的最大合理井底流壓piwfmax，由式(20)計算直井生產井的最小合理井底流壓piwfmin；

5)在保證注聚水平井外、內管注入聚合物的推進速度vξ1i和vξ2i相同的條件下，由式(6)計算各直井生產井的產量Qoi和Qii；

6)由(10)計算注聚水平井外管注入量Qho，由式(11)計算注聚水平井內管注入量Qhi；

7)由式(7)計算波及效率Ea；

8)判斷Eamax

9)piwf=piwf+1，判斷式(17)，若成立，則轉至步驟(5)，反之，則轉至步驟(10)；

10)由式(12)可知，波及效率Eamax=max{Ea}，與此對應的注聚水平井外管注入量Qho、內管注入量Qhi為最優外管注入量Qoo、內管注入量Qio。

3 計算實例

根據上述水平井與直井混合井網水平井雙管注聚合理配注量的優化方法，計算遼河油田水驅高246塊高2-蓮H16水平井雙管注水時的合理配注量。

如圖1所示，高246塊高2-蓮H16水平注水井周圍有10口與其連通的直井生產井。含油面積A為1.31km2，平均油層有效厚度h為38.1m，油層平均滲透率k為0.588μm2，孔隙度φ為0.218，油層中部深度Hm為1 560m，高2-蓮H16水平注水井外管水平段長度2lo為120.98m，內管水平段長度2li為222.12m。

圖1 高246塊高2-蓮H16井位圖Fig.1 Well location map of G2-H16 in G246 Block

應用上述優化方法計算的遼河油田水驅高246塊高2-蓮H16水平井注水量與波及效率的關系如圖2所示。

圖2 高246塊高2-蓮H16井注水量與波及效率關系Fig.2 Relationship between injection volume and sweep efficiency of G2-H16 in G246 Block

圖2表明，當外管注水量為27m3/d，內管注水量為46m3/d時，波及效率較高。遼河油田水驅高246塊高2-蓮H16水平井的外管實際日配注量為28m3/d，內管實際日配注量為45m3/d。因此，根據上述方法計算出的水平井的合理配注量與油田實際配注量基本吻合，可以認為本文給出的水平井注聚-直井采油混合井網水平井雙管注聚合理配注量的優化方法是正確的。

4 結 語

建立了水平井注聚-直井采油混合井網水平井雙管注聚合理配注量優化的數學模型，給出了這種混合井網水平井外、內管合理配注量的優化方法。對遼河油田水驅高246塊高2-蓮H16水平井的應用分析表明，根據上述方法計算出的水平井的合理配注量與油田實際配注量基本吻合，驗證了本文數學模型及優化方法的正確性。研究結果為水平井與直井混合井網的滲流理論提供了重要支持，為今后研究更符合油藏中滲流實際情況的水平井數學物理模型提供了參考。

[1] JOSHI S D. A review of horizontal well and drain hole technology [C]// SPE16868.[S.l.]:[s.n.], 1988: 226-230.

[2] 郎兆新, 張麗華, 程林松, 等. 水平井與直井聯合開采問題——五點法面積井網[J]. 石油大學學報(自然科學版), 1993, 17(6)：50-55. LANG Zhaoxin, ZHANG Lihua, CHENG Linsong, et al. Combined production of vertical and horizontal wells—A 5-spot well pattern[J]. Journal of the University of Petroleum, China, 1993, 17(6): 50-55.

[3] 曲德斌, 葛家理, 王德民. 水平井與直井聯合面積布井的開發理論研究(一)——一般的五點面積井網[J]. 石油勘探與開發, 1995, 22(1)：35-38. QU Debin, GE Jiali, WANG Demin. Combination production of vertical and horizontal wells—5-spot well pattern [J]. Petroleum Exploration and Development, 1995, 22(1): 35-38.

[4] 胡軍暉, 郡捷年. 水平井產能的影響因素[J]. 石油勘探與開發, 1994, 21(5)：44-52. HU Junhui, YAN Jienian. The factors affecting the productivity of horizontal well [J]. Petroleum Exploration and Development, 1994, 21(5): 44-52.

[5] 尹洪軍, 呂杭, 蔡明, 等. 非均質氣藏水平井三維滲流的產能預測方法[J]. 數學的實踐與認識, 2011, 41(16)：59-64. YIN Hongjun, LYU Hang, CAI Ming, et al. Productivity prediction method of 3D seepage horizontal wells in heterogeneity gas reservoir [J]. Mathematics in Practice and Theory, 2011, 41(16): 59-64.

[6] 于金萍. 水平井與直井聯合布井不規則井網水井配注量的確定方法研究[D]. 大慶：東北石油大學, 2010. YU Jingping. Study on Determination Approaches of Water Injection Rate of Skewed Pattern with a Combination of Horizontal and Vertical Wells [D]. Daqing: Northeast Petroleum University, 2010.

[7] 鐘會影, 蔡明, 楊釗, 等. 考慮基巖向井筒供氣的氣藏壓裂水平井產能預測模型[J]. 數學的實踐與認識, 2014, 44(18)：132-137. ZHONG Huiying, CAI Ming, YANG Zhao, et al. Productivity prediction model for fractured horizontal wells considering the gas rate of matrix flowing into well bore [J]. Mathematics in Practice and Theory, 2014, 44(18): 132-137.

[8] TSUNEKI I, TATSUNERI D, YOJI N. Stability of oil-water emulsion with mobile surface charge[J]. Colloids and Surfaces A: Physicochemical and Engineering Aspects, 2006, 29(1): 128-141.

[9] JONG Y K, MYUNG G S, JONG D. Zeta potential of nanobubbles generated by ultrasonication in aqueous alkyl polyglycoside solutions[J]. Journal of Colloid and Interface Science, 2000, 23(3): 285-291.

[10]王德民, 程杰成, 吳軍政, 等. 聚合物驅油技術在大慶油田的應用[J]. 石油學報, 2005, 26(1)：74-78. WANG Demin, CHENG Jiecheng, WU Junzheng, et al. Application of polymer flooding technology in Daqing Oilfield[J]. Acta Petrolei Sinica, 2005, 26(1):74-78.

[11]ABOLRAZL E, ELHEM G, TORAJ M. Separation of water-in-oil emulsion using microfiltration[J]. Desalination, 2005, 185(3): 371-382.

[12]孫煥泉. 勝利油田三次采油技術的實踐與認識[J]. 石油勘探與開發, 2006, 33(3)：262-266. SUN Huanquan. Practice and understanding on tertiary recovery in Shengli Oilfield[J]. Petroleum Exploration and Development, 2006, 33(3): 262-266.

[13]袁淋, 李曉平, 張璐, 等. 水平井穩態產能公式對比與分析[J]. 巖性油氣藏, 2013, 25(6)：127-132. YUAN Lin, LI Xiaoping, ZHANG Lu, et al. Analysis and correlation of steady-state productivity formulas for horizontal wells[J]. Lithologic Reservoirs, 2013, 25(6): 127-132.

[14]王大為, 李曉平. 水平井產能分析理論研究進展[J]. 巖性油氣藏, 2011, 23(2)：118-123. WANG Dawei, LI Xiaoping. Advances in deliverability analysis of horizontal well[J]. Lithologic Reservoirs, 2011, 23(2): 118-123.

[15]于雷, 車飛, 許磊明, 等. 水平井產能與人工裂縫參數相關性模擬實驗[J]. 斷塊油氣田, 2015, 22(2)：221-223. YU Lei, CHE Fei, XU Leiming, et al. Correlation simulation experiment of productivity and artificial fracture parameters for horizontal well[J]. Fault-Block Oil & Gas Fie1d, 2015, 22(2): 221-223.

Optimization method of rational polymer injection allocation rate for horizontal well dual pipe polymer injection

SUN Zhe1, SUN Wei1, ZENG Sheng2, WAN Qingsheng3

(1. Key Laboratory of Enhanced Oil Recovery of Ministry Education, Northeast Petroleum University, Daqing, Heilongjiang 163318, China; 2. China Petroleum Engineering Company Limited Beijing Company, Beijing 100000, China; 3. Storage and Transportation Sales Branch, Daqing Oilfield Corp PetroChina, Daqing, Heilongjiang 163000, China)

In order to study the problems of horizontal well polymer injection allocation rate, considering the maximum sweep efficiency as the object function, taking the bottom hole flowing pressure of the injection well and the production wells as the constraint conditions, the mathematical model of the optimization of the horizontal well dual pipe rational polymer injection allocation rate is established under the combined well pattern of horizontal well for polymer injection and vertical wells for oil extraction. The rational injection allocation rate of the G2-H16 horizontal well in waterflood Gao 246 Block of Liaohe Oilfield is given. The results indicate that the calculated rational injection allocation rate of the horizontal well is nearly identical to actual data in the oilfield, which proves the correctness of the method.

seepage mechanics; horizontal and vertical wells; horizontal well dual pipe polymer injection; rational polymer injection allocation rate; optimization method

1008-1534(2016)04-0319-05

2016-03-17;

2016-04-29；責任編輯：馮 民

國家自然科學基金(51574086);東北石油大學研究生創新科研項目(YJSCX2015-010NEPU)

孫 哲(1988—)，女，山東高密人，博士研究生，主要從事提高采收率技術方面的研究。

E-mail:sunzhe1988dqyt@126.com

TE357.46

A

10.7535/hbgykj.2016yx04010

孫 哲，孫 巍，曾 晟, 等.水平井雙管注聚合理配注量的優化方法[J].河北工業科技,2016,33(4):319-323. SUN Zhe, SUN Wei, ZENG Sheng, et al.Optimization method of rational polymer injection allocation rate for horizontal well dual pipe polymer injection[J].Hebei Journal of Industrial Science and Technology,2016,33(4):319-323.