CN108843492B - 通过无人机进行风机偏航角测算方法及系统 - Google Patents

通过无人机进行风机偏航角测算方法及系统 Download PDF

Info

Publication number
CN108843492B
CN108843492B CN201810627916.8A CN201810627916A CN108843492B CN 108843492 B CN108843492 B CN 108843492B CN 201810627916 A CN201810627916 A CN 201810627916A CN 108843492 B CN108843492 B CN 108843492B
Authority
CN
China
Prior art keywords
point
unmanned aerial
aerial vehicle
blades
calculating
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
Application number
CN201810627916.8A
Other languages
English (en)
Other versions
CN108843492A (zh
Inventor
刘迅
陈小明
尚黎民
叶华
柯严
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shanghai Clobotics Technology Co ltd
Original Assignee
Shanghai Clobotics Technology Co ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Shanghai Clobotics Technology Co ltd filed Critical Shanghai Clobotics Technology Co ltd
Priority to CN201810627916.8A priority Critical patent/CN108843492B/zh
Publication of CN108843492A publication Critical patent/CN108843492A/zh
Application granted granted Critical
Publication of CN108843492B publication Critical patent/CN108843492B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F03—MACHINES OR ENGINES FOR LIQUIDS; WIND, SPRING, OR WEIGHT MOTORS; PRODUCING MECHANICAL POWER OR A REACTIVE PROPULSIVE THRUST, NOT OTHERWISE PROVIDED FOR
    • F03D—WIND MOTORS
    • F03D7/00—Controlling wind motors 
    • F03D7/02—Controlling wind motors  the wind motors having rotation axis substantially parallel to the air flow entering the rotor
    • F03D7/0204—Controlling wind motors  the wind motors having rotation axis substantially parallel to the air flow entering the rotor for orientation in relation to wind direction
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C11/00—Photogrammetry or videogrammetry, e.g. stereogrammetry; Photographic surveying
    • G01C11/02—Picture taking arrangements specially adapted for photogrammetry or photographic surveying, e.g. controlling overlapping of pictures
    • G—PHYSICS
    • G05—CONTROLLING; REGULATING
    • G05D—SYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
    • G05D1/00—Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
    • G05D1/10—Simultaneous control of position or course in three dimensions
    • G05D1/101—Simultaneous control of position or course in three dimensions specially adapted for aircraft
    • 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
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00—Energy generation through renewable energy sources
    • Y02E10/70—Wind energy
    • Y02E10/72—Wind turbines with rotation axis in wind direction

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • General Physics & Mathematics (AREA)
  • Chemical & Material Sciences (AREA)
  • Mechanical Engineering (AREA)
  • General Engineering & Computer Science (AREA)
  • Multimedia (AREA)
  • Combustion & Propulsion (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Sustainable Energy (AREA)
  • Sustainable Development (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Automation & Control Theory (AREA)
  • Navigation (AREA)

Abstract

本发明提供了一种通过无人机进行风机偏航角测算方法及系统，包括如下步骤：控制无人机以风塔高度绕风机飞行，通过图像传感器采集叶轮的视频流；对视频流中的叶片进行检测，对三个叶片进行实时跟踪，并实时计算三个叶片的相对位置及重叠度；当检测到两个叶片完全重叠时，认定此时无人机飞行到风轮平面β上，读取此时位置传感器获取的点P1的位置信息；根据点P1的位置信息计算与点P1以风塔呈轴对称分别的点P2的第一位置信息；根据点P1的位置信息、点P2的第一位置信息以及地球质心计算出风轮平面β，进而计算出风轮平面的偏航角。本发明能够通过无人机测算出风轮平面的偏航角，从而能够便于风机进行建模分析，为实现风机的综合检测提供了便利。

Description

通过无人机进行风机偏航角测算方法及系统
技术领域
本发明涉及风机检测，具体地，涉及一种通过无人机进行风机偏航角测算方法及系统。
背景技术
叶片作为风电系统关键部件，数量多且型面结构复杂，检测难度大。而传统的检测方法是采用专用工装对叶片进行检测，效率低，影响了叶片的检测精度和效率。近年来，随着三坐标测量的普及，建立三坐标测量对叶片的跟踪检测就成为一种趋势，偏航角和方位角是三坐标测量中两个重要参数。
当风连续不断改变方向，风电机叶片轮无法及时跟踪风向，使得风轮转轴与风向不平行，叶轮处于偏航状态，同时作用在叶片上的空气动力载荷也使得桨翼处于运动状态。因此，精确测量偏航角和转角对于叶片以及风电系统性能的研究和检测至关重要。
发明内容
针对现有技术中的缺陷，本发明的目的是提供一种通过无人机进行风机偏航角测算方法及系统。
根据本发明提供的通过无人机进行风机偏航角测算方法，所述风机包括风塔和设置在风塔顶端的叶轮，所述叶轮包括三个沿周向均匀分布的叶片，包括如下步骤：
步骤S1：控制无人机以风塔高度绕风机飞行，当无人机在飞行过程中，通过图像传感器采集所述叶轮的视频流；
步骤S2：对所述视频流中的叶片进行检测，当检测到风机的三个叶片时，对三个叶片进行实时跟踪，并实时计算三个叶片的相对位置及重叠度；
步骤S3：当检测到两个叶片完全重叠时，认定此时无人机飞行到风轮平面β上，读取此时位置传感器获取的点P1的位置信息；
步骤S4：根据点P1的位置信息计算与点P1以风塔呈轴对称分布的点P2的第一位置信息；
步骤S5：根据点P1的位置信息、点P2的第一位置信息以及地球质心计算出风轮平面β，进而计算出所述风轮平面的偏航角。
优选地，在步骤S3和步骤S4之间还包括如下步骤：
-让无人机继续飞行，当再次检测到两个叶片完全重叠时，读取此时位置传感器获取的点P2的第二位置信息，通过点P2的第二位置信息对点P2的第一位置信息进行验证。
优选地，所述检测到两个叶片完全重叠时的点P1通过如下方式计算得出：
P1＝P[min(τ)] (1-1)
Figure GDA0002311521400000021
其中，τ为二值图像流ti中目标行数累加值，P为无人机的实时位置，目标行数累加值τ为根据ti(x,y)的值累加生成，根据式(1-2)，当ti(x,y)＝1时累加一次；
当目标行数累加值τ为最小时，确定所述两个叶片完全重叠；
其中x表示二值图像流ti的x轴坐标值；y表示二值图像流ti的y轴坐标值。
优选地，还包括如下点P1、P2的位置验证步骤：
步骤M1：将点P1的位置信息、P2的第一位置信息转换至地球坐标系(Xe,Ye,Ze)，点P1、P2的位置信息通过GPS模块，采用经度、纬度、高度表示，转换计算公式为：
Figure GDA0002311521400000022
N为纬度B处的卯酉圈曲率半径，E为地球第一偏心率，
Figure GDA0002311521400000023
E＝(a2-b2)/a2,a为地球长半径，b为地球短半径，B为位置信息中的纬度，L为位置信息中的风塔经度，H为位置信息中的风塔高度；
步骤M2：验证点P2、P1在地球坐标的下的位置关系，即
Figure GDA0002311521400000024
Figure GDA0002311521400000025
其中
Figure GDA0002311521400000026
为点P2与点P1之间的直线距离，
Figure GDA0002311521400000027
为点P1离风轮中心距离，
Figure GDA0002311521400000028
为点P2离风轮中心距离；
步骤M3：计算精确度ratio,判断精确度ratio是否满足98％<ratio<102％；
Figure GDA0002311521400000029
优选地，还包括如下转角γ的测算步骤：
步骤N1：让无人机位于风机正前方，距离风塔底端设定距离的点P0位置以垂直升空至风塔高度，获取点PT的位置信息；
步骤N2：读取图像传感器于点PT采用的叶轮图像，将所述叶轮图像去除噪声；
步骤N3：对叶轮图像进行边缘检测，检测出目标叶片信息，通过角点检测计算叶尖点坐标，计算三个叶片几何中心点即为风轮中心在图像中的坐标Pwind centre；
步骤N4：连接叶尖点坐标与Pwind centre坐标确定一目标直线，进一步的计算该目标直线的斜率,即可得到该叶片的转角大小；
其中点PT表示无人机升空至所述风塔高度时的无人机的位置点。
优选地，所述无人机上搭载位置传感器、图像传感器和机载电脑；
所述位置传感器、所述图像传感器连接所述机载电脑；所述位置传感器用于实时读取无人机位置信息，图像传感器用于拍摄风机叶片生成风机叶片图像，机载电脑用于无人机位置信息和风机叶片图像的处理。
优选地，所述步骤S5具体为：
-将点P1的位置信息、点P2的第一位置信息转换至地球坐标系，再根据转换至地球坐标系中点P1、点P2及地球质心计算出风轮平面β，进而可得到在地球坐标系中，风轮平面β在地球坐标系(Xe,Ye,Ze)中的方向向量，进一步计算所述方向向量与东北天坐标系中Y轴之间形成的偏航角。
本发明提供的通过无人机进行风机偏航角测算系统，用于实现所述的通过无人机进行风机偏航角测算方法，包括：
飞行控制模块，用于控制无人机以风塔高度绕风机飞行，当无人机在飞行过程中，通过图像传感器采集叶片的视频流；
重叠度计算模块，用于对所述视频流中的叶片进行检测，当检测到风机的三个叶片时，对三个叶片进行实时跟踪，并实时计算三个叶片的相对位置及重叠度；
位置信息获取模块，用于当检测到两个叶片完全重叠时，认定此时无人机飞行到风轮平面β上，读取此时位置传感器获取的点P1的位置信息；
位置信息计算模块，用于根据点P1的位置信息计算与点P1以风塔呈轴对称分布的点P2的第一位置信息；
偏航角计算模块，用于根据点P1的位置信息、点P2的第一位置信息以及地球质心计算出风轮平面β，进而计算出所述风轮平面的偏航角。
与现有技术相比，本发明具有如下的有益效果：
本发明能够通过无人机测算风轮平面的方向向量与东北天坐标系Y轴之间形成的偏航角，从而能够根据风向调整风轮转轴，便于提高了风机的发电效率；本发明能够通过无人机测算出风轮平面的偏航角，从而能够便于风机进行建模分析，为实现风机的综合检测提供了便利；本发明中在无人机一次飞行过程中，不仅能够实现偏航角的测量，而且能够实现叶片转角的大小测量，提高了测量效果。
附图说明
通过阅读参照以下附图对非限制性实施例所作的详细描述，本发明的其它特征、目的和优点将会变得更明显：
图1为本发明中通过无人机进行风机偏航角测算方法的步骤流程图；
图2为本发明中通过无人机进行风机偏航角测算方法的原理示意图；
图3为本发明中风轮平面β的示意图；
图4为本发明中偏航角示意图；
图5为本发明中通过视觉检测方法计算转角γ的示意图；
图6为本发明中通过无人机进行风机偏航角测算系统的模块示意图。
图中：
1为第一平面δ；
2为飞行轨迹曲线s；
3为风轮平面β；
4为直线l；
5为点P1；
6为点P2。
具体实施方式
下面结合具体实施例对本发明进行详细说明。以下实施例将有助于本领域的技术人员进一步理解本发明，但不以任何形式限制本发明。应当指出的是，对本领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干变形和改进。这些都属于本发明的保护范围。
针对传统叶片检测的难度与技术上不足和确定偏航角和转角对于风电行业研究的重要性，本发明提出一种通过无人机进行风机偏航角测算方法及系统。
图1为本发明中通过无人机进行风机偏航角测算方法的步骤流程图，如图1所示，本发明提供的通过无人机进行风机偏航角测算方法，所述风机包括风塔和设置在风塔顶端的叶轮，所述叶轮包括三个沿周向均匀分布的叶片，包括如下步骤：
步骤S1：控制无人机以风塔高度绕风机飞行，当无人机在飞行过程中，通过图像传感器采集所述叶轮的视频流；
步骤S2：对所述视频流中的叶片进行检测，当检测到风机的三个叶片时，对三个叶片进行实时跟踪，并实时计算三个叶片的相对位置及重叠度；
步骤S3：当检测到两个叶片完全重叠时，认定此时无人机飞行到风轮平面β上，读取此时位置传感器获取的点P1的位置信息；
步骤S4：根据点P1的位置信息计算与点P1以风塔呈轴对称分布的点P2的第一位置信息；
步骤S5：根据点P1的位置信息、点P2的第一位置信息以及地球质心计算出风轮平面β，8在地球坐标系(Xe,Ye,Ze)中的方向向量，进一步的计算所述方向向量与东北天坐标系(ENU)中Y轴之间形成的偏航角。
在本实施例中，三个叶片中任意相邻的两叶片之间的角度为120度。
在本实施例中，在步骤S3和步骤S4之间还包括如下步骤：
-让无人机继续飞行，当再次检测到两个叶片完全重叠时，读取此时位置传感器获取的点P2的第二位置信息，通过点P2的第二位置信息对点P2的第一位置信息进行验证，从而提高了算法的效率。
所述无人机上搭载位置传感器、图像传感器和机载电脑；所述位置传感器、所述图像传感器连接所述机载电脑；
当所述无人机绕风机飞行时，所述位置传感器用于实时读取无人机位置信息，图像传感器用于拍摄风机叶片生成风机叶片图像，机载电脑用于无人机位置信息和风机叶片图像的处理。
在本实施例中，所述图像传感器工业相机。
根据叶片在不同视觉角度呈现不同姿态来精确估计P1，P2，结合地球质点三个不共线的位置点确定风轮平面β，即可求得偏航角a。同理转角的测量是先定位PT，同时读取PT对应图像，应用视觉图像检测出叶片姿态的方位角度。
图2为本发明中通过无人机进行风机偏航角测算方法的原理示意图，如图2所述，无人机绕风机轮毂飞行一周形成第一平面δ和飞行轨迹曲线s，如图2所示，第一平面δ与风轮平面β相交于直线l；直线l与飞行轨迹曲线s相交于点P1、P2。
由于点P1、P2在风轮平面β上，因此在确定点P1、P2后配合地球质心便能够确定风轮平面β。
所述无人机在绕风机轮毂飞行时，所述图像传感器采集叶片的视频流，所述位置传感器采集所述视频流对应的位置信息。
由于现在的水平轴大型风力发电机组多采用三叶片形式，根据平面视角遮挡原理，当无人机正好位于点P1或点P2时，会有一个叶片被其他两个叶片遮挡，所述图像传感器检测到风机叶片图像为两个叶片时，进一步的根据点P1、P2的位置特殊性，应用视觉跟踪方法能够将点P1、P2标定出来。
无人机在飞行时实时读取所述图像传感器拍摄到视频流fi，并对图像视频流fi进行预处理生成仅含叶片目标的二值图像流ti。
当无人机接近点P1或点P2时，三叶片中有两叶片接近重叠或者一叶片被部分遮挡，当其重叠率达到最大或者只能检测到两个叶片时，图像传感器检测到二值图像流ti近似为一条斜方向窄带，当无人机位于点P1或P2时，所述窄带的宽度为最小，即二值图像流ti中目标行数累加值τ为最小。
P1＝P[min(τ)] (1-1)
Figure GDA0002311521400000061
其中，τ为二值图像流ti中目标行数累加值，P为无人机的实时位置，P1为感兴趣位置，fi表示图像传感器采集的视频图像流，τ为根据ti(x,y)的值累加生成，根据式(1-2)，当ti(x,y)＝1时累加一次；
其中x表示二值图像流ti的x轴坐标值；y表示二值图像流ti的y轴坐标值。
因直线l与飞行轨迹曲线s相交于点P1、P2，即点P1、P2相对于轮毂存在对称关系，当先计算出点P1位置，便可计算出点P2大致位置，再借助无人机去验证点P2，从而进一步提高算法的效率。
当进行P0、P1位置验证时包括如下步骤：
步骤M1：将点P0、P1、P2的位置信息转换至地球坐标系(Xe,Ye,Ze)；在本实施例中，所述位置传感器为GPS模块，点P0、P1、P2的位置信息通过GPS模块，采用经度、纬度、高度表示；
转换计算公式为：
Figure GDA0002311521400000071
N为纬度B处的卯酉圈曲率半径，E为地球第一偏心率，
Figure GDA0002311521400000072
E＝(a2-b2)/a2,a为地球长半径，b为地球短半径，B为位置信息中的纬度，L为位置信息中的风塔经度，H为位置信息中的风塔高度；
步骤M2：验证点P2、P1在地球坐标的下的位置关系，即
Figure GDA0002311521400000073
Figure GDA0002311521400000074
其中
Figure GDA0002311521400000075
为P2，P1点直线距离，
Figure GDA0002311521400000076
为P1离风轮中心距离，
Figure GDA0002311521400000077
为P2离风轮中心距离；
步骤M3：计算精确度ratio,判断精确度ratio是否满足98％<ratio<102％；
Figure GDA0002311521400000078
图3为本发明中风轮平面β的示意图，图4为本发明中偏航角示意图，如图3、4所示，基于点P1、P2以及地球质心计算出风轮平面β，即可得到风机偏航在地球坐标系(Xe,Ye,Ze)中的方向向量，进一步的计算所述方向向量与Ye轴的偏航角；
为后续图像识别技术计算出精确的转角，应当在图像传感器位于叶轮正前方时读取叶轮的信息，所以点PT唯一。根据已经确定出的点P0、P1的位置信息，进一步的选择无人机拍摄路径曲线s中PT点，结合点P0的位置信息和风塔高度即可得到PT位置，然后读取PT点位置图像传感器拍摄的图像信息。
图5为本发明中通过视觉检测方法计算转角γ的示意图，如图5所示，本发明提供的通过无人机进行风机偏航角测算方法，还包括如下转角γ的测算步骤：
步骤A1：让无人机位于风机正前方，距离风塔底端设定距离的点P0位置以垂直升空至风塔高度，获取点PT的位置信息；
步骤A2：读取图像传感器于点PT采用的叶轮图像，将所述叶轮图像去除噪声；
步骤A3：对叶轮图像进行边缘检测，检测出目标叶片信息，通过角点检测计算叶尖点坐标，计算三个叶片几何中心点即为风轮中心在图像中的坐标Pwind centre；
步骤A4：连接叶尖点坐标与Pwind centre坐标确定一目标直线，进一步的计算该目标直线的斜率,即可得到该叶片的转角大小；
其中点PT表示无人机升空至所述风塔高度时的无人机的位置点。
图6为本发明中通过无人机进行风机偏航角测算系统的模块示意图，如图6所示，本发明提供的通过无人机进行风机偏航角测算系统100，包括：
飞行控制模块101，用于控制无人机以风塔高度绕风机飞行，当无人机在飞行过程中，通过图像传感器采集叶片的视频流；
重叠度计算模块102，用于对所述视频流中的叶片进行检测，当检测到风机的三个叶片时，对三个叶片进行实时跟踪，并实时计算三个叶片的相对位置及重叠度；
位置信息获取模块103，用于当检测到两个叶片完全重叠时，认定此时无人机飞行到风轮平面β上，读取此时位置传感器获取的点P1的位置信息；
位置信息计算模块104，用于根据点P1的位置信息计算与点P1以风塔呈轴对称分布的点P2的第一位置信息；
偏航角计算模块105，用于根据点P1的位置信息、点P2的第一位置信息以及地球质心计算出风轮平面β，进而计算出所述风轮平面的偏航角。
当实现本发明提供的通过无人机进行风机偏航角测算方法时，能够同时实现风机偏航角和转角γ的测量，包括如下步骤
步骤S1：记录位于风塔底端点P0位置的位置信息；
步骤S2：让无人机以风塔高度绕风机自动飞行一周形成飞行轨迹曲线s，飞行轨迹曲线s半径为叶片半径的整数倍，根据P0点位置和风塔高度确定飞行高度。
步骤S3：无人机在飞行过程中，通过图像传感器采集的视频流进行叶片检测，当检测到风机的三个叶片时，对三个叶片进行实时跟踪，并实时计算三个叶片的相对位置及重叠度；
步骤S4：当检测到两个叶片完全重叠时，认定此时无人机飞行到风轮平面β上，读取此时位置传感器获取的点P1的位置信息；
步骤S5：根据点P1的位置信息计算与点P1以风塔呈轴对称分别的点P2的第一位置信息；
步骤S6：让无人机继续飞行，当再次检测到两个叶片完全重叠时，读取此时位置传感器获取的点P2的第二位置信息，通过P2的第二位置信息对点P2的第一位置信息进行验证。
先在风塔下记录GPS点P0，并已知风塔的高度和叶片的长度；控制无人机飞起，高度为风塔高度，绕风轮中心点距以两倍叶片半径距离飞行；飞行到风机左侧，并计算风机两个叶片重叠率，判定是否位置在风机平面上，根据式(1-1)，得到该点的位置，即为点P1；确定P1点，机载电脑根据位置对称性，得到估计P2点估计位置，无人机飞行到风机右侧平面估计位置后继续飞行，即可获取P2点；通过P1、P2及地球质心求取风机偏航角。读取点P1(lat:40.17208455248887,lon:107.27228840493933)；读取点P2(lat:40.175627519375055,lon:107.27312725768884)；读取Pwind centre(lat:40.17421597362377,lon:107.27278682293934)；由公式1-3将P1、P2、Pwind centre点的经纬度位置信息转换成直角坐标系信息。得到
Figure GDA0002311521400000091
Figure GDA0002311521400000092
据公式1-4得ratio＝99.9975％，满足98％<ratio<102％。
本发明能够通过无人机测算风轮平面的方向向量与东北天坐标系Y轴之间形成的偏航角，从而能够根据风向调整风轮转轴，便于提高了风机的发电效率；本发明能够通过无人机测算出风轮平面的偏航角，从而能够便于风机进行建模分析，为实现风机的综合检测提供了便利；本发明中在无人机一次飞行过程中，不仅能够实现偏航角的测量，而且能够实现叶片转角的大小测量，提高了测量效果。
以上对本发明的具体实施例进行了描述。需要理解的是，本发明并不局限于上述特定实施方式，本领域技术人员可以在权利要求的范围内做出各种变形或修改，这并不影响本发明的实质内容。

Claims (8)

1.一种通过无人机进行风机偏航角测算方法，所述风机包括风塔和设置在风塔顶端的叶轮，所述叶轮包括轮毂和三个沿轮毂周向均匀分布的叶片，其特征在于，包括如下步骤：
步骤S1：控制无人机以风塔高度绕风机飞行，当无人机在飞行过程中，通过图像传感器采集所述叶轮的视频流；
步骤S2：对所述视频流中的叶片进行检测，当检测到风机的三个叶片时，对三个叶片进行实时跟踪，并实时计算三个叶片的相对位置及重叠度；
步骤S3：当检测到两个叶片完全重叠时，认定此时无人机飞行到风轮平面β上，读取此时位置传感器获取的点P1的位置信息；
步骤S4：根据点P1的位置信息计算与点P1以风塔呈轴对称分布的点P2的第一位置信息；
步骤S5：根据点P1的位置信息、点P2的第一位置信息以及地球质心计算出风轮平面β，进而计算出所述风轮平面的偏航角。
2.根据权利要求1所述的通过无人机进行风机偏航角测算方法，其特征在于，在步骤S3和步骤S4之间还包括如下步骤：
-让无人机继续飞行，当再次检测到两个叶片完全重叠时，读取此时位置传感器获取的点P2的第二位置信息，通过点P2的第二位置信息对点P2的第一位置信息进行验证。
3.根据权利要求1所述的通过无人机进行风机偏航角测算方法，其特征在于，所述检测到两个叶片完全重叠时的点P1通过如下方式计算得出：
P1＝P[min(τ)] (1-1)
Figure FDA0002311521390000011
其中，τ为二值图像流ti中目标行数累加值，P为无人机的实时位置，目标行数累加值τ为根据ti(x,y)的值累加生成，根据式(1-2)，当ti(x,y)＝1时累加一次；
当目标行数累加值τ为最小时，确定所述两个叶片完全重叠；
其中x表示二值图像流ti的x轴坐标值；y表示二值图像流ti的y轴坐标值。
4.根据权利要求1所述的通过无人机进行风机偏航角测算方法，其特征在于，还包括如下点P1、P2的位置验证步骤：
步骤M1：将点P1的位置信息、P2的第一位置信息转换至地球坐标系(Xe,Ye,Ze)，点P1、P2的位置信息通过GPS模块，采用经度、纬度、高度表示，转换计算公式为：
Figure FDA0002311521390000021
N为纬度B处的卯酉圈曲率半径，E为地球第一偏心率，
Figure FDA0002311521390000022
E＝(a2-b2)/a2,a为地球长半径，b为地球短半径，B为位置信息中的纬度，L为位置信息中的风塔经度，H为位置信息中的风塔高度；
步骤M2：验证点P2、P1在地球坐标的下的位置关系，即
Figure FDA0002311521390000023
Figure FDA0002311521390000024
其中
Figure FDA0002311521390000025
为点P2与点P1之间的直线距离，
Figure FDA0002311521390000026
为点P1离风轮中心距离，
Figure FDA0002311521390000027
为点P2离风轮中心距离；
步骤M3：计算精确度ratio,判断精确度ratio是否满足98％<ratio<102％；
Figure FDA0002311521390000028
5.根据权利要求1所述的通过无人机进行风机偏航角测算方法，其特征在于，还包括如下转角γ的测算步骤：
步骤N1：让无人机位于风机正前方，距离风塔底端设定距离的点P0位置以垂直升空至风塔高度，获取点PT的位置信息；
步骤N2：读取图像传感器于点PT采用的叶轮图像，将所述叶轮图像去除噪声；
步骤N3：对叶轮图像进行边缘检测，检测出目标叶片信息，通过角点检测计算叶尖点坐标，计算三个叶片几何中心点即为风轮中心在图像中的坐标Pwind centre；
步骤N4：连接叶尖点坐标与Pwind centre坐标确定一目标直线，进一步的计算该目标直线的斜率,即可得到该叶片的转角大小；
其中点PT表示无人机升空至所述风塔高度时的无人机的位置点。
6.根据权利要求1所述的通过无人机进行风机偏航角测算方法，其特征在于，所述无人机上搭载位置传感器、图像传感器和机载电脑；
所述位置传感器、所述图像传感器连接所述机载电脑；所述位置传感器用于实时读取无人机位置信息，图像传感器用于拍摄风机叶片生成风机叶片图像，机载电脑用于无人机位置信息和风机叶片图像的处理。
7.根据权利要求1所述的通过无人机进行风机偏航角测算方法，其特征在于，所述步骤S5具体为：
-将点P1的位置信息、点P2的第一位置信息转换至地球坐标系，再根据转换至地球坐标系中点P1、点P2及地球质心计算出风轮平面β，进而可得到在地球坐标系中，风轮平面β在地球坐标系(Xe,Ye,Ze)中的方向向量，进一步计算所述方向向量与东北天坐标系中Y轴之间形成的偏航角。
8.一种通过无人机进行风机偏航角测算系统，其特征在于，用于实现权利要求1至7任一项所述的通过无人机进行风机偏航角测算方法，包括：
飞行控制模块，用于控制无人机以风塔高度绕风机飞行，当无人机在飞行过程中，通过图像传感器采集叶片的视频流；
重叠度计算模块，用于对所述视频流中的叶片进行检测，当检测到风机的三个叶片时，对三个叶片进行实时跟踪，并实时计算三个叶片的相对位置及重叠度；
位置信息获取模块，用于当检测到两个叶片完全重叠时，认定此时无人机飞行到风轮平面β上，读取此时位置传感器获取的点P1的位置信息；
位置信息计算模块，用于根据点P1的位置信息计算与点P1以风塔呈轴对称分布的点P2的第一位置信息；
偏航角计算模块，用于根据点P1的位置信息、点P2的第一位置信息以及地球质心计算出风轮平面β，进而计算出所述风轮平面的偏航角。
CN201810627916.8A 2018-06-19 2018-06-19 通过无人机进行风机偏航角测算方法及系统 Active CN108843492B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810627916.8A CN108843492B (zh) 2018-06-19 2018-06-19 通过无人机进行风机偏航角测算方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810627916.8A CN108843492B (zh) 2018-06-19 2018-06-19 通过无人机进行风机偏航角测算方法及系统

Publications (2)

Publication Number Publication Date
CN108843492A CN108843492A (zh) 2018-11-20
CN108843492B true CN108843492B (zh) 2020-04-21

Family

ID=64202144

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810627916.8A Active CN108843492B (zh) 2018-06-19 2018-06-19 通过无人机进行风机偏航角测算方法及系统

Country Status (1)

Country Link
CN (1) CN108843492B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110687560A (zh) * 2019-10-12 2020-01-14 中国电建集团贵州工程有限公司 一种叶轮组合快速定位方法
CN115143056B (zh) * 2022-07-21 2023-06-16 中国长江三峡集团有限公司 风力发电机的停机姿态参数测量方法及装置
CN115480589B (zh) * 2022-09-06 2023-07-25 中科云尚(南京)智能技术有限公司 基于无人机的风机巡检航线生成方法及系统
CN117212077B (zh) * 2023-11-08 2024-02-06 云南滇能智慧能源有限公司 一种风电机的风轮故障监测方法、装置、设备及存储介质

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4589633B2 (ja) * 2004-01-29 2010-12-01 富士重工業株式会社 水平軸風車及び水平軸風車の制御方法
DK2599993T3 (en) * 2011-12-01 2016-06-06 Siemens Ag A method for determining krøjningsvinklen of a component of a wind turbine
US9562515B2 (en) * 2013-02-01 2017-02-07 General Electric Company Method and apparatus for wind turbine noise reduction
CN104948384B (zh) * 2015-06-16 2019-01-29 新疆金风科技股份有限公司 风力发电机组偏航角度校准方法和系统
CN107084723B (zh) * 2017-05-12 2019-07-02 中国人民解放军91550部队 一种海洋环境下的水下航行体运动轨迹估计方法
CN107884781B (zh) * 2017-11-07 2020-06-02 北京电子工程总体研究所 一种双无人机跟踪测距方法
CN107956652A (zh) * 2017-11-14 2018-04-24 青岛华创风能有限公司 一种风力发电机组偏航测量方法及装置

Also Published As

Publication number Publication date
CN108843492A (zh) 2018-11-20

Similar Documents

Publication Publication Date Title
CN108733079B (zh) 通过无人机对风机进行自动巡检飞行路径确定方法及系统
CN108843492B (zh) 通过无人机进行风机偏航角测算方法及系统
CN108894933B (zh) 通过无人机对风机叶尖跟踪检测时跟丢再捕捉方法及系统
CN108915959B (zh) 通过无人机对风机叶尖区域绕行检测方法及系统
CN108869197B (zh) 通过无人机进行风机高度精确测量方法及系统
CN110631593B (zh) 一种用于自动驾驶场景的多传感器融合定位方法
CN108869196B (zh) 通过无人机对风机后侧叶根区域检测方法及系统
US10703479B2 (en) Unmanned aerial vehicle, control systems for unmanned aerial vehicle and control method thereof
CN105484938B (zh) 风力发电机组的偏航控制方法及装置
CN109031312B (zh) 适用于烟囱内部作业的飞行平台定位装置和定位方法
CN110475966A (zh) 确定风力涡轮机的转子平面的定向
CN108204879B (zh) 一种转动惯量的光学测量方法及系统
CN101532841A (zh) 基于地标捕获跟踪的飞行器导航定位方法
CN106885918B (zh) 一种面向多旋翼飞行器的多信息融合实时风速估计方法
CN110426178B (zh) 一种基于尾流示踪的风场测量方法及系统
CN106155081A (zh) 一种旋翼无人机大范围目标监测和精确定位方法
CN113885580A (zh) 基于无人机实现自动化巡检风机的路径规划方法及系统
CN108225273B (zh) 一种基于传感器先验知识的实时跑道检测方法
CN104504255A (zh) 一种螺旋翼升力和阻力力矩的确定方法
Guo et al. Detecting and positioning of wind turbine blade tips for uav-based automatic inspection
CN112577606A (zh) 一种双无人机搭载主动热成像的风机叶片巡检方法
CN110007318B (zh) 风场干扰下基于卡尔曼滤波的单无人机判断gps欺骗的方法
Wildmann et al. Towards vertical wind and turbulent flux estimation with multicopter uncrewed aircraft systems
CN106408589A (zh) 基于车载俯视相机的车辆运动测量方法
CN109764864A (zh) 一种基于颜色识别的室内无人机位姿获取方法及系统

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant