基于小波與灰色模型的滑坡時間預(yù)測預(yù)報(DOC 63頁).doc

摘 要我國是世界上滑坡災(zāi)害最為嚴重的國家之一?；聻?zāi)害往往會造成嚴重的經(jīng)濟損失，威脅人民生命安全。據(jù)統(tǒng)計，我國有70多個城市，460多個縣收到滑坡災(zāi)害的威脅及危害，每年平均至少造成15億23億元的經(jīng)濟損失，由于近年來自然災(zāi)害增多，滑坡也相應(yīng)增加。如果我們能夠知道滑坡發(fā)生的時間，積極采取應(yīng)對措施，就可以將滑坡危害降到最低，減少不必要的人員傷亡和經(jīng)濟損失。因此對滑坡時間進行預(yù)測預(yù)報研究具有十分重要的意義。本文基于以上目的，針對滑坡預(yù)測預(yù)報中存在的問題進行時間預(yù)測預(yù)報研究。本文主要通過對小波理論和灰色模型的介紹和研究，將小波降噪和灰色模型應(yīng)用于滑坡時間預(yù)測預(yù)報。主要研究內(nèi)容如下：（1）對滑坡時間預(yù)測預(yù)報及其發(fā)展趨勢做了簡單的介紹。介紹了滑坡監(jiān)測技術(shù)和滑坡預(yù)測模型。（2）時間預(yù)測預(yù)報之前需要對觀測的位移時間數(shù)據(jù)進行處理，去除噪聲，以提高預(yù)測精度。本文選用小波的方法對滑坡觀測的位移數(shù)據(jù)進行去噪處理。在此之前介紹了小波理論，小波函數(shù)和小波閾值的選取規(guī)則。后面通過選用不同的閾值來對滑坡數(shù)據(jù)進行降噪處理，并利用信噪比來比較降噪效果。（3）簡單介紹灰色系統(tǒng)模型，利用matlab編寫相應(yīng)的程序。并利用灰色模型對已經(jīng)降噪處理的數(shù)據(jù)進行時間預(yù)測，并分析結(jié)果。（4）把小波閾值降噪方法和灰色模型應(yīng)用到新灘滑坡，對觀測數(shù)據(jù)進行降噪處理并預(yù)測時間，并將結(jié)果與真實值進行對比分析。關(guān)鍵詞：滑坡 小波分析 小波閾值 灰色模型 降噪 時間預(yù)測AbstractChina is one of the countries which have the most serious landslide calamity in the world。Landslide calamity can always cause serious economic losses, which threatens the safety of peoples life.According to the statistics, about 70 cities and 460 counties in our country have been threatened and damaged by the landslide hazard. As a result, the economic usually suffers from a loss between at least 1.5 billion and 2.3 billion every year. As the natural disasters increase recently, the landslides also increase accordingly.If we know that the landslide occurred at certain time and take actions, then the landslide hazard can be reduced at most, besides we can reduce unnecessary casualties and economic losses.Consequently, research of forecasting Landslide-time is significant.Based on the purpose above, in this paper, for the problems in the prediction of landslide, we will make a detailed research. According to the introduction and study of wavelet theory and Grey model, wavelet noise reduction and Grey model can be applied to landslide prediction.The main research contents are as follows:（1）A succinct introduction of the landslide-time prediction and development tendency, including the landslide monitoring technology and landslide prediction model.（2）Before theTime prediction, In order to improve prediction accuracy, it is necessary to process observed displacement time data and remove noise. In this paper, the technique of wavelet is applied to remove the noise in the observed landslide displacement data.Wavelet theory, wavelet function and the selection rules of wavelet threshold are proposed before, Then conduct noise reduction for landslide data by different threshold values, besides, through the signal to noise ratio to compare the result of noise reduction.（3）Proposing Grey model succinctly, programming by MATLAB and using Grey model to estimate the time of the noise reductive data, then analyzing the results.（4）Apply the noise reduction method for wavelet threshold and Grey model to the new beach landslide, conduct noise reduction for the observed data and estimate the time, then compare the results to the true value.Keyword: Landslide;wavelet analysis;wavelet threshold;Grey model;wavelet denosing ;Time prediction目 錄第一章 緒論11.1 論文的研究背景11. 2滑坡災(zāi)害時間預(yù)測預(yù)報及其發(fā)展趨勢21.3 滑坡預(yù)測預(yù)報研究目的和意義41.4論文的主要內(nèi)容5第二章 滑坡監(jiān)測技術(shù)與預(yù)測模型62.1 滑坡監(jiān)測技術(shù)62.2 滑坡預(yù)報模型6第三章 小波理論83.1小波變換83.2 小波變換原理與分類93.2.1 小波變換原理93.2.2 小波變換分類93.3 多分辨率分析133.4 Mallat算法153.4.1 Mallat算法綜述153.4.2 Mallat分解算法163.4.3 Mallat合成算法173.5 常用的小波函數(shù)183.6 小波變換用于信號降噪的原理223.6.1 基本降噪模型與原理223.7 非線性小波變換閾值法253.7.1 算法實現(xiàn)253.7.2 閾值的選取263.8 小波降噪效果評價指標28第四章 滑坡監(jiān)測灰色模型304.1 灰色系統(tǒng)介紹304.2 滑坡監(jiān)測數(shù)據(jù)的生成314.2.1 累加生成(AGO)324.2.2 累減生成（IAGO）334.2.3 均值生成344.2.4 級比生成344.3 滑坡監(jiān)測灰色建模方法354.3.1 灰色模型的一般形式（GM(n,h)）354.3.2 GM(1.1)模型364.4 GM(1,1)模型的精度檢驗394.4.1 相對誤差檢驗法394.4.2 后驗差檢驗法394.4.3 關(guān)聯(lián)度檢驗法41第五章 小波降噪與灰色模型預(yù)測的實例分析425.1 新灘滑坡簡介425.1.1 滑坡復(fù)活的地質(zhì)環(huán)境425.1.2 新灘滑坡的基本特征425.2 小波數(shù)據(jù)處理435.2.1 小波閾值法降噪445.2.2 去噪結(jié)果分析495.3 灰色預(yù)測預(yù)報505.3.1 滑坡灰色模型預(yù)測505.3.2 時間預(yù)測結(jié)果分析52結(jié)論與展望53參考文獻54致謝5759長安大學碩士論文第一章 緒論1.1 論文的研究背景滑坡是斜坡上的巖土體沿某一界面發(fā)生剪切破壞向坡下運動的現(xiàn)象（DB 50/1432003，地質(zhì)災(zāi)害防治工程勘察規(guī)范）。是一種突發(fā)的而且具有嚴重危害的地質(zhì)災(zāi)害。由于滑坡發(fā)生的頻率較大、范圍比較廣，其災(zāi)害性已經(jīng)非常嚴重，是僅次于地震的地質(zhì)災(zāi)害?；露喟l(fā)生在山區(qū)，我國是一個多山的國家，山區(qū)面積占到全國陸地面積的69%，因而我國也是滑坡災(zāi)害極為頻繁、損失程度極大的國家。在過去的二十多年，我國相繼發(fā)生了一系列重大的滑坡災(zāi)害，1998年，全國共發(fā)生滑坡、崩塌和泥石流達18萬處，較大規(guī)模的有447處，造成1157人死亡，1萬多人受傷，50多萬間房屋被毀壞，經(jīng)濟損失270億元人民幣；二零零年6-7月，陜西省南部的安康等地連續(xù)受到5次特大暴雨襲擊，誘發(fā)了規(guī)模不等的滑坡、崩塌、泥石流等災(zāi)害超過4000處，造成213人死亡，并使多處橋梁、路基被毀，5條國、省干線公路和49條縣鄉(xiāng)級公路中斷，襄渝鐵路中斷行車達7天之久1。在國外，滑坡災(zāi)害一樣嚴重而且廣泛。1938年加拿大Montagneuse河谷發(fā)生m的歷史大滑坡，堵塞河流形成1.5長的水庫，至1952年被湖相沉積填滿，1988年才潰決；1952年冬天美國洛杉磯滑坡造成750萬美元的損失；1963年意大利瓦伊昂水庫發(fā)生體積約為m的大滑坡，致使當時世界上最高的雙曲拱壩失效，涌浪奪去了該壩下游2600多人的生命。1970年，秘魯因地震引起的Huascaran山區(qū)的一起山崩，造成18000人死亡。1977年瑞典維塔滑坡， 150多間房屋收到破壞或重創(chuàng)，9人死亡，直接經(jīng)濟損失約為1 5000萬瑞典克朗；1999年10月，因大雨墨西哥普遍引發(fā)滑坡，造成死亡、失蹤約600人，20多萬人無家可歸。事實上，許多國家、地區(qū)，如日本的新潟縣及其臨縣，英國的南威爾士，俄羅斯的高加索與黑海沿岸，美國加州、新澤西、德克薩斯，肯尼亞的中部地區(qū)，法國南部阿爾卑斯，意大利中部等，均為滑坡多發(fā)地區(qū)或曾經(jīng)報道過大型滑坡災(zāi)害32 34。通過上面的敘述可以發(fā)現(xiàn)，滑坡災(zāi)害危害非常大，如果我們可以預(yù)先知道其發(fā)生的地點、時間、強度和影響，就可以預(yù)先進行防范，減少損失，因此對滑坡災(zāi)害預(yù)測預(yù)報進行深入、系統(tǒng)的研究是十分必要的，該研究具有重要的學術(shù)價值和實用意義。20世紀60年代滑坡預(yù)測預(yù)報研究是才開始發(fā)展起來，現(xiàn)在已經(jīng)成為滑坡研究中的一個熱門課題。因為滑坡問題比較復(fù)雜，現(xiàn)在滑坡時間預(yù)測預(yù)報還是一個世界性的科學難題。在滑坡研究領(lǐng)域，國內(nèi)外專家學者已取得了一些成績，也有預(yù)報成功的例子，如長江三峽的新灘滑坡等。但這些成功預(yù)報大多是通過監(jiān)測工作實現(xiàn)的。但是并沒有完全解決該用怎樣的理論，建立何種理論模型對滑坡進行預(yù)測的問題。甚至我們不知道是否有這個模型，是否能夠?qū)Ρ浪?、滑坡進行有效的時間預(yù)報，還是一個處于探索研究中的問題33。從日本學者齋藤迪孝(M. Satio)提出的滑坡時間預(yù)報經(jīng)驗公式至今四十多年的研究歷史，可以把滑坡預(yù)測預(yù)報研究歷史劃分為四個階段2：1、現(xiàn)象、經(jīng)驗預(yù)測預(yù)報階段(20世紀60-70年代)60年代日本學者齋藤迪孝(M. Satio)提出一個預(yù)報滑坡的經(jīng)驗公式及圖解，即著名的“齋藤法”。根據(jù)1969年智利Chuquicamata礦滑坡監(jiān)測時間-位移曲線Hoek(1977)，提出了與齋藤類似的外延線預(yù)測預(yù)報法。F.O.Jones（1961）、P.C.Stevenson（1977）、T.H.Nilsen,et al.(1970)、E.Fussganger(1976)、G.Guidicini（1976）等先后對滑坡預(yù)報的統(tǒng)計學或經(jīng)驗法方法進行研究。2、統(tǒng)計、灰色預(yù)測預(yù)報階段(20世紀80年代) 應(yīng)用統(tǒng)計數(shù)學理論建立滑坡變形參數(shù)和時間之間相互關(guān)系來進行滑坡預(yù)測預(yù)報，其主要有回歸預(yù)測模型、時間序列分析預(yù)測模型及灰色系統(tǒng)模型等三大類。3、非線性理論預(yù)測預(yù)報階段(20世紀90年代)這一階段的特點有多種預(yù)報方法的綜合研究與應(yīng)用；廣泛的現(xiàn)代數(shù)理科學新理論應(yīng)用于滑坡預(yù)報理論研究；滑坡預(yù)報的技術(shù)手段得到前所未有的發(fā)展。4、系統(tǒng)綜合和實時跟蹤動態(tài)預(yù)測預(yù)報階段(20世紀90年代末期-至今)。1. 2滑坡災(zāi)害時間預(yù)測預(yù)報及其發(fā)展趨勢滑坡的時間預(yù)測預(yù)報分為三個階段，即長期，短期，臨滑三個階段32：1、長期預(yù)報 滑坡災(zāi)害長期預(yù)報是以某地區(qū)或某一滑坡對象在今后可能活躍的年份范圍為預(yù)測目標，其結(jié)果是近似的。長期預(yù)測需要的資料有：歷史滑坡及其活動性材料、易滑坡地區(qū)的地質(zhì)情況。例如:一些地方的區(qū)域性滑坡災(zāi)害大部分情況與特殊的氣象情況有聯(lián)系。大區(qū)域的滑坡活動很多時候是由長時間下雨與特大暴雨引發(fā)的。因此，通過研究滑坡災(zāi)害活動與氣候情況之間的關(guān)系，了解滑坡旋回周期和氣候變化周期之間的對應(yīng)關(guān)系對滑坡災(zāi)害長期時間預(yù)測是有很大的作用。另外，地震也很容易產(chǎn)生新的滑坡并會使老滑坡復(fù)活，通過地震活動與區(qū)域構(gòu)造應(yīng)力場的周期性與規(guī)律性的特點，可以大概預(yù)測沿活動構(gòu)造帶的滑坡活動性的特點，然而這是趨勢性的、超長周期的預(yù)測。自然界許多的大型滑坡，特別是一些千萬方的特大型滑坡，不單單和外界條件有關(guān)，很大程度上是受地質(zhì)條件與內(nèi)動力因素影響，這種滑坡的發(fā)生、發(fā)展和演化具有自身演化的旋回周期性。因此研究這些滑坡災(zāi)害活動的歷史對研究滑坡長期預(yù)測預(yù)報有很大的意義。2、短期預(yù)報 滑坡短期時間預(yù)報是指依據(jù)滑坡影響因素對滑坡作用的強度與隨時間的變化規(guī)律，在掌握滑坡動態(tài)發(fā)展趨勢的基礎(chǔ)上，分析推斷滑坡短時間內(nèi)的演化情況，預(yù)測可能發(fā)生滑坡的時間，一般是預(yù)測的月份。該預(yù)測的兩個重要途徑為:滑坡自身運動的內(nèi)在規(guī)律和活動性趨勢分析?；禄顒有杂|發(fā)因素的規(guī)律性與滑坡活動性之間關(guān)系的研究途徑。3、臨滑預(yù)測預(yù)報 滑坡災(zāi)害臨滑預(yù)測與預(yù)報是不盡相同的，滑坡時間預(yù)測不僅是一個科學問題，而且是一個實際問題，臨滑預(yù)報還有管理上的特征。但是滑坡時間臨滑預(yù)測與預(yù)報又是一體的?；骂A(yù)測正確的例子顯示，只有持續(xù)的監(jiān)測滑坡才能夠避免滑坡的發(fā)生。但是實際上，我們不可能不間斷的去監(jiān)測所有對我們有危害的滑坡，因為我們沒有足夠資本。因此，從這個方面來說，我們無法避免滑坡造成的損失，但是如果我們能夠在滑坡發(fā)生前能夠通過滑坡的一些信息預(yù)測出滑坡發(fā)生的具體時間，就可以減輕滑坡造成的危害?；卵芯砍跗?，日本學者日本學者齋藤迪孝(M. Satio)研究認為土質(zhì)類滑坡運動分為三個階段：初始蠕變階段，等速蠕變階段，加速位移階段。我們可以利用第三階段的滑坡位移監(jiān)測數(shù)據(jù)根據(jù)曲線圖解法對滑坡進行臨滑預(yù)報。但是，我們并不知道是不是所有的滑坡都適合這個規(guī)律。因此，滑坡時間臨滑預(yù)測預(yù)報需要不間斷的位移監(jiān)測數(shù)據(jù)，而且觀測數(shù)據(jù)必須有好的精度，數(shù)據(jù)也應(yīng)該足夠多。只有這樣，才能更準確的預(yù)測時間。另外，開展位移場以外的物理監(jiān)測，如摩擦熱、聲發(fā)射等監(jiān)測，也是擴大臨滑預(yù)測預(yù)報的重要途徑。目前滑坡預(yù)測預(yù)報已經(jīng)取得了巨大的成就。但是因為一些不可控制的，隨即的因素影響，滑坡時間預(yù)測預(yù)報現(xiàn)在仍然是一個世界性的難題，很難準確的預(yù)報滑坡發(fā)生的時間。現(xiàn)在絕大部分的滑坡都不可能做到發(fā)生前進行預(yù)報，主要是因為以下兩方面的原因：一是由于造成滑坡災(zāi)害的因素非常復(fù)雜，有眾多的不確定因素，滑坡的本身破壞方式，變形機制與變形過程多種多樣。主要造成我們很難準確預(yù)報滑坡。另外，由于受到各種條件的限制，絕大部分的滑坡都沒有進行長時間的監(jiān)測，也沒有采取防治措施，因此造成很大部分滑坡突然發(fā)生。還有部分人覺得滑坡時間預(yù)測非常難，而滑坡空間預(yù)測比較簡單，因此導致滑坡時間預(yù)測的研究被忽視。以上都是現(xiàn)階段我們國家滑坡預(yù)測中存在的問題。近三十多年來，國內(nèi)外專家學者對滑坡的研究做了大量的工作，獲得了一定的成果，但是現(xiàn)有的研究成果還是難以實現(xiàn)滑坡災(zāi)害的滑動時間，活動強度，危害程度等方面進行全面，快速，準確的預(yù)報，因此防災(zāi)決策部門對滑坡預(yù)報有進一步更高的要求。目前對于滑坡災(zāi)害時間預(yù)測預(yù)報還處在初級階段，需要有更大的發(fā)展。普遍適用性和外推適應(yīng)性程度低是現(xiàn)在已經(jīng)有的滑坡時間預(yù)測預(yù)報理論和方法存在的主要問題。因而在以后的滑坡時間預(yù)測預(yù)報的研究中將會向以下幾個方面發(fā)展：(1)滑坡動態(tài)變形的系統(tǒng)監(jiān)測和綜合參數(shù)預(yù)測預(yù)報模型研究。(2)基于非線性科學理論的滑坡時間預(yù)測預(yù)報技術(shù)研究。 (3)人類工程活動與滑坡活動之間的相互作用關(guān)系研究。 (4)滑坡動態(tài)變形過程中地表水和地下水的作用機理研究。(5)滑坡變形破壞作用機理研究。 (6)滑坡時間預(yù)測預(yù)報的組合預(yù)測模型研究。 (7)滑坡變形破壞的計算機仿真模擬預(yù)測研究。(8)新理論新方法在滑坡預(yù)測預(yù)報工作中的應(yīng)用研究。1.3 滑坡預(yù)測預(yù)報研究目的和意義滑坡災(zāi)害與土地資源利用、能源開發(fā)、建筑物穩(wěn)定性、環(huán)境保護、生命與財產(chǎn)保障等具有密切的關(guān)系?；骂A(yù)測預(yù)報，不管是空間預(yù)測預(yù)報、時間預(yù)測預(yù)報，還是災(zāi)情損失預(yù)測預(yù)報，都是工程地質(zhì)和環(huán)境地質(zhì)研究領(lǐng)域內(nèi)的世界難題，開展滑坡災(zāi)害預(yù)測預(yù)報的目的可概括為以下幾點：(1)為重大國民經(jīng)濟建設(shè)的規(guī)劃提供宏觀災(zāi)害的預(yù)防，減災(zāi)和避災(zāi)的科學依據(jù)。 (2)為具體建筑物和建設(shè)場地的布置，選擇無滑坡且不易形成新的滑坡的地段，使得工程能夠正常進行。 (3)預(yù)測滑坡防治的必要性和可能的措施，預(yù)測不加固斜坡而使用建筑場地的可能性。 (4)及時預(yù)報滑坡災(zāi)害活動和工程或建筑物的危險位移量，以此來預(yù)測滑坡發(fā)生的時間，避免滑坡產(chǎn)生大位移所造成的危害。減少或防止人員傷亡和經(jīng)濟損失。 (5)為滑坡多發(fā)危害性大的地區(qū)制定相關(guān)政策法規(guī)提供科學的理論依據(jù)。確定相關(guān)的防治措施。1.4論文的主要內(nèi)容近年來，小波理論在地球科學領(lǐng)路得到了廣泛的應(yīng)用，滑坡過程是一個復(fù)雜的過程，我們?nèi)绻軌蚶眯〔ɡ碚搧韺骂A(yù)測預(yù)報進行研究，將有利于我們了解滑坡過程的本質(zhì)。本文將采用小波理論和灰色模型結(jié)合的方式來對滑坡數(shù)據(jù)進行處理預(yù)測預(yù)報。我們知道小波在信號處理方面已經(jīng)有很好的成績，特別是在信號去噪方面，所以我們將利用小波分析的信號降噪的方法來對滑坡體位移監(jiān)測曲線進行降噪處理，在此基礎(chǔ)上再利用灰色預(yù)測模型來預(yù)測滑坡發(fā)生的時間。本文的主要內(nèi)容有下面幾個方面：1、滑坡監(jiān)測技術(shù)與預(yù)報的模型介紹。2、小波理論與灰色理論介紹。3、利用小波理論在滑坡時間預(yù)測預(yù)報前對滑坡體位移時間曲線的降噪。小波降噪的閾值選取規(guī)則，通過選取不同閾值來對滑坡時間位移曲線進行處理并通過信噪比來比較結(jié)果的優(yōu)劣。4、建立灰色模型，并把小波降噪后的數(shù)據(jù)代入到灰色模型中，以新灘滑坡為例進行時間預(yù)測預(yù)報輸出預(yù)測預(yù)報結(jié)果。第二章 滑坡監(jiān)測技術(shù)與預(yù)測模型2.1 滑坡監(jiān)測技術(shù) 從不同監(jiān)測方法獲得的觀測量的屬性來看，監(jiān)測內(nèi)容可以分為四大類：變形監(jiān)測分為地表變形監(jiān)測和深部位移監(jiān)測。地表變形監(jiān)測以常規(guī)測量為主要監(jiān)測技術(shù)（如大地測量、GPS測量等），滑坡體位移監(jiān)測是滑坡體最直觀的信息。地下水監(jiān)測主要監(jiān)測地下水位（或地下水壓力）、孔隙水壓力和地下水水質(zhì)等。應(yīng)力應(yīng)變監(jiān)測主要監(jiān)測地應(yīng)力以及地應(yīng)力變化引起的巖土聲發(fā)射現(xiàn)象、邊坡錨固后的錨桿和錨索的應(yīng)力應(yīng)變等。誘發(fā)因素監(jiān)測主要監(jiān)測氣象、地震、人類工程活動等30?；麦w地表變形監(jiān)測一般采用常規(guī)大地測量方法、全球定位系統(tǒng)（GPS）方法、測量機器人技術(shù)、三維激光掃描技術(shù)、地面數(shù)字攝影測量方法、遙感（RS）方法等。滑坡體深部位移監(jiān)測主要通過對埋設(shè)在地質(zhì)鉆孔中的測斜管、多點位移計、滑動變形計以及埋設(shè)在豎井和平硐中的伸縮儀等儀器設(shè)備的變位情況來獲取資料。另外還有滑坡地下水觀測，巖石邊坡聲發(fā)射監(jiān)測，邊坡錨固力監(jiān)測，地溫和地應(yīng)力監(jiān)測和誘發(fā)因素大氣降雨、氣溫、相關(guān)人類活動監(jiān)測和宏觀地質(zhì)巡查等。滑坡監(jiān)測剖面可以分為重要監(jiān)測剖面和一般監(jiān)測剖面?；卤O(jiān)測剖面的選擇應(yīng)該結(jié)合地質(zhì)結(jié)構(gòu)、成因機制、變形特征，監(jiān)測點則布置在滑動面的兩側(cè)或者滑動面的上下盤面上，并分析該剖面上全部監(jiān)測點的功能并予以綜合，建立該剖面在水平方向和豎直方向上代表滑坡體的變形塊體范圍及其組合。滑坡變形監(jiān)測網(wǎng)和監(jiān)測點應(yīng)分為滿足不同的布設(shè)及要求。其他專用監(jiān)測儀器，如伸縮儀、測縫計、水平位錯計、液體靜力水準儀等，根據(jù)勘探平硐、豎井所在位置，適量布設(shè)在滑動面兩側(cè)或上下。邊坡錨固監(jiān)測應(yīng)依據(jù)規(guī)范和設(shè)計要求按一定比例安裝監(jiān)測儀器進行錨固力變化觀測?；卤O(jiān)測的監(jiān)測網(wǎng)點、水準點、基巖標埋設(shè)以及多點位移計埋設(shè)安裝等，必須滿足相關(guān)規(guī)范和設(shè)計的要求。用于滑坡監(jiān)測的各類儀器設(shè)備必須經(jīng)過檢定合格后才可以使用。針對不同的監(jiān)測儀器設(shè)備，原始數(shù)據(jù)采集可以分別設(shè)計為人工采集、半自動采集、自動采集三種方式來完成。2.2 滑坡預(yù)報模型滑坡預(yù)測預(yù)報研究的核心和關(guān)鍵問題就是建立一種與滑坡體相適應(yīng)，能夠準確預(yù)報，操作性強的預(yù)測預(yù)報模型。國內(nèi)外許多的專家學者經(jīng)過一段時間的研究，從不同的角度提出了各種各樣的預(yù)報模型。綜合目前國內(nèi)外已有的預(yù)報方法可以分為以下三類。1.確定性預(yù)報模型確定性預(yù)報模型是把與滑坡以及環(huán)境有關(guān)的各類參數(shù)用測定的量予以數(shù)值化，并利用明確的函數(shù)來表達其數(shù)學關(guān)系。該模型的預(yù)測反映了滑坡的物理性質(zhì)，更適合于斜坡或者山體滑坡單體的預(yù)測。確定性方法有齋藤迪孝方法、有限元，邊界元，離散元，耦合方法，還有極限平衡分析法。最先發(fā)展起來的齋藤迪孝方法和與它相關(guān)的一些方法預(yù)測的精度不高，計算的時間也屬于大概的時間，并不準確，這些方法適合滑坡中短期、臨滑預(yù)報。有限元，邊界元，離散元的運用可以模擬分析滑坡穩(wěn)定性以及滑坡位移變化，因此來對滑坡進行預(yù)報。考慮滑坡各種影響因素的極限分析法具有物理概念清晰，計算也很簡單的特點，這種方法適合滑坡長期預(yù)報。2.非確定性預(yù)報模型非確定性預(yù)報模型主要的重點是對已經(jīng)存在的滑坡和與其相關(guān)的外界因素，地質(zhì)環(huán)境條件進行調(diào)查，掌握這些統(tǒng)計規(guī)律，并不注重滑坡機理嚴格的數(shù)學表達。非確定性模型比較適合區(qū)域性的土地的利用和國土開發(fā)規(guī)劃。建立在因果分析和統(tǒng)計分析基礎(chǔ)上的各種預(yù)報模型都可以歸屬此類。例如移動平均法、灰色模型、指數(shù)平滑、非線性預(yù)測模型、時間序列分析預(yù)測模型、回歸分析法、趨勢疊加法、生長曲線法、卡爾曼濾波法、動態(tài)跟蹤法等多種方法。這些模型和滑坡觀測的數(shù)據(jù)有關(guān)系，只要滑坡觀測的位移數(shù)據(jù)足夠多，分布規(guī)律就能夠很好的進行滑坡預(yù)測預(yù)報。非確定性模型適用于滑坡中短期預(yù)報。3.工程地質(zhì)類比模型 工程地質(zhì)類比模型是比較被測對象近似的參照對象，類比推理被測對象以后的發(fā)展情況，是介于確定性模型與非確定性模型之間的一種方法。一些直觀的評判方法、黃金分割法、綜合信息預(yù)報模型以及力學圖解法都屬于此種模型。黃金分割法簡單適用，精度也高，多適用于滑坡的中長期預(yù)報。綜合信息預(yù)報模型適用于臨滑預(yù)報，可以識別滑坡所處的變形階段，便于現(xiàn)場操作。力學圖解法可用做滑坡的判據(jù)，并能判定滑坡的破壞型式。模糊評判法、神經(jīng)網(wǎng)絡(luò)法也可以歸入此類。第三章 小波理論3.1小波變換為了克服短時傅里葉變換的單一分辨率缺陷，發(fā)展了小波變換。小波變換是一種工具，它把數(shù)據(jù)、函數(shù)或算子分割成不同頻率的成分，然后再用分解的方法去研究對應(yīng)尺度下的成分。小波變換繼承并發(fā)展了短時傅里葉變換的局部化思想，并且克服了其窗口大小與形狀固定不變的缺點。它不但可以同時從時頻域觀測信號的局部特征，而且時間分辨率和頻率分辨率都是可以變化的，是非常好的信號處理方法。由于小波變換分析時、頻局域化的卓越效果而譽滿全球，因此它成為分析非平穩(wěn)信號的重要工具。因此，小波變換成為信號處理、信息獲取與處理等多領(lǐng)域首選的數(shù)學分析工具，在國內(nèi)外形成一次又一次的研究高潮。它與分形理論、神經(jīng)網(wǎng)絡(luò)聯(lián)合互補，對當今自然科學產(chǎn)生巨大的沖擊力和推動作用67。1981年，法國物理學家Morlet在分析地質(zhì)數(shù)據(jù)時基于群論首先提出了小波變換這一概念。1985年，法國數(shù)學家Meyer首先提出光滑的小波正交基，后被稱為Meyer基，對小波理論做出了巨大貢獻。1988年，數(shù)學家Daubechies提出了具有緊支集的光滑正交小波基Daubechies基，為小波應(yīng)用研究增添了催化劑。后來，信號分析專家Mallat提出了多分辨分析的概念，給出了構(gòu)造正交小波基的一般方法。Mallat受金字塔算法的啟發(fā)，以多分辨分析為基礎(chǔ)提出了著名的快速小波算法Mallat算法，這是小波理論突破性的成果，Mallat算法的提出，宣告了小波從理論研究走向?qū)拸V的應(yīng)用研究8。小波變換是傅里葉變換發(fā)展史上的里程碑式的進展910111213。它是一種信號的時間尺度（時間頻率）分析方法，它具有多分辨分析（Multiresolution Analysis）的特點，而且在時頻兩域都具有表征信號局部特征的能力，是一種窗口大?。创翱诿娣e）固定不變，而其形狀是可以改變的。即在低頻部分具有較高的頻率分辨率與較低的時間分辨率，而在高頻部分具有較高的時間分辨率與較低的頻率分辨率，被譽為“數(shù)學顯微鏡”，它能根據(jù)觀察的對象自動“調(diào)焦”，對信號進行多尺度細化分析，以取得令人滿意的最佳效果；正是這種特性，使小波變換具有對信號的“自適應(yīng)性”，它能根據(jù)被分析的對象自動調(diào)整有關(guān)參數(shù)，很適合于探測正常信號中夾帶的瞬態(tài)反?，F(xiàn)象并展示其成分，可見利用連續(xù)小波變換進行動態(tài)系統(tǒng)故障檢測與診斷具有良好的效果。3.2 小波變換原理與分類3.2.1 小波變換原理小波變換從誕生開始，已經(jīng)有了幾種改型，其中常見的三種形式如下。（1） 經(jīng)典的小波變換。 （3.1）式中，基本小波滿足允許條件和非零條件。（2） 改正的小波變換。 （3.2）（3） 進一步改正的小波變換。 （3.3）式中，即為經(jīng)典小波變換式（3.1）；即為改正的小波變換式（3.2）。3.2.2 小波變換分類小波變換有很多種形式，所有形式都來源于基本公式。有連續(xù)小波變換和離散小波變換兩種。在離散小波變換中我們進一步做如下定義：有冗余的離散系統(tǒng)（框架）和正交小波基及其他小波基。3.2.2.1 連續(xù)小波變換一維信號的連續(xù)小波變換，就是利用一個小波函數(shù)將信號變換到時間尺度域上。它的時頻局部化是變化的，在高頻處時間分辨率高，低頻處時間分辨率低。其定義如下：設(shè)（其中代表平方可積的實數(shù)空間，即能量有限的信號空間），它的傅里葉變換是，則當滿足下面的條件： （3.4）時，則是一母小波或基本小波。把經(jīng)伸縮、平移就得到一個小波序列： （3.5）式中，代表伸縮因子；代表平移因子。對于任意函數(shù)，其連續(xù)小波變換定義為： （3.6）小波變換的逆變換（重構(gòu)公式）為： （3.7）因為在小波變換中起著觀測窗的作用，所以還應(yīng)該滿足一般函數(shù)的約束條件： （3.8）故是R中的一個連續(xù)函數(shù)。為了使重構(gòu)公式成立，在原點必定為零，即： （3.9）由式（3.9）可以發(fā)現(xiàn)小波函數(shù)必然具有振蕩性。為了使信號重構(gòu)的實現(xiàn)在數(shù)值上是穩(wěn)定的，不僅要完全重構(gòu)條件外，還需要小波的傅里葉變換滿足以下的穩(wěn)定性條件： （3.10）從穩(wěn)定性條件出發(fā)，可以定義一個對偶小波： （3.11）連續(xù)小波變換具有如下性質(zhì)：（1）線性性：各個分量的小波變換之和與一個多分量信號的小波變換相等。設(shè)，則 （3.12）（2）平移不變性：若的小波變換為，則的小波變換為 （3.13）（3）伸縮共變性：若的小波變換為，則的小波變換為 （3.14）式中，。（4）自相似性：相對于不同的平移參數(shù)和不同的尺度參數(shù)的連續(xù)小波變換之間是自相似的。（5）冗余性：連續(xù)小波變換中存在信息表述的冗余度（redundancy）。盡管冗余可以提高信號重建時計算的穩(wěn)定性，但是它增加了小波變換的計算時間和分析、解釋小波變換結(jié)果的困難。小波變換的冗余性事實上也是自相似性的直接反映，它主要表現(xiàn)在以下兩個方面：（1）由連續(xù)小波變換恢復(fù)原信號的重構(gòu)分式不是唯一的。也就是說，信號的小波變換與小波重構(gòu)不存在一一對應(yīng)關(guān)系。（2）小波變換的核函數(shù)即小波函數(shù)存在許多可能的選擇，它們可以是非正交小波、正交小波、雙正交小波，甚至允許是彼此線性相關(guān)的。3.2.2.2 離散小波變換如果在每個可能的尺度下都去計算小波系數(shù)，計算量將相當大，會產(chǎn)生驚人的數(shù)據(jù)量。在實際應(yīng)用過程中，特別是在計算機上實現(xiàn)的時候，連續(xù)小波必須加以離散化。這一離散化只是針對連續(xù)的尺度參數(shù)和連續(xù)平移參數(shù)的，而不是針對時間變量的。連續(xù)小波變換中小波函數(shù)的定義為式（3.5），在離散化時，通?？傁拗浦蝗≌?，令，可見在離散小波變換中，相容性條件為 （3.15）一般最常用的辦法就是將尺度參數(shù)和平移參數(shù)分別令，。對應(yīng)的離散小波函數(shù)： （3.16）對應(yīng)的離散小波系數(shù)為 （3.17）其重構(gòu)公式為 （3.18）式中，是一個與信號無關(guān)的常數(shù)。為了使小波變換具有可變化的時間分辨率與頻率分辨率，適應(yīng)待分析信號的非平穩(wěn)性，通常需要改變尺度參數(shù)和平移參數(shù)的大小，使得小波變換具有“變焦距”的功能。在實際應(yīng)用中主要采用二進制動態(tài)采樣網(wǎng)絡(luò)。經(jīng)這種離散化后的小波和相應(yīng)小波變換，稱之為二進小波（Dyadic Wavelet）和二進小波變換。此時，取，。式（3.16）變?yōu)?（3.19）二進小波對信號的分析具有變焦距的功能。如果有一個放大倍數(shù)為，它對應(yīng)的為觀測到信號的某部分內(nèi)容。二進小波不同于連續(xù)小波的離散化，它只是對尺度參數(shù)進行離散化，而對時間域上的平移參數(shù)保持連續(xù)變化。二進小波不破壞信號在時間域上的平移不變量，這也正是二進小波和正交小波基相比所具有的獨特的優(yōu)點。設(shè)函數(shù)，假設(shè)存在兩個常數(shù)、，且使得穩(wěn)定性條件幾乎全部滿足，即 （3.20）則為一個二進小波。若，則稱為最穩(wěn)定條件。而函數(shù)序列叫做的二進小波變換，其中 （3.21） 式（3-21）對應(yīng)的逆變換為 （3.22）3.3 多分辨率分析多分辨率分析（Multi-Resolution Analysis，簡稱MRA）的概念是1988年S.Mallat在構(gòu)造正交小波基時，將計算機視覺領(lǐng)域內(nèi)的多分辨率思想率先引入小波變換，在空間的概念上形象地說明了小波變換的多分辨特性，并使用多分辨分析將此前Meyer等提出的各種具體小波基的構(gòu)造法統(tǒng)一起來，在Burt and Adelson圖像分解和重構(gòu)的塔式算法的啟發(fā)下，基于它的多分辨框架，提出了塔式多分辨率分解與重構(gòu)算法7，給出了正交小波的構(gòu)造方法以及正交小波的快速算法Mallat快速小波分解和重構(gòu)算法81415。對于多分辨率分析的理解，可以用一個三層的分解進行說明，其小波分解樹如圖3.1所示。SA1D1A2D2D3A3 圖3.1 三層多分辨分析樹結(jié)構(gòu)圖 上圖中，A表示低頻，D表示高頻，末尾的序號數(shù)表示小波分解的層數(shù)（也即尺度數(shù)）。從圖上可以看出，多分辨分析僅僅是對低頻的部分進行進一步分解，而不予以考慮高頻的部分。信號S的分解關(guān)系式為：。需要強調(diào)的是，這里僅僅是用一個層分解來進行說明的，如果要進行進一步的分解，則可以把低頻部分A3分解成低頻部分A4和高頻部分D4，如繼續(xù)再分解，則以此類推。多分辨分析的最終目的是為了構(gòu)造一個在頻率上高度逼近L2(R)空間的正交小波基（或正交小波包基），這些頻率分辨率不同的正交小波基相當于帶寬各異的帶通濾波器。從三層多分辨分析樹型結(jié)構(gòu)圖3.1可以看出，多分辨分析只對低頻空間部分進行進一步的分解，使頻率的分辨率變得越來越高。 空間L2(R)中的多分辨率分析是指L2(R)中滿足如下條件的一個空間序列： (3.23)且有 (3.24) (3.25) (3.26)即空間序列中的每一個空間均為空間序列中心V0的尺度變換。多分辨分析的另一重要特性：空間V0在整數(shù)平移下的不變性，即 (3.27)上式蘊涵：若，則，有。另外，我們常常要求存在函數(shù)，使 是V0上的一個規(guī)范正交基 (3.28)其中，規(guī)定對于，。由此，式（3.26）和式（3.28）蘊涵：，是Vj上的一個規(guī)范正交基。我們通常稱為多分辨分析的“尺度函數(shù)”。多分辨分析的基本思想是：當一組閉子空間滿足式（3.23）式（3.28）時，存在一個一個上的規(guī)范正交基，其中。對于，有 （3.29）（為上的正交投影）。而且小波有顯式的構(gòu)造公式。 對于，定義為在上的正交補空間，于是有 （3.30）及 （3.31）（事實上，假設(shè)，則）。于是，當時，有 （3.32）其中，全部所有的子空間都是正交的。根據(jù)式（3.24）和式（3.25），有 （3.33）上式將分解為一系列相互正交的子空間，而且空間繼承了尺度特性，即 (3.34)式（3.29）可等價地表達為：對于固定的，構(gòu)成了上的一組規(guī)范正交基。再根據(jù)式（3.33）和式（3.34）可知：整個小波族是上的一個規(guī)范正交基。3.4 Mallat算法81415Mallat算法在小波變換里占據(jù)非常重要的地位，相當于傅里葉變換中的快速傅里葉算法。Mallat在1987年提出了多分辨分析的概念，通過將計算機的多尺度分析思想引入到小波分析中而提出的，統(tǒng)一了在此之前的所有具體正交小波基的構(gòu)造，并且提出了相應(yīng)的分解和重構(gòu)的快速算法，即Mallat算法。3.4.1 Mallat算法綜述多分辨分析是構(gòu)造正交小波基的有力工具，下面從多分辨分析出發(fā)，我們可以得到正交小波變換和反變換的快速算法。設(shè)是中的正交小波基，則對于任意信號，有下面的無窮級數(shù)展開式： （3.35）由于一般的不具有初等解析表達式，直接用公式求展開系數(shù)非常不方便。取任意的，設(shè)為空間到空間的正交投影算子，為空間到空間的正交投影算子，則有： （3.36）其中，小波系數(shù)實際上就是離散二進網(wǎng)格上的小波變換，這樣的多分辨分析和小波分析就聯(lián)系起來了。由得到：，即下面的式子成立： （3.37）3.4.2 Mallat分解算法Mallat分解算法是計算尺度系數(shù)和小波系數(shù)的算法，為了簡化分析，假設(shè)、都是實函數(shù)，在數(shù)學推導過程中忽略積分式中的復(fù)共軛。由二尺度方程得到： （3.38）從而有： （3.39）同理可以推出： （3.40）上面式（3.39）和（3.40）就是Mallat算法的分解公式。當時，充分逼近，因此，任取，可以取充分小的，使得。由于是一個雙向無限的空間帶，在數(shù)學上為方便起見，我們可以認為，則有下面的分解關(guān)系?？臻g上的分解關(guān)系為： （3.41）信號上的分解關(guān)系為（）： （3.42）Mallat算法的分解過程可以用一個金字塔逐次分解的圖來表示，令，系數(shù)的塔式分解圖如下圖3.2所示：C0DnD3D1C1CnCn-1C2D2圖3.2 Mallat算法塔式分解圖3.4.3 Mallat合成算法Mallat合成算法就是根據(jù)圖3.2中的系數(shù)Cn、Dn、D2 、D1推導出C0。由上面的推導得到： （3.43）將上式兩邊同時與作內(nèi)積得： （3.44）即為： （3.45）還因為： （3.46） （3.47）將式（3.46）和式（3.47）帶入（3.45）于是得到Mallat算法的合成公式： （3.48）3.5 常用的小波函數(shù)1、Haar小波函數(shù)81516Haar小波函數(shù)是在小波分析中最早用到的一個具有緊支撐的正交小波函數(shù)，A.Haar于1990年提出的一種正交函數(shù)系。Haar小波的構(gòu)造過程為：給出尺度函數(shù)： （3.49）計算： （3.50）構(gòu)造小波： （3.51）Haar小波函數(shù)具有正交性、雙正交性、緊支撐性（支撐長度為1）、對稱性。它既可以用于連續(xù)小波變換，也可以用于離散小波變換。其消失矩階數(shù)為1，濾波器長度為2。圖3.3 Haar小波的尺度函數(shù)和小波函數(shù)圖2、Meyer小波 Meyer小波是在頻率域中定義的。尺度函數(shù)： （3.52）小波函數(shù)； （3.53）Meyer小波不是緊支撐的，但是收斂得很快。圖3.4 Meyer小波的尺度函數(shù)和小波函數(shù)圖3、Daubechies(dbN)小波81516Daubechie函數(shù)是世界著名小波分析學者Inrid Daubechie對尺度為2的整數(shù)冪（即）構(gòu)造的小波函數(shù)。Haar小波是其尺度為1的特例，其他的小波沒有明確的表達式。此小波函數(shù)具有正交性、雙正交性、緊支撐性，既可以進行連續(xù)小波變換也可以進行離散小波變換，其支撐寬度為2N-1，濾波器長度為2N，小波函數(shù)的消失矩階數(shù)為N。大多數(shù)dbN不具有對稱性，該小波還具有正則性。圖3-5 db6小波的尺度函數(shù)和小波函數(shù)圖4、Symlets小波它是對dbN小波函數(shù)的一種改進，性質(zhì)與dbN小波相似，但更近似于對稱小波。5、Morlet小波該小波常用的是復(fù)值Morlet小波： （3.54）其傅里葉變換為： （3.55）由式（3.55）可以看出，當時，滿足容許條件，因此Morlet是容許小波。Morlet小波是一種復(fù)值小波，能夠提取信號中的幅值和相位信息，其時域、頻域都具有很好的局域性，常用于復(fù)數(shù)信號的分解及時頻分析中，Morlet小波在推廣到N維時具有很好的角度選擇性。另外，Morlet小波沒有尺度函數(shù)，而且是正交分解的。6、墨西哥帽(Mexican Hat)小波墨西哥帽(Mexican Hat)小波函數(shù)是Gauss函數(shù)的二階導數(shù)，由于其小波函數(shù)的圖形像墨西哥帽的截面，所以稱這個函數(shù)為墨西哥帽函數(shù)。墨西哥帽小波函數(shù)為 （3.56）其小波函數(shù)在時域與頻域都具有很好的局部性，具有對稱性，它適用于連續(xù)的小波變換，它的尺度函數(shù)不存在，因此不具有正交性。圖3-6 Morlet和Mexican Hat小波的小波函數(shù)圖7、雙正交小波（biorthogonal小波）雙正交小波（biorthogonal小波）的引入是為了解決對稱性與精確信號重構(gòu)