PVA凝胶内基质扩散系数测定及其传质性能的优化

管清坤, 杨宏. PVA凝胶内基质扩散系数测定及其传质性能的优化[J]. 环境工程学报, 2017, 11(3): 1375-1382. doi: 10.12030/j.cjee.201512042
引用本文: 管清坤, 杨宏. PVA凝胶内基质扩散系数测定及其传质性能的优化[J]. 环境工程学报, 2017, 11(3): 1375-1382. doi: 10.12030/j.cjee.201512042
GUAN Qingkun, YANG Hong. Determination of diffusion coefficient of PVA gel and optimization of its mass transfer performance[J]. Chinese Journal of Environmental Engineering, 2017, 11(3): 1375-1382. doi: 10.12030/j.cjee.201512042
Citation: GUAN Qingkun, YANG Hong. Determination of diffusion coefficient of PVA gel and optimization of its mass transfer performance[J]. Chinese Journal of Environmental Engineering, 2017, 11(3): 1375-1382. doi: 10.12030/j.cjee.201512042

PVA凝胶内基质扩散系数测定及其传质性能的优化

  • 基金项目:

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

  • 中图分类号: X703

Determination of diffusion coefficient of PVA gel and optimization of its mass transfer performance

  • Fund Project:
  • 摘要: 凝胶的传质性能是影响其处理效能的主要因素之一,为了提高PVA凝胶内基质的扩散能力,对PVA凝胶的传质性能进行了研究,通过单因素实验确定添加剂活性炭和CaCO3对传质性能的影响程度,采用2因素5水平的响应面分析法(response surface methodology,RSM)优化了PVA凝胶的传质性能。研究结果表明:添加适量的活性炭和CaCO3可有效改善PVA凝胶的传质性能;在添加剂最优添加量下,NH4+-N、NO2--N和NO3--N在凝胶中的有效扩散系数分别为1.363 7×10-9、1.085 0×10-9和1.019 9×10-9 m2·s-1。
  • 污泥未经处理随意排放堆置,会造成严重的环境污染问题。国际上污泥主要有土地利用、卫生填埋、焚烧和投海等4种处置方式[1]。其中,填埋处置对技术指标要求相对宽松、运行成本低,是现阶段我国污泥处置的主要方式,且为简易的单独填埋,即污泥经过脱水消化后,直接倾倒于事先设置好的填埋坑中,并采用膜或土覆盖进行封场。由于我国污水处理厂对污泥处理的重视度不高,技术资金投入力度也不够,导致污泥的含水率高、物理力学性质差,不仅达不到市政污泥的填埋标准,而且造成填埋场库容的日益紧张,更严重的是会埋下安全隐患[2],如深圳下坪垃圾填埋场和山西太原垃圾填埋场均发生过填埋体的滑坡事故。为此,在《城镇污水处理厂污泥处理处置技术指南》[3]的国家规范中对填埋污泥的各项指标做出了明确规定。与此同时,我国的污泥产量也在逐年增加,目前,国内上海老港、成都长安、深圳下坪、杭州天子岭等填埋库区库容已经出现严重不足。因此,污泥填埋场内坑体加固与库内污泥深度脱水减量成为目前多数填埋场所面临的问题。

    现阶段常用机械压滤方式对污泥进行深度脱水。从机械脱水原理来看,机械压滤的过程实质上就是污泥的排水固结过程,即在总应力作用下孔隙水不断被排出的过程。孙政等[4]对污水处理厂脱水污泥的固结特性进行了研究,发现污泥的固结规律与一般黏土差别较大,超孔隙水压力的消散较慢。朱婧等[5]对污泥、淤泥、粘土的压缩特性进行了对比研究,认为污泥与淤泥的固结不同,在外力荷载下其固结过程可以分多个阶段。王鹏等[6]采用纤维加筋技术,研究了不同掺量下加筋污泥的固结压缩特性。范惜辉等[7]选用普通硅酸盐水泥和硫铝酸盐水泥作为固化材料,研究了固化污泥在不同应力下的压缩、渗透规律。机械压滤技术一般是先采用化学药剂预调质,使污泥颗粒的结合水释放出来之后,再其进行深度脱水,将湿基含水率降至60%以下。采用药剂真空预压法处理污泥也是如此,调质改性后的污泥与工程废浆类似,在真空预压过程中存在流固的两相转变,并在大部分时间里处于弹塑性状态,此时需要采用土力学中的固结理论进行分析[8]。武亚军课题组[9][10]对于无机药剂调质过的新鲜污泥的真空固结特性进行了研究,由于暂存库区污泥与新鲜污泥性质不同,固结特性也必然有差异,而目前关于这方面的研究并未见有所报道。此外,FeCl3是比较常用的一种调质药剂,而芬顿试剂在污水处理中应用较多,但不常用于污泥调质,因此,一方面为了对新鲜污泥与暂存库污泥进行对比,另一方面为了对FeCl3和芬顿试剂的调质效果进行对比,本研究采用土力学中的固结实验对分别采用2种不同药剂调质过的填埋污泥的压缩固结特性进行了研究,研究结果可为机械压滤和真空预压处理填埋污泥的工程实践提供参考。

    实验选用的药剂分别为FeCl3·6H2O、FeSO4、浓硫酸,以上药剂均为分析纯(AR)。实验所需H2O2通过40%的双氧水颗粒(昌乐鑫富强商贸有限公司)按浓度比例添加。实验仪器主要包括中压固结仪和电子天平等。

    对暂存库区填埋污泥与新鲜污泥的各项物理指标进行了测试,其中比重采用比重瓶法测试;密度采用环刀法测试;含水率采用低温烘干法测试;有机物采用灼失量法测试。结果表明,填埋污泥与新鲜污泥的含水率分别为74.1%和82.17%,有机质含量分别为40.9%和64.9%,比重分别为1.87和1.57,密度分别为1.2 g·cm−3和1.02 g·cm−3。由此可见,填埋污泥具有比新鲜污泥含水率低、有机物含量低、比重和密度大等特点。

    固结实验的药剂调质方案中氯化铁的添加量分别为0%、10%、20%、30%和40%;芬顿试剂的添加方案如表1所示。装入烧杯中置于常温下放置24 h,待污泥与药剂充分反应后,再均匀装填入固结仪,每个实验组别设置2组平行实验。由于污泥含水率较高,初级固结应力较大时容易发生冒浆,选取初级固结应力为3.125 kPa,加荷比为1,将最大固结应力增加至400 kPa。根据《土工试验方法标准》(GB/T50123-1999),加载过程中按规定时间记录百分表读数,由于污泥稳定达到稳定标准时间较长,每级加载48 h。第1级固结应力p设置为3.125 kPa,之后按6.25、12.5、25、50、100 kPa依次加载,以沉降量小于0.005 mm·h−1为沉降稳定的标准。

    表 1  污泥固结实验芬顿试剂调质方案
    Table 1.  Consolidation test plan of sludge conditioned by Fenton reagent
    编号Fe2+/%H2O2/%H2O2/Fe2+
    1441
    2461.5
    3482
    4881.5
    58122
    68163
      注:添加量表示占污泥干基的质量比。
     | Show Table
    DownLoad: CSV

    添加药剂之后污泥的孔隙比e (指污泥中孔隙体积与固体体积的比值,初始孔隙比e0采用含水率和比重进行换算,压缩过程中的孔隙比根据压缩量测试)会发生较大的变化,不同种类的药剂添加量与初始孔隙比的关系如图1所示。由图1可知,经过药剂调质改性后,e均有不同程度的增大。采用FeCl3调质后(图1(a)),污泥的初始孔隙比e0变化明显,从原始污泥的4.098上升至6.681,但随着药剂掺量的增加,污泥的孔隙变化较为平缓,最终达到7.244。采用芬顿改性后(图1(b)),当Fe2+的掺量为4%时,污泥孔隙比随着H2O2掺量的增加变化明显,由4.802上升至7.092;当Fe2+的掺量为8%时,污泥孔隙比随着H2O2掺量的增加变化较为缓慢,最终达到4.908。这是由于在药剂调质过程中产生了大量气体,这些气体不能完全从污泥中排出,而是积存分布在污泥内部,导致污泥的空隙变多,从而使得孔隙比增大。

    图 1  不同药剂种类的添加量与初始孔隙比的关系
    Figure 1.  Relationship between sludge initial porosity and the dosage of different agents

    压缩实验每级荷载的加载周期为48 h,不同FeCl3添加量下改性污泥的孔隙比e与荷载p的关系如图2(a)所示。污泥初始孔隙比为4.098,略大于常规的软黏土,经过药剂调质后,污泥的初始孔隙比随着药剂添加量的增加逐渐变大,当药剂添加量为40%时,孔隙比达到7.244。不同芬顿配比掺量下改性污泥的e-p关系如图2(b)所示。由图2(b)可知,对比2种药剂调质后的污泥发现,在初级荷载作用下,样品的孔隙比迅速减小。通过对固结应力为100 kPa时的孔隙比变化量进行了分析,发现调质污泥的压缩量基本均达到总压缩量的70%以上。这是因为在前期压缩过程中,调质污泥较原始污泥颗粒间的空隙总量更多,颗粒间没有形成骨架,强度较低,在较低应力作用下,孔隙水排出顺畅,压缩量大,孔隙比减小幅度大。经过3.125、6.25、12.5、25、50 kPa荷载作用下,芬顿改性污泥的沉降量较大,孔隙被大幅压缩;当荷载大于50 kPa时,污泥沉降速率逐渐减慢,沉降幅度逐渐减小,污泥孔隙比被压缩幅度也逐渐减小。由于原始污泥中有机质含量较高,存在大量具有一定承载力的微生物残体和胶结絮状有机物,通过添加FeCl3与芬顿试剂可以一定程度上破坏微生物残体和胞外聚合物,减少了有机物的含量,样品更容易发生固结压缩。

    图 2  调质污泥的孔隙比-压力曲线
    Figure 2.  e-p curve of conditioned sludge

    将调质污泥的孔隙比e与固结应力p之间的关系可以绘制成半对数坐标曲线 (e-lgp),如图3所示。由图3(a)可知,孔隙比e与固结压力lgp之间呈明显的线性关系,这一结果与常规淤泥类似。填埋污泥的压缩指数为0.64,调质污泥的压缩指数在0.776~0.795,跟新鲜脱水污泥差别较大,且与常规淤泥在数值上也较为接近[1, 5]。由图3(a)可知,污泥初始孔隙比的拟合值要略大于实验实测值。这是由于污泥的机械脱水和长期填埋类似于加卸载过程,压缩之后产生不仅存在塑形变形,而且也会发生一定程度的回弹。污泥的实际孔隙比和理论孔隙比的差值在一定程度上反映了不可恢复的塑形变形。同时,重塑制样及拟合精度也会对该结果产生一定影响。由于污泥中含有凝胶状结构,颗粒接触点处有一定的胶结力,能承受一定的压力而变形较小,使得在初期加荷阶段曲线平缓。此外,一般的原状土由于前期固结应力的存在会发生自重应力下的固结。其压缩曲线会出现屈服应力的折点,污泥的e-lgp曲线近似为一条直线,由此可知,调质污泥不存在应力屈服点,属于欠固结土。不同芬顿配比掺量下改性污泥的e-lgp曲线如图3(b)所示。污泥孔隙比随固结应力增大基本呈线性减小,压缩指数Cc为0.444~0.591,整体上小于原始污泥和经FeCl3调质后的污泥,和常规淤泥土较为接近,但仍属于高压缩性土。

    图 3  调质污泥的孔隙比-压力曲线
    Figure 3.  e-lgp curve of conditioned sludge

    固结系数Cv是表示孔隙水压力消散快慢的物理量,固结系数越大,固结速度越快,反之越慢。采用时间平方根法可得到调质污泥固结系数Cv与固结应力p之间的关系。图4(a)为采用FeCl3在各级压力下的固结系数变化结果。由图4(a)可知,在初级压力下,调质污泥的固结系数在10−3 cm2·s−1数量级变化,随着固结应力的增大,污泥的固结系数逐渐减小。此外,随着FeCl3掺量的增大,固结系数也越大,且在前几级固结应力下固结系数的减小幅度也越来越明显。由各条固结系数曲线关系可以说明在每一级固结应力下,随着FeCl3添加量的增加,污泥的固结系数增大,即FeCl3掺量越多,固结过程中孔隙水压力消散越快,这一点与新鲜脱水污泥固结系数的变化规律一致[9]

    图 4  调质污泥固结系数与固结应力的关系
    Figure 4.  Relation between consolidation coefficient and consolidation stress of conditioned sludge

    对比芬顿调质的实验结果(图4(b))可知:当Fe2+的添加量为4%时,样品的固结系数随着H2O2添加量的增加而增大,当H2O2的掺量为8%时达到最大;当Fe2+的添加量为8%时,样品的固结系数随着H2O2掺量的增大呈现先增大后减小的趋势。这是由于当H2O2添加量过多时,不仅不能分解产生更多的羟基自由基,反而会使最初产生的羟基自由基发生泯灭[11]。就初级固结应力下的固结系数而言,芬顿试剂改性后初级固结应力下Cvmax=9.88×10−3 cm2·s−1,当固结应力增大到400 kPa时,Cv=1.85×10−3 cm2·s−1;经过40%的FeCl3调质后Cvmax=2.91×10−3 cm2·s−1,随着固结应力的增大,Cv减小至4.98×10−4 cm2·s−1。因此,当Fe2+添加量为4%、H2O2掺量为8%时,在固结应力作用下污泥的孔压消散最快。

    污泥与淤泥、黏土最大的区别是污泥的固体物质中存在40%~60%的有机物,这些有机物大多数是生物处理过程中的微生物残体[12]。因此,污泥中的水分赋存状态非常特殊,除了具有孔隙水、表面结合水以外,存在絮凝体内部的结合水和细胞颗粒内部的细胞水(或称为生物水)[13]。这些水赋存于可以承载一定压力的有机物絮体中,这使污泥中水分难以快速排出,因此,孔隙水压力消散时间非常漫长[14]。添加药剂在一定程度上使得微生物残体胞内水以及有机絮体中的结合水释放,从而大大缩短了固结时间。

    图5所示,通过固结系数可以推演出污泥在各级固结压力下的渗透系数k。由图5(a)可知:k和固结应力的规律与固结系数Cv和固结应力的规律相似,受固结应力影响较大;在0~25 kPa阶段,污泥的渗透系数下降明显,渗透性变差,这是因为大孔隙被压缩成小孔隙或密闭孔隙,孔隙比迅速减小导致排水困难。经过试剂调质后,长期填埋污泥的渗透系数增大,初级固结应力下的k从10−7 cm·s−1数量级增大到10−6 cm·s−1数量级,随着压力的增大,k减小为10−8 cm·s−1数量级;当FeCl3的掺量为40%时,样品在初级固结应力下的k=4.439×10−6 cm·s−1,随着固结应力的增加,k减小至3.796×10−8 cm·s−1;采用芬顿试剂调质的污泥在初级固结应力下的k=6.48×10−6 cm·s−1 (图5(b)),随着固结应力的增大k则下降至9.94×10−8 cm·s−1。若以固结系数和渗透系数作为污泥固结效果好坏的指标,芬顿试剂的效果更佳。

    图 5  调质污泥渗透系数和固结应力的关系
    Figure 5.  Relation between permeability coefficient and consolidation stress of conditioned sludge

    此外,土体渗透性与其孔隙比密切相关。有研究[5]表明,土体孔隙比e与lgk存在一定的关系。图6为在不同FeCl3掺量下调质污泥渗透系数与孔隙比关系曲线。由图6可知,随着孔隙比的减小,渗透系数也逐渐减小,反之,渗透性增大,e与lgk之间的线性关系近似成立。在一定孔隙比范围内,相同孔隙比下10%添加量的污泥渗透系数一直小于同样孔隙比的其他掺量污泥,20%、30%和40%添加量下的污泥在孔隙比为5~7时渗透系数较为接近,但是随着孔隙比减小,实验组污泥的渗透系数出现差异,且随着药剂掺量减小而递减,添加量为20%的实验组渗透系数接近于10%添加量的实验组。

    图 6  经FeCl3调质污泥渗透系数与孔隙比关系
    Figure 6.  Relation between permeability coefficient and porosity of sludge conditioned by FeCl3

    对比调质污泥和原始污泥可以发现,在重合的孔隙比区间内,相同孔隙比下对照组的渗透系数要大于添加药剂的实验组。这是因为与天然细粒土一样,污泥由于初期的加药絮凝和板框压滤,其初始状态的结构也很复杂,一旦扰动,原有的过水通道的形状、大小及其分布都会改变,故渗透系数也不同。这一点与普通的性质相近,相同孔隙比时扰动土样的渗透系数通常小于原状土样[15]。实际加药时由于搅拌分散以及药剂的作用导致污泥颗粒分散变小,絮状结构一定程度上被破坏,使得调质后的污泥在相同孔隙比下的k小于原始污泥。

    1)污泥经试剂调质后能在较短时间内排水固结稳定;调质污泥在低荷载水平下沉降量较大,在高荷载水平下沉降逐渐平稳,孔隙变化不大;经过FeCl3调质后的污泥压缩性增大,压缩指数由0.64增大至0.776~0.795。

    2)在初级固结应力下,调质污泥的固结系数在10−3 cm2·s−1数量级内变化。添加FeCl3的实验组Cv,max=2.91×10−3 cm2·s−1;芬顿调质实验组Cv,max=9.88×10−3 cm2·s−1。比阻和固结系数并不是简单呈负相关性,两者之间的定量关系还需要进一步研究。

    3)渗透系数受固结应力影响较大。当FeCl3的掺量为40%时,样品在初级固结应力下的渗透系数为4.439×10−6 cm·s−1,在400 kPa下,渗透系数减小为3.796×10−8 cm·s−1;采用芬顿试剂调质的污泥在初级固结应力下,k=6.48×10−6 cm·s−1,在400 kPa下,k=9.94×10−8 cm·s−1

    4)在芬顿试剂最小添加量时(4% Fe2++4% H2O2)的调质效果均比FeCl3最大添加量40%时的调质效果要好,因此,建议在工程实践中采用芬顿试剂进行调质污泥。

  • [1] 李丽娟,马贵平,赵林果. 固定化酶载体研究进展[J]. 中国生物工程杂志,2015,35(11):105-113
    [2] 赵欣欣,孙玲,董玉玮,等. 固定化微生物技术及其在污水处理中的应用[J]. 水处理技术,2015,41(7):17-20
    [3] 张业健,叶海仁,郑向勇,等. 固定化包埋技术在水处理领域的应用进展[J]. 工业水处理,2011,31(1):9-12
    [4] ALI M,OSHIKI M,RATHNAYAKE L,et al. Rapid and successful start-up of anammox process by immobilizing the minimal quantity of biomass in PVA-SA gel beads[J]. Water Research,2015,79:147-157
    [5] 方定,谌婕妤,吴丁山,等. 微生物包埋固定化技术及其在废水处理中的应用[J]. 广东化工,2015,42(17):133-134
    [6] SHULER M L,KARGI F. Bioprocess Engineering[M]. New York:Prentice Hall,2002
    [7] JANJAROEN D,LING Fangqiong,MONROY G,et al. Roles of ionic strength and biofilm roughness on adhesion kinetics of Escherichia coli onto groundwater biofilm grown on PVC surfaces[J]. Water Research,2013,47(7):2531-2542
    [8] 茆云汉,王建龙. 聚乙烯醇固定化微生物新方法的研究[J]. 环境科学学报,2013,33(2):370-376
    [9] CHEN Kuocheng,LIN Yingfeng. Immobilization of microorganisms with phosphorylated polyvinyl alcohol (PVA)gel[J]. Enzyme and Microbial Technology,1994,16(1):79-83
    [10] SONG S H,CHOI S S,PARK K,et al. Novel hybrid immobilization of microorganisms and its applications to biological denitrification[J]. Enzyme and Microbial Technology,2005,37(6):567-573
    [11] SEKER D C,ZAIN N A M. Response surface optimization of glucose production from liquid pineapple waste using immobilized invertase in PVA-alginate-sulfate beads[J]. Separation and Purification Technology,2014,133:48-54
    [12] 黄真真,陈桂秋,曾光明,等. 固定化微生物技术及其处理废水机制的研究进展[J]. 环境污染与防治,2015,37(10):77-85
    [13] 王杏佳,徐丽,高静,等. PVA固定化微生物凝胶小球的制备及性能研究[J]. 塑料工业,2014,42(2):89-92
    [14] TAKEI T,IKEDA K,IJIMA H,et al. Fabrication of poly (vinyl alcohol)hydrogel beads crosslinked using sodium sulfate for microorganism immobilization[J]. Process Biochemistry,2011,46(2):566-571
    [15] KIMURA Y,ITOKAWA H,NOTO K,et al. Stability of autotrophic nitrogen removal system under four non-steady operations[J]. Bioresource Technology,2013,137:196-201
    [16] BAE H,YANG H,CHOI M,et al. Optimization of the mechanical strength of PVA/alginate gel beads and their effects on the ammonia-oxidizing activity[J]. Desalination and Water Treatment,2015,53(9):2412-2420
    [17] WU K Y A,WISECARVER K D. Cell immobilization using PVA crosslinked with boric acid[J]. Biotechnology and Bioengineering,1992,39(4):447-449
    [18] ARAVINTHAN V,TAKIZAWA S,FUJITA K,et al. Factors affecting nitrogen removal from domestic wastewater using immobilized bacteria[J]. Water Science and Technology,1998,38(1):193-201
    [19] 李峰,吕锡武,严伟. 聚乙烯醇作为固定化细胞包埋剂的研究[J]. 中国给水排水,2000,16(12):14-17
    [20] 钟成华,刘鹏,张文东,等. 卡拉胶、膨润土改性PVA包埋小球[J]. 环境工程学报,2013,7(8):2837-2843
    [21] 李志荣. 包埋颗粒内硝化菌增殖衰亡规律及硝化特性研究[D]. 上海:上海交通大学,2012
    [22] 吴亚杰,李辉军,戴昕,等. 无机材料对固定化微生物凝胶球性能的影响[J]. 环境工程学报,2014,8(8):3289-3293
    [23] 郑佩,陈芳艳,唐玉斌,等. 固定化菌藻微球的制备、表征及其对富营养化湖水的修复[J]. 环境工程学报,2014,8(5):1999-2005
    [24] 晋凯迪,于鲁冀,陈慧敏,等. 亚硝化细菌吸附-包埋固定化方法的优化及其脱氮性能[J]. 水处理技术,2015,41(7):57-60
    [25] REDDY L V A,WEE Y J,YUN J S,et al. Optimization of alkaline protease production by batch culture of Bacillus sp. RKY3 through Plackett-Burman and response surface methodological approaches[J]. Bioresource Technology,2008,99(7):2242-2249
  • 加载中
计量
  • 文章访问数:  3504
  • HTML全文浏览数:  3011
  • PDF下载数:  407
  • 施引文献:  0
出版历程
  • 收稿日期:  2016-01-29
  • 刊出日期:  2017-03-10
管清坤, 杨宏. PVA凝胶内基质扩散系数测定及其传质性能的优化[J]. 环境工程学报, 2017, 11(3): 1375-1382. doi: 10.12030/j.cjee.201512042
引用本文: 管清坤, 杨宏. PVA凝胶内基质扩散系数测定及其传质性能的优化[J]. 环境工程学报, 2017, 11(3): 1375-1382. doi: 10.12030/j.cjee.201512042
GUAN Qingkun, YANG Hong. Determination of diffusion coefficient of PVA gel and optimization of its mass transfer performance[J]. Chinese Journal of Environmental Engineering, 2017, 11(3): 1375-1382. doi: 10.12030/j.cjee.201512042
Citation: GUAN Qingkun, YANG Hong. Determination of diffusion coefficient of PVA gel and optimization of its mass transfer performance[J]. Chinese Journal of Environmental Engineering, 2017, 11(3): 1375-1382. doi: 10.12030/j.cjee.201512042

PVA凝胶内基质扩散系数测定及其传质性能的优化

  • 1. 北京工业大学建筑工程学院, 北京市水质科学与水环境恢复重点实验室, 北京 100124
基金项目:

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

摘要: 凝胶的传质性能是影响其处理效能的主要因素之一,为了提高PVA凝胶内基质的扩散能力,对PVA凝胶的传质性能进行了研究,通过单因素实验确定添加剂活性炭和CaCO3对传质性能的影响程度,采用2因素5水平的响应面分析法(response surface methodology,RSM)优化了PVA凝胶的传质性能。研究结果表明:添加适量的活性炭和CaCO3可有效改善PVA凝胶的传质性能;在添加剂最优添加量下,NH4+-N、NO2--N和NO3--N在凝胶中的有效扩散系数分别为1.363 7×10-9、1.085 0×10-9和1.019 9×10-9 m2·s-1。

English Abstract

参考文献 (25)

返回顶部

目录

/

返回文章
返回