亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        海水淡化用汽液分離元件研究進展

        2012-04-09 10:25:51張曉晨胡孔誠
        化學工業(yè)與工程 2012年2期
        關鍵詞:汽液流板氣速

        張曉晨,張 銘,胡孔誠

        (1.青島科技大學化學與分子工程學院,山東 青島,266042; 2.國家海洋局天津海水淡化與綜合利用研究所,天津 300192;3.國家海水利用工程技術研究中心,天津 300192)

        世界人口的持續(xù)增加和經濟的不斷發(fā)展,使得水資源短缺成為全球性問題。海水淡化作為一種淡水資源增量技術,可以增加淡水供應量,已成為世界各國解決淡水短缺的重要途徑之一。截至2010年,全球淡化水裝機規(guī)模已達6 520萬t/d,解決了2億多人的飲水問題[1]。

        目前,主流的商業(yè)化海水淡化技術包括:多級閃蒸(MSF)、低溫多效(MED)和反滲透(RO)3種。前2種均為熱法海水淡化技術,其本質是在一定條件下使海水蒸發(fā)產生蒸汽,蒸汽冷凝得到淡水的過程。海水蒸發(fā)時由于沸騰和液滴碰撞會使蒸汽中夾帶有一定量的海水液滴,這些液滴如不能有效去除將會影響產水水質。為了去除蒸汽中夾帶的海水小液滴,保證淡化水水質,在熱法海水淡化裝置中通常會設置汽液分離元件,用于分離蒸汽中夾帶的海水液滴。圖1a)和圖1b)分別給出了多級閃蒸和低溫多效2種主流熱法海水淡化系統(tǒng)中汽液分離元件的示意圖。

        圖1 海水淡化裝置中的汽液分離元件Fig.1 Gas-liquid separation in desalination plant

        目前用于熱法海水淡化工程中的汽液分離元件主要可分為折流板式(也稱葉片式)和絲網式2大類,其分離性能主要受操作工況和元件結構的影響,操作工況主要包括:操作氣速、液滴粒徑、液滴夾帶量等。折板式汽液分離元件的結構參數(shù)包括:折板間距、折板形式、角度、長度、安裝方式及折鉤設置情況等;絲網式汽液分離元件的結構參數(shù)則包括絲網層數(shù)、厚度、編制方式、網絲直徑等。本論文將就這兩類廣泛使用的汽液分離元件的相關研究進展和案例進行綜述,以期對相關研究和工程應用提供借鑒。

        1 折流板式汽液分離元件

        折流板式汽液分離元件是熱法海水淡化工程中應用最為廣泛的1種汽液分離元件,其具有分離效率高、流動阻力小、結構簡單、抗污堵能力強、操作彈性大等特點[2]。

        1.1 分離原理

        折流板式汽液分離元件由大量平行的金屬折流板組成,氣流在彎折流道中急速轉向,夾帶液滴在較大慣性力作用下與折流板碰撞被分離[3]。不同形式的折流板結構如圖2所示。

        圖2 折流板式汽液分離元件結構圖Fig.2 Structure of wave-plate mist eliminator

        1.2 研究進展

        折流板式汽液分離元件的性能主要受折流板形式、操作氣速、板間距等因素的影響。為了獲得最適宜的折流板式汽液分離元件結構,國內外研究者采用試驗和數(shù)值模擬的方法對折流板汽液分離過程開展了眾多研究工作。

        1.2.1試驗研究

        早在1939年Houghton等[4]就對折流板式汽液分離元件進行了研究,他們指出,折流板式汽液分離元件的分離效率隨氣速的增加呈現(xiàn)先增大后減小的趨勢,并且認為產生這一現(xiàn)象的主要原因是二次夾帶。

        顏學升[5]和楊柳等[6]對V形折流板進行了研究,考察了操作條件和結構參數(shù)對分離效率的影響,結果表明,液滴粒徑增加使液滴弛豫時間增加,有利于液滴分離;而折板角度增加,使氣流流線趨于直線形;板間距增加,則增加了液滴顆粒的流通面積,都不利于液滴分離。

        進一步的研究發(fā)現(xiàn),在折流板上設置折鉤(見圖2e)可起到導流的作用,能有效降低二次夾帶的發(fā)生,提高分離效率[9-10]。

        李嘉等[11]對帶鉤折流板和普通折流板的性能進行了對比試驗,研究結果表明,在相同的氣速、折流板板間距和角度的情況下,帶鉤折流板的分離性能要優(yōu)于普通折流板,原因是折鉤處液膜的在表面張力作用下厚度變薄,減少了因氣流對液膜剪切作用引起的二次夾帶;另一方面,折鉤處產生的旋渦會對小粒徑的液滴產生吸引作用,使之進入折鉤被分離,從而提高了分離效率。

        王為術等[12]也針對帶鉤折流板式汽液分離元件進行了試驗研究,結果表明,氣速是影響分離效率的主要因素,分離效率隨氣速的增加呈先增大后減小的趨勢,氣速超過5 m/s后分離效率急劇降低,這主要是由于氣速越大,氣流對液滴的攜帶力越強,在折流板上來不及排出而被氣流攜帶走的液體就越多。此外,試驗結果還表明相同氣速條件下,提高載液量可以增加氣流中液滴相互碰撞的幾率,進而提高分離效率。

        肖立春等[13]的研究認為,折角為72 °,板間距為20 mm時,帶鉤折流板分離性能最佳,而折流板的長度對分離效率影響并不明顯。

        在大量試驗研究基礎上,人們也設計了一系列不同結構形式的流板。

        鐘利春[14]設計了一種表面具有凹槽的波紋型折流板,這種結構一方面可以增加液滴同折流板壁碰撞的幾率和強度,使液滴易于被折流板攔截;另一方面被攔截的液滴可在凹槽內匯聚,并沿凹槽流出,極大的減少了二次夾帶的發(fā)生,從而提高了分離效率。

        李明慧[15]設計了一種兩級折流板元件,第1級主要起分離作用,而第2級除了分離作用外,還可以起到部分導流的作用,減少葉片下游的二次夾帶,從而達到提高分離效率和降低壓降的效果。

        徐月等[16]通過優(yōu)化多層折流板的布置,使流體在通道內交錯流動,實現(xiàn)了折流板出口處氣流的均勻分布,提高了裝置的分離效率。

        阮國嶺等[17]針對海水淡化設備規(guī)模大,不易清洗的特點,開發(fā)了一種用于熱法海水淡化的汽液分離元件,折流板被固定于開口槽架上,這種結構使折流板便于拆卸,易于更換和清理。

        1.2.2數(shù)值模擬研究

        近年來,計算流體力學(CFD)的方法開始應用于汽液分離元件的研究,研究者借助于數(shù)值模擬的方法,對汽液分離過程進行了一系列的模擬研究。

        Bürkholz[18]在1989年提出分離效率會受到液滴粒徑、折流板間距和角度等因素的影響,之后的研究者們對這些影響因素的作用機理進行了大量的研究。

        徐淑君等[19]基于RNGk-ε模型,采用汽液兩相雙向耦合方法對單級V形折流板式汽液分離元件進行了研究。研究顯示液滴粒徑對分離效率的影響非常大,進氣速度為3 m/s時,粒徑為50 μm液滴的去除率可達100%,而粒徑小于20 μm的液滴去除率還不到35%。這主要是由于小粒徑液滴受Saffman升力作用明顯,顆粒運動軌跡不穩(wěn)定,不易被去除。

        趙毅等[20]采用k-ε湍流模型對折線形折流板式汽液分離元件[見圖2c)]的分離性能進行了研究,板間距為20 mm、氣速到達6 m/s時,折流板對粒徑20 μm以上液滴的分離效率近100%。同時,模擬結果顯示,分離效率還隨折流板偏轉角增加而增加,偏轉角在51°左右時,折流板對大于30 μm的液滴分離效率可以達到100%。

        孫志春等[21]采用雷諾應力模型,氣-液雙向耦合算法,研究了V型折流板的分離性能,結果同樣顯示,折流板間距為20 mm時對應的分離效率最高。

        人們通過對STDk-ε模型(標準k-ε湍流模型)、低Rek-ε模型、RSM(Reynolds Stress Model,雷諾應力模型)、SSTk-ω(剪應力傳遞湍流模型)和RNGk-ε(修正的k-ε雙方程模型)5種常用湍流模型的對比研究[19, 22-25],并且與試驗數(shù)據(jù)進行比較,研究結果表明RNGk-ε湍流模型能較好地處理強旋及彎曲流動的情況,而低Rek-ε模型則較適于對近壁面流動情況進行模擬。

        對如圖2d)所示的折流板內流場分布的模擬研究結果顯示[26]:液滴主要在第1折彎處的左壁面和第2拐彎處的右壁面被分離。這主要是由于離心力的作用,液滴由內彎側向外彎側聚集,大顆粒液滴主要集中在外彎側,尤其在第1個拐彎處,液滴的聚集現(xiàn)象更加明顯。

        孫志春等[21]將冷態(tài)試驗結果與CFD模擬結果進行了對比,發(fā)現(xiàn)液滴粒徑分布為10 μm,氣速在5~7 m/s時的數(shù)值模擬結果與試驗結果非常接近,作者認為這是由于10 μm液滴粒徑接近亞微觀粒子,與計算模型中的粒子吻合度較好所致。

        2 絲網式汽液分離元件

        絲網是另外一種應用較廣的汽液分離元件,具有體積小、分離效率高、質量輕及安裝操作維修方便等優(yōu)點。

        2.1 分離原理

        絲網式汽液分離元件的結構如圖3所示,其分離原理主要有3種:直接攔截、慣性撞擊和擴散攔截。直接攔截就是氣體流過絲網結構時,氣流中大于絲網孔徑的液滴被攔截下來;慣性撞擊是利用液滴慣性比氣流大的特點,當夾帶有液滴的氣流流經絲網時,氣體順著絲網結構改變方向,而液滴由于慣性作用撞擊到絲網上被分離出來;擴散攔截是針對粒徑小于0.1 μm的小液滴,這種液滴以不規(guī)則形式沿著流體流線運動,碰到絲網被富集分離下來[27]。

        2.2 研究進展

        絲網式汽液分離元件廣泛應用于化工、食品、制冷和海水淡化等領域,影響其分離效率和壓降的因素主要包括:絲網厚度、操作氣速和編制方式等。

        圖3 絲網式汽液分離元件結構圖Fig.3 Structure of wire mesh eliminator

        2.2.1試驗研究

        19世紀50和60年代,York和Poppele[28-29]的論文中對絲網式汽液分離元件進行了研究,并且歸納了一個計算最高操作氣速的半經驗公式,用以對絲網的設計進行指導。之后,Carpenter等[30]對絲網式汽液分離元件進行了較為系統(tǒng)的研究,并對絲網的層數(shù)、網絲直徑及比表面積等結構參數(shù)進行考察,得到分離效率的計算公式。

        Hisham等[31]在研究中對操作氣速、絲網厚度和編制方式等影響因素進行了全面的考察,結果顯示:絲網式汽液分離元件的壓降隨氣速增加呈線性上升的趨勢;分離效率在氣速為1~6 m/s的范圍內隨氣速的增大而增加,高氣速下分離效率變化趨于平緩;網絲的厚度對分離效率影響不大,但容易造成壓降的增加;密織的絲網雖然有利于液滴的分離但同時也容易引發(fā)二次夾帶。

        李柏松等[32]以水-空氣體系為研究對象,對水平安裝的絲網式汽液分離元件進行了研究,考察了液滴粒徑和液滴濃度對絲網分離效率的影響。研究結果表明:分離效率隨液滴粒徑的增大而增大,液滴粒徑為10 μm時,對應的分離效率為90%,而粒徑增大到70 μm時,分離效率接近100%;分離效率也隨液滴濃度的增大而增大,這主要是由于液滴濃度的增大,使得液滴之間容易發(fā)生碰撞而凝集成較大液滴,從而容易被分離。

        Capps[33]對絲網式汽液分離元件的結構參數(shù)進行了研究。結果表明,網絲直徑80~280 μm、厚度65~150 mm,比表面積為250~650 m2/m3,空隙率在97~99%之間的絲網對10 μm左右液滴的分離效率可達到99%。

        Bürkholz[34]針對絲網液滴二次夾帶問題行了研究,認為氣速控制在4~5 m/s可以減少夾帶造成的效率下降。同時還根據(jù)試驗結果推導出了一個簡單的近似公式,通過絲網結構來確定可分離的液滴粒徑范圍。

        Lerner[35]開發(fā)了一種用于多級閃蒸(MSF)海水淡化裝置中的絲網式汽液分離元件,其體積較小,但有效通量可以達到傳統(tǒng)絲網的2~3倍。

        基于試驗研究結果,人們也設計開發(fā)出了一系列新型的絲網式汽液分離元件。洪放[36]發(fā)明了一種可調的絲網式汽液分離元件,在進氣側安裝了可調節(jié)風門,使用過程中可根據(jù)需要對風門面積進行調節(jié),使系統(tǒng)在最適宜氣速下工作,以保證其達到最高的分離效率。唐曉東等[37]則發(fā)明了一種汽液分離裝置,利用特殊絲網使液滴速度降低以提高絲網對液滴的捕捉收集能力,并通過氣體和液滴透過絲網柔性層的多次折流提高其分離效率。陳長林[38]設計了一種波形絲網除沫器,其由經模壓后的絲網和絲網格柵組成,這種結構可以使分離下來的液滴迅速匯流到低處,聚集成大液滴后排出,可有效減少液滴在絲網上的滯留時間,避免發(fā)生二次夾帶降低分離效率。王玉潔等[39]將絲網與板式聚結填料相結合,先通過板式聚結填料進行初步分離,再經由絲網進行二次分離,不僅減少了二次夾帶的發(fā)生,提高了分離效率,還可以有效減少對設備的腐蝕,延長其使用壽命。

        2.2.2數(shù)值模擬研究

        研究者同樣利用CFD模擬的方法對絲網式汽液分離元件進行了研究。

        Rahimi等[40]利用FLUENT對絲網式汽液分離元件的分離性能進行了計算。結果顯示:當氣速為6 m/s時,分離效率接近100%,但當氣速繼續(xù)增加時,本已凝結在網絲上的液滴又重新被高速氣流夾帶,造成分離效率的下降。絲網網格的密度同樣會影響分離效率,較大的網格使液滴不能很好的和網絲接觸,以致無法凝聚在絲網上,導致分離效率較低。

        張李偉等[41]選取20層隨機交錯排列的絲網組成的絲網汽液分離器模型進行數(shù)值模擬,并利用經驗公式的計算結果進行驗證。結果表明,在分離效率大于80%情況下的擬合情況較好,而對于分離效率較低的情況,計算值與試驗值仍存在較大的偏差。

        由于絲網式汽液分離元件的內部結構復雜,流動情況多變,CFD模擬過程中對其流場進行的簡化,會對模擬結果的準確性造成很大的影響,所以對于絲網式汽液分離元件的CFD模擬研究尚不成熟,仍需研究者對數(shù)值模型進行進一步的優(yōu)化。

        3 工程實例

        折流板式和絲網式汽液分離元件各有優(yōu)缺,前者具有通量大、壓降小、易清洗、不宜堵塞等優(yōu)點,但成本較高;后者分離效率高、安裝簡單、成本較低,但操作彈性小,易被堵塞。因此在實際工程中需綜合考慮海水淡化工藝流程、操作溫度、原水水質等情況進行選擇。表1給出了部分國內外熱法海水淡化工程中所用汽液分離元件。

        表1 國內外部分多效蒸餾海水淡化裝置用汽液分離元件Table 1 Application of gas-liquid separation to MSF

        4 發(fā)展方向

        折流板和絲網作為兩種主要海水淡化工程用汽液分離元件形式,已有多年的應用,今后仍將是主要應用形式。未來海水淡化用汽液分離元件的研究將主要集中在以下幾方面:

        1)開發(fā)新型汽液分離元件,在保證裝置分離效率的前提下,不斷降低操作壓降。

        2)將玻璃鋼、工程塑料、鋁板和新型新材料應用于折板流式和絲網式汽液分離元件,以不斷降低汽液分離元件成本。

        3)完善數(shù)值模擬模型,提高CFD模擬結果的準確性。

        4)通過研究建立不同形式汽液分離元件的數(shù)學模型,用于指導工程設計。

        參考文獻:

        [1]IDA Desalination Yearbook 2010-2011[R].UK:Media Analytics Ltd,2010

        [2]MOORE M J,SIEVERDING C H.Two-phase steam flow in turbines and separators[M].Washington:Hemisphere Publishing Corporation,1976

        [3]任相軍,王振波,金有海.氣液分離技術設備進展[J].過濾與分離,2008,18(2):43-47

        [4]HOUGHTON H G,RASFORD W H.Measurements on eliminators and the development of a new type for use at high gas velocities[J].Transactions of the American Chemical Engineers,1939,35:427-433

        [5]顏學升,王助良,吳健衛(wèi).脫硫除霧設備性能研究[J].農機化研究,2003(7):136-138

        [6]楊柳,王世和,王小明.脫硫除霧器除霧特性的研究[J].中國動力工程學報,2005,25(2):289-292

        [7]MONAT J P,MCNULTY K J,MICHELSON I S,etal.Accurate evaluation of chevron mist eliminators[J].Chemical Engineering Progress,1986,82:32-39

        [8]AZZOPARDI B J,SANAULLAH K S.Re-entrainment in wave-plate mist eliminators [J].Chemical Engineering Science,2002,57: 3 557-3 563

        [9]MCNULTY K J,MONAT J P,HANSEN O V.Performance of commercial chevron mist eliminators[J].Chemical Engineering Progress,1987,83:48-55

        [10]JAMES P W,WANG Y,AZZOPARDI B J,etal.The role of drainage channels in the performance of wave-plate mist eliminators[J].Chemical Engineering Research & Design,2003,81:639-648

        [11]李嘉,黃素逸,王曉墨,等.波形板分離器的冷態(tài)試驗試驗研究[J].華中科技大學學報:自然科學版,2008,36(1):112-114

        [12]王為術,徐維暉,羅毓珊,等.帶鉤波形板分離器分離性能試驗研究[J].華北水利水電學院學報,2010,31(2):52-54

        [13]肖立春,丁志江,李強,等.汽水分離器分離效率的冷態(tài)試驗試驗研究[J].熱能動力工程,2010,25(2):177-179

        [14]鐘利春.平行波形除霧器:中國,20062000289.8[P].2007-01-24

        [15]李明慧.波形折流板式除霧器:中國, 2006 2016 5207.5[P].2008-01-09

        [16]徐月,付嫚,覃亮,等.一種疊式折流板脫水器:中國,201020132060.6[P].2010-11-24

        [17]呂慶春,阮國嶺.用于蒸餾法海水淡化捕沫裝置:中國,200520025714.4[P].2006-05-24

        [18]BURKHOLZ A.Droplet Separation[M].New York:VCH publishers,1989

        [19]徐淑君,姚征,朱懿淵.波紋板除霧器兩相流動的數(shù)值模擬與分析[J].上海理工大學學報,2007,29(3):275-280

        [20]趙毅,華偉,王亞君,等.濕式煙氣脫硫塔中折線型擋板除霧器分離效率的數(shù)值模擬[J].中國動力工程學報,2005,25(2):293-297

        [21]孫志春,郭永紅,肖海平,等.鼓泡脫硫塔除霧器除霧特性數(shù)值研究及實驗驗證[J].中國電機工程學報,2010,30(8):68-75

        [22]WANG W,DAVIES G A.CFD studies of separation of mists from gases using vane-type separators [J].Transactions of the Institute of Chemical Engineers,Part A,1996,74:232-238

        [23]GILLANDT I,RIEHLE C,F(xiàn)RITSCHING U.Gas-particle flow in a comparison of measurements and simulations[J].Forschung Im Ingenieur Wesen-engineering Research,1996,62(11/12):315

        [24]郭鵬宇,楊震.除霧器通道內二維兩相流場的數(shù)值模擬[J].電力科學與工程,2005,(2):31-33

        [25]樊水沖,楊學忠,郭勇.波紋板除霧器流場與效率的數(shù)值計算[J].煉油與化工,2006,17(2):5-9

        [26]GALLETTI C, BRUNAZZI E, TOGNOTTI L.A numerical model for gas flow and droplet motion in wave-plate mist eliminators with drainage channels [J].Chemical Engineering Science,2008,63:5 639-5 652

        [27]史永紅.絲網氣液分離器分離性能和壓力損失分析[J].石油化工設備,2006,35(31):35-39

        [28]YORK O H.Performance of wire-mesh demisters [J].Chemical Engineering Progress,1954,50(8):421-424

        [29]YORK O H,POPPELE E W.Wire mesh mist eliminators[J].Chemical Engineering Progress,1963,59 (6):45-50

        [30]CARPENTER C L,OTHMER D F.Entrainment removal by a wire-mesh separator[J].AIC E Journal,1955,1:549-557

        [31]HISHAM T,EL DEDDOUKY,IMAD M ALAFFIQI,etal.Performance of wire mesh mist eliminator [J].Chemical Engineering and Processing,2000,39 (2):129-139

        [32]李柏松.天然氣過濾器氣液分離性能的試驗試驗研究[J].天然氣工業(yè),2007,27(10):123-125

        [33]CAPPS R W.Properly specify wire-mesh mist eliminators [J].Chemical Engineering Progress,1994,74(4):49-55

        [34]BURKHOLZ A.Die Beschreibung der partikelabscheidung durch tr?gheitskr?fte mit hilfe einer dimensionsanalytisch abgeleiteten kennzahl [J].Chem Eng Techol,1986,58:548-556

        [35]LENRNER B J.High-tech mist elimination in multi-stage evaporators[J].Plant Operations Prog,1986,(5):52-56

        [36]洪放.可調絲網除沫器:中國,95224980.4[P].1996-07-17

        [37]唐曉東,王志紅,諸林.一種天然氣脫硫脫水凈化裝置:中國,200610020946.X[P].2006-10-11

        [38]陳長林.波形絲網除沫器:中國,97236744.6[P].1999-01-27

        [39]王玉潔,李海平,王靜,等.一種捕沫分離裝置:中國,200910260116.8[P].2010-06-09

        [40]RAHIM I R,ABBASPOUR D.Determination of pressure drop in wire mesh mist eliminator by CFD[J].Chemical Engineering and Processing,2008(47):1 504-1 508

        [41]張李偉,孫海鷗,蘊蕾.船用絲網除霧器分離效率計算[J].船舶工程,2007,29(5):1-8

        猜你喜歡
        汽液流板氣速
        汽液通道結構對超薄熱管性能的影響分析
        結構參數(shù)對交錯百葉折流板管殼式換熱器性能影響的研究
        新型錐形孔折流板管殼式換熱器數(shù)值模擬研究
        水管內汽液兩相流壓降計算研究
        濕法煙氣脫硫吸收塔阻力特性實驗研究
        新型折板除霧器的流場和壓降數(shù)值模擬
        螺旋折流板換熱器在高溫高壓環(huán)境中的應用研究
        當代化工(2015年12期)2015-10-27 01:42:29
        高壓環(huán)境下垂直管內汽液攪混流相界面特征
        環(huán)己烯、1,2-環(huán)氧環(huán)己烷和2-環(huán)己烯-1-酮二元和三元系的等壓汽液平衡
        化工科技(2014年6期)2014-06-09 06:43:26
        D120 mm流化床冷模實驗研究
        化工科技(2014年5期)2014-06-09 05:17:22
        在线免费午夜视频一区二区| 中文精品久久久久人妻不卡| 粗大猛烈进出白浆视频| 在线综合亚洲欧洲综合网站 | 午夜精品久久久久成人| 亚洲V在线激情| 蜜臀精品一区二区三区| 青青草在线免费观看在线| 日韩久久无码免费看A| 五月激情婷婷丁香| 亚洲国产精品国自产拍av在线| 日本五十路人妻在线一区二区| 中国老熟女重囗味hdxx| 亚洲av无码专区亚洲av桃| 国产内射视频在线播放| 日韩有码在线观看视频| 久久国产劲爆∧v内射-百度| 亚洲乱亚洲乱少妇无码99p | 精品国产av一区二区三区| 99久久人人爽亚洲精品美女| 日本一区二区啪啪视频| 亚洲免费人成网站在线观看| 精品一区中文字幕在线观看| 中国丰满熟妇xxxx性| 精品三级久久久久久久电影| av天堂吧手机版在线观看| 久久精品国产亚洲av蜜点| 99精品国产在热久久| 国产小屁孩cao大人| 一区二区久久精品66国产精品| 澳门蜜桃av成人av| 人妻丝袜av中文系列先锋影音| 一本久到久久亚洲综合| 国产免费人成网站在线播放| 久久久精品国产亚洲av网深田| 色噜噜狠狠一区二区三区果冻| 野外三级国产在线观看| 中文字幕高清视频婷婷| 国产午夜福利精品一区二区三区| 国产精品毛片无遮挡高清| 区二区三区亚洲精品无|