云南省滑坡災(zāi)害分形特征研究與分形評價模型探討

史興旺，管新邦，邊 筠

(煤炭科學技術(shù)研究院有限公司安全分院，北京 100013)

分形理論在分析復雜的非線性系統(tǒng)方面取得了良好的效果[1]。例如Sivakumar B.[2]論證了運用多重分形理論框架解釋不同氣候帶降雨特征的可行性，Kotowski P.[3]提出了一套完整的方法來測量鐵基合金表面斷裂的分形維數(shù)，并證明所測斷裂面的分形維數(shù)不依賴于該研究對象的材料，Zhang Y等[4]將分形理論運用到城市地理中，來反映和模擬城市邊界形態(tài)。對于滑坡災(zāi)害來說，因其受到多種因素的綜合作用和影響，滑坡地表變形演化過程往往具有復雜的非線性特征[5]，并且隨著荷載值增大整體呈現(xiàn)非線性增大的趨勢[6]。這為適用于解決復雜非線性問題的分形理論找到了合適的切入點。近年來，分形理論在滑坡研究領(lǐng)域也受到了部分學者的關(guān)注和認可。例如Katz O.等[7]通過模型實驗探討了邊坡破壞類型和破壞頻率之間的分形關(guān)系，Li C等[8]分析了滑坡發(fā)生和降雨之間的冪指數(shù)關(guān)系。此外，針對滑坡位移趨勢研究，鄭明新等[9]分析了滑坡演化過程中滑坡變形曲線的分維特征，并首次提出了滑坡大滑分維預(yù)報公式。樊曉一等[10]對不同穩(wěn)定狀態(tài)下地震滑坡的地形剖面線的多重分形維數(shù)演化特征進行了研究，并表明相對穩(wěn)定滑坡的多重分形維數(shù)演化特征具有明顯的有序性。

綜合前人研究，本文以分形理論為基礎(chǔ)，對研究區(qū)滑坡的分形特征進行研究，獲得不為整數(shù)的分形維數(shù)，進一步得到基于分維數(shù)的滑坡災(zāi)害對不同影響因素的相關(guān)性，與邏輯回歸模型、證據(jù)權(quán)模型、證據(jù)可信度函數(shù)模型和信息量模型所得到的相關(guān)性具有較高的一致性?；碌膸缀涡螒B(tài)及空間分布具有復雜的不確定性，不同的滑坡類型有不同的形態(tài)特征，幾何形態(tài)特征和空間分布對于滑坡的研究有著重要的意義，分形理論可對建立滑坡災(zāi)害易發(fā)性評價模型提供新的科學模型。

1 研究方法

世界是非線性的，自然界中分形無處不在。分形理論是用分數(shù)維度的視角和數(shù)學方法描述和研究客觀事物，該理論認為局部特征形態(tài)與整體特征形態(tài)具有自相似的特征。分形理論突破了從傳統(tǒng)的點線面和時空的研究方法，是一門以不規(guī)則幾何形態(tài)為研究對象的新幾何學。分形理論已經(jīng)應(yīng)用到圖像處理、自然科學、管理學、工程技術(shù)、經(jīng)濟學和文化藝術(shù)等各個方面。分維是分形理論中的基本參數(shù)和重要原則。研究滑坡的分維值可以了解其復雜程度和演化規(guī)律?；禄顒訌姸鹊纳邔⒊尸F(xiàn)出其分維值的升高，分維值較高時，則面臨較高的產(chǎn)生滑坡的可能性?；庐a(chǎn)生后其分維值降低，可以根據(jù)某個時期內(nèi)斜坡系統(tǒng)分維值的變化及其升幅的大小來預(yù)測滑坡產(chǎn)生的數(shù)量[11]。

本文采用遙感數(shù)據(jù)為高分一號遙感數(shù)據(jù)，解譯識別的滑坡災(zāi)害數(shù)據(jù)庫分析。通過對遙感圖像進行正射校正、圖像融合、圖像增強處理和圖像鑲嵌與裁剪以及配合大量的野外調(diào)查，特別是小比例尺的地質(zhì)災(zāi)害詳查完成解譯。結(jié)合云南省圖件資料的整理和對云南省滑坡編錄資料的分析，利用盒維數(shù)法對研究區(qū)的滑坡災(zāi)害分形特征進行分析研究。

(1)常維分形可用公式表示為：

N=Cr-D.

(1)

其中，N為特征尺度為r的物體數(shù)目，r為特征尺度，C為常數(shù)，D為分維數(shù)。

其定義為：設(shè)A屬于RN，在歐氏距離下，統(tǒng)計用邊長1/2n的盒子個數(shù)Nn(A)，則有：

(2)

當n增大時，統(tǒng)計計算Nn(A)，隨著一組系列的r1，r2，……rn，得到一組N(r1)，N(r2)，……N(rn)，以點[log(r)，logN(r)]為坐標作雙對數(shù)圖，利用最小二乘法，擬合成一條直線得出斜率，其絕對值即為分形維數(shù)。其直線公式為：

lnNn(A)=a+b1/ln2n.

(3)

其中a為常數(shù)，b為分形維數(shù)。

(2)變維分形其函數(shù)關(guān)系表示為：

D=g(r).

(4)

對于N與r之間的函數(shù)關(guān)系N=f(r)，令f(r)=Cr-D，則有

(5)

有研究表明，任一參數(shù)A隨重現(xiàn)期的變化規(guī)律均可以轉(zhuǎn)化為變維分形的形式。若將常維分形看作1階分形，若其中的常維分形為D=C′/rD′，則可以對1階分形再進行一次分形，以此類推可得到N階分形[12]。其表達式為：

N1=C1/rD1.

(6)

其中D1=C2/rD2；D2為常數(shù)。

N2=C1/rD1.

(7)

……

N階分形，其中Dn-1=Cn/rDn；Dn為常數(shù)。

Nn=Cn-1/rDn-1.

(8)

由此可見，對于任意階的變維分形，都可以轉(zhuǎn)化為常維分形的形式。多階分形可以通過階數(shù)累計處理，衍生出累計和變化分形。具體步驟為：

1)將原始數(shù)對(Ni，ri)按照從小到大的序列投影到雙對數(shù)坐標系中，計算相鄰兩點間在雙對數(shù)坐標系中的Di,i+1。Di,i+1的求解公式為：

Di,i+1=ln(Ni/Ni+1)ln(ri/ri+1).

(9)

2)對(N1,N2,L,Ni,L,Nn)進行累計和計算。

{S1i}={N1,N1+N2,N1+N2+N3,…};i=1,2,3,…,n.

(10)

{S2i}={S11,S11+S12,S11+S12+S13,…};ji=1,2,3,…,n.

(11)

{S3i}={S21,S21+S22,S21+S22+S23,…};i=1,2,3,…,n.

(12)

{Sni}={S(n-1)1,S(n-1)1+S(n-1)2,S(n-1)1+S(n-1)2+S(n-1)3,…};i=1,2,3,…,n.

(13)

其中，S1，S2，S3，…Sn為相應(yīng)分形階數(shù)的累計和序列。

3)根據(jù)以上公式，由S1i可以得到在雙對數(shù)坐標系中的n-1條斜率不等的線段，這些斜率的相反數(shù)則為一階累計和變維分形維數(shù)，同理可以得到高于一階累計和變維分形維數(shù)。隨著分形階數(shù)的增加，雙對數(shù)坐標系中的各階累計和形成的數(shù)據(jù)點最終會呈現(xiàn)出很好的線性關(guān)系。此時可以用線性回歸模型計算該線性關(guān)系的參數(shù)，得到其常分維數(shù)D。

2 常維分形特征

常維分形是指分維數(shù)D不隨特征尺度r的變化而變化，二者之間不存在某種函數(shù)關(guān)系。在不考慮滑坡地質(zhì)災(zāi)害分布地區(qū)的地貌地形特征，把滑坡看作是二維平面上的點，運用盒維數(shù)法來研究常維分形。盒維數(shù)法是一個常用而又簡單的方法，在操作過程中利用公式不斷縮小盒子的尺寸，建立尺寸為2n的灰度圖像，與災(zāi)害點圖像覆蓋，得到各尺寸盒子與點的相交程度，然后在對數(shù)坐標系中建立相交盒子個數(shù)N與盒子尺寸r的散點圖，采用最小二乘法進行線性擬合，驗證其擬合程度的相關(guān)系數(shù)，所得擬合直線的斜率即為待求圖形的分形維數(shù)。

在滑坡地質(zhì)災(zāi)害的空間分布關(guān)系中，同樣可以用盒維數(shù)法進行分形特征的研究。圖1為滑坡空間地質(zhì)災(zāi)害的雙對數(shù)模型，線性回歸方程為：y=-0.8432x+6.8836，R2=0.9906。其擬合程度接近于1，說明云南省的滑坡地質(zhì)災(zāi)害在空間分布上具有明顯的分形特征。可以看出研究區(qū)的滑坡分布受制于一定的分形維數(shù)，又因為滑坡產(chǎn)生后其分維值會降低，我們可以求解不同時期的研究區(qū)的滑坡分維值，統(tǒng)計其變化趨勢特征，預(yù)測研究區(qū)滑坡是否會發(fā)生，同時可以根據(jù)分維降低的幅度，預(yù)測滑坡數(shù)量和規(guī)模。

圖1 云南省滑坡點空間分形特征Fig.1 Yunnan landslide spatial fractal features

云南省的水分、地形與構(gòu)造等條件的組合對滑坡災(zāi)害的頻繁發(fā)生提供了便利，從而表現(xiàn)出滑坡空間維數(shù)值D=0.8432，高于中國大陸及三大分區(qū)滑坡空間分維值[13]，說明云南省滑坡災(zāi)害在全國范圍來講屬于高發(fā)區(qū)。

基于“過去和現(xiàn)在是打開將來之門的鑰匙”理論，過去和現(xiàn)在的滑坡的空間分布，將來也將是最有可能發(fā)生滑坡的地區(qū)[14]。因此，在進行滑坡易發(fā)區(qū)評價時，可以運用滑坡的分形公式進行擬合，計算滑坡周圍不同距離內(nèi)發(fā)生滑坡的概率。因此，在進行滑坡相關(guān)分析中，應(yīng)研究不同區(qū)域的分形特征，對比不同地區(qū)滑坡災(zāi)害的分維值，建立分維值分區(qū)，對分維值大的因素要重點考慮。

3 滑坡影響因子與空間分布變維分形特征

研究滑坡的分布與其影響因子的分形關(guān)系，可以得到滑坡對不同影響因素的相關(guān)性，進而為易發(fā)性的評價和分區(qū)提供基礎(chǔ)資料。然而在這類相關(guān)研究中，嚴格滿足常維分形的事件是沒有的，這時常維分形則常常不能滿足工作需要。在分形中，分維數(shù)D不是常數(shù)，而是與特征尺度r之間存在著某種函數(shù)關(guān)系時，即為變維分形。

3.1 滑坡面積與坡度的分形特征

坡度是控制滑坡發(fā)育的重要因素之一。按照坡度步長為10°，對滑坡面積和與坡度區(qū)間進行統(tǒng)計，見表1所示?？梢钥闯鲭S著坡度的增加，圖形出現(xiàn)先增后減的趨勢，在20°到30°之間出現(xiàn)峰值。對數(shù)據(jù)進行雙對數(shù)變換處理，發(fā)現(xiàn)并不呈現(xiàn)出常維分形的特征；對數(shù)據(jù)進行高階累計和變換處理，可以看出滑坡面積的空間分布與坡度在4階累積和曲線上呈現(xiàn)很好的線性關(guān)系，見圖2。對4階累計和變換數(shù)據(jù)進行線性回歸擬合，R2=0.9968，擬合相關(guān)系數(shù)接近于1，這表明滑坡面積空間分布與坡度呈4階累積和變維分形關(guān)系，其分維值D4=2.1384。

表1 坡度區(qū)間與滑坡面積關(guān)系統(tǒng)計表Table 1 Interval slope and landslide area relationship tables

圖2 滑坡面積與坡度分維數(shù)分形特征Fig.2 Landslide area and slope fractal dimension characteristic curve

3.2 滑坡面積與高程的分形特征

高程的增加會影響滑坡地質(zhì)災(zāi)害的發(fā)生概率和規(guī)模大小。按照高程步長為500 m，對滑坡面積和與高程區(qū)間進行分等級統(tǒng)計，大于3 500 m為一個等級，見表2所示。可以看出，隨著高程的增加滑坡面積與高程之間并不存在明顯的線性關(guān)系。對滑坡面積與高程之間進行高階累計和變維分形變換，在3階累計和與特征尺度r的雙對數(shù)模型中具有很好的線性關(guān)系，見圖3。表明滑坡面積與高程之間存在著3階累計和變維分形關(guān)系，其分維值D3=3.5059，相關(guān)系數(shù)R2=0.9961，擬合相關(guān)系數(shù)接近于1。

表2 高程區(qū)間與滑坡面積關(guān)系統(tǒng)計表Table 2 Elevation range and landslide area relationship tables

圖3 滑坡面積與高程分維數(shù)分形特征Fig.3 Landslide area and elevation fractal dimension characteristic curve

3.3 滑坡面積與斷裂緩沖區(qū)之間的分形特征

在ArcGIS中對斷裂緩沖區(qū)設(shè)計步長為1 km、緩沖半徑為10 km的緩沖區(qū)，統(tǒng)計各緩沖區(qū)內(nèi)滑坡地質(zhì)災(zāi)害的面積，見表3所示?？梢钥闯鲈陔p對數(shù)關(guān)系中，滑坡面積與斷裂緩沖區(qū)之間不存在常維分形特征。在高階累計和變維分形變換中，滑坡面積與斷裂緩沖區(qū)在2階累計和雙對數(shù)模型中的線性關(guān)系要好于3階累計和模型的線性關(guān)系，見圖4。其線性回歸曲線擬合程度R2=0.9992，擬合相關(guān)系數(shù)接近于1。表明滑坡面積與斷裂緩沖區(qū)存在2階累計和變維分形關(guān)系，其分維值D2=1.6005。

表3 斷裂緩沖區(qū)間與滑坡面積關(guān)系統(tǒng)計表Table 3 Fracture buffer range and landslide area relationship tables

圖4 滑坡面積與斷裂緩沖區(qū)分維數(shù)分形特征Fig.4 Landslide area and fracture buffer range fractal dimension characteristic curve

3.4 滑坡面積與巖性之間的分形特征

根據(jù)巖性的軟硬程度屬性，運用數(shù)字代碼建立滑坡面積與巖性之間的關(guān)系，見表4所示?？梢钥闯鲈诟唠A累計和變維分形中，滑坡面積與巖性在2階累計和雙對數(shù)坐標系中存在著最優(yōu)的線性關(guān)系，見圖5，其擬合相關(guān)系數(shù)R2=0.9899，擬合相關(guān)系數(shù)接近于1。表明滑坡面積與巖性之間存著在2階累計和變維分形，其分維數(shù)D2=2.3876。

表4 巖性軟硬等級與滑坡面積關(guān)系統(tǒng)計表Table 4 Lithology grade and landslide area relationship tables

圖5 滑坡面積與巖性分維數(shù)分形特征Fig.5 Landslide area and lithology grade fractal dimension characteristic curve

4 結(jié)果分析

以上研究結(jié)果表明，云南省滑坡地質(zhì)災(zāi)害在空間上具有明顯的分形特征。云南省滑坡地質(zhì)災(zāi)害的面積空間分布與坡度、高程、斷裂緩沖區(qū)和巖性等環(huán)境因子呈現(xiàn)不同階的變維分形特征：與斷裂緩沖區(qū)和巖性呈現(xiàn)2階累計和變維分形，與高程呈現(xiàn)3階累計和變維分形，與坡度呈4階累計和變維分形。分維值分別為坡度D4=2.1384，高程D3=3.5059，斷裂緩沖區(qū)D2=1.6005，巖性D2=2.3876。

為了研究滑坡災(zāi)害的分布與各環(huán)境因子之間的變維分形的關(guān)系，分別采用了邏輯回歸模型[15-17]、證據(jù)權(quán)模型[18-19]、證據(jù)可信度函數(shù)模型[20-21]和信息量模型[22]對研究區(qū)的滑坡災(zāi)害與影響因素的關(guān)系進行了統(tǒng)計分析，分別得到了回歸系數(shù)、C值、Bel值和信息量值。這些量值分別代表了這4種評價模型所構(gòu)建的影響因素在滑坡災(zāi)害事故發(fā)生中的權(quán)重，其量值的大小即為影響因素權(quán)重值的大小。對比發(fā)現(xiàn)分形維數(shù)與以上四種方法得到的權(quán)重值在評價滑坡災(zāi)害分布上表現(xiàn)出共同的相關(guān)性，即分維值越大環(huán)境因子對滑坡地質(zhì)災(zāi)害的影響越大，見表5。因此推斷分維數(shù)可以用于表征滑坡地質(zhì)災(zāi)害與環(huán)境因子之間在分形上的相關(guān)程度，分維數(shù)越大，則環(huán)境因子對滑坡地質(zhì)災(zāi)害的影響越大，其敏感度排序為：高程>地層巖性>坡度>斷裂緩沖區(qū)。

表5 不同評價模型的影響因素敏感性分級Table 5 Different factors evaluation model sensitivity rating

可見分形理論得出的滑坡災(zāi)害對不同影響因素的敏感性與邏輯回歸模型、證據(jù)權(quán)模型、證據(jù)可信度函數(shù)模型、信息量模型所得到的結(jié)論具有很高的一致性，可以依此建立基于分形理論的滑坡災(zāi)害易發(fā)性評價模型。但是，由于分維模型的建立及分維數(shù)的求取與影響因子分類分級數(shù)目、分級方式和提取單元有很大的關(guān)系。不同精度的數(shù)據(jù)、不同的分級方法求取的計算結(jié)果可能會有較大差異，不同的區(qū)域由于環(huán)境特點不同，其結(jié)果也可能有很大差別。因此在建立分形評價模型時要嚴謹考慮這些問題。

5 結(jié)論

(1)云南省滑坡地質(zhì)災(zāi)害的面積空間分布與坡度、高程、斷裂緩沖區(qū)和巖性等環(huán)境因子呈現(xiàn)不同階的變維分形特征，與斷裂緩沖區(qū)和巖性呈2階累計和變維分形，與高程呈3階累計和變維分形，與坡度呈4階累計和變維分形。

(2)分維數(shù)可以反映滑坡地質(zhì)災(zāi)害與環(huán)境因子之間在分形上的復雜程度，其特征是分維數(shù)越大，則環(huán)境因子對滑坡地質(zhì)災(zāi)害的影響越大?；碌挠绊懸蜃用舾卸鹊呐判蚪Y(jié)果為：高程>地層巖性>坡度>斷裂緩沖區(qū)。

(3)分形理論得出的滑坡災(zāi)害對不同影響因素的敏感性與邏輯回歸模型、證據(jù)權(quán)模型、證據(jù)可信度函數(shù)模型和信息量模型所得到的結(jié)論具有很高的一致性，可以依此建立基于分形理論的滑坡災(zāi)害易發(fā)性評價模型。