[发明专利]一种基于双层规划的共享自行车区域投放规划方法

权利要求书

1.一种基于双层规划的共享自行车区域投放规划方法，其特征在于，包括以下步骤：

1)获取拟提供服务的共享自行车系统的规划基础资料，包括备选共享自行车区域投放点资料、交通小区资料、交通小区和备选共享自行车区域投放点距离资料、交通小区之间的共享自行车出行需求资料、建设运营成本资料、路网相关资料；

2)根据获取得到的基础资料，构建各条出行路径的广义出行成本函数，其中广义出行成本由步行时间和骑行时间两部分构成；

构建各条出行路径的广义出行成本函数，包含如下步骤：

2.1)建立共享自行车出行路径：共享自行车出行路径是指用户使用共享自行车从起点到终点的单向出行过程假设有k条路径，每条路径共包括5个阶段：a、从起点交通小区r步行至借车点m；b、在借车点m借车；c、使用共享自行车从借车点m骑行至另一还车点n；d、在还车点n还车；e、从还车点n步行至终点交通小区s；

2.2)计算步行时间：步行时间是指用户从起点步行至借车点m或从还车点n步行至终点所用的时间；假设交通小区出行需求均匀分布，因此共享自行车出行的起点和终点能够看成起讫点交通小区的形心；步行时间计算公式为：

式中，分别为起点交通小区r的形心至借车点m和还车点n至终点交通小区s的形心的步行时间；s_rm、s_ns分别为起点交通小区r的形心至借车点m和还车点n至终点交通小区s的形心的距离；v_w是用户平均步行速度；

2.3)计算骑行时间：骑行时间是指用户从借车点m骑行共享自行车到还车点n的时间，骑行时间与借还点距离以及路段流量有关，骑行时间计算公式为：

式中，为借车点m骑行共享自行车到还车点n的骑行时间；c为单位租赁费用；τ为旅行时间价值；s_mn是借车点m和还车点n之间的距离；v_f是在自由流状态下借车点m和还车点n之间的共享自行车的速度；k_j是借车点m和还车点n之间的共享自行车阻塞密度；ξ为速度折减系数，若设置了共享自行车专用道ξ＝1，否则0＜ξ＜1；q_mn为m和n之间的共享自行车流量，计算公式为为OD对(r,s)之间的出行路径k上的共享自行车出行量；

2.4)计算各出行路径的广义出行成本函数：用户在OD对(r,s)间第k条路径使用共享自行车出行的成本包括步行时间成本和骑行时间成本，计算公式为：

式中，为交通小区r到交通小区s的广义出行成本；分别为起点交通小区r的形心至借车点m和还车点n至终点交通小区s的形心的步行时间；为借车点m骑行共享自行车到还车点n的骑行时间；

3)分别以共享自行车网络系统最优和用户均衡为上层模型和下层模型，建立基于双层规划的共享自行车区域投放规划模型，包含如下步骤：

3.1)对模型的建立提出如下假设：

①模型中的投放点为区域投放的虚拟投放点，该点的投放量是这一区域投放量之和；

②共享自行车流是单向的，即总是从借车点到还车点，任意投放点都具备借车和还车功能，视其在出行路径中的位置而定；

③各交通小区的需求或吸引发生点是该小区共享自行车发生量和吸引量最集中的点，但是该点并不一定是小区的形心，不考虑小区内的交通出行，各小区的需求量是根据数据获得的需求量，因此建立的模型是一个确定需求条件下的静态模型；

3.2)以系统用户出行成本最小为目标函数建立上层模型，需要满足的约束条件为：①共享自行车系统建设资金不超过投资资金上限；②共享自行车系统日常运营管理费用不超过系统运营资金上限；③共享自行车区域投放点的数量不超过相应范围；

Subject to:

式中，q_rm为起点r至借车点m借车的客流量；q_ns为还车点n至终点s借车的客流量；q_mn为m和n之间的共享自行车流量，计算公式为为OD对(r,s)之间的出行路径k上的共享自行车出行量；分别为起点交通小区r的形心至借车点m和还车点n至终点交通小区s的形心的步行时间；为借车点m骑行共享自行车到还车点n的骑行时间；b_x为投放的共享自行车数量；y_x为决策变量，当投放点被选中时y_x＝1，否则y_x＝0；c_b为共享自行车单价；C_f为共享自行车系统建设的资金上限；c_o为投放点每小时的运营费用；C_v为系统运营资金的上限；为投放点数量下限，取值不小于规划区域内交通小区的个数；为投放点数量上限；

3.3)下层模型为共享自行车系统用户均衡模型，需要满足的约束条件为：①OD对(r,s)之间各路径上的共享自行车流量之和等于该OD对之间的共享自行车需求量；②路径上的共享自行车流量为非负值；

Subject to:

式中，q_mn为m和n之间的共享自行车流量；为借车点m骑行共享自行车到还车点n的骑行时间；q_rs为OD对(r,s)之间的共享自行车出行量；为OD对(r,s)之间的出行路径k上的共享自行车出行量；

4)对共享自行车区域投放规划模型进行求解，得到各区域投放点的位置以及共享自行车配备方案；其中，求解共享自行车区域投放规划模型的方法为启发式算法，选择结合遗传算法的混合粒子群算法进行求解，求解过程如下：

①将下层模型转化为其一阶条件，公式如下；用下层模型的一阶条件代替下层模型，得到一个与双层规划模型等价的单层模型；

式中，M、ε均为正数，M＞＞ε；为表征流量大小的变量，当否则为OD对(r,s)之间的出行路径k上的共享自行车出行量；为交通小区r到交通小区s的广义出行成本；u_rs为OD对(r,s)之间的最短出行时间；

②令迭代次数t＝0，利用随机函数生成初始值和初始速度以及决策变量y_x的初始值当投放点被选中时y_x＝1，否则y_x＝0；

③采用粒子群算法更新的位置和速度采用遗传算法中的选择、交叉、变异策略更新y_x的位置

④根据算法的收敛准则判断结果是否满足迭代停止条件，若满足，则得到问题的最终优化结果；若不满足则返回到步骤③继续求解；

⑤根据模型的求解结果，得到共享自行车区域投放点的位置以及需要投放的共享自行车数量。

2.根据权利要求1所述的一种基于双层规划的共享自行车区域投放规划方法，其特征在于，在步骤1)中，所述备选共享自行车区域投放点资料包括拟提供共享自行车的服务范围、服务范围内的共享自行车客流产生情况和共享自行车客流产生点位置，以此确定共享自行车区域投放备选投放点位置，共享自行车客流产生情况以及产生点位置通过交通调查获得；所述交通小区资料包括拟提供共享自行车的服务范围内交通小区的划分情况以及交通小区内共享自行车投放区域的用地性质；所述交通小区和备选共享自行车区域投放点距离资料包括交通小区形心到共享自行车区域投放点的距离以及各共享自行车区域投放点之间的距离，从共享自行车系统规划方案获得；所述交通小区之间的共享自行车出行需求资料是指各交通小区时间的共享自行车出行OD量数据，从历史数据或者交通调查获得；所述建设运营成本资料包括共享自行车单价、运营维护费用、租赁费用和旅行时间价值，通过市场调查和共享自行车系统规划方案获得；所述路网相关资料包括拟提供共享自行车服务范围内的路径情况、道路阻抗函数、共享自行车骑行速度和用户步行速度，从共享自行车系统规划方案获得。