[发明专利]一种基于多层次分区的公共自行车高峰期调度方法

权利要求书

1.一种基于多层次分区的公共自行车高峰期调度方法，其特征在于，包括以下步骤：

1)选定研究范围，采集基础资料，包括规划资料、各租赁点的设施资料、需求资料，其中，所述规划资料包括规划租赁点位置、各租赁点规划配备车辆数和租赁点用地性质类型，该租赁点用地性质类型主要分为交通枢纽、公交中途站、公建商业类、居住类和休闲类五类；所述各租赁点的设施资料指实际锁车桩数；所述需求资料指该租赁点早晚高峰短时段内公共自行车相对准确的借还量，一般记借车量为公共自行车出行发生量，为正值，而还车量为出行吸引量，为负值；

2)确定上层调度区域，分析各租赁点锁车桩数和借还需求关系，通过各点间车辆流动关联度将周围租赁点纳入调度区域，直至上层调度区域内为公共自行车高峰期主要流动范围；

3)在调度区域内进行小区划分，根据各租赁点高峰期运营数据，结合聚类分析和人工调整方法进行小区划分和编号；

4)确定区间调度方案：根据步骤3)的小区划分结果，统计各小区高峰期短时段内总需求量以及小区各租赁点中最早的时间要求，代入调度模型中计算得到区间调度方案；

5)确定区内固定调度方案：根据各小区区内租赁点的分布和道路条件，利用最短路模型进行区内固定调度路线的确定，并结合小区内各租赁点的短时段需求量确定区内调度方案；

6)结合区间和区内调度方案制定高峰期整体方案。

2.根据权利要求1所述的一种基于多层次分区的公共自行车高峰期调度方法，其特征在于：在步骤2)中，分析确定上层调度区域的主要步骤包括：

2.1)判断各租赁点的高峰期每10分钟需求量和其锁车桩数量的差异，若差值超过锁车桩数量的一半，则将其作为供不应求的调度需求点；

2.2)通过公共自行车借还刷卡数据或预约数据，得到调度需求点主要自行车流向租赁点或自行车吸引租赁点，将其纳入上层调度区域内；

2.3)判断调度区域内租赁点互为主要吸引、发生租赁点，即流入和流出调度区域的公共自行车较少时，停止纳入新租赁点；

2.4)判断调度区域内租赁点个数，较多则剔除部分租借关联度较低的租赁点，最终确定调度区域范围；

在步骤3)中，小区划分的主要步骤包括：

3.1)输入步骤2)得到的调度区域内各租赁点高峰期每10分钟的需求数据，利用SPSS软件进行聚类分析，将高峰期借还特性相近的租赁点归为一类；

3.2)根据聚类结果，结合设计人员经验判断，将聚类相同且距离不大于300米的公共自行车租赁点划分为一个小区；

3.3)判断小区内租赁点个数N，若N>6，则剔除部分距离较远的租赁点另划分为一区，否则确定小区范围及内部租赁点，对小区进行编号a∈{a₁,a₂,...,a_m}，同时对小区内部租赁点进行标号b_a∈{b_a1,b_a2,...,b_an}；

在步骤4)中，确定区间调度方案的主要步骤包括：

4.1)统计各小区高峰期每10分钟的需求数据，即小区内所有租赁点每10分钟需求数据之和，将此数据作为小区需求数据；

4.2)统计各小区内租赁点的最早的时间窗要求，将其作为相应小区的时间窗；

4.3)采取设置虚拟小区的形式，将分区按相同顺序多次编号，按高峰期2小时计算，每10分钟统计一次需求量，最多来回调度12次，则以小区个数M为一个周期，循环12个周期编号，即a小区及其虚拟点编号集合为A＝{a,a+M,a+2M,...,a+11M}；

4.4)绘制高峰期期间所有小区的各时段时间窗需求表，其中应包括所有小区编号、形心坐标、锁车桩数、需求量和时间区间；

4.5)根据调度的实际情况，考虑运输车辆的固定成本、行驶成本以及客户的满意度，建立目标函数：

minZ=C0Σj=1nΣk=1mx0jk+C1Σi=0nΣj=0nΣk=1mdijxijk+C2P[Σi=1nΣk=1mmax(ti-li,0)+Σi=1nΣk=1mmax(ei-ti,0)]]]>

其中，表达式等号右边第一部分为固定成本，其取值与人工费用、车辆折旧费用、车辆其它固定费用有关；第二部分为行驶成本，其取值与车辆单位油耗、租赁点间直线距离、城市道路直线系数有关；第三部分为惩罚成本，其取值与运输车辆到达服务点的实际时间、租赁点规定的服务时间、乘客等待时间有关；

4.6)根据既定目标函数以及实际调度情况，确定模型的约束条件如下：

①运输车辆必须从中心车场出发

②运输车辆必须回到中心车场

③运输车辆经过某一小区时，其调度量加上运输车辆原有的载运量不小于0且不大于运输车辆的最大载运能力0≤ω_ijk-q_j≤Q；其中，ω_ijk≥0、i,j∈Nk∈{1,2,...,m}；

④进入某小区的车辆，必须要从该小区离开

⑤每个小区只有一辆车服务

⑥任一小区的需求量不大于调度车辆的容量0≤|q_i|y_ik≤Q；

⑦到达小区j的时刻t_j等于到达上一个小区i并开始服务的时刻与该点服务时间与两点间行程时间t_ij之和t_j＝t_i+t_ij+t_uii,j＝0,1,...,n，t_ij≠0；

⑧到达各小区的上一时间窗的时间必须小于到达该小区下一时间窗的时间t_i＜t_i+M；

上式中：

N：小区编号{0,1,2，…,n}，0代表中心车场；

M：实际小区个数；

K：所需运输车辆数{1,2，…,m}；

d_ij：小区i到小区j的直接距离，假设d_ij＝d_ji；

q_i：小区i的调度量，调出为负，调入为正；

ω_ijk：运输车辆k从租赁点i到租赁点j车上的载运量；

l_i：小区i的最晚服务时间；

e_i：小区i的最早服务时间；

t_i：运输车辆到达小区i的时间；

t_j：运输车辆到达小区j的时间；

t_ij：运输车辆从小区i到达小区j的时间；

t_ui：运输车辆在小区i的装卸时间；

P：时间惩罚系数，与乘客的等待时间成正相关关系；

C₀：运输车辆的固定成本；

C₁：运输车辆单位里程的成本费用；

C₂：等待时间的成本费用；

4.7)采用遗传算法对模型进行求解，主要步骤如下：

①染色体编码

遗传算法在进行搜索前，先将解空间的解数据表示成遗传空间的基因串结构数据，这种从表现型到基因型的映射称为编码；适合本次研究调度模型的遗传编码方式是符号编码，具体表示为g＝(w₁,w₂,...,w_N)，其中w表示每个租赁点被哪辆车服务，w_i＝(m-num,s-num)表示第i个租赁点接受车辆编号为m的运输车辆的配送服务；

②生成初始种群

系统通过染色体的编码方式来随机生成一系列染色体，作为最初的初始种群；初始种群的数量视调度区域的规模大小来确定，其取值在50-200之间；

③计算适应度值

本调度模型中，求的是目标函数的最小值，因此将函数值的倒数作为个体的适应度值，函数值越小的个体，适应度值越大，个体就越优；

f_i＝1/Z

④选择操作

选择操中从旧群体中以一定概率选择优良个体组成新的种群，以繁殖得到下一代个体，个体被选中的概率跟适应度值有关，个体适应度值越高，被选中的概率越大；本次研究选择轮盘赌法，即基于适应度比例的选择策略，个体ⁱ被选中的概率为其中，f_i为个体i的适应值；N为种群个体数目；

⑤交叉操作

从种群中随机选择两个个体，通过两个染色体的交换组合，把父串的优秀特征遗传给子串，从而得到新的优秀个体；

⑥变异操作

从种群中随机选取一个个体，选择个体中的一点进行变异以产生更优秀的个体；

⑦判断是否满足终止条件；停止规则一般为最大进化代数；或算法在连续进化一定代数后，解的适应度没有明显改进；

4.8)整理求解得到的调度顺序和各租赁点应调度的车辆数，制定小区间调度路线方案；

在步骤5)中，确定区内固定调度路线的主要步骤包括：

5.1)采集小区内各租赁点周边道路高峰期交通量数据，以及各租赁点间的连接道路长度；分别对两类数据进行分级赋值，数据最大值赋5，最小值赋1，中间数值根据比例关系赋值，保留两位有效数字；

5.2)将各租赁点间每条连接道路的交通量等级值和道路长度等级值求和，作为该条连接道路的阻抗

式中：Q_ijk——小区内租赁点i至租赁点j间第k条连接道路的交通量等级值；

L_ijk——小区内租赁点i至租赁点j间第k条连接道路的长度等级值；

5.3)比较小区内租赁点i至租赁点j所有连接道路的阻抗值，取最小者作为租赁点i和租赁点j间的路径，其阻抗

5.4)为确保区内服务点有且仅有一次均被服务到，建立最优路径模型，代入道路阻抗值求得小区内固定最优路径：

s.t.x_ij+x_ji≤1

Σi=1MXij=1]]>

Σj=1MXij=1]]>

i≠j

式中：Ω——小区内调度总阻抗值；

5.5)根据步骤5.4)得到的最优调度顺序，结合步骤5.3)中最小阻抗值道路，确定区内调度路线；

5.6)结合各租赁点的需求数据制定区内调度方案。

3.根据权利要求1所述的一种基于多层次分区的公共自行车高峰期调度方法，其特征在于：在步骤6)中，将区内调度方案衔接于区间调度之后，形成从小区到租赁点的完整调度方案。