1.一种基于电能质量评估的牵引供电系统故障风险预警方法，其步骤为：a、数据采集与预处理

安装在牵引供电系统中的数据采集装置，采集得到以下数据矩阵：a1、牵引变压器一次侧的A相电压总谐波畸变率矩阵DA：

其中， 为第i天第j个采样时刻的A相电压谐波畸变率，i＝1，2，…，I；j＝1，2，…，J；I为采样的总天数，J为每一天的采样总次数；

a2、牵引变压器一次侧B相的电压总谐波畸变率矩阵DB：

其中， 表示第i天第j个采样时刻的B相电压谐波畸变率；

a3、牵引变压器一次侧C相的电压总谐波畸变率矩阵DC：

其中， 表示第i天第j个采样时刻的C相电压谐波畸变率；

a4、牵引变压器一次侧电压不平衡度矩阵Ε：

其中， 表示第i天第j个采样时刻的电压不平衡度值；

a5、牵引变压器二次侧母线电压有效值矩阵V：

其中， 表示第i天第j个采样时刻的牵引侧母线电压有效值；

b、计算电能质量数据的概率分布

b1、计算谐波畸变率数据的概率分布b1.1、数据处理装置算出A相电压总谐波畸变率矩阵DA中第j列的平均值，也即I天中的第j个采样时刻的A相电压总谐波畸变率的平均值 同时算出A相电压总谐波畸变率矩阵DA中第j列的方差，也即I天中的第j个采样时刻的A相电压总谐波畸变率的方差进而得到第j个采样时刻的A相电压总谐波畸变率 的概率密度函数其中exp表示指数函数运算；

b1.2、数据处理装置同时算出B相电压总谐波畸变率矩阵DB中第j列的平均值，也即I天中的第j个采样时刻的B相电压总谐波畸变率的平均值 并算出B相电压总谐波畸变率矩阵DB中第j列的方差，也即I天中的第j个采样时刻的B相电压总谐波畸变率的方差进而得到第j个采样时刻的B相电压总谐波畸变率 的概率密度函数b1.3、算出C相电压总谐波畸变率矩阵DC中第j列的平均值，也即I天中的第j个采样时刻的C相电压总谐波畸变率的平均值 同时算出C相电压总谐波畸变率矩阵DC中第j列的方差，也即I天中的第j个采样时刻的C相电压总谐波畸变率的方差进而得到第j个采样时刻的C相电压总谐波畸变率 的概率密度函数b2、计算电压不平衡度数据的概率分布数据处理装置算出电压不平衡度矩阵E中第j列的平均值，也即I天中的第j个采样时刻的电压不平衡度平均值 同时算出电压不平衡度矩阵E中第j列的方差，也即I天中的第j个采样时刻的电压不平衡度的方差进而得到第j个采样时刻电压不平衡度ej的概率密度函数f( ej )，b3、计算电压偏差数据的概率分布

数据处理装置算出电压偏差矩阵V中第j列的平均值，也即I天中的第j个采样时刻的电压偏差数据平均值 同时算出电压偏差矩阵V中第j列的方差，也即I天中的第j个采样时刻的电压偏差数据的方差进而得到第j个采样时刻电压偏差数据vj的概率密度函数f( vj )，c、计算电能质量因素的风险值

c1、计算谐波风险值

c1.1由下式得到第j个采样时刻A相电压谐波畸变率 的严重度函数

进而得到j个采样时刻A相电压谐波的风险值 其中表示对 进行微分运算， 表示从负无穷大到正无穷大的无穷积分；

c1.2由下式得到第j个采样时刻B相电压谐波畸变率 的严重度函数

进而得到j个采样时刻B相电压谐波的风险值 其中表示对 进行微分运算；

c1.3由下式得到第j个采样时刻C相电压谐波畸变率 的严重度函数

进而得到第j个采样时刻的C相电压谐波风险值 其中表示对 进行微分运算；

c2、计算电压不平衡风险值

由下式得到第j个采样时刻电压不平衡风险度ej严重度函数SE(ej)，

进而将电压不平衡风险严重度函数SE(ej)与电压不平衡度ej的概率密度函数f(ej)相乘，再计算其无穷积分，即得到j个采样时刻电压不平衡风险值 其中dej表示对ej进行微分运算；

c3、计算电压偏差风险值

由下式得到第j个采样时刻电压偏差vj严重度函数SV(vj)，

进而将电压偏差风险严重度函数SV(vj)与电压偏差vj的概率密度函数f(vj)相乘，再计算其无穷积分，即得到j个采样时刻电压偏差风险值 其中dvj表示对vj进行微分运算；

d、故障风险预警

d1、如第j个采样时刻A相电压谐波的风险值 小于1，则判定系统A相不会出现谐波故障，不给出谐波风险预警信息；否则，判定系统A相出现谐波故障的可能性大，给出系统在j时刻可能出现A相谐波故障的风险预警信息；

如第j个采样时刻B相电压谐波的风险值 小于1，则判定系统B相不会出现谐波故障，不给出谐波风险预警信息；否则，判定系统B相出现谐波故障的可能性大，给出系统在j时刻可能出现B相谐波故障的风险预警信息；

如第j个采样时刻C相电压谐波的风险值 小于1，则判定系统C相不会出现谐波故障，不给出谐波风险预警信息；否则，判定系统C相出现谐波故障的可能性大，给出系统在j时刻可能出现C相谐波故障的风险预警信息；

d2、电压不平衡风险预警

如第j个采样时刻电压不平衡风险值 小于1，则判定系统不会出现电压不平衡故障，不给出电压不平衡风险预警信息；否则，判定系统出现电压不平衡的可能性大，给出系统在j时刻可能出现电压不平衡故障的风险预警信息；

d3、电压偏差风险预警

如第j个采样时刻电压偏差风险值 小于1，则判定系统不会出现电压偏差故障，不给出电压偏差风险预警信息；否则，判定系统出现电压偏差的可能性大，给出系统在j时刻可能出现电压偏差故障的风险预警信息。