1、2023 年第 6 期何可沁,等:空间大数据支持下城市安全生产风险评估方法引文格式:何可沁,程南炜,邓敏 空间大数据支持下城市安全生产风险评估方法J 测绘通报,2023(6):167-171 DOI:10 13474/j cnki 11-2246 2023 0188空间大数据支持下城市安全生产风险评估方法何可沁,程南炜,邓敏(中南大学地球科学与信息物理学院,湖南 长沙 410083)摘要:现有城市安全生产风险评估研究大多基于统计调查数据微观地评估单个企业或园区,没有充分考虑安全生产事故空间效应,难以适应于城市尺度安全生产风险评估。本文基于 H-V 风险评估框架综合分析区域致灾因子危险性、区域目
2、标暴露度和区域救援适应性等多个方面,构建城市安全生产风险评估体系,在多源地理空间数据支撑下量化分析城市安全生产风险指标,揭示城市安全生产风险分布格局,为城市安全生产防控与治理提供决策知识服务。以昆山市为研究区域,实现空间大数据支持下的城市安全生产风险评估。试验结果表明,昆山市安全生产高风险区域主要集中于玉山镇东部和南部,其中玉山镇南部致灾因子危险性较高,玉山镇东部致灾因子危险性与区域目标暴露度均较高。关键词:地理大数据;城市安全生产;风险评估中图分类号:P208文献标识码:A文章编号:0494-0911(2023)06-0167-05Urban safety production risk a
3、ssessment supported by spatial big dataHE Keqin,CHENG Nanwei,DENG Min(School of Geosciences and Info-Physics,Central South University,Changsha 410083,China)Abstract:Most of the existing studies on urban safety production risk assessment are based on statistical survey data to assess a singleenterpri
4、se or park,and without sufficiently considering the spatial effect of safety production accidents,which is difficult to adapt to theurban scale safety production risk assessment Based on the H-V risk assessment framework,this paper constructs the urban safetyproduction risk assessment system by comp
5、rehensively analyzing regional hazard factors,regional target exposure and regional rescueadaptation With the support of multi-source geospatial data,the urban safety production risk index is quantified based on spatialanalysis to reveal the distribution pattern of urban safety production risk,which
6、 could provide decision-making knowledge service for theprevention and control of urban production safety Taking Kunshan as the study area,the urban safety production risk assessmentsupported by spatial big data is realized The experimental results show that the high-risk areas of safety production
7、in Kunshan aremainly concentrated in the east and south of Yushan Town,and the risk of disaster factors is higher in the south of Yushan Town,whilethe risk of disaster factors and regional target exposure are higher in the east of Yushan TownKey words:geographic big data;urban safety production;risk
8、 assessment随着我国工业化与城镇化进程的快速推进,集约分布、迅速发展的生产企业给社会经济增长注入了强劲活力,同时也给城市安全运行带来了新挑战1-2。当前,我国生产企业数量增长快速,存在工作场地密集、地下空间交错和燃气管线复杂等特点,传统性和非传统性风险因素持续增加3,城市安全生产风险随之不断积累,生产事故概率和危害逐步增大。近年来,重大生产企业事故表明,当前城市安全生产风险管控能力滞后于城市生产建设速度4。生产企业类事故呈现事件突发、危害连锁的时空特征5,不仅对生产企业自身造成重大损失,而且给人民生命财产安全和城市公共服务设施保障带来严重威胁6。2018 年党中央、国务院印发了 关于
9、推进城市安全发展的意见,提出“强化安全风险管控,对城市安全风险进行全面辨识评估”的总体要求7-8。对城市开展大范围安全生产风险评估,构建系统、科学、可泛化的城市安全生产风险评估体系,发现生产企业高危隐患区域,对减少或避免重特大安全事故发生、确保城市安全生产及维护社会稳定具有重大现实意义9。目前,国内外专家学者从不同角度对风险进行描述10-11。文献 12提出风险是可能发生的损失程度与不确定性;文献 13 认为风险是危险事件发生概率;文献 14指出风险是致灾因子、承灾体及761收稿日期:2022-06-05作者简介:何可沁(1997),女,硕士,主要研究方向为时空数据挖掘和风险评估。E-mail
10、:1575198842 qq com测绘通报2023 年第 6 期其脆弱性协同作用的效果,即风险=致灾因子承灾体脆弱性;联合国国际减灾战略(ISDR)组织提出风险是由人为或自然因素等致灾因子作用于承灾体脆弱性而引发破坏事件的概率,即风险=致灾因子脆弱性;国际减轻自然灾害十年(IDNDR)组织认为风险是指某区域特定时段内某个或系列事件造成破坏的严重程度15-16,即风险=事件发生的概率脆弱性。综上,现有风险评估多基于危险性(H)和脆弱性(V)维度构建风险评估指标,为本文城市安全生产风险评估提供了坚实理论框架。安全生产风险是指在生产过程中由安全隐患可能引发的、造成大面积生命财产安全遭受重大损失的灾
11、难性社会后果,通常包括火灾和爆炸等安全生产事故。现有城市安全生产风险研究主要基于企业点要素,即关注具体的、单一的危险源评价。如文献 15采用层次分析法建立企业综合环境风险指数,对松木岛化工园区内部企业个体进行分级评估;文献 17 收集了某工业园区中 5 家企业重大危险源分布及地理环境数据,综合个人风险、社会风险和潜在生命损失评估企业风险;文献 18 针对钢铁企业建立风险评估指标体系,运用层次分析法与熵权法对钢铁企业进行风险评估与评级。这类研究对数据精度和粒度要求高,调查数据内容涉及广且获取困难,较少顾及危险性的空间叠加效应和周围经济社会的脆弱程度,难以支持城市尺度的安全生产风险评估需求。随着位
12、置感知技术与数据库技术的快速发展,多源泛在地理大数据为全面认知城市“画像”提供了丰富数据源。此外,利用空间分析方法探索企业、人口、医疗、消防等城市要素分布格局与配置关系成为空间评估与服务决策的常用手段19-20。因此,本文基于 H-V 风险评估框架构建融合致灾因子危险性、区域目标暴露度和区域救援适应性的城市安全生产风险指标体系;在此基础上,协同多源地理空间大数据和空间分析技术精细量化大范围区域安全生产风险,为城市安全生产风险重点管控和企业分布格局优化等服务提供辅助决策知识。1研究区概况与数据预处理1.1研究区概况昆山市位于江苏省东南部,地处上海市和苏州市之间,内部河网密布,地势平坦,空气质量良
13、好。近年来,昆山市受上海“浦东效应”和本地人才计划等政策的大力支持,吸引了大量生产企业在此落地,经济发展快速。截至 2022 年,昆山市已连续17 年蝉联全国经济发展百强县首位。然而,生产企业的快速发展同时也带来了诸多安全生产问题,火灾和爆炸等潜在生产事故极大地威胁着人民生命财产安全与社会稳定。当前,昆山市正处于城市安全风险评估的重要阶段,因此本文选取昆山市作为研究区域,从安全生产的角度对城市安全进行评估。1.2数据源和数据预处理研究数据如图 1 所示,包括:通过 OSM 获取的道路和行政区划等基础地理矢量数据;从高德地图获取的昆山市医院、派出所和消防队位置信息;从 WorldPop 获取的昆
14、山市人口分布信息;从市委市政府联合各部门实施的风险评估方案中获取的昆山市火灾、爆炸等企业生产危险事件清单,以及操作不当、设备故障等企业生产隐患事件清单。数据预处理主要包括冗余处理、空间坐标系统一、数据裁剪与合并、数据归一化和空间数据网格化等步骤。具体的,仅保留一条原始数据中重复企业数据行;将企业、危险事件与隐患事件进行空间匹配,使之具有空间位置信息;运用 ArcGIS 软件,将数据统一至 WGS-84 坐标系下,通过昆山市行政区划空间信息,对企业数据、POI 数据、道路数据和人口信息进行裁剪提取;将昆山市划分为 1 km 1 km的规则网格,并作为研究单元。图 1研究数据概况8612023 年
15、第 6 期何可沁,等:空间大数据支持下城市安全生产风险评估方法2研究方法现有城市安全生产风险评估多基于单个生产企业或园区设计,对城市尺度风险分布格局和空间效应认识不足,难以有效支持大范围风险评估任务。本文将生产事故空间效应嵌入致灾因子危险性,进一步将承载体脆弱性细化为区域目标暴露度和救援适应性,构建城市安全生产风险评估指标体系。在此基础上,基于 H-V 风险评价框架计算城市安全生产风险值,综合指标体系如图 2 所示。图 2综合指标体系2.1城市安全生产风险指标体系构建H-V 风险评价框架分为致灾因子危险性和承载体脆弱性。其中,承载体脆弱性可进一步细分为区域目标暴露度和区域救援适应性21。2.1
16、.1致灾因子危险性致灾因子危险性指可能会危及城市安全生产的灾害源因素,本文选取生产企业自身危险性和空间分布效应两个测度。企业自身危险性由网格内企业自身火灾、爆炸等危险及操作不当、设备故障等隐患的联合分布衡量,公式为Risk(a)=j=1danj(a)danjhidj(a)hidj(1)式中,danj(a)和danj分别为网格单元 a 第 j 类危险等级和研究区域第 j 类危险均值;hidj(a)和hidj分别为网格单元a第j类隐患等级和研究区域第j类隐患均值。考虑到生产企业聚集分布导致的潜在危险和空间叠加,网格单元 a 空间分布危险性由自身及周围网格单元生产企业空间密度的邻近加权描述,即Vul
17、(a)=C(a)area(a)+W(a,b)C(b)area(b)(2)式中,C(a)和 area(a)分别为网格单元 a 内生产企业数量和面积;C(b)和 area(b)分别为邻近网格单元 b内生产企业数量和面积;W(a,d)为网格单元a和b的标准正态分布距离权重。在此基础上,网格a的致灾因子危险性表示为d(a)=K(Vul(a)K(Risk(a)(3)式中,K()为数据归一化。2.1.2区域目标暴露度区域目标暴露度指危险环境中居民生命财产和事物的暴露程度,是安全生产事故风险存在的必要条件21,本文选取人口暴露度和社会经济暴露度两个因素。人口暴露度由网格单元 a 内人口密度衡量,公式为vul
18、_pop(a)=pop(a)/area(a)(4)式中,pop(a)和 area(a)分别为网格单元 a 内的人口数量和面积。社会经济暴露度由网格单元 a 内基础设施密度衡量,公式为vul_gdp(a)=C_poi(a)/area(a)(5)式中,C_poi(a)和 area(a)分别为网格单元 a 内基础设施数量和面积。网格 a 的区域目标暴露度表示为c(a)=w1K(vul_gdp(a)+w2K(vul_pop(a)(6)式中,w1和 w2表示熵权法计算得到的指标权重。2.1.3区域救援适应性区域救援适应性指受灾区域得到救援的速度或可能性大小,保证了承载体免受或少受事故影响的程度21。本文
19、选取医疗救援可达性、消防救援可达性和道路通达性 3 个因素。网格单元 a 的医疗救援可达性和消防救援可达性分别由 a 中心距综合医院和消防站最短路网距离衡量,公式为Dist_h(a)=min(net_dis(a,hos)(7)Dist_f(a)=min(net_dis(a,fir)(8)式 中,min()为 取 最 小 值;net_dis(a,hos)和net_dis(a,fir)分别为网格单元 a 中心距所有医院和消防站的路网距离数组。道路通达性由网格单元 a 内道路密度衡量,公式为Ral(a)=R_l(a)/area(a)(9)式中,R_l(a)和 area(a)分别表示网格单元 a 的路
20、网总长度和面积。网格 a 区域救援适应性表示为p(a)=(w3K(Dist_h(a)+w4K(Dist_f(a)Ral(a)(10)式中,w3和 w4表示熵权法计算得到的指标权重。2.2城市安全生产风险评估H-V 风险评估框架认为风险是灾害危险性和承961测绘通报2023 年第 6 期灾体脆弱性等风险指标的联合分布,即风险(R)=危险性(H)脆弱性(V)。基于本文建立的城市安全生产风险指标体系,城市安全生产风险由致灾因子危险性、区域目标暴露度和救援适应性空间联合分布衡量,即R(a)=d(a)c(a)p(a)(11)式中,d(a)、c(a)和 p(a)分别为致灾因子危险性、区域目标暴露度和救援适
21、应性。当 d(a)、c(a)和 p(a)同时增大或减小时,R(a)随之增大或减小;当 d(a)、c(a)和 p(a)变化不一致时,R(a)取决于三者间相对变化速率;当 d(a)或c(a)=0,即不存在潜在灾害或无人员和设施财产分布时,R(a)=0。3结果与分析本文将昆山市划分为 1 km1 km 的规则网格,利用多源空间大数据进行昆山市安全生产风险评估,结果如图 3 所示。图 3昆山市安全生产风险评估结果图 3(a)表示综合昆山市企业自身火灾、爆炸等自身危险性与空间分布危险性的致灾因子危险性分布。昆山市安全生产致灾高危区域主要位于玉山镇南部、东部地区和经济技术开发区西部地区。大量的化学化工企业
22、,金属制造业等高风险企业在此聚集,而且企业聚集所带来的风险性呈现出“1+12”的空间叠加效应,威胁周围企业及其他城市要素安全,从致灾因子角度需要对 3 个地区尤其是黑色高危险地区进行重点关注。图 3(b)表示综合昆山市人口暴露度和社会经济暴露度的区域目标暴露度评价结果。昆山市区域目标暴露度高值主要位于玉山镇东部和经济技术开发区,陆家镇、张浦镇、千灯镇也呈现小范围高值集中分布。可以看出,玉山镇和经济技术开发区作为昆山市经济较为发达地区,城市基础设施与人口分布相对稠密,一旦发生企业事故灾害,将严重威胁周边人民及社会财产安全。图 3(c)表示昆山市区域救援适应性分布,即综合道路通达性、消防救援可达性
23、和医疗救援可达性的企业事故灾害应对能力。黑色表示应急能力较强,白色表示应急能力较弱。图中表明昆山市风险应急能力具有显著异质性,其中,玉山镇对城市风险应急能力较好,而经济技术开发区风险应急能力较弱。图 3(d)表示融合致灾因子危险性、区域目标暴露度和区域救援适应性的昆山市企业安全生产风险评估结果,风险最高区域主要位于玉山镇东部和南部。一方面,昆山市城市居民、基础设施和社会财产大量聚集于玉山镇东部,即使该地应急能力较强,但是当前应急能力难以有效应对致灾因子危险性和目标暴露度。因此,该地除了加强企业安全生产工作外,还需进一步提高应急救援能力;另一方面,虽然玉山镇南部并未聚集大量居民与财产,但是该地企
24、业自身和空间效应产生的致灾因子风险极高,亟须加强消防、医疗和道路通达性,以应对高风险突发企业灾难事故。综合图 3,昆山市人口和社会暴露度分布较广,而致灾因子危险性较高区域、区域救援能力较低区域和安全生产风险较高区域具有相似分布,可知本文方法能够有效顾及致灾因子、承载体和系统环境特征,且以致灾因子危险性和城市环境救援适应性为主要影响因素,为主动发现城市安全生产高风险区域、降低灾害发生概率甚至遏制灾害发生提供先验认知知识。0712023 年第 6 期何可沁,等:空间大数据支持下城市安全生产风险评估方法4结语本文顾及灾害空间效应与城市设施空间分布,基于 H-V 风险评估框架,构建了致灾因子危险性、区
25、域目标暴露度和区域救援适应性的城市安全生产风险评估方法,然后利用多源空间大数据量化评估了昆山市安全生产风险分布格局。可以发现,昆山市安全生产高风险分布于人民与社会财产较为集中的玉山镇东部和南部;同时这些区域致灾因子危险性较高且救援适应性较弱,存在连锁增强的潜在危险,需要重点优化与防控治理。基于 H-V 风险评估框架和空间分析技术设计了城市安全生产风险评估方法,在多源地理空间数据支撑下适用于大范围研究区安全生产风险评估任务。然而,基于规则网格作为地理分析单元,忽略了城市空间天然结构。此外,对企业自身危险性数据依赖性较高,如何基于事故发生后果危险性设计致灾因子危险性需要进一步验证研究。参考文献:1
26、张超,徐凤娇,王皖,等 城市公共安全风险分析标准研究 J 标准科学,2020(6):19-24 2乔庆梅 我国安全生产监督管理问题探析 J 中国软科学,2006(6):20-30 3黄健 城市安全生产风险评估体系构建及存在问题 J 化工进展,2019,38(S1):288-292 4解读 危险化学品生产、储存装置个人可接受风险标准和社会可接受风险标准(试行)J 劳动保护,2014(8):44-46 5胡馨升 化工园区事故风险评估方法及应用研究 D 北京:中国矿业大学(北京),2018 6朱晓莉 化工园区区域评价方法的选择及风险分析 D 天津:天津理工大学,2016 7中华人民共和国国家发展和改
27、革委员会 基于风险检验的基础方法:SY/T 67142008 S 北京:石油工业出版社,2008 8中共中央办公厅 国务院办公厅印发关于推进城市安全发展的意见 J 中国安全生产科学技术,2018,14(1):48 9ZHANG S C Discussion into safety risk evaluation andaccident investigation methods J Safety Health Environment,2015(12):7-10 10 LUO Y Risk analysis and safety assessment M 2nded Beijing:Chemic
28、al Industry Press,2010 11 Shenzhen Emergency Management Office Strengtheningtheriskmanagementandconsolidatingthesafetyfoundationto carry out the urban public safety riskassessment and publish the public safety white paper J China Emergency Management,2014(10):27-33 12 WILSON R,CROUCH E ARisk assessm
29、ent andcomparisons:an introduction J Science,1987,236(4799):267 13 JOHNSTON R J,GREGORY D R,PRATT G,et al Thedictionary of human geographyM S l :Wiley-Blackwell,2009 14 孙爱军 工业园区事故风险评价研究 D 天津:南开大学,2011 15 田园 基于层次分析法的化工企业环境风险评价研究 D 大连:大连理工大学,2015 16 王国栋,颜爱华,侯蕊芳,等 城市安全风险评估方法体系研究及实践 J 中国安全生产科学技术,2019,15(7):46-50 17 张洪姣,侯蕊芳 工业园区重大危险源风险评价实例 J 中国公共安全(学术版),2019(1):33-36 18 贾秋璇,董宪伟 钢铁企业安全风险评估指标体系研究 J 华北理工大学学报(自然科学版),2022,44(1):117-123 19 蒋雯京,程春梅,张艳蓓,等 基于 GIS/AHP 集成的浙江省洪涝灾害风险评估 J 测绘通报,2019(2):125-130 20 华一新 空间决策支持的技术框架J 测绘通报,2015(2):1-4 21 葛全胜,邹铭,郑景云 中国自然灾害风险综合评估初步研究 M 北京:科学出版社,2008(责任编辑:马保卫)171