水平井二維穩(wěn)定滲流場分析及應(yīng)用

玄建

中國石化勝利油田分公司勘探開發(fā)研究院，山東 東營 257000

引 言

隨著鉆井技術(shù)的不斷進(jìn)步，水平井在油氣田開發(fā)中得到廣泛應(yīng)用，尤其是在低滲透、頁巖油氣藏開發(fā)中發(fā)揮了較大作用[1-4]。很多學(xué)者研究水平井時(shí)，重點(diǎn)研究水平井產(chǎn)量問題。Joshi通過垂直及水平面內(nèi)的二維滲流問題，推導(dǎo)出均質(zhì)油藏水平井穩(wěn)態(tài)產(chǎn)能方程[5]。也有學(xué)者對水平井二維滲流進(jìn)行了深入研究：李曉軍等通過保角變換研究斷塊油藏水平井與直井組合井網(wǎng)滲流場[6]，張強(qiáng)等按照橢球體模型分析穩(wěn)定滲流場[7]，袁淋等重點(diǎn)研究水平井及其酸化措施后的產(chǎn)能問題[8-10]?？傮w來看，目前，國內(nèi)外的研究成果公式極其復(fù)雜，難以推廣應(yīng)用。在薄層油藏開發(fā)中，為提高開發(fā)效果，應(yīng)用水平井開發(fā)是一種常用方法[11-12]，這種情況可簡化為水平井二維滲流，因此，研究水平井二維穩(wěn)定滲流具有重要意義。本次研究通過儒可夫斯基變換，把水平井二維穩(wěn)定滲流轉(zhuǎn)換為直井徑向二維穩(wěn)定滲流，開展水平井二維穩(wěn)定滲流場中等勢線、流線以及滲流速度研究，并據(jù)此開展水平段長度優(yōu)化、水平段方向優(yōu)化、水平井垂直注水井一注一采井網(wǎng)優(yōu)化，并給出油氣田開發(fā)過程中部署水平井的具體應(yīng)用方法。

1 水平井二維穩(wěn)定滲流場分析

水平井二維穩(wěn)定滲流場是滲流力學(xué)分析中的重要內(nèi)容，研究水平井二維穩(wěn)定滲流場具有重要意義。由于水平井滲流場較復(fù)雜，因此應(yīng)用保角變換分析。保角變換是一種化難為易的方法，將邊界形狀比較復(fù)雜的滲流場轉(zhuǎn)換為較為簡單的滲流場。根據(jù)儒可夫斯基公式變換

式中：w-保角變換后變量，m；

c-橢圓簇半焦距，m；

z-保角變換前變量，m。

z平面內(nèi)的同心圓簇變換為w平面內(nèi)以（–c，0）、（c，0）為焦點(diǎn)的橢圓簇，特別地，當(dāng)z平面內(nèi)的同心圓半徑逼近c(diǎn)時(shí)，w平面內(nèi)的橢圓逼近從（–c，0）到（c，0）的一條割線，由于水平井可視為無限導(dǎo)流管，因此這條割線視為二維滲流場中的水平井，那么，滲流場中的等勢線為共焦點(diǎn)的橢圓簇

式中：x-橫坐標(biāo)，m；

y-縱坐標(biāo)，m；

a-橢圓半長軸，m；

b-橢圓半短軸，m。

相應(yīng)地，z平面內(nèi)通過圓心的射線簇變換為w平面內(nèi)以（–c，0）、（c，0）為焦點(diǎn)的雙曲線簇，即水平井滲流場中的流線簇

式中：a′-雙曲線半實(shí)軸，m；

b′-雙曲線半虛軸，m。

根據(jù)保角變換，求得水平井產(chǎn)量公式

式中：Q-水平井產(chǎn)量，m3/s；

pe-油藏邊界壓力，Pa；

pc-井底壓力，Pa；

μ-流體黏度，mPa·s；

K-滲透率，D；

h-油藏厚度，m；

ae-油藏邊界橢圓半長軸，m；

be-油藏邊界橢圓半短軸，m。

在z平面內(nèi)，任一點(diǎn)滲流速度為

式中：vz-z平面內(nèi)滲流速度，m/s；

θ-流線夾角，rad；

r-半徑，m；

dlz-z平面內(nèi)等勢線上兩點(diǎn)reiθ、rei(θ+dθ)的距離，m。

在w平面內(nèi)，任一點(diǎn)滲流速度為

式中：vw-w平面內(nèi)滲流速度，m/s；

dlw-dlz在w平面內(nèi)的變換，m。

由滲流速度公式可以看出，當(dāng)r等于1時(shí)，代表水平井井壁上的滲流速度，即流體流入水平井時(shí)，在水平井中點(diǎn)（θ=π/2 或θ=3π/2）滲流速度最小，愈靠近水平井兩端（θ=0或θ=π）滲流速度愈大（圖1）。

圖1 水平井井壁不同位置上的滲流速度圖Fig.1 The map of seepage velocity on the borehole wall of horizontal well

2 應(yīng)用算例

2.1 水平段長度優(yōu)化

某常規(guī)薄層圓形油藏面積約3.1 km2（半徑1 km），具有恒定邊界壓力20 MPa，厚度3.0 m，單儲系數(shù)15.0×104t/（km2·m），滲透率0.3 D，地下原油黏度5.0 mPa·s。在油藏中心部署一口水平井，生產(chǎn)壓差2.0 MPa，水平井產(chǎn)量根據(jù)水平井二維穩(wěn)定滲流公式計(jì)算得到。水平井日產(chǎn)油量和日產(chǎn)油量增長率隨水平段長度變化情況見圖2。從圖2可見，當(dāng)水平段生產(chǎn)井段（如：射孔完井時(shí)的射孔井段）小于100 m時(shí)，水平井日產(chǎn)油量增長率較大，當(dāng)水平段生產(chǎn)井段（射孔井段）超過100 m后，水平井日產(chǎn)油量增長率小于0.5%。當(dāng)油藏半徑不同（2，4，8 km），油藏其他參數(shù)均相同時(shí)，水平井日產(chǎn)油量增長率與半徑為1 km的油藏一致，即當(dāng)水平段生產(chǎn)井段（射孔井段）超過100 m后，水平井日產(chǎn)油量增長率小于0.5%，因此該類油藏部署水平井時(shí)，建議水平段生產(chǎn)井段（射孔井段）不宜超過100 m。

圖2 水平井日產(chǎn)油量和日產(chǎn)油量增長率隨水平段長度變化圖Fig.2 The change of daily output and its increasing rate with the length of horizontal section

2.2 水平段方向優(yōu)化

油氣田開發(fā)過程中部署水平井時(shí)通常根據(jù)數(shù)值模擬及油藏工程方法優(yōu)化水平段方向[13-14]，而根據(jù)水平井二維穩(wěn)定滲流場優(yōu)化水平段方向也同樣是一種重要的方法。根據(jù)水平井二維穩(wěn)定滲流速度公式計(jì)算得到水平井二維穩(wěn)定滲流場中的等勢線和流線分布（圖3）。滲流場中任意r處，在水平段的延伸方向（θ=0、θ=π）滲流速度最大，滲流速度可簡化為

在水平段的中垂線方向（θ=π/2、θ=3π/2）滲流速度最小，滲流速度可簡化為

因此，水平井二維穩(wěn)定滲流場中，當(dāng)r相等時(shí)，水平段延伸方向的滲流速度始終大于中垂線上的滲流速度，因此，與一口直井的二維徑向穩(wěn)定滲流不同，在水平井二維穩(wěn)定滲流場中存在主流線，且主流線為水平段延伸方向（θ=0、θ=π）的直線，因此在部署水平井時(shí)需要考慮水平段的方向。

圖3 水平井二維穩(wěn)定滲流場中等勢線及流線分布圖Fig.3 The map of the equipotential lines and the streamlines in the two-dimensional steady seepage field of horizontal well

墾71塊東營組油藏是具有氣頂和較強(qiáng)邊水能量的構(gòu)造油藏，地質(zhì)儲量642.0×104t，含油面積4.8km2。開發(fā)初期，氣竄和邊水錐進(jìn)現(xiàn)象較為普遍，為減緩開發(fā)過程中的氣竄和邊水錐進(jìn)，計(jì)劃在該油藏腰部部署水平井開發(fā)，那么氣竄方向和邊水錐進(jìn)方向必須避開水平井滲流場中的主流線。由于構(gòu)造油藏中的油氣界面、油水界面平行于構(gòu)造線，因此部署水平井時(shí)，水平段方向應(yīng)與構(gòu)造線平行，才能實(shí)現(xiàn)該類油藏水平井開發(fā)的最佳效果。2016年，墾71塊東營組油藏在腰部部署6口水平井，方案實(shí)施后，綜合含水由96.2%降至81.7%，開發(fā)效果較好。

2.3 水平井垂直注水井一注一采井網(wǎng)優(yōu)化

水平井采油、垂直注水井注水是常見的一注一采井網(wǎng)，在開發(fā)方案設(shè)計(jì)中常遇到井位優(yōu)化問題，很多學(xué)者的研究較為復(fù)雜[15-20]，使用不便。筆者在水平井二維穩(wěn)定滲流場分析基礎(chǔ)上，通過平移、旋轉(zhuǎn)及儒可夫斯基變換，得到垂直注水井水平井形成的一注一采二維穩(wěn)定滲流場流線分布圖（圖4）。

圖4 水平井垂直注水井不同方位組合一注一采流線分布圖Fig.4 The streamlines of different patterns with one horizontal well and one vertical well

圖4中，垂直注水井的位置不同，則流線分布不同，但主流線均在以（–c，0）、（c，0）為焦點(diǎn)的雙曲線上，對判斷水平井來水方向和見水位置具有一定指導(dǎo)意義。特別地，當(dāng)垂直注水井在水平井中垂線上時(shí)（θ=π/2 或θ=3π/2），主流線即為中垂線；當(dāng)垂直注水井在水平井延伸方向上時(shí)（θ=0或θ=π），主流線與水平井方向一致，當(dāng)垂直注水井在任一位置時(shí)（θ=θ0），主流線在雙曲線上。根據(jù)推導(dǎo)的水平井滲流場滲流速度公式，結(jié)合水平井垂直注水井不同方位組合形成的一注一采流線分布圖，得出水平井中垂線上是部署垂直注水井的最優(yōu)位置。當(dāng)垂直注水井部署在中垂線上時(shí)，與其它位置的垂直注水井相比，水驅(qū)波及面積相對較大，主流線滲流速度相對較小，水驅(qū)前緣推進(jìn)相對較慢，能夠延長水平井無水采油時(shí)間，提高經(jīng)濟(jì)效益。

3 結(jié) 語

（1）在水平井二維穩(wěn)定滲流場分析基礎(chǔ)上，得到水平井二維穩(wěn)定滲流過程中滲流速度參數(shù)公式。特別地，在水平井井壁上，水平井中點(diǎn)處滲流速度最小，愈靠近水平井兩端滲流速度愈大。

（2）常規(guī)油藏水平井生產(chǎn)井段（射孔井段）不宜超過100 m，實(shí)際應(yīng)用中，應(yīng)當(dāng)根據(jù)油藏的具體參數(shù)優(yōu)化得到最佳的水平段長度；水平段方向應(yīng)于構(gòu)造線平行。

（3）水平井與垂直注水井形成的一注一采二維穩(wěn)定滲流場中的流線分布圖，主流線在共焦點(diǎn)的雙曲線上，對判斷水平井來水方向和見水位置具有一定指導(dǎo)意義。分析得出垂直注水井在水平井的中垂線上時(shí)，水驅(qū)波及面積相對較大，開發(fā)效果較好，是最佳部署位置。

[1]王歡，廖新維，趙曉亮，等.超低滲透油藏分段多簇壓裂水平井產(chǎn)能影響因素與滲流規(guī)律 以鄂爾多斯盆地長8超低滲透油藏為例[J].油氣地質(zhì)與采收率，2014，21（6）：107-110.doi：10.3969/j.issn.1009-9603.2014.06.027 WANG Huan,LIAO Xinwei,ZHAO Xiaoliang,et al.A study on productivity and flow regimes of segmented multi-cluster fractured horizontal well in ultra-low perme-ability reservoir:A case of Chang 8 ultra-low permeability reservoir in Ordos Basin[J].Petroleum Geology and Recovery Efficiency,2014,21(6):107–110.doi:10.3969/j.-issn.1009-9603.2014.06.027

[2]徐慶巖，于靖之，蔣文文，等.特低滲透油藏壓裂水平井與直井聯(lián)合開采井網(wǎng)產(chǎn)能預(yù)測新方法[J].科學(xué)技術(shù)與工程，2014，14（1）：30-34.doi：10.3969/j.issn.1671-1815.2014.18.007 XU Qingyan, YU Jingzhi, JIANG Wenwen, et al. A new method of predicting productivity of combined fractured horizontal well and vertical well for ultra-low permeability reservoir[J]. Science Technology and Engineering, 2014,14(1): 30–34. doi: 10.3969/j.issn.1671-1815.2014.18.007

[3]卞曉冰，張景臣，張士誠.低/特低滲透油藏壓裂水平井布縫模式研究[J].科學(xué)技術(shù)與工程，2015，15（11）：58-62.doi：10.3969/j.issn.1671-1815.2015.11.010 BIAN Xiaobing,ZHANG Jingchen,ZHANG Shicheng.Research of fracture distribution pattern for fractured horizontal wells in low&ultra-low permeability reservoirs[J].Science Technology and Engineering,2015,15(11):58–62.doi:10.3969/j.issn.1671-1815.2015.11.-010

[4]郭小哲，周長沙.頁巖氣藏壓裂水平井滲流數(shù)值模型的建立[J].西南石油大學(xué)學(xué)報(bào)（自然科學(xué)版），2014，36（5）：90-96.doi：10.11885/j.issn.1674-5086.2013.09.-11.07 GUO Xiaozhe,ZHOU Changsha.Seepage numerical model for fractured horizontal well in shale gas reservoir[J].Journal of Southwest Petroleum University(Science&Technology Edition),2014,36(5):90–96.doi:10.-11885/j.issn.1674-5086.2013.09.11.07

[5]鄧森，王怒濤，饒慶鋒，等.考慮摩阻損失的水平井長度優(yōu)化[J].大慶石油地質(zhì)與開發(fā)，2017，36（1）：61-64.doi：10.3969/j.issn.1000-3754.2017.01.010 DENG Sen,WANG Nutao,RAO Qingfeng,et al.Optimization of the horizontal well length considering the frictionloss[J].Petroleum Geology and Oilfield Development in Daqing,2017,36(1):61–64.doi:10.3969/j.issn.1000-3754.2017.01.010

[6]李曉軍，張琪，路智勇.斷塊油藏水平井與直井組合井網(wǎng)滲流解析解及其應(yīng)用[J].中國石油大學(xué)學(xué)報(bào)（自然科學(xué)版），2010，34（3）：85-88.doi：10.3969/j.issn.1673-5005.2010.03.018 LI Xiaojun,ZHANG Qi,LU Zhiyong.Analytic solution of horizontal-vertical composed well pattern seepage in fault-block reservoir and its application[J].Journal of China University of Petroleum(Edition of Natural Science),2010,34(3):85–88.doi:10.3969/j.issn.1673-5005.2010.-03.018

[7]張強(qiáng)，劉生國，張俊松，等.考慮流體壓縮性的應(yīng)力敏感油藏水平井穩(wěn)定滲流模型[J].巖性油氣藏，2012，24（4）：121-124.doi：10.3969/j.issn.1673-8926.-2012.04.024 ZHANG Qiang, LIU Shengguo, ZHANG Junsong, et al.A steady seepage flow model of horizontal wells based on fluid compressibility and stress sensitivity[J]. Lithologic Reservoirs, 2012, 24(4): 121–124. doi: 10.3969/j.issn.-1673-8926.2012.04.024

[8]袁淋，王朝明，李曉平，等.致密砂巖氣藏氣水同產(chǎn)水平井產(chǎn)能公式推導(dǎo)及應(yīng)用[J].巖性油氣藏，2016，28（3）：121-126.doi：10.3969/j.issn.1673-8926.2016.03.016 YUAN Lin,WANG Chaoming,LI Xiaoping,et al.Deduction and application of productivity formula for horizontal well with producing water in tight sandstone gas reservoir[J].Lithologic Reservoirs,2016,28(3):121–126.doi:10.3969/j.issn.1673-8926.2016.03.016

[9]袁淋，李曉平.水平井酸化后產(chǎn)能公式推導(dǎo)及應(yīng)用[J].巖性油氣藏，2014，26（1）：127-130.doi：10.3969/j.-issn.1673-8926.2014.01.022 YUAN Lin,LI Xiaoping.Deduction and application of productivity formula of horizontal wells after acidizing[J].Lithologic Reservoirs,2014,26(1):127–130.doi:10.-3969/j.issn.1673-8926.2014.01.022

[10]袁淋，李曉平，張璐，等.水平井穩(wěn)態(tài)產(chǎn)能公式對比與分析[J].巖性油氣藏，2013，25（6）：127-132.doi：10.-3969/j.issn.1673-8926.2013.06.024 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.doi:10.3969/j.issn.1673-8926.2013.06.024

[11]李英明，任海洋.水平井技術(shù)在洼70塊超稠油薄層油藏的應(yīng)用[J].斷塊油氣田，2012，19（S1）：40-43.doi：10.6056/dkyqt2012z1011 LI Yingming,REN Haiyang.Application of horizontal well in thin reservoir with super-heavy oil in Block Wa70[J].Fault-Block Oil&Gas Field,2012,19(S1):40–43.doi:10.6056/dkyqt2012z1011

[12]聶彬，劉月田.特低豐度薄層油藏水平井開發(fā)效果評價(jià)體系研究[J].油氣井測試，2013，22（5）：14-17.doi：10.3969/j.issn.1004-4388.2013.05.004 NIE Bin,LIU Yuetian.Research on horizontal well development effect assessment in super-low abundance thin reservoir[J].Well Testing,2013,22(5):14–17.doi:10.-3969/j.issn.1004-4388.2013.05.004

[13]祝效華，諶柯宇，毛宇，等.水平井段脈沖攜巖數(shù)值模擬研究[J].西南石油大學(xué)學(xué)報(bào)（自然科學(xué)版），2016，38（6）：164-171.doi：10.11885/j.issn.1674-5086.2015.-01.14.01 ZHU Xiaohua,CHEN Keyu,MAO Yu,et al.Numerical simulation of cuttings pulse carrying in a horizontal well[J].Journal of Southwest Petroleum University(Science&Technology Edition),2016,38(6):164–171.doi:10.11885/j.issn.1674-5086.2015.01.14.01

[14]曹麗娜，李曉平，羅誠，等.裂縫性低滲氣藏水平井不穩(wěn)定產(chǎn)量遞減探討[J].西南石油大學(xué)學(xué)報(bào)（自然科學(xué)版），2017，39（3）：103-110.doi：10.11885/j.issn.1674-5086.2016.03.21.03 CAO Lina,LI Xiaoping,LUO Cheng,et al.Investigation of unstable production decline in fractured lowpermeability gas reservoir horizontal wells[J].Journal of Southwest Petroleum University(Science&Technology Edition),2017,39(3):103–11.doi:10.11885/j.issn.1674-5086.2016.03.21.03

[15]沈瑞，李兆國，段寶江，等.王窯油田西部低滲透油藏正方形反九點(diǎn)井網(wǎng)加密調(diào)整物理模擬[J].油氣地質(zhì)與采收率，2014，21（3）：89-91.doi：10.13673/j.cnki.cn37-1359/te.2014.03.024 SHEN Rui, LI Zhaoguo, DUAN Baojiang, et al. Physical simulation of square inverted nine-spot pattern infilling and adjustment for low permeability reservoirs in western Wangyao Oilfield[J]. Petroleum Geology and Recovery Efficiency, 2014, 21(3): 89–91. doi: 10.13673/j.cnki.-cn37-1359/te.2014.03.024

[16]杜殿發(fā)，鄭洋，趙艷武，等.反九點(diǎn)聯(lián)合井網(wǎng)水平井蒸汽驅(qū)產(chǎn)能公式[J].特種油氣藏，2016，23（4）：90-92.doi：10.3969/j.issn.1006-6535.2016.04.020 DU Dianfa,ZHENG Yang,ZHAO Yanwu,et al.Horizontal well productivity equation of steam flooding with inverted nine-spot combination well pattern[J].Special Oil&Gas Reservoirs,2016,23(4):90–92.doi:10.3969/-j.issn.1006-6535.2016.04.020

[17]趙春森，肖丹鳳，宋文玲，等.水平井與直井交錯(cuò)井網(wǎng)優(yōu)化方法[J].石油勘探與開發(fā)，2005，32（1）：119-122.doi：10.3969/j.issn.1008-4800.2016.07.098 ZHAO Chunsen,XIAO Danfeng,SONG Wenling,et al.Optimization of the staggered well pattern of mixed horizontal and vertical wells[J].Petroleum Exploration and Development,2005,32(1):119–122.doi:10.3969/j.issn.-1008-4800.2016.07.098

[18]蔣廷學(xué)，卞曉冰，王海濤，等.深層頁巖氣水平井體積壓裂技術(shù)[J].天然氣工業(yè)，2017，37（1）：90-96.doi：10.3787/j.issn.1000-0976.2017.01.011 JIANG Tingxue,BIAN Xiaobing,WANG Haitao,et al.Volume fracturing of deep shale gas horizontal wells[J].Natural Gas Industry,2017,37(1):90–96.doi:10.3787/j.-issn.1000-0976.2017.01.011

[19]趙金洲，陳曦宇，李勇明，等.水平井分段多簇壓裂模擬分析及射孔優(yōu)化[J].石油勘探與開發(fā)，2017，44（1）：117-124.doi：10.11698/ped.2017.01.14 ZHAO Jinzhou,CHEN Xiyu,LI Yongming,et al.Numerical simulation of multi-stage fracturing and optimization of perforation in a horizontal well[J].Petroleum Exploration and Development,2017,44(1):117–124.doi:10.-11698/ped.2017.01.14

[20]田鴻照，孫野.直井與水平井組合吞吐轉(zhuǎn)汽驅(qū)操作參數(shù)優(yōu)選[J].巖性油氣藏，2013，25（3）：127-130.doi：10.3969/j.issn.1673-8926.2013.03.023 TIAN Hongzhao,SUN Ye.Operation parameters optimization of steam flooding with the combination of vertical and horizontal wells[J].Lithologic Reservoirs,2013,25(3):127–130.doi:10.3969/j.issn.1673-8926.2013.03.-023