[发明专利]一种对北斗连续运行参考站覆盖区域构建应变场的方法

权利要求书

1.一种对北斗连续运行参考站覆盖区域构建应变场的方法，其特征在于：包括如下步骤：

步骤一：获取北斗连续运行参考站观测数据进行预处理，并利用高精度长基线数据处理软件解算北斗数据，对粗差进行踢除，获得所有站点的位移场；

步骤二：根据北斗连续运行参考站位置分布，构建站点的Delaunay三角网；

步骤三：用坐标方式构建每个三角形顶点的索引；

步骤四：利用步骤一解算得到的每个站点的位移值给对应三角形顶点赋值，建立位移场和顶点坐标的二元一次函数关系；

式中，(x,y)代表三角形的顶点位置；(u,v)代表北斗观测站点获取的位移值；

步骤五：在任一三角形中，根据步骤四公式，将三个站点的位移和站点坐标建立函数关系：

步骤六：将所有站点的位置信息和位移场信息带入步骤五建立的函数关系中，解算得到系数a₁,a₂,a₃和b₁,b₂,b₃；

步骤七：利用步骤六解算得到的参数，解算得到北斗连续运行参考站覆盖范围内的主应变率；

步骤八：计算最大主应变率和最小主应变率；

步骤九：计算最大剪切应变率和面膨胀率。

2.根据权利要求1所述的一种对北斗连续运行参考站覆盖区域构建应变场的方法，其特征在于：步骤七中的主应变率计算方法如下：

其中，ε_x、ε_y表示三角形单元内x方向和y方向的主应变率，ε_xy代表三角形单元在oxy平面内的剪切应变率。

3.根据权利要求2所述的一种对北斗连续运行参考站覆盖区域构建应变场的方法，其特征在于：步骤八中最大主应变率和最小主应变率计算方法如下：

其中，ε₁表示研究区域内最大主应变率，ε₂表示最小主应变率，为ε₁的方向，ε₂方向同理也可求出。

4.根据权利要求3所述的一种对北斗连续运行参考站覆盖区域构建应变场的方法，其特征在于：步骤九中最大剪切应变率和面膨胀率计算方法如下：

△＝ε₁+ε₂＝ε_e+ε_n

其中，γ_max表示区域最大剪切应变率，△表示面膨胀率。