基于GNSS-RTK技术的风电机组机舱位移和轨迹研究_刘智超.pdf

2023.01.DQGY 19CHINA ELECTRICAL EQUIPMENT INDUSTRYPRODVCT AND TECHNIC产品与技术基准站的信号，通过对所收到的信号进行实时处理，给出移动站的三维坐标，如图1所示，RTK技术可以在很短的时间内获得厘米级的定位精度8-10。图1 GNSS-RTK技术拓扑图基于GNSS-RTK技术的风电机组机舱位移测量的现场实施方案如图2所示，移动站通过RTK定位技术获得经纬度信息传递给风电机组主控，主控通过经纬度转位移计算模块以及前期的标定参数实时计算机舱偏航坐标系下X方向和Y方向的数据。0 引言风力发电机组塔架振动频率处于低频段，通过获得机舱（塔顶）低频绝对位移数据，可以对机组叶轮不平衡、机组稳定性评估、塔架载荷反演等方面提供有效数据支撑，无论是对机组稳定性控制还是优化仿真和技术研究具有重要的价值和意义。常规的测量手段很难获得有效、稳定和低成本的测量方案1-4，而GNSS-RTK定位技术利用全球导航卫星系统的RTK载波相位差分技术，实时处理基准站和移动站所得载波相位的差分结果，能够获得动态的厘米级定位精度5-7。1 技术说明1.1 宏观性能测试RTK（Real Time Kinematic）是以载波相位观测值进行实时动态相对定位的技术，其原理是将位于基准站上的接收机观测的卫星数据，通过数据通信链（如无线电台）实时发送出去，而位于附近的移动站（流动站）接收机在对卫星观测的同时，也接收来自刘智超 张鹏 姜禹含 齐男 李永亮 李满达吐（蒙东协合开鲁风力发电有限公司）基于GNSS-RTK技术的风电机组机舱位移和轨迹研究摘要：获得风力发电机组机舱低频绝对位移数据，可以对机组叶轮不平衡、机组稳定性评估、塔架载荷反演等方面提供有效数据支撑，无论是从机组稳定性控制还是优化仿真和技术研究具有重要的价值和意义。通过采用GNSS-RTK技术、特定标定方法、经纬度转位移方法和相关尺寸链解算公式计算获得了机舱绝对位移结果。通过现场实测数据与Bladed仿真论证了GNSS-RTK技术在风电机组上能够达到厘米级的绝对定位精度，通过分析工况数据、机组运行状态以及机舱定位定向数据，机舱运行和待机时的位移结果和工况相关性明显，说明了测量的机组绝对位移数据的有效性和准确性。关键词：风电机组；GNSS-RTK；机舱绝对位移电台天线电台已知位置基准站(GNSS接收机)流动站（待测点）1020km电器工业202301设计排版.indd 272023.2.6 11:18:53 AM2023.01.DQGY 20CHINA ELECTRICAL EQUIPMENT INDUSTRYPRODVCT AND TECHNIC产品与技术图2 实施方案示意图1.2 坐标转换说明机舱绝对位移计算原理图如图3所示。O点为塔底坐标系原点，对应塔底塔筒圆心，O1对应塔顶塔筒的圆心，忽略塔筒制造安装和重力影响，O点将与O1点重合，实际情况由于机舱重心不在塔筒轴线上，O1通常会向机舱前部倾斜一段距离，O2点为机组受到风载作用，机舱运动到达的随机位置对应的坐标点。点A、点B和点C位于测风支架之上，A和C位置分别安装两只天线，A为主天线，C为副天线，两天线可以测量AC指向，可用于向量解析坐标数据，B位置为测风支架中心，AO1指O1状态下A点位置，其他坐标以此类推11。图3 坐标转换关系图A坐标的变化即为测量点的轨迹变化，OO2的变化即为机舱的绝对位移变化，计算关系如式（1）所示：（1）地球椭球的第一偏心率公式如式（2）所示：（2）地球椭球的卯酉圈曲率半径公式如式（3）所示：（3）AO2在塔筒坐标系下距离坐标（Xi,Yi）计算公式如式（4）和式（5）所示：（4）（5）AB距离（测点和测风支架中心距离）通过现场实测获得，O1Bo1尺寸通过机组设计图纸获得，通过坐标向量转换关系可以获得O1点距离坐标O1（X1，Y1），O2（X2，Y2）点坐标同理可以获得。2 标定说明2.1 安装部署安装部署机组基本参数信息如表1所示，测量天线安装实物图如图4所示，右边为主天线，左边为副天线。表1 测试机组基本信息信息数值塔高120m叶轮直径140m容量2.5MW叶片长度69m额定转速11.1rpm图4 现场移动站和基站天线实物图Co2Bo2Ao2Bo1Ao1Co1O2O1OsinYiXi电器工业202301设计排版.indd 282023.2.6 11:18:55 AM2023.01.DQGY 21CHINA ELECTRICAL EQUIPMENT INDUSTRYPRODVCT AND TECHNIC产品与技术2.2 安装部署卫星定位定向系统安装在机舱顶部测风支架上，可输出航向角、经度、纬度等，通过标定与换算，可实现塔顶绝对位移的测量。标定的意义在于通过机舱在无风或小风工况下，机舱偏航若干周绘制出定位天线的轨迹，进而计算出塔架实际轴线经纬度坐标，作为后续卫星定位数据处理过程中的参考坐标。绘制出标定过程机组经纬度变化时序图，如图4所示。以经度为横坐标，纬度为纵坐标，绘制出经纬度椭圆，如图5所示。图5 标定全过程经纬度变化情况通过上述椭圆拟合和迭代求得塔顶中心点O点最小误差的经纬度，如表2所示。如前所述，O点可以理解为塔架无内应力（塔底无弯矩载荷）状态下塔顶所在的位置，卫星定位系统可同时测量出机组所处的海拔高度。表2中的数据将直接应用于算法程序中。表2 塔顶中心点O点经纬度参数数值纬度35.346347238245720经度114.2086539510335机舱定位天线海拔高度191.86m通过机舱设计图纸和现场安装实测得到O1Bo1尺寸和Ao1Bo1尺寸，通过式（1）获得OAo1测量天线轨迹和OO1机舱轨迹，如图6所示。（a）测量天线轨迹图 （b）机舱位移轨迹图图6 测量天线和机舱轨迹图（a）机舱绝对位移幅值图 （b）机舱绝对位移分布图图7 机舱绝对位移幅值和分布直方图如图7所示，机舱在无风、小风条件下连续偏航，通过GNSS-RTK测量的机舱位移均值为0.138m，最大值0.22m，最小值0.07m，均方根0.028，可以观察出，偏航过程中机舱存在晃动，晃动幅度大多数数据在0.05m，这应该是风载和偏航系统内部激励力等因素导致的。通过实测数据得到机舱由于重力作用向前移动了0.138m，通过机组Bladed静态仿真得到机舱由于重力原因会向前移动0.126m，距离实测值仅仅相差1.2cm，可见通过该测量方案可以得到机舱位移数据，精度可达到厘米级别。3 标定说明3.1 安装部署机组塔顶绝对位移是通过卫星主天线经纬度坐标向位移坐标转化，以及相关尺寸链解算而得到。通过分析工况数据、机组运行状态以及机舱定位定向数据电器工业202301设计排版.indd 292023.2.6 11:18:56 AM2023.01.DQGY 22CHINA ELECTRICAL EQUIPMENT INDUSTRYPRODVCT AND TECHNIC产品与技术进而分析机组定位数据的有效性和稳定性等信息。如图8所示，分别展示了卫星测量天线轨迹、机舱绝对位移轨迹，机舱航向是以机舱机尾指向机头为X正方向，机舱右侧指向左侧为Y方向正方向（右手坐标系）。待机时机舱偏航坐标Y方向位移与风速之间的关系。图9 发电时机舱X方向位移与风速之间的关系图10 待机时机舱X方向位移与风速之间的关系可以得到以下规律和结论：1）机组发电时，随着风速增加，在9m/s风速时，机舱X方向位移达到最大值，均值为-0.561m，0.1和0.9分位数分别为-0.674m和-0.374m，随着风速递增，机舱位移最大值减小，这和机组控制策略和出力方式有关。2）机组待机时，机舱X方向位移变化不大，基本在0.08m左右前后微微晃动，按照标定时机舱前倾0.14m来计算，风载使得机舱后移了0.06m。3）机组待机时，机舱Y方向位移变化不大，基本在0左右晃动。4）机组发电时，机舱Y方向位移在随着风速从5m/s上升，机舱晃动位置在0.1m左右，即来风方向是叶轮的右前方，推断为机组对风偏差或地形等原因导致。（a）卫星测量天线轨迹图 （b）机舱轨迹图（c）机舱轨迹图-机舱方向 （d）机舱轨迹坐标时序图图8 机舱运动情况下数据分析机舱在风载作用下，机舱X方向位移基本为负值，包括小风和停机在内，机舱运动范围在-1.00.3m之间。其中在数据中可以明显看出一次停机过程，机舱从-0.6m急速移动到0.3m后振荡稳定下来。通过以上分析可以得出GNSS-RTK能够很好测量机舱绝对位移，可以达到厘米级精度。3.2 相关性分析当机组机舱正对风向时，机舱受到向后的推力，尤其是机组并网发电时，机舱前后晃动的幅度很大，如图9和图10所示，展示了现场机组发电和待机时机舱偏航坐标X方向位移与风速之间的关系，蓝色散点为实测值，三角标记分别为0.1和0.9分位点，红色圆点为均值。如图11和图12所示，展示了现场机组发电和电器工业202301设计排版.indd 302023.2.6 11:18:58 AM2023.01.DQGY 23CHINA ELECTRICAL EQUIPMENT INDUSTRYPRODVCT AND TECHNIC产品与技术图11 发电时机舱Y方向位移与风速之间的关系图12 待机时机舱Y方向位移与风速之间的关系4 结束语本文提出一种基于GNSS-RTK技术的风电机组机舱位移和轨迹的研究方法，阐明了标定的基本方法和应用案例，通过仿真和实测数据的维度论证了技术方案的有效性和准确性，可以得出以下结论。1）建立测点位置和机舱中心坐标尺寸链关系，通过向量转换能够将经纬度信息转化为坐标信息，进而转换为机舱绝对位移信息，通过小风偏航测试活动就能标定实现。2）机组小风偏航过程中，机舱中心位置变化的Bladed仿真结果和实测结果偏差仅为1.2cm，完全能够实现厘米级的绝对定位精度。3）通过分析工况数据、机组运行状态以及机舱定位定向数据，机舱运行和待机时的位移结果和工况相关性明显，说明了测量的机组绝对位移数据的有效性和准确性。参考文献1 王潇,俞佳男,唐波,等.基于GNSS的风电机组低频振动位移计算方法研究J.机电工程,2021,38(8):996-1002，1057.2 刘国祥.多兆瓦风电机组独立变桨系统控制策略的研究与应用D.乌鲁木齐：新疆农业大学,2013.3 巫发明,杨从新,王青,等.大型风电机组塔筒动态载荷识别研究与验证J.振动与冲击,2020,39(20):197-203.4 吕鹏.时变载荷作用下大型风力机塔筒动态特性分析D.沈阳：沈阳工业大学,2014.5 罗虹,赵政权.高度角对GNSS多系统组合短基线RTK影响J.导航定位学报,2022,10(4):145-152.6 程建华,崔珈浩,贾春.GNSS-RTK/UWB紧组合定位模型性能分析与提升C.第十三届中国卫星导航年会论文集S03卫星导航系统与增强,2022:123-130.7 黄会宝,高强.基于低成本GNSS接收机的多天线RTK定位研究J.地理空间信息,2022,20(1):94-97，9.8 付宁波.GNSS-RTK测量精度与重复性精度检测分析J.江西测绘,2021(4):5-7.9 Olivart i Llop J.M.,MorenoSalinas D.,Snchez J.Full Real-Time Positioning and Attitude System Based on GNSS-RTK TechnologyJ.Sustainability,2020,12(23).10 Xiong Chunbao,Niu Yanbo,Wang Zhilin,Yuan Linlin.Dynamic monitoring of a super high-rise structure based on GNSS-RTK technique combining CEEMDAN and wavelet threshold analysisJ.European Journal of Environmental and Civil Engineering,2019,25(10).11 于亮亮.确定风力发电机组的塔顶的振动位移的方法和设备:中国，201810474365.6 P.2019-11-26.12 杨红雷.融合多类卫星GNSS/SLR数据的精密轨道确定及大地测量参数解算D.济南：山东大学,2021.（收稿日期：2022-11-03）电器工业202301设计排版.indd 312023.2.6 11:18:59 AM