Page 98 - 《华中农业大学学报（自然科学版）》2024年第1期

92 华 中 农 业 大 学 学 报 第 43 卷

居民的游憩特征数据，利用公园的平均游憩时长除 公园的单位面积游憩强度的最小值， t max 为所有综合
以公园面积，并进行标准化得到公园的单位面积游 公园的单位面积游憩强度的最大值。
憩强度标准化值，之后构建公园的单位面积游憩强 基于以上内容构建改进多出行模式的两步移动
度指数对其服务供给规模进行修正 ［16］ ，具体计算公
搜索模型，分为以下 2 个步骤：第一步，计算公园的
式如下： 服务能力，由于居民到达综合公园采用的交通方式
S j = s j(1 + T j) （1）
不同，因此相同时间阈值下，公园的服务范围会因
t j - t min
T j = （2） 为居民所采用的交通方式而产生改变，本研究根据
t max - t min
式（1）~（2）中， S j 为公园 j 的服务供给规模， s j 为 前面识别的采用步行、骑行以及车行方式到达公园
公园 j 的面积， T j 为公园 j 的单位面积游憩强度标准 的用户比例乘以累计的经过高斯函数衰减之后的
化值， t j 为公园 j的单位面积游憩强度，用公园的平均 用户数量，并进行叠加以计算公园的综合服务能
游憩时长与公园面积的比值来衡量， t min 为所有综合 力值：
S j
R j = （3）
)
)
)
B 1∑ G( d k，j，1 ，d 0 × P k，1 + B 2∑ G( d k，j，2 ，d 0 × P k，2 + B 3∑ G( d k，j，3 ，d 0 × P k，3
k ∈{d k，j，1 ≤ d 0 } k ∈{d k，j，2 ≤ d 0 } k ∈{d k，j，3 ≤ d 0 }
式（6）中， I i，k，l 表示居住单元网格 i 公园绿地可
式（3）中， R j 表示公园 j的服务能力值； B 1 、 B 2 、 B 3
分别表示步行、骑行以及车行的叠加权重； d k，j，1 、 达性与居住人数的 Moran’s I（莫兰指数）， Moran’s
d k，j，2 、 d k，j，3 分别表示居民点 k到公园绿地 j以步行、骑 I 大于 0 或者小于 0 分别表示公园绿地可达性和居
行、车行方式出行的时间； d 0 表示时间阈值， P k，1 、 住人数存在正、负空间关联性，等于 0 则表示不存在
P k，2 、 P k，3 分别表示步行、骑行以及车行在阈值范围内 空间关联性 ［27］ ； Y i，k 表示居住网格单元 i 的公园绿
--
居民点 k的需求规模，利用居民点 k的居住人数表示。 地可达性， Y j，l 表示居住网格单元 j 的居住人数；Y k
--
G( d k，j，1 ，d 0) 、G( d k，j，2 ，d 0) 、G( d k，j，3 ，d 0) 为高斯函数， 表示公园绿地可达性的均值，Y l 表示居住人数的均
值； σ k 为公园绿地可达性的方差， σ l 为居住人数的方
以G( d k，j，1 ，d 0) 为例，计算公式如下：
差； W i，j 为居住网格单元 i 和居住网格单元 j 之间的
2 ( )
ï ï ì - × d k，j，1 2 - 1 空 间 权 重 矩 阵 ，在 空 间 权 重 矩 阵 的 构 建 上 ，以
1
ï ï e
2
d 0
G( d k，j，1 ，d 0) = í - e ，d k，j，1 ≤ d 0 （4）
ï ï - 1 Queen 相邻定义邻近关系，阶数为 1，即 2 个居住网
ï ï î 0 1 - e 2 ， d k，j，1 > d 0
格之间若有共同边，则 W i，j 为 1，否则为 0； n 为居住
第二步，计算公园的可达性：
网格的数量。
A i = B 1∑ G( d i，j，1 ，d 0 )× R j +
3）洛伦兹曲线和基尼系数法。洛伦兹曲线与基
i ∈{ d i，j，1 ≤ d 0 }
B 2∑ G( d i，j，2 ，d 0 )× R j +
i ∈{ d i，j，2 ≤ d 0 } 尼系数常被用于定量测度社会收入分配的公平程
B 3∑ G( d i，j，3 ，d 0 )× R j （5） 度，由于收入分配的公平程度在本质上与公园绿地
i ∈{ d i，j，3 ≤ d 0 }
式（5）中， A i 为居住点 i的公园绿地可达性，其值 资源分配较为相似，这 2种方法也在环境公平领域被
越大表示公园绿地的可达性越高。 广泛应用 ［28］ 。本研究基于洛伦兹曲线与基尼系数构
2）双变量局部空间自相关。基于两步移动搜索 建公园绿地供需平衡模型，对中心城区整体公园绿
法计算出公园绿地可达性本质是基于居民和公园之 地分配的供需平衡进行测度，其中洛伦兹曲线定义
间的距离以及公园的覆盖范围计算出的加权人均公 为居住点 k 的累计人口比例 P k 与其累计获取的绿地
园绿地面积。因此，本研究将在 2种时间阈值下计算 资源比例 C k 之间的函数关系，基尼系数的计算公式
出的公园绿地可达性作为居住网格单元的绿地供 如下：
n
给，居住网格单元的居住人数作为居住点的绿地需 G = 1 - ∑ k = 1 (P k - P k - 1) (C k + C k - 1 ) （7）
k
求，基于双变量空间自相关 ［26］ 探讨居住网格单元的 ∑ i = 1 A i r i
C k = （8）
n
绿地供给与绿地需求之间的匹配程度，得到中心城 ∑ i = 1 A i r i
区综合公园的供需平衡关系，计算公式如下： 式（7）~（8）中： G为中心城区的公园绿地供需平
-- - 衡指数， n为中心城区内的居住点数， k为居住点按照
Y i，k - Y k n Y j，l - Y l
I i，k，l = × ∑ j = 1 W i，j × （6）
σ k σ l 其公园绿地可达性按升序排序后的序号； A i 为居住