基于转矩实时观测叶轮疲劳寿命的在线估计装置及方法
【专利摘要】基于转矩实时观测叶轮疲劳寿命的在线估计装置及方法，涉及双馈型风电机组的智能化健康状态评估技术。为了解决现有的技术不适用于评估叶轮疲劳寿命的问题。将软件模拟分析得到的叶轮实际转动惯量转换为叶轮等效转动惯量；由基于旋转坐标变换或基于电磁功率计算得到的发电机电磁转矩，获到发电机的输入转矩；结合转速传感器得到的发电机转速，得到叶轮气动转矩，通过雨流量法根据叶轮气动转矩得到叶轮载荷谱，根据叶轮载荷谱和叶轮参考载荷谱，评估叶轮疲劳寿命。本方法简单易行，不通过额外传感器，降低了成本，提高了可靠性。本发明适用于叶轮疲劳寿命的在线估计。
【专利说明】
基于转矩实时观测叶轮疲劳寿命的在线估计装置及方法
技术领域
[0001] 本发明涉及双馈型风电机组的智能化健康状态评估技术。
【背景技术】
[0002] 风电机组长期工作于野外环境并负荷运行，如何对机组状态进行监测与分析是提 高机组智能化水平的一个重要方面。其中，叶轮是机组的重要组成部分，机组工作时，风能 将在叶轮上产生极大的气动转矩，该转矩产生的机械能通过发电机转化为电能送往电网， 该气动转矩最终也会传递到塔顶，造成塔顶载荷。叶轮气动转矩反应了叶轮及塔顶载荷的 变化，而载荷的变化将引起材料的疲劳损伤，因此根据叶轮转矩监测可对叶轮疲劳寿命做 出评估。
[0003] 由于叶轮载荷对机组的可靠运行极为重要，现有的风电机组在设计时，都会对叶 轮载荷进行分析与优化，以满足使用年限标准。但很多分析是在实验假想情况下进行的，难 以反映真实风机运行状况。因此，研究叶轮的实时载荷是非常有意义的。叶轮的载荷分析包 括多方向上的各种应力及扭矩的分析，叶轮的气动转矩是叶轮最大的一种转矩。
[0004] 现有的对叶轮气动转矩载荷进行分析，一般有下列方法。
[0005] (1)通过软件建模方法进行模拟分析。通过Matlab或者专门的风电分析软件，如 Blade，可对不同风况(风速、风向）下叶片的受力和转矩进行分析。同时，结合风速分布曲 线，可近似得到载荷循环次数，从而验证叶轮疲劳载荷设计是否达到要求。此外，通过软件 模拟，还可以较为准确的计算不同风况下的最大载荷和极限载荷，为机组在极端风况下的 运行提供依据。但整体而言，这种方法不适用于机组疲劳寿命分析，机组微观地址的差异， 以及控制上差异，可能造成机组实际疲劳载荷谱与理论分析存在较大差异，难以用于评价 机组真实疲劳寿命。
[0006] (2)采用器件直接测量载荷。相对软件模拟方法，这种方法直接通过传感器测零部 件应力，能直观反映器件的应力和扭矩等信息，具有极强的实用性，便于得到部件载荷谱实 现疲劳寿命分析。其中，传感器通常选用阻抗随机械形变而变化的应变片，把应变片粘贴在 被测零部件表面，通过测量阻抗实现对应力的测量。但这种方法存在传感器自身疲劳寿命 问题，以及容易受到环境因素，例如温度的影响，同时存在安装困难等问题，增加了硬件成 本，可靠性不高。
【发明内容】
[0007] 本发明的目的是为了解决现有的技术不适用于评估叶轮疲劳寿命的问题，从而提 供基于转矩实时观测叶轮疲劳寿命的在线估计装置及方法。
[0008] 本发明所述第一种基于转矩实时观测叶轮疲劳寿命的在线估计装置，该装置包括 电压传感器、电流传感器、两相静止坐标变换器、第一转矩观测器、电机输入转矩计算器、叶 轮转矩计算器、叶轮载荷谱生成器、叶轮疲劳寿命估计器和转速传感器；
[0009] 所述电压传感器感应定子侧的三相电压和转子侧的三相电压，所述电流传感器感 应定子侧的三相电流和转子侧的三相电流；
[0010]电压传感器的输出端与两相静止坐标变换器的电压输入端相连，电流传感器的输 出端与两相静止坐标变换器的电流输入端相连；
[0011 ]两相静止坐标变换器的输出端与第一转矩观测器的电信号输入端相连；
[0012]第一转矩观测器的输出端与电机输入转矩计算器的电磁转矩输入端相连；
[0013 ]电机输入转矩计算器的输出端与叶轮转矩计算器的电机输入转矩输入端相连；
[0014] 转速传感器感应发电机转速，转速传感器的三个输出端分别与第一转矩观测器的 发电机转速输入端、电机输入转矩计算器的发电机转速输入端和叶轮转矩计算器的发电机 转速输入端相连；
[0015] 叶轮转矩计算器的输出端与叶轮载荷谱生成器的输入端相连；
[0016] 叶轮载荷谱生成器的输出端与叶轮疲劳寿命估计器的输入端相连。
[0017] 本发明所述第二种基于转矩实时观测叶轮疲劳寿命的在线估计装置，该装置包括 电压传感器、电流传感器、旋转坐标变换器、锁相环、光电编码器、积分器、减法器、第二矩阵 观测器、电机输入转矩计算器、叶轮转矩计算器、叶轮载荷谱生成器和叶轮疲劳寿命估计 器；
[0018] 所述电压传感器感应定子侧的三相电压，所述电流传感器感应定子侧的三相电流 和转子侧的三相电流；
[0019 ]电流传感器的输出端与旋转坐标变换器的电流输入端相连；
[0020] 电压传感器的定子三相电压输出端与锁相环的输入端相连；
[0021] 锁相环的定子电压旋转角度输出端分别与旋转坐标变换器的定子电压旋转角度 输入端和减法器的定子电压旋转角度输入端相连；
[0022 ]光电编码器米集电机转子机械转速，光电编码器的转子机械转速输出端与积分器 的输入端相连；
[0023] 积分器的转子机械角度输出端与减法器的转子机械角度输入端相连；
[0024] 减法器的转子同步旋转角度输出端与旋转坐标变换器的转子同步旋转角度输入 端相连；
[0025] 旋转坐标变换器输出端与第二转矩观测器的电信号输入端相连；
[0026] 第二转矩观测器的输出端与电机输入转矩计算器的电磁转矩输入端相连；
[0027] 电机输入转矩计算器的输出端与叶轮转矩计算器的电机输入转矩输入端相连；
[0028] 转速传感器感应发电机转速，转速传感器的两个输出端分别与电机输入转矩计算 器的发电机转速输入端和叶轮转矩计算器的发电机转速输入端相连；
[0029] 叶轮转矩计算器的输出端与叶轮载荷谱生成器的输入端相连；
[0030] 叶轮载荷谱生成器的输出端与叶轮疲劳寿命估计器的输入端相连。
[0031] 根据第一种基于转矩实时观测叶轮疲劳寿命的在线估计装置的在线估计方法，该 方法包括以下步骤：
[0032] 步骤一、电压传感器感应定子侧的三相电压和转子侧的三相电压，电流传感器感 应定子侧的三相电流和转子侧的三相电流;转速传感器感应发电机转速《 g;
[0033]步骤二、两相静止坐标变换器将定子侧的三相电压转换为静止坐标系下的定子侧 电压Uscc和定子侧电压Usf!，将定子侧的三相电流转换为静止坐标系下的定子侧电流iscc和定 子侧电流isP，将转子侧的三相电压转换为静止坐标系下的转子侧电压Ura和转子侧电压Urf!， 将转子侧的三相电流转换为静止坐标系下的转子侧电流和转子侧电流ire;
[0034] 步骤三、第一转矩观测器结合发电机转速C0g，将定子侧电压Usa、定子侧电压Use、定 子侧电流iSCI、定子侧电流i sp、转子侧电压Ura、转子侧电压Urf!、转子侧电流ira和转子侧电流 irf!转换为发电机电磁转矩Tan;
[0035] 步骤四、电机输入转矩计算器结合发电机转速C0g，将发电机电磁转矩转换为发 电机输入转矩Tg;
[0036] 步骤五、叶轮转矩计算器结合发电机转速C0g，将发电机输入转矩Tg转换为叶轮气 动转矩T a;
[0037] 步骤六、叶轮载荷谱生成器通过雨流量法根据叶轮气动转矩Ta得到叶轮载荷谱；
[0038] 步骤七、将步骤六得到的叶轮载荷谱和预设叶轮参考载荷谱进行比较，得到比较 结果，该结果作为叶轮疲劳寿命评估值。
[0039] 根据第二种基于转矩实时观测叶轮疲劳寿命的在线估计装置的在线估计方法，该 方法包括以下步骤：
[0040] 步骤一、电压传感器感应定子侧的三相电压，电流传感器感应定子侧的三相电流 和转子侧的三相电流;转速传感器感应发电机转速 wg;光电编码器米集电机转子机械转速 CO r ;
[0041] 步骤二、锁相环将定子侧的三相电压转换为定子电压旋转角度0S;
[0042] 步骤三、积分器将电机转子机械转速c〇r转换为发电机转子机械角度；
[0043] 步骤四、减法器结合定子电压旋转角度0S，将发电机转子机械角度转换为转子同 步旋转角度I;
[0044] 步骤五、旋转坐标变换器结合定子电压旋转角度0S和转子同步旋转角度0r，分别将 定子侧的三相电流转换为旋转坐标系下的定子侧电流ids和定子侧电流i qs，将转子侧的三 相电流转换为旋转坐标系下的转子侧电流idr和转子侧电流iqr;
[0045] 步骤六、第二转矩观测器将定子侧电流ids、定子侧电流iqs、转子侧电流idr和转子 侧电流转换为发电机电磁转矩1^;
[0046] 步骤七、电机输入转矩计算器结合发电机转速c〇g，将发电机电磁转矩Tem转换为发 电机输入转矩T g;
[0047] 步骤八、叶轮转矩计算器结合发电机转速《 g，将发电机输入转矩Tg转换为叶轮气 动转矩Ta;
[0048] 步骤九、叶轮载荷谱生成器通过雨流量法根据叶轮气动转矩Ta得到叶轮载荷谱；
[0049] 步骤十、将步骤九得到的叶轮载荷谱和预设叶轮参考载荷谱进行比较，得到比较 结果，该结果作为叶轮疲劳寿命评估值。
[0050] 本发明的原理为：
[0051] 风机，齿轮及电机之间的通过轴承进行连接，忽略轴系的旋转阻尼和弹性作用，并 将风力机和齿轮箱的转动惯量折算到电机侧，可得到如下运动方程，
(I)
[0053]其中Tg为发电机输入转矩，Ta为叶轮气动转矩，J%为叶轮等效转动惯量，为发 电机转速。
[0054] 可见，如果能够获得发电机的输入转矩、发电机转速以及叶轮等效转动惯量，则气 动转矩可以通过计算方法得到。发电机转速是风电机组正常监测参数，获取较为容易。等效 转动惯量可根据下列方法得到：
[0055] J，g = Ja/ng ⑵
[0056] 即叶轮等效转动惯量为叶轮实际转动惯量Ja除以齿轮变速箱的变速比ng。叶
[0057] 轮实际转动惯量可根据现有技术的有限元分析软件得到。
[0058]发电机输入转矩可通过在线测量方法得到。类似公式（1)，忽略摩擦阻尼影响，电 机运动方程有，
(3)
[0060] 其中Jg为发电机转动惯量，可由厂家提供也可由现有技术的常规实验获得，p为发 电机极对数。为发电机电磁转矩，可通过直接计算获得。
[0061] 对于风电场大量采用的双馈型发电机，其电磁转矩可在同步旋转坐标系上，由定 转子电流计算转矩。其转矩方程为
[0062] Tem = pLm(iqsidr-idsiqr) (4)
[0063] 其中，p为发电机极对数，Lm为发电机互感参数，ids和iqs分别为定子电流在定子侧 电压旋转坐标系下的d轴分量和q轴分量，idr和iqr分别为转子电流在转子侧旋转坐标系下 的d轴分量和q轴分量。为获得同步旋转坐标的旋转角，需要通过锁相环获得定子电压旋转 角度，以及通过转速传感器获得发电机转速。
[0064]此方法需要坐标变换，同时需要准确知道发电机互感参数，因此本发明还提出采 用电磁功率的方法计算发电机电磁转矩，该方法检测精度高、可靠性高。
[0065]对于双馈电机，其机械功率Pm与定子电磁功率匕和转子的电磁功率Pr之和应平衡， 即Pm=Ps+Pr，再根据功率与转矩关系可得发电机电磁转矩，具体计算方法为，
[0067 ] Pr-irctUrct-^ lrJ3UrJ3
[0068] Ps = isaUsa+isPUsP (5)
[0069] 其中ira，irfi，ura，u rf!及isa，isfi，usa，usf!分别是转子侧电流在两相静止坐标系下的a 轴分量和0轴分量、转子侧电压在两相静止坐标系下的^^由分量和0轴分量、定子侧电流在两 相静止坐标系下的a轴分量和0轴分量、定子侧电压在两相静止坐标系下的a轴分量和0轴分 量。得到叶轮气动转矩后，可对其载荷循环进行统计。本发明采用现有技术的雨流量法实现 统计。确定每个雨滴循环的起点和终点，可得到该载荷循环的幅值，即载荷变化幅值。风电 机组现场运行时，在线统计更新不同幅值下载荷的循环次数，这样即可得到叶轮气动转矩 的载荷谱。
[0070] 将得到的叶轮载荷谱与根据叶轮理论及实验获得的参考载荷谱进行比较，采用现 有技术的Miner累积损伤理论评估方法，最终可估算叶轮疲劳寿命。
[0071] 本发明直接根据监测发电机组运行时测量的电气参数，通过发电机电磁转矩与叶 轮气动转矩之间的关系，得到叶轮气动转矩，并通过雨流量法对叶轮载荷循环进行记录，以 实现叶轮疲劳寿命在线估计，实现机组可智能化运行。通过计算发电机输入转矩，逆推叶轮 气动转矩，实现不直接测量叶轮转矩值，而通过等效计算的方法得到叶轮转矩。
[0072] 有益效果:与现有技术的常规方法不同的是，本方法不直接测量叶轮转矩，而是通 过等效计算的方法，在已知叶轮和电机的转动惯量基础上，从发电机的电磁转矩直接逆推 得到叶轮转矩，从而避免了直接测量转矩时应变片安装困难，易受环境影响问题，同时降低 了成本，提高了可靠性。获得叶轮转矩后，通过雨流量法的载荷循环次数统计，得到载荷谱， 从而实现叶轮疲劳寿命在线估计。针对发电机的电磁转矩计算问题，提出了采用了基于旋 转坐标变换和基于电磁功率计算的两种方法，实现了冗余观测。本方法简单易行，不通过额 外传感器，电压传感器、电流传感器和转速传感器都是对风电机组进行常规检测用的器件， 本发明只根据风电机组常规运行时监测的数据，通过风电机组模型实现叶轮气动转矩测 量，提高了风组的自检自查能力和装置的智能化水平。
[0073] 本发明适用于叶轮疲劳寿命的在线估计。
【附图说明】
[0074] 图1是【具体实施方式】三所述的基于转矩实时观测叶轮疲劳寿命的在线估计方法的 简化框图；
[0075] 图2是【具体实施方式】一中的测量发电机电磁转矩的原理示意图；
[0076] 图3是【具体实施方式】一中的基于转矩实时观测叶轮疲劳寿命的在线估计装置的结 构示意图；
[0077] 图4是【具体实施方式】二中的基于转矩实时观测叶轮疲劳寿命的在线估计装置的结 构示意图；
[0078] 图5是【具体实施方式】四中的雨流量法进行载荷循环统计在示意图；
[0079] 图6是【具体实施方式】四中的发电机电磁转矩实际值与计算值对比图；
[0080] 图7是【具体实施方式】四中的发电机输入转矩实际值与计算值对比图。
【具体实施方式】
【具体实施方式】 [0081] 一：结合图1至图3具体说明本实施方式，本实施方式所述的基于转 矩实时观测叶轮疲劳寿命的在线估计装置，该装置包括电压传感器1、电流传感器2、两相静 止坐标变换器3、第一转矩观测器9、电机输入转矩计算器11、叶轮转矩计算器12、叶轮载荷 谱生成器13、叶轮疲劳寿命估计器14和转速传感器15;
[0082] 所述电压传感器1感应定子侧的三相电压和转子侧的三相电压，所述电流传感器2 感应定子侧的三相电流和转子侧的三相电流；
[0083]电压传感器1的输出端与两相静止坐标变换器3的电压输入端相连，电流传感器2 的输出端与两相静止坐标变换器3的电流输入端相连；
[0084] 两相静止坐标变换器3的输出端与第一转矩观测器9的电信号输入端相连；
[0085] 第一转矩观测器9的输出端与电机输入转矩计算器11的电磁转矩输入端相连；
[0086] 电机输入转矩计算器11的输出端与叶轮转矩计算器12的电机输入转矩输入端相 连；
[0087] 转速传感器15感应发电机转速，转速传感器15的三个输出端分别与第一转矩观测 器9的发电机转速输入端、电机输入转矩计算器11的发电机转速输入端和叶轮转矩计算器 12的发电机转速输入端相连；
[0088] 叶轮转矩计算器12的输出端与叶轮载荷谱生成器13的输入端相连；
[0089] 叶轮载荷谱生成器13的输出端与叶轮疲劳寿命估计器14的输入端相连。
[0090] 图2是测量发电机电磁转矩的结构示意图，只需要获得静止坐标系a-0轴上转子侧 电流与电压值和定子侧电流与电压值，即：^、;[1^、111^、111^和;[；^、;[；^、11；^、11；^，以此按公式(5) 计算功率和发电机电磁转矩即可。因两相静止坐标变换器输出电压包含大量开关谐波，因 此Ura、Urf需要进行低通滤波处理。图2中，abc/邱为两相静止坐标变换器，DFIG为双馈感应 电机，U sa、Usb、Usc为定子侧的三相电压，Ura、Urb、Urc为转子侧的三相电压，isa、isb、isc为定子 侦_三相电流，ira、irb、ire为转子侧的三相电流。
【具体实施方式】 [0091] 二：结合图2和图4具体说明本实施方式，本实施方式所述的基于转 矩实时观测叶轮疲劳寿命的在线估计装置，该装置包括电压传感器1、电流传感器2、旋转坐 标变换器4、锁相环5、光电编码器6、积分器7、减法器8、第二矩阵观测器10、电机输入转矩计 算器11、叶轮转矩计算器12、叶轮载荷谱生成器13和叶轮疲劳寿命估计器14;
[0092] 所述电压传感器1感应定子侧的三相电压，所述电流传感器2感应定子侧的三相电 流和转子侧的三相电流；
[0093 ]电流传感器2的输出端与旋转坐标变换器4的电流输入端相连；
[0094] 电压传感器1的定子三相电压输出端与锁相环5的输入端相连；
[0095] 锁相环5的定子电压旋转角度输出端分别与旋转坐标变换器4的定子电压旋转角 度输入端和减法器8的定子电压旋转角度输入端相连；
[0096] 光电编码器6米集电机转子机械转速，光电编码器6的转子机械转速输出端与积分 器7的输入端相连；
[0097] 积分器7的转子机械角度输出端与减法器8的转子机械角度输入端相连；
[0098] 减法器8的转子同步旋转角度输出端与旋转坐标变换器4的转子同步旋转角度输 入端相连；
[0099 ]旋转坐标变换器4输出端与第二转矩观测器10的电信号输入端相连；
[0100] 第二转矩观测器10的输出端与电机输入转矩计算器11的电磁转矩输入端相连；
[0101] 电机输入转矩计算器11的输出端与叶轮转矩计算器12的电机输入转矩输入端相 连；
[0102] 转速传感器15感应发电机转速，转速传感器15的两个输出端分别与电机输入转矩 计算器11的发电机转速输入端和叶轮转矩计算器12的发电机转速输入端相连；
[0103] 叶轮转矩计算器12的输出端与叶轮载荷谱生成器13的输入端相连；
[0104] 叶轮载荷谱生成器13的输出端与叶轮疲劳寿命估计器14的输入端相连。
[0105]图2是测量发电机电磁转矩的结构示意图。首先对定子侧三相电压usa、Usb、u sc进行 采样，并通过数字锁相环技术，得到定子电压旋转角度9S，同时，通过光电编码器，获得电机 转子机械转速，其积分为发电机转子机械角度，再与定子电压旋转角度9 S做差，即得到转 子同步旋转角度I。根据定子电压旋转角度9S、转子同步旋转角度I、定子侧的三相电流和 转子侧的三相电流，进行旋转坐标变换，分别得到定子在同步旋转坐标d-q轴上的电流及和 转子在同步旋转坐标d-q轴上的电流ids、iqs、idr和iqr，再根据公式(4)即可求得发电机电磁 转矩。abc/dq为旋转坐标变换器。可以同时采取【具体实施方式】一和【具体实施方式】二的两种 发电机电磁转矩的测量装置，以提高检测可靠性，如图2所示。
【具体实施方式】 [0106] 三:结合图1具体说明本实施方式，根据一所述的基于 转矩实时观测叶轮疲劳寿命的在线估计装置的在线估计方法，
[0107] 该方法包括以下步骤：
[0108] 步骤一、电压传感器1感应定子侧的三相电压和转子侧的三相电压，电流传感器2 感应定子侧的三相电流和转子侧的三相电流;转速传感器15感应发电机转速cog;
[0109] 步骤二、两相静止坐标变换器3将定子侧的三相电压转换为静止坐标系下的定子 侧电压Usa和定子侧电压Use，将定子侧的三相电流转换为静止坐标系下的定子侧电流isa和 定子侧电流ise，将转子侧的三相电压转换为静止坐标系下的转子侧电压Ura和转子侧电压 Urf!，将转子侧的三相电流转换为静止坐标系下的转子侧电流ira和转子侧电流irfi;
[0110] 步骤三、第一转矩观测器9结合发电机转速^，将定子侧电压Usa、定子侧电压Use、 定子侧电流isec、定子侧电流i sp、转子侧电压Urcc、转子侧电压Urf!、转子侧电流ira和转子侧电 流irf!转换为发电机电磁转矩1^;
[0111] 步骤四、电机输入转矩计算器11结合发电机转速CO g，将发电机电磁转矩转换为 发电机输入转矩Tg;
[0112] 步骤五、叶轮转矩计算器12结合发电机转速^^，将发电机输入转矩Tg转换为叶轮 气动转矩Ta;
[0113] 步骤六、叶轮载荷谱生成器13通过雨流量法根据叶轮气动转矩Ta得到叶轮载荷 谱；
[0114] 步骤七、将步骤六得到的叶轮载荷谱和预设叶轮参考载荷谱进行比较，得到比较 结果，该结果作为叶轮疲劳寿命评估值。
[0115] 图1是基于转矩实时观测叶轮疲劳寿命的在线估计方法的简化框图。首先，根据软 件模拟分析得到的叶轮实际转动惯量，通过公式(2)转换为高速轴电机侧为叶轮等效转动 惯量J^;根据发电机电磁转矩T em，通过公式(3)计算得到发电机的输入转矩Tg;最后，结合 转速传感器得到的发电机转速，通过式(1)得到叶轮气动转矩。
[0116] 通过该方法，避免直接使用应变式应力传感器测量方法带来的传感器安装及易受 环境影响问题，通过机组模型（即本发明所用的公式)和风电机组的电气测量参数直接对叶 轮转矩进行观测，提高了整个装置的可靠性。得到叶轮气动转矩后，通过雨流量法进行载荷 循环统计，即可得到载荷谱，同叶片的参考载荷谱曲线对照，即可在线估计叶轮的疲劳寿 命。
【具体实施方式】 [0117] 四：结合图1、图5至图7具体说明本实施方式，根据二 所述的基于转矩实时观测叶轮疲劳寿命的在线估计装置的在线估计方法，
[0118] 该方法包括以下步骤：
[0119] 步骤一、电压传感器1感应定子侧的三相电压，电流传感器2感应定子侧的三相电 流和转子侧的三相电流;转速传感器15感应发电机转速〇 g;光电编码器6米集电机转子机 械转速w r ;
[0120] 步骤二、锁相环5将定子侧的三相电压转换为定子电压旋转角度0S;
[0121] 步骤三、积分器7将电机转子机械转速cor转换为发电机转子机械角度；
[0122] 步骤四、减法器8结合定子电压旋转角度0S，将发电机转子机械角度转换为转子同 步旋转角度I;
[0123] 步骤五、旋转坐标变换器4结合定子电压旋转角度0S和转子同步旋转角度0r，分别 将定子侧的三相电流转换为旋转坐标系下的定子侧电流ids和定子侧电流i qs，将转子侧的 三相电流转换为旋转坐标系下的转子侧电流idr和转子侧电流iqr;
[0124] 步骤六、第二转矩观测器10将定子侧电流ids、定子侧电流iqs、转子侧电流i dr和转 子侧电流iqr转换为发电机电磁转矩Tan;
[0125] 步骤七、电机输入转矩计算器11结合发电机转速C〇g，将发电机电磁转矩Tem转换为 发电机输入转矩T g;
[0126] 步骤八、叶轮转矩计算器12结合发电机转速co g，将发电机输入转矩Tg转换为叶轮 气动转矩Ta;
[0127] 步骤九、叶轮载荷谱生成器13通过雨流量法根据叶轮气动转矩Ta得到叶轮载荷 谱；
[0128] 步骤十、将步骤九得到的叶轮载荷谱和预设叶轮参考载荷谱进行比较，得到比较 结果，该结果作为叶轮疲劳寿命评估值。
[0129] 可以同时采取【具体实施方式】三和【具体实施方式】四的两种发电机电磁转矩的获得 方法，以提高检测可靠性。
[0130]图1是基于转矩实时观测叶轮疲劳寿命的在线估计方法的简化框图。首先，根据软 件模拟分析得到的叶轮实际转动惯量，通过公式2转换为高速轴电机侧为叶轮等效转动惯 量根据发电机电磁转矩!-，通过公式3计算得到发电机的输入转矩T g;最后，结合转速 传感器得到的发电机转速，通过式1得到叶轮气动转矩。
[0131] 通过该方法，避免直接使用应变式应力传感器测量方法带来的传感器安装及易受 环境影响问题，通过机组模型（即本发明所用的公式)和风电机组的电气测量参数直接对叶 轮转矩进行观测，提高了整个装置的可靠性。得到叶轮气动转矩后，通过雨流量法进行载荷 循环统计，即可得到载荷谱，同叶片的参考载荷谱曲线对照(参考载荷谱曲线由厂家提供或 者仿真模拟计算)即可在线估计叶轮的疲劳寿命。
[0132] 图6为采用【具体实施方式】三和【具体实施方式】四所述的基于转矩实时观测叶轮疲劳 寿命的在线估计方法得到的发电机电磁转矩计算值与实际值对比图，其中A为采用具体实 施方式三的方法得到的曲线，B为采用【具体实施方式】四的方法得到的曲线，C为发电机电磁 转矩实际值。图7是采用本实施方式得到的发电机输入转矩值计算值与实际值对比图，图中 虚线的曲线为计算值，实线的曲线为实际值，可以看到采用本实施方式得到的计算值与实 际值仅存在微小差异。
[0133] 图5是雨流量法进行载荷循环统计在示意图，在得到叶轮气动转矩!^后通过雨流 量法对载荷循环进行统计，雨滴法如图5所示，横轴为载荷，即计算得到的叶轮气动转矩T a， 纵轴为时间。雨滴从波谷A处滴落，当到达B时，沿B落下，直到降落到与E点平行，因 E点同为 波谷且幅值比A处更小。同样，B点雨滴流到C，沿C下落到C'，并下落到E最后落下。C处雨滴落 到D并落下，同样止于与E平行处。D处雨滴下落到E，在此过程中遇到雨滴C'，停止，C-D-C'构 成一个载荷循环。类似的，F-G-F'构成一个载荷循环。记录载荷循环的幅值，并统计每个幅 值的循环次数，即可得到载荷谱。
【主权项】
1. 基于转矩实时观测叶轮疲劳寿命的在线估计装置，其特征在于，该装置包括电压传 感器（1 )、电流传感器(2)、两相静止坐标变换器(3)、第一转矩观测器(9)、电机输入转矩计 算器(11)、叶轮转矩计算器(12)、叶轮载荷谱生成器（13)、叶轮疲劳寿命估计器(14)和转速 传感器(15); 所述电压传感器（1)感应定子侧的三相电压和转子侧的三相电压，所述电流传感器(2) 感应定子侧的三相电流和转子侧的三相电流； 电压传感器（1)的输出端与两相静止坐标变换器（3)的电压输入端相连，电流传感器 (2)的输出端与两相静止坐标变换器(3)的电流输入端相连； 两相静止坐标变换器(3)的输出端与第一转矩观测器(9)的电信号输入端相连； 第一转矩观测器(9)的输出端与电机输入转矩计算器(11)的电磁转矩输入端相连； 电机输入转矩计算器（11)的输出端与叶轮转矩计算器（12)的电机输入转矩输入端相 连； 转速传感器（15)感应发电机转速，转速传感器（15)的三个输出端分别与第一转矩观测 器(9)的发电机转速输入端、电机输入转矩计算器（11)的发电机转速输入端和叶轮转矩计 算器(12)的发电机转速输入端相连； 叶轮转矩计算器(12)的输出端与叶轮载荷谱生成器(13)的输入端相连； 叶轮载荷谱生成器(13)的输出端与叶轮疲劳寿命估计器(14)的输入端相连。2. 基于转矩实时观测叶轮疲劳寿命的在线估计装置，其特征在于，该装置包括电压传 感器（1)、电流传感器(2)、旋转坐标变换器(4)、锁相环(5)、光电编码器(6)、积分器(7)、减 法器(8)、第二矩阵观测器(10)、电机输入转矩计算器（11)、叶轮转矩计算器（12)、叶轮载荷 谱生成器(13)和叶轮疲劳寿命估计器(14); 所述电压传感器（1)感应定子侧的三相电压，所述电流传感器(2)感应定子侧的三相电 流和转子侧的三相电流； 电流传感器(2)的输出端与旋转坐标变换器(4)的电流输入端相连； 电压传感器(1)的定子三相电压输出端与锁相环(5)的输入端相连； 锁相环(5)的定子电压旋转角度输出端分别与旋转坐标变换器(4)的定子电压旋转角 度输入端和减法器(8)的定子电压旋转角度输入端相连； 光电编码器(6)米集电机转子机械转速，光电编码器(6)的转子机械转速输出端与积分 器(7)的输入端相连； 积分器(7)的转子机械角度输出端与减法器(8)的转子机械角度输入端相连； 减法器(8)的转子同步旋转角度输出端与旋转坐标变换器(4)的转子同步旋转角度输 入端相连； 旋转坐标变换器(4)输出端与第二转矩观测器(10)的电信号输入端相连； 第二转矩观测器(10)的输出端与电机输入转矩计算器(11)的电磁转矩输入端相连； 电机输入转矩计算器（11)的输出端与叶轮转矩计算器（12)的电机输入转矩输入端相 连； 转速传感器（15)感应发电机转速，转速传感器（15)的两个输出端分别与电机输入转矩 计算器(11)的发电机转速输入端和叶轮转矩计算器(12)的发电机转速输入端相连； 叶轮转矩计算器(12)的输出端与叶轮载荷谱生成器(13)的输入端相连； 叶轮载荷谱生成器(13)的输出端与叶轮疲劳寿命估计器(14)的输入端相连。3. 根据权利要求1所述的基于转矩实时观测叶轮疲劳寿命的在线估计装置的在线估计 方法，其特征在于，该方法包括以下步骤： 步骤一、电压传感器（1)感应定子侧的三相电压和转子侧的三相电压，电流传感器（2) 感应定子侧的三相电流和转子侧的三相电流;转速传感器(15)感应发电机转速c〇g; 步骤二、两相静止坐标变换器(3)将定子侧的三相电压转换为静止坐标系下的定子侧 电压Uscc和定子侧电压Usfi，将定子侧的三相电流转换为静止坐标系下的定子侧电流iscc和定 子侧电流id，将转子侧的三相电压转换为静止坐标系下的转子侧电压Ura和转子侧电压Urf!， 将转子侧的三相电流转换为静止坐标系下的转子侧电流和转子侧电流irfi; 步骤三、第一转矩观测器(9)结合发电机转速COg，将定子侧电压Usa、定子侧电压Use、定 子侧电流i SCi、定子侧电流i sf!、转子侧电压Ura、转子侧电压Urfi、转子侧电流ira和转子侧电流 irf!转换为发电机电磁转矩Te*; 步骤四、电机输入转矩计算器（11)结合发电机转速《8，将发电机电磁转矩Te3m转换为发 电机输入转矩Tg; 步骤五、叶轮转矩计算器（12)结合发电机转速ω g，将发电机输入转矩Tg转换为叶轮气 动转矩Ta; 步骤六、叶轮载荷谱生成器(13)通过雨流量法根据叶轮气动转矩Ta得到叶轮载荷谱； 步骤七、将步骤六得到的叶轮载荷谱和预设叶轮参考载荷谱进行比较，得到比较结果， 该结果作为叶轮疲劳寿命评估值。4. 根据权利要求2所述的基于转矩实时观测叶轮疲劳寿命的在线估计装置的在线估计 方法，其特征在于，该方法包括以下步骤： 步骤一、电压传感器（1)感应定子侧的三相电压，电流传感器(2)感应定子侧的三相电 流和转子侧的三相电流;转速传感器（15)感应发电机转速ω g;光电编码器(6)采集电机转 子机械转速ω r ; 步骤二、锁相环(5)将定子侧的三相电压转换为定子电压旋转角度0S; 步骤三、积分器(7)将电机转子机械转速换为发电机转子机械角度； 步骤四、减法器(8)结合定子电压旋转角度03，将发电机转子机械角度转换为转子同步 旋转角度I; 步骤五、旋转坐标变换器(4)结合定子电压旋转角度03和转子同步旋转角度0r，分别将 定子侧的三相电流转换为旋转坐标系下的定子侧电流ids和定子侧电流iqs，将转子侧的三 相电流转换为旋转坐标系下的转子侧电流idr和转子侧电流i qr; 步骤六、第二转矩观测器（10)将定子侧电流ids、定子侧电流iqs、转子侧电流idr和转子 侧电流iqr转换为发电机电磁转矩Te*; 步骤七、电机输入转矩计算器（11)结合发电机转速COg，将发电机电磁转矩Tem转换为发 电机输入转矩Tg; 步骤八、叶轮转矩计算器（12)结合发电机转速ω g，将发电机输入转矩Tg转换为叶轮气 动转矩Ta; 步骤九、叶轮载荷谱生成器(13)通过雨流量法根据叶轮气动转矩Ta得到叶轮载荷谱； 步骤十、将步骤九得到的叶轮载荷谱和预设叶轮参考载荷谱进行比较，得到比较结果， 该结果作为叶轮疲劳寿命评估值。
【文档编号】H02P21/14GK105958886SQ201610320171
【公开日】2016年9月21日
【申请日】2016年5月13日
【发明人】杨荣峰, 苏勋文, 张学广, 徐殿国, 张磊, 孙运涛
【申请人】国家电网公司, 国网山东省电力公司电力科学研究院, 许继集团有限公司, 哈尔滨工业大学, 黑龙江科技大学, 许昌许继风电科技有限公司