应用拟合推估法建立我国大陆地壳水平运动速度场的研究

曾安敏，刘光明，秦显平，唐颖哲

(1．信息工程大学地理空间信息学院，河南郑州450052;2．地理信息工程国家重点实验室，陕西西安710054;3．西安测绘研究所，陕西 西安710054)

一、引 言

中国大陆除遭受印度板块的强烈碰撞以及菲律宾板块、太平洋板块的俯冲影响外，板块内形变也非常复杂和不均匀[1-3]。在具有如此复杂的地壳运动背景的大陆构建大地测量坐标框架，必须顾及板块的整体运动和局部运动。我国已建立了2000中国坐标系CGCS2000[4-6]，2008年8月正式启用相应 的 参 考 框 架 CTRF2000[6-7]。 但 由 于 建 立CTRF2000时的资料截止到2002年，重复观测时间较短，绝大部分框架点未赋予点位速度场信息，无法反映我国坐标框架的动态变化，不能完整地表述我国坐标系统的现势性[8-9]。由于板块运动引起的坐标变化是不容忽视的，速度场数据成为现代大地测量不可或缺的基础性数据。地壳运动的描述往往采用欧拉矢量法和统计拟合的方法。基于实测速度信息可以采用最小二乘准则求得欧拉矢量的估值，常用的统计拟合方法有多面函数法、有限元插值法、拟合推估法等。近10年来，国内不少学者在2002年前利用“网络工程”数据陆续开展我国大陆板块运动相关研究工作[10-12]，但由于使用数据少，其成果的实用价值有一定的局限性。2009年前利用“网络工程”数据研究了最新的速度场，有学者利用局部欧拉矢量法建立了3度格网平均值[13]，但用于建立局部欧拉矢量的点及其分布是其关键。也有学者利用有限元法建立了速度模型[14]，但面元的确定是关键。……