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

    黃河下游丁壩群對河勢控制效果研究

    2022-07-26 00:51:10竹,舟*,瑩,
    大連理工大學(xué)學(xué)報 2022年4期
    關(guān)鍵詞:工程

    彭 秀 竹, 徐 向 舟*, 趙 瑩, 高 璐

    ( 1.大連理工大學(xué) 水利工程學(xué)院, 遼寧 大連 116024;2.山東省水利科學(xué)研究院, 山東 濟南 250013 )

    0 引 言

    黃河下游河道不斷遷徙變化,在一定條件下就會決溢泛濫,改走新道.2019年國家將黃河流域生態(tài)保護和高質(zhì)量發(fā)展放在戰(zhàn)略高度,明確要求實施河道和灘區(qū)綜合提升治理工程.隨著小浪底水庫調(diào)水調(diào)沙的運用,下游河道的水沙條件發(fā)生了變化,原有的工程布置已經(jīng)不能滿足當(dāng)前河道治理情況,尤其是部分靠溜工程的上首坍塌嚴重,致使主溜超出工程控制范圍,這對防汛搶險十分不利,需要加快開展游蕩型河段控導(dǎo)工程的續(xù)建.然而河道整治工程的設(shè)計按照原控導(dǎo)工程確定是否合理,目前還缺乏必要的技術(shù)論證.鑒于黃河下游本身水沙特性的復(fù)雜性以及相關(guān)研究的缺乏,有必要開展整治工程對黃河下游游蕩型河道影響的研究.

    丁壩的布設(shè)和結(jié)構(gòu)形式的選擇直接關(guān)系到河道的治理效果.已有研究通過模型試驗和數(shù)值計算等方法,分析了丁壩對河道中流體特性的影響,相關(guān)成果大多聚焦在速度分布[1]、自由水面形態(tài)和紊動特性[2-3]等方面,比如已有許多學(xué)者根據(jù)模型試驗和數(shù)值模擬分析了單丁壩挑角[4]、長度[5]和壩型[6]等對流場及沖刷的影響.以往對丁壩的研究多是針對單丁壩某一具體問題深入分析,很少有把丁壩群當(dāng)作整體研究對象來考慮.事實上,單丁壩對河道整治有局限性,現(xiàn)實大多采用丁壩群系統(tǒng)配合作用,在避免產(chǎn)生新水害的同時,還可以達到河道整治的目的.加之彎道處壩體作用下的水流、泥沙和河床的相互作用機理復(fù)雜,因此仍需對河灣處壩群周圍水流特征和河床變形進行更加深入的分析,從而驗證丁壩群在緩流落淤、保護岸坡等方面的效果,以尋求具有顯著防洪效益的護岸方式.此外,以往對丁壩周圍流動特性的研究多集中于順直河道[5],而黃河下游為游蕩型河道,河灣段水沙作用機理復(fù)雜,在討論丁壩群效用機制時難以直接應(yīng)用現(xiàn)有成果.鑒于此,本文運用MIKE21軟件,構(gòu)建黃河下游渾水條件下典型游蕩型河道的二維水動力數(shù)值模型,對比分析在不同挑角、布設(shè)排數(shù)和結(jié)構(gòu)形式情況下,丁壩群對游蕩型河道平面河勢的控制效果,為丁壩設(shè)計施工提供參考,以期為黃河下游工程布局提供科學(xué)思路.

    1 區(qū)域概況

    全長約275 km的孟津至高村河段,位處黃河下游.該河段是典型的游蕩型河段,河勢擺動劇烈,易沖易淤,灘槽泥沙交換過程復(fù)雜.而老君堂河段整個位于游蕩型河段中,在山東省東明縣黃河南灘內(nèi).1974年開始在該河段修建老君堂控導(dǎo)工程,上接左岸周營口工程來溜,下送溜至榆林工程.從老君堂河勢變化發(fā)展來看,如果工程對灘岸灘尖在大溜作用下,淘刷嚴重,繼續(xù)坍塌,灘嘴將掉落,主流有進一步下滑的趨勢,屆時工程下首丁壩將會進一步出險和右岸灘岸坍塌加?。疄榉乐购觿堇^續(xù)下滑,使老君堂工程發(fā)揮應(yīng)有的迎溜送溜作用,維持目前河勢,計劃修建老君堂下延工程(圖1).

    圖1 研究區(qū)域Fig.1 Study area

    2 數(shù)值模擬

    2.1 模型計算范圍和數(shù)據(jù)來源

    本次數(shù)學(xué)模型的計算運用MIKE21軟件,建立基于非正交網(wǎng)格的老君堂附近典型天然河道二維水動力數(shù)值模型.計算范圍上起馬寨斷面,下至楊小寨,涵蓋整個老君堂河段,河段總長12 km,如圖1所示.模型應(yīng)用的是2017年黃河下游干流河道地形圖數(shù)據(jù),并在此基礎(chǔ)上添加了黃河水文勘測局對黃河下游老君堂河段實測的7個汛前斷面數(shù)據(jù).河道中原有控導(dǎo)工程地形采用不過流邊界形式處理.利用網(wǎng)格生成器對計算河段進行非結(jié)構(gòu)化、結(jié)構(gòu)化網(wǎng)格劃分,大部分采用三角網(wǎng)格,為了保持工程前后邊界的順直,在擬建工程的區(qū)域采用四邊形網(wǎng)格,并且進行網(wǎng)格加密處理.模擬河段共剖分了12 207個網(wǎng)格點.

    2.2 模擬試驗方案

    模型的邊界條件按照河道治理設(shè)計的水文部分,進口邊界條件選擇進口總流量,設(shè)計流量4 000 m3/s,設(shè)計含沙量30 kg/m3[7],出口邊界水面高程為63 m.本文采用控制變量法,以布設(shè)排數(shù)為變量設(shè)計了1、3、5排3種方案;以挑角為變量設(shè)計了30°、45°、75° 3種方案.以原有的整治工程設(shè)計方案作為制訂試驗方案的基本依據(jù),均為下挑實體丁壩,丁壩布設(shè)在河灣凹岸(圖1圈畫部分),尺寸按照原控導(dǎo)工程設(shè)置,壩長約100 m,壩間距約1倍壩長,壩頂高程為4 000 m3/s流量水位+0.5 m.在此基礎(chǔ)上,確定最佳布設(shè)方案,選取合理且相同的平面布設(shè)方式,對實體丁壩和透水丁壩進行數(shù)值模擬.其中透水丁壩透空率選擇為30%,結(jié)構(gòu)型式見圖2,其他設(shè)置與實體丁壩相同.同時為了定量分析實體丁壩與透水丁壩對凹岸的緩流效果,引入流速減小貢獻率這一概念,定義如下:

    圖2 透水丁壩立面示意圖Fig.2 Elevation schematic of permeable spur-dike

    K=(|vq-vz|/vz)×100%

    (1)

    式中:vq為無壩時測點流速,vz為設(shè)置丁壩后的測點流速.

    采用相對標(biāo)準偏差對丁壩壩區(qū)上下游斷面流速進行均勻性分析.相對標(biāo)準偏差計算公式如下:

    (2)

    (3)

    3 試驗結(jié)果

    3.1 實體丁壩不同挑角對流場的影響

    挑角與有效壩長直接相關(guān),其影響河道自由水面和流態(tài)穩(wěn)定.因此合理選擇挑角對矯正河勢和維護工程安全起著關(guān)鍵作用.本文模擬了單排實體丁壩不同挑角下河勢變化情況,選取實體丁壩壩軸線上下游20 m處典型斷面作為水位和流速的提取斷面,具體位置見圖3.實體丁壩周圍河勢變化情況如表1所示.從表1可以看出,丁壩上游斷面水位升高,隨著挑角的增大,水位雍高均值增大.當(dāng)丁壩挑角為75°時,水位雍高均值達到最大值1.388 m,同時流速均勻性指標(biāo)值Cv也達到3種角度中的最大值0.418.Cv越大,流速分布的均勻性越低.證明此挑角工況下,斷面流速變幅最大,水流動能和勢能之間轉(zhuǎn)化頻繁.在丁壩下游均發(fā)生跌水現(xiàn)象.下游Cv隨著挑角的增大逐漸減小,與上游變化正好相反.挑角為30°時,下游斷面流速變化相對活躍.由此可見,河灣處水沙機理復(fù)雜,在不同的斷面處,丁壩周圍的流速變幅隨挑角變化的規(guī)律不可一概而論.從水流本身特性來看,丁壩的存在致使原有河道過水面積突然減小,打破了水流的暢行,使水流流速重新分布.由于丁壩的束窄作用,流速在丁壩斷面處陡增.水流受到慣性作用影響,繞過壩頭后過流斷面進一步收縮,流速繼續(xù)增大.同時流速變化必然引起水流勢能與動能的轉(zhuǎn)化[8],水流過壩時勢能轉(zhuǎn)變?yōu)閯幽埽憩F(xiàn)為壩后跌水.綜合上下游斷面的流速均勻性指標(biāo)值來看,在下挑角為45°時,區(qū)域內(nèi)流速最為穩(wěn)定,水流的動勢能轉(zhuǎn)化較弱,從而對壩基安全構(gòu)成的威脅也降至最?。?/p>

    圖3 典型斷面所在位置Fig.3 Positions of the typical segments

    表1 不同挑角下實體丁壩上下游流態(tài)變化情況Tab.1 Change of the flow regime in the upstream and downstream of the solid spur-dike with different oriented angles

    3.2 實體丁壩不同布設(shè)排數(shù)下的挑流效果

    單丁壩對凹岸的防護作用有限,丁壩群的配置較單丁壩在挑離凹岸高速流團、增加凹岸低流區(qū)方面有更為優(yōu)異的工作效果.本文模擬的是45°挑角下不同布設(shè)排數(shù)下的挑流情況.由圖4可知,在沒有布設(shè)丁壩時,流線整體靠近凹岸,容易造成凹沖凸淤,河灣迅速向下移發(fā)展,使河道更加彎曲.在單排丁壩的作用下,散亂在凹岸的流線逐漸集中在河道中央,在壩前發(fā)生明顯轉(zhuǎn)向,繞過壩頭而下,部分沿壩頭轉(zhuǎn)向的水流逆流而上形成一定的回流區(qū).主流繞過壩頭一定距離后流線又開始向凹岸擴散.隨著丁壩排數(shù)的增加,擴散的主流線會重新聚集起來,歸順于主河道,從而減輕對凹岸的沖刷.

    在河灣段設(shè)置丁壩后,其挑流作用發(fā)揮,將沖擊凹岸的高速水流挑離至河道中心區(qū)域沖刷河道,整體高速流往凸岸偏移,靠近凸岸的淤積泥沙得以被沖刷.圖5是丁壩作用區(qū)的流速云圖,其中流速指沿水深平均值.從中可以直觀地看出,在沒有設(shè)置丁壩之前,紅色區(qū)塊的高速流團直擊凹岸,嚴重威脅河岸安全.在設(shè)置丁壩之后,高速流區(qū)(紅色部分)向河道中心偏移,并且隨著排數(shù)的增加,偏移的效果愈加明顯,同時靠近凹岸的低流區(qū)面積也隨之增加,說明壩后保護區(qū)的范圍隨排數(shù)增加呈擴大趨勢.另外值得一提的是,丁壩排數(shù)雖然增加,但是沖擊壩頭的最大流速卻降低.在單排丁壩作用下,沖擊壩頭最大流速為7 m/s,3排丁壩的情況下最大流速為6.5 m/s.可以反映出,較單丁壩而言,丁壩群體作用可以延長工程使用年限,減輕水流對壩頭的沖刷,工程整體工作更具安全性.

    (a) 排數(shù)為0

    3.3 實體丁壩與透水丁壩應(yīng)用模擬及對比分析

    實體丁壩與透水丁壩對于河道斷面水位h的影響效果相似,程度不同.模型按照挑角為45°,排數(shù)p=5布設(shè),水文要素的提取斷面布置在丁壩群上下游各20 m處,見圖6.如圖7所示,橫坐標(biāo)l為距右岸距離,丁壩上游的橫向水面線雍高明顯,凹岸水位雍高程度略大于凸岸,整體橫向水面線呈現(xiàn)凹岸高凸岸低的走勢;水流繞過壩體后出現(xiàn)跌水現(xiàn)象,下游橫向水面線降低,壩后水深較淺,尤其是對于實體丁壩,壩后水位降低更加明顯.可以看出相對于實體丁壩,透水丁壩在補充壩后水量這方面效果顯著,透水丁壩壩后水深較實體丁壩高0.8 m左右.出現(xiàn)上述水流現(xiàn)象的原因可能是:在凹岸設(shè)置丁壩之后,局部水面條件被改變,阻斷了部分水流過流,出現(xiàn)滯留雍阻現(xiàn)象,導(dǎo)致丁壩上游水位升高;另一方面,丁壩不同的結(jié)構(gòu)形式會對橫斷面上水面線的變化造成一定程度的影響.由于壩身的透水性,上游水流能通過壩間空隙對壩后水量進行持續(xù)補充,普遍增加了水流過壩后的凹岸水深,從而有效改善了壩后水深[9],水流從上游往下更為平緩地過渡.本研究結(jié)果符合文獻[9]中物理模型試驗所示現(xiàn)象.圖8展示的是丁壩上下游斷面上的流速分布情況,其中某點流速是指沿水深的平均值.在上游斷面處,無論是布置實體丁壩還是透水丁壩,流速同原河道相比均減小,并且在右岸流速減小的幅度更大,在遠離丁壩的一側(cè)流速逐漸一致.在下游斷面的壩根、壩身附近,流速均在1 m/s以下.壩頭附近,流速陡增,實體丁壩相鄰點流速最大增幅為69.7%,透水丁壩為60%.總體來看,透水丁壩周圍的流場更為平緩,水流流態(tài)也更好,有利于工程壽命的延長.

    (a) 排數(shù)為0

    圖6 典型斷面和測點所在位置Fig.6 Positions of the typical segments and measuring points

    (a) 上游斷面

    (a) 上游斷面

    丁壩的存在不僅影響河勢,同時也影響著含沙量的重新分配.圖9展示了渾水條件下兩類丁壩在工作時刻,壩群周圍的懸沙濃度ρ分布情況,其中懸沙濃度代表的是沿水深的平均值.壩區(qū)的含沙量明顯低于主河道內(nèi)的含沙量.從壩區(qū)至主流區(qū),實體丁壩泥沙濃度最大變化幅度為35.5%,透水丁壩為30.3%.同原河道相比,設(shè)置兩種類型丁壩后,壩后泥沙均得到很大程度的淤積,實體丁壩壩后含沙量從約30 kg/m3下降到最低 24.72 kg/m3,透水丁壩則下降至24.08 kg/m3.而凸岸和主河道淤積的泥沙被沖刷,懸沙濃度增加,這種懸沙濃度的重新分配情況表征了布設(shè)丁壩之后,河灣段不再處于凹沖凸淤、河岸下移的不利情勢.根據(jù)丁壩調(diào)整的流路,河灣段會逐漸恢復(fù)到順直的有利狀態(tài).泥沙濃度變化趨勢與已有的一項插板透水丁壩分流分沙模擬研究中的變化趨勢一致[10],進一步驗證了模擬結(jié)果的可靠性.

    圖9 丁壩群附近懸沙濃度分布Fig.9 Distribution of the sediment concentration near the spur-dike group

    4 結(jié)果討論

    4.1 挑角和排數(shù)對實體丁壩導(dǎo)流影響機理分析

    單丁壩的河道整治效果有局限性,在實際治理中丁壩大多以丁壩群或者雙丁壩的形式出現(xiàn).從本文中的單丁壩與壩群的挑流效果對比圖(圖4、5)可以直觀看出,丁壩群發(fā)揮的作用并非是單丁壩作用的簡單累加.在單丁壩的挑流作用下,原本貼近河岸的水流被束窄,集中到河道中央.隨著過壩后過流面積增大,水流以拋物線形式擴散,再次流落到河岸邊繼續(xù)沖刷.而當(dāng)丁壩群作為系統(tǒng)工作時,水流先被第一座丁壩挑離開河岸偏向河道中間一定范圍,流線還未擴散至河岸又在第二座丁壩的作用下被二次挑流,水流再次向河道中心聚攏,循環(huán)往復(fù)直到尾壩發(fā)揮完作用.由于前排的丁壩往往為后排的丁壩提供“掩護”,合理分擔(dān)水勢,所以丁壩群壩頭受沖刷強度不如單個丁壩作用時強烈.丁壩的挑角關(guān)系到壩面的迎沖強度,也間接決定了河床的束窄程度.已有研究經(jīng)過模型試驗,通過對不同挑角下丁壩附近水流結(jié)構(gòu)分析認為:挑角越大,丁壩對水流的挑流作用越強烈,在壩頭附近的流速變幅越大,同時壩頭旋渦強度也越大,且在30°挑角時水流最為穩(wěn)定[11].本文分析發(fā)現(xiàn)在不同的斷面位置,流速分布的均勻性隨挑角增大的變化規(guī)律不一致.在上游斷面,流速均勻性隨挑角的增大而減小,30°最穩(wěn)定;在下游斷面隨挑角的增大而增大,75°最穩(wěn)定.可以從以下兩方面來解釋產(chǎn)生差異的原因:(1)所述的模型試驗[11]是基于順直河道,而本文研究的是黃河下游游蕩河道,河灣處水流流態(tài)和水力條件十分復(fù)雜.(2)模型試驗觀察的是限于壩頭處的流速變幅,而本文是針對丁壩上下游作用區(qū).不過結(jié)合實際工程經(jīng)驗,建議丁壩的挑角均設(shè)計為下挑45°為宜.

    4.2 丁壩群透水性對緩流落淤的影響

    圖10展示的是應(yīng)用實體丁壩和透水丁壩后凹岸的流速狀態(tài),測點處的流速是指沿水深平均值,具體測點位置見圖6.從圖10中可以看出丁壩發(fā)揮挑流作用,致使凹岸流速大幅度降低,起到護灘保堤的作用,降低了高速水流在離心力作用下對凹岸造成的強烈沖刷.實體丁壩對河段凹岸流速減小平均貢獻率為85%,透水丁壩為84%,兩者對凹岸都具有很顯著的緩流效果,實體丁壩的緩流效果略好一些.不過通過改變透水丁壩的透水率和透水孔數(shù),可以接近甚至超過應(yīng)用實體丁壩時所達到的緩流效果.而且兩者在促進壩后落淤方面的效果也很接近,沉積的泥沙不僅為壩體提供支撐,也在一定程度上填補了凹岸的形態(tài)空缺,逐漸“裁灣”,抑制河道凹岸繼續(xù)下滑趨勢.通過對比分析兩者的模擬試驗結(jié)果發(fā)現(xiàn),在黃河下游河灣段應(yīng)用透水丁壩群的緩流落淤效果與實體丁壩群基本相同,并未因部分透水而使其治理效果減弱.從理論上來說,只要透水丁壩的透水率不是太大,兩者的治理效果也不應(yīng)該相差過大.究其原因是整治工程的控導(dǎo)作用主要取決于水流動力條件和工程的邊界條件.前者主要包括流量、入流角等因子,后者則包括工程的長度、平面形式、消能大小等[12].由于控制變量,在模型設(shè)置中,實體丁壩和透水丁壩有差異的僅是平面形式,導(dǎo)致消能大小有所不同.實體丁壩通過水流撞擊密閉平面消能,透水丁壩通過透水孔減少工程前大河流量,其作用同樣降低水流動能.因此可推論,只要透水率不是大到使大河流量銳減的程度,兩者的治理作用應(yīng)該是相近的.反觀透水丁壩因為它的透水性,能維持水沙的縱向交換,更加符合柔性的生態(tài)調(diào)節(jié)原則,且不影響洪水過流,更加貼近自然調(diào)節(jié)的管理理念.

    圖10 凹岸測點流速Fig.10 Velocity at measuring points near the concave bank

    4.3 需要注意和改進的問題

    自20世紀70年代以來,透水丁壩因其良好的適應(yīng)性和功能性已經(jīng)廣泛應(yīng)用于河道治理、工程防護等諸多領(lǐng)域[13].其主要形式有水力插板透水丁壩、井柱網(wǎng)格透水丁壩、鋼筋混凝土透水樁壩和井柱鋼絲籠壩等.其中鋼筋混凝土透水樁壩在近幾十年黃河下游河道整治工程中的應(yīng)用逐漸增多,如河南段東安工程和韋灘工程等.然而鋼筋混凝土透水樁壩存在易斷樁、易沖毀、施工質(zhì)量難以把控等缺陷.例如2021年4月東安控導(dǎo)單排透水樁壩在經(jīng)過水流頂沖后,工程中間沖毀80 m以上.東安工程的水毀問題說明目前亟須尋找可靠的設(shè)計方法滿足工程安全性的需求.透水丁壩在設(shè)計和施工技術(shù)方面存在相關(guān)規(guī)范較少的問題,為了使其設(shè)計和施工達到標(biāo)準化,需要制定透水丁壩的工程設(shè)計指南以及施工技術(shù)規(guī)范,并且對透水丁壩的水流特性、沖淤特性和穩(wěn)定性計算等進行更加深入的研究.另一方面,為了深入貫徹落實生態(tài)文明建設(shè)的理念,加強流域生態(tài)保護,迫切需要對河道工程進行仿生態(tài)改造結(jié)構(gòu)形式與新材料的相關(guān)研究,在滿足防洪需求壩型的基礎(chǔ)上,分析有利于水生生物繁衍生息的河道工程材料和結(jié)構(gòu).并且要求合理選擇壩址,利用沖淤環(huán)境,為水生動植物創(chuàng)造良好的生境條件.今后研究將向構(gòu)建生態(tài)友好型透水丁壩發(fā)展,尋找這種兼顧防洪安全與生態(tài)友好的建設(shè)措施.

    5 結(jié) 論

    (1)挑角直接影響丁壩周圍自由水面及斷面流速分布的均勻性.在丁壩上游斷面,挑角為75°時,水位變化最大,壩前雍高1.388 m,此時流速分布的均勻性也最差.但在下游斷面卻是挑角為30°時流速分布的均勻性最差.由于河灣處水沙機理復(fù)雜,綜合上下游斷面的流速均勻性指標(biāo)值來看,在下挑角為45°時,區(qū)域內(nèi)流速分布最為穩(wěn)定,水面的流態(tài)較緩,對岸基的動態(tài)壓力相對較小,對河灣能起到更為理想的防護作用.

    (2)較單丁壩而言,壩群具有更優(yōu)的送溜挑流效果,壩群的存在能將沖擊凹岸的高速水流挑離至河道中心,保護岸線.丁壩群體作用可以延長工程使用年限,減輕水流對壩頭的沖刷.單丁壩作用下沖擊壩頭流速最大可達7 m/s,而群體作用下最大流速為6.5 m/s.因此適宜的排數(shù)不僅能有效控制河勢,也有利于提高壩群整體的安全性能.

    (3)在黃河下游河灣段應(yīng)用透水丁壩群的緩流落淤效果與實體丁壩群相差不大.實體丁壩壩后懸沙濃度從30 kg/m3下降到最低24.72 kg/m3,透水丁壩則下降至24.08 kg/m3;兩者對河段凹岸流速減小平均貢獻率均高達80%以上,具有很強的護岸效果.此外透水丁壩因為其獨具的透水性能有利于水生生物繁衍生息,兼顧了防洪安全和生態(tài)友好,具有廣闊的發(fā)展前景.

    猜你喜歡
    工程
    《工程爆破》》征稿征訂
    工程爆破(2022年3期)2022-07-26 01:58:56
    《工程爆破》征稿簡則
    工程爆破(2022年2期)2022-06-17 14:13:56
    子午工程
    太空探索(2016年6期)2016-07-10 12:09:06
    工程
    工程
    工程
    工程
    工程
    工程
    工程
    秋霞在线观看毛片| 夫妻午夜视频| av播播在线观看一区| 男的添女的下面高潮视频| 夜夜看夜夜爽夜夜摸| 人体艺术视频欧美日本| a级一级毛片免费在线观看| 午夜免费观看性视频| 国产一级毛片在线| 午夜激情福利司机影院| 国产 精品1| 日日摸夜夜添夜夜爱| xxx大片免费视频| 成年人免费黄色播放视频 | 国产熟女欧美一区二区| 欧美区成人在线视频| 亚洲成色77777| 在线精品无人区一区二区三| 麻豆成人av视频| 国产伦理片在线播放av一区| 亚洲,欧美,日韩| 黑人巨大精品欧美一区二区蜜桃 | 欧美日韩精品成人综合77777| 亚洲一区二区三区欧美精品| 最新的欧美精品一区二区| 伊人亚洲综合成人网| 久久久午夜欧美精品| 欧美精品一区二区大全| 欧美三级亚洲精品| 亚洲成人av在线免费| www.av在线官网国产| 精品酒店卫生间| 久久久国产精品麻豆| 中文字幕人妻丝袜制服| 在线看a的网站| 黑丝袜美女国产一区| 免费观看av网站的网址| 99视频精品全部免费 在线| 久久久久久久久久久免费av| 岛国毛片在线播放| 免费少妇av软件| 国产爽快片一区二区三区| 老女人水多毛片| 日韩电影二区| 国产黄片美女视频| 欧美亚洲 丝袜 人妻 在线| 亚洲国产av新网站| 一级a做视频免费观看| 国产免费视频播放在线视频| 精品一区二区三区视频在线| 久久久久久久大尺度免费视频| 免费少妇av软件| 亚洲欧美日韩东京热| 永久免费av网站大全| 久久久久久久久大av| 人人妻人人添人人爽欧美一区卜| 亚洲精品国产av成人精品| 国产精品一二三区在线看| 久久久久久久亚洲中文字幕| 国产一区二区三区综合在线观看 | 精品亚洲成国产av| 在线精品无人区一区二区三| 午夜福利网站1000一区二区三区| 日本与韩国留学比较| 中文精品一卡2卡3卡4更新| 亚洲,一卡二卡三卡| 日韩一区二区三区影片| 免费人成在线观看视频色| 精品国产一区二区三区久久久樱花| 五月伊人婷婷丁香| 国产视频首页在线观看| 精品一区二区三区视频在线| av在线老鸭窝| 成人亚洲欧美一区二区av| 精华霜和精华液先用哪个| 欧美少妇被猛烈插入视频| 午夜老司机福利剧场| 亚洲精品国产av蜜桃| 伦精品一区二区三区| 国内少妇人妻偷人精品xxx网站| 人人妻人人澡人人看| 色视频在线一区二区三区| 简卡轻食公司| 精品亚洲成a人片在线观看| 五月开心婷婷网| 自拍偷自拍亚洲精品老妇| 亚洲国产欧美在线一区| 嫩草影院入口| 日本黄色片子视频| 一级毛片aaaaaa免费看小| 免费观看在线日韩| 最近最新中文字幕免费大全7| 在线观看免费高清a一片| 三级国产精品片| av视频免费观看在线观看| 男女免费视频国产| 成人国产av品久久久| 精品人妻偷拍中文字幕| 我的女老师完整版在线观看| 高清不卡的av网站| 欧美丝袜亚洲另类| 国产伦理片在线播放av一区| 久久久亚洲精品成人影院| 自拍欧美九色日韩亚洲蝌蚪91 | 永久免费av网站大全| 黑人高潮一二区| 卡戴珊不雅视频在线播放| 嘟嘟电影网在线观看| 久久久欧美国产精品| 寂寞人妻少妇视频99o| 一本—道久久a久久精品蜜桃钙片| 精品少妇久久久久久888优播| 欧美 亚洲 国产 日韩一| 99热全是精品| 在线观看人妻少妇| 国产精品99久久久久久久久| 精品国产露脸久久av麻豆| 一区二区三区精品91| 一级毛片 在线播放| 久久综合国产亚洲精品| 成人无遮挡网站| 欧美丝袜亚洲另类| 久久精品熟女亚洲av麻豆精品| 日韩熟女老妇一区二区性免费视频| 婷婷色麻豆天堂久久| 18禁裸乳无遮挡动漫免费视频| 免费观看性生交大片5| 欧美日韩av久久| 特大巨黑吊av在线直播| 国产深夜福利视频在线观看| 又黄又爽又刺激的免费视频.| 日本黄色日本黄色录像| 你懂的网址亚洲精品在线观看| 伊人亚洲综合成人网| 亚洲伊人久久精品综合| 色哟哟·www| 亚洲,一卡二卡三卡| 国产黄色视频一区二区在线观看| 又爽又黄a免费视频| 国产毛片在线视频| 人体艺术视频欧美日本| 少妇被粗大猛烈的视频| 欧美国产精品一级二级三级 | 岛国毛片在线播放| 22中文网久久字幕| 极品人妻少妇av视频| 久久综合国产亚洲精品| 五月天丁香电影| 91成人精品电影| 精品酒店卫生间| 汤姆久久久久久久影院中文字幕| 一本—道久久a久久精品蜜桃钙片| 黄色怎么调成土黄色| 国产精品嫩草影院av在线观看| 热re99久久精品国产66热6| 精品一品国产午夜福利视频| 国产黄片美女视频| 国产真实伦视频高清在线观看| 2018国产大陆天天弄谢| 国模一区二区三区四区视频| 国产 精品1| 国产成人精品一,二区| 久久久久久久久久久免费av| av在线观看视频网站免费| xxx大片免费视频| 亚洲高清免费不卡视频| 欧美日韩av久久| 国产日韩欧美在线精品| www.av在线官网国产| av播播在线观看一区| 久久久精品94久久精品| 国产精品嫩草影院av在线观看| 伊人久久国产一区二区| 国产精品99久久久久久久久| 亚洲精品乱码久久久久久按摩| 亚洲精品色激情综合| 精品少妇内射三级| 欧美 日韩 精品 国产| 日韩精品免费视频一区二区三区 | 一本—道久久a久久精品蜜桃钙片| 日本爱情动作片www.在线观看| 久久久久国产精品人妻一区二区| videos熟女内射| 免费大片黄手机在线观看| 美女大奶头黄色视频| 爱豆传媒免费全集在线观看| 六月丁香七月| 精品99又大又爽又粗少妇毛片| 国产伦理片在线播放av一区| 男女边摸边吃奶| 99久久中文字幕三级久久日本| 久久 成人 亚洲| 欧美3d第一页| 少妇的逼水好多| 中文天堂在线官网| 亚洲av男天堂| 青青草视频在线视频观看| 麻豆成人午夜福利视频| 国产91av在线免费观看| 国产精品99久久久久久久久| 如日韩欧美国产精品一区二区三区 | 51国产日韩欧美| 性色av一级| 国产精品国产av在线观看| 一级毛片黄色毛片免费观看视频| 一本—道久久a久久精品蜜桃钙片| 男男h啪啪无遮挡| 亚洲人成网站在线播| 一个人看视频在线观看www免费| 91午夜精品亚洲一区二区三区| 国产亚洲欧美精品永久| 色婷婷久久久亚洲欧美| av黄色大香蕉| 成人免费观看视频高清| 日本色播在线视频| 亚洲精品久久午夜乱码| 极品教师在线视频| 内射极品少妇av片p| 蜜臀久久99精品久久宅男| 亚洲国产欧美在线一区| 日本欧美国产在线视频| 国产一区有黄有色的免费视频| 国产有黄有色有爽视频| 精品久久久久久电影网| 黄色毛片三级朝国网站 | 国产精品成人在线| 黑人猛操日本美女一级片| 51国产日韩欧美| av视频免费观看在线观看| 三级国产精品欧美在线观看| 日本猛色少妇xxxxx猛交久久| 高清欧美精品videossex| 一级a做视频免费观看| 日韩大片免费观看网站| 伊人久久国产一区二区| 国产精品国产av在线观看| 女人久久www免费人成看片| 日本wwww免费看| videossex国产| 黑人巨大精品欧美一区二区蜜桃 | 久久久久久久国产电影| 久久久国产欧美日韩av| 亚洲成色77777| 一区二区av电影网| 卡戴珊不雅视频在线播放| 熟妇人妻不卡中文字幕| 亚洲欧美日韩东京热| 黄片无遮挡物在线观看| 美女xxoo啪啪120秒动态图| 欧美丝袜亚洲另类| 日日摸夜夜添夜夜添av毛片| 国产精品福利在线免费观看| 国产精品女同一区二区软件| 国产视频内射| 亚洲无线观看免费| 婷婷色麻豆天堂久久| 亚洲精品国产色婷婷电影| 日韩在线高清观看一区二区三区| 国产白丝娇喘喷水9色精品| 精品国产乱码久久久久久小说| 国精品久久久久久国模美| 久久久久久久国产电影| 亚洲经典国产精华液单| 99久久中文字幕三级久久日本| 高清欧美精品videossex| 欧美一级a爱片免费观看看| 麻豆精品久久久久久蜜桃| 亚洲精品乱久久久久久| 我要看黄色一级片免费的| 国产乱人偷精品视频| 国产有黄有色有爽视频| 亚洲国产精品专区欧美| 国语对白做爰xxxⅹ性视频网站| 久久精品国产亚洲av天美| 国产伦理片在线播放av一区| 国产真实伦视频高清在线观看| 国产在视频线精品| 国产伦精品一区二区三区四那| 蜜臀久久99精品久久宅男| 综合色丁香网| 久久精品国产a三级三级三级| 人体艺术视频欧美日本| 久久久久视频综合| 麻豆成人午夜福利视频| 欧美高清成人免费视频www| av免费观看日本| 99久久综合免费| 精品久久久噜噜| 日日啪夜夜撸| 天美传媒精品一区二区| 亚洲图色成人| 老司机影院毛片| 一级毛片 在线播放| av不卡在线播放| 日日摸夜夜添夜夜添av毛片| 哪个播放器可以免费观看大片| 六月丁香七月| 99热这里只有是精品50| 久久热精品热| 噜噜噜噜噜久久久久久91| 91精品国产九色| 亚洲中文av在线| 国产精品人妻久久久影院| 丝袜脚勾引网站| 91精品一卡2卡3卡4卡| 国产精品福利在线免费观看| 99九九线精品视频在线观看视频| 青春草国产在线视频| 九草在线视频观看| 国产精品国产三级国产av玫瑰| 国产淫片久久久久久久久| 久久毛片免费看一区二区三区| 婷婷色综合大香蕉| 日韩欧美精品免费久久| 国产午夜精品一二区理论片| 精品久久久久久久久亚洲| 欧美日韩国产mv在线观看视频| 午夜影院在线不卡| 只有这里有精品99| 亚洲欧美日韩卡通动漫| av在线app专区| 欧美三级亚洲精品| 国产精品一区二区在线不卡| 纯流量卡能插随身wifi吗| 夫妻性生交免费视频一级片| 汤姆久久久久久久影院中文字幕| 自拍欧美九色日韩亚洲蝌蚪91 | 蜜桃久久精品国产亚洲av| 波野结衣二区三区在线| 亚洲精品色激情综合| 免费黄色在线免费观看| 在线精品无人区一区二区三| 日韩制服骚丝袜av| 男人和女人高潮做爰伦理| 国产精品不卡视频一区二区| av福利片在线| 日本欧美国产在线视频| 国产精品不卡视频一区二区| 亚洲熟女精品中文字幕| 少妇猛男粗大的猛烈进出视频| 国产爽快片一区二区三区| 国产免费一级a男人的天堂| 99九九线精品视频在线观看视频| 三上悠亚av全集在线观看 | 国产免费一区二区三区四区乱码| 一级爰片在线观看| 亚洲精品自拍成人| 久久韩国三级中文字幕| 国产精品麻豆人妻色哟哟久久| 午夜日本视频在线| 国产高清有码在线观看视频| 九九爱精品视频在线观看| 午夜福利在线观看免费完整高清在| 日韩强制内射视频| 国产精品伦人一区二区| av天堂中文字幕网| 国产精品三级大全| 性色av一级| 日本免费在线观看一区| 欧美国产精品一级二级三级 | 成人黄色视频免费在线看| 国语对白做爰xxxⅹ性视频网站| 另类亚洲欧美激情| 99热国产这里只有精品6| 色婷婷久久久亚洲欧美| 午夜福利,免费看| 亚州av有码| 免费看日本二区| 在线看a的网站| 亚洲伊人久久精品综合| 久久韩国三级中文字幕| 日韩欧美一区视频在线观看 | 久久99热6这里只有精品| 国产黄片美女视频| 亚洲欧洲日产国产| 中文字幕免费在线视频6| 欧美精品一区二区大全| 91在线精品国自产拍蜜月| 中文字幕久久专区| 婷婷色综合大香蕉| 亚洲欧美中文字幕日韩二区| 亚洲久久久国产精品| 18禁裸乳无遮挡动漫免费视频| 国产亚洲一区二区精品| 九九在线视频观看精品| av天堂久久9| 老女人水多毛片| 国产亚洲5aaaaa淫片| 六月丁香七月| 欧美高清成人免费视频www| 日本vs欧美在线观看视频 | 午夜日本视频在线| 精品少妇内射三级| 少妇被粗大的猛进出69影院 | 国产精品99久久99久久久不卡 | 亚洲第一av免费看| 深夜a级毛片| 我要看黄色一级片免费的| 内射极品少妇av片p| 亚洲av福利一区| 妹子高潮喷水视频| 久久99一区二区三区| 免费观看性生交大片5| 免费少妇av软件| 婷婷色综合大香蕉| 99久国产av精品国产电影| 国产毛片在线视频| xxx大片免费视频| 在线亚洲精品国产二区图片欧美 | 日产精品乱码卡一卡2卡三| 啦啦啦视频在线资源免费观看| h视频一区二区三区| 少妇人妻久久综合中文| 成人二区视频| 亚洲精品aⅴ在线观看| 亚洲高清免费不卡视频| kizo精华| 国产成人精品婷婷| 久久精品国产亚洲av天美| 久久99蜜桃精品久久| 亚洲欧美中文字幕日韩二区| 三级经典国产精品| 成人特级av手机在线观看| 亚洲第一av免费看| 一区二区av电影网| 国产亚洲一区二区精品| 永久免费av网站大全| 成人无遮挡网站| 国产精品.久久久| 青春草亚洲视频在线观看| 久久婷婷青草| 一本—道久久a久久精品蜜桃钙片| 中文字幕制服av| 日产精品乱码卡一卡2卡三| 男人爽女人下面视频在线观看| 国产男女超爽视频在线观看| 亚洲色图综合在线观看| 中文资源天堂在线| 高清av免费在线| 老女人水多毛片| 亚洲美女视频黄频| 免费看av在线观看网站| 亚洲精品亚洲一区二区| 中文乱码字字幕精品一区二区三区| 夫妻性生交免费视频一级片| 久久综合国产亚洲精品| av视频免费观看在线观看| 在线观看三级黄色| 精品99又大又爽又粗少妇毛片| 一级爰片在线观看| 在线天堂最新版资源| 日日爽夜夜爽网站| 日日摸夜夜添夜夜添av毛片| 一边亲一边摸免费视频| 又爽又黄a免费视频| 日本午夜av视频| 麻豆精品久久久久久蜜桃| 日韩一本色道免费dvd| 我要看黄色一级片免费的| 国产黄片美女视频| 久久精品久久精品一区二区三区| 久久亚洲国产成人精品v| 欧美国产精品一级二级三级 | 一本色道久久久久久精品综合| 内地一区二区视频在线| 国产午夜精品久久久久久一区二区三区| 久久久久国产精品人妻一区二区| 亚洲国产精品一区二区三区在线| 六月丁香七月| 精品酒店卫生间| 日本av手机在线免费观看| 久久久欧美国产精品| 日本猛色少妇xxxxx猛交久久| 国产在线一区二区三区精| 哪个播放器可以免费观看大片| 99热这里只有是精品在线观看| 国产永久视频网站| 欧美 亚洲 国产 日韩一| 欧美区成人在线视频| 美女内射精品一级片tv| 亚洲在久久综合| 人人澡人人妻人| 中国美白少妇内射xxxbb| 国产视频首页在线观看| 午夜视频国产福利| 国产精品三级大全| 久久精品熟女亚洲av麻豆精品| 自拍欧美九色日韩亚洲蝌蚪91 | 亚洲精品色激情综合| 日本91视频免费播放| av免费在线看不卡| 亚洲精品乱码久久久v下载方式| 哪个播放器可以免费观看大片| 妹子高潮喷水视频| 只有这里有精品99| 99热6这里只有精品| 亚洲av中文av极速乱| 国产中年淑女户外野战色| a级毛片在线看网站| 麻豆乱淫一区二区| 男女啪啪激烈高潮av片| 国产老妇伦熟女老妇高清| 高清黄色对白视频在线免费看 | 全区人妻精品视频| 中文字幕精品免费在线观看视频 | 香蕉精品网在线| 黄片无遮挡物在线观看| 久久人人爽人人片av| 夜夜看夜夜爽夜夜摸| 亚洲综合精品二区| 亚洲国产成人一精品久久久| 欧美精品一区二区大全| 老司机影院成人| 在线看a的网站| 王馨瑶露胸无遮挡在线观看| 亚洲成色77777| 亚洲欧美清纯卡通| 熟女电影av网| 我要看黄色一级片免费的| 亚洲第一区二区三区不卡| 亚洲四区av| 国产成人精品久久久久久| 大香蕉97超碰在线| 国产av码专区亚洲av| 中文字幕精品免费在线观看视频 | 一级av片app| av福利片在线| 九九爱精品视频在线观看| 亚洲美女视频黄频| 久久久久久久久久久丰满| 一级毛片 在线播放| 国产日韩欧美在线精品| 国产男人的电影天堂91| 日韩精品免费视频一区二区三区 | 亚洲不卡免费看| 菩萨蛮人人尽说江南好唐韦庄| 老熟女久久久| 国产伦精品一区二区三区视频9| 亚洲图色成人| 我的女老师完整版在线观看| 又黄又爽又刺激的免费视频.| 亚洲国产欧美在线一区| 一个人免费看片子| 各种免费的搞黄视频| 久久久久久久亚洲中文字幕| 一本色道久久久久久精品综合| 精品一区二区三卡| 国国产精品蜜臀av免费| 中文字幕久久专区| 国产男女内射视频| 亚洲人成网站在线观看播放| 亚洲精品国产成人久久av| 熟女av电影| 亚洲精品久久午夜乱码| 伦理电影大哥的女人| 一级毛片我不卡| 亚洲国产精品一区二区三区在线| av女优亚洲男人天堂| 一级,二级,三级黄色视频| 三级国产精品片| 如日韩欧美国产精品一区二区三区 | 久久精品熟女亚洲av麻豆精品| 日韩三级伦理在线观看| 美女中出高潮动态图| 免费久久久久久久精品成人欧美视频 | 赤兔流量卡办理| 亚洲国产精品一区三区| 99re6热这里在线精品视频| 啦啦啦视频在线资源免费观看| 黄色配什么色好看| 亚洲国产最新在线播放| 精品熟女少妇av免费看| 男人爽女人下面视频在线观看| 国产黄色免费在线视频| 一级av片app| 亚洲怡红院男人天堂| 国产色爽女视频免费观看| 啦啦啦在线观看免费高清www| 亚洲人成网站在线观看播放| 天美传媒精品一区二区| 成人亚洲精品一区在线观看| 免费观看的影片在线观看| 纵有疾风起免费观看全集完整版| 狂野欧美激情性bbbbbb| 国产黄频视频在线观看| 久久久久久人妻| 大码成人一级视频| 久久久精品免费免费高清| 在线观看美女被高潮喷水网站| 大香蕉久久网| 国产亚洲精品久久久com| 久久热精品热| 成人二区视频| 国产成人91sexporn| 黄色欧美视频在线观看| 久久ye,这里只有精品| 99国产精品免费福利视频| 精品国产国语对白av| 人妻夜夜爽99麻豆av| 插阴视频在线观看视频| 欧美精品一区二区大全| 在线观看www视频免费| 日韩精品有码人妻一区| 国产免费又黄又爽又色| 久久影院123| 国产白丝娇喘喷水9色精品| 在现免费观看毛片| 免费人妻精品一区二区三区视频| 久热久热在线精品观看| 中国三级夫妇交换| 少妇的逼水好多| 一二三四中文在线观看免费高清| 色网站视频免费| 久久久久国产精品人妻一区二区| 观看av在线不卡|