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

    HS算法在生物生態(tài)組合污水處理中的應用

    2015-04-14 12:28:06娥,王
    計算機工程與應用 2015年1期
    關鍵詞:搜索算法污水處理向量

    李 娥,王 波

    重慶大學 計算機學院,重慶 400044

    1 引言

    生物生態(tài)組合污水處理[1]是將傳統(tǒng)生物處理系統(tǒng)和生態(tài)處理系統(tǒng)聯(lián)合起來對污水進行處理,因地制宜,充分發(fā)揮各自的優(yōu)勢,形成投資省,能耗低,運行穩(wěn)定,管理方便的污水處理系統(tǒng)。吳敏、陳秀榮等[2-3]提出的生物生態(tài)組合工藝優(yōu)于常規(guī)二級生物處理效果,系統(tǒng)工程造價和能耗顯著降低,達到了生物段靈活、高效和生態(tài)段低耗、穩(wěn)定的最佳組合,然而均未給出合理分配生物段生態(tài)段污水處理負荷的解決方法。眾所周知,氧氣是微生物的主要動力能源,污水處理過程中曝氣量的多少直接影響著生物段中污染物的降解量和污水處理質(zhì)量。曝氣時間過長會導致污泥老化,曝氣時間過短會使生物段生化反應不充分,最終超過生態(tài)段污水處理能力導致出水超標。因此掌握合理的曝氣時間,不僅是污水處理質(zhì)量的保證,也成為污水處理負荷分配的關鍵。

    污水處理是一個高度并行非線性的動態(tài)處理系統(tǒng),國際水質(zhì)協(xié)會提出的Activated Sludge Models(ASM)[4]系列活性污泥模型受到大力推廣,但該系列模型結(jié)構過于復雜,待定參數(shù)過多,因此國內(nèi)外眾多學者將神經(jīng)元網(wǎng)絡建模應用到污水處理過程中[5-7]。徑向基(Radial Basis Function,RBF)神經(jīng)網(wǎng)絡以其訓練速度快,能以任意精度逼近任意連續(xù)函數(shù)等優(yōu)點得到廣泛的應用。近幾年,采用遺傳算法、退火算法和粒子群算法等智能算法優(yōu)化神經(jīng)網(wǎng)絡參數(shù)的方法被提出,然而遺傳算法和退火算法計算復雜度較高,且在搜索最優(yōu)參數(shù)時較為耗時,粒子群雖然具有較快的收斂速度,但也易陷入局部極小值。和聲搜索(Harmony search,HS)算法[8-9]是韓國學者Geem等人受音樂家創(chuàng)作優(yōu)美和聲的啟發(fā),提出的一種啟發(fā)式全局搜索算法,具有較好的尋優(yōu)精度和跳出局部最優(yōu)的能力。目前,該方法已在多維多極值函數(shù)優(yōu)化、土坡穩(wěn)定分析等問題中得到了廣泛應用。為了更好地解決優(yōu)化問題,針對基本和聲搜索算法中記憶庫取值概率(Harmony Memory Considering Rate,HMCR)、音調(diào)調(diào)節(jié)概率(Pitch Adjusting Rate,PAR)和音調(diào)微調(diào)帶寬(Bandwidth,BW)都是采用常數(shù)值這一缺點。Mahdavi等[10]提出了基于動態(tài)PAR和BW的改進和聲搜索算法;高立群等[11]將之與粒子群搜索算法相結(jié)合,提出了自適應和聲粒子群搜索算法;陳瑩珍等[12]將混沌機制作用在種群的初始化和算法停滯時的變異中,提出了混沌自適應和聲搜索算法;杜文莉等[13]將和聲搜索算法與混沌的粒子群算法相融合,該算法收斂速度快,穩(wěn)定性高,已經(jīng)成功應用于重油熱解模型的參數(shù)估計。

    本文在前人研究成果的基礎上,將混沌優(yōu)化思想融合到整個和聲搜索算法中,提出了改進混沌和聲搜索算法(Improved Chaotic Harmony Search,ICHS)。利用ICHS算法來優(yōu)化RBF神經(jīng)網(wǎng)絡參數(shù),建立ICHS-RBF智能模型,并將該模型應用到生物生態(tài)組合污水處理過程中。通過該模型實時有效地選擇合適的處理工況,以曝氣時間為指標進行污水處理負荷分配,預測污水出水水質(zhì),并在LabVIEW(Laboratory Virtual Instrument Engineering Workbench)平臺上集成該模型的負荷分配策略和出水水質(zhì)預測功能,對污水處理過程進行動態(tài)模擬。

    2 改進混沌和聲搜索算法

    2.1 基本和聲搜索算法

    和聲搜索算法是最近幾年提出的啟發(fā)式全局搜索算法,其靈感來源于音樂家創(chuàng)作優(yōu)美樂曲的過程。音樂家們通過反復調(diào)整樂隊中各個樂器的音調(diào)以最終獲得由聽眾的感受(美學評價)所決定的最優(yōu)美動聽的樂曲,和聲搜索算法就是尋找對目標函數(shù)適應度最高的全局最優(yōu)解。和聲搜索算法可以概括為三個簡單的步驟:初始化;產(chǎn)生新和聲;更新和聲記憶庫(Harmony Memory,HM)。首先產(chǎn)生HMS(Harmony Memory Size)個初始和聲向量存儲至HM中。然后對和聲向量的每一個分量以概率HMCR在HM中進行搜索,以概率1-HMCR在HM外的解空間中隨機產(chǎn)生。對來自于HM中的分量,以概率PAR決定是否需要對其進行量為BW的變異微調(diào)。最后,將新的和聲向量與HM中的最差解進行比較,若新解優(yōu)于最差解,則替換之,更新和聲記憶庫。如此迭代循環(huán),直到滿足終止條件。

    2.2 混沌優(yōu)化算法

    混沌運動具有很強的隨機性,遍歷性和規(guī)律性,能夠在一個特定區(qū)域非重復地到達所有的狀態(tài),其看似雜亂的變化過程,其實隱藏著精細的內(nèi)在規(guī)律。式(1)是典型的混沌Logistic映射:

    其中μ為混沌吸引子,j為迭代次數(shù)。當μ為4時,系統(tǒng)處于完全混沌狀態(tài)?;诙屋d波的混沌優(yōu)化算法步驟如下:

    步驟1混沌變量的初始值記為x=(x1,x2,…,xm),根據(jù)式(2)將混沌變量映射到對應的Logistic變量空間中,記為y=(y1,y2,…,ym)。

    步驟2采用μ=4時的式(1)對y進行混沌映射,同時計算yj的適應度值fk(yj)。如果fk+1(yj+1)優(yōu)于fk(yj),則放棄yj,否則放棄yj+1,k繼續(xù)迭代。

    步驟3若迭代若干步后,yj始終保持適應度值fk(yj)最優(yōu),則以yj為初始值按式(1)再次進行混沌映射,即第二次載波,產(chǎn)生新的混沌和聲向量序列。

    步驟4如果滿足終止條件則停止搜索,輸出最優(yōu)yj和f(yj)的值。

    2.3 改進混沌和聲搜索算法

    (1)混沌優(yōu)化策略

    混沌算法優(yōu)化和聲搜索算法策略主要表現(xiàn)在三個方面:其一,初始化HM。從混沌和聲向量序列中取HMS個和聲向量初始化HM,而非在解空間中隨機產(chǎn)生。其二,當以1-HMCR的概率在HM外獲取新解時,同樣在當前混沌和聲向量序列中取值而非隨機產(chǎn)生新解。其三,采用基于二次載波的混沌算法優(yōu)化HM中的和聲向量,產(chǎn)生潛在的最優(yōu)候選解。從當前混沌和聲向量序列中取解,并與HM中的和聲向量進行比較。如果該和聲優(yōu)于HM中最差和聲,則代替之;若差于最差和聲向量或與HM中某和聲向量相同,則不作任何操作,繼續(xù)取下一個和聲進行比較,直至某和聲向量在若干次迭代中均保持為最優(yōu)解,則以該和聲向量為初始值,并以量BW進行變異調(diào)節(jié),產(chǎn)生新的混沌和聲向量序列。

    混沌運動可以使搜索遍布整個解空間,保證解的多樣性,以免陷入局部最優(yōu);采用基于二次載波的混沌優(yōu)化算法將潛在的最優(yōu)和聲向量選入HM中,加快了算法的收斂速度。

    (2)記憶庫取值概率HMCR

    HMCR決定了新的和聲向量是繼承HM中的優(yōu)良解信息還是從解空間中隨機產(chǎn)生。在算法運行前期,需要較小的HMCR以便從全局范圍中搜索潛在的最優(yōu)和聲向量;在算法運行中后期,已經(jīng)達到最優(yōu)解鄰域,需要較大的HMCR以便進行局部搜索。因此采用文獻[12]中的方法,確定HMCR的動態(tài)變化如下式所示,其中gn表示當前迭代次數(shù),NI為最大迭代次數(shù)。Hmax和Hmin分別代表HMCR的最大最小值。

    (3)音調(diào)調(diào)節(jié)概率PAR和微調(diào)帶寬BW

    PAR代表了來自HM中的新和聲是否需要進行微調(diào)變異的概率。在算法搜索前期,較小的值有助于保留原始和聲的信息,盡快找到最優(yōu)解;在算法搜索后期,較大的值促使新的和聲進行微調(diào)變異,促進解的多樣性,幫助算法跳出局部最優(yōu)。采用文獻[10]中的方法,PAR的變化如下所示:

    音調(diào)微調(diào)帶寬BW的值隨著迭代次數(shù)以指數(shù)趨勢下降,其變化如式所示:

    2.4 ICHS優(yōu)化RBF神經(jīng)網(wǎng)絡算法流程

    采用ICHS算法優(yōu)化RBF神經(jīng)網(wǎng)絡各項參數(shù)的基本思想是將網(wǎng)絡學習過程看成是在最優(yōu)解空間中搜索最優(yōu)和聲向量的過程,和聲向量是由RBF網(wǎng)絡的基函數(shù)中心C、寬度σ和權值ω組成。衡量和聲向量是否最優(yōu)的標準,稱為適應度函數(shù)。在神經(jīng)網(wǎng)絡優(yōu)化問題中,適應度函數(shù)需考量網(wǎng)絡逼近精度和網(wǎng)絡復雜度兩方面的性能,使兩者的綜合指標達到最小。網(wǎng)絡的逼近精度是指網(wǎng)絡期望輸出和實際輸出差值的平方和,網(wǎng)絡復雜度則主要由隱藏層神經(jīng)元個數(shù)決定。因此確定適應度函數(shù)為:

    其中,P為網(wǎng)絡期望輸出,T為實際輸出,R為訓練樣本的個數(shù),s為隱藏層神經(jīng)元個數(shù)。可見,當標準誤差較小,隱藏層神經(jīng)元個數(shù)較少時,適應度函數(shù)值則越小。

    算法的終止條件有兩種可能:一是算法的迭代次數(shù)大于規(guī)定的最大迭代次數(shù),另外一種則是當適應度值小于給定的ε時,即網(wǎng)絡的逼近精度和網(wǎng)絡復雜度均達到某種較優(yōu)狀態(tài),算法終止。

    ICHS優(yōu)化RBF神經(jīng)網(wǎng)絡的算法描述如下:

    步驟1初始化。設置各類參數(shù)的初始值,算法的最大迭代次數(shù)NI以及初始化和聲記憶庫,確定網(wǎng)絡的拓撲結(jié)構。

    步驟2將和聲向量轉(zhuǎn)化為神經(jīng)網(wǎng)絡各參數(shù),計算適應度值。

    步驟3利用混沌優(yōu)化算法選擇適應度值最優(yōu)的和聲向量,產(chǎn)生混沌和聲向量序列。

    步驟4通過動態(tài)HMCR、PAR和BW產(chǎn)生新的和聲。

    步驟5將新的和聲向量與HM中的和聲向量作比較,更新和聲記憶庫。

    步驟6重復步驟2~5,直到達到終止條件。

    ICHS優(yōu)化RBF神經(jīng)網(wǎng)絡的流程如圖1所示。

    圖1 ICHS優(yōu)化RBF神經(jīng)網(wǎng)絡的算法流程圖

    3 改進和聲搜索算法在污水處理中的應用

    以重慶某采用序批式生物膜法(SBBR)[15]和人工濕地相結(jié)合的小城鎮(zhèn)污水處理廠為實驗對象,以其歷史數(shù)據(jù)實時建立基于ICHS-RBF的工況選擇智能模型和出水水質(zhì)預測模型,并對該模型進行測試,同時動態(tài)模擬該污水廠的污水處理過程。

    3.1 ICHS-RBF污水處理智能模型

    本文建立兩種不同功能的模型:工況選擇模型和出水預測模型。將生物生態(tài)組合污水處理工藝分為兩種工況:生態(tài)污水處理工況和生物生態(tài)協(xié)同污水處理工況。先利用工況選擇模型對污水進行分類,若僅采用生態(tài)處理方法即可達到排放標準,則選擇生態(tài)污水處理工況,無需再進行污水處理負荷的分配;若污水的處理負荷較大,僅靠生態(tài)處理方法已無法完成對污水的凈化,則選擇用生物生態(tài)協(xié)同污水處理工況,具體流程如圖2所示。

    圖2 生物生態(tài)協(xié)同污水處理流程圖

    工況選擇模型根據(jù)進水水質(zhì)參數(shù)確定合理的處理工況,分別對生態(tài)污水處理工況和生物生態(tài)協(xié)同污水處理工況編號為0和1,以污水流速、污水溫度、COD濃度、TN濃度和TP濃度作為輸入節(jié)點,工況編號作為輸出節(jié)點。

    出水預測模型是利用污水處理過程中進出水水質(zhì)參數(shù)之間的非線性映射關系,預測出水水質(zhì),同時通過確定生物生態(tài)協(xié)同工況中生物段的曝氣時間對污水處理負荷進行合理分配,盡可能節(jié)省能源。為了提高出水預測精度,采用分步預測的方法,即在生物處理段和生態(tài)處理段各自建立網(wǎng)絡模型。當選擇工況0時,曝氣時間為0,污水進水參數(shù)直接作為生態(tài)段的輸入?yún)?shù)值;當選擇工況1時,生物段的出水參數(shù)除曝氣時間外,其余作為生態(tài)段的輸入?yún)?shù)值,如圖3所示。

    3.2 數(shù)據(jù)預處理

    實驗數(shù)據(jù)分別從不同季節(jié)中的平時水質(zhì)數(shù)據(jù)和節(jié)假日期間水質(zhì)數(shù)據(jù)中隨機抽取而來,共計96組,隨機取出6組作為測試數(shù)據(jù),其余作為訓練數(shù)據(jù)。訓練網(wǎng)絡時,首先應對樣本數(shù)據(jù)進行預處理,剔除異常數(shù)據(jù),并按公式(7)進行歸一化處理。

    圖3 出水預測網(wǎng)絡拓撲圖

    X、Y分別是處理前后的數(shù)據(jù)值;Xmax、Xmin,Ymax、Ymin分別是處理前后樣本的最大值和最小值。預測數(shù)據(jù)同樣采用樣本數(shù)據(jù)的映射結(jié)構進行歸一化處理,預測結(jié)果進行反歸一化操作。

    3.3 結(jié)果分析

    完成ICHS算法對RBF神經(jīng)網(wǎng)絡的參數(shù)尋優(yōu),最終工況選擇模型的工況序號結(jié)果是1、1、1、0、1、0,生物段和生態(tài)段的出水預測結(jié)果如表1所示。第4、6組數(shù)據(jù)采用工況0便可達到污水排放標準;其他均采用生物生態(tài)協(xié)同污水處理工況,曝氣時間預測值與實測值很接近,分別曝氣1.5 h、3.2 h、1.9 h和2.4 h后,排水到生態(tài)段作進一步的處理,出水水質(zhì)滿足一級排放B標準,實現(xiàn)了利用曝氣時間來進行污水處理負荷的分配,同時出水水質(zhì)參數(shù)的預測值與實測值也相差不大。其中,TN濃度的預測效果最佳,效果比較穩(wěn)定,相對誤差維持在1%~4%之間;TP濃度的預測效果最差,抖動較大,最大相對誤差達到6.3%;COD和NH4+-N預測效果良好。

    3.4 基于LabVIEW的污水處理過程動態(tài)模擬

    近年來,基于LabVIEW的模擬運行平臺[16-17]在污水處理領域得到了廣泛的應用,它不僅方便對污水廠進行整體設計,形象地演示污水處理過程,進行大量數(shù)值計算,而且還可實現(xiàn)污水處理過程的實時控制。本系統(tǒng)利用MATLAB Script節(jié)點調(diào)用ICHS-RBF智能模型MATLAB程序,待模型訓練成功后,輸入測試數(shù)據(jù)即可顯示出水水質(zhì)的預測結(jié)果,同時根據(jù)工況選擇模型給出的工況編號和對應的曝氣時間,將自動調(diào)用相應的工況子程序來動態(tài)演示污水處理的具體流程。圖4為工況1時該廠的污水處理流程的模擬演示。

    表1 出水水質(zhì)預測值與實測值比較

    圖4 工況1污水處理過程模擬

    4 結(jié)束語

    本文以改進的混沌和聲搜索算法對RBF神經(jīng)網(wǎng)絡進行優(yōu)化并應用于污水處理領域,建立ICHS-RBF智能模型來確定合理的污水處理工況和曝氣時間,達到合理分配生物段和生態(tài)段污水處理負荷的目的;出水水質(zhì)預測的最高相對誤差為6.3%,最低相對誤差為0.4%,預測效果較好;基于該模型的生物生態(tài)組合污水處理過程模擬系統(tǒng)具有直觀性好,人機交互方便等特點,具有一定的實用價值。

    [1]苑宏英,周琪.生物-生態(tài)聯(lián)合污水處理系統(tǒng)[J].環(huán)境污染治理技術與設備,2005,6(10):19-23.

    [2]吳敏,楊健,馬運才,等.生物-生態(tài)過濾工藝處理城鎮(zhèn)污水[J].中國給水排水,2002,18(4):37-39.

    [3]陳秀榮,周琪,張杰.生物-生態(tài)協(xié)同工藝深度處理城市污水[J].哈爾濱工業(yè)大學學報,2008,40(4):650-655.

    [4]Henze M,Gujer W,Mino T,et al.Activated Sludge Models ASM1,ASM2,ASM2d and ASM3.IWA task group on mathematical modeling for design and operation of biological wastewater treatment,2000.

    [5]Nasr M S,Moustafa M A E,Seif H A E,et al.Application of artification neuralnetwork for the prediction of EL-AGAMY wastewater treatment plant performance-EGYPT[J].Alexandria Engineering Journal,2012,51:37-43.

    [6]韓建秋,喬俊飛.改進Elman網(wǎng)絡的污水處理過程建模[J].計算機與應用化學,2011,28(7):847-850.

    [7]吳凡松,彭永臻,王維斌,等.RBF網(wǎng)絡用于活性污泥系統(tǒng)仿真研究[J].中國積水排水,2003,19(9):48-49.

    [8]Geem Z W,Kim J H,Loganathan G V.A new heuristic optimization algorithm:harmony search[J].Simulation,2001,76(2):60-68.

    [9]雍龍泉.和聲搜索算法研究進展[J].計算機系統(tǒng)應用,2011,20(7):244-248.

    [10]Mahdavi M,F(xiàn)esanghary M,Damangir E.An improved harmony search algorithm for solving optimization problem[J].Applied Mathematics and Computation,2007,188(2):1567-1597.

    [11]高立群,葛延峰,孔芝,等.自適應和聲粒子群搜索算法[J].控制與決策,2010,25(7):1101-1104.

    [12]陳瑩珍,高岳林.混沌的自適應和聲搜索算法[J].太原理工大學學報,2011,42(2):141-144.

    [13]杜文莉,張海龍,錢鋒.融合和聲搜索的混沌粒子群優(yōu)化算法及工業(yè)應用[J].清華大學學報:自然科學版,2012,52(3):325-330.

    [14]劉俊商,秀芹,盧建剛,等.采用中心聚類與PSO的RBF網(wǎng)絡設計方法[J].計算機工程與應用,2009,45(36):212-215.

    [15]姚敬博,李亞峰,王勛,等.SBBR工藝脫氮除磷理論及影響因素[J].工業(yè)安全與環(huán)保,2009,35(8):18-19.

    [16]喬俊飛,韓紅桂,張穎.基于LabVIEW的污水處理過程動態(tài)模擬研究[J].儀器儀表學報,2008,29(4):879-882.

    [17]Huaigang Z,Simeng F.Application of virtual instrument in sewage multi-parameter online monitoring system[C]//Proceedings of IEEE International Workshop on Intelligent Systems and Applications,2009:1788-1791.

    猜你喜歡
    搜索算法污水處理向量
    我國鄉(xiāng)鎮(zhèn)污水處理模式的探討
    向量的分解
    太原市61個村要建污水處理設施嚴禁直排入河
    改進的和聲搜索算法求解凸二次規(guī)劃及線性規(guī)劃
    聚焦“向量與三角”創(chuàng)新題
    涌浪機在污水處理中的應用
    向量垂直在解析幾何中的應用
    向量五種“變身” 玩轉(zhuǎn)圓錐曲線
    基于汽車接力的潮流轉(zhuǎn)移快速搜索算法
    基于逐維改進的自適應步長布谷鳥搜索算法
    久99久视频精品免费| 自拍偷自拍亚洲精品老妇| 中文资源天堂在线| 色av中文字幕| 久久久久国内视频| 亚洲性夜色夜夜综合| 热99在线观看视频| 一进一出抽搐gif免费好疼| 一本综合久久免费| 他把我摸到了高潮在线观看| 成人特级黄色片久久久久久久| 伦理电影大哥的女人| 最新中文字幕久久久久| 两性午夜刺激爽爽歪歪视频在线观看| 欧美国产日韩亚洲一区| 欧美成狂野欧美在线观看| 色哟哟·www| 一级a爱片免费观看的视频| 欧美区成人在线视频| 日日摸夜夜添夜夜添小说| 夜夜夜夜夜久久久久| 91久久精品国产一区二区成人| 国产精品免费一区二区三区在线| 成人毛片a级毛片在线播放| 一本综合久久免费| 别揉我奶头 嗯啊视频| 亚洲最大成人手机在线| 天堂√8在线中文| 美女大奶头视频| 麻豆成人av在线观看| 高潮久久久久久久久久久不卡| 三级毛片av免费| 国产精品女同一区二区软件 | 99热精品在线国产| 成人性生交大片免费视频hd| 欧美一区二区亚洲| 亚洲片人在线观看| 国产免费一级a男人的天堂| 国产亚洲精品av在线| 久久精品国产亚洲av天美| 精品熟女少妇八av免费久了| 久久久久国内视频| 亚洲aⅴ乱码一区二区在线播放| 九九热线精品视视频播放| 国产成人福利小说| 久久精品影院6| 热99在线观看视频| 午夜福利在线在线| 精品一区二区三区人妻视频| 久久国产精品影院| 精品福利观看| 免费看光身美女| 欧美绝顶高潮抽搐喷水| 国产探花在线观看一区二区| 亚洲第一电影网av| 老司机福利观看| 免费在线观看亚洲国产| av女优亚洲男人天堂| 日韩欧美精品免费久久 | 真人做人爱边吃奶动态| 高潮久久久久久久久久久不卡| 国产av一区在线观看免费| 国产成人啪精品午夜网站| 在线看三级毛片| 一本久久中文字幕| 亚洲精品影视一区二区三区av| 日本a在线网址| 一级毛片久久久久久久久女| 久久久久国产精品人妻aⅴ院| 成人无遮挡网站| 久久久久久大精品| 亚洲欧美日韩无卡精品| 国产三级中文精品| 久久精品91蜜桃| 黄色一级大片看看| 亚洲成人久久爱视频| 男女之事视频高清在线观看| 欧美一区二区国产精品久久精品| 黄色女人牲交| 美女免费视频网站| 欧美日韩黄片免| 我的女老师完整版在线观看| 亚洲精品在线观看二区| 两个人的视频大全免费| 亚洲中文字幕日韩| 国产精品国产高清国产av| 韩国av一区二区三区四区| 免费在线观看日本一区| 欧美国产日韩亚洲一区| 老司机午夜福利在线观看视频| 麻豆成人午夜福利视频| 俺也久久电影网| 欧美zozozo另类| 免费看美女性在线毛片视频| 美女高潮喷水抽搐中文字幕| 久久人人精品亚洲av| 国产乱人视频| 身体一侧抽搐| 亚洲中文字幕一区二区三区有码在线看| 欧美日韩瑟瑟在线播放| 亚洲精品在线美女| 在线天堂最新版资源| 亚洲专区国产一区二区| 成人国产综合亚洲| 国产亚洲精品综合一区在线观看| 亚洲天堂国产精品一区在线| 高清在线国产一区| 亚洲第一电影网av| 欧美国产日韩亚洲一区| 亚洲天堂国产精品一区在线| 色综合欧美亚洲国产小说| 老司机午夜福利在线观看视频| 国产成人欧美在线观看| 国产精品久久久久久久久免 | 亚洲av免费高清在线观看| 日韩欧美精品免费久久 | 91九色精品人成在线观看| 亚洲在线观看片| 中文字幕人妻熟人妻熟丝袜美| 精品久久久久久成人av| 中文字幕免费在线视频6| 一级av片app| 久久草成人影院| 国产69精品久久久久777片| 午夜激情欧美在线| 亚洲人与动物交配视频| 久久久国产成人免费| 午夜精品久久久久久毛片777| 最后的刺客免费高清国语| 中国美女看黄片| 色播亚洲综合网| 成人午夜高清在线视频| 久久6这里有精品| 午夜激情福利司机影院| 日本一本二区三区精品| 亚洲精品一卡2卡三卡4卡5卡| 人人妻,人人澡人人爽秒播| 好男人电影高清在线观看| 日本黄色视频三级网站网址| 老熟妇乱子伦视频在线观看| 成人高潮视频无遮挡免费网站| 欧美激情国产日韩精品一区| 免费人成在线观看视频色| 久久国产精品人妻蜜桃| 高潮久久久久久久久久久不卡| 国产精品野战在线观看| 欧美潮喷喷水| 亚洲中文日韩欧美视频| 精品久久久久久久人妻蜜臀av| 国产欧美日韩精品一区二区| 中文字幕久久专区| 脱女人内裤的视频| 免费在线观看日本一区| 波多野结衣巨乳人妻| 国产精品不卡视频一区二区 | www日本黄色视频网| 3wmmmm亚洲av在线观看| 中文资源天堂在线| 亚洲午夜理论影院| 欧美3d第一页| 欧洲精品卡2卡3卡4卡5卡区| 欧美国产日韩亚洲一区| 亚洲人成网站高清观看| 亚洲av中文字字幕乱码综合| 色综合婷婷激情| 久久精品国产亚洲av涩爱 | 国产精品野战在线观看| 少妇的逼水好多| netflix在线观看网站| 特级一级黄色大片| 伦理电影大哥的女人| 久久精品91蜜桃| 91午夜精品亚洲一区二区三区 | 日本与韩国留学比较| 深夜a级毛片| 婷婷六月久久综合丁香| 国产高清视频在线观看网站| 综合色av麻豆| 狠狠狠狠99中文字幕| 日韩欧美精品免费久久 | 亚洲精品在线美女| 日本五十路高清| 在线观看一区二区三区| 色吧在线观看| 午夜福利免费观看在线| 狠狠狠狠99中文字幕| 桃红色精品国产亚洲av| 亚洲人成伊人成综合网2020| 如何舔出高潮| 午夜久久久久精精品| 看片在线看免费视频| 国产白丝娇喘喷水9色精品| 身体一侧抽搐| 美女免费视频网站| 国产精品免费一区二区三区在线| 久久精品综合一区二区三区| 日本免费a在线| 毛片一级片免费看久久久久 | 久久婷婷人人爽人人干人人爱| 男女床上黄色一级片免费看| 嫩草影院精品99| 久久久久久久精品吃奶| 国产人妻一区二区三区在| 久久草成人影院| 日本一本二区三区精品| 国产熟女xx| 变态另类成人亚洲欧美熟女| 欧美又色又爽又黄视频| 久久亚洲精品不卡| 色哟哟·www| 国产成人影院久久av| 国产亚洲欧美98| 内地一区二区视频在线| 少妇的逼水好多| 波野结衣二区三区在线| 又粗又爽又猛毛片免费看| 男女做爰动态图高潮gif福利片| 国产精品三级大全| 国产精品久久久久久人妻精品电影| 麻豆成人午夜福利视频| 亚洲不卡免费看| 欧美高清成人免费视频www| 国产黄片美女视频| 日韩免费av在线播放| 最新中文字幕久久久久| 天堂动漫精品| av视频在线观看入口| 精品一区二区三区视频在线| 欧美日韩亚洲国产一区二区在线观看| 看片在线看免费视频| 国产老妇女一区| 国产三级在线视频| 少妇的逼好多水| 国模一区二区三区四区视频| 两人在一起打扑克的视频| 成人美女网站在线观看视频| 精品不卡国产一区二区三区| 中文字幕人成人乱码亚洲影| 国产国拍精品亚洲av在线观看| 51午夜福利影视在线观看| 日本一二三区视频观看| 成人亚洲精品av一区二区| 在线国产一区二区在线| 欧美区成人在线视频| av国产免费在线观看| 变态另类丝袜制服| 色精品久久人妻99蜜桃| 美女高潮的动态| 成人欧美大片| 在线观看66精品国产| 亚洲真实伦在线观看| 宅男免费午夜| 毛片女人毛片| 内地一区二区视频在线| 国产亚洲精品综合一区在线观看| 国产老妇女一区| 一级a爱片免费观看的视频| 九色成人免费人妻av| 啪啪无遮挡十八禁网站| 在线免费观看不下载黄p国产 | 中文字幕精品亚洲无线码一区| 亚洲av成人精品一区久久| 人妻制服诱惑在线中文字幕| 男人舔奶头视频| 中文字幕精品亚洲无线码一区| 久久久久精品国产欧美久久久| 久久久久免费精品人妻一区二区| 动漫黄色视频在线观看| 午夜福利在线在线| 亚洲av.av天堂| 91在线精品国自产拍蜜月| 18禁黄网站禁片午夜丰满| 直男gayav资源| 国产成+人综合+亚洲专区| 欧美+亚洲+日韩+国产| 草草在线视频免费看| 最近视频中文字幕2019在线8| 日韩亚洲欧美综合| 别揉我奶头 嗯啊视频| 国产伦在线观看视频一区| 国内精品美女久久久久久| 最新中文字幕久久久久| 国内毛片毛片毛片毛片毛片| 美女 人体艺术 gogo| 又黄又爽又刺激的免费视频.| 日韩精品中文字幕看吧| 国产三级中文精品| 天堂√8在线中文| 日韩高清综合在线| 精品99又大又爽又粗少妇毛片 | 看片在线看免费视频| 国产视频一区二区在线看| www日本黄色视频网| 深爱激情五月婷婷| 欧美色视频一区免费| 99久国产av精品| 级片在线观看| 日韩av在线大香蕉| 国产精品av视频在线免费观看| 中文字幕高清在线视频| 色噜噜av男人的天堂激情| av在线观看视频网站免费| 精品一区二区免费观看| 嫩草影院新地址| 中出人妻视频一区二区| 成人永久免费在线观看视频| 久久精品人妻少妇| 有码 亚洲区| 一本一本综合久久| 最近最新中文字幕大全电影3| 亚洲色图av天堂| 少妇被粗大猛烈的视频| 性色av乱码一区二区三区2| 美女 人体艺术 gogo| 国产高清三级在线| 网址你懂的国产日韩在线| 搡老岳熟女国产| 亚洲av电影在线进入| 欧美xxxx性猛交bbbb| 成人欧美大片| 免费在线观看影片大全网站| 欧美黄色片欧美黄色片| 久久久久国内视频| 校园春色视频在线观看| 日日摸夜夜添夜夜添av毛片 | 国产亚洲精品久久久久久毛片| 婷婷六月久久综合丁香| 欧美成人性av电影在线观看| 色视频www国产| 小蜜桃在线观看免费完整版高清| 91久久精品电影网| 色吧在线观看| 午夜福利视频1000在线观看| 精品久久久久久久久久免费视频| 精品国产亚洲在线| 久久中文看片网| 久久性视频一级片| www.熟女人妻精品国产| 婷婷精品国产亚洲av| 九九久久精品国产亚洲av麻豆| 色精品久久人妻99蜜桃| 亚洲久久久久久中文字幕| 亚洲欧美日韩东京热| 欧美bdsm另类| 女人十人毛片免费观看3o分钟| 中文字幕免费在线视频6| 99久久九九国产精品国产免费| 热99在线观看视频| 成年女人永久免费观看视频| 简卡轻食公司| 99久久精品一区二区三区| 日韩欧美 国产精品| 午夜亚洲福利在线播放| 日本一本二区三区精品| 亚洲专区国产一区二区| 草草在线视频免费看| 国产精品久久视频播放| 内地一区二区视频在线| 国产乱人伦免费视频| 久久这里只有精品中国| 日韩av在线大香蕉| 精品99又大又爽又粗少妇毛片 | 国产亚洲精品av在线| 免费看a级黄色片| 露出奶头的视频| 精品无人区乱码1区二区| 午夜视频国产福利| 国产精品久久久久久亚洲av鲁大| 午夜福利免费观看在线| 美女高潮的动态| 国产伦一二天堂av在线观看| 欧美又色又爽又黄视频| 亚洲狠狠婷婷综合久久图片| 三级毛片av免费| 国产在线男女| 国产精品亚洲一级av第二区| 精品乱码久久久久久99久播| 久久久久国内视频| 高潮久久久久久久久久久不卡| 亚洲 国产 在线| 欧美一区二区亚洲| 男女那种视频在线观看| 在线观看av片永久免费下载| 久久久久久久午夜电影| 欧美午夜高清在线| 一区二区三区四区激情视频 | 久久精品人妻少妇| 色综合站精品国产| 99热这里只有是精品在线观看 | 99在线人妻在线中文字幕| 欧美xxxx性猛交bbbb| 欧美日韩国产亚洲二区| 精华霜和精华液先用哪个| 亚洲七黄色美女视频| 国产成人啪精品午夜网站| 国产精品一区二区性色av| 欧美日韩中文字幕国产精品一区二区三区| 免费在线观看亚洲国产| 黄色配什么色好看| 丰满的人妻完整版| 成人性生交大片免费视频hd| netflix在线观看网站| 免费看a级黄色片| 国产精品久久久久久人妻精品电影| 99久久成人亚洲精品观看| 小说图片视频综合网站| 草草在线视频免费看| 给我免费播放毛片高清在线观看| 夜夜躁狠狠躁天天躁| 成人毛片a级毛片在线播放| 一级毛片久久久久久久久女| 极品教师在线视频| 精品国产亚洲在线| 极品教师在线免费播放| 99久久成人亚洲精品观看| 我的老师免费观看完整版| 俄罗斯特黄特色一大片| 中文字幕免费在线视频6| 国产午夜福利久久久久久| 国产欧美日韩精品亚洲av| 国内精品久久久久久久电影| 黄色一级大片看看| 国产精品,欧美在线| 国产精华一区二区三区| 国产激情偷乱视频一区二区| 亚洲精品成人久久久久久| 男人的好看免费观看在线视频| 亚洲av熟女| 又黄又爽又刺激的免费视频.| 国产精品影院久久| 男女视频在线观看网站免费| 人妻丰满熟妇av一区二区三区| 成人av一区二区三区在线看| 欧美日韩乱码在线| 老女人水多毛片| av在线观看视频网站免费| 一边摸一边抽搐一进一小说| 免费一级毛片在线播放高清视频| 日韩欧美国产一区二区入口| 99国产精品一区二区蜜桃av| 成人美女网站在线观看视频| 午夜福利18| 美女cb高潮喷水在线观看| 一本综合久久免费| 99热只有精品国产| 欧美丝袜亚洲另类 | 欧美精品啪啪一区二区三区| 国产欧美日韩一区二区三| av女优亚洲男人天堂| 亚洲无线观看免费| 国产主播在线观看一区二区| 麻豆国产av国片精品| 久久久久精品国产欧美久久久| 变态另类成人亚洲欧美熟女| a级毛片免费高清观看在线播放| 欧美极品一区二区三区四区| 九九在线视频观看精品| 欧美日韩中文字幕国产精品一区二区三区| 如何舔出高潮| 日韩中文字幕欧美一区二区| 亚洲国产欧美人成| 成人鲁丝片一二三区免费| 久久精品91蜜桃| 久久久精品欧美日韩精品| 黄色丝袜av网址大全| 波多野结衣高清无吗| 欧美xxxx黑人xx丫x性爽| 看黄色毛片网站| av在线蜜桃| 一区二区三区四区激情视频 | 国产精品一区二区免费欧美| 亚洲自拍偷在线| 一本久久中文字幕| 两性午夜刺激爽爽歪歪视频在线观看| 婷婷六月久久综合丁香| 狂野欧美白嫩少妇大欣赏| 免费av不卡在线播放| 啪啪无遮挡十八禁网站| 国内精品久久久久精免费| av中文乱码字幕在线| 久久久久九九精品影院| 内射极品少妇av片p| 国产91精品成人一区二区三区| 非洲黑人性xxxx精品又粗又长| 麻豆成人av在线观看| 免费观看人在逋| 国产色婷婷99| 国产高清激情床上av| xxxwww97欧美| 日韩中文字幕欧美一区二区| 国产伦人伦偷精品视频| 久久精品国产亚洲av天美| av中文乱码字幕在线| 男人和女人高潮做爰伦理| 少妇高潮的动态图| 久久久精品欧美日韩精品| 日韩欧美在线二视频| 国产三级在线视频| 能在线免费观看的黄片| 村上凉子中文字幕在线| 精品无人区乱码1区二区| 国内少妇人妻偷人精品xxx网站| 久久久国产成人免费| 亚洲一区二区三区色噜噜| 69av精品久久久久久| 欧美日韩国产亚洲二区| 有码 亚洲区| av黄色大香蕉| www日本黄色视频网| 无人区码免费观看不卡| 九九在线视频观看精品| 岛国在线免费视频观看| 亚洲午夜理论影院| 十八禁网站免费在线| 日本三级黄在线观看| 亚洲精品一区av在线观看| 亚洲最大成人手机在线| 成人高潮视频无遮挡免费网站| 色综合站精品国产| 高潮久久久久久久久久久不卡| 久久婷婷人人爽人人干人人爱| 久久久国产成人精品二区| 大型黄色视频在线免费观看| 性色avwww在线观看| 国产精品乱码一区二三区的特点| 丰满乱子伦码专区| 精品久久久久久久久久免费视频| 午夜福利欧美成人| 久久99热6这里只有精品| 一进一出抽搐gif免费好疼| 高清毛片免费观看视频网站| 国内精品一区二区在线观看| 90打野战视频偷拍视频| 国产色婷婷99| 国产91精品成人一区二区三区| 午夜福利成人在线免费观看| 日本五十路高清| 精品人妻熟女av久视频| 一区二区三区激情视频| 亚洲五月天丁香| 国产成+人综合+亚洲专区| 国产亚洲欧美在线一区二区| 国产aⅴ精品一区二区三区波| 久久精品国产自在天天线| 男人狂女人下面高潮的视频| 亚洲国产精品sss在线观看| 国产美女午夜福利| 一区二区三区高清视频在线| 毛片女人毛片| 亚洲中文字幕日韩| 俺也久久电影网| 亚洲精品一卡2卡三卡4卡5卡| 精品免费久久久久久久清纯| netflix在线观看网站| 露出奶头的视频| 色综合婷婷激情| 亚洲av日韩精品久久久久久密| 亚洲在线自拍视频| 亚洲自偷自拍三级| 乱码一卡2卡4卡精品| 亚洲综合色惰| 直男gayav资源| 久久精品国产亚洲av香蕉五月| 精品日产1卡2卡| 91av网一区二区| 性欧美人与动物交配| 自拍偷自拍亚洲精品老妇| 亚洲欧美激情综合另类| 婷婷精品国产亚洲av在线| 国产成人av教育| 亚洲 国产 在线| 可以在线观看的亚洲视频| 久久精品91蜜桃| 一本久久中文字幕| 亚洲成人免费电影在线观看| 俺也久久电影网| 国产精品久久久久久久电影| 18美女黄网站色大片免费观看| 免费av不卡在线播放| 色精品久久人妻99蜜桃| 成人国产综合亚洲| 亚洲国产精品成人综合色| 少妇高潮的动态图| 欧美成人a在线观看| 三级男女做爰猛烈吃奶摸视频| 亚洲精品影视一区二区三区av| 欧美不卡视频在线免费观看| 欧美成人性av电影在线观看| 日韩国内少妇激情av| 日本a在线网址| 亚洲av五月六月丁香网| 精品久久久久久久末码| 国产精品亚洲一级av第二区| 日韩欧美精品免费久久 | 日本免费a在线| 日日摸夜夜添夜夜添av毛片 | 国产av不卡久久| 亚洲av中文字字幕乱码综合| 国产一区二区三区在线臀色熟女| 国产精品不卡视频一区二区 | 午夜日韩欧美国产| 99国产综合亚洲精品| 国产精品乱码一区二三区的特点| 亚洲av第一区精品v没综合| 看十八女毛片水多多多| 成人无遮挡网站| 变态另类成人亚洲欧美熟女| 午夜免费激情av| 午夜两性在线视频| 国产亚洲精品久久久com| 每晚都被弄得嗷嗷叫到高潮| 国产美女午夜福利| 精品人妻视频免费看| 性欧美人与动物交配| 午夜福利欧美成人|