《1 前言》
1 前言
在水环境质量评价中, 文献[1]首次采用了可拓评价方法。在该方法中, 待评物元的量值为点, 水质标准的量值为区间, 评价方法建立在点与区间的距离——距的概念基础上。该方法也同时应用于围岩稳定性分类[2]。但在实际应用中, 由于测量、计算所带来的数据误差、信息不完全带来的数据不完全、人为主观因素带来的数据不准确等原因, 待评物元的量值可能为一区间, 这时要考虑区间与区间的距离, 并在区间与区间之距——IR距概念[3]的基础上, 给出相应的可拓评价方法。
《2 区间距》
2 区间距
在实域上起作用是点和区间, 而点x与点y的距离在经典数学中已有定义, 即
ρ(x,y)=|x−y|。
为了建立关联函数, 还要定义点与区间距离——距的概念[4,5,6,7]。
经典数学中的有限开区间 (a, b) 、闭区间[a, b]、半开半闭区间 (a, b]和[a, b) 统称为区间, 用符号〈a, b〉表示, 也就是说, 区间〈a, b〉可以包含a (或b) , 也可以不包含a (或b) 。设x为实域 (-∞, +∞) 上的任一点, 〈a, b〉为实域上任一有限区间, 则称
ρ(x,⟨a,b⟩)=∣∣x−a+b2∣∣−12(b−a)
此处距的概念与经典数学中距离稍有不同。当 〈a, b〉中不含x时, ρ (x, 〈a, b〉) 与经典数学中点与区间距离d概念相同, 即: ρ (x, 〈a, b〉) = d, d是离x最近区间端点与x的距离。当x点在〈a, b〉之内时, 经典数学中认为点与区间距离d = 0;而这里的距为负值, 负值的不同表示x在区间〈a, b〉内位置不同。例如, 设给定区间〈3, 5〉, 则ρ (4, 〈3, 5〉) =|4-4|-1= -1。
为了建立区间上的关联函数, 文献[3]引入了区间与区间之距——IR距和区间关于区间的位值——IR位值的概念:
设 〈x, y〉 与 〈a, b〉 为实数域上的任意有限区间, 则称
ρI(⟨x,y⟩,⟨a,b⟩)=12(ρ(x,⟨a,b⟩)+ρ(y,⟨a,b⟩))
为区间 〈x, y〉 与区间〈a, b〉的区间距, 简称IR距。当〈x, y〉为一单点集时, IR距就是点与区间之距。为了区别, 点x与区间〈a, b〉之距用ρ (x, 〈a, b〉) 表示, 区间〈x, y〉 与区间〈a, b〉的区间距用ρI (〈x, y〉, 〈a, b〉) 表示。
设X1=〈a, b〉, X2=〈c, d〉且X1⊂ X2, 则
DI(⟨x,y⟩,X1,X2)={−1,[x,y]⊂[a,b]ρI(⟨x,y⟩,X2)−ρI(⟨x,y⟩,X1),其它
称为区间〈x, y〉关于区间 X1 , X2的区间位值, 简称IR位值。当 〈x, y〉为一单点集时, IR位值就是点关于区间的位值。为了区别, 点关于区间X1 , X2的位值用D (x, X1 , X2) 表示, 区间〈x, y〉关于区间X1 , X2的位区间位值用DI (〈x, y〉, X1 , X2) 表示。
关于区间论域上的可拓集理论以及建立在IR距上的关联函数理论参见文献[3]。
《3 区间可拓评价方法》
3 区间可拓评价方法
R1=(N1,C,V1)=⎡⎣⎢⎢⎢⎢⎢N1C1C2⋮Cn⟨a11,b11⟩⟨a21,b21⟩⋮⟨an1,bn1⟩⎤⎦⎥⎥⎥⎥⎥‚R2=(N2,C,V2)=⎡⎣⎢⎢⎢⎢⎢N2C1C2⋮Cn⟨a12,b12⟩⟨a22,b22⟩⋮⟨an2,bn2⟩⎤⎦⎥⎥⎥⎥⎥‚⋯⋯Rm=(Nm,C,Vm)=⎡⎣⎢⎢⎢⎢⎢NmC1C2⋮Cn⟨a1m,b1m⟩⟨a2m,b2m⟩⋮⟨anm,bnm⟩⎤⎦⎥⎥⎥⎥⎥
为m个同征 (C1, C2, …, Cn) 物元[6,7], 则称
R0=[NCN1V1N2V2⋯⋯NmVm]=⎡⎣⎢⎢⎢⎢⎢⎢⎢NC1C2⋮CnN1⟨a11,b11⟩⟨a21,b21⟩⋮⟨an1,bn1⟩N2⟨a12,b12⟩⟨a22,b22⟩⋮⟨an2,bn2⟩⋯⋯⋯⋯N⟨a1m,b1m⟩⟨a2m,b2m⟩⋮⟨amm,bnm⟩⎤⎦⎥⎥⎥⎥⎥⎥⎥
是同征物元R1, R2, …, Rm的同征物元体, 其中Nj表示所划分的第j个评价类别, Ci表示第i个评价指标, Vij = 〈aij, bij〉分别为Nj关于指标Ci所规定的量值范围, 即各类别关于对应的评价指标所取的数据范围经典域。令
RP=(P,C,VP)=⎡⎣⎢⎢⎢⎢P,C1,C2,⋮Cn,V1PV2P⋮VnP⎤⎦⎥⎥⎥⎥=⎡⎣⎢⎢⎢⎢⎢P‚C1,C2,⋮Cn,⟨a1P,b1P⟩⟨a2P,b2P⟩⋮⟨anP,bnP⟩⎤⎦⎥⎥⎥⎥⎥‚
其中, P表示类别的全体, ViP = 〈aiP, biP〉为P关于Ci 所取的量值的范围, 即P的节域, 且Vij⊂ViP (i = 1, 2, …, n, j = 1, 2, …, m) 。
步2 确定待评物元 对待评的事物q, 把所检测得到的数据或分析结果用物元
Rq=(q,C,Vq)=⎡⎣⎢⎢⎢⎢⎢q,C1,C2,⋮Cn,V1qV2q⋮Vnq⎤⎦⎥⎥⎥⎥⎥=⎡⎣⎢⎢⎢⎢⎢q‚C1,C2,⋮Cn,⟨c1q,d1q⟩⟨c2q,d2q⟩⋮⟨cnq,dnq⟩⎤⎦⎥⎥⎥⎥⎥
表示, 称为事物q的待评物元, 式中q表示某事物, 〈ciq, diq〉为q关于Ci的量值, 即待评事物检测所得的具体数据。
步3 确定权系数 确定指标Ci的权系数为αi, 且
Kj(Viq)=ρI(⟨ciq,diq⟩,⟨aij,bij⟩)DI(⟨ciq,diq⟩,⟨aij,bij⟩,⟨aiP,biP⟩)
其中Kj (Viq) 为待评事物q对于评价指标Ci的量值〈ciq, diq〉关于等级j的关联度。
Kj(q)=∑i=1nαiKj(Viq)。
步6 等级评定 若kj0(q)=maxj‰∈{1,2,⋯,m}Kj(q), 则评定q属于等级j0。令
K¯¯¯j(q)=Kj(q)−minjKj(q)maxjKj(q)−minjKj(q),j∗=∑j=1mjK¯¯¯j(q)/∑j=1mK¯¯¯j(q),
《4 确定权重的可拓方法》
4 确定权重的可拓方法
设Viq=⟨ciq,diq⟩‚Vij=⟨aij,bij⟩,rij(Viq,Vij)=12(rij(ciq,Vij)+rij(diq,Vij)),其中 rij(x,Vij)=⎧⎩⎨⎪⎪⎪⎪⎪⎪⎪⎪2(x−aij)bij−aij, x≤aij+bij22(bij−x)bij−aij, x≥aij+bij2i=1,2,⋯‚n, j=1,2,⋯,m, (1)则rijmax(Viq,Vijmax)=maxj{rij(Viq,Vij)}。 (2)
如果指标i的数据落入的类别越大, 该指标应赋以越大的权重, 则取
ri=⎧⎩⎨⎪⎪⎪⎪⎪⎪⎪⎪jmax×(1+rijmax(Viq,Vijmax)), 当rijmax(Viq,Vijmax)≥−0.5时;jmax×0.5, 当rijmax(Viq,Vijmax)<−0.5时。 (3)
如果指标i的数据落入的类别越大, 该指标应赋以越小的权重, 则取
ri=⎧⎩⎨⎪⎪⎪⎪⎪⎪⎪⎪(m−jmax+1)×(1+rijmax(Viq,Vijmax)), 当rijmax(Viq,Vijmax)≥−0.5时;(m−jmax+1)×0.5, 当rijmax(Viq,Vijmax)<−0.5时。 (4)
αi=ri/∑i=1nri。 (5)
《5 实例计算》
5 实例计算
为了便于比较, 本文仍引用文献[1,8]中的同一实例, 即汾河上游、中游、下游3个断面的水质评价。文献[8]用的是灰关联评价方法。文献[1]用的是可拓评价方法。
文献[8]将水质标准取为一个定值, 带有一定的局限性, 也不符合标准实际。文献[1]弥补了这一缺陷, 将水质标准取为一个区间。汾河上游的兰村、中游的石滩和下游的河津3个断面设立了水质监测站, 1981—1984年每年分别在丰水期、平水期、枯水期取样监测。文献[1,8]中, 待评物元的量值取为各断面实测元素浓度的平均值, 该平均值只能反映监测数据的一个数字特征, 就像水质标准取为一个定值一样, 不能很好地反映其实际问题。本文将待评物元的量值取为各断面实测元素浓度的区间值。
R=⎡⎣⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢qDOBOD酚CN−HgAsCr+6兰村⟨9.1‚ 9.5⟩⟨1‚ 1.5⟩⟨0.0008‚ 0.003⟩⟨0.0009‚ 0.0015⟩⟨0.0007‚ 0.00089⟩⟨0.002‚ 0.004⟩⟨0.0055‚ 0.0068⟩石滩⟨8.2‚ 8.8⟩⟨3.9‚ 4.45⟩⟨0.006‚ 0.008⟩⟨0.002‚ 0.004⟩⟨0.0005‚ 0.0007⟩⟨0.068‚ 0.076⟩⟨0.01‚ 0.02⟩河津⟨6.4‚ 6.7⟩⟨7‚ 7.5⟩⟨0.01‚ 0.02⟩⟨0.008‚ 0.009⟩⟨0.00045‚ 0.00055⟩⟨0.038‚ 0.045⟩⟨0.00065‚ 0.00078⟩⎤⎦⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥
R0=⎡⎣⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢NDOBOD酚CN−HgAsCr+61级⟨9‚ 10⟩⟨1‚ 3⟩⟨0‚ 0.001⟩⟨0‚ 0.005⟩⟨0‚ 0.00005⟩⟨0‚ 0.05⟩⟨0‚ 0.01⟩2级⟨6‚ 9⟩⟨3‚ 4⟩⟨0.001‚ 0.002⟩⟨0.005‚ 0.05⟩⟨0‚ 0.00005⟩⟨0‚ 0.05⟩⟨0.01‚ 0.05⟩3级⟨5‚ 6⟩⟨4‚ 5⟩⟨0.002‚ 0.005⟩⟨0.05‚ 0.2⟩⟨0.00005‚ 0.0001⟩⟨0‚ 0.05⟩⟨0.01‚ 0.05⟩‚⟨0.01‚⋯0.05⟩4级⟨3‚ 5⟩⟨5‚ 6⟩⟨0.005‚ 0.01⟩⟨0.05‚ 0.2⟩⟨0.0001‚ 0.001⟩⟨0.05‚ 0.1⟩⟨0.05‚ 0.1⟩5级⟨2‚ 3⟩⟨6‚ 10⟩⟨0.01‚ 0.1⟩⟨0.05‚ 0.2⟩⟨0.001‚ 0.001⟩⟨0.05‚ 0.1⟩⎤⎦⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥
RP=(P,C,VP)=⎡⎣⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢PDOBOD酚CN−HgAsCr+6⟨1‚ 10⟩⟨0‚ 10⟩⟨0‚ 0.1⟩⟨0‚ 0.2⟩⟨0‚ 0.001⟩⟨0‚ 0.1⟩⟨0‚ 0.1⟩⎤⎦⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥。
按式 (1) , (2) , (4) , (5) 计算其权重为:
W7×3=⎡⎣⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢0.200.160.070.220.040.090.230.200.150.120.310.070.070.090.190.060.040.260.070.140.25⎤⎦⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥。
按评价方法步3、步4、步5中的计算公式计算, 其评价结果见表1。
表1 区域水环境质量的区间可拓评价结果
Table 1 The interval extension assessment result of water environmental quality
《表1》
Kj (q) | N1 | N2 | N3 | N4 | N5 | max | j0 | j* |
兰村 | 0.050 | -0.452 | -0.608 | -0.718 | -0.858 | 0.050 | 1 | 1.78 |
石滩 | -0.228 | -0.147 | -0.498 | -0.499 | -0.588 | -0.147 | 2 | 2.01 |
河津 | -0.239 | -0.133 | -0.472 | -0.485 | -0.496 | -0.133 | 2 | 1.68 |
表1中:N1数据分别表示汾河上游、中游、下游3个断面水质属于1级的关联度;同理, N2, N3, N4, N5下面的数据分别表示3个断面水质属于对应2级、3级、4级与5级的关联度;j0数据分别表示3个断面水质最大关联度所对应的级别;j*数据分别表示3个断面水质的级别变量特征值。从表1可以看出, 兰村水质属1级 (j0 = 1) 偏2级 (j* = 1.78) , 石滩水质属2级 (j0=2, j*=2.01) , 河津水质属2级 (j0=2, j* =1.68) 。
文献[1]采用可拓评价方法, 评价结果为:兰村水质属1级 (j0 = 1) 偏2级 (j* = 1.76) , 石滩水质属2级 (j0 = 2, j* = 2.01) , 河津水质属2级 (j0 = 2, j* = 1.84) 。
文献[8]采用灰关联评价方法, 评价结果为:兰村水质属1级 (关联度为0.862) , 石滩水质属2级 (关联度为0.538) , 河津水质属2级 (关联度为0.627) 。
从以上可以看出, 本文区间可拓评价方法与文献[1,8]评价的结果是一致的。
《6 结语》
6 结语
解决了水质等级标度为区间同时实测水质也为区间时进行水质评价的问题。所提供的方法简单、易行、合理, 具有一定的实际意义, 并对待评物元的量值为区间的问题都适用, 具有普适性。
所给的权重公式是用简单关联函数构造的。该公式计算简单, 易于应用。不过, 从rij计算出ri时要根据实际情况选择合理公式。