流泥石流危险性评价已经成为当前国内外灾害科学研究的热点之一,也是泥石流预测预报和防灾减灾的重要内容。国外对泥石流危险性评价始于20世纪70年代末期,日本在1977年首次对泥石流危险性开展了研究提出了一系列泥石流危险性判定的方法,到90年代初,各种数学模型如线性回归、层次分析法、灰色理论、遗传神经网络、模糊数学等方法开始在这一研究中得以广泛应用和发展
泥石流危险性是根据泥石流危险因子评价的,刘希林采用灰色系统理论中的灰色统计方法,根据专家通信调查总结出2项主要危险因子以及8项次要危险因子,并得到广泛应用。目前这些广泛应用的评价方法简单易行,但多是将刘希林提出的10个评价因子直接或者经过简单转换之后代入相应的公式得到泥石流危险度,未能充分考虑各评价因子与危险度的相关性,并且在某些情况下,获取以上全部评价指标是较为困难的,因此需要对以上指标进行筛选,以实现利用较少的指标进行危险性评价的目的,通径分析方法则能较好地解决上述问题,考虑了各影响因子与危险度之间的相关性,根据相关程度的强弱对其进行筛选,使得根据较少的评价因子实现危险性评价成为可能。
本文以大西沟泥石流为例,利用通径分析方法,对以上因子进行进一步的筛选,选择了与危险度相关程度较强的7个评价因子,并利用筛选后的评价因子以可拓学理论为基础确定泥石流危险性等级,为在仅有少量评价指标或不易获得泥石流全部基本特征指标的情况下进行泥石流危险性评价提供了一种新的思路,同时为利用通径分析与可拓理论进行泥石流危险性评价提供了参考。
通径分析是美国数量遗传学家Sewall Wright于1921年提出来的一种多元统计技术,它是简单相关分析的继续,在多元回归的基础上将相关系数加以分解,通过直接通径、间接通径及总通径系数分别表示某一变量对因变量的直接作用效果、通过其他变量对因变量的间接作用效果和综合作用效果[2,8,9]。
对于一般的多元线性回归分析,设有自变量X1,X2,,Xk和因变量Y,
利用最小二乘法求出式(4)各自变量线性回归系数的求解模型,在此基础上,进行一定的数量变换,则可得出如下各简单相关系数的分解方程:
以上为通径分析的基本模型,每一个自变量与因变量之间的简单相关系数可被分解为直接通径效果部分和间接通径效果部分,其中rij为各自变量的简单相关系数,riy为自变量与因变量之间的简单相关系数,PiY为直接通径,表示Xi对Y的直接影响效应,rijPjY为间接通径,表示Xi通过Xj对因变量Y的间接影响效应。
笔者以大西沟地区30条泥石流沟的实际调查资料为基础,对刘希林模型所提出的流域面积X1、主沟长度X2、流域最大相对高差X3、流域切割密度X4、主沟床弯曲系数X5、泥沙补给段长度比X6、24 h最大降雨量X7和流域内人口密度X8等8项次要危险因子以及一次泥石流最大冲出量X9、泥石流发生频率X102个主要危险因子进行通径分析。
将利用刘希林模型计算出的泥石流危险度作为因变量Y,某地区30条泥石流沟的基本参数与危险度见基于通径分析与可拓学的公路泥石流危险性评价。利用SPSS16.0软件进行分析计算,得到各变量与因变量之间的相关系数,见中的计算结果可知,10个评价因子中对泥石流危险度影响最主要的几项从大到小分别是X2(主沟长度)、X9(一次泥石流最大冲出量)、X1(流域面积)、X6(泥沙补给段长度比)、X3(流域最大相对高差)、X5(主沟床弯曲系数),其简单相关系数分别达到0.720、0.707、0.641、0.602、0.573、0.230,相对而言,其他因子对其影响较小,如人口密度与危险度的简单相关系数仅为0.019,远不及上述几项因子。此外,24 h最大降雨量是诱发泥石流暴发的一个不可忽略的重要因素,但由于上文中作为基础数据的30条泥石流沟处于相近的区域范围内,其降雨量取值相同,因此导致在上述计算分析过程中并未体现出24 h最大降雨量对泥石流危险度的影响,所以在进行危险性评价时,可以主要依据于上述分析筛选出的6项因子,同时考虑24 h最大降雨量这一因子,使得在保证评价结果正确性的前提下,以减少评价因子的数量,简化评价过程。
在以上确定的评价因子的基础上,利用可拓学方法对其进行危险性评价[10,11,12]。
由集合P、它的n个特征Ci以及将标准量值范围加以拓展了的取值范围Vpi=pi,bpi组成的物元称为节域物元,可表示为:
对于待判别的泥石流沟,将调查结果用物元表示,称为泥石流危险性评价的待评物元,表示为:
根据可拓集合的关联函数[13],可得到待评泥石流沟物元关于各类等级的关联度,第i条泥石流沟的第j个因子危险性类别的关联度为:
式中:Wi为各评价指标的权系数,可通过专家打分法、简单关联函数法、按评价因子的灰关联排序等方法来确定[2,11,1]。若Kt0(p)=maxKt(p),则评定p属于等级t0。
大西沟泥石流流域面积大约为0.67 km2,高程1030~1300 m,泥石流主沟曲线 km,直线 km,泥沙补给段长度比大约为51%,一次最大冲出量大约为1.93104m3,当地24 h最大降雨量为217 mm。该泥石流沟的三维立体图见。
本文通径分析结果显示,在目前泥石流危险性评价时常用的10个因子中,主沟长度、一次泥石流最大冲出量、流域面积、泥沙补给段长度比、流域内最大相对高差、主沟床弯曲系数其简单相关系数分别达到0.720、0.707、0.641、0.602、0.573、0.230,同时由于24 h最大降雨量对泥石流的启动具有至关重要的作用,因此本文最终选取以上7个因子并利用可拓学方法进行泥石流危险性评价,评价方法和结果都具有较好的可信度。
与其他评价模型相比,本文中评价因子的筛选过程考虑了各项因子与泥石流危险度相关程度的强弱,使得评价结果更加合理可信,同时减少了评价因子的数量,特别是减少了泥石流发生频率这类不易获取的因子,使得泥石流危险性评价更加简化易行。在对各项因子进行筛选后将可拓学方法应用到泥石流危险性评价中,评价过程合理,能较真实地反映实际泥石流的危险等级,可为类似工程的危险性评价提供可行的思路。
泥石流灾害及其孕灾环境具有高阶模糊性，现有的模型和方法在表示和分析泥石流灾害时存在明显的尺度效应。首先提出了我国泥石流危险性评价的尺度模型，通过尺度分析的方法分析了孕灾因子的适宜尺度区间。然后基于区间II-型模糊集表示孕灾因子的高阶模糊隶属度，以此描述泥石流及其孕灾环境的复杂性；以区间权重描述各孕灾因子的重要性，并提出了适宜尺度区间上的危险性评价方法。试验表明，这种方法能够克服许多单尺度分析方法的缺陷，在泥石流危险性评价过程中能充分利用孕灾因子在适宜尺度区间上的可能状态，能顾及专家的对孕灾因子重要性评价的不确定性，评价结果比较客观。耀世注册=登录首页？
