曲网格克里金网格化方法及其在重力测量近区地形改正中的应用

㊀第３９卷第４期物㊀探㊀与㊀化㊀探
Ｖｏｌ．３９，Ｎｏ．４㊀㊀２０１５年８月
ＧＥＯＰＨＹＳＩＣＡＬ＆ＧＥＯＣＨＥＭＩＣＡＬＥＸＰＬＯＲＡＴＩＯＮ
Ａｕｇ．，２０１５㊀
ｄｏｉ：１０．１１７２０／ｗｔｙｈｔ．２０１５．４．３１唐小平，刘宽厚，耿涛，等．曲网格克里金网格化方法及其在重力测量近区地形改正中的应用［Ｊ］．物探与化探，２０１５，３９（４）：８４８－８５４．ｈｔｔｐ：／／ｄｏｉ．ｏｒｇ／１０．１１７２０／ｗｔｙｈｔ．２０１５．４．３１
ＴａｎｇＸＰ，ＬｉｕＫＨ，ＧｅｎｇＴ，ｅｔａｌ．Ｃｕｒｖｅｄ⁃ｇｒｉｄｋｒｉｇｉｎｇｇｒｉｄｄｉｎｇｍｅｔｈｏｄａｎｄｉｔｓａｐｐｌｉｃａｔｉｏｎｔｏｎｅａｒ⁃ｓｔａｔｉｏｎｔｅｒｒａｉｎｃｏｒｒｅｃｔｉｏｎｉｎｇｒａｖｉｔｙｍｅａｓｕｒｅｍｅｎｔ［Ｊ］．Ｇｅ⁃ｏｐｈｙｓｉｃａｌａｎｄＧｅｏｃｈｅｍｉｃａｌＥｘｐｌｏｒａｔｉｏｎ，２０１５，３９（４）：８４８－８５４．ｈｔｔｐ：／／ｄｏｉ．ｏｒｇ／１０．１１７２０／ｗｔｙｈｔ．２０１５．４．３１
曲网格克里金网格化方法
及其在重力测量近区地形改正中的应用
唐小平１，２，刘宽厚１，耿涛１，冯治汉１
（１．中国地质调查局西安地质调查中心，陕西西安㊀７１００５４；２．国土资源部岩浆作用与找矿重点实验室，陕西西安㊀７１００７４）
㊀㊀收稿日期：２０１４⁃０６⁃２５；修回日期：２０１５⁃０３⁃２７
㊀㊀基金项目：国家青年基金项目（４１４０４１１０）；中国地质调查局项目（１２１２０１１２２０６００㊁１２１２０１１３０３１６００）；陕西省面上基金项目（２０１５ＪＭ４１３３）
摘要：文中将研究一种高精度㊁强稳定性与灵活性的近区地改值计算方法㊂首先将克里金网格化技术与曲网格方法相结合，研究曲网格克里金网格化方法（简写为ＣＧＫＧ⁃Ｍ）及其扩边技术，进而实现数据的圆域网格化；再将该方法应用于圆域近区地改值的计算中，实现任意方位㊁任意环数的圆域近区地改值计算；最后进行了模型与实际地形的测试和误差分析㊂研究结果表明，ＣＧＫＧ⁃Ｍ近区地改算法具有方法简单㊁计算精度高㊁灵活性与稳定性强㊁计算速度快等特点，是与ＧＴＣＳ⁃１型近区地改仪数据采集格式相匹配的一种高效㊁快捷的近区地改值计算方法㊂关键词：曲网格克里金网格化法；圆域；近区地形改正
中图分类号：Ｐ６３１㊀㊀㊀文献标识码：Ａ㊀㊀㊀文章编号：１０００－８９１８（２０１５）０４－０８４８－０７
㊀㊀近区地形改正是重力勘探的一项重要工作，传统近区地改值的计算方法主要有基于扇区（圆域）与网格（方域）两种算法［１－２］，但均有所不足㊂基于网格化的方域算法在计算时，首先将矩形框区域网格化并求出总地改值，再减去四个角部分地改值，从
而得出实际地改值㊂当地形特别复杂或地形数据较少时，这种算法精度较低，容易引起畸变㊂圆域（扇区）算法是将计算区域按方位和环数分成若干个扇区，再求取每个扇区的平均地改值，最后计算合成后的总地改值㊂该方法在求取扇区平均高差，要求每个扇区参与计算平均高差的点数不少于５个，计算时若扇区划分太小或数据量过小或分布不均匀时，则难以满足上述要求；若扇区太大，测点周围的地形描述精度将会整体下降，从而影响地改精度㊂总的来讲，方域算法在消去远端四角地形影响时存有一定误差，圆域算法则在方法灵活性与地形数据分布等方面受到限制㊂故而急需研究一种高精度㊁强稳定性与灵活性的近区地改值计算方法，以克服传统
近区地改中存在的诸多问题㊂
近区地改算法最早是圆域地改算法，由Ｈａｍｍｅｒ于１９３９年提出［３］，并沿用至今㊂随着计算机技术的发展，上世纪６０年代开始出现基于计算机软件开发的方域网格化算法［４］，并在此后的发展中，该方法在重力近区㊁中区地改中得到了很大的发展和应用［５－６］㊂一直以来，地形数据精度及计算方法对重力地形改正的影响受到了广泛关注，冯治汉［７］使用了直立长方体法，并采用不同网格间距的高程数据计算北祁连地区的重力调查地形改正值，并系统地讨论了地改精度，得出越小网格距的高程数据对地形的模拟程度越高；张俊等［８］为了得出高程数据网格间距对表面积分法㊁直立长方体法和平均高程直立长方体积法计算中区地形改正值精度的影响，对选择的４５０测点，采用不同网格间距的高程数据，以中区分段的方式比较了三种方法的地改精度，得出表面积分法受高程数据网格化间隔的影响较小，直方体受到的影响最大，同时指出影响中区地
㊀４期唐小平等：曲网格克里金网格化方法及其在重力测量近区地形改正中的应用
改精度的最大影响地段为５０ ５００ｍ㊂，随着现在光学与测量技术的发展，在近区地改方面，为了解决地形测量的精度问题，近年来人们开始研制基于现代光学测量仪器的近区地形改正仪器系统㊂Ｓｃｈｉａｖｏｎｅ
Ｄ等对摄影制图㊁数字化地图和高精度的激光扫描数据进行了比较，得出了激光扫描数据精度最高㊁地形重建时间最短的结论［９］，不久又研制出了基于激光扫描技术的近区地改仪［１０］，但该仪器价格昂贵，比较笨重；刘宽厚㊁耿涛等研究和实验了基于激光测距技术下的近区地改，研制出了重力测量近区地形改正系统［１１］，该系统具有使用简单㊁仪器轻便㊁能适应复杂地形等特点，主要针对０ ２０ｍ内的近区地改㊂本算法即为该系统的配套算法之一㊂
曲网格技术是基于复杂地形与内部结构建模的一类数值建模方法㊂１９８５年，Ｔｈｏｍｐｓｏｎ等［１２］系统地阐述了曲网格技术在各类方程㊁各种模型中使用的基本理论；随后，这一技术被广泛的应用于地球物理的各个领域［１３－１４］，并在此基础上发展出了正交网格等更有效的计算网格［１５－１６］㊂
地球物理中的网格化技术很多，如线性插值法㊁多元二次函数㊁反插值法以及克里金网格化法等［１７］㊂克里金网格化法最早由南非采矿工程师ＫｒｉｇｅＤＧ提出，故命名为 克里金 法［１７］，后经Ｍａｔｈｅｒｏｎ［１８－１９］发展而逐步完善，是一种空间内插法㊂文中通过引入曲网格数据计算转化原理，先将圆域测量数据按方位与距离展开成方域，再在方域中应用克里金网格化法进行网格化计算，最后应用曲网格数据计算反转化原理还原到圆域中，计算出相应点高差，并采用圆域地改值计算方法计算出相应的地改值㊂
１㊀方法原理
１．１㊀圆的曲网格计算与极坐标系下克里金法在坐标转换中，极坐标系（θ，ｒ，ｈ）与直角坐标系（ｘ，ｙ，ｚ）之间的转换关系为
ｘ＝ｒｃｏｓθ，ｙ＝ｒｓｉｎθ，ｚ＝ｈ㊂（１）对式（１）引入曲网格计算，将极坐标转换为计算坐标（ε，η，ｚ）㊂设计算区域为：εɪ（０，１），ηɪ（０，１），得
ε＝θ２π，η＝ｒ－ｒ１
ｒ２－ｒ１
，ｚ＝ｈ㊂（２）
将上式进行变形，从而得到（θ，ｒ，ｈ）的反算公式θ（ε）＝２πε，ｒ（η）＝ｒ１＋（ｒ２－ｒ１）η，
ｈ＝ｚ㊂（３）将式（３）代入式（１）简化得到的泛函关系
ｘ（ε，η）＝ｒ（η）ｃｏｓ（ε），ｙ（ε，η）＝ｒ（η）ｓｉｎ（ε），
ｈ＝ｚ㊂（４）将式（４）代入克里金网格化方法中，求取曲网格克里金网格化的基本计算公式㊂常规克里金方法
Ｚ∗＝ｍ（ｕ）＋ðｎ（ｕ）ａ＝１［λａ（ｕ）－ｍ（ｕ）］，（５）则在曲网格下的计算公式即表示为
Ｚ∗＝ｍ（ｕ）＋ðｎ（ｕ）ａ＝１［λａ（ｕ）－ｍ（ｕ）］＝ｍ［ｕ（ｘ，ｙ，ｚ）］＋ðｎ（ｕ）ａ＝１｛λａ［ｕ（ｘ，ｙ，ｚ）］－
ｍ［ｕ（ｘ，ｙ，ｚ）］｝
＝ｍ｛ｕ［ｘ（ε，η），ｙ（ε，η），ｚ］｝＋
ðｎ（ｕ）ａ＝１｛λａ［ｕ（ｘ（ε，η），ｙ（ε，η），ｚ）］－
ｍ｛ｕ［ｘ（ε，η），ｙ（ε，η），ｚ）］｝㊂（６）式（６）即为最终曲网格克里金网格化法，该式实际控制坐标为（ε，η，ｚ），通过式（６）可以计算出任意给定的（ε，η）点上的坐标高差值ｚ，而通过式（３）可以求出（ε，η，ｚ）所对应的极坐标（θ，ｒ，ｈ），从而实现了曲网格克里金网格化方法在圆域网格化中的应用㊂
１．２㊀曲网格网格化方法在近区地改中的应用首先说明ＧＴＣＳ⁃１型仪器的数据格式，数据以重力测点为圆心，按方位和环数分布于测点四周（图１ａ），（０，０）为测点坐标，数据测量点数位于其周围，稀疏分布在半径２０ｍ的圆形内㊂将ＧＴＣＳ⁃１的测量数据转换为（θ，ｒ，ｈ）后，应用式（２），计算坐标系中的（ε，η，ｚ）
ε＝θ２π，η＝ｒ－ｒ１
ｒ２－ｒ１
，ｚ＝ｈ㊂（７）其中：θ表示方位角度；ｒ为平距；ｈ为高差；ｒ１㊁ｒ２表示地改外半径和内半径；在０ ２０ｍ近区地改中，ｒ２＝２０ｍ，ｒ１＝０．５ｍ，０ɤθɤ２π，０ɤεɤ１，０ɤηɤ１㊂经过式（２）可将（θ，ｒ，ｈ）转化为曲网格的计算域（ε，η，ｚ），代入式（６）中，即可实现任意（ε，η）下的ｚ值求取，进而调用反算公式（３），即可得到相应的（θ，ｒ，ｈ）㊂具体应用过程如图１ｂ 图１ｄ所示，图１ｂ为将ＧＴＣＳ⁃１的测量数据转换为（θ，ｒ，ｈ）后，按式（７）展布于计算空间之中，将图２ｂ代入式（６）中，对各节点的高差进行网格化（图１ｃ），最后通过式（４）即可回到物理空间，成标准的圆环分布（图
㊃９４８㊃
物㊀探㊀与㊀化㊀探３９卷㊀
１ｄ）㊂假如这些点刚好为该扇区的重心位置，则可计算出所在扇区的地改值，将所有扇区地改值叠加，即可得出最终地改值㊂
另外，从图１ｃ中可以看出，与普通网格化方法相似，可以自由选择η（对应环数）和ε（对应方位
数）的间隔，而在图１ｄ也会得到相应的空间坐标，从而可实现任意方位㊁任意环数的近区地改值计算，至于地形测量精度与网格参数选取之间的关系，则不在本文研究范围之内
㊂
ａ ＧＴＣＳ⁃１型近区地改仪原始数据采集高差特征点位；ｂ ａ中各数据点转化为计算坐标系点位；ｃ ｂ中需要网格化出的目标点位，调用式（６），采用１６方位㊁２０环进行计算，网格化出其高差；ｄ ｃ中各点反算后得到的直角坐标系点位
图１㊀曲网格克里金网格化方法（ＣＧＫＧ⁃Ｍ）在近区地改（０ ２０ｍ）中的应用示意
１．３㊀周期与对称扩边技术
扩边技术是网格化中重要的一项工作，常用的扩边技术有常数扩边㊁余弦扩边等㊂文中根据测量数据具有周期性和关于圆心对称的特征，采用周期扩边与对称扩边技术相结合，实现了对原始数据的实数据扩边，从而大大提高了网格化精度㊂图２为对称扩边与周期扩边技术示意图，Ａ区为实测数据区，模拟ＧＴＣＳ⁃１型近区地改仪野外数据采集点位分布，按１６方位㊁单方位２０点测量，每个方位控制２２．５ʎ，测量距离大于２０ｍ，转换到计算域中（Ａ区所示），由Ｃ可知，ＧＴＣＳ⁃１型近区地改仪的采集数据是以测点为圆心的圆内，而圆具有周期性和圆点对称性的特点，从而可以使用原始测量的数据对网格化区域进行扩边㊂如图Ｂ所示，根据周期性特征，第１方位的左边可使用第
１６㊁１５㊁ 等方位进行扩边，而第１６方位的右边可使用第１㊁２㊁ 等方位进行扩边；图２ｃ给出圆心对称扩边结果，如第１方位与第９方位关于圆心对称，所以第１和９方位可互为中心对称扩边数据，以此类推，可实现各方位关于圆点对称的数据扩边㊂
完成上述扩边后，对Ａ区上部的扩边可自由选择常数据补０或余弦扩边等扩边方法，这里不多做阐述㊂
２㊀模型测试
２．１㊀测试模型一：圆锥体模型
选用圆锥体模型对本方法进行测试，锥体斜面倾角３０ʎ，水平半径３０ｍ，采用１６方位㊁１０环对理论数据抽稀，以模拟该地形下的ＧＴＣＳ⁃１型近区地改仪数据采集情况㊂再应用曲网格克里金网格化方法（ＣＧＫＧ⁃Ｍ）对抽稀数据进行地形重构，并与理论模型进行比较，结果如图３㊁图４所示㊂图３给出了理论模型与ＣＧＫＧ⁃Ｍ重构地形的高差比
㊃
０５８㊃
㊀４期唐小平等：
曲网格克里金网格化方法及其在重力测量近区地形改正中的应用注：１ １６数值为数据测量方位；Ａ 网格数据区；Ｂ㊁Ｄ 周期扩边区域；Ｃ 中心对称扩边区
图２㊀
周期与对称扩边技术示意
注：等值线为理论模型高差等值线；黑色圆点为ＣＧＫＧ⁃Ｍ网出的高差点，呈圆环状分布
图３㊀凹陷模型等高线及其曲网格克里金网格化技术重构
较，从图中可以看出ＣＧＫＧ⁃Ｍ重构地形与理论模型差距很小，几乎完全吻合㊂为了进一步说明二者之间的差异，图４给出了ＣＧＫＧ⁃Ｍ重构地形与理论模型之间的剖面比较，图４ａ为二者的高差，圆点代表理论模型， ＋ 字地表ＣＧＫＧ⁃Ｍ重构地形，二者几乎完全重合，将该剖面的理论地形高与
ＣＧＫＧ⁃Ｍ重构地形高相减，从而得出二者的实际差距δｈ＝Ｈ－ｈ，图４ｂ给出了δｈ的分布情况，从图
４ｂ中可以看出δｈ很小，在０值附近成直线分布，均方根误差为０．００００４６ｍ，说明了本方法精度较高㊂
㊃
１５８㊃
物㊀探㊀与㊀化㊀探３９卷㊀
２．２㊀测点数据测试
随机选择野外实测模型进行数据处理分析，以测试本方法的实际应用情况㊂选择一处中间为
沟，三面为台地的模型作为本次试验的理论模型，采用１６方位㊁５ ８点测量模式，
先将测得数据应
ａ ＣＧＫＧ⁃Ｍ网出的高差；ｂ ＣＧＫＧ⁃Ｍ计算值与理论值之间的差异
图４㊀
曲网格克里金网格化法高差点重构与理论值之间的比较
注：等值线为理论模型高差等值线；黑色圆点为ＣＧＫＧ⁃Ｍ网出的高差点，呈圆环状分布
图５㊀复杂模型等高线及其曲网格克里金网格化技术重构
用Ｓｕｒｆｅｒ进行高差网格化处理，选用参数如下：克里金法为网格化方法，ｘ㊁ｙ坐标均为－２０ ２０ｍ，网格间距１ｍ㊂将网格化成果画出等值线，间隔
０．２５ｍ，再采用本文所述ＣＧＫＧ⁃Ｍ方法进行网格化计算，结果如图５中黑色圆点所示㊂从图中可以看出，黑色圆点所注明的高差与等高线大多能对应，只有较远的地方和个别地方有所差异，这是由于：网格化方位与测量方位不完全重合，以及远区数据量较少的缘故，但总体上讲这两种方法地形重构效果相当，说明了本方法精度可与常用的Ｓｕｒｆｅｒ网格化方法媲美㊂
３㊀误差分析
前面２．１㊁２．２节给出了两种不同模型的数据测试结果，进行了关于高差的相关误差分析，所以这里不再作有关高差的误差分析了，主要进行地改值的误差分析㊂选用圆锥体模型，选择锥体倾角变化范围为０ʎ ７５ʎ，水平半径为２０ｍ，分别采用锥体理论地改值计算与基于ＣＧＫＧ⁃Ｍ方法的圆域地改值计算［１－２］：
㊃
２５８㊃
㊀４期唐小平等：曲网格克里金网格化方法及其在重力测量近区地形改正中的应用锥体理论公式
Δｇ＝
２πＧρＲ
ｎ
（１－ｃｏｓｉ）㊀㊀扇形公式
ΔｇＴ＝
２πＧρＲ
ｎ
（Ｒｍ＋１－Ｒｍ＋Ｒ２ｍ
＋Δｈ２－Ｒ２ｍ＋１
＋Δｈ２）单位均为１０－８ｍ／ｓ２，计算结果如图６所示㊂
从图６ａ中可以看出ＣＧＫＧ⁃Ｍ计算值与理论
值地改值几乎重合，只有当锥体倾角变得很大时稍稍产生差异，而这个差异与理论地改值之间的相对误差，如图６ｂ所示，随着锥体倾角的增加，下
降趋势明显，当倾角大于５０ʎ又有所上升，最终收于０．２％左右㊂这说明了本文所采用的基于ＣＧＫＧ⁃Ｍ的近区地改值计算方法精度较高，稳定性较强，另外，ＧＴＣＳ⁃１型近区地改仪最大测量仰角为７０ʎ左右，所以本方法完全适应该仪器的数据处理范围，无需担心发散现象㊂
图６ｃ㊁６ｄ为文献［１－２］所示传统平面克里金网格方法计算的地改值及其与理论值之间的相对误差分析，从图中可以看出传统平面克里金网格方法随着坡度增大
，地改值与相对误差均逐渐增大，这是由于随着地形角度增加，平面网格的四角对地改值的影响逐渐增大，该方法稳定性较差㊂
ａ ＣＧＫＧ⁃Ｍ近区地形改正计算值与真值对比；ｂ ＣＧＫＧ⁃Ｍ计算值与近区地改真值之间的相对误差；ｃ 地改真值与传统圆域算法下的近区地改值对比；ｄ 传统圆域近区地改与真值之间的误差分析
图６㊀曲网格克里金网格化法（ＣＫＧ⁃Ｍ）与传统方法近区地改计算值及其与理论地改值之间的误差分析
４㊀结论
文中首先详细介绍了曲网格克里金网格化方法（ＣＧＫＧ⁃Ｍ）的原理与推导过程和ＧＴＣＳ⁃１近区地改仪的数据特点，再系统介绍了该方法的周期扩边与对称扩边，进行有关模型的地形重构测试和高差精度对比及其与ｓｕｒｆｅｒ进行实际数据的网格化地形重构比较，最后详细介绍了基于本方法
的圆域近区地改值计算，并进行了误差分析，总结出该方法具有精度高㊁稳定性强㊁使用方便㊁可与ＧＴＣＳ⁃１型近区地改仪配合使用等特点，具有广泛的应用前景㊂除了在近区地改中被应用外，该方法还可在中㊁远区地改及其他相关圆域网格化计算中使用㊂同时，对开发变网格圆域网格化方法具有重要意义㊂
㊃
３５８㊃
物㊀探㊀与㊀化㊀探３９卷㊀
参考文献：
［１］㊀ＤＺ０００４⁃９１重力调查技术规范（１ʒ５００００）［Ｓ］．１９９１，中华人民共和国国土资源部．
［２］㊀ＤＺ／Ｔ０１７１⁃１９９７大比列尺重力勘察规范［Ｓ］．１９９７，中华人民共和国国土资源部．
［３］㊀ＨａｍｍｅｒＳ．Ｔｅｒａｉｎｃｏｒｒｅｃｔｉｏｎｓｆｏｒｇｒａｖｉｍｅｔｅｒｓｔａｔｉｏｎｓ［Ｊ］．Ｇｅｏｐｈｙｓ⁃ｉｃｓ，１９３９，４：１８４－１９４．
［４］㊀ＮａｇｙＤ．Ｔｈｅｇｒａｖｉｔａｔｉｏｎａｌａｔｔｒａｃｔｉｏｎｏｆａｒｉｇｈｔｒｅｃｔａｎｇｕｌａｒｐｒｉｓｍ［Ｊ］．Ｇｅｏｐｈｙｓｉｃｓ，１９６６，３１：３６２－３７１．
［５］㊀ＦｏｒｂｅｒｇＲ．ＧｒａｖｉｔｙｆｉｅｌｄｔｅｒｒａｉｎｅｆｆｅｃｔｃｏｍｐｕｔａｔｉｏｎｓｂｙＦＦＴ［Ｊ］．ＢｕｌｌＧｅｏｄ，１９８５，５９：３４２－３６０．
［６］㊀ＬｉＹＣ，ＳｉｄｅｒｉｓＭ．Ｉｍｐｒｏｖｅｄｇｒａｖｉｍｅｔｒｉｃｔｅｒｒａｉｎｃｏｒｒｅｃｔｉｏｎ［Ｊ］．Ｇｅ⁃ｏｐｈｙｓｉｃａｌＪｏｕｒｎａｌｏｆＩｎｔｅｒｎａｔｉｏｎａｌ，１９９４，１１９：７４０－７５２．［７］㊀冯治汉．区域重力调查中的中区地形改正方法及精度［Ｊ］．物探与化探，２００７，３１（５）：４５５－４５８．
［８］㊀张俊，张宝松，邸兵叶，等．高程数据网格间距对重力中区地形改正精度的影响［Ｊ］．物探与化探，２０１４，３８（１）：１５７－１６１．［９］㊀ＳｃｈｉａｖｏｎｅＤ，ＣａｐｏｌｏｎｇｏＤ，ＬｏｄｄｏＭ．Ｈｉｇｈｒｅｓｏｌｕｔｉｏｎｄｅｍｓｆｏｒｎｅａｒｓｔａｔｉｏｎｔｅｒｒａｉｎｃｏｒｒｅｃｔｉｏｎｉｎｇｒａｖｉｍｅｔｅｒｙ［Ｒ］．ＥＧＭ２００７Ｉｎ⁃ｔｅｒｎａｔｉｏｎａｌｗｏｒｋｓｈｏｐ，２００７．
［１０］ＳｃｈｉａｖｏｎｅＤ，ＣａｐｏｌｏｎｇｏＤ，ＬｏｄｄｏＭ．Ｎｅａｒ⁃ｓｔａｔｉｏｎｔｏｐｏｇｒａｐｈｉｃｍａｓｅｓｃｏｒｒｅｃｔｉｏｎｆｏｒｈｉｇｈ⁃ａｃｃｕｒａｃｙｇｒａｖｉｍｅｔｒｉｃ［Ｊ］．ＧｅｏｐｈｙｓｉｃｓＰｒｏｄｐｅｃｔｉｎｇ，２００９，５７（５）：７３９－７５１．［１１］刘宽厚，耿涛，杨怀英，等．基于便携激光测距仪的重力测量近区地形改正系统［Ｊ］．物探与化探，２０１２，３６（３）：４０３－４０８．［１２］ＴｈｏｍｐｓｏｎＪＦ，ＷａｒｓｉＺＵＡ，ＭａｓｔｉｎＣＷ．ＮｕｍｅｒｉｃａｌＧｒｉｄＧｅｎｅｒａ⁃ｔｉｏｎ：ＦｏｕｎｄａｔｉｏｎｓａｎｄＡｐｐｌｉｃａｔｉｏｎｓ［Ｐ］．ＮｏｒｔｈＨｏｌｌａｎｄ，ＩＳＢＮ：０４４４００９８５Ｘ，１９８５．
［１３］ＦｏｍｂｅｒｇＢ．Ｔｈｅｐｓｅｕｄｏｓｐｅｃｔｒａｌｍｅｔｈｏｄ：ａｃｃｕｒａｔｅｒｅｐｒｅｓｅｎｔａｔｉｏｎｏｆｉｎｔｅｒｆａｃｅｓｉｎｅｌａｓｔｉｃｗａｖｅｃａｌｃｕｌａｔｉｏｎｓ［Ｊ］．Ｇｅｏｐｈｙｓｉｃｓ，１９８８，５３：６２５－６３７．
［１４］ＮｉｅｌｓｅｎＰ，ＦｌｅｍｍｉｎｇＩ，ＢｅｒｇＰ，ｅｔａｌ．Ｕｓｉｎｇｔｈｅｐｓｅｕｄｏ⁃ｓｐｅｃｔｒａｌｍｅｔｈｏｄｏｎｃｕｒｖｅｄｇｒｉｄｓｆｏｒ２Ｄｅｌａｓｔｉｃｆｏｒｗａｒｄｍｏｄｅｌｉｎｇ［Ｊ］．Ｇｅｏ⁃ｐｈｙｓｉｃａｌＰｒｏｓｐｅｃｔｉｎｇ，１９９５，４３：３６９－３９５．
［１５］蒋丽丽，孙建国．基于Ｐｏｉｓｓｏｎ方程的曲网格生成技术［Ｊ］．世界地质，２００８，２７（３）：２９８－３０５．
［１６］孙建国，蒋丽丽．用于起伏地表条件下的地球物理场数值模拟的正交曲网格生成技术［Ｊ］．石油地球物理勘探，２００９，４４（４）：４９４－５００．
［１７］郭良辉，孟小红，郭志宏，等．地球物理不规则分布数据的空间网格化法［Ｊ］．物探与化探，２００５，２９（５）：４３８－４４２．［１８］ＭａｔｈｅｒｏｎＧ．Ｐｒｉｎｃｉｐｌｅｓｏｆｇｅｏｓｔａｔｉｓｔｉｃｓ［Ｊ］．ＥｃｏｎｏｍｉｃＧｅｏｌｏｇｙ，１９６３，５８：１２４６－１２６６．
［１９］ＭａｔｈｅｒｏｎＧ．Ｔｈｅｉｎｔｒｉｎｓｉｃｒａｎｄｏｍｆｕｎｃｔｉｏｎｓ，ａｎｄｔｈｅｉｒａｐｐｌｉｃａｔｉｏｎｓ［Ｊ］．Ａｄｖ．Ａｐｐｌ．Ｐｒｏｂ．，１９７３，５：４３９－４６８．
Ｃｕｒｖｅｄ⁃ｇｒｉｄｋｒｉｇｉｎｇｇｒｉｄｄｉｎｇｍｅｔｈｏｄａｎｄｉｔｓａｐｐｌｉｃａｔｉｏｎｔｏｎｅａｒ⁃ｓｔａｔｉｏｎｔｅｒｒａｉｎｃｏｒｒｅｃｔｉｏｎ
ｉｎｇｒａｖｉｔｙｍｅａｓｕｒｅｍｅｎｔ
ＴＡＮＧＸｉａｏ⁃Ｐｉｎｇ１，２，ＬＩＵＫｕａｎ⁃Ｈｏｕ１，ＧＥＮＧＴａｏ１，ＦＥＮＧＺｈｉ⁃Ｈａｎ１
（１．Ｘｉ＇ａｎＣｅｎｔｅｒｏｆＧｅｏｌｏｇｙＳｕｒｖｅｙ，ＣｈｉｎａＧｅｏｌｏｇｙＳｕｒｖｅｙ，Ｘｉ＇ａｎ㊀７１００５４，Ｃｈｉｎａ；２．ＫｅｙｌａｂｏｒａｔｏｒｙｆｏｒｔｈｅＳｔｕｄｙｏｆＦｏｃｕｓｅｄＭａｇｍａｔｉｓｍａｎｄＧｉａｎｔｏｒｅＤｅｐｏｓｉｔｓ，ＭＬＲ，Ｘｉ＇ａｎ㊀７１００７４，Ｃｈｉｎａ）
Ａｂｓｔｒａｃｔ：Ｉｎｔｈｉｓｐａｐｅｒ，ａｈｉｇｈ⁃ａｃｃｕｒａｃｙａｎｄｓｔｒｏｎｇｓｔａｂｉｌｉｔｙｎｅａｒ⁃ｚｏｎｅｔｅｒｒａｉｎｃｏｒｒｅｃｔｉｏｎｍｅｔｈｏｄｉｓｄｉｓｃｕｓｓｅｄ．Ｆｉｒｓｔｌｙ，ｃｏｍｂｉｎｉｎｇＫｒｉｇ⁃ｉｎｇｇｒｉｄｄｉｎｇｍｅｔｈｏｄｗｉｔｈｃｕｒｖｅｄ⁃ｇｒｉｄｃｏｍｐｕｔａｔｉｏｎｔｅｃｈｎｏｌｏｇｙ，ａＣｕｒｖｅｄ⁃ＧｒｉｄＫｒｉｇｉｎｇＧｒｉｄｄｉｎｇＭｅｔｈｏｄ（ｓｈｏｒｔｆｏｒＣＧＫＧ⁃Ｍ）ａｎｄｉｔｓｅｘ⁃ｔｅｎｄｉｎｇｅｄｇｅｔｅｃｈｎｉｑｕｅｓａｒｅｉｍｐｌｅｍｅｎｔｅｄｔｏｒｅａｌｉｚｅｄａｔａ⁃ｇｒｉｄｄｉｎｇｉｎｃｉｒｃｌｅａｒｅａ；Ｎｅｘｔ，ｔｈｉｓｍｅｔｈｏｄｉｓｕｓｅｄｔｏｃｏｍｐｕｔｅｎｅａｒ⁃ｚｏｎｅｔｅｒｒａｉｎｃｏｒｒｅｃｔｉｏｎｖａｌｕｅａｎｄｒｅａｌｉｚｅｔｈｅｇｒａｖｉｔｙｍｅａｓｕｒｅｍｅｎｔｎｅａｒ⁃ｚｏｎｅｔｅｒｒａｉｎｃｏｒｒｅｃｔｉｏｎｖａｌｕｅ＇ｓｃｏｍｐｕｔａｔｉｏｎｗｈｉｃｈｃｏｎｔａｉｎｓａｎｙａｎｇｌｅｓａｎｄｃｉｒｃｌｅｓ．Ｆｉｎａｌｌｙ，ｓｏｍｅｍｏｄｅｌｉｎｇｔｅｓｔｓ，ｒｅａｌｄａｔａｃｏｍｐｕｔａｔｉｏｎａｎｄｅｒｒｏｒａｎａｌｙｓｉｓｏｆｔｈｉｓｍｅｔｈｏｄａｒｅｃａｒｒｉｅｄｏｕｔ．ＴｈｅｒｅｓｅａｒｃｈｒｅｓｕｌｔｈａｓｓｈｏｗｎｔｈａｔＣＧＫＧ⁃Ｍｎｅａｒ⁃ｚｏｎｅｔｅｒｒａｉｎｃｏｒｒｅｃｔｉｏｎｖａｌｕｅｃｏｍｐｕｔａｔｉｏｎｍｅｔｈｏｄｈａｓｓｏｍｅｃｈａｒａｃｔｅｒｉｓｔｉｃｓ，ｓｕｃｈａｓｈｉｇｈｃｏｍｐｕｔａ⁃ｔｉｏｎａｃｃｕｒａｃｙ，ｓｔｒｏｎｇｓｔａｂｉｌｉｔｙａｎｄｆｌｅｘｉｂｉｌｉｔｙ．Ｉｔｉｓａｓｉｍｐｌｅ，ｅｆｆｉｃｉｅｎｔａｎｄｆａｓｔｍｅｔｈｏｄｆｏｒｎｅａｒ⁃ｚｏｎｅｔｅｒｒａｉｎｃｏｒｒｅｃｔｉｏｎｖａｌｕｅｃｏｍｐｕｔａｔｉｏｎｗｈｉｃｈｍａｔｃｈｅｓｔｈｅｍｅａｓｕｒｅｍｅｎｔｄａｔａｆｏｒｍｏｆＧＴＣＳ⁃１ｇｒａｖｉｔｙｍｅａｓｕｒｅｍｅｎｔｎｅａｒ⁃ｓｔａｔｉｏｎｃｏｒｒｅｃｔｉｏｎｉｎｓｔｒｕｍｅｎｔ．
Ｋｅｙｗｏｒｄｓ：ｃｕｒｖｅｄ⁃ｇｒｉｄｋｒｉｇｉｎｇｇｒｉｄｄｉｎｇｍｅｔｈｏｄ；ｃｉｒｃｌｅｒｅｇｉｏｎ；ｎｅａｒｓｔａｔｉｏｎｔｅｒｒａｉｎｃｏｒｒｅｃｔｉｏｎ
作者简介：唐小平（１９８２－），男，硕士，助理研究员，主要从事矿产综合地球物理研究工作㊂Ｅ⁃ｍａｉｌ：ｘｉａｏｐｉｎｇｔａｎｇ＠１２６．ｃｏｍ㊃４５８㊃。