姜周曙， 翁翔彬， 王 劍， 雷淳正

(杭州電子科技大學自動化學院，浙江杭州310018)

淡水短缺已成為阻礙社會經(jīng)濟發(fā)展的重要瓶頸，海水淡化技術(shù)是解決淡水資源短缺的主要途徑之一。經(jīng)過50多年的發(fā)展，從20世紀90年代開始，以能耗低、建造周期短、設(shè)備安裝簡便等諸多優(yōu)點集一身的反滲透海水淡化技術(shù)被公認為是21世紀海水淡化的主導技術(shù)［1］，目前在國內(nèi)外已得到廣泛應用。

隨著反滲透技術(shù)的大力發(fā)展，全球各地的反滲透海水淡化工程的規(guī)模也越來越大，日產(chǎn)幾十萬噸級的海水淡化工廠不斷出現(xiàn)。以色列阿什凱隆的哈代拉海水淡化廠日產(chǎn)淡水33萬t，為現(xiàn)在世界上運營中最大反滲透海水淡化廠。在我國，浙江舟山六橫島日產(chǎn)10萬t的海水淡化工程即將竣工;河北曹妃甸阿科凌海水淡化項目完成一期日產(chǎn)5萬t的工程，全部工程預期日產(chǎn)淡水100萬t，如若建成將成為世界上最大的反滲透工程。針對如此大規(guī)模復雜的海水淡化工程，對其進行故障分析是保證系統(tǒng)正常運行的一個重要前提。文中將從反滲透海水淡化工程的工藝流程入手，以“脫鹽率與產(chǎn)水量上升”故障為例，建立故障樹并進行定性定量分析，找出該故障的成因和防止其發(fā)生的方法。

1 大型反滲透海水淡化工程工藝流程分析

根據(jù)相關(guān)資料統(tǒng)計，海水淡化工程出現(xiàn)80%的故障都源于對海水原水的預處理不當，因此研究海水淡化工藝流程尤其是預處理階段的步驟極其重要。給水預處理對反滲透安全運行至關(guān)重要。膜材料具有一定的化學穩(wěn)定性，余氯、pH值、細菌、溫度、某些化學物質(zhì)等對穩(wěn)定性有很大影響，反滲透系統(tǒng)對給水水質(zhì)要求較高，通常要求SDI＜3.0(5.0)。因此反滲透系統(tǒng)之前通常配有較為完善的預處理系統(tǒng)，預處理的作用是為了降低污染指數(shù)，使系統(tǒng)經(jīng)過更長時間運行后再進行清洗;減少硬度和有機物的沉淀，確保膜免受機械和化學損傷，以使膜有良好的性能和足夠長的使用壽命。下面以浙江舟山六橫島日產(chǎn)10萬t反滲透海水淡化工程為例，介紹流程工藝以及預處理環(huán)節(jié)。

典型的大規(guī)模反滲透海水淡化工程包括取水、預處理、反滲透、后處理和供水等處理單元。設(shè)在海邊的深水泵將海水原水送入預處理單元。海水的一級預處理單元通過加入FeCl3、骨膠并使用隔板絮凝等處理環(huán)節(jié)，有效去除海水中的懸浮物、有機物、無機膠體等雜質(zhì)，在二級預處理單元中海水依次通過多介質(zhì)機械過濾器和保安過濾器以保證出水水質(zhì)達到反滲透水質(zhì)要求，兩級預處理后的海水經(jīng)高壓泵和能量回收裝置提升到反滲透所需壓力后，送入反滲透單元進行淡化，淡化后的產(chǎn)品水經(jīng)過后處理調(diào)質(zhì)后，供水廠調(diào)度。反滲透后的高壓濃海水經(jīng)能量回收后，變?yōu)榈蛪簼夂Ｋ?，?jīng)合理處置后排放。整個工程工藝流程如圖1所示。

圖1 反滲透海水淡化系統(tǒng)工藝流程Fig.1 Reverse osmosis desalination process

2 反滲透海水淡化故障樹

故障樹分析(Fault Tree Analysis，F(xiàn)TA)是一種將系統(tǒng)故障形成的原因由總體至部分按樹枝狀逐級細化的分析方法，因而是對復雜動態(tài)系統(tǒng)的設(shè)計、工廠試驗或現(xiàn)場發(fā)現(xiàn)失效形式進行可靠性分析的工具，其目的是判明故障，確定故障的原因、影響和發(fā)生概率［2］。

FTA最早由美國貝爾實驗室的沃森和莫恩斯提出，其后波音公司哈斯爾等人將其運用到了宇航領(lǐng)域，使飛機的設(shè)計有了重要的改進。而后美國原子能委員會利用FTA分析了核電站可能發(fā)生的各種事故的概率，并由此得出核能是非常安全的能源結(jié)論。隨著FTA方法的逐漸完整和日益廣泛的應用，已被公認為是可靠性、安全性分析的一種快速高效的方法［3］。

2.1 建立故障樹

故障樹模型是描述診斷對象結(jié)構(gòu)、功能和關(guān)系的一種定性因果模型［4］。以出現(xiàn)的故障為出發(fā)點，逐級找出導致發(fā)生故障的原因，并用相應的符號將故障現(xiàn)象和各級故障原因連接起來，形成一個因果關(guān)系圖。故障樹建造的完善程度將直接影響專家系統(tǒng)的診斷推理效率及正確性［5］。建樹過程是一種反復的過程，要在廣泛吸取系統(tǒng)設(shè)計、使用過程的知識經(jīng)驗并深入分析系統(tǒng)的實際情況的基礎(chǔ)上不斷修改完善。故障樹建造過程的實質(zhì)就是尋找出所研究的系統(tǒng)故障和導致系統(tǒng)故障的諸多因素之間的邏輯關(guān)系，并把這種關(guān)系用故障樹的圖形符號表示出來［6］，故障樹中常用事件與邏輯門如表1所示。

脫鹽率和產(chǎn)水量是反滲透海水淡化工程中兩個重要的指標，淡化水合格與否完全取決于這兩個指標是否達標。因此，一旦脫鹽率和產(chǎn)水量有異常變化，就要對整個系統(tǒng)進行排查，確定故障位置，采取相應的糾正措施［7-8］。文中選取“脫鹽率下降且產(chǎn)水量上升”為頂事件的例子，自上而下繪制導致該故障發(fā)生的故障樹并對其進行分析(見圖2)。

表1 故障樹事件符號與邏輯門符號釋義Tab.1 Explanation of logic and event signs

圖2 “脫鹽率下降且產(chǎn)水量上升”故障樹Fig.2 Fault tree of desalting rate reduce and water production rise

2.2 故障樹定性分析

為了找出導致頂事件發(fā)生的所有可能的故障模式，對故障樹進行定性分析［9］。一般的定性分析方法有求解故障樹的最小割集、最小徑集以及結(jié)構(gòu)重要度。

2.2.1 最小割集 在故障樹中，幾個底事件的集合同時發(fā)生時引起頂事件的發(fā)生，這個集合便稱為割集。而在這個割集中若任意去掉其中任意一個底事件，就不再是割集，即不能再引起頂事件的發(fā)生，則這個割集被稱之為最小割集(Minimal Cut Sets，MCS)。對于故障樹來說，它的全部最小割集的完整集合代表了頂事件發(fā)生的所有可能性。

2.2.2 最小徑集 在故障樹中，幾個事件的集合同時不發(fā)生時不會引起頂事件的發(fā)生，這個集合便稱為徑集。而在這個徑集中若任意去掉其中任意一個底事件，就不再是徑集，即要引起頂事件的發(fā)生，則這個徑集被稱之為最小徑集(Minimal Path Sets，MPS)。對于故障樹來說，它的全部最小徑集的完整集合代表了避免頂事件發(fā)生的各種安全模式。

2.2.3 結(jié)構(gòu)重要度 在故障樹中，每個底事件對故障頂事件有不同的影響，各底事件相對應的重要程度也有大小之分。因此，重要度分析就是一種研究底事件對故障頂事件的影響程度，并確定其重要性的分析方法［9］，有助于有效地改進系統(tǒng)的可靠度。結(jié)構(gòu)重要度是在不考慮各底事件的發(fā)生概率，或者說假設(shè)在各底事件的發(fā)生概率都是一樣的情

2.3 故障樹定量分析

故障樹定量分析方法［10］前提要求出各底事件發(fā)生的概率，從而利用最小割集與最小徑集來計算頂事件發(fā)生的概率。

1)利用最小割集計算頂事件發(fā)生概率:況下，從故障樹的邏輯關(guān)系對各底事件對故障頂事件的影響程度進行分析。這是一種定性的重要度分析［7］。

求解底事件的結(jié)構(gòu)重要系數(shù)是求解結(jié)構(gòu)重要度的一種重要的方法。當某個底事件的狀態(tài)從正常狀態(tài)0變?yōu)楣收蠣顟B(tài)1，而其他底事件的狀態(tài)保持不變時，則頂事件也從正常狀態(tài)0變?yōu)楣收蠣顟B(tài)1時如公式1，2所示:

其中，φ(0i，x)和φ(1i，x)分別代表底事件Xi的0狀態(tài)與1狀態(tài)下頂事件的狀態(tài)。2個公式表示這個底事件的狀態(tài)變化對頂事件的發(fā)生發(fā)揮了作用。在有m個底事件的故障樹中，兩種狀態(tài)的互不相容的組合數(shù)共有2m個。當把第Xi個底事件作為變化對象時，其余(m－1)個底事件狀態(tài)所對應保持不變的對照組共有2m－1個組合。在這2m－1個對照組中有多少對照組使頂事件的狀態(tài)發(fā)生變化的比值就是該事件Xi的結(jié)構(gòu)重要系數(shù)，用公式3表示:

2)利用最小徑集計算頂事件發(fā)生概率:

式(5)中，r，s分別表示最小徑集的序數(shù);Xi∈kr∪ks表示第i個底事件Xi或?qū)儆诘趓個最小徑集，或?qū)儆诘趕個最小徑集。1≤r＜s≤k表示r，s的取值范圍。

在實際情況中，要合理運用以上兩個計算頂事件發(fā)生概率的公式，若故障樹最小割集較少，則用公式(4)來計算，若最小徑集較少，則用公式(5)來計算。以上公式均是考慮底事件相互獨立。若底事件相關(guān)聯(lián)，則必須考慮相容事件和相斥事件的概率計算問題。另在故障樹分析中，有些很復雜很龐大的故障樹要精確求出頂事件的發(fā)生概率非常困難，這需要采用簡單、準確、快速的近似算法來替代最小割集最小徑集法，一般有求近似區(qū)間法，首項近似法以及獨立近似法等，在這里不做展開。

2.4 “脫鹽率下降且產(chǎn)水量上升”故障分析與維護

通過上述故障樹分析方法，對圖2故障樹的邏輯關(guān)系進行分析計算，“脫鹽率下降且產(chǎn)水量上升”故障T存在19個最小割集:K1={X2}，K2={X3}，X19}，K19={X1，X20}。19 個最小割集表示了造成系統(tǒng)“脫鹽率下降且產(chǎn)水量上升”的19種故障模式。同時存在2 個最小徑集個最小徑集表示了2 種保證系統(tǒng)正常運行的安全模式。并得出底事件的結(jié)構(gòu)重要度順序為

將如圖2所示的故障樹的底事件，按結(jié)構(gòu)重要度排列之后不難發(fā)現(xiàn)，誘發(fā)薄膜氧化所引發(fā)的故障底事件對頂事件“脫鹽率下降且產(chǎn)水量上升”故障的影響最大。在實際反滲透海水淡化系統(tǒng)運行中最需要防范以及注意這些故障的發(fā)生。

底事件的結(jié)構(gòu)重要度為防止故障的發(fā)生提供了參考。實際上，要對誘發(fā)“脫鹽率下降且產(chǎn)水量上升”故障的其他故障進行分析，相關(guān)故障若無法恢復，就會帶來嚴重的經(jīng)濟損失，諸如此類也應給予高度重視。

將“脫鹽率下降且產(chǎn)水量上升”故障對系統(tǒng)的影響嚴重程度分為3級。其中，系統(tǒng)組件不可恢復性損害故障為嚴重型，用“★★★”表示;系統(tǒng)中部分組件造成不可恢復性損害的故障為中度型，用“★★”表示;只需加以清洗、消毒等方式便能消去的故障為輕微型，用“★”表示;若不及時處理將會對系統(tǒng)中部分組件造成不可恢復性損害的故障用“★★☆”表示［11］。“脫鹽率下降且產(chǎn)水量上升”故障分析如表2所示。

表2 “脫鹽率下降且產(chǎn)水量上升”故障分析Tab.2 Fault analysis of decreases of desalination rate and water output

從表2的故障危害程度可以看到雖然引起“‘O’形圈泄漏”、“中心管損壞”、“中心疊片破裂”等這些中間故障事件的故障的底事件結(jié)構(gòu)重要度不高，但所造成的危害卻是非常嚴重、不易恢復的。所以實際的維護故障預防不但要考慮FTA的底事件的結(jié)構(gòu)重要度，還要根據(jù)故障可能會造成生產(chǎn)和經(jīng)濟上的損失而引起相關(guān)重視。

3 結(jié)語

通過對反滲透海水淡化工程工藝流程的研究，繪制了以“脫鹽率下降且產(chǎn)水量上升”故障為頂事件的故障樹，并分別提出了對故障樹進行定性定量分析的方法。在對繪制的故障樹進行分析后，求出最小割集、最小徑集和結(jié)構(gòu)重要度。將引起該故障的各種故障模式一一列出，并發(fā)現(xiàn)了最重要的事件為引起薄膜氧化的故障(X1，X2，X3)。

通過對故障的深入研究對引發(fā)“脫鹽率下降且產(chǎn)水量上升”故障的故障中間事件提出了相關(guān)的處理方法;并發(fā)現(xiàn)“膜組件機械損壞”(C2)等故障對反滲透海水淡化系統(tǒng)影響較大，且因此在海水淡化生產(chǎn)時應對上述故障加以密切關(guān)注，并給予定期的檢查與維護，保證反滲透海水淡化系統(tǒng)正常運行。故障樹分析法為及時檢查、排除故障并采取有效措施預防故障的發(fā)生提供了全面而可靠的信息，為反滲透海水淡化故障診斷專家系統(tǒng)的開發(fā)奠定了基礎(chǔ)。

［1］高從堦.海水淡化淺談［J］.化工技術(shù)經(jīng)濟，2003，21(11):24-26.GAO Cong-jie.Seawater desalination［J］.Chemical Techno-Economic，2003，21(11):24-26.(in Chinese)

［2］朱繼洲.故障樹原理和應用［M］.西安:西安交通大學出版社，1989.

［3］Fariz AbdulRahman，Athi Varuttamaseni，Michael Kintner，et al.Application of fault tree analysis for customer reliability assessment of a distribution power system［J］.Reliability Engineering and System Safety，2013，111:76-85.

［4］Michael Anthony，Robert Arno，Neal Dowling，et al.Reliability analysis for power to fire pump using fault tree and RBD［J］.Industry Applications，2013，49(2):997-1003.

［5］王小軍，周昌盛，呂漢興.一種基于關(guān)系數(shù)據(jù)庫的知識表示和推理方法［J］.控制工程，2005(12):8-12.WANG Xiao-jun，ZHOU Chang-sheng，LV Han-xing.Know ledge representation and reasoning method based on relation database［J］.Control Engineering of China，2005(12):8-12.(in Chinese)

［6］趙之剛，張海勇，張鋒.基于專家系統(tǒng)的艦載通信裝備故障診斷研究［J］.現(xiàn)代防御技術(shù)，2007(12):41-45.ZHAO Zi-gang，ZHANG Hai-yong，ZHANG Feng.Research on fault diagnosis expert system of warship communication equipment［J］.Modern Defence Technology，2007(12):41-45.(in Chinese)

［7］馮逸仙.反滲透水處理系統(tǒng)工程［M］.北京:中國電力出版社，2005.

［8］Karim Bourouni.Availability assessment of a reverse osmosis plant:comparison between reliability block diagram and fault tree analysis methods［J］.Desalination，2013，313:66-76.

［9］馮肇瑞，崔國璋.安全系統(tǒng)工程［M］.北京:冶金工業(yè)出版社，1987.

［10］姜周曙，胡亞才，屠傳經(jīng).FTA在溴化鋰吸收式制冷機故障診斷中的應用［J］.浙江大學學報:工學版，2001，31(2):209-213.JING Zhou-shu，HU Ya-cai，TU Chuan-jing.Application of FTA in the fault diagnose of lithium bromide-water absorption chillers［J］.Journal of Zhejiang University:Engineering Science，2001，31(2):209-213.(in Chinese)

［11］李朋.反滲透海水淡化故障診斷專家系統(tǒng)的研究與開發(fā)［D］.杭州:杭州電子科技大學，2011.