“情景-应对”型理论体系的发展及其在地震灾害应急管理中的应用探讨

陈波; 王芳; 肖本夫

doi:10.11899/zzfy20210401

“情景-应对”型理论体系的发展及其在地震灾害应急管理中的应用探讨

doi: 10.11899/zzfy20210401

1.
中国地震局地球物理研究所, 北京100081
2.
四川省地震局, 成都610041

基金项目: 国家重点研发计划项目 (2019YFC1509402)；中国地震局地球物理研究所基本科研业务费专项(DQJB19A0132)

详细信息

作者简介:
陈波，男，生于1987年。博士，副研究员。主要从事地震工程学与城市减灾技术研究工作。E-mail：chenbo@cea-igp.ac.cn

计量
- 文章访问数: 1092
- HTML全文浏览量: 334
- PDF下载量: 138
- 被引次数: 0
出版历程
- 收稿日期: 2021-10-18
- 刊出日期: 2021-12-31

The Development of “Scenario-response” Theoretical System and Its Application in Earthquake Disaster Emergency Management

1.
Institute of Geophysics, China Earthquake Administration, Beijing 100081, China
2.
Sichuan Earthquake Agency, Chengdu 610041, China

摘要

摘要: “情景-应对”型理论体系基于事件静、动态情景分析生成应急决策方案，已成为重大突发事件应急决策的基本范式。近年来，“情景-应对”研究理论的迅速发展极大地拓宽了其应用范围，逐渐被引入地震应急中，在一定程度上能够弥补传统“预测-应对”模式的不足，为震后高效行动及充分应对提供有效支撑。本文从情景识别、情景推演及应对处置方案生成等方面详细介绍了“情景-应对”理论体系发展中涉及的关键技术，阐述了其在地震灾害应急管理中的应用。通过对已有研究的回顾，分析讨论了该领域的最新进展及需进一步解决的问题。
- 情景-应对 /
- 情景识别 /
- 情景推演 /
- 地震灾害 /
- 应急管理
Abstract: The “scenario-response” theoretical system, which generates emergency decision-making schemes based on the static and dynamic scenario analysis of unexpected incidents, has become the basic paradigm of emergency decision-making for major emergencies. In recent years, the rapid development of the “scenario-response” theory has greatly broadened its scope of application, and has been gradually applicable to earthquake emergency. To some extent, it could make up for the shortcomings of the traditional “prediction-response” model, and provide effective support for efficient and adequate post-earthquake actions. This paper introduces the key technologies involved in the development of the “scenario-response” from the aspects of scenario recognition, scenario deduction, and emergency decision-making schemes generation in detail, and then elaborates its application in earthquake disaster emergency management. Through the review of previous studies, the new progress in the field and the existing problems are analyzed and discussed.
- Scenario-response /
- Scenario recognition /
- Scenario deduction /
- Earthquake disaster /
- Emergency management

HTML全文

. 引言

静乐井水位是山西乃至华北地区映震能力较强的地下流体观测井之一，井水位异常能客观反映区域应力场的变化，全球7级以上地震均能记录到水震波。1983年始测以来出现的几次高值异常与山西北部至晋冀蒙交界地区强震活动时间具有准同步性，在1989年大同-阳高M_s6.1级地震、1991年忻州M_s5.1级地震、1991年大同-阳高M_s5.8级地震、1998年张北M_s6.2级地震、1999年大同-阳高M_s5.6级等中强地震发生前有较好的异常显示。因此将静乐井水位高值异常作为晋冀蒙交界地区地震趋势预测的重要判据（张淑亮等，1997；车用太等，1999；王爱英等，1998）。

2017年8月静乐井水位又出现大于多年平均值的高值异常，到目前为止异常仍在持续。按照以往震例时间指示意义，目前已接近发震时间，未来晋冀蒙交界区震情形势如何发展？目前的异常是否真实反映了构造活动加剧所引起的区域应力场和深部地球物理场的改变？是否与中强震孕育过程有关？静乐井水位目前的异常特征不足以解析异常区的应力状态，从而也直接影响对晋冀蒙交界地区地震形势的认识和未来地震趋势的判定。

多年的观测表明，影响静乐井水位动态变化的主要干扰因素为短时强降雨与河流水体的荷载作用。无论是短时间的强降雨还是洪水期河流的荷载作用，均可在短时间内造成井水位的快速变化（张淑亮等，2005；车用太等，2004；鱼金子等，1994；张昭栋等，1989；张昭栋等，1990）。目前对这2种干扰因素的识别与排除基本上以定性分析为主，定量分析方法较少。特别是对于同时受2种因素影响的静乐井水位而言，在异常分析与判定工作中，仅采用与同期降雨量对比的方法分析井水位变化与降雨量的关系，而河流水位变化对井水位的影响分析工作至今尚未开展。为此，本文利用数值模拟与含水层垂向应力反演的方法，对降雨与河流荷载作用下的静乐井区应力进行定量计算，并从区域应力场的角度分析2种方法计算结果的差异，以期为地震趋势判定提供较可靠的预测依据。

1. 井区概况

静乐流体观测井位于吕梁山断块隆起区的静乐向斜，位于东碾河断裂带上，地处东碾河南岸，与东碾河相距约130m，其北部和南部均为山区（见图 1）。静乐井深362.92m，井区为河谷地区，主要发育有下古生界灰岩、白云岩喀斯特含水层及第四系冲洪积砂砾石含水层，2个含水层间有较强的水力联系。观测层岩性时代为中奥陶系灰岩裂隙溶洞水，其顶板埋深为46.3m，单位涌水量84.5L/（s·m），渗透系数为2.41—11.5m/d。观测层具有厚度大、透水性强的特点。含水层岩石溶洞发育，静乐井在230—290m处通过破碎带，地下水由北西经井区向东南方向排泄，井区处于强径流区，地下水补给主要来自北、南山区基岩裸露区的大气降水补给和东碾河河水沿断裂破碎带入渗补给。

图 1 井区地质构造

Figure 1. Well area tectonic map

下载: 全尺寸图片幻灯片

2. 强降雨与河流荷载作用下静乐井水位动态特征

多年的观测表明，大于15mm的集中降雨将引起静乐井水位的快速上升（张淑亮等，2005），可能由于降雨在地表形成荷载，使得含水层受力状态发生变化，引起水井水位上升。东碾河水位变化是影响静乐井水位变化的又一重要因素。当大面积洪水作用通过井区时，井水位迅速大幅上升，几乎没有时间滞后，洪水期水位上升，洪水过后水位下降，水位升降幅度取决于洪水深度。

图 2所示为静乐井水位与井区降雨量和东碾河流量的对比。由图 2可知，201年8月20日左右单日降雨量为82.8mm时，东碾河流量由原来的0.777m³/s增至9.8m³/s，增加了9.023m³/s；静乐井水位由9.074m升至7.747m，水位上升幅度达1.327m。由此可知，降雨与河流荷载作用是影响静乐井水位动态变化的主要因素，但降雨结束后井水位并未迅速回落至上升前的水平。张淑亮等（1997）曾对静乐井水位高值变化与降雨的关系进行了较系统的研究，认为降雨量的增大虽可在短时间内造成井水位升降变化，但不能改变水位的年变动态，只有在区域应力场增强的背景条件下，降雨引起的附加力源才可能造成井水位年动态的异常变化，降雨荷载是区域应力场变化的诱发因素，也可能是在2017年雨季结束后井水位仍处于高值的主要原因。

图 2 静乐井水位与降雨量和东碾河流量的对比

Figure 2. Comparison of Jingle well level with rainfall and Dongnian river discharge

下载: 全尺寸图片幻灯片

3. 降雨与河流荷载作用下井区应力数值模拟分析

利用井区周围地壳结构、不同地层介质弹性参数建立三维地质模型，预测降雨与河流荷载作用下的井区应力。由于本文基于数值模拟与含水层垂向应力反演的方法进行定量计算，因此，依据静乐井井区地层岩性特征和井-含水层特征参数选取计算参数。因静乐井位于河谷地区，故计算时以变质岩系中大理岩力学性质指标经验数据为依据。

3.1 三维地质模型建立

根据静乐井区地层结构、河流位置及研究目标，设定有限元尺寸为1km×1km×360m，河流宽76m。采用SOLID185单元，模型采用由面拉伸成体的建模方式，取井孔周围围岩弹性模量为5×10¹⁰MPa，泊松比为0.25，密度为2600kg/m³。

网格划分的目的是对模型实现离散化，并用适当数量的网格单元得到最精确的解。本文选用Workblench中的自动网格划分工具，以满足模型要求。网格疏密程度直接影响计算结果的精度，但网格加密将增加CPU计算时间，且需更大的存储空间。为保证模型网格密度和计算精度，本文选择的网格尺寸为1.5m，同时适当增加井区附近区域的网格密度。其他参数均采用默认值，可满足模型的需求。

载荷和约束是ANSYS软件求解计算的边界条件，由所选单元的自由度形势定义。由于本文主要分析河流与降雨荷载的影响，因此选择垂向应力加载。为防止模型移动，模型底端与侧面采用固定约束，限制发生移动和转动。

考虑数值模拟结果的多样性及模型造成的误差，建模时选择不同的边界条件，并反复测试计算，以可反映降雨与最大河流荷载作用对井区的最大影响范围作为模型的边界条件。由于降雨与河流荷载的作用主要表现为垂向应力，因此模型初始地应力只考虑重力。因本文通过数值模拟与反演的方法分析静乐井水位在降雨与河流荷载作用下的垂向应力在量级上是否具有相似性，进而确定影响水位变化的主要因素，对计算精度的要求相对较低，故选取物理力学参数时以研究区内地层岩性的平均值作为模型计算参数。

3.2 结果分析

由于2017年是近年来静乐井水位年变幅最大的年份，也是降雨与河流荷载影响最大的年份，故本文以2017年8月洪水期东碾河水位、静乐井区降雨量和静乐井水位为研究对象，计算降雨和河流荷载作用下的静乐井井区的垂向位移和垂向应力。

（1）单降雨荷载的影响

以2017年8月降雨量作为加载体，将降雨量等效为荷载，计算得到静乐井垂向位移和垂向应力如图 3所示。由图 3可知，降雨荷载对静乐井产生的垂向位移约为0.6mm，垂向应力约为12hPa。

图 3 单降雨荷载作用下静乐井垂向位移和垂向应力

Figure 3. Vertical displacement and vertical stress of Jingle well under single rainfall load

下载: 全尺寸图片幻灯片

（2）单河流荷载的影响

以洪水期东碾河水位最大变化量作为加载体，并将河水水位等效为荷载计算得到静乐井垂向位移和垂向应力如图 3所示。由图 3可知，洪水期河流荷载对静乐井产生的垂向位移约为0.3mm，垂向应力约为15hPa。

图 4 单河流荷载作用下静乐井垂向位移和垂向应力

Figure 4. Vertical displacement and vertical stress of Jingle well under single river load

下载: 全尺寸图片幻灯片

（3）降雨与河流荷载的综合影响

由于静乐井水位变化是在因降雨与河流荷载的共同作用，单河流和单降雨的影响仅能反映井水位的部分变化，因此，需考虑二者的共同作用。以2017年8月东碾河水位最大变化量和降雨量作为共同加载体，并将二者等效为荷载，计算得到静乐井垂向位移和垂向应力如图 5所示。由图 5可知，降雨和河流荷载共同作用对静乐井产生的垂向位移约为0.7mm，垂向应力约为30hPa，与单河流和单降雨影响的二者之和基本一致，表明降雨与河流荷载作用是影响静乐井水位变化的主要因素。

图 5 河流与降雨荷载共同作用下静乐井垂向位移和垂向应力

Figure 5. Vertical displacement and vertical stress of Jingle well under combined action of river and rainfall loads

下载: 全尺寸图片幻灯片

4. 静乐井-含水层垂向应力反演

数值模拟结果表明，降雨和河流荷载是影响静乐井水位变化的主要因素。为进一步验证数值模拟结果能否真实反映井区实际应力变化，还需将其与实际水位观测值所反映的井区应力进行对比。在由井水位变化反演井区应力场变化方面，目前取得了一些进展，如黄辅琼等（2004）利用华北地区40多个深井水位动态变化资料，探讨研究了华北地区现今构造应力场状态；孙小龙等（2011）运用小波分析法提取出能反映水位多年动态变化的趋势信息，反演出华北地区多年构造应力场的变化特征；丁风和等（2015）利用部分含水层介质参数、井水位变化量与含水层垂向应力变化量的关系式，定量分析了张渤带地区构造应力场的动态变化过程。

本文根据静乐井-含水层介质参数及井水位变化量等数据，反演井水位变化引起的含水层垂向应力。应力反演所用数据时段为2016—2017年（根据东碾河水位与流量始测时间确定）。反演前首先对观测数据进行预处理（丁风和等，2015），去除长周期干扰因素、降雨、气压等年际变化对水位的影响；然后采用去趋势法和傅里叶周期法，剔除地下水开采和降水补给等长周期干扰对水位的影响。

4.1 反演所需各类参数计算

（1）气压系数与潮汐效率计算

利用卡尔曼滤波方法消除地球固体潮对所选时段水位、气压整点值数据的影响。利用滤波后的日均值数据，采用高阶差分和一元线性回归等方法获得静乐井气压系数。利用维尼迪科夫潮汐调和分析方法获得静乐井水位M₂波潮汐效率。气压系数和潮汐效率计算结果如表 1所示。

表 1 静乐井-含水层多种参数计算结果

Table 1. Calculation results of various parameters of Jingle well-aquifer

时间	气压系数/（mm·hPa⁻¹）	潮汐效率/（mm·10⁻⁹）	体积压缩系数/pa⁻¹	孔隙度/%	固体体积压缩系数/pa⁻¹	水位变化量/m	含水层垂向应力/hPa
2016-01	4.618	0.191	69.285	1.806	5.117	0.013	1.4019
2016-02	2.522	0.330	9.766	2.442	0.751	0.007	0.5703
2016-03	3.267	0.308	167.63	1.008	6.112	−0.004	−0.8453
2016-04	5.997	0.180	960.22	0.434	18.121	−0.014	−5.9814
2016-05	4.328	0.342	530.25	0.516	10.970	0.005	1.8874
2016-06	5.032	0.264	4926.8	0.053	10.884	−0.0009	−3.2224
2016-07	5.742	0.316	221.25	0.657	6.260	−0.003	−0.8655
2016-08
9.897	0.238	592.66	0.203	5.651	−0.017	−15.5497
2016-09	3.960	0.167	170.08	1.215	8.059	−0.006	−0.9238
2016-10	4.642	0.148	90.599	1.785	6.621	−0.007	−0.8333
2016-11	2.459	0.269	162.6	1.366	6.875	0.002	0.2584
2016-12	3.767	0.147	1709.5	0.106	6.9188	0.008	14.8921
2017-01	4.007	0.154	123.43	1.52	7.348	0.007	0.8916
2017-02	2.861	0.214	254.43	1.089	9.405	0.001	0.1522
2017-03	2.784	0.227	1454.1	0.071	3.449	0.001	1.3134
2017-04	3.808	0.292	148.14	0.997	5.684	−0.007	−1.3365
2017-05	6.093	0.206	600.47	0.494	12.952	0.011	4.2036
2017-06	4.756	0.285	84.983	1.233	4.318	0.005	0.8134
2017-07	6.225	0.249	114.4	1.032	5.173	−0.004	−0.7851
2017-08	18.400	0.191	4233.5	0.018	36.619	−0.03	−323.3068
2017-09	15.148	0.164	751.41	0.388	14.207	−0.015	−7.5280
2017-10	6.703	0.269	3324.2	0.055	8.1277	−0.012	−42.1179
2017-11	3.473	0.161	813.31	0.676	20.424	0.010	2.9897
2017-12	—	0.246	—	—	—	—	—

下载: 导出CSV

| 显示表格

（2）静乐井-含水层介质参数计算

不排水状态下，井水位气压系数和潮汐效率可分别表示为（Bredehoeft，1967；李春洪等，1990；张昭栋等，1993）：

$$ {B_{\rm{p}}} = \frac{{n\beta }}{{\alpha + n\beta }} $$

(1)

$$ {B_g} = - \frac{{1 - n}}{{\rho g\left[ {\left({1 - n} \right)\alpha + \beta } \right]}} $$

(2)

由式（1）、式（2）可得：

$$ {B_g} = - \frac{{1 - n}}{{\rho gng\left[ {\frac{{\left({1 - n} \right)\left({1 - {B_{\rm{p}}}} \right)}}{{{B_{\rm{p}}}}} + 1} \right]}} $$

(3)

式中B_P为井水位气压系数；B_g为M₂波潮汐效率；n 为含水层孔隙度；α为含水层固体骨架体积压缩系数；β为含水层内水的体积压缩系数；ρ为水的密度；g为重力加速度（ρg=0.098hPa/mm）。因此，n、β可依据式（3）与表 1中的潮汐效率和气压系数得到。根据式（1）或式（2）可求α。

4.2 静乐井-含水层垂向应力反演

利用表 1数据，根据水平层状含水层模式下，含水层介质参数（孔隙度、水和固体骨架体积压缩系数等）、井水位变化量与含水层垂向应力变化量的关系式（张昭栋等，1987）等来计算含水层垂向应力，含水层垂向应力计算如下：

$$ \Delta {\sigma _2} = \frac{{ - 2\beta \rho g}}{{n\beta + \left({1 - n} \right)/E}} \bullet \Delta H $$

(4)

式中Δσ_z为含水层垂向应力变化量；E 为含水层固体骨架的杨氏模量（E =1/α）；△ H 为剔除地下水开采、降雨和气压影响后的含水层应力变化引起的压力水头变化量，即井水位变化量。当井-含水层系统所受应力增强，即Δσ_z >0时，井水位上升，水位埋深H变小，其变化量ΔH < 0；当井-含水层系统所受应力减弱，即Δσ_z < 0时，井水位下降，水位埋深H变大，其变化量ΔH>0。利用式（4）与表 1中有关参数计算静乐井2016—2017年含水层垂向应力，结果如表 1、图 6所示。可知，静乐井含水层垂向应力高值时段与降雨量高值时段、东碾河流量高值时段较吻合，再次表明降雨与河流荷载作用是造成含水层垂向应力高值变化的主要因素，井水位的大幅突升变化与垂向应力增大密切相关。

图 6 东碾河流量、静乐井-含水层垂向应力与降雨量的对比

Figure 6. Flow rate of Dongnian river, vertical stress of Jingle well-aquifer and rainfall comparison chart

下载: 全尺寸图片幻灯片

5. 静乐井水位异常分析

5.1 模型预测值与井水位反演值对比分析

由数值模拟结果可知，无论单河流荷载还是单降雨荷载效应均对洪水期静乐井水位造成了一定影响，二者造成的垂向位移累计0.9mm，垂向应力累计27hPa。河流和降雨综合影响模型结果也得到类似结果（垂向位移累计0.7mm，垂向应力累计30hPa）。由含水层垂向应力反演结果可知，洪水期含水层的应力显著增强，由正常时段的4.2hPa增至79hPa左右。2种方法计算结果均表明静乐井水位上升变化与河流荷载效应和降雨荷载效应密切相关。但2种方法计算结果存在一定差异，降雨与河流荷载作用下井区垂向应力数值模拟结果小于含水层垂向应力反演值，相差39hPa。

由于数值模拟计算得到的垂向应力是研究区自身重力、降雨与河流荷载作用下产生的垂向应力之和，而由井水位反演的垂向应力不仅包含上述因素产生的垂向应力，可能还含有一些因构造活动引起的垂向应力变化。事实上2种方法的计算结果的确存在一定差异。虽然2种计算方法均受建模条件和参数选取的影响而存在计算误差，但根据以往的研究，误差值一般不会大于计算值，因此2种计算结果的差异性是客观存在的。按照表 1中静乐井水位平均气压系数4.3521mm/hPa可推算出2种方法垂向应力计算结果差值引起的井水位变化幅度约为16.97cm，大于正常时段水位平均月变化幅度0.9cm的水平。表明影响静乐井水位快速上升的因素除降雨与河流荷载作用外，可能还受其他因素的影响。

5.2 水位高值异常与构造活动的关系

为进一步探讨静乐井水位2017年8月以来的高值异常除受降雨与河流荷载影响外是否受构造活动的影响，本文从可反映区域应力状态变化的参数和测量值（缺震和b值、地震矩释放、GPS基线等）及同一构造区其他前兆测项变化特征与静乐井水位高值异常进行对比分析。

大量研究结果表明，b值与其对应区域应力状态、地壳破裂强度有关（刘雁冰等，2017；朱艾斓等，2005；吴小平，1990）。1983年以来山西北部至晋冀蒙交界区3级以上地震缺震和b值曲线有4次异常过程（见图 7），第1次异常出现在1988年，异常持续过程中在山西地震带出现一组中强地震活动，先后发生了大同-阳高6.1级地震、侯马4.9级和忻州5.1级地震、大同-阳高5.8级地震等，其中大同-阳高6.1级地震被认为是华北地区在经历唐山大地4期张淑亮等：基于数值模拟与含水层垂向应力反演的静乐井水位异常分析863震后进入新活跃幕的标志，是华北应力场趋于增强的结果；第2次异常出现在1995年初，异常持续过程中先后发生了1996年5月包头6.4级地震、1998年1月张北6.2级地震、1999年3月张北5.6级地震、1999年11月大同-阳高5.6级地震等中强地震活动，既是1989年开始的华北新一轮活跃幕的继续，又是华北应力场增强的结果；第3次异常始于2007年，这次异常过程中在山西北部至晋冀蒙交界区发生了几次4级以上中等地震，部分学者认为，这组中等地震的发生是受汶川地震的影响，鄂尔多斯块体与华北平原块体相对挤压和扭错增强，导致山西带形变场与构造应力场由原来的构造拉张转为构造挤压，是汶川地震延迟触发的结果（刘峡等, 2013, 朱艾斓等，2010；刘瑞春等，2014）；2017年又出现同步异常变化，目前异常仍在发展，且缺震和b值的几次异常与静乐井水位高值异常较为同步，因此，可认为静乐井水位高值异常反映的是华北区域应力场增强的过程，是山西地震带地震活动增强的重要判据（张淑亮等，1997）。2017年静乐井水位高值异常也可能与区域应力场增强有一定关系。

图 7 静乐井水位月均值与山西地震带3级地震缺震和b值的对比

Figure 7. Comparison chart of monthly average water level of jingle well with magnitude 3 earthquake deficiency and b value in Shanxi seismic belt

下载: 全尺寸图片幻灯片

由跨过静乐井的山西岢岚-山西太原GPS基线时序曲线图 8可知，2016年以来基线由原来的下降趋势（基线缩短）变为缓慢上升（基线拉长），应力状态在2016年发生了变化，由静乐井所在区域地震矩加速释放时图 9可知，2017年井区存在明显的加速释放低值区。GPS基线和地震矩释放均表明2017年静乐井所在区域应力场在增强。

图 8 山西岢岚-山西太原GPS基线时序曲线

Figure 8. GPS baseline timing curve of Kelan, Shanxi-Taiyuan, Shanxi

下载: 全尺寸图片幻灯片

图 9 静乐井所在区域地震矩加速释放时空图

Figure 9. Accelerated release diagram of seismic moment in jingle well area

下载: 全尺寸图片幻灯片

静乐井构造上位于山西西部吕梁山隆起区，位于同一构造区距静乐井较近的前兆测点有宁武钻孔应变（60km）和神池钻孔应变（95km），在静乐井水位突升前后这2个测点的应变观测值也出现大幅度反向变化，所以说静乐井水位高值变化不是独立事件（见图 10），可能反映了所在构造区应力场的改变。

图 10 静乐井水位与同一构造区其它前兆测项对比图

Figure 10. Comparison between Jingle well water level and other precursor measurements in the same structure area

下载: 全尺寸图片幻灯片

综上所述，静乐井水位的高值异常除与降雨、河流荷载作用有关外，井区构造活动增强可能也是一种因素。

6. 结论与认识

通过数值模拟、井-含水层垂向应力反演、区域应力场变化特征等对静乐井水位2017年以来出现的高值异常进行初步探讨，得出以下结论：

（1）根据井区周围地壳结构、不同地层介质弹性参数建立三维地质模型，并以2017年洪水期东碾河水位最大变化量和降雨量作为共同加载体，模拟得到静乐井垂向位移约为0.7mm，垂向应力约为30hPa。

（2）利用静乐井含水层孔隙度、固体骨架体积压缩系数和水的体积压缩系数等参数，根据水平层状含水层模式下含水层介质参数、井水位变化量与含水层垂向应力变化量的关系式，得到静乐井-含水层垂向应力约为79hPa。

（3）数值模拟与井-含水层垂向应力反演结果均表明，静乐井水位高值异常与同时段降雨量增多、河流荷载效应增强关系密切，但2种结果存在一定差异，含水层垂向应力反演值大于降雨与河流荷载效应，相差39hPa，这种差异可能反映静乐井水位除受降雨与河流荷载作用外，还可能受其他因素的影响。

（4）通过与同一构造区距静乐井较近的前兆测点宁武钻孔应变和神池钻孔应变对比分析可知，静乐井水位突升前后2个测点应变观测值出现大幅度反向变化，表明静乐井水位高值变化不是独立事件，可能反映了所在构造区应力场发生了改变。

（5）可反映区域应力场变化特征的几种参数对2种方法计算结果产生差异的原因表明，2017年静乐井水位高值异常期间山西地震带3级地震缺震和b值、穿过静乐井的GPS基线和地震矩释放均存在显著的异常，推测静乐井水位高值异常除受降雨与河流荷载作用的影响外，也可能受构造活动增强的影响。

图 1 突发事件情景演化的动态贝叶斯网络

Figure 1. Dynamic Bayesian network for the scenario evolution of emergent events

下载: 全尺寸图片幻灯片

图 2 突发事件应急决策方案生成基本流程

Figure 2. Flow chart of generating emergency decision plan for emergency events

下载: 全尺寸图片幻灯片

图 3 地震灾害情景链

Figure 3. Earthquake disaster scenario chains

下载: 全尺寸图片幻灯片

图 4 生命线应急任务分解结构

Figure 4. Breakdown structure for the lifeline emergency task

下载: 全尺寸图片幻灯片

表 1 部分情景条件概率取值

Table 1. Condition probabilities of some scenarios

情景名称	取值	发生概率
建筑物倒塌	极少倒塌	0.350 8
	大半倒塌	0.223 4
	绝大多数倒塌	0.425 7
通信中断	没有中断	0
通信中断	中断	1
交通中断	没有中断	0
交通中断	中断	1
供水、供电、供气设施破坏	设施破坏但短时间能恢复	0.331 8
供水、供电、供气设施破坏	瘫痪	0.668 2
危险品泄漏	没有泄漏	0.662 6
危险品泄漏	泄漏	0.337 4
火灾	极少（10次以下）	0.786 1
火灾	次数多/发生大规模火灾	0.213 9
人员死亡	在可接受死亡人数范围内	0.389 7
人员死亡	大于可接受死亡人数	0.610 3
生活物资缺乏程度	不缺乏	0.129 7
生活物资缺乏程度	缺乏	0.870 3
社会大众恐慌程度	无	0
社会大众恐慌程度	恐慌严重	1
社会公共安全事件严重程度	无	0.456 0
社会公共安全事件严重程度	哄抢、抢劫等	0.544 0
疫情种类	无疫情或传染性小	0.576 5
疫情种类	死亡率高的传染性疾病	0.423 5

下载: 导出CSV

参考文献(60)

[1]	蔡友军, 林均岐, 刘金龙等, 2015. 基于贝叶斯模型的地震直接经济损失快速评估方法研究. 地震工程与工程振动, 35(2): 144—150 Cai Y. J. , Lin J. Q. , Liu J. L. , et al. , 2015. Study on the rapid evaluation method for direct economic loss of earthquake based on Bayesian model. Earthquake Engineering and Engineering Dynamics, 35(2): 144—150. (in Chinese)
[2]	曹海峰, 2018. 非常规突发事件应急预案研究. 北京: 社会科学文献出版社, 126—131 Cao H. F. , 2018. Research on emergency plan of unconventional emergencies with scenario construction method. Beijing: Social Sciences Academic Press (China), 126—131. (in Chinese)
[3]	陈波, 徐超, 温增平, 2015. 三维激光扫描技术在震害调查中的应用. 震灾防御技术, 10(1): 87—94 doi: 10.11899/zzfy20150109 Chen B. , Xu C. , Wen Z. P. , 2015. The applications of 3D laser scanning technology in seismic damage survey. Technology for Earthquake Disaster Prevention, 10(1): 87—94. (in Chinese) doi: 10.11899/zzfy20150109
[4]	陈晨, 俞政, 张新梅, 2014. 考虑应急活动的非常规突发事件情景刻画研究. 中国安全科学学报, 24(10): 163—169 Chen C. , Yu Z. , Zhang X. M. , 2014. Study on characterization of unconventional emergency scenarios considering emergency activities. China Safety Science Journal, 24(10): 163—169. (in Chinese)
[5]	陈迎欣, 李烨, 魏薇, 2018. 云计算环境下基于任务-资源匹配的地震应急资源分配研究. 科技管理研究, 38(19): 202—208 doi: 10.3969/j.issn.1000-7695.2018.19.029 Chen Y. X. , Li Y. , Wei W. , 2018. Research on allocation of earthquake emergency resource based on task-resources matching in cloud computing environment. Science and Technology Management Research, 38(19): 202—208. (in Chinese) doi: 10.3969/j.issn.1000-7695.2018.19.029
[6]	程希萌, 沈占锋, 邢廷炎等, 2016. 基于高分遥感影像的地震受灾建筑物提取与倒损情况快速评估. 自然灾害学报, 25(3): 22—31 Cheng X. M. , Shen Z. F. , Xing T. Y. , et al. , 2016. Damaged building extraction and rapid assessment for earthquake disasters based on high-resolution remote sensing images. Journal of Natural Disasters, 25(3): 22—31. (in Chinese)
[7]	冯海成, 冯翰, 2012. 国际救援队抗灾体适能的调查报告——及对初中生抗灾体适能培养的建议. 中国应急救援, (1): 14—16.
[8]	巩前胜, 2018. “情景—应对”型应急决策中情景识别关键技术研究. 西安: 西安科技大学. Gong Q. S. , 2018. Research on the key techniques of scenario recognition in “Scenario-Response” emergency decision-making. Xi’an: Xi’an University of Science and Technology. (in Chinese)
[9]	姜卉, 黄钧, 2009. 罕见重大突发事件应急实时决策中的情景演变. 华中科技大学学报(社会科学版), 23(1): 104—108 Jiang H. , Huang J. , 2009. The study on the issues of scenario evolvement in real-time decision making of infrequent fatal emergencies. Journal of Huazhong University of Science and Technology (Social Science Edition), 23(1): 104—108. (in Chinese)
[10]	姜卉, 2012. 应急实时决策中的情景表达及情景间关系研究. 电子科技大学学报(社科版), 14(1): 48—52 Jiang H. , 2012. Expression of scenarios and the relationship among scenarios in emergency real-Time decision-Making. Journal of UESTC (Social Sciences Edition), 14(1): 48—52. (in Chinese)
[11]	姜卉, 侯建盛, 2012. 基于情景推演与网络图技术的应急处置模式研究. 数学的实践与认识, 42(24): 21—32 doi: 10.3969/j.issn.1000-0984.2012.24.003 Jiang H. , Hou J. S. , 2012. The study on the emergency disposition mode based on scenario derivation and network chart. Mathematics in Practice and Theory, 42(24): 21—32. (in Chinese) doi: 10.3969/j.issn.1000-0984.2012.24.003
[12]	李晋, 叶涵, 刘智等, 2017. 西南地区承灾体数据采集与动态更新系统. 震灾防御技术, 12(4): 858—869 doi: 10.11899/zzfy20170414 Li J. , Ye H. , Liu Z. , et al. , 2017. Data collection and dynamic update system for hazard-affect buildings in southwestern China. Technology for Earthquake Disaster Prevention, 12(4): 858—869. (in Chinese) doi: 10.11899/zzfy20170414
[13]	李姜, 张合, 刘志辉等, 2021. 基于遥感影像的建筑抗震能力分析——以张家口万全区为例. 震灾防御技术, 16(2): 371—380 doi: 10.11899/zzfy20210216 Li J. , Zhang H. , Liu Z. H. , et al. , 2021. Analysis of seismic capacity of buildings based on remote sensing image——taking Zhangjiakou Wanquan district as an example. Technology for Earthquake Disaster Prevention, 16(2): 371—380. (in Chinese) doi: 10.11899/zzfy20210216
[14]	李勇建, 乔晓娇, 孙晓晨等, 2015. 基于系统动力学的突发事件演化模型. 系统工程学报, 30(3): 306—318 Li Y. J. , Qiao X. J. , Sun X. C. , et al. , 2015. Modeling the evolution of emergency based on system dynamics. Journal of Systems Engineering, 30(3): 306—318. (in Chinese)
[15]	刘铁民, 2010. 危机型突发事件应对与挑战. 中国安全生产科学技术, 6(1): 8—12 Liu T. M. , 2010. Crisis-type emergency response and challenges. Journal of Safety Science and Technology, 6(1): 8—12. (in Chinese)
[16]	卢弋, 陈霖, 冯伟, 2021. 基于案例推理的城市轨道交通应急预警决策. 交通工程, 21(1): 74—79, 85 Lu Y. , Chen L. , Feng W. , 2021. Research on decision-making of urban rail transit emergency warning based on case reasoning. Journal of Transportation Engineering, 21(1): 74—79, 85. (in Chinese)
[17]	马亚杰, 李永义, 韩秀丽, 2007. 基于人工神经网络的地震经济损失评估. 世界地震工程, 23(1): 146—150 doi: 10.3969/j.issn.1007-6069.2007.01.028 Ma Y. J. , Li Y. Y. , Han X. L. , 2007. Earthquake economic loss estimation based on the ANN. World Earthquake Engineering, 23(1): 146—150. (in Chinese) doi: 10.3969/j.issn.1007-6069.2007.01.028
[18]	马祖军, 谢自莉, 2012. 基于贝叶斯网络的城市地震次生灾害演化机理分析. 灾害学, 27(4): 1—5, 24 doi: 10.3969/j.issn.1000-811X.2012.04.001 Ma Z. J. , Xie Z. L. , 2012. Evolution mechanism of earthquake-induced urban disasters based on bayesian networks. Journal of Catastrophology, 27(4): 1—5, 24. (in Chinese) doi: 10.3969/j.issn.1000-811X.2012.04.001
[19]	穆宣社, 2015. 基于地理空间大数据的应急指挥辅助决策平台研究. 测绘通报, (6): 93—96 Mu X. S. , 2015. Research on emergency command decision support platform based on geospatial big data. Bulletin of Surveying and Mapping, (6): 93—96. (in Chinese)
[20]	齐珉, 齐文华, 苏桂武, 2020. 基于新浪微博的2017年四川九寨沟7.0级地震舆情情感分析. 华北地震科学, 38(1): 57—63 Qi M. , Qi W. H. , Su G. W. , 2020.2017 Sichuan Jiuzhaigou M7.0 earthquake sentiment analysis based on Sina Weibo. North China Earthquake Sciences, 38(1): 57—63. (in Chinese)
[21]	郄子君, 2018. 基于关键承灾体的区域复杂灾害情景建模研究. 大连: 大连理工大学. Qie Z. J., 2018. Elements-at-risk as a framework for modeling the complex disaster scenario of hazard-affected regions. Dalian: Dalian University of Technology. (in Chinese)
[22]	郄子君, 荣莉莉, 2020. 面向灾害情景推演的区域模型构建方法研究. 管理评论, 32(10): 276—292 Qie Z. J. , Rong L. L. , 2020. A construction method of hazard-affected region for disaster scenario evolution. Management Review, 32(10): 276—292. (in Chinese)
[23]	盛勇, 孙庆云, 王永明, 2015. 突发事件情景演化及关键要素提取方法. 中国安全生产科学技术, 11(1): 17—21 Sheng Y. , Sun Q. Y. , Wang Y. M. , 2015. Emergency scenario evolution and extraction method of key elements. Journal of Safety Science and Technology, 11(1): 17—21. (in Chinese)
[24]	苏桂武, Rodgers J. , 田青等, 2020. 参与式渭南地震情景构建: 中国地震风险对策的行动研究示范. 地震地质, 42(6): 1446—1473 doi: 10.3969/j.issn.0253-4967.2020.06.012 Su G. W. , Rodgers J. , Tian Q. , et al. , 2020. The participatory construction of a seismic scenario for Weinan city: a pilot action research to address the improvement of earthquake disaster risk reduction in China. Seismology and Geology, 42(6): 1446—1473. (in Chinese) doi: 10.3969/j.issn.0253-4967.2020.06.012
[25]	汤文宇, 2007. CBR的应用研究. 南京: 南京邮电大学.
[26]	王浩林, 范小晶, 2018. 基于情景匹配的海上船舶溢油事故危险性快速评估. 环境工程, 36(S1): 848—853 Wang H. L. , Fan X. J. , 2018. Timely risk assessment of shipping oil spill based on situation matching. Environmental Engineering, 36(S1): 848—853. (in Chinese)
[27]	王慧彦, 王建飞, 张敬军, 2016. 基于情景构建的地震巨灾准备计划框架. 自然灾害学报, 25(6): 112—119 Wang H. Y. , Wang J. F. , Zhang J. J. , 2016. Framework of earthquake catastrophe preparedness planning based on scenarios. Journal of Natural Disasters, 25(6): 112—119. (in Chinese)
[28]	王慧彦, 王建飞, 刘晓静, 2019. 基于情景的巨灾应对协同方法研究——以地震为例. 中国社会公共安全研究报告, (2): 79—97.
[29]	王建飞, 2015. 基于情景平衡的地震灾害应急方案形成研究. 焦作: 河南理工大学. Wang J. F., 2015. Study on the earthquake emergency plan based on scenario balance. Jiaozuo: Henan Polytechnic University. (in Chinese)
[30]	王建飞, 张景发, 2018. 基于情景平衡的地震应急方案生成模型. 安全与环境学报, 18(5): 1923—1930 Wang J. F. , Zhang J. F. , 2018. Model for preparing an earthquake emergency plan based on the background scenario balance. Journal of Safety and Environment, 18(5): 1923—1930. (in Chinese)
[31]	王循庆, 李勇建, 孙华丽, 2014. 基于随机Petri网的群体性突发事件情景演变模型. 管理评论, 26(8): 53—62, 116 Wang X. Q. , Li Y. J. , Sun H. L. , 2014. Scenario evolvement model of unexpected incidents involving mass participation based on stochastic Petri nets. Management Review, 26(8): 53—62, 116. (in Chinese)
[32]	王艳东, 李昊, 王腾等, 2016. 基于社交媒体的突发事件应急信息挖掘与分析. 武汉大学学报·信息科学版, 41(3): 290—297 Wang Y. D. , Li H. , Wang T. , et al. , 2016. The mining and analysis of emergency information in sudden events based on social media. Geomatics and Information Science of Wuhan University, 41(3): 290—297. (in Chinese)
[33]	武旭鹏, 夏登友, 李健行, 2014. 非常规突发事件情景描述方法研究. 中国安全科学学报, 24(4): 159—165 Wu X. P. , Xia D. Y. , Li J. H. , 2014. Study on method for describing unconventional emergency scenario. China Safety Science Journal, 24(4): 159—165. (in Chinese)
[34]	夏登友, 2015. 基于“情景—应对”的非常规突发灾害事故应急决策技术研究. 北京: 北京理工大学. Xia D. Y., 2015. Research on emergency decision-making for unconventional disasters and accidents based on scenario-response. Beijing: Beijing Institute of Technology. (in Chinese)
[35]	许冲, 田颖颖, 沈玲玲等, 2018.2015年尼泊尔廓尔喀M_W7.8地震滑坡数据库. 地震地质, 40(5): 1115—1128 Xu C. , Tian Y. Y. , Shen L. L. , et al. , 2018. Database of landslides triggered by 2015 Gorkha (Nepal) M_W7.8 earthquake. Seismology and Geology, 40(5): 1115—1128. (in Chinese)
[36]	杨保华, 方志耕, 刘思峰等, 2012. 基于GERTS网络的非常规突发事件情景推演共力耦合模型. 系统工程理论与实践, 32(5): 963—970 doi: 10.3969/j.issn.1000-6788.2012.05.007 Yang B. H. , Fang Z. G. , Liu S. F. , et al. , 2012. Model of co-coupling in unconventional incidents based on GERTS network. Systems Engineering-Theory & Practice, 32(5): 963—970. (in Chinese) doi: 10.3969/j.issn.1000-6788.2012.05.007
[37]	杨继君, 佘廉, 2016. 面向多灾点需求的应急资源调度博弈模型及优化. 中国管理科学, 24(8): 154—163 Yang J. J. , She L. , 2016. Game model and optimization based on resource requirements of multiple crisis locations. Chinese Journal of Management Science, 24(8): 154—163. (in Chinese)
[38]	杨青, 杨帆, 2012. 基于元胞自动机的突发传染病事件演化模型. 系统工程学报, 27(6): 727—738 doi: 10.3969/j.issn.1000-5781.2012.06.002 Yang Q. , Yang F. , 2012. Emergency epidemics spread model using cellular automata. Journal of Systems Engineering, 27(6): 727—738. (in Chinese) doi: 10.3969/j.issn.1000-5781.2012.06.002
[39]	于山, 王海霞, 马亚杰, 2005. 三层BP神经网络地震灾害人员伤亡预测模型. 地震工程与工程振动, 25(6): 113—117 doi: 10.3969/j.issn.1000-1301.2005.06.019 Yu S. , Wang H. X. , Ma Y. J. , 2005. Three-layer BP network model for estimation of casualties in an earthquake. Earthquake Engineering and Engineering Dynamics, 25(6): 113—117. (in Chinese) doi: 10.3969/j.issn.1000-1301.2005.06.019
[40]	袁晓芳, 田水承, 王莉, 2011. 基于PSR与贝叶斯网络的非常规突发事件情景分析. 中国安全科学学报, 21(1): 169—176 doi: 10.3969/j.issn.1003-3033.2011.01.028 Yuan X. F. , Tian S. C. , Wang L. , 2011. Scenario analysis of unconventional emergency based on PSR model and Bayesian networks. China Safety Science Journal, 21(1): 169—176. (in Chinese) doi: 10.3969/j.issn.1003-3033.2011.01.028
[41]	张磊, 王延章, 陈雪龙, 2016. 基于知识元的非常规突发事件情景模糊推演方法. 系统工程学报, 31(6): 729—738 Zhang L. , Wang Y. Z. , Chen X. L. , 2016. Fuzzy inference method for unconventional events scenarios based on knowledge unit. Journal of Systems Engineering, 31(6): 729—738. (in Chinese)
[42]	张明红, 佘廉, 2016. 基于情景的突发事件演化模型研究——以青岛“11.22”事故为例. 情报杂志, 35(5): 65—71 Zhang M. H. , She L. , 2016. A case study on the evolutionary model of emergency events based on scenario: the Qingdao "11.22" accident. Journal of Intelligence, 35(5): 65—71. (in Chinese)
[43]	张明红, 2016. 基于案例的非常规突发事件情景推理方法研究. 武汉: 华中科技大学. Zhang M. H., 2016. Study on unconventional emergency scenario reasoning method the cases-based. Wuhan: Huazhong University of Science and Technology. (in Chinese)
[44]	张文娟, 2019. 基于移动通信大数据的地震灾区人口伤亡获取与评估系统设计. 地震工程学报, 41(4): 1066—1071 doi: 10.3969/j.issn.1000-0844.2019.04.1066 Zhang W. J. , 2019. Design of the population casualty acquisition and evaluation system in earthquake disaster areas based on mobile communication big data. China Earthquake Engineering Journal, 41(4): 1066—1071. (in Chinese) doi: 10.3969/j.issn.1000-0844.2019.04.1066
[45]	张莹, 郭红梅, 尹文刚等, 2019. 基于SIFT特征与SVM分类的地震灾情图像信息异常检测方法. 地震研究, 42(2): 265—272 doi: 10.3969/j.issn.1000-0666.2019.02.015 Zhang Y. , Guo H. M. , Yin W. G. , et al. , 2019. Detection method of earthquake disaster image anomaly based on SIFT feature and SVM classification. Journal of Seismological Research, 42(2): 265—272. (in Chinese) doi: 10.3969/j.issn.1000-0666.2019.02.015
[46]	张英菊, 2016. 案例推理技术在环境群体性事件应急决策中的应用研究. 安全与环境工程, 23(1): 94—99. Zhang Y. J. , 2016. Research on application of CBR in emergency decision-making for environmental massive events. Safety and Environmental Engineering, 23(1): 94—99. (in Chinese)
[47]	张媛, 王金萍, 赖俊彦等, 2020. 互联网+社区地震应急桌面演练的设计与实施——以青海省海西州为例. 震灾防御技术, 15(2): 411—418 Zhang Y. , Wang J. P. , Lai J. Y. , et al. , 2020. Design and implementation of internet + community earthquake emergency tabletop exercises——a case study of Haixi prefecture in Qinghai province. Technology for Earthquake Disaster Prevention, 15(2): 411—418. (in Chinese)
[48]	仲秋雁, 郭艳敏, 王宁等, 2012. 基于知识元的非常规突发事件情景模型研究. 情报科学, 30(1): 115—120 Zhong Q. Y. , Guo Y. M. , Wang N. , et al. , 2012. Research on unconventional emergency scenario model based on knowledge element. Information Science, 30(1): 115—120. (in Chinese)
[49]	朱达邈, 2017. 地震灾评推演训练中主要情景的客观度优化研究. 北京: 中国科学院大学. Zhu D. M., 2017. Research on objective optimization of main scenes in earthquake disaster assessment deduction training. Beijing: University of Chinese Academy of Sciences. (in Chinese)
[50]	Cosgrave J. , 1996. Decision making in emergencies. Disaster Prevention and Management, 5(4): 28—35. doi: 10.1108/09653569610127424
[51]	Fahey L. , 1999. Competitor scenarios: projecting a rival’s marketplace strategy. Competitive Intelligence Review, 10(2): 65—85. doi: 10.1002/(SICI)1520-6386(199932)10:2<65::AID-CIR9>3.0.CO;2-Z
[52]	Gadomski A. M. , Bologna S. , Costanzo G. D. , et al. , 2001. Towards intelligent decision support systems for emergency managers: the IDA approach. International Journal of Risk Assessment and Management, 2(3—4): 224—242.
[53]	Gilbert A. L. , 2000. Using multiple scenario analysis to map the competitive futurescape: a practice-based perspective. Competitive Intelligence Review, 11(2): 12—19. doi: 10.1002/(SICI)1520-6386(200032)11:2<12::AID-CIR4>3.0.CO;2-E
[54]	Kahn H., Wiener A. J., 1967. The year 2000: a framework for speculation on the next thirty-three years. New York: MacMillan and Co Limited.
[55]	Kowalski-Trakofler K. M. , Vaught C. , Scharf T. , 2003. Judgment and decision making under stress: an overview for emergency managers. International Journal of Emergency Management, 1(3): 278—289. doi: 10.1504/IJEM.2003.003297
[56]	Li X. L. , Li Z. Q. , Yang J. S. , et al. , 2020. Seismic vulnerability comparison between rural Weinan and other rural areas in Western China. International Journal of Disaster Risk Reduction, 48: 101576. doi: 10.1016/j.ijdrr.2020.101576
[57]	Liu Y. H. , So E. , Li Z. Q. , et al. , 2020. Scenario-based seismic vulnerability and hazard analyses to help direct disaster risk reduction in rural Weinan, China. International Journal of Disaster Risk Reduction, 48: 101577. doi: 10.1016/j.ijdrr.2020.101577
[58]	Rodgers J. , Su G. W. , Qi W. H. , et al. , 2020. Creating an earthquake scenario in China: a case study in Weinan City, Shaanxi Province. International Journal of Disaster Risk Reduction, 42: 101305. doi: 10.1016/j.ijdrr.2019.101305
[59]	Wei Y. M. , Liang Q. M. , Fan Y. , 2006. A scenario analysis of energy requirements and Energy Intensity for China’s rapidly developing society in the year 2020. Technological Forecasting and Social Change, 73(4): 405—421. doi: 10.1016/j.techfore.2004.12.003
[60]	Xu C., Dai F. C., Xu X. W., et al., 2012. GIS-based support vector machine modeling of earthquake-triggered landslide susceptibility in the Jianjiang River watershed, China. Geomorphology, 145—146: 70—80.