CN106233115B - 估计飞行器发动机的物理参数的测量值是否正常的方法 - Google Patents
估计飞行器发动机的物理参数的测量值是否正常的方法 Download PDFInfo
- Publication number
- CN106233115B CN106233115B CN201580016936.6A CN201580016936A CN106233115B CN 106233115 B CN106233115 B CN 106233115B CN 201580016936 A CN201580016936 A CN 201580016936A CN 106233115 B CN106233115 B CN 106233115B
- Authority
- CN
- China
- Prior art keywords
- value
- exec
- variance
- running parameter
- parameter
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 39
- 230000005856 abnormality Effects 0.000 claims abstract description 17
- 238000005259 measurement Methods 0.000 claims abstract description 17
- 238000012545 processing Methods 0.000 claims description 21
- 230000002159 abnormal effect Effects 0.000 claims description 12
- 238000013500 data storage Methods 0.000 claims description 12
- 238000004590 computer program Methods 0.000 claims description 8
- 239000000446 fuel Substances 0.000 claims description 8
- 238000003860 storage Methods 0.000 claims description 8
- 230000006353 environmental stress Effects 0.000 claims description 3
- 239000012530 fluid Substances 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 claims 1
- 238000002485 combustion reaction Methods 0.000 claims 1
- 239000000463 material Substances 0.000 claims 1
- 238000012360 testing method Methods 0.000 description 29
- 230000006870 function Effects 0.000 description 6
- 238000010586 diagram Methods 0.000 description 4
- 238000009826 distribution Methods 0.000 description 4
- 230000015572 biosynthetic process Effects 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 241000208340 Araliaceae Species 0.000 description 2
- 241000196324 Embryophyta Species 0.000 description 2
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 2
- 235000003140 Panax quinquefolius Nutrition 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 230000002349 favourable effect Effects 0.000 description 2
- 230000014509 gene expression Effects 0.000 description 2
- 235000008434 ginseng Nutrition 0.000 description 2
- 238000002347 injection Methods 0.000 description 2
- 239000007924 injection Substances 0.000 description 2
- 238000005304 joining Methods 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 230000000454 anti-cipatory effect Effects 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000002950 deficient Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000008450 motivation Effects 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 230000035882 stress Effects 0.000 description 1
- 230000001960 triggered effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M15/00—Testing of engines
- G01M15/14—Testing gas-turbine engines or jet-propulsion engines
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B23/00—Testing or monitoring of control systems or parts thereof
- G05B23/02—Electric testing or monitoring
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F11/00—Error detection; Error correction; Monitoring
-
- G—PHYSICS
- G07—CHECKING-DEVICES
- G07C—TIME OR ATTENDANCE REGISTERS; REGISTERING OR INDICATING THE WORKING OF MACHINES; GENERATING RANDOM NUMBERS; VOTING OR LOTTERY APPARATUS; ARRANGEMENTS, SYSTEMS OR APPARATUS FOR CHECKING NOT PROVIDED FOR ELSEWHERE
- G07C11/00—Arrangements, systems or apparatus for checking, e.g. the occurrence of a condition, not provided for elsewhere
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Chemical & Material Sciences (AREA)
- Combustion & Propulsion (AREA)
- Theoretical Computer Science (AREA)
- Quality & Reliability (AREA)
- General Engineering & Computer Science (AREA)
- Automation & Control Theory (AREA)
- Testing And Monitoring For Control Systems (AREA)
- Testing Or Calibration Of Command Recording Devices (AREA)
- Complex Calculations (AREA)
Abstract
本发明涉及一种用于对针对装置(1)的工作参数的值(xexec)的由传感器(20)测量的所述装置(1)的物理参数的值(yexec)的特征是否正常进行估计的方法。所述方法的特征在于包括实施以下步骤:(a)基于与所述多个(xi;yi)i∈[1,n]对有关的回归模型计算物理参数的估计值(yexec);(b)计算相关的残差(resexec);(c)基于与所述多个(xi;yi)i∈[1,n]对有关的方差模型计算物理参数的方差的估计值(varexec);(d)基于所测量的值(yexec)计算异常分数(Zscore),计算估计值(varexec),并计算对于所述多个(xi;yi)i∈[1,n]对的平均残差值(mean);(e)将所测量的值(yexec)的所述异常分数(Zscore)与标准偏差阈值(σ)相比较;并且如果所述异常分数(Zscore)大于所述阈值(σ),(f)在接口装置(33)上将所述测量值识别为异常。
Description
技术领域
本发明涉及航空设备测试领域。
更确切地,本发明涉及一种用于对装置、特别是飞行器发动机的物理参数的测量值是否正常进行估计的方法。
背景技术
测试台是用于在受控条件下的测试中对测试机器(通常为飞行器发动机)的性能进行测量以观察测试机器行为的平台。在这样的测试中,测试台和/或测试机器配备的传感器采集数据,数据被传送到指定的数据库,该数据在下文中也被称为测试数据库,也称为测试库。简单起见,传感器被称为测试台传感器,由此包括在测试范围内的发动机配备的传感器。
测试机器通常是所要测试的开发的原型(所采集的数据因此由研究部门使用以对机器进行改进和开发)或者所要检查规格和可靠性的最终产品(所采集的数据因此由质量团队使用)。替代性地,测试机器可以是整个发动机,或者是发动机的用于部分测试的组件。
然而,由于测试台和/或测试机器的一个或更多个传感器的缺陷经常发生,采集的数据具有异常值或反常值。则这些在“失效”采集过程中获取的数据对测试数据库造成“污染”。
这对库的使用者(特别是研究部门)造成问题,该使用者使用所存储的全部数据来特别进行性能比较。这些数据还能够用于重新调节测试台上的发动机的数值模型。
另外,在测试过程中,快速检测任何有缺陷的传感器是很重要的,以至于能够根据异常的严重性停止测试。事实上,测试是非常昂贵的,并且因此对测试、特别是测试的输出的优化很重要。
法国专利FR2965915描述了一种用于实时监视测试台的传感器的示例性方法,该方法使得能够就传感器故障发出信号,但是这并不能够检测测试中的测量质量的漂移。另外,已知的方法不可避免地与测试台关联且无法在测试库中存储的测量值中后验地识别异常测量值。
因而,期望可靠、高效且可重复地控制与发动机有关的参数的测量,从而简单地识别测量值集中的异常测量值。
发明内容
根据第一方面,本发明提供了一种用于对针对所述装置的工作参数的值的由传感器测量的装置的物理参数的值是否正常进行估计的方法,存储在数据库中的多个对中的每一对限定对于所述工作参数的值的所述物理参数的参考值,所述数据库存储在数据存储装置上,所述方法的特征在于包括由数据处理装置实施以下步骤:
(a)基于与所述多个对有关的回归模型计算对于所述工作参数的值的所述物理参数的估计值;
(b)计算相关的残差;
(c)基于与所述多个对有关的方差模型计算对于所述工作参数的值的所述物理参数的方差的估计值;
(d)基于所述残差、所述方差的估计值和对于所述多个对的平均残差值计算所测量的值的异常分数;
(e)将所测量的值的所述异常分数与标准偏差阈值相比较;
(f)如果所述异常分数高于所述阈值,则在接口装置上告知所述测量值为异常。
基于工作参数对方差进行估计使得(根据使用背景)不再需要测量不确定性中的方差。关于回归模型建立的置信区间(Zscore的公差σ)因此比恒定方差更加真实,并且异常点的检测完全更加精确。告警模型因此更加高效。
根据其他的有利且非限制性特征:
·该方法包括处理所述数据库的所述多个对的在先阶段,该在先阶段包括由数据处理装置实施以下步骤:
(a0)通过基于所述工作参数的值x并根据所述(xi;yi)i∈【1,n】对的集对所述物理参数的值y回归建模来确定与所述多个(xi;yi)i∈【1,n】对有关的所述回归模型,其中,xi;yi表明存储在数据库中的对的值;
(a1)对于所述对中的每一个,计算所述物理参数的估计值和相关的残差;
(a2)计算所述残差的平均值;
(a3)在尺寸为w的滑动窗口上计算残差方差的值的集,每个值与一对的所述工作参数的值有关;
(a4)通过基于所述工作参数的值x并根据所述{xj;varj}j∈【1,n-w+1】的集对所述残差方差的值var回归建模来确定与所述多个对有关的所述方差模型,其中,varj表明所计算的残差方差值,并且xj表明相关的工作参数的值;
·所述在先阶段包括基于所确定的方差模型确定与所确定的回归模型有关的置信区域并将所述置信区域显示在所述接口装置上的步骤(a5);
·所述置信区域由公式为的上边界和公式为的下边界限定,其中,f表示所述回归模型,并且g表示所述方差模型;
·所述异常分数(Zscore)通过公式获取,其中,resexec是与所述物理参数的测量值有关的所述残差,是所述方差的估计值,并且mean是对于所述多个对的所述平均残差值;
·所述装置是飞行器发动机;
·所述物理参数从与所述发动机有关的压强、内部温度、流体流动速率、噪声水平和燃料密度中选择;
·所述工作参数从与所述发动机有关的速率、燃料流动速率、燃料温度、环境压强和环境温度中选择;
·所述发动机被置于包括传感器的测试台上,所述步骤(e)包括如果告知所述测量值异常则将所述测试台停机;
·步骤(e)包括如果未告知所述测量值异常则将由所述物理参数的所测量的值和所述工作参数的值形成的对添加到对的所述数据库中。
根据第二方面,本发明涉及一种用于对针对所述装置的工作参数的值的由传感器测量的装置的物理参数的值是否正常进行估计的设备,包括数据处理装置、数据存储装置和接口装置,所述数据存储装置件多个对存储在数据库中,每个对限定对于所述工作参数的值xi的所述物理参数的参考值,所述设备的特征在于,所述数据处理装置被配置为实施:
-用于基于与所述多个对有关的回归模型计算对于所述工作参数的值的所述物理参数的估计值的模块;
-用于计算相关的残差的模块;
-用于基于与所述多个对有关的方差模型计算对于所述工作参数的值的所述物理参数的方差的估计值的模块;
-用于基于所述残差、所述方差的估计值和对于所述多个对的平均残差值计算所测量的值的异常分数的模块;
-用于将所测量的值的所述异常分数与标准偏差阈值相比较的模块;
-用于若所述异常分数高于所述阈值则在接口装置上发送告警信号以告知所述测量值为异常的模块。
根据其他有利且非限制性特征:
·所述数据处理模块(31)还被配置为实施:
-用于通过基于所述工作参数的值x并根据所述(xi;yi)i∈【1,n】对的集对所述物理参数的值y回归建模来确定与所述多个对有关的所述回归模型的模块,其中,xi;yi表明存储在数据库中的对的值;
-用于对于所述对中的每一个计算所述物理参数的估计值和相关的残差的模块;
-用于计算所述残差的平均值的模块;
-用于在尺寸为w的滑动窗口上计算残差方差的值的集的模块,每个值与一对的所述工作参数的值有关;
-用于通过基于所述工作参数的值x并根据所述{xj;varj}j∈【1,n-w+1】的集对所述残差方差的值var回归建模来确定与所述多个对有关的所述方差模型的模块,其中,varj表明所计算的残差方差值,并且xj表明相关的工作参数的值。
根据第三方面,本发明涉及一种***,该***包括:
-测试台,包括传感器并且适用于容纳装置;
-根据第二方面所述的用于对针对所述装置的工作参数的值的由传感器测量的所述装置的物理参数的值是否正常进行估计的设备。
根据第四和第五方面,本发明涉及一种计算机程序产品和可由计算设备读取的存储装置。该计算机程序产品包括代码指令,该代码指令用于实现根据本发明的第一方面的用于对对针对所述装置的工作参数的值的由传感器测量的装置的物理参数的值是否正常进行估计的方法。所述存储装置上的计算机程序产品包括代码指令,所述代码指令用于实现根据第一方面所述的用于对针对所述装置的工作参数的值的由传感器测量的装置的物理参数的值是否正常进行估计的方法。
附图说明
通过阅读依照一优选实施例的说明,本发明的进一步地特征和优点将变得明显。该说明参照附图给出,在附图中:
图1表示实施根据本发明的方法的一示例性环境;
图2a-2b示出了根据本发明的示例性方法的两个阶段的步骤;
图3a-3f表示在根据本发明的方法的各步骤过程中获取的数据接口上的显示区域。
具体实施方式
参见图1,本发明是用于对针对发动机1的工作参数的值xexec的由受控装置1的物理参数的传感器20测量的值yexec是否正常进行估计的方法。受控装置1特别地是飞行器发动机1,通常是涡轮机的整机或部分,特别地是双流涡轮机。
尽管本说明是针对(在飞行或测试台中的)涡轮机的一个优选示例进行的,然而可以理解,本发明并不仅限于此,并且“装置(device)”指的是接收至少一个工作参数作为“输入”的任何物理***,有必要通过传感器对该至少一个工作参数进行监视,并且直接对该涡轮机的行为进行建模是很复杂的。测试机器这一含义是当装置在例如测试台上进行测试(而非实际使用)时赋予的。
本方法因此能够由任何推进装置来实施,该推进装置是或不是陆地的。
除了推进***之外,将很容易理解,本发明能够具有航空方面的其他应用,特别是起落架***,或者任何其他需要对本身工作进行监视的复杂设备装置。
在另一领域中,例如在铁路领域中,能够对火车故障***或车辆的其他***件进行监视。
本发明还涵盖其他的制造或实验***。
本方法可应用于任何的测量监视(包括在装置的使用寿命内),但是优选地,测试测量的示例如下:发动机1(或另一装置)被置于连接有传感器20的测试台2中。测试台2用于模拟发动机1在实际条件下的工作。本发明的目的是验证在测试过程中得到的测量值是否有效。因此,测量值是否正常与该测量值的有效性直接相关。声明为无效的测量值被认为是异常的。
工作参数(值x)是“解释性”或“预测性”变量,这与测量的物理参数(值y)相对,测量的物理参数是“用于解释”或“用于预测”的变量。换言之,工作参数的值x是原因,而物理参数的值y是结果。
更准确地,工作参数是与使用者能够作用或者环境具有影响的发动机1有关的受控值。换言之,工作参数是输入命令。在本说明的下文中,发动机速率(即发动机1的转子每单位时间进行的旋转数量,当时间单位为分钟时经常被直接称为“发动机rpm(每分钟转速)”)将被采用,但是将会理解,能够监视诸如注入发动机1的燃料流动速率、注入发动机1的燃料温度、发动机1周围的环境压强以及发动机1周围的环境温度之类的许多其他的工作参数。该参数是针对发动机选择的输入参数。
物理参数代表发动机1响应于被施加工作参数的预期行为的特征物理量,针对该特征物理量来适配传感器20。可以理解,可以提供多个适配于各物理量的传感器20。特别地,该物理参数能够是发动机1的一点处的压强、发动机1的一点处的内部温度、发动机1的一点处的流体流动速率、发动机1产生的噪声水平、发动机1的燃料密度等。这些本领域的技术人员将根据测试目的来选择要测量的物理量的类型。在本说明的下文中,压强传感器20测量的压强将被作为示例采用。
(x;y)对指明点采集:对于工作参数的值x,物理参数的值y由传感器20测量。
数据库(该数据库构成前述的测试数据库)中存储有多个(xi;yi)i∈【1,n】对,每一个对限定了针对工作参数的值xi的物理参数的参考值yi,该数据库自身存储在数据存储装置32中。“参考”值意味着这些值是在稳定后的测试阶段中被获取的,即能够被认为是正常的。可能的异常值已经从库中移除。
数据存储装置32(通常为硬盘)和数据处理装置31(例如处理器)是设备3中连接到传感器20(如图1中所示)处或直接整合到测试台2中的设备。设备3(或另一设备)配备有接口装置(诸如键盘和屏幕)以与数据库交互,并且特别用于显示结果(参见下文)。
(xexec;yexec)对指代“被监视”的测量,即试图针对该测量值来确定测量值yexec是否正常。
监视的测量值能够是实时(特别是在发动机的使用寿命中工作过程中)获取的测量值或甚至是已经存储在数据库中并待用的延时测量值(xexec;yexec)对。
用于实时查询并分析测试数据库的***能够由数据处理装置31实施。
学习阶段
本方法包括两个阶段。第一阶段为学习阶段,且第二阶段为运行阶段。优选地,预先实施学习阶段以便建立模型,该模型将在下文中描述(并且可能将模型存储在数据存储装置32上);并且之后在每个新的测量值采集中实施运行阶段。运行阶段使得可以对测量的值yexec是否正常进行估计,这是本发明的目的。学习阶段能够时不时地恢复以更新模型。
替代性地,很有可能不先行学习,并且可能在每次实施运行阶段时实时确定模型。
本说明中的下文中将对两个阶段进行说明。
学习阶段能够被视为仅对库中的数据进行处理的步骤集(即独立于(xexec;yexec)对)。
参见图2a(使用速率/压强参数的示例),学习阶段始于步骤(a0):根据(xi;yi)i∈【1,n】对的值的集将y回归建模为x的函数来确定与所述多个(xi;yi)i∈【1,n】对有关的回归模型。该回归模型将被用在运行阶段中。
回归指代本领域的技术人员公知的分析一个变量(本申请中为y)与一个或更多个其他变量(本申请中为x)的关系的一套统计方法。换言之,步骤(a0)在于确定函数f,该函数f使得对于给定的关联类型,基于值xi的值yi最佳近似(apporoximer au mieux)。关联类型同样可知为线性、多项式、指数、对数、回归等。
所使用的关联类型有利地根据曲线形状来选择,并且能够通过使确定系数最大化的优化(例如专利申请FR2939928中所述的方式)来自动完成。
图3a示出了由n个(xi;yi)i∈【1,n】对形成的散点图以及通过回归获取的模型。
在第二步骤(a1)中,学***均值:
一旦建立回归模型,数据处理装置31确定在运行阶段中使用的另一模型:方差模型。
同样与所述多个(xi;yi)i∈【1,n】对有关的该模型通常在两个连续步骤中进行计算。在步骤(a3)中,对残差方差值(varj)j∈【1,n-w+1】集进行计算,该组值将使得步骤(a4)中的回归能够实施以获取方差模型。
通过已知方式,方差被作为平方残差的平均值来计算。利用尺寸为w的“滑动窗口”(w是算法的预定参数,优选地,选择为具有足够的尺寸,以避免置信管(tube deconfiance)扩大,参见下文。举例来说,能够采用样本数量的10%。)来获取多个方差值。更准确地,对于每个方差值varj计算(xi;yi)i∈【1,n】对集中的子集(指数j至j+w-1)。特别地,
类似于步骤(a0),在步骤(a4)中通过根据x对var回归建模来确定与所述多个(xi;yi)i∈【1,n】对有关的方差模型。换言之,在此步骤中,,数据处理装置31确定函数g以使得对于j∈【1,n-w+1】有
图3c示出了由n-w+1个{xj;varj}j∈【1,n-w+1】对形成的散点图和通过回归获取的模型。
在可选步骤(a5)中,学习阶段能够包括对与回归模型有关的“置信管”(即置信区域)进行限定。该管包括上边界和下边界。在这两个边界之间,关于(x;y)对的点被认为是正常的并且不在外部,正如在下文中能够看出的。
置信管由相对于回归模型的幅度偏差或者标准偏差阈值σ来限定(参见下文)。如之前的解释,该管限定了作为方差的函数的置信区间,该置信区间增强结果的真实性:假否和假真的数量大幅减少。
上边界的等式由公式给定,并且下边界的等式由公式给定。
置信管优选地被添加到由n个(xi;yi)i∈【1,n】对形成的散点图和由(图3a中的类型的)回归获取的模型的图示中,以便生成图3d。如下所述,该管使得能够对运行阶段的结果进行预测,并且以非常直观的方式(且以对于非本领域的技术人员容易理解的方式)说明测量值是否正常。
将会注意到,点(xi;yi)i∈【1,n】的分布随参数x并不是均等的。在工作参数为速率的示例中,位于高速率的测量值事实上远多于处于低速率的测量值。这样的不均等分布能够引起如果参数w的选择不佳,则置信管在其开始处扩大,如图3e中所示。
优选地,参数w被选得足够大,以使得函数g递增。为此,步骤(a3)和之后的步骤(a4)能够随着w值的增大迭代重复直至满足该条件。
运行阶段
如前文中的解释,学习阶段代表用于加快运行阶段(对应于根据本发明的本方法的核心)的准备工作。学习阶段能够与运行阶段“同时”迭代进行。在这一部分的说明中,将引用全部之前所述的相关公式。
该阶段使得可以对针对工作参数的值xexec的物理参数的测量值yexec是否正常进行估计。
该阶段由图2b示出。如果学习阶段预先实施,则如该图中所示,能够从数据存储装置32中加载模型。
在步骤(a)中,物理参数的估计值针对工作参数的值xexec基于回归模型来计算,该回归模型与所述多个(xi;yi)i∈【1,n】有关(该模型可以在学习阶段确定)。该计算如前文中由公式解释。
在接下来的步骤(b)中,相关残差Resexec的计算为:
在步骤(c)(该步骤可以与步骤(a)和步骤(b)中的一个或两个一起执行)中,物理参数的方差的估计值针对工作参数的值xexec基于方差模型来计算,该方差模型与所述多个(xi;yi)i∈【1,n】对有关(该模型可以在学习阶段确定)。该计算如前文中由公式解释。
针对所述多个(xi;yi)i∈【1,n】对的残差、方差的估计值和残差平均值mean(可选地在学习阶段中计算)使得数据处理装置能够在步骤(d)中计算测量值yexec的异常分数Zscore(也称为置信分数)。优选地,,异常分数由公式给定。该分数越高,测量值越有可能异常。
在步骤(e)中,测量值yexec的异常分数(Zscore)与先前所述的阈值σ相比较(由标准偏差的数量表达的阈值,例如为三个至六个标准偏差)。
如果异常分数Zscore大于所述阈值σ,则在步骤(f)中在接口33上告知该测量值异常。如果这是在测试台2上的实时测试,则能够触发告警,并且测试台2被停止(当前测试无效,并且在测试台被重新使用之前应当对异常进行分析)。(xexec;yexec)对则不被添加到数据库中(或者在等待检查时被删除)。
在相反的情况下,测量值被认为是正常的,并且(xexec;yexec)对联合数据库的参考值。有可能恢复学习阶段以更新模型。
优选地,(例如用于研究部门的)综合报告则自动下发。
将要注意的是,异常测量值将被表示为在前述置信管之外,例如如图3f中所示。
特别地,如果Zscore>σ,则(等于置信管的直径的一半,这意味着相关的点在置信管之外。)
事实上,残差符合平均值为mean的正态律。测试(与Zscore相比较)在于确定当前的观测是否有可能来自相同的律。
如果绘制出正态律的分布,能够看出,超过六个标准偏差时,位置在分布尾部。观察到该标准偏差的正常事件的概率的量级为10-9。
设备和***
用于实施之前所述(对针对装置1的工作参数的值xexec的由传感器20测量的装置1的物理参数的值yexec是否正常进行估计的)方法的设备3(图1中所示)包括数据处理装置31、数据存储装置32和接口装置33。
数据存储装置32在数据库中存储多个(xi;yi)i∈【1,n】对,每个(xi;yi)i∈【1,n】限定对于工作参数的值xi的物理参数的参考值yi。
数据处理装置31配置为实现:
-基于与所述多个(xi;yi)i∈【1,n】对有关的回归模型计算针对工作参数的值xexec的物理参数的估计值的模块;
-计算有关的残差resexec的模块;
-基于与所述多个(xi;yi)i∈【1,n】对有关的方差模型计算针对工作参数的值xexec的物理参数的方差的估计值的模块;
-基于残差resexec、估计方差值和所述多个(xi;yi)i∈【1,n】对的平均残差值计算测量值yexec的异常分数Zscore(Z分数)的模块;
-将测量值yexec的异常分数Zscore与标准偏差阈值σ相比较的模块;
-若异常分数Zscore高于所述阈值σ则在接口装置33上发送告警信号以告知测量值异常的模块。
如果设备3还实施学习阶段,则数据处理装置31还配置为实现:
-通过基于{xi;yi}i∈【1,n】集中的值x对值y回归建模来确定与所述多个(xi;yi)i∈【1,n】对有关的所述回归模型的模块(该模块可以配置为在数据装置32中存储回归模型以待将来使用);
-针对每一对计算物理参数的估计值和彼此相关的残差resi的模块;
-在尺寸为w的滑动窗口上计算残差方差的(varj)j∈【1,n-w+1】值集的模块;
-通过基于{xj;varj}j∈【1,n-w+1】集中的x对var回归建模来确定与所述多个(xi;yi)i∈【1,n】对有关的所述方差模型的模块(该模块可以配置为在数据装置32中存储方差模型以待将来使用)。
设备3优选地被解释为包括在***中,该***还包括测试台2,测试台2包括对物理参数的值yexec进行测量的传感器20并且适于容纳装置1。
计算机程序产品
根据第四和第五方面,本发明涉及一种计算机程序产品和可由计算设备读取的存储装置。该计算机程序产品包括代码指令,该代码指令用于(在数据处理装置31、特别是设备3的数据处理装置上)实现根据本发明的第一方面的用于对针对所述装置1的工作参数的值xexec的由传感器20测量的装置1的物理参数的值yexec是否正常进行估计的方法。该存储装置例如是该设备3的数据存储装置32,并且该存储装置上具有该计算机程序产品。
Claims (14)
1.一种用于对针对装置(1)的工作参数的值xexec而言由传感器(20)测量的所述装置(1)的物理参数的值yexec是否正常进行估计的方法,存储在数据库中的多个(xi;yi)i∈【1,n】对中的每一对限定对于所述工作参数的值xi的所述物理参数的参考值yi,所述数据库存储在数据存储装置(32)上,所述方法的特征在于,包括:由数据处理装置(31)实施以下步骤:
(a)基于与所述多个(xi;yi)i∈【1,n】对有关的回归模型计算对于所述工作参数的值xexec的所述物理参数的估计值
(b)计算相关的残差resexec;
(c)基于与所述多个(xi;yi)i∈【1,n】对有关的方差模型计算对于所述工作参数的值xexec的所述物理参数的方差的估计值
(d)基于所述残差resexec、所述方差的估计值和对于所述多个(xi;yi)i∈【1,n】对的平均残差值mean计算所测量的值yexec的异常分数Zscore;
(e)将所测量的值yexec的所述异常分数Zscore与标准偏差阈值σ相比较;
(f)如果所述异常分数Zscore高于所述标准偏差阈值σ,则在接口装置(33)上告知所述测量值为异常。
2.根据权利要求1所述的方法,包括处理所述数据库的所述多个(xi;yi)i∈【1,n】对的在先阶段,该在先阶段包括由所述数据处理装置(31)实施以下步骤:
(a0)通过基于所述工作参数的值x并根据所述(xi;yi)i∈【1,n】对的集对所述物理参数的值y回归建模来确定与所述多个(xi;yi)i∈【1,n】对有关的所述回归模型,其中,xi;yi表明存储在数据库中的对的值;
(a1)对于所述(xi;yi)i∈【1,n】对中的每一个,计算所述物理参数的估计值和相关的残差resi;
(a2)计算所述残差resi的平均值mean;
(a3)在尺寸为w的滑动窗口上计算残差方差的值(varj)j∈【1,n-w+1】的集,每个值与一(xi;yi)i∈【1,n】对的所述工作参数的值xj有关;
(a4)通过基于所述工作参数的值x并根据{xj;varj}j∈【1,n-w+1】的集对所述残差方差的值var回归建模来确定与所述多个(xi;yi)i∈【1,n】对有关的所述方差模型,其中,varj指代所计算的残差方差值,并且xj指代相关的工作参数的值。
3.根据权利要求2所述的方法,其中,所述在先阶段包括基于所确定的方差模型确定与所确定的回归模型有关的置信区域并将所述置信区域显示在所述接口装置(33)上的步骤(a5)。
4.根据权利要求3所述的方法,其中,所述置信区域由公式的上边界和公式为的下边界限定,其中,f表示所述回归模型,并且g表示所述方差模型。
5.根据权利要求1至4中任一项所述的方法,其中,所述异常分数Zscore通过公式获取,其中,resexec是与所述物理参数的测量值有关的所述残差,是所述方差的估计值,并且mean是所述多个(xi;yi)i∈【1,n】对的所述平均残差值。
6.根据权利要求1至4中任一项所述的方法,其中,所述装置(1)是飞行器发动机。
7.根据权利要求6所述的方法,其中,所述物理参数从与所述发动机(1)有关的压强、内部温度、流体流动速率、噪声水平和燃料密度中选择。
8.根据权利要求6所述的方法,其中,所述工作参数从与所述发动机(1)有关的速率、燃料流动速率、燃料温度、环境压强和环境温度中选择。
9.根据权利要求8所述的方法,其中,所述发动机(1)被置于包括传感器(20)的测试台(2)上,所述步骤(f)包括如果告知所述测量值异常则将所述测试台(2)停机。
10.根据权利要求1至4中任一项所述的方法,其中,所述步骤(f)包括如果未告知所述测量值异常则将(xexec;yexec)对添加到(xi;yi)i∈【1,n】对的所述数据库中,所述(xexec;yexec)对由所述物理参数的所测量的值和所述工作参数的值形成。
11.一种用于对针对装置(1)的工作参数的值xexec而言由传感器(20)测量的所述装置(1)的物理参数的值yexec是否正常进行估计的设备(3),包括数据处理装置(31)、数据存储装置(32)和接口装置(33),所述数据存储装置将多个(xi;yi)i∈【1,n】对存储在数据库中,每个(xi;yi)i∈【1,n】对限定对于所述工作参数的值xi的所述物理参数的参考值yi,所述设备(3)的特征在于,所述数据处理装置(31)被配置为实施:
-用于基于与所述多个(xi;yi)i∈【1,n】对有关的回归模型计算对于所述工作参数的值xexec的所述物理参数的估计值的模块;
-用于计算相关的残差resexec的模块;
-用于基于与所述多个(xi;yi)i∈【1,n】对有关的方差模型计算对于所述工作参数的值xexec的所述物理参数的方差的估计值的模块;
-用于基于所述残差resexec、所述方差的估计值和对于所述多个(xi;yi)i∈【1,n】对的平均残差值mean计算所测量的值yexec的异常分数Zscore的模块;
-用于将所测量的值yexec的所述异常分数Zscore与标准偏差阈值σ相比较的模块;
-用于若所述异常分数Zscore高于所述标准偏差阈值σ则在接口装置(33)上发送告警信号以告知所述测量值为异常的模块。
12.根据权利要求11所述的设备,其中,所述数据处理装置(31)还被配置为实施:
-用于通过基于所述工作参数的值x并根据所述(xi;yi)i∈【1,n】对的集对所述物理参数的值y回归建模来确定与所述多个(xi;yi)i∈【1,n】对有关的所述回归模型的模块,其中,xi;yi指代存储在数据库中的一个对的值;
-用于对于所述(xi;yi)i∈【1,n】对中的每一个计算所述物理参数的估计值和相关的残差resi的模块;
-用于计算所述残差resi的平均值mean的模块;
-用于在尺寸为w的滑动窗口上计算残差方差的值(varj)j∈【1,n-w+1】的集的模块,每个值与一(xi;yi)i∈【1,n】对的所述工作参数的值xj有关;
-用于通过基于所述工作参数的值x并根据{xj;varj}j∈【1,n-w+1】的集对所述残差方差的值var回归建模来确定与所述多个(xi;yi)i∈【1,n】对有关的所述方差模型的模块,其中,varj指代所计算的残差方差值,并且xj指代相关的工作参数的值。
13.一种***,包括:
-测试台(2),包括传感器(20)并且适用于容纳装置(1);
-根据权利要求11和12中任一项所述的用于对针对所述装置(1)的工作参数的值xexec而言由传感器(20)测量的所述装置(1)的物理参数的值yexec是否正常进行估计的设备(3)。
14.一种存储装置,可由计算设备读取,所述存储装置上的计算机程序产品包括代码指令,所述代码指令用于实现根据权利要求1至9中任一项所述的用于对针对所述装置(1)的工作参数的值xexec而言由传感器(20)测量的装置(1)的物理参数的值yexec是否正常进行估计的方法。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR1452650 | 2014-03-27 | ||
FR1452650A FR3019295B1 (fr) | 2014-03-27 | 2014-03-27 | Procede d'estimation du caractere normal ou non d'une valeur mesuree d'un parametre physique d'un moteur d'aeronef |
PCT/FR2015/050785 WO2015145085A1 (fr) | 2014-03-27 | 2015-03-26 | Procédé d'estimation du caractère normal ou non d'une valeur mesurée d'un paramètre physique d'un moteur d'aéronef |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106233115A CN106233115A (zh) | 2016-12-14 |
CN106233115B true CN106233115B (zh) | 2018-02-23 |
Family
ID=50780782
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201580016936.6A Active CN106233115B (zh) | 2014-03-27 | 2015-03-26 | 估计飞行器发动机的物理参数的测量值是否正常的方法 |
Country Status (8)
Country | Link |
---|---|
US (1) | US10060831B2 (zh) |
EP (1) | EP3123139B1 (zh) |
CN (1) | CN106233115B (zh) |
BR (1) | BR112016022187B1 (zh) |
CA (1) | CA2943397C (zh) |
FR (1) | FR3019295B1 (zh) |
RU (1) | RU2686252C2 (zh) |
WO (1) | WO2015145085A1 (zh) |
Families Citing this family (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102016202370A1 (de) * | 2016-02-17 | 2017-08-17 | MTU Aero Engines AG | Verfahren zum Ermitteln eines Einflusses eines Innen-Prüfstands auf eine im Innen-Prüfstand betriebene Gasturbine |
CN106053087A (zh) * | 2016-05-12 | 2016-10-26 | 大工(青岛)新能源材料技术研究院有限公司 | 一种燃气内燃机测评方法 |
CN106200621A (zh) * | 2016-08-17 | 2016-12-07 | 中国航空工业集团公司西安飞行自动控制研究所 | 一种飞控***正常工作判据的获取方法 |
SE541872C2 (en) * | 2016-09-22 | 2020-01-02 | Scania Cv Ab | Method and system for predicting the fuel consumption for a vehicle |
US10642852B2 (en) | 2016-09-26 | 2020-05-05 | Splunk Inc. | Storing and querying metrics data |
US10587635B2 (en) * | 2017-03-31 | 2020-03-10 | The Boeing Company | On-board networked anomaly detection (ONAD) modules |
FR3074590B1 (fr) * | 2017-12-04 | 2023-03-17 | Soc Air France | Methode de prediction d'une anomalie de fonctionnement d'un ou plusieurs equipements d'un ensemble |
FR3075359B1 (fr) * | 2017-12-18 | 2019-11-22 | Safran Aircraft Engines | Detection et accomodation de circuits ouverts intermittents sur un capteur moteur d'une turbomachine d'aeronef |
CN108734371A (zh) * | 2018-02-12 | 2018-11-02 | 阿里巴巴集团控股有限公司 | 一种针对风控指令的处理方法、装置及设备 |
CN108563739B (zh) * | 2018-04-11 | 2021-05-04 | 平安科技(深圳)有限公司 | 天气数据获取方法及装置、计算机装置及可读存储介质 |
FR3095424A1 (fr) * | 2019-04-23 | 2020-10-30 | Safran | Système et procédé de surveillance d’un moteur d’aéronef |
CN112580903A (zh) * | 2019-09-27 | 2021-03-30 | 华晨宝马汽车有限公司 | 用于评估发动机质量稳定性的方法和装置及存储介质 |
US11615656B2 (en) * | 2019-12-20 | 2023-03-28 | Pratt & Whitney Canada Corp. | Method and system for diagnosing an engine or an aircraft |
CN111967508A (zh) * | 2020-07-31 | 2020-11-20 | 复旦大学 | 基于显著图的时间序列异常点检测方法 |
CN112179522A (zh) * | 2020-09-14 | 2021-01-05 | 中国空气动力研究与发展中心超高速空气动力研究所 | 一种飞行器表面温度传感器测试***及测试方法 |
CN112378670B (zh) * | 2020-11-10 | 2021-10-15 | 北京航空航天大学 | 一种基于改进粒子滤波的火箭发动机故障检测方法 |
CN116391131A (zh) * | 2020-12-17 | 2023-07-04 | 浙江吉利控股集团有限公司 | 一种电池温度预测方法及*** |
Family Cites Families (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3593886B2 (ja) * | 1998-06-17 | 2004-11-24 | 日産自動車株式会社 | 車両用能動型振動制御装置 |
RU2282940C2 (ru) * | 2001-12-03 | 2006-08-27 | Юнайтид Текнолоджиз Копэрейшн | Комплексная система для контроля за состоянием изделия из удаленного пункта |
US7464533B2 (en) * | 2003-01-28 | 2008-12-16 | General Electric Company | Apparatus for operating gas turbine engines |
US7072797B2 (en) * | 2003-08-29 | 2006-07-04 | Honeywell International, Inc. | Trending system and method using monotonic regression |
US7062370B2 (en) * | 2004-03-30 | 2006-06-13 | Honeywell International Inc. | Model-based detection, diagnosis of turbine engine faults |
DE102005054735B4 (de) * | 2005-11-17 | 2019-07-04 | Robert Bosch Gmbh | Verfahren und Vorrichtung zum Betreiben einer Brennkraftmaschine |
US20090030752A1 (en) * | 2007-07-27 | 2009-01-29 | General Electric Company | Fleet anomaly detection method |
US20090043447A1 (en) * | 2007-08-07 | 2009-02-12 | General Electric Company | Systems and Methods for Model-Based Sensor Fault Detection and Isolation |
US20090112519A1 (en) * | 2007-10-31 | 2009-04-30 | United Technologies Corporation | Foreign object/domestic object damage assessment |
JP2010025024A (ja) * | 2008-07-22 | 2010-02-04 | Toyota Motor Corp | 温度センサの異常検出装置 |
FR2939928B1 (fr) | 2008-12-15 | 2012-08-03 | Snecma | Standardisation de donnees utilisees pour la surveillance d'un moteur d'aeronef |
DE102009002682B4 (de) * | 2009-04-28 | 2012-09-20 | Airbus Operations Gmbh | Einrichtung und Verfahren zur Residuenauswertung eines Residuums zur Erkennung von Systemfehlern im Systemverhalten eines Systems eines Flugzeugs |
FR2965915B1 (fr) | 2010-10-11 | 2013-08-02 | Snecma | Systeme de surveillance d'un banc d'essai de moteur d'aeronef |
US8219356B2 (en) * | 2010-12-23 | 2012-07-10 | General Electric Company | System and method for detecting anomalies in wind turbines |
FR2983528B1 (fr) * | 2011-12-05 | 2014-01-17 | Snecma | Methode de surveillance d'une chaine de mesure d'un turboreacteur |
-
2014
- 2014-03-27 FR FR1452650A patent/FR3019295B1/fr active Active
-
2015
- 2015-03-26 WO PCT/FR2015/050785 patent/WO2015145085A1/fr active Application Filing
- 2015-03-26 RU RU2016142123A patent/RU2686252C2/ru active
- 2015-03-26 CA CA2943397A patent/CA2943397C/fr active Active
- 2015-03-26 US US15/129,279 patent/US10060831B2/en active Active
- 2015-03-26 CN CN201580016936.6A patent/CN106233115B/zh active Active
- 2015-03-26 BR BR112016022187-7A patent/BR112016022187B1/pt active IP Right Grant
- 2015-03-26 EP EP15717037.4A patent/EP3123139B1/fr active Active
Also Published As
Publication number | Publication date |
---|---|
EP3123139B1 (fr) | 2018-05-23 |
BR112016022187A2 (zh) | 2017-08-15 |
EP3123139A1 (fr) | 2017-02-01 |
CA2943397A1 (fr) | 2015-10-01 |
US20170176292A1 (en) | 2017-06-22 |
US10060831B2 (en) | 2018-08-28 |
RU2016142123A3 (zh) | 2018-09-24 |
RU2686252C2 (ru) | 2019-04-24 |
CA2943397C (fr) | 2019-10-15 |
FR3019295B1 (fr) | 2016-03-18 |
WO2015145085A1 (fr) | 2015-10-01 |
FR3019295A1 (fr) | 2015-10-02 |
BR112016022187B1 (pt) | 2022-09-13 |
RU2016142123A (ru) | 2018-04-28 |
CN106233115A (zh) | 2016-12-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106233115B (zh) | 估计飞行器发动机的物理参数的测量值是否正常的方法 | |
EP3441947B1 (en) | System and method for remaining useful life determination | |
JP7340265B2 (ja) | 異常検出装置、異常検出方法、およびプログラム | |
CN112508105B (zh) | 一种采油机故障检测与检索方法 | |
US9779495B2 (en) | Anomaly diagnosis method and apparatus | |
EP2384971B1 (en) | Method of determining a maneuver performed by an aircraft | |
US11061390B2 (en) | System fault isolation and ambiguity resolution | |
US20130274964A1 (en) | Flight data monitoring and validation | |
CN110414154B (zh) | 一种带有双测点的风机部件温度异常检测和报警方法 | |
US20160369777A1 (en) | System and method for detecting anomaly conditions of sensor attached devices | |
US20090276136A1 (en) | Method for calculating confidence on prediction in fault diagnosis systems | |
US20170331844A1 (en) | Systems and methods for assessing airframe health | |
US20200143292A1 (en) | Signature enhancement for deviation measurement-based classification of a detected anomaly in an industrial asset | |
CN108896299A (zh) | 一种齿轮箱故障检测方法 | |
JP7468724B2 (ja) | 配管診断装置、配管診断方法、及びプログラム | |
US9607451B2 (en) | Method and a system for merging health indicators of a device | |
CN107077134A (zh) | 用于对监控飞机发动机的***进行验证的工具 | |
CN103775832A (zh) | 基于瞬变流反问题方法的输油管道漏失检测的装置 | |
KR20170053692A (ko) | 커널 회귀 모델의 앙상블을 위한 장치 및 방법 | |
JP6714498B2 (ja) | 設備診断装置及び設備診断方法 | |
CN103217172A (zh) | 一种卡尔曼滤波传感器信息融合的故障检测方法 | |
US7366639B2 (en) | Methods for establishing alerts and/or alert limits for monitoring mechanical devices | |
WO2018155536A1 (ja) | 位置決定装置、それを備えた位置決定システム、及び位置決定方法並びに位置決定プログラム | |
CN106528940B (zh) | 基于数理模型对列车车轴性能进行评价和预警的方法 | |
Dzakowic et al. | Advanced Techniques for the verification and validation of prognostics & health management capabilities |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | 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 |