最新消息: USBMI致力于为网友们分享Windows、安卓、IOS等主流手机系统相关的资讯以及评测、同时提供相关教程、应用、软件下载等服务。

灯泡贯流式双调节水轮发电机组调速系统建模与仿真

IT圈 admin 33浏览 0评论

2024年6月13日发(作者:仲睿才)

第33卷第4期 

2010年8月 

水电站机电技术 

V01.33 No.4 

Mechanical&Electrical Technique of Hydropower Station 

Aug.2010 23 

灯泡贯流式双调节水轮发电机组调速系统建模与仿真 

彭天波 

(湖北省电力试验研究院,湖北武汉430077) 

摘要:以灯泡贯流式双调节水轮发电机组为例,首先进行调速系统参数现场测试,并经分析、整理后得出了水轮机 

调节系统各环节的重要参数,之后分别采用最小二乘法(LS)和基因遗传算法(GA)对原动机及调节系统建模和参数 

辨识,得到非线性”自定义”模型,并与电网稳定计算通用模型PSASP进行了比较分析,从而建立了水轮机调节系统 

在并网状态下的模型和数据,为电力系统仿真研究提供与实际系统状况相吻合的水轮机调节系统模型,满足了电力 

系统稳定计算的要求。 

关键词:灯泡贯流式机组;参数测试;调速系统建模;仿真 

中图分类号:TV734.4 文献标识码:B 文章编号:1672—5387(2010)04—0023—05 

某航电枢纽工程1号机组为灯泡贯流式水轮发 

蜗壳未充水,调速器处于自动运行方式;切除人 

频率给定为额定值;机组 

电机组,由东芝水电设备航州)有限公司生产,其主 

工死区;置增益为整定值,

模拟并网(断路器位置信号短接)。 

要设备参数如下: 

水轮机型号GZ(T ̄5003)一WP-691;最大功率 

1.1.2试验方法 

16.980 MW;额定功率15.432 MW;最大水头8.4 m;设 

计水头4.7m;最小水头15m;设i十流量 ̄<366.47m3/s。发 

由调速器测试仪给调速器提供额定的机频输入 

信号,以开度给定将导叶接力器调整到50%行程附 

近。然后,将调速器处于自动运行方式,升高或降低 

频率使接力器全关或全开,调整频率信号值,使之按 

电机型号SFWG15—84—7300;额定容量16.67 MVA; 

额定功率15 MW;额定电压10.5 kV;额定电流916.4 A; 

额定频率50 Hz;额定转速71.4 r/min;飞逸转速 

个方向逐次升高或降低,在导叶接力器每次变化 

240 r/min;功率因素0.9(滞后);转动惯量3 000 t・m ; 

惯性时间常数 6.13(额定水头)/2.06(最大水头)s。 

调速器调节参数调整范围如下: 

稳定后,用调速器测试仪记录该次信号频率值及相 

应的接力器行程、浆叶行程,分别绘制频率升高或降 

低的静态特性曲线。每条曲线在接力器行程10% 

比例增益 :0.5—20;积分增益Ki:0.025~10; 

微分增益Kd:0~40;永态转差率b。调整范围: 

90%的范围内,测点不少于1O点。 

设定bt=5%,Td=2 s,Tn=0,bp=6%,分三个水头 

4 in,额定水头4.7 in,最小水头1.5 m) 

0~10%;人工转速死区调整范围:0~±0-3 Hz;电气开 

(最大水头8.

限调整范围:0~100%;“频率给定”(数字给定)调整 

进行试验。 

范围:45~55 Hz;“功率给定”(数字给定)调整范围: 

0 ̄IO0%。 

1.1.3试验结果 

试验结果见表1,额定水头下的试验曲线见图1。 

表1调速器静特性试验结果 

调速器液压控制系统参数如下: 

型号:PZWST150—6.3一Stars;额定工作压力: 

6.3 MPa;主配压阀直径:150 mm;主配压阀最大行 

程:15 nl/n;滤芯精度:10/J,m(0120D010BN/HC);差 

压发讯器发讯压差:0.2 MPa。 

1调速器建模参数测试条件、测试方法与测 

试结果 

1.1调速器静特性试验(测量印、jj(、 ̄L等参数) 

1.1.1试验条件 

收稿日期:2010—05-19 

作者简介:彭天波(1966一),男,硕士,高级工程师,从事发电厂经济运 

行与控制工作。 

24 水电站机电技术 第33卷 

验组号标称偏标枥称一定差标偏一标称一^^宴测实冽 ~ 一 呲 砩一i ~ 8孵阳 骗8孤加加如r‘ 坫;8,0‘ ∞8}宝 

\ 

\ \ 。} 

 ’I 

lll l l\ -_ 

\ \●_ 

. 。 

)0  10 O0 20 0O 3O 00 40 D0 60 00 8D 00 7O 00 e0 0 0 q0 00 1 C 

图1额定水头4.7ITI调速器静态特性 

1.2频率扰动试验(测量 、 、丁n等参数) 

1.2.1试验条件 

蜗壳未充水,调速器处于自动运行方式;切除人 

工死区;置增益为整定值,bp=0,频率给定为额定值; 

机组模拟并网(断路器位置信号短接)。分别在空载、 

并网一次调频模式下进行。 

1.2.2试验方法 

由调速器测试仪给调速器提供额定的机频输人 

信号,以开度给定将导叶接力器调整到5%行程附 

近。断开主接反馈,用调速器测试仪改变机组输入频 

率,记录频率改变后、PID输出、导叶接力器行程变 

化过程曲线。试验茁凋速器参数设定为 =0,bt=20%, 

8 s,Tn=0.2 s和bD=0,bt=20%,Td=8 S,Tn=0 s两 

种参数组合下进行。 

1.2.3试验结果 

测试结果汇总见表2。以第4组bp=0,bt=20%, 

8 S,Tn=0.2 S为代表的试验录波图见图2所示。 

表2调速器控制参数测试结果 

试 bt(%) Yd(s) 7h(s) 

注:表中标定偏差=(K-S3/S,K为实测值,s为标称值。 

1 679.. .. ii iii  ii i 

0 0。 

图2第4组频率扰动试验 

1.3接力器关闭规律试验(测量主接最短关闭时间 

丁s、最短开机时间 ) 

测量主接最短关闭时间 、最短开机时间 

结果见表3。 

表3接力器关闭规律试验结果 

测 

簇 

1.4测量人工死区试验 

1.4.1试验条件 

蜗壳未充水,调速器处于自动运行方式。 

1.4.2试验方法 

用阶跃频率信号法校验调速器的转速死区,把 

测试仪输出的频率信号整定为50 Hz,于额定频率 

基础上每步0.01 Hz(或0.005 Hz)施加正(或负)的 

阶跃频率偏差信号,并逐一录波;直到接力器开始产 

生与此信号对应的改变,此时的频率累积偏差信号 

即为人工死区。 

1.4.3试验结果 

经过测试调速系统人工死区为0.05 Hz。 

1.5频率及开度给定扰动试验(测量接力器反应时 

间 、桨叶反应时间 ) 

1.5.1试验条件 

蜗壳未充水,调速器处于自动运行方式;切除人 

工死区;置增益为整定值,频率给定为额定值;机组 

模拟并网(断路器位置信号短接)。断开主接反馈。 

1.5.2试验方法 

由调速器测试仪给调速器提供额定的机频输入 

信号,以开度给定将导叶接力器行程调整到拐点以 

上,给定频率或开度扰动,用调速器测试仪记录频率 

或开度给定变化后导叶接力器行程、桨叶行程变化 

过程曲线。 

1.5-3试验结果 

试验结果见表4、表5。 

第4期 彭天波:灯泡贯流式双调节水轮发电机组调速系统建模与仿真 25 

表4接力器反应时间丁y实测值 

表5桨叶反应时间 实测值 

1.00 

1.00 

接力器开度一3% 

接力器开度+3% 

1.6并网频率扰动试验(测量动态响应时间) 

1.6.1试验条件 

机组带负荷稳定运行,调速器处于自动运行方 

式。 

1.6.2试验方法 

机组带70%额定负荷(10 Mw)运行,用调速器 

测试仪改变输入调速器的频率,模拟改变电网频率。 

把测试仪输出的额定频率信号接人调速器的机频信 

号输入端,分别于额定网频基础上施加正负阶跃偏 

差频率信号(从0.1 Hz、一0.1 Hz;0.2 Hz、一0.2 Hz;到 

0.3 Hz、一0.3 Hz,每个模拟网频信号持续60 s以上 

(以接力器稳定时间为准),频率变化范围49.70— 

50.30 Hz);用调速器测试仪分别对模拟网频频率、 

有功功率和导叶开度、桨叶角度信号进行录波。 

试验实际水头5.4 m,仿真水头6.0 m,开度模 

式,bp=4%,bt=20%,Td=2s,Tn=8s。 

1.6.3试验结果 

试验数据应用于调速系统建模,其中-I-0.2 Hz 

的频率扰动试验示波图见图3、图4所示 

 ●

一 ●●● 

。● 

曲驯:一棚蛆臻奉 一■叶主接 一紫 主接 一膏功功奉 一 嘲氟奉 

图3+0.2Hz频率扰动试验 

|| 

__‘ 岍一 r  :

。 

。 

00 8 70 1740 2610 34 80 43 50 52 20 ̄43∞69∞78 30 87 

图4—0.2Hz频率扰动试验 

2水轮机调速系统建模 

现场试验完成以后,通过对试验数据进行分析、 

整理,得出了水轮机调节系统各环节的重要参数。分 

别采用最小二乘法(LS)和基因遗传算法(GA)对原 

动机及调节系统建模和参数辨识,得到非线性“自定 

义”模型,并与电网稳定计算通用模型PSASP进行 

了比较分析。 

2.1自定义调速器模型、参数及仿真 

水轮机调节系统在并网状态下的数学模型(自 

定义模型)如图5一图7所示。 

图5调速器导叶电气部分(控制器)模型框图 

图6调速器导叶液压系统模型框图 

H 

图7调速器桨叶控制系统模型框图 

水电站机电技术 第33卷 

模型中符号含义和参数见表6。 

表6自定义调速器模型符号含义与参数 

2.2自定义机组一引水道模型、参数及仿真 

水轮机机组一引水道在并网状态下的数学模型 

(自定义模型)如图10所示。 

图1O机组一引水道模型框图 

模型中符号含义和参数见表7,仿真效果见图 

l1、图12所示。 

表7自定义机组一引水道模型符号含义与参数 

名称 

初始导叶开度 

初始功率 

初始桨叶开度 

最小二乘参数辨 

识算法(LS) 

基因遗传参数辨 

识算法(GA) 

符号 

rg 

Pg 

Cg 

数值 单位 

P.u. 

P.u. 

P.u 

传递函数△尸(s)/Ay(s)= 

(一5.00s+1.69)/(5.38s 2+2.00s+1.oo) 

传 递 函 数 △ s)/ △ 

(-4.98s+1.68)/(1.03s 2+3.47s+1。oo) 

s)= 

得到模型后进行仿真计算。调速器导叶开度仿真 

结果见图8,调速器桨叶开度仿真结果曲线见图9。 

图11机组一引水道LS法辨识效果图 

时闻 

图8调速器模型导叶开度仿真与实测过程比较 

图l2机组一引水道GA法辨识效果图 

2.3自定义模型仿真输出与实测动态过程 匕较 

采用Ls及GA法所得自定义模型仿真输出与 

实测动态过程比较见图13、图14所示。 

图9调速器模型桨叶开度仿真与实测过程比较 

第4期 彭天波:灯泡贯流式双调节水轮发电机组调速系统建模与仿真 

表8 PSASP模型符号含义与参数 

27 

《 

避 

l 

100 

80 

篡60 

40 

图13 Ls法自定义模型仿真输出与实测动态过程比较 

图14 GA法自定义模型仿真输出与实测动态过程比较 图16 IJs法自定义模型与PSASP模型的动态过程比较 

2.4水轮机调节系统在并网状态下的自定义模型与 

PSASP模型的比较 

PSASP模型见图15所示,其模型符号含义及参 

数见表8,本次在1号机组建立的自定义模型与 

PSASP模型的动态过程比较见图16、图17。 

图17 GA法自定义模型与PSASP模型的动态过程比较 

3结论与建议 

根据该电站1号机组的现场实测参数和数据, 

建立了1号机组水轮机调节系统在并网状态下的数 

学模型,并进行参数辨识。仿真表明所建模型能较好 

地反映系统在并网状态下对频率变化的动态响应, 

并与电网稳定计算通用模型PSASP进行了比较,可 

以看出,自定义模型更能反映原动机及调节系统实 

际情况,在电力系统稳定计算中采用效果显著。 

图15 PSASP模型框图 

2024年6月13日发(作者:仲睿才)

第33卷第4期 

2010年8月 

水电站机电技术 

V01.33 No.4 

Mechanical&Electrical Technique of Hydropower Station 

Aug.2010 23 

灯泡贯流式双调节水轮发电机组调速系统建模与仿真 

彭天波 

(湖北省电力试验研究院,湖北武汉430077) 

摘要:以灯泡贯流式双调节水轮发电机组为例,首先进行调速系统参数现场测试,并经分析、整理后得出了水轮机 

调节系统各环节的重要参数,之后分别采用最小二乘法(LS)和基因遗传算法(GA)对原动机及调节系统建模和参数 

辨识,得到非线性”自定义”模型,并与电网稳定计算通用模型PSASP进行了比较分析,从而建立了水轮机调节系统 

在并网状态下的模型和数据,为电力系统仿真研究提供与实际系统状况相吻合的水轮机调节系统模型,满足了电力 

系统稳定计算的要求。 

关键词:灯泡贯流式机组;参数测试;调速系统建模;仿真 

中图分类号:TV734.4 文献标识码:B 文章编号:1672—5387(2010)04—0023—05 

某航电枢纽工程1号机组为灯泡贯流式水轮发 

蜗壳未充水,调速器处于自动运行方式;切除人 

频率给定为额定值;机组 

电机组,由东芝水电设备航州)有限公司生产,其主 

工死区;置增益为整定值,

模拟并网(断路器位置信号短接)。 

要设备参数如下: 

水轮机型号GZ(T ̄5003)一WP-691;最大功率 

1.1.2试验方法 

16.980 MW;额定功率15.432 MW;最大水头8.4 m;设 

计水头4.7m;最小水头15m;设i十流量 ̄<366.47m3/s。发 

由调速器测试仪给调速器提供额定的机频输入 

信号,以开度给定将导叶接力器调整到50%行程附 

近。然后,将调速器处于自动运行方式,升高或降低 

频率使接力器全关或全开,调整频率信号值,使之按 

电机型号SFWG15—84—7300;额定容量16.67 MVA; 

额定功率15 MW;额定电压10.5 kV;额定电流916.4 A; 

额定频率50 Hz;额定转速71.4 r/min;飞逸转速 

个方向逐次升高或降低,在导叶接力器每次变化 

240 r/min;功率因素0.9(滞后);转动惯量3 000 t・m ; 

惯性时间常数 6.13(额定水头)/2.06(最大水头)s。 

调速器调节参数调整范围如下: 

稳定后,用调速器测试仪记录该次信号频率值及相 

应的接力器行程、浆叶行程,分别绘制频率升高或降 

低的静态特性曲线。每条曲线在接力器行程10% 

比例增益 :0.5—20;积分增益Ki:0.025~10; 

微分增益Kd:0~40;永态转差率b。调整范围: 

90%的范围内,测点不少于1O点。 

设定bt=5%,Td=2 s,Tn=0,bp=6%,分三个水头 

4 in,额定水头4.7 in,最小水头1.5 m) 

0~10%;人工转速死区调整范围:0~±0-3 Hz;电气开 

(最大水头8.

限调整范围:0~100%;“频率给定”(数字给定)调整 

进行试验。 

范围:45~55 Hz;“功率给定”(数字给定)调整范围: 

0 ̄IO0%。 

1.1.3试验结果 

试验结果见表1,额定水头下的试验曲线见图1。 

表1调速器静特性试验结果 

调速器液压控制系统参数如下: 

型号:PZWST150—6.3一Stars;额定工作压力: 

6.3 MPa;主配压阀直径:150 mm;主配压阀最大行 

程:15 nl/n;滤芯精度:10/J,m(0120D010BN/HC);差 

压发讯器发讯压差:0.2 MPa。 

1调速器建模参数测试条件、测试方法与测 

试结果 

1.1调速器静特性试验(测量印、jj(、 ̄L等参数) 

1.1.1试验条件 

收稿日期:2010—05-19 

作者简介:彭天波(1966一),男,硕士,高级工程师,从事发电厂经济运 

行与控制工作。 

24 水电站机电技术 第33卷 

验组号标称偏标枥称一定差标偏一标称一^^宴测实冽 ~ 一 呲 砩一i ~ 8孵阳 骗8孤加加如r‘ 坫;8,0‘ ∞8}宝 

\ 

\ \ 。} 

 ’I 

lll l l\ -_ 

\ \●_ 

. 。 

)0  10 O0 20 0O 3O 00 40 D0 60 00 8D 00 7O 00 e0 0 0 q0 00 1 C 

图1额定水头4.7ITI调速器静态特性 

1.2频率扰动试验(测量 、 、丁n等参数) 

1.2.1试验条件 

蜗壳未充水,调速器处于自动运行方式;切除人 

工死区;置增益为整定值,bp=0,频率给定为额定值; 

机组模拟并网(断路器位置信号短接)。分别在空载、 

并网一次调频模式下进行。 

1.2.2试验方法 

由调速器测试仪给调速器提供额定的机频输人 

信号,以开度给定将导叶接力器调整到5%行程附 

近。断开主接反馈,用调速器测试仪改变机组输入频 

率,记录频率改变后、PID输出、导叶接力器行程变 

化过程曲线。试验茁凋速器参数设定为 =0,bt=20%, 

8 s,Tn=0.2 s和bD=0,bt=20%,Td=8 S,Tn=0 s两 

种参数组合下进行。 

1.2.3试验结果 

测试结果汇总见表2。以第4组bp=0,bt=20%, 

8 S,Tn=0.2 S为代表的试验录波图见图2所示。 

表2调速器控制参数测试结果 

试 bt(%) Yd(s) 7h(s) 

注:表中标定偏差=(K-S3/S,K为实测值,s为标称值。 

1 679.. .. ii iii  ii i 

0 0。 

图2第4组频率扰动试验 

1.3接力器关闭规律试验(测量主接最短关闭时间 

丁s、最短开机时间 ) 

测量主接最短关闭时间 、最短开机时间 

结果见表3。 

表3接力器关闭规律试验结果 

测 

簇 

1.4测量人工死区试验 

1.4.1试验条件 

蜗壳未充水,调速器处于自动运行方式。 

1.4.2试验方法 

用阶跃频率信号法校验调速器的转速死区,把 

测试仪输出的频率信号整定为50 Hz,于额定频率 

基础上每步0.01 Hz(或0.005 Hz)施加正(或负)的 

阶跃频率偏差信号,并逐一录波;直到接力器开始产 

生与此信号对应的改变,此时的频率累积偏差信号 

即为人工死区。 

1.4.3试验结果 

经过测试调速系统人工死区为0.05 Hz。 

1.5频率及开度给定扰动试验(测量接力器反应时 

间 、桨叶反应时间 ) 

1.5.1试验条件 

蜗壳未充水,调速器处于自动运行方式;切除人 

工死区;置增益为整定值,频率给定为额定值;机组 

模拟并网(断路器位置信号短接)。断开主接反馈。 

1.5.2试验方法 

由调速器测试仪给调速器提供额定的机频输入 

信号,以开度给定将导叶接力器行程调整到拐点以 

上,给定频率或开度扰动,用调速器测试仪记录频率 

或开度给定变化后导叶接力器行程、桨叶行程变化 

过程曲线。 

1.5-3试验结果 

试验结果见表4、表5。 

第4期 彭天波:灯泡贯流式双调节水轮发电机组调速系统建模与仿真 25 

表4接力器反应时间丁y实测值 

表5桨叶反应时间 实测值 

1.00 

1.00 

接力器开度一3% 

接力器开度+3% 

1.6并网频率扰动试验(测量动态响应时间) 

1.6.1试验条件 

机组带负荷稳定运行,调速器处于自动运行方 

式。 

1.6.2试验方法 

机组带70%额定负荷(10 Mw)运行,用调速器 

测试仪改变输入调速器的频率,模拟改变电网频率。 

把测试仪输出的额定频率信号接人调速器的机频信 

号输入端,分别于额定网频基础上施加正负阶跃偏 

差频率信号(从0.1 Hz、一0.1 Hz;0.2 Hz、一0.2 Hz;到 

0.3 Hz、一0.3 Hz,每个模拟网频信号持续60 s以上 

(以接力器稳定时间为准),频率变化范围49.70— 

50.30 Hz);用调速器测试仪分别对模拟网频频率、 

有功功率和导叶开度、桨叶角度信号进行录波。 

试验实际水头5.4 m,仿真水头6.0 m,开度模 

式,bp=4%,bt=20%,Td=2s,Tn=8s。 

1.6.3试验结果 

试验数据应用于调速系统建模,其中-I-0.2 Hz 

的频率扰动试验示波图见图3、图4所示 

 ●

一 ●●● 

。● 

曲驯:一棚蛆臻奉 一■叶主接 一紫 主接 一膏功功奉 一 嘲氟奉 

图3+0.2Hz频率扰动试验 

|| 

__‘ 岍一 r  :

。 

。 

00 8 70 1740 2610 34 80 43 50 52 20 ̄43∞69∞78 30 87 

图4—0.2Hz频率扰动试验 

2水轮机调速系统建模 

现场试验完成以后,通过对试验数据进行分析、 

整理,得出了水轮机调节系统各环节的重要参数。分 

别采用最小二乘法(LS)和基因遗传算法(GA)对原 

动机及调节系统建模和参数辨识,得到非线性“自定 

义”模型,并与电网稳定计算通用模型PSASP进行 

了比较分析。 

2.1自定义调速器模型、参数及仿真 

水轮机调节系统在并网状态下的数学模型(自 

定义模型)如图5一图7所示。 

图5调速器导叶电气部分(控制器)模型框图 

图6调速器导叶液压系统模型框图 

H 

图7调速器桨叶控制系统模型框图 

水电站机电技术 第33卷 

模型中符号含义和参数见表6。 

表6自定义调速器模型符号含义与参数 

2.2自定义机组一引水道模型、参数及仿真 

水轮机机组一引水道在并网状态下的数学模型 

(自定义模型)如图10所示。 

图1O机组一引水道模型框图 

模型中符号含义和参数见表7,仿真效果见图 

l1、图12所示。 

表7自定义机组一引水道模型符号含义与参数 

名称 

初始导叶开度 

初始功率 

初始桨叶开度 

最小二乘参数辨 

识算法(LS) 

基因遗传参数辨 

识算法(GA) 

符号 

rg 

Pg 

Cg 

数值 单位 

P.u. 

P.u. 

P.u 

传递函数△尸(s)/Ay(s)= 

(一5.00s+1.69)/(5.38s 2+2.00s+1.oo) 

传 递 函 数 △ s)/ △ 

(-4.98s+1.68)/(1.03s 2+3.47s+1。oo) 

s)= 

得到模型后进行仿真计算。调速器导叶开度仿真 

结果见图8,调速器桨叶开度仿真结果曲线见图9。 

图11机组一引水道LS法辨识效果图 

时闻 

图8调速器模型导叶开度仿真与实测过程比较 

图l2机组一引水道GA法辨识效果图 

2.3自定义模型仿真输出与实测动态过程 匕较 

采用Ls及GA法所得自定义模型仿真输出与 

实测动态过程比较见图13、图14所示。 

图9调速器模型桨叶开度仿真与实测过程比较 

第4期 彭天波:灯泡贯流式双调节水轮发电机组调速系统建模与仿真 

表8 PSASP模型符号含义与参数 

27 

《 

避 

l 

100 

80 

篡60 

40 

图13 Ls法自定义模型仿真输出与实测动态过程比较 

图14 GA法自定义模型仿真输出与实测动态过程比较 图16 IJs法自定义模型与PSASP模型的动态过程比较 

2.4水轮机调节系统在并网状态下的自定义模型与 

PSASP模型的比较 

PSASP模型见图15所示,其模型符号含义及参 

数见表8,本次在1号机组建立的自定义模型与 

PSASP模型的动态过程比较见图16、图17。 

图17 GA法自定义模型与PSASP模型的动态过程比较 

3结论与建议 

根据该电站1号机组的现场实测参数和数据, 

建立了1号机组水轮机调节系统在并网状态下的数 

学模型,并进行参数辨识。仿真表明所建模型能较好 

地反映系统在并网状态下对频率变化的动态响应, 

并与电网稳定计算通用模型PSASP进行了比较,可 

以看出,自定义模型更能反映原动机及调节系统实 

际情况,在电力系统稳定计算中采用效果显著。 

图15 PSASP模型框图 

发布评论

评论列表 (0)

  1. 暂无评论