• <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欧美777| av欧美777| 亚洲熟妇中文字幕五十中出| 嫩草影院精品99| 老汉色av国产亚洲站长工具| av片东京热男人的天堂| 欧美黑人巨大hd| 男人舔女人下体高潮全视频| 欧美+亚洲+日韩+国产| 99在线视频只有这里精品首页| 欧美性猛交黑人性爽| 草草在线视频免费看| 日韩欧美一区二区三区在线观看| 欧美黄色片欧美黄色片| 国产精品九九99| 日韩欧美免费精品| 母亲3免费完整高清在线观看| 国产成人影院久久av| 大香蕉久久成人网| 1024视频免费在线观看| 麻豆av在线久日| 老鸭窝网址在线观看| 女生性感内裤真人,穿戴方法视频| 窝窝影院91人妻| 亚洲欧美精品综合久久99| 国产又黄又爽又无遮挡在线| 日本一本二区三区精品| 亚洲第一欧美日韩一区二区三区| 色在线成人网| 每晚都被弄得嗷嗷叫到高潮| 999久久久国产精品视频| 国产亚洲精品av在线| 国产精品久久久av美女十八| 999精品在线视频| 国产精品影院久久| 日韩精品免费视频一区二区三区| 久久婷婷成人综合色麻豆| 久久这里只有精品19| av福利片在线| 日本a在线网址| 久久人妻福利社区极品人妻图片| 夜夜爽天天搞| 成人18禁在线播放| 变态另类成人亚洲欧美熟女| 久久中文看片网| 俺也久久电影网| 在线播放国产精品三级| 一边摸一边做爽爽视频免费| 一边摸一边做爽爽视频免费| 成人特级黄色片久久久久久久| 久久香蕉精品热| 亚洲精品一卡2卡三卡4卡5卡| 久久欧美精品欧美久久欧美| а√天堂www在线а√下载| 欧美黄色淫秽网站| 亚洲欧美精品综合一区二区三区| 欧美黑人欧美精品刺激| 非洲黑人性xxxx精品又粗又长| 精品国内亚洲2022精品成人| 好看av亚洲va欧美ⅴa在| 欧美绝顶高潮抽搐喷水| 91九色精品人成在线观看| 日韩国内少妇激情av| 精品一区二区三区视频在线观看免费| 美女高潮到喷水免费观看| 欧美乱妇无乱码| 精品一区二区三区av网在线观看| 嫩草影视91久久| 亚洲久久久国产精品| 51午夜福利影视在线观看| 久久精品成人免费网站| 两性夫妻黄色片| 女人被狂操c到高潮| 国产成人啪精品午夜网站| 欧美久久黑人一区二区| 曰老女人黄片| 久久婷婷人人爽人人干人人爱| 夜夜爽天天搞| 亚洲欧美日韩高清在线视频| 亚洲专区中文字幕在线| 久久亚洲精品不卡| 变态另类丝袜制服| 天天添夜夜摸| 亚洲黑人精品在线| 亚洲avbb在线观看| 亚洲 欧美一区二区三区| 亚洲九九香蕉| 大型黄色视频在线免费观看| 亚洲av成人不卡在线观看播放网| 91麻豆精品激情在线观看国产| 日韩高清综合在线| 亚洲精华国产精华精| 日韩一卡2卡3卡4卡2021年| 脱女人内裤的视频| 神马国产精品三级电影在线观看 | av有码第一页| 搡老岳熟女国产| 亚洲中文日韩欧美视频| 日韩国内少妇激情av| 色哟哟哟哟哟哟| 国产精品久久久av美女十八| 精品国产美女av久久久久小说| 天天躁夜夜躁狠狠躁躁| 一级a爱片免费观看的视频| www日本黄色视频网| 久久天堂一区二区三区四区| 亚洲精品粉嫩美女一区| АⅤ资源中文在线天堂| 在线观看免费日韩欧美大片| www.精华液| 日韩大尺度精品在线看网址| 欧美 亚洲 国产 日韩一| 丝袜在线中文字幕| 久热这里只有精品99| а√天堂www在线а√下载| 欧美色欧美亚洲另类二区| 高潮久久久久久久久久久不卡| 激情在线观看视频在线高清| 日韩免费av在线播放| 国产高清激情床上av| 亚洲七黄色美女视频| 欧美在线一区亚洲| 麻豆成人午夜福利视频| 国产av一区在线观看免费| 99久久99久久久精品蜜桃| 亚洲精品在线美女| 国产成人欧美| 91成人精品电影| 俄罗斯特黄特色一大片| 在线av久久热| 午夜免费鲁丝| 欧美激情久久久久久爽电影| 国产精品久久久人人做人人爽| 午夜福利在线观看吧| 我的亚洲天堂| 欧美日本亚洲视频在线播放| 日本a在线网址| 国产人伦9x9x在线观看| 久久久久亚洲av毛片大全| 制服丝袜大香蕉在线| 午夜福利成人在线免费观看| 久久精品91蜜桃| 黄色视频不卡| 午夜久久久在线观看| 久久精品91蜜桃| 国产精品爽爽va在线观看网站 | 国产亚洲精品综合一区在线观看 | av视频在线观看入口| 制服人妻中文乱码| 色在线成人网| 1024手机看黄色片| 欧美乱码精品一区二区三区| 18禁国产床啪视频网站| 免费在线观看完整版高清| a级毛片在线看网站| or卡值多少钱| 亚洲熟妇中文字幕五十中出| 久久九九热精品免费| 极品教师在线免费播放| 亚洲人成伊人成综合网2020| 久久久久免费精品人妻一区二区 | 天堂√8在线中文| 别揉我奶头~嗯~啊~动态视频| 久久午夜综合久久蜜桃| 亚洲av中文字字幕乱码综合 | 变态另类成人亚洲欧美熟女| 男男h啪啪无遮挡| 国产av一区在线观看免费| 很黄的视频免费| 国产熟女午夜一区二区三区| 长腿黑丝高跟| 国产精品一区二区精品视频观看| 免费看十八禁软件| 巨乳人妻的诱惑在线观看| 最新美女视频免费是黄的| av有码第一页| 久久久久久免费高清国产稀缺| www国产在线视频色| 在线观看日韩欧美| 国产精华一区二区三区| 国产成人精品久久二区二区91| 亚洲人成77777在线视频| 欧美黄色片欧美黄色片| 亚洲中文字幕日韩| 黑人欧美特级aaaaaa片| 18禁裸乳无遮挡免费网站照片 | 久久久久久久久免费视频了| 一边摸一边抽搐一进一小说| 宅男免费午夜| 很黄的视频免费| 亚洲av成人一区二区三| 大型av网站在线播放| 美女免费视频网站| 亚洲真实伦在线观看| 老熟妇仑乱视频hdxx| 成人av一区二区三区在线看| netflix在线观看网站| 久久天躁狠狠躁夜夜2o2o| 国产精品综合久久久久久久免费| 亚洲成av人片免费观看| 波多野结衣高清作品| 男人舔女人的私密视频| 国产免费男女视频| 国产激情欧美一区二区| 国产精品亚洲av一区麻豆| 亚洲成av人片免费观看| 国产成人精品无人区| 在线天堂中文资源库| 久久婷婷人人爽人人干人人爱| 好男人电影高清在线观看| av天堂在线播放| 国产亚洲欧美精品永久| 久久中文字幕一级| 久久久久亚洲av毛片大全| 国产aⅴ精品一区二区三区波| 国产欧美日韩一区二区精品| 香蕉丝袜av| 99热只有精品国产| 一区二区日韩欧美中文字幕| 老熟妇仑乱视频hdxx| 亚洲五月色婷婷综合| 久久久久免费精品人妻一区二区 | 国产精品日韩av在线免费观看| 在线av久久热| 成年人黄色毛片网站| 久久久久久国产a免费观看| 亚洲精品久久成人aⅴ小说| 国产精品香港三级国产av潘金莲| 无限看片的www在线观看| 一夜夜www| 久久国产精品人妻蜜桃| 人妻久久中文字幕网| 久久草成人影院| 一级a爱视频在线免费观看| 日韩欧美国产在线观看| 亚洲午夜精品一区,二区,三区| 男人舔奶头视频| 久久精品影院6| tocl精华| 天堂影院成人在线观看| 午夜两性在线视频| 国产精品久久久久久精品电影 | 亚洲黑人精品在线| 久久久久久久午夜电影| a级毛片a级免费在线| 久久久国产精品麻豆| 哪里可以看免费的av片| 男女那种视频在线观看| 人人澡人人妻人| 搞女人的毛片| 色老头精品视频在线观看| 夜夜看夜夜爽夜夜摸| 91成人精品电影| 男女做爰动态图高潮gif福利片| 国产av不卡久久| 日本精品一区二区三区蜜桃| 国产精品久久久久久亚洲av鲁大| 黄色 视频免费看| www日本黄色视频网| 日韩欧美国产一区二区入口| 神马国产精品三级电影在线观看 | 精品国产国语对白av| 一级黄色大片毛片| 一级a爱片免费观看的视频| 国产成年人精品一区二区| 日韩欧美免费精品| 1024视频免费在线观看| 久久天躁狠狠躁夜夜2o2o| 两人在一起打扑克的视频| 日日爽夜夜爽网站| 999久久久精品免费观看国产| 久热这里只有精品99| 亚洲aⅴ乱码一区二区在线播放 | 欧美中文综合在线视频| 国产精品亚洲美女久久久| 伊人久久大香线蕉亚洲五| 啦啦啦 在线观看视频| 国产真实乱freesex| 精品国产乱子伦一区二区三区| 黑人巨大精品欧美一区二区mp4| 国产三级黄色录像| 精品久久久久久久久久久久久 | av福利片在线| 亚洲精品中文字幕在线视频| 母亲3免费完整高清在线观看| 9191精品国产免费久久| 久久香蕉国产精品| 国产又黄又爽又无遮挡在线| 久久狼人影院| 国产亚洲欧美98| 久久精品91蜜桃| 变态另类丝袜制服| 欧美黑人巨大hd| 久99久视频精品免费| 69av精品久久久久久| av福利片在线| 99久久99久久久精品蜜桃| 国产爱豆传媒在线观看 | 国产不卡一卡二| 国产精品av久久久久免费| 午夜激情福利司机影院| 丁香六月欧美| 日韩免费av在线播放| 国产精品日韩av在线免费观看| 麻豆国产av国片精品| 国产精品 欧美亚洲| 久久精品国产综合久久久| 999久久久精品免费观看国产| 自线自在国产av| 国产一卡二卡三卡精品| 中文字幕人妻丝袜一区二区| 亚洲av成人av| 亚洲三区欧美一区| 精品日产1卡2卡| 黄色丝袜av网址大全| 国产伦在线观看视频一区| 黄色 视频免费看| 久久中文看片网| 久久久久久久精品吃奶| 中文字幕最新亚洲高清| 久久精品aⅴ一区二区三区四区| 老司机午夜十八禁免费视频| ponron亚洲| 高清毛片免费观看视频网站| 国产精品99久久99久久久不卡| 日韩欧美免费精品| 国产成人精品久久二区二区免费| 19禁男女啪啪无遮挡网站| 国产av一区在线观看免费| 一a级毛片在线观看| 桃色一区二区三区在线观看| 亚洲精品在线观看二区| 国产极品粉嫩免费观看在线| 老熟妇仑乱视频hdxx| 男人的好看免费观看在线视频 | 变态另类成人亚洲欧美熟女| 天堂√8在线中文| 日本撒尿小便嘘嘘汇集6| 麻豆久久精品国产亚洲av| 丝袜人妻中文字幕| 国产精品免费视频内射| 18禁裸乳无遮挡免费网站照片 | 少妇裸体淫交视频免费看高清 | 亚洲色图av天堂| 免费观看人在逋| 亚洲天堂国产精品一区在线| 97人妻精品一区二区三区麻豆 | 韩国精品一区二区三区| 国产精品久久久av美女十八| 母亲3免费完整高清在线观看| 日本三级黄在线观看| 18禁观看日本| АⅤ资源中文在线天堂| 欧美黄色淫秽网站| 手机成人av网站| 国产精品久久久av美女十八| 午夜福利在线观看吧| √禁漫天堂资源中文www| 亚洲精品国产区一区二| 69av精品久久久久久| 超碰成人久久| 两个人视频免费观看高清| 亚洲自偷自拍图片 自拍| 叶爱在线成人免费视频播放| 熟妇人妻久久中文字幕3abv| 日韩欧美免费精品| 超碰成人久久| 69av精品久久久久久| 757午夜福利合集在线观看| 两个人视频免费观看高清| 日韩国内少妇激情av| 成在线人永久免费视频| 久久久久久久久久黄片| 可以免费在线观看a视频的电影网站| 久久久久亚洲av毛片大全| 午夜老司机福利片| 日韩欧美国产一区二区入口| 十分钟在线观看高清视频www| 亚洲专区中文字幕在线| 午夜精品久久久久久毛片777| 91九色精品人成在线观看| avwww免费| tocl精华| 精品久久久久久久久久久久久 | 很黄的视频免费| 欧美乱码精品一区二区三区| 久久久国产欧美日韩av| 1024视频免费在线观看| 女性被躁到高潮视频| 欧美激情 高清一区二区三区| 99在线视频只有这里精品首页| 亚洲国产日韩欧美精品在线观看 | 亚洲av美国av| 精品欧美国产一区二区三| 麻豆一二三区av精品| 国产av在哪里看| 国产精品久久久久久人妻精品电影| 国产成人精品久久二区二区免费| 免费观看人在逋| 欧美绝顶高潮抽搐喷水| 亚洲精品中文字幕一二三四区| 亚洲va日本ⅴa欧美va伊人久久| 日本a在线网址| 少妇熟女aⅴ在线视频| 日本a在线网址| 又黄又爽又免费观看的视频| 久久亚洲真实| 黑人巨大精品欧美一区二区mp4| 禁无遮挡网站| 久久久久久九九精品二区国产 | 热re99久久国产66热| 亚洲av日韩精品久久久久久密| 中文资源天堂在线| 听说在线观看完整版免费高清| 亚洲一区二区三区不卡视频| 亚洲av第一区精品v没综合| 亚洲午夜精品一区,二区,三区| 免费搜索国产男女视频| 中文亚洲av片在线观看爽| 青草久久国产| 不卡一级毛片| 日韩有码中文字幕| 国产亚洲精品第一综合不卡| 欧美一级毛片孕妇| 日本 av在线| 成人18禁高潮啪啪吃奶动态图| 亚洲成av片中文字幕在线观看| 九色国产91popny在线| 变态另类丝袜制服| 午夜日韩欧美国产| 香蕉丝袜av| 女性被躁到高潮视频| 亚洲国产欧美一区二区综合| 亚洲欧美精品综合一区二区三区| 1024视频免费在线观看| 一区二区三区高清视频在线| 操出白浆在线播放| av中文乱码字幕在线| 少妇的丰满在线观看| 日韩成人在线观看一区二区三区| 性欧美人与动物交配| 无遮挡黄片免费观看| 法律面前人人平等表现在哪些方面| 男女做爰动态图高潮gif福利片| 国内久久婷婷六月综合欲色啪| 亚洲第一青青草原| 国产精品野战在线观看| 听说在线观看完整版免费高清| 不卡av一区二区三区| 久久亚洲真实| 亚洲国产高清在线一区二区三 | 国产精品1区2区在线观看.| 亚洲精华国产精华精| 精品国产亚洲在线| 国产日本99.免费观看| 亚洲自拍偷在线| 欧美成人免费av一区二区三区| 一区二区日韩欧美中文字幕| 欧美一级a爱片免费观看看 | 国产成人精品久久二区二区91| 真人做人爱边吃奶动态| 日韩视频一区二区在线观看| 搡老熟女国产l中国老女人| 19禁男女啪啪无遮挡网站| 色尼玛亚洲综合影院| 身体一侧抽搐| 国产亚洲欧美在线一区二区| 丰满人妻熟妇乱又伦精品不卡| 国产男靠女视频免费网站| 99久久无色码亚洲精品果冻| 9191精品国产免费久久| 哪里可以看免费的av片| 久久人妻av系列| 香蕉av资源在线| 色综合欧美亚洲国产小说| 亚洲全国av大片| 两性夫妻黄色片| 国产精品1区2区在线观看.| 性色av乱码一区二区三区2| 波多野结衣巨乳人妻| 久久国产精品男人的天堂亚洲| 人人妻人人澡人人看| 国产午夜福利久久久久久| 日本 欧美在线| 一本大道久久a久久精品| 变态另类丝袜制服| 一本精品99久久精品77| 久久精品人妻少妇| 在线永久观看黄色视频| 亚洲人成77777在线视频| 国产av一区在线观看免费| 亚洲精品中文字幕在线视频| 亚洲精品美女久久久久99蜜臀| 亚洲色图 男人天堂 中文字幕| 国产成人精品久久二区二区免费| 宅男免费午夜| 久久久精品国产亚洲av高清涩受| 中文字幕高清在线视频| 国产成人系列免费观看| 男女视频在线观看网站免费 | 18禁国产床啪视频网站| 国产视频内射| 老汉色∧v一级毛片| 欧美性猛交黑人性爽| 欧洲精品卡2卡3卡4卡5卡区| 好男人电影高清在线观看| 人成视频在线观看免费观看| 中文字幕av电影在线播放| 婷婷六月久久综合丁香| 亚洲欧美激情综合另类| 亚洲成人久久性| 久久亚洲真实| 亚洲一区二区三区不卡视频| 亚洲精品粉嫩美女一区| 午夜久久久在线观看| 美女大奶头视频| www国产在线视频色| 1024香蕉在线观看| 国产精品一区二区免费欧美| 久久香蕉国产精品| 女人爽到高潮嗷嗷叫在线视频| 久久中文字幕一级| 18禁观看日本| 亚洲专区字幕在线| 亚洲成av片中文字幕在线观看| 丰满人妻熟妇乱又伦精品不卡| 国产主播在线观看一区二区| svipshipincom国产片| 国产蜜桃级精品一区二区三区| 成年女人毛片免费观看观看9| 亚洲精品一区av在线观看| 99热6这里只有精品| 免费在线观看影片大全网站| 欧美zozozo另类| 人人妻人人看人人澡| 嫩草影院精品99| 精品国产国语对白av| 精品国产乱码久久久久久男人| 亚洲自偷自拍图片 自拍| av视频在线观看入口| 老司机在亚洲福利影院| 亚洲五月天丁香| 侵犯人妻中文字幕一二三四区| 最近最新免费中文字幕在线| 亚洲电影在线观看av| 桃红色精品国产亚洲av| 亚洲专区中文字幕在线| 欧美国产日韩亚洲一区| 女生性感内裤真人,穿戴方法视频| 女性被躁到高潮视频| 最近在线观看免费完整版| 美女午夜性视频免费| 国产精品一区二区免费欧美| 少妇裸体淫交视频免费看高清 | 亚洲av五月六月丁香网| 亚洲成人久久爱视频| 男女下面进入的视频免费午夜 | 亚洲成人国产一区在线观看| 人人妻人人看人人澡| 久久精品91无色码中文字幕| 女人高潮潮喷娇喘18禁视频| 亚洲第一电影网av| 日韩视频一区二区在线观看| 亚洲人成网站在线播放欧美日韩| 欧美色欧美亚洲另类二区| 香蕉av资源在线| 国产伦在线观看视频一区| 国产av又大| 一区福利在线观看| 真人做人爱边吃奶动态| www.www免费av| 亚洲国产毛片av蜜桃av| 国产精品电影一区二区三区| 亚洲国产精品成人综合色| 日韩精品青青久久久久久| 在线播放国产精品三级| 91成年电影在线观看| 在线观看舔阴道视频| 大型av网站在线播放| 国产精品免费视频内射| 1024视频免费在线观看| 亚洲五月色婷婷综合| 一区二区三区精品91| 妹子高潮喷水视频| 亚洲色图 男人天堂 中文字幕| 老熟妇仑乱视频hdxx| 免费在线观看黄色视频的| 男女做爰动态图高潮gif福利片| 90打野战视频偷拍视频| 91在线观看av| 国产精品 国内视频| 国产一级毛片七仙女欲春2 | 欧美乱码精品一区二区三区| 日韩欧美免费精品| 亚洲色图av天堂| 久久精品成人免费网站| 一二三四在线观看免费中文在| 午夜免费激情av| 日本免费一区二区三区高清不卡| 日韩一卡2卡3卡4卡2021年| 97人妻精品一区二区三区麻豆 | 久久国产精品影院| 一本一本综合久久| 伊人久久大香线蕉亚洲五| 在线av久久热| 国产高清激情床上av| 高潮久久久久久久久久久不卡| 国产午夜精品久久久久久| 好看av亚洲va欧美ⅴa在| 精品免费久久久久久久清纯| 在线观看免费视频日本深夜| 99在线视频只有这里精品首页|