數(shù)學(xué)建模論文1.doc

河流流域水量調(diào)度技術(shù)方法應(yīng)用研究摘要本文首先分析了題目條件與表中數(shù)據(jù)，根據(jù)表中數(shù)據(jù)和所給的條件建立了有關(guān)串并聯(lián)、混聯(lián)和優(yōu)化調(diào)度的幾種模型，并對其進(jìn)行誤差分析，逐步改進(jìn)，找出比較優(yōu)秀的模型。運(yùn)用數(shù)據(jù)擬合，最小乘法等常用方法對數(shù)據(jù)進(jìn)行分析和建立模型，同時還利用MATLAB 的強(qiáng)大圖形處理功能，進(jìn)行數(shù)據(jù)的可視化，根據(jù)可視化的結(jié)果做出分析和判斷。并運(yùn)用回歸分析對所建立模型進(jìn)行有效檢驗，最后進(jìn)行模型的評價和推廣。同時提出了一些實用性的建議.問題1：我們根據(jù)題目的要求和圖形的分析建立了串并聯(lián)的模型，模型的建立有利于水資源的調(diào)度.問題2：在分析水庫群聯(lián)接方式及特點的基礎(chǔ)上，根據(jù)各水庫庫容大小、調(diào)節(jié)程度和水文上的差別等因素，建立模型為. 在對模型改進(jìn)的基礎(chǔ)上建立了模型II。對模型進(jìn)行了合理的理論證明和推導(dǎo)，所給出的理論證明結(jié)果大約為.，然后借助于.數(shù)學(xué)算法和MATLAB和lingo等軟件，對附件中所提供的數(shù)據(jù)進(jìn)行了篩選，去除異常數(shù)據(jù),對殘缺數(shù)據(jù)進(jìn)行適當(dāng)補(bǔ)充,并從中隨機(jī)抽取了3組數(shù)據(jù)（每組8個采樣）對理論結(jié)果進(jìn)行了數(shù)據(jù)模擬，結(jié)果顯示，理論結(jié)果與數(shù)據(jù)模擬結(jié)果吻合。問題3：我們建立在對問題二深入的分析，建立.關(guān)鍵字：1 問題重述（1）針對撫河流域水資源調(diào)度問題，建立干流串聯(lián)或混聯(lián)模型及支流并聯(lián)模型。（2）根據(jù)劃分區(qū)域分析典型的控制斷面，利用分析得到的影響因素、相關(guān)信息，應(yīng)用干流、支流兩個模型，以1月作為開始調(diào)度期，建立合理的數(shù)學(xué)模型，實現(xiàn)水庫群的聯(lián)合調(diào)度解決方案。（3）對于旱季（7月-9月），給出水量調(diào)度整體的指導(dǎo)性意見2 問題分析（1） 問題一：這是一個類似物理電路中各用電器的連接的問題.題目中主要研究的是：針對撫河流域水資源調(diào)度問題，搞清楚各水庫的關(guān)系有利于在干旱時期能給予需要合理調(diào)配用水，充分調(diào)動各方面水利資源，動態(tài)分配用水量. 要分析個水庫或河段之間的串并聯(lián)之間的關(guān)系,關(guān)機(jī)那是看,它們是否有水力上的直接聯(lián)系.如果個河段或水庫是布置在同一條河流，形如梯級的水庫群。各水庫有各自的開發(fā)目標(biāo)，同時各庫的徑流之間通過人工（天然）河渠形成上下聯(lián)系.,則為串聯(lián)方式.如果各水庫是位于相鄰的幾條干支流或不同河流上的一排水庫，它們有各自的集水面積，故無水力上的聯(lián)系。當(dāng)為同一目標(biāo)共同工作時，才有水利上的聯(lián)系，則為并聯(lián)方式.（2）問題二： 本問題是一個優(yōu)化問題。目的是要求在滿足題目要求先滿足生活用水，后滿足農(nóng)業(yè)用水，其次是工業(yè)用水，兼顧生態(tài)環(huán)境，進(jìn)行區(qū)域性的限制用水，并啟動相應(yīng)水庫調(diào)度措施，進(jìn)而實現(xiàn)對相應(yīng)河段的水量進(jìn)行實時調(diào)控、均勻分配水量資源。取得優(yōu)先分配原則和用水總量最小.為了便于分析，所以把整個流域分化成了六個區(qū)域。區(qū)域S1s1:控制斷面1以上區(qū)域S2s2:控制斷面1-3及2以上區(qū)域S3s3：控制斷面3-4區(qū)域S4s4：控制斷面4-9及6、7、8 以上區(qū)域S5s5：控制斷面5以上區(qū)域S6s6：控制斷面6以上 通過對降雨量、水庫水量損失、水庫的調(diào)節(jié)能力和需求量等因素的分析與限制，本題就變成了一個帶約束條件的優(yōu)化問題。（2） 問題三：我們通過對問題二建立了.模型，分析了.3 模型的假設(shè)（1）假設(shè)每月按三十天計算（2）4 定義與符號說明（1）為各測試站點同期降水量（mm）（2）P流域平均降水量（mm）（3）n測站數(shù)（4）計算時段 內(nèi)庫區(qū)水面蒸發(fā)強(qiáng)度，以水層深度（mm）計；（5）計算時段 內(nèi)庫區(qū)陸面蒸發(fā)強(qiáng)度，以水層深度（mm）計；（6） 計算時段 內(nèi)水庫平均水面面積（km2）（7）水面蒸發(fā)皿折算系數(shù)，一般為0.650.80。5 模型的建立與求解5.1我們通過對題目的理解與圖形的分析得到了問題一的兩個模型如下圖所示.（圖一、圖二）5.1.1 干流串聯(lián)模型圖二 從圖一可看出：四個控制段面把A河分為三個河段分別為S1、S2、S3。我們把這三個部分是干流上的串聯(lián)方式.由此我們建立串聯(lián)模型如圖一. 干流串聯(lián)方式最主要的特點是上下河段的水力聯(lián)系，主要表現(xiàn)在以下幾方面：（1）位于上游的、較高調(diào)節(jié)程度的河段對天然來水起了調(diào)節(jié)作用，改變了下游河段的流量在時間上的變化過程。（2）串聯(lián)的河段的天然入流往往同步性較好，因而上游河段枯水期調(diào)節(jié)而提高的下泄流量及洪水期減少的下瀉流量，多半對減輕下游河段的調(diào)節(jié)任務(wù)有利。（3）當(dāng)串聯(lián)的各河段為同一用水部門服務(wù)時，它們的工作便具有互為補(bǔ)充的關(guān)系。5.1.2 支流并聯(lián)模型圖 二從圖二中我們可以看出：大型水庫1、中型水庫4、中型水庫5、C河、D河，這五的關(guān)系就屬于支流上的并聯(lián)方式，由此我們建立支流上的并聯(lián)模型如圖二.支流并聯(lián)最主要的特點是各庫間與河流間的水利聯(lián)系，主要特點如下：（1）并聯(lián)水庫群和與河流的天然入流往往同步性較差，因而水庫間“水文補(bǔ)償”的調(diào)節(jié)作用更明顯。（2）當(dāng)并聯(lián)水庫群與河流為同一開發(fā)目標(biāo)服務(wù)時，它們的工作便具有互為補(bǔ)充的關(guān)系。5.2要實現(xiàn)水庫群的聯(lián)合調(diào)度解決方案，供水量與需水量必須滿足必須滿足一定的約束條件。 5.2.1 流域降雨量的計算聚雨量對優(yōu)化的影響.流域平均降水量的計算方法主要有：算術(shù)平均法、加權(quán)平均法、泰森多邊形法、等雨量線法等。在面積較大的流域，最好用泰森多邊形法，計算流域的平均降水量：小流域常用加權(quán)平均法；在平地上可用算術(shù)平均法和等雨量線法。.我們通過對上表的分析計算出了各月的聚雨量，見圖三 圖三.我們所知道的降雨量計算方法中算術(shù)平均法，即：對于地形起伏不大，降雨分布均勻，測試站布設(shè)合理或較多的情況下，算術(shù)平均法簡單、而且也能獲得滿意的結(jié)果。式中：為各測試站點同期降水量（mm）P流域平均降水量（mm）n測站數(shù)故我們以下采用算術(shù)平均法計算河流流域各控制斷面的聚水量。結(jié)果分析：根據(jù)數(shù)據(jù)表的數(shù)據(jù)得出多年平均降水量為1705.382mm。降水的年內(nèi)、年際分配極不均勻。年內(nèi)降水主要集中在46月份，占全年的47.22%以上.（太少）5.2.2 水庫的水量損失水庫建成蓄水后，因改變河流天然狀況及庫內(nèi)外水力條件而引起額外的水量損失，主要包括蒸發(fā)損失和滲透損失。 水庫蓄水后，使庫區(qū)形成廣闊水面，原有的陸面蒸發(fā)變?yōu)樗嬲舭l(fā)。由于流入水庫的徑流資料是根據(jù)建庫前壩址附近觀測資料整編得出，其中已計入陸面蒸發(fā)部分。因此，計算時段t（年、月）水庫的蒸發(fā)損失是指由陸面面積變?yōu)樗婷娣e所增加的額外蒸發(fā)量 （以m3計），即 （21） 式中：計算時段 內(nèi)庫區(qū)水面蒸發(fā)強(qiáng)度，以水層深度（mm）計； 計算時段 內(nèi)庫區(qū)陸面蒸發(fā)強(qiáng)度，以水層深度（mm）計； 計算時段 內(nèi)水庫平均水面面積（km2）； 1000單位換算系數(shù)，1 mmkm2=106/103 m3=103 m3。 水庫水面蒸發(fā)可根據(jù)水庫附近蒸發(fā)站或氣象站蒸發(fā)資料折算成自然水面蒸發(fā)，即 （22）式中：水面蒸發(fā)皿實測水面蒸發(fā)(mm)； 水面蒸發(fā)皿折算系數(shù)，一般為0.650.80。陸面蒸發(fā)，尚無較成熟的計算方法，在水庫設(shè)計中常采用多年平均降雨量和多年平均徑流深之差，作為陸面蒸發(fā)的估算值。有我國各個地方的觀測資料，由蒸發(fā)皿實測的水面蒸發(fā)量為1506mm，蒸發(fā)皿折算系數(shù)K=0.8，流域多年平均年降水量H= 15348.44mm。解: （1）陸面蒸發(fā)量：降雨量 - 徑流量 （2）水面蒸發(fā)量： （3）水庫的各月額外蒸發(fā)損失量:如下表：月蒸發(fā)量（降水量-徑流量）（億/m3）月份1月2月3月4月5月6月7月8月9月10月11月12月年蒸發(fā)量s1:控制面1以上1.462.74.44.865.965.532.1862.581.671.361.10.834.5s2:控制面 1-3及2以上-1.5-2-4-6.5-7.9-10-5.56-2.8-2.6-2.4-2-2-49s3:控制面3-4-2.9-4-8-13-16-19-11.2-5.5-4.7-4.4-3-3-95s4:控制面4-9及6、7、8以上-7.2-11-21-32-39-45-25.3-13-11-9.3-8-7-228s5:控制面5以上2.23-3-5-7.7-8.7-9.8-5.07-2.9-2.5-1.9-2-1-51s6:控制面6以上0.36-0-0-0.7-0.8-0.9-0.51-0.3-0.2-0.2-0-0-4.7另一方面，建庫之后，由于水庫蓄水，水位抬高，水壓力的增大改變了庫區(qū)周圍地下水的流動狀態(tài)，因而產(chǎn)生了水庫的滲漏損失。水庫的滲漏損失主要包括下面幾個方面：（l）通過能透水的壩身(如土壩、堆石壩等) 的滲漏，以及閘門、水輪機(jī)等的漏水；（2）通過壩基及繞壩兩翼的滲漏；（3）通過庫底、庫周流向較低的透水層的滲漏。一般可按滲漏理論的達(dá)西公式估算滲漏的損失量。計算時所需的數(shù)據(jù)(如滲漏系數(shù)、滲徑長度等)必須根據(jù)庫區(qū)及壩址的水文地質(zhì)、地形、水工建筑物的型式等條件來決定，而這些地質(zhì)條件及滲流運(yùn)動均較復(fù)雜，往往難以用理論計算的方法獲得較好的成果。因此，在生產(chǎn)實際中，常根據(jù)水文地質(zhì)情況，定出一些經(jīng)驗性的數(shù)據(jù)，作為初步估算滲漏損失的依據(jù)。若以一年或一月的滲漏損失相當(dāng)于水庫蓄水容積的一定百分?jǐn)?shù)來估算時，則采用如下數(shù)值：（l）水文地質(zhì)條件優(yōu)良(指庫床為不滲水層，地下水面與庫面接近)，010%/年或01%/月。（2）透水性條件中等，10%20%/年或1%1.5%/月。（3）水文地質(zhì)條件較差，20%40%/年或1.5%3%/月。在水庫運(yùn)行的最初幾年，滲漏損失往往較大(大于上述經(jīng)驗數(shù)據(jù))，因為初蓄時，為了濕潤土壤及抬高地下水位需要額外損失水量。水庫運(yùn)行多年之后，因為庫床泥沙顆粒間的空隙逐漸被水內(nèi)細(xì)泥或粘土淤塞，滲漏系數(shù)變小，同時庫岸四周地下水位逐漸抬高，滲漏量減少?，F(xiàn)根據(jù)四水庫的水文地質(zhì)條件得出其透水性條件為每月滲透1.25%。得出以下表格：水庫名稱正常蓄水位(m)有效庫容(億m3)死水位(m)死庫容(億m3)調(diào)節(jié)庫容(億m3)1月份初始水位（m）一月份初庫容量(億m3)平均每月滲透量(億m3)大型水庫111404.49710761.1363.36110851.43320.17915大型水庫29704.649361.742.99392.080.26中型水庫38370.6688220.1010.5678240.1720.0215中型水庫47600.8147200.1120.7027230.1650.0206255.2.3 水位和流量 原始數(shù)據(jù)給出的表格中只是包含了：大型水庫一：9761007；大型水庫二：836871；中型水庫三：760.5794；中型水庫四：626670；但是未能給出各大水庫的正常水位對應(yīng)的流量。于是運(yùn)用了二次項擬合，得出了以下的數(shù)據(jù)和表格：水位與流量圖水位與流量的關(guān)系大型水庫1 水位（m）976977978979980981982983984985流量 （m3/s）2762.41161872763875276808381020水位（m）986987988989990995998100010051007流量（m3/s）1220142016301870211034504520523071608000水位（m）1010101510201025103010351040104510501055流量（m3/s）9502 12155 15128 18421 22034 25967 30220 34793 39686 44899 水位（m）1060106510701075108010851090109511001105流量（m3/s）50433 56286 62460 68953 75767 82901 90355 98129 106223 114637 水位（m）1110111511201125113011351140流量（m3/s）123371 132425 141800 151494 161509 171844 182498 大型水庫2水位（m）836836.5837837.5838838.5839840842845流量（m3/s）01025355580105180420900水位（m）847850852855857860862865871877流量（m3/s）12501850225057706870867010030121301662023865 水位（m）883889895901907913919925931937流量（m3/s）31403 39970 49568 60195 71852 84539 98255 113002 128778 145583 水位（m）943流量（m3/s）163419 中型水庫3水位（m）760.5761763765767769771773775777流量（m3/s）2626.965.3104153211277353445564水位（m）779781783785787788789790794798流量（m3/s）70384910001160132014101500159019502443 水位（m）802806810814818822826830834838流量（m3/s）2943 3489 4083 4722 5409 6142 6921 7747 8620 9539 水位（m）842流量（m3/s）10505 中型水庫4水位（m）626626.5627628629630631632633634流量（m3/s）621372234306769351210150018002110水位（m）635636637639640642650660670680流量（m3/s）2420274030803810420050008400123001650021772 水位（m）690700710720730740750760770780流量（m3/s）27283 33245 39659 46524 53840 61609 69828 78499 87622 97196 水位（m）790流量（m3/s）107222 5.2.4 徑流調(diào)節(jié) （不是很清楚想干嘛，說的不明白）眾所周知，河川徑流在一年之內(nèi)或者在年際之間的豐枯變化都是很大的。我國河流年內(nèi)洪水季的水量往往要占全年來水總量的70%80%。河川徑流的劇烈變化，給人類帶來很多不利的后果，如汛期大洪水容易造成災(zāi)害，而枯水期水少，不能滿足興利需要。因此，無論是為了消除或減輕洪水災(zāi)害，還是為了滿足興利需要，都要求采取措施，對天然徑流進(jìn)行控制和調(diào)節(jié)。圖四為年調(diào)節(jié)示意圖。圖上表明，只需一部分多余水量將水庫蓄滿（圖中橫線所示），其余的多余水量（斜線部分），只能由溢洪道棄掉。圖中豎影線部分表示由水庫放出的水量，以補(bǔ)充枯水季天然水量的不足，其總水量相當(dāng)于水庫的調(diào)節(jié)庫容。 水庫的興利庫容能夠蓄納設(shè)計枯水年豐水期的全部余水量時，稱為完全年調(diào)節(jié)；若興利庫容相對較小，不足以蓄納設(shè)計枯水年豐水期的全部余水量而產(chǎn)生棄水時，稱為不完全年調(diào)節(jié)，或季調(diào)節(jié)。這是規(guī)劃設(shè)計中劃分水庫調(diào)節(jié)性能所采用的界定。必須指出，從水庫實際運(yùn)行看，這種劃分是相對的，完全年調(diào)節(jié)遇到比設(shè)計枯水年徑流量更豐的年份，就不可能達(dá)到完全年調(diào)節(jié)。供水 蓄水 棄水 供水Q(m3/S) （月）圖四 q 5.2.5生活用水、工農(nóng)業(yè)等用水所占比列（用餅形圖畫出來接著還要分析）5.2.6 水庫的調(diào)節(jié)能里要分析（要用數(shù)據(jù)和圖形說明）5.27 模型求解模型一： 接下來根據(jù)已知條件建立聯(lián)合調(diào)度的數(shù)學(xué)模型，找出最準(zhǔn)則與目標(biāo)函數(shù)來是想水庫群的聯(lián)合調(diào)度解決方案。研究庫群聯(lián)合優(yōu)化調(diào)度必先確定最優(yōu)準(zhǔn)則。這里根據(jù)州河流域四庫水文特性及水資源利用現(xiàn)狀，確定最優(yōu)準(zhǔn)則為：四水庫控制并利用水資源量最大；四水庫內(nèi)水面蒸發(fā)量最小；四大水庫自起始調(diào)度（一月）月起，河流各月平均水面面積分別記作：S1，S2，S3；S4，S5，S6。各月平均水面蒸發(fā)系數(shù)與蒸發(fā)量分別記作：K1，K2，K6。E1，E2，E6；于橋水庫各月出庫水量分別記作：QO1，QO2，QO4。則由最優(yōu)準(zhǔn)則可寫出求解四庫水量損失的數(shù)學(xué)表達(dá)式，亦即目標(biāo)函數(shù)如下：Min【（S1+S2+S3+S6）KE/10+QOj】1，2，6, j=1,2,3,4約束方程（1）蓄水量約束PwQj 1，2，3,4；1，2，12式中，w為j水庫第個月份的庫蓄水量；P、Q分別為j水庫第個月份的答應(yīng)最小、最大庫蓄水量。（2）需供水量約束P+Uw 1，2，3；1，2，12式中，U為水庫第個月份需供水量，其它同上。（3）起調(diào)蓄量約束W111085; w21939; w31824; w41723; 式中，w1為水庫第1個月份庫蓄水量，均為起調(diào)蓄量。先確定四庫聯(lián)合調(diào)度的決策變量為：起調(diào)蓄量w1、入庫水量QIw、需供水量Uw、蒸發(fā)系數(shù)K、蒸發(fā)量E，則綜合目標(biāo)函數(shù)及約束方程，可獲得流域四庫聯(lián)合調(diào)度的優(yōu)化數(shù)學(xué)模型如下：MINEQMin【（S1+S2+S3）KE/10+QO3i】St w1w13W+UwW式中，1，2，3；1，2，12。由州河流域四水庫水文資料，可確定四水庫平均每年的起調(diào)蓄量、入庫水量、需供水量、蒸發(fā)系數(shù)及蒸發(fā)量等參數(shù)。根據(jù)已知資料，可確定大型水庫一、大型水庫二、中型水庫三、中型水庫四起調(diào)蓄量分別為： 1.4332、939.0、0.172、0.165（108m3）。運(yùn)用輪庫尋優(yōu)算法，可算得實現(xiàn)聯(lián)合調(diào)度時似水庫平均每年最小水量損失為81915.23(104m3)，最小蒸發(fā)量10960.25 (104m3)，于橋水庫最小放水量70954.98(104m3)。三庫優(yōu)化調(diào)度方案中蓄量、水位及水庫放水量(表.4)。 州河流域盡管有較強(qiáng)的產(chǎn)流能力，卻因流域面積中型水庫三和中型水庫四兩水庫庫容相對較小，聯(lián)合調(diào)度的難度較大。這里，從徑流調(diào)節(jié)及流域水資源聯(lián)合調(diào)度角度出發(fā)，運(yùn)用數(shù)學(xué)模型與輪庫尋優(yōu)算法，編寫計算機(jī)程序?qū)α饔騼?nèi)三庫進(jìn)行尋優(yōu)計算，使流域水資源損失減至最小，也為庫群聯(lián)合調(diào)度提供了實用的計算方法。模型二：設(shè)有4個水庫(大型水庫1、大型水庫2、中型水庫3、中型水庫4)，有15個區(qū)域記為1、2、3、4、5、6、7、8、9、10、11、12、13、14、15。區(qū)域地段需水量1月2月3月4月5月6月區(qū)域1控制斷面1以上0.0619370.0619370.0619370.0619370.0619370.061937區(qū)域1控制斷面1控制斷面30.2383610.2383610.2383610.2383610.2383610.238361區(qū)域1控制斷面3控制斷面40.4737610.4737610.4737610.4737610.4737610.473761區(qū)域1控制斷面4控制斷面90.19030.19030.19030.19030.19030.1903區(qū)域1控制斷面2以上0.08850.08850.08850.08850.08850.0885區(qū)域1控制斷面2河口0.0748010.0748010.0748010.0748010.0748010.074801區(qū)域1控制斷面5以上0.261210.261210.261210.261210.261210.26121區(qū)域1控制斷面5河口0.1172930.1172930.1172930.1172930.1172930.117293區(qū)域1控制斷面6以上0.0584280.0584280.0584280.0584280.0584280.058428區(qū)域1控制斷面6河口0.0332160.0332160.0332160.0332160.0332160.033216區(qū)域1控制斷面7河口0.2542680.2542680.2542680.2542680.2542680.254268區(qū)域1控制斷面8河口0.3357840.3357840.3357840.3357840.3357840.335784區(qū)域1灌溉區(qū)10.2153130.2153130.2153130.2153130.2153130.215313區(qū)域1灌溉區(qū)20.2110140.2110140.2110140.2110140.2110140.211014區(qū)域1灌溉區(qū)30.2328670.2328670.2328670.2328670.2328670.232867區(qū)域地段需水量7月8月9月10月11月12月區(qū)域1控制斷面1以上0.258190.258190.258190.0619370.0619370.061937區(qū)域1控制斷面1控制斷面30.721250.721250.721250.2383610.2383610.238361區(qū)域1控制斷面3控制斷面41.6563831.6563831.6563830.4737610.4737610.473761區(qū)域1控制斷面4控制斷面90.7287670.7287670.7287670.19030.19030.1903區(qū)域1控制斷面2以上0.3071670.3071670.3071670.08850.08850.0885區(qū)域1控制斷面2河口0.2695970.2695970.2695970.0748010.0748010.074801區(qū)域1控制斷面5以上1.0030371.0030371.0030370.261210.261210.26121區(qū)域1控制斷面5河口0.449120.449120.449120.1172930.1172930.117293區(qū)域1控制斷面6以上0.216050.216050.216050.0584280.0584280.058428區(qū)域1控制斷面6河口0.1226870.1226870.1226870.0332160.0332160.033216區(qū)域1控制斷面7河口1.0558631.0558631.0558630.2542680.2542680.254268區(qū)域1控制斷面8河口1.3943131.3943131.3943130.3357840.3357840.335784區(qū)域1灌溉區(qū)11.167061.167061.167060.2153130.2153130.215313區(qū)域1灌溉區(qū)21.155291.155291.155290.2110140.2110140.211014區(qū)域1灌溉區(qū)31.1977331.1977331.1977330.2328670.2328670.232867水庫入庫量1月2月3月4月5月6月大型水庫10.946081.062721.0964161.2545282.28614412.49344大型水庫21.7889981.8348771.7921092.3734944.7179589.400147中型水庫30.246240.2410560.1399680.2306882.877126.2208中型水庫40.142560.8786880.3188160.622081.477445.9616水庫入庫量7月8月9月10月11月12月大型水庫111.508483.991685.339522.1720961.9206721.194912大型水庫215.171244.173124.224962.151361.87921.757376中型水庫34.66561.373761.03680.362880.51840.274752中型水庫48.29440.6481.632960.414720.959041.01088