各水准点的高程；
[0042] 平差计算模型为：误差方程V = BX-L，V改正数矩阵，B系数矩阵，X未知数矩阵， L常数项矩阵；
[0043] 精度统计：统计平差单位权中误差与最弱点中误差；
[0044] Q23 :重力数据处理
[0045] Q231 :资料整理分析：
[0046] 相对重力测量数据预处理：外业资料汇集与整理，包括点位编号、汇总坐标与高 程，整理相对重力联测资料等内容；数据处理计算，包括仪器读数换算成毫伽值，进行气压 改正、仪器高改正、潮汐改正后得到初步观测值，经零漂改正后得到最后观测值，然后进行 重力段差计算、精度计算等；
[0047] 绝对重力测量资料的整理：全面收集和整理已有绝对重力测量成果，建立绝对测 量成果档案，整理与分析绝对重力测量资料；
[0048] Q232 :平差方案确定：
[0049] 原则：平差计算时，基本网采用"弱基准"，平差时不固定任何重力点；
[0050] 权：在平差计算过程中，绝对重力观测量、飞机联测相对重力观测量和所有长基线 标定相对重力观测量都始终采用先验权；其它相对重力观测量采用抗差估计重新定权；
[0051] 仪器参数：采取由多至少的原则，逐次舍去各仪器振幅较小或振幅与中误差相 当或相位中误差较大的周期误差，只顾及影响显著的周期项，以便达到仪器参数的合理匹 配；
[0052] Q233 :平差计算：
[0053] 观测量：绝对观测量和相对观测量；
[0054] 误差计算方程为：V = AX-L ;
[0055] 绝对重力观测量误差方程为R -g,°，&为i点的平差重力值，乂为i点的绝对 重力仪测定的观测重力值：
[0056] 相对重力观测误差计算方程为：
[0057]
[0058] gl、g]分别为测站i、j点平差后的重力值，g RZl、gRz]分别为测站i、j点经过四项改 正的相对联测最后观测值，Ri、Rj仪器在测站i、j点的观测读数，Ck重力仪的M次多项式 格值函数的K次格值改正因子，X n、Yn周期误差的参数，T n周期误差的周期；
[0059] 精度评定：单位权中误差m。的计算公式为
t观测量的总个数，t 必要观测量的总个数；
[0060] Q3 :多网并置资源的应用：将上述步骤得到的数据进行建立全国各地区厘米级高 精度似大地水准面模型、大陆速度场模型、全国地面沉降动态模型，应用于现代大地基准的 建立与维护。
[0061] 作为上述技术方案的优选，所述Q2步骤中对GNSS数据处理采用美国麻省理工学 院研制的高精度数据处理软件GAMIT/GL0BK或者BERNESE软件。
[0062] 作为上述技术方案的优选，所述Q21、Q22与Q23步骤不分前后，同时进行。
[0063] 作为上述技术方案的优选，所述Q22步骤中水准数据处理采用的步骤流程为：
[0064] 根据二等水准测量规范，采集水准网数据；
[0065] 进行标尺改正；初步计算测站的概略高程值；
[0066] 然后进行水准观测数据的质量控制，剔除粗差，控制系统误差影响；
[0067] 利用自由网平差、拟稳平差方法进行平差，评定计算；最后生成水准网成果及数据 处理报告。
[0068] 作为上述技术方案的优选，所述Q23步骤中重力数据处理采用的步骤流程为：
[0069] 按照《国家重力控制测量规范》（GB/T20256-2006)，采集重力网数据；
[0070] 整理分析观测数据；
[0071] 进行数据质量控制，剔除粗差，控制系统误差影响；
[0072] 进行气压、仪器高、潮汐、零漂等改正；确定合理有效的平差方案，对重力网进行平 差，评定计算；最后生成重力网成果及数据处理报告。
[0073] 本发明实施例提供的一种新型大地基准建网方法，与传统的技术相比，结构简单， 通过GNSS、水准、重力网并置，可以建立高精度似大地水准面、大陆速度场模型、地面沉降模 型，精度提高到厘米级，实现现代大地基准建立与维护，进一步丰富地质灾害与地理国情监 测手段，达到增效的目的。
【附图说明】
[0074] 为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现 有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本 发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可 以根据这些附图获得其他的附图。
[0075] 图1为本发明实施例1的一种新型大地基准建网方法的步骤流程示意图。
[0076] 图2为本发明实施例1的一种新型大地基准建网方法的Q22步骤中水准数据处理 的步骤流程示意图。
[0077] 图3为本发明实施例1的一种新型大地基准建网方法的Q23步骤中重力数据处理 采用的步骤流程示意图。
[0078] 图4为本发明实施例1的一种新型大地基准建网方法的Ql步骤中GNSS水准重力 观测墩的结构示意图。
[0079] 图5为本发明实施例2的一种新型大地基准建网方法应用的高程控制模型的结构 示意图。
[0080] 图6为本发明实施例3的一种新型大地基准建网方法的似大地水准面模型的结构 示意图。
[0081] 图7为本发明实施例4的一种新型大地基准建网方法的全国地面沉降动态模型的 结构示意图。
【具体实施方式】
[0082] 下面将结合本发明的附图，对本发明的技术方案进行清楚、完整地描述，显然，所 描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例， 本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发 明保护的范围。
[0083] 实施例1
[0084] 如图1-4所示，本发明实施例提供的一种新型大地基准建网方法，如下步骤：
[0085] GNSS水准重力观测墩建设
[0086] 如图4所示，按照GNSS水准点观测墩建设要求，在符合条件站址建立观测墩。按 照《国家重力控制测量规范》要求，重力观测墩建设在GNSS水准点观测墩附近，且重力观测 墩基座与GNSS观测墩墩体分离建设。
[0087] GNSS观测按照《全球定位系统（GPS)测量规范》（GB/T 18314-2009)进行观测，水 准测量按照《国家一、二等水准测量规范》（GB 12897-2006)测量。重力测量按照《国家重 力控制测量规范》进行观测数据采集。
[0088] GNSS、水准、重力网数据处理
[0089] GNSS数据处理
[0090] 1)软件
[0091] 采用美国麻省理工学院研制的高精度数据处理软件GAMIT/GL0BK或者Bemese软 件。
[0092] 2)数据预处理
[0093] 数据整理：依据外业观测手簿，将同一天的观测数据放在一起，数据格式为Rinex 格式，并进行数据正确性的检验。
[0094] 国家CORS站数据收集：国家GPS连续运行站是我国大地测量高精度的观测站，处 理国家GPS连续运行站及周围IGS站观测数据。
[0095] 数据标准化：使用随机软件标准化，形成观测数据文件和广播星历文件。
[0096] 天线高归算：按照天线结构，天线高统一采用观测值归算。
[0097] 3)基线解算
[0098] 先验坐标获取：差分模式获取GPS观测站的先验坐标，其坐标可以达到0.1 m以内 的精度。
[0099] 参数设置：设置高度角、对流层、观测值类型等其它参数。
[0100] 基线解算：以年积日为单位，进行基线解算。
[0101] 基线重复性检验：基线重复性按下式计算
[0102]
[0103] 其中，η为同一基线的观测时段总数