CN102545762B - 直驱式波浪发电系统的控制方法 - Google Patents
直驱式波浪发电系统的控制方法 Download PDFInfo
- Publication number
- CN102545762B CN102545762B CN201110429361.4A CN201110429361A CN102545762B CN 102545762 B CN102545762 B CN 102545762B CN 201110429361 A CN201110429361 A CN 201110429361A CN 102545762 B CN102545762 B CN 102545762B
- Authority
- CN
- China
- Prior art keywords
- msub
- mrow
- permanent magnet
- stator
- coordinate
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 16
- 238000010248 power generation Methods 0.000 title claims abstract description 13
- 230000009466 transformation Effects 0.000 claims abstract description 14
- 238000004804 winding Methods 0.000 claims description 10
- 230000001360 synchronised effect Effects 0.000 claims description 9
- 238000013016 damping Methods 0.000 claims description 7
- 238000006073 displacement reaction Methods 0.000 claims description 6
- 239000013535 sea water Substances 0.000 claims description 6
- 230000001131 transforming effect Effects 0.000 claims description 3
- 238000013178 mathematical model Methods 0.000 claims description 2
- 238000003306 harvesting Methods 0.000 claims 1
- 230000004044 response Effects 0.000 abstract description 3
- 238000011217 control strategy Methods 0.000 abstract description 2
- 230000004907 flux Effects 0.000 description 6
- 238000010586 diagram Methods 0.000 description 4
- 230000001788 irregular Effects 0.000 description 3
- 239000011159 matrix material Substances 0.000 description 3
- 230000005284 excitation Effects 0.000 description 2
- 230000001052 transient effect Effects 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000007599 discharging Methods 0.000 description 1
- 238000004870 electrical engineering Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P80/00—Climate change mitigation technologies for sector-wide applications
- Y02P80/10—Efficient use of energy, e.g. using compressed air or pressurized fluid as energy carrier
Landscapes
- Control Of Eletrric Generators (AREA)
Abstract
本发明涉及一种直驱式波浪发电系统的控制方法,包括以下步骤:1)通过位置传感器和速度传感器测量的到发电机转子的位置和速度;2)结合其他的参数计算得到当前的定子电流的控制目标值;3)同时测量发电机定子电流的瞬时值,通过坐标变换将a-b-c坐标下的定子电流变换到d-q坐标下;4)将定子电流的实测值和控制目标值比较后输入至比例积分控制器,比例积分控制器输出d-q坐标中定子电压的控制目标值,再通过坐标逆变换,将d-q坐标下的定子电压转换到a-b-c坐标;5)将a-b-c坐标下定子电压的目标值输入至直线永磁发电机中,对变换器进行控制,输出目标定子电压,使得定子电流与控制目标值相等。所述的直驱式波浪发电系统由阿基米德浮子和直线永磁发电机组成。该控制策略响应速度快,效率更高,能适用于全海况条件下。
Description
技术领域
本发明属于电气工程领域电力系统及其自动化研究方向,公布了一种能适用于全海况(规则波浪和不规则波浪,以及线性波浪和非线性波浪)条件下的直驱式波浪发电系统的最优控制方法。
背景技术
直驱式波浪发电系统(WEC)通常由浮子和直线永磁发电机组成,其结构简单,效率较高。基于阿基米德浮子(Archimedes Wave Swing,AWS)的波浪发电系统是第一个直接驱动式的波浪发电系统。对其进行最优控制目的是为了最大限度地从波浪中获取能量,例如通过将海水压入或排出AWS内的圆柱体,改变圆柱体内空气的体积和压力,使得AWS能与波浪共振达到控制目的,但这一方法响应速度慢,难以跟踪波浪的变化,并且控制的精度也不高;又例如在对线性发电机转子上的受力分析的基础上,将WEC等效于一个串联谐振电路,采用对反应力的控制来最大限度的获取波浪中的能量;再例如通过对线性永磁发电机的电气控制来调整AWS的阻尼达到控制的目的,这种方法使用起来比较方便和精确。但是这些方法都是在假定波浪是规则波浪的前提下得到的,而在实际工况下,作用于AWS的波浪不可能是规则的。亟需提供一种控制方法来解决问题。
发明内容
发明目的:本发明所要解决的技术问题是针对实际工况下作用于AWS的波浪是非规则的前提下提供一种响应速度快、易实现地控制方法更有效地获取波浪中的能量。
技术方案:本发明通过如下技术手段加以实现:
一种直驱式波浪发电系统的控制方法,包括以下步骤
1)通过位置传感器和速度传感器测量得到发电机转子的位置和速度;
2)浮子的水动力阻尼系数参数、发电机动子速度、永磁动子的磁链、波长等参数计算得到当前的定子电流的控制目标值;
3)同时测量发电机定子电流的瞬时值,通过坐标变换将a-b-c坐标下的定子电流通过派克变换变换到d-q坐标下;
4)将定子电流的实测值和控制目标值比较后输入至比例积分控制器,比例积分控制器输出d-q坐标中定子电压的控制目标值,再通过坐标逆变换,将d-q坐标下的定子电压转换回a-b-c坐标;
5)将a-b-c坐标下定子电压的目标值输入至发电机中,对变换器进行控制,输出 目标定子电压,使得定子电流与控制目标值相等。
所述的直驱式波浪发电系统由阿基米德浮子和直线永磁发电机组成。
所述的阿基米德浮子的数学模型为:
所述的直线永磁发电机在d-q坐标下的简化模型为:
定子电流为:
id=0
其中x为AWS浮子和直线永磁发电机动子的位移、y为AWS浮子和直线永磁发电机动子的速度、mtot为所有运动部件的质量总和、Fg为来自于直线永磁发电机作用于动子上的力、βwb为AWS的水动力阻尼系数、ks为AWS的弹性系数、Fw为波浪作用在阿基米德浮子上力的总和、Ls为同步电感、ω为直线永磁发电机的电角速度、Xs为同步电抗、R为单相绕组电阻、u为电压、i为电流、s为直线永磁发电机定子侧下标、d为d轴侧下标、q为q轴侧下标、ψPM为永磁动子的磁链。
所述阿基米德浮子与波浪共振,且来自于线性永磁发电机的阻力等于来自于海水的阻力。
所述的定子电流为q轴定子电流。
有益效果:本发明能应用于全海况条件下,通过控制使得AWS于波浪共振,更加有效地获取波浪中的能量,提高了WEC系统的效率。
附图说明
图1为本发明基于AWS波浪发电系统结构的示意图;
图2为本发明单质量块模型的示意图;
图3为本发明最优控制的控制框图。
具体实施方式
为了从波浪中获取最大能量,必须采用控制手段使得AWS与波浪共振。本控制方法相应速度快、能量获取效率高,且能适用于全海况条件下。
由于直线永磁发电机是由AWS直接驱动的,AWS和直线永磁发电机的动子可以用如图2所示的单质量块模型来描述。由于波浪发电系统通常是通过电力电子元件与电网相连,因此,从电网侧是感受不到AWS的惯性,所以,用于电力系统分析的波浪发电系统模型中,AWS的运动方程可以采用如下简化模型:
其中,x和y分别为AWS浮子和直线永磁发电机动子的位移和运动速度;mtot为所有运动部件的质量总和,包括AWS浮子、直线永磁发电机的动子、AWS上方的水,以及其他和浮子一同运动的部分;Fg来自于直线永磁发电机作用于动子上的力;βwb为AWS的水动力阻尼系数;ks为AWS的弹性系数,Fw为波浪作用在AWS浮子上力的总和,是整个波浪发电系统激励。
直线永磁发电机本质上是永磁发电机,与传统的旋转永磁发电机主要区别是发电机动子的运动规律不同。传统的旋转电机旋转速度和旋转方向保持不变,而直线永磁电机的动子(Translator)是进行来回往复运动,其速度的大小和方向是变化的。根据直线永磁发电机动子的运动特点,分为动子运动速度大于零和小于零两种情况来讨论直线永磁发电的模型。假定直线永磁发电机是对称的,并且当发电机的动子位移为零时,匝链在a相绕组上的磁链等于零。a-b-c坐标下的直线永磁发电机的模型可以写为。
1)y>0
采用与同步发电机相同的正方向定义,定子磁链方程可以写为:
ψs=-Liabc+ψPM_abc (3)
其中,ψs为定子磁链向量,ψs={ψa,ψb,ψc]T;iabc为定子电流向量,iabc=[ia,ib,ic]T;L为 电感矩阵;ψPM_abc为匝链在定子绕组上的动子磁链。L和ψPM_abc可以写为:
λ为直线永磁发电机的极距;ψPM为永磁动子的磁链;Lss为定子绕组的自感;M为定子绕组之间的互感,由于直线永磁发电的动子要比定子要长的多,通常在其运动过程中能够覆盖整个定子,所以Lss和M都保持恒定。
a-b-c坐标下定子电压方程为:
其中,us_abc为定子端电压向量;R为定子电阻矩阵,R=diag(Rs,Rs,Rs),Rs为定子绕组的电阻。
将式(3)代入式(4),定子电压方程可以写为:
参照传统旋转发电机的派克变换,选择如下的变换将直线永磁发电机在a-b-c坐标下的模型变换为d-q坐标下的模型。坐标变换如下:
fdq0=Dfabc (6)
其中, 需要指出,在稳态情况下,用于同步发电机建模d-q坐标的旋转速度和旋转方向保持不变,而这里的d-q坐标是固定在直线发电机的动子上,并且随着直线发电机的动子做来回往复运动,因此,d-q坐标运动速度的大小和方向是随时间变化的。
电压方程(5)左右两边同时乘以变换矩阵D可得:
其中:
将A和S代入式(7),可得:
其中,ω为直线永磁发电机的电角速度,ω=2πv/λ;Ls=Lss-M。
2)y<0
当直线永磁发电机的动子向相反的方向运动,运动到任意一个位置时,定子的速度大小,以及匝链在定子绕组上的磁链与动子正方向运动时是相等的,只是速度的方向相反,所以在定子上感应出的电压和电流的方向大小相同,而方向相反,所以在a-b-c坐标中具有如下关系:
fabc,v>0(x)=-fabc,v<0(x) (10)
另外由式(6)可以看出,坐标变换只是与位置x有关,所以式(6)在动子到达的任何位置都适用,式(10)左右两边同时乘以D可得:
fdq0,v>0(x)=-Dfabc,v<0(x)=-fdq0,v<0(x) (11)
将-uds,-uqs,-iqs,-ids和-ω代入到式(8)和式(9)中可得:
综合式(8),(9),(12)和(13),直线永磁发电机在d-q坐标中的模型可以写为:
其中,Xs为直线永磁发电机的同步电抗,Xs=|ω|Ls。
3、直线永磁发电机的实用模型
在电力系统动态分析计算中,由于发电机定子的动态过程相对比较短,通常忽略发电机的定子暂态,即令
由于是永磁发电机,定子直轴磁链,ψds和定子交轴磁链,ψqs可以写为:
ψds=Lsids+ψPM (17)
ψqs=Lsiqs (18)
由此可得:
则实用模型的定于电压万程变为
uds=-Rsids+Xsiqs (20)
uqs=-Rsiqs-Xsids+ωψPM (21)
功率方程为:
P=udsids+uqsiqs (22)
Q=uqsids-udsiqs (23)
在上述公式中:
x AWS浮子和直线永磁发电机动子的位移。
y AWS浮子和直线永磁发电机动子的速度。
mtot 所有运动部件的质量总和,包括AWS浮子、直线永磁发电机的动子、AWS上方的水,以及其他和浮子一同运动的部分。
Fg 来自于直线永磁发电机作用于动子上的力。
βwb AWS的水动力阻尼系数。
ks AWS的弹性系数。
Fw 波浪作用在AWS浮子上力的总和,是整个波浪发电系统激励。
Ls 同步电感。
ω直线永磁发电机的电角速度,ω=2πv/λ;λ表示波长。
Xs 同步电抗。
R 单相绕组电阻。
u 电压。
i 电流。
s 直线永磁发电机定子侧下标。
d,q d轴侧、q轴侧下标。
ψPM 永磁动子的磁链。
在实际工况下,考虑到波浪是非规则的:
则来自于直线永磁发电机作用于动子的力可以表示为:
其中Eq为直线永磁发电机的电动势:Eq=ωψPM,ω为动子的角速度:ω=2π/λ。
则(6)式可以表示为:
将(5)、(7)代入(2)可得:
上式中,等式右边为系统的输入量,y和x为状态变量;iq是控制变量;mtot,π,ψPM,λ,βwb和ks是常量,因此这是一个线性系统,可以分成两个子系统:
其中:x=x1+x2,v=v1+v2。
本发明是通过对线性永磁发电机的q轴定子电流的控制来达到最大限度获取波浪能量的目的的,为了达到这一控制目标要满足以下两个条件:
(1)AWS与波浪共振;
(2)来自于线性永磁发电机的阻力等于来自于海水的阻力。
对于子系统1,假设:
x1=X1sin(ω1t) (34)
v1=V1cos(ω1t) (35)
Eq1=E1cos(ω1t) (3
6)
则上述两个条件可以表示为:
通过对iq11和iq12的控制实现子系统1的阻尼控制,那么子系统1的最优控制器的电流设 定值为:
同样对第二个子系统可以得到:
则q轴定子电流值可以表示为:
为了减少定子中的能量损耗,d轴定子电流值为0,即:id_ref=0。
控制框图如图3所示。本控制策略通过对背靠背变换器中发电机侧变换器的控制来实现。通过位置传感器和速度传感器测量的到发电机转子的位置和速度,进而结合其他的参数计算得到当前的定子电流的控制目标值。同时测量发电机定子电流的瞬时值,通过坐标变换将a-b-c坐标下的定子电流变换到d-q坐标下。然后将定子电流的实测值和控制目标值比较后输入至比例积分控制器,比例积分控制器输出d-q坐标中定子电压的控制目标值,再通过坐标逆变换,将d-q坐标下的定子电压转换到a-b-c坐标。最后将a-b-c坐标下定子电压的目标值输入至直线永磁发电机中,对变换器进行控制,输出目标定子电压,进而实现对定子电流的控制,使得定子电流与控制目标值相等,获取最大的波浪能。
本发明能应用于全海况条件下,通过控制使得AWS于波浪共振,更加有效地获取波浪中的能量,提高了WEC系统的效率。
Claims (4)
1.一种直驱式波浪发电系统的控制方法,其特征在于:包括以下步骤:
1)通过位置传感器和速度传感器测量得到阿基米德浮子的位移和速度;
2)基于直驱式波浪发电系统的参数,包括运动部件的质量总和mtot、阿基米德浮子的弹性系数ks、阿基米德浮子的水动力阻尼系数参数βwb、永磁发电机动子的磁链ψPM和直线永磁发电机的极距λ,以及测量得到的阿基米德浮子的位移x和速度v,计算得到当前的定子电流的控制目标值为和id=0,这时阿基米德浮子与波浪共振,且来自于线性永磁发电机的阻力等于来自于海水的阻力;
3)测量发电机定子电流的瞬时值,通过坐标变换将三相a-b-c坐标下的定子电流通过派克变换变换到d-q坐标下;
4)将定子电流的实测值和控制目标值比较后输入至比例积分控制器,比例积分控制器输出d-q坐标中定子电压的控制目标值,再通过坐标逆变换,将d-q坐标下的定子电压通过转换回到a-b-c坐标;
5)将a-b-c坐标下定子电压的目标值输入至直线永磁发电机的端口,使得定子电流与控制目标值相等。
2.根据权利要求1所述的直驱式波浪发电系统的控制方法,其特征在于:所述的直驱式波浪发电系统由阿基米德浮子和直线永磁发电机组成。
3.根据权利要求1所述的直驱式波浪发电系统的控制方法,其特征在于:
所述的阿基米德浮子的数学模型为:
所述的直线永磁发电机在d-q坐标下的模型为:
其中,x、v、mtot、ks、βwb、ψPM、同上,Fg为来自于直线永磁发电机作用于动子上的力、Fw为波浪作用在阿基米德浮子上力的总和、Ls为同步电感、ω为直线永磁发电机的电角速度、Xs为同步电抗、R为单相绕组电阻、u为电压、i为电流、s为直线 永磁发电机定子侧下标、d为d轴侧下标、q为q轴侧下标,
阿基米德浮子模型和直线永磁发电机模型,共同构成了直驱式波浪发电系统的模型。
4.根据权利要求1所述的直驱式波浪发电系统的控制方法,其特征在于:所述控制是使得阿基米德浮子与波浪共振,且来自于线性永磁发电机的阻力等于来自于海水的阻力,从而获取最大的波浪能。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110429361.4A CN102545762B (zh) | 2011-12-20 | 2011-12-20 | 直驱式波浪发电系统的控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110429361.4A CN102545762B (zh) | 2011-12-20 | 2011-12-20 | 直驱式波浪发电系统的控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102545762A CN102545762A (zh) | 2012-07-04 |
CN102545762B true CN102545762B (zh) | 2015-06-24 |
Family
ID=46351863
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201110429361.4A Active CN102545762B (zh) | 2011-12-20 | 2011-12-20 | 直驱式波浪发电系统的控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102545762B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102966486A (zh) * | 2012-11-27 | 2013-03-13 | 天津大学 | 浮子式海浪发电装置的变特性曲线控制方法 |
JP6046182B2 (ja) * | 2015-02-27 | 2016-12-14 | ファナック株式会社 | 振動を抑制する機能を備えたモータ制御装置 |
CN106953566A (zh) * | 2017-05-03 | 2017-07-14 | 广东工业大学 | 一种基于模糊pid控制器的频率匹配的方法与装置 |
CN107104616B (zh) * | 2017-07-10 | 2019-07-19 | 广东工业大学 | 一种直驱式波浪发电系统分层鲁棒性控制方法及装置 |
CN107800340A (zh) * | 2017-11-23 | 2018-03-13 | 广东工业大学 | 一种波浪发电输出最大功率的控制方法及相关装置 |
CN113335467A (zh) * | 2021-05-28 | 2021-09-03 | 哈尔滨工业大学(深圳) | 一种基于波浪能和太阳能组合发电的海洋浮式平台 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101141110A (zh) * | 2007-09-07 | 2008-03-12 | 浙江大学 | 变速恒频双馈异步风力发电机转子电流无延时控制方法 |
CN101247104A (zh) * | 2008-03-27 | 2008-08-20 | 无锡市飞象电子有限公司 | 定子磁链定向的交流电机旋转速度和转子位置推测方法 |
CN101604953A (zh) * | 2009-07-13 | 2009-12-16 | 山西合创电力科技有限公司 | 永磁风力发电机矢量及直接转矩控制综合方法 |
CN101771380A (zh) * | 2010-02-11 | 2010-07-07 | 哈尔滨工业大学(威海) | 直接转矩控制的变频器空间矢量调制方法 |
CN101902192A (zh) * | 2010-07-15 | 2010-12-01 | 福州大学 | 混合式步进电机直接自控制方法 |
-
2011
- 2011-12-20 CN CN201110429361.4A patent/CN102545762B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101141110A (zh) * | 2007-09-07 | 2008-03-12 | 浙江大学 | 变速恒频双馈异步风力发电机转子电流无延时控制方法 |
CN101247104A (zh) * | 2008-03-27 | 2008-08-20 | 无锡市飞象电子有限公司 | 定子磁链定向的交流电机旋转速度和转子位置推测方法 |
CN101604953A (zh) * | 2009-07-13 | 2009-12-16 | 山西合创电力科技有限公司 | 永磁风力发电机矢量及直接转矩控制综合方法 |
CN101771380A (zh) * | 2010-02-11 | 2010-07-07 | 哈尔滨工业大学(威海) | 直接转矩控制的变频器空间矢量调制方法 |
CN101902192A (zh) * | 2010-07-15 | 2010-12-01 | 福州大学 | 混合式步进电机直接自控制方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102545762A (zh) | 2012-07-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102545762B (zh) | 直驱式波浪发电系统的控制方法 | |
CN102005843B (zh) | 无刷双馈风力发电机及无刷双馈风力发电系统控制方法 | |
CN103187919B (zh) | 一种永磁同步电机弱磁调速的系统和方法 | |
CN111914442B (zh) | 一种磁链可调永磁同步磁阻电机的建模方法 | |
CN103501146B (zh) | 无刷直流电机驱动系统的换相转矩脉动抑制方法及系统 | |
CN105024608B (zh) | 一种不平衡电网下基于矩阵变换器的pmsg比例积分谐振控制方法 | |
CN103872951A (zh) | 基于滑模磁链观测器的永磁同步电机转矩控制方法 | |
CN101286724A (zh) | 一种永磁同步电机直接转矩控制方法 | |
CN104300861B (zh) | 一种三相永磁同步电机的控制方法 | |
CN103338002A (zh) | 一种永磁同步电机永磁体磁通和交轴电感的辨识方法 | |
CN101764567A (zh) | 永磁同步风力发电机复合矢量控制方法 | |
CN101557193B (zh) | 非正弦反电动势表面式交流永磁电机矢量控制方法 | |
CN103346728A (zh) | 双馈风力发电机定转子电流传感器故障判断方法 | |
CN206640528U (zh) | 永磁同步电机无权值模型预测转矩控制系统 | |
CN109861605B (zh) | 一种永磁同步电机无差拍转矩预测控制方法 | |
CN108418485B (zh) | 一种隐极式混合励磁电机恒功率损耗模型预测控制方法 | |
CN202696533U (zh) | 一种变速永磁交流发电机系统 | |
CN102355192B (zh) | 双馈风力发电机无功功率的控制方法 | |
CN107070346B (zh) | 一种永磁磁阻型双转子电机的解耦控制方法 | |
CN103986381A (zh) | 海浪发电系统的微网构建最优化功率因数复合控制方法 | |
CN102032242A (zh) | 基于电参量信息融合的液压系统运行状态在线监测方法 | |
CN102522942A (zh) | 双馈风力发电机励磁控制方法 | |
CN102332861A (zh) | 双馈风力发电机有功功率的控制方法 | |
CN203377608U (zh) | 直驱永磁同步风力发电系统的电机侧变流器控制系统 | |
CN103066904B (zh) | 一种永磁风力发电机最大功率跟踪控制方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |