地下水允许开采量确定的风险分析.docVIP

地下水允许开采量确定的风险分析 (1) 　　而风险分析(Risk analysis)则相应包括三个相互联系的部分，即：风险辩识、风险估算和风险评价。风险辩识即为分析系统在特定的决策条件下发生什么样的事故，对地下水供水工程来说，就是分析在某一决策条件下运行供水工程会产生什么环境问题；风险估算则着眼于定量地描述事故发生的概率，即计算出事故发生的可能性有多大；风险评价的主要任务是阐述事故将会产生什么样的社会、经济和环境负效应，即对事故所产生的后果进行评价。 2 风险分析的研究方法 　　目前，用于风险及不确定性分析的方法很多，如：重现期法、直接积分法、蒙特卡罗法、均值一次二阶距法、加强一次二阶距法和二阶距法等。本文所采用的风险分析方法是目前被广泛应用于科研和生产实际工作中的蒙特卡罗方法，该方法的优点在于无需以显式给出风险因子与研究变量之间的关系，而是通过从实测资料的统计分析中，求得其经验关系[1]。 　　蒙特卡罗(Monte-Carlo，简称为M-C)方法，又称统计试验方法，是人工产生和利用随机数方法的总称。它是一类通过对有关的随机变量或随机过程的随机抽样，来求解数学、物理和工程技术问题近似解的数值方法。具体来说，就是对所要求解的问题，构造一种随机变量或随机过程，使其某一数值特征(如数学期望)为所求问题的解，然后对所构造的随机变量或过程进行抽样，并由得到的样本算出相应的参数值，作为所求问题的近似解。 　　用蒙特卡罗方法模拟产生满足一定分布的水文地质参数，首先要用某种特定的方法产生该种分布的随机数，这一过程称为随机抽样。随机变量的分布有多种(如：均匀分布、正态分布、皮尔逊Ⅲ型、F-分布等)，不同分布对应的随机数序列也不同。但就随机数的产生而言，最基本的随机变量是在区间[0，1]上服从均匀分布的随机变量。服从其它分布随机变量的随机数可由[0，1]上均匀分布的随机数变换产生。 　　产生随机数的方法大致可分三类，即随机数表方法、物理方法和数学方法。目前，最为常用的是数学方法。用数学方法产生随机数是通过数学递推式运算实现的，而由此产生的数值序列到一定长度以后，或退化为零，或周而复始地出现周期现象。因此，由数学方法产生的随机数并不是真正的随机数。但是，只要其周期充分长，而且能够通过[0，1]上均匀分布变量所需满足的各种统计检验，如分布的均匀性、抽样的随机性、独立性以及前后的一致性等，则从统计模拟的角度来看，就可以把它们当作随机数使用。为了和真正的随机数相区别，通常把用数学方法产生的随机数称为“伪随机数”。用数学方法产生的“伪随机数”的特点是速度快，占用计算机的内存小，对所模拟的问题可以进行复查，一般又有较好的概率统计性质。 　　用数学方法产生伪随机数的方法有多种，包括迭代取中法、移位法和同余法，而最为常用的是同余法。同余法包括乘同余法、加同余法、混合同余法等[2]。 　　本次采用混合同余法，其递推公式为： xi+1=λxi+c (2) γi+1=xi+1/m (mod m) (3) 　　其中，x,λ,c和m都取非负整数。当c=0时，上式即为乘同余法的递推公式。 　　按上述方法由计算机产生的伪随机数序列为[0，1]上的均匀随机数，在实际问题中尚需变换为给定或已知分布的随机样本值，即进行随机抽样。由于本文研究的水文地质参数大都服从均匀分布或正态分布，因此只需介绍正态随机变量的抽样。 　　正态随机变量的随机抽样方法有：哈斯汀(Hasting)有理逼近方法、统计近似方法和变换方法[3]。本次采用统计近似方法。产生均值为μ、方差为σ2的正态分布随机数Y的计算公式为： (4) 其中n足够大。通常取n=12，其近似程度已是相当好了，此时有： (5) 其中RNi为0到1之间的均匀分布的随机数[2]。 3 地下水允许开采量确定的风险分析 　　对地下水供水水源地的持续运行来说，风险Pf就是指供水工程承担的供水外来负荷L(Loading)大于工程本身的承载能力R(Resistance或Capacity)的概率： Pf=P(LR) (6) 式中：L与R均可为包括众多水文地质参数的函数。对地下水供水水源地来说，L主要为地下水的开采量，R则为地下水的允许开采量。因此，允许开采量的确定是非常重要，其数值的准确与否直接关系着水源地能否持续运行。允许开采量的计算方法有很多，但每一种都包含着一些不确定性因素。例如在最简单的通过计算地下水的补给量来确定允许开采量的方法中，要涉及到降水入渗补给量、河道渗漏补给量、灌溉水回渗量、其它含水层的越流补给量、地下水的侧向径流补给量等。 　　本文作者长期在山东省济宁市开展地下水资源评价与此同时管理工作，掌握了大量的地质及地下水的动态资料等，且深感过去的资源评价仅局限于确定性模型研究的不足。为此，选定济宁市作为研究实例。对于济宁