快速检索        
  武汉大学学报·信息科学版  2016, Vol. 41 Issue (4): 535-540

文章信息

袁鹏, 孙宏飞, 秦昌威, 张理想
YUAN Peng, SUN Hongfei, QIN Changwei, ZHANG Lixiang
安徽CORS参考站三维速度场分析
Analysis of Anhui CORS Reference Stations 3D Velocity Field
武汉大学学报·信息科学版, 2016, 41(4): 535-540
Geomatics and Information Science of Wuhan University, 2016, 41(4): 535-540
http://dx.doi.org/10.13203/j.whugis20140552

文章历史

收稿日期: 2015-03-02

安徽CORS参考站三维速度场分析
袁鹏1 , 孙宏飞2, 秦昌威3,4, 张理想2    
1. 武汉大学卫星导航定位技术研究中心, 湖北 武汉, 430079;
2. 安徽省基础测绘信息中心, 安徽 合肥, 230031;
3. 武汉大学测绘学院, 湖北 武汉, 430079;
4. 湖北省国土测绘院, 湖北 武汉, 430010
摘要: 利用GAMIT/GLOBK软件处理安徽连续运行参考站(Continuously Operating Reference Stations,CORS)2011-11~2013-09近2 a的数据,获取了参考站的三维速度场,并分析了安徽省地壳水平运动和陆地垂直运动。结果表明,安徽CORS参考站在ITRF2008框架下的平均水平运动速率为34.44 mm/a,方向为E18.01°S;相对于欧亚板块的平均水平运动速率为6.13 mm/a,方向为E7.31°S;在安徽CORS整体无净旋转基准下,淮北平原的参考站存在1~4 mm/a的水平相对运动。郯庐断裂带在嘉山至泗县一带表现为左旋运动,嘉山以南运动较弱。在垂直方向上,安徽CORS参考站大致以淮河为界表现出南升北降的态势。淮河以南平均隆升速率为2.72 mm/a,淮河以北平均沉降速率为2.98 mm/a。该地区的地面沉降可能与过量抽取地下水和矿产开发等人为因素有关。
关键词: 连续运行参考站网络     速度场     地壳运动     地面沉降    
Analysis of Anhui CORS Reference Stations 3D Velocity Field
YUAN Peng1, SUN Hongfei2, QIN Changwei3,4, ZHANG Lixiang2    
1. Research Center of GNSS, Wuhan University, Wuhan 430079, China;
2. Surveying and Mapping Geomatics Center of Anhui, Hefei 230031, China;
3. School of Geodesy and Geomatics, Wuhan University, Wuhan 430079, China;
4. Hubei Institute of Land Surveying and Mapping, Wuhan 430010, China
First author: YUAN Peng, PhD candidate, specializes in geodesy and geodynamics. E-mail: yuanpeng@whu.edu.cn
Foundation support: The National Natural Science Foundation of China, No. 41374033; the National High Technology Research and Development Program of China (863 Program), No. 2012AA12A209.
Abstract: The most important causes of horizontal reference stations are human factors such as excessive pumping of groundwater and immoderate mineral exploitation. GPS data from Anhui Operating Reference Stations(AHCORS) between 2011-11 and 2013-9 was continuously processed the using GAMIT/GLOBK software. Then, based on these results, the 3D velocities of the reference stations were determined, and crustal horizontal movement and land vertical movement was analyzed in detail. It was found that the mean horizontal velocity of AHCORS reference stations was 34.44 mm/a in the direction of E18.01°S under ITRF2008 reference frame and 6.13 mm/a in the direction of E7.31°S with regard to the Eurasia Plate. Stations located in Huaibei Plain have moved in horizontal with 1~4 mm/a relative to No-Net-Rotation of the whole AHCORS. Tancheng-Lujiang Fault showed sinistral strike slip feature between Jiashan and Sixian, with very weak movement in the south of Jiashan. As for the vertical direction, AHCORS reference stations uplift in the south of the Huai River, but subside in the north of it. The average uplift rate of the stations located in the south of the Huai River is 2.72 mm/a, while the average subside rate of the stations located in the north of the Huai River is 2.98 mm/a.
Key words: CORS     velocity     crust motion     land subsidence    

全球定位系统(GPS)观测可提供高精度、大范围和准实时的地壳运动定量数据,已广泛应用于地壳水平运动研究[1, 2]等领域。通过长期观测并采用精密的数据处理策略,利用GPS技术也能监测陆地垂直运动[3]。连续运行参考站网络(Continuously Operating Reference Stations,CORS)是一种重要的地理空间信息基础设施,它能够提供长期、连续、稳定的GPS观测数据。安徽省卫星定位综合服务系统(AHCORS)于2011年建成,包含均匀分布全省的50余个参考站,为研究安徽省现今三维运动状态积累了大量基础数据。

安徽省地处华北、扬子断块和秦岭大别断褶区3个大地构造单元的接壤地带,郯庐断裂带斜穿而过,长江、淮河横贯其中,是一个地质构造比较复杂的地区,地质灾害较为严重。尤其是郯庐断裂带的活动与淮北平原的地面沉降引人关注,不少学者已对此展开了研究[4, 5, 6, 7, 8, 9]。本文获取了AHCORS参考站的三维速度场,更为精确、完整地描述和分析了安徽省地壳水平运动状态,并第一次定量地给出了由GPS技术揭示的淮北平原地面沉降状况。

1 数据处理

本文采用了AHCORS 50个参考站2011年11月~2013年9月每隔7天一组的观测数据进行处理。此外,为了获得AHCORS参考站在全球框架下的速度场,本文还解算了全球73个国际全球卫星导航定位系统服务(International GNSS Service,IGS)站的历史同期观测数据。

采用GAMIT/GLOBK(Ver.10.5)软件[10]进行数据处理,主要分为三步。

(1)利用GAMIT获得AHCORS参考站和国内及周边IGS站的单日松弛解。处理策略主要包括:①同时解算卫星轨道和测站坐标;②使用双差观测值,并根据验后相位残差重新定权;③使用igs08.atx绝对天线相位中心改正模型[11];④对流层延迟改正选用维也纳映射函数(Vienna Mapping Function 1,VMF1)[12],每2 h估计1个天顶湿延迟参数,每天估计两对水平梯度参数;⑤顾及二、三阶电离层延迟的影响,垂直总电子含量(Vertical Total Electron Content,VTEC)值来自于欧洲定轨中心(Center for Orbit Determination in Europe,CODE),地磁场模型采用第11代国际地磁参考场(International Geomagnetic Reference Field 11,IGRF11)[13];⑥顾及固体潮、极潮、海潮的影响,海潮改正采用有限元解2004(Finite Element Solutions 2004,FES2004)格网模型[14];⑦使用IGS精密星历和地球自转参数,并给予适当约束。

(2)按照与步骤(1)相同的策略解算全球73个IGS站,获得单日松弛解。

(3)利用GLOBK综合各网单日松弛解,并以IGS站在国际地球参考框架2008(International Terrestrial Reference Frame 2008,ITRF2008)下的坐标和速度为基准,估算AHCORS参考站的速度。

2 速度场分析 2.1 水平速度场

表 1给出了ITRF2008框架下AHCORS参考站的三维速度估值和中误差统计。表 1表明,AHCORS参考站在ITRF2008框架下N、E方向速度的精度较高,平均值为0.2 mm/a,明显优于U方向。但GPS测站运动存在有色噪声,采用有色噪声假设得到的速度不确定度可达到白噪声假设的4~5倍[15]。因此,本文采用白噪声假设得到的速度不确定度可能过于乐观。参考站平均水平运动速率为34.44 mm/a,方向为E18.01°S,与文献[1]一致。

表 1 ITRF2008框架下AHCORS的三维速度估值及中误差统计/(mm·5a-1) Tab. 1 Horizontal and Vertical Velocities and Their Uncertainties of AHCORS in ITRF2008/(mm·5a-1)
VEVNVUσEσNσU
最大值34.83-7.598.220.240.230.86
最小值31.15-15.21-13.500.190.190.65
平均值32.73-10.660.780.210.200.73

任何运动都是相对的,描述地壳运动必须借助确定的参考基准。不同参考基准下的速度具有不同含义,能够突出不同的运动特征。ITRF框架下的速度是相对于整个地球的水平板块运动无净旋转条件而言的。它可以分为三部分[1, 16],分别研究区域内部相对运动,研究区域整体相对于欧亚板块的刚性运动以及欧亚板块相对于ITRF框架的刚性运动。相对于欧亚板块的速度场扣除了欧亚板块相对于ITRF框架的刚性运动。区域整体无净旋转速度场进一步扣除了研究区域整体相对于欧亚板块的刚性运动,仅包含区域内部的相对运动。

为分析AHCORS参考站相对于欧亚板块的运动状态,采用ITRF2008板块运动模型[17],对ITRF2008框架下的速度进行转换。结果表明,AHCORS参考站相对于欧亚板块的平均运动速率为6.13 mm/a,方向为E7.31°S,与文献[1]一致。这反映出安徽省整体上具有在统一的大陆动力学环境下连续变形的特性。

为进一步分析AHCORS参考站的水平差异性运动,求取了参考站在AHCORS整体无净旋转基准下的速度场,如图 1所示。本文中,误差椭圆代表 95%置信度。具体步骤为:基于AHCORS参考站在ITRF2008框架下的实测速度,采用最小二乘法估计参考站整体相对于ITRF2008框架的欧拉矢量,然后将实测速度与欧拉矢量模型反算速度作差即可。

图 1可以看出,淮河是安徽省一条明显的水平速度分界线。淮河以南的大部分参考站水平相对运动速度很小。淮河以北除萧砀地区外的参考站,均存在1~4 mm/a的相对运动。淮北平原断裂构造发育,新构造和地震活动频繁,固镇地区(BBGZ)于1979年发生5.0级地震[4]。位于淮河构造变形带和郯庐断裂带中段(潍坊-嘉山)与南段(嘉山-广济)过渡区域的嘉山(CZMG)、五河(BBWH)、泗县(SZSX)、灵璧(SZLB)和固镇(BBGZ)等站表现出左旋运动特征,与文献[5]一致。郯庐断裂带南段大部分参考站速度很小,仅合肥(CHJU)和舒城(LASC)存在约1 mm/a的压性运动。总体来说,在安徽境内,郯庐断裂带仅嘉山地区较为活跃,嘉山以南则相对较弱。这证实了郯庐断裂带第四纪活动具有分段特征,且中、南段运动特征差异明显[5, 6]的推论。因此,郯庐断裂带中、南段过渡区域的地壳活动应予以特别关注。

图 1 AHCORS参考站在整体无净旋转基准下的水平速度场 Fig. 1 Horizontal Velocity Field of AHCORS Relative to No-net-rotation of the Whole AHCORS
2.2 垂直速度场

图 2展示了AHCORS参考站的垂直速度场及其等值线。图 2表明,安徽省陆地垂直运动的总体态势为南升北降,淮河是升降运动的分界线。淮河以南表现出隆升趋势,参考站垂直运动速度的平均值为2.72 mm/a。其中,江淮丘陵和黄山地区的参考站隆升较为明显,约为2~5 mm/a,而位于大别山区和沿江平原的参考站则较弱。淮北平原的大部分参考站表现出沉降趋势,平均沉降速度为2.98 mm/a,其西北部地区沉降速率较大,平均值为6.03 mm/a,最大值出现在宿州砀山(SZDS),达到13.50 mm/a。

图 2 AHCORS参考站垂直速度场和等值线 Fig. 2 Vertical Velocity Field and Its Contours of AHCORS

文献[7, 8]采用水准测量的方法研究了安徽省陆地垂直运动,研究得出的整体升降趋势与本文基本一致,数值也较为接近,但局部地区存在差异。在淮北平原西北部,本文获得的沉降速率比文献[7, 8]大2~6 mm/a;在江淮丘陵,本文获得的隆升速率与文献[7]一致,与文献[8]相反。这些差异可能是多方面原因造成的。技术上,本文采用GPS技术,而文献[7, 8]采用的是水准测量方法,易受误差积累的影响。空间上,AHCORS参考站分布均匀,而文献[7, 8]的水准点沿国家一、二等水准路线布设,分布不均,难以全面地反映安徽省陆地垂直运动的细节信息。时间上,本文的研究时段为2011~2013年,文献[7, 8]的研究时段分别为20世纪50~80年代和1951~1999年,而淮北平原20世纪80年代中后期才开始发生较为严重的地面沉降。此外,文献[9]发现江淮丘陵存在着伴随郯庐断裂带的张压交替运动而产生的周期性升降运动,这也可能是本文与文献[8]得到的该地区升降趋势相反的原因之一。但这种周期性升降运动是否存在还有待更长时间的观测来证实。

本文得到的AHCORS参考站沉降速率见表 2。结合表 2图 2来看,除安庆潜山(AQQS)外,所有沉降的参考站都位于淮北平原。沉降速率较大的参考站主要分布在阜阳、亳州、宿州砀山和蚌埠固镇等地。

表 2 AHCORS参考站沉降速率及其中误差/(mm·5a-1) Tab. 2 Subside Velocities and Their Uncertainties of AHCORS in ITRF2008 /(mm·5a-1)
测站所在地区沉降速率中误差
SZDS宿州砀山13.500.68
FYLQ阜阳临泉10.590.65
BBGZ蚌埠固镇7.210.68
FYJS阜阳界首6.210.67
BZBZ亳州市5.630.81
BZGY亳州涡阳5.050.72
FYFY阜阳市4.080.74
FYYS阜阳颍上3.620.73
SZSX宿州泗县3.570.69
FYFN阜阳阜南2.890.7
BZLX亳州利辛2.690.69
AQQS安庆潜山0.410.8

淮北平原主要由淮河和黄河冲积物组成,松散层厚度可达数百m。该区人口占全省的43%,但水资源仅占全省的约20%,城乡工农业生产和生活用水主要来自于地下水[18]。由于过度抽取地下水,淮北平原发生了明显的地面沉降。淮北平原西部的阜阳和亳州的孔隙承压水处于同一承压水盆地中,由于开采强度大、补给条件差,承压水水位持续下降,造成粘土层压密释水,从而引发严重的地面沉降[19]。淮北平原最北端的宿州砀山水文地质条件与平原西部类似,浅层地下水开采率超过50%[18, 19],也形成了超采漏斗,产生了地面沉降。蚌埠固镇矿产资源丰富,矿产资源过量开采可能是引发该地区地面沉降的重要因素。而且,该地区地壳活动较强,不排除构造因素造成沉降。

值得注意的是,安庆望江站(AQWJ)的上升速率达到8.22±0.72 mm/a。该站地处宿松-枞阳断裂东南侧,在新构造上属于长江下游安庆弱上升区[20],近期曾发生2005年九江Ms5.7和2011年安庆Ms4.8等多次4级以上地震[21]。但周边测站却没有明显的上升趋势,因此怀疑该站的上升可能与该站当地地壳异常活动有关。

3 结语

本文采用GAMIT/GLOBK软件处理了AHCORS参考站2011~2013年的数据以及全球73个IGS站的历史同期数据,确定了AHCORS参考站在ITRF2008框架、相对于欧亚板块和研究区域整体无净旋转基准下的三种水平速度场及垂直速度场,并进行了分析。结果表明,AHCORS参考站在ITRF2008框架下的平均水平运动速率为34.44 mm/a,方向为E18.01°S;在ITRF2008板块运动模型下,相对于欧亚板块的平均水平运动速率为6.13 mm/a,方向为E7.31°S;在AHCORS整体无净旋转基准下,淮北平原的参考站存在1~4 mm/a的相对运动。郯庐断裂带中、南段运动特征差异明显,在嘉山至泗县一带表现为左旋运动,在嘉山以南运动较弱。在垂直方向上,AHCORS参考站以淮河为界,表现出南升北降的态势。淮河以南平均隆升速率为2.72 mm/a。其中,江淮丘陵地区隆升明显,约为2~5 mm/a。淮北平原平均沉降速率为2.98 mm/a,其西北部地面沉降较为严重,平均沉降速率达到6.03 mm/a。这主要是由于过量抽取地下水和开发矿产等人为因素引发的。

参考文献
[1] Wang Qi. GPS Determination of Present-day Crustal Deformation in China:Velocity Field and Tectonic Implication[D]. Wuhan:Wuhan University, 2004(王琪. 用GPS监测中国大陆现今地壳运动:变形速度场与构造解释[D]. 武汉:武汉大学, 2004)
[2] Qu Wei, Zhang Qin, Wang Qingliang, et al. Present Tectonic Motion Feature of Shanxi Graben Based on the GPS Monitoring Data[J]. Geomatics and Information Science of Wuhan University, 2012,37(02):195-198,246,259(瞿伟,张勤,王庆良,等. 利用GPS观测资料分析山西地堑系现今地壳构造活动特征[J]. 武汉大学学报·信息科学版, 2012,37(02):195-198,246,259)
[3] 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):1401-1404(曾波, 张彦芬, 姜卫平,等. 山西CORS网基准站速度场分析[J]. 武汉大学学报·信息科学版, 2012, 37(12):1401-1404)
[4] Liu Dongwang, Liu Zemin, Shen Xiaoqi, et al. The Characteristics of Modern Structure Stress Field in the Huaihe Structure Deformation Belt, Anhui and Its Neighborhood Block[J].Earthquake Research in China, 2004, 20(4):364-371(刘东旺, 刘泽民, 沈小七,等. 安徽淮河构造变形带及邻近块体现代构造应力场特征[J].中国地震, 2004, 20(4):364-371)
[5] Liu Xiaoxia, Jiang Zaisen, Wu Yanqiang, et al. Motion and Deformation State of the Tancheng-Lujiang Fault Zone Derived from GPS Data[J]. Earthquakes, 2012, 32(4):1-10(刘晓霞, 江在森, 武艳强,等. 利用GPS资料研究郯庐带现今运动及变形状态[J].地震, 2012, 32(4):1-10)
[6] Shi Wei, Zhang Yueqiao, Dong Shuwen, et al. Quaternary Activity and Segmentation Behavior of the Middle Portion of the Tan-Lu Fault Zone[J].Acta Geoscientia Sinica, 2003, 24(1):11-18(施炜,张岳桥,董树文,等.郯庐断裂带中段第四纪活动及其分段特征[J].地球学报, 2003, 24(1):11-18)
[7] Feng Haojian, Gu Dansheng, Zhang Li, et al. The Research on the Earth Crust Vertical Movement Characteristic and Mechanism in Eastern China[J].Acta Geodaetica et Cartographica Sinica, 1998,27(1):16-23(冯浩鉴, 顾旦生, 张莉,等. 中国东部地区地壳垂直运动规律及其机制研究[J].测绘学报, 1998,27(1):16-23)
[8] Dong Hongwen, Gu Dansheng, Li Guozhi, et al. Research on Vertical Recent Crustal Movement of the Mainland of China[J]. Acta Geodaetica et Cartographica Sinica, 2002,31(2):100-103(董鸿闻,顾旦生,李国智,等. 中国大陆现今地壳垂直运动研究[J].测绘学报, 2002,31(2):100-103)
[9] Yang Guohua, Han Yueping. Recent Crustal Vertical Movement in Jiangsu-Shandong-Anhui Area and Characteristics of Activity in Tancheng-Lujiang Fault Zone[J]. Crustal Deformation and Earthquake, 1997, 17(01):39-44(杨国华,韩月萍.苏鲁皖地区现今地壳垂直运动及郯庐断裂带的活动特征[J]. 地壳形变与地震, 1997, 17(01):39-44)
[10] MIT. GAMIT/GLOBK 10.5[OL]. http://www-gpsg.mit.edu/-simon/gtgk/docs.htm, 2015
[11] IGS. IGSMAIL-6354 Upcoming Switch to IGS08/igs08.atx[OL]. http://igscb.jpl.nasa. gov/pipermail/igsmail/2011/006346.html, 2015
[12] Boehm J,Werl B, Schuh H. Troposphere Mapping Functions for GPS and Very Long Baseline Interferometry from European Centre for Medium-range Weather Forecasts Operational Analysis Data[J]. Journal of Geophysical Research:Solid Earth, 2006, 111(B02406):1-9
[13] Petrie E J, King M A, Moore P, et al. Higher-order Ionospheric Effects on the GPS Reference Frame and Velocities[J]. Journal of Geophysical Research:Solid Earth,2010, 115(B03417):1-8
[14] Lyard F, Lefevre F, Letellier T, et al. Modeling the Global Ocean Tides:Modern Insights from FES2004[J]. Ocean Dyn, 2006,56(5-6):394-415
[15] Santamaría-Gómez A, Bouin M N, Collilieux X, et al. Correlated Errors in GPS Position Time Series:Implications for Velocity Estimates[J]. Journal of Geophysical Research:Solid Earth, 2011, 116:B01405
[16] Jiang Zaisen, Liu Jingnan. The Method in Establishing Strain Field and Velocity Field of Crustal Movement Using Least Squares Collocation[J]. Chinese J. Geophys, 2010, 53(05):1109-1117(江在森, 刘经南. 应用最小二乘配置建立地壳运动速度场与应变场的方法[J].地球物理学报, 2010, 53(05):1109-1117)
[17] Altamimi Z, Métivier L, Collilieux X, et al. ITRF2008 Plate Motion Model[J]. Journal of Geophysical Research:Solid Earth, 2012, 117:B07402
[18] Zhao Hui, Chen Xiaofeng. Huaibei Plain Groundwater Resources Development and Protection Research[J]. Express Water Resources & Hydropower Information, 2010, 31(2):20-22(赵晖,陈小凤.淮北平原地下水资源开发利用及保护研究[J].水利水电快报,2010,31(2):20-22)
[19] Tao Yuezan, Xi Daoying. Groundwater Evaluation and Pumping-limit Planning of Anhui Huaibei Plain[J]. Geotechnical Investigation & Surveying, 2006,1:30-33,22(陶月赞,席道瑛.安徽淮北平原地下水超采评价与限采规划[J].工程勘察, 2006,1:30-33,22)
[20] Zhai Hongtao, Zheng Yingping, Li Guang, et al. Study on Recent Active Age of Susong-Zongyang Fault and Its Seismic Hazard in the Future[J]. Journal of Disaster Prevention and Mitigation Engineering, 2010,05:567-571(翟洪涛,郑颖平,李光,等. 宿松-枞阳断裂最新活动时代及未来地震危险性研究[J]. 防灾减灾工程学报,2010,05:567-571)
[21] China Earthquake Networks. History Search[OL]. http://www.ceic.ac.cn/history,2015(中国地震台网. 历史查询[OL]. http://www.ceic.ac.cn/history, 2015)