郑州市行政边界通过全国地理信息资源目录服务系统获取，以《郑州市城市更新专项规划》中的地图为基准，对边界进行调整，工业GDP、登记注册企业数、综合能源消费量、GDP、人口数据源于《郑州统计年鉴2023》，DEM数据、植被净初级生产力(Net Primary Production, NPP)下载于地理空间数据云，路网数据采用开放街道地图(Open street map, OSM)道路数据集，污染源企业源于国家级环保备案数据库，避难所、物资库、敏感点、监控点位、企业废水排放量、历史事故率等数据源于郑州市本地综合性环保信息平台，应急演练数据为郑州市突发环境事件应急演练专家打分均值，高中及大学以上文化程度人口数据源于《郑州市第七次全国人口普查公报(第四号)》。

按照刘云熹等提出的突发环境事件应急能力理论框架 [5] ，结合专家访谈与文献分析法，从风险源与敏感性评估、应急准备与响应能力、社会支持与恢复力3个维度出发，构建包含12个指标的突发环境事件应急能力评估指标体系，具体见 表1 ；参考余光辉等研究成果 [3] ，选用坡度和NPP来表征自然因素，选用工业GDP、登记注册企业数、综合能源消费量来表征社会经济因素。

注：参考《企业突发环境事件风险分级方法》，将污染源企业划分为石油化工、有机化学制造，有色金属冶炼、电镀，危险废物处置、医疗垃圾清理，制药、农药生产，一般工业等五类，行业风险系数依次为1.8、1.6、2.0、1.5、1。

$R_i=\alpha \times \frac{S_i}{S_M}+\beta \times \frac{H_i}{H_M}+\gamma \times \frac{D_i}{D_m}$ (1)

式中 $S_i$为污染源密度， $S_M$为 $S_i$最大值， $H_i$为近5年事故发生率， $H_M$为 $H_i$最大值， $D_i$为最近敏感点距离， $D_m$为 $D_i$最小值，权重 $\alpha$、 $\beta$、 $\gamma$分别为0.4、0.3、0.3。

本文将熵值法和层次分析法结合，从而构建突发环境事件应急能力评估指标体系，即分别通过熵值法和层次分析法计算出指标权重，然后对其进行加权融合，该方法可较大程度避免由于主观因素而产生的结果误差 [10] ，计算步骤如下：

(1) 设有 $n$个研究单元， $m$个指标， $x_{ij}$为研究单元 $i$的第 $j$个指标， $x_{ij}$标准化后得到 $y_{ij}$，计算第 $j$项指标的熵值 $e_j$：

$e_j=-\frac{1}{\ln \left(n\right)}{\displaystyle \underset{i}{\sum }{y}_{ij}\ln \left(y_{ij}\right)}$ (2)

(2) 计算客观权重 $w_j$：

$w_j=\frac{1-e_j}{{\displaystyle \underset{j}{\sum }1-e_j}}$ (3)

(3) 通过专家打分构建判断矩阵，采用特征向量法计算客观权重 $w_a$。

(4) 采用乘法合成法平衡主客观权重，求得最终权重 $w_z$：

$w_z=\frac{w_j{w}_a}{{\displaystyle \sum w_j{w}_a}}$ (4)

(5) 计算各区县得分：

$H_i={\displaystyle \underset{j}{\sum }\left(w_z{y}_{ij}\right)}$ (5)

地理探测器是探测空间分异性、以及揭示其背后影响因素的一种统计学方法 [11] ，共包含4个探测器，其中因子探测可用于探测某因子 $X$在多大程度上解释属性 $Y$的空间分异，交互探测用于识别不同因子之间的交互作用，即评估因子 $X_1$和 $X_2$共同作用时是否会增加对因变量 $Y$的解释力，因子探测公式如下：

$q=1-\frac{{\displaystyle \underset{h=1}{\overset{L}{\sum }}{N}_h{\sigma }_h^2}}{N{\sigma }^2}=1-\frac{SSW}{SST}$ (6)

$SSW={\displaystyle \underset{h=1}{\overset{L}{\sum }}{N}_h{\sigma }_h^2},SST=N{\sigma }^2$ (7)

式中 $L$为变量 $Y$的分层，下标 $h$表示第 $h$层， $N$为全区的单元数， ${\sigma }^2$为全区 $Y$值的方差， $SSW$为层内方差之和， $SST$为全区总方差， $q$值越大表示自变量 $X$对 $Y$的解释力越强，反之则越弱。

交互探测通过比较 $q\left(X_1\right)$、 $q\left(X_2\right)$与 $q\left(X_1\cap X_2\right)$的关系，进而将交互类型划分为双因子增强和非线性增强。

$q\left(X_1\cap X_2\right)>M\left(q\left(X_1\right),q\left(X_2\right)\right)$ (8)

$q\left(X_1\cap X_2\right)>q\left(X_1\right)+q\left(X_2\right)$ (9)

式中 $q\left(X_1\right)$、 $q\left(X_2\right)$表示两种因子 $X_1$和 $X_2$对 $Y$的 $q$值， $q\left(X_1\cap X_2\right)$表示交互时的 $q$值， $M\left(q\left(X_1\right),q\left(X_2\right)\right)$表示 $q\left(X_1\right)$和 $q\left(X_2\right)$中的最大值。

根据各开发区、区(县、市)的突发环境事件应急能力得分情况，采用Jenks自然断点法将其划分为低水平、中水平和高水平，并进行可视化表达。如 图1 所示，郑州市突发环境事件应急能力整体呈现出由中心向外围逐渐增强的圈层结构，分异特征明显，其中，巩义、新密等西部县级市和以金水区、郑东新区和经开区等为代表的东北部区域应急能力较强，得分均在54以上，显示出较高的风险防控和应急响应水平；相较之下，位于郑州市中部的荥阳市、中原区、二七区、管城区和新郑市，应急能力则相对较弱，得分一般在38至41之间，反映出这些区域在应急准备、资源调配及社会支持方面存在一定的不足；以中牟县、航空港区等为代表的东部区县和位于北部的上街区、惠济区和高新区以及位于西南的登封市的应急能力处于中等水平，得分在42至53分之间，表明这些区县虽显示出一定的应急响应潜力，但仍需进一步加强风险防控和应急准备。

进一步研究各开发区、区(县、市)的应急能力指标得分情况，以便揭示其在应对突发环境事件时的优势与不足，如 表2 所示，从风险源与敏感性评估来看，中牟县表现最优，除敏感点密度外、其余各项指标得分均在8分以上，航空港区各项指标得分也相对均衡，显示出较强的风险防控能力；相较之下，管城区、二七区和荥阳市在环境风险指数、敏感点密度和黄河排放量等指标上具有明显劣势，这些地区工业活动密集、人类社会活动频繁，突发环境事件发生概率较高，需进一步加强环境风险源企业的监管和防控。

在应急准备与响应能力方面，新密市和巩义市整体较为优秀，各项指标分值均在5分以上，表明这些区域在应急准备方面较为充分，能够在突发环境事件发生时迅速调配资源，为应急救援提供有力保障；与之对应，管城区、中原区和航空港区的缺点则更为突出，在避难设施密度、应急物资储备、监控网络密度及应急演练成效四项指标上的得分均低于3分，反映出这些地区在应急响应能力方面尚需加强。

就社会支持与恢复力而言，金水区在路网可达性、经济发展水平和人力资源素质方面占据显著优势，指标得分均在6分以上，在突发环境事件发生后能够迅速调动资源、恢复社会秩序，为应急管理和恢复工作提供有力支持；对比而言，登封市、新密市和中牟县经济发展水平较低、路网可达性较差，这在一定程度上限制了其社会支持与恢复力的发挥，需加强基础设施建设，提升经济发展水平，优化人口素质，以提高突发环境事件应对能力。

将各影响因素按Jenks自然断点法分为5类、确保类间差异最大化，而后导入地理探测器之中，求得因子探测的 $q$值，如 表3 所示，从 $q$值来看，登记注册企业数(0.489)和坡度(0.473)是影响郑州市环境应急能力的主导因素，登记注册企业数作为社会经济活动的重要指标，其高值区域往往伴随着较高的工业活动和潜在的环境风险，坡度则揭示了自然地理条件对应急能力分布的制约作用，坡度较大的区域可能面临更多的地质灾害和地形限制，不利于应急物资的快速调配和人员的及时疏散；工业GDP (0.301)、NPP (0.179)和综合能源消费量(0.195)也对环境应急能力产生了一定影响，工业GDP的增长往往伴随着工业活动的扩张，增加了环境风险；NPP作为自然生产力指标，其变化可能影响生态系统的稳定性和恢复力，进而影响环境应急的社会支持与恢复力；综合能源消费量则反映了能源消耗水平和污染排放强度，是影响环境应急能力的重要因素之一。

交互探测结果显示，相较于单因子 $q$值，不同影响因素的双因子 $q$值均得到显著提高，说明影响因素的交互作用有效增强了对环境应急能力的解释，如 表4 所示，NPP与工业GDP、NPP与登记注册企业数、NPP与综合能源消费量、工业GDP与综合能源消费量等4类影响因素的交互作用为非线性增强，其余影响因素为双因子增强，其中交互作用最强的为工业GDP与登记注册企业数(0.839)，这一结果反映了社会经济活动的密集程度与环境应急能力之间的紧密联系，工业GDP的增长往往伴随着工业活动的扩张和登记注册企业数的增加，这不仅带来了潜在的环境风险，也对环境应急管理和响应能力提出了更高的要求；交互作用最弱的为坡度与工业GDP (0.549)，这可能是由于坡度作为自然地理条件，主要影响应急物资的快速调配和人员的及时疏散，而工业GDP则更多地反映了社会经济活动的规模和强度，两者虽然都对环境应急能力产生影响，但其作用机理相对独立，因此交互作用较弱。

注：右上角标N的交互类型为非线性增强，其余为双因子增强。

本研究以郑州市16个开发区、区(县、市)为研究对象，从风险源与敏感性评估、应急准备与响应能力、社会支持与恢复力三方面出发，通过熵值法和层次分析法构建突发环境事件应急能力评估指标体系，然后对郑州市应急能力空间格局和能力特征进行分析，最后研究了不同影响因素的作用机理，主要结论如下：

(1) 郑州市突发环境事件应急能力整体呈由中心向外围逐渐增强的圈层结构特征，西部巩义、新密等县级市及东北部金水区、郑东新区等区域综合得分最高，展现出较强的风险防控与应急响应水平；中部二七区、管城区等核心区域能力薄弱，应急资源调配与社会支持存在明显短板；东部中牟县、航空港区及北部上街区等处于中等水平，需进一步强化风险源管控与应急准备。

(2) 郑州市各区域环境能力特征差异显著，中牟县与航空港区在污染源风险值、环境风险指数等方面较为占据优势，管城区、二七区在敏感点密度与黄河排放量等指标上处于劣势地位，新密市、巩义市在避难设施密度、物资储备等指标上表现优秀，金水区在路网可达性与经济发展水平上表现优异，登封市、新密市受基础设施与人口素质的限制较为明显。

(3) 郑州市环境应急能力受自然因素与社会经济因素的交互驱动，因子探测结果表明，登记注册企业数与坡度是影响环境应急能力的主导因素，反映出背后的社会经济活动强度与地形制约效应；交互探测结果显示，工业GDP与登记注册企业数交互效果最强，凸显了经济活动密集区环境风险与应急需求的协同增强。

本文着眼于郑州市突发环境事件应急能力空间格局、区域特征及影响因素的研究，研究结果表明，郑州市不同区域的突发环境事件应急能力差异显著，在风险源与敏感性评估、应急准备与响应能力、社会支持与恢复力3个维度上各有优劣，这可能会导致突发环境事件发生时部分地区难以有效应对。鉴于此，本文提出一系列提升策略，以期为提升郑州市环境应急管理水平提供部分参考：

(1) 重点强化敏感点防护与快速响应：针对金水区、中原区等人口稠密区存在的敏感点防护薄弱、应急资源调度效率不足等问题，建议生态环境局联合街道办事处建立“1公里应急响应圈”，即以敏感点为中心划定辐射范围，配置微型应急物资储备站(含防毒面具、应急通讯设备)，由社区网格员负责日常维护。

(2) 优化企业风险联防联控机制：航空港区、经开区等工业集聚区的核心问题在于工业企业集聚导致风险源叠加、跨部门协同处置能力不足，建议应急管理局组织实施企业分级管控，按照行业风险系数强制要求高风险企业定期开展应急演练，检查企业应急预案，收集危险品库存数据，整合周边救援力量。

(3) 加强生态屏障与风险预警能力：荥阳市、登封市等由于地形条件的限制容易出现监测盲区，引发次生突发环境事件的概率较大，建议自然资源规划局协同水利局开展“生态修复 + 风险预防”工程，即在河流源头区种植净水植物(如芦苇、香蒲)，同时设置生物毒性预警浮标，在水质异常时自动触发报警，为及时开展后续处置工作争取时间。