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

    八字腦灘群航道整治數(shù)學(xué)模型計(jì)算與分析

    2016-05-25 00:37:00倪志輝吳立春
    關(guān)鍵詞:八字河段淤積

    倪志輝,郭 毅,吳立春

    (1. 重慶交通大學(xué) 水利水運(yùn)工程教育部重點(diǎn)實(shí)驗(yàn)室 國(guó)家內(nèi)河航道整治工程技術(shù)研究中心,重慶 400074;2. 重慶交通大學(xué) 西南水運(yùn)工程科學(xué)研究所,重慶 400016;3. 重慶第二師范學(xué)院,重慶 400067 )

    八字腦灘群航道整治數(shù)學(xué)模型計(jì)算與分析

    倪志輝1,2,郭 毅1,吳立春3

    (1. 重慶交通大學(xué) 水利水運(yùn)工程教育部重點(diǎn)實(shí)驗(yàn)室 國(guó)家內(nèi)河航道整治工程技術(shù)研究中心,重慶 400074;2. 重慶交通大學(xué) 西南水運(yùn)工程科學(xué)研究所,重慶 400016;3. 重慶第二師范學(xué)院,重慶 400067 )

    嘉陵江八字腦灘群位于鳳儀航電樞紐庫(kù)尾,由于馬回電站發(fā)電泄水口與船閘出口分散布置,導(dǎo)致八字腦灘段右汊槽窄水淺,左汊又極其彎曲和狹窄,不具備通航條件;同時(shí),大石塘灘段的沙卵石江心洲將嘉陵江水分為左右兩汊,左汊較順直,但枯期水深不足,右汊彎曲半徑較小,不能滿足高等級(jí)航道建設(shè)標(biāo)準(zhǔn)。針對(duì)嘉陵江八字腦灘群河段的礙航特征,應(yīng)用數(shù)學(xué)模型對(duì)兩種不同的整治方案進(jìn)行水流條件以及泥沙淤積情況計(jì)算分析與比較。研究結(jié)果表明:該數(shù)學(xué)模型可較好地模擬河段水流的運(yùn)動(dòng)及泥沙的淤積情況;從整治效果看,方案一相對(duì)較優(yōu),泥沙回淤量較少,航槽相對(duì)穩(wěn)定,可作為該灘整治方案,供設(shè)計(jì)參考。

    航道工程;數(shù)學(xué)模型;泥沙淤積;灘群;通航條件

    航道整治是河道治理的一個(gè)部分,是一種用整治建筑物調(diào)整和控制水流,穩(wěn)定有利河勢(shì),以改善航道航行條件的工程措施,其中使用最廣泛的是:丁壩、潛壩、潛堤和挖槽等方式[1]。中國(guó)的航道整治有著悠久的歷史。早在1565年,潘季馴就提出了“以堤束水,以水攻沙”的整治原則。21世紀(jì)以來(lái),隨著水運(yùn)事業(yè)的發(fā)展,中國(guó)數(shù)量眾多的復(fù)雜灘群都進(jìn)行了航道整治工程。例如,曹民雄等[2]在對(duì)長(zhǎng)江上游和尚巖灘群灘性分析的基礎(chǔ)上,采取了1∶400的正態(tài)模型研究整治方案,并取得較好效果;楊禮生等[3]針對(duì)復(fù)雜的泉港灘群航段,開發(fā)了河工模型,對(duì)該河段航道整治工程進(jìn)行了方案優(yōu)化;費(fèi)曉昕等[4]建立了河工模型,用于研究整治樞紐壩下變動(dòng)回水區(qū)的航道,并提出合理可行的整治方案。同樣,國(guó)內(nèi)外其他學(xué)者[5-12]在航道整治方面也做出了許多研究。

    嘉陵江梯級(jí)開發(fā)從20世紀(jì)80年代末蓬安馬回梯級(jí)開工建設(shè)起,歷時(shí)20多年,目前已接近尾聲。由于開發(fā)建設(shè)時(shí)間長(zhǎng),各梯級(jí)投資主體多樣,再加上各時(shí)期開發(fā)項(xiàng)目的側(cè)重點(diǎn)不同,特別是早期開發(fā)的梯級(jí)將航運(yùn)擺在較次要的位置,便形成了各航運(yùn)梯級(jí)以正常擋水位銜接的格局。嘉陵江八字腦灘群位于鳳儀樞紐的回水末端,馬回樞紐船閘下游。電站泄水發(fā)電、水位消落時(shí),在梯級(jí)庫(kù)尾段便有部分航段的部分灘險(xiǎn)礙航。因此筆者通過建立二維水沙數(shù)學(xué)模型模擬該灘群附近航道情況,并對(duì)部分礙航灘險(xiǎn)進(jìn)行整治,以滿足航運(yùn)需要。

    1 河段礙航機(jī)理分析

    嘉陵江流域地理位置為東經(jīng)102°~109°和北緯29°40′~34°30′,地勢(shì)由西北向東南傾斜,地形復(fù)雜。八字腦灘群(含八字腦、大石塘、滿天星)位于鳳儀航電樞紐庫(kù)尾,該灘群上鄰馬回樞紐船閘出口,灘段總長(zhǎng)約5.4 km。河段徑流主要來(lái)源于降水,年內(nèi)變化與降水基本一致,年內(nèi)年際變化均較大。根據(jù)金溪水文站實(shí)測(cè)資料統(tǒng)計(jì)分析,本河段汛期為5—10月,占全年水量的79.8%,其中7—9月占全年水量的51.9%,非汛期11月—翌年4月僅占全年水量的20.2%。

    圖1 八字腦灘群河段河勢(shì)Fig.1 Diagram of Bazinao Shoal Group section regime

    八字腦灘群平面形態(tài)十分復(fù)雜,上段(八字腦段)寬窄相間,最寬處達(dá)到約820 m,最窄處僅250 m左右,形成1號(hào)和2號(hào)沙卵石江心洲。上游1號(hào)江心洲長(zhǎng)約350 m,寬約430 m,高程在279.5~285.3 m之間;下游2號(hào)江心洲長(zhǎng)約1 200 m,寬約360 m,高程在279.5~282.2 m之間。水流被兩個(gè)江心洲切割分離成多股,形成了“周家浩”等多個(gè)支汊。由于船閘出口即河道右汊——周家浩槽窄水淺,加之馬回電站發(fā)電泄水口與船閘出口分散布置(電站泄水在上游,船閘出口在下游),造成主流棄右浩而沿左汊下泄,而左汊又極其彎曲和狹窄,根本不具備通航條件,使該灘枯水期船舶無(wú)法正常航行,幾乎到了斷航的境地。

    大石塘灘位于3號(hào)江心洲上端,該洲長(zhǎng)約1 900 m,寬約640 m,高程在279.5~285.9 m之間。高大的沙卵石江心洲將嘉陵江水一分為二,形成左右兩汊。左汊較順直,但枯期水深不足,右汊為目前航槽,但該槽彎曲半徑較小,僅約380 m,不能滿足高等級(jí)航道建設(shè)標(biāo)準(zhǔn),航行條件較差。

    2 二維水沙數(shù)學(xué)模型及驗(yàn)證

    2.1 模型控制方程

    2.1.1 二維水流運(yùn)動(dòng)方程

    1)水流連續(xù)方程

    (1)

    2)ξ方向動(dòng)量方程

    (2)

    3)η方向動(dòng)量方程

    (3)

    式中:ξ,η分別為正交曲線坐標(biāo)系中的兩個(gè)正交曲線坐標(biāo);u,v分別為沿ξ,η方向的流速;h為水深;H為水位;g為重力加速度;Cξ,Cη為正交曲線坐標(biāo)系中的拉梅系數(shù)。

    2.1.2 泥沙運(yùn)動(dòng)方程

    1) 二維懸移質(zhì)不平衡輸沙基本方程

    SL)

    (4)

    2) 河床變形方程

    (5)

    2.2 模型計(jì)算范圍及網(wǎng)格剖分

    模擬中數(shù)學(xué)模型的入口邊界選定在馬回樞紐出水口上游約2.4 km的長(zhǎng)生溝處,出口邊界選定在正源鎮(zhèn)下碼頭,計(jì)算區(qū)域總長(zhǎng)約8.5 km。

    圖2 模型網(wǎng)格Fig.2 Grid graph of model

    模型網(wǎng)格主要由六節(jié)點(diǎn)的三角形網(wǎng)格組成(圖2),共有56 805個(gè)網(wǎng)格節(jié)點(diǎn),27 984個(gè)網(wǎng)格單元。網(wǎng)格單元尺寸順?biāo)鞣较蚣s為30~60 m,沿河寬方向間距約為30~50 m,并在馬回樞紐下引航道附近區(qū)域進(jìn)行局部加密,尺寸約為5~10 m。采用上游流量、下游水位作為模型的邊界條件。模型參數(shù)中紊動(dòng)擴(kuò)散系數(shù)采用25 m2/s;糙率采用0.025~0.035。

    2.3 模型驗(yàn)證

    2.3.1 水位驗(yàn)證

    模型驗(yàn)證所采用的資料為2014年3月實(shí)測(cè)的流速流向資料(Q=291 m3/s)及水位資料。圖3對(duì)平面二維數(shù)模的計(jì)算水位與原型實(shí)測(cè)水位進(jìn)行了比較,二者符合程度較高。從圖3的水位偏差值可見,一般偏差在±0.08 m以內(nèi),個(gè)別水尺最大相差也未超過±0.1 m,可見數(shù)模計(jì)算與原型實(shí)測(cè)資料吻合較好。

    圖3 水位驗(yàn)證Fig.3 Water level verification

    2.3.2 流速驗(yàn)證

    應(yīng)用建立的水流數(shù)學(xué)模型,對(duì)擬定的驗(yàn)證流量進(jìn)行了流場(chǎng)數(shù)值模擬。并對(duì)數(shù)模計(jì)算流速與原型實(shí)測(cè)資料進(jìn)行了對(duì)比。從圖4可見,兩者流速大小和分布以及最大值、最小值的位置均較為一致。各測(cè)點(diǎn)流速的計(jì)算值與實(shí)測(cè)值之間的差值大多在±0.01 m/s內(nèi),極個(gè)別誤差較大的也不超過±0.02 m/s,相對(duì)偏差均在±10%以內(nèi)。

    圖4 Q=291 m3/s實(shí)測(cè)流速驗(yàn)證Fig.4 Measured flow velocity in verification Q=291 m3/s

    2.3.3 河床沖淤驗(yàn)證

    本次計(jì)算以2010年4月實(shí)測(cè)地形建立模型,按2010年4月—2014年3月之間的水沙過程進(jìn)行計(jì)算,對(duì)照2014年3月地形檢驗(yàn)河床沖淤的相似性。經(jīng)過對(duì)沖淤時(shí)間比尺、輸沙率比尺的反復(fù)調(diào)整,最終達(dá)到與原型河床沖淤相似的要求。實(shí)測(cè)以及計(jì)算驗(yàn)證時(shí)段內(nèi)嘉陵江河段(8.5 km)的沖刷量和淤積量見表1。

    表1 沖淤量驗(yàn)證Table 1 The verified scale about scour and silting amount

    由表1可知,實(shí)測(cè)驗(yàn)證時(shí)段內(nèi)嘉陵江河段(8.5 km)的沖刷量為17.13×104m3,淤積量為14.53×104m3;計(jì)算相同河段同一時(shí)期內(nèi)的的沖刷量為11.95×104m3,淤積量為13.47×104m3。實(shí)測(cè)與計(jì)算的淤積量相差不大,但沖刷量卻存在一定差別,究其原因主要是在該河段存在人為采砂的影響。

    3 整治原則及方案

    3.1 整治原則

    根據(jù)八字腦復(fù)雜灘群河段的水沙特性、淺灘成因、演變規(guī)律及整治要求,為適應(yīng)嘉陵江航道建設(shè)大幅度提高航道等級(jí)的需要,擬定其整治思路為在該河段歷次整治的基礎(chǔ)上,采用疏浚、筑壩相結(jié)合措施,優(yōu)化整治線型,完善丁壩布置,加固穩(wěn)定邊灘,縮窄整治線寬度和提高整治水位,以增強(qiáng)輸沙能力,減少泥沙回淤,維持航槽穩(wěn)定。

    1)八字腦灘群航道整治宜采用疏浚與筑壩相結(jié)合措施,并根據(jù)各淺灘成因及演變規(guī)律各有側(cè)重。

    2)在該河段歷次整治的基礎(chǔ)上,完善整治建筑物布置,加大整治力度,延長(zhǎng)沖刷時(shí)間,增加航槽流速,減少泥沙回淤。

    3)依托主導(dǎo)河岸,利用河道深槽,順應(yīng)中、枯水河勢(shì),優(yōu)化整治線布置,構(gòu)成平順的整治線型。

    4)挖槽軸線沿河道深泓布置,并盡量與中、枯水水流動(dòng)力軸線保持一致,以利航槽的穩(wěn)定。

    3.2 整治標(biāo)準(zhǔn)

    八字腦灘群航道整治模型范圍上起馬回電站出水口上游2.4 km,下至正源鎮(zhèn)下碼頭,全長(zhǎng)約8.5 km,河段擬按內(nèi)河Ⅳ級(jí)航道進(jìn)行整治,航道尺度為為1.9 m×50 m×480 m(水深×槽寬×彎曲半徑)。

    測(cè)量期間在該灘設(shè)置了6把基本水尺,從上游至下游6把基本水尺的設(shè)計(jì)水位均為279.50 m,相應(yīng)設(shè)計(jì)流量為268 m3/s。

    根據(jù)理論、計(jì)算、結(jié)合嘉陵江優(yōu)良河段分析和整治實(shí)踐經(jīng)驗(yàn)等,綜合確定該灘的整治水位為設(shè)計(jì)水位上1.2 m,整治線寬度為160 m,相應(yīng)整治流量為1 100 m3/s。

    3.3 整治方案

    3.3.1 方案一

    方案一整治線走向和位置依托左側(cè)河岸,利用比較堅(jiān)實(shí)的河岸、硬角、磯頭作為整治線的控制點(diǎn),沿河道深泓布置成微彎平順的連續(xù)曲線,兩組反向曲線間以直線聯(lián)接。進(jìn)口段整治線上接馬回樞紐船閘下引航道,與主流方向斜向相交,逐漸向左岸過渡,大石塘灘靠左布置;出口段與左岸深槽平順相接。具體布置及整治措施如圖5。

    在進(jìn)口段補(bǔ)建馬回船閘下引航道調(diào)順段,并在引航道口門末端布置1#炸礁,同時(shí)在左岸筑1#堵壩一座,1#護(hù)岸一座,以封堵支汊,集中水流,穩(wěn)定航槽;在大石塘灘右汊筑1#鎖壩一座及2#、3#護(hù)岸兩座,以調(diào)節(jié)兩汊分流比,保護(hù)邊灘,維持新開航槽的穩(wěn)定;在出口段筑2#洲尾順壩,以束水歸槽,維護(hù)航槽穩(wěn)定;在水深不足1.9 m的淺區(qū)布置1#、2#、3#挖槽3處,以滿足航道尺度要求。

    圖5 方案一布置Fig.5 Arrangement plan of scenario 1

    3.3.2 方案二

    方案二八字腦灘整治線前段布置同方案一,大石塘及滿天星灘將整治線布置于原右汊主航槽,并沿河道深泓布置成微彎平順的連續(xù)曲線。具體布置及整治措施如圖6。

    在進(jìn)口段補(bǔ)建馬回船閘下引航道調(diào)順段,并在引航道口門末端布置1#炸礁,同時(shí)在左岸筑1#堵壩一座,1#護(hù)岸一座,以封堵支汊,集中水流,穩(wěn)定航槽;在大石塘灘上游筑1#順壩一座,右汊進(jìn)口筑2#、3#護(hù)岸兩座,以增加主汊流量,調(diào)順?biāo)?,保護(hù)邊灘;在滿天星灘段右側(cè)筑順壩一座,并切除左側(cè)凸岸1#邊灘,同時(shí)拆除原有丁壩(即1#拆壩),以增大航槽彎曲半徑和拓寬;在水深不足1.9 m的淺區(qū)布置1#、2#挖槽兩處,以滿足航道尺度要求。

    圖6 方案二布置Fig.6 Arrangement plan of scenario 2

    4 數(shù)模計(jì)算結(jié)果分析與比較

    4.1 通航水流條件

    經(jīng)數(shù)模計(jì)算兩個(gè)方案實(shí)施后,水流均平緩,流態(tài)均良好,航槽流速、比降也相差不大,且都在船舶航行的允許范圍內(nèi),具體見表2。

    在Q=268 m3/s時(shí),兩個(gè)方案較天然情況下的流速都有增加,但變化不大,水面比降在0.07‰~0.38‰之間。Q=1 100 m3/s時(shí),方案二的航槽流速增幅較方案一更大,特別是CS21斷面較天然情況增大89%,航槽內(nèi)橫向流速過大,彎曲半徑也相對(duì)偏小,可能導(dǎo)致船舶上行難度增加。Q=4 000 m3/s時(shí),方案一和方案二的航槽流速和水面比降與天然情況相比都有所降低,能較好的改善船舶航行的流態(tài)環(huán)境。

    表2 方案流速、比降對(duì)比Table 2 The comparison of velocity and gradient ratio of 2 scenarios

    4.2 泥沙沖淤情況

    嘉陵江八字腦灘群河段無(wú)推移質(zhì)泥沙實(shí)測(cè)資料,作為現(xiàn)場(chǎng)查勘取樣參考,共取7個(gè)基坑試樣。經(jīng)分析,八字腦灘群7個(gè)試坑平均中值粒徑為37.7 mm,最大粒徑為220 mm。其水文資料根據(jù)金溪水文站2000—2005年日流量資料,得到6年平均流量為549 m3/s;其中2003年平均流量為595 m3/s,與6年平均流量接近,故以此作為本次計(jì)算的代表年,并對(duì)其進(jìn)行水沙過程概化和計(jì)算,其計(jì)算成果如下:

    方案一整治后挖槽普遍淤積,淤積總量約0.38萬(wàn)m3,淤積厚度在0.04~0.38 m之間。其中八字腦灘挖槽上段左側(cè)沖刷,沖刷深度在0.11~0.34 m;挖槽下段右側(cè)淤積,淤積厚度在0.07~0.31 m;大石塘灘左槽疏浚區(qū)以淤積為主,其淤積厚度在0.04~0.38 m之間。各挖槽淤積情況見表3。

    表3 方案一挖槽淤積情況Table 3 The scale about channeling sedimentation of scenario 1

    方案二實(shí)施后挖槽內(nèi)泥沙淤積量為0.43萬(wàn)m3,挖槽沖淤變化基本同方案一。其中大石塘灘受彎道水流影響,右側(cè)河床即凸岸淤積較多,淤積厚度在0.04~0.32 m;滿天星灘受彎道水流及水庫(kù)壅水的雙重影響,挖槽普遍淤積,淤積厚度在0.15~0.37 m。各挖槽淤積情況見表4。

    表4 方案二挖槽淤積情況Table 4 The scale about channeling sedimentation of scenario 2

    因此,綜合表3和表4比較分析得出兩個(gè)結(jié)論:①兩個(gè)方案最小航槽寬度均為41 m,可滿足Ⅳ級(jí)航道35 m的最小航寬要求;②一個(gè)水文年后,兩個(gè)方案在疏浚區(qū)均有一定的泥沙回淤,但方案二的挖槽內(nèi)泥沙淤積量較方案一略為更大些,其中在八字腦灘兩個(gè)方案的回淤量及淤積分布相近,但大石塘灘左、右兩槽淤積量相差約20%。

    4.3 航道穩(wěn)定性分析

    對(duì)于多個(gè)淺段組成的淺灘群挖槽的穩(wěn)定來(lái)說,要求輸沙能力保持沿程遞增的態(tài)勢(shì),即下游淺灘的輸沙能力應(yīng)大于上游淺灘;淺灘的輸沙能力應(yīng)大于深槽。只有這樣,航道挖槽的穩(wěn)定才能得以維持。

    筆者采用沙莫夫推移質(zhì)輸沙率公式,并結(jié)合試驗(yàn)資料計(jì)算整治流量下各段的輸沙率:

    (6)

    式中:UC為起動(dòng)流速;中值粒徑D以mm計(jì)。

    計(jì)算整治后各淺區(qū)、深槽的推移質(zhì)輸沙率,分析整治工程實(shí)施后的航道穩(wěn)定。輸沙率結(jié)果見表5和表6。

    表5 方案一工程河段推移質(zhì)輸沙率(Q=1 100 m3/s)Table 5 The scale about bedload transport rate of scenario 1

    表6 方案二工程河段推移質(zhì)輸沙率(Q=1 100 m3/s)Table 6 The scale about bedload transport rate of scenario 2

    由上可知,從整治效果看(航槽穩(wěn)定),兩方案整治后該灘下段的輸沙能力略大于上段,淺區(qū)的輸沙能力明顯大于深槽,均有利于挖槽的穩(wěn)定。但在滿天星灘的局部航段,方案二受彎道水流和水庫(kù)壅水的雙重影響,造成局部挖槽淤積嚴(yán)重,可能影響航道穩(wěn)定。

    5 結(jié) 論

    1)建立的數(shù)學(xué)模型及驗(yàn)證計(jì)算的水面線、流速分布及河床沖淤等與原型實(shí)測(cè)資料基本吻合,可用于進(jìn)行航道整治模擬計(jì)算。

    2)水流條件計(jì)算表明:兩個(gè)方案整治后船閘下引航道與八字腦右槽平順銜接,水流平緩,流態(tài)良好,通航水流條件明顯改善。但右槽(方案二)在整治流量下航槽局部段橫向流速增大,可能導(dǎo)致船舶上行難度增加。

    3)河床沖淤計(jì)算表明:一個(gè)水文年后,兩個(gè)方案疏浚區(qū)均有一定的泥沙回淤,但淤積后的航道寬度均能滿足Ⅳ級(jí)航道的最小航寬。但方案二在大石塘灘航段受彎道水流及水庫(kù)壅水的雙重影響,致使泥沙回淤顯著,可能會(huì)破壞挖槽的穩(wěn)定性。

    4)從整治效果看(通航水流條件、航槽穩(wěn)定等),方案一相對(duì)較優(yōu),泥沙回淤量較少,航槽相對(duì)穩(wěn)定,可作為該灘整治方案。

    [1] 鞠俊.河口航道整治工程水動(dòng)力及河床變形響應(yīng)數(shù)值模擬研究[D].南京:南京水利科學(xué)研究院,2008. JU Jun.ComputationalRsearchontheHydrodynamicandBedMorphologicResponseofWaterwayRegulationEngineeringatEstuarine[D].Nanjing:Nanjing Hydraulic Research Institute,2008.

    [2] 曹民雄,蔡國(guó)正,王秀紅.長(zhǎng)江上游和尚巖灘群航道整治研究[J].水利水運(yùn)工程學(xué)報(bào),2010(3):1-9. CAO Minxiong,CAI Guozheng,WANG Xiuhong.Waterway regulation of Monk Rock rapids group in the upper reaches of the Yangtze River[J].Hydro-ScienceandEngneering,2010(3):1-9.

    [3] 楊禮生,呂升奇,劉貴平.贛江泉港灘群Ⅲ級(jí)航道整治工程方案對(duì)比及整治效果分析[J].水運(yùn)工程,2012(8):102-106. YANG Lisheng,LV Shengqi,LIU Guiping.Scheme comparison and effect analysis on waterway engineering of Qunguang reach downstream the Ganjiang river[J].Port&WaterwayEngineering,2012(8):102-106.

    [4] 費(fèi)曉昕,張緒進(jìn),倪志輝,等.北江白石窯水利樞紐壩下變動(dòng)回水區(qū)航道整治試驗(yàn)研究[J].水運(yùn)工程,2013(4):146-150. FEI Xiaoxin,ZHANG Xujin,NI Zhihui,et al.Model experiment on fluctuating backwater area under the dam of Beijiang Baishiyao water conservancy[J].Port&WaterwayEngineering,2013(4):146-150.

    [5] RAJARANTNAM N,NWATCHUKWU B A.Flow near groin-like structures[J].JournalofHydraulicEngineering,1983,109(3):463- 480.

    [6] SURIAN N.Downstream variation in grain size along an Alpine river:analysis of controls and processes[J].Geomorphology,2002,43(1/2):137-149.

    [7] 倪志輝,易靜,張緒進(jìn).感潮河網(wǎng)區(qū)航道整治工程潮流數(shù)值模擬[J].水運(yùn)工程,2013(6):108-113. NI Zhihui,YI Jing,ZHANG Xujing.Tidal flow numerical simulation study on regulation project in tidal river network area[J].Port&WaterwayEngineering,2013(6):108-113.

    [8] 王秀紅,曹民雄,馬愛興,等.烏江沙陀電站變動(dòng)回水區(qū)航道整治二維水流數(shù)學(xué)模型研究[J].水運(yùn)工程,2012(11):156-160. WANG Xiuhong,CAO Minxiong,MA Aixing,et al.Fluctuating backwater area waterway regulation of Wujiang Shatuo hydropower station by 2D mathematical model[J].Port&WaterwayEngineering,2012(11):156-160.

    [9] 陳界仁,陳國(guó)祥.贛江航道整治工程泥沙數(shù)學(xué)模型[J].河海大學(xué)學(xué)報(bào)(自然科學(xué)版),2003,31(2):224-228. CHEN Jieren,CHEN Guoxiang.2-D sediment mathematical model for Ganjiang river waterway project[J].JournalofHohaiUniversity(NaturalSciences),2003,31(2):224-228.

    [10] 王義安,陸永軍.二維泥沙數(shù)學(xué)模型在松花江五股流航道整治中的應(yīng)用[J].水利水運(yùn)工程學(xué)報(bào),2002(2):61-65. WANG Yi’an,LU Yongjun.Application of 2-D sediment mathematical model in training Wuguliu reach channel in Songhuajiang river[J].Hydro-ScienceandEngineering,2002(2):61-65.

    [11] 倪志輝,劉亞博,張緒進(jìn),等.北江白石窯樞紐壩下航道穩(wěn)定性初步分析[J].水運(yùn)工程,2014(6):90-95. NI Zhihui,LIU Yabo,ZHANG Xujing.Analysis of waterway stability under the dam of Beijiang Baishiyao water conservancy[J].Port&WaterwayEngineering,2014(6):90-95.

    [12] 倪志輝,劉亞博,吳立春,等.岷江老木孔航電樞紐庫(kù)區(qū)泥沙淤積計(jì)算與分析[J].長(zhǎng)江流域資源與環(huán)境.2015,24 (1):149-155. NI Zhihui,LIU Yabo,WU Lichun.Calculation and analysis of the sediment deposition in Laomukong navigation-power junction on the Minjiang river[J].ResourcesandEnvironmentintheYangtzeBasin,2015,24 (1):149-155.

    Numerical Model of Channel Regulation for Bazinao Shoal Group

    NI Zhihui1 2,GUO Yi1,WU Lichun3

    (1. Key Laboratory of Hydraulic and Waterway Engineering of the Ministry of Education,National Engineering Research Center for Inland Waterway Regulation,Chongqing Jiaotong University,Chongqing 400074, P. R. China;2. Southwestern Research Institute of Water Transportation Engineering,Chongqing Jiaotong University, Chongqing 400016, P. R. China; 3.Chongqing University of Education, Chongqing 400067, P. R. China)

    Jialing River Bazinao group shoal section is located in the tail of Fengyi navigation power junction reservoir. Due to the separation of the water discharge outlet of Mahui power station and vessel lock exit, the right branch of Bazinaotan is of narrow channel of shallow water and the left branch is extremely curved and narrow, which can not accommodate navigation. Further in Dashi shoal section the sand-pebble bar in water splits Jialing river into two branches right and left of which the left branch is relatively straight but is of insufficient water depth in dry season while the right branch is of smaller curve radius and does not meet requirement for high class channel construction. Aiming at these navigation impeding characteristic problems of Bazinaotan shoal group reach of Jialing river, a mathematical model was set up to calculate, analyze and compare the two different treatment schemes for current and sediment conditions. The results show that this mathematical model can well simulate reach flow and sediment accumulation situation. As to the end result of regulation and rectification, scheme one is superior for such advantages of less siltation, more stable channel and thus can be taken as implementation scheme for reference of design.

    waterway engineering; mathematical model;sediment deposition; shoal group; navigation condition

    10.3969/j.issn.1674-0696.2016.02.16

    2015-02-02;

    2015-06-01

    國(guó)家自然科學(xué)基金項(xiàng)目(41306078);內(nèi)河航道整治技術(shù)交通行業(yè)重點(diǎn)實(shí)驗(yàn)室項(xiàng)目(NHHD-201514);國(guó)家內(nèi)河航道整治工程技術(shù)研究中心水利水運(yùn)工程教育部重點(diǎn)實(shí)驗(yàn)室項(xiàng)目(SLK2016B03)

    倪志輝(1980—),男,湖南衡陽(yáng)人,副研究員,博士,主要從事河流海岸水動(dòng)力學(xué)、環(huán)境數(shù)值模擬方面的研究。E-mail:benny251@163.com。

    O319.56

    A

    1674-0696(2016)02-070-05

    猜你喜歡
    八字河段淤積
    短短八字,道盡一生
    意林(2023年18期)2023-11-05 10:45:10
    淤積性皮炎知多少
    妊娠期肝內(nèi)膽汁淤積癥
    肝博士(2020年5期)2021-01-18 02:50:28
    SL流量計(jì)在特殊河段的應(yīng)用——以河源水文站為例
    石瀘高速公路(紅河段)正式通車
    石油瀝青(2020年1期)2020-05-25 06:54:10
    5類人跟泡腳“八字不合”
    淤積與浚疏:清朝時(shí)期福州西湖的治理史
    膽汁淤積性肝病問題解答
    肝博士(2015年2期)2015-02-27 10:49:51
    淺析黃河河段維修養(yǎng)護(hù)施工工藝
    河南科技(2014年14期)2014-02-27 14:11:44
    防凌期寧蒙河段流量演進(jìn)的分期分河段混合算法研究
    国语对白做爰xxxⅹ性视频网站| 在线看a的网站| 久久久久久久久久久久大奶| 老鸭窝网址在线观看| 亚洲欧洲国产日韩| 国产精品久久久久久av不卡| 黄色 视频免费看| 1024香蕉在线观看| 中文字幕人妻丝袜制服| 爱豆传媒免费全集在线观看| 久久国内精品自在自线图片| 精品亚洲成国产av| 亚洲第一区二区三区不卡| 亚洲av成人精品一二三区| 国产日韩一区二区三区精品不卡| 国产成人aa在线观看| 午夜福利视频在线观看免费| 狠狠婷婷综合久久久久久88av| 久久这里只有精品19| 看免费av毛片| 一二三四在线观看免费中文在| 男人添女人高潮全过程视频| 肉色欧美久久久久久久蜜桃| av免费在线看不卡| 欧美 亚洲 国产 日韩一| 伊人久久国产一区二区| 国产精品av久久久久免费| 少妇人妻 视频| 免费观看无遮挡的男女| 亚洲av在线观看美女高潮| 九草在线视频观看| 亚洲天堂av无毛| 美女xxoo啪啪120秒动态图| 一二三四在线观看免费中文在| 久久婷婷青草| 最新的欧美精品一区二区| 欧美中文综合在线视频| 亚洲一区中文字幕在线| 国产色婷婷99| 天天躁夜夜躁狠狠久久av| 一边亲一边摸免费视频| 亚洲欧美日韩另类电影网站| 午夜福利一区二区在线看| 国产成人精品无人区| 国产一区亚洲一区在线观看| 亚洲欧美中文字幕日韩二区| 亚洲精品久久成人aⅴ小说| 满18在线观看网站| 国产精品欧美亚洲77777| 下体分泌物呈黄色| 成年动漫av网址| 国产成人精品无人区| 高清av免费在线| 免费观看无遮挡的男女| 亚洲五月色婷婷综合| 老汉色∧v一级毛片| 自线自在国产av| 免费在线观看完整版高清| 欧美最新免费一区二区三区| av免费观看日本| 精品第一国产精品| 欧美日韩亚洲高清精品| 夫妻性生交免费视频一级片| 日韩中文字幕视频在线看片| 一二三四在线观看免费中文在| 国产av一区二区精品久久| 交换朋友夫妻互换小说| 欧美国产精品一级二级三级| 久久精品夜色国产| 久久精品国产亚洲av涩爱| 国产 一区精品| 黄片播放在线免费| 日韩制服丝袜自拍偷拍| 亚洲欧美精品自产自拍| 伦理电影免费视频| 国产一级毛片在线| 黑丝袜美女国产一区| 国产97色在线日韩免费| 欧美日韩精品网址| 午夜福利,免费看| 国产 一区精品| 啦啦啦在线免费观看视频4| 国产欧美日韩一区二区三区在线| 色网站视频免费| 最近最新中文字幕免费大全7| 国产在视频线精品| 久久这里有精品视频免费| 日韩av不卡免费在线播放| av一本久久久久| 少妇人妻精品综合一区二区| 日韩,欧美,国产一区二区三区| 少妇被粗大猛烈的视频| 我要看黄色一级片免费的| 国产成人精品一,二区| 亚洲第一区二区三区不卡| 国产男人的电影天堂91| 老司机影院成人| 精品一品国产午夜福利视频| 国产精品无大码| 日本av手机在线免费观看| 日韩中文字幕欧美一区二区 | 亚洲av电影在线进入| 大话2 男鬼变身卡| 国产精品国产三级专区第一集| 日韩制服骚丝袜av| 中国国产av一级| 高清欧美精品videossex| 国产成人欧美| 另类亚洲欧美激情| 久久女婷五月综合色啪小说| 满18在线观看网站| 日韩精品有码人妻一区| 精品第一国产精品| 不卡av一区二区三区| 香蕉精品网在线| 午夜福利乱码中文字幕| 欧美日韩一区二区视频在线观看视频在线| 欧美日韩一级在线毛片| 99久国产av精品国产电影| 国产日韩欧美在线精品| 亚洲人成电影观看| 精品少妇内射三级| 日韩精品免费视频一区二区三区| 如日韩欧美国产精品一区二区三区| 高清不卡的av网站| 欧美97在线视频| 91午夜精品亚洲一区二区三区| 制服人妻中文乱码| 精品亚洲乱码少妇综合久久| 久久免费观看电影| 美女福利国产在线| 午夜免费观看性视频| 只有这里有精品99| 咕卡用的链子| 国产又色又爽无遮挡免| 久久这里只有精品19| www日本在线高清视频| 国产av一区二区精品久久| 中文欧美无线码| 国产成人免费无遮挡视频| 十八禁高潮呻吟视频| 久久精品国产亚洲av高清一级| 国产精品欧美亚洲77777| 丰满饥渴人妻一区二区三| 亚洲精品美女久久av网站| 性少妇av在线| 国产精品国产三级国产专区5o| 黄色配什么色好看| 国产精品免费大片| 熟妇人妻不卡中文字幕| 亚洲精品视频女| 久久久久久久精品精品| 亚洲少妇的诱惑av| 人成视频在线观看免费观看| 一二三四中文在线观看免费高清| 午夜免费观看性视频| 国产精品一国产av| 亚洲精品日本国产第一区| 国产毛片在线视频| 亚洲精品乱久久久久久| 久久久久久免费高清国产稀缺| 成人免费观看视频高清| 可以免费在线观看a视频的电影网站 | 在现免费观看毛片| 黄色配什么色好看| 久久久精品94久久精品| 久久精品国产a三级三级三级| 亚洲av日韩在线播放| 国产人伦9x9x在线观看 | 亚洲国产欧美在线一区| 欧美日韩亚洲高清精品| 天天躁夜夜躁狠狠久久av| 啦啦啦中文免费视频观看日本| 天天躁狠狠躁夜夜躁狠狠躁| 国产熟女午夜一区二区三区| 曰老女人黄片| 久久女婷五月综合色啪小说| 精品第一国产精品| 精品国产一区二区三区四区第35| 国产欧美日韩综合在线一区二区| 一本大道久久a久久精品| 午夜免费男女啪啪视频观看| 亚洲一级一片aⅴ在线观看| 免费黄色在线免费观看| 丝袜美腿诱惑在线| 久久青草综合色| 91成人精品电影| 2022亚洲国产成人精品| 中文字幕精品免费在线观看视频| 日韩,欧美,国产一区二区三区| 可以免费在线观看a视频的电影网站 | 乱人伦中国视频| 天堂中文最新版在线下载| 亚洲情色 制服丝袜| 亚洲精品一二三| 精品福利永久在线观看| 日本av手机在线免费观看| 欧美精品高潮呻吟av久久| 久久99蜜桃精品久久| 性色avwww在线观看| 日韩视频在线欧美| 精品久久久久久电影网| 妹子高潮喷水视频| av网站免费在线观看视频| av在线播放精品| 美女中出高潮动态图| 波多野结衣一区麻豆| 亚洲四区av| 熟女av电影| av免费在线看不卡| 999久久久国产精品视频| 边亲边吃奶的免费视频| 亚洲综合精品二区| 精品酒店卫生间| 有码 亚洲区| 中文字幕亚洲精品专区| 国产黄色免费在线视频| www日本在线高清视频| 亚洲国产日韩一区二区| 大片电影免费在线观看免费| 国产精品蜜桃在线观看| 国产极品天堂在线| 曰老女人黄片| 肉色欧美久久久久久久蜜桃| 大片免费播放器 马上看| 国产不卡av网站在线观看| 久久这里只有精品19| 一级黄片播放器| 国产av码专区亚洲av| 久久人人爽av亚洲精品天堂| 国产av一区二区精品久久| 亚洲伊人久久精品综合| 18禁动态无遮挡网站| 一级片免费观看大全| 嫩草影院入口| 九草在线视频观看| 青青草视频在线视频观看| 国产xxxxx性猛交| 欧美日韩一级在线毛片| 精品国产超薄肉色丝袜足j| 天堂中文最新版在线下载| 亚洲,一卡二卡三卡| 精品一品国产午夜福利视频| 成年女人毛片免费观看观看9 | 天美传媒精品一区二区| 伊人久久大香线蕉亚洲五| 不卡视频在线观看欧美| 日韩人妻精品一区2区三区| 欧美日韩亚洲国产一区二区在线观看 | 日日撸夜夜添| 亚洲国产成人一精品久久久| 91aial.com中文字幕在线观看| 亚洲精品日本国产第一区| 伦理电影大哥的女人| 精品一区在线观看国产| 男人添女人高潮全过程视频| 一本—道久久a久久精品蜜桃钙片| 另类精品久久| 这个男人来自地球电影免费观看 | 日本91视频免费播放| 男人操女人黄网站| 一区二区三区激情视频| 午夜日本视频在线| 五月天丁香电影| 国产成人精品在线电影| 啦啦啦视频在线资源免费观看| 热99久久久久精品小说推荐| 国产精品免费大片| 精品国产一区二区三区四区第35| 各种免费的搞黄视频| 成年av动漫网址| 国产精品成人在线| 亚洲国产色片| 寂寞人妻少妇视频99o| 日日爽夜夜爽网站| 亚洲精品美女久久av网站| 亚洲精品国产色婷婷电影| 久久久国产欧美日韩av| 极品人妻少妇av视频| 亚洲av综合色区一区| 天天操日日干夜夜撸| 一区二区三区乱码不卡18| 久久久久精品人妻al黑| 日韩成人av中文字幕在线观看| 在线观看免费视频网站a站| 美女国产视频在线观看| 久久精品国产a三级三级三级| 久久久亚洲精品成人影院| 在线观看国产h片| 一二三四中文在线观看免费高清| 久久久国产一区二区| 日韩一区二区三区影片| 99国产精品免费福利视频| 日韩不卡一区二区三区视频在线| 美女福利国产在线| 国产色婷婷99| 久久久亚洲精品成人影院| 各种免费的搞黄视频| 日韩 亚洲 欧美在线| 欧美激情极品国产一区二区三区| 日日撸夜夜添| a 毛片基地| 精品一品国产午夜福利视频| 亚洲美女黄色视频免费看| av国产久精品久网站免费入址| 日韩视频在线欧美| 亚洲精品久久久久久婷婷小说| 99久久综合免费| 啦啦啦啦在线视频资源| 国产av精品麻豆| 一边亲一边摸免费视频| 三上悠亚av全集在线观看| h视频一区二区三区| 一区在线观看完整版| 一级毛片我不卡| 国产精品香港三级国产av潘金莲 | 丝袜脚勾引网站| 欧美人与善性xxx| 亚洲精品乱久久久久久| 美女国产高潮福利片在线看| 国产成人一区二区在线| 少妇被粗大猛烈的视频| 亚洲国产最新在线播放| 亚洲人成网站在线观看播放| freevideosex欧美| 高清欧美精品videossex| 久久久久久久亚洲中文字幕| 母亲3免费完整高清在线观看 | 麻豆乱淫一区二区| www.av在线官网国产| 日韩一卡2卡3卡4卡2021年| 激情视频va一区二区三区| 在线天堂最新版资源| 欧美另类一区| 丝袜在线中文字幕| 国产一区二区激情短视频 | 国产一区亚洲一区在线观看| 在现免费观看毛片| 国产精品蜜桃在线观看| 波野结衣二区三区在线| 午夜久久久在线观看| 亚洲图色成人| 2022亚洲国产成人精品| 亚洲av免费高清在线观看| 精品少妇黑人巨大在线播放| 性少妇av在线| 欧美日韩精品网址| 成人免费观看视频高清| www.熟女人妻精品国产| 久久精品国产亚洲av天美| 人成视频在线观看免费观看| 欧美国产精品一级二级三级| 天堂俺去俺来也www色官网| 欧美激情极品国产一区二区三区| 精品一区二区三区四区五区乱码 | 狠狠精品人妻久久久久久综合| av.在线天堂| 久久热在线av| 99久久人妻综合| 久久久久久久久久久久大奶| 欧美日韩视频高清一区二区三区二| 亚洲图色成人| 视频在线观看一区二区三区| 毛片一级片免费看久久久久| 欧美av亚洲av综合av国产av | 国产 一区精品| 成年人午夜在线观看视频| 人人妻人人爽人人添夜夜欢视频| 日韩大片免费观看网站| 久久久精品免费免费高清| 久久人人爽av亚洲精品天堂| 精品国产露脸久久av麻豆| 亚洲美女视频黄频| 91久久精品国产一区二区三区| 精品久久久久久电影网| 极品人妻少妇av视频| 日韩一卡2卡3卡4卡2021年| 美女国产视频在线观看| 久久热在线av| 国产高清不卡午夜福利| 国产一区有黄有色的免费视频| 99九九在线精品视频| 天天躁狠狠躁夜夜躁狠狠躁| 男女下面插进去视频免费观看| 777米奇影视久久| 亚洲欧美一区二区三区久久| 极品少妇高潮喷水抽搐| 午夜免费男女啪啪视频观看| 亚洲欧美精品自产自拍| 视频区图区小说| 高清不卡的av网站| 99久久人妻综合| 亚洲成人av在线免费| 少妇被粗大猛烈的视频| 中文字幕av电影在线播放| 看非洲黑人一级黄片| 一区二区av电影网| 国产男人的电影天堂91| 成人影院久久| 精品少妇久久久久久888优播| 九色亚洲精品在线播放| 五月天丁香电影| 免费日韩欧美在线观看| 欧美成人午夜精品| 国产精品蜜桃在线观看| 亚洲久久久国产精品| 这个男人来自地球电影免费观看 | 校园人妻丝袜中文字幕| 国产一区二区三区综合在线观看| 日本爱情动作片www.在线观看| 国产一区二区三区综合在线观看| 午夜日韩欧美国产| 精品国产乱码久久久久久男人| 中文字幕精品免费在线观看视频| 男女下面插进去视频免费观看| 黄色视频在线播放观看不卡| 亚洲激情五月婷婷啪啪| 久久久久国产精品人妻一区二区| 大香蕉久久成人网| 飞空精品影院首页| 国产国语露脸激情在线看| www.自偷自拍.com| 一本—道久久a久久精品蜜桃钙片| 一个人免费看片子| 91国产中文字幕| 青青草视频在线视频观看| 亚洲婷婷狠狠爱综合网| 各种免费的搞黄视频| 国产一区二区 视频在线| 色播在线永久视频| 丰满饥渴人妻一区二区三| 欧美黄色片欧美黄色片| 国产黄频视频在线观看| 青春草亚洲视频在线观看| 亚洲国产成人一精品久久久| 搡女人真爽免费视频火全软件| 日日啪夜夜爽| 精品午夜福利在线看| 亚洲精品久久成人aⅴ小说| 丰满乱子伦码专区| 深夜精品福利| 91精品伊人久久大香线蕉| 国产av码专区亚洲av| 亚洲av免费高清在线观看| 亚洲一码二码三码区别大吗| 国产 精品1| 只有这里有精品99| 自线自在国产av| 日本vs欧美在线观看视频| 久久久a久久爽久久v久久| 在线观看三级黄色| 亚洲视频免费观看视频| xxxhd国产人妻xxx| 久久久久久久久久人人人人人人| 男人爽女人下面视频在线观看| 两个人看的免费小视频| 日韩不卡一区二区三区视频在线| 狠狠精品人妻久久久久久综合| 只有这里有精品99| 日韩 亚洲 欧美在线| 免费高清在线观看视频在线观看| 999精品在线视频| 日韩伦理黄色片| 九九爱精品视频在线观看| 免费不卡的大黄色大毛片视频在线观看| av在线老鸭窝| 国产高清国产精品国产三级| 成年人免费黄色播放视频| 亚洲国产欧美网| 黄频高清免费视频| 少妇 在线观看| 香蕉国产在线看| 黑人巨大精品欧美一区二区蜜桃| av又黄又爽大尺度在线免费看| 成年动漫av网址| 搡老乐熟女国产| 不卡视频在线观看欧美| av国产久精品久网站免费入址| 午夜福利在线观看免费完整高清在| 国产成人aa在线观看| 男人操女人黄网站| 在线观看人妻少妇| 一边亲一边摸免费视频| 深夜精品福利| av电影中文网址| 欧美精品一区二区大全| 欧美亚洲日本最大视频资源| 精品久久久久久电影网| 久久精品aⅴ一区二区三区四区 | 国产精品久久久久久久久免| 国产成人免费观看mmmm| 一区二区日韩欧美中文字幕| 亚洲色图综合在线观看| 国产精品av久久久久免费| 欧美日韩av久久| 国产一区二区三区av在线| 女人高潮潮喷娇喘18禁视频| 男女边吃奶边做爰视频| 日本黄色日本黄色录像| av电影中文网址| 视频在线观看一区二区三区| 久久久久久久国产电影| 在线观看三级黄色| 人人妻人人爽人人添夜夜欢视频| 午夜日韩欧美国产| 交换朋友夫妻互换小说| 18禁动态无遮挡网站| 99热网站在线观看| 亚洲视频免费观看视频| 校园人妻丝袜中文字幕| 美女脱内裤让男人舔精品视频| 欧美日韩av久久| 制服诱惑二区| www.自偷自拍.com| xxx大片免费视频| 欧美+日韩+精品| 一级毛片我不卡| a级毛片在线看网站| 欧美成人精品欧美一级黄| av免费在线看不卡| 亚洲精华国产精华液的使用体验| 国产1区2区3区精品| 大话2 男鬼变身卡| 精品国产超薄肉色丝袜足j| 亚洲伊人久久精品综合| 欧美变态另类bdsm刘玥| 亚洲成人手机| 男女边摸边吃奶| 国产欧美日韩一区二区三区在线| 精品久久蜜臀av无| 亚洲成av片中文字幕在线观看 | 大香蕉久久成人网| av线在线观看网站| 97在线视频观看| 九草在线视频观看| a级毛片黄视频| 99香蕉大伊视频| 美女午夜性视频免费| 精品少妇久久久久久888优播| 亚洲第一区二区三区不卡| 热re99久久精品国产66热6| 久久久久久伊人网av| 另类亚洲欧美激情| 午夜日韩欧美国产| 18在线观看网站| av线在线观看网站| 成人影院久久| 国产av码专区亚洲av| 天堂俺去俺来也www色官网| 少妇人妻久久综合中文| 久久亚洲国产成人精品v| 久久午夜福利片| 国产极品天堂在线| 亚洲少妇的诱惑av| 青春草亚洲视频在线观看| 精品国产露脸久久av麻豆| 国产精品一二三区在线看| 日本vs欧美在线观看视频| 亚洲欧美成人精品一区二区| 一级毛片我不卡| 中文欧美无线码| 国产在线一区二区三区精| 国产高清国产精品国产三级| 久久国产精品大桥未久av| 麻豆精品久久久久久蜜桃| 日日撸夜夜添| 久热这里只有精品99| 黄色怎么调成土黄色| 黑丝袜美女国产一区| 国产日韩欧美视频二区| 亚洲av.av天堂| 熟女av电影| 不卡av一区二区三区| 日本av免费视频播放| 一级爰片在线观看| 国产成人免费观看mmmm| 久久久久久久国产电影| 国产精品成人在线| 下体分泌物呈黄色| 色94色欧美一区二区| 男女国产视频网站| 中文字幕亚洲精品专区| 美女大奶头黄色视频| 国产日韩一区二区三区精品不卡| 国产亚洲精品第一综合不卡| 视频区图区小说| 中文字幕亚洲精品专区| av免费观看日本| 天堂俺去俺来也www色官网| 最新的欧美精品一区二区| 在线 av 中文字幕| 亚洲精品,欧美精品| 亚洲国产毛片av蜜桃av| 中文字幕最新亚洲高清| 中文字幕人妻丝袜制服| 亚洲精品在线美女| 又粗又硬又长又爽又黄的视频| 国产白丝娇喘喷水9色精品| 国产精品免费视频内射| 七月丁香在线播放| 伦精品一区二区三区| 女性被躁到高潮视频| 国产欧美日韩综合在线一区二区| 美女午夜性视频免费| 大片免费播放器 马上看| 最新的欧美精品一区二区| 久久久a久久爽久久v久久| 久久热在线av| 国产av精品麻豆| 妹子高潮喷水视频| 大码成人一级视频| 成人亚洲精品一区在线观看| 九色亚洲精品在线播放| 一本—道久久a久久精品蜜桃钙片| 久久精品国产亚洲av涩爱|