石羊河流域地下水

1、石羊河流域地下水時(shí)空動(dòng)態(tài)研究答 辯 人：朱小燕 指導(dǎo)老師：馬金輝副教授地理信息系統(tǒng)專(zhuān)業(yè)碩士研究生學(xué)位論文答辯 2008年6月31號(hào)1一、緒論 1.研究背景及意義 石羊河流域地處西北內(nèi)陸干旱區(qū)，是我國(guó)典型的資源型缺水地區(qū)，人均水資源占有量?jī)H775立方米，耕地畝均水資源占有量?jī)H280立方米，是河西走廊最早開(kāi)采利用地下水的地區(qū)。 近二十年來(lái)，隨著流域內(nèi)經(jīng)濟(jì)社會(huì)的發(fā)展和人口的不斷增加，水資源供需矛盾日趨突出，水資源開(kāi)發(fā)利用嚴(yán)重超載，現(xiàn)已成為河西走廊人口密度最大、水資源最短缺、供需矛盾最突出、生態(tài)與環(huán)境惡化最嚴(yán)重的流域。 2二、研究?jī)?nèi)容國(guó)家政府高度重視多學(xué)科共同關(guān)注為流域水文地質(zhì)研究提供基礎(chǔ)資料，也為當(dāng)

2、地地下水資源管理、地下水開(kāi)發(fā)利用、工農(nóng)業(yè)供水、水利工程建設(shè)以及流域生態(tài)恢復(fù)建設(shè)等提供可靠的決策依據(jù) 32.地下水研究歷程 公元前7世紀(jì)公元1800年：感性認(rèn)識(shí)、現(xiàn)象描述和主觀推測(cè)。18001940年：水文學(xué)理論取得了突破性進(jìn)展，地下水水文學(xué)的學(xué)科地位得以確立。1940年后：以計(jì)算機(jī)技術(shù)、遙感、GIS和數(shù)學(xué)模擬方法等諸多方法的引入和應(yīng)用為主要標(biāo)志，地下水研究進(jìn)入了的快速發(fā)展階段。4二、研究?jī)?nèi)容3.研究現(xiàn)狀 地下水受自然因素與人為因素共同影響，具有強(qiáng)烈的時(shí)空變化特點(diǎn)，其動(dòng)態(tài)過(guò)程存在地區(qū)上的差異性以及隨時(shí)間變化的特點(diǎn)，特別是容易受到外界條件影響的淺層地下水，這種地區(qū)上的特征及多年變化比較明顯。 現(xiàn)階

3、段，數(shù)學(xué)模型法是研究地下水變化的主要方法，此類(lèi)方法大體分為確定性模型和隨機(jī)模型：5二、研究?jī)?nèi)容（1）確定性模型 由地下水運(yùn)動(dòng)微分方程、控制方程、初始條件和邊界條件構(gòu)成的數(shù)學(xué)模型，常用解析法、物理法、數(shù)值法進(jìn)行模擬求解。（2）隨機(jī)模型 地下水動(dòng)態(tài)受多種因素的影響，包括氣象、水文、地質(zhì)、人為活動(dòng)等，可以用概率統(tǒng)計(jì)分析方法找出這些不確定性因素的規(guī)律，從而建立相應(yīng)的隨機(jī)模型。6二、研究?jī)?nèi)容模型比較：確定性模型 1）能比較準(zhǔn)確的預(yù)測(cè)地下水動(dòng)態(tài)變化2）基于許多假設(shè)條件基礎(chǔ)之上建立起來(lái)的，很難真實(shí)地模擬實(shí)際復(fù)雜的水文地質(zhì)條件。隨機(jī)模型1）無(wú)法表示出系統(tǒng)各個(gè)要素之間的動(dòng)力學(xué)關(guān)系2）對(duì)于缺乏地下水地質(zhì)資料或者含

4、水層參數(shù)資料不足，卻積累一定數(shù)量的水文氣象、開(kāi)采量和地下水動(dòng)態(tài)資料時(shí)，實(shí)用而有效。720世紀(jì)80年代開(kāi)始 熱衷于探討石羊河流域地下水資源的綜合開(kāi)發(fā)利用狀況，地下水動(dòng)態(tài)變化特征、水資源轉(zhuǎn)化等問(wèn)題。步入90年代后 引入先進(jìn)技術(shù)遙感、同位素技術(shù)； 地下水系統(tǒng)管理模式、水資源合理利用途徑、地下水來(lái)源及年齡、流域環(huán)境變化、水資源承載能力與可持續(xù)發(fā)展研究等方面作了相關(guān)研究。4.石羊河流域水資源研究現(xiàn)狀8二、研究?jī)?nèi)容2000年至今 1） 遙感、同位素技術(shù)被廣泛應(yīng)用于對(duì)探測(cè)地下水空間分布、地下水來(lái)源、流路、年齡等研究中； 2） 數(shù)值模擬軟件、一般統(tǒng)計(jì)學(xué)方法、地統(tǒng)計(jì)學(xué)方法被廣泛應(yīng)用于模擬地下水系統(tǒng)的動(dòng)態(tài)過(guò)程以及

5、時(shí)空分布變化，并進(jìn)一步探討影響地下水動(dòng)態(tài)變化的因子； 3）通過(guò)數(shù)學(xué)模型方法建立水資源評(píng)價(jià)指標(biāo)，來(lái)評(píng)價(jià)石羊河流域水資源承載能力、水資源脆弱性以及水資源危機(jī)等。9二、研究?jī)?nèi)容代表性的文獻(xiàn)有： 肖篤寧（2003） 地統(tǒng)計(jì)學(xué)和GIS結(jié)合，定量研究了民勤綠洲地下水埋深和礦化度的空間變異特征及其與土地利用變化的關(guān)系； 張若玲（2006） 利用水平衡方法、同位素技術(shù)、地球化學(xué)技術(shù)分析總結(jié)了整個(gè)地區(qū)的地下水動(dòng)態(tài)分布特征、同位素特征以及水化學(xué)分布規(guī)律； 孫艷偉（2007） 采用綜合指數(shù)法、灰色關(guān)聯(lián)投影法和基于MATLAB的人工神經(jīng)網(wǎng)絡(luò)，并結(jié)合地GIS對(duì)流域內(nèi)各個(gè)分區(qū)的地下水脆弱性進(jìn)行評(píng)價(jià)； Zhenmei Ma

6、（2008） 通過(guò)統(tǒng)計(jì)學(xué)方法分析了石羊河流域的徑流量變化趨勢(shì)和改變點(diǎn)，評(píng)估了氣候變化和人類(lèi)活動(dòng)對(duì)流域徑流的影響。10二、研究?jī)?nèi)容 小結(jié)： 大致包括以下研究?jī)?nèi)容： （1）地下水文要素的動(dòng)態(tài)變化特征及時(shí)空分布； （2）地下水資源開(kāi)發(fā)利用對(duì)生態(tài)環(huán)境的影響問(wèn)題 （3）流域水資源分布特征及水資源轉(zhuǎn)化規(guī)律的研究； （4）地下水資源管理模式優(yōu)化及可持續(xù)利用評(píng)價(jià) （5）氣候變化和人類(lèi)活動(dòng)對(duì)水資源變化的影響。115本文研究?jī)?nèi)容及方法時(shí)域動(dòng)態(tài)變化特征一般統(tǒng)計(jì)學(xué)方法空間分布特征地統(tǒng)計(jì)學(xué)方法綠洲地下潛水水系空間分布基于數(shù)字地下水位模型的地表徑流模擬12312二、研究區(qū)概況13 石羊河流域水資源比較貧乏，屬資源性缺水區(qū)

7、。從90年代末期開(kāi)始，地下水就已經(jīng)逐步成為石羊河流域的主要供水水源，年平均開(kāi)采量占全流域水資源利用總量的50%以上。 14三、石羊河流域地下水位年際變化特征1.原理 根據(jù)地下水水位數(shù)據(jù)的特征，一個(gè)非穩(wěn)定的地下水時(shí)間序列的組合模型，可表達(dá)為: h(t)=T(t)+P(t)+R(t) 其中:h(t)為地下水水位；T(t)為趨勢(shì)函數(shù)；P(t)為周期函數(shù)；R(t)為隨機(jī)特征函數(shù)；t為時(shí)間。151.1趨勢(shì)及突變點(diǎn)檢驗(yàn)趨勢(shì)檢驗(yàn)Mann-Kendall檢驗(yàn) 原假設(shè)H0:時(shí)間序列數(shù)據(jù)(x1，xn)是n個(gè)獨(dú)立的、隨機(jī)變量同分布的樣本 ；備擇假設(shè)H1是雙邊檢驗(yàn):對(duì)于所有的k，jn且kj，xk和xj的分布是不相同的

8、，Mann-Kendall檢驗(yàn)的統(tǒng)計(jì)變量S被定義為： 式中，分別為j，k年相應(yīng)的測(cè)量值。且kj，16當(dāng)n10時(shí)，標(biāo)準(zhǔn)的正態(tài)統(tǒng)計(jì)變量Z可用 下式計(jì)算：17突變點(diǎn)檢驗(yàn)Pettitt 檢驗(yàn) 水文序列變點(diǎn)檢測(cè)與識(shí)別是上世紀(jì)70年代以來(lái)統(tǒng)計(jì)技術(shù)方面的研究熱點(diǎn)，主要用以確定水文序列中突然出現(xiàn)變化的點(diǎn)，將時(shí)間序列劃分為差異明顯的各子時(shí)間段。 18則最顯著的突變點(diǎn)為:突變點(diǎn)的顯著性水平可以通過(guò)下式計(jì)算得出: 191.2時(shí)間序列模型 一個(gè)非穩(wěn)定的地下水時(shí)間序列的組合模型是由趨勢(shì)函數(shù)T(t)、周期函數(shù)P（t）和隨機(jī)特征函數(shù)R(t)線性組合而成，模型的每一項(xiàng) 由以下函數(shù)構(gòu)建：（1）趨勢(shì)擬合 20（2）周期性分析 分

9、離趨勢(shì)項(xiàng)后，對(duì)剩余的序列進(jìn)行周期分析，主要采用諧波分析法識(shí)別和提取周期項(xiàng)P(t)。（3）自相關(guān)特征分析 消除趨勢(shì)項(xiàng)和近似周期項(xiàng)后的剩余序列為平穩(wěn)隨機(jī)系列項(xiàng)，可直接對(duì)其用自回歸模型求解 R(t)=0+1R(t-1)+ 2R(t-2)+pR(t-p) 212.數(shù)據(jù)說(shuō)明223結(jié)果分析3.1地下水位年際變化特征及原因分析變化趨勢(shì)檢驗(yàn)結(jié)果23 根據(jù)上述六個(gè)站點(diǎn)地下水位隨時(shí)間變化過(guò)程，可將其劃分為以下四種模式： （1）單調(diào)下降型 （2）波動(dòng)衰減型 （3）中末期加速型 （4）末期趨穩(wěn)型變化模式分類(lèi)24原因分析25 自1984年至2005年，分別由2.02和2.22億立方米，減少到0.78和1.07億立方米,

10、平均每年都以0.05億立方米的速度逐年減少，整體也呈現(xiàn)逐年下降趨勢(shì)。 263.2地下水位階段變化特征及原因分析突變點(diǎn)檢測(cè)結(jié)果27分階段地下水位下降速度和下降加速度變化圖28變點(diǎn)年原因分析29地下水位年際變化的模擬 R(t)=0+1R(t-1)+ 2R(t-2)+pR(t-p) 30擬合結(jié)果31四、石羊河流域地下水位空間分布特征1.原理 1.1數(shù)字地下水位模型的概念 本文所建立的數(shù)字地下水位模型（Digital Watertable Model，簡(jiǎn)稱(chēng)DWTM)，采用規(guī)則格網(wǎng)數(shù)據(jù)結(jié)構(gòu)，該格網(wǎng)有多個(gè)大小相同且互相相鄰的單元格組成，每一單元格的行列編號(hào)表示該單元格所對(duì)應(yīng)位置的X、Y坐標(biāo)，而單元格的值表

11、則示該X、Y坐標(biāo)點(diǎn)的地下水位。321.2協(xié)同克立格空間插值方法 協(xié)同克里格法利用多個(gè)變量之間的互相關(guān)性，用這些變量對(duì)其中一個(gè)變量進(jìn)行局部估計(jì)的方法，但是只有當(dāng)兩個(gè)變量之間相關(guān)性較大時(shí)，建立和計(jì)算交叉協(xié)方差函數(shù)和變異函數(shù)才有意義，協(xié)同克里格估計(jì)才能夠保證精確和有效。其定義為： 設(shè)點(diǎn)處第個(gè)變量為待估變量，則其估計(jì)量是K個(gè)協(xié)同區(qū)域化變量樣本點(diǎn)數(shù)據(jù)的線性組合：332 數(shù)據(jù)說(shuō)明343結(jié)果分析（1）地下水位的空間分布特征3536空間分布特征： 地下水位的逐級(jí)遞減方向與地表水系方向相一致，在山前地帶地表水系與地下水聯(lián)系緊密，地下水位相對(duì)較高，但在河流下游地區(qū)以及無(wú)地表徑流的地區(qū)地下水位則較低。 地下水位的變化存在明顯的分帶特征，民勤紅崖山水庫(kù)前的地下水位變化帶最為明顯，這是由于受阿拉古山和紅崖山的阻隔，民勤紅崖山水庫(kù)前后分屬于兩個(gè)不同的盆地，致使其前后差異明顯。 地下水位在2001年的空間分布狀態(tài)與1984年相比，民勤縣東北部地區(qū)的地下水位低值區(qū)，逐漸向西南方向移動(dòng)，地下水位低值區(qū)的面積有所擴(kuò)大，這與民勤為緩解當(dāng)?shù)厮Y源供需矛盾，大量開(kāi)采地下水的人為活動(dòng)是密不可分的。37（2）地下水位下降幅度的空間分布特征38（3）地下水位空間異質(zhì)性的逐年變化特征39（4）地下水