一种风暴潮灾害风险精细化评估的计算方法

文档序号:10489538阅读:224来源:国知局
一种风暴潮灾害风险精细化评估的计算方法【专利摘要】本发明涉及一种风暴潮灾害风险精细化评估的计算方法,包括对受风暴潮灾害影响的地区建立不同的分类网格单元;以海岸线为基准计算不同的网格单元质心到海岸线的距离;利用高精度的数字高程地形,计算每个网格内的平均高程;根据经验统计公式计算风暴潮最大增水,根据平均高程和最大增水,获取每个格网单元的平均深度;根据研究区内不同岸线类型分别赋不同的易损性,计算不同网格单元的最大可能性损失。本发明降低了目前风暴潮灾害风险评估的局限性,可以在更细的空间单元上反映损失可能性差异,提高评估结果的精度。【专利说明】一种风暴潮灾害风险精细化评估的计算方法
技术领域
[0001]本发明涉及沿海灾害风险评估的
技术领域
,特别涉及针对受风暴潮灾害影响的小尺度区域的损失率计算方法。【
背景技术
】[0002]风暴潮灾害不仅居我国海洋灾害之首,而且已成为威胁我国沿海经济发展的最严重的自然灾害之一。目如国内的风暴潮灾害综合风险评估技术体系还不成熟,风暴潮灾害风险评估结果不足以满足实际应用需求。[0003]风暴潮灾害发生的可能损失一直是灾害风险评估的难题。国外研究机构商业模型保险公司等针对风暴潮风险评估及损失综合评估开发了一系列模型,用于保险费率的厘定以及灾后理赔定损。美国联邦应急管理署开发的多灾种损失评估软件MH-HAZUS模型在洪水模块中考虑了风暴潮引发沿岸淹没的影响,可以评估风暴潮及其带来的洪水产生的直接破坏、间接损失评估及次生灾害影响评估,为灾后的风暴潮应对提供科学的决策支持。EQECAT巨灾风险模型公司根据亨伯河口和泰晤士河河口等欧洲地区实际情况,采用经验统计方法,开发了风暴潮灾情评估模型(石先武,谭骏,国志兴,等.风暴潮灾害风险评估研究综述.地球科学进展.2013,28(8):866-874.)。[0004]随着灾害保险特别是巨灾保险和巨灾再保险市场的发展,传统的风暴潮灾害损失经验统计模型满足不了应用的需求,具体表现在风暴潮风险评估的系统性研究还不够全面和深入,风暴潮灾害风险评估的研究成果还不能满足我国风暴潮风险管理的需求。[0005]在评价单元上,从全球和全国尺度上灾害损失过程主要以市、县行政区作为评价单元,也有的以土地利用或者城市用地单元开展区域综合自然灾害风险评估。(张斌,赵前胜,姜瑜君.区域承灾体脆弱性指标体系与精细量化模型研究.灾害学.2010,25(2):36-40;辜智慧,徐伟,葛怡,等.基于城市用地单元的自然灾害风险评估概念模型.中国安全科学学报·2012,ν·22(04):110-115·)。[0006]在评价计算方法上,对各类建(构)筑物、基础设施结构采取分类清单易损性分析。但由于需要收集的资料类型和数量比较多,工作量大,并且易损性存在地域差别,难以进行有效比较。[0007]风暴潮灾害风险评估结果向定量化、区域综合化、管理空间化的方向发展,现有的风暴潮灾害风险评估不能定量表现各区域之间的具体差别,因而不能满足风暴潮灾害风险管理的要求。【
发明内容】[0008]针对上述技术不足,本发明的目的提供一种风暴潮灾害风险精细化评估方法。[0009]本发明解决其技术问题所采用的技术方案是:一种风暴潮灾害风险精细化评估的计算方法,包括以下步骤:[0010]1)将待评估地区的风暴潮影响区域矢量图层进行网格单元分割;[0011]2)根据风向、风速得到该区域的风暴潮暴增水指数f1;[0012]3)将各网格单元的平均高程作为各网格单元的风暴潮暴露度指数D;平均高程低于心的网格单元执行下一步骤;[0013]4)根据海岸线类型得到风暴潮浪对岸线的易损性指数f3;[0014]5)根据各网格单元距海岸线的最近距离得到风暴潮距离衰减指数f2;[0015]6)根据风暴潮暴增水指数fi、风暴潮暴露度指数D、风暴潮距离衰减指数f2和易损性指数f3得到各单元的损失率;损失率值用于表示风暴潮灾害的风险。[0016]所述根据风向、风速得到该区域的风暴潮暴增水指数通过以下公式得到:[0017]fi(F,t,9)=AFtsinB[0018]其中,fi为风暴潮暴增水指数,λ为系数,F为风速,Θ为风向与岸线的角度,t为大风持续时间。[0019]所述根据各网格单元距海岸线的最近距离得到风暴潮距离衰减指数5通过以下公式得到:[0020][0021]其中,d为某网格单元质心坐标到矢量图层中海岸线的最短距离,k为调节系数。[0022]所述根据海岸线类型得到风暴潮浪对岸线的易损性指数f3通过以下公式得到:<0.8-1.0,z为建设围堤0.6-0,8,z为养殖围堤[0023]f3(z)t-0:.4-0.6,z为码1-汴线0.2-0.4,z为砂砾岸线<0-0.2,z为其他岸线[0024]所述损失率1^=(;1^1(?,1:,9)-〇)*1^((1)*1^3(2)。[0025]本发明具有以下有益效果及优点:[0026]1.本发明主要针对沿海受风暴潮严重影响的地区,考虑了暴风风向、暴风风速大小、暴风持续时间等因素,同时考虑了地形、岸线类型、海岸线走向等诸多其他因素(如图1海水淹没示意图所示,表示风暴潮损失与海水淹没损失与高程、距岸线距离的关系。d为某单元网格A中心到海岸线的距离,D为某单元网格A内的平均高程,Tmax为风暴潮增水的最高增水水位),提出了一套新的风暴潮灾害损失的计算方法,可以在沿海单元不同地形高低、不同岸线类型,不同岸线走向情况下,该方法可以计算出风暴潮影响的不同区域范围的可能性损失差异,也可以在不同的风速和风向下,为沿海区域内到海岸线不同距离的网格单元计算出不同的损失率。[0027]2.针对现有风险评估技术在较小空间尺度并未实现的研究现状,本发明的"精细化评估"的含义,即体现了比国内外目前现有的风暴潮灾害风险在评估技术方法上可行,在空间尺度上分辨率更高。[0028]3.本发明综合考虑风、地形、岸线类型等多种因素,其优点在于实现了风暴潮灾害风险的微观尺度的空间模拟,在空间上较好地反映了风暴潮损害的可能性损失差异,精度较高。为沿海风暴潮灾害频发的地区灾害保险的实施提供了基础,同时为灾害保险费率的计算提供了依据,也为今后该地区的风暴潮灾害风险区域规划等提供决策依据。【附图说明】[0029]图1是风暴潮灾害损失与影响因素示意图;[0030]图2是本发明的方法流程图;[0031]图3是沿海岸线生成受风暴潮影响区域的单元格网(IOOmX100m)图;[0032]图4a是西北风向、东北岸线走向之间夹角的示意图;[0033]图4b是东北风向、东北岸线走向之间夹角的不意图;[0034]图5是在不同风速下的风暴潮淹没范围与可能性损失风险程度(IOOmX100m)。【具体实施方式】[0035]下面结合实施例对本发明做进一步的详细说明。[0036]风暴潮是我国海洋灾害造成沿海地区损失最严重的一种灾害。本发明以自然灾害风险评估为基础,提出了微观尺度的风暴潮灾害风险评估方法,利用历史灾情和孕灾环境中致灾因子与灾害发生的关系大小,并以龙口市沿岸风暴潮灾害风险为例进行了定量化评估研究。本发明以地形、天气等环境参数计算风暴潮灾害可能损失,计算灾害损失率,确定承灾体灾害风险的可能损失率,从而为风暴潮灾害风险评估提供技术支撑,为沿海灾害财产保险的基准费率制定提供了计算依据,同时也为其他洪水淹没损失的计算提供了思路借签。[0037]如图2所示,本发明提供一种针对受风暴潮灾害影响的小尺度区域的损失率计算方法。其中,首先对受风暴潮灾害影响的地区建立不同的分类网格单元,以此网格单元作为受风暴潮影响差异的分类基础,其次,以海岸线为基准计算不同的网格单元质心到海岸线的距离,再次,利用高精度的数字高程地形,计算每个网格内的平均高程;接着,根据经验统计公式计算风暴潮最大增水,根据平均高程和最大增水,获取每个格网单元的平均深度;然后,根据研究区内不同岸线类型分别赋不同的易损性,最后在以上基础上,根据建立的模型计算不同网格单元的最大可能性损失。这种方法降低了目前风暴潮灾害风险评估的局限性,可以在更细的空间单元上反映风暴潮灾害损失可能性差异,提高了评估结果的精度。[0038]①基础地理信息入库和处理。利用地理信息系统软件(arcgis10.0)软件将评估区域的地形图,海岸线矢量图层统一到平面直角坐标系下,导入数据库;根据实地调查或历史资料,沿海岸线确定受风暴潮影响的地区,生成淹没面积的polygon(多边形)图层;对受风暴潮灾害影响的区域进行单元分割,采用fishnet(渔网)命令以100mXIOOm网格单元或者以所需网格单元进行分割。如图3所示单元格网与海岸线在空间上的相对位置,每个单元具有中心经炜度坐标等信息,经②-⑦步骤处理后,每个网格单元属性表中包括对应岸段的最大增水、单元格平均高程、单元中心距海岸线的距离、淹没水深、易损性等信息。以下步骤对每个网格单元进行处理。[0039]②计算暴风增水指数h。在海水比较浅的莱州湾地区,风的动力作用远大于气压差的作用,对增水量值的估计,主要考虑风向α、风速F和风的持续时间T。风暴潮各阶段的增水量的评估采用以下经验公式进行:[0040]Tmax=AFtsinB(1)[0041]Tmax为风暴潮增水指数,即为复相关系数,F为风速,Θ为风向与岸线的角度,t为大风持续时间。λ根据验潮站历史水位与风速的关系建立统计模型计算可得,具体方法参见文献:[0042][1]山义昌,张芹,杨付津,等.莱州湾风暴潮的形成与增水量评估.山东气象.2006(04):4-5。[0043][2]尹延英,徐振山,张胜平,等.莱州湾夏营站温带风暴潮预报方案编制[J].海洋预报.1996(02):58-62。[0044][3]徐振山.莱州湾温带风暴潮的分析及其预报.山东水利科技.1995(2):6_8。[0045][4]张胜平,徐振山,王军,等.莱州湾温带风暴潮预报研究.海洋预报.2002,64-72(5)。[0046]根据上述文献得计算公式[0047]式中,fmax是该段岸线最大增水的估计值,fmar是该段岸线最大增水的均值。[0048]在本实例中,λ取值如下式:[0049]在东南风阶段,Tmaxl=0.005Ftsin9(3)[0050]在东北风阶段,Tmax2=0.0lFtsinQ(4)[0051]其中,Tmaxl、Tmax2为增水高度值(m);F为风速(m/s);t为大风持续时间(h);θ为风向与岸线之间的角度。当Θ=90°时,sin0=1,大风的增水作用最大;当大风风向偏离90°时,sin0〈l,当风向转到大于180°时,大风没有增水作用为0,即风暴方向由海向岸有增水作用,反之由岸向海方向起减水作用。[0052]风暴潮增水量的大小,不仅和诱发的天气系统的路径、强度、风速等因素有关,还和风向与岸线走向之间的角度Θ有关。因此,风暴潮最高水位Tmax不仅受到风速的影响,而且受风向α和岸线地形走势邱勺影响,最大夹角Θ计算可用以下公式计算。[0053]θ=α-β,〇。<α<360。,0。<β<360。(5)[0054]图4a和图4b是风向、东北岸线走向之间夹角的示意图,本发明以岸线β不变为例说明风向α变化时Θ的计算方式。以正北方向N为0轴起,α为风向与0轴之间的角度,β为海岸线与〇轴之间的角度,Θ计算公式适用于任意风向与岸线走向的组合。[0055]因此,利用上述经验统计公式即可算出每段岸线的最大增水量Tmax,根据评估区域的历史资料建立起风力大小、风向与增水量之间的统计经验公式,公式为Tmax=AFtsinΘ,确定经验参数λ(λ值一般从实测历史风暴潮值中反推求出,取值大致在〇.000459~0.001之间),根据风力大小、风向等计算出风暴潮的最大增水位。[0056]③计算暴露度指数D:暴露性是指人、财、物直接暴露于风暴潮洪水中。每个网格单元内承险体暴露度与沿岸位置与地形有关系,可用承险体所处位置的平均高程D代表。如果不考虑风直接致灾的情形,在风暴潮水位一定的情况下,在最高潮位之上的承险体不受损失,而最高潮位之下的承险体则遭受淹没的损失。[0057]各单元所在位置的平均高程可根据高精度的数字高程地形(DEM)提取。高精度的DEM可通过当地国土测绘局购买或无人机飞行实地扫描之后通过pix4dmapper提取。本步骤通过arcgis利用RasterCalculator(栅格计算器)工具将高于f1的高程值做掩膜处理,然后对低于fi的网格单元利用zonalstatisticsastable(以表格显示分区统计)工具计算得到包含每个网格的高程的平均值属性表,之后按照公共字段链接到图3所示的格网,得到各网格单元的平均高程。[0058]④根据②和③计算各单元的淹没指数。计算依据按照低于风暴潮洪水水位的区域参与计算,高于风暴潮洪水水位的区域不参与计算。本步骤可直接在Arcgis中计算(见表1),在单元网格内的属性表中,计算淹没指数公式为Tmax-D。[0059]⑤计算各网格单元的距离衰减指数f2。由于距海岸线的距离d不一样,其损失程度也不一致。各单元受风暴潮灾害风浪的作用范围有能量差异,根据野外调查,在一定的距离范围内,一般遵循距离衰减原则,离岸线越近,受破坏越大,离岸线越远,破坏性越小。因此,以受风暴潮影响的最大距离为基准,根据各单元格内到海岸线的最短距离d分段计算距离衰减指数:[0060](G)[0061]可根据实际调查确定当地最大风暴潮海水淹没距离,本实例取最大海水淹没距离为SOOnuk为地面粗糙度调节系数,由于沿岸建筑、植被覆盖等地面粗糙程度会使不同海水淹没距离会有所区别,计算依据可参考曼宁糙率系数进行分级(按相对大小在0~1之间)确定。[0062]本步骤通过arcgis获取各单元的质心坐标,通过arcgis的〃最近距离〃工具计算出各质心坐标到海岸线的最短距离d。[0063]⑥计算风暴潮浪对岸线的易损性指数f3。易损性是指承灾体易于受到致灾洪水的破坏、伤害或损坏的特性,反映了各类承灾体对洪灾的承受能力。风暴海浪对各种承灾体的破坏作用是不同的,损失率由于承灾体不同,其损失程度不一样。岸线可分为基岩岸线、建设围堤、码头岸线、砂砾质岸线、养殖围堤和其他类型岸线,根据风暴潮对各类型岸线的破坏程度和作用强度,我们采取一种简单的分类方法,即根据岸线地貌类型、经济价值等对研究区的岸线易损性进行赋值:^0.8-1.0,Z为建设围堤0.6-0.8,z为养殖围堤[0064]f3(z.)=^0.4-0,6,z为码.1-片:线(7)0.2_0.4,z为砂砾岸线、0-0.2,z为其他岸线[0065]本步骤在arcgis中对各岸线类型作缓冲区分析,以风暴潮海水最大淹没距离为半径,以buffer缓冲区工具按照800m执行缓冲区分析,并以f3(z)类型赋属性值,将各类岸线缓冲区的易损性值用Arcgistoolbox分析工具箱下Identify工具赋给各网格单元。[0066]⑦根据以下公式计算各单元的损失率[0067]在以上基础上,针对风暴潮的直接损失,我们提出了单一网格内直接损失率的计算方法:[0068]LR=(fi(F,t,9)-D)*f2(d)*f3(z)(8)[0069]其中,LR(L〇SSRate)为风暴潮的可能性直接损失率,为风暴潮增水指数,D为风暴潮暴露度指数,f2(d)为风暴潮距尚装减指数,f3(Z)为风暴潮浪对岸线的易损性指数。[0070]表1网格属性表损失率计算[0072]在属性表中,公式(3)-(8)通过vb或python语言按照字段值进行分类计算。[0073]计算结果如图5所示在风速分别为20m/秒和25m/秒的情况下,且风暴持续36小时,计算的风暴淹没的范围和可能性损失大小,由于实例中岸线为砂砾质岸线,取其易损性值为0.3。将结果按自然断点法分为6类,图中值越大,损失率越大。从图中可以看出,可能性损失呈现规律性的分布,离海岸线越近,可能性损失值越大;在沿海岸线方向上,受所在单元高程对海水淹没的影响,风险值规律不明显。但在风速20m/秒和25m/秒两种情况下,一方面,后者比前者淹没的范围大,另一方面,靠近岸线的同一网格单元内可能性损失值后者也比前者大。【主权项】1.一种风暴潮灾害风险精细化评估的计算方法,其特征在于包括W下步骤:1)将待评估地区的风暴潮影响区域矢量图层进行网格单元分割;2)根据风向、风速得到该区域的风暴潮暴增水指数fi;3)将各网格单元的平均高程作为各网格单元的风暴潮暴露度指数D;平均高程低于fi的网格单元执行下一步骤;4)根据海岸线类型得到风暴潮浪对岸线的易损性指数f3;5)根据各网格单元距海岸线的最近距离得到风暴潮距离衰减指数f2;6)根据风暴潮暴增水指数fi、风暴潮暴露度指数D、风暴潮距离衰减指数f2和易损性指数f3得到各单元的损失率;损失率值用于表示风暴潮灾害的风险。2.根据权利要求1所述的一种风暴潮灾害风险精细化评估的计算方法,其特征在于所述根据风向、风速得到该区域的风暴潮暴增水指数fi通过W下公式得到:fi(F,t,白)=AFtsin白其中,fi为风暴潮暴增水指数,λ为系数,F为风速,Θ为风向与岸线的角度,t为大风持续时间。3.根据权利要求1所述的一种风暴潮灾害风险精细化评估的计算方法,其特征在于所述根据各网格单元距海岸线的最近距离得到风暴潮距离衰减指数f2,通过W下公式得到:其中,d为某网格单元质屯、坐标到矢量图层中海岸线的最短距离,k为调节系数。4.根据权利要求1所述的一种风暴潮灾害风险精细化评估的计算方法,其特征在于所述根据海岸线类型得到风暴潮浪对岸线的易损性指数f3通过W下公式得到:5.根据权利要求1所述的一种风暴潮灾害风险精细化评估的计算方法,其特征在于所述损失率:LR=(fi(F,t,目)-D)*f2(d)*f3(z)。【文档编号】G06Q10/06GK105844427SQ201610230264【公开日】2016年8月10日【申请日】2016年4月14日【发明人】于良巨,施平,侯西勇,邢前国【申请人】中国科学院烟台海岸带研究所
网友询问留言 已有0条留言
  • 还没有人留言评论。精彩留言会获得点赞!
1