最新消息: USBMI致力于为网友们分享Windows、安卓、IOS等主流手机系统相关的资讯以及评测、同时提供相关教程、应用、软件下载等服务。

基于集对分析的污水处理工艺设计优化

IT圈 admin 25浏览 0评论

2024年6月5日发(作者:赫丹山)

维普资讯

第34卷第11期 

湖南大学学报(自然科学版) 

V01.34.No.11 

NOV.2 0 0 7 

2 0 0 7年1 1月 

Journal of Hunan University(Natural Sciences) 

文章编号:1000.2472(2007)11—0070.06 

基于集x,3 析的污水处理.T。艺设计优化 

蒋 茹 ,曾光明,李晓东,杜春艳,李建兵,宋琳玲 

(湖南大学环境科学与工程学院,湖南长沙410082) 

摘 要:针对污水处理工艺设计多目标、多指标,以及确定性与多种不确定性共存的特 

点,将集对分析理论用于污水处理工艺设计的优化决策,拓展了结构函数一同异反联系度 

的概念;在确定权重时引入熵权修正广泛应用的层次分析法的主观权重,一定程度上解决了 

信息丢失问题;其独特的稳定性分析通过基序分析其稳定区域,判定扩展序的存在,保证了 

决策的准确性和可靠性.方法的有效性通过实例研究得到了验证.实例考察了某污水处理厂 

的4种候选方案,即 /o法、三沟式氧化沟、厌氧一单沟式氧化沟和SBR法.结果表明,对该 

污水处理厂而言,厌氧一单沟式氧化沟为最佳方案,综合效益最高.所得排序稳定,不存在扩 

展序. 

关键词:污水处理;工艺设计;优化;集对分析;熵权 

中图分类号:X32、021 文献标识码:A 

Wastewater Treatment Process Design Optimization 

Based on Set Pair Analysis 

JIANG Ru ,ZENG Guang—ming,LI Xiao-dong,DU Chun—yan,LI Jian—bing,SONG Lin—ling 

(College of Environmental Science and Engineering,Hunan Univ,Changsha,Hunan 410082,China) 

Abstract:Considering the characteristics of multi—objective,multi—index and coexistence of certainty and un— 

certainties,set pair analysis theory was introduced into wastewater treatment process design optimization.The 

concept of structural function,that is,identity—discrepancy—contrary connection degree was developed.Further— 

more,a subjective and objective integrated approach to determine weights was applied to avoid information loss. 

The unique stability analysis for basic ranking to determine other extra ranking could ensure the veracity and sta— 

bility of decision results.The effectiveness of the approach was verified through case study.For some municipal 

wastewater tratmente plants,four alternatives( /o,triple oxidation ditch,anaerobic single oxidation ditch 

and SBR)were evaluated and compared.The results showed that anaerobic single oxidation ditch WaS the opti— 

mal scheme and could obtain the maximum general benefits.The resultant ranking was stable and there was no 

extra ranking. 

Key words:wastewater treatment;process design;optimization;set pair analysis;ent1 opy weight 

污水处理工艺设计直接影响着项目建设能否达 

到预期的目标以及投人运营后的效果,如污水处理 

的经济性、可靠性和稳定性.目前采用的方法有层次 

分析法(AFIP) 卜引、灰色关联法[引、灰色局势决策 

法 和模糊综合评价法 等,但还存在以下不足: 

AHP常因对标度把握不准而无法构造良好的判断 

÷收稿日期:2007-02_o2 

基金项目:国家自然科学基金(70171055、50179011)、国家杰出青年科学基金(50225926)、高等学校博士学科点专项科研基金 

(20020532017)资助项目;2000年教育部高等学校优秀青年教师教学科研奖励计划资助项目 

作者简介:蒋茹(1979一),女,浙江金华人,湖南大学博士研究生 

t通讯联系人.E—mail:hunan_jiangru@163.com 

维普资讯

第11期 蒋 茹等:基于集对分析的污水处理工艺设计优化 71 

矩阵;灰色方法考虑到灰色性的特点,但指标值分布 

过于离散时,由于白化函数的范围窄,可能导致最终 

优化排序的错误;模糊综合评判法由于基本运算规 

对环境的影响等 

2工艺设计优化决策的特点 

2.1 多目标性与目标的相互矛盾性 

则强调取小取大,过分突出极值的作用,忽略了各项 

指标的综合效应,且模糊隶属度的确定存在较大的 

随机性,这些都可能导致优化结果的不确定性. 

集对分析是我国学者赵克勤于1989年提出的 

污水处理工艺优化是一个涉及经济、技术、环 

境、社会文化等方面的多目标、多层次决策过程.而 

种新的处理不确定性信息的系统分析方法 一 . 

近十年来,随着理论的不断完善和人们对环境不确 

定性问题研究的日益重视,已成功地用于城市规划、 

质量评价、污染预报、资源开发利用等方面.集对分 

析的核心思想是把确定性和不确定性作为一个互相 

联系、互相制约、互相渗透,在一定条件下可互相转 

化的确定一不确定系统来处理.它用联系度来统一 

处理模糊、随机、信息不完全所导致的系统不确定 

性.针对污水处理工艺优化决策的多目标、多指标特 

点-s ,以及确定性和多种不确定性并存的现象,笔者 

应用集对分析法对污水处理工艺设计优化的问题进 

行了研究. 

1工艺方案优化的评价指标体系 

1.1经济指标 

主要表现为工程建设投资和运行费用.前者包 

括土地费用、直接建设安装工程费用、设备费用、管 

理费用和财务费用;后者包括电费、人工费、药剂费、 

维护费、污泥处理以及折旧费等.经济指标可用诸如 

单位水量处理费用等综合指标来表示. 

1.2工艺指标 

主要包括单位时间处理水量、进水水质、处理效 

率、出水水质等与处理工艺相关的指标. 

1.3可靠性指标 

可靠性指标包括所选用工艺技术的新颖性及先 

进性、技术的成熟程度、系统对水量水质变化的适应 

性、设备的耐久性、设备维修的难易程度等. 

1.4管理指标 

主要有操作技术的复杂程度、操作劳动的繁简 

与强度、控制与调节的简便程度、设备和器材的标准 

化程度、药剂的通用性程度及普及性、运输量以及施 

工期的长短等. 

1.5资源占有及环境指标 

包括占地面积、占地的价值级别、劳动指标、副 

产品及甲烷气的回收利用、污水及污泥的利用程度、 

电耗及其他能耗、污泥产量、药剂耗量、二次污染及 

这些目标有时具有不同重要性甚至是相互矛盾的. 

尤其是近年来,对出水水质和处理水量等提出了较 

高的要求,但能达到要求的此类工艺要求在资金和 

占地方面大投人.另一方面,在我国目前经济发展的 

情况以及城市用地日益紧张的形势下,要求污水处 

理厂的建设占用土地面积尽可能地小、投资尽可能 

地少.因此协调好这些相互矛盾的多目标是我们在 

实际工作中需要解决的一个主要任务. 

2.2多层次和多指标 

从文献[3]中的递阶层次结构可看出,由于评价 

系统的复杂性,工艺设计优化呈现的是一个多层次 

结构,分别由目标层、准则层和指标层构成. 

2.3定性指标与定量指标相结合 

在评价指标中既有定量指标,又有定性指标,前 

者可通过企业的统计数据得到,后者需要利用专家 

法得到具体的指标值.这些指标又可分为越大越好 

的效益型指标和越小越好的成本型指标. 

2.4评价指标的相互交错性 

在影响选择的各因素之间并非相互独立,它们 

之间存在依赖性和反馈性,相互依存、相互影响.如 

果忽略它们的差异性以及之间的关系则会得到与实 

际情况不太吻合的排序结果. 

2.5不确定性 

污水处理工艺优化涉及的因素多且复杂,相当 

部分评价因素不能定量描述,具有不确定性(模糊 

性、灰色性、未确知性).同时处理系统各种行为的输 

人、输出也不是确定的数值.如果在决策过程中无视 

上述不确定性,将会损失大量的信息,导致决策的不 

准确和不科学.此外,各指标的影响程度也不相同, 

其权重的确定存在较大的主观不确定性. 

综上所述,污水处理工艺优化的过程是定性与 

定量、确定性与多种不确定性同时存在的多目标决 

策问题.而集对分析理论用联系度来统一处理模糊、 

随机、信息不完全所导致的综合不确定性,并作为一 

个确定一不确定系统来研究.这是将集对分析理论 

引人污水处理工艺优化决策过程的原因. 

维普资讯

72 湖南大学学报(自然科学版) 

合评价水平的不变 ¨.而污水处理工艺优选涉及的 

3现有优化策略存在的主要问题 

利用集对分析中的重要概念——联系度,来构 

建工艺设计优化的效用函数.根据不同的定义,分别 

形成两种不同的优化策略. 

优化策略一 :根据给定的各个方案,确定出一 

个理想方案,将各待选方案与理想方案构筑一个集 

指标,很多是相互影响、相互联系的,因此采用各待 

选方案指标值与理想方案指标值同一度的线性加 

权,显然不能反映实际. 

2)指标的权重是相互影响的,指标权重的赋值 

与指标的量化相比具有更大的不确定性和随意性. 

目前普遍采用的赋权方法是AHP法.AHP法是一 

种定性与定量相结合的决策分析技术,将问题分解 

成若干层次,由专家和决策者对所列指标通过两两 

对,计算备选方案各指标与理想方案各对应指标的 

同一度、差异度和对立度,并计算集对的联系度,寻 

找与理想方案最贴近的那个方案;或者单一地利用 

联系度中的同一度分量,结合各指标的权重计算待 

选方案与理想方案的加权同一度,根据所得的同一 

度数值大小确定方案的优劣次序. 

优化策略二 们:将待评价的评价指标量与评价 

标准构筑一个集对,对于一个具体决策来说,设有 

N个评价指标,其中有s,F。,F:,P个指标分别处 

于I,Ⅱ,Ⅲ,Ⅳ级标准,则该决策的联系度表达式为 

= +

争。+ :+ . 。,f:分别表示差异不确 

C 

定度, 表示对立度标记.设a= ,b。= ,b:= 

,c=丽/-

n 

则a,b (,z=1,2),c依次表示为同一 

度、差异不确定度、对立度.原联系度表达式可简写 

为 :a+b。 。+b2 i2+ .根据集对分析理论,同 

度、对立度是相对确定的,而差异度则相对不确 

定,a,b,c三者是对同一问题不同侧面的全面刻 

画,依据a,b,c三者大小关系及定量分析,可判断 

实际决策目标状况.a值越大,b,c值越小,工艺级 

别越优先.另外,根据a,b,c的具体数值,可对污水 

处理厂工艺状况进行分级,即以评价因子的数量指 

标相对于决策评价标准的达标、超标及其所占比例, 

确定污水处理厂最佳工艺.但是在实际工作中,某项 

评价指标的评价标准没有确切的定义,由于地区差 

异性,很难实现统一. 

总的来说,现有的应用主要存在以下不足: 

1)线性加权对于各评价指标相互联系的污水处 

理工艺方案优化不合理.对于文献中的对各方案的 

同一度进行线性加权,一般只适用于各评价指标相 

互独立的场合,若各评价指标间不独立,“和”的结果 

必然是信息的重复,难以反映客观实际;特别是线性 

加权法可使各评价指标间得到线性地补偿,任一指 

标值的减少,都可以用另一些指标的增量来维持综 

比较重要程度而逐层进行判断评分,利用计算判断 

矩阵的特征向量确定下层指标对上层指标的贡献程 

度,从而得到最底层指标对总目标而言重要性的赋 

权结果 .但AHP确定权重存在较大的主观随意 

性,当采用专家咨询方式时,容易产生循环而不满足 

传递性公理,导致标度把握不准并丢失部分信 

. 

3)我们认为联系度表达式原则上是一个结构函 

数,直接将联系度作为一确定的数值来进行各种运 

算是值得商榷的.若只单一的利用联系度中的同一 

度分量,受理想点的限制,尚未考虑对立度,具有一 

定的主观性和片面性,实质为一种理想点法. 

4基于集对分析的污水处理工艺优化多属 

性决策模型 

多属性决策问题一般可记为G=(S,E,W, 

D).方案集S={S。,S2,…,S },指标集E={e。,ez, 

,e },指标权重集W={ 。, :,…, },决策矩 

阵D:(d ,) ,dh为方案S 关于指标e 的属性 

值.决策矩阵经归一化处理,所得的归一化矩阵记为 

,, 

P=(P ,) .对效益型指标,有Ph=盖,d = Ⅱ

一 

max{d。,,d:,,…,d };对成本型指标,有Pk,= 

dml- 

d =min{dl ,d2,,…,d ,}. 

Ⅱ 

4.1熵权法确定指标权重 

指标权重以数量形式反映各指标在综合评价中 

的重要程度.若同一指标下各方案的指标值相差越 

大,则反映该指标在综合评价中比较重要.若某项指 

标值全部相等,则该指标在综合评价中不起作用.信 

息论中,信息熵是信息量的度量,其基本出发点是把 

获得的信息看作用以消除不确定性的东西.故信息 

量的大小可用被消除的不确定性的多少来表示.某 

项指标值变异程度越大,该指标提供的信息量越大, 

则该指标的权重也越大,反之亦然.因此可以根据各 

维普资讯

第11期 蒋茹等:基于集对分析的污水处理工艺设计优化 73 

指标提供的信息量大小来确定指标权重. 

利用熵权和普遍应用的AHP法确定权重,其出 

发点和理论基础均不同,前者从各指标提供的信息 

量多少出发,依据的是信息论;后者是由决策者分析 

6 ==-——cz ,——c ,=: — ! — :: ! ; 

(4) 

因此,对于效益型指标,集对{P “ }的联系度为 

评价,人为主观色彩较重,理论基础是决策论.本文 

采用主观与客观相结合的方法,在AHP法基础上, 

运用熵值理论对各专家评估水平赋予权重,确定各 

指标的组合权重.科学的确定权重方法应采用两者 

的算术平均值_1引. 

为确定各指标的权重,先计算各指标的熵值 

E,=一(Inm) ∑P lnPh (1) 

式中只,=Pk /∑/ ,若 =o'贝IJ ,=0. 

其次,计算各指标的偏差度D,=1一E,和客观 

加权因子 ,=D,/∑Di.设用AHP法得到的各指 

标的权重为 (r=1,2,…,7z),可用 ,进一步进行 

修正,即组合权重为 

硼 =( ,+ )/2 (2) 

4.2 基于集对分析的多属性决策模型n 

设最优方案和最劣方案分别记为U,V,它们 

可根据系统目标和客观条件,从方案集的内部或外 

部选定.U,V的选定涉及到指标类型,记诸方案关 

于P 的比较区间为[ ,,“,].最优方案和最劣方案 

关于P,的属性值分别为 ,, ,.当 属于效益型指 

标时,有“,=“,, = ,;当P 属于成本型指标 

时,有一U,=73,,一73 =U,.最优方案和最劣方案构成 

了诸方案的比较区间[V,U]. 

1)在P,的比较区间中确定集对{P }的联 

系度. 

对效益型指标, r“ 十 ,_73 一,U南十 一73 分别表示在  

数值上P 与“,,P 与 ,的接近程度.其值越大,接 

近程度越大,在确定区间[ ,“,]内,P 与 的接 

近恰说明了Pb与“,的远离.故定义集对{P “,} 

的同一度和对立度分别为 

户h 

一 

“h一

“ + 

’ 

p “, , ,^、 

h 

j 

由归一化条件a+b+c=1,定义差异度为 

P^, .(“,一Pb)(P 一 ).

,,

} _ 十——— _ -十 

. 

“,+ 

(5)

)户^ 、J 

 

对成本型指标,可确定集对{P ,}的联系度为 

“, , 

(“,一Pb)(户 一 ,)..

ult,

I 十— 十 

 

k ・ 

鸟 (6)o  

式(5),(6)中的同一度和对立度分别说明了方案 

关于指标e,的趋优程度和趋劣程度.两式关于同一 

度和对立度互为对称,且具有相同的差异度,说明了 

刻划的对立性和一致性.当P =“,或 ,时,差异度 

最小为零,表示刻划的不确定性最小;当Pb= 

时,差异度最大为2 /(“,+ ,),表示 

刻划的不确定性最大(但不超过0.5). 

将以上两式统一记为 

{d.. }=ab+6h +c^ (7) 

2)在比较区间[V,U]中计算集对{ ,U}的 

联系度 

{ .u}=a^+bk +c (8) 

式中 =∑ n打,bk=∑ , =∑ ,. 

r=l r=1 ,=l 

3) 与最优方案U的贴近度 

=— 一

. 

(9) 

a^1.C^ 

显然, 最大所对应的方案为最优方案,此时 

得到的排序为基序. 

4.3 稳定性分析 

当 >0,对a 正向修正,表示对方案 接近理 

想方案U持肯定意见,且 越接近1,修正作用越 

强.反之,当i<0,对 正向修正,表示对方案 接 

近理想方案U持否定意见,且 越接近一1,修正作 

用越强 引.因此,可通过分析 在[一1,1]内的变 

化,对方案排序的稳定性进行分析.为清楚表达,记 

i为△ ,则△ ∈[一1,1].假定对于整个方案集,△ 

相同.设 基序名次为1,记原 为r ;当△ ≠0 

时,其值变为 .设基序中有r:<r ,则根据式(9) 

可得 

维普资讯

74 湖南大学学报(自然科学版) 

cpa^一cka >0 (10) 

表2 各工艺设计方案的联系度及相对贴近度 

Tab.2 Weights of attributes by AItP and entropy method 

1)当0≤Ai≤1时,式(9)变为 

一(b 

n^+kAi)+ 

: 

(11)

为维持基序r:<r ,需满足不等式Cpa 一cka > 

(ckbp—cp )Ai.当ckbp一 6 ≤0时,有△ ∈[0, 

1];当ckbp—cp ≤0时,有0≤Ai≤ 

min{ , }. 

2)当一1≤△ <0时,式(9)变为 

一n^+(毛 bc^一 ^  )Ai.’ (12) 

同理,为维持基序 h<r ,有当akb 一apb ≥ 

0时,有△ ∈[一1,0);当akb 一apbl<0时,有 

max

【-1, }≤Ai<o0 .

取以上两种情况的并集即可得到基序排序 

<r:的稳定区域.以上进行的是局部分析,为了保 

证全序的稳定性,只需依次两两比较分析确定△ 的 

范围,则所有△ 取值区间的交集即为排序结果的稳 

定区域. 

4.4 实例研究 

拟建一座1O万m3/d的污水处理厂,二级生化 

处理,并要求脱氮除磷.有4个方案: /o法(S。)、三 

沟式氧化沟(S )、厌氧一单沟式氧化沟(S。)和SBR 

法(S ).所选定的评价指标共8个,分别为项目投 

资、经营成本、占地面积、氮磷去除效果、污泥处理效 

果、运行稳定性、工艺成熟性以及操作的难易程度. 

各方案的具体指标值可参见文献[3]. 

各指标权重、各方案的联系度及相对贴近度计 

算结果分别见表1和表2. 

表1 各指标的权重值 

Tab.1 Detailed infomat|0n 0f attribⅡtes 

根据表2中各方案的相对贴近度rJ,得到在相 

对确定条件下的排序结果为r >r >r。>r .因 

相对贴近度越大,其所对应的方案越优,所以方案优 

劣顺利分别为S,,S:,S ,S ,即厌氧一单沟式氧化沟、 

三沟式氧化沟、 /o法、SBR法,此为基序.该排序 

结果与文献[3]所采用的多层次灰色关联法得到的 

排序结果基本一致,只是在第2和第3位序上存在 

偏差. 

为了保证本文所得结果的稳定性,以基序为基 

础进行稳定性分析,其结果见表3.排序相邻方案的 

稳定性定量分析计算表明,得该排序结果的稳定区 

域为△ ∈[一1,1],这与i本身的定义区域是一致 

的.也就是说,不存在扩展序,通过本方法得到的排 

序结果是稳定的,其排序结果优于文献[3]的结果. 

表3排序稳定性分析 

Tab.3 Connection degrees and relative closeness 

degrees of each alternative 

5 结 论 

本文引入集对分析理论解决污水处理工艺方案 

优化的多属性决策问题,概念清晰,计算简便,结果 

稳定,便于编程实现.同时,为了减少决策者因受专 

业知识的限制和偏好等主观因素的影响,充分利用 

方案所反映的固有信息,在确定权重时引入熵权修 

正层次分析法主观权重的思想,综合考虑主、客观因 

素的影响,使得决策结果更趋向合理和科学. 

参考文献 

[1] ELLIS K V,TANG S L.Wastewater treatment optimization 

维普资讯

第11期 蒋茹等:基于集对分析的污水处理工艺设计优化 75 

model for developing world.I:Model development[J].Journal 

of Environmental Engineering Division,ASCE,1991,117:501 

581. 

[2]ELLIS K V,TANG S L.Wastewater treatment optimization 

model for development world.II:Model testing[J].Journal of 

Environmentla Engineering Division,AsCE,1994,120:610— 

624. 

[3]ZENG G M,JIANG R,HUANG G H,et a1.Optimization of 

wastewater treatment laternatives selection by hierarchy grey re・ 

lational analysis[J].Journal of Environmentla Management, 

2007,82:250—259. 

[4]王玲玲,曾光明,黄国和,等.小区域污水处理的多目标规划设 

计探讨[J].环境工程,2004,22(6):88—89. 

WANG LL,ZENGGM,HUANGGH,eta1.Study onmulti— 

objective programming for wastewater treatment in small region 

[J].Environmentla Engineering,2004,22(6):88—89.(In Chi・ 

neae) 

[5]凌猛,杭世王君.城市污水处理厂工艺方案模糊决策方法的应 

用[J].给水排水,1998,24(3):6—9. 

LIN M,HANG S J.The application of fuzzy decision—making 

approach in the process altemativea of municipal wastewater 

traetment plnat[J].Water and Wastewater Engineering,1998, 

24(3):6~9.(In Chinese) 

[6]赵克勤,宣爱理.集对论——一种新的不确定性理论、方法与应 

用[J].系统工程,1996,14(1):18—23. 

ZHAO K Q,XUAN A L.Application of lf novel uncertainty theo— 

ry nad method——Set pair analysis[J].Systems Engineering, 

1996,14(1):18—23.(InChinese) 

[7]JIANGY L,XUCF,YAOY,et a1.A new approachfor repr. 

aeenting and processing uncertainty knowledge[C]//The 2003 

IEEE Intemational Conference on Information Reuse and Inte. 

gration.USA.2003:466—470. 

[8]ANNELIES J B,HEINZ A P,RALF 0,et a1.Indiactors for 

the sustainabiilty assessment of wastewater treatment systems 

[J].Urban Water,2002(4):153—161. 

[9]王国平,王洪光.集对分析用于污水处理厂的综合评价[J].江 

苏环境科技,2002,15(1):16—18. 

WANG G P,WANG H G.Comprehensive assessment in sewage 

treatment plant by integration and contrast analysis[J].Jiangsu 

Environmental Science and Technology,2002,15(1):16—18. 

(In Chinese) 

[1O]李凡修,梅平,陈武.多元集对模型在污水处理厂改造决策中应 

用[J].环境科学与技术,2004,27(6):92—94. 

LI FX,MEI P,CHENW.Set pairanalysismethodwithmulti・ 

factor connection number for plnaning of WTP reformation[J]. 

EnvironmentlaScience and Technology,2004,27(6):92—94. 

(In Chinese) 

[11]李香云,刘理才.一种集成的科技实力综合评价方法[J].科学 

学研究,2004,22(6):85—87. 

LI X Y.LIU L C.A literature review on the study and evalu・ 

ation of scientific and technological strength[J].Studies in cSience 

of cSience,2004,22(6):85—87.(In Chinese) 

[12]张佳林,尉晓君.基于AHP法的地下水体综合评价[J].湖南 

大学学报:自然科学版,2006,33(5):136—140. 

ZHANG J L,WEI X J.Comprehensive assessment of groundwa— 

tel"based on AHPEJ].Journal of Hunan University:Natural cSi. 

encea,2006,33(5):136—140.(InChinese) 

[13]方创琳,毛汉英.区域发展规划指标体系建立方法探讨[J].地 

理学报,1999,54(5):410—419. 

FANG C L.MA0 H Y.A system of indicators for regional de— 

velopment planning[J].Acta Geographiac Siniac,1999,54(5): 

410—419.(In Chinese) 

[14]张文泉,张世英,江立勤.基于熵的决策评价模型及应用[J]. 

系统工程学报,1995,10(3):69—75. 

ZHANG W Q,ZHANG S Y,JIANG L Q.Application of deci・ 

sion—making model based on entropy[J].Journal of Systems 

Engineering,1995,10(3):69—75.(InChinese) 

[15]张斌.多目标系统决策的模糊集对分析方法[J].系统工程理 

论与实践,1997(12):108—114. 

ZHANGB.Fuzzyset pairanalysisformulti.objcetive systemdeci— 

sion.making[J].Systems Engineering—Theory&Practice,1997 

(12):108—114.(InChinese) 

2024年6月5日发(作者:赫丹山)

维普资讯

第34卷第11期 

湖南大学学报(自然科学版) 

V01.34.No.11 

NOV.2 0 0 7 

2 0 0 7年1 1月 

Journal of Hunan University(Natural Sciences) 

文章编号:1000.2472(2007)11—0070.06 

基于集x,3 析的污水处理.T。艺设计优化 

蒋 茹 ,曾光明,李晓东,杜春艳,李建兵,宋琳玲 

(湖南大学环境科学与工程学院,湖南长沙410082) 

摘 要:针对污水处理工艺设计多目标、多指标,以及确定性与多种不确定性共存的特 

点,将集对分析理论用于污水处理工艺设计的优化决策,拓展了结构函数一同异反联系度 

的概念;在确定权重时引入熵权修正广泛应用的层次分析法的主观权重,一定程度上解决了 

信息丢失问题;其独特的稳定性分析通过基序分析其稳定区域,判定扩展序的存在,保证了 

决策的准确性和可靠性.方法的有效性通过实例研究得到了验证.实例考察了某污水处理厂 

的4种候选方案,即 /o法、三沟式氧化沟、厌氧一单沟式氧化沟和SBR法.结果表明,对该 

污水处理厂而言,厌氧一单沟式氧化沟为最佳方案,综合效益最高.所得排序稳定,不存在扩 

展序. 

关键词:污水处理;工艺设计;优化;集对分析;熵权 

中图分类号:X32、021 文献标识码:A 

Wastewater Treatment Process Design Optimization 

Based on Set Pair Analysis 

JIANG Ru ,ZENG Guang—ming,LI Xiao-dong,DU Chun—yan,LI Jian—bing,SONG Lin—ling 

(College of Environmental Science and Engineering,Hunan Univ,Changsha,Hunan 410082,China) 

Abstract:Considering the characteristics of multi—objective,multi—index and coexistence of certainty and un— 

certainties,set pair analysis theory was introduced into wastewater treatment process design optimization.The 

concept of structural function,that is,identity—discrepancy—contrary connection degree was developed.Further— 

more,a subjective and objective integrated approach to determine weights was applied to avoid information loss. 

The unique stability analysis for basic ranking to determine other extra ranking could ensure the veracity and sta— 

bility of decision results.The effectiveness of the approach was verified through case study.For some municipal 

wastewater tratmente plants,four alternatives( /o,triple oxidation ditch,anaerobic single oxidation ditch 

and SBR)were evaluated and compared.The results showed that anaerobic single oxidation ditch WaS the opti— 

mal scheme and could obtain the maximum general benefits.The resultant ranking was stable and there was no 

extra ranking. 

Key words:wastewater treatment;process design;optimization;set pair analysis;ent1 opy weight 

污水处理工艺设计直接影响着项目建设能否达 

到预期的目标以及投人运营后的效果,如污水处理 

的经济性、可靠性和稳定性.目前采用的方法有层次 

分析法(AFIP) 卜引、灰色关联法[引、灰色局势决策 

法 和模糊综合评价法 等,但还存在以下不足: 

AHP常因对标度把握不准而无法构造良好的判断 

÷收稿日期:2007-02_o2 

基金项目:国家自然科学基金(70171055、50179011)、国家杰出青年科学基金(50225926)、高等学校博士学科点专项科研基金 

(20020532017)资助项目;2000年教育部高等学校优秀青年教师教学科研奖励计划资助项目 

作者简介:蒋茹(1979一),女,浙江金华人,湖南大学博士研究生 

t通讯联系人.E—mail:hunan_jiangru@163.com 

维普资讯

第11期 蒋 茹等:基于集对分析的污水处理工艺设计优化 71 

矩阵;灰色方法考虑到灰色性的特点,但指标值分布 

过于离散时,由于白化函数的范围窄,可能导致最终 

优化排序的错误;模糊综合评判法由于基本运算规 

对环境的影响等 

2工艺设计优化决策的特点 

2.1 多目标性与目标的相互矛盾性 

则强调取小取大,过分突出极值的作用,忽略了各项 

指标的综合效应,且模糊隶属度的确定存在较大的 

随机性,这些都可能导致优化结果的不确定性. 

集对分析是我国学者赵克勤于1989年提出的 

污水处理工艺优化是一个涉及经济、技术、环 

境、社会文化等方面的多目标、多层次决策过程.而 

种新的处理不确定性信息的系统分析方法 一 . 

近十年来,随着理论的不断完善和人们对环境不确 

定性问题研究的日益重视,已成功地用于城市规划、 

质量评价、污染预报、资源开发利用等方面.集对分 

析的核心思想是把确定性和不确定性作为一个互相 

联系、互相制约、互相渗透,在一定条件下可互相转 

化的确定一不确定系统来处理.它用联系度来统一 

处理模糊、随机、信息不完全所导致的系统不确定 

性.针对污水处理工艺优化决策的多目标、多指标特 

点-s ,以及确定性和多种不确定性并存的现象,笔者 

应用集对分析法对污水处理工艺设计优化的问题进 

行了研究. 

1工艺方案优化的评价指标体系 

1.1经济指标 

主要表现为工程建设投资和运行费用.前者包 

括土地费用、直接建设安装工程费用、设备费用、管 

理费用和财务费用;后者包括电费、人工费、药剂费、 

维护费、污泥处理以及折旧费等.经济指标可用诸如 

单位水量处理费用等综合指标来表示. 

1.2工艺指标 

主要包括单位时间处理水量、进水水质、处理效 

率、出水水质等与处理工艺相关的指标. 

1.3可靠性指标 

可靠性指标包括所选用工艺技术的新颖性及先 

进性、技术的成熟程度、系统对水量水质变化的适应 

性、设备的耐久性、设备维修的难易程度等. 

1.4管理指标 

主要有操作技术的复杂程度、操作劳动的繁简 

与强度、控制与调节的简便程度、设备和器材的标准 

化程度、药剂的通用性程度及普及性、运输量以及施 

工期的长短等. 

1.5资源占有及环境指标 

包括占地面积、占地的价值级别、劳动指标、副 

产品及甲烷气的回收利用、污水及污泥的利用程度、 

电耗及其他能耗、污泥产量、药剂耗量、二次污染及 

这些目标有时具有不同重要性甚至是相互矛盾的. 

尤其是近年来,对出水水质和处理水量等提出了较 

高的要求,但能达到要求的此类工艺要求在资金和 

占地方面大投人.另一方面,在我国目前经济发展的 

情况以及城市用地日益紧张的形势下,要求污水处 

理厂的建设占用土地面积尽可能地小、投资尽可能 

地少.因此协调好这些相互矛盾的多目标是我们在 

实际工作中需要解决的一个主要任务. 

2.2多层次和多指标 

从文献[3]中的递阶层次结构可看出,由于评价 

系统的复杂性,工艺设计优化呈现的是一个多层次 

结构,分别由目标层、准则层和指标层构成. 

2.3定性指标与定量指标相结合 

在评价指标中既有定量指标,又有定性指标,前 

者可通过企业的统计数据得到,后者需要利用专家 

法得到具体的指标值.这些指标又可分为越大越好 

的效益型指标和越小越好的成本型指标. 

2.4评价指标的相互交错性 

在影响选择的各因素之间并非相互独立,它们 

之间存在依赖性和反馈性,相互依存、相互影响.如 

果忽略它们的差异性以及之间的关系则会得到与实 

际情况不太吻合的排序结果. 

2.5不确定性 

污水处理工艺优化涉及的因素多且复杂,相当 

部分评价因素不能定量描述,具有不确定性(模糊 

性、灰色性、未确知性).同时处理系统各种行为的输 

人、输出也不是确定的数值.如果在决策过程中无视 

上述不确定性,将会损失大量的信息,导致决策的不 

准确和不科学.此外,各指标的影响程度也不相同, 

其权重的确定存在较大的主观不确定性. 

综上所述,污水处理工艺优化的过程是定性与 

定量、确定性与多种不确定性同时存在的多目标决 

策问题.而集对分析理论用联系度来统一处理模糊、 

随机、信息不完全所导致的综合不确定性,并作为一 

个确定一不确定系统来研究.这是将集对分析理论 

引人污水处理工艺优化决策过程的原因. 

维普资讯

72 湖南大学学报(自然科学版) 

合评价水平的不变 ¨.而污水处理工艺优选涉及的 

3现有优化策略存在的主要问题 

利用集对分析中的重要概念——联系度,来构 

建工艺设计优化的效用函数.根据不同的定义,分别 

形成两种不同的优化策略. 

优化策略一 :根据给定的各个方案,确定出一 

个理想方案,将各待选方案与理想方案构筑一个集 

指标,很多是相互影响、相互联系的,因此采用各待 

选方案指标值与理想方案指标值同一度的线性加 

权,显然不能反映实际. 

2)指标的权重是相互影响的,指标权重的赋值 

与指标的量化相比具有更大的不确定性和随意性. 

目前普遍采用的赋权方法是AHP法.AHP法是一 

种定性与定量相结合的决策分析技术,将问题分解 

成若干层次,由专家和决策者对所列指标通过两两 

对,计算备选方案各指标与理想方案各对应指标的 

同一度、差异度和对立度,并计算集对的联系度,寻 

找与理想方案最贴近的那个方案;或者单一地利用 

联系度中的同一度分量,结合各指标的权重计算待 

选方案与理想方案的加权同一度,根据所得的同一 

度数值大小确定方案的优劣次序. 

优化策略二 们:将待评价的评价指标量与评价 

标准构筑一个集对,对于一个具体决策来说,设有 

N个评价指标,其中有s,F。,F:,P个指标分别处 

于I,Ⅱ,Ⅲ,Ⅳ级标准,则该决策的联系度表达式为 

= +

争。+ :+ . 。,f:分别表示差异不确 

C 

定度, 表示对立度标记.设a= ,b。= ,b:= 

,c=丽/-

n 

则a,b (,z=1,2),c依次表示为同一 

度、差异不确定度、对立度.原联系度表达式可简写 

为 :a+b。 。+b2 i2+ .根据集对分析理论,同 

度、对立度是相对确定的,而差异度则相对不确 

定,a,b,c三者是对同一问题不同侧面的全面刻 

画,依据a,b,c三者大小关系及定量分析,可判断 

实际决策目标状况.a值越大,b,c值越小,工艺级 

别越优先.另外,根据a,b,c的具体数值,可对污水 

处理厂工艺状况进行分级,即以评价因子的数量指 

标相对于决策评价标准的达标、超标及其所占比例, 

确定污水处理厂最佳工艺.但是在实际工作中,某项 

评价指标的评价标准没有确切的定义,由于地区差 

异性,很难实现统一. 

总的来说,现有的应用主要存在以下不足: 

1)线性加权对于各评价指标相互联系的污水处 

理工艺方案优化不合理.对于文献中的对各方案的 

同一度进行线性加权,一般只适用于各评价指标相 

互独立的场合,若各评价指标间不独立,“和”的结果 

必然是信息的重复,难以反映客观实际;特别是线性 

加权法可使各评价指标间得到线性地补偿,任一指 

标值的减少,都可以用另一些指标的增量来维持综 

比较重要程度而逐层进行判断评分,利用计算判断 

矩阵的特征向量确定下层指标对上层指标的贡献程 

度,从而得到最底层指标对总目标而言重要性的赋 

权结果 .但AHP确定权重存在较大的主观随意 

性,当采用专家咨询方式时,容易产生循环而不满足 

传递性公理,导致标度把握不准并丢失部分信 

. 

3)我们认为联系度表达式原则上是一个结构函 

数,直接将联系度作为一确定的数值来进行各种运 

算是值得商榷的.若只单一的利用联系度中的同一 

度分量,受理想点的限制,尚未考虑对立度,具有一 

定的主观性和片面性,实质为一种理想点法. 

4基于集对分析的污水处理工艺优化多属 

性决策模型 

多属性决策问题一般可记为G=(S,E,W, 

D).方案集S={S。,S2,…,S },指标集E={e。,ez, 

,e },指标权重集W={ 。, :,…, },决策矩 

阵D:(d ,) ,dh为方案S 关于指标e 的属性 

值.决策矩阵经归一化处理,所得的归一化矩阵记为 

,, 

P=(P ,) .对效益型指标,有Ph=盖,d = Ⅱ

一 

max{d。,,d:,,…,d };对成本型指标,有Pk,= 

dml- 

d =min{dl ,d2,,…,d ,}. 

Ⅱ 

4.1熵权法确定指标权重 

指标权重以数量形式反映各指标在综合评价中 

的重要程度.若同一指标下各方案的指标值相差越 

大,则反映该指标在综合评价中比较重要.若某项指 

标值全部相等,则该指标在综合评价中不起作用.信 

息论中,信息熵是信息量的度量,其基本出发点是把 

获得的信息看作用以消除不确定性的东西.故信息 

量的大小可用被消除的不确定性的多少来表示.某 

项指标值变异程度越大,该指标提供的信息量越大, 

则该指标的权重也越大,反之亦然.因此可以根据各 

维普资讯

第11期 蒋茹等:基于集对分析的污水处理工艺设计优化 73 

指标提供的信息量大小来确定指标权重. 

利用熵权和普遍应用的AHP法确定权重,其出 

发点和理论基础均不同,前者从各指标提供的信息 

量多少出发,依据的是信息论;后者是由决策者分析 

6 ==-——cz ,——c ,=: — ! — :: ! ; 

(4) 

因此,对于效益型指标,集对{P “ }的联系度为 

评价,人为主观色彩较重,理论基础是决策论.本文 

采用主观与客观相结合的方法,在AHP法基础上, 

运用熵值理论对各专家评估水平赋予权重,确定各 

指标的组合权重.科学的确定权重方法应采用两者 

的算术平均值_1引. 

为确定各指标的权重,先计算各指标的熵值 

E,=一(Inm) ∑P lnPh (1) 

式中只,=Pk /∑/ ,若 =o'贝IJ ,=0. 

其次,计算各指标的偏差度D,=1一E,和客观 

加权因子 ,=D,/∑Di.设用AHP法得到的各指 

标的权重为 (r=1,2,…,7z),可用 ,进一步进行 

修正,即组合权重为 

硼 =( ,+ )/2 (2) 

4.2 基于集对分析的多属性决策模型n 

设最优方案和最劣方案分别记为U,V,它们 

可根据系统目标和客观条件,从方案集的内部或外 

部选定.U,V的选定涉及到指标类型,记诸方案关 

于P 的比较区间为[ ,,“,].最优方案和最劣方案 

关于P,的属性值分别为 ,, ,.当 属于效益型指 

标时,有“,=“,, = ,;当P 属于成本型指标 

时,有一U,=73,,一73 =U,.最优方案和最劣方案构成 

了诸方案的比较区间[V,U]. 

1)在P,的比较区间中确定集对{P }的联 

系度. 

对效益型指标, r“ 十 ,_73 一,U南十 一73 分别表示在  

数值上P 与“,,P 与 ,的接近程度.其值越大,接 

近程度越大,在确定区间[ ,“,]内,P 与 的接 

近恰说明了Pb与“,的远离.故定义集对{P “,} 

的同一度和对立度分别为 

户h 

一 

“h一

“ + 

’ 

p “, , ,^、 

h 

j 

由归一化条件a+b+c=1,定义差异度为 

P^, .(“,一Pb)(P 一 ).

,,

} _ 十——— _ -十 

. 

“,+ 

(5)

)户^ 、J 

 

对成本型指标,可确定集对{P ,}的联系度为 

“, , 

(“,一Pb)(户 一 ,)..

ult,

I 十— 十 

 

k ・ 

鸟 (6)o  

式(5),(6)中的同一度和对立度分别说明了方案 

关于指标e,的趋优程度和趋劣程度.两式关于同一 

度和对立度互为对称,且具有相同的差异度,说明了 

刻划的对立性和一致性.当P =“,或 ,时,差异度 

最小为零,表示刻划的不确定性最小;当Pb= 

时,差异度最大为2 /(“,+ ,),表示 

刻划的不确定性最大(但不超过0.5). 

将以上两式统一记为 

{d.. }=ab+6h +c^ (7) 

2)在比较区间[V,U]中计算集对{ ,U}的 

联系度 

{ .u}=a^+bk +c (8) 

式中 =∑ n打,bk=∑ , =∑ ,. 

r=l r=1 ,=l 

3) 与最优方案U的贴近度 

=— 一

. 

(9) 

a^1.C^ 

显然, 最大所对应的方案为最优方案,此时 

得到的排序为基序. 

4.3 稳定性分析 

当 >0,对a 正向修正,表示对方案 接近理 

想方案U持肯定意见,且 越接近1,修正作用越 

强.反之,当i<0,对 正向修正,表示对方案 接 

近理想方案U持否定意见,且 越接近一1,修正作 

用越强 引.因此,可通过分析 在[一1,1]内的变 

化,对方案排序的稳定性进行分析.为清楚表达,记 

i为△ ,则△ ∈[一1,1].假定对于整个方案集,△ 

相同.设 基序名次为1,记原 为r ;当△ ≠0 

时,其值变为 .设基序中有r:<r ,则根据式(9) 

可得 

维普资讯

74 湖南大学学报(自然科学版) 

cpa^一cka >0 (10) 

表2 各工艺设计方案的联系度及相对贴近度 

Tab.2 Weights of attributes by AItP and entropy method 

1)当0≤Ai≤1时,式(9)变为 

一(b 

n^+kAi)+ 

: 

(11)

为维持基序r:<r ,需满足不等式Cpa 一cka > 

(ckbp—cp )Ai.当ckbp一 6 ≤0时,有△ ∈[0, 

1];当ckbp—cp ≤0时,有0≤Ai≤ 

min{ , }. 

2)当一1≤△ <0时,式(9)变为 

一n^+(毛 bc^一 ^  )Ai.’ (12) 

同理,为维持基序 h<r ,有当akb 一apb ≥ 

0时,有△ ∈[一1,0);当akb 一apbl<0时,有 

max

【-1, }≤Ai<o0 .

取以上两种情况的并集即可得到基序排序 

<r:的稳定区域.以上进行的是局部分析,为了保 

证全序的稳定性,只需依次两两比较分析确定△ 的 

范围,则所有△ 取值区间的交集即为排序结果的稳 

定区域. 

4.4 实例研究 

拟建一座1O万m3/d的污水处理厂,二级生化 

处理,并要求脱氮除磷.有4个方案: /o法(S。)、三 

沟式氧化沟(S )、厌氧一单沟式氧化沟(S。)和SBR 

法(S ).所选定的评价指标共8个,分别为项目投 

资、经营成本、占地面积、氮磷去除效果、污泥处理效 

果、运行稳定性、工艺成熟性以及操作的难易程度. 

各方案的具体指标值可参见文献[3]. 

各指标权重、各方案的联系度及相对贴近度计 

算结果分别见表1和表2. 

表1 各指标的权重值 

Tab.1 Detailed infomat|0n 0f attribⅡtes 

根据表2中各方案的相对贴近度rJ,得到在相 

对确定条件下的排序结果为r >r >r。>r .因 

相对贴近度越大,其所对应的方案越优,所以方案优 

劣顺利分别为S,,S:,S ,S ,即厌氧一单沟式氧化沟、 

三沟式氧化沟、 /o法、SBR法,此为基序.该排序 

结果与文献[3]所采用的多层次灰色关联法得到的 

排序结果基本一致,只是在第2和第3位序上存在 

偏差. 

为了保证本文所得结果的稳定性,以基序为基 

础进行稳定性分析,其结果见表3.排序相邻方案的 

稳定性定量分析计算表明,得该排序结果的稳定区 

域为△ ∈[一1,1],这与i本身的定义区域是一致 

的.也就是说,不存在扩展序,通过本方法得到的排 

序结果是稳定的,其排序结果优于文献[3]的结果. 

表3排序稳定性分析 

Tab.3 Connection degrees and relative closeness 

degrees of each alternative 

5 结 论 

本文引入集对分析理论解决污水处理工艺方案 

优化的多属性决策问题,概念清晰,计算简便,结果 

稳定,便于编程实现.同时,为了减少决策者因受专 

业知识的限制和偏好等主观因素的影响,充分利用 

方案所反映的固有信息,在确定权重时引入熵权修 

正层次分析法主观权重的思想,综合考虑主、客观因 

素的影响,使得决策结果更趋向合理和科学. 

参考文献 

[1] ELLIS K V,TANG S L.Wastewater treatment optimization 

维普资讯

第11期 蒋茹等:基于集对分析的污水处理工艺设计优化 75 

model for developing world.I:Model development[J].Journal 

of Environmental Engineering Division,ASCE,1991,117:501 

581. 

[2]ELLIS K V,TANG S L.Wastewater treatment optimization 

model for development world.II:Model testing[J].Journal of 

Environmentla Engineering Division,AsCE,1994,120:610— 

624. 

[3]ZENG G M,JIANG R,HUANG G H,et a1.Optimization of 

wastewater treatment laternatives selection by hierarchy grey re・ 

lational analysis[J].Journal of Environmentla Management, 

2007,82:250—259. 

[4]王玲玲,曾光明,黄国和,等.小区域污水处理的多目标规划设 

计探讨[J].环境工程,2004,22(6):88—89. 

WANG LL,ZENGGM,HUANGGH,eta1.Study onmulti— 

objective programming for wastewater treatment in small region 

[J].Environmentla Engineering,2004,22(6):88—89.(In Chi・ 

neae) 

[5]凌猛,杭世王君.城市污水处理厂工艺方案模糊决策方法的应 

用[J].给水排水,1998,24(3):6—9. 

LIN M,HANG S J.The application of fuzzy decision—making 

approach in the process altemativea of municipal wastewater 

traetment plnat[J].Water and Wastewater Engineering,1998, 

24(3):6~9.(In Chinese) 

[6]赵克勤,宣爱理.集对论——一种新的不确定性理论、方法与应 

用[J].系统工程,1996,14(1):18—23. 

ZHAO K Q,XUAN A L.Application of lf novel uncertainty theo— 

ry nad method——Set pair analysis[J].Systems Engineering, 

1996,14(1):18—23.(InChinese) 

[7]JIANGY L,XUCF,YAOY,et a1.A new approachfor repr. 

aeenting and processing uncertainty knowledge[C]//The 2003 

IEEE Intemational Conference on Information Reuse and Inte. 

gration.USA.2003:466—470. 

[8]ANNELIES J B,HEINZ A P,RALF 0,et a1.Indiactors for 

the sustainabiilty assessment of wastewater treatment systems 

[J].Urban Water,2002(4):153—161. 

[9]王国平,王洪光.集对分析用于污水处理厂的综合评价[J].江 

苏环境科技,2002,15(1):16—18. 

WANG G P,WANG H G.Comprehensive assessment in sewage 

treatment plant by integration and contrast analysis[J].Jiangsu 

Environmental Science and Technology,2002,15(1):16—18. 

(In Chinese) 

[1O]李凡修,梅平,陈武.多元集对模型在污水处理厂改造决策中应 

用[J].环境科学与技术,2004,27(6):92—94. 

LI FX,MEI P,CHENW.Set pairanalysismethodwithmulti・ 

factor connection number for plnaning of WTP reformation[J]. 

EnvironmentlaScience and Technology,2004,27(6):92—94. 

(In Chinese) 

[11]李香云,刘理才.一种集成的科技实力综合评价方法[J].科学 

学研究,2004,22(6):85—87. 

LI X Y.LIU L C.A literature review on the study and evalu・ 

ation of scientific and technological strength[J].Studies in cSience 

of cSience,2004,22(6):85—87.(In Chinese) 

[12]张佳林,尉晓君.基于AHP法的地下水体综合评价[J].湖南 

大学学报:自然科学版,2006,33(5):136—140. 

ZHANG J L,WEI X J.Comprehensive assessment of groundwa— 

tel"based on AHPEJ].Journal of Hunan University:Natural cSi. 

encea,2006,33(5):136—140.(InChinese) 

[13]方创琳,毛汉英.区域发展规划指标体系建立方法探讨[J].地 

理学报,1999,54(5):410—419. 

FANG C L.MA0 H Y.A system of indicators for regional de— 

velopment planning[J].Acta Geographiac Siniac,1999,54(5): 

410—419.(In Chinese) 

[14]张文泉,张世英,江立勤.基于熵的决策评价模型及应用[J]. 

系统工程学报,1995,10(3):69—75. 

ZHANG W Q,ZHANG S Y,JIANG L Q.Application of deci・ 

sion—making model based on entropy[J].Journal of Systems 

Engineering,1995,10(3):69—75.(InChinese) 

[15]张斌.多目标系统决策的模糊集对分析方法[J].系统工程理 

论与实践,1997(12):108—114. 

ZHANGB.Fuzzyset pairanalysisformulti.objcetive systemdeci— 

sion.making[J].Systems Engineering—Theory&Practice,1997 

(12):108—114.(InChinese) 

发布评论

评论列表 (0)

  1. 暂无评论