劉道信
(北京交通大學(xué)中國產(chǎn)業(yè)安全研究中心,北京 100044)
可動凝膠型聚合物驅(qū)油方面,國內(nèi)、外已開展了大量的室內(nèi)實驗機理研究,但在數(shù)值模擬方面研究較少。本文在考慮了驅(qū)油重要因素[1-3]的基礎(chǔ)上,建立了三維三相十組分可動凝膠型聚合物驅(qū)數(shù)學(xué)模型,解法采用有限體積法(FVM),形成了數(shù)值模擬軟件[4]。并針對大慶采油二廠先導(dǎo)試驗區(qū)開展了機理性研究,包括:可動凝膠型聚合物注入體積、段塞組合、殘余阻力系數(shù)及注入不可動凝膠的對比等,從理論上預(yù)測了注入?yún)?shù)的影響因素,設(shè)計了注采合理參數(shù)。
可動凝膠型聚合物驅(qū)油機理數(shù)學(xué)模型包括:交聯(lián)反應(yīng)過程中的組分產(chǎn)生和損耗、化學(xué)組分的吸附和滯留、一二價離子的質(zhì)量交換、凝膠的降解、流動性能,以及化學(xué)組分對水相粘度和滲透率下降系數(shù)等。
聚合物和交聯(lián)劑不斷損耗,逐漸生成凝膠。在成膠時間達到之前,溶液的流動類似于未交聯(lián)的聚合物。在成膠時間達到以后,則隨濃度的變化形成弱凝膠體系。可動凝膠型聚合劑的特點是有凝膠的強度和粘度,但又能隨壓力變化而流動。
(1)
式中:Ri為單位時間內(nèi),由于交聯(lián)反應(yīng)消耗或產(chǎn)生的i組分質(zhì)量濃度,i=p,r,g(聚合物、交聯(lián)劑、凝膠);C為相濃度;Kp,Kr為化學(xué)反應(yīng)速度系數(shù),利用實驗成膠時間得出[2];指數(shù)m,n為反應(yīng)階數(shù)(常數(shù),m≥1,n≥1);右式中首個負號表示質(zhì)量消耗。
生成的可動凝膠能夠影響水相的流動特性[6],表現(xiàn)為流動阻力增加??赏ㄟ^下述方程描述:
Rkg=Rkg(Cg,v)
(2)
式中:Rkg為水相滲透率降低系數(shù);Cg為可動凝膠濃度;v為可動凝膠的流動速度。
這組參數(shù)一般直接用給定不同含鹽量下的Rkg~Cgg表格插值計算,聚合物和交聯(lián)劑降低系數(shù)描述方法與凝膠類似。
由于動態(tài)描述凝膠在巖石表面上的吸附比較困難,為近似遵循Langmuir等溫吸附理論[7]的靜態(tài)吸附,并且考慮聚合物的吸附與鹽度的關(guān)系為可逆的,與濃度的關(guān)系為不可逆的??捎孟率矫枋觯?/p>
(3)
式中:qadg為吸附量;a1,d1為平衡吸附常數(shù),其值由實驗室測定;C為相濃度。
對于不符合Langmuir吸附規(guī)律的情況,可用實驗室給定的不同含鹽量下的聚合物吸附曲線插值計算吸附量[8],聚合物和交聯(lián)劑吸附量描述方法與凝膠類似。
(4)
式中:
先導(dǎo)試驗區(qū)為一個五點井組,油層深度1089.3m,地層溫度45℃,原始地層壓力12.85MPa,飽和壓力10.28MPa,地層原油粘度10.2MPa·s,地層水礦化度7445mg/L,二價離子含量20mg/L,平均有效厚度14.55m,平均滲透率627.9×10-3μm2地質(zhì)儲量34.1×104t。該模擬井組平面上劃分為420 個非結(jié)構(gòu)性網(wǎng)格(X方向20個,Y方向21個),縱向上網(wǎng)格6層,總網(wǎng)格數(shù)為2520個。注入化學(xué)劑溶液參數(shù)見表1。
表1 凝膠型聚合物性質(zhì)
為了選擇合適的可動凝膠型聚合物注入體積,開展了溶液注入體積分別為0.38PV、0.48PV和0.57PV的開發(fā)效果預(yù)測(圖1)。設(shè)定條件:注入凝膠聚合物溶液濃度為6000mg/L,注入凝膠聚合物5個月基本為不可動。之后,可動凝膠聚合物溶液濃度減為1000mg/L,年產(chǎn)液速度11%,后期為凝膠驅(qū)油,含水和累積產(chǎn)油預(yù)測曲線見圖1。首先,可以看出,凝膠型聚合物驅(qū)可較大幅度增加采油量及減小含水率,效果比繼續(xù)水驅(qū)效果好得多,井組可累積增油42.0×103m3左右,含水率最低降至60%。同時,可以看出,注入體積越大越累積產(chǎn)油越多,結(jié)合單井增油量及噸化學(xué)劑增油量分析,建議初注入體積為0.48PV。
圖1 可動凝膠聚合物不同用量開發(fā)效果對比
在可動凝膠型聚合物驅(qū)總用量為480mg/L×PV的基礎(chǔ)上,選擇了三種段塞組合進行開發(fā)效果預(yù)測。一種是可動凝膠型聚合物溶液濃度為1000mg/L,注入段塞為0.48PV;第二種由三個段塞組成,可動凝膠型聚合物溶液濃度為1500mg/L,注入段塞為0.08 PV+可動凝膠型聚合物溶液濃度為1000mg/L,注入段塞為0.3 PV+可動凝膠型聚合物溶液濃度為600mg/L,注入段塞為0.1PV;第三種由兩個段塞組成,可動凝膠型聚合物溶液濃度為1500mg/L,注入段塞為0.1 PV+可動凝膠型聚合物溶液濃度為800mg/L,注入段塞為0.4125PV。預(yù)測結(jié)果見表2。結(jié)合增產(chǎn)及噸化學(xué)劑增油量來看,采用兩個段塞可動凝膠型聚合物溶液濃度為1500mg/L,注入段塞為0.1 PV+可動凝膠型聚合物溶液濃度為800mg/L,注入段塞為0.4125PV時最佳。
表2 可動凝膠型聚合物溶液段塞組合開發(fā)效果對比
為了研究不同殘余阻力系數(shù)對凝膠驅(qū)的影響,依據(jù)現(xiàn)有配方聚合物溶液濃度4000mg/L時的殘余阻力系數(shù)42.2,分別計算了殘余阻力系的0.8、1.2、1.4倍時的開發(fā)效果預(yù)測。預(yù)測條件是不可動凝膠型聚合物4個月后,進行可動凝膠聚合物驅(qū),結(jié)果見表3。由計算結(jié)果可見看出,殘余阻力系數(shù)在一定范圍越大,增加采收率值越高,但增油量及聚合物效率的增長幅度在大于1.2倍后減小。結(jié)合實際油藏的井口注入壓力及采油速度的影響,目前殘余阻力系數(shù)比較合理。
表3 不同殘余阻力系數(shù)開發(fā)效果對比
在可動凝膠型聚合物驅(qū)油前,可注入高濃度的交聯(lián)系劑,使可動凝膠體系的成膠粘度和成膠時間加強,達到油藏深部調(diào)剖的作用。因此,開展了高濃度凝膠注入時間的研究,選擇注入高濃度凝膠劑溶液為2個月、3個月、4個月和5個月;注入2個月累積注入0.024PV,注入3個月累積注入0.028PV,注入4個月累積注入0.032PV,注入5個月累積注入0.038PV。預(yù)測結(jié)果見圖2。由計算結(jié)果可見看出,注入調(diào)剖時間越長越好,但是注入4個月與5個月的累積產(chǎn)油量也接近,變化幅度變小。綜合考慮認為:在可動凝膠型聚合物驅(qū)油前,可考慮注入4個月的高濃度凝膠型聚合物作為油藏深部調(diào)剖。
圖2 不同凝膠劑注入時間開發(fā)效果對比
1)高濃度可動凝膠型聚合物體系的成膠粘度強和成膠時間長,可用于油藏的深度調(diào)剖;低濃度可動凝膠型聚合物體系配制成弱凝膠體系,可用于整體區(qū)塊的驅(qū)油工作,是一種具有很好應(yīng)用前景的提高采收率的手段。
2)通過可動凝膠型聚合物體系的重要化學(xué)驅(qū)油機理研究,建立了數(shù)學(xué)模型,解法采用有限體積法(FVM),編制了三維三相九組分可動凝膠型聚合物驅(qū)數(shù)值模擬器。該模擬器可用于凝膠型聚合物體系的數(shù)值模擬實驗研究。
3)結(jié)合勝利油田試驗區(qū),凝膠型聚合物驅(qū)可較大幅度增加采油量及減小含水率,井組可累積增油42.0×103m3左右,含水率最低降至60%,提高采收率達14.2%。同時,對可動凝膠型聚合物的注入體積、段塞組合、殘余阻力系數(shù)及注入不可動凝膠等參數(shù)進行了優(yōu)化分析,指導(dǎo)了開發(fā)參數(shù)的制定。
[1] Smith J E.The Transition Pressure:A Quick Method for Quantifying Polyacrylamide Gel Strength,SPE #18739[C].Presented at the SPE International Symposium on Oilfield Chemistry in Houston,TX,February 8-10,1989.
[2] Todd B J,Willhite G P,Green D W.A Mathematical Model of In-Situ Gelation of Polyacrylamide by a Redox Process,SPE/DOE #20215[C].Presented at the SPE/DOE Seventh Symposium on Enhanced Oil Recovery,Tulsa,OK,April 22-25,1990.
[3] 袁士義.聚合物地下交聯(lián)調(diào)剖數(shù)學(xué)模型[J].石油學(xué)報,1991,12(1):49-59.
[4] Mack J C,Smith J E.In- Depth Colloidal Dispersion Gels Improve Oil Recovery Efficiency,SPE/DOE #27780[C].
Presented at the SPE/DOE Ninth Symposium on Improved Oil Recovery[J].Tulsa,OK,April 17-20,1994.
[5] Gao,H.W.and M.M.Chang.A Three-Dimensional,Three-Phase Simulator for Permeability Modification Treatments Using Gelled Polymers[R].DOE Report NIPER-388,March 1990.NTIS Order No.DE90000227.
[6] Camilleri,D.,Engelsen,S.,Lake,L.W.,et al.Description of an Improved Compositional Micellar/Polymer Simulator[R].SPE #13967,SPERE,November 1987.
[7] 李幅塏.黑油和組分模型的應(yīng)用[M].北京:石油工業(yè)出版社,1996.
[8] 付天郁,曹鳳,邵振波.聚合物驅(qū)控制程度的計算方法及應(yīng)用[J].大慶石油地質(zhì)與開發(fā),2004,23(3):81-82.