论文部分内容阅读
地壳运动及其形变的运动学及动力学研究一直是地球科学领域的热点和难点课题。近几十年来,以GPS为代表的高精度空间大地测量技术的崛起,开辟了研究地壳运动与板块构造变化的新途径,尤其是在地壳弹塑性形变研究方面取得较大进展,但对其数学物理机制的研究尚缺乏系统性,该领域的研究需进一步深入。本文针对地壳形变分析中传统地壳刚体旋转运动模型难以顾及板内形变的问题,根据地壳形变的物理机制,利用模型补偿思想,采用不同数学方法描述板内形变;针对现有弹性地壳运动形变分析模型对块体内部形变假设过强的问题,提出两种优化的弹性地壳形变分析模型。利用“中国大陆构造环境监测网络(CMONOC)”在环渤海地区的GPS观测速度场资料,研究了环渤海区域整体及各主要块体的地壳运动及形变规律,具体内容包括:1、讨论了大地测量观测数据与地质、地球物理资料建立地壳运动及其形变分析模型的原理,分析了研究地壳形变的必要性和科学意义;梳理了空间大地测量资料在地壳运动及形变分析领域的主要研究成果和存在的主要问题。2、阐述了大地测量及地质和地球物理两种技术手段进行地壳运动及形变建模分析的参考框架理论,并讨论分析了ITRF参考框架的建立和更新方法。3、探讨了地壳弹性运动模型的提出依据及建模理论,建立了喜马拉雅块体和祁连山块体的整体旋转与均匀应变模型(REHSM)和整体旋转与线性应变(RELSM)两种弹性运动模型。研究结果显示,两种模型结果与实际运动的差异均较小,明显优于刚体地壳运动模型,且对块体形变分析结果与地质结论符合较好。4、鉴于地壳整体旋转与均匀应变模型可以表达块体内部主要变形,本文提出利用“非参数”描述偏离整体旋转与均匀应变模型的不规则形变,从而建立地壳形变的半参数分析模型。对祁连山块体的GPS速度场数据及环渤海区域“陆态网”数据进行拟合分析,结果表明,新模型有效的改善了整体旋转与均匀应变模型。由于新模型并不假设块体内部形变是均匀变化的,也不是线性变化的,实际上是优化了两种弹性运动形变分析模型。5、将块体内部偏离整体旋转与均匀应变的不规则形变看作是对整体旋转与均匀应变模型所确定的主体运动的干扰信号,建立了地壳形变分析的最小二乘配置模型,考虑到地壳形变均有明显的区域性和方向性差异特点,建模时采用分区分方向确定协方差函数相关参数。分析结果表明,最小二乘配置模型的分析结果优于块体的整体旋转与均匀应变模型,表明本文方法进一步精化了整体旋转与均匀应变模型。6、针对GPS速度场数据可能包含粗差影响地壳形变分析结果的问题,提出一种利用小波变换剔除速度场噪声从而优化第5条中地壳形变分析的最小二乘配置模型。对中国大陆6个次级板块的速度场数据进行建模分析结果显示,采用小波变换对块体运动速率进行预处理,再进行相关参数的求解,得到的参数解更加符合实际,处理结果的精度要高于最小二乘配置方法。7、建立了环渤海区域及邻区的小波辅助的最小二乘配置法优化的整体旋转与均匀应变模型,绘制了环渤海区域地壳形变、应变的空间分布图像;所得应变参数表明:环渤海区域整体上呈现0.24×10-9/a的NW-SE22.1°的双向趋势性扩张运动,主压应变轴方向基本为NEE-SWW,主张应变轴基本为NNW-SSE,且主压应变轴方向为NE47.70°89.74°,这与地球物理方法得到该区的主压应变轴的优势方向是大体一致的,表明环渤海区域现今地壳运动是相对稳定的。