CN108877201B - 一种基于交通指挥指数的警力资源优化方法 - Google Patents

一种基于交通指挥指数的警力资源优化方法 Download PDF

Info

Publication number
CN108877201B
CN108877201B CN201810878089.XA CN201810878089A CN108877201B CN 108877201 B CN108877201 B CN 108877201B CN 201810878089 A CN201810878089 A CN 201810878089A CN 108877201 B CN108877201 B CN 108877201B
Authority
CN
China
Prior art keywords
traffic
event
police
representing
events
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
CN201810878089.XA
Other languages
English (en)
Other versions
CN108877201A (zh
Inventor
缪明月
周芳
李强
张军
李国军
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Capital International Airport Public Security Bureau
CAPITAL UNIVERSITY OF ECONOMICS AND BUSINESS
Original Assignee
Beijing Capital International Airport Public Security Bureau
CAPITAL UNIVERSITY OF ECONOMICS AND BUSINESS
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Beijing Capital International Airport Public Security Bureau, CAPITAL UNIVERSITY OF ECONOMICS AND BUSINESS filed Critical Beijing Capital International Airport Public Security Bureau
Priority to CN201810878089.XA priority Critical patent/CN108877201B/zh
Publication of CN108877201A publication Critical patent/CN108877201A/zh
Application granted granted Critical
Publication of CN108877201B publication Critical patent/CN108877201B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G08SIGNALLING
    • G08GTRAFFIC CONTROL SYSTEMS
    • G08G1/00Traffic control systems for road vehicles
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/04Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/40Business processes related to the transportation industry

Landscapes

  • Business, Economics & Management (AREA)
  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Economics (AREA)
  • General Physics & Mathematics (AREA)
  • Human Resources & Organizations (AREA)
  • Strategic Management (AREA)
  • Theoretical Computer Science (AREA)
  • Marketing (AREA)
  • Tourism & Hospitality (AREA)
  • General Business, Economics & Management (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Game Theory and Decision Science (AREA)
  • Development Economics (AREA)
  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明涉及一种基于交通指挥指数的警力资源优化方法,该方法通过量化参数的方式求解出表征警情状况的交通指挥指数值,并将排队论与泊松分布结合应用于交通警察优化调度中,使交通警察数量充足,以确保快速响应在指定到达时间内发生的各类交通事件,通过使用泊松分布来计算概率大于0.95时发生交通事件的次数,并将求得的值归入排队模型不等式中,从而进一步导出最小巡逻交警组数,本发明首次将故障树分析与交通指挥指数计算模型相结合,通过修改事件发生率将排队模型与泊松分布相结合,考虑交通运行规律,提出了交通指挥指数计算模型可广泛应用于警务巡逻和警备人员的优化调度中。

Description

一种基于交通指挥指数的警力资源优化方法
技术领域
本发明涉及交通领域,更具体的说是涉及一种基于交通指挥指数的警力资源优化方法。
背景技术
随着中国经济的突飞猛进以及现代交通的日益发达,道路上的警情数量不断增多,警情的发生给广大群众以及交警指挥部门造成极大的困扰。因此,为了处理各种交通警情,充足的警力以及对警力资源的合理部署是不可或缺的。现有的交通警力资源部署主要是交警指挥部门依靠专家经验或者是所管辖区内的道路状况来实现的,这种方法掺杂主观性因素较多,而交通事件的发生随机性比较强,因此可能会导致一些路段警力不足,不能对发生的警情及时处理。
为了实现传统警务模式向大数据背景下的现代化勤务模式的转型升级,应用最简单的描述方法来精准描述交通的运行态势,迫切需要研究运用交通指挥指数。交通指挥指数由交通拥堵指数、交通事件指数、交通执法指数构成。运用交通指挥指数可以实现上勤方案的动态调整,但是由于交通指挥指数建立在多指标评价体系之上,各指标性质不同,目前将交通指挥指数运用于警力资源部署的方案并未取得很好的应用效果。
因此,如何提供一种警力资源分配合理、部署过程可量化、科学直观的基于交通指挥指数的警力资源部署方法是本领域技术人员亟需解决的问题。
发明内容
本发明所要解决的技术问题是提供一种警力资源分配合理、部署过程可量化、科学直观的基于交通指挥指数的警力资源部署方法。
本发明解决上述技术问题的技术方案如下:一种基于交通指挥指数的警力资源优化方法,该方法包括以下步骤:
将目标区域自由划分为多个基本反应格,并按照24小时制划分为24个时间段;
构建FAT树图,并计算FAT树图中顶事件发生的概率,计算公式为:
pT=pa+p1c+p1e=pa+pcprc+pepre
Figure GDA0002506339680000021
Figure GDA0002506339680000022
Figure GDA0002506339680000023
其中,a、c、e分别表示交通事件、交通拥堵和交通违章事件;pT表示FAT树图中顶事件发生的概率,即交通警察出动的概率;p1c,p1e是FAT树图中中间事件发生的概率;pa、pc、pe分别表示交通事件、交通拥堵和交通违章事件发生的概率;prc、pre是FAT树图中基础事件发生的概率,即对应事件引发派遣警察的概率;ga,gc,ge分别表示交通事件、交通拥堵和交通违章事件的关键重要度,其作为对应标准化数据的系数;故障树分析是一种描述事件因果关系的有方向的“树”,是系统安全工程中的重要的分析方法之一,它能对各种系统的危险性进行识别评价,既适用于定性分析,又能进行定量分析,具有简明、形象化的特点,体现了以系统工程方法研究安全问题的系统性、准确性和预测性;
通过定期交通监控统计每天中第i基本反应格第j时段内交通事件、交通拥堵和交通违章事件的处理人员的平均数量,并对统计结果进行数据归一化处理,数据归一化处理的公式为:
Figure GDA0002506339680000031
Figure GDA0002506339680000032
Figure GDA0002506339680000033
其中,Naij、Ncij和Neij分别表示每天中第i基本反应格第j时段内交通事件、交通拥堵和交通违章事件的处理人员的平均数量;
Figure GDA0002506339680000034
Figure GDA0002506339680000035
分别表示每天中第i基本反应格第j时段内交通事件、交通拥堵和交通违章事件的处理人员的平均数量标准值;Namin和Namax分别表示Naij的最小值和最大值;Ncmin和Ncmax分别表示Ncij的最小值和最大值;Nemin和Nemax分别表示Neij的最小值和最大值;f(a)为min-max标准化函数;a为式中的Naij、Ncij或Neij
计算每天中第i基本反应格第j时段内的交通指挥指数值,计算公式为;
Figure GDA0002506339680000036
Figure GDA0002506339680000037
其中,go表示除交通事件、交通拥堵和交通违章事件以外的其他事件关键重要度;N* oij表示每天中第i基本反应格第j时段内其他事件的处理人员的平均数量标准值;wij表示在每天中第i基本反应格第j时段内处理交通事件所需时间;w* ij表示在每天中第i基本反应格第j时段内处理交通事件所需时间的标准值;Iij表示第i基本反应格第j时段内交通指挥指数值;
利用地理信息分析软件在目标区域的地图中标记出多个所述基本反应格,并分别对交通事件、交通拥堵和交通违章事件进行位置标定,根据第i基本反应格第j时段内交通指挥指数值Iij计算质心坐标,确定目标区域内最佳巡更点的位置;
利用排队论计算事件发生后当事人在接受事件处理前的平均等待时间Wq,计算公式为:
Figure GDA0002506339680000041
其中,c表示巡逻交警组数,λ表示一个时段交通事件的发生次数,μ表示巡逻小组每小时处理的事件数,ρ=λ/(cμ)表示事件处理强度,Lq表示等候处理的排队队长,k∈[0,c-1],P0表示没有事件正在被处理的概率;
设置当事人的最短平均等待时间h,即应满足Wq≤h,求得巡逻交警组数,比如设置当事人的平均等待时间不超过15min,即应满足Wq≤0.25;
用泊松分布计算在一个时段内发生x次交通事件的累计概率,并计算一个时段内累计概率大于0.95时发生交通事件的平均次数;
将求得的一个时段内累计概率大于0.95时发生交通事件的次数作为新的λ的值,代入事件发生后事件当事人在接受事件处理前的平均等待时间的计算公式中,进一步求得最小巡逻交警组数;
根据处于巡逻岗和备勤岗的交警的实际工作制度,在求得的最小巡逻交警组数的基础上建立用于警力调度的整数规划模型,进一步求得道路巡逻交警的人数和办公室待命交警的人数;
根据交通警察的实际工作数据并利用建立的用于警力调度的整数规划模型,完成警务人员中道路巡逻交警与办公室待命交警的调动工作。
本发明的有益效果是:通过科学量化参数的方式求解出表征警情状况的交通指挥指数值,并将排队论与泊松分布结合应用于交通警察优化调度中,使交通警察数量充足,以确保快速响应在指定到达时间内发生的事件,通过使用泊松分布来计算概率大于0.95时发生交通事件的次数,并将求得的值归入排队模型不等式中,从而进一步导出最小巡逻交警组数,本发明首次将排队模型与泊松分布相结合,考虑交通运行规律,本发明首次将故障树分析与交通指挥指数计算模型相结合,将排队模型与泊松分布相结合,考虑交通运行规律,提出了交通指挥指数计算模型,可广泛应用于警务巡逻和警备人员的优化调度中。
在上述技术方案的基础上,本发明还可以做如下改进。
进一步,对统计结果进行数据归一化处理的条件为:
N(i+1)(j+1)>Nij,N* (i+1)(j+1)>N*ij
f(kNij)=f(Nij)
Figure GDA0002506339680000051
其中,f(Nij)表示Nij的min-max标准化函数;Nij表示Naij、Ncij或Neij;N* ij表示
Figure GDA0002506339680000052
Figure GDA0002506339680000053
进一步,所述每天中第i基本反应单元第j时段内处理交通事件所需时间wij的计算公式为:
Figure GDA0002506339680000054
Figure GDA0002506339680000055
其中,dij表示从事件发生地到交通驻地的距离;(xij,yij)是事件发生地的坐标位置;(x0,y0)是交通驻地的坐标位置;vij表示处理人员从交通驻地到事件发生地的平均速度。
进一步,所述每天中第i基本反应单元第j时段内处理交通事件所需时间wij的计算公式为:
Figure GDA0002506339680000056
Figure GDA0002506339680000057
其中,dij表示从事件发生地到交通驻地的距离;(xij,yij)是事件发生地的坐标位置;(x0,y0)是交通驻地的坐标位置;vij表示处理人员从交通驻地到事件发生地的平均速度。
进一步,处理人员从交通驻地到事件发生地的平均速度vij的计算公式为:
Figure GDA0002506339680000061
其中,vfij表示夜间处理人员从交通驻地到事件发生地的平均速度,rVfij表示除夜间的其他时段处理人员从交通驻地到事件发生地的平均速度;r表示夜间与其他时段处理人员从交通驻地到事件发生地的平均速度的比率。
进一步,根据第i基本反应单元第j时段内交通指挥指数值Iij计算质心坐标,计算公式为:
Figure GDA0002506339680000062
Figure GDA0002506339680000063
其中,Ihj表示第h个基本反应区域第j小时的交通指挥调度指数,其值等于该区域所有点的交通指挥调度指数之和;xh表示质心坐标中的横坐标,yh表示质心坐标中的纵坐标。
进一步,所述整数规划模型的基本公式如下:
Figure GDA0002506339680000064
其中,n∈Z+,fj表示道路巡逻交警的人数,xj表示办公室待命交警的人数。
进一步,一个时段内发生x次交通事件的累计概率的计算公式为:
Figure GDA0002506339680000071
其中,系数λ表示一个时段内交通事件的平均发生率,x表示一个时段内交通事件发生次数。
进一步,交通警察的实际工作数据包括交警的工作时间、休息时间、不同岗位人员换岗规则。
采用上述进一步方案的有益效果是在根据量化参数后得到的数据部署警力资源的过程中也考虑到了交警在实际工作中的作息时间及不同岗位人员的换岗规则,在保证交警的工作规律不被打乱的前提下合理实施本方法。
进一步,所述地理信息分析软件为ArcGIS地理信息系统软件。ArcGIS地理信息系统软件为用户提供一个可伸缩的、全面的GIS平台,该软件包含了许多的可编程组件,从细粒度的对象(例如单个的几何对象)到粗粒度的对象(例如与现有ArcMap文档交互的地图对象)涉及面极广,这些对象为开发者集成了全面的GIS功能,用户应用功能全面的GIS工作站来编辑地理数据集,建立数据编辑和质量控制的工作流,创建地图和分析模型并将这些工作和方法记录成文档,方便管理人员对目标区域巡更点的查看和管理。
附图说明
图1为本发明一种基于交通指挥指数的警力资源优化方法的方法流程图;
图2为本发明一种基于交通指挥指数的警力资源优化方法中FAT树图的结构示意图;
图3为本发明一种基于交通指挥指数的警力资源优化方法中利用泊松分布反映在等待时间为15分钟以内时不同的概率下巡逻交警数在不同时段内变化情况的折线图;
图4为本发明一种基于交通指挥指数的警力资源优化方法中利用排队论反映不同的等待时间下巡逻交警数在各时间段内变化情况的折线图;
图5为本发明实施例中反映首都国际机场路附近平均1天内各时段的交通事件起数的折线图;
图6为本发明实施例中工作日北京首都国际机场交通支队管界内拥堵指数的变化趋势图;
图7为本发明实施例中周末北京首都国际机场交通支队管界内拥堵指数的变化趋势图。
具体实施方式
以下结合附图对本发明的原理和特征进行描述,所举实例只用于解释本发明,并非用于限定本发明的范围。
如图1所示,一种基于交通指挥指数的警力资源优化方法,该方法包括以下步骤:
S1:将目标区域自由划分为多个基本反应格,并按照24小时制划分为24个时间段;
S2:参见附图2,构建FAT树图,并计算FAT树图中顶事件发生的概率,计算公式为:
pT=pa+p1c+p1e=pa+pcprc+pepre
Figure GDA0002506339680000081
Figure GDA0002506339680000082
Figure GDA0002506339680000083
其中,a、c、e分别表示交通事件、交通拥堵和交通违章事件;pT表示FAT树图中顶事件发生的概率;p1c,p1e是FAT树图中中间事件发生的概率;pa、pc、pe分别表示交通事件、交通拥堵和交通违章事件发生的概率;prc、pre是FAT树图中基础事件发生的概率;ga,gc,ge分别表示交通事件、交通拥堵和交通违章事件的关键重要度,其作为对应标准化数据的系数;
S3:通过定期交通监控统计每天中第i基本反应格第j时段内交通事件、交通拥堵和交通违章事件的处理人员的平均数量,并对统计结果进行数据归一化处理,数据归一化处理的公式为:
Figure GDA0002506339680000091
Figure GDA0002506339680000092
Figure GDA0002506339680000093
其中,Naij、Ncij和Neij分别表示每天中第i基本反应格第j时段内交通事件、交通拥堵和交通违章事件的处理人员的平均数量;
Figure GDA0002506339680000094
Figure GDA0002506339680000095
分别表示每天中第i基本反应格第j时段内交通事件、交通拥堵和交通违章事件的处理人员的平均数量标准值;Namin和Namax分别表示Naij的最小值和最大值;Ncmin和Ncmax分别表示Ncij的最小值和最大值;Nemin和Nemax分别表示Neij的最小值和最大值;f(a)为min-max标准化函数;a为式中的Naij、Ncij或Neij
S4:计算每天中第i基本反应格第j时段内的交通指挥指数值,计算公式为:
Figure GDA0002506339680000096
Figure GDA0002506339680000097
其中,go表示除交通事件、交通拥堵和交通违章事件以外的其他事件关键重要度;N* oij表示每天中第i基本反应格第j时段内其他事件的处理人员的平均数量标准值;wij表示在每天中第i基本反应格第j时段内处理交通事件所需时间;w* ij表示在每天中第i基本反应格第j时段内处理交通事件所需时间的标准值;Iij表示第i基本反应格第j时段内交通指挥指数值;
利用地理信息分析软件在目标区域的地图中标记出多个所述基本反应格,并分别对交通事件、交通拥堵和交通违章事件进行位置标定,根据第i基本反应格第j时段内交通指挥指数值Iij计算质心坐标,确定目标区域内最佳巡更点的位置;
S5:利用排队论计算事件发生后当事人在接受事件处理前的平均等待时间Wq,计算公式为:
Figure GDA0002506339680000101
其中,c表示巡逻交警组数,λ表示一个时段交通事件的发生次数,μ表示巡逻小组每小时处理的事件数,ρ=λ/(cμ)表示事件处理强度,Lq表示等候处理的排队队长,k∈[0,c-1],P0表示没有事件正在被处理的概率;
S6:设置当事人的最短平均等待时间h,即应满足Wq≤h,求得巡逻交警组数;
S7:用泊松分布计算在一个时段内发生x次交通事件的累计概率,并计算一个时段内累计概率大于0.95时发生交通事件的平均次数;
S8:将求得的一个时段内累计概率大于0.95时发生交通事件的次数作为新的λ的值,代入事件发生后事件当事人在接受事件处理前的平均等待时间的计算公式中,进一步求得最小巡逻交警组数;
S9:根据处于巡逻岗和备勤岗的交警的实际工作制度,在求得的最小巡逻交警组数的基础上建立用于警力调度的整数规划模型,进一步求得道路巡逻交警的人数和办公室待命交警的人数;
S10:根据交通警察的实际工作数据并利用建立的用于警力调度的整数规划模型,完成警务人员中道路巡逻交警与办公室待命交警的调动工作。
本实施中,对统计结果进行数据归一化处理的条件为:
N(i+1)(j+1)>Nij,N* (i+1)(j+1)>N*ij
f(kNij)=f(Nij)
Figure GDA0002506339680000111
其中,f(Nij)表示Nij的min-max标准化函数;Nij表示Naij、Ncij或Neij;N* ij表示
Figure GDA0002506339680000112
Figure GDA0002506339680000113
具体地,每天中第i基本反应单元第j时段内处理交通事件所需时间wij的计算公式为:
Figure GDA0002506339680000114
Figure GDA0002506339680000115
其中,dij表示从事件发生地到交通驻地的距离;(xij,yij)是事件发生地的坐标位置;(x0,y0)是交通驻地的坐标位置;vij表示处理人员从交通驻地到事件发生地的平均速度。
具体地,处理人员从交通驻地到事件发生地的平均速度vij的计算公式为:
Figure GDA0002506339680000116
其中,vfij表示夜间处理人员从交通驻地到事件发生地的平均速度,rVfij表示除夜间的其他时段处理人员从交通驻地到事件发生地的平均速度;r表示夜间与其他时段处理人员从交通驻地到事件发生地的平均速度的比率。
具体地,根据第i基本反应单元第j时段内交通指挥指数值Iij计算质心坐标,计算公式为:
Figure GDA0002506339680000121
Figure GDA0002506339680000122
其中,Ihj表示第h个基本反应区域第j小时的交通指挥调度指数,其值等于该区域所有点的交通指挥调度指数之和;xh表示质心坐标中的横坐标,yh表示质心坐标中的纵坐标。
具体地,整数规划模型的基本公式如下:
Figure GDA0002506339680000123
其中,n∈Z+,fj表示道路巡逻交警的人数,xj表示办公室待命交警的人数。
上述实施例中,一个时段内发生x次交通事件的累计概率的计算公式为:
Figure GDA0002506339680000124
其中,系数λ表示一个时段内交通事件的平均发生率,x表示一个时段内交通事件发生次数。
具体地,本实施例中交通警察的实际工作数据包括交警的工作时间、休息时间、不同岗位人员换岗规则。
具体地,地理信息分析软件为ArcGIS地理信息系统软件。ArcGIS地理信息系统软件为用户提供一个可伸缩的、全面的GIS平台,该软件包含了许多的可编程组件,从细粒度的对象(例如单个的几何对象)到粗粒度的对象(例如与现有ArcMap文档交互的地图对象)涉及面极广,这些对象为开发者集成了全面的GIS功能,用户应用功能全面的GIS工作站来编辑地理数据集,建立数据编辑和质量控制的工作流,创建地图和分析模型并将这些工作和方法记录成文档,方便管理人员对目标区域巡更点的查看和管理。
参见附图3和图4,当时间小于8时,泊松分布与排队论中的Wq的差别相对较大,随着时间的增加,差异变得越来越小。对于泊松分布,当累积概率小于或等于0.7时所需的交通警察数量明显大于累积概率小于或等于0.6时的交通警察数量。随着时间的增大,这一点没有任何变化。当然,我们的重点不是泊松分布本身,而是泊松分布和排队模型之间的差异。在道路上巡逻的警察总数和警察待命的办公室至少应该足以应付不超过95%的事件,为确定警察巡逻道路的数量和警备人员的数量,有必要考虑目前警员的总数与执法和事件处理的分配比例。
参见附图5,在经过对首都国际机场路附近415天的现场统计,得到平均1天即24小时内各时段的交通事件起数,并绘制成折线图,通过对上述时限内现场交通事件进行观察分析,发现一般情况下,某小时时段内,有现场执法记录的人员多则4人,少则1人,甚至没有记录,对工作日及周末的实测数据分别进行整理,具体结果参见下表1和表2。
表1工作日统计结果数据表
Figure GDA0002506339680000141
表2周末统计结果数据表
Figure GDA0002506339680000151
根据高德地图对于拥堵指数的定义是vfj/vj,vfj即在第j小时的最高车速,一般被认为是发生在车流量很小的夜间。vj是第j小时的平均车速。北京首都国际机场交通支队管界内拥堵指数在周末与工作日的变化趋势参见附图6和附图7,其中附图6为工作日北京首都国际机场交通支队管界内拥堵指数的变化趋势,图7为周末北京首都国际机场交通支队管界内拥堵指数的变化趋势。
由图6和图7可以看出,在北京首都国际机场存在明显的早晚高峰。考虑到拥堵因素以及便于计算的目的,对于计算备勤警力人数:从早7点至晚7点,以五分钟到达现场为标准,即wq≤0.08。其他时段,以十五分钟到达现场为标准,即wq≤0.25。备勤警力加上巡逻警力可以满足处理95%的各类交通事件,且满足到达现场时间要求。这样,整数规划模型的各个约束参数如下表3所示:
表3整数规划模型的关键参数统计表
Figure GDA0002506339680000161
考虑到基本条件和交通警察工作习惯,对于整数规划模型,模型建立的主要参考交警实际工作数据如下:
(1)警察的工作时间可分为两个时期:日班时间为8:30~17:30,夜班时间为17:30~8:30。
(2)每个交通警察巡逻队在第一天开始只在白班工作,然后在第二天工作24小时(白班和夜班),第三天和第四天休息。
(3)巡逻交警岗位与备勤交警岗位是相对独立的。
(4)备勤警力(主要工作时间都在办公室)每周有两天休息
(5)备勤警力一周须有一晚上值班时间。
无论工作日抑或周末,无论白班或夜班,道路上巡逻的交警数量最少需要4人,可就是最少需要16名巡逻交警,按照每组2人建立8组巡逻交警,分别是Team 1~Team 8。具体分配结果参见下表4:
表4巡逻警察组排班表
Figure GDA0002506339680000171
根据上述数据,整数规划模型建立如下:
Figure GDA0002506339680000172
约束条件s.t为:
xd1+xd4+xd5+xd6+xd7≥10
xd1+xd2+xd5+xd6+xd7≥10
xd1+xd2+xd3+xd6+xd7≥10
xd1+xd2+xd3+xd4+xd7≥10
xd1+xd2+xd3+xd4+xd5≥10
xd6+xd2+xd3+xd4+xd5≥6
xd7+xd2+xd3+xd4+xd5≥6
Figure GDA0002506339680000173
xdn≥0,xdn∈Z
其中,minG是最少的备勤警力数,xdn(n=1...7)分别表示从周一到周日每天配备的备勤警力数。从周一到周日每天配备的备勤警力数的统计结果可参见下表5:
表5一周内备勤警力数统计表
Figure GDA0002506339680000174
通过上表及相关计算,得到备勤警力目标值为38,加上18名巡逻警力,总共需要54名交通警察。备勤警力日常主要在办公室工作,其中32人每周一开始工作,其他6人每周三开始工作,确保备勤警力一周有2天休息。综合来看,一方面警力安排满足了覆盖95%交通事件处理所需警力的需求,另一方面,也达到了有关警察到达现场的时间要求。
本实施例提供的一种基于交通指挥指数的警力资源优化方法,通过量化参数的方式求解出表征警情状况的交通指挥指数值,并将排队论与泊松分布结合应用于交通警察优化调度中,使交通警察数量充足,以确保快速响应在指定到达时间内发生的事件,通过使用泊松分布来改变计算待机警察的最小数目时的平均事件率的值,将其归入排队模型不等式,从而进一步导出警察待命的最小数,本发明首次将排队模型与泊松分布相结合,考虑交通运行规律,可广泛应用于警务巡逻和警备人员的优化调度中。
以上所述仅为本发明的较佳实施例,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (8)

1.一种基于交通指挥指数的警力资源优化方法,其特征在于,包括以下步骤:
将目标区域自由划分为多个基本反应格,并按照24小时制划分为24个时间段;
构建FAT树图,并计算FAT树图中顶事件发生的概率,计算公式为:
pT=pa+p1c+p1e=pa+pcprc+pepre
Figure FDA0002506339670000011
Figure FDA0002506339670000012
Figure FDA0002506339670000013
其中,a、c、e分别表示交通事件、交通拥堵和交通违章事件;pT表示FAT树图中顶事件发生的概率;p1c,p1e是FAT树图中中间事件发生的概率;pa、pc、pe分别表示交通事件、交通拥堵和交通违章事件发生的概率;prc、pre是FAT树图中基础事件发生的概率;ga,gc,ge分别表示交通事件、交通拥堵和交通违章事件的关键重要度,其作为对应标准化数据的系数;
通过定期交通监控统计每天中第i基本反应格第j时段内交通事件、交通拥堵和交通违章事件的处理人员的平均数量,并对统计结果进行数据归一化处理,数据归一化处理的公式为:
Figure FDA0002506339670000014
Figure FDA0002506339670000015
Figure FDA0002506339670000016
其中,Naij、Ncij和Neij分别表示每天中第i基本反应格第j时段内交通事件、交通拥堵和交通违章事件的处理人员的平均数量;
Figure FDA0002506339670000024
Figure FDA0002506339670000025
分别表示每天中第i基本反应格第j时段内交通事件、交通拥堵和交通违章事件的处理人员的平均数量标准值;Namin和Namax分别表示Naij的最小值和最大值;Ncmin和Ncmax分别表示Ncij的最小值和最大值;Nemin和Nemax分别表示Neij的最小值和最大值;f(a)为min-max标准化函数,a为式中的Naij、Ncij或Neij
计算每天中第i基本反应格第j时段内的交通指挥指数值,计算公式为:
Figure FDA0002506339670000021
Figure FDA0002506339670000022
其中,go表示除交通事件、交通拥堵和交通违章事件以外的其他事件的关键重要度;N* oij表示每天中第i基本反应格第j时段内其他事件的处理人员的平均数量标准值;wij表示在每天中第i基本反应格第j时段内处理交通事件所需时间;w* ij表示在每天中第i基本反应格第j时段内处理交通事件所需时间的标准值;Iij表示第i基本反应格第j时段内交通指挥指数值;
利用地理信息分析软件在目标区域的地图中标记出多个所述基本反应格,并分别对交通事件、交通拥堵和交通违章事件进行位置标定,根据第i基本反应格第j时段内交通指挥指数值Iij计算质心坐标,确定目标区域内最佳巡更点的位置;
利用排队论计算事件发生后当事人在接受事件处理前的平均等待时间Wq,计算公式为:
Figure FDA0002506339670000023
其中,c表示巡逻交警组数,λ表示一个时段交通事件的发生次数,μ表示巡逻小组每小时处理的事件数,ρ=λ/(cμ)表示事件处理强度,Lq表示等候处理的排队队长,k∈[0,c-1],P0表示没有事件正在被处理的概率;
设置当事人的最短平均等待时间h,即应满足Wq≤h,求得巡逻交警组数;
用泊松分布计算在一个时段内发生x次交通事件的累计概率,并计算一个时段内累计概率大于0.95时发生交通事件的平均次数;
将求得的一个时段内累计概率大于0.95时发生交通事件的次数作为新的λ的值,代入事件发生后事件当事人在接受事件处理前的平均等待时间的计算公式中,进一步求得最小巡逻交警组数;
根据处于巡逻岗和备勤岗的交警的实际工作制度,在求得的最小巡逻交警组数的基础上建立用于警力调度的整数规划模型,进一步求得道路巡逻交警的人数和办公室待命交警的人数;
根据交通警察的实际工作数据并利用建立的用于警力调度的整数规划模型,完成道路巡逻交警与办公室待命交警的调动工作。
2.根据权利要求1所述一种基于交通指挥指数的警力资源优化方法,其特征在于,对统计结果进行数据归一化处理的条件为:
N(i+1)(j+1)>Nij,N* (i+1)(j+1)>N*ij
f(kNij)=f(Nij)
Figure FDA0002506339670000031
其中,f(Nij)表示Nij的min-max标准化函数;Nij表示Naij、Ncij或Neij;N* ij表示
Figure FDA0002506339670000032
Figure FDA0002506339670000033
3.根据权利要求1所述一种基于交通指挥指数的警力资源优化方法,其特征在于,所述每天中第i基本反应单元第j时段内处理交通事件所需时间wij的计算公式为:
Figure FDA0002506339670000041
Figure FDA0002506339670000042
其中,dij表示从事件发生地到交通驻地的距离;(xij,yij)是事件发生地的坐标位置;(x0,y0)是交通驻地的坐标位置;vij表示处理人员从交通驻地到事件发生地的平均速度。
4.根据权利要求1所述一种基于交通指挥指数的警力资源优化方法,其特征在于,处理人员从交通驻地到事件发生地的平均速度vij的计算公式为:
Figure FDA0002506339670000043
其中,vfij表示夜间处理人员从交通驻地到事件发生地的平均速度,rVfij表示除夜间的其他时段处理人员从交通驻地到事件发生地的平均速度;r表示夜间与其他时段处理人员从交通驻地到事件发生地的平均速度的比率。
5.根据权利要求1所述一种基于交通指挥指数的警力资源优化方法,其特征在于,根据第i基本反应单元第j时段内交通指挥指数值Iij计算质心坐标,计算公式为:
Figure FDA0002506339670000044
Figure FDA0002506339670000045
其中,Ihj表示第h个基本反应区域第j小时的交通指挥调度指数,其值等于该区域所有点的交通指挥调度指数之和;xh表示质心坐标中的横坐标,yh表示质心坐标中的纵坐标。
6.根据权利要求1所述一种基于交通指挥指数的警力资源优化方法,其特征在于,所述整数规划模型的基本公式如下:
Figure FDA0002506339670000051
其中,n∈Z+,fj表示道路巡逻交警的人数,xj表示办公室待命交警的人数。
7.根据权利要求1-6任一项所述一种基于交通指挥指数的警力资源优化方法,其特征在于,交通警察的实际工作数据包括交警的工作时间、休息时间、不同岗位人员换岗规则。
8.根据权利要求1-6任一项所述一种基于交通指挥指数的警力资源优化方法,其特征在于,所述地理信息分析软件为ArcGIS地理信息系统软件。
CN201810878089.XA 2018-08-03 2018-08-03 一种基于交通指挥指数的警力资源优化方法 Expired - Fee Related CN108877201B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810878089.XA CN108877201B (zh) 2018-08-03 2018-08-03 一种基于交通指挥指数的警力资源优化方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810878089.XA CN108877201B (zh) 2018-08-03 2018-08-03 一种基于交通指挥指数的警力资源优化方法

Publications (2)

Publication Number Publication Date
CN108877201A CN108877201A (zh) 2018-11-23
CN108877201B true CN108877201B (zh) 2020-10-27

Family

ID=64307380

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810878089.XA Expired - Fee Related CN108877201B (zh) 2018-08-03 2018-08-03 一种基于交通指挥指数的警力资源优化方法

Country Status (1)

Country Link
CN (1) CN108877201B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111160537B (zh) * 2020-01-03 2022-08-19 南京邮电大学 一种基于ann的路口交通警力资源调度系统

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6973384B2 (en) * 2001-12-06 2005-12-06 Bellsouth Intellectual Property Corporation Automated location-intelligent traffic notification service systems and methods
CN102044011B (zh) * 2009-10-20 2013-07-31 北京交通大学 一种警力资源调度的方法及系统
CN107944727A (zh) * 2017-12-04 2018-04-20 杨川 一种基于历史安全度模型的警力调度融合方法
CN108269000A (zh) * 2017-12-22 2018-07-10 武汉烽火众智数字技术有限责任公司 基于警情大数据时空分析的智能警力部署方法及系统

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
面向交通事件处置的交通警力资源调度方法的研究;马威;《中国优秀硕士学位论文全文数据库 社会科学Ⅰ辑》;20150915(第09期);全文 *

Also Published As

Publication number Publication date
CN108877201A (zh) 2018-11-23

Similar Documents

Publication Publication Date Title
Savas Simulation and cost-effectiveness analysis of New York's emergency ambulance service
Southworth Regional evacuation modeling: A state-of-the-art review
Khattak et al. Modeling traffic incident duration using quantile regression
CN110705747A (zh) 一种基于大数据的智慧公共交通云脑体系
CN108428341A (zh) 一种基于人机融合的交通应急管理调度方法及系统
Salicrú et al. Timetable-based operation in urban transport: Run-time optimisation and improvements in the operating process
US6058370A (en) Method of forecasting ambulance service demand
CN107564270A (zh) 一种智能公交调度运营方法
Zaki et al. A simulation model for the analysis and management of an emergency service system
CN113902237A (zh) 一种现场作业风险管控系统
CN112883310A (zh) 一种基于大数据的多层次综合性社区服务平台及设备
CN114548656B (zh) 一种安全应急管理指挥调度系统
CN110309952B (zh) 一种基于通勤模型的城市就业空间布局优化辅助系统
CN110097226A (zh) 一种资源调度方法、系统、电子设备及可读存储介质
CN108877201B (zh) 一种基于交通指挥指数的警力资源优化方法
CN116362420B (zh) 一种应急救援指挥中心管理方法
CN114529438A (zh) 一种救援现场风险评估与资源管理系统
CN115470940A (zh) 一种基于bim技术的建筑工程运维管理系统
CN111832929B (zh) 一种机场值机的动态调度方法及系统
CN109426945A (zh) 一种面向移动互联网的交通秩序管理方法及系统
Greasley A redesign of a road traffic accident reporting system using business process simulation
CN112085341A (zh) 一种适用于电力生产作业全过程风险管控的主站系统
Kwak et al. Police patrol beat design: Allocation of effort and evaluation of expected performance
Miao et al. Novel Method for Optimizing Traffic Police Manpower by Standby and Patrol
Mabert Planning police communications support

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20201027

Termination date: 20210803