劉慧 劉德安

摘要：利用四川省宜賓市的1953—2010年的春季降水量數(shù)據(jù)資料，分別應(yīng)用正態(tài)分布、Gamma分布、對(duì)數(shù)正態(tài)分布和P-Ⅲ型分布對(duì)該地區(qū)的降水分布特征進(jìn)行研究。在此基礎(chǔ)上，利用離差平方和最小準(zhǔn)則進(jìn)行4種分布函數(shù)的擬合優(yōu)度檢驗(yàn)。結(jié)果表明，4種分布函數(shù)均能擬合宜賓市春季降水分布;通過(guò)擬合優(yōu)度檢驗(yàn)表明，對(duì)數(shù)正態(tài)分布對(duì)宜賓市春季降水量的擬合效果最好，Gamma分布次之，P-Ⅲ型分布和正態(tài)分布的擬合效果相對(duì)較差。

關(guān)鍵詞：降水量擬合;擬合檢驗(yàn);對(duì)比分析

中圖分類號(hào)：P468? ? ? ? ?文獻(xiàn)標(biāo)識(shí)碼：A

文章編號(hào)：0439-8114（2019）07-0047-05

Abstract： By utilizing spring precipitation data of Yibin city in Sichuan province from 1953 to 2010， normal distribution， Gamma distribution， logarithmic normal distribution， P-Ⅲ distribution were applied to study the region's precipitation distribution characteristics. Based on this， the minimum of sum square variation method was used for goodness of fit test of the four distribution functions. The result shows that the four distribution functions can all fit with precipitation distribution of spring in Yibin city. The goodness of fit tests show that logarithmic normal distribution is best fit with the precipitation distribution of spring in Yibin city， Gamma distribution is the second best. The fitting effect of P-III distribution and normal distribution is relatively poor.

Key words： precipitation totals fitness; fitness test; comparative analysis

降水是地面觀測(cè)中最基本的氣象要素之一，也是決定某地氣候特征的關(guān)鍵因子。降水不僅直接影響土壤水分平衡，也間接影響太陽(yáng)輻射、氣溫、濕度等天氣變量，同時(shí)降雨量也是農(nóng)業(yè)生產(chǎn)和農(nóng)業(yè)生態(tài)系統(tǒng)的重要影響因子。而某地春季的降水量更是對(duì)農(nóng)作物的春耕有著不可小視的作用，關(guān)乎到作物的初生長(zhǎng)。

研究降水量及其分布規(guī)律不僅有重要的氣候?qū)W意義，而且對(duì)農(nóng)業(yè)生產(chǎn)、水資源管理與應(yīng)用以及防災(zāi)、減災(zāi)有重要的現(xiàn)實(shí)意義。因此，對(duì)春季降水量概率分布的研究是很有必要的，這就需要一個(gè)能夠較好擬合本地春季降水量的概率分布函數(shù)，以此對(duì)本地的春季降水有所把握和合理利用。同時(shí)，在氣候統(tǒng)計(jì)診斷和預(yù)報(bào)中，經(jīng)常在正態(tài)假定的前提下進(jìn)行研究工作，但如果氣候變量不符合正態(tài)假定的前提，則將會(huì)給預(yù)報(bào)帶來(lái)一定的誤差。從這個(gè)方面來(lái)講，對(duì)各種氣象要素所服從的各種非正態(tài)分布的研究有著重要意義。

關(guān)于降水量概率分布的研究，目前國(guó)內(nèi)外已經(jīng)取得了很多研究成果。在國(guó)內(nèi)，科學(xué)工作者對(duì)各氣象要素所服從分布情況開展了大量的研究工作。在正態(tài)分布模型的研究中，胡文東等[1]指出寧夏大多數(shù)氣象站的年降水量服從正態(tài)分布，而季度降水量及逐月降水量服從正態(tài)分布的氣象站相對(duì)較少;周雁翎等[2]在對(duì)粵西降水量分布的研究中指出，除了各站夏季降水量均服從正態(tài)分布外，粵西各站月、春、秋、冬和年降水量大都不服從正態(tài)分布。對(duì)于多雨地區(qū)，情況則不同，吳慧[3]在對(duì)海南省降水量正態(tài)分布特征的研究中指出，在0.01的信度檢驗(yàn)下，海南省各站春、夏、秋季降水量基本符合正態(tài)分布，冬季僅北部和東部地區(qū)符合正態(tài)分布，月降水量的正態(tài)性也隨季節(jié)而變化且具有局地性。由于大多的天氣氣候極值往往出現(xiàn)于非正態(tài)時(shí)間序列中，僅用正態(tài)分布去模擬降水序列可能產(chǎn)生較大誤差。從這個(gè)意義上也反映出，認(rèn)識(shí)各種氣象要素多服從的各種非正態(tài)分布很有意義。

各國(guó)科學(xué)工作者都致力于非正態(tài)性分布的研究。如Entekhabi等[4]提出了用Gamma分布描述土壤水分和徑流。自2002年以來(lái)，有學(xué)者先后利用 Gamma分布，WEIBULL分布和負(fù)指數(shù)分布擬合了降水量、降水日數(shù)、徑流和蒸發(fā)等各種地表參量的地理空間分布函數(shù)，取得了較好的效果[5-7]。段莉珠等[8]根據(jù)潞西氣象站實(shí)測(cè)降水量資料，驗(yàn)證了該站降水量為P-Ⅲ型分布;吳樟榮[9]得出P-Ⅲ型分布能很好地?cái)M合廣東省登陸臺(tái)風(fēng)的年最大風(fēng)速資料;周建康等[10]通過(guò)對(duì)南京站降水量的統(tǒng)計(jì)分析，驗(yàn)證了南京站降水量雨型為P-Ⅲ型分布。丁裕國(guó)[11]曾用Γ分布模式研究過(guò)北京市等5個(gè)站的月降水量的統(tǒng)計(jì)特征;在國(guó)外，Mooley[12]曾用Γ概率分布模式對(duì)亞洲夏季月降水量進(jìn)行模擬。李偉光等[13]對(duì)海南省?？谑胁煌瑫r(shí)間尺度降水量進(jìn)行了分布擬合，指出Gamma分布函數(shù)對(duì)?？谑胁煌瑫r(shí)間尺度降水量均有較好擬合效果;石強(qiáng)等[14]用Gamma分布擬合青島市各月和年降水量，并建立了適用于青島市降水量擬合的二參數(shù)Gamma分布。由此看出，非正態(tài)分布所涉獵的研究范圍十分廣泛，因此非正態(tài)分布的研究對(duì)做好氣候的統(tǒng)計(jì)診斷和預(yù)報(bào)具有十分重要的意義。

根據(jù)前人的研究成果，丁裕國(guó)等[15]指出，正態(tài)分布常用于擬合各地月平均氣溫、年平均氣溫、部分地區(qū)的日（時(shí)）平均氣溫、多雨地區(qū)的年降水量;各種偏態(tài)分布常用于擬合地區(qū)逐日（時(shí)）平均或最高（低）氣溫、日（月）降水量、無(wú)雨期長(zhǎng)度、相對(duì)濕度、能見度、短歷時(shí)降水強(qiáng)度、徑流、土壤水分等。由此可以看出，正態(tài)分布和各種非正態(tài)分布都能夠?qū)σ欢ǖ臍庀蟋F(xiàn)象進(jìn)行擬合，且均具有十分重要的意義。本研究分別利用正態(tài)分布以及P-Ⅲ型分布、Gamma分布、對(duì)數(shù)正態(tài)分布3種偏態(tài)分布對(duì)四川省宜賓市春季（3—5月）的降水量資料進(jìn)行了擬合，并做出了對(duì)比分析，旨在尋求對(duì)四川省宜賓市春季降水量擬合更好的分布函數(shù)，以期進(jìn)一步提高本地氣候統(tǒng)計(jì)診斷和預(yù)報(bào)的準(zhǔn)確度。

1? 資料與方法

1.1? 資料來(lái)源與研究區(qū)概況

數(shù)據(jù)資料取自四川省宜賓市的春季降水量資料，為1953—2010年的春季（3—5月）累積降水量數(shù)據(jù)。宜賓市位于四川省南部，處于川、滇、黔三省結(jié)合部的金沙江、岷江、長(zhǎng)江匯流地帶。地跨北緯27°50′—29°16′、東經(jīng)103°36′—105°20′。宜賓市屬中亞熱帶濕潤(rùn)季風(fēng)氣候，低丘、河谷兼有南亞熱帶的氣候?qū)傩?具有氣候溫和、熱量豐足、雨量充沛、光照適宜、無(wú)霜期長(zhǎng)、冬暖春早、四季分明的特點(diǎn);年平均氣溫18 ℃左右;年平均降水量1 050～1 618 mm，5—10月為雨季，降水量占全年的81.7%，主汛期為7—9月，降雨量更集中，占全年總降雨量的51%;年平均日照時(shí)間為1 000～1 130 h，無(wú)霜期為334～360 d;年平均風(fēng)速僅為1.23 m/s，多為西北風(fēng)和東北風(fēng)，靜風(fēng)頻率較大，高達(dá)34%～53%，風(fēng)速小。

1.2? 研究方法

1.2.1? 降水概率分布模型及其參數(shù)估計(jì)? 首先假設(shè)降水量概率分布分別服從正態(tài)分布、Gamma分布、對(duì)數(shù)正態(tài)分布和P-Ⅲ型分布。在進(jìn)行降水量的概率分布擬合時(shí)，對(duì)正態(tài)分布的參數(shù)估計(jì)方法采用矩法估計(jì)（根據(jù)樣本矩直接估計(jì)兩個(gè)參數(shù)），對(duì)于非正態(tài)分布的Gamma分布、對(duì)數(shù)正態(tài)分布和P-Ⅲ型分布選用極大似然估計(jì)的方法確定函數(shù)的參數(shù)。極大似然估計(jì)方法是參數(shù)估計(jì)中最重要、也是應(yīng)用最廣泛的方法之一，也就是選取使得似然函數(shù)達(dá)到最大的數(shù)值作為參數(shù)的估計(jì)值。

1）正態(tài)分布及其參數(shù)估計(jì)。正態(tài)分布作為最常見且應(yīng)用最廣泛的隨機(jī)變量的分布，在隨機(jī)變量的分布中居于中心地位。其在概率論與數(shù)理統(tǒng)計(jì)中無(wú)論在理論研究還是實(shí)際應(yīng)用上都占有極其重要的地位，并且數(shù)理統(tǒng)計(jì)中許多重要問(wèn)題的解決都是以正態(tài)分布為基礎(chǔ)的。在氣候統(tǒng)計(jì)診斷和預(yù)測(cè)中，大多數(shù)診斷方法和預(yù)測(cè)模型都是建立在假設(shè)氣候變量是正態(tài)分布的前提下進(jìn)行的。因此，對(duì)于氣候變量是否為正態(tài)分布的檢驗(yàn)不僅很有必要而且很重要。以往對(duì)氣象諸要素的統(tǒng)計(jì)特性研究表明，月平均氣溫、多雨地區(qū)的年降水量一般近似服從正態(tài)分布。本研究亦采用正態(tài)分布進(jìn)行宜賓市春季降水量的分布擬合，其概率密度函數(shù)為：

1.2.2? 擬合假設(shè)檢驗(yàn)方法? 在降水量概率分布擬合之后，還要進(jìn)行擬合優(yōu)度的檢驗(yàn)。本研究采用K-S檢驗(yàn)法來(lái)完成。K-S檢驗(yàn)全稱是Kolmogorov-Smirnov檢驗(yàn)，亦稱擬合優(yōu)度檢驗(yàn)法。它是用來(lái)檢驗(yàn)給定的一組樣本數(shù)據(jù)的分布是否符合某一特定分布的方法，是檢驗(yàn)小樣本擬合效果的常用方法。它的檢驗(yàn)方法是以樣本數(shù)據(jù)的累計(jì)頻數(shù)分布與特定理論分布比較，若兩者間的差距很小，則推論該樣本取自某特定分布族。

以Gamma函數(shù)為例，K-S檢驗(yàn)進(jìn)行擬合效果的檢驗(yàn)方法為設(shè)x1，x2，x3，…，xn為來(lái)自一個(gè)總體具有分布函數(shù)Γ（x）的隨機(jī)樣本，且設(shè)該樣本是按從小到大的順序排列的有序樣本，于是可以應(yīng)用有序樣本求得Γ（x）的置信界限值，有序樣本的分布函數(shù)[15]：

1.2.3? 對(duì)比分析方法? 對(duì)以上4種降水量概率分布函數(shù)的擬合效果進(jìn)行對(duì)比分析，比較哪種函數(shù)對(duì)宜賓市的春季降水量有更好的擬合效果。擬合優(yōu)度評(píng)價(jià)是選擇概率分布函數(shù)的重要指標(biāo)，比較常用的擬合優(yōu)度評(píng)價(jià)方法有直觀的圖形分析法、Genest-Rivest方法、AIC信息準(zhǔn)則、離差平方和最小準(zhǔn)則（OLS）。本研究采用離差平方和最小準(zhǔn)則法來(lái)對(duì)函數(shù)進(jìn)行定量的擬合優(yōu)度檢驗(yàn)，以確定最優(yōu)的分布函數(shù)。OLS的計(jì)算公式為：

2? 應(yīng)用

以四川省宜賓市為例，采用宜賓市的春季降水量數(shù)據(jù)資料，進(jìn)行春季降水量的擬合、擬合檢驗(yàn)和對(duì)比分析，以期得到最優(yōu)的擬合宜賓市春季降水量的分布函數(shù)。

2.1? 各概率分布函數(shù)參數(shù)、概率值的確定

根據(jù)前面的討論，首先假設(shè)降雨量是服從正態(tài)分布、Gamma分布、對(duì)數(shù)正態(tài)分布和P-Ⅲ型分布的。其中，正態(tài)分布的參數(shù)可以按照公式（2）計(jì)算，Gamma分布、對(duì)數(shù)正態(tài)分布和P-Ⅲ型分布的參數(shù)采用極大似然估計(jì)法求得，可分別按照公式（4）、（6）、（8）進(jìn)行估計(jì)。各概率分布函數(shù)參數(shù)的計(jì)算結(jié)果如表1所示。

將表1中相應(yīng)的參數(shù)分別代入公式（1）、（3）、（5）、（7），即可求得各自的概率密度。從所得的計(jì)算結(jié)果可以直觀地看到各降水量區(qū)間出現(xiàn)的可能性，為減災(zāi)防災(zāi)工程提供相應(yīng)的科學(xué)依據(jù)。當(dāng)降水量值大于100 mm時(shí)，通過(guò)各分布函數(shù)計(jì)算的概率值得到其出現(xiàn)的最大概率僅為5%左右，說(shuō)明宜賓市春季降水量超過(guò)100 mm的概率不是很大;而其春季降水量低于30 mm出現(xiàn)的概率也較低，從各分布計(jì)算的概率值來(lái)看，最高也只是5%左右。

2.2? 擬合檢驗(yàn)

根據(jù)宜賓市春季降水量數(shù)據(jù)，利用K-S檢驗(yàn)方法計(jì)算的結(jié)果如表2所示。查表可知，在n=58， α=0.05時(shí)的K-S檢驗(yàn)的檢驗(yàn)統(tǒng)計(jì)量值Dnα為? 0.175 19。由表2可見，對(duì)于4種概率分布函數(shù)，均存在Dnα>Dn，說(shuō)明文中擬合的降水量理論分布與實(shí)際頻率分布差異較小，可以通過(guò)假設(shè)，即宜賓市春季降水量服從正態(tài)分布、Gamma分布、對(duì)數(shù)正態(tài)分布和P-Ⅲ型分布。為檢驗(yàn)各自概率分布函數(shù)的擬合效果，分別畫出正態(tài)分布、Gamma分布、對(duì)數(shù)正態(tài)分布和P-Ⅲ分布對(duì)宜賓市春季降水量的擬合曲線，如圖1、圖2、圖3和圖4所示。從圖1可以看出，正態(tài)分布在40～60 mm降水量區(qū)段的擬合效果較好，而在其他降水量區(qū)段的擬合效果相對(duì)較差;從圖2、圖3可以看出，Gamma分布和對(duì)數(shù)正態(tài)分布在30～90 mm降水量區(qū)段的擬合效果較好，在90 mm降水量以上區(qū)段的擬合效果相對(duì)較差;從圖4可以看出，P-Ⅲ型分布在40～90 mm降水量區(qū)段的擬合效果較好，而在90～120 mm降水量區(qū)段的擬合效果較差?？傮w而言，對(duì)數(shù)正態(tài)分布和Gamma分布的整體擬合效果較好，P-Ⅲ型分布和正態(tài)分布的擬合效果相對(duì)較差，此為最優(yōu)擬合分布的判斷提供了一定的參考價(jià)值。

2.3? 對(duì)比分析

為了準(zhǔn)確判定最優(yōu)的擬合分布函數(shù)，通過(guò)離差平方和最小準(zhǔn)則法進(jìn)行宜賓市春季降水量擬合最佳的分布函數(shù)分析。當(dāng)OLS值最小，則其對(duì)應(yīng)的分布就為擬合效果最佳的分布。由表3可知，對(duì)數(shù)正態(tài)分布的OLS值最小，說(shuō)明對(duì)數(shù)正態(tài)分布對(duì)春季降水量的擬合效果最好，與實(shí)際宜賓市的春季降雨情況最接近。

為更加直觀地體現(xiàn)這一點(diǎn)，作出各個(gè)降水量概率分布函數(shù)擬合效果的對(duì)比分析圖。從圖5可以看出，正態(tài)分布概率理論曲線和經(jīng)驗(yàn)分布曲線擬合的效果最差;對(duì)數(shù)正態(tài)分布、P-Ⅲ型分布及Gamma分布的頻率理論曲線和經(jīng)驗(yàn)分布曲線都擬合的較好，和OLS結(jié)果值相吻合。因此可知，對(duì)數(shù)正態(tài)分布模擬的降雨和實(shí)際的基本一致，可以用來(lái)模擬實(shí)際的降雨情況，這個(gè)分布函數(shù)模型是合理的。

3? 結(jié)論

本研究根據(jù)四川省宜賓市的春季降水量資料，利用矩法估計(jì)、極大似然估計(jì)的參數(shù)估計(jì)方法，分別采用正態(tài)分布以及Gamma分布、對(duì)數(shù)正態(tài)分布和P-Ⅲ型分布對(duì)春季降水量分布進(jìn)行擬合。并且利用K-S檢驗(yàn)法對(duì)各分布函數(shù)進(jìn)行檢驗(yàn)，同時(shí)利用離差平方和最小準(zhǔn)則（OLS）對(duì)正態(tài)分布和3種非正態(tài)分布（對(duì)數(shù)正態(tài)分布、Gamma分布和P-Ⅲ型分布）的擬合效果做了對(duì)比分析，由計(jì)算結(jié)果可以得到以下主要結(jié)論。

1）正態(tài)分布、Gamma分布、對(duì)數(shù)正態(tài)分布和P-Ⅲ型分布均可以對(duì)宜賓市的春季降水量分布進(jìn)行擬合。

2）計(jì)算結(jié)果顯示，對(duì)數(shù)正態(tài)分布的OLS值最小，即對(duì)數(shù)正態(tài)分布對(duì)宜賓市春季降水量的擬合效果最好;Gamma分布次之，P-Ⅲ型分布再次，正態(tài)分布最差。

3）通過(guò)本次研究，可以為宜賓市的氣候統(tǒng)計(jì)診斷和預(yù)報(bào)工作提供一定的參考，為本地春季農(nóng)耕提供理論上的支撐，具有理論意義和業(yè)務(wù)實(shí)用價(jià)值。

參考文獻(xiàn)：

[1] 胡文東，陳曉光，李艷春，等.寧夏月、季、年降水量正態(tài)性分析[J].中國(guó)沙漠，2006，26（6）：963-968.

[2] 周雁翎，張瓊雄，張? 練，等.粵西降水量的分布特征及正態(tài)分布檢驗(yàn)[J].廣東氣象，2007，29（3）：32-34.

[3] 吳? 慧.海南省降水量的正態(tài)分布特征及正態(tài)化變換[J].廣東氣象，2005（2）：12-13.

[4] ENTEKHABI，D，EAGLESON P S. Land surface hydrology parameterization for AGCM including subgrid scale spatial variability[J].J Climate，1989，2：816-830.

[5] 劉晶淼，丁裕國(guó)，周秀驥，等.基于降水氣候強(qiáng)迫的一種地表徑流估計(jì)方法[J].氣象學(xué)報(bào)，2004，62（3）：285-293.

[6] 王紀(jì)軍，劉晶淼，丁裕國(guó).用于非均勻地表通量估計(jì)的一種積分算子[J].高原氣象，2004，23（5）：605-611.

[7] 劉晶淼，丁裕國(guó)，李? 云，等.一種度量區(qū)域旱澇災(zāi)害嚴(yán)重性的指標(biāo)[J].自然災(zāi)害學(xué)報(bào)，2004，13（5）：16-19.

[8] 段莉珠，番朝成.德宏潞西氣象站降水量統(tǒng)計(jì)分析[J].云南水力發(fā)電，2009，25（1）：8-10.

[9] 吳樟榮.基于P-Ⅲ型分布函數(shù)的臺(tái)風(fēng)風(fēng)險(xiǎn)水平分析[J].浙江師范大學(xué)學(xué)報(bào)，2002，25（2）：166-168.

[10] 周建康，唐運(yùn)憶，徐志俠.南京站降水量的統(tǒng)計(jì)分析[J].水文，2003，23（6）：35-38.

[11] 丁裕國(guó).降水量概率分布的一種間接模式[J].南京氣象學(xué)院學(xué)報(bào)，1987，4（10）：407-415.

[12] MOOLEY D A. Gamma probability model for asian summer monsoon monthly rainfall[J].Monthly weather review，1973，101：160-170.

[13] 李偉光，朱乃海，陳匯林，等.海口不同時(shí)間尺度降水量統(tǒng)計(jì)擬合[J].安徽農(nóng)學(xué)通報(bào)，2007，13（12）：189-190，227.

[14] 石? 強(qiáng)，林滋新.青島地區(qū)降水的Gamma分布模式[J].山東氣象，1987（2）：15-19.

[15] 丁裕國(guó)，申紅艷，江志紅，等.氣候概率分布理論及其應(yīng)用新進(jìn)展[J].氣象科技，2009，37（3）：257-262.