第55卷第1O期 电测与仪表 Electrical Measurement&Instrumentation VoI.55 No.10 2018年5月25日 May.25,2018 基于滑模控制的直驱式波浪发电系统MPPT控制策略冰 林巧梅 ,杨俊华 ,蔡浩然 ,杨金明 (1.广东工业大学自动化学院,广州510006;2.华南理工大学电力学院,广州510641) 摘要:为提高直驱型波浪发电系统的功率捕获和转换能力,通过分析系统运动方程,建立其数学模型。依据波 浪入射频率和幅值,构造系统最优功率输出条件,计算q轴电流参考值,应用滑模变结构控制算法实时跟踪参 考电流。针对发电系统中含有的未知干扰,引入鲁棒控制,使发电机状态亦能保持实时跟踪最优参考值。仿真 结果表明,不论是否有干扰存在,所提控制方案都能实现最优参考电流跟踪和最大功率捕获跟踪,系统鲁棒性 提高。 关键词:直驱式波浪发电系统;永磁直线发电机;最大功率点追踪;滑模鲁棒控制 中图分类号:TM619 文献标识码:A 文章编号:1001—1390(2018)10-0090-06 MPPT algorithm for direct—drive wave power generation system based on sliding mode control Lin Qiaomei ,Yang Junhua ,Cai Haoran ,Yang Jinming ( .Faculty of Automation,Guangdong University f oTechnology,Guangzhou 51 0006,China. 2.School of Electric Power,South China University of Technology,Guangzhou 510641,China) Abstract:In order to improve the power absorption and conversion capability of the direct—drive wave power genera- tion system,the motion equation of the system is analyzed and its mathematical model of the wave power system is es— tablished.According to the incident frequency and amplitude of the wave,the optimal power output condition of the system is constructed and the q-axis current reference value is calculated.The sliding mode control algorithm is imple- mented to track the reference current in real time.Aiming at the uncertain disturbance in wave power generation sys— tem,the robust control is added,SO that the generator state can still keep track of the optimal reference in real time. The simulation results show that the proposed robust control strategy can keep track the optimal reference current in re— al time and maintain the maximum energy absorption despite the presence of uncertain perturbation,which improves the robustness of the controlled system. Keywords:direct—drive wave power generation,PMLG,MPPT,sliding mode robust control 0 引 言 如何提高波浪发电系统的功率捕获和能量转换 效率,是波浪发电技术的关键,当波浪转换装置 (WEC)的频率与波浪频率相匹配时,可发生共振现 象,此时波浪力与速度同相,波浪转换装置可从波浪 中获取最大能量 。稳态条件下,具有固定质量和 清洁无污染的可再生能源中,风能、太阳能随机 性较大,而单位小时内波浪能变化幅度小、功率密度 大¨ 。据相关数据显示,波浪能平均密度达2 kW/m ~3 kW/m ,是风能密度(0.5 kW/m )的4~6倍,全 球海洋能发电量可达l0 Tw 2j。我国海岸线漫长,波 几何结构的波浪捕获装置的固有频率,不易改变,但 通过调节能量转换装置(FrO)控制发电机电磁力,改 变系统运动幅值和相位,可改变WEC频率,达到系 统共振 。 浪能蕴藏量高达0.15亿kW,可商业利用的装机容 量亦高达0.3亿kW~0.5亿kWl3 J,开发利用波浪能 可缓解我国能源匮乏的窘境。 基金项目:国家自然科学基金资助项目(513770265);广东省科技计划 为获取波浪发电系统的最大功率,系统须满足 一项目(2016B090912006);广东省自然科学基金项目(2015A030313487);广东 定的幅值和相位条件 J。幅值条件受时变海浪的 省教育部产学研合作专项资金(2013B090500089) 一影响,难以达到最优的幅值状态,所以幅值优化控制 90一 第55卷第l0期 电测与仪表 Electrical Measurement&Instrumentation VOI.55 No.10 May.25,2018 2018年5月25日 策略存在本质上的不足;通过调节WEC的谐振频率 实现相位优化的控制方式应用广泛,典型的控制策 略有反应式控制策略和锁存控制策略。文献[9]应 用反应式控制策略扩大波浪发电系统的有效频率捕 捉范围,收集更多的规则海浪频谱,控制简单且易于 工程实现,但针对非谐振式的波浪发电系统,该控制 策略依赖于对未来海浪信息的预测。文献[10]提出 了相位锁存法,通过合理锁存和释放浮子,强迫浮子 速度与入射波浪力同相,从而达到系统的最大波能 捕获,该方法的关键和难点在于如何确定WEC的最 佳锁存和释放时间。 ’ 实际上,波浪发电系统是一个非线性强耦合系 统,为此,文献[11]提出了非线性内模控制法,但其 控制效果取决于对系统模型估算的准确性,而估算 误差在所难免,因此难以完全抑制系统的不确定性 干扰。为实现波浪发电系统的最大波能捕获,文献 [12]提出了模糊控制策略,该模糊控制系统以WEC 的结构阻尼为控制参数,根据波浪幅值和系统累积 捕获的平均功率这两类历史数据,采用遗传算法估 算出系统阻尼,经模糊控制器调节得到期望的系统 阻尼参数,进而使系统达到最优运行。但该控制方 法仅考虑波浪幅值和系统累积捕获的平均功率作为 优化参数,其有效性有待进一步考证。借鉴风力发 电中最大功率跟踪控制策略¨ ,文献[14]采用了变 特性曲线控制策略,实验表明该策略不仅可提高功 率输出,且能避免电磁转矩波动。但变特性曲线中 线性上升区合适的斜率值难以确定,斜率太大或太 小都会使电磁转矩不能到达额定值。 为提高直驱型波浪发电系统MPPT控制策略的 效率,通过分析系统运动方程,建立系统数学模型。 以时变波浪周期和幅值为输入信号,在满足最优输 出功率条件下,在线计算q轴最佳参考电流。当系统 不含不确定性干扰时,采用等效滑模控制项实时跟 踪参考电流;当系统含不确定性干扰时,加入鲁棒控 制项,抑制干扰,使发电机状态实时保持跟踪最优参 考值,提高控制系统鲁棒性。 1 直驱式波浪转换装置的数学模型 1.1结构与原理 直驱式波浪转换装置如图1所示。 海浪能捕获装置主要由浮子构成,能量转换装 置主要由永磁直线发电机和变换器构成 。浮子随 海浪振荡,与之连接的动子直接驱动永磁同步发电 机做直线往复运动发电。 图1 直驱式波浪发电系统结构图 Fig.1 Structure diagram of the direct—drive wave power generation system 1.2动力模型 据牛顿第二定律,可得单一垂直自由度运动的 直驱式波浪发电装置的运动方程为: Mx(t)=Fh+F 一G (1) 式中 为浮子与动子等运动装置的总质量; 为动子位移; 为海水压强作用在浮子上的垂直分 力;G为运动装置的重力。 当海浪与浮子运动幅值都较小,则 可用如下 方程表示: Fh=Fd+F +Fh (2) 式中 表示人射波浪作用在浮子上的力;F 表 示辐射力;F 表示静浮力。 定动子位移 :0,F =0,则有: Fh =一pgSx=一Kx (3) 式中g为重力加速度;p为海水密度;S为浮子 的有效横截面积; =pgS。 针对规则的波浪,F 可表示为: F =一 互一 (4) 式中/4表示由F 产生的附加质量;B表示由F 产生的阻尼系数; 为动子速度。 波浪与浮子的特性直接影响人射波浪的作用效 果,即: Fd=A互 + + (5) 式中 为海浪往复运动的位移。 联立上述各方程,可得波浪发电动力学方程: (M+A) + +K 一F 一G=Fd (6) 式中 的平衡位置只受常数项重力G的影响, ~91— 第55卷第l0期 电测与仪表 ElectricalMeasurement&Instrumentation Vo1.55 No.1O May.25,2018 2018年5月25日 考虑重力后,选取新平衡位置 H ,简化式(6): (M+ )互+ +K 一F =Fd (7) 1.3永磁直线电机数学模型 为建立永磁直线发电机数学模型,作如下基本 假定:(1)动子及永磁体上均无阻尼绕组;(2)铁芯不 饱和;(3)不计磁滞和涡流损耗。 依照电动机惯例,定子磁链方程为: 十 (8) t,pq=Lqiq 式中 为永磁体转子磁链; 、 、L 、L 、 、 为d、q轴定子磁链、电感和电流。 永磁同步发电机的定子电压方程为: 砜+ + lUq +警+ 式中R为定子相电阻。 电机的反电磁力方程为: F :3盯 (1O) 式中下为极距。 由式(10)可知,永磁电机的反电磁力由定子d、q 轴电流分量决定。为减少损耗,提高控制效率,采用 i =0的矢量控制方法,定子电流完全决定反电磁力 的大小。此时电机产生的反电磁力可表示为: F = (11) 2直驱式波浪转换装置系统分析 当入射频率与波浪发电装置的固有频率匹配 时,发生共振,系统可捕捉到最大波浪能 引。从式 (7)可知,只有F 是可控制的力,若使系统发生共 振,可控制反电磁力使系统接近最优功率提取条件。 为实时调节反电磁力F ,构造其最优值为速度 与位移的线性组合: F =一R 一K (12) 式中R 表示产生有功功率的阻尼系数;K 表示 产生无功功率的弹性系数。 将式(12)代人式(7)中,有: ( +4) +(曰+B ) +(K+ ) =Fd (13) 此为典型的二阶线性微分方程,可表示为RLC 串联电路形式,等效模拟电路如图2所示,入射波浪 力 用等效电压源替代,动子速度 模拟电流,,反 电磁力F 模拟负载两端的电压。 一92一 M+A 1 I吖丫 一1 B J ~一 V、,\ 一 卜Ⅲ Rg {1 一 图2 直驱式波浪发电系统等效电路图 Fig.2 Equivalent circuit diagram of a direct—drive wave power generation system 据最大传输定理可知,当系统发生谐振,且负载 电阻R =B时,负载可获得最大功率,对应电机可从 波浪中捕获最大功率。实现最大功率传输条件为: R =B (14) ∞(M+A): ∞(+)=—— (15)(15) 式中to为等效电路频率。 由式(15)可得模拟电路的谐振频率 。: otO:^√M+A r ̄ +K s .,.(16) 而反电磁力中的弹性系数为: Ks= (M+A)一K (17) 将式(14)和式(17)代人式(11)可得到q轴参 考电流i : = { +[too(M+A)一 2)(18) 根据参考电流,设计直驱式波浪发电系统机侧 变换器最大波能捕获模型,如图3所示。 图3 最大波能捕获控制模型 Fig.3 Energy—maximizing absorption control mode 由图3可知,根据最优状态下计算的发电机q轴 参考电流i。 ,令d轴参考电流i =0,可到d、q轴 电流偏差信号,输入到滑模控制器,得到电压控制信 号 ,u ,再经过坐标变换得到两相静止电压u ,u。, 施加到SVPWM发生器产生驱动机侧变换器工作的 信号。由此可见,只要发电机的d、q轴电流能够实 第55卷第l0期 电测与仪表 Electrical neasurement&Instrumentation V01.55 NO.10 May.25,2018 2018年5月25日 时跟踪参考输入,就能保证能量转换装置(发电机) 获得最大波能。 3 直驱式波浪转换装置控制策略 滑模控制律可由等效控制项u 和切换鲁棒控制 项u 组成。在系统没有不确定性干扰时,等效控制 项u 。跟踪参考电流;当系统含有未知干扰时,切换到 鲁棒控制项 。 ,克服系统扰动,提高系统鲁棒性,保 证系统运行不离开滑模面。 取L :L。=L,由永磁同步发电机定子电压方程 转化可得: r, “d . . 一 一 (19) l【 一 一 一 定义d、q轴电流误差分别为 =0一i ,e。=t—q— q;控制量为 d,“q。 不考虑不确定性干扰情况下,设计d、q轴电流 滑模控制器切换函数: ed (20) tsq= q 取 = :0,有: =一 + =o l = =t—q— Uq+ ¨ + 0 (21) 则d、q轴电流滑模控制器的等效控制项分 别为: Jf u铂=u = +砜 (22) Uoq:“ : + + + 当系统函数含未知干扰d(t)时,为保证系统运 行状态不离开滑模面,设计切换鲁棒控制项: =ssgn(S) (23) 此时,d、q轴滑模控制器的输出由等效控制项和切 换控制项组成,则d,q轴电流滑模控制器输出分别为: fl =M d=Meq+Mq+ s = w=——: q+ d+s q+Rid+ssggn( 5)5) JI “ = +“ = + ¨ 'iT …gn s) (24) 式中 为严格的正实数。 d、q轴控制器输出满足Lyapunov稳定性,稳定 性证明如F,选择Lyapunov函数: : 1 s (25) [一 等 + + ·[_ 】≤0 (26) 满足Lyapunov稳定性条件。 vq = 口·s= · q q— 一+ + + J~ 。一 。 l 一 一 一 一争 一 ] 。 qESg(S) 一 + + + L[_ 一 ]≤0 (27) 满足Lyapunov稳定性条件。 4仿真实验与分析 ‘ 为验证所述系统结构和控制策略的有效性,基于 Maflab/Simulink环境,搭建小型直驱式波浪发电系统仿 真模型,仿真参数为,定子电感L=8.2 mH,定子电阻尺 =6.48 Q,永磁体磁链 =0.147 Wb,极距 =0.05 in, 运动部件质量M=1 000 kg,附加质量 =2(DO kg,阻尼 系数B:6 000 N,弹力系数K=8 000 N。 波浪力变化如图4所示,在时间5 S之前、5 S到 10 S之间、10 S之后,波浪频率分别为3 rad/s、6 rad/s、9 rad/s,波浪幅值分别为1 m、2 In、3m。 0 2 4 删 6 8 1O 12 14 16 18 2O t/s 图4 波浪力变化曲线图 Fig.4 Changing curves of the wave force 在不含未知干扰d(t)的情况下,滑模控制器的 输出只含有等效控制项u 从图5中可知,等效控制下,即便入射频率和幅 值不断变化,浮子速度仍能保持与波浪力同相,满足 最大波能捕获条件;没有任何控制时,浮子速度与波 浪力无法同相,亦无法捕获最大波能。 一93— 第55卷第1O期 电测与仪表 Electrical Measurement&Instrumentation Vo1.55 NO.10 May.25,2018 2018年5月25日 6 墨 4i n i至 4 石4 2 2 。 一一一胁_三~ / = c一 ~ 0 2 4 6 8 1O 12 14 16 18 2O /s (a)等效控制下浮子速度与波浪力波形 { 喜 t (b)无效控制下浮子速度与波浪力波形 图5 浮子速度与波浪力的仿真波形图 Fig.5 Simulation waveforms of buoy’S velocity and wave excitation force ,s (a)等效控制下电流跟踪曲线 t/s (b)无效控制下电流跟踪曲线 图6正常情况下电流跟踪曲线 Fig.6 Current tracking CuYVe under normal case 从图6中可知,等效控制下,发电机能实时跟踪 最佳参考电流,获得最大功率跟踪捕获;没有控制 时,不能实现最佳参考电流和最大波能跟踪。 由图7可知,等效控制下,WEC捕获到的平均功 率大概为6 kW,是没有控制时平均功率的60倍,由 此证明了所提控制策略的有效性。 当系统含有未知扰动d(t)时,为克服扰动,提高 控制系统的鲁棒性,滑模控制器输出由等效控制项 和切换鲁棒控制项构成,即 :u 。+ 。 ...——94...—— 一一~t~( r 士f, 三jW 二_… …卅㈨ Ⅲ…㈨ 0 . . ~川…卅川 …m川 ……川川 一 ; H . . .一, Ⅲ …… ¨ 附㈨ 一…m㈨ Ⅲ…㈨ Ⅲ…㈨㈨ T,●}t } } Re fer≤≤ 0/一< …:e、n ce cu rr/e、n t/\、 8 8 ㈣ 《 、一一0. 5 R al e me current l 0 2 4 6 8 lU l2 l4 16 l8 20 t (b)无效控制下电流跟踪曲线 图8扰动情况下电流跟踪曲线 Fig.8 Current tracking curve under disturbance conditions 图9比较了扰动情况下,加入鲁棒控制的平均功 率捕获情况。显而易见,没有加入鲁棒控制时,平均 功率为3 kwh,相比正常情况,功率吸收减少50%; 加入鲁棒控制后,平均功率为4.8 kWh,虽然相比正 常情况,功率吸收减少了20%,但相比没有加入鲁棒 控制,功率吸收增加了60%。 t}s 图9扰动情况下不同控制算法控制下平均功率捕获 曲线图 Fig.9 Average power absorption curve using different control algorithms under disturbance conditions 第55卷第10期 2018年5月25日 电测与仪表 Electrical Measurement&InstrlHnentation V0I.55 No.10 May.25,2018 5 结束语 为提高直驱型波浪发电系统的功率捕获效率和 转换能力,通过分析系统运动方程,结合波浪入射频 率和幅值,在满足最优功率输出条件下,在线计算g 轴电流参考值,采用等效滑模控制算法实时跟踪参 考电流。针对含未知干扰的波浪发电系统,加入鲁 棒控制,实现最优参考值实时跟踪。仿真结果表明, 所提控制方案在正常和扰动情况下,都可维持系统 最大功率捕获跟踪,能量转化效率提高。 参考文献 [1]MeAnhur S,Brekken T.K.A.Oceanwave power data generationfor grid integration studies[C]//2010 IEEE Power and Energy Society General Meeting.Minneapolis,MN.US:IEEE,July 2010:1-6. [2]J.Falnes.A review of wave—energy extraction[J].Marine Struc- tures,2007,20(4):185-201. [3]Koivusaari R.(2010).Available:http://www.aw—energy.com/ blog.html [4]Tedeschi E,Molinas M.Impact of control strategies on the rating of electric power take off for wave eneryg conversion[C]//2010 IEEE International Symposium on Industrila Electronics.Bari,Italy:IEEE, - 2010:2406-2411. [5]GreenhowM,Rosen JH,ReedM.Contorl strategiesforthe clam wave energy device[J].Applied ocean research,1984,6(4):197-206. [6]康庆,肖曦,聂赞相,等.直驱型海浪发电系统输出功率优化控 制策略[J].电力系统自动化,2013,37(3):24-29. [7]Shek J K H,Macpherson D E,Mueller M A,et a1.Reaction force control of a linear electrical generator for direct drive wave eneryg con— version[J].IET Renewable Power Generation,2007,1(1):17-24. [8]游亚戈,李伟,刘伟民,等.海洋能发电技术的发展现状与前景 [J].电力系统自动化,2010,34(14):1—12. [9]D.Valerio,P.Beirao,J.sa da Costa.Optimisation of wave energy extraction with the Archimedes Wave Swing[J].Ocean Engineering, 2007,34(17-18):2330-2344. [10]A.Babarit,G.Duclos,A.H.Cle ment.Comparison of latching control strategies for a heaving wave eneryg device in random sea[J]. Appliedocean research,2004,26(5):227-238. [1 1]A.F.de O.Fale ̄o.Modelling and control of oscillating—body wave energy conveners with hydraulic power take—off nad gas accumulator [J].Ocean Engineering,2007,(34):2021-2032. [12]M.A.Jama,A.Assi,Addy.W,et a1.Fuzzy Logic Based reactive Controller for Heaving Wave Energy Conveners[C]//2012 IEEE In International Conference on Renewable Energies ofr Developing Coun- tires.Proc:IEEE,2012:1-5. [13]Xia C L,Geng Q,Gu X,et a1.Input-output Feedback Linearization and Speed Control of a Surface Permanent magnet Synchronous Wind eGnerator with the Bostchopper Convener[J].IEEE Transactions on IndustirM Electornics,2012,59(9):3489-3500. [14]Xi Xiao,Xuanrui Huang,Qing Kang.A Hill Climbing Method— Based Maximum Power Point Tracking Strategy for Direct—Drive Wave Energy Conveners[J].IEEE Trnasactions on Industiral Electornics, 2016,63(1):257-267. [15]Drew B,Plummer A R,Sahinkaya M N.A review of wave eneryg eonve ̄er technology[J].Proceedings of the Institution of Mechanical Engineers,Part A:Journal of Power and Energy,2009,223(8): 887-9o2. 作者简介: 林巧梅(1991一),女,硕士研究生,研究方向为波浪 发电最大功率捕获。Email:571089407@qq.corn 杨俊华(1965一),男,博士,教授,主要研究方向为电机电器及其控 制,风力发电机组的设计与控制。Email:yly93@gdut.edu.cn 蔡浩然(1993一),男,硕士研究生,研究方向为交流电机控制、混沌控 制。Email:475536415@qq.corn 杨金明(1962一),男,博士,教授,研究方向为能源与动力系统中电力 电子及其控制技术。 收稿日期:2017-06—15;修回日期:2017 ̄8-22 (王艳丽编发) 一95—