盛程

（1．上海發(fā)電設(shè)備成套設(shè)計研究院，上海200240；2．華北電力大學(xué)核熱工安全與標準化研究所，北京102206）

自然循環(huán)CHF影響因素的灰色關(guān)聯(lián)度分析

盛程1，2

（1．上海發(fā)電設(shè)備成套設(shè)計研究院，上海200240；2．華北電力大學(xué)核熱工安全與標準化研究所，北京102206）

臨界熱流密度（CHF）對反應(yīng)堆安全起著至關(guān)重要的作用。自然循環(huán)條件下的CHF發(fā)生受到很多因素影響，目前對其特征的把握尚不完善。基于實驗獲得的常壓自然循環(huán)實驗數(shù)據(jù)，采用反映系統(tǒng)序列間正、負相關(guān)性的灰色關(guān)聯(lián)度量化模型，分析了加熱段入口焓值、流速和出口干度對自然循環(huán)CHF的影響。通過關(guān)聯(lián)度計算得到，入口焓值和出口干度與CHF呈負相關(guān)，工質(zhì)流速與CHF呈正相關(guān)。

自然循環(huán)；臨界熱流密度；灰色關(guān)聯(lián)度；影響因素

0 引言

由于受到傳熱惡化等的限制，工業(yè)設(shè)備中兩相沸騰的熱流密度難以得到進一步提高［1］。預(yù)測及預(yù)防換熱設(shè)備中出現(xiàn)的傳熱惡化，掌握關(guān)鍵因素的變化規(guī)律和影響程度，一直是熱工水力學(xué)領(lǐng)域的研究熱點。核電站發(fā)生例如全廠斷電（SBO）的嚴重事故后，堆芯余熱將完全依賴自然循環(huán)導(dǎo)出。但是由于一回路冷卻水流量較低，并且與強迫循環(huán)相比，自然循環(huán)會出現(xiàn)靜態(tài)流量漂移現(xiàn)象，平均流量在漂移點后將急劇下降，更容易發(fā)生傳熱惡化。CHF本身是一個復(fù)雜的非線性發(fā)展過程，而自然循環(huán)下的熱工水力現(xiàn)象不確定性更大，使得定性分析的難度增大。具有正、負相關(guān)性的灰色關(guān)聯(lián)分析方法是基于行為因子序列的微觀或宏觀幾何接近，以分析和確定因子之間影響程度的一種定量分析工具［2］?；趯嶒灁?shù)據(jù)，對影響CHF的入口焓值、工質(zhì)流速和出口干度這3個因素進行灰色關(guān)聯(lián)度計算，為進一步分析自然循環(huán)條件下CHF的發(fā)生機理奠定基礎(chǔ)。

1 實驗裝置和方法

如圖1所示，實驗回路主要包括實驗段、預(yù)熱器、穩(wěn)壓器和冷凝器等設(shè)備。實驗平臺高3.3 m，寬2 m。加熱通道的總長度為1 000 mm，橫截面為40 mm×5 mm的矩形。此外，在垂直加熱板的側(cè)面等間隔插入12根加熱管，單根功率均為2.5 kW。在實驗段出口附近距離內(nèi)壁面5 mm位置插入高溫?zé)犭娕?。實驗段正面為石英玻璃視窗，其余各面均包覆保溫材料，整個實驗段的熱效率為60%。

實驗壓力為常壓，工質(zhì)為去離子水。實驗過程中，在回路內(nèi)兩相自然循環(huán)流動形成后，以0.25 kW/min的速率逐漸提高實驗段功率，并開啟數(shù)據(jù)采集軟件對數(shù)據(jù)進行實時記錄。通過設(shè)定金屬保護溫度，當(dāng)加熱段壁溫或升溫速率達到設(shè)定的安全值時，自動切斷所有功率，以保護實驗設(shè)備安全。此時保存各項數(shù)據(jù)，記錄發(fā)生CHF時的功率，并計算出對應(yīng)的熱流密度值。

圖1 自然循環(huán)實驗回路

2 灰色關(guān)聯(lián)度模型

灰色關(guān)聯(lián)度可以定量反映系統(tǒng)變化過程中各因素之間的關(guān)聯(lián)程度及相對變化特征［2］。為了表現(xiàn)主行為與各因素間的正、負相關(guān)聯(lián)系，采用改進的灰色關(guān)聯(lián)度模型［3］。其構(gòu)建思想是通過不同序列間各時段的斜率比來反映相互間的相關(guān)程度。用變異系數(shù)表示2個序列斜率比值的分散程度。斜率比值越接近1，表示關(guān)聯(lián)性越好［4］。拓廣變異系數(shù)，定義廣義變異系數(shù)，考察斜率比的值是否在1附近。基于以上思路建立的模型［3］為

式中：X0為參考序列，系統(tǒng)主行為序列，X0=（x0（1），x0（2），…，x0（n））；Xi為比較序列，系統(tǒng)相關(guān)因素序列，在j時刻上為Xj（j），Xi=（xi（1），xi（2），…，xi（n））（i=1，2，…，m）；δ為變異系數(shù)；ξ為廣義變異系數(shù)；γ（X0，Xi）為X0和Xi的關(guān)聯(lián)度。具體定義見文獻［4］。

3 計算方法

3.1 數(shù)據(jù)預(yù)處理

在實際應(yīng)用中，由于各數(shù)據(jù)序列的量綱和數(shù)量級存在差異，要得到正確的結(jié)論并不容易。所以在關(guān)聯(lián)分析前需要對各序列數(shù)據(jù)進行無量綱化。通?？衫镁禑o量綱法進行處理。即

上述模型需將序列｛xi（j）｝始點零化，參考序列和各比較序列始點零化像［5］，分別為

3.2 計算流程

灰色關(guān)聯(lián)度的計算流程按文獻［4］執(zhí)行。

4 計算分析

4.1 影響因素灰色關(guān)聯(lián)度計算

采用作者在本課題中獲得的數(shù)據(jù)［6］，見表1，進行灰色關(guān)聯(lián)度分析。

利用表1的實驗數(shù)據(jù)，根據(jù)3.2節(jié)中文獻［4］采用的計算流程，在應(yīng)用公式（2）~（3）對數(shù)據(jù)進行預(yù)處理后，基于公式（1），得到如下所示的計算結(jié)果

表1 臨界熱流密度及相關(guān)因素值

4.2 灰色關(guān)聯(lián)度分析

根據(jù)4.1節(jié)計算結(jié)果可以得到：在自然循環(huán)條件下，加熱段入口焓值與qCHF呈負相關(guān)性，即隨著入口焓值的增加，qCHF逐漸減小；反映自然循環(huán)流速的入口雷諾數(shù)與qCHF呈正相關(guān)性，即隨著工質(zhì)流速的增加，qCHF增加；加熱段出口干度與qCHF呈負相關(guān)性，即隨著出口干度的增加，qCHF減小。根據(jù)此結(jié)論可進行如下分析。

入口焓值的增大，即一定壓力下入口過冷度的減小使通道內(nèi)工質(zhì)的溫度水平升高，液體在較小的加熱量下就能發(fā)生汽化。同時加熱段內(nèi)的流體密度差減小，自然對流效應(yīng)的減弱使通道的換熱能力變差。這兩方面共同作用，使CHF隨著入口焓值的增大而減小。

隨著流速增大，加熱表面的汽泡脫離頻率增大，脫離后汽泡的運動速度也加快，蒸汽很快流出通道，因此不易出現(xiàn)環(huán)狀流型和汽泡壅塞。此外，隨著更多的過冷液體進入加熱通道，加熱面能夠獲得更好的冷卻，壁面過熱度減小，蒸汽產(chǎn)生速率下降，因此在較大的加熱條件下CHF才會發(fā)生。

隨著加熱段出口干度的增加，通道出口處流體的焓值增大，液體蒸發(fā)所需的熱量減小，因而傳熱惡化發(fā)生時需要的加熱量減小。當(dāng)出口干度逐漸增大時，意味著蒸汽所占的比重增大，而液體所占的比重減小。當(dāng)出口干度較大，同時通道出口處的流型為彈狀或環(huán)狀流時，加熱量不需要太大就可以在較短時間內(nèi)將覆蓋加熱面的液膜層蒸干，并引起溫度突升。

5 結(jié)語

應(yīng)用一種可以體現(xiàn)序列間正、負相關(guān)性的改進型灰色關(guān)聯(lián)度模型，結(jié)合實驗數(shù)據(jù)，得到了自然循環(huán)條件下CHF相關(guān)因素中加熱段入口焓值、工質(zhì)流速和出口干度的影響程度量化值，并進行了分析。隨著實驗段入口焓值和出口干度的增加，臨界熱流密度均減小，體現(xiàn)了與CHF的負相關(guān)性；而隨著工質(zhì)流速的增加，臨界熱流密度增加，體現(xiàn)了與CHF的正相關(guān)性。通過合理應(yīng)用反映正負相關(guān)性的灰色關(guān)聯(lián)度模型，可以為衡量自然循環(huán)CHF相關(guān)因素的影響程度提供一種量化方法。

［1］盛程，周濤，李精精，等.窄通道自然循環(huán)臨界熱流密度的非線性分析［J］.原子能科學(xué)技術(shù)，2012，46（11）：1 330-1 335.

［2］鄧聚龍.灰理論基礎(chǔ)［M］.武漢：華中科技大學(xué)出版社，2002.

［3］曹明霞.灰色關(guān)聯(lián)度模型正負性問題的研究及其改進［J］.系統(tǒng)工程與電子技術(shù)，2008，30（6）：1 086-1 088.

［4］劉平，周濤，黃彥平.基于灰色關(guān)聯(lián)度的圓管內(nèi)流動CHF影響因素分析［C］.2010年第五屆綠色財富（中國）論壇，北京，2010.

［5］劉思峰，黨耀國，方志耕.灰色系統(tǒng)理論及其應(yīng)用［M］.北京：科學(xué)出版社，2004.

［6］盛程，周濤，張蕾，等.窄矩形通道自然循環(huán)流動停滯與臨界熱流密度研究［J］.核科學(xué)與工程，2013，33（1）：65-75.

Gray Correlation Degree Analysis of Influencing Factors in Natural Circulation

SHENG Cheng1，2
（1.Shanghai Power Equipment Research Institute，Shanghai 200240，China；2.Institute of nuclear thermal-hydraulic safety and standardization，North China Electric Power University，Beijing 102206，China）

Critical heat flux（CHF）plays a vital role in the safety of reactor.CHF in rectangular channels under natural circulation are affected by many factors whose features are still not grasped perfectly.Grey correlation degree calculation model which reflects the positive and negative relation of sequences is applied to analyze the influences of inlet enthalpy value，flow velocity，and outlet dryness on CHF under natural circulation conditions.Results showed that there were negative correlations between inlet enthalpy value or outlet dryness and CHF，while a positive correlation existed between flow velocity and CHF.

natural circulation；critical heat flux；gray correlation degree；influencing factor

TK124

A

1007－9904（2015）02－0015－03

2015-01-07

盛程（1985），男，工程師，主要從事反應(yīng)堆熱工水力方面的研究。

國家自然科學(xué)基金項目（50976033）資助。