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

    低溶解氧下氨氧化過(guò)程神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)控制模型

    2017-02-22 07:23:27馮紅利劉秀紅黃斯婷楊玉兵周薛揚(yáng)北京工業(yè)大學(xué)北京市水質(zhì)科學(xué)與水環(huán)境恢復(fù)工程重點(diǎn)實(shí)驗(yàn)室北京000中國(guó)人民大學(xué)環(huán)境學(xué)院北京0087
    中國(guó)環(huán)境科學(xué) 2017年1期
    關(guān)鍵詞:控制參數(shù)硝化氨氮

    馮紅利,劉秀紅,楊 慶*,黃斯婷,崔 斌,周 桐,楊玉兵,周薛揚(yáng)(.北京工業(yè)大學(xué),北京市水質(zhì)科學(xué)與水環(huán)境恢復(fù)工程重點(diǎn)實(shí)驗(yàn)室,北京 000;.中國(guó)人民大學(xué)環(huán)境學(xué)院,北京 0087)

    低溶解氧下氨氧化過(guò)程神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)控制模型

    馮紅利1,劉秀紅2,楊 慶1*,黃斯婷1,崔 斌1,周 桐1,楊玉兵1,周薛揚(yáng)1(1.北京工業(yè)大學(xué),北京市水質(zhì)科學(xué)與水環(huán)境恢復(fù)工程重點(diǎn)實(shí)驗(yàn)室,北京 100022;2.中國(guó)人民大學(xué)環(huán)境學(xué)院,北京 100872)

    在低溶解氧(DO)狀態(tài)下,以城市生活污水為研究對(duì)象,將神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)的方法應(yīng)用到污水處理過(guò)程中,建立了基于神經(jīng)網(wǎng)絡(luò)的氨氧化過(guò)程預(yù)測(cè)控制模型,預(yù)測(cè)并控制污水處理氨氧化過(guò)程. 該模型分為兩部分,一是根據(jù)在線pH值變化預(yù)測(cè)氨氧化結(jié)束時(shí)間,其相關(guān)系數(shù)R值為0.9985;二是根據(jù)在線pH值實(shí)時(shí)預(yù)測(cè)氨氮濃度,R值為0.9083.試驗(yàn)結(jié)果表明該模型預(yù)測(cè)精度高、可控性好,具有較好的適應(yīng)性和穩(wěn)定性,對(duì)實(shí)現(xiàn)并穩(wěn)定短程硝化以及促進(jìn)主流工藝(厭氧氨氧化)有重要的指導(dǎo)和借鑒意義.

    低溶解氧;氨氧化過(guò)程;神經(jīng)網(wǎng)絡(luò);模型;pH

    污水生物脫氮一直是國(guó)內(nèi)外研究的熱點(diǎn).目前國(guó)內(nèi)外對(duì)污水生物脫氮主要集中在新工藝和機(jī)理研究.新工藝如短程硝化反硝化工藝

    [1-4,12-13,20]、同步硝化反硝化工藝、厭氧氨氧化工藝[21],而這些工藝都直接或間接地涉及到氨氧化過(guò)程[17-19],因此研究氨氧化過(guò)程對(duì)這些工藝的實(shí)際工程應(yīng)用具有重要的指導(dǎo)與借鑒意義.城市污水的水質(zhì)、水量與環(huán)境條件等均隨時(shí)間變化,同時(shí),城市污水氨氧化過(guò)程具有非線性、時(shí)變性和不確定性等特點(diǎn),因此,難于建立精準(zhǔn)的數(shù)學(xué)模型;即使建立了數(shù)學(xué)模型, 結(jié)構(gòu)也較為復(fù)雜,難于實(shí)際工程應(yīng)用.為有利于模型的實(shí)際應(yīng)用,不但應(yīng)使模型盡量簡(jiǎn)化,控制算法簡(jiǎn)單,易于實(shí)現(xiàn),而且在系統(tǒng)不確定性因素的影響下能夠保持良好的性能,滿足實(shí)時(shí)控制的需求.神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)控制是把神經(jīng)網(wǎng)絡(luò)與預(yù)測(cè)控制相結(jié)合,將預(yù)測(cè)控制的滾動(dòng)優(yōu)化思想與神經(jīng)網(wǎng)絡(luò)精確描述非線性和不確定性動(dòng)態(tài)過(guò)程的特性有機(jī)結(jié)合[5-7].近年來(lái),許多學(xué)者提出了多種基于神經(jīng)網(wǎng)絡(luò)的預(yù)測(cè)控制模型

    [8-11],部分學(xué)者將神經(jīng)網(wǎng)絡(luò)用于進(jìn)、出水指標(biāo)的預(yù)測(cè)[9],例如將LM-BP神經(jīng)網(wǎng)絡(luò)技術(shù)應(yīng)用于水質(zhì)預(yù)測(cè)[11].但關(guān)于生因此,本研究以城市生活污水為研究對(duì)象,在物脫氮過(guò)程的神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)模型未見(jiàn)報(bào)道.

    低溶解氧(0.4~0.6mg/L)條件下,考察了氨氧化過(guò)程建立神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)控制模型的可行性,并對(duì)不同條件下,模型預(yù)測(cè)的準(zhǔn)確性等進(jìn)行了研究.

    1 材料與方法

    1.1 污泥與污水

    試驗(yàn)所用種泥取自穩(wěn)定運(yùn)行2年的大型中試短程 SBR內(nèi)活性污泥[14],污水取自北京市某家屬小區(qū)實(shí)際生活污水.進(jìn)水水質(zhì)和控制參數(shù)如表1.

    表1 進(jìn)水水質(zhì)與控制參數(shù)Table 1 Influent quality & control parameters

    1.2 試驗(yàn)裝置與運(yùn)行模式

    圖1 SBR試驗(yàn)裝置示意Fig.1 Schematic diagram of SBR

    小試SBR,材質(zhì)為有機(jī)玻璃,有效容積19.5L,高 1100mm,內(nèi)徑 150mm,試驗(yàn)裝置如圖 1所示.SBR進(jìn)水由液位計(jì)控制,每天運(yùn)行一個(gè)周期,整個(gè)反應(yīng)過(guò)程中,由加熱棒和溫控系統(tǒng)控制反應(yīng)溫度在 25 ,℃好氧反應(yīng)階段由空氣壓縮機(jī)曝氣,控制系統(tǒng)內(nèi)DO在0.4~0.6mg/L;缺氧階段投加乙醇作為反硝化碳源,整個(gè)反應(yīng)過(guò)程一直進(jìn)行攪拌,以保持系統(tǒng)內(nèi)泥水混合狀態(tài),好氧缺氧反應(yīng)時(shí)間按照實(shí)時(shí)控制策略[15]運(yùn)行,反應(yīng)結(jié)束后沉淀30min,排水12L.

    1.3 模型建立過(guò)程

    模型建立整體過(guò)程分為 3個(gè)部分:1.數(shù)據(jù)采集:通過(guò)在線儀表將實(shí)時(shí)采集數(shù)據(jù);2.數(shù)據(jù)處理:將第一步采集到的數(shù)據(jù)進(jìn)行預(yù)處理,從中選擇足夠的樣本和合適的變量參數(shù);3.模型建立:將處理后的數(shù)據(jù)進(jìn)行建模,首先選擇合適的模型,進(jìn)行訓(xùn)練、校正和測(cè)試,達(dá)到精度要求后,進(jìn)行在線預(yù)測(cè)、過(guò)程監(jiān)測(cè)及傳感器監(jiān)測(cè).

    1.3.1 預(yù)測(cè)模型變量參數(shù)的選擇 污水生物處理過(guò)程控制參數(shù)一般分為直接參數(shù)和間接參數(shù)兩種,直接參數(shù)是指通過(guò)在線傳感器直接獲得的COD、NH4+-N、NO2--N、NO3

    --N等污染物的濃度值,但由于這些直接參數(shù)的傳感器價(jià)格昂貴、存在滯后性等原因,一直以來(lái)應(yīng)用得不是很多.間接參數(shù)是指一些與污染物濃度存在一定關(guān)系,可以間接反映系統(tǒng)內(nèi)反應(yīng)過(guò)程的控制參數(shù),如,DO、pH和ORP等.間接參數(shù)由于檢測(cè)方便準(zhǔn)確穩(wěn)定,且經(jīng)濟(jì)實(shí)用,應(yīng)用比較廣泛[16].在好氧硝化階段,DO和 pH值可以指示氨氧化反應(yīng)終點(diǎn);當(dāng)硝化結(jié)束后,停止曝氣進(jìn)入缺氧階段,ORP和pH值曲線上的變化點(diǎn)可準(zhǔn)確指示反硝化終點(diǎn)[10]. pH值在氨氧化結(jié)束和反硝化結(jié)束時(shí)都會(huì)出現(xiàn)明顯的變化點(diǎn),與DO、ORP等控制參數(shù)相比,選取pH值作為短程深度脫氮過(guò)程的控制參數(shù)既可節(jié)省數(shù)據(jù)存儲(chǔ)的空間,又可減少控制器運(yùn)算次數(shù),使控制策略得以簡(jiǎn)化.本試驗(yàn)在低溶解氧(0.4~0.6mg/L)下實(shí)現(xiàn)穩(wěn)定的短程,好氧階段 DO維持恒定,因此,選取pH值作為BP神經(jīng)網(wǎng)絡(luò)的預(yù)測(cè)參數(shù)指標(biāo),用于預(yù)測(cè)氨氧化過(guò)程中氨氮濃度的變化. 1.3.2 神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)模型的建立 為預(yù)測(cè)氨氧化結(jié)束時(shí)間和氨氮濃度,本研究建立的 3層 BP神經(jīng)網(wǎng)絡(luò)如圖2所示,其輸入層包括在線監(jiān)測(cè)的pH值和好氧反應(yīng)過(guò)程間隔 30~60min實(shí)際測(cè)定的氨氮濃度,輸出層為預(yù)測(cè)的pH值和氨氮濃度.

    圖2 神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)Fig.2 Structure of neural network

    神經(jīng)網(wǎng)絡(luò)訓(xùn)練函數(shù)如下:

    激活函數(shù):激活函數(shù)不管對(duì)于識(shí)別率或收斂速度都有顯著的影響.在逼近高次曲線時(shí),S形函數(shù)精度比線性函數(shù)要高得多,隱含層的激活函數(shù)采取S函數(shù):f(x) = 1/(1+exp(-x)) (1)

    訓(xùn)練函數(shù)的優(yōu)化:traingdx為自適應(yīng) lr動(dòng)量梯度下降法,與其它訓(xùn)練函數(shù)(trainm, trainbfg, trainrp和traingd等)相比,該訓(xùn)練函數(shù)具有精度較高,運(yùn)行速度快,且內(nèi)存占用小的優(yōu)點(diǎn).

    1.BP神經(jīng)網(wǎng)絡(luò)的設(shè)計(jì):典型的BP神經(jīng)網(wǎng)絡(luò)為3層結(jié)構(gòu)(輸入層、隱含層和輸出層).輸入層接受算法輸入,在本模型中對(duì)應(yīng)為X.輸入層節(jié)點(diǎn)的個(gè)數(shù)和X的維數(shù)對(duì)應(yīng),本研究中,輸入X為標(biāo)量,輸入層節(jié)點(diǎn)個(gè)數(shù)為 1.隱含層結(jié)點(diǎn)個(gè)數(shù)可以自行設(shè)定,在本實(shí)例中建議設(shè)定為4個(gè).輸出層結(jié)點(diǎn)個(gè)數(shù)和輸出Y的維數(shù)對(duì)應(yīng),本研究中,輸出y為標(biāo)量,輸出層節(jié)點(diǎn)個(gè)數(shù)為1.

    2.算法基本思路:

    模型結(jié)構(gòu)和訓(xùn)練函數(shù)確定后,給定組樣本:即穩(wěn)定運(yùn)行的217批pH值和對(duì)應(yīng)NH4+-N數(shù)據(jù),對(duì)模型進(jìn)行訓(xùn)練,具體過(guò)程如下:

    1) 神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)為1-4-1型,選0.1作為最大容許誤差,并給出隱層到輸出層的權(quán)矩陣,輸入層到隱層的權(quán)矩陣,將其初始化為某一小的隨機(jī)權(quán)矩陣.

    2) 依次輸入樣本,計(jì)算網(wǎng)絡(luò)的實(shí)際輸出.

    3) 計(jì)算誤差.

    4) 計(jì)算各層權(quán)值的變化量.

    5) 調(diào)節(jié)權(quán)值.

    6) 若誤差在容許范圍內(nèi),則轉(zhuǎn)2),否則轉(zhuǎn)7).

    7) 若達(dá)到訓(xùn)練時(shí)間,結(jié)束;否則,轉(zhuǎn) 2)進(jìn)入下一輪的計(jì)算,若直到規(guī)定的次數(shù)(2500次)時(shí)還沒(méi)有使輸出層達(dá)到期望的輸出,則強(qiáng)行退出,終止程序.

    3.基于在線pH值預(yù)測(cè)氨氧化結(jié)束時(shí)間:在低DO(0.4~0.6mg/L)狀態(tài)下,建立一個(gè)3層的BP神經(jīng)網(wǎng)絡(luò)模型(包含 4個(gè)隱含層),利用輸出后的誤差來(lái)估計(jì)輸出層的直接前導(dǎo)層的誤差,再用這個(gè)誤差估計(jì)更前一層的誤差,如此一層一層的反傳下去,就獲得了所有其他各層的誤差估計(jì).神經(jīng)網(wǎng)絡(luò)在外界輸入樣本的刺激下不斷改變網(wǎng)絡(luò)的連接權(quán)值,以使網(wǎng)絡(luò)的輸出不斷地接近期望的輸出,權(quán)值不斷調(diào)整的過(guò)程,也就是網(wǎng)絡(luò)的學(xué)習(xí)訓(xùn)練過(guò)程.此過(guò)程一直進(jìn)行到網(wǎng)絡(luò)輸出的誤差減少到可接受的程度或進(jìn)行到預(yù)先設(shè)定的學(xué)習(xí)次數(shù)為止.采用穩(wěn)定運(yùn)行的217批pH值數(shù)據(jù)訓(xùn)練神經(jīng)網(wǎng)絡(luò)模型,通過(guò)滾動(dòng)預(yù)測(cè)的方法預(yù)測(cè)氨氧化結(jié)束時(shí)間.具體如下:一個(gè)SBR好氧反應(yīng)周期中,前50個(gè)點(diǎn)只進(jìn)行數(shù)據(jù)采集,不進(jìn)行預(yù)測(cè),從第51個(gè)pH值開(kāi)始預(yù)測(cè),用第1~50個(gè)pH值預(yù)測(cè)第51個(gè)pH值,將預(yù)測(cè)出的第51個(gè)當(dāng)作已知pH值,用第2~52個(gè)pH值預(yù)測(cè)第52個(gè)pH值,將預(yù)測(cè)的第52個(gè)作為已知pH值,用第3~52個(gè)pH值預(yù)測(cè)第53個(gè)pH值,以此類推,滾動(dòng)預(yù)測(cè),用前50個(gè)實(shí)際在線監(jiān)測(cè)的pH值預(yù)測(cè)后20個(gè)pH值值.在第51個(gè)點(diǎn)時(shí)顯示了預(yù)測(cè)的51~70個(gè)pH值值,根據(jù)預(yù)測(cè)的pH值變化曲線的 ‘氨谷’點(diǎn),可預(yù)測(cè)氨氧化結(jié)束時(shí)間.

    基于在線 pH值實(shí)時(shí)預(yù)測(cè)氨氮濃度:在低DO(0.4~0.6mg/L)狀態(tài)下,建立一個(gè)3層的BP神經(jīng)網(wǎng)絡(luò)模型(包含 4個(gè)隱含層),采用穩(wěn)定運(yùn)行的217批pH值和對(duì)應(yīng)NH4+-N數(shù)據(jù)訓(xùn)練神經(jīng)網(wǎng)絡(luò)模型,其中pH值為在線連續(xù)監(jiān)測(cè),氨氮為實(shí)際測(cè)定不同時(shí)間點(diǎn)的氨氮濃度.根據(jù)在線監(jiān)測(cè) pH值實(shí)時(shí)預(yù)測(cè)NH4+-N的具體預(yù)測(cè)方法:根據(jù)pH值與對(duì)應(yīng)NH4+-N濃度構(gòu)建模型,根據(jù)pH值預(yù)測(cè)氨氮結(jié)束時(shí)間模型中預(yù)測(cè)的pH值和所建立的pH值和氨氮模型,來(lái)預(yù)測(cè)對(duì)應(yīng)的下一時(shí)刻N(yùn)H4+-N值.

    1.4 分析方法

    每周進(jìn)行3-5次周期性檢測(cè)分析.每個(gè)周期根據(jù)pH值、DO(WTW3420,DO/pH值便攜式測(cè)定儀)變化關(guān)鍵點(diǎn)取樣.COD采用COD快速測(cè)定儀(連華,5B-3(C))測(cè)定.NH4+-N、NO2--N、NO3

    --N采用流動(dòng)注射儀(QC8500Series2)進(jìn)行測(cè)定,污泥濃度(MLSS)用重量法測(cè)定.

    2 結(jié)果與討論

    2.1 模型訓(xùn)練結(jié)果

    圖3 基于在線pH值預(yù)測(cè)氨氧化結(jié)束時(shí)間模型訓(xùn)練結(jié)果Fig.3 Train result of ammonia oxidizedtermination prediction model based on on-line pH

    圖4 基于在線pH值實(shí)時(shí)預(yù)測(cè)氨氮濃度模型訓(xùn)練結(jié)果Fig.4 Train result of real time ammonia concentration prediction model based on on-line pH

    實(shí)現(xiàn)與穩(wěn)定短程硝化均需準(zhǔn)確控制氨氧化結(jié)束時(shí)間,防止亞硝酸鹽氧化過(guò)程.采用實(shí)時(shí)控制方式,利用pH值曲線上的變化點(diǎn)確定氨氧化結(jié)束時(shí)間的控制方式中,只有當(dāng)氨氧化徹底結(jié)束,且亞硝酸鹽氧化開(kāi)始后,才能出現(xiàn)變化點(diǎn),易導(dǎo)致過(guò)度曝氣.基于在線pH值預(yù)測(cè)氨氧化結(jié)束時(shí)間及氨氮濃度,能夠更有效防止氨氮過(guò)度氧化.

    2.1.1 基于在線pH值預(yù)測(cè)氨氧化結(jié)束時(shí)間 圖3給出了基于在線pH值預(yù)測(cè)氨氧化結(jié)束時(shí)間模型的訓(xùn)練結(jié)果.由圖可知整體相關(guān)系數(shù)(R)值為0.99853,說(shuō)明該模型的相關(guān)性很好,擬合度高.均方誤差(MSE)可評(píng)價(jià)數(shù)據(jù)的變化程度,是衡量“平均誤差”的一種較方便的方法,本模型MSE值僅為 0.00030108,說(shuō)明預(yù)測(cè)模型描述實(shí)驗(yàn)數(shù)據(jù)的精確度很高.

    2.1.2 基于在線pH值實(shí)時(shí)預(yù)測(cè)氨氮濃度 圖4給出了基于在線pH值實(shí)時(shí)預(yù)測(cè)氨氮濃度模型的訓(xùn)練結(jié)果,其整體R值為0.9083,說(shuō)明該模型的相關(guān)性較好,擬合度高.但由于pH值為實(shí)時(shí)在線監(jiān)測(cè)的數(shù)據(jù),數(shù)據(jù)采集時(shí)間間隔 1min,而氨氮為30~60min不同時(shí)間點(diǎn),關(guān)鍵反應(yīng)階段氨氮濃度數(shù)據(jù),氨氮樣本數(shù)據(jù)相對(duì)較少,這使得該模型精確度略有下降.雖然該模型精確度低于基于 pH值預(yù)測(cè)氨氧化結(jié)束時(shí)間的模型,但整體仍具有較好的預(yù)測(cè)效果.

    2.2 不同運(yùn)行條件,氨氧化結(jié)束時(shí)間預(yù)測(cè)效果

    為進(jìn)一步考察不同運(yùn)行條件下,氨氧化結(jié)束時(shí)間預(yù)測(cè)效果,本研究針對(duì)溫度和污泥濃度下氨氧化結(jié)束時(shí)間預(yù)測(cè)效果進(jìn)行了研究.

    2.2.1 不同溫度條件,氨氧化結(jié)束時(shí)間預(yù)測(cè)效果 圖 5給出了不同溫度下氨氧化結(jié)束時(shí)間預(yù)測(cè)效果,在污泥濃度(MLSS)為2800mg/L時(shí),溫度分別為21℃、26℃、31 ,℃氨氮去除率均在99%以上.由圖可知,隨著溫度的升高,氨氧化速率逐漸提高.模型預(yù)測(cè)誤差精度分別為2.78%、2.08%、3.61%以內(nèi),預(yù)測(cè)精度均較高,效果顯著.

    2.2.2 不同 MLSS下氨氧化結(jié)束時(shí)間預(yù)測(cè)效果 圖6給出了不同污泥濃度(MLSS)下氨氧化結(jié)束時(shí)間預(yù)測(cè)效果圖,在 26℃時(shí),MLSS分別為1400和2800mg/L時(shí),氨氮去除率均在99%以上.隨著污泥濃度的提高,氨氧化速率也隨之加快,預(yù)測(cè)誤差精度分別為0%和3.13%.該模型的預(yù)測(cè)精度高,預(yù)測(cè)效果好.

    圖5 不同溫度下氨氧化結(jié)束時(shí)間預(yù)測(cè)效果Fig.5 Predict result of ammonia oxidation termination at ddifferent temperature

    圖6 不同MLSS下氨氧化結(jié)束時(shí)間預(yù)測(cè)效果Fig.6 Predict result of ammonia oxidation termination at ddifferent MLSS

    2.2.3 長(zhǎng)期運(yùn)行 SBR,預(yù)測(cè)氨氧化結(jié)束時(shí)間效果 為進(jìn)一步確定模型預(yù)測(cè)準(zhǔn)確性,采用建立的模型對(duì)運(yùn)行160天的SBR氨氧化結(jié)束時(shí)間進(jìn)行預(yù)測(cè).圖7比較了實(shí)際測(cè)定值與模型預(yù)測(cè)結(jié)果.氨氧化結(jié)束時(shí)間預(yù)測(cè)值與實(shí)際值基本相同,預(yù)測(cè)精度在80%以上.

    圖7 實(shí)際測(cè)定值與模型預(yù)測(cè)效果比較Fig.7 Compare of real and predict results of ammonia oxidation termination

    基于上述研究結(jié)果,氨氮結(jié)束時(shí)間預(yù)測(cè)模型具有較好的適應(yīng)性,在不同溫度和污泥濃度條件下,預(yù)測(cè)精度仍較高;同時(shí),從長(zhǎng)期運(yùn)行 SBR預(yù)測(cè)結(jié)果看,盡管水質(zhì)及運(yùn)行控制參數(shù)有所變化,氨氧化結(jié)束預(yù)測(cè)時(shí)間有所不同,但模型預(yù)測(cè)精度仍較高,說(shuō)明模型具有良好的適應(yīng)性、準(zhǔn)確性和穩(wěn)定性.基于BP神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)氨氧化結(jié)束時(shí)間,為實(shí)現(xiàn)與穩(wěn)定短程硝化提供了一種新的方法.

    2.3 基于在線pH值實(shí)時(shí)預(yù)測(cè)氨氮濃度效果

    城市污水自養(yǎng)脫氮工藝已成為國(guó)內(nèi)外學(xué)者廣泛研究熱點(diǎn).兩段式短程厭氧氨氧化工藝,不但要實(shí)現(xiàn)短程,且要求僅有約一半的氨氮氧化為亞硝氮,即實(shí)現(xiàn)半短程硝化,這是是實(shí)現(xiàn)主流厭氧氨氧化的關(guān)鍵和難點(diǎn).圖 8A比較了模型預(yù)測(cè)氨氮濃度和實(shí)際測(cè)定氨氮濃度,結(jié)果顯示模型可準(zhǔn)確預(yù)測(cè)氨氮濃度,與實(shí)際氨氮濃度基本相同.基于實(shí)現(xiàn)半短程的條件,需將曝氣時(shí)間控制在H點(diǎn),控制NH4+-N/NO2--N為1左右.圖8B給出了采用模型預(yù)測(cè)半短程的效果,模型預(yù)測(cè)過(guò)程中同時(shí)取樣測(cè)定水質(zhì)情況,結(jié)果顯示模型能夠較準(zhǔn)確地預(yù)測(cè)半短程硝化曝氣時(shí)間,預(yù)測(cè)氨氮濃度與實(shí)測(cè)氨氮濃度基本相同.

    圖8 基于在線pH值實(shí)時(shí)預(yù)測(cè)氨氮效果Fig.8 Result of real-time predict ammonia based on on-line pH

    3 結(jié)論

    3.1 在低DO條件下,采用城市污水,在SBR反應(yīng)器中,選取 pH值作為預(yù)測(cè)參數(shù)指標(biāo)建立了三層 BP神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)模型.氨氧化過(guò)程結(jié)束時(shí)間預(yù)測(cè)控制模型動(dòng)態(tài)性能和自適應(yīng)性能較好,模型R值高達(dá)0.9985.在不同溫度和污泥濃度條件下,預(yù)測(cè)精度仍較高;同時(shí),從長(zhǎng)期運(yùn)行 SBR預(yù)測(cè)結(jié)果看,盡管水質(zhì)及運(yùn)行控制參數(shù)有所變化,氨氧化結(jié)束預(yù)測(cè)時(shí)間有所不同,但模型預(yù)測(cè)精度仍較高,說(shuō)明模型具有良好的適應(yīng)性、準(zhǔn)確性和穩(wěn)定性.基于 BP神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)氨氧化結(jié)束時(shí)間,為實(shí)現(xiàn)與穩(wěn)定短程硝化提供了一種新的方法.

    3.2 基于pH值預(yù)測(cè)反應(yīng)過(guò)程中氨氮濃度模型,雖然受氨氮量和進(jìn)水氨氮濃度波動(dòng)過(guò)大的影響,但該模型的R 值仍可達(dá)0.9083.模型可準(zhǔn)確預(yù)測(cè)氨氮濃度,與實(shí)際氨氮濃度基本相同.基于實(shí)現(xiàn)半短程的條件,模型能夠較準(zhǔn)確地預(yù)測(cè)半短程硝化曝氣時(shí)間,這為實(shí)現(xiàn)城市污水兩段式短程厭氧氨氧化工藝提供了一種新的控制方法和途徑.

    [1] Hellinga C, Schellen A A J C, Mulder J W, et al. The SHARON process: An innovative method for nitrogen removal fromammonium-rich wastewater [J]. Water Science & Technology, 1998,37(9):135-142.

    [2] Chung J W, Bae W. Nitrite reduction by a mixed culture under conditions relevant to shortcut biological nitrogen removal [J]. Biodegradation, 2002,13(3):163-170.

    [3] Peng Y Z, Zhu G B. Biological nitrogen removal with nitrification and denitrification via nitrite pathway [J]. Applied Microbiology &Biotechnology, 2006,73(1):15-26.

    [4] 李澤兵,李 軍,李 妍,等.短程硝化反硝化技術(shù)研究進(jìn)展 [J].給水排水, 2011,37(4):163-168.

    [5] Aslan S, Miller L, Dahab M. Ammonium oxidation via nitrite accumulation under limited oxygen concentration in sequencing batch reactors [J]. Bioresource Technology, 2009,100(2):59-664.

    [6] 許 超,陳治剛,邵惠鶴.預(yù)測(cè)控制技術(shù)及應(yīng)用發(fā)展綜述 [J]. 化工自動(dòng)化及儀表, 2002,29(3):1-10.

    [7] 陳增強(qiáng),袁著祉,張 燕.基于神經(jīng)網(wǎng)絡(luò)的非線性預(yù)測(cè)控制綜述[J]. 控制工程, 2002,9(4):7-11.

    [8] 段向軍.基于神經(jīng)網(wǎng)絡(luò)的預(yù)測(cè)控制方法研究 [D]. 黑龍江:大慶石油學(xué)院, 2005.

    [9] 韓紅貴,甄博然,喬俊飛.動(dòng)態(tài)結(jié)構(gòu)優(yōu)化神經(jīng)網(wǎng)絡(luò)及其在溶解氧控制中的應(yīng)用 [J]. 信息與控制, 2010,39(3):354-360.

    [10] 管 秋,王萬(wàn)良,徐新黎,等.基于神經(jīng)網(wǎng)絡(luò)的污水處理指標(biāo)軟測(cè)量研究 [J]. 環(huán)境污染與防治, 2006,28(2):156-158.

    [11] Wang W L, Ren M, Guan Q. Soft-sensing method for wastewater treatment based on BP neural network [C]//Proc. of the 4th World Congress on Intelligent Control and Automation. 上海:同濟(jì)大學(xué)出版社, 2002:2330-2332.

    [12] 田建平,曹東衛(wèi),李海楠.LM-BP神經(jīng)網(wǎng)絡(luò)在于橋水庫(kù)水質(zhì)預(yù)測(cè)中的應(yīng)用 [J]. 水利信息化, 2010,3:31-34.

    [13] Li J P, Elliott D. Nielsen M, et al. Long-term partial nitrification in an intermittently aerated sequencing batch reactor (SBR) treating ammonium-rich wastewater under controlled oxygenlimited conditions [J]. Biochemical Engineering Journal, 2011, 55(3):215-222.

    [14] Liu G Q, Wang J M. Long-term low DO enriches and shifts nitrifier community in activated sludge [J]. Environmental Science Technology, 2013,47(10):5109-5117.

    [15] 顧聲波,王淑瑩,彭永臻.短程深度脫氮中試工藝的低溫啟動(dòng)和維持 [J]. 環(huán)境科學(xué), 2013,34(8):3164-3170.

    [16] Yang Q, Peng Y Z, Liu X H, et al. Nitrogen removal via nitrite form municipal wastewater at low temperatures using real-time control to optimize nitrifying communities [J]. Environment Science Technology, 2007,41(23):8159-8164.

    [17] 楊 慶,彭永臻,王淑瑩,等.SBR法短程深度脫氮過(guò)程分析與控制模式的確立 [J]. 環(huán)境科學(xué), 2009,30(4):1084-1089.

    [18] Ruiz G, Jeison D, Chamy R. Nitrification with high nitrite accumulation for the treatment of wastewater with high ammonia concentration [J]. Water Research, 2003,37(1):1371-1377.

    [19] Xue Y, Yang F L, Liu S T, et al. The influence of control factors on the start up and operation for partial nitrification in membrane bioreactor [J]. Bioresource Technology, 2009,100(3):1055-1060.

    [20] Munoz C, Rojas D, Candia O, et al. Supervisory control system to enhance partial nitrification in an activated sludge reactor [J]. Chemical Engineering Journal, 2009,145(3):453-460.

    [21] 王淑瑩,黃惠珺,郭建華,等. DO對(duì)SBR短程硝化系統(tǒng)的短期和長(zhǎng)期影響 [J]. 北京工業(yè)大學(xué)學(xué)報(bào), 2010,36(8):1104-1110.

    [22] Du Y, Peng Y Z, Cao S B, et al. Advanced nitrogen removal with simultaneous Anammox and denitrification in sequencing batch reactor [J]. Bioresource Technology, 2014,162(6):316-322.

    Neural network prediction and control model for ammonia oxidizing process under low DO concentration.

    FENG Hong-li1, LIU Xiu-hong2, YANG Qing1*,HUANG Si-ting1, CUI Bin1, ZHOU Tong1, YANG Yu-bing1, ZHOU Xue-yang1
    (1.Key Laboratory of Beijing Water Quality Science and Water Environment Recovery Engineering, Beijing University of Technology, Beijing 100022;2.School of Environment & Natural Resources, Renmin University of China, Beijing 100872). China Environmental Science, 2017,37(1):139~145

    Under low dissolved oxygen (DO) concentration, the neural network prediction method was applied in SBR for treating domestic wastewater. The neural network control model was built to predict and control the ammonia oxidizing process. The model was divided into two parts. In the first part with the correlation coefficient (R value) of 0.9985, the end of ammonia oxidization was predicted based on the on-line pH variations. In the second part with R value of 0.9083, the ammonia concentration was real-time predicted based on the on-line pH variations. The results showed that the model with high prediction accuracy, good controllability, better adaptability and stability, can not only benefit for achieving and stabilizing short-cut, but also promote the application of anaerobic ammonium oxidation for treating domestic wastewater.

    low dissolved oxygen;ammonia oxidizing process;neural network;model;pH

    X703

    A

    1000-6923(2017)01-0139-07

    馮紅利(1990-),男,河北保定人,北京工業(yè)大學(xué)碩士研究生,研究方向?yàn)槲鬯锩摰捌溥^(guò)程控制.

    2016-03-31

    國(guó)家自然科學(xué)基金(51508561);北京市委組織部青年拔尖團(tuán)隊(duì);北京市優(yōu)秀人才培養(yǎng)資助計(jì)劃

    * 責(zé)任作者, 副教授, yangqing@bjut.edu.cn

    猜你喜歡
    控制參數(shù)硝化氨氮
    懸浮物對(duì)水質(zhì)氨氮測(cè)定的影響
    化工管理(2022年14期)2022-12-02 11:43:52
    高超聲速飛行器滑模控制參數(shù)整定方法設(shè)計(jì)*
    改進(jìn)型T-S模糊神經(jīng)網(wǎng)絡(luò)的出水氨氮預(yù)測(cè)
    云南化工(2021年8期)2021-12-21 06:37:36
    Birkhoff系統(tǒng)穩(wěn)定性的動(dòng)力學(xué)控制1)
    MBBR中進(jìn)水有機(jī)負(fù)荷對(duì)短程硝化反硝化的影響
    基于PI與準(zhǔn)PR調(diào)節(jié)的并網(wǎng)逆變器控制參數(shù)設(shè)計(jì)
    黑龍江電力(2017年1期)2017-05-17 04:25:08
    氧化絮凝技術(shù)處理高鹽高氨氮廢水的實(shí)驗(yàn)研究
    厭氧氨氧化與反硝化耦合脫氮除碳研究Ⅰ:
    間位芳綸生產(chǎn)廢水氨氮的強(qiáng)化處理及工程實(shí)踐
    海水反硝化和厭氧氨氧化速率同步測(cè)定的15N示蹤法及其應(yīng)用
    国产中年淑女户外野战色| 男人添女人高潮全过程视频| 成人无遮挡网站| 在线看a的网站| 直男gayav资源| 免费观看无遮挡的男女| 2021天堂中文幕一二区在线观| 亚洲在久久综合| 亚洲久久久久久中文字幕| 国产精品爽爽va在线观看网站| 80岁老熟妇乱子伦牲交| 国产高清国产精品国产三级 | 亚洲精品影视一区二区三区av| 在线观看av片永久免费下载| 性色avwww在线观看| 免费观看av网站的网址| 精品人妻一区二区三区麻豆| 日本爱情动作片www.在线观看| 国产午夜精品久久久久久一区二区三区| 免费看光身美女| 麻豆精品久久久久久蜜桃| 久久97久久精品| 久久精品国产自在天天线| 三级男女做爰猛烈吃奶摸视频| 99久久人妻综合| 亚洲四区av| 一级毛片电影观看| 日本色播在线视频| 黄色配什么色好看| 国产乱人偷精品视频| 简卡轻食公司| 深爱激情五月婷婷| 少妇被粗大猛烈的视频| 欧美三级亚洲精品| 久久久精品免费免费高清| 亚洲欧美日韩另类电影网站 | 成人欧美大片| 亚洲av免费在线观看| 麻豆久久精品国产亚洲av| 成人免费观看视频高清| 夜夜爽夜夜爽视频| 欧美一区二区亚洲| 精品国产三级普通话版| 三级国产精品欧美在线观看| 国产色婷婷99| 一级av片app| 国产亚洲av嫩草精品影院| 有码 亚洲区| 国产精品福利在线免费观看| 国产av码专区亚洲av| 久久久精品免费免费高清| 亚洲成人一二三区av| 欧美xxxx黑人xx丫x性爽| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | freevideosex欧美| 亚洲欧美日韩东京热| 国产久久久一区二区三区| 最近2019中文字幕mv第一页| 亚洲,一卡二卡三卡| 国产成人精品婷婷| 日本黄大片高清| tube8黄色片| 97精品久久久久久久久久精品| 99久久中文字幕三级久久日本| 伦精品一区二区三区| 久久久a久久爽久久v久久| 中文字幕制服av| 在现免费观看毛片| 大片电影免费在线观看免费| 国产精品一及| 全区人妻精品视频| 真实男女啪啪啪动态图| 久久精品国产a三级三级三级| 欧美 日韩 精品 国产| 狂野欧美激情性bbbbbb| 高清欧美精品videossex| 国模一区二区三区四区视频| 国产av国产精品国产| 亚洲av一区综合| 国产黄a三级三级三级人| 成人特级av手机在线观看| 草草在线视频免费看| 神马国产精品三级电影在线观看| av在线观看视频网站免费| 在线观看av片永久免费下载| 婷婷色麻豆天堂久久| 午夜福利高清视频| 国产毛片a区久久久久| 小蜜桃在线观看免费完整版高清| 男人添女人高潮全过程视频| 男女无遮挡免费网站观看| a级一级毛片免费在线观看| 免费av观看视频| 91精品一卡2卡3卡4卡| 麻豆成人av视频| 欧美激情久久久久久爽电影| .国产精品久久| 永久网站在线| 亚洲真实伦在线观看| 免费播放大片免费观看视频在线观看| 99久久精品热视频| 大码成人一级视频| 欧美性感艳星| 六月丁香七月| 免费观看的影片在线观看| 国产淫片久久久久久久久| 亚洲人与动物交配视频| 三级国产精品片| 午夜精品国产一区二区电影 | 麻豆乱淫一区二区| 蜜桃亚洲精品一区二区三区| 在线观看人妻少妇| 欧美潮喷喷水| 人妻一区二区av| 国产一区二区三区综合在线观看 | 18禁裸乳无遮挡动漫免费视频 | 男人和女人高潮做爰伦理| 国产人妻一区二区三区在| 一级黄片播放器| 成人鲁丝片一二三区免费| 亚洲成人一二三区av| 久久久久久久国产电影| 亚洲精品一区蜜桃| 涩涩av久久男人的天堂| 欧美一级a爱片免费观看看| 建设人人有责人人尽责人人享有的 | 日本av手机在线免费观看| 欧美少妇被猛烈插入视频| 九九久久精品国产亚洲av麻豆| 日韩制服骚丝袜av| 久久99热这里只有精品18| 成人特级av手机在线观看| 在线观看国产h片| 色视频在线一区二区三区| 校园人妻丝袜中文字幕| 18禁裸乳无遮挡免费网站照片| 久久久久精品性色| 嘟嘟电影网在线观看| 中文字幕免费在线视频6| 久久精品国产亚洲av天美| 精品久久久久久久人妻蜜臀av| www.av在线官网国产| 成年版毛片免费区| 亚洲成人av在线免费| 亚洲av日韩在线播放| 中文欧美无线码| 男女下面进入的视频免费午夜| 草草在线视频免费看| xxx大片免费视频| 国产亚洲最大av| 国产精品久久久久久av不卡| 麻豆成人午夜福利视频| 久久久久精品性色| 最近的中文字幕免费完整| 超碰97精品在线观看| 99热这里只有精品一区| 免费看日本二区| 内地一区二区视频在线| 久久久a久久爽久久v久久| 亚洲三级黄色毛片| 国产午夜精品一二区理论片| 中文字幕久久专区| 免费av观看视频| 免费大片黄手机在线观看| 国产老妇女一区| 男的添女的下面高潮视频| 免费高清在线观看视频在线观看| 国产 一区精品| 国产伦精品一区二区三区视频9| 亚洲av成人精品一二三区| 自拍欧美九色日韩亚洲蝌蚪91 | 一级黄片播放器| 最近2019中文字幕mv第一页| 亚洲欧美中文字幕日韩二区| 日本黄色片子视频| 国产精品麻豆人妻色哟哟久久| 人妻一区二区av| 国产欧美日韩精品一区二区| 久久久a久久爽久久v久久| 老女人水多毛片| 成人午夜精彩视频在线观看| 国产毛片在线视频| www.av在线官网国产| 国产淫语在线视频| 国产又色又爽无遮挡免| 大香蕉久久网| av.在线天堂| 欧美日韩在线观看h| 亚洲av福利一区| 一级毛片电影观看| 久久久久久久精品精品| 在线观看免费高清a一片| 亚洲四区av| 久久99热6这里只有精品| 最近最新中文字幕大全电影3| 亚洲人成网站在线观看播放| 特级一级黄色大片| 99久久九九国产精品国产免费| 精品人妻偷拍中文字幕| 亚洲伊人久久精品综合| 精品一区在线观看国产| 日韩免费高清中文字幕av| 建设人人有责人人尽责人人享有的 | 九草在线视频观看| 国产精品一区二区性色av| 精品国产三级普通话版| 欧美性猛交╳xxx乱大交人| 亚洲国产最新在线播放| 久久久久性生活片| 观看免费一级毛片| 91精品一卡2卡3卡4卡| 亚洲色图av天堂| 最近中文字幕2019免费版| 2021天堂中文幕一二区在线观| 各种免费的搞黄视频| 亚洲一级一片aⅴ在线观看| 男女边吃奶边做爰视频| 一级毛片电影观看| 国产免费一区二区三区四区乱码| 久久精品国产亚洲网站| 最近最新中文字幕大全电影3| 欧美日本视频| 国产av码专区亚洲av| 2021天堂中文幕一二区在线观| 亚洲成人中文字幕在线播放| 免费不卡的大黄色大毛片视频在线观看| 丰满少妇做爰视频| 国产亚洲av片在线观看秒播厂| 欧美潮喷喷水| 麻豆久久精品国产亚洲av| 777米奇影视久久| 午夜激情福利司机影院| 国产又色又爽无遮挡免| 欧美一区二区亚洲| 国产精品蜜桃在线观看| 亚洲成色77777| 一本色道久久久久久精品综合| 在现免费观看毛片| 成年女人在线观看亚洲视频 | 亚洲va在线va天堂va国产| 成人美女网站在线观看视频| 亚洲国产高清在线一区二区三| xxx大片免费视频| 精品人妻熟女av久视频| 亚洲精华国产精华液的使用体验| 日本与韩国留学比较| 国产综合精华液| 欧美 日韩 精品 国产| 伦理电影大哥的女人| 一本色道久久久久久精品综合| 人妻少妇偷人精品九色| 亚洲精品一二三| 亚洲综合色惰| 高清毛片免费看| 亚洲av中文av极速乱| 午夜福利视频1000在线观看| 精品酒店卫生间| 国产精品爽爽va在线观看网站| 久久久a久久爽久久v久久| 97人妻精品一区二区三区麻豆| 乱系列少妇在线播放| 国产精品一区二区在线观看99| 美女xxoo啪啪120秒动态图| 成人美女网站在线观看视频| 我的女老师完整版在线观看| 成年av动漫网址| 亚洲av国产av综合av卡| 亚洲国产色片| 你懂的网址亚洲精品在线观看| av在线播放精品| 成人午夜精彩视频在线观看| 一边亲一边摸免费视频| 国产免费又黄又爽又色| 爱豆传媒免费全集在线观看| av在线app专区| 黄片无遮挡物在线观看| 午夜激情福利司机影院| 天天躁夜夜躁狠狠久久av| 精品久久久久久久末码| 国产毛片在线视频| 成人特级av手机在线观看| tube8黄色片| 成人国产麻豆网| 啦啦啦啦在线视频资源| 神马国产精品三级电影在线观看| 亚洲欧美一区二区三区国产| 日本-黄色视频高清免费观看| 久久6这里有精品| 国产成人午夜福利电影在线观看| 亚州av有码| 精品99又大又爽又粗少妇毛片| 久久久午夜欧美精品| 男人和女人高潮做爰伦理| 久久女婷五月综合色啪小说 | 国产午夜精品一二区理论片| 国国产精品蜜臀av免费| 青春草视频在线免费观看| 国产精品秋霞免费鲁丝片| 国产69精品久久久久777片| 在线观看av片永久免费下载| 黄色视频在线播放观看不卡| 我的老师免费观看完整版| 精品久久久精品久久久| 欧美最新免费一区二区三区| 国产男女内射视频| 在线观看人妻少妇| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 午夜福利网站1000一区二区三区| 国内精品宾馆在线| 在线天堂最新版资源| 性插视频无遮挡在线免费观看| 国产成人freesex在线| 久久影院123| 亚洲av中文字字幕乱码综合| 中国三级夫妇交换| 在线免费观看不下载黄p国产| 国产毛片在线视频| 舔av片在线| 老女人水多毛片| 身体一侧抽搐| 日韩欧美精品免费久久| av天堂中文字幕网| 国产黄片视频在线免费观看| 色5月婷婷丁香| 日韩强制内射视频| 中文资源天堂在线| 在线观看国产h片| 赤兔流量卡办理| 亚洲av福利一区| 免费黄频网站在线观看国产| 亚洲av.av天堂| 直男gayav资源| 国内揄拍国产精品人妻在线| 亚洲精品乱码久久久v下载方式| 国产真实伦视频高清在线观看| 麻豆乱淫一区二区| 日本与韩国留学比较| 免费在线观看成人毛片| 亚洲精品久久久久久婷婷小说| 最近2019中文字幕mv第一页| 国产爱豆传媒在线观看| 少妇 在线观看| 天天躁夜夜躁狠狠久久av| 亚洲国产色片| 交换朋友夫妻互换小说| 久久久久精品性色| 亚洲精品自拍成人| 在线免费观看不下载黄p国产| 久久热精品热| 麻豆乱淫一区二区| 日韩欧美精品v在线| 国产乱来视频区| 日日啪夜夜撸| 99热这里只有是精品在线观看| 亚洲精品亚洲一区二区| 久久精品国产亚洲网站| 国产精品无大码| 热re99久久精品国产66热6| 亚洲精品视频女| 水蜜桃什么品种好| 久久精品人妻少妇| 99热全是精品| 国产精品福利在线免费观看| 日本黄色片子视频| av在线蜜桃| 国产男人的电影天堂91| 美女主播在线视频| 国产高清有码在线观看视频| 欧美高清性xxxxhd video| 成人高潮视频无遮挡免费网站| 大话2 男鬼变身卡| av一本久久久久| 麻豆乱淫一区二区| 九草在线视频观看| 97在线视频观看| 久久精品国产a三级三级三级| 国产永久视频网站| 久久午夜福利片| 麻豆成人av视频| 久久精品综合一区二区三区| 熟妇人妻不卡中文字幕| 国产免费一级a男人的天堂| 69av精品久久久久久| 精品视频人人做人人爽| www.色视频.com| 国产亚洲一区二区精品| 国产精品久久久久久久电影| 免费观看性生交大片5| 国产成人精品一,二区| 精品人妻视频免费看| 久久久久网色| 99久久九九国产精品国产免费| 久久久久久久午夜电影| 国产伦理片在线播放av一区| 一区二区三区四区激情视频| 亚洲成人av在线免费| 插阴视频在线观看视频| 午夜福利在线观看免费完整高清在| 五月天丁香电影| 久久久久久九九精品二区国产| 午夜爱爱视频在线播放| 亚洲精品,欧美精品| 国产亚洲一区二区精品| 在线观看一区二区三区| 老女人水多毛片| xxx大片免费视频| 秋霞伦理黄片| 国产综合精华液| 欧美激情久久久久久爽电影| 午夜老司机福利剧场| 九九在线视频观看精品| 简卡轻食公司| av播播在线观看一区| 国产精品蜜桃在线观看| 国产成人免费观看mmmm| av在线老鸭窝| 丰满乱子伦码专区| 天天一区二区日本电影三级| 性插视频无遮挡在线免费观看| 国产av码专区亚洲av| 日韩亚洲欧美综合| 亚洲一区二区三区欧美精品 | 亚洲欧美日韩东京热| 毛片一级片免费看久久久久| 高清午夜精品一区二区三区| 亚洲成人久久爱视频| 一区二区三区精品91| 晚上一个人看的免费电影| 视频区图区小说| 老女人水多毛片| 国产av不卡久久| 国产精品人妻久久久久久| 在线观看美女被高潮喷水网站| 国产成人精品福利久久| 99久久精品国产国产毛片| 交换朋友夫妻互换小说| 热re99久久精品国产66热6| 在线观看一区二区三区| 99热这里只有是精品50| 久久久久久国产a免费观看| kizo精华| 在线观看免费高清a一片| 在线天堂最新版资源| 丝袜脚勾引网站| 2021天堂中文幕一二区在线观| 中文字幕亚洲精品专区| 亚洲精品色激情综合| 波多野结衣巨乳人妻| 午夜免费鲁丝| 波野结衣二区三区在线| 下体分泌物呈黄色| 中文欧美无线码| 高清毛片免费看| 国产一区亚洲一区在线观看| 在线播放无遮挡| 丝袜美腿在线中文| 国产伦理片在线播放av一区| 中文天堂在线官网| 精品少妇久久久久久888优播| 别揉我奶头 嗯啊视频| 国产女主播在线喷水免费视频网站| 最近中文字幕高清免费大全6| 午夜精品一区二区三区免费看| 亚洲精品日韩av片在线观看| 亚洲国产色片| 伦理电影大哥的女人| 偷拍熟女少妇极品色| 成人美女网站在线观看视频| 国内揄拍国产精品人妻在线| 美女高潮的动态| 交换朋友夫妻互换小说| 美女内射精品一级片tv| 国产精品国产av在线观看| 在线观看三级黄色| 亚洲图色成人| 国产一区二区三区av在线| 久久精品夜色国产| 亚洲欧美精品自产自拍| 日韩强制内射视频| 亚洲人成网站在线播| 午夜福利在线在线| 在线观看免费高清a一片| 亚洲av成人精品一区久久| 一级片'在线观看视频| 亚洲性久久影院| 卡戴珊不雅视频在线播放| 亚洲欧美一区二区三区黑人 | 超碰av人人做人人爽久久| 一个人看的www免费观看视频| 久久久色成人| 亚洲熟女精品中文字幕| 国产日韩欧美亚洲二区| 久久久久久久午夜电影| 亚洲色图av天堂| 91精品伊人久久大香线蕉| 高清视频免费观看一区二区| 欧美少妇被猛烈插入视频| 亚洲av成人精品一区久久| 最近最新中文字幕免费大全7| 如何舔出高潮| 亚洲精品日韩av片在线观看| 在线免费观看不下载黄p国产| 精品人妻偷拍中文字幕| 男男h啪啪无遮挡| 人妻制服诱惑在线中文字幕| 少妇裸体淫交视频免费看高清| kizo精华| 国产国拍精品亚洲av在线观看| 中文字幕免费在线视频6| 亚洲天堂av无毛| 麻豆乱淫一区二区| 亚洲精品自拍成人| 男人爽女人下面视频在线观看| 99精国产麻豆久久婷婷| 国产高清有码在线观看视频| 成人欧美大片| 有码 亚洲区| 秋霞伦理黄片| 男女下面进入的视频免费午夜| 国产一级毛片在线| 久久国产乱子免费精品| 嘟嘟电影网在线观看| 亚洲成人精品中文字幕电影| 国产69精品久久久久777片| 国产精品无大码| 一本一本综合久久| 麻豆成人午夜福利视频| 国产毛片a区久久久久| 我要看日韩黄色一级片| 久久这里有精品视频免费| 亚州av有码| 如何舔出高潮| 成人黄色视频免费在线看| 久久精品国产a三级三级三级| 亚洲av.av天堂| 大又大粗又爽又黄少妇毛片口| 久久久精品欧美日韩精品| 看免费成人av毛片| 久久韩国三级中文字幕| 亚洲成色77777| 国产亚洲av嫩草精品影院| 日韩一本色道免费dvd| 久久精品久久久久久久性| 3wmmmm亚洲av在线观看| 精品人妻一区二区三区麻豆| 少妇人妻久久综合中文| 国产成人a区在线观看| 欧美少妇被猛烈插入视频| 国产伦在线观看视频一区| 天天躁日日操中文字幕| 亚洲第一区二区三区不卡| 亚洲人与动物交配视频| 亚洲欧美日韩无卡精品| 精品少妇久久久久久888优播| 一级毛片aaaaaa免费看小| 精品午夜福利在线看| 日韩电影二区| 色哟哟·www| 看黄色毛片网站| 黄色怎么调成土黄色| 欧美+日韩+精品| 观看美女的网站| 特级一级黄色大片| 久久午夜福利片| 两个人的视频大全免费| 国产精品国产三级国产专区5o| 欧美日韩在线观看h| 国产精品久久久久久久久免| 亚洲精品亚洲一区二区| 又爽又黄无遮挡网站| 国产女主播在线喷水免费视频网站| 麻豆精品久久久久久蜜桃| 美女主播在线视频| 久久久久性生活片| 日韩欧美精品v在线| 亚洲精品国产av蜜桃| 亚洲欧洲国产日韩| 亚洲国产成人一精品久久久| 日韩精品有码人妻一区| 日本三级黄在线观看| 亚洲av免费在线观看| 欧美高清成人免费视频www| 国产免费视频播放在线视频| 大话2 男鬼变身卡| 色视频www国产| 菩萨蛮人人尽说江南好唐韦庄| a级一级毛片免费在线观看| 久久精品久久精品一区二区三区| 夜夜看夜夜爽夜夜摸| a级毛色黄片| 日韩不卡一区二区三区视频在线| 青春草视频在线免费观看| 综合色丁香网| 精品国产一区二区三区久久久樱花 | 日日摸夜夜添夜夜添av毛片| 国产成人免费无遮挡视频| 午夜福利网站1000一区二区三区| 欧美丝袜亚洲另类| 国产真实伦视频高清在线观看| 国产日韩欧美亚洲二区| 久久久久久久午夜电影| 日韩国内少妇激情av| 97热精品久久久久久| 两个人的视频大全免费| 嘟嘟电影网在线观看| 99热网站在线观看| 国产精品国产av在线观看| 国产免费视频播放在线视频| 麻豆久久精品国产亚洲av| 久久精品久久久久久久性| 免费电影在线观看免费观看| 熟女av电影| 一本久久精品| 久久99精品国语久久久| 别揉我奶头 嗯啊视频|