[发明专利]一种基于GNSS的地表变形连续监测方法

摘要

本发明公开了一种基于GNSS的地表变形连续监测方法，包括以下步骤：步骤一、建立GNSS观测站；步骤二、获取观测基线的长度变形量；步骤三、计算观测基线的应变分量；步骤四、计算观测基线的应变分量的平均值；步骤五、计算待监测区域的应变参数；步骤六、计算观测基线的剪切斜应变参数；步骤七、重复步骤二至步骤六，对待监测区域进行连续监测。本发明结构简单、设计合理，通过检测GNSS观测站的位移量，得到观测基线的长度变形量，由观测基线的长度变形量计算待监测区域的应变参数和每一条观测基线的剪切斜应变参数，对被监测区域的应力积累变化提供有效参考，使用操作方便，实用性强，使用效果好，便于推广使用。

主权项

1.一种基于GNSS的地表变形连续监测方法，其特征在于，包括以下步骤：步骤一、建立GNSS观测站：在待监测区域的地表表面建立n个GNSS观测站，n为大于3的整数，任意三个GNSS观测站不共线，任意两个GNSS观测站的连线为观测基线，观测基线的条数共i条，步骤二、获取观测基线的长度变形量：n个GNSS观测站分别实时获取其所在位置处的位移量，并将获取到的n个GNSS观测站所在位置处的位移量分别实时传输给计算机(2)，计算机(2)根据接收到的n个GNSS观测站所在位置处的位移量计算得到i条观测基线的长度变形量，并将i条观测基线的长度变形量存储在观测基线长度变形量库e中，e＝{e₁,e₂,...,e_j,...,e_i}，其中，e_j表示第j条观测基线的长度变形量，j＝1、2、...、i；步骤三、计算观测基线的应变分量：计算机(2)根据公式计算和其中，表示第j条观测基线在x轴方向上的应变分量，表示第j条观测基线在y轴方向上的应变分量，表示第j条观测基线在xy坐标轴角平均线方向上的应变分量，θ_j表示第j条观测基线的方向角，κ_j表示第j条观测基线的应变分量系数；步骤四、计算观测基线的应变分量的平均值：计算机(2)根据公式计算和其中，表示i条观测基线在x轴方向上的应变分量的平均值，表示i条观测基线在y轴方向上的应变分量的平均值，表示i条观测基线在xy坐标轴角平均线方向上的应变分量的平均值；步骤五、计算待监测区域的应变参数：计算机(2)根据公式计算待监测区域的应变参数，待监测区域的应变参数包括第一剪应变参数r₁、第二剪应变参数r₂和面应变参数ε_α；步骤六、计算观测基线的剪切斜应变参数：计算机(2)根据步骤五中的第一剪应变参数r₁和第二剪应变参数r₂，获得i条观测基线的剪切斜应变参数，其中，i条观测基线中的第j条观测基线的剪切斜应变参数r(θ_j)的计算公式为：r(θ_j)＝r₁sin2θ_j+r₂cos2θ_j；步骤七、重复步骤二至步骤六，对待监测区域进行连续监测。