石 琴，雒 儀

(甘肅省平?jīng)鏊乃Y源勘測局，甘肅 平?jīng)?744000)

1 測驗(yàn)系統(tǒng)簡介

1.1 水文軟件發(fā)展

20世紀(jì)80年代初，水文資料一直采用手工記載、計(jì)算整理，工作量大、錯(cuò)誤率高，直接影響整編效率。20世紀(jì)90年代后期，隨著計(jì)算機(jī)技術(shù)的引入，形成了較為成熟的DOS版水文測算軟件，但無可視化窗口，且無法對文字、圖形進(jìn)行處理[1]。進(jìn)入21世紀(jì)，隨著水文測驗(yàn)儀器設(shè)備自動(dòng)化水平的提高，出現(xiàn)基于Visual Basic.NET編程語言的水文測算軟件，可以處理基本水文資料[2]。但隨著水文行業(yè)規(guī)范的修訂完善、計(jì)算機(jī)更新?lián)Q代周期縮短和開發(fā)者后期缺乏維護(hù)等原因，此軟件已不滿足當(dāng)前水文規(guī)范的要求。依據(jù)計(jì)算機(jī)技術(shù)發(fā)展，甘肅省水文局采用先進(jìn)的Java編程語言，研發(fā)的《甘肅省水文測驗(yàn)系統(tǒng)》(簡稱測驗(yàn)系統(tǒng))軟件，滿足當(dāng)前水文測驗(yàn)整編規(guī)范與計(jì)算機(jī)操作系統(tǒng)，并在一定時(shí)間內(nèi)保持領(lǐng)先水平。本文通過測驗(yàn)系統(tǒng)在靈臺(tái)水文站在站資料計(jì)算整編中的應(yīng)用，分析其對于基層測站的適用性，并介紹使用中的技術(shù)要點(diǎn)。

1.2 測驗(yàn)系統(tǒng)特性

測驗(yàn)系統(tǒng)集成了水文數(shù)據(jù)分析、圖形顯示、文檔讀寫和人機(jī)交互，各項(xiàng)功能之間通過My SQL數(shù)據(jù)庫鏈接[3]；具有數(shù)據(jù)錄入、分析計(jì)算、圖形處理、成果導(dǎo)出等功能。其中分析計(jì)算是以水文行業(yè)規(guī)范為標(biāo)準(zhǔn)編譯的資料分析計(jì)算模塊，主要包含自動(dòng)監(jiān)測資料分析、流量分析計(jì)算、水位流量關(guān)系曲線等水文資料處理相關(guān)功能。原始資料由軟件界面以人工輸入或文本導(dǎo)入等形式取得，保存于數(shù)據(jù)庫分類表中，通過測驗(yàn)系統(tǒng)分析計(jì)算后形成整編成果，可導(dǎo)出電子文檔或打印輸出。各類成果數(shù)據(jù)均符合水文行業(yè)技術(shù)規(guī)范標(biāo)準(zhǔn)[4]和《水文資料整編規(guī)范》(SL 247—2012)的要求。

2 靈臺(tái)水文站簡介

靈臺(tái)水文站位于平?jīng)鍪徐`臺(tái)縣中臺(tái)鎮(zhèn)，地理位置為東經(jīng)107°36′15″，北緯35°04′13″，地處黃河流域涇河水系達(dá)溪河，集水面積為1 500 km2，距河口為33.5 km，屬于國家重要水文站、區(qū)域代表站、中央報(bào)汛站、省界斷面站和國家水質(zhì)監(jiān)測站。隸屬于甘肅省平?jīng)鏊乃Y源勘測局管理，管轄雨量站5處，水位站3處，定點(diǎn)洪水調(diào)查河段1處，水質(zhì)調(diào)查斷面1處，測驗(yàn)任務(wù)包含測站斷面和地形測量、蒸降、水位、流量、泥沙、暴雨洪水調(diào)查等項(xiàng)目。

3 測驗(yàn)系統(tǒng)應(yīng)用

測驗(yàn)系統(tǒng)采用純組件UI，精簡平面風(fēng)格。界面左側(cè)為控制區(qū)域，右側(cè)為計(jì)算與顯示區(qū)域。左中部為程序各主要項(xiàng)目的選擇按鈕，其中包括“測站信息”、“斷面數(shù)據(jù)”、“流量測驗(yàn)”、“流量成果”等項(xiàng)目。左上部為各主要項(xiàng)目計(jì)算控制區(qū)域，包括提示說明以及各計(jì)算控制按鈕。左下部為其他按鈕，包括快捷操作、聲明、設(shè)置等。軟件最下方為狀態(tài)條，包括軟件版本、提示欄。測驗(yàn)系統(tǒng)軟件界面見圖1所示。

3.1.1測站信息

測站信息模塊是測驗(yàn)系統(tǒng)數(shù)據(jù)庫的基本組成體系，在使用軟件前需先建立測站信息，信息包括：測站站碼、站名、水系信息與各類計(jì)算系數(shù)。其中至少錄入“序號(hào)”、“測站編碼”、“站名”3個(gè)信息，該測站才會(huì)被程序保存與調(diào)用，且“測站編碼”與“站名”必須是唯一值。測站信息錄入界面見圖2所示。

圖2 測站信息錄入界面示意

3.1.2制表信息

當(dāng)對資料成果進(jìn)行打印輸出時(shí)，會(huì)涉及到各類成果表的制表信息和輔助文本。測驗(yàn)系統(tǒng)支持對成果表中各項(xiàng)輔助說明文本的填寫，分為兩部分。一處位于“制表信息”模塊中，主要為測算人員、輔助文本等全局性信息；另一處位于測驗(yàn)系統(tǒng)各項(xiàng)計(jì)算功能中，主要為測算時(shí)間、設(shè)備型號(hào)等獨(dú)立信息，其中流速儀法計(jì)算制表信息見圖3所示。

圖3 流速儀法計(jì)算制表信息

3.2 自動(dòng)監(jiān)測資料分析

靈臺(tái)水文站自動(dòng)監(jiān)測設(shè)備主要集中在水位和降水量測驗(yàn)項(xiàng)目，所采用的儀器分別是SEBA PULS型雷達(dá)水位計(jì)和JH-2型雨量計(jì)，2款設(shè)備均以5 min為最小單位依次采集數(shù)據(jù)，單個(gè)站點(diǎn)全年可觀測105 120次，人工無法處理。測驗(yàn)系統(tǒng)可實(shí)現(xiàn)自動(dòng)監(jiān)測設(shè)備資料的精簡摘錄，剔除資料中存在的重復(fù)值、特異值、無效值、亂碼等問題，資料滿足與《水文資料整編系統(tǒng)》HDP5.0的對接，從而實(shí)現(xiàn)自動(dòng)監(jiān)測資料整編[5]。自動(dòng)監(jiān)測資料可通過“宏電RTU前置平臺(tái)”遠(yuǎn)程下載，首先確保采集的原始資料未改動(dòng)，在操作界面選擇對應(yīng)的測站和需要分析的年份，單擊“載入”將原始資料添加至測驗(yàn)系統(tǒng)，此時(shí)程序提示“原始資料加載完成”，最后經(jīng)“分析”可得到處理完成的數(shù)據(jù)，再通過《水文資料整編系統(tǒng)》HDP5.0整編后可得逐日降水量表、降水量摘錄表、逐日平均水位表、洪水水位摘錄表等成果資料。

以降水量資料整編為例，選擇2020年度主汛期JH-2型雨量計(jì)資料與人工記錄進(jìn)行比測分析，靈臺(tái)水文站降水量總樣本為64個(gè)，各月總量差值為0.4～2.0 mm，比值在1.5%～3.0%區(qū)間；最大日量差為0.1～2.5 mm；降水量日數(shù)差為0～2 d，主要因素是JH-2型雨量計(jì)的最小感應(yīng)量為0.2 mm，無法識(shí)別0.1 mm降水。綜上分析，JH-2型雨量計(jì)資料誤差符合標(biāo)準(zhǔn)，通過測驗(yàn)系統(tǒng)處理后的資料達(dá)到整編規(guī)范要求。靈臺(tái)水文站主汛期降水量對照見圖4。

圖4 靈臺(tái)水文站主汛期降水量對照示意

3.3 流量資料分析

測驗(yàn)系統(tǒng)分別支持對流速儀法、浮標(biāo)法、比降面積法的流量分析計(jì)算。在靈臺(tái)水文站流量測驗(yàn)項(xiàng)目中，一般在水位低于996.32 m時(shí)屬于低水位區(qū)間，此時(shí)采用流速儀法涉水測流。測驗(yàn)系統(tǒng)通過分析法計(jì)算流速儀采集的數(shù)據(jù)，計(jì)算內(nèi)容主要包括測點(diǎn)流速、垂線平均流速、部分平均流速、部分流量、斷面流量等；當(dāng)水位位于996.33～997.21 m時(shí)屬于中高水位區(qū)間，此時(shí)水面流速波動(dòng)較為平穩(wěn)，適合采用雷達(dá)波流速儀進(jìn)行測算；當(dāng)水位高于997.21 m時(shí)屬于高水位，此時(shí)水面波動(dòng)較為明顯，受測驗(yàn)條件與河流特性的制約，采用浮標(biāo)法進(jìn)行流量測驗(yàn)。測驗(yàn)系統(tǒng)采用內(nèi)插法對水邊加以確定，分析過水?dāng)嗝婷娣e，根據(jù)實(shí)測水面流速計(jì)算斷面流量[6]；當(dāng)汛后進(jìn)行定點(diǎn)洪水調(diào)查時(shí)，測驗(yàn)系統(tǒng)可通過實(shí)測河段的水面比降、大斷面等數(shù)據(jù)，采用曼寧公式來推算歷史最大流量。本小結(jié)主要以流速儀法流量計(jì)算為例，簡要闡述測驗(yàn)系統(tǒng)在流量分析計(jì)算中的應(yīng)用。

3.3.1流速儀法計(jì)算

流速儀法流量計(jì)算模塊能夠處理測深垂線單點(diǎn)法與多點(diǎn)法，并支持任意的測速垂線布設(shè)，可通過大斷面與水位自動(dòng)提取過水?dāng)嗝娣秶ＶС帜媪?、停滯、分流等特殊水情，適用范圍廣，計(jì)算后可得《暢流期流速儀法實(shí)測流量記載表》[7]。由于靈臺(tái)水文站斷面控制較好，在測算中可以借用大斷面來布設(shè)測流斷面，先將測站大斷面資料填寫至“斷面數(shù)據(jù)”，返回“流速儀法”進(jìn)行操作。

1) 單擊“斷面數(shù)據(jù)”選擇所需要借用的斷面名稱(見圖5)。

圖5 測驗(yàn)系統(tǒng)“流速儀法”模塊計(jì)算選項(xiàng)

2) 選擇“垂線測次”，下拉可選“一點(diǎn)法”、“二點(diǎn)法”與“精測法”。“垂線測次”指斷面自動(dòng)布設(shè)時(shí)垂線的測點(diǎn)數(shù)，具體規(guī)則見圖6。

圖6 測驗(yàn)系統(tǒng)“垂線測次”參數(shù)設(shè)置選項(xiàng)

3) 填寫“開始水位”與“結(jié)束水位”，其為斷面布設(shè)的截止高程，也就是水面所對應(yīng)的水位。“開始水位”必須符合所選擇斷面的特性，不得大于所選擇斷面的“最高”高程或小于“最低”高程，否則系統(tǒng)將不予布設(shè)斷面。

4) 以上數(shù)據(jù)填寫完成后單擊“布設(shè)”，右側(cè)計(jì)算表格中將會(huì)顯示提取的斷面數(shù)據(jù)以供使用，將會(huì)自動(dòng)生成“號(hào)數(shù)”、“起點(diǎn)距”、“測流水位”、“河底高程”、“水深”、“相對”、“測點(diǎn)深”(見圖7)。手工將實(shí)測“N”、“T”或“測點(diǎn)流速”填寫進(jìn)行后續(xù)流量計(jì)算與成果表打印輸出。

圖7 測驗(yàn)系統(tǒng)“斷面布設(shè)”成果輸出界面

3.4 水位流量曲線分析

水位流量關(guān)系曲線分析適用于河道水情報(bào)汛中臨時(shí)曲線的率定，借用實(shí)時(shí)測驗(yàn)的流量數(shù)據(jù)可快速率定關(guān)系曲線，在特殊情況下測站無法開展流量測驗(yàn)或時(shí)間緊迫時(shí)可根據(jù)已有的實(shí)測數(shù)據(jù)建立水位流量相關(guān)關(guān)系，通過水位流量關(guān)系示意推求瞬時(shí)流量，滿足汛期水情隨測隨報(bào)的要求[8]。介于水情報(bào)汛對及時(shí)性的要求，在水位流量關(guān)系曲線率定時(shí)，其原始數(shù)據(jù)只需錄入實(shí)測水位與流量值，測驗(yàn)系統(tǒng)通過三項(xiàng)式回歸方程對各數(shù)值進(jìn)行擬合，并動(dòng)態(tài)繪制成光滑的曲線[9]。此時(shí)即可通過關(guān)系曲線圖直觀的查閱不同水位下的流量數(shù)值，支持高水延長(見圖8)。

圖8 靈臺(tái)水文站水位流量關(guān)系曲線示意

4 應(yīng)用成果

參與完成靈臺(tái)站年度SEBA PULS型雷達(dá)水位計(jì)水位資料整編；完成靈臺(tái)站及附屬雨量站(坷臺(tái)、百里、楊新莊、杜家川、李家灣)JH-2型雨量計(jì)降水量資料整編；參與完成年度流量測驗(yàn)與整編任務(wù)，共測得流量成果78份；其中流速儀法流量成果43份，浮標(biāo)法流量成果35份；參與完成年度定點(diǎn)洪水調(diào)查任務(wù)，共測得澗河、蒲河、南河洪水調(diào)查成果3份；參與完成年度報(bào)汛任務(wù)，尤其是在應(yīng)對“8.17”特大暴雨洪水中該系統(tǒng)發(fā)揮了重要作用，達(dá)到了水情報(bào)汛期的及時(shí)性與精確度雙層要求。

5 結(jié)語

靈臺(tái)水文站于2020年度使用《甘肅省水文測驗(yàn)系統(tǒng)》，參與測站日常測驗(yàn)與水文資料在站整編工作，在降水量、水位、流量、報(bào)汛等測驗(yàn)項(xiàng)目中成效顯著，較好的完成了年度測驗(yàn)任務(wù)。測驗(yàn)系統(tǒng)操作簡單、界面友好、功能完善，滿足基層水文工作需求，實(shí)現(xiàn)了水文資料在站整編的數(shù)字化與自動(dòng)化，提高了在站整編工作效率和質(zhì)量。各項(xiàng)計(jì)算結(jié)果符合水文行業(yè)“四舍六入，逢五奇進(jìn)偶舍”的運(yùn)算規(guī)則[10]，符合《水文年鑒匯編刊印規(guī)范》(SL 460-2009)規(guī)定要求，可在甘肅省部分水文站及省外相關(guān)單位推廣使用。