买专利卖专利找龙图腾，真高效！ 查专利查商标用IPTOP,全免费！专利年费监控用IP管家,真方便！

龙图腾网恭喜玉溪矿业有限公司申请的专利减小矿区测量高程异常差值的方法获国家发明授权专利权，本发明授权专利权由国家知识产权局授予，授权公告号为：CN115406401B 。

龙图腾网通过国家知识产权局官网在2025-03-07发布的发明授权授权公告中获悉：该发明授权的专利申请号/专利号为：202211139312.1，技术领域涉及：G01C5/00；该发明授权减小矿区测量高程异常差值的方法是由俞艳波;李光全;李小松;杨双龙;李琦;苏海华;普庆红设计研发完成，并于2022-09-19向国家知识产权局提交的专利申请。

本减小矿区测量高程异常差值的方法在说明书摘要公布了：本申请公开了一种减小矿区测量高程异常差值的方法，包括以下步骤：步骤S1：从观测重力异常Δg中去除模型重力异常ΔgGM和地形重力效应δNT得到残差重力异常，将残差重力异常代入Stokes公式，按下式计算残差大地水准面δNr：步骤S2：在缺乏重力资料的情况下，根据高程异常的波谱结构特点分解大地水准面高N；步骤S3：在未知点上，由EGM2008地球重力场模型计算出未知点的模型高程异常，在所得结果上加上未知点剩余高程异常ξC，得到未知点最终的高程异常值，求得未知点的正常高。通过实践，基于“移去—恢复”法的多面函数拟合方法在大红山铜矿矿区的拟合效果较好，拟合后正常高精度可达到±2cm的要求。

本发明授权减小矿区测量高程异常差值的方法在权利要求书中公布了：1.一种减小矿区测量高程异常差值的方法，其特征在于，包括以下步骤：步骤S1：重力异常观测值分解为：由位模型计算得到的模型重力异常，地形重力效应和残差重力异常，不考虑地形的影响，以上分解可以分成两部分，即由位模型计算得到的模型重力异常部分和残余部分，这两种分解方式可以分别表示为：Δg＝ΔgGM+δΔgT+δΔgr其中，Δg为观测重力异常，ΔgGM为模型重力异常，δΔgT为局部地形改正值，δΔgr为地形质量对计算点的引力；N＝NGM+δNT+δNr其中，N为重力异常观测值，NGM为模型重力异常，δNT为地形重力效应，δNr为残差重力异常；步骤S11：移去过程，按下式将模型重力异常ΔgGM和地形重力效应δNT从观测重力异常Δg中去除得到残差重力异常δΔgr:δΔgr＝Δg-ΔgGM-δgT式中，模型重力异常ΔgGM采用位系数法按下式计算得到: 上式中，ψ，λ为计算点的纬度和经度，GM是地心引力常数，R为地球平均半径，为完全规格化位系数，为完全规格化缔合Legendre函数，N为位模型球谐展开的最大阶数；a为拟合系数；δΔgr为地形质量对计算点的引力，仅考虑相对于布格片的局部地形起伏的影响，即局部地形校正δΔgTC，考虑相对于大地水准面的地形起伏的影响，计算公式为： 其中，δNtotal为大地水准面地形起伏影响值，Gρ为计算点的地心距离，γ为平均正常重力，hp为计算点的地形高，步骤S12：将残差重力异常代入Stokes公式，按下式计算残差大地水准面δNr： 其中，R为地球平均半径，γ为平均正常重力，Δgr为地形质量对计算点的引力，Sψ为计算点的纬度；步骤S13：进行大地水准面的恢复，分别计算NGM和δNT，其中NGM由位系数计算得到： 式中，ρ为计算点的地心距离，γ为平均正常重力；δNT按下式计算： 其中，δNT为局部地形影响，ρ为地壳密度，hp为计算点的地形高，h为流动点的地形高，l为计算点至流动点的距离，G为地心距离；步骤S2：在缺乏重力资料的情况下，根据高程异常的波谱结构特点，将大地水准面高N分解为：ξ＝ξGM+ξΔG+ξT其中，ξΔG为中波部分，称为残差高程异常，ξT为短波部分地形高；ξGM为长波部分，称为模型高程异常，采用EGM2008重力场模型按下式计算ξGM： 其中，ψ，λ为计算点的纬度和经度，GM是地心引力常数，R为地球平均半径，为完全规格化位系数，为完全规格化缔合Legendre函数，N为位模型球谐展开的最大阶数；a为拟合系数；ρ为计算点的地心距离，γ为平均正常重力，ξΔG表示中波部分，称为残差高程异常，通过求解剩余重力异常的边值得到： 其中，R为地球平均半径，γ为平均正常重力，为计算点的纬度、经度及地心距离值，为高程异常；为点的纬度；ξT——表示短波部分，通过求解地形高正得到： 其中，ρ为计算点的地心距离，γ为平均正常重力，hp为计算点的地形高，l为计算点至流动点的距离，h为流动点的地形高，G为点的地心距离；在没有数字高程模型数据情况下，将ξΔG和ξT部分合并后用数学模型逼近的方法表征，将GPS点的高程异常分为两部分求解，即：ξ＝ξGM+ξC式中，ξGM——重力场求得的高程异常；ξC——剩余高程异常；步骤S3：将高程异常移去EGM2008模型大地水准面高的剩余值认为是剩余高程异常值ξC；基于EGM2008模型的高程异常值，并以GPS大地高和三等水准测量得到的正常高求差，得到各已知点的真实高程异常值，剩余高程异常值由模型高程异常值与真实高程异常值求差得到，并认为是各已知点的真实剩余高程异常值；结合各已知点的CGCS2000坐标值，将剩余高程异常值作为已知值，分别代入到加权平均模型、多项式曲线模型、平面模型、二次曲面模型、三次曲面模型、移动曲面模型和多面函数模型中，建立起以CGCS2000坐标值为参数，求解各未知点的剩余高程异常值函数模型；各未知点基于EGM2008模型，计算出各点的模型高程异常，最后得到各未知点的高程异常值；计算公式如下： 式中，ξxi,yi——未知点的高程异常值；ξCxi,yi——未知点基于函数模型的剩余高程异常插值；ξGMxi,yi——未知点基于EGM2008的模型高程异常；Ffitxi,yi——剩余高程异常值插值函数；FEGM2008xi,yi——EGM2008模型高程异常插值函数；求得未知点的正常高。

如需购买、转让、实施、许可或投资类似专利技术，可联系本专利的申请人或专利权人玉溪矿业有限公司，其通讯地址为：653100 云南省玉溪市红塔区珊瑚路84号；或者联系龙图腾网官方客服，联系龙图腾网可拨打电话0551-65771310或微信搜索“龙图腾网”。