一种以风力发电机组叶片运行轨迹计算净空距离方法与流程

1.本发明涉及风力发电测控技术领域，尤其是一种以风力发电机组叶片运行轨迹计算净空距离方法。


背景技术：

2.随着风力发电的大力发展，大家在追求快速批量发展过程中，风力发电机组的安全保障变得越来越重要。同时，随着2010年前大力发展的风电场，目前随着时间、机组性能、寿命都在缩减。风电机组扫塔、倒塔、叶片断裂等事故也是频繁出现。另外，目前三北平原地区风力发电机组装机容量已趋于饱和状态，现在风力发电机组安装正逐步向海上和山地区域发展，从而推动了高塔筒、长叶片技术发展。通过以gps监测的叶片轨迹，计算出任何时刻叶片和塔筒的净空距离。并集合风力发电机组的状态数据，完成数据关联分析和数学建模，在风力发电机组运行过程中，可提前预知风力发电机组的叶片净空距离。将叶片净空距离反馈至风力发电机组控制系统，通过闭环控制，一旦叶片净空距离小于安全值，通过控制并干预当前状态，避免叶片扫塔，保证叶片和塔筒的安全，从而避免事故发生。
3.目前，风电行业叶片净空距离均采用监测系统直接测量获得，主要分为以下三种：1、机舱安装激光雷达或则毫米波雷达方案：通过在机舱底部安装激光激光或则毫米波发生和接受装置，通过特定光束，一般为3束光。通过光反射原理，一旦反射距离发生改变，认为当前叶片到塔筒距离为当前光束到塔筒距离。该方案缺点：只有三个点位测量，数据误差大，无法实时感知当前叶片净空距离。同时，方案受气象天气影响较大，在多雾天气，数据基本不可用。
4.2、机舱安装摄像头方案：通过在机舱安装摄像头，拍摄叶片运行视频，通过机器学习的方法来推算当前叶片到塔筒的净空距离。该方案缺点：需要大量经验图片，容易漏选。摄像机摄像范围有限，只能拍摄固定区域的叶片情况。同时整套方案受外界天气影响较大，无法有效识别当前图像。另外，由于机器学习算法本身的缺陷，存在一定的误差和误报，无法精确预知当前叶片和塔筒净空距离。
5.3、塔筒相应高度安装雷达测距装置：通过在塔筒上安装导轨，导轨上测距装置根据机组当前偏航状态进行跟踪，并实时测量叶片净空距离。该方法安装复杂，同时安装过程可能会损坏塔筒，可能会带来新的安全隐患。同时该方案成本价格昂贵。


技术实现要素：

6.本发明目的在于：针对上述问题，提出了一种以风力发电机叶片运行轨迹计算叶片净空方法，克服了现有技术无法精准获取叶片净空距离，也无法动态高效的的对机组安全性能进行防护的问题。
7.本发明是通过下述方案来实现的：一种以风力发电机组叶片运行轨迹计算净空距离方法，包括以下步骤：通过叶片轨迹计算出当前叶尖三维坐标；
建立以地球圆心为原点坐标系；叶尖三维坐标换算，同时建立塔筒结构模型；最终获取以塔筒圆心为原点的坐标系统的叶尖三维坐标；计算出过轮毂到塔筒顶部圆最短点投影到塔筒底内点并并行过轮毂到塔筒顶部圆最短点动态直线方程；计算叶片净空距离；通过机器学习方法构建叶片净空距离的预测模型；据机组自身及环境状态、叶片净空距离，实时判断叶片塔筒安全。
8.在构建叶片净空距离的预测模型之后设置有叶片净空距离预测步骤。
9.叶片轨迹计算出当前叶尖三维坐标，借助已有风力发电机组叶片轨迹测量系统给出的叶尖定位位置数据；由于系统给出原始数据为经纬度和海拔，使用换算算法计算出当前叶尖的三维坐标。
10.构建叶尖三维坐标具体步骤为，通过测量塔筒底部经纬度海拔，计算出塔筒底部圆心的坐标为；同时计算得出的叶尖的得坐标为x,y,z；通过两个坐标计算，叶根轨迹的三维坐标分别为（）。
11.构建最短点动态直线方程步骤为：设塔筒上面连接机舱的部分的半径为r，塔基半径为r，为偏航角度，上塔筒顶部到塔筒底部的垂直距离为h，则塔筒上方距叶片最近点坐标为；则塔筒塔基距叶片最近的点为；现设叶片和塔筒发生碰撞的区域为从塔基往上的h米高度，那么以高度为h的与叶片最近的点在塔筒上的坐标点为，那么该点到塔基的投影为，与原点的向量为，该截面在塔基的投影线为:因此叶片到该截面的距离可以抽象为该点在塔基的投影点到该截面到塔基的投影线的距离。因此,则叶片到塔筒的最短距离d的长度；。
12.计算叶片净空距离具体为，根据叶尖当前位置坐标投影到塔底平面点的坐标，通过点到直接垂线计算方式，计算出当前的最短距离，从而计算出当前叶片净空距离。
13.数据建模具体步骤为，根据已计算出的叶片净空距离及主控私有通讯协议获取的风力发电机组状态数据，通过机器学习方法对上述数据进行分析建模和训练优化，得到基于风力发电机组状态数据的叶片净空距离的预测模型。
14.据机组状态和叶片净空距离，实时判断叶片塔筒安全的具体步骤为：通过计算给出的当前叶片净空距离和预知的叶片净空距离；通过通讯，数据接入主控系统，主控系统根据当前和预知叶片净空距离及风机当前运行状态，并结合规则抽象法判断当前机组状态当前安全状态；一旦叶片净空距离处于危险状态及等级，机组通过降低转速、变桨等降载控制策略避免叶片扫塔，从而保证叶片和塔筒安全，提高机组安全性能；最后，通过当前和预知叶片净空距离及风机当前运行状态结合规则抽象法发现当前机组安全状态，若机组回归安全状态，则脱离降载控制，正常发电，保障机组的发电量。
15.所述降载控制策略，通过当前和预知叶片净空距离及风机当前运行状态，机组危险状态及等级等数据结合风力发电机组主控控制逻辑和风力发电机组运行专家经验法得到的先进的降载控制策略，在保证机组安全运行的情况下，极大的保留机组的发电能力。
16.所述规则抽象法包括风力发电机组主控控制逻辑、风力发电机组运行状态和风力发电机组运行专家经验法一种或三种任意组合。
17.综上所述，由于采用了上述技术方案，本发明的有益效果是：1、本发明可以通过采取并结合规则抽象法判断当前机组状态当前安全状态；一旦叶片净空距离处于危险状态及等级，机组通过降低转速、变桨等降载控制策略避免叶片扫塔，从而保证叶片和塔筒安全，提高机组安全性能；通过当前和预知叶片净空距离及风机当前运行状态结合规则抽象法发现当前机组安全状态，若机组回归安全状态，则脱离降载控制，正常发电，保障机组的发电量，整体上提升安全防护性能。
18.2、本发明通过实时感知当前净空距离，净空距离数据进入风力发电机组控制系统。并通过风力发电机组当前状态和算法结果预知下一时刻是否属于危险状态，一旦到达危险极限，风力发电机组控制系统向风机发出控制指令，或停机或偏航，让风力发电机组叶片避免扫塔，保证风力发电机组叶片塔筒处于安全状态。另外，由于高精度叶片净空距离的计算，通过数据分析可算出风力发电机组叶片扫塔的极限状态，并集合叶片本身载荷和相关技术指标，可为风力发电机组叶片降载提供数据支撑。
19.3、本发明采用使获得的叶片净空距离更准确。
附图说明
20.图1是以风力发电机组叶片运行轨迹计算叶片净空距离方法步骤图；图2是以风力发电机组叶片运行轨迹计算叶片净空距离方法塔筒坐标图。
具体实施方式
21.本说明书中公开的所有特征，或公开的所有方法或过程中的步骤，除了互相排斥的特征和/或步骤以外，均可以以任何方式组合。
22.本说明书（包括任何附加权利要求、摘要）中公开的任一特征，除非特别叙述，均可被其他等效或具有类似目的的替代特征加以替换。即，除非特别叙述，每个特征只是一系列
等效或类似特征中的一个例子而已。
23.在本发明的描述中，需要理解的是，术语“上”、“下”、“左”、“右”等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本发明和简化描述，而不是指示或暗示所指的设备或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本发明的限制。
24.此外，术语“第一”、“第二”等仅用于描述目的，而不能理解为指示或暗示相对重要性或隐含指明所指示的技术特征的数量。由此，限定有“第一”、“第二”等的特征可以明示或隐含地包括一个或多个该特征。
25.实施例1如图1~图2所示，一种以风力发电机组叶片运行轨迹计算净空距离方法，其包括以下步骤：通过已有的风力发电机组叶片轨迹监测系统获取叶尖的gps定位信息数据。
26.通过当前gps定位信息，使用经纬度、海拔转换算法，计算出当前叶尖以地球球心为原点的三维坐标。
27.查询塔筒设计图纸，建立塔筒的结构模型，为一个正圆锥形。
28.使用rtk技术测量出塔筒底部圆心的定位数据，并转换成平面坐标系中的三维坐标，并建立以塔筒底部圆心为原点的坐标系。
29.完成叶尖已有的三维坐标转换，最终获取以塔筒圆心为原点的坐标系统的叶尖三维坐标。
30.根据点到直线原理，并集合当前机组偏航角度和塔筒参数，计算出过轮毂到塔筒顶部圆最短点投影到塔筒底内点并并行过轮毂到塔筒顶部圆最短点动态直线方程。
31.点到直线最短距离为垂线，计算当前叶尖点到当前直线的距离。
32.集合已经计算出的叶片净空距离和机组状态数据，完成数据建模，建立数学模型。
33.根据当前机组状态和数据模型，可预测叶片净空距离。
34.叶片净空距离实时闭环到主控系统中，根据机组状态和叶片净空距离，实时判断，一旦超过安全阈值，及时干预机组控制，保证叶片和塔筒安全。
35.其中通过叶片轨迹计算出当前叶尖三维坐标，需要借助已有风力发电机组叶片轨迹测量系统给出的叶尖定位位置数据；由于系统给出原始数据为经纬度和海拔，需要使用换算算法计算出当前叶尖的三维坐标。
36.本方法中计算出当前叶尖的三维坐标的具体步骤如下；该算法通过地理坐标的定义规则，计算公式如下：第一步：查询相关参数ra_wgs84=6378137.0//地球长轴半径，单位：米f_wgs84=1/298.257223563//扁率rb_wgs84=(ra_wgs84*(1-f_wgs84))//短轴半径，单位:米e_wgs84=(0.081819190842621494334802451753867)//第一偏心率e2_wgs84=(0.0066943799901413169961372335400465)//第一偏心率的平方第二步：换算当前系统测量的经纬度为：latitude（经度）、longitude（纬度）、height（海拔高
度）cosphi = cos(rad(latitude))coslamda = cos(rad(longitude))sinphi = sin(rad(latitude))sinlamda = sin(rad(longitude))n = ra_wgs84 / sqrt(1-e2_wgs84 * sinphi * sinphi)，其中cosphi、coslamda、sinphi、sinlamda和n 均为自定义的中间常量；则当前叶尖的三维坐标可表示为：x = (n height) * cosphi * coslamday = (n height) * cosphi * sinlamda z = (n * (1.0
ꢀ‑ꢀ
e2_wgs84) height) * sinphi其中建立塔筒结构模型，需要利用塔筒机械结构，塔筒属于锥形结构，并查询塔筒相关参数，包括塔筒顶部小圆半径，塔筒底部大圆半径，塔筒高度。绘制出塔筒结构模型图。
37.步骤如下：1、塔筒结构属于锥形结构，查询塔筒相关参数，包括：顶部小圆半径:r，底部大圆半径;r ，塔筒中心轴高度为：h；建立塔筒锥形结构模型如下如图1所示；其中在建立以塔筒底部圆心为原点坐标系，需要使用rtk技术，测量塔筒底部圆心经纬度和海拔，并通过坐标系换算，计算出当前塔筒圆心三维坐标，并以该点为原点绘制出塔筒的平面坐标系；其中叶尖三维坐标换算，由于直接测量后换算的叶尖三维坐标，是以地球球心为原点球坐标下的三维坐标；使用相关算法，完成坐标系转换成，最终转换成以塔筒底部圆心为原点的塔筒坐标系中的三维坐标；其中具体步骤为：通过测量塔筒底部经纬度海拔，并通过前面公式，计算出塔筒底部圆心的坐标为；通过前面公式，计算得出的叶尖的得坐标为x,y,z。通过两个坐标计算，叶根轨迹的三维坐标分别为（）。
38.其中建立最短距离直线动态方程，是根据点到直线原理，根据机组偏航角度值和塔筒相关参数，计算出过轮毂到塔筒顶部圆最短点投影到塔筒底内点并并行过轮毂到塔筒顶部圆最短点动态直线方程；其具体步骤为如下：如图2所示，塔筒一般为上面半径较小，下面半径较大的圆台，那么设塔筒上面连接机舱的部分的半径为r，塔基半径为r，为偏航角度（从主控状态数据表中查询），上塔筒顶部到塔筒底部的垂直距离（塔筒长度）为h，则塔筒上方距叶片最近点（轮毂中心到底部塔筒圆周上最近点，图中a点）坐标为；则塔筒塔基距叶片最近的点（图中c点）为，现假设叶片和塔筒发生碰撞的区域为从塔基往上的h米高度，那么设以高度为h的与叶片最近的点为塔筒上的坐标点（图中h0点），则图中h0点的坐标可由h0、h1和c点及a0、a1和c点购成的相似三角形得到，为
，那么该点到塔基的投影（h1点）坐标点位为，与原点的向量为，那么，可以h1与原点的向量为法向量，h0和h1点位为截面坐标点，形成截面，那么该截面在塔基的投影线（图中l）为:因此叶片到该截面的距离可以抽象为该点在塔基的投影点到该截面到塔基的投影线的距离。因此,则叶片到塔筒的最短距离d的长度；上述公式即为最短距离直线动态方程。
39.其中计算叶片净空距离，根据点直线最短距离原理，根据叶尖当前位置坐标投影到塔底平面点的坐标，通过点到直接垂线计算方式，计算出当前的最短距离，从而计算出当前叶片净空距离；其中数据建模具体步骤为，根据已计算出的叶片净空距离及主控私有通讯协议获取的风力发电机组状态数据，通过专家经验法、聚类方法、深度学习和多重线性分析方法等机器学习方法对上述数据进行分析建模和训练优化，得到基于风力发电机组状态数据的叶片净空距离的预测模型；其中叶片净空距离的预测模型，可依据历史运行的数据和风力发电机组实时状态数据，提前预知未来的一段时间的叶片净空距离；其中闭环控制保证叶片塔筒安全步骤为：通过计算给出的当前叶片净空距离和预知的叶片净空距离；通过通讯，数据接入主控系统，主控系统根据当前和预知叶片净空距离及风机当前运行状态，并结合规则抽象法判断当前机组状态当前安全状态；一旦叶片净空距离处于危险状态及等级，机组通过降低转速、变桨等降载控制策略避免叶片扫塔，从而保证叶片和塔筒安全，提高机组安全性能；最后，通过当前和预知叶片净空距离及风机当前运行状态结合规则抽象法发现当前机组安全状态，若机组回归安全状态，则脱离降载控制，正常发电，保障机组的发电量。
40.所述降载控制策略，则通过上述数据(包括当前和预知叶片净空距离及风机当前运行状态，机组危险状态及等级等数据)结合风力发电机组主控控制逻辑和风力发电机组运行专家经验法得到的先进的降载控制策略，在保证机组安全运行的情况下，极大的保留机组的发电能力。
41.所述规则抽象法包括风力发电机组主控控制逻辑、风力发电机组运行状态和风力发电机组运行专家经验法一种或三种任意组合。
42.以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。