唐海華，李 琪，黃瓅瑤，周 超，羅 斌

(長(zhǎng)江勘測(cè)規(guī)劃設(shè)計(jì)研究有限責(zé)任公司，武漢 430010)

1 研究背景

流域洪水調(diào)控計(jì)算本質(zhì)上是深度運(yùn)用計(jì)算機(jī)技術(shù)，模擬流域洪水在其防洪體系下的傳播演進(jìn)與調(diào)度控制過(guò)程。我國(guó)流域的防洪體系主要包括水庫(kù)、河段、測(cè)站、堤防、蓄滯洪區(qū)、水閘、泵站等水利對(duì)象[1-2]，重點(diǎn)涉及河道洪水一維演進(jìn)、水庫(kù)洪水調(diào)度、水庫(kù)回水模擬、水庫(kù)潰壩模擬、堤防潰決模擬、蓄滯洪區(qū)分洪計(jì)算、泄洪閘門分配計(jì)算、洪水二維淹沒(méi)模擬、洪災(zāi)損失計(jì)算和洪水風(fēng)險(xiǎn)評(píng)估等業(yè)務(wù)計(jì)算[3-7]。充分運(yùn)用各類專業(yè)模型對(duì)水利工程對(duì)象群開展聯(lián)合調(diào)度控制，對(duì)提升流域防洪能力、降低洪水淹沒(méi)風(fēng)險(xiǎn)、減小洪水災(zāi)害損失具有重要支撐作用。

目前，流域各層級(jí)防汛調(diào)度部門主要依靠防汛調(diào)度決策指揮系統(tǒng)開展洪水調(diào)控工作。此類系統(tǒng)主要基于固定的范圍、對(duì)象、流程和模型進(jìn)行設(shè)計(jì)建設(shè)[8-9]，其應(yīng)用功能模塊、防洪控制對(duì)象、調(diào)度計(jì)算任務(wù)、業(yè)務(wù)數(shù)據(jù)流程、設(shè)施控制方式等大多按定制化方式開發(fā)實(shí)現(xiàn)，各類專業(yè)模型分散在不同的軟件功能中，原本連續(xù)的洪水調(diào)控任務(wù)被模塊化分割，缺乏整體性、系統(tǒng)性和一致性，數(shù)據(jù)銜接與信息處理需大量人工干預(yù)，流程繁瑣；同時(shí)，流域及各類水利工程的洪水調(diào)度方案通常需要定期修編，對(duì)應(yīng)的系統(tǒng)功能及專業(yè)模型必須配套開發(fā)調(diào)整才能滿足新方案要求，工作量大，效率較低；此外，當(dāng)遭遇超標(biāo)準(zhǔn)洪水時(shí)，由于超標(biāo)準(zhǔn)洪水具有偶發(fā)性、不確定性、復(fù)雜性和特殊性，調(diào)度人員通常需要根據(jù)洪水的演變態(tài)勢(shì)在線開展多方案調(diào)控計(jì)算模擬，才能科學(xué)下達(dá)防汛決策和調(diào)控指令，計(jì)算范圍、對(duì)象和流程會(huì)隨著洪水影響程度不同而動(dòng)態(tài)變化，專業(yè)模型、參數(shù)類型與數(shù)據(jù)規(guī)模也較為龐大，無(wú)法提前預(yù)置計(jì)算對(duì)象和業(yè)務(wù)流程[10]。

綜上，按定制化方式建設(shè)防汛調(diào)度決策支持系統(tǒng)很難支撐變化場(chǎng)景下的不確定性洪水調(diào)控需求，應(yīng)對(duì)能力明顯不足。組態(tài)化、組件化和流程化技術(shù)是應(yīng)對(duì)系統(tǒng)研發(fā)不確定性需求的重要手段。馬振宇等[11]將組態(tài)設(shè)計(jì)思想引入管理信息系統(tǒng)結(jié)構(gòu)及其業(yè)務(wù)工作流，使系統(tǒng)可有效適應(yīng)業(yè)務(wù)流程重組，大幅減少了需求的不確定性風(fēng)險(xiǎn)；唐海華等[12]充分應(yīng)用組件化技術(shù)提出了水庫(kù)群調(diào)度的多模型一體化集成架構(gòu)，實(shí)現(xiàn)了模型算法的通用集成；萬(wàn)書鵬等[13]基于組態(tài)化和流程化模式開展了調(diào)控系統(tǒng)的微服務(wù)編排技術(shù)研究，實(shí)現(xiàn)了系統(tǒng)業(yè)務(wù)場(chǎng)景的靈活組裝，滿足了各級(jí)調(diào)度控制中心業(yè)務(wù)場(chǎng)景的功能組合需求。本文充分借鑒上述思想，深入分析其與洪水調(diào)控計(jì)算的融合應(yīng)用方法，重點(diǎn)研究不同洪水調(diào)控場(chǎng)景下差異化計(jì)算需求的一體化動(dòng)態(tài)耦合技術(shù)。

2 關(guān)鍵思路

為解決洪水調(diào)控模擬過(guò)程的定制化和模塊分割問(wèn)題，實(shí)現(xiàn)流域洪水的河道傳播演進(jìn)、水庫(kù)調(diào)度決策、庫(kù)區(qū)回水線推求、泄洪控制計(jì)劃擬定、堤壩潰決與蓄滯洪區(qū)分洪、洪水淹沒(méi)分析、洪災(zāi)損失統(tǒng)計(jì)及洪水風(fēng)險(xiǎn)等級(jí)評(píng)估等全過(guò)程一體化模擬，并支持按需自定義組織不同規(guī)模、不同類型的洪水調(diào)控計(jì)算任務(wù)，本文將組件化、組態(tài)化和流程引擎技術(shù)引入流域洪水調(diào)度控制領(lǐng)域。

組件化技術(shù)是模塊化思想的演進(jìn)和延伸，主要通過(guò)深層次解耦與隔離不需要關(guān)注的部分來(lái)強(qiáng)化模塊角色的可轉(zhuǎn)換性；組態(tài)化技術(shù)的精髓在于可通過(guò)“搭積木”的方式來(lái)按需配置自己想實(shí)現(xiàn)的不同模擬計(jì)算功能，而不依賴編寫計(jì)算機(jī)程序代碼；流程引擎技術(shù)的核心則是將不同的目標(biāo)任務(wù)分解為節(jié)點(diǎn)、流向和流程，然后根據(jù)不同的角色、分工、條件來(lái)決定信息流和業(yè)務(wù)流的傳遞方向與轉(zhuǎn)換邏輯，從而實(shí)現(xiàn)多項(xiàng)任務(wù)的靈活銜接與自動(dòng)關(guān)聯(lián)。

將以上技術(shù)與洪水調(diào)控計(jì)算的水利對(duì)象、業(yè)務(wù)邏輯和模型算法深度融合，即可構(gòu)建一套以“流程節(jié)點(diǎn)組件分解”與“調(diào)控計(jì)算組態(tài)耦合”相銜接的技術(shù)方案，按照“先分后合”的方式實(shí)現(xiàn)各類流域洪水調(diào)控任務(wù)的按需組態(tài)耦合，以動(dòng)態(tài)配置搭建替代程序定制開發(fā)，有效應(yīng)對(duì)不同洪水調(diào)控場(chǎng)景下的不確定性計(jì)算需求。

3 技術(shù)方案

3.1 流程節(jié)點(diǎn)組件分解

針對(duì)流域防洪體系各類水利對(duì)象、模型算法和方案參數(shù)，按對(duì)象庫(kù)、規(guī)則庫(kù)和節(jié)點(diǎn)庫(kù)3類進(jìn)行分解和管理，為調(diào)控計(jì)算組態(tài)耦合提供關(guān)鍵要素支撐，如圖1所示。

圖1 流程節(jié)點(diǎn)組件分解示意圖Fig.1 Schematic diagram of process node decomposition

(1)防洪體系數(shù)字化建模：梳理流域內(nèi)所有防洪體系資料，按照不同計(jì)算對(duì)象類型，分別進(jìn)行數(shù)字化建模，分類定義并量化當(dāng)前流域內(nèi)所有對(duì)象的基礎(chǔ)屬性、設(shè)計(jì)參數(shù)和特征指標(biāo)，構(gòu)建防洪體系對(duì)象庫(kù)。

(2)模型參數(shù)規(guī)范化定義：針對(duì)洪水調(diào)控涉及的各類計(jì)算模型，按對(duì)象庫(kù)分類方式分別梳理輸入輸出參數(shù)，并剔除重復(fù)變量，統(tǒng)一制定數(shù)據(jù)結(jié)構(gòu)標(biāo)準(zhǔn)，構(gòu)建所有模型的參數(shù)庫(kù)。

(3)模型算法組件化封裝：將洪水調(diào)控涉及的各類計(jì)算模型與具體水利對(duì)象解耦，統(tǒng)一采用參數(shù)庫(kù)中的抽象定義作為輸入輸出接口，按組件化方式實(shí)現(xiàn)算法開發(fā)，提升模型自身的通用性和移植性，構(gòu)建支撐洪水調(diào)控計(jì)算的組件庫(kù)。

(4)模型組件節(jié)點(diǎn)化構(gòu)建：針對(duì)所有模型算法組件，根據(jù)其接口需求逐一與參數(shù)庫(kù)關(guān)聯(lián)，構(gòu)成具有“輸入-計(jì)算-輸出”標(biāo)準(zhǔn)化結(jié)構(gòu)的流程節(jié)點(diǎn)，最終將組件庫(kù)和參數(shù)庫(kù)全部封裝為面向調(diào)控任務(wù)的節(jié)點(diǎn)庫(kù)。

(5)對(duì)象耦合校驗(yàn)邏輯制定：由于節(jié)點(diǎn)庫(kù)中每個(gè)節(jié)點(diǎn)已全部與水利對(duì)象解耦，當(dāng)創(chuàng)建流程節(jié)點(diǎn)開展實(shí)例計(jì)算時(shí)，必須與對(duì)象庫(kù)中的某一類或幾類水利對(duì)象進(jìn)行耦合。因此，針對(duì)任一類水利對(duì)象，必須判斷節(jié)點(diǎn)與對(duì)象耦合的有效性，只有節(jié)點(diǎn)內(nèi)封裝的模型組件接口與該類對(duì)象的屬性參數(shù)存在關(guān)聯(lián)性才為有效耦合。

定義對(duì)象庫(kù)的類型集合OT,類型總數(shù)m；節(jié)點(diǎn)庫(kù)的節(jié)點(diǎn)集合ND，節(jié)點(diǎn)總數(shù)n。可制定節(jié)點(diǎn)對(duì)象耦合有效性集合C為

(1)

式中C(i,j)為0-1變量，表示節(jié)點(diǎn)NDi與對(duì)象OTj的耦合有效性，且i∈(1,2，…，n),j∈(1,2，…，m)。

針對(duì)集合C中C(i,j)=1的元素，構(gòu)建對(duì)應(yīng)的節(jié)點(diǎn)對(duì)象適配集合A為

(2)

式中A(i,j) 表示節(jié)點(diǎn)NDi與對(duì)象OTj的接口參數(shù)適配關(guān)系，具體表達(dá)式根據(jù)不同節(jié)點(diǎn)和對(duì)象類型分別定義。若C(i,j)=0，則對(duì)應(yīng)A(i,j)無(wú)效，按置空處理。

(6)節(jié)點(diǎn)流向校驗(yàn)邏輯制定：節(jié)點(diǎn)庫(kù)中的每個(gè)節(jié)點(diǎn)在開展組合計(jì)算流程搭建時(shí)，需與其他節(jié)點(diǎn)進(jìn)行連接。若將節(jié)點(diǎn)A連接至節(jié)點(diǎn)B時(shí)，則產(chǎn)生節(jié)點(diǎn)流向A-B，其中A為上級(jí)節(jié)點(diǎn)，B為下級(jí)節(jié)點(diǎn)。此時(shí)，必須判斷節(jié)點(diǎn)流向的有效性，其中上級(jí)節(jié)點(diǎn)的輸出參數(shù)中至少存在一項(xiàng)與下級(jí)節(jié)點(diǎn)的輸入?yún)?shù)關(guān)聯(lián)，節(jié)點(diǎn)流向才有效，否則無(wú)效。

構(gòu)建節(jié)點(diǎn)流向有效性集合V為

(3)

式中V(i,j)為0-1變量，表示節(jié)點(diǎn)流向NDi-NDj的有效性，且i∈(1,2，…，n),j∈(1,2，…，n)。

針對(duì)集合V中V(i,j)=1的元素，構(gòu)建對(duì)應(yīng)的節(jié)點(diǎn)參數(shù)關(guān)聯(lián)集合L為

(4)

式中L(i,j)表示節(jié)點(diǎn)流向NDi-NDj的參數(shù)關(guān)聯(lián)邏輯，具體表達(dá)式根據(jù)不同節(jié)點(diǎn)類型分別定義。若V(i,j)=0，則對(duì)應(yīng)L(i,j)無(wú)效，按置空處理。

式(1)—式(4)共同構(gòu)成了所有節(jié)點(diǎn)、對(duì)象和流向的邏輯校驗(yàn)規(guī)則庫(kù)。

圖2 調(diào)控計(jì)算組態(tài)耦合流程示意圖Fig.2 Flowchart of configuration coupling for control calculation

3.2 調(diào)控計(jì)算組態(tài)耦合

根據(jù)對(duì)象庫(kù)、規(guī)則庫(kù)和節(jié)點(diǎn)庫(kù)成果，針對(duì)不同調(diào)控計(jì)算任務(wù)的業(yè)務(wù)邏輯，依次組織對(duì)應(yīng)的對(duì)象、規(guī)則和節(jié)點(diǎn)，即可按流程化方式實(shí)現(xiàn)該任務(wù)的組態(tài)耦合，如圖2所示。具體過(guò)程為：?jiǎn)?dòng)流程，從節(jié)點(diǎn)庫(kù)中選擇一個(gè)節(jié)點(diǎn)作為流程起點(diǎn)；為起始節(jié)點(diǎn)耦合計(jì)算對(duì)象，并根據(jù)規(guī)則庫(kù)校驗(yàn)當(dāng)前節(jié)點(diǎn)與對(duì)象耦合是否有效，若無(wú)效，則重選對(duì)象類型；根據(jù)組合需求創(chuàng)建下級(jí)節(jié)點(diǎn)，并根據(jù)規(guī)則庫(kù)校驗(yàn)上下級(jí)節(jié)點(diǎn)流向是否有效，若無(wú)效，則重選下級(jí)節(jié)點(diǎn)；為下級(jí)節(jié)點(diǎn)耦合計(jì)算對(duì)象，并校驗(yàn)有效性；根據(jù)組合計(jì)算需求，創(chuàng)建新的下級(jí)流程節(jié)點(diǎn)并完成對(duì)應(yīng)的節(jié)點(diǎn)流向和對(duì)象耦合校驗(yàn)，若下級(jí)需同時(shí)開展多項(xiàng)任務(wù)，則創(chuàng)建多個(gè)并聯(lián)子節(jié)點(diǎn)，依次循環(huán)，直至節(jié)點(diǎn)創(chuàng)建完畢；將上述所有流程節(jié)點(diǎn)的輸入信息進(jìn)行歸納集中，各上、下級(jí)節(jié)點(diǎn)之間存在關(guān)聯(lián)的參數(shù)，根據(jù)式(4)的邏輯規(guī)則確定，其余則由用戶輸入確認(rèn)；執(zhí)行計(jì)算任務(wù)，按流程順序逐節(jié)點(diǎn)完成計(jì)算和信息傳遞；輸出所有流程節(jié)點(diǎn)的對(duì)象計(jì)算結(jié)果，流程結(jié)束。

4 實(shí)例應(yīng)用

本文技術(shù)已成功應(yīng)用于長(zhǎng)江流域控制性水工程防災(zāi)聯(lián)合調(diào)度示范系統(tǒng)，實(shí)現(xiàn)了水庫(kù)群防洪調(diào)度、河道洪水演進(jìn)、庫(kù)區(qū)回水分析等一體化流域洪水調(diào)控計(jì)算。以長(zhǎng)江中上游為例，其防洪體系主要涵蓋：梨園、阿海、三峽等16座水庫(kù)；梨園—阿海、阿?！鸢矘颉⑷龒{—宜昌等19個(gè)河段；李莊、朱沱、宜昌等8個(gè)測(cè)站，以及三峽庫(kù)區(qū)140個(gè)河道斷面對(duì)象。防洪體系拓?fù)浣Y(jié)構(gòu)概化如圖3所示。

圖3 防洪體系拓?fù)浣Y(jié)構(gòu)示意圖Fig.3 Topological structure of flood control system

首先，建立水庫(kù)、河段、測(cè)站、河道斷面等4種對(duì)象類型，分別記為RV、RR、ST、RS；然后針對(duì)各類型的實(shí)體對(duì)象進(jìn)行數(shù)字化建模，構(gòu)成對(duì)象庫(kù)。

針對(duì)4種對(duì)象分類定義輸入輸出及屬性參數(shù)，構(gòu)成參數(shù)庫(kù)。其中水庫(kù)參數(shù)包括入庫(kù)流量、出庫(kù)流量、壩上水位、水位庫(kù)容曲線、泄流能力曲線、下游水位流量曲線、最高水位約束、最低水位約束、最大出庫(kù)約束、最小出庫(kù)約束、下游安全泄量、回水?dāng)嗝婕?；河段參?shù)包括上邊界編碼、下邊界編碼、區(qū)間匯入站編碼集、河道斷面編碼集等；測(cè)站參數(shù)主要包括雨量、水位、流量、關(guān)聯(lián)斷面編碼、警戒水位、保證水位、警戒流量、保證流量、水位流量曲線等；河道斷面參數(shù)包括起點(diǎn)距-高程數(shù)組、下斷面間距、糙率、水位、流量、移民水位、土地征用水位等。

將河道洪水一維演進(jìn)、水庫(kù)洪水調(diào)度、水庫(kù)回水模擬3個(gè)模型按統(tǒng)一接口梳理輸入輸出參數(shù)，然后采用“輸入接口+算法+輸出接口”的格式進(jìn)行通用化封裝，形成模型組件庫(kù)。

將各模型組件的輸入、輸出接口逐一與參數(shù)庫(kù)關(guān)聯(lián)，可封裝為河道洪水一維演進(jìn)、水庫(kù)洪水調(diào)度、水庫(kù)回水模擬共3個(gè)流程節(jié)點(diǎn)，分別記為FR、RD、BW，形成節(jié)點(diǎn)庫(kù)。

令對(duì)象庫(kù)中4種對(duì)象類型集合為OT，節(jié)點(diǎn)庫(kù)中3個(gè)節(jié)點(diǎn)集合為ND，則：

OT={RV,RR,ST,RS} ；

(5)

ND={FR,RD,BW} 。

(6)

將式(5)—式(6)代入式(1)，得4×3的節(jié)點(diǎn)與對(duì)象類型耦合有效性集合C，即

(7)

式(7)中C(i,j)=1的耦合項(xiàng)繼續(xù)構(gòu)建節(jié)點(diǎn)NDi與對(duì)象OTj的接口參數(shù)適配關(guān)系集合A，用于確定對(duì)象庫(kù)中參與節(jié)點(diǎn)計(jì)算的實(shí)例對(duì)象編碼，以及參數(shù)庫(kù)中與節(jié)點(diǎn)計(jì)算相關(guān)聯(lián)的對(duì)象參數(shù)。

將式(6)代入式(3)，得3×3的節(jié)點(diǎn)流向有效性集合為

(8)

式(8)中V(i,j)=1的節(jié)點(diǎn)繼續(xù)構(gòu)建NDi對(duì)應(yīng)NDj的參數(shù)關(guān)聯(lián)集合L，用于確定NDi的輸出參數(shù)與NDj的輸入?yún)?shù)之間的關(guān)聯(lián)方式。

至此，本例所有對(duì)象庫(kù)、節(jié)點(diǎn)庫(kù)和規(guī)則庫(kù)全部創(chuàng)建完成，根據(jù)計(jì)算任務(wù)分解需求，其搭建步驟為：

(1)啟動(dòng)流程，從節(jié)點(diǎn)庫(kù)中選擇RD為流程起點(diǎn)，根據(jù)式(7)得有效耦合對(duì)象類型為RV，選擇龍頭水庫(kù)梨園。

(2)根據(jù)式(8)得RD的有效節(jié)點(diǎn)流向?yàn)镕R和BW。創(chuàng)建二級(jí)節(jié)點(diǎn)FR，其有效耦合對(duì)象類型為RR，選擇梨園—阿海河段。

(3)按照上述方式，根據(jù)圖3的拓?fù)潢P(guān)系依次創(chuàng)建下級(jí)節(jié)點(diǎn)并耦合計(jì)算對(duì)象。其中三峽水庫(kù)調(diào)度完成后需同時(shí)開展回水計(jì)算和下游河道演算，因此包含BW和FR兩個(gè)子節(jié)點(diǎn)，分別耦合三峽水庫(kù)和三峽—宜昌河段。

(4)根據(jù)式(2)歸集上述各級(jí)流程節(jié)點(diǎn)的輸入信息。其中RD節(jié)點(diǎn)包括各水庫(kù)的入庫(kù)流量、壩上水位、水位庫(kù)容曲線、泄流能力曲線、最高水位約束、最低水位約束、最大出庫(kù)約束、最小出庫(kù)約束；RR節(jié)點(diǎn)包括各河段的上邊界流量、下邊界水位、區(qū)間匯入流量，以及河道斷面集合的水位、起點(diǎn)距-高程數(shù)組、斷面間距和糙率；BW節(jié)點(diǎn)包括各水庫(kù)的入庫(kù)流量、出庫(kù)流量、壩上水位、下游水位流量曲線，以及回水?dāng)嗝婕系乃?、起點(diǎn)距-高程數(shù)組、斷面間距和糙率。

(5)根據(jù)式(4)將各級(jí)節(jié)點(diǎn)的上級(jí)輸出參數(shù)與下級(jí)輸入?yún)?shù)進(jìn)行關(guān)聯(lián)，如梨園—阿海河段的下邊界流量可直接關(guān)聯(lián)阿海水庫(kù)的入庫(kù)流量、阿海水庫(kù)的出庫(kù)流量則直接關(guān)聯(lián)阿海—金安橋河段的上邊界流量，其余節(jié)點(diǎn)依此類推。對(duì)于無(wú)法關(guān)聯(lián)的參數(shù)，直接從數(shù)據(jù)庫(kù)自動(dòng)提取，若需用戶交互則開放至界面輸入。

(6)執(zhí)行流程任務(wù)，各級(jí)節(jié)點(diǎn)按流程順序依次計(jì)算，并自動(dòng)往下級(jí)節(jié)點(diǎn)流轉(zhuǎn)信息，計(jì)算完成后，輸出所有對(duì)象計(jì)算結(jié)果，流程結(jié)束。

5 結(jié) 語(yǔ)

根據(jù)本文提出的流域洪水調(diào)控計(jì)算組態(tài)耦合技術(shù)成果，開發(fā)配套軟件操作界面，可有效解決不確定性洪水差異化調(diào)控需求的快速響應(yīng)問(wèn)題，支撐防汛決策支持系統(tǒng)的敏捷開發(fā)搭建。若任務(wù)需求發(fā)生變化，根據(jù)新的業(yè)務(wù)邏輯靈活調(diào)整流程節(jié)點(diǎn)和計(jì)算對(duì)象，即可完成動(dòng)態(tài)適應(yīng)。目前，本文技術(shù)僅在流域洪水調(diào)控領(lǐng)域進(jìn)行了應(yīng)用驗(yàn)證，未來(lái)還可進(jìn)一步凝練概化，增強(qiáng)通用性與普適性，拓展應(yīng)用于水資源配置、引調(diào)水工程水量調(diào)度、水生態(tài)調(diào)度、水環(huán)境水質(zhì)模擬、水污染擴(kuò)散模擬、泥沙調(diào)度模擬、水電站群發(fā)電調(diào)度、電網(wǎng)水風(fēng)光互補(bǔ)調(diào)節(jié)等類似業(yè)務(wù)，在水利、水電行業(yè)具有廣闊的應(yīng)用前景。