基于InSAR技术的张家口定点形变观测场地形变监测.pdf

高晨，孙丽娜，马栋，等.基于 InSAR 技术的张家口定点形变观测场地形变监测 J.华北地震科学，2023,41（3）：101-107.doi:10.3969/j.issn.10031375.2023.03.015.GAO Chen，SUN Lina，MA Dong，et al.Deformation Detection and Analysis of Fixed Point Deformation Observation at ZhangjiakouBased on InSAR TechnologyJ.North China Earthquake Sciences，2023,41（3）：101-107.doi:10.3969/j.issn.10031375.2023.03.015.基于 InSAR 技术的张家口定点形变观测场地形变监测高 晨1,2，孙丽娜1,2*，马 栋1,2，屈 曼1,2，侯晓真1,2，刘思宇1,2（1.河北红山巨厚沉积与地震灾害国家野外科学观测研究站,河北 邢台 054000；2.河北省地震局,石家庄 050021）摘要：针对张家口定点形变观测场地伸缩仪 NS 向观测曲线在 2020 年 6 月至 12 月出现加速上升的异常变化，以张家口定点形变观测场地为研究区域，基于 Sentinel-1 SAR 降轨数据并采用小基线集技术对场地进行形变监测，并获取了 20182020 年该区垂向年平均形变速率和形变时间序列。研究结果表明：张家口断裂两侧地表在垂向的相对运动速率小于 1 mm/a，断层活动性较弱；张家口断裂在 NNW-SSE 向拉张应力作用下，使得断层南盘地表相对北盘出现较为明显的沉降运动，导致张家口伸缩仪 NS 向观测出现张性加速趋势变化异常现象。关键词：伸缩仪；定点形变观测场地；SBAS-InSAR；张家口断裂；异常中图分类号：P315.725 文献标志码：A 文章编号：10031375（2023）03010107doi:10.3969/j.issn.10031375.2023.03.015 0 引言合 成 孔 径 雷 达 干 涉 测 量（InterferometricSynthetic Aperture Radar，简称 InSAR）技术是近年来新兴的一项空间大地测量技术，被广泛应用于地壳形变监测，具有全天候、无接触、覆盖范围广、高分辨率、高精度、强时效和具有穿透性等特点1-2。其中，SBAS-InSAR（Small Baseline Subset，简称 SBAS）技术是由 Berardino 和 Lanari 等人于 2002 年提出的一种形变时间时序分析技术3，该技术能够有效降低时空失相关、大气延迟的影响，对地表进行长时间序列的缓慢形变监测4。2020 年 6 月至 12 月，张家口定点形变观测场地伸缩仪 NS 向观测曲线开始出现大幅度加速上升的异常，呈现 NS 向加速拉张变化（图 1a）。经过河北省地震局河北地震台、中国地震台网中心和张家口地震监测中心站联合进行现场异常核实，发现SS-Y 型伸缩仪运行状态正常且周围不存在明显的地面荷载变化。将 20172020 年张家口伸缩仪NS 向观测数据与降雨量、气温和气压分别做相关性分析（图 1b1d），相关性系数分别为0.008 7、0.209 4、0.168 2。因此，认为张家口伸缩仪 NS 向异常与气温、气压、降水和地面荷载等因素不具有较强的相关性，经讨论将该项异常定为异常测项进行跟踪。张家口水管倾斜仪 NS 向观测曲线形态在2020 年 6 月 16 日至 2020 年 9 月 3 日期间同步出现大幅度南倾趋势变化（2019 年观测数据受干扰较为明显，未使用）（图 2）。由于地震前连续观测的形变测项（倾斜仪、伸缩仪等）会出现中短期前兆异常现象，且地震活动在地质上属于断裂运动的表现5-6。综上所述，张家口伸缩仪 NS 向观测异常可能与观测仪器附近活动断裂两侧地表在垂直方向上的差异性构造运动有关。本文拟使用 SBAS-InSAR 技术对张家口定点形变 观 测 场 地 的 垂 直 形 变 进 行 监 测，采 用 降 轨Sentinel-1 SAR 影像数据，分析研究区形变空间分布特征，揭示伸缩仪 NS 向形变异常的原因。收稿日期：2022-09-23基金项目：中国地震局地震科技星火计划攻关项目（XH200304）；河北地震科技星火计划项目（DZ20190424079）；河北省重点研发计划项目（20375404D）。第一作者简介：高晨（1989），男，工程师，主要从事 InSAR 和 GPS 形变监测与地球动力学研究，E-mail：gaoch_*通讯作者简介：孙丽娜（1982），女，高级工程师，主要从事地震活动性研究，E-mail：第 41 卷第 3 期华 北 地 震 科 学Vol.41，No.32023 年 7 月NORTH CHINA EARTHQUAKE SCIENCESJul.，2023 1 观测背景张家口定点形变观测场地地处阴山 EW 向构造带与山西构造带北段、张家口-渤海地震带（以下简称张渤带）西北段的交汇部位，附近的主要断裂为张家口断裂7。张家口断裂是张渤带西部一条重要断裂，是山西断陷带与燕山块体的构造分界之一8-9，控制了张家口-宣化盆地的北部边界，受到NWW 向张渤带和 NNE 向山西断陷带的复合作用。其北侧中生代地层抬升，形成低山丘陵；南侧下降，接受第四纪沉积多构成山间盆地8。张家口断裂全长约 70 km，由 NWW 向和近 EW 向两组多段断层组成的呈正断为主兼左旋走滑断裂，垂直滑动 速 率 约 为 1.33 mm/a，左 旋 走 滑 速 率 为 1.00.6 mm/a9-12。断裂由西、中、东 3 段组成：西段较短，长约 10 km，位于两条 NE 向断裂之间，沿河谷展布，隐伏于地下，构成低丘与高台的构造分界；中部主体段落起自万全镇西北，向东偏南延至西望山北东，长达 39 km，由 NWW 向和近 EW 向 2 组多条长短不等的次级断层斜列组成，为山盆之间的构造边界；东段由 NWW 向和 NW 向 2 条断裂组成，长约21 km，构造形迹不很清楚，止于其东的 NE 向断裂（图 3）10。张家口断裂的最新活动时代为晚更新世中晚期至全新世10,13。张家口地震监测中心站 SS-Y 伸缩仪为分钟值数字化记录方式，其基线长度：NS 向基线长 28.62 m，方 位 角 0；EW 向 基 线 长 28.22 m，方 位 角 90；NW 向基线长 29.81 m，方位角 45。仪器分辨率优于 11010，漂移量 1106/a，标定精度优于 1%。固体潮记录清晰，年变形态明显。63036912应变量20172018201920202021年份6303691250050应变量20172018201920202021年份 a 伸缩仪 NS 向观测曲线c 伸缩仪 NS 向观测与气温对比图1071076303691250403020100应变量降雨量/mm气压/Pa20172018201920202021年份63036912950920910940930900应变量20172018201920202021年份b 伸缩仪 NS 向观测与降雨量对比图d 伸缩仪 NS 向观测与气压对比图107107NS 向降雨量NS 向气压NS 向气温气温/图 1 张家口伸缩仪 NS 向观测曲线及其与气象三要素对比图Fig.1 The observation curve of tensometer in the NS direction at Zhangjiakou seismic station and it comparison with the three meteorologicalelements 100010020030012345678910 11 12倾斜量/()月份103201720182020 图 2 水管倾斜仪 NS 向观测曲线年变对比图Fig.2 The annual variation contrast chart of the water pipe tiltmeterin the NS direction102华 北 地 震 科 学41 卷 2 技术方法与数据处理2.1SBAS-InSAR 技术基本原理SBAS-InSAR 技术是 Berardino 和 Lanari 等3人在 2002 年提出的一种经典的基于多主影像的 InSAR时间序列分析方法。该技术主要使用具有较短时-空基线的影像对产生干涉图提高相干性，然后对选取的干涉对进行差分干涉处理并进行相位解缠，根据相干像元相位和观测时间的关系，采用奇异值分解（SVD）方法将多个小基线集合数据联合起来求解，求得影像序列间地表形变速度的最小范数最小二乘解2,14。2.2SBAS 时序分析本 文 使 用 欧 空 局（European Space Agency，ESA）获取的 Sentinel-1 共计 38 景降轨 SAR 数据，时间跨度从 20182020 年，成像模式为干涉宽幅模式（IW），详细数据参数见表 1。根据研究区范围（114.7 115.43E），（40.58 41.04N），对 SLC（Single Look Complex）影像进行裁剪。本文采用欧空局发布的对应的 AUX_POEORB 精密轨道参数文件；DEM 数据采用美国宇航局（NASA）提供的空间分 辨 率 为 30 m 的 SRTM DEM 数 据，用 于 针 对TOPS（Terrain Observation with Progressive）模式的增强谱分集（Enhanced Spectral Diversity，ESD）配准、模拟地形相位及地形相位去除15。表 1 研究区 Sentinel-1 卫星影像参数Table 1 Satellite image parameters of Sentinel-1 in study area 卫星 轨道方向 轨道号 开始获取日期 截止获取日期 极化方式 影像数S1BD472018-01-022 020-12-29VV38 本文选取 2018 年 9 月 11 日数据作为主影像，然后按照设定时间阈值为 300 d，空间基线阈值为200 m，并剔除干涉效果较差的干涉对来降低时空去相干的影响，最终共筛选出 87 个干涉对（图 4）。图 4 中的黄色点代表主影像，绿色点代表从影像，灰色线代表干涉像对。其中，最大时间基线为 264 d，最长空间基线为 118.19 m。对所有干涉像对进行干涉处理，包括干涉图生成，干涉图去平、滤波和相位解缠，并把所有像对配准到主影像上。其中，多视处理采用方位向视数为 1，距离向视数为 4；通过分析相干系数图，研究区内非植被覆盖区的相干系数均大于 0.2，因此将解缠相关系数阈值设置为 0.2，并采用 Delaunay MCF 方法对阈值大于 0.2 的像元点进行解缠；采用 SRTM DEM 数据模拟地形相位，以去除地形相位误差。在解缠后的干涉图像上选取远离形变区相对稳定的点作为 GCP 点，并采用多项式方法进行轨道精炼和重去平。第一次 SBAS 反演，先估算残余地形，并通过二次解缠优化全部干涉图；第二次 SBAS 反演，通过时空高通滤波和低通滤波对大气影响进行估计并去除大气相位误差，把 N4100N4050N4040N4030N050 km050 km11440E11500E11520E影像覆盖范围研究区特征点台站剖面线断层F1:张家口断裂;F2:万全断裂;F3:洋河断裂;F4:洗马林-水泉断裂 图 3 张家口定点形变观测场地附近断层分布图Fig.3 Fault distribution near the fixed point deformation observation site at Zhangjiakou Station 15010050050100150空间基线/m20180101201810022019070320200402 20210101日期 图 4 影像干涉对组合图Fig.4 Image Interference Pair Composition3 期高晨，等：基于 InSAR 技术的张家口定点形变观测场地形变监测103 所有 LOS 向形变序列和年平均形变速率转换到地理坐标系下。最后根据投影关系转换为垂直向年平均形变速率图，投影转换公式为：vUvLOScos(1)式中：vU表示垂直向年平均形变速率矢量；vLOS表示 LOS 向年平均形变速率矢量。2.3监测结果验证由于研究区内缺少外部地面实测数据（水准或者 GNSS），故选取相同数据的 PS-InSAR 垂向形变测量值结果16（图 5b）与 SBAS-InSAR 垂向形变测量值结果（图 5a）进行交叉互检来评价 SBAS-InSAR 垂向形变反演精度及可靠性17。考虑到研究区包括山地和平地地貌，在研究区范围内随机选取 10 个高相干点进行 SBAS-InSAR 年平均垂向形变速率精度定量评价（表 2）。研究以 SBAS-InSAR 反演得到的年平均垂向形变速率为参考，PS-InSAR 反演得到的年平均垂向形变速率为对照，采用均方根误差来评估两者反演数据结果的一致性：RMSE=vuuuuuutNi=1y(i)x(i)2N(2)式中：y 表示 SBAS-InSAR 年平均垂向形变速率值；x 表示 PS-InSAR 年平均垂向形变速率值；N 表示垂向形变速率值个数，i=1,2，N。高相干点的年平均垂向形变速率测量差异性指标：均方根误差为 1.15 mm/a，间接验证了 SBAS-InSAR 垂向形变反演的定量精度，且 SBAS-InSAR技术在研究区范围内形变监测可达 mm 级。3 地表沉降监测结果使用 SBAS-InSAR 技术对张家口定点形变观测场地的地表微小形变进行监测，获得了 20182020年期间垂直方向上的地表形变监测结果（图 5a）。地表形变监测结果显示，张家口定点形变观测场地附近区域垂直方向上未发现大范围显著的地表沉降或者隆升。张家口定点形变观测场地位于张家口断裂中段附近，因此本文重点研究张家口断裂中段。为了进一步定量分析张家口断裂中段的活动特征，绘制了 4 条垂直于断层走向的剖面线（图 3），即图 3 中的 A-A1D-D1，剖面线长度分别为 24 km、23 km、35 km 和 35 km。剖面线 A-A1 和B-B1 位于张家口台形变观测点西侧，剖面线 C-C1 和D-D1 位于张家口台形变观测点东侧，断层两侧地表年平均垂向形变速率几乎没有差异（图 6），剖面 表 2 高相干点垂向年平均形变速率测量值 mm/aTable 2 Measured value of vertical annual average deformation rate of high coherence point mm/a 点号SBAS-InSAR垂直形变速率PS-InSAR垂直形变速率点号SBAS-InSAR垂直形变速率PS-InSAR垂直形变速率11.348 0040.684 30260.113 2750.198 00220.796 6711.032 78670.792 1720.085 82230.475 3420.457 12480.616 3692.120 09240.480 8070.428 11191.333 2951.159 70650.173 4811.766 369100.419 7940.712 606 4100N4140N11500E11520E图例观测山河断裂速率/(mma1)6.1210.77010 kmF1:张家口断裂；F2:万全断裂;F3:洋河断裂a 研究区垂直形变场(SBAS-InSAR)4100N4140N11500E11520E图例观测山河断裂速率/(mma1)6.2911.64010 kmF1:张家口断裂；F2:万全断裂;F3:洋河断裂b 研究区垂直形变场(PS-InSAR)图 5 研究区垂直形变场Fig.5 Vertical deformation field of the study area104华 北 地 震 科 学41 卷 B-B1 断层两侧地表垂向形变速率与其他剖面略有差异，其南盘距离断层约 5 km 以外沉降速率变化较大，存在一个明显的下降趋势变化（图 6b），经光学影像核实此处为居民区，形变信息为非构造信号。各剖面位置断层两侧地表垂直方向上相对速率分别约为 0.38 mm/a、0.16 mm/a、0.8 mm/a、0.46 mm/a（图 6a6d），这与周月玲的结果相似10。分别在张家口定点形变观测场地的北盘和南盘分别选取 4 个特征点（图 3），S1 和 S2 均位于张家口市区西侧八角台山，S3 位于张家口市区东侧东山，S4 位于翟家庄村北基岩出露点，4 个点距断层直线距离分别约为 4 km、6 km、3.4 km、10 km；N1和 N3 位于佛爷山，N2 位于安家沟，N4 位于石人山，4 个点距离断层直线距离分别约为 10 km、6.5 km、3 km、5 km。以上特征点均分布在基岩出露区域且相干系数大于 0.9，得到垂直方向上累计形变时间序列 结 果（图 7a 7b）。由 于 Sentinel-1 卫 星 在2020 年 3 月 16 日至 7 月 26 日期间缺少覆盖研究区的降轨数据，因此断层北盘和南盘的特征点在这个期间垂直方向上形变量数值较大，在图 7a7b 中表现为折线斜率较大。图 7a 显示，张家口断裂北盘在 20182020 年期间地表在垂直方向上一直呈现隆升运动趋势。图 7b 显示，张家口断裂南盘在2018 年 19 月期间地表在垂直方向上形变量不明显，2018 年 10 月至 2020 年 7 月前运动趋势与张家口断裂北盘基本保持一致，表现为隆升运动；2020年 712 月期间，张家口断裂南盘地表运动趋势发生较为明显的转折，表现为沉降运动，该时间段南盘地表运动趋势与 2018 年和 2019 年同期的隆升运动相反，且张家口断裂南盘地表运动趋势发生转折的时间与张家口台伸缩仪 NS 向地表呈现加速拉张变化（图 1a）和水管倾斜仪 NS 向出现大幅度向南倾斜（图 2）异常现象的发生和持续时间基本一致。4 结果分析与讨论根据前人已有研究结果，张家口-宣化盆地现今的区域构造应力场以 NEE-SWW 向挤压、NNW-SSE 向拉张为主要特征11，张家口断裂作为张家口盆地北缘边界的主控断裂，在第四纪期间持续活动，其中部段落活动强烈8；根据张家口至张北1998、2002 和 2006 年各一期、2013 年 3 期精密水准观测数据资料分析结果显示，水准测线整体在张家口断裂附近形成垂直高梯度带区，并以张家口断裂为界表现为北升南降趋势，且 2013 年以后断层两侧 101234115051015速率/(mma1)距离断层的距离/kma A-A1 两侧地表垂向形变速率NS32101234201510505101520速率/(mma1)距离断层的距离/kmc C-C1 两侧地表垂向形变速率NS10123151050511速率/(mma1)距离断层的距离/kmb B-B1 两侧地表垂向形变速率NS3210123201510505101520速率/(mma1)距离断层的距离/kmd D-D1 两侧地表垂向形变速率图 6 张家口断裂剖面线形变示意图NS 图 6 张家口断裂剖面线形变示意图Fig.6 Deformation profiles of the Zhangjiakou fault3 期高晨，等：基于 InSAR 技术的张家口定点形变观测场地形变监测105 垂直形变速率增大，同时张家口断裂的水平运动弱于垂直运动9,18。根 据 上 述 对 于 张 家 口 断 裂 的 研 究 结 果，2020 年 612 月，张家口伸缩仪 NS 向观测曲线形态与往年同期呈张性趋势变化基本一致，但是曲线上升速率显著高于往年同期，出现加速张性变化趋势异常，张家口水管倾斜仪 NS 向观测曲线形态同时出现加速向南倾斜变化。针对上述异常，通过与利用 SBAS-InSAR 技术获取的张家口断裂两侧特征点垂直方向上累计形变时间序列曲线进行对比，发现在 2020 年 612 月，张家口南盘地表相对于北盘呈较为明显沉降运动趋势，与 2018 年和 2019 年同时间段南盘地表运动趋势相反。根据张家口断裂附近现今区域应力场主要特征并结合特征点累计形变时间序列，该断层在 2020 年 612 月间可能受NNW-SSE 向拉张应力作用，使断层南盘地表相对于北盘在垂直方向上出现较为明显的沉降运动，因此导致了张家口伸缩仪 NS 向观测曲线在该时段出现张性变化加速趋势变化和水管仪 NS 向观测曲线出现加速向南倾斜趋势变化的异常。5 结论基于欧空局提供的 Sentinel-1 降轨 SAR 影像数据，使用 SBAS-InSAR 技术对张家口定点形变观测场地地表形变进行监测研究，得到 20182020 年期间该场地及附近区域地表在垂直方向上形变结果。1）张家口断裂两侧地表在垂直方向的相对运动速率小于 1 mm/a，断层活动性较弱。2）张家口断裂在在 2020 年 612 月间可能受NNW-SSE 向区域拉张应力作用，使断层南盘地表相对于北盘在垂直方向出现较为明显的沉降运动趋势，导致张家口伸缩仪 NS 向观测出现张性加速趋势变化异常。致谢感谢欧空局提供的 Sentinel-1 雷达数据。2018-01-022018-02-072018-03-032018-03-272018-04-202018-05-142018-06-072018-07-012018-07-252018-08-182018-09-112018-10-052018-10-292018-11-222018-12-162019-01-092019-02-022019-02-262019-03-222019-04-152019-05-092019-06-022019-08-132019-09-182019-10-242019-11-172019-12-112020-01-282020-02-212020-03-162020-07-022020-07-262020-08-192020-09-122020-10-062020-10-302020-11-232020-12-29543210累计形变量/mmN1N2N3N4日期a 北盘垂向形变时间序列曲线2018-01-022018-02-072018-03-032018-03-272018-04-202018-05-142018-06-072018-07-012018-07-252018-08-182018-09-112018-10-052018-10-292018-11-222018-12-162019-01-092019-02-022019-02-262019-03-222019-04-152019-05-092019-06-022019-08-132019-09-182019-10-242019-11-172019-12-112020-01-282020-02-212020-03-162020-07-022020-07-262020-08-192020-09-122020-10-062020-10-302020-11-232020-12-295432101234累计形变量/mmN1N2N3N4日期b 南盘垂向形变时间序列曲线 图 7 特征点垂向形变时间序列曲线Fig.7 Vertical deformation time series curves of characteristic points106华 北 地 震 科 学41 卷 参考文献：朱建军,李志伟,胡俊.InSAR 变形监测方法与研究进展 J.测绘学报,2017,46(10):1717-1733.1赵强,苏琴.基于 InSAR 技术的冕宁跨断层场地形变监测与分析 J.地震研究,2018,41(3):361-367.2Berardino P,Fornaro G,Lanari R,et al.A new algorithm for surface deformation monitoring based on small baseline differential SARinterferogramsJ.IEEE Transactions on Geoscience and Remote Sensing,2002,40(11):2375-2383.3李达,邓喀中,高晓雄,等.基于 SBAS-InSAR 的矿区地表沉降监测与分析 J.武汉大学学报:信息科学版,2018,43(10):1531-1537.4车兆宏,谢觉民,张晶,等.张北 6.2 级地震地形变前兆演化特征 J.地震,1999,19(4):315-322.5吴中海,周春景,冯卉,等.青海玉树地区活动断裂与地震 J.地质通报,2014,33(4):419-469.6马栋,陈建国,郑爽,等.荷载对钻孔应变观测影响的实例分析 J.地震研究,2014,37(1):79-85.7杨歧焱,彭远黔,周月玲,等.综合地质-地球物理方法在山间盆地断裂探测中的应用以张家口断裂为例 J.地球物理学进展,2016,31(6):2451-2457.8孟宪纲,占伟,朱文武,等.晋冀蒙交界地区地壳形变运动和地震活动性研究 J.地球物理学进展,2017,32(5):1907-1914.9周月玲,尤惠川.张家口断裂第四纪构造变形与活动性研究 J.震灾防御技术,2010,5(2):157-166.10张红艳,谢富仁,荆振杰.京西北盆岭构造区现代构造应力场的非均匀特征 J.地球物理学报,2009,52(12):3061-3071.11陈长云,贺建明,李腊月,等.基于跨断层和 GPS 资料综合分析张渤带运动特征 J.华南地震,2016,36(3):17-28.12尤惠川,邵翠茹,杨歧焱,等.张家口断裂发现第四纪晚期活动证据 J.震灾防御技术,2008,3(4):474-477.13冯文凯,顿佳伟,易小宇,等.基于 SBAS-InSAR 技术的金沙江流域沃达村巨型老滑坡形变分析 J.工程地质学报,2020,28(2):384-393.14吴文豪.哨兵雷达卫星 TOPS 模式干涉处理研究 D.武汉:武汉大学,2016.15高胜,曾琪明,焦健,等.永久散射体雷达干涉研究综述 J.遥感技术与应用,2016,31(1):86-94.16何海英,陈彩芬,陈富龙,等.张家口明长城景观廊道 Sentinel-1 影像 SBAS 形变监测示范研究 J.国土资源遥感,2021,33(1):205-213.17孟宪纲,陈聚忠,张俊青,等.张家口至张北垂直形变异常及分析 J.大地测量与地球动力学,2015,35(1):13-15.18DeformationDetectionandAnalysisofFixedPointDeformationObservationatZhangjiakouBasedonInSARTechnologyGAO Chen1,2,SUN Lina1,2*,MA Dong1,2,QU Man1,2,HOU Xiaozhen1,2,LIU Siyu1,2（1.Hebei Hongshan National Observatory on Thick Sediments and Seismic Hazards,Xingtai 054000,China;2.Hebei Earthquake Agency,Shijiazhung 050021,China）Abstract:The NS direction of tensometer of Zhangjiakou fixed-point deformation observation site has acceleratedrising anomaly from June to December 2020.In view of this anomaly,we take Zhangjiakou fixed-point deformationobservation site as the research area,based on sentinel-1 SAR descending imaging and small baseline subset(SBAS)technology for InSAR time series analysis.we obtains the vertical annual average deformation rate and deformationtime series in this area from 2018 to 2020.The results show that the vertical relative movement rate of the surface onboth sides of Zhangjiakou fault is less than 1 mm/a,and the fault activity is weak.Under the action of NNW-SSEtensile stress,the southern side surface of Zhangjiakou fault appears obvious subsidence movement relative to thenorthern side,this is the reason why the observation curves of the NS component of tensometer in Zhangjiakou seismicstation showed the abnormal change of tensile acceleration trend.Keywords:tensometer；fixed point deformation observation site；SBAS-InSAR；Zhangjiakou fault；Abnormal3 期高晨，等：基于 InSAR 技术的张家口定点形变观测场地形变监测107