基于WOA優(yōu)化模糊PID的設施智能水肥系統(tǒng)設計與實現(xiàn)

趙澤能，許敏界，華 珊，李雙偉，徐志福，韓愷源，陳貴才

（1.浙江科技學院生物與化學工程學院，浙江 杭州 310023；2.浙江省農(nóng)業(yè)科學院農(nóng)業(yè)裝備研究所，浙江 杭州 310021）

0 引 言

水肥一體化是建立在滴灌系統(tǒng)基礎上的現(xiàn)代農(nóng)業(yè)技術，旨在實現(xiàn)節(jié)水、節(jié)肥、高產(chǎn)、高效的目標。以肥調(diào)水、以水促肥為核心理念，在現(xiàn)代農(nóng)業(yè)科學領域受到廣泛關注，并對農(nóng)業(yè)的可持續(xù)發(fā)展產(chǎn)生著重要的影響[1,2]。在水肥一體化系統(tǒng)中，精確配肥至關重要。通過控制肥液濃度在最佳范圍內(nèi)，可促進作物根系的生長發(fā)育，并增強土壤中礦物質(zhì)的吸收能力[3]，例如不同水肥的調(diào)控對水稻生長特征表現(xiàn)就較為明顯[4]。隨著5G、物聯(lián)網(wǎng)、現(xiàn)代工業(yè)控制和無線通信等技術應用于農(nóng)業(yè)生產(chǎn)中，現(xiàn)代農(nóng)業(yè)逐漸向智慧農(nóng)業(yè)轉(zhuǎn)型[5]。

水肥一體化系統(tǒng)的時變性、非線性和滯后性等問題會對水肥灌溉的運行效率產(chǎn)生不利影響，因此制定高效的肥水灌溉控制策略至關重要。目前國內(nèi)外主要采用PID 以及模糊PID的控制策略實現(xiàn)對系統(tǒng)的精確控制[6,7]。王正等[8]提出一種改進Smith預估補償?shù)哪：齈ID 控制，引入Smith預估器并對其進行結(jié)構(gòu)優(yōu)化，解決被控對象模型不適配引起的決策失誤問題，該方法相較于傳統(tǒng)模糊PID具有更好的控制品質(zhì)和更短的響應時間。謝佩軍[9]等提出改進離散灰色預測的新型變論域模糊PID 控制方法，融合超前控制的改進型離散灰色預測，增加了系統(tǒng)的適應性，解決了大時滯問題，提高了智能灌溉系統(tǒng)的肥液控制效果。宋卓研[10]等利用POS 和BP 神經(jīng)網(wǎng)絡算法優(yōu)化PID 控制參數(shù)，引入預測補償和提前干預，有效解決了非線性、時變性和滯后性等問題。許景輝[11]等設計出基于人群搜索算法優(yōu)化PID的智能灌溉控制策略，提高了系統(tǒng)性能，并將人群搜索算法與多種算法進行對比，體現(xiàn)算法優(yōu)越性，但該方法仍舊停留在仿真階段。吳琦[12]等利用蟻群算法優(yōu)化變論域伸縮因子，根據(jù)系統(tǒng)環(huán)境變化實時改變伸縮因子，實現(xiàn)控制pH在7.0。以上研究為基于本地端的控制，朱德蘭[13]等將自整定模糊PID 算法嵌入遠程服務終端設計出遠程調(diào)節(jié)EC的水肥控制系統(tǒng)。綜上所述，智能算法優(yōu)化PID一定程度上解決了水肥系統(tǒng)時變性、自適應能力差等問題，但當前研究大多還處于試驗仿真階段，遠程控制一定程度上提高了工作效率，但與智能算法結(jié)合的研究相對較少。

本文結(jié)合智能算法與遠程控制技術，設計了基于WOA 優(yōu)化模糊PID的設施智能水肥系統(tǒng)。利用鯨魚優(yōu)化算法收斂速度快、自適應強等特點優(yōu)化模糊PID參數(shù)，結(jié)合農(nóng)業(yè)物聯(lián)網(wǎng)遠程監(jiān)測和控制技術，實現(xiàn)對水肥EC值的遠程精確調(diào)控。

1 水肥系統(tǒng)架構(gòu)設計

1.1 硬件結(jié)構(gòu)設計

水肥系統(tǒng)結(jié)構(gòu)圖如圖1所示，該水肥系統(tǒng)主要由恒壓供水模塊、母液補給模塊、出肥模塊、灌溉區(qū)以及控制中心組成。恒壓供水模塊主要由蓄水池、疊片過濾器、流量計、電磁閥、離心變頻泵組成。該模塊主要為系統(tǒng)提供純凈的水源以及動力。母液補給模塊由母液桶、酸液桶、疊片過濾器、電動球閥、流量計以及文丘里吸肥器組成。母液桶和酸液或母液桶中裝有攪拌機與液位傳感器，電動球閥通過RS-485 總線收發(fā)信號并智能調(diào)節(jié)閥門開度以控制流量。流量計通過RS-485 總線傳遞母液流量數(shù)據(jù)。文丘里吸肥器利用流速產(chǎn)生壓力差的原理吸收母液至混肥桶。出肥模塊由混肥桶、電磁閥、EC傳感器、pH 傳感器、變頻離心泵、流量計組成。混肥桶具有同時出肥和回流功能。EC傳感器和pH 傳感器分別檢測EC以及pH值。變頻離心泵為吸肥提供動力。控制中心主要由上位機、PLC控制單元組成。上位機負責將獲取的傳感器數(shù)據(jù)傳輸給遠程服務器并將反饋的決策值給PLC控制控制單元。PLC執(zhí)行指令完成相應動作，實現(xiàn)精準控制水肥EC值。

圖1 水肥系統(tǒng)結(jié)構(gòu)圖Fig.1 Structure diagram for water and fertilizer system

1.2 遠程控制系統(tǒng)

遠程控制系統(tǒng)主要分為遠程端和本地端兩部分，如圖2所示，本地端設備主要由觸摸屏、PLC、壓力傳感器、EC傳感器、流量計、pH 傳感器、電動球閥、電磁閥、攪拌機、離心泵以及本地端服務器組成。遠程端主要由云服務器、用戶終端以及開發(fā)者服務器組成。組態(tài)屏通過RS-485 總線與PLC 進行通訊，PLC控制電動球閥、電磁閥、攪拌機以及離心泵執(zhí)行指令并接受傳感器數(shù)據(jù)。本地端服務器通過RS-232 總線接收來自PLC 的信號并通過4G 信號與遠程端進行數(shù)據(jù)交互。用戶終端對接云服務器發(fā)送和接受數(shù)據(jù)便捷控制本地端設備運行，用戶終端主要包括PC 機以及移動端。開發(fā)者服務器通過SDK提供的API接口對接云服務器，將傳感器獲得的數(shù)據(jù)進行存儲以及算法處理，將接收的本地端EC傳感器數(shù)據(jù)通過基于WOA的模糊PID 算法處理。PLC 將遠程端算法輸出信號通過RS-485總線輸入給電動球閥，電動球閥智能調(diào)節(jié)開度控制注入管道的母液量進而實時控制系統(tǒng)出肥口溶液的EC值。

圖2 遠程控制結(jié)構(gòu)圖Fig.2 Remote control structure diagram

2 基于WOA優(yōu)化模糊PID的控制模型設計

2.1 水肥電導率分析模型

電導率的控制過程滿足質(zhì)量守恒原理，當系統(tǒng)趨于動態(tài)穩(wěn)定時，流入混肥桶的水和流入混肥桶中的母液的質(zhì)量和等于混肥桶中的肥液和流出混肥桶的肥液的質(zhì)量和，同時考慮流入混肥桶中的母液總量由電磁閥調(diào)節(jié)開度來實現(xiàn)，因此表達式為：

式中：V(t)為混肥桶中溶液體積，L；q1(t)為流入混肥桶的水流量，L/s；C1為流入混肥桶的水濃度，mol/L；q100%為電磁閥開度為100%時的流量，L/s；wi(t)為電磁閥開度信號；Cmi為流入混肥桶的i號母液濃度，mol/L；q2(t)為流出混肥桶的肥液流量，L/s；C2(t)為流出混肥桶的肥液濃度，mol/L。

在實際應用中，系統(tǒng)達到動態(tài)平衡，混肥桶中的V(t)恒定可視為常數(shù)VT，各流入母液吸肥管道負壓相等，為了便于分析，假設各母液管道控制電磁閥開度相等，母液濃度相等，又由于濃度與EC成正比關系，可將C1替換成流入混肥桶的電導率E1，將Cmi替換成流入各混肥桶母液電導率Em，將C2(t)替換成流出混肥桶的肥液電導率E2(t)，wi(t)替換成各路相等的電磁閥開度信號w(t)，式（1）可表示為：

從實際混肥特點看，可認為EC控制模型為典型一階滯后模型[14]，由于清水的EC值遠小于母液EC值可忽略，此時表達式可為：

將式（3）拉氏變換后得：

式中：E(s)和w(s)分別為E2(t)和w(t)拉氏變換后的表達方式。

由式（4）可知，控制響應特征為一階線性系統(tǒng)，通過試驗數(shù)據(jù)測得，當q100%=0.121 L/s，Em=5 mS/cm，VT=60 L，系統(tǒng)運行時出口壓力為0.16 MPa 時，相應的q2=1.1 L/s。實測滯后時間為0.57 s。將上述數(shù)據(jù)代入式（4）中的得到EC的近似傳遞函數(shù)為：

2.2 模糊PID控制

模糊PID 控制主要分為模糊控制和PID 控制兩部分，如圖3 所示。系統(tǒng)的誤差e和誤差變化率ec作為模糊控制器和PID控制器的輸入，模糊控制器通過將輸入信號模糊化、規(guī)則庫模糊推理以及解模糊化處理，輸出PID 參數(shù)修正量Kp1、Ki1、Kd1、并作為PID控制器的輸入，實時修正PID參數(shù)值。

圖3 模糊PID模型框圖Fig.3 Block diagram for fuzzy PID model

PID參數(shù)的調(diào)整公式為：

式中：Kp、Ki、Kd分別為修正后的比例、積分、微分系數(shù)，Kp0、Ki0、Kd0分別為PID修正前的比例、積分、微分參數(shù)。

模糊控制輸入的檢測值給定EC值偏差e、偏差變化率ec以及輸出的修正量Kp1、Ki1、Kd1通過模糊化分別轉(zhuǎn)換為模糊語言變量E、ΔE/Δt和KP、KI、KD。對應的模糊語言變量分別在其論域上定義7 個模糊子集{NB,NM,NS,ZO,PS,PM,PB}，表示負大、負中、負小、零、正小、正中、正大，量化論域都為{-6,-4,-2,0,2,4,6}。這些模糊子集采用三角形隸屬函數(shù)模糊化，Mamdani推理機制，重心解模糊化的方法。模糊控制規(guī)則如表1所示。

表1 模糊控制規(guī)則表Tab.1 Fuzzy control rules

2.3 WOA優(yōu)化模糊PID控制模型

2.3.1 WOA優(yōu)化算法

WOA 優(yōu)化算法是一種基于仿生學思想的優(yōu)化算法，由Mirjalili 等人于2016 年提出[15]。其靈感來源于座頭鯨的捕食行為，通過更新鯨魚位置尋找最優(yōu)解，具體的尋優(yōu)過程可以分為以下步驟：圍獵、氣泡網(wǎng)攻擊、隨機搜索機制。

（1）圍獵。鯨魚接近包圍獵物過程中發(fā)現(xiàn)最佳鯨魚位置值X*(t)，其他個體位置為Xi(t)，第i個鯨魚在t+1 代的搜索位置為Xi(t+1)，則其他鯨魚在第i代向最優(yōu)位置靠近的迭代公式為：

式中：D表示當前位置與最優(yōu)位置之間的距離；A和C為隨機參數(shù)；r為[0,1]的隨機值；t為當前迭代次數(shù)；Tmax為最大迭代次數(shù)。

（2）氣泡網(wǎng)攻擊。該階段模擬鯨魚發(fā)動氣泡捕食，分別有收縮包圍和螺旋式包圍兩種方式，以50%的概率更新位置，收縮包圍的數(shù)學模型和圍獵行為基本相似，只是將式（9）中A的范圍縮小至[-1,1]，其他公式保持不變。螺旋式包圍的數(shù)學模型如下：

式中：b表示對數(shù)螺旋形狀常數(shù)；l為[-1,1]隨機數(shù)。

（3）隨機搜索機制。收縮包圍的捕食行為的數(shù)學模型中∣A∣＜1，而當∣A∣＞1 時，說明鯨魚個體距最優(yōu)位置鯨魚較遠，此時當前鯨魚會選擇一個隨機鯨魚位置靠近，擴大搜索范圍，提高全局搜索能力，公式如下：

式中：Xrand(t)表示當前隨機鯨魚位置；D1表示當前隨機鯨魚位置與當前最優(yōu)鯨魚位置間距離。

2.3.2 WOA優(yōu)化模糊PID

模型如圖4 所示，通過WOA 優(yōu)化算法優(yōu)化模糊PID 的量化因子Ke、Kec以及比例因子Kp、Ki、Kd，在設定范圍內(nèi)找出滿足系統(tǒng)要求的最佳量化因子和比例因子，讓系統(tǒng)調(diào)控在具備專家經(jīng)驗的同時賦予其自學習能力，提高系統(tǒng)的自適應能力和魯棒性。

圖4 基于WOA優(yōu)化模糊PID控制框圖Fig.4 Block diagram for fuzzy PID control based on WOA optimization

WOA 優(yōu)化模糊PID 參數(shù)的流程圖如圖5 所示，具體步驟如下：

圖5 WOA優(yōu)化模糊PID參數(shù)流程Fig.5 The flow of WOA optimization for fuzzy PID Parameters

（1）首先初始化鯨魚種群。

（2）計算鯨魚個體的適應度，適應度函數(shù)選擇絕對積分性能指標（ITAE）作為適應度函數(shù)：

式中：JITAE表示絕對積分性能指標函數(shù)值；e(t)為系統(tǒng)誤差；t為積分時間。

（3）根據(jù)適應度值判斷是否滿足終止條件（滿足最小適應度值或達到最大迭代次數(shù)），若滿足將最佳位置參數(shù)代入系統(tǒng)退出循環(huán)，否則進入步驟（4）。

（4）利用鯨魚算法機制更新鯨魚位置并計算適應度值，將當前適應度值與歷史最佳適應度值進行比較，并替換成最優(yōu)位置參數(shù)。

（5）再次進入步驟（3）進行判斷，直到找到最佳位置參數(shù)。

圖6 為基于WOA 優(yōu)化模糊PID 的水肥系統(tǒng)控制流程圖，該流程步驟如下：

圖6 基于WOA優(yōu)化模糊PID的水肥系統(tǒng)控制流程圖Fig.6 Control flow chart for water and fertilizer system based on WOA optimization fuzzy PID

（1）用戶設定混肥EC值。

（2）獲得設定值與實際檢測值進行比較，獲得誤差e(t)以及誤差變化率de(t)。

（3）根據(jù)輸入的誤差以及誤差變化率，利用WOA 算法對PID參數(shù)進行優(yōu)化，獲得修正后的PID參數(shù)。

（4）將PID控制模塊的輸出信號傳遞至電動球閥調(diào)節(jié)閥門開度，控制水肥EC值。

（5）判斷肥液EC值是否到達指定范圍，若否，反饋當前水肥溶液EC值，再次回到步驟（2），直到滿足要求出肥。

3 試驗結(jié)果分析

3.1 控制模型仿真與性能對比

為了驗證基于WOA 優(yōu)化模糊PID 水肥系統(tǒng)的有效性，利用MATLAB中的Simulink對其進行建模仿真如圖7所示。

通過不斷的調(diào)試，選擇PID 參數(shù)Kp0=4，Ki0=0.483，Kd0=2.99，采用設定EC值為1 mS/cm的階躍信號作為輸入，分別對3種控制方法進行仿真實驗，結(jié)果如圖8所示。

圖8 控制模型仿真曲線Fig.8 Simulation curve for control models

由曲線可知，基于WOA 優(yōu)化模糊PID 與模糊PID、PID 控制相比響應速度更快，穩(wěn)定性更好，超調(diào)量明顯更小，穩(wěn)態(tài)精度較高，系統(tǒng)響應曲線平滑，魯棒性較好。各控制模型系統(tǒng)性能對比如表2所示。

表2 各控制模型系統(tǒng)性能對比Tab.2 Comparison of system performance for control models

由表2可知，基于WOA 優(yōu)化的模糊PID 系統(tǒng)較傳統(tǒng)的PID控制系統(tǒng)，超調(diào)量僅為PID 控制的2.7%，調(diào)節(jié)時間縮短了86.5%，系統(tǒng)的響應速度提高了5.4%，穩(wěn)態(tài)誤差降低了99.8%；較模糊PID 控制系統(tǒng)，超調(diào)量僅為模糊PID 控制的7.5%，調(diào)節(jié)時間縮短了85.4%，系統(tǒng)的響應速度保持基本不變，穩(wěn)態(tài)誤差降低了97.3%。

3.2 系統(tǒng)測試

系統(tǒng)裝置如圖9所示。于2023年3月在浙江省農(nóng)業(yè)科學院桑園科研創(chuàng)新實驗基地進行實地測試，母液選擇水溶性復合肥（N-P2O5-K2O，氮∶磷∶鉀=17∶17∶17）配置成的電導率為5 mS/cm的溶液。

圖9 水肥系統(tǒng)測試裝置Fig.9 Testing device for water and fertilizer system

試驗結(jié)果如表3 所示，由表3 可知基于WOA 的遠程模糊PID 控制的水肥系統(tǒng)較PID 控制和模糊PID 控制能夠較為準確的調(diào)節(jié)EC至設定值范圍附近，且EC波動幅度也明顯較小，隨著EC設定值的增加，EC波動幅度減小，系統(tǒng)的穩(wěn)態(tài)精度明顯提高，系統(tǒng)適應性較強，魯棒性較好。與傳統(tǒng)水肥灌溉方式相比，系統(tǒng)施肥用量減少了18%，節(jié)水12%。

表3 水肥EC值檢測對比表 mS/cmTab.3 Comparison table of EC values for water and fertilizer

4 結(jié) 論

本文將基于WOA 優(yōu)化模糊PID 的控制算法應用至水肥一體化系統(tǒng)中，配合遠程控制系統(tǒng)對水肥溶液EC值遠程調(diào)控，解決了水肥控制系統(tǒng)的非線性和時變性等問題，實現(xiàn)了水肥一體化系統(tǒng)的精確控制和穩(wěn)定運行，具有一定的實用價值。

（1）將農(nóng)業(yè)物聯(lián)網(wǎng)、自動化控制和無線通信等先進信息技術集成應用于水肥一體化系統(tǒng)中，實現(xiàn)了對水肥EC的精確控制與系統(tǒng)的高效運行。

（2）在滿足環(huán)境要求的情況下，根據(jù)水肥混肥原理建立簡化的水肥混肥數(shù)學模型，為EC的精確調(diào)控提供了理論基礎。MATLAB 仿真對比試驗結(jié)果表明該系統(tǒng)較傳統(tǒng)的PID控制模型，系統(tǒng)的超調(diào)量僅為PID 控制的2.7%，調(diào)節(jié)時間縮短了86.5%，穩(wěn)態(tài)誤差降低了99.8%。經(jīng)過系統(tǒng)測試可見，該水肥系統(tǒng)能夠較為準確的調(diào)節(jié)EC至設定值范圍附近，且EC波動幅度明顯較小，系統(tǒng)的穩(wěn)態(tài)精度明顯提高，系統(tǒng)適應性較強，魯棒性較好。與傳統(tǒng)水肥灌溉方式相比，系統(tǒng)施肥用量減少了18%，節(jié)水12%。