首页 / 专利库 / 电脑编程 / 算法 / 基于卡尔曼滤波器的输电线路故障定位方法

基于卡尔曼滤波器的输电线路故障定位方法

阅读:506发布:2023-05-28

专利汇可以提供基于卡尔曼滤波器的输电线路故障定位方法专利检索,专利查询,专利分析的服务。并且基于卡尔曼 滤波器 的输电线路故障 定位 方法,给定输电线模型和测量值;引入故障状态α;通过基于卡尔曼滤波器EnKF的动态状态估计DSE过程,得到输电线路故障类型和故障定位的初步估计值;将附加信息作为额外的约束,通过附加信息重新运行动态状态估计DS,以执行增强的卡尔曼滤波器EnKF;测量冗余的增加,得到对输电线路故障 位置 准确估计值。本 发明 方法是一种两步骤的集成卡尔曼滤波器(EnKF)方法,在不需要推导雅可比矩阵、降低了实现 算法 的难度的同时,充分利用了模型非线性和误差统计的全面演化,使故障位置预估更加精确,数据更加稳定。,下面是基于卡尔曼滤波器的输电线路故障定位方法专利的具体信息内容。

1.基于卡尔曼滤波器的输电线路故障定位方法,其特征在于包括以下步骤:
步骤1:给定输电线模型和测量值;
步骤2:引入故障状态α;
步骤3:通过基于卡尔曼滤波器EnKF的动态状态估计DSE过程,得到输电线路故障类型和故障定位的初步估计值;
步骤4:将附加信息作为额外的约束,通过附加信息重新运行动态状态估计DS,以执行增强的卡尔曼滤波器EnKF;
步骤5:测量冗余的增加,得到对输电线路故障位置准确估计值。
2.根据权利要求1所述基于卡尔曼滤波器的输电线路故障定位方法,其特征在于:
所述步骤1中,给先定输电线模型和测量值:
xk=f(xk-1,uk-1)+wk-1    (1);
zk=h(xk,uk)+vk-1      (2);
其中,公式(1)(2)中k是步长数;
公式(1)中,xk是k时刻输电线路的状态,xk-1是k-1时刻输电线路的状态,uk-1是k-1时刻系统的控制量,f(xk-1,uk-1)是前一时刻状态xk-1、前一时刻控制量uk-1和xk之间的映射关系,wk-1表示对状态量xk处理和测量过程伴随的随机噪声,其均值为零,方差为Qd;
公式(2)中,zk是k时刻的测量值,h(xk,uk)是xk,uk和zk之间的映射关系,vk-1表示对测量值zk处理和测量过程伴随的随机噪声,其均值为零,方差为Rd。
3.根据权利要求1所述基于卡尔曼滤波器的输电线路故障定位方法,其特征在于:
所述步骤2中,当输电线路发生故障时,从故障到终端的距离被视为故障期间不变的新状态α,输电线模型变成:
公式(3)中,αk是k时刻的故障状态,xk是k时刻的输电线路状态,f(xk-1,uk-1,ak-1)是xk-1,uk-1,αk-1和xk之间的映射关系,wa,k-1是k-1时刻的过程和测量产生的噪声;
当引入故障状态α作为状态时,输电线模型的非线性会增加,故定义一个扩展的状态向量为y=(x,α)T;将公式(1)和公式(2)重写为:
yk=f(yk-1,uk-1)+wk-1    (4);
zk=h(yk,uk)+vk-1    (5);
公式(4),公式(5)中k是步长数;
公式(4)中的yk是k时刻输电线路的状态,yk-1是k-1时刻输电线路的状态,f(yk-1,uk-1)是yk-1,uk-1和yk之间的映射关系,uk-1是k-1时刻对系统的控制量,wk-1表示时刻的过程和测量产生的噪声;
公式(5)中的zk是k时刻的测量值,h(yk,uk)是yk,uk和zk之间的映射关系,vk-1表示k-1过程和测量产生的噪声。
4.根据权利要求1所述基于卡尔曼滤波器的输电线路故障定位方法,其特征在于:
步骤3中,基于EnKF的DSE使用多次的预测校正过程进行分析,如下所示:
EnKF预测:
公式(6)中, 表示输电线路在k时刻的状态且由预测步骤计算得到,uk-1是k-1时刻对系统的控制量, 表示 uk-1和 之间的映射关系, 表示k-1时刻的过程和
观测产生的噪声且有校正步骤计算得到;
公式(7)中, 表示时刻的测量值且由预测步骤计算得到,uk是k时刻对系统的控制量,表示 uk和 之间的映射关系, 表示k-1时刻的过程和观测产生的噪声且有
校正步骤计算得到;
公式(8)中, 表示输电线路在k时刻的状态且由预测步骤计算得到,N是样本数, 是yk的共轭矩阵;
公式(9)中, 表示时刻的测量值且由预测步骤计算得到,N是样本数, 是zk的共轭矩阵;
EnKF校正:
公式(10)(11)(12)(13)中,N是样本数, 是误差协方差矩阵的共轭矩阵, 是雅克比矩阵的转置矩阵, 表示输电线路在k时刻的状态且由预测步骤计算得到, 表示时刻的测量值且由预测步骤计算得到,Kk是卡尔曼增益;
其中,Pk是误差协方差矩阵,Hk是雅克比矩阵,Kk是卡尔曼增益,上标i-表示在预测步骤计算状态和计算观测值。
5.根据权利要求4所述基于卡尔曼滤波器的输电线路故障定位方法,其特征在于:
基于EnKF的动态状态估计DSE过程,由权利要求4基于EnKF的DSE预测校正过程分析基础上,已知故障前和故障后的多段输电线路模型,在出现故障区段的连接点处稍有不同,无故障段连接点的三相电压电流总是相等的;
公式(14)中, 表示三相输电线上任意一点的左右两端随时间变化
的电压值;
公式(15)中, 表示三相输电线任意一点的左右两端随时间变化的
电流值;
公式(16)中, 表示三相电路故障点左右两端随时间变化的电压值;
将公式(14)-(16)包含在基于EnKF的DSE预测校正过程中,将得到故障定位 的一个近似值,以及故障点处的三相电压的估计值;根据故障点处的三相电压的估计值,能够对故障类型进行分类。
6.根据权利要求5所述基于卡尔曼滤波器的输电线路故障定位方法,其特征在于:
故障阻抗越小,出现故障相电压的变化越大;
故障类型是通过检查三相电压来确定的,若B和C相电压等于标称值,而A相电压很小,那么这是一个A-G单相接地短路故障;若C相电压等于标称值,而A相和B相的电压值接近标称值的一半,那么这是一个A-B的相间短路故障;若C相电压等于标称值,但是相A和B电压非常小,接近零,那么它是一个AB-G两相接地的短路故障。
7.根据权利要求1所述基于卡尔曼滤波器的输电线路故障定位方法,其特征在于:
步骤4中,利用步骤3所得的故障类型,导出了故障点处的附加电流关系式;
如果第一步确定故障是A-G故障,则附加的电流关系是:
公式(17)中, 分别表示故障线路的b相的任意一点左右两端随时间变化
的电流值,即B相的电流值为标称值;
公式(18)中, 分别表示故障线路的c相的任意一点左右两端随时间变化
的电流值,即C相的电流值为标称值;
当重新运行增强的基于EnKF的DSE时,上述第一步所得到的故障类型和故障点距离的估计值作为额外的约束,作为增加状态估计过程的测量冗余;
如果是AB-G两相接地的短路故障,电流关系为:
公式(19)中, 分别表示故障线路的c相的任意一点左右两端随时间变化
的电流值,即C相的电流为标称值,且A,B两相的电压值接近零;
如果是A-B相间短路故障,电流关系为:
公式(20)中, 分别表示故障线路的c相的任意一点左右两端随时间变化
的电流值,即C相的电流为标称值,且A,B两相的电压值约为标称值的一半;
如果是三相相间短路故障,三相的电压值均接近零。
8.根据权利要求1所述基于卡尔曼滤波器的输电线路故障定位方法,其特征在于:
步骤5中,在重新运行增强EnKF的DSE过程时,由于较好的初始状态和增加的测量冗余,实现了对故障位置的更精确的估计,执行增强的EnKF过程,能够给出更精确的故障定位估计值。

说明书全文

基于卡尔曼滤波器的输电线路故障定位方法

技术领域

[0001] 本发明涉及输电线路故障诊断技术领域,具体涉及一种基于卡尔曼滤波器的输电线路故障定位方法。

背景技术

[0002] 现有输电线路故障诊断技术中,准确地指出输电线路上的故障,能够最大限度地节省人和成本。当前输电线路故障测距技术主要分为三大类:单端阻抗法、双端阻抗法和行波法。
[0003] 单端阻抗法和双端阻抗法,均利用线路端子计算故障点到终端的阻抗从而确定故障位置,这两种方法性能受故障路径电阻、负载以及源参数的影响;行波法对波到达所需时间进行监测,以确定故障的位置,需要极高的采样率,且当故障发生点接近总线时无法保证结果精确度。
[0004] 此外,基于DSE(动态状态估计)的输电线路故障定位方法,该方法使用加权最小二乘法WLS,将故障距离视为DSE(动态状态估计)处理过程的状态,但是DSE(动态状态估计)需要准确的故障类型,并且雅可比矩阵结构复杂,会提升计算难度,同时将加权最小二乘法WLS用于DSE(动态状态估计),可能导致非线性状态估计问题。
[0005] 针对上述方法的不足,如何在保证输电线路故障测距方法足够精确的前提下,使得故障定位方法步骤和计算尽量简单,是一个亟待解决的技术问题。

发明内容

[0006] 为解决上述技术问题,本发明提供一种基于卡尔曼滤波器的输电线路故障定位方法,该方法是一种两步骤的集成卡尔曼滤波器(EnKF)方法,在不需要推导雅可比矩阵、降低了实现算法的难度的同时,充分利用了模型非线性和误差统计的全面演化,使故障位置预估更加精确,数据更加稳定。
[0007] 本发明采取的技术方案为:
[0008] 基于卡尔曼滤波器的输电线路故障定位方法,包括以下步骤:
[0009] 步骤1:给定输电线模型和测量值;
[0010] 步骤2:引入故障状态α;
[0011] 步骤3:通过基于卡尔曼滤波器EnKF的动态状态估计DSE过程,得到输电线路故障类型和故障定位的初步估计值;
[0012] 步骤4:将附加信息作为额外的约束,通过附加信息重新运行动态状态估计DS,以执行增强的卡尔曼滤波器EnKF;
[0013] 步骤5:测量冗余的增加,得到对输电线路故障位置准确估计值。
[0014] 上述方法具有高精确度的同时,满足了输电线路故障定位过程对简短步骤的需求。
[0015] 所述步骤1中,
[0016] 给先定输电线模型和测量值:
[0017] xk=f(xk-1,uk-1)+wk-1  (1);
[0018] zk=h(xk,uk)+vk-1  (2);
[0019] 其中,公式(1)(2)中k是步长数;
[0020] 公式(1)中,xk是k时刻输电线路的状态,xk-1是k-1时刻输电线路的状态,uk-1是k-1时刻系统的控制量,f(xk-1,uk-1)是前一时刻状态xk-1、前一时刻控制量uk-1和xk之间的映射关系,wk-1表示对状态量xk处理和测量过程伴随的随机噪声,其均值为零,方差为Qd;
[0021] 公式(2)中,zk是k时刻的测量值,h(xk,uk)是xk,uk和zk之间的映射关系,vk-1表示对测量值zk处理和测量过程伴随的随机噪声,其均值为零,方差为Rd。
[0022] 状态集合表示为 的N×n矩阵;
[0023] 其中,n是状态维数,N是样本数。
[0024] 公式(1),公式(2)就是输电线路动态评估的线性模型微分公式,以k为参数就是离散化。
[0025] 所述步骤2中,
[0026] 当输电线路发生故障时,从故障到终端的距离被视为故障期间不变的新状态α,新输电线模型变成:
[0027]
[0028] 公式(3)中,αk是k时刻的故障状态,xk是k时刻的输电线路状态,f(xk-1,uk-1,ak-1)是xk-1,uk-1,αk-1和xk之间的映射关系,wa,k-1是k-1时刻的过程和测量产生的噪声。
[0029] 当引入故障状态α作为状态时,输电线模型的非线性会增加,故定义一个扩展的状态向量为y=(x,α)T。将公式(1)和(2)重写为:
[0030] yk=f(yk-1,uk-1)+wk-1  (4);
[0031] zk=h(yk,uk)+vk-1  (5);
[0032] 公式(4),(5)中是k步长数;
[0033] 公式(1)中的yk是k时刻输电线路的状态,yk-1是k-1时刻输电线路的状态,f(yk-1,uk-1)是yk-1,uk-1和yk之间的映射关系,uk-1是k-1时刻对系统的控制量,wk-1表示时刻的过程和测量产生的噪声;
[0034] 公式(2)中的zk是k时刻的测量值,h(yk,uk)是yk,uk和zk之间的映射关系,vk-1表示k-1过程和测量产生的噪声。
[0035] 步骤3中,基于EnKF的DSE使用多次的预测校正过程进行分析,如下所示:
[0036] EnKF预测:
[0037]
[0038] 公式(6)中, 表示输电线路在k时刻的状态且由预测步骤计算得到,uk-1是k-1时刻对系统的控制量, 表示 uk-1和 之间的映射关系, 表示k-1时刻的过程和观测产生的噪声且有校正步骤计算得到;
[0039]
[0040] 公式(7)中, 表示时刻的测量值且由预测步骤计算得到,uk是k时刻对系统的控制量, 表示 uk和 之间的映射关系, 表示k-1时刻的过程和观测产生的噪声且有校正步骤计算得到;
[0041]
[0042] 公式(8)中, 表示输电线路在k时刻的状态且由预测步骤计算得到,N是样本数,是yk的共轭矩阵。
[0043]
[0044] 公式(9)中, 表示时刻的测量值且由预测步骤计算得到,N是样本数, 是zk的共轭矩阵。
[0045] EnKF校正:
[0046]
[0047]
[0048]
[0049]
[0050] 公式(10)、(11)、(12)、(13)中,N是样本数, 是误差协方差矩阵的共轭矩阵,是雅克比矩阵的转置矩阵, 表示输电线路在k时刻的状态且由预测步骤计算得到, 表示时刻的测量值且由预测步骤计算得到,Kk是卡尔曼增益。
[0051] 其中,Pk是误差协方差矩阵,Hk是雅克比矩阵,Kk是卡尔曼增益,上标i-表示在预测步骤计算状态和计算观测值。注意,Jacobian矩阵Hk不需要显式地从非线性公式组导出。
[0052] 基于EnKF的动态状态估计DSE过程,由上述基于EnKF的DSE预测校正过程分析基础上,已知故障前和故障后的多段输电线路模型,在出现故障区段的连接点处稍有不同,无故障段连接点的三相电压电流总是相等的。
[0053]
[0054] 公式(14)中, 表示三相输电线上任意一点的左右两端随时间变化的电压值。
[0055]
[0056] 公式(15)中, 表示三相输电线任意一点的左右两端随时间变化的电流值。
[0057] 其中,i≠m,上述电压和电流关系表现了基尔霍夫的电压和电流定律。然而,在故障点只有式(14)所示的电压关系存在,如果故障类型未知,那么当前关系仍然不清楚;
[0058]
[0059] 公式(16)中, 表示三相电路故障点左右两端随时间变化的电压值。
[0060] 将公式(14)-(16)包含在基于EnKF的DSE预测校正过程中,将得到故障定位 的一个近似值,以及故障点处的三相电压的估计值;
[0061] 根据故障点处的三相电压的估计值,能够对故障类型进行分类:
[0062] 通常,故障阻抗越小,出现故障相电压的变化越大;
[0063] 故障类型是通过检查三相电压来确定的,若B和C相电压等于标称值,而A相电压很小,那么这是一个A-G单相接地短路故障;若C相电压等于标称值,而A相和B相的电压值接近标称值的一半,那么这是一个A-B的相间短路故障;若C相电压等于标称值,但是相A和B电压非常小,接近零,,那么它是一个AB-G两相接地的短路故障。
[0064] 如图3(a)、图3(b)、图3(c)、图3(d)所示不同故障类型所对应的电压曲线:图3(a)为单相接地故障电压曲线图;图3(b)为两相相间短路电压曲线图;图3(c)为两相相间短路且接地电压曲线图;图3(d)为三相相间短路且接地电压曲线图。
[0065] 通过步骤3得到了故障位置的近似值 和故障类型;故障位置的近似值 作为增强的EnKF状态估计过程的初始猜测。
[0066] 步骤4中,利用步骤3所得的故障类型,导出了故障点处的附加电流关系式;
[0067] 如果第一步确定故障是A-G故障,则附加的电流关系是:
[0068]
[0069] 公式(17)中, 分别表示故障线路的b相的任意一点左右两端随时间变化的电流值,即B相的电流值为标称值。
[0070]
[0071] 公式(18)中, 分别表示故障线路的c相的任意一点左右两端随时间变化的电流值,即C相的电流值为标称值。
[0072] 当重新运行增强的基于EnKF的DSE时,上述第一步所得到的故障类型和故障点距离的估计值作为额外的约束,作为增加状态估计过程的测量冗余。
[0073] 如果是AB-G两相接地的短路故障,电流关系为:
[0074]
[0075] 公式(19)中, 分别表示故障线路的c相的任意一点左右两端随时间变化的电流值,即C相的电流为标称值,且A,B两相的电压值接近零。
[0076] 如果是A-B相间短路故障,电流关系为:
[0077]
[0078] 公式(20)中, 分别表示故障线路的c相的任意一点左右两端随时间变化的电流值,即C相的电流为标称值,且A,B两相的电压值约为标称值的一半。
[0079] 如果是三相相间短路故障,三相的电压值均接近零。
[0080] 步骤5中,在重新运行增强EnKF的DSE过程时,由于较好的初始状态和增加的测量冗余,实现了对故障位置的更精确的估计,执行增强的EnKF过程,能够给出更精确的故障定位估计值。
[0081] 较好的初始状态指的是输电线路的中间段,由于中间段距离输电线路上任意一点的距离都相对较短,所以距离故障位置也相对较近,所以被称为较好的初始状态。
[0082] 本发明一种基于卡尔曼滤波器的输电线路故障定位方法,不需要预先了解故障位置和故障类型,也能在半周期内准确定位输电线路故障。
[0083] 本发明一种基于卡尔曼滤波器的输电线路故障定位方法,包括运行基于ENKF的方法,使用附加信息重新运行DSE过程来执行增强的ENKF两大步骤。前者得到故障位置的初步估计值和故障类型;后者将得到更精确的故障定位。
[0084] 本发明一种基于卡尔曼滤波器的输电线路故障定位方法,技术效果如下:
[0085] 1:设计了一种简单的两步骤故障定位方法,重复运行DSE过程可以得到更加精确的故障位置的估计值。
[0086] 2:该方法不需要准确预知故障的位置,故障的类型,也不需要导出雅可比矩阵,只需一个半周期的数据窗,降低了方法的计算难度。附图说明
[0087] 图1为本发明基于EnKF方法流程图
[0088] 图2为本发明方法的系统图。
[0089] 图3(a)为单相接地故障电压曲线图;
[0090] 图3(b)为两相相间短路电压曲线图;
[0091] 图3(c)为两相相间短路且接地电压曲线图;
[0092] 图3(d)为三相相间短路且接地电压曲线图。

具体实施方式

[0093] 基于卡尔曼滤波器的输电线路故障定位方法,包括以下步骤:
[0094] 步骤1:给定输电线模型和测量值;
[0095] 步骤2:引入故障状态α;
[0096] 步骤3:通过基于卡尔曼滤波器EnKF的动态状态估计DSE过程,得到输电线路故障类型和故障定位的初步估计值;
[0097] 步骤4:将附加信息作为额外的约束,通过附加信息重新运行动态状态估计DS,以执行增强的卡尔曼滤波器EnKF;
[0098] 步骤5:测量冗余的增加,得到对输电线路故障位置准确估计值。
[0099] 实施例
[0100] 本实施例以500kV、200km传输线上的两组故障场景为研究对象,对基于EnKF的故障定位方法,单端阻抗法,双端阻抗法在输电线路上每20km段的故障定位误差进行比较。上述方法主要包括:运行基于EnKF的方法和使用附加信息重新运行DSE过程来执行增强的EnKF两大步骤。
[0101] 基于EnKF的方法包括给定输电线路模型,给定输电线路故障状态,运行基于EnKF的DSE过程3个步骤;基于EnKF的DSE由反复的EnKF预测和EnKF校正进行总结;基于EnKF的方法能够给出故障类型和故障定位的初步估计值。而执行增强的EnKF过程将给出更精确的故障定位估计值。
[0102] 作为本发明的一种优选方案,基于EnKF方法的流程图如图1所示,具体实施步骤如下:
[0103] (1):给定输电线的模型和测量值,本实施例训练过程中,相关参数的类型及初值如下表1:
[0104] 表1输电线路及等效系统参数表
[0105] l 200-km V 500-kvz1,Line (0.01273+j0.3520)Ω/km c1,Line 12.74nF/km
z0,Line (0.3864+j1.5557)Ω/km c0,Line 7.751nF/km
z1,SL (1.1547+j11.5473)Ω/km z1,SR 2z1,SL
z0,SL (2.8868+j20.2078)Ω/km z0,SR 2z0,SL
AngSL 0° AngSR -30°
[0106] (2):输入故障状态,本实施例中输电线路的故障状态设定如下:
[0107] 第一种情况是单相(A相)接地短路故障,故障阻抗为15Ω;
[0108] 第二种是相间(B-C)短路故障,故障阻抗为0.01Ω;
[0109] (3):运行基于EnKF的DSE过程。
[0110] 基于EnKF的方法具体运行步骤如下,流程图如图1所示:
[0111] (1)、给定输电线模型和测量值:
[0112] xk=f(xk-1,uk-1)+wk-1  (1)
[0113] zk=h(xk,uk)+vk-1  (2)
[0114] 式中,x是状态,u是输入,z是测量值,k是步长数,w,v是过程,Qd、Rd、f和h的零均值和协方差观测噪声是映射或函数。集合表示为x=[x1...xN]形式的n×N矩阵,其中n是状态维数,N是样本数。公式(1),(2)是基于DSE的输电线路故障定位中给出的线性模型微分公式的离散化。
[0115] (2)、当输电线路发生故障时,将从故障到终端的距离视为故障期间不变的新状态α。
[0116] 新模型变成:
[0117]
[0118] 当引入α作为状态时,模型的非线性会增加,故定义一个扩展的状态向量为y=(x,α)T。将公式(1)和(2)重写为:
[0119] yk=f(yk-1,uk-1)+wk-1  (4)
[0120] zk=h(yk,uk)+vk-1  (5)
[0121] (3)、将基于ENKF的DSE使用多次的预测校正过程进行总结,如下所示:
[0122] ENKF预测:
[0123]
[0124]
[0125]
[0126]
[0127] ENKF校正:
[0128]
[0129]
[0130]
[0131]
[0132] 式中,Pk是误差协方差矩阵,Hk是雅克比矩阵,Kk是卡尔曼增益,上标i-表示在预测步骤计算状态和计算观测值。Jacobian矩阵Hk不需要显式地从非线性公式组导出。
[0133] 执行增强的EnKF的具体实施步骤如下:
[0134] (1)根据判断出的故障类型给出附加的电流关系作为附加信息;
[0135] (2)重新运行基于EnKF的DSE过程以得到更精确的故障定位估计值。
[0136] 实施例最终测试结果如下:
[0137] 表2 A相接地短路故障情形比较表
[0138]
[0139]
[0140] 由表2得到:EnKF方法在用于发生单相接地短路故障的故障点距离的估计时,平均误差最小,因此对比传统的单端阻抗和双端阻抗法,EnKF更适用于发生单相接地短路故障的故障点距离的估计。
[0141] 表3 B-C相间短路故障情形比较表
[0142]
[0143] 由表3得到:EnKF方法在用于发生双相相间短路故障的故障点距离的估计时,平均误差最小,因此对比传统的单端阻抗和双端阻抗法,EnKF更适用于发生双相相间短路故障的故障点距离的估计。
高效检索全球专利

专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台,是提供专利分析,专利查询,专利检索等数据服务功能的知识产权数据服务商。

我们的产品包含105个国家的1.26亿组数据,免费查、免费专利分析。

申请试用

分析报告

专利汇分析报告产品可以对行业情报数据进行梳理分析,涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度,系统通过AI智能系统对图表进行解读,只需1分钟,一键生成行业专利分析报告。

申请试用

QQ群二维码
意见反馈