【文献阅读笔记-418 CEUS】基于地理加权可加性引力模型和机器学习的城市内部通勤流边界效应的表示-解释框架

01

论文基本信息

论文题目

A representation-explanation framework for flow-based border effects on intra-urban commuting using geographically weighted additive gravity model and machine learning

论文作者

第一作者：Junyang Zhao

通讯作者：Yongxi Gong

第一单位：哈尔滨工业大学-深圳国际设计学院

期刊名称及分区

JCR分区： 区域与城市规划/地理学  Q1   

中科院分区：  地球科学1区TOP（大类） 地理学1区TOP（小类）

IF=8.9

发表时间及周期

Received： 10 June 2026

Revised： 30 July 2026

Accepted：28 August 2026

周期： 约 79天

关键词

城市内部通勤；边界效应；表征-解释框架；地理加权加性引力模型

Highlights

• 一种表示-解释框架量化了基于流的边界效应。

• 提出了基于流的地理加权可加重力模型（FB-GWAGM）。

• 一个加性重力模型将边界厚度表示为等效公里数。

• 行政边界在空间上异质性超过自然边界。

• SHAP 属性将边界厚度表示为公里；交通连通性占主导地位。

DOI: https://doi.org/10.1016/j.compenvurbsys.2026.102516

摘要

边界效应会抑制人员和货物的跨边界流动，正是这种抑制作用推动了边界两侧建成环境和经济社会条件的逐渐分化。因此，从流动中衡量边界效应比从下游结果（如房价差异或服务供给缺口）中推断它们更接近边界效应本身。为了应对基于流动的边界效应本质上不可观测且无法系统表征的挑战，本研究提出了一种表征-解释框架。在表征阶段，采用将每次边界穿越视为物理距离的加性增量来处理的基于流动的地理加权加性重力模型（FB-GWAGM）。在解释阶段，采用具有 SHapley Additive exPlanations 分析的可靠性加权随机森林模型来识别影响边界厚度的建成环境和经济社会因素。 基于模型结果，识别并量化了边界厚度的空间异质性，并在同一度量标准下分解了建成环境和经济社会因素对此异质性的贡献，从而实现了不同边界类型和因素维度的直接比较。这些发现表明，该表征-解释框架为量化基于流量的边界效应提供了一种可复制且具有物理可解释性的方法，深化了对边界如何塑造城市空间互动的理解，并揭示了边界厚度异质性的社会经济依存性。

03

研究背景及研究问题

3.1 研究背景

• 城市通勤是城市内部最频繁、最具刚性的空间移动形式之一，其形成受土地利用多样性、职住平衡、交通可达性、房价、地方服务设施以及距离摩擦等因素影响。

• 与这些连续渐变因素不同，行政边界和自然屏障（如山体、河流）会在边界位置产生更为离散、突出的互动断裂，对跨界通勤形成额外阻滞。

• 既有研究多通过房价差异、公共服务差距、价格空间相关性削弱等边界效应的“后果”来测量边界，但论文强调：边界效应首先作用于跨界流动，建成环境差异只是长期累积后的结果，因此从流动中测量更贴近边界效应本体。

• 另一方面，边界效应无法从原始流量数据中直接读取，它只是嵌入在OD通勤流中的隐含阻抗项，必须借助模型估计才能恢复。

• 广州同时具有明显的制度性边界（区级行政边界）与自然性边界（山体与河网），为考察多类型边界如何共同塑造城市内部通勤提供了典型场景。

3.2 原文提出的三个研究缺口

• 基于流量的边界效应不可直接观测，难以同时把它作为“通勤流的解释变量”和“被建成环境/社会经济因素解释的因变量”纳入同一分析链条。

• 边界效应具有明显空间异质性，但既有方法多为全局估计，难以揭示不同OD对之间的差异。

• 边界效应本质上是出行阻抗的一部分，应以距离等价值表达；既有乘性重力模型通常给出无量纲系数，难以用公里解释，也难以用同一单位量化各因素贡献。

3.3 研究问题（RQs）

RQ1： 如何从通勤流中表征边界效应，并在统一框架下进一步解释其成因？

RQ2： 边界厚度的空间异质性表现为何？哪些建成环境与社会经济因素与这种异质性相关？

RQ3： 如何将边界厚度及其影响因素的贡献统一表示为距离等价值单位？

04

数据源与方法

1.数据源

• 通勤流数据：采用中国移动广州2023年6月工作日匿名手机信令数据。识别18—65岁工作年龄人口的居住地与工作地，在100米网格尺度提取547,156条通勤OD记录，再聚合到PSU尺度得到24,076个正流OD；剔除流量≤3的弱联系后，最终用于建模的样本为18,672个OD对。

• 边界与解释变量数据：行政边界来自行政区边界；自然边界包括由30 m SRTM DEM识别的山体屏障，以及国家基础地理信息库提取的五级及以上河流和水体。解释变量还整合了OSM道路、Amap POI/公交/地铁/路径规划、企查查企业数据及安居客/贝壳房价数据。

  

2.研究区域

研究对象为广州，全市常住人口约1870万，包含11个行政区。分析基本单元为182个街道/镇级次行政单元（subdistrict, PSU），这是中国城市交通规划、公共服务配置和土地利用管理的重要基层尺度。

3.研究方法

• 方法框架：第一阶段使用FB-GWAGM，将边界穿越次数作为对欧氏距离的加性阻抗增量，并在四维OD空间内进行地理加权局地拟合，得到以公里表示的行政与自然边界厚度；带宽由AICc与CV-MSE联合确定，最终取K=3000。第二阶段以局地边界厚度为因变量，构建可靠性加权随机森林并结合SHAP解释不同建成环境、社会经济和自然障碍强度与边界厚度异质性的关联。

• 
  

• 模型验证：论文还设置了2×2消融实验（全局/局地 × 乘性/加性）并采用五折空间分块交叉验证，对表征能力进行比较。

  

05

研究结果

• 模型表征效果显著优于全局模型：带宽联合评估后选取 K=3000。作为局地拟合初值的全局乘性模型，R2=0.621；进一步采用局地加性模型后，平均局地 R2=0.732，范围为0.234—0.888，说明广州通勤中的边界效应具有明显空间异质性，单一全局参数不足以刻画。

  
  

• 消融实验验证“局地+加性”组合最优：五折空间分块交叉验证中，M4（即FB-GWAGM，局地加性）取得最高平均测试 R2=0.609 和最低平均RMSE=204.59；相较之下，全局加性模型M2在某一保留区块上出现严重失稳，均值 R2为负，表明单一全局加性参数集空间迁移性不足。

  

• 边界厚度具有可接受的估计可靠性：行政边界厚度与自然边界厚度的局地显著比例分别为98.51%和84.19%，平均局地 R2为0.732，因此二者均进入第二阶段解释模型。附录中的上限敏感性分析和固定γ敏感性分析表明：中心分布与总体空间排序较稳定，尤其是行政边界厚度结果较稳健。

  

• 两类边界厚度的平均水平不同：局地估计显示，行政边界厚度均值为5.861 km/次穿越，自然边界厚度均值为3.010 km/次穿越。这意味着在平均意义上，穿越一次行政边界相当于额外增加约5.861公里出行阻抗，而穿越一次自然边界相当于增加约3.010公里。

• 空间格局清晰分化：行政边界厚度高值主要出现在从化与花都/白云/增城交界以及白云—黄埔边界，局地值最高可达接近30 km；而天河、海珠、荔湾、越秀等核心城区之间边界厚度最低，仅约2.83—4.01 km。自然边界厚度高值则集中在白云山—王子山山体周边，显示陡坡与稀疏穿越设施叠加时边界阻抗最强；河流边界多表现为中低值，南部平坦地区接近零。

  

• 第二阶段解释模型表现较好：加权调参后的随机森林在测试集上，行政边界模型和自然边界模型的 R2分别为0.748和0.735。

• 关键关联因素存在显著差异：

  对行政边界厚度，SHAP显示最重要变量是BRD_WR（56.8%），其次为PORRD（12.1%）、ED_R（5.5%）、CDI（3.9%） 和 ED_W（3.9%）。

  

  对自然边界厚度，最重要变量是BRD_WR（38.7%），其次是BARRIER_SEV（27.9%） 和 PORRD（9.0%），再是 ED_W（3.7%） 与 ED_R（3.4%）。

  

• 由于第一阶段以公里度量边界厚度，SHAP贡献也能直接以公里表示。例如，行政边界模型中 BRD_WR 的低值与高值分别对应约 +9.256 km 与 −4.217 km 的平均贡献；自然边界模型中 BRD_WR 低值与高值对应约 +3.722 km 与 −2.001 km，体现出清晰的非线性阈值特征。

  

06

主要讨论内容

• 行政边界比自然边界更“厚”：论文认为，这主要源于两类边界的机制差异。行政边界的阻滞机制是制度性的，区级治理在规划权、公共服务供给和土地利用管理上的分割，会在边界两侧形成服务覆盖、基础设施投入和空间组织上的不连续性，因此其平均阻抗高于自然边界。

• 行政边界厚度并非仅由“是否有边界”决定，而与周边城市化条件密切相关：虽然所有区界在制度上具有同等行政属性，但局地结果呈现出明显的核心—外围梯度：中心城区边界厚度较低，外围边界较高。论文据此指出，行政边界效应还与跨界交通供给、道路网络结构、企业多样性与功能分化等可调节的建成环境条件相关。换言之，边界本身不变，但周边通达性和功能互补性不同，会对应不同的跨界渗透水平。

• 自然边界的平均厚度较低，但局地差异更集中于特定屏障位置：自然边界的作用机制主要是物理性的，山体陡坡、水体宽度等在特定地段形成显著阻断，而非像行政边界那样沿整条边界均匀发挥作用。因此，自然边界呈现“均值较低、局地极值较高”的特征，特别是在山体周边。

• 自然边界并非只受自然条件支配，也与周边可达性环境共同关联：在自然边界模型中，BARRIER_SEV的重要性仅次于BRD_WR，说明地形/水系强度确实关键；但公交连接和道路结构仍然占有较大权重，表明自然边界的“可穿越性”还取决于桥梁、道路、轨道和公交网络如何将两侧空间重新连通。

  

• 政策启示强调“差异化边界治理”而非简单消除边界：论文明确指出，行政边界与自然边界都具有治理、生态和空间组织功能，因此政策目标不是取消边界，而是在必要通勤联系被过度阻碍的边界段提升通透性。对高厚度自然边界，应优先改善跨障碍联通设施，同时兼顾生态保护；对行政边界，则应加强相邻辖区在交通供给、道路连续性和互补型土地利用上的协同。广州的两个河流案例进一步说明：同为河流边界，不同厚度对应不同规划响应——外围从化的流溪河更需要审慎提升必要跨河连接，而中心城区珠江段则更应重在维护和协调既有高通透网络。

07

结论

本文基于广州18,672个街道尺度通勤OD对，构建了一个面向基于流量的边界效应的“表征—解释”两阶段框架。第一阶段用FB-GWAGM从通勤流中提取以公里表示的局地边界厚度，第二阶段以可靠性加权随机森林+SHAP解释其空间异质性。结果表明，广州通勤中的行政边界厚度（均值5.861 km）高于自然边界厚度（均值3.010 km），且二者都存在显著空间异质性；行政边界厚度更多与交通供给、道路结构和企业多样性等城市条件相关，自然边界厚度则同时与自然障碍强度和交通连通性相关。论文的核心贡献在于：不仅把边界效应从隐含流量中恢复出来，还将其及其影响因素统一转化为距离等价值单位。作者也指出，研究受限于单一时间窗口和低频流量样本剔除，未来可向多城市和纵向时序拓展。

08

亮点汇总

• 提出“表征—解释”一体化框架：把边界效应同时作为通勤流中的隐含阻抗项和需要解释的结果变量，形成完整分析链条。

• 从“流量”而非“下游结果”测量边界效应：相较于房价差异、服务差距等后果变量，更贴近边界效应的直接作用过程。

• 采用加成重力模型表达边界穿越阻抗：将每次边界穿越视为对物理距离的额外增加，使边界厚度具有明确的公里意义。

• 引入地理加权局地估计：在四维OD空间中进行局地拟合，成功识别行政边界与自然边界的空间异质性。

• 第二阶段使用可靠性加权RF与SHAP：考虑第一阶段估计不确定性，并揭示变量的非线性、阈值和方向性贡献。

• 实现“边界厚度—因素贡献”同单位比较：不仅边界厚度以公里表示，SHAP贡献也能以公里解释，增强了跨边界类型、跨因素维度的可比性。

• 发现广州通勤边界效应具有差异化机制：行政边界更厚、呈核心—外围梯度；自然边界平均更薄但在山体等局部区域高度集中。

• 提出差异化规划启示：不主张简单消除边界，而是依据不同边界段的厚度和功能背景，实施有针对性的跨界联通与协调治理。