Study on dynamic response characteristics of ground fissure sites in Songzhuang, Beijing
-
摘要:
地裂缝灾害对北京局部地区人居安全和城乡规划建设造成了较大威胁。为分析北京地裂缝活动对发育场地动力响应的影响特征、规律和范围，本文选取发育较典型的通州区宋庄地裂缝为研究对象，基于野外调查，得到了地裂缝现场发育特征和影响范围，并通过地脉动测试手段，获取了场地卓越频率等重要参数；利用3种动力分析方法，对地裂缝场地的动力响应特性和规律进行了研究。结果表明：北京宋庄地裂缝场地的卓越频率和卓越周期与测点选取的位置无明显关系；地裂缝场地有着显著的“上盘效应”，位于地裂缝上盘的傅里叶卓越频率、反应谱卓越周期和Arias烈度峰值总大于下盘，且在距离地裂缝较近的区域，地脉动测试幅值存在着较为突出的“放大效应”；地裂缝对场地动力响应的影响范围宽度为距地裂缝上下盘各12 m。研究结果可对北京宋庄地裂缝发育场地的规划建设和工程结构抗震设防提供一定的参考。
Abstract:Abstracts: Ground fissure disasters pose significant threats to the safety of residential areas and urban-rural planning in localized areas of Beijing. To analyze the impact of ground fissure activity on the dynamic response characteristics, patterns, and scope of affected sites, this study focuses on the Songzhuang ground fissure in Tongzhou district , a typical development area. Based on the field investigations, the study captured the development characteristics and influence scope of ground fissures and obtain key parameters such as the predominant frequency of the site through microtremor testing. Three dynamic analysis methods were applied to investigate the dynamic response characteristics and patterns of the ground fissure site. The results show that there is no obvious relationship between the excellence frequency and excellence period of the ground fissure site in Songzhuang, Beijing and the location of the measurement points. The ground fissure site exhibits a pronounced "hanging wall effect", with the Fourier Superior Frequency, the excellent cycle of the response spectra, and the peak intensity of Arias located in the hanging wall of the ground fissure generally exceeding those on the footwall. Additionally, the amplitude of the ground pulsation test is greater than that of the footwall in the area nearer to the ground fissure. The results of this study can provide valuable references for planning, construction, and seismic fortification of structures in the Songzhuang ground fissure development area.
-
Keywords:
- ground fissure /
- microtremor /
- dynamic response /
- amplification effect /
- extent of impact
-
0. 引 言
在三峡库区中有着众多失稳模式各异的危岩体，这些临江的危岩体一旦失稳破坏，可能引发灾害性涌浪，严重威胁航道和沿岸居民的生命财产安全，例如可能发生滑移式破坏箭穿洞危岩体[1]、可能发生压溃式破坏的龙门寨危岩体[2]、可能发生倾倒式破坏的棺木岭危岩体[3]、可能发生坠落式破坏的望霞危岩体[4]等。目前对滑坡产生涌浪灾害的研究较为丰富，但对于危岩体崩塌产生涌浪灾害研究零散、不系统，快速准确评估危岩体崩塌产生涌浪灾害对库区防灾减灾具有重要意义。
目前，包含滑坡涌浪灾害在内的涌浪研究中，主要通过公式法、概化物理试验模型法、缩尺物理模型试验法和数值模拟等方法开展[5]。如牟英华等[6]选取的美国土木工程学会建议的推算公式法对文昌阁危岩体的涌浪进行分析，代云霞等[7]采用潘家铮公式法和守恒法的理论计算方法对巫山县某崩塌体崩滑速度、入江体积、初始涌浪高度及涌浪的传播衰减、爬高进行计算，对比理论计算结果和调查数据；杨渠锋等 [8]采用概化物理实验模型法， 通过不同堆放组合下的矩形滑坡体, 模拟了三峡库区典型的陡岩崩塌发生、涌浪形成及传播的整个试验过程,得出了初始涌浪传播的波高衰减系数计算公式和陡岩崩塌涌浪爬高的计算公式;易臻彦[9]通过 FLOW-3D 进行了危岩崩塌涌浪数值模拟，得到了危岩入水运动过程和涌浪高度变化情况。表1总结了目前常用的涌浪评价方法的优缺点，其中公式法因其技术门槛低、使用便捷且经济的优势，在实际工程中应用广泛。
表 1 涌浪研究方法及优劣Table 1. Surge research methods and advantages and disadvantages方法 公式法 概化物理试验模型法 缩尺物理模型试验法 数值模拟法 准确度 中 中-高 高 高 时间 少 少 非常多 多 经费 低 低 非常高 中等 技术门槛 低 中等 高 高 由于试验规模和成本限制，大多涌浪公式研究是根据涌浪的三阶段，即涌浪产生、涌浪传播和涌浪爬高，对涌浪的全过程分阶段进行的研究，少有适用于涌浪全过程的公式计算体系。且大多研究关注滑坡涌浪，并少有考虑危岩体失稳模式不同所带来的影响。危岩体的失稳模式不同时，其涌浪特征也会不同，包括入水点的不同[10]、水体运动特征不同[1]、涌浪波类型不同[11]等。基于此，本文对现有的阶段性涌浪公式进行了系统化梳理，建立了考虑危岩体不同失稳模式的涌浪公式计算体系，并在此基础上编制了危岩涌浪计算引擎。以不同失稳模式的典型危岩体为例，进行了涌浪危险性分析，同时验证了公式计算体系的正确性。这对公式法在危岩体涌浪计算的运用具有重要意义，有利于库区崩塌涌浪防灾减灾。
1. 涌浪公式法计算体系
崩滑体涌浪灾害主要可分为三个阶段：涌浪产生、涌浪传播和涌浪爬高[12]。当崩滑体失稳入江时，所激起的涌浪，即为涌浪的产生阶段；如图1所示，在涌浪产生后，涌浪通常以环的形式向外传播，当抵达对岸后发生反射和干涉，并随着传播距离的增加，逐渐由环状传播转变为平行传播[13]；当涌浪传播至库岸时，则会在沿岸产生爬高。涌浪在产生和传播阶段时，主要影响航道安全，而爬高阶段则主要影响沿岸居民安全。
1.1 涌浪产生阶段
崩滑体的失稳模式将直接影响涌浪的产生阶段，危岩体的失稳模式主要有滑移式、坠落式、压溃式和倾倒式，他们产生的涌浪特征也有所不同，主要体现在①入水点的空间位置不同。入水点不同，形成的原始涌浪形态就有差异。从岸边激发的涌浪为新月形，而从离岸较远处激发的涌浪呈环形或山包形[10]。②动力差异导致形成的涌浪特征有差异。不同的岩体运动方式，激发的涌浪水体运动特征不一样。例如滑移式入水造成的水体涌浪以水平方向运动为主，而坠落垂直入水造成的水体涌浪以垂直方向的运动为主[1]。③水波类型的差异。不同岩体不同的冲击方式和不同的水深条件下，形成的涌浪波类型不同，如孤立波、椭圆余弦波以及非线性过渡波等[11]。而不同波型在传播过程中的衰减亦有差异。因此，在涌浪的产生阶段，应先根据危岩体的形态特征判断其失稳模式，从而选择公式计算危岩体的最大首浪高度。由于对滑坡涌浪的研究成果已经非常深入且丰富，因此本文对于滑移式危岩体不再进行讨论。
坠落式危岩体后部发育有近乎垂直的结构面，下部临空条件良好[14]，一旦危岩体失稳，在下落过程中基本不会受阻挡，危岩体主要发生垂直运动。在坠落式危岩冲击水体时，首先在水面激起“弹坑”，而后水从各个方向向弹坑处汇聚，当速度较大时可能在弹坑位置形成射流，水体被挤压向四周扩散，此时涌浪波以同心圆的形式向外传播，如图2所示。
对于坠落式危岩，美国土木工程协会推荐方法[15]和潘家铮法[16]均通过理论推导，采用垂直下落的箱型模型进行假设，在假设问题的几何参数和动力学参数已知的情况下，得到了理论解，同时进行了试验验证了理论公式的准确性，两者采取的模型与坠落式危岩体的运动模式一致，可以选取这两者方法推算的公式对坠落式危岩体涌浪的最大波幅进行评估。
压溃式危岩体的主要特点为存在破碎的基座。当基座在上方岩体的重力作用下被压裂，被裂隙切割的破碎危岩体随即失稳破坏，发生解体下沉运动，随后发生碎屑流运动[17]。水体受到碎屑体的推挤而向四周涌开，碎屑体动能持续传递给水体，在周围形成涌浪[18]，如图3所示。
对于压溃式危岩体，其运动模式与颗粒柱体崩塌下坠−滑动的复合运动的破坏模式较为类似。可采用张全等[17]和Huang[19]基于物理试验推导得出的计算公式，对压溃式危岩的入水速度和涌浪最大波幅进行估算。
倾倒式危岩体主控结构面倾角变化较大，且主控结构面下端部潜存于陡崖或陡坡岩体内，在重力作用下通常围绕主控结构面的下端部或下端部与临空面的交点旋转倾倒破坏。在旋转过程中主要受重力作用，在倾倒过程中将持续推挤水体向前运动，水体形成渐进的强烈飞溅浪花，随后形成涌浪波向前传播，如图4所示。
目前对于倾倒式危岩体所适用的涌浪公式计算方法研究较为罕见，Noda[15]的研究表明，对于一定范围内的一维或二维问题，当岩石块的体积相对于水体通道内水体积较小，并且满足其他理想化假设时，可以借助已知速度的箱型模型来预测波浪的生成情况。因此，可借鉴Maciel等[20]的简化方式，通过截取与其倾倒方向相垂直的剖面，将问题简化为三角形块体下落，并采用美国土木工程协会推荐的Noda理论推导的计算方法对于该危岩体的涌浪最大波幅进行近似计算。表2列举了不同失稳模式下对应的最大首浪计算公式及适用条件。
表 2 最大首浪公式Table 2. Part of propagation wave formula公式 来源 适用模式 适用条件 $ V = \dfrac{v}{{\sqrt {g{h_w}} }} $ $ \dfrac{{{H_s}}}{{{h_w}}} $ $ \dfrac{{{H_{\max }}}}{{{H_s}}} $ Noda[15]
（1970）坠落式
倾倒式h落＞hw $ \dfrac{{{H_{\max }}}}{{{h_w}}} = 1.17\dfrac{v}{{\sqrt {g{h_w}} }} $ 潘家铮[16]
（1980）坠落式 0.5<Fr数<2 $ \dfrac{{{H_{\max }}}}{{{h_w}}} = 0.159{a^{ - 0.641}}{\left( {\dfrac{v}{{\sqrt {g{h_w}} }}} \right)^{1.6}}{\left( {\dfrac{{{h_1}}}{{{h_w}}}} \right)^{0.839}}{\left( {\dfrac{d}{{{h_w}}}} \right)^{0.23}}{\left( {\dfrac{{{h_1}}}{{{h_2}}}} \right)^{0.294}} $
$ v = 0.636{a^{0.172}}{\left( {g{h_1}} \right)^{0.5}}{\left( {\dfrac{{{h_1}}}{{{h_w}}}} \right)^{ - 0.007}}{\left( {\dfrac{{{h_2}}}{{{h_w}}}} \right)^{0.156}}{\left( {\dfrac{{{h_1}}}{{{h_2}}}} \right)^{ - 0.003}} $张全等[17]
（2021）压溃式 2<a<16
30 m<h2<120 m式中，v为入水最大速度，hw为水深，Hs为滑坡体厚度，由V与$ \dfrac{{{H_s}}}{{{h_w}}} $确认波浪特性，Hmax为最大首浪高度，由V与$ \dfrac{{{H_{\max }}}}{{{H_s}}} $确认涌浪高度，h落为坠落高度；a为柱体高宽比，h1为危岩体高度，h2基座高度，d为危岩体宽度；Fr为相对Froude数，可根据$ {F_r} = \dfrac{v}{{\sqrt {g{h_w}} }} $进行计算。 1.2 涌浪传播阶段
涌浪在传播过程中将对航道产生严重威胁，可通过传播过程中的涌浪高度对航道危险性进行评价。根据王家山滑坡的相似模型试验[21]，涌浪的传播阶段可更细致地分为环状传播阶段和平行传播阶段。在这两个传播阶段中，涌浪波幅的衰减规律也有所不同，因此对涌浪的传播阶段的计算应采用两个适用范围相重叠的计算公式展开。如环状传播阶段可采用Huber and Hager[13]、Panizzo et al. [22]、Heller et al. [23]、Mohammed and Fritz [24]、Heller and Spinneken [25]物理试验推导的波幅计算公式；平行传播阶段可采用水科院经验公式[26]、潘家铮理论推导公式[16]、Huber and Hager物理试验推导公式[13]及殷坤龙和汪洋物理试验推导公式[27]等，其中Huber and Hager推导的两传播阶段计算公式、殷坤龙和汪洋等推导的两传播阶段计算公式可组合使用；根据Froude数范围选取公式组合进行计算，以表3所示的公式组合为例，对于0.54<Fr<0.9的崩塌体，可采用Heller and Spinneken [25]和殷坤龙和汪洋的公式，对环状传播和平行传播的涌浪进行计算。
表 3 部分传播浪公式Table 3. Part of propagation wave formula公式 来源 范围 适用条件 $ \dfrac{{{H_P}}}{{{h_w}}} = 1.67\sin (\alpha ){\cos ^2}\left(\dfrac{2}{3}\gamma \right){\left(\dfrac{\rho }{{{\rho _w}}}\right)^{\tfrac{1}{4}}}{\left[ {\dfrac{{{V_s}}}{{(d{h_w}^2)}}} \right]^{\tfrac{1}{2}}}{\left(\dfrac{r}{{{h_w}}}\right)^{ - \tfrac{2}{3}}} $ Huber and Hager[13](1997) 近场环状传播 高倾角入水 $ \dfrac{{{H_P}}}{{{h_w}}} = \dfrac{3}{4}{\left(\dfrac{v}{{\sqrt {g{h_w}} }}{\left(\dfrac{s}{{{h_w}}}\right)^{0.5}}{\left(\dfrac{{\rho {V_s}}}{{{\rho _w}d{h_w}^2}}\right)^{0.25}}{\left(\cos (\dfrac{6}{7}\alpha )\right)^{0.5}}{\left(\dfrac{x}{{{h_w}}}\right)^{ - \tfrac{1}{3}}}\right)^{\tfrac{4}{5}}} $ 沿程平行传播 $ \dfrac{{{H_P}}}{{{h_w}}} = 1.47\dfrac{{{H_{\max }}}}{{{h_w}}}{\left(\dfrac{x}{{{h_w}}}\right)^{ - 0.5}} $ 殷坤龙和汪洋[27](2008) 近场环状传播 0.063<Fr数<0.9
5°<α<45°$ \dfrac{{{H_P}}}{{{h_w}}} = \dfrac{{{H_{\max }}}}{{{h_w}}}{e^{ - 0.4\left(\tfrac{x}{{{h_w}}}\right)0.35}} $ 沿程平行传播 $ \dfrac{{{H_P}}}{{{h_w}}} = 2.75F_r^{0.67}S{M^{0.6}}{\left(\dfrac{r}{{{h_w}}}\right)^{ - 1}}{f_\gamma } $
$ {f_\gamma } = {\cos ^{2\left(1 + e^{ - 0.2\left(\tfrac{r}{{{h_w}}}\right)}\right)}}\left(\dfrac{2}{3}\lambda \right) $Heller and Spinneken[25]
(2015)近场环状传播 0.54<Fr数<2.47
α=45°式中，HP为传播浪高度，$\varTheta $为入水角度，©为滑动方向与径向夹角，$\lambda $为危岩体密度，$\lambda_w $为水密度，Vs为体积，r为径向距离，x为距危岩体距离；S为滑坡相对厚度，可根据$ S = \dfrac{{{H_s}}}{{{h_w}}} $进行计算，M为危岩相对质量，可根据$ M = \dfrac{m}{{{\rho _s}wh_w^2}} $进行计算，m为危岩体质量。 环状传播区大多可简化为以崩滑体入水点为半径、河道宽度为半径的半圆形区域，在涌浪的整个影响区域中占比较小，因此在计算危岩体涌浪时，也可仅采用平行传播浪的计算公式，再由传播浪的数据参与涌浪爬高的计算。
1.3 涌浪爬高阶段
通过涌浪传播过程中的高度可对航道内的危险情况进行预测，而对沿岸的城集镇和居民点而言，则需对涌浪爬高进行评估。涌浪爬坡高度受微地貌、首浪高度等等因素的影响，其规律很难总结。查阅文献，Synolakis [28]、Müller [29]及殷坤龙和汪洋[27]物理试验推导提出的公式可进行涌浪爬高计算。根据涌浪的传播规律，针对环状传播区和平行传播区，殷坤龙和汪洋[27]将涌浪的爬高也分为了正对岸爬高和沿程爬高，两者公式可结合使用，如表4所示。
表 4 部分爬高浪公式Table 4. Part of climbing wave formula公式 来源 适用范围 适用条件 $ \dfrac{{{H_R}}}{{{h_w}}} = 2.831{(\cot \beta )^{\tfrac{1}{2}}}{\left(\dfrac{{{H_P}}}{{{h_w}}}\right)^{\tfrac{5}{4}}} $ Synolakis [28]
(1987)正对岸爬高 坡比1:19.85 $ \dfrac{{{H_R}}}{{{h_w}}} = 1.25{\left(\dfrac{{90}}{\beta }\right)^{0.2}}{\left(\dfrac{{{H_P}}}{{{h_w}}}\right)^{1.25}}{\left(\dfrac{{{H_P}}}{L}\right)^{ - 0.15}} $ Müller [29]
(1995)正对岸爬高 坡比1:1、1:3或坡度90° $ \dfrac{{{H_R}}}{{{h_w}}} = 2.3\dfrac{{{H_P}}}{{{h_w}}}{\left(\dfrac{{90}}{\beta }\right)^{0.2}} $ 殷坤龙和汪洋[27](2008) 正对岸爬高 0.063<Fr数<0.9
5<α<45$ \dfrac{{{H_R}}}{{{h_w}}} = \left((2.3{\left(\dfrac{{90}}{\beta }\right)^{0.2}} - 1)\cos \delta + 1\right)\dfrac{{{H_P}}}{{{h_w}}} $ 沿程爬高 式中，HR为爬高，®为岸坡坡角，L为坝前波长；™为爬坡方位角，根据河道宽度B和计算点与滑坡的水平距离xs计算$ \cos \delta = \dfrac{B}{{\sqrt {{B^2} + x_s^2} }} $。 选用正对岸爬高公式计算危岩体入江位置正对岸的涌浪爬高进行计算，选用沿程爬高公式对涌浪影响范围内的爬高进行计算，同时可以沿河道平行传播的浪高分别计算该位置的正对岸爬高，以此Synolakis 和Müller 提出的公式可近似计算沿程的爬高值。
2. 考虑失稳模式的涌浪公式计算体系和计算引擎
根据涌浪全过程规律和现有的阶段式公式，可组合形成考虑失稳模式的涌浪公式计算体系，并以此开发考虑失稳模式的涌浪公式计算引擎，计算引擎的流程图如图5所示，根据失稳模式选取最大首浪高度计算公式计算首浪高度，根据危岩体实际参数范围选取传播浪高度计算公式计算传播浪高度，根据危岩体所处河道信息选取爬高计算公式计算涌浪爬坡高度，引擎页面如图6所示。
为了快速评价涌浪灾害，依照危岩体涌浪公式计算体系，本研究团队基于Java Spring技术体系开发了危岩体涌浪计算引擎。该引擎采用MVC分层架构和Spring多层框架，核心开发语言为Java，前台使用JavaScript，并通过cesium WebGL引擎生成二维/三维涌浪模型，引擎支持浏览器/服务器应用模式，数据存储采用MySQL关系型数据库。
该计算引擎利用三维地图技术进行涌浪危岩模型的计算与展示，并具备涌浪危害范围分析以及河道最大波高剖面线绘制的功能。引擎运用特定的算法和公式进行涌浪参数的计算，允许用户根据不同情况选择公式组合并输入参数，以进行动态参数计算。同时，引擎支持细粒度的数据处理，以提高数据处理的精细度和准确性，并具备快速响应的能力。
用户可以模拟不同的崩塌场景，并自定义涌浪计算的范围和水位条件。引擎能够生成包含计算结果的报告，报告中通过三维模型和图表等形式展示数据，以展示危岩崩塌可能引发的涌浪波高分布及其影响范围。
3. 实例计算
3.1 压溃式危岩体涌浪计算
数值模拟作为一种重要的科学研究方法，在物理规律清晰、模型合理的情况下，数值模拟的结果在大多数情况下能够较好地反映现实情况，具有一定的参考价值。基于考虑失稳模式的危岩体涌浪公式计算体系和计算引擎，以压溃式的典型危岩体进行涌浪计算和危险性分析，与数值模拟计算结果进行对比来验证公式计算体系的可行性，计算过程遵循图5所示流程。
三峡库区龙门寨危岩体位于中国长江支流大宁河水域小三峡中的龙门峡右岸[2]。岩体位于著名的大宁河小三峡旅游景区，距巫山县城2.9公里。危岩体上游处因贯穿裂缝发育而与母体山体分离，形成孤立的柱状结构，平均高差为190 m，平均危岩横宽和厚度约40 m，主崩方向为250°，危岩体方量约为30.4万m3，危岩体下有厚约12 m 的泥灰岩劣化基座, 泥灰岩平台下部为碎裂结构泥质灰岩陡崖，垂向劈理发育，受长期水位变化的影响,其岩体强度逐年降低，危岩体可能发生压溃式破坏，地质剖面图如图7。龙门峡处河谷宽约200 m， 河床高程约85 m，涌浪计算参数如表5。
![]()
图 7 龙门寨危岩体剖面图（郑嘉豪等，2020）Figure 7. Profile of Longmenzhai perilous rock mass (Zheng Jiahao et al., 2020)根据危岩体的失稳模式和龙门寨危岩体参数范围，选用张全等[17]的计算公式，对其入水速度和涌浪最大波幅进行估算。由于是高倾角入水且存在支流，两岸并非完全平行，因此选择水科院[26]提出的公式进行传播浪高度计算，危岩体坡度接近90°，最后选用Müller [29]推导的公式计算典型位置的涌浪爬高，涌浪的计算范围为危岩体沿河道上下2 km。
表 5 危岩体计算参数Table 5. Calculation parameters of dangerous rock mass符号 g hw v h1 h2 d a Vs 参数 重力加速度（m/s2） 水深（m） 入水速度（m/s） 高度（m） 基座高度（m） 宽度（m） 柱体高宽比 体积（万m³） 145m 9.8 65 27.1 190 12 40 4.75 30.4 175m 95 25.6 以175 m水位为例，图8展示了考虑失稳模式的涌浪公式计算体系和郑嘉豪等[2]对箭穿洞危岩体涌浪数值模拟的涌浪计算结果，图8（a）展示了涌浪沿河道传播中的波幅情况；图8（b）展示了按照《滑坡涌浪危险性评估》规范[30]计算的涌浪危险程度分区及沿岸居民点位置的爬高情况，数值模拟得最大涌浪高度为11.6 m，公式体系计算得最大涌浪高度为13.9 m，对比计算结果，由张全等推导的公式计算出的结果较数值模拟偏高，但整体误差不超过20%；由河道波幅传播图可知，水科院公式计算出传播浪较数值模拟结果偏低，但整体衰减规律较为吻合，涌浪的衰减由快到慢，500 m以内涌浪急速衰减，500 m后涌浪的衰减逐渐变缓。综合来看，公式体系计算结果显示危岩体的造浪特征和传播规律与数值模拟结果较为吻合，因此这套公式体系可用来粗略对龙门寨危岩体的涌浪进行计算，同时也验证了涌浪公式法计算体系的可行性。将验证过的公式计算体系应用到典型坠落式危岩体和倾倒式危岩体失稳破坏产生的涌浪计算与分析之中。
![]()
图 8 压溃式涌浪公式计算体系计算结果Figure 8. The calculation results of the collapse surge formula calculation system3.2 坠落式危岩体涌浪计算
2020年7月10日，渔峡口危岩失稳坠落，渔峡口危岩位于清江右岸，距渔峡口镇约1 km处。距离上游水布垭镇直线距离16 km，下游长阳县城68 km。危岩分布高程239 m～272 m，总体为陡崖，坡面近直立，危岩体所在岸坡相对高差100 m，危岩体被裂隙切割，下部并无阻碍，结构面贯通后危岩体随之坠落，为典型的坠落式危岩体，渔峡口码头监控记录下了渔峡口危岩体的崩塌—造浪和传播—爬高过程，如图9所示，图9（a）为危岩体失稳入水时刻，而后激起高约26 m的水舌（图9（b）），形成约70 m的最大射流（图9（c）），并产生涌浪沿着河道进行传播（图9（d）），最大涌浪波幅约为10.4 m。
![]()
图 9 渔峡口危岩体涌浪过程Figure 9. Surge process of Yuxiakou perilous rock mass（a） t = 0 s；（b） t = 4 s；（c） t =8 s；（d） t =14 s由于渔峡口危岩的失稳模式为坠落式，坠落高度大于水深，近乎垂直入水，且两岸并非完全平行，应优先选用美国土木工程学会推荐首浪高度公式[15]和水科院传播公式[26]的公式体系，根据影像提取出表6所示的计算参数进行计算。
表 6 危岩体计算参数Table 6. Calculation parameters of dangerous rock mass符号 g hw v HS Vs 参数 重力加速度
（m/s2）水深
（m）入水速度
（m/s）平均厚度
（m）体积
（万m³）渔峡口危岩体 9.8 63 30 3 0.01 如图10所示，计算得出渔峡口危岩体入水的最大涌浪高度为2.61 m，在传播200 m后降低至0.16 m，进入到低风险范围，与视频影像中传播至渔船位置，涌浪波对渔船未有较大影响相接近。
![]()
图 10 坠落式涌浪公式计算体系计算结果Figure 10. The calculation results of the falling surge formula calculation system3.3 倾倒式危岩体涌浪计算
2022年1月8日，一块状危岩体从Capitólio旅游区的悬崖上失稳折断，倾倒在湖面上漂泊的游船船员身上，同时形成涌浪，造成10人死亡，30人受伤。
由Maciel等[20]，将危岩体进行模型假设为高为30 m，宽为4 m，厚为4 m的箱型模型。同时假设倾倒危岩体在接触水面时为已知速度9.42±0.94 m/s的箱型模型垂直下落，如图11。Maciel等随即应用美国土木工程协会推荐法[15]对此次危岩造成的最大首浪波幅进行了计算，最终得到最大波幅为3.75±0.06 m。
![]()
图 11 Capitólio危岩体涌浪过程（Maciel，2023）：（a）危岩体原始尺寸；（b）发生倾倒及概化尺寸；（c）形成飞溅及入水速度Figure 11. Surge process of Capitólio perilous rock mass (Maciel, 2023)由于此处水深为4m，可忽略水深对于涌浪的影响，则选用水科院传播公式[26]的公式组合，如图11所示的计算参数进行计算。
如图12所示，巴西Capitólio危岩体入水的最大涌浪高度为3.75 m，传播200 m后降低至0.26 m，进入到低风险范围。
![]()
图 12 倾倒式涌浪公式计算体系计算结果Figure 12. Calculation results of toppling surge formula calculation system4. 结论与建议
4.1 结论
（1）针对危岩崩滑体涌浪现象的复杂性和多样性，本文首先对现有的涌浪计算公式进行了系统性的筛选和分析。基于危岩崩滑体涌浪的三个阶段以及失稳模式不同，成功地建立了考虑不同失稳模式的危岩体涌浪公式计算体系。
（2）为了快速评价涌浪灾害，依照崩滑体涌浪公式计算体系，结合计算机技术开发崩滑体涌浪计算引擎。基于三维地图实现涌浪崩塌模型计算、崩塌涌浪展示、涌浪危害范围、河道最大波高剖面线等功能，可大大提高对危岩体崩塌产生的涌浪进行分析效率与准确度，有利于库区涌浪防灾减灾。
（3）运用涌浪公式法计算体系和涌浪计算引擎对于典型压溃式危岩体龙门寨危岩体进行涌浪分析，发现175 m时最大首浪高度为13.9 m，传播至2 km处传播浪高度为1.75 m，码头处爬高值为2.91 m，而数值模拟的计算结果分别为11.9 m、2 m、2.1 m，误差在20%以内，验证了涌浪公式法计算体系和涌浪计算引擎的可行性，并运用相同方法对于典型坠落式危岩体和倾倒式危岩体进行了涌浪分析，两者危岩体涌浪传播200 m后都进入低风险区域，体现了该计算体系在不同失稳模式下的应用情况。
4.2 建议
（1）涌浪公式法计算体系能够快速计算涌浪，但在实际应用中仍需根据具体情况进行修正和调整。因为涌浪的产生和传播受到多种因素的影响，如地形、水深、速度等[31]，这些因素的变化都可能对涌浪的计算结果产生影响。因此，在利用该计算体系进行涌浪预测时，需要充分考虑实际情况选取公式，以提高预测结果的准确性和可靠性。
（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 宋庄地裂缝浅表部地层剖面图（据赵龙，2018修改）
Figure 3. Stratigraphic section of the shallow strata of Songzhuang ground fissure (modified from Zhao Long, 2018)
![]()
图 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 双埠头村地裂缝上盘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 沟渠村地裂缝上盘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
![]()
图 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 双埠头村地裂缝上盘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 沟渠村地裂缝上盘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
![]()
图 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
![]()
图 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 选取地震仪的主要技术参数
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
-
[1] 何红前. 渭河盆地地裂缝成因机理研究[D]. 西安:长安大学,2011. [HE Hongqian. Study on genetic mechanism of ground fissures in Weihe Basin[D]. Xi’an:Changan University,2011. (in Chinese with English abstract)] HE Hongqian. Study on genetic mechanism of ground fissures in Weihe Basin［D］. Xi’an: Changan University, 2011. （in Chinese with English abstract）
[2] 王景明,地裂缝及其灾害的理论与应用. 地裂缝及其灾害的理论与应用[M]. 西安:陕西科学技术出版社,2000. [WANG Jingming. Theory of ground fissures hazards and its application[M]. Xi’an:Shaanxi Science & Technology Press,2000. (in Chinese)] WANG Jingming. Theory of ground fissures hazards and its application［M］. Xi’an: Shaanxi Science & Technology Press, 2000. （in Chinese）
[3] 赵忠海. 北京地区地裂缝灾害的分布特征及其成因探讨[J]. 地质灾害与环境保护,2006,17(3):75 − 78. [ZHAO Zhonghai. Discussion on the distribution characteristics and genetic type of the land crack in Beijing[J]. Journal of Geological Hazards and Environment Preservation,2006,17(3):75 − 78. (in Chinese with English abstract)] DOI: 10.3969/j.issn.1006-4362.2006.03.018 ZHAO Zhonghai. Discussion on the distribution characteristics and genetic type of the land crack in Beijing［J］. Journal of Geological Hazards and Environment Preservation, 2006, 17（3）: 75 − 78. （in Chinese with English abstract） DOI: 10.3969/j.issn.1006-4362.2006.03.018
[4] 耿大玉,李忠生. 中美两国的地裂缝灾害[J]. 地震学报,2000,22(4):433 − 441. [GENG Dayu,LI Zhongsheng. Ground fissure hazards in the United States and China[J]. Acta Seismologica Sinica,2000,22(4):433 − 441. (in Chinese with English abstract)] GENG Dayu, LI Zhongsheng. Ground fissure hazards in the United States and China［J］. Acta Seismologica Sinica, 2000, 22（4）: 433 − 441. （in Chinese with English abstract）
[5] Michael C. Carpenter. Earth fissures and subsidence complicate development of desert water resources [J]. U. S. Geological Survey,Tucson,Arizona,1999,pp. 65 − 78.
[6] AYRES J E. Two Clovis fluteld points from southern Arizona[J]. Kiva,1970,35(3):121 − 124. DOI: 10.1080/00231940.1970.11760562
[7] Heidke,James M. The Earliest Tucson Basin Pottery[J]. Archaeology in Tucson,1997,11(3):9 − 10.
[8] Heilen,Michael P. An Archaeological Theory of Landscapes [D]. Ph. D. dissertation,Department of Anthropology,University of Arizona,Tucson,2005.
[9] HUCKELL B B. The paleo-Indian and archaic occupation of the Tucson Basin:an overview[J]. Kiva,1984,49(3/4):133 − 145.
[10] KREITLER C W. Fault control of subsidence,Houston,texasa[J]. Groundwater,1977,15(3):203 − 214. DOI: 10.1111/j.1745-6584.1977.tb03165.x
[11] 赵龙,北京市平原区高精度(1∶5万)地面沉降和地裂缝调查[R],北京市,北京市地质工程勘察院,2018-12-01. [Zhao L,A high-precision (1∶50,000) survey of ground settlement and ground fissures in the plain area of Beijing. [R],Beijing,Beijing Geological Engineering Survey Institute,2018-12-01 (in Chinese)] Zhao L, A high-precision （1∶50, 000） survey of ground settlement and ground fissures in the plain area of Beijing. ［R］, Beijing, Beijing Geological Engineering Survey Institute, 2018-12-01 （in Chinese）
[12] 贾三满,郭萌. 从高丽营探槽分析黄庄—高丽营断裂与地裂缝的关系[J]. 城市地质,2007,2(4):24 − 28. [JIA Sanman,GUO Meng. The relation between huangzhung-gaoliying fault and by gaoliying trench and earth fissure[J]. City Geology,2007,2(4):24 − 28. (in Chinese with English abstract)] JIA Sanman, GUO Meng. The relation between huangzhung-gaoliying fault and by gaoliying trench and earth fissure［J］. City Geology, 2007, 2（4）: 24 − 28. （in Chinese with English abstract）
[13] 彭建兵,卢全中,邓亚虹,等. 国网智能电网研究院工程场地地裂缝勘察评价报告[R]. 西安:长安大学工程设计研究院,2010. [Peng J,Lu Q,Deng Y,et al,Report on the Investigation and Evaluation of Ground Fissures at the State Grid Smart Grid Research Institute Project Site [R],Xi’an,Chang'an University Engineering Design Institute,2010 (in Chinese)] Peng J, Lu Q, Deng Y, et al, Report on the Investigation and Evaluation of Ground Fissures at the State Grid Smart Grid Research Institute Project Site ［R］, Xi’an, Chang'an University Engineering Design Institute, 2010 （in Chinese）
[14] 刘德成,靳小平,周自梁. 灰色关联度分析法在北京通州区地裂缝灾害危险性评价中的应用[J]. 中国地质灾害与防治学报,2010,21(3):73 − 76. [LIU Decheng,JIN Xiaoping,ZHOU Ziliang. Application of grey degree of association evaluation method in ground fissures in Tongzhou District of Beijing[J]. The Chinese Journal of Geological Hazard and Control,2010,21(3):73 − 76. (in Chinese with English abstract)] LIU Decheng, JIN Xiaoping, ZHOU Ziliang. Application of grey degree of association evaluation method in ground fissures in Tongzhou District of Beijing［J］. The Chinese Journal of Geological Hazard and Control, 2010, 21（3）: 73 − 76. （in Chinese with English abstract）
[15] 姜媛,王荣,田芳,等. 北京地区地面沉降与地裂缝关系研究[J]. 城市地质,2014,9(4):6 − 10. [JIANG Yuan,WANG Rong,TIAN Fang,et al. Study on the relationship between land subsidence and ground fissures in Beijing[J]. Urban Geology,2014,9(4):6 − 10. (in Chinese with English abstract)] JIANG Yuan, WANG Rong, TIAN Fang, et al. Study on the relationship between land subsidence and ground fissures in Beijing［J］. Urban Geology, 2014, 9（4）: 6 − 10. （in Chinese with English abstract）
[16] 卢全中,彭建兵,邓亚虹,等. 北京北七家—高丽营地裂缝破坏特征及影响带宽度[J]. 工程勘察,2014,42(6):5 − 11. [LU Quanzhong,PENG Jianbing,DENG Yahong,et al. Failure characteristics and influence width of Beiqijia-Gaoliying ground fissure in Beijing[J]. Geotechnical Investigation & Surveying,2014,42(6):5 − 11. (in Chinese with English abstract)] LU Quanzhong, PENG Jianbing, DENG Yahong, et al. Failure characteristics and influence width of Beiqijia-Gaoliying ground fissure in Beijing［J］. Geotechnical Investigation & Surveying, 2014, 42（6）: 5 − 11. （in Chinese with English abstract）
[17] 刘方翠,祁生文,彭建兵,等. 北京市地裂缝分布与发育规律[J]. 工程地质学报,2016,24(6):1269 − 1277. [LIU Fangcui,QI Shengwen,PENG Jianbing,et al. Characters of the ground fissures developing in Beijing[J]. Journal of Engineering Geology,2016,24(6):1269 − 1277. (in Chinese with English abstract)] LIU Fangcui, QI Shengwen, PENG Jianbing, et al. Characters of the ground fissures developing in Beijing［J］. Journal of Engineering Geology, 2016, 24（6）: 1269 − 1277. （in Chinese with English abstract）
[18] 周永恒,杨肖肖,丰成君,等. 北京平原区黄庄—高丽营断裂(房山—涞水段)第四纪活动特征的浅层综合探测证据[J]. 地球学报,2021,42(5):677 − 689. [ZHOU Yongheng,YANG Xiaoxiao,FENG Chengjun,et al. Evidence of shallow synthetic exploration of quaternary activity characteristics along Fangshan–Laishui section of Huangzhuang–gaoliying fault in Beijing Plain[J]. Acta Geoscientica Sinica,2021,42(5):677 − 689. (in Chinese with English abstract)] DOI: 10.3975/cagsb.2020.012502 ZHOU Yongheng, YANG Xiaoxiao, FENG Chengjun, et al. Evidence of shallow synthetic exploration of quaternary activity characteristics along Fangshan–Laishui section of Huangzhuang–gaoliying fault in Beijing Plain［J］. Acta Geoscientica Sinica, 2021, 42（5）: 677 − 689. （in Chinese with English abstract） DOI: 10.3975/cagsb.2020.012502
[19] 关金环,高明亮,宫辉力. 首都国际机场区域差异性沉降原因探讨[J]. 测绘科学,2021,46(9):67 − 75. [GUAN Jinhuan,GAO Mingliang,GONG Huili. Discussion on the causes of regional differential settlement of Beijing-capital international airport[J]. Science of Surveying and Mapping,2021,46(9):67 − 75. (in Chinese with English abstract)] GUAN Jinhuan, GAO Mingliang, GONG Huili. Discussion on the causes of regional differential settlement of Beijing-capital international airport［J］. Science of Surveying and Mapping, 2021, 46（9）: 67 − 75. （in Chinese with English abstract）
[20] 孟振江,彭建兵,李超,等. 耦合型地裂缝活动特征与成因机制模拟研究——以北京宋庄地裂缝为例[J]. 水文地质工程地质,2023,50(3):138 − 148. [MENG Zhenjiang,PENG Jianbing,LI Chao,et al. A simulation study of the activity characteristics and genetic mechanism of coupled ground fissures:exemplified by the songzhuang ground fissure in Beijing[J]. Hydrogeology & Engineering Geology,2023,50(3):138 − 148. (in Chinese with English abstract)] MENG Zhenjiang, PENG Jianbing, LI Chao, et al. A simulation study of the activity characteristics and genetic mechanism of coupled ground fissures: exemplified by the songzhuang ground fissure in Beijing［J］. Hydrogeology & Engineering Geology, 2023, 50（3）: 138 − 148. （in Chinese with English abstract）
[21] 任雅哲,丰成君,戚帮申,等. 北京顺义断裂活动对首都机场地裂缝影响定量研究[J]. 地质力学学报,2023,29(5):685 − 703. [REN Yazhe,FENG Chengjun,QI Bangshen,et al. Quantitative research of the impact of Shunyi fault activity on the ground fissures in the Beijing Capital International Airport,China[J]. Journal of Geomechanics,2023,29(5):685 − 703. (in Chinese with English abstract)] DOI: 10.12090/j.issn.1006-6616.2023063 REN Yazhe, FENG Chengjun, QI Bangshen, et al. Quantitative research of the impact of Shunyi fault activity on the ground fissures in the Beijing Capital International Airport, China［J］. Journal of Geomechanics, 2023, 29（5）: 685 − 703. （in Chinese with English abstract） DOI: 10.12090/j.issn.1006-6616.2023063
[22] 林通. 地层结构和场地特征对地脉动测试的影响[D]. 上海:同济大学,2008. [LIN Tong. Influence of stratum structure and site characteristics on ground pulsation test[D]. Shanghai:Tongji University,2008. (in Chinese with English abstract)] LIN Tong. Influence of stratum structure and site characteristics on ground pulsation test［D］. Shanghai: Tongji University, 2008. （in Chinese with English abstract）
[23] 吴志坚,王兰民,孙军杰,等. 青藏高原典型多年冻土场地地面脉动特征[J]. 岩石力学与工程学报,2008,27(11):2316 − 2323. [WU Zhijian,WANG Lanmin,SUN Junjie,et al. Characteristics of microseism at typical permafrost sites in Qinghai—Tibet Plateau[J]. Chinese Journal of Rock Mechanics and Engineering,2008,27(11):2316 − 2323. (in Chinese with English abstract)] WU Zhijian, WANG Lanmin, SUN Junjie, et al. Characteristics of microseism at typical permafrost sites in Qinghai—Tibet Plateau［J］. Chinese Journal of Rock Mechanics and Engineering, 2008, 27（11）: 2316 − 2323. （in Chinese with English abstract）
[24] 薛捷,王郑睿,王蕾. 用地脉动测试的方法对地裂缝场地的评价[J]. 国土资源导刊,2014,11(10):70 − 73. [XUE Jie,WANG Zhengrui,WANG Lei. Evaluation of ground fissure site by the method of ground pulsation test[J]. Land & Resources Herald,2014,11(10):70 − 73. (in Chinese)] XUE Jie, WANG Zhengrui, WANG Lei. Evaluation of ground fissure site by the method of ground pulsation test［J］. Land & Resources Herald, 2014, 11（10）: 70 − 73. （in Chinese）
[25] 慕焕东. 地裂缝场地地震放大效应研究[D]. 西安:长安大学,2014. [MU Huandong. Study on seismic amplification effect of ground fissure site[D]. Xi’an:Changan University,2014. (in Chinese with English abstract)] MU Huandong. Study on seismic amplification effect of ground fissure site［D］. Xi’an: Changan University, 2014. （in Chinese with English abstract）
[26] 张磊刚,邓亚虹,薛捷,等. 地裂缝场地地脉动响应特征初探——以西安F6地裂缝为例[J]. 工程地质学报,2017,25(3):784 − 793. [ZHANG Leigang,DENG Yahong,XUE Jie,et al. Preliminary research on dynamic response characteristics of microtremor at ground fissure sites:Illustrated with F6 ground fissure in Xi’an[J]. Journal of Engineering Geology,2017,25(3):784 − 793. (in Chinese with English abstract)] ZHANG Leigang, DENG Yahong, XUE Jie, et al. Preliminary research on dynamic response characteristics of microtremor at ground fissure sites: Illustrated with F6 ground fissure in Xi’an［J］. Journal of Engineering Geology, 2017, 25（3）: 784 − 793. （in Chinese with English abstract）
[27] 王晗,邓亚虹,慕焕东,等. 西安典型地裂缝场地地脉动测试及地震响应特征分析[J]. 中国地质灾害与防治学报,2022,33(4):55 − 64. [WANG Han,DENG Yahong,MU Huandong,et al. Ground pulsation tests and analysis on seismic response of typical ground fissure sites in Xi’an[J]. The Chinese Journal of Geological Hazard and Control,2022,33(4):55 − 64. (in Chinese with English abstract)] WANG Han, DENG Yahong, MU Huandong, et al. Ground pulsation tests and analysis on seismic response of typical ground fissure sites in Xi’an［J］. The Chinese Journal of Geological Hazard and Control, 2022, 33（4）: 55 − 64. （in Chinese with English abstract）
[28] 崔思颖,邓亚虹,曹歌,等. 河北平原典型地裂缝场地动力响应特征研究[J]. 煤田地质与勘探,2022,50(6):81 − 93. [CUI Siying,DENG Yahong,CAO Ge,et al. Research on dynamic response characteristics of typical ground fissure sites in Hebei Plain[J]. Coal Geology & Exploration,2022,50(6):81 − 93. (in Chinese with English abstract)] CUI Siying, DENG Yahong, CAO Ge, et al. Research on dynamic response characteristics of typical ground fissure sites in Hebei Plain［J］. Coal Geology & Exploration, 2022, 50（6）: 81 − 93. （in Chinese with English abstract）
[29] CHANG J,DENG Y H,XUAN Y,et al. The dynamic response of sites with earth fissures as revealed by microtremor analysis:A case study in the Linfen Basin,China[J]. Soil Dynamics and Earthquake Engineering,2020,132:106076. DOI: 10.1016/j.soildyn.2020.106076
[30] CHANG J,DENG Y H,MU H D,et al. A microtremor study to reveal the dynamic response of earth fissure site:the case study in Fenwei Basins,China[J]. Environmental Earth Sciences,2022,81(3):80. DOI: 10.1007/s12665-022-10217-y
[31] XUAN Y,DENG Y H,HE J,et al. Microtremor-based analysis of the dynamic response characteristics of earth-fissured sites in the Datong basin,China[J]. Earthquake Engineering and Engineering Vibration,2021,20(3):567 − 582. DOI: 10.1007/s11803-021-2039-2
[32] XUAN Y,DENG Y H,CHANG J,et al. Dynamic characteristics of fault-controlled earth fissure sites in Jingyang County,Weihe Basin,China[J]. Acta Geophysica,2022,70(4):1687 − 1700. DOI: 10.1007/s11600-022-00833-4
[33] 杨涛,宫辉力,赵文吉,等. 北京顺义区地裂缝分布特征及成因分析[J]. 自然灾害学报,2010,19(6):100 − 106. [YANG Tao,GONG Huili,ZHAO Wenji,et al. Distribution characteristics and cause analysis of ground fissures in Shunyi District of Beijing[J]. Journal of Natural Disasters,2010,19(6):100 − 106. (in Chinese with English abstract)] YANG Tao, GONG Huili, ZHAO Wenji, et al. Distribution characteristics and cause analysis of ground fissures in Shunyi District of Beijing［J］. Journal of Natural Disasters, 2010, 19（6）: 100 − 106. （in Chinese with English abstract）
[34] 赵龙,罗勇,李玉梅,等. 北京平原区地裂缝受灾体形态特征及影响因素[J]. 水文地质工程地质,2019,46(6):156 − 164. [ZHAO Long,LUO Yong,LI Yumei,et al. Characteristics of disaster-affected bodies and influence factors for earth fissure in Beijing Plain[J]. Hydrogeology & Engineering Geology,2019,46(6):156 − 164. (in Chinese with English abstract)] ZHAO Long, LUO Yong, LI Yumei, et al. Characteristics of disaster-affected bodies and influence factors for earth fissure in Beijing Plain［J］. Hydrogeology & Engineering Geology, 2019, 46（6）: 156 − 164. （in Chinese with English abstract）
[35] 赵龙,李玉梅,崔文君,等. 北京宋庄地裂缝灾害特征及影响因素分析[J]. 工程地质学报,2018,26(6):1600 − 1610. [ZHAO Long,LI Yumei,CUI Wenjun,et al. Disaster characteristics and influence factors for ground fissures at songzhuang village in Beijing[J]. Journal of Engineering Geology,2018,26(6):1600 − 1610. (in Chinese with English abstract)] ZHAO Long, LI Yumei, CUI Wenjun, et al. Disaster characteristics and influence factors for ground fissures at songzhuang village in Beijing［J］. Journal of Engineering Geology, 2018, 26（6）: 1600 − 1610. （in Chinese with English abstract）
[36] 赵晨曦,陈柘舟,刘建凯,等. 基于野外实地调查的北京宋庄地裂缝经济损失评估[J]. 上海国土资源,2023,44(2):113 − 118. [ZHAO Chenxi,CHEN Zhezhou,LIU Jiankai,et al. Economic loss assessment of ground fissures in Songzhuang,Beijing based on field investigation[J]. Shanghai Land & Resources,2023,44(2):113 − 118. (in Chinese with English abstract)] ZHAO Chenxi, CHEN Zhezhou, LIU Jiankai, et al. Economic loss assessment of ground fissures in Songzhuang, Beijing based on field investigation［J］. Shanghai Land & Resources, 2023, 44（2）: 113 − 118. （in Chinese with English abstract）
计量
- 文章访问数: 31
- HTML全文浏览量: 7
- PDF下载量: 18
邮件订阅
RSS