张 勇, 梁冠军, 黄展华

(滁州职业技术学院 土木工程系,安徽 滁州 239000)

自2008年7月1日开始，我国正式采用2000国家大地坐标系，且此后的各类测绘地理信息成果一律采用2000国家大地坐标系.而我们国家的CGCS2000框架目前只能提供相应于ITRF97框架、2000历元的静态坐标.因此不利于科研和工程中应用最为普及的GPS精确定位(ITRF2005框架和当前历元)，有进行转换的必要性.故须以国家能获得的高分辨率速度场资料建立我国地壳运动速度场模型，并在此基础上进行框架转换及历元归算，方可将测量成果纳入我国2000大地坐标系.

基于上述原因，通过插值或数值模拟方法求得具有足够点位密度且分布均匀的速度场已经成为迫切需要解决的问题.国内外不少专家学者已经采用多种方法建立了我国大陆区域地壳运动速度场.刘经南等采用多面函数法对中国区域水平运动进行了数字模拟[1]，杨少敏等用双三次样条函数建立了中国大陆构造形变场[2]，蒋志浩采用有限元插值方法(三角形三次多项式插值算法)建立了中国地壳运动速度场模型[3].从数学的角度发现插值点与已知点的关系是这些模型的优点，但是没有充分考虑板块固有的物理及力学意义，因此是有缺陷的.

20世纪中叶地质学家D.G.Krig首先提出了Kriging插值方法，该方法用变异函数为工具，是一种线形无偏且计算准确的空间内插方法，得到了越来越广泛的应用.大量实践证明，该插值法建立起来的模型能够很好……