强电离放电脱除NO的化学反应动力学模拟

李萍, 依成武, 李文浩, 李珏, 王晶. 强电离放电脱除NO的化学反应动力学模拟[J]. 环境工程学报, 2017, 11(9): 5078-5083. doi: 10.12030/j.cjee.201609250
引用本文: 李萍, 依成武, 李文浩, 李珏, 王晶. 强电离放电脱除NO的化学反应动力学模拟[J]. 环境工程学报, 2017, 11(9): 5078-5083. doi: 10.12030/j.cjee.201609250
LI Ping, YI Chengwu, LI Wenhao, LI Jue, WANG Jing. Chemical reaction kinetics simulation of NO removal by strong ionization discharge[J]. Chinese Journal of Environmental Engineering, 2017, 11(9): 5078-5083. doi: 10.12030/j.cjee.201609250
Citation: LI Ping, YI Chengwu, LI Wenhao, LI Jue, WANG Jing. Chemical reaction kinetics simulation of NO removal by strong ionization discharge[J]. Chinese Journal of Environmental Engineering, 2017, 11(9): 5078-5083. doi: 10.12030/j.cjee.201609250

强电离放电脱除NO的化学反应动力学模拟

  • 基金项目:

    国家自然科学基金资助项目(51278229)

    江苏省六大人才峰会基金资助项目(JNHB-018)

  • 中图分类号: X511

Chemical reaction kinetics simulation of NO removal by strong ionization discharge

  • Fund Project:
  • 摘要: 强电离放电脱除NO是集化学、物理学和环境工程学为一体的交叉学科,其脱除过程本质上是等离子体中进行的一系列微观过程,难以通过实验进行直接研究。故通过建立强电离放电脱除模拟烟气中NO的化学模型并运用MATLAB对其进行数值模拟以探讨NO的去除机理。数值模拟结果表明,在800 μs内产生大量·O和·OH,而仅产生少量·N,由此可知NO的去除主要是通过与·O和·OH的氧化反应生成HNO3实现的,在脱除过程中NO先被氧化成NO2,再进一步被转化成HNO3。同时,通过实验验证了所建立模型的合理性。实验结果表明:随着NO初始浓度的增加,NO脱除效率逐渐下降;随着含氧量的增加,NO脱除效率显著增大。
  • 河流自净是一个涉及物理、化学和生物的复杂过程,是河流在一定空间内恢复其洁净状态的现象[1-2]。河流自净能力的恢复是城市生态环境建设和景观保护的重要环节,而目前城市内河流普遍采用“三面光”的梯形硬质化渠道,河水流速快,沉降性能低,改变了原有自然生态本底和水文特征,削弱了河流的自净能力。目前,河流只监测基本的水文参数和水质参数,同时,监测河流健康状况的方法对监测员的技术要求高且不能做到在线实时监测。因此,迫切需要一种在线监测河流水质参数和自净能力的方法。

    荧光和紫外光谱技术因其具有灵敏度高、用量少、测量简单、不消耗化学试剂等优点[3],近年来,被广泛应用于河流、湖泊、海洋等自然水体中污染物的监测[4-5]以及污水处理厂的过程控制[6-7]、工业废水中特定污染物的鉴别[8-9]。三维激发发射矩阵(3D-EEM)光谱,被称为“荧光指纹”,被广泛应用于检测废水、表征河流中溶解性有机物(DOM)[10]。紫外可见光谱分析中特定波长254 nm处的吸光度值(UV254)可作为总有机碳(TOC)和溶解性有机碳(DOC)的替代参数[11-12]。河流净化过程包括稀释、沉淀、曝气等多种化学与生物机制,可以采用数学模型进行评价[13],KARRASCH等[14]从浮游微生物的胞外酶角度分析得出,工业废水使微生物耐受性增强,降解能力提高,赵长森等[15]采用生物学指数与水生物指示环境结合的方法评价水样污染程度、生态系统稳定性与河流及水库的健康程度。

    河流水质与自净能力的传统测定方法及参数选取较为复杂,而从河流微生物的生理状态的角度分析河流的自净能力鲜有研究。本研究将人工净化与河流自净功能的协同作用发挥出来,以渭河流域西安段的河流及污水处理厂为考察对象,采用三维荧光光谱、紫外光谱联用呼吸图谱的方法,考察了不同性质的河流及污水处理厂各处理单元中微生物与有机物之间的作用关系,探讨了光谱法与呼吸图谱法联用表征河流状态及自净能力的可行性,以期得到河流水质和健康状况的综合评判方法。

    2018年4—6月,对陕西省西安市境内的河流进行采样,包括皂河(A1~A6)、太平河(A7)、灞河(A8)3条纳污河,泾河(N1)、渭河(N2~N4)2条天然水体。为了对比自然净化和人工净化的异同,按照工艺处理单元顺序,对污水处理厂WWTP1分别采样,W1~W7分别为进水(格栅后端)、曝气沉砂池、氧化沟泥水混合物、氧化沟沉淀30 min后的上清液、氧化沟沉淀2 h后的上清液终沉池出水和最终排水。其中,通过收集氧化沟不同沉淀时间的上清液,可模拟得到终沉池沉淀过程中的样品。

    在渭河的众多支流中,皂河、太平河和灞河属于“三面光”设计类型的渠道式城市纳污河,其中,皂河[16]具有最大的泄洪和纳污能力,全长35.8 km,集水面积283 km2,接纳西安市城区60%的生活污水、工业废水及3个污水处理厂的出水。太平河属于皂河的支流,接纳西安市西部的污废水。灞河[17]接纳西安市东部的污废水,在流入渭河前设有人工湿地。泾河是渭河最大的支流,接纳陕西北部的污废水。

    采用有机玻璃采样器于水下0.5 m处收集得到样品,将水引到无菌聚乙烯瓶中,然后通过冰袋运输至实验室,放入冰箱4 °C冷藏保存,分析前,将水样自然升温至25 °C,将收集的样品混合均匀,量取300 mL进行呼吸图谱的测定,另外100 mL水样通过0.45 µm滤膜过滤,以除去大尺寸的悬浮固体,用于光谱测定。所有样品的检测分析均在采样结束后2~3 d内完成。

    采用日立F-7000型荧光分光光度计进行三维荧光光谱检测。检测条件为:采用氙弧灯为激发光源,激发波长Ex=200~400 nm,发射波长Em=200~600 nm,狭缝宽度与扫描间隔均为5 nm,扫描速度为2 000 nm·min−1,响应时间为0.5 s,灵敏度为中等,光倍增管电压为700 V,采用超纯水(18.3 Ω)作为空白水样,以消除水的拉曼散射。采用752N紫外分光度计于波长254 nm处测量UV254

    呼吸图谱采用序批式呼吸计量法[18],于西安绿标水环境公司提供的BM300分析平台进行测定,分别获得现场呼吸速率OURS、内源呼吸速率OURe和总呼吸速率OURT

    根据污染物负荷及断面功能属性将全部采样断面分为3类:纳污河(A1~A8);自然水体(N1~N4);污水处理厂(W1~W7)。

    代表性断面的三维荧光图谱见图1。根据CHEN等[19]的三维荧光矩阵图五区划分法,识别出上述3类断面的5个特征峰(图1(b)),分别为类色氨酸T峰(Ex/Em=275 nm/340 nm)、类酪氨酸S峰(Ex/Em=225 nm/340 nm)、腐殖质C峰(Ex/Em=(310~320) nm /(380~410) nm)、富里酸A峰(Ex/Em=(240~260) nm/(380~400) nm),3类水样表现出峰位置及荧光强度的差异。纳污河自上游至下游各峰的最大荧光强度呈下降趋势,皂河源头及上游、污水厂进水有机物含量极高,且类蛋白峰在水样中占优势,类腐殖酸荧光强度相对较低,这是由于这些断面是河流接纳污水的源头,有机污染程度高,微生物含量高;而皂河下游、渭河、污水厂二沉池及出水中类腐殖酸占优势,这是因为这些断面经过污水厂的强化生物作用及河流自净作用后,有机物含量小,微生物繁殖速率慢,这与HENDERSON等[20]的描述一致,T峰反映的是不稳定易降解有机物,在废水中占主导地位,与废水微生物活性相关,与BOD之间的相关性较强,腐殖质C峰、富里酸A峰为难降解有机物,在天然水中占优势。

    图 1  代表性断面的三维荧光图谱
    Figure 1.  Three-dimensional fluorescence spectra of representative sections

    图2所示,以传统水质参数COD值作为参照,可以看到UV254FT与COD变化趋势基本一致,其中FT代表T峰的最大荧光强度。纳污河自上游至下游污染程度逐渐降低,自然水体污染程度低且稳定,污水处理厂水样的COD和UV254值从进水至出水在氧化沟工艺阶段出现极大值,而FT在进水出现极大值,FT与BOD呈显著正相关性,这表示可生物降解的有机物,FT可用于监测污水处理厂工艺处理过程中有机物的去除效果。

    图 2  3类水样的不同有机污染指标对比
    Figure 2.  Comparison of organic pollution indicators of three types of water bodies

    表1所示,FTFC分别代表T峰和C峰的最大荧光强度,在皂河源头A1流入渭河N4的过程中,FTFC逐渐减小,对应的污染物的去除率分别为69.0%、49.2%,而污水处理厂从进水W1至出水W7过程中,T峰与C峰对应的污染物去除率分别为68.0%、33.0%,河流与污水厂的T峰去除率基本相同,河流中微生物去除难降解有机物的能力高于污水处理厂,这说明河流中微生物群落与污水处理厂有所不同,且河流微生物更容易降解难降解有机物。

    表 1  3类水体的不同断面的水质、光谱及呼吸图谱特征参数
    Table 1.  Characteristic parameters of water quality, spectra and respirogram of different sections of three types of water bodies
    水样类型采样断面编号COD/(mg·L−1)光谱特征参数呼吸图谱特征参数/(mg·(L·h)−1)
    UV254FTFCHIXBIXOURSOUReOURenOURT
    纳污河A1120.200.257 8933 7100.490.946.212.041.5910.27
    A289.500.196 1222 9670.580.992.430.870.724.53
    A373.000.184 5472 8840.610.972.990.910.746.50
    A445.000.163 7352 7320.680.941.460.621.873.05
    A554.000.163 7562 7150.671.000.650.320.291.50
    A639.000.143 2322 2850.671.041.840.880.845.54
    A732.000.132 8852 4160.681.000.880.500.513.07
    A825.000.112 4741 6070.640.971.700.310.474.81
    自然水体N118.50 0.061 2677850.630.931.200.490.445.10
    N223.00 0.072 1531 2970.501.050.510.880.814.23
    N322.50 0.112 7221 8090.641.080.390.580.602.25
    N426.50 0.082 4481 8830.551.031.030.700.765.82
    污水处理厂W1132.000.259 9994 0470.460.954.851.491.416.26
    W2127.00 0.258 8684 0500.460.964.181.211.078.56
    W3161.00 0.384 8854 0060.721.0121.013.5825.8424.43
    W4145.00 0.33 9733 2590.680.963.971.3914.2016.81
    W5123.00 0.253 6544 6690.730.950.770.870.982.59
    W632.00 0.143 5742 7820.681.032.480.670.648.22
    W722.500.123 2002 7120.711.001.910.860.888.81
     | Show Table
    DownLoad: CSV

    皂河源头A1的腐殖化指数HIXb[21]较小(表1),这说明DOM较不稳定,易于生物降解;河流断面A4~A8、N1、N3的HIXb要高于其他断面,同时这些断面的FT较小,其中A4、A5、A7、N3均为接纳污水处理厂排放水的河流断面,可见污水厂排放水中DOM腐殖化程度更高,DOM更稳定,易生物降解的物质较少,A6、A8、N1断面的水质较好,难降解物质占优势;而污水厂进入生化处理阶段后的断面A3~A7的HIXb高于河流,且没有明显的下降趋势,这说明污水厂的微生物对去除难降解物质的能力匮乏,从DOM的腐殖化程度的角度分析证实了上述结论。

    本研究的全部断面的自生源指数BIX[22]为0.93~1.08,差异较小,DOM具有较强的自生源特征,是生物细菌活动产生的。

    总呼吸速率OURT反映的是基质不受限制条件下微生物最大的呼吸速率,其值越高,说明微生物的降解有机物的潜能越高。微生物的现场呼吸速率OURS反映的是采样时水样中微生物的现场活性,OURs越高,说明微生物数量及活性较高。纳污河的OURs普遍比自然水体高(表1),说明纳污河中微生物数量大且活性较高,其接纳了大量人类活动产生的生活污水,含大量类蛋白、脂肪等有机物,致使微生物大量繁殖,而纳污河断面中A5和A7由于接纳了污水厂出水,有机物浓度被稀释,微生物现场活性较弱;污水厂的W3点为氧化沟泥水混合物,其OURS极高,这与污水厂生物处理阶段活性污泥含量高的结果一致,此阶段微生物大量繁殖,降解有机物的速率极高,与自然净化的慢速过程形成对比。

    通过OUR=OURT′、FT= FT′,将OURT-FT(图3(a))平面划分为4个区域(T1~T4);通过OUR= OURS′、FT= FT′,将OURS-FT(图3(b))平面划分为4个区域(S1~S4)。其中,OURT′、OURS′、FT′的定量根据河流接纳的水样性质的不同而有所不同。

    图 3  河流断面的FT与OURS、OURT联用表征河流污染程度与自净能力
    Figure 3.  FT of river sections combined with OURS and OURT for characterizing river pollution degree and self-purification ability

    T1表示OURT> OURT′、FT<FT′,为超量潜在自净能力区域,包含断面A6、A8、N1、N2、N4,这说明断面污染程度较低,但微生物的潜在降解能力较高,自净能力较高,此时河流进入自净过程的后续阶段。同时,这种生物活性的改善可能归因于外部环境条件的改善,如更好的供氧,因为生物活性比需要降解的污染物更充分,这表明水生生态系统的健康。T2表示OURT> OURT′、FT>FT′,为受污染区域,包含断面A1~A3,这说明河流受到严重有机污染,虽然微生物活性极强且大量繁殖,但有机物处于超饱和状态,可能超过了微生物的降解能力与河流自净容量,一段时间后,会形成黑臭水样。T3表示OURT< OURT′、FT<FT′,为人工净化完成区域,包含断面A4、A5、A7、N3,这些断面接纳污水厂出水,污染物被稀释,虽然水样表观上恢复了原本干净的状态,但仍含有较多复杂不易降解的有机物,须汇入河流,经微生物长期降解才能恢复水样原本健康的状态。T4表示OURT< OURT′、FT>FT′,为无法判定区域,本研究无断面出现此情况,污染程度大而微生物呼吸速率小的情况出现的概率非常小,也不符合自净理论。

    S1表示OURS> OURS′、FT<FT′,为无法判定区域,本研究中无断面出现在此区域,同时,有机物含量小而微生物现场耗氧速率较强的情况出现的概率非常小,也不符合自净理论;S2表示OURS> OURS′、FT>FT′,为受污染区域,包含断面A1~A3,与上述T2区域描述基本一致;S3表示OURS< OURS′、FT<FT′,为受基质限制的区域,包含除A1~A3以外的其他断面,说明微生物现场活性受基质限制而无法生存;S4表示OURS< OURS′、FT>FT′,为无法判定区域,说明有机物含量大而微生物现场耗氧速率较弱的情况出现的概率也非常小,同样不符合自净理论。

    综上,结合2种不同的区域划分方法,可以判定河流的状态和自净能力:A1~A3为受污染断面,污染程度可能超过了自净容量;A4、A5、A7和N3为人工净化完成的断面,仍有大量难降解有机物须经过水样自净完成净化过程;A6、A8、N1、N2和N4为进行到水样自净过程的后续阶段的断面,生物活性受基质含量限制,但微生物具有超量潜在自净能力,可能归因于外部环境的改善,如溶解氧的升高或更适宜微生物生存的温度。

    1) 3类水样(纳污河、自然水体与污水处理厂各处理单元)表现出不同的荧光强度与荧光特征值,纳污河及污水处理厂前处理单元普遍类蛋白最大荧光强度高于自然水体,自然水体及污水处理厂后处理单元的DOM的腐殖化程度高于纳污河,类蛋白最大荧光强度可用于区分水体类型及判定有机污染程度,腐殖化程度可用于判定DOM的难易降解程度。

    2)类蛋白T峰(Ex/Em=275 nm/340 nm)的最大荧光强度FT可作为反映污水处理过程中有机污染程度及微生物量的指标,纳污河自上游至下游、污水厂自进水至出水的FT逐渐减小,河流微生物群落与污水处理厂的活性污泥有所不同,纳污河微生物去除难降解有机物(C峰)的能力高于活性污泥。

    3) OURS通常用于表征微生物现场活性,OURT用于表征微生物降解有机物的潜能,自然水体的有机污染程度及微生物现场活性均较低,纳污河与污水处理厂生物处理单元的微生物现场活性与潜能颇高,微生物呼吸速率与水体有机污染程度密切相关。

    4)采用呼吸图谱与紫外光谱、三维荧光光谱联用,以OURS-FT,OURT-FT这2种区域划分方式为依据,研究了河流微生物与有机污染之间的作用关系,建立了判定河流的污染状态和自净能力的定量指标,为城镇两极分化条件下的河流生态管理提供参考。

  • [1] GAFFNEY J S, MARLEY N A. The impacts of combustion emissions on air quality and climate from coal to biofuels and beyond[J]. Atmospheric Environment, 2009, 43(1):23-36
    [2] 赵静, 严金英, 邱婧伟, 等. NaClO2/NaClO复合吸收液同时脱硫脱硝[J]. 环境工程学报, 2012, 6(10):3684-3688
    [3] MOK Y S, NAM C M, CHO M H, et al. Decomposition of volatile organic compounds and nitric oxide by non-thermal plasma discharge processes[J]. IEEE Transactions on Plasma Science, 2002, 30(1):408-416
    [4] CHANG J S, URASHIMA K, TONG Y X, et al. Simultaneous removal of NOx and SO2 from coal boiler flue gases by DC corona discharge ammonia radical shower systems:Pilot plant tests[J]. Journal of Electrostatics, 2003, 57(3):313-323
    [5] KIM Y, KANG W S, PARK J M, et al. Experimental and numerical analysis of streamers in pulsed corona and dielectric barrier discharge[J]. IEEE Transactions on Plasma Science, 2004, 32(1):18-24
    [6] DONG L M, LAN S, YANG J X, et al. Plasma chemical reaction for nitric oxide and sulfur dioxide removal in corona discharge reactor[C]. Electrical Insulation and Dielectric Phenomena, 2003 Report. Conference on LEEE. 2003:440-443
    [7] RAJANIKANTH B S, SUSHMA B R. Injection of N-radicals into diesel engine exhaust treated by plasma for improved NOx removal:A feasibility study[J]. Plasma Science & Technology, 2006, 8(2):202-206
    [8] MORENO S H, PACHECO M P, PACHECO-SOTELO J, et al. Modeling and experimental study on nitric oxide treatment using dielectric barrier discharge[J]. IEEE Transactions on Plasma Science, 2007, 35(5):1533-1540
    [9] 张文豪, 杜传进, 王伟. 低温等离子体技术净化柴油机尾气NOx的化学反应动力学模拟研究[J]. 拖拉机与农用运输车, 2008, 35(1):79-81
    [10] 孙永明, 夏文虎, 张勤. 低温等离子体技术净化船舶柴油机尾气的化学反应动力学模拟[J]. 上海海事大学学报, 2015, 36(2):79-83
    [11] SATHIAMOORTHY G, KALYANA S, FINNEY W C,et al. Chemical reaction kinetics and reactor modeling of NOx removal in a pulsed streamer corona discharge reactor[J]. Industrial & Engineering Chemistry Research, 1999, 38(5):1844-1855
    [12] WANG Z, ZHOU J, ZHU Y, et al. Simultaneous removal of NOx, SO2 and Hg in nitrogen flow in a narrow reactor by ozone injection:Experimental results[J]. Fuel Processing Technology, 2007, 88(8):817-823
    [13] MÄTZING H. Chemical kinetics of flue gas cleaning by irradiation with electrons[J]. Advances in Chemical Physics, 2007,80:315-402
    [14] TOKUNAGA O, SUZUKI N. Radiation chemical reactions in NOx and SO2, removals from flue gas[J]. Radiation Physics & Chemistry, 1984, 24(1):145-165
    [15] HERRON J T, GREEN D S. Chemical kinetics database and predictive schemes for non-thermal humid air plasma chemistry. Part Ⅱ. Neutral species reactions[J]. Plasma Chemistry & Plasma Processing, 2001, 21(3):459-481
    [16] 李文绚, 金保侠. 气体动力学计算方法[M]. 北京:机械工业出版社, 1990:14-32
  • 期刊类型引用(5)

    1. 马瑞,陈莎,樊守彬,焦玉方,王东洋. 北京市空气质量特征与典型区域污染源贡献分析. 环境科学学报. 2024(03): 317-326 . 百度学术
    2. 周雨玲,马得玺. 工业区大气污染源排放与环境健康风险评估. 皮革制作与环保科技. 2024(12): 82-84 . 百度学术
    3. 赵丽娟,黄梓朦,徐志钮,张硕,支九英. 基于石墨烯涂层高灵敏SPR-PCF大气污染浓度检测传感器及特性研究. 半导体光电. 2024(06): 870-878 . 百度学术
    4. 蔡芸,韩惠敏. 均等化试点政策与超大城市公共卫生系统韧性评价. 中国卫生事业管理. 2023(10): 730-735+750 . 百度学术
    5. 姜宠. 基于虚拟现实的智慧图书馆绿化环境优化设计系统. 信息与电脑(理论版). 2023(23): 91-93 . 百度学术

    其他类型引用(4)

  • 加载中
    Created with Highcharts 5.0.7访问量Chart context menu近一年内文章摘要浏览量、全文浏览量、PDF下载量统计信息摘要浏览量全文浏览量PDF下载量2024-052024-062024-072024-082024-092024-102024-112024-122025-012025-022025-032025-0400.250.50.7511.25Highcharts.com
    Created with Highcharts 5.0.7Chart context menu访问类别分布DOWNLOAD: 18.4 %DOWNLOAD: 18.4 %FULLTEXT: 63.5 %FULLTEXT: 63.5 %META: 18.1 %META: 18.1 %DOWNLOADFULLTEXTMETAHighcharts.com
    Created with Highcharts 5.0.7Chart context menu访问地区分布其他: 99.1 %其他: 99.1 %XX: 0.7 %XX: 0.7 %天津: 0.2 %天津: 0.2 %其他XX天津Highcharts.com
计量
  • 文章访问数:  1850
  • HTML全文浏览数:  1448
  • PDF下载数:  493
  • 施引文献:  9
出版历程
  • 收稿日期:  2016-12-19
  • 刊出日期:  2017-08-26
李萍, 依成武, 李文浩, 李珏, 王晶. 强电离放电脱除NO的化学反应动力学模拟[J]. 环境工程学报, 2017, 11(9): 5078-5083. doi: 10.12030/j.cjee.201609250
引用本文: 李萍, 依成武, 李文浩, 李珏, 王晶. 强电离放电脱除NO的化学反应动力学模拟[J]. 环境工程学报, 2017, 11(9): 5078-5083. doi: 10.12030/j.cjee.201609250
LI Ping, YI Chengwu, LI Wenhao, LI Jue, WANG Jing. Chemical reaction kinetics simulation of NO removal by strong ionization discharge[J]. Chinese Journal of Environmental Engineering, 2017, 11(9): 5078-5083. doi: 10.12030/j.cjee.201609250
Citation: LI Ping, YI Chengwu, LI Wenhao, LI Jue, WANG Jing. Chemical reaction kinetics simulation of NO removal by strong ionization discharge[J]. Chinese Journal of Environmental Engineering, 2017, 11(9): 5078-5083. doi: 10.12030/j.cjee.201609250

强电离放电脱除NO的化学反应动力学模拟

  • 1. 江苏大学环境与安全工程学院, 镇江 212013
基金项目:

国家自然科学基金资助项目(51278229)

江苏省六大人才峰会基金资助项目(JNHB-018)

摘要: 强电离放电脱除NO是集化学、物理学和环境工程学为一体的交叉学科,其脱除过程本质上是等离子体中进行的一系列微观过程,难以通过实验进行直接研究。故通过建立强电离放电脱除模拟烟气中NO的化学模型并运用MATLAB对其进行数值模拟以探讨NO的去除机理。数值模拟结果表明,在800 μs内产生大量·O和·OH,而仅产生少量·N,由此可知NO的去除主要是通过与·O和·OH的氧化反应生成HNO3实现的,在脱除过程中NO先被氧化成NO2,再进一步被转化成HNO3。同时,通过实验验证了所建立模型的合理性。实验结果表明:随着NO初始浓度的增加,NO脱除效率逐渐下降;随着含氧量的增加,NO脱除效率显著增大。

English Abstract

参考文献 (16)

返回顶部

目录

/

返回文章
返回