川滇菱形块体应变率场解算及其地壳水平形变趋势分析.docVIP

川滇菱形块体应变率场解算及其地壳水平形变趋势分析 　　摘要：基于普通Kriging插值理论建立了川滇地区均匀格网节点上的速率值，然后结合弹性力学几何方程与Taylor级数展开解算每个图形单元的平均应变分量，从而最终得到了川滇地区地壳应变率场分布。结果显示，川滇块体总体处于近南北向的拉张状态，局部地壳压缩趋势明显。最大剪应变分布、单位体积地震能量累积速率分布和区内重大走滑断裂分布一致性良好，能够清晰的显示出青藏高原内部物质运移的趋势。文中应变率场解算结果显示青藏高原东缘物质的逃逸形式与区域乃至高原整体长期的地质活动特征之间具有继承关系。 　　关键词：Kriging理论，泰勒级数，川滇菱形块体，应变率，GPS速度场 　　 　　 Strain Field Calculation and Horizontal Crustal Deformation Analysis of Sichuan-Yunnan Block 　　 　　Abstract：Based on ordinary Kriging interpolation theory, we interpolated the scattered GPS velocity data of Sichuan-Yunnan block and its adjacent areas to grid point values by the Kriging method. We both calculated the strain components of each unit and the strain trends of the block based on Taylor series expansion and strain component expressions in elastic mechanics. The results show that the block is in a N-S expanded state but the local compression trend is obvious. Distributions of maximum shear strain, seismic energy and main slip fracture are very consistent. The result of the strain rate in this study suggests that the contemporary tectonic of the Sichuan-Yunnan block inherits the long term geological deformation. 　　Keyword: Kriging theory, Taylor series, Sichuan-Yunnan block, Strain rate, GPS velocity field; 　　中图分类号：TN823文献标识码： A 文章编号： 　　1 引言 　　随着GPS观测资料的不断积累，利用GPS观测数据计算地应变率的方法在地壳动力学的研究当中的应用越来越广泛。但是对于同一个研究区域，不同的学者基于几乎相同的观测资料，其解算结果往往差异明显。如何有效减小空间分布不均的离散数据对应变率解算带来的误差，目前还没有严密、完备的数学模型。利用GPS观测资料计算地应变率的本质是定义在点函数上的数值微分，加之GPS观测资料是在地球表面以离散点的形式表示。Savage等（2001）在球坐标系下推导了计??图形单元平均应变参数的表达式，但是这种方法在进行地应变率的解算时需要根据离散的观测数据建立一系列图形单元，通过解算每个图形单元内的应变参数，最终获得研究区域整体的应变场。但是实际GPS测站分布由于受到多种观测因素的制约，其空间分布往往极不规则，直接根据实测资料建立满足要求的子单元往往不太现实。那么如何在同时顾及应变解算精度及其分辨率的情况下，获取研究区域的应变率场呢？研究这普遍应用反距离加权的最小二乘法与最小二乘配置来建立均匀分布的形变场[1,2]。反距离加权的最小二乘法的解算结果往往受观测点密集程度的影响很大。姚道荣[3]等研究结果表明，当空间变量不满足二阶平稳假设时，最小二乘配置就难以应用。并且在相同的观测噪声下，数据分辨率对最小二乘配置估值的影响远大于Kriging插值结果。地质统计学中的Kriging插值方法是专门应对空间分布不均匀的离散数据插值问题的一种行之有效的方法。该方法可以根据空间数据的结构性与变异性，实现无偏、最优的插值。 　　川滇菱形块体位于青藏高原东缘，是高原物质东流的主要通道。由于来受到青藏高原物质东南方向的顺时针旋转挤出以及阿萨姆构