• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    基于水文水動力學(xué)耦合的洪水預(yù)報模型研究及應(yīng)用

    2017-02-27 23:36:42楊甜甜梁國華何斌李瑩芹
    南水北調(diào)與水利科技 2017年1期

    楊甜甜+梁國華+何斌+李瑩芹

    摘要:為了提高下墊面變化劇烈流域的洪水預(yù)報精度,在傳統(tǒng)流域水文模型的基礎(chǔ)上耦合水動力學(xué)模型,建立水文水動力耦合洪水預(yù)報模型。首先利用水文模型獲得某一斷面的流量過程作為水動力學(xué)模型的邊界條件;之后利用一維水動力學(xué)模型進行河道洪水演進計算,推求流域出口斷面的流量過程;最后以煙臺市外夾河流域為例進行驗證。結(jié)果表明,所建水文水動力耦合模型模擬的產(chǎn)流合格率較高,流量過程與實測值吻合,在一定程度上彌補了集總式水文模型不能考慮河道內(nèi)復(fù)雜水流運動的不足,因此對具有復(fù)雜水文、水力條件的流域的洪水預(yù)報具有重要的指導(dǎo)意義。

    關(guān)鍵詞:大伙房流域模型;一維水動力學(xué)模型;水文水動力學(xué)耦合模型;洪水預(yù)報

    中圖分類號:TV122;P338文獻標(biāo)志碼:A文章編號:1672-1683(2017)01-0072-07

    Abstract:In order to improve the accuracy of flood forecasting in basins with dramatic changes in the base course,this paper proposes a coupled hydrological-hydrodynamic flood forecasting model which couples the traditional hydrological model with the one-dimensional hydrodynamic model.First,the flow process of a section was obtained by the hydrological model and was used as the boundary condition of the one-dimensional hydrodynamic model.Then the one-dimensional hydrodynamic model was applied to simulate the flood routing in the river channel.The proposed method was tested on the Waijia River basin in Yantai City,China.The results indicate that the coupled model can simulate runoff yield with a high qualification rate,and generate flow processes that tally with the measured values.It can be concluded that this coupled hydrological-hydrodynamic model may make a significant contribution to flood forecasting in basins with complex hydrologic and hydraulic conditions thanks to its ability to model water movement in the basin.

    Key words:DHF model;one-dimensional hydrodynamic model;coupled hydrological-hydrodynamic model;flood forecasting

    流域洪水預(yù)報作為洪水調(diào)度的基礎(chǔ),是防洪減災(zāi)的一項重要的非工程措施。20世紀60年代以來,世界各國都投入大量人力物力研制了適用于本國流域特點的洪水預(yù)報模型[1]。通常,在模型中將流域徑流過程分為產(chǎn)流和匯流兩部分。相對產(chǎn)流而言,匯流過程不僅受產(chǎn)流過程的影響而且還受流域地形、河網(wǎng)、河道等多種因素的影響,其規(guī)律多變不易描述。傳統(tǒng)的集總式流域水文模型通常采用馬斯京根法或其它水文學(xué)方法進行演進計算。但是,由于僅能根據(jù)水量平衡方程和槽蓄方程進行流量計算,模擬結(jié)果往往不甚理想。尤其是在丘陵、平原混合區(qū)域,河道內(nèi)水流運動受地形條件、上游來水、回水頂托等因素的影響,傳統(tǒng)的洪水預(yù)報模型中河網(wǎng)匯流部分無法描述這些影響因素而導(dǎo)致流域出口斷面的洪水預(yù)報精度降低[2-7]。

    國內(nèi)外學(xué)者針對該問題進行了一些嘗試和研究[8-9],在水文學(xué)方法的實際應(yīng)用中加入河道水流運動物理機制的描述,以求盡可能地反映流域匯流真實情況,通常借助計算機等工具求解圣維南方程組,得到河道內(nèi)不同斷面的流量和水位,進而補充傳統(tǒng)洪水預(yù)報無法得到的信息,提高預(yù)報精度。然而,我國幅員遼闊、流域各有特點,洪水預(yù)報模型在實際應(yīng)用上仍然存在著許多問題。因此,在充分考慮流域產(chǎn)流、匯流特點,研制適合流域?qū)嶋H情況的水文水動力學(xué)耦合的洪水預(yù)報模型,顯得尤為重要。本文以煙臺市外夾河流域為研究背景,在傳統(tǒng)洪水預(yù)報的匯流部分加入計算水力學(xué)的求解方法,采用一維水動力學(xué)方法進行模擬,結(jié)合水文學(xué)模型構(gòu)建水文水動力學(xué)耦合模型以彌補水文學(xué)模型不能考慮河道內(nèi)復(fù)雜水流運動的不足。

    1 耦合洪水預(yù)報模型的構(gòu)建

    1.1 水文模型

    降雨徑流是流域洪水的重要成因,水文模型是模擬流域上游區(qū)域以及中下游區(qū)間產(chǎn)、匯流過程的重要工具,在我國應(yīng)用較廣泛的主要是基于蓄滿產(chǎn)流的新安江流域水文模型和基于超滲產(chǎn)流的大伙房流域模型[10](DHF)。其中,DHF模型于1973年由遼寧省大伙房水庫管理局提出,由產(chǎn)流和匯流兩部分模型組成,是一種集總式概念性水文模型。其中,產(chǎn)流模型采用超滲產(chǎn)流機理進行計算,以拋物線描述流域表層蓄水容量的空間分配曲線,以霍爾登下滲曲線為基礎(chǔ)亦采用拋物線描述雙層下滲率分配曲線;匯流模型采用據(jù)多年預(yù)報經(jīng)驗總結(jié)得到的變強度、變匯流速度的經(jīng)驗單位線進行計算。DHF模型共16個產(chǎn)匯流參數(shù),需優(yōu)選確定的只有6個參數(shù)。整個模型結(jié)構(gòu)合理,數(shù)學(xué)方程簡單,便于計算運用,模型的模擬精度和檢驗精度較高,已在東北的桓仁、碧流河、石佛寺等10多個流域的洪水預(yù)報中得到良好的應(yīng)用[11-13]。

    根據(jù)研究流域的特點,選擇適用于北方干旱半干旱地區(qū)的DHF作為本文耦合模型中的水文模型部分來計算研究流域上游斷面的出流,將其結(jié)果作為水動力學(xué)模型的輸入條件。

    1.2 水動力學(xué)模型

    一維水動力學(xué)模型根據(jù)離散方式的不同有多種求解方法,包括特征線法、有線差分法、有限體積法等。其中有限差分法的Preissmann四點偏心隱格式在江河的洪水演進計算中應(yīng)用最為廣泛,且其計算精度、計算效率、穩(wěn)定性完全滿足本文模型需求。因此,采用Preissmann四點偏心隱格式對描述天然河道一維水流運動的圣維南方程組離散,使用追趕法對離散后方程組進行求解。根據(jù)不同的邊界條件,可設(shè)不同遞推關(guān)系,依次求得河道內(nèi)各斷面的水位、流量值[14-16]。

    水動力學(xué)模型的對象往往是一個大型、復(fù)雜的非恒定水流系統(tǒng),為了盡可能真實地模擬實際河道中的水流運動狀態(tài),在應(yīng)用時需要對天然河道的水力學(xué)參數(shù)及模型的邊界條件、初始條件、參數(shù)條件進行明確。

    (1)天然河道地形復(fù)雜,斷面形狀不規(guī)則,給水力學(xué)參數(shù)的計算帶來困難,在實際應(yīng)用中通常將其進行概化,使其在基本符合斷面形狀變化的基礎(chǔ)上又能降低模型復(fù)雜度提高計算效率。

    (2)上游邊界條件是非恒定流發(fā)生過程中,首斷面應(yīng)滿足的水力條件。它常以流量過程線或水位過程線來表示,而下邊界條件通常以水位過程線或水位流量關(guān)系曲線來表示。

    (3)除特殊條件的處理外,模型要想穩(wěn)定計算,必須給予能反映河道初始狀態(tài)的水位、流量值,可以是恒定狀態(tài)的,也可以是非恒定狀態(tài)的。

    (4)模型中的待定參數(shù)主要有反映固體邊界作用的糙率系數(shù),差分格式的權(quán)系數(shù)θ,時間步長和空間步長。河道的糙率因影響因素眾多,很難準(zhǔn)確求得,模型可在河道實際情況基礎(chǔ)上給出糙率初值,然后根據(jù)計算結(jié)果進行優(yōu)化;權(quán)系數(shù)直接影響模型穩(wěn)定性和計算精度,一般選取大于等于0.6的強穩(wěn)定狀態(tài);時間步長根據(jù)計算需求及穩(wěn)定性要求可取30 s、60 s等;空間步長取兩斷面間距,為不等距步長,大約為100~500 m,其中在河道斷面形狀變化較劇烈和彎道河段,距離步長較小,順直河道距離步長較大,對斷面間距較大者根據(jù)實際需求進行插值。

    1.3 耦合模型構(gòu)建方法

    集總式水文模型專注于降雨、徑流的產(chǎn)生過程和規(guī)律,將整個水文系統(tǒng)看作一個整體而進行研究,不考慮系統(tǒng)內(nèi)部的不均勻性和河道水流運動規(guī)律。在流域下墊面均一、形狀規(guī)則、區(qū)域內(nèi)河網(wǎng)較簡單、匯流規(guī)律較單一的情況下,單獨使用集總式水文模型進行洪水預(yù)報不僅可以降低模型復(fù)雜程度而且也能取得較為理想的結(jié)果。但是,大多數(shù)流域的下墊面不會一成不變,所以忽略水流運動規(guī)律的集總式水文模型無法滿足工程需要,其預(yù)報精度會受地形、河網(wǎng)等影響。水動力學(xué)模型著重于河道水流運動規(guī)律的詳實描述,可以反映河道內(nèi)細微差別。因此,兩者耦合使用正好可以彌補傳統(tǒng)集總式水文模型的不足。

    耦合模型的建立重點在其耦合方法上,本文選擇連接的方式進行耦合,見圖1。采用大伙房流域水文模型進行河道上邊界斷面和各支流控制流域出口斷面的流量過程的預(yù)測;將其計算結(jié)果作為上邊界條件和內(nèi)邊界條件直接輸入到一維水動力學(xué)模型中;一維水力學(xué)模型的數(shù)據(jù)不反饋輸入到水文學(xué)模型,二者是單向輸入的。

    2 耦合模型的應(yīng)用

    2.1 研究流域

    煙臺市外夾河流域(圖2)面積1 072 km2,平均長度60 km,寬度15.5 km,形狀系數(shù)0.26,上游寬、下游窄,呈“梨”型;河道干流全長64 km,干流平均比降1.32‰。流域為典型的丘陵、平原混合區(qū)域,其上游譚家莊水文站以上區(qū)域和各支流上游區(qū)域主要為山地、丘陵區(qū);而譚家莊站與福山站之間以及主干河道兩側(cè)則多以平原、洼地為主,因此不同下墊面條件導(dǎo)致區(qū)域內(nèi)匯流時間有很大差異。同時,福山水文站距離入海口僅有8 km,河道水流受海潮頂托回水現(xiàn)象嚴重,導(dǎo)致出口斷面——福山水文站的洪水預(yù)報精度降低;另外,外夾河河道彎曲度大、河床多有突起,極易壅水漫灘威脅兩岸安全。

    本文選取外夾河福山站以上流域建立大伙房流域模型和水文水動力耦合模型,用于驗證水文水動力學(xué)耦合模型的合理性和對洪水預(yù)報精度的改善。研究區(qū)域內(nèi)有福山站、譚家莊站、古現(xiàn)村站、桃村站、清泉埠站、鐵口站6個雨量站,均有長系列實測降雨資料;應(yīng)用大伙房流域模型和水文移置法進行研究區(qū)域的降雨徑流模擬計算,得到譚家莊水文站的流量過程和區(qū)域中譚家莊水文站與福山水文站間間的各支流流量過程,作為水動力學(xué)模型的上邊界條件和內(nèi)邊界條件。水動力學(xué)模型計算河段選擇外夾河譚家莊水文站至福山水文站河段,河段劃分121個斷面,5個特殊結(jié)點,全長約27.9 km。譚家莊站、福山站是外夾河上比較大的水文站,有實測水位、流量資料,本文利用福山站來檢驗耦合模型的可靠性。

    2.2 模型條件處理

    水文水動力學(xué)耦合模型中其產(chǎn)匯流特性及水流特征不光遵循特定理論,還有一定的經(jīng)驗性,單純地注重數(shù)學(xué)模型的處理不足以保證模型模擬得到良好的效果,還必須對研究區(qū)域及河段的特征信息進行補充并加以經(jīng)驗性的處理,即需對實際工程情況進行特殊處理。

    (1)區(qū)間入流的處理。

    研究區(qū)域中譚家莊水文站與福山水文站之間河段還有沐浴河、友誼河、東風(fēng)河、芹河、橫河5條支流入流,但各支流流域控制面積較小,河流短,汛期洪量小,并且缺少各支流實測斷面資料及水文資料。為降低耦合模型中一維水動力學(xué)模型的復(fù)雜程度,將其簡化處理不做河網(wǎng)計算,而是作為主干河道旁側(cè)入流。因此,需要對各支流流域進行降雨徑流模擬,出流直接耦合于水動力學(xué)模型。但是,由于外夾河各支流缺乏實測水文資料,并且各支流流量較小,其入流對主干河道水量、峰現(xiàn)時間等影響不大,因此不采用水文模型對各支流控制區(qū)域進行流量模擬,而是采用水文移置的方法進行處理:直接根據(jù)各支流控制面積和譚家莊以上區(qū)域的面積比,移用譚家莊水文站以上區(qū)域大伙房流域模型計算結(jié)果。其中,東風(fēng)河與芹河入流點距離較近,為了降低模型復(fù)雜程度將其簡化為一條支流來考慮。各支流控制流域面積與譚家莊以上區(qū)域面積比見表1。

    (2)河道斷面處理。

    本研究河段上游屬于天然狀態(tài)河槽,斷面窄而深,基本無灘地、河心淤積帶的存在,而下游處于城鎮(zhèn)地帶,人工修筑堤防并進行過大規(guī)模的河道治理,斷面形狀變化不劇烈。在本模型中,將上游斷面的尖銳突起部分修整為規(guī)則的梯形或復(fù)式斷面;而下游斷面較為規(guī)則,將其概化為矩形、梯形或復(fù)式斷面。由于研究河段的流量集中在主槽內(nèi),灘地的行洪流量較小,所以在本文計算中暫不考慮復(fù)式斷面的水流差異。

    (3)模型邊界條件處理。

    福山水文站是研究河段外夾河與鄰近清洋河匯合之前的一個控制站,同時也是水動力學(xué)模型的下邊界,距入??? km。由于清洋河上游為門樓水庫,汛期受水庫泄洪及潮水頂托影響,存在死水及回水情況,增加了福山斷面邊界條件處理的復(fù)雜性。研究比較下列方案:A方案給定福山站水位流量關(guān)系;B方案給定福山水位過程;C方案給定福山流量過程。試算結(jié)果表明,福山邊界條件的給定與模型計算的準(zhǔn)確性及穩(wěn)定性密切相關(guān)。若要強制給定福山水位過程或流量過程將致使水位與流量不協(xié)調(diào),從而使計算失穩(wěn)??紤]到本模型側(cè)重于干流福山站流量、水位的精確預(yù)報,通過多種方案的比較最后選定的方案A。

    (4)模型初始條件處理。

    模型中的水動力學(xué)部分采用恒定初始狀態(tài),即以任一初始水位,推算出各計算斷面在足夠長時間恒定流量(采用場次洪水初始時刻實測流量值)輸入的情況下到達穩(wěn)定狀態(tài)下的河道內(nèi)各斷面處水位和流量,并以此值作為模型計算的初始值,進行非恒定流的計算。

    (5)模型參數(shù)條件處理。

    水動力學(xué)模型的空間步長采用相鄰斷面間距長度,為保證模型計算的穩(wěn)定性及空間步長的合理性,對斷面間距較大的按500 m的距離利用相鄰斷面進行插值。時間步長根據(jù)計算需求及穩(wěn)定性要求取60 s,權(quán)重θ取0.8,既保證模型精度又計算穩(wěn)定。根據(jù)研究河段現(xiàn)有資料情況以及精度、計算速度等的要求,研究河段上游基本處于天然河道狀態(tài),糙率取0.05;下游處于城鎮(zhèn)地帶,河床、岸壁的粗糙度小,糙率取0.028。

    2.3 模型結(jié)果與分析

    (1)耦合模型模擬結(jié)果。

    選用流域出口斷面——福山水文站1976年-1990年期間資料完整性較好、有代表性的16場洪水進行模擬計算,從而驗證所采用的水文水動力耦合模型能否有效地描述流域產(chǎn)匯流特性、能否真實反映河道的水流情況,模型參數(shù)及地形、斷面的概化是否合理,以及對洪水預(yù)報精度是否有所改善。模擬結(jié)果如表2所示。

    從耦合模型的產(chǎn)流模擬結(jié)果可知:a.其產(chǎn)流模擬僅有750813場次不合格,其余場次洪水徑流深預(yù)報相對誤差均在20%以內(nèi),合格率達94%;b.760629、790731、850726三個場次洪水總徑流量預(yù)報偏大,原因為流域內(nèi)塘壩、洼地等的截流作用未納入水文模型中,導(dǎo)致在極度干旱情況下的產(chǎn)流量偏大;c.大部分場次洪水產(chǎn)流模擬值偏小,因為譚家莊上游區(qū)域內(nèi)植被情況較好,下游區(qū)域為城鎮(zhèn),不透水面積大,而模型構(gòu)建時對上、下游區(qū)域進行均一處理未考慮其下墊面的變化。

    從耦合模型的匯流模擬結(jié)果可知:a.有15場洪水的洪峰流量和峰現(xiàn)時間的預(yù)報誤差在允許范圍內(nèi),合格率為94%;b.確定性系數(shù)有10場在0.7以上,只有730821場次洪水在0.5以下,流量過程擬合較好;c.從圖3所示的4場典型洪水過程模擬結(jié)果看,耦合模型對于洪水流量過程整體上擬合良好,所確定的參數(shù)基本能反映流域產(chǎn)匯流特性以及河道水流特性;但是,對于峰前、峰后小洪水以及復(fù)峰型洪水反映不靈敏,如760818場次洪水的大洪峰前的小峰起峰不明顯,850818場次洪水的退水?dāng)M合不好,分析原因是主干各支流入流進行了概化處理,未考慮其產(chǎn)匯流差異性。

    綜上所述,耦合模型產(chǎn)、匯流模擬合格率均為94%,且福山水文站模擬流量過程與實測值擬合良好,說明所建立的水文水動力耦合模型能夠較好地反映流域的降雨-徑流響應(yīng)過程。

    (2)與大伙房模型結(jié)果對比分析。

    單獨使用大伙房模型對福山水文站以上流域進行場次洪水模擬,其結(jié)果如表3所示。與水文水動力耦合模型模擬結(jié)果進行對比分析可知:a.大伙房模型模擬的產(chǎn)流結(jié)果合格率為81%,低于耦合模型的94%;匯流模擬結(jié)果合格率只有56%,遠低于耦合模型的94%;b.大伙房模型模擬的洪水過程確定性系數(shù)只有7場在0.7以上,最高只有0.85,低于耦合模型的0.97;另有820825場次和900724場次洪水確定性系數(shù)過小,說明大伙房模型對洪水過程的模擬偏差較大,不能很好地反映流域匯流特性;c.大伙房模型模擬的峰現(xiàn)時間大部分都有1個時段的差值,其中730831場次洪水峰現(xiàn)時間差為4個時段,預(yù)報不合格。d.從圖3所示的4場典型洪水過程模擬結(jié)果可以看出,大伙房模型模擬的洪水流量過程與實際洪水?dāng)M合效果不是很好,退水較快。

    綜上所述,單純地使用大伙房流域水文模型進行流域洪水模擬往往會影響匯流的準(zhǔn)確性,而水動力學(xué)模型的加入使得匯流更加貼合流域?qū)嶋H情況。

    (3)1990年以后洪水模擬。

    煙臺市外夾夾河自1990年起修建了大量的橡膠壩、攔河閘以及跨河大橋等工程,致使流域的產(chǎn)、匯流條件有較大的改變。但是,由于缺乏閘壩參數(shù)及運行資料,耦合模型建立過程中并未考慮河道中的人工建筑物。本文利用建立的耦合模型對1990年以后的9場洪水進行模擬計算,分析其是否對流域產(chǎn)匯流有影響,從而驗證本文所建耦合模型的適用性。

    分析表4的模擬結(jié)果可知:a.1990年以后洪水的產(chǎn)、匯流模擬結(jié)果的合格率分別為66.7%、88.9%,較1990年前洪水的模擬合格率有所降低;b.大部分場次洪水產(chǎn)流模擬值較實際值偏大,而1990年以前洪水模擬值大部分偏小。雖然在實際工程運用中,當(dāng)汛期遇大洪水時各橡膠壩均將其壩頂高程降到最低,但其對洪水依然有攔蓄作用,對產(chǎn)流產(chǎn)生一定影響;c.不合格場次洪水的產(chǎn)流量、洪峰流量值均較小,說明橡膠壩、攔河閘、跨河大橋等工程的修建會會大大降低小洪水模擬的精確性,而對大洪水影響作用不明顯。綜上所述,雖然河道人工建筑物的修建改變了流域的產(chǎn)匯流條件使得耦合模型無法很好地反映小洪水的產(chǎn)匯流特性,但其對大洪水依然適用。

    3 結(jié)論

    本文建立了復(fù)雜下墊面條件的水文預(yù)報和天然河道復(fù)雜地形和水流條件的水動力學(xué)計算的耦合洪水預(yù)報模型,并利用煙臺市外夾河流域的實測水文資料進行模型的應(yīng)用驗證,結(jié)果如下。

    (1)本文所建的水文水動力耦合模型模擬的洪水流量過程與實測過程擬合良好,模型計算穩(wěn)定、快速,可用于實際洪水的預(yù)報。

    (2)通過與大伙房模型的模擬結(jié)果對比發(fā)現(xiàn),耦合模型能夠更好地解決單純使用集總式水文模型而不能反映河道水流運動的弊端,大大提高了流域洪水預(yù)報模擬精度,更具有通用性。

    (3)限于數(shù)據(jù)資料的缺乏,模型實際應(yīng)用過程中采取了一些概化的方法,同時未考慮流域內(nèi)小水庫、小型塘壩和洼地等的調(diào)蓄作用和橡膠壩、攔河閘、跨河大橋等建筑物對河道水流運動的影響。因此,對于小洪水受以上因素影響較大,導(dǎo)致本文所建模型的模擬與實測值有所偏差,而其更適用于大洪水的模擬。

    參考文獻(References):

    [1] 徐宗學(xué).水文模型:回顧與展望[J].北京師范大學(xué)學(xué)報:自然科學(xué)版,2010,46(3):278-289.(XU Zong-xue.Hydrological models:past,present and future[J].Journal of Beijing Normal University:Natural Science,2010,46(3):278-289.(in Chinese))

    [2] 常露,劉開磊,姚成,等.復(fù)雜河道洪水預(yù)報系統(tǒng)研究-以淮河王家壩至小柳巷區(qū)間流域為例[J].湖泊科學(xué),2013,25(3):422-427.(CHANG Lu,LIU Kai-lei,YAO Cheng,et al.Real-time flood forecasting system for complicated river channels:a case study from Wangjiaba to Xiaoliuxiang section in the Huaihe River basin[J].Lake Science,2013,25(3):422-427.(in Chinese))

    [3] 徐興亞,方紅衛(wèi),張岳峰,等.基于粒子濾波同化的河道洪水實時概率預(yù)報模型與應(yīng)用[J].水科學(xué)進展,2015,3.(XU Xing-ya,F(xiàn)ANG Hong-wei,ZHANG Yue-feng,et al.A real-time probabilistic channel flood forecasting model and application based on particle filters[J].Advances in Water Science,2015,3.(in Chinese))

    [4] 薛天柱.洪水演算方法在白龍江河段洪水預(yù)報中的應(yīng)用研究[D].蘭州:蘭州大學(xué),2012.(XUE Tian-zhu.Method in Bailong River Flood Forecasting[D].Lanzhou:Lanzhou University,2012.(in Chinese))

    [5] 趙人俊.流域水文模型[M].北京:水利電力出版社,1983.(ZHAO Ren-jun.Hydrologic models of watersheds[M].Beijing:Water Conservancy and Electric Power Press,1983.(in Chinese))

    [6] 程亮,王宗志,胡四一,等.強烈下滲條件下天然河道洪水演進模擬方法[J].中國科學(xué):地球科學(xué),2015,45(2):207-215.(CHENG Liang,WANG Zong-zhi,HU Si-yi,et al.Flood routing model incorporating intensive streambed infiltration[J].Scientia Sinica Terrae,2015,45(2):207-215.(in Chinese))

    [7] 包紅軍,李致家,王莉莉.具有行蓄洪區(qū)的復(fù)雜水系實時洪水預(yù)報研究[J].水力發(fā)電學(xué)報,2009(4):5-12.(BAO Hong-jun,LI Zhi-jia,WANG Li-li.Real-time flood forecasting for complex watershed with flood diversion and flood retarding areas[J].Journal of Hydroelectric Engineering,2009(4):5-12.(in Chinese))

    [8] 王船海,李光熾.實用河網(wǎng)水流計算[M].南京:河海大學(xué)出版社,1997.(WANG Chuan-hai,LI Guang-zhi.Practical guide on river network flow calculation[M].Nanjing:Hohai University Press,1997.(in Chinese))

    [9] 楊國錄.河流數(shù)學(xué)模型[M].北京:海洋出版社,1993.(YANG Guo-lu.Mathematical model of river[M].Beijing:China Ocean Press,1993.(in Chinese))

    [10] 大連理工大學(xué),國家防汛抗旱總指揮部辦公室.水庫防洪預(yù)報調(diào)度方法及應(yīng)用[M].北京:中國水利水電出版社,1996:43-52.(Dalian University of Technology,The Office of Flood Control and Drought Relief Headquarters.Method and application of reservoir flood forecasting operation[M].Beijing:China Water & Power Press,1996.(in Chinese))

    [11] 王猛,彭勇,梁國華.大伙房模型在石佛寺流域洪水預(yù)報方案中的應(yīng)用[J].南水北調(diào)與水利科技,2012,10(2):93-97.(WANG Meng,PENG Yong,LIANG Guo-hua.Application of DHF model in flood forecasting for Shifosi Basin[J].South-to-North Water Transfers and Water Science & Technology,2012,10(2):93-97.(in Chinese))

    [12] 馬玉軍.碧流河水庫流域洪水預(yù)報方案及設(shè)計洪水過程修勻研究[D].大連:大連理工大學(xué),2013.(MA Yu-jun.Study on flood forecasting scheme and graduation of design flood hydrograph of Biliuhe Reservoir[D].Dalian:Dalian University of Technology,2013.(in Chinese))

    [13] 李娜.基于知識的桓仁水庫實時洪水預(yù)報方法應(yīng)用研究[D].大連:大連理工大大學(xué),2008.(LI Na.Study and application on real time flood forecast method of Huanren Reservoir based on knowledge[D].Dalian:Dalian University of Technology,2008.(in Chinese))

    [14] 杜獻梅,趙克玉,馮玉權(quán).沁河中游復(fù)式斷面洪水演進模擬[J].西北農(nóng)林科技大學(xué)學(xué)報:自然科學(xué)版,2013(2):223-227.(DU Xian-mei,ZHAO Ke-yu,F(xiàn)ENG Yu-quan.Simulation of flood routing with compound cross-section in middle reach of the Qin River[J].Journal of Northwest A & F University:Natural Science Edition,2013(2):223-227.(in Chinese))

    [15] 李云良,張奇,姚靜,等.鄱陽湖湖泊流域系統(tǒng)水文水動力聯(lián)合模擬[J].湖泊科學(xué),2013,25(2):227-335.(LI Yun-liang,ZHANG Qi,YAO Jing,et al.Integrated simulation of hydrological and hydrodynamic processes for Lake Poyang-catchment System [J].Journal of Lake Sciences,2013,25(2):227-335.(in Chinese))

    [16] 楊洵,梁國華,周惠成.基于MIKE11的太子河觀-葠河段水文水動力模型研究[J].水電能源科學(xué),2010,28(11):84-87.(YANG Xun,LIANG Guo-hua,ZHOU Hui-cheng.Study on hydrology and hydrodynamic model in Guanyinge-Shenwo section of Taizihe River based on MIKE11[J].Water Resources and Power,2010,28(11):84-87.(in Chinese))

    如何舔出高潮| 我要看日韩黄色一级片| av女优亚洲男人天堂| 欧洲精品卡2卡3卡4卡5卡区| 欧美+日韩+精品| 老熟妇仑乱视频hdxx| 最好的美女福利视频网| 最近视频中文字幕2019在线8| 日韩亚洲欧美综合| 一本精品99久久精品77| 在线观看66精品国产| 五月伊人婷婷丁香| 久久久久久国产a免费观看| 97热精品久久久久久| 亚洲男人的天堂狠狠| 国产伦精品一区二区三区四那| 精品国产亚洲在线| 99在线人妻在线中文字幕| 国产精华一区二区三区| 国产精品一区二区免费欧美| 悠悠久久av| 香蕉av资源在线| 夜夜爽天天搞| 日韩欧美国产在线观看| 成年女人看的毛片在线观看| 国产国拍精品亚洲av在线观看| 美女cb高潮喷水在线观看| 亚洲男人的天堂狠狠| 国内精品一区二区在线观看| 国内揄拍国产精品人妻在线| 色精品久久人妻99蜜桃| 国产成人a区在线观看| 精品不卡国产一区二区三区| 亚洲av.av天堂| 欧美黑人巨大hd| 99热这里只有精品一区| 欧美日韩综合久久久久久 | 免费在线观看影片大全网站| 亚洲熟妇中文字幕五十中出| 18禁黄网站禁片午夜丰满| 99精品在免费线老司机午夜| 亚洲 欧美 日韩 在线 免费| 国产国拍精品亚洲av在线观看| 久久九九热精品免费| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 女生性感内裤真人,穿戴方法视频| 男人和女人高潮做爰伦理| av在线老鸭窝| .国产精品久久| 少妇丰满av| or卡值多少钱| 精品日产1卡2卡| 真实男女啪啪啪动态图| 成人国产综合亚洲| 又粗又爽又猛毛片免费看| 搡老岳熟女国产| 久久久精品欧美日韩精品| 一边摸一边抽搐一进一小说| 三级毛片av免费| 内地一区二区视频在线| 日日夜夜操网爽| 免费无遮挡裸体视频| 午夜免费成人在线视频| 国产男靠女视频免费网站| 日本免费一区二区三区高清不卡| 亚洲真实伦在线观看| 午夜福利欧美成人| www.色视频.com| 色综合亚洲欧美另类图片| 国产精品一区二区三区四区久久| a级毛片a级免费在线| 欧美激情在线99| 国产淫片久久久久久久久 | 色综合婷婷激情| 伦理电影大哥的女人| 12—13女人毛片做爰片一| 国产精品av视频在线免费观看| 男女之事视频高清在线观看| 日本在线视频免费播放| 一级黄色大片毛片| 久久久久久久久久黄片| 免费在线观看亚洲国产| 少妇高潮的动态图| 我的女老师完整版在线观看| 亚洲成人久久爱视频| 小蜜桃在线观看免费完整版高清| 97人妻精品一区二区三区麻豆| 精品免费久久久久久久清纯| 99久久成人亚洲精品观看| 九九在线视频观看精品| 国产三级黄色录像| 色视频www国产| 天天一区二区日本电影三级| 国产免费av片在线观看野外av| 欧美成人a在线观看| 丰满人妻熟妇乱又伦精品不卡| 亚洲精品一卡2卡三卡4卡5卡| 村上凉子中文字幕在线| 色综合站精品国产| 亚洲成人久久爱视频| 午夜福利成人在线免费观看| 国产在视频线在精品| 国内精品一区二区在线观看| 亚洲国产精品999在线| 97超视频在线观看视频| 美女xxoo啪啪120秒动态图 | 国产精品98久久久久久宅男小说| 欧美日韩黄片免| 久久国产乱子伦精品免费另类| 日韩欧美一区二区三区在线观看| 麻豆国产97在线/欧美| 精品免费久久久久久久清纯| 嫁个100分男人电影在线观看| 美女 人体艺术 gogo| 如何舔出高潮| 亚洲中文日韩欧美视频| 在线天堂最新版资源| 97人妻精品一区二区三区麻豆| 69av精品久久久久久| 亚洲国产精品久久男人天堂| 久久久久久九九精品二区国产| 亚洲无线在线观看| 婷婷丁香在线五月| 日本黄色视频三级网站网址| 噜噜噜噜噜久久久久久91| 国产69精品久久久久777片| 国产毛片a区久久久久| 九九久久精品国产亚洲av麻豆| 十八禁人妻一区二区| 国产色婷婷99| 亚洲精品日韩av片在线观看| 99riav亚洲国产免费| 午夜精品一区二区三区免费看| 中文字幕av成人在线电影| 免费看美女性在线毛片视频| 免费黄网站久久成人精品 | 国内精品久久久久久久电影| 久久精品国产亚洲av涩爱 | 国产欧美日韩精品一区二区| 成人三级黄色视频| 午夜福利视频1000在线观看| 日日摸夜夜添夜夜添av毛片 | 亚洲黑人精品在线| 成人欧美大片| 国内揄拍国产精品人妻在线| 国产美女午夜福利| 亚洲av二区三区四区| 88av欧美| 91麻豆av在线| 久久久久久九九精品二区国产| 少妇裸体淫交视频免费看高清| 精品99又大又爽又粗少妇毛片 | 亚州av有码| 亚洲人与动物交配视频| 亚洲经典国产精华液单 | 久久精品国产亚洲av涩爱 | 国产亚洲精品综合一区在线观看| 国产精品av视频在线免费观看| 一本一本综合久久| 成人国产综合亚洲| АⅤ资源中文在线天堂| 色综合欧美亚洲国产小说| 九九在线视频观看精品| 欧美性猛交黑人性爽| 老司机午夜福利在线观看视频| 久久国产精品影院| 国内毛片毛片毛片毛片毛片| 午夜精品一区二区三区免费看| 人人妻人人澡欧美一区二区| 给我免费播放毛片高清在线观看| 亚洲欧美日韩东京热| 日韩欧美三级三区| 亚洲精品456在线播放app | 亚洲自偷自拍三级| 国产精品亚洲av一区麻豆| h日本视频在线播放| 少妇的逼好多水| 男人和女人高潮做爰伦理| 国产精品永久免费网站| 成人国产一区最新在线观看| 天堂动漫精品| 国产日本99.免费观看| av福利片在线观看| 制服丝袜大香蕉在线| 又爽又黄a免费视频| 好男人电影高清在线观看| 69av精品久久久久久| 亚洲欧美日韩高清在线视频| 男女下面进入的视频免费午夜| 1024手机看黄色片| 亚洲色图av天堂| 国产精品三级大全| 亚洲人成网站在线播| 亚洲第一区二区三区不卡| 18禁在线播放成人免费| 亚洲三级黄色毛片| 他把我摸到了高潮在线观看| 国产69精品久久久久777片| 女同久久另类99精品国产91| 久久精品综合一区二区三区| 国产av一区在线观看免费| 91字幕亚洲| 久久久久久久精品吃奶| 国产精品人妻久久久久久| 久久婷婷人人爽人人干人人爱| 亚洲人成电影免费在线| 小说图片视频综合网站| 一个人免费在线观看的高清视频| 男女那种视频在线观看| 一区福利在线观看| 性欧美人与动物交配| 亚洲人成网站在线播放欧美日韩| 国产高清有码在线观看视频| 中文字幕av在线有码专区| 久久婷婷人人爽人人干人人爱| 国产麻豆成人av免费视频| 亚洲经典国产精华液单 | 少妇高潮的动态图| 国产三级黄色录像| 99热这里只有是精品50| 中文字幕人成人乱码亚洲影| 一级黄片播放器| 欧美最新免费一区二区三区 | 天堂网av新在线| av国产免费在线观看| 久9热在线精品视频| 久久精品国产亚洲av涩爱 | 日本熟妇午夜| 成人特级av手机在线观看| 亚洲人成网站在线播| 又粗又爽又猛毛片免费看| 成人无遮挡网站| 在线观看舔阴道视频| 国产私拍福利视频在线观看| 校园春色视频在线观看| 日韩成人在线观看一区二区三区| 三级毛片av免费| 小说图片视频综合网站| 亚洲aⅴ乱码一区二区在线播放| 啦啦啦观看免费观看视频高清| 久久精品国产自在天天线| av女优亚洲男人天堂| h日本视频在线播放| 国产一区二区三区在线臀色熟女| 亚洲欧美日韩高清在线视频| av欧美777| 午夜激情欧美在线| 色哟哟哟哟哟哟| 亚洲aⅴ乱码一区二区在线播放| 少妇的逼水好多| 91麻豆精品激情在线观看国产| 99热这里只有是精品在线观看 | 99热6这里只有精品| 欧美日韩综合久久久久久 | 久久久久九九精品影院| 国产黄片美女视频| 亚洲,欧美精品.| 99热6这里只有精品| 欧美国产日韩亚洲一区| 两人在一起打扑克的视频| 亚洲午夜理论影院| 又粗又爽又猛毛片免费看| 国产精品精品国产色婷婷| 蜜桃久久精品国产亚洲av| 欧美日韩综合久久久久久 | 1000部很黄的大片| 丰满乱子伦码专区| 99久久精品国产亚洲精品| 小蜜桃在线观看免费完整版高清| 无遮挡黄片免费观看| 免费看美女性在线毛片视频| x7x7x7水蜜桃| 国产伦在线观看视频一区| 搡老岳熟女国产| 亚洲黑人精品在线| 别揉我奶头 嗯啊视频| 99国产精品一区二区三区| 深夜a级毛片| 亚洲欧美日韩卡通动漫| 直男gayav资源| 欧美区成人在线视频| 精品午夜福利在线看| 国产精品美女特级片免费视频播放器| 国产精品人妻久久久久久| 午夜福利欧美成人| 国产亚洲精品久久久com| 人妻丰满熟妇av一区二区三区| 国产乱人视频| 亚洲av第一区精品v没综合| 亚洲天堂国产精品一区在线| 久久久久久九九精品二区国产| 久久久久国内视频| 网址你懂的国产日韩在线| 99国产综合亚洲精品| 国产精品野战在线观看| 国产高清视频在线观看网站| 久久人妻av系列| 草草在线视频免费看| 国产精品亚洲av一区麻豆| 在线观看av片永久免费下载| 麻豆国产av国片精品| 在线观看66精品国产| 岛国在线免费视频观看| 欧美日韩瑟瑟在线播放| 日日干狠狠操夜夜爽| 国产精品人妻久久久久久| 床上黄色一级片| 国产精品一区二区免费欧美| 亚州av有码| 欧美国产日韩亚洲一区| 久久久成人免费电影| 如何舔出高潮| 国产不卡一卡二| 国产亚洲欧美98| 成人国产一区最新在线观看| 热99re8久久精品国产| 999久久久精品免费观看国产| 在线观看午夜福利视频| 欧美激情久久久久久爽电影| 久久久久久久精品吃奶| 欧美不卡视频在线免费观看| 搡老妇女老女人老熟妇| 亚洲欧美日韩东京热| 日韩欧美免费精品| 欧美日韩国产亚洲二区| 极品教师在线免费播放| 国产伦一二天堂av在线观看| 国产精品精品国产色婷婷| 波多野结衣巨乳人妻| 亚洲精品粉嫩美女一区| 99精品在免费线老司机午夜| 国产成年人精品一区二区| 精品人妻一区二区三区麻豆 | 在线a可以看的网站| 亚洲黑人精品在线| 国产精品一区二区三区四区免费观看 | 丁香六月欧美| 国产精品三级大全| 成人精品一区二区免费| 久久久精品欧美日韩精品| 小蜜桃在线观看免费完整版高清| 亚洲av电影在线进入| 看片在线看免费视频| 天堂√8在线中文| 啦啦啦观看免费观看视频高清| 亚洲av美国av| 国产在线精品亚洲第一网站| 长腿黑丝高跟| 波多野结衣巨乳人妻| 国产一区二区激情短视频| 免费观看精品视频网站| 一二三四社区在线视频社区8| 久久久久久久亚洲中文字幕 | 搡老熟女国产l中国老女人| 国产精品爽爽va在线观看网站| 真人做人爱边吃奶动态| 亚洲av电影在线进入| 久久久久久久久中文| 成人特级av手机在线观看| 国产在线男女| 嫩草影院精品99| 午夜视频国产福利| 99国产综合亚洲精品| 搡老熟女国产l中国老女人| 最近视频中文字幕2019在线8| 99久久99久久久精品蜜桃| 中文字幕熟女人妻在线| 亚洲一区二区三区不卡视频| 一区二区三区四区激情视频 | 男人和女人高潮做爰伦理| 欧美日韩黄片免| 国产精品嫩草影院av在线观看 | 日韩成人在线观看一区二区三区| 老司机午夜福利在线观看视频| 97人妻精品一区二区三区麻豆| 亚洲精品456在线播放app | 啪啪无遮挡十八禁网站| 欧美日韩综合久久久久久 | 日本撒尿小便嘘嘘汇集6| 婷婷精品国产亚洲av在线| 午夜激情福利司机影院| 深夜a级毛片| 99久国产av精品| 色噜噜av男人的天堂激情| 精品日产1卡2卡| 一区福利在线观看| 国产aⅴ精品一区二区三区波| 久久久久久久久久黄片| 99热精品在线国产| 舔av片在线| 欧美+日韩+精品| 亚洲七黄色美女视频| 欧美精品啪啪一区二区三区| 午夜福利在线观看吧| 麻豆成人av在线观看| 日韩高清综合在线| 国产免费一级a男人的天堂| 久久久国产成人免费| 成人毛片a级毛片在线播放| 一本一本综合久久| 9191精品国产免费久久| 18禁在线播放成人免费| 黄色一级大片看看| 久久午夜福利片| 一区二区三区激情视频| 国产美女午夜福利| 日韩欧美精品v在线| 免费看日本二区| eeuss影院久久| 亚洲人成网站在线播| 精品一区二区免费观看| 看免费av毛片| 午夜福利免费观看在线| 日本熟妇午夜| www.www免费av| 日本黄色视频三级网站网址| 精品人妻熟女av久视频| 村上凉子中文字幕在线| 亚洲精品乱码久久久v下载方式| 欧美午夜高清在线| 夜夜躁狠狠躁天天躁| 精品久久国产蜜桃| 欧美+日韩+精品| 久久久久久久精品吃奶| 欧美三级亚洲精品| 亚洲人与动物交配视频| 美女大奶头视频| 亚洲美女搞黄在线观看 | 在线观看一区二区三区| 久久国产精品人妻蜜桃| 欧美日韩中文字幕国产精品一区二区三区| 熟妇人妻久久中文字幕3abv| 一区二区三区免费毛片| 国产高潮美女av| 97超级碰碰碰精品色视频在线观看| 亚洲乱码一区二区免费版| 精品午夜福利视频在线观看一区| 一进一出好大好爽视频| 99热精品在线国产| 少妇裸体淫交视频免费看高清| 白带黄色成豆腐渣| 亚洲av免费高清在线观看| 精品人妻偷拍中文字幕| 免费看日本二区| 黄色女人牲交| 日韩大尺度精品在线看网址| 嫩草影院精品99| 黄色配什么色好看| 国产成人影院久久av| 成人av一区二区三区在线看| 国产69精品久久久久777片| 国内揄拍国产精品人妻在线| 99热这里只有精品一区| 男女做爰动态图高潮gif福利片| 最近最新中文字幕大全电影3| 日韩欧美三级三区| 99热这里只有是精品50| 宅男免费午夜| www.熟女人妻精品国产| 午夜福利欧美成人| 一级av片app| 亚洲精品在线观看二区| 好男人在线观看高清免费视频| 欧美区成人在线视频| 欧美成人a在线观看| 欧美激情久久久久久爽电影| 制服丝袜大香蕉在线| 在线观看舔阴道视频| 成人特级av手机在线观看| 欧美潮喷喷水| 精品午夜福利在线看| 国产aⅴ精品一区二区三区波| 国产精品av视频在线免费观看| 亚洲av不卡在线观看| 久久国产乱子伦精品免费另类| 欧美色欧美亚洲另类二区| 久久久久久九九精品二区国产| 好男人在线观看高清免费视频| 成人特级黄色片久久久久久久| 最近最新免费中文字幕在线| 成年版毛片免费区| 免费搜索国产男女视频| 亚洲欧美日韩东京热| 69人妻影院| 伊人久久精品亚洲午夜| 日韩欧美免费精品| 伊人久久精品亚洲午夜| 一个人免费在线观看电影| 国产伦在线观看视频一区| 99国产综合亚洲精品| 欧美+日韩+精品| 欧美在线一区亚洲| 国语自产精品视频在线第100页| 欧美黑人巨大hd| 女同久久另类99精品国产91| 一本久久中文字幕| 国产午夜福利久久久久久| 成年人黄色毛片网站| av在线蜜桃| 国产真实乱freesex| 老熟妇仑乱视频hdxx| 51国产日韩欧美| 国产乱人视频| 性欧美人与动物交配| 国产午夜精品久久久久久一区二区三区 | av专区在线播放| 男女床上黄色一级片免费看| 欧美潮喷喷水| 亚洲狠狠婷婷综合久久图片| 午夜日韩欧美国产| 一个人免费在线观看的高清视频| 啦啦啦观看免费观看视频高清| 美女免费视频网站| 亚洲人成伊人成综合网2020| 久久人人精品亚洲av| 午夜老司机福利剧场| 一本综合久久免费| 乱人视频在线观看| 国产精品人妻久久久久久| 久久久久久大精品| 亚洲成av人片在线播放无| 我要看日韩黄色一级片| 搡老熟女国产l中国老女人| 极品教师在线免费播放| 搞女人的毛片| 一个人免费在线观看电影| 日本熟妇午夜| 一区二区三区四区激情视频 | 久久人人精品亚洲av| 在线看三级毛片| av在线观看视频网站免费| 女人被狂操c到高潮| 自拍偷自拍亚洲精品老妇| 两人在一起打扑克的视频| 亚洲国产精品999在线| 中出人妻视频一区二区| 日韩中字成人| 亚洲av电影不卡..在线观看| 最新中文字幕久久久久| 欧美国产日韩亚洲一区| 亚州av有码| 啦啦啦观看免费观看视频高清| 男女床上黄色一级片免费看| 成人欧美大片| 亚洲狠狠婷婷综合久久图片| 国产伦在线观看视频一区| а√天堂www在线а√下载| 黄色配什么色好看| 久久天躁狠狠躁夜夜2o2o| 小说图片视频综合网站| 亚洲成av人片免费观看| 精品人妻偷拍中文字幕| xxxwww97欧美| 成人三级黄色视频| 日本三级黄在线观看| 丰满人妻一区二区三区视频av| 久久久精品大字幕| 丰满乱子伦码专区| 在线a可以看的网站| 99热这里只有是精品在线观看 | 9191精品国产免费久久| 男插女下体视频免费在线播放| 亚洲成人精品中文字幕电影| 国内精品久久久久精免费| 久久久久久久精品吃奶| 又粗又爽又猛毛片免费看| 级片在线观看| 一个人免费在线观看电影| 亚洲成人久久性| 欧洲精品卡2卡3卡4卡5卡区| 亚洲狠狠婷婷综合久久图片| 亚洲电影在线观看av| 一级av片app| 国产黄a三级三级三级人| 黄色日韩在线| 日本撒尿小便嘘嘘汇集6| 国产成年人精品一区二区| 757午夜福利合集在线观看| 日韩 亚洲 欧美在线| 香蕉av资源在线| 特级一级黄色大片| 成人性生交大片免费视频hd| 黄色女人牲交| 成人永久免费在线观看视频| av女优亚洲男人天堂| 国产蜜桃级精品一区二区三区| 久久久久免费精品人妻一区二区| 久久久色成人| 精品国产三级普通话版| 国内精品久久久久精免费| 啪啪无遮挡十八禁网站| 精品人妻1区二区| 欧美一级a爱片免费观看看| 亚洲国产精品成人综合色| 此物有八面人人有两片| 日本一本二区三区精品| 国内毛片毛片毛片毛片毛片| 日韩欧美 国产精品| 国产av一区在线观看免费| 观看美女的网站| 精品人妻偷拍中文字幕| 国产aⅴ精品一区二区三区波| 国产高清有码在线观看视频| 亚洲av成人不卡在线观看播放网| 国产又黄又爽又无遮挡在线| 国产黄色小视频在线观看| 大型黄色视频在线免费观看| 观看免费一级毛片| 69av精品久久久久久| av中文乱码字幕在线| 看免费av毛片| 在线观看午夜福利视频| 五月伊人婷婷丁香| 99久久99久久久精品蜜桃| 日韩亚洲欧美综合|