长距离水工隧洞基于GNSS网坐标协因数阵的隧洞横向贯通误差预计法.docxVIP

附录 A （规范性） 基于GNSS网坐标协因数阵的隧洞横向贯通误差预计法 A.1 基于GNSS网坐标协因数阵严密计算隧洞横向贯通误差的原理 设隧洞工程独立坐标系与隧洞贯通面的相对关系如图A.1所示，将隧洞工程独立坐标系顺时针旋转角，使得旋转后的轴与隧洞贯通面垂直，旋转后的坐标系为。 图A.1 隧洞独立坐标系与隧洞贯通面相对位置关系示意图 由图A.1中各点的几何关系可得贯通点P从进洞推算的坐标与从出洞推算坐标的坐标差表达式为： （A.1） 式中： 、——分别为由洞口点A、D推算贯通点P的坐标；——分别为AB、DE边的方位角。 不考虑联系测量中定向角以及洞内导线测量边长的测量误差（因为本节仅考虑洞外控制测量误差引起的横向贯通误差），可得贯通面的横向贯通中误差的微分表达式为： （A.2） 顾及到下式（A.3）及式（A.4）： （A.3） （A. 4） 整理可得到横向贯通中误差的详细微分表达式： （A.5） 其中，及如下式（A.6）及式（A.7）所示： （A. 6） （A.7） 由误差传播定律可知，横向贯通中误差为下式： （A.8） 式中： ——洞外控制网的验后单位权中误差； ——坐标协因数矩阵。 通过上面一系列的推导可知，当隧洞洞外控制网的和可知时，其横向贯通中误差就可以按照上式（A.8）严密计算得到。但是上文中提到的一系列商用软件一般都不提供验后坐标协因数阵，如何才能得到这两个必要的参数，就成了严密计算GNSS网隧洞横向贯通误差的关键所在。 A.2 一点一方向平差中单位权中误差和坐标协因数阵的计算方法 前文已经提到，由于长距离水工隧洞对施工精度和贯通误差的严格要求，需要建立工程独立控制网，采用GNSS建立隧洞工程独立控制网的步骤如下： （1）结合具体的隧洞工程独立控制网的建网需求，挑选出一个适宜的地面固定点及与其相关的一条特定方向。 （2）根据前面介绍的方法，选择合适的中央子午线与工程投影面大地高，进行GNSS工程控制网一点一方向平差。平差时先在坐标系中进行三维控制网无约束平差，再将三维无约束平差结果转换至高斯平面，最后在高斯平面上通过平移旋转即可得到工程独立坐标系中的成果。 设固定地面起算点的平面坐标为（），其同一坐标系中的特定起算方向方位角为，二者均为隧洞独立坐标系中的成果；该点对应的三维无约束平差后转换得到的高斯平面直角坐标为（），由三维无约束平差成果转换至高斯平面直角坐标系，经反算得到的该方向坐标方位角为，则洞外控制网一点一方向平差后的高斯平面直角坐标经平移旋转后在独立坐标系中的坐标为： （A.9） 式中： 、——GNSS控制点在隧洞独立坐标系中的平面坐标； 、——经过三维无约束平差后转换得到的高斯平面直角坐标系中的坐标。 基于上述一系列的平差转换思想，根据方差-协方差传播定律，可以得到变换后坐标的协方差阵为： （A.10） 式中： ——经三维无约束平差转换，计算得到的高斯平面坐标中的相应方差阵； ——分块对角阵，其主对角线上的矩阵子块为： （A.11） 通过工程独立控制网一点一方向平差，可以直接得到隧洞独立控制网的验后单位权中误差，又通过式（A.9）和式（A.10）可以计算得到隧洞独立坐标系下的协因数矩阵，由此便可通过式（A.8）严密计算得到GNSS误差引起的隧洞横向贯通中误差。除了上述方法，隧洞独立坐标系协因数阵还可通过三维和二维坐标系中的一系列转化得到，下一节即详细介绍该方法。 A.3 GNSS网坐标协因数阵计算的另外一种方法 GNSS控制网基线所在的坐标系为WGS84空间直角坐标系，GNSS测量得到的基线向量经平差、转换可得到大地坐标系坐标；再经过高斯投影可得到高斯平面直角坐标系坐标；最后经转换可得隧洞独立坐标系坐标。隧洞洞外GNSS控制网的坐标协因数阵亦可通过上述一系列三维和二维坐标转化得到，下面介绍其详细计算过程。 （1）GNSS控制网空间直角坐标协因数阵的计算 洞外GNSS控制网以空间直角坐标系中基线向量为观测值，其误差方程可表达如下： （A.12） 式中： 等——近似坐标； 等——近似坐标相应改正数； ——基线分量残差。 将所有独立基线向量的观测方程组合成矩阵形式即为： （A.13） 式中： ——残差项； ——系数阵； ——改正数项； ——常数项。 使用间接平差原理可分别得到验后单位权中误差以及空间直角坐标协因数阵： （A.14） 式中： ——多余观测的个数； ——相应权矩阵。 （2）GNSS控制网大地坐标协因数阵的计算 GNSS控制网空间直角坐标与大地坐标的微分关系式如下所示： （A.15） 式中： L、B——大地经度、纬度； h——大地高； M、N——分别为子午圈、卯酉圈曲率半径。 若忽略GNSS控制网中各点间的相关性，并令为分块矩阵