19卷 2年 第 1第4期 7月 安全与环境工程 Safety and Environmental Engineering Vo1.19 NO.4 Ju1. 2012 基于变维分形理论的滑坡空间分布特征研究 侯 威,张 茜,楼蓉蓉,彭少波 (中国地质大学工程学院,武汉430074) 摘要:为了对滑坡的空间分布特征进行定量的研究,以香溪河流域高阳一官庄坪段为例,基于变维分形理论对 滑坡空间分布的分形模型进行了优化,定量地得到了各影响因子与滑坡空问分布之间的敏感性关系,同时对滑坡 空间分布的成因进行了定性分析,结果与利用分形理论得到的定量敏感性关系相一致,表明优化的分形模型具有 可靠性。 关键词:滑坡;空间分布;分形特征;变维分形理论;敏感性 中图分类号:X43;P642.21 文献标识码:A 文章编号:1671-1556(2012)04—0017—05 Research on the Characteristics of the Spatial Distribution of Landslides Based on Variable Dimension Fractal Rong—rong,PENG Shao—bo HoU Wei。ZHANG Qian,L0U (Faculty of Engineering,China University of Geosicences,Wuhan 430074,China) Abstract:In order tO make the quantitative research on the characteristics of the spatial distribution of landslides,this paper takes the Gaoyang—Guanzhuangping section of Xiangxi River basin as an example and applies variable dimension fractal theory to the optimization of the fractal mode1.The paper qualitatively analyses the relationship between the sensitivity of spatial distribution and impact factors,and the cause of the spatial distribution of landslides.The results show the accordance of two methods,indicating that the optimized fractal model is reliable. Key words:landslide;spatial distribution;dimension characteristcs;variable dimension fractal theory;sensitivity 0 引 言 目前,对滑坡空间分布特征的研究大多还是定 这种方法来定量地反映滑坡空间分布与各影响因子 之间的敏感性关系。 性研究,但近年来,一些学者开始将变维分形理论应 用到滑坡空间分布特征的定量研究中。如薛天放 1 研究模型 据申维等[4 研究,一般分形统计模型为 y—CX ‘。’ (1) 等口 ]、王志旺L3 通过对滑坡在各尺度下的分布面积 来研究滑坡空间分布的分形特征,使滑坡的空问分 布特征定量化,但是滑坡的分布面积并不能完全反 式中:C为常数;F(D)为一个参数为D的函数,它 与D具有相互对应的关系,所以可以将研究分维数 D转变为研究分维数函数F(D)。 已有研究 q 主要通过对滑坡在各尺度下的分 布面积来研究滑坡空间分布的分形特征,但这并不 映出滑坡的空间分布规律,因为在研究区内不同尺 度的分布面积是不一样的。本文的目的在于找到一 个新的模型来更加准确地反映滑坡空间分布与各影 响因素的敏感性关系,为此以香溪河流域高阳一官 庄坪段为例,基于变维分形理论,用各尺度下滑坡出 现的比例来研究滑坡空间分布的分形特征,试图用 收稿日期:2011一O7—21 修回日期:2011—1O一25 能完全反映出滑坡的空问分布规律,因此本文尝试 通过对各尺度下滑坡出现的比例来研究滑坡的空间 作者简介:侯威(1987一),男,硕士研究生,主要研究方向为岩土工程与稳定性分析。E—mail:hot—hou@hotmail.corn l8 安全与环境工程 第19卷 分布特征。于是得到下面的关系式: P—N /N (2) 式中:P为滑坡在某尺度下出现的比例;N 为滑坡 在某尺度下的面积(km );N 为研究区中某尺度下 的总面积(km )。 其中N 、N 很显然具有变维分形特征,于是套 用一般分形统计模型y—CX ,可以得到: N 一C r l 。 (3) N,一C,rF2‘。 (4) 根据上述公式,可以得到 P—C /C,rFl‘D’F2‘D (5) 设F(D)===F (D)一F (D),于是有 P—Cr 。 (6) 但F(D)在很多时候是一个变量,于是就需要进行 变维分形,将F(D)转换成常分维数。根据何昱华 的 研究,可通过一种累积和的系列变换,使变换后的数据 呈很好的线性关系,这样就可以将变换后的数据用常 维分形来处理。该方法主要有以下几步变换: (1)将原始数据点(N ,/'i)(其中i取1到//'的整 数)经过变换,绘制成双对数曲线。如果得到的曲线 不能很好地拟合成一条直线,就需进行累加变换。对 此可将得到的序列{N }={N , ,…, }进行累 加,可得出如下一阶、二阶以至于 阶的累积和序列: {S }一{Nl,N1+N2,…,N +N2+…+N ) { }==={S ,S +S ,…,S +S +…+S } {S }一{S( 1).,S( 1)_+S( 1)2,…,S( 1)'+ S( 一】)1l+…+S( 】】} (2)将各阶累积和的数据点画在双对数坐标系 中,拟合成一条直线,找出拟合度最高的那组序列的 数据点,于是就可以确定累积和阶次与常分维数。 最后就可以通过最终得到的常分维数来进行分形特 征研究。 2研究区滑坡空间分布的分形特征研究 2.1研究区概况 2.1.1地形地貌条件 研究区内地貌上属鄂西构造侵蚀、剥蚀中低山 区,近南北流向的香溪河是区内最低侵蚀基准面,该 河两岸峰峦叠嶂,河谷形态受岩性的控制,但基本上 呈不对称“V”字型峡谷。香溪河右岸属逆向坡,坡 角随岩性的不同有明显变化,软岩地段一般为20。 左右,厚层砂岩地段在4O。左右;香溪河左岸为顺向 坡,地形坡度大体与岩层倾角一致或稍缓,一般为 4O。~50。。香溪河两岸近东西向支沟较为密集,多 常年流水,均属横向深切沟,两侧沟坡较陡,多呈 “V”字。 2.1.2地层岩性 研究区香溪河流域高阳一官庄坪段主要出露地 层包括:三叠系下统嘉陵江组(T )、中统巴东组 (T 6)和上统沙镇溪组(T。5),侏罗系下统香溪组 (J )、聂家山组(J n)、下沙溪庙组和上沙溪庙组 地层以及各种成因类型的第四系地层。 2.1.3 区域地质构造 在区域地质构造上,研究区属秭归盆地东侧边 缘与黄陵背斜西翼的交接部位,系典型的单斜构造, 岩层总体走向近南北,倾角一般为45。~50。,倾向 西,此构造形式总体控制了香溪河谷的发育及形态, 并构成了以产出于厚层砂岩中的陡倾单面斜坡为主 体的顺向坡。 2.1.4水文地质条件 根据本区地下水的赋存介质,可将地下水划分 为赋存于第四系覆盖层中的子L隙水和赋存于基岩裂 隙中的基岩裂隙水。孔隙水主要接受降雨入渗补给 和部分生活用水的补给,其水位埋深较浅;基岩裂隙 水主要接受降水入渗和上层孔隙水补给,主要以泉 的形式排泄。 2.2滑坡空间分布与各影响因素的分段变维分形 特征 2.2.1 滑坡空间分布与地层的分段变维分形特征 根据研究区内的地质图和野外地质调查得到的 滑坡分布图,运用ArcGIS可得到滑坡的地层分布 图(见图1),并统计出各地层在研究区内的分布面 积和研究区内各地层的滑坡分布面积;然后根据这 两组分布面积,可计算得出滑坡在研究区各地层内 河 图1 滑坡的地层分布图 Fig.1 Strata distribution of landslides 第4期 候威等:基于变维分形理论的滑坡空间分布特征研究 所占的比例,并将所占的比例进行从大到小的排序,用 r来编号,依次取r===1,2,…,6,分别代表不同的地层, 见表1。值得注意的是,没有滑坡的地层,不参加计算, 因为在这些地层中很难产生滑坡,分布特征已经十分 系,相关系数R一0.999 43,于是得到滑坡的空间分 布与地层呈二阶累积和变维分形关系,常分维数 D一1.505 30,见图4。 0 明显,无须再进行变维分形特征研究。 表1滑坡分布与地层的关系表 Table 1 Strata and the distr.hution of landslides 以滑坡在各地层中出现的比例N(r)的对数为 纵坐标,以地层编号r的对数为横坐标作图,可得到 N(r)一r的双对数曲线,见图2。 一1 1 一——1 一1 警一2 一2 —2 0 O.2 O.4 0.6 0-8 log r 图2 滑坡分布比例与地层的原始分段分维序列图 Fig.2 Relationship between the primary dimension series of strata and the distribution ratio of landslides 很显然,所得曲线明显不是一条直线,说明滑坡 分布与地层呈变维分形关系,于是对其进行一阶累 积和变维分形变换,可得到图3。 ~0.7 0.8 0.9 1.0 1.1 l-2 —1.3 图3滑坡分布比例与地层的一阶累积和分段分维序列图 Fig.3 Relationship between the first phase accumulation dimen- sion series of strata and the distribution ratio of landslides 进行一阶累积和变维分形变换后,曲线仍然不 是一条直线,于是进行二阶累积和变维分形变换,在 进行二阶累积和变维分形变换后,数据点呈线性关 0.2 一O.4 0.6 三一0I8 1.0 —1.2 1.4 图4滑坡分布比例与地层的二阶累积和分段分维序列图 Fig.4 Relationship between the second-phase accumulation di— mension series of strala and the distribution ratio of landslides 2.2.2 滑坡空间分布与高程的分段变维分形特征 研究区内高程为80~l 747 1TI,野外调查时三 峡库区水位为160 In,滑坡分布在160~560 In的高 程范围内,于是在160~560 rn中每50 m分段,运 用AcrGIS可得到滑坡的高程分布图(见图5),并对 研究区DEM滑坡空问分布图进行统计分析,可得 到滑坡分布与高程的关系表(见表2)。 香滑一..一I 豳 一_.局一局一向高高一溪坡程程程程程程程程程可 12 2 2 3 3 4 4 5 51 6o m以下 0~260m O~31Om O~360m O~410m 0 ̄460m 0~51Om O~560m 0 m以上 图5滑坡的高程分布图 Fig.5 Elevation distribution of landslides 表2滑坡分布与高程的关系表 Table 2 Elevation and distribution of landslides 局一尚~向一局l 61 6l 620 安全与环境工程 第19卷 在二阶累积和变维分形变换后,数据点呈线性 关系,相关系数R一0.999 65,于是得到滑坡的空间 分布与高程的二阶累积和变维分形关系,常分维数 D—1.449 36,见图6。 0 0 关系,相关系数R一0.999 53,于是可得到滑坡的空 间分布与坡度的二阶累积和变维分形关系,常分维 数D一1.494 31,见图8。 —0 O 0 0 一一溪坡~。:。。香滑 m加如∞如 I .__._ 1 o 0.2 0.4 o.6 0.8 1.0 logr 图6 滑坡分布比例与高程的二阶累积和分段分维序列图 Fig.6 Relationship between the second—phase accumulation dimension series of elevation and the distr|bution ratio of Iandslides 2.2.3 滑坡空间分布与坡度的分段变维分形特征 根据研究区DEM,用ArcGIS提取坡度,对坡 度分级,每隔10。分段,可得到滑坡的坡度分布图 (见图7),并进行统计分析,可得到滑坡分布与坡度 的关系表(见表3)。 图7滑坡的坡度分布图 Fig.7 Slope distribution of landslides 表3滑坡分布与坡度的关系表 Table 3 Slope and distribution of landslides 在二阶累积和变维分形变换后,数据点呈线性 logr 河 。。0。口。图8滑坡分布比例与坡度的二阶累积和分段分维序列图 Fig.8 Relationship between the second—phase accumulation dimension series of slope and the distribution ratio of landslides 2.2.4滑坡空间分布与坡向的分段变维分形特征 根据研究区DEM,采用ArcGIS提取坡向,对 坡向分级,每隔45。分段,分成N、NE、E、SE、S、SW、 W、NW 8个方向,可得到滑坡的坡向分布图(见图 9),并进行统计分析,可得到滑坡分布与坡向的关系 表(见表4)。 溪坡 __ll ____一一一 , 河 图9 滑坡的坡向分布图 Fig.9 Aspect distribution of landslides 表4滑坡分布与坡向的关系表 Table 4 Aspect and distribution of landslides 以 口。上 Ⅳ 第4期 侯威等:基于变维分形理论的滑坡空间分布特征研究 通过一阶累积和变维分形变换后,得到的数据 点呈线性关系,相关系数R一0.994 21,于是可得到 滑坡的空间分布与坡向呈一阶累积和变维分形关 系,常分维数D=0.858 80,见图1O,无须再进行二 阶累积和变维分形变换。 一O. 一0. 一0. 一1. 一1. 一1. 一1. 一1. 一1. 一1. 1ogr 图lO 滑坡分布比例与坡向的一阶累积和分段分维序列图 Fig.10 Relationship between the first-phase accumulation dimension series of aspect and the distribution ratio of landslides 3研究区滑坡空间分布的成因分析 3.1 坡向对滑坡空间分布规律的影响 河流附近的滑坡绝大多数都是朝向河流或者沟 谷方向,在香溪河流域,由于香溪河大体上是南北走 向,岩层倾向大体是西向,两组主要节理倾向为东向 和南向,而这两组节理的走向正好与香溪河、冲沟以 及长江的走向相同,可见这两组节理决定了香溪河、 冲沟以及长江的走向,所以也间接地决定了滑坡的 坡向,因此坡向造成滑坡的分布很有规律,大多都朝 向河流或冲沟。 3.2 地层对滑坡空间分布规律的影响 在J。一 地层中产生滑坡的比例最高,为 5.759 4 ,在J 和J。X地层中出现滑坡的比例也较 高,分别为4.743 6 和4.339 8 ,而从这三组地层的 岩性来看,多为黏土质粉砂岩,较软弱,且在这三组地 层之上多为石英砂岩,之下多为灰岩白云岩,都较为 坚硬,所以滑坡集中分布在这三组地层中。 3.3坡度对滑坡空间分布规律的影响 据研究,由于滑坡过程是由陡坡变为缓坡的能量 释放过程,所以滑动坡体的总体坡度较周围山体平缓, 有的甚至成为平坦地形或凹地【6]。在三峡库区,大多 数滑坡也都是缓倾角滑坡,在香溪河流域也不例外,在 坡度为10。~2O。的区域内出现滑坡的比例最高,为 3.347 ,在坡度为20。~3Oo、10。以下的区域内出现滑坡 的比例也较高,分别为2.693 4%和2.320 5%。 3.4 高程对滑坡空间分布规律的影响 在香溪河流域,高程与到河流的距离是相互联 系的,高程越低,离河流的距离也越近。三峡大坝蓄 水后,水位在135~175 m之间波动,研究区内高程 在210 m以下出现滑坡的比例最高,高程在210~ 310 m区域内出现滑坡的比例也较高。由于水的作 用,离河流越近,滑坡的分布越密集。 4结论 (1)研究区滑坡的空间分布与坡向具有一阶累 积和变维分形分布特征,与地层、坡度、高程均具有二 阶累积和变维分形分布特征,表明滑坡的空间分布对 坡向的敏感性最高,要明显大于岩性、坡度与高程;滑 坡空间分布与地层、坡度、高程的常分维数D分别为 1.505 30、1.494 31、1.449 36,表明滑坡的空间分布对 地层的敏感性最高,坡度次之,高程较低。 (2)通过对研究区滑坡空间分布的成因进行定 性分析可以发现,滑坡空间分布对坡向、地层、坡度、 高程的敏感性关系与用变维分形理论得到的敏感性 关系是相对应的,说明用变维分形理论来研究滑坡 的空间分布特征是可行、可靠的,且可以定量地反映 滑坡的空间分布与各影响因素之间的敏感性关系。 (3)利用此次定量研究得出的各影响因子与滑 坡空间分布之间的敏感性关系,可为基于GIS的滑 坡空间分布预测打下基础。即在采用GIS预测滑 坡的空间分布时,可通过确定各影响因子的权重,从 而找到一套滑坡空间分布预测方法,并通过某一区 域内的滑坡资料,提取定量化的滑坡空间分布特征, 然后预测出整个研究区域内的滑坡空间分布。 (4)本次对滑坡空间分布特征的分形研究是将 分形理论以及变维分形理论应用在地质灾害方面的 一种尝试和对前人研究滑坡空间分布的分形特征的 一种优化,相信该研究在地质灾害研究中将会有广 阔的应用前景。 参考文献: [1]薛天放,杨庆,栾茂田.基于GIS技术的滑坡空问分布的分形特 征研究[J].岩土力学,2007,28(2):347—351,354. [2]薛天放.基于地理信息系统技术(GIS)的滑坡空问分布的分形特 征研究[D].大连:大连理工大学,硕士学位论文,2005. [3]王志旺.基于GIS技术的区域滑坡分形特征分析与危险性评价 [D].武汉:中国地质大学,博士学位论文,2010. [4]申维,赵鹏大.分形统计模型的理论研究及其在地质学中的应用 [J].地球科学,1998,33(2):234--243. [5]付昱华.变换形成的分形与海洋环境数据分析预测[J].海洋通 报,2000,19(1):79—88. [6]唐辉明.工程地质学基础[M].北京:化学工业出版社,2008:98.