石 朋,侯爰冰,馬欣欣,陳 喜,張志才
(1.河海大學(xué)水文水資源與水利工程科學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室,江蘇 南京 210098;
2.河海大學(xué)水文水資源學(xué)院,江蘇 南京 210098)
“南石(石漠化)北沙(沙漠化)”一直是制約我國(guó)西部地區(qū)可持續(xù)發(fā)展的兩大生態(tài)環(huán)境問題。西南喀斯特地區(qū)面積約54萬km2,包括以貴州為中心的廣西、云南 、四川、湖南、湖北 、廣東 、江西、重慶八省(市),是全球三大喀斯特集中分布區(qū)中連片裸露碳酸鹽巖面積最大的區(qū)域。該區(qū)內(nèi)山地面積大,降水豐沛,土層淺薄,可溶巖成土速率緩慢,水土流失危險(xiǎn)度高,屬非地帶性的脆弱生態(tài)帶[1]。同時(shí),該區(qū)內(nèi)人口達(dá)1億多,貧困人口相對(duì)集中,人地矛盾非常突出。坡地植被一旦被破壞,土壤侵蝕作用將會(huì)加劇,往往導(dǎo)致薄土層全部流失,造成嚴(yán)重的石漠化,水分、養(yǎng)分調(diào)蓄能力迅速降低。目前的石漠化土地面積已達(dá) 17萬 km2,占該區(qū)喀斯特面積的 31%,近2000萬人的生存環(huán)境受到嚴(yán)重威脅。更嚴(yán)重的是,石漠化面積還在不斷擴(kuò)大,不僅使土地生產(chǎn)力下降,地表植被覆蓋率銳減,系統(tǒng)水源涵養(yǎng)能力削弱,地表水源枯竭,更造成土地資源喪失,糧食減產(chǎn)。另外,該區(qū)地處長(zhǎng)江和珠江兩大水系的上游,一旦生態(tài)屏障崩潰,將嚴(yán)重危及中國(guó)半壁江山的生態(tài)安全。因此,西南喀斯特流域水循環(huán)研究一直是國(guó)家的科研重點(diǎn)。
由可溶巖組成的非均質(zhì)含水介質(zhì),構(gòu)成了巖溶喀斯特流域所特有的二元三維空間結(jié)構(gòu)系統(tǒng)。該系統(tǒng)在水流侵蝕和溶蝕的雙重作用下,由地表、地下 2個(gè)子系統(tǒng)構(gòu)成,其中,地表子系統(tǒng)包括溶溝、溶隙、峰叢、天窗、漏斗及落水洞等;地下子系統(tǒng)則主要由溶洞、管道、伏流及地下河等組成。地表成為喀斯特徑流中的發(fā)生場(chǎng)和分配場(chǎng),地下則成為喀斯特徑流中的輸送場(chǎng)和調(diào)蓄場(chǎng),而最終會(huì)通過地表河或地下河以較大的落差匯入江河。在巖溶系統(tǒng)的整個(gè)發(fā)展過程中,3種結(jié)構(gòu)的巖溶含水層或含水體系得到了普遍認(rèn)可,分別是:①細(xì)小結(jié)構(gòu)體系,以封閉的空隙占優(yōu)勢(shì);②粗大結(jié)構(gòu)體系,以基本上未封閉的空隙占優(yōu)勢(shì);③復(fù)活體系。因此,喀斯特流域在結(jié)構(gòu)上有4個(gè)基本特性[2],即不均一的雙重含水介質(zhì)結(jié)構(gòu)、二元流場(chǎng)形態(tài)結(jié)構(gòu)、三維空間地域結(jié)構(gòu)以及功能上的耗散結(jié)構(gòu)。
長(zhǎng)期以來,人們對(duì)喀斯特流域不同于常態(tài)的水文和地貌問題做了大量的研究工作,但大多只側(cè)重于單方面的研究。水文問題和地貌問題之間有著十分密切的聯(lián)系,從長(zhǎng)的時(shí)間尺度來看,喀斯特水是喀斯特地貌演化的直接動(dòng)力;而從短的時(shí)間尺度來看,地貌形態(tài)幾乎是穩(wěn)定的,流域內(nèi)不同的地貌組合是導(dǎo)致區(qū)內(nèi)特殊水文過程的根本動(dòng)因。近年來,國(guó)內(nèi)外的一些學(xué)者也認(rèn)識(shí)到這一點(diǎn),并注重于水文地貌系統(tǒng)的研究,把喀斯特流域的水文地貌看作是一個(gè)相互作用、相互影響、互為因果的有機(jī)系統(tǒng),不但從水的角度出發(fā)考慮地貌問題,同時(shí)也從地貌的角度出發(fā)研究水文現(xiàn)象,進(jìn)而促進(jìn)了水文學(xué)和地貌學(xué)的交叉學(xué)科——水文地貌學(xué)的興起。
整個(gè)喀斯特流域地貌大系統(tǒng)中除了有構(gòu)造地貌系統(tǒng)、巖石地貌系統(tǒng)和氣候地貌系統(tǒng)外,還有水文地貌系統(tǒng)。水文與地貌這兩大要素之間相互制約、相互影響,一方面在不同的水文條件下,發(fā)育了不同的喀斯特地貌;另一方面,不同的喀斯特地貌也制約著不同的水文過程。
在喀斯特地貌系統(tǒng)中,不同的水文過程作用往往使喀斯特系統(tǒng)在時(shí)空上表現(xiàn)出不同的水文地貌效應(yīng)(如以水流侵蝕為主的地貌結(jié)構(gòu)或以水流溶蝕為主的地貌結(jié)構(gòu));與此相應(yīng)的,喀斯特地貌的多樣性和差異性必然也會(huì)導(dǎo)致流域水文過程不同。一般來說,喀斯特流域地貌結(jié)構(gòu)對(duì)水體過程的作用可歸納為3類:蓄水作用、滯水作用和導(dǎo)水作用。表層溶蝕裂隙的入滲強(qiáng)度大、厚度小,往往表現(xiàn)為滯水作用;在峰林平原、峰林盆地中,土層厚度較大,具有較強(qiáng)的蓄水作用;而諸如地表、地下河及谷地等,則往往表現(xiàn)為導(dǎo)水作用。
在我國(guó),西南地區(qū)的喀斯特巖溶介質(zhì)比北方地區(qū)的具有更強(qiáng)的非均質(zhì)性和各向異性。大氣降水和地表水通過巖溶地表滲入地下轉(zhuǎn)化為地下水,對(duì)巖溶地下結(jié)構(gòu)不斷進(jìn)行改造,從而逐漸形成巖溶通道(管道、洞穴等),并有序構(gòu)成巖溶地下河系[3]。在這樣的地下河系統(tǒng)中,水流以管道流為主,巖溶管道或地下河構(gòu)成地下徑流的主干。在2個(gè)管道流之間,地下水在層理、節(jié)理和溶孔中運(yùn)動(dòng),此時(shí)的運(yùn)動(dòng)狀態(tài)與一般的滲流相似。管道流由于流動(dòng)阻力小,流速非常大。據(jù)有關(guān)資料[4]統(tǒng)計(jì),一般情況下地下河發(fā)育區(qū)的巖溶水流速枯水期為 8.64~17.28km/d,平水期為17.28~ 43.20km/d,洪水期為 43.20~ 129.60km/d,極端流速可達(dá)172.80km/d,可見在管道流或地下河發(fā)育區(qū)的地下水運(yùn)動(dòng)已不符合達(dá)西定律。此外,巖溶地下河系統(tǒng)中還包含巖溶裂隙、巖溶裂縫和巖溶孔隙等多種介質(zhì)體,形成了一個(gè)多重的復(fù)合體系,具有高度的非均質(zhì)性,使得地下河系統(tǒng)中的水流運(yùn)動(dòng)常呈現(xiàn)出達(dá)西流與非達(dá)西流并存的現(xiàn)象[5],這也為如何準(zhǔn)確描述巖溶地下水的運(yùn)動(dòng)規(guī)律增加了困難。
過去對(duì)巖溶喀斯特流域的水循環(huán)過程往往偏重于定性研究,直到 20世紀(jì)90年代以后才陸續(xù)出現(xiàn)一些針對(duì)巖溶流域水循環(huán)過程定量研究的水文模型。我國(guó)西南巖溶流域具有特殊的二元三維地質(zhì)結(jié)構(gòu)的極其復(fù)雜的地下河系統(tǒng)。在長(zhǎng)期風(fēng)化侵蝕作用下,喀斯特流域中大量的裂隙、溶溝、漏斗、天窗和落水洞使得其水文循環(huán)極其復(fù)雜。降雨時(shí),地面徑流通過這些大型管道迅速進(jìn)入地下成為地下徑流,使一次雨洪過程中地下徑流占較大比重,這與非巖溶流域的特征不同,加之多重復(fù)合體系的多種水流運(yùn)動(dòng)介質(zhì)體及達(dá)西流與非達(dá)西流并存的現(xiàn)象,給巖溶流域建立水文模型帶來相當(dāng)?shù)睦щy。
目前在巖溶喀斯特地區(qū)建立的水文模型大體可概括為以下3類:
第1類通過對(duì)在非巖溶流域上成功使用的概念性集總式流域水文模型進(jìn)行適度改造,使其在一定程度上能夠適用于巖溶地區(qū)。如袁道先[6]在水箱模型的基礎(chǔ)上開發(fā)出的桂林“丫吉模型”是最早在西南巖溶地區(qū)成功建立的流域水文模型之一。此類模型[7-12]雖然在模型構(gòu)建時(shí)對(duì)流域進(jìn)行了一定的單元?jiǎng)澐?但由于對(duì)水文過程、輸入變量及邊界條件的概化,對(duì)流域空間變異性的忽略,以及對(duì)喀斯特流域產(chǎn)匯流空間分布的不明確,使得該類模型從嚴(yán)格意義上說仍屬于黑箱模型,且在實(shí)際應(yīng)用中存在諸多問題。當(dāng)將此類集總式模型應(yīng)用于地形地貌空間變異性高度發(fā)育的喀斯特流域時(shí),會(huì)發(fā)現(xiàn)模型與客觀實(shí)際嚴(yán)重偏離,從而影響模型的計(jì)算精度和可信程度。
第2類是結(jié)合灰色理論、分形方法及人工神經(jīng)網(wǎng)絡(luò)等現(xiàn)代數(shù)學(xué)方法所建立起來的數(shù)理統(tǒng)計(jì)模型[13-14]。束龍倉(cāng)等[15]將門限自回歸模型與小波BP神經(jīng)網(wǎng)絡(luò)相結(jié)合應(yīng)用于后寨河流域出口流量的預(yù)測(cè);石朋等[16-17]基于線性矩法和MK等方法分析了貴州巖溶流域徑流過程的年度和季節(jié)變化、洪枯水變化規(guī)律及其演變的空間特征;陳宏峰等[18]利用BP人工神經(jīng)網(wǎng)絡(luò)建立了湖南洛塔河地區(qū)表層巖溶泉的徑流預(yù)測(cè)模型等。以數(shù)學(xué)方法為基礎(chǔ)的數(shù)理統(tǒng)計(jì)模型往往需要大量長(zhǎng)系列的歷史觀測(cè)資料,然而,我國(guó)西南巖溶流域處于經(jīng)濟(jì)發(fā)展相對(duì)落后的地區(qū),除少數(shù)試驗(yàn)流域近年來在國(guó)家重大研究項(xiàng)目的資助下開展了一些連續(xù)觀測(cè)外,大部分地區(qū)資料匱乏甚至無資料,從而使統(tǒng)計(jì)模型在理論上難以獲得支撐。由于此類方法缺乏必要的物理基礎(chǔ),所建模型用于模擬尚可,用于預(yù)測(cè)則存在較大的不確定性,因此在巖溶水文循環(huán)研究中很難有進(jìn)一步發(fā)展的空間。
第3類是非巖溶流域上適用的分布式/半分布式水文模型,嘗試將其應(yīng)用于巖溶流域。與集總式模型不同,分布式模型能夠同時(shí)考慮流域水文過程和植被、土壤、地形地貌等地理因素的時(shí)空變異性及其對(duì)水文過程的影響,一些固有的自然地理要素往往被作為模型要素加以考慮,因此,分布式模型能夠較好地描述流域內(nèi)細(xì)部的水文過程。任啟偉[19]在SWAT模型基礎(chǔ)上改建的刁江半分布式模型可視為分布式模型在我國(guó)西南喀斯特流域上的有益嘗試。張志才等[20]在DHSVM模型基礎(chǔ)上所建立的巖溶試驗(yàn)流域分布式模型以及蒙?;ǖ萚21]基于落水洞所構(gòu)建的一個(gè)半分布式水文模型也有一定的計(jì)算精度。在建模理論上分布式水文模型是適用于復(fù)雜條件下的西南喀斯特流域的,但也存在一些亟待解決的問題:①巖溶產(chǎn)流模式的確定。由于表層巖溶帶的存在使得巖溶流域的產(chǎn)流模式明顯不同于非巖溶流域的產(chǎn)流模式[22-24],表層巖溶帶裂隙率高,滲透性好,地表徑流只出現(xiàn)在連續(xù)降雨且雨強(qiáng)極大的降雨過程中[25]。因此,需要明確表層巖溶帶在降水入滲補(bǔ)給和產(chǎn)流中的作用[26-27]。②復(fù)雜下墊面條件的模型描述。落水洞等巖溶洼地往往是地表水快速匯集進(jìn)入地下水系的連接處[28]。分布式水文模型根據(jù)DEM采用填洼法進(jìn)行柵格匯流演算,不能將該模型直接照搬到喀斯特流域,而應(yīng)在喀斯特地貌識(shí)別基礎(chǔ)上建立地表與地下水水系的連接關(guān)系。③巖溶徑流成分及匯流路徑的確定。目前柵格型匯流主要根據(jù)地形坡度進(jìn)行計(jì)算,而產(chǎn)生于表層巖溶帶的水流不僅受地形影響,還取決于巖溶裂隙帶滲透系數(shù)的各向異性。因此,喀斯特產(chǎn)匯流受喀斯特地貌、土壤、喀斯特裂隙及地表和地下水系控制,目前尚無可直接套用或引用的模式,需要在詳細(xì)調(diào)查、巖溶水動(dòng)態(tài)加密監(jiān)測(cè)基礎(chǔ)上開展深入研究[29]。④下墊面空間異質(zhì)性導(dǎo)致的尺度效應(yīng)問題?;谛〕叨冉⒌奈锢砟P图皵?shù)值解在流域尺度上由于維數(shù)災(zāi)、參數(shù)災(zāi),增加了模型的不確定性,降低了模型的適用性。需要關(guān)注從基本尺度得到的基本規(guī)律如何應(yīng)用到流域及區(qū)域尺度上。就目前而言,分布式模型在喀斯特地區(qū)的應(yīng)用尚不多見,仍需要更多的實(shí)踐檢驗(yàn)。
針對(duì)西南喀斯特流域下墊面條件的高度空間異質(zhì)性,眾多學(xué)者開始關(guān)注不同類型下墊面條件的水文特性,從而開展了大量的室內(nèi)外試驗(yàn)。此類試驗(yàn)大體上可歸納為2類:一是對(duì)局地水文地質(zhì)參數(shù)的求證;二是對(duì)地表以下水流路徑的探測(cè)。
在水文地質(zhì)參數(shù)的求證方面,章程等[30-31]通過野外溶蝕試片法證明了不同土地利用類型下溶蝕量存在顯著差異性,同時(shí)還利用多參數(shù)自動(dòng)記錄儀探討了巖溶泉的水位、水溫等與降雨的響應(yīng)關(guān)系;張志才等[32]利用時(shí)域反射儀測(cè)定山體不同位置在不同植被覆蓋條件下的土壤含水率,分析了喀斯特峰叢山體土壤水的空間分布規(guī)律及其影響因素;劉延惠等[33]在巖溶流域設(shè)置徑流小區(qū),通過實(shí)地觀測(cè)分析不同植被類型下巖溶地表徑流的差異;王臘春等[34]通過脈沖試驗(yàn)和降雨天然脈沖分析,論述了后寨河流域巖溶水水文特征與流域貯水結(jié)構(gòu)、地質(zhì)地貌結(jié)構(gòu)的相互關(guān)系,以及裂隙流與管道流的轉(zhuǎn)化關(guān)系。魯程鵬等[35]利用電導(dǎo)率與溶質(zhì)濃度的正相關(guān)關(guān)系,在示蹤試驗(yàn)基礎(chǔ)上提出一套估算地下水流速及含水層滲透系數(shù)的方法;鄧振平等[36]采用弱酸性紅A熒光燃料類試劑作為示蹤材料,分析了巖溶落水洞和地下河出口的連通性,以及地下分水嶺的分布特征。上述試驗(yàn)研究對(duì)非均一含水介質(zhì)巖溶水運(yùn)動(dòng)研究及巖溶水文模擬具有一定的參考意義,但考慮到喀斯特流域的特殊情況,這些試驗(yàn)的結(jié)論是否具有代表性,能否在流域尺度上加以推廣值得商榷。
在地表以下水流路徑的探測(cè)方面,目前國(guó)內(nèi)尚少有此方面的研究性論文,現(xiàn)有的方法仍是以水文地質(zhì)鉆探方法為主。此類方法工程量較大,在一般的研究流域難以實(shí)施。近年來,連續(xù)電導(dǎo)率剖面儀及地質(zhì)雷達(dá)的出現(xiàn)在一定程度上提高了效率,但同時(shí)仍存在一些問題。以RH-4型連續(xù)電導(dǎo)率剖面儀為例,其在地表以下 10~500m能夠精細(xì)地反演地下情況,但自地表至地表以下10m為其探測(cè)盲區(qū),而這10m范圍恰恰是西南喀斯特流域中表層巖溶帶強(qiáng)烈發(fā)育的區(qū)段。再如地質(zhì)雷達(dá),其優(yōu)點(diǎn)是靈活便攜,且能在各種復(fù)雜地形條件下正常使用,但缺點(diǎn)是易受干擾,外界微弱的電磁信號(hào)即會(huì)影響其工作結(jié)果。此外,地表土層厚度、土壤含水量等均會(huì)影響其探測(cè)結(jié)果,而對(duì)于探測(cè)結(jié)果的反演也是見仁見智,并無嚴(yán)格統(tǒng)一的標(biāo)準(zhǔn),這樣一來,對(duì)一些細(xì)小的或由破碎巖石構(gòu)成的地下水流通道,往往難以被識(shí)別出來。
巖溶流域特殊的地質(zhì)結(jié)構(gòu)使得其下墊面條件不同于非巖溶流域。表層巖溶帶的存在使其包氣帶結(jié)構(gòu)更為復(fù)雜。那么,如何在此類特殊包氣帶條件下探討其產(chǎn)流機(jī)制和產(chǎn)流模式,將是未來巖溶流域水循環(huán)研究的一個(gè)重點(diǎn)。盡管已有的研究表明,巖溶流域雖土層淺薄,但仍屬蓄滿產(chǎn)流模式,但此種蓄滿不僅僅是土壤的蓄滿,同時(shí)需要考慮表層巖溶帶的蓄滿,而表層巖溶帶內(nèi)的水流運(yùn)移儲(chǔ)存方式又不同于土壤,那么,如何構(gòu)建出這種“土壤+表層巖溶帶”的復(fù)合蓄滿產(chǎn)流模式將是一個(gè)難點(diǎn)。
巖溶地下河系統(tǒng)中的水流運(yùn)動(dòng)常呈現(xiàn)出達(dá)西流與非達(dá)西流并存且地表河與地下河交替出現(xiàn)的情況。對(duì)于這種復(fù)雜的水流情態(tài),如何通過合理的數(shù)學(xué)物理方法加以描述將是巖溶環(huán)境下匯流過程研究的主要難點(diǎn)。在此類環(huán)境下,除非在某些重點(diǎn)區(qū)域進(jìn)行鉆井等強(qiáng)化觀測(cè),對(duì)于大部分地區(qū),由于難以確定邊界條件,傳統(tǒng)的水動(dòng)力學(xué)方法顯然不再適用,如單純使用線性水庫(kù)進(jìn)行演算,則又回歸到黑箱模型中。因此,如何在其中找到一個(gè)合適的切入點(diǎn)顯得尤為重要。目前,將分形技術(shù)用于模擬復(fù)雜的地下河系統(tǒng)與管道裂隙系統(tǒng)并概化出理想地下河系統(tǒng)響應(yīng)函數(shù)是一種有益的嘗試。
由于水文模型能夠給出流域降雨徑流的變化過程,因此,在喀斯特流域開展流域水文模型研究仍將是今后的主要研究?jī)?nèi)容之一。但水文模型的研究方法需要與研究側(cè)重點(diǎn)相適應(yīng),如研究重點(diǎn)在于考慮較大尺度上流域或區(qū)域的水資源量問題,則宜采用結(jié)構(gòu)簡(jiǎn)單又不完全失其物理意義的概念性集總式水文模型;當(dāng)以面積較小的試驗(yàn)流域?yàn)檠芯繉?duì)象,重點(diǎn)關(guān)注巖溶水循環(huán)過程的具體細(xì)節(jié)時(shí),需要構(gòu)建分布式巖溶流域水文模型,通過劃分網(wǎng)格或代表性單元,盡可能真實(shí)地反映降雨徑流的時(shí)空變化??紤]到目前的具體情況,采用一些半分布式模型并加以改進(jìn),也不失為一種折中的方法。
地貌水文學(xué)研究在非巖溶流域中一直是地貌學(xué)家和水文學(xué)家都十分關(guān)注的領(lǐng)域,而在巖溶流域中,此方面的研究更值得重點(diǎn)關(guān)注。不同的巖溶地貌類型及特殊地質(zhì)單元,諸如峰叢、峰林、天窗、落水洞等,對(duì)產(chǎn)流和匯流都會(huì)造成巨大的影響。如果能夠建立某種模式,通過地形地貌來確定不同巖溶類型下的產(chǎn)匯流特征,對(duì)于巖溶環(huán)境下的水循環(huán)研究將具有劃時(shí)代的意義。
喀斯特空間變異性及地表-地下雙重結(jié)構(gòu)對(duì)水文-生態(tài)鏈接具有重要控制作用。近年來,人類非理性的活動(dòng)對(duì)西南喀斯特山地生態(tài)水文格局的破壞十分嚴(yán)重,地表生態(tài)水文過程已經(jīng)喪失了協(xié)調(diào)發(fā)展的趨勢(shì)。因此,選擇西南喀斯特典型流域,針對(duì)地貌、植被高度空間異質(zhì)性和地表-地下雙重結(jié)構(gòu)水文地質(zhì)特征研究流域尺度的水循環(huán)動(dòng)力過程對(duì)生態(tài)環(huán)境的影響,對(duì)喀斯特流域生態(tài)系統(tǒng)的恢復(fù)和重建具有重要意義,勢(shì)必將成為未來重點(diǎn)研究?jī)?nèi)容之一。
盡管現(xiàn)有的地質(zhì)雷達(dá)、連續(xù)電導(dǎo)率剖面儀等設(shè)備應(yīng)用于巖溶流域時(shí)還存在諸多問題,但弄清巖溶地下管道的具體位置及其連接情況,對(duì)于研究巖溶地下水的運(yùn)移具有十分重要的意義。因此,在上述儀器設(shè)備基礎(chǔ)上,進(jìn)一步開發(fā)適用于巖溶地下勘探的新技術(shù)新方法顯得尤為重要。隨著3S技術(shù)的不斷進(jìn)步,將3S技術(shù)應(yīng)用到巖溶流域特殊地質(zhì)結(jié)構(gòu)的識(shí)別,可能是今后一個(gè)值得關(guān)注的發(fā)展方向。
西南巖溶流域地處我國(guó)長(zhǎng)江和珠江兩大水系的上游,具有重要的生態(tài)和環(huán)境意義,巖溶流域水循環(huán)研究對(duì)于該地區(qū)的環(huán)境保護(hù)和可持續(xù)發(fā)展具有十分重要的意義。在總結(jié)現(xiàn)有國(guó)內(nèi)外學(xué)者對(duì)我國(guó)西南巖溶流域水循環(huán)相關(guān)研究的基礎(chǔ)上,指出了目前研究中存在的一些問題,并對(duì)未來可能的研究方向提出了一些淺見。
:
[1]李陽兵,王世杰,容麗,等.西南喀斯特山地石漠化及生態(tài)恢復(fù)研究展望[J].生態(tài)學(xué)雜志,2004,23(6):84-88.
[2]梁虹,楊明德.喀斯特流域水文地貌系統(tǒng)及其識(shí)別方法探討[J].中國(guó)巖溶,1994,13(1):1-9.
[3]郭純青.中國(guó)巖溶地下河系及其水資源[J].水文地質(zhì)工程地質(zhì),2001(5):43-45.
[4]朱遠(yuǎn)峰.中國(guó)巖溶水系統(tǒng)和巖溶水資源研究進(jìn)展[C]//工程地質(zhì)、水文地質(zhì)、環(huán)境地質(zhì)論文集.北京:地震出版社,1993:194-199.
[5]郭琳,陳植華.巖溶地區(qū)地下河系統(tǒng)水資源定量評(píng)價(jià)的問題與出路[J].中國(guó)巖溶,2006,25(1):1-6.
[6]袁道先.中國(guó)南方裸露型巖溶峰叢山區(qū)巖溶水系統(tǒng)及其數(shù)學(xué)模型的研究:以桂林丫吉村為例[M].南寧:廣西師范大學(xué)出版社,1995.
[7]盧衛(wèi)中.巖溶流域水文模型初步研究[J].水文,1995(2):29-33.
[8]王臘春,史運(yùn)良.非閉合流域巖溶水模型[J].水科學(xué)進(jìn)展,1995,6(4):318-323.
[9]史運(yùn)良,王臘春.巖溶水概念分布模型[J].南京大學(xué)學(xué)報(bào):地理版,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ò)和遺傳算法的濟(jì)南市區(qū)巖溶地下水預(yù)報(bào)模型研究[J].山東農(nóng)業(yè)大學(xué)學(xué)報(bào),2004,35(3):436-441.
[15]束龍倉(cāng),董貴明,陶玉飛,地下河天窗水位變化分析及預(yù)測(cè)[J].水利學(xué)報(bào),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年來的徑流過程時(shí)空演變特征及驅(qū)動(dòng)因素探討[J].地球與環(huán)境,2009,37(3):207-213.
[18]陳宏峰,朱明秋,夏日元.湖南洛塔干河豬場(chǎng)表層巖溶泉BP人工神經(jīng)網(wǎng)絡(luò)分析[J].中國(guó)巖溶,2005,24(4):300-304.
[19]任啟偉.基于改進(jìn)SWAT模型的西南巖溶流域水量評(píng)價(jià)方法研究[D].武漢:中國(guó)地質(zhì)大學(xué),2006.
[20]張志才,陳喜,石朋,等.喀斯特流域分布式水文模型及植被生態(tài)水文效應(yīng)[J].水科學(xué)進(jìn)展,2009,20(6):806-811.
[21]蒙?;?王臘春,蘇維詞,等.基于落水洞的巖溶半分布式水文模型的構(gòu)建及其應(yīng)用[J].地理科學(xué),2009,29(4):550-554.
[22]賀中華,梁虹,黃法蘇,等.基于DEM 的喀斯特流域地貌類型的識(shí)別[J].大地測(cè)量與地球動(dòng)力學(xué),2008,28(3):46-53.
[23]蔣忠誠(chéng),袁道先.表層巖溶帶的巖溶動(dòng)力學(xué)特征及其環(huán)境和資源意義[J].地球?qū)W報(bào),1999,20(3):302-308.
[24]鄒勝章,張文慧,梁小平,等.表層巖溶帶調(diào)蓄系數(shù)定量計(jì)算:以湘西洛塔趙家灣為例[J].水文地質(zhì)工程地質(zhì),2005(4):37-42.
[25]彭韜,王世杰,張信寶,等.喀斯特坡地地表徑流系數(shù)監(jiān)測(cè)初報(bào)[J].地球與環(huán)境,2008,36(2):125-129.
[26]姜光輝,吳吉春,郭芳,等.森林覆蓋的喀斯特地區(qū)表層巖溶帶的產(chǎn)流閾值[J].水科學(xué)進(jìn)展,2008,19(1):72-77.
[27]姜光輝,郭芳.我國(guó)西南巖溶區(qū)表層巖溶帶的水文動(dòng)態(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]覃小群,蔣忠誠(chéng).表層巖溶帶及其水循環(huán)的研究進(jìn)展與發(fā)展方向[J].中國(guó)巖溶,2005,24(3):250-254.
[30]章程,謝運(yùn)球,呂勇,等.不同土地利用方式對(duì)巖溶作用的影響:以廣西弄拉峰叢洼地巖溶系統(tǒng)為例[J].地理學(xué)報(bào),2006,61(11):1181-1188.
[31]章程,袁道先,李強(qiáng),等.典型表層巖溶動(dòng)力系統(tǒng)的環(huán)境敏感性研究:以廣西弄拉峰叢洼地巖溶系統(tǒng)為例[J].水文地質(zhì)工程地質(zhì),2005(1):78-80.
[32]張志才,陳喜,石朋,等.貴州喀斯特峰叢山體土壤水分布特征及其影響因素[J].長(zhǎng)江流域資源與環(huán)境,2008,17(5):803-807.
[33]劉延惠,張喜,崔迎春,等.貴州開陽喀斯特山地幾種不同植被類型的地表徑流研究[J].貴州林業(yè)科技,2005,33(2):8-10.
[34]王臘春,許有鵬,張立峰,等.貴州普定后寨地下河流域巖溶水特征研究[J].地理科學(xué),2000,20(6):557-562.
[35]魯程鵬,束龍倉(cāng),苑利波,等.基于示蹤試驗(yàn)求解巖溶含水層水文地質(zhì)參數(shù)[J].吉林大學(xué)學(xué)報(bào):地球科學(xué)版,2009,39(4):717-721.
[36]鄧振平,周小紅,何師意,等.西南巖溶石山地區(qū)巖溶地下水示蹤試驗(yàn)與分析[J].中國(guó)巖溶,2007,26(2):163-169.