[0029] 结合附图1‑6,对本发明技术方案作进一步描述:
[0030] 一、系统总体功能的设计
[0031] 本发明一种油气浓度大数据智能检测系统由由测量参数采集平台和多点油气浓度预测子系统两部分组成,测量参数采集平台实现对油气浓度和影响油气浓度的参数进行检测,多点油气浓度预测子系统实现对测量参数进行处理与预测;测量参数采集平台包括测量参数的检测节点、网关节点、现场监控端、云平台和手机App组成,检测节点构建成CAN总线网络通信来实现检测节点和网关节点之间的CAN现场总线网络通信;检测节点将检测的测量参数通过网关节点发送给现场监控端和云平台,网关节点、云平台、现场监控端和手机App之间实现测量参数和相关控制信息的双向传输。测量参数采集平台见图1所示。
[0032] 二、检测节点的设计
[0033] 采用大量基于CAN总线通信网络的检测节点作为测量参数感知终端,检测节点通过CAN总线通信网络实现与网关节点之间的信息相互交互。检测节点包括采集油气浓度、温度和风速的传感器和对应的信号调理电路、MSP430单片机和CAN总线网络通信的CAN总线接口;检测节点的软件主要实现CAN总线网络通信和测量参数的采集与预处理。软件采用C语言程序设计,兼容程度高,大大提高了软件设计开发的工作效率,增强了程序代码的可靠性、可读性和可移植性。检测节点结构见图4。
[0034] 三、网关节点设计
[0035] 网关节点包括CAN接口、NB‑IoT模块、MSP430单片机和RS232接口,网关节点包括CAN接口实现与检测节点之间通信的CAN现场总线通信网络,NB‑IoT模块实现网关节点与云平台之间的数据双向交互,RS232接口连接现场监控端,实现网关与现场监控端之间的信息交互。网关节点结构见图5。
[0036] 四、现场监控端的软件设计
[0037] 现场监控端是一台工业控制计算机,现场监控端主要实现对测量参数进行采集与处理和预测,实现与检测节点与网关节点的信息交互,现场监控端主要功能为通信参数设置、数据分析与数据管理和通过多点油气浓度预测子系统对测量参数进行智能处理与预测,该管理软件选择了Microsoft Visual++6.0作为开发工具,调用系统的Mscomm通信控件来设计通讯程序,现场监控端软件功能见图5;多点油气浓度预测子系统结构如图2所示,多点油气浓度预测子系统由参数检测模块、检测参数融合模型和二元联系数的DRNN神经网络预测模型组成,多点油气浓度预测子系统的设计过程如下:
[0038] 参数检测模块包括多个参数检测模型,多个油气浓度传感器、多个温度传感器、多个风速传感器分别感知被检测环境的油气浓度值、温度值和风速值分别作为对应的参数检测模块的多个参数检测模型的输入,3个参数检测模块的多个参数检测模型输出分别作为对应的检测参数融合模型的输入;参数检测模型的结构功能见图3所示;
[0039] 1、参数检测模型设计
[0040] 参数检测模型由带时滞单元的Adaline神经网络模型、EMD经验模态分解模型、GM(1,1)灰色预测模型、多个NARX神经网络预测模型、2个按拍延迟线TDL、2个ARIMA预测模型和二元联系数的小波神经网络模型组成;
[0041] (1)、带时滞单元的Adaline神经网络模型设计
[0042] 参数检测传感器输出作为带时滞单元的Adaline神经网络模型的输入,带时滞单元的Adaline神经网络模型输出作为EMD经验模态分解模型的输入;带时滞单元的Adaline神经网络模型由2个按拍延迟线TDL和Adaline神经网络组成,参数传感器输出作为对应的按拍延迟线TDL的输入,该按拍延迟线TDL的输出作为Adaline神经网络的输入,Adaline神经网络的输出作为对应的按拍延迟线TDL的输入,该按拍延迟线TDL的输出为带时滞单元的Adaline神经网络模型的输出;Adaline神经网络模型的自适应线性单元(Adaptive Linear Element)是早期的神经网络模型之一,该模型的输入信号可写成向量的形式:X(K)=[x0T(K),x1(K),…xn(K)] ,每一组输入信号对应有一组权值向量相对应表示为:W(K)=[k0(K),k1(K),…k(K)],x0(K)等于负1时是Adaline神经网络模型的偏置值决定神经元的兴奋或抑制状态,可根据Adaline神经网络模型的输入向量和权值向量定义网络输出为:
[0043]
[0044] 在Adaline神经网络模型中,有一特殊输入即理想响应输出d(K),把它送入Adaline神经网络模型中,然后通过网络的输出y(K)进行比较,将差值送到学习算法机制中,以调整权向量直到获得最佳权向量,y(K)与d(K)趋向一致,权向量的调整过程即为网络的学习过程,学习算法是学习过程的核心部分,Adaline神经网络模型的权值优化搜索算法采用LMS算法最小二乘法。
[0045] (2)、EMD经验模态分解模型设计
[0046] 带时滞单元的Adaline神经网络模型输出作为EMD经验模态分解模型的输入,EMD经验模态分解模型输出的测量参数低频趋势值作为GM(1,1)灰色预测模型的输入;EMD经验模态分解模型输出的多个测量参数高频趋势值分别作为对应的多个NARX神经网络预测模型的输入;EMD经验模态分解是一种自适应信号筛选方法,具有计算简单、直观、基于经验和自适应的特点。它能将存在于测量参数信号中不同特征的趋势逐级筛选出来,得到多个高频波动部分(IMF)和低频趋势部分。EMD经验模态分解出来的IMF分量包含了测量参数信号从高到低不同频率段的成分,每个频率段包含的频率分辨率都随信号本身变化,具有自适应多分辨分析特性。使用EMD经验模态分解的目的就是为了更准确地提取故障信息。IMF分量必须同时满足两个条件:①在待分解参数测量信号中,极值点的数目与过零点的数目相等,或最多相差一个;②在任一时间上,由局部极大值和局部极小值定义的包络均值为零。EMD经验模态分解方法针对带时滞单元的Adaline神经网络模型输出值信号的“筛分”过程步骤如下:
[0047] (a)带时滞单元的Adaline神经网络模型输出值信号所有的局部极值点,然后用三次样条线将左右的局部极大值点连接起来形成上包络线。
[0048] (b)在用三次样条线将带时滞单元的Adaline神经网络模型输出值的局部极小值点连接起来形成下包络线,上、下包络线应该包络所有的数据点。
[0049] (c)上、下包络线的平均值记为m1(t),求出:
[0050] x(t)‑m1(t)=h1(t) (2)
[0051] x(t)为带时滞单元的Adaline神经网络模型输出值原始信号,如果h1(t)是一个IMF,那么h1(t)就是x(t)的第一个IMF分量。记c1(t)=h1k(t),则c1(t)为信号x(t)的第一个满足IMF条件的分量。
[0052] (d)将c1(t)从x(t)中分离出来,得到:
[0053] r1(t)=x(t)‑c1(t) (3)
[0054] 将r1(t)作为原始数据重复步骤(a)‑步骤(c),得到x(t)的第2个满足IMF条件的分量c2。重复循环n次,得到信号x(t)的n个满足IMF条件的分量。这样通过EMD经验模态分解模型把带时滞单元的Adaline神经网络模型输出分解为低频趋势部分和多个高频波动部分,EMD经验分解模型如图3所示。
[0055] (3)、GM(1,1)灰色预测模型设计
[0056] EMD经验模态分解模型输出的测量参数低频趋势值作为GM(1,1)灰色预测模型的输入,GM(1,1)灰色预测模型输出和多个NARX神经网络预测模型输出分别作为二元联系数的小波神经网络模型的对应输入;;GM(1,1)灰色预测方法较传统的统计预测方法有着较多的优点,它不需要确定预测变量是否服从正态分布,不需要大的样本统计量,不需要根据测量参数低频趋势值输入变量的变化而随时改变预测模型,通过累加生成技术,建立统一的微分方程模型,累加测量参数低频趋势原始值还原后得出预测结果,微分方程模型具有更高的预测精度。建立GM(1,1)灰色预测模型的实质是对低频趋势值原始数据作一次累加生成,使生成数列呈现一定规律,通过建立微分方程模型,求得拟合曲线,用以对参数测量低频趋势值进行预测。
[0057] (4)、多个NARX神经网络预测模型设计
[0058] EMD经验模态分解模型输出的多个测量参数高频趋势值分别作为对应的多个NARX神经网络预测模型的输入,GM(1,1)灰色预测模型输出和多个NARX神经网络预测模型输出分别作为二元联系数的小波神经网络模型的对应输入;NARX神经网络预测模型是一种带输出反馈连接的动态递归神经网络,在拓扑连接关系上可等效为有输入时延的BP神经网络加上输出到输入的时延反馈连接,其结构由输入层、时延层、隐层和输出层构成,其中输入层节点用于信号输入,时延层节点用于输入信号和输出反馈信号的时间延迟,隐层节点利用激活函数对时延后的信号做非线性运算,输出层节点则用于将隐层输出做线性加权获得最终网络输出。NARX神经网络第i个隐层节点的输出hi为:
[0059]
[0060] NARX神经网络第j个输出层节点输出oj为:
[0061]
[0062] 本发明专利的NARX神经网络包括输入层、时延层、隐层和输出层。
[0063] (5)、ARIMA预测模型设计
[0064] 2个按拍延迟线TDL输出分别作为对应的ARIMA预测模型的输入,2个ARIMA预测模型输出作为二元联系数的小波神经网络模型的对应输入;ARIMA(Auto regressive Integrated Moving Average)预测模型是自回归积分滑动平均模型,它将自回归模型(Autoregressive,AR)和滑动平均模型(Moving Average,MA)有机地组合起来,使之成为一种综合的预测方法。作为有效的现代数据处理方法之一,它被誉为时间序列预测方法中最复杂最高级的模型,在实际应用中,由于输入原始数据序列往往表现出一定的趋势或循环特征,不满足ARMA模型对时间序列的平稳性要求,而取差分是消除数据趋势性的一种方便和有效的方法。基于差分后的数据序列建立的模型称为ARIMA预测模型,记为{Xt}‑ARIMA(p,d,q),其中p、q称为模型的阶,d表示差分的次数。显然,当d为0时,ARIMA预测模型为ARMA模型,其定义为:
[0065] xt=b1xt‑1+…+bpxt‑p+εt+a1εt‑1+…+aqεt‑q (6)
[0066] {xt}为要预测的二元联系数的小波神经网络输出的压力确定值a和波动值b的数2
据序列,{εt}~WN(0,σ)。ARIMA预测模型建立主要包括模型的识别、参数估计和模型诊断。
模型识别主要包括时间序列的预处理和模型参数的初步定阶;模型定阶完成之后需要通过时间序列观察值并结合p,d,q值来对模型中的未知参数进行估计;模型的诊断主要是针对整个模型的显著性检验和模型中参数的显著性检验。通常模型的建立是个不断优化的过程,模型优化常用的为AIC和BIC准则,即最小信息量准则其值越小,模型越合适,BIC准则是针对AIC准则对大样本序列的不足所做的改进。
[0067] (6)、二元联系数的小波神经网络模型设计
[0068] GM(1,1)灰色预测模型输出和多个NARX神经网络预测模型输出分别作为二元联系数的小波神经网络模型的对应输入,二元联系数的小波神经网络模型输出测量参数的确定值c和波动值d构成测量参数的二元联系数为c+di,测量参数的确定值c和波动值d分别作为对应的2个按拍延迟线TDL的对应输入和二元联系数的小波神经网络模型的对应输入,2个按拍延迟线TDL输出分别作为对应的ARIMA预测模型的输入,2个ARIMA预测模型输出作为二元联系数的小波神经网络模型的对应输入,二元联系数的小波神经网络模型是输出为二元联系数的小波神经网络模型,小波神经网络模型WNN(Wavelet Neural Networks)是在小波理论基础上,结合人工神经网络而提出的一种前馈型网络。它是以小波函数为神经元的激励函数,小波的伸缩、平移因子以及连接权重,在对误差能量函数的优化过程中被自适应调整。设小波神经网络模型的输入信号可以表示为输入的一维向量xi(i=1,2,…,n),输出信号表示为yk(k=1,2,…,m),小波神经网络模型输出层输出值的计算公式为:
[0069]
[0070] 公式中ωij输入层i节点和隐含层j节点间的连接权值, 为小波基函数,bj为小波基函数的平移因子,aj小波基函数的伸缩因子,ωjk为隐含层j节点和输出层k节点间的连接权值。本专利中的小波神经网络模型的权值和阈值的修正算法采用梯度修正法来更新网络权值和小波基函数参数,从而使小波神经网络输出不断逼近期望输出。二元联系数的小波神经网络模型输出作为被测量参数的二元联系数值。
[0071] 2、检测参数融合模型设计
[0072] ①、构建参数测量传感器的时间序列二元联系数阵列
[0073] 一段时间多个参数检测模型输出的二元联系数构成参数测量传感器的时间序列二元联系数阵列,设有n个参数测量传感器和m个时刻的nm个参数测量传感器的二元联系数值构成n行和m列的参数测量传感器的时间序列二元联系数阵列,设同一参数测量传感器的不同时刻二元联系数值为Aij(t),Aij(t+1),…,Aij(m),则所有参数测量传感器的时间序列二元联系数阵列为:
[0074]
[0075] ②、计算参数测量传感器的时间序列二元联系数值的距离融合权重
[0076] 同一时刻所有参数测量传感器的二元联系数的平均值构成二元联系数阵列的正理想值,二元联系数阵列的正理想值为:
[0077] 同一时刻所有参数测量传感器的二元联系数值与对应的二元联系数阵列的正理想值贴近度最小的二元联系数值构成二元联系数阵列的负理想值,二元联系数阵列的负理想值为:
[0078]
[0079] 每个参数测量传感器的时间序列二元联系数值的正理想值贴近度是每个参数测量传感器的时间序列二元联系数与二元联系数阵列的正理想值的贴近度为:
[0080]
[0081] 每个参数测量传感器的时间序列二元联系数值的负理想值贴近度是每个参数测量传感器的时间序列二元联系数值与二元联系数阵列的负理想值的贴近度为:
[0082]
[0083] 每个参数测量传感器的时间序列二元联系数值的正理想值贴近度除以该参数测量传感器的时间序列二元联系数值的负理想值贴近度与该参数测量传感器的时间序列二元联系数值的正理想值贴近度的和得到的商为每个参数测量传感器的时间序列二元联系数值的距离相对贴近度,公式为:
[0084]
[0085] 通过(16)公式计算可以知道,每个参数测量传感器的时间序列二元联系数值的距离相对贴近度越大,则该参数测量传感器的时间序列二元联系数值与对应的正理想值就越接近,否则该参数测量传感器的时间序列二元联系数值与对应的正理想值就越远,根据这个原理确定每个参数测量传感器的时间序列二元联系数值的距离相对贴近度除以所有参数测量传感器的时间序列二元联系数值的距离相对贴近度的和得到的商为每个参数测量传感器的时间序列二元联系数值的距离贴近度融合权重为:
[0086]
[0087] ③、计算参数测量传感器的时间序列二元联系数值的灰色关联度融合权重[0088] 每个参数测量传感器的时间序列二元联系数值与二元联系数阵列正理想值的灰色关联度为:
[0089]
[0090] 通过计算每个参数测量传感器的时间序列二元联系数值与二元联系数阵列正理想值的灰色关联度,可以构建每个参数测量传感器的时间序列二元联系数值的灰色关联度矩阵:
[0091]
[0092] 根据公式(16)可以得到每个参数测量传感器的时间序列二元联系数值与时间序列参数测量传感器的二元联系数阵列正理想值之间的灰色关联度,如下公式所示:
[0093]
[0094] 同理,每个参数测量传感器的时间序列二元联系数值与时间序列参数测量传感器的二元联系数阵列负理想值的灰色关联度为,定义如下公式:
[0095]
[0096] 同理,通过计算每个参数测量传感器的时间序列二元联系数值与二元联系数阵列负理想值的灰色关联度,可以构建每个参数测量传感器的时间序列二元联系数值的灰色关联度矩阵:
[0097]
[0098] 根据公式(19)可以得到每个参数测量传感器的时间序列二元联系数值与二元联系数阵列的负理想值之间的灰色关联度,如下公式所示:
[0099]
[0100] 每个参数测量传感器的时间序列二元联系数值与二元联系数阵列的正理想值之间的灰色关联度除以该参数测量传感器的时间序列二元联系数值与二元联系数阵列的正理想值之间的灰色关联度相加该参数测量传感器的时间序列二元联系数值与二元联系数阵列的负理想值之间的灰色关联度的和得到的商为该参数测量传感器的时间序列二元联系数值的灰色关联相对贴近度为:
[0101]
[0102] 通过(21)公式计算可以知道,每个参数测量传感器的时间序列二元联系数值的灰色关联相对贴近度越大,则该参数测量传感器的时间序列二元联系数值与对应的二元联系数阵列的正理想值的形状相似度越大,否则该参数测量传感器的时间序列二元联系数值与对应的二元联系数阵列的正理想值的形状相似度越小,根据这个原理确定每个参数测量传感器的时间序列二元联系数值的灰色关联相对贴近度除以所有参数测量传感器的时间序列二元联系数值的灰色关联相对贴近度的和得到的商为该参数测量传感器的时间序列二元联系数值的灰色关联度融合权重为:
[0103]
[0104] ④、计算多个参数测量传感器的时间序列二元联系数值的融合值
[0105] 根据每个参数测量传感器的时间序列二元联系数值的距离融合权重αi和灰色关联度融合权重βi,求取均方根组合权重γi,显然γi与αi、βi和都应尽可能接近,根据最小相对信息熵原理有:
[0106]
[0107] 用拉格朗日乘子法解上述优化问题得:
[0108]
[0109] 根据公式(24)可以知每个参数测量传感器的时间序列二元联系数值的距离融合权重和灰色关联度融合权重乘积的均方根占所有参数测量传感器的时间序列二元联系数值的距离融合权重和灰色关联度融合权重乘积的均方根和的比为该参数测量传感器的时间序列二元联系数值融合的均方根组合权重。
[0110] 根据参数测量传感器的时间序列二元联系数值的距离融合权重αi和灰色关联度融合权重βi进行线性组合得到该参数测量传感器的时间序列二元联系数值融合的线性组合权重θi,公式为:
[0111] θi=ααi+ββi (25)
[0112] 根据公式(27)和公式(28)得到该参数测量传感器的时间序列二元联系数值融合的二元联系数融合权重为wi:
[0113] wi=[min(θi,γi),max(θi,γi)] (26)
[0114] 从公式(26)可以知每个参数测量传感器的时间序列二元联系数值的均方根组合权重、线性组合权重重按照从小到大排序构成该参数测量传感器的时间序列二元联系数值的组合融合权重。根据同一时刻每个参数测量传感器的时间序列二元联系数值与该参数测量传感器的时间序列二元联系数值的组合融合权重的积相加得到的和为所有参数测量传感器的时间序列二元联系数融合值为:
[0115]
[0116] 3、二元联系数的DRNN神经网络预测模型设计
[0117] 3个检测参数融合模型输出作为二元联系数的DRNN神经网络预测模型的对应输入,二元联系数的DRNN神经网络预测模型输出油气浓度的确定值a和波动值b构成油气浓度的二元联系数为a+bi,油气浓度的确定值a和波动值b分别作为二元联系数的DRNN神经网络预测模型的对应输入;二元联系数的DRNN神经网络预测模型为输出和输入都是二元联系数的DRNN神经网络模型,二元联系数的DRNN神经网络预测模型为3层网络结构,其隐层为回归层,在本专利二元联系数的DRNN神经网络预测模型中,设I=[I1(t),I2(t),…,In(t)]为网络输入向量,其中Ii(t)为二元联系数的DRNN神经网络预测模型第i个神经元t时刻的输入,回归层第j个神经元的输出为Xj(t),Sj(t)为第j个回归神经元输入总和,f(·)为S的函数,则O(t)为二元联系数的DRNN神经网络预测模型输出。则二元联系数的DRNN神经网络预测模型输出为:
[0118]
[0119] 二元联系数的DRNN神经网络预测模型输出作为油气浓度的二元联系数值。
[0120] 五、一种油气浓度大数据检测系统的设计举例
[0121] 根据油气浓度检测系统的实际状况,系统布置测量参数采集平台的检测节点、网关节点和现场监控端的平面布置安装图,其中检测节点的传感器根据检测的需要均衡布置在油气浓度挥发的各个方位,通过该系统实现对油气浓度参数进行采集与预测。
[0122] 本发明方案所公开的技术手段不仅限于上述实施方式所公开的技术手段,还包括由以上技术特征任意组合所组成的技术方案。应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰这些改进和润饰也视为本发明的保护范围。