水资源系统中的不确定性及风险分析方法X
第26卷 第2期Vol.26 No.2干旱区地理 2003年6月Jun. 2003ARID LAND GEOGRAPHY
文章编号:1000-6060(2003)02-0116-06
水资源系统中的不确定性及风险分析方法
左其亭1,2 吴泽宁1 赵 伟3
(1郑州大学环境与水利学院,郑州 450002; 2中国科学院新疆生态与地理研究所,乌鲁木齐 830011;
3新疆水利规划设计局,乌鲁木齐 830001)
Ξ
摘 要:水资源系统中广泛存在着不确定性,其对系统的影响很大,甚至带来灾害性风险(如洪水、干旱),是风险产生的根本原因,也是水资源系统研究遇到的难点问题之一。(包括随机性、模糊性、灰色性及未确知性)及研究方法,处理方法入手,介绍基于概率统计学(处理随机性)的风险计算模型,()、基于灰数
(处理灰性)和基于未确知数(处理未确知性)。
关键词:水资源系统 不确定性 中图分类号:TV213.91 引言
水,是生命之源,是人类及一切生物赖以生存的一种不可缺少的宝贵资源。从水的基本属性来看,水资源系统是客观存在于自然界的天然系统。然而,它与人类社会经济发展有着千丝万缕的联系。人类发展依赖水资源,又对水资源有越来越大的影响,特别是随着社会经济的发展,出现的水问题(如干旱、洪水、水污染)日益突出。
随着社会实践的需要和研究方法的发展,人们已经看到了把水资源看成一个系统来研究的重要性。在人类社会早期,人们只能通过简单的劳动(如堤防、沟渠),来改善水资源利用;随着社会的进一步发展,对水资源的需求量增加,开发利用水资源的力度加大,综合开发利用水资源的工程措施和水资源系统分析与管理的思想逐步形成;到了20世纪中期,由于经济迅速发展和人口剧烈膨胀,水资源供需紧张和生态环境恶化的局面越来越严峻。因此,需要把水资源纳入到一个系统,来研究合理开发、利用和优化配置水资源问题。
水资源系统是一个十分复杂的不确定性系统,广泛存在着随机性、模糊性、灰色性、未确知性。正是由于不确定性的广泛存在,水资源系统在开发、利
Ξ收稿日期:2002-11-06; 修回日期:2003-02-28
用及获得经济效益的同时,也存在一定的风险(如干
旱、洪水、水污染、经济亏损、生态环境破坏)。那么,如何评价和降低风险?如何进行水资源系统风险规划与管理?解决这些问题的基础是要定量研究不确定性问题以及寻找风险分析计算方法。
2 水资源不确定性系统
2.1 客观存在的不确定性类型
1,2〕
到目前为止,人们已认识到四种不确定性〔:(1)由于条件提供的不充分和偶然因素的干
扰
,使几种人们已经知道的确定结果的出现呈现偶然性,但在某次试验中不能预料哪一个结果发生。这种不确定性即为“随机性”;
(2)由于事物的复杂性,事物的界线不分明,对其概念不能给出确定的描述和确切的评定标准。这种不确定性即为“模糊性”;
(3)由于事物的复杂性、信道上噪音干扰和接收系统能力(含人的辩识能力)的限制,人们只知系统的部分信息或信息量所呈现的大致范围。这种部分已知、部分未知的不确定性即称为“灰色性”;(4)纯主观上的、认识上的不确定性称为“未确知性”。与灰色性相比,它具有较多的信息量,不但
基金项目:国家重点基础研究发展规划项目(G1999043608)及中国科学院知识创新工程重要方向项目(KZCX3-SW-327)
作者简介:左其亭,1967年生,男,河南固始县人,工学博士。主要研究方向为水文水环境与水资源管理。已发表论文66篇,主持和参
加课题12项。E-mail:[email protected]
2期 左其亭等:水资源系统中的不确定性及风险分析方法 117
知道信息量的取值范围,还知道所求量在该区间的分布状态。
处理各种不确定性已有了各自的数学方法。处理随机性的数学方法是随机理论与概率统计;处理模糊性的数学方法是模糊数学;处理灰色性的数学方法是灰色数学;处理未确知性的数学方法是未确知数学。但在实践中常常遇到在同一个系统中几种不确定性同时出现或交叉出现的情况。因此,研究综合处理不确定性的数学方法──不确定性数学是非常必要的。不确定性数学是以概率统计、模糊数学、灰色数学、未确知数学为主体,四方面密切相关、互相渗透组成的数学体系〔1,3〕。2.2 水资源系统中的不确定性
从水资源系统的输入,性一面。个方面:;其二是系统内部结构本身存在的不确定性;其三是系统输出存在的不确定性。下面就来简单说明水资源系统中存在的四种不确定性。
(1)水资源系统中存在的随机性从系统的输入与输出关系来看,由于系统环境的变化可能会影响系统功能的变化,这种变化常常表现为随机性。比如,流域降雨径流量的大小会随着降雨的随机变化而变化,从而表现为随机性。从系统结构及内部状态来看,由于水结构的复杂性,使得人们对水文规律的认识和参数的获得常常带有随机性。比如,对地下水水文地质参数的试验,常常是随机地选择钻探点和抽水试验点,得到的水文地质参数当然也存在一定的随机性。
(2)水资源系统中存在的模糊性水资源系统是一个复杂的系统,在系统输入与输出、系统结构与状态等方面,很多有关的概念界限不分明。比如,我们常说的“洪水季节”与“枯水季节”“、稳定流”与“非稳定流”、“含水层”与“隔水层”等等,都是界限不分明的模糊概念。
(3)水资源系统中存在的灰色性从水质点的运移机理、流速状态以及介质结构来看,人们不可能对其完全清清楚楚,即永远是灰的。
从建立的系统模型来看,由于人类认识能力的限制、测试手段的限制,使得我们获取的资料或精度
不高或残缺不全。当然,在这种条件下建立的系统
模型也是灰的。
从系统的功能来看,由于对系统的结构、内部状态以及系统边界都不可能完全清楚,因此,对系统功能的了解和认识也不完全清楚和精确。比如,对降雨径流量的计算,由于种种原因,常常计算结果误差较大,即所得的结论也应该是灰的。
(4)水资源系统中存在的未确知性
由于水资源系统是一个庞大的复杂系统,对系、,一;,,,,导水通道又在何处等。虽,但如果要准确做出这一判断,必然要耗费很大的人力、物力。也就是说,我们总可以在一定程度上做得到。但实际中不需要我们完完全全了解,只需要通过某些手段较准确地估计或计算出隔水层的渗透量即可。这就是未确知性。
2.3 水资源不确定性系统的概念
从以上分析可以看出,水资源系统广泛存在着四种不确定性,即随机性、模糊性、灰色性、未确知性。作为一个客观的水资源系统,确切地说
,是一个未确知性的灰色—模糊—随机系统,称之为“水资源
〔4~8〕不确定性系统”。
在研究水资源系统时,不可避免地会遇到不确定性的处理问题。我们应该正确理解水资源系统中存在的不确定性,不能一遇到有关不确定性问题,就把它概化成简单的确定性问题,而应该对不确定性问题采取合理的处理方法,使之更切合实际。
3 水资源系统不确定性带来的风险问
题
由于客观世界的复杂性和人类认识的局限性,人类所有活动和决策都不可避免地伴随着不确定性的影响,因而也不可避免地冒一定风险。如洪水灾害、干旱灾害、地震灾害、台风灾害、战争以及股市下跌、经济突然滑坡等等。
关于风险的定义有许多种,可以从不同的角度来对风险进行定义。比如,从探讨风险与损失之间的联系入手,定义风险是产生损失的可能性;从某种
干 旱 区 地 理 26卷11 8
程度的损失产生的概率入手,定义风险是产生某种程度损失的机会。
无论是哪种定义都承认风险是由不确定性因素产生的。正是由于客观世界和人类活动中存在大量的不确定性因素,致使人们不能准确地预测未来事件发生的后果及其可能性大小。从这个意义上讲,所谓风险是指发生损失的不确定性。
风险具有三个基本特征:客观性、损失性和不确定性。“客观性”是指不论人们是否意识到,风险都是客观存在的,如地震、风暴、洪水等。“损失性”是指有可能带来损失的事件才称为有风险,即风险总是与损失相联系的,如洪水灾害、干旱灾害、地震灾害等。“不确定性”,,而,识别、评估不确定性因素带来的风险,并在此基础上优化组合各种风险管理技术,作出风险决策,包括水资源系统风险规划与管理,对风险实施有效的控制和妥善处理风险所致损失的后果,以期以最少的投入获得最大的安全保障〔9~11〕。
或 E(x)=
∫xf(x)dx
A
(4)
式中:xi,p(xi)—离散型风险变量及相应的概率;f(x)—连续型风险变量的密度函数。
标准差计算式:
2
σ=D(x)=(5) E(x-)
则,风险度(即变异系数)计算模型方程式:
σ
(6) FD=()
Ex
4.2 4 水资源系统风险分析计算模型
在风险分析中,对风险出现的可能性及大小有
两种度量方法:一种是基于风险出现的概率,采用“风险率”度量;一种是基于风险的变异性测度,采用“风险度”度量。4.1 基于概率统计的风险计算模型
(1)风险率计算根据事件A的性质,有两种确定方法:
(如洪水),一种是,如果事件A为“失事事件”
则事件A的风险率为:
(1) FP=P(A)
其中,FP为风险率,P(A)为事件A的概率。以下类同。
(如可供另一种是,如果事件A为“安全事件”
水量),则事件A的风险率计算式为:
(2) FP=1-P(A)
(2)风险度计算期望值计算式:
(3) E(x)=∑xip(xi)
x∈A
i
“,对“风
(),而是用〔0,()。
(设模糊事件A的模糊概率为P(A),则模糊风险率的计算也有两种情况:
(如洪一种是,如果模糊事件A为“失事事件”
水),则A的模糊风险率为:
(7) FP(A)=P(A)
(如另一种是,如果模糊事件A为“安全事件”
可供水量),则A的模糊风险率为:
(8) FP(A)=1-P(A)
其中,P(A)的计算式如下: P(A)=
μ(x)dp=∫
U
A
E(μA(x))
(9)
或
P(A)=ΣμA(xi)p(xi)
x∈A
i
(10)
其中:U—是模糊子集;μA(x)—是A的隶属度函数;p(xi)—是xi对应的概率。
(2)模糊风险度计算模糊事件A的期望值:
xμ(x)dp
E(A)==
P(A)
μ(x)dp∫
xμA(x)dp
A
A
(11)
U
或
E(A)=方差定义为:σ(A)=
2
x∈A
ΣxiμA(xi)p(xi)
p(A)
(12)
P(A)
2
∫
U
(x-E(A))2μA(x)dx
(13)
22
=E(A)-(E(A))
其中,E(A)=
P(A)
U
x2μA(x)dp
2期 左其亭等:水资源系统中的不确定性及风险分析方法 119
A的模糊风险度定义为:FD(A)=
2
()
∧
(如另一种是,如果灰色事件A为“安全事件”
∧
E(A)
安全供水量),则A的灰色风险率为
:
∧
∧
=
()
(x-E(A))μ(x)dx
xμ(x)dx2
A
FP(A)=1-P(A)
(3)灰色风险度计算
∧
(18)
A
灰色事件A的期望值:
xdp=
(x)+μ(x)
(x+(xP(A)
=
E(A)
2
2
μ(x)+μ(x)
-
(14)
∧
E(A)
4.3 基于灰数的风险计算模型
(1)灰集合与灰数的定义
所谓G是论域U上的一个灰子集,从U到闭区间〔0,1〕G:U→0,u)0〕和
μG:U→〔0,1〕,u|→μG(u)∈〔0,1〕μG分别称为G的上隶属其中,G≥μG;G与
函数和下隶属函数;G(u)与μG(u)分别称为元素
u相对于G的上隶属度和下隶属度。
∧
(19)
P(A)
∧
x∈A
∑
xi
∧
(x)()
p(xi)
2
∧
E(A)=(20)
P(A)
方差定义为:σ(A)=
2
∧
∧
P(A)
U∧
(x-E(A))2
∧
∧
(x)+μ(x)
dp
2
(21)
设论域U=R,则称灰集合(x)
μ(x),(x)∈〔 G, x∈R, 0,1〕
μ(x)
为灰数,并简记作G。
(2)灰色风险率计算
∧
∧
22
=E(A)-(E(A))
其中,E(A2)=
∧
∧
∧
P(A)
U
x
2
(x)+μ(x)
dp(22)
2
A的灰色风险度定义为:
∧
灰色事件A的灰色概率P(A)可以表达为:
∧
FD(A)=
2∧
∧
P(A)=
∫
U
=E(
∧
μ(x)+(x)
=dp
2
μ(x)+(x)
2
)
E(A)
(15)(16)
=
∧
μ(xi)+(xi)
或P(A)=p(xi)
2∧
x∈A
i
(x)+μ(x)x
dx
U
(x-E(A))2
∧
(x)+μ(x)
dx
2
∧
P(A)
=
2
∧
∧
其中,p(xi)—是xi对应的概率;其它符号同前。
∧
∧
2
∧
(23)
E(A)
设灰色事件A的灰色概率为P(A),则灰色风险率的计算也有两种情况:
∧
4.4 基于未确知数的风险计算模型
(1)一般未确知数的定义
(如洪一种是,如果灰色事件A为“失事事件”
∧
水灾害),则A的灰色风险率为
∧
∧
FP(A)=P(A)(17)
对于区间数〔a,b〕,若函数F(x)满足下列条
)内为单调非减函数、件:①F(x)在(-∞,+∞右
连续,至多有有限个间断点。②当x
干 旱 区 地 理 26卷12 0
则称〔a,b〕与F(x)构成一个未确知数,记作
{〔a,b〕,F(x)}〔a,b〕为取值区间,F(x)为区间〔a,b〕上的
主观可信度分布函数。
(2)未确知风险率计算
根据一般未确知数的定义可知,真值落到(x1,
x2)的可信度为F(x2)-F(x1)。这种可信度可以
参考文献〔References〕
〔1〕 王清印,倪天智,刘开第等.不确定性数学及其研究方向〔A〕.
灰色系统学术论文集〔C〕.开封:河南大学出版社,1993.151~
154.〔WangQing2yin,NiTian2zhi,LiuKai2dietal.Uncertaintymathematics.
ProceedingsofGreySystemStudy.
Kaifeng:
Henanuniversitypress.1993.151~〕
定义为未确知风险率,即
FP(A)=F(x2)-F(x1)
其中,A∈〔x1,x2〕
(3)未确知风险度计算
(24)
〔2〕 左其亭,.〕.:〔ZuoQi2ting,Ma
2ofgroundwatersystemandHydrogeologyandEngineeringGeolo2,21(5):41~45.〕
未确知事件A的期望值:
E(x)=ix,xi1((26)
3〕 刘开弟,吴和琴,王念鹏等.不确定性信息数学处理及应用
〔M〕.北京:科学出版社,1999,1~20.〔LiuKai2di,WuHe2qin,
WangNian2pengetal.mathematicalprocessofuncertaintyinfor2mationanditsapplications.Beijing:sciencepress.1999,1~20.〕
或 E(xφ1
x2
式中,xi,φ(xi)—离散型风险变量及相应的分φ(xi)=F(xi)-F(xi-1);φ(x)—布密度函数。
连续型风险变量的密度函数。
标准差计算式:σ=
D(x)=
E(x-x)
2
〔4〕 左其亭,王中根.现代水文学〔M〕.郑州:黄河水利出版社,
2002.17~20.〔ZuoQi2ting,WangZhong2gen.Modernhydrolo2gy.Zhengzhou:YellowRiverHydrologicalPress.2002.17~22.〕
(27)
〔5〕 夏军.灰色系统水文学〔M〕.武汉:华中理工大学出版社,1999.
1~23.〔XiaJun.Greysystemhydrology.Wuhan:Huazhongu2niversityofscienceandtechnologypress.1999.1~23.〕
则,风险度(即变异系数)计算模型方程式:
(28) FD=()
Ex
〔6〕 王清印,左其亭,王峰松等.灰色数学基础〔M〕.武汉:华中理工
大学出版社,1996.7~34.〔WangQing2yin,ZuoQi2ting,Wang
Feng2songetal.Basicgreymathematics.Wuhan:Huazhonguni2versityofscienceandtechnologypress.1996,7~34.〕
5 结束语
(1)本文仅探讨两方面内容:一是针对水资源
〔7〕 刘绍英,庞彦军,苑乐仁等.未确知数学及其应用〔M〕.保定:河
北大学出版社,1994.1~8.〔LiuShao2ying,PangYian2jun,
YuanYue2renetal.Unascertaintymathematicsanditsapplica2tions.Baoding:Hebeiuniversitypress.1994.1~8.〕
系统的特点,论证其存在的不确定性类型,包括模糊性、随机性、灰色性及未确知性,提出水资源不确定性系统的概念;二是从不确定性因素的数学处理方
法入手,介绍基于概率统计学(处理随机性)的风险计算模型,提出基于模糊隶属度(处理模糊性)、基于灰数(处理灰性)和基于未确知数(处理未确知性)的风险计算模型。
(2)因篇幅所限,本文对含有多种不确定性的复杂系统的风险计算模型没有进一步推导。关于这一点,可以参照不确定数学的研究方法,按照上述定义风险的思路可以进行推导,这里不再多叙。
(3)因篇幅所限,本文没有对各模型的应用列举实例。这不影响风险分析计算模型的应用。相关成果可参见文〔12〕或其它研究成果。
〔8〕 ZbigniewWKundzewicz.NewUncertaintyConceptsinHydrolo2
gyandWaterResources〔M〕.CambridgeUniversityPress.2001.1~23.
〔9〕 王丽萍,傅湘编著.洪灾风险及经济分析〔M〕.武汉:武汉水利
电力大学出版社,1999.1~71.〔WangLi2ping,FuXiang.
Floodriskandeconomicanalysis.Wuhan:Wuhanuniversityofhydraulicandelectricengineeringpress.1999.1~71.〕
〔10〕张海波,张捷斌.新疆流域水资源分配模拟系统设计〔J〕.干旱
区地理,2001,24(3):251~254.〔ZhangHai2Bo,ZhangJie2bin.
ModelingsystemdesignforwaterresourceallocationofriverbasininXinjiang.AridLandGeography,2001,24(3):251~254.〕
〔11〕李新,杨德刚.塔里木河水资源利用的效益与生态损失〔J〕.
干
旱区地理,2001,24(4):327~331.〔LiXin,YangDe2gang,
BenefitandecologicallossofwaterutilizationinTarimriver.Arid
2期 左其亭等:水资源系统中的不确定性及风险分析方法 121
LandGeography,2001.24(4):327~331.〕
Fuzzyriskcomputationmodelanditsapplications.JournalofZhengzhouUniversityofTechnology.2001,22(3):78~80.〕
〔12〕左其亭,吴泽宁.模糊风险计算模型及其应用研究〔J〕.郑州工
业大学学报,2001,22(3):78~80.〔ZuoQi2ting,WuZe2ning.
UncertaintiesinWaterResourcesSystemand
RiskAnalysisZUOQi2ting1,2 3
(1SchoolofEnvironment&,,hou450002,China;
2XinjiangofScience,Urumqi830011,China;
3andManagementBureau,Urumqi830001,China)
Abstract:Onekeyinwaterresourcessystemresearchtodayishowtoquantifyuncertainty,whichisinherentineverywaterresourcessystem.Uncertaintiesaretheoriginofriskinhydrologyandwaterre2sources.Inthispaper,thetypesofuncertaintiesinwaterresourcessystemarecomprehensivelyconcluded,suchasrandom,fuzzy,greyandunascertainty,andtheconceptofuncertainsystemofwaterresourcesisputsfor2ward.Fromthepointsofmathematicalproceedingofuncertaintyfactors,theriskanalysismodelbasedonprob2abilitytheory(relatedtothefieldofrandom)isintroduced,andnewriskanalysismodelsbasedonfuzzysubor2
dinativedegree(relatedtothefieldoffuzzy),greynumber(relatedtothefieldofgrey)andunascertainedmathematics(relatedtothefieldofunascertainty)areputforward.Itlayafoundationforriskplanningandmanagementofwaterresources.
Keyword:WaterResources
System,Uncertainty,Risk,RiskAnalysis