精細搜索策略應用于質量交換網絡綜合

作者：楊嶺崔國民周志強肖媛來源：《化工學報》日期：2022-11-03人氣：1100

隨著能源與環(huán)保問題日益嚴峻，我國提出了“碳達峰·碳中和”的戰(zhàn)略目標[1]，這對企業(yè)降低排放的能力提出新的要求，使得發(fā)展新的節(jié)能減排和清潔生產技術成為當下的研究熱點。質量交換網絡(mass exchange network, MEN)[2]作為一項廣泛應用于化工原料的提純、萃取、吸附和污染物去除技術，可顯著降低廢水、廢氣的產生和質量分離劑(mass separating agent, MSA)的消耗，近年來受到越來越多的關注。

質量交換網絡的概念由El-Halwagi等[3]于1989年提出，是在換熱網絡研究的基礎上發(fā)展而來的，因此許多學者[4-6]將換熱網絡的夾點法與數(shù)學規(guī)劃法運用到質量交換網絡上。El-Halwagi等[3]率先使用熱力學夾點法設計質量交換網絡，在濃度區(qū)間表中確定夾點，以此確定具有最大質量交換量的初始網絡，并根據(jù)初始網絡改進為成本最低的網絡。Hallale等[7-8]提出了一種新的y-y*復合曲線圖，用于處理涉及夾點同一側存在多個MSA的問題，并指出質量交換單元數(shù)量最小并不一定使成本最低，網絡設計時應充分考慮傳質驅動力。Gadalla[9]將夾點分析原理與質量交換網絡聯(lián)系起來，繪制富流股濃度與等效貧液流濃度的對比圖，所得圖用于分析當前網絡的性能，并設計新的網絡結構。Yanwarizal等[10]通過構建濃度與質量負荷圖同步確定最小MSA用量和貧富流股間的傳質負荷。

此外，數(shù)學規(guī)劃法在質量交換網絡綜合上也有著廣泛的應用。El-Halwagi等[11]使用兩階段合成方法，第一階段使用線性規(guī)劃確定外部MSA的最低成本和夾點，第二階段使用混合整數(shù)線性規(guī)劃減少質量交換器的數(shù)量，從而獲得綜合費用最低的網絡結構。Isafiade等[12]提出了一種基于區(qū)間的混合整數(shù)非線性規(guī)劃(mixed integer and non-linear programming, MINLP)超結構模型，可同時利用超結構與夾點技術的優(yōu)點，此外，還規(guī)定了分流必須以等濃度混合，雖然簡化了模型的復雜度，但一定程度上限制了模型的求解域。李紹軍等[13]建立了基于分級超結構的無分流質量交換網絡MINLP模型，并采用Alopex算法與遺傳算法結合的混合算法進行求解，可同時權衡操作費用和投資費用。謝會等[14]建立了一種基于區(qū)間的無分流超結構模型，其中過程貧流股與外部貧流股具有同等匹配機會，且允許重復匹配產生，并將列隊競爭算法應用于該模型。都健等[15]以傳質濃度差作為優(yōu)化變量，并考慮了多組分系統(tǒng)，利用改進的遺傳算法進行求解，取得了較好的優(yōu)化結果。侯創(chuàng)等[16]在數(shù)學模型中引入了取整函數(shù)，并通過LINGO求解器求解，解決了理論塔板數(shù)和實際塔板數(shù)存在差異的問題，提高了塔板的利用效率的同時降低了年綜合費用。

啟發(fā)式算法與MINLP問題的適應性和強大的搜索能力在許多研究中得以驗證，但其難以兼顧全局搜索和局部搜索，一定程度上限制了啟發(fā)式算法的求解能力。鑒于此，提出一種新的應用于質量交換網絡的精細搜索策略，分別利用了啟發(fā)式方法簡單有效、求解域大的特點和確定性方法精度高、收斂快的特點，建立兩種精細搜索方式，分別為帶有個體回代與分化的高精度強制進化隨機游走算法(random walk algorithm with compulsive evolution, RWCE)，以及坐標輪換法與黃金分割法結合的方法。最后通過兩個算例驗證了兩種精細搜索方法的有效性。

1 數(shù)學模型

1.1 問題描述

根據(jù)目標組分的濃度不同，將過程流股分為貧流股和富流股。為了使富流股中的組分濃度滿足出口濃度要求，通過質量分離器將富流股中所含目標組分通過吸收、解吸、吸附、萃取、離子交換等方式轉移到貧流股中。質量交換網絡優(yōu)化的目的是在滿足出口濃度要求以及傳質可行性的前提下，獲得最具經濟性的網絡結構。

如圖1所示，貧富流股之間的連接線表示傳質單元，并假設為逆流布置。第ir股富流股和第il股貧流股進行傳質操作，富流股從進口濃度y $i r i n$ 降低至出口濃度y $i r o u t$ ，貧流股從進口濃度x $i l i n$ 升高至出口濃度x $i l o u t$ 。M表示質量交換量，可由富流股的進出口濃度差與流量Gir 的乘積求得，亦可用貧流股濃度差與流量Lil 求得。

圖1

圖1 質量交換網絡示意圖

Fig.1 Sketch map of mass exchange network

貧富流股之間的傳質存在相平衡關系，如式(1)所示，對于稀溶液，在所涉及濃度范圍內，可假設相平衡關系是線性的，且與富流股中其他可溶性物質無關。

$y_{i r}^{*} = m_{i r, i l} x_{i l} + b_{i r, i l}$ (1)

式中，xil 為貧流股的濃度，此時能與貧流股傳質的富流股的最小濃度為Δy $i r *$ ；m和b為與貧富流股的性質有關的常數(shù)。圖2表示傳質過程中富流股的濃度變化以及貧流股的平衡濃度變化。其中Δy表示富流股濃度差；Δy*表示貧流股的平衡濃度差；Δyleft和Δyright分別為傳質單元兩端的傳質驅動力。

圖2

圖2 貧富流股的相平衡關系

Fig.2 Phase equilibrium of rich and lean stream

1.2 目標函數(shù)

目標函數(shù)是獲得質量交換網絡的最低年度綜合費用(total annual cost, TAC)，包括質量交換器的投資費用C1和外部貧流股的操作費用C2，如式(2)所示。操作費用取決于外部貧流股的流量，而對于不同工況的分離方式，投資費用計算方式不盡相同。當使用板式塔時，通過增加塔板數(shù)來增加質量交換器的傳質能力，此時的投資費用由塔板數(shù)決定，而使用填料塔時，投資費用由質量交換器的質量決定。操作費用如式(3)所示，投資費用具體計算方式在算例驗證部分詳細介紹。

$m i n T A C = C_{1} + C_{2}$ (2) $C_{2} = \sum_{i l = 1}^{N L} C_{0} L_{i l}, i l \in N L$ (3)

1.3 有分流節(jié)點非結構模型

以往的研究中常采用分級超結構(stage-wise superstructure, SWS)模型[17]合成HEN、MEN，而本文基于換熱網絡優(yōu)化的思想，將有分流節(jié)點非結構模型(nodes-based non-structural model, NNM)[18]應用于質量交換網絡的優(yōu)化，該模型以節(jié)點表示傳質單元潛在的位置，具有更強的靈活性和更寬廣的求解域。模型如圖3所示，結構中有NR股富流股和NL股貧流股，箭頭表示流股流動方向，所有質量交換單元均為逆流布置。富/貧流股上分別設置NRG和NLG個分流組，富/貧流股經過每個分流組時，分裂為NRS和NLS個分流，該組所有的分流在末端進行非等濃度混合后流入下一分流組。每個分流分別包含NRD和NLD個節(jié)點，因此所有的富流股節(jié)點數(shù)為NR×NRG×NRS×NRD，貧流股總節(jié)點數(shù)則為NL×NLG×NLS×NLD。初始狀態(tài)下，整個結構為空結構，優(yōu)化開始后隨機在富流股和貧流股上抽取兩個節(jié)點，連接形成質量交換單元。

圖3

圖3 有分流節(jié)點非結構模型

Fig.3 Schematic of the nodes-based non-structural model with streams split

1.4 約束條件

(1) 傳質可行性約束

$\begin{array}{l} y_{i r, j r, k r, l r}^{i n} > m_{i r, i l} x_{i l, j l, k l, l l}^{o u t} + b_{i r, i l} \\ i r \in N R, j r \in N R G, k r \in N R S, l r \in N R D \\ i l \in N L, j l \in N L G, k l \in N L S, l l \in N L D \end{array}$ (4) $y_{i r, j r, k r, l r}^{o u t} > m_{i r, i l} x_{i l, j l, k l, l l}^{i n} + b_{i r, i l}$ (5)

(2) 傳質單元質量守恒

$G_{i r, j r, k r} (y_{i r, j r, k r, l r}^{i n} - y_{i r, j r, k r, l r}^{o u t}) = L_{i l, j l, k l} (x_{i l, j l, k l, l l}^{o u t} - x_{i l, j l, k l, l l}^{i n})$ (6)

(3) 流股質量守恒

$G_{i r} (y_{i r}^{i n} - y_{i r}^{o u t}) = \sum_{i r = 1}^{N R} \sum_{j r = 1}^{N R G} \sum_{k r = 1}^{N R S} \sum_{l r = 1}^{N R D} [G_{i r, j r, k r} (y_{i r, j r, k r, l r}^{i n} - y_{i r, j r, k r, l r}^{o u t}) Z_{i r, j r, k r, l r}]$ (7) $L_{i l} (x_{i l}^{o u t} - x_{i l}^{i n}) = \sum_{i l = 1}^{N L} \sum_{j l = 1}^{N L G} \sum_{k l = 1}^{N L S} \sum_{l l = 1}^{N L D} [L_{i l, j l, k l} (x_{i l, j l, k l, l l}^{o u t} - x_{i l, j l, k l, l l}^{i n}) Z_{i l, j l, k l, l l}]$ (8)

(4) 節(jié)點出口濃度約束

$y_{i r, j r, k r, l r}^{i n} \geq y_{i r, j r, k r, l r}^{o u t}$ (9) $x_{i l, j l, k l, l l}^{o u t} \geq x_{i l, j l, k l, l l}^{i n}$ (10)

(5) 流股出口濃度約束

$y_{i r}^{o u t} \leq y_{i r}^{t a r g e t}$ (11) $x_{i l}^{o u t} \leq x_{i l}^{t a r g e t}$ (12)

(6) 流量約束

$\sum_{k r = 1}^{N R S} S P R_{i r, j r, k r} G_{i r, j r, k r} = G_{i r} \leq G_{i r, m a x}$ (13) $\sum_{k l = 1}^{N L S} S P L_{i l, j l, k l} L_{i l, j l, k l} = L_{i l} \leq L_{i l, m a x}$ (14)

2 強制進化隨機游走算法

本文基于節(jié)點非結構模型，采用RWCE算法[19]同步優(yōu)化質量交換網絡，該算法已在換熱網絡問題上取得了許多優(yōu)秀的成果[20-21]，驗證了該算法強大的全局搜索能力。RWCE算法的基本原理為：種群中的個體獨立進化，每個個體從空結構開始，以目標函數(shù)減小為進化方向，通過隨機生成傳質單元和隨機增減傳質量、分流比及流量實現(xiàn)整型變量和連續(xù)變量的同步進化，同時，有一定的概率接受目標函數(shù)較高的結構，具備跳出局部最優(yōu)解的能力。與以塔板數(shù)為優(yōu)化變量的算法不同，本算法以傳質量作為優(yōu)化變量，并根據(jù)質量衡算反求節(jié)點的進出口濃度，可避免迭代計算造成的優(yōu)化效率低下的問題，顯著提高優(yōu)化效率。

2.1 初始化

讀取算例的輸入?yún)?shù)，設定網絡的結構參數(shù)，包括富/貧流股的分流組數(shù)、分流數(shù)、節(jié)點數(shù)。設置種群規(guī)模Np，初始狀態(tài)下的所有個體均為空結構，并且結構中所有節(jié)點對應的變量信息全部歸零。

2.2 連續(xù)變量游走

需要進行游走操作的變量包括傳質量、分流比和貧流股流量。依次對所有節(jié)點做循環(huán)，若節(jié)點上存在傳質單元，則該傳質單元的連續(xù)變量有γ的概率參與游走。當判定要參與游走時，則分流比游走概率為β，傳質量游走概率為1-β，流量游走概率為δ。分流比、傳質量、貧流股流量分別按照式(15)~式(18)游走。游走公式中1-2r確定隨機的游走正負方向，rΔL確定隨機的游走步長，此公式保證步長是以0為中心的正態(tài)分布，即大概率以小步長搜索，小概率以大步長搜索。

$S P R_{i r, j r, k r}^{'} = S P R_{i r, j r, k r} + (1 - 2 r_{1}) r_{2} Δ L_{S P}$ (15) $S P L_{i l, j l, k l}^{'} = S P L_{i l, j l, k l} + (1 - 2 r_{3}) r_{4} Δ L_{S P}$ (16) $M_{i l, j l, k l, l l}^{'} = M_{i l, j l, k l, l l} + (1 - 2 r_{5}) r_{6} Δ L_{M}$ (17) $L_{i l}^{'} = L_{i l} + (1 - 2 r_{7}) r_{8} Δ L_{L}$ (18)

經過游走后，若一些傳質單元的傳質量、分流比或流量小于預設的下限Mmin、SPRmin、SPLmin、Lmin，認為這些傳質單元已無存在的必要，因此將這種傳質單元清除，傳質量、分流比、流量以及連接關系全部清零。

2.3 新單元生成

每次迭代有一定概率ε選中任意兩個貧富流股的節(jié)點，連接形成一個新的傳質單元，并賦予此傳質單元一個初始的傳質量。若該傳質單元所在的分流上沒有其他傳質單元則需要賦予其初始分流比，若該傳質單元所在的貧流股上沒有其他傳質單元，則還需賦予其初始流量，如式(19)~式(22)所示。

$M_{i r, j r, k r, l r}^{_{n e w}} = M_{m a x} r_{9}$ (19) $S P R_{i r, j r, k r}^{n e w} = \{\begin{array}{l} S P R_{m a x} r_{10}, i f S P R_{i r, j r, k r} = 0 \\ S P R_{i r, j r, k r}, e l s e \end{array}$ (20) $S P L_{i l, j l, k l}^{n e w} = \{\begin{array}{l} S P L_{m a x} r_{11}, i f S P L_{i l, j l, k l} = 0 \\ S P L_{i l, j l, k l}, e l s e \end{array}$ (21) $L_{i l}^{n e w} = \{\begin{array}{l} L_{m a x} r_{12}, i f L_{i l} = 0 \\ L_{i l}, e l s e \end{array}$ (22)

2.4 選擇與變異

如式(23)所示，當生成階段結束后，計算費用F1，并與游走前的費用F0進行比較，選擇費用較低的結構保留至下一次迭代中。然而，為了防止優(yōu)化陷入局部最優(yōu)，即便F1<F0，也仍有一定概率ζ接受費用較差的結構，避免算法早熟。

$S_{i t + 1} = \{\begin{array}{l} S_{i t} i f F_{1} \leq F_{0} \\ S_{i t} e l s e i f r_{13} < ζ \\ S_{i t - 1} e l s e \end{array}$ (23)

式中，向量 Sit 表示第it次迭代的網絡結構，包含所有節(jié)點的全部變量信息。

3 精細搜索策略的提出

3.1 全局搜索與局部搜索的矛盾

啟發(fā)式方法是一類具有貪婪性的隨機搜索方法，這決定了啟發(fā)式方法難以兼顧全局搜索和局部搜索。以RWCE算法優(yōu)化廢水脫酚算例為例進行說明，設定兩組不同精度的優(yōu)化參數(shù)，對比觀察二者所得結果的特點和差異，流股數(shù)據(jù)和投資費用計算公式見文獻[12]。保持模型結構參數(shù)相同，設定不同的優(yōu)化參數(shù)，參數(shù)A設定較大的步長、進化概率、生成概率以及接受差解概率，而參數(shù)B則將以上參數(shù)設定為較小值，經過相同次數(shù)的迭代后對比兩套參數(shù)所得的結果。當以參數(shù)A進行優(yōu)化時，在第4.205×108次迭代得到費用為338816 USD·a-1的結果并收斂于此，該結果的結構如圖4所示。從圖中可知，內部貧流股S1和S2的流量分別為4.9262 kg·s-1和2.2539 kg·s-1，而二者的流量限制分別為5.0000 kg·s-1和3.0000 kg·s-1，由于使用內部貧流股不會帶來費用的增加，因此應盡可能多地使用內部貧流股。由此可知，該結構的連續(xù)變量仍有優(yōu)化的空間，但受限于基礎算法的固定優(yōu)化參數(shù)，難以得到精度更高的結果，算法往往收斂于非局部最優(yōu)解。因此，出于保持全局搜索能力的目的而設定的固定參數(shù)不利于求解更精確的網絡結構，即局部搜索能力不足。

圖4

圖4 參數(shù)A得到的結構

Fig.4 The structure obtained with parameter A

當以參數(shù)B進行優(yōu)化時，如圖5所示，得到了一個包含4個傳質單元的結構。圖6為兩套參數(shù)優(yōu)化過程的TAC曲線，可以看出使用參數(shù)B時TAC下降頻次較少，且下降幅度較小，很快便收斂于368794 USD·a-1，此費用較前者更高。實驗結果說明，全程的高精度的優(yōu)化對全局優(yōu)化來說是不利的，不僅優(yōu)化效率降低，而且一旦陷入局部最優(yōu)解便難以跳出，導致算法的持續(xù)優(yōu)化能力不足。

圖5

圖5 參數(shù)B得到的結構

Fig.5 The structure obtained with parameter B

圖6

圖6 不同優(yōu)化參數(shù)的TAC迭代曲線

Fig.6 TAC iteration curves with different optimization parameters

分析可知，產生以上兩個實驗現(xiàn)象的原因如下：首先，MINLP問題的局部最優(yōu)解眾多，全局搜索和局部搜索難以兼顧，無論傾向于哪種搜索都會錯過一些更好的解；其次，整個優(yōu)化過程對精度的要求是不同的，優(yōu)化初期可以搜索到的可行解數(shù)量巨大，對算法全局搜索能力的需求更大，而優(yōu)化后期很難在全局范圍內搜索到更好的解，但在局部范圍內仍可以搜索到更好的解，因此對算法局部搜索能力的需求更大。鑒于此，可采用前期以較大的步長、進化概率、生成概率，以及接受差解概率進行全局搜索的方式，獲得一個初始解，再通過精度更高的方法對該解做進一步的優(yōu)化，以彌補初始解由于精度不足而未被充分優(yōu)化的問題。

3.2 精細搜索策略

提出兩種應用于質量交換網絡優(yōu)化的精細搜索方法。方法1的核心思想是：將基礎RWCE得到的解復制生成大量的新個體，每個新個體均具有相同的初始結構。以該初始結構為中心，在鄰域內以更精細化的隨機游走進行高精度搜索，并將所得到的更優(yōu)解作為初始解，再次復制給全部個體進行鄰域內的高精度搜索，以此實現(xiàn)更強的局部搜索能力。方法2的核心思想是：借助確定性方法的高精度和快速收斂的特性，對連續(xù)變量進行優(yōu)化。使用坐標輪換法[22]對基礎RWCE算法得到的結構進行搜索，該方法將多維的MEN問題轉化為一維優(yōu)化問題，并通過黃金分割法求解極小值點，可降低求解難度，同時也擁有較快的求解速度。

精細搜索是以基礎算法獲得的結構為起點的局部搜索，而局部搜索能獲得更好解的前提是基礎算法得到的結構處于較優(yōu)的解空間內。為了保證基礎RWCE算法能獲得具有優(yōu)化潛力的結構，本文通過算法參數(shù)的設置來強化初次優(yōu)化時的全局搜索能力，包括較大的游走步長、傳質單元進化概率、生成概率和接受差解概率，保證算法在全局范圍內獲得一個具有優(yōu)化潛力的結構，而非收斂于局部最優(yōu)解，在此基礎之上，再進行精細搜索便可得到更優(yōu)解。兩種精細搜索方法步驟如下。

方法1

（1) 初始化階段：與初始結構設置相同的模型參數(shù)，將結構中所有節(jié)點的傳質量、濃度、分流比和流量全部置零。

(2) 個體分化階段：讀取基礎RWCE算法所得最優(yōu)解的結構數(shù)據(jù)，并賦值給種群中的所有個體，將該結構的TAC記錄為F0。

(3) 迭代階段：種群中每個個體使用高精度的RWCE算法依次優(yōu)化，當?shù)螖?shù)達到ITf,max時，對比種群最優(yōu)費用Fp與初始費用F0，如果Fp小于F0，則將種群最優(yōu)個體的結構 Sp記為歷史最優(yōu)結構 SB，否則將初始結構 S0記為歷史最優(yōu)結構，最后將迭代次數(shù)歸零。

(4) 回代階段：將歷史最優(yōu)解回代到步驟(2)中，結構數(shù)據(jù)賦值給種群中的全部個體，作為下一輪迭代的初始結構。

值得注意的是，此方法是基于基礎RWCE實現(xiàn)的，保留了RWCE算法的游走、生成、選擇與變異等機制，在以局部搜索為主的同時也并未喪失全局搜索能力，因此不僅可以對連續(xù)變量做進一步優(yōu)化，也具有一定的結構變異能力。此外，方法1是基于隨機搜索的，收斂速度比較慢，雖然精度較基礎算法有所提高，但得到的解仍不是數(shù)學意義上的局部最優(yōu)，因此在對求解速度和精度有較高要求時，該方法存在一定的不足。操作流程如圖7所示。

圖7

圖7 基于啟發(fā)式算法的精細搜索策略流程

Fig.7 Flow chart of fine search strategy based on heuristic algorithm

方法2

MEN問題包含傳質量、分流比、流量三種連續(xù)變量，它們共同組成n維向量 x。本方法對優(yōu)化問題進行降維處理，依次沿著 x 的各個分量方向進行一維搜索，同時保持剩余n-1個分量不變。式(24)~式(26)表示以 xk-1為起點，沿著 ek 分量方向進行一維搜索，求取該方向上的極小值點和最佳步長λB，并將該極小值點作為下一分量的優(yōu)化起點，對剩下的分量采取相同的操作。當所有分量都優(yōu)化過一輪之后，再次從第一個分量開始優(yōu)化，直至某一輪優(yōu)化前后的目標函數(shù)值之差小于η時終止迭代。

$f (x_{k}) = m i n f (x_{k - 1} + λ e_{k})$ (24) $x_{k} = x_{k - 1} + λ_{B} e_{k}$ (25) $e_{k} = (0_{(1)}, 0_{(2)}, \dots, 1_{(k)}, \dots, 0_{(n)})^{T}$ (26)

其中極小值點求取關鍵在于λ的取值，本文采用的是黃金分割法求取每個分量的最佳步長。在初始搜索區(qū)間[a,b]內取兩個點α1和α2，其值如式(27)、式(28)所示，計算f (α1)和f (α2)。

$α_{1} = b - 0.618 (b - a)$ (27) $α_{2} = a + 0.618 (b - a)$ (28)

若f (α1)<f (α2)，則舍棄區(qū)間(α2,b]，并將區(qū)間[a,α2]作為新搜索區(qū)間[a1,b1]。若f (α1)>f (α2)，則舍棄區(qū)間[a,α1)，并將區(qū)間[α1,b]作為新搜索區(qū)間[a1,b1]。重復上述步驟，當搜索區(qū)間長度縮短至預先設定長度θ時，取區(qū)間中點α*作為極小值點，終止迭代，α*與起始點之差即為最佳步長。

基于啟發(fā)式算法的方法1是一種通用的精細搜索方法，可有效處理不同工況的算例，同時保持了一定的結構變異能力。而處理以填料塔為傳質單元的算例時，由于傳質單元投資費用是關于傳質量的連續(xù)函數(shù)，因此基于確定性方法的方法2具有較高的精度和較快的收斂速度。但方法2只能對連續(xù)變量進行優(yōu)化，并且只適用于在一個較優(yōu)解的基礎上做小范圍的調整，不適用于全局優(yōu)化。

本文求解環(huán)境為Intel(R) Xeon(R) Gold 6226R CPU，主頻2.9 GHz，64 GB RAM，采用Fortran90編程。

4 算例驗證

4.1 算例1 焦爐氣脫硫

本算例取自于文獻[3]，流股數(shù)據(jù)如表1所示，主要目的是去除焦爐氣中的硫化氫，以減少燃燒產生的SO2。過程貧流股S1采用稀氨水作為吸收劑，外部貧流股S2采用冷凍甲醇作為吸收劑。傳質單元采用篩板塔，投資費用如式(29)~式(34)所示。當相平衡方程為線性時，塔板數(shù)的近似計算通常使用Kremser方程求得，但由于方程形式和對數(shù)項的存在，計算過程中可能會產生分母為0的情況，這導致了公式的非連續(xù)性。因此本文采用Fraser等[23]提出的塔板數(shù)公式(30)，通過Underwood近似將Kremser方程轉化為連續(xù)函數(shù)，以此降低計算難度，其中n取1/3，每塊塔板的單價為4552 USD·a-1。貧流股總操作費用由單位操作費用乘以流量求得。RWCE算法參數(shù)為γ=0.3，β=0.2，δ=0.01，ε=0.002，ζ=0.01，ΔLSP=0.01，ΔLM=3×10-5，ΔLL=0.01，Mmax=0.002，Lmax=0.002，Mmin=10-5，SPRmin=0.01，SPLmin=0.01，Lmin=0.02。

表1 算例1的流股數(shù)據(jù)

Table 1 Stream data of case 1

流股	最大流量/(kg·s-1)	入口濃度/(kg H2S·kg-1)	目標濃度/(kg H2S·kg-1)	m	b	C0
R1	0.9000	0.07000	0.00030
R2	0.1000	0.05100	0.00010
S1	2.3000	0.00060	0.03100	1.45	0	117360
S2	∞	0.00020	0.00350	0.26	0	176040

新窗口打開| 下載CSV

$C_{1} = \sum_{i r = 1}^{N R} \sum_{j r = 1}^{N R G} \sum_{k r = 1}^{N R S} 4552 N_{i r, j r, k r}$ (29) $N_{i r, j r, k r} = {(\frac{Δ y_{i r, j r, k r}^{n} + Δ y_{i r, j r, k r}^{* n}}{Δ y_{r i g h t, i r, j r, k r}^{n} + Δ y_{l e f t, i r, j r, k r}^{n}})}^{1 / n}$ (30) $Δ y_{i r, j r, k r} = y_{i r, j r, k r}^{i n} - y_{i r, j r, k r}^{o u t}$ (31) $Δ y_{i r, j r, k r}^{*} = (m_{i r, i l} x_{i l, j l, k l}^{o u t} + b_{i r, i l}) - (m_{i r, i l} x_{i l, j l, k l}^{i n} + b_{i r, i l})$ (32) $Δ y_{l e f t, i r, j r, k r} = y_{i r, j r, k r}^{i n} - (m_{i r, i l} x_{i l, j l, k l}^{o u t} + b_{i r, i l})$ (33) $Δ y_{r i g h t, i r, j r, k r} = y_{i r, j r, k r}^{o u t} - (m_{i r, i l} x_{i l, j l, k l}^{i n} + b_{i r, i l})$ (34)

使用基礎RWCE優(yōu)化此算例，花費時長13278 s得到TAC=408996 USD·a-1的結果，包含4個傳質單元，總計22塊塔板，并且S1流股分為了兩股。此結果已低于現(xiàn)有文獻中的結果，但仍有繼續(xù)優(yōu)化的空間，因此使用精細搜索策略對該結果做進一步優(yōu)化。采用方法1經過25441 s獲得了TAC=407308 USD·a-1的結果，而采用方法2經過24 s獲得TAC=407801 USD·a-1的結果，方法1獲得了費用更低的結果。方法1的整個優(yōu)化過程的TAC曲線如圖8所示。基礎RWCE算法前期使TAC不斷下降，隨后TAC長期停滯，因此在5.000×108次迭代后使用方法1進行精細搜索，引入精細搜索策略后TAC快速下降，并在接下來的優(yōu)化過程中多次下降。方法1所得結果相比基礎RWCE算法所得結果TAC下降了1688 USD·a-1，相比文獻中的最優(yōu)結果TAC下降了3257 USD·a-1。精細搜索前后的MEN結構如圖9所示，相比基礎RWCE所得結果，S1的流量下降0.0027 kg·s-1，S2的流量增加了0.0181 kg·s-1，貧流股操作費用上升了2862 USD·a-1，但傳質單元R1-S1的塔板數(shù)減少了1塊，節(jié)省了4552 USD·a-1的投資費用，因此使TAC進一步下降。表2將本文結果與文獻中的結果進行了對比，得到了目前最優(yōu)的結果，由此可見，精細搜索策略可以實現(xiàn)對已有結構的深度優(yōu)化，提升了算法的局部搜索能力。

圖8

圖8 算例1優(yōu)化全過程的TAC迭代曲線

Fig.8 TAC iterative curve of the whole optimization process of case 1

圖9

圖9 算例1精細搜索前后的網絡結構

Fig.9 Network structure of case 1 before and after fine search

表2 算例1的結果對比

Table 2 Result comparison of case 1

文獻	單元數(shù)	總塔板數(shù)	TAC/(USD·a-1)
[24]	4	—	530471
[25]	5	—	469968
[7]	5	—	431613
[26]	4	25	429700
[27]	4	21	422293
[28]	4	24	420545
[15]	4	20	412500
[16]	4	22	411166
[13]	6	19	410971
[14]	6	19	410565
本文	4	21	407308

新窗口打開| 下載CSV

4.2 算例2 空氣除氨

本算例取自于文獻[29]，富流股由5股空氣組成，主要目的是從空氣中去除氨氣。S1、S2為過程貧流股，S3為外部貧流股，傳質單元采用填料塔，費用計算方法如式(35)~式(37)所示，流股數(shù)據(jù)如表3所示。

$C_{1} = \sum_{i r = 1}^{N R} \sum_{j r = 1}^{N R G} \sum_{k r = 1}^{N R S} (618 M A S S_{i r, j r, k r}^{0.66} \times 0.225)$ (35) $M A S S_{i r, j r, k r} = \frac{M_{i r, j r, k r}}{0.02 L M C D_{i r, j r, k r}}$ (36) $L M C D_{i r, j r, k r} = \frac{Δ y_{l e f t} - Δ y_{r i g h t}}{l n (Δ y_{l e f t} / Δ y_{r i g h t})}$ (37)

式中，618為單位高度投資費用；MASSir,jr,kr 表示設備質量；0.225為年度化系數(shù)，年操作時長為8150 h；0.02為總傳質系數(shù)；LMCDir,jr,kr 表示對數(shù)平均濃度差；Δyleft和Δyright分別表示傳質單元左右兩端的傳質驅動力。RWCE算法參數(shù)為γ=0.3，β=0.2，δ=0.1，ε=0.005，ζ=0.001，ΔLSP=0.03，ΔLM=3×10-4，ΔLL=0.1，Mmax=0.002，Lmax=1，Mmin=5×10-5，SPRmin=0.01，SPLmin=0.01，Lmin=0.02。

表3 算例2的流股數(shù)據(jù)

Table 3 Stream data of case 2

流股	最大流量/(kg·s-1)	入口濃度/(kg NH3·kg-1)	目標濃度/(kg NH3·kg-1)	m	b	C0
R1	2.0000	0.00500	0.00100
R2	4.0000	0.00500	0.00250
R3	3.5000	0.01100	0.00250
R4	1.5000	0.01000	0.00500
R5	0.5000	0.00800	0.00250
S1	1.8000	0.00170	0.00710	1.2	0	0
S2	1.0000	0.00250	0.00850	1	0	0
S3	∞	0.00000	0.01700	0.5	0	0.001

新窗口打開| 下載CSV

高志輝[28]采用自適應模擬退火遺傳算法，得到一個包含2次分流、9個傳質單元的結構，年綜合費用為120717 USD·a-1，但從給出的單元流股信息表中發(fā)現(xiàn)，5號單元的富端濃度不符合傳質可行性約束，因此結果存在誤差。本文采用基礎RWCE算法在1.845×108次迭代時得到TAC=129067 USD·a-1的解，用時21977 s，并且直到5.000×108次迭代也未能實現(xiàn)費用下降，因此使用精細搜索策略對該結構進一步優(yōu)化。采用方法1進行精細搜索時，經過1.751×108次迭代，用時19263 s，獲得了128065 USD·a-1的結果。采用方法2進行精細搜索用時51 s，經過1255次迭代，TAC下降至127807 USD·a-1，如圖10所示。在本算例中，方法2的速度明顯快于方法1，并且費用也有小幅下降，表現(xiàn)出更好的適用性。本文最優(yōu)結果相比基礎程序得到的結果下降了1260 USD·a-1，相比文獻[25]中的結果下降了2093 USD·a-1。精細搜索前后的網絡結構如圖11所示，由于本算例中使用填料塔作為傳質單元，因此圖中圓圈中的數(shù)字表示傳質單元的序號而非塔板數(shù)。

圖10

圖10 算例2基礎RWCE與精細搜索優(yōu)化過程的TAC曲線

Fig.10 TAC curves of basic RWCE and fine search optimization process of case 2

圖11

圖11 算例2精細搜索前后的結構

Fig.11 Network structure of case 2 before and after fine search

圖11兩個結構間的差異主要在貧流股S2和S3的流量，以及與之連接的傳質單元的傳質量。經過精細搜索后，S2的流量增加了0.0449 kg·s-1，與之相連的8號單元傳質量增加了0.00058 kg·s-1，與8號單元耦合的9號單元傳質量減少了0.00058 kg·s-1，與9號單元在同一組分流的1號單元傳質量增加了0.00028 kg·s-1，與1號單元耦合的2號單元減少了0.00028 kg·s-1，這導致了外部貧流股S3承擔的傳質負荷降低，因此所需流量也降低，使得操作費用降低了469 USD·a-1。此外，精細搜索后傳質量有了更好的分配方式，使得結構中設備總投資降低790 USD·a-1。與文獻結果的對比列于表4中，驗證了基于確定性方法的精細搜索策略對提升算法局部搜索能力的有效性。

表4 算例2的結果對比

Table 4 Result comparison of case 2

文獻	單元數(shù)	操作費用/(USD·a-1)	投資費用/(USD·a-1)	TAC/(USD·a-1)
[30]	8	85203	48797	134000
[12]	7	82410	50913	133323
[25]	9	81301	48599	129900
[28]	8	73350	47367	120717
本文	9	77508	50299	127807

新窗口打開| 下載CSV

5 結論

（1）低精度的優(yōu)化參數(shù)雖然全局搜索能力較強，但會導致連續(xù)變量難以得到充分優(yōu)化，無法收斂至局部最優(yōu)解。而高精度的優(yōu)化參數(shù)全局搜索能力較弱，陷入局部最優(yōu)解陷阱后便難以跳出，導致算法收斂于較差的結果。

（2）提出兩種用于對初始解深度優(yōu)化的精細搜索策略。方法1將最優(yōu)結構進行不斷回代和分化，并以高精度的RWCE算法進行局部搜索，具有適用范圍廣、優(yōu)化效果好的優(yōu)點，且保留了一定的全局搜索能力。方法2采用的坐標輪換法對每個變量的最優(yōu)解逐個求解，而最優(yōu)解通過黃金分割法求得，具有速度快、精度高的優(yōu)點，適合使用填料塔的實例。

（3）將精細搜索策略應用于焦爐氣脫硫和空氣除氨算例中，分別得到了407308 USD·a-1和127807 USD·a-1的優(yōu)化結果，與現(xiàn)有文獻結果相比經濟性提升，驗證了該策略的有效性，為費用最小質量交換網絡問題提供了一種新的行之有效的優(yōu)化方法。

符號說明

C0，C1，C2	分別為操作費用系數(shù)、投資費用、操作費用，USD·a-1
e	單位向量
F0，F1	分別為游走前后的年綜合費用，USD·a-1
IT	迭代次數(shù)
L，G	分別為貧、富流股流量，kg·s-1
LMCD	對數(shù)平均濃度差，kg·kg-1
ΔL	游走步長
M	傳質量，kg
MASS	傳質設備質量，kg
m，b	相平衡常數(shù)
N	塔板數(shù)
NL，NR	分別為貧、富流股的流股數(shù)
NLD，NRD	分別為貧、富流股的節(jié)點數(shù)
NLG，NRG	分別為貧、富流股的分流組數(shù)
NLS，NRS	分別為貧、富流股的分流數(shù)
Np	種群規(guī)模
r	區(qū)間(0,1)內的隨機數(shù)
S	質量交換網絡結構
SPL，SPR	分別為貧、富流股的分流比
x	連續(xù)變量組成的n維向量
x，y	分別為貧、富流股濃度，kg·kg-1
Z	0/1變量
α	區(qū)間邊界點
β	分流比進化概率
γ	進化概率
δ	流量進化概率
ε	新單元生成概率
ζ	接受差解概率
η	收斂誤差
θ	收斂區(qū)間長度
λ	一維搜索步長
上角標
in	流股進口
new	新傳質單元
out	流股出口
target	目標濃度
*	平衡濃度
下角標
B	歷史最優(yōu)
il，ir	分別為貧、富流股編號
it	迭代次數(shù)
jl，jr	分別為貧、富流股分流組編號
kl，kr	分別為貧、富流股分流編號
L	流量
ll，lr	分別為貧、富流股節(jié)點編號
M	傳質量
max	上限
min	下限
p	種群最優(yōu)
SP	分流比

關鍵字：優(yōu)秀論文

上一篇：基于知識數(shù)據(jù)化表達的制冷空調系統(tǒng)故障診斷方法
下一篇：深度融合特征提取網絡及其在化工過程軟測量中的應用

欄目分類

熱門排行

推薦信息

期刊知識