一种基于生产台账的涉嫌非法运输危化品企业的辨识方法与流程

本发明涉及货物运输领域，具体涉及一种基于生产台账的涉嫌非法运输危化品企业的辨识方法。

背景技术：

危化品是指具有爆炸、易燃、毒害、感染、腐蚀、放射性等特性，在运输、储存、生产、经营、使用和处置中，易造成人身伤亡、财产损毁和环境污染而需要特别防护的物质及物品。随着国民经济的持续快速发展，危化品运输量逐年上升，在这其中，非法运输量占比不轻。

目前有不少学者对危化品的运输进行了研究，但大多数都是研究为危化品的合法运输，例如对合法运输的安全风险进行评估评价，而辨识危化品企业是否存在非法营运这一块的研究少之又少，危化品企业非法营运危化品对公共安全存在非常大的隐患。

技术实现要素：

本发明要解决的技术问题是提供一种基于企业的台账数据，建立模型，预测某时期企业的生产量、运输量，判断出某危化品企业是否涉嫌非法营运原材料、成品的辨识方法。

为了解决上述技术问题，本发明采用以下技术方案：

本发明提供一种基于生产台账的涉嫌非法运输危化品企业的辨识方法，包括如下步骤：

步骤一、基于危化品企业台账记录的原材料入库量、出库量、库存量，危化品成品生产量、入库量、出库量、库存量、运输量，建立时间序列模型进行多时间段预测，从而判断企业是否涉嫌非法运输危化品原材料；

步骤二、基于危化品企业台账数据以各主要原材料、燃料、水等为自变量，以危化成品生产量为因变量，建立多元回归模型进行预测，从而判断企业是否涉嫌非法运输危化成品；

步骤三、通过视频监控车辆，比对台账记录的合法车辆的运输情况及非企业车辆运输情况，判断是否涉嫌非法运输行为；

所述步骤一中的时间序列模型如下：

即：t期估计值＝α*(t期实际值)+(1-α)*(t-1期估计值)；

其中，s为危化成品生产量、y为运输量，α取值一般在0.3-0.5之间，将t时期根据数据变化节点分为多时间段，聚类后再预测能使预测结果；

所述步骤二中多元回归模型如下：

ht＝b0+b1x1+b2x2+…+bkxk+∈

其中，ht为t时期预测的威化成品生产量，b0，b1，b2,...,bk是模型的参数，∈为误差项；

所述辨识方法还包括应用加权平均法融合步骤一中时间序列模型和步骤二中多元回归模型预测的危化品生产量，得到目标时期成品生产量最终预测值；

所述加权平均融合法的计算公式如下：

其中，lt为t时期危化成品生产量的最终预测值。

进一步的，所述步骤一还包括通过将t时期运用时间序列模型预测的值比对企业台账记录的原材料及成品的各数据，计算累计差，判断两者比对的累计差是否在容忍区间里，从而判断该企业是否涉嫌非法运输原材料及成品，比对的累计差公式如下:

其中，∑δ为预测值与台账数据s的累计差。

进一步的，所述步骤二还包括通过t时期运用多元回归模型预测的生产值比对企业台账记录的危化品成品生产量，计算累计差，判断两者比对的累计差是否在某容忍区间值，从而判断企业是否涉嫌非法运营危化品，比对的累计差公式如下：

其中，∑φ为预测值ht与台账数据h的累计差。

进一步的，所述加权平均融合法通过t时期运用最终预测值比对台长数据，计算累计差，比对的累计差公式如下：

∑γ＝lt-l

其中，∑γ为最终预测值lt与台账数据l的累计差。

进一步的，所述步骤三还包括观察台账记录数据与预测数据在t时期的库存量与运输量是否符合总库存量计算等式，总库存计算等式如下：

kt1＝kt2+kt3+kt4

其中，kt1为t时期总库存量，kt2为t时期剩余库存量，kt3为t时期合法运输量，kt4为t时期非法运输量。

进一步的，所述步骤一的时间序列模型使用python进行多时间段预测。

进一步的，所述步骤一中数据聚类后：

a、当每组数据量都小于20时，第一个数据取前三项的均值；

b、当每组数据量大于20时，所有数据正常计算。

进一步的，所述辨识方法还包括基于车辆运输记录长期累计运输量明显差别，出入时间不对，从而判断企业涉嫌非法运输危化品。

进一步的，所述辨识方法还包括基于视频监控的车辆数据，存在企业未登记的车辆频繁出入，从而判断企业涉嫌非法运输危化品。

本发明的有益效果在于：

1、本发明方法从企业源头出发，基于危化品企业台账，危化品在时间上有着明显的供应需求规律，可通过建立时间序列模型预测危化品原材料入库量、出库量、消耗量、库存量，危化成品生产量、运输量、入库量、出库量、库存量；同时建立多元回归模型算法，选取相关关键性变量为自变量，以危化品的成品生产量为因变量进行预测某时期危化成品生产量；比对某时间段内危化成品的入库量、出库量、库存量、运输量，运用加权平均方法融合两种预测方法预测的危化成品生产量，将台账的记录数据与之对比分析，计算某时期预测值与记录数据的累计差大小，从而辨识危化品企业是否涉嫌非法营运危化成品。进一步分析原材料，同理将时间序列模型预测的各主要原材料入库量、出库量、库存量与台账记录数据比对分析，辨识危化品企业是否涉嫌非法运输危化品原材料。判断企业是否涉嫌危化品的非法运营。

2、本发明方法基于视频监控的车辆数据，比对台账记录的企业登记的合法车辆运输情况和视频监控的车辆运输情况，若长期与车辆运输记录在运输量、出入时间有明显差别，可判定企业涉嫌非法运输危化品；若长期有企业未登记的车辆频繁出入，可判断企业涉嫌非法运输危化品。根据运输量，基于剩余库存思想，判断企业是否涉嫌危化品的非法运营。

附图说明

为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。

图1是本发明方法基于生产台账的涉嫌非法运输危化品企业的辨识方法的流程图。

具体实施方式

为了更好地阐述该发明的内容，下面通过具体实施例对本发明进一步的验证。特在此说明，实施例只是为更直接地描述本发明，它们只是本发明的一部分，不能对本发明构成任何限制。

本发明实施例提供一种基于生产台账的涉嫌非法运输危化品企业的辨识方法，具体包括如下步骤：

步骤一：第一层面，根据得到的企业台账数据，包含某时间段内各原材料入库量(r1)、出库量(c1)、库存量(k1)、消耗量(x)，危化成品生产量(s)、入库量(r2)、出库量(c2)、运输量(y)、库存量(k2)，生产成品时所需燃料用量、电量、水量，车辆出入运输情况等。选取的危化品原材料关键变量为a、b、c，时间段长为58天，其入库量、出库量、库存量整理如下表1，基于库存量、出库量和入库量进行初步分析判断是否有涉嫌非法营运的可能；

表1企业台账各原材料数据

步骤二：运用时间序列模型以时间为轴训练数据，预测目标时期的危化成品生产量。手动将数据聚类，聚类后每组数据量都小于20，因此第一个数值取前三项的均值：

s1＝y1+y2+y3

由：

借助python进行预测，经调试，α取0.4效果较好，得到预测结果如表3所示，python代码如下所示：

建立多元回归模型，以各主要原材料(a、b、c)消耗量、水(d)消耗量等为自变量，消耗量情况见下表2，以危化成品生产量为因变量，预测目标时期危化成品生产量。

由：

ht＝b0+b1x1+b2x′2+…+bkxk+∈

借助excel回归模型，得到回归表达式，预测结果如表3：

ht＝1.822273a+0.325682d–27.1818

表2各主要原材料及水的消耗量

步骤三：将时间序列模型预测的危化成品生产量与多元回归模型预测的危化成品生产量进行融合得目标时期危化成品生产量最终预测值。

由：

得到目标时期危化成品生产量最终预测值，结果如下表3：

表3各方法预测危化成品量比较

步骤四：将各方法预测的值与台账数据比对，根据两者的累计差是否在容忍区间内，判断危化企业是否涉嫌危化品的运输。

由：

累计差＝预测值–台账值

得每个预测值与台账数据的差值，进行累加，可得：

一次指数平滑模型累计差为

多元回归模型累计差为∑δ＝ht–h＝103.76

加权平均融合法累计差为∑δ＝lt–l＝69.08

假定容忍值为30，即∑δ不在容忍区间内，则可判定该危化企业涉嫌危化成品非法运输，可采取侦查手段进行查证；

步骤五：同理比对用时间序列(一次指数平滑法)预测的关键原材料与台账数据，进一步具体分析，若是原材料累计差不在容忍区间，则判定该企业涉嫌运输危化品原材料，可采取侦查手段进行查证；

步骤六：第二层面，利用视频监控该企业合法车辆运输情况，长期比对企业台账记录的企业车辆运输情况，例如根据车牌识别出某车辆在某时间内的出入情况、运输量情况是否和台账记录一致，判断企业合法运输车辆是否存在非法运输；根据车牌辨别某经常性出入车辆是否为该企业的运输车辆，是否存在外来车辆非法运输；

表4车辆运输情况比较

从表4可以看出比对台账记录的运输情况和视频监控的车辆运输情况中，车牌为651和438的企业登记合法车辆出入时间不一致，涉嫌非法运输；车牌u569的车辆频繁出入企业周边，且在台账记录的合法车辆出入时间内出现，涉嫌非法运输危化品。

根据运输量，基于剩余库存思想，即总库存量＝剩余库存量+合法运输量+非法运输量，计算企业的危化品原材料及成品的运输情况，以危化成品为例，如下表5：

表5危化成品出入库量与运输量

即：

非法运输量＝库存总量–剩余库存量–合法运输量

＝31633.09–31035＝598.08

步骤七：基于企业台账，若企业源头生产与车辆两方面都有或有其一某种情况不符合，则可判定该危化品企业涉嫌危化品的运营，由以上操作计算，建议采取侦查手段进行查证，辨识过程结束。

本文中所描述的具体实施例仅仅是对本发明精神作举例说明。本发明所述技术领域的技术人员可以所描述的具体实施例做各种各样的修改或补充或采用类似的方式替代，但并不会偏离本发明的精神或者超越所附权利要求书所定义的范围。