ISSN 1003-8035 CN 11-2852/P
  • 中国科技核心期刊
  • CSCD收录期刊
  • Caj-cd规范获奖期刊
  • Scopus 收录期刊
  • DOAJ 收录期刊
  • GeoRef收录期刊
欢迎扫码关注“i环境微平台”

基于改进时序InSAR技术的东莞地面沉降时空演变特征

戴真印, 刘岳霖, 张丽平, 张贤

戴真印,刘岳霖,张丽平,等. 基于改进时序InSAR技术的东莞地面沉降时空演变特征[J]. 中国地质灾害与防治学报,2023,34(1): 58-67. DOI: 10.16031/j.cnki.issn.1003-8035.202112028
引用本文: 戴真印,刘岳霖,张丽平,等. 基于改进时序InSAR技术的东莞地面沉降时空演变特征[J]. 中国地质灾害与防治学报,2023,34(1): 58-67. DOI: 10.16031/j.cnki.issn.1003-8035.202112028
DAI Zhenyin,LIU Yuelin,ZHANG Liping,et al. Spatial-temporal evolution characteristics of land subsidence in Dongguan City based on improved InSAR technology[J]. The Chinese Journal of Geological Hazard and Control,2023,34(1): 58-67. DOI: 10.16031/j.cnki.issn.1003-8035.202112028
Citation: DAI Zhenyin,LIU Yuelin,ZHANG Liping,et al. Spatial-temporal evolution characteristics of land subsidence in Dongguan City based on improved InSAR technology[J]. The Chinese Journal of Geological Hazard and Control,2023,34(1): 58-67. DOI: 10.16031/j.cnki.issn.1003-8035.202112028

基于改进时序InSAR技术的东莞地面沉降时空演变特征

基金项目: 广东省地勘事业发展基金项目(2019202)
详细信息
    作者简介:

    戴真印(1978-),男,安徽安庆人,本科,高级工程师,主要研究方向为工程地质与地质灾害。E-mail:zhydai@126.com

    通讯作者:

    刘岳霖(1990-),男,内蒙古赤峰人,硕士,工程师,主要从事地质环境调查评价、地质灾害防治工作。E-mail:liuyuelin@email.cugb.edu.cn

  • 中图分类号: P642.26

Spatial-temporal evolution characteristics of land subsidence in Dongguan City based on improved InSAR technology

  • 摘要: 东莞市是珠三角城市群和粤港澳大湾区的重要节点城市,深厚欠固结软弱土及其诱发的地面沉降已成为湾区内代表性的区域地质灾害问题,影响城市地质环境安全。为研究东莞市地面沉降发育规律及时空演变特征,采用改进时序InSAR技术对覆盖东莞地区的137景Sentinel-1 SLC SAR影像数据进行处理,分析了2015年6月至2020年6月地表形变动态演化规律。结果表明:(1)全域内地表沉降变形整体较稳定,沉降发育区占市域面积的34.6%,变形严重发育区主要集中在麻涌镇、道滘镇、洪梅镇、中堂镇、沙田镇及滨海湾新区;(2)大部分沉降变形点处于缓慢发展变形阶段,年平均沉降速率在20 mm/a以内,累计沉降量在1000 mm以内;(3)结合形变监测和现场调查,认为地面沉降与深厚软土发育和人类工程活动的耦合作用有很强的相关性。证明该方法能较好地识别和反映城市复杂形态区地面沉降发育的时空演化特征,为灾害预警、减避及治理提供技术支持。
    Abstract: Dongguan City is an important city of Pearl River Delta urban agglomeration and Guangdong-Hong Kong-Macao Greater Bay Area. Deep unconsolidated soft soil and its land subsidence have become a representative regional geological disaster in the Bay area, affecting the safety of urban geological environment. In order to study the development features and spatial-temporal evolution characteristics of land subsidence in Dongguan City, 137 sentinel-1 SLC SAR images covering whole Dongguan City were processed by improved InSAR technology, and the dynamic evolution characteristics of land deformation from June 2015 to June 2020 was analyzed. The results show that: (1) The land surface subsidence and deformation are stable in the whole region, and the subsidence developing areas account for 34.6% of the total urban area. The serious subsidence areas are mainly concentrated in Mayong Town, Daojiao Town, Hongmei Town, Zhongtang Town, Shatian Town and Binhai Bay New Area. (2) Most of the subsidence points are in the slow developing stage, the annual average subsidence rate is within 20 mm/yr, and the accumulated settlement is less than 1000 mm. (3) Combined with deformation monitor results and field investigation, land subsidence hazard has a great correlation with the coupling effect of deep soft soil development and human engineering activities. This method can better identify and reflect the temporal and spatial evolution characteristics of soft land subsidence development in urban area, and provide technical support for disaster early warning, mitigation and management.
  • 地裂缝为常见的地质灾害之一,一般指在内外应力协同作用下,岩土体的连续性被破坏,并发生变形破裂,进而在浅表层形成的具一定规模的变形带或裂缝的灾害现象,其变形破坏常呈现为水平剪切、垂直剪切与拉张破坏[13]。地裂缝在许多国家都有广泛发育,在美国、中国、墨西哥、澳大利亚以及埃塞俄比亚、肯尼亚等地均出现了大规模的地裂缝灾害,对当地居民生产生活造成严重影响,产生了难以计量的破坏和损失[410]

    调查发现,自上世纪60年代起,北京平原区勘察到地裂缝超过40条,至今已发育有7个地裂缝带,其中以高丽营、顺义和通州地区所发育的地裂缝最为典型,对周围建筑和公共设施造成了严重破坏,引起了人们的高度关注[11]。在北京地裂缝的发育特征和成因机理研究方面:贾三满等认为高丽营地裂缝主要受断裂的控制,为基底断裂活动延伸至地表导致的,是断裂构造和地面沉降共同作用的结果[12];彭建兵等通过对顺义区土沟村地裂缝进行研究,认为该地裂缝与下伏断裂相连,并呈现出同沉积断层的活动特征[13];刘德成等认为通州地裂缝是非构造性地裂缝,主要受地貌形态所控制,并分布于河流两岸,顺河堤方向展布[14];姜媛等在对大量资料分析的基础上,得到了高丽营地裂缝和当地地面沉降在时间、空间和活动方式上的关联,认为二者在加剧活动时间上具有同步性,在活动特征上具有相似性,并指出地面沉降在断裂两侧引起的不均匀沉降是高丽营地裂缝形成的主要因素[15];卢全中等对北七家至高丽营段的地裂缝进行研究,确定其上、下盘影响带总宽度为43 m,地裂缝对房屋与道路的破坏形式有3种[16];刘方翠等认为是区域构造活动及差异沉降变形机制的共同作用导致了北京平原区顺义地裂缝、高丽营地裂缝、羊房地裂缝和北小营地裂缝的快速形成发育[17];周永恒等通过地球物理探测技术与联合钻孔剖面方法揭示了黄庄-高丽营断裂(房山-涞水段)具典型的正断活动特征[18];关金环等发现北京首都国际机场发生的差异性沉降主要是由顺义-良乡断裂所发育的顺义地裂缝造成的[19];孟振江等以北京地裂缝为原型研究了耦合型地裂缝的发育活动特征和成因机制,揭示了不同位错量和水位下降量引发的地层位移场和应力场的变化特征[20];任雅哲等依据断层位错理论定量分析了顺义断裂蠕滑活动对机场地裂缝形成的影响,认为顺义断裂上盘中、深层地下水抽采引起的地面差异沉降是导致机场地裂缝加剧的主要因素[21]。以上研究在一定程度上揭示了北京地裂缝的成因与影响因素,为后续的深入研究提供了重要的研究基础和理论参考。

    作为特殊的工程场地,地裂缝场地对动力的响应明显且产生的影响较大。相关研究表明,傅里叶频率与场地条件、土层结构完整性及力学性质等紧密相关,而地脉动频率及卓越周期都可以反映出试验场地土体结构的变化[2223]。在利用地脉动测试方法研究方面:薛捷等通过地脉动测试的方法对西安地裂缝场地进行了研究,从而为西安市抗震设防工作提供理论依据[24];慕焕东等使用FLAC3D数值分析软件对西安地裂缝场地开展研究,发现地裂缝场地具有明显的放大效应,且上盘效应显著[25];张磊刚等应用地脉动测试手段对西安F6地裂缝场地进行研究,探明了场地因地裂缝而产生的放大效应和峰值变化规律[26];王晗等通过地脉动现场测试方法,得到了西安地区典型地裂缝的场地放大效应和变化规律[27];崔思颖等选取河北平原具有代表性的地裂缝进行地脉动测试揭示了动力响应特征[28];Chang J和Xuan Y等利用地脉动测试技术获得了山西临汾和大同两处盆地的典型地裂缝场地的数据并开展研究,研究结果表明可以通过地脉动信号来分析地裂缝场地的动力响应特征,并揭示了其变化规律与放大效应[2932]

    在对北京地区的地裂缝研究方面,前人所做研究大多集中于对地裂缝破坏现象的描述和产状勘察,以及从发育场地地貌和地质构造角度对其成因机理进行研究,在利用地脉动测试北京地裂缝发育场地的动力响应研究方面还属于空白,因该方法的应用还未被普及,加上北京多个地裂缝段现场的干扰较大,故尚未开展专门的地脉动现场测试[33]。鉴于此,本文的研究对象选取为北京通州区宋庄地裂缝场地,采用现场踏勘、地脉动测试和动力分析等方法,开展针对性的宋庄地裂缝场地的动力响应研究,以探究地裂缝的存在对场地动力响应的“峰值”及其“卓越频率”所造成的影响,并在此基础上揭示地裂缝发育场地的动力响应规律,以期为宋庄地区场地工程建设和防震减灾工程的设计和规划提供理论依据,同时为其他地裂缝发育区开展此类研究提供借鉴与参考。

    北京平原区地处太行山隆起与燕山隆起和华北断陷盆地相交汇的地带,自新生代以来,该地区主要处于伸展变形的构造环境,地质构造复杂且独特,断裂构造发育。近年来,随着人类工程活动的日益加剧,引发了多起地裂缝灾害[33]

    通州区地裂缝于上世纪60年代开始零星出现,唐山大地震发生后,郎府、西集、宋庄等地都出现了大量的地裂缝,其中以宋庄地裂缝最为典型。宋庄镇地处北京迭断陷与大兴迭隆起相交界的区域,其下伏具有活动迹象的早、中更新世南苑-通县断裂,该断裂带长约110 km,走向NE35~50°,倾向NW,倾角60°~75°,NW盘下降,SE盘上升,为正断层(图1)。断层西北侧震旦亚界埋深400~1000余米,上覆有巨厚的新近系、第四系沉积,而东南侧寒武系及中上元古界埋深仅60~300米,上覆有第四系及少量新近系沉积;对断层附近钻孔岩芯观察,具有明显挤压现象,应属压型断层。宋庄地裂缝于2010年后活动尤为剧烈[34]。经野外调查发现,宋庄地裂缝沿南苑-通县断裂展布,南起小中河在地表呈近直线延伸,经双埠头村、沟渠村、大庞村至平家疃村后向北东方向延伸,长度约8.7 km,整体走向NE31°,倾向NW。该地裂缝以水平拉张破坏为主,地表水平拉张量在2.15~56.7 cm之间,垂直位错量较小,地裂缝影响变形带最宽处达到400 m(图2)。

    图  1  宋庄地裂缝与周边断裂分布图
    Figure  1.  Distribution map of mainly fractures and location of Songzhuang ground fissures
    图  2  宋庄地裂缝致灾示意图
    Figure  2.  Schematic diagram of disaster caused by songzhuang ground fissures

    宋庄地裂缝的地层剖面如图3示,地表以下由粉砂、细砂、粉土、粉质黏土互层组成[35]。在浅表层,地裂缝近直立,地裂缝以偏拉张开裂为主,张开量为1~15 cm;在地下埋深50 m以下的地层开始逐渐出现位错现象,埋深53.2 m处的钻孔岩芯显示,在浅黑灰色砂质粘土层中,发育不连续构造面,构造面倾角64°,且两侧地层颜色不同。总体上分析发现地层上下盘位错量随埋深逐渐增大,地裂缝NW盘地层相对下降,西北侧第四系地层沉积厚度大于东南侧。在宋庄地裂缝影响范围内,致灾最为严重的是双埠头村至大庞村段,据统计,2016年、2018年和2022年地裂缝造成的直接经济损失总值分别为4623.24万元、7033.49万元和10751.14万元,可见地裂缝的影响和破坏在不断加剧[36]

    图  3  宋庄地裂缝浅表部地层剖面图(据赵龙,2018修改)
    Figure  3.  Stratigraphic section of the shallow strata of Songzhuang ground fissure (modified from Zhao Long, 2018)

    地脉动通常是指地表面振动周期在0.05~10s、振幅在千分之几至几微米的微幅振动,是由环境振动在地球表面所产生的复杂随机振动,包括地球内部应力、地震等自然因素及交通等人类活动所引起的振动,并经过不同传波介质和不同场地环境后综合作用于地球表面的结果。目前地脉动研究大致分为两类:一类是常时微动,一般是指周期小于1s的微动,侧重其振幅、周期性与观测点场地分类和振动特性研究;另一类是长波微动,一般是指周期大于1s的微动,研究的重点是由微动提取面波,进而解释推断观测点地下横波波速1s的常时微动,可用以研究地裂缝场地动力响应特性。

    地脉动测试要求的频段范围在0.1~15Hz之间,本次进行地脉动测试所使用的设备是高灵敏度伺服型速度网络地震仪,可满足测试精度要求,主要规格参数见表1

    表  1  选取地震仪的主要技术参数
    Table  1.  Main technical parameters of the Seismic Monitor
    仪器型号 CV-374AV
    动态范围 136 db
    灵敏度 1000 mv/(cm/s)
    通频带 0.1~100 Hz
    测量范围 ±0.02 m/s
    线性 0.03%
    采样频率 100 Hz
    AD分辨频率 24 bit
    下载: 导出CSV 
    | 显示表格

    高灵敏度伺服型速度网络地震仪尺寸为180*120*100 mm,灵敏度较高,在进行地脉动测试时,严格按照操作规程测试,确保测试数据的精确性。为避开不利环境因素和人类活动对地脉动监测结果的影响,选择天气条件良好情况下进行地脉动测试,测试工作选择在凌晨1点至5点期间。

    选择在北京通州区宋庄镇发育较为典型的双埠头村地裂缝段进行地脉动测试,现场布设测点具体原则为:在与主裂缝走向垂直的方向布设2条测线,结合地裂缝的影响带宽度范围,布设长度为60 m的测线,其中,上、下盘长度各设30 m。将地裂缝作为测线的中点,向远离地裂缝的方向进行测点的布设。以上盘为例,首先,在距地裂缝1.5 m处布设第一个测点,再沿测线以3 m为间距布设测点,共布设5个;接着按每隔5 m布设一个测点,共布设3个。因此,所布设测点与地裂缝间的水平距离由近到远设成1.5、3、6、9、12、15、20、25和30 m。下盘测点与上盘对称布设。

    针对地脉动观测各个测点的数据,在稳定波段内以10秒为单位进行波段的截取并通过Pwave32软件将其转换导出为加速度时程曲线,之后应用SeismoSignal软件对数据进行预处理,以得到反应谱曲线、傅里叶谱曲线以及Arias烈度曲线,通过频谱特征分析,进而获得地裂缝场地地脉动响应的加速度幅值、卓越频率、影响范围等场地动力响应特征;通过利用3种从不同角度出发的动力分析方法,对采集并经处理的地脉动数据进行分析,进一步提炼地脉动产生的场地动力特性。

    通过傅里叶谱分析,可以将地脉动加速度时程曲线转换成傅里叶谱曲线,从而得到场地卓越频率和傅里叶谱峰值,通过对不同测点卓越频率和峰值的对比,可得到场地不同区域地脉动响应的分布特征,从而获得地裂缝场地动力响应特征。

    图4图5为双埠头村地裂缝测线上盘和下盘在三个方向上的傅里叶谱图(其中X、Y、Z分别为垂直地裂缝走向方向、平行地裂缝走向方向及和垂直地面方向,下文同),图中各测点的傅里叶谱峰值显著,呈现出多峰的特点。主频带宽度较大,而各测点的卓越频率相对集中,因此其波速的变化没有明显的规律。傅里叶谱图反映出的多峰、谱面积大以及谐波丰富等特征,这与双埠头村地裂缝发育场地中较为松软的粘性土层结构特征相一致。通过对双埠头村地裂缝上下盘的傅里叶谱统计分析,得到该区域地裂缝场地的卓越频率在2.10~4.30 Hz之间浮动。而场地卓越频率的大小和测点距地裂缝的水平距离间的关系并不显著,均值约为3.30 Hz。

    图  4  双埠头村地裂缝上盘不同测点XYZ方向上傅里叶谱
    Figure  4.  Fourier spectra in the XYZ direction at different measurement points on the hanging wall of the ground fissures in Shuangbutou village
    图  5  双埠头村地裂缝下盘不同测点XYZ方向上傅里叶谱
    Figure  5.  Fourier spectra in the XYZ direction at different measurement points on the footwall of the ground fissures in Shuangbutou village

    图6图7为双埠头村地裂缝各测点的傅里叶卓越峰值图。从图中可以看出,在距离地裂缝区域较近的位置,地脉动响应由于地裂缝的存在而出现了显著的放大效应,傅里叶谱峰值与该距离呈负相关关系,距地裂缝的距离越大,峰值越低,在距离较远处曲线趋于平缓,其临界值约在距地裂缝12 m处,傅里叶谱峰值大致保持稳定不变。在距离地裂缝12 m以内的区域,振动的幅值和能量密度相对较高,而在距离地裂缝12米以外的区域,地表振动的幅值和能量密度则相对较低,这种现象表明靠近地裂缝的区域地表振动会比远离地裂缝的区域更强烈。

    图  6  双埠头村地裂缝上盘XYZ方向上傅里叶卓越峰值
    Figure  6.  The peak value of Fourier spectrum of the hanging wall in XYZ directions in Shuangbutou village
    图  7  双埠头村地裂缝下盘XYZ方向上傅里叶卓越峰值
    Figure  7.  The peak value of Fourier spectrum of the footwall in XYZ directions in Shuangbutou village

    图8图9分别为沟渠村地裂缝各测点X、Y、Z三个方向上傅里叶谱峰值与距地裂缝的距离变化关系的散点及曲线图,发现距离地裂缝越远的位置则傅里叶谱峰值越小,直至在较远的位置,曲线逐渐接近平缓;最终在距地裂缝约12 m处,曲线基本保持不变,傅里叶谱峰值趋于一致。通过统计分析,得出沟渠村地裂缝的场地卓越频率在2.30~4.30 Hz范围内变化,卓越频率的高低与测点所在的位置之间没有明显关系,其均值约为3.40 Hz。

    图  8  沟渠村地裂缝上盘XYZ方向上傅里叶卓越峰值
    Figure  8.  The peak value of Fourier spectrum of the hanging wall in XYZ directions in Gouqu village
    图  9  沟渠村地裂缝下盘XYZ方向上傅里叶卓越峰值
    Figure  9.  The peak value of Fourier spectrum of the footwall in XYZ directions in Gouqu village

    通过反应谱分析,可以得到特征周期和相应的反应谱峰值,通过对不同测点特征周期和反应谱峰值的对比,可以得到场地不同区域地脉动响应的分布特征,从而获得地裂缝场地动力响应特征。

    经数据分析得到,双埠头村地裂缝各测点在同一方向上的反应谱谱型均为“单峰型”,且具有主峰突出、频带窄和谱面积小的特点;反应谱卓越周期集中在0.1 s~0.2 s之间,部分卓越周期可达0.3s,此外卓越周期值与地裂缝和测点间的水平距离无明显关系(图1011)。

    图  10  双埠头村地裂缝上盘不同测点XYZ方向上反应谱
    Figure  10.  Response spectrum in the XYZ direction at different measurement points on the hanging wall of the ground fissures in Shuangbutou Village
    图  11  双埠头村地裂缝下盘不同测点XYZ方向上反应谱
    Figure  11.  Response spectrum in the XYZ direction at different measurement points on the footwall of the ground fissures in Shuangbutou Village

    图12图13可得,在距离地裂缝较近处,因地裂缝的存在地脉动响应表现出显著的放大效应,随着与地裂缝距离的增大,地脉动的响应加速度基本不变,在距地裂缝12m处趋于平稳,反应谱峰值基本一致。

    图  12  双埠头村地裂缝上盘XYZ方向上反应谱加速度峰值变化曲线
    Figure  12.  Response spectrum acceleration peak value variation curves in XYZ direction on the hanging wall of the ground fissures in Shuangboutou village
    图  13  双埠头村地裂缝下盘XYZ方向上反应谱加速度峰值变化曲线
    Figure  13.  Variation curves of peak values of reaction spectral acceleration in XYZ direction on the footwall of the ground fissures in Shuangbutou Village

    图14图15可以看出,沟渠村地裂缝各测点在X、Y、Z三个方向的上反应谱加速峰值随与地裂缝间的垂直距离变化,距地裂缝较近的区域的地脉动响应表现出显著的放大效应,且距离地裂缝越远,反应谱峰值越小,直至曲线接近稳定不变,反应谱峰值基本一致。

    图  14  沟渠村地裂缝上盘XYZ方向上反应谱加速度峰值变化曲线
    Figure  14.  Variation Curve of Peak Acceleration in the XYZ Directions of the Response Spectrum on the Hanging Wall of the Ground Fissure in Gouqu Village
    图  15  沟渠村地裂缝下盘XYZ方向上反应谱加速度峰值变化曲线
    Figure  15.  Variation Curve of Peak Acceleration in the XYZ Directions of the Response Spectrum on the Footwall of the Ground Fissure in Gouqu Village

    通过对比双埠头村地裂缝和沟渠村地裂缝反应谱曲线,发现布设在宋庄地裂缝上下盘的测点,在与地裂缝同等距离时,上盘反应谱的峰值均大于下盘,这表明地裂缝的存在使得地裂缝两侧场地的地脉动响应发生了差异变化,存在显著“上盘效应”,且上盘峰值大于下盘。

    Arias烈度是一种目前较为成熟的评判场地地震动响应强弱、反映场地抗震能力和设防标准的方法,通过分析Arias烈度曲线的变化特征进而了解地裂缝场地不同位置的抗震特性,可从场地抗震烈度方面揭示地裂缝对场地的放大效应。

    图16图17分别为双埠头村地裂缝上盘与下盘不同测点的Arias烈度曲线图,图中注释栏“MP”表示监测点。结果表明,Arias烈度呈现出随时间逐渐增长的趋势,但各个测点的增长幅度有所不同,且有着不同的最大值。

    图  16  双埠头村地裂缝上盘不同测点XYZ方向上Arias烈度曲线
    Figure  16.  Arias intensity profiles in the XYZ direction at different measurement points on the hanging wall of the ground fissure in Shuangbutou village
    图  17  双埠头村地裂缝下盘不同测点XYZ方向上Arias烈度曲线
    Figure  17.  Arias intensity profiles in the XYZ direction at different measurement points on the footwall of the ground fissure in Shuangbutou village

    为反映地裂缝对场地不同位置烈度值产生的影响,将双埠头村地裂缝各测点Arias烈度的峰值沿垂直地裂缝走向方向展开,得到上下盘各测点在X、Y、Z三个方向上的Arias烈度峰值随与地裂缝距离变化的散点数据曲线图(图1819)。可见,场地不同位置的烈度值受地裂缝的影响程度也不同。Arias烈度在距离地裂缝较近的区域,其峰值较大,且该烈度值随测点与地裂缝间的距离增大而减小并逐渐趋于稳定,其临界距离为12 m,而临界区域以外的范围受放大效应的影响较小。

    图  18  双埠头村地裂缝上盘XYZ方向上Arias烈度峰值随距离变化曲线
    Figure  18.  Variation curves of peak Arias intensity with distance in the XYZ direction on the hanging wall of the ground fissures in Shuangbutou village
    图  19  双埠头村地裂缝下盘XYZ方向上Arias烈度峰值随距离变化曲线
    Figure  19.  Variation curves of peak Arias intensity with distance in the XYZ direction on the footwall of the ground fissures in Shuangbutou village

    (1)宋庄地裂缝延伸8.7 km,地表破坏以拉张开裂为主,水平拉张量在2.15~56.7 cm之间;其中双埠头村地裂缝的傅里叶卓越频率主要集中在2.10~4.30 Hz之间,反应谱的卓越周期集中在0.1~0.2 s之间;沟渠村地裂缝的傅里叶卓越频率集中在2.30~4.30 Hz之间,反映了地裂缝不同区段活动性的变化。

    (2)通过反应谱分析法、傅里叶谱分析法分析得到,测点的卓越频率无明显的规律性变化;Arias烈度峰值在靠近地裂缝处较大,且随着距离的增大而逐渐减小,在距地裂缝较远的区域趋于平缓。

    (3)宋庄地裂缝场地动力响应在反应谱、傅里叶谱和Arias烈度上都存在明显的“放大效应”,即靠近地裂缝的区域大于远离地裂缝的区域,且该临界距离为距地裂缝12 m。

    (4)宋庄地裂缝场地有着较为显著的“上盘效应”,即距地裂缝相同距离情况下,位于上盘的傅里叶卓越频率、反应谱卓越周期和Arias烈度峰值总大于下盘。

  • 图  1   研究区位置及范围

    Figure  1.   Location of the study area

    图  2   研究区及影像覆盖范围

    Figure  2.   Study area and image coverage

    图  3   干涉对小基线集网络连接图

    Figure  3.   Small baseline set network connection diagram of interference pair

    图  4   三次线性回归后的高程残差示意图

    Figure  4.   Schematic diagram of elevation residual after three linear regression

    图  5   改进的时序InSAR技术流程图

    Figure  5.   Flow chart of improved time-series InSAR

    图  6   东莞市2015—2020年间各年年平均形变速率

    Figure  6.   Annual average subsiding rate between 2015 and 2020

    图  7   东莞市主要沉降变形区分布

    Figure  7.   The distribution of serious subsidence area

    图  8   长安镇典型监测区的时间序列形变特征

    Figure  8.   Subsiding characteristics of time-series in typical monitoring area of Chang’an Town

    图  9   长安镇调查点地面沉降变形特征

    Figure  9.   Characteristics of land subsidence deformation in Chang’an Town

    图  10   重点区软土厚度与InSAR解译累计沉降量对比图

    Figure  10.   Comparison of soft soil thickness and InSAR interpretation accumulated settlement in main investigation areas

    表  1   Sentinel-1卫星数据参数表

    Table  1   Parameters of Sentinel-1 satellite

    Sentinel-1参数值
    轨道高度/km700
    重返周期/d12
    入射角/(°)29~46
    分辨率5 m×20 m
    幅宽/km250
    极化方式HH+HV、VH+VV、HH、VV
    影像时间2015年6月15日—2020年6月12日
    影像数量/景137
    下载: 导出CSV

    表  2   东莞市沉降速率统计表

    Table  2   Statistics data of annual average subsiding rate

    发育程度年均形变速率/(mm·a−1分布面积/km2
    2015年2016年2017年2018年2019年2020年平均值
    弱发育0~−102444.272453.042444.702409.442413.012429.432432.31
    中等发育−20~−1013.706.4111.7645.3542.8326.5924.44
    −30~−201.460.532.474.173.082.662.40
    强发育−40~−300.510.110.800.800.790.820.64
    −40及以上0.160.010.370.340.390.600.31
    下载: 导出CSV
  • [1] 葛伟丽,李元杰,张春明,等. 基于InSAR技术的内蒙古巴彦淖尔市地面沉降演化特征及成因分析[J]. 水文地质工程地质,2022,49(4):198 − 206. [GE Weili,LI Yuanjie,ZHANG Chunming,et al. An attribution analysis of land subsidence features in the city of Bayannur in InnerMongolia based on InSAR[J]. Hydrogeology & Engineering Geology,2022,49(4):198 − 206. (in Chinese with English abstract)
    [2] 王双,严学新,揭江,等. 珠三角平原区软土分布与地面沉降相关性分析[J]. 上海国土资源,2019,40(2):75 − 79. [WANG Shuang,YAN Xuexin,JIE Jiang,et al. Correlation analysis between soft soil distribution and land subsidence in the Pearl River Delta Plain[J]. Shanghai Land & Resources,2019,40(2):75 − 79. (in Chinese with English abstract) DOI: 10.3969/j.issn.2095-1329.2019.02.015
    [3] 姚佳明,姚鑫,陈剑,等. 基于InSAR技术的缓倾煤层开采诱发顺层岩体地表变形模式研究[J]. 水文地质工程地质,2020,47(3):135 − 146. [YAO Jiaming,YAO Xin,CHEN Jian,et al. A study of deformation mode and formation mechanism of a bedding landslide induced by mining of gently inclined coal seam based on InSAR technology[J]. Hydrogeology & Engineering Geology,2020,47(3):135 − 146. (in Chinese with English abstract)
    [4] 张凯翔. 基于“3S”技术的地质灾害监测预警系统在我国应用现状[J]. 中国地质灾害与防治学报,2020,31(6):1 − 11. [ZHANG Kaixiang. Review on geological disaster monitoring and early warning system based on “3S” technology in China[J]. The Chinese Journal of Geological Hazard and Control,2020,31(6):1 − 11. (in Chinese with English abstract) DOI: 10.16031/j.cnki.issn.1003-8035.2020.06.01
    [5] 潘光永,陶秋香,陈洋,等. 基于SBAS-InSAR的山东济阳矿区沉降监测与分析[J]. 中国地质灾害与防治学报,2020,31(4):100 − 106. [PAN Guangyong,TAO Qiuxiang,CHEN Yang,et al. Monitoring and analysis of sedimentation in Jiyang mining area of Shandong Province based on SBAS-InSAR[J]. The Chinese Journal of Geological Hazard and Control,2020,31(4):100 − 106. (in Chinese with English abstract)
    [6]

    FERRETTI A,PRATI C,ROCCA F. Nonlinear subsidence rate estimation using permanent scatterers in differential SAR interferometry[J]. IEEE Transactions on Geoscience and Remote Sensing,2000,38(5):2202 − 2212. DOI: 10.1109/36.868878

    [7]

    FERRETTI A, PRATI C, ROCCA F L. Permanent scatterers in SAR interferometry[C]//Remote Sensing. Proc SPIE 3869, SAR Image Analysis, Modeling, and Techniques II, Florence, Italy. 1999, 3869: 139 − 145.

    [8]

    BERARDINO P,FORNARO G,LANARI R,et al. A new algorithm for surface deformation monitoring based on small baseline differential SAR interferograms[J]. IEEE Transactions on Geoscience and Remote Sensing,2002,40(11):2375 − 2383. DOI: 10.1109/TGRS.2002.803792

    [9] 姜兆英,于胜文,陶秋香. StaMPS-MTI技术在地面沉降监测中的应用[J]. 西南交通大学学报,2017,52(2):295 − 302. [JIANG Zhaoying,YU Shengwen,TAO Qiuxiang. Application of StaMPS-MTI technology in monitoring ground subsidence[J]. Journal of Southwest Jiaotong University,2017,52(2):295 − 302. (in Chinese with English abstract) DOI: 10.3969/j.issn.0258-2724.2017.02.012
    [10] 孔祥如,罗勇,刘贺,等. PS-InSAR技术在北京通州区地面沉降监测中的应用[J]. 城市地质,2021,16(1):25 − 31. [KONG Xiangru,LUO Yong,LIU He,et al. Application of PS-InSAR technology in the land subsidence survey in Tongzhou District,Beijing[J]. Urban Geology,2021,16(1):25 − 31. (in Chinese with English abstract) DOI: 10.3969/j.issn.1007-1903.2021.01.004
    [11] 晏霞,刘媛媛,赵振宇. 利用时序InSAR技术监测南水进京后北京平原地区的地面沉降[J]. 地球物理学进展,2021,36(6):2351 − 2361. [YAN Xia,LIU Yuanyuan,ZHAO Zhenyu. Land subsidence monitoring after the start of the South to North Water Transfer in the Beijing plain based on multi-temporal InSAR[J]. Progress in Geophysics,2021,36(6):2351 − 2361. (in Chinese with English abstract)
    [12] 狄桂栓. 基于InSAR技术的黄河三角洲区域地表形变浅析[J]. 地理空间信息,2020,18(9):106 − 109. [DI Guishuan. Primary investigation of surface deformation in the Yellow River Delta based on InSAR technology[J]. Geospatial Information,2020,18(9):106 − 109. (in Chinese with English abstract) DOI: 10.3969/j.issn.1672-4623.2020.09.027
    [13] 张宏雪,曾润强,孙萍萍,等. 基于SBAS-InSAR的延安新区地表沉降监测与分析[J]. 兰州大学学报(自然科学版),2021,57(2):185 − 193. [ZHANG Hongxue,ZENG Runqiang,SUN Pingping,et al. Surface subsidence monitoring and analysis in Yan'an New District based on SBAS-InSAR[J]. Journal of Lanzhou University (Natural Sciences),2021,57(2):185 − 193. (in Chinese with English abstract)
    [14] 蒲川豪,许强,蒋亚楠,等. 延安新区地面沉降分布及影响因素的时序InSAR监测分析[J]. 武汉大学学报(信息科学版),2020,45(11):1728 − 1738. [PU Chuanhao,XU Qiang,JIANG Yanan,et al. Analysis of land subsidence distribution and influencing factors in Yan'an new district based on time series InSAR[J]. Geomatics and Information Science of Wuhan University,2020,45(11):1728 − 1738. (in Chinese with English abstract) DOI: 10.13203/j.whugis20190372
    [15] 冉培廉,李少达,杨晓霞,等. 基于SBAS-InSAR技术的西安市地面沉降监测[J]. 河南理工大学学报(自然科学版),2021,40(3):66 − 74. [RAN Peilian,LI Shaoda,YANG Xiaoxia,et al. Monitoring of Xi'an City land subsidence based on SBAS-InSAR[J]. Journal of Henan Polytechnic University (Natural Science),2021,40(3):66 − 74. (in Chinese with English abstract)
    [16] 范军,左小清,李涛,等. PS-InSAR和SBAS-InSAR技术对昆明主城区地面沉降监测的对比分析[J]. 测绘工程,2018,27(6):50 − 58. [FAN Jun,ZUO Xiaoqing,LI Tao,et al. Analysis and comparison of PS-InSAR and SBAS-InSAR for ground subsidence monitoring in the main city of Kunming[J]. Engineering of Surveying and Mapping,2018,27(6):50 − 58. (in Chinese with English abstract) DOI: 10.19349/j.cnki.issn1006-7949.2018.06.010
    [17] 周书东, 张益. 东莞市工程地质特征及常见工程问题[J]. 现代隧道技术, 2020, 57(增刊1): 78 − 88

    ZHOU Shudong, ZHANG Yi. Engineering geological characteristics and common engineering issues in Dongguan City[J]. Modern Tunnelling Technology, 2020, 57(Sup 1): 78 − 88. (in Chinese with English abstract)

    [18] 黄建龙,刘亦农,曾伟国. 粤港澳大湾区地质特点与地质环境保护策略分析[J]. 人民珠江,2019,40(9):103 − 109. [HUANG Janlong,LIU Yinong,ZENG Weiguo. Analysis of geological characteristics and geological environment protection in Guangdong-Hong Kong-Macao Greater Bay Area[J]. Pearl River,2019,40(9):103 − 109. (in Chinese with English abstract) DOI: 10.3969/j.issn.1001-9235.2019.09.015
    [19] 杨魁,杨建兵,江冰茹. Sentinel-1卫星综述[J]. 城市勘测,2015(2):24 − 27. [YANG Kui,YANG Jianbing,JIANG Bingru. Sentinel-1 satellite overview[J]. Urban Geotechnical Investigation & Surveying,2015(2):24 − 27. (in Chinese with English abstract) DOI: 10.3969/j.issn.1672-8262.2015.02.006
    [20] 刘琦,岳国森,丁孝兵,等. 佛山地铁沿线时序InSAR形变时空特征分析[J]. 武汉大学学报(信息科学版),2019,44(7):1099 − 1106. [LIU Qi,YUE Guosen,DING Xiaobing,et al. Temporal and spatial characteristics analysis of deformation along Foshan subway using time series InSAR[J]. Geomatics and Information Science of Wuhan University,2019,44(7):1099 − 1106. (in Chinese with English abstract) DOI: 10.13203/j.whugis20190025
    [21] 伍素贞,谢荣安,谢文珠,等. 时序InSAR 2018—2019年广州市南沙区形变监测与分析[J]. 工程勘察,2020,48(8):48 − 52. [WU Suzhen,XIE Rongan,XIE Wenzhu,et al. Ground subsidence monitoring in Nansha District by using Sentinel 1A/B SAR images[J]. Geotechnical Investigation & Surveying,2020,48(8):48 − 52. (in Chinese with English abstract)
    [22] 梁景才. 雷州半岛地区InSAR时序形变监测与分析[J]. 地矿测绘,2020,36(4):21 − 26. [LIANG Jingcai. InSAR time series deformation monitoring and analysis in Leizhou peninsula[J]. Surveying and Mapping of Geology and Mineral Resources,2020,36(4):21 − 26. (in Chinese with English abstract) DOI: 10.3969/j.issn.1007-9394.2020.04.006
    [23] 谢思梅,吴曼乔,钟惠润. 短基线集技术(SBAS)在南沙自贸区地面沉降监测中的应用研究[J]. 信息技术与信息化,2020(8):136 − 138. [XIE Simei,WU Manqiao,ZHONG Huirun. Application of short baseline set technique (SBAS) in land subsidence monitoring in Nansha free trade zone[J]. Information Technology and Informatization,2020(8):136 − 138. (in Chinese) DOI: 10.3969/j.issn.1672-9528.2020.08.044
    [24] 贾会会,张海清,李克达,等. 融合分布式散射体时序InSAR技术在矿区形变调查中的应用[J]. 吉林大学学报(地球科学版),2022,52(1):202 − 213. [JIA Huihui, ZHANG Haiqing, LI Keda, et al. Application of fusion distributed scatterer time-series InSAR technique in deformation survey of mining area[J]. Journal of Jilin University (Earth Science Edition),2022,52(1):202 − 213. (in Chinese with English abstract)
图(10)  /  表(2)
计量
  • 文章访问数:  3161
  • HTML全文浏览量:  1225
  • PDF下载量:  408
  • 被引次数: 0
出版历程
  • 收稿日期:  2021-12-23
  • 修回日期:  2022-05-29
  • 网络出版日期:  2022-11-29
  • 刊出日期:  2023-02-24

目录

/

返回文章
返回