• <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
    9191精品国产免费久久| 在线看a的网站| 亚洲中文日韩欧美视频| 日本91视频免费播放| 汤姆久久久久久久影院中文字幕| 蜜桃在线观看..| 国产成人免费观看mmmm| 国产精品香港三级国产av潘金莲 | 欧美日韩视频高清一区二区三区二| 一区福利在线观看| 丰满迷人的少妇在线观看| 性色av乱码一区二区三区2| 美女中出高潮动态图| 曰老女人黄片| av天堂久久9| 国产精品偷伦视频观看了| 日韩,欧美,国产一区二区三区| 激情视频va一区二区三区| 久9热在线精品视频| 国产在线观看jvid| www日本在线高清视频| 桃花免费在线播放| 久久久久久久大尺度免费视频| 日本猛色少妇xxxxx猛交久久| 国产国语露脸激情在线看| 超色免费av| 欧美日本中文国产一区发布| 亚洲少妇的诱惑av| 伊人久久大香线蕉亚洲五| 欧美日韩黄片免| 91精品伊人久久大香线蕉| 自线自在国产av| 欧美久久黑人一区二区| 777米奇影视久久| 女性生殖器流出的白浆| 欧美日韩国产mv在线观看视频| 51午夜福利影视在线观看| 成年动漫av网址| 少妇粗大呻吟视频| 亚洲欧美激情在线| 国产亚洲欧美精品永久| 精品国产乱码久久久久久小说| 色精品久久人妻99蜜桃| 国产精品一区二区在线观看99| 男女午夜视频在线观看| 亚洲欧美激情在线| xxxhd国产人妻xxx| 69精品国产乱码久久久| 男女之事视频高清在线观看 | 国产视频一区二区在线看| 69精品国产乱码久久久| 成人亚洲欧美一区二区av| 久热这里只有精品99| 久久这里只有精品19| 99国产精品一区二区蜜桃av | 欧美变态另类bdsm刘玥| 日韩制服丝袜自拍偷拍| 国产麻豆69| 人成视频在线观看免费观看| 久久人人爽人人片av| 亚洲九九香蕉| 黄色一级大片看看| 天天躁夜夜躁狠狠躁躁| 在现免费观看毛片| 免费在线观看日本一区| 成人国语在线视频| 99国产综合亚洲精品| 亚洲精品一二三| 亚洲精品一区蜜桃| 日日爽夜夜爽网站| 色综合欧美亚洲国产小说| 丝袜人妻中文字幕| 国产午夜精品一二区理论片| 国产亚洲av高清不卡| svipshipincom国产片| 黑人猛操日本美女一级片| 久久久久视频综合| 国产精品国产三级专区第一集| a级片在线免费高清观看视频| 麻豆av在线久日| 国产伦人伦偷精品视频| 80岁老熟妇乱子伦牲交| 99九九在线精品视频| 一级毛片 在线播放| 制服诱惑二区| 国产精品成人在线| 91国产中文字幕| 精品国产超薄肉色丝袜足j| 亚洲 欧美一区二区三区| 成人三级做爰电影| 狠狠婷婷综合久久久久久88av| 亚洲成av片中文字幕在线观看| 亚洲视频免费观看视频| 人人妻,人人澡人人爽秒播 | 91九色精品人成在线观看| 精品久久久久久电影网| 国产免费现黄频在线看| 亚洲精品美女久久av网站| 久久毛片免费看一区二区三区| 不卡av一区二区三区| 18禁黄网站禁片午夜丰满| 亚洲激情五月婷婷啪啪| 永久免费av网站大全| 亚洲国产成人一精品久久久| 色视频在线一区二区三区| 一边亲一边摸免费视频| 国产成人一区二区三区免费视频网站 | 一二三四在线观看免费中文在| 日韩 亚洲 欧美在线| 日本一区二区免费在线视频| 国产欧美日韩一区二区三 | 久久人妻熟女aⅴ| 男人舔女人的私密视频| 一本综合久久免费| 久久99热这里只频精品6学生| 成人三级做爰电影| 精品一区二区三区四区五区乱码 | 啦啦啦在线观看免费高清www| 日本色播在线视频| 精品福利永久在线观看| 99精国产麻豆久久婷婷| 久久久久精品人妻al黑| 精品福利观看| 欧美精品av麻豆av| 日韩大码丰满熟妇| 国产精品成人在线| 十八禁网站网址无遮挡| 女性被躁到高潮视频| 色婷婷av一区二区三区视频| 超色免费av| 久久99精品国语久久久| 色94色欧美一区二区| 国产又色又爽无遮挡免| 午夜免费观看性视频| 亚洲精品久久成人aⅴ小说| 午夜精品国产一区二区电影| 中文字幕最新亚洲高清| 少妇人妻 视频| 中文欧美无线码| 久久99精品国语久久久| 日韩av不卡免费在线播放| 国产成人系列免费观看| 久久久久久人人人人人| 国产黄色视频一区二区在线观看| 亚洲成人手机| 日本欧美国产在线视频| 国产精品久久久久成人av| 亚洲欧美一区二区三区国产| 久久影院123| 一级毛片 在线播放| 国产黄频视频在线观看| 亚洲中文av在线| 91精品三级在线观看| 人体艺术视频欧美日本| 男女之事视频高清在线观看 | 丁香六月天网| 久久久久久免费高清国产稀缺| 亚洲专区中文字幕在线| 中文字幕亚洲精品专区| 99九九在线精品视频| 国产av国产精品国产| 啦啦啦啦在线视频资源| 国产视频首页在线观看| 国产激情久久老熟女| 国产麻豆69| 国产麻豆69| 欧美中文综合在线视频| 国产精品人妻久久久影院| 后天国语完整版免费观看| 国产一卡二卡三卡精品| 免费少妇av软件| 国产在线一区二区三区精| 免费在线观看视频国产中文字幕亚洲 | 美女视频免费永久观看网站| 男人爽女人下面视频在线观看| www.熟女人妻精品国产| 亚洲av电影在线观看一区二区三区| 9色porny在线观看| 最近手机中文字幕大全| 亚洲成人免费电影在线观看 | 午夜福利免费观看在线| 国产免费现黄频在线看| 母亲3免费完整高清在线观看| 人人妻人人澡人人看| 国产片内射在线| 在线观看免费高清a一片| 亚洲综合色网址| e午夜精品久久久久久久| 最黄视频免费看| 观看av在线不卡| 狂野欧美激情性xxxx| 2018国产大陆天天弄谢| 十八禁人妻一区二区| 嫁个100分男人电影在线观看 | svipshipincom国产片| 日韩av免费高清视频| 国产精品一区二区免费欧美 | 制服人妻中文乱码| 老司机深夜福利视频在线观看 | 国产1区2区3区精品| 欧美日韩国产mv在线观看视频| 亚洲自偷自拍图片 自拍| 人妻人人澡人人爽人人| 亚洲一区中文字幕在线| 精品少妇一区二区三区视频日本电影| 欧美+亚洲+日韩+国产| 黑人巨大精品欧美一区二区蜜桃| 亚洲伊人色综图| 一区二区三区激情视频| 丝袜在线中文字幕| 精品亚洲乱码少妇综合久久| 考比视频在线观看| 黄色 视频免费看| 亚洲精品成人av观看孕妇| 另类精品久久| 天天躁夜夜躁狠狠躁躁| tube8黄色片| 赤兔流量卡办理| 国产成人一区二区在线| 亚洲精品一二三| 亚洲av电影在线进入| 大型av网站在线播放| 伊人亚洲综合成人网| 国产成人影院久久av| 欧美亚洲 丝袜 人妻 在线| 欧美日韩视频精品一区| 免费在线观看黄色视频的| 国产成人免费观看mmmm| 久久久久视频综合| 国产精品欧美亚洲77777| 久久久精品94久久精品| 国产成人av激情在线播放| 午夜福利视频精品| 一本大道久久a久久精品| 水蜜桃什么品种好| 亚洲成人免费电影在线观看 | 搡老岳熟女国产| 欧美av亚洲av综合av国产av| 国产不卡av网站在线观看| 成人免费观看视频高清| 成人午夜精彩视频在线观看| av在线app专区| 少妇精品久久久久久久| 欧美黑人欧美精品刺激| 久久精品国产亚洲av高清一级| 精品一区二区三区四区五区乱码 | 婷婷色综合大香蕉| 首页视频小说图片口味搜索 | 亚洲成av片中文字幕在线观看| videosex国产| 亚洲久久久国产精品| av线在线观看网站| 老熟女久久久| 亚洲精品一二三| 成人亚洲欧美一区二区av| 亚洲成人国产一区在线观看 | 国产免费福利视频在线观看| 超色免费av| 狂野欧美激情性xxxx| 高潮久久久久久久久久久不卡| 国产在视频线精品| 久久人人爽人人片av| 国产精品香港三级国产av潘金莲 | 亚洲情色 制服丝袜| 一级毛片我不卡| 免费在线观看视频国产中文字幕亚洲 | 波多野结衣一区麻豆| 99精品久久久久人妻精品| √禁漫天堂资源中文www| 欧美性长视频在线观看| 免费黄频网站在线观看国产| 宅男免费午夜| 宅男免费午夜| 狠狠婷婷综合久久久久久88av| 18在线观看网站| 丰满人妻熟妇乱又伦精品不卡| 一本久久精品| 人人澡人人妻人| 久久国产精品人妻蜜桃| 咕卡用的链子| 国产亚洲av高清不卡| 一个人免费看片子| 19禁男女啪啪无遮挡网站| 美女大奶头黄色视频| 亚洲人成电影观看| 久久久国产一区二区| 热re99久久国产66热| 丝袜脚勾引网站| 少妇粗大呻吟视频| 亚洲第一av免费看| 一本一本久久a久久精品综合妖精| 国产av精品麻豆| 久久性视频一级片| 2018国产大陆天天弄谢| 亚洲av男天堂| 日韩av免费高清视频| 国产1区2区3区精品| 97人妻天天添夜夜摸| 纯流量卡能插随身wifi吗| h视频一区二区三区| 黄色一级大片看看| 亚洲激情五月婷婷啪啪| 日韩一卡2卡3卡4卡2021年| 成人影院久久| 咕卡用的链子| 999精品在线视频| 男女无遮挡免费网站观看| 看十八女毛片水多多多| 久久天躁狠狠躁夜夜2o2o | 亚洲欧美中文字幕日韩二区| 午夜免费成人在线视频| 久久久久精品人妻al黑| 电影成人av| 2018国产大陆天天弄谢| 99久久99久久久精品蜜桃| 婷婷色av中文字幕| 亚洲国产日韩一区二区| 99精国产麻豆久久婷婷| 婷婷色麻豆天堂久久| 亚洲av国产av综合av卡| 熟女少妇亚洲综合色aaa.| 美国免费a级毛片| 欧美成人午夜精品| 日韩欧美一区视频在线观看| 亚洲国产精品999| 男女边摸边吃奶| 操美女的视频在线观看| 国产免费视频播放在线视频| 欧美 日韩 精品 国产| 久久av网站| 满18在线观看网站| 一区在线观看完整版| 日韩中文字幕视频在线看片| 制服人妻中文乱码| 中文乱码字字幕精品一区二区三区| 亚洲伊人色综图| 人体艺术视频欧美日本| 最黄视频免费看| 一边亲一边摸免费视频| 日韩免费高清中文字幕av| 纯流量卡能插随身wifi吗| 欧美黑人欧美精品刺激| 大型av网站在线播放| 在线av久久热| 国产精品久久久久成人av| 搡老岳熟女国产| 国产伦人伦偷精品视频| 国产成人欧美在线观看 | 少妇人妻久久综合中文| 久久热在线av| 国产视频一区二区在线看| 91精品三级在线观看| 久9热在线精品视频| √禁漫天堂资源中文www| 欧美性长视频在线观看| 国产爽快片一区二区三区| 亚洲成av片中文字幕在线观看| 一区二区三区精品91| 欧美国产精品一级二级三级| 亚洲午夜精品一区,二区,三区| 三上悠亚av全集在线观看| 免费少妇av软件| 免费观看av网站的网址| 成人国语在线视频| 精品国产国语对白av| 日韩制服丝袜自拍偷拍| 男女免费视频国产| 18在线观看网站| 天堂8中文在线网| 丰满饥渴人妻一区二区三| 久久毛片免费看一区二区三区| 十八禁人妻一区二区| 黄色a级毛片大全视频| 国产精品久久久久久精品古装| 成人午夜精彩视频在线观看| 人成视频在线观看免费观看| 精品欧美一区二区三区在线| 精品少妇内射三级| 日韩制服丝袜自拍偷拍| 亚洲av日韩精品久久久久久密 | 亚洲国产欧美在线一区| 9191精品国产免费久久| 老熟女久久久| 午夜福利在线免费观看网站| 一区二区av电影网| 久久精品久久精品一区二区三区| 亚洲中文字幕日韩| 久久久久精品人妻al黑| 免费在线观看黄色视频的| 亚洲情色 制服丝袜| 麻豆乱淫一区二区| 中文字幕人妻丝袜一区二区| 少妇人妻 视频| 精品国产国语对白av| 国产精品一二三区在线看| av视频免费观看在线观看| 日韩大码丰满熟妇| 男女国产视频网站| 国产亚洲一区二区精品| 亚洲男人天堂网一区| 亚洲欧洲国产日韩| 国精品久久久久久国模美| 国产成人欧美| 性高湖久久久久久久久免费观看| 免费不卡黄色视频| 免费日韩欧美在线观看| 久久亚洲国产成人精品v| 两性夫妻黄色片| 啦啦啦在线免费观看视频4| 国产成人欧美| 九草在线视频观看| 日本vs欧美在线观看视频| 一本久久精品| 在线观看免费视频网站a站| 免费高清在线观看视频在线观看| 各种免费的搞黄视频| 国语对白做爰xxxⅹ性视频网站| 亚洲第一av免费看| 夫妻性生交免费视频一级片| 久久99精品国语久久久| 高清不卡的av网站| 美女视频免费永久观看网站| 十八禁网站网址无遮挡| 波多野结衣一区麻豆| 制服人妻中文乱码| 丝瓜视频免费看黄片| av网站免费在线观看视频| 日本av手机在线免费观看| 国产野战对白在线观看| 美女脱内裤让男人舔精品视频| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲欧美精品综合一区二区三区| 欧美黄色片欧美黄色片| av福利片在线| 欧美人与性动交α欧美软件| 国产精品久久久av美女十八| 国产精品国产三级专区第一集| 午夜福利视频精品| 国产xxxxx性猛交| 男女边摸边吃奶| av在线播放精品| 在线观看人妻少妇| 肉色欧美久久久久久久蜜桃| av网站在线播放免费| 欧美日韩黄片免| 性色av乱码一区二区三区2| 久久热在线av| 亚洲图色成人| 视频区欧美日本亚洲| av又黄又爽大尺度在线免费看| 老司机在亚洲福利影院| 另类亚洲欧美激情| 在线观看免费午夜福利视频| 丰满少妇做爰视频| 国产又爽黄色视频| 青草久久国产| 女人爽到高潮嗷嗷叫在线视频| 精品一区二区三卡| 欧美日韩视频精品一区| 免费一级毛片在线播放高清视频 | 咕卡用的链子| 精品一区二区三卡| 久久性视频一级片| 精品卡一卡二卡四卡免费| 在线观看人妻少妇| 视频在线观看一区二区三区| 精品第一国产精品| 91字幕亚洲| 18在线观看网站| 9191精品国产免费久久| 国产精品久久久久成人av| 美女扒开内裤让男人捅视频| 嫩草影视91久久| 欧美日韩成人在线一区二区| 久久久久网色| 色94色欧美一区二区| 成人影院久久| 欧美在线黄色| 99香蕉大伊视频| a 毛片基地| 久久女婷五月综合色啪小说| 搡老乐熟女国产| 日韩视频在线欧美| 少妇人妻 视频| cao死你这个sao货| 久久ye,这里只有精品| 亚洲国产毛片av蜜桃av| 精品国产乱码久久久久久小说| 免费日韩欧美在线观看| 亚洲人成77777在线视频| 久久青草综合色| 国产欧美日韩一区二区三 | 精品人妻熟女毛片av久久网站| 免费久久久久久久精品成人欧美视频| 少妇粗大呻吟视频| 男女边摸边吃奶| 欧美性长视频在线观看| 欧美激情极品国产一区二区三区| 黑人猛操日本美女一级片| 精品亚洲成国产av| 亚洲欧美成人综合另类久久久| 妹子高潮喷水视频| 麻豆av在线久日| a级毛片黄视频| 亚洲九九香蕉| 99久久综合免费| 宅男免费午夜| 日韩大片免费观看网站| 久久久久久免费高清国产稀缺| 男的添女的下面高潮视频| 咕卡用的链子| 亚洲中文av在线| 亚洲情色 制服丝袜| 国产午夜精品一二区理论片| 两人在一起打扑克的视频| 亚洲中文字幕日韩| 交换朋友夫妻互换小说| 王馨瑶露胸无遮挡在线观看| 久久久久精品人妻al黑| 少妇人妻久久综合中文| 热99国产精品久久久久久7| 精品人妻在线不人妻| 老司机午夜十八禁免费视频| 精品久久久久久电影网| 亚洲欧美激情在线| 欧美精品人与动牲交sv欧美| 高清不卡的av网站| 一本一本久久a久久精品综合妖精| 超色免费av| 欧美人与性动交α欧美软件| 欧美97在线视频| 日本wwww免费看| av电影中文网址| 国产欧美亚洲国产| 国产成人av教育| 一区二区三区激情视频| 亚洲国产日韩一区二区| 丰满人妻熟妇乱又伦精品不卡| 国产亚洲午夜精品一区二区久久| av线在线观看网站| 99国产精品一区二区蜜桃av | 精品少妇黑人巨大在线播放| 国产欧美日韩综合在线一区二区| a 毛片基地| 五月开心婷婷网| av天堂久久9| 亚洲国产毛片av蜜桃av| 久久九九热精品免费| 两个人看的免费小视频| 欧美变态另类bdsm刘玥| 久久久久视频综合| 久久久精品区二区三区| 国产精品久久久久久人妻精品电影 | 亚洲,欧美,日韩| 国产视频首页在线观看| 成人午夜精彩视频在线观看| 欧美日韩福利视频一区二区| 在线观看www视频免费| 人人妻人人澡人人爽人人夜夜| 真人做人爱边吃奶动态| 看免费成人av毛片| 日韩伦理黄色片| 欧美精品啪啪一区二区三区 | 女人精品久久久久毛片| 亚洲自偷自拍图片 自拍| 成在线人永久免费视频| 亚洲,一卡二卡三卡| 日韩电影二区| 91麻豆av在线| 国产精品久久久久久精品古装| 欧美日韩成人在线一区二区| 精品卡一卡二卡四卡免费| 久久久国产精品麻豆| 亚洲,一卡二卡三卡| 久久精品熟女亚洲av麻豆精品| 一本色道久久久久久精品综合| 另类精品久久| 亚洲美女黄色视频免费看| 天天躁日日躁夜夜躁夜夜| 免费高清在线观看日韩| 亚洲av电影在线观看一区二区三区| 国产91精品成人一区二区三区 | 国产亚洲午夜精品一区二区久久| 亚洲欧美清纯卡通| 久久毛片免费看一区二区三区| 婷婷色综合大香蕉| 国产老妇伦熟女老妇高清| 性色av乱码一区二区三区2| 免费久久久久久久精品成人欧美视频| 男人添女人高潮全过程视频| xxx大片免费视频| 妹子高潮喷水视频| 国产一区二区三区av在线| 一本综合久久免费| 亚洲第一青青草原| 欧美黄色淫秽网站| 韩国高清视频一区二区三区| 婷婷色麻豆天堂久久| 丁香六月欧美| 国产高清不卡午夜福利| 欧美日韩国产mv在线观看视频| 国产又色又爽无遮挡免| 成人影院久久| 无遮挡黄片免费观看| 久久久久视频综合| 少妇精品久久久久久久| 欧美亚洲 丝袜 人妻 在线| 国产精品国产三级国产专区5o| 国产一区二区激情短视频 | 中国美女看黄片| 色播在线永久视频| 99国产综合亚洲精品| 一区二区三区四区激情视频| 黄色一级大片看看| 少妇猛男粗大的猛烈进出视频|