• <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竹,舟*,瑩,
    大連理工大學學報 2022年4期
    關鍵詞:工程

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

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

    0 引 言

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

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

    1 區(qū)域概況

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

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

    2 數值模擬

    2.1 模型計算范圍和數據來源

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

    2.2 模擬試驗方案

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

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

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

    (1)

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

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

    (2)

    (3)

    3 試驗結果

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

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

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

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

    (a) 排數為0

    3.3 實體丁壩與透水丁壩應用模擬及對比分析

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

    (a) 排數為0

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

    (a) 上游斷面

    (a) 上游斷面

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

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

    4 結果討論

    4.1 挑角和排數對實體丁壩導流影響機理分析

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

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

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

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

    4.3 需要注意和改進的問題

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

    5 結 論

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

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

    (3)在黃河下游河灣段應用透水丁壩群的緩流落淤效果與實體丁壩群相差不大.實體丁壩壩后懸沙濃度從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
    工程
    工程
    工程
    工程
    工程
    工程
    工程
    久久精品aⅴ一区二区三区四区| 欧美日韩亚洲国产一区二区在线观看| 99久久综合精品五月天人人| 怎么达到女性高潮| 99精品在免费线老司机午夜| 欧美国产精品va在线观看不卡| 久热爱精品视频在线9| 精品久久久久久久久久免费视频| 亚洲第一电影网av| 国内精品久久久久精免费| 日韩欧美国产在线观看| 国产在线精品亚洲第一网站| 国产亚洲精品一区二区www| 一本大道久久a久久精品| av天堂久久9| 国产成人av教育| 国产麻豆成人av免费视频| 免费观看人在逋| 亚洲va日本ⅴa欧美va伊人久久| 午夜日韩欧美国产| 久久久久久免费高清国产稀缺| 嫩草影视91久久| 免费少妇av软件| 91麻豆精品激情在线观看国产| 日韩三级视频一区二区三区| 欧美日韩亚洲综合一区二区三区_| 欧美激情高清一区二区三区| 91国产中文字幕| 精品免费久久久久久久清纯| 国产国语露脸激情在线看| 国产麻豆69| 亚洲一区高清亚洲精品| 亚洲 欧美一区二区三区| 好看av亚洲va欧美ⅴa在| 免费无遮挡裸体视频| 悠悠久久av| 无遮挡黄片免费观看| 国产精品,欧美在线| 女生性感内裤真人,穿戴方法视频| 国产伦人伦偷精品视频| 欧美乱妇无乱码| 丝袜人妻中文字幕| 国产精品秋霞免费鲁丝片| 成熟少妇高潮喷水视频| 男男h啪啪无遮挡| 色婷婷久久久亚洲欧美| 91字幕亚洲| 午夜a级毛片| 啦啦啦 在线观看视频| 欧美日韩亚洲综合一区二区三区_| 丝袜人妻中文字幕| 熟女少妇亚洲综合色aaa.| 国产精品电影一区二区三区| 曰老女人黄片| 不卡一级毛片| 久久亚洲精品不卡| 国产精品久久久人人做人人爽| 美女高潮喷水抽搐中文字幕| 国产高清有码在线观看视频 | 不卡一级毛片| 欧美黑人精品巨大| 长腿黑丝高跟| 人成视频在线观看免费观看| 亚洲第一欧美日韩一区二区三区| 一级毛片女人18水好多| 一级黄色大片毛片| 美女大奶头视频| 亚洲av美国av| 亚洲国产欧美一区二区综合| 大码成人一级视频| 一级毛片精品| 亚洲中文字幕日韩| 一区二区日韩欧美中文字幕| 国产av精品麻豆| 在线观看午夜福利视频| 女人精品久久久久毛片| 在线播放国产精品三级| 亚洲aⅴ乱码一区二区在线播放 | 黑人巨大精品欧美一区二区mp4| 黄色 视频免费看| 亚洲自拍偷在线| 欧美老熟妇乱子伦牲交| 亚洲av第一区精品v没综合| 在线观看日韩欧美| 9色porny在线观看| 色婷婷久久久亚洲欧美| 亚洲自拍偷在线| 亚洲久久久国产精品| 久久人妻熟女aⅴ| 中文字幕人妻熟女乱码| 最新美女视频免费是黄的| 黑人操中国人逼视频| 精品不卡国产一区二区三区| 男女下面进入的视频免费午夜 | 美女高潮到喷水免费观看| 国产av又大| 最近最新中文字幕大全电影3 | 不卡一级毛片| 欧美日韩中文字幕国产精品一区二区三区 | 欧美激情 高清一区二区三区| 一级,二级,三级黄色视频| 午夜福利成人在线免费观看| 国产精品免费视频内射| 岛国视频午夜一区免费看| 这个男人来自地球电影免费观看| 精品一区二区三区av网在线观看| 国语自产精品视频在线第100页| 一区二区三区高清视频在线| 在线免费观看的www视频| 叶爱在线成人免费视频播放| 99在线视频只有这里精品首页| 午夜久久久在线观看| 色综合婷婷激情| 中文亚洲av片在线观看爽| 在线视频色国产色| 免费高清在线观看日韩| 国产国语露脸激情在线看| 成人18禁高潮啪啪吃奶动态图| 国产亚洲精品av在线| 亚洲色图av天堂| 乱人伦中国视频| 99热只有精品国产| 亚洲精品久久国产高清桃花| 午夜成年电影在线免费观看| 国产免费男女视频| 婷婷精品国产亚洲av在线| 亚洲熟妇中文字幕五十中出| 99国产精品一区二区蜜桃av| 国产97色在线日韩免费| 99久久99久久久精品蜜桃| 中文亚洲av片在线观看爽| 男人舔女人下体高潮全视频| 高清黄色对白视频在线免费看| 99香蕉大伊视频| 欧美中文综合在线视频| 国产亚洲精品一区二区www| 成人av一区二区三区在线看| 免费高清在线观看日韩| 首页视频小说图片口味搜索| 又黄又爽又免费观看的视频| 欧美午夜高清在线| av网站免费在线观看视频| 国产区一区二久久| 久久精品国产综合久久久| 不卡一级毛片| 变态另类丝袜制服| 操出白浆在线播放| 看免费av毛片| 亚洲国产高清在线一区二区三 | 精品国产超薄肉色丝袜足j| 国产精品久久久av美女十八| 亚洲人成77777在线视频| 如日韩欧美国产精品一区二区三区| 亚洲avbb在线观看| 久久精品国产综合久久久| 色哟哟哟哟哟哟| 两个人免费观看高清视频| 叶爱在线成人免费视频播放| 精品国产一区二区久久| 在线国产一区二区在线| 麻豆久久精品国产亚洲av| 男女做爰动态图高潮gif福利片 | 国产麻豆成人av免费视频| 两个人看的免费小视频| 精品一区二区三区视频在线观看免费| 国产精品影院久久| 露出奶头的视频| 亚洲一卡2卡3卡4卡5卡精品中文| 免费一级毛片在线播放高清视频 | 亚洲 欧美 日韩 在线 免费| 免费少妇av软件| 免费看美女性在线毛片视频| 国产成人av教育| 窝窝影院91人妻| 久久婷婷成人综合色麻豆| 亚洲成人国产一区在线观看| 欧美成人一区二区免费高清观看 | 91麻豆av在线| 亚洲色图av天堂| 波多野结衣高清无吗| 国产精品免费一区二区三区在线| 中国美女看黄片| 国产成人啪精品午夜网站| 少妇熟女aⅴ在线视频| 在线视频色国产色| 免费看美女性在线毛片视频| 女人精品久久久久毛片| 久久精品国产综合久久久| 91字幕亚洲| 国产一区二区在线av高清观看| 亚洲成av人片免费观看| 国产免费男女视频| АⅤ资源中文在线天堂| 老司机午夜福利在线观看视频| 欧美日本视频| 美女午夜性视频免费| 久久中文字幕人妻熟女| 欧美一级毛片孕妇| 久久精品亚洲熟妇少妇任你| 嫩草影视91久久| 在线天堂中文资源库| 亚洲精华国产精华精| 少妇的丰满在线观看| 正在播放国产对白刺激| 熟妇人妻久久中文字幕3abv| 亚洲成人精品中文字幕电影| 一级,二级,三级黄色视频| 成人18禁在线播放| 一级a爱视频在线免费观看| 国产主播在线观看一区二区| 色精品久久人妻99蜜桃| 久久伊人香网站| 国产精品一区二区精品视频观看| 国产精品免费一区二区三区在线| 国产色视频综合| 9色porny在线观看| 亚洲伊人色综图| 一个人免费在线观看的高清视频| 免费在线观看亚洲国产| 咕卡用的链子| 免费看a级黄色片| 嫩草影院精品99| 后天国语完整版免费观看| 欧美日韩乱码在线| 免费久久久久久久精品成人欧美视频| 又黄又爽又免费观看的视频| 男女做爰动态图高潮gif福利片 | aaaaa片日本免费| 色综合婷婷激情| 国产野战对白在线观看| av超薄肉色丝袜交足视频| 久久伊人香网站| 黄网站色视频无遮挡免费观看| 日本精品一区二区三区蜜桃| 国产熟女午夜一区二区三区| 国产激情欧美一区二区| 99国产极品粉嫩在线观看| 俄罗斯特黄特色一大片| 老司机福利观看| 一边摸一边抽搐一进一出视频| 日韩 欧美 亚洲 中文字幕| 亚洲黑人精品在线| 国产aⅴ精品一区二区三区波| 日韩中文字幕欧美一区二区| 久久国产精品男人的天堂亚洲| or卡值多少钱| 精品久久久久久成人av| 国产99白浆流出| 国产一区二区激情短视频| bbb黄色大片| 99国产精品99久久久久| 亚洲美女黄片视频| 一个人观看的视频www高清免费观看 | 亚洲欧美日韩无卡精品| 亚洲黑人精品在线| 夜夜躁狠狠躁天天躁| 黄色 视频免费看| 亚洲人成77777在线视频| 久久天躁狠狠躁夜夜2o2o| 国内久久婷婷六月综合欲色啪| 国产亚洲精品久久久久久毛片| 国产亚洲精品久久久久久毛片| 欧美黄色片欧美黄色片| 精品熟女少妇八av免费久了| 日本一区二区免费在线视频| 中文字幕另类日韩欧美亚洲嫩草| 久久这里只有精品19| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲专区国产一区二区| 18禁国产床啪视频网站| 日本vs欧美在线观看视频| 久久精品亚洲精品国产色婷小说| 久久婷婷成人综合色麻豆| 淫秽高清视频在线观看| 美女扒开内裤让男人捅视频| 国产成人啪精品午夜网站| 在线观看66精品国产| 深夜精品福利| 亚洲欧美激情综合另类| 成年版毛片免费区| 99国产极品粉嫩在线观看| 亚洲av熟女| 亚洲精华国产精华精| 欧美久久黑人一区二区| 91九色精品人成在线观看| 搡老岳熟女国产| 国产成人精品久久二区二区91| 女警被强在线播放| 亚洲色图综合在线观看| 久久人人精品亚洲av| 亚洲七黄色美女视频| 免费在线观看黄色视频的| 伦理电影免费视频| 手机成人av网站| 自线自在国产av| 天堂√8在线中文| 国产区一区二久久| 激情在线观看视频在线高清| 久久精品国产清高在天天线| 69精品国产乱码久久久| 在线观看免费视频日本深夜| 亚洲精品中文字幕一二三四区| 婷婷丁香在线五月| 脱女人内裤的视频| 欧美黄色片欧美黄色片| 熟女少妇亚洲综合色aaa.| 色哟哟哟哟哟哟| 国产精品自产拍在线观看55亚洲| 日日爽夜夜爽网站| 亚洲,欧美精品.| 午夜老司机福利片| 女人被躁到高潮嗷嗷叫费观| 日韩视频一区二区在线观看| 国产精品野战在线观看| 91麻豆精品激情在线观看国产| 国产精品98久久久久久宅男小说| 69av精品久久久久久| 精品国产国语对白av| 国产精品二区激情视频| 亚洲一区高清亚洲精品| 村上凉子中文字幕在线| 一级作爱视频免费观看| 国产精品一区二区精品视频观看| 狠狠狠狠99中文字幕| 91大片在线观看| 日韩欧美一区视频在线观看| 在线观看日韩欧美| 成人特级黄色片久久久久久久| 88av欧美| 精品一区二区三区av网在线观看| 国产免费av片在线观看野外av| 99精品久久久久人妻精品| 一区二区三区激情视频| 久久中文字幕人妻熟女| 麻豆一二三区av精品| 亚洲第一av免费看| 婷婷丁香在线五月| 精品福利观看| 亚洲黑人精品在线| 欧美绝顶高潮抽搐喷水| 啦啦啦 在线观看视频| 人人澡人人妻人| av天堂在线播放| 欧美日韩一级在线毛片| 757午夜福利合集在线观看| 俄罗斯特黄特色一大片| 久久久久久久久中文| 丁香欧美五月| 欧美一级毛片孕妇| 午夜视频精品福利| 午夜精品久久久久久毛片777| 中文亚洲av片在线观看爽| 日本免费a在线| 精品国产国语对白av| 男女做爰动态图高潮gif福利片 | 香蕉久久夜色| 亚洲性夜色夜夜综合| 国产片内射在线| www.自偷自拍.com| 亚洲人成电影观看| 亚洲国产高清在线一区二区三 | 91精品三级在线观看| 三级毛片av免费| 日本免费一区二区三区高清不卡 | 人人澡人人妻人| 亚洲第一青青草原| 免费不卡黄色视频| 午夜影院日韩av| 男人的好看免费观看在线视频 | 午夜久久久在线观看| 国产一级毛片七仙女欲春2 | 无人区码免费观看不卡| 女人被躁到高潮嗷嗷叫费观| 国产精品久久久久久人妻精品电影| 国产免费男女视频| 久久精品影院6| 啦啦啦观看免费观看视频高清 | 别揉我奶头~嗯~啊~动态视频| 亚洲三区欧美一区| 国产高清有码在线观看视频 | 他把我摸到了高潮在线观看| 亚洲免费av在线视频| 国产精品影院久久| 亚洲成av人片免费观看| 麻豆久久精品国产亚洲av| 国产精品久久视频播放| 成人国产综合亚洲| 母亲3免费完整高清在线观看| 日日夜夜操网爽| 国产一区二区在线av高清观看| 国产激情久久老熟女| 99re在线观看精品视频| 最近最新免费中文字幕在线| 精品不卡国产一区二区三区| 国产精品九九99| 亚洲色图综合在线观看| 黄色视频不卡| av天堂久久9| 丰满的人妻完整版| 可以免费在线观看a视频的电影网站| 无遮挡黄片免费观看| 人成视频在线观看免费观看| 精品少妇一区二区三区视频日本电影| 久久久久九九精品影院| 国产乱人伦免费视频| 99国产精品一区二区蜜桃av| 亚洲国产日韩欧美精品在线观看 | 多毛熟女@视频| 久久久久亚洲av毛片大全| 精品国产一区二区久久| 波多野结衣av一区二区av| tocl精华| 91字幕亚洲| 99久久精品国产亚洲精品| 妹子高潮喷水视频| 精品熟女少妇八av免费久了| 国产精品免费视频内射| 欧美成狂野欧美在线观看| 欧美亚洲日本最大视频资源| 在线观看午夜福利视频| 免费在线观看影片大全网站| cao死你这个sao货| 亚洲欧美激情在线| 叶爱在线成人免费视频播放| 久久 成人 亚洲| 久久久精品欧美日韩精品| 天堂影院成人在线观看| 长腿黑丝高跟| 老司机午夜福利在线观看视频| 国产精品香港三级国产av潘金莲| 亚洲第一电影网av| 久久狼人影院| 久久国产精品人妻蜜桃| 国产高清激情床上av| 91大片在线观看| 国产精品二区激情视频| 欧美久久黑人一区二区| 亚洲av日韩精品久久久久久密| 欧美久久黑人一区二区| 久久精品人人爽人人爽视色| 国产激情久久老熟女| 在线观看免费视频网站a站| svipshipincom国产片| 欧美成狂野欧美在线观看| 波多野结衣一区麻豆| 欧美不卡视频在线免费观看 | 精品卡一卡二卡四卡免费| 久久中文字幕人妻熟女| 美女国产高潮福利片在线看| 亚洲七黄色美女视频| 变态另类成人亚洲欧美熟女 | 美女国产高潮福利片在线看| 亚洲中文av在线| 国内精品久久久久精免费| 一a级毛片在线观看| 热99re8久久精品国产| 国产午夜福利久久久久久| 国产欧美日韩综合在线一区二区| 国产男靠女视频免费网站| 后天国语完整版免费观看| 久久欧美精品欧美久久欧美| 一本久久中文字幕| 亚洲全国av大片| 女人被狂操c到高潮| 韩国精品一区二区三区| 精品久久蜜臀av无| 99久久久亚洲精品蜜臀av| 日韩高清综合在线| 色婷婷久久久亚洲欧美| 青草久久国产| 韩国精品一区二区三区| 在线观看免费午夜福利视频| 久久国产乱子伦精品免费另类| 十八禁人妻一区二区| 一级黄色大片毛片| 母亲3免费完整高清在线观看| 欧美日韩瑟瑟在线播放| www.自偷自拍.com| 久99久视频精品免费| 亚洲欧美日韩无卡精品| 欧美大码av| 婷婷六月久久综合丁香| 日日爽夜夜爽网站| 欧美日韩乱码在线| 一本综合久久免费| 一个人免费在线观看的高清视频| 脱女人内裤的视频| 亚洲av第一区精品v没综合| 亚洲人成网站在线播放欧美日韩| 麻豆av在线久日| 国产欧美日韩一区二区三区在线| 国产成人欧美在线观看| 亚洲自拍偷在线| 伦理电影免费视频| 一区二区日韩欧美中文字幕| 成人欧美大片| 亚洲七黄色美女视频| 久久精品亚洲精品国产色婷小说| a在线观看视频网站| 日本免费一区二区三区高清不卡 | 色哟哟哟哟哟哟| 这个男人来自地球电影免费观看| 欧美+亚洲+日韩+国产| 欧美色欧美亚洲另类二区 | 亚洲免费av在线视频| 好男人电影高清在线观看| 欧美成人免费av一区二区三区| 1024香蕉在线观看| 久久午夜综合久久蜜桃| 国产精品自产拍在线观看55亚洲| 午夜激情av网站| 色播亚洲综合网| 每晚都被弄得嗷嗷叫到高潮| 一级片免费观看大全| 妹子高潮喷水视频| 国产成人欧美在线观看| 午夜福利18| 亚洲欧美精品综合久久99| 午夜福利18| 99在线视频只有这里精品首页| 成人特级黄色片久久久久久久| 男女之事视频高清在线观看| 男男h啪啪无遮挡| 精品高清国产在线一区| 国产精品1区2区在线观看.| 在线天堂中文资源库| 精品久久蜜臀av无| 亚洲国产精品sss在线观看| 欧美在线一区亚洲| 久久久久久久久免费视频了| 欧美黑人精品巨大| 很黄的视频免费| 伊人久久大香线蕉亚洲五| av超薄肉色丝袜交足视频| 757午夜福利合集在线观看| 一边摸一边做爽爽视频免费| 精品电影一区二区在线| 久久久久久大精品| 三级毛片av免费| 国产成+人综合+亚洲专区| 亚洲人成77777在线视频| 午夜福利成人在线免费观看| 美女国产高潮福利片在线看| 免费观看人在逋| 精品一区二区三区四区五区乱码| 桃红色精品国产亚洲av| 国产精品98久久久久久宅男小说| 国产亚洲欧美在线一区二区| 精品第一国产精品| 免费看a级黄色片| 亚洲人成77777在线视频| 成人国语在线视频| 夜夜看夜夜爽夜夜摸| 亚洲av第一区精品v没综合| 日韩欧美国产在线观看| 在线观看免费视频网站a站| 欧美日韩一级在线毛片| 欧美日韩黄片免| 亚洲av电影在线进入| 波多野结衣巨乳人妻| 国产在线观看jvid| 桃色一区二区三区在线观看| 亚洲欧美日韩另类电影网站| 午夜a级毛片| 桃色一区二区三区在线观看| 亚洲专区国产一区二区| 最近最新免费中文字幕在线| 亚洲欧美精品综合一区二区三区| 日韩一卡2卡3卡4卡2021年| 一进一出好大好爽视频| 啦啦啦观看免费观看视频高清 | 国内久久婷婷六月综合欲色啪| 欧美日韩中文字幕国产精品一区二区三区 | 久久国产精品人妻蜜桃| 亚洲久久久国产精品| 又黄又爽又免费观看的视频| 亚洲aⅴ乱码一区二区在线播放 | 精品久久久久久久久久免费视频| 国产免费av片在线观看野外av| 在线观看舔阴道视频| 香蕉国产在线看| 日韩免费av在线播放| 亚洲国产精品久久男人天堂| 黄色视频,在线免费观看| 搞女人的毛片| 后天国语完整版免费观看| 亚洲三区欧美一区| 91麻豆av在线| 成人免费观看视频高清| 亚洲av片天天在线观看| 黄片大片在线免费观看| 国产真人三级小视频在线观看| 亚洲avbb在线观看| 成在线人永久免费视频| 亚洲精品中文字幕在线视频| 国产精品秋霞免费鲁丝片| 精品乱码久久久久久99久播| 岛国在线观看网站| 这个男人来自地球电影免费观看| 亚洲国产毛片av蜜桃av| 国产99久久九九免费精品| 琪琪午夜伦伦电影理论片6080| 久久香蕉精品热| 久久中文看片网| 又黄又粗又硬又大视频| 国产熟女午夜一区二区三区| 亚洲 国产 在线| av天堂久久9| 两个人视频免费观看高清| 成人精品一区二区免费| a在线观看视频网站| 国产xxxxx性猛交| 大码成人一级视频| 亚洲精品国产精品久久久不卡| 亚洲欧美精品综合一区二区三区|