高级检索

ISSN1001-3806CN51-1125/TN 网站地图

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

大型风电叶片动态摄影测量网络优化

丰伟 董明利 孙鹏

引用本文:
Citation:

大型风电叶片动态摄影测量网络优化

    作者简介: 丰伟(1993-),男,硕士研究生,现主要从事摄影测量的研究.
    通讯作者: 董明利, dongml@bistu.edu.cn
  • 基金项目:

    北京市教委科技计划重点资助项目 KZ201711232029

  • 中图分类号: TM315;P234

Dynamic photogrammetry network optimization for large wind turbine blades

    Corresponding author: DONG Mingli, dongml@bistu.edu.cn ;
  • CLC number: TM315;P234

  • 摘要: 在大型风电叶片动态摄影测量中,为了对相机的站位进行优化,采用一种变异操作改进型遗传算法作为摄影测量网络优化方法,通过光线束前方交会的误差传递建立测量误差模型,以空间坐标测量误差的标准差为网络优化的目标,同时根据被测风电叶片几何结构和实际环境确定了相应的约束条件进行仿真实验,得到了最优的相机站位。结果表明,在以叶片长度为40m的风机为被测物的仿真实验中,最优站位的空间坐标测量误差标准差为2.7mm;通过对叶片长度为3.5m的风机模型进行实测实验验证,最优站位的相对测量误差为0.009%,最大误差为0.617mm。该研究为风电叶片摄影测量的网络优化提供了参考。
  • Figure 1.  Camera stations and large wind turbine

    Figure 2.  Flow chart of network optimization algorithm

    Figure 3.  The angle between the optical axis of the camera and the plane of the wind wheel

    Figure 4.  Initial set of stations

    Figure 5.  Variation of the standard deviation of the optimal station and the mean standard deviation of all stations

    Figure 6.  Variation of the X-coordinate of the optimal station with iterations

    Figure 7.  Variation of the Z-coordinate of the optimal station with iterations

    Figure 8.  Experimental site and the turbine to be measured

    Table 1.  The standard deviation and coordinates of the optimal station in iteration

    iterations the spatial standard deviation of the optimal station/mm X-coordinate/ mm Z-coordinate/ mm
    1 4.100 8×104 2.1×105
    10 3.027 6.353×104 1.467×105
    20 2.884 6.536×104 1.418×105
    30 2.749 7.258×104 1.406×105
    40 2.741 7.551×104 1.421×105
    50 2.741 7.551×104 1.421×105
    60 2.741 7.370×104 1.409×105
    70 2.741 7.352×104 1.408×105
    80 2.740 7.352×104 1.408×105
    90 2.740 7.352×104 1.408×105
    100 2.740 7.352×104 1.408×105
    下载: 导出CSV

    Table 2.  Comparison of measurement results between optimal and stochastic stations

    station number maximum error/mm average error/mm root mean square/mm relative error/%
    1(optimal station) 1.392 0.617 0.836 0.009
    2 2.241 1.273 1.562 0.018
    3 3.277 2.852 2.872 0.041
    4 5.639 3.057 3.568 0.051
    5 2.715 2.556 2.564 0.037
    6 3.415 2.909 2.931 0.042
    7 2.685 2.036 2.097 0.029
    8 3.415 2.523 2.640 0.036
    9 4.382 3.142 3.257 0.045
    10 5.149 3.127 3.496 0.045
    下载: 导出CSV
  • [1]

    XU J J, CHEN Y N. The installed capacity of wind power in China is the first in the world for four consecutive years. The problems and countermeasures for the development of wind power industry to be solved[J]. China Economic Weekly, 2016(16): 80-81(in Chin-ese).
    [2]

    GRIFFITH D T, MAYES R L, HUNTER P S. Excitation methods for a 60kW vertical axis wind turbine[M].New York, USA: Springer, 2011: 329-338.
    [3]

    OSGOOD R, BIR G, MUTHA H, et al. Full-scale modal wind turbine tests: Comparing shaker excitation with wind excitation[M].New York, USA: Springer, 2011: 113-124.
    [4]

    KONG C, BANG J, SUGIYAMA Y. Structural investigation of composite wind turbine blade considering various load cases and fatigue life[J]. Energy, 2005, 30(11): 2101-2114.
    [5]

    TAMURA Y, MATSUI M, PAGNINI L C, et al. Measurement of wind-induced response of buildings using RTK-GPS[J]. Journal of Wind Engineering and Industrial Aerodynamics, 2002, 90(12/15): 1783-1793.
    [6]

    BREUER P, CHMIELEWSKI T, GÓRSKI P, et al. Application of GPS technology to measurement of displacement of high-rise structures due to weak winds[J]. Journal of Wind Engineering and Industrial Aerodynamics, 2002, 90(3): 223-230. doi: 10.1016/S0167-6105(01)00221-5
    [7]

    NICKITOPOULOU A, PROTOPSALTI K, STIROS S. Monitoring dynamic and quasi-static deformations of large flexible engineering structures with GPS: Accuracy, limitations and promises[J]. Engineering Structures, 2006, 28(10): 1471-1482. doi: 10.1016/j.engstruct.2006.02.001
    [8]

    CHEN Y, NI Y Q, YE X W, et al. Structural health monitoring of wind turbine blade using fiber Bragg grating sensors and fiber optic rotary joint[J]. Proceedings of the SPIE, 2012, 8345: 86-91.
    [9]

    SCHROEDER K, ECKE W, APITZ J, et al. A fibre Bragg grating sensor system monitors operational load in a wind turbine rotor blade[J]. Measurement Science and Technology, 2006, 17(5): 1167-1172. doi: 10.1088/0957-0233/17/5/S39
    [10]

    FRASER C S. Network design considerations for non-topographic photogrammetry[J]. Photogrammetric Engineering & Remote Sensing, 1984, 50(8): 1115-1126.
    [11]

    OLAGUE G. Automated photogrammetric network design using genetic algorithms[J]. Photogrammetric Engineering & Remote Sensing, 2002, 68(5): 423-432.
    [12]

    QIAO Y J, WANG H R, ZHAO Y J. Study on binocular vision measurement network layout for large curved surface parts[J]. Ch-inese Journal of Scientific Instrumen, 2015, 36(4): 913-918(in Chinese).
    [13]

    CHEN J Y, DONG M L, WANG J, et al. Photogrammetric network design for large-scale trough concentrator surface measurement[J]. Renewable Energy Resources, 2016, 34(3): 353-359(in Chin-ese).
    [14]

    LIU Y, XIAO Sh D, ZHANG R, et al. Initial estimation of digital image correlated deformation based on genetic algorithms[J]. Laser Technology, 2020, 44(1): 130-135(in Chinese).
  • [1] 曾亮维蔡元静谭础深黄佐华 . 斯托克斯光偏振态测量系统的优化. 激光技术, 2017, 41(1): 74-78. doi: 10.7510/jgjs.issn.1001-3806.2017.01.016
    [2] 魏继锋彭勇高学燕周山 . 简易型大直径激光束功率能量在线测量装置. 激光技术, 2011, 35(3): 322-325. doi: 10.3969/j.issn.1001-3806.2011.03.010
    [3] 龙龙李宗峰 . 基于激光位移传感器的3维位置测量算法研究. 激光技术, 2017, 41(4): 531-536. doi: 10.7510/jgjs.issn.1001-3806.2017.04.015
    [4] 范宜艳赵斌马国鹭 . 基于光学标靶与测距仪的隐藏区域坐标测量. 激光技术, 2014, 38(6): 723-728. doi: 10.7510/jgjs.issn.1001-3806.2014.06.001
    [5] 杜志广颜树华林存宝王国超魏春华 . 双路信号相位同步测量系统设计与实现. 激光技术, 2016, 40(3): 315-319. doi: 10.7510/jgjs.issn.1001-3806.2016.03.003
    [6] 肖兴维马国鹭曾国英陆野 . 正交视觉与倾角仪组合空间位姿测量方法研究. 激光技术, 2020, 44(3): 278-282. doi: 10.7510/jgjs.issn.1001-3806.2020.03.002
    [7] 朱红伟叶会英 . 光反馈自混合干涉系统反馈水平的研究与测量. 激光技术, 2010, 34(6): 847-850. doi: 10.3969/j.issn.1001-3806.2010.06.034
    [8] 张德斌宋余华王全胜杜亚清张新兴张豪 . 激光发散角测量的误差分析. 激光技术, 2016, 40(6): 926-929. doi: 10.7510/jgjs.issn.1001-3806.2016.06.031
    [9] 陈曼龙 . 机器视觉螺纹测量的误差分析. 激光技术, 2014, 38(1): 109-113. doi: 10.7510/jgjs.issn.1001-3806.2014.01.024
    [10] 梁金辉赵冬娥董娟 . 光电靶的设计与改进. 激光技术, 2008, 32(5): 456-459.
    [11] 彭敦云宋连科栗开婷郭文静 . 偏光干涉法测量液晶的双折射率. 激光技术, 2014, 38(3): 422-424. doi: 10.7510/jgjs.issn.1001-3806.2014.03.030
    [12] 刘顺涛骆华芬陈雪梅徐静 . 结构光测量系统的标定方法综述. 激光技术, 2015, 39(2): 252-258. doi: 10.7510/jgjs.issn.1001-3806.2015.02.023
    [13] 肖长江张景超魏勇李兴元胡学良 . 基于激光视觉原理测量玻璃中气泡的尺寸. 激光技术, 2015, 39(3): 391-394. doi: 10.7510/jgjs.issn.1001-3806.2015.03.024
    [14] 聂雪莹项飞荻黄欣刘劲松杨振刚王可嘉 . 金属平板的太赫兹雷达散射截面测量. 激光技术, 2016, 40(5): 676-681. doi: 10.7510/jgjs.issn.1001-3806.2016.05.012
    [15] 金岩华叶会英孙瑶 . 基于相位解卷的位移测量系统. 激光技术, 2009, 33(4): 359-361,365. doi: 10.3969/j.issn.1001-3806.2009.04.007
    [16] 杨初平刘建斌谭穗妍翁嘉文 . 应用频率积分相位解调测量径向畸变. 激光技术, 2014, 38(3): 402-405. doi: 10.7510/jgjs.issn.1001-3806.2014.03.026
    [17] 周策策李杏华 . 基于机器视觉的螺纹参量测量系统. 激光技术, 2016, 40(5): 643-647. doi: 10.7510/jgjs.issn.1001-3806.2016.05.006
    [18] 王国超颜树华高雷谢学东田震 . 光栅干涉位移测量技术发展综述. 激光技术, 2010, 34(5): 661-664,716. doi: 10.3969/j.issn.1001-3806.2010.O5.023
    [19] 李建民王蕴芬田野牛振凤刘伟东韩冰刘钰马艳丽 . 基于远场干涉测量棱镜内气泡直径. 激光技术, 2010, 34(1): 67-70. doi: 10.3969/j.issn.1001-3806.2010.01.019
    [20] 胡林亭史德民李佩军任成才 . 激光监测系统测量精度的检测方法. 激光技术, 2008, 32(6): 670-672.
  • 加载中
图(8) / 表(2)
计量
  • 文章访问数:  435
  • HTML全文浏览量:  361
  • PDF下载量:  1
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-03-12
  • 录用日期:  2020-04-07
  • 刊出日期:  2021-01-25

大型风电叶片动态摄影测量网络优化

    通讯作者: 董明利, dongml@bistu.edu.cn
    作者简介: 丰伟(1993-),男,硕士研究生,现主要从事摄影测量的研究
  • 1. 北京信息科技大学 光电测试技术及仪器教育部重点实验室,北京 100192
  • 2. 北京信息科技大学 机电系统测控北京市重点实验室, 北京 100192
基金项目:  北京市教委科技计划重点资助项目 KZ201711232029

摘要: 在大型风电叶片动态摄影测量中,为了对相机的站位进行优化,采用一种变异操作改进型遗传算法作为摄影测量网络优化方法,通过光线束前方交会的误差传递建立测量误差模型,以空间坐标测量误差的标准差为网络优化的目标,同时根据被测风电叶片几何结构和实际环境确定了相应的约束条件进行仿真实验,得到了最优的相机站位。结果表明,在以叶片长度为40m的风机为被测物的仿真实验中,最优站位的空间坐标测量误差标准差为2.7mm;通过对叶片长度为3.5m的风机模型进行实测实验验证,最优站位的相对测量误差为0.009%,最大误差为0.617mm。该研究为风电叶片摄影测量的网络优化提供了参考。

English Abstract

    • 风能作为一种绿色环保的可再生能源,是世界各国能源发展的重要方向。在过去十几年,我国的风电装机量已位居世界第一[1]。风电叶片作为风机的重要部件,其空气动力学特性对风机的发电效率和使用寿命有重要影响。为使风机高效安全运行,需要测量运行过程中叶片的运动和形变数据来获得其空气动力学特性。

      风电叶片测量方法很多,但这些测量方法多为接触式测量,包括通过在叶片上布设应变片、全球定位系统(global positioning system,GPS)传感器或光纤传感器的方式测量叶片的运动和形变[2-9]。上述测量方法技术成熟,测量精度高,但接触式测量方法需要布设大量的传感器与电缆,使测量过程复杂并增加了测量成本。同时,风机多运行于恶劣环境中,叶片表面的传感器易受环境影响,降低测量精度。摄影测量是一种非接触测量技术,具有精度高、操作简单和成本低等优点,将摄影测量技术应用于风电叶片的测量有着很好的发展前景。

      在摄影测量中,为了获得叶片运动的高精度3-D测量数据,需要建立由多台相机组成的摄影测量网络。由于风电叶片位于几十米的高空中,尺寸巨大,而测量相机位于地面,站位数量少、约束条件多,不同站位的测量精度差别很大。为提高风电叶片的测量精度,需要对摄影测量网络进行优化。

      网络优化作为摄影测量领域重要的研究内容,在国内外都展开了广泛研究。1984年,加拿大的FRASER已经提出了非地形摄影测量的网络优化问题,并分析了网络优化中的约束条件[10]。2002年,OLAGUE提出了使用遗传算法实现自动摄影测量网络优化[11]。相比国外,国内相关研究起步较晚,多为应用方面的研究及算法改进。2015年,哈尔滨理工大学的QIAO等人采用改进的多目标并列选择遗传算法对大尺寸曲面零件双目视觉测量网络的规划进行了研究[12],兼顾分辨率与覆盖率两个优化目标,有效地解决了双目视觉网络规划问题。2016年,北京信息科技大学的CHEN等人采用遗传算法针对槽式聚光器面形摄影测量进行网络规划[13],得到了最优的相机布局,提高了测量精度。上述优化方法测量时可采用的相机站位多,站位布设灵活,但这些方法都只能应用于尺寸不超过10m的静止物体的静态摄影测量,对于相机站位只有2个,站位的约束条件多, 被测物尺寸巨大且处于运动状态的大型风电叶片动态摄影测量并不适用,因此,需要针对风电叶片动态摄影测量中相机站位少、约束条件多等因素,设计其网络优化方法。

      本文中研究了大型风电叶片动态摄影测量的网络优化方法,该方法通过3-D重构过程的误差传递建立误差模型,以空间坐标测量误差的标准差作为优化目标,测量时的相机站位为优化变量,通过一种变异改进型遗传算法求风电叶片摄影测量中相机的最优站位[14]

    • 大型风电叶片的摄影测量系统如图 1所示。

      Figure 1.  Camera stations and large wind turbine

      图 1中,A点代表风机的旋转中心,O点为风机塔柱底部,C, D点代表相机站位,C, D两点关于YOZ平面对称, B点代表相机站位C, DZ轴上的投影。在拍摄时,左右两台相机对准风机的旋转中心, 使风机的旋转中心位于图像正中央。图中,φ1, φ2为绕Y轴方向的旋转角,ω代表相机站位的仰角,由于并未涉及到Z方向的旋转,绕Z方向的旋转角κ=0。由此可以推出C, D站位的旋转矩阵RC, RD分别为:

      $ \left\{ { \begin{array}{l} {\mathit{\boldsymbol{R}}_C} = \left[ {\begin{array}{*{20}{c}} {\cos \kappa }&{ - \sin \kappa }&0\\ {\sin \kappa }&{\cos \kappa }&0\\ 0&0&1 \end{array}} \right] \times \\ {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} \left[ {\begin{array}{*{20}{c}} 1&0&0\\ 0&{\cos \omega }&{\sin \omega }\\ 0&{ - \sin \omega }&{\cos \omega } \end{array}} \right]\left[ {\begin{array}{*{20}{c}} {\cos {\varphi _1}}&0&{ - \sin {\varphi _1}}\\ 0&1&0\\ {\sin {\varphi _1}}&0&{\cos {\varphi _1}} \end{array}} \right] = \\ {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} \left[ {\begin{array}{*{20}{l}} {{R_{C,11}}}&{{R_{C,12}}}&{{R_{C,13}}}\\ {{R_{C,21}}}&{{R_{C,22}}}&{{R_{C,23}}}\\ {{R_{C,31}}}&{{R_{C,32}}}&{{R_{C,33}}} \end{array}} \right]\\ {\mathit{\boldsymbol{R}}_D} = \left[ {\begin{array}{*{20}{c}} {\cos \kappa }&{ - \sin \kappa }&0\\ {\sin \kappa }&{\cos \kappa }&0\\ 0&0&1 \end{array}} \right]\left[ {\begin{array}{*{20}{c}} 1&0&0\\ 0&{\cos \omega }&{\sin \omega }\\ 0&{ - \sin \omega }&{\cos \omega } \end{array}} \right] \times \\ {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} \left[ {\begin{array}{*{20}{c}} {\cos {\varphi _2}}&0&{ - \sin {\varphi _2}}\\ 0&1&0\\ {\sin {\varphi _2}}&0&{\cos {\varphi _2}} \end{array}} \right] = \left[ {\begin{array}{*{20}{c}} {{R_{D,11}}}&{{R_{D,12}}}&{{R_{D,13}}}\\ {{R_{D,21}}}&{{R_{D,2}}}&{{R_{D,23}}}\\ {{R_{D,31}}}&{{R_{D,22}}}&{{R_{D,33}}} \end{array}} \right] \end{array} } \right. $

      (1)

      结合图 1与(1)式,C, D两点相机的成像过程可由下面的共线方程组得到:

      $ \begin{array}{l} \begin{array}{*{20}{c}} {{f_{{x_1}}}:{x_1} = }\\ { - f\frac{{{R_{C,11}}\left[ {X - \left( { - {X_C}} \right)} \right] + {R_{C,12}}\left( {Y - {Y_C}} \right) + {R_{C,13}}\left( {Z - {Z_C}} \right)}}{{{R_{C,31}}\left[ {X - \left( { - {X_C}} \right)} \right] + {R_{C,32}}\left( {Y - {Y_C}} \right) + {R_{C,33}}\left( {Z - {Z_C}} \right)}},}\\ {{f_{{y_1}}}:{y_1} = }\\ { - f\frac{{{R_{C,21}}\left[ {X - \left( { - {X_C}} \right)} \right] + {R_{C,22}}\left( {Y - {Y_C}} \right) + {R_{C,23}}\left( {Z - {Z_C}} \right)}}{{{R_{C,31}}\left[ {X - \left( { - {X_C}} \right)} \right] + {R_{C,32}}\left( {Y - {Y_C}} \right) + {R_{C,33}}\left( {Z - {Z_C}} \right)}},}\\ {{f_{{x_2}}}:{x_2} = }\\ { - f\frac{{{R_{D,11}}\left[ {X - \left( { - {X_C}} \right)} \right] + {R_{D,12}}\left( {Y - {Y_C}} \right) + {R_{D,13}}\left( {Z - {Z_C}} \right)}}{{{R_{D,31}}\left[ {X - \left( { - {X_C}} \right)} \right] + {R_{D,32}}\left( {Y - {Y_C}} \right) + {R_{D,33}}\left( {Z - {Z_C}} \right)}},} \end{array}\\ \begin{array}{*{20}{c}} {{f_{{y_2}}}:{y_2} = }\\ { - f\frac{{{R_{D,21}}[X - ( - {X_C})] + {R_{D,22}}(Y - {Y_C}) + {R_{D,23}}(Z - {Z_C})}}{{{R_{D,31}}[X - ( - {X_C})] + {R_{D,32}}(Y - {Y_C}) + {R_{D,33}}(Z - {Z_C})}}} \end{array} \end{array} $

      (2)

      式中,(x1y1)代表C点相机拍摄的被测点像面坐标,(x2y2)代表D点相机拍摄的被测点像面坐标,f为相机主距,XYZ代表被测点的空间坐标,(-XCYCZC)和(XCYCZC)分别代表C点和D点相机位置的空间坐标,RC, ijRD, ij (i, j=1, 2, 3)代表(1)式中两台相机相机坐标系与空间坐标系相互转换的旋转矩阵中的元素。将共线方程线性化,可以利用光线束平差求解目标点的最优坐标值。对(2)式进行1阶泰勒展开并转化为矩阵形式, 可得到下式:

      $ \begin{align} & \mathit{\boldsymbol{F}}=\left[ \begin{matrix} {{h}_{{{x}_{1}}}} \\ {{h}_{{{y}_{1}}}} \\ {{h}_{{{x}_{2}}}} \\ {{h}_{{{y}_{2}}}} \\ \end{matrix} \right]=\left[ \begin{matrix} {{x}_{1}}-{{f}_{{{x}_{1}}}}({{X}_{0}},{{Y}_{0}},{{Z}_{0}}) \\ {{y}_{1}}-{{f}_{{{y}_{1}}}}({{X}_{0}},{{Y}_{0}},{{Z}_{0}}) \\ {{x}_{2}}-{{f}_{{{x}_{2}}}}({{X}_{0}},{{Y}_{0}},{{Z}_{0}}) \\ {{y}_{2}}-{{f}_{{{y}_{2}}}}({{X}_{0}},{{Y}_{0}},{{Z}_{0}}) \\ \end{matrix} \right]= \\ & {{\left[ \begin{matrix} \frac{\text{d}{{f}_{{{x}_{1}}}}}{~\text{d}X} & \frac{\text{d}{{f}_{{{x}_{1}}}}}{~\text{d}Y} & \frac{\text{d}{{f}_{{{x}_{1}}}}}{~\text{d}Z} \\ \frac{\text{d}{{f}_{{{y}_{1}}}}}{~\text{d}X} & \frac{\text{d}{{f}_{{{y}_{1}}}}}{~\text{d}Y} & \frac{\text{d}{{f}_{{{y}_{1}}}}}{~\text{d}Z} \\ \frac{\text{d}{{f}_{{{x}_{2}}}}}{~\text{d}X} & \frac{\text{d}{{f}_{{{x}_{2}}}}}{~\text{d}Y} & \frac{\text{d}{{f}_{{{x}_{2}}}}}{~\text{d}Z} \\ \frac{\text{d}{{f}_{{{y}_{2}}}}}{~\text{d}X} & \frac{\text{d}{{f}_{{{y}_{2}}}}}{~\text{d}Y} & \frac{\text{d}{{f}_{{{y}_{2}}}}}{~\text{d}Z} \\ \end{matrix} \right]}_{\begin{smallmatrix} X={{X}_{0}} \\ Y={{Y}_{0}} \\ Z={{Z}_{0}} \end{smallmatrix}}}\left[ \begin{matrix} {{\delta }_{{{X}_{0}}}} \\ {{\delta }_{{{Y}_{0}}}} \\ {{\delta }_{{{Z}_{0}}}} \\ \end{matrix} \right]=\mathit{\boldsymbol{J}}\left[ \begin{array}{*{35}{l}} {{\delta }_{{{X}_{0}}}} \\ {{\delta }_{{{Y}_{0}}}} \\ {{\delta }_{{{Z}_{0}}}} \\ \end{array} \right] \\ \end{align} $

      (3)

      式中,hx1hy1hx2hy2为像面坐标的测量值和计算值的差值,F为差值矩阵,(X0Y0Z0)为泰勒展开时被测点空间坐标的初始值,其中dfxk/dX,dfxk/dY,dfxk/dZ,dfyk/dX,dfyk/dY,dfyk/dZ(k=1, 2)为(2)式中各式对空间坐标的偏导数,J为偏导数矩阵,δX0δY0δZ0为空间坐标的修正值,在迭代过程中不断修正初始值。将上式化简变形可得到空间坐标修正值:

      $ \left[ \begin{matrix} {{\delta }_{{{X}_{0}}}} \\ {{\delta }_{{{Y}_{0}}}} \\ {{\delta }_{{{Z}_{0}}}} \\ \end{matrix} \right]={{({{{\mathit{\boldsymbol{J}}}}^{\text{T}}}{\mathit{\boldsymbol{J}}})}^{-1}}{{{\mathit{\boldsymbol{J}}}}^{\text{T}}}{\mathit{\boldsymbol{F}}} $

      (4)

      将(4)式代入最小二乘迭代,不断修正空间坐标。随着迭代次数增加,坐标值不断接近真实值。

    • 空间坐标测量误差的标准差可以由(4)式通过协方差传播率推导得到。

      $ {{{\mathit{\boldsymbol{D}}}}_{\delta }}={{\left( {{{\mathit{\boldsymbol{J}}}}^{\text{T}}}{\mathit{\boldsymbol{J}}} \right)}^{-1}}{{{\mathit{\boldsymbol{J}}}}^{\text{T}}}{{{\mathit{\boldsymbol{D}}}}_{{\mathit{\boldsymbol{F}}}}}{{\left[ {{\left( {{{\mathit{\boldsymbol{J}}}}^{\text{T}}}{\mathit{\boldsymbol{J}}} \right)}^{-1}}{{{\mathit{\boldsymbol{J}}}}^{\text{T}}} \right]}^{\text{T}}} $

      (5)

      式中,Dδ是空间坐标修正值的协方差矩阵, DF是像面坐标的协方差矩阵。在理想状态下,各像面点相互独立,所有像面点坐标误差水平相同,则像面点的协方差矩阵DF为像面点定位误差标准差的平方σF2与单位矩阵E的乘积,如下式所示:

      $ {{{\mathit{\boldsymbol{D}}}}_{{\mathit{\boldsymbol{F}}}}}={{\sigma }_{{\mathit{\boldsymbol{F}}}}}^{2}{\mathit{\boldsymbol{E}}} $

      (6)

      将(6)式代入(5)式,得到空间坐标修正值的协方差矩阵。

      $ {{{\mathit{\boldsymbol{D}}}}_{\delta }}={{\sigma }_{{\mathit{\boldsymbol{F}}}}}^{2}{{({{{\mathit{\boldsymbol{J}}}}^{\text{T}}}{\mathit{\boldsymbol{J}}})}^{-1}} $

      (7)

      Dδ的对角线元素分别为σX2σY2σZ2,其中σXσYσZ为被测点在XYZ 3个方向上误差的标准差。将像面定点误差标准差代入(5)式,求不同被测点的Dδ矩阵,并对Dδ的对角线元素的平方根求均值,记作在当前相机站位下空间坐标测量误差的标准差,如下式所示:

      $ \sigma =\sum\limits_{i=1}^{n}{\frac{\sqrt{{{{\mathit{\boldsymbol{D}}}}_{\delta ,i}}(1,1)}+\sqrt{{{{\mathit{\boldsymbol{D}}}}_{\delta ,i}}(2,2)}+\sqrt{{{{\mathit{\boldsymbol{D}}}}_{\delta ,i}}(3,3)}}{3n}} $

      (8)

      式中,n为被测点数。

    • 本文中提出一种变异改进型遗传算法作为网络优化方法。在遗传算法中,种群经过多次迭代优化后,可以得到适应度较高的次优站位,该站位与最优站位在空间位置上十分接近。在对该站位的变异操作中,传统的变异方法通过随机替换站位中的参量来实现变异,这种方法随机性太强,反而有可能劣化当前的次优站位。本文中通过改进变异方法,检索当前站位相邻区域,提高了局部搜索能力,可以更快地检索到与次优站位接近的最优站位,进而提高收敛速度。具体流程如图 2所示。

      Figure 2.  Flow chart of network optimization algorithm

      (1) 初始化种群:根据实际情况建立约束条件,选择多组满足约束条件的相机站位,采用实值编码建立初始站位集合(初始种群)。

      (2) 计算适应度:根据(8)式计算站位集合中各站位的空间坐标测量误差的标准差,以标准差的倒数作为适应度,如下式所示:

      $ S=\frac{1}{\sigma } $

      (9)

      式中, S为适应度,σ为标准差。

      (4) 选择复制:根据站位集合中各站位的适应度选择站位进行复制,适应度高的站位被选择的概率高,被选择的站位继续进行交叉变异操作。

      (5) 交叉:随机将种群中序号接近的站位XZ坐标进行交换,例如将P1P2两站位交叉得到新的站位P3P4

      $ \begin{array}{*{35}{l}} {{P}_{1}}({{X}_{1}},{{Z}_{1}}) \\ {{P}_{2}}({{X}_{2}},{{Z}_{2}}) \\ \end{array}\Rightarrow \begin{array}{*{35}{l}} {{P}_{3}}({{X}_{1}},{{Z}_{2}}) \\ {{P}_{4}}({{X}_{2}},{{Z}_{1}}) \\ \end{array} $

      (10)

      (6) 变异:这里的变异指的是将站位坐标在相邻区域内进行随机位移。这种变异方法局部搜索能力更强,收敛更快,如下式所示:

      $ {{P}_{\text{ old }}}(X,Z)\Rightarrow {{P}_{\text{ new }}}(X\pm a,Z)\text{or}{{P}_{\text{ new }}}(X,Z\pm a) $

      (11)

      式中, a为随机位移量, PnewPold是新、旧站位。

      (7) 建立新种群:将经过交叉变异之后的站位重新计算适应度,根据适应度选取站位生成新的站位集合,继续迭代。

    • 图 1所示,在仿真实验中,以轮毂高度80m,叶片长度40m,风轮直径80m的1.5MW风机作为被测物,将两台相机对称地布设在风机正前方的水平地面上,在风机叶片上均匀地布设被测点。

    • 仿真实验过程中,应对站位集合中的站位设立约束条件,使求得的站位具有合理性和可行性。

    • 在实验过程中为了使风机叶片上所有的被测点都能被相机采集到,必须对相机镜头光心到风轮中心的距离做出限制。

      $ L\ge \frac{fD}{d} $

      (12)

      式中,L是相机光心到风机旋转中心的距离,D是风轮直径,d是相机拍摄图像的最小尺寸,f是相机主距。例如在仿真实验时,大型风机的叶片长度为40m,风轮直径为80m,仿真相机拍摄的图片尺寸为4872pixel×3248pixel,像素尺寸为7.4μm,主距为20mm,那么L应大于67m。

    • 在摄影测量中,回光反射目标点具有可视角度范围,为了得到被测点明亮清晰的图像,应对被测点所在平面与相机光轴的夹角做出相应的约束。该约束在风电叶片测量中体现为相机光轴与叶片旋转平面的夹角大于30°。

      图 3所示,角αβ代表相机与风轮上极限位置处被测点T(XT, YT, ZT)的连线与风轮所在平面的夹角,可视角度约束条件可通过下式体现:

      Figure 3.  The angle between the optical axis of the camera and the plane of the wind wheel

      $ \left\{ \begin{array}{*{35}{l}} \tan \alpha =\frac{CE}{TE}\ge \tan {{30}{{}^\circ }} \\ \tan \beta =\frac{DF}{TF}\ge \tan {{30}{{}^\circ }} \\ \end{array} \right. $

      (13)
    • 在仿真实验中采用空间坐标测量误差的标准差作为优化目标,需要确定像面点定位误差的标准差。摄影测量目标点定位精度一般在1/20pixel。这里采用的相机像素尺寸为7.4μm,则像面点定位误差的标准差为0.37μm,再结合(7)式和(8)式就可以得到相机位于该站位时空间坐标测量误差的标准差。

    • 在网络优化的仿真实验中,选取28对均匀散布的对称站位作为初始站位集合,如图 4所示。

      Figure 4.  Initial set of stations

      图中每一对关于Y-O-Z平面对称的三角形代表站位集合中的一组站位,每组2个站位Z坐标相等,X坐标互为相反数,Y坐标代表高度,由于相机位于地面,默认Y坐标为0。阴影区域为不满足第2.1节中所提到约束条件的站位区域。将初始站位集合代入第1.3节中所述的网络优化算法,迭代一定次数得出最优站位。最优站位的空间坐标测量误差标准差和所有站位标准差均值的变化曲线如图 5所示。

      Figure 5.  Variation of the standard deviation of the optimal station and the mean standard deviation of all stations

      图 5可知,随着迭代次数增加,最优站位的空间坐标标准差不断减少并最终稳定于2.7mm。种群中最优站位的XZ坐标随迭代次数增加的变化如图 6图 7所示。

      Figure 6.  Variation of the X-coordinate of the optimal station with iterations

      Figure 7.  Variation of the Z-coordinate of the optimal station with iterations

      在前30次迭代时,当前种群的最优站位的XZ坐标变化明显,迭代次数超过70次后,相机站位的X坐标收敛于73.5m,Z坐标收敛于140.8m。仿真实验中网络优化的结果如表 1所示。

      Table 1.  The standard deviation and coordinates of the optimal station in iteration

      iterations the spatial standard deviation of the optimal station/mm X-coordinate/ mm Z-coordinate/ mm
      1 4.100 8×104 2.1×105
      10 3.027 6.353×104 1.467×105
      20 2.884 6.536×104 1.418×105
      30 2.749 7.258×104 1.406×105
      40 2.741 7.551×104 1.421×105
      50 2.741 7.551×104 1.421×105
      60 2.741 7.370×104 1.409×105
      70 2.741 7.352×104 1.408×105
      80 2.740 7.352×104 1.408×105
      90 2.740 7.352×104 1.408×105
      100 2.740 7.352×104 1.408×105

      表 1中数据可知,在仿真实验中,最优的一组相机站位坐标为(-7.352×104, 0, 1.408×105)与(7.352×104, 0,1.408×105)(单位为mm),相机位于该组站位时,空间坐标测量误差的标准差为2.7mm。

    • 在实测实验中,采用叶片长度为3.5m的风机模型为被测物,在上面均匀布设被测点,如图 8所示。

      Figure 8.  Experimental site and the turbine to be measured

      使用V-STARS摄影测量系统测量被测点的空间坐标作为后续相机站位优化实验中被测点的坐标真值。实测实验中采用两台联合视觉技术(allied vision technologies,AVT)高精度工业相机拍摄图像,其分辨率为4872pixel×3248pixel,像素尺寸为7.4μm,镜头焦距为20mm。

      将两台相机布设完毕,触发相机同时拍摄,使用动态摄影测量软件对拍摄的图片进行处理,得到每个目标点的空间坐标值。

    • 先将仿真得到的最优站位等比缩放,将相机布设在最优站位进行拍摄测量,再选择多个随机站位拍摄测量,计算相邻目标点间的距离,与V-STARS系统拍摄得到的真值对比,得到的数据如表 2所示。

      Table 2.  Comparison of measurement results between optimal and stochastic stations

      station number maximum error/mm average error/mm root mean square/mm relative error/%
      1(optimal station) 1.392 0.617 0.836 0.009
      2 2.241 1.273 1.562 0.018
      3 3.277 2.852 2.872 0.041
      4 5.639 3.057 3.568 0.051
      5 2.715 2.556 2.564 0.037
      6 3.415 2.909 2.931 0.042
      7 2.685 2.036 2.097 0.029
      8 3.415 2.523 2.640 0.036
      9 4.382 3.142 3.257 0.045
      10 5.149 3.127 3.496 0.045

      表 2中,站位1为仿真得到的最优站位,其余站位为随机站位。由表中数据可知,当相机位于站位1时,测量误差小于所有的随机站位,证明当相机位于网络优化后的最优站位时,其摄影测量的精度最高。

    • 针对大型风电叶片动态摄影测量的网络优化问题,本文中以摄影测量3维重构时产生的空间坐标测量误差标准差为优化目标,根据摄影测量中的约束条件对站位做出限制,设计变异操作改进型遗传算法为基础的网络优化算法,通过该算法在仿真实验中得到了相机的最优站位。在实测实验中,当相机位于最优站位时,其测量精度优于其它站位,为风电叶片摄影测量的网络优化提供了理论技术支持。

参考文献 (14)

目录

    /

    返回文章
    返回