趙詠梅，趙 夢(mèng)，焦樹林，李 靜，曾友華，趙宗權(quán)

(1.畢節(jié)市人民政府防汛抗旱指揮部辦公室,貴州 畢節(jié) 551700；2.貴州師范大學(xué) 地理與環(huán)境科學(xué)學(xué)院, 貴州 貴陽(yáng) 550001)

1 引言

水是生命之源，是人類生存和維持良好自然生態(tài)環(huán)境的自然資源，對(duì)維持社會(huì)經(jīng)濟(jì)發(fā)展具有不可替代的作用[1，2]。水資源在很大程度上已經(jīng)成為制約中國(guó)經(jīng)濟(jì)轉(zhuǎn)型以及增長(zhǎng)的關(guān)鍵要素，在過(guò)去40年的高速經(jīng)濟(jì)增長(zhǎng)過(guò)程中，水資源污染和浪費(fèi)等現(xiàn)象極為嚴(yán)重，90%的地下水遭到不同程度的污染，各類污染事件不斷發(fā)生[3]。所以在國(guó)民經(jīng)濟(jì)發(fā)展和社會(huì)進(jìn)步中對(duì)水資源的治理與協(xié)調(diào)至關(guān)重要。

據(jù)統(tǒng)計(jì)：畢節(jié)市多年平均水資源總量約為1.12×1010m3，根據(jù)2014年畢節(jié)市水資源公報(bào)，畢節(jié)市地下水資源量為4.09×109m3，其中長(zhǎng)江流域?yàn)?.93×109m3，珠江流域?yàn)?.65×108m3，全市大于1×104m的河流有193條，2014年畢節(jié)市地表水資源總量為1.38×1010m3，產(chǎn)水規(guī)模為5.14×105m3/km2[4]。而就貴州省畢節(jié)市關(guān)于水資源的研究大多集中在對(duì)暴雨的分析[5～9]、水資源生態(tài)承載力探究[10]、水資源保護(hù)與評(píng)價(jià)[11，12]、降水分析[13]，以及地下水水文地球化學(xué)特征[14]等方面，但在不同頻率條件下的降水量的研究目前還未有報(bào)道。

鑒于此，本文以畢節(jié)市為研究對(duì)象(圖1)，通過(guò)降水經(jīng)驗(yàn)頻率計(jì)算得到畢節(jié)氣象站近50年(1961～2015年)降水的經(jīng)驗(yàn)頻率曲線，選用皮爾遜Ⅲ型理論曲線進(jìn)行多年不同頻率下降水的擬合，得出一條最適合畢節(jié)站的降水頻率理論擬合曲線。通過(guò)這條擬合曲線推導(dǎo)畢節(jié)地區(qū)不同頻率條件下的降水量，以期為該地區(qū)水資源的更好利用提供理論依據(jù)，為相關(guān)水利部門管理水資源提供指導(dǎo)及類似的研究區(qū)水資源管理提供參考和借鑒。

圖1 研究區(qū)分布

2 材料與方法

2.1 研究區(qū)概況

畢節(jié)市為貴州省下轄地級(jí)市，位于貴州西北部，貴州金三角地區(qū)之一，屬烏蒙山腹地，川、滇、黔之交界帶，長(zhǎng)江珠江之屏障，西鄰云南，北接四川，是烏江、北盤江、赤水河發(fā)源地，是一個(gè)多民族聚居、歷史文化燦爛、資源富集的地方。國(guó)家“西電東送”的重要能源基地，國(guó)家新型能源化工基地，國(guó)家新能源汽車高新技術(shù)產(chǎn)業(yè)化基地，國(guó)家生物醫(yī)藥產(chǎn)業(yè)基地，現(xiàn)代山地高效生態(tài)農(nóng)業(yè)、新能源、新型建材、以大數(shù)據(jù)為核心的服務(wù)外包和呼叫中心等多種新興產(chǎn)業(yè)的集聚地，全國(guó)唯一以“開發(fā)扶貧、生態(tài)建設(shè)”為主題的試驗(yàn)區(qū)。川、滇、黔、渝結(jié)合部區(qū)域性中心城市，西南地區(qū)區(qū)域性重要綜合交通樞紐，珠三角連接西南地區(qū)、長(zhǎng)三角連接?xùn)|盟地區(qū)的重要通道。礦產(chǎn)資源儲(chǔ)量巨大，已探明礦產(chǎn)60多種，磷礦儲(chǔ)量名列全國(guó)前茅，煤炭遠(yuǎn)景儲(chǔ)量居于長(zhǎng)江以南之冠，鐵礦、銅礦、鉛礦、稀土礦儲(chǔ)量處在貴州第一；生物資源多樣，動(dòng)植物資源2800多種，有全國(guó)、全省如馬鈴薯之鄉(xiāng)、白蒜之鄉(xiāng)等眾多“地理標(biāo)志”；水能資源豐富，河湖水系縱橫交錯(cuò)，流域面積100 km2以上的河流有80條，是古夜郎政治經(jīng)濟(jì)文化中心之一，中國(guó)南方古人類文化發(fā)祥地，景色秀麗，被譽(yù)為“洞天湖地、花海鶴鄉(xiāng)、避暑天堂”；氣候涼爽，是避暑旅游城市觀測(cè)點(diǎn)，也是三省紅都，為長(zhǎng)江以南最后一塊革命根據(jù)地，厚重的歷史文化在全國(guó)具有唯一性。

2.1.1 地形地貌

畢節(jié)市位于貴州省西北部，地處烏蒙山區(qū)。從地理區(qū)位上看，東接貴陽(yáng)市、遵義市，南與安順市、六盤水市相連，西與云南昭通市、曲靖市相鄰，北與四川瀘州市接壤。從地形上來(lái)看，該區(qū)位于云貴高原的東部斜坡上，地勢(shì)西部高東部低，該區(qū)域以山地為主，境內(nèi)河流眾多，而且?guī)r溶中多發(fā)育有地下暗河，畢節(jié)市地質(zhì)構(gòu)造復(fù)雜，褶皺層與斷裂帶較多，屬于典型的喀斯特巖溶地貌。主要的山脈有大婁山、烏蒙山、老王山，分別位于畢節(jié)市的北部、西部、西南部。巖溶地貌類型多種多樣，東部多為石峰林、谷地、峰叢、緩丘和洼地，中部多為峰叢、槽谷以及丘陵洼地，西部主要是高原、巖溶、緩丘和盆地。從地勢(shì)上來(lái)看，畢節(jié)市自西向東呈三級(jí)階梯狀，第一階梯為高原地帶，平均海拔為2000～2400 m，包括威寧縣、赫章縣西部；第二階梯為山地，平均海拔為1400～1800 m，包括赫章縣東部、七星關(guān)區(qū)、大方縣、黔西縣、納雍縣以及織金縣西部；第三級(jí)階梯為低山丘陵帶，平均海拔為1000～1400 m，包括金沙和黔西兩縣及織金縣東部[4]。

2.1.2 河流水文

烏蒙山位于赫章縣城西部和威寧縣境內(nèi)，其山系為白水河、牛欄江、烏江和北盤江的分水嶺；而大婁山脈西起赫章縣東部，向東經(jīng)過(guò)七星關(guān)區(qū)、大方縣北部，之后進(jìn)人金沙縣，一直延伸到遵義市境內(nèi)，其山脈在畢節(jié)市境內(nèi)為烏江水系和長(zhǎng)江水系的分水嶺；老王山山脈呈西北-東南走向，西北處與烏蒙山東支相連接，向東延伸到六盤水市的水城、六枝一帶，為烏江上游三岔河與珠江水系北盤江的分水嶺[4](表1)。

表1 畢節(jié)市水系河流

數(shù)據(jù)來(lái)源：2014年畢節(jié)市水資源公報(bào)

2.1.3 社會(huì)經(jīng)濟(jì)發(fā)展

根據(jù)統(tǒng)計(jì)年鑒可知：國(guó)內(nèi)生產(chǎn)總值方面，七星關(guān)區(qū)和金沙縣約占畢節(jié)市生產(chǎn)總值的48%。七星關(guān)區(qū)經(jīng)濟(jì)發(fā)展水平處各縣的領(lǐng)先地位，常住人口處于全市領(lǐng)先地位。威寧縣土地面積最大并且常住人口最多，經(jīng)濟(jì)發(fā)展水平靠前；從常住人口來(lái)看，大多集中在威寧縣、七星關(guān)區(qū)和織金縣；從土地面積來(lái)看，大多集中在威寧縣、大方縣以及七星關(guān)區(qū)(表2)。

表2 2015年畢節(jié)市各縣經(jīng)濟(jì)狀況

資料來(lái)源：2016年畢節(jié)市統(tǒng)計(jì)年鑒

綜上可知，畢節(jié)市的社會(huì)經(jīng)濟(jì)發(fā)展并不平衡，地區(qū)間發(fā)展存在明顯的差異。但畢節(jié)市已形成了以七星關(guān)區(qū)為核心的發(fā)展格局，城鎮(zhèn)人口大部分集中在七星關(guān)區(qū)和威寧縣。人們趨向經(jīng)濟(jì)更為發(fā)達(dá)的中心區(qū)域集中居住，而土地面積的大小對(duì)于常住人口的影響并不大。

2.2 數(shù)據(jù)處理與分析

本文選取1961～2015畢節(jié)水文站年降水量統(tǒng)計(jì)數(shù)據(jù)，通過(guò)降水經(jīng)驗(yàn)頻率計(jì)算得到畢節(jié)站近50年(1961～2015)降水的經(jīng)驗(yàn)頻率曲線(圖2、3)，利用spss軟件做出畢節(jié)站1961～2015年降水趨勢(shì)圖(圖2)。

2.2.1 畢節(jié)站降水變化特征

對(duì)畢節(jié)站近50年(1961～2015)年降水量進(jìn)行統(tǒng)計(jì)并利用spss軟件分析得出，年降水量的線性趨勢(shì)方程為：y=3711.7-1.4x，年降水量總體上呈鋸齒狀下降的態(tài)勢(shì)，1961～2015年的平均降水量為888.37 mm，年降水量最小值出現(xiàn)在1989年，為615.0 mm，年降水量最多的年份為1983年，為1285.44 mm。

圖2 畢節(jié)站降水變化趨勢(shì)

2.2.2 頻率計(jì)算原理

2.2.2.1 隨機(jī)變量

實(shí)踐表明，偶然現(xiàn)象雖有它不確定的一面，但觀測(cè)大量同類的偶然現(xiàn)象，還是可以發(fā)現(xiàn)它的統(tǒng)計(jì)規(guī)律。由于這種規(guī)律是和它的出現(xiàn)機(jī)會(huì)分不開，因此在數(shù)學(xué)上稱這種偶然現(xiàn)象為隨機(jī)現(xiàn)象。隨機(jī)現(xiàn)象可以取得不同的數(shù)值，這些數(shù)值可以用一個(gè)變量來(lái)表示，這種變量就叫做隨機(jī)變量。隨機(jī)變量的各個(gè)數(shù)值都對(duì)應(yīng)一定的出現(xiàn)機(jī)會(huì)。水文計(jì)算中的流量和雨量等水文資料都是隨機(jī)變量，如畢節(jié)氣象站1999年的年降雨量792.9 mm就是一個(gè)隨機(jī)變量。

2.2.2.2 總體和樣本

數(shù)理統(tǒng)計(jì)中，把隨機(jī)變量的全體，亦即包括整體的系列，稱為總體。從總體中任意抽取一部分系列，稱為樣本。許多水文現(xiàn)象的總體是無(wú)限的，實(shí)際上也無(wú)法取得。例如畢節(jié)氣象站的年降雨量，其總體應(yīng)包括從古至今以至延長(zhǎng)到未來(lái)無(wú)限年代的所有資料。而現(xiàn)在所觀測(cè)到的幾十年的年降雨量資料，只不過(guò)是總體中的一部分，是一個(gè)很有限的樣本。

2.2.2.3 頻率

在數(shù)理統(tǒng)計(jì)中，頻率是指隨即變量某值在樣本中的出現(xiàn)機(jī)會(huì)(即概率或機(jī)率)。在水文計(jì)算中，頻率則表示等于或大于某一變量可能出現(xiàn)的機(jī)會(huì)；具體地說(shuō)，等于或大于某一變量的出現(xiàn)次數(shù)與總次數(shù)的比之，就稱為這個(gè)變量的頻率。頻率通常以百分?jǐn)?shù)表示。例如用畢節(jié)氣象站55年的實(shí)測(cè)降水量資料，將這55年的資料按大小順序排列，則等于或大于第一個(gè)最大雨量值在55年中只有一次，等于或大于第二大的雨量值在55年中有兩次，依次類推，則等于或大于最后一個(gè)最小的雨量值在55年中共有55次。將這些等于或大于的次數(shù)分別除以總次數(shù)，所得的百分?jǐn)?shù)便是各個(gè)年雨量值所對(duì)應(yīng)的頻率，如表3所示。頻率常用符號(hào)P來(lái)表示，P=5%表示等于或大于某變量的出現(xiàn)機(jī)會(huì)為5%。

表3 畢節(jié)站降水頻率計(jì)算

續(xù)表3

年份年降水量序號(hào)按大小排序模比系數(shù)kiki-1(ki-1)2(ki-1)3pCvCs1998839.91 38842.10 0.95 -0.05 0.00 0.00 67.861999852.69 39841.03 0.95 -0.05 0.00 0.00 69.6420001067.26 40839.91 0.95 -0.05 0.00 0.00 71.432001920.44 41829.09 0.93 -0.07 0.00 0.00 73.212002688.71 42801.09 0.90 -0.10 0.01 0.00 75.002003767.44 43783.08 0.88 -0.12 0.01 0.00 76.792004854.43 44781.49 0.88 -0.12 0.01 0.00 78.572005767.04 45767.44 0.86 -0.14 0.02 0.00 80.362006960.65 46767.04 0.86 -0.14 0.02 0.00 82.142007895.72 47742.14 0.84 -0.16 0.03 0.00 83.9320081114.86 48735.40 0.83 -0.17 0.03 -0.01 85.712009699.16 49720.33 0.81 -0.19 0.04 -0.01 87.502010742.14 50699.16 0.79 -0.21 0.05 -0.01 89.292011665.84 51688.71 0.78 -0.22 0.05 -0.01 91.072012829.09 52665.84 0.75 -0.25 0.06 -0.02 92.862013783.08 53655.27 0.74 -0.26 0.07 -0.02 94.6420141014.28 54645.08 0.73 -0.27 0.08 -0.02 96.4320151040.93 55615.50 0.69 -0.31 0.09 -0.03 98.21平均值888.37

注：以上數(shù)據(jù)來(lái)自于貴州省統(tǒng)計(jì)氣象統(tǒng)計(jì)年鑒

2.2.2.4 重現(xiàn)期

頻率是一個(gè)抽象的統(tǒng)計(jì)術(shù)語(yǔ)，不便于理解。因此，水文計(jì)算中經(jīng)常用重現(xiàn)期來(lái)表示隨機(jī)變量的出現(xiàn)機(jī)會(huì)。所謂重現(xiàn)期是指在許多的實(shí)驗(yàn)里，某一事件重復(fù)出現(xiàn)的時(shí)間間隔的平均數(shù)，即平均的重現(xiàn)間隔期。重現(xiàn)期常用符號(hào)T來(lái)表示。頻率與重現(xiàn)期的關(guān)系有兩種表示法：

當(dāng)研究暴雨洪水問(wèn)題時(shí)，一般p<50%，采用

(1)

式(1)中：T表示重現(xiàn)期以年計(jì)；P表示頻率以小數(shù)或百分?jǐn)?shù)計(jì)。例如P=1%，即平均 100 年可以出現(xiàn)1次，亦即重現(xiàn)期為100年，或稱100年一遇。

當(dāng)研究枯水問(wèn)題時(shí)，一般p>50%，采用：

(2)

2.2.3 經(jīng)驗(yàn)頻率計(jì)算

根據(jù)實(shí)測(cè)水文資料，按從大到小的次序進(jìn)行排列，然后用經(jīng)驗(yàn)頻率公式計(jì)算出來(lái)的頻率，成為經(jīng)驗(yàn)頻率。以水文變量為縱坐標(biāo)，以經(jīng)驗(yàn)頻率為橫坐標(biāo)，點(diǎn)繪經(jīng)驗(yàn)頻率點(diǎn)據(jù)，根據(jù)點(diǎn)群分布趨勢(shì)，繪出一條平滑的曲線，稱為頻率曲線。經(jīng)驗(yàn)頻率曲線是根據(jù)實(shí)測(cè)資料(樣本)直接繪成的。

2.2.3.1 經(jīng)驗(yàn)頻率公式

計(jì)算經(jīng)驗(yàn)頻率的公式有好幾種，目前我國(guó)常用計(jì)算經(jīng)驗(yàn)頻率公式為：

(3)

式(3)中：p等于或大于某變量的經(jīng)驗(yàn)頻率；m為系列按遞減次序排列時(shí)，各隨即變量的順序號(hào)；n表示觀測(cè)資料的總項(xiàng)數(shù)。

2.2.3.2 經(jīng)驗(yàn)頻率曲線

用畢節(jié)站1961～2015年實(shí)測(cè)年的降雨量資料，按表格式計(jì)算其經(jīng)驗(yàn)頻率。表3中第(9)欄經(jīng)驗(yàn)頻率是按式(3)計(jì)算得出，以表3中第(4)欄的年降雨量為縱坐標(biāo)，第(9)欄相應(yīng)的經(jīng)驗(yàn)頻率P為橫坐標(biāo)，在等分格紙上點(diǎn)繪經(jīng)驗(yàn)頻率點(diǎn)據(jù)，然后通過(guò)點(diǎn)群中間目估繪出一條平滑曲線，即為該站多年降水量經(jīng)驗(yàn)頻率曲線(圖3)。

圖3 畢節(jié)氣象站1961～2015降水經(jīng)驗(yàn)頻率分布

2.2.3.3 統(tǒng)計(jì)參數(shù)

水文系列同各種隨機(jī)變量的分布情況是不相同的。有的變量數(shù)值大，有的變量數(shù)值??；有的變量比較集中，有的則比較分散；有的變量呈對(duì)稱分布，有的則不對(duì)稱，等等。在數(shù)理統(tǒng)計(jì)中，水文學(xué)中常用的統(tǒng)計(jì)參數(shù)為均值，離差系數(shù)和偏態(tài)系數(shù)等。

均值，即算術(shù)平均數(shù)，用表示，計(jì)算如下：

(4)

離差系數(shù)：

(5)

偏態(tài)系數(shù)也叫做偏差系數(shù)，以符號(hào)Cs表示。它是反映變量在均值兩邊是否對(duì)稱或不對(duì)稱程度的一個(gè)系數(shù)。

(6)

2.2.4 配線法

(7)

表4 頻率曲線選配計(jì)算

續(xù)表4

頻率(p)第一次配線第二次配線X=888.37X=888.37Cv=0.15Cv=0.15Cs=0.5CvCs=1CvKpxpkpXp30.000 1.08957.111.08955.7640.000 1.04920.671.03918.9850.000 1.00886.821.00885.0460.000 0.96853.160.96851.5370.000 0.92817.380.92816.1575.000 0.90797.660.90796.7780.000 0.87775.790.87775.3685.000 0.84750.40.84750.6390.000 0.81718.630.81719.8895.000 0.76671.870.76675.0197.000 0.72641.720.73646.3399.000 0.66585.250.67593.1499.900 0.55489.830.57504.8799.990 0.46412.620.49434.9699.999 0.39346.550.42376.24頻率(p)第三次配線第四次配線X=888.37X=888.37Cv=0.15Cv=0.15Cs=1.5CvCs=2CvKpxpkpXp0.001 1.74 1544.581.771574.630.010 1.63 1449.111.661471.220.033 1.57 1395.371.591413.430.050 1.55 1376.461.571393.170.100 1.51 1343.181.531357.610.200 1.47 1308.421.491320.610.333 1.44 1281.691.451292.260.500 1.42 1259.691.4312691.000 1.37 1220.221.381227.42.000 1.33 1177.841.331182.973.333 1.29 1144.241.291147.925.000 1.26 1115.741.261118.3210.000 1.20 1062.031.201062.8520.000 1.12 998.821.12998.1425.000 1.10 975.321.10974.2630.000 1.07 954.461.07953.1340.000 1.03 917.371.03915.7650.000 0.99 883.380.99881.7260.000 0.96 850.020.96848.5370.000 0.92 815.040.92813.9675.000 0.90 795.980.90795.2380.000 0.87 775.010.87774.7285.000 0.85 750.920.85751.2790.000 0.81 721.130.81722.4595.000 0.76 678.030.77681.1297.000 0.73 650.720.74655.1699.000 0.68 600.560.68608.0199.900 0.58 518.850.60532.6999.990 0.51 455.550.54475.7799.999 0.45 403.430.48429.94

3 結(jié)果分析

圖4 畢節(jié)氣象站1961-2015年降水經(jīng)驗(yàn)頻率分布

4 結(jié)語(yǔ)

開發(fā)和改造河流，興建和管理農(nóng)田水利工程都必須掌握河流長(zhǎng)期水文形式和近期的確切水情。只有通過(guò)通過(guò)水文預(yù)報(bào)是預(yù)測(cè)面臨期間的為工程施工、運(yùn)營(yíng)、管理提供科學(xué)依據(jù)，而水利計(jì)算成果關(guān)系到工程的規(guī)模和工程效益，必須加以重視。