[发明专利]一种中观层面事故多发道路判别的方法

权利要求书

1.一种中观层面事故多发道路判别的方法，其特征在于，步骤如下：

步骤1：获取道路上所有路段及交叉口的几何设计、交通运行及事故数据；利用街景地图得到道路几何设计数据，基于线圈检测设备获取交通流量数据，根据浮动车数据提取路段的运行速度，根据交通事故的严重程度将事故分为物损事故、伤亡事故；

步骤2：沿道路方向，将道路划分为中观研究单元；根据横断面设计参数、道路运行速度、道路范围，将道路划分为中观研究单元，即中观单元内的横断面、运行特征差异性较小；

步骤3：与道路相交方向，两侧均选用350m作为影响道路运行状况的路网范围，采用介度中心度量化路网的集聚性，将路网划分为方格，不规则方格，混合型三种路网形态；

步骤4：构建中观安全分析样本数据集；提取各中观单元的道路几何特征变量、交通运行特征变量及总事故数、物损及伤亡事故数，构建交通安全分析的样本数据集；

步骤5：建立随机效应负二项安全分析模型；

步骤5.1：针对总事故数构建随机效应负二项交通安全模型；假设路段交通事故数服从负二项分布，y_ij～Negbin(θ_ij，r)，模型方程为y_ii为主干道i第j个中观单元的观测事故数，θ_ij是y_ij的期望值，离散系数r服从Gamma分布(0.001，0.001)，X_ii为自变量，β为估计系数，为主干道随机效应，服从正态分布(0，1/a)，其中a为精度参数，服从Gamma分布(0.001，0.001)；

步骤5.2：针对物损及伤亡事故数构建二元因变量随机效应负二项交通安全模型；假设物损及伤亡事故服从负二项分布，y_ij1～Negbin(θ_ij1，r₁)，y_ij2～Negbin(θ_ij2，r₂)，模型方程为其中y_ij1某中观单元的物损事故数，y_ij2为中观单元的伤亡事故数，θ_ij1是y_ij1的期望值，θ_ij2是y_ij2的期望值；是随机效应项；u_ij1和u_ij2是误差项，服从正态分布u_ij1～N(0，1/τ)，τ是精确度参数，τ～gamma(0.001，0.001)；X_ij1和X_ij2是自变量；β和δ_m为待估计系数；

使用全贝叶斯方法估计交通安全模型参数；首先为参数设置一个特定的先验分布，再与观测数据结合得出后验分布，并通过马尔科夫链蒙特卡洛法完成参数估计；理论框架为其中y为已发生的事故数，θ为事故数期望，L(y|θ)为似然函数，π(θ)为先验分布，π(θ|y)为在给定y条件下θ的后验分布，即将要发生的事故数期望，∫L(y|θ)π(θ)dθ为观测数据的边缘概率分布；

步骤6：计算安全可提高空间；安全可提高空间为贝叶斯估计事故数与类似地点平均事故期望的差值，即其中，为地点ij事故类型k的安全可提高空间，k＝0表示总事故，k＝1表示物损事故，k＝2表示伤亡事故；表示地点ij的事故超出了类似地点，并且由该地点相关特征所引起，可采取改善措施将超出的事故降低到平均值；将样本数据代入计算公式，计算所有地点总体、物损、伤亡事故的安全可提高空间，从而判别事故多发道路。