gps数据在庐断裂带中二段运动及变形分析中的应用.docxVIP

gps数据在庐断裂带中二段运动及变形分析中的应用 桐带及围区构造分析 断裂带是东亚大陆上一些北北东向的重大扰动系的主干断裂带。它通常起北起黑龙江省，南起长江。它通常向北东移动，穿过中国大陆东部。 郯庐断裂带是由多组呈斜列分布的多条断裂所组成的断裂带, 龚伟等根据各段的地质结构特征, 将该断裂带北起黑龙江边的肇兴县, 南达长江北岸的广济县, 分为三段, 即: 北段(肇兴—昌图), 中段(昌图—嘉山), 南段 (嘉山—广济)。 郯庐断裂带中南段在第四纪以来活断层最发育, 历史上曾发生过多次中强以上地震, 是我国目前主要地震活动区域之一。 2008年5月12日汶川8.0级地震发生后, 地下流体多测项震后效应在郯庐断裂带中南段部分地区一直持续存在, 同时速度结构研究结果表明该区域存在发生中强以上地震的介质速度结构条件, 因此加强对郯庐带现今地壳变形特征的认识, 对郯庐带未来地震危险性的判定有着重要的意义。 利用GPS观测资料研究地壳形变特征主要包括对GPS资料进行分析和对位移场/速度场及其衍生场进行分析。 通过对速度场及其衍生场进行分析, 研究块体的运动特征、 块体的变形特征、 大范围低频应变场分布特征等已经广泛被采用。 而利用GPS观测资料直接分析断裂带的运动及变形特征却为数不多, 特别是对华北地区主要断裂带变形特征的研究。 本文在前人研究的基础上, 首先基于块体整体运动及均匀应变模型, 对中国地壳运动观测网络GPS区域网1999—2001年、 2001—2004年、 2004—2007年及2007—2009年四期速度场数据进行处理, 分析了郯庐带的运动及变形的动态演化特征, 分析汶川地震前后郯庐带变形特征的差异特性; 并通过计算1999—2004年, 2004—2007年及2007—2009年三期扣除渤海湾区域整体刚性运动的GPS速度场, 分析了环渤海区域周边地块间的相对运动状态, 进一步认识整个环渤海区域的构造运动背景。 1 gps数据的选择和数据处理方法 1.1 gps速度场协方差建模及模型参数 刘晓霞等利用GPS观测资料对中国大陆常用的三个描述块体运动及变形的模型, 即李延兴等提出的块体整体旋转与线性应变模型(RELSM), 块体整体旋转与均匀应变模型(REHSM)和刚体旋转模型(RRM)进行了分析。 文中采用拟合残差的离散程度作为判定指标, 结果显示在华北地区REHSM的拟合残差的离散度最小, 为三种模型中的最优模型, 能较真实的反映华北地区的运动及变形方式。 本文选用中国地壳运动观测网络GPS区域网1999—2001年、 2001—2004年、 2004—2007年及2007—2009年四期相对于欧亚大陆的GPS速度场结果, 首先利用这四期的资料, 基于块体运动与均匀应变模型(REHSM)对郯庐带中南段的运动及变形特征进行研究, 而后采用1999—2004年、 2004—2007年、 2007—2009年三期GPS速度场, 扣除环渤海区域的整体刚性运动, 研究GPS速度场所揭示的环渤海地区的内部变形特征。 因华北地区非构造因素噪声较显著, 所以首先利用最小二乘配置方法对四期GPS速度场进行滤波, 剔除非构造因素噪声。 江在森等对最小二乘配置滤波及速度场协方差矩阵的构建理论进行了详细的讨论, 并给出中国大陆高斯型函数模型参数K=0.001658。 鉴于本文研究的区域为华北地区, 比中国大陆区域小, 所以利用华北地区1999—2004年期相对于欧亚大陆的GPS速度场(图1, 一倍误差椭圆)对高斯型函数模型参数K进行了重新计算。 图1给出了中国大陆地壳运动观测网络1999—2004年GPS速度场(华北地区346个站), 将华北地区实测的扣除整体刚性运动的速度作为信号计算得到的协方差与高斯函数拟合的结果如图2所示。 图2显示1999—2004年期的高斯型函数参数C=0.659,K=0.002212, 本文即根据此高斯型函数参数建立最小二乘配置滤波的协方差阵, 分别对1999—2001年、 2001—2004年、 2004—2007年、 2007—2009年四期GPS速度场进行滤波, 然后将滤波后的GPS速度场资料作为REHSM模型的输入, 求解各块体的模型参数, 然后沿边界带取9个代表点, 基于REHSM反算这9个边界带的速度, 最后取两相邻地块模拟的速度的差矢量, 进而分析边界带的运动及变形状况。 表1为四期的华北平原地块、 鲁东黄海地块的模型参数及精度情况。 虽然华北地区的变形量很小, 但表1显示基于REHSM模型解算的欧拉参数及应变参数依旧比解算的参数误差高2个数量级, 由此可见REHSM模型的模型参数是可信的, 由REHSM模型计算的边界点的速度在一定程度上能反映边界带运动的真实情况。 郯庐带北段两侧块体内的GPS站较