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

    基于正交試驗(yàn)的翼片式斜板沉沙池截除率研究

    2020-12-04 03:35:26趙法鑫楊力行鄭祖國(guó)
    關(guān)鍵詞:沉沙池斜板傾角

    趙法鑫,李 琳,楊力行,鄭祖國(guó)

    (1.新疆農(nóng)業(yè)大學(xué) 水利與土木工程學(xué)院,新疆 烏魯木齊 830052;2.浙江水利水電學(xué)院 水利與環(huán)境工程學(xué)院,浙江 杭州 310018)

    1 研究背景

    翼片式斜板沉沙池是是由日本學(xué)者丹保憲仁在一系列試驗(yàn)后提出的一種新型水沙分離技術(shù),又被稱為迷宮式沉沙池,其結(jié)構(gòu)是在斜板上布置與水流方向垂直的翼片(如圖1所示),水流流經(jīng)翼片形成三個(gè)區(qū)域:主流區(qū)、渦流區(qū)、環(huán)流區(qū)。泥沙顆粒在渦流區(qū)的作用下被帶入翼片槽內(nèi),槽內(nèi)環(huán)流區(qū)流速遠(yuǎn)小于主流區(qū),為泥沙沉降提供了良好的水力條件。在渦流和環(huán)流的共同作用下實(shí)現(xiàn)水沙的動(dòng)態(tài)分離。李繼震[1]研究了橫向流翼片式沉淀池內(nèi)沿深度方向濁度的變化規(guī)律,結(jié)合水平方向除濁公式,提出了包含水平、豎直兩個(gè)方向的除濁數(shù)學(xué)模式。He Cheng等[2]通過(guò)CFD數(shù)值模擬對(duì)翼片式斜板沉沙池內(nèi)的流場(chǎng)進(jìn)行了研究,分析發(fā)現(xiàn)水流在流經(jīng)翼片后形成的渦流區(qū)可將主流區(qū)的泥沙顆粒送入環(huán)流區(qū),泥沙顆粒在環(huán)流區(qū)內(nèi)增加了碰撞幾率有助于沉降,并通過(guò)模型試驗(yàn)驗(yàn)證了在沉沙池內(nèi)布置翼片式斜板可有效的提高截除率。胥標(biāo)[3]采用湍流雷諾模型及離散相模型對(duì)翼片式斜板沉沙池內(nèi)水沙兩相流進(jìn)行了模擬,分別就翼片格數(shù)、翼片寬度、斜板寬度和斜板傾角各單因素對(duì)截除率的影響做出了研究,得到了沉降效果最優(yōu)時(shí)的結(jié)構(gòu)參數(shù)。倪建華[4]將翼片式斜板應(yīng)用于含沙量較高的長(zhǎng)江以及嘉陵江的取水頭部,對(duì)粗沙細(xì)沙和粗細(xì)沙在不同進(jìn)水流速下進(jìn)行除沙效率及結(jié)構(gòu)參數(shù)優(yōu)化試驗(yàn),得出對(duì)除沙效率影響較高的因素為斷面流速、翼片斜板尺寸、傾角和翼片格數(shù)等。趙法鑫等[5]通過(guò)對(duì)三種不同沉沙池的各種表面負(fù)荷下的截除率進(jìn)行試驗(yàn),證明了翼片式斜板沉沙池截除率的優(yōu)越性,且在翼片式斜板沉沙池中已沉降的泥沙不易受水流影響再次懸浮流出沉沙池。陳濱[6]在鎮(zhèn)江市金山水廠翼片斜板沉淀池內(nèi)進(jìn)行原型試驗(yàn),分析了進(jìn)水均勻性、絮凝效果和翼片格數(shù)對(duì)去濁效果的影響程度。

    圖1 翼片式斜板流區(qū)示意圖

    前人通過(guò)試驗(yàn)和數(shù)值模擬的方法揭示了翼片式斜板沉沙池的泥沙沉降機(jī)理,對(duì)相關(guān)結(jié)構(gòu)參數(shù)如翼片高度、翼片間距、翼片傾角等開展了單因素試驗(yàn),定性分析了各因素對(duì)截除率的影響。由于影響截除率的因素較多,且各個(gè)影響因素對(duì)截除率的敏感性和顯著性尚缺乏相應(yīng)研究成果,實(shí)際工程應(yīng)用時(shí)翼片式斜板沉沙池的結(jié)構(gòu)優(yōu)化過(guò)程缺乏依據(jù)。同時(shí),由于前人對(duì)翼片式斜板沉沙池的截除率與各結(jié)構(gòu)參數(shù)最優(yōu)組合間的相關(guān)關(guān)系尚未開展系統(tǒng)研究,使其應(yīng)用于工程實(shí)際時(shí)缺乏設(shè)計(jì)依據(jù)。因此,本文擬基于正交試驗(yàn),通過(guò)極差分析得出影響截除率的敏感因素,借助方差分析法判斷各因素對(duì)截除率影響的顯著性。采用投影尋蹤回歸(PPR)方法建立穩(wěn)定可靠的截除率與表面負(fù)荷、翼片高度及翼片間距的預(yù)測(cè)模型,獲得各因素對(duì)截除率的貢獻(xiàn)權(quán)重,和不同表面負(fù)荷時(shí)的截除率最優(yōu)時(shí)的結(jié)構(gòu)參數(shù)組合。研究成果可望為保證較高截除率時(shí)的翼片式沉沙池結(jié)構(gòu)設(shè)計(jì)和優(yōu)化提供依據(jù)和方法。

    2 截除率正交試驗(yàn)簡(jiǎn)介

    2.1 試驗(yàn)裝置為了方便觀察試驗(yàn)過(guò)程中的水流與顆粒運(yùn)動(dòng)的運(yùn)動(dòng)情況,本試驗(yàn)沉沙池采用透明有機(jī)玻璃制成,翼片式斜板采用黑色PVC塑料板制成。沉沙池主要尺寸及構(gòu)成如圖2所示。應(yīng)用沉沙池池寬(D)對(duì)翼片間距和翼片高度進(jìn)行無(wú)量綱化處理,試驗(yàn)中翼片間距分別為D/3、D/2、D/3,翼片高度分別為D/3、D/2、2D/3。翼片傾角分別取90°、75°、60°。

    圖2 翼片式斜板沉沙池模型系統(tǒng)試驗(yàn)布置圖(單位:mm)

    2.2 試驗(yàn)材料選擇試驗(yàn)選用核桃沙作為模型沙進(jìn)行試驗(yàn)。核桃沙粒徑范圍0.8~1.2 mm,堆積密度為0.81 t/m3。核桃沙呈黃褐色,能清晰地觀察到其在沉沙池內(nèi)及翼片槽內(nèi)的運(yùn)動(dòng)軌跡,且易于被攔截不會(huì)再次進(jìn)入循環(huán)系統(tǒng)。前人研究結(jié)果表明核桃沙屬于輕質(zhì)沙,其沉降特性與天然沙接近,使用核桃沙進(jìn)行試驗(yàn)可以獲得可靠的試驗(yàn)結(jié)果[7-8]。

    2.3 試驗(yàn)方案試驗(yàn)選取翼片間距(A)、翼片傾角(B)、翼片高度(C)三個(gè)影響因素,每個(gè)因素下設(shè)計(jì)三個(gè)水平變量。本次試驗(yàn)所選因素及水平設(shè)計(jì)如表1所示。試驗(yàn)選用了三因素三水平試驗(yàn)方案,故正交表應(yīng)選擇L9(34),并將最后一列設(shè)為空列即誤差列,以便方差分析。其試驗(yàn)方案見(jiàn)表2。

    表1 試驗(yàn)因素及水平設(shè)計(jì)

    表2 正交試驗(yàn)設(shè)計(jì)方案

    筆者通過(guò)前期試驗(yàn)發(fā)現(xiàn),在表面負(fù)荷小于16 m3/(h·m2)時(shí)翼片式斜板沉沙池截除率相比普通斜板沉沙池優(yōu)勢(shì)并不明顯,且采用翼片式斜板反而會(huì)增加建設(shè)成本,故本試驗(yàn)選取輸沙率為0.455 g/s、表面負(fù)荷分別為:16 m3/(h·m2),20 m3/(h·m2),24 m3/(h·m2),28 m3/(h·m2),32 m3/(h·m2)情況下進(jìn)行截除率試驗(yàn)。

    2.4 試驗(yàn)步驟圖2為翼片式斜板沉沙池模型系統(tǒng)試驗(yàn)布置圖。試驗(yàn)開始時(shí)從進(jìn)水槽上方投放模型沙,在進(jìn)流紊動(dòng)作用下核桃沙與槽內(nèi)水流充分混合,模型沙在沉沙池內(nèi)以懸移運(yùn)動(dòng)為主,經(jīng)翼片式斜板區(qū)域后,由出水擋板頂溢流入集水槽再經(jīng)出水口流出。出水口下方放置回水箱,回水箱上方布置篩網(wǎng)(篩網(wǎng)孔徑100目)攔截隨水流溢出的模型沙,保證隨水流溢出的模型沙不會(huì)再次進(jìn)入沉沙池中?;厮?jīng)進(jìn)水管再次進(jìn)入進(jìn)水槽,實(shí)現(xiàn)整個(gè)系統(tǒng)的自循環(huán)。

    3 試驗(yàn)結(jié)果分析

    3.1 截除率的各影響因素敏感性與顯著性分析對(duì)表2所列各方案下的翼片式斜板沉沙池截除率進(jìn)行試驗(yàn),每組試驗(yàn)結(jié)束后收集篩網(wǎng)上的核桃沙,用烘干機(jī)烘干后使用電子秤(精度為0.1 g)稱其質(zhì)量記為M1,加沙總質(zhì)量記為M2,(M2-M1)/M2記為各組沉沙池的截除率。對(duì)試驗(yàn)所選用的各表面負(fù)荷下正交試驗(yàn)結(jié)果進(jìn)行極差分析。因素的極差大小反應(yīng)了該因素水平變動(dòng)時(shí),試驗(yàn)指標(biāo)即截除率的變動(dòng)幅度。極差越大說(shuō)明該因素對(duì)截除率的影響也越大,故可以獲得各因素敏感性及最優(yōu)參數(shù)組合,結(jié)果見(jiàn)表3,表中Kij表示第i因素第j水平所對(duì)應(yīng)的試驗(yàn)指標(biāo)和的平均。從表中可以看出,各試驗(yàn)工況下空列的極差值均小于其他因素的極差,表明在試驗(yàn)過(guò)程中試驗(yàn)誤差對(duì)截除率影響小于其他三項(xiàng),即該組試驗(yàn)數(shù)據(jù)有效,可進(jìn)行方差分析,繼而獲得翼片間距、傾角、高度對(duì)截除率的影響。

    由表3還可知,在表面負(fù)荷16 m3/(h·m2)、20 m3/(h·m2)、24 m3/(h·m2)時(shí)三者極差分析結(jié)果一致。表面負(fù)荷一定的情況下,因素翼片間距(A)極差最大,其次是翼片高度(C),最后是翼片傾角(B)。截除率對(duì)各因素敏感性由大到小為A-C-B,即表面負(fù)荷16 m3/(h·m2)、20 m3/(h·m2)、24 m3/(h·m2)時(shí)翼片間距(A)對(duì)截除率起著主要影響作用,翼片式斜板最優(yōu)結(jié)構(gòu)組合為翼片間距取D/2,翼片傾角取60°,翼片高度取D/2。

    表面負(fù)荷28m3/(h·m2)及32 m3/(h·m2)時(shí)兩者極差分析結(jié)果一致。因素翼片高度(C)極差最大,其次是翼片間距(A),最后是翼片傾角(B)。截除率對(duì)各因素敏感性由大到小為C-A-B。翼片式斜板最優(yōu)結(jié)構(gòu)組合為翼片間距取D/3,翼片傾角取90°,翼片高度D/2。

    為了更加直觀的分析各因素對(duì)截除率的影響,根據(jù)表3可做出各因素對(duì)截除率影響的直觀分析圖,如圖3所示。從圖中可以看出,表面負(fù)荷16 m3/(h·m2)、20 m3/(h·m2)、24 m3/(h·m2)時(shí)截除率隨著翼片間距、翼片高度的增大先增大后減小,隨著翼片傾角傾斜程度的增加而增加,且隨著翼片間距的變化最為明顯,即截除率對(duì)翼片間距最為敏感。28 m3/(h·m2)及32 m3/(h·m2)時(shí),截除率隨著翼片間距的增大而減小,隨著翼片夾角傾角的變大而減小,隨著翼片高度的增加截除率呈現(xiàn)先增大后急劇變小的趨勢(shì)。翼片高度對(duì)截除率的變化影響最大,即截除率對(duì)翼片高度最為敏感。

    表3 截除率極差分析結(jié)果

    圖3 截除率影響因素直觀圖

    為了得到各因素對(duì)截除率的影響程度是否顯著,對(duì)正交試驗(yàn)結(jié)果進(jìn)行方差分析,在方差分析中將因素的顯著性水平劃分為四個(gè)等級(jí)(1)當(dāng)F≥F0.01為高度顯著影響,記為“***”,(2)當(dāng)F0.01>F≥F0.05為顯著影響,記為“**”,(3)當(dāng)F0.05>F≥F0.1為影響不十分顯著,記為“*”,(4)當(dāng)F<F0.1為無(wú)顯著影響,記為“×”。分析結(jié)果見(jiàn)表4。由表4可以看出表面負(fù)荷16 m3/(h·m2)時(shí)翼片間距為顯著影響因素,傾角、高度為無(wú)顯著影響因素。20 m3/(h·m2)時(shí)翼片間距為顯著影響因素,高度為不十分顯著因素,傾角為無(wú)顯著影響因素。24 m3/(h·m2)時(shí)翼片間距為不十分顯著因素,傾角、高度為無(wú)顯著影響因素。28 m3/(h·m2)時(shí)翼片高度為顯著影響因素,間距、傾角為無(wú)顯著影響因素。32 m3/(h·m2)時(shí)翼片高度為高度顯著影響因素,間距為不十分顯著因素,傾角為無(wú)顯著影響因素。

    3.2 基于投影尋蹤回歸模型(PPR)的截除率仿真計(jì)算及分析投影尋蹤回歸的原理是將試驗(yàn)得到的高維數(shù)據(jù)利用計(jì)算機(jī)將其投影到低維子空間上,借此實(shí)現(xiàn)“高維數(shù)據(jù)低維化”的目的,通過(guò)對(duì)低維子空間數(shù)據(jù)的分析,找出能夠反映高維原始數(shù)據(jù)內(nèi)部特征的非正態(tài)信息和規(guī)律并進(jìn)行回歸分析,從而實(shí)現(xiàn)對(duì)高維數(shù)據(jù)的分析[9]。PPR模型表達(dá)方式及詳細(xì)計(jì)算方法可見(jiàn)參考文獻(xiàn)[10-11],此處不再贅述。

    通過(guò)上文正交試驗(yàn)結(jié)果極差分析和方差分析可以看出,影響截除率的主要因素為翼片間距與翼片高度,翼片傾角是非顯著性因素,為提高運(yùn)算效率及精度在使用PPR建模時(shí)不考慮翼片傾角,而是將翼片間距、翼片高度、表面負(fù)荷作為建模參數(shù)。除45組建模樣本外,事先預(yù)留了18組檢驗(yàn)樣本。穩(wěn)定性分析時(shí)首先采用PPR程序?qū)颖具M(jìn)行建模運(yùn)算,并基于此模型對(duì)建模樣本進(jìn)行“還原擬合”計(jì)算,取“還原擬合”值與建模樣本實(shí)測(cè)值做差,一般相對(duì)誤差|δ|<10%作為合格率判斷標(biāo)準(zhǔn)[12],借此確定模型的合格率。之后將預(yù)留檢驗(yàn)樣本數(shù)據(jù)帶入PPR模型,計(jì)算預(yù)留檢驗(yàn)樣本合格率,當(dāng)其也滿足合格標(biāo)準(zhǔn)時(shí),認(rèn)為此PPR模型具有較好的穩(wěn)定性與可靠度。

    采用鄭祖國(guó)等[13-14]自主開發(fā)的PPR程序?qū)σ砥叫卑宄辽吵亟爻蔬M(jìn)行建模分析。運(yùn)行時(shí)指定5個(gè)投影參數(shù)分別為:P=3,Q=1,S=0.1,M=5,Mu=3。其中P為影響因素個(gè)數(shù);Q為因變量個(gè)數(shù);S為光滑系數(shù),決定了模型的靈敏度,S越小則模型越靈敏,取值范圍為0<S<1;M為嶺函數(shù)上限個(gè)數(shù);Mu為嶺函數(shù)最優(yōu)個(gè)數(shù)。M與Mu決定了模型尋找數(shù)據(jù)內(nèi)在結(jié)構(gòu)的精密程度。通過(guò)PPR建模運(yùn)算可得到模型的投影方向如式(1)所列,模型所需嶺函數(shù)如圖4所示。同時(shí)可得到相應(yīng)嶺函數(shù)的權(quán)重系數(shù)β如式(2)所示。

    表4 截除率方差分析結(jié)果

    圖4 翼片式斜板沉沙池截除率PPR模型嶺函數(shù)圖

    此外,用PPR程序?qū)σ砥叫卑宄辽吵亟爻蔬M(jìn)行計(jì)算,得到各影響因子對(duì)翼片式沉沙池截除率貢獻(xiàn)權(quán)重系數(shù)為:表面負(fù)荷=1,翼片間距=0.686,翼片高度=0.527。

    依據(jù)上述PPR建模過(guò)程,對(duì)截除率建模樣本與其預(yù)留檢驗(yàn)樣本進(jìn)行預(yù)測(cè),并與實(shí)測(cè)值進(jìn)行比較,如圖5所示。由圖5可知,建模樣本與檢驗(yàn)樣本的截除率預(yù)測(cè)值與實(shí)測(cè)值偏差較小,當(dāng)相對(duì)誤差|δ|<6%時(shí),建模樣本與檢驗(yàn)樣本的合格率均為100%,當(dāng)相對(duì)誤差|δ|<5%時(shí)建模樣本與檢驗(yàn)樣本的合格率分別為93.33%和94.44%,平均相對(duì)誤差分別為2.03%和1.92%。表明PPR模型具有較好的可靠性,故可根據(jù)此模型對(duì)未進(jìn)行試驗(yàn)的工況進(jìn)行預(yù)測(cè),分析各種工況下翼片式斜板沉沙池的最優(yōu)結(jié)構(gòu)參數(shù)組合,為實(shí)際應(yīng)用提供參考。

    例如,借助此模型分別對(duì)表面負(fù)荷取16,17,…,31,32 m3(/h·m2);翼片間距分別取A/D為0.333,0.417,0.5,0.583,0.667;翼片高度分別取C/D為0.33,0.417,0.5,0.5833,0.667時(shí)的不同試驗(yàn)組合進(jìn)行仿真運(yùn)算。因數(shù)據(jù)較多,以表面負(fù)荷22 m3(/h·m2)與30 m3(/h·m2)為例,繪制翼片間距(A)/池寬(D)與翼片高度(C)/池寬(D)對(duì)截除率影響的等值線圖即圖6。通過(guò)等值線圖可以直觀的觀察到在某一表面負(fù)荷下翼片式斜板最優(yōu)的間距及高度組合情況。由圖6可以看出,不同表面負(fù)荷下截除率都隨著翼片間距占比的增加呈現(xiàn)先增加后減小的趨勢(shì),最優(yōu)區(qū)間在翼片間距占池寬0.5~0.6附近,不同表面負(fù)荷下截除率也隨著翼片高度占比的增加呈現(xiàn)先增加后減小的趨勢(shì),最優(yōu)區(qū)間在翼片高度占池寬0.5附近。

    圖5 截除率PPR建模計(jì)算精度分析

    圖6 截除率與A/D、C/D的關(guān)系

    4 結(jié)論

    選用模型沙(核桃沙),基于正交試驗(yàn)在多種表面負(fù)荷下對(duì)不同結(jié)構(gòu)參數(shù)的翼片式斜板沉沙池的截除率進(jìn)行了系統(tǒng)試驗(yàn),采用極差分析法與方差分析法對(duì)試驗(yàn)結(jié)果進(jìn)行敏感性與顯著性分析,并建立了有較好可靠性的PPR運(yùn)算模型,借助此模型對(duì)未能進(jìn)行試驗(yàn)的工況進(jìn)行預(yù)測(cè),主要得到以下結(jié)論:

    (1)不同表面負(fù)荷下,截除率的敏感因素與顯著性因素是不同的。在表面負(fù)荷小于等于24 m3/(h·m2)時(shí),截除率的敏感因素是翼片間距,主要顯著性因素也為翼片間距。當(dāng)表面負(fù)荷大于24 m3/(h·m2)時(shí),截除率的敏感性因素為翼片高度,主要顯著性因素也為翼片高度。翼片傾角不是截除率的敏感因素與顯著性因素。

    (2)結(jié)合敏感性分析、顯著性分析及PPR仿真計(jì)算分析,可知在表面負(fù)荷16~32 m3/(h·m2)范圍內(nèi)時(shí),翼片間距與池寬的比值取0.35~0.55時(shí)截除率較優(yōu),翼片高度與池寬的比值取0.45~0.55時(shí)截除率較優(yōu),翼片傾角在表面負(fù)荷小于等于24 m3/(h·m2)時(shí)可取60°,表面負(fù)荷大于24 m3/(h·m2)可取90°。

    (3)將PPR程序用于翼片式斜板沉沙池截除率高維建模中,建立了穩(wěn)定可靠的計(jì)算模型,能夠客觀的挖掘翼片式斜板沉沙池截除率高維數(shù)據(jù)內(nèi)在的本質(zhì)結(jié)構(gòu),確定各影響因子對(duì)截除率貢獻(xiàn)權(quán)重系數(shù)由大到小為:表面負(fù)荷,翼片間距,翼片高度。所采用的方法能對(duì)不同結(jié)構(gòu)參數(shù)的翼片式斜板截除率進(jìn)行預(yù)測(cè)且有著較好的可靠性,可以給出不同表面負(fù)荷下截除率最優(yōu)結(jié)構(gòu)參數(shù)建議,為生產(chǎn)設(shè)計(jì)提供理論支撐。

    猜你喜歡
    沉沙池斜板傾角
    一種用于塑料粒子均勻冷卻的冷卻裝置
    淺談水廠浮沉池斜板更新改造
    地球軸傾角的改斜歸正
    激光傾角儀在CT引導(dǎo)下經(jīng)皮肺穿刺活檢中的應(yīng)用
    車輪外傾角和前束角匹配研究
    北京汽車(2021年3期)2021-07-17 08:25:42
    斜板沉淀池在煉鋼污水處理中的應(yīng)用分析
    滲透破壞沖淤式沉沙池技術(shù)研究
    引黃濟(jì)青沉沙池運(yùn)行存在問(wèn)題及對(duì)策
    山東水利(2017年5期)2017-04-07 06:29:51
    自排式沉沙池上游渠道水力學(xué)特性研究
    新型翼片式斜板沉淀池的數(shù)值模擬與優(yōu)化研究
    人人妻人人澡人人爽人人夜夜 | 亚洲成人中文字幕在线播放| 中文字幕免费在线视频6| 久久久精品欧美日韩精品| 男的添女的下面高潮视频| 亚洲成人av在线免费| 精品不卡国产一区二区三区| 亚洲三级黄色毛片| a级毛色黄片| 国产精品久久久久久av不卡| 亚洲精品456在线播放app| 亚洲精品影视一区二区三区av| 国内揄拍国产精品人妻在线| 欧美变态另类bdsm刘玥| 国产大屁股一区二区在线视频| 亚洲图色成人| 91久久精品国产一区二区成人| 亚洲av福利一区| av播播在线观看一区| av免费观看日本| 99视频精品全部免费 在线| 可以在线观看毛片的网站| 大又大粗又爽又黄少妇毛片口| 国产精品伦人一区二区| 性色avwww在线观看| av一本久久久久| 少妇丰满av| 一级毛片黄色毛片免费观看视频| 丝袜喷水一区| 精品一区二区三区视频在线| 欧美性感艳星| 寂寞人妻少妇视频99o| 在线观看人妻少妇| 国产v大片淫在线免费观看| 最近中文字幕高清免费大全6| 日韩制服骚丝袜av| 两个人的视频大全免费| 国产又色又爽无遮挡免| 国产精品美女特级片免费视频播放器| 2018国产大陆天天弄谢| 97人妻精品一区二区三区麻豆| 免费大片黄手机在线观看| 秋霞伦理黄片| 亚洲精品国产av蜜桃| 非洲黑人性xxxx精品又粗又长| 国产免费视频播放在线视频 | 欧美3d第一页| 最后的刺客免费高清国语| 色哟哟·www| 18禁在线播放成人免费| 久久99热这里只频精品6学生| 久久久久久久久大av| 又粗又硬又长又爽又黄的视频| 97超视频在线观看视频| 免费观看a级毛片全部| 亚洲怡红院男人天堂| 亚洲av成人精品一二三区| 久久热精品热| 国产精品1区2区在线观看.| 成人美女网站在线观看视频| 午夜精品一区二区三区免费看| 亚洲成人一二三区av| 韩国av在线不卡| 亚洲综合精品二区| 日韩亚洲欧美综合| 国产在视频线精品| 国产亚洲一区二区精品| 尤物成人国产欧美一区二区三区| 国产在线男女| 少妇高潮的动态图| 全区人妻精品视频| 在线观看av片永久免费下载| 欧美区成人在线视频| h日本视频在线播放| 精品国产露脸久久av麻豆 | 国产精品久久久久久精品电影| 有码 亚洲区| 免费看美女性在线毛片视频| videossex国产| 激情 狠狠 欧美| 国产精品国产三级国产av玫瑰| 亚洲国产精品sss在线观看| 美女国产视频在线观看| 亚洲av福利一区| 一边亲一边摸免费视频| 一级毛片电影观看| 91久久精品国产一区二区成人| 婷婷色综合www| 午夜福利视频1000在线观看| 免费黄色在线免费观看| 日韩欧美国产在线观看| 亚洲不卡免费看| 亚洲av成人精品一区久久| 亚洲国产日韩欧美精品在线观看| 精品酒店卫生间| 国产男女超爽视频在线观看| 亚洲精品日韩在线中文字幕| 国内少妇人妻偷人精品xxx网站| 一区二区三区高清视频在线| 亚洲欧美精品专区久久| 蜜臀久久99精品久久宅男| 国产伦精品一区二区三区视频9| 国产成年人精品一区二区| 网址你懂的国产日韩在线| 免费看光身美女| 日韩欧美一区视频在线观看 | 99热全是精品| 日韩一区二区三区影片| 成年免费大片在线观看| 亚洲av电影不卡..在线观看| 欧美一级a爱片免费观看看| 色5月婷婷丁香| 欧美三级亚洲精品| 国产午夜精品久久久久久一区二区三区| 亚洲精品日韩av片在线观看| 精品一区二区免费观看| 久久久久精品久久久久真实原创| 亚洲自偷自拍三级| 狂野欧美激情性xxxx在线观看| 亚洲精品久久午夜乱码| 久久午夜福利片| 好男人在线观看高清免费视频| 欧美成人精品欧美一级黄| 久久草成人影院| 欧美高清性xxxxhd video| 乱码一卡2卡4卡精品| 久久精品国产亚洲av天美| 在线观看人妻少妇| 久久草成人影院| 人妻少妇偷人精品九色| 18禁裸乳无遮挡免费网站照片| 一区二区三区四区激情视频| 舔av片在线| 久久精品综合一区二区三区| 爱豆传媒免费全集在线观看| 97超视频在线观看视频| 亚洲av成人精品一二三区| 国产 一区精品| 亚洲av国产av综合av卡| 中文字幕av成人在线电影| 天堂中文最新版在线下载 | 校园人妻丝袜中文字幕| 亚洲综合色惰| 少妇被粗大猛烈的视频| 日韩成人伦理影院| av在线亚洲专区| 一级黄片播放器| 精品一区二区三区视频在线| 能在线免费看毛片的网站| 日日摸夜夜添夜夜爱| 91久久精品国产一区二区三区| 国产亚洲91精品色在线| av又黄又爽大尺度在线免费看| 日韩不卡一区二区三区视频在线| 亚洲成人中文字幕在线播放| 国产午夜精品久久久久久一区二区三区| 91aial.com中文字幕在线观看| 日韩精品青青久久久久久| 国产免费福利视频在线观看| 国产av码专区亚洲av| 日日啪夜夜撸| 纵有疾风起免费观看全集完整版 | 亚洲精品成人久久久久久| 3wmmmm亚洲av在线观看| 亚洲欧美精品专区久久| 午夜福利在线在线| 亚洲激情五月婷婷啪啪| 久久99蜜桃精品久久| 国产亚洲5aaaaa淫片| 日韩av在线免费看完整版不卡| 人妻夜夜爽99麻豆av| 中文天堂在线官网| 纵有疾风起免费观看全集完整版 | 久久精品久久久久久噜噜老黄| 亚洲伊人久久精品综合| 草草在线视频免费看| 日韩中字成人| 欧美日韩精品成人综合77777| 免费观看a级毛片全部| 中文欧美无线码| 国产av在哪里看| 国产三级在线视频| 五月伊人婷婷丁香| 免费观看a级毛片全部| 日本与韩国留学比较| 插逼视频在线观看| 亚洲最大成人中文| 观看免费一级毛片| 亚洲av中文字字幕乱码综合| 精品国产一区二区三区久久久樱花 | 性插视频无遮挡在线免费观看| 男插女下体视频免费在线播放| 美女高潮的动态| 天堂av国产一区二区熟女人妻| 成人亚洲精品一区在线观看 | 日日撸夜夜添| 美女xxoo啪啪120秒动态图| 精品人妻熟女av久视频| 国产精品蜜桃在线观看| 看黄色毛片网站| 夜夜看夜夜爽夜夜摸| 精品人妻视频免费看| 精品国产三级普通话版| 国产午夜精品久久久久久一区二区三区| 日韩欧美精品v在线| 亚洲一级一片aⅴ在线观看| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 免费观看精品视频网站| 赤兔流量卡办理| 特级一级黄色大片| 三级经典国产精品| 欧美三级亚洲精品| 久久久亚洲精品成人影院| 国产极品天堂在线| 日本色播在线视频| 国产一区亚洲一区在线观看| 中文精品一卡2卡3卡4更新| 国产精品日韩av在线免费观看| 免费观看性生交大片5| 国模一区二区三区四区视频| 2021少妇久久久久久久久久久| 国产精品一二三区在线看| 久热久热在线精品观看| 国产久久久一区二区三区| 80岁老熟妇乱子伦牲交| 高清视频免费观看一区二区 | 国产一区二区亚洲精品在线观看| 国产精品一区二区三区四区久久| 高清av免费在线| 免费看美女性在线毛片视频| 欧美zozozo另类| 国产精品熟女久久久久浪| 免费黄色在线免费观看| 男人狂女人下面高潮的视频| 亚洲人成网站在线播| 欧美97在线视频| 中文资源天堂在线| 欧美激情国产日韩精品一区| av.在线天堂| 国产高潮美女av| 国产 一区精品| 国产av不卡久久| 成人漫画全彩无遮挡| 18禁在线无遮挡免费观看视频| 好男人在线观看高清免费视频| 免费在线观看成人毛片| 成年av动漫网址| 日韩电影二区| 久久鲁丝午夜福利片| 亚洲电影在线观看av| 亚洲av日韩在线播放| 亚洲欧美精品专区久久| 99热网站在线观看| 一级毛片电影观看| av卡一久久| 亚洲欧美中文字幕日韩二区| 亚洲精品456在线播放app| 国产 一区 欧美 日韩| 美女被艹到高潮喷水动态| 久久久久久久久久成人| 九九在线视频观看精品| 婷婷色麻豆天堂久久| 日韩制服骚丝袜av| 最近视频中文字幕2019在线8| 性插视频无遮挡在线免费观看| 成人性生交大片免费视频hd| 精品久久久久久久久久久久久| 80岁老熟妇乱子伦牲交| 亚洲av二区三区四区| 亚洲真实伦在线观看| 26uuu在线亚洲综合色| 床上黄色一级片| 国产精品久久久久久av不卡| 午夜爱爱视频在线播放| 亚洲久久久久久中文字幕| 国产色婷婷99| 不卡视频在线观看欧美| 国产成人精品久久久久久| 美女大奶头视频| 精品久久久久久久久av| 日韩av在线大香蕉| 亚洲人与动物交配视频| 黄色欧美视频在线观看| 一个人看视频在线观看www免费| 大香蕉久久网| 汤姆久久久久久久影院中文字幕 | 精品久久久久久久久av| 亚洲aⅴ乱码一区二区在线播放| 男人爽女人下面视频在线观看| 18禁动态无遮挡网站| 国产色婷婷99| 亚洲最大成人中文| 永久免费av网站大全| 亚洲成人中文字幕在线播放| 亚洲国产欧美在线一区| 能在线免费观看的黄片| 国产欧美日韩精品一区二区| 国产黄频视频在线观看| 美女黄网站色视频| 两个人视频免费观看高清| 亚洲国产色片| 久久久久久久久久久丰满| 中文欧美无线码| 精品国产三级普通话版| 麻豆国产97在线/欧美| 国产精品久久视频播放| 国产v大片淫在线免费观看| 黄色一级大片看看| 免费播放大片免费观看视频在线观看| 国产成人a∨麻豆精品| 亚洲美女搞黄在线观看| 久久鲁丝午夜福利片| 久久国产乱子免费精品| 校园人妻丝袜中文字幕| 人妻系列 视频| 波多野结衣巨乳人妻| 亚洲av国产av综合av卡| 亚洲最大成人av| 免费在线观看成人毛片| 亚洲最大成人中文| 在线观看一区二区三区| 插阴视频在线观看视频| 丝袜喷水一区| ponron亚洲| 欧美精品国产亚洲| 非洲黑人性xxxx精品又粗又长| 国产高清不卡午夜福利| 97在线视频观看| 建设人人有责人人尽责人人享有的 | 亚洲最大成人手机在线| 久久精品国产鲁丝片午夜精品| 国产成人freesex在线| 中文欧美无线码| 激情五月婷婷亚洲| 人妻系列 视频| 成人午夜精彩视频在线观看| 午夜日本视频在线| 日韩欧美三级三区| av女优亚洲男人天堂| 黄色欧美视频在线观看| 人妻夜夜爽99麻豆av| 精品久久国产蜜桃| 少妇人妻精品综合一区二区| 国产精品麻豆人妻色哟哟久久 | 淫秽高清视频在线观看| 中文乱码字字幕精品一区二区三区 | 十八禁网站网址无遮挡 | 欧美xxⅹ黑人| 日韩不卡一区二区三区视频在线| 午夜亚洲福利在线播放| 亚洲精品成人av观看孕妇| av专区在线播放| 国产爱豆传媒在线观看| 亚洲av福利一区| 亚洲最大成人手机在线| 观看免费一级毛片| 亚洲成人av在线免费| 亚洲欧美成人综合另类久久久| 国产精品久久久久久av不卡| 亚洲国产成人一精品久久久| 搞女人的毛片| av在线观看视频网站免费| 国产精品爽爽va在线观看网站| 最近2019中文字幕mv第一页| 国产在线男女| 国产乱人视频| 久久久久久久午夜电影| 熟女人妻精品中文字幕| 啦啦啦韩国在线观看视频| 国产免费福利视频在线观看| 亚洲av成人av| 午夜精品在线福利| 成人美女网站在线观看视频| 国产精品久久久久久久电影| 国产精品久久久久久精品电影| 六月丁香七月| 国产成人精品一,二区| 精品熟女少妇av免费看| 99热网站在线观看| 一区二区三区乱码不卡18| 亚洲天堂国产精品一区在线| 国产有黄有色有爽视频| 国产精品国产三级国产av玫瑰| 18禁裸乳无遮挡免费网站照片| 国产亚洲av片在线观看秒播厂 | 3wmmmm亚洲av在线观看| 久久国内精品自在自线图片| 人人妻人人看人人澡| 26uuu在线亚洲综合色| 国产人妻一区二区三区在| 久久久久久久久中文| 我的女老师完整版在线观看| 免费黄色在线免费观看| 久久精品国产亚洲av涩爱| 国产成人福利小说| 亚洲人成网站在线播| 久久久久久九九精品二区国产| 精品一区在线观看国产| 秋霞在线观看毛片| 特级一级黄色大片| 少妇被粗大猛烈的视频| 在线观看免费高清a一片| av线在线观看网站| 色网站视频免费| 97超视频在线观看视频| 中文乱码字字幕精品一区二区三区 | 日韩成人伦理影院| 亚洲国产精品sss在线观看| 欧美zozozo另类| av黄色大香蕉| 永久免费av网站大全| 久久久精品欧美日韩精品| 久久精品久久久久久久性| 国产成人freesex在线| 亚洲欧洲国产日韩| 国产在视频线精品| 可以在线观看毛片的网站| 亚洲精品久久久久久婷婷小说| 一个人免费在线观看电影| 精品久久久久久成人av| 久久精品人妻少妇| 日韩成人av中文字幕在线观看| 国产在线一区二区三区精| 国产成人精品久久久久久| a级毛片免费高清观看在线播放| 麻豆精品久久久久久蜜桃| 成人高潮视频无遮挡免费网站| 国产极品天堂在线| 最近中文字幕高清免费大全6| 日本三级黄在线观看| 欧美精品一区二区大全| 99九九线精品视频在线观看视频| 国产精品蜜桃在线观看| 插逼视频在线观看| 97在线视频观看| 一级av片app| 免费观看av网站的网址| or卡值多少钱| 亚洲精品国产av蜜桃| 国产午夜福利久久久久久| 国产精品爽爽va在线观看网站| 亚洲av不卡在线观看| 免费观看av网站的网址| 激情五月婷婷亚洲| 久久久久久久国产电影| 天堂中文最新版在线下载 | 日本-黄色视频高清免费观看| 色综合站精品国产| 嘟嘟电影网在线观看| 永久网站在线| or卡值多少钱| 中文字幕人妻熟人妻熟丝袜美| 亚洲精品中文字幕在线视频 | 精品亚洲乱码少妇综合久久| 有码 亚洲区| 久久久久久久久久久免费av| 亚州av有码| 国产淫片久久久久久久久| 美女被艹到高潮喷水动态| 亚洲成人精品中文字幕电影| 亚洲国产最新在线播放| 成年av动漫网址| 人妻制服诱惑在线中文字幕| 国产精品福利在线免费观看| av免费在线看不卡| videossex国产| 欧美极品一区二区三区四区| 日本-黄色视频高清免费观看| 久久久久久国产a免费观看| 色视频www国产| 欧美人与善性xxx| av国产免费在线观看| 麻豆成人午夜福利视频| 九草在线视频观看| 又粗又硬又长又爽又黄的视频| 51国产日韩欧美| 亚洲电影在线观看av| 在线观看美女被高潮喷水网站| 国产亚洲av片在线观看秒播厂 | 91午夜精品亚洲一区二区三区| 在线观看美女被高潮喷水网站| 中文字幕久久专区| av国产免费在线观看| 三级毛片av免费| 国产精品一区二区在线观看99 | 天天躁日日操中文字幕| 少妇的逼水好多| 26uuu在线亚洲综合色| 美女内射精品一级片tv| 久久久久久久亚洲中文字幕| 亚洲激情五月婷婷啪啪| 国产综合精华液| 免费观看精品视频网站| 在线观看av片永久免费下载| 免费大片18禁| 国产 一区 欧美 日韩| 不卡视频在线观看欧美| 欧美日韩在线观看h| 夫妻午夜视频| 非洲黑人性xxxx精品又粗又长| 午夜精品一区二区三区免费看| 1000部很黄的大片| 热99在线观看视频| 69人妻影院| 成人鲁丝片一二三区免费| 欧美 日韩 精品 国产| 男女那种视频在线观看| 九九在线视频观看精品| 26uuu在线亚洲综合色| 国精品久久久久久国模美| 啦啦啦啦在线视频资源| 一个人看视频在线观看www免费| 男女啪啪激烈高潮av片| 97在线视频观看| 午夜福利视频1000在线观看| 国产日韩欧美在线精品| 免费大片黄手机在线观看| 天天躁夜夜躁狠狠久久av| 黄片无遮挡物在线观看| 国产精品女同一区二区软件| 狂野欧美白嫩少妇大欣赏| 久久精品熟女亚洲av麻豆精品 | 久久久欧美国产精品| 3wmmmm亚洲av在线观看| 亚洲国产色片| 黄色日韩在线| av在线天堂中文字幕| 丝袜喷水一区| 精品久久久久久久久亚洲| 日本欧美国产在线视频| 国产精品一区二区三区四区免费观看| 男女边吃奶边做爰视频| 又爽又黄a免费视频| 亚洲乱码一区二区免费版| 国产精品熟女久久久久浪| 色综合色国产| 亚洲aⅴ乱码一区二区在线播放| 啦啦啦韩国在线观看视频| 日韩不卡一区二区三区视频在线| 精品一区二区三区视频在线| h日本视频在线播放| 91久久精品国产一区二区三区| 尤物成人国产欧美一区二区三区| 亚洲精品自拍成人| 亚洲欧美中文字幕日韩二区| 大陆偷拍与自拍| 91精品伊人久久大香线蕉| 国产伦精品一区二区三区四那| 久久99热这里只有精品18| 久热久热在线精品观看| 只有这里有精品99| 国产麻豆成人av免费视频| 寂寞人妻少妇视频99o| 日本一二三区视频观看| 国产白丝娇喘喷水9色精品| 一夜夜www| 麻豆成人av视频| 一级二级三级毛片免费看| 成人午夜精彩视频在线观看| 亚洲国产欧美人成| 久久精品熟女亚洲av麻豆精品 | 啦啦啦中文免费视频观看日本| 亚洲欧美成人精品一区二区| 日韩av免费高清视频| 尤物成人国产欧美一区二区三区| 亚洲精品日韩av片在线观看| 乱系列少妇在线播放| 亚洲精品aⅴ在线观看| av女优亚洲男人天堂| 国产午夜福利久久久久久| 午夜免费激情av| 最近中文字幕2019免费版| 一区二区三区免费毛片| 亚洲最大成人中文| 99九九线精品视频在线观看视频| 在线播放无遮挡| 国产精品日韩av在线免费观看| 国产淫语在线视频| 久久久国产一区二区| 一级毛片aaaaaa免费看小| 国产av不卡久久| 免费观看精品视频网站| 国产精品三级大全| 美女被艹到高潮喷水动态| 精品一区在线观看国产| 久久久精品免费免费高清| 亚洲欧洲日产国产| 国产伦理片在线播放av一区| 久久久精品欧美日韩精品| kizo精华| 成人午夜高清在线视频| 网址你懂的国产日韩在线| 一区二区三区乱码不卡18| 免费观看av网站的网址| 久久久久网色| 在线观看av片永久免费下载| 三级毛片av免费| 免费看a级黄色片| 成年免费大片在线观看| 久久国产乱子免费精品| 高清毛片免费看| 黑人高潮一二区| 久久精品熟女亚洲av麻豆精品 | 亚洲av一区综合| 少妇高潮的动态图| 欧美激情在线99| 国产精品99久久久久久久久| 综合色av麻豆| 一个人看视频在线观看www免费| 国产精品爽爽va在线观看网站| 天堂av国产一区二区熟女人妻| 97在线视频观看| 成人午夜精彩视频在线观看| 麻豆国产97在线/欧美| 日韩在线高清观看一区二区三区|