， ， ，東亞(中國石油大學(xué)(華東) 化學(xué)工程學(xué)院， 山東 青島 266580)

隨著陸地石油資源的日漸枯竭，新能源還不能滿足人類社會的需求，各國紛紛把目光投向深海[1-2]。半潛式海洋鉆井平臺是深海油氣開采的重要設(shè)備，動力定位是其核心技術(shù)，當(dāng)海洋平臺受到風(fēng)、浪、流的干擾力和力矩的情況下，動力定位系統(tǒng)產(chǎn)生抵抗干擾的反力和反力矩，可維持平臺在期望位置或期望軌跡[3]。

由于海洋環(huán)境的復(fù)雜性，對海洋平臺動力定位系統(tǒng)操作精度的要求也嚴(yán)苛。為增強系統(tǒng)的魯棒性，近年來一些魯棒性較強的控制方法被提出，其中滑模變結(jié)構(gòu)控制具有較強的魯棒性的優(yōu)點得到廣泛運用[4]。文獻[5-6]首次將滑?？刂茟?yīng)用于動力定位系統(tǒng)，相比于傳統(tǒng)PID和擴展卡爾曼濾波器控制，取得了良好的控制效果。但是滑?？刂拼嬖谝欢ǖ亩墩駟栴}，影響了實際應(yīng)用。本研究基于文獻[7-8]，利用一種低通濾波器來實現(xiàn)虛擬控制輸入，并將其與滑?？刂扑惴ńY(jié)合，實現(xiàn)滑?？刂扑惴憫?yīng)迅速、對參數(shù)變化及擾動不靈敏[14]，有效削弱控制器輸出的抖振，改進現(xiàn)有海洋平臺滑模動力定位方法。

滑?？刂坡稍O(shè)計中采用等效切換控制，控制輸入由等效控制量和切換魯棒控制量組成，等效控制保證系統(tǒng)的狀態(tài)在滑模面上，切換控制保證系統(tǒng)的狀態(tài)在有限時間內(nèi)到達滑模面并沿滑模面漸進收斂到零點。在滑?？刂破骱推脚_之間加入濾波器，可有效抑制控制輸入的抖振，使得實際輸出的力矩信號更為平滑，從物理意義上降低平臺推進系統(tǒng)的能耗和磨損。仿真結(jié)果表明，相比于傳統(tǒng)的PD控制，基于低通濾波器的等效切換控制效果良好，響應(yīng)速度更快，控制精度更高，也更為穩(wěn)定。

1 問題描述

海洋平臺三自由度運動模型[9]為：

(1)

其中，質(zhì)量矩陣M和阻尼矩陣D分別如下：

(2)

(3)

其中，J(φ)代表慣性坐標(biāo)系和載體坐標(biāo)系速度向量轉(zhuǎn)換矩陣：

(4)

圖1 低通濾波器結(jié)構(gòu)原理圖Fig.1 The diagram of low pass filter

2 滑?？刂破髟O(shè)計

2.1 低通濾波器工作原理

基于低通濾波器的滑模控制方法將濾波器接入控制器與控制對象海洋平臺中間，低通濾波器對控制器輸出進行濾波，可抑制滑模控制結(jié)構(gòu)本身所帶來的抖振，使實際輸出的控制力矩更為平滑，在實際工程應(yīng)用中有重要價值。其控制系統(tǒng)結(jié)構(gòu)圖如圖1[14]所示。

圖1中，u(t)為虛擬控制輸入；τ(t)為實際控制輸入?；？刂破鞲鶕?jù)海洋平臺期望和實際位置的偏差進行計算，從而構(gòu)成完整的閉環(huán)系統(tǒng)。設(shè)計低通濾波器：

(5)

其中：λ>0,λ∈R。s代表微分函數(shù)。

由圖1和式(5)可得虛擬控制輸入和實際控制輸入關(guān)系為：

(6)

2.2 滑?？刂破髟O(shè)計及穩(wěn)定性證明

根據(jù)式(1)推導(dǎo)出：

(7)

將式(7)代入到式(1)中可得：

(8)

將式(8)代入濾波器導(dǎo)出公式(6)可得：

(9)

其中，定義

(10)

式(9)可寫為：

(11)

將(11)展開可以提取出：

(12)

海洋平臺的位置跟蹤誤差為實際位置向量與期望位置向量之差：

e=η-ηd。

(13)

設(shè)計滑模函數(shù)如下：

(14)

其中：Λi=diag(λi1,λi2,…,λin),λij>0,i=1,2,j=1,2…,n。

對式(14)進行求導(dǎo)可得：

(15)

將式(11)代入可得：

(16)

(17)

設(shè)定切換控制量：

(18)

其中，K∈R為控制器參數(shù)。由式(17)和式(18)可以得到虛擬控制輸入量為：

u=ueq+usw。

(19)

其中：ueq代表等效虛擬控制輸入量；usw代表切換虛擬控制輸入量。

為保證系統(tǒng)穩(wěn)定性，需要將設(shè)計的切換控制量代入Lyapunov函數(shù)進行穩(wěn)定性證明。

定理1對于如表達式(8)、(12)所示的系統(tǒng)，如果控制器設(shè)計為(19)，則跟蹤誤差e將會漸進收斂到0。

證明：選取如下Lyapunov函數(shù)：

(20)

對其進行求導(dǎo)：

(21)

3 仿真結(jié)果分析

用一座第六代半潛式海洋鉆井平臺作為實驗仿真對象?？刂茖ο髤?shù)為：長114.07 m，寬78.68 m，排水量51 624 t，吃水深度19 m[12]。根據(jù)海洋平臺三自由度運動模型方程，可得平臺質(zhì)量矩陣和阻尼矩陣如下[10]：

(22)

環(huán)境干擾采用高頻和低頻的復(fù)合形式[13]，形式如下：

(23)

3.1 定位仿真

圖2為平臺在海洋平面上三個自由度縱蕩、橫蕩、艏搖的定位響應(yīng)圖像，可以看出，平臺定位響應(yīng)迅速，到達指定位置后運行穩(wěn)定，除縱蕩方向，橫蕩和艏搖方向幾乎沒有超調(diào)。

3.2 跟蹤仿真

圖3為兩種控制方法下平臺在海洋平面上三個自由度縱蕩、橫蕩、艏搖的跟蹤響應(yīng)圖像，由圖像對比可知，基于低通濾波器的等效切換滑?？刂祈憫?yīng)速度略快于PD控制，并且控制精度要比PD控制更高，傳統(tǒng)PD控制在位置跟蹤中出現(xiàn)了平頂現(xiàn)象，無法獲得高精度的控制效果。圖4為基于低通濾波器等效切換滑?？刂坪蛡鹘y(tǒng)PD控制下的平臺在x和y二維平面內(nèi)軌跡運動圖像，可以看出，基于低通濾波器的滑?？刂菩Ч?，精度更高。

圖2 平臺縱蕩、橫蕩、艏搖定位圖像Fig.2 Surge,sway and yaw of positioning images

圖3 平臺縱蕩、橫蕩、艏搖跟蹤圖像Fig.3 Surge,sway and yaw of tracking images

圖4 平臺二維平面跟蹤圖像Fig.4 Tracking image in the xy plane

圖5～7分別為平臺在x、y、z方向濾波前后的控制輸入，經(jīng)過濾波器濾波后，平臺實際的輸出力矩中噪聲信號被過濾，實際控制輸入力矩曲線光滑穩(wěn)定，能有效減少推進裝置不必要的能量消耗，從而達到低耗高效，延長動力裝置的使用壽命。

圖5 平臺縱蕩控制輸入與實際力矩輸出對比Fig.5 Comparison of control input and practical torque in the direction of surge

圖6 平臺橫蕩控制輸入與實際力矩輸出對比Fig.6 Comparison of control input and practical torque in the direction of sway

圖7 平臺艏搖控制輸入與實際力矩輸出對比Fig.7 Comparison of control input and practical torque in the direction of yaw

4 結(jié) 論

針對半潛式海洋鉆井平臺的動力定位問題，在滑?？刂频幕A(chǔ)上，將等效切換控制方法與濾波器結(jié)合，設(shè)計出基于低通濾波器的等效切換滑?？刂破?。與傳統(tǒng)PD控制方法進行對比，仿真結(jié)果表明基于低通濾波器等效切換控制方法的控制品質(zhì)及響應(yīng)特性較為優(yōu)良，平臺在定位和跟蹤仿真中響應(yīng)迅速，運行穩(wěn)定；加入低通濾波器之后，控制器輸出抖振被有效抑制，控制系統(tǒng)的實際力矩輸出曲線更為平滑，減少海洋平臺推進器的能耗和推進系統(tǒng)內(nèi)部的磨損。

參考文獻：

[1]金秋，張國忠.世界海洋油氣開發(fā)現(xiàn)狀及前景展望[J].國際石油經(jīng)濟，2005，13(3):43-44.

JIN Qiu，ZHANG Guozhong.The current situation and prospect of offshore oil and gas development in the world[J].International Petroleum Economics，2005，13(3):43-44.

[2]張大剛.深海油田的開發(fā)—當(dāng)前國際應(yīng)用及發(fā)展趨勢[J].中國造船，2005，46(4):41-46.

ZHANG Dagang.The development of deep water oilfield:International application and development trends in current situation[J].China Shipbuilding，2005，46(4):41-46.

[3]趙志高，楊建民，王磊,等.動力定位系統(tǒng)發(fā)展?fàn)顩r及研究方法[J].海洋工程，2002，20(1):91-97.

ZHAO Zhigao，YANG Jianmin，WANG Lei,et al.The development and research of dynamic positioning system[J].Oceanographic Engineering，2002，20(1):91-97.

[4]謝彬.深水半潛式鉆井平臺設(shè)計與建造技術(shù)[M].北京：石油工業(yè)出版社，2013:196-197.

[5]TANNURI E A，AGOSTINHOA C，MORISHITA H M，et al.Dynamic positioning systems:An experimental analysis of sliding mode control[J].Control Engineering Practice，2010，18(10):1121-1132.

[6]BESSA W M，DUTRA M S，KREUZER E.Dynamic positioning of underwater robotic vehicles with thruster dynamics compensation[J].International Journal of Advanced Robotic Systems，2013，10:1-8.

[7]李嬌嬌，施小成，孟羽澤,等.基于非線性濾波器的船舶動力定位系統(tǒng)[J].應(yīng)用科技，2017，44(2):23-28.

LI Jiaojiao，SHI Xiaocheng，MENG Yuze，et al.The dynamic positioning system of ships based on nonlinear filter[J].Applied Technology，2017，44(2):23-28.

[8]BALCHEN J G，JENSSEN N A，MATHISEN E，et al.Dynamic positioning of floating vessles based on Kalman Filtering and optimal control [C]∥IEEE Conference on Decision and Control Including the Symposium on adaptive Processes.1980,107(3):852-864.

[9]邊信黔，付明玉，王元慧.船舶動力定位[M].北京：科學(xué)出版社，2011:51-53.

[10]賈欣樂，楊鹽生.船舶運動數(shù)學(xué)模型-機理建模與辨識建模[M].大連：大連海事大學(xué)出版社，1999:106-108.

[11] 鄭大鐘.線性系統(tǒng)理論[M].2版.北京：清華大學(xué)出版社，2002:218-230.

[12]WEI L，PANG Y J.Research on dynamic simulation of DP for a deep water semi-submersible platform[C]∥Proceedings of 2012 International Conference on Mechanical Engineering and Material Science.2012:642-645.

[13]楊楊.動力定位船舶非線性自適應(yīng)控制研究[D].大連：大連海事大學(xué)，2013:26-29.

[14]劉金琨.滑模變結(jié)構(gòu)控制MATLAB仿真[M].北京：清華大學(xué)出版社，2015:217-222.