• <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久久久久免费| 日韩熟女老妇一区二区性免费视频| 女人久久www免费人成看片| 国产亚洲精品久久久久5区| 满18在线观看网站| 美女大奶头黄色视频| 母亲3免费完整高清在线观看| 只有这里有精品99| 久久精品亚洲av国产电影网| 日本色播在线视频| 99热网站在线观看| 久久午夜综合久久蜜桃| 亚洲精品日本国产第一区| 丁香六月天网| 男女免费视频国产| 欧美中文综合在线视频| 亚洲精品成人av观看孕妇| 国产淫语在线视频| 色精品久久人妻99蜜桃| 国产亚洲欧美精品永久| 黄色视频在线播放观看不卡| 精品国产一区二区三区四区第35| 国产成人精品无人区| 中文字幕色久视频| 曰老女人黄片| 国产一区有黄有色的免费视频| 日韩熟女老妇一区二区性免费视频| 婷婷色综合大香蕉| 丝袜美足系列| 精品亚洲乱码少妇综合久久| 欧美日韩亚洲综合一区二区三区_| 人人妻人人澡人人看| 中文字幕av电影在线播放| 蜜桃国产av成人99| 性高湖久久久久久久久免费观看| 少妇人妻久久综合中文| 国产精品久久久久成人av| 亚洲精品国产av蜜桃| 国产一区二区三区综合在线观看| 欧美人与性动交α欧美精品济南到| 久久国产精品人妻蜜桃| 性色av一级| 午夜福利,免费看| 只有这里有精品99| 国产成人免费无遮挡视频| 悠悠久久av| 精品少妇黑人巨大在线播放| 久久精品久久久久久噜噜老黄| 脱女人内裤的视频| 欧美日韩亚洲高清精品| 一二三四社区在线视频社区8| 1024视频免费在线观看| 亚洲国产av新网站| 免费不卡黄色视频| 亚洲色图综合在线观看| 2021少妇久久久久久久久久久| 一区二区三区激情视频| 国产日韩欧美亚洲二区| 亚洲精品一二三| 国产成人精品久久久久久| 老鸭窝网址在线观看| 亚洲av国产av综合av卡| 婷婷丁香在线五月| 免费日韩欧美在线观看| av视频免费观看在线观看| 1024视频免费在线观看| 亚洲,一卡二卡三卡| 亚洲国产成人一精品久久久| 国产高清不卡午夜福利| 亚洲精品国产一区二区精华液| a 毛片基地| 成人午夜精彩视频在线观看| 免费看十八禁软件| 中文字幕最新亚洲高清| 中文字幕亚洲精品专区| 久久中文字幕一级| 国产亚洲午夜精品一区二区久久| 热re99久久精品国产66热6| 免费在线观看完整版高清| 一区二区三区精品91| 青青草视频在线视频观看| 亚洲精品中文字幕在线视频| 精品久久久久久久毛片微露脸 | 精品一区二区三区四区五区乱码 | 99re6热这里在线精品视频| 国产一区二区三区av在线| 亚洲精品久久午夜乱码| 男女午夜视频在线观看| 色精品久久人妻99蜜桃| 丰满迷人的少妇在线观看| 热re99久久精品国产66热6| 欧美日本中文国产一区发布| 免费不卡黄色视频| 国产精品免费视频内射| 免费在线观看完整版高清| 男人操女人黄网站| 最近手机中文字幕大全| 午夜免费男女啪啪视频观看| 亚洲美女黄色视频免费看| 一区二区三区四区激情视频| 91成人精品电影| 亚洲七黄色美女视频| 一级片'在线观看视频| 欧美精品亚洲一区二区| 亚洲国产日韩一区二区| 久久久久久亚洲精品国产蜜桃av| 人人妻人人爽人人添夜夜欢视频| 两性夫妻黄色片| 久久综合国产亚洲精品| 国产爽快片一区二区三区| 亚洲成色77777| 十分钟在线观看高清视频www| 热99国产精品久久久久久7| 在线观看免费视频网站a站| 18在线观看网站| av天堂久久9| 婷婷成人精品国产| 免费观看人在逋| 高清av免费在线| 秋霞在线观看毛片| 夫妻性生交免费视频一级片| 午夜久久久在线观看| 国产精品国产三级国产专区5o| 99国产精品免费福利视频| 免费在线观看黄色视频的| 国产免费一区二区三区四区乱码| 成年动漫av网址| 最新的欧美精品一区二区| 欧美日韩视频高清一区二区三区二| 中文字幕制服av| 国产av精品麻豆| 2018国产大陆天天弄谢| 精品国产超薄肉色丝袜足j| 大香蕉久久网| 啦啦啦在线观看免费高清www| 国产精品一区二区精品视频观看| 国产精品二区激情视频| xxx大片免费视频| 久久国产精品影院| 色网站视频免费| 国产男人的电影天堂91| 久久久久久久大尺度免费视频| 亚洲人成77777在线视频| 一级片'在线观看视频| 亚洲成人免费电影在线观看 | 一区二区三区乱码不卡18| 亚洲第一av免费看| 亚洲一卡2卡3卡4卡5卡精品中文| 国产精品熟女久久久久浪| 国产91精品成人一区二区三区 | 日韩大码丰满熟妇| 亚洲欧美一区二区三区久久| 免费看不卡的av| 美女主播在线视频| 国产精品一区二区精品视频观看| 黑丝袜美女国产一区| 一级片'在线观看视频| 亚洲精品在线美女| 满18在线观看网站| 精品高清国产在线一区| 国产午夜精品一二区理论片| 午夜老司机福利片| 女性被躁到高潮视频| 亚洲国产精品999| 人人妻人人澡人人看| 五月天丁香电影| 人人妻人人爽人人添夜夜欢视频| 女人精品久久久久毛片| 一二三四在线观看免费中文在| www.av在线官网国产| 亚洲成国产人片在线观看| 97精品久久久久久久久久精品| 日本黄色日本黄色录像| 久久久久久亚洲精品国产蜜桃av| 欧美日韩黄片免| 国产有黄有色有爽视频| 脱女人内裤的视频| 日本欧美视频一区| bbb黄色大片| 老司机影院成人| 亚洲,欧美,日韩| 七月丁香在线播放| 久久人人爽人人片av| 成人亚洲精品一区在线观看| 久久影院123| 丝瓜视频免费看黄片| 欧美在线一区亚洲| 啦啦啦在线免费观看视频4| 在线观看一区二区三区激情| 国产欧美日韩综合在线一区二区| 亚洲五月色婷婷综合| 亚洲av日韩精品久久久久久密 | videos熟女内射| 视频在线观看一区二区三区| 国语对白做爰xxxⅹ性视频网站| 另类亚洲欧美激情| 亚洲国产av新网站| 国产精品久久久久久精品古装| 欧美日韩亚洲综合一区二区三区_| 久久ye,这里只有精品| 国产在线一区二区三区精| 免费看十八禁软件| 亚洲欧美一区二区三区久久| 亚洲成人国产一区在线观看 | 午夜免费男女啪啪视频观看| 99热网站在线观看| 欧美 日韩 精品 国产| 亚洲欧美精品自产自拍| 国产1区2区3区精品| 精品国产一区二区三区四区第35| 亚洲一卡2卡3卡4卡5卡精品中文| 又紧又爽又黄一区二区| 波野结衣二区三区在线| 成人影院久久| 日韩电影二区| 女人精品久久久久毛片| 国产精品久久久久久精品古装| 国产成人精品久久久久久| 亚洲精品久久午夜乱码| 久久久欧美国产精品| 日韩制服骚丝袜av| 亚洲av综合色区一区| 满18在线观看网站| 精品久久久精品久久久| 五月开心婷婷网| 国产精品三级大全| 少妇精品久久久久久久| 日本av免费视频播放| 在线 av 中文字幕| 久久人妻福利社区极品人妻图片 | 美女脱内裤让男人舔精品视频| 欧美大码av| 91成人精品电影| 欧美日韩视频精品一区| 亚洲欧美精品综合一区二区三区| 欧美激情高清一区二区三区| 99国产精品99久久久久| 99国产精品免费福利视频| 午夜激情久久久久久久| 国产野战对白在线观看| 亚洲av成人精品一二三区| 亚洲精品第二区| 大陆偷拍与自拍| 久久国产精品人妻蜜桃| 国产av一区二区精品久久| 国产97色在线日韩免费| 嫩草影视91久久| 汤姆久久久久久久影院中文字幕| 国产片特级美女逼逼视频| 午夜福利免费观看在线| 啦啦啦中文免费视频观看日本| 午夜免费观看性视频| 亚洲国产av影院在线观看| 丰满饥渴人妻一区二区三| 日本a在线网址| av国产精品久久久久影院| 久久久精品免费免费高清| 菩萨蛮人人尽说江南好唐韦庄| 日韩一区二区三区影片| 男女边摸边吃奶| 一本综合久久免费| 极品人妻少妇av视频| 亚洲五月色婷婷综合| 黄频高清免费视频| 亚洲欧美一区二区三区久久| 亚洲av日韩精品久久久久久密 | 精品国产一区二区三区四区第35| 99久久精品国产亚洲精品| 另类精品久久| 欧美日韩成人在线一区二区| 欧美精品啪啪一区二区三区 | 精品久久久久久久毛片微露脸 | 最近手机中文字幕大全| 岛国毛片在线播放| 久久久久久久精品精品| 亚洲av美国av| 制服人妻中文乱码| 亚洲视频免费观看视频| 777米奇影视久久| 亚洲伊人色综图| 久久精品亚洲av国产电影网| 狂野欧美激情性bbbbbb| 高清黄色对白视频在线免费看| 亚洲国产欧美网| 99国产精品免费福利视频| 丰满人妻熟妇乱又伦精品不卡| 91成人精品电影| 久久精品久久久久久噜噜老黄| 建设人人有责人人尽责人人享有的| 高清黄色对白视频在线免费看| av国产久精品久网站免费入址| 亚洲精品美女久久久久99蜜臀 | 岛国毛片在线播放| av在线播放精品| 宅男免费午夜| 国产激情久久老熟女| 国产野战对白在线观看| 黄色片一级片一级黄色片| 亚洲av综合色区一区| 麻豆av在线久日| 欧美精品啪啪一区二区三区 | 亚洲一区二区三区欧美精品| 老司机亚洲免费影院| 久久女婷五月综合色啪小说| 美女福利国产在线| 国产精品.久久久| 青青草视频在线视频观看| 国产精品国产av在线观看| 91麻豆av在线| 青春草视频在线免费观看| 国产成人91sexporn| 在线亚洲精品国产二区图片欧美| 美女主播在线视频| 纯流量卡能插随身wifi吗| 在线观看免费日韩欧美大片| 一级a爱视频在线免费观看| 国产成人精品久久久久久| 一级毛片我不卡| av有码第一页| 真人做人爱边吃奶动态| 成人亚洲精品一区在线观看| 丝袜美足系列| 天天添夜夜摸| 国产成人a∨麻豆精品| 亚洲国产日韩一区二区| 日日夜夜操网爽| 亚洲成av片中文字幕在线观看| 久久久久精品国产欧美久久久 | www.999成人在线观看| 青春草亚洲视频在线观看| 在线观看免费视频网站a站| 亚洲国产欧美在线一区| 日韩视频在线欧美| 又粗又硬又长又爽又黄的视频| 欧美av亚洲av综合av国产av| 国产成人欧美在线观看 | 国产激情久久老熟女| 日本欧美视频一区| 一边亲一边摸免费视频| 80岁老熟妇乱子伦牲交| 久久ye,这里只有精品| 赤兔流量卡办理| 夫妻午夜视频| 国产欧美日韩一区二区三 | 黄色毛片三级朝国网站| 成人免费观看视频高清| 亚洲国产欧美网| 色94色欧美一区二区| 看免费av毛片| 女性生殖器流出的白浆| 热re99久久国产66热| 国产精品.久久久| 啦啦啦视频在线资源免费观看| 久久精品国产亚洲av涩爱| 91老司机精品| 欧美激情极品国产一区二区三区| 国产成人a∨麻豆精品| 五月开心婷婷网| 国产1区2区3区精品| 国产三级黄色录像| 午夜日韩欧美国产| 香蕉丝袜av| 一本综合久久免费| 在线看a的网站| 不卡av一区二区三区| 少妇裸体淫交视频免费看高清 | 天堂8中文在线网| 亚洲av美国av| 国精品久久久久久国模美| 91麻豆精品激情在线观看国产 | 美女国产高潮福利片在线看| 天天躁夜夜躁狠狠久久av| 欧美精品亚洲一区二区| 亚洲精品久久久久久婷婷小说| 99热全是精品| 两性夫妻黄色片| 韩国高清视频一区二区三区| 亚洲一卡2卡3卡4卡5卡精品中文| 国产深夜福利视频在线观看| 桃花免费在线播放| 美女午夜性视频免费| 大片电影免费在线观看免费| 国产午夜精品一二区理论片| 国产成人av激情在线播放| www.999成人在线观看| 脱女人内裤的视频| 亚洲人成77777在线视频| 91九色精品人成在线观看| 叶爱在线成人免费视频播放| 久久精品亚洲熟妇少妇任你| 老鸭窝网址在线观看| 亚洲专区国产一区二区| 日本av免费视频播放| 欧美成人精品欧美一级黄| av福利片在线| 亚洲少妇的诱惑av| 久久久久国产一级毛片高清牌| 亚洲国产看品久久| 国产高清国产精品国产三级| 捣出白浆h1v1| 国产成人精品久久二区二区免费| 免费女性裸体啪啪无遮挡网站| 老司机亚洲免费影院| 男的添女的下面高潮视频| 亚洲精品一卡2卡三卡4卡5卡 | 亚洲av成人精品一二三区| 日本五十路高清| 久久久久视频综合| 少妇被粗大的猛进出69影院| 人人妻人人添人人爽欧美一区卜| 国产精品.久久久| 午夜91福利影院| 男女之事视频高清在线观看 | 国产精品av久久久久免费| 中文字幕精品免费在线观看视频| 男的添女的下面高潮视频| 亚洲欧美日韩另类电影网站| 国产一区二区三区av在线| 亚洲精品国产av蜜桃| 少妇人妻久久综合中文| 青春草视频在线免费观看| 久久久久精品国产欧美久久久 | 亚洲九九香蕉| 久久久久网色| 大香蕉久久成人网| 国产精品久久久av美女十八| www.精华液| 熟女av电影| 看免费av毛片| 欧美日韩亚洲综合一区二区三区_| 在现免费观看毛片| 欧美亚洲 丝袜 人妻 在线| 亚洲精品中文字幕在线视频| 日韩一卡2卡3卡4卡2021年| 午夜影院在线不卡| 国产高清国产精品国产三级| 国产精品九九99| 美国免费a级毛片| 国产福利在线免费观看视频| 亚洲欧美激情在线| 亚洲人成电影免费在线| 91成人精品电影| 美女扒开内裤让男人捅视频| 亚洲,欧美,日韩| 美女福利国产在线| 国产福利在线免费观看视频| 韩国高清视频一区二区三区| 国产精品秋霞免费鲁丝片| 国产成人系列免费观看| 搡老乐熟女国产| 老汉色av国产亚洲站长工具| av电影中文网址| 伊人亚洲综合成人网| 午夜两性在线视频| 波多野结衣一区麻豆| 久久精品国产a三级三级三级| 亚洲成国产人片在线观看| 人妻 亚洲 视频| 欧美久久黑人一区二区| 老司机影院毛片| 中文乱码字字幕精品一区二区三区| 久久精品成人免费网站| 亚洲成av片中文字幕在线观看| 人人澡人人妻人| 久久女婷五月综合色啪小说| 亚洲午夜精品一区,二区,三区| 在线观看免费日韩欧美大片| 国产在线一区二区三区精| 又粗又硬又长又爽又黄的视频| 18禁国产床啪视频网站| 国产成人精品久久二区二区免费| 国产成人精品无人区| 天堂俺去俺来也www色官网| 国产欧美日韩综合在线一区二区| 久久国产精品男人的天堂亚洲| 中文字幕最新亚洲高清| 90打野战视频偷拍视频| 黄色视频在线播放观看不卡| 99精品久久久久人妻精品| 午夜福利视频在线观看免费| 国产视频首页在线观看| 久久精品人人爽人人爽视色| 老司机影院成人| 国产色视频综合| 精品亚洲乱码少妇综合久久| 久久久精品94久久精品| av有码第一页| 亚洲国产精品成人久久小说| 国产熟女午夜一区二区三区| 国产在线观看jvid| 青青草视频在线视频观看| 色网站视频免费| 一二三四在线观看免费中文在| 美女主播在线视频| 美女午夜性视频免费| 久久人人爽av亚洲精品天堂| 99精国产麻豆久久婷婷| 秋霞在线观看毛片| 极品人妻少妇av视频| 一区二区日韩欧美中文字幕| 黄色毛片三级朝国网站| 国产成人精品久久二区二区91| 激情视频va一区二区三区| 性少妇av在线| 男人添女人高潮全过程视频| 我的亚洲天堂| 狠狠婷婷综合久久久久久88av| 一边摸一边做爽爽视频免费| 日本vs欧美在线观看视频| 考比视频在线观看| 亚洲成人手机| 老鸭窝网址在线观看| 久久久久久久大尺度免费视频| 国产精品人妻久久久影院| 熟女av电影| 咕卡用的链子| 国产黄色免费在线视频| 国产深夜福利视频在线观看| 久久久精品94久久精品| 亚洲av电影在线进入| 97人妻天天添夜夜摸| 免费观看a级毛片全部| 少妇的丰满在线观看| 午夜福利,免费看| av电影中文网址| 久久99精品国语久久久| 女人高潮潮喷娇喘18禁视频| 亚洲精品国产色婷婷电影| 欧美激情极品国产一区二区三区| 手机成人av网站| 中文字幕人妻熟女乱码| 99久久人妻综合| 一级毛片 在线播放| 国产免费又黄又爽又色| 夫妻午夜视频| 尾随美女入室| 免费在线观看日本一区| 久久天堂一区二区三区四区| 亚洲av欧美aⅴ国产| 久久久久国产精品人妻一区二区| 精品久久久久久久毛片微露脸 | 中文字幕另类日韩欧美亚洲嫩草| 久久青草综合色| 深夜精品福利| 国精品久久久久久国模美| 成人亚洲精品一区在线观看| 在现免费观看毛片| www.av在线官网国产| 欧美人与性动交α欧美软件| av网站免费在线观看视频| 999精品在线视频| 最新在线观看一区二区三区 | 国产成人av教育| 亚洲国产毛片av蜜桃av| 夫妻午夜视频| 久久精品国产亚洲av涩爱| 丰满少妇做爰视频| 国产亚洲精品久久久久5区| 国产精品一区二区精品视频观看| 亚洲熟女精品中文字幕| 亚洲国产精品国产精品| 精品久久蜜臀av无| 日韩 欧美 亚洲 中文字幕| 国产高清国产精品国产三级| 精品人妻熟女毛片av久久网站| 精品少妇久久久久久888优播| 美女高潮到喷水免费观看| 人妻人人澡人人爽人人| 欧美精品人与动牲交sv欧美| 亚洲精品第二区| 制服诱惑二区| 18禁观看日本| 日韩一卡2卡3卡4卡2021年| 国产xxxxx性猛交| 丰满人妻熟妇乱又伦精品不卡| 深夜精品福利| 91精品伊人久久大香线蕉| 一级毛片黄色毛片免费观看视频| 曰老女人黄片| 免费高清在线观看视频在线观看| 欧美成人精品欧美一级黄| 久久久精品国产亚洲av高清涩受| 国产精品久久久av美女十八| www.熟女人妻精品国产| 国产又色又爽无遮挡免| 免费在线观看日本一区| 99国产综合亚洲精品| 手机成人av网站| 国产黄频视频在线观看| 三上悠亚av全集在线观看| 爱豆传媒免费全集在线观看| 80岁老熟妇乱子伦牲交| 色视频在线一区二区三区| 国产男女超爽视频在线观看| 男人爽女人下面视频在线观看| 欧美亚洲 丝袜 人妻 在线| 欧美成人精品欧美一级黄| 日韩一卡2卡3卡4卡2021年| 熟女av电影| 久久久久网色| 极品人妻少妇av视频| 91精品三级在线观看| 国产黄色视频一区二区在线观看| 亚洲av男天堂| 亚洲国产日韩一区二区| 一级黄片播放器| 777米奇影视久久| 波野结衣二区三区在线| 欧美日韩黄片免| svipshipincom国产片| 狠狠婷婷综合久久久久久88av| 纯流量卡能插随身wifi吗| 在线观看国产h片| 久久精品国产亚洲av高清一级|