陕西韩城台巨幅地倾斜异常的成因分析

杨小林, 王军, 王希彬, 王新, 窦玛丽, 危自根

杨小林,王军,王希彬,王新,窦玛丽,危自根. 2018. 陕西韩城台巨幅地倾斜异常的成因分析. 地震学报,40(6):760−773. doi:10.11939/jass.20180025. DOI: 10.11939/jass.20180025
引用本文: 杨小林,王军,王希彬,王新,窦玛丽,危自根. 2018. 陕西韩城台巨幅地倾斜异常的成因分析. 地震学报,40(6):760−773. doi:10.11939/jass.20180025. DOI: 10.11939/jass.20180025
Yang X L,Wang J,Wang X B,Wang X,Dou M L,Wei Z G. 2018. What causes the remarkable tilt anomalies at the Hancheng geodynamic observatory in Shaanxi Province? Acta Seismologica Sinica40(6):760−773. doi:10.11939/jass.20180025. DOI: 10.11939/jass.20180025
Citation: Yang X L,Wang J,Wang X B,Wang X,Dou M L,Wei Z G. 2018. What causes the remarkable tilt anomalies at the Hancheng geodynamic observatory in Shaanxi Province? Acta Seismologica Sinica40(6):760−773. doi:10.11939/jass.20180025. DOI: 10.11939/jass.20180025

陕西韩城台巨幅地倾斜异常的成因分析

基金项目: 国家重点研发计划“深地资源勘查开采”重点专项(2016YFC0600402)和中国地震局2017年度震情跟踪课题(2017010206)联合资助
详细信息
    通讯作者:

    杨小林: e-mail:yang-xiaolin123@163.com

  • 中图分类号: P315.72+5

What causes the remarkable tilt anomalies at the Hancheng geodynamic observatory in Shaanxi Province?

  • 摘要: 巨幅地倾斜异常既可能是地震前兆异常信息,也有可能是仪器问题或环境干扰所致的异常信号. 有效地厘清其性质,对地震前兆异常的及时识别与科学判定至关重要. 陕西韩城台金属水平摆EW分量自2010年以来连续两次出现巨幅东倾和西倾异常,幅度分别约达140″和180″,该巨幅地倾斜异常产生的根源至今尚未被厘清. 鉴于此,本研究依据韩城台所在区域的水文、构造和地震活动等特征,提出并分析了地下水动力变化、韩城断裂慢滑移和地壳应力场变动等3种可能的成因机制. 结果表明,第一种成因机制难以有效地解释巨幅地倾斜异常,第二和第三种成因机制则具有一定的可能性,但证据还不够充分. 因此,更可靠的物理解释尚需更多的观测和更深入的研究. 尽管本文未能给出该巨幅异常的真正成因,但所采取的分析方法可为今后巨幅地倾斜异常性质的判定工作提供有益的参考.
    Abstract: The remarkable tilt anomalies could be the earthquake precursors, but may also be caused by instrumental factors and environmental disturbances. Thus, the question arises on how to distinguish the earthquake precursors from the non-tectonic factors, which is very important to effectively and reasonably detect earthquake precursors. Since 2010, two remarkable tilt anomalies have been recorded by metallic horizontal pendulums in E-W component at the Hancheng observatory in Shannxi Province, and the amount of east- and west-ward tilt approxi-mately reach up to 140″ and 180″, respectively, but these two remarkable tilt anomalies have not been reasonably and clearly interpreted till now. Here, we propose and compare three different causal mechanisms possibly responsible for these anomalous phenomena according to the regional hydrological, tectonic and seismicity characteristics, i.e. ① hydrodynamics-induced surface tilt, ② a long-term slow slip event on the northeastern segment of the Hancheng fault, and ③ variations of the regional tectonic stress field during the anomalous period. We then theoretically calculated the poroelastic deformation and the fault slip amount, and finally inversed the focal mechanism solutions of 85 earthquakes (2.0≤ML≤4.8) that occurred between 2008 and 2015 with the aim of determining the regional stress field changes (35°N—36°N, 110°E—111°E) in the crust. Our results show that the first possibility can be shown unlikely, but it is difficult to rule out the second and the third possibility according to the current eviden-ce. To further prove and confirm the causal relationship between deformation of tectonic origin and the anomalies, more comprehensive tilt and crustal deformation measurements are necessary in the Hancheng region in the future, furthermore, more intensive researches are also needed to reveal and determine the causal mechanisms of these anomalies. Unfortunately, we fail to find the real causal mechanism, but the approaches used in this study could be helpful to investigate the causal origin of remarkable anomalies recorded by tiltmeters in the near future.
  • 强震除了会造成房屋倒塌、道路中断等直接损失,还会引起火灾、水灾、滑坡、泥石流等次生灾害.空间对地观测技术的发展为地震调查提供了经济、快速、安全的评估方式 (李德仁,2009):在震前预案阶段,遥感图像为风险评估提供参考;在灾后应急阶段,遥感图像用于震害灾情快速评估和地震次生灾害的动态过程监测;在损失评估和重建阶段,遥感图像辅助震害建筑物、道路等的灾情损失调查、宏观烈度的快速评估以及灾区重建选址等工作.尤其是合成孔径雷达 (synthetic aperture radar,简写为SAR) 和激光雷达 (light detection and ranging,简写为LiDAR) 等新型遥感技术的不断发展,为震后灾情快速评估与应急救援工作增添了新的技术途径.合成孔径雷达由于不受气候、天气情况影响,能够记录地物侧面、表面粗糙度和角反射信息,而激光雷达能够获取震害地物目标表面的高精度三维坐标.空间探测技术的不断进步,使得遥感朝多空间分辨率、多时间分辨率和多光谱分辨率的方向发展.遥感影像资源的多样性丰富了遥感震害评估的手段,但同时也增加了地物目标震害识别方法的复杂度.需综合分析目标震害遥感机理,才能有针对性地选择遥感数据种类和处理方法,提高震害评估的效率与精度.

    震害遥感特征主要指震害在遥感图像上的表现现象.来源于不同平台、不同传感器上的图像特征差异很大,这是由于光谱反射、粗糙度和散射特性等不同所致.遥感图像特征是以典型图像的解剖、原理图示、表格分析等方式对图像的色调、图形、纹理、地物关系等进行描述,这方面的研究成果已有很多. Brunner等 (2010)在分析高分辨率光学图像和SAR图像中震害建筑物特征的基础上,通过两者特征的相互融合,实现了震害建筑物的评估;赵福军 (2010)研究了光学图像中的震害建筑物特征,并以对象为单元构建了震害特征规则集;刘金玉等 (2013)从SAR图像成像机理和目标的后向散射特征出发,细致分析了建筑物的震害特征,为寻求目标检测方法提供了思路;焦其松 (2015)以北川地震地面LiDAR数据为例,详细分析了不同震害类型建筑物的三维特征表现,并提出了基于形状分析的震害建筑物提取方法.前人的研究主要集中于单一遥感数据源的特征分析或两种数据源的特征融合,且针对特定的研究区域,所采用的特征集合并未完全覆盖多源遥感图像中的特征,应用的特征具有较大的局限性,因此需细致梳理多源遥感图像中不同震害建筑物类型的特征表现,掌握震害地物目标在图像中的成像机理与特征,为构建震害自动化提取流程和分析方法提供参考.

    鉴于此,本文拟在介绍震害遥感分类分级的基础上,重点细致地介绍震害建构筑物在光学、合成孔径雷达和激光雷达图像中各自的表现特征,并联合多源遥感影像综合对比分析建筑物震害特征,从而达到梳理震害提取特征的目的,为地震应急与评估过程中不同源遥感图像的表现优势与应用方法提供参考.

    地震发生之后,由于建筑物的外观结构、建造材料等多有差异,其破坏形式也各不相同.长期以来,地震震害分类分级主要通过人工实地勘测完成,目前已形成一套完整的震害破坏等级判定体系.自20世纪60年代以来,遥感技术由于其具有宏观、综合等特点而被引入震害信息识别中,经过数十年的研究和应用,也已形成了一套初步的遥感震害分类分级方案.

    按照地震现场工作的调查规范 (中华人民共和国国家质量监督检验检疫总局, 中国国家标准化管理委员会,2011),建筑物震害可划分为5个破坏级别:基本完好,轻微破坏,中等破坏,严重破坏,毁坏.宏观地面调查能够比较全面地考察建筑物的承重结构类型、承重构件 (墙体等) 是否出现裂缝、是否被破坏等情况,而在遥感影像上判断建筑物的震害等级主要是根据建筑物顶部的破坏情况.遥感震害调查一般将建筑物震害类型分为基本完好、中等破坏、毁坏等3类 (王栋梁,2007张磊,2008).宏观地面调查分类与遥感震害调查分类的对应关系列于表 1.

    表  1  宏观地面调查与遥感震害调查的建筑物震害类型对应表
    Table  1.  The correspondence of seismic damage types between macroscopic ground investigation and remote sensing investigation
    序号 宏观地面调查分类 遥感震害调查分类
    1 基本完好
    2 轻微破坏 基本完好
    3 中等破坏
    4 严重破坏 中等破坏
    5 毁坏 毁坏
    下载: 导出CSV 
    | 显示表格

    在遥感震害评估工作中,多源遥感数据源主要为光学影像、SAR图像和LiDAR图像,震害建筑物目标在不同源遥感图像上所反映的特征差异较大.为有效地综合利用这些特征,需对多源遥感图像震害特征进行梳理分析.本节将按照多源遥感图像类型,重点分析不同震害程度建筑物的图像特征.

    为了能够快速准确地检测出建筑物的倒塌损毁情况,需要获取光学遥感影像上建筑物及其震害的图像特征,并将文字描述性特征转化为特征参量规则集,有针对性地选择监督分类或面向对象分析等方法来实现震害建筑物的识别 (Mitomi et al,2000Turker,San,2004Zhu et al,2008).

    图 1给出了基本完好、中等破坏、毁坏等3类震害等级的建筑物影像,主要从灰度特征、形状特征、纹理特征、上下文特征等方面对不同等级震害在影像中的特征进行对比分析.

    图  1  汶川地震都江堰地区震后光学影像
    (a) 基本完好建筑物;(b) 中等破坏建筑物;(c) 毁坏建筑物
    Figure  1.  Optical images of Dujiangyan area after the Wenchuan earthquake
    (a) Intact building; (b) Moderately damaged building; (c) Destroyed building

    1) 灰度特征.基本完好建筑物表面的灰度分布比较均匀,建筑物顶部亮度差异较大;建筑物倒塌后,因其顶部破碎使灰度分布变得复杂 (柳稼航,2003).

    2) 形状特征.基本完好建筑物的几何形状规整,以直角矩形为主; 倒塌建筑物顶部出现断裂或几何形状不完整,或者由于建筑物顶部部分倒塌使得边界模糊不清,或者出现大的断裂或者部分坍落、边角坠落现象.

    3) 纹理特征.基本完好建筑物的顶部具有明显的纹理特征,纹理通常与边界线有一致的方向性, 部分大型厂房和学校楼房的顶部有明暗相间的纹理;建筑物的倒塌,破坏了原来影像中建筑物在结构和空间布局上的规律性,在影像上表现为轮廓消失、图案零乱、纹理不规则等特征 (赵福军,2010).

    4) 上下文特征.建筑物周围的主要地物是街道两旁的树木和街道等; 建筑物毁坏后,建筑物周边的道路上布满大量的瓦砾或墙体碎屑,造成道路的堵塞或隔断.

    SAR通过主动发射微波脉冲,接收目标物体对雷达波束的后向散射信号,所接收的信号同时以振幅和相位的形式保存 (张红等,2009).成像的SAR图像包含了丰富的地物目标信号,这些信号能够较好地体现成像目标的多维特征 (Matsuoka,Yamazaki, 2002, 2004).由于震害地物自身结构比较复杂,SAR图像丰富的信息量为研究震害目标的成像机理提供了有利的条件.下面主要从基本完好建筑物、中等破坏建筑物、毁坏建筑物等3类震害建筑物类型出发,分析其在SAR图像中的特征及机理.

    在中低分辨率图像中,单个像素所反映的是单个或多个建筑及其部件、背景的回波综合叠加结果,仅从图像中很难辨别单个建筑的结构和状态.建筑物在图像上一般以区域的形式进行描述.

    在高分辨率图像中,所呈现的是建筑物的外形及表面结构.规则分布的强散射点、角反射亮线、建筑阴影暗区构成了高分辨率SAR图像中建筑物的特征 (邵芸等,2008).一般来说,完好的建筑在SAR图像上包含叠掩区、角反射区、屋顶区和阴影区等部分,如图 2所示.从雷达图像特征上来看,城区完好的建筑物呈现比较规则的排列,彼此之间的空间关系符合楼群特征.由于SAR的侧视成像特点,在近距,迎入射方向,墙面的叠掩区强度通常高于地表的叠掩区强度;稍远处,由于墙体与地表的二面角反射效应使得房屋一侧的后向散射强度高而呈现“1”或“L”形强回波特征,强回波位置仍然符合房屋的排列关系;再远处,屋顶较弱的单次反射;在背向雷达入射方向,一般可以看到较清晰的矩形状阴影区 (Gamba et al,2007).

    图  2  完好建筑物SAR图像成像示例
    A为叠掩区,B为角反射区,C为屋顶反射区,D为阴影区,θ为入射角
    Figure  2.  SAR imaging of intact buildings
    A represents the layover area, B represents the angular reflection area, C represents the roof reflection area, D represents the shadow area, and θ represents incident angle

    图 3为不同分辨率SAR图像中的建筑 (区).在分辨率为20 m的ENVISAT ASAR图像 (图 3a) 中,建筑表现为亮点.由于图像分辨率较低,在建筑密集的中心城区,建筑的后向散射基本混合在一起,成为亮区,很难分辨单个建筑.在分辨率为8 m的RADARSAT-1图像 (图 3b) 中,可见整齐排列的建筑.在分辨率为3 m的TerraSAR-X图像 (图 3c) 中,单个建筑已可分辨,且部分建筑的叠掩、二次散射、阴影等特征已可见.在分辨率为1 m的COSMO-SkyMed SAR图像 (图 3d) 中,单个建筑的图像更为清晰,甚至可分辨其部分结构,例如由房顶的建筑结构所形成的亮线/点以及建筑的进一步细节信息 (如窗台角反射产生的亮点) 也可显示出来.

    图  3  不同分辨率SAR图像中的建筑 (区)
    (a) 分辨率为20 m的ENVISAT ASAR图像; (b) 分辨率为8 m的RADARSAT-1图像; (c) 分辨率为3 m的TerraSAR-X图像;(d) 分辨率为1 m的COSMO-SkyMed图像
    Figure  3.  The building samples in different resolution SAR images
    (a) 20 m-resolution ENVISAT ASAR image; (b) 8 m-resolution RADARSAT-1 image; (c) 3 m-resolution TerraSAR-X image; (d) 1 m-resolution COSMO-SkyMed image

    SAR图像中基本完好建筑物的主要特征为 (安立强,2011):①灰度分布特征.朝向雷达传感器一侧的建筑物墙面与地面互相垂直形成二面角反射器,回波较强,形成一条明显的亮条纹,叠掩效应、二面角反射等共同形成建筑物的主线条,使得建筑物在图像上表现为高亮度,建筑物的侧面及水平屋顶导致图像上形成黑色区;②几何特征.建筑物具有规则的几何外形,通常是长方形或“L”形;③高分辨率SAR图像上,可见建筑物的一些细节.

    中等破坏建筑物的倒塌部分由于失去了原有的几何特性而形成零散的斑点或阴影,建筑物未倒塌部分在SAR图像上依然呈现规则的几何形状,倒塌部分表面结构杂乱 (刘金玉等,2013).在这种情况下,建筑物形成的二面角反射受到破坏,亮条纹在几何形状上变得不规则,在一片区域内均匀散射,其成像原理如图 4所示.

    图  4  建筑物震害散射特征示意图
    Figure  4.  Schematic diagram of scatter characteristics of earthquake-damaged buildings

    建筑物屋顶破坏但结构完整的中等破坏建筑物,如图 5c所示,其墙面与地面之间的二面角反射未被破坏,依然形成一条亮色条纹,但屋顶因破坏产生的参差结构形成了多次散射,亮度与破坏前相比增强,叠掩区与周围环境的反差降低,如图 6所示.

    图  5  汶川地震都江堰地区中等破坏建筑物光学图像与SAR图像
    (a) 震前谷歌图像;(b) 震后航片;(c) 震后COSMO-SkyMed图像
    Figure  5.  Optical image and SAR image of moderately damaged buildings in Dujiangyan area after the Wenchuan earthquake
    (a) Google image before the earthquake; (b) Aerial image after the earthquake; (c) COSMO-SkyMed image after the earthquake
    图  6  屋顶破坏建筑物光学图像 (a) 与SAR图像 (b).
    红色框表示建筑物破坏区域
    Figure  6.  Optical image (a) and SAR image (b) of roof-damaged building
    The damaged buildings are delineated by red squares

    建筑物结构和屋顶均未被破坏,只是一侧墙体倒塌,破坏一侧原有的二面角反射器消失,其形成的亮条纹在图像中消失不见,叠掩区的反射特征也发生变化,亮度与周围环境反差降低;而另一侧墙体完好,仍能形成阴影.

    建筑物受到地震毁坏时,二面角反射器受到破坏,建筑物破坏形成的瓦砾、水泥结构等碎块形成若干个小的反射个体,SAR图像中建筑物原有的几何结构模糊甚至消失,后向散射方式变为任意方向的多次散射 (赵凌君,2009),如图 7所示.毁坏建筑物在图像上的主要识别规律为图像亮度大,在城市建筑物区某位置可形成一团无规律的亮斑 (图 8).

    图  7  毁坏建筑物震害散射特征示意图
    Figure  7.  Schematic diagram of scatter characteristics of destroyed buildings
    图  8  毁坏建筑物 (红圈内) SAR图像
    Figure  8.  SAR image of destroyed buildings delineated by red circles

    与传统遥感数据相比,点云数据是对建筑物完整三维结构的记录,可以直观地表达目标的三维特征;通过点云数据可以直接量测建筑物的间距、三维坐标、方位角等信息,从而计算建筑物长度、面积和体积等信息,为建筑物地震破坏模式的分析和数值模拟提供重要数据基础,有利于推动地震工程领域的研究 (焦其松,2015).

    建筑物的破坏程度主要通过高度及形状结构的变化来反映,破坏特征在对应的点云数据中表现为点云空间几何形态的变化,包括高度改变、体积减小、建筑边界轮廓变化等,反映在建筑物震害上主要有以下4种类型.

    1) 建筑结构下部和端部的破坏较为严重.地震测区内砖混结构房屋下部的破坏普遍较上部严重得多,裂缝分布多,损毁较为严重.另外,很多房屋的端角部和整体建筑的端部均发生相对严重的破坏,甚至倒塌,如图 9所示.可见,破坏处已失去原来的平直表面,角边呈不规则状,表面粗糙,棱角状,原来的窗户边沿也显示为多边形状.

    图  9  结构破坏建筑物的LiDAR三维几何特征示意图
    (a) 顶面倾斜的建筑物; (b) 底层平行倒塌的建筑物
    Figure  9.  Schematic diagram of LiDAR three-dimensional geometry for buildings with structural destruction
    (a) A building with inclined top; (b) A building with parallel collapsed bottom

    2) 建筑物或局部倒塌.很多年代久远的房屋并未经过抗震设计,房屋的整体性差; 而砖混结构、砌体结构房屋的薄弱层由于应力集中,往往成为最先破坏的地方,最终导致整个房屋倒塌或局部倒塌,其LiDAR三维几何图如图 10所示.可以看到,倒塌建筑物整体呈不规则外形,缺乏整齐的平坦外观,仅外形显示建筑物的存在.

    图  10  倒塌建筑物的LiDAR三维几何示意图
    (a) 完全垮塌建筑物 (杂乱无章); (b) 完全垮塌建筑物 (似有外形)
    Figure  10.  Schematic diagram of LiDAR three-dimensional geometry for collapsed buildings
    (a) Completely collapsed building (disorderly and unsystematic); (b) Completely collapsed building (seems to have a shape)

    3) 墙体的震害破坏.大部分墙体在较强地震作用下均会产生裂缝,主要表现为:墙体呈严重的斜裂缝或“X”形裂缝,大多开裂处的砖块被剪断,窗下墙出现“X”形裂缝,因底层塌陷而导致的房屋倾斜或部分倒塌等,如图 11所示.可见,墙体倒塌,底部有不规则的堆积物,内墙或内框架暴露.

    图  11  墙体破坏建筑物的的LiDAR三维几何示意图
    (a) 建筑物墙体开裂;(b) 建筑物中部截断
    Figure  11.  Schematic diagram of LiDAR three-dimensional geometry for wall-damaged buildings
    (a) Building with wall crack; (b) Middle-cut building

    4) 大跨度厂房破坏.钢屋架、大跨度的厂房由于跨度大且没有足够的支撑,梁与承重墙之间的连接性差,在地震力作用下易发生坍塌破坏,如图 12所示.可以看到,建筑物大跨度破坏范围大,常见高的梁柱独立支撑,呈规则状垮落.

    图  12  大跨度破坏建筑物的LiDAR三维几何示意图
    Figure  12.  Schematic diagram of LiDAR three-dimensional geometry for a large-span damaged building

    将所研究的建筑物用理想化的几何模型来表达,即将建筑物的轮廓简化为长方体,按照屋顶形状将建筑物几何模型分为平顶建筑物与尖顶建筑物两种模型.通过高分辨率光学影像解译与实地调查图像的对比分析,确定建筑物图像的位置及几何形状,再据此提取SAR图像中的建筑物单体信息,对SAR图像沿建筑物墙面以一定间隔作剖面分析,并统计分析其后向散射特征,从而获取不同类型震害建筑物图像特征,为震害信息提取提供理论依据.

    假定建筑物是孤立的,墙面和周围的地面均平滑, 其SAR图像特征如图 13所示.

    图  13  平顶建筑物SAR图像特征示意图 (金鼎坚,2012)
    (a) θ<arctan (h/w); (b) θ>arctan (h/w)θ表示雷达入射角,w表示建筑物侧面宽度,h表示建筑物高度,g表示地面后向散射,c表示角反射,f表示墙面后向散射,r表示屋顶后向散射,s表示阴影.第二行表示地距图像上雷达回波的组成,条带高度表示回波的相对强度
    Figure  13.  Schematic diagram of SAR image of flat-top buildings
    (a) θ < arctan (h/w); (b) θ > arctan (h/w) θ represents incident angle of radar, w represents building side width, h represents building height, g represents ground backscattering, c represents angle reflection, f represents wall backscattering, r represents roof backscattering, s represents shadow. The second line represents the composition of radar echoes on the ground image, and the height of the strip indicates the relative intensity of the echo

    1) 当θ<arctan (h/w) 时 (图 13a),沿距离向由近及远,图像上分别为地面散射、叠掩、角反射、阴影等,其中叠掩包括地面回波、墙面回波和屋顶散射叠加在一起形成的叠掩 (g+f+r),地面和墙面的回波组成的叠掩如图中g+f所示.

    图 14为都江堰地区一个多层平顶建筑物的图像及实地调查照片,其宽度w小于高度h,满足条件θ<arctan (h/w),可简化为图 13a的几何模型.

    参照光学图像 (图 14a) 及实地调查照片 (图 14c) 可知,该建筑物属于完好建筑物.对SAR图像沿建筑物墙面以一定间隔作剖面分析,如图 14b中红色条带所示,根据条带方向求得平均剖面值,所得的沿距离向的后向散射强度如图 15所示.剖面线中沿距离向由近及远分别为地面的地面散射 (A)、叠掩和角反射 (B)、阴影 (C) 和地面散射 (D) 部分.

    图  14  多层平顶建筑物的光学图像
    (a)、SAR图像 (b) 与实地调查照片 (c) 图 (b) 中红色箭头表示视线方向,红色条带表示剖面线位置
    Figure  14.  Optical image (a), SAR image (b) and field survey image (c) of a multilayer flat-top building
    In Fig.(b), red arrow indicates direction of sight, and the red stripe represents the location of the profile line
    图  15  θ<arctan (h/w) 时的建筑物剖面线
    A为屋前地面散射部分;B为叠掩 (由地面回波、墙面回波和屋顶回波叠加形成) 部分,且存在地面和墙面形成的角反射效应;C为阴影部分;D为屋后地面散射部分
    Figure  15.  Building profile line on the condition of θ < arctan (h/w)
    A represents ground scattering in front of house; B represent layerover (formed by ground echo, wall echo and roof echo), and there is an angle reflection effect formed by ground and wall; C represents shadow; D represents ground scattering behind the house

    2) 当θ>arctan (h/w) 时 (图 13b),沿着距离向由近及远,图像上分别为地面的后向散射、叠掩、角反射、屋顶后向散射、阴影和地面散射等部分.此时地面回波、墙面回波和屋顶回波叠加在一起形成叠掩 (g+f+r);屋顶的后向散射除了一部分叠掩外,还另有一部分能被传感器单独接收,此时建筑物的阴影最少.

    图 16为都江堰一个多层平顶建筑物图像及实地调查照片,其高度h小于宽度w,满足条件θ>arctan (h/w),可简化为图 13b的几何模型.

    图  16  平顶建筑物的光学图像 (a)、SAR图像 (b) 与实地调查照片 (c)
    图 (b) 中红色箭头表示视线方向,红色条带表示剖面线位置
    Figure  16.  Optical image (a), SAR image (b) and field survey image (c) of the flat-top building
    In Fig.(b), red arrow indicates direction of sight, and the red stripe represents the location of the profile line

    参照光学图像 (图 16a) 和实地调查照片 (图 16c),该建筑物属于完好建筑物.对SAR图像沿建筑物墙面以一定间隔作剖面分析,如图 16b中红色条带所示,再根据条带方向求得平均剖面值,所得的沿距离向的后向散射强度如图 17所示.剖面线中沿距离向由近及远分别为地面散射、叠掩、角反射、屋顶散射、阴影和地面散射部分.

    图  17  θ>arctan (h/w) 时的建筑物剖面线
    A为屋前地面散射部分;B为叠掩 (由地面回波、墙面回波和屋顶回波叠加形成) 部分,且存在地面和墙面形成的角反射效应;C为屋顶散射部分;D为阴影部分;E为屋后地面散射部分
    Figure  17.  Building profile on the condition of θ > arctan (h/w)
    A represents ground scattering in front of house; Brepresents layerover (formed by ground echo, wall echo and roof echo), and there is an angle reflection effect formed by ground and wall; C represents roof scattering; D represents shadow; E represents ground scattering behind the house

    尖顶建筑物的散射特征与平顶建筑不同,除了叠掩、角反射和阴影特征之外,还有屋顶的单次反射效应.假设尖顶建筑物的理论模型为等腰三角形,倾斜屋面与平面形成的角度相同,尖顶建筑物示意图如图 18所示.以雷达波入射角θ与屋顶倾斜角度α的大小关系作为分类标准,分θα (图 18a) 和θα (图 18b) 两种情况予以分析.

    图  18  尖顶建筑物SAR图像特征示意图. (a) θα;(b) θα
    Figure  18.  Schematic diagram of SAR image of a spire building. (a) θ < α; (b) θ > α

    1) 当θα时 (图 18a),雷达波无法照射到背面屋顶,只有正面屋顶的回波与地面回波和墙面回波组成叠掩.

    2) 当θα时 (图 18b),屋顶的散射包括两部分,即来自正面屋顶的散射和来自背面屋顶的散射,因此在角反射效应之前的叠掩包括3部分:地面、正面屋顶和背面屋顶的回波组成的叠掩 (g+rf+rb),地面、墙面和背面屋顶的散射组成的叠掩 (g+f+rb),地面和墙面的散射组成的叠掩 (g+f).

    图 19为都江堰一个尖顶建筑物的图像和实地调查照片,其高度h小于宽度w,满足条件θα,可简化为图 18b的几何模型.

    图  19  尖顶建筑物的光学图像 (a)、SAR图像 (b) 和实地调查照片 (c)
    图 (b) 中,红色箭头表示视线方向,红色条带表示剖面线位置
    Figure  19.  Optical image (a), SAR image (b) and field survey image (c) of the spire building
    In Fig.(b), red arrow indicates direction of sight, and the red stripe represents the location of the profile line

    参照光学图像 (图 19a) 和实地调查照片 (图 19c) 可知,该建筑物属于完好建筑物.对SAR图像沿建筑物墙面以一定间隔作剖面分析,如图 19b中红色条带所示,再根据条带方向求得平均剖面值,所得的沿距离向的后向散射强度图如图 20所示.剖面线中沿距离向从左到右分别为地面散射、叠掩、角反射、屋顶散射、阴影和屋后地面散射部分.

    图  20  θα时建筑物剖面图
    A为地面散射部分; B为叠掩 (包括地面回波、墙面回波和正面屋顶回波叠加组成) 部分,且存在地面和墙面形成的角反射效应;C为屋顶散射部分; D为阴影部分; E为屋后地面散射部分
    Figure  20.  Building profile on the condition of θ > α
    A represents ground scattering; B represents layerover (formed by ground echo, wall echo and roof echo), and there is an angle reflection effect formed by ground and wall, C represents roof scattering, D represents shadow, E represents ground scattering behind the house

    完全倒塌的建筑物成为瓦砾堆,失去几何走向特征,规则的叠掩、角反射形成的亮线和阴影等特征消失,其后向散射强度一般略高于阴影.由于建筑材料中的瓦砾、钢筋在局部构成小的角反射器,形成较强回波,致使其表现为许多高亮的点. 图 21图 22中 (a) 和 (b) 分别为完全倒塌建筑物的光学图像和SAR图像,对其剖面分析分别如图 21c图 22c所示,图中出现无明显分布规律的峰值,这是由建筑物毁坏过程中小角反射器形成的强回波所造成的,对比两者可知,其剖面并无明显规律.

    图  21  损毁建筑物的光学图像 (a)、SAR图像 (b) 和SAR图像剖面线分析图 (c)
    图 (c) 中AA′为垂直于损毁建筑物屋顶的剖面
    Figure  21.  Optical image (a), SAR image (b) and SAR image profile (c) of the destroyed building
    In Fig.(c), AA′ is a profile perpendicular to the roof of the destroyed building
    图  22  损毁建筑物光学图像 (a)、SAR图像 (b) 和SAR图像剖面线分析图 (c)
    图 (c) 中AA′为垂直于损毁建筑物屋顶的剖面
    Figure  22.  Optical image (a), SAR image (b) and SAR image profile (c) of the destroyed building
    In Fig.(c), AA′ is a profile perpendicular to the roof of the destroyed building

    基于遥感图像的震害信息快速判识对于地震应急与救援决策具有重要的意义.本文通过分析建筑物震害在光学影像中的表现特征、在SAR图像中的散射机理以及在LiDAR图像中的三维特征,针对相同建筑物不同源遥感数据的联合分析,总结了不同震害程度建筑物在3类遥感影像中的表现特征.

    地震发生后,能够获取到大批量多传感器、多源遥感影像数据,不同类型的遥感影像数据对震害建筑物的表现形式不同,反映建筑物震害破坏形式的侧重点亦有所差异.震后光学遥感图像获取的色彩信息符合人眼色觉原理,便于人工判读,但是光学影像获取时易受到云雾天气的影响,而SAR图像能克服恶劣天气的影响,穿透云雾,且随着SAR图像分辨率的提高,其对建筑物表面细节信息的反映更精细. LiDAR图像能获取建筑物的三维信息,其中:机载LiDAR设备由于飞行高度较高,视场角不受地面障碍物的干扰,但其成本较高且不易获取建筑物的侧面信息;地面LiDAR能获取建筑物的侧面信息,但其视场角受障碍物遮挡影响较大且时效性不强.因此需要有效地利用多源遥感影像数据,如机载LiDAR图像与地面LiDAR图像的融合、LiDAR图像与高分辨率光学图像的融合以及光学图像与SAR图像的融合,通过不同类型影像特征的融合,深入开展多源遥感图像震害特征研究,充分发挥每种传感器类型影像的特长,相互补充,进而提高信息识别的能力与精度.

    多源遥感震害机理与特征的分析突破了传统震害分析仅限于图像灰度特征的不足,将遥感成像机理引入到震害分析中,从而更深入地挖掘遥感震害特征,辅助遥感震害信息的判识.未来可根据震害类型的不同,选择合适的提取方法与模型,同时通过充分挖掘3类遥感图像中建筑物震害的特征表现,进一步基于震害成像机理构建多源数据融合识别模型,从而为震害分析提供有效的理论支撑,也有助于多源数据融合的信息识别方法的研究,使遥感技术更好地为震害应急与评估服务.

  • 图  1   研究区域及台站概况

    (a) 区域构造背景、台站分布、历史强震及2008—2016年ML≥2.0地震分布;(b) 韩城台、龙门水文站、韩城断裂NE段和黄河的地理位置;(c) 韩城台及韩城断裂剖面图;(d) NS和EW分量金属水平摆;(e) 韩城台及周边地区与黄河之间地下水动力示意图

    Figure  1.   Map view of the investigation region and stations containing regional tectonic setting,earthquake events,hydrology and station locations

    (a) Regional tectonic setting,the distribution of ML≥2.0 events during the period from 2008 to 2016 and significant historical earthquakes;(b) Distribution of the Longmen hydrological station,Hancheng station,NE segment of Hancheng fault and Yellow River;(c) Picture showing the Hancheng station and Hancheng fault plane;(d) Picture of NS and EW components of the metallic horizontal pendulums;(e) Conceptual cartoon illustrating the groundwater hydrodynamics between Hancheng station and Yellow River

    图  2   韩城台地倾斜观测曲线

    (a) NS分量变化;(b) EW分量巨幅异常变化;(c) 地倾斜的矢量变化

    Figure  2.   The secular ground tilt recorded by Hancheng station

    (a) Tilt variations in NS component;(b) Remarkable tilt anomalies in EW component;(c) The hodograph showing the changes of tilt vector

    图  3   经验模态分解方法提取的地倾斜周年变化信号(a)和观测室气温变化(b)

    Figure  3.   Annual variations of tilt retrieved by EMD (a) and room temperature changes recorded at Hancheng observatory (b)

    图  4   1977年黄河特大洪峰引起的巨幅地倾斜

    (a) 龙门水文站记录的黄河特大洪峰;(b) 韩城台记录到的SE向巨幅地倾斜

    Figure  4.   Remarkable ground tilt induced by the great flood peak of the Yellow River in 1977

    (a) The 1977 great flood peak of the Yellow River recorded at Longmen hydrological station;(b) SE-ward tilt recorded at Hancheng station

    图  5   黄河水位(a)、S34井水位(b)及运城气象站(c,d)观测的降雨变化

    Figure  5.   Hydrology and rainfall information of Hancheng station and its adjacent regions

    (a) The water level of Yellow River recorded at Longmen hydrological station since 2007;(b) Changes in groundwater level in S34 well;(c) Time series of cumulative annual and daily precipitation observed at Yuncheng meteorological station;(d) Cumulative daily and detrended time series from 2001 to 2017

    图  6   断裂慢滑移引起的EW分量的理论倾斜场

    Figure  6.   Theoretical ground tilt in EW component induced by slow slip event on fault

    图  7   巨幅地倾斜异常前后研究区内85次2.0≤ML≤4.8地震的震源机制解及应力场反演结果

    (a) 2008年10月10日至2010年2月26日研究区内28次地震的震源机制解;(b) 2010年3月14日至2015年9月15日研究区内57次地震的震源机制解;(c) 2008年10月10日至2010年2月26日研究区内应力张量;(d) 2010年3月14日至2015年9月15日研究区内应力张量

    Figure  7.   Focal mechanism solutions for the 85 earthquakes with 2.0≤ML≤4.8 and stress field inversion in the studied area before and after the remarkable tilt anomalies

    (a) Focal mechanism solutions for the 28 earthquakes from Octorber 10,2018 to February 26,2010;(b) Focal mechanism solutions for the 57 earthquakes from March 14,2010 to September 15,2015;(c) The principal stress axes for the period from October 10,2008 to February 26,2010;(d) The principal stress axes for the period from March 14,2010 to September 15,2015

    表  1   韩城台概况

    Table  1   General information of Hancheng station

    台基
    岩性
    海拔
    /m
    观测室情况 金属水平摆仪器参数 仪器支墩
    材质及尺寸
    覆盖层
    厚度/m
    室内气温
    年变幅/℃
    室内气温
    日变幅/℃
    相对
    湿度
    仪器
    类型
    记录
    方式
    折合摆长
    /mm
    使用周期
    /s
    奥陶系
    灰岩
    460 0 17—18 ≤0.5 90% JB 光记录 NS分量:25.2
    EW分量:25.4
    18—19 混凝土,
    高0.7 m
    下载: 导出CSV
  • 薄万举. 2010. 形变异常与干扰关系的再认识[J]. 大地测量与地球动力学,30(1):5–8.

    Bo W J. 2010. Study on relation between crustal deformation anomaly and disturbances[J]. Journal of Geodesy and Geodyna-mics,30(1):5–8 (in Chinese).

    高伯贤,高雪. 2011. 韩城矿区南部奥陶系灰岩岩溶水特征[J]. 陕西煤炭,30(1):43–45.

    Gao B X,Gao X. 2011. Characteristics of Ordovician limestone karst water in southern Hancheng mining area[J]. Shaanxi Coal,30(1):43–45 (in Chinese).

    郭平战. 2015. 韩城矿区岩溶形成机理及其水文地质特征[J]. 地下水,37(5):54–57.

    Guo P Z. 2015. Karst formation mechanism of karst in Hancheng mining area and its hydrogeological characteristics[J]. Ground Water,37(5):54–57 (in Chinese).

    何毅. 2012. 近60年来渭河流域气候变化研究[D]. 杨凌: 西北农林科技大学: 16–37.

    He Y. 2012. Climate Change of Wei River Basin in Last 60 Years[D]. Yangling: Northwest A&F University: 16–37 (in Chinese).

    黄辅琼,陈颙,白长清,张晶,晏锐,杨明波,兰从欣,张晓东,江在森. 2005. 八宝山断层的变形行为与降雨及地下水的关系[J]. 地震学报,27(6):637–646.

    Huang F Q,Chen Y,Bai C Q,Zhang J,Yan R,Yang M B,Lan C X,Zhang X D,Jiang Z S. 2005. The correlation of deformation behavior with precipitation and groundwater of the Babaoshan fault in Beijing[J]. Acta Seisomogical Sinica,27(6):637–646 (in Chinese).

    扈桂让,李自红,闫小兵,赵晋泉,曾金艳,郭瑾. 2017. 韩城断裂晚第四纪活动性研究[J]. 地震地质,39(1):206–217.

    Hu G R,Li Z H,Yan X B,Zhao J Q,Zeng J Y,Guo J. 2017. The study of Late Quaternary activity of Hancheng fault[J]. Seismology and Geology,39(1):206–217 (in Chinese).

    陆一锋,徐鸣洁,王良书,米宁,李华,于大勇. 2012. 鄂尔多斯东南缘地区的地壳结构[J]. 科学通报,57(1):59–64.

    Lu Y F,Xu M J,Wang L S,Mi N,Li H,Yu D Y. 2011. Crustal structure of the southeastern margin of the Ordos block[J]. Chinese Science Bulletin,56(35):3854–3859. doi: 10.1007/s11434-011-4847-7

    牛安福,顾国华,曹景鹏,张凌空,闫伟,赵静,吉平. 2013. 芦山MS7.0地震前远、近场形变时空演化特征研究[J]. 地震学报,35(5):670–680.

    Niu A F,Gu G H,Cao J P,Zhang L K,Yan W,Zhao J,Ji P. 2013. On the preseismic deformation changes prior to the Lushan MS7.0 earthquake[J]. Acta Seismologica Sinica,35(5):670–680 (in Chinese).

    水利部黄河水利委员会. 2017. 水情信息[EB/OL]. [2018–01–21]. http://61.163.88.227:8006/hwsq.aspx.

    Yellow River Conservancy Commission of the Ministry of Water Resources. 2017. The information of water[EB/OL]. [2018–01–21]. http://61.163.88.227:8006/hwsq.aspx (in Chinese).

    吴富春,张鸿福,景北科,段锋,张义民. 1999. 陕西几例特殊的无震异常及其成因分析[J]. 西北地震学报,21(3):268–273.

    Wu F C,Zhang H F,Jing B K,Duan F,Zhang Y M. 1999. Study on some special aseismic precursory anomalies observed from Shaanxi seismic network and their causes[J]. Northwestern Seismological Journal,21(3):268–273 (in Chinese).

    田中豊. 1972. 地殼変動連続観測(長期変動と地震前駆変動)[C]//地震予知研究シンポジウム. 东京: 东京大学出版社: 35–44.

    Tanaka Y. 1972. Long-term crustal deformation and precursory variation in seismicity[C]//Proceedings of the Earthquake Prediction Research Symposium 1972. Tokyo: University of Tokyo Press: 35–44 (in Japanese).

    里嘉千茂. 1989. 伝播性歪に関する数値実験[J]. 測地学会誌,35(1):27–36.

    Sato K. 1989. Numerical experiments on strain migration[J]. Journal of the Geodetic Society of Japan,35(1):27–36 (in Japanese).

    内藤宏人,吉川澄夫. 1999. 地殻変動解析支援プログラムMICAP-Gの開発[J]. 地震,52(1):101–103.

    Naito H,Yoshikawa S. 1999. A program to assist crustal deformation analysis[J]. Journal of the Geodetic Society of Japan,52(1):101–103 (in Japanese).

    竹本修三, 和田安男, 伊藤潔, 福田洋一, 森井亙, 百瀬秀夫, 中村光邦. 2004. 地殻ひずみの観測に及ぼす局所的日照変化の影響—花山と立山観測室のデータ比較[R].京都大学防災研究所年報, 47(B): 725–734.

    Takemoto S, Wada Y, Ito K, Fukuda Y, Morii W, Momose H, Nakamura M. 2004. Effect of local sunshine changes upon crustal strain observations: Comparison of strain data obtained at Kwasan and Tateyama stations[R]. Disaster Prevention Research Institute Annuals, Kyoto University, 47(B): 725–734 (in Japanese).

    Agnew D C. 1986. Strainmeters and tiltmeters[J]. Rev Geophys,24(3):579–624. doi: 10.1029/RG024i003p00579

    Bilham R G,Beavan R J. 1979. Strains and tilts on crustal blocks[J]. Tectonophysics,52(1/2/3/4):121–138. doi: 10.1016/0040-1951(79)90216-6

    Braitenberg C,Nagy I. 2014. Illustrating the superposition of signals recorded by the Grotta Gigante pendulums with musical analogues[J]. Acta Carsol,43(1):139–147.

    Brimich L,Bednárik M,Bezák V,Kohút I,Bán D,Eper-Pápai I,Mentes G. 2016. Extensometric observation of Earth tides and local tectonic processes at the Vyhne station,Slovakia[J]. Contrib Geophys Geod,46(2):75–90. doi: 10.1515/congeo-2016-0006

    Bykov V G,Trofimenko S V. 2016. Slow strain waves in blocky geological media from GPS and seismological observations on the Amurian plate[J]. Nonlin Process Geophys,23(6):467–475. doi: 10.5194/npg-23-467-2016

    Caniven Y,Dominguez S,Soliva R,Peyret M,Cattin R,Maerten F. 2017. Relationships between along-fault heterogeneous normal stress and fault slip patterns during the seismic cycle:Insights from a strike-slip fault laboratory model[J]. Earth Planet Sci Lett,480:147–157. doi: 10.1016/j.jpgl.2017.10.009

    Cicerone R D,Ebel J E,Britton J. 2009. A systematic compilation of earthquake precursors[J]. Tectonophysics,476(3/4):371–396.

    Detournay E, Cheng A H D. 1993. Fundamental of Poroelasticity in Comprehensive Rock Engineering: Principles, Practice & Projects, Vol.2[M]. Oxford: Pergamon Press: 127–128.

    Dragoni M,Bonafede M,Boschi E. 1984. On the interpretation of slow ground deformation precursory to the 1976 Friuli earthquake[J]. Pure Appl Geophys,122(6):781–792.

    Eper-Pápai I,Mentes G,Kis M,Koppán A. 2014. Comparison of two extensometric stations in Hungary[J]. J Geodyn,80:3–11. doi: 10.1016/j.jog.2014.02.007

    Evans K,Wyatt F. 1984. Water table effects on the measurement of earth strain[J]. Tectonophysics,108(3/4):323–337.

    Fréchet J,Rivera L. 2012. Horizontal pendulum development and the legacy of Ernst von Rebeur-Paschwitz[J]. J Seismol,16(2):315–343. doi: 10.1007/s10950-011-9272-5

    Fukuyama E. 2015. Dynamic faulting on a conjugate fault system detected by near-fault tilt measurements[J]. Earth Planets Space,67:38. doi: 10.1186/s40623-015-0207-1

    Gershenzon N I,Bykov V G,Bambakidis G. 2009. Strain waves,earthquakes,slow earthquakes,and afterslip in the framework of the Frenkel-Kontorova model[J]. Phys Rev E,79(5):056601. doi: 10.1103/PhysRevE.79.056601

    Goulty N R. 1976. Strainmeters and tiltmeters in geophysics[J]. Tectonophysics,34(3):245–256.

    Hao M,Wang Q L,Cui D X,Liu L W,Zhou L. 2016. Present-day crustal vertical motion around the Ordos block constrained by precise leveling and GPS data[J]. Surv Geophys,37(5):923–936. doi: 10.1007/s10712-016-9375-1

    Harada M,Furuzawa T,Teraishi M,Ohya F. 2003. Temporal and spatial correlations of the strain field in tectonic active region,southern Kyusyu,Japan[J]. J Geodyn,35:471–481. doi: 10.1016/S0264-3707(03)00008-5

    Hardebeck J L,Shearer P M. 2002. A new method for determining first-motion focal mechanisms[J]. Bull Seismol Soc Am,92(6):2264–2276. doi: 10.1785/0120010200

    Harrison J C,Herbst K. 1977. Thermoelastic strains and tilts revised[J]. Geophys Res Lett,4(11):535–537. doi: 10.1029/GL004i011p00535

    Hisz D B,Murdoch L C,Germanovich L N. 2013. A portable borehole extensometer and tiltmeter for characterizing aquifers[J]. Water Resour Res,49(12):7900–7910. doi: 10.1002/wrcr.20500

    Huang N E,Shen Z,Long S R,Wu M C,Shih H H,Zheng Q,Yen N C,Tung C C,Liu H H. 1998. The empirical mode decomposition and the Hilbert spectrum for nonlinear and non-stationary time series analysis[J]. Proc Roy Soc A Math Phys Eng Sci,454(1971):903–995. doi: 10.1098/rspa.1998.0193

    Ishii H,Sato T,Takagi A. 1980. Characteristics of strain migration in the northeastern Japan Arc (II):Amplitude characteris-tics[J]. J Geod Soc Japan,26(1):17–25.

    Jahr T,Jentzsch G,Gebauer A,Lau T. 2008. Deformation,seismicity,and fluids:Results of the 2004/2005 water injection experiment at the KTB/Germany[J]. J Geophys Res,113(B11):B11410. doi: 10.1029/2008JB005610

    Kartvelishvili K Z. 2010. Investigation of deformational processes in Tbilisi underground earth-tidal laboratory[J]. J Georgian Geophys Soc,Phys Atmos,Ocean and Space Plasma,14(B):197–201.

    Kawai K,Sekine S,Fuji N R,Geller R J. 2009. Waveform inversion for D″ structure beneath northern Asia using Hi-net tiltmeter data[J]. Geophys Res Lett,36(20):L20314. doi: 10.1029/2009GL039651

    Kimura T,Obara K,Kimura H,Hirose H. 2011. Automated detection of slow slip events within the Nankai subduction zone[J]. Geophys Res Lett,38(1):L01311. doi: 10.1029/2010GL045899

    Linde A T,Gladwin M T,Johnston M J S,Gwyther R L,Bilham R G. 1996. A slow earthquake sequence on the San Andreas fault[J]. Nature,383(6595):65–68. doi: 10.1038/383065a0

    Lindsey E O,FialkoY,Bock Y,Sandwell D T,Bilham R. 2014. Localized and distributed creep along the southern San Andreas fault[J]. J Geophys Res,119(10):7909–7922. doi: 10.1002/2014JB011275

    Martínez-Garzón P,Kwiatek G,Ickrath M,Bohnhoff M. 2013. MSATSI:A MATLAB package for stress inversion combining solid classic methodology,a new simplified user-handling and a visualization tool[J]. Seismol Res Lett,85(4):896–904.

    Mentes G. 2008. Observation of recent tectonic movements by extensometers in the Pannonian Basin[J]. J Geodyn,45(4/5):169–177.

    Mentes G. 2017. The role of recent tectonics and hydrological processes in the evolution of recurring landslides on the Danube’s high bank in Dunaföldvár,Hungary[J]. J Geodyn,290:200–210.

    Michelson A A. 1914. Preliminary results of measurements of the rigidity of the earth[J]. Astrophys J,39:105–138. doi: 10.1086/142058

    NOAA. 2017. Climate data online[EB/OL]. [2017−10−12]. https://www.ncdc.noaa.gov/cdo-web/.

    Okada Y. 1985. Surface deformation due to shear and tensile faults in a half-space[J]. Bull Seismol Soc Am,74(5):1135–1154.

    Rikitake T. 1987. Earthquake precursors in Japan:Precursor time and detectability[J]. Tectonophysics,136(3/4):265–282.

    Rikitake T. 1988. Earthquake prediction:An empirical approach[J]. Tectonophysics,148(3/4):195–210.

    Schuite J,Longuevergne L,Bour O,Burbey T J,Boudin F,Lavenant N,Davy P. 2017. Understanding the hydromechanical behavior of a fault zone from transient surface tilt and fluid pressure observations at hourly time scales[J]. Water Resour Res,53(12):10558–10582. doi: 10.1002/2017WR020588

    Sgrigna V,D'ambrosio C,Yanovskaya T B. 2002. Numerical modeling of preseismic slow movements of crustal blocks caused by quasi-horizontal tectonic forces[J]. Phys Earth Planet Int,129(3/4):313–324.

    Stein R S,Barka A A,Dieterich J H. 1997. Progressive failure on the North Anatolian fault since 1939 by earthquake stress trig-gering[J]. Geophys J Int,128(3):594–604. doi: 10.1111/gji.1997.128.issue-3

    Takemoto S. 1991. Some problems on detection of earthquake precursors by means of continuous monitoring of crustal strains and tilts[J]. J Geophys Res,96(B6):10377–10390. doi: 10.1029/91JB00239

    Timofeev V Y,Ardyukov D G,Boyko E V,Gribanova E I,Semibalamut V M,Timofeev A V,Yaroshevich A V. 2012. Strain and displacement rates during a large earthquake in the South Baikal region[J]. Russ Geol Geophys,53(8):798–816. doi: 10.1016/j.rgg.2012.06.007

    Timofeev V Y,Masalsky O K,Ardyukov D G,Timofeev A V. 2015. Local deformation and rheological parameters by measurements in Talaya station gallery (Baikal region)[J]. Geodyn Tectonophys,6(2):241–253. doi: 10.5800/GT-2015-6-2-0179

    Tsai V C. 2011. A model for seasonal changes in GPS positions and seismic wave speeds due to thermoelastic and hydrologic variations[J]. J Geophys Res,116(B4):B04404. doi: 10.1029/2010JB008156

    Wang H F. 2000. Theory of Linear Poroelasticity[M]. Princeton: Princeton University Press: 265–266.

    Yamazaki K. 2013. An attempt to correct strain data measured with vault-housed extensometers under variations in temperature[J]. Tectonophysics,599:89–96. doi: 10.1016/j.tecto.2013.04.001

    Zadro M,Braitenberg C. 1999. Measurements and interpretations of tilt-strain gauges in seismically active areas[J]. Earth Sci Rev,47(3):151–187.

  • 期刊类型引用(3)

    1. 刘同振,高原. 青藏高原东北缘地壳地震各向异性研究进展. 中国地震. 2023(02): 225-242 . 百度学术
    2. 蔡光耀,王未来,吴建平,房立华. 鄂尔多斯及邻区基于程函方程的面波层析成像. 地球物理学报. 2021(04): 1215-1226 . 百度学术
    3. 莘海亮,曾宪伟,康敏,高级. 海原弧形构造区地壳三维精细速度结构成像. 地球物理学报. 2020(03): 897-914 . 百度学术

    其他类型引用(4)

图(7)  /  表(1)
计量
  • 文章访问数:  1939
  • HTML全文浏览量:  1115
  • PDF下载量:  110
  • 被引次数: 7
出版历程
  • 收稿日期:  2018-02-23
  • 修回日期:  2018-04-18
  • 网络出版日期:  2018-08-19
  • 发布日期:  2018-10-31

目录

/

返回文章
返回