四川省CORS网地壳运动速度场模型研究

随着空间测量技术的不断发展，连续运行观测站(CORS)的建立不仅提高了相对定位的精度，也为地球物理研究提供了便利条件。GPS具有高精度、全天候和高机动性等特点，使通过获得大范围地壳运动时间序列研究某一地区速度场成为可能[1]。四川盆地位于构造活动比较活跃的几个断裂带上，如龙门山断裂带、鲜水河断裂带、安宁河-则木河断裂带等[2]，近年来地震频发。本文利用GAMIT/GLOBK软件解算2006～2015年的观测数据，得到四川省CORS站在ITRF2008框架下的水平速度场和垂直速度场，并将其转换到欧亚框架下，为揭示四川省的地壳运动规律提供依据。

1 CORS网数据处理

本文选择覆盖四川地区的50个CORS站作为研究速度场的基准站。为了提高数据解算的精度，并入13个均匀分布在我国及周边地区的IGS站组成区域网[3]。先用GAMIT进行基线解算得到单日松弛解，然后用GLOBK进行滤波平差，最终获得四川省CORS站的速度场模型。

陆游常常自觉不自觉地以一个越人之眼，比较着、衡量着梁益地区的各个方面，而他的比较与衡量中总带着一些政治、经济、文化各方面的优越感。这种与生俱来的优越感，让陆游在远离政治文化中心的梁益地区生活以及书写时，时常具有一种有意无意地居高临下式的审视与判断，将原本可以更加客观的比较衡量，变成带有明显偏爱色彩的俯视性考量。

1.1 IGS站选取

引入IGS站可以对整个四川CORS网起到控制作用[4]。本文选择我国11个IGS站(BJFS、CHAN、SHAO、WUHN、HKWS、KUNM、LHAZ、URUM、GUAO、TWIF、NCKU)以及周边2个站(SUWN、ULAB)。

1.2 获取速度场方案

主要考虑以下几个因素：电离层的折射影响用LC观测值；截止高度角为10°；海潮文件改正采用GAMIT提供的模型。速度场获取步骤如下：

由表2可知，无论是N方向或E方向，还是整体水平方向，位于断裂带或者断裂带周边CORS站的运动速率都高于平原地区非断裂带CORS站的速率，但在垂直方向上不能确定断裂带上CORS站的运动速率就一定大于非断裂带上的运动速率。

人民满意的教育，只有开头，没有结尾；广州好教育，只有起点，没有终点。在屈哨兵局长的带领下，一群对教育用心、用情、用力的教育人正在路上，也永远在路上；对好教育的探索与追寻正在路上，也永远在路上。我们相信，有他们的努力，古老而又年轻的广州将迎来教育大发展、民生更幸福的美好明天！

由GPS连续观测站获得速度场，一般都是根据连续观测站得到的坐标时间序列，估计由于同震位移或者仪器更换产生的阶跃以及可能存在的周期性运动，从而构建模型求得速度场。从时间序列角度，两点就能确定一条直线，即2 d的数据就能求得一个速度。经计算认为3 a是求稳定速率的最低时间长度要求，且数据连续性要好[6]，因此本文选用10 a的数据长度得到的速度是可靠的。表1给出了ITRF2008框架下四川省CORS站的三维速度估值及中误差统计。

图2给出四川省CORS站水平方向速度场模型。结合表1看出，四川省CORS站在东西方向的最大运动速率为50.71 mm/a，最小为31.6 mm/a，平均为37.95 mm/a；南北方向的最大运动速率为-0.9 mm/a，最小为-19.62 mm/a，平均为-11.05 mm/a；水平方向平均运动速率为38.72 mm/a，优势方向为S73.9°E，与青藏高原物质整体沿川滇活动地块东南向运动这一结论相符[7]。

2 四川省CORS站速度场分析

图1 CHDU站坐标时间序列图 Fig.1 Coordinate time series of CHDU station

2)由于GLOBK的主要原理是卡尔曼滤波，将单时段获得的h-file作为GLOBK的输入并获得测站的原始坐标时间序列，通过分析找到异常点予以剔除，获得干净的时间序列[5]。

表1 ITRF2008框架下四川CORS的三维速度估值及中误差统计/mm·a-1 Tab.1 Three dimensional velocity estimation and error estimation of Sichuan CORS under ITRF2008 framework

VEVNVUSESNSU最大值50．710-0．90038．6200．8481．1981．438最小值31．600-19．620-9．5700．4530．5300．470平均值37．955-11．0513．4290．6110．6940．653

NFC手机支付是一种把手机变成虚拟卡的卡模拟技术，是一种近场高频无线通信技术，允许设备间点对点非接触式数据交互，由非接射频RFID技术演变而来。NFC设备靠近非接读卡器后NFC设备进入工作状态，NFC设备与读卡器根据APDU协议建立交互，读卡器将敏感数据发给NFC设备的安全载体进行处理，NFC设备将响应数据发给读卡器。

由于四川盆地位于几个地质构造活动比较活跃的断裂带上，为了对比断裂周围的CORS站运动速度和平原非断裂带上CORS站运动速率大小，给出龙门山断裂带、鲜水河断裂带、安宁河-则木河以及平原地带的CORS站结果，如表2所示。

表2 平原非断裂带与断裂带CORS站三维速度估值/mm·a-1 Tab.2 Three dimensional velocity estimation of nonfault zone and fault zone CORS station in plain

地貌VEVNVUVH成都平原35．653-9．1402．90536．806龙门山36．873-9．4822．73538．073鲜水河41．858-12．7088．16343．745安宁河-则木河38．527-16．362-0．91756．889

1)用GAMIT进行处理，得到单日松驰解，检查基线重复率和post_nrms是否符合精度要求。

2.1 水平速度场分析

3)把所有单日解得到的h-file合并成一个单独的H-file，然后运行glred/glorg得到所有CORS站新的坐标时间序列，然后运行scgps程序获得速度场。图1给出CHDU站的坐标时间序列。

由表1可知，CORS网的速度场解算精度都比较高，N、E、U方向上的均方根误差均值优于0.7 mm，水平方向整体运动速度较垂直方向要快，说明四川省的地壳运动在水平方向运动较为活跃。

图2 四川CORS网ITRF2008框架下水平方向速度场模型 Fig.2 Horizontal direction velocity field model of Sichuan CORS network ITRF2008 framework

NNR-NUVEL-1A模型作为ITRF的国际标准模型，它是基于地质和地球物理数Ma的资料建立的[8]，如表3所示。通过NNR-NUVEL-1A模型将ITRF2008框架下水平方向速度转换到欧亚框架下[9]，表4给出了欧亚框架下的四川省CORS站的三维速度估值及中误差统计。

表3 NNR-NUVEL1-1A 模型 Tab.3 NNR-NUVEL1-1A model

ωx/rad·Ma-1ωy/rad·Ma-1ωz/rad·Ma-1经度/(°)纬度/(°)w/(°)·Ma-1-0．0009836-0．00239840．0031558-112．350．60．234

表4 欧亚框架下四川省CORS站的三维速度估值及中误差统计 Tab.4 Three dimensional velocity estimation and error estimation of CORS stations in Sichuan province under the framework of Eurasia

水平速度VEVNSESN最大值28．50126．501 7．0590．8481．482最小值7．8077．465-10．4870．4240．630平均值14．39813．758-1．4650．6091．001

对比表1和表4看出，在欧亚框架下四川省CORS站在东西向和南北向的测站速度均小于ITRF2008框架下的速度，尤其是东西向测站速度平均减小24 mm/a。欧亚框架下测站水平方向的平均速度为14.40 mm/a，较ITRF2008框架下的速度降低24.32 mm/a。图3为欧亚框架下四川省CORS站的水平方向速度场模型。

图3 四川CORS网欧亚框架下水平速度场模型 Fig.3 Horizontal velocity field model under the framework of Sichuan

由图3看出，所有测站水平方向的速度方向具有较好的一致性，优势方向为S77.78°E，在龙门山断裂带、鲜水河断裂带和安宁河-则木河断裂带上的测站水平方向运动速度平均值大于平原非断裂带测站。

2.2 垂直速度场分析

图4给出四川省CORS站垂直方向速度场模型。结合表1，相较于水平方向，所有测站的垂向构造运动速度相对比较微弱，平均为3.43 mm/a。通过图1给出的坐标时间序列也可以看出，在垂直方向非构造干扰比较显著，且具有明显的周期性变化。四川省CORS站的垂向速度在大部分区域表现为隆升趋势，有较少的一部分测站区域表现为下降状态。从整体上看，川西地区的测站表现为隆升状态而且北部CORS站上升的平均速度要大于其他地区，出现下沉现象的测站基本上位于东部平原地区。

关于四川盆地地壳隆升问题，国内学者普遍认为与青藏高原物质的南东向滑移受阻密切相关，但具体的滑移受阻机制仍然有着不同的观点。而对于地表下沉则可能与地下水过度流失以及煤矿、石油等自然资源的过度采伐有关。通过图4可以看出，地表下降的主要测站有ANYU、CAXI、GANL、JIGE、QLAI、XICH、YBIN、ZTON、CHDU、KAIY，主要位于东部平原地区，如成都、宜宾、邛崃、西昌，这些地方人口密度较大，且工农业发达，在城市建设中地面工程垂向破坏较大。例如，成都的下沉速度为0.23 mm/a，下沉原因主要包括：近年来随着地铁的修建以及房地产建设的迅速发展，使得垂向地表破坏较大；其次，在施工建设中需要抽取地下水，过度采伐也是地表下沉的原因之一；地质构造方面，成都平原处于挤压状态，成都西侧受逆冲作用的影响呈现抬升趋势，东侧呈现沉降趋势。

图4 四川CORS网ITRF2008框架下垂直方向速度场模型 Fig.4 Vertical velocity field model under the framework of ITRF2008 in Sichuan

3 应变场分析

速度场虽然能够直观地反映CORS站区域运动特性，但其与参考基准的选取密不可分，且速度场无法反映一些精细的形变特征。而应变场则与参考基准无关，且区域应变场直接反映局部应力应变的空间分布特性[9]。

一年一度的全国水利厅局长会议是水利工作的航塔。翻读近几年的工作报告，“民生水利”这一关键词频频提及，成为社会各界关注的热点。

图5 四川地区地壳主应变率 Fig.5 Crustal principal strain rate in Sichuan area

由图6和图7知，四川地区大部分区域地壳膨胀变形，其中较大的区域为龙门山断裂带和鲜水河断裂带南部以及四川西南部，最大膨胀率达到3.14×10-11/a；断裂带区域膨胀率为负，代表压缩变形，最小量为-8.59×10-11/a；最大剪切应变率在鲜水河断裂带的北西段和南东段附近达到峰值，幅度最大为8.33×10-11/a，另一高值区位于安宁河-则木河断裂带北部及龙门山断裂带南部。

图6 四川地区地壳面膨胀率 Fig.6 The expansion rate of crustal surface in Sichuan area

图7 四川地区地壳最大剪切应变率 Fig.7 Maximum shear strain rate of crust in Sichuan area

4 结 语

1)ITRF2008框架下四川省CORS站的东西方向平均速度为37.95 mm/a，欧亚框架下东西方向速度为13.76 mm/a；ITRF2008框架下南北方向平均运动速度为-11.05 mm/a，欧亚框架下南北方向速度为-1.46 mm/a。欧亚框架下的速度更准确，且水平方向的速度大于垂直方向。

2)ITRF2008框架下四川省CORS站的水平方向平均运动速率为38.72 mm/a，优势方向为S73.9°E；欧亚框架下水平方向的平均运动速率为14.40 mm/a，优势方向为S77.78°E，较ITRF2008框架下的速度降低了24.32 mm/a。CORS站整体向南东方向运动。

图5为四川区域地壳主应变率。可见，四川地区最小主应变率表现为主压应变，为-5.4×10-11/a；最大主应变率表现为主张应变，为7.31×1011/a。总体来看，川西地区和鲜水河断裂东北部主应变率相对较大，以主张应变为主；东部平原地区主应变率较小，以主压应变为主。三大主要断裂带上主应变率相对较大，以主张应变为主；东部平原地区主应变率较小，以主压应变为主。

3)四川省CORS站垂直方向平均运动速度为3.43 mm/a，大部分测站区域表现为隆升状态，位于东部平原的少数CORS站处于下降趋势。北部测站的隆升平均速度大于其他地区。

2.2 治疗前后两组Try-2、CER及AMY水平对比治疗后研究组Try-2、AMY水平均明显低于对照组，而CER水平明显高于对照组，差异均有统计学意义（P均＜0.05），见表1。

（二）新课导入课堂导入是整个教学过程中的一个有机组成部分,也是一个重要的教学环节。在语文课堂教学中,新课的导入成功与否，关键要看是否能够激发学生的求知欲。怎样导入新课，应该是：原则上求“趣”，形式上求“新”，内容上求“巧”，方式上求“活”，语言上求“精”。教师制作的微课要根据新课知识点设计新颖而实用的问题，以吸引学生的注意力，为新课的讲解做好铺垫。如人教版四年级上册的《梅花魂》，要紧紧抓住祖父对那幅梅花的热爱，以及梅花的高贵品格，微课设计上可以用反问形式激发学生的好奇心，不断指引他们往下走。

4)四川地区地壳平均最大主应变率为7.31×1011/a，最小主应变率为-5.4×10-11/a，且大部分川西区域地壳膨胀变形，面膨胀率低值区域和最大剪切应变率的高值区域分布在鲜水河以及其他断裂带附近。

致谢：数据资料由陆态网络及四川省地震局提供，在此表示感谢。

参考文献

[1] 曾波,张彦芬,姜卫平，等.山西CORS网基准站速度场分析[J].武汉大学学报：信息科学版,2012,37(12):1 401-1 404(Zeng Bo, Zhang Yanfen, Jiang Weiping, et al.Velocity Analysis of Shanxi CORS Fiducial Station[J].Geomatics and Information Science of Wuhan University,2012,37(12)：1 401-1 404)

[2] 张荣斗,黄丁发,廖华，等.基于GPS的四川地区现时地壳形变分析[J]，测绘科学,2008,33(6):30-33(Zhang Rongdou, Huang Dingfa, Liao Hua, et al. Analysis and Research of Sichuan Regional Crustal Deformation Based on GPS[J].Science of Surveying and Mapping,2008,33(6):30-33)

[3] 李婧. “陆态网络”基准站坐标时间序列变化特性分析[D].郑州：信息工程大学,2013(Li Jing. Analysis of the ‘CMONC Fiducial Stations’ Coordinate Time Series Variation Characteristics[D].Zhengzhou: Information Engineering University,2013

[4] 鄂栋臣,詹必伟,姜卫平，等.应用GAMIT/GLOBK软件进行高精度GPS数据处理[J].极地研究,2005,17(3):173-182(E Dongchen, Zhan Biwei, Jiang Weiping, et al. High-Precision GPS Data Procession by GAMIT/GLOBK[J].Chinese Journal of Polar Research ,2005,17(3):173-182)

[5] 王庆良,崔笃信,王文萍,等.川西地区现今垂直地壳运动研究[J].中国科学:地球科学,2008,38(5):598-610(Wang Qingliang,Cui Duxin,Wang Wenping, et al. Research on Vertical Crustal Movement in Western Sichuan[J]. Science China Earth Science,2008,38(5)：598-610)

[6] 金双根,朱文耀.亚洲板块运动欧拉参数的确定及其板内形变的分析[J].大地构造与成矿学,2002,26(2):138-142(Jin Shuanggen, Zhu Wenyao. Intra-Plate Deformation and Construction of Plate Motion Euler Parameter in the Asia Region[J]. Geotectonic et Metallogenia,2002,26(2)：138-142)

[7] 曹里,杨凯钧.多应变参数维度的川滇应变率场构造及分析[J].测绘科学,41(10):53-58(Cao Li，Yang Kaijun. The Structure of Strain and Characteristics Analysis of Sichuan-Yunnan Area under Multi-Strain Parameter Dimension[J]. Science of Surveying and Mapping, 41(10):53-58)

[8] 孙宏飞,余建平,张理想.利用GAMIT/GLOBK解算安徽CORS基准站速度场[C].2014年GPS、大地专业委员会学术年会,2014(Sun Hongfei,Yu Jianping,Zhang Lixiang. Using GAMIT / GLOBK to Solve Anhui Speedway of CORS Base Station[C]. 2014 GPS, the Earth Professional Committee Academic Annual, 2014)

[9] 姚宜斌.利用高精度复测GPS网进行中国大陆区域地壳运动特征分析[J].地球物理学进展,2008,23(4):1 030-1 037(Yao Yibin. Analysis on the Characteristics of Crustal Movement in Mainland China by Using High Accuracy Recosted GPS Network[J].Progress in Geophysics. 2008,23(4):1 030-1 037)