基于EDF與IRS聯合研究的ORCHIDEE 2.0，為燃煤電廠飛灰排放控制提供可視化、可量化的決策工具

關鍵詞
electrostatic precipitator,fly ash emission,SO3 injection,back-corona,ORCHIDEE,靜電除塵,煙氣治理

在燃煤電廠中，靜電除塵器（ESP）依然是控制飛灰排放的主力裝備。為了滿足日益嚴格的排放標準，ESP 通常要長期運行在 99.5% 以上的收塵效率，這意味著從煤場到煙囪的每一個環節——煤質、燃燒、煙氣流場、電氣運行、振打維護乃至 SO3 調質——都可能對排放水平產生放大效應。如何讓運行人員在控制室中“看見”這些復雜耦合效應，并在幾分鐘內完成“如果……會怎樣”的工況推演，正是 ORCHIDEE 軟件試圖解決的問題。

這項研究由法國電力集團 EDF 的研發中心與運營電廠（Cordemais、Le Havre）以及意大利 IRS 公司聯合完成[8][11]。ORCHIDEE（法文全稱意為“燃煤電廠靜電除塵效率優化”）并非傳統意義上的廠商經驗公式軟件，而是立足于物理模擬的工業化工具：在不依賴制造商經驗數據庫的前提下，利用中等復雜度的物理建模對靜電除塵過程進行快速求解，為運行、檢修和技術改造提供定量支撐。

靜電除塵是一個高度耦合的多物理場過程。氣體流動、電場分布、離子產生與遷移、顆粒充電與遷移、積灰與振打排灰、飛灰電阻率與反電暈（back-corona）等環節相互影響。工業界常用的 Deutsch 方程雖然簡單易用，卻將這些細節全部“平均化”；而高校與研究機構開發的全物理模型，又往往計算量龐大、難以用于值班室的快速決策。ORCHIDEE 的策略，是在這兩種極端之間找到工程上可接受的“中間層”：保留關鍵物理機制，用參數化與插值算法替代全細節求解，把復雜模型“預計算”，將現場查詢環節壓縮到幾分鐘內即可給出結果。

其核心方法分為兩步：首先，基于 IRS 的 SPES 物理模型[2–7]，對代表 ESP 運行狀態的幾個主導參數（如煙氣流速、施加電壓、顆粒粒徑分布）進行系統的變分分析，覆蓋機組典型及極端工況；然后，通過多維非線性插值算法，把這些預計算結果封裝成標準 PC 上可快速調用的性能“地圖”，運行人員只需輸入當下煤質、負荷、氧量、煙溫、電場運行狀態等參數，即可即時得到分室、分電場的收塵效率和總排放濃度。對中心區電場，插值預測與全模型的偏差通常小于 10%，邊界區域也控制在 20% 以內，滿足工程診斷精度要求。

在最初版本基礎上，ORCHIDEE 2.0 進一步向實際工況靠攏，新增了預除塵模塊、SO3 注入調質模塊、積灰層動態模塊、反電暈模塊以及故障模擬模塊等，使其從“單純算效率”的工具，演進為可模擬完整煙氣治理鏈條的工業決策平臺。

在預除塵環節，ORCHIDEE 2.0 引入了兩種模型，兼顧理論性與實用性。一類是基于 Zenz 理論的旋風除塵模型，通過停留時間與切割粒徑的相關關系，將標準幾何條件下的旋風器效率與粒徑分布關聯起來，適用于滿足嚴格設計規則的規范旋風器；另一類是“固定捕集率”模型，允許用戶直接按粒徑段輸入來自制造商手冊或現場測試的效率，用以覆蓋大量非標準結構的預除塵設備。這一模塊讓運行人員可以量化“加裝/改造預除塵”對后端 ESP 效率與排放濃度的影響，為整體煙氣治理方案優化提供數據基礎。

飛灰電阻率是決定 ESP 成敗的關鍵參數，低硫煤、超潔凈煤和高堿金屬煤在工業實踐中屢屢因為反電暈或電場失效而導致排放超標。針對這一痛點，ORCHIDEE 2.0 采用并改進了 Bickelhaupt 提出的電阻率預測模型[12][13]，在原始相關式基礎上處理了插值區間的非連續性問題，擴大了適用溫度與成分范圍，并將在線或設計 SO3 注入量直接并入燃燒生成 SO3 的總濃度中統一計算。在 130 ℃ 工況下，模型模擬的電阻率–SO3 濃度關系與 Bickelhaupt 報告中的結果高度吻合，趨勢和數值誤差均處于工程可接受范圍[13]。

更具工程價值的是將 SO3 調質的效果貫穿到整體 ESP 模擬中。以法國 Le Havre 4 號機組為例，該 600 MW 機組配套兩箱四電場 Lurgi ESP，原煤為南非產低硫煤，未經 SO3 調質時飛灰電阻率約 2×10^11 Ω·cm，若直接入 ESP，排放遠超 50 mg/Nm3 法規限值。運行方在入口 130 ℃ 煙溫下注入約 15 ppm SO3 后，ORCHIDEE 預測煙囪排放約 28 mg/Nm3（干基 6%O2），而現場四小時連續測試的平均值為 32 mg/Nm3，誤差處于儀表精度與工況波動范圍內[11]。類似地，在英國 Cottam 1 號 500 MW 機組上，南非煤未經調質時模擬排放約 260 mg/Nm3，通過約 11.6 ppm SO3 注入，模擬值降至 36 mg/Nm3，而實測結果為 38 mg/Nm3，驗證了 ORCHIDEE 在不同鍋爐、不同 ESP 結構和不同煤質條件下的適用性。

傳統 ESP 模型多止步于“顆粒進入電場前”的過程，而真正困擾運行人員的，往往是已經沉積在極板上的積灰層行為：何時開始滑落？何時被風蝕再飛揚？振打太快還是太慢對排放究竟有多大影響？針對這一“黑箱”環節，EDF 與 IRS 借助歐盟 ABRICOS 工業試驗項目[14]，在工業級試驗 ESP 上開展了系統試驗與質量平衡測量，從而首次獲得了可量化的灰層生長、滑移、再飛揚與振打瞬態排放數據。在此基礎上，ORCHIDEE 2.0 構建了灰層動力學模型：一方面，通過顆粒沉積通量、風蝕與離子風導致的再飛揚通量，以及沿板面下滑通量的連續質量守恒，計算灰層厚度的緩慢演化；另一方面，在振打瞬間引入“塊體脫落–部分直接落入灰斗–部分被再捕集–部分重返氣相”的瞬態質量轉移模型[15]，從而把振打制度（頻率、順序、停歇時間）與排放濃度建立起定量聯系。

更復雜的情況是高電阻灰條件下的反電暈。ORCHIDEE 2.0 在模型中顯式引入“灰層電壓降”與“反電暈電流”兩大效應[9][10][16]：一方面，沉積層自身的高電阻會在其內部形成明顯電壓降，壓縮有效放電間隙電壓，導致正常負輝光放電減弱；另一方面，當灰層電壓升高到一定閾值時，會在層內微孔中觸發正電暈放電，產生向陽極遷移的正離子流。雖然這一正離子流與主負離子流在電流上疊加，但兩者的空間電荷部分抵消，顯著降低了顆粒的凈充電率和遷移速度，整體收塵效率下降。ORCHIDEE 通過等效電路形式對輝光電流與反電暈電流–電壓關系進行非線性求解，再反推對顆粒充電與效率的影響，使工程人員不再停留于“看到電流上升、效率反降”的經驗判斷，而是可以在仿真中提前識別“反電暈窗口”，用配煤或 SO3 調質將工況保持在安全區間。

這些模塊之間并非簡單串聯，而是通過一個非線性主循環進行耦合求解：預除塵影響入口粉塵濃度與粒徑分布，進而改變 ESP 內沉積與再飛揚通量；灰層厚度又影響電阻率有效路徑與電壓分配，反過來改變放電電流與反電暈強度；SO3 注入和煤質變化改變本征電阻率，對上述過程同時施加“基線”修正；電氣與機械故障（如極線結灰、極板變形、絕緣子泄漏）則通過改變局部電場與流場條件，引發全局效率變化。ORCHIDEE 2.0 將這些反饋路徑納入統一迭代框架中，直至粉塵濃度與電氣參數收斂，從而給出穩定的收塵效率和分場 V–I 曲線。

在實際電廠應用中，這種多模塊耦合的價值體現在兩個層面。一是對運行調整的指導：例如在 Cottam 1 號機組上，通過將振打電機停歇時間加倍（慢振打）或減半（快振打）的現場試驗顯示，偏離“標準振打”兩側都會造成排放上升；ORCHIDEE 模擬值與實測數據的趨勢和量級高度一致，說明模型能夠抓住灰層動力學與排放之間的本質聯系。二是對故障診斷的輔助：在 Cordemais 5 號機組上，研究人員現場采集了各電場分段的電壓–電流曲線，并在檢修停機后進行內部檢查。結果顯示，某一分段存在嚴重極線結渣，其 V–I 曲線整體“右移下沉”，對應電流抑制特征。ORCHIDEE 的故障庫中引入“極線/極板結灰”“極線–極板偏心”“絕緣子表面泄漏”等故障模型，通過擬合實際 V–I 曲線，能夠反推最可能的故障類型與嚴重程度，為維護人員提供更精確的檢修優先級依據。

從行業視角來看，ORCHIDEE 的意義在于把原本散落在科研論文中的物理模型、復雜的數值求解與電廠一線的實際決策場景打通。對運行部門而言，它是一個可以快速回答“如果提高負荷、如果更換煤源、如果某電場退出、如果調整 SO3 注入量、如果修改振打周期”會對排放產生何種影響的推演平臺，有助于在不降負荷、不大修設備的前提下，通過過程優化維持甚至提高 ESP 收塵效率。對檢修與技術部門而言，它提供了一種“先在仿真中試錯”的工具，可以在實施前預估更換內部構件、改造電源、增加電場、調整煙道布置等方案的減排效果與性價比，避免投入“無感知收益”的改造項目。

隨著我國燃煤機組向靈活性調峰、超低排放與高比例低硫煤摻燒方向發展，類似 ORCHIDEE 這樣兼具物理真實性和工程易用性的 ESP 仿真工具，將愈發成為電廠煙氣治理和運維團隊的標配“第二雙眼睛”。在更加嚴格的排放約束與更加復雜的工況波動面前，誰能更早把物理機理轉化為日常運行可用的“按鈕”和“曲線”，誰就能在環保達標與經濟性之間掌握更大的主動權。

參考文獻
[1] Salvi J. A methodology to simplify the ESP models[C]//Proceedings of VIII ICESP. Birmingham, USA, 2001.
[2] Gallimberti I. The simulation of corona discharges under practical precipitator conditions[C]//Symposium on the Transfer and Utilisation of Particulate Control Technology. EPRI, New Orleans, Vol.1, 1986.
[3] Bellagamba B, Mattachini F, Gallimberti I, Turri R, Gazzani A, Tromboni U. A mathematical model for simulation of large scale electrostatic precipitators[C]//5th International Conference on Electrostatic Precipitation. Washington D.C., 1993.
[4] Gallimberti I, Gazzani A, Tromboni U, Lami E, Mattachini F, Trebbi G. Physical simulation of the particle migration in ESP, Part I – Model description[C]//VI ICESP. Budapest, 1996.
[5] Gallimberti I, Gazzani A, Tromboni U, Lami E, Mattachini F, Trebbi G. Physical simulation of the particle migration in ESP, Part II – Application results[C]//VI ICESP. Budapest, 1996.
[6] Lami E, Mattachini F, Gallimberti I, Turri R, Tromboni U. A numerical procedure for computing the voltage–current characteristics in ESP configuration[J]. Journal of Electrostatics, 1995, 34: 385–399.
[7] Gallimberti I. Recent advancements in the physical modelling of electrostatic precipitators[J]. Journal of Electrostatics, 1998, 43: 219–247.
[8] Arrondel V, Bacchiega G, Gallimberti I. ESP modelling: from university to industrial application[C]//VIII ICESP. Birmingham, USA, 2001.
[9] Gallimberti I. Detailed mass balance in electrostatic precipitators under industrial operating conditions[C]//IX ICESP. Kruger, South Africa, 2004.
[10] Bacchiega G, Gallimberti I, Arrondel V, Raizer Ph, Lecointre J, Hamlil M. Static and dynamic back-corona characteristics[C]//IX ICESP. Kruger, South Africa, 2004.
[11] Arrondel V, Salvi J, Gallimberti I, Bacchiega G. ORCHIDEE: efficiency optimisation of coal ash collection in electrostatic precipitators[C]//IX ICESP. Kruger, South Africa, 2004.
[12] Bickelhaupt R. A technique for predicting fly ash resistivity[R]. EPA-600/7-79-204, 1979.
[13] Bickelhaupt R. Fly ash resistivity prediction improvement with emphasis on sulphur trioxide[R]. EPA-600/7-86/010, 1986.
[14] Bacchiega G, Gallimberti I, Sani E, Sala R, Arrondel V, Hamlil M, Christensen E. Experimental study of the mass balance in a pilot industrial ESP[J]. Journal of Electrostatics, 2005.
[15] Arrondel V, Caraman N, Hamlil M, Bacchiega G, Gallimberti I. Development of an industrial model of rapping – effect on the collecting efficiency[C]//X ICESP. Cairns, Australia, 2006.
[16] Bacchiega G, Gallimberti I, Arrondel V, Caraman N, Hamlil M. Back-corona model for prediction of ESP efficiency and voltage-current characteristics[C]//X ICESP. Cairns, Australia, 2006.

獲取更多靜電除塵相關專業論文，請訪問 https://isesp.org/conference-papers/