CN102985884B - 用于更新模型预测控制器中模型的方法和系统 - Google Patents

用于更新模型预测控制器中模型的方法和系统 Download PDF

Info

Publication number
CN102985884B
CN102985884B CN201180019994.6A CN201180019994A CN102985884B CN 102985884 B CN102985884 B CN 102985884B CN 201180019994 A CN201180019994 A CN 201180019994A CN 102985884 B CN102985884 B CN 102985884B
Authority
CN
China
Prior art keywords
model
mpc
processing factory
mpm
performance grade
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.)
Expired - Fee Related
Application number
CN201180019994.6A
Other languages
English (en)
Other versions
CN102985884A (zh
Inventor
A.巴德维
N.库巴
S.布哈特
N.N.南多拉
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
ABB Schweiz AG
Original Assignee
ABB Research Ltd Switzerland
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by ABB Research Ltd Switzerland filed Critical ABB Research Ltd Switzerland
Publication of CN102985884A publication Critical patent/CN102985884A/zh
Application granted granted Critical
Publication of CN102985884B publication Critical patent/CN102985884B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B13/00Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
    • G05B13/02Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
    • G05B13/04Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
    • G05B13/048Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators using a predictor
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B13/00Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
    • G05B13/02Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
    • G05B13/04Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
    • G05B13/041Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a variable is automatically adjusted to optimise the performance
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B13/00Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
    • G05B13/02Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
    • G05B13/04Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
    • G05B13/042Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a parameter or coefficient is automatically adjusted to optimise the performance

Landscapes

  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Evolutionary Computation (AREA)
  • Medical Informatics (AREA)
  • Software Systems (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Feedback Control In General (AREA)
  • Testing And Monitoring For Control Systems (AREA)

Abstract

本发明涉及用于更新模型预测控制器中模型的方法。该方法包括评估加工厂的操作性能等级与期望性能等级的偏差。通过更新模型预测控制器中的模型针对模型工厂不匹配诊断模型预测控制。诊断模型预测控制器包括确定与模型工厂不匹配相关的模型预测误差,量化模型工厂不匹配,并更新所述模型预测控制器中的模型。本发明还涉及用于根据本发明方法更新模型预测控制器中模型的系统。

Description

用于更新模型预测控制器中模型的方法和系统
技术领域
本发明涉及模型预测控制(MPC),并且更具体地说,涉及用于更新模型预测控制器中模型的方法和系统。
背景技术
在过程控制工业中,采用先进过程控制(APC)以便降低操作成本,获得高生产率,始终如一地保持质量,还有其它类似原因。APC允许从过程控制工业的目前操作方案变换到改进且更大生产率的操作方案,并且还适应过程控制工业中涉及的过程的操作和设计约束。
一般而言,多变量APC在多变量预测控制器中实现了更流行的先进多变量控制方案,称为模型预测控制(MPC)。MPC使用加工厂中涉及的过程的数学模型,以便预测过程的未来动态行为,从而提供其工厂的过程和操作的最优操控变量。由此,可以理解,模型的准确性是有效且成功实现MPC的关键要素。
实际上,工厂动态特性改变,从而引起模型与工厂之间的不匹配。该不匹配称为模型工厂不匹配(MPM)。MPM导致工厂动态特性的不准确预测。使用模型受所述MPM影响的APC将降级MPC及其总体控制性能,这因此改变产品质量并引起经济损失。
为了消除控制器的性能降级,在检测到差的控制器性能时更新模型变得很重要。可通过良好确立的MPC性能监视来检测差的控制器性能。例如,简单方法可以是分析预测误差,预测误差是模型预测与真实输出之间的差。在检测到差的控制器性能之后,可识别并使用确立的诊断技术诊断造成差的控制器性能的原因,诸如差的模型(即MPM)、未测量的干扰以及约束饱和等。
当前,在识别或检测到差的模型时,MPM由模型诊断,并在重新识别模型之后诊断。重新识别模型涉及设计扰动信号、判定和/或考虑扰动期间工厂的操作条件、选择适当模型并估计模型参数。这要求专业程度高并且相当耗时。还有,更重要的是,它涉及较长的扰动周期,由此在扰动周期期间制造大量或大数量质量低的产品,通常称为不合格产品。
因而,存在对于显著降低扰动周期并更有效地降低和/或消除MPM的替换方法的需要。
发明目的
本发明的一个目的是降低诊断MPM期间的扰动周期以便降低和/或消除MPM。
本发明的另一目的是更有效地降低MPM。
发明内容
本发明旨在提供用于更新加工厂的模型预测控制器中模型的方法和系统,以便降低和/或消除由MPM引起的、加工厂的操作性能等级与期望性能等级的偏差。
相应地,本发明提供了用于模型预测控制器中的模型更新的方法。该方法包括评估加工厂的操作性能等级与期望性能等级的偏差的步骤。通过更新模型预测控制器中的模型及其MPC针对模型工厂不匹配(MPM)诊断模型预测控制(MPC)。诊断所述MPC包括确定与MPM相关的模型预测误差。然后,量化MPM并更新所述模型预测控制器中的模型。
相应地,本发明还提供了用于根据本发明方法更新模型预测控制器中模型的系统。该系统包括用于监视MPC性能的性能监视单元、用于检测造成加工厂的操作性能等级与期望性能等级的偏差的原因的检测单元以及用于生成特征上仅足以计算m+1个非零滞后相关系数的白噪声信号的白噪声生成单元,其中m是MPC的操控变量数。提供估计器单元用于通过计算加工厂的控制器的输出值与MPC的现有模型的输出值之间的差来估计模型预测误差。还有,提供量化器单元用于量化MPM。另外,该系统还具有更新器单元用于基于量化的MPM提供更新模型,并用于更新模型预测控制器中加工厂的模型。
附图说明
参考附图,附图中:
图1描绘了加工厂中的典型模型预测控制设置;以及
图2示出了根据本发明的模型预测控制器中的模型更新的框架。
具体实施方式
参考图1和2作为非详尽性示范实施例的例证描述本发明。
在图1中,加工厂(100)由分布式控制系统DCS(110)中的初级等级控制器、诸如PID控制器等控制。提供了模型预测控制(120),模型预测控制(120)使用加工厂中涉及的过程的数学模型,下文称为模型,以便预测过程的未来动态行为,相应地提供其工厂的过程和操作的最优操控变量。提供性能监视单元(130)用于监视MPC(120)的性能。另外,检测单元(140)检测造成加工厂的操作性能等级与期望性能等级的偏差的原因。
设想本文上面提到的偏差已经由作为模型与加工厂之间的不匹配的模型工厂不匹配(MPM)造成或由于模型工厂不匹配(MPM)而出现,还相对于加工厂的开环条件进一步说明本发明,该设想本质上不是约束性的,而纯粹是示范性的,以便更好地理解。
针对MPM诊断MPC变成强制性的。本文进一步参考图2来描述本发明。本文提到的诊断还包含相对于校正所述偏差及其MPM来实现MPC(220)的模型(240)中的必要校正。
相应地,离线采取MPC(220),其中向控制器的设置点值是具有预定值的恒定值,或者是紧接在离线采取MPC(220)之前从MPC(220)演进的。在后一情况下MPC(220)被视为在线,其中控制器的设置点是MPC(220)的操控变量。提供白噪声生成单元(230)用于生成白噪声信号,所述白噪声信号在特征上仅足以计算模型预测误差与MPC的操控变量之间的(m+1个,其中m是MPC的操控变量数)非零滞后相关系数。白噪声信号被添加到设置点,如下式所示:
这里,SPi(t)是在时刻t当MPC处于“离线”模式时的第i个设置点值,Wi(t)是在时刻t的第i个白噪声序列中的值,并且(SPi)DCS,t是要到DCS(210)的第i个设置点值。应该注意,在MPC(220)“在线”的情况下,SPi(t)由MPC(220)计算,并且因此是MPC(220)的操控变量。再者,在MPC(220)“在线”模式下,
(SPi)DCS,t值还被传到现有模型(240),现有模型(240)的输出由表示,它们是受控变量的预测值。来自工厂(200)的受控变量的实际值由y表示。可通过在线测量或通过其它适当手段、诸如定期实验室分析获得来自工厂(200)的受控变量值。测量值y与模型(240)预测的值之间的差称为模型预测误差或模型残差,并且由表示。注意,在图2中,为了标示简单,(SPi)DCS,t表示。
量化器单元(260)计算之间的滞后相关系数。这些系数的值用于计算现有模型与当前工厂之间的“间隙”,即量化MPM。
考虑加工厂的单输入和单输出开环条件,为了理解上简单,说明增益不匹配的如下情况。这绝不是约束性的,而是纯粹示范性的,并且相对本发明非详尽的。类似地,本发明也适用于多输入和多输出以及闭环条件,并且可共同扩展地应用。
对于第一阶时延系统,该多项式表示如下:
其中y(k)是工厂的输出;是模型的输出;且e(k)是模型误差;
a和b是工厂的参数;am和bm是模型的参数;
k是采样时刻;
tdp是工厂的时延;tdm是模型的时延。
仅考虑增益不匹配(即a=am且tdP=tdm):
要注意,即使在所有三个参数、即a、b和td中都存在不匹配,也可应用本发明的方法。为了简化原因,仅考虑一个参数、即增益或b和bm中的不匹配。
现在,给出在两个时间序列x1与x2之间的滞后m的相关系数:
然后,根据等式1,
在滞后0,
在滞后1,
在滞后tdp+1,
在滞后tdp+2,
在滞后tdp+n,
因而,在滞后tdp+1观察到第一非零相关系数。此外,该相关显示出指数式衰减。
现在,假设是在滞后tdp+1在之间观察到的相关系数。我们然后可写:
这里,b-bm是量化的MPM。
更新器单元(270)基于量化的MPM将现有模型(240)更新成更新模型(280),并更新模型预测控制器中MPC(220)的模型。MPC(220)然后切换到在线模式。
本发明还可扩展成以与本文之前在说明书中所描述的类似方式更新子模型。由于MPM而需要更新的子模型的识别可通过已知技术(诸如应用部分相关分析的技术或任何其它适当技术或方法)来完成。这样,在子模型级更新模型还减少了更新MPM的模型的工作,并变得更有效。
因而,视情况而定,本发明可应用于更新整个模型或子模型。
此说明书中描述的示范实施例相对于本发明的范围实质上不是约束性的。将满足相同功能性和方法并且在本说明书中未明确提到或阐述的其它修改、改变、变化被视为也在本发明的范围内。在本发明范围的意义内,使用单数术语将包含其复数形式,并且反之亦然。
本发明在化学和石油化学、水泥、制浆和纸以及医药工业中发现了广泛用途,这里举出了几个示例。一些其它应用包含发电、供应链管理和行为健康等。

Claims (7)

1.一种用于更新模型预测控制器中模型的方法,所述方法包括:
a.评估加工厂的操作性能等级与期望性能等级的偏差;
b.针对模型工厂不匹配MPM诊断模型预测控制MPC,用于更新模型预测控制器中的模型及其MPC;其中诊断所述MPC包括:
a.确定与MPM相关的模型预测误差,其中确定与MPM相关的模型预测误差包括:
i.向所述加工厂中的所述控制器的设置点值添加白噪声,
ii.向所述加工厂中的所述控制器和所述MPC的现有模型提供具有添加的白噪声的设置点值,
iii.通过计算所述加工厂的输出值与所述MPC的所述现有模型的输出值之间的差来估计所述模型预测误差,
b.使用所述模型预测误差与所述MPC的每一个操控变量之间的非零滞后相关系数量化所述MPM,所述MPC的所述操控变量是所述加工厂中的所述控制器的所述具有添加的白噪声的设置点值;以及
c.基于所量化的MPM来更新所述模型预测控制器中的所述模型以便校正所述模型;
d.基于所述MPC的所校正的模型来校正加工厂的操作性能等级与期望性能等级的所述偏差。
2.如权利要求1所述的方法,其中评估所述加工厂的所述操作性能等级与所述期望性能等级的偏差包含监视所述MPC的性能;以及检测造成所述偏差的原因,所述加工厂的所述操作性能等级与所述期望性能等级的所述偏差对应于MPM。
3.如权利要求1所述的方法,其中诊断所述MPC包含通过针对所述MPM校正所述MPC的模型及其所述MPC来校正所述加工厂的所述操作性能等级与所述期望性能等级的所述偏差。
4.如权利要求1-3中任一项所述的方法,其中离线或在线执行所述模型更新。
5.如权利要求1-3中任一项所述的方法,其中所述加工厂中涉及的过程是开环或闭环。
6.如权利要求1-3中任一项所述的方法,其中所述模型包含所述加工厂的整个模型和/或其子模型。
7.一种用于根据权利要求1-6中任一项所述的方法更新模型预测控制器中模型的系统,所述系统包括:
性能监视单元,用于监视所述MPC的性能;
检测单元,用于检测所述加工厂的所述操作性能等级与所述期望性能等级的偏差的原因;
白噪声生成单元,用于生成特征上仅足以计算非零滞后相关系数的白噪声信号;
估计器单元,用于通过计算所述加工厂的输出值与所述MPC的现有模型的输出值之间的差来估计所述模型预测误差;
量化器单元,用于量化所述MPM;以及
更新器单元,用于基于量化的MPM提供更新模型并更新所述模型预测控制器中所述加工厂的模型。
CN201180019994.6A 2010-04-19 2011-04-19 用于更新模型预测控制器中模型的方法和系统 Expired - Fee Related CN102985884B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
IN1086/CHE/2010 2010-04-19
IN1086CH2010 2010-04-19
PCT/IB2011/000845 WO2011132050A1 (en) 2010-04-19 2011-04-19 A method and system for updating a model in a model predictive controller

Publications (2)

Publication Number Publication Date
CN102985884A CN102985884A (zh) 2013-03-20
CN102985884B true CN102985884B (zh) 2015-11-25

Family

ID=44120184

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201180019994.6A Expired - Fee Related CN102985884B (zh) 2010-04-19 2011-04-19 用于更新模型预测控制器中模型的方法和系统

Country Status (5)

Country Link
US (1) US9581980B2 (zh)
EP (1) EP2561411B1 (zh)
CN (1) CN102985884B (zh)
ES (1) ES2842967T3 (zh)
WO (1) WO2011132050A1 (zh)

Families Citing this family (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2013087973A1 (en) * 2011-12-16 2013-06-20 Metso Automation Oy Method of tuning a process controller
DK177915B1 (en) * 2013-05-28 2015-01-05 Core As Process control method
CN104698976B (zh) * 2014-12-23 2017-06-16 南京工业大学 一种预测控制模型性能下降的深度诊断方法
WO2017055948A1 (en) * 2015-09-28 2017-04-06 Koninklijke Philips N.V. Methods and systems for controlling gas flow using a proportional flow valve
JP6443311B2 (ja) * 2015-11-30 2018-12-26 オムロン株式会社 制御装置、制御プログラムおよび記録媒体
US10459428B2 (en) 2016-03-08 2019-10-29 Honeywell Limited Optimal closed-loop input design for identification of flat-sheet process models
US10358983B2 (en) * 2016-04-19 2019-07-23 General Electric Company Asset degradation model baselinening system and method
US11449046B2 (en) 2016-09-16 2022-09-20 Honeywell Limited Model-plant mismatch detection with support vector machine for cross-directional process behavior monitoring
US10809674B2 (en) * 2016-09-16 2020-10-20 Honeywell Limited Model-plant mismatch detection using model parameter data clustering for paper machines or other systems
US10174456B2 (en) 2016-09-23 2019-01-08 Honeywell Limited Technique to improve paper machine cross-directional model predictive control performance by creating a measurement profile reference trajectory
US10358771B2 (en) 2016-09-23 2019-07-23 Honeywell Limited Method of designing model predictive control for cross directional flat sheet manufacturing processes to guarantee spatial robustness and to prevent actuator picketing
US10309059B2 (en) * 2016-09-23 2019-06-04 Honeywell International Inc. Method of designing model predictive control for cross directional flat sheet manufacturing processes to guarantee temporal robust stability and performance
WO2020010629A1 (zh) * 2018-07-13 2020-01-16 西门子股份公司 过程控制预测模型调整方法、装置和过程控制器
EP3671373A1 (de) 2018-12-19 2020-06-24 Bystronic Laser AG Prognoseverfahren zur prognose eines bearbeitungsergebnisses für ein laserbearbeitungsverfahren
US11525375B2 (en) 2020-04-09 2022-12-13 General Electric Company Modeling and control of gas cycle power plant operation with variant control profile
CN113760625B (zh) * 2021-06-30 2024-02-02 浙江中控技术股份有限公司 一种模型预测控制器性能的评估方法及监控系统
CN115309044B (zh) * 2022-07-26 2024-05-17 福建工程学院 一种基于模型预测控制的机械臂角速度控制方法

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1742241A (zh) * 2002-12-09 2006-03-01 搭篷技术公司 具有可变动态特性的过程的自适应控制的一种系统和方法

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3578886B2 (ja) * 1997-05-02 2004-10-20 東京エレクトロン株式会社 プロセス制御システムとそのプロセスデータ転送制御方法
US7187989B2 (en) * 2003-12-22 2007-03-06 Fakhruddin T Attarwala Use of core process models in model predictive controller
US7096085B2 (en) * 2004-05-28 2006-08-22 Applied Materials Process control by distinguishing a white noise component of a process variance
US7536232B2 (en) * 2004-08-27 2009-05-19 Alstom Technology Ltd Model predictive control of air pollution control processes
US7451004B2 (en) * 2005-09-30 2008-11-11 Fisher-Rosemount Systems, Inc. On-line adaptive model predictive control in a process control system
US20080243289A1 (en) 2007-03-28 2008-10-02 Honeywell International, Inc. Model maintenance architecture for advanced process control
JP2009181392A (ja) * 2008-01-31 2009-08-13 Omron Corp モデル予測制御方法およびモデル予測制御装置
CN101925866B (zh) * 2008-01-31 2016-06-01 费希尔-罗斯蒙特系统公司 具有用来补偿模型失配的调节的鲁棒的自适应模型预测控制器

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1742241A (zh) * 2002-12-09 2006-03-01 搭篷技术公司 具有可变动态特性的过程的自适应控制的一种系统和方法

Also Published As

Publication number Publication date
EP2561411A1 (en) 2013-02-27
CN102985884A (zh) 2013-03-20
EP2561411B1 (en) 2020-11-04
ES2842967T3 (es) 2021-07-15
US20130041482A1 (en) 2013-02-14
WO2011132050A1 (en) 2011-10-27
US9581980B2 (en) 2017-02-28

Similar Documents

Publication Publication Date Title
CN102985884B (zh) 用于更新模型预测控制器中模型的方法和系统
US9405286B2 (en) Use of predictors in process control systems with wireless or intermittent process measurements
US9436174B2 (en) Kalman filters in process control systems
US10061286B2 (en) Model predictive control using wireless process signals
US8095227B2 (en) Apparatuses, systems, and methods utilizing adaptive control
US9298176B2 (en) Compensating for setpoint changes in a non-periodically updated controller
US8756039B2 (en) Rapid process model identification and generation
US20080243289A1 (en) Model maintenance architecture for advanced process control
JP2002049409A (ja) プロセス制御システムにおける適応推定モデル
US20080274240A1 (en) Adaptive controller and expert system food processing
EP2884354A1 (en) Model-based predictive controller with steady-state model adaptation
JP7517774B2 (ja) 非周期的に更新されるコントローラにおける速度に基づく制御、プロセスを制御する方法、プロセスコントローラ
EP3270240B1 (en) Adaptive control techniques for ph control or control of other industrial processes
CN109642685A (zh) 使用定位器模型诊断气动控制阀的方法与装置
US20160378101A1 (en) Methods, apparatuses and systems for controlling a valve based on a combination of a characteristic curve for the valve and a proportional, integral and derivative signal valve
CN104615161A (zh) 一种基于模型预测的双回路水箱液位控制方法
CN111610779B (zh) 一种数据驱动的非线性系统执行器故障因子的辨识方法
DK177915B1 (en) Process control method
JP7521712B1 (ja) 制御装置、制御方法、及びプログラム
US20240149544A1 (en) System, apparatus and method for programmable logic control of pressure in a tire curing press bladder
JP4379698B2 (ja) シート状製品を製造する方法およびシート状製品製造装置
Alvarez et al. Bioprocess modelling for learning model predictive control (L-MPC)
RU2242040C1 (ru) Адаптивная цифровая система управления нестационарными объектами со связанными параметрами

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20191128

Address after: Baden, Switzerland

Patentee after: ABB Switzerland Co.,Ltd.

Address before: Zurich

Patentee before: ABB RESEARCH Ltd.

CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20151125