肖惠仁，姚禹歌，張 縵，呂俊復(fù)，張 揚(yáng)，祁傳西，陶世健，胡昌華，王振山

(1.清華大學(xué) 能源與動(dòng)力工程系，北京 100084; 2.陜西彬長(zhǎng)新民塬發(fā)電有限公司，陜西 西安 710018; 3.青島嘉能海諾電力設(shè)備有限公司,山東 青島 266300)

過(guò)去的20 a，我國(guó)經(jīng)濟(jì)水平高速增長(zhǎng)，帶動(dòng)了能源行業(yè)的快速崛起，人們生活質(zhì)量得到很大提升，與此同時(shí)環(huán)境污染問題也變得愈發(fā)不可忽視。在嚴(yán)峻的環(huán)境壓力下，我國(guó)推行了極其嚴(yán)格的污染物超低排放標(biāo)準(zhǔn)。煤燃燒發(fā)電由于其污染物排放量大，面臨著嚴(yán)峻考驗(yàn)。相比于其他煤燃燒發(fā)電技術(shù)，循環(huán)流化床(CFB)作為一種可大規(guī)模商業(yè)化推廣的清潔燃燒技術(shù)，擁有著煤種適應(yīng)性強(qiáng)、污染物排放低、負(fù)荷調(diào)節(jié)簡(jiǎn)單等特點(diǎn)，對(duì)于我國(guó)的能源清潔高效利用有著重要意義[1-5]。

灰控閥是CFB鍋爐的一個(gè)重要部件，主要應(yīng)用于爐膛排渣以及回料器到外置式換熱器的物料分配[6-9]，其運(yùn)行可靠性及耐用性直接影響到鍋爐的運(yùn)行情況以及電廠效益。CFB鍋爐的灰控閥有氣動(dòng)閥和機(jī)械閥兩類，目前市面上使用較多的氣動(dòng)閥有L閥、U閥、V閥、J閥等[10-14]。理論上，氣動(dòng)閥能更好的控制和調(diào)節(jié)物料流量，但研究發(fā)現(xiàn)各股流化風(fēng)的作用非線性甚至非單調(diào)，實(shí)際調(diào)節(jié)效果不甚理想[15-16]。錐形閥作為使用最廣的一種機(jī)械閥，具有實(shí)際運(yùn)行效果好、物料流量控制穩(wěn)定、調(diào)節(jié)簡(jiǎn)單等優(yōu)點(diǎn)，在大容量循環(huán)流化床鍋爐得到廣泛應(yīng)用[17-19]。錐形閥可以布置在外置式換熱床與回料器之間，調(diào)節(jié)從回料器進(jìn)入外置式換熱床的熱灰流量，從而控制外置式換熱床的換熱量，同時(shí)也對(duì)爐膛內(nèi)部溫度起到調(diào)控作用;錐形閥也常被用于排渣系統(tǒng)，控制排渣流量。

錐形閥的調(diào)節(jié)性能對(duì)CFB鍋爐的設(shè)計(jì)和運(yùn)行十分重要。目前，針對(duì)錐形閥的相關(guān)研究較少，其設(shè)計(jì)、運(yùn)行多是經(jīng)驗(yàn)性的，缺乏可靠的理論與研究數(shù)據(jù)。筆者設(shè)計(jì)并搭建了冷態(tài)實(shí)驗(yàn)臺(tái)，測(cè)量不同錐形閥開度、流化風(fēng)速、初始床料高度以及閥體標(biāo)稱尺寸條件下的閥流量，分析以上幾個(gè)影響因素對(duì)錐形閥內(nèi)部流動(dòng)的影響，探究錐形閥的流量特性，以期為錐形閥的設(shè)計(jì)與運(yùn)行控制提供依據(jù)。

1 孔口流動(dòng)系統(tǒng)

圖1為一種典型的錐形閥結(jié)構(gòu)示意圖，包括閥體、閥芯和電動(dòng)驅(qū)動(dòng)機(jī)構(gòu)。閥體為中空結(jié)構(gòu)，上部開有清掃孔通道，防止錐形閥堵塞結(jié)渣等問題出現(xiàn)。在實(shí)際工業(yè)應(yīng)用中，通過(guò)電動(dòng)驅(qū)動(dòng)結(jié)構(gòu)來(lái)調(diào)節(jié)閥芯的位置，從而調(diào)節(jié)錐形閥內(nèi)部的通流面積，進(jìn)而達(dá)到調(diào)節(jié)錐形閥物料流量的目的。實(shí)際實(shí)驗(yàn)使用的錐形閥如圖1(c)所示，僅包括閥體和閥芯部分，而不包括清掃孔和電動(dòng)驅(qū)動(dòng)機(jī)構(gòu)等對(duì)錐形閥流動(dòng)特性無(wú)影響的配件。

顆粒流過(guò)錐形閥，本質(zhì)上是顆粒通過(guò)孔口的流動(dòng)。關(guān)于顆粒通過(guò)孔口系統(tǒng)的流動(dòng)，已經(jīng)有很多相關(guān)研究。秦霽光等[20]只考慮重力作用，對(duì)顆粒自由下落流過(guò)小孔進(jìn)行了分析，結(jié)合實(shí)驗(yàn)結(jié)果，提出孔口流率Gp為

Gp=0.24g0.5(d0-dp)2.5ρp

(1)

式中，g為重力加速度，g=9.81 m/s2;d0為孔口流通直徑;dp為顆粒平均粒徑;ρp為顆粒密度。

郭天民等[21]分析了豎直下降管中的顆?？卓诹鞒銮闆r，將顆粒料柱視為連續(xù)介質(zhì)，利用伯努利方程，得到孔口流率為

(2)

式中，Cd為常數(shù)，主要考慮物料由孔口流出時(shí)的收縮現(xiàn)象;ε為顆?？紫堵?H為豎直管的高度。

BEVERLOO等[22]與ZHANG[23]分析了孔口兩側(cè)無(wú)壓差、順重力情況下顆?？卓诹鲃?dòng)情況。BEVERLOO等[22]發(fā)現(xiàn)孔口表面自由應(yīng)力對(duì)顆粒通過(guò)孔口系統(tǒng)會(huì)有影響，并對(duì)顆粒通過(guò)孔口的流動(dòng)邊界進(jìn)行了修正，認(rèn)為在顆粒流過(guò)孔口時(shí)，孔口最外側(cè)周邊區(qū)域無(wú)顆粒流過(guò)，且無(wú)顆粒流經(jīng)區(qū)域大小與顆粒直徑有關(guān)，因此增加了一個(gè)與顆粒形狀相關(guān)的修正系數(shù)K，認(rèn)為有效流通直徑為d0-Kdp。ZHANG[23]在觀測(cè)孔口附近顆粒流動(dòng)時(shí)，發(fā)現(xiàn)孔口附近氣泡的不穩(wěn)定性會(huì)對(duì)顆粒流率有影響，因此增加了一項(xiàng)與孔口附近顆粒間摩擦力有關(guān)的修正因子τ:

Gp=Cd(1-ε)(d0-Kdp)2.5ρpg0.5τ0.5

(3)

SPINK等[24]發(fā)現(xiàn)細(xì)顆粒與粗顆粒在重力作用下的流動(dòng)規(guī)律明顯不同。對(duì)于粗顆粒，顆粒對(duì)流體助力影響不敏感，一般不出現(xiàn)下料困難問題，顆粒下料流率隨著顆粒尺寸的增大而緩慢減小;但對(duì)于細(xì)顆粒，隨著顆粒粒徑減小，下料流率先上升到最大之后再逐步降低。

CREWDSON等[25]認(rèn)為細(xì)顆粒下料時(shí)，受到氣固流體動(dòng)力學(xué)作用力，這種間隙流體作用力相當(dāng)于在顆粒上加了一個(gè)容重，并以垂直壓力梯度dp/dr的形式來(lái)體現(xiàn)這一作用力，即

(4)

式中，θ為料倉(cāng)錐角。

由于忽略了顆粒間黏性力，因此式(4)只適用于黏聚性較小的粉體顆粒。式(4)也可用于孔口兩側(cè)有壓差時(shí)的情況。

目前，使用最廣泛的孔口流率計(jì)算公式為L(zhǎng)EUNG公式[26]，可適用于實(shí)際工業(yè)過(guò)程中大多數(shù)的顆?？卓谙到y(tǒng)流動(dòng)情況，且計(jì)算誤差較小。LEUNG公式視氣固混合物為無(wú)黏性流體，由托里拆利定理及力平衡方程導(dǎo)出:

(5)

式中，Gp0為無(wú)壓差時(shí)的顆?？卓诹髀?Δp0為孔口兩側(cè)壓差。

DAVIDSON等[27]研究了錐形料倉(cāng)的粉體下料情況，并提出了Hourglass理論，Hourglass理論適用于具有一定黏聚性的粉體:

(6)

式中，Kp為與顆粒性質(zhì)相關(guān)的系數(shù)。

前人關(guān)于顆?？卓谙到y(tǒng)的研究基本都是針對(duì)顆粒垂直下落及等流通截面面積的情況，但實(shí)際情況中，包括錐形閥實(shí)際使用時(shí)，往往是屬于顆粒傾斜下落流過(guò)變截面通道，本文將通過(guò)對(duì)錐形閥顆粒流動(dòng)進(jìn)行研究，簡(jiǎn)單探討顆粒傾斜下落通過(guò)變截面通道時(shí)的流動(dòng)特性。

在前人提出的經(jīng)驗(yàn)公式中，易見通流面積對(duì)閥流量影響作用較大。當(dāng)實(shí)驗(yàn)工況和顆粒性質(zhì)給定時(shí)，閥流量只與閥的通流面積有關(guān)，而對(duì)于錐形閥，閥的通流面積由閥芯位置(閥芯位移)來(lái)調(diào)節(jié)和決定。定義閥開度為閥芯相對(duì)于閥閉合時(shí)移動(dòng)的位移x，通流面積定義為最小通流截面面積，對(duì)應(yīng)的流通長(zhǎng)度為d。以閥體標(biāo)稱尺寸為30 mm的錐形閥為例，計(jì)算得到錐形閥位移x同閥通流面積S之間的函數(shù)關(guān)系為

S=(83.54-3d)d

(7)

對(duì)S(x)函數(shù)進(jìn)行線性擬合，得到對(duì)應(yīng)的擬合函數(shù)為

S=20.758x

(8)

在x∈(0,25 mm)時(shí)，平均誤差為4%，因此在本文考慮和計(jì)算的范圍內(nèi)，可以認(rèn)為S與x之間近似呈線性關(guān)系。

2 實(shí)驗(yàn)系統(tǒng)

實(shí)驗(yàn)所用物料為玻璃微珠，其顆粒球形度為0.85。通過(guò)馬爾文粒度分析，得到床料顆粒的顆粒粒度分布如圖2所示，平均粒徑d32為171.2 μm，中位粒徑d50為150.87 μm。

圖2 床料顆粒粒徑分布情況Fig.2 Particle size distribution of bed material

通過(guò)實(shí)驗(yàn)測(cè)量得到了床料顆粒的真實(shí)密度和堆積密度，分別為2.468 8 g/cm3和1.534 5 g/cm3，從而可知床料顆粒的密堆積體積分?jǐn)?shù)為0.621 6，空隙率為0.378 4，屬于Geldart-B類粒子。

為確定返料器內(nèi)部物料的流態(tài)，需得到臨界流化風(fēng)速，并將之與返料器內(nèi)部風(fēng)速進(jìn)行對(duì)比[28]。筆者采用關(guān)聯(lián)式計(jì)算與實(shí)驗(yàn)法來(lái)得到臨界流化風(fēng)速umf，計(jì)算所用關(guān)聯(lián)式[29]為

(9)

對(duì)于實(shí)驗(yàn)所用的平均粒徑為171.2 μm的玻璃微珠，計(jì)算得到其臨界流化風(fēng)速為0.074 m/s。實(shí)驗(yàn)通過(guò)下行法得到其臨界流化風(fēng)速為0.07～0.09 m/s。實(shí)際返料器內(nèi)部流化風(fēng)速為0.2 m/s，大于臨界流化風(fēng)速，返料器內(nèi)部物料處于流化狀態(tài)，從而可將返料器簡(jiǎn)化為一個(gè)帶錐形閥出口的流化床進(jìn)行實(shí)驗(yàn)研究。

如圖3所示，實(shí)驗(yàn)裝置由長(zhǎng)方體流化床、帶布風(fēng)板的風(fēng)室、錐形閥、風(fēng)機(jī)、給料器和料桶組成。流化床截面積為500 mm×100 mm，高度為1 000 mm。流化床頂部開口出風(fēng)，流化床底部為風(fēng)室，風(fēng)室中布置了兩塊布風(fēng)板，保證風(fēng)能被均勻給入到流化床中，流化風(fēng)的風(fēng)量調(diào)節(jié)通過(guò)安裝在風(fēng)機(jī)上的變頻器來(lái)實(shí)現(xiàn)。流化床前墻安裝有進(jìn)料漏斗，通過(guò)給料器和進(jìn)料漏斗來(lái)實(shí)現(xiàn)穩(wěn)定持續(xù)給料。實(shí)驗(yàn)中通過(guò)變頻器控制給料機(jī)的給料流量，使得給料流量與錐形閥出料流量大致相等，保持流化床內(nèi)部物料量平衡。流化床后墻安裝有錐形閥，錐形閥底部出料，出料口下方放置有料筒，通過(guò)測(cè)量流化穩(wěn)定后一定時(shí)間內(nèi)從錐形閥進(jìn)入料桶的物料質(zhì)量，計(jì)算得到錐形閥的物料流量。流化床開有8排6列共48個(gè)測(cè)量孔，可通過(guò)U型測(cè)壓管及光纖來(lái)測(cè)量流化床內(nèi)部的壓力及顆粒體積分?jǐn)?shù)。

圖3 錐形閥流量測(cè)量實(shí)驗(yàn)臺(tái)示意Fig.3 Schematic of experimental bench with cone valve

在搭建的實(shí)驗(yàn)裝置中，錐形閥開度x、流化風(fēng)速u、初始床料高度h、閥體標(biāo)稱尺寸D會(huì)對(duì)閥流量造成影響，通過(guò)測(cè)量不同工況(x,u,h,D)下錐形閥流量G、床內(nèi)壓力和顆粒體積分?jǐn)?shù)分布，來(lái)對(duì)錐形閥的流量特性進(jìn)行實(shí)驗(yàn)性研究。實(shí)驗(yàn)中所使用錐形閥標(biāo)稱尺寸為20，25，30，35 mm，初始床料高度分別為0.34，0.40，0.50 m;流化風(fēng)速調(diào)節(jié)范圍為0.15～0.60 m/s;給料機(jī)流量與閥出口流量保持基本相等，以維持實(shí)驗(yàn)中流化床床存量的動(dòng)態(tài)平衡及錐形閥入口壓力的大致穩(wěn)定。

在實(shí)驗(yàn)前，對(duì)給料器進(jìn)行了標(biāo)定，獲得了給料器不同頻率下的給料速率，以便能根據(jù)錐形閥物料流出速率的大小調(diào)節(jié)給料器頻率，從而保持流化床內(nèi)部物料量平衡。從圖4所示的標(biāo)定結(jié)果容易看出給料流量與給料器頻率呈現(xiàn)出較好的線性關(guān)系。

圖4 給料器流量隨頻率變化曲線Fig.4 Relationship between the feeding rate and frequency

3 實(shí)驗(yàn)結(jié)果分析與討論

3.1 閥流量隨開度x變化關(guān)系

圖5為測(cè)量得到的錐形閥流量G隨開度x的變化情況及擬合曲線。在閥開度為x∈(4 mm,15 mm)時(shí)，流量隨開度增加呈現(xiàn)出了很好地線性關(guān)系，當(dāng)x增大到x=xm后，繼續(xù)增大閥開度x，發(fā)現(xiàn)流量G基本處于穩(wěn)定狀態(tài)，此時(shí)閥流量達(dá)到當(dāng)前工況下的最大值Gm，定義流量最大值Gm為飽和流量，定義Gm對(duì)應(yīng)的最小開度xm為滿開度。

圖5 錐形閥流量隨開度變化關(guān)系Fig.5 Relationship between the flow rate and opening

結(jié)合文獻(xiàn)調(diào)研中顆粒孔口系統(tǒng)的相關(guān)研究，對(duì)錐形閥內(nèi)部流動(dòng)進(jìn)行分析，發(fā)現(xiàn)閥流量由床內(nèi)物料流動(dòng)情況和閥內(nèi)流通面積決定，而這2者又可表征為錐形閥進(jìn)料能力和錐形閥排料能力，2者中的較小值決定了閥的排料能力，也就是閥流量。從式(7)可以看出開度影響錐形閥內(nèi)部通流面積，開度增大時(shí)，閥排料能力隨之增加，但錐形閥進(jìn)料能力并不隨開度變化而變化。因此在增大開度的過(guò)程中，會(huì)存在滿開度，使得閥排料能力等于閥進(jìn)料能力，繼續(xù)增大開度，閥通流面積和排料能力增大，但閥進(jìn)料能力沒有增加，閥流量并不增加。

從對(duì)顆?？卓谙到y(tǒng)的文獻(xiàn)調(diào)研中可以看到，閥流量與孔口流通面積大致呈線性關(guān)系，從本文給出的閥通流面積S與開度x關(guān)系式(7),(8)中，可以看到S與x之間近似呈線性關(guān)系，這與實(shí)驗(yàn)得到流量與開度之間呈現(xiàn)線性關(guān)系相吻合。

3.2 閥流量隨流化風(fēng)速u變化關(guān)系

由圖6可知，隨著流化風(fēng)速的增加，錐形閥流量大致呈增加的趨勢(shì)，但增長(zhǎng)較慢。本實(shí)驗(yàn)使用的物料為80目(80～400 μm)粒徑的玻璃微珠，其對(duì)應(yīng)的臨界流化風(fēng)速umf為0.08 m/s，在u=0.2 m/s=2.5umf時(shí)，平均閥流量為279 g/s，而u=0.4 m/s=5umf時(shí)，平均閥流量為285 g/s，僅增長(zhǎng)了約2%，因此認(rèn)為當(dāng)流化風(fēng)速處于大于臨界流化風(fēng)速的一定范圍內(nèi)，其對(duì)閥流量影響較小。此外測(cè)量了不同流化風(fēng)速下的床壓，發(fā)現(xiàn)流化風(fēng)速改變對(duì)床壓幾乎沒有影響。

圖6 閥流量隨流化風(fēng)速變化關(guān)系Fig.6 Relationship between the flow rate and fluidized air velocity

3.3 閥流量隨床層高度h變化關(guān)系

由圖7可知，隨著初始床料高度的增加，閥流量有一個(gè)小幅度的增加，對(duì)于閥體標(biāo)稱尺寸D=20 mm，流化風(fēng)速u=0.2 m/s下的工況，h=34 mm對(duì)應(yīng)的飽和流量Gm=305.36 g/s，h=40 cm對(duì)應(yīng)的Gm=331.3 g/s，h=50 cm對(duì)應(yīng)的Gm=354.94 g/s;不同初始床料高度h對(duì)應(yīng)的滿開度xm也不同，隨著h的增大，xm有小幅度增大，但增長(zhǎng)幅度小于閥流量。初始床料高度對(duì)于閥流量的影響主要體現(xiàn)為對(duì)床壓的影響。初始床料高度增加，床壓變大，使得閥流量增大。

圖7 閥流量隨初始床料高度變化關(guān)系Fig.7 Relationship between the flow rate and initial height of bed material

3.4 閥流量隨閥體標(biāo)稱尺寸D變化關(guān)系

圖8給出了不同閥體標(biāo)稱尺寸下，錐形閥流量與開度的關(guān)系，可以看出不同標(biāo)稱尺寸下，閥流量隨開度變化規(guī)律一致。對(duì)圖8數(shù)據(jù)進(jìn)行分析處理，可以得到不同閥體標(biāo)稱尺寸工況對(duì)應(yīng)的飽和流量和飽和開度，如圖9所示，可以看到隨著標(biāo)稱尺寸增加，閥流量與對(duì)應(yīng)的飽和流量、滿開度也隨之增大。閥體標(biāo)稱尺寸D影響閥入口截面面積和閥內(nèi)流通區(qū)域，進(jìn)而影響閥進(jìn)料能力和閥出料能力，這也解釋了圖9中閥體標(biāo)稱尺寸D增大時(shí)，閥流量增大，同時(shí)飽和流量也隨之增大。

圖8 不同閥體標(biāo)稱尺寸下閥流量與開度關(guān)系Fig.8 Relationship between the flow rate and opening under different nominal size of valve

圖9 飽和流量與滿開度隨閥體標(biāo)稱尺寸變化關(guān)系Fig.9 Relationship between saturation flow,full-opening and the nominal size of valve

根據(jù)孔口系統(tǒng)物料流量隨通流半徑或通流面積之間的關(guān)系，易知閥流量與標(biāo)稱尺寸之間也應(yīng)呈與孔口系統(tǒng)類似的冪次關(guān)系，同時(shí)D=0時(shí)，應(yīng)有Gm=0。對(duì)實(shí)驗(yàn)數(shù)據(jù)進(jìn)行擬合，得到在本文實(shí)驗(yàn)范圍內(nèi)，閥體標(biāo)稱尺寸D與飽和流量Gm、滿開度xm之間的關(guān)系為

(9)

在顆?？卓谙到y(tǒng)流動(dòng)的文獻(xiàn)調(diào)研中，提到了流量G與尺寸D之間呈2.0～2.5次方關(guān)系，而從本文擬合得到的式(9)，可以看出Gm與D之間呈現(xiàn)2.133 2次方關(guān)系，與前人的相關(guān)研究結(jié)果吻合較好[19-25]。

4 結(jié) 論

(1)隨開度x的增加，閥流量先呈現(xiàn)出一小段弧形上升，之后線性上升到飽和流量點(diǎn)[xm,Gm]，繼續(xù)增大開度，錐形閥流量保持Gm不變。

(2)在流化風(fēng)速u>2umf的一定風(fēng)速范圍內(nèi)進(jìn)行實(shí)驗(yàn)，發(fā)現(xiàn)流化風(fēng)速u對(duì)閥流量影響較小，流化風(fēng)速?gòu)?.5umf增加到5.0umf，流量?jī)H增加2%，同時(shí)發(fā)現(xiàn)在不同流化風(fēng)速下，流化床內(nèi)部壓力分布幾乎不發(fā)生改變。

(3)改變初始床料高度h進(jìn)行實(shí)驗(yàn)，發(fā)現(xiàn)初始床料高度h對(duì)床壓有正向作用，隨著h的增大，飽和流量Gm及對(duì)應(yīng)的滿開度xm也有小幅度增加，同時(shí)Gm的增加幅度大于xm。

(4)不同閥體標(biāo)稱尺寸下，錐形閥的流量變化規(guī)律相似，利用前人擬合的顆粒孔口系統(tǒng)流量經(jīng)驗(yàn)公式形式來(lái)對(duì)本文的實(shí)驗(yàn)結(jié)果進(jìn)行擬合，發(fā)現(xiàn)飽和流量Gm與閥體標(biāo)稱尺寸s之間大致呈2.1次方關(guān)系，與文獻(xiàn)調(diào)研中前人的相關(guān)研究結(jié)果大致吻合。