普子恒 趙 偉 熊鈺瑤 史星濤 方春華

（1.三峽大學(xué) 電氣與新能源學(xué)院, 湖北 宜昌 443002；2.湖北省輸電線路工程技術(shù)研究中心, 湖北 宜昌443002）

隨著社會經(jīng)濟的發(fā)展,用電量逐年增長,城市電纜化率越來越高；然而由于電纜破損、絕緣老化等原因可能發(fā)生火災(zāi),造成巨大事故損失.據(jù)統(tǒng)計,2013年至2016年因電纜燃燒引發(fā)的火災(zāi)約占電氣火災(zāi)數(shù)量的60%,并有逐漸增長的趨勢[1].為降低電纜火災(zāi)發(fā)生率,減小火災(zāi)事故影響,需對電力電纜燃燒特性進(jìn)行分析,獲得電纜火災(zāi)蔓延、煙氣發(fā)展等規(guī)律,并采取相應(yīng)的火災(zāi)防護策略.

電力電纜燃燒試驗具有一定的破壞性,存在成本高、不可逆的缺點,增加了研究難度；近年來流體力學(xué)理論發(fā)展迅速,電纜燃燒模擬方法也得到了進(jìn)一步的完善.通過大量的研究驗證[2-5],電纜燃燒數(shù)值模擬已得到了廣泛應(yīng)用[6-9].模擬電力電纜燃燒的過程主要分為固相熱解與氣相反應(yīng)兩個部分,固相熱解模擬電纜受熱后可燃成分的形成,氣相反應(yīng)模擬可燃成分在空氣進(jìn)行燃燒反應(yīng)釋放熱量的過程.彭玉輝等對狹長通道電纜火災(zāi)進(jìn)行了研究,采用了簡單熱解即設(shè)置恒定熱釋放火源,氣相反應(yīng)選擇乙烯氣相反應(yīng),重點分析了煙氣運動規(guī)律以及火災(zāi)探測器在電纜火災(zāi)中的適用性[10]；芬蘭國家技術(shù)研究中心Matala等[11]為研究核電站的消防系統(tǒng)配置,對核電站電纜燃燒進(jìn)行了模擬,模型的氣相反應(yīng)定義了聚氯乙烯氣相反應(yīng)；高明旭[12]在研究綜合管廊火災(zāi)安全性中,設(shè)置了簡單熱解中的t2模型,氣相反應(yīng)定義了聚氯乙烯氣相反應(yīng).王明年等[13]對綜合管廊內(nèi)的電纜設(shè)置了簡單熱解,通過仿真模擬得到火災(zāi)功率的發(fā)展與典型區(qū)域的溫度變化,并研究了電纜倉內(nèi)結(jié)構(gòu)參數(shù)對火災(zāi)發(fā)展的影響.針對復(fù)合結(jié)構(gòu)的電力電纜燃燒模擬,多數(shù)學(xué)者對固相熱解采用了簡單熱解的方式,而氣相反應(yīng)采用了基礎(chǔ)的甲烷、乙烯氣相反應(yīng)或者電纜護套材料聚氯乙烯的氣相反應(yīng)；且沒有過多地開展不同設(shè)置下電纜燃燒模擬的分析與對比工作,僅將研究重點放在了電纜燃燒模擬后的結(jié)果分析與實際應(yīng)用[14-16].由于電力電纜具有多層結(jié)構(gòu),且不同層所采用的材料均不同,實際燃燒較為復(fù)雜,上述仿真模擬與試驗測量值存在一定差異,需改進(jìn)燃燒模擬的方法,才能得到更準(zhǔn)確的電纜火災(zāi)蔓延和煙氣發(fā)展規(guī)律.

本文首先對簡單熱解與復(fù)雜熱解進(jìn)行了仿真模擬,分析了相同受熱條件下的固相質(zhì)量熱解損失率差異.其次,基于電纜燃燒產(chǎn)物含量對聚氯乙烯反應(yīng)方程進(jìn)行修正,并與聚氯乙烯氣相反應(yīng)以及基礎(chǔ)的甲烷、乙烯氣相反應(yīng)進(jìn)行了模擬,分析了相同熱釋放速率條件下的溫度場分布規(guī)律.最后,根據(jù)試驗建立不同的電纜燃燒模型進(jìn)行模擬驗證.

1 燃燒模型

1.1 電力電纜燃燒模擬基本原理

電力電纜具有多層結(jié)構(gòu),由護套層、絕緣層和金屬導(dǎo)體3部分組成.電力電纜的燃燒過程相較于一般的電線電纜更為復(fù)雜,主要包括兩個過程：(1)在熱源的作用下,電力電纜存在明顯的溫度差異,最外側(cè)護套層相較于內(nèi)部溫升更快,開始發(fā)生熱解反應(yīng),且較為劇烈；而內(nèi)部的絕緣層通過熱傳導(dǎo)的方式傳遞熱量,熱解反應(yīng)較為緩慢.熱解產(chǎn)生不同的揮發(fā)性氣體,既包含可燃成分,也包含非可燃成分,同時還會產(chǎn)生碳煙顆粒.(2)熱解產(chǎn)生的可燃性氣體在初期主要由護套層產(chǎn)生,氣體進(jìn)入氣相與固體表面的空氣接觸,進(jìn)行氣相氧化反應(yīng),釋放熱量,形成燃燒[17]；發(fā)展階段的氣體將由護套層與絕緣層共同產(chǎn)生,燃燒更為劇烈；燃燒后期的氣體主要由可燃物燃燒后所形成的焦炭產(chǎn)物熱解產(chǎn)生.依據(jù)實際燃燒機理,電纜的燃燒模擬劃分為固相熱解與氣相反應(yīng),是兩者耦合所形成的,燃燒示意圖如圖1所示.

圖1 電纜燃燒示意圖

1.2 固相熱解

電力電纜的護套層一般采用含鹵族元素的材料,例如聚氯乙烯,起到防腐蝕、防磨損的保護作用；絕緣層采用聚乙烯、交聯(lián)聚乙烯等材料,起到絕緣的作用.溫度的不同和材料的差異性導(dǎo)致電纜各部分產(chǎn)生了多種熱解反應(yīng),每種熱解反應(yīng)速率與產(chǎn)物依據(jù)材料自身的熱解屬性與外界條件共同決定.當(dāng)電纜燃燒模擬場景的邊界條件涉及到熱傳導(dǎo)和熱解反應(yīng)時,可以分為簡單熱解與復(fù)雜熱解方式進(jìn)行相關(guān)參數(shù)的設(shè)定.

簡單熱解的方式需設(shè)置電纜各個結(jié)構(gòu)層的傳熱系數(shù)、密度以及比熱容等熱物理參數(shù),并應(yīng)用一維熱傳導(dǎo)方程計算內(nèi)部的熱量傳遞,方程見式(1)；對于材料的熱解反應(yīng)可進(jìn)行簡化,將電纜看作等厚度薄板進(jìn)行研究,不涉及固體內(nèi)部熱解的具體過程,只根據(jù)電力電纜整體燃燒熱解性能設(shè)置單位面積熱釋放速率與燃點.

式中：ρs為密度；cs為比熱容為內(nèi)熱源項.

復(fù)雜熱解的方式除設(shè)置電纜各個結(jié)構(gòu)層的傳熱系數(shù)、密度與比熱容等基本的熱物理參數(shù)用于導(dǎo)熱計算外,針對電纜不同結(jié)構(gòu)熱解情況還應(yīng)設(shè)置各自的熱解動力學(xué)參數(shù),并使用Arrhenius公式計算反應(yīng)速率：

式中：Wi為第i個熱解反應(yīng)速率；Ai為第i個反應(yīng)指前因子；Ei為第i個反應(yīng)活化能；R為通用氣體常數(shù),其值為8.314 kJ/(kmol·K)；Ts為反應(yīng)時的溫度(K)；N為反應(yīng)級數(shù)；ρi為當(dāng)前固相種類的質(zhì)量濃度；ρ0為邊界處固相表面的初始密度；XO2為模擬反應(yīng)速率受局部氧濃度影響.

1.3 氣相反應(yīng)

電纜熱解產(chǎn)生的可燃?xì)怏w種類繁多,與空氣中的氧氣發(fā)生多聯(lián)式化學(xué)反應(yīng),所形成的產(chǎn)物具有多態(tài)、數(shù)量多的特點,模擬過程中為實現(xiàn)對氣體物種跟蹤,對所有氣體物種求解質(zhì)量輸運方程將大大增加計算成本,且極易導(dǎo)致計算的不收斂.可燃?xì)怏w與氧氣反應(yīng)釋放熱量的動力學(xué)機理目前尚不完善,為簡化模擬過程,避免復(fù)雜化學(xué)反應(yīng)的計算,采用混合控制燃燒模型.模型將反應(yīng)分成了燃料、空氣與產(chǎn)物3部分來描述氣體燃料和氧氣的湍流燃燒,每個部分均作為一種“集總物種”,模擬中將每個集總物種看作一個整體進(jìn)行輸運并求解跟蹤,見式(3)：

式中：燃料化學(xué)式中B為除碳、氫與氧之外的其他元素,一般為氮或鹵族元素.空氣集總物種的化學(xué)方程式系數(shù)為：

vs與vCO的化學(xué)方程式系數(shù)通過式(5)由該燃料燃燒生成的CO與碳煙顆粒的量產(chǎn)率推出：

式中：WF為燃料相對分子質(zhì)量；Wn為燃燒產(chǎn)物的相對分子質(zhì)量；yn為燃燒產(chǎn)物的量產(chǎn)率.

在模擬場景中,電纜中不同材料熱解的多種可燃?xì)怏w均由氣相反應(yīng)中的燃料所替代；然而不同材料的燃燒熱各有不同,為使不同種類材料燃燒釋放相對應(yīng)的熱量,依據(jù)材料的燃燒熱與定義的氣相反應(yīng)的燃燒熱之間的差異,調(diào)整氣體燃料質(zhì)量產(chǎn)生率.即材料質(zhì)量損失率仍由熱解公式確定,但熱解單位時間所產(chǎn)生的氣體燃料質(zhì)量由材料燃燒熱與氣相反應(yīng)燃燒熱的比值決定：

式中：為材料熱解進(jìn)入氣相反應(yīng)的燃料；mf為熱解公式計算得出的燃料質(zhì)量；Hf為材料的燃燒熱；Hc為氣相反應(yīng)的燃燒熱.

從理論上分析,氣相反應(yīng)設(shè)置不同方程所建立的不同燃燒模型在模擬過程中,雖然不同反應(yīng)方程的燃燒熱有所不同,熱解產(chǎn)生的氣體也由反應(yīng)方程的燃料所替代,但燃料在進(jìn)入氣相之前會將所定義反應(yīng)方程的燃燒熱與電纜材料燃燒熱代入公式(6)進(jìn)行計算.如此一來,根據(jù)計算結(jié)果就能得到產(chǎn)生相同熱量的燃料質(zhì)量數(shù)值,從而保證了不同氣相反應(yīng)的燃燒模型在模擬同一電纜燃燒過程中產(chǎn)生的熱量相同.

2 固相熱解模擬與分析

2.1 固相熱解的模擬設(shè)置

為研究不同方式設(shè)置的固相熱解對模擬電纜燃燒的影響,將復(fù)雜熱解與文獻(xiàn)使用較多的簡單熱解[12-16]在相同受熱條件下進(jìn)行質(zhì)量熱解損失率的對比分析.選取典型電力電纜進(jìn)行分析,外護套為聚氯乙烯材料,絕緣層為聚乙烯材料.建立電纜熱解樣本,樣本厚度設(shè)定1 mm,形成熱薄層,便于測試固相模型的熱解情況.按照簡單熱解與復(fù)雜熱解分別設(shè)置參數(shù),加熱溫度設(shè)定為700℃,計算區(qū)域為0.8 m×0.8 m×0.4 m的長方體空間,網(wǎng)格尺寸為2 cm×2 cm×1 cm,總網(wǎng)格數(shù)為64 000個,熱解模擬時間為60 s.

2.2 固相熱解質(zhì)量損失對比分析

在相同的溫度加熱條件下,固相熱解設(shè)置的兩種方式包括簡單熱解與復(fù)雜熱解模擬的質(zhì)量損失情況,如圖2所示.

圖2 不同固相熱解設(shè)定方式的電纜質(zhì)量損失率

在初期的恒定受熱條件下,曲線變化見圖2右上角.復(fù)雜熱解設(shè)置的電纜模型在溫度上升的過程中,依據(jù)Arrhenius公式進(jìn)行熱解,出現(xiàn)了持續(xù)的質(zhì)量損失；然而簡單熱解由于溫度未達(dá)到設(shè)定的燃點,沒有出現(xiàn)變化,直到5.4 s后,溫度達(dá)到了燃點,質(zhì)量損失出現(xiàn)大幅度的提升.在熱解的中期,溫度達(dá)到穩(wěn)定,簡單熱解在12 s左右達(dá)到了0.45 g/s的質(zhì)量損失峰值；復(fù)雜熱解在10 s左右達(dá)到0.3 g/s的質(zhì)量損失峰值后,表面質(zhì)量密度下降,熱解速率逐漸減小.在熱解的后期,由于電纜質(zhì)量減少,簡單熱解質(zhì)量損失率下降快,在30 s左右結(jié)束熱解；復(fù)雜熱解設(shè)定了絕緣層以及護套層的材料熱解參數(shù),在30 s后主要進(jìn)行的是絕緣層的熱解,熱解質(zhì)量損失在0.15 g/s左右,直至50 s后,僅剩下少部分的電纜樣本,熱解速率逐漸下降,在60 s降為0.01 g/s.

在受熱的條件下,固相的質(zhì)量損失通過式(6)轉(zhuǎn)換至燃料,并進(jìn)入到氣相中燃燒形成熱量,即固相的質(zhì)量損失率的大小代表了熱釋放速率的數(shù)值.在固相熱解模擬中,可以看出,簡單熱解雖然通過燃點與單位面積熱釋放速率的設(shè)定可以在熱解前期與復(fù)雜熱解保持一致,但隨著溫度的升高與熱解的深入進(jìn)行,兩種方式的質(zhì)量損失率即熱釋放速率出現(xiàn)了顯著的差異.相對而言,復(fù)雜熱解基于電纜的不同結(jié)構(gòu)設(shè)置了不同的熱解參數(shù),可能更能反應(yīng)出真實情況下的電纜熱解情況.

3 氣相反應(yīng)模擬與分析

3.1 基于電纜燃燒產(chǎn)物含量修正系數(shù)的氣相反應(yīng)

電力電纜主要燃燒物為護套層、絕緣層,且兩者材料大多不一致,而眾多學(xué)者采用的聚氯乙烯氣相反應(yīng)或基礎(chǔ)的甲烷、乙烯氣相反應(yīng)[18]并沒有考慮到絕緣層的燃燒,需對其進(jìn)行修正.針對典型電力電纜,外護套為聚氯乙烯材料,絕緣層為聚乙烯材料,對其進(jìn)行氣相反應(yīng)的修正.首先根據(jù)聚氯乙烯的燃燒產(chǎn)物量產(chǎn)率,按照式(5)計算,配出相應(yīng)的聚氯乙烯氣相反應(yīng),量產(chǎn)率見表1,反應(yīng)方程見表2中的2號.

表1 燃燒產(chǎn)物量產(chǎn)率

在電力電纜燃燒過程中,絕緣層的聚乙烯材料燃燒也會產(chǎn)生二氧化碳、一氧化碳以及碳煙等產(chǎn)物,而鹵代物氯化氫則是由聚氯乙烯產(chǎn)生,可假設(shè)電纜燃燒是聚氯乙烯與聚乙烯按照比例燃燒的過程,根據(jù)表1中典型電纜燃燒產(chǎn)物氯化氫的量產(chǎn)率,即可建立關(guān)于燃燒比例的方程：

式中：X為護套層燃燒的比例；WF1、WF2為護套層、絕緣層材料的相對分子質(zhì)量；WHCl為燃燒產(chǎn)物氯化氫的相對分子質(zhì)量；yHCl為電纜燃燒產(chǎn)物氯化氫量產(chǎn)率；v1(HCl)為護套層氣相反應(yīng)方程氯化氫的系數(shù).

通過式(7)可計算出聚氯乙烯的比例X為0.4,則可得到電纜燃燒的氣相燃料分子式為C2H3.6Cl0.4,相對分子質(zhì)量為41.8.最后根據(jù)表1中典型電纜的碳煙量產(chǎn)率與一氧化碳量產(chǎn)率,按照式(5)計算出相應(yīng)的系數(shù),配平氣相反應(yīng)方程,見表2中的1號.

表2 氣相反應(yīng)方程

3.2 氣相反應(yīng)的模擬設(shè)置

為研究不同氣相反應(yīng)對模擬電纜燃燒的影響,將基于電纜燃燒產(chǎn)物含量修正的氣相反應(yīng)與基礎(chǔ)的甲烷、乙烯氣相反應(yīng)以及聚氯乙烯氣相反應(yīng)在相同熱釋放速率條件下進(jìn)行溫度場的對比分析.具體氣相反應(yīng)方程見表2.

計算區(qū)域為0.8 m×0.8 m×0.8 m的正方體空間,采用均勻網(wǎng)格,網(wǎng)格尺寸為2 cm,總網(wǎng)格數(shù)為6.4×104個.空間底面中心設(shè)置0.4 m×0.4 m的燃料入口,模擬電纜熱解產(chǎn)生可燃?xì)怏w被不同氣體燃料替代后釋放出相同熱量的情況.計算區(qū)域除底面外其余為開放邊界在計算區(qū)域,同時在Y=0.4 m處設(shè)置溫度監(jiān)測面,記錄氣體燃料燃燒溫度,如圖3所示.

圖3 模擬空間示意圖

3.3 燃料溫度場分布對比

圖4為4種氣體燃料依據(jù)定義的氣相反應(yīng)燃燒產(chǎn)生8 k W熱釋放率的溫度場,由于不同氣體燃料燃燒熱不同,形成8 k W熱釋放率所需燃料分別為：500.0 mg/s的C2H3.6Cl0.4、487.8 mg/s的C2H3Cl、160.0 mg/s的CH4與169.8 mg/s的C2H4.

圖4 不同氣體燃料燃燒產(chǎn)生8 kW熱釋放率的溫度場

圖中燃燒功率較小,氣體流動平穩(wěn),波動較小,形成的溫度等值線圍繞著燃燒中心區(qū)域,該區(qū)域氧氣充足,燃料充分反應(yīng),在4種不同氣體燃料中C2H3Cl形成的溫度場中心區(qū)域溫度最高,達(dá)到560℃,C2H3.6Cl0.4中心區(qū)域溫度僅為406℃.在Z軸方向上,C2H3.6Cl0.4的100℃等溫線位于0.7 m高度以上,125℃至400℃等溫線較為分散.C2H3Cl燃料的100℃等溫線位于0.6 m高度左右,由于4種燃料的熱釋放速率相同,C2H3Cl中心區(qū)域溫度高,因此125℃至560℃等溫線分布范圍會更小.CH4與C2H4燃料燃燒形成的等溫線較為一致,其125℃至500℃的等溫線分布比C2H3.6Cl0.4更為分散,比C2H3Cl更為集中；由于C2H4燃料中心區(qū)域相較于CH4溫度更高,在200℃以上的等溫線分布范圍更小.通過網(wǎng)格單元的數(shù)據(jù)計算,超過50℃區(qū)域的相同位置點,C2H3.6Cl0.4氣相反應(yīng)與C2H3Cl、CH4、C2H4氣相反應(yīng)平均相差了12.33℃、8.83℃、11.06℃,中心區(qū)域最高溫度減小了154℃、108℃、135℃.

圖5為4種氣體燃料燃燒產(chǎn)生40 k W熱釋放率的溫度場.此時由于計算區(qū)域熱釋放速率大,溫度高,中心區(qū)域位置有所提高.

圖5 不同氣體燃料燃燒產(chǎn)生40 k W熱釋放率的溫度場

4種燃料中,C2H3Cl中心區(qū)域溫度仍為最高,達(dá)到了928℃,超過400℃的等溫線更為集中.CH4、C2H4的中心溫度分別達(dá)到了832℃、897℃,整體的等溫線僅出現(xiàn)了較小的差異,差異體現(xiàn)為：300℃與400℃的溫度范圍內(nèi),CH4覆蓋范圍更大；在800℃以上的中心區(qū)域,C2H4溫度更高.C2H3.6Cl0.4隨著熱釋放速率的增加,其中心區(qū)域溫度僅為523℃,但300℃至400℃的中低溫區(qū)域覆蓋范圍變得更大.通過網(wǎng)格單元的數(shù)據(jù)計算,超過100℃區(qū)域的相同位置點,C2H3.6Cl0.4氣相反應(yīng)與C2H3Cl、CH4、C2H4氣相反應(yīng)平均相差90.18℃、79.89℃、83.75℃,中心區(qū)域最高溫度減小了405℃、309℃、374℃.

3.4 溫度場分布差異分析

燃燒模型基于混合限制,所定義集總物種的快速化學(xué)反應(yīng)和燃料的平均化學(xué)源項使用渦流耗散概念(EDC)建模[19].假設(shè)計算單元內(nèi)的局部氣體存在兩種狀態(tài)：未混合或完全混合,混合部分將進(jìn)行反應(yīng).當(dāng)氣體燃料初始為非預(yù)混狀態(tài),燃料質(zhì)量變化速率為：

計算單元內(nèi)混合區(qū)中的燃料成分變化為：

通過式(8)與式(9)可得到基于EDC模型的燃料質(zhì)量變化速率：

單位體積內(nèi)的燃料質(zhì)量變化速率決定了單位體積內(nèi)的熱釋放速率：

式中為單位體積熱釋放速率；Δh為燃料的燃燒熱.

4種燃料消耗單位的質(zhì)量所需空氣質(zhì)量各不同,燃料本質(zhì)是與空氣中的氧氣反應(yīng),為便于分析,考慮燃料與所需氧氣質(zhì)量的關(guān)系,參數(shù)見表3.

表3 不同燃料消耗氧氣參數(shù)表

在燃料過量的燃燒中心區(qū)域,即式(10)最小值取到第2種情況下,s越小,單位體積內(nèi)燃料質(zhì)量變化速率更快.C2H3Cl的s僅為0.783,分別是C2H3.6Cl0.4、CH4、C2H4的0.43倍、0.20倍以及0.23倍,單位體積內(nèi)燃料質(zhì)量變化速率則是另外3種的2.34倍、5.11倍與4.38倍；表2中所列出的燃燒熱中C2H3Cl是C2H3.6Cl0.4、CH4、C2H4燃燒熱的1.03倍、0.33倍以及0.35倍,盡管C2H3Cl的燃燒熱小于CH4、C2H4的燃燒熱,根據(jù)式(10)溫度場中心區(qū)域的單位體積熱釋放速率,C2H3Cl仍將大于其余3種的單位體積熱釋放速率,因此C2H3Cl燃料中心區(qū)域溫度比其他3種高,并且隨著燃料流量增多,中心區(qū)域溫度差距更為明顯.在整體熱釋放率設(shè)定一致的條件下,C2H3Cl在燃燒中心區(qū)域具有更多的熱釋放速率,Z軸方向上較高位置的熱釋放率會偏小,即出現(xiàn)溫度高的等溫線分布集中,分布在高處的等溫線溫度較低的現(xiàn)象.C2H3.6Cl0.4單位體積熱釋放速率最低,故在圖中表現(xiàn)為等溫線的分散分布.

4 試驗仿真算例

通過上述的分析,可以看出不同的固相熱解設(shè)定方式在相同受熱條件下,熱解質(zhì)量損失率即熱釋放速率存在差異性；不同的氣相反應(yīng)在相同熱釋放速率條件下,形成的溫度場分布也有所不同.因此,根據(jù)實際的電纜燃燒試驗,將復(fù)雜熱解與基于電纜燃燒產(chǎn)物含量系數(shù)修正的1號氣相反應(yīng)相結(jié)合,建立電纜燃燒模型；同時將復(fù)雜熱解與表2中的2、3、4號氣相反應(yīng)、簡單熱解與表2中的1、2號與3、4號氣相反應(yīng)相結(jié)合,建立8種電纜燃燒模型,進(jìn)行了仿真對比的驗證.

4.1 試驗介紹

電纜燃燒試驗選取美國CHRISTIFIRE項目中電纜水平燃燒試驗[20],其數(shù)據(jù)全面可靠,適于進(jìn)行仿真模型的對比驗證.電纜結(jié)構(gòu)參數(shù)見表4,電纜燃燒試驗如圖6所示.

表4 電纜結(jié)構(gòu)參數(shù)表

圖6 電纜燃燒試驗

試驗采用3個托盤,間距30 cm,每個托盤長2.4 m,寬0.45 m.托盤上排列著試驗電纜.在托盤下方20 cm處放置一個0.3 m×0.3 m的燃燒器,用于點燃電纜,即試驗所設(shè)置的火源,火源的功率約為40 k W,試驗過程中火源持續(xù)燃燒,直至試驗結(jié)束關(guān)閉火源.

4.2 不同模型的電纜燃燒仿真對比

根據(jù)試驗設(shè)置建立仿真模型,由于試驗中整體電纜燃燒的功率發(fā)展具有代表性,相較于其他物理量能夠更直觀反映仿真的準(zhǔn)確性,8種不同的電纜燃燒模型的燃燒功率變化與試驗的燃燒功率變化對比如圖7所示.

圖7 不同仿真模型與試驗燃燒功率對比圖

圖7中,8種電纜燃燒模型的燃燒功率曲線發(fā)展不同.采用復(fù)雜熱解與簡單熱解的燃燒模型,兩者之間整體釋放的總熱量存在差異,但在采用了相同的固相熱解方式的條件下,釋放的總熱量基本相同.

在固相采用簡單熱解方式的條件下,2號為C2H3Cl燃料的氣相反應(yīng),由于其在相同熱釋放速率形成的溫度場中心區(qū)域溫度最高,電纜在受熱條件下熱解產(chǎn)生燃料速率最快,燃料反應(yīng)的熱量一部分提高了環(huán)境溫度,一部分再次反饋給電纜,導(dǎo)致電纜熱解面積增大與熱解速率增快,曲線呈現(xiàn)出了燃燒功率的快速發(fā)展.但是在溫度達(dá)到設(shè)定的燃點后,熱解速率并不會隨著溫度的提供繼續(xù)增加,導(dǎo)致不同的氣相反應(yīng)形成的功率峰值差距不大.1號的C2H3.6Cl0.4功率曲線發(fā)展最慢,3號的CH4與4號的C2H4氣相反應(yīng)形成的溫度場相似,因此燃燒曲線除發(fā)展速率4號的C2H4略快于3號的CH4之外,其余較為一致.2號的C2H3Cl氣相反應(yīng)大約在1 200 k W,而1號的C2H3.6Cl0.4與3、4號的CH4、C2H4燃料氣相反應(yīng)均為1 100 k W.

在固相采用復(fù)雜熱解方式的條件下,由于熱解反應(yīng)受到溫度的影響較為顯著,溫度的不斷提高同時也促進(jìn)熱解速率的不斷增快.3號的C2H3Cl燃料氣相反應(yīng)燃燒功率峰值最高,且與其他3種氣相反應(yīng)拉開了明顯的差距.功率曲線在800 s就達(dá)到1 900 k W,在1 200 s左右基本燃燒完畢.3號的CH4與4號的C2H4兩種燃料之間的燃料熱與單位質(zhì)量消耗的氧氣質(zhì)量等系數(shù)較為接近,根據(jù)式(10)計算CH4燃料單位體積熱釋放速率小于C2H4燃料,造成了燃燒功率曲線上升階段出現(xiàn)較晚,增長速率更小,即CH4燃料在上升階段釋放了更少的熱量,在整體熱量大致相同的情況下,功率曲線衰減階段也更慢.此外,也可以看出采用復(fù)雜熱解相較于簡單熱解,更容易受到氣相反應(yīng)的影響.

復(fù)雜熱解與基于電纜燃燒產(chǎn)物含量修正系數(shù)的氣相反應(yīng)相結(jié)合的電纜燃燒模型,即1號的C2H3.6Cl0.4燃料氣相反應(yīng)與復(fù)雜熱解的仿真燃燒功率與試驗功率發(fā)展較為吻合.仿真與試驗的上升階段接近一致,試驗中電纜燃燒功率峰值為770 k W,仿真中燃燒功率高于試驗,達(dá)到820 k W；在衰減階段仿真燃燒功率下降速率更快,在2 300 s燃燒功率降為0,試驗曲線衰減過程在1 400 s后燃燒功率出現(xiàn)下降速率先增快后減慢的情況,下降曲線較為平緩,直至3 200 s左右燃燒完畢.仿真與試驗的燃燒功率誤差主要來自于燃燒后期功率的下降階段,整體誤差約為16.4%.造成誤差的主要原因是仿真中對電纜結(jié)構(gòu)與固相熱解步驟進(jìn)行了簡化,試驗中電纜實際燃燒進(jìn)行多種復(fù)雜熱解；在功率上升階段,主要是由電纜護套與絕緣材料燃燒產(chǎn)生大量熱量；在功率下降階段,電纜燃燒生成的產(chǎn)物例如殘留的碳繼續(xù)熱解反應(yīng)產(chǎn)生少量熱量.

5 結(jié) 語

本文提出了復(fù)雜熱解與基于電纜燃燒產(chǎn)物含量系數(shù)修正的氣相反應(yīng)相結(jié)合的電纜燃燒模型,分別對固相熱解與氣相反應(yīng)進(jìn)行了深入分析,并選取了試驗進(jìn)行驗證對比,得到了以下結(jié)論：

1)固相熱解的復(fù)雜熱解與簡單熱解在相同的受熱條件下,質(zhì)量損失率即熱釋放速率各不相同.針對由不同材料所組成、具有多層結(jié)構(gòu)的電力電纜而言,復(fù)雜熱解能隨著外界溫度的變化,對不同層的熱解速率分別進(jìn)行實時的調(diào)整,更符合實際的熱解情況.

2)不同的氣相反應(yīng)在相同條件下形成的溫度場具有差異性,氣相反應(yīng)中燃料單位質(zhì)量消耗的氧氣質(zhì)量以及燃燒熱參數(shù)將影響單位體積熱釋放速率.消耗單位質(zhì)量燃料所需的氧氣質(zhì)量越小,燃燒熱越大,則單位體積熱釋放速率更高,形成的溫度場高溫區(qū)域更為集中,數(shù)值更高.

3)復(fù)雜熱解與基于電纜燃燒產(chǎn)物含量系數(shù)修正的氣相反應(yīng)相結(jié)合的電纜燃燒模型可對多層結(jié)構(gòu)的電力電纜燃燒進(jìn)行仿真模擬,仿真結(jié)果更符合試驗功率.兩者誤差約為16.4%,誤差主要來自電纜燃燒的后期,沒有模擬實際燃燒形成的產(chǎn)物例如碳化物等的再次燃燒.

電纜熱解過程模擬仍需深入研究,以便建立更完善的仿真模型,進(jìn)一步提高仿真的準(zhǔn)確性.