動態熱條件下原油罐儲維溫過程主要影響因素分析

孫 巍，劉玉多，成慶林，趙立新，王志華

(1.東北石油大學 提高油氣采收率教育部重點實驗室，黑龍江 大慶 163318； 2.黑龍江省石油石化多相介質處理及污染防治重點實驗室，黑龍江 大慶 163318)

國家“十四五”現代能源體系規劃提出的關于實施重點行業領域節能降碳行動，要求加強工業領域節能和能效提升，推動重點行業節能改造，實現能耗“雙控”向碳排放總量和強度“雙控”轉變[1-2]。油田企業一直是能源消耗和碳排放大戶，每年用于油品維溫以及外輸的熱量占到系統耗熱量的85%以上，導致能源消耗量大，碳排放量增加[3-4]。目前，國內近49%的油田企業采用加熱盤管對具有高黏、易凝特性的原油進行維溫儲存，而儲罐邊界受外界環境波動較大，需根據不同季節變化對儲罐內部維溫條件進行相應地調節[5-8]。因此，有必要針對不同季節氣候對應儲罐維溫條件下的儲罐邊界油溫變化情況進行主要影響因素分析。

以往的學者都是在忽略盤管內部強制對流的基礎上，將盤管作為定熱源開展對儲罐內部油流狀態的研究。周凱等[9]、ZHAO等[10]將盤管的熱邊界條件作為恒溫邊界，在建立儲罐管式加熱二維理論模型的基礎上，研究了儲罐加熱過程中內部原油溫度場、流場分布情況。Sterling等[11-12]通過電加熱器對家用熱水儲罐進行加熱模擬，發現受加熱源影響，罐內流體形成密度差出現分層，在罐底區域出現了溫度較低的加熱死角。Magazinovi[13]、Ramadan[14]等采用數值模擬的方法對不同加熱管數下的水箱內部溫度場分布情況進行分析，并開展了不同加熱管數對水箱內部水溫變化、吸收熱能以及傳熱速率影響程度的研究，將儲罐所處動態熱環境作為定解條件處理來研究儲罐維溫過程中熱量傳遞過程。Kocijel等[15]在研究不同儲罐尺寸、加熱源距邊界位置對罐內流體溫度分布影響的基礎上，以壁面邊界傳熱系數變化為依據，分析了不同季節下儲罐邊界熱損失情況。孫巍等[16-18]通過建立外部動態熱環境以及原油變物性模型，從不同液位高度、盤管加熱結構以及外部邊界條件等角度，開展了由儲罐內部原油向外界傳熱規律影響的研究。王曉東[19]以原油溫度場的不均勻程度和加熱效率作為加熱效果的評價指標，定性地分析了不同加熱溫度下儲罐邊界熱量損失情況，并對不同類型、不用液位下的儲罐維溫溫度進行了優化調整。

筆者在綜合考慮外界動態熱環境、原油變物性參數、盤管動態熱源以及空間結構的基礎上，建立大型含蠟原油罐儲維溫過程三維理論模型，開展了對儲罐維溫過程中內部原油溫度場、流場分布規律的研究，提出了盤管周向效應對罐內含蠟原油對流傳熱的影響規律，劃分了外界環境影響區、內部維溫熱流影響區以及過渡區域；采用多元非線性回歸的方法，建立不同季節下的邊界區域主控影響因素模型，實現定量表征罐邊界溫度與動態熱環境、原油變物性參數之間的影響機制，為油庫生產管理實現降低能源消耗、減少碳排放的目標提供理論支持。

1 考慮動態熱環境下的原油罐儲維溫過程三維理論模型

1.1 儲罐維溫物理模型

原油罐儲維溫過程物理模型如圖1所示。從圖1可以看出，在原油儲罐維溫過程中，盤管內部加熱介質以強制對流的形式向管壁進行傳熱，極大地提高了盤管的能效利用率，促進了熱量向罐內油品的傳遞，盤管附近原油經溫升后密度減小，與其他區域油品形成密度差，在罐內形成自然對流。為便于研究含蠟原油罐儲維溫過程熱量傳遞特性，對模型作出如下簡化：①忽略原油內部的相變過程；②忽略儲罐內部構件對原油維溫效果的影響；③忽略外界風速對儲罐邊界油溫的影響。

1.2 控制方程

在已知原油物理參數的基礎上，通過建立相應的控制方程，對原油儲罐維溫過程中內部傳熱與流動規律進行研究，采用LES大渦模擬方法更好地對大尺度空間下儲罐內部形成的大渦流結構進行求解[20-25]。

①質量守恒方程：

(1)

式中：ρoil為油品密度，kg/m3；x、y、z分別為x、y、z方向坐標，m；u、v、w分別為x、y、z方向油流速度，m/s。

圖1 原油罐儲維溫過程物理模型圖Fig.1 Physical model diagram of temperature-maintaining process in crude oil tank storage(a)Global physical model of storage tank maintenance temperature；(b)Local physical model diagram of heating coil

②動量守恒方程：

(2)

(3)

(4)

式中：p為油品靜水壓力，Pa；μoil為油品動力黏度，Pa·s。

③能量守恒方程：

(5)

(6)

④大渦運動方程：

濾波函數：

(7)

式中：D為流體區域；x為空間坐標；G(x，x′)為濾波函數，決定了所求解的渦的尺度。

通過濾波函數處理瞬時狀態的Navier-Stokes方程及連續方程，得到大渦模擬的控制方程為：

(8)

(9)

(10)

根據Smagorinsky的基本模型，假定亞格子尺度應力為：

(11)

(12)

1.3 原油變物性模型建立

原油物性測量值與擬合值對比情況見圖2。由圖2可知，隨著原油溫度的升高，油品密度、導熱系數逐漸降低呈線性變化規律，黏度呈指數變化規律，受原油溫度變化影響較大，而原油比熱容呈現出先升高后降低的變化趨勢。通過測得的數據進行擬合，得到各參數與溫度之間的數學關系，如式(13)～(16)所示。

油品密度：

圖2 原油物性測量值與擬合值對比Fig.2 Comparison of the measured and fitted values of crude oil physical properties(a)Crude oil density；(b)Crude oil thermal conductivity；(c)Crude oil viscosity；(d)Crude oil specific heat

(13)

油品導熱系數：

λoil=-0.00007Toil+0.1366

(14)

油品比熱容：

(15)

油品黏度：

μoil=37.805e-0.097Toil

(16)

1.4 動態熱環境

隨著一天中光照角度的不斷變化，儲罐罐頂、罐壁在同一時刻下所接受的太陽輻射熱有所不同，導致儲罐邊界油溫呈震蕩波動式變化，因此，有必要針對儲罐不同邊界位置建立相應的外界動態熱環境理論模型[26-28]。

罐頂邊界條件：

①無太陽照射：

(17)

②有太陽照射：

(18)

罐壁邊界條件：

①無太陽照射：

r=rtank，0≤x≤htank，

(19)

②有太陽照射：

(20)

罐底邊界條件：

(21)

將等溫度場和零速度場作為原油儲罐維溫過程數值模擬的初始條件：t=0，T=T0(T0為原油初始溫度，℃)，u=v=w=0。

2 原油儲罐維溫過程的數值模擬與實驗驗證

采用尺寸為φ160 mm×8、距罐底1 m位置的4組加熱盤管，罐壁選用導熱系數為0.035 W/(m·℃)、厚度為0.08 m的保溫材料對罐內原油進行維溫，儲罐基礎設計參數如表1所示，儲罐動態邊界條件設計參數如表2所示。分別對盤管局部尺寸為0.08～0.02 m進行了網格劃分，對不同盤管尺寸下的儲罐內部熱流傳遞進行研究，結果如表3和圖3所示。0.04 m以下盤管網格尺寸原油溫度場模擬結果誤差小于0.6%，故最終選取0.04 m的盤管局部網格尺寸進行網格劃分，最終得到罐內原油網格數量為12640980，盤管內部熱水域網格數量為4581410的三維儲罐維溫模型。

表1 100000 m3浮頂儲罐基礎設計參數Table 1 Design parameters of 100000 cubic meters floating roof storage tank foundation

表2 儲罐熱邊界參數Table 2 Thermal boundary parameters of storage tanks

表3 局部網格尺寸誤差對比Table 3 Comparison of local grid size errors

圖3 儲罐模型網格劃分Fig.3 Grid division of storage tank model(a)Global grid division of storage tank；(b)Local grid division of coil

通過瞬時N-S方程計算原油儲罐內部的大尺度湍流運動，由于所采用的熱邊界條件、原油物性參數分別隨時間、油品溫度呈動態變化，故需使用UDF編程得以實現；壓力-速度耦合方面，運用速度進口和壓力出口來模擬熱水流經盤管過程。以1×105 m3大型原油儲罐為研究對象，運用VITO MTT系統對儲罐內部油溫變化的時實數據監測，通過與數值模擬結果進行分析對比，檢驗數值求解過程的準確性。罐內原油溫度場各測點實測數據與數值模擬結果相對誤差分析對比結果如表4所示。由表4可知，數值模擬結果符合罐內原油溫度場實際分布情況，而靠近罐頂區域油品受大氣溫度、太陽輻射影響較大，傳熱形式復雜多樣，導致19.58 m測點相對誤差最大為13.51%，測量值與數值模擬結果出現較大偏離。

表4 罐內原油溫度場實測結果與模擬結果對比Table 4 Comparison between measured and simulation results of crude oil temperature field in tank

3 原油罐儲維溫過程的傳熱與流動規律

含蠟原油罐儲維溫過程中，受盤管內部強制對流換熱影響，促進盤管區域油品熱量的吸收，導致油品密度變小、黏度降低，罐內原油形成密度差，從而引起油品間發生自然對流，促進了罐內油品的混合，儲罐內部油品運動跡線如圖4所示。由圖4可知，在大尺度空間下，吸收盤管熱量后的原油，沿罐壁位置開始上浮，形成了由盤管→罐壁→罐中心→罐底的整體大型渦流結構，盤管進出口間隔區域依托罐內大渦流結構，在黏滯力的作用下形成了流速為0.001～0.005 m/s的小型湍渦結構。

圖4 罐內油品速度跡線圖Fig.4 Velocity trace diagram of oil in tank

儲罐維溫過程中，罐壁位置形成了流速為0.012～0.015 m/s的湍渦結構，并由兩側逐步向儲罐中心區域發展，罐內油品混合程度提高，儲罐總體油溫上升，儲罐軸向位置油品熱流分布如圖5所示。由圖5可知，儲罐頂部直接與外部動態熱環境接觸，熱量交換頻繁，罐頂區域油品熱量耗散嚴重，罐壁位置受保溫層作用，油溫變化小，接近儲罐平均油溫，罐底位置距大渦流結構較遠，受罐內自然對流換熱強度弱，導致冷流體聚集在罐底中心位置形成低溫區。儲罐軸向不同位置油品溫度場、流場分布情況如圖6所示。由圖6可知，a區域受無盤管加熱區域形成不完全的小型渦流結構影響，罐內整體大渦流結構發生偏移，且a區域盤管周向效應最弱，導致罐底低溫區域厚度最大；而由儲罐中心a區向c、e區過渡中，兩股渦旋不斷靠近，盤管周向效應增強，罐底形成的低溫區域減小。

儲罐作為大尺度空間，罐內熱流由罐壁向儲罐內部油品進行換熱，其徑向位置油品熱流分布如圖7所示。由圖7可知，隨著熱流向罐中心接近，熱量逐漸散失，在儲罐中心出現不規則的低溫區域，邊界位置距離儲罐中心較遠，熱流傳遞較慢，油品溫度低于平均。儲罐徑向不同位置油品溫度場、流場分布情況如圖8所示。由圖8可以看出，a區域油品受盤管加熱影響，罐底對應盤管位置油品溫度高于其他區域；由a區到d區逐漸由罐內大渦結構下方轉變為上方區域，油品間對流傳熱強度提高，并在無盤管區域出現局部湍渦；由d區至e區，渦流強度明顯減弱，且e區受外界環境影響，熱量散失嚴重。

圖6 儲罐軸向各位置油品溫度場、流場云圖Fig.6 Nephograms of oil temperature field and flow field at each axial position of storage tank y/m：(a)0；(b)15；(c)30；(d)-15；(e)-30

圖7 儲罐徑向位置油品云圖Fig.7 Nephogram of oils in the radial position of storage tank(a)Radial temperature field nephogram；(b)Vector diagram of radial flow field

4 原油罐儲維溫過程主要影響因素分析

罐內油品溫度受動態熱條件、原油變物性影響呈現波動式變化規律，因此采用多元非線性回歸的方式，建立原油溫度場與各影響因素之間的數學關

圖8 儲罐徑向各位置油品溫度場、流場云圖Fig.8 Nephograms of oil temperature field and flow field at each radial position of storage tank z/m：(a)0；(b)5；(c)10；(d)15；(e)20

系，確定各因素對油品溫度影響的權重，如式(22)～(25)所示。

①數據進行歸一化處理：

(22)

式中：X*為歸一化后的樣本數據；max(x)為樣本數據的最大值；min(x)為樣本數據的最小值。

②建立熱量傳遞過程多元非線性回歸模型：

(23)

③給定回歸模型初始值，進行迭代求解：

(24)

式中：Ti為儲罐邊界油品溫度的模擬值；T′為儲罐邊界油品溫度的回歸值；SStot為總誤差，對于總誤差來源主要有2個：影響因素多樣，導致輸出值Ti多樣；隨機誤差SSres。顯然，隨機誤差SSres越小，R2越接近于1。

④建立原油儲罐邊界熱量傳遞主控影響因素數學模型：

(25)

儲罐全年用以保證內部原油安全運行的維溫條件隨季節轉變而發生變化，使得外部、內部因素在罐邊界形成的主要影響區域存在差異，因此，有必要對動態熱環境下儲罐維溫過程邊界油溫變化的主要影響因素進行研究，儲罐在不同季節下的維溫條件如表5所示。由表5可知，受季節變化影響，夏季外界環境溫度較高，儲罐內外溫差小，熱量散失少，內部所需維溫條件較低；而冬季儲罐邊界位置受外部低溫環境影響較大，熱量散失嚴重，有必要提高儲罐內部維溫條件，以確保儲罐安全穩定運行。

表5 不同季節儲罐維溫條件Table 5 Temperature-maintaining conditions of storage tanks in different seasons

4.1 罐頂油溫變化影響因素定量分析

由于外界動態熱環境溫度對罐頂位置油溫變化影響較大，因此選取罐頂位置有太陽輻射時間段進行數據分析，回歸結果如表6～表9所示。由表6～表9可知，在擬合方程中各變量次冪相同情況下，自變量自身系數可以代表對因變量邊界油溫變化的影響程度，通過計算得到各影響因素權重占比，由此劃分出外部、內部因素主要影響區以及過渡區域。

表6 春季罐頂位置油溫變化擬合方程Table 6 Fitting equation of oil temperature change at tank top in spring

表7 夏季罐頂位置油溫變化擬合方程Table 7 Fitting equation of oil temperature change at tank top in summer

根據罐頂位置各因素對邊界油溫變化的影響權重占比將罐頂位置劃分不同影響區域進行研究：Ⅰ(內部維溫熱流影響區)、Ⅱ(過渡區)、Ⅲ(外界環境影響區)，如圖9所示。由圖9(b)可以看出，夏季對區域Ⅰ油溫變化影響較大的因素為原油密度、比熱，權重占比高達80%以上，受大尺度空間自然對流影響較強；區域Ⅲ形成了以太陽輻射傳熱為主、與外界環境溫差形成的自然對流為輔，最大厚度為1.79 m的外界環境影響區。由夏季(見圖9(b))向冬季(見圖9(d))轉變的過程中，大氣溫度成為影響罐頂邊界油溫變化的主導因素，且受冬季罐內維溫熱流影響，使得罐頂多數油品處于受內、外因素綜合作用下，最大厚度為1.68 m的過渡區域Ⅱ。

表8 秋季罐頂位置油溫變化擬合方程Table 8 Fitting equation of oil temperature change at tank top in autumn

表9 冬季罐頂位置油溫變化擬合方程Table 9 Fitting equation of oil temperature change at tank top in winter

圖9 罐頂位置各因素對邊界油溫變化的影響權重占比Fig.9 Weight proportion of various factors of tank top position on boundary oil temperature change(a)Spring；(b)Summer；(c)Autumn；(d)Winter

4.2 罐壁油溫變化影響因素定量分析

罐壁位置受保溫層作用，外部動態熱環境對其內部油品溫度變化小，故選取整天數據進行研究，回歸方程如表10～表13所示。由表10～表13可計算得到各因素對罐壁區域油溫變化的影響權重，實現對罐壁位置油品內、外主要影響區域的定量分析。

表10 春季罐壁位置油溫變化擬合方程Table 10 Fitting equation of oil temperature change at tank wall position in spring

表11 夏季罐壁位置油溫變化擬合方程Table 11 Fitting equation of oil temperature change at tank wall position in summer

表12 秋季罐壁位置油溫變化擬合方程Table 12 Fitting equation of oil temperature change at tank wall position in autumn

表13 冬季罐壁位置油溫變化擬合方程Table 13 Fitting equation of oil temperature change at tank wall position in winter

罐壁位置各因素對邊界油溫變化的影響權重占比如圖10所示。由圖10可以看出，罐壁位置受保溫層影響，外部邊界條件對儲罐內部油溫變化影響較弱，罐內形成自然對流促進罐壁處油品熱量傳遞。其中，原油密度對其油溫變化影響權重高達50%以上，且油品黏度更容易受盤管升溫熱流的影響。而夏季(見圖10(b))儲罐所處外界大氣溫度接近于內部油品溫度以及受罐壁保溫層作用，儲罐內外幾乎不發生熱量交換。

4.3 罐底油溫變化影響因素定量分析

罐底位置直接與土壤接觸，受外界大氣、太陽輻射熱影響小，相關回歸方程如表14～表17所示。由表14～表17可知，同次冪擬合方程下的各影響因素自身系數可得到對罐底邊界油溫變化的權重占比，從而實現罐底油溫變化主要影響因素的定量分析。

根據罐底位置各因素對邊界油溫變化的影響權重占比，劃分出罐底位置外部、內部主要影響區域以及過渡區域，如圖11所示。由圖11可以看出，罐底位置油溫變化主要受內部原油物性影響較大，其中原油密度權重占比近50%，形成以自然對流傳熱為主的維溫熱流影響區；隨著外界土壤溫度降低，罐內外溫差變大，熱量散失嚴重，外部因素影響占比逐漸變大，邊界形成以導熱為主，厚度為0.13 m的外界土壤影響區Ⅰ。

圖10 罐壁位置各因素對邊界油溫變化的影響權重占比Fig.10 Weight proportion of various factors of tank wall position on boundary oil temperature change(a)Spring；(b)Summer；(c)Autumn；(d)Winter

表14 春季罐底位置油溫變化擬合方程Table 14 Fitting equation of oil temperature change at tank bottom in spring

表15 夏季罐底位置油溫變化擬合方程Table 15 Fitting equation of oil temperature change at tank bottom in summer

表16 秋季罐底位置油溫變化擬合方程Table 16 Fitting equation of oil temperature change at tank bottom in autumn

表17 冬季罐底位置油溫變化擬合方程Table 17 Fitting equation of oil temperature change at tank bottom in winter

圖11 罐底位置各因素對邊界油溫變化的影響權重占比Fig.11 Weight proportion of various factors of tank bottom position on boundary oil temperature change(a)Spring；(b)Summer；(c)Autumn；(d)Winter

5 結 論

(1)在大尺度空間儲罐內部，罐壁位置形成流速為0.012～0.015 m/s的湍渦結構，形成了由盤管→罐壁→罐中心→罐底的整體大型渦流結構，受無加熱盤管區域小型湍渦結構影響，使得大渦結構發生偏移。罐底位置受罐內自然對流換熱強度弱，導致冷流體聚集在罐底中心位置形成低溫區，且隨著兩側盤管距離漸近，盤管周向效應增強，罐底低溫區域明顯得到改善。

(2)罐頂油溫受季節變化影響較大，形成了以太陽輻射傳熱、與外界環境溫差形成的自然對流的最大厚度為1.79 m的外界環境影響區；受冬季罐內維溫熱流影響，大部分油品處在內、外因素綜合作用下，最大厚度為1.68 m的過渡區域；隨著外界土壤溫度降低，罐底位置形成了以導熱為主，厚度為0.13 m的外界土壤影響區；罐壁受保溫層影響，邊界形成了以原油自然對流為主的內部維溫熱流影響區，且原油各物性參數對其油溫變化影響比重波動較小。