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

    碼頭群密度及非均勻度對近岸區(qū)洪水動力特性的影響研究

    2017-12-05 05:56:49
    水利水電快報 2017年11期
    關(guān)鍵詞:均勻度水流碼頭

    (1.長江勘測規(guī)劃設(shè)計研究有限責(zé)任公司 江河整治公司,湖北 武漢 430010;2. 鄂州市河道堤防管理處,湖北 鄂州 436000)

    碼頭群密度及非均勻度對近岸區(qū)洪水動力特性的影響研究

    江磊1蘆偉宏2胡春燕1侯衛(wèi)國1陳前海1陳正兵1

    (1.長江勘測規(guī)劃設(shè)計研究有限責(zé)任公司江河整治公司,湖北武漢430010;2.鄂州市河道堤防管理處,湖北鄂州436000)

    目前,碼頭群疊加對近岸水動力變化疊加影響的認識尚不清楚,研究碼頭群的疊加效應(yīng)對于指導(dǎo)工程實踐和優(yōu)化港區(qū)布局具有重要意義。采用資料分析和數(shù)值模擬相結(jié)合的方式,以南京河段新濟洲左汊左岸七壩港區(qū)為研究對象,分析了量化碼頭群數(shù)量和碼頭群排列方式的指標,選擇不同的量化指標模擬分析了碼頭群密度及碼頭群非均勻度對近岸洪水動力特性的疊加影響。結(jié)果表明,水位流速變化幅度及影響范圍與碼頭群密度正相關(guān),近岸水位壅高降低及流速增大減小的位置與碼頭群密度無關(guān),但幅度隨碼頭群密度的增大而增大;碼頭群非均勻度的影響與碼頭群密度不同,水位流速變化幅度及影響范圍隨非均勻度先增大后減小,近岸水位壅高降低與流速增大減小位置隨碼頭群非均勻度先增大后減小。碼頭群密度對水流特性的影響較碼頭群非均勻度大,在實踐過程中應(yīng)控制碼頭群的密度。

    碼頭群密度;非均勻度;洪水流量;疊加影響

    河道內(nèi)碼頭修建占用過水面積、增大局部阻力,導(dǎo)致局部水位壅高、流速增大,影響防洪安全并危及堤防穩(wěn)定。由于高樁碼頭對各類型的地基具有較強的適應(yīng)性,因而成為我國廣泛使用的碼頭結(jié)構(gòu)形式之一[1]。隨著國民經(jīng)濟發(fā)展,城市河道人類活動加劇,岸線開發(fā)利用程度增加,大量涉河工程的興建導(dǎo)致河道內(nèi)水位流速較天然河道更為多變,防洪形勢更加復(fù)雜[2]。城市局部河段內(nèi)由于水深條件好,通常碼頭連片布置,較短的岸線內(nèi)存在數(shù)個碼頭,形成碼頭群。碼頭群的存在可能會導(dǎo)致水流運動條件下碼頭間疊加或遮蔽效應(yīng)的產(chǎn)生,使其對水流運動特性的影響較單個碼頭更為明顯、復(fù)雜,也進一步加劇了防洪的安全隱患[3]。因此研究碼頭群的影響在豐富研究體系方面具有一定理論意義,而在保障防洪安全方面具有較大實際意義。

    大量研究者對單個碼頭或單塊樁群對水流結(jié)構(gòu)的影響進行了深入研究[4]。單個碼頭對水流結(jié)構(gòu)的影響主要表現(xiàn)為碼頭上游水位壅高、兩側(cè)流速增大,但對水位流速影響范圍較為有限[5-6]。單塊樁群對水流影響主要與樁群形狀、樁群密度、樁群寬度及長度有關(guān),樁群密度越大,樁群寬度則越寬;樁群長度越長,其對水流的影響越大[7-9]。目前少有研究涉及到碼頭群或樁群塊聯(lián)合體對水流的影響。碼頭群對水流的影響的研究主要涉及以下兩個問題:①如何量化碼頭群,應(yīng)采用何種指標對碼頭數(shù)量、排列方式進行定量描述?②不同量化指標條件下碼頭群對水流結(jié)構(gòu),尤其是對洪水水流結(jié)構(gòu)有何影響?

    本文以南京河段新濟洲左汊七壩港區(qū)為研究對象,采用平面二維水流數(shù)學(xué)模型研究了碼頭群對洪水水流結(jié)構(gòu)的影響,研究內(nèi)容主要包括:①碼頭群密度及排列方式量化指標初探;②碼頭群密度對洪水動力特性的影響;③碼頭群排列方式對洪水動力特性的影響;④碼頭群密度及非均勻度對洪水動力特性影響的比較分析。

    1 數(shù)值模型建立及研究工況

    1.1 基本控制方程

    本文通過二維水流數(shù)學(xué)模型模擬研究碼頭群對區(qū)域水動力特征的疊加影響,區(qū)域水動力特征變量(水位、流速、近岸水位流速等)通過求解水深平均的平面二維水流運動方程組得到??刂品匠探M的表達形式為

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

    (1)

    (2)水流動量方程:

    (2)

    (3)

    式中,u、v分別為垂線平均流速在x、y方向上的分量,m/s;h為總水深,h=η+d,d為靜止水深,η為水位,m;g為重力加速度,m/s2;τsx、τsy分別為表面切應(yīng)力分量,N/m2;τbx、τby分別為底部切應(yīng)力分量,N/m2;ρ為水的密度,ρ0為參考水密度,kg/m3;f為Coriolis參數(shù)(Ω為地球自傳角速度,Φ為地理緯度),f=2ΩsinΦ;fv和fu分別為地球自傳引起的加速度,m/s2;Sxx、Sxy、Syx、Syy為輻射應(yīng)力分量,N/m2;Txx、Txy、Tyx、Tyy為水平粘性應(yīng)力項,N/m2。

    Tij橫向上包括粘滯摩擦、動蕩摩擦和差別平流項。對這幾項的估測使用的是一種基于水深平均流速梯度的渦粘性公式:

    (4)

    式中,vs為渦粘系數(shù),m2/s。

    1.2 方程離散及求解

    采用三角網(wǎng)格對計算區(qū)域進行離散,對重點研究區(qū)域進行加密處理。在離散區(qū)域的三角形控制單元上,首先對方程組(1)~(3)進行積分,將方程離散為3部分,即時間積分項、空間積分項和源項積分項。時間積分項和源項積分項可直接利用顯式方法求解??臻g積分項的求解是關(guān)鍵,首先利用高斯公式將面積分轉(zhuǎn)化為線積分,線積分即為單元邊處的界面數(shù)值通量。界面數(shù)值通量有多種算法,本文模型采用求解法向一維Riemann問題的方法來計算界面數(shù)值通量。大量理論及實踐研究表明此方法在求解河道水流運動時具有較高的精度和穩(wěn)定性。

    1.3 模型率定與驗證

    本文所建模型范圍上起彭心洲洲頭、下至大勝關(guān),長度約為55 km。采用2011年8月實測的水位和流速分布驗證數(shù)學(xué)模型,驗證內(nèi)容包括新濟洲左右汊水面線及左右汊內(nèi)的流速分布。驗證結(jié)果表明,模擬計算水面線與實測水面線吻合度較高,模擬計算流速大小及分布形態(tài)與實測流速大小及分布形態(tài)較吻合,誤差滿足計算精度的要求,說明模型能較好地模擬新濟洲河段的水流運動。

    1.4 碼頭阻水特征分析

    防洪設(shè)計洪水條件下,高樁碼頭的阻水建筑物主要包括前沿平臺:橫梁、縱梁及樁基,平臺、橫梁、縱梁合稱為上部結(jié)構(gòu),樁基稱為下部結(jié)構(gòu),下部結(jié)構(gòu)是上部結(jié)構(gòu)的支撐建筑物。數(shù)學(xué)模型中以樁基為基本單元將碼頭概化為樁群,樁群中單個樁的下部為圓柱結(jié)構(gòu),上部為立方體結(jié)構(gòu)。大量統(tǒng)計分析結(jié)果表明碼頭各細部結(jié)構(gòu)占用過水面積的比率A平∶A橫∶A縱∶A樁為1∶0.379∶2.5∶8.533。由此可見,碼頭的主要阻水建筑物為下部樁基。模型中樁群的阻水效果主要通過增加局部網(wǎng)格的阻力實現(xiàn),阻力大小等于水流中障礙物的繞流阻力,其計算表達式為

    F=0.5ργCDAeV2

    (5)

    式中,F(xiàn)為單個樁在水流中的繞流阻力,N;γ為樁的形狀系數(shù),實際計算過程中根據(jù)不同的形狀取值;CD為拖拽力系數(shù);Ae為樁的有效阻水面積,m2;V為水流的平均流速,m/s。

    1.5 碼頭群量化特征分析

    港區(qū)內(nèi)碼頭群通常順岸布置、呈帶狀分布,因此碼頭群的特征主要體現(xiàn)在沿河道的縱向上。目前尚無較好的量化指標來綜合反映碼頭群的縱向特征。碼頭群縱向特征的量化主要涉及兩個問題:①如何在固定長度岸線內(nèi)量化碼頭數(shù)量;②固定碼頭數(shù)量條件下,如何量化分析碼頭排列方式差異。針對上述兩個問題,在分析長江南京河段左右岸碼頭分布特性的基礎(chǔ)上,本文提出了如下兩個指標。

    (1)碼頭群密度系數(shù)σ。碼頭群數(shù)量的多少對水流運動特性的影響可能存在明顯差異。碼頭群密度系數(shù)σ主要用來定量描述河段內(nèi)碼頭群數(shù)量,其定義為岸段內(nèi)碼頭平臺前沿線長度累加之和與岸線總長度的比率,表達式為

    (6)

    式中,li為第i個碼頭的平臺前沿線長度,m;L為岸線的總長度,m。由上式可知,密度系數(shù)σ克服了天然河道內(nèi)碼頭零碎、尺度大小不一,難以統(tǒng)計量化等缺點,是較好的天然河道碼頭建設(shè)密度量化指標。σ的取值范圍為[0,1],當(dāng)研究岸段內(nèi)無碼頭時,σ為0;當(dāng)研究岸段碼頭無間距全排列時,σ為1。

    (2) 碼頭群非均勻度系數(shù)ε。在河道某一岸段內(nèi)碼頭群數(shù)量不變的條件下,碼頭群順岸存在多種排列方式,不同排列方式對水流特性可能存在不同的影響。碼頭群非均勻度系數(shù)ε主要用來定量描述河段碼頭群排列方式的差異,其定義為岸段內(nèi)相鄰兩碼頭的最小間距dmin與最大間距dmax的比率。由于在dmin及dmax確定的條件下,dmin及dmax相對于水流流向存在兩種位置關(guān)系,即dmin在上游、dmax在下游和dmin在下游、dmax在上游。這就導(dǎo)致對于同一個ε值,依然存在兩種不同的碼頭群排列方式。因此附加定義當(dāng)dmin位于上游時,ε值取為正;dmin位于下游時,ε值取為負,具體表達式為

    (7)

    式中,dmin為岸段內(nèi)相鄰碼頭的最小間距,m;dmax為岸段內(nèi)相鄰碼頭的最大間距,m。由上式可知,ε的取值范圍為[-1,1]。當(dāng)岸段內(nèi)碼頭均勻排列時,ε為1;當(dāng)岸段內(nèi)相鄰碼頭搭界時,ε為0;當(dāng)岸段碼頭接近于均勻排列且dmin位于下游時,ε→-1。

    1.6 研究工況

    本文選取南京河段新濟洲左汊左岸七壩港區(qū)為研究對象,選定岸段岸線長度約為3 500 m。現(xiàn)狀條件下共有4個碼頭,碼頭密度系數(shù)σ為0.246,碼頭非均勻度系數(shù)ε為0.199。設(shè)定兩種典型情況分別將研究碼頭群密度、碼頭群非均勻程度對水流結(jié)構(gòu)的影響。典型情況1為碼頭群密度影響水流研究。在現(xiàn)狀碼頭密度條件下,分別將密度減小和增大1倍,形成3種工況,各工況碼頭均勻排列,非均勻度系數(shù)均為1.0(見表1)。典型情況2為碼頭群非均勻度影響水流研究。在均勻排列碼頭條件下,調(diào)整碼頭位置形成3種工況,各工況碼頭數(shù)量相同,密度系數(shù)均為0.212(見表2)。

    表1 碼頭密度影響水流的計算工況

    表2 碼頭非均勻系數(shù)影響水流的計算工況

    2 計算結(jié)果分析與討論

    2.1 碼頭群密度影響

    (1) 局部水位流速變化。碼頭興建后,由于占用了過水面積且阻力增大,水位變化主要發(fā)生在上游區(qū)域,以壅高為主;流速變化則較多發(fā)生在下游區(qū)域,主要表現(xiàn)為流速的減小。表3為局部區(qū)域內(nèi)水位最大壅高△Hmax及流速最大減小△vmax隨碼頭群密度σ的變化。由表3可知,△Hmax及△vmax隨σ的增大而增大,表明隨著碼頭群密度的增加,水位壅高及流速減小產(chǎn)生了明顯的疊加效應(yīng)。△Hmax及△vmax隨σ的增加總體呈線性增長趨勢,因此采用線性模型對數(shù)據(jù)進行擬合:△Hmax與σ的線性斜率系數(shù)k=11.64,回歸系數(shù)R2=0.97;△vmax與σ的線性斜率系數(shù)k=8.16,回歸系數(shù)R2=0.89。由此可證明碼頭群密度對水位壅高幅值的影響更明顯。

    表3 水位最大壅高及流速最大減小隨碼頭密度的變化

    圖1為局部區(qū)域水位壅高1.0 mm等值線隨碼頭群密度σ的變化圖,圖2為局部區(qū)域流速減小2.0 cm/s等值線隨碼頭群密度σ的變化圖。由圖1可知,水位壅高影響的區(qū)域大小S1及流速減小影響的區(qū)域大小S2亦隨碼頭群密度σ的增大而增大:當(dāng)σ=0.123時,S1=0.99×104m2;當(dāng)σ=0.246時,S1=21.87×104m2;當(dāng)σ=0.492時,S1=206.53×104m2。當(dāng)σ=0.123時,S2=15.89×104m2;當(dāng)σ=0.246時,S2=57.40×104m2;當(dāng)σ=0.492時,S2=62.35×104m2。由此可見,碼頭群密度對水位壅高范圍的影響較流速減小范圍的影響更明顯。

    圖1 不同密度條件下水位壅高1 mm等值線

    圖2 不同密度條件下流速減小2 cm/s等值線

    (2) 近岸水位流速變化。圖3為近岸水位變化△H隨碼頭群密度σ的變化圖,圖4為近岸流速△v變化隨碼頭群密度σ的變化圖(圖中x為自上游向下游的距離)。從圖3~4可知,近岸水位流速的變化幅度小于局部區(qū)域水位流速的變化幅度。碼頭密度對最大水位壅高(△Hmax)和降低(△Hmin)、最大流速增大(△vmax)和減小(△vmin)位置無明顯影響,但對其幅值影響顯著。當(dāng)σ=0.123時,△Hmax=0.55 mm,△Hmin=-0.75 mm;當(dāng)σ=0.246時,△Hmax=1.00 mm,△Hmin=-2.40 mm;當(dāng)σ=0.492時,△Hmax=1.80 mm,△Hmin=-4.65 mm。當(dāng)σ=0.123時,△vmax=0.80 cm/s,△vmin=-0.60 cm/s;當(dāng)σ=0.246時,△vmax=2.00 cm/s,△vmin=-4.00 cm/s;當(dāng)σ=0.492時,△vmax=4.50 cm/s,△vmin=-8.00 cm/s。由此可見,隨著碼頭密度增大,近岸水位流速基本呈線性增長趨勢。

    圖3 不同碼頭密度條件下近岸水位變化

    圖4 不同碼頭密度條件下近岸流速變化

    2.2 碼頭群非均勻度影響

    (1) 水位流速變化。碼頭群非均勻度變化主要改變碼頭的排列方式,非均勻度值絕對值越小,碼頭群排列越不規(guī)則。表4為局部區(qū)域內(nèi)水位最大壅高△Hmax及流速最大減小△vmax隨碼頭群非均勻度ε的變化。從表4可知,△Hmax及△vmax隨ε的增大呈先增大后減小的變化趨勢,說明水位壅高及流速減小與碼頭群非均勻度并非是單一的關(guān)系?!鱄max及△vmax隨ε增加無明顯的變化趨勢,采用線性模型對數(shù)據(jù)進行擬合:△Hmax與σ的線性斜率系數(shù)k=0.71,回歸系數(shù)R2=0.35;△vmax與ε的線性斜率系數(shù)k=0.31,回歸系數(shù)R2=0.12;由此可見,碼頭群非均勻度對水位壅高及流速減小影響不明顯。

    圖5為局部區(qū)域水位壅高1.0 mm等值線隨碼頭群非均勻度ε的變化,圖6為局部區(qū)域流速減小2.0 cm/s等值線隨碼頭群非均勻度ε的變化圖。從圖5~6可知,水位壅高影響的區(qū)域大小S1隨碼頭群非均勻度ε的增大呈先增大后減小的變化趨勢,當(dāng)ε=-0.334時,S1=3.06×104m2;當(dāng)ε=0.334時,S1=17.15×104m2;當(dāng)ε=1.000時,S1=3.22×104m2。流速減小影響的區(qū)域大小S2隨碼頭群非均勻度ε的增大基本保持不變,當(dāng)ε=-0.334時,S2=43.7×104m2;當(dāng)ε=0.334時,S2=47.40×104m2;當(dāng)ε=1.000時,S2=42.35×104m2。由此可見,碼頭非均勻度對水位壅高范圍的影響較流速減小范圍的影響更明顯。

    表4 水位最大壅高及流速最大減小隨碼頭非均勻度的變化

    圖5 不同非均勻度條件下水位壅高1 mm等值線

    圖7 不同非均勻度條件下近岸水位變化

    圖8 不同非均勻度條件下近岸流速變化

    (2) 近岸水位流速變化。圖7為近岸水位變化△H隨碼頭群非均勻度ε的變化圖,圖8為近岸流速△v變化隨碼頭群非均勻度ε的變化圖(圖中x為自上游向下游的距離)。從圖7~8可知,近岸水位流速的變化幅度小于局部區(qū)域流速的變化幅度。近岸水位的變化趨勢依然是上游壅高下游減小,最大水位壅高(△Hmax)隨ε的增加呈現(xiàn)先增大后減小的趨勢,最大水位降低(△Hmin)隨ε的增加呈現(xiàn)先減小后增大的趨勢。當(dāng)ε=-0.334時,△Hmax=0.70 mm,△Hmin=-2.20 mm;當(dāng)ε=0.334時,△Hmax=1.10 mm,△Hmin=-1.90 mm;當(dāng)ε=1.000時,△Hmax=0.80 mm,△Hmin=-2.1 mm。近岸流速的變化較為復(fù)雜,ε既影響流速變化的位置也影響流速變化的幅度。最大流速增大(△vmax)隨ε的增加先減小后增大,最大流速減小(△vmin)隨ε基本保持不變位置。當(dāng)ε=-0.334時,△vmax=0.60 cm/s,△vmin=-3.20 cm/s;當(dāng)ε=0.334時,△vmax=1.40 cm/s,△vmin=-3.15 cm/s;當(dāng)ε=1.000時,△vmax=0.60 cm/s,△vmin=-3.20 cm/s。由此可見,碼頭非均勻度對近岸水位的影響較流速的影響更明顯。

    相鄰碼頭最小間距dmin出現(xiàn)的位置可能是碼頭非均勻度影響呈先增大后減小的重要原因,對于ε=-0.334和ε=0.334兩種工況,盡管非均勻度絕對值大小相同,但兩者的差別主要體現(xiàn)在:當(dāng)ε=0.334時,dmin位于上游;當(dāng)ε=-0.334時,dmin位于下游,上游碼頭群可能對下游碼頭群的影響產(chǎn)生了明顯的屏蔽效果。為驗證此設(shè)想,將碼頭群密度影響和碼頭群非均勻度影響的計算結(jié)果進行統(tǒng)一分析,統(tǒng)計各工況上游區(qū)域碼頭最小間距dmin,由于工況ε=-0.334條件下dmin位于下游區(qū)域,因此用最大間距dmax代替上游區(qū)域的最小間距dmin。對上游區(qū)相鄰碼頭最小間距dmin與局部最大水位壅高△Hmax的關(guān)系進行統(tǒng)計,并畫圖進行回歸分析(見圖9)。從圖可知,隨dmin增大,△Hmax逐漸減小,兩者呈現(xiàn)明顯的負相關(guān)關(guān)系,采用冪指數(shù)函數(shù)進行回歸,回歸系數(shù)R2=0.815 1。這也進一步說明碼頭群非均勻度對水流的影響主要與上游區(qū)相鄰碼頭的最小間距有關(guān)。

    圖9 上游區(qū)最小碼頭間距與最大水位壅高關(guān)系

    3 結(jié) 論

    (1) 洪水動力特性變化與碼頭群密度正相關(guān),局部水位最大壅高及影響范圍、局部流速最大減小及影響范圍、近岸水位流速增大減小均隨碼頭密度的增大而增大。

    (2) 當(dāng)相鄰碼頭最小間距位于上游時,洪水動力特性變化與碼頭非均勻度負相關(guān),局部水位最大壅高及影響范圍、局部流速最大減小及影響范圍、近岸水位流速增大減小均隨碼頭群非均勻度的增大而增大。

    (3) 當(dāng)相鄰碼頭最小間距位于下游時,局部水位最大壅高及影響范圍、局部流速最大減小及影響范圍、近岸水位流速增大減小均隨碼頭群非均勻度的增大保持不變。

    (4) 碼頭群密度對洪水動力特性的影響較碼頭群非均勻度大,在實踐過程中應(yīng)更注意控制岸線內(nèi)的碼頭密度。

    [1] 黃本勝,程香菊,袁麗蓉,等. 碼頭樁群對河道行洪與流場影響的三維數(shù)值模擬[J]. 水動力學(xué)研究與進展A輯,2010,25(1):99-105.

    [2] 吳飛. 高樁碼頭樁群對河道水流特性影響的數(shù)值模擬研究[J]. 水資源與水工程研究,2008年,2008,19(1):89-93.

    [3] 李彬,劉博,鄭國棟. 碼頭群連體擴建對河道水流影響的數(shù)值模擬[J]. 廣東水利水電,2016(6):1-3.

    [4] 李文文,黃本勝,候杰. 高樁碼頭樁群水流特性的試驗研究[J],新疆農(nóng)業(yè)大學(xué)學(xué)報(自然科學(xué)版),2004,27(4):78-81.

    [5] 王春華. 高樁碼頭樁群對流場及底床影響試驗研究[D]. 南京:南京水利科學(xué)研究院, 2010:30-35.

    [6] 李光熾,周晶晏,張貴壽. 高樁碼頭對河道流場影響的數(shù)值模擬[J]. 河海大學(xué)學(xué)報(自然科學(xué)版),2004,32(2):216-220.

    [7] Stephen E. Coleman. Clearwater local scour at complex piers[J],Journal of Hydraulic Engineering,2005,131(4):330-334.

    [8] Siniscalchi F, Nikora V I, Aberle J. Plant patch hydrodynamics in streams: mean flow, turbulence, and drag forces.[J] Water Resources Research, 2012, 48 (1): 1513-1519.

    [9] Rominger J, Nepf H M. Flow adjustment and interior flow associated with a rectangular porous obstruction.[J] Journal of Fluid Mechanics, 2011, 680: 636-659.

    (編輯:李曉濛)

    2017-09-15

    國家重點研發(fā)計劃課題(2016YFC0402306)

    江 磊,男,長江勘測規(guī)劃設(shè)計研究有限責(zé)任公司江河整治公司,工程師.

    1006-0081(2017)11-0047-06

    TV131.2

    A

    猜你喜歡
    均勻度水流碼頭
    低播量下雜交稻產(chǎn)量形成對種植均勻度的響應(yīng)
    作物研究(2023年2期)2023-05-28 13:44:14
    全自動化碼頭來了
    哪股水流噴得更遠
    能俘獲光的水流
    我只知身在水中,不覺水流
    文苑(2020年6期)2020-06-22 08:41:56
    均勻度控制不佳可致肉種雞晚產(chǎn)
    前往碼頭
    在碼頭上釣魚
    錦綸長絲染色均勻度判色新方法
    復(fù)方丹參片中冰片的含量均勻度研究
    中成藥(2014年10期)2014-02-28 22:29:24
    色播在线永久视频| 亚洲精华国产精华精| 国内毛片毛片毛片毛片毛片| 男女免费视频国产| 国产精品久久视频播放| 在线观看免费午夜福利视频| 色婷婷av一区二区三区视频| 可以免费在线观看a视频的电影网站| 国产97色在线日韩免费| 自拍欧美九色日韩亚洲蝌蚪91| 99久久综合精品五月天人人| 99久久综合精品五月天人人| 久久ye,这里只有精品| 国产精品一区二区在线不卡| 精品国内亚洲2022精品成人 | 成年版毛片免费区| 九色亚洲精品在线播放| 欧美亚洲日本最大视频资源| 免费看a级黄色片| 另类亚洲欧美激情| 天堂中文最新版在线下载| 国产国语露脸激情在线看| 久久午夜综合久久蜜桃| 欧美人与性动交α欧美软件| 热re99久久精品国产66热6| 精品一品国产午夜福利视频| av网站免费在线观看视频| 少妇的丰满在线观看| 精品国内亚洲2022精品成人 | 69av精品久久久久久| 少妇的丰满在线观看| 美国免费a级毛片| 悠悠久久av| 巨乳人妻的诱惑在线观看| 午夜福利视频在线观看免费| 亚洲五月色婷婷综合| 最新在线观看一区二区三区| 国产av一区二区精品久久| 亚洲专区国产一区二区| 亚洲在线自拍视频| 中文字幕另类日韩欧美亚洲嫩草| 777久久人妻少妇嫩草av网站| 丁香六月欧美| 中文字幕色久视频| 亚洲人成77777在线视频| 亚洲av成人av| 久久久久久久国产电影| 看黄色毛片网站| 国产男女内射视频| 免费在线观看黄色视频的| 两性夫妻黄色片| 麻豆国产av国片精品| 久久国产精品男人的天堂亚洲| 高清欧美精品videossex| 成年女人毛片免费观看观看9 | 免费久久久久久久精品成人欧美视频| 久久久久久久午夜电影 | 激情在线观看视频在线高清 | av中文乱码字幕在线| 欧美国产精品va在线观看不卡| 黄色怎么调成土黄色| 欧美在线一区亚洲| av在线播放免费不卡| 在线十欧美十亚洲十日本专区| 欧美精品av麻豆av| av天堂久久9| 777米奇影视久久| 欧美日韩亚洲国产一区二区在线观看 | 欧美精品亚洲一区二区| 日本vs欧美在线观看视频| 亚洲av成人av| 十八禁高潮呻吟视频| 女性被躁到高潮视频| 国产精品一区二区在线不卡| 女警被强在线播放| 岛国毛片在线播放| 国产欧美日韩一区二区三| 午夜91福利影院| 亚洲精品一卡2卡三卡4卡5卡| 国产aⅴ精品一区二区三区波| 久久国产精品大桥未久av| 不卡一级毛片| 久久精品亚洲熟妇少妇任你| 狂野欧美激情性xxxx| 国产精品成人在线| 久久久久久久久久久久大奶| 18在线观看网站| 色精品久久人妻99蜜桃| av在线播放免费不卡| 国产精品二区激情视频| 亚洲第一青青草原| 午夜福利,免费看| 国产精品香港三级国产av潘金莲| 一区福利在线观看| 久久国产精品大桥未久av| 国产av又大| 欧美国产精品一级二级三级| 美国免费a级毛片| 女人被躁到高潮嗷嗷叫费观| 久久狼人影院| 一二三四社区在线视频社区8| 亚洲精品自拍成人| 亚洲av欧美aⅴ国产| 国产精品秋霞免费鲁丝片| 免费在线观看亚洲国产| 每晚都被弄得嗷嗷叫到高潮| 午夜两性在线视频| 另类亚洲欧美激情| 丝袜人妻中文字幕| e午夜精品久久久久久久| 看黄色毛片网站| 美女视频免费永久观看网站| 80岁老熟妇乱子伦牲交| 午夜亚洲福利在线播放| 国产成人免费观看mmmm| 成熟少妇高潮喷水视频| 亚洲精品av麻豆狂野| 欧美乱码精品一区二区三区| 亚洲五月色婷婷综合| 黄色女人牲交| 久久99一区二区三区| 操美女的视频在线观看| 亚洲精品久久午夜乱码| 又紧又爽又黄一区二区| av视频免费观看在线观看| 精品国产一区二区三区久久久樱花| 18禁裸乳无遮挡免费网站照片 | 女同久久另类99精品国产91| 涩涩av久久男人的天堂| 一本一本久久a久久精品综合妖精| 亚洲国产精品一区二区三区在线| 在线观看免费午夜福利视频| 丰满迷人的少妇在线观看| 一边摸一边抽搐一进一出视频| 日本五十路高清| 午夜福利免费观看在线| 久久香蕉激情| 夜夜躁狠狠躁天天躁| 亚洲欧美精品综合一区二区三区| 午夜福利欧美成人| 国产男女超爽视频在线观看| 精品福利永久在线观看| 老汉色av国产亚洲站长工具| 母亲3免费完整高清在线观看| 精品一区二区三区四区五区乱码| 高潮久久久久久久久久久不卡| 国产精品国产高清国产av | 欧美亚洲 丝袜 人妻 在线| 精品午夜福利视频在线观看一区| 国产成人免费无遮挡视频| 在线观看免费视频日本深夜| 美女午夜性视频免费| 欧美日韩亚洲综合一区二区三区_| 一区在线观看完整版| 国产日韩一区二区三区精品不卡| 一二三四在线观看免费中文在| 少妇被粗大的猛进出69影院| 亚洲一码二码三码区别大吗| 久久久久视频综合| 看免费av毛片| 麻豆成人av在线观看| 午夜福利欧美成人| 中文字幕另类日韩欧美亚洲嫩草| 精品熟女少妇八av免费久了| 国产精品 欧美亚洲| 在线国产一区二区在线| 国产亚洲精品久久久久久毛片 | 99热网站在线观看| 91老司机精品| 亚洲色图av天堂| 久久国产精品影院| 亚洲美女黄片视频| 深夜精品福利| 久久久久久久精品吃奶| √禁漫天堂资源中文www| 国产亚洲精品一区二区www | 久久香蕉国产精品| av电影中文网址| 水蜜桃什么品种好| 十八禁高潮呻吟视频| 视频区欧美日本亚洲| 中国美女看黄片| 欧美另类亚洲清纯唯美| 高清欧美精品videossex| 99久久人妻综合| 99在线人妻在线中文字幕 | 欧美国产精品va在线观看不卡| 国产精品影院久久| 亚洲av日韩精品久久久久久密| 久久性视频一级片| 欧美精品av麻豆av| 水蜜桃什么品种好| 搡老岳熟女国产| 国产亚洲精品久久久久5区| aaaaa片日本免费| 一级毛片高清免费大全| 一进一出抽搐gif免费好疼 | 人妻一区二区av| 午夜福利在线观看吧| 99久久综合精品五月天人人| 成人永久免费在线观看视频| av中文乱码字幕在线| 亚洲中文av在线| 国产熟女午夜一区二区三区| 国产成人av激情在线播放| 黄网站色视频无遮挡免费观看| 丰满的人妻完整版| a级毛片在线看网站| 丰满人妻熟妇乱又伦精品不卡| 老司机亚洲免费影院| 国产成人系列免费观看| 俄罗斯特黄特色一大片| 两人在一起打扑克的视频| 一本综合久久免费| tube8黄色片| 狠狠狠狠99中文字幕| 一级片'在线观看视频| 欧美不卡视频在线免费观看 | 午夜视频精品福利| 99在线人妻在线中文字幕 | 黑人巨大精品欧美一区二区蜜桃| 五月开心婷婷网| 人人妻人人澡人人看| 亚洲国产欧美网| 一进一出抽搐gif免费好疼 | 中国美女看黄片| 午夜久久久在线观看| 一边摸一边抽搐一进一出视频| 一区二区日韩欧美中文字幕| 久久久久国内视频| 免费av中文字幕在线| 国产亚洲av高清不卡| 人人妻人人爽人人添夜夜欢视频| 一区在线观看完整版| 99国产精品99久久久久| 91成人精品电影| 性色av乱码一区二区三区2| av网站在线播放免费| av免费在线观看网站| 亚洲第一青青草原| 久久这里只有精品19| 女人爽到高潮嗷嗷叫在线视频| 国产亚洲精品第一综合不卡| 一夜夜www| 夜夜夜夜夜久久久久| 一本大道久久a久久精品| 啦啦啦免费观看视频1| av超薄肉色丝袜交足视频| 多毛熟女@视频| 亚洲欧美精品综合一区二区三区| 18禁国产床啪视频网站| 中文欧美无线码| 黄色怎么调成土黄色| 伊人久久大香线蕉亚洲五| 香蕉丝袜av| 老司机午夜福利在线观看视频| 国产精品乱码一区二三区的特点 | 无限看片的www在线观看| 一a级毛片在线观看| 大香蕉久久成人网| 欧美大码av| 9191精品国产免费久久| av线在线观看网站| 90打野战视频偷拍视频| 精品熟女少妇八av免费久了| 成人精品一区二区免费| 乱人伦中国视频| 国产有黄有色有爽视频| 欧美一级毛片孕妇| 国产乱人伦免费视频| 热99久久久久精品小说推荐| 中文欧美无线码| 啦啦啦在线免费观看视频4| 久久人人97超碰香蕉20202| 黑丝袜美女国产一区| 国产精品亚洲av一区麻豆| 老熟妇乱子伦视频在线观看| 午夜福利一区二区在线看| 男人的好看免费观看在线视频 | 久久久久国产一级毛片高清牌| 亚洲专区国产一区二区| 人成视频在线观看免费观看| 国产精品乱码一区二三区的特点 | 日韩欧美三级三区| 国产片内射在线| 国产一区有黄有色的免费视频| 久久香蕉激情| 亚洲三区欧美一区| 国产成人精品久久二区二区91| 亚洲av片天天在线观看| 日韩熟女老妇一区二区性免费视频| 最近最新中文字幕大全电影3 | 国产极品粉嫩免费观看在线| 在线观看免费高清a一片| 午夜精品国产一区二区电影| 欧美黄色片欧美黄色片| av网站免费在线观看视频| www.熟女人妻精品国产| 亚洲精品国产区一区二| 日本欧美视频一区| 啪啪无遮挡十八禁网站| xxx96com| 又紧又爽又黄一区二区| 亚洲国产毛片av蜜桃av| 大型av网站在线播放| 很黄的视频免费| 免费在线观看视频国产中文字幕亚洲| 人妻 亚洲 视频| 亚洲一卡2卡3卡4卡5卡精品中文| 啦啦啦在线免费观看视频4| 老汉色∧v一级毛片| 天天添夜夜摸| 免费av中文字幕在线| 男人舔女人的私密视频| 欧美日韩福利视频一区二区| 麻豆乱淫一区二区| 国产黄色免费在线视频| 亚洲欧美日韩另类电影网站| 久久亚洲真实| 欧美 亚洲 国产 日韩一| 五月开心婷婷网| 黄色a级毛片大全视频| av欧美777| 女人精品久久久久毛片| 亚洲美女黄片视频| 亚洲成国产人片在线观看| 久久精品国产a三级三级三级| 久久久国产成人免费| 女同久久另类99精品国产91| 人成视频在线观看免费观看| 无限看片的www在线观看| 成人亚洲精品一区在线观看| 中文字幕高清在线视频| 久久精品人人爽人人爽视色| 看黄色毛片网站| 国产欧美亚洲国产| 两个人看的免费小视频| 亚洲男人天堂网一区| av一本久久久久| 久久国产精品男人的天堂亚洲| 久久久久久久久久久久大奶| 九色亚洲精品在线播放| 视频在线观看一区二区三区| 亚洲成人手机| 日韩视频一区二区在线观看| 亚洲欧美一区二区三区久久| 大陆偷拍与自拍| 一级毛片精品| 成人亚洲精品一区在线观看| 亚洲av片天天在线观看| 欧美日韩亚洲综合一区二区三区_| 日日爽夜夜爽网站| 捣出白浆h1v1| 啦啦啦视频在线资源免费观看| 亚洲第一青青草原| www.精华液| 国产不卡一卡二| 国产亚洲欧美98| 亚洲精品在线观看二区| 视频区图区小说| 老熟妇仑乱视频hdxx| 一区在线观看完整版| 精品国产亚洲在线| 久9热在线精品视频| 在线观看www视频免费| a级毛片在线看网站| 女性被躁到高潮视频| 精品国产乱码久久久久久男人| 亚洲 欧美一区二区三区| 午夜福利欧美成人| 久久香蕉激情| 亚洲五月色婷婷综合| 亚洲男人天堂网一区| 热99re8久久精品国产| 麻豆成人av在线观看| 国产成人精品无人区| 亚洲欧美日韩高清在线视频| 欧美在线黄色| 中文字幕人妻熟女乱码| 色老头精品视频在线观看| 99精国产麻豆久久婷婷| 国产成人av激情在线播放| 亚洲国产毛片av蜜桃av| 欧美色视频一区免费| 不卡av一区二区三区| 亚洲精品美女久久av网站| 女人被狂操c到高潮| 少妇粗大呻吟视频| 亚洲av欧美aⅴ国产| 成人18禁高潮啪啪吃奶动态图| 极品少妇高潮喷水抽搐| 桃红色精品国产亚洲av| 精品久久久精品久久久| 久久久精品免费免费高清| 最新在线观看一区二区三区| 妹子高潮喷水视频| 亚洲aⅴ乱码一区二区在线播放 | 日韩三级视频一区二区三区| www.熟女人妻精品国产| 亚洲成人免费电影在线观看| 国产主播在线观看一区二区| 午夜老司机福利片| 亚洲va日本ⅴa欧美va伊人久久| 日本五十路高清| 久久精品亚洲av国产电影网| 欧美国产精品va在线观看不卡| 国产欧美亚洲国产| 大片电影免费在线观看免费| 国产亚洲一区二区精品| 大香蕉久久网| 国产成人欧美在线观看 | 亚洲五月天丁香| 国精品久久久久久国模美| 国产精品久久久av美女十八| 亚洲九九香蕉| 高潮久久久久久久久久久不卡| 亚洲少妇的诱惑av| 欧美成狂野欧美在线观看| 国产成+人综合+亚洲专区| 人妻一区二区av| 成人三级做爰电影| 日韩欧美国产一区二区入口| 亚洲aⅴ乱码一区二区在线播放 | 18禁裸乳无遮挡免费网站照片 | 精品国产一区二区三区四区第35| 大码成人一级视频| 久久亚洲精品不卡| 淫妇啪啪啪对白视频| 国产午夜精品久久久久久| 啦啦啦视频在线资源免费观看| 亚洲av电影在线进入| 欧美日韩国产mv在线观看视频| 国产有黄有色有爽视频| 99国产精品一区二区三区| 国产精品久久久久久人妻精品电影| 97人妻天天添夜夜摸| 国产成人影院久久av| 妹子高潮喷水视频| 99香蕉大伊视频| 亚洲精品中文字幕在线视频| 三级毛片av免费| 精品久久久久久,| 精品久久蜜臀av无| 在线永久观看黄色视频| 91大片在线观看| 一二三四社区在线视频社区8| 男男h啪啪无遮挡| 国产精品av久久久久免费| 免费在线观看视频国产中文字幕亚洲| 国内久久婷婷六月综合欲色啪| 亚洲专区中文字幕在线| 国产精品美女特级片免费视频播放器 | 操美女的视频在线观看| 99re在线观看精品视频| 欧美大码av| 亚洲精品国产精品久久久不卡| 一级毛片精品| 色94色欧美一区二区| av不卡在线播放| 亚洲av熟女| 国产乱人伦免费视频| 亚洲av熟女| 中文亚洲av片在线观看爽 | 精品亚洲成国产av| www.精华液| 在线免费观看的www视频| av国产精品久久久久影院| 午夜福利免费观看在线| av不卡在线播放| 最新的欧美精品一区二区| 999精品在线视频| 日本撒尿小便嘘嘘汇集6| 波多野结衣av一区二区av| 国产区一区二久久| tube8黄色片| 久久久久久亚洲精品国产蜜桃av| 久久精品aⅴ一区二区三区四区| 两个人看的免费小视频| 日韩熟女老妇一区二区性免费视频| 99国产精品99久久久久| 男女免费视频国产| 在线观看免费视频日本深夜| 午夜老司机福利片| 村上凉子中文字幕在线| svipshipincom国产片| 欧美黄色片欧美黄色片| 女人爽到高潮嗷嗷叫在线视频| 欧美在线一区亚洲| 老司机亚洲免费影院| 亚洲五月天丁香| 欧美日韩亚洲高清精品| 亚洲人成伊人成综合网2020| 欧美在线黄色| 黑人欧美特级aaaaaa片| 丝袜美腿诱惑在线| 亚洲欧美精品综合一区二区三区| 日韩有码中文字幕| 国产xxxxx性猛交| 高清在线国产一区| 人人妻人人爽人人添夜夜欢视频| 啦啦啦视频在线资源免费观看| 日本a在线网址| 无限看片的www在线观看| 最新的欧美精品一区二区| 午夜精品国产一区二区电影| 人妻久久中文字幕网| 99热只有精品国产| 一区二区三区国产精品乱码| 亚洲aⅴ乱码一区二区在线播放 | 1024香蕉在线观看| 首页视频小说图片口味搜索| 久久精品国产清高在天天线| 我的亚洲天堂| 国产精品av久久久久免费| 变态另类成人亚洲欧美熟女 | 日本a在线网址| av免费在线观看网站| 免费在线观看黄色视频的| 99re在线观看精品视频| 黑人欧美特级aaaaaa片| 免费观看人在逋| 9热在线视频观看99| 欧美精品一区二区免费开放| 精品亚洲成a人片在线观看| 成人影院久久| 国产精品av久久久久免费| 国产精品98久久久久久宅男小说| 黄色片一级片一级黄色片| 国产又爽黄色视频| 黄色怎么调成土黄色| 日韩欧美免费精品| 亚洲欧美激情在线| 免费观看a级毛片全部| 丁香欧美五月| 国产在线精品亚洲第一网站| 黄频高清免费视频| 女同久久另类99精品国产91| 99热网站在线观看| 国产不卡一卡二| 国产激情欧美一区二区| 天堂中文最新版在线下载| 欧美性长视频在线观看| av超薄肉色丝袜交足视频| 美女国产高潮福利片在线看| 视频在线观看一区二区三区| 国产亚洲一区二区精品| 国产欧美日韩综合在线一区二区| 好男人电影高清在线观看| 操美女的视频在线观看| 国产日韩欧美亚洲二区| 精品欧美一区二区三区在线| 国产又色又爽无遮挡免费看| 日韩欧美在线二视频 | 国产精品欧美亚洲77777| 少妇被粗大的猛进出69影院| 久久人妻av系列| 中文字幕人妻丝袜一区二区| xxxhd国产人妻xxx| 19禁男女啪啪无遮挡网站| 最新美女视频免费是黄的| 两个人免费观看高清视频| 久久人妻av系列| 18禁国产床啪视频网站| 亚洲欧洲精品一区二区精品久久久| 欧美黑人精品巨大| 国产成人啪精品午夜网站| 久久狼人影院| www.熟女人妻精品国产| av视频免费观看在线观看| 免费av中文字幕在线| 一进一出好大好爽视频| 色综合欧美亚洲国产小说| 久久精品熟女亚洲av麻豆精品| 校园春色视频在线观看| 亚洲片人在线观看| 亚洲精品乱久久久久久| 欧美av亚洲av综合av国产av| 大片电影免费在线观看免费| 国产淫语在线视频| 黄色丝袜av网址大全| 波多野结衣av一区二区av| 一级片'在线观看视频| 91成人精品电影| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲片人在线观看| 中出人妻视频一区二区| 黄色女人牲交| 久久久久久久久免费视频了| 一二三四社区在线视频社区8| 黄色怎么调成土黄色| 国产男女超爽视频在线观看| 国产精品一区二区精品视频观看| 午夜精品在线福利| 国产成人欧美| 国产男女内射视频| 国产午夜精品久久久久久| 女性被躁到高潮视频| 久久这里只有精品19| 无遮挡黄片免费观看| 亚洲 欧美一区二区三区| 大型av网站在线播放| 亚洲精品在线观看二区| 欧美大码av| 久久99一区二区三区| 欧美+亚洲+日韩+国产| 丰满迷人的少妇在线观看| 成人免费观看视频高清| 一区二区三区激情视频| 亚洲成国产人片在线观看| 亚洲中文av在线| 一级a爱视频在线免费观看| 十分钟在线观看高清视频www| 日韩三级视频一区二区三区| 国产精品久久久久成人av| 国产成人av教育|