水利工程論文-淺析數(shù)據(jù)庫技術(shù)在海河流域水資源綜合規(guī)劃中的應(yīng)用論文關(guān)鍵詞：數(shù)據(jù)庫技術(shù)；水資源規(guī)劃；供需分析；海河流域論文摘要：海河流域水資源綜合規(guī)劃采用水資源系統(tǒng)配置模擬模型對全流域進(jìn)行了長系列的供需分析，生成了近2萬行供需平衡數(shù)據(jù)。如果采用常規(guī)手工統(tǒng)計(jì)分析方法進(jìn)行供需分析，難于適應(yīng)規(guī)劃的進(jìn)度和精度要求。介紹如何利用Foxpro數(shù)據(jù)庫系統(tǒng)提供的強(qiáng)大數(shù)據(jù)處理功能，實(shí)現(xiàn)供需平衡數(shù)據(jù)的統(tǒng)計(jì)分析和輸出，并給出了實(shí)用的程序片段，不僅對提高水資源規(guī)劃工作效率有所幫助，而且其思路和方法對其他專業(yè)規(guī)劃也具有一定的參考價(jià)值。1引言海河水利委員會在水資源綜合規(guī)劃中委托中國水科院水資源所開發(fā)研制了海河流域水資源系統(tǒng)配置模型,配置模型以水資源三級區(qū)套地市作為基本計(jì)算分區(qū)，同時(shí)又根據(jù)規(guī)劃精度的要求對分區(qū)作了進(jìn)一步細(xì)化，共劃分107個(gè)計(jì)算單元。利用配置模型對全流域進(jìn)行了長系列的供需分析生成了基準(zhǔn)年、2010、2020和2030年4個(gè)規(guī)劃水平年107個(gè)計(jì)算單元45年的供需平衡數(shù)據(jù)，不僅數(shù)據(jù)量大，而且方案反復(fù)調(diào)整,十分煩瑣，若采用常規(guī)手工統(tǒng)計(jì)分析方法，對計(jì)算單元供需平衡數(shù)據(jù)按不同分區(qū)進(jìn)行各種頻率的供需分析，難于適應(yīng)規(guī)劃的進(jìn)度和精度要求。運(yùn)用Foxpro數(shù)據(jù)庫系統(tǒng)提供的強(qiáng)大數(shù)據(jù)處理功能基本實(shí)現(xiàn)了不同水資源分區(qū)供需分析長系列成果生成、均值統(tǒng)計(jì)、各種保證率成果確定、大型調(diào)水工程不同保證率供水量統(tǒng)計(jì)及最終成果表編制等工作的程序化,極大地提高了工作效率。2建立基礎(chǔ)數(shù)據(jù)庫水資源配置模型運(yùn)行后導(dǎo)出的供需平衡數(shù)據(jù)為包含107個(gè)計(jì)算單元的1956-2000年的系列數(shù)據(jù)文件數(shù)據(jù)格式既可選擇Excel表格，也可選擇以Tab為分隔的數(shù)據(jù)文件表格基本格式見表1(僅列部分字段)。Foxpro數(shù)據(jù)庫系統(tǒng)具有直接導(dǎo)入以上兩種數(shù)據(jù)格式文件并生成相應(yīng)數(shù)據(jù)庫的功能。以Tba分隔數(shù)據(jù)文件為例，如水資源配置模型導(dǎo)出文件名為balancedat則導(dǎo)入并建立數(shù)據(jù)庫的程序?yàn)椋簎se單元供需分析dbfdelelallpackappendfrombalanceDattypedelimitedwithtab其中，“單元供需分析dbf”為事先建立其字段為與balancedat相一致的數(shù)據(jù)庫文件。按水資源綜合規(guī)劃要求供需分析成果以省級行政區(qū)套水資源二級區(qū)(以下簡稱省套二級區(qū))為統(tǒng)計(jì)單元保證率主要根據(jù)統(tǒng)計(jì)單元的降水頻率確定。因此需建立一個(gè)以省套二級區(qū)為統(tǒng)計(jì)單元的長系列降水量數(shù)據(jù)庫“降水排序dbf”字段包括：年份、二級區(qū)碼、省級碼、降水量和降水排序(降水量降序排列號)5個(gè)字段，其中年份為1956-2000年。為了簡化統(tǒng)計(jì)分析的程序代碼，構(gòu)建了一個(gè)貯存省級行政區(qū)及二級區(qū)組合碼的數(shù)據(jù)庫“省套二級區(qū)編碼dbf，數(shù)據(jù)庫記錄包括了海河流域省級行政區(qū)與二級區(qū)的所有可能組合海河流域有l(wèi)6個(gè)省套二級區(qū)，因此共有l(wèi)6個(gè)記錄。字段包括：wt_code(二級區(qū)碼)、pry_code(省級碼)、wt_na(二級區(qū)名稱)、plv_na(省級區(qū)名稱)4個(gè)字段。3數(shù)據(jù)庫記錄編碼供需平衡數(shù)據(jù)的統(tǒng)計(jì)分析，主要以計(jì)算單元編碼作為統(tǒng)計(jì)區(qū)域分類的線索，選擇合適的編碼方式和規(guī)則尤為重要。以字母和數(shù)據(jù)構(gòu)成的編碼具有簡潔、便于重新組合的優(yōu)點(diǎn)。在海河流域水資源配置模型中，計(jì)算單元編碼由1個(gè)字母和l5個(gè)數(shù)字組成，其中前3位代表水資源二級區(qū)，如C02為海河北系；前7位代表水資源三級區(qū)，如C010100為灤河山區(qū)；89位為省級行政區(qū)，如l3為河北省；而第l0位以后的編碼為地市級行政區(qū)等信息。對計(jì)算單元供需平衡數(shù)據(jù)按不同分區(qū)要求進(jìn)行統(tǒng)計(jì)分析的操作，實(shí)際上是對數(shù)據(jù)庫記錄(行)進(jìn)行分類、合并、排序、提取的過程。為了達(dá)到分類合并的目的，在“單元供需分析dbf”數(shù)據(jù)庫中增加了1個(gè)“組合碼”字段，用于存儲重新分類統(tǒng)計(jì)的組合編碼。4供需平衡數(shù)據(jù)的統(tǒng)計(jì)分析為了生成省套二級區(qū)的供需分析長系列成果，并進(jìn)行各種頻率的供需分析，首先需對“單元供需分析dbf”數(shù)據(jù)庫中的每個(gè)記錄進(jìn)行分類，如哪些記錄屬于海河北系河北省1956年的數(shù)據(jù)?通過提取單元編碼中代表二級區(qū)及省級行政區(qū)的編碼，并與年份進(jìn)行組合，即可形成不同年份省套二級區(qū)的組合碼，將其存人數(shù)據(jù)庫的組合碼字段，為同類合并(相加)創(chuàng)造了條件。數(shù)據(jù)庫系統(tǒng)對已索引或排序的數(shù)據(jù)庫提供了方便的同類合并功能，以下程序片段，輕而易舉地完成了省套二級區(qū)供需分析長系列數(shù)據(jù)的整理工作。use單元供需分析dbfcopyalltotempdbfusetempdbfreplall組合碼withyearl+substr(code0，1，3)+substr(code0，8，2)1indexon組合碼toch_codetotalon組合碼to省套二級區(qū)系列dbf其中：tempdbf為臨時(shí)數(shù)據(jù)庫，是“單元供需分析dbf”數(shù)據(jù)庫的復(fù)制文件對其進(jìn)行操作可避免改變原有數(shù)據(jù)庫：yearl和code0分別為年份和單元編碼字段。在數(shù)據(jù)庫系統(tǒng)的工作區(qū)1打開先前建立的“省套二級區(qū)編碼dbf”，并保持打開狀態(tài)。在工作區(qū)2反復(fù)打開或關(guān)閉用于查找不同降水頻率對應(yīng)年份的“降水排序dbf”、用于提取不同保證率供需平衡數(shù)據(jù)的“省套二級區(qū)系列dbf”及存儲統(tǒng)計(jì)分析結(jié)果的“省套二級區(qū)供需分析dbf”，并進(jìn)行數(shù)據(jù)的統(tǒng)計(jì)、復(fù)制、添加、分區(qū)名稱及保證率的填寫等操作。以工作區(qū)1第1個(gè)記錄(第一個(gè)省套二級區(qū))為起始，提取省套二級區(qū)信息，利用dowhilenoteofO循環(huán)語句。在兩個(gè)工作區(qū)之間來回變換，進(jìn)行重復(fù)的統(tǒng)計(jì)分析操作，直到完成最后一個(gè)省套二級區(qū)的統(tǒng)計(jì)分析。程序片段(僅以50保證率和平均值統(tǒng)計(jì)為例)如下：sele1&選擇工作區(qū)1use省套二級區(qū)編碼dbfdowhilenoteof0&循環(huán)開始prv_nal=prv_na&提取要統(tǒng)計(jì)的省級區(qū)名wt_na1=wt_na&提取要統(tǒng)計(jì)的二級區(qū)名prv_codel=pry_code&提取要統(tǒng)計(jì)的省級區(qū)編碼wt_codel=wt_code&提取要統(tǒng)計(jì)的二級區(qū)編碼sele2&選擇工作區(qū)2use降水排頻Dbf&在排頻文件中找到年份locatefor省級碼=prv_codeland二級區(qū)碼=wt_codeland降水排序=23year50=年份use省套二級區(qū)系列Dbf&在數(shù)據(jù)庫中提取對應(yīng)省套二級區(qū)和頻率的記錄Locateforsubstr(code0，1，3)=wt_code1andsubstr(code0，8，2)=prv_codelandyear0=year50Copynext1totmp50dbfaverageall城鎮(zhèn)需水to城鎮(zhèn)需水1forsubstr(code0，1，3)=wt_codelandsubstr(code0，8，2)=pry_co