楊矞琦,郭繼香,王 翔,王詩棋,劉小雨,吳 楊,朱博圓,楊小輝,王育杉
(1.中國石油大學(北京)重質(zhì)油國家重點實驗室克拉瑪依分室,新疆克拉瑪依 834000;2.中國石油大學(北京)非常規(guī)油氣科學技術研究院,北京 102249;3.中國石油化工股份有限公司西北油田分公司,新疆巴音郭楞蒙古自治州 841600)
重油冷采主要采用降黏劑,使用水溶性好的表面活性劑配制成一定濃度的水溶液,添加到重油中,形成O/W型乳狀液?;瘜W降黏劑室內(nèi)降黏效果與實際井筒降黏效果差距較大,根本原因是對加入表面活性劑后復雜流體井筒舉升流動規(guī)律認識不清[1-3]。在垂直上升的液-液流動中,觀察到的流動型態(tài)一般可分為6種主要類型,即分散流、泡沫流、段塞流、彈狀流、中心環(huán)狀流和擾動流[4-8]。一般認為,影響油水兩相流型的因素很多,包括流體性質(zhì)(如密度、黏度、表面或界面張力、相互溶解度等)、管道性質(zhì)(如幾何形狀、尺寸和管道方向)以及工況(如溫度、壓力、流向、流量、流量比等)。目前對于油水兩相流動型態(tài)的影響規(guī)律研究普遍采用油水界面張力在20~34 mN/m,形成的油-水多相流的分類、流型圖、流型轉(zhuǎn)變機制等建立在物性簡單、密度較低的油相基礎上[9-13]。筆者采用自主研發(fā)的垂直管流裝置,實驗研究加入表面活性劑對重油水兩相垂直管流流動型態(tài)的影響機制,油水界面張力小于1 mN/m,建立加入表面活性劑條件下重油水流動型態(tài)圖版,分析油水乳狀液在垂直管流過程中的形成規(guī)律,推導流型轉(zhuǎn)變界限方程。
材料:油相為透明環(huán)烷油(30 ℃時黏度為581 mPa·s,密度為1.889 g/cm3),一種主要由環(huán)烷烴組成的石油組分,包括環(huán)戊烷、環(huán)己烷及其同系物,以及質(zhì)量分數(shù)為4.894%的氯,濟南丹龍化工有限公司;原油樣品,塔河油田;水相為自來水(密度為1.03 g/cm3);亞甲基藍,北京博恩試劑有限公司;曲拉通-100表面活性劑,北京博恩試劑有限公司;SDG-2水溶性復合降黏劑,實驗室自制[14]。
儀器:E5637-C型界面張力計,法國桑切斯技術;垂直管流裝置(實驗室自制),海安石油科研儀器有限公司[15]。
采用懸滴法,將裝有油樣的樣品池放置在油浴內(nèi),油浴鍋內(nèi)硅油的界面低于樣品池中油樣的界面,以便能夠通過攝像技術觀察到水滴的界面形變,通過界面張力計測定油/水界面張力。油/水界面張力測定的溫度為30 ℃,每次測量時間為3 h,直到界面張力穩(wěn)定,測量至少重復3次。
在長2.00 m、內(nèi)徑0.01 m的不銹鋼管內(nèi)進行油水兩相垂直向上流動測量,管長、直徑比(L/D)為100,垂直管流裝置流程圖如圖1所示。Bergman等[16]提出完全發(fā)展的流動是指由于流體顆粒和管壁之間的剪切應力引起的黏性效應,當流體穿過一定長度(定義為進水特征長度)的直管時,流體產(chǎn)生完全發(fā)展的速度分布。因為本文中在整個實驗條件范圍內(nèi)油和水都處于層流狀態(tài)。在層流的情況下,進水特征長度[17]計算公式為
Lh/D=0.05Re.
(1)
式中,Lh為進水特征長度,m;Re為雷諾數(shù);D為管道的直徑,m。
計算得到的對應于水和油雷諾數(shù)最大值的Lh/D分別為0.143和75,均低于實驗的L/D值200,說明在本文中操作條件下,可以達到完全發(fā)展流動。
圖1 垂直管流裝置流程Fig.1 Schematic diagram of experimental setup
實驗步驟:將油通過注入泵A泵入循環(huán)泵中,將循環(huán)泵中的油通過伺服電機驅(qū)動活塞移動以一定的流速注入垂直管道,同時利用注入泵B將容器B中的表面活性劑水溶液通過注入點泵入垂直管道中的上升油流中,注水管段直徑0.003 m,位于垂直上升管的底部,注入點的幾何形狀見圖2。
當循環(huán)泵中的活塞從一側(cè)移動到另一側(cè)時,一組實驗過程結束。入口持水率Yw由0至90%可計算如下:
Yw=Qw/(Qw+Qo).
(2)
式中,Qw和Qo分別為水和油的流量,m3/s。
圖2 將水注入主要測試管段的注入點Fig.2 Injection point for water into oil in the main test loop
油、水流量測定及實驗條件:在固定的入口持水率、溫度和壓力條件下,當流動達到穩(wěn)定狀態(tài)時,通過高速攝像和電導探針技術對流動型態(tài)進行觀察和記錄,每個流動條件持續(xù)10 min。
根據(jù)實際井筒乳化降黏開采過程中,向原油中加入乳化降黏劑后,油水界面張力一般小于1 mN/m,這樣有利于水包油乳狀液的形成,達到降黏減阻的目的。為考察油水界面張力小于1 mN/m時的重油油水兩相流動規(guī)律,配置質(zhì)量分數(shù)為0.5%的曲拉通-100表面活性劑水溶液,使油水界面張力為0.752 mN/m,固定流速0.03 m/s,控制持水率為10%~80%,常溫常壓下,觀察加入表面活性劑對不同持水率油水兩相在垂直管道中流動型態(tài)的影響,并與不加表面活性劑的油水兩相流動特征進行對比,結果見圖3、4。
圖3 加表面活性劑流速為0.03 m/s不同持水率下油水兩相流動型態(tài)Fig.3 Flow pattern of oil-water two-phase flow at a flow rate of 0.03 m/s with surfactant
圖4 不加表面活性劑速度為0.03 m/s時油水兩相流動型態(tài)Fig.4 Oil-water two-phase flow patterns at a velocity of 0.03 m/s without surfactant
從圖3、4看出,環(huán)烷油在可視釜中表現(xiàn)為淡黃色,亞甲基藍染色的水呈深藍色,水與環(huán)烷油主流從注入點(底部)流入,可視釜(頂部)流出。當持水率從10%逐漸增加到80%的過程中,主要觀察到8種流動型態(tài),分別為油包水分散流(DW/O)、油包水泡狀流(BW/O)、油包水彈狀流(SW/O)、油包水蠕狀流(CEW/O)、環(huán)狀流(AF)、擾動流(Churn)、水包油泡狀流(BO/W)和水包油分散流(DO/W)。其中加入表面活性劑后的擾動流(Churn Flow)(圖3)與常規(guī)油水兩相流動的擾動流[18](圖4)相差較大,表現(xiàn)為水相穿過油相呈帶狀分布并在延伸方向上波動的剪切流型。出現(xiàn)該流型的原因是在油水界面張力非常小的情況下,局部持水率急劇增多使水相紊流程度加劇,環(huán)狀流無法繼續(xù)維持而發(fā)生徑向波動,導致流動型態(tài)呈卷曲條帶狀。加入表面活性劑溶液后,油水界面張力急劇降低,水相易與油相形成水膜,表面活性劑水溶液與油相形成泡沫狀,各個流動型態(tài)也不再是傳統(tǒng)油水兩相流動型態(tài),表現(xiàn)為兩種或兩種以上流動型態(tài)并存,如持水率10%時表現(xiàn)為油包水泡狀流和分散流并存,但以泡狀流為主;持水率30%時,油包水彈狀流和泡狀流并存,彈狀流為主;持水率40%時,油包水蠕狀流、彈狀流、泡狀流并存,蠕狀流為主。持水率增大到60%時,水相受油相徑向阻力增大,由于加入表面活性劑水溶液與油相界面張力降低,水相無法繼續(xù)保持環(huán)狀流,在油相中形成擾動流。繼續(xù)增大入口持水率至70%時,相反轉(zhuǎn)發(fā)生,油為連續(xù)相的流動轉(zhuǎn)變?yōu)橐运疄檫B續(xù)相的流動,出現(xiàn)水包油泡狀流和水包油分散流。當入口持水率大于30%運行一段時間后,整個體系微藍且透光性變差,說明形成了水包油型乳狀液。而不加表面活性劑的情況下,主要以單一流態(tài)為主,環(huán)狀流并未出現(xiàn),且當持水率增大至80%時,仍然為油包水流動,相轉(zhuǎn)換未發(fā)生,這與不加表活劑差別較大,說明加入表活劑后油水界面張力降低,使油水界面更易被拉長變形,且油水界面更易破碎,有利于油水相轉(zhuǎn)換的發(fā)生。
實際生產(chǎn)過程中會出現(xiàn)油井產(chǎn)量忽高忽低,地層周期性出水的現(xiàn)象,持水率大小影響油水兩相流動型態(tài)與乳狀液形成過程,流速大小影響油水混合物在垂直井筒中流動摩阻與O/W乳狀液形成快慢,對于重油而言,持水率越大、流速越快,原油與水更易形成水包油乳狀液,能更快起到減阻效果。因此,研究不同流速下,不同含水范圍加入降黏劑后油水兩相垂直管流流態(tài)特征對實際降黏生產(chǎn)指導具有重要意義。
在常溫常壓下,配置質(zhì)量分數(shù)為0.5%的曲拉通-100表面活性劑水溶液,改變持水率10%~70%,考察流速分別為0.03、0.06、0.09、0.15 m/s時油水兩相流動型態(tài)特征,結果見圖5~7。
圖5 流速為0.06 m/s時油水兩相流動型態(tài)Fig.5 Oil-water two-phase flow pattern at velocity of 0.06 m/s
圖6 流速為0.09 m/s時油水兩相流動型態(tài)Fig.6 Oil-water two-phase flow pattern at velocity of 0.09 m/s
圖7 流速為0.15 m/s時油水兩相流動型態(tài)Fig.7 Oil-water two-phase flow pattern at velocity of 0.15 m/s
從圖5~7看出:流速為0.06 m/s,持水率為20%~70%時,出現(xiàn)了泡狀流、彈狀流、環(huán)狀流、擾動流;流速為0.09 m/s,持水率為20%、30%時為泡狀流,持水率為40%時為彈狀流,持水率為50%時為環(huán)狀流,大于60%為擾動流;流速為0.15 m/s,持水率20%~30%為泡狀流,持水率大于40%為擾動流。說明流速越高,持水率越大,越易形成擾動流,擾動流下,長條環(huán)狀水流被打破分散,使水相與油相混合程度更高,更易形成水包油乳狀液,使原油黏度降低。因此在高產(chǎn)量超稠油降黏開采過程中可適當下調(diào)摻入綜合持水率,實現(xiàn)低成本降黏開采。
根據(jù)加入表面活性劑后不同流速、不同持水率下油水兩相垂直管流流態(tài)特征實驗結果,建立常溫常壓加劑后形成乳狀液之前油水兩相流動型態(tài),結果見圖8。
圖8 加劑條件下油水兩相流動型態(tài)Fig.8 Oil-water two-phase flow pattern with surfactant
從圖8看出,混合流速越低,越容易出現(xiàn)彈狀流、蠕狀流等大粒徑水泡流,隨流速增加,彈狀流與蠕狀流范圍逐漸變窄甚至消失,泡狀流、環(huán)狀流和擾動流等分散度更高的流動型態(tài)所占比例逐漸增加。由于界面張力的極大降低,且水相密度低于油相,油相對水相的浮力作用使水相的表觀速度大于油相的表觀速度,油水界面處存在相對速度差,因此在低界面張力下水相更易被拉長,形成長條狀的彈狀流、蠕狀流以及環(huán)狀流。油水混合速度越大,水相更易在油相中分散,形成細小泡狀流和擾動流,油包水流動更易轉(zhuǎn)變?yōu)樗土鲃印?/p>
在垂直管流裝置中,控制油相流速0.12 m/s,使油相在管道中循環(huán)流動,通過注入泵將0.5%表面活性劑水溶液注入油流參與循環(huán),使整個體系持水率為50%,觀察加入表面活性劑溶液后,在垂直管道中O/W乳狀液的形成過程見圖9。
從圖9看出,50%持水率的表面活性劑水溶液加入油相中后,先以連續(xù)的擾動流流動型態(tài)穿過油相;到2 min時,連續(xù)的水相逐漸分散、溶解,油珠穿插在連續(xù)的水相中,形成更明顯的攪拌狀流動;3 min時,水相繼續(xù)分散,水相與油相形成油水界面膜,此時泡狀流出現(xiàn)(水泡放大10倍圖像見圖10(a));5 min時,部分水相與油相形成水包油乳狀液,透光性變差,水相與形成的水包油乳狀液以環(huán)狀流的形式在垂直管道中上升流動;7 min時,更多的水包油乳狀液形成,少部分水相分散在水包油乳狀液中;8 min完全形成O/W乳狀液(10倍放大圖見圖10(b))。
圖9 垂直管道加劑條件下O/W乳狀液形成過程Fig.9 Process of O/W emulsion formation with surfactant at uptake
圖10 形成乳狀液過程放大圖像Fig.10 Zooming images of emulsion forming process
從圖10看出,泡狀流時,水相在連續(xù)的油相中形成一個一個的水泡,說明表面活性劑降低油水表面自由能,使水泡穩(wěn)定且均勻的分散在油相中;形成水包油乳狀液后,小油滴被水包裹,液滴尺寸明顯降低,分散度增加,穩(wěn)定的存在在連續(xù)相中。
由于實際原油不可視,可視觀察法難以得到實際原油油水兩相流動型態(tài),為探索能夠測量原油油水兩相流動型態(tài)的方法,采用4組平行電阻探針,連接電阻率測量儀,固定溫度50 ℃,在混合流速為0.06 m/s下,模擬油管76 mm管徑中23.68 m3/d流量,通過電腦實時測量并記錄不同持水率(30%、50%、70%)時,塔河原油在加入水溶性復合降黏劑SDG-2的條件下,水包油乳狀液形成過程中電阻值變化曲線,并分析該曲線所對應油水兩相流動型態(tài),結果見圖11~13。
從圖11看出,持水率為30%時,加劑后油水兩相流動型態(tài)隨時間變化規(guī)律為:彈狀流—環(huán)狀流—泡狀流—O/W乳狀液流動。彈狀流流型下,降黏劑水溶液以大段塞的形式存在于油相中,間歇式上升,電阻探針信號變化明顯,呈周期式波動;當達到環(huán)狀流型時,降黏劑水溶液穿過油相中心,電阻探針信號波動明顯,表現(xiàn)為一段低值波動,一段高值波動;繼續(xù)循環(huán)運行一段時間,油水兩相流動型態(tài)轉(zhuǎn)變?yōu)榕轄盍?,電阻探針信號在小幅度范圍?nèi)振蕩,且振蕩頻率密集,每一個水泡接觸到探針,都會使信號向高頻值波動;當形成水包油乳狀液后,油滴均勻分散在水相中,水相為連續(xù)相,此時4組電阻探針顯示的電阻值會有幾秒鐘的持續(xù)時間,且保持在高電阻值附近,不同時域存在電阻值的負脈沖波動,每次的負脈沖波動都說明有油滴穿過探針。從加入降黏劑水溶液到完全形成O/W乳狀液,共用時900 s。
圖11 30%持水率下形成O/W乳狀液過程探針信號Fig.11 Probe signal of O/W emulsion formed at a water holdup of 30%
圖12 50%持水率時原油與水形成乳狀液過程探針信號Fig.12 Probe signal of O/W emulsion formed at a water holdup of 50%
從圖12看出,50%持水率條件下,4組探針信號顯示整個體系形成乳狀液過程為:環(huán)狀流—泡狀流—水包油乳狀液。持水率從30%增大到50%,不存在彈狀流,這是因為持水率增大,由于降黏劑降低油水界面張力作用,水相更易被油相“拉長”,成為長條形而非圓球形,從而形成中心是水相、接近管壁處為油相的水中心環(huán)狀流,流動一段時間后,環(huán)狀流流型逐漸被打破,水相以水滴的形式分散在油相中,形成泡狀流,此時油水混合程度加大,繼續(xù)運行一段時間,整個體系形成O/W型乳狀液,從加入降黏劑到完全形成乳狀液共用時540 s。
圖13 70%持水率時原油與水形成乳狀液探針信號Fig.13 Probe signal of O/W emulsion formed at a water holdup of 70%
從圖13看出,70%持水率條件下,四組探針信號顯示加劑后油水兩相流動型態(tài)為:擾動流—水包油乳狀液。擾動流型實際可看作是泡狀流、蠕狀流、環(huán)狀流的混合流型,因此擾動流的探針信號是泡狀流、蠕狀流、環(huán)狀流探針信號的集合,表現(xiàn)為高電阻值、低電阻值上下波動幅度大,且呈現(xiàn)一定持續(xù)時間。70%持水率下,由于持水率增大,水相受油相徑向阻力增大,由于加入降黏劑的水溶液與油相界面張力降低,水相無法繼續(xù)保持環(huán)狀流,轉(zhuǎn)變?yōu)椴▌觿×业臄_動流,此流型下油水兩相更易形成O/W型乳狀液,整個過程共用時280 s。
從圖11~13看出,基于電阻探針法測量得到的油水乳狀液形成過程與透明可視觀察法得到的油水乳狀液形成過程基本相符,說明對于不可視的原油,采用電阻探針法測量油水兩相流動規(guī)律具有一定可行性,為原油-水兩相流動規(guī)律的在線監(jiān)測方法提供了理論依據(jù)。
在油水兩相流動過程中,當油包水流動轉(zhuǎn)變?yōu)樗土鲃?,即相轉(zhuǎn)化發(fā)生,容易導致流體流變特性的突變;同時,液相與管道壁面的接觸方式及何種液相接觸管壁都會隨流型的變化而變化,從而對管道的腐蝕程度造成不同程度的影響,并且在反相點或其附近會產(chǎn)生動量、傳熱傳質(zhì)和壓降的急劇變化。因此,相轉(zhuǎn)換是油水管道設計和運行中需要考慮的主要因素之一?;谄屏髂P蚚6,9-10,19],在運動波理論的基礎上通過分析分散相在連續(xù)相中的濃度分布隨持水率的變化關系,推導出油包水流型向水包油流型轉(zhuǎn)換的邊界方程,該方法已獲得國內(nèi)外學者的認可。本文中采用數(shù)值劃分流型的方法,考慮到輕油水和摻表面活性劑條件下重油水兩相流動在持水率0~100%之間發(fā)生了相轉(zhuǎn)換現(xiàn)象,對常溫常壓下這兩種流動體系發(fā)生相轉(zhuǎn)換時的流型轉(zhuǎn)換界限方程進行分析。
根據(jù)漂移流模型可得相分布因數(shù)Co與持水率Yw和混合速度um之間的關系為
uso/Yo=Coum+u∝(1-Yo)N,
(3)
令u∝(1-Yo)N=X,則
uso/Yo=Coum+X.
(4)
式中,uso為油相表觀速度;Yo為持油率;u∝為單一油滴的極限上升速度;um為混合上升速度;Co為相分布因數(shù);N為與粒徑有關的指數(shù)。
單一油滴的極限上升速度u∝可通過式(4)獲得,這樣每次固定一個Yo值,做相應的uso/Yo與um的相關關系圖,并對相關關系圖進行線性擬合就可得到Co及X,則
N=ln(X/u∝)/lnYo.
(5)
用這種方法可以得到Co與Yo,N與Yo的對應關系,并能看出,若要得到每一組Co、N與Yo的相關關系,需要在固定持水率Yo的情況下考察uso/Yo與um的相關關系。
通過對油水兩相流動型態(tài)特征及持水率實驗結果進行分析,將實驗中獲得的uso、um、u∝、Yo分別代入式(3)和式(4),進行線性回歸后,得到Co及N隨持水率Yw的變化關系,結果見圖14。
圖14 相分布因數(shù)Co指數(shù)N隨持水率Yw的變化Fig.14 Changes of Co and N with Yw
從圖14看出,當0.45 摻表面活性劑條件下重油水兩相: Co=0.747+1.89Yw-6.99Yw2+9.65Yw3- 3.49Yw4. (6) N=5.445-16.473Yw+87.539Yw2-156.783Yw3+ 93.542Yw4. (7) 把漂移模型式(3)變形為 um(1-YoCo)=uso+u∝YoYwN. (8) 再由持水率Yw的定義得 Kw=1-uso/um=YoCo+(u∝YoYwN)/um. (9) 當重油水兩相Yw約為 0.60時,Co=0.896 8,N=5.14,代入式(9),設置初始值并對方程擬合后,得到油包水流型向過渡流型轉(zhuǎn)化的邊界方程。 摻表面活性劑條件下重油水兩相: Kw=0.448+0.014 7um-1× ((ρo-ρw)gσρo-1)1/2(1-(-1+(1+ 4Mσ))1/2(2Mσ)-1)1/2. (10) 式中,ρo和ρw分別為油和水的密度;M為表征流體黏度的無量綱特性參數(shù);σ為油水界面張力;g為重力加速度。 式(10)定義了摻表面活性劑重油水兩相流型由水包油向過渡流型轉(zhuǎn)換的邊界,若給定流體性質(zhì),則u∝已知,就可以在Yw~um坐標系下按式(10)繪出曲線。當um趨于無窮大時,對摻表面活性劑下重油水兩相流動,Yw的漸近值分別為Kw∝=0.456,Kw∝=0.448。當摻表活劑重油水兩相Yw約為0.65時,流型開始向水包油轉(zhuǎn)化,此時重油水Co=1.371、N=5.65,設定初始值并對方程進行擬合后得到過渡流型向水包油流型轉(zhuǎn)化的邊界方程。 摻表面活性劑重油水: Kw=0.548+0.014 1um-1((ρo-ρw)gσρo-1)1/2(1- (-1+(1+4Mσ))1/2(2Mσ)-1)1/2. (11) 當um趨于無窮大時,Yw的漸近值分別為Kw∝=0.537和0.548。綜合式(10)和式(11)得到垂直上升管中油水兩相流流型轉(zhuǎn)換分界線,將前文通過觀察到的油水兩相流動型態(tài)實驗結果與模型得到的流型對比,結果如圖15所示(其中紅線為水包油流型,藍線為油包水流型)。 圖15 實驗與模擬得到的流型對比Fig.15 Comparison of experimental and simulated flow patterns 從圖16看出:對于摻表面活性劑重油水兩相流動,通過實驗觀察到的油包水泡狀流、油包水彈狀流、油包水段塞流幾乎全部都落在了油包水流型界限方程內(nèi),而水包油分散流、水包油泡狀流、水包油彈狀流和水包油段塞流也基本落在水包油流型界限方程內(nèi),吻合度達90%以上。油包水數(shù)值邊界和水包油數(shù)值邊界中間區(qū)域為過渡流型,其中擾動流(Churn)即為油包水流型向水包油流型轉(zhuǎn)換的過渡流型。說明流型界限方程能夠較好的將本文研究的不同類型油水流動型態(tài)進行劃分。基于該方法可對其余條件下流型轉(zhuǎn)換邊界進行數(shù)值辨識和劃分。 (1)加入表面活性劑后,由于油水界面張力的急劇降低,水相在油相中的分散度增加,油水界面更易發(fā)生變形,使油水泡狀流轉(zhuǎn)變?yōu)殚L條狀的彈狀流或蠕狀流,油水流動型態(tài)往往表現(xiàn)為2種或2種以上流態(tài)共存,如泡狀流和彈狀流、環(huán)狀流和擾動流共存。 (2)基于可視觀察和電阻探針法探究了加劑后油水乳狀液的形成規(guī)律,發(fā)現(xiàn)油水乳化過程為:大液滴的彈狀流→環(huán)狀流→擾動流→泡狀流→油水乳狀液,油水混合流速越大、持水率越高,擾動流越易發(fā)生,形成乳狀液所需時間越短。 (3)根據(jù)漂移流模型,建立了加劑后重油水兩相流動持水率預測和流動型態(tài)轉(zhuǎn)變界限方程。流型界限方程能夠較好地將實驗研究的不同類型油水流動型態(tài)進行劃分,基于該方法可對其余條件下流型轉(zhuǎn)換邊界進行數(shù)值辨識和劃分。4 結 論