Page 121 - 2024年第55卷第10期
P. 121

(4)补偿技术层。
                  考虑到传统的全球导航卫星系统定位技术无法满足定位精度或精度保证率的要求,本研究在数学
              模型中定义了一套补偿技术来解决这一问题。补偿技术层定义了在悬崖边壁、倒悬体遮挡等多种复杂
              环境下压实施工质量监控的补偿内容。压实参数监控精度( PA)和监控精度保证率(GA)在约束条件
                                                                       a                      a
              层中描述如式(5)所示:
                                                     PA= PA
                                                    {   a    GNSS PA c                                 (5)
                                                     GA= GA
                                                        a    GNSS GA c
              式中:PA       和 PA分别是全球导航卫星系统技术和补偿技术的监测精度;GA                                和 GA分别是全球导
                       GNSS    c                                                     GNSS    c
              航卫星系统技术和补偿技术的监测精度保证率;表示改进的卡尔曼滤波算法。补偿技术层是对以往
              研究的区别和改进。
                  ( 5)控制参数计算公式层。
                  根据定位技术获得的监测数据,定义了控制参数的计算公式,如式(6)所示:
                                          D = Z (m,n,k) - Z (m,n,k - 1)
                                            a
                                                                       2
                                                       2
                                                                2
                                          S= (x- x)+ (y- y)(z - z) ?(t - t)
                                           a  槡   1  2     1  2    1  2    2  1
                                          R = N(P(x,y,z,t) =(m,n,k))                                    (6)
                                            a
                                          PA= (x,y,z) - (x,y,z)
                                             a     m   m   m    t   t  t
                                          GA= N(PA≥PA)?N(ALL) × 100 %
                                                          s
                                                    a
                                             a
              式中:Z(m,n,k)为点(m,n,k)在碾压单元坝面上的高度;x、y、z为碾压机位置的空间坐标;t
              为记录时间;N()是一个计数函数。详细的计算原理和公式参考了文献[30];P()是一种算法,可以
              获得单位曲面点,然后计算压实遍数;( x、y,z)(x,y,z)为监控空间坐标和实际空间坐标。
                                                    m
                                                        m
                                                                    t
                                                                        t
                                                            m
                                                                 t
                  该数学模型考虑了在悬崖边壁、倒悬体遮挡等多种复杂环境下压实施工质量监控的环境特征,如
              在数学模型中定义的环境属性层,弥补了已有研究中模型缺少环境特征的不足                                     [31 - 34] 。特别地,在该模
              型中,首次提出了耦合 GNSS、UWB及 RTS的智能监控 PCT,以解决因全球导航卫星系统技术在悬崖
              边壁、倒悬体遮挡等多种复杂环境下信号不足,大坝施工质量监控技术的性能精度降低的问题,从而
              更好地适应高山峡谷地区大坝施工碾压监控高精度的需求。
              3 融合 GNSS?UWB?RTS碾压混凝土坝碾压施工智能监控 PCT

              3.1 卡尔曼滤波算法 卡尔曼滤波通过输入的观测值与系统估计值对系统的状态进行最优估计,是一
              种线性最小方差估计的方法,具有能够实时输入数据并进行预测、计算量较小以及便于计算机编程等
              诸多优点。基于碾压机运动模型的卡尔曼滤波过程如下。
                  碾压机能够根据实时接收的卫星信号计算出实时的三维坐标 x、y、z,其中 t为时间,时间间隔
                                                                                 t
                                                                          t
                                                                              t
              为 Δ t。同时,碾压机的三维速度为 v、v、v,碾压机的三维加速度为 a、a、a,由于碾压机在作
                                                       zt
                                                                                     yt
                                               xt
                                                   yt
                                                                                         zt
                                                                                 xt
              业过程中近似为匀速直线运动,且时间间隔 Δ t较短,因此忽略加速度对碾压机位置的影响。因此定
              义状态矩阵为:
                                               X = [ x y z v      v   v]  T                             (7)
                                                t    t  t   t  xt  yt  zt
                  系统在 Δ t的时间内状态转移过程为:
                                                              I       0
                                                               3 × 3
                                              X = AX_ ,A =                                              (8)
                                                      t - 1
                                                t
                                                               0   Δ t·I
                                                                        3 × 3
                  系统的先验估计协方差 P 为:
                                         t
                                                                T
                                                     P= AP_ A + Q                                       (9)
                                                       t
                                                             t - 1
                  系统的先验估计给出基于前一时刻的系统状态,以及碾压机的运动模型得出的系统当前时刻状态的
                                                                                                   2
                                                                                              —   1 5 1 —
   116   117   118   119   120   121   122   123   124   125   126