999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

纏繞式變頻電磁水處理器電磁頻率對抑垢效果的影響

2015-10-15 05:35:37王建國李雨通鄧麗娟
化工學報 2015年3期
關鍵詞:信號實驗

王建國,李雨通,鄧麗娟

?

纏繞式變頻電磁水處理器電磁頻率對抑垢效果的影響

王建國,李雨通,鄧麗娟

(東北電力大學自動化工程學院,吉林省吉林市132012)

針對纏繞式變頻電磁水處理系統,利用ANSYS仿真軟件對處理腔內水溶液的磁感應強度和感生電流進行分析,探究水體內磁感應強度以及感生電流做功與激勵信號頻率的內在關系,并利用電磁抑垢效果在線評價實驗臺,通過改變纏繞式處理腔中電磁場的頻率大小分析其對特定硬水循環換熱實驗管道抑垢效果的影響規律,進而尋求理想的抑垢頻率。實驗結果表明,電磁頻率在分界頻率附近時水處理系統能到達較長的污垢誘導期以及較好的抑垢率,抑垢效果最佳。

纏繞式處理腔;頻率;磁感應強度;感生電流;抑垢效果;對流;傳熱;數值分析

引 言

電力、鋼鐵、冶金等工業領域大型設備的正常運行都離不開循環水冷卻系統。碳酸鈣垢是水冷應用中最常見的一種垢型,由于污垢的熱阻相對較大,換熱面結垢通常會降低其換熱效率,同時流通面積的減小會增加泵耗。結垢也會增加補水、排污水以及清理換熱面的費用,嚴重時會造成爆管等事故[1]。相對于傳統的化學除垢抑垢方法,電磁抑垢法以其操作簡單、綠色環保、價格低廉等優點廣泛應用于工業和生活水冷系統[2-6]。但是考慮到實驗條件的限制以及水體的多樣性,國內外所報道的研究成果存在差異性[7-12]。針對纏繞式變頻水處理系統的研究往往是在特定水質參數和實驗工況下進行,所以得到的結論也不盡相同。姜德寧等[13]對不同硬度的溶液樣本經過處理后,認為所加激磁電壓幅值越高阻垢效果越好,當磁感應強度達到0.9 mT時結垢量最小;金貞花等[14]對硬度在300 mg·L-1的硬水溶液施加交變電信號進行處理,認為激磁線圈內的電流越大阻垢效果越顯著;陳璨等[15]認為頻率為1.2 kHz左右時抑垢效果最好。而邢曉凱等[16]在自行設計的實驗條件下,認為激磁頻率為700 Hz時抗垢效果最佳;費繼友等[17]還有些學者通過自制的掃頻裝置,認為掃頻可以對不同水質達到一定的抑垢效果。

本研究針對電磁能量分布比較均勻[18]的纏繞式處理腔,從腔體水溶液內感生電流做功值和磁感應強度兩個方面開展研究,根據實驗臺實際條件繞制合理的線圈匝數,在自主研制的電磁抑垢效果在線評價系統的基礎上得出不同頻率下特定硬水循環流經換熱實驗管道時的抑垢效果及其變化規律,進而對纏繞式變頻電磁水處理技術的發展起到指導性作用。

1 纏繞式變頻電磁水處理系統電磁場分析

由于污垢的形成受到很多未知因素的影響,目前變頻電磁水處理的原理尚不明確,磁抑垢機理一直存在爭議[19],公認度比較高的有“洛倫茲力理論”[20]、“氫鍵斷裂理論”[20]、“磁滯效應理論”[21]、“極化作用理論”[22]等。綜合這些理論,可以將纏繞式變頻電磁水處理工作的主要原理概括為:電磁信號發生裝置對緊密纏繞在輸水管上的線圈施加一定頻率和幅值的方波電壓,從而使線圈內部產生與水流方向大致平行的交變磁場,交變磁場會在水體內部感生出與線圈纏繞方向平行的感生電場,輸水管道內的水體在交變磁場和電場的作用下內部離子和分子的物理化學性質發生某些改變,進而影響晶核的形成以及沉淀顆粒的生長,從而達到抑垢的效果。但是,在抑垢除垢的過程中感生電場和感生磁場的哪一個起到關鍵作用并無定論。韓勇等[23]認為感生電場對水體做功越多抑垢效果越好;而金貞花等[14]認為激磁線圈內的電流越大,磁感應強度越大,阻垢效果越顯著。本研究利用ANSYS有限元軟件探究感生電流做功值和磁感應強度大小與激勵線圈結構和信號頻率的內在關系,實現最優抑垢效果的分析。

ANSYS有限元軟件包是一個多用途的有限元法計算機設計程序,可以用來求解結構、流體、電力、電磁及碰撞等問題[24]。本研究基于纏繞式變頻電磁水處理系統的物理模型,利用ANSYS軟件的PLANE53、CIRCU124、POST1和POST26等單元對處理腔進行2D瞬態電磁仿真分析。選取10 mmol·L-1碳酸鈣過飽和溶液進行模型仿真,溶液pH為9.0,電導率為2.400×103μS·cm-1,溫度為27℃。激磁線圈為單層密繞的螺線管,線圈材料為漆包銅線,銅線的電阻率在外界溫度25℃時為1.7×10-8Ω·m,銅線直徑為2 mm,匝數為50。輸水管道采用PVC材料,壁厚2.5 mm,內徑為10 mm,管長0.25m。纏繞式變頻電磁水處理器2D模型如圖1所示,2D仿真模型如圖2所示。

圖1 纏繞式變頻電磁水處理器2D模型

圖2 纏繞式變頻電磁水處理器2D仿真模型

1.1 電磁頻率對水處理系統內感生電流的影響

由于纏繞式變頻電磁水處理系統水體內的感生電流方向與激勵線圈內電流方向平行,也呈環狀分布,水體感生電流密度分布如圖3所示,藍色部位為感生電流分布,粉色部位為激勵線圈。

圖3 水體感生電流密度分布

可以將處理腔內溶液等效為與激勵線圈同軸的單匝線圈進行處理,并達到很好的精度[23]。經ANSYS對所建模型的仿真,如圖4線圈電流信號波形和圖5感生電流信號波形,可知兩種信號都呈周期性指數變化。

圖4 線圈電流信號波形

圖5 感生電流信號波形

可根據線圈互感的電路方程推導出水體內感生電流的數學表達式[式(1)、式(2)]以及感生電流做功的計算公式[式(3)]。

(1)

(2)

式中,1、2分別為激勵線圈電感、水體電感, H;1、2分別為激勵線圈電阻水體電阻,Ω;為處理腔內水容積與循環水總體積的比值;0為實驗周期,s;為激勵線圈的電壓幅值,V。

根據做功公式[式(3)]可知,當參數以及激勵線圈結構一定時,感生電流做功量只與激勵電壓以及激勵信號頻率有關,當信號頻率確定時做功量隨激勵信號電壓幅值增大而增大。

當激勵線圈電壓幅值為1 V時,取參量為 1×10-12,線圈參量11的比值為0.667 ms,代入式(3),觀察電磁頻率對感生電流做功量的影響曲線,如圖6所示。可以看出,當激勵電壓幅值和溶液參數以及激勵線圈結構一定時,感生電流做功量只取決于激勵信號頻率的大小,當頻率在0~525 Hz附近時做功值隨頻率增加迅速增大,而當頻率大于525 Hz后做功值增長很不明顯。

圖6 電磁頻率對感生電流做功的影響

通過以上的討論可知,一味地增大激磁信號頻率并不能很明顯地增加做功量,雖然通過式(3)的函數表達式可以證明是關于的單調遞增函數,但是當高于某一個分界頻率值時,隨著頻率的繼續增大,的變化卻很緩慢,并且當趨近無窮大時收斂于某一個值。

由于式(3)很復雜,直接通過數學推導很難找到分界頻率與線圈時間常數的函數關系,所以選取幾組不同的值,找到其對應的分界頻率點進行擬合,找到對應的函數關系,擬合曲線如圖7所示。通過圖7可以知道臨界頻率0是關于的反比例函數,比例系數為0.35,所以0僅由決定,即激勵線圈是由電阻和電感決定。由于纏繞式變頻電磁水處理器結構參數的限制,分界頻率值往往處于低頻段內,所以纏繞式變頻電磁水處理更適合在低頻段工作。因為過高的頻率不僅不會使有顯著的提升,而且還會使纏繞式變頻電磁水處理系統設計復雜化,特別是造成線圈嚴重發熱以及電磁輻射等弊端,提高系統的造價,影響系統電路的精確度。

圖7 線圈時間常數與臨界頻率的關系

之所以會出現分界頻率,是由于溶液中的離子分布可以看作是由若干個離子氛組成的,當離子氛的電平衡遭到破壞,中心離子要從先前的離子氛中移出時,就會對離子溶液造成大的擾動作用。而纏繞式變頻電磁水處理裝置工作時,其產生的交變電場作用于碳酸鈣水溶液,當電磁頻率較低時離子氛的中心離子可以從離子氛中移出,致使大范圍的帶電離子重新排列,隨著頻率的增加離子氛的中心離子能夠更快速地從離子氛中移出,而當頻率高至迫使中心離子恰好移出其離子氛時對溶液的擾動可以達到最理想的效果,但是隨著頻率的繼續增加,將會導致有些離子氛的中心離子只能在其離子氛內做往返運動,擾動作用不會明顯提升。所以確實會存在一個分界頻率對溶液的做功值影響最明顯,擾動作用也最強。

1.2 電磁頻率對水處理系統內磁感應強度的影響

不同線圈匝數條件下,改變激勵信號頻率,觀察水體中心磁感應強度的大小,可以得到表1。由表可知,在激勵信號電壓幅值和線圈匝數不變的前提下,水體中心磁感應強度隨激勵信號頻率的增大也呈現遞減趨勢。

表1 水體中心磁感應強度與線圈匝數以及信號頻率的關系

磁感應強度計算公式為

式中,為磁感應強度,Gs;為磁場強度,A·m-1;為線圈匝數;為電流,A;e為有效磁路長度,m;為磁導率,H·m-1。

由式(5)和式(6)可知,在線圈的激勵電壓幅值一定時磁感應強度的大小與線圈內有效電流直接相關,然而有效電流由線圈阻抗決定,當線圈匝數一定時,信號頻率增大,阻抗會大幅度增大,所以會使磁感應強度減小。

2 實 驗

2.1 實驗系統

電磁抑垢效果在線評價實驗平臺結構如圖8 所示。

圖8 電磁抑垢效果在線評價實驗臺

1—water cooling system; 2—air cooling system; 3—water tank; 4—circulating water pump; 5—drainage system; 6—pipeline of A side; 7—pipeline of B side; 8—spiral winding processing chamber; 9—electromagnetic signal generator; 10—thermostatic water bath; 11—U-shaped tubes; 12—water quality parameters on-line instrumentation; 13—inlet temperature of two pipes; 14—outlet temperature of two pipes; 15—wall temperature of two pipes; 16—bath temperature

實驗系統采用雙路對比的實驗方法,研究相同水質、相同工況、加電磁處理和不加電磁處理的污垢熱阻及水質參數變換特性。主要由工質循環系統、冷卻循環系統、信號采集系統和變頻電磁水處理裝置4部分組成。工質循環系統由恒溫水浴箱、工質循環管道組成。恒溫水浴箱底部安裝有加熱器,用來保持箱內水溫恒定至設定溫度。通過變頻器控制循環水泵,獲得實驗需要的流速,同時方便進行不同流速的實驗。冷卻循環系統由空冷系統和水冷系統組成。在水冷系統工作的基礎上,空冷系統通過PID智能溫度控制器和變頻器控制風冷冷卻器的風機轉速來調節風量大小,從而控制循環工質的入口溫度,使入口溫度保持恒定。信號采集系統由工控機、智能控制模擬量前端、水質參數在線儀表及DS18B20測溫元件組成。變頻電磁水處理裝置由變頻電磁信號發生器和纏繞式電磁處理腔組成。變頻信號發生器產生可變頻率的大功率方波信號,通過纏繞電極式處理腔作用于流過的工質,從而達到抑垢效果。

2.2 污垢熱阻監測及抑垢率計算

基于污垢熱阻在線監測技術評價電磁場的抑垢效果。

污垢熱阻[25]的定義如下

式中,f為污垢熱阻,m2·K·W-1;wf1wf2wf3為3個測點的壁溫,即管壁與污垢間的界面實測溫度,℃;s為污垢與流體的界面溫度,℃;為熱通量,W·m-2。

s可通過式(8)求得

式中,fo、fi分別為實驗管段出、入口處溫度,℃;為實驗管道長度,m;為實驗管道內徑,m;為Stanton數。

式中,c為介質比熱容,J·kg·K-1;為介質密度,kg·m-3。

綜合以上分析可知,只要測出工藝流體在研究管段的出入口溫度、流速和管壁溫度,利用上述關系就可以得到污垢熱阻。相應地,電磁場的抑垢率據化工行業標準得到[26]。

抑垢率公式如下

2.3 實驗條件

根據實驗臺現有條件,參考有關文獻對磁記憶時間的描述,初步設定處理腔長度為25 cm。結合線圈匝數和信號頻率對水體感生電流做功以及磁感應強度的影響,本應選擇大線徑多匝的纏繞方式,但是考慮到線圈過流發熱,在保證磁感應強度和感生電流盡量大的前提下,本實驗采用細線徑多匝的方案進行實驗研究。實驗介質為去離子水與無水Na2CO3、CaCl2按摩爾比1:1配制的溶液,鈣硬度為1000mg·L-1。

實驗中,下水箱的循環介質溫度控制在27℃,循環介質的流速設置為0.4m·s-1,恒溫水浴溫度控制在50℃。將銅質漆包線纏繞在加磁管道B管入口管段,呈螺旋管狀,線圈兩端接上自制變頻電磁裝置,實驗中裝置輸出波形為方波信號,輸出電流和頻率大小根據具體實驗進行調整。實驗周期為7~10 d,在線監測平臺每分鐘自動采集1次實驗參數,在線自動計算污垢熱阻并將數據存盤。

3 頻率對抑垢效果影響的實驗結果與分析

基于電磁場頻率對水處理系統的感生電流做功以及磁感應強度的影響,通過對線圈的電阻和電感的測量,計算感生電流做功的分界頻率為1.108 kHz。所以,基于電磁抑垢效果在線評價技術,在保證線圈所加交變方波電壓峰值不變的條件下,調節信號發生裝置的信號頻率在0~10 kHz范圍內以0.5 kHz的步長遞增,使水處理系統在不同頻率下工作,進而得到不同頻率電磁場條件下的污垢熱阻變化曲線,如圖9所示。一方面,通過曲線的變化可以觀察不同頻率電磁場條件下污垢誘導期(污垢誘導期指的是從溶液與表面接觸的時刻起到檢測到析晶沉積物的這段時間),來間接反映電磁頻率對污垢形成的抑制作用;另一方面也可以通過實驗數據計算不同頻率電磁場條件下抑垢率的變化,如圖10所示,直觀地反映了電磁頻率對污垢形成的影響。

圖9 不同頻率電磁場對污垢熱阻的影響

圖10 不同頻率電磁場作用下的抑垢率

污垢誘導期越長,說明模擬換熱器表面生成的垢越少,加磁處理具有一定的抑垢效果。從圖9不同頻率電磁場條件下的污垢熱阻變化曲線可以觀察到頻率為1 kHz以及0.5 kHz下加磁處理后換熱器表面污垢的誘導期到實驗結束時還未結束,頻率為5 kHz和10 kHz下加磁處理后換熱器表面污垢誘導期分別為40 h和60 h左右。從以上污垢誘導期的長短可以看出,加磁處理下的污垢誘導期比未加磁處理長,在頻率為1 kHz、0.5 kHz加磁處理下的污垢誘導期較其他頻率下長。大量的實驗研究證明,在1.108 kHz分界頻率附近的頻率點所得的模擬換熱器表面的污垢誘導期都在100 h以上,抑垢效果最明顯,證明了當電磁頻率在分界頻率附近時抑垢效果最明顯的理論。

從圖10不同頻率電磁場作用下的抑垢率曲線可以直觀地觀察到,1 kHz處的抑垢率最大,可達到90%左右,抑垢效果最好,也說明了當電磁頻率在分界頻率附近時抑垢效果最明顯。

4 結 論

(1)利用ANSYS仿真軟件對纏繞式變頻電磁水處理建立2D模型,并利用電路模型等效,計算出感生電流做功的表達式。當激勵信號電壓幅值一定時,處理腔內感應電流做功的大小與激勵信號頻率有關,并存在一個分界頻率,當激勵頻率在分界頻率以下時感應電流做功隨頻率的增大做功值大幅增加,在分界頻率以上感應電流做功變化甚微,此分界頻率大小在左右。

(2)由ANSYS仿真結果可知,纏繞式變頻電磁水處理器的線圈匝數一定時,處理腔中心磁感應強度隨激勵信號的頻率增加,磁感應強度呈現遞減趨勢,所以頻率在低頻段更適合,保持在20 kHz之內為優。

(3)本實驗工況下,分界頻率為1.108 kHz,當保證激勵信號電壓幅值一定時,改變激勵信號頻率在0~10 kHz內變化,動態實驗結果表明,在分界電磁頻率附近頻率點1 kHz時,模擬換熱面的污垢誘導期長達100 h,抑垢率也在90%左右,電磁抑垢效果理想,說明纏繞式變頻電磁水處理器在分界頻率附近工作時能達到較好的抑垢效果。

References

[1] Lipusa L C, A?ko B, Hamlerb A. Electromagnets for high-flow water processing [J]., 2011, 50 (9): 952-958

[2] Wang Jianguo (王建國), He Fang (何芳), Di Hao (邸昊). Correlation analysis of magnetic field and conductivity, pH value in electromagnetic restraint of scale formation [J].(化工學報), 2012, 63 (5): 1468-1473

[3] Shelver W L, Kamp L M, Church J L.Measurement of triclosan in water using a magnetic particle enzyme immunoassay [J]., 2007, 55 (10): 3758-3763

[4] Xiao Zengli (肖曾利), Pu Chunsheng (蒲春生). Review of technology and application on magnetized water treatment for anti-scaling [J].(斷塊油氣田), 2010, 17 (1): 121-125

[5] Chen Bowu (陳博武). Status and trends of the fouling monitoring technology of cooling water [J].(工業水處理), 2006, 26 (5): 4-8

[6] He Jun (何俊), Zhao Zongze (趙宗澤), Li Yuehua (李躍華). Research status and progress in scale inhibition and removal by physical technology [J].(工業水處理), 2010, 30 (9): 5-9

[7] Cheolho Bai. Effect of electronic anti-fouling treatment on fouling mitigation with circulating cooling-tower water [J]., 2001, 28 (5): 671-680

[8] Chang M C, Tai C Y. Effect of the magnetic field on the growth rate of aragonite and the precipitation of CaCO3[J]., 2010, 164 (1): 1-9

[9] Zhang Xuefeng (張雪峰), Niu Yonghong (牛永紅). Research progress on application of magnetic treatment industrial water technology [J].(化工進展), 2011, 30: 831-833

[10] Alimi Fathi, Tlili Mohamed, Gabrielli Claude. Effect of a magnetic water treatment on homogeneous and heterogeneous precipitation of calcium carbonate [J]., 2006, 40: 1941-1950

[11] Botello-Zubiate M E, Alvarez A, Martinez-Villafane A. Influence of magnetic water treatment on the calcium carbonate phase formation and the electrochemical corrosion behavior of carbon steel [J]., 2004, 369 (1/2): 256-259

[12] Song Liu. A study of fouling in a heat exchanger with an application of an electronic anti-fouling technology[D]. Michigan USA: Drexel University, 1999

[13] Jiang Dening (姜德寧), Zhou Kun (周昆). Experimental study on anti-scale of pipes based on alternating magnetic field [J].(天津工程師范學院學報), 2008, 18 (1): 21-25

[14] Jin Zhenhua (金貞花), Chen Yongchang (陳永昌), Xing Xiaokai (邢曉凱). Investigation on the high voltage electrostatic anti-fouling technology for heat transfer enhancement [J].(工程熱物理學報), 2006, 27 (6): 1035-1037

[15] Chen Can (陳璨), Liu Runhua (劉潤華). The impact of electromagnetic frequency on CaCO3solubility quantitative research [J].(科學技術與工程), 2009, 9 (15): 4567-4570

[16] Xing Xiaokai(邢曉凱), Ma Chongfang (馬重芳), ChenYongchang(陳永昌). Experimental study of electromagnetic anti-fouling technology to enhanced heat transfer [J].(石油大學學報), 2005, 29 (1): 79-83

[17] Fei Jiyou (費繼友), Li Yuquan (李玉泉), Bai Xin (白鑫). The electromagnetic frequency conversion water anti-scale technology and experimental research [J].(化工自動化及儀表), 2010, 38 (2): 157-161

[18] Chen Xiaozhuan (陳小磚), Liu Jianhua (柳建華), Ma Chongfang (馬重芳). Experimental research on mechanism of electromagnetic anti-fouling technology on cooling water flowing through absorber [J].(流體機械), 2011, 11 (39): 67-70

[19] Sung H Lee. A study of physical water treatment technology to mitigate the mineral fouling in a heat exchanger[D]. USA: Thesis of Drexel University, 2002

[20] Ding Jie (丁潔),Dong Lingyun (董凌云). Magnetic treatment water anti-scale and the influence in the anti-scaling induction period//9th National Chemical Engineering Science and Technology[C]. Qingdao, China, 1998: 216-220

[21] Xie Qifen (謝綺芬). The anti-scale mechanism and design of magnetic water device [J].(水處理技術), 1985, 11 (2): 49-52

[22] Donsldson J D, Grimes S. Lifting the scales from our pipes [J]., 1988, 117 (18): 43-46

[23] Han Yong (韓勇), Zhao Yongping (趙永平), Chai Xin (柴鑫). Analytical modeling of induced current for DC pulse water processor system and optimizing of excitation coil [J].(電機與控制學報), 2011, 9 (15): 31-38

[24] Yan Zhaowen (閻照文). ANSYS10.0 Electromagnetic Analysis Technology and Engineering Example Explanation (ANSYS10.0電磁分析技術及實例詳解) [M]. Beijing: China Water Power Press, 2006

[25] Yang Shanrang (楊善讓), Xu Zhiming (徐志明), Sun Lingfang (孫靈芳). Heat Exchanger Fouling and Countermeasures (換熱設備污垢與對策) [M]. 2nd ed. Beijing: Science Press, 2004: 1-5

[26] HG/T 2160—1991: Cooling water dynamic simulation test method[S]. Beijing: Chemical Industry Press, 1996

Influence of electromagnetic frequency on scale inhibition for spiral winding variable frequencyelectromagnetic water processor

WANG Jianguo, LI Yutong, DENG Lijuan

School of Automation EngineeringNortheast Dianli UniversityJilinJilinChina

For the spiral winding variable frequency electromagnetic water processor, magnetic induction intensity of water in processing cavity and induced current were analyzed with ANSYS simulation software. The inner relationship of magnetic induction intensity and induced current work with excitation signal frequency was studied. Using an online evaluation test bench for electromagnetic scale inhibition, the influence of spiral winding processing cavity frequency on scale inhibition of hard water circulation heat transfer experiment pipe was observed, and the ideal scale inhibition frequency was determined. Under this experimental condition, the frequency near the line frequency had longer scale induction period and stronger scale inhibition effect, and anti-scaling efficiency was the best.

spiral windingprocessing cavity; frequency; magnetic induction intensity; induced current; anti-scaling effect; convection; heat transfer; numerical analysis

2014-09-16.

LI Yutong, 302533570@qq.com

10.11949/j.issn.0438-1157.20141385

TQ 085

A

0438—1157(2015)03—0972—07

國家自然科學基金項目(51176028);吉林省重點科技 攻關項目(20140204006SF)。

2014-09-16收到初稿,2014-11-28收到修改稿。

聯系人:李雨通。第一作者:王建國(1963—),男,教授。

supported by the National Natural Science Foundation of China (51176028) and the Key Scientific and Technological Project of Jilin Province of China (20140204006SF).

猜你喜歡
信號實驗
記一次有趣的實驗
微型實驗里看“燃燒”
信號
鴨綠江(2021年35期)2021-04-19 12:24:18
完形填空二則
做個怪怪長實驗
孩子停止長個的信號
NO與NO2相互轉化實驗的改進
實踐十號上的19項實驗
太空探索(2016年5期)2016-07-12 15:17:55
基于LabVIEW的力加載信號采集與PID控制
一種基于極大似然估計的信號盲抽取算法
主站蜘蛛池模板: 国产一区二区三区夜色| 999国产精品永久免费视频精品久久| 免费黄色国产视频| 福利在线不卡| 99精品福利视频| 中字无码精油按摩中出视频| 波多野结衣一区二区三区四区| 国产一级毛片高清完整视频版| 噜噜噜久久| 午夜精品久久久久久久无码软件| 欧美狠狠干| 狠狠久久综合伊人不卡| 国产精品99久久久| 成人年鲁鲁在线观看视频| 伊人久久精品亚洲午夜| 亚洲高清中文字幕在线看不卡| 欧美精品啪啪一区二区三区| 国产一级毛片网站| 久久精品这里只有精99品| 激情亚洲天堂| 毛片在线看网站| 26uuu国产精品视频| 欧亚日韩Av| 秋霞午夜国产精品成人片| 国产成a人片在线播放| 制服丝袜 91视频| 狠狠v日韩v欧美v| 欧美福利在线| 亚洲综合精品第一页| 99热亚洲精品6码| 精品福利网| 国产精品午夜电影| 国产精品免费入口视频| 伊大人香蕉久久网欧美| 欧美亚洲一二三区| 中国美女**毛片录像在线| 国产69精品久久久久孕妇大杂乱 | 欧美人人干| 亚洲成a人片在线观看88| 国产精品视频免费网站| 97在线观看视频免费| 久久www视频| 亚洲狠狠婷婷综合久久久久| 日本人妻丰满熟妇区| 99热这里只有免费国产精品| 88av在线播放| 成人va亚洲va欧美天堂| 欧美日在线观看| 日韩AV手机在线观看蜜芽| 午夜毛片免费看| 精品久久综合1区2区3区激情| 国产特级毛片| 呦女亚洲一区精品| 99免费视频观看| 日韩av电影一区二区三区四区| 一级毛片免费不卡在线| 国产女人在线视频| 成人国产精品一级毛片天堂| 亚洲无码高清视频在线观看 | 免费在线视频a| 国产欧美视频在线| 日韩A∨精品日韩精品无码| 日韩免费毛片| 狠狠色噜噜狠狠狠狠奇米777| 欧美日韩免费在线视频| 四虎影视8848永久精品| 国产精品毛片一区| 国产精品无码AV片在线观看播放| 人妻无码AⅤ中文字| 国产精品自在在线午夜区app| 国内精品手机在线观看视频| 精品一區二區久久久久久久網站| 久久婷婷六月| 国产区福利小视频在线观看尤物| 综合色88| 亚洲大学生视频在线播放| 国产呦精品一区二区三区网站| 四虎永久免费地址| 啦啦啦网站在线观看a毛片| 亚洲美女操| 久久国产精品嫖妓| 亚洲综合在线网|