專利名稱:一種基于水功能區(qū)的水量水質(zhì)調(diào)控方法
技術(shù)領(lǐng)域:
本發(fā)明涉及一種基于水功能區(qū)的水量水質(zhì)調(diào)控方法,是一種水資源調(diào)控的方法,是一種對數(shù)千平方千米及以上的區(qū)域或流域的河流水功能區(qū)的水量、水質(zhì)進(jìn)行分析和水環(huán)境控制的方法。
背景技術(shù):
對生態(tài)環(huán)境的分析和調(diào)控需要將自然環(huán)境和人類活動的諸多因素進(jìn)行整合,剔除影響較小的因素,確定關(guān)鍵因素,建立與實際河流狀態(tài)最近似的模型。并在模型的應(yīng)用中不斷的擴(kuò)充和修正模型,加入更多的因素,使模型越發(fā)接近實際的河流和水環(huán)境狀態(tài)。早期的流域級的大范圍水量水質(zhì)的分析由于受到當(dāng)時科研水平和手段的限制,只能尋找數(shù)個關(guān)鍵性的因素建立模型。這樣的模型只能近似的模擬水量與污染物入河量狀況,與實際相差較大?,F(xiàn)有的水量水質(zhì)調(diào)控方法對于污染控制的分析缺乏對用水導(dǎo)致水量和污染負(fù)荷排放對水功能區(qū)達(dá)標(biāo)狀況的雙重影響因素的考慮,難以從技術(shù)角度分析合理的水量水質(zhì)聯(lián)合調(diào)控措施,以實現(xiàn)水功能區(qū)達(dá)標(biāo)要求?,F(xiàn)有河流納污能力計算方法主要是以河流設(shè)計流量分析計算水功能區(qū)的納污能力,通過對水文歷史資料分析,通過頻率分析得出給定的偏枯保證率(通常為75%或90% )條件下的河流斷面流量,以該流量為基礎(chǔ),考慮河流自凈能力設(shè)計分析可以容許水功能區(qū)不超標(biāo)的污染負(fù)荷入河量。沒有考慮到人工用水和工程調(diào)度(人工的水源分配)與自然水循環(huán)之間的動態(tài)關(guān)系,不能反映人工取水、水利工程調(diào)控等因素改變徑流過程進(jìn)而影響納污能力的普遍性實際情況。
發(fā)明內(nèi)容
為了克服現(xiàn)有技術(shù)的問題,本發(fā)明提出了一種基于水功能區(qū)的水量水質(zhì)調(diào)控方法,所述的方法綜合分析用水與工程調(diào)度對河道流量過程的影響,建立調(diào)控模型,并提出了基于用水調(diào)控的水功能區(qū)動態(tài)納污能力的分析計算方式,為水功能區(qū)污染負(fù)荷入河控制量分析提供基礎(chǔ),實現(xiàn)流域水質(zhì)水量總量控制技術(shù)。本發(fā)明所述的方法對水量調(diào)控和水質(zhì)目標(biāo)控制分別按照水資源分區(qū)和河流水功能區(qū)進(jìn)行控制,并建立二者之間的水量傳輸排放關(guān)系。本發(fā)明的目的是這樣實現(xiàn)的一種基于水功能區(qū)的水量水質(zhì)調(diào)控方法,所述方法的步驟如下構(gòu)建水量調(diào)控模擬模型的步驟用于以流域水資源分區(qū)和行政分區(qū)嵌套形成水功能區(qū)的計算單元,并將水功能區(qū)的主要控制工程節(jié)點和水功能區(qū)控制斷面作為基本信息模擬單元;分析所述水功能區(qū)范圍內(nèi)的水量平衡信息,建立計算單元與水功能區(qū)地表徑流匯流與用水退水關(guān)系,計算用水影響下的水功能區(qū)控制斷面過流;所述的水功能區(qū)斷面的流量計算公式如下Q = W±+W*-U+S式中Q-水功能區(qū)控制斷面過流量,
W上-上游來水量,W本-水功能區(qū)所在單元地表徑流量,U-水功能區(qū)所在單元地表取用水量,S-水功能區(qū)所在單元用水后退入水功能區(qū)斷面水量;計算水功能區(qū)控制斷面流量的步驟用于使用水量調(diào)控模擬模型計算水功能區(qū)的斷面流量;計算入河污染負(fù)荷量的步驟用于對入河的污染源進(jìn)行量化的估算,包括對工業(yè)點源產(chǎn)生量估算和城鎮(zhèn)生活點源產(chǎn)生量估算;基于水功能區(qū)的納污能力計算的步驟用于根據(jù)水量調(diào)控模擬模型的流量計算結(jié)果,采用一維水質(zhì)模擬公式計算用水影響條件下的水功能區(qū)納污能力
M = 31. 5 (Cs-C0 exp (-kL/u)) exp (kx/u) Qr,式中M_污染物納污能力(t/a),Cs-水功能區(qū)下斷面水質(zhì)目標(biāo)控制因子濃度(mg/L),Ctl-水功能區(qū)上斷面水質(zhì)目標(biāo)控制因子濃度(mg/L),k-污染物綜合衰減系數(shù)(1/s),L-水功能區(qū)長度(m),X-排污口下游斷面距水功能區(qū)下斷面縱向距離(m),u-設(shè)計流量下的平均流速(m/s),Qr-設(shè)計流量(m3/s);分析水功能區(qū)達(dá)標(biāo)狀況的步驟用于使用水功能區(qū)達(dá)標(biāo)分析模型,通過污染負(fù)荷產(chǎn)生入河量與水功能區(qū)納污能力對比分析水功能區(qū)達(dá)標(biāo)狀況,如果若水功能區(qū)達(dá)標(biāo)率不能實現(xiàn),則回到“構(gòu)建水量調(diào)控模擬模型的步驟”調(diào)整數(shù)據(jù)量,重新進(jìn)行計算,直至實現(xiàn)水功能區(qū)達(dá)標(biāo)率,則結(jié)束過程。本發(fā)明產(chǎn)生的有益效果是本發(fā)明通過構(gòu)建流域級水功能區(qū)的水量調(diào)控模擬模型,考慮人工用水、工程調(diào)度與自然水循環(huán)關(guān)系,同時構(gòu)建水功能區(qū)納污能力分析模型,使用納污能力分析模型以及水量調(diào)控模擬模型為依據(jù)分析計算水功能區(qū)的納污能力,根據(jù)分析計算結(jié)果對流域的水量和水質(zhì)進(jìn)行調(diào)控。與現(xiàn)有技術(shù)相比,本發(fā)明充分考慮人工用水、工程調(diào)度與自然水循環(huán)三方面的因素,對流域水質(zhì)的納污能力分析更加準(zhǔn)確,使水量調(diào)控更加合理。
下面結(jié)合附圖和實施例對本發(fā)明作進(jìn)一步說明。圖I是本發(fā)明的實施例一所述方法的流程圖;圖2是本發(fā)明實施例一所述方法中水量調(diào)控模擬模型的水源用戶配置關(guān)系示意圖;圖3是本發(fā)明實施例一所述方法的水量調(diào)控模擬模型的計算流程示意圖;圖4是計算單元與水功能區(qū)關(guān)系示意圖;圖5是本發(fā)明的實施例二所述水量調(diào)控模擬模型的計算基本流程圖。
具體實施例方式實施例一本實施例是一種基于水功能區(qū)的水量水質(zhì)調(diào)控方法。本實施例采用了一種水資源分區(qū)和水功能區(qū)嵌套計算的新形式,實現(xiàn)了水量調(diào)控模擬和河流納污能力計算的融合。本實施例的關(guān)鍵在于一是整體水量水質(zhì)調(diào)控的技術(shù)設(shè)計;二是基于水功能區(qū)的水量調(diào)控模擬模型;三是基于用水的動態(tài)納污能力計算;四是水功能區(qū)達(dá)標(biāo)狀況分析模型。其中,基于水功能區(qū)的水量調(diào)控模擬模型計算是本實施例的重點所在。本實施例所述方法中涉及水資源分區(qū)、水功能區(qū)及計算單元等不同性質(zhì)而又互相關(guān)聯(lián)的信息承載體。其中,計算單元為流域劃分中的水資源分區(qū)與行政區(qū)單元嵌套形成的區(qū)域,承載區(qū)域水資源量以及供用水等信息;水功能區(qū)為表征具有水質(zhì)目標(biāo)要求的河流分段劃分。 所述的水資源分區(qū)是指針對水資源規(guī)劃提出的流域劃分,水資源分區(qū)是一種類似行政區(qū)域的面狀結(jié)構(gòu)。水功能區(qū)是以河流為劃分對象,科學(xué)合理地在相應(yīng)水域劃定具有特定功能、滿足水資源開發(fā)利用和保護(hù)要求的河流分段劃分。水功能區(qū)是以河流為主的線狀信息結(jié)構(gòu)。水量調(diào)控模擬模型中的計算單元是水資源分區(qū)和行政分區(qū)通過地理位置上相互嵌套關(guān)系形成的區(qū)域,模擬中每個計算單元包括不同用戶的用水需求、污水處理水平、小型工程調(diào)蓄與供水能力等信息。本實施例所述方法根據(jù)對于現(xiàn)有納污能力計算局限性的分析,對水質(zhì)水量的總量調(diào)控分析涉及對水量過程、水質(zhì)模擬過程的聯(lián)合分析。上述各部分工作通過方案設(shè)置和多層次反饋實現(xiàn)對水質(zhì)水量聯(lián)合調(diào)控的模擬分析。本實施例所述方法包括水量調(diào)控模擬模型、入河污染負(fù)荷預(yù)測模型、納污能力分析模型和水功能區(qū)達(dá)標(biāo)分析模型。模型間的相互關(guān)系為水量調(diào)控模擬模型計算用水與水功能區(qū)控制斷面流量;入河污染負(fù)荷預(yù)測模型分析入河污染負(fù)荷量;納污能力分析模型根據(jù)水量模擬調(diào)控的流量結(jié)果計算特定用水情景下水功能區(qū)納污能力,通過納污能力分析模型對比動態(tài)水量條件下的水功能區(qū)納污能力;水功能區(qū)達(dá)標(biāo)分析模型通過污染負(fù)荷產(chǎn)生入河量與水功能區(qū)納污能力對比分析水功能區(qū)達(dá)標(biāo)狀況。在上述模型順序計算得出水功能區(qū)達(dá)標(biāo)狀況后,若水功能區(qū)達(dá)標(biāo)率不能實現(xiàn),可以進(jìn)一步根據(jù)輸入信息對用水量、工程調(diào)度方案、污水處理和污染控制措施進(jìn)行反饋調(diào)整,從源用水、治污、調(diào)度等環(huán)節(jié)實現(xiàn)水量水質(zhì)的聯(lián)合調(diào)控,提出總量控制目標(biāo),達(dá)到污染控制目的。本實施例所述方法的步驟如下,如圖I所示構(gòu)建水量調(diào)控模擬模型的步驟用于以流域水資源分區(qū)和行政分區(qū)嵌套形成水功能區(qū)的計算單元,并將水功能區(qū)的主要控制工程節(jié)點和水功能區(qū)控制斷面作為基本信息模擬單元;分析所述水功能區(qū)范圍內(nèi)的水量平衡信息,建立計算單元與水功能區(qū)地表徑流匯流與用水退水關(guān)系,計算用水影響下的水功能區(qū)控制斷面過流;所述的水功能區(qū)斷面的流量計算公式如下Q = W±+W*-U+S式中Q-水功能區(qū)控制斷面過流量,W上-上游來水量,
W$-水功能區(qū)所在單元地表徑流量,U-水功能區(qū)所在單元地表取用水量,S-水功能區(qū)所在 單元用水后退入水功能區(qū)斷面水量。水量調(diào)控模擬模型現(xiàn)有關(guān)于水量調(diào)控模擬模型,在國際上開發(fā)和應(yīng)用比較成熟的有SWAT模型。針對地下水有M0DFL0W模型。SWAT模型是一種分布式水文模型,同時考慮了用水對水文過程的影響。郭正鑫(2009)就北京溫榆河為例,使用工具ARCGIS處理獲得的數(shù)據(jù),建立起能夠描述河網(wǎng)閘壩水系統(tǒng)的SWAT模型來實現(xiàn)對流域水量水質(zhì)的聯(lián)合調(diào)控。但SWAT模型的應(yīng)用需要以分布式水文過程為基礎(chǔ),需要大量基礎(chǔ)信息資料,包括對區(qū)域水文氣象、地形地貌、土地利用等信息數(shù)據(jù)的處理輸入,同時還需要應(yīng)用一定長度歷史水文數(shù)據(jù)對模型進(jìn)行驗證和參數(shù)率定??紤]與現(xiàn)有的供用水統(tǒng)計信息資料相一致,本實施例提出的水量調(diào)控模擬模型不采用分布式模型方式,主要采用概念化處理系統(tǒng)水量過程的方式實現(xiàn)模擬,同時采用模塊化分割方法實現(xiàn)不同類別水量的模擬。水量調(diào)控模擬模型主要針對用水、河道水量調(diào)度過程進(jìn)行模擬,分析不同用水和調(diào)度方案對水平衡過程的影響效果。水量調(diào)控模型需要對現(xiàn)狀和未來的水資源供需分析進(jìn)行模擬計算,通過情景方案設(shè)置和反饋協(xié)調(diào)河道外生活、生產(chǎn)和生態(tài)用水與河道內(nèi)生態(tài)、發(fā)電以及航運等用水,提出流域水資源在不同地區(qū)和行業(yè)的合理配置方案,并給出不同用水和工程調(diào)度情景下各區(qū)域的廢污水排放狀況和重要節(jié)點斷面的水量過流狀況,并給出各個水功能分區(qū)節(jié)點的過流狀況。根據(jù)上述要求,本實施例所述的水量調(diào)控模擬模型的主要設(shè)計原則如下(I)確定模擬規(guī)模與精度根據(jù)水量總量調(diào)控需求,設(shè)置基本的供用水計算單元和斷面節(jié)點。為配合水功能區(qū)的納污能力計算,采用水資源分區(qū)套行政區(qū)作為水量供需平衡計算單元,主要控制節(jié)點為大型以上蓄水工程、重要引提水工程、重要省界斷面以及水功能區(qū)節(jié)點斷面。(2)流域水循環(huán)與水量供用耗排二元水循環(huán)關(guān)系為基礎(chǔ)根據(jù)流域天然水系與計算單元的供用耗排水關(guān)系繪制描述水量關(guān)系的水資源網(wǎng)絡(luò)節(jié)點圖。將計算單元、主要工程、控制性節(jié)點以及供用耗排水等系統(tǒng)元素采用概化的“點”、“線”元素表達(dá),反映流域水循環(huán)與水資源供、用、耗、排過程,作為模擬計算的基礎(chǔ)。(3)約束條件在模型中設(shè)置反映各分區(qū)水資源配置中的工程和非工程約束,如供水渠道、管道的過水能力約束,不同來水情況下分水指標(biāo)的約束等。對于發(fā)電工程考慮其發(fā)電調(diào)度需求以及與供水、河道內(nèi)水量調(diào)度需求的關(guān)系,設(shè)置河道內(nèi)水量需求約束條件。(4)模塊化處理模塊化處理即是根據(jù)計算過程中各部分功能將系統(tǒng)整體劃分為相互關(guān)聯(lián)但又具備一定獨立性的功能模塊,針對各模塊的過程進(jìn)行詳細(xì)設(shè)計的方法。模型模擬中根據(jù)各類水源的利用和處理定義了多個相對封閉和獨立的子系統(tǒng),同時對每一類子系統(tǒng)又規(guī)范了其計算范圍和框架。通過模塊化分割可以對各個子系統(tǒng)功能作封閉式處理,并通過各模塊有機(jī)整合實現(xiàn)系統(tǒng)功能。通過模塊化分割處理可以深入細(xì)致地處理各功能塊內(nèi)部的細(xì)節(jié)流程,使系統(tǒng)框架更清晰、可擴(kuò)展性更強。計算水功能區(qū)控制斷面流量的步驟用于使用水量調(diào)控模擬模型計算水功能區(qū)的斷面流量。本步驟使用水功能區(qū)斷面的流量計算公式,采用不同的用水情景時,調(diào)整水功能區(qū)所在單元取用水量U,計算出水功能區(qū)控制斷面在不同用水情景下的過流量,分析出水量變化條件下的水質(zhì)狀況模擬和控制目標(biāo),實現(xiàn)了人工用水與水循環(huán)過程動態(tài)關(guān)系的模擬。水量調(diào)控模擬模型通過設(shè)置供用水計算單元和節(jié)點進(jìn)行水量平衡計算。為了配合水功能區(qū)的納污能力計算,采用水資源分區(qū)嵌套行政區(qū)作為計算單元。計算單元是系統(tǒng)模擬的中心環(huán)節(jié),水資源供用耗排的分析都是在計算單元內(nèi)部完成的。節(jié)點包含工程節(jié)點和控制節(jié)點兩類工程節(jié)點特指在系統(tǒng)圖上單列的蓄引提工程,這類工程在模擬計算中單獨參與計算,可以對受水單元或者其他工程供水,同時也可以按照工程運行目標(biāo)對發(fā)電航運和生態(tài)與環(huán)境等河道內(nèi)用水需求進(jìn)行水量調(diào)控??刂乒?jié)點指有水量或水質(zhì)控制要求的河道或渠道斷面,它具有與工程節(jié)點相同的各種水力關(guān)系,但一般不具備調(diào)蓄能力。將計算單元、主要工程節(jié)點、控制節(jié)點以及供用耗排水等系統(tǒng)元素采用概化的“點”、“線”元素表達(dá),繪制描述流域水量關(guān)系的水資源系統(tǒng)網(wǎng)絡(luò)節(jié)點圖,反映流域水循環(huán)與水資源供用耗排過程,以此作為模擬計算的基礎(chǔ)。水量調(diào)控模型以行政區(qū)域和流域片區(qū)結(jié)合的計算單元為基本用水單位模擬水量 過程,可以較為準(zhǔn)確的反映天然水循環(huán)過程,從整體上對流域水循環(huán)過程進(jìn)行仿真模擬。相對與現(xiàn)有的水質(zhì)模擬主要基于河流層面信息進(jìn)行污染負(fù)荷入河過程的分析,降低了模擬的計算步長和計算單元的尺度,擴(kuò)大了可模擬的范圍,可以更好的反映較大范圍內(nèi)區(qū)域的水循環(huán)、供用水過程條件下的區(qū)域整體水量過程和相應(yīng)的納污能力指標(biāo)。水量調(diào)控模型主要以模擬方式完成水量供需平衡計算、用戶用耗排水和河道水量平衡過程。以一個計算單元為例,供需平衡計算公式為
權(quán)利要求
1.一種基于水功能區(qū)的水量水質(zhì)調(diào)控方法,其特征在于所述方法的步驟如下 構(gòu)建水量調(diào)控模擬模型的步驟用于以流域水資源分區(qū)和行政分區(qū)嵌套形成水功能區(qū)的計算單元,并將水功能區(qū)的主要控制工程節(jié)點和水功能區(qū)控制斷面作為基本信息模擬單元;分析所述水功能區(qū)范圍內(nèi)的水量平衡信息,建立計算單元與水功能區(qū)地表徑流匯流與用水退水關(guān)系,計算用水影響下的水功能區(qū)控制斷面過流;所述的水功能區(qū)斷面的流量計算公式如下Q = W 上+W 本-U+S 式中Q_水功能區(qū)控制斷面過流量, W上-上游來水量, W*-水功能區(qū)所在單元地表徑流量, U-水功能區(qū)所在單元地表取用水量, S-水功能區(qū)所在單元用水后退入水功能區(qū)斷面水量; 計算水功能區(qū)控制斷面流量的步驟用于使用水量調(diào)控模擬模型計算水功能區(qū)的斷面流量; 計算入河污染負(fù)荷量的步驟用于對入河的污染源進(jìn)行量化的估算,包括對工業(yè)點源產(chǎn)生量估算和城鎮(zhèn)生活點源產(chǎn)生量估算; 基于水功能區(qū)的納污能力計算的步驟用于根據(jù)水量調(diào)控模擬模型的流量計算結(jié)果,采用一維水質(zhì)模擬公式計算用水影響條件下的水功能區(qū)納污能力 M = 31. 5 (Cs-C0 exp (-kL/u)) exp (kx/u) Qr, 式中M_污染物納污能力(t/a), Cs-水功能區(qū)下斷面水質(zhì)目標(biāo)控制因子濃度(mg/L), C0-水功能區(qū)上斷面水質(zhì)目標(biāo)控制因子濃度(mg/L), k-污染物綜合衰減系數(shù)(1/s), L-水功能區(qū)長度(m), X-排污口下游斷面距水功能區(qū)下斷面縱向距離(m), U-設(shè)計流量下的平均流速(m/s), Qr-設(shè)計流量(m3/s); 分析水功能區(qū)達(dá)標(biāo)狀況的步驟用于使用水功能區(qū)達(dá)標(biāo)分析模型,通過污染負(fù)荷產(chǎn)生入河量與水功能區(qū)納污能力對比分析水功能區(qū)達(dá)標(biāo)狀況,如果若水功能區(qū)達(dá)標(biāo)率不能實現(xiàn),則回到“構(gòu)建水量調(diào)控模擬模型的步驟”調(diào)整數(shù)據(jù)量,重新進(jìn)行計算,直至實現(xiàn)水功能區(qū)達(dá)標(biāo)率,則結(jié)束過程。
2.根據(jù)權(quán)利要求I所述的方法,其特征在于,所述的計算水功能區(qū)控制斷面流量的步驟中用于根據(jù)模塊化分割方法水量調(diào)控模擬模型的計算包括的子步驟如下 計算單元本地徑流用水模擬包括第一階段計算本地徑流的利用及用戶間配置,第二階段計算河網(wǎng)入流統(tǒng)計及供水計算,第三階段計算網(wǎng)超蓄水量排放; 計算單元非常規(guī)水源利用模擬非常規(guī)水源主要包括雨水利用、微咸水和海水利用; 地下水利用模擬對于淺層地下水,將每個計算單元視為只有一個地下水庫,其補給水量為地下水可利用量,對每個單元輸入地下水可開采量的年過程; 地表工程供水模擬對于大型水庫、引水工程等,以本地徑流計算后得出的水量作為入庫水量,根據(jù)水庫入流、上時段存蓄水量結(jié)合水庫調(diào)度規(guī)則得出可供水量,對用戶進(jìn)行水量配置,并計算超過蓄水能力后的工程下泄流量; 控制斷面水平衡模擬根據(jù)計算單元對水庫節(jié)點、水功能區(qū)控制斷面等的徑流匯流量、單元用水耗水后的排水量,通過水量平衡計算各斷面的過流量; 時段結(jié)束判斷判斷計算的時段是否到達(dá)設(shè)定的末尾,如果到達(dá)則結(jié)束過程,如果未達(dá)到則回到“計算單元本地徑流用水模擬”的子步驟,繼續(xù)下移時段的計算。
3.根據(jù)權(quán)利要求I或2所述的方法,其特征在于,所述的水量水質(zhì)聯(lián)合分析調(diào)控的步驟中的分析方法包括兩種分析方式 第一種分析方式是在確定用水總量條件下模擬各水功能區(qū)斷面的徑流過程,得出各區(qū)域納污能力總量;第二種分析方式是通過對比納污能力計算結(jié)果與實際排污預(yù)測量,調(diào)整用水控制量實現(xiàn)水功能區(qū)斷面水質(zhì)達(dá)標(biāo)目標(biāo)。
全文摘要
本發(fā)明涉及一種基于水功能區(qū)的水量水質(zhì)調(diào)控方法,所述方法的步驟包括構(gòu)建水量調(diào)控模擬模型的步驟;計算水功能區(qū)控制斷面流量的步驟;計算入河污染負(fù)荷量的步驟;基于水功能區(qū)的納污能力計算的步驟;分析水功能區(qū)達(dá)標(biāo)狀況的步驟。本發(fā)明通過構(gòu)建流域級水功能區(qū)的水量調(diào)控模擬模型,考慮人工用水、工程調(diào)度與自然水循環(huán)關(guān)系,同時構(gòu)建水功能區(qū)納污能力分析模型,使用納污能力分析模型以水量調(diào)控模擬模型為依據(jù)分析計算水功能區(qū)的納污能力,根據(jù)分析計算結(jié)果對流域的水量和水質(zhì)進(jìn)行調(diào)控。與現(xiàn)有技術(shù)相比,本發(fā)明充分考慮人工用水、工程調(diào)度與自然水循環(huán)三方面的因素,對流域水質(zhì)的納污能力分析更加準(zhǔn)確,使水量調(diào)控更加合理。
文檔編號G06F19/00GK102750448SQ20121019101
公開日2012年10月24日 申請日期2012年6月11日 優(yōu)先權(quán)日2012年6月11日
發(fā)明者游進(jìn)軍, 牛存穩(wěn), 薛小妮, 褚俊英, 賀華翔, 賈仰文, 賈玲, 魏娜 申請人:中國水利水電科學(xué)研究院