左 信，殷衛(wèi)兵，譚壯壯

（1.中國石油大學（北京），北京 102249；2.中國石油集團海洋工程有限公司，北京 100028）

進入21世紀，過程工業(yè)的規(guī)模不斷擴大，隨之而來的安全問題得到了人們的廣泛關注。石油化工等過程工業(yè)中，任何一個小的失誤都有可能帶來極其慘重的后果。為了降低事故發(fā)生的概率、減少事故發(fā)生帶來的損失，包括石油化工等越來越多的過程工業(yè)采用了安全儀表系統(tǒng)SIS（safety instrumented system）。SIS一般由傳感器、邏輯控制器和執(zhí)行器3部分組成，通常獨立于基本過程控制系統(tǒng)BPCS（basic process control system），專門用于處理緊急安全事故[1]。SIS系統(tǒng)的性能由多個性能指標衡量，應使用多目標優(yōu)化方法對SIS進行優(yōu)化。但是，傳統(tǒng)的多目標優(yōu)化方法在對SIS進行優(yōu)化時，存在最優(yōu)解分布不合理的問題，這是SIS優(yōu)化問題的特殊性導致的，因此本文提出一種改進的多目標優(yōu)化方法對SIS進行優(yōu)化。

1 安全儀表系統(tǒng)的多目標優(yōu)化

近年來，對于安全儀表系統(tǒng)的研究主要集中于可靠性模型、共因失效分析、冗余結(jié)構(gòu)優(yōu)化等方面[2-4]。衡量SIS系統(tǒng)的主要性能指標有3個：安全完整性水平SIL（safety integrity level）、誤跳車率STR（spurious trip rate）以及生命周期成本LCC（life circle cost）。SIL是安全儀表系統(tǒng)最重要的性能指標，它表征了安全儀表系統(tǒng)所能達到的可靠性等級，根據(jù)IEC 61508的要求[1]，SIL分為4級，劃分依據(jù)為平均要求時失效概率PFDavg（average probability of dangerous failure on demand）和系統(tǒng)結(jié)構(gòu)，表1為SIL等級與PFDavg的對應關系[1]。SIS系統(tǒng)運行過程中，實際生產(chǎn)過程并不存在安全隱患，而由于SIS自身原因?qū)е铝隋e誤跳車，即SIS系統(tǒng)進行了誤動作，帶來了不必要的跳車損失，這種事件發(fā)生的概率可以用STR來描述。另外，LCC也是SIS系統(tǒng)在設計和運行過程中應該考慮的一個重要指標。

表1 低要求操作模式下的安全完整性水平Tab.1 Request of safety integrity level

由前述可知，一個好的SIS系統(tǒng)需要在多個方面表現(xiàn)良好。然而PFDavg、STR和LCC是存在矛盾的3個指標，即某一指標的提高往往會帶來其它指標的下降，所以，對于SIS系統(tǒng)的多目標優(yōu)化成為了安全儀表系統(tǒng)的一個研究領域[3-5]。在以往的研究中[3，5]，主要是對SIS系統(tǒng)的PFDavg、STR和LCC對于冗余結(jié)構(gòu)MooN（M-out-of-N）和測試周期TI（Test Interval）的計算模型進行多目標優(yōu)化。但是，在實際的設計要求中，往往對于可靠性的要求是系統(tǒng)的SIL達到某一等級。因此可以認為在SIS系統(tǒng)的多目標優(yōu)化中，應該在達到SIL要求的前提下，即PF Davg達到某一SIL等級的要求時，使得STR和LCC盡可能小。然而，原有的多目標遺傳算法選擇解集的一個原則是使解集中的解均勻的分布在Pareto前沿上，而并非更多的達到某一值，所以針對這一具體問題的特殊要求，筆者通過改進原有的多目標遺傳算法，在解集中解的數(shù)量不變的情況下得到更加符合實際應用要求的解。

2 安全儀表系統(tǒng)的多目標優(yōu)化模型

安全儀表系統(tǒng)的3個主要指標為：PFDavg、STR和LCC。這3個指標主要由冗余結(jié)構(gòu)MooN和測試周期TI決定，因此建立這3個指標對于MooN和TI的計算模型[2，5]。

2.1 安全儀表系統(tǒng)的PFDavg模型

系統(tǒng)整體的PFDavg為

式中：P，PS，PL，PA分別為整個系統(tǒng)、傳感器子系統(tǒng)、邏輯控制器子系統(tǒng)、執(zhí)行器子系統(tǒng)的平均要求時失效率。其中，每個子系統(tǒng)有獨立的MooN冗余結(jié)構(gòu)。

對于每個子系統(tǒng)，P由獨立失效部分PI和共因失效部分PC組成，即：

式中：獨立失效部分PI表示子系統(tǒng)的失效僅僅是某一個通道的失效而產(chǎn)生的，在MooN系統(tǒng)中只有當獨立失效的通道數(shù)大于N減去M的差時子系統(tǒng)才會失效；共因失效部分PC表示通道的失效是由于共同原因造成的，這種失效一旦發(fā)生往往意味著整個系統(tǒng)的失效。

2.2 安全儀表系統(tǒng)的STR模型

STR的一般計算模型：

式中：S，SS，SL，SA分別為整個系統(tǒng)、傳感器子系統(tǒng)、邏輯控制器子系統(tǒng)、執(zhí)行器子系統(tǒng)的誤跳車率。對于每個子系統(tǒng)，S由獨立失效部分SI和共因失效部分SC組成，即：

式中：獨立失效部分SI表示誤跳車僅是某一個通道的安全失效而引起的，在MooN系統(tǒng)中只有當獨立安全失效的通道數(shù)大于或等于M時發(fā)生誤跳車；共因失效部分SC表示誤跳車是由共同原因造成的，這種失效一旦發(fā)生往往意味著誤跳車的發(fā)生。

2.3 安全儀表系統(tǒng)的LCC模型

生命周期成本L計算模型為

其中：Cc為固定成本，元；Co為運行成本，元；Ci為風險成本，元。固定成本Cc包括系統(tǒng)的設計、購買、安裝、啟動等由系統(tǒng)結(jié)構(gòu)決定的費用；運行成本Co包括檢測成本、維護成本，主要有系統(tǒng)結(jié)構(gòu)和測試周期決定；風險成本Ci包括由危險事故和誤跳車產(chǎn)生的費用，主要有系統(tǒng)的PFDavg和STR決定。

以上3個模型的具體計算模型可見文獻[5]。

3 改進的安全儀表系統(tǒng)的多目標優(yōu)化算法

3.1 安全儀表系統(tǒng)多目標優(yōu)化問題的數(shù)學描述

安全儀表系統(tǒng)的 3個主要性能指標 PFDavg、STR和LCC，是由冗余結(jié)構(gòu)MooN和測試周期TI決定的。因此，可以將安全儀表系統(tǒng)的優(yōu)化問題描述為以MooN和TI為決策變量，PFDavg、STR和LCC為目標函數(shù)的多目標優(yōu)化問題，如式（6）所示：

式中：下標i表示第i個子系統(tǒng)；T為測試周期，月；Tmin和Tmax分別為測試周期下限和上限，月；N和M表示冗余結(jié)構(gòu)為N取M結(jié)構(gòu)；Nmax為最大冗余度。

3.2 非支配排序遺傳算法

PFDavg、STR和LCC是安全儀表系統(tǒng)的3個主要性能指標，單一的對其中某個優(yōu)化并不能為設計安全儀表系統(tǒng)帶來較高的參考價值。因此，安全儀表系統(tǒng)有多個指標需要同時進行優(yōu)化。目前，常用的多目標優(yōu)化方法有2種：第1種是將多個目標函數(shù)加權(quán)后相加得到一個目標函數(shù)，再采用單目標優(yōu)化的方法對其進行優(yōu)化，這種方法最終得到一個單一的解，這個解往往在某一方面很好，但在其它方面可能很差，這樣的解沒有太多的參考價值；第2種方法是采用多目標優(yōu)化算法，其最終產(chǎn)生一組位于Pareto最優(yōu)前沿上的解，這組解多數(shù)情況下含有多個解，彼此之間不互相支配，可以供決策者選擇。

在多目標優(yōu)化中，優(yōu)化目標是一個由多個目標函數(shù)組成的函數(shù)向量。假設所有的優(yōu)化目標函數(shù)都是要最小化，對于2個可行解x1和x2，所有的目標函數(shù)fl都滿足fl（x1）≤fl（x2），則稱解x1支配解x2。如果一個解集中的所有解都不被解空間中的其他解支配，則該解集為Pareto最優(yōu)前沿。多目標優(yōu)化的最終目標就是求取問題的Pareto最優(yōu)前沿。由于冗余分配問題被證明為NP-hard問題[6]，因此，文中采用了非支配排序遺傳算法NSGA-II[7]，它是一種高效的多目標優(yōu)化算法[8]。

NSGA-II的核心思想是將當前世代的所有個體按照適應度進行排序，排序靠前的個體獲得更高的參與遺傳操作的概率，而適應度由群體分級（Rank）和計算擁擠距離（Crowding Distance）決定。分級機制的原理是：種群中不被其它個體支配的個體獲得Rank為1的分級，只被Rank為1的個體支配的個體獲得Rank為2的分級，只被Rank為1和2的個體支配的個體獲得Rank為3的分級，以此類推。擁擠距離為個體在目標函數(shù)空間中前后相鄰的個體之間的距離，如圖1所示，為一個含有2個目標f1和f2的多目標優(yōu)化問題，圖中x1，x2，x3和x4為遺傳算法中的4個個體，其中x2的擁擠距離為d1，即x2鄰近個體x1和x3之間的距離，同樣x3的擁擠距離為d2。在NSGA-II中，Rank值低的個體會有更高的幾率參與遺傳操作，同Rank值的個體中擁擠距離大的個體優(yōu)先參加遺傳操作，由此可見NSGA-II算法選取解的策略是讓最終的解集盡可能均勻分布在Pareto最優(yōu)前沿上。

利用改進的多目標遺傳算法對安全儀表系統(tǒng)進行優(yōu)化的基本步驟如下：

步驟1建立以PFDavg、STR、LCC為優(yōu)化目標，Ni，Mi，Ti為決策變量的目標函數(shù)，如式（6）所示。

步驟2生成初始種群，并為隨機生成每個個體的｛Ni，Mi，Ti｝。

步驟3計算個體的適應度｛P，S，L｝。計算個體的Rank和擁擠距離值，以進行非支配排序。

步驟4判斷遺傳代次是否等于最大代次，若滿足進入步驟5；否則，遺傳代次加1，并執(zhí)行交叉、變異以及選擇等遺傳算子，生成新的種群，重復步驟3。

步驟5選擇種群中Rank=1的個體，得到Prato最優(yōu)前沿。

圖1 NSGA-II算法擁擠距離示意圖Fig.1 Crowding Distance of NSGA-II

3.3 NSGA-II算法應用于安全儀表系統(tǒng)優(yōu)化的缺陷

NSGA-II算法以及其它非精確算法得到的多目標最優(yōu)解，并不是完整的Pareto解集。如NSGA-II算法等基于遺傳算法的多目標優(yōu)化算法，其最優(yōu)解集中解的個數(shù)最大為種群中個體的數(shù)量，當實際Pareto解集中解的數(shù)量大于個體數(shù)量時，算法將無法得到完整的Pareto解集。如圖2所示為一個2目標優(yōu)化問題的Pareto解集和NSGA-II產(chǎn)生的最優(yōu)解集。

圖中實心圓為NSGA-II產(chǎn)生的解，空心圓為Pareto解集中NSGA-II未產(chǎn)生的解。從圖中可以看出Pareto解集中有18個解，而NSGA-II只產(chǎn)生了9個解，這是由于設定優(yōu)化參數(shù)時將種群數(shù)量設定成了9。上述問題可以將種群數(shù)量設定成18以上來解決。但是，實際中的問題的Pareto解集往往有上千甚至上萬個解，而單靠增加NSGA-II的種群數(shù)量是不可行的，因為當種群數(shù)量增大時，計算的復雜度將會變得很大。因此，多目標優(yōu)化算法所產(chǎn)生的解往往不能產(chǎn)生完整的Pareto解集，只能產(chǎn)生一定數(shù)量的解來描述Pareto解集。如前所述，NSGA-II采用秩和擁擠距離結(jié)合的方法來產(chǎn)生最優(yōu)解集，其目的就在于保持種群的多樣性，同時使最終的解均勻的分布在Pareto前沿上，這樣來達到用優(yōu)化最優(yōu)解盡可能完整的描述Pareto最優(yōu)解集的目的。

采用擁擠距離的排序策略，可以很好的使NSGA-II產(chǎn)生的最優(yōu)解均勻的分布在Pareto最優(yōu)前沿上，這對于大部分多目標優(yōu)化問題來說是很好的結(jié)果。但是，SIS系統(tǒng)存在其特殊性。在實際的設計中，常常要求的是SIS系統(tǒng)的SIL等級，即對于可靠性的要求是系統(tǒng)的SIL達到某一等級，而并非PFDavg最小。因此可以認為在SIS系統(tǒng)的多目標優(yōu)化中，應該在達到SIL要求的前提下，即PFDavg達到某一SIL等級的要求時，使得STR和LCC盡可能小。然而，原有的多目標遺傳算法選擇解集的一個原則是使解集中的解均勻的分布在Pareto前沿上，而并非更多的達到某一值，所以針對這一具體問題的特殊要求，本章通過改進原有的多目標遺傳算法，在解集中解的數(shù)量不變的情況下得到更加符合實際應用要求的解。

3.4 改進的安全儀表系統(tǒng)多目標優(yōu)化算法

NSGA-II是一種基于遺傳算法的多目標優(yōu)化算法，因此，由NSGA-II產(chǎn)生的解集中解的數(shù)量不會超過其種群數(shù)量，換言之，大多數(shù)情況下所產(chǎn)生的解集是Pareto最優(yōu)前沿的子集。在安全儀表系統(tǒng)的多目標優(yōu)化問題中，Pareto最優(yōu)前沿中一般有上千個解，但是如果把NSGA-II的種群數(shù)量選為1000以上，算法的運算復雜度將會大的無法接受，事實上一般種群數(shù)量都在100左右。如何在有限的種群數(shù)量下盡可能完整的表現(xiàn)Pareto最優(yōu)前沿成為了多目標遺傳算法存在的一個問題。如前所述，NSGA-II使用群體分級機制和計算擁擠距離相結(jié)合的方法來解決這一問題，使優(yōu)化結(jié)果在保持非支配的條件下均勻的分布在Pareto最優(yōu)前沿上。

這種策略在大多數(shù)問題中是適合的，但是在安全儀表系統(tǒng)優(yōu)化中，往往希望得到的是PFDavg略微優(yōu)于某一SIL等級要求的水平，而其它2個指標更好的解，換言之，在多數(shù)設計過程中更加關注的是系統(tǒng)的SIL等級而非實際的PFDavg值。如表2所示，2個解都達到了SIL3的要求，第2個解的PFDavg比SIL3低很多，但是還達不到SIL4；第1個解PFDavg只比SIL3的要求低一點，但卻具有優(yōu)于第2個解的STR和LCC。實際上，PFDavg、STR和LCC這3個指標大多數(shù)時候是矛盾的，所以PFDavg較好往往意味著另2個解較差，這一點可以從表2的對比看出。因此，應當在優(yōu)化的過程中使得具有略低于某一SIL等級要求（即10-1，10-2，10-3或10-4）的PFDavg的解更容易遺傳到后代，即剛剛達到某一SIL等級對PFDavg要求但STR和LCC更好的解。

為達到此目的，在擁擠距離中加入PFDavg與SIL要求值的距離的項，使得個體按距離排列的過程中同時按照PFDavg與SIL要求值的距離排列。如下式所示：

表2 低要求操作模式下的安全完整性水平Tab.2 Request of safety integrity level

圖2 Pareto解集和NSGA-II產(chǎn)生的最優(yōu)解集Fig.2 Pareto solution set and NSGA-II optimal solution set

式（7）中：D（x）為原NSGA-II中解x的擁擠距離的計算公式，xpre與 xnext為與解 x鄰近的解；式（8）中：D′（x）為改進算法中解x的擁擠距離計算公式，其中z為SIL等級的邊界值，如表1所示。加入這個距離后，遺傳算法會向著具有更接近SIL等級要求的方向優(yōu)化，這樣可以得到更多的有利于決策的解。用以上獲得的值代替原NSGA-II算法中步驟4中的Distance值，將使得PFDavg略優(yōu)于某一SIL等級要求的個體更容易參與遺傳操作。

改進的多目標優(yōu)化算法能夠?qū)⒂邢薜慕饧性诟鱾€SIL等級附近，獲得更有利于設計人員作出決策的解集，對于安全儀表系統(tǒng)的優(yōu)化有很大價值。

4 應用實例

為對比提出的改進優(yōu)化方法與原優(yōu)化方法的不同，以文獻[5]中的例子為例，同時應用提出的改進算法和原算法對安全儀表系統(tǒng)冗余分配問題進行優(yōu)化。例子中的安全儀表系統(tǒng)為壓力保護系統(tǒng)，其示意圖如圖3所示[3]。

如圖3中，PT為傳感器，PLC為邏輯控制器，F(xiàn)C為最終執(zhí)行機構(gòu)，M/N代表MooN冗余結(jié)構(gòu)。例子中壓力保護系統(tǒng)用于保護化學反應器，使當中的壓力不致過高。傳感器子系統(tǒng)用于測量化學反應器中的壓力，傳送入邏輯控制子系統(tǒng)；邏輯控制子系統(tǒng)通過預設的程序?qū)ο到y(tǒng)狀態(tài)進行判斷，當化學反應器中的壓力高于預設的安全壓力時發(fā)出動作信號；執(zhí)行器子系統(tǒng)接受邏輯控制子系統(tǒng)的控制信號并產(chǎn)生相應的動作。本文使用文獻[5]中的模型作為安全儀表系統(tǒng)的計算模型，并采用相同的實驗數(shù)據(jù)進行實驗。

圖3 壓力保護系統(tǒng)示意圖Fig.3 Schematic diagram of pressure protection system

運用提出的改進優(yōu)化方法對上述問題進行優(yōu)化，并與原有的NSGA-II優(yōu)化算法[7]進行比較。優(yōu)化結(jié)果如圖4和圖5所示。

圖4 兩種優(yōu)化方法的優(yōu)化結(jié)果三維圖Fig.4 Optimization results of the two methods

圖5 兩種優(yōu)化算法的優(yōu)化結(jié)果二維圖Fig.5 Optimization results.using two algorithms

圖4所示為原有的方法和改進方法產(chǎn)生的解在以PFDavg、STR和LCC為坐標軸的三維圖中的分布，從圖中可以看出2種方法都產(chǎn)生了非支配的解集；原有的方法產(chǎn)生解均勻的分布在三維空間中，而改進方法產(chǎn)生的解集中在SIL等級要求的PFDavg的邊緣附近。

圖5為改進的算法和原有的NSGA-II算法的比較，所有的最優(yōu)解按照PFDavg進行降序排列。從圖中可以看出2種方法產(chǎn)生的編號為1到編號為34的解，均滿足SIL2的要求且LCC相差不多，但是改進方法的STR優(yōu)于原方法；2種方法產(chǎn)生的編號為37到編號為54的解，均滿足SIL3的要求，但是改進方法的STR和LCC優(yōu)于原方法；2種方法產(chǎn)生的編號為59到編號為90的解，均滿足SIL4的要求，但是改進方法的STR和LCC優(yōu)于原方法。

因此，可以得出原有方法所產(chǎn)生解的PFDavg均勻的分布在10-5到10-2的區(qū)間內(nèi)，改進方法所產(chǎn)生解的PFDavg大多分布在略微優(yōu)于SIL等級要求的PFDavg值的附近；同處于一個SIL等級的解中，改進方法所產(chǎn)生的解在STR和LCC 2個方面大多優(yōu)于原有的方法。

5 結(jié)語

優(yōu)良的安全儀表系統(tǒng)應在安全完整性水平（SIL），誤跳車率（STR）以及生命周期成本（LCC）3個方面達到綜合最優(yōu)化要求，通過采用改進的安全儀表系統(tǒng)多目標優(yōu)化算法，得到了更好的滿足安全完整性水平偏好綜合要求的設計方案。改進的安全儀表系統(tǒng)多目標優(yōu)化算法體現(xiàn)出良好的優(yōu)越性：首先，改進的安全儀表系統(tǒng)多目標優(yōu)化算法能夠使得到的解具有略微優(yōu)于某一SIL等級的PFDavg值；其次，改進算法產(chǎn)生的解與原有NSGA-II方法產(chǎn)生的解相比較，在具有相同SIL等級時具有更優(yōu)的STR和LCC值；第三，改進算法更加適應安全儀表系統(tǒng)的多目標優(yōu)化問題，能夠為設計人員提供更有意義的參考。

[1] IEC 61508-1998.Functional Safety of Electrical/Electronic/ Programmable Electronic Safety Related Systems[S].1998.

[2] 王麗君，左信，潘高峰.安全儀表系統(tǒng)的不同運行階段實時可靠性分析[J].計算機與應用化學，2012，29（11）：1331-1334.

[3] Torres-Echeverría A C，Martorell S，Thompson H A.Design optimization of a safety-instrumented system based on RAMS+ C addressing IEC 61508 requirements and diverse redundancy [J].Reliability Engineering and System Safety，2009，94（2）：162-179.

[4] Torres-Echeverría A C，Martorell S，Thompson H A.Multiobjective optimization of de sign and testing of safety instrumented systems with MooN voting architectures using a genetic algorithm[J].Reliability Engineering and System Safety，2012，106（10）：45-60.

[5] 左信，譚壯壯.安全儀表系統(tǒng)冗余分配的多目標優(yōu)化[J].中國安全科學學報，2013，23（11）：37-42.

[6] Maw-Sheng Chern.On the computational complexity of reliability redundancy allocation in a series system[J].Operations Research Letters，1992，11（5）：309-315.

[7] Kalyanmoy Deb，Samir Agrawal，Amrit Pratap，et al.A fast and elitistmulti-objective genetic algorithm：NSGA-II[J].IEEE Transactions on Evolutionary Computation，2002（2）：182-197.

[8] Carlos A Coello.Evolutionary multi-objective optimization：a historical view of the field[J].Computational Intelligence Magazine，IEEE，2006（1）：28-36.