Scientia Geographica Sinica  2012 , 32 (7): 885-891 https://doi.org/10.13249/j.cnki.sgs.2012.07.885

Orginal Article

岷江上游潜在性泥石流堰塞湖危害及判识

柳金峰12, 游勇12, 陈兴长3

1.中国科学院山地灾害与地表过程重点实验室, 四川 成都 610041
2.中国科学院水利部成都山地灾害与环境研究所, 四川 成都 610041
3.西南科技大学环境与资源学院, 四川 绵阳 621010

Potential Damages and Identification of Debris Flow Barrier Lakes in theUpper Min River Drainage

LIU Jin-feng12, YOU Yong12, CHEN Xing-chang3

1. Key Laboratory of Mountain Surface Process and Hazards, Chinese Academy of Sciences, Chengdu, Sichuan 610041, China
2. Institute of Mountain Hazards and Environment, Chinese Academy of Sciences, Chengdu, Sichuan 610041, China
3. School of Environment and Resources, Southwest University of Science and Technology, Mianyang, Sichuan 621010, China

中图分类号:  P642.23

文献标识码:  A

文章编号:  1000-0690(2012)07-0885-07

收稿日期: 2011-05-3

修回日期:  2011-09-12

网络出版日期:  2012-07-20

版权声明:  2012 《地理科学》编辑部 本文是开放获取期刊文献,在以下情况下可以自由使用:学术研究、学术交流、科研教学等,但不允许用于商业目的.

基金资助:  国家科技支撑计划课题(2009BAK56B05)、国家自然科学基金项目(40971014)资助

作者简介:

作者简介:柳金峰(1979-),男,安徽泾县人,助理研究员,主要从事山地灾害实验与防治工程研究。E-mail:liujf@imde.ac.cn

展开

摘要

汶川地震造成泥石流形成条件的改变,其次生灾害堰塞湖的危害已开始显现。如何对震后潜在性泥石流堰塞湖进行判识,成为迫切需要解决的问题。选取岷江上游映秀至汶川段为研究区,通过分析震后泥石流形成条件的变化、典型泥石流堰塞湖的危害及松散物质储量,选取潜在性泥石流堰塞湖的判识指标,利用模糊物元可拓模型,建立潜在性泥石流堰塞湖的综合判识模式。通过判识,研究区形成泥石流堰塞湖可能性高的一级支沟有17条,主要集中分布在映秀镇至草坡乡段,此段将是今后受堰塞湖危害的高危地段。

关键词: 岷江上游 ; 泥石流 ; 堰塞湖 ; 潜在危害 ; 判识

Abstract

The Wenchuan earthquake caused numerous landslides and collapses that provided abundant loose solid materials for future mobilization as debris flows. The formation conditions for debris flows were changed consequently. Debris flows will be very active for a long time in the affected area. Their secondary disaster, barrier lake, is beginning to appear. Therefore, it is urgent to identify those debris flow gullies which are potential to form barrier lakes following the environmental changes caused by the earthquake. This paper selected the upper Min River from Yinxiu to Wenchuan as study area, analyzed the changed formation conditions of debris flows and the damages of a typical debris flow barrier lake, and interpreted the loose solid materials from TM satellite images and aerial photographs. Then, the identification factors were selected and classified based on the factor independence analysis. At last, the fuzzy matter-element extension theory was used to construct a model to identify potential barrier lakes. The potential future increase in debris flow activity relates directly to the great increase in loose materials. The interpretation results indicated that the area of the landslides and collapses caused by the earthquake is 39.40 km2, accounting for 6.97% of the total study area. Assuming 55% of the landslide and collapse area were their deposition area and the average deposition depth was 10 m, the volume of the loose solid materials can be estimated as 2.13×108m3 which was sufficient for debris flow formation. Based on the factor independence analysis, the volume of the loose solid materials per unit area, the discharge ratio between the debris flow gully and the main river, the slope gradient of debris flow channel, the main river width, and the included angle between the debris flow gully and the main river were determined as identification factors. According to the multi-factor comprehensive identification, among the 55 first-grade gullies, 17 gullies have high probabilities for potential debris flow barrier lakes, 15 have moderate probabilities, and 23 gullies have low probabilities. The most likely debris flow gullies for forming barrier lakes mainly distribute in the middle-to-lower reaches of the study area (Yingxiu Town to Caopo Town), and those with moderate and low probabilities are mainly distributed in the upper reaches of the study area (Caopo Town to Wenchuan County). In addition, among the 17 potential debris flow gullies with high probabilities, 4 gullies already yielded barrier lakes during the 2008-2010 rainy seasons. The remaining debris flow gullies will produce barrier lakes in the future with high probability.

Keywords: Upper Min River ; debris flow ; barrier lakes ; potential damages ; identification

0

PDF (698KB) 元数据 多维度评价 相关文章 收藏文章

本文引用格式 导出 EndNote Ris Bibtex

柳金峰, 游勇, 陈兴长. 岷江上游潜在性泥石流堰塞湖危害及判识[J]. , 2012, 32(7): 885-891 https://doi.org/10.13249/j.cnki.sgs.2012.07.885

LIU Jin-feng, YOU Yong, CHEN Xing-chang. Potential Damages and Identification of Debris Flow Barrier Lakes in theUpper Min River Drainage[J]. Scientia Geographica Sinica, 2012, 32(7): 885-891 https://doi.org/10.13249/j.cnki.sgs.2012.07.885

引 言

汶川地震诱发的大量崩塌、滑坡提供了丰富的松散物源,泥石流灾害会在未来相当长一段时间内处于活跃期[1,2]。其次生灾害堰塞湖的危害已开始显现,如岷江上游震后2008~2009年雨季只形成了磨子沟泥石流堰塞湖,但是到2010年就形成4个堰塞湖,且造成了较大危害,如2010年8月14日,映秀镇红椿沟暴发特大泥石流,泥石流直接冲入岷江,形成了长450 m、宽300 m、厚度达30 m的堰塞体,堵塞导致岷江改道,洪水涌入重建的映秀新城。震后泥石流的活跃导致堵河形成堰塞湖的可能性大大增加,因此,如何对震后潜在性泥石流堰塞湖进行判识,以保障灾区人民生命财产的安全,成为迫切需要解决的问题。本文选取岷江上游映秀至汶川段为研究区,通过分析震后泥石流形成条件的变化、典型泥石流堰塞湖的危害,选取潜在泥石流堰塞湖的判识指标,利用模糊物元可拓模型,建立潜在泥石流堰塞湖的综合判识模式。

1 背景概况

1.1 研究区概况

本文选取岷江上游映秀-汶川段为研究区域。研究区域位于四川省阿坝州汶川县境内,集水区面积565 km2,沿岷江的一级支沟共有55条(图1)。

图1   研究区地质、降水及地震前后泥石流沟分布

Fig.1   The geology, rainfall of the study area and debris flow distribution before and after earthquake

研究区位于青藏高原的东部边缘地带,属第一阶梯青藏高原向第二阶梯四川盆地及边缘山地过渡地带。地貌类型复杂,以高山峡谷为主,相对高度一般为1 000~4 000 m[3]。研究区域土地主要集中分布在海拔1 000~4 000 m的高程范围内,占总面积的93.0 %,其中海拔≥2 000 m高程的面积占总面积的77.7%。对研究区域山坡坡度进行分级统计可知,流域内≥25°的陡坡地占87.2%,可见研究区的山坡坡度主要以陡坡为主。

研究区位于四川盆地边沿的秦岭纬向构造带、龙门山北东向构造带与马尔康北西向构造带间的三角形地块内,在整体上属松潘-甘孜地槽褶皱系的北段,龙门山华夏系、岷江经向系、岷江旋扭构造系三大构造体系在区内均有展布。区内褶皱和断裂发育,尤其是深大断裂发育,如茂汶断裂、岷江断裂、米亚罗断裂等。出露的地层从元古界到新生界都有,岩性主要以陆相碎岩、火山碎岩、碳酸盐岩和区域变质岩为主(图1)。

研究区处于海陆季风气候向高原季风气候的过渡带,受水气来源和地形的影响,岷江上游各地的降水有较大差异(图1),干湿季节明显,全区多年平均降雨量在750 mm左右[4]。汶川县的漩口是岷江上游的降水中心,年降水量达1 327.6 mm。区内降水在季节上分配不均匀,降水主要集中在下半年(5~10月),占全年降水量的80%~90%,且降水强度大,一日最大降水量为79.9 mm。

综上,研究区的地形、地质和降水条件有利于泥石流的形成,该区域历来就是泥石流灾害多发的区域[4]。据震前统计,研究区共分布有泥石流沟15条(图1)。地震造成研究区物源的剧增,泥石流活动已经进入活跃期,并会持续很长一段时间,高度活跃的泥石流进入岷江后,不仅会抬升、挤压主河床,同时堵塞主河形成堰塞湖的可能性将会大大增加。

1.2 研究区堰塞湖危害概况

据野外考察,研究区震后在2008~2010年的雨季中,有22条沟暴发了泥石流,其中只有两条为震前老泥石流沟,其余20条是震后新增的泥石流沟(图1)。在已暴发的泥石流沟中,有4条(红椿沟、小磨子沟、关山沟和磨子沟)形成了堰塞湖,并造成了一定的危害。下面以岷江右岸的一级支沟关山沟为例,介绍泥石流堰塞湖及其危害情况。

关山沟位于汶川县银杏乡境内,该沟距离映秀镇城区约5.4 km,距“5.12”汶川地震的宏观震中牛圈沟的直线距离为6.7 km(图1)。关山沟在汶川地震前并未暴发过泥石流,沟内植被条件较好,松散物源少。震后,松散物源剧增,地震引发的崩塌、滑坡的总方量约122×104m3,如此巨大储量的松散物源,只要遇到一定的降水条件,非常容易形成泥石流[5]。据野外考察,震后2008~2010年雨季,该沟共计暴发了十余场泥石流。2008年雨季后,泥石流在沟口形成了一个长约80 m,宽约150 m的堆积扇,泥石流并未造成岷江堵塞,而是将岷江逼向对岸,将岷江挤压成只有约5 m宽的河道(图2a,拍摄于2008年10月);2009年雨季后,沟口的泥石流堆积扇不断扩大,扇长约230 m,宽约250 m,泥石流堵断岷江,但是很快就溃决了,并未造成较大危害(图2b,拍摄于2009年12月);2010年雨季后,关山沟泥石流继续暴发,堵塞岷江,形成堰塞湖(图2c,拍摄于2010年11月),回水淹没距沟口上游1.6 km岷江右岸的银杏乡东界脑村(图1),97户新建的灾民重建房全部被淹,平均淹没深度约1.5 m(图2d,拍摄于2010年11月)。

图2   关山沟泥石流堰塞湖及其危害

Fig.2   The debris flow barrier lake of Guanshan gully and its damages

2 震后泥石流形成条件的变化

震后泥石流形成条件的最大变化主要由松散物质条件的变化所决定。松散物质丰富的沟道,在较低降雨条件下很容易起动形成泥石流;而那些松散物质条件相对较少的沟道,无论是形成泥石流的规模还是暴发频率,都不及松散物质丰富的沟道。松散物质条件是震后泥石流形成的关键因素[6]

本文采用TM卫星影像面上解译、典型支沟航片解译、典型灾害点现场调查三者相结合的方法对研究区震后松散物质进行了解译。数据主要采用震后2008年雨季前CBERS-2和TM卫星影像、1∶25 000和1∶5 000航空影像(图3)。

图3   研究区崩塌和滑坡分布

Fig.3   The distribution of landslides induced by earthquake

结果表明,研究区地震直接引发的崩塌、滑坡的总地表面积为39.40 km2,占总面积的6.97%。假定崩塌、滑坡总地表面积的55%为其堆积区面积,且平均堆积厚度为10 m,可估算得到研究区松散物质的总方量为2.13×108m3。由图3可知,地震直接引发的松散物质从下游映秀镇向汶川县城递减,这与距离发震断层北川映秀断裂带的远近有关,据发震断层越近,其产生的崩塌、滑坡越多。充足的松散物源使泥石流暴发的频率和规模大大增加,加之岷江两侧山坡产生的大量松散物质直接进入主河,造成主河床抬升和变窄,泥石流直接进入主河后,更容易造成堵塞,形成堰塞湖,从而造成更大危害。

3 潜在性泥石流堰塞湖多因子综合判识

3.1 判识指标

本文采用多因子综合判识的方法,对可能造成堵河形成堰塞湖的泥石流沟进行判识。从泥石流本身的活动强度和沟口处主河的条件出发,选取单位面积松散物源方量、支主流量比、泥石流沟床比降、主河宽度、主支夹角5个判识指标。结合目前现有文献研究[7~9],指标的分级见表1表1中支主流量比指标是支沟泥石流峰值流量(P=5%)与主河流量的比值,泥石流流量采用配方法计算[10],其中堵塞系数的取法采用震后修改值[11]。主河流量采用雨季(5~10)月平均最小流量,数据来源于汶川县威州水文站观测数据(2003~2008年)。

表1   潜在性泥石流堰塞湖判识指标分级

Table 1   The classification of the indices for potential debris flow barrier lakes identification

泥石流堰塞湖发生可能性
单位面积松散物源
方量c1(×104m3/km2
<1010~50>50
支主流量比c2<0.50.5~1>1
泥石流沟床比降c3(‰)<100100~300>300
主河宽度c4(m)>10050~100<50
主支夹角c5(°)<3030~90>90

新窗口打开

指标选取后,需对各个指标进行因子的独立性分析[12]。本文采用SPSS 18.0软件的非参数检验工具,将选取的因子分为5组,进行多个独立样本的检验。通过软件提供的三种检验方法:Kruskal Wallis方法、中值方法和Jonckheere-Terpstra 方法检验后,其渐近显著性值均小于0.05,表明在95%置信水平下,选取的5个指标之间有显著性差异,即五个指标是相互独立的。

由于各个指标的作用不同,因此,在潜在性泥石流堰塞湖判识之前,需要确定各个指标的权重。本文采用层次分析法确定各个指标的权重,首先构造判断矩阵如下:

Z=Cij=c1c2c3c4c5c112355c2121235c31312123c415131212c5151513121

通过一致性检验,确定单位面积松散物源方量、支主流量比、泥石流沟床比降、主河宽度、主支夹角5个因子的权重依次为0.45、0.26、0.15、0.09、0.05。

3.2 判识方法

本文采用模糊物元可拓性理论对潜在性泥石流堰塞湖进行多因子综合判识。模糊物元模型是通过引进物元,并对其进行变换和运算,从而从定性和定量两个角度去研究解决问题的规律和方法;可拓集合和关联函数是可拓学定量化工具,可拓方法的理论基础是物元理论和可拓集合理论[13~17]。研究区潜在性泥石流堰塞湖判识的计算步骤如下:

1) 物元分析。物元分析是用来处理在某些条件下,用通常的方法无法达到预期目标不相容问题规律的一种分析方法。在物元分析中,物元的表达形式为:

R=(M,C,X) (1)

式中,R为物元,为潜在性泥石流堰塞湖判识;M为所描述的事物,为所需评价的潜在性泥石流堰塞湖可能性的高低;C为所描述的事物的特征,为选取的指标体系;X为量值,为指标的量化值。

2) 物元矩阵。根据选取的5个潜在泥石流堰塞湖判识的指标,构建物元矩阵如下:

R=Mc1x1c2x2c3x3c4x4c5x5(2)

3) 节域对象物元矩阵。根据表1中选取的5个指标的取值范围,构建节域对象物元矩阵如下:

R=MPc1[ap1,bp1]c2[ap2,bp2]c3[ap3,bp3]c4[ap4,bp4]c5[ap5,bp5](3)

式中,MP为由标准事物加上可转化为标准的事物组成的节域对象;xpi为节域对象关于特征ci的量值范围,xpi=[api,bpi]。

经典域对象物元矩阵可表示为:

R=MBc1[aB1,bB1]c2[aB2,bB2]c3[aB3,bB3]c4[aB4,bB4]c5[aB5,bB5](4)

式中,MB 为标准对象;xBi=[aBi,bBi],表示标准对象MB关于特征ci的量值范围。显然有xBixpi(i=1,2,…,5)。根据表1中五个选取指标的分级范围,分别构建潜在性泥石流堰塞湖可能性高、中、低三种情况下的经典域对象物元矩阵。

4) 关联函数。在物元评价中,关联函数使解决不相容问题的结果量化。若区间x0=[a,b],x1=[c,d],且x0x1,则关联度为:

kj(xi)=-ρ(xi,xij)xijxixijρ(xi,xij)ρ(xi,xpi)-ρ(xi,xij)xixij(5)

式中,

ρ(xi,xij)=xi-(aij+bij)/2-(bij-aij)/2(6)

ρ(xi,xpi)=xi-(api+bpi)/2-(bpi-api)/2(7)

5) 计算待评价单元综合关联度,并确定其所属等级。根据上式得到的关联度值,待评价单元关于潜在性泥石流堰塞湖可能性等级j 的综合关联度按下式计算:

αj=i=15wijkj(xi)(8)

式中,αj为综合关联度,wij为评价指标的权重。

最后通过比较待评价单元各个等级的综合关联度值,其值最大的等级即为其所属的潜在性泥石流堰塞湖可能性等级。

3.3 判识结果

本文以岷江上游映秀-汶川段55条一级支沟为评价单元,根据上述判识指标和标准物源模型,利用Matlab7.0软件编程进行计算,对研究区的潜在泥石流堰塞湖进行判识(图4表2)。通过判识,研究区55条一级支沟中,形成泥石流堰塞湖可能性高的有17条,占总数的30.9%;可能性中等的有15条,占总数的27.3%;可能性低的有23条,占总数的41.8%(表2)。

图4   潜在泥石流堰塞湖判识结果

Fig.4   The identification results of the potential debris flow barrier lakes

表2   研究区典型一级支沟特征值及潜在泥石流堰塞湖判识结果

Table 2   The characteristic value of first-grade gullies and identification results of the potential debris flow barrier lakes

编号沟名流域面积主沟长度相对高差主沟比降松散物方量支主流量比主河宽度主支夹角判识
结果
(km2)(km)(m)(‰)(×104 m3)(m)(°)
M1郭主铺沟5.705.272093345980.436097
M2万村沟2.853.49164542500.236683
M3新桥沟15.336.2022922722451.356776
M4七盘沟52.6515.5130701539384.3874101
M5磨刀溪沟1.752.691741582340.2510088
M6板桥沟27.1811.1930082103632.387992
M7板子沟54.7317.1540351765811.755296
M8小茅坪沟3.573.291768483160.325677
M9木瓜圆沟3.743.172032509160.4382101
M10大溪沟15.767.6328632737791.6682107
M11马埝坪沟43.3815.22376218514212.188287
M12苏村沟6.035.0023333731630.5812143
M13簇头沟21.898.94293326211081.6289126
M14高店子沟7.355.012434389860.72108101
M15瓦窑坪沟1.203.091825532130.1188100
M16新店沟2.774.002069439500.266086
M17小坪下沟2.042.801615502750.22127102
M18窰子沟7.144.6720603621040.8180101
M19华溪沟10.416.461968239550.8411352
M20大水沟2.142.961664478800.158958
M21石马山7#沟0.671.88133960680.058056
M22石马山6#沟0.491.361170747130.047597
M23石马山5#沟0.361.261144791110.0369112
M24桃关沟50.6914.28300016116442.917573
M25石马山4#沟1.762.201638649350.157563
M26石马山3#沟0.671.911497706150.058099
M27石马山2#沟0.612.201643673180.055992
M28佛堂坝沟33.869.87247919624702.636586
M29石马山1#沟1.302.441659612420.118080
M30彻底关沟16.517.91236223910491.1410076
M31高家沟3.593.2118324857670.3452110
M32罗圈湾沟28.059.62244819820542.3170114
M33张坪岩沟0.671.501437863210.093593
M34银杏坪沟7.004.7220373455050.664878
M35苏坡店上沟0.421.47120773870.04120118
M36苏坡店沟2.843.291633397420.2212084
M37野牛沟24.3810.57301620116551.6510591
M38一碗水沟3.463.361437233900.327245
M39磨子沟7.155.5326523847690.595877
M40下坪沟1.482.662076702470.1410890
M41上坪沟0.682.021634686170.079098
M42马家坪沟2.292.5617826272670.247080
M43太平驿沟26.199.1822902009122.288577
M44清水沟3.854.11159631130.4210895
M45清水店沟39.0413.56320118216602.809770
M46关山沟2.042.7118716211220.566585
M47小磨子沟2.742.8117625023190.536092
M48窝氹沟0.641.541000563150.116079
M49无音寺沟5.344.4716562892060.753875
M50白家林沟0.631.781083548160.094082
M51老街沟0.481.59108761460.064095
M52老街南沟0.511.751091559120.064084
M53白家林南沟0.661.571011559100.154093
M54烧房沟0.651.901117529260.097092
M55红椿沟5.203.7712892742130.797088

新窗口打开

从判识结果可知,研究区形成泥石流堰塞湖可能性高的沟道主要集中在北川映秀断层20 km以内的范围,已经形成堰塞湖的4条泥石流沟(红椿沟、小磨子沟、关山沟和磨子沟)都在此范围内,其主要原因是震后泥石流的所需的松散物源丰富,同时地震造成地表植被破坏,植被拦截降雨的作用减弱,有利于坡面的产、汇流,加之震后沟道的堵塞不仅仅只是沟道弯曲及局部堵塞情况,而是多级的堵塞,多级溃决后的级联效应,增大了泥石流的规模。形成条件的变化导致震后泥石流活动强度、规模的急剧增大,使得泥石流的形成由降雨控制型逐步转为松散土体控制型[18],从而加大了泥石流堰塞湖形成的可能性。

4 结 论

1) 汶川地震诱发的大量崩塌、滑坡为泥石流的形成提供了丰富的松散物质来源,泥石流会在相当长一段时间内处于活跃期,形成泥石流堰塞湖的可能性加大,目前泥石流堰塞湖的危害性已经显现,因此,在今后对泥石流及其次生灾害堰塞湖的危害需加大重视。

2) 震后泥石流的形成条件发生变化,主要表现在松散固体物质的剧增上。通过遥感解译和估算,研究区55条一级支沟内松散堆积物方量约为2.13×108 m3。松散固体物质在研究区内的分布是不均匀的,其主要原因是由于距发震断层北川映秀断裂带远近不同而造成的,距离发震断层越近,其对发震断层的敏感性越大。

3) 通过对潜在泥石流堰塞湖的判识,研究区形成潜在性泥石流堰塞湖可能性高的有17条,主要集中分布在映秀镇至草坡乡段,已经形成堰塞湖的4条泥石流沟都分布在此段,此段将是今后受堰塞湖危害的高危地段,今后应加强此段泥石流沟的防治及监测预警工作。

The authors have declared that no competing interests exist.


参考文献

[1] 崔鹏,韦方强,何思明,.

5·12汶川地震诱发的山地灾害及减灾措施

[J].山地学报,2008,26(3): 280~282.

[本文引用: 1]     

[2] 吴树仁,石菊松,姚鑫,.

四川汶川地震地质灾害活动强度分析评价

[J].地质通报,2008,27(11): 1900~1906.

[本文引用: 1]     

[3] 游勇,柳金峰.

汶川8级地震对岷江上游泥石流灾害防治的影响

[J].四川大学学报(工程科学版),2009,41(增刊):16~22.

[本文引用: 1]     

[4] 常晓军,丁俊,魏伦武,.

岷江上游地质灾害发育分布规律初探

[J].沉积与特提斯地质,2007,27(1):103~108.

[本文引用: 2]     

[5] 柳金峰,游勇.

汶川地震触发潜在性泥石流研究:以岷江上游关山沟为例

[J].四川大学学报(工程科学版),2009,41(增刊):70~75.

[本文引用: 1]     

[6] 胡凯衡,游勇,庄建琦,.

北川地震重灾区泥石流特征与减灾对策

[J].地理科学,2010,30(4):566~570.

[本文引用: 1]     

[7] 徐永年,匡尚富.

泥石流入汇的危险性判别指标

[J].自然灾害学报,2002,11(3):33~38.

[本文引用: 1]     

[8] 郭志学,曹叔尤,刘兴年,.

泥石流堵江影响因素试验研究

[J].水利学报,2004,11:39~45.

[9] 张金山,沈兴菊,谢洪.

泥石流堵河影响因素研究:以岷江上游为例

[J].灾害学,2007,22(2):82~86.

[本文引用: 1]     

[10] 周必凡,李德基,罗德富,.泥石流防治指南[M].北京:科学出版社,1991.

[本文引用: 1]     

[11] 游勇,柳金峰,陈兴长.

“5·12”汶川地震后北川苏保河流域泥石流危害及特征

[J].山地学报, 2010, 28(3): 358~366.

[本文引用: 1]     

[12] 刘希林,陈宜娟.

泥石流风险区划方法及其应用——以四川西部地区为例

[J].地理科学,2010,30(4):558~565.

[本文引用: 1]     

[13] 郝航程,朱方海.

基于可拓学理论的边坡潜在破坏模式识别方法

[J]. 地下空间与工程学报,2007,3(4):698~702.

[本文引用: 1]     

[14] 付静,刘国东,张贞,.

模糊物元模型在土壤侵蚀等级划分中的应用

[J].水土保持通报, 2008, 28(3): 138~140.

[15] 黄辉玲,吴次芳.

基于可拓学的生态市建设评价——以哈尔滨市为例

[J]. 地理科学,2009,29(5):651~657.

[16] 李如忠,汪明武,金菊良.

地下水环境风险的模糊多指标分析方法

[J].地理科学,2010,30(2):229~235.

[17] 徐辉,陈又星.

模糊聚类模型及其在国际市场划分中的应用研究

[J].地理科学,2010,30(3):350~354.

[本文引用: 1]     

[18] 崔鹏,庄建琦,陈兴长,.

汶川地震区震后泥石流活动特征与防治对策

[J].四川大学学报(工程科学版),2010,42(5):10~19.

[本文引用: 1]     

/