董鴻霖，朱小良

(東南大學(xué) 能源與環(huán)境學(xué)院，南京 210096)

我國大氣污染形勢較為嚴(yán)峻，電力行業(yè)燃燒煤炭產(chǎn)生的SO2、NOx、顆粒物(PM)等為主要的污染來源之一，為貫徹落實《中華人民共和國環(huán)境保護法》、《大氣污染防治行動計劃》等法律法規(guī)，環(huán)境保護部發(fā)布了一系列的標(biāo)準(zhǔn)和規(guī)范，使得燃煤電廠必須準(zhǔn)確控制污染物的排放。對氣體污染物排放量的準(zhǔn)確測量是以對煙氣成分、溫度、壓力、濕度、流量等的準(zhǔn)確測量為基礎(chǔ)的[1]。就目前的技術(shù)手段而言，利用煙氣連續(xù)排放監(jiān)測系統(tǒng)(CEMS)是掌握燃煤電廠氣體污染物排放狀況的重要手段，可通過對煙氣參數(shù)的測量計算得到氣體污染物的排放量[2]。CEMS由一系列子系統(tǒng)組成，其中，準(zhǔn)確測量流量是CEMS有效工作的基礎(chǔ)。CEMS通過點測量或線測量的方式，采用壓差法、熱平衡法或超聲波法得到煙氣截面的平均流速，這就要求所選取的測量截面必須具有代表性。因此，根據(jù)GB/T 16157—1996 《固定污染源排氣中顆粒物測定與氣態(tài)污染物采樣方法》的要求，CEMS所選取的截面流速需要相對平穩(wěn)，應(yīng)避開彎頭、斷面等急劇變化的部位；上游長度應(yīng)大于或等于4倍煙道直徑、下游長度應(yīng)大于或等于2倍煙道直徑。然而現(xiàn)實中火電廠在脫硫后的引風(fēng)機至煙囪進口部分的矩形尾部煙道存在90°的彎道且直管段長度僅為煙道直徑的1～3倍。這就導(dǎo)致煙氣在尾部煙道不能充分發(fā)展，煙道內(nèi)流場情況較為復(fù)雜，存在高、中、低速區(qū)，以及回流區(qū)、二次流等。CEMS大多采用基于速度-面積法的截面插入多點式風(fēng)量測量裝置[3]，其測點布置方式存在測點布置過多的問題，對于大型煙道至少需要30個以上，且對于帶有90°彎的尾部煙道，有部分測點布置于流速波動較大的回流區(qū)導(dǎo)致測量數(shù)據(jù)精度較低。

筆者基于流體力學(xué)中的自動模型化(簡稱自?；?理論，在充分分析煙道內(nèi)流動狀況的基礎(chǔ)上找到最優(yōu)的流速測點布置方案，為燃煤電廠的污染物排放控制打下基礎(chǔ)。

1 煙道內(nèi)流動狀況理論分析

尾部煙道內(nèi)的煙氣為黏性流體，其流動過程符合流體力學(xué)中的相關(guān)概念。以往帶有90°彎的矩形尾部煙道測點布置大多直接采用速度-面積法中對數(shù)-線性法或?qū)?shù)-切比雪夫法的測點布置方案，沒有充分分析煙道內(nèi)煙氣的真實流動情況，測點部分布置于流動狀況復(fù)雜的區(qū)域?qū)е略摬糠譁y點數(shù)據(jù)波動較大，甚至出現(xiàn)讀數(shù)為零的情況，影響流量的準(zhǔn)確測量。因此，有必要對尾部煙道的流動情況進行分析并指導(dǎo)測點的布置。

1.1 自模化

自?；噶黧w雷諾數(shù)在某區(qū)間范圍內(nèi)變化時失去影響作用或作用減小的一種現(xiàn)象，該區(qū)間稱為自?；瘏^(qū)[4]。黏性流體存在兩個自?；瘏^(qū)，在層流范圍內(nèi)流體的流速場彼此相似，層流范圍內(nèi)的雷諾數(shù)范圍稱為第一自?；瘏^(qū)。當(dāng)雷諾數(shù)過大且超過一定界限時，黏性力影響相對減弱，此時即使提高雷諾數(shù)，流動現(xiàn)象和流動性能也不會有太大變化；雖然雷諾數(shù)不同但有相同的黏性效果[5]，該雷諾數(shù)區(qū)間稱為第二自?；瘏^(qū)。對于大型燃煤電廠，其脫硫后至煙囪進口部分的尾部煙道流速一般在5～25 m/s，對于矩形煙道，異形斷面管道特征尺寸de為：

(1)

式中：A為管道截面積；S為管道潤濕周長。

雷諾數(shù)的表達式為：

(2)

式中：Re為雷諾數(shù)；u為流體平均速度；v為流體運動黏度。對于大型火電廠尾部煙道，不同工況下其煙氣雷諾數(shù)一般大于106。

流體進入第二自?；瘏^(qū)的雷諾數(shù)臨界值(Re臨)[5]見表1(d為管徑，mm；K當(dāng)為當(dāng)量粗糙度，蒸汽管道的K當(dāng)為0.1～0.5 mm，風(fēng)道的K當(dāng)為0.8 mm)。由表1可以看出：大型火電廠尾部煙道內(nèi)煙氣雷諾數(shù)遠大于臨界值，煙氣處在第二自?；瘏^(qū)，即理論上可以認(rèn)為煙道內(nèi)煙氣在此流速范圍內(nèi)形成的流場相似。

表1 進入第二自模化區(qū)的雷諾數(shù)臨界值

1.2 中值定理

根據(jù)自?；碚?，火電廠尾部煙道在不同的進口條件下，煙道內(nèi)流場相似。測量尾部煙道截面平均流速最理想的測量方法是找到一個或多個固定點，在不同的進口煙氣流速條件下該固定點測得的流速能夠代表整個截面上的平均流速。顯然，在煙道尾部任一截面存在高、中、低速區(qū)，速度介于0(近壁面、低速區(qū))至大于截面平均流速(高速區(qū))的最大值。

由于速度在空間和時間上的變化是連續(xù)可微的，根據(jù)中值定理，在某一定截面上一定存在一條或多條連續(xù)閉合曲線，該曲線上所有點的速度均為截面上的平均流速。根據(jù)自?；碚?，不同進口條件下同一截面上對應(yīng)的等速線分布應(yīng)該具有一定的規(guī)律。

1.3 二次流

二次流是由主流引起的伴隨流動，常見的二次流有Taylor渦流、Gortler渦流及Dean渦流[6]。Dean渦流是由于流體在有彎折的管道中流動，垂直于流體流動方向的壓力梯度不平衡，并在離心力的作用下，在流體的截面上下產(chǎn)生與旋向相反的漩渦。當(dāng)雷諾數(shù)較低時，漩渦上下分布對稱且穩(wěn)定，雷諾數(shù)過大則會破壞該對稱性。

對于帶有90°彎的矩形煙道，煙氣經(jīng)過彎折的過程中會產(chǎn)生Dean渦流并影響下游流動，由于Dean渦流的存在，煙氣流線會偏離來流方向，使用煙道常用的煙氣測速裝置會產(chǎn)生誤差。對于熱線風(fēng)速儀，此時測得的煙氣流速大小為煙氣的實際流速而非流速在來流方向上的分速度；對于畢托管等測量裝置，截面上二次流會使得測速裝置難以正對測點速度方向而導(dǎo)致測量誤差。故在選擇測量截面的時候，必須考慮截面上二次流的影響。

2 模型的建立及仿真

2.1 仿真模型

大型火電廠尾部煙道進口流速及其截面特征尺寸在大于一定條件下，其尾部煙道內(nèi)流體雷諾數(shù)一定位于第二自模化區(qū)，為探究一般性規(guī)律，建立和某600 MW火電廠尾部煙道具有相似尺寸及形狀的尾部煙道及煙囪的簡化模型(見圖1)。該尾部煙道長度為12 m，其中有4 m為彎折部分，從脫硫后引風(fēng)機到彎折部分的長度為8 m，煙道截面為4 m×4 m的矩形；煙囪高度60 m，煙道設(shè)置距離地面高8 m，該煙道存在90°彎且彎折后尾部煙道長度僅為8 m，顯然無法滿足CEMS截面選取的要求。

圖1 煙道-煙囪結(jié)構(gòu)示意圖

2.2 邊界條件

采用3D建模軟件完成對煙道及煙囪的幾何建模后，綜合考慮計算效率及計算精度，采用結(jié)構(gòu)化網(wǎng)格為主、非結(jié)構(gòu)化網(wǎng)格為輔的網(wǎng)格劃分方式，其中煙道部分主要為結(jié)構(gòu)化網(wǎng)格。利用計算流體力學(xué)(CFD)軟件進行求解計算，考慮到一般尾部煙道進口流速為5～25 m/s，馬赫數(shù)小于0.1，對應(yīng)可壓縮系數(shù)小于0.003，故假設(shè)其為不可壓縮流體。為驗證本文研究對象及研究范圍適用于自?；碚?，設(shè)置進口流體速度分別為4 m/s、10 m/s、13 m/s、19 m/s，湍流強度為5%，溫度為50 ℃，密度為1.16 kg/m3，動力黏度為1.72×10-5Pa·s，采用標(biāo)準(zhǔn)k-ε湍流模型。求解選用壓力求解器，標(biāo)準(zhǔn)SIMPLE算法，其中壓力、動量空間離散采用二階迎風(fēng)格式，求解殘差余項設(shè)置為10-5。

2.3 仿真結(jié)果分析

該測點布置方案主要位于尾部煙道，故重點分析尾部煙道內(nèi)流動情況。圖2為進口流速為10 m/s時尾部煙道內(nèi)的流線分布情況。

圖2 尾部煙道流線圖

將正對來流方向右側(cè)(左側(cè))的煙道壁面稱為右側(cè)(左側(cè))壁面，由圖2可以看出：煙氣以10 m/s的速度從尾部煙道進口均勻進入尾部煙道，在經(jīng)過彎折部分后，煙道內(nèi)流速分布不再均勻；靠近煙道左側(cè)出現(xiàn)了速度遠大于進口設(shè)定流速的高速區(qū)，經(jīng)過彎折之前靠近煙道右側(cè)部分的流體在慣性力和離心力的共同作用下，不再沿著右側(cè)壁面流動，導(dǎo)致煙氣流過90°彎后出現(xiàn)了流動混亂的低速區(qū)及回流區(qū)；由于速度是溫度和壓力的連續(xù)函數(shù)，故在低速區(qū)和高速區(qū)中間部分為流速接近進口設(shè)定流速的中速區(qū)。

在該模型中，上游與下游的距離比為2∶1時對應(yīng)的煙道截面距離煙囪進口2.67 m。將4 m/s、10 m/s、13 m/s、19 m/s進口流速條件下距離煙囪6 m、4 m、2.67 m、1.5 m的截面上速度分布云圖進行比較，結(jié)果見圖3～圖6。

圖3 距離煙囪6 m截面速度云圖

圖4 距離煙囪4 m截面速度云圖

圖5 距離煙囪2.67 m截面速度云圖

圖6 距離煙囪1.5 m截面速度云圖

由圖3～圖6可以看出：在不同進口流速條件下，距離煙囪不同位置截面上的速度云圖分布相似，均存在高、中、低速區(qū)及回流區(qū)，不同速度區(qū)域隨著流體的發(fā)展呈現(xiàn)高速區(qū)范圍逐漸減小，中速區(qū)及低速區(qū)范圍逐漸增大的趨勢；中速區(qū)范圍增大意味著隨著流體的發(fā)展，中速區(qū)截面上速度接近進口流速設(shè)定值的空間范圍在增大，速度分布梯度減小，在該截面的中速區(qū)上測量的誤差就會降低。在4種進口流速條件下流體均處于自?；瘏^(qū)，流場內(nèi)流動規(guī)律相似，因此自?；碚撨m用于尾部煙道在進口流速大于5 m/s的情況，即在5～25 m/s進口流速的條件下，煙道內(nèi)部流場分布相似，對于一定范圍內(nèi)不同進口流速條件下分析測量截面的位置，可以在單一進口流速條件下重點進行分析。

3 測點布置

3.1 測點布置截面的選取

選取10 m/s的進口流速條件，提取距離煙囪進口6 m、5 m、4 m、2.67 m、2 m、1.5 m、1.3 m、1.2 m、1.1 m、1 m、0.9 m、0.8 m截面上的流線并進行比對分析，結(jié)果見圖7。

圖7 距離煙囪進口不同距離截面流線對比圖

由圖7可以看出：隨著尾部煙道流動的發(fā)展，截面流線分布也在不斷變化。在距離煙囪進口5～6 m，截面上高速區(qū)和中速區(qū)流線橫向分布較為均勻，意味著該范圍截面上流速存在橫向分速度。隨著流動的發(fā)展，在距離煙囪進口1.2～4 m都可以觀察到存在較為混亂的二次流，尤其是距離煙囪進口2.67 m的截面上，流線分布極其復(fù)雜，不僅回流區(qū)存在二次流，位于高速區(qū)的上下兩側(cè)同樣存在二次流，當(dāng)距離煙囪進口小于1.5 m時，高速區(qū)上下兩側(cè)的二次流逐漸消失，當(dāng)距離煙囪進口約1.2 m時低速回流區(qū)內(nèi)的二次流也逐漸消失，考慮到1.1 m時截面上的二次流剛消失，而0.9 m時位于中速區(qū)上側(cè)的區(qū)域偏離主流方向范圍略大，在實際的流速測量中會產(chǎn)生較大偏差，對比之下選擇距離煙囪進口1 m的截面進行測點布置。

3.2 普適性測點的布置研究

針對已經(jīng)選定的煙道截面，其上測點布置必須具有普適性，即在不同進口煙氣流動情況下都有很好的精確度，該布置方式在不同邊界條件下所得到的流速在數(shù)值上應(yīng)靠近所設(shè)定的進口流速。

將仿真后距離煙囪進口1 m截面的流線圖進行對比，結(jié)果見圖8。

圖8 不同進口流速下距離煙囪1 m截面流線分布

由圖8可以看出：不同進口流速下距離煙囪進口1 m處截面上的流線分布變化不大，高速區(qū)、中速區(qū)、低速區(qū)及回流區(qū)的分布情況不隨進口流速的變化而變化且分布相對均勻，故可以認(rèn)定在該平面上一定存在可以代表進口流速的最優(yōu)測點，截面中部區(qū)域為中速區(qū)且其內(nèi)包含有進口流速設(shè)定值；中部區(qū)域的流線相當(dāng)密集且整個平面的流線由中部區(qū)域呈向外發(fā)散狀，說明空間內(nèi)此區(qū)域的主流方向受到其他如二次流等的干擾很小，偏離主流方向的分速度很小，所以不同進口流速條件下該測點一定位于流線圖中中部區(qū)域內(nèi)。將截面上進口流速設(shè)定值的速度等值線提取繪制于同一坐標(biāo)軸中，結(jié)果見圖9。

圖9 不同流速下等值線

由圖9可以看出：不同進口流速條件下距離煙囪進口1 m截面上對應(yīng)等速線分布情況相似，為了更好地判斷普適性測點的選取位置，重點觀察高度為0.5～3.5 m，距離煙道左側(cè)2.0～2.2 m內(nèi)煙氣等速線分布情況(見圖10)。

圖10 高度為0.5～3.5 m，距離煙道左側(cè)2.0～2.2 m內(nèi)煙氣等速線分布

由圖10可以看出：隨著進口流速設(shè)定值的改變，等速線在高度0.9～1.3 m、對應(yīng)距離煙道左側(cè)2.04～2.07 m；高度3.0～3.3 m、對應(yīng)距離煙道左側(cè)2.02～2.04 m的等速線分布較為密集，即在該范圍內(nèi)的任意一條等速線上的點所仿真出的速度都和初始設(shè)定值極為接近。于是重點觀察高度為0.9～1.3 m、3.0～3.3 m的等速線分布情況(見圖11、圖12)。

圖11 高度為0.9～1.3 m內(nèi)煙氣等速線分布

圖12 高度為3.0～3.3 m內(nèi)煙氣等速線分布

由圖11、圖12可以看出：高度1.13 m、距離煙道左側(cè)2.053 m(記為A點)，高度3.16 m、距離煙道左側(cè)2.024 m(記為B點)處的兩個點上等速線更加密集，故考慮選擇這兩個點為該截面上的普適性測點，其在截面上的具體位置見圖13。

圖13 最優(yōu)測點位置

由圖13可以看出：A點位置較為理想，附近的流線呈密集的發(fā)散狀說明該點附近流線沿水平方向發(fā)展且來流方向以外的流動干擾很小，使用設(shè)備(如畢托管等)進行流速測量時可以準(zhǔn)確對準(zhǔn)來流方向，具有較高的精確度；而B點位于煙道的上側(cè)，附近流場呈現(xiàn)較為均勻的向煙道左側(cè)流動的流線，說明在該點上除主流方向的流動外，煙氣還有自中速區(qū)向高速區(qū)流動的分速度。

3.3 測點精確度分析

4 m/s、10 m/s、13 m/s、19 m/s進口條件下對應(yīng)A、B點的絕對速度VA、VB及垂直于截面方向分速度VAw、VBw見表2。未參與測點分析的進口流速為7 m/s、16 m/s、22 m/s，此時A、B點流速數(shù)據(jù)見表3。

表2 進口流速設(shè)定值下A、B兩點流速 m/s

表3 未分析的進口流速設(shè)定值下A、B兩點流速 m/s

以進口流速設(shè)定值為橫坐標(biāo)，相對誤差為縱坐標(biāo)，將VA與VAw的相對誤差、VB與VBw相對誤差、VA和VB均值與進口流速設(shè)定值相對誤差繪于同一坐標(biāo)系(見圖14)。由表2、表3及圖14可以看出：不同流速下選定的A、B測點都能很好地代表進口流速設(shè)定值，其中A點流線相較于B點流線更好地垂直于截面，而VB與VBw的相對誤差小于0.25%，可見B點流線同樣近似垂直于選定截面，VA和VB的均值與進口流速設(shè)定值之間的相對誤差全部小于0.3%，可見在充分考慮截面選取、速度場分布及截面流線分布情況下，選取的A、B點在不同進口流速條件下其測量結(jié)果皆具有極高的精確度。

圖14 不同進口流速下測點相對誤差

4 結(jié)語

通過對大型燃煤電廠尾部煙道進行理論分析，選取具有代表性的尾部煙道結(jié)構(gòu)進行內(nèi)部流場研究，最終選定合適的速度測量截面和測點布置方案。具體結(jié)論如下：

(1)對于燃煤電廠大尺寸尾部煙道，其在不同工況下內(nèi)部流體雷諾數(shù)位于第二自?；瘏^(qū)，其流場分布滿足自?；煌r下當(dāng)內(nèi)部流場趨于穩(wěn)定，即工況保持穩(wěn)定時，其內(nèi)部流場相似。

(2)煙氣在經(jīng)過彎折后，其截面上會產(chǎn)生二次流使得流體速度偏離測量截面垂直方向，截面上二次流現(xiàn)象隨著流動的發(fā)展會逐漸消失?，F(xiàn)實中對于測量截面的選取應(yīng)充分考慮二次流的影響。

(3)充分考慮自模化理論，通過將多個等速線繪制于同一坐標(biāo)系中的辦法找到等速線分布最為密集的兩個點，其在不同工況下均有很好的適用性。根據(jù)仿真結(jié)果，進口流速在4～22 m/s時，兩個測點的速度均值與設(shè)定值之間的相對誤差均小于0.3%，且流速方向與測量截面夾角近乎垂直，在實際測量過程中可以極大地減小由于測速儀器布置時無法正對來流方向帶來的部分測量誤差。