一种机械密封智能调节系统，包括：监测装置，实时测量与密封相关的信号；控制装置，包括状态判断模块、可行域分析模块、评估模块和寻优模块，其中，所述状态判断模块解读由监测装置得到的信号以估计密封的状态，所述可行域分析模块依据当前状态判断调节装置的调节作用可以使密封达到的状态，所述评估模块对密封的状态评估价值，所述寻优模块从密封调节的可行域中尽可能地寻找得到较高价值的调节方式，据此对调节装置发出指令；调节装置，根据控制装置给出的指令主动地对密封施加作用，本发明在实时监测并分析密封状态的基础上，依据用户要求的综合性能目标，寻找最佳的且通过调节可达的状态，并通过调节装置将密封调节到这一状态。
控制装置，包括状态判断模块、可行域分析模块、评估模块和寻优模块，其中，所述状态判断模块解读由监测装置得到的信号以估计密封的状态，所述可行域分析模块依据当前状态判断调节装置的调节作用可以使密封达到的状态，所述评估模块对密封的状态评估价值，所述寻优模块从密封调节的可行域中尽可能地寻找得到具有最高的安全和/或经济效益的调节方式，据此对调节装置发出指令；
调节装置，根据控制装置给出的指令主动地对密封施加作用；
所述密封采用静环浮动式设计，静环(1)经由弹簧(4)浮动支承在静环座(3)上，且与静环座(3)间由副密封(5)阻隔泄漏，可在施加电压时改变密封端面锥高δ的压电叠堆环(6)镶嵌在密封静环(1)上，锥高增大会增强密封的流体静压效应，增大静环与动环的间隙，从而减轻或消除接触，但会增大泄漏量；
所述监测装置采用静环背部对称布置的两个微型声发射传感器(7)，测量来自密封摩擦副的声发射信号；
所述状态判断模块依据声发射传感器测得的信号判断密封状态，具体方式如下：
(1)根据预先测定的材料声发射谱线库，由时间长度远小于转轴周期的一小段声发射信号得出接触剧烈程度和接触压力分布，通过实验测定，标定由声发射信号各频带的功率到两环轴向中心距离和两环夹角的函数；
(2)对比两个微型声发射传感器(7)测得的信号的相位得到声源中主要成分的大致位置，由这一位置随转轴旋转的变化，及两环轴向中心距离和两环夹角随转轴旋转的变化，得出动环偏摆、静环偏斜力矩和轴向压紧程度；
所述可行域分析模块，在调节装置可达的范围内预设若干个调节装置锥高节点，分析各节点处密封将会处于的状态，并通过它们插值得到可达状态的曲线；
所述评估模块，首先由给定的密封状态得到预期的磨损率和泄漏率，然后将它们折算成密封环磨损和物料泄漏带来的经济损失率；
所述寻优模块，在由前述的可行域分析模块得出的可行域中进行一维寻优，以得到可以使经济损失率尽可能小的目标状态，并分析达成目标状态所需的锥高，调节装置据此改变压电叠堆环(6)两端的电压；
或者，密封采用静环浮动式设计，静环(1)端面上的三个节流孔(8)分别被导到独立受控的三个流体源(10)上，每个流体源的压强都可在一定范围内独立地改变，调节这些流体源的压力能够影响端面上的流体压力分布，从而实现对密封状态的调整，其能够有效地对两环整体的靠近或远离以及静环的倾斜进行补偿，但对于动环倾斜、密封端面的磨损、密封支承零部件老化、流体介质含磨粒，则无法直接补偿；
所述监测装置采用四个测量密封环运动的电涡流位移传感器(9)，其中三个周向均布在静环的背面，一个布置在动环的背面；
所述状态判断模块中，由已测得的密封环运动，判断密封的动环偏摆状况和静环浮动支承状况：动环背面传感器测得的即为动环的偏摆运动；静环背面所测得的运动中，滤除与动环背面传感器测得的一致的成分，即为静环浮动支承的状况，最终用四个参量表示密封状态；
所述可行域分析模块，预设三个流体源压强极限设定值的2 3 ＝8种组合，采用机械密封数值模型分析每种组合下密封将会处于的状态，并由它们围成调节可达状态的大致边界；
所述评估模块，由给定的密封状态得出预期的磨损率，在考虑此磨损率影响的情况下计算未来的预期泄漏率变化，并进一步推算这样的未来泄漏变化条件下冷却系统将工作机降到安全温度以下所需的时间；
所述寻优模块，采用模拟退火算法，在得出的可行域中随机地变换产生候选状态，利用评估模块计算其时间需求，并以一定的接受准则决定是否用新的候选状态进行下一步迭代，最终找到使时间需求尽可能小的目标状态，并在找到后分析达成目标状态所需的受控流体源压强，发送给调节装置。
2.根据权利要求1所述机械密封智能调节系统，其特征在于，所述监测装置主要包括传感器、信号预处理装置、采集卡、采集软件及必要的作为载体的终端设备。
3.根据权利要求1所述机械密封智能调节系统，其特征在于，所述与密封相关的信号包括上下游压力、膜内压力、转速、密封环温度、密封环位移、密封环间电容以及摩擦副声发射，所述调节装置对密封施加的作用为力、压强变化、形状变化以及加热/冷却中的一种或者多种组合。
4.根据权利要求1所述机械密封智能调节系统，其特征在于，所述状态判断模块对监测装置测得的信号进行滤波并按照标定关系得出被测物理量的值；所述评估模块首先由给定的密封状态计算其性能，然后将性能参数按照实际应用场景的需求综合起来，从而将每个密封状态映射为一个价值量。
5.一种基于权利要求1所述机械密封智能调节系统的机械密封智能调节方法，其特征在于，实时测量与密封相关的信号，执行控制算法，生成对密封施加作用的指令，进行密封调节，其中，所述控制算法依次包括状态判断环节、可行域分析环节、评估环节和寻优环节，所述状态判断环节解读由监测装置得到的信号以估计密封的状态，所述可行域分析环节依据当前状态判断调节装置的调节作用可以使密封达到的状态，所述评估环节对密封的状态评估价值，所述寻优环节从密封调节的可行域中尽可能地寻找得到较高价值的调节方式，据此发出指令。
6.根据权利要求5所述机械密封智能调节方法，其特征在于，所述与密封相关的信号包括上下游压力、膜内压力、转速、密封环温度、密封环位移、密封环间电容以及摩擦副声发射，所述对密封施加的作用为力、压强变化、形状变化以及加热/冷却中的一种或者多种组合。
7.根据权利要求5所述机械密封智能调节方法，其特征在于，所述状态判断环节中，对测得的信号进行滤波并按照标定关系得出被测物理量的值，当传感器的测量并不能提供确定密封状态的全部信息时，状态判断算法给出密封状态的概率分布，并依据调节后的信号变化作进一步判断。
8.根据权利要求5所述机械密封智能调节方法，其特征在于，所述评估环节中，首先由给定的密封状态计算其性能，然后将性能参数按照实际应用场景的需求综合起来，从而将每个密封状态映射为一个价值量。
9.根据权利要求5所述机械密封智能调节方法，其特征在于，所述状态判断、可行域分析和评估均依据密封的特性模型设计精确或近似计算的方法，所述密封的特性模型指特定状态的密封会带来怎样的性能和产生怎样的可监测要素，及未来状态将如何转变。
一种机械密封智能调节系统与方法
技术领域
[0001] 本发明属于流体密封领域和自动控制领域，特别涉及一种机械密封智能调节系统与方法。
背景技术
[0002] 机械密封是旋转机械设备中的一种常见轴端密封形式。其配对的动环和静环相对旋转，形成具有一定稳定性的摩擦副，在限制泄漏的同时减轻甚至消除接触。然而，这种来自于对密封的润滑、动力学、传热与变形的精妙设计的稳定性有时仍然无法在较强的异常因素的影响下维持摩擦副的稳定，导致密封性能下降，寿命缩短，甚至发生无预兆的失效。
[0003] 针对这一问题，人们提出引入在密封中引入外加的调节装置，在实时监测的基础上，采用反馈控制的方法尝试将密封维持在某种理想状态附近。这种方法的不足之处在于，由于机械密封机理的复杂性和结构的紧凑性，调节装置在机械密封系统中的调节能力往往是有限的，这包括两个方面：一是许多异常因素是无法(或至少很难)被一种或多种特定的调节装置的功能所刚好抵消的，更何况实践上不同种类的调节装置在占用空间上可能存在干涉；二是调节装置的调节强度是有限的，尤其是在一些高参数的机械密封中(通常来说，对可靠性有极高要求的也正是高参数的密封)，密封本身的机理带来的稳定性也相当强大，调节装置所能带来的改变可能并不占据支配地位。这两点意味着，如果按反馈控制的方法设定一个确定的目标状态，这个目标状态一般是不可达的，然而作为非线性系统的密封的状态与性能的关系是复杂的，即：并非状态越接近某个高性能的目标状态，性能就越佳。进一步说，在不同的具体应用场合下，用户对密封的各方面性能(如泄漏、磨损量、扭矩等)可能有着不同的需要，应当从安全性、经济性等多个角度去综合评判对性能的需要。
发明内容
[0004] 为了克服上述现有技术的缺点，本发明的目的在于提供一种机械密封智能调节系统与方法，其在实时监测并分析密封状态的基础上，依据用户要求的综合性能目标，寻找最佳的且通过调节可达的状态，并通过调节装置将密封调节到这一状态。
[0005] 为了实现上述目的，本发明采用的技术方案是：
[0006] 一种机械密封智能调节系统，包括：
[0007] 监测装置，实时测量与密封相关的信号；
[0008] 控制装置，包括状态判断模块、可行域分析模块、评估模块和寻优模块，其中，所述状态判断模块解读由监测装置得到的信号以估计密封的状态，所述可行域分析模块依据当前状态判断调节装置的调节作用可以使密封达到的状态，所述评估模块对密封的状态评估价值，所述寻优模块从密封调节的可行域中尽可能地寻找得到具有最高的安全和/或经济效益的调节方式，据此对调节装置发出指令；
[0009] 调节装置，根据控制装置给出的指令主动地对密封施加作用。
[0010] 所述监测装置主要包括传感器、信号预处理装置、采集卡、采集软件及必要的作为载体的终端设备。
[0011] 所述与密封相关的信号包括上下游压力、膜内压力、转速、密封环温度、密封环位移、密封环间电容以及摩擦副声发射，所述调节装置对密封施加的作用为力、压强变化、形状变化以及加热/冷却中的一种或者多种组合。
[0012] 所述状态判断模块对监测装置测得的信号进行滤波并按照标定关系得出被测物理量的值；所述评估模块首先由给定的密封状态计算其性能，然后将性能参数按照实际应用场景的需求综合起来，从而将每个密封状态映射为一个价值量。
[0013] 本发明还提供了一种机械密封智能调节方法，实时测量与密封相关的信号，执行控制算法，生成对密封施加作用的指令，进行密封调节，其中，所述控制算法依次包括状态判断环节、可行域分析环节、评估环节和寻优环节，所述状态判断环节解读由监测装置得到的信号以估计密封的状态，所述可行域分析环节依据当前状态判断调节装置的调节作用可以使密封达到的状态，所述评估环节对密封的状态评估价值，所述寻优环节从密封调节的可行域中尽可能地寻找得到较高价值的调节方式，据此发出指令。
[0014] 所述状态判断环节中，对测得的信号进行滤波并按照标定关系得出被测物理量的值，当传感器的测量并不能提供确定密封状态的全部信息时，状态判断算法给出密封状态的概率分布，并依据调节后的信号变化作进一步判断。
[0015] 所述状态判断、可行域分析和评估均依据密封的特性模型设计精确或近似计算的方法，所述密封的特性模型指特定状态的密封会带来怎样的性能和产生怎样的可监测要素，及未来状态将如何转变。
[0016] 与现有技术相比，本发明的有益效果是：
[0017] 针对机械密封主动调节中线性反馈算法由于调节机构的能力限制和密封系统的非线性特征而不适用的问题，根据具体的应用场合制定最优化策略。
附图说明
[0018] 图1是本发明系统三个装置的典型组合方式，必要时组合方式也可进行调整。在它们的作用下，密封可以尽可能地达成具体应用场景下的最大价值，达到精准调控的目的。
[0019] 图2是本发明实施例1中采用静环浮动式设计的密封结构示意图。
[0020] 图3是本发明实施例1中锥高调节的示意图。
[0021] 图4是本发明实施例2中采用静环浮动式设计的密封结构示意图。
[0022] 图5是图3中的静环端面示意图。
具体实施方式
[0023] 为了使本发明的发明目的、技术方案及技术效果更加清楚明白，以下结合附图对本发明的具体实施例进行描述。应当理解，此处所描述的具体实施例仅用以解释本发明，并不用于限定本发明。
[0024] 参考图1，本发明为一种机械密封智能调节系统，其在实时监测并分析密封状态的基础上，依据用户要求的综合性能目标，寻找最佳的且通过调节可达的状态，并通过调节装置将密封调节到这一状态，其包括：
[0025] (1)监测装置
[0026] 监测装置主要包括传感器、信号预处理装置(如放大器、信号转换器等)、采集卡、采集软件，及必要的作为载体的终端设备(如工控机)。监测装置实时地测量与密封相关的多种信号，可用的通常包括上下游压力、膜内压力、转速、密封环温度、密封环位移、密封环间电容、摩擦副声发射等。
[0027] (2)控制装置
[0028] 以软件和运行它的终端设备为控制算法载体。区别于传统的反馈控制，本发明控制算法主要包括以下环节：状态判断、可行域分析、评估和寻优。
[0029] 状态判断环节解读由监测装置得到的信号以估计密封的状态。其对测得的信号进行滤波并按照标定关系得出被测物理量的值。对部分信号种类来说，还需要依据相应的密封模型进行推算(比如，由密封环中某点测得的温度，依据密封的传热分析来推算密封端面摩擦剧烈程度)。有时传感器的测量并不能提供确定密封状态的全部信息，此时状态判断算法给出密封状态的概率分布，并依据调节后的信号变化作进一步判断。
[0030] 可行域分析环节依据当前状态判断调节装置的调节作用可以使密封达到的状态。
[0031] 评估环节对密封的状态评估价值。其首先由给定的密封状态计算其性能，然后将性能参数按照实际应用场景的需求综合起来。总的说来这一环节将每个密封状态映射为一个价值量。
[0032] 寻优环节从密封调节的可行域中尽可能地寻找得到较高价值的调节方式，据此对调节装置发出指令。
[0033] 其中，状态判断、可行域分析和评估均依据密封的特性模型——即特定状态的密封会带来怎样的性能和产生怎样的可监测要素，及未来状态将如何转变——设计精确或近似计算的方法。
[0034] (3)调节装置
[0035] 调节装置根据控制算法给出的指令主动地对密封施加作用(这种作用可能是各种形式的，如力、压强变化、形状变化、加热/冷却等，当然也包括了多种形式作用的组合)。
[0036] 以下是本发明的两个具体实施例。
[0037] 实施例1：
[0038] 某机械密封用于密封无危险、无污染的物料，希望在避免密封环的磨损以避免更换密封环的成本的同时，尽可能减少泄漏以提高物料产率，最终实现经济效益的最大化。
[0039] 请参见图2，密封采用静环浮动式设计(静环1经由弹簧4浮动支承在静环座3上，且与静环座3间由副密封5阻隔泄漏)。将压电叠堆环6镶嵌在密封静环1上，其可在施加电压时改变密封端面的锥高δ(见图3)：锥高(本例仅涉及收敛锥)增大会增强密封的流体静压效应，增大静环与动环的间隙，从而减轻或消除接触，但会增大泄漏量。
[0040] 采用静环背部对称布置的两个微型声发射传感器7测量来自密封摩擦副的声发射信号。
[0041] 状态判断算法依据声发射传感器测得的信号判断密封状态。具体方式如下：
[0042] (1)根据预先测定的材料声发射谱线库(特定材料配副接触摩擦引发的声发射信号出现的若干频带，以及摩擦剧烈程度不同导致的各频带比例变化)，由一小段(时间长度远小于转轴周期的)声发射信号得出接触剧烈程度和接触压力分布(当接触更多地源于两环的压紧时，接触压力的分布会更为分散；若接触更多地源于两环的偏斜，则接触压力的分布更集中)。通过实验测定，可以标定由声发射信号各频带的功率到两环轴向中心距离和两环夹角的函数；
[0043] (2)对比两个声发射传感器测得的信号的相位来得到声源中主要成分的大致位置。由这一位置随转轴旋转的变化，及(1)得出的两环轴向中心距离及两环夹角随转轴旋转的变化，得出动环偏摆、静环偏斜力矩和轴向压紧程度。
[0044] 可行性分析算法在调节装置可达的范围内预设若干个调节装置锥高节点，分析各节点处密封将会处于的状态，并通过它们插值得到可达状态的曲线。采用已经比较成熟的机械密封数值模型即可实现这样的分析。
[0045] 评估算法首先由给定的密封状态得到预期的磨损率(与接触压力的分布和转速有关)和泄漏率(由动环偏摆、静环偏斜力矩和轴向压紧程度计算出)，然后将它们折算成密封环磨损和物料泄漏带来的经济损失率。
[0046] 寻优算法在由前述的可行性算法得出的可行域中进行一维寻优，以得到可以使经济损失率尽可能小的目标状态(由算法中生成的候选目标状态计算经济损失率的过程通过前述评估算法完成)，并分析达成目标状态所需的锥高，调节装置据此改变压电叠堆环6两端的电压。考虑到问题的非线性，为避免收敛到局部最优解，采用模拟退火算法进行寻优：
[0047] (1)针对可行域为一维这一特征，将可行域变换为一个状态参量p；
[0048] (2)取可行的初始状态参量p(0)，由前述评估算法计算其相应的损失率q(0)；
[0049] (3)令n＝0；
[0050] (4)基于当前状态参量，即p(n)，随机变换出p’(n+1)＝p(n)+e(q(n))，其中e是一个随机量，其服从以k1·q(n)为标准差、0为均值的高斯分布(但从分布中截去使p’(n+1)不可行的部分)；
[0051] (5)由前述评估算法计算p’(n+1)所对应的损失率q’(n+1)；
[0052] (6)若q’(n+1)≤q(n)则接受新解，即令p(n+1)＝p(n)，相应地q(n+1)＝q’(n+1)；否则，以P＝exp(-(q’(n+1)-q(n))/(k2·q(n)))的概率接受新解。若未接受新解，则令p(n+
1)＝p(n)，相应地q(n+1)＝q(n)；
[0053] (7)判断是否满足结束条件(n达到上限或连续多次未接受新解)，若结束则将此p(n+1)作为结果；否则，n累加1，跳转到(4)。
[0054] 实施例2：
[0055] 某用于冷却系统的密封，需要在发生故障时做出应急调节，尽可能保证工作机温度降到足够低之前冷却系统保持正常工作。
[0056] 请参见图4，密封采用静环浮动式设计，图5为其静环端面示意图，静环1端面上的三个节流孔8分别被导到独立受控的三个流体源10上，每个流体源的压强都可在一定范围内独立地改变。调节这些流体源的压力可以影响端面上的流体压力分布，从而实现对密封状态的调整。其可以有效地对两环整体的靠近或远离以及静环的倾斜进行补偿，但对于其他问题(如动环倾斜、密封端面的磨损、密封支承零部件老化、流体介质含磨粒等)则无法直接补偿，但依然可能通过合理的调节来间接地减少其不利影响。
[0057] 布置四个电涡流位移传感器9以测量密封环的运动。其中三个(9-1,2,3)周向均布在静环的背面，一个(9-4)布置在动环的背面。
[0058] 状态判断算法由已测得的密封环运动，判断密封的动环偏摆状况和静环浮动支承状况：动环背面传感器测得的即为动环的偏摆运动；静环背面所测得的运动中，滤除与动环背面传感器测得的一致的成分(设静环背面测得运动为Us，动环背面测得运动为Ur，定义一种经标准化的范数||·||并寻找min(||Us-α·Ur||)，给出范数最小化后的Us-α·Ur)，即为静环浮动支承的状况(其进一步分解为分为轴向和两个角向共三个分量)，最终用四个参量表示密封状态。
[0059] 可行性分析算法预设了三个流体源压强极限设定值的2^3＝8种组合，采用机械密封数值模型分析每种组合下密封将会处于的状态，并由它们围成调节可达状态的大致边界。
[0060] 评估算法由给定的密封状态得出预期的磨损率，在考虑此磨损率影响的情况下计算未来的预期泄漏率变化，并进一步推算这样的未来泄漏变化条件下冷却系统将工作机降到安全温度以下所需的时间。
[0061] 寻优算法采用模拟退火算法，在由前述的可行性算法得出的可行域中随机地变换产生候选状态，用前述评估算法计算其时间需求，并以一定的接受准则决定是否用新的候选状态进行下一步迭代。最终找到使时间需求尽可能小的目标状态，并在找到后分析达成目标状态所需的受控流体源压强，发送给调节装置。具体做法与实施例1的基本类似：
[0062] (1)将前述的表示密封状态的四个参量标准化，定义为状态向量p；
[0063] (2)取可行的初始状态向量p(0)，由前述评估算法计算其相应的时间需求q(0)；
[0064] (3)令n＝0；
[0065] (4)基于当前状态参量，即p(n)，随机变换出p’(n+1)＝p(n)+e(q(n))，其中e是一个随机向量，其服从以k1·q(n)为标准差、0为均值的四维高斯分布(但从分布中截去使p’(n+1)不可行的部分)；
[0066] (5)由前述评估算法计算p’(n+1)所对应的时间需求q’(n+1)；
[0067] (6)若q’(n+1)≤q(n)则接受新解，即令p(n+1)＝p(n)，相应地q(n+1)＝q’(n+1)；否则，以P＝exp(-(q’(n+1)-q(n))/(k2·q(n)))的概率接受新解。若未接受新解，则令p(n+
1)＝p(n)，相应地q(n+1)＝q(n)；
[0068] (7)判断是否满足结束条件(n达到上限或连续多次未接受新解)，若结束则将此p(n+1)作为结果；否则，n累加1，跳转到(4)。
