《Modelling》：Data-Driven Modelling of Structural Interdependencies in Urban Flood-Resilience Support Systems: A Panel Granger Network Approach
编辑推荐：
城市洪涝韧性取决于社会保护、生态条件、经济资源与基础设施服务之间的相互作用。本研究构建面板Granger网络框架，考察2011—2022年京津冀地区13个城市15项韧性支撑指标之间的时序预测关系。年度序列经方向一致性处理、一阶差分与标准化。采用一阶滞后的Dum
城市洪涝韧性取决于社会保护、生态条件、经济资源与基础设施服务之间的相互作用。本研究构建面板Granger网络框架，考察2011—2022年京津冀地区13个城市15项韧性支撑指标之间的时序预测关系。年度序列经方向一致性处理、一阶差分与标准化。采用一阶滞后的Dumitrescu–Hurlin（D-H）检验，并结合名义Benjamini–Hochberg（BH）筛选，在210个有序指标对中识别出26条候选有向链接；采用双向固定效应（TWFE）模型估计带符号的边权重。在保留的链接中，22条（84.62%）连接不同维度，网络在单一弱连通分量内包含三个混合维度社区。四指标Borda排序将总抚养比、医疗保险参保率、R&D经费占GDP比重与污水处理率列于前四位；排水管道密度与卫生技术人员并列第五。留一城市分析保留基线链接的80.8%–96.2%，而留一年份分析仅保留23.1%–57.7%。定向移除实验进一步刻画有向可达性与网络连通性的变化。上述发现建立了跨维度预测组织的经验描述，并揭示所选网络对时间样本构成的敏感性。该框架整合了基于面板的链接筛选、网络刻画与设定诊断，用于城市洪涝韧性支撑指标的结构分析。
研究背景方面，城市洪涝受极端水文气候、土地开发、基础设施服务与社会经济条件交互作用影响，强降雨可超出排水系统负荷、破坏交通与基础设施并影响基本公共服务可及性，且后果会跨越基础设施与制度边界传播，因此理解支撑城市功能连续性的各项条件之间的关系对识别韧性规划中的协调节点具有重要意义。城市洪涝韧性涵盖吸收冲击、维持基本功能、恢复服务与适应恢复过程的能力。既有研究主要沿三条路径展开：一是构建多维指标体系并采用熵权法、层次分析法等多准则技术进行综合测度；二是考察城市与城市群尺度韧性的时空分异；三是通过灰色关联分析、障碍度模型、可解释机器学习与地理探测器识别关联因素。这些工作为韧性状况测度及其空间分异与潜在决定因素提供了经验基础，但指标本身的时间组织问题尚未充分展开——社会保护、公共卫生服务、生态空间、财政资源与基础设施投资可能共同演化并存在滞后预测关系，此类关系可在组分层面表征其方向、连通性与在多维支撑系统中的位置，从而与总量韧性评估形成互补。复杂网络分析为表征组分层面关系提供了工具，已有研究将其用于高密度沿海城市韧性量化与城市内涝风险评估，高维向量自回归模型的网络Granger因果统计方法亦有发展。基于此，本研究分析京津冀13个城市2011—2022年社会、自然、经济与基础设施四个维度15项城市洪涝韧性支撑指标，回答三个问题：哪些滞后预测关系在所选指标间出现并以何种有向网络组织；该网络呈现何种跨维度格局、社区与结构位置；其链接与排序在替代设定与定向扰动下是否稳定。
研究人员开展的研究为：先对负向指标（城市人口密度S1、城镇登记失业率S2、总抚养比S5）作方向一致性变换，再作一阶差分与标准化；采用一阶滞后D-H异质面板Granger非因果检验筛选有序指标对，以名义BH调整q值≤0.10为保留标准；对保留链接以TWFE模型估计带符号描述性权重；进而通过网络连通性、跨维度关系、贪婪模块度最大化社区检测与四指标（PageRank、加权出强度、加权入强度、有向介数）Borda排序刻画结构，并以节点与边移除实验、留一城市、留一年份、阈值变动与权重设定变动评估稳定性。主要结论为：基线网络含15个节点与26条有向候选链接，密度0.1238，正权重10条、负权重16条；构成单一弱连通分量与7个强连通分量（最大含9个节点），有向路径覆盖约59.52%的有序节点对，二值有向全局效率0.2819；总绝对权重2.7253，自然连通度0.0377；跨维度链接占84.62%，绝对权重占比84.96%，形成三个混合维度社区（加权模块度0.4348）；Borda排序前四位为S5、S4、E3、I3，I1与S3并列第五；S5具最高PageRank（0.2639）、加权出强度（0.4681）与加权入强度（0.6249），E3具最高有向介数（0.2592）；节点移除的标准化AUC为相对最大弱连通分量（LWCC）0.3622、有向全局效率0.1727，随机均值分别为0.4211与0.2550；留一城市保留80.8%–96.2%基线链接，留一年份仅保留23.1%–57.7%。其意义在于建立了从面板链接筛选到网络刻画与设定诊断的一体化经验工作流，为城市韧性支撑条件的组分层面结构分析提供可复现描述，并明确揭示结论对时间样本构成的敏感性。
关键技术方法（不超过250字）：研究采用京津冀13城市2011—2022年面板数据，指标主要来自《中国城市统计年鉴》、各市统计公报、《河北经济年鉴》（2012—2019）与《河北统计年鉴》（2020—2023），年降水量作为外部气象协变量纳入补充条件设定，个别缺失值按相邻年份线性插值。方法上依次为：方向一致性变换、一阶差分与全样本标准化；ADF–Fisher与Pesaran CIPS面板平稳性诊断；一阶滞后Dumitrescu–Hurlin异质面板Granger非因果检验与名义BH多重检验筛选；双向固定效应模型估计带符号边权重并辅以Mean-Group估计；复杂网络分析（连通性、全局效率、自然连通度、贪婪模块度社区检测、四指标Borda排序）；以及节点与边定向移除、2000次随机移除与度保持零模型对比、留一城市与留一年份再估计等敏感性实验。
研究结果方面，3.1网络构建结果：以一阶滞后D-H统计量、上尾p值约定与名义BH调整q≤0.10，从210个有序对中保留26条候选链接并附TWFE带符号权重。3.1.1网络组成与整体结构：网络密度0.1238，单一弱连通分量，7个强连通分量（最大9节点），有向路径覆盖约59.52%有序节点对，二值有向全局效率0.2819，总绝对权重2.7253，自然连通度0.0377（二阶近似0.0373）。3.1.2社区结构与跨维度关系：绝对权重无向投影上贪婪模块度最大化得到三个混合维度社区（加权模块度0.4348），跨维度链接22/26（84.62%），其绝对权重占84.96%；最大源—目标权重份额依次为基础设施→社会（17.25%）、社会→基础设施（15.64%）、社会→经济（14.89%）、基础设施→自然（14.52%）、经济→社会（13.47%）。3.2关键节点与结构角色：Borda排序前四位为S5、S4、E3、I3，I1与S3并列第五，其后为E2、N1、E5、N2；S5具最高PageRank、加权出强度与加权入强度（0.2639、0.4681、0.6249），有向介数0.2234；E3有向介数最高（0.2592），S3为0.2198，S4为0.0577；I3加权出、入强度为0.4673与0.1784，I1为0.2520与0.2506，N1加权入强度0.3958。3.3拓扑敏感性：Borda定向节点移除下相对LWCC与有向全局效率的标准化AUC为0.3622与0.1727，随机均值0.4211与0.2550，六项预设比较中BH调整q值分别为0.1810与0.0990；按绝对权重降序移除边的AUC为0.6462与0.3704，按加权边介数移除为0.6128与0.3428，随机均值0.6345与0.3966，六项比较均未达BH调整5%阈值，说明连通性与有向可达性对不同移除规则响应不同。3.4替代设定下稳定性：q≤0.05与q≤0.20分别保留18与34条链接，剔除北京、天津得27条（与基线共享23条，Jaccard=0.767）；留一城市保留21–28条链接、保留80.8%–96.2%基线边（Jaccard=0.724–0.893），留一年份保留7–26条、23.1%–57.7%（Jaccard=0.214–0.484）；Borda排序留一城市相关系数0.763–0.961、Top-5重叠3/5–5/5，留一年份为?0.256–0.728与1/5–4/5；q≤0.01、0.05、0.20下秩相关为0.845、0.859、0.714（Top-5重叠均3/5），仅正权重网络秩相关0.2595、重叠2/5，Mean-Group权重秩相关0.7841、重叠4/5；剔除I3后保留21条链接且均与基线共享，共同14项指标Borda秩相关0.958、重叠4/5。诊断显示210个回归中79个残差Pesaran CD检验名义BH调整q≤0.05，年度截面去均值得19条链接（11条共享基线），26个条件源系数中仅2个满足名义BH调整q≤0.10，两种合成完全零过程下单向程序在96.9%与98.3%的800次重复中至少选中一条边。
讨论部分指出，本研究的核心方法贡献在于对城市韧性支撑系统内纵向关系作组分层面表征，弥补综合指数、城市排序与时空评估难以揭示指标间预测关联方向与组织方式的不足；经验上得到稀疏图、强跨维度表征、三个混合维度社区与差异化指标拓扑位置，并与京津冀区域已有综合韧性得分、空间分异与驱动归因研究形成互补。结构启示上，跨维度链接集中揭示社会保护、公共卫生、经济资源、生态空间与城市基础设施之间的反复出现的界面；第一个混合维度社区联结N2、I2、E1、N1、S3、I3与E5，第二个社区联结E4、S1、S2与I1，第三个社区联结S4、S5、E3与E2；S5、S4、E3、I3居Borda前列，E3与S3具突出有向介数，提示潜在桥接位置；基础设施→社会、社会→基础设施、社会→经济与基础设施→自然的权重份额提示应针对排水与环境服务同社会经济及公共卫生条件的时序协调开展定向分析，未来可结合淹没记录、基础设施运行数据与服务可及性测度检验关联在洪涝事件与恢复期是否持续。移除实验提供拓扑视角：定向节点移除下有向全局效率较早下降而弱连通性较缓，效率比较达名义BH调整q=0.0990，但六项移除对比均未达调整5%阈值，故轨迹仅作描述性证据与后续验证基准。局限与未来研究方面：年度城市指标适于追踪中期演变，事件尺度应用需更高频的降雨、淹没、排水运行与恢复数据及刻画短历时强度与空间差异的降雨指标；推断受短时间维度、截面依赖与双变量D-H筛选设计影响，基线D-H筛选采用城市特定双变量回归而年份固定效应仅纳入TWFE加权与条件敏感性模型，残差依赖检验、年度去均值、年份剔除、条件TWFE与合成零模拟均显示链接选择存在实质变异，26个条件源系数仅2个满足名义BH调整q≤0.10，完全零实验在设定生成过程下频繁误选，故名义BH调整q值仅作本设计下的筛选标准，保留链接解释为候选滞后预测关联，更强推断校准需显式处理共同冲击与短面板的方法；指标系统受可比年度城市序列可得性限制，纳入排水能力、泵站运行、淹没暴露与应急响应可及性可更直接表征洪涝过程，空间显式的基础设施与人口数据可支持城市内部比较；若干图测度依赖绝对权重或无向投影，仅概括带符号有向网络的特定侧面，节点与边移除实验刻画指定扰动下的估计图结构变化，融合物理基础设施、空间约束、运行能力与恢复过程的多层扩展可支持跨统计与工程表征的分析。结论部分译文：基于2011—2022年京津冀13城市面板，本研究为15项城市洪涝韧性支撑指标构建了有向、加权、带符号的候选关联网络。一阶滞后D-H筛选在210个可能有序链接中保留26条，其中22条（84.62%）连接不同维度；所得网络包含一个弱连通分量与三个混合维度社区，记录了所选指标系统内跨维度滞后预测关联占主导地位。Borda分析识别出横跨社会保护、公共卫生、经济资源与基础设施的突出位置：总抚养比、医疗保险参保率、R&D经费占GDP比重与污水处理率居前四位，排水管道密度与每万人卫生技术人员并列第五；结合强度与介数测度，这些排序刻画了所选指标如何贡献于筛选图的连通性与方向性组织。网络实验同时量化拓扑响应与设定敏感性：定向节点移除下相对LWCC与有向全局效率的标准化AUC为0.3622与0.1727，随机均值0.4211与0.2550；高绝对权重边移除对应AUC为0.6462与0.3704，加权边介数移除为0.6128与0.3428，随机均值0.6345与0.3966；六项AUC对比均未达BH调整5%阈值。剔除一个城市时再估计保留80.8%–96.2%基线链接，剔除一个年份时保留23.1%–57.7%，表明时间抽样是网络稳定性的重要决定因素。总体而言，本研究建立了推导组分层面候选关联、考察跨维度网络结构并评估推断关系模式稳定性的集成分析工作流，其主要贡献是对城市韧性支撑指标时间组织的系统网络描述及明确的设定诊断。
