人工湿地填料磷去除效果的影响因素分析.pdf

农业环境科学学报 2 0 0 8，2 7(2)：7 4 8 7 5 2 J o u r n a l o f Ag r o E n v i r o n me n t S c i e n c e 人工湿地填料磷去除效果的影响因素分析 汤显强，李金 中2，刘学功 2，黄岁棵(1 南开大学环境科学与工程学院，天津 3 0 0 0 7 1；2 天津市水利科学研究所，天津 3 0 0 0 6 1)摘要：为分析人工湿地填料磷去除效果的影响因素，选取页岩、粗砾石、铁矿石、麦饭石等不同填料进行等温磷吸附实验和室内 小试磷去除研究。结果表明，填料类型和粒径影响填料磷吸附性能。上述填料磷吸附性能依次为页岩 铁矿石麦饭石 粗砾石且 页岩最大磷吸附量随粒径增加而减小。水力负荷(H L A)变化对不同填料磷去除效率影响效果不同，当 HL A从 2 0 0 mm d 降至 1 4 2 8 m i l 1 d-m 时 粗砾石构建湿地总磷去除率增加大约 1 0；麦饭石因为吸附渐进饱和总磷去除率呈现出不断下降趋势；铁矿石可 能受微生物影响导致 F e 等金属元素释放，总磷去除率从大约 6(H L A=2 0 0 m m d )上升至 3 5(H L A=1 4 2 8 mm d。)左右；页岩受 H L A变化影响较小，总磷去除率基本稳定在 8 0 以上。综合上述实验结果发现填料种类、粒径及水力负荷均影响湿地填料磷去除 效果。关键词：湿地；填料；磷去除；水力负荷 中图分类号：X 7 0 3 1 文献标识码：A 文章编号：1 6 7 2 2 0 4 3(2 0 0 8)0 2 0 7 4 8 0 5 I n fl u e n c i n g F a c t o r s o f Ph o s p h o r u s Re mo v a l P e r f o r ma n c e o f Co n s t r u c t e d W e t l a n d s S u b s t r a t e s T ANG Xi a n q i a n g ，L I J i n-z h o n g：，L I U X u e-g o n g：，HUANG S u i l i a n g (1 C o U e g e o f E n v i ron m e n t a l S c i e n c e a n d E n g i n e e ri n g，N a n k a i U n i v e r s i t y，T i a n j i n 3 0 0 0 7 1，C h i n a；2 T i a n j i n H y d r a u l i c S c i e n t i f i c R e s e a r c h I n s t i t u t e，T i a n j i n 3 0 0 0 6 1，C h i n a)Ab s t r a c t：I n o r d e r t o a n aly z e t h e f a c t o r s t h a t i n fl u e n c e d p h o s p h o r u s r e mo v a l p e r f o r ma n c e o f c o n s t r u c t e d w e t l a n d s s u b s t r a t e s，s h ale，c o u r s e g r a v e l，i r o n s t o n e a n d h o rnb l e n d e w e r e u s e d t o c o n d u c t t h e p h o s p h o r u s i s o t h e rm al a d s o r p t i o n e x p e ri me n t an d i n d o o r p i l o t s c ale p h o s p h o r u s r e mo v al r e s e a r c h P h o s p h o r u s i s o t h e rm a l a d s o rpt i o n e x p e ri me n t a l r e s u l t s i n d i c a t e d t h a t s u b s t r a t e t y p e s a n d gra i n s i z e s c o u l d b o t h i n fl u e n c e p h o s p h o r u s r e mo v a l c a p a c i t i e s T h e r e wa s a n o r d e r i n r a n k f o r p h o s p h o r u s a d s o rp t i o n c a p a c i t i e s o f the f o u r s u b s t r a t e s me n t i o n e d a b o v e：s h ale i ron s t o n e h o rnb l e n d e c o u r s e g r a v e l a n d t h e ma x i mu m p h o s p h o r u s a d s o rpt i o n c a p a c i t y o f s h a l e i n c r e a s e d a s the gra i n s i z e o f the ma t e r i al d e c r e a s e d T P D mo v a l e ffic i e n c y o f i n d o o r-p i l o t s c ale c o n s t r u c t e d w e tl a n d s wi t h d i f f e r e n t s u b s a t e s v a ri e d wi t h h y d r a u l i c l o a d i n g r a t e f HLA 1 Wh e n HL A wa s d e c l i n e d f r o m 2 0 0 mm d t o 1 4 2 8 mm d ，T P r e mo v al r a t e o f c o u pe g r a v e l b e d c o n s t r u c t e d we tl a n d i n c r e a s e d a b o u t 1 0 wh i l e h o r n b l e n d e b e d a p p e a r e d d e p r e s s i v e t I le n d s i n c e the s u b s t r a t e a c h i e v e d s a t u r a t e d a d s o r p t i o n g r a d u all y I r o n s t o n e tha t wa s m a y b e a ff e c t e d b y t h e a c t i o n o f mi c roo r g a n i s m s l e d t o t h e i n c r e a s e o f T P r e m o v a l r a t e s fr o m a b o u t 6(H L A=2 0 0 m m d )t o a p p r o x i ma t e l y 3 5(1 4 2 8 m m d-1)，o w i n g t o t h e r e l e a s e o f e l e m e n t F e S h ale b e d w a s h a r d l y affe c t e d b y t h e c h a n g e o f H L A，a n d the T P r e m o v a l r a t e m a i n t a i n e d a b o v e 8 0 d u ri n g the e x p e rime n t al p e ri o d T h e s e r e s u l t s i mp l i e d t h a t s u b s t r a t e t y p e s，g r a i n s i z e a n d h y d r a u l i c l o a d i n g r a t e we r e the f a c t o r s tha t i nfl u e n c e d p h o s p h o r u s r e mo v a l p e rfo r ma n c e o f c o n s t ruc t e d we t l a n d s s u b s t r a t e s K e y w o r d s：c o n s t ruc t e d w e tl a n d；s u b s t r a t e；p h o s p h o rus r e mo v al；h y d r a u l i c l o a d i n g r a t e(H L A)收稿 日期：2 0 0 7 0 3-2 3 基金项目：天津市科委重大基金项 目“天津市外环河水质改善及水资源保护工程”；国家 自 然科学基金(5 0 4 7 9 0 3 4)；教育部留学归国服务人员启动 基金；南开大学“引进人才科研启动项 目”；现代水利科技创新项 目X D S 2 0 0 7 0 5 作者简介：汤显强(1 9 8 1 一)，男，博士研究生，主要研究方 向为水污染控制。E ma i l：a n b i t i o n m a i l n a n k a i e d u c a 通讯联 系人：黄岁棵E m a i l：s l h u a n g n a n k ai e d u a n 维普资讯 http:/ 第 2 7卷第 2 期 农 业 环 境 科 学 学 报 7 4 9 近年来，人工湿地废水生态处理技术在 国内外被 广泛应用于处理城市生活污水_ 1 1，农业 面源污染嘲、富 营养化水体I3、4 1 等也具备较好的磷去除效果。关于人工 湿地磷去除机理及其国内外应用现状的总结详见文 献 5 。人工 湿地除磷 主要靠 湿地填 料对 磷 的吸 附沉 降 5 9 1。分析人工湿地填料磷去除性能的影响因素对 选择具有较高磷去除性能的填料和提高人工湿地除 磷效率意义重大。人工湿地填料磷去除性能研究大多 集中于除磷填料的选择，即通过等温磷吸附实验 比较 不同填料 的兰格缪尔最大磷吸附量 2 1。此法选用填 料的量较少且粒径与实际工程所用填料不尽相 同，填 料磷去除性能可能与实际处理工程不符。为科学合理 选择人工湿地磷去除填料，本文通过填料等温磷 吸附 实验、室内小试填料磷去除实验综合比较不同湿地填 料磷去除效果，开展湿地填料磷去除性能的影响因素 分析。1 实验材料及方法 1 1实验材料及填料部分物理性质的测定 考虑到填料来源及成本等因素，选用天津蓟县 的 页岩、粗砾石、铁矿石、麦饭石作为实验材料。测定填 料干容重、孔隙率、水力渗透系数、颗粒粒径分布等部 分物理性质见表 1。其中填料干容重和孔隙率采用标 准土壤科学方法测定_ l 3 1，水力渗透系数依据达西渗透 原理测定 ll 2 1，根据 0 2 5 2 0 mm共 7种孔径筛 网的筛 分结果确定填料的颗粒粒径分布。1 2等温吸附实验 准确称取 1 0 5下干燥 2 h粒径为 0 5 1 0 m m 的页岩、粗砾石、铁矿石、麦饭石 1 g于 1 0 0 mL具塞 离心管 中，分 别准确加入含磷 浓度为 0 0、1 0、2 5、5 0、1 0 0、2 5 0、5 0 0、1 0 0 0、2 0 0 0 mg L 的 K H 2 P O 4 溶 液 2 0 m L。各离心管 中加入 1滴 甲苯防止微生物活 动。在水浴恒温振荡器上连续振荡 2 4 h(2 5，1 7 5 r m i n )，振荡完毕后在 5 0 0 0 r m i n 下离心 1 5 mi n，取上清液用钼锑抗比色法测定磷含量，溶液 中磷含量 减小量即填料的磷吸附量。以溶液平衡吸附浓度对填 料平衡吸附量作图即填料等温磷吸附曲线。为确定填料颗粒粒径对磷吸附效果的影响，选用 粒径 分别 为 0 0 0 5 m m、0 5 1 0 m m、1 0 2 0 mm的 页岩测定不同粒径下页岩磷吸附量的变化，测定步骤 同上。1 3小试实验装置简介及运行 1 3 1 小试实验装置简介 采用室内无植物复合垂直潜流填料柱测定 比较 填料磷去除效果，实验装置共 4套，分别装填页岩、麦 饭石、粗砾石和铁矿石。每套装置(图 1)由污水配给 箱、下行柱和上行柱(按污水流向)等构成，下行柱和 上行柱的内径、柱高相同，具体是 1 5 e m和 1 2 0 e m，材 质为有机玻璃。进水管位于下行柱 1 0 0 e m处，出水管 位于上行柱 9 0 e m处，比进水管低 1 0 e m以保证污水 自重力流。下行柱与上行柱填料分别装填至与进水管 口和出水管 口等高的位置。1 3 2 室内小试运行 系统 自2 0 0 6年 3月中旬于室内正常运行，温度 控制为 2 0 2 5。设计水力负荷(H L A)为 2 0 0 m m d 污水配给箱进水闸 图 1人工湿地小试实验装置简图 Fi g ur e 1 S c h ema t i c d i a g r a m o f pi l o t s c a l e c on s t r u c t e d we t l a n d e x p e r i me n t a l a p p ar a t u s 表 1填料部分物理性质 Ta b l e 1 S o me p hy s i c a l pr op e r t i e s o f t he s ubs t r a t e i n t h e c o n s t r uc t e d we t l a n ds 维普资讯 http:/ 7 5 0 汤显强等：人工湿地填料磷去除效果的影响因素分析 2 0 0 8 年 3 月(各装置 日进污水约 7 L)，为研究 H L A变化对湿地 填料除磷效率 的影响，一个月后，H L A降至 1 4 2 8 m i l l d 一 1(各装置 日进污水约 5 L)。取卫津河河水作为 小试污水水源(贮于室内污水箱，容积为 3 m s)，主要 水 质参数 C O D：1 1 5 1 4 0 mg L ，T N：7 8 9 5 m g L ，T P：1 4 1 6 m g L 一，p H：8 3 2 8 6 4。每 隔 2 d在 进 水 口，出水口采集水样，测定水样中磷浓度变化。2 实验结果与讨论 2 1湿地填料的物理性质 从表 1 看出，页岩和粗砾石的孔隙率较高，铁矿 石与粗砾石的干容重较大。水力渗透系数 可用来 比较不同填料对湿地水流的滞留能力，值越低，填 料与污水的接触效果越好t l 2】。粗砾石的水力渗透系数 最大，其次是铁矿石，麦饭石与页岩较低，分别为 6 和8 c m -。d 6 o d 。为填料颗粒的均匀性参数，其值越 低，颗粒均匀度越高，比表面积越大。粗砾石与麦饭石 颗粒具有较高的均匀度，铁矿石与页岩均匀性较差。2 2填料类型、粒径对磷吸附性能的影响 图 2为页岩、粗砾石、铁矿石、麦饭石的等温磷吸 附曲线。填料类型影响填料磷吸附效果，不同填料磷 吸附效果为页岩 铁矿石 麦饭石 粗砾石。图 3表 明填料粒径影响填料最大磷吸附量，当页岩颗粒粒径 从 0 0 0 5 mm、0 5 1 0 m m增加至 1 0 2 0 mi l l 时，其 亲 謦 塑 登 溶液平衡浓度 C e mg L-i 图 2 填料对磷的等温吸附曲线(0 5 m md 1 mm)F i g u r e 2 C u r v e s o f i s o t h e r m a d s o r p t i o n o f p h o s p h o r u s b y t h e s u b s t r a t e f 0 5 mmd 1 m m)謦 谪 莲 譬 溶液平衡浓度 C e m g L-图 3不同粒径页岩等温磷吸附曲线 F i g u r e 3 C u r v e s o f i s o t h e r m a d s o rpt i o n o f p h o s p h o r u s b y s h a l e w i t hd i ff e ren t g r a i n s i z e s 最大磷 吸附量分别从 9 5 3 8 m g P k g 一 页岩、6 1 9 7 m g P k g-l 页岩下降至 3 8 3 8 m g P k g-l 页岩，随着填料 粒径的增大，页岩最大磷吸附量逐渐变小。国外学者 S e o 等研究不同牡蛎壳填料的磷吸附性能时，也发现 填料粒径增加会导致填料最大磷吸附量减小。填料 粒径增加使颗粒比表面积减小，引起单位质量填料表 面磷吸附点位减少可能是填料最大磷吸附量降低的 重要原 因 1 3 1 4】。为分析填料类型对磷吸附性能的影响，分别采用 L i n e a r(线性)、F r e u n d l i c h(弗兰德里希)和 l a n g mu i r (兰格缪尔)吸附模型模拟拟合图2 4 种填料磷吸附 过程。L i n e ar，F r e u n d l i c h，l a n g m u i r 方程如下：J N=p C+q,N=K C ，-等。式中：为填料单位质量吸附量，mg P k g t 填料；C为 溶液 中磷平衡浓度，mg L t；队q分别为 L i n e ar方程的 斜率 与截距，L k g 一，m g P k g t；b填 料的饱 和吸附量，mg P k g-填料；，k为常数。从表 2中拟合参数 z 看，L i n e ar拟合结果相对 较差，不适合模拟上述 4 种填料等温磷吸附曲线。页 岩、铁矿石与麦饭石 3 种填料的l a n g m u i r 拟合结果优 于 F r e u n d l i c h 拟合结果。国外的 D r i z o，S a k a d e v a n 等报 道过类似现象，与 F r e u n d l i c h拟合结果相 比，页岩、铁 表 2填料等温吸附曲线拟合参数表 T a b l e 2 T h e s i mu l a t i o n p a r a me t e r s o f i s o t h e rm a d s o rpti o n c u rve s b y the s u b s t r a t e s 瑚 姗伽姗抛 o 维普资讯 http:/ 第 2 7 卷第 2 期 农 业 环 境 科 学 学 报 7 5 l 矿渣、角闪石的磷吸附过程更适合用 l a n g m u i r 吸附等 温线描述 8,1 1】，其给出的解释是可能矿物材料(页岩、铁矿渣、角闪石)颗粒表面均匀度较高，比较符合 l a n g m u i r 单分子层吸附假设。4种材料中粗砾石磷吸 附过程最符合 F r e u n d l i c h吸附模 型，拟合相关度达到 0 9 9 2 4，这说明粗砾石 的磷 吸附量随着吸附平衡浓度 的增加呈指数下降趋势，磷 吸附过程不适合 l a n g m u i r 吸附模型。2 _ 3水力负荷变化对 室内小试填料磷 去除性能的影 响 图 4为水力负荷为(HL A)2 0 0 mm d 和 1 4 2 8 m m d 时页岩、粗砾石、麦饭石 和铁矿石 T P去除率变化 曲线(竖线为 H L A变化分界线)。一一页岩 一 铁矿石 一 粗砾石 一 麦饭石-、I1 -“t 一0 进水浓度 rI r H -_ _ -，-一 HL A=2 0 o mm d HLA=1 4 28 mm d 。o、b、。一 o。，o_ o_ o 一 r TVT t 一 1 一 ，、，、一下 。j：I|：C 嫩二 二 3 2 5 4 2 4 1 0 4 2 O 4 2 8 5 6 5 一l 4 日期 月一 日 图 4 室内小试人工湿地填料 T P去除率 F i g u r e 4 TP r e mov a l r a t e s o f i nd o o r p i l ot s c a l e c on s t r u c t e d we t l a nd 水力负荷(H L A)变化影响填料 T P 去除效率。当 H L A从 2 0 0 m m d 降至 1 4 2 8 mm d 时 粗砾石、铁 矿石 T P去除率呈现上升趋势。水力负荷下降导致水 力停留时间增加 污染物与填料接触时间延长 粗砾 石、铁矿石 T P去除率相应出现上升趋势，其中粗砾 石 T P去除率上升较稳 定 增加大约 1 0：铁矿石填 料 T P去除率上升较快，从大约 6(H L A=2 0 0 m m d )上升至 3 5 左右(H L A=1 4 2 8 m m d )。受水力负荷降 低的影响 F e 与磷 的结合效果不断改善可能是铁矿 石 T P去除率大幅上升的另外一个可能原因。根据测 定进 出水 p H值发现(图 5)，铁矿 石填料 出水 p H平 均值从 8 2 6(H L A=2 0 0 m m d )降低至 8 0 2(H L A=1 4 2 8 m m d 一)，说 明随着 HL A降低，填料 与污水、微 生物等接触时间增长 铁矿石可能受微生物及其他 因 素作用，其与磷的结合效果不断改善I l l I，T P去除 I生能 不断提高。但碱性环境下，乙酸提取实验表明 与铁 磷、铝磷相 比，钙磷是磷 的主要结合形态 l 1 3 与铁相 8 4尊LA=200mmd-。毛：。_0-。-。-o，。：叶 ：。_o。J 比，钙在铁矿石中的含量一般较低。因此铁矿石填料 T P去除率不断上升是否还受其他因素影响尚需进一 步研究。水力负荷(HL A)从 2 0 0 mm d 降至 1 4 2 8 mm d 后，页岩 T P去除率保持相对稳定，麦饭石不断下降。页岩填料磷吸附量较大(图 4)，当进水磷浓度较低且 稳定时，填料磷 吸附稳定且不易饱和。因此页岩填料 T P去除率保持相对较稳定，并处于 8 0 以上。麦饭石 T P去除率随着 H L A的降低而降低的原因可能是填 料吸附已达到饱和。处理初期，麦饭石填料 由于颗粒 粒径较小且分布均匀(d 。，d 6 o d 。值较小)，填料具有 较好的表面磷吸附效果，磷去除I生能较好。但根据表 2麦饭石等温磷吸附曲线的拟合情况看 麦饭石最大 磷 吸附量较小且符合 L a n g m u i r 吸附模 式。L a n g m u i r 吸附为表面单分子层吸附，具有饱 和吸附量，因此 当 麦饭石填料吸附达到饱和后，表面不再提供新 的磷吸 附点位。此时降低水力负荷并不能提高填料 的 T P去 除效率。3 结论 本文选择页岩、粗砾石、麦饭石与铁矿石作为实 验材料。通过填料等温磷吸附实验、室 内小试填料磷 去除实验分析人工湿地填料磷去除性能的影响因素。填料等温磷 吸附实验结果表明湿地填料类型和填料 粒径均影响填料磷去除效果。水力负荷是影响湿地填 料磷去除效果的另一因素。当 H L A从 2 0 0 m m d 降 至 1 4 2 8 mm d 时，粗砾石、铁矿石构建湿地 T P去除 率呈现不 同程度上升，页岩 T P去除率保持相对稳 定，麦饭石接近饱和吸附后 T P去除率随 HL A减小不 但没有增加，反而不断下降。一_ r I k 爱*d II 8 7 6 5 4 3 2 如 加 m 0 踌(I LL 维普资讯 http:/ 7 5 2 汤显强等：人工湿地填料磷去除效果的影响因素分析 2 0 0 8 年 3 月 本文仅 限于填料磷去除性能 比较和影响因素分 析，没有考虑微生物、湿地植物、污水类型等对湿地填 料磷去除性能的影响，基于本文实验结果，将选择页 岩为湿地磷去除填料进一步研究湿地植物、微生物等 因素对人工湿地磷去除性能的影响。参考文献：【1】D r i z o A，C o m e a u Y，F o r g e t C，e t a1 P h o s p h o r u s s a t u r a ti o n p o t e n t i a l：A p a r a me t e r f o r e s ti ma t i n g t h e l o n g e v i t y o f c o n s t r u c t e d w e t l a n d s y s t e ms J En v i r o n Sc i Te c h n o l，2 0 0 2，3 6：4 6 4 2 4 64 8 【2】R i v e r a R，Wa r r e n A，C u r d s C R，e t a1 T h e印p l i c a t i o n o f t h e r o o t z o n e me tho d f o r the t r e a t me n t and i l s e o f h i g h s n g t h a b a t t o i r w a s t e i n Me x i c o J Wa t S i c T e c h，1 9 9 7，3 5：2 7 1-2 7 8 【3】Wo r r al l P，P e b e r d y K J，Mi l l e t M C C o n s t r u c t e d w e t l and s and n a t u r e c o n s e r v a t i o n 忉 Wa r S i c T e c h，1 9 9 7，3 5：2 0 5 2 1 3 【4】朱彬刑 用水生植物净化富营养化水体的研究进展叨 上海环境 科学，2 0 0 2，2 1(9)：5 64 5 6 7 5】汤显强，黄岁棵 人工湿地去污机理及 国内外应用现状 J】_ 水处理技 术，2 0 0 7，3 3(2)：9 1 3 【6】G r e e n w a y M，Wo o l i e y A C o n s tr u c t ed w e t l a n d s i n Q u e e n s l and：p e r f o r-ma l i c e e ffic i e n c y a n d n u t r i e n t b i o a c c u mu l a t io n J】E c o l E n g，1 9 9 9，1 2：3 9 -5 5 7】T ann e r C C，S u k i a s J P S，U p s d e U M P R e l a t i o n s h i p s b e t w e e n l o a d i n g r a t e s an d po l l u t a n t r e mo v a l d u rin g ma t u r a t i o n o f g r a v e l-b e d c o n s t r u c t ed w e t land s J E n v i r o n Q ，1 9 9 8，2 7：4 4 8-4 5 8 【8】D i a z O A，R e d d y K R，Mo o r e P A S o l u b il i t y o f i n o r g ani c P i n s t r e a r n w a-t e r a s i n fl u e n c ed b y p H and C a c o n c e n a tl o n 叨 Wa t e r Re s，1 9 9 4，2 8：1 7 5 5 1 7 6 3 【9】P ant Hl(Redd y KR，L e mo nE P h osp h o rus ret e n ti o n c a p a c i tyo f r o o t b e d med i ao f s u b s u r f a c eflo w c o n s t r u c t ed we tl an d s Ec o l o g i c a l En g i n e e r-i n g,2 0 01，1 7：3 4 5-3 5 5 【1 o】谭洪新，周琪 湿地填料的磷吸附特性及潜流人工湿地除磷效果 研究叨 农业环境科学学报，2 0 0 5，2 4(2)：3 5 3 3 5 6 【1 1】S a k a d e v an B a v o r H J P h o s p h a t e a d s o r p t i o n c h a r a c t e ri s t i c s o f s o il s，s l a g an d z e o l i t e t o b e u s ed a s s ub s t r a t e s i n c o n s tr u c t ed we tl a n d s y s t e ms J】Wa t e r R e s，1 9 9 8，3 2，2：3 9 3-3 9 9 【1 2】D ri z o A，e t a1 P h y s i c al-c h e m i c al s c r e e n i n g o f p h o s p h a t e rem o v i n g s I l b s t r a t e s fo r u s i n g i n c o n s t r u c t ed w e tl and s y s t e m s叨Wate r Re s，1 9 9 9，33，1 7：3 5 9 6-3 6 0 2【1 3】P r o c h a s k a C Z o u b o u H s AL R e m o v a l o f p h o s p h a t e s b y p il o t v e r t i c a l-fl o w co n s t ru c t e d we tl a n d s u s i n g a mi x t u re o f s an d an d d o l o mi t e 8 s s u b s t r a t e J】E c o l o g i c a l En g i nee r i n g,2 0 0 6，2 6：2 9 3-3 0 3 f 1 4 S e oD C，C h o J S，L e eH J e t a1 P h o s p h o rusret e n t i o n c a p a c it y o f fi l t e r medi a for e s t i ma t i n g the l o n g e v i ty o f con s t r u c t e d we tl an d Wate r Re s，2 0 o 5 3 2 4-4 5-2 4 57 维普资讯 http:/