Application of 3D Laser Scanning and Quantitative Morphology Analysis Method to Bedrock Fault Surface in Paleo-seismic Research in Bedrock Area
-
摘要: 由于对第四纪地层的严重依赖，传统古地震探槽研究方法在基岩区难以发挥作用，导致无法获取基岩区断层的强震活动历史。本研究以山西地堑系的交城断裂为目标断裂，以断裂北段2处基岩断层面为研究对象，通过三维激光扫描技术获取基岩断层面高精度形貌，基于变差函数法结合滑动窗口操作量化断层表面形貌特征，开展在基岩区提取断裂古地震信息的实例研究。结果显示，2处基岩断层面的形貌在高度上具有明显的分段特征，指示了断层面在地震事件作用下的分段出露过程。这种断层面形貌分段特征可以用来识别古地震事件和同震位移量。在思西村基岩断层面上，识别出由老到新的3次古地震事件，同震倾滑位移量依次为2.0 m、1.9 m和2.3 m，在上兰镇基岩断层面上，识别出由老到新的3次古地震事件，同震倾滑位移量依次为1.4 m、2.5 m和2.0 m，指示了交城断裂北、中段具有产生同震位移量大于2 m、震级大于7.5级的破裂型地震的能力。上述研究成果表明，基于三维激光扫描和形貌量化分析方法开展基岩断层面古地震研究，可以准确而高效地识别古地震事件次数和同震位移量，扩展古地震的研究对象，拓宽古地震的研究空间。在未来的研究中，可以适时地开展宇宙成因核素测年以测定断层面的暴露年龄，获得发震年代，给予地震序列年龄框架。Abstract: Due to the heavy reliance on Quaternary strata, traditional trench research methods are difficult to apply in the bedrock area, resulting in the inability to obtain strong seismic activity history on faults in the bedrock area. In this study, the Jiaocheng Fault in the Shanxi Rift is selected as the target fault, and the two bedrock fault surfaces in the northern segment are taken as the study object to carry out case studies. The variogram combined with the sliding window operation quantifies the morphologic features of the fault surface and extracts the paleo-seismic information. The results show that the morphology of the two bedrock fault surfaces has segmentation in height, indicating the exposure process of the fault surfaces under the action of periodic seismic events. Therefore, paleo-earthquake events and co-seismic slips can be identified by the morphological segments. Accordingly, three paleo-earthquake events were identified on the bedrock fault of Sixi village, and the co-seismic dip slips were 2.0 m, 1.9 m and 2.3 m in sequence. The co-seismic dip slips of the three paleoearthquakes are 1.4m, 2.5 m and 2.0 m in sequence. The results indicate that the northern segment of the northern and middle segments of Jiaocheng fault can produce ruptures with a co-seismic slip greater than 2 m and a magnitude greater than M7.5. The above results also show that the paleo-seismic research on bedrock fault surfaces based on t-LiDAR and morphological quantitative analysis can accurately and efficiently identify the number of paleo-earthquake and co-seismic slips, increase the objects of paleo-earthquakes, and expand research space of paleo-earthquakes. In future research, cosmogenic nuclide dating can be carried out timely to determine the exposure age of the fault surface, obtain the earthquake ages, and provide an age frame to the earthquake sequences.1) 2为区分变量，此处用D0表示同震位移，与刘静等（1996）中同震位移D含义一致。
-
引言
在活动断裂研究中，许多地区处于基岩区。目前最常用的探槽技术难以在基岩区开展断层的古地震研究工作（McCalpin，2009）。如何识别基岩区古地震事件、获取同震位移量是当前活动断裂研究工作中一项极为重要的问题（Goodall等，2021）。基岩断层面在一次次地震事件作用下会分段出露，由于暴露时间的先后差异，会在断层表面形成差异性的形貌风化条带（Giaccio等，2003）。因此，通过开展高精度基岩断层面形貌的定量分析，识别断层表面的形貌风化条带，是解决这一问题的有效方法（Wiatr等，2015）。三维激光扫描仪是一种高精度形貌扫描测量仪器，可以快速获取基岩断层表面精细形貌数据（Jones等，2009）。基于分形理论的变差函数法与滑动窗口运算结合，可以有效地量化断层表面的形貌特征（He等，2016）。将三维激光扫描技术与断层面形貌分析方法结合，可能成为基岩区开展古地震研究的有效途径（Zou等，2021；Giaccio等，2003）。
通常，一条断层在空间上会交替展布于基岩区和沉积区。传统的古地震研究方法主要是通过探槽技术分析第四纪沉积物的断错特征和年代（Sieh，1978，1984；Young等，2002；Galli等，2008；Ran等，2010；王林等，2017），其研究对象主要是有第四系覆盖的沉积区断层，重点关注断错的第四纪地层，而对展布于基岩区的断层，仅记录其空间展布和倾角、倾向等几何特征，没有充分揭示其所蕴含的地震活动信息。这在一定程度上限制了活动断裂的研究对象和空间尺度，导致难以获取断裂完整的地震活动历史（Zreda和Noll，1998；Palumbo等，2004）。特别是，基岩断层面是一种能够长期保存的活动构造地貌，其抗侵蚀能力强且不存在沉积间断等问题，能够连续的记录古地震事件，借助于断层面形貌定量分析方法，可以快速地获得较精准的同震位移量（Zou等，2020）。现阶段，有必要将基于三维激光扫描技术的断层面古地震研究与沉积区探槽结果结合，获取断裂更加完整的古地震历史，更好地评价断裂的活动特征。前人在-山西地堑系南部的霍山断裂和罗云山断裂上，开展了基岩断层面形貌分析的试验性工作，获得了断层面蕴含的古地震信息（He等，2016；Zou等，2020）。本文选用美国Trimble公司生产的三维激光扫描仪，对位于山西地堑系交城断裂北段的2处基岩断层面开展形貌扫描与定量分析工作，以识别古地震事件的期次和同震位移量，探索三维激光扫描技术与断层面形貌定量分析方法在基岩区古地震信息提取中的应用，丰富基岩断层面古地震的研究实例。
1. 断层面形貌定量分析法
三维激光扫描仪，也称地面光探测和测距仪器（陆基LiDAR），是一种高精度表面形貌测量仪器，作为一种非接触式测量系统，其已在地质灾害检测、同震变形的微地貌测量和震后建筑物破坏形式调查等方面取得了广泛应用（李峰等，2008；Jones等，2009）。该仪器的基本原理是，通过发射激光束至物体表面并接收其返回信号，测量扫描仪器中心与被测物体之间的距离，从而实现物体表面几何形貌的精细测量（Wiatr等，2015）。
本文基于陆基LiDAR扫描获取断层表面形貌的具体方法如下（图1）：（1）使用型号为Trimble GX 3-D的陆基LiDAR对断层表面进行激光测量，获取断层表面高精度、高分辨率形貌数据。相邻扫描点之间的距离随扫描仪与断层面之间的距离变化，最远点距扫描仪中心距离控制在15 m内，以保证获取的断层面形貌点云数据分辨率达2 mm。通过对保存完好、没有植被和沉积物覆盖的断层表面区域进行扫描，可获得断层表面的原始三维点云数据；（2）在扫描点相对位置不变的情况下，使用Realworks Survey Advanced 6.1软件移动、旋转具有三维空间坐标（x，y，z）的点云数据，将原来倾斜的断层面点云数据整体平整，变换后形貌数据的x轴为断层面走向，y轴为断面倾向，z轴为断层表面起伏方向；（3）使用Surfer 12软件，采用最大自然领域法，将已置平但不规则的点云数据插值为单元格大小为2 mm×2 mm的DEM数据；（4）利用Global Mapper 17软件将断层表面的DEM数据裁剪成规则的矩形，以开展断层面形貌的定量分析工作。数据采集的更多细节，参见前人在山西霍山和罗云山地区开展的断层面形貌分析工作（He等，2016；Zou等，2020）。
图 1 基岩断层面高精度形貌的扫描与处理Figure 1. Scanning and processing of high-precision morphology of bedrock fault surfaces自然的断层表面具有自相似或自相仿的形态特征，可以采用分形方法研究，并量化表达为分形维数D（Power等，1991；Sagy等，2009）。本研究通过各向同性变差函数与滑动窗口运算相结合，计算断层面的分形维数值D。变差函数法是处理DEM高程场的一种十分有效的方法（Chase，1992；Davies等，1999），其本质是描述高程场的平均差异r (t) 随点距t的变化关系：
$$ r(t)=c t^{K} $$ (1) 式中，c为常数。对于二维面状分形体，分形维数 D 与分形指数 K具有以下关系（Bi等，2012）：
$$ D=\dfrac{(3-K)}{2} $$ (2) 式中，分形维数 D 代表自然表面的复杂度（Klinkenberg，1992；Sung等，2004），理论值在2～3之间，是反映表面形貌特征的定量化参数（Mandelbrot，1967；Elliot，1989）。
根据Sung等（1998）提出的准则，选择尺寸大小为 66 mm×66 mm、130 mm×130 mm 和 258 mm×258 mm 的滑动窗口，采用变差函数法对基岩断层表面的分形维数值进行计算（图2）。窗口尺寸为W，沿断层表面每次滑动W/2，由左至右、由上至下遍历整个DEM 数据场，获得分形维数值在基岩断层表面的栅格分布图。为更好地显示断层表面形貌沿高度的变化，以行作为单元计算同一高度上分形维数的平均值和标准差，获得分形维数值的散点变化图。最后，采用Student’s t-test检验法，对沿高度方向上的分形维数值进行检验分段，划分形貌条带（Zou等，2020）。
图 2 基岩断层面高精度形貌的量化分析Figure 2. Quantitative analysis of high-precision morphology of bedrock fault surfaces埋藏于地下的基岩断层表面主要受构造滑动控制，具有初始的形态（Sagy等，2007）。一旦断层面在地震事件作用下出露于地表，风化作用则成为影响其表面形貌的主要因素。区域风化速率基本一致的情况下，断层表面的风化形貌特征主要与出露时间相关（Wiatr等，2015）。在一次次地震事件作用下，断层面会分段出露，先、后出露的断层表面因经历不同的暴露时长而显示出不同的形貌风化特征，从而沿高度方向上断层表面形貌呈现分段特征（He等，2016）。这种形貌特征一旦被定量表达，则可以用于恢复基岩断层面的出露过程，识别古地震期次和同震位移量（图3；Zou等，2020）。
图 3 基岩断层面的出露过程和相应的形貌特征Figure 3. The exposure process of the bedrock fault surface and the corresponding morphological characteristics2. 研究实例
山西地堑系中部的交城断裂，位于吕梁山山前，是太原盆地的活动边界断裂。断裂北起泥屯镇向西南延伸至汾阳市，总长约125 km，整体走向NE-SW，倾向SE，南北端分别被石岭关隆起和灵石隆起截止（谢新生等，2008）。根据其几何形状和活动性，断裂分为北、中、南3段，南段无全新世活动证据。先前的古地震探槽研究揭示了断裂北段和中段的3次全新世古地震事件，2段在发震时间上具有同步性（图4；谢新生等，2008；郭慧等，2012）。其中，断裂北段多处出露保存完好的基岩断层面为开展陆基LiDAR扫描和断层面定量形貌分析，提供了良好的研究对象。本研究选择位于北段的思西村和上兰镇2处基岩断层面开展研究，2个断层面相距约22公里，岩性均为灰岩，研究点附近均没有分支断层，保证了记录的古地震事件的完整性和同震位移量的准确性。
图 4 目标断裂与研究点的位置和基本信息Figure 4. Location and basic information of the aimed fault and study points2.1 交城断裂—思西村
思西村基岩断层整体走向NNE-SSW，地表延伸数百米，断层面倾向SEE，倾角60°。选择的断层面扫描区域远离建筑物等人为活动痕迹，并与山麓发育的小沟壑保持有一定距离（图5（a））。野外调查显示基岩断层面呈现出3个明显的风化条带：（1）断层面下部（风化条带3），基岩断层面相对完整，表面可观察到微小裂隙；（2）断层面中部（风化条带2），基岩断层面开始呈块状，表面伴有侵蚀坑和较大的裂缝；（3）断层面上部（风化条带1），基岩断层面被茂密的植被覆盖，表面严重破碎（图5（b））。在整个断层面上，均可观察到细小地衣残体的不均匀分布特征，未观察到碳酸岩层理。使用陆基LiDAR 扫描获取三维点云数据，并按前文所述流程处理生成2 mm×2 mm的DEM数据（图5（c）），开展进一步的量化分析。
采用变差函数法结合滑动窗口操作，对思西村基岩断层面的高精度形貌数据进行定量化分析。定量化形貌分形结果呈现出明显的3分段特性，且66 mm×66 mm、130 mm×130 mm 和 258 mm×258 mm 3个滑动窗口尺寸下的分形结果表现出较好的一致性。沿着断层面高度方向上，在2.3～2.8 m和4.7～5 m处，断层面分形维数存在2个缓变区间，这2个缓变区间将断层表面形貌分成3段。段1位于断层面上部，分段高度H上= 2.0 m，分形维数D上=2.3218±0.0485（66 mm）、D上=2.2979±0.0458（130 mm）、D上=2.2842±0.0339（258 mm）。段2位于断层面中部，分段高度H中= 1.9 m，分形维数D中=2.4118±0.331（66 mm）、D中=2.3876±0.0262（130 mm）、D中=2.3771±0.0361（258 mm）。段3位于断层面下部，分段高度H下= 2.3 m，分形维数D下=2.4303±0.0506（66 mm）、D下=2.4107±0.0500（130 mm）、D下=2.3889±0.0441（258 mm）（图6（a）、表1）。
表 1 思西村基岩断层面的分形维数与分段高度Table 1. Characteristic fractal value （D）and heights of morphological segments on Sixi village bedrock fault surface断层面 滑动窗口尺寸 分形维数（D） 标准差（σ） 下段 中段 上段 下段 中段 上段 思西村
断层面66 mm×66 mm 2.4303 2.4118 2.3218 0.0506 0.0331 0.0485 130 mm×130 mm 2.4107 2.3876 2.2979 0.0500 0.0262 0.0458 258 mm×258 mm 2.3889 2.3771 2.2842 0.0441 0.0361 0.0339 分段高度H/m / 2.3 1.9 2.0 / / / 
图 6 思西村基岩断层面形貌的定量分析和古地震事件识别结果Figure 6. Quantitative morphology analysis and paleo-earthquake identification results for Sixi village bedrock fault surface基岩断层面的分段形貌特征指示了断层面在周期性地震事件作用下分带出露的过程。先出露的断层面和后出露的断层面，经历不同的风化时间会形成形貌差异条带。因此，可以根据形貌分段的数量和高度，识别古地震事件的期次和同震位移量。据此，思西村基岩断层面的3个形貌分段（H上、H中和H下）可能指示了3次古地震事件（E1、E2 和 E3），同震倾滑位移量分别为 2.0 m、1.9 m 和 2.3 m，根据实测的断层面倾角60°，垂直同震滑动位移量分别为1.74 m、1.65 m和1.99 m。断层破裂长度L和同震位移D0是估计古地震震级的常用参数（McCalpin，2009）。然而，准确的古地震地表破裂长度通常不易获得，同震位移量不仅反映了地震能量的释放，而且与地震的震级呈正相关（Wells等，1994）。本研究中通过基岩断层的形貌分析获得的同震位移量可用于估计震级。根据刘静等（1996）提出的华北正断层震级M和位移D0之间的经验公式，估算3次古地震事件的震级分别为MS7.5、MS7.5和MS7.6（图6（b））。探槽研究工作揭示的断裂北段和中段（总长105 km）全新世以来具有同步地震破裂性，使用华北正断层的震级M和位移D0、长度L之间的经验公式（刘静等，1996），估算得3次古地震事件的震级分别为MS7.8、MS7.8和MS7.9。
2.2 交城断裂—上兰镇
上兰镇基岩断层位于交城断裂北段的中部，断层坎高约10 m，倾向SE，倾角70°。基岩断层坎7.5 m以上部分逐渐后退，7.5 m以下保存完好且可以进行扫描用于古地震事件分析（图7（a））。上兰镇调查点远离河流和洪积扇等地貌，因此断层坎底部没有相关的侵蚀破坏或沉积掩埋，确保了断层面的完整性。选择形貌扫描区域时，避开受人为活动影响（如植被）的基岩断面，以保证采集的形貌数据的可用性。由于风化程度和暴露时间不同的影响，在断层面不同高度处观察到3个明显的形貌特征（图7（b））：（1）断层面下部（风化条带3）相对完整，表面可观察到小的张性裂缝和擦痕；（2）断层面中部（风化条带2）表面形貌开始变得粗糙，可见凹坑；（3）断层面上部（风化条带1）呈现出被风化碎屑覆盖的特征。同思西村一样，使用 陆基LiDAR 扫描保存完好的区域并获得点云数据，通过插值得到单元大小为2 mm×2 mm的DEM数据（图7（c））。值得注意的是，在上兰镇基岩断层面未发现地衣残体和碳酸岩层理，断面底部观测到擦痕（图7（d）），擦痕方向近垂直，指示了断层面受构造活动出露、具有正断运动的特征。在断层表面观测到局部存在张性裂隙，这可能会对断层面形貌的量化分析造成一定影响，但由于是局部现象，可通过对同一高度方向上分形维数取拟合均值的方法进行消除。
采用变差函数法结合滑动窗口操作，对上兰镇基岩断层面的高精度形貌数据进行定量化分析。上兰镇基岩断层面形貌在3种滑动窗口尺寸下同样呈现出3分段特征。沿断层面高度方向上，在2～3 m和5.5～6.1 m处，断层面分形维数存在2两个缓变区间，这2两个缓变区间将断层表面形貌分成3段。段1位于断层面上部，分段高度H上= 1.4 m，分形维数D上=2.2361±0.0419（66 mm）、D上=2.2122±0.0458（130 mm）、D上=2.1685±0.0174（258 mm）。段2位于断层面中部，分段高度H中= 2.5m，分形维数D中=2.5331±0.815（66 mm）、D中=2.5204±0.0653（130 mm）、D中=2.5025±0.0502（258 mm）。段3位于断层面下部，分段高度H下= 2.0 m，分形维数D下=2.7105±0.0266（66 mm）、D下=2.6810±0.0441（130 mm）、D下=2.6596±0.0165（258 mm）（图8（a）、表2）。值得注意的是，在张裂隙发育的地方，断层表面的分形维数值会呈现出低值情况，造成分形维数的波动，但是这种微小的波动不构成独立的形貌分段，对形貌分段结果没有实质影响，仅反映了断层面的局部形貌特征。
表 2 上兰镇基岩断层面的分形维数和分段高度Table 2. Characteristic fractal value （D） and heights of morphological segments on Shanglanzhen bedrock fault surface断层面 滑动窗口尺寸 分形维数（D） 标准差（σ） 下段 中段 上段 下段 中段 上段 上兰镇
断层面66 mm×66 mm 2.7105 2.5331 2.2361 0.0266 0.0815 0.0419 130 mm ×130 mm 2.6810 2.5204 2.2122 0.0441 0.0653 0.0458 258 mm×258 mm 2.6596 2.5025 2.1685 0.0165 0.0502 0.0174 分段高度H/m / 2.0 2.5 1.4 / / / 
图 8 上兰镇基岩断层面形貌定量分析和古地震事件识别结果Figure 8. Quantitative morphology and paleo-earthquake identification results of Shanglanzhen bedrock fault surface上兰镇基岩断层面的分段形貌特征指示了断层面的分带出露过程，这种分带出露过程与周期性的地震事件有关。因此，可以根据形貌分段的数量和高度来识别古地震事件的期次和同震位移量。据此，上兰镇基岩断层面的3个形貌分段（H上、H中和H下）可能指示了3次古地震事件（E1、E2 和 E3），同震倾滑位移量分别为 1.4 m、2.5 m 和 2.0 m，根据实测的断层面倾角70°，垂直同震滑动位移量分别为1.32 m、2.35 m和1.88 m。根据华北正断层的震级M和同震位移D0之间的经验公式（刘静等，1996），估算3次古地震事件的震级分别为MS7.4、MS7.6和MS7.5。在上兰镇基岩断层剖面中，7.5 m处存在明显的坡折，在此高度以上断层面开始后退，这可能指示了较早的地震事件，但是由于时间久远难以通过形貌学方法识别（图8（b））。根据本研究获得的3次古地震事件的同震位移量D0，结合交城断裂北、中段的断裂长度，使用华北正断层的震级M和位移D0、长度L之间的经验公式（刘静等，1996），估算得3次古地震事件的震级分别在MS7.7、MS7.9和MS7.8。
上兰镇和思西村基岩断层面的形貌分形结果在分段性上是一致的，均呈现出3分段特征，且分段特征均通过了统计学方法Student’s t-test的检验（图9）。根据形貌分段结果，思西村基岩断层面形貌特征揭示出由老到新的3次古地震事件，同震倾滑位移量依次为2.0 m、1.9 m和2.3 m；上兰镇基岩断层面揭示出由老到新的3次古地震事件，同震倾滑位移量依次为1.4 m、2.5 m和2.0 m。但二者形貌特征在细节上存在一定的差异（图8），主要表现为以下2点：一是在相邻2个形貌分段的过渡带宽度上，思西村断层面过渡带较窄，而上兰镇断层面过渡带较宽，这可能指示了相较于思西村断层面，上兰镇基岩断层面在间震期受弱水动力条件的影响程度更强；二是在分段内的分形维数集中度上，思西村断层面较上兰镇断层面更加离散，显示出段内形貌的不均一性，这可能与前文所述的思西村断层面上存在随机分布的碳化地衣残体有关。
图 9 思西村和上兰镇调查点的Student’s t-test分段检验结果Figure 9. Morphology segmentation results of Sixi village and Shanglanzhen survey points based on Student's t-test method谢新生等（2008）通过对4个探槽的研究确定了交城断裂北段和中段经历了3次全新世的古地震事件，相应的地震事件年龄为～8.36 ka、5.32～6.1 ka和3.06～3.53 ka。郭慧等（2012）通过对3个探槽的研究得出了类似结论，揭示了交城断裂北段和中段3次全新世的古地震事件，确定地震事件的年龄分别为8.53～8.56 ka、～5.91 ka、3.06～3.74 ka。因此，交城断裂北、中段的全新世地震活动在期次和年代范围方面得到了较好确定，根据毛凤英等（1995）提出的事件窗口分析法，将3次古地震事件的年龄分别限定为8.36～8.56 ka、5.32～6.14 ka和3.06～3.74 ka（图10）。
图 10 交城断裂简图与全新世以来古地震信息Figure 10. Schematic map and paleoearthquake information of the JCF since the Holocene本研究通过基岩断层面的定量形貌分析，揭示了3次古地震事件及对应的同震位移量。从地震事件的数量（3次）和地震强度（M＞7.5的破裂型地震）来看，本研究与已有的探槽古地震研究具有较好的一致性。2个基岩断层露头附近（2 km以内）都有前人开挖的探槽，且2个基岩断层点与开挖的探槽位于同一断层段上，附近没有分支断层迹象。基于上述分析，本研究揭示的3次古地震可能指示了全新世的3次古地震事件，但更为直接的证据需要通过绝对定年手段（如宇宙成因核素36Cl暴露测年法）予以进一步确认。
将本研究结果与已有研究相结合，可以获取更为完整的交城断裂古地震信息，有助于更为合理的区域地震危险性评价。考虑到交城断裂北段和中段，全新世以来发生了复发间隔在～2.6 kyr的3次破裂型地震事件，具有超过3 kyr的离逝时间，能够产生＞2 m的同震位移，具备震级M＞7.5的发震能力，应对交城断裂的北、中段活动性和周边地区的地震危险性予以密切关注。
3. 结论与讨论
陆基LiDAR扫描与高精度断层面形貌定量分析方法，为基岩区开展古地震研究工作提供了新思路。以山西地堑系交城断裂北段的2处基岩断层面为研究实例，开展了断层面形貌学定量分析。结果表明，2处基岩断层面形貌在高度上均具有分段特征，指示了断层面在地震事件作用下分带出露的过程。通过识别断层表面形貌分段和分段个数来确定古地震事件和相应的同震位移量，在思西村基岩断层面识别出由老到新的3次古地震事件，同震倾滑位移量依次为2.0 m、1.9 m和2.3 m，在上兰镇基岩断层面识别出由老到新的3次古地震事件，同震倾滑位移量依次为1.4 m、2.5 m和2.0 m。基岩断层面形貌学结果结合已有探槽研究显示，交城断裂北、中段能够产生同震位移量＞2 m、震级M＞7.5的破裂型地震，对交城断裂北、中段的地震活动性和周边地区的地震危险性应当予以密切关注。
研究表明，定量的形貌学分析方法应用于基岩断层面，可以准确、高效地识别古地震事件次数和同震位移量，扩展古地震研究对象，拓宽古地震研究空间。值得关注的是，运用基岩断层面形貌学方法开展基岩区活动断裂研究，需要完整地获取断层表面的精细形貌特征，因此对于断层表面的完整性要求较高。为此，选址时需要通过野外调查，保证研究点远离沟谷侵蚀域，避开断层坎顶部的破碎带和剥蚀坡，选择基岩断层下部完整的部分进行研究。同时，需要注意所选断层面的自然属性，选择的基岩露头最好远离农田、大坝等人为建筑，尽可能排除人类活动干扰。并不是每个沿断层走向的基岩面露头点都适合开展断层面形貌分析工作，因此，本研究方法只能获得典型点位的同震位移测量，难以获取沿断裂走向上的活动性变化信息。
本研究通过断层擦痕等构造形迹观测和断错地貌特征分析，论证了基岩断层面的构造属性，基于高精度形貌量化分析技术识别了断层面的地震活动历史。需要注意的是，这一研究仅对断层表面的形貌特征进行了静态分段，其研究基础在于形貌分段之间暴露时间的差异。气候变化会造成基岩断层表面的风化速率改变，进而对断层表面形貌的变化快慢产生影响。因此，未来的工作中，有必要考虑山西地区全新世以来的风化速率这一变化量，开展断层面形貌随暴露时间和风化速率动态演化过程的研究。此外，没有对基岩断层面露头开展绝对定年工作，缺乏地震事件的年代框架也是本研究的局限。今后的研究中，应适时地开展 36Cl 等宇宙成因核素测年以获得绝对年龄，限定地震事件的发震年代，以更好地与沉积区探槽古地震的研究成果综合对比。这是将来的研究中需要着力解决的问题，也是我们未来工作的重点。
致谢 感谢审稿专家和期刊编辑的宝贵修改意见！
-
图 1 基岩断层面高精度形貌的扫描与处理
Figure 1. Scanning and processing of high-precision morphology of bedrock fault surfaces
图 2 基岩断层面高精度形貌的量化分析
Figure 2. Quantitative analysis of high-precision morphology of bedrock fault surfaces
图 3 基岩断层面的出露过程和相应的形貌特征
Figure 3. The exposure process of the bedrock fault surface and the corresponding morphological characteristics
图 4 目标断裂与研究点的位置和基本信息
Figure 4. Location and basic information of the aimed fault and study points
图 6 思西村基岩断层面形貌的定量分析和古地震事件识别结果
Figure 6. Quantitative morphology analysis and paleo-earthquake identification results for Sixi village bedrock fault surface
图 8 上兰镇基岩断层面形貌定量分析和古地震事件识别结果
Figure 8. Quantitative morphology and paleo-earthquake identification results of Shanglanzhen bedrock fault surface
图 9 思西村和上兰镇调查点的Student’s t-test分段检验结果
Figure 9. Morphology segmentation results of Sixi village and Shanglanzhen survey points based on Student's t-test method
图 10 交城断裂简图与全新世以来古地震信息
Figure 10. Schematic map and paleoearthquake information of the JCF since the Holocene
表 1 思西村基岩断层面的分形维数与分段高度
Table 1. Characteristic fractal value （D）and heights of morphological segments on Sixi village bedrock fault surface
断层面 滑动窗口尺寸 分形维数（D） 标准差（σ） 下段 中段 上段 下段 中段 上段 思西村
断层面66 mm×66 mm 2.4303 2.4118 2.3218 0.0506 0.0331 0.0485 130 mm×130 mm 2.4107 2.3876 2.2979 0.0500 0.0262 0.0458 258 mm×258 mm 2.3889 2.3771 2.2842 0.0441 0.0361 0.0339 分段高度H/m / 2.3 1.9 2.0 / / / 
下载: 导出CSV
表 2 上兰镇基岩断层面的分形维数和分段高度
Table 2. Characteristic fractal value （D） and heights of morphological segments on Shanglanzhen bedrock fault surface
断层面 滑动窗口尺寸 分形维数（D） 标准差（σ） 下段 中段 上段 下段 中段 上段 上兰镇
断层面66 mm×66 mm 2.7105 2.5331 2.2361 0.0266 0.0815 0.0419 130 mm ×130 mm 2.6810 2.5204 2.2122 0.0441 0.0653 0.0458 258 mm×258 mm 2.6596 2.5025 2.1685 0.0165 0.0502 0.0174 分段高度H/m / 2.0 2.5 1.4 / / /
下载: 导出CSV
-
[1] 郭慧, 江娃利, 谢新生, 2012. 山西交城断裂北端及中段3个大型探槽全新世断错现象分析. 地震地质, 34（1）: 76—92 doi: 10.3969/j.issn.0253-4967.2012.01.008Guo H. , Jiang W. L. , Xie X. S. , 2012. Analysis of Holocene faulting phenomena revealed in the three trenches along the northern and central Jiaocheng Fault, Shanxi. Seismology and Geology, 34(1): 76—92. (in Chinese) doi: 10.3969/j.issn.0253-4967.2012.01.008 [2] 李峰, 徐锡伟, 陈桂华等, 2008. 高精度测量方法在汶川MS 8.0地震地表破裂带考察中的应用. 地震地质, 30（4）: 1065—1075 doi: 10.3969/j.issn.0253-4967.2008.04.022Li F. , Xu X. W. , Chen G. H. , et al. , 2008. The application of different high-precision survey methods in the investigation of the MS 8.0 Wenchuan earthquake surface ruptures. Seismology and Geology, 30(4): 1065—1075. (in Chinese) doi: 10.3969/j.issn.0253-4967.2008.04.022 [3] 刘静, 汪良谋, 1996. 震级与地震破裂参数统计关系式系数的讨论. 地震地质, 18（3）: 225—236Liu J. , Wang L. M. , 1996. Discussion on coefficients of statistical regressions between magnitude and earthquake rupture parameters. Seismology and Geology, 18(3): 225—236. (in Chinese) [4] 毛凤英, 张培震, 1995. 古地震研究的逐次限定方法与新疆北部主要断裂带的古地震研究. 活动断裂研究, （4）: 153—164Mao F. Y. , Zhang P. Z. , 1995. Progressive constraining method in paleoseismic study and paleoearthquakes along the major active faults in northern Xinjiang. Research on Active Fault, (4): 153—164. (in Chinese) [5] 王林, 田勤俭, 李德文, 等, 2017. 蔚广盆地南缘断裂带唐山口段山前断层活动性研究. 震灾防御技术, 12（2）: 319—328 doi: 10.11899/zzfy20170208Wang L. , Tian Q. J. , Li D. W. , et al. , 2017. The activity of the piedmont fault along the Tangshankou segment of Yuguang basin southern marginal fault. Technology for Earthquake Disaster Prevention, 12(2): 319—328. (in Chinese) doi: 10.11899/zzfy20170208 [6] 谢新生, 江娃利, 孙昌斌等, 2008. 山西交城断裂带多个大探槽全新世古地震活动对比研究. 地震地质, 30（2）: 412—430 doi: 10.3969/j.issn.0253-4967.2008.02.007Xie X. S. , Jiang W. L. , Sun C. B. , et al. , 2008. Comparison study on Holocene paleoseismic activities among multi-trenches along the Jiaocheng Fault zone, Shanxi. Seismology and Geology, 30(2): 412—430. (in Chinese) doi: 10.3969/j.issn.0253-4967.2008.02.007 [7] Bi L. S., He H. L., Wei Z. Y., et al., 2012. Fractal properties of landforms in the ordos block and surrounding areas, China. Geomorphology, 175—176: 151—162. [8] Chase C. G. , 1992. Fluvial landsculpting and the fractal dimension of topography. Geomorphology, 5(1—2): 39—57. doi: 10.1016/0169-555X(92)90057-U [9] Davies S. , Hall P. , 1999. Fractal analysis of surface roughness by using spatial data. Journal of the Royal Statistical Society: Series B (Statistical Methodology), 61(1): 3—37. doi: 10.1111/1467-9868.00160 [10] Elliot J. K. , 1989. An investigation of the change in surface roughness through time on the foreland of Austre Okstindbreen, North Norway. Computers & Geosciences, 15(2): 209—217. [11] Galli P. , Galadini F. , Pantosti D. , 2008. Twenty years of paleoseismology in Italy. Earth-Science Reviews, 88(1—2): 89—117. doi: 10.1016/j.earscirev.2008.01.001 [12] Giaccio B. , Galadini F. , Sposato A. , et al. , 2003. Image processing and roughness analysis of exposed bedrock fault planes as a tool for paleoseismological analysis: results from the Campo Felice fault (central Apennines, Italy). Geomorphology, 49(3—4): 281—301. doi: 10.1016/S0169-555X(02)00191-5 [13] Goodall H. J. , Gregory L. C. , Wedmore L. N. J. , et al. , 2021. Determining histories of slip on normal faults with bedrock scarps using cosmogenic nuclide exposure data. Tectonics, 40(3): e2020TC006457. [14] He H. L. , Wei Z. Y. , Densmore A. , 2016. Quantitative morphology of bedrock fault surfaces and identification of paleo-earthquakes. Tectonophysics, 693: 22—31. doi: 10.1016/j.tecto.2016.09.032 [15] Jones R. R. , Kokkalas S. , McCaffrey K. J. W. , 2009. Quantitative analysis and visualization of nonplanar fault surfaces using terrestrial laser scanning (LiDAR)—The Arkitsa fault, Central Greece, as a case study. Geosphere, 5(6): 465—482. [16] Klinkenberg B., 1992. Fractals and morphometric measures: is there a relationship? Geomorphology, 5(1—2): 5—20. [17] Mandelbrot B. , 1967. How long is the coast of Britain? Statistical self-similarity and fractional dimension. Science, 156(3775): 636—638. doi: 10.1126/science.156.3775.636 [18] McCalpin J., 2009. Paleoseismology. 2nd ed. Burlington: Academic Press. [19] Palumbo L. , Benedetti L. , Bourlès D. , et al. , 2004. Slip history of the Magnola fault (Apennines, Central Italy) from 36Cl surface exposure dating: evidence for strong earthquakes over the Holocene. Earth and Planetary Science Letters, 225(1—2): 163—176. doi: 10.1016/j.jpgl.2004.06.012 [20] Power W. L. , Tullis T. E. , 1991. Euclidean and fractal models for the description of rock surface roughness. Journal of Geophysical Research: Solid Earth, 96(B1): 415—424. doi: 10.1029/90JB02107 [21] Ran Y. K. , Chen L. C. , Chen J. , et al. , 2010. Paleoseismic evidence and repeat time of large earthquakes at three sites along the Longmenshan fault zone. Tectonophysics, 491(1—4): 141—153. doi: 10.1016/j.tecto.2010.01.009 [22] Sagy A. , Brodsky E. E. , Axen G. J. , 2007. Evolution of fault-surface roughness with slip. Geology, 35(3): 283—286. doi: 10.1130/G23235A.1 [23] Sagy A. , Brodsky E. E. , 2009. Geometric and rheological asperities in an exposed fault zone. Journal of Geophysical Research: Solid Earth, 114(B2): B02301. [24] Sieh K. E. , 1978. Prehistoric large earthquakes produced by slip on the San Andreas Fault at Pallett Creek, California. Journal of Geophysical Research: Solid Earth, 83(B8): 3907—3939. doi: 10.1029/JB083iB08p03907 [25] Sieh K. E. , 1984. Lateral offsets and revised dates of large prehistoric earthquakes at Pallett Creek, southern California. Journal of Geophysical Research: Solid Earth, 89(B9): 7641—7670. doi: 10.1029/JB089iB09p07641 [26] Sung Q. C., Chen Y. C., Chao P. C., 1998. Spatial variation of fractal parameters and its geological implications. Terrestrial, Atmospheric and Oceanic Sciences, 9(4): 655—672. [27] Sung Q. C. , Chen Y. C. , 2004. Self-affinity dimensions of topography and its implications in morphotectonics: an example from Taiwan. Geomorphology, 62(3—4): 181—198. doi: 10.1016/j.geomorph.2004.02.012 [28] Wells D. L. , Coppersmith K. J. , 1994. New empirical relationships among magnitude, rupture length, rupture width, rupture area, and surface displacement. Bulletin of the Seismological Society of America, 84(4): 974—1002. [29] Wiatr T. , Papanikolaou I. , Fernández-Steeger T. , et al. , 2015. Bedrock fault scarp history: insight from t-LiDAR backscatter behaviour and analysis of structure changes. Geomorphology, 228: 421—431. doi: 10.1016/j.geomorph.2014.09.021 [30] Young J. J. , Arrowsmith J. R. , Colini L. , et al. , 2002. Three-dimensional excavation and recent rupture history along the Cholame segment of the San Andreas fault. Bulletin of the Seismological Society of America, 92(7): 2670—2688. doi: 10.1785/0120000604 [31] Zou J. J. , He H. L. , Yokoyama Y. , et al. , 2020. Seismic history of a bedrock fault scarp using quantitative morphology together with multiple dating methods: a case study of the Luoyunshan piedmont fault, southwestern Shanxi Rift, China. Tectonophysics, 788: 228473. doi: 10.1016/j.tecto.2020.228473 [32] Zou J. J., He H. L., Yokoyama Y., et al., 2021. Identification of Paleoearthquakes and Coseismic Slips on a Normal Fault Using High-Precision Quantitative Morphology: application to the Jiaocheng Fault in the Shanxi Rift, China. Lithosphere, (Special 2): 2550879. [33] Zreda M. , Noll J. S. , 1998. Ages of prehistoric earthquakes revealed by cosmogenic chlorine-36 in a bedrock fault scarp at Hebgen Lake. Science, 282(5391): 1097—1099. doi: 10.1126/science.282.5391.1097 期刊类型引用(3)
1. 张军龙，徐岳仁，李文巧，陈峰. 中强地震发震构造标志浅析——以2023年积石山M_S6.2地震为例. 地震. 2024(01): 226-234 . 
百度学术2. 陈陆望，郑忻，张杰，张苗，胡永胜，郑剑，殷晓曦. 基于断层量化隐伏型煤矿基岩含水层水化学空间演化机制研究. 地质学报. 2024(07): 2245-2257 . 百度学术3. 高雅真. 面向高分辨率激光扫描点云的图像色彩界面交互配色设计研究. 激光杂志. 2024(08): 159-163 . 百度学术其他类型引用(1)
-
下载:
点击查看大图
计量
- 文章访问数: 350
- HTML全文浏览量: 70
- PDF下载量: 29
- 被引次数: 4
