1、第 卷 第 期 年 月人 民 长 江 ,收稿日期:基金项目:国家自然科学基金项目();贵州省科技计划资助项目(黔科合平台人才 号)作者简介:胡 贝,女,硕士研究生,研究方向为矿山安全与灾害防治。:通信作者:李希建,男,教授,博士生导师,博士,研究方向为煤矿瓦斯灾害防治、煤层气开发。:文章编号:()引用本文:胡贝,李希建,代芳瑞,等 基于组合赋权和云模型的危岩体危险性评价 人民长江,():基于组合赋权和云模型的危岩体危险性评价胡贝,李 希 建,代 芳 瑞,汪 圣 伟,(贵州大学 矿业学院,贵州 贵阳;复杂地质矿山开采安全技术工程中心,贵州 贵阳;贵州大学 瓦斯灾害防治与煤层气开发研究所,贵州 贵
2、阳)摘要:为更加合理地判别危岩体的危险性,以具有丹霞地貌特征的危岩体为研究对象,从危岩体不稳定性、危破坏性和量观文化价值 个方面着手,选取 个评价指标构建危岩体危险性评价指标体系;结合层次分析法()、基于指标相关性的指标权重确定方法()和云模型(),构建组合赋权 云模型危岩体危险评价模型,并选取 组测点进行模型可靠性验证。结果表明:各评价指标权重排序为主控结构面特征 危岩规模 边坡坡度 科学文化价值 景观价值 失稳后致灾可能性 坡形 水文条件;基于组合赋权 云模型得到的各组危岩体危险性判别结果与实际测量结果基本吻合,评价结果准确性高,模型具有可靠性。关 键 词:危岩体;危险性;组合赋权;云模型
3、;确定度中图法分类号:文献标志码:引 言不同工程类型的岩体受构造运动、自然位置和人类活动的影响,逐渐发育成危岩体,随之造成山体滑坡、崩塌等灾害,从而给人员安全及工程建设带来极大威胁。因此,科学准确地判定危岩体危险性,并采取针对性的预防措施变得尤为迫切。学者们在分析危岩体危险性方面已取得了一定进展。唐红梅等基于 的方法构建模糊综合评价模型,并通过为各指标赋加权重获得了崩塌体危险性评估结果。刘晓然等运用可变模糊集理论建立了崩塌危险性评价指标的相对隶属函数,采用离差最大化方法计算各指标的权重。陈宙翔等借助无人机勘查公路边坡危岩崩塌路段信息,进而利用 软件模拟危岩体崩落运动轨迹以评价公路边坡危岩崩塌的
4、危险性。徐伟等通过对郭达山后山危岩进行岩体分类与质量评定实现危岩危险性的划分。张硕等根据坠落式危岩体的稳定性极限状态特征实现崩塌危险性评估。高永才等通过岩石力学实验,采用静力解析法判定景区危岩体危险性。高相波等基于逆向工程建模方法利用 有限差分数值计算软件进行危岩体危险性分析。黎尤等采用三维模拟软件 对汶川县三官庙村崩塌部位危岩崩塌灾害的运动特征进行模拟,分析评价了潜在危岩区的危险性,并得出其危险影响范围。唐红梅等在层次分析法的基础上,利用 技术成功地对巫山县进行了危险性分区。上述研究采用了不同方法对危岩体的危险性进行分析,但危岩体的危险性影响因素多且难以量化,故本文引入云模型,通过定性概念与
5、定量计算之间的任意转换,实现危岩危险性等级的判别。危岩体受流水、风化和重力等外力因素的共同作用,发育形成丹霞地貌。丹霞地貌以红色陡崖坡为主,坡面多呈顶斜、身陡、麓缓状态,受自身的不稳定以 人 民 长 江 年及自然环境、人类活动的影响,易出现崖体破裂、崩塌、滑坡和泥石流等灾害现象,严重制约着人类活动及工程建设。鉴于此,为判别危岩体危险性等级,笔者以具有丹霞特征的危岩体为研究对象,分析了影响危岩危险性的主要因素,构建了危岩体危险性评价指标体系,再采用层次分析法(,)、基于指标相关性的指标权重确定方法(,)和云模型(,)建立危岩体危险性等级评价模型,并结合实例验证其准确性。危岩体危险性因素分析危岩体
6、危险性主要与危岩体不稳定性、危岩体破坏性和危岩体价值 个因素有关。危岩体不稳定性()坡型。坡面平滑且不存在临空状态,则难形成危岩。若坡面粗糙且临空环境好,则易造成滑塌或坠落。()坡度。危岩的坡度越大,其稳定性越差;坡度越小,其稳定性就相对越大。因此,危岩的稳定性与坡度呈负相关。()主控结构面特征。主控结构面特征主要包含结构面粗糙度、填充情况、宽度和角度。其中结构面的角度与危岩稳定性有着更直接的关系。为了避免各影响因素之间的相互影响作用,选取主控结构面的倾角作为评价危岩危险性的主要指标之一。()岩性。单一均质岩体,结构稳定,难造成破坏。而由多种相似岩层组成的岩体,稳定性差,易受扰动发生破坏。()
7、水文条件。自然工况下,大部分危岩及边坡处于基本稳定状态。但在暴雨工况下,危岩及边坡易出现失稳破坏的现象。危岩体破坏性()危岩规模。危岩体规模主要以危岩体方量来判别,危岩体的规模与危岩体危险性呈正相关关系,即规模越大,造成的危害性就越大。()失稳后致灾可能性。危岩经过运动与承灾体发生碰撞,对承灾体造成损害,失稳后能否与承灾体发生接触是能否成灾的关键,以此用危岩体与承灾体之间的水平来衡量致灾可能性。危岩体价值()科学文化价值。危岩为普通岩体,价值不大;危岩科学文化价值较小;危岩科学文化价值较大;危岩稀有,科学文化价值很大。()景观价值。危岩为普通岩体,无较大价值;危岩具有景观价值较小;危岩具有景观
8、价值较大;危岩景观价值很大,且稀有。等级划分基于上述对危岩体危险性的影响因素分析,以丹霞地貌为例构建危岩体危险性评价指标体系及等级划分表(见表)。由于各因素都是定性阐述,难以定量计算,参照文献采用的量化方法可将各个影响因素分为四级,分别赋值、,赋值数值越高,说明丹霞地貌区域的危岩体危险性越大。将危岩体的危险性等级分为级危险性大,即危岩体灾害发育程度为强发育;级危险性中等,即危岩体灾害发育程度为中等发育;级危险性小,危岩体灾害发育程度为弱发育。再根据量化分值进行划分:级(分),级(,分),级(分)。表 危岩体危险性影响因素等级划分 评价指标危险性等级标准目标层准则层因素层级级级危岩体危险性不稳定
9、性坡形凸曲线、缓陡缓、洞穴型直坡面、陡缓陡、负地形曲线阶梯曲线评价边坡坡度 岩性组合极软岩含软弱夹层组合、均质软岩 较软岩均质坚硬岩或较坚硬岩主控结构面角度 水文条件经年有水流经年有水湿润或雨季有少量水经年无水危破坏性危岩规模 失稳后致灾可能性水平距离 水平距离 水平距离 自身价值科学文化价值极大科学文化价值,稀有具有一定的科学文化价值普通岩体,无科学文化价值景观价值极大景观价值;代表性价值具有一定的景观价值无景观价值 组合赋权确定指标权重 层次分析法根据危岩体危险性 个影响因素指标构造初始评判矩阵,求出各指标所占权重。主观权重向量式为?()式中:代表各指标主观权重;?代表判断矩阵行乘积的 次
10、方根。第 期 胡 贝,等:基于组合赋权和云模型的危岩体危险性评价 法赋权算法 法是一种由指标间的对比强度和指标间的冲突性来衡量指标客观权重的量化方法。依据危岩体危险性评价指标,利用()对其加以改进,详细步骤如下:()构建初始评判矩阵。设有 个评价对象,个评价指标,初始评判矩阵见式()。()式中:表示第 个评价对象第 个指标所对应的评分数据 ,;,()。()标准化矩阵。根据 方法标准化处理 得:()式中:?;?,(?)。()计算变异系数。()式中:表示第 个指标的变异系数。()计算独立性系数。根据标准化矩阵 经过皮尔逊系数计算法,得到相关系数矩阵 (),并算得独立性系数:()()()综合性系数。
11、()()确定客观权重。()组合权重为更加全面地考虑指标重要度,引入拉格朗日乘数法对主观赋权 和客观赋权 进行融合,获得更加真实的权重值,进而计算得到各指标的组合权重 (,),即 ()式中:代表评价指标的组合权重,代表主观权重,代表客观权重。构建组合赋权 云模型危岩体危险评价模型 云模型的生成利用 软件,通过正态云发生器来实现危岩体危险性的定性定量转换,算法步骤如下:()生成以 为期望值及 为方差的随机数;()生成以 为期望值及 为方差的随机数;()计算危岩体危险性评价指标确定度值:()()式中:表示确定度。云模型中,云的特征值借助期望、熵 和超熵来反映。期望 表示云的重心位置,熵 表示定性概念
12、的模糊程度,超熵 表示熵值的熵。则危岩体危险性的评价指标隶属于不同等级的数字特征值可依据式()求得。()式中:,分别代表同一等级标准的上限值、下限值;取值为 。根据公式()算得危岩体危险性评价指标云模型的数字特征值,如表 所列。表 云模型的数字特征值 危险等级(,)(,)(,)(,)在已知云模型的数字特征值的基础上,通过 软件生成的 朵正态云,且分别对应一个危险性等级,如图 所示。横轴表示危岩体评价指标值,纵轴表示确定度,因此可以很直观地表示评价指标评分分值、状态等级与确定度三者之间的隶属程度关系。根据正态云发生器的算法,分别计算危岩体危险性的各评价指标数据 隶属于云 的确定度 (,;,)。再
13、结合各评价指标的组合权重,算出评价对象的综合确定度 表达式为()人 民 长 江 年式中:代表综合确定度,代表确定度,代表各指标组合权重。危岩体危险性等级 判别模型为 ,()图 危岩体隶属于各危险等级的云模型 评判流程组合赋权 云模型危岩体危险评价模型的评价过程为:按照危岩体危险性等级划分标准求出相应等级的云数值特征值;由 软件生成危岩体隶属于各危险等级的云模型,并根据样本数据按照正态云发生器的算法计算各指标对应于各级别的确定度;结合各指标的组合权重,计算综合确定度;依据最大隶属度原则判别危岩体的危险性等级。具体流程如图 所示。图 危岩体危险性判别流程 工程应用 工程概况选取广东省韶关市丹霞山景
14、区内的丹霞地貌作为实际工程案例,验证所建模型的可行性。丹霞山地处广东省韶关市仁化县,是以丹霞地貌景观为主的风景区,该景区岩体多由红色砂烁岩构成,且因景区内部具有赤壁丹崖等特色而被称之为丹霞地貌。鉴于丹霞地貌坡面陡峭,岩性力学强度较差,垂直节理较发育,易发育较多危岩,选择该区域进行研究有一定的现实意义。根据现场勘查结果对不同测试地点的危险性评价指标进行量化评分,详情见表。表 景观危岩体危险性评分 测点编号 组合权重的确定危岩体危险性评价指标的组合权重见表。表 危岩体危险性评价指标权重 指标主观权重客观权重组合权重坡形 边坡坡度 岩性组合 主控结构面角度 水文条件 危岩规模 失稳后致灾可能性 科学
15、文化价值 景观价值 由表 可以看出,各评价指标权重排序为 主控结构面特征 危岩规模 边坡坡度 科学文化价值 景观价值 失稳后致灾可能性 坡形 水文条件。由各指标权重结果可以看出,主控结构面特征是影响危岩体危险性的重要因素,且主控结构面的倾角越大的景区地点,危岩体越易发育,表明此处的危岩体越不稳定,危险性越大。危险性等级评价结果以 测点为例说明计算确定度的过程。根据表 的云模型数字特征值,并将案例中 测点对应的 个指标数据代入到公式(),算得每个指标隶属于各危险等级的确定度,见表。通过将表 中的各指标组合权重和表 中 测点的各指标 第 期 胡 贝,等:基于组合赋权和云模型的危岩体危险性评价确定度
16、代入到公式()中,算得综合确定度 (,),再由公式()可知()()(),因此 测点的危险性隶属于级的可能性最大。同理依次得到 组测点的综合确定度,得出各片区危岩体危险性等级判定结果,见表。表 测点的各指标确定度 确定度()()()表 危岩体危险性判别结果 测点编号综合确定度()()()判别结果实测结果 由表 可知,选取的 组测点中有 组的危险性等级为级,危险性大,有 组的危险性等级为级,危险性中等。根据所选地点危险性评价结果可知,该景区危岩体发育较多,发生危岩体破坏的可能性也随之增大,因此针对危险性大的危岩应采取主动防护的手段加以治理。通过对此表中危险性等级,可以看出云模型得到的各危岩体危险性判别结果与文献中的实际测量结果基本吻合,验证了该模型的准确性。值得注意的是,测点 和 的危险性等级与实测结果不一致,这主要是由于丹霞地貌所处地势复杂,但该模型评价指标级别还不够精细所致。未来,应通过调查不同地区的危岩体,对危岩体危险性评价指标级别进行更加精细地划分。结 论本研究从影响危岩体危险性的危岩体不稳定性、危岩体破坏性和危岩体价值 个主要方面,选取 个评价指标构建危岩体危险性评价指标体系。采