吳海君，杜慶軍，宮汝祥

（1.中國(guó)石油大學(xué)（華東），山東 青島266580；2.中海油田服務(wù)股份有限公司油田生產(chǎn)研究院，天津300450）

0 引言

水平井多元熱流體吞吐技術(shù)已在渤海稠油開發(fā)實(shí)踐中取得了較好的增產(chǎn)效果［1-4］。由于注氣量較大且油藏非均質(zhì)性嚴(yán)重，特別是大孔道的存在，使得井間氣竄現(xiàn)象嚴(yán)重，影響了注入井和鄰井的動(dòng)態(tài)［5-9］，降低了熱能利用率和多元熱流體的開發(fā)效果，干擾了油井的正常生產(chǎn)［10-13］。必須采取關(guān)井以及減小注氣強(qiáng)度等措施，以便減小氣竄的影響。在后續(xù)生產(chǎn)過程中，也需對(duì)不同竄流程度的氣竄現(xiàn)象采取不同的調(diào)整措施［14-16］。因此，需要建立不同氣竄程度的表征方法，從而對(duì)井間竄流程度進(jìn)行分類研究，為后續(xù)調(diào)整方式的制定提供指導(dǎo)。

本文利用數(shù)值模擬方法，建立了存在大孔道的水平井多元熱流體井間竄流基礎(chǔ)模型，提出用井間竄流系數(shù)、 竄流距離和形狀因子3 種參數(shù)來定量表征氣竄程度；在此基礎(chǔ)上，討論了大孔道滲透率和注入量對(duì)3種表征方法的影響。

1 基礎(chǔ)模型的建立

南堡×油田x 井生產(chǎn)層位為Nm05 小層，垂直深度943 m，生產(chǎn)初期地層壓力8 MPa，水平段長(zhǎng)215 m。在多輪次吞吐生產(chǎn)期間，該井與鄰近的y 井之間存在明顯的氣竄特征。在注氣期間，x 井的油壓和套壓迅速下降，y 井的油壓由0 上升至2.5 MPa；在下泵生產(chǎn)期，x井地層流壓下降較快，由初始的9.4 MPa 下降到3.0 MPa，產(chǎn)能下降迅速。根據(jù)y 井和x 井的近井地層性質(zhì)，分析認(rèn)為，兩井間存在大孔道。以x 井實(shí)際油藏特征為依據(jù)，建立了水平井多元熱流體吞吐竄流模型，其井位及滲透率分布見圖1。

模型基本參數(shù)為：油藏溫度56 ℃，地下原油黏度665 mPa·s，油層厚度6 m，初始含油飽和度0.68，平面滲透率4 000×10-3μm2；在I，J 和K 方向上劃分的網(wǎng)格數(shù)為51×25×8，網(wǎng)格尺寸為25 m×15 m×0.75 m，水平段長(zhǎng)215 m。

采用吞吐井交替開井開發(fā)方式，4 口邊井以120 m3/d 的排液速度同時(shí)生產(chǎn)2 個(gè)月后，自中心井開始注入多元熱流體。其中熱水注入溫度為240 ℃，注入速度為165 m3/d，地面標(biāo)準(zhǔn)狀況下，氣體注入速度為57 600 m3/d，注入量為3 300 m3。

圖1 模型井位及滲透率分布

2 竄流程度的表征方法

4 口邊井在生產(chǎn)過程中壓力逐漸降低，產(chǎn)油量逐漸減小，當(dāng)中心井注入多元熱流體后，由于氣體的擴(kuò)散和井間的壓力差，注入氣會(huì)向井周圍擴(kuò)散。為了便于研究井間氣體的竄流程度，以其中一口井（3#井）的日產(chǎn)氣和累計(jì)產(chǎn)氣曲線為研究對(duì)象，提出竄流系數(shù)的概念，定量表征井間氣體的竄流程度。利用氣體驅(qū)替前緣場(chǎng)圖，提出了竄流前緣距離和形狀因子2 個(gè)參數(shù)，對(duì)竄流程度進(jìn)行表征。

2.1 井間竄流系數(shù)

圖2為3#井日產(chǎn)氣量和累計(jì)產(chǎn)氣量隨生產(chǎn)時(shí)間的變化?？梢钥闯觯涸谥行木膺^程中，3#井日產(chǎn)氣量逐漸增大；當(dāng)中心井注入結(jié)束時(shí)，3#井的日產(chǎn)氣量達(dá)到最大；在之后的生產(chǎn)過程中，日產(chǎn)氣量逐漸減小。

在注入量一定的條件下，竄流越嚴(yán)重，則鄰井日產(chǎn)氣量的峰值越大，生產(chǎn)結(jié)束時(shí)的累計(jì)產(chǎn)氣量也越大，因此，引入竄流系數(shù)的概念來表征井間竄流程度的大小。其定義為，在一個(gè)生產(chǎn)周期內(nèi)，鄰井最大日產(chǎn)氣量和鄰井累計(jì)產(chǎn)氣量的乘積與累計(jì)注氣量平方的比值。其表達(dá)式為

式中：n 為井間竄流系數(shù)；Qmax為鄰井最大日產(chǎn)氣量，m3；Qc為鄰井累計(jì)產(chǎn)氣量，104m3；Qi為中心井注氣總量，104m3。

數(shù)值模擬計(jì)算得出，在基礎(chǔ)方案下，3#井日產(chǎn)氣峰值為7 364.4 m3，累計(jì)產(chǎn)氣量為25.34×104m3，累計(jì)注氣量為111.8×104m3。用式（1）計(jì)算得到基礎(chǔ)方案的竄流系數(shù)為14.43。

圖2 3#井日產(chǎn)氣量和累計(jì)產(chǎn)氣量隨時(shí)間的變化

2.2 竄流前緣距離和形狀因子

圖3為高滲帶滲透率為平均滲透率的8 倍、 注入結(jié)束時(shí)的含氣飽和度分布。竄流前緣距離定義為注氣結(jié)束時(shí)注入井與竄流前緣間的最大距離。圖中L 即為竄流前緣距離。當(dāng)高滲帶滲透率為平均滲透率的8 倍時(shí)，L 值為158.94 m。

圖3 竄流前緣距離示意

為了研究氣體擴(kuò)散前緣形狀的不規(guī)則程度，引入形狀因子表征圖形規(guī)則程度，形狀因子越大，圖形越規(guī)則。形狀因子定義為，注氣結(jié)束時(shí)氣體擴(kuò)散前緣所包含的面積與周長(zhǎng)平方的比值。即：

式中：m 為形狀因子；s 為面積，m2；c 為周長(zhǎng)，m。

當(dāng)高滲帶滲透率倍數(shù)為8 倍時(shí)，周長(zhǎng)為695.13 m，面積為17 115.15 m2，根據(jù)式（2）計(jì)算得到形狀因子為0.035 42。

3 大孔道滲透率對(duì)表征參數(shù)的影響

分別討論了井間大孔道滲透率依次為平均滲透率的4，8，12，16 和20 倍時(shí)3#井產(chǎn)氣量的變化，并計(jì)算得到相應(yīng)的竄流系數(shù)、竄流前緣距離和形狀因子。

3.1 對(duì)竄流系數(shù)的影響

不同大孔道滲透率倍數(shù)下，3#井日產(chǎn)氣量的變化如圖4所示。不同滲透率時(shí)竄流系數(shù)見表1。由圖可知，隨著高滲條帶滲透率的增大，3#井日產(chǎn)氣量峰值增加，且出現(xiàn)的時(shí)間更早；竄流系數(shù)也隨高滲條帶滲透率的增加而迅速增加。其主要原因在于，隨著高滲帶滲透率的增大，更多的氣體會(huì)以更快的速度沿高滲條帶突破，造成鄰井產(chǎn)氣量增加，使注入氣體利用率降低。

圖4 不同高滲帶滲透率倍數(shù)下3#井日產(chǎn)氣量

表1 不同高滲帶滲透率倍數(shù)下的竄流系數(shù)

3.2 對(duì)竄流前緣距離的影響

圖5為高滲帶滲透率倍數(shù)分別為4，8，12，16 和20 倍時(shí)注入20 d 后含氣飽和度分布，表2為不同滲透率倍數(shù)下注入結(jié)束時(shí)的注入井與竄流前緣的距離L值。由圖可以看出，大孔道滲透率越大，氣體沿高滲條帶向生產(chǎn)井竄流程度越嚴(yán)重，在相同時(shí)刻，竄流前緣與注入井間距離越大。

圖5 不同滲透率倍數(shù)下注入20 d 后含氣飽和度分布

表2 不同滲透率倍數(shù)下的竄流前緣距離

為了解竄流前緣距離隨注入時(shí)間的變化，分別模擬了滲透率倍數(shù)為8 時(shí)注入4，8，12，16，20 d 時(shí)含氣飽和度分布。模擬結(jié)果表明，當(dāng)注入時(shí)間較短時(shí)，無明顯氣竄，但隨著注入時(shí)間的增加，注入量增加，氣體擴(kuò)散范圍增大，竄流前緣距離也隨之增大。

利用數(shù)值模擬技術(shù)，綜合模擬研究了不同大孔道滲透率倍數(shù)在不同時(shí)間竄流前緣距離的變化，結(jié)果見圖6。由圖可以看出：注入時(shí)間越長(zhǎng)，相同滲透率倍數(shù)條件下的竄流前緣距離越大；滲透率倍數(shù)越大，竄流前緣距離隨時(shí)間的增加幅度越大。

圖6 注入時(shí)間對(duì)竄流前緣距離的影響

3.3 對(duì)竄流形狀因子的影響

利用數(shù)值模擬技術(shù)，還研究了大孔道不同滲透率倍數(shù)下在注入結(jié)束時(shí)竄流形狀因子的變化，結(jié)果見表3。由表可知，隨著滲透率的增大，形狀因子逐漸降低，竄流前緣不規(guī)則程度增加。這主要是由于，高滲帶滲透率倍數(shù)越大，氣體沿高滲帶竄進(jìn)程度越明顯，氣體擴(kuò)散范圍越不規(guī)則。

表3 不同滲透率倍數(shù)下的形狀因子統(tǒng)計(jì)

4 結(jié)論

1）提出了井間竄流系數(shù)、竄流距離和形狀因子3個(gè)參數(shù)來定量表征竄流程度的大小。

2）隨著大孔道滲透率的增大，竄流系數(shù)和竄流距離均增大，形狀因子減小。

3）隨著注入量的增大，竄流距離逐漸增大；注入量相同時(shí)，大孔道滲透率越大，形狀因子越小，竄流程度越大。

［1］黃穎輝，劉東，羅義科.海上多元熱流體吞吐先導(dǎo)試驗(yàn)井生產(chǎn)規(guī)律研究［J］.特種油氣藏，2013，20（2）：84-86.

［2］林濤，孫永濤，劉海濤，等.CO2，N2與蒸汽混合增效作用［J］.斷塊油氣田，2013，20（2）：246-247.

［3］唐曉旭，馬躍，孫永濤.海上稠油多元熱流體吞吐工藝研究及現(xiàn)場(chǎng)試驗(yàn)［J］.中國(guó)海上油氣，2011，23（3）：185-188.

［4］劉小鴻，張風(fēng)義，黃凱，等.南堡35-2 海上稠油油田熱采初探［J］.油氣藏評(píng)價(jià)與開發(fā)，2011，1（1）：61-64.

［5］林濤，孫永濤，孫玉豹，等.多元熱流體返出氣增產(chǎn)技術(shù)研究［J］.斷塊油氣田，2013，20（1）：126-128.

［6］何應(yīng)付，高慧梅，周錫生.改善特低滲透油藏注二氧化碳驅(qū)油效果的方法［J］.斷塊油氣田，2011，18（4）：512-515.

［7］楊兵，李敬松，祁成祥，等.海上稠油油藏多元熱流體吞吐開采技術(shù)優(yōu)化研究［J］.石油地質(zhì)與工程，2012，26（1）：54-56.

［8］楊光璐，唐震，王中元.邊底水稠油油藏蒸汽驅(qū)注采系統(tǒng)優(yōu)化研究［J］.特種油氣藏，2006，13（2）：51-55.

［9］邵紅云，陳雷，李林輝，等.勝利油田注聚區(qū)塊深調(diào)防竄技術(shù)研究與應(yīng)用［J］.斷塊油氣田，2005，12（1）：81-83.

［10］王建波，高云叢，王科戰(zhàn).腰英臺(tái)特低滲透油藏CO2驅(qū)油井見氣規(guī)律研究［J］.斷塊油氣田，2013，20（1）：118-122.

［11］陳民鋒，郎兆新，莫小國(guó).超稠油油藏蒸汽吞吐參數(shù)優(yōu)選及合理開發(fā)界限的確定［J］.石油大學(xué)學(xué)報(bào)：自然科學(xué)版，2002，26（1）：39-42.

［12］鄭家朋，東曉虎，劉慧卿，等.稠油油藏注蒸汽開發(fā)汽竄特征研究［J］.特種油氣藏，2012，19（6）：72-75.

［13］曹嫣鑌，劉冬青，張仲平，等.勝利油田超稠油蒸汽驅(qū)汽竄控制技術(shù)［J］.石油勘探與開發(fā)，2012，39（6）：739-743.

［14］張紅玲，劉慧卿，王晗，等.蒸汽吞吐汽竄調(diào)剖參數(shù)優(yōu)化設(shè)計(jì)研究［J］.石油學(xué)報(bào)，2007，28（2）：105-108.

［15］牛淑芳.注蒸汽氣化氮碳泡沫封堵汽竄技術(shù)及應(yīng)用［J］.油氣地質(zhì)與采收率，2005，12（1）：78-80.

［16］耿朝暉.北過蒸汽驅(qū)汽竄控制技術(shù)初探［J］.中外能源，2012，17（7）：49-52.