周学均，陈小强，谢 磊

（1.中电国际胡布发电有限公司,巴基斯坦 卡拉奇74200；2.中电华创电力技术研究有限公司，江苏 苏州215123；3.浙江大学 智能系统与控制研究所，浙江 杭州310027）

现代化工业过程规模庞大，有数百至数千个高度耦合的控制回路，为了使生产效率最大化，保证产品质量，所有关键回路必须处于或接近最佳状态下运行[1]。振荡是过程控制系统性能下降的常见原因[2]。Z震荡可以在一个控制回路中产生并传播到其他单元，形成厂级振荡[3]。这些震荡有可能造成能源和材料的浪费，产品质量波动，甚至威胁控制系统的安全性。因此，对工业过程回路进行准确的因果分析，快速查明故障根源并采取相应的纠正措施，是保证安全生产和提高经济效益的关键。然而，由于工业控制系统中采集的数据往往被非平稳因素、噪声、周期项等干扰，对振荡源进行定位是非常有挑战的任务，因此大量的因果关系分析方法被用于震荡源的定位研究。

文献[4]使用互相关分析法来进行因果分析。方法虽然简单，但不适用于非线性的情况。文献[5]将主成分分析与多变量格兰杰因果关系结合起来进行厂级振荡的定位。格兰杰因果方法是一个广泛应用的方法，但对变量独立性和线性具有较严苛的要求。除了这些时域的因果分析方法，谱格兰杰、偏有向相干[6]、有向传递函数[7]等一些频域的方法也被开发出来。然而，谱格兰杰是基于傅里叶变换的方法，因此只能适用于线性平稳过程。……