-
在气候变暖的背景下,全球极端天气灾害发生的频率和造成的经济损失成倍增加,威胁着人们的生命财产安全,制约着社会经济的可持续发展。就灾害发生的频率、造成的经济损失和受影响的人口数而言,洪涝灾害是全球范围内最严重的自然灾害之一。在中国,过去20年的洪涝灾害影响了9亿人,占全世界受洪灾影响人口的55%[1]。由于社会经济的快速发展,中国洪涝灾害造成的直接经济损失呈现上升趋势[2-3]。因此,对中国洪涝灾害进行研究具有重要的现实意义和紧迫性。面对日益严重和频繁的气象灾害,除了通过减少温室气体排放来减缓气候变暖的速度,通过改进预报方法降低极端事件的不确定性外,如何采取和制定适应性战略和政策来减少灾害损失已成为一个重要的前沿课题[4]。国内外从2014年开始陆续颁布了一系列应对包括洪水灾害在内的灾害适应战略和政策[5-7],旨在制定和实施综合防灾减灾规划,建设具有风险抵御能力和可持续发展的韧性社会。为实现上述适应战略和政策,洪涝灾害综合经济损失评估是一项重要的基础性工作[8],可以为洪涝灾害的补偿和恢复机制提供数据支持。
洪涝灾害的经济影响包括直接经济损失和间接经济损失。直接经济损失包括农林牧渔业、工业信息交通运输业、水利设施等经济损失,由水利部公布。间接经济损失是指因直接损失导致的生产能力下降而造成的供需失衡通过产业链上下游关联传播引发的进一步的GDP下降[9-10]。间接经济损失评估是制定减灾政策的重要内容,可以为确定防洪基础设施的投入力度和灾后恢复重建定量模拟提供依据,也是适应性评估的重要组成部分,可以确保掌握适应措施的全部经济效益。目前,中国洪涝灾害间接经济损失评估的模型机制尚不完善,评估时长多以年、月为单位[11-12],评估区域多为单一的省份、市[13-16]而不是多区域整体评估,纳入的直接损失数据还不够全面。总的来说,洪涝灾害间接经济损失评估目前还没有统一可靠的方法。基于此,本文拟建立一套高精度、完善合理的中国洪涝灾害间接经济损失评估方法,在将洪涝灾害的直接影响转化为相应部门的生产能力损失的基础上,利用最新开发的基于复杂网络与适应性主体的环境经济系统CLUES模型,逐日高精度的模拟洪涝灾害导致的生产能力下降在产业链中传播引发的间接经济损失。以2017年中国洪涝灾害为例,展示模型方法的作用机理和评估效果,了解中国各省份各部门在面对洪涝灾害冲击时的脆弱性和关键风险节点,并对产业链进行风险预警。帮助制定适应性战略,减少洪涝灾害造成的间接经济损失。
-
洪涝灾害的直接冲击以及通过产业链传播的间接冲击造成的影响持续时间多在一年之内,因此需要使用中短期的高精度环境经济系统模型模拟。本文利用bob手机在线登陆新近研发的基于复杂网络和适应性主体的环境经济系统模型(Climate-Resilient and Low-Carbon Unfolding Economic Scenarios,以下简称CLUES),高时空精度地评估洪涝灾害的间接经济损失。模型采用面向对象的模拟方式,在虚拟世界[36]中设定生产、消费和运输主体以复杂而现实的方式进行互动,模型所刻画的主体行为均存在微观基础[37-39]。本文中模型的网络结构数据是中国省级多区域投入产出表(MRIO):“中国经济体”是研究的实例;各省份—部门是生产主体;各省份的最终消费是消费主体;生产主体之间、生产主体和消费主体之间的运输链条就是运输主体,其长度是按照省份间距离和陆路运输速度估算的运输天数。
将外生的洪涝灾害直接损失数据输入时,各主体在一个给定的复杂网络之中进行动态演化与互动,最终CLUES模型会输出所有地区—部门365天逐日的非均衡的GDP,与没有外生冲击时的稳态GDP对比,即可得出相应冲击通过产业链造成的间接经济损失。各主体的行为如下所述。
-
生产主体可以生产产品,向供应网络中相连接的其他主体送出产品与订单,如果遇到外部冲击,则可以表现出一定的适应性行为,如补充存货、调整上游供应商的订单份额、利用闲置生产能力、调整生产技术、灾后重建等。
1. 生产产品:基于列昂惕夫生产函数
在遭受洪灾影响的短期时间内,各生产主体不会发生投入要素的替代行为,因此生产量是基于列昂惕夫生产函数得出,受订单总量、生产能力、原材料供应三个方面限制
$$ {X^a}{{ = {\rm{min}}}}\left\{ {{O^{{{{\rm{tot}}}}}}{{, }}{X^{{{{\rm{cap}}}}}}{{, {\rm{min}}}}\left\{ {{X^{s'}}} \right\}{{, {\rm{min}}}}\left\{ {{X^{r',s'}}} \right\}} \right\} $$ (1) 其中,
$ {{X}}^{{a}} $ 为实际的生产量;${{O}}^{{{\rm{tot}}}}$ 代表总订单;${{X}}^{{}{{\rm{cap}}}}$ 为生产能力;r为地区;s为部门;$ {{X}}^{{s'}} $ 、$ {{X}}^{{r',s'}} $ 为不同产品库存约束下的可能生产水平。$$ {O^{{\rm{tot}}}} = \sum {{O^{ \leftarrow \left\langle {{p}} \right\rangle }}} (r',s') + \sum {{O^{ \leftarrow \left\langle {{c}} \right\rangle }}} (r') + {O^{{E}}} $$ (2) 其中,
${O^{ \leftarrow \left\langle {{p}} \right\rangle }}$ 为其他生产主体的订单;${O^{ \leftarrow \left\langle {{c}} \right\rangle }}$ 为来自消费主体的订单;${{O}}^{{←{{E}}}}$ 为出口订单。$$ {X^{{\rm{cap}}}} = \alpha (1 - \theta )\overline X $$ (3) 其中,
${{X}}^{{{\rm{cap}}}}$ 为生产能力;α为产能过剩比例;θ为相对于事件前水平的生产能力损失比例;$\overline X$ 为事前的生产能力。如果s′是同质产品
$$ {I^{{{{R}}}}}(s') = n(s')X_{t - 1}^aa(s') $$ (4) $$ {X^{s{'}}} = \left\{ {\begin{aligned} & {\frac{{I(s{'})}}{{a(s{'})}} \quad \quad \quad \quad \quad\quad \;{\rm{if }}\;\;I(s{'}) \geqslant {{\psi }}{I^R}(s{'})} \\ & {(\frac{{I(s{{'}})}}{{\psi {I^R}(s{'})}})\frac{{I(s{{'}})}}{{a(s{{'}})}} \;\quad\quad{\rm{if}} \;\;I(s{{'}}) < {{\psi }}{I^R}(s{{'}})} \end{aligned}} \right. $$ (5) 如果s′是异质产品
$$ {I^R}(r',s') = n(s')X_{t - 1}^aa(r',s') $$ (6) $$ {X^{r{'},s{'}}} = \left\{ {\begin{aligned} & {\frac{{I(r{'},s{'})}}{{a(r{'},s{'})}} \quad \quad\quad\quad\quad \quad \;\;{\rm{if }}\;\;I(r{'},s{'}) \geqslant \psi {I^R}(r{'},s{'})} \\ & {(\frac{{I(r{'},s{'})}}{{{{\psi }}{I^R}(r{'},s{'})}})\frac{{I(r{'},s{'})}}{{a(r{'},s{'})}} \quad\;\;{\rm{if }}\;\;I(r{'},s{'}) < \psi {I^R}(r{'},s{'})} \end{aligned}} \right. $$ (7) 其中,
${{I}}^{{{{R}}}}$ 为所需的库存水平;n为不同产品的目标使用天数;${{X}}^{{{{a}}}}$ 为实际的生产量;${a}$ 为单位产量的原材料投入需求(消耗);${{X}}^{{{\rm{s}}}}$ 为不同产品库存约束下的可能生产水平;I为当前的库存水平;$\psi $ 是所需库存水平的比率,低于该比率,生产主体将仅使用部分库存,以便在模拟步骤之间顺利生产(默认情况下,$\psi $ 的值为0)。2. 送出产品:基于习惯—适应性准则
产品的分配方式是基于习惯—适应准则。当实际生产量等于订单总量时,产品按照订单量进行分配;当实际生产量小于订单总量时,产品按照各类订单的份额进行分配。具体如下
$$ \begin{aligned} & \left\{ {\begin{aligned} & {{Z^{ \to \left\langle {{{{p}}}} \right\rangle }}(r',s') = {O^{ \leftarrow \left\langle {{{{p}}}} \right\rangle }}(r',s')} \\ & {{Z^{ \to \left\langle {{c}} \right\rangle }}(r')\begin{aligned} {} \end{aligned} = {O^{ \leftarrow \left\langle {{{{c}}}} \right\rangle }}(r')\begin{aligned} {}&{}&{}&\quad\quad\quad\quad\quad\;\;\;{{{{\rm{if}}}}\begin{aligned} \;{{X^{{a}}} = {O^{{{{\rm{tot}}}}}}}&{} \end{aligned}} \end{aligned}} \\ & {{Z^ \to }^{{{{E}}}} = {O^{{{{E}}}}}} \end{aligned} } \right. \\ & \left\{ {\begin{aligned} & {{Z^{ \to \left\langle {{{{p}}}} \right\rangle }}(r',s') = {X^{{a}}} \times \frac{{{O^{ \leftarrow \left\langle {{p}} \right\rangle }}(r',s')}}{{{O^{{\rm{tot}}}}}}} \\ & {{Z^{ \to \left\langle {{{{c}}}} \right\rangle }}(r')\begin{aligned} {} \end{aligned} = {X^{{a}}} \times \frac{{{O^{ \leftarrow \left\langle {{c}} \right\rangle }}(r')}}{{{O^{{{{\rm{tot}}}}}}}}\begin{aligned} {}&{}&\quad\quad\quad\;{{{{\rm{if}}}}\begin{aligned} \;{{X^a} < {O^{{\rm{tot}}}}}&{} \end{aligned}}&{} \end{aligned}} \\ & {{Z^ \to }^{{{{E}}}} \begin{aligned} {} \end{aligned} = {X^{{a}}} \times \frac{{{O^{{{{E}}}}}}}{{{O^{{{{\rm{tot}}}}}}}}} \end{aligned} } \right. \\ \end{aligned} $$ (8) 其中,
${Z^{ \to \left\langle {{{{p}}}} \right\rangle }}$ 为送到生产主体的产品数量;${Z^{ \to \left\langle {{{{c}}}} \right\rangle }}$ 为送到消费主体的产品数量;${Z^{ \to \left\langle {{{{E}}}} \right\rangle }}$ 为出口产品的数量。3. 送出订单:从而向目标水平补充存货
$$ {I_t}(s'){{ = }}{I_{t - 1}}(s') - X_{t - 1}^a{a_{t - 1}}(s') + \sum {{Z^{ \leftarrow \left\langle p \right\rangle }}} (r',s') $$ (9) 其中,
${{I}}_{{t}}\left({s'}\right)$ 为生产主体在t期关于$ {{s}}^{{'}} $ 的库存;${I}_{t-1}\left({s'}\right)$ 为生产主体在t−1期关于$ {{s}}^{{'}} $ 的库存;$ {{X}}_{{t-1}}^{{a}} $ 为t−1期实际的生产量;${\alpha }_{t-1}\left({s'}\right)$ 为t−1期的单位产量的$ {{s}}^{{'}} $ 消耗量;${{Z^{ \leftarrow \left\langle p \right\rangle }}}$ 为${r}^{{'}} $ 生产主体提供的$ {{s}}^{{'}} $ 产品量。如果s′是同质产品
$$ {I^{{{{T}}}}}(s') = n(s') \times {\rm{min}} \left\{ {{O^{{\rm{tot}}}},{X^{{\rm{cap}}}}} \right\}a(s') $$ (10) $$ {O^{ \to \left\langle {{{{p}}}} \right\rangle }}(s'){{ = }}a(s'){X^a} + \left\{{I^T}(s') - I(s')\right\}\frac{{\Delta t}}{{{\tau _I}}} $$ (11) 其中,
${{I}}^{{{{T}}}}\left({s'}\right)$ 为对于$ {s'} $ 的目标库存量;${O^{ \to \left\langle {{{{p}}}} \right\rangle }}$ 为向不同生产主体发送的订单;${\Delta }{t} $ 为每个时间步长的长度(取1天);${{\tau _I}}$ 为调整到目标库存水平需要的时间,取180天。4. 调整上游供应商:而缓解中间品短缺
根据如下方法将这期给不同区域供应商的订单份额进行更为现实的调整:如果某供应商上期提供的原材料大于给它的订单量,那么这期就增加给其的订单份额;反之,则减少。具体如下
$$ o_t^{ \to \left\langle p \right\rangle }(r',s') = o_{t - 1}^{ \to \left\langle p \right\rangle }(r',s') + \left(\frac{{{Z^{ \to \left\langle p \right\rangle }}(r,s|\left\langle p \right\rangle (r',s')}}{{\displaystyle\sum\limits_{r'} {{Z^{ \to \left\langle p \right\rangle }}(r,s|\left\langle p \right\rangle (r',s')} }} - \frac{{O_{t - 1}^{ \to \left\langle p \right\rangle }(r',s')}}{{O_{t - 1}^{ \to \left\langle p \right\rangle }(s')}}\right)\frac{{\Delta t}}{{{\tau _0}}} $$ (12) 其中,
$o^{ \to \left\langle p \right\rangle }$ 为给生产主体的订单份额;${Z^{ \to \left\langle p \right\rangle }}$ 为送到生产主体的产品量;${{\tau _0}}$ 为调整到目标订单需要的时间,取180天。向区域r′中的供应商订购
$$ O_{}^{ \to \left\langle p \right\rangle }(r',s') = {O^{ \to \left\langle p \right\rangle }}(s') \times {o^{ \to \left\langle p \right\rangle }}(r',s') $$ (13) 5. 利用闲置生产能力:从而保证产品供给
$$ {\alpha _{t + 1}} = \left\{ {\begin{array}{*{20}{l}} {{\alpha _t} + ({\alpha ^{\max }} - {\alpha _t})(\frac{{{O^{{{{\rm{tot}}}}}}{{ - }}{X^{{a}}}}}{{{O^{{{{\rm{tot}}}}}}}})\frac{{\Delta t}}{{{\tau _a}}} \;\;\;\;\;\;\quad {\rm{if}} \;{X^{{a}}}{{ < }}{O^{{{{\rm{tot}}}}}}} \\ {{\alpha _t} - ({\alpha _t} - 1)\frac{{\Delta t}}{{{\tau _a}}}\quad \quad \quad \quad \quad \quad \quad {\rm{if}} \;{X^{{a}}}{{ = }}{O^{{{{\rm{tot}}}}}}} \end{array}} \right. $$ (14) 其中,
${\alpha ^{\max }}$ 为最大可能的生产过剩比例,取1.2;${\tau _a}$ 为调整到最大生产能力需要的时间,取180天。6. 调整生产技术(从而适应中间品短缺)
对于开放经济体,假设s′是同质产品,如果满足下式,说明s′存在缺失
$$ {O_{t - 1}}^{ \to \left\langle {{p}} \right\rangle }(s') > {Z^{ \to \left\langle {{p}} \right\rangle }}(r,s|\left\langle p \right\rangle (s')) $$ (15) 在国内不能满足s′产品的供应时,生产主体将s′产品的订单转移到国外
$$ {a_t}(s'){{ = }}{a_{t - 1}}(s') - \left(\frac{{O_{t - 1}^{ \to \left\langle p \right\rangle }(s') - {Z^{ \to \left\langle p \right\rangle }}(r,s|\left\langle p \right\rangle (s')}}{{O_{t - 1}^{ \to \left\langle p \right\rangle }(s')}}\right){a_{t - 1}}(s') \times \frac{{\Delta t}}{{{\tau _A}^ \downarrow }} $$ (16) $$ {i_t} = {i_{t - 1}} + \sum\limits_{{{s' = 1}}}^s {{a_{{t}}}\left( {s'} \right)} $$ (17) 其中,
${{\tau _A} }$ 为技术适应需要的时间,取180天。当s′的订单和收货量满足下式,表明s′不存在缺失
$$ {O_{t - 1}}^{ \to \left\langle p \right\rangle }(s') = {Z^{ \to \left\langle p \right\rangle }}(r,s|\left\langle p \right\rangle (s')) $$ (18) 在国内能满足s′产品的供应时,生产主体应该将s′的订单从进口退回
$$ {a_t}(s'){{ = }}{a_{t - 1}}(s') + \left(\frac{{\bar a(s') - {a_{t - 1}}(s')}}{{\bar a(s')}}\right)(\bar a(s') - {a_{t - 1}}(s')) \times \frac{{\Delta t}}{{{{\tau _A}}^ \uparrow }} $$ (19) $$ {i_t} = {i_{t - 1}} + \sum\limits_{{{s' = 1}}}^s {{a_{{t}}}\left( {s'} \right)} $$ (20) 对于封闭经济,
$ {{a}}_{{t}}\left({s}{'}\right) $ 没有理由偏离$\overline a\left(s\right)$ 。7. 灾后重建(逐渐恢复生产能力)
当遭受洪灾的冲击导致生产能力降低时,生产主体会通过灾后重建慢慢恢复受损的生产能力,直到生产能力损失为0。具体如下
$$ {\theta _{t + 1}} = \left(1 - \frac{{\Delta t}}{{{\tau _\theta }}}\right){\theta _t} $$ (21) 其中,
${\tau _\theta }$ 为生产能力恢复需要的时间,取180天。8. 记忆:存储往期变量以指导行为
生产主体存储当前模拟期间的相关状态变量,这些变量将在未来的计算中使用。
-
消费主体也可以根据现实的供需状况,调整不同供应商所占的份额
1. 准备送出订单
根据如下方式将这期给不同地区生产者的订单份额进行更为现实的调整:如果某生产者上期提供的商品大于给它的订单量,那么这期就增加给其的订单份额;反之,则减少。如下所示
$$ {o}_{t}^{\to \langle p\rangle }({r}^{\prime },{s}^{\prime })={o}_{t-1}^{\to \langle p\rangle }({r}^{\prime },{s}^{\prime })+\left(\frac{{Z}^{\to \langle C\rangle }(r,s|\langle p\rangle ({r}^{\prime },{s}^{\prime })}{{\displaystyle \sum _{{r}^{\prime }}{Z}^{\to \langle C\rangle }(r,s|\langle p\rangle ({r}^{\prime },{s}^{\prime })}}-\frac{{O}_{t-1}^{\to \langle p\rangle }({r}^{\prime },{s}^{\prime })}{{O}_{t-1}^{\to \langle p\rangle }({s}^{\prime })}\right)\frac{\Delta t}{{\tau }_{0}} $$ (22) 其中,
${{Z}^{\to \langle C\rangle }}$ 为生产主体送到消费主体的产品量。向区域r′中的生产主体订购
$$ O_{}^{ \to \left\langle p \right\rangle }(r',s'){{ = }}{O^{ \to \left\langle p \right\rangle }}(s') \times{o^{ \to \left\langle p \right\rangle }}(r',s') $$ (23) 如果s是异质的,则无须计算份额。
2. 记忆:存储往期变量以指导行为
消费主体存储当前模拟期间的相关状态变量,这些变量将在未来的计算中使用。
-
运输主体模拟了生产主体之间、生产主体和消费主体之间实际的货物运输过程,是连接各主体的运输链条,其总长度是各主体之间的运输天数。各主体之间运输天数通过估算得到,涉及到的31个省及直辖市之间的距离来自高德地图,运输速度取35千米/小时[38]
1. 装载产品:把生产主体生产的产品装载到运输链条上
$$ {T_1} = \left\langle {{{{P}}}} \right\rangle {Z^{ \to \left\langle C \right\rangle }}(r') $$ (24) 2. 运输产品:每一天向前运输一步
$$ {T_{l + 1}} = {T_l} $$ (25) 3. 卸载产品:把运输链条上的产品卸载给相应的生产、消费主体
$$ \left\langle {P/C} \right\rangle {Z^{ \leftarrow \left( P \right)}}(r',s') = {T_l} $$ (26) -
本文选择2017年的全国洪涝灾害作为研究对象的三个原因:2017年全国洪涝灾害直接经济损失为2 142.53亿元,较接近历年来洪涝灾害直接经济损失的平均值;2017年中国发生洪涝灾害的区域分布较广,除天津外其他省份都受到了不同程度的直接影响;2017年中国洪涝灾害发生的流域具有代表性,长江中游区域,洞庭湖、松花江、珠江等流域等都是历年经常发生洪灾的地区。本文的研究对象为中国省级水平,进行多区域的分析,与单个省份的分析相比更符合现实,因为企业的产业链不会仅局限于本省,而是会遍布全国各个省份。
2017年中国31个省份(缺少港、澳、台地区数据)和42个部门间的多区域投入产出表[40]数据来自CEADs数据库;洪涝灾害直接损失数据来自2017年《水旱灾害公报》《全国洪涝灾情》《中国气象灾害年鉴》;生产能力(提供服务)损失比例估算数据来自《中国统计年鉴》《中国工业统计年鉴》《中国固定资产投资建设总额》。全国人均居住面积采用2016年住建部公布的数据。
-
1. 将全国范围的直接损失分配到各省份
《全国洪涝灾情》仅公布了全国范围的农牧渔业、工业交通运输业的经济损失,没有细分到各省,需要根据各省相应行业的实物损失占比将损失进行分配。
1)农林牧渔业损失分配方式
各省份农业受灾情况有相应的统计数据,而农业和林、牧、渔的受灾程度比较相同,根据各省份农业的损失占全国总损失的比例,将全国的农林牧渔经济损失分配到各省份。
$$ {\rm{Agr}}_i^{\rm{dis}} ={\rm {Agr}}_i^{\rm{affect}} \times 10\% + {\rm{Agr}}_i^{\rm{damage}} \times 30\% + {\rm{Agr}}_i^{\rm{fail}} \times 80\% $$ (27) $$ {\rm{loss}}_i^{\rm{AFAF}} = \dfrac{{{\rm{Agr}}_i^{\rm{dis}}}}{{\displaystyle\sum\limits_{i = 1}^{31} {{\rm{Agr}}_i^{\rm{dis}}} }} \times {\rm{loss}}_T^{\rm{AFAF}} $$ (28) 其中,
${{\rm{Agr}}}_{{i}}^{{\rm{dis}}}$ 为i省农业受损面积;${{\rm{Agr}}}_{{i}}^{{\rm{affect}}}$ 为i省农业受灾面积;${{\rm{Agr}}}_{{i}}^{{\rm{damage}}}$ 为i省农业成灾面积;${{\rm{Agr}}}_{{i}}^{{\rm{fail}}}$ 为i省农业绝收面积;${{\rm{loss}}}_{{i}}^{{\rm{AFAF}}}$ 为i省农林牧渔经济损失;${{\rm{loss}}}_{{T}}^{{\rm{AFAF}}}$ 为全国农林牧渔经济损失。2)工业交通运输业分配方式
工业交通运输业直接经济损失涵盖了采矿业、交通运输业、电力业和通信业的直接经济损失。首先以采矿业、交通运输业、电力业和通信业四个行业对气象灾害敏感度得分为权重,将工业交通运输业直接经济损失进行分配,采矿业权重为4.93%,交通运输业权重为51.41%,电力业权重为42.43%,通信业权重1.23%[41]。然后,将各个行业的全国总损失分配到各省份。分配方式同上述的农林牧渔业,依据相应行业各省份实物受损占比进行分配:采矿业以各省停产工矿企业个数占全国停产工矿企业总数占比为依据;电力行业以各省供电线路中断次数占全国供电线路中断总次数比例为依据;通信业以各省通信中断次数占全国通信中断总次数为依据。
将交通运输业进一步划分为铁路运输业、道路运输业、水上运输业和航空运输业(占比很低,忽略不计)四个子行业,分别计算各省份子行业的损失占比,汇总得到各省份交通运输业受损占比
$$ {\rm{Tra}}_i^{\rm{dis}} = \frac{{{\rm{Int}}_i^{\rm{rail}}}}{{{\rm{Int}}_T^{\rm{rail}}}} \times {\rm{PF}}_i^{\rm{rail}} + \frac{{{\rm{{Int}}}_i^{\rm{road}}}}{{{\rm{Int}}_T^{\rm{road}}}} \times {\rm{PF}}_i^{\rm{road}} + \frac{{{\rm{{Int}}}_i^{\rm{water}}}}{{{\rm{Int}}_T^{\rm{water}}}} \times {\rm{PF}}_i^{\rm{water}} $$ (29) $$ {\rm{loss}}_i^{\rm{Tra}} = \dfrac{{{\rm{Tra}}_i^{\rm{dis}}}}{{\displaystyle\sum\limits_{i = 1}^{31} {{\rm{Tra}}_i^{\rm{dis}}} }} \times {\rm{loss}}_T^{\rm{Tra}} $$ (30) 其中,
${{{\rm{Tra}}}}_{{i}}^{{\rm{dis}}}$ 为i省交通运输业受损占比;${{\rm{{Int}}}}_{{i}}^{{\rm{rail}}}$ ,${{\rm{{Int}}}}_{{i}}^{{\rm{road}}}$ ,${{{\rm{Int}}}}_{{i}}^{{\rm{water}}}$ 为i省铁路,公路中断次数,港口关停次数;${{\rm{Int}}}_{{T}}^{{\rm{rail}}}$ ,${{\rm{Int}}}_{{T}}^{{\rm{road}}}$ ,${{\rm{Int}}}_{{T}}^{{\rm{water}}}$ 为全国铁路、公路中断次数、港口关停次数;${{\rm{PF}}}_{{i}}^{{\rm{rail}}}$ 为i省铁路货运量占全国铁路货运量比例;${{\rm{PF}}}_{{i}}^{{\rm{road}}}$ 为i省公路货运量占全国公路货运量比例;${{\rm{PF}}}_{{i}}^{{\rm{water}}}$ 为i省水路货运量占全国水路货运量比例;${{\rm{loss}}}_{{i}}^{{\rm{Tra}}}$ 为i省交通运输业经济损失;${{\rm{loss}}}_{{T}}^{{\rm{AFAF}}}$ 为全国交通运输业经济损失。部分省份的应急管理厅可能会公布本省份当年的分行业详细的洪涝灾害损失情况,本文根据各省实际公布的数据对从全国分配到各省份的数据进行修正。
2. 将各省直接损失数值转换为MRIO表相应省份—部门生产(提供服务)能力损失比例
根据图1中洪涝灾害直接影响的类型,本文采用以下两种方法转化其直接影响:
1)直接经济损失转换为部门生产(提供服务)能力损失比例方式
产品产值类经济损失
$$ \theta _i^j = \frac{{{\rm{loss}}_i^j}}{{{\rm{value}}_i^j}} $$ (31) 资本存量类经济损失
$$ \theta _i^j = \frac{{{\rm{loss}}_i^j}}{{{\rm{capital}}_i^j}} $$ (32) 其中,
${\theta _i^j}$ 为i省j部门的生产能力(提供服务)损失比例;${{\rm{loss}}}_{{i}}^{{j}}$ 为i省j部门的直接经济损失;${{\rm{value}}}_{{i}}^{{j}}$ 为i省j部门的总产值;$ {{K}}_{{i}}^{{j}} $ 为i省j部门的固定资本存量。本文采用永续盘存法对上述部门的固定资本存量进行估计[42]
$$ K_{it}^j = K_{i(t - 1)}^j(1 - D) + I_{i(t - 1)}^j/{P_{it}} $$ (33) 其中,
$ {{K}}_{{it}}^{{j}} $ 为i省的j部门t期的固定资本存量;$ {{K}}_{{i(t-1)}}^{{j}} $ 为i省的j部门t−1期的固定资本存量;D为资本折旧率,取9.2%[43];$ {{}{I}}_{{i(t-1)}}^{{j}} $ 为i省的j部门t−1期建设投资;以2003年作为基期,以10%作为分母去除2003年各行业的固定资产建设投资数值;$ {{P}}_{{it}} $ 为i省t期的固定资产价格指数。2)房屋倒塌和人口受灾转换为部门生产(提供服务)能力损失比例方式
洪涝灾害受灾人口是指在洪涝灾害中生产生活遭受损失的人口数量(同一人遭受一次以上洪涝灾害时,只统计一次,不重复计)。受到洪灾影响的人可能在一段时间(几天或者一周)内,在上班工作的时候迟到几个小时,本文的损失数据是以每一日为步长输入的,因此将受灾人口处理为受到洪涝灾害的影响人口在当天全天不能工作。对于房屋倒塌,本文认为房屋所有者会对受损房屋进行处理修缮,房屋倒塌当天不能参加工作。因此把人口受灾比例和房屋受损比例之和当作一般制造业(MRIO中6~24部门)和建筑业的生产能力损失比例。
$$ {\rm{Pop}}_i^{\rm{dis}} = \frac{{{\rm{pop}}_i^{{\rm{affect}}}}}{{{\rm{pop}}_i^T}} $$ (34) $$ {\rm{Hou}}_i^{\rm{dis}} = \frac{{{\rm{Are{a}}^{{\rm{Ave}}}} \times {\rm{House}}_i^{\rm{Ave}} \times ({\rm{Hou}}_i^{\rm{damage}} \times 10\% + {\rm{Hou}}_i^{\rm{collapse}} \times 80\% )}}{{{\rm{Hou}}_i^T}} $$ (35) $$ \theta _i^{{\rm{man\& con}}} = {\rm{Pop}}_i^{\rm{dis}} + {\rm{Hou}}_i^{\rm{dis}} $$ (36) 其中,
${{\rm{Pop}}}_{{i}}^{{\rm{dis}}}$ 为i省人口受灾比例;${{\rm{pop}}}_{{i}}^{{\rm{affect}}}$ 为i省受灾人口;${{\rm{pop}}}_{{i}}^{{T}}$ 为i省年末总人口;${{\rm{Hou}}}_{{i}}^{{\rm{dis}}}$ 为i省房屋损坏比例;${{\rm{Area}}}^{{\rm{Ave}}}$ 为全国人均居住面积;${{\rm{House}}}_{{i}}^{{\rm{Ave}}}$ 为i省平均每户人数;${{\rm{Hou}}}_{{i}}^{{\rm{damage}}}$ 为i省房屋受损数;${{\rm{Hou}}}_{{i}}^{{\rm{collapse}}}$ 为i省房屋倒塌数;${{\rm{Hou}}}_{{i}}^{{T}}$ 为i省房屋建筑总面积;${\theta _i^{{\rm{man\& con}}}}$ 为i省一般制造业和建筑业部门的生产能力损失比例。3. 将各省份年度损失分配到洪涝灾害发生的具体日期
上述计算得到的各项数据为全年的总损失情况,是对每一天受到的洪涝灾害损失汇总得到的。CLUES模型要求以每一天为步长对各省份—部门进行生产(提供服务)能力损失比例进行输入,本文根据每一场洪涝灾害造成对应日期的人口受灾数量、直接经济损失以及应急响应等级、洪水预警等级、紧急状态将洪涝灾害的年度总损失分配到具体的日期。
4. 将各省份—部门生产(提供)服务能力损失比例纳入CLUES模型进行模拟
1)用中国省级多区域投入产出表为中国经济体赋值;
2)初始化世界和主体,根据已有变量推出其他模拟需要的变量;
3)洪灾冲击发生,输入各省份—部门当天的生产(提供服务)能力损失比例;
4)连接各主体,使其按上文给定顺序行动;
5)下一天开始,重复3)和4),总共模拟365天(从2017年5月1日—2018年4月30日)。
-
图2显示了利用CLUES模型模拟的2017年中国洪涝灾害直接损失通过产业链传播造成GDP的间接损失。横轴是模拟的日期,2017年5月1日是第一天,共365天。黑色虚线表示洪涝灾害冲击的时间范围,从5月9日开始到11月7日结束。将稳态的GDP值减去CLUES模型模拟的在2017年洪涝灾害冲击下的GDP值,得到相应的损失值。每日损失区域表示每日的GDP损失,可以清晰地再现2017洪涝灾害的过程。最大单日GDP损失26.5亿元出现在7月1日,主要是长江2017年第1号洪水导致湖南、广西、江西、四川等省遭受较大直接损失带来的冲击。其次是6月24日的16.7亿元,受江西、浙江、贵州等地洪涝灾害的冲击。其他单日GDP损失较大的日期包括8月23日,是由于台风天鸽在广东省登陆引起的洪涝灾害带来的冲击。以及10月5日,是由于发生在湖北的洪涝灾害带来的冲击。尽管7月15日吉林发生了严重的洪涝灾害,但该冲击并未造成较大的单日GDP损失。累积损失区域表示累计的GDP损失,2017年洪涝灾害通过产业链传播的间接经济损失是621.6亿元,占直接经济损失2142.5亿元的29%。
图3展示了各省每日GDP损失曲线,以及洪涝灾害发生时间和单次占全年比例(仅标注洪灾直接损失超过10亿元的省份)。图3中被标注的曲线都有明显的波峰,这些波峰出现的时间与各省份洪涝灾害直接冲击时间高度一致,而曲线其他位置相对平缓,这表明这些省份的间接损失主要是由该省份洪涝灾害直接冲击造成的,而其他省的洪涝灾害通过省际产业链传递造成的损失相对较小。这可能是由于运输成本等原因,企业对产业链上下游大部分进行了本地化。而图3中未标注省份(海南、青海和新疆除外)的曲线没有明显的波峰且曲线波幅较大,说明这些省份的间接损失主要是由于其他省份的洪涝灾害通过产业链传播造成的,而不是本省洪涝灾害导致的,因为这些省份受洪灾的影响较小,直接损失很低几乎为零。这表明,洪涝灾害的直接冲击会外溢到其他地区,可能是因为这些省份与受洪涝灾害冲击较大的省份之间产业链连接紧密。海南、青海和新疆等地洪涝灾害造成的直接经济损失虽然很小,但受灾人口比例较高,导致其每日GDP损失曲线出现明显峰值。
各省间接经济损失出现时间与该省份遭受直接冲击的时候非常一致,表明产业链效应是立即发生的,且在间接经济损失在出现后的一周之内就降至很低的水平,表明洪涝灾害的直接冲击对各省份经济体来说影响整体较小。各省份损失曲线的高点多集中在7月初,是由于2017年长江1号洪水造成的生产能力下降通过产业链传播造成,个别省份损失曲线的高点出现在该省份受洪涝灾害严重影响的时间。
将MRIO表中42部门按照行业类型汇总为了六大部门类,包括农林牧渔业(1)、采矿业(2~5)、制造业(6~23)、电力燃气水供应业(24~26)、建筑业(27)和服务业(28~42)。图4展示了全国和各省上述部门类的逐日和累积损失,各部门类的损失主要发生在6月下旬、7月上旬、8月下旬和10月初,是因为部分省份在这些时期遭受了较重的洪涝灾害冲击。制造业的损失占到所有部门损失的50.6%,因为人口受灾会导致制造业停产,而因为产业链的密切联系会遭受到进一步的级联损失。其次是服务业,占25.6%,服务业作为末端行业,各行业的直接冲击会通过产业链传导对其造成影响。而受直接损失冲击严重的农业、采矿业等间接经济损失占比很低,因为这些行业多处于生产链的初始端,不易受到产业链的级联影响。各省份各部门损失占比与全国的情况相似,且直接经济损失越大的省份,其制造业损失占比也相对较高;直接经济损失越小的省份,服务业损失占比越高。
图5横轴代表模拟的365天,纵轴是按照1302个省份—部门损失值的降序排序,颜色变化代表GDP损失,颜色轴上限为13万元。本文使用热力图显示了2017年中国洪涝灾害的经济足迹。从图5可以发现,自2017年长江1号洪水发生后,前200个省份—部门每天的GDP损失都超过13万元,直到模拟结束,这表明产业链对这些省份—部门的影响将持续很长一段时间。这200个省份—部门包括受洪涝灾害直接冲击较严重的湖南、江西、湖北、广东、四川等省的制造业部门,以及受洪涝灾害冲击较小的上海、北京等地区的金融业、山西的煤炭采选产品;第201~700个省份—部门,损失值在7—10月份较大,而在洪涝灾害直接冲击结束后,变得很小,表明产业链的影响在几个月之内逐渐退。而701~1 302个省份—部门,损失值非常低,且很快接近于零,表明产业链冲击的时间非常短。从时间轴上看,各省份—部门损失较多的日期集中在洪涝灾害冲击较严重的7月初、8月下旬及10月上旬。
图6展示了模拟第62天(2017年7月4日)每个地区每种产品的稀缺性指数,纵轴是按MRIO表中顺序排序的42个部门,每个部门代表一种产品,横轴是按间接经济损失占GDP比例降序进行排序的31个地区。稀缺性指数被定义为(产品需求—产品供应)/(产品需求)。产品供给是这一天运到这个地区的该产品,产品需求是这一天所有企业生产所需的这种产品加上消费者的需求。每个单元格的颜色代表该地区部门的产品稀缺指数。如果这个指标大于0,就是供不应求。本文计算了31个地区42种产品365天的稀缺指数,可以绘制365张与图6相同的图片,由于篇幅的限制,本文选择了一天内所有地区—部门的稀缺指数之和最大的进行展示。
图6显示,在7月4日这一天,湖南省的制造业稀缺指数都超过了0.1,表明超过10%的需求得不到满足,进而导致生产能力损失和附加值损失。广西省的金属和非金属开采、非金属制品和金属冶炼、专用、通信设备、金属制品修理,建筑业等稀缺指数都较大,接近0.1。因为湖南和广西两个省份在7月初遭受了当年“长江1号洪水”的冲击,造成的直接损失通过产业链进行的传播、引发了上述产品的短缺。江西省的采矿业、金属制品和电力行业稀缺指数较高,接近0.08,这些行业的稀缺是其GDP损失的主要原因。从各种产品上来看,各地区的专用设备稀缺指数都比较高,这说明该产品容易受到产业链的影响,进而导致间接经济损失。因此,适应性资金应该优先像该部门进行分配,比如适当增加其库存,提高产能过剩的比例等,来减少直接冲击通过产业链传播造成的间接经济损失。
通过经济足迹和稀缺指数图,本文逐日高精度展示了洪涝灾害通过产业链传播造成的经济损失情况以及损失产生的原因。本文的研究结果可以通过预测具体时间、具体地区的产品短缺情况,为供应链进行预警,从而尽可能实现洪灾通过产业链传播引发的间接经济损失的避免。
根据表1中2017年洪涝灾害冲击通过产业链传播造成的间接经济损失在各省的分布情况和大小, 31个省份分为四个部分:经济欠发达地区,如青海、宁夏、西藏,这些省份的间接损失都在2亿元以下;二是受洪涝灾后冲击较小的省份,如天津、北京、上海,或者甘肃、辽宁、黑龙江等受到影响但制造业不发达地区,这些省份的间接损失均在10亿元以下;三是受洪涝灾害影响较大的省份,如吉林、广西、陕西等,或者制造业发达的地区,如江苏、河南、重庆等,这些省份的间接经济损失在10~30亿元之间;四是受洪涝灾害冲击严重的湖南、广东、江西、湖北等地,其间接经济损失超过了40亿元,因为企业产业链上下游企业大部分还是在本省,受到洪涝灾害冲击严重的省份也是受间接经济影响最大的省份。
表 1各省份因洪涝灾害导致的经济损失和比例
单位:亿元 省份 直接经济损失 间接经济损失 间接经济损失/直接经济损失 总经济损失 湖南 524.42 127.84 0.24 652.26 吉林 368.92 21.44 0.06 390.36 广东 314.38 45.89 0.15 360.27 江西 106.74 64.60 0.61 171.34 广西 130.75 28.04 0.21 158.79 湖北 99.53 58.73 0.59 158.26 陕西 98.04 17.25 0.18 115.29 四川 66.83 16.84 0.25 83.67 辽宁 69.03 9.95 0.14 78.98 浙江 46.96 22.94 0.49 69.90 云南 52.22 14.56 0.28 66.78 贵州 40.31 17.93 0.44 58.24 甘肃 49.55 4.91 0.10 54.46 重庆 15.94 29.05 1.82 44.99 黑龙江 29.28 7.93 0.27 37.21 福建 23.02 9.67 0.42 32.69 山东 14.06 18.45 1.31 32.51 河南 10.17 19.68 1.94 29.85 安徽 11.68 17.79 1.52 29.47 内蒙古 17.29 9.33 0.54 26.62 江苏 4.54 14.35 3.16 18.89 西藏 16.58 1.34 0.08 17.92 河北 8.32 9.25 1.11 17.57 山西 5.76 8.02 1.39 13.78 新疆 7.01 4.37 0.62 11.38 海南 6.01 4.32 0.72 10.33 北京 1.24 5.26 4.24 6.50 上海 0.12 6.21 51.73 6.33 天津 0.00 4.06 / 4.06 青海 2.42 0.75 0.31 3.17 宁夏 1.41 0.84 0.59 2.25 根据表1展示的洪涝灾害造成的各省2017年总经济损失以及间接和直接经济损失的比值(各省直接经济损失来自《水旱灾害公报》,各省间接经济损失值由笔者估算得到,总经济损失是直接和间接经济损失之和)。本文用该比值衡量各省份受洪涝灾害产业链传播影响的严重程度。间接经济损失与直接经济损失比值最高的是上海,其次是北京、江苏、河南、重庆,因为这些省份受到洪涝灾害直接冲击较小,或者这些省份与其他省份的产业链联系密切;吉林、西藏、甘肃、辽宁这些省份的比值都很小,因为这些省份受到洪涝灾害直接冲击相对较大,也可能是因为这些省份制造业不发达,产业链较短;湖南、广东、广西等受洪涝灾害直接冲击较严重,其间接经济损失约占直接经济损失的20%。
本文的结果是针对2017年特定洪涝灾害冲击情景,主要是发生在长江中游区域性大洪水,洞庭湖水系的流域性大洪水,以及台风登陆导致的珠江大洪水。中国洪水还可能发生在黄河流域和长江下游,不同洪涝灾害情景造成的间接经济损失情况与本文的结果可能存在差异,需要进一步的针对性分析。
Dynamic Loss Assessment of Flood Disaster Industrial Chain in China
-
摘要:随着水文现象强度频率的增加和社会经济的持续发展,中国洪涝灾害造成的直接经济损失呈上升趋势,为制定防灾减灾、适应战略和灾后补偿恢复政策,全面评估洪涝灾害的经济损失是一项重要工作。基于复杂网络和适应性主体的环境经济系统模型,在多区域评估框架下,逐日高精度模拟了2017年中国长江中游区域性洪水、洞庭湖流域性洪水、松花江洪水和台风天鸽风暴潮洪水等造成的直接损失通过产业链传播引发的间接经济损失,展示了各省份—部门受产业链冲击和产品短缺情况,体现了模型的供应链预警能力。结果显示,2017年间接经济损失为621.6亿元,占直接经济损失的29%,80%以上来自没有直接经济损失的制造业和服务业;洪涝灾害直接影响会通过产业链扩散到未受洪灾影响的省份,但间接经济损失大多发生在省内,这可能与各企业本地化的产业链选择有关;制造业发达地区单位直接冲击造成的间接经济损失更高,可能与其资本密集度高以及产业间联系密切的特定经济结构有关。建议各级政府加大非工程措施适应性资金投入,增强社会公众和承载体抗灾能力;建立洪水保险和风险补偿机制,促进灾区尽快恢复重建;关键省份—行业要进行库存储备和产能冗余,减少洪涝灾害造成的直接经济损失和间接经济损失。Abstract:With the increase in the intensity and frequency of hydrological phenomena and the rapid development of social economy, direct economic losses caused by floods in China show an upward trend. In order to formulate disaster prevention and mitigation plans, adaptive strategies, and post-disaster compensation and recovery policies, it is important to comprehensively evaluate the economic losses caused by floods. Under the framework of multi-regional assessment, the environmental economic system model based on complex network and adaptive agents was used to simulate the daily indirect economic losses caused by regional floods in China transmitted through the industrial chain, such as floods in the middle reaches of the Yangtze River, floods in Dongting Lake Basin, floods in Songhua River, and the storm surge of Typhoon Tiange in 2017. The model shows the impact of the industrial chain and product shortage of all sectors in the regions concerned and reflects the early warning ability of the supply chain. The results show that the indirect economic loss in 2017 was 62.16 billion yuan, accounting for 29% of the direct economic loss, and more than 80% of the loss came from the manufacturing and service industries, which suffered no direct economic loss. The direct impact of flood disasters also spread to unaffected provinces through the supply chain, but most of the indirect economic losses mostly occurred in this local province, which may be related to the localized supply chain choice of enterprises. The indirect losses caused by work units’ direct impact in areas with developed manufacturing industries are higher, which may be related to their specific economic structures with high capital intensity and close inter-industry ties. Governments at all levels should increase their investment of adaptation funds in non-engineering measures to improve disaster resistance of the public and the carrier, set up flood insurance and risk compensation mechanisms to promote the recovery and reconstruction of disaster areas as soon as possible. Key regions and industries should carry out inventory reserve and capacity redundancy, to reduce the direct and indirect economic losses caused by floods.
-
表 1各省份因洪涝灾害导致的经济损失和比例
单位:亿元 省份 直接经济损失 间接经济损失 间接经济损失/直接经济损失 总经济损失 湖南 524.42 127.84 0.24 652.26 吉林 368.92 21.44 0.06 390.36 广东 314.38 45.89 0.15 360.27 江西 106.74 64.60 0.61 171.34 广西 130.75 28.04 0.21 158.79 湖北 99.53 58.73 0.59 158.26 陕西 98.04 17.25 0.18 115.29 四川 66.83 16.84 0.25 83.67 辽宁 69.03 9.95 0.14 78.98 浙江 46.96 22.94 0.49 69.90 云南 52.22 14.56 0.28 66.78 贵州 40.31 17.93 0.44 58.24 甘肃 49.55 4.91 0.10 54.46 重庆 15.94 29.05 1.82 44.99 黑龙江 29.28 7.93 0.27 37.21 福建 23.02 9.67 0.42 32.69 山东 14.06 18.45 1.31 32.51 河南 10.17 19.68 1.94 29.85 安徽 11.68 17.79 1.52 29.47 内蒙古 17.29 9.33 0.54 26.62 江苏 4.54 14.35 3.16 18.89 西藏 16.58 1.34 0.08 17.92 河北 8.32 9.25 1.11 17.57 山西 5.76 8.02 1.39 13.78 新疆 7.01 4.37 0.62 11.38 海南 6.01 4.32 0.72 10.33 北京 1.24 5.26 4.24 6.50 上海 0.12 6.21 51.73 6.33 天津 0.00 4.06 / 4.06 青海 2.42 0.75 0.31 3.17 宁夏 1.41 0.84 0.59 2.25 -
[1] United Nations Office for Disaster Risk Reduction, Centre for Research on the Epidemiology of Disasters. The human cost of disasters: an overview of the last 20 years (2000—2019) [EB/OL]. (2020-10-12)[2022-02-06]. https://www.undrr.org/publication/human-cost-disasters-overview-last-20-years-2000-2019. [2] 中国水旱灾害防御公报编写组. 《中国水旱灾害防御公报2020》概要 [J]. 中国防汛抗旱, 2021, 31(11): 26-32.doi:10.16867/j.issn.1673-9264.2021323 [3] WILLNER S N, OTTO C, LEVERMANN A. Global economic response to river floods [J]. Nature Climate Change, 2018, 8(7): 594-598.doi:10.1038/s41558-018-0173-2 [4] 李宁, 张正涛, 陈曦, 等. 论自然灾害经济损失评估研究的重要性 [J]. 地理科学进展, 2017, 36(2): 256-63. [5] 中共中央, 国务院. 关于推进防灾减灾救灾体制机制改革的意见 [EB/OL ]. (2017-01-10) [2022-02-06]. http://www.gov.cn/zhengce/2017-01/10/content_5158595.htm. [6] 国家发展改革委, 住房城乡建设部. 关于印发气候适应型城市建设试点工作的通知 [EB/OL ]. (2017-02-21)[2022-02-06]. http://www.gov.cn/xinwen/2017-02/25/content_5170863.htm#allContent. [7] United Nations International Strategy for Disaster Reduction . Sendai framework for disaster risk reduction 2015—2030 [EB/OL]. (2015-03-18)[2022-01-20].https://www.undrr.org/publication/sendai-framework-disaster-risk-reduction-2015-2030. [8] 王宝华, 付强, 谢永刚, 等. 国内外洪水灾害经济损失评估方法综述 [J]. 灾害学, 2007(3): 95-99.doi:10.3969/j.issn.1000-811X.2007.03.022 [9] 吴先华, 徐中兵, 袁迎蕾, 等. 台风灾害的关联经济损失评估: 以江苏省为例 [J]. 灾害学, 2014, 29(2): 77-83. [10] 李宁, 吴吉东, 崔维佳. 基于ARIO模型的汶川地震灾后恢复重建期模拟 [J]. 自然灾害学报, 2012, 21(2): 68-75. [11] 李廉水, 蔡洋, 谭玲. 基于动态CGE模型的中国暴雨洪涝灾害综合经济损失评估研究 [J]. 河海大学学报(哲学社会科学版), 2020, 22(1): 28-36+106. [12] 张鹏, 李宁. ARIO模型对灾害引起的部门经济损失数据完备程度的敏感性分析: 以湖南为例 [J]. 灾害学, 2014, 29(04): 105-10.doi:10.3969/j.issn.1000-811X.2014.04.020 [13] 宁思雨, 黄晶, 汪志强, 等. 基于投入产出法的洪涝灾害间接经济损失评估: 以湖北省为例 [J]. 地理科学进展, 2020, 39(3): 420-432.doi:10.18306/dlkxjz.2020.03.007 [14] 张鹏, 李宁, 吴吉东, 等. 基于投入产出模型的区域洪涝灾害间接经济损失评估 [J]. 长江流域资源与环境, 2012, 21(6): 773-779. [15] 张正涛, 李宁, 冯介玲, 等. 从重建资金与效率角度定量评估灾后经济恢复力的变化——以武汉市“武汉市“市“率角度定量暴雨洪涝灾害为例 [J]. 灾害学, 2018, 33(4): 211-216.doi:10.3969/j.issn.1000-811X.2018.04.036 [16] 张鹏, 李宁, 刘雪琴, 等. 基于投入产出模型的洪涝灾害间接经济损失定量分析 [J]. 北京师范大学学报(自然科学版), 2012, 48(4): 425-431. [17] JONKMAN S N, BOČKARJOVA M, KOK M, et al. Integrated hydrodynamic and economic modelling of flood damage in the Netherlands [J]. Ecological Economics, 2008, 66(1): 77-90.doi:10.1016/j.ecolecon.2007.12.022 [18] VEEN A V D, LOGTMEIJER C. Economic hotspots: visualizing vulnerability to flooding [J]. Natural hazards, 2005, 36(1): 65-80. [19] 路琮, 魏一鸣, 范英, 等. 灾害对国民经济影响的定量分析模型及其应用 [J]. 自然灾害学报, 2002(3): 15-20.doi:10.3969/j.issn.1004-4574.2002.03.003 [20] OKUYAMA Y. Critical review of methodologies on disaster impact estimation [R]. UN Assessment on the Economics of Disaster Risk Reduction, 2008. [21] BORGOMEO E, VADHEIM B, WOLDEYES F , et al. The distributional and multi-sectoral impacts of rainfall shocks: evidence from computable general equilibrium modelling for the Awash Basin, Ethiopia [J]. Ecological Economics, 2018, 146: 621-632.doi:10.1016/j.ecolecon.2017.11.038 [22] JOSHI S R, VIELLE M, BABONNEAU F, et al. Physical and economic consequences of sea-level rise: a coupled GIS and CGE analysis under uncertainties [J]. Environmental and Resource Economics, 2015, 65(4): 813-839. [23] PAUW K, THURLOW J, BACHU M, et al. The economic costs of extreme weather events: a hydrometeorological CGE analysis for Malawi [J]. Environment and Development Economics, 2011, 16(2): 177-198.doi:10.1017/S1355770X10000471 [24] GREENBERG M R, LAHR M, MANTELL N. Understanding the economic costs and benefits of catastrophes and their aftermath: a review and suggestions for the US federal government [J]. Risk Analysis: An International Journal, 2007, 27(1): 83-96.doi:10.1111/j.1539-6924.2006.00861.x [25] HALLEGATTE S, RANGER N, MESTRE O, et al. Assessing climate change impacts, sea level rise and storm surge risk in port cities: a case study on copenhagen [J]. Climatic change, 2011, 104(1): 113-37.doi:10.1007/s10584-010-9978-3 [26] LI J, CRAWFORD-BROWN D, SYDDALL M, et al. Modeling imbalanced economic recovery following a natural disaster using input-output analysis [J]. Risk Anal, 2013, 33(10): 1908-23. [27] KOKS E E, BOCKARJOVA M, DE MOEL H, et al. Integrated direct and indirect flood risk modeling: development and sensitivity analysis [J]. Risk Anal, 2015, 35(5): 882-900.doi:10.1111/risa.12300 [28] MENDOZA-TINOCO D, HU Y, ZENG Z, et al. Flood footprint Assessment: A Multiregional Case of 2009 Central European Floods [J]. Risk Anal, 2020, 40(8): 1612-31.doi:10.1111/risa.13497 [29] ZENG Z, GUAN D, STEENGE A E, et al. Flood footprint assessment: a new approach for flood-induced indirect economic impact measurement and post-flood recovery [J/OL]. Journal of Hydrology, 2019, 579: 124204. https://www.sciencedirect.com/science/article/abs/pii/S0022169419309394. [30] MENDOZA-TINOCO D, GUAN D, ZENG Z, et al. Flood footprint of the 2007 floods in the UK: the case of the yorkshire and the humber region [J]. Journal of Cleaner Production, 2017, 168: 655-67.doi:10.1016/j.jclepro.2017.09.016 [31] VERSCHUUR J, KOKS E E, HALL J W. Observed impacts of the COVID-19 pandemic on global trade [J]. Nature Human Behaviour, 2021, 5(3): 305-7.doi:10.1038/s41562-021-01060-5 [32] TEZUKA S, TAKIGUCHI H, KAZAMA S, et al. Estimation of the effects of climate change on flood-triggered economic losses in Japan [J]. International Journal of Disaster Risk Reduction, 2014, 9: 58-67.doi:10.1016/j.ijdrr.2014.03.004 [33] WARD P J, JONGMAN B, WEILAND F S, et al. Assessing flood risk at the global scale: model setup, results, and sensitivity [J/OL]. Environmental Research Letters, 2013, 8(4): 044019. https://iopscience.iop.org/article/10.1088/1748-9326/8/4/044019/meta. [34] 张正涛, 崔鹏, 李宁, 等. 武汉市"2016.07. 06"暴雨洪涝灾害跨区域经济波及效应评估研究 [J]. 气候变化研究进展, 2020, 16(4): 433-41. [35] 曹玮. 洪涝灾害的经济影响与防灾减灾能力评估研究 [D]. 长沙: 湖南大学, 2013. [36] FARMER J D, FOLEY D. The economy needs agent-based modelling [J]. Nature, 2009, 460(7256): 685-686.doi:10.1038/460685a [37] INOUE H, TODO Y. Firm-level propagation of shocks through supply-chain networks [J]. Nature Sustainability, 2019, 2(9): 841-847.doi:10.1038/s41893-019-0351-x [38] OTTO C, WILLNER S N, WENZ L, et al. Modeling loss-propagation in the global supply network: the dynamic agent-based model acclimate [J]. Journal of Economic Dynamics and Control, 2017, 83: 232-269.doi:10.1016/j.jedc.2017.08.001 [39] HALLEGATTE S. An adaptive regional input‐output model and its application to the assessment of the economic cost of Katrina [J]. Risk Analysis: An International Journal, 2008, 28(3): 779-799.doi:10.1111/j.1539-6924.2008.01046.x [40] ZHENG H, ZHANG Z, WEI W, et al. Regional determinants of China’s consumption-based emissions in the economic transition [J/OL]. Environmental Research Letters, 2020, 15(7): 074001.https://iopscience.iop.org/article/10.1088/1748-9326/ab794f/meta. [41] 张钛仁, 宋善允, 田翠英, 等. 中国行业气象服务效益评估方法与分析研究 [J]. 气象软科学, 2007(4): 11. [42] 李宾. 中国资本存量估算的比较分析 [J]. 数量经济技术经济研究, 2011, 28(12): 21-36+54. [43] 金戈. 中国基础设施资本存量估算 [J]. 经济研究, 2012, 47(4): 4-14+100.