你的位置:水利新闻 >> 资讯 >> 水利科技资讯 >> 详细内容 在线投稿

响应面法优化光催化/臭氧氧化深度处理炼油废水

排行榜 收藏 打印 发给朋友 举报 来源: 环境工程学报 发布者:cjk3d
热度80票 浏览37次 时间:2015年4月01日 19:34

b{^6? x0作者:倪可 王利平 肖春宝 李新颖水利新闻1Em:i4?"j

E fEA.a-p5m o+Z0摘要:利用响应面方法(RSM)对光催化/臭氧氧化深度处理炼油废水工艺进行优化,考察了臭氧通量、光催化剂投加量、初始pH 和反应时间对于处理效果的影响,提出采用该工艺的数学模型及优化后的工艺参数。结果表明,各影响因子对COD 去除率影响顺序为反应时间>光催化剂投加量>初始pH >臭氧通量,方程的F 值为11.54,相关系数为0.9537,调整相关系数为0.915,说明数学模型可以较好地模拟真实的反应曲面。优化得到最佳的工艺参数: 臭氧通量1.05 L/min、光催化剂投加量0.33 g/L、初始pH 7.51、反应时间96.95 min,在该条件下,对COD 去除率为97.88%,与预测值99.49% 接近。采用95%处理水和5%新鲜水混合,水质达到了循环冷凝水的补充水水质指标要求。水利新闻0v;ap#r`g$~-C7w)T

水利新闻M zQ1@8Ci F

关键词:炼油废水 响应面法 工艺优化 光催化 臭氧氧化

^2|[mk;|K k&S3[0

3pWIF7]dZ ]#I0中图分类号X703 文献标识码A 文章编号1673-9108(2014)12-5349-07水利新闻TFd)k[:Z

{+Z l*{ Y-r2Qu;F0Optimization of photocatalytic/ozonation treatment process for oil-refinery wastewater with response surface methodology

F'[+vkK0qc0 水利新闻N4w0F9rM4y

Ni Ke1 Wang Liping1 Xiao Chunbao2 Li Xinying1

&a ps y Y`D;]Y T0

2G }^^1I0Abstract:Response surface methodology was employed to optimize the photocatalytic/ozonation treatment process for oil-refinery wastewater.The effects of O3 flux,catalyst dosing quantity,initial pH and reaction time on COD removal were investigated.What’s more,the mathematical model and optimized parameters of process were also achieved.The results showed that the effects of different factors on COD removal were as follows: reaction time > catalyst dosing quantity > initial pH > reaction time.The mathematical model could greatly response the practical process with the equation F-value of 11.54,correlation coefficient of 0.9537 and adjusted correlation coefficient of 0.915.The optimum process parameters of oil-refinery wastewater treatment were O3 flux of 1.05 L/min,catalyst dosing quantity of 0.33 g/L,initial pH of 7.51,reaction time of 96.95 min.Under the optimum conditions,the actual COD removal rate was 97.88% which closed to the predicted value of 99.49%.The mixed water with the ratio of treated water and fresh water of 95∶ 5 reached the circulation condensate makeup water quality standards.

}4v?`G0rXQ0 水利新闻}(oR `Pp7sh-Z

Key words:oil-refinery wastewater; response surface methodology; process optimization; phototcatalysis; ozonation水利新闻j.WX2n^

水利新闻(os ^ E1B*nt1ecl

炼油废水的深度处理是指对于石化企业二级污水处理后水,用生化法、离子交换法、超滤法和化学氧化法等进行再次处理,处理后水用于回用,如作为循环冷却水的补充水等,不但可以减轻对于环境的污染,也可以为炼油企业节约用水,增加经济效益[1]。高级氧化技术可有效氧化降解高浓度有机污染物,也常常作为废水深化处理的技术手段[2,3]。光催化是高级氧化技术的一种,处理废水的先进性已被公认,但是由于TiO2纳米颗粒局限性: TiO2粉末易凝聚,不易从溶液中分离,无法回收再利用[4]。在处理工业有机废水方面研究较少,如何能将光催化技术运用于产业化处理过程是如今的研究热点[5]。臭氧氧化法作为另一种较为成熟的高级氧化技术,近年来在废水处理领域运用较为广泛[6]。考虑将光催化和臭氧催化氧化技术进行联用,对于炼油废水进行深度处理,以期达到回用的目的。

n0]y!o/x$fs ap0 水利新闻:khJ#r"|'{op-r

响应面方法(RSM)以统计学为基础,用于设计实验,评估单个变量以及多个变量之间的相互作用,可在有限的实验次数下优化实验操作的参数[7-9]。同传统实验方法相比,RSM 不易造成时间和资源的浪费,还能对实验过程进行分析、结果进行预测[10]。RSM 已在国外制浆废水处理[11]、饮用水净化[12]和重金属处理[13]等领域得到了广泛的应用。RSM 作为一种先进的研究方法也逐渐引起了国内研究学者的重视,运用该方法优化深度处理废水回用尚属首次。

._^ QsixQ&V0 水利新闻w c GjRt.{H_

本研究采用RSM 的中心组合设计方法(CCD)进行实验设计,利用Design Expert 8.0.5 软件(Stat-Ease 公司,美国)对实验结果进行了分析,并建立了以COD(COD 去除率反映了体系氧化降解程度)为响应值的二次多项式。对深度处理工艺条件进行优化,预测了实验结果,并在优化条件下对炼油废水进行了深度处理实验,评价了水质,对用于循环冷凝水补充水具有参考价值。水利新闻$`ji#W%HN!u%`

k,\T6jh/v5R01 材料和方法水利新闻~1w+V j&UqZ0q

7]? cRO2fto1d01.1 实验用水水利新闻9WG+K7HTS

q:b%@k8qO_0炼油废水取自安庆石化综合污水处理厂二级处理后的外排水,水质结果如下: pH 为6.46,SS 为23mg/L,石油类为123.6 mg/L,COD 为174 mg/L,NH3-N 为41.9 mg/L,硫化物为12.9 mg/L。

'ZI JCz6O0 水利新闻X&Nz#dmr^C

参考同济大学[14]、刘守新等[15]的制备方法,并进行改进,采用溶胶-凝胶法在改性球形活性炭颗粒表面负载TiO2晶体,在N2保护下程序升温煅烧制得TiO2/AC 光催化剂。该种催化剂物理特征及其吸附性能平均指标如下: 颗粒大小为1.88 mm,装填密度为577 g/L,碘吸附值为1 092 mg/g,孔容积为0.86 cm3/g,比表面积为1 377m2/g,负载比例为21.49%,TiO2晶体粒径为11.2 nm。

IfMe/Fy0

5w [/A8iT01.2 实验装置

K;nVo lY&C0

9\'} q9XC x!E~/_7b0zO0本实验采用自制的光反应器及其配套装置,见图1。光源为紫外光光源,单一紫外波长为254 nm。紫外灯光(30 W)距离水面距离3.5 cm保持不变。水利新闻3a7vh6s\)L

水利新闻yiX-X"u-\!Nx

1.3 实验方法

yN)SA)B*V o,C0 水利新闻"L0rY?n hA

用c (NaOH)1.00 mol/L 或c (H2SO4)1.00mol/L 调节废水初始pH,投加光催化剂,开启紫外光反应装置,泵入臭氧,对炼油废水进行氧化处理,反应一定时间后,将处理后溶液进行离心分离取上清液,ρ(COD)采用重铬酸钾法测定,计算COD 的去除率。利用RSM 进行实验设计、参数优化和结果预测。水利新闻.B,n$J(Q~hyC)O @8h

水利新闻R[3\5kVl8}uk&B"d

1.4 实验设计

R*a%b])TF6O&L3~0 水利新闻!Hl1G4k5kv

实验选取臭氧通量、光催化剂投加量、初始pH和反应时间为实验因素,按照CCD 法设计了一个4因素5 水平的实验方案,各因素水平和编码见表1。水利新闻3wLr&lCR[

:a @[D J.|0以炼油废水COD 去除率(y)为响应值,用多项式回归分析对实验数据进行拟合,可以得到二次多项式模型[16,17],其模型为:水利新闻1nV kH y}R TF

水利新闻0{8M~ dQ)Y

式中: y 为响应值; xi 、xj为实验因素; β0为常系数;βi为线性系数; βii为二次项系数; βij为交互项系数;ε为随机误差。

kP#Q%UU+j-~0 水利新闻1\Om/@3fj M

二项式模型拟合质量的优劣是由决定系数(R2)所决定的,使用二项式模型和方差分析(ANOVA)对数据进行拟合和分析,以获得自变量和响应变量之间的关系式。水利新闻7mW2}~ Q

m s%J(z6|02 结果与分析

X$YK j0n]0 水利新闻%a&c d)h V$LZ!qz$R9@c%|

2.1 回归方程与数据分析水利新闻4~:S]+xv-lL/b

&xl u,?2U RW0一共进行30 组实验,每次实验都按照设计方案进行,然后测定ρ(COD),最后计算得到COD 去除率,实验设计方案和结果见表2。水利新闻]LB iW#s5G k&T%}C

1KZV@C o8TeY0实验序号如1、4、15、22、25 和26 的实验因素的条件是相同的,主要是为了验证实验Pure error(纯误差)[18]。利用Design Expert 软件对表2 的实验数据进行多元回归拟合,最后得到臭氧通量、催化剂投加量、初始pH 和反应时间与炼油废水COD 去除率之间的二次多项式回归方程:水利新闻^p:[ dsY ~

水利新闻v3l4\Ae~U

使用方差分析(ANOVA)来检查二次多项式模型对实验结果进行模拟的充分性和显著性,方差分析结果见表3。水利新闻DVvC3\f

h3?-D0i7L9{!m0表3 为方差分析结果。模型方程统计显著性是由F 值确定,F 值越大,则表明方程的显著性越强[19]。P 值小于0.0500,表明模型因素项具有显著性,P 值大于0.1000,认为模型因素项是非显著性的。由表3 可知,二次多项式模型的F 为11.54,远大于1,P < 0.0001,说明模型具有较好的回归效果和较强的显著性。ρ(O3)、ρ(催化剂)、初始pH 和反应时间的F 值分别为0.066、20.18、20.15 和55.51,因此各因素对处理效果影响的显著性顺序为反应时间>催化剂投加量>初始pH >臭氧通量。

"JN9s1_"\|]4S(J0

])Ea w g z;C0相关系数R2 和Radj2是检验模型可信度和准确性的重要指标,R2 和Radj2越靠近1,表明模型越能有效反映实验的数据,R2和Radj2越靠近0,表明模型越不能有效反映实验的数据。该回归方程的相关系数(R2)为0.9537,调整相关系数(Radj2)为0.9150,说明该回归方程能较好地模拟真实的曲面。

"R?9U@0xG5K!Eb&z0

M D}O#|0图2 ~图4 为残差图,残差是Design Expert 软件预测值与实际值之间的误差。内学生化残差用来表征标准偏差偏离实际、预测响应值的程度,在图形上表现为数据点是否呈现线性分布。外学生化残差是用于考虑各个响应值数据相对于拟合的回归方程是否为异常点,大部分学生化残差分布在± 3.5 范围之内。由图2 可见,实验数据点呈现了线性分布,而且表明了回归模型拟合得较好,该直线上的数据点不存在任何问题,进一步证实了预测值与实际值较接近。由图3、图4 可见,数据点随机分布,没有任何趋势,由于实验操作中的误差原因,图4 中的第3 个数据点超出范围,在后续的数据处理中将其舍去。

yaLYUfrK0

XK$@Uw _0图5 为废水COD 去除率的实测值与预测值关系图。由图可见,实验数据点基本分布在直线上,深度处理后炼油废水COD 去除率的实际值和预测值的皮尔逊相关系数(pearson correlation coefficient)为0.981,说明模型预测值与实测值之间存在良好的线性关系,模型具有较高的准确度。水利新闻mTDR d|;ctm@

水利新闻)i yo8^;Unw` a0p

2.2 因素效应分析

xMi,oy6S0

-LY9Aerv'Q"Zj5L0为了更好地考察臭氧通量、光催化剂投加量、初始pH 和反应时间交互作用对COD 去除效果的影响,绘制了反应三维曲面图和等高线图。

LSbo8em-s0 水利新闻cu^:FT Nja J

当初始pH 7,反应时间70 min,臭氧通量与光催化投加量的交互作用见图6。水利新闻k,OX A%Q9p

水利新闻 eI6u7eQV

当pH 为7,光催化剂投加量为0.3 g/L,臭氧通量与反应时间的交互作用见图7。

,O ^,Y6sm~._5K4u0 水利新闻#B!W7Le Yi

当臭氧通量为1.5 L/min,反应时间为70 min时,光催化剂投加量与pH 的交互作用见图8。水利新闻"U'kcd+T;V"L0c

水利新闻'TH@C!hK e

当臭氧通量为1.5 L/min,初始pH 为7 时,光催化剂投加量与反应时间的交互作用见图9。水利新闻+Y%x L-C$}by

+U ?y*I-Kf)q0由上述的三维反应曲面图及其对应的等高线图可知,反应过程中始终存在光催化和臭氧氧化这两个同时进行的反应单元,其中的TiO2/AC 催化剂既是光催化剂也是臭氧催化氧化的催化剂。当催化剂的投加量较少,反应时间较长时,光催化单元占主导地位,作为载体活性炭大量吸附溶液中的污染分子,创造出一个高浓度的污染物环境,通过传质平衡将污染分子转移到TiO2分子表面,污染物分子与TiO2光反应产生的大量的OH 进行反应,通入的O3因为催化剂投入量的不足,主要是在紫外光照射下,形成少量的OH 参与催化氧化过程,延长反应时间增加催化氧化的时间可以显著提高COD 的去除效率。当催化剂投加量较大,反应时间较短时,此时的臭氧催化氧化单元占主导地位,在很短的时间内具有大比表面积的载体将污染物分子与O3一同吸附至催化剂的表面,增加了O3与污染物质的接触效率,同时,由于高浓度的O3环境,使更多的O3分解成为氧化性能更强的OH,在很短的时间内完成催化氧化反应。如果一方面增加催化剂的投加量,一方面增加反应时间,COD 的去除率提高明显,最终可以达到将近100%。

LLy B'@Pk0

B&`/\ l e1j02.3 工艺优化

b!Q~ C B5Rx1B0 水利新闻Q ~!}mA gKd

优化目的是寻找最佳反应条件使得响应值最大化,即寻找处理效果最好的工艺条件。根据响应面二次多项式回归方程,使用Design Expert软件,求解COD 去除率的最大值,得到优化条件见表4。水利新闻NA%MD-^ CK]3^

XY U gkuxyi+d02.4 验证实验水利新闻V,D1b$a6{ ka{

!s j5nX N0采用优化方案,对炼油废水进行深度处理,进行了5 组平行实验,各个指标取平均值,并将实验结果与《安庆石化回用于循环冷却水补水的水质标准》进行了对比,实验结果见表5。

3h^_}8lUL7LA0 水利新闻 uv"I)c)Mw_

由表5 可知,对炼油废水进行深度处理后COD的去除率为97.88%,与预测结果99.49% 接近,说明RSM 对光催化/臭氧氧化深度处理炼油废水工艺过程优化结果是可行的。处理后水质除了腐蚀速率外均达到了《安庆石化回用于循环冷却水补水的水质标准》中指标的要求。对光催化/臭氧氧化深度处理后的95% 出水和5% 新鲜水混合,进行静态阻垢和旋转挂片实验,结果如下: 最大阻垢率为83.1%,最小阻垢率为70.8%,平均阻垢率为76.95%。最大腐蚀速率为0.07 mm/a,最小腐蚀速率为0.036 mm/a,平均腐蚀速率为0.053 mm/a。水质达到了《安庆石化回用于循环冷却水补水的水质标准》中腐蚀速率的要求。说明该水源用于石化回用循环冷却水补水是可行的,具有明显的经济效益和环境效益。水利新闻eE I G,_S

水利新闻CH/u8}(C5B3x

3 结论水利新闻g,v4nO*@8Mg S

水利新闻*j1p!I s:g$\,{2F$\z

(1)采用光催化/臭氧氧化深度处理炼油废水,利用Design Expert 软件建立了反应的数学模型,进行了方差分析,其中数学模型的F < 0.0001,R2 和R2adj分别为0.9537 和0.9150 模型是显著的,回归方程可以很好地模拟真实的反应曲面。建立的响应面模型影响因素及影响因素的显著性顺序为: 反应时间>催化剂投加量>初始pH >臭氧通量。

4PQ)UkU(]S0

:z8q+?ug(l0(2)采用响应面方法(RSM)对反应工艺进行优化,得到光催化/臭氧氧化深度处理炼油废水的最优条件: 臭氧通量1.05 L/min,光催化剂投加量0.33 g/L,初始pH 7.51 和反应时间96.95 min。水利新闻v J7l Yh5A \.E

水利新闻.A1?P7hU0d-hI

(3)绘制了因素间的三维曲面图和等高线图,存在同时进行的光催化和臭氧催化氧化单元,反应单元的主导地位由TiO2/AC 光催化剂的投加量决定。水利新闻I#u8Q!Np5u%zz4e

水利新闻;z pi"j;f!POCQ

(4)在最优条件下对炼油废水进行了深度处理,COD 的去除率为97.88%,与预测值接近,95%出水和5%新鲜水混合后,用于石化回用循环冷却水补水是可行的,具有经济和环境效益,同时为后续中试实验提供了参考。

nDf)} MX0 水利新闻P4_Dm3F][!d

参考文献

TE SA!nT1e0

Z0A;Xi7zf p2xk0[1]苟钰娴,杨凤林,石利军.一体式A/O 摇动床工艺处理石化废水.环境工程学报,2007,1(2):34-38水利新闻3k%}-p \4V6W9V"l0rL

m%PUbY:ob0Ge Y.X.,Yang F.L.,Shi L.J.Treatment of petrochemical wastewater using integrative A/O swimbed process.Chinese Journal of Environmental Engineering,2007,1(2):34-38(in Chinese)水利新闻5Zj'c4K)s%E[;tI

水利新闻*S|7Y!\_:pW

[2]Zhu X.B.,Tian J.P.,Liu R.,et al.Optimization of Fenton and electro-Fenton oxidation of biologically treated coking wastewater using response surface methodology.Separation and Purification Technology,2011,81 (3):444-450

j%S| O(d S}E]0

-Yv:K6O F2s |$t%W}J0[3]张冬梅,刘蕾,褚衍洋.电-Fenton 氧化降解2,4 -二氯酚.环境科学研究,2012,25(9):1041-1046

DJ'\2Z&tm4D0

'bD H Uik-L?0Zhang D.M.,Liu L.,Zhu Y.Y.Degradation of 2,4 -dichlorophenol by electro-Fenton oxidation process.Research of Environmental Sciences,2012,25 (9):1041-1046 (in Chinese)

0x0rRG-\_ ?o0

.F X$jU~K8X ?0[4]Arnan I.,Dona R.J.M.,Rendon E.T.,et al.TiO2 activation by using activated carbon as a support: Part Ⅱ.Applied Catalysis B: Environmental,2003,44(2):153-160水利新闻"C']pW5D.J:P

水利新闻"U(Id ^&s.dc

[5]钱伯章,朱建芳.炼油废水处理技术新进展.化工环保,2009,29(2):99-104

}Tzx9F&o a t*It0

+A`H;y k8d'D$k&r!a0Qian B.Z.,Zhu J.F.Development of technologies for petrochemical wastewater treatment.Environmental Protection of Chemical Industry,2009,29 (2):99-104 (in Chinese)

_c5| j+tIIrGf[0 水利新闻:N'V1n*I8h;Ces!U&n

[6]蔡少卿,戴启洲,王佳裕.非均相催化臭氧处理高浓度制药废水的研究.环境科学学报,2011,31(7):1440-1449

|5O$ZB.~p J0 水利新闻'g;FJ~[$w[Y

Cai S.Q.,Dai Q.Z.,Wang J.Y.Heterogeneous catalytic ozone oxidation for pharmaceutical wastewater treatment.Acta Scientiae Circumstantiae,2011,31 (7):1440-1449(in Chinese)水利新闻M#bK6u%l&M

"t~U ] W(}L0[7] Alaton I.A.,Tureli G.,Hanci T.O.Treatment of azodye production wastewaters using photo-Fenton-like advanced oxidation processes: Optimization by response surface methodology.Journal of Photochemistry and Photobiology A: Chemistry,2009,202(2-3):142-153水利新闻bs)[jbu"O8v

水利新闻Y/mH0iT2}:o mP

[8]Zhong K.,Wang Q.Optimization of ultrasonic extraction of polysaccharides from dried longan pulp using response surface methodology.Carbohydrate Polymers,2010,80 (1):19-25

)a"h9t3]F`;F0e,X `%J0 水利新闻%NYd#W;a Ig Q

[9] Jain M.,Garg V.K.,Kadirvelu K.Investigation of Cr (VI)adsorption onto chemically treated Helianthus annuus: Optimization using response surface methodology.Bioresource Technology,2011,102(2):600-605水利新闻y-v4D1i{!c"V [R@

水利新闻(mO'znMO'Q

[10]崔旭峰,杜尔登.响应面方法优化光催化降解抗生素类污染物过程研究.环境科学与管理,2010,35 (7):70-74,99水利新闻UN2EGT2Q

水利新闻 {6iVX-b}PU

Cui X.F.,Du E.D.Optimization of photocatalytic degradation of antibiotic pollutant by response surface method.Environmental Science and Management,2010,102(2):600-605(in Chinese)

#K3H$F e6J0 水利新闻*E)p!a:?@ez

[11]Wang J.P.,Chen Y.Z.,Wang Y.,et al.Optimization of the coagulation flocculation process for pulp mill wastewater treatment using a combination of uniform design and response surface methodology.Water Research,2011,45(17):5633-5640

!k%I1y}&v N.|0 水利新闻}(`9EIki A

[12]Hasa H.A.,Abdullah S.R.,Kamarudin S.K,et al.Response surface methodology for optimization of simultaneous COD,NH + 4 -N and Mn2 + removal from drinking water by biological aerated filter.Desalination,2011,275(1-3):50-61

%j.Yg+r@0pc-G\0 水利新闻J {vz7c\p'[s

[13]Bhatti M.S.,Reddy A.S.,Thukral A.K.Electrocoagulation removal of Cr(VI)from simulated wastewater using response surface methodology.Journal of Hazardous Materials,2009,172(2-3):839-846

-{s8[nJ~5N,r0 水利新闻 b|,zQ_TG.Mxw

[14]胡中华,王晓静.一种活性炭负载二氧化钛光催化剂的制备方法: 中国,CN101244383A,2008-08-20水利新闻9ud6|!t7WeV"l T_#x

5fjvQ mP'o7])@{(T0[15]刘守新,陈曦.TiO2/活性炭负载型光催化剂的溶胶-凝胶法合成及表征.催化学报,2008,29(1):19-24

)~kMzL0R y8^0 水利新闻kR&h-f \

Liu S.X.,Chen X.Sol-gel preparation and characterization of activated carbon supported TiO2 photocatalyst.Chinese Journal of Catalysis,2008,29(1):19-24(in Chinese)

+fr,rG5H"vF0 水利新闻a7nSoLB?&j

[16]Neseli S.,Yaldiz S.,Turkes E.Optimization of tool geometry parameters for turning operations based on the response surface methodology.Measurement,2011,44(3):580-587水利新闻6e9FV pY#VSG.w

\j,y.JWI0[17]Bezerra M.A.,Santelli R.E.,Oliverra E.P.,et al.Re-sponse surface methodology (RSM)as a tool for optimization in analytical chemistry.Talanta,2008,76 (5):965-977水利新闻JwK.p!Bd$M

)e~BG z7j E0[18] Zhang Z.M.,Zheng H.L.Optimization for decolorization of azo dye acid green 20 by ultrasound and H2O2 using response surface methodology.Journal of Hazardous Materials,2009,172(2-3):1388-1393

#k;U j8^-a.}(s1c0

0a!t${-F4g0[19]Singh P.,Shera S.S.,Banik J.,et al.Optimization ofcultural conditions using response surface methodology versus artificial neural network and modeling of L-glutaminase production by Bacillus cereus MTCC 1305.Bioresource Technology,2013,137(4):261-269水利新闻 gu Q,l-Ffo| XHU

水利新闻h;Ax TP7h5ll'e

作者简介: 倪可(1989-),男,硕士研究生,研究方向: 污水处理理论及技术。水利新闻&~/jNf8Y&gXi

顶:7 踩:10
对本篇资讯内容的质量打分:
当前平均分:-1.06 (34次打分)
【已经有29人表态】
10票
感动
3票
路过
4票
高兴
4票
难过
2票
搞笑
3票
愤怒
1票
无聊
2票
同情
下一篇:市政泵站雨天排放对黄浦江、苏州河水质的影响上一篇:多级跌水有机填料型人工湿地处理农村生活污水

广告投放区


水资讯网接受会员自助投稿


2015年河海大学举办EFDC_Explorer7.2培训

广告投放区

广告投放区