基于改進蟻獅算法的間歇反應釜優化方法研究

汪會　潘海鵬　張益波

摘? 要：針對蟻獅算法易陷入局部極值點和收斂速度慢的問題，提出一種基于自適應邊界調節策略和分段搜索策略的改進型蟻獅算法。通過引入自適應調節因子對螞蟻圍繞蟻獅游走時的范圍進行自適應改變，提高了算法的全局尋優能力和尋優精度。融合粒子群算法的記憶保存思想和柯西變異算子對螞蟻位置進行分段搜索，豐富了螞蟻種群的多樣性，解決了算法易陷入局部最優的問題。將改進的算法應用于乙酸乙酯生產過程的間歇式反應釜，并與蟻獅算法進行比較，實驗表明改進后算法收斂速度更快，對間歇式反應釜生產過程中的反應溫度、產物濃度和反應時間等參數的優化效果明顯。

關鍵詞：改進蟻獅算法;間歇式反應釜;自適應邊界調節策略;分段搜索策略;參數優化

中圖分類號：TP273? ? ?文獻標識碼：A

Research on Optimization Method of Batch Reactor

based on Improved Ant Lion Algorithm

WANG Hui， PAN Haipeng， ZHANG Yibo

（School of Mechanical and Automatic Control， Zhejiang Sci-Tech University， Hangzhou 310018， China）

1119549054@qq.com; pan@zstu.edu.cn; zhangy41@163.com

Abstract： This paper proposes an improved ant lion algorithm based on adaptive boundaries and segmentation search strategy to address the problems that ant-lion algorithm tends to fall into local extremum and slow convergence speed. The adaptive adjustment factor is introduced to adaptively change the range of ants around the ant lion， which improves global search ability and search accuracy of the algorithm. The idea of memory preservation of particle swarm optimization and Cauchy mutation operator are combined to search the ant's position in segments， which enriches the diversity of ant population and solves the problem that the algorithm is easy to fall into local optimum. The improved algorithm is applied to the batch reactor of ethyl acetate production process. Compared with the traditional ant lion algorithm， experimental results show that the improved algorithm converges faster， and it has a better optimization effect on parameters， such as reaction temperature， product concentration and reaction time in the batch reaction production process.

Keywords： improved ant lion algorithm; batch reactor; adaptive boundary adjustment strategy; segmented search

strategy; parameter optimization

1? ?引言（Introduction）

間歇生產過程是將有限量的物料按規定的加工順序，在一個或多個設備中加工，達到一定的反應程度后一次性取出全部反應物料，獲得有限量產品的加工過程[1-3]，被廣泛應用于要求多品種、小批量、優質高產等的生產領域。與連續過程相比，間歇生產過程因其工作點不穩定、有一定的運行時間和重復運行的特點，它的優化方法研究受到廣泛關注[4]。

乙酸乙酯生產是一個典型的間歇生產過程，含有三個生產單元，分別是反應單元、中和單元和精餾單元[5]。反應單元是乙酸乙酯生產的第一個工段，反應生成的乙酸乙酯濃度直接影響后面精餾單元提純的產品質量，且反應單元和中和單元耗時較長，而精餾單元的耗時較短，因此生產過程中各單元之間存在節拍失衡問題，因此對反應單元的產物濃度和反應時間的優化十分必要。

由乙酸乙酯反應單元的特性看出，該單元的模型復雜，求解難度較大，因而利用數值求解是更優的選擇，如無須梯度信息的迭代動態規劃[6]、遺傳算法[7]和粒子群算法[8]等。李宏光等[9]采用迭代動態規劃算法，應用于終端時刻固定的間歇生產過程的動態優化問題中，但迭代動態規劃算法需對狀態變量進行離散化，計算量很大。趙博等[10]同時考慮費用和環境影響，利用遺傳算法求解間歇化工過程的多目標優化問題，遺傳算法的計算量由種群大小決定，即由控制變量決定，當遇到一些控制變量較多的過程時，遺傳算法的計算工程也會十分復雜。陳偉等[11]針對序列二次規劃算法求解含復雜約束的間歇生產過程時全局搜索弱的問題，采用粒子群算法增強了其對含復雜約束的優化問題的求解能力，但粒子群算法容易產生早熟收斂，且局部尋優能力不足。

蟻獅優化算法（Ant Lion Optimizer， ALO）[12]是Mirjalili于2015 年提出的一種仿生智能算法，具有調節參數少、尋優精度高等優點，成為近幾年進化計算領域的研究熱點之一，已被廣泛應用于變壓器設計、電力系統、無線電傳感器以及無人機航線規劃等領域。徐欽帥等[13]為了提高無線傳感器網絡的節點覆蓋率，引入邊界收縮因子和判斷機制而提出混合變異的改進蟻獅算法（Mixed Strategy based ALO， MS-ALO），提高了網絡覆蓋率，并使節點分布更加均勻。劉景森等[14]基于自適應邊界算子和動態比例系數提出了優選策略的自適應蟻獅優化算法（Preferred Strategy Self-adaptive ALO， PSALO），取得了較好的尋優效果，但他們只在測試函數上進行了仿真，沒有從工程角度進行相關實驗。

本文以間歇式反應生產過程為對象，以產品質量、反應時間和能耗為指標，針對ALO收斂速度慢且易陷入局部極值問題，引入自適應邊界調節策略和分段搜索策略[15]，提出一種基于自適應邊界和柯西變異的改進蟻獅優化（Ant Lion Optimization with Self-adaptive and Cauchy Mutation， CSALO）算法。通過對螞蟻圍繞蟻獅隨機游走的范圍進行自適應調整，隨機改變螞蟻游走范圍的大小，增強了算法的全局尋優能力和尋優精度。在螞蟻搜索可行解過程中引入分段搜索策略，搜索的前半段參考粒子群算法的記憶保存策略更新螞蟻位置，增加算法的多樣性，提高算法的局部搜索能力;算法的后半段加入柯西變異算子，提高了算法尋到全局最優解的概率，增強了算法跳出局部極值點的能力，并加快了算法尋優速度。最后將CSALO應用于乙酸乙酯生產的間歇式反應單元的參數優化問題，通過五組仿真實驗證明，與蟻獅算法相比，本文提出的方法收斂速度更快，在保證反應時間較短的情況下，提高了乙酸乙酯的濃度，降低了生產能耗，優化效果明顯。

2? ?蟻獅優化算法（Ant lion optimization algorithm）

蟻獅是一種肉食性昆蟲，以昆蟲為食，以其獨特的捕食方式而得名。蟻獅利用它尖銳的下顎在沙子里設下陷阱，等到覓食的螞蟻掉入陷阱后，蟻獅就會將其吞食。捕食成功后，蟻獅挖下一個陷阱繼續等待捕捉螞蟻，如圖1所示。

Mirjalili根據其習性提出了蟻獅優化算法[16]，步驟如下：

（1）螞蟻在尋找食物時，按照式（1）進行隨機游走：

其中，是螞蟻的位置，是計算螞蟻游走位置的累計和，是最大迭代次數，是一個代數式，如式（2）所示：

其中，是隨機游走步長，是一個隨機函數，如式（3）所示：

其中，是在[0，1]均勻分布的隨機數。

（2）在每次優化過程中，都會更新螞蟻的位置，但如果使用式（1）更新螞蟻位置，則無法保證螞蟻是在求解空間內進行隨機搜索，因此Mirjalili對式（1）進行了歸一化處理，如式（4）所示：

其中，是第維變量在第次迭代時的標準化位置，和分別是第維變量隨機游走步長的最小值和最大值，和是第維變量在第次迭代中隨機游走的最小值和最大值，如式（5）、式（6）所示：

其中，是第代中所有螞蟻的位置最小值，是第代中所有螞蟻位置的最大值，、是第代中第只和第只蟻獅的位置。

（3）蟻獅依靠陷阱捕捉螞蟻，這里的“陷阱”實際是螞蟻隨機游走范圍不斷縮小，如式（7）、式（8）所示：

其中，，是當前迭代次數，是最大迭代次數，是一個常數。

（4）蟻獅捕食成功后，將自己的位置轉移至捕食到螞蟻的位置，如式（9）所示：

其中，是第代中第只螞蟻的位置，是優化目標函數。

（5）在每次迭代過程中，都會將適應度值最好的蟻獅保存下來作為精英蟻獅，這會影響螞蟻的游走，需更新螞蟻的位置，如式（10）所示：

其中，是在第代中圍繞由輪盤賭方式選中的蟻獅游走的螞蟻位置，是第次迭代過程中圍繞適應度最優的精英蟻獅游走的螞蟻位置。

3? 改進型蟻獅優化算法（Improved ant lion optimization algorithm）

3.1? ?自適應邊界調節策略

在蟻獅優化算法的式（7）、式（8）中，螞蟻圍繞蟻獅游走的活動范圍僅受當前迭代次數和最大迭代次數的影響，這就導致在當前迭代次數下的所有螞蟻的游走范圍是一樣的，使螞蟻種群的多樣性下降，影響算法的尋優精度。對此，本文借鑒文獻[14]的自適應邊界策略，并在此基礎上在螞蟻的游走步長里融入了迭代次數的影響，提出改進自適應邊界調節策略，對式（7）、式（8）中的進行如下改進：

其中，是[0，1]內均勻分布的隨機數。

3.2? ?基于柯西變異的分段搜索策略

根據式（4）可知，螞蟻的位置由自身位置、蟻獅位置和所設陷阱位置決定，因此蟻獅算法具有較強的局部搜索能力，但易被局部極值點“迷惑”。為了提高算法跳出局部極值點的能力，并加快算法的收斂速度，對算法進行分段化處理。

算法搜索的前半段，借鑒粒子群算法中對每代最優解進行記憶保存的思想，對式（4）進行如下改進：

其中，是第只螞蟻在第代的位置，是第代螞蟻的最優位置，是第代中第只螞蟻的位置，和是第代中任意兩只螞蟻的位置，和是[0，1]內均勻分布的隨機數，是慣性權重，是學習因子，是螞蟻的種群大小，是螞蟻位置更新的步長函數，如式（15）所示：

算法的后半段，利用柯西變異“幫助”算法跳出局部極值點，縮短算法的尋優時間。柯西分布的概率密度函數如式（16）所示：

其中，是位置參數，是比例參數，當，時是標準柯西分布，如式（17）所示：

與高斯分布相比，標準柯西分布在原點處的波峰較矮，在兩端的曲線坡度較緩，取值范圍較廣，因此柯西變異的擾動能力更強。將柯西變異引入螞蟻位置更新公式中，可以以較短的時間來搜索相鄰區間，不受局部極值點的“迷惑”，提升算法的全局尋優性。因此，對式（4）進行如下改進：

其中，是標準柯西分布，是[0，1]內均勻分布的隨機數。算法搜索的前半段和后半段由概率決定，

，是最大迭代次數。

4? ?仿真實例（Simulation example）

4.1? ?乙酸乙酯生產過程間歇式反應釜描述

設乙酸乙酯間歇生產過程反應單元發生如下反應：。是乙醇和乙酸混合溶液，是乙酸乙酯，是乙醚。反應過程的微分方程描述[17]，如式（19）所示：

其中，、表示、的濃度。

反應過程的溫度是控制變量，控制目標函數是：

其中，和是數量級因子，用于統一量綱，，。代表反應過程中的溫度，可行域范圍為[600，960]。代表產出符合質量要求的乙酸乙酯所需時間。合格產物的濃度是0.47，因此，控制優化的性能指標是在處獲得最大值，并在時得出反應過程的反應溫度的控制軌跡。

4.2? ?仿真結果

為了驗證算法CSALO的有效性和可行性，將CSALO與蟻獅算法（ALO）、具有自適應邊界的蟻獅算法（SALO）和基于柯西變異分段搜索的蟻獅算法（CALO）在相同的條件下進行仿真對比分析。所用軟件信息：Windows 10、MATLAB R2017a。參數初始化：螞蟻和蟻獅種群大小均為40，最大迭代次數為500。運行5 次，得到目標函數值和產物的濃度值，如表1、表2所示。目標函數優化曲線和反應過程的反應溫度控制軌跡如圖2、圖3所示。

表1、表2的尋優結果表明，在目標函數值的仿真結果上，CSALO的平均值和最劣值最高，CALO的最優值最高，而SALO的最優值、最劣值和平均值均最差。綜合來看，CSALO算法的優化效果最好，SALO算法的優化效果最差，SALO和ALO算法有陷入局部最優的問題，CALO算法雖然克服了局部最優問題，但算法的穩定性較差。圖2、圖3表明，CSALO算法收斂速度最快，尋優所需時間最短。

由上述分析可知，CSALO的尋優精度和收斂速度明顯更優，算法求解的穩定性更好，表現出較強的全局搜索能力、局部開發能力和跳出局部極值能力。這主要是因為在搜索螞蟻位置時加入了柯西變異，有利于算法從局部極值中跳出，螞蟻游走過程中的自適應邊界和分段搜索策略有效擴大了螞蟻種群多樣性，提高了算法找到理論最優值的概率。

綜上所述，針對乙酸乙酯間歇反應單元的參數優化問題，CSALO算法在間歇型乙酸乙酯生產過程中的產物濃度、反應溫度和生產效率的平衡點的優化效果顯著。比起蟻獅算法，在工程應用上，CSALO算法仍然具有較好的尋優精度和收斂速度，在保證反應時間較短的情況下，有效提高了乙酸乙酯的濃度，降低了能耗，表現出較好的尋優性能。

5? ?結論（Conclusion）

反應釜是化工生產的第一工段，對產品的質量起到決定性作用。為了有效地提高反應釜產物的濃度和降低能耗，本文提出了基于蟻獅算法的優化方案，并針對蟻獅算法存在的問題提出了改進型蟻獅算法。仿真實驗表明，改進的蟻獅算法與蟻獅算法相比提高了乙酸乙酯的濃度，降低了能耗，優化效果更好。

就化工生產而言，除了反應釜，還有中和單元和精餾單元，而針對它們的優化方案還有待研究。

參考文獻（References）

[1] 陸寧云，王福利，高福榮，等.間歇過程的統計建模與在線監測[J].自動化學報，2006，32（3）：400-410.

[2] LEE J， LEE K S， LEE J H， et al. An online batch span minimization and quality control strategy for batch and semi-batch processes[J]. Control Engineering Practice， 2001， 9（8）：901-909.

[3] 賈立，曹魯明，邱銘森.基于建模誤差PDF形狀的間歇過程數據驅動模型[J].儀器儀表學報，2012，33（7）：1505-1512.

[4] 陳治綱，許超，邵惠鶴.間歇過程優化與先進控制綜述[J].化工自動化及儀表，2003，30（3）：1-6.

[5] 楊嵐.基于PSO的間歇型乙酸乙酯生產線的優化方法研究[D].杭州：浙江理工大學，2018.

[6] LUUS R. Optimal control of batch reactors by iterative dynamic programming[J]. Journal of Process Control， 1994， 4（4）：218-226.

[7] NGUANG S K， CHEN L Z， CHEN X D. Optimisation of fed-batch culture of hybridoma cells using genetic algorithms[J]. ISA Transactions， 2001， 40（4）：381-389.

[8] VENTER G， JAROSLAW S S. Particle Swarm Optimization[J]. AIAA Journal， 2003， 41（8）：129-132.

[9] 李宏光，劉驥鵬，黃靜雯.自適應變步長迭代動態規劃方法及其在間歇過程優化中的應用[J].控制與決策，2015，30（11）：

2048-2054.

[10] 趙博，袁希鋼，羅祎青.考慮環境影響的間歇過程多目標最優化設計[J].化工進展，2007，26（1）：113-118.

[11] 陳偉，賈立.間歇過程PSO-SQP混合優化算法研究[J].儀器儀表學報，2016，37（2）：339-347.

[12] SEYEDALI M， PRADEEP J， SHAHRZAD S. Multi-objective ant lion optimizer： A multi-objective optimization algorithm for solving engineering problems[J]. Applied Intelligence， 2017， 46（1）：79-95.

[13] 徐欽帥，何慶，魏康園.改進蟻獅算法的無線傳感器網絡覆蓋優化[J].傳感技術學報，2019，32（2）：18.

[14] 劉景森，霍宇，李煜.優選策略的自適應蟻獅優化算法[J].模式識別與人工智能，2020，33（2）：121-132.

[15] 何慶，林杰，徐航.混合柯西變異和均勻分布的蝗蟲優化算法[J/OL].控制與決策，2020，36（7）：1-10.https：//doi.org/10.13195/j.kzyjc.2019.1609.

[16] SEYEDALI M. The ant lion optimizer[J]. Advances in Engineering Software， 2015， 83（1）：80-98.

[17] ZHANG J. Modelling and multi-objective optimal control of batch processes using recurrent neuro-fuzzy networks[J]. International Journal of Automation and Computing， 2006， 20（1）：1-7.

作者簡介：

汪? 會（1995-），女，碩士生.研究領域：間歇型乙酸乙酯生產線的優化.

潘海鵬（1965-），男，碩士，教授.研究領域：系統建模和控制，智能檢測與控制，工業過程控制與綜合自動化.

張益波（1980-），男，博士，副教授.研究領域：非線性系統控制，多智能體控制.