張在孝 黃啟玉 張 汛 單錦旭 朱祥瑞 張 雨
(1. 中海石油(中國)有限公司深圳分公司 廣東深圳 518054; 2. 中國石油大學(北京)城市油氣輸配北京市重點實驗室 北京 102249)
流花16-2油田海底原油管道在不保溫輸送輕質原油的過程中,在海水與管道壁面的大溫差強制對流換熱及輕質原油的高析蠟點的作用下,海底原油管道的結蠟問題較為嚴重。海底管道蠟沉積一直是原油管道流動保障迫切關心的問題[1-2],一般而言,C17—C55的正構烷烴被認為是石蠟[3]。蠟在原油中的溶解度是溫度的函數,當油溫低于析蠟點時,蠟分子就會從原油中析出[4]。當累積析蠟量質量分數達到1%~2%時,蠟晶便可在原油中形成膠凝結構,將原油束縛其中[5-6]。海底管道環境溫度常低于輸送原油的析蠟點,因而在原油的輸送過程中,由于管壁處蠟分子溶解度的降低,在蠟分子徑向濃度梯度的作用下,蠟沉積物則容易在管壁處形成[7]。
蠟沉積機理主要包括分子擴散、剪切彌散、布朗擴散及重力沉降,其中分子擴散機理被大多數學者所接受[8-9]。影響蠟沉積物沉積速率的影響因素較多[10-14],包括油流溫度、環境溫度、油流與管壁溫度差、流速、原油組成以及沉積時間。一些學者以局部管道為對象開展研究,認為隨著沉積時間的增長,管壁結蠟層厚度的增加會減小原油的對外散熱能力,導致蠟沉積速率減小,蠟沉積厚度緩慢增大[14-16]。然而對于管道全線,尤其是大溫差海底管道全線,局部管段蠟沉積速率的減小無法代表管道全線的結蠟情況。因此,本文利用室內環道實驗裝置對流花16-2油田原油進行蠟沉積規律的探討,根據實驗結果建立流花16-2油田海底管道蠟沉積預測模型,模擬研究沉積時間對管道沿線蠟沉積速率、蠟沉積厚度以及單位時間管道全線總結蠟量的影響,并開展理論分析,為今后海底管道清管提供一定的理論指導。
1.1.1流花16-2油田原油基礎物性
流花16-2油田海底管道輸送原油的基礎物性參數情況:20 ℃時原油密度810.0 kg/m3,凝點0 ℃,析蠟點28.7 ℃,含蠟量5.84%,其溶解度系數曲線和黏溫曲線如圖1所示。

圖1 流花16-2油田管輸原油溶解度系數曲線以及黏溫曲線
1.1.2環道實驗裝置
利用室內小型環道對流花16-2油田原油蠟沉積規律開展研究,實驗裝置如圖2所示。該裝置主要由實驗管路系統、恒溫水浴系統、數據采集系統和空氣吹掃系統組成。原油經螺桿泵后在測試管道沉積,測試段管壁溫度有恒溫水浴系統控制,管壁溫度為海水環境溫度。實驗結束后,先用空氣吹掃系統清除管內剩余原油,再將測試段拆卸進行升溫加熱,使管壁沉積物融化后掉落至燒杯進行稱重。

圖2 流花16-2油田管輸原油蠟沉積環道實驗裝置示意圖
海底管道蠟沉積模型的通用模型為[17]
(1)

在不同溫度、流速、油壁溫差以及沉積時間的條件下,流花16-2油田原油蠟沉積速率實驗結果如表1所示。

表1 流花16-2油田原油蠟沉積速率實驗結果
將流花16-2油田原油蠟沉積速率的實驗參數代入式(1)并對實驗結果進行擬合,建立流花16-2油田海底管道蠟沉積速率模型為
(2)
流花16-2油田海底管道共計23.1 km,管道外徑φ254 mm,壁厚9.55 mm,保溫海底管道沿線環境最低溫度7.9 ℃,0~5.8 km段的設計總傳熱系數69 W/(m2·℃),5.8~23.1 km段的設計總傳熱系數104 W/(m2·℃),設計輸量5 406 m3/d,啟輸溫度112 ℃。
根據流花16-2油田海底管道基礎參數和原油基礎物性參數計算管道運行參數的變化,從而預測管道沿線蠟沉積速率及蠟層厚度的變化。
1) 管壁處徑向溫度梯度。
根據熱平衡關系,不同管段的徑向溫度梯度可表示為[18]
(3)
其中
(4)

2) 管道沿線原油黏度變化及管壁剪切率。
由于輸送油品含蠟量低可近似為牛頓流體,油品黏度隨溫度變化不明顯,因此可用Arrhenius方程表示[19]:
μ=AeEa/(RT)
(5)
式(5)中:Ea為黏性流動活化能,J/mol;A為指前因子或頻率因子,Pa·s;R為氣體常數,取8.314 J/(mol·K);T為絕對溫度,K。
牛頓流體紊流狀態下的管壁剪切率為[18]
(6)
3) 管道沿線總傳熱系數變化。
在管壁蠟沉積的影響下,蠟沉積后管道沿線的總傳熱系數可表示為[20]
(7)
式(7)中:Kw為蠟沉積后管道沿線總傳熱系數,W/(m2·℃);λw為結蠟層的導熱系數,W/(m2·℃);δ為結蠟層厚度,m。
4) 單位時間管道全線總結蠟量變化。
管道全線結蠟日增加量可表示為
(8)
式(8)中:Wall為單位時間管道全線總結蠟量,m3;δi,j為沉積j天后管道第i段結蠟層厚度,m;li,j為沉積j天后管道第i段管段長度,m;Δt為時間,d;N為管道全線劃分的總段數。
2.2.1預測結果
流花16-2油田海底管道在設計時按2個分段管線取不同的總傳熱系數,然而實際上該系數是連續變化的。因此,本文在預測該油田管道沿線蠟沉積分布(圖3)時,選取3個管道初始總傳熱系數K1、K2、K3,值分別為69、90和104 W/(m2·℃),以使結果更接近實際情況。

圖3 流花16-2油田海底管道不同沉積時間內管道沿線的結蠟情況分布
由圖3可知,隨著沉積時間的增長,流花16-2油田海底管道在3種總傳熱系數條件下管道沿線蠟沉積速率、結蠟厚度均呈現相似的變化規律:①蠟沉積速率和沉積厚度均呈現出先增加到峰值,后逐漸降低的趨勢;②蠟沉積速率的峰值位置隨沉積時間增長逐漸向管道末端移動,峰值大小逐漸降低;③蠟沉積厚度的峰值位置隨沉積時間增長逐漸向管道末端移動;④隨著沉積時間的增長,峰值位置前的管段蠟沉積速率變低,峰值位置后管段蠟沉積速率卻逐漸升高;⑤管道全線最大蠟層厚度逐漸增大,峰值位置前管段的蠟層厚度增加幅度要小于后管段。
然而,管道蠟沉積速率峰值大小的降低并不能代表管段全線蠟沉積速率的降低。因此本文提出了單位時間管道全線總結蠟量的概念,用于表述管道全線蠟沉積速率的大小。如圖4所示,單位時間管道全線總結蠟量隨著沉積時間的增長而逐漸增大。這一現象表明,管壁的結蠟保溫作用并不能有效減緩管道蠟沉積問題,相反,管道全線每天的結蠟增加量逐漸增多。為了保障管道正常運行,通常以2 mm的蠟沉積厚度作為清管周期的確定方法[16]。由圖3可知,流花16-2油田海底管道在運行3~5 d后管道全線的最大蠟沉積厚度都接近2 mm,因此推薦將管道全線最大蠟層厚度2 mm或管道全線總結蠟量10 m3的沉積時間作為清管周期(3~5 d),以保障管道的安全運行。

圖4 流花16-2油田海底管道單位時間管道全線總結蠟量隨沉積時間變化
2.2.2預測結果分析
由式(3)、(4)和(7)可知,在蠟沉積的過程中,沉積在管壁的蠟層會減小管段的總傳熱系數,提高管段油溫,降低單位管長溫降,從而使該管段蠟沉積速率減小。管道環境溫度、輸送油品溫度及兩者的溫差均影響著管道的蠟沉積。在流花16-2油田海底管道深水不保溫、強制對流換熱的條件下,進入管道的油流在很短的距離即可達到環境溫度,而這段距離即為管壁結蠟的主要區域。根據圖5可知,蠟沉積后管道沿線總傳熱系數Kw的最小值位于管道沿線最大蠟層厚度處,當最大蠟層處管段總傳熱系數減小,后管段的油壁溫差就會增大,因此蠟沉積速率峰值位置逐漸向管道末端移動。隨著沉積時間的增加,管道沿線蠟層厚度逐漸增大,導致蠟沉積后管道沿線總傳熱系數隨沉積時間逐漸下降,從而導致蠟沉積速率峰值大小逐漸降低。

圖5 流花16-2油田海底管道沿線總傳熱系數和溫降變化
根據單位時間管道全線總結蠟量的定義,式(8)還可以表示為

(9)
根據式(9),單位時間管道全線總結蠟量與管道蠟沉積速率及里程圍成的面積成正比。如圖6所示,在1 d和30 d的沉積時間內,管道沿線蠟沉積速率與里程圍成的面積可大致分為Ⅰ、Ⅱ、Ⅲ3個區域,Ⅲ區的面積(SⅢ)大于Ⅰ區的面積(SⅠ)表明30 d沉積時間的單位時間管道全線總結蠟量大于1 d沉積時間的。在海底管道不保溫輸送過程中,管壁結蠟保溫作用使得管道沿線總傳熱系數減小,進而導致管道沿線蠟沉積速率峰值的后移,使得管道主要結蠟區域變寬,蠟沉積速率與距離圍成的面積變大,即單位時間管道全線總結蠟量隨著沉積時間逐漸變大。同樣,隨著管道總傳熱系數的增大,一方面海水與管道壁面強制對流換熱加劇,導致管道結蠟位置向管道前端移動,管道沿程最大結蠟速率峰值增大;另一方面使管道沿線溫度梯度增大,導致油流溫度快速接近環境溫度,即主要結蠟區域變窄。因此,管道總傳熱系數越大,單位時間管道全線總結蠟量越小。

圖6 流花16-2油田沉積時間1 d和30 d的管道沿線蠟沉積速率分布
1) 流花16-2油田海底管道結蠟速率較快,在3~5 d的運行時間內,管道沿線最大蠟層厚度即可達到2 mm左右。隨著沉積時間的延長,管壁結蠟層的保溫作用使管道沿線總傳熱系數減小,管道蠟沉積速率峰值降低,且位置向管道末端方向移動。
2) 在深海不保溫、海水與管道間強制對流換熱的條件下,管壁的結蠟保溫作用使管道沿線最大蠟層厚度位置處的總傳熱系數迅速降低,導致后管段油壁溫差的升高,隨著沉積時間的增長,蠟沉積區域逐漸向后管段遷移,擴大管道全線蠟沉積區域,增加單位時間管道全線總結蠟量。
3) 鑒于流花16-2油田海底管道結蠟速率快、單位時間管道全線總結蠟量隨運行時間的增加而逐漸升高的特性,推薦以管道沿線最高蠟層厚度2 mm或管道全線總結蠟量10 m3的沉積時間,即3~5 d作為清管周期,以保障海底管道的輸送安全。