一种基于虚拟加速度的gnss变形监测方法及系统的制作方法
【技术领域】
[OOOU 本发明设及变形监测领域，特别设及一种基于虚拟加速度的GNSS变形监测方法 及系统。
【背景技术】
[0002] 变形监测主要是利用各种仪器和手段对被监测对象的变形进行持续观测，获得被 监测对象的动态形变数据，通过对动态形变数据的统计分析，建立被监测对象变化情况的 数学预测模型，对被监测对象未来变化趋势做出合理预测，当变形超过某个特定允许值，贝U 认为有可能是灾害事故发生的前兆。由于变形监测是对被监测对象空间位置及内部形态随 时间变化特征进行连续监测，故对监测信息的实时性、连续性和准确性要求较高。
[0003] GNSS(全球卫星导航定位系统）具有速度快、全天候观测和测点间无需通视等优 点，其数据采集、传输、处理、分析、显示和存储等自动化程度高，已成为高精度同步变形监 测的主要技术手段之一。GNSS主要包括PPP技术（精密单点定位技术）和RTK技术（实时 动态测量定位技术）。PPP技术使用单台测量设备即可完成定位的精密解算，但在实践中必 须采用延时较长的精密星历，难实现较高实时性要求。而RTK技术可直接使用广播星历，即 可实现较高实时性要求。
[0004] RTK技术是基于高精度载波相位的差分技术，其在监测区布设多个接收站，并通过 网络将数据传至控制中屯、，选择稳定区接收站作为基准站，形变区接收站作为流动站，通过 组基线的方式解算流动站的位置，从而推算出被监测对象的变形情况。实践表明，RTK技术 中基于基线解算流动站的位置可很好的消除卫星轨道误差、电离层误差和对流层误差等站 间公共误差，但接收机观测噪声误差不易消除。
【发明内容】
[0005] 本发明所要解决的技术问题是提供一种基于虚拟加速度的GNSS变形监测方法及 系统，其能有效削弱RTK技术中实时变形监测中接收机观测噪声误差，及事后变形监测中 接收机观测噪声误差，适用于大型建筑物的长期微小自然形变的变形监测。
[0006] 本发明解决上述技术问题的技术方案如下：
[0007] 一种基于虚拟加速度的GNSS变形监测方法，包括如下步骤：
[000引步骤1，虚拟流动站准静态条件的加速度，对加速度W第一采样间隔进行采样并存 储。
[0009] 步骤2,提取对流动站W第二采样间隔开始采样时流动站的速度和位置，并将所述 速度和位置分别作为初始速度和初始位置；设置预设零偏。
[0010] 步骤3,提取初始速度、初始位置、预设零偏及所述初始速度和初始位置相应的加 速度，并对所述加速度进行零偏修正，计算单位第一采样间隔后流动站的末速度和末位置。
[0011] 步骤4,将末速度和末位置分别作为初始速度和初始位置，重复执行步骤3,至第 一采样间隔的累积时长与第二采样间隔相等时，提取此刻计算所得的末速度和末位置分别 作为基于加速度解算所得流动站的第一速度和第一位置。
[0012] 步骤5,利用所述第一速度、第一位置和预设零偏，结合基于基线解算所得流动站 的第二速度和第二位置，并将所述第二速度修正到预定取值范围内，建立卡尔曼滤波模型 求解流动站的运动状态改变量，进一步求解实时修正所得流动站的第=速度、第=位置和 修正零偏。
[0013] 步骤6,将所述第=速度、第=位置和修正零偏分别作为初始速度、初始位置和预 设零偏，重复执行步骤3至5,求解后续实时修正所得流动站的第=速度、第=位置和修正 零偏，至全部实时修正完成；并W全部第=位置为基础采用固定区间最优平滑算法，求解事 后修正所得流动站的第四位置。
[0014] 本发明的有益效果是：虚拟流动站准静态条件的加速度，基于加速度解算流动站 的第一速度和第一位置，利用第一速度、第一位置和预设零偏，结合RTK技术中基于基线解 算流动站的第二速度和第二位置，并将所述第二速度修正到预定取值范围内，建立卡尔曼 滤波模型求解流动站的运动状态改变量，进一步求解实时修正流动站的第=速度、第=位 置和修正零偏；当全部实时修正完成，W全部第=位置为基础采用固定区间最优平滑算法， 求解事后修正所得流动站的第四位置；其能有效削弱RTK技术中实时变形监测中接收机观 测噪声误差，及事后变形监测中接收机观测噪声误差，适用于大型建筑物的长期微小自然 形变的变形监测。
[0015] 在上述技术方案的基础上，本发明还可W做如下改进。
[0016] 进一步，所述加速度的表达式为：
[0017] a=n"+n,t
[0018] 其中，所述a为加速度，rv为有色噪声，n，t为白噪声。
[0019] 所述对加速度W第一采样间隔进行采样并存储的【具体实施方式】为；W第一采样间 隔采样流动站在ECEF坐标系下X轴、y轴和Z轴=个方向的加速度，并将同一时刻X轴、y 轴和Z轴=个方向的加速度列为一行，按采样先后顺序逐行W文本形式存储。
[0020] 进一步，所述步骤2中预设零偏为0。
[0021] 进一步，所述步骤3的【具体实施方式】包括如下步骤：
[0022] 步骤Al,提取所述初始速度和初始位置对应时刻的第一加速度，及单位第一采样 间隔后的第二加速度，对所述第一加速度根据如下第一公式进行零偏修正得修正第一加速 度，对所述第二加速度根据如下第二公式进行零偏修正得修正第二加速度。
[0023] 所述第一公式如下：
[0024]
[0025] 所述第二公式如下：
[0026] 媒1=(2,+1-护
[0027]其中，所述< 为修正第一加速度，a,为第一加速度，<1为修正第二加速度，aw为 第二加速度，ba为预设零偏。
[002引步骤A2,利用所述修正第一加速度和修正第二加速度，根据如下第S公式和第四 公式分别计算单位第一采样间隔后流动站的末速度和末位置。
[0029]所述第S公式如下：
[0033] 其中，所述At为第一采样间隔，Vt为初始速度，VW为末速度，St为初始位置，SW 为末位置。
[0034] 进一步，所述步骤5的【具体实施方式】包括如下步骤：
[0035] 步骤BI,将基准站位置固定，利用基准站在EWF坐标系下的坐标，根据基线解算 出的向量推算出基于基线解算所得流动站的第二速度和第二位置。
[0036] 步骤B2,利用所述第一速度、第一位置、第二速度和第二位置建立流动站在ECEF坐标系下的系统模型方程和观测模型方程，并将所述第二速度置零，建立卡尔曼滤波模型 求解流动站的运动状态改变量；其中所述运动状态改变量包括流动站的位置改变量、速度 改变量和零偏改变量。
[0037] 所述系统模型方程为：
[00%]
[0039] 其中，所述AP为位置改变量，A/，为位置改变量导数，AV为速度改变量，Al)为速 度改变量导数，Ab为零偏改变量，Ai为零偏改变量导数，n为系统模型白噪声，q为系统模 型驱动白噪声，-!为连续马尔科夫过程状态系数，I为3阶单位矩阵，心为主对角线上 元素全部为的3阶对角矩阵。
[0040] 所述观测模型方程为：
[0041]
[0042] 其中，所述pa为第一位置，PS为第二位置，Va为第一速度，VS为第二速度且为零，nP 为位置观测噪声，riy为速度观测噪声。
[0043] 步骤B3,利用所述流动站的运动状态改变量、第一速度、第一位置和预设零偏，根 据如下第五公式求解实时修正所得流动站的第=速度、第=位置和修正零偏。
[0044] 所述第五公式如下：
[0045]
[0046] 其中，所述P为第S位置，V为第S速度，b为修正零偏。
[0047] 本发明的另一技术方案如下：
[0048] 一种基于虚拟加速度的GNSS变形监测系统，包括加速度生成模块、初始状态提取 模块、基于加速度中间状态生成模块、基于加速度解算模块、实时修正模块和事后修正模 块。
[0049] 所述加速度生成模块，其用于虚拟流动站准静态条件的加速度，对加速度W第一 采样间隔进行采样并存储。
[0050] 所述初始状态提取模块，其用于提取对流动站W第二采样间隔开始采样时流动站 的速度和位置，并将所述速度和位置分别作为初始速度和初始位置；设置预设零偏。
[0051] 所述基于加速度中间状态生成模块，其用于提取初始速度、初始位置、预设零偏及 所述初始速度和初始位置相应的加速度，并对所述加速度进行零偏修正，计算单位第一采 样间隔后流动站的末速度和末位置。
[0052] 所述基于加速度解算模块，其用于将末速度和末位置分别作为初始速度和初始位 置，重复驱动基于加速度中间状态生成模块工作，至第一采样间隔的累积时长与第二采样 间隔相等时，提取此刻计算所得的末速度和末位置分别作为基于加速度解算所得流动站的 第一速度和第一位置。
[0化3] 所述实时修正模块，其用于利用所述第一速度、第一位置和预设零偏，结合基于基 线解算所得流动站的第二速度和第二位置，并将所述第二速度修正到预定取值范围内，建 立卡尔曼滤波模