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

    混合多目標骨干粒子群優(yōu)化算法在污水處理過程優(yōu)化控制中的應用

    2017-10-13 06:09:58周紅標喬俊飛
    化工學報 2017年9期
    關鍵詞:引導者設定值污水處理

    周紅標,喬俊飛

    ?

    混合多目標骨干粒子群優(yōu)化算法在污水處理過程優(yōu)化控制中的應用

    周紅標1,2,3,喬俊飛1,2

    (1北京工業(yè)大學信息學部,北京 100124;2計算智能和智能系統北京市重點實驗室,北京 100124;3淮陰工學院自動化學院,江蘇淮安 223003)

    通過對污水生化處理過程的分析,選取能耗和罰款最低為優(yōu)化目標,建立污水生化處理過程多目標優(yōu)化控制模型。為了提高Pareto最優(yōu)解集的收斂性和多樣性,提出一種基于Pareto支配和分解的混合多目標骨干粒子群優(yōu)化算法(HBBMOPSO)。該方法采用帶自適應懲罰因子的分解方法選取個體引導者,采用Pareto支配和擁擠距離法維護外部檔案和選取全局引導者。此外,采用精英學習策略增強粒子跳出局部Pareto前沿的能力。最后,將HBBMOPSO與自組織模糊神經網絡預測模型和自組織控制器相結合,實現污水生化處理過程溶解氧和硝態(tài)氮設定值的動態(tài)尋優(yōu)、智能決策和底層跟蹤控制。利用國際基準仿真平臺BSM1進行實驗驗證,結果表明所提HBBMOPSO方法在保證出水水質參數達標的前提下,能夠有效降低污水處理過程的能耗。

    污水;優(yōu)化;過程控制;粒子群;分解

    引 言

    活性污泥法是城市污水處理廠普遍采用的污水生物處理方法[1-2]。該方法通過向廢水中連續(xù)充入空氣,經過一定反應時間后,因好氧微生物繁殖而形成活性污泥。利用活性污泥的生物凝聚、吸附和氧化作用,分解去除污水中的有機污染物[3]。為了滿足出水排放標準,降低罰款(在國際基準模型中,一般利用出水水質來表征罰款),污水處理廠經常處于滿負荷運行,即通過鼓風機和回流泵保持好氧區(qū)溶解氧(dissolved oxygen,O)濃度和厭氧區(qū)硝態(tài)氮(nitrate nitrogen,NO)濃度處于較高水平[4]。但是,鼓風機和回流泵的運轉需要大量的能量供給,這不可避免地增加了運行成本。同時,從生化反應機理來看,只有合適的O和NO才能保證硝化和反硝化反應的順利進行[5]。因此,應根據實際運行情況對O和NO的設定值進行動態(tài)尋優(yōu),構建以降低能耗(energy consumption,EC)和出水水質(effluent quality,EQ)為目標的多目標優(yōu)化控制策略,從而提高污水處理效果和降低運行成本。

    在過去的十幾年里,單目標優(yōu)化控制在污水處理過程中已取得了豐碩成果[6-8]。但是,單目標優(yōu)化控制主要關注能耗指標,容易導致出水水質參數超標,增加運行成本。隨后,有研究者通過權重因子構建損失函數將包含能耗和出水水質的多目標優(yōu)化問題轉換成單目標優(yōu)化問題[9-11]。但是,這種方法不僅權重因子難以確定,而且難以取得能耗和出水水質之間的最佳平衡。多目標優(yōu)化算法能夠克服上述缺點。因此,污水處理過程的多目標優(yōu)化控制受到了廣泛關注[12-18]。文獻[13]利用集成有多目標優(yōu)化算法的交互式軟件IND-NIMBUS構建污水處理過程全流程的運行優(yōu)化模型。文獻[14]采用NSGAII算法建立溫室氣體排放、操作成本和出水污染物濃度的多目標優(yōu)化模型。文獻[15]采用動態(tài)多目標優(yōu)化算法建立開發(fā)成本和污水處理質量的兩目標優(yōu)化模型。文獻[16-18]采用NSGAII算法建立能耗和出水水質的兩目標優(yōu)化模型。通過上述分析發(fā)現,在污水處理過程中采用多目標優(yōu)化控制不僅能夠保證出水水質參數達標,而且能夠有效降低能耗和運行成本。但是上述研究對多目標優(yōu)化算法的討論和改進還存在不足,同時NSGAII算法運行成本較高,解的收斂性和多樣性有待提高,能耗的降低還有一定的空間。

    模擬鳥類覓食行為的粒子群優(yōu)化(particle swarm optimization,PSO)算法[19]具有收斂速度快、概念簡明和易于實現的優(yōu)點,在多目標優(yōu)化領域應用更加普遍[20-22]。但是傳統粒子群優(yōu)化算法在決策空間中探索和開發(fā)能力的平衡嚴重依賴于對飛行參數的調控[22],限制了其在實際工程中的應用。文獻[23]提出的骨干粒子群優(yōu)化(barebones particle swarm optimization,BBPSO)算法取消了速度更新公式,采用關于粒子個體引導者和全局引導者的高斯采樣完成粒子位置更新,無須額外的參數設置,更加適合實際工程應用,因此有研究者將其擴展用于解決多目標優(yōu)化問題[24-25]。

    Pareto解集的收斂性和多樣性是多目標優(yōu)化算法面臨的兩個主要問題。目前,一般通過Pareto支配關系選取個體引導者以促進種群收斂,通過密度估計策略(如自適應網格[20]、擁擠距離[26]等)來保持種群多樣性。文獻[27]提出的MOEA/D方法采用分解策略將多目標優(yōu)化問題分解成一組單目標優(yōu)化問題,具有更低的計算復雜度、更好的收斂性和多樣性。因此,有研究者將分解方法集成到多目標粒子群算法當中[28]。但是,同時考慮綜合利用分解和Pareto支配的多目標粒子群優(yōu)化算法的研究相對較少,尤其是在多目標骨干粒子群優(yōu)化中未見相關報道。因此,本文同時利用Pareto支配和帶自適應懲罰因子的分解方法來構建混合多目標骨干粒子群優(yōu)化(hybrid barebones multi-objective particle swarm optimization,HBBMOPSO)算法,以期提高Pareto解集的收斂性和多樣性,從而提供給決策者更高質量的候選解。

    本文在分析污水處理過程控制變量、操作變量、優(yōu)化目標和出水水質性能指標的基礎上,首先利用前期研究的自組織模糊神經網絡[29](self-organizing fuzzy neural network,SOFNN)建立EC、EQ和出水水質參數的目標函數;然后利用本文所提的HBBMOPSO算法對O和NO設定值進行動態(tài)尋優(yōu),并進一步利用智能決策系統選取偏好解作為當前優(yōu)化周期的優(yōu)化設定值;最后利用前期研究的SOFNN控制器[30]實現對優(yōu)化設定值的底層跟蹤控制。所有實驗均采用活性污泥污水處理國際基準仿真平臺(Benchmark Simulation Model No.1, BSM1)進行驗證。

    1 污水處理過程多目標優(yōu)化問題

    BSM1是由國際水質協會(IWAQ)和歐盟科學技術與合作組織(COST)合作開發(fā),可用于公平評價不同控制策略和優(yōu)化方法。BSM1采用典型的前置反硝化脫氮A/O工藝[31],由生化反應池和二沉池組成,如圖1所示。生化反應池包含5個單元,前2個單元是缺氧區(qū),主要完成反硝化反應,后3個單元是好氧區(qū),主要完成硝化反應。生化反應池第3、4單元的氧氣轉換系數要求保持在恒定值240 d?1。由于入水流量和組分濃度呈現強非線性關系,O和NO控制器的目標就是分別通過調節(jié)第5單元的曝氣量La5和內回流量a來控制第5單元的溶解氧濃度(O,5)和第2單元的硝態(tài)氮濃度(NO,2)[30-32]。多目標優(yōu)化的目的就是分別通過動態(tài)調整O,5和NO,2的設定值[O,sp,NO,sp],實現EC和EQ之間的最佳平衡。

    在污水處理過程中,O,5和NO,2設定值不僅影響EQ,而且與EC有密切關系。因此,為了取得EC和EQ之間的最佳平衡,利用HBBMOPSO算法來處理這對相互沖突的目標函數。

    由于曝氣能耗(AE)和泵送能耗(PE)占總能耗的70%以上[18],因此優(yōu)化問題的EC定義為AE和PE之和,即

    EC=AE+PE (1)

    按照BSM1機理模型的定義,AE和PE為[31]

    (3)

    式中,V和Lai分別為第個單元的體積和氧氣轉換系數;O,sat為溶解氧飽和濃度;為優(yōu)化周期;a、r和w分別為內回流量、污泥回流量和污泥排放量。

    EQ表示向受納水體排放出水污染物所需支付的罰款,按照BSM1定義,EQ的公式為[18,32]

    式中,SS、COD、NO、NKj和BOD5分別為固體懸浮物濃度、化學需氧量、硝態(tài)氮濃度、凱氏氮濃度和5日生化需氧量;e是上清液排出量。

    出水水質參數的約束條件為BSM1基準中給出的達標限定值,具體為[18]

    其中,出水總氮tot為出水硝態(tài)氮NO和出水凱氏氮NKj之和

    綜上所述,污水處理過程約束多目標優(yōu)化問題描述如下

    其中,=[1,2]=[O,sp,NO,sp]為優(yōu)化向量;f()為出水水質參數與優(yōu)化向量之間的關系,=1,2,…,5;和是第個優(yōu)化變量的上下限,=1,2。

    由式(1)~式(4)可以發(fā)現,EC主要與操作變量La5和a有關,EQ主要與5種出水水質參數有關。但是,EC、EQ與O,5、NO,2這兩個優(yōu)化變量之間并沒有明確的數學關系,同時構成約束條件的關鍵出水水質參數也不能在線測量。因此,首先需要利用數據驅動思想建立EC、EQ和出水水質參數的精確軟測量模型,作為優(yōu)化目標函數;然后,采用多目標優(yōu)化算法對O,5和NO,2進行動態(tài)尋優(yōu),獲取一組等同優(yōu)秀的Pareto解集;進一步,采用智能決策系統從Pareto解集選取偏好解作為優(yōu)化設定值;最后,設計具有控制精度高和穩(wěn)定性好的自適應控制器以跟蹤優(yōu)化設定值,從而實現多目標優(yōu)化控制。本文構建的污水處理過程多目標優(yōu)化控制整體架構如圖1所示,其中自組織預測模型和自組織控制器的設計過程見文獻[29-30],本文主要介紹多目標優(yōu)化算法的設計。

    BSM1模型自帶晴好天氣、陰雨天氣和暴雨天氣3種不同工況下入水流量和組分濃度的數據文件,采樣間隔為15 min[31],來源于某實際污水處理廠連續(xù)兩周的運行操作數據。圖2所示為晴好天氣下入水流量和入水S、B,H、NH濃度,能夠反映污水處理過程強非線性、不確定性嚴重和強耦合的非線性特點[31-32]。

    2 混合多目標骨干粒子群優(yōu)化算法

    下面首先闡述分解方法、引導者選擇、外部檔案維護、粒子位置更新和精英學習策略,然后給出HBBMOPSO算法的整體流程,最后給出污水生化處理過程優(yōu)化控制的整體流程。

    2.1 分解方法

    MOEA/D采用分解方法將多目標優(yōu)化問題分解為一系列單目標優(yōu)化子問題,然后采用進化算法協同解決這些單目標子問題,主要方法有:加權和法、Tchebycheff法和基于懲罰的邊界交集法(penalty-based boundary intersection,PBI)。研究表明,PBI方法具有更大的優(yōu)勢[27]。因此本文將PBI方法引入到BBMOPSO中,用于更新粒子的個體引導者,構建基于分解和支配相結合的混合多目標骨干粒子群算法(HBBMOPSO)。

    在PBI方法中[27],單目標優(yōu)化子問題定義為

    其中

    式中,>0為懲罰參數,=(1,2,…,w)T為預定義的均勻分布的權向量集。PBI方法如圖3所示,為理想參考點,為代表搜索方向的直線,為解()在上的投影,1為*和之間的距離,2為()和之間的垂直距離。

    在PBI方法中,通常取為固定值5.0[27-28]。然而,最近的研究表明[33],值越小收斂速度越快,值越大多樣性越好。文獻[33]采用值隨迭代次數線性增加策略,平衡了收斂性和多樣性。本文提出了一個新的值非線性調整策略,公式如下

    其中,的初始值為5,min1,max10;step=(max?min)/max。ΔEntropy為算法進化過程Pareto熵的變化量,由Pareto解映射到平行格坐標系統后得到,具體見文獻[21]。用于區(qū)分進化狀態(tài)的閾值由式(9)計算得到[21]

    (9)

    其中,為目標個數,為時刻檔案中解的個數。

    圖4給出了本文算法在第2天第1個優(yōu)化周期獲得的Pareto熵及其變化量,圖5給出了對應的值變化情況??梢钥闯觯軌螂S著進化狀態(tài)非線性調整,起到加快收斂速度和提高多樣性的目的。

    2.2 引導者選擇

    引導者選擇包括個體引導者pbest和全局引導者gbest的更新:pbest采用分解的方法從外部檔案選擇,并且考慮約束條件的影響;gbest隨機選自外部檔案中具有較大擁擠距離的前10%的解。

    pbest更新:對于當前粒子,根據解的可行性和權向量從外部檔案中為其選擇一個具有較低聚合函數值的解作為其個體引導者,具體流程見算法1。其中,對于粒子,其約束違背程度定義為

    其中,g(),=1,2,…,為粒子的個不等式約束條件;h(),=+1,+2,…,為粒子的?個等式約束條件。

    Algorithm 1:Update of pbest

    for= 1 to

    set pbest=1;

    for= 2 to ||

    if pbest&Aare feasible

    if(pbest|w,*) >(A|w,*)

    pbest=A;

    end if

    end if

    if pbestis feasible&Ais infeasible

    pbestremain;

    end if

    if pbestis infeasible&Ais feasible

    pbest=A;

    end if

    if pbest&Aare infeasible

    pbest=arg min cv(pbest,A);

    end if

    end for

    end for

    return pbest

    2.3 粒子位置更新

    粒子位置更新公式為[24]

    其中,μ() = [×pbest()+(1?)×gbest()]/2 和σ()=|pbest()?gbest()|,為[0,1]之間均勻分布的隨機數。采用pbest和gbest的隨機權重和作為高斯采樣的中心,擴大了粒子的搜索范圍;以50%的交叉概率繼承全局引導者gbest的信息,算法更加注重對gbest的開發(fā),提高了收斂速度。

    2.4 外部檔案維護

    HBBMOPSO算法采用有限容量的外部檔案來存儲非支配解,需要采用合適的外部檔案維護算法來提高解的多樣性和收斂性。目前常用的有自適應網格法[20]和擁擠距離法[26],由于擁擠距離法無須設置額外參數,更加適合工程應用,本文采用擁擠距離法[26,33]來維護外部檔案,具體流程見算法2。

    假設為外部檔案,為產生的新種群,a為外部檔案的最大容量。在ArchiveUpdate算法中,CheckDominance(S,A)函數用于評價SA之間的Pareto支配關系。當S支配A時,該函數返回1;否則,當A支配SSA互不支配時,返回?1;最終,算法將返回更新后的外部檔案。

    Algorithm 2:Update of External Archive

    for= 1 to ||

    for= 1 to ||

    flag = CheckDominance(S,A);

    if flag == 1

    markAas a dominated solution;

    else if flag == ?1

    markSas a dominated solution; break;

    end if

    end for

    delete the marked dominated solutions from;

    if Sis not marked as a dominated solution

    addSto;

    if || >N

    compute the crowding distance;

    delete the most crowded one;

    end if

    end if

    end for

    return

    2.5 精英學習策略

    為了克服粒子快速收斂特性帶來的易陷入局部Pareto前沿的缺陷,MOPSO一般都帶有粒子變異操作。但是,粒子變異算子難以彌補粒子搜索潛在的脆弱性[34],對于復雜的實際問題,種群容易陷入局部Pareto前沿。因此,設計了精英學習策略(elitism learning strategy,ELS)實現外部檔案中非支配解的進化搜索[34-35],具體流程見算法3。

    在算法3中,模擬二進制交叉(simulated binary crossover,SBX)操作用于精英解之間交換有用的基因片段,多項式變異(polynomial mutation,PM)操作通過產生一個小的擾動以對局部區(qū)域進行搜索。首先從外部檔案中選擇一定數目適應度值fitness(a,)較大的非支配解形成精英子集,本文的數目取||的一半。對于外部檔案中每一個解A(=1,2,…, ||),在[1, ||]隨機產生一個整數。然后,對AE構成的父代解執(zhí)行SBX操作產生子代解。最終從子代解1和2中隨機選取一個執(zhí)行PM操作,形成新的解S

    Algorithm 3:Elitism Learning

    for= 1 to ||

    generate a random integerin [1, ||];

    {1,2} = SBX(A,E)

    generate a random integerin [1, 2];

    S= PM(C);

    end for

    return

    2.6 HBBMOPSO算法整體流程

    下面給出HBBMOPSO算法的整體流程,首先對粒子位置和參考點進行初始化,然后算法進入主循環(huán)進行迭代,直至滿足停止條件,具體見算法4。

    Algorithm 4: HBBMOPSO

    initialize ev=0,=null, and={1,2,…,p};

    generatewell-distributed weighted vectors;

    for= 1 to

    randomly generate the positionxofp;

    evaluate the objective values ofp;

    set pbest=pand initialize*;

    end for

    update external archive (Algorithm 2);

    ev= ev +;

    while ev≤ max_ev

    environmental detecting using ΔEntropy

    adaptively adjustusing Eq. (8)

    for=1 to

    select gbest(Algorithm 1);

    update the positionxofpusing Eq. (11);

    evaluate the objective values forp;

    update pbest;

    update the reference point* in Eq. (7);

    end for

    ev= ev +;

    update external archive (Algorithm 2);

    perform ELS (Algorithm 3);

    evaluate new solutions’ objective values;

    update the reference point* in Eq. (7);

    update external archive (Algorithm 2);

    ev= ev +;

    end while

    output

    2.7 優(yōu)化控制整體流程

    污水處理過程多目標優(yōu)化控制整體流程如下。

    (1)利用BSM1基準離線產生500組建模數據,輸入為O,5、NO,2和入水水質參數,輸出為EC、EQ和出水水質參數;利用SOFNN[29]預測模型建立優(yōu)化算法的目標函數。

    其中,=[1,2,…,x]為SOFNN的輸入向量,為輸入變量數;=[1j,2j,…,c]和=[1j,2j,…,δ]分別為規(guī)則層第個神經元對應的隸屬函數層神經元的高斯函數中心和寬度;φ為規(guī)則層第個神經元歸一化后的輸出;w為規(guī)則層第個神經元與輸出層第個神經元之間的連接權;=1,2,…,,為規(guī)則層神經元數;=1,2,…,,為輸出變量數。SOFNN采用二階算法優(yōu)化中心、寬度和權值,采用奇異值分解優(yōu)化網絡結構,構建了結構緊湊、學習速度快和預測精度高的網絡模型[29]。

    (2)利用所提的HBBMOPSO算法(Algorithm 4)對EC、EQ和出水水質參數的目標函數進行優(yōu)化,得到一組Pareto最優(yōu)解集。

    (3)利用模糊隸屬函數法設計智能決策系統,從包含眾多等同優(yōu)秀解的Pareto解集中選取一個偏好解,以確定當前優(yōu)化周期內O,5和NO,2的底層跟蹤控制設定值。

    對于第個目標函數F,外部儲備集中非支配解的滿意度定義為

    其中,max和min分別是第個目標函數F的最大值和最小值。歸一化的非支配解的滿意度為

    (14)

    其中,為目標個數,|A|為外部儲備集元素個數。選取μ的最大值對應的解為偏好解。

    (4)利用基于互信息的自組織模糊神經網絡控制器對[O,sp,NO,sp]優(yōu)化設定值進行跟蹤控制,控制器在復雜工況的自適應能力有一定程度提高[30]。若BSM1連續(xù)14天數據仿真完畢,則結束;否則,返回步驟(2)進行下一個優(yōu)化周期的優(yōu)化控制。

    3 實驗驗證

    3.1 參數設置

    所有實驗均在BSM1基準仿真平臺上進行,控制器采樣周期為45 s,優(yōu)化周期為2 h,O,5和NO,2的優(yōu)化設定值范圍分別為0.4~3和0.2~2 mg·L?1。對于HBBMOPSO算法,種群規(guī)模為50,外部檔案規(guī)模a為50,目標函數最大評價次數ev為10000(即最大迭代次數max=200),交叉概率c為0.9,變異概率m為1/,SBX和PM的分布性指數c和m均為20。

    3.2 結果和分析

    圖6給出了分別利用SOFNN和FNN預測模型得到的能耗EC和出水水質EQ的預測結果。可以看出,SOFNN的預測精度要比FNN更高,為多目標優(yōu)化算法提供了更精確的目標函數。圖7給出了第2 d第1個優(yōu)化周期HBBMOPSO算法在固定、線性遞增和自適應3種情況下的逼近Pareto前沿??梢钥闯?,自適應獲取的Pareto解收斂性和多樣性更好。圖8給出了HBBMOPSO與NSGAII、dMOPSO、BBMOPSO算法獲得的第2 d第1個優(yōu)化周期的Pareto最優(yōu)解集以及決策系統決策出的偏好解。從圖7可以看出,本文算法獲取的Pareto前沿的收斂性和分布性更好,獲取的優(yōu)化設定值具有更低的EC和EQ。

    圖9給出了7 d的設定值優(yōu)化情況和底層SOFNN自組織控制器跟蹤效果。從圖9可以看出,在優(yōu)化控制過程中,O,sp和NO,sp隨著入水條件和組分濃度的變化進行動態(tài)調節(jié),實現了出水水質和能耗這對沖突目標之間的動態(tài)平衡;同時,SOFNN控制器的控制效果要明顯優(yōu)于PID控制,尤其是硝態(tài)氮的控制精度有明顯的提升,實現了復雜工況下的快速高精度的跟蹤控制。

    圖10給出了多目標優(yōu)化控制過程出水水質參數變化情況。從圖10可以看出,相比于PID恒定值控制,優(yōu)化控制方法下,NH濃度沒有觀察到明顯上升趨勢,tot濃度總體呈現下降趨勢。另外,BOD5、COD和TSS基本維持不變,表明優(yōu)化控制在有效降低能耗的同時并沒有引起出水水質參數超標而造成罰款上升的情況。從污水處理過程機理分析可知,tot和NH是一對具有相互競爭關系的出水水質參數,通過多目標優(yōu)化算法能夠取得這兩者之間的最佳平衡。

    表1給出了晴好天氣工況下不同優(yōu)化控制策略獲得的AE、PE、EC和EQ的具體情況。其中,influent表示平均入水水質參數;PID表示恒定值控制,即O,5和NO,2分別取2和1 mg·L?1;SOOC表示單目標優(yōu)化控制;Hopfield表示采用文獻[9]介紹的Hopfield神經網絡實現的多目標優(yōu)化控制;DDAOC表示采用文獻[11]介紹的自適應動態(tài)規(guī)劃方法實現的多目標優(yōu)化控制;NSGAII表示采用文獻[18]介紹的神經網絡和多目標優(yōu)化算法相結合的優(yōu)化控制;dMOPSO表示采用文獻[28]介紹的基于分解的多目標優(yōu)化算法,BBMOPSO表示采用文獻[24]介紹的多目標骨干粒子群算法,這兩者的參數均按照原文設置。從表1可以看出,相比PID恒定值控制,在晴好天氣工況下,HBBMOPSO方法的EC下降9.60%,并且EQ上升幅度較小。實驗結果表明HBBMOPSO方法不僅能夠保證出水水質參數達標,而且能夠有效降低能耗。同時,與dMOPSO和BBMOPSO相比,優(yōu)化效果有較大程度的提升,表明了帶有自適應懲罰參數的混合多目標骨干粒子群算法能夠獲取更高質量的解。此外,NSGAII方法的EQ上升幅度較大,主要是由于出水tot較大引起。單目標優(yōu)化SOOC方法的EC下降較小,節(jié)能效果不太明顯。對表進一步分析,可以看出,HBBMOPSO方法下的NH略有上升,tot有較大程度下降,其他水質參數變化幅度不大,這與圖7水質參數變化情況相吻合。需要指出的是,本文所有對比方法均能夠保證平均出水水質參數符合達標限定值。從出水污染物去除率的角度分析,采用HBBMOPSO優(yōu)化控制方法,相比入水水質參數,出水NH、tot、BOD5、COD和TSS的去除率分別達到了92.30%、72.61%、96.29%、71.66%和93.74%。

    表1 晴好天氣不同優(yōu)化控制策略的能耗和出水水質比較

    ① Results are listed in original papers.

    為了進一步驗證HBBMOPSO優(yōu)化方法在復雜天氣工況下的適應能力,利用BSM1基準模型中的陰雨和暴雨兩種工況的數據文件進行了實驗驗證,表2給出了不同優(yōu)化控制策略的對比結果。從表2可以看出,相比PID恒定值控制,在陰雨和暴雨天氣工況下,HBBMOPSO的EC分別下降了9.29%和8.77%,并且EQ沒有明顯上升。相比NSGAII、dMOPSO和BBMOPSO,HBBMOPSO優(yōu)化方法的能耗更低,表明HBBMOPSO在復雜工況下也能夠取得良好的節(jié)能效果。

    表2 陰雨和暴雨天氣不同優(yōu)化控制策略的EC和EQ比較

    ①Results are listed in original papers.

    4 結 論

    針對活性污泥法污水處理過程能耗和罰款最低的多目標優(yōu)化問題,提出了混合多目標骨干粒子群優(yōu)化(HBBMOPSO)算法。該方法將Pareto支配和分解策略相結合,利用分解策略選取個體引導者,利用擁擠距離和Pareto支配關系維護外部檔案和選取全局引導者。同時,設計了分解方法中懲罰參數的自適應調整策略,提高了算法收斂速度和非支配解的質量。此外,設計了外部檔案精英學習策略,提高了種群跳出局部Pareto前沿的能力。將HBBMOPSO算法與自組織模糊神經網絡預測模型、自組織控制器相結合,實現了污水處理過程多目標優(yōu)化運行控制。

    與單目標優(yōu)化、NSGAII、dMOPSO和BBMOPSO等方法相比,HBBMOPSO方法的尋優(yōu)能力更強,克服了傳統粒子群優(yōu)化算法對飛行參數敏感的不足,獲取的Pareto最優(yōu)解的收斂性和多樣性均有較大程度的改善,為智能決策系統提供了一組質量更高的候選解。利用智能決策系統選取的偏好解在保證出水水質參數達標的前提下,取得了更低的能耗。此外,實驗結果還表明,自組織模糊神經網絡構建的能耗和出水水質優(yōu)化目標模型具有較高的建模精度,自組織控制器能夠較快地跟蹤優(yōu)化設定值,具有較好的穩(wěn)定性、魯棒性和控制精度,能夠滿足多目標優(yōu)化控制的需求。

    下一步的研究工作是充分挖掘污水處理過程蘊含的知識信息,將其與過程數據相結合,構建基于知識和數據混合驅動的污水生化處理過程多目標優(yōu)化運行控制策略。

    References

    [1] WAN J F, GU J, ZHAO Q,. COD capture: a feasible option towards energy self-sufficient domestic wastewater treatment[J]. Scientific Reports, 2016, 6(4): 1-9.

    [2] OTURAN M A, AARON J J. Advanced oxidation processes in water/wastewater treatment: principles and applications. A review[J]. Critical Reviews in Environmental Science and Technology, 2014, 44(23): 2577-2641.

    [3] SANTíN I, PEDRET C, VILANOVA R,. Removing violations of the effluent pollution in a wastewater treatment process[J]. Chemical Engineering Journal, 2015, 279(11): 207-219.

    [4] JUDD S J. The status of industrial and municipal effluent treatment with membrane bioreactor technology[J]. Chemical Engineering Journal, 2016, 305(12): 37-45.

    [5] ?MAND L, CARLSSON B. Optimal aeration control in a nitrifying activated sludge process[J]. Water Research, 2012, 46(7): 2101-2110.

    [6] HREIZ R, LATIFI M A, ROCHE N. Optimal design and operation of activated sludge processes: state-of-the-art[J]. Chemical Engineering Journal, 2015, 281(12): 900-920.

    [7] OSTACE G S, BAEZA J A, GUERRERO J,. Development and economic assessment of different WWTP control strategies for optimal simultaneous removal of carbon, nitrogen and phosphorus[J]. Computers and Chemical Engineering, 2013, 53(6): 164-177.

    [8] SANTIN I, PEDRET C, VILANOVA R. Applying variable dissolved oxygen set point in a two level hierarchical control structure to a wastewater treatment process[J]. Journal of Process Control, 2015, 28(4): 40-55.

    [9] 韓廣, 喬俊飛, 韓紅桂, 等. 基于Hopfield神經網絡的污水處理過程優(yōu)化控制[J]. 控制與決策, 2014, 29(11): 2085-2088. HAN G, QIAO J F, HAN H G,. Optimal control for wastewater treatment process based on Hopfield neural network[J]. Control and Decision, 2014, 29(11): 2085-2088.

    [10] MACHADO V C, GABRIEL D, LAFUENTE J,. Cost and effluent quality controllers design based on the relative gain array for a nutrient removal WWTP[J]. Water Research, 2009, 43(20): 5129-5141.

    [11] QIAO J F, BO Y C, CHAI W,. Adaptive optimal control for a wastewater treatment plant based on a data-driven method[J]. Water Science and Technology, 2013, 67(10): 2314-2320.

    [12] DAI H L, CHEN W L, LU X W. The application of multi-objective optimization method for activated sludge process: a review[J]. Water Science and Technology, 2016, 73(2): 223-235.

    [13] HAKANEN J, SAHLSTEDT K, MIETTINEN K. Wastewater treatment plant design and operation under multiple conflicting objective functions[J]. Environmental Modelling and Software, 2013, 46(4): 240-249.

    [14] SWEETAPPLE C, FU G T, BUTLER D. Multi-objective optimisation of wastewater treatment plant control to reduce greenhouse gas emissions[J]. Water Research, 2014, 55(2): 52-62.

    [15] HREIZ R, ROCHE N, BENYAHIA B,. Multi-objective optimal control of small-size wastewater treatment plants[J]. Chemical Engineering Research and Design, 2015, 102(7): 345-353.

    [16] CHEN W L, LU X W, YAO C H. Optimal strategies evaluated by multi-objective optimization method for improving the performance of a novel cycle operating activated sludge process[J]. Chemical Engineering Journal, 2015, 260(9): 492-502.

    [17] ZHANG R, XIE W M, YU H Q,. Optimizing municipal wastewater treatment plants using an improved multi-objective optimization method[J]. Bioresource Technology, 2014, 157(2): 161-165.

    [18] QIAO J F, ZHANG W. Dynamic multi-objective optimization control for wastewater treatment process[J]. Neural Computing and Applications, 2016, 28(10): 1-11.

    [19] KENNEDY J, EBERHART R C. Particle swarm optimization[C]// Proceeding of IEEE International Conference on Neural Networks. 1995: 1942-1948.

    [20] COELLO C A C, PULIDO G T, LECHUGA M S. Handling multiple objectives with particle swarm optimization[J]. IEEE Transactions on Evolutionary Computation, 2004, 8(3): 256-279.

    [21] HU W, YEN G G. Adaptive multiobjective particle swarm optimization based on parallel cell coordinate system[J]. IEEE Transactions on Evolutionary Computation, 2015, 19(1): 1-18.

    [22] TRIPATHI P K, BANDYOPADHYAY S, PAL S K. Multi-objective particle swarm optimization with time variant inertia and acceleration coefficients[J]. Information Sciences, 2007, 177(22): 5033-5049.

    [23] KENNEDY J. Bare bones particle swarms[C]// Proceedings of IEEE Swarm Intelligence Symposium. 2003: 80-87.

    [24] ZHANG Y, GONG D W, DING Z H. A bare-bones multi-objective particle swarm optimization algorithm for environmental/economic dispatch[J]. Information Sciences, 2012, 192(6): 213-227.

    [25] ZHANG Y, GONG D W, GENG N,. Hybrid bare-bones PSO for dynamic economic dispatch with valve-point effects[J]. Applied Soft Computing, 2014, 18(5): 248-260.

    [26] DEB K, PRATAP A, AGARWAL S,. A fast and elitist multiobjective genetic algorithm: NSGA-II[J]. IEEE Transactions on Evolutionary Computation, 2002, 6(2): 182-197.

    [27] ZHANG Q F, LI H. MOEA/D: a multiobjective evolutionary algorithm based on decomposition[J]. IEEE Transactions on Evolutionary Computation, 2007, 11(6): 712-731.

    [28] ZAPOTECAS MARTíNEZ S, COELLO COELLO C A. A multi-objective particle swarm optimizer based on decomposition[C]// Proceedings of the 13th Annual Conference on Genetic and Evolutionary Computation. 2011: 69-76.

    [29] 喬俊飛, 周紅標. 基于自組織模糊神經網絡的出水總磷預測[J]. 控制理論與應用, 2017, 34(2): 224-232. QIAO J F, ZHOU H B. Prediction of effluent total phosphorus based on self-organizing fuzzy neural network[J]. Control Theory and Applications, 2017, 34(2): 224-232.

    [30] 周紅標. 基于自組織模糊神經網絡的污水處理過程溶解氧控制[J]. 化工學報, 2017, 68(4): 1516-1524. ZHOU H B. Dissolved oxygen control of the wastewater treatment process using self-organizing fuzzy neural network[J]. CIESC Journal, 2017, 68(4): 1516-1524.

    [31] JEPPSSON U, PONS M N. The COST benchmark simulation model-current state and future perspective[J]. Control Engineering Practice, 2004, 12(3): 299-304.

    [32] HAN H G, QIAN H H, QIAO J F. Nonlinear multiobjective model-predictive control scheme for wastewater treatment process[J]. Journal of Process Control, 2014, 24(3): 47-59.

    [33] YANG S X, JIANG S Y, JIANG Y. Improving the multiobjective evolutionary algorithm based on decomposition with new penalty schemes[J]. Soft Computing, 2016, 21(2): 1-15.

    [34] LIN Q Z, LI J Q, DU Z H,. A novel multi-objective particle swarm optimization with multiple search strategies[J]. European Journal of Operational Research, 2015, 247(3): 732-744.

    [35] LIN Q Z, CHEN J Y, ZHAN Z H,. A hybrid evolutionary immune algorithm for multiobjective optimization problems[J]. IEEE Transactions on Evolutionary Computation, 2016, 20(5): 711-729.

    Optimal control of wastewater treatment process using hybrid multi-objective barebones particle swarm optimization algorithm

    ZHOU Hongbiao1,2,3, QIAO Junfei1,2

    (1Faculty of Information Technology,Beijing University of Technology, Beijing 100124, China;2Beijing Key Laboratory of Computational Intelligence and Intelligent System, Beijing 100124, China;3Faculty of Automation, Huaiyin Institute of Technology, Huai’an 223003, Jiangsu, China)

    Through analysis of biological wastewater treatment process (WWTP), a multi-objective optimal control strategy was developed with targets of minimizing both energy consumption and amercement. A hybrid multi-objective barebones particle swarm optimization (HBBMOPSO) algorithm based on Pareto dominance and decomposition was proposed to improve convergence and diversity of optimized set of Pareto solutions. In HBBMOPSO, selection of personal leaders was determined from self-adaptive penalty factor decomposition while maintenance of external dossiers and selection of global leaders were determined from dominance and crowded distance. Furthermore, elitism learning strategy was adopted to facilitate particle escaping from local Pareto fronts. Finally, HBBMOPSO was combined with self-organizing fuzzy nerve network modeler and controller to realize dynamic optimization, intelligent decision, and background monitoring on dissolved oxygen and nitrate nitrogen in biological WWTP. Experimental study on international standardized simulator platform BSM1 showed that HBBMOPSO method can effectively reduce energy consumption under the premise of ensuring effluent to meet quality standard.

    wastewater; optimization; process control; particle swarm; decomposition

    10.11949/j.issn.0438-1157.20170583

    TP 273

    A

    0438—1157(2017)09—3511—11

    2017-05-09收到初稿,2017-06-05收到修改稿。

    喬俊飛。

    周紅標(1980—),男,博士研究生。

    國家自然科學基金重點項目(61533002)。

    2017-05-09.

    Prof. QIAO Junfei, junfeiq@bjut.edu.cn

    supported by the State Key Program of National Natural Science Foundation of China (61533002).

    猜你喜歡
    引導者設定值污水處理
    我國鄉(xiāng)鎮(zhèn)污水處理模式的探討
    冷軋鍍鋅光整機延伸率控制模式的解析與優(yōu)化
    湖北農機化(2021年7期)2021-12-07 17:18:46
    太原市61個村要建污水處理設施嚴禁直排入河
    引導者 傳播者 擔當者——新年寄語《人大建設》
    人大建設(2019年1期)2019-07-13 06:00:30
    涌浪機在污水處理中的應用
    西氣東輸二線人機交互界面的防呆系統
    自動化儀表(2015年5期)2015-06-15 19:01:34
    目標設定值傳感器的原理與檢修
    對工頻耐壓試驗跳閘電流設定值問題的探討
    三亞加快污水處理工程建設
    數學教師如何當好“引導者”
    亚洲成人av在线免费| 久久久久久人人人人人| 咕卡用的链子| 丰满乱子伦码专区| 一区二区av电影网| 婷婷色综合大香蕉| 国产伦人伦偷精品视频| a级片在线免费高清观看视频| 人人妻,人人澡人人爽秒播 | 99香蕉大伊视频| 免费高清在线观看视频在线观看| 老汉色av国产亚洲站长工具| 少妇精品久久久久久久| a 毛片基地| 精品卡一卡二卡四卡免费| 久久久久国产一级毛片高清牌| 母亲3免费完整高清在线观看| 欧美黑人欧美精品刺激| 精品亚洲成国产av| 九九爱精品视频在线观看| 精品少妇一区二区三区视频日本电影 | a级毛片在线看网站| 99精国产麻豆久久婷婷| 国产欧美亚洲国产| 欧美黑人精品巨大| 午夜激情av网站| 一本久久精品| 亚洲精品av麻豆狂野| 国产男人的电影天堂91| 一本—道久久a久久精品蜜桃钙片| 亚洲天堂av无毛| 在线观看免费日韩欧美大片| 亚洲欧美一区二区三区黑人| 人体艺术视频欧美日本| 日韩熟女老妇一区二区性免费视频| 日本黄色日本黄色录像| 国产av一区二区精品久久| 国产免费又黄又爽又色| 成年av动漫网址| 国产精品国产三级国产专区5o| 91精品伊人久久大香线蕉| 国产日韩一区二区三区精品不卡| 中文字幕色久视频| 久久久精品国产亚洲av高清涩受| 国产精品成人在线| 欧美在线一区亚洲| 在线天堂最新版资源| 看十八女毛片水多多多| 亚洲,欧美,日韩| 亚洲精品久久成人aⅴ小说| 国产爽快片一区二区三区| 超碰97精品在线观看| 丝袜在线中文字幕| 久久久久久免费高清国产稀缺| 国产精品免费视频内射| 侵犯人妻中文字幕一二三四区| 最近2019中文字幕mv第一页| 老司机亚洲免费影院| av不卡在线播放| 五月天丁香电影| 超碰成人久久| 老司机深夜福利视频在线观看 | 亚洲在久久综合| 十八禁网站网址无遮挡| 一本一本久久a久久精品综合妖精| 午夜福利一区二区在线看| 80岁老熟妇乱子伦牲交| 国产午夜精品一二区理论片| av视频免费观看在线观看| avwww免费| avwww免费| 蜜桃在线观看..| 国产又色又爽无遮挡免| 国产亚洲av片在线观看秒播厂| 嫩草影视91久久| 中文字幕制服av| 伊人久久大香线蕉亚洲五| 国产无遮挡羞羞视频在线观看| 久久久久久免费高清国产稀缺| 国产精品蜜桃在线观看| 丝瓜视频免费看黄片| 亚洲精品国产av蜜桃| 在线亚洲精品国产二区图片欧美| 悠悠久久av| 天天躁夜夜躁狠狠躁躁| 丝瓜视频免费看黄片| 国产成人系列免费观看| 久久精品aⅴ一区二区三区四区| 欧美在线黄色| 搡老岳熟女国产| 亚洲第一av免费看| 国产日韩欧美在线精品| 亚洲成人免费av在线播放| tube8黄色片| 日韩精品免费视频一区二区三区| 国产精品国产av在线观看| 国产成人午夜福利电影在线观看| 一本一本久久a久久精品综合妖精| 看免费成人av毛片| 精品酒店卫生间| 色网站视频免费| 午夜福利在线免费观看网站| 国产高清不卡午夜福利| 免费高清在线观看日韩| 免费观看av网站的网址| 久久免费观看电影| 午夜av观看不卡| 国产精品秋霞免费鲁丝片| 日韩中文字幕视频在线看片| 高清av免费在线| 9色porny在线观看| 啦啦啦在线免费观看视频4| 午夜日韩欧美国产| 丰满少妇做爰视频| 免费在线观看视频国产中文字幕亚洲 | 国产精品免费视频内射| 青春草国产在线视频| 久久久精品区二区三区| 国产精品三级大全| 亚洲一级一片aⅴ在线观看| 国产97色在线日韩免费| xxx大片免费视频| 伦理电影免费视频| 国产成人av激情在线播放| 最近中文字幕高清免费大全6| 亚洲国产中文字幕在线视频| 国产精品三级大全| 91老司机精品| 久久 成人 亚洲| av国产久精品久网站免费入址| 精品久久久精品久久久| 男人爽女人下面视频在线观看| 最近的中文字幕免费完整| 另类精品久久| 九九爱精品视频在线观看| 国产xxxxx性猛交| 国产午夜精品一二区理论片| 亚洲国产日韩一区二区| 国产有黄有色有爽视频| 男女边摸边吃奶| 69精品国产乱码久久久| 精品一区二区三区四区五区乱码 | 亚洲在久久综合| 久久久久久久久久久免费av| 男女之事视频高清在线观看 | 亚洲精品一区蜜桃| 国产又爽黄色视频| 免费在线观看完整版高清| 蜜桃在线观看..| 国产精品.久久久| 狠狠婷婷综合久久久久久88av| 免费高清在线观看视频在线观看| 亚洲国产av影院在线观看| 一边亲一边摸免费视频| 亚洲情色 制服丝袜| 卡戴珊不雅视频在线播放| 一本久久精品| 99热全是精品| 极品少妇高潮喷水抽搐| 国产av一区二区精品久久| 深夜精品福利| 午夜福利免费观看在线| 男的添女的下面高潮视频| 最黄视频免费看| 成人国产av品久久久| 日本午夜av视频| 亚洲av成人精品一二三区| 777米奇影视久久| 两个人看的免费小视频| 热re99久久国产66热| 日本欧美国产在线视频| 老司机深夜福利视频在线观看 | 日韩一本色道免费dvd| 日韩av不卡免费在线播放| 亚洲av国产av综合av卡| 婷婷色综合大香蕉| 美女扒开内裤让男人捅视频| 日日爽夜夜爽网站| 中文乱码字字幕精品一区二区三区| 亚洲五月色婷婷综合| 亚洲国产成人一精品久久久| 精品国产一区二区三区四区第35| 国产探花极品一区二区| 国产高清不卡午夜福利| 一级,二级,三级黄色视频| 少妇人妻久久综合中文| av网站在线播放免费| 久久国产亚洲av麻豆专区| 美女国产高潮福利片在线看| 国产欧美日韩一区二区三区在线| 成年动漫av网址| 欧美人与善性xxx| 最近中文字幕2019免费版| 在线观看免费日韩欧美大片| 亚洲国产欧美在线一区| 两个人看的免费小视频| 一本大道久久a久久精品| 曰老女人黄片| 亚洲熟女毛片儿| 97在线人人人人妻| 赤兔流量卡办理| 国产在线免费精品| 久久天躁狠狠躁夜夜2o2o | av网站免费在线观看视频| 国产精品一国产av| e午夜精品久久久久久久| 欧美日韩视频高清一区二区三区二| 高清不卡的av网站| 欧美日韩精品网址| 亚洲图色成人| 高清视频免费观看一区二区| 免费女性裸体啪啪无遮挡网站| 欧美精品一区二区免费开放| 极品人妻少妇av视频| 亚洲精品国产一区二区精华液| 免费观看性生交大片5| av在线播放精品| 国产精品久久久久成人av| 麻豆av在线久日| 成人手机av| 中文精品一卡2卡3卡4更新| 国产一区亚洲一区在线观看| xxxhd国产人妻xxx| 亚洲自偷自拍图片 自拍| 日韩中文字幕欧美一区二区 | 国产亚洲最大av| 欧美另类一区| 观看美女的网站| av电影中文网址| 老汉色∧v一级毛片| 国产亚洲av片在线观看秒播厂| 久久精品国产亚洲av涩爱| 婷婷色麻豆天堂久久| 精品一区二区免费观看| 多毛熟女@视频| 我的亚洲天堂| 超碰成人久久| 国产在线一区二区三区精| 一本一本久久a久久精品综合妖精| 久久99一区二区三区| 别揉我奶头~嗯~啊~动态视频 | 欧美变态另类bdsm刘玥| 国产在线免费精品| 亚洲av中文av极速乱| 中文字幕色久视频| 我的亚洲天堂| 一边摸一边抽搐一进一出视频| 国产日韩欧美亚洲二区| 日韩,欧美,国产一区二区三区| 日韩大码丰满熟妇| 国产欧美日韩综合在线一区二区| 可以免费在线观看a视频的电影网站 | 欧美国产精品一级二级三级| 又黄又粗又硬又大视频| 久久久久久久国产电影| 亚洲国产日韩一区二区| 久久性视频一级片| 亚洲国产欧美在线一区| 中文字幕高清在线视频| 999久久久国产精品视频| 80岁老熟妇乱子伦牲交| 美女大奶头黄色视频| 久久这里只有精品19| 一二三四中文在线观看免费高清| 国产 精品1| 欧美精品高潮呻吟av久久| 黄频高清免费视频| 少妇人妻 视频| 丁香六月欧美| 97人妻天天添夜夜摸| 九九爱精品视频在线观看| 高清视频免费观看一区二区| 日韩av不卡免费在线播放| 国产一区亚洲一区在线观看| 亚洲国产最新在线播放| 欧美 日韩 精品 国产| 欧美人与善性xxx| 亚洲精品,欧美精品| 成人影院久久| av女优亚洲男人天堂| 涩涩av久久男人的天堂| 国产男女超爽视频在线观看| 哪个播放器可以免费观看大片| 中国国产av一级| 少妇人妻精品综合一区二区| 日韩欧美一区视频在线观看| 天天操日日干夜夜撸| 777米奇影视久久| 亚洲精品国产av成人精品| 免费黄频网站在线观看国产| 99精国产麻豆久久婷婷| 成人国产麻豆网| av卡一久久| 精品亚洲乱码少妇综合久久| 欧美日本中文国产一区发布| 哪个播放器可以免费观看大片| 国产精品女同一区二区软件| 国产精品国产三级国产专区5o| 亚洲成人国产一区在线观看 | 午夜福利网站1000一区二区三区| 亚洲精品乱久久久久久| 亚洲成人手机| 国产亚洲av高清不卡| 成年动漫av网址| 电影成人av| 久久久精品94久久精品| 国产精品一区二区精品视频观看| 91成人精品电影| 亚洲第一区二区三区不卡| 成年女人毛片免费观看观看9 | 国产淫语在线视频| 男女下面插进去视频免费观看| 一区二区av电影网| 国产成人a∨麻豆精品| 久久精品国产a三级三级三级| 亚洲欧美成人精品一区二区| 男女边摸边吃奶| 久久影院123| 下体分泌物呈黄色| 另类亚洲欧美激情| 久久99精品国语久久久| 亚洲国产欧美在线一区| 国产不卡av网站在线观看| 最近最新中文字幕大全免费视频 | av免费观看日本| 啦啦啦视频在线资源免费观看| 亚洲美女搞黄在线观看| 久久女婷五月综合色啪小说| 性少妇av在线| 精品卡一卡二卡四卡免费| 韩国精品一区二区三区| 亚洲男人天堂网一区| 男人操女人黄网站| 高清不卡的av网站| 黑人巨大精品欧美一区二区蜜桃| 黄色一级大片看看| av福利片在线| av在线播放精品| 亚洲国产精品国产精品| 大香蕉久久网| 亚洲免费av在线视频| 大片免费播放器 马上看| videos熟女内射| 国产伦人伦偷精品视频| 日韩av免费高清视频| 777久久人妻少妇嫩草av网站| 欧美变态另类bdsm刘玥| 免费人妻精品一区二区三区视频| 久久综合国产亚洲精品| 久久久久久人人人人人| 国产免费现黄频在线看| 菩萨蛮人人尽说江南好唐韦庄| 日日摸夜夜添夜夜爱| 亚洲伊人色综图| 亚洲av电影在线进入| 精品一区二区免费观看| 热99国产精品久久久久久7| 国产免费又黄又爽又色| 一级片免费观看大全| 国产精品一区二区在线不卡| 国产亚洲一区二区精品| 十八禁网站网址无遮挡| 久久久精品国产亚洲av高清涩受| 99久久综合免费| 涩涩av久久男人的天堂| 亚洲欧美一区二区三区黑人| 美女扒开内裤让男人捅视频| 99热全是精品| 亚洲精品美女久久av网站| 两性夫妻黄色片| 亚洲成人手机| 国产一区有黄有色的免费视频| 久久久久久久久久久久大奶| 黑人猛操日本美女一级片| 99热全是精品| 成年人免费黄色播放视频| 成人国产麻豆网| 亚洲自偷自拍图片 自拍| 久久久精品免费免费高清| 国产精品 国内视频| 亚洲精品国产色婷婷电影| 免费看av在线观看网站| 岛国毛片在线播放| 久久精品久久久久久噜噜老黄| 又黄又粗又硬又大视频| 午夜福利乱码中文字幕| 亚洲天堂av无毛| 亚洲国产最新在线播放| 日韩电影二区| 欧美精品一区二区大全| 国产亚洲欧美精品永久| 99久国产av精品国产电影| 秋霞伦理黄片| 黄片无遮挡物在线观看| 中文字幕高清在线视频| 波多野结衣av一区二区av| 亚洲国产欧美一区二区综合| 国产精品一区二区在线不卡| 久久久久久人人人人人| 欧美日韩一区二区视频在线观看视频在线| av网站在线播放免费| 久久精品久久久久久久性| 在线天堂中文资源库| 中国三级夫妇交换| 亚洲精品日本国产第一区| 国产色婷婷99| av在线观看视频网站免费| 嫩草影院入口| 亚洲精品久久成人aⅴ小说| 哪个播放器可以免费观看大片| 亚洲国产欧美网| 一区福利在线观看| 男人操女人黄网站| 美国免费a级毛片| 欧美变态另类bdsm刘玥| av福利片在线| 热99国产精品久久久久久7| 少妇被粗大的猛进出69影院| 如日韩欧美国产精品一区二区三区| 一二三四中文在线观看免费高清| 波多野结衣av一区二区av| 国产黄色免费在线视频| 男男h啪啪无遮挡| 国产色婷婷99| 欧美日韩视频精品一区| 色播在线永久视频| 秋霞伦理黄片| 成人黄色视频免费在线看| 中文字幕精品免费在线观看视频| 国产xxxxx性猛交| av免费观看日本| 一边摸一边抽搐一进一出视频| 欧美激情 高清一区二区三区| 亚洲欧美成人综合另类久久久| av电影中文网址| 欧美黑人欧美精品刺激| 最近最新中文字幕免费大全7| 777久久人妻少妇嫩草av网站| 中文字幕人妻丝袜一区二区 | 在线免费观看不下载黄p国产| 看免费成人av毛片| 欧美中文综合在线视频| 久久久精品区二区三区| 亚洲av电影在线进入| 成人漫画全彩无遮挡| 精品亚洲成a人片在线观看| 中文天堂在线官网| 一区二区三区四区激情视频| 日韩视频在线欧美| 日韩av不卡免费在线播放| 毛片一级片免费看久久久久| 日韩 欧美 亚洲 中文字幕| 国产福利在线免费观看视频| 国产亚洲最大av| 九九爱精品视频在线观看| 亚洲激情五月婷婷啪啪| av国产久精品久网站免费入址| 午夜福利乱码中文字幕| 精品视频人人做人人爽| 欧美精品一区二区免费开放| 午夜91福利影院| 午夜福利一区二区在线看| 久久人妻熟女aⅴ| 一级毛片 在线播放| 国产日韩欧美视频二区| 99国产综合亚洲精品| svipshipincom国产片| 十八禁高潮呻吟视频| 亚洲人成电影观看| 高清欧美精品videossex| 亚洲精品一二三| 亚洲一级一片aⅴ在线观看| 国产精品久久久久久精品古装| 久久久久久久久久久免费av| 午夜久久久在线观看| 国产深夜福利视频在线观看| 久久天躁狠狠躁夜夜2o2o | 久久久久精品国产欧美久久久 | 精品人妻熟女毛片av久久网站| 2018国产大陆天天弄谢| 熟妇人妻不卡中文字幕| 这个男人来自地球电影免费观看 | 国产探花极品一区二区| 高清av免费在线| 老熟女久久久| 免费人妻精品一区二区三区视频| 久久精品国产综合久久久| 欧美激情 高清一区二区三区| 少妇的丰满在线观看| 国产日韩欧美视频二区| 一级毛片黄色毛片免费观看视频| 综合色丁香网| 无遮挡黄片免费观看| 亚洲天堂av无毛| 日韩 亚洲 欧美在线| 在线天堂中文资源库| av天堂久久9| 久久国产精品男人的天堂亚洲| 制服诱惑二区| 黄色视频在线播放观看不卡| 丝袜美腿诱惑在线| 亚洲色图综合在线观看| 考比视频在线观看| 视频区图区小说| 性色av一级| 永久免费av网站大全| 男男h啪啪无遮挡| 亚洲,欧美精品.| av又黄又爽大尺度在线免费看| 久久久久久人人人人人| 午夜日本视频在线| 99久国产av精品国产电影| 久久国产精品大桥未久av| 91老司机精品| 成人午夜精彩视频在线观看| 久久久久精品久久久久真实原创| 热re99久久精品国产66热6| 高清黄色对白视频在线免费看| 亚洲精品av麻豆狂野| 亚洲欧美成人精品一区二区| 欧美日韩亚洲国产一区二区在线观看 | 国产日韩欧美视频二区| 永久免费av网站大全| 欧美乱码精品一区二区三区| 久久久国产精品麻豆| 日本爱情动作片www.在线观看| 久久久欧美国产精品| 大码成人一级视频| 伊人久久国产一区二区| 如日韩欧美国产精品一区二区三区| 精品少妇一区二区三区视频日本电影 | 岛国毛片在线播放| 国产一区二区 视频在线| kizo精华| 中文精品一卡2卡3卡4更新| 新久久久久国产一级毛片| 色精品久久人妻99蜜桃| 中文乱码字字幕精品一区二区三区| 国产av精品麻豆| 丁香六月欧美| www.av在线官网国产| 综合色丁香网| 男女床上黄色一级片免费看| av网站在线播放免费| 亚洲欧美色中文字幕在线| 亚洲免费av在线视频| 伊人久久国产一区二区| 中文字幕人妻熟女乱码| 国产精品久久久久久久久免| 午夜av观看不卡| 精品国产超薄肉色丝袜足j| 亚洲一区二区三区欧美精品| 国产片特级美女逼逼视频| 最近最新中文字幕大全免费视频 | 少妇猛男粗大的猛烈进出视频| 国产在线视频一区二区| 亚洲情色 制服丝袜| 1024香蕉在线观看| 中国国产av一级| 国产又色又爽无遮挡免| 欧美少妇被猛烈插入视频| 青春草视频在线免费观看| 成人三级做爰电影| 欧美精品高潮呻吟av久久| 国产成人欧美在线观看 | 操出白浆在线播放| 18禁动态无遮挡网站| 青草久久国产| 侵犯人妻中文字幕一二三四区| 亚洲七黄色美女视频| 亚洲国产最新在线播放| 免费高清在线观看日韩| 美女福利国产在线| 色网站视频免费| 亚洲欧美一区二区三区久久| 最近中文字幕2019免费版| 欧美乱码精品一区二区三区| 人人妻人人澡人人爽人人夜夜| 国产成人免费无遮挡视频| 美女高潮到喷水免费观看| 久久综合国产亚洲精品| 18禁国产床啪视频网站| 在线观看国产h片| 高清视频免费观看一区二区| 两个人看的免费小视频| 黄片小视频在线播放| 成人影院久久| 熟女av电影| 人人妻人人添人人爽欧美一区卜| 欧美人与善性xxx| 日本爱情动作片www.在线观看| 捣出白浆h1v1| 国精品久久久久久国模美| 国产欧美亚洲国产| 最近中文字幕高清免费大全6| 亚洲精品美女久久av网站| 一二三四中文在线观看免费高清| 最近中文字幕高清免费大全6| 亚洲人成77777在线视频| 国产av国产精品国产| 一本一本久久a久久精品综合妖精| 色视频在线一区二区三区| 国产在线视频一区二区| 久久国产亚洲av麻豆专区| 久久久久久久久久久免费av| 午夜福利影视在线免费观看| 色网站视频免费| 无限看片的www在线观看| 免费观看av网站的网址| 超碰成人久久| 最近最新中文字幕大全免费视频 | 婷婷色综合www| 大话2 男鬼变身卡| 九色亚洲精品在线播放| 男女边摸边吃奶|