石 朋,侯爰冰,馬欣欣,陳 喜,張志才

(1.河海大學(xué)水文水資源與水利工程科學(xué)國家重點實驗室,江蘇 南京 210098;

2.河海大學(xué)水文水資源學(xué)院,江蘇 南京 210098)

“南石(石漠化)北沙(沙漠化)”一直是制約我國西部地區(qū)可持續(xù)發(fā)展的兩大生態(tài)環(huán)境問題。西南喀斯特地區(qū)面積約54萬km2,包括以貴州為中心的廣西、云南 、四川、湖南、湖北 、廣東 、江西、重慶八省(市),是全球三大喀斯特集中分布區(qū)中連片裸露碳酸鹽巖面積最大的區(qū)域。該區(qū)內(nèi)山地面積大,降水豐沛,土層淺薄,可溶巖成土速率緩慢,水土流失危險度高,屬非地帶性的脆弱生態(tài)帶[1]。同時,該區(qū)內(nèi)人口達1億多,貧困人口相對集中,人地矛盾非常突出。坡地植被一旦被破壞,土壤侵蝕作用將會加劇,往往導(dǎo)致薄土層全部流失,造成嚴(yán)重的石漠化,水分、養(yǎng)分調(diào)蓄能力迅速降低。目前的石漠化土地面積已達 17萬 km2,占該區(qū)喀斯特面積的 31%,近2000萬人的生存環(huán)境受到嚴(yán)重威脅。更嚴(yán)重的是,石漠化面積還在不斷擴大,不僅使土地生產(chǎn)力下降,地表植被覆蓋率銳減,系統(tǒng)水源涵養(yǎng)能力削弱,地表水源枯竭,更造成土地資源喪失,糧食減產(chǎn)。另外,該區(qū)地處長江和珠江兩大水系的上游,一旦生態(tài)屏障崩潰,將嚴(yán)重危及中國半壁江山的生態(tài)安全。因此,西南喀斯特流域水循環(huán)研究一直是國家的科研重點。

1 西南喀斯特流域的水文特性

1.1 喀斯特流域的空間結(jié)構(gòu)特征

由可溶巖組成的非均質(zhì)含水介質(zhì),構(gòu)成了巖溶喀斯特流域所特有的二元三維空間結(jié)構(gòu)系統(tǒng)。該系統(tǒng)在水流侵蝕和溶蝕的雙重作用下,由地表、地下 2個子系統(tǒng)構(gòu)成,其中,地表子系統(tǒng)包括溶溝、溶隙、峰叢、天窗、漏斗及落水洞等;地下子系統(tǒng)則主要由溶洞、管道、伏流及地下河等組成。地表成為喀斯特徑流中的發(fā)生場和分配場,地下則成為喀斯特徑流中的輸送場和調(diào)蓄場,而最終會通過地表河或地下河以較大的落差匯入江河。在巖溶系統(tǒng)的整個發(fā)展過程中,3種結(jié)構(gòu)的巖溶含水層或含水體系得到了普遍認(rèn)可,分別是:①細小結(jié)構(gòu)體系,以封閉的空隙占優(yōu)勢;②粗大結(jié)構(gòu)體系,以基本上未封閉的空隙占優(yōu)勢;③復(fù)活體系。因此,喀斯特流域在結(jié)構(gòu)上有4個基本特性[2],即不均一的雙重含水介質(zhì)結(jié)構(gòu)、二元流場形態(tài)結(jié)構(gòu)、三維空間地域結(jié)構(gòu)以及功能上的耗散結(jié)構(gòu)。

1.2 喀斯特流域的地貌水文效應(yīng)

長期以來,人們對喀斯特流域不同于常態(tài)的水文和地貌問題做了大量的研究工作,但大多只側(cè)重于單方面的研究。水文問題和地貌問題之間有著十分密切的聯(lián)系,從長的時間尺度來看,喀斯特水是喀斯特地貌演化的直接動力;而從短的時間尺度來看,地貌形態(tài)幾乎是穩(wěn)定的,流域內(nèi)不同的地貌組合是導(dǎo)致區(qū)內(nèi)特殊水文過程的根本動因。近年來,國內(nèi)外的一些學(xué)者也認(rèn)識到這一點,并注重于水文地貌系統(tǒng)的研究,把喀斯特流域的水文地貌看作是一個相互作用、相互影響、互為因果的有機系統(tǒng),不但從水的角度出發(fā)考慮地貌問題,同時也從地貌的角度出發(fā)研究水文現(xiàn)象,進而促進了水文學(xué)和地貌學(xué)的交叉學(xué)科——水文地貌學(xué)的興起。

整個喀斯特流域地貌大系統(tǒng)中除了有構(gòu)造地貌系統(tǒng)、巖石地貌系統(tǒng)和氣候地貌系統(tǒng)外,還有水文地貌系統(tǒng)。水文與地貌這兩大要素之間相互制約、相互影響,一方面在不同的水文條件下,發(fā)育了不同的喀斯特地貌;另一方面,不同的喀斯特地貌也制約著不同的水文過程。

在喀斯特地貌系統(tǒng)中,不同的水文過程作用往往使喀斯特系統(tǒng)在時空上表現(xiàn)出不同的水文地貌效應(yīng)(如以水流侵蝕為主的地貌結(jié)構(gòu)或以水流溶蝕為主的地貌結(jié)構(gòu));與此相應(yīng)的,喀斯特地貌的多樣性和差異性必然也會導(dǎo)致流域水文過程不同。一般來說,喀斯特流域地貌結(jié)構(gòu)對水體過程的作用可歸納為3類:蓄水作用、滯水作用和導(dǎo)水作用。表層溶蝕裂隙的入滲強度大、厚度小,往往表現(xiàn)為滯水作用;在峰林平原、峰林盆地中,土層厚度較大,具有較強的蓄水作用;而諸如地表、地下河及谷地等,則往往表現(xiàn)為導(dǎo)水作用。

1.3 喀斯特流域的地下水運動特征

在我國,西南地區(qū)的喀斯特巖溶介質(zhì)比北方地區(qū)的具有更強的非均質(zhì)性和各向異性。大氣降水和地表水通過巖溶地表滲入地下轉(zhuǎn)化為地下水,對巖溶地下結(jié)構(gòu)不斷進行改造,從而逐漸形成巖溶通道(管道、洞穴等),并有序構(gòu)成巖溶地下河系[3]。在這樣的地下河系統(tǒng)中,水流以管道流為主,巖溶管道或地下河構(gòu)成地下徑流的主干。在2個管道流之間,地下水在層理、節(jié)理和溶孔中運動,此時的運動狀態(tài)與一般的滲流相似。管道流由于流動阻力小,流速非常大。據(jù)有關(guān)資料[4]統(tǒng)計,一般情況下地下河發(fā)育區(qū)的巖溶水流速枯水期為 8.64～17.28km/d,平水期為17.28～ 43.20km/d,洪水期為 43.20～ 129.60km/d,極端流速可達172.80km/d,可見在管道流或地下河發(fā)育區(qū)的地下水運動已不符合達西定律。此外,巖溶地下河系統(tǒng)中還包含巖溶裂隙、巖溶裂縫和巖溶孔隙等多種介質(zhì)體,形成了一個多重的復(fù)合體系,具有高度的非均質(zhì)性,使得地下河系統(tǒng)中的水流運動常呈現(xiàn)出達西流與非達西流并存的現(xiàn)象[5],這也為如何準(zhǔn)確描述巖溶地下水的運動規(guī)律增加了困難。

2 喀斯特流域水循環(huán)的研究方法

2.1 流域水文模型方法

過去對巖溶喀斯特流域的水循環(huán)過程往往偏重于定性研究,直到 20世紀(jì)90年代以后才陸續(xù)出現(xiàn)一些針對巖溶流域水循環(huán)過程定量研究的水文模型。我國西南巖溶流域具有特殊的二元三維地質(zhì)結(jié)構(gòu)的極其復(fù)雜的地下河系統(tǒng)。在長期風(fēng)化侵蝕作用下,喀斯特流域中大量的裂隙、溶溝、漏斗、天窗和落水洞使得其水文循環(huán)極其復(fù)雜。降雨時,地面徑流通過這些大型管道迅速進入地下成為地下徑流,使一次雨洪過程中地下徑流占較大比重,這與非巖溶流域的特征不同,加之多重復(fù)合體系的多種水流運動介質(zhì)體及達西流與非達西流并存的現(xiàn)象,給巖溶流域建立水文模型帶來相當(dāng)?shù)睦щy。

目前在巖溶喀斯特地區(qū)建立的水文模型大體可概括為以下3類:

第1類通過對在非巖溶流域上成功使用的概念性集總式流域水文模型進行適度改造,使其在一定程度上能夠適用于巖溶地區(qū)。如袁道先[6]在水箱模型的基礎(chǔ)上開發(fā)出的桂林“丫吉模型”是最早在西南巖溶地區(qū)成功建立的流域水文模型之一。此類模型[7-12]雖然在模型構(gòu)建時對流域進行了一定的單元劃分,但由于對水文過程、輸入變量及邊界條件的概化,對流域空間變異性的忽略,以及對喀斯特流域產(chǎn)匯流空間分布的不明確,使得該類模型從嚴(yán)格意義上說仍屬于黑箱模型,且在實際應(yīng)用中存在諸多問題。當(dāng)將此類集總式模型應(yīng)用于地形地貌空間變異性高度發(fā)育的喀斯特流域時,會發(fā)現(xiàn)模型與客觀實際嚴(yán)重偏離,從而影響模型的計算精度和可信程度。

第2類是結(jié)合灰色理論、分形方法及人工神經(jīng)網(wǎng)絡(luò)等現(xiàn)代數(shù)學(xué)方法所建立起來的數(shù)理統(tǒng)計模型[13-14]。束龍倉等[15]將門限自回歸模型與小波BP神經(jīng)網(wǎng)絡(luò)相結(jié)合應(yīng)用于后寨河流域出口流量的預(yù)測;石朋等[16-17]基于線性矩法和MK等方法分析了貴州巖溶流域徑流過程的年度和季節(jié)變化、洪枯水變化規(guī)律及其演變的空間特征;陳宏峰等[18]利用BP人工神經(jīng)網(wǎng)絡(luò)建立了湖南洛塔河地區(qū)表層巖溶泉的徑流預(yù)測模型等。以數(shù)學(xué)方法為基礎(chǔ)的數(shù)理統(tǒng)計模型往往需要大量長系列的歷史觀測資料,然而,我國西南巖溶流域處于經(jīng)濟發(fā)展相對落后的地區(qū),除少數(shù)試驗流域近年來在國家重大研究項目的資助下開展了一些連續(xù)觀測外,大部分地區(qū)資料匱乏甚至無資料,從而使統(tǒng)計模型在理論上難以獲得支撐。由于此類方法缺乏必要的物理基礎(chǔ),所建模型用于模擬尚可,用于預(yù)測則存在較大的不確定性,因此在巖溶水文循環(huán)研究中很難有進一步發(fā)展的空間。

第3類是非巖溶流域上適用的分布式/半分布式水文模型,嘗試將其應(yīng)用于巖溶流域。與集總式模型不同,分布式模型能夠同時考慮流域水文過程和植被、土壤、地形地貌等地理因素的時空變異性及其對水文過程的影響,一些固有的自然地理要素往往被作為模型要素加以考慮,因此,分布式模型能夠較好地描述流域內(nèi)細部的水文過程。任啟偉[19]在SWAT模型基礎(chǔ)上改建的刁江半分布式模型可視為分布式模型在我國西南喀斯特流域上的有益嘗試。張志才等[20]在DHSVM模型基礎(chǔ)上所建立的巖溶試驗流域分布式模型以及蒙海花等[21]基于落水洞所構(gòu)建的一個半分布式水文模型也有一定的計算精度。在建模理論上分布式水文模型是適用于復(fù)雜條件下的西南喀斯特流域的,但也存在一些亟待解決的問題:①巖溶產(chǎn)流模式的確定。由于表層巖溶帶的存在使得巖溶流域的產(chǎn)流模式明顯不同于非巖溶流域的產(chǎn)流模式[22-24],表層巖溶帶裂隙率高,滲透性好,地表徑流只出現(xiàn)在連續(xù)降雨且雨強極大的降雨過程中[25]。因此,需要明確表層巖溶帶在降水入滲補給和產(chǎn)流中的作用[26-27]。②復(fù)雜下墊面條件的模型描述。落水洞等巖溶洼地往往是地表水快速匯集進入地下水系的連接處[28]。分布式水文模型根據(jù)DEM采用填洼法進行柵格匯流演算,不能將該模型直接照搬到喀斯特流域,而應(yīng)在喀斯特地貌識別基礎(chǔ)上建立地表與地下水水系的連接關(guān)系。③巖溶徑流成分及匯流路徑的確定。目前柵格型匯流主要根據(jù)地形坡度進行計算,而產(chǎn)生于表層巖溶帶的水流不僅受地形影響,還取決于巖溶裂隙帶滲透系數(shù)的各向異性。因此,喀斯特產(chǎn)匯流受喀斯特地貌、土壤、喀斯特裂隙及地表和地下水系控制,目前尚無可直接套用或引用的模式,需要在詳細調(diào)查、巖溶水動態(tài)加密監(jiān)測基礎(chǔ)上開展深入研究[29]。④下墊面空間異質(zhì)性導(dǎo)致的尺度效應(yīng)問題。基于小尺度建立的物理模型及數(shù)值解在流域尺度上由于維數(shù)災(zāi)、參數(shù)災(zāi),增加了模型的不確定性,降低了模型的適用性。需要關(guān)注從基本尺度得到的基本規(guī)律如何應(yīng)用到流域及區(qū)域尺度上。就目前而言,分布式模型在喀斯特地區(qū)的應(yīng)用尚不多見,仍需要更多的實踐檢驗。

2.2 室內(nèi)外試驗方法

針對西南喀斯特流域下墊面條件的高度空間異質(zhì)性,眾多學(xué)者開始關(guān)注不同類型下墊面條件的水文特性,從而開展了大量的室內(nèi)外試驗。此類試驗大體上可歸納為2類:一是對局地水文地質(zhì)參數(shù)的求證;二是對地表以下水流路徑的探測。

在水文地質(zhì)參數(shù)的求證方面,章程等[30-31]通過野外溶蝕試片法證明了不同土地利用類型下溶蝕量存在顯著差異性,同時還利用多參數(shù)自動記錄儀探討了巖溶泉的水位、水溫等與降雨的響應(yīng)關(guān)系;張志才等[32]利用時域反射儀測定山體不同位置在不同植被覆蓋條件下的土壤含水率,分析了喀斯特峰叢山體土壤水的空間分布規(guī)律及其影響因素;劉延惠等[33]在巖溶流域設(shè)置徑流小區(qū),通過實地觀測分析不同植被類型下巖溶地表徑流的差異;王臘春等[34]通過脈沖試驗和降雨天然脈沖分析,論述了后寨河流域巖溶水水文特征與流域貯水結(jié)構(gòu)、地質(zhì)地貌結(jié)構(gòu)的相互關(guān)系,以及裂隙流與管道流的轉(zhuǎn)化關(guān)系。魯程鵬等[35]利用電導(dǎo)率與溶質(zhì)濃度的正相關(guān)關(guān)系,在示蹤試驗基礎(chǔ)上提出一套估算地下水流速及含水層滲透系數(shù)的方法;鄧振平等[36]采用弱酸性紅A熒光燃料類試劑作為示蹤材料,分析了巖溶落水洞和地下河出口的連通性,以及地下分水嶺的分布特征。上述試驗研究對非均一含水介質(zhì)巖溶水運動研究及巖溶水文模擬具有一定的參考意義,但考慮到喀斯特流域的特殊情況,這些試驗的結(jié)論是否具有代表性,能否在流域尺度上加以推廣值得商榷。

在地表以下水流路徑的探測方面,目前國內(nèi)尚少有此方面的研究性論文,現(xiàn)有的方法仍是以水文地質(zhì)鉆探方法為主。此類方法工程量較大,在一般的研究流域難以實施。近年來,連續(xù)電導(dǎo)率剖面儀及地質(zhì)雷達的出現(xiàn)在一定程度上提高了效率,但同時仍存在一些問題。以RH-4型連續(xù)電導(dǎo)率剖面儀為例,其在地表以下 10～500m能夠精細地反演地下情況,但自地表至地表以下10m為其探測盲區(qū),而這10m范圍恰恰是西南喀斯特流域中表層巖溶帶強烈發(fā)育的區(qū)段。再如地質(zhì)雷達,其優(yōu)點是靈活便攜,且能在各種復(fù)雜地形條件下正常使用,但缺點是易受干擾,外界微弱的電磁信號即會影響其工作結(jié)果。此外,地表土層厚度、土壤含水量等均會影響其探測結(jié)果,而對于探測結(jié)果的反演也是見仁見智,并無嚴(yán)格統(tǒng)一的標(biāo)準(zhǔn),這樣一來,對一些細小的或由破碎巖石構(gòu)成的地下水流通道,往往難以被識別出來。

3 喀斯特水循環(huán)研究的發(fā)展方向

3.1 巖溶環(huán)境下產(chǎn)流過程的研究

巖溶流域特殊的地質(zhì)結(jié)構(gòu)使得其下墊面條件不同于非巖溶流域。表層巖溶帶的存在使其包氣帶結(jié)構(gòu)更為復(fù)雜。那么,如何在此類特殊包氣帶條件下探討其產(chǎn)流機制和產(chǎn)流模式,將是未來巖溶流域水循環(huán)研究的一個重點。盡管已有的研究表明,巖溶流域雖土層淺薄,但仍屬蓄滿產(chǎn)流模式,但此種蓄滿不僅僅是土壤的蓄滿,同時需要考慮表層巖溶帶的蓄滿,而表層巖溶帶內(nèi)的水流運移儲存方式又不同于土壤,那么,如何構(gòu)建出這種“土壤+表層巖溶帶”的復(fù)合蓄滿產(chǎn)流模式將是一個難點。

3.2 巖溶環(huán)境下匯流過程的研究

巖溶地下河系統(tǒng)中的水流運動常呈現(xiàn)出達西流與非達西流并存且地表河與地下河交替出現(xiàn)的情況。對于這種復(fù)雜的水流情態(tài),如何通過合理的數(shù)學(xué)物理方法加以描述將是巖溶環(huán)境下匯流過程研究的主要難點。在此類環(huán)境下,除非在某些重點區(qū)域進行鉆井等強化觀測,對于大部分地區(qū),由于難以確定邊界條件,傳統(tǒng)的水動力學(xué)方法顯然不再適用,如單純使用線性水庫進行演算,則又回歸到黑箱模型中。因此,如何在其中找到一個合適的切入點顯得尤為重要。目前,將分形技術(shù)用于模擬復(fù)雜的地下河系統(tǒng)與管道裂隙系統(tǒng)并概化出理想地下河系統(tǒng)響應(yīng)函數(shù)是一種有益的嘗試。

3.3 喀斯特流域水循環(huán)過程的量化研究

由于水文模型能夠給出流域降雨徑流的變化過程,因此,在喀斯特流域開展流域水文模型研究仍將是今后的主要研究內(nèi)容之一。但水文模型的研究方法需要與研究側(cè)重點相適應(yīng),如研究重點在于考慮較大尺度上流域或區(qū)域的水資源量問題,則宜采用結(jié)構(gòu)簡單又不完全失其物理意義的概念性集總式水文模型;當(dāng)以面積較小的試驗流域為研究對象,重點關(guān)注巖溶水循環(huán)過程的具體細節(jié)時,需要構(gòu)建分布式巖溶流域水文模型,通過劃分網(wǎng)格或代表性單元,盡可能真實地反映降雨徑流的時空變化?？紤]到目前的具體情況,采用一些半分布式模型并加以改進,也不失為一種折中的方法。

3.4 喀斯特流域的地貌水文學(xué)研究

地貌水文學(xué)研究在非巖溶流域中一直是地貌學(xué)家和水文學(xué)家都十分關(guān)注的領(lǐng)域,而在巖溶流域中,此方面的研究更值得重點關(guān)注。不同的巖溶地貌類型及特殊地質(zhì)單元,諸如峰叢、峰林、天窗、落水洞等,對產(chǎn)流和匯流都會造成巨大的影響。如果能夠建立某種模式,通過地形地貌來確定不同巖溶類型下的產(chǎn)匯流特征,對于巖溶環(huán)境下的水循環(huán)研究將具有劃時代的意義。

3.5 喀斯特流域的生態(tài)系統(tǒng)研究

喀斯特空間變異性及地表-地下雙重結(jié)構(gòu)對水文-生態(tài)鏈接具有重要控制作用。近年來,人類非理性的活動對西南喀斯特山地生態(tài)水文格局的破壞十分嚴(yán)重,地表生態(tài)水文過程已經(jīng)喪失了協(xié)調(diào)發(fā)展的趨勢。因此,選擇西南喀斯特典型流域,針對地貌、植被高度空間異質(zhì)性和地表-地下雙重結(jié)構(gòu)水文地質(zhì)特征研究流域尺度的水循環(huán)動力過程對生態(tài)環(huán)境的影響,對喀斯特流域生態(tài)系統(tǒng)的恢復(fù)和重建具有重要意義,勢必將成為未來重點研究內(nèi)容之一。

3.6 新技術(shù)新方法的應(yīng)用研究

盡管現(xiàn)有的地質(zhì)雷達、連續(xù)電導(dǎo)率剖面儀等設(shè)備應(yīng)用于巖溶流域時還存在諸多問題,但弄清巖溶地下管道的具體位置及其連接情況,對于研究巖溶地下水的運移具有十分重要的意義。因此,在上述儀器設(shè)備基礎(chǔ)上,進一步開發(fā)適用于巖溶地下勘探的新技術(shù)新方法顯得尤為重要。隨著3S技術(shù)的不斷進步,將3S技術(shù)應(yīng)用到巖溶流域特殊地質(zhì)結(jié)構(gòu)的識別,可能是今后一個值得關(guān)注的發(fā)展方向。

4 結(jié) 語

西南巖溶流域地處我國長江和珠江兩大水系的上游,具有重要的生態(tài)和環(huán)境意義,巖溶流域水循環(huán)研究對于該地區(qū)的環(huán)境保護和可持續(xù)發(fā)展具有十分重要的意義。在總結(jié)現(xiàn)有國內(nèi)外學(xué)者對我國西南巖溶流域水循環(huán)相關(guān)研究的基礎(chǔ)上,指出了目前研究中存在的一些問題,并對未來可能的研究方向提出了一些淺見。

:

[1]李陽兵,王世杰,容麗,等.西南喀斯特山地石漠化及生態(tài)恢復(fù)研究展望[J].生態(tài)學(xué)雜志,2004,23(6):84-88.

[2]梁虹,楊明德.喀斯特流域水文地貌系統(tǒng)及其識別方法探討[J].中國巖溶,1994,13(1):1-9.

[3]郭純青.中國巖溶地下河系及其水資源[J].水文地質(zhì)工程地質(zhì),2001(5):43-45.

[4]朱遠峰.中國巖溶水系統(tǒng)和巖溶水資源研究進展[C]//工程地質(zhì)、水文地質(zhì)、環(huán)境地質(zhì)論文集.北京:地震出版社,1993:194-199.

[5]郭琳,陳植華.巖溶地區(qū)地下河系統(tǒng)水資源定量評價的問題與出路[J].中國巖溶,2006,25(1):1-6.

[6]袁道先.中國南方裸露型巖溶峰叢山區(qū)巖溶水系統(tǒng)及其數(shù)學(xué)模型的研究:以桂林丫吉村為例[M].南寧:廣西師范大學(xué)出版社,1995.

[7]盧衛(wèi)中.巖溶流域水文模型初步研究[J].水文,1995(2):29-33.

[8]王臘春,史運良.非閉合流域巖溶水模型[J].水科學(xué)進展,1995,6(4):318-323.

[9]史運良,王臘春.巖溶水概念分布模型[J].南京大學(xué)學(xué)報:地理版,1992(13):134-161.

[10]程根偉.新安江巖溶水文模型[J].水電能源科學(xué),1991,9(2):139-144.

[11]FLEURY P,PLAGNES V,BAKALOWICZ M.Modelling of the functioning ofkarstaquiferswith a reservoirmodel:applicationto fortaine de Vaucluse(South of France)[J].Journal of Hydrology,2007,345:38-49.

[12]RIMMER A,SALINGAR Y.Modelling precipitation-streamflow processes in karstbasin:the case of the Jordan River sources,Israel[J].Journal of Hydrology,2006,331:524-542.

[13]COMTE V B,JOUDE H,ROESCH A,et al.Time series analyses for karst river interactions assessment:case of the Coulazou River(Southern France)[J].Journal of Hydrology,2008,349:98-114.

[14]張保祥.基于神經(jīng)網(wǎng)絡(luò)和遺傳算法的濟南市區(qū)巖溶地下水預(yù)報模型研究[J].山東農(nóng)業(yè)大學(xué)學(xué)報,2004,35(3):436-441.

[15]束龍倉,董貴明,陶玉飛,地下河天窗水位變化分析及預(yù)測[J].水利學(xué)報,2009,40(5):529-534.

[16]SHI Peng,CHEN Xi,QU Si-min,et al.Regional frequency analysis of low flow based on L-moments:case study in karst area,southwest China[J].Journal of Hydrologic Engineering,2010,15(5):370-377.

[17]楊濤,陳喜,謝海文,等.貴州近50年來的徑流過程時空演變特征及驅(qū)動因素探討[J].地球與環(huán)境,2009,37(3):207-213.

[18]陳宏峰,朱明秋,夏日元.湖南洛塔干河豬場表層巖溶泉BP人工神經(jīng)網(wǎng)絡(luò)分析[J].中國巖溶,2005,24(4):300-304.

[19]任啟偉.基于改進SWAT模型的西南巖溶流域水量評價方法研究[D].武漢:中國地質(zhì)大學(xué),2006.

[20]張志才,陳喜,石朋,等.喀斯特流域分布式水文模型及植被生態(tài)水文效應(yīng)[J].水科學(xué)進展,2009,20(6):806-811.

[21]蒙?；?王臘春,蘇維詞,等.基于落水洞的巖溶半分布式水文模型的構(gòu)建及其應(yīng)用[J].地理科學(xué),2009,29(4):550-554.

[22]賀中華,梁虹,黃法蘇,等.基于DEM 的喀斯特流域地貌類型的識別[J].大地測量與地球動力學(xué),2008,28(3):46-53.

[23]蔣忠誠,袁道先.表層巖溶帶的巖溶動力學(xué)特征及其環(huán)境和資源意義[J].地球?qū)W報,1999,20(3):302-308.

[24]鄒勝章,張文慧,梁小平,等.表層巖溶帶調(diào)蓄系數(shù)定量計算:以湘西洛塔趙家灣為例[J].水文地質(zhì)工程地質(zhì),2005(4):37-42.

[25]彭韜,王世杰,張信寶,等.喀斯特坡地地表徑流系數(shù)監(jiān)測初報[J].地球與環(huán)境,2008,36(2):125-129.

[26]姜光輝,吳吉春,郭芳,等.森林覆蓋的喀斯特地區(qū)表層巖溶帶的產(chǎn)流閾值[J].水科學(xué)進展,2008,19(1):72-77.

[27]姜光輝,郭芳.我國西南巖溶區(qū)表層巖溶帶的水文動態(tài)分析[J].水文地質(zhì)工程地質(zhì),2009(5):89-93.

[28]PERRIN J,JEANNIN P,ZWAHLEN F.Epikarst storage in a karst aquifer:a conceptual model based on isotopic data,Milandre test site,Switzerland[J].Journal ofHydrology,2003,279:106-124.

[29]覃小群,蔣忠誠.表層巖溶帶及其水循環(huán)的研究進展與發(fā)展方向[J].中國巖溶,2005,24(3):250-254.

[30]章程,謝運球,呂勇,等.不同土地利用方式對巖溶作用的影響:以廣西弄拉峰叢洼地巖溶系統(tǒng)為例[J].地理學(xué)報,2006,61(11):1181-1188.

[31]章程,袁道先,李強,等.典型表層巖溶動力系統(tǒng)的環(huán)境敏感性研究:以廣西弄拉峰叢洼地巖溶系統(tǒng)為例[J].水文地質(zhì)工程地質(zhì),2005(1):78-80.

[32]張志才,陳喜,石朋,等.貴州喀斯特峰叢山體土壤水分布特征及其影響因素[J].長江流域資源與環(huán)境,2008,17(5):803-807.

[33]劉延惠,張喜,崔迎春,等.貴州開陽喀斯特山地幾種不同植被類型的地表徑流研究[J].貴州林業(yè)科技,2005,33(2):8-10.

[34]王臘春,許有鵬,張立峰,等.貴州普定后寨地下河流域巖溶水特征研究[J].地理科學(xué),2000,20(6):557-562.

[35]魯程鵬,束龍倉,苑利波,等.基于示蹤試驗求解巖溶含水層水文地質(zhì)參數(shù)[J].吉林大學(xué)學(xué)報:地球科學(xué)版,2009,39(4):717-721.

[36]鄧振平,周小紅,何師意,等.西南巖溶石山地區(qū)巖溶地下水示蹤試驗與分析[J].中國巖溶,2007,26(2):163-169.