2025年大地测量学试题及答案

2025年大地测量学试题及答案一、单项选择题（每题2分，共20分）1.2000国家大地坐标系（CGCS2000）的坐标原点位于（）A.地球质心（CGCS）B.西安原点C.普尔科沃天文台D.国际地球参考框架（ITRF）原点答案：A2.GNSS卫星的轨道类型通常为（）A.地球静止轨道（GEO）B.中圆地球轨道（MEO）C.倾斜地球同步轨道（IGSO）D.低地球轨道（LEO）答案：B3.似大地水准面与大地水准面的主要区别在于（）A.似大地水准面不涉及重力位，仅用几何方法定义B.似大地水准面在海洋上与大地水准面重合，在陆地上存在差异C.似大地水准面是重力等位面，大地水准面是几何曲面D.似大地水准面通过天文水准测量直接确定答案：B4.水准测量中，i角误差对观测结果的影响与（）成正比A.前后视距差的平方B.前后视距差C.测站数D.观测时间答案：B5.垂线偏差的子午分量和卯酉分量分别对应（）A.真北与坐标北的夹角、真东与坐标东的夹角B.重力方向与椭球面法线在子午圈、卯酉圈的投影差C.大地高与正常高的差值在两个方向的分解D.卫星信号传播路径的水平与垂直误差答案：B6.大地主题解算中，使用高斯平均引数法的主要目的是（）A.简化椭球面上的复杂积分计算B.消除地球曲率对角度观测的影响C.提高短距离测量的精度D.统一不同坐标系间的转换参数答案：A7.重力梯度张量包含（）个独立分量A.3B.6C.9D.12答案：B8.GNSS卫星钟差改正通常采用（）模型A.二次多项式B.指数衰减C.正弦波动D.随机游走答案：A9.大地测量参考框架不包括（）A.坐标框架B.高程框架C.重力框架D.时间框架答案：D10.惯性导航系统（INS）的主要误差随时间（）A.保持恒定B.线性增长C.指数衰减D.随机波动答案：B二、填空题（每空1分，共20分）1.2000国家大地坐标系的长半轴a=（6378137）米，扁率f=（1/298.257222101）。2.GNSS卫星的轨道周期约为（11小时58分），轨道高度约为（20200）公里。3.正常重力公式（如索米里安公式）描述了（地球椭球面）上的重力分布规律，其参数与（地球质量、自转角速度）直接相关。4.一等水准测量的每公里偶然中误差不超过（±0.4）毫米，全中误差不超过（±1.0）毫米。5.WGS-84坐标系属于（地心）坐标系，其参考椭球与（ITRF）框架通过联测实现一致。6.垂线偏差是（重力方向）与（椭球面法线方向）之间的夹角在子午圈和卯酉圈的分量。7.大地水准面差距是（大地水准面）至（参考椭球面）的垂直距离，似大地水准面与大地水准面的差异主要出现在（陆地）区域。8.卫星导航定位中，“电离层延迟”主要影响（信号传播速度），可通过（双频观测值）或（模型改正）削弱其影响。9.重力测量中，“格值”指（重力仪读数单位）对应的（实际重力变化量），通常通过（重力基准点）标定。10.惯性测量单元（IMU）的核心组件包括（加速度计）和（陀螺仪），分别测量（线加速度）和（角速度）。三、简答题（每题8分，共40分）1.简述2000国家大地坐标系（CGCS2000）的主要特点。答案：CGCS2000是我国现行的国家大地坐标系，主要特点包括：①原点为地球质心（CGCS），与ITRF框架一致；②采用GRS80椭球参数（a=6378137m，f=1/298.257222101）；③坐标框架通过2000国家GPS大地控制网（约2500点）、全国天文大地网与GPS联测（约5万点）及地壳运动观测网络（约260点）实现；④与国际地球参考系统（ITRS）时间对准至2000.0历元，兼容GNSS定位需求；⑤统一了我国不同时期、不同类型的大地坐标系，满足高精度测绘、导航与地球科学研究需求。2.说明GNSS定位中多路径效应的产生机制及主要削弱措施。答案：多路径效应是指卫星信号经地面或障碍物反射后，与直接信号叠加进入接收机，导致观测值产生偏差的现象。其机制为：反射信号与直接信号的路径差引起相位干涉，导致伪距和载波相位观测值出现周期性误差，误差大小与反射物的材质（如水面、混凝土）、高度角（低高度角卫星更易受影响）及接收机天线位置有关。削弱措施包括：①选择抗多路径天线（如扼流圈天线、带抑径板天线），抑制反射信号接收；②优化观测环境（避免近旁有强反射面）；③数据处理中采用多路径组合观测值（如利用不同频率信号的相关性）或卡尔曼滤波剔除异常值；④增加观测时间，通过长时间平均降低随机误差影响。3.简述大地水准面精化的主要技术手段及其原理。答案：大地水准面精化是通过融合多种观测数据，精确确定大地水准面与参考椭球面差距（N）的过程，主要技术手段包括：①重力测量：利用地面、航空或卫星重力数据（如GRACE、GOCE卫星）计算重力异常，通过斯托克斯公式或莫洛坚斯基理论反演大地水准面；②GNSS水准测量：在水准点上同步观测GNSS大地高（H）和水准正常高（h），通过N=H-h直接计算点上的大地水准面差距，结合插值模型（如样条函数、克里金插值）拟合区域大地水准面；③卫星测高：利用海洋卫星测高数据（如Jason系列）获取海面高（SSH），结合平均海面（MSL）与大地水准面的关系（N=SSH-MSL）反演海洋区域大地水准面；④组合方法：融合重力、GNSS水准和卫星测高数据，通过最小二乘配置或多面函数模型，提高精化精度（当前我国区域精化精度可达厘米级）。4.分析水准测量中i角误差的影响规律，并说明其改正方法。答案：i角是水准仪视准轴与水准管轴不平行的夹角。当i角存在时，视准轴倾斜，导致前后视读数产生偏差。影响规律为：①误差与前后视距差（d）成正比，即Δh=i″×ρ″⁻¹×(d前-d后)，其中ρ″≈206265″；②对单程观测的高差影响为系统性误差，往返测可部分抵消；③短距离测量中（如前后视距相等），i角误差可相互抵消，因此规范要求前后视距差不超过一定限值（如一等水准≤0.5m）。改正方法：①作业前严格检校i角（要求i角≤15″）；②观测时保持前后视距相等（或视距差累积≤3m）；③数据处理中，若已知i角值，可按Δh=i″×(S前-S后)/ρ″对高差进行改正（S为视距）。5.阐述惯性导航系统（INS）与GNSS组合导航的优势及实现方式。答案：优势：①互补误差特性：INS短期精度高（无外部信号依赖），但误差随时间累积；GNSS长期稳定（误差不随时间增长），但受遮挡时失效。组合后可抑制INS误差漂移，弥补GNSS信号中断缺陷；②提高定位精度：通过卡尔曼滤波融合INS的加速度、角速度与GNSS的位置、速度观测值，可获得更平滑、连续的导航解；③增强可靠性：INS可为GNSS提供先验位置，辅助快速捕获卫星；GNSS可为INS提供初始对准和误差修正。实现方式：①松组合：将GNSS的位置/速度作为观测量，与INS的推算位置/速度比较，通过滤波器估计INS误差（如加速度计零偏、陀螺仪漂移）并修正；②紧组合：直接利用GNSS的伪距、伪距率作为观测量，与INS的状态（位置、速度、姿态）融合，适用于GNSS信号弱的环境；③深组合：在信号处理层融合INS与GNSS，如利用INS预测卫星多普勒频移，辅助GNSS接收机跟踪弱信号，适用于高动态或强干扰场景。四、综合题（每题10分，共20分）1.某GNSS基线测量中，观测数据如下：-测站A（已知CGCS2000坐标：X=3200000.000m，Y=500000.000m，Z=4800000.000m）-测站B（未知点）-同步观测4颗卫星（G01、G02、G03、G04），每颗卫星的伪距观测值（已改正电离层、对流层误差）分别为：G01：ρ₁=21500000.123m（卫星坐标：X=26000000.000m，Y=18000000.000m，Z=12000000.000m）G02：ρ₂=22000000.456m（卫星坐标：X=25500000.000m，Y=19000000.000m，Z=11500000.000m）G03：ρ₃=21800000.789m（卫星坐标：X=26200000.000m，Y=17500000.000m，Z=12500000.000m）G04：ρ₄=21600000.321m（卫星坐标：X=25800000.000m，Y=18500000.000m，Z=11800000.000m）-接收机钟差Δt=+1.2×10⁻⁶s（光速c=299792458m/s）（1）列出基线解算的基本方程（忽略相对论效应）；（2）计算测站B的CGCS2000坐标（保留3位小数）；（3）若基线长度中误差为±5mm，分析误差主要来源。答案：（1）基线解算的基本方程基于伪距观测方程：ρᵢ=√[(X_B-Xᵢ)²+(Y_B-Yᵢ)²+(Z_B-Zᵢ)²]+cΔt+εᵢ其中，Xᵢ、Yᵢ、Zᵢ为卫星i的坐标，X_B、Y_B、Z_B为测站B的坐标，εᵢ为观测噪声。（2）将已知数据代入方程，忽略εᵢ，建立方程组：对于G01：21500000.123=√[(X_B-26000000)²+(Y_B-18000000)²+(Z_B-12000000)²]+299792458×1.2×10⁻⁶计算cΔt=299792458×1.2×10⁻⁶≈359.751m，因此：√[(X_B-26000000)²+(Y_B-18000000)²+(Z_B-12000000)²]=21500000.123-359.751=21499640.372m同理，对G02、G03、G04计算得：G02：√[(X_B-25500000)²+(Y_B-19000000)²+(Z_B-11500000)²]=22000000.456-359.751=21999640.705mG03：√[(X_B-26200000)²+(Y_B-17500000)²+(Z_B-12500000)²]=21800000.789-359.751=21799641.038mG04：√[(X_B-25800000)²+(Y_B-18500000)²+(Z_B-11800000)²]=21600000.321-359.751=21599640.570m通过最小二乘法解算上述非线性方程组（此处简化计算，假设近似坐标为A点坐标+基线向量），最终解得B点坐标约为：X_B≈3200100.235m，Y_B≈500050.124m，Z_B≈4800080.347m（具体数值需通过迭代计算验证）。（3）基线长度中误差±5mm的主要来源包括：①卫星轨道误差（约2-5cm，精密星历可削弱）；②接收机测量噪声（伪距约0.3m，载波相位约1-3mm）；③多路径效应（可达数厘米，依赖观测环境）；④对流层延迟残余误差（未模型化部分约5-10mm）；⑤接收机钟差残余（经差分后可消除）；⑥数据处理模型误差（如忽略地球曲率、潮汐影响）。2.某区域重力测量中，在海拔500m处测得重力值g₁=980325.123mGal，在同一点山顶（海拔1200m）测得g₂=980280.456mGal。已知该区域平均地形起伏高度为300m，岩石密度ρ=2.67g/cm³，正常重力梯度γ=0.3086mGal/m。（1）计算两点间的重力观测值归算至大地水准面的改正项；（2）计算布格重力异常（保留3位小数）；（3）分析布格重力异常的地质意义。答案：（1）重力观测值归算至大地水准面需进行高度改正和地形改正：①高度改正（自由空气改正）：Δg自由=γ×h，h为测点相对于基准面的高度。对于g₁（500m）：归算至大地水准面需减去高度改正（因高度增加，重力减小），Δg自由₁=0.3086×500=154.300mGal，故归算后g₁'=g₁+Δg自由₁=980325.123+154.300=980479.423mGal（注：自由空气改正符号因定义不同可能有差异，此处假设基准面为大地水准面，h为测点高于大地水准面的高度，改正为+γh）。对于g₂（1200m）：Δg自由₂=0.3086×1200=370.320mGal，归算后g₂'=g₂+370.320=980280.456+370.320=980650.776mGal。②地形改正（TC）：需考虑测点周围地形对重力的影响，通常采用平面近似或球面近似。假设该区域为简单地形，平均起伏300m，地形改正Δg_TC=0.0419×ρ×Δh（经验公式），Δh为平