一种适用于多机电力系统调频分析的融合建模方法及装置与流程

本发明涉及电网调度通信领域，特别涉及一种适用于多机电力系统调频分析的融合建模方法及装置。

背景技术：

电力系统的频率稳定是保证电力系统的安全稳定运行的必要条件。目前对区域电力系统的研究多用于潮流分析，此类模型对机组的动态特性进行了很多简化；对于调频的动态特性分析的数学模型多局限于一到两种发电机组，这与实际电力系统运行情况相差较远；还有一些数学模型由于未将暂态过程和动态过程区分，致使运算速度十分缓慢；在考虑风电机组这类新能源发电机组时，一般未将有功功率和无功功率解耦，不利于有效的分析。并且，现有数学模型不仅鲜有考虑气象参数的随机性对调频的影响，而且更是鲜有将机理和数据相结合的适用于多机电力系统调频分析的融合模型。

技术实现要素：

本发明的主要目的是提出一种适用于多机电力系统调频分析的融合建模方法，旨在克服以上问题。

为实现上述目的，本发明提出的一种适用于多机电力系统调频分析的融合建模方法，包括：

s10向空冷机组模型输入空冷给定负荷信号、迎面风速信号、环境温度信号和上一次调频的系统频率的偏差量δω，向火电机组模型输入火电给定负荷信号和上一次调频的系统频率的偏差量δω，向水电机组模型输入水电给定负荷信号和上一次调频的系统频率的偏差量δω，向风电机组模型输入风速信号；

s20空冷机组模型输出的有功功率火电机组模型输出的有功功率水电机组模型输出的有功功率和风电机组模型输出的有功功率分别乘以各自的功率份额系数α1、α2、α3和α4之后叠加，获得电力系统中机组发电量的标幺值α1+α2+α3+α4＝1；

s30利用电力系统中机组发电量的标幺值和电网中的用户用电量的标幺值根据转子方程求得系统频率的偏差量δω和实际的系统频率ωm，其中ta∑表示等效转子时间常数，β∑表示等效摩擦系数，*表示取标幺值，表示同步发电机的转子转速的标幺值，也即系统频率的标幺值；

s40将所述系统频率的偏差量δω，用反馈的方式作为下一次调频的系统频率的偏差量δω；

s50采用实际数据对以上各机组模型进行验证及修正，生成适用于多机电力系统高频分析的模型。

优选地，以上所述各机组模型除风电机组模型外的其他机组模型还根据频率偏差参与调频，该调频方法为：设定转子方程中的积分初值为1，则转子方程输出为频率ωm，频率ωm减去频率参考值ωpref，将频率的偏差量用反馈的方式作为机组的输入，即可实现一次调频。

优选地，所述水电机组模型采用的是基于电力系统建模与仿真软件的区域电力系统模型中自带的水轮机模块，以给定负荷和频率偏差量为输入、输出为水电机组发出的有功功率；

所述空冷机组模型采用的是基于电力系统建模与仿真软件的区域电力系统模型中自带的空冷机组模块，以环境温度、迎面风速、给定负荷和频率偏差量为输入，输出空冷机组发出的有功功率；

所述火电机组模型采用的是基于电力系统建模与仿真软件的区域电力系统模型中自带的火电机组模块，以给定负荷和频率偏差量为输入，输出火电机组发出的有功功率；

所述风电机组模型均是采用的是基于电力系统建模与仿真软件的区域电力系统模型中自带的风电机组模块，以风速为输入，输出风电机组的有功功率。

优选地，所述火电机组发出的功率是由高压缸、两个高压加热器、中压缸和低压缸一起提供时，功率份额系数表示各部分发的功率占总功率的比例,为火电机组输出功率的标幺值，所述火电机组模型的建模方法为：

上一次调频的系统频率的偏差量δω除以调差系数δ1后为火电给定负荷信号减去所述后同时作为高压调节阀的传递函数和中压调节阀的传递函数的输入，高压调节阀的传递函数的输出与高压主汽阀的通过系数相乘后作为高压缸容积方程的输入，高压缸容积方程的输出作为第一高压加热器的输入；

中压调节阀的传递函数的输出、中压主汽阀的通过系数和再热容积方程的输出同时相乘后作为中压缸容积方程的输入，

所述中压缸容积方程的输入同时以负反馈的形式与第一高压加热器的输出作差，所述差作为再热容积方程的输入，

中压缸容积方程的输出利用第二高压加热器加热后，再利用除氧器进行加热，除氧器的输出作为低压缸容积方程的输入；

高压缸容积方程的输出、第一高压加热器的输出、中压缸容积方程的输出、第二高压加热器的输出和低压缸容积方程的输出分别乘以各自的功率份额系数c1、c2、c3、c4和c5后进行相加，相加后得到的信号为火电机组模型输出的有功功率c1+c2+c3+c4+c5＝1。

优选地，所述高压调节阀和中压调节阀还设有限幅环节，限定输出量的范围。

优选地，所述空冷机组模型的建模方法为：

上一次调频的系统频率的偏差量δω除以调差系数δ2后为空冷给定负荷信号减去所述后同时作为高压调节阀的传递函数和中压调节阀的传递函数的输入，高压调节阀的传递函数的输出与高压主汽阀的通过系数相乘后作为高压缸容积方程的输入，高压缸容积方程的输出作为第一高压加热器的输入；

中压调节阀的传递函数的输出、中压主汽阀的通过系数和再热容积方程的输出同时相乘后作为中压缸容积方程的输入，

所述中压缸容积方程的输入同时以负反馈的形式与第一高压加热器的输出作差，所述差作为再热容积方程的输入，

中压缸容积方程的输出利用第二高压加热器加热后，再利用除氧器进行加热，除氧器的输出作为低压缸容积方程的输入；

低压缸容积方程的输出、迎面风速信号和环境温度信号同时作为空冷系统方程的输入，空冷系统方程输出低压缸出口焓值，低压缸出口焓乘以低压缸输出量，即可得到低压缸输出的有功功率，所述空冷系统方程是运用matlab内的系统函数求出低压缸出口焓值；

高压缸容积方程的输出、第一高压加热器的输出、中压缸容积方程的输出、第二高压加热器的输出、低压缸发出的功率分别乘以各自的功率份额系数b1、b2、b3、b4和b5后进行相加，相加后得到的信号为空冷机组模型输出的有功功率b1+b2+b3+b4+b5＝1。

优选地，所述第一高压加热器、第二高压加热器和除氧器均采用回热抽汽模型来建模模拟，所述回热抽汽模型的建模方法包括：

回热抽汽模型的输入乘以1+α为d1，所述d1与d4相减为回热抽汽模型的输出，

回热抽汽模型输出的信号乘以后为d3，d3和d6相减后为d2，d2作为符号函数sign的输入，d2还作为函数的输入，抽汽逆止门的输入为常数1，抽汽逆止门的输出、符号函数sign的输出和函数的输出同时相乘后为d5，d5作为抽气惯性模型的输入，抽气惯性模型的输出乘以后为d6，d5乘以α后为d4；

所述抽汽逆止门的传递函数为：

tchf为抽汽逆止阀伺服阀时间常数，tcqf为抽汽逆止阀油动机时间常数，其中的积分器初值为1，sign表示符号函数，当输入为正数，结果为1，当输入为负数，结果为-1，当输入为0，结果为0，

函数表示对输入信号取绝对值，再开平方，

抽汽惯性模型为惯性环节，

第一高压加热器、第二高压加热器和除氧器中的传递函数分别为：

以上三个式子中，tgj1、tgj2、tcy分别为第一高压加热器、第二高压加热器、除氧器的回热抽汽惯性时间常数，比例系数α和β，在第一高压加热器中，分别等于αgj1和βgj1，在第二高压加热器中，分别等于αgj2和βgj2，在除氧器中，分别等于αcy和βcy。

优选地，所述高压主汽阀的通过系数和中压主汽阀的通过系数均为常数1；所述高压调节阀的传递函数和中压调节阀的传递函数均为：

thf为滑阀的惯性时间常数，tvh为油动机的惯性时间常数；

所述高压缸容积方程为

所述中压缸容积方程为

所述低压缸容积方程为

所述再热容积方程为

上式中，th为高压缸的惯性时间常数，tm为中压缸的惯性时间常数，tl为低压缸的惯性时间常数，trh再热环节的惯性时间常数，所有积分器初始值为1。

优选地，所述风电机组模型的建模方法为：

风速信号、风电机组的桨距角控制系统输出的桨距角和传动轴输出的角速度wwt同时作为风轮机的输入，其中所述风电机组指的是双馈异步风机，所述风轮机以风速、传动轴的角速度和桨距角为输入量，输出风轮机的机械转矩，所述传动轴以风轮机的机械转矩和风机转子转速为输入，输出传动轴的角速度和风机的机械转矩所述桨距角控制系统以风机转子转速和测量的风机输出功率为输入，输出桨距角；

根据风电机组的传动轴输出的机械转矩与函数f(u)输出的电磁转矩利用公式求得风轮机转子转速所述f(u)为根据风轮机转子转速计算电磁转矩的函数；

所述风轮机转子转速作为函数f(u)的输入；

所述风轮机转子转速和风轮机的输出作为传动轴的输入；

所述风轮机转子转速与传动轴的输出的机械转矩相乘后为风电机组模型输出的有功功率

所述风电机组模型输出的有功功率和所述风轮机转子转速同时作为桨距角控制系统的输入。

优选地，所述s30中转子方程求得系统频率的偏差量δω和实际的系统频率ωm的方法为：

积分器是指将两个输入之差加以累计，累积方法可以是从零递增，积分器的初始值为0。累积方法也可以是从某一设定值递减；由于模型中物理参量都取标幺值，因此从某一设定值递减时，积分器的初始值为1。

若积分器的初始值为0，则转子方程的输出为频率的偏差量δω，所述频率的偏差量δω加上频率参考值ωpref得到实际的系统频率ωm；

若积分器的初始值为1，则转子方程的输出为实际的系统频率ωm，所述实际的系统频率ωm减去频率参考值ωpref得到频率的偏差量δω。

用反馈的方式作为下一次调频的系统频率的偏差量δω，既可得到适用于多机电力系统调频分析的物理模型。

进一步，利用实际数据对模型进行校验和修正，得到最终的适用于多机电力系统调频分析的融合模型。

本发明的优点在于，本发明的模型着眼于区域电力系统一次调频的动态过程，不考虑无功和电压调节的暂态过程。与现有的考虑暂态过程的区域电力系统模型相比，仿真速度提高了几十倍。由于受到仿真速度的限制，传统模型一般是只考虑一种机组的单机无穷大系统。而使用本发明建模方法建立的模型的区域电力系统，全面包含风电、水电、火电和空冷机组模型，仿真效果更贴近实际。此外，风电机组和空冷机组的运行都受气象因素的影响，这在本发明的模型中加以考虑，对于研究气象因素影响调频过程有参考价值。另外，本发明利用实际数据对最终的物理模型进行校核，实现了将机理和数据相融合的建模策略，进一步提高了模型精度和实际应用可靠性。使用本发明的建模方法建立的模型研究区域电力系统的频率特性，给电网调度部门提供参考，保证电力系统的安全稳定运行。

附图说明

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

图1为本发明所述的适用于多机电力系统调频分析的模型的融合建模方法实现步骤；

图2为本发明所述的适用于多机电力系统调频分析的机理模型示意图；

图3为所述的火电机组模型的原理示意图；

图4为所述的空冷机组模型的原理示意图；

图5为所述的回热抽汽模型的原理示意图；

图6为所述的风电机组机组模型的原理示意图；

图7为具体实施方式七所述的将机理及数据的相融合的建模方法示意图，

本发明目的的实现、功能特点及优点将结合实施例，参照附图做进一步说明。

具体实施方式

下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明的一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

需要说明，若本发明实施例中有涉及方向性指示(诸如上、下、左、右、前、后……)，则该方向性指示仅用于解释在某一特定姿态(如附图所示)下各部件之间的相对位置关系、运动情况等，如果该特定姿态发生改变时，则该方向性指示也相应地随之改变。

另外，若本发明实施例中有涉及“第一”、“第二”等的描述，则该“第一”、“第二”等的描述仅用于描述目的，而不能理解为指示或暗示其相对重要性或者隐含指明所指示的技术特征的数量。由此，限定有“第一”、“第二”的特征可以明示或者隐含地包括至少一个该特征。另外，各个实施例之间的技术方案可以相互结合，但是必须是以本领域普通技术人员能够实现为基础，当技术方案的结合出现相互矛盾或无法实现时应当认为这种技术方案的结合不存在，也不在本发明要求的保护范围之内。

具体实施方式一：结合图1和图2说明本实施方式，本实施方式所述的适用于多机电力系统调频分析的机理与数据相融合的建模方法，它包括如下步骤：

步骤一：向空冷机组模型输入空冷给定负荷信号、迎面风速信号、环境温度信号和上一次调频的系统频率的偏差量δω，向火电机组模型输入火电给定负荷信号和上一次调频的系统频率的偏差量δω，向水电机组模型输入水电给定负荷信号和上一次调频的系统频率的偏差量δω，向风电机组模型输入风速信号；

步骤二：空冷机组模型输出的有功功率火电机组模型输出的有功功率水电机组模型输出的有功功率和风电机组模型输出的有功功率分别乘以各自的功率份额系数α1、α2、α3和α4之后叠加，获得电力系统中机组发电量的标幺值α1+α2+α3+α4＝1；

步骤三：利用电力系统中机组发电量的标幺值和电网中的用户用电量的标幺值根据转子方程求得系统频率的偏差量δω和实际的系统频率ωm，其中所述转子方程中ta∑表示等效转子时间常数，β∑表示等效摩擦系数，*表示取标幺值；

步骤四：将所述系统频率的偏差量δω，用反馈的方式作为下一次调频的系统频率的偏差量δω，适用于多机电力系统调频分析的模型建立完成。

步骤五：为了进一步提升整体模型的准确性及精度，使用实际数据对各个模型进行验证及修正，适用于多机电力系统调频分析的模型建立完成。

各机组模型均是为了计算各机组向系统提供的有功功率，同时除风电机组模型外的其他机组还根据频率偏差参与调频。

本模型中的水电机组采用的是matlab/simulink(电力系统建模与仿真软件)自带的水轮机模块。运用matlab/simulink软件，建立了包含风电、水电、火电和空冷机组的区域电力系统模型。模型总框图如图2所示。空冷机组以环境温度、迎面风速、给定负荷和频率偏差量为输入，输出为空冷机组发出的有功功率。火电机组以给定负荷和频率偏差量为输入，输出为火电机组发出的有功功率。水电机组以给定负荷和频率偏差量为输入，输出为水电机组发出的有功功率。风电机组以风速为输入，输出为风电机组的有功功率。将空冷机组、火电机组、水电机组和风电机组发出的有功功率(标幺值)分别表示为再用这些机组的有功功率分别乘以各自的功率份额系数α1、α2、α3、α4，再作和，即可得到区域电力系统所有机组提供的总的有功功率。功率份额系数可用机组的发电容量除以系统总容量得到，因此四个功率份额系数的和应为100％。根据转子方程(1)，即可计算出系统的频率。

式中，表示同步发电机的转子转速的标幺值，也即系统频率的标幺值，ta∑表示等效转子时间常数，本系统中除风机外均为同步发电机，将所有同步发电机等效为一台同步发电机，再用功率份额系数乘以各自的转子时间常数，即可求出等效转子时间常数ta∑；

β∑表示等效摩擦系数，与ta∑的求取方法是一样的，也是根据功率份额系数来求取；

分别为系统中机组发电量和用户用电量的标幺值；

由上式可知，转子方程是惯性环节，起延迟的作用，转子方程是将功率的不平衡经过惯性环节转换为频率的偏差量。将转子方程中的积分初值取0，则转子方程的输出为频率的偏差量δω，频率的偏差量加上频率的参考值，就可以得到实际的频率值ωm。若转子方程中的积分初值为1，则转子方程输出为频率ωm，频率ωm减去频率参考值ωpref，将频率的偏差量用反馈的方式作为机组的输入，即可实现一次调频。

具体实施方式二：结合图3说明本实施方式，本实施方式是对具体实施方式一所述的适用于多机电力系统调频分析的模型的机理建模方法的进一步限定，所述火电机组模型的建模方法为：

上一次调频的系统频率的偏差量δω除以调差系数δ1后为火电给定负荷信号减去所述后同时作为高压调节阀的传递函数和中压调节阀的传递函数的输入，高压调节阀的传递函数的输出与高压主汽阀的通过系数相乘后作为高压缸容积方程的输入，高压缸容积方程的输出作为第一高压加热器的输入；

中压调节阀的传递函数的输出、中压主汽阀的通过系数和再热容积方程的输出同时相乘后作为中压缸容积方程的输入，

所述中压缸容积方程的输入同时以负反馈的形式与第一高压加热器的输出作差，所述差作为再热容积方程的输入，

中压缸容积方程的输出利用第二高压加热器加热后，再利用除氧器进行加热，除氧器的输出作为低压缸容积方程的输入；

高压缸容积方程的输出、第一高压加热器的输出、中压缸容积方程的输出、第二高压加热器的输出和低压缸容积方程的输出分别乘以各自的功率份额系数c1、c2、c3、c4和c5后进行相加，相加后得到的信号为火电机组模型输出的有功功率c1+c2+c3+c4+c5＝1。

火电机组发出的功率是由高压缸、两个高压加热器、中压缸和低压缸一起提供，功率份额系数表示各部分发的功率占总功率的比例,为火电机组输出功率的标幺值，此外高压调节阀和中压调节阀还要加上限幅环节，限定输出量的范围。

具体实施方式三：结合图4说明本实施方式，本实施方式是对具体实施方式一所述的适用于多机电力系统调频分析的模型的机理建模方法的进一步限定，所述空冷机组模型的建模方法为：

上一次调频的系统频率的偏差量δω除以调差系数δ2后为空冷给定负荷信号减去所述后同时作为高压调节阀的传递函数和中压调节阀的传递函数的输入，高压调节阀的传递函数的输出与高压主汽阀的通过系数相乘后作为高压缸容积方程的输入，高压缸容积方程的输出作为第一高压加热器的输入；

中压调节阀的传递函数的输出、中压主汽阀的通过系数和再热容积方程的输出同时相乘后作为中压缸容积方程的输入，

所述中压缸容积方程的输入同时以负反馈的形式与第一高压加热器的输出作差，所述差作为再热容积方程的输入，

中压缸容积方程的输出利用第二高压加热器加热后，再利用除氧器进行加热，除氧器的输出作为低压缸容积方程的输入；

低压缸容积方程的输出、迎面风速信号和环境温度信号同时作为空冷系统方程的输入，空冷系统方程输出低压缸出口焓值，低压缸出口焓乘以低压缸输出量，即可得到低压缸输出的有功功率；

高压缸容积方程的输出、第一高压加热器的输出、中压缸容积方程的输出、第二高压加热器的输出、低压缸发出的功率分别乘以各自的功率份额系数b1、b2、b3、b4和b5后进行相加，相加后得到的信号为空冷机组模型输出的有功功率b1+b2+b3+b4+b5＝1。

空冷机组与火电机组相比，区别仅在低压缸出口处。所以空冷机组模型与火电机组模型相比多了一个空冷系统方程。空冷系统方程以低压缸输出量、环境温度、迎面风速为输入，输出低压缸出口焓值，低压缸出口焓乘以低压缸输出量，即可得到低压缸输出的有功功率。空冷系统方程内部是运用s函数求出低压缸出口焓值。

具体实施方式四：结合图5说明本实施方式，本实施方式是对具体实施方式二或三所述的适用于多机电力系统调频分析的模型的机理建模方法的进一步限定，

所述第一高压加热器、第二高压加热器和除氧器均采用回热抽汽模型来建模模拟，所述回热抽汽模型的建模方法包括如下步骤：

回热抽汽模型的输入乘以1+α为d1，所述d1与d4相减为回热抽汽模型的输出，

回热抽汽模型输出的信号乘以后为d3，d3和d6相减后为d2，d2作为符号函数sign的输入，d2还作为函数的输入，抽汽逆止门的输入为常数1，抽汽逆止门的输出、符号函数sign的输出和函数的输出同时相乘后为d5，d5作为抽气惯性模型的输入，抽气惯性模型的输出乘以后为d6，d5乘以α后为d4；

所述抽汽逆止门的传递函数为

tchf为抽汽逆止阀伺服阀时间常数，tcqf为抽汽逆止阀油动机时间常数，其中的积分器初值为1，sign表示符号函数，当输入为正数，结果为1，当输入为负数，结果为-1，当输入为0，结果为0。

函数表示对输入信号取绝对值，再开平方。

抽汽惯性模型为惯性环节，

第一高压加热器、第二高压加热器和除氧器中的传递函数分别为：

上面三个式子中，tgj1、tgj2、tcy分别为第一高压加热器、第二高压加热器、除氧器的回热抽汽惯性时间常数。

比例系数α和β，在第一高压加热器中，分别等于αgj1和βgj1，在第二高压加热器中，分别等于αgj2和βgj2，在除氧器中，分别等于αcy和βcy。

具体实施方式五：本实施方式是对具体实施方式二或三所述的适用于多机电力系统调频分析的模型的机理建模方法的进一步限定，

所述高压主汽阀的通过系数和中压主汽阀的通过系数均为常数1；所述高压调节阀的传递函数和中压调节阀的传递函数均为：

thf为滑阀的惯性时间常数，tvh为油动机的惯性时间常数；

所述高压缸容积方程为

所述中压缸容积方程为

所述低压缸容积方程为

所述再热容积方程为

式中，th为高压缸的惯性时间常数，tm为中压缸的惯性时间常数，tl为低压缸的惯性时间常数，trh再热环节的惯性时间常数，所有积分器初始值为1。

具体实施方式六：结合图6说明本实施方式，本实施方式是对具体实施方式一所述的适用于多机电力系统调频分析的模型的机理建模方法的进一步限定，所述风电机组模型的建模方法为：

风速信号、风电机组的桨距角控制系统输出的桨距角和传动轴输出的角速度wwt同时作为风轮机的输入；

根据风电机组的传动轴输出的机械转矩与函数f(u)输出的电磁转矩利用公式求得风轮机转子转速所述f(u)为根据风轮机转子转速计算电磁转矩的函数；

所述风轮机转子转速作为函数f(u)的输入；

所述风轮机转子转速和风轮机的输出作为传动轴的输入；

所述风轮机转子转速与传动轴的输出的机械转矩相乘后为风电机组模型输出的有功功率

所述风电机组模型输出的有功功率和所述风轮机转子转速同时作为桨距角控制系统的输入。

所述风电机组指的是双馈异步风机。图6中pitch是桨距角，是机械转矩，是电磁转矩，是根据风机转子转速计算得到。f(u)为根据风机转子转速计算电磁转矩的函数。图6中风轮机、传动轴、桨距角控制系统三个模块均采用传统方法建模，机械转矩乘以转子转速即可得到风机发出的功率

风轮机以风速、传动轴的角速度和桨距角为输入量，输出风轮机的机械转矩；传动轴以风轮机的机械转矩和风机转子转速为输入，输出传动轴的角速度和风机的机械转矩桨距角控制系统以风机转子转速和测量的风机输出功率为输入，输出桨距角pitch。

具体实施方式七：结合图7，本实施方式是对具体实施方式一所述的适用于多机电力系统调频分析的模型的机理建模方法的进一步限定，

所述步骤三中，根据转子方程求得系统频率的偏差量δω和实际的系统频率ωm的方法为：

若积分器的初始值为0，则转子方程的输出为频率的偏差量δω，所述频率的偏差量δω加上频率参考值ωpref得到实际的系统频率ωm；

若积分器的初始值为1，则转子方程的输出为实际的系统频率ωm，所述实际的系统频率ωm减去频率参考值ωpref得到频率的偏差量δω。

转子方程的求解本质就是一个积分器赋初值求解过程，积分器是指将两个输入之差加以累计，累积方法可以是从零递增，积分器的初始值为0。累积方法也可以是从某一设定值递减；由于模型中物理参量都取标幺值，因此从某一设定值递减时，积分器的初始值为1。

具体实施方式七：结合图7说明本实施方式，本实施方式是对具体实施方式一所述的适用于多机电力系统调频分析的模型的融合建模方法的进一步限定。如下图7所示，将实际数据输入上述机理模型，根据模型输出与实际数据的偏差对模型进行校正，直至模型输出与实际数据输出基本相同为止。

以上所述仅为本发明的优选实施例，并非因此限制本发明的专利范围，凡是在本发明的发明构思下，利用本发明说明书及附图内容所作的等效结构变换，或直接/间接运用在其他相关的技术领域均包括在本发明的专利保护范围内。