一种基于单频GNSS定位技术监测高铁大桥动态挠度的方法技术

本发明专利技术提出了一种基于单频GNSS定位技术监测高铁大桥动态挠度的方法，步骤如下：S1、接收机实时获取基站和桥梁监测点的GNSS单频载波相位观测值；S2、利用E‑ADOP和GDOP的因子双重约束方法选择待测卫星；S3、基于球冠谐分析建立对流层实时模型，并解算对流层延迟；S4、建立多路径实时模型，解算多路径误差改正数，并根据误差改正数更新卫地距离数；S5、基于对流层延迟和更新后的卫地距离计算桥梁监测点的GNSS单频载波相位双差观测量，并解算模糊度；S6、将模糊度回代实时计算桥梁监测点的三维坐标值，获得桥梁动态挠度值。本发明专利技术方法操作和计算较为简单，可以自动实时监测高铁大桥的动态挠度，通过建模不断改正测量数据，提高桥梁动态监测的精准度。

【技术实现步骤摘要】
一种基于单频GNSS定位技术监测高铁大桥动态挠度的方法
本专利技术涉及一种实时监测高铁大桥动态挠度的方法，属于桥梁健康监测

技术介绍
随着经济的发展，铁路桥梁越来越普及，在运营过程中铁路桥梁不可避免的会出现老化和损伤，桥梁的承载能力和安全性降低，因此从20世纪90年起，国内外开始了桥梁健康监测方面的研究。在铁路桥梁中，高速铁路由于其运行时速较高，在250到350公里内，且高铁桥梁绝大多数建在野外，跨度大都在100米以上，所以高铁桥梁结构的动力性能需要更高的刚度要求，要有较强的抗弯、抗扭刚度，以防止桥梁出现较大挠度变形和振幅，高铁桥梁的监测系统及内容与一般桥梁的存在一定差异，高铁桥梁的安全性监测更加至关重要。主梁挠度是进行桥梁刚度监测、安全性、适用性评价的最重要的指标之一，对评价高铁桥梁安全性和适用性有重要意义；其中动挠度则直接反映移动活载作用下梁体的挠度动态变化，通过在桥梁荷载试验中获取的动挠度时程曲线,可以计算出反映桥梁动力响应的冲击系数，了解桥梁周期性变位规律，判断桥梁结构的整体性及薄弱部位。目前,国内外在桥梁的动挠度监测方面做了大量研究,动态挠度监测主要方法有以下几种：机器人监测法、激光投射法、连通管法、GPS法等。机器人监测法是由马达驱动全站仪，按照设定的程序依次自动扫描桥梁测点，得到桥梁的变形情况，此方法能够较快速的进行大量程和高精度监测，但成本较高，受环境影响大，不能进行多点同步监测；激光投射式挠度监测法是利用激光的准直性，在桥梁的测点上安装发射源或目标靶，通过计算激光投影的坐标来推算桥梁挠度的方法，此方法的成本较高，且容易受尘埃、水汽等物体的干扰导致数据失真；连通管挠度监测法是利用连通管内液面水平的原理，沿梁体布设连通管道，通过测点处连通管内液位的变化，换算得到桥梁挠度的变化，此方法原理简单、成本较低、精度高，但是只能人工测量，实时性太差，采集的数据无法及时处理；GPS监测法利用卫星与设置在桥梁上监测点的接收机来测量解算桥梁动态挠度，不受通视条件、距离影响，可全天候自动实时获取测点三维绝对坐标，但之前一直存在接收机设备动态监测精度低、成本高的问题。随着近些年高采样频率GNSS接收机的出现、实时差分法等数据处理方法的完善、BDS/伽利略/GPS多系统技术迅速发展以及接收机设备价格下降，GPS监测法越来越多的运用在桥梁监测中，现阶段GNSS技术监测桥梁更多应用于公路桥梁，因为高铁大桥的刚度要求，变形小，相应的研究应用还比较少。
技术实现思路
针对目前高铁大桥健康监测技术不够成熟、成本较高的问题，本专利技术提出了一种基于单频GNSS定位技术监测高铁大桥动态挠度的方法，在高铁大桥的监测点设置接收机，结合卫星实时获取基站和桥梁监测点的GNSS单频载波相位观测值，通过对流层延迟和多路径误差改正数消除观测值的一部分误差，最终利用改正后的观测值实时解算高铁大桥监测点的三维坐标。为解决上述技术问题，本专利技术采用了如下技术手段：一种基于单频GNSS定位技术监测高铁大桥动态挠度的方法，具体包括以下步骤：S1、利用设置在高铁大桥上的接收机实时获取基站和桥梁监测点的GNSS单频载波相位观测值；S2、利用E-ADOP和GDOP的因子双重约束方法选择待测卫星；S3、基于球冠谐分析建立对流层实时模型，并解算对流层延迟；S4、建立多路径实时模型，解算多路径误差改正数，并根据误差改正数更新卫地距离数；S5、基于对流层延迟和更新后的卫地距离计算桥梁监测点的GNSS单频载波相位双差观测量，并解算模糊度；S6、将模糊度回代实时计算桥梁监测点的三维坐标值，获得桥梁动态挠度值。进一步的，步骤S2的具体操作如下：S21、计算GDOP值，公式如下：其中，tr(.)表示对矩阵主对角线元素求和，A为由某种特征测量值建立的线性方程组的系数矩阵，Y＝AX，Y表示已知基站位置的向量，X为监测点位置的未知向量。当ATA为非奇异矩阵，即上述方程数大于未知数个数，则采用最小二乘法计算X，公式如下：X＝(ATA)-1ATY(2)S22、计算E-ADOP值，公式如下：其中，m代表卫星的个数，ωs＝sin(θs)2，θs为卫星s的仰角，s＝1,2,…,m，Cr为伪距代数余子式矩阵，Ch为代数余子式矩阵的阶。S23、GDOP值与卫星星座构成的多面体有关，GDOP对定位误差有放大作用，所以GDOP值越小，定位精度越好；根据E-ADOP数值可以判断卫星数和所选卫星中高高度角卫星的比例，卫星数增多、高高度角卫星比例增大会提高定位效果，所以E-ADOP值越小，定位成功率越高，精度越高。计算不同卫星组合对应的GDOP和E-ADOP值，选择使GDOP和E-ADOP值最小的卫星作为待测卫星。进一步的，步骤S3中对流层延迟的计算公式如下：其中，r为地心距离，δ为球冠余纬坐标，λ为球冠经度坐标，a为地球平均半径，a＝6378.137km，H为建模参考面高度,nk为非整阶勒让德函数序列，k是根nk的序号，k和l为自然数，且k≥l，和为规格化球冠谐系数，为非整阶的Schmidt正交规格化缔合勒让德函数。用n＝nk(m)表示n的实根，即n＝nk(m)，有：其中，F为超几何函数，为规格化因子，的表达式如下：其中，Γ表示伽马函数。进一步的，步骤S4的具体操作如下：S41、设接收机天线为R，接受机接受到的信号波形为在x＝x0处对信号产生多路径效应，有：其中，θ为卫星高度角,β为波形在x0处的切线与水平方向的夹角,h为天线相位中心到反射面的垂直距离，A为波形振幅，ω为角速度，为x0处的相位。S42、令根据当前历元时刻的波形和卫星高度角计算x0：S43、计算x0处产生的多路径信号的距离，公式如下：其中，S44、计算x0处波形到直达信号的垂直距离，公式如下：S45、计算卫星发射的直达信号的距离，公式如下：S46、计算多路径效应误差改正数，公式如下：S47、根据多路径误差改正数更新卫地距离数ρ：其中，为卫地距离的实测值。进一步的，步骤S5中的GNSS单频载波相位双差观测量的计算公式如下：其中，为t1时刻接收机i、j和卫星p、q对应的双差观测值，为t1时刻接收机i、j和卫星p、q对应的卫地距离二次差值，f为信号频率，c为光速，为t1时刻接收机i、j和卫星p、q对应的整周模糊度二次差值，为t1时刻接收机i、j和卫星p、q对应的电离层延迟改正数二次差值，为t1时刻接收机i、j和卫星p、q对应的对流层延迟改正数二次差值。采用以上技术手段后可以获得以下优势：本专利技术提出了一种基于单频GNSS定位技术监测高铁大桥动态挠度的方法，在高铁大桥上的监测点安装接收机，自动接收卫星发送的信号，采用双重因子约束法筛选所有卫星，提高监测信号的精度，针对所选卫星的观测值，分别分析对流层延迟和多路径效应，不断改正观测数据，根据基站点坐标和改正后的观测数据计算高铁大桥上监测点的三维坐标，监测点的高程坐标值变化即大桥的动态挠度。与传统监测方法相比，本专利技术方法无需太多人工参与，操作较为简单，抗干扰能力强，不易受周围环境影像，整个监测过程的计算量不大，监测精度尤其是动态监测的精度比较高，可以实时、高效监测高铁大桥的健康情况。附图说明图1为本专利技术一种基于单频GNSS定位技术监测高铁大桥动态挠度的方法的步骤流程图。图2为本专利技术一种本文档来自技高网...

【技术保护点】
1.一种基于单频GNSS定位技术监测高铁大桥动态挠度的方法，其特征在于，包括以下步骤：S1、利用设置在高铁大桥上的接收机实时获取基站和桥梁监测点的GNSS单频载波相位观测值；S2、利用E‑ADOP和GDOP的因子双重约束方法选择待测卫星；S3、基于球冠谐分析建立对流层实时模型，并解算对流层延迟；S4、建立多路径实时模型，解算多路径误差改正数，并根据误差改正数更新卫地距离数；S5、基于对流层延迟和更新后的卫地距离计算桥梁监测点的GNSS单频载波相位双差观测量，并解算模糊度；S6、将模糊度回代实时计算桥梁监测点的三维坐标值，获得桥梁动态挠度值。

【技术特征摘要】
1.一种基于单频GNSS定位技术监测高铁大桥动态挠度的方法，其特征在于，包括以下步骤：S1、利用设置在高铁大桥上的接收机实时获取基站和桥梁监测点的GNSS单频载波相位观测值；S2、利用E-ADOP和GDOP的因子双重约束方法选择待测卫星；S3、基于球冠谐分析建立对流层实时模型，并解算对流层延迟；S4、建立多路径实时模型，解算多路径误差改正数，并根据误差改正数更新卫地距离数；S5、基于对流层延迟和更新后的卫地距离计算桥梁监测点的GNSS单频载波相位双差观测量，并解算模糊度；S6、将模糊度回代实时计算桥梁监测点的三维坐标值，获得桥梁动态挠度值。2.根据权利要求1所述的一种基于单频GNSS定位技术监测高铁大桥动态挠度的方法，其特征在于，步骤S2的具体操作如下：S21、计算GDOP值，公式如下：其中，tr(.)表示对矩阵主对角线元素求和，A为由某种特征测量值建立的线性方程组的系数矩阵，Y＝AX，Y表示已知基站位置的向量，X为监测点位置的未知向量；当ATA为非奇异矩阵，即上述方程数大于未知数个数，则采用最小二乘法计算X，公式如下：X＝(ATA)-1ATYS22、计算E-ADOP值，公式如下：其中，m代表卫星的个数，ωs＝sin(θs)2，θs为卫星s的仰角，s＝1,2,…,m，Cr为伪距代数余子式矩阵，Ch为代数余子式矩阵的阶；S23、计算不同卫星组合对应的GDOP和E-ADOP值，选择使GDOP和E-ADOP值最小的卫星作为待测卫星。3.根据权利要求1所述的一种基于单频GNSS定位技术监测高铁大桥动态挠度的方法，其特征在于，步骤S3中对流层延迟的计算公式如下：其中，r为地心距离，δ为球冠余纬坐标，λ为球冠经度坐标，a为地球平均半径，a＝637...

【专利技术属性】
技术研发人员：张秋昭，段文再，吴来义，戴新军，岳青，刘华，
申请(专利权)人：中国矿业大学，中铁大桥南京桥隧诊治有限公司，
类型：发明
国别省市：江苏,32