Risk Assessment of Earthquake Disaster Impact on Industrial Chain
-
摘要: 保障产业链安全稳定、具备抗灾韧性,成为我国经济高质量发展的重要基础,因此基于已收集整理的全球大震巨灾对产业链冲击影响资料,梳理分析地震灾害对不同行业及产业链的影响,尤其是针对地震巨灾对产业链影响的特征及要素进行研究。基于全球地震巨灾对产业链的冲击影响研究,通过量化方式给出产业链地震灾害风险评估水平及不同产业链震后恢复能力,并提出产业链地震安全保障举措建议,为研究我国产业链地震敏感性提供参考。Abstract: Ensuring the safety, stability, and resilience of the industrial chain is crucial for the high-quality development of China's economy. Therefore, based on the collected and organized data on the impact of global earthquake disasters on the industrial chain, this study examines the impact of earthquake disasters on the different industry and industrial chain, especially the characteristics and elements of the impact of earthquake disasters on the industrial chain. Drawing from the examination of the impact of global earthquake disasters on the industrial chain, this study quantitatively provides the level of earthquake risk assessment for the industrial chain, along with the post-earthquake recovery capabilities of different industrial chains. Additionally, this study presents measures and suggestions for ensuring earthquake safety within the industrial chain, providing a basis for investigating the earthquake sensitivity of China's industrial chain in China.
-
Key words:
- Industrial chain /
- Earthquake disasters /
- Natural disasters /
- Risk assessment /
- Resilience
-
引言
据河北省地震台网测定,2018年2月12日18时31分在河北省廊坊市永清县(39.37°N,116.67°E)发生4.3级地震,震源深度20km。此次地震未造成人员伤亡,在震中区域造成个别房屋墙皮脱落,河北省廊坊市、固安市震感较强,北京市、天津市普遍有感。地震发生后,中国地震局工程力学研究所对永清4.3级地震强震动记录进行了远程和现场记录回收。强震动记录是地震能量在地表或建筑结构上的真实反映,是工程地震、地震学研究的基础数据资料(朱升初等,2018),本文对此次地震构造背景及强震动记录进行了分析,对永清4.3级地震烈度采用建立虚拟台站方法进行了校正,为河北地区烈度调查、地震区划等提供基础资料。
1. 发震构造背景
永清4.3级地震区域构造主要以NE、NW向断裂为主,地震发生在廊固凹陷内(图 1),廊固凹陷是冀中坳陷北部的一个古近系湖盆,整个凹陷大体呈NE走向,南北长约90km,东西宽20—40km,面积约2600km2。该凹陷是受区域伸展并伴有走滑作用形成的典型箕状凹陷,整个凹陷表现为北断南超、西断东超的构造格局。由于受NE向断裂构造控制,表现出整体不均匀下降的箕状断陷特征。该凹陷内曾于1621年发生廊坊永清$ 5\frac{1}{2} $级地震(国家地震局震害防御司,1995)。
图 1 永清4.3级地震构造背景Figure 1. The seismotectonic background map of the Yongqing MS4.3 earthquake(1)施庄断裂(2)南口山前断裂(3)南口-孙河断裂(4)小汤山-东北旺断裂(5)黄庄-高丽营断裂 (6)顺义-前门-良乡断裂(7)南苑-通县断裂(8)新夏垫断裂(9)夏垫断裂(10)大兴凸起东缘断裂 (11)蓟县山前断裂(12)宝坻断裂(13)桐柏断裂(14)河西务断裂(15)天津断裂(16)大城东断裂 (17)沧东断裂(18)海河断裂(19)牛驼镇凸起东缘断裂(20)徐水断裂(21)东垒子-涞水断裂 (22)永定河断裂(23)徐水南断裂本次地震发生在廊固凹陷东侧控制断裂河西务断裂带附近,该断裂北起河西务西北,向西南经大王务,延伸至永清东码头镇,全长近50km,走向NE,倾向SE,为上陡下缓的铲形正断层,断层产状与震源机制参数基本一致(王晓山等,2018),所以永清4.3级地震发震构造可能为河西务断裂。周月玲等(2018)通过浅层地震勘探和钻孔联合剖面探测,结合年代样品测试,对河西务断裂活动性进行了综合研究,揭示断裂上断点埋深约150m或以浅,第四系底界面垂直断错为20—45m,断裂最新活动时代为晚更新世早期。
2. 观测台站及记录获取
2.1 强震动台站分布
本次地震强震动数据主要由中国地震局工程力学研究所提供,共获取3个分量自由场记录74组,常规处理后得到222条3分向加速度记录,获取记录的强震动台站均为土层台,记录台站分布如图 2所示,由图 2可知,触发台站主要分布在震中北部,南部无台站触发,这可能是因为首都圈台站较密集,而河北省境内强震动台站稀疏,或与仪器型号不同且很多台站老旧未改造有关。触发台站的震中距为29.4—185.5km,其中震中距<50km的记录有5组,震中距为50—100km的有29组,震中距为100—150km的有37组,震中距为150—200km的有3组,距震中最近的台站为长子营台(11ZZY),震中距为29.4km。虽然此次地震获取的强震动记录较全面,但由于震级较小,未造成大的灾害损失,且震中附近台站相对稀疏,能供工程使用的强震动记录较少。
2.2 强震动记录获取
在强震动记录使用分析前,首先对波形数据进行消除基线漂移,即原始强震动记录通过减去事前记录平均值的方法进行零线调整(吕俊强等,2015),通过分析计算得到震中距R<50km的5个台站强震动记录基本信息及相关参数,见表 1。其中距震中最近的长子营台记录到的最大峰值加速度为7.9cm/s2,觅子店台记录到的东西向最大峰值加速度为36.1cm/s2、南北向最大峰值加速度为40.4cm/s2,红星中学台记录到的垂直方向最大峰值加速度为10.8cm/s2,5个台站加速度时程曲线如图 3。由表 1可知,离震中较近的长子营台(震中距29.4km)较离震中较远的觅子店台(震中距38.4km)和于家务台(震中距36.8km)各分量加速度小,影响地震动强度的因素除震中距外,还包括仪器型号、场地特征条件和地质构造、场地土层结构、覆盖层厚度、地形及断裂分布等(王海云等,2010;温瑞智等,2013;郭秋娜等,2018)。本文根据温瑞智等(2015)提出的适用于我国强震动台站场地划分经验曲线,采用H/V谱比法对本次地震触发台站进行分类,分类结果见图 4。由图 4可知,记录到的本次地震台站大部分为Ⅰ类场地(51.35%),其次为Ⅱ类场地(35.14%)和Ⅲ类场地(13.51%),其中,觅子店台和红星中学台为Ⅰ类场地,于家务台和张家湾中学台为Ⅱ类场地,长子营台为Ⅲ类场地。长子营台与觅子店台和于家务台记录差异较大,这可能与场地类型有关。
表 1 强震动记录(震中距<50km)及相关参数Table 1. Strong motion records(epicenter distance < 50km) and related parameters台站名称 台站代码 仪器型号 震中距/km PGA/cm·s-2 EW SN UD 觅子店 11MZD REFTEK/SLJ-100 38.4 36.1 40.4 8.5 红星中学 11YHZ REFTEK/SLJ-100 43.4 -5.4 -4.6 10.8 于家务 11YJW ETNA/ES-T 36.8 -21.7 -18.5 -8.1 张家湾中学 11ZJW REFTEK/SLJ-100 47.9 -5.8 4.2 -4.4 长子营 11ZZY REFTEK/SLJ-100 29.4 5.0 -7.9 -3.5 3. 谱分析
加速度反应谱是目前国内外工程抗震设计的重要依据,反应谱形状、幅值及卓越周期的差异体现了场地条件、震级对地面运动反应谱的影响(温瑞智等,2013;徐钦等,2017)。为避免地震时建筑结构发生共振,确保结构安全,有必要对加速度反应谱成分进行分析。本文选取震中距<50km的5个典型强震动台站进行频谱分析,阻尼比取0.05,得到本次地震EW、SN、UD向加速度反应谱(图 5)。
由图 5可知,觅子店台记录到了本次地震EW向和SN向加速度反应谱最大值,红星中学台记录到了本次地震SN向加速度反应谱最大值。觅子店台与长子营台在EW和SN向差别最大,觅子店台与张家湾中学台在SN向差别最大。5个典型台站记录到的加速度反应谱高频成分丰富,谱卓越周期主要为0.1—0.3s,0.3s后迅速下降。根据现场调查,震区主要房屋类型为土木结构、木房顶砌体结构、预制板砌体结构、钢结构简易房屋4种,本次地震除少数土木结构房屋原有裂缝扩大、墙皮开裂外,其他类型房屋未发现新开裂、屋顶掉瓦、墙皮掉灰现象,中小城市主要建筑物结构自振周期为0.3—1s(王文才等,2018),因此本次地震对该震区自振周期范围内的建筑物破坏影响较小。通过对比不同震中距台站各方向反应谱特征周期发现,震中距对特征周期的影响不明显。
傅氏谱能反映地震波在频域上的幅值特性,地震动记录中包含多种振动频率,不同地区的主要振动频率不同,结构自身具有不同频率、幅频响应,因而在不同频率振动作用下,结构振动不同。为避免地震作用下结构发生共振现象,确保结构安全,有必要对地震频率成分进行分析(冉志杰等,2012)。
本文计算了永清4.3级地震震中距50km以内的5个典型台站强震动记录傅氏谱(图 6),由图 6可知典型台站傅氏谱谱型主要为单峰值型,强震动记录EW分量傅氏谱最大峰值为3.5—8.5Hz,SN分量傅氏谱最大峰值为2.9—6.2Hz,UD分量傅氏谱最大峰值为2.1—9.8Hz。总体来看,UD分量频率范围分布较宽,EW分量傅氏谱最大峰值大于SN分量和UD分量傅氏谱峰值。一般而言,随着震中距的增大,地震动记录长周期(低频)越显著,但本次地震5个典型台站特征不明显,可能与仪器型号及场地特征条件有关。另外,EW分量傅氏谱最大峰值随着震中距的增大呈先减小后增大的趋势,而SN分量和UD分量无该特征。
本文以东西向为正方向,统计峰值加速度反应谱与永清4.3级地震触发台站同震中所处方位角和震中距的关系(图 7)。由图 7可知,峰值加速度反应谱在同一震中距下整体随着方位角的增大而增大,达70°左右后下降,这可能与附近河西务断裂构造走向有关;峰值加速度反应谱值随着震中距的增加而减小。
4. 地震烈度分析
地震烈度是衡量地震破坏和地震动强弱程度的重要指标。Wald等(1999)认为当烈度较小时,用峰值加速度计算烈度较合适,因为人对加速度较敏感。随着强震动观测技术的发展,由于以各观测点峰值加速度作为依据可真实反映地震发生时各点地面运动强度(田秀丰等,2013),基于地震动记录快速确定地震仪器烈度,越来越受到研究者的重视(崔建文等,2016)。目前国内外基于强震动记录资料的地震烈度计算方法较多,不同方法利用不同的地震动参数,包括时间、频度和强度等,建立地震烈度与地震动参数之间的关系。本文采用美国ShakeMap烈度计算方法(李俊等,2010),计算永清4.3级地震事件记录台站烈度,并采用克里格插值方法进行网格化插值,绘制等值线图(图 8)。Wald等(1999)的研究认为,当烈度低于Ⅴ度时峰值加速度PGA与烈度相关性明显高于峰值速度PGV,而当烈度为Ⅴ—Ⅶ度时,PGA和PGV与烈度相关程度相当,考虑本次地震烈度较低,故采用峰值加速度计算烈度。由图 9可知,本次地震仪器烈度等震线长轴沿NW走向,且极震区位于震中东北部,而不是震中附近,但由现场宏观烈度调查(图 8)与震源机制解(王晓山等,2018)可知该地震等震线长轴沿NE向展布,这种差异现象的出现可能与台站分布不均匀有关,震中南部无强震动记录,所以图 8中震中北部附近出现数据空白现象,不利于获取实际烈度长轴方向真实性。
可靠快速生成地震烈度可为人员伤亡估计、应急救援决策和工程抢险修复决策提供重要依据,为解决由于缺少地震动观测值出现空白区域现象及烈度长轴方向和实际方向不一致的问题,本文首先利用永清4.3级地震中获取的所有台站加速度峰值,按地震动峰值衰减模型(王国新,2001)确定该地区地震动衰减关系,然后在整个地震区内建立空间随机假设台站进行补点插值,得到新的地震烈度分布图。地震动峰值衰减模型(王国新,2001)如下:
$$ \lg Y = {{\rm{C}}_1} + {{\rm{C}}_2}M + {{\rm{C}}_3}M\lg \left({R + {R_0}} \right) + {{\rm{C}}_4}\lg \left({R + {R_0}} \right) \pm \varepsilon $$ (1) 式中,Y为PGA(峰值加速度),M为震级,R为震中距,C1、C2、C3、C4为待求常数,R0为近场饱和因子,ε为统计分析误差项。采用通用全局优化法拟合得到地震动峰值加速度随震级和震中距的变化规律,考虑本次地震震级较小,故选取震中距<100km的Ⅰ类场地强震动数据进行拟合,其中水平向按EW、NS向加速度峰值矢量拟合回归分析,拟合曲线如图 9,计算得到的衰减关系式为:
$$ \lg Y = 1.6882 + 0.3926M - 0.1546M\lg \left({R + 5} \right) + 0.7000\lg \left({R + 5} \right) \pm 0.2287 $$ (2) 峰值加速度随震中距的变化如图 10所示,由图 10可知峰值加速度随着震中距的增加而减小。同一震中距不同台站测得的峰值加速度离散性较大。通过拟合各台站实际PGA值,与俞言祥等(2013)地震动衰减关系进行对比,可知20km范围内2条曲线吻合较好,而20km外拟合曲线呈衰减相对缓慢状态,这可能与本次地震震级较小、记录台站均为土层台、远场有一定放大作用有关(王海云,2011)。由观测值得到的预测方程虽存在一定离散性,但整体趋势与实际结果保持一致。
在震区空白区域随机产生N个空间坐标的假设台站,按式(2)计算水平峰值加速度,按美国ShakeMap烈度计算方法转换成仪器烈度,绘制增加假设台站的烈度分布图(图 11),由图 11可知,增加假设台站的仪器烈度分布图与现场烈度调查分布图极震区方向基本一致,进一步验证了河西务断裂为本次地震发震构造的可能,解决了台站分布不均匀导致地震烈度影响场计算出现缺值现象的问题,为以后缺少台站记录震区提供准确快速制作烈度分布图的思路,为震害调查和地震应急救援提供重要依据。
5. 结论
本文介绍了河北永清4.3级地震发震构造背景,并对国家强震动台网中心在本次地震中获取的加速度记录特征进行了分析,为河北省地震动特征和震害调查提供了一定参考,结论如下:
(1) 河北永清4.3级地震在廊固凹陷东侧控制断裂河西务断裂带附近,地震烈度极震区走向与该断裂展布方向一致,结合震源机制解初步判定本次地震发震构造为河西务断裂。
(2) 74个强震动台站获取到强震动加速度记录,震中距位于29.4—185.5km,其中觅子店台记录到的本次地震东西向最大峰值加速度为36.1cm/s2,垂直方向最大加速度为40.4cm/s2,于家务台记录到的南北向最大峰值加速度为18.5cm/s2。采用H/V谱比法对触发台站进行分类,可知Ⅰ类场地居多。
(3) 通过对5个典型台站反应谱分析发现永清4.3级地震加速度反应谱高频成分丰富,谱卓越周期主要为0.1—0.3s。傅氏谱谱型主要为单峰值型,UD分量频率范围分布较宽,EW分量傅氏谱最大峰值大于SN分量和UD分量。
(4) 通过分析峰值加速度反应谱与震中距、方位角的关系,发现峰值加速度反应谱整体随着方位角的增大而增大,达70°左右后下降,这可能与附近河西务断裂构造走向有关,有待进一步研究。其次,峰值加速度反应谱值整体随着震中距的增加而减小。
(5) 由本次地震仪器烈度分布图可知,台站分布不均匀性直接影响烈度分布的真实性,为避免此类问题,拟合出该地区地震动衰减关系,随机产生空间假设台站进行弥补,为以后缺少台站记录震区提供准确快速制作烈度分布图的思路,为震害调查和地震应急救援提供重要依据。由于收集到的强震动数据有限,文中仅采用本次地震获取的峰值加速度拟合当地地震动衰减关系,拟合结果可能存在误差,今后将收集更多的强震动数据进行拟合,以校正该地区地震动衰减关系。
-
表 1 国民经济行业分类及代码
Table 1. Classification and codes of the national economic industries
代码 行业 简称 A 农、林、牧、渔业 农林牧渔 B 采矿业 采矿 C 制造业 制造 D 电力、燃气及水的生产和供应业 水电气暖 E 建筑业 建筑 F 交通运输、仓储和邮政业 运输仓储 G 信息传输、计算机服务和软件业 信息服务 H 批发和零售业 批发零售 I 住宿和餐饮业 住宿餐饮 J 金融业 金融 K 房地产业 房地产 L 租赁和商务服务业 租赁商务 M 科学研究、技术服务和地质勘查业 科技服务 N 水利、环境和公共设施管理业 公共设施 O 居民服务和其他服务 其他服务 P 教育 教育 Q 卫生、社会保障和社会福利业 卫生社保 R 文化、体育和娱乐业 文体娱乐 S 公共管理和社会组织 公共管理 T 国际组织 国际组织 表 2 国民经济行业产业链分类及代码
Table 2. Classification and codes of the national economic industry chain
代码 名称 代码 名称 代码 名称 a 农业 i 造纸 q IT b 林业 j 煤炭 r Internet c 畜牧渔 k 机械制造 s 教育 d 农工贸 l 汽车 t 文体娱乐 e 肉类 m 钢铁 u 旅游 f 蔬菜 n 电信 v 金融 g 医药 o 服装 w 媒介 h 化工 p 高新技术 表 3 全球地震巨灾研究震例
Table 3. Some case examples of global earthquake catastrophe
编号 地震日期/(年-月-日) 纬度 经度 震中地名 深度/km 震级M/级 震中烈度 海啸/m 死亡人数/万人 损失/亿美元 1 1868-08-13 18.60ºS 71.00ºW 智利阿里卡近海 25 8½ 11 18.0 2.50 3.0 2 1908-12-28 38.15ºN 15.68ºE 意大利西西里岛 10 7¼ 11 13.0 8.22 1.2 3 1923-09-01 35.10ºN 139.50ºE 日本横滨 35 7.9 — 13.0 14.28 6.00 4 1939-01-25 36.20ºS 72.20ºW 智利奇廉 60 7.7 10 — 3.00 9.2 5 1960-02-29 30.50ºN 9.50ºW 摩洛哥阿加迪尔 — 6.8 10 — 1.31 1.2 6 1970-05-31 9.25ºS 78.84ºW 秘鲁卡斯马近海 73 7.9 10 0.4 6.68 5.3 7 1972-12-23 12.35ºN 86.13ºW 尼加拉瓜马拉瓜 7 6.2 — — 1.00 8.0 8 1976-02-04 15.30ºN 89.15ºW 危地马拉城 12 7.5 9 0.1 2.30 21.5 9 1976-07-27 39.61ºN 117.89ºE 中国河北唐山 16 7.8 11 — 24.20 122.3 10 1988-12-07 40.93ºN 44.11ºE 亚美尼亚斯皮塔克 5 6.8 10 — 2.50 162.0 11 1990-06-20 37.00ºN 49.19ºE 伊朗加兹温 19 7.9 7 — 4.00 80.0 12 1999-08-17 40.77ºN 30.00ºE 土耳其伊斯坦布尔 13 7.6 10 2.5 1.71 200.0 13 2001-01-26 23.39ºN 70.23ºE 印度普杰 16 7.7 10 — 2.00 26.2 14 2004-12-26 3.30ºN 95.98ºE 印尼锡默卢岛 30 9.1 10 50.9 22.79 100.0 15 2008-05-12 30.95ºN 103.40ºE 中国四川汶川 19 8.0 11 — 8.77 1 207.3 16 2010-01-12 18.44ºN 72.57ºW 海地太子港 13 7.3 10 3.0 22.26 80.0 17 2011-03-11 38.32ºN 142.37ºE 日本宫城海域 32 9.0 12 24.0 2.78 1234.0 18 2023/02/06 37.15ºN 36.95ºE 土耳其南部 20 7.8 11 — 4.50 841.0 表 4 全球地震巨灾对行业的冲击影响统计结果
Table 4. Statistical results of the impact of some global earthquake catastrophes on the industry
序号 地震日期/
(年-月-日)震中地名 震级M/级 行业代码 A B C D E F G H I J K L M N O P Q R S T 1 1868-08-13 智利阿里卡近海 8½ ● — — — ● ● — ● ● ● ● — — — — ● — — ● — 2 1908-12-28 意大利西西里岛 7¼ ● — — — ● ● — ● ● ● ● — — — — ● — — ● — 3 1923-09-01 日本横滨 7.9 ● — — — ● ● — ● ● ● ● — — — — ● ● — ● — 4 1939-01-25 智利奇廉 7.7 ● — — — ● ● — ● ● ● ● — — — — ● ● — ● — 5 1960-02-29 摩洛哥阿加迪尔 6.8 — — — ● ● ● — ● ● ● ● — — — — ● ● ● ● — 6 1970-05-31 秘鲁卡斯马近海 7.9 ● — — ● ● ● — ● ● ● ● — — ● — ● ● ● ● — 7 1972-12-23 尼加拉瓜马拉瓜 6.2 — — — ● ● — — ● ● ● ● — — — — ● ● ● ● — 8 1976-02-04 危地马拉城 7.5 ● — — ● ● ● — ● ● ● ● — — ● ● ● ● ● ● — 9 1976-07-27 中国河北唐山 7.8 ● ● ● ● ● ● — ● ● ● ● ● ● ● ● ● ● ● ● — 10 1988-12-07 亚美尼亚斯皮塔克 6.8 — — — ● ● ● — ● ● ● ● ● — ● ● ● ● ● ● ● 11 1990-06-20 伊朗加兹温 7.9 ● — ● ● ● ● — ● ● ● ● ● ● ● ● ● ● ● ● ● 12 1999-08-17 土耳其伊斯坦布尔 7.6 ● — ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● 13 2001-01-26 印度普杰 7.7 ● — ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● 14 2004-12-26 印尼锡默卢岛 9.1 ● — ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● 15 2008-05-12 中国四川汶川 8.0 ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● 16 2010-01-12 海地太子港 7.3 ● — ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● 17 2011-03-11 日本宫城海域 9.0 ● — ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● 18 2023-02-06 土耳其南部 7.8 ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● ● 合计 15 3 9 14 18 17 7 18 18 18 18 10 9 12 11 18 16 14 18 9 表 5 影响等级划分
Table 5. Impact level classification
严重影响 中度影响 轻度影响 轻微影响 基本无影响 I>0.8 0.8≥I>0.6 0.6≥I>0.4 0.4≥I>0.2 0.2≥I 表 6 全球地震巨灾对产业链的冲击影响统计结果
Table 6. Statistical results of the impact of some global earthquake catastrophes on the industry chain
序号 地震时间/
(年-月)震中地名 震级M/级 产业链代码 a b c d e f g h i j k l m n o p q r s t u v w 1 1868-08 智利阿里卡近海 8½ ● — ● ● ● ● — — — — — — — — — — — — ● — — ● — 2 1908-12 意大利西西里岛 7¼ ● — ● ● ● ● — — — — — — — — — — — — ● — — ● — 3 1923-09 日本横滨 7.9 ● — ● ● ● ● ● — — — — — — — — — — — ● — — ● — 4 1939-01 智利奇廉 7.7 ● — — ● ● ● — — — — — — — ● — — — — ● — — ● — 5 1960-02 摩洛哥阿加迪尔 6.8 ● — — ● ● ● ● — — — — — — ● — — — — ● ● — ● — 6 1970-05 秘鲁卡斯马近海 7.9 ● — ● ● ● ● ● — — — — — — ● — — — — ● ● — ● — 7 1972-12 尼加拉瓜马拉瓜 6.2 ● — — ● ● ● ● — — — — — — ● — — — — ● ● — ● — 8 1976-02 危地马拉城 7.5 ● — ● ● ● ● ● — — — — — — ● — — — — ● ● — ● — 9 1976-07 中国河北唐山 7.8 ● ● — ● ● ● ● ● ● ● ● — ● ● ● — — — ● ● — ● — 10 1988-12 亚美尼亚斯皮塔克 6.8 ● — — ● ● ● ● — — — — — — ● — — — — ● ● — ● — 11 1990-06 伊朗加兹温 7.9 ● — — ● ● ● ● — — — ● — — ● — — — — ● ● — ● — 12 1999-08 土耳其伊斯坦布尔 7.6 ● — ● ● ● ● ● — — — ● — — ● ● — ● ● ● ● ● ● ● 13 2001-01 印度普杰 7.7 ● — — ● ● ● ● — — — ● — — ● ● ● ● ● ● ● ● ● ● 14 2004-12 印尼锡默卢岛 9.1 ● — ● ● ● ● ● — — — ● — — ● ● ● ● ● ● ● ● ● ● 15 2008-05 中国四川汶川 8.0 ● ● ● ● ● ● ● ● — — ● ● ● ● ● ● ● ● ● ● ● ● ● 16 2010-01 海地太子港 7.3 — — ● ● ● ● ● — — — ● — — ● ● ● ● ● ● ● ● ● ● 17 2011-03 日本宫城海域 9.0 ● — ● ● ● ● ● ● — — ● ● ● ● ● ● ● ● ● ● ● ● ● 18 2023-02 土耳其南部 7.8 ● — ● ● ● ● ● ● — ● ● ● ● ● ● ● ● ● ● ● ● ● ● 合计 17 2 11 18 18 18 15 4 1 2 9 3 4 15 8 6 7 7 18 14 7 18 7 表 7 全球地震巨灾对产业链结构的影响统计结果
Table 7. Statistical results of the impact of some global earthquake catastrophes on structure of industrial chain
序号 地震日期/
(年-月-日)震中地名 震级M/级 生产链 供应链 销售链 代理链 管理链 1 1868-08-13 智利阿里卡近海 8½ ● ● ● — — 2 1908-12-28 意大利西西里岛 7¼ ● ● ● — — 3 1923-09-01 日本横滨 7.9 ● ● ● — — 4 1939-01-25 智利奇廉 7.7 ● ● ● — — 5 1960-02-29 摩洛哥阿加迪尔 6.8 ● ● — — — 6 1970-05-31 秘鲁卡斯马近海 7.9 ● ● ● — — 7 1972-12-23 尼加拉瓜马拉瓜 6.2 ● ● — — — 8 1976-02-04 危地马拉城 7.5 ● ● ● — ● 9 1976-07-27 中国河北唐山 7.8 ● ● ● — ● 10 1988-12-07 亚美尼亚斯皮塔克 6.8 ● ● ● — — 11 1990-06-20 伊朗加兹温 7.9 ● ● ● ● ● 12 1999-08-17 土耳其伊斯坦布尔 7.6 ● ● ● ● ● 13 2001-01-26 印度普杰 7.7 ● ● ● ● ● 14 2004-12-26 印尼锡默卢岛 9.1 ● ● ● ● ● 15 2008-05-12 中国四川汶川 8.0 ● ● ● ● ● 16 2010-01-12 海地太子港 7.3 ● ● ● ● ● 17 2011-03-11 日本宫城海域 9.0 ● ● ● ● ● 18 2023-02-06 土耳其南部 7.8 ● ● ● ● ● 合计 18 18 16 8 10 表 8 全球地震巨灾对生产链的影响统计结果
Table 8. Statistical results of the impact of some global earthquake catastrophes on the production chain
序号 地震日期/
(年-月-日)震中地名 震级M/级 开采捕捞 种植养殖 研发培育 制作加工 上市产品 1 1868-08-13 智利阿里卡近海 8½ ● ● — — — 2 1908-12-28 意大利西西里岛 7¼ ● ● — — — 3 1923-09-01 日本横滨 7.9 ● ● — ● ● 4 1939-01-25 智利奇廉 7.7 — ● — — ● 5 1960-02-29 摩洛哥阿加迪尔 6.8 — ● — — ● 6 1970-05-31 秘鲁卡斯马近海 7.9 ● ● — — ● 7 1972-12-23 尼加拉瓜马拉瓜 6.2 — ● — — ● 8 1976-02-04 危地马拉城 7.5 ● ● — — ● 9 1976-07-27 中国河北唐山 7.8 ● ● ● ● ● 10 1988-12-07 亚美尼亚斯皮塔克 6.8 — ● ● ● ● 11 1990-06-20 伊朗加兹温 7.9 — ● — ● ● 12 1999-08-17 土耳其伊斯坦布尔 7.6 ● ● ● ● ● 13 2001-01-26 印度普杰 7.7 — ● — ● ● 14 2004-12-26 印尼锡默卢岛 9.1 ● ● ● ● ● 15 2008-05-12 中国四川汶川 8.0 ● ● ● ● ● 16 2010-01-12 海地太子港 7.3 ● ● — ● ● 17 2011-03-11 日本宫城海域 9.0 ● ● ● ● ● 18 2023-02-06 土耳其南部 7.8 ● ● ● ● ● 合计 12 18 7 11 16 表 9 全球地震巨灾对供应链的影响统计结果
Table 9. Statistical results of the impact of some global earthquake catastrophes on the supply chain
序号 地震日期/
(年-月-日)震中地名 震级M/级 采购 运输 存储 配送 1 1868-08-13 智利阿里卡近海 8½ — ● — — 2 1908-12-28 意大利西西里岛 7¼ ● ● — — 3 1923-09-01 日本横滨 7.9 ● ● ● — 4 1939-01-25 智利奇廉 7.7 — ● ● — 5 1960-02-29 摩洛哥阿加迪尔 6.8 — ● ● — 6 1970-05-31 秘鲁卡斯马近海 7.9 ● ● — — 7 1972-12-23 尼加拉瓜马拉瓜 6.2 — ● — — 8 1976-02-04 危地马拉城 7.5 — ● ● — 9 1976-07-27 中国河北唐山 7.8 ● ● ● ● 10 1988-12-07 亚美尼亚斯皮塔克 6.8 — ● ● ● 11 1990-06-20 伊朗加兹温 7.9 — ● ● ● 12 1999-08-17 土耳其伊斯坦布尔 7.6 ● ● ● ● 13 2001-01-26 印度普杰 7.7 ● ● — ● 14 2004-12-26 印尼锡默卢岛 9.1 ● ● ● ● 15 2008-05-12 中国四川汶川 8.0 ● ● ● ● 16 2010-01-12 海地太子港 7.3 ● ● ● ● 17 2011-03-11 日本宫城海域 9.0 ● ● ● ● 18 2023-02-06 土耳其南部 7.8 ● ● ● ● 合计 11 18 13 10 表 10 全球地震巨灾对销售链的影响统计结果
Table 10. Statistical results of the impact of some global earthquake catastrophes on the sale chain
序号 地震日期/
(年-月-日)震中地名 震级M/级 总销 分销 批发 零售 1 1868-08-13 智利阿里卡近海 8½ — — — ● 2 1908-12-28 意大利西西里岛 7¼ — — — ● 3 1923-09-01 日本横滨 7.9 — — ● ● 4 1939-01-25 智利奇廉 7.7 — — ● ● 5 1960-02-29 摩洛哥阿加迪尔 6.8 — — — ● 6 1970-05-31 秘鲁卡斯马近海 7.9 — — ● ● 7 1972-12-23 尼加拉瓜马拉瓜 6.2 — — — ● 8 1976-02-04 危地马拉城 7.5 — ● ● ● 9 1976-07-27 中国河北唐山 7.8 ● ● ● ● 10 1988-12-07 亚美尼亚斯皮塔克 6.8 — ● ● ● 11 1990-06-20 伊朗加兹温 7.9 — ● ● ● 12 1999-08-17 土耳其伊斯坦布尔 7.6 ● ● ● ● 13 2001-01-26 印度普杰 7.7 — ● ● ● 14 2004-12-26 印尼锡默卢岛 9.1 ● ● ● ● 15 2008-05-12 中国四川汶川 8.0 ● ● ● ● 16 2010-01-12 海地太子港 7.3 ● ● ● ● 17 2011-03-11 日本宫城海域 9.0 ● ● ● ● 18 2023-02-06 土耳其南部 7.8 ● ● ● ● 合计 7 11 14 18 表 11 风险等级划分
Table 11. Risk level classification
高度风险 较高风险 中度风险 较低风险 低度风险 Ri>0.8 0.8≥Ri>0.6 0.6≥Ri>0.4 0.4≥Ri>0.2 0.2≥Ri 表 12 产业链地震灾害风险评估
Table 12. Earthquake disaster risk assessment of industrial chains
排
序受冲击的产业链I0 相关行业I1 相关行业I2 相关行业I3 风险
指数风险
评估名称 数量 影响因子 名称 影响因子 名称 影响因子 名称 影响因子 1 金融 18 1.000 金融 1.000 房地产 1.000 公共管理 1.000 1.000 高 2 教育 18 1.000 教育 1.000 建筑 1.000 水电气暖 0.778 0.944 高 3 农工贸 18 1.000 运输仓储 0.944 批发零售 1.000 水电气暖 0.778 0.919 高 4 蔬菜 18 1.000 农林牧渔 0.833 批发零售 1.000 运输仓储 0.944 0.911 高 5 肉类 18 1.000 农林牧渔 0.833 批发零售 1.000 运输仓储 0.944 0.911 高 6 农业 17 0.944 农林牧渔 0.833 运输仓储 0.944 水电气暖 0.778 0.805 高 7 医药 15 0.833 卫生社保 0.889 公共管理 1.000 批发零售 1.000 0.792 较高 8 文体娱乐 14 0.778 文体娱乐 0.778 住宿餐饮 1.000 水电气暖 0.778 0.657 较高 9 畜牧渔 11 0.611 农林牧渔 0.833 运输仓储 0.944 水电气暖 0.778 0.521 中 10 电信 15 0.833 信息服务 0.389 水电气暖 0.778 公共设施 0.667 0.479 中 11 旅游 7 0.389 住宿餐饮 1.000 建筑 1.000 公共管理 1.000 0.389 较低 12 制造 9 0.500 制造 0.500 建筑 1.000 运输仓储 0.944 0.381 较低 13 服装 8 0.444 制造 0.500 运输仓储 0.944 批发零售 1.000 0.337 较低 14 媒介 7 0.389 信息服务 0.389 房地产 1.000 公共管理 1.000 0.282 较低 15 Internet 7 0.389 信息服务 0.389 公共管理 1.000 水电气暖 0.778 0.260 较低 16 IT 7 0.389 信息服务 0.389 科技服务 0.500 水电气暖 0.778 0.202 较低 17 高新技术 6 0.333 科技服务 0.500 制造 0.500 水电气暖 0.778 0.190 低 18 化工 4 0.222 运输仓储 0.944 制造 0.500 批发零售 1.000 0.183 低 19 钢铁 4 0.222 采矿 0.167 运输仓储 0.944 批发零售 1.000 0.135 低 20 汽车 3 0.167 制造 0.500 运输仓储 0.944 批发零售 1.000 0.126 低 21 林业 2 0.111 农林牧渔 0.833 采矿 0.167 批发零售 1.000 0.075 低 22 煤炭 2 0.111 采矿 0.167 运输仓储 0.944 批发零售 1.000 0.068 低 23 造纸 1 0.056 制造 0.500 农林牧渔 0.833 批发零售 1.000 0.040 低 表 13 恢复等级划分
Table 13. Recovery level classification
强 较强 中 较弱 弱 CR>0.8 0.8≥CR>0.6 0.6≥CR>0.4 0.4≥CR>0.2 0.2≥CR 表 14 产业链震灾恢复能力评估
Table 14. Evaluation of the disaster recovery capability of the industrial chains
排序 产业链 自主恢复
比率S外部干预
比率H恢复时间/d 恢复指数 综合
指数恢复
能力应急t 正常T 应急恢复 正常恢复 1 Internet 0.9 0.1 2 15 1.000 1.000 1.000 强 2 电信 0.8 0.2 2 15 0.966 0.939 0.952 强 3 医药 0.3 0.7 2 15 0.793 0.636 0.715 较强 4 肉类 0.4 0.6 3 20 0.552 0.523 0.537 中 5 蔬菜 0.4 0.6 3 20 0.552 0.523 0.537 中 6 服装 0.3 0.7 3 20 0.529 0.477 0.503 中 7 IT 1.0 0.0 10 20 0.207 0.795 0.501 中 8 农业 0.9 0.1 10 20 0.200 0.750 0.475 中 9 畜牧渔 0.9 0.1 10 20 0.200 0.750 0.475 中 10 林业 1.0 0.0 10 30 0.207 0.530 0.369 较弱 11 造纸 1.0 0.0 10 30 0.207 0.530 0.369 较弱 12 汽车 1.0 0.0 10 30 0.207 0.530 0.369 较弱 13 高新技术 0.9 0.1 10 30 0.200 0.500 0.350 较弱 14 机械制造 0.9 0.1 10 30 0.200 0.500 0.350 较弱 15 钢铁 0.9 0.1 10 30 0.200 0.500 0.350 较弱 16 煤炭 1.0 0.0 15 30 0.138 0.530 0.334 较弱 17 化工 0.8 0.2 10 30 0.193 0.470 0.331 较弱 18 农工贸 0.9 0.1 15 30 0.133 0.500 0.317 较弱 19 媒介 0.8 0.2 15 30 0.129 0.470 0.299 较弱 20 文体娱乐 0.7 0.3 20 45 0.093 0.293 0.193 弱 21 金融 0.3 0.7 10 60 0.159 0.159 0.159 弱 22 教育 0.4 0.6 15 60 0.110 0.174 0.142 弱 23 旅游 0.7 0.3 30 90 0.062 0.146 0.104 弱 -
陈春慧,2013. 突发事件对汽车产业供应链影响研究−−以日本地震对广州汽车产业影响为例. 北方经贸,(3):60−61. 陈德富,杜义飞,李仕明等,2010. 汶川地震灾后重建中发展产业承载能力的机制与途径研究−−以德阳市磷化工产业为例. 电子科技大学学报(社科版),12(5):11−15,28.Chen D. F., Du Y. F., Li S. M., et al, 2010. Research on mechanism and path of developing industrial carrying capacity in post-disaster reconstruction-a case study of phosphorus chemical industry in Deyang City. Journal of University of Electronic Science and Technology of China (Social Sciences Edition), 12(5): 11−15,28.(in Chinese) 陈和,王荣祥,2011. 日本大地震对广东经济影响研究:基于产业链的视角. 战略决策研究,2(3):32−36.Chen H., Wang R.X. 2011. A research on the economic impact of Japan’s earthquake on Guangdong economy from the perspective of industry chain. Journal of Strategy and Decision-Making, 2(3): 32−36.(in Chinese) 杜兴,2015. 地震灾害间接经济损失评估方法研究. 西安:西安建筑科技大学.Du X.,2015. Research on assessment methods of earthquake-caused indirect economic loss. Xi’an:Xi’an University of Architecture and Technology. (in Chinese) 方美燕,吴海民,2008. 汶川震灾后龙门山经济圈的产业链重塑. 经济管理,(17):80−83.Fang M. Y., Wu H. M. 2008. A study on rebuilding and adjustments of industrial chain of post-disaster Longmenshan economic circle. Economic Management, (17): 80−83.(in Chinese) 刚健华,何青,2011. 日本大地震的溢出性经济影响与中国对策. 经济理论与经济管理,(6):69−76.Gang J. H., He Q. 2011. Economic consequences of the 2011 Tohoku Pacific earthquake. Economic Theory and Business Management, (6): 69−76.(in Chinese) 郭文博,2011. 浅析企业供应链全球化的弊端及对策−−基于日本地震对汽车产业影响的反思. 经济视角,(3):65,67. 郭增建,陈鑫连,1991. 城市地震对策. 北京:地震出版社. 黄汉权,2020-12-31(03). 多措并举增强产业链供应链自主可控能力. 中国经济时报. 景冰冰,施唯,2023. 2023年土耳其 M S7.8地震灾区损失与保险研究. 震灾防御技术,18(3):495−504.Jing B. B., Shi W. 2023. Study on loss and earthquake insurance of Turkey M S7.8 earthquake in 2023. Technology for Earthquake Disaster Prevention, 18(3): 495−504.(in Chinese) 刘贵富,赵英才,2006. 产业链的分类研究. 学术交流,(8):102−106.Liu G. F., Zhao Y. C. 2006. Exploring the category of industrial links. Academic Exchange, (8): 102−106.(in Chinese) 刘佳,2019. 震后汶川体育旅游产业发展研究. 成都:西华大学.Liu J.,2019. Research on Wenchuan sports tourism industry development after earthquake. Chengdu:Xihua University. (in Chinese) 潘泓旭,2011. 日本地震对我国IT行业的影响及解决措施. 现代商业,(14):59. 彭燕,2010. 地震灾后安县农业主导产业选择及发展研究. 雅安:四川农业大学.Peng Y.,2010. Research on development and selection of agricultural leading industries in Anxian after earthquake. Ya’an:Sichuan Agricultural University. (in Chinese) 沈铭辉,李淑娟,2011. 日本大地震对日本和世界经济的影响. 武汉金融,(5):40−42. 盛朝迅,2021. 新发展格局下推动产业链供应链安全稳定发展的思路与策略. 改革,(2):1−13.Sheng C.X. 2021. Thoughts and strategies for promoting the safe and stable development of industrial chain and supply chain under the new development pattern. Reform, (2): 1−13.(in Chinese) 史培军,2016. 灾害风险科学. 北京:北京师范大学出版社.Shi P. J.,2016. Disaster risk science. Beijing:Beijing Normal University Publishing Group. (in Chinese) 宋志平,张国民,刘杰等,2011. 全球地震灾害信息目录(9999B.C.~2010A.D.). 北京:地震出版社.Song Z. P., Zhang G. M., Liu J., et al, 2011. Disaster information catalog of global earthquakes (9999B.C.~2010A.D.). Beijing:Seismological Press. (in Chinese) 王彬彬,2012. 汶川地震灾后产业重建中的能力建设. 中国集体经济,(10):47−48. 王俊,陈柳钦,2011. 日本强震及核危机对全球经济的影响. 日本问题研究,25(2):1−10.Wang J., Chen L. Q. 2011. Impact of Japan’s earthquake and the nuclear crisis on the global economy. Japanese Research, 25(2): 1−10.(in Chinese) 魏来,陈宏,张洁,2011. 灾变性事件引起的消费品产业链价格波及效应研究. 管理评论,23(10):25−30.Wei L., Chen H., Zhang J. 2011. A study of price ripple effect of catastrophic events on the consumable industry chain. Management Review, 23(10): 25−30.(in Chinese) 徐玖平,杨春燕,2008. 四川汶川特大地震灾后重建的产业集群调整分析. 中国人口·资源与环境,18(6):142−151Xu J. P., Yang C. Y. 2008. Industrial cluster adjustment in post Wenchuan earthquake. China Population, Resources and Environment,, 18(6): 142−151. (in Chinese) 杨方,张义斌,葛旭波,2011. 东日本大地震及核危机对能源电力的影响. 中国电力,44(6):1−6.Yang F., Zhang Y. B., Ge X. B. 2011. Impact of Japan earthquake and nuclear crisis on energy and power industry. Electric Power, 44(6): 1−6.(in Chinese) 叶利扬,2013. 大地震后日本电子企业战略调整. 沈阳:辽宁大学.Ye L. Y.,2013. International strategic adjustment of Japanese electronics enterprises after the earthquake. Shenyang:Liaoning University. (in Chinese) 尹训飞,2019. 制造业产业链安全现状与对策. 中国工业和信息化,(7):54−59. doi: 10.3969/j.issn.1674-9138.2019.07.008 袁丽,邱桂华,2011. 日本大地震对中国及世界经济影响分析. 现代商贸工业,23(6):103−104. doi: 10.3969/j.issn.1672-3198.2011.06.061 张辉,申震宇,贾越等,2020. 供应链安全可控研究. 物流科技,43(11):143−145.Zhang H., Shen Z. Y., Jia Y., et al, 2020. Study on supply chain safety and controllability. Logistics Sci-Tech, 43(11): 143−145.(in Chinese) 张升,李兆焱,张思宇等,2023. 2023年土耳其7.8级地震灾害特征. 世界地震工程,39(3):45−55.Zhang S., Li Z. Y., Zhang S. Y., et al, 2023. Disaster characteristics of Turkey M 7.8 earthquake in 2023. World Earthquake Engineering, 39(3): 45−55.(in Chinese) 张义博,2021. 产业链安全内涵与评价体系. 中国经贸导刊,(10):55−59. 张玉来,2011. 从震灾看全球产业链上的“日本元素”−−以半导体产业为核心. 日本研究,(3):30−36. 张肇诚,陈琪福,郑大林,2013. 震例总结研究探讨. 北京:地震出版社. 赵福军,2019. 在构建开放经济格局中提升产业链安全. 国际商务财会,(4):3−6,11. -