PID參數對動力定位系統定位精度的影響

張嫦利， 王 磊， 李 博， 徐勝文

(上海交通大學 海洋工程國家重點實驗室， 上海 200240)

?

PID參數對動力定位系統定位精度的影響

張嫦利， 王 磊， 李 博， 徐勝文

(上海交通大學 海洋工程國家重點實驗室， 上海 200240)

海洋開發促進了動力定位技術在海洋結構物上的應用。動力定位系統能使海洋結構物保持在目標位置和艏向角上。PID控制器在動力定位系統中應用廣泛，其參數的選擇顯著影響動力定位系統的定位精度。為了提高動力定位系統定位精度，借助一個時域模擬軟件對PID控制器的參數即比例增益系數Kp、積分增益系數Kd和微分增益系數Ki進行了研究。結果顯示,動力定位系統的定位精度受到控制器參數選擇的影響。對于研究的動力定位駁船，通過試湊法找到滿足定位要求的一組控制參數后，通過控制變量法輕微增加Kp、Kd或Ki會使得定位精度增加。試湊法和控制變量法選取控制器參數的方法將對工程中操作動力定位系統具有一定的參考價值。

動力定位； PID控制； 時域模擬； 定位精度

0 引 言

隨著海洋資源需求的激增，諸如鉆井、鋪管、海洋調查、潛水等海上作業越來越廣泛[1]。海洋結構物經常需要保持在目標位置和特定的艏向角上以滿足作業要求[2]。現今用于海洋結構物定位的主流方式有錨泊定位、動力定位和錨泊輔助動力定位。其中，動力定位系統最基本的功能是能夠借助推進器使浮體保持在特定位置和艏向角上[3]。

動力定位系統主要包括位置測量系統、控制系統、推力系統和動力系統。其中，控制系統接收位置測量系統發出的位置和艏向角等信息，通過由外部環境計算出的擾動力和測得的位置，計算出保持目標位置和艏向角需要的總推力和總力矩，并借助推力分配策略進行分配。目前控制方法有PID控制、H∞魯棒控制、模糊控制、神經網絡控制、遺傳算法控制等[4-5]。

PID控制因為控制簡單、魯棒性好，是實際生產中最常用的一種二階線性控制方法[6]。PID控制方法由PID參數的選擇所確定，這對整個動力定位系統產生影響。衡量動力定位系統的標準主要有功率消耗和定位精度。動力定位船舶的定位精度受到PID參數的直接影響[7]。PID控制參數的選擇有多種方法，主要有試湊法、經驗數據法和擴充臨界比例度法[8]。

本文的研究對象是某裝有動力定位系統的駁船，主要目的是研究PID參數對動力定位系統定位精度的影響。本文中借助的軟件是上海交通大學構建的動力定位試驗模擬程序，用試湊法整定出可行的一組參數。一旦發現可行的參數，就可以通過控制變量法進行參數優化。對時域模擬結果進行統計分析，具體分析參數對定位精度的影響，同時給出了參數選擇的建議。

1 3自由度運動學和動力學方程

1.1 坐標系定義

1.2 3自由度運動學方程

動力定位系統一般只考慮水平面3自由度運動，即：縱蕩、橫蕩和艏搖。運動方程為：

(1)

1.3 3自由度動力學方程

由于動力定位系統只考慮浮體的低頻運動，因此只需要研究低頻運動的方程。低頻運動的3自由度運動學方程可寫為[9]：

(2)

2 PID控制原理

PID控制器是根據系統的誤差，利用比例、積分、微分控制計算得到控制信息的一種控制方法。比例(P)控制的輸出量和輸入的誤差信號成正比關系；積分(I)控制的輸出量與輸入誤差信號的積分成正比關系；微分(D)控制中，控制器的輸出和輸入信號的微分成正比關系[10-11]。

動力定位系統的控制系統要實現濾波和反饋控制兩大功能。動力定位系統不考慮高頻運動，因此需要濾去位置測量系統測到的高頻信號和噪聲。控制系統還需要計算出所需的總推力和力矩，并分配給各個推進器。這里的總推力和力矩指的是抵制橫蕩、縱蕩的力和抵制艏搖的力矩。PID控制器的目標是通過只補償觀測器估算的低頻運動，使得海洋結構物保持在特定的位置和艏向角上[12]。

本文借助的時域模擬程序是在Simulink平臺上建模。具體包括位置設定模塊、控制系統模塊、推力系統模塊、波浪載荷模塊、動力學方程模塊和狀態觀測器模塊[9]。控制系統模塊中選擇了PID 控制算法和Kalman濾波，使得控制系統形成了位置信息的后反饋系統加上風前饋系統。

動力定位系統的PID控制方法如下式所示[13]：

(3)

(4)

(5)

式中:Kp、Ki、Kd分別是比例、積分和微分增益系數;Δx、Δy、Δψ分別是大地坐標系中x、y方向和艏向上實際位置和目標位置的偏差;Tx、Ty、Mψ分別是控制系統計算出的x、y方向的推力和艏向的力矩;TWfx、TWfy、TWfψ分別是風前饋在x、y方向的力和艏向上的力矩。比例增益系數的作用是提供與位置偏差成比例的回復力；積分增益系數能夠抵消緩慢變化和靜態的平均力，微分增益系數提供的是阻尼力。

PID控制器的工作原理如圖2所示[14]。

3 數學模型

3.1 駁船參數

本文的研究對象是裝有動力定位系統的駁船，三維模型如圖3所示。

圖3 裝有動力定位系統的駁船的三維模型

駁船的主要尺度如下：長度LPP=150.00 m，寬度B=50.00 m,吃水T=10.00 m,重心高度(距基線)KG=10.00 m,排水量V=73 750 m3,橫搖慣性半徑Kxx=20.00 m,縱搖慣性半徑Kyy=39.00 m,艏搖慣性半徑Kzz=39.00 m。

此外，駁船采用6個全方位螺旋槳，螺旋槳直徑3.6 m，最大推力800 kN。具體布置如圖3所示。

3.2 海洋環境載荷

本文未考慮風載荷作用，由于此種情況下浪、流同向作用時為最惡劣的環境條件，故取浪、流同向作用作為計算的環境條件，其具體參數:環境力方向180°,有義波高HS=5.27 m,譜峰周期TP=13.4 s,譜峰因子γ=3.3,波浪譜Jonswap,流速VC=1 m/s-1。數值模擬中，環境載荷方向取180°，即在大地坐標系中沿著x軸的負方向。

4 PID控制器參數的整定

應用經典的試湊法對PID的參數進行反復調整，直至找到滿足定位要求的一組參數。動力定位船舶允許的位移是水深2%～4%[15]。駁船實際水深超過100 m。因此可取參數應保證駁船的位移在2～4 m。

設定目標艏向角30°，應用時域模擬程序對每組選定參數進行12 600 s的數值模擬，通過定位精度看參數是否可取。通過試湊法，發現Kp= 0.001 5、Kd= 0.06、Ki= 1×10-5這組參數滿足要求[9]。這組參數下，縱蕩、橫蕩、平均偏移半徑、艏搖的標準差分別是1.242、 0.619 2、1.387 8 m和0.968 5°。駁船的定位表現如圖4和圖5所示。在進行數據分析時，去除了首尾各1 800 s以保證數據的穩定性。

圖4 縱蕩和橫蕩的時歷曲線

圖5 艏搖的時歷曲線

5 數值模擬與分析

5.1 數值模擬

本文借助上海交通大學編寫的時域模擬程序，駁船的目標艏向角ψsy有30°和60°。

數值模擬的目的是比較不同參數下駁船定位效果。計算得到駁船的縱蕩、橫蕩、平均偏移半徑、艏搖的標準差。更高的定位精度代表著更優的定位性能，即將海洋結構物保持在目標位置和艏向上的能力。本文應用上述程序對每個情況進行12 600 s的數值模擬。

數值模擬中使用的PID參數如表1所示。用試湊法整定出可行參數后，運用控制變量法進行參數優化。對于每個目標艏向角，每次只改變Kp,Kd和Ki中的1個。對于每個參數選定3組數，因此每個艏向角對應9組參數。

5.2 計算結果和分析

對每組參數下的駁船進行數值模擬，直接得到位移(縱蕩、橫蕩、艏搖)的時歷。通過統計計算得到位移(縱蕩、橫蕩、平均偏移半徑、艏搖)的標準差，如表2所示。

從表2可以看出，當Kp,Kd和Ki中有任一個微量增加時，駁船的位移減少。這種現象對Kp最明顯。對每個參數以最小值為參照點，通過比較增量的百分比來研究。可以看出，在目標艏向角為30°時，當Kp增加時橫蕩的減小值達到19.65%；在目標艏向角為60°時，當Kp增加時橫蕩的減小值達到18.64%，平均偏移半徑在動力定位的定位精度中權重較大；當目標艏向角為30°時，分別微增Kp,Kd和Ki，平均偏移半徑的減小量分別達到12% , 5.7% 和 3.4%。對于目標艏向角60°來說，相應的值是15.3% , 3.2%和3.0%。位移的趨勢圖如圖6～8所示。

表1 數值模擬中控制器選擇的PID參數

表2 位移的標準差

圖6 Kp的變化對位移的影響

圖7 Ki的變化對位移的影響

(1) 由圖6可以看到,給驗證可行的比例系數Kp增加或減少一個微量值，定位精度有明顯變化。具體表現在，微量減小Kp會增大駁船的位移，微量增大Kp會減小駁船的位移。從PID控制原理上來說，一般情況下，增大Kp會加快系統的反應速度，對減小靜差有利[12]。靜差是控制系統準確性的重要指標，指的是被控制量的穩定值和給定值的差。在動力定位系統中，增大Kp有利于減小位置和艏向角的靜差，即位移。需要注意的是，這樣的趨勢只在一定范圍內成立，因為過大的比例系數會使系統產生較大的超調和振蕩，從而使穩定性變差。

圖8 Kd的變化對位移的影響

(2) 由圖7和圖8可以發現，給驗證可行的Ki和Kd增加或減少一個微量值，定位精度有變化，變化程度沒有改變Kp引起的大。具體表現在，微量減小Ki或Kd會增大駁船的位移，微量增大Ki或Kd會減小駁船的位移。一般情況下，增大Ki會加大積分作用，有利于增大穩定性，但會增大超調。增大Kd能加快系統響應，減小超調、增大穩定性，但會降低系統抗干擾的能力。因此，Ki和Kd對定位精度的有利影響趨勢也是在一定范圍內成立的。

(3) 由表2和圖6～8可知，改變Kp、Kd和Ki對橫蕩的影響最明顯。在目標艏向角30°時，分別增大Kp、Kd和Ki，橫蕩的減小量達到19.65%，9.54%和7.75%。在目標艏向角為60°時，分別增大Kp、Kd和Ki，橫蕩的減小量達到18.64%，1.53%和2.21%。這與推進器的布置方式有關。

6 結 論

分析了PID參數對某動力定位駁船的定位精度的影響，并提出了相關建議，得到了以下結論：

(1) 分別增大比例系數Kp、 微分系數Kd和積分系數Ki都能減小位移，提高駁船的定位精度。這個結論在一定范圍內保持成立。因此，在工程實際或模型

試驗中，通過試湊法得到滿足定位要求的一組參數后，即可微量增大某一PID控制系數來優化定位效果。

(2) 一定范圍內，增大比例系數Kp對改善動力定位系統的定位精度效果最為明顯。

(3) 本文研究的駁船的橫蕩受到PID參數的影響最大。工程實際或模型試驗中，對于具有相似布置方式的推進器的船舶，可以借鑒相關參考。

(4) 本文只采用了控制變量法對PID參數進行優化，工程實際中可以考慮同時增大PID參數來嘗試獲得更高的定位精度。

[1] Jorrit-Jan Serraris.Time domain analysis for DP simulations[C]∥AE, Hawaii,2009.

[2] 楊 歡，王 磊，李 欣.錨泊輔助動力定位研究與進展[J].實驗室研究與探索,2012,31(4):88-92.

[3] Asgeir J Soresen. A survey of dynamic positioning control systems[J]. Annual Reviews in Control, 2011,35:123-136.

[4] 趙之韻.船舶動力定位系統控制器設計[D].大連：大連海事大學,2009.

[5] 鄒 剛.基于PID 控制的錨泊輔助動力定位系統研究[D].上海：上海交通大學,2014.

[6] 阿依布·阿布都拉，王彩霞.PID控制在冶金和熱處理中的應用[J].Hot Working Technology, 2013，42(6): 226-229.

[7] Xu Shengwen, Wang Xuefeng, Wang Lei,etal. Dynamic positioning with roll-pitch motion control for a semi-submersible[C]∥OPE, Alaska, 2013:973 -978.

[8] 金 奇，鄧志杰.PID控制原理及參數整定方法[J].重慶工學院學報,2008,22(5):91-95.

[9] 李 博.動力定位系統的環境力前饋研究[D].上海：上海交通大學,2013.

[10] 鄒 剛，王 磊，劉翰林. 錨泊輔助動力定位系統艏向控制模式研究[J].海洋工程，2013,31(5):11-17.

[11] 楊巨慶，黃 健，段麗華. PID控制技術與應用[J]. 哈爾濱師范大學自然科學學報，2004,20(2)：76-79.

[12] Dong T Nguyen, Asgeir J S?nsen. Set point chasing for thruster-assisted position mooring[J]. IEEE Journal of Oceanic Engineering, 2009, 34:548-558.

[13] 劉翰林.旁靠多浮體動力定位系統數值模擬與研究[D] .上海：上海交通大學,2013.

[14] 周 利.半潛式鉆井平臺動力定位推力系統理論設計及模型試驗研究[D] .上海：上海交通大學,2009.

[15] Wang Lei, Sun Pan, Wang Liang. Research on time domain simulation of dynamic positioning for a deep water semi-submersible platform[J].Scientia Sinica Phys, Mech & Astron,2011, 41(2) :123-131.

Effects of Selection of PID Coefficients on Positioning Accuracy of DP System

ZHANGChang-li,WANGLei,LIBo,XUSheng-wen

(State Key Laboratory of Ocean Engineering, Shanghai Jiaotong University, Shanghai 200240, China)

Offshore exploration has promoted the development of dynamically positioned (DP) marine structures. With dynamic positioning system, ocean constructions can be kept at a desired position and heading. PID controller is widely used in the DP system (DPS). Furthermore, the positioning ability is essentially affected by the selection of PID coefficients. To enhance the positioning accuracy, the influence of the parameters of PID control, which include proportional gain, derivative gain and integral gain are investigated by a dynamic positioning simulation program. The obtained results indicate that the positioning accuracy is influenced strongly by these coefficients. As for the studied barge, a set of parameters is obtained by trial-and-error method, the positioning accuracy of the DPS increases as increasingKp,KiorKdslightly using control variable method. Additionally, a guidance that combines trial-and-error method and control variable method on the selection of the coefficients is given, the method may be useful for the operator of DPS.

dynamic positioning; PID control; time domain simulation; positioning accuracy

2014-05-07

國家自然科學基金項目(51179103)

張嫦利(1989-)，女，江蘇南通人，碩士生，主要從事動力定位研究。Tel.:18818212284; E-mail：zhangchangli@sjtu.edu.cn

王 磊(1971-)，男，安徽蚌埠人，副教授，主要從事動力定位系統的研究。Tel.:021-62932025; E-mail：wanglei@sjtu.edu.cn

P 751

A

1006-7167(2015)03-0008-05