资源描述:
第 l 期 2 0 1 6年 1月 锅炉制造 BOI L ER MANUFACTURI NG No . 1 J a n . 2 01 6 不确定度原理在火电厂管道热效率测算 中的应用 曾大海, 石奇光, 吴 文, 张静秋, 李瑞雪 上海电力学院能源与机械工程学院, 上海 2 0 0 0 9 0 摘要 火电厂全面节能管理及其机组调峰 , 开始重视管道系统热效率对 电厂煤耗的影响。依托国家能源局 新颁的 火力发电厂能量平衡导则 第 3部分 热平衡 D L / T 6 0 6 . 3 2 0 1 4 行业标准中关于火力发电厂管道 热效率的测算方法, 应用测量不确定度原理, 研究管道热效率不确定度的评定方法, 并以 3 5 0 M W 火电机组 D C S在线监测数据为例, 用常规焓差法定量计算和定性分析了影响管道热效率不确定度的主要因素, 对准确 评价火力发电厂供电煤耗率及建立统一的火电机组管道热效率测试导则提供参考。 关键词 火电厂; 管道热效率; 不确定度 ; 分析原理 中图分类号 T K 2 1 2 文献标 识码 A 文章编号 C N 2 31 2 4 9 2 0 1 5 0 5 0 0 1 0 0 7 Ap p l i c a t i o n o f Un c e r t a i n t y Effi c i e n c y Te s t p r i nc i p l e i n Pi p i n g Th e r m a l f l0 r Po we r Pl a n t z e n g D a h a i , S h i Q i g u a n g ,W u W e n , Z h a n g J i n g q i u , L i R u i x u e C o l l e g e o f E n e r g y a n d M e c h a n i c a l E n g i n e e r i n g , S h a n g h a i U n i v e r s i t y o f E l e c t r i c P o w e r , S h a n g h a i 2 0 0 0 9 0 , C h i n a Abs t r a c t Fo r c o mp r e h e n s i v e e n e r g y s a v i n g ma na g e me n t o f c o a lfi r e d p o we r pl a n t a n d l o a d r e g u l a t i o n,t h e i nfl u e n c e o f pi pe s y s t e m t h e r ma l e f f i c i e n c y o f t h e p o we r p l a n t c o a l c o ns ump t i o n h a d a t t r a c t e d t h e a t t e n t i o n .I n t h i s p a p e r ,t h e me a s u r i n g a n d c alc u l a t i o n me t h o d s o f p o w e r p l a n t p i p e l i n e t h e r - ma l e f f i c i e n c y h a s b e e n u s e d,b a s e d o n“G u i d e f o r e n e r gy b a l a n c e o f t h e r ma l p o we r p l a n t p a r t 3 t h e r m a l e n e r g y b a l a n c e ” D L / T 6 0 6 . 3 2 0 1 4 w h i c h w a s p r o mu l g a t e d b y t h e N a t i o n a l E n e r gy A d mi n i s t r a t i o n.Th e p rinc i p l e o f me a s u r e me n t un c e rta i nt y h a s b e e n r e r e nc e d.Un c e rta i n t y e v a l ua t i o n a n d a naly s i s me t h o d o f p i p i n g t h e r ma l e f f i c i e n c y i s r e s e a r c h e d. T a k i n g t h e 3 5 0MW po we r p l a n t for e x a mp l e,t h e d a t a o f DC S i s u s e d for c a l c u l a t i o n a n d a n a l y s i s t o e s t i ma t e t h e f a c t o r s o f t h e u n c e r t a i n - t y o f p i p e l i n e t h e r mal e ffic i e n c y b y e n t h a l p y d i f f e r e n c e me t h o d . T h i s c a n p r o v i d e t h e r e f e r e n c e for a c c u r a t e a s s e s s me n t o f t h e s u pp l y c o a l c o n s u mp t i o n o f p o we r p l a n t a nd t h e e s t a b l i s h me n t o f u n i fie d p i p i n g t h e r ma l e f f i c i e n c y t e s t g u i d e l i n e o f p o we r p l a n t . Ke y wo r dsPo we r p l a n t ;Pi p i n g t he rm a l e ffi c i e n c y;Un c e rta i n t y;An aly s i s t h e o r y 0 引 言 在火力发电厂进行全厂热经济性评价 时, 管 道热力系统是必要的组成部分。当前 , 火 电企业 正全面深化节能管理及其经常性的调峰运行 , 更 加重视管道系统热效率对电厂煤耗的影响。以往 的工程实 际中一般 都将管道热效率 取为经验值 99 % lL l J ,事实上 , 在不 同负荷状态下 管道热效 率 是变化的 2 - 4 1 。不同容量 、 一次再热或两次再 热 不同循环 的火 电机组管道热效率亦不同。火电厂 管道热效率测算的准确性对电厂煤耗率的评价影 响很大 , 管道效率变化 1 %, 机组发 电煤耗变化约 收稿 日期 2 0 1 5 1 01 1 作者简介 曾大海 1 9 8 7 一 , 男, 四川成都人, 硕士, 研究方向为火电厂热力系统节能。 学兔兔 w w w .x u e t u t u .c o m 第 1期 曾大海, 等 不确定度原理在火电厂管道热效率测算中的应用 ‘1 1 3 g c e / k W h 。同时 , 相对于火力发 电厂 的锅炉和 汽轮机已有工程实用统一的热效率实验准则, 然 而, 火力发电厂管道热力系统未有统一的热效率 测试标准。 本文依托新颁的 火力发 电厂能量平衡 导则 第3部分 热平衡 D L / T 6 0 6 . 3 2 0 1 4 以下简 称 导则 行业标准中关于火力发电厂管道热效 率的测算方法 , 应用测量不确定度 相关原理 , 建立了火电厂管道热效率不确 定度评定分析模 型, 给出了管道热效率不确定度 计算的完整流程 , 并 以 3 5 0 MW 火 电机组 D C S 分 散控制系统 在线 数据测试管道热效率 为例 , 阐 述了管道热效率相关参数的采集及其数据处理方 法, 对正、 反平衡管道热效率测算结果进行不确定 度评定, 探讨了降低管道热效率不确定度和提高 管道热效率的技术措施, 为火力发电厂开展管道 热效率在线试验和准确评价机组供电煤耗率提供 参考, 为研究建立统一的火电机组管道热效率测 试导则亦有现实意义 。 1 管道热效率与全厂热经济性 的关 系 火力发电厂热经济性指标主要有全厂热效率 全厂热耗率 H R 和供电煤耗率 b 6, 与管道 热效率 叼 有着密切联系 , 计算方法为 叼 印 7 7 叼 1 职 k J / k W . h 叼c p 2 b g 标准煤 / k w h , , 一 A / 3 式中, 叼 印 为全厂热效率, %; 叼 为锅炉热效 率, %; 叼 管道热效率, %; 叼 为汽轮发电机组热 效率 , %; A为厂用电率 , %。 2 管道热效率的不确定度 管道热效率测量的不确定度是表征管道热效 率测量值分散性的一个量化数值, 这个数值带有 “4 - ” 号 , 它决定了管道热效率真值的范围。 2 . 1 管道热效率测量不确定度评定方法 在工程实践 中不确定度分析代替误差分析 , 对测试结果的评判 已得 到较为广泛的应用 。 1 9 9 5年 , 国际标准化组织 I S O 起 草了测量不确 定度表示指南 G u i d e t o t h e E x p r e s s i o n o f U n c e r - t a i n t y i n M e a s u r e m e n t 简称 G U M 卜 。我国于 1 9 9 9年发 布了 测量 不确定度评定 与表示 J J F 1 0 5 91 9 9 9 , 并 于 2 0 1 2年 进行 了修 订 , 使不确定度的分析规范化 。 不确定的评定综合了A、 类不确定度及不 确定度的合成 , 并通过扩展不确定度表示最终 的 测量结果。 在管道热效率的不确定度评定中, A类不确 定度是通过多次测量 以减少直接测量 参数 的误 差 , 其最终结果是通过统计学 中的标准差 s 进行 表征; B类不确定度是测量误差 限 △与测量参数 的误差统计学分布取值的比值, 是采用非统计学 方法所得到的。 合成不确定度 “ 卵 是管道热效率模 型中 可直接澳 0 量参数 的灵敏度与对应的标准不确定度 共 同表征的不确定度评定最终结果的重要参数 。 其中参数的灵敏度C 是测量模型对该参数的偏导 数 , 测量参数的标准不确定度是包 含 A、 类不 确定度 的综合参数。 扩 展 不 确 定 度 是 管 道 合 成 不 确 定 度 卵 与包含因子 k 的乘积, 是管道热效率可能 值包含的半区间, 一般而言, k 2 , 当 k 取其它值 时, 需要特别说明 ; 管道热效率测量不确定度 的最 终表达形式为 [ 7 一 , 叼 ]。 2 . 2 管道热效率测量不确定度分析模型 火 电厂管道热效率是指汽轮机的热耗量占锅 炉热负荷的百分 比, 可通过正平衡法或反平衡法 计算所得 , 两种计算方法所得结果_致。火 电厂 机组管道热力系统示意图如图 1 所示。 图 1 具有一次再热火 电机组 管道热 力系统 示意图 一 高温过热器 出 E l 蒸汽焓 , k J / k g ; D 一高 温过热器出口蒸汽流量 , t / h ; 一 主蒸汽流量 机 侧 , t/ h ; 一主蒸汽焓 机侧 , k J / k g ; D 一再热 蒸汽流量 , t / h ; 一再热蒸汽冷段 机侧 焓 , k J / 学兔兔 w w w .x u e t u t u .c o m 1 2 锅炉制造 总第 2 5 5期 k g ; 一再热蒸汽冷段 炉侧 焓, k J / k g ; 二 . 一 再热蒸汽热段 炉侧 焓, k J / k g ; 二 一再热蒸汽热 段 机侧 焓 , k J / k g ; D 一补水流量 , t / h ; 一补 充水焓 , k J / k g ; D m 一给水流量 , t / h ; 一给水焓 , k J / k g ; D 一锅炉侧给水流量 , t / , h; .6 一锅炉侧 给水焓 , k J / k g ; D 一厂用辅助蒸汽流量 , t / h ; 。 一 厂用辅助蒸汽焓 , k J / k g ; 一厂用辅助蒸汽返 回 水系数 , %; 一厂用辅助蒸汽返 回水焓 , k J / k g ; D 一主蒸汽泄漏 流量 , t / h ; D p ~主 给水泄 漏 流 量 , t / h 。 依托 导则 , 火电厂管道热效率不确定度测 量模型如式 4 、 5 所示。 正平衡管道热效率 7 p 1 0 0 % 厂 , P 6 , D 6 , T o , P o , D 0 , , p 二 , ⋯ , ,p 4 Y b 反平衡管道热效率 An 叼 g d 1一 1 0 0 % , P 6 , D6 , , P 0 , D 0 ,⋯ , ‰ , p f w 5 Y b △Q Q 1 Q 2十Q 3Q 4十Q 5Q 6 6 式中 7 7 鲥一管 道 热效 率 , %;Q 一锅 炉 热 负荷 入量 , 且充分考虑 了对测量结果 的不确定度产生 输 出热 量 Q 。 ,k J;Q 。 一汽 轮 机 热 耗 ,k J; 影响的分量。由此可见 , 当采用焓差法时 , 正 、 反 A Q 一机组管道热损失, k J; Q 一新蒸汽管道 平衡管道热效率不确定性的影响因素是一致的。 散热损失 , k J; Q 一带热量的工质泄漏热损失 , 焓值在管道热效率计算 中是重要参数 , 在实 k J; Q 一再热蒸汽 管道散热损失 , k J;Q 一给 际测量中无法直接在测量仪表显示 , 拟合过热蒸 水管道热损失 , k J; Q 一厂用辅助系统热损失 , 汽与给水 、 补水 比焓与温度、 压力的关系是进行管 k J; Q 一锅炉连续排污热损失, .j} t, 。 道热效率不确定度评定分析的基础, 各焓值拟合 所建立的模型涵盖了影响测量结果的全部输 表达式如表 1所示 。 表 1 管道热效率不确定评定所需比焓拟合公式 注 ①过热蒸汽 比焓表达式 中, 温度 T的单位为 K, 给水 比焓 、 补水 比焓表达式 中温度 t 的单位为℃ ; ②式中压力 P的单位 为 MP a ; ③式 中 拟合公式 的误差是与基于 I A P WSI F 9 7和 I F C 6 7标准 的计算值对 比。 1 将式 4 与表 1过热蒸汽、 给水和补水 比 率不确定度分析模型。 焓表达式结合得到可测量参数的管道正平衡热效 一 1 0 0 % 7 式中 D 。, D 6 , D ,h , T o , T b , ‰ , ⋯, P 。 , P 6 , r , O ; 合成标准不确定度 为 为由测量仪表直接测量。 式 7中 各 参 数 分 别 测 量 , 相 互 独 立 , 善 摹 砜 学兔兔 w w w .x u e t u t u .c o m 第 1期 曾大海 , 等 不确定度原理在火电厂管道热效率测算中的应用 。l 3 [ O D o 『 [ 『 [ ] 中 D 。, D 6 , D r , ,M , ‰ ,⋯, P 。, , 响 分别为可测参数 , D , D , , , o o o [ “ 。 】 『 0 p b u ] [ “ 蜘 ] 8 , ⋯ , P o, P , 包含 A、 B不确定度 的标准 不确定度 ] , 即 “ 。 。 . . . 厢 , 9 式 7 各测参数求偏导得到 D 。, D , D , ⋯ , 的灵敏度 。 D 0, 巩 , ⋯ , 。 ‰ , 如 1 0 最终, 正平衡管道热效率测量不确定度为 1 叼p 2u 叼 p。 2 同理, 将式 5 与表 1 过热蒸汽、 给水和补 水 比焓表达式结合得到可测量参数的管道反平衡 热效率不确定度分析模 型, 求得测量模 型下 的各 参数灵敏度和合成标准不确定度 “ 叼 。 最终 , 反平衡 管道热效 率测量 不确定 度为 叩 1 7 2u 叩 。 2 . 3 管道热效率测量不确定度评定流程 依据测量不确定度原理 , 管道热效率测量不 确定度评定过程如图 2 。 图2 测量不确定度评定过程 3 管道热效率测量不确定度案例分析 本文选取某 电厂亚临界压力一次 中间再热 , 3 5 0 MW 煤 、 煤气混烧 的机组 , 匹配 日本三菱重工 公 司设计制造 的两缸 T C 4 F 2 8型汽轮 机, 分析测 量不确定度在管道热效率 的应用。在计算值 中不 确定度 /Z 包含了A类、 B类不确定度。分析 前提是不考虑锅炉连续排污热损失以及给水管道 散热损失 。 3 . 1 数据 采集 进行管道热效率测量 的不确定度评定 , 对测 试数据进行处理是基础 , 进行 不确定度分析计算 是核心。 3 . 1 . 1 测量仪 器及精度 . 火电厂管道热效率测量的不确定度评定涉及 的测量参数主要有发电机功率 、 压力 、 温度和流量 四大类 , 测量仪器及精度见表 2所示 。 3 . 1 . 2 D C S在线数据采集 火电厂 D C S 分散 控制系统 具有强大 的数 据采集功能 。在不具备完整 、 规范化的管道热效 率测试条件 的前提下 , 运用 D C S在线采集数据进 行不确定的评定具有一定的现实意义。 火 电厂数据采集点繁多 , 实时更新快 , 机组负 荷及相关测点参数变化迅速 , 因此 , 对机组 的工况 进行稳定性判定是基础 , 选取稳定工况段数据进 表 2测量仪器及精度 学兔兔 w w w .x u e t u t u .c o m 1 4 锅炉制造 总第 2 5 5期 行工程定量分析才有可比性。基于 导则 进行 对应的稳定参数分别如表 3 和图3 所示。 机组工况判稳, 判定参数及要求, 归一化稳定工况 2 O 1. 5 1 0 O. 5 赫 0 0 一 0 5 1. O 一 1 5 2.O 0 1 0 0 2 0 0 3 0 0 40 0 5 00 6 00 7 0 0 8 0 0 9 0 0 数 搠 个数 图 3 稳定工况稳态数据段 注 ①图中所示 为稳态工况下负荷 、 主蒸汽压力 、 一 级抽汽压力数值归一化后所成图形 ②温度 变化范 围在 5 ℃ 内, 稳定工况下数值归一化后波动较大 , 属于正常状况。 选取机组夏季 2 3 0 02 4 0 0期 间, 3 4 5 MW 负荷稳定工况段的测量参数 , 采用焓差法进行管 道正、 反平衡热效率计算及其不确定度 的分析分 别如表 4 、 5 。 表 3 机组稳定工况判定参数及要求 川 注 A 、 A ⋯ 、 A 分别为每类参数对应的参数最大值、 参数最小值、 参数平均值 表 4 管道热效率计算结果 注 以上计算值为稳定工况下稳定 负荷段对应 的多组测量值的平均值 计算所得 。 计算得到管道正 、 反平衡热效率为9 7 . 5 3 6 % 。 其 中主蒸汽管道散热损失为 0 . 0 8 5 % , 再热蒸 汽 管道散热损失为 0 . 1 4 % , 带热量工质泄漏损失为 0 . 0 5 9 % , 厂用辅助蒸汽热量损失为 2 . 1 8 %。 学兔兔 w w w .x u e t u t u .c o m 第 1期 曾大海 , 等 不确定度原理在火电厂管道热效率测算中的应用 1 5 计算表明, 管道正平衡热效率为 9 7 . 5 3 6 4 - 0 . 0 6 1 %; 不确定度影响最大的为温度测量值, 主蒸汽温度 机侧 、 主蒸汽温度 炉侧 、 再热冷 段蒸汽温度 炉侧 、 再热热段蒸汽温度 炉侧 ; 其次是压力值的测量, 主蒸汽压力 机侧 、 主蒸 汽压力 炉侧 、 再热冷段蒸汽压力 炉侧 ; 流量 的测量对管道效率的影响普遍偏小, 可忽略不计。 管道反平衡热效率为 9 7 . 5 3 6 4- 0 . 0 6 1 2 %; 不确 定度影响因素与管道正平衡热效率不确定度影响 因素相 同。 管道正、 反平衡热效率进行测量不确定度的 评定结果亦几乎一致, 综合精确度, 最终管道正、 反平衡热效率为 9 7 . 5 3 60 . 0 6 1 %。该工况下 锅炉热效率为 9 2 %, 汽轮发电机组热效率 4 5 %, 厂用 电 率 5 % , 则 全 厂 热 效 率 为 4 0 . 3 8 0 0 . 0 2 5 % , 全 厂 热 耗 为 8 9 1 5 . 3 4-5 .6k J / k W h ,供 电煤耗 率为 3 2 0 . 6 4 4-0 . 2 0 g e e / k W h , 按机组年平均运行 4 0 0 0小时 , 标准煤 5 0 0元每 吨计算 , 全年将 产生 2 8万元 的费用 偏 差 。 若 管 道 效 率 取 9 9 %, 全 厂 热 效 率 为 4 0 . 9 8 6 %, 全厂热耗为 8 7 8 3 . 5 k J / k W h , 供 电煤耗率为 3 1 5 . 9 g c e / k W h 。对 比管道热效 率取值 9 9 %和实际计算值 9 7 . 5 3 6 4-0 . 0 6 1 % , 全厂热效率、 全厂热耗和供电煤耗率分别偏差 [ 0 . 5 8 1 , 0 . 6 3 1 ] %、 [ 1 2 6 . 3 , 1 3 7 . 4 ] k J / k W h 和 [ 4 . 5 4 , 4 . 9 4] g c e / k W h 。 4 结 论 1 依托 D L / T 6 0 6 . 3 2 0 1 4 , 采用常规的焓差 法建立 了火 电厂管道正、 反平衡热效率 的测量不 确定度分析模型, 拟合了适用于亚临界火 电机组 的给水焓值, 以及补水焓值表达式; 给出了管道热 效率不确定度评定 的完整流程和计算方法 , 分析 管道热效率不确定度 的影响因素。 2 计算分析表明, 在不考虑给水管道散热损 失 、 锅炉连续排污热损失的条件下 , 管道正 、 反平 衡热效率为 9 7 . 5 3 6 4- 0 . 0 6 1 %, 全厂供电煤耗 率为 3 2 0 . 6 4 0 . 2 0 g c e / k W h , 不考虑不确 定度将产生 2 8万元每年的费用偏差 ; 与管道效率 学兔兔 w w w .x u e t u t u .c o m 1 6 锅炉制造 总第 2 5 5期 取 9 9 %相 比, 全厂供 电煤耗偏差达 [ 4 . 5 4, 4 . 9 4] g c e / k W h , 引起 评估 偏 差达 3 1 7 . 8万 元 至 3 4 5 . 8万元每年。 3 不确定度分析表明, 对管道热效率不确定 度影响最大的因素是温度测量值 , 其次是压力值 的测量 , 降低其不确定度应增加温度值的测量 准 确度 ; 同时 , 增设高压缸排汽温度 、 压力 , 中压缸入 口温度 、 压力 的测量点 , 加强火电机组工质泄漏点 的相关测量 , 可更全面 的评价火 电厂各管道热损 失 , 进而全面评价管道效率及其不确定度。 4 开展管道热效率不确定度 的评定对研究 建立统一、 规范化的火力发电厂管道热效率的试 验准则 , 加强火力发电厂热力管道安装 、 检修对管 道热力性能的评价和验收具有一定 的工程实际意 义 ; 当管道散热采用传热学等不同测算方法时 , 其 正 、 反平衡管道热效率不确定度 的影响 因素是不 同的。 参考文献 郑体宽. 热力发电厂[ M] . 北京 中国电力出版社 , 2 o o 8 . 王文欢, 刘祺福等. 不同负荷下燃煤电厂管道效率 变化特性 的研 究 [ J ] . 华东 电力 , 2 0 0 9 , 3 7 6 1 O 451 0 48. 石奇光, 仝宁等. 发电厂管道热效率及行业标准的 探讨[ J ] . 华东电力, 2 0 0 4, 3 2 1 4 4 4 7 . 石奇光, 薛玉兰, 马庆. 关于发电厂管道热效率的 反平衡算法及其分析 [ J ] . 华东工业大学学报, 1 9 9 7 , 1 9 3 9 91 0 5 . 王志国, 马一太, 卢苇. 不确定度分析原理在锅炉 热效率测算 中的应用 [ J ] . 中国电机工程学报, [ 6 ] [ 7 ] [ 8 ] [ 9 ] 2 0 0 5 , 2 5 3 1 2 51 2 9 . 郭振宇, 吴正勇等. 汽轮机性能试验不确定的评定 的研究[ J ] . 动力工程学报, 2 0 1 0 , 3 0 1 0 7 4 9 7 5 4. 阎维平 , 杜海玲. 锅炉热效率测试的不确定度评定 [ J ] . 锅炉技术, 2 0 1 0 , 4 l 5 1 01 7 . Ahme t Eg e,Ha el Me h me t Sa h i n. De t e r mi na t i o n o f un c e r t a i nt i e s i n e n e r g y a n d e x e r g y a na l y s i s o f a p o we r p l a n t [ J ] . E n e r g y C o n v e r s i o n a n d Ma n a g e m e n t , 2 0 1 4,8 5 3 9 94 0 6 . C l a a s W a g n e r ,Ki m.H. E s b e n s e n . A s y s t e ma t i c a p pr o a c h t o a s s e s s i n g me a s u r e me nt u n c e r t a i n t y f o r C02 e mi s s i o n s f r o m c o a l fi r e d p o w e r p l a n t s Mi s s i n g c o n t ri b u t i o n s f r o m t h e T h e o r y of S a m p l i n g T O S [ J ] . C h e m i c a l E n g i n e e r i n g R e s e a r c h a n d D e s i g n, 2 01 1. 8 9 1 5 721 5 8 6. 吕当振, 段学农等. 各因素测量不确定度对锅炉热 效率测试质量的影响[ J ] . 热力发电, 2 0 1 4 , 4 3 5 5 45 8. G u i d e t o t h e E x p r e s s i o n o f Un c e r t a i n t y i n me a s u r e me n t G U M [ S ] . 曾翔君 , 骆一萍等. 测量不确定度导论[ M] . 西安 西安交通大学出版社 , 2 0 1 1 . J J F 1 0 5 9 . 1 2 0 1 2 , 测量不确定度评定与表示[ S ] . D E / T 6 0 6 . 32 0 1 4, 火力发电厂能量平衡导则 第 三部分 热平衡[ s ] . 方超. 高精度的过热水蒸汽焓函数简单实用数 学模型的探讨 [ J ] .中国电机工程学报, 1 9 8 9 , 9 1 3 7 4 7 . 张子慧. 热工测量与自动控制[ M] .北京 中国建 筑工业出版社, 2 0 0 7 . 藏迎诃 迎趋 藏 迎 、 广 告 储 感 1● 1j 1 J 1J 1 J m U rL r;r L rL rL 1 J] J 1J 1 j、 j l 2 3 4 5 学兔兔 w w w .x u e t u t u .c o m
展开阅读全文