地震破裂带附近同震位移和变形模式对于深入理解地震破裂过程、断层行为及活动断层与地形地貌关系等至关重要。文中提出一种迭代最近点(ICP)算法, 利用地震前后地形点云进行差分确定断层近场三维同震地表位移。在川西大凉山断裂带交际河断层上选取2期SfM地形点云叠加同震位移场模拟震前-震后点云, 测试ICP方法获取同震位移场的精度。该方法在网格边长>50m的条件下可准确恢复同震位移场的方向和幅度, 水平和垂直精度分别为20~10cm, 这与地形点云定位精度相当。随着点云密度和网格窗口尺寸减小, 该方法恢复同震位移场等的精度将降低。通过分析树木生长、房屋建设、河流侵蚀等地形变化对恢复位移场的潜在影响, 研究结果表明扩大网格尺寸可使震前-震后点云具有足够地形结构进行匹配, 以减小局部地形变化对恢复位移场的影响, 网格窗口尺寸是在具有足够地形结构的大尺度和具有更精细分辨率的小尺度之间的权衡。文中所述的ICP方法利用震前-震后高精度点云可获取地震破裂带附近精细的三维地表位移场, 为浅层断层滑动和破裂带变形提供了新的约束, 有助于研究地震破裂过程和断层生长过程。
2025年1月7日9时5分, 西藏日喀则市定日县发生6.8级地震, 震源深度10km, 地震最高烈度达Ⅸ度。文中通过高分2号遥感数据开展震前与震后影像解译, 揭示发震断裂登么错断裂的同震地表破裂在几何样式上具有南、北段分异性, 北段长约3km, 表现为几何结构相对简单、破裂宽度较窄的“集中型破裂”样式, 具有连续性和继承性; 南段长约12km, 表现为几何结构较为复杂、破裂宽度较宽的“弥散型破裂”样式, 具有断续性和新生性。震中区附近28个行政村的房屋倒塌数量和面积统计显示, 震中附近3个行政乡的受灾程度由高到低依次为长所乡、措果乡、曲洛乡。将受灾行政村按地质地理环境分区, 受灾面积由大到小依次为近微观震中区、靠近地表破裂带的湖泊区、远离微观震中和地表破裂的近基岩山区。同震地表破裂解译表明, 登么错湖以北和以南存在2个未破裂的断层区段, 结合地震由南向北单侧破裂的过程及震后余震分布的特点, 认为南段未破裂区的地震危险性可能更高。在区域尺度上, 曲洛断裂、登么错断裂、郭加断裂、定结断裂这4条断裂组成的正断层对本次地震均有震后响应, 考虑近期中强地震能量释放的特点, 曲洛断裂和定结断裂的地震危险性值得关注。最后, 依据此次地震特点, 文中提出在高寒高海拔、地质地貌类型多样、断层变形样式弥散、少数民族聚居地区开展震害防御的对策与建议。
随着断层面形貌测量、 宇宙成因核素测年等技术的发展, 基岩断层面成为在基岩区古地震研究的良好研究对象。但是, 受沟谷侵蚀、 沉积覆盖和人为改造等非构造因素影响, 并非所有出露的基岩断层面都能完整记录和保存古地震信息。因此, 必须在基岩断层面上选择合适的研究点。传统的野外测量方法耗时费力, 难以掌握基岩断层面的整体信息, 即使通过高分辨率卫星影像也无法获得基岩断层面的精细结构。小型无人机航测(sUAV)采用低空摄影测量技术, 可快速获取基岩断层面高精度的地表三维结构, 为筛选目标工作点提供便利。文中以山西地堑系六棱山北麓断裂的马家窑基岩断坎为例, 采用小型无人机航测结合SfM摄影测量技术, 获得了基岩断层面的高精度三维地形数据, 通过精细地貌解译和断层坎剖面分析识别了沿断层坎分布的崩塌破坏、 沉积埋藏和侵蚀剥露等非构造因素, 分析了侵蚀、 埋藏和构造这3种典型的基岩断层面的形态特征, 选取了适合开展古地震研究的基岩断层面点, 展现了小型无人机在基岩区古地震研究工作中发挥的重要应用潜力。
同震地表破裂带的空间展布及形变特征是地球深部断层活动在地表的直观地貌表现, 不但记录着地震破裂和断层运动的信息, 还反映了区域应力和地壳运动状况。因此, 开展震后地震地表破裂带调查对于了解发震断层的构造活动尤为重要。高精度地形观测技术可以获取前所未有的高时空分辨率的地球表面特征, 为辨别历史地震地表破裂遗迹、 提取地表同震位移、 活动构造地质填图等提供高质量数据。文中选取富蕴1931年地震地表破裂带作为研究区, 利用SfM(Structural from Motion)摄影测量技术生成分辨率为1m的数字高程模型(DEM), 详细识别地表破裂并测量冲沟的右旋位移。基于地表破裂的几何及构造地貌特征, 将富蕴地震地表破裂带由北向南分为S1、 S2、 S3、 S4 4段, 其间以挤压隆起或拉分盆地相连接。沿破裂带共获得194组最新冲沟的右旋水平位移, 得到1931年同震位移的平均值为(5.06±0.13)m。同震位移局部缺失或突变的区域与几何阶区的位置也有良好的对应关系。以上结果填补了对富蕴地震地表破裂精细形态研究的空白, 也进一步展示高分辨率的地形数据在活动构造研究中良好的应用价值。
研究地震地表破裂的几何结构与同震位移分布的空间关联性, 对于理解断层面几何形貌与地震活动的相关性具有重要意义。 文中运用滑动窗口方法, 对28个震例的地表破裂非规则几何结构特征(走向变化、 阶区大小和粗糙度)和同震位移分别进行了量化和均值计算, 并计算了二者之间的相关性, 最后分析讨论了走滑型和倾滑型震例的相关性之间的差异性及成因。 研究结果表明: 1)相对于倾滑型, 走滑型地震地表破裂的几何结构与同震位移的相关性更高, 表明走滑型断层(破裂)的几何特征在分段中具有更高的参考价值; 2)在走滑型震例中, 地表破裂的几何结构与同震位移之间存在负相关关系, 即地表破裂上某段的阶区、 走向变化、 粗糙度的特征值越大, 其对应的同震位移量越小; 3)走滑型地震的地表破裂非规则几何结构特征中, 与同震位移分布的相关程度由高到低依次为阶区大小、 走向变化、 粗糙度。
探槽揭露的地质构造信息十分珍贵,是正确认识断层活动性和古地震信息的关键资料,但是探槽本身无法长久妥善保存,因此完整精确地记录探槽揭露的地质构造信息一直受到人们的重视。传统的探槽记录方法主要有素描和照片拼接2种,但前者缺乏色彩信息且包含过多主观解释内容;后者在照片几何拼接过程中会出现接边、畸变等问题。文中以甘肃省天桥沟-黄羊川断裂的柳条河探槽为例,介绍1种完整记录探槽高精度的三维地质构造信息的技术方法和工作流程。该方法基于三维图像建模和正射镶嵌技术,首先在野外对清理好的探槽建立控制点或比例尺,用数码相机采集探槽所有剖面及周边微地貌照片,然后在室内用图像建模软件进行探槽的三维图像建模,最后生成探槽的三维图像模型及其剖面的正射镶嵌图。结果表明:探槽的三维图像建模既可以实现高分辨率正射照片镶嵌的自动化,又可以高效、低成本地获得高精度的探槽三维图像模型,不仅有利于探槽信息的存储,更可以直观地获取更多地质构造三维信息。
构造地貌能够记录构造运动的长期积累,并能通过相关参数定量表达。但是,由于地貌还受到气候和岩性等因素的影响,所以如何利用相关参数来表达地貌与构造之间的关系是一个研究热点。文中利用水平分辨率为90m的SRTM数据,采用元分维模型,通过变差函数法计算了滇西南地区的二维分形参数。并对比了在不同岩性不同气候条件下分形参数的分布特征。结果显示,二维分形参数与岩性、气候因素相关性不强。而分维值与构造活动因素相关性很好,在构造活动强烈区域分维值较低,在构造活动不强烈区域分维值较高。这说明分维值能在一定程度上反映构造活动的强弱,为研究区域构造活动提供了一种新的手段。
地震断层滑移过程中的摩擦升温能够指示断层的动态摩擦性质,为我们认识地震的动力学过程及能量分配提供一种途径;断层带摩擦生热的残留标志也能作为识别地震活动的依据。镜质体反射率是煤和石油天然气领域广泛使用的地质温度计,在断层岩和断层力学研究中也具有潜在的应用价值。文中对采自汶川地震断裂带的多种断层岩样品开展了镜质体反射率研究,揭示了断层活动对断层岩镜质体反射率的影响;从断层角砾带向断层核部,镜质体反射率呈现出增大的趋势,指示着地震过程中断层的摩擦生热过程。八角庙露头样品的镜质体反射率明显高于深溪沟露头,可能暗示着2个露头上断层活动性质的差异。研究还表明,对断层剖面的镜质体反射率开展系统性研究有助于识别以往地震的滑动带以及限定滑动带的宽度。对深溪沟露头上断层滑动面样品开展的测试表明,2条宽度约2mm的局部黑色滑动带内存在镜质体反射率的异常高值;数值计算表明汶川地震过程中该露头附近断层滑动面上的宏观摩擦升温可能<162℃,同震阶段的动态摩擦系数较低,这些估计与高速摩擦实验的结果总体上是一致的。
大凉山断裂带是中国大陆大型强震断裂带鲜水河-小江断裂系的重要组成部分,其中段的普雄断裂是断裂带中最长的1条次级断裂。该次级断裂的古地震活动情况对评价该区域地震活动水平和建立地震灾害预防措施具有重要意义。近年开展的活动断层填图和古地震探槽研究表明普雄断裂是1条晚更新世以来活动强烈的略向W倾的高角度左旋走滑断层。沿断裂开挖的2个探槽分别揭露了2次和3次古地震事件,分别发生在8206 BC-1172 AD、1084-1549 AD和17434-7557 BC、1577-959 BC、927-1360 AD。结合探槽古地震结果和历史地震记载,进行建模分析得到断裂的离逝时间约为0.7ka,与距今倒数第2次事件的时间间隔约为2.3ka。同时,根据震级与地表破裂长度的经验公式估算古地震事件震级在7级以上。
探槽技术虽然已经发展成为古地震研究的重要手段, 结合定年技术可以识别出古地震事件与地震重复间隔, 但是仍然存在较大的不确定性和局限性。例如, 目视解译难以区分黄土沉积内部的细层理, 严重影响了古地震事件发生位置或时间的判定。如何提高古地震研究的精度和准确性, 降低古地震事件判定的不确定性, 是目前面临的一个迫切问题。山西洞峪沟剖面位于山西临汾盆地东北角, 横跨霍山山前断裂带, 不仅揭示了较好的黄土沉积序列, 还揭露了明显的地层错断事件。因此, 该剖面是一个开展高精度探槽古地震研究, 降低判定古地震事件不确定性的理想场所。根据高精度的粒度与磁化率变化曲线, 结合目视分层解译结果, 对洞峪沟黄土剖面进行了精细分层, 界定了各层的厚度和边界。依据细分层和断层两盘地层的对应关系, 将u6层沉积以来的3次断错事件的发生位置和时间作了再限定, 它们分别发生在u5-7、u4以及u2的顶部, 对应断层上盘埋深7.1m、4.7m与2.9m。根据释光测年结果以及断层上盘地层的平均沉积速率, 推测3次断错事件发生的时间分别在 (48.1±1.5)~(43.2±2.5)ka BP接近45.8ka BP、(35.0±2.4)~(30.6±1.3)ka BP接近32.8kaBP、(26.4±0.8)~(20.9±0.7)ka BP接近23.3ka BP。根据3次黄土-古土壤沉积旋回的厚度差, 判定3次地震的同震垂直位移分别为0.5m、0.4m和1.3m, 累计位移2.2m。依据高精度粒度与磁化率变化曲线的地层分层方法, 为有效降低黄土地区古地震研究的不确定性提供了一种很好的参考。
地形数据的质量(精度和分辨率)影响着地球科学的研究水平。LiDAR测量是目前获取高分辨率地形数据的有效技术方法之一, 但是其高昂的测量成本和相对复杂的后期数据处理限制了LiDAR技术的大众化应用。近年来, 一种被称为SfM(Structure from Motion)的适合大众化使用的新的高精度3维地形数据获取技术开始引起人们的注意。这种新型数字摄影测量技术可以利用高效的图像特征匹配算法从多视角照片中提取重叠区域的3维地形数据。由于SfM技术仅需要目标物体的照片, 而且对相机拍摄位置、图像尺度及拍摄焦距没有要求, 因此利用简单测量平台采集地面照片就可以获取高质量的3维地形数据。与LiDAR技术相比, 大大降低了获取高精度数据的成本, 使得高精度3维地形数据的使用大众化。文中介绍了SfM技术的基本原理和流程, 展示了SfM技术获取高精度3维地形数据的简单而有效的特性, 特别适合于植被稀少的区域。文中利用近千米高空拍摄的、具有约70%重叠度的一套随LiDAR飞行采集的数字航空照片生成具有真彩色的高密度SfM点云数据, 点密度高达25.5个/m2, 可生成分辨率0.2m的DEM(数字高程模型)。对比相同区域的LiDAR点云数据, 统计分析表明58.3%的LiDAR数据与SfM数据的垂直偏差<0.1m, 88.3%的LiDAR数据的垂直偏差<0.2m; 而且发现不同地貌的SfM数据精度存在差异, 平缓地形的SfM数据精度高于陡峭地形的SfM数据精度。文中还介绍了以氦气球作为拍摄平台的SfM测量系统, 可以快捷地获取高精度的3D地形数据和正射影像, 比目前常用的差分GPS测量具有更高的效率和数据精度。
大凉山断裂带是川滇菱形活动块体东边界的重要组成部分, 其活动习性和滑动速率对于了解青藏高原东南缘的块体运动和构造变形具有重要的理论意义, 同时对于该地区地震危险性评价和地震中长期预测有着重要的应用价值。大凉山断裂带由6条分支断裂构成, 竹马断裂是其北段的一条分支。文中在详细野外调查的基础上, 通过典型构造地貌的高精度GPS测量和断错时间的约束, 得出竹马断裂是一条以左旋走滑为主的全新世活动断裂, 全新世以来的水平滑动速率为1.5~3.1mm/a; 利用古地震探槽揭示出竹马断裂的2次古地震事件时间分别为(50.3±5.7)~30ka BP和30~(17.4±1.2)ka BP。考虑到断裂带北段还存在一条与竹马断裂平行的分支断裂——公益海断裂, 大凉山断裂带北段的滑动速率与南段基本一致。鲜水河-小江断裂系中段在大凉山断裂带上分配的滑动速率与并行的安宁河断裂带和则木河断裂带也大体相当, 断裂系中段的滑动速率之和与南北两段的速率大致吻合。大凉山断裂带的存在使断裂系在几何上成为一个完整的弧形构造, 也弥补了中段滑动速率的亏损, 使各段的滑动保持协调一致。此外, 沿竹马断裂大量发育的冲洪积台地, 根据沉积特征和测年结果, 认为是新仙女木事件后末次冰消期的冰川融水形成的冰水堆积, 类似的地貌面可能广泛分布于青藏高原东南缘。
基岩断层面形貌特征的定量分析方法是研究断层活动历史和识别古地震的有效手段, 是对探槽技术的有力补充, 可以弥补在基岩区开展古地震研究技术手段的严重不足。文中以山西霍山山前断裂上的断层面(崖)为例, 采用各向同性变差函数法计算了断层面形貌的2D分维值, 并通过沿断层走向逐条平均的方式获得了2D分维值在断层倾向上的分布特征, 结果显示断层面形貌在倾向上具有显著的分带性特征, 由此获得的每个分带的特征分维值随断层面高度的增加呈阶跃式增加。这种断层面形貌在倾向上的显著分带性特征、特征分维值的阶跃式增加反映了断层面的出露方式不是渐进式, 而是不连续的间歇式; 这种间歇式分带性出露方式与断层的周期性地震活动相关, 每个分带的特征值与每个分带的出露时间相关, 即与地震发生时间相关。因此: 1)可以通过每一个断层面分带的特征分维值估计该分带出露的时间; 2)断层面分带的宽度可以用来估计每次破裂地震的同震位移量。霍山山前断裂带上3个基岩断层面形貌分带揭示出3次地震事件; 3个2D特征分维值不仅显示3个分带自上而下由老变新的阶跃特征, 而且3个阶跃之间大致相等的特征分维值差显示了3次地震发生间隔是大致相等的; 断层面分带的宽度反映出破裂地震的倾向同震位移量大约为3.5m。此外, 在阶跃式分带之间还存在宽0.5~1m的较窄分带, 这些窄分带的特征分维值是随断层面高度增加而逐渐增大的, 显示在间震期基岩断层面被逐渐剥露的特征。
机载LiDAR技术为描绘活动构造相关构造地貌和最新的地表形变提供更精确的基础数据。如何将LiDAR新技术、新数据应用于活动构造填图和活动断层地震危险性评价等方面,是今后活动构造研究领域的一个重要的发展方向。文中以新疆天山北麓的独山子背斜-逆冲断裂带为试验区,开展了基于LiDAR数据的活动构造填图实验研究。首先,采用机载LiDAR技术进行数据采集,获得点云密度为6.6个/m2、平均点间距为0.39m的LiDAR原始数据;其次,利用试验区内12个测量精度可达mm级的GPS静态测量点评估LiDAR的相对垂直精度为0.12m、均方差值为0.078m;最后,对密度为6.4个/m2的地面点云数据进行DEM最佳分辨率评估,利用反距离权重算法获得0.5m分辨率的数据高程模型(DEM)。该分辨率的DEM数据足以完成独山子背斜-逆冲断裂带的精细构造地貌特征的确定以及高精度的空间解译。文中仅使用DEM可视化工具从不同虚拟的视角、不同色度或其他处理方式来识别微构造地貌、划分地貌面和确定断层位置等,宏观上获得与前人通过航片解译和野外调查一致的断裂分布特征,微观上较前者具有更高的精细程度。此外,数据采集、数据质量检验、数据处理及数据应用等技术和方法适用于其他能够获得LiDAR地形数据的活动断裂研究工作。
以甘孜-玉树断裂带东南段的地质地貌为研究对象,在遥感解译的基础上,通过对典型地区的详细野外调查和探槽研究对该段晚第四纪活动性进行研究。在断裂沿线的生康乡、仁果乡、错阿乡、日阿乡进行了断错地貌分析和晚第四纪滑动速率计算, 生康区的水平滑动速率为(7.6±0.5)mm/a, 垂直滑动速率为(1.1±0.1)mm/a; 仁果区的水平滑动速率为(8.0±0.3)mm/a,垂直滑动速率为(1.1±0.1)mm/a; 错阿区的水平滑动速率为(10.3±0.4)mm/a; 日阿区的水平滑动速率为(10.8±0.8)mm/a, 垂直滑动速率为(1.1±0.1)mm/a。在仁果乡和错阿乡进行了探槽研究,两处探槽都揭示了多次古地震事件,虽然揭露的断层构造样式有所不同,但总体上都是以走滑为主兼有一定的逆冲分量。综合古地震事件和滑动速率分析表明,甘孜-玉树断裂带东南段晚第四纪尤其是全新世以来活动剧烈。