发明名称 基于三点模型的柴油机故障预测方法
摘要 本发明的目的在于提供基于三点模型的柴油机故障预测方法,首先采集柴油机运行状态数据并进行数值滤波,获得无干扰时序数列。继而拟合时序状态参数变化趋势,建立柴油机三点预测模型,据此前向外推,获得预测数据,从而实现对柴油机未来工作状态的趋势预测。根据本发明的故障诊断结果,可以在设备劣化初期准确预测出柴油机未来劣化趋势,指导工作人员对柴油机视情维修;本发明基于三点模型预测柴油机故障,具有不需要样本数据全部统计特征、运算简捷、使用方便、结果准确以及尤其适用于近期预测等优点。
申请公布号 CN104295373A 申请公布日期 2015.01.21
申请号 CN201410525303.5 申请日期 2014.10.08
申请人 哈尔滨工程大学 发明人 李文辉;董佳莹;费景洲;王金鑫;王忠巍
分类号 F02B77/08(2006.01)I 主分类号 F02B77/08(2006.01)I
代理机构 代理人
主权项 基于三点模型的柴油机故障预测方法,其特征是:(1)等时采集检测柴油机运行参数,所采集的运行参数包括:功率、油耗率、转速、烟度、噪声、滑油压力、排气温度、增压后进气压力;(2)将获得的柴油机运行参数进行中位值滤波,将含有噪声的运行参数转换为无干扰信号,对于每一项参数,按照时间顺序,由远及近分别构成建模原序列;(3)将所获取的无干扰信号作为建模数据,获取三点坐标:三点坐标的获取步骤如下:a、基于三点法的预测模型包括三种:二次抛物线型、直线型与指数曲线型、三次抛物线型,从中选择一种预测模型,建立柴油机运行参数三点模型;b、确定三点坐标:对于柴油机的一项运行参数,三点纵坐标的确定应满足如下原则:1)原序列总项数N应当是奇数,此时N=n,当总项数为偶数时,则去除第一项检测数据,使N奇数化,此时N‑1=n,n为实际序列总项数,同时,n不小于9;2)若实际序列总项数9≤n&lt;15,则在保持检测数据排列顺序不变的前提下,抽取实际序列的前3项、正中3项、后3项检测数据,分别构成初期、中期、近期子序列;每个子序列算出一个加权平均数,得到的三个加权平均数就是需要确定的三点的纵坐标;对于每一子序列,检测数据对应权数的选择由远及近分别取1、2、3;3)若实际序列总项数n≥15,则在保持检测数据排列顺序不变的前提下,抽取原序列的前5项、正中5项、后5项检测数据,分别构成初期、中期、近期子序列。每个子序列算出一个加权平均数,得到的三个加权平均数就是需要确定的三点的纵坐标;对于每一子序列,检测数据对应权数由远及近依次选取1、2、3、4、5;在满足以上三个原则的前提下,将初期、中期、近期三点的坐标分别设为M<sub>1</sub>(t<sub>1</sub>,R),M<sub>2</sub>(t<sub>2</sub>,S),M<sub>3</sub>(t<sub>3</sub>,T)其中,M<sub>1</sub>,M<sub>2</sub>,M<sub>3</sub>分别代表初期、中期、近期三点;t<sub>1</sub>,t<sub>2</sub>,t<sub>3</sub>分别为初期、中期、近期三点横坐标;R,S,T分别为初期、中期、近期三点纵坐标;设正中项数<img file="FDA0000582702850000021.GIF" wi="260" he="141" />实际序列各项检测值为y<sub>1</sub>、y<sub>2</sub>、…y<sub>d</sub>…y<sub>n</sub>;因此对于五项加权平均,三点的纵坐标为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>R</mi><mo>=</mo><mfrac><mn>1</mn><mn>15</mn></mfrac><mrow><mo>(</mo><msub><mi>y</mi><mn>1</mn></msub><mo>+</mo><msub><mrow><mn>2</mn><mi>y</mi></mrow><mn>2</mn></msub><mo>+</mo><msub><mrow><mn>3</mn><mi>y</mi></mrow><mn>3</mn></msub><mo>+</mo><msub><mrow><mn>4</mn><mi>y</mi></mrow><mn>4</mn></msub><mo>+</mo><msub><mrow><mn>5</mn><mi>y</mi></mrow><mn>5</mn></msub><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000582702850000022.GIF" wi="696" he="121" /></maths><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>S</mi><mo>=</mo><mfrac><mn>1</mn><mn>15</mn></mfrac><mrow><mo>(</mo><msub><mi>y</mi><mrow><mi>d</mi><mo>-</mo><mn>2</mn></mrow></msub><mo>+</mo><mn>2</mn><msub><mi>y</mi><mrow><mi>d</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>+</mo><mn>3</mn><msub><mi>y</mi><mi>d</mi></msub><mo>+</mo><msub><mrow><mn>4</mn><mi>y</mi></mrow><mrow><mi>d</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>+</mo><msub><mrow><mn>5</mn><mi>y</mi></mrow><mrow><mi>d</mi><mo>+</mo><mn>2</mn></mrow></msub><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000582702850000023.GIF" wi="873" he="127" /></maths><maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>T</mi><mo>=</mo><mfrac><mn>1</mn><mn>15</mn></mfrac><mrow><mo>(</mo><msub><mi>y</mi><mrow><mi>n</mi><mo>-</mo><mn>4</mn></mrow></msub><mo>+</mo><msub><mrow><mn>2</mn><mi>y</mi></mrow><mrow><mi>n</mi><mo>-</mo><mn>3</mn></mrow></msub><mo>+</mo><msub><mrow><mn>3</mn><mi>y</mi></mrow><mrow><mi>n</mi><mo>-</mo><mn>2</mn></mrow></msub><mo>+</mo><msub><mrow><mn>4</mn><mi>y</mi></mrow><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>+</mo><msub><mrow><mn>5</mn><mi>y</mi></mrow><mi>n</mi></msub><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000582702850000024.GIF" wi="867" he="128" /></maths>对于三项加权平均,三点的纵坐标为:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mi>R</mi><mo>=</mo><mfrac><mn>1</mn><mn>6</mn></mfrac><mrow><mo>(</mo><msub><mi>y</mi><mn>1</mn></msub><mo>+</mo><mn>2</mn><msub><mi>y</mi><mn>2</mn></msub><mo>+</mo><mn>3</mn><msub><mi>y</mi><mn>3</mn></msub><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000582702850000025.GIF" wi="463" he="128" /></maths><maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mi>S</mi><mo>=</mo><mfrac><mn>1</mn><mn>6</mn></mfrac><mrow><mo>(</mo><msub><mi>y</mi><mrow><mi>d</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>+</mo><mn>2</mn><msub><mi>y</mi><mi>d</mi></msub><mo>+</mo><msub><mrow><mn>3</mn><mi>y</mi></mrow><mrow><mi>d</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000582702850000026.GIF" wi="517" he="125" /></maths><maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><mi>T</mi><mo>=</mo><mfrac><mn>1</mn><mn>6</mn></mfrac><mrow><mo>(</mo><msub><mi>y</mi><mrow><mi>n</mi><mo>-</mo><mn>2</mn></mrow></msub><mo>+</mo><mn>2</mn><msub><mi>y</mi><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>+</mo><mn>3</mn><msub><mi>y</mi><mi>n</mi></msub><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000582702850000027.GIF" wi="527" he="127" /></maths>同样的,三点的横坐标也取各项检测值对应时间点的加权平均数;对于五项加权平均,三点横坐标为:<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><msub><mi>t</mi><mn>1</mn></msub><mo>=</mo><mfrac><mn>1</mn><mn>15</mn></mfrac><mrow><mo>(</mo><mn>1</mn><mo>+</mo><mn>2</mn><mo>&times;</mo><mn>2</mn><mo>+</mo><mn>3</mn><mo>&times;</mo><mn>3</mn><mo>+</mo><mn>4</mn><mo>&times;</mo><mn>4</mn><mo>+</mo><mn>5</mn><mo>&times;</mo><mn>5</mn><mo>)</mo></mrow><mo>=</mo><mfrac><mn>11</mn><mn>3</mn></mfrac><mo>=</mo><mn>3</mn><mo>+</mo><mfrac><mn>2</mn><mn>3</mn></mfrac></mrow>]]></math><img file="FDA0000582702850000028.GIF" wi="1034" he="118" /></maths><maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><msub><mi>t</mi><mn>2</mn></msub><mo>=</mo><mi>d</mi><mo>+</mo><mfrac><mn>2</mn><mn>3</mn></mfrac><mo>=</mo><mfrac><mrow><mi>n</mi><mo>+</mo><mn>1</mn></mrow><mn>2</mn></mfrac><mo>+</mo><mfrac><mn>2</mn><mn>3</mn></mfrac><mo>=</mo><mfrac><mrow><mn>3</mn><mi>n</mi><mo>+</mo><mn>7</mn></mrow><mn>6</mn></mfrac></mrow>]]></math><img file="FDA0000582702850000029.GIF" wi="647" he="133" /></maths><maths num="0009" id="cmaths0009"><math><![CDATA[<mrow><msub><mi>t</mi><mn>3</mn></msub><mo>=</mo><mrow><mo>(</mo><mi>n</mi><mo>-</mo><mn>2</mn><mo>)</mo></mrow><mo>+</mo><mfrac><mn>2</mn><mn>3</mn></mfrac><mo>=</mo><mi>n</mi><mo>-</mo><mfrac><mn>4</mn><mn>3</mn></mfrac></mrow>]]></math><img file="FDA00005827028500000210.GIF" wi="490" he="126" /></maths>对于五项加权平均,三点的坐标为:<maths num="0010" id="cmaths0010"><math><![CDATA[<mrow><msub><mi>M</mi><mn>1</mn></msub><mrow><mo>(</mo><mfrac><mn>11</mn><mn>3</mn></mfrac><mo>,</mo><mi>R</mi><mo>)</mo></mrow><msub><mi>M</mi><mn>2</mn></msub><mrow><mo>(</mo><mfrac><mrow><mn>3</mn><mi>n</mi><mo>+</mo><mn>7</mn></mrow><mn>6</mn></mfrac><mo>,</mo><mi>S</mi><mo>)</mo></mrow><mo>,</mo><msub><mi>M</mi><mn>3</mn></msub><mrow><mo>(</mo><mi>n</mi><mo>-</mo><mfrac><mn>4</mn><mn>3</mn></mfrac><mo>,</mo><mi>T</mi><mo>)</mo></mrow></mrow>]]></math><img file="FDA00005827028500000211.GIF" wi="819" he="133" /></maths>对于三项加权平均,三点的横坐标为:<maths num="0011" id="cmaths0011"><math><![CDATA[<mrow><msub><mi>t</mi><mn>1</mn></msub><mo>=</mo><mfrac><mn>1</mn><mn>6</mn></mfrac><mrow><mo>(</mo><mn>1</mn><mo>+</mo><mn>2</mn><mo>&times;</mo><mn>2</mn><mo>+</mo><mn>3</mn><mo>&times;</mo><mn>3</mn><mo>)</mo></mrow><mo>=</mo><mfrac><mn>7</mn><mn>3</mn></mfrac><mo>=</mo><mn>2</mn><mo>+</mo><mfrac><mn>1</mn><mn>3</mn></mfrac></mrow>]]></math><img file="FDA00005827028500000212.GIF" wi="690" he="126" /></maths><maths num="0012" id="cmaths0012"><math><![CDATA[<mrow><msub><mi>t</mi><mn>2</mn></msub><mo>=</mo><mi>d</mi><mo>+</mo><mfrac><mn>1</mn><mn>3</mn></mfrac><mo>=</mo><mfrac><mrow><mi>n</mi><mo>+</mo><mn>1</mn></mrow><mn>2</mn></mfrac><mo>+</mo><mfrac><mn>1</mn><mn>3</mn></mfrac><mo>=</mo><mfrac><mrow><mn>3</mn><mi>n</mi><mo>+</mo><mn>5</mn></mrow><mn>6</mn></mfrac></mrow>]]></math><img file="FDA00005827028500000213.GIF" wi="612" he="131" /></maths><maths num="0013" id="cmaths0013"><math><![CDATA[<mrow><msub><mi>t</mi><mn>3</mn></msub><mo>=</mo><mrow><mo>(</mo><mi>n</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>+</mo><mfrac><mn>1</mn><mn>3</mn></mfrac><mo>=</mo><mi>n</mi><mo>-</mo><mfrac><mn>2</mn><mn>3</mn></mfrac></mrow>]]></math><img file="FDA00005827028500000214.GIF" wi="455" he="127" /></maths>对于三项加权平均,三点的坐标为:<maths num="0014" id="cmaths0014"><math><![CDATA[<mrow><msub><mi>M</mi><mn>1</mn></msub><mrow><mo>(</mo><mfrac><mn>7</mn><mn>3</mn></mfrac><mo>,</mo><mi>R</mi><mo>)</mo></mrow><mo>,</mo><msub><mi>M</mi><mn>2</mn></msub><mrow><mo>(</mo><mfrac><mrow><mn>3</mn><mi>n</mi><mo>+</mo><mn>5</mn></mrow><mn>6</mn></mfrac><mo>,</mo><mi>S</mi><mo>)</mo></mrow><mo>,</mo><msub><mi>M</mi><mn>3</mn></msub><mrow><mo>(</mo><mi>n</mi><mo>-</mo><mfrac><mn>2</mn><mn>3</mn></mfrac><mo>,</mo><mi>T</mi><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000582702850000031.GIF" wi="800" he="130" /></maths>(4)计算趋势线模型参数估计值,并生成预测数据:对于二次抛物线型,其预测趋势线模型为<img file="FDA0000582702850000032.GIF" wi="335" he="76" />对于直线型和指数曲线型,其预测趋势线模型分别为<img file="FDA0000582702850000033.GIF" wi="210" he="66" />和<img file="FDA0000582702850000034.GIF" wi="183" he="77" />对于三次抛物线型,其预测趋势线模型为<img file="FDA0000582702850000035.GIF" wi="442" he="81" />将步骤(3)中三点坐标带入,求得各预测模型的参数估计值,完成建模;在此基础上,生成预测数据;(5)返回步骤(3)建立模型,直至获得所有所需柴油机运行状态参数预测值。
地址 150001 黑龙江省哈尔滨市南岗区南通大街145号哈尔滨工程大学科技处知识产权办公室