基于層次分析法的德州地面沉降易發(fā)區(qū)的可拓學(xué)劃分應(yīng)用

張永偉,李紅霞,商婷婷，賈超,梁浩,楊亞賓

(1.山東省地質(zhì)環(huán)境監(jiān)測總站，山東 濟(jì)南　250014；2.山東大學(xué)，山東 濟(jì)南　250061；3.濟(jì)南市育賢小學(xué)，山東 濟(jì)南　250024；4.山東省魯北地質(zhì)工程勘察院，山東 德州　253015)

?

技術(shù)方法

基于層次分析法的德州地面沉降易發(fā)區(qū)的可拓學(xué)劃分應(yīng)用

張永偉1，2,李紅霞3,商婷婷1，賈超2,梁浩1,楊亞賓4

(1.山東省地質(zhì)環(huán)境監(jiān)測總站，山東 濟(jì)南250014；2.山東大學(xué)，山東 濟(jì)南250061；3.濟(jì)南市育賢小學(xué)，山東 濟(jì)南250024；4.山東省魯北地質(zhì)工程勘察院，山東 德州253015)

摘要：基于物元模型的可拓學(xué)理論分析，提出華北平原德州地面沉降易發(fā)區(qū)的可拓學(xué)劃分方法。利用地面沉降易發(fā)性分級、分類標(biāo)準(zhǔn)和影響因子分析，構(gòu)造出經(jīng)典域物元和節(jié)域物元，應(yīng)用物元和可拓集合中的關(guān)聯(lián)函數(shù)，建立了易發(fā)性等級綜合評判的可拓評價模型，通過基于層次分析法的可拓學(xué)評價分析，采用正方形等間距(2km×2km)剖分了135個網(wǎng)格，再將剖分單元轉(zhuǎn)化成面元，經(jīng)過等值差分，得到了地面沉降易發(fā)區(qū)劃結(jié)果。計算結(jié)果顯示，可拓學(xué)方法能夠?qū)崿F(xiàn)定量化、多因子評價地面沉降易發(fā)性等級問題，從而科學(xué)合理地指導(dǎo)地面沉降綜合分區(qū)防治。

關(guān)鍵詞：地面沉降；可拓學(xué)；物元模型；易發(fā)區(qū)劃；層次分析法

引文格式：張永偉,李紅霞,賈超,等.基于層次分析法的德州地面沉降易發(fā)區(qū)的可拓學(xué)劃分應(yīng)用[J].山東國土資源，2016,32(7)：58-64.ZHANG Yongwei，LI Hongxia，SHANG Tingting，etc. Application of Extenincs Theory for Partition of Land Subsidence Susceptible Area of Dezhou City Based on Analytic Hierarchy Process[J].Shandong Land and Resources, 2016,32(7)：58-64.

地面沉降是一種致使地面高程降低的緩變性地質(zhì)災(zāi)害[1]。基于層次分析的可拓學(xué)是較為科學(xué)的定量化評價方法，采用這種方法對華北平原地面沉降最為嚴(yán)重的地區(qū)之一德州進(jìn)行易發(fā)區(qū)劃分，對于該區(qū)的地面沉降綜合防治具有一定的指導(dǎo)作用。

1研究區(qū)概況

德州市位于黃河下游沖積平原區(qū)，1965年深層地下水開采以來，深層地下水位降落漏斗迅速在縱、橫向上拓展，漏斗中心水位埋深已超過140m。20世紀(jì)80年代，德州市范圍內(nèi)產(chǎn)生了不同程度的地面沉降。德州是華北平原地面沉降出現(xiàn)最早、影響最大和最為嚴(yán)重的地區(qū)之一，已與河北省形成連片沉降區(qū)[2]。

德州市地面沉降各監(jiān)測點(diǎn)以每年16.9～62.5mm(1991—2010年)的速率下沉。沉降中心位于城區(qū)的國棉一廠院內(nèi)(D62)，19年累計沉降量為-1186.9mm，年均沉降量為-62.5mm。城區(qū)邊緣的木材公司院內(nèi)(D9)累計沉降量(19年)為-636.9mm，年均沉降量(19年平均值)-33.5mm。外圍(YD1)累計沉降量為-321.7mm，年均沉降量(19年平均值)-16.9mm(表1)。由此看出，監(jiān)測區(qū)各觀測點(diǎn)存在明顯的不均勻沉降。隨著觀測點(diǎn)不均勻沉降的發(fā)生，地面點(diǎn)高差會越來越大[3]。

2可拓學(xué)模型

2.1可拓學(xué)簡介

可拓學(xué)(extenics)是一門橫斷學(xué)科和交叉學(xué)科，由廣東工業(yè)大學(xué)蔡文為首的中國學(xué)者創(chuàng)立，是研究事物拓展的可能性和開拓創(chuàng)新的規(guī)律與方法[4]。目前，可拓學(xué)應(yīng)用于工業(yè)[5]、醫(yī)學(xué)、農(nóng)業(yè)、土地等級評價等領(lǐng)域，在地質(zhì)學(xué)上應(yīng)用如圍巖質(zhì)量分類[6]、泥石流危險度評價、地質(zhì)災(zāi)害脆弱性評價和損毀山體影響度評價等[7]。

假設(shè)事物的等級為N，它關(guān)于特征c的量值為v，以此三元組作為表述事物的基本元，為物元R構(gòu)成的三要素，R=(N，c，v)。

可拓學(xué)應(yīng)用的基本思路為：將評價目標(biāo)根據(jù)數(shù)據(jù)值進(jìn)行等級劃分，構(gòu)建可拓模型，再將評價指標(biāo)依次代入各等級的集合中進(jìn)行多指標(biāo)評定，根據(jù)評定結(jié)果與各等級集合的關(guān)聯(lián)度大小進(jìn)行比較來確定評定等級。

表1　德州市1991—2010年各監(jiān)測點(diǎn)高程及沉降量統(tǒng)計

2.2可拓學(xué)模型建立

2.2.1確定經(jīng)典域

(1)

式中：N0j表示所劃分的第j(j=1，2，3，…，m)個等級，ci(i=1，2，3，…，n)表示易發(fā)性等級N0j的因素，Voji=﹤aoji，boji﹥?yōu)镹0j關(guān)于因素ci所確定的量值范圍——經(jīng)典域。

2.2.2確定節(jié)域

(2)

式中：P為等級的全體Vpi=﹤api，bpi﹥,為P關(guān)于因素i所取量值的范圍，即P的節(jié)域。

2.2.3確定待評物元

(3)

式中：p為待評價目標(biāo)，ci為影響易發(fā)等級的因素，vi為p關(guān)于因素ci的量值。

2.2.4確定評價指標(biāo)關(guān)于各等級的關(guān)聯(lián)度

各單項評價指標(biāo)Vi關(guān)于各等級j的關(guān)聯(lián)度為

(4)

2.2.5待評物元關(guān)于各類別等級的關(guān)聯(lián)度

待評目標(biāo)p關(guān)于評價等級j的關(guān)聯(lián)度為：

(5)

2.2.6可拓評判等級

(6)

(7)

式中：j*為P的級別變量特征值，從j*中可看出偏向另一類的程度。

3地面沉降易發(fā)性等級可拓評價

3.1評價體系的建立

3.1.1評價因子的確定

在綜合分析地面沉降影響因素的基礎(chǔ)上，參考有關(guān)研究成果[8-10]，選擇地面沉降狀況因子和地質(zhì)環(huán)境條件因子(包括地下水開采狀況)作為評價因子。地面沉降狀況因子包括累計沉降量、沉降速率和預(yù)測沉降量[11-14]。地質(zhì)環(huán)境條件因子包括水位允許降深、新生界厚度、軟土層厚度和地面高程。

3.1.2評價指標(biāo)體系

(1)體系的層次構(gòu)成。評價體系由評價目標(biāo)層、影響因子層(準(zhǔn)則層)、要素屬性層和要素指標(biāo)層構(gòu)成見表2。

(2)評價指標(biāo)標(biāo)準(zhǔn)：易發(fā)程度高為Ⅰ級；易發(fā)程度中等為Ⅱ級；易發(fā)程度低為Ⅲ級。

3.2確定評價區(qū)并制作基礎(chǔ)網(wǎng)格

該次評價范圍確定為德州市德城區(qū)，面積為539km2。該次評價采用正方形等間距網(wǎng)格進(jìn)行剖分，網(wǎng)格大小為2km×2km，共剖分為135個網(wǎng)格。然后將剖分單元轉(zhuǎn)化成面元，并對每一剖分單元進(jìn)行連續(xù)編號，將單元編號作為各面元的“ID”屬性賦值。將單元的各種指標(biāo)值分別賦予各單元的屬性值，建立單元屬性庫。

3.3單元數(shù)據(jù)的獲取

根據(jù)野外實測值、觀測值、地形圖、地質(zhì)圖量算值等方式獲取單元數(shù)據(jù)。

3.4權(quán)值的確定

權(quán)值的確定方法采用層次分析法(Analytic Hierarchy Process，簡稱AHP)[15]。層次分析定權(quán)法判斷矩陣標(biāo)度分級依據(jù)表3。權(quán)值計算采取自上而下(目標(biāo)層A—因子層B—屬性層C)的方法。

表2　評價指標(biāo)體系與指標(biāo)等級劃分

表3　層次分析定權(quán)法的判斷矩陣標(biāo)度分級及其含義

3.4.1構(gòu)建目標(biāo)層A的綜合判斷矩陣

為實現(xiàn)B層對目標(biāo)層A層的描述，對B層指標(biāo)進(jìn)行相對重要性兩兩比較，確定地面沉降狀況(B1)和地質(zhì)環(huán)境條件(B2)的組合群對總目標(biāo)的權(quán)重關(guān)系。目標(biāo)層A的判斷矩陣為：A={bij│i，j=1-2}2×2在上述工作的基礎(chǔ)上，構(gòu)建的地面沉降危害性評價權(quán)重綜合判斷矩陣見表4。

表4　地面沉降危害性評價權(quán)重綜合判斷矩陣

3.4.2構(gòu)建2個因子層B的判斷矩陣

分別根據(jù)因子層(B1、B2)的要求構(gòu)建B層和C層的判斷矩陣(表5、表6)。

表5　地面沉降狀況因子相對重要性判斷矩陣

表6　地質(zhì)環(huán)境狀況因子相對重要性判斷矩陣

3.4.3層次排序及求解權(quán)向量

根據(jù)判斷矩陣，利用線性代數(shù)，求出矩陣的最大特征根所對應(yīng)的特征向量。所求特征向量即為各評價因素的重要性排序，歸一化后，也就是權(quán)數(shù)分配。

對于一個正向量W=(W1，W2，…，Wn)T，其標(biāo)準(zhǔn)化向量W為：

(8)

式中：Ws(Ws1，Ws2，…，Wsn)T為同一層次相應(yīng)因子對于上一層次某個指標(biāo)相對重要性的排序權(quán)值。

目前常見的計算因子排序權(quán)重向量的方法主要有方根法、和積法、特征根法、最小二乘法和對數(shù)最小二乘法等。鑒于方根法計算方便，該次采用方根法近似求解，并用Excel進(jìn)行計算[16]。計算結(jié)果見表7。

表7　矩陣權(quán)值計算成果

3.4.4判斷矩陣一致性檢驗

單層排序權(quán)值是否合理，需進(jìn)行一致性檢驗，公式為：

(9)

(10)

式中：CR為一致性比例；RI為平均一致性指標(biāo)，其數(shù)據(jù)如表8所示。

CI為一致性指標(biāo)。當(dāng)CR<0.1時，達(dá)到一致性；否則，需重新調(diào)整。判斷矩陣一致性檢驗結(jié)果如表9所示，CR值均<0.1，判斷矩陣均具有滿意的一致性。

表8　層次分析法的平均隨機(jī)一致性指標(biāo)值

表9　矩陣一致性檢驗結(jié)果

3.4.5確定層次總排序權(quán)值

目標(biāo)層A包含2個因子B1和B2，其權(quán)值分別為a1和a2；B層包含7個因子C1，C2，C3，...，Ck，...，C7；它們對因子Bj的層次排序權(quán)值分別為b1j，b2j，...，bnj(當(dāng)Ck與Bj無關(guān)時，bnj=0)。C層總排序權(quán)值為二層權(quán)值乘積代數(shù)和。計算公式為：

(11)

計算結(jié)果見表10。

表10　評價體系權(quán)值計算結(jié)果

3.5可拓評價

根據(jù)式(6)計算出每個單元網(wǎng)格的K值，確定每個單位網(wǎng)格的可拓評判等級，分高、中、低三級評價其易發(fā)程度，最后根據(jù)所有單元網(wǎng)格可拓等級，并綜合考慮地質(zhì)環(huán)境條件和地面沉降狀況綜合確定地面沉降高、中、低易發(fā)區(qū)。評價過程簡述如下。

3.5.1物元的構(gòu)造

將易發(fā)程度等級Ⅰ,Ⅱ,Ⅲ分別記為N01,N02,N03;可以構(gòu)造出易發(fā)程度各等級的經(jīng)典物元和節(jié)域物元。

3.5.2確定待評物元

Ri為第i個待評的單元格，vij為第i個待評單元格第j個預(yù)測指標(biāo)實際值。按公式(4)計算其評價指標(biāo)C1關(guān)于各易發(fā)程度等級的關(guān)聯(lián)度。按照公式(5)計算出待評物元關(guān)于各類別等級的關(guān)聯(lián)度，評價單元(沉降中心D62監(jiān)測點(diǎn)所在單元格)計算結(jié)果見表11。

表11　關(guān)聯(lián)度計算結(jié)果

表中K01(P)=0.184；K02(P)=-0.579;K03(P)=-0.727，經(jīng)比較分析，D62所在單元格易發(fā)程度屬于Ⅰ級，即D62所在單元為地面沉降高易發(fā)區(qū)，其他單元按照同樣的方法計算確定其易發(fā)級別。根據(jù)所有單元格的評價結(jié)果，再將各剖分單元轉(zhuǎn)化成面元，最后經(jīng)過等值差分，結(jié)合地質(zhì)環(huán)境條件等綜合劃分出不同級別的地面沉降地質(zhì)災(zāi)害易發(fā)區(qū)。

4地面沉降易發(fā)區(qū)劃

根據(jù)地面沉降地質(zhì)災(zāi)害現(xiàn)狀、發(fā)育程度、危害程度、發(fā)展趨勢、地質(zhì)環(huán)境條件，結(jié)合地下水開采等評價因子，按照基于層次分析的可拓學(xué)綜合評價等級大小進(jìn)行易發(fā)區(qū)劃分，地面沉降分為高、中和低3類易發(fā)區(qū)。評價區(qū)易發(fā)程度分區(qū)結(jié)果如圖1所示。

圖1　地面沉降易發(fā)程度分區(qū)圖

4.1高易發(fā)區(qū)(Ⅰ)

位于德州市中心地帶，面積約200km2。該區(qū)深層地下水位埋深多大于100m，是深層地下水的集中開采區(qū)和地面沉降重點(diǎn)區(qū)域。地面沉降累計沉降量局部大于1000mm，1991—2010年年平均沉降速率為35～62.5mm。該區(qū)深層地下水主要為工業(yè)、第三產(chǎn)業(yè)和城市生活用水。

該區(qū)應(yīng)嚴(yán)格執(zhí)行深層地下水壓采方案，調(diào)整產(chǎn)業(yè)結(jié)構(gòu)；充分利用南水北調(diào)、黃河水源補(bǔ)充水資源；實施生態(tài)補(bǔ)償原則和用途管制，實行階梯水價；對第Ⅲ,Ⅳ含水層組實施人工回灌；對地下水、地面沉降進(jìn)行長期監(jiān)測，加強(qiáng)沉降風(fēng)險評估，建立地面沉降監(jiān)測預(yù)警信息平臺。

4.2中易發(fā)區(qū)(Ⅱ)

位于高易發(fā)區(qū)東南側(cè)，面積139km2，是工業(yè)和經(jīng)濟(jì)較為發(fā)達(dá)的地段。該區(qū)深層地下水位埋深90～110m，2010年地面沉降累計量500～1000mm，1991—2010年年平均沉降速率為20～35mm，是深層地下水開采較集中區(qū)，是地面沉降和深層地下水降落漏斗的發(fā)展地段。該區(qū)深層地下水主要為工業(yè)、農(nóng)業(yè)、第三產(chǎn)業(yè)和城鄉(xiāng)生活使用。

該區(qū)宜作為后備水源或應(yīng)急水源，實行控制性開采，盡快建立分層標(biāo)等地面沉降自動監(jiān)測體系。

4.3低易發(fā)區(qū)(Ⅲ)

位于研究區(qū)東南邊緣地段以及東部京滬高鐵沿線一帶，面積200km2，是德城區(qū)農(nóng)業(yè)和第三產(chǎn)業(yè)未來發(fā)展的重要地段。該區(qū)深層地下水位埋深50～90m，2010年地面沉降累計沉降量300～450mm，1991—2010年年平均沉降速率為10～20mm，深層地下水開采比較分散，是地面沉降和深層地下水降落漏斗的周邊地帶。該區(qū)深層地下水主要滿足鄉(xiāng)鎮(zhèn)人民的生活用水，其次是部分小企業(yè)的工業(yè)用水。

該區(qū)深層地下水實行控制開采，開發(fā)利用淺層地下水，適量引用地表水；提倡農(nóng)業(yè)節(jié)水灌溉；可在基巖埋深小的地段建立1組地面沉降監(jiān)測基巖標(biāo)。

5結(jié)論

可拓學(xué)地面沉降易發(fā)性評價具有定量化、多因子、綜合性的特點(diǎn)。它引入了描述客觀事物性質(zhì)變化的關(guān)聯(lián)函數(shù)工具，對實變函數(shù)中距離的概念加以拓展，為質(zhì)、量兼顧解決2大類7小類地面沉降易發(fā)性評價因子間的矛盾問題提供了定量手段。

該次地面沉降物元模型的建立，實現(xiàn)了在沉降易發(fā)區(qū)等級劃分時，將三元組(易發(fā)區(qū)評價等級N，關(guān)于評價因子c的量值v)放在一個統(tǒng)一體中進(jìn)行綜合評價。

綜合考慮地質(zhì)環(huán)境條件、地下水開采規(guī)劃所預(yù)測的累計沉降量等值線變化趨勢與可拓學(xué)劃分的地面沉降易發(fā)性分區(qū)形態(tài)具有一定的相似性，反映出地面沉降具有區(qū)域緩變性和受地下水開采影響明顯的特點(diǎn)。這一相似性規(guī)律需要進(jìn)一步研究。

綜上，基于層次分析法的可拓學(xué)理論進(jìn)行地面沉降易發(fā)區(qū)劃，方法是可行的，結(jié)果基本符合當(dāng)?shù)貙嶋H情況。

參考文獻(xiàn)：

[1]林良俊，孫曉明，柳富田.新時期沿海地區(qū)國土開發(fā)面臨的形勢和主要環(huán)境地質(zhì)問題及對策研究[J].水文地質(zhì)工程地質(zhì)，2010，37，(3):1-3.

[2]孟慶峰，戴魯旗.德州市地面沉降現(xiàn)狀及防治對策[J].山東國土資源，2008，24(11):8-9.

[3]張永偉.華北平原德州地面沉降成生機(jī)理、監(jiān)測預(yù)警與可控性研究[D].濟(jì)南：山東大學(xué)，2014.

[4]蔡文.物元模型及應(yīng)用[M].北京:科學(xué)技術(shù)文獻(xiàn)出版社,1994.

[5]XU Y, ZHU QX. A new extension theory-based production operation method in industrial process[J]. Chin J Chem Eng., 2013, 21(1):44-54.

[6]賈超、肖樹芳、劉寧，可拓學(xué)理論在洞室?guī)r體質(zhì)量評價中的應(yīng)用[J]. 巖石力學(xué)與工程學(xué)報，2003，(5)：751-756.

[7]ZHANG YW, LI SC & MENG FQ. Application of extenics theory for evaluating effect degree of damaged mountains based on analytic hierarchy process[J]. Environ Earth Sci, 2014,71:4463-4471. DOI 10.1007/s12665-013-2840-y.

[8]賈德旺，岳躍華，葉進(jìn)霞，等.濟(jì)寧城區(qū)地面沉降演變規(guī)律及經(jīng)濟(jì)損失估算方法[J].山東國土資源，2006，22(11):38-41.

[9]陳杰，朱國榮，顧阿明,等.Biot固結(jié)理論在地面沉降計算中的應(yīng)用[J].水文地質(zhì)工程地質(zhì)，2003，(2):28-31.

[10]龔士良.集對分析及在城市地而沉降研究中的應(yīng)用[J].水文地質(zhì)工程地質(zhì),1998，(5)：33-40.

[11]李斌，劉峰．遺傳神經(jīng)網(wǎng)絡(luò)在地面沉降趨勢預(yù)測中的應(yīng)用．山東國土資源，2013，29(8)：47-49．

[12]郭小萌，宮輝力，朱峰,等.基于矩陣的GM(1,1)模型預(yù)測地面沉降——以北京市平原區(qū)為例[J].水文地質(zhì)工程地質(zhì),2013,40(06):101-105.

[13]楊玉堂，劉清德.聊城市主要地質(zhì)災(zāi)害現(xiàn)狀及防治工作進(jìn)展[J].山東國土資源，2008，24(1-2):20-21.

[14]王小剛，鄒祖光，王秀芹,等.東營市城區(qū)地面沉降影響因素[J].山東國土資源，2006,22(5):50-53.

[15]王哲，易發(fā)成.基于層次分析法的綿陽市地質(zhì)災(zāi)害易發(fā)性評價[J].自然災(zāi)害學(xué)報，2009，18(1)：14-23.

[16]許紹雙，Excel在層次分析法中的應(yīng)用[J]中國管理信息化，2006,11(9):17-19.

收稿日期：2015-12-17；

修訂日期：2016-02-20；編輯：王敏

基金項目：山東省地面沉降監(jiān)測與防治研究(工作項目編碼：1212011220182)、華北平原德州地面沉降監(jiān)測預(yù)警系統(tǒng)建設(shè)(魯國土資發(fā)[2012]112號)

作者簡介：張永偉(1972—)，男，山東泗水人，研究員，主要從事水工環(huán)及地質(zhì)災(zāi)害防治等方面的工作；E-mail：zyw526@126.com

中圖分類號：P642.26

文獻(xiàn)標(biāo)識碼：B

Application of Extenincs Theory for Partition of Land Subsidence Susceptible Area of Dezhou City Based on Analytic Hierarchy Process

ZHANG Yongwei1,2，LI Hongxia3，SHANG Tingting1, JIA Chao2，LIANG Hao1，YANG Yabin3

(1. Shandong Monitoring Center of Geo-Environment,Jinan,250014,China；2. Shandong University, Jinan,250061, China. 3. Jinan Yuxian primary school,Jinan,250024,China ;4. Lubei Geo-engineering Exploration Institute，Dezhou, 253015,China)

Abstract:An evaluating method is proposed to the partition of susceptible area of Dezhou city in Huabei plain based on the matter element model as deduced from the theories of extenics. Using the susceptible degree, classification criteria and influencing factors of land subsidence, the elements in classical and sectorized field are defined, whereby an extension-evaluating model is developed for comprehensive evaluation of the classification of susceptible degree of land subsidence by the correlation function of matter elements and extension sets. The analytic hierarchy process methodology for calculating weight aggregation of evaluating factors is proposed. All region are divided into 135 grids using the square(2km×2km), then the split unit is transformed into surface element after equivalent difference. The evaluating results are thus obtained through such an evaluation. The calculation results show that the present method is effective and applicable to forecast the classification of susceptible degree of land subsidence on multiple-factor quantitatively, which lends the support for control planning for land subsidence.

Key words:Land subsidence; extenics; matter element model; partition of susceptible area; analytic hierarchy process