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

    基于線性自抗擾的城市污水脫氮控制仿真

    2016-05-11 02:15:03魏偉王小藝王藩劉載文北京工商大學(xué)計算機與信息工程學(xué)院北京100048
    化工學(xué)報 2016年3期
    關(guān)鍵詞:動力學(xué)模型控制污水

    魏偉,王小藝,王藩,劉載文(北京工商大學(xué)計算機與信息工程學(xué)院,北京 100048)

    ?

    基于線性自抗擾的城市污水脫氮控制仿真

    魏偉,王小藝,王藩,劉載文
    (北京工商大學(xué)計算機與信息工程學(xué)院,北京 100048)

    摘要:根據(jù)A/O工藝的脫氮原理,基于城市污水處理的1號基準動力學(xué)模型,分析影響脫氮過程的兩個重要因素:硝態(tài)氮和溶解氧,利用線性自抗擾實現(xiàn)城市污水處理過程的脫氮控制?;谙鯌B(tài)氮和溶解氧分別設(shè)計線性自抗擾控制,既避免對復(fù)雜數(shù)學(xué)模型的依賴,又克服影響污水出水水質(zhì)的波動因素,可獲得較好的污水處理效果。動態(tài)仿真結(jié)果表明,線性自抗擾具有更好的控制性能。基于線性自抗擾的脫氮控制有利于提高污水處理的脫氮效率,改善出水水質(zhì)。

    關(guān)鍵詞:污水;動力學(xué)模型;脫氮;線性自抗擾;控制;動態(tài)仿真

    2015-12-22收到初稿,2015-12-29收到修改稿。

    聯(lián)系人及第一作者:魏偉(1984—),男,博士,副教授。

    引 言

    在我國廢水中,化學(xué)需氧量和氨氮排放量為廢水中主要污染物,城市污水化學(xué)需氧量排放量為889.8 t,氨氮排放量為131.4萬噸[1]?;瘜W(xué)需氧量中含有大量有機氮,使得有機氮和氨氮占據(jù)城市污水的較大比例。在城市污水中,氮元素主要以4種形式存在:有機氮、氨氮、硝酸鹽和亞硝酸鹽。氮是維持生態(tài)系統(tǒng)營養(yǎng)物質(zhì)循環(huán)的重要元素,但是污水中氮含量過高又會使水體富營養(yǎng)化[2],進一步污染水體。因此,控制城市污水中的氨氮、有機氮和總氮濃度成為污水處理的重要指標。目前,廣泛使用的方法是生物脫氮,主要有3個基本過程:氨解、硝化和反硝化。氨解將有機氮轉(zhuǎn)化為氨氮;硝化將氨氮轉(zhuǎn)化為亞硝酸鹽和硝酸鹽;反硝化將亞硝酸鹽和硝酸鹽轉(zhuǎn)化為氮氣,進而實現(xiàn)生物脫氮[3-5]。其中,硝化反應(yīng)需氧量較大,反硝化反應(yīng)在缺氧環(huán)境下發(fā)生,需要合理控制污泥年齡及外加碳源[6]。

    在污水處理脫氮工藝中,應(yīng)用較多的是基于活性污泥法的A/O工藝[7]。A/O工藝由缺氧池、好氧池及沉淀池串聯(lián)而成。缺氧池的溶解氧濃度(dissolved oxygen,DO)一般小于0.2 mg·L?1,好氧池的DO濃度一般為2~4 mg·L?1。好氧池末端與缺氧池相連,形成內(nèi)循環(huán);沉淀池末端、入水端口與缺氧池相連,形成外循環(huán)。內(nèi)、外循環(huán)實現(xiàn)對污泥的循環(huán)利用以及污水循環(huán)處理。好氧池的DO濃度、內(nèi)循環(huán)流量、外循環(huán)流量以及外加碳源量都會影響A/O工藝的脫氮效率[8]。為研究基于A/O工藝的脫氮控制,選取國際水協(xié)會發(fā)布的基于污水處理廠工藝過程的1號基準仿真模型(Benchmark Simulation Model No.1,BSM1)。BSM1為研究污水處理控制提供了一個標準平臺,便于不同控制方案以及策略相互比較。本文研究基于BSM1的脫氮控制,利用氧轉(zhuǎn)移系數(shù)和內(nèi)循環(huán)流量控制DO濃度及硝態(tài)氮濃度。

    目前,用于污水處理的控制算法主要有PI控制、神經(jīng)網(wǎng)絡(luò)控制等。PI控制因結(jié)構(gòu)簡單、容易實現(xiàn)被最早用于污水處理領(lǐng)域[9-10]。然而,污水出水水質(zhì)受到污水來源(組分)的不確定性以及污水處理過程的非線性(主要是污水流量的波動、水質(zhì)組分的變化、組分間的相互作用以及天氣)等干擾因素的影響,PI控制因?qū)Ψ蔷€性擾動的抑制能力較差而使控制精度降低。神經(jīng)網(wǎng)絡(luò)控制通過大量輸入輸出數(shù)據(jù)訓(xùn)練來保證控制精度。文獻[11]介紹了一種通過輸入輸出數(shù)據(jù)辨識控制器參數(shù)的方法,稱為虛擬參考反饋調(diào)節(jié)(virtual reference feedback tuning,VRFT),并將其應(yīng)用于BSM1的控制中。VRFT通過輸入和輸出數(shù)據(jù)的反饋進行辨識,對控制器進行在線或者離線調(diào)節(jié),適應(yīng)于雙變量或者多變量非線性系統(tǒng)控制,但是VRFT與神經(jīng)網(wǎng)絡(luò)控制方法相似,控制精度取決于訓(xùn)練數(shù)據(jù)量的大小以及訓(xùn)練數(shù)據(jù)和實際數(shù)據(jù)的相關(guān)性。在城市污水處理過程中,影響出水水質(zhì)的不確定因素較多,VRFT的性能受到一定的限制。

    自抗擾控制是一種不依賴于對象和擾動精確模型的控制算法,它通過實時估計和補償影響控制性能的各種擾動而獲得良好的控制效果[12-14]。然而,自抗擾控制有十幾個可調(diào)參數(shù),不便整定,線性自抗擾控制(linear active disturbance rejection control,LADRC)將可調(diào)參數(shù)減少為控制器帶寬和觀測器帶寬,具有更為明確的物理意義且控制性能同樣優(yōu)異[15-19]。本文亦采用LADRC進行城市污水脫氮控制的仿真研究。

    1 BSM1模型脫氮分析

    BSM1由生化反應(yīng)池和二沉池組成(如圖1所示),其中5個生化反應(yīng)池由2個缺氧池和3個好氧池組成,使用活性污泥1號模型模擬每個生化反應(yīng)池中的生物進程[20]。二沉池共10層,使用二次指數(shù)沉淀速率模型模擬二沉池的沉淀過程[21-22]。BSM1通過內(nèi)循環(huán)和外循環(huán)回流使系統(tǒng)形成閉環(huán)。BSM1包含城市污水的13種基本組分,主要有溶解氧、含氮物質(zhì)、活性細菌等。

    圖1 BSM1工藝結(jié)構(gòu)Fig.1 Process structure of BSM1

    在MATLAB/Simulink平臺上,BSM1通過M語言和C語言實現(xiàn),每個模塊使用S函數(shù)封裝。除了生化反應(yīng)池、二沉池和進出污水管線外,需要在1號生化反應(yīng)池進水口之前加入延時模塊,避免在仿真過程中出現(xiàn)數(shù)據(jù)環(huán)。

    測試模型前,模型的流量、曝氣量等需依據(jù)實際情況進行假設(shè)和限制。模型測試分為兩個過程:靜態(tài)測試和動態(tài)測試。與文獻[9]的標準值比對,測試結(jié)果在合理的誤差范圍內(nèi),BSM1符合污水處理廠的標準模型要求。選取14 d污水數(shù)據(jù)(數(shù)據(jù)來源于http://www.benchmarkWWTP.org/)作為BSM1的輸入,圖2(a)為污水流量,圖2(b)為含氮水質(zhì)組分的濃度。

    仿真時間設(shè)為14 d,圖3(a)為1、2號生化反應(yīng)池的DO濃度,圖3(b)為3~5號生化反應(yīng)池的DO濃度。對比A/O工藝中對DO濃度范圍的限定,需控制5號生化反應(yīng)池的DO濃度,使BSM1缺氧池和好氧池的DO濃度符合A/O工藝要求。在缺氧池中,異養(yǎng)菌通過反硝化作用,將硝態(tài)氮轉(zhuǎn)化為氮氣;同時,含氮有機污染物通過氨解,將氨基轉(zhuǎn)換為游離態(tài)的氨氮進入好氧池中,通過自養(yǎng)菌的硝化作用,將氨氮轉(zhuǎn)化成硝態(tài)氮。BSM1有內(nèi)、外回流循環(huán)系統(tǒng),將硝化反應(yīng)生成的硝態(tài)氮送回缺氧池;同時培養(yǎng)出利于反硝化反應(yīng)的污泥。通過硝化反應(yīng)和反硝化反應(yīng),可實現(xiàn)脫氮。為提高脫氮效率,需合理控制好氧池中的DO濃度和缺氧池中的硝態(tài)氮濃度。好氧池中,若DO濃度過低,會抑制硝化反應(yīng);相反,過高的DO濃度又會增加曝氣設(shè)備的能耗,并且過量的DO經(jīng)內(nèi)、外循環(huán)進入缺氧池會抑制反硝化反應(yīng)。為降低能耗,需控制好氧池末端5號生化反應(yīng)池的DO濃度。硝態(tài)氮濃度是連接硝化和反硝化的中間產(chǎn)物,在好氧池中,硝態(tài)氮濃度不斷累積,控制缺氧池末端硝態(tài)氮濃度可提高反硝化反應(yīng)效率。由于硝態(tài)氮在好氧池中難以轉(zhuǎn)化和去除,故控制缺氧池末端硝態(tài)氮濃度對降低出水總氮含量有重要意義[11,23-24]。采取的控制方案見表1。

    圖2 污水流量及含氮水質(zhì)組分濃度Fig.2 Flow of waste water and concentration of nitrogen in waste water

    圖3 5個生化反應(yīng)池的溶解氧濃度Fig.3 Concentration of dissolved oxygen in 5 biochemical reaction pools

    表1 基于BSM1的污水處理控制方案Table 1 Programmes for process control of waste water based on BSM1

    2 線性自抗擾控制設(shè)計

    2.1 線性自抗擾控制

    三階線性自抗擾控制結(jié)構(gòu)如圖4所示。

    圖4 三階線性自抗擾控制結(jié)構(gòu)Fig.4 Structure of three-order LADRC

    圖中虛線框部分為線性自抗擾控制器,主要包括控制律部分和線性擴張狀態(tài)觀測器部分。除給定信號yr外,LADRC的輸入僅為控制量u和被控對象輸出y。線性擴張狀態(tài)觀測器依據(jù)u和y獲得狀態(tài)估計z1和z2以及總擾動的估計z3;控制量u由給定信號yr、狀態(tài)估計z1和z2、總擾動的估計z3以及參數(shù)b0組成。

    因此,線性自抗擾控制器在設(shè)計時并未用到精確的模型信息,這對提高閉環(huán)系統(tǒng)的魯棒性具有極其重要的作用。

    2.2 硝態(tài)氮的線性自抗擾控制設(shè)計

    生化反應(yīng)池中水質(zhì)組分的物料平衡方程為

    其中,k表示生化反應(yīng)池序號,當k= a,表示為內(nèi)循環(huán),當k =0,表示為流入污水,當k= r,表示為外循環(huán);Z表示水質(zhì)組分濃度,V表示生化反應(yīng)池體積,Q表示流出流量,i表示水質(zhì)組分序號,r為水質(zhì)組分反應(yīng)速率。2號生化反應(yīng)池中SNO,2的物料平衡方程為[9]

    1號生化反應(yīng)池中SNO,1的物料平衡方程為

    對式(2)求導(dǎo),再將式(5)代入,得到

    令SNO,2= x1,Qa= u(t),式(7)可改寫為

    其中,狀態(tài)x3表示未知擾動f( t ),χ(t )為未知擾動的變化率,式(8)可寫為

    其中

    式(9)的狀態(tài)擴張觀測器為

    于是,基于參數(shù)ωo的擴張狀態(tài)觀測器為

    PD控制器設(shè)計為

    其中,r1為SNO,2的期望值;取kd= 2ωc,kp=;取b為其附近的某一值b0。

    2.3 溶解氧的線性自抗擾控制設(shè)計

    5號生化反應(yīng)池DO濃度物料平衡方程[9]為

    其中,KLa為曝氣設(shè)備爆氣后的氧轉(zhuǎn)移系數(shù),S?為DO的未飽和系數(shù)。5號生化反應(yīng)池DO模型可視為一個一階系統(tǒng)

    圖5 受控的DO5和SNO,2響應(yīng)曲線Fig.5 Response curves of controlled DO5and SNO,2

    PD控制器為

    其中r2為DO的期望值;取kd=ωc;取b為其附近的某一值b0。

    3 仿真研究

    對溶解氧的LADRC參數(shù)設(shè)置為:b0= 1,ωo= 600,ωc= 400。通過仿真,溶解氧7~14 d跟蹤曲線如圖5(a)所示,DO5很好地跟蹤了期望值2 mg·L?1。對硝態(tài)氮的LADRC參數(shù)設(shè)置為:b0= 0.04,ωo= 600,ωc= 400,通過仿真,硝態(tài)氮7~14 d跟蹤曲線如圖5(b)所示,SNO,2較好地跟蹤了期望值1 mg·L?1。

    仿真結(jié)果表明,污水流量及污水組分發(fā)生變化時,LADRC能較好地克服這些干擾因素,保證控制性能。

    表2對比了LADRC、VRFT與PI控制的性能指標以及水質(zhì)能耗參數(shù)。LADRC控制下溶解氧和硝態(tài)氮的平均誤差、最大誤差和標準差都遠小于VRFT和PI控制,控制精度較高;采用LADRC的出水水質(zhì)也優(yōu)于其他兩種控制方法。表2中各指標定義為[9]:

    誤差絕對值積分

    最大誤差

    表2 LADRC、VRFT和PI控制性能對比Table 2 Comparisons among LADRC, VRFT and PI

    標準差

    出水水質(zhì)

    成本指數(shù):

    其中,AE為曝氣池電力消耗,PE為水泵能耗,SP為污泥能耗,ME為混合耗能,其計算式分別為

    表3 污水處理前后的水質(zhì)對比Table 3 Comparisons of effluent quality before and after control

    表4 不同天氣條件下污水處理后的水質(zhì)對比Table 4 Comparisons of effluent quality under different weather

    表3給出了污水處理前、后的水質(zhì)變化對比。①是污水未處理時的水質(zhì)指標;②為沒有施加控制時,污水處理后的出水水質(zhì)指標;③為經(jīng)LADRC脫氮控制后的污水出水水質(zhì)指標。數(shù)據(jù)②③基本達到了我國城鎮(zhèn)污水處理廠污染物排放一級A類標準[25]。數(shù)據(jù)③與數(shù)據(jù)②比較,經(jīng)脫氮控制后,氨氮濃度進一步下降58%,凱氏定氮下降42%,總氮含量上升6%,生化需氧量(biochemical oxygen demand,BOD)與化學(xué)需氧量(chemical oxygen demand,COD)基本保持不變,綜合水質(zhì)提高11%。

    在BSM1中,總氮為凱氏定氮、硝態(tài)氮濃度之和,凱氏定氮為氨氮和有機氮之和。從表3氨氮、凱氏定氮和總氮濃度可以得出,經(jīng)脫氮控制后,污水處理過程中氨解和硝化反應(yīng)得到促進,但是由于A/O工藝的缺氧池前置,使得硝化反應(yīng)在污水處理生化反應(yīng)的末端,部分硝態(tài)氮隨污水直接排出,導(dǎo)致在綜合水質(zhì)提高的情況下,總氮含量沒有降低。

    總之,污水經(jīng)脫氮處理后,在保持總氮含量變化不大的前提下,有機氮和氨氮濃度進一步降低,脫氮效率明顯提高,水質(zhì)得到改善。

    以上利用晴天的污水數(shù)據(jù)進行仿真,為驗證LADRC的抗擾能力,選取雨天和暴雨天的污水數(shù)據(jù)進行仿真。保持LADRC的控制參數(shù)不變,受控的DO5和SNO,2響應(yīng)曲線如圖6所示。

    顯然,即使天氣不同,污水流量及組分發(fā)生變化,圖6(a)所示受控溶解氧濃度幾乎保持不變,圖6(b)所示受控硝態(tài)氮濃度在第11~13天有變化外,其他時間亦幾乎保持不變。表4給出了幾種天氣情況下,經(jīng)LADRC脫氮控制后的污水出水水質(zhì)指標對比。

    表中數(shù)據(jù)基本達到了我國城鎮(zhèn)污水處理廠污染物排放一級A類標準[25]。可見,線性自抗擾控制能夠保證系統(tǒng)對各種不確定因素具有良好的魯棒性,以保證出水水質(zhì)。

    4 結(jié) 論

    A/O為城市污水脫氮處理較為常見的工藝。在脫氮控制過程中,由于污水自身的復(fù)雜性和波動性,其建模和控制難度很大。本文基于BSM1模型,研究基于A/O工藝的脫氮控制。仿真結(jié)果表明,LADRC控制性能較好,結(jié)構(gòu)簡單,易于在污水處理等復(fù)雜被控對象上實現(xiàn),是一種較為實用的污水處理控制方法,可有效改善污水的出水水質(zhì)。

    References

    [1] 國家環(huán)境保護總局. 2013年中國環(huán)境狀況公報污染物排放情況[R]. 北京, 2014.

    Ministry of Environmental Protection of the People’s Republic of China. The bulletin of discharge of pollutants in China in 2013 [R]. Beijing, 2014.

    [2] 鄭興燦, 李亞新. 污水除磷脫氮技術(shù) [M]. 北京: 中國建筑工業(yè)出版社, 1998: 221.

    ZHENG X C, LI Y X. The Technology of Wastewater Nitrogen and Phosphorus Removal [M]. Beijing: China Building Industry Press, 1998: 221.

    [3] 張景紅, 胡海濤. 城市污水氮源與脫氮技術(shù)的發(fā)展 [J]. 化學(xué)工程與裝備, 2010, 10 (2): 172-173. DOI: 10.3969/j.issn.1003-0735.2010. 02.059.

    ZHANG J H, HU H T. The development of technology of the nitrogen source and nitrogen removal [J]. Chemical Engineering & Equipment, 2010, 10 (2): 172-173. DOI: 10.3969/j.issn.1003-0735.2010.02.059.

    [4] 周雪, 劉娜. 淺談新型除磷脫氮技術(shù) [J]. 科學(xué)咨詢, 2014, 8 (31): 50-51. DOI: 10.3969/j.issn.1671-4822.2014.24.031.

    ZHOU X, LIU N. Discussion on the new nitrogen and phosphorus removal technology [J]. The Scientific Advisory, 2014, 8 (31): 50-51. DOI: 10.3969/j.issn.1671-4822.2014.24.031.

    [5] KASHIF R, DAE H A, DAE S L. Simultaneous organic carbon and nitrogen removal in an anoxic-oxic activated sludge system under various operating conditions [J]. Bioresource Technology, 2014, 162 (6): 373-378.

    [6] 劉陽, 白瑩, 周娜. 城市污水脫氮除磷中的化學(xué)問題 [J]. 環(huán)境保護科學(xué), 2008, 34 (3): 60-62. DOI: 10.3969/j.issn.1004-6216.2008.03.019.

    LIU Y, BAI Y, ZHOU N. Chemical problems during nitrogen and phosphate removal of municipal wastewater treatment [J]. Environmental Protection Science, 2008, 34 (3): 60-62. DOI: 10.3969/j.issn.1004-6216.2008.03.019.

    [7] 馬勇, 王淑瑩, 王曉蓮, 等. A/O脫氮工藝的控制策略和應(yīng)用性研究 [J]. 環(huán)境污染治理技術(shù)與設(shè)備, 2003, 4 (10): 18-22. DOI: 10.3969/j.issn.1673-9108.2003.10.004.

    MA Y, WANG S Y, WANG X L, et al. Control strategies and application research of A/O nitrogen removal process [J]. Techniques and Equipment for Environmental Pollution Control, 2003, 4 (10): 18-22. DOI: 10.3969/j.issn. 1673-9108.2003.10.004.

    [8] 彭永臻, 王曉蓮, 王淑瑩. A/O脫氮工藝影響因素及其控制策略的研究 [J]. 哈爾濱工業(yè)大學(xué)學(xué)報, 2005, 37 (8): 1053-1057. DOI: 10.3321/j.issn:0367-6234.2005.08.011.

    PENG Y Z, WANG X L, WANG S Y. Study on influence factors and control strategies A/O nitrogen removal process [J]. Journal of Harbin Institute of Technology, 2005, 37 (8): 1053-1057. DOI: 10.3321/j.issn:0367-6234.2005.08.011.

    [9] ALEX J, BENEDETTI L, COPP J, et al. Benchmark Simulation Model No.1 (BSM1) [R]. Sweden: Dept. of Industrial Electrical Engineering and Automation, Lund University, 2008.

    [10] VRECKO D, HVALA N, KOCIJAN J. Wastewater treatment benchmark: what can be achieved with simple control [C]// Proc. ICA 2001: vol. 1. Malm?, Sweden, 2001: 221-228.

    [11] DAVID J R, XAVIER F A, JEPPSSON U. Application of multivariate virtual reference feedback tuning for wastewater treatment plant control [J]. Control Engineering Practice, 2012, 20 (5): 499-510.

    [12] 韓京清. 自抗擾控制器及其應(yīng)用 [J]. 控制與決策, 1998, 13 (1): 19-23. DOI: 10.13195/j.cd.1998.01.19.hanjq.004.

    HAN J Q. Auto disturbances rejection controller and it’s applications [J]. Control and Decision, 1998, 13 (1): 19-23. DOI: 10.13195/j.cd.1998. 01.19.hanjq.004.

    [13] 韓京清. 從PID技術(shù)到自抗擾技術(shù) [J]. 控制工程, 2002, 9 (3): 13-18. DOI: 10.14107/j.cnki.kzgc.2002.03.003.

    HAN J Q. From PID to active disturbance rejection control [J]. Control Engineering, 2002, 9 (3): 13-18. DOI: 10.14107/j.cnki.kzgc. 2002.03.003.

    [14] 黃一, 張文革. 自抗擾控制器的發(fā)展 [J].控制理論與應(yīng)用, 2002, 19 (4): 486-492. DOI: 10.3969/j.issn.1000-8152.2002.04.001.

    HUANG Y, ZHANG W G. Development of active disturbance rejection controller [J]. Control Theory and Applications, 2002, 19 (4): 486-492. DOI: 10.3969/j.issn.1000-8152.2002.04.001.

    [15] 高志強. 自抗擾控制思想探究 [J]. 控制理論與應(yīng)用, 2013, 30 (12): 1498-1510. DOI: 10.7641/CTA.2013.31087.

    GAO Z Q. On the foundation of active disturbance rejection control [J]. Control Theory & Applications, 2013, 30 (12): 1498-1510. DOI: 10.7641/CTA.2013.31087.

    [16] Gao Z Q. Scaling and bandwidth-parameterization based controller tuning [C]// Proceedings of the American Control Conference. Denver, 2003: 4989-4996.

    [17] 蘇思賢, 楊慧中. 一類多變量系統(tǒng)的自抗擾非線性動態(tài)解耦控制[J]. 化工學(xué)報, 2010, 61 (8): 1949-1954.

    SU S X, YANG H Z. Active-disturbance-rejection dynamic nonlinear decoupling control for a class of multivariable systems [J]. CIESC Journal, 2010, 61 (8): 1949-1954.

    [18] Huang C E, LI D H, XUE Y L. Active disturbance rejection control for the ALSTOM gasifier benchmark problem [J]. Control Engineering Practice, 2013, 21 (4): 556-564.

    [19] 張園, 孫明瑋, 陳增強. 強制循環(huán)蒸發(fā)系統(tǒng)線性自抗擾解耦控制的魯棒設(shè)計 [J]. 化工學(xué)報, 2015, 66 (S2): 263-270. DOI: 10.11949/ j.issn.0438-157.20141914.

    ZHANG Y, SUN M W, CHEN Z Q. Robust design of linear active disturbance rejection decoupling control for forced-circulationevaporation system [J]. CIESC Journal, 2015, 66 (S2): 263-270. DOI: 10.11949/j.issn.0438-157. 20141914.

    [20] HENZE M, GRADY C P L, GUJER W, et al. Activated sludge model No.1: IAWPRC scientific and technical report No.1 [R]. London: IAW-PRC, 1987.

    [21] TAKACS I, PATRY G, NOLASCO D, et al. A dynamic model of the clarification thickening process [J]. Water Research, 1991, 25 (10): 1263-1271.

    [22] The European Cooperation in the Field of Scientific and Technical Research. Action 624: optimal management of wastewater systems [EB/OL]. http://www.ensic.inpl-nancy. fr/COSTWWTP

    [23] 陳文亮, 姚重華, 呂錫武. 活性污泥工藝的多目標優(yōu)化分析 [J].環(huán)境科學(xué)學(xué)報, 2013, 33 (7): 1918-1925. DOI: 10.13671/j.hjkxxb.2013. 07.025.

    CHEN W L, YAO C H, Lü X W. Analysis of activated sludge process by multi-objective optimization [J]. Acta Scientiae Circumstantiae, 2013, 33 (7): 1918-1925. DOI: 10.13671/j.hjkxxb.2013.07.025.

    [24] SHEN W H, CHEN X Q, PONS M N. Model predictive control for wastewater treatment process with feedforward compensation [J]. Chemical Engineering Journal, 2009, 155 (1/2): 161-174.

    [25] 國家環(huán)境保護總局. 城鎮(zhèn)污水處理廠污染物排放標準: GB 18918—2002 [S]. 北京: 中國標準出版社, 2002.

    Ministry of Environmental Protection of the People’s Republic of China. Discharge standard of pollutants for municipal wastewater treatment plant: GB 18918—2002 [S]. Beijing: Standards Press of China, 2002.

    研究論文

    Received date: 2015-12-22.

    Foundation item: supported by the National Natural Science Foundation of China (61403006), the Importation and Development of High-Caliber Talents Project of Beijing Municipal Institutions (YETP1449, CIT&TCD201404031), the Project of Scientific and Technological Innovation Platform (PXM2015_ 014213_000063) and the Major Project of Beijing Municipal Education Commission Science and Technology Development Plans (KZ201510011011).

    Nitrogen removal in wastewater treatment processes based on linear active disturbance rejection control and its dynamic simulation

    WEI Wei, WANG Xiaoyi, WANG Fan, LIU Zaiwen
    (School of Computer and Information Engineering, Beijing Technology and Business University, Beijing 100048, China)

    Abstract:Based on A/O nitrogen removal process and the kinetic model, Benchmark Simulation Model No.1, two significant factors, nitrate nitrogen and dissolved oxygen, affecting nitrogen removal in wastewater treatment process are considered. Linear active disturbance rejection control (LADRC) is adopted to realize the removal of nitrogen. The control laws designed for nitrate nitrogen and dissolved oxygen do not rely on faithful mathematical models, which make the effluent quality of wastewater be robust enough to uncertainties and fluctuant factors, and nice performance is obtained. Dynamic simulation results show that LADRC has better performance. It can improve the efficiency of nitrogen removal and the effluent quality of wastewater.

    Key words:wastewater; kinetic model; nitrogen removal; LADRC; control; dynamic simulation

    DOI:10.11949/j.issn.0438-1157.20151947

    中圖分類號:TP 273

    文獻標志碼:A

    文章編號:0438—1157(2016)03—1032—08

    基金項目:國家自然科學(xué)基金項目(61403006);北京市屬高等學(xué)校高層次人才引進與培養(yǎng)計劃項目(YETP1449,CIT&TCD201404031);科研基地建設(shè)—科技創(chuàng)新平臺—現(xiàn)代商科特色項目(PXM2015_014213_000063);北京市教委科技計劃重點項目(KZ201510011011)。

    Corresponding author:WEI Wei, weiweizdh@126.com

    猜你喜歡
    動力學(xué)模型控制污水
    什么是核污水
    處理污水
    污水的自述
    小讀者(2019年20期)2020-01-04 02:13:58
    污水罐坍塌
    南方周末(2019-12-05)2019-12-05 11:17:43
    研究模糊控制下車輛的側(cè)傾穩(wěn)定性仿真分析
    房地產(chǎn)投資系統(tǒng)動力學(xué)模型分析
    社會轉(zhuǎn)型期行政權(quán)控制的路徑探索
    科技視界(2016年21期)2016-10-17 17:06:18
    會計預(yù)算控制現(xiàn)狀及方法
    淺談高層建筑沉降監(jiān)測關(guān)鍵環(huán)節(jié)控制
    科技視界(2016年20期)2016-09-29 13:07:14
    保險公司財務(wù)風(fēng)險管理及控制研究
    cao死你这个sao货| 久久性视频一级片| 欧美黑人巨大hd| 国产精品久久视频播放| 在线观看日韩欧美| 久久草成人影院| 99热这里只有精品一区 | 亚洲人成电影免费在线| 中出人妻视频一区二区| 欧美大码av| 看免费av毛片| 亚洲最大成人中文| 每晚都被弄得嗷嗷叫到高潮| 精品不卡国产一区二区三区| 亚洲真实伦在线观看| 欧美中文日本在线观看视频| 亚洲精品乱码久久久v下载方式 | 超碰成人久久| 无限看片的www在线观看| 又黄又爽又免费观看的视频| 亚洲熟女毛片儿| 精品一区二区三区视频在线观看免费| 久久中文看片网| 天天添夜夜摸| tocl精华| 精品99又大又爽又粗少妇毛片 | 久久午夜综合久久蜜桃| 午夜成年电影在线免费观看| 在线观看舔阴道视频| 欧美日韩精品网址| 精品人妻1区二区| 久久国产乱子伦精品免费另类| 国产一区二区在线观看日韩 | 亚洲电影在线观看av| 日本免费一区二区三区高清不卡| 久久精品夜夜夜夜夜久久蜜豆| 亚洲av美国av| 狠狠狠狠99中文字幕| 国产成人影院久久av| 国产淫片久久久久久久久 | 国产高潮美女av| 悠悠久久av| 亚洲成人精品中文字幕电影| www.999成人在线观看| 精品久久久久久久毛片微露脸| 男女之事视频高清在线观看| 日韩有码中文字幕| 蜜桃久久精品国产亚洲av| 午夜久久久久精精品| 窝窝影院91人妻| 国产精品亚洲美女久久久| 又紧又爽又黄一区二区| 男女午夜视频在线观看| 97超视频在线观看视频| 国产精品爽爽va在线观看网站| 搡老妇女老女人老熟妇| 亚洲成人久久爱视频| 91麻豆av在线| 成人一区二区视频在线观看| 波多野结衣巨乳人妻| 中文字幕av在线有码专区| 国产精品一区二区三区四区免费观看 | 午夜福利高清视频| 日本 av在线| 亚洲精品色激情综合| av在线天堂中文字幕| 国产黄色小视频在线观看| 免费观看的影片在线观看| 国内精品久久久久久久电影| 国内久久婷婷六月综合欲色啪| 看免费av毛片| 色老头精品视频在线观看| 国产精品综合久久久久久久免费| 久久精品国产亚洲av香蕉五月| 久久久久久久午夜电影| 一本久久中文字幕| www国产在线视频色| 亚洲在线自拍视频| 99热精品在线国产| 一二三四在线观看免费中文在| 午夜福利在线在线| 深夜精品福利| 亚洲国产欧美一区二区综合| 欧美黄色淫秽网站| 中文字幕av在线有码专区| 日本成人三级电影网站| 午夜福利在线在线| 亚洲中文av在线| 亚洲aⅴ乱码一区二区在线播放| 别揉我奶头~嗯~啊~动态视频| 国产精品亚洲av一区麻豆| 日本一本二区三区精品| 国语自产精品视频在线第100页| 国产精品久久久久久亚洲av鲁大| 亚洲熟妇中文字幕五十中出| а√天堂www在线а√下载| 无人区码免费观看不卡| 夜夜躁狠狠躁天天躁| 色综合婷婷激情| 欧美日韩综合久久久久久 | 精品国产美女av久久久久小说| 国产亚洲欧美在线一区二区| 中亚洲国语对白在线视频| 精品无人区乱码1区二区| 国产综合懂色| 亚洲欧美精品综合一区二区三区| 欧美+亚洲+日韩+国产| 精品国产超薄肉色丝袜足j| 国产精品久久久久久亚洲av鲁大| 亚洲国产欧美人成| 在线观看免费午夜福利视频| 亚洲成人久久爱视频| 亚洲国产精品久久男人天堂| 国产成人精品久久二区二区免费| 国产黄片美女视频| 色哟哟哟哟哟哟| 国产av在哪里看| 热99在线观看视频| 18禁观看日本| 国产主播在线观看一区二区| 在线永久观看黄色视频| 九色成人免费人妻av| 国产精品久久久久久人妻精品电影| 别揉我奶头~嗯~啊~动态视频| 中文字幕最新亚洲高清| 日韩国内少妇激情av| 午夜视频精品福利| 婷婷六月久久综合丁香| av在线天堂中文字幕| 精品一区二区三区视频在线 | 最近在线观看免费完整版| 熟女少妇亚洲综合色aaa.| 久久精品国产综合久久久| 搡老熟女国产l中国老女人| 色综合婷婷激情| 亚洲成人久久性| 热99在线观看视频| 视频区欧美日本亚洲| 俺也久久电影网| 精品久久久久久成人av| 国产伦一二天堂av在线观看| 12—13女人毛片做爰片一| 亚洲,欧美精品.| 夜夜爽天天搞| 俺也久久电影网| 美女高潮的动态| 美女高潮的动态| 日日摸夜夜添夜夜添小说| 身体一侧抽搐| 99久久精品热视频| 男女之事视频高清在线观看| 国产亚洲欧美98| 99国产精品99久久久久| 国内精品美女久久久久久| 男女做爰动态图高潮gif福利片| 精品久久久久久久毛片微露脸| 1024手机看黄色片| 在线播放国产精品三级| 搞女人的毛片| 久久久久久人人人人人| 国产精品一及| 97超视频在线观看视频| 久久99热这里只有精品18| 偷拍熟女少妇极品色| 亚洲黑人精品在线| 欧美3d第一页| av女优亚洲男人天堂 | 亚洲av成人精品一区久久| 白带黄色成豆腐渣| 九九在线视频观看精品| 午夜福利18| 午夜影院日韩av| 亚洲18禁久久av| 99久久无色码亚洲精品果冻| 久久草成人影院| 亚洲无线观看免费| 日韩欧美国产一区二区入口| 男人舔奶头视频| 长腿黑丝高跟| 嫁个100分男人电影在线观看| 免费在线观看成人毛片| 一级毛片女人18水好多| 国产欧美日韩精品一区二区| 精品不卡国产一区二区三区| 18禁黄网站禁片免费观看直播| 国产精品爽爽va在线观看网站| 桃红色精品国产亚洲av| 99热6这里只有精品| 精品久久久久久久人妻蜜臀av| 午夜福利在线观看免费完整高清在 | 97超级碰碰碰精品色视频在线观看| 日韩大尺度精品在线看网址| 成人特级黄色片久久久久久久| 男女做爰动态图高潮gif福利片| 久久亚洲真实| 又爽又黄无遮挡网站| 在线观看舔阴道视频| 午夜福利在线在线| 精品久久久久久久毛片微露脸| 国产成人一区二区三区免费视频网站| 啦啦啦韩国在线观看视频| 国产99白浆流出| 国产精品自产拍在线观看55亚洲| 狂野欧美白嫩少妇大欣赏| 亚洲av成人一区二区三| 亚洲电影在线观看av| 男女床上黄色一级片免费看| 性色avwww在线观看| www.熟女人妻精品国产| 久久国产精品影院| 久久99热这里只有精品18| 欧美成人一区二区免费高清观看 | 男女视频在线观看网站免费| 老司机午夜福利在线观看视频| 啦啦啦观看免费观看视频高清| 精品一区二区三区视频在线观看免费| 深夜精品福利| 1024手机看黄色片| 国产成人欧美在线观看| 亚洲国产日韩欧美精品在线观看 | 麻豆一二三区av精品| 国产午夜福利久久久久久| 国产野战对白在线观看| 亚洲九九香蕉| 日韩av在线大香蕉| 精品久久久久久成人av| av福利片在线观看| 国产精品久久久人人做人人爽| 国产欧美日韩一区二区精品| 欧美最黄视频在线播放免费| 日韩人妻高清精品专区| 亚洲美女视频黄频| 久久久久久久精品吃奶| 色av中文字幕| 欧美激情久久久久久爽电影| avwww免费| 国产精品 国内视频| 久久午夜综合久久蜜桃| 欧美三级亚洲精品| 中文字幕人妻丝袜一区二区| tocl精华| 看片在线看免费视频| 韩国av一区二区三区四区| 97碰自拍视频| 国产精品av久久久久免费| 国产日本99.免费观看| 欧美乱码精品一区二区三区| av天堂在线播放| 18禁美女被吸乳视频| 亚洲国产中文字幕在线视频| 老司机在亚洲福利影院| 国产精品久久久久久精品电影| 美女cb高潮喷水在线观看 | 国产精品久久电影中文字幕| 亚洲欧美日韩东京热| 美女免费视频网站| 成人午夜高清在线视频| 欧美中文综合在线视频| 亚洲欧美精品综合一区二区三区| 国产欧美日韩一区二区三| 男插女下体视频免费在线播放| 亚洲欧洲精品一区二区精品久久久| 国产亚洲av高清不卡| 香蕉av资源在线| 俄罗斯特黄特色一大片| 国产精品亚洲av一区麻豆| 欧美精品啪啪一区二区三区| 精品国产三级普通话版| 国产成人影院久久av| 免费看日本二区| 成人高潮视频无遮挡免费网站| 亚洲国产精品sss在线观看| 精品乱码久久久久久99久播| 久99久视频精品免费| 国产欧美日韩精品一区二区| 日韩中文字幕欧美一区二区| 亚洲av成人av| 99热6这里只有精品| 在线十欧美十亚洲十日本专区| 国产精品久久久久久精品电影| 国产亚洲av嫩草精品影院| 亚洲精品456在线播放app | 亚洲av熟女| 欧洲精品卡2卡3卡4卡5卡区| 日本成人三级电影网站| 国产精品免费一区二区三区在线| 叶爱在线成人免费视频播放| 成年人黄色毛片网站| 99热精品在线国产| 欧美精品啪啪一区二区三区| 一进一出抽搐gif免费好疼| 欧美日韩中文字幕国产精品一区二区三区| 国产美女午夜福利| 国产熟女xx| 亚洲成人中文字幕在线播放| 午夜福利免费观看在线| 美女高潮的动态| 亚洲国产精品合色在线| 亚洲激情在线av| 男女之事视频高清在线观看| 欧美国产日韩亚洲一区| 夜夜看夜夜爽夜夜摸| 亚洲18禁久久av| 日韩有码中文字幕| 欧美乱码精品一区二区三区| 草草在线视频免费看| 成人无遮挡网站| av黄色大香蕉| 色综合婷婷激情| www日本黄色视频网| 久久天躁狠狠躁夜夜2o2o| 亚洲人成网站高清观看| 国产99白浆流出| 国产精品爽爽va在线观看网站| 国产单亲对白刺激| 精品不卡国产一区二区三区| 国产亚洲欧美在线一区二区| 久久草成人影院| 日本熟妇午夜| 久久这里只有精品中国| 两性夫妻黄色片| 麻豆一二三区av精品| 色在线成人网| 国产高清有码在线观看视频| 久久久久精品国产欧美久久久| 欧美日韩黄片免| 国产精品九九99| 天天躁日日操中文字幕| 中文字幕久久专区| 亚洲精品一区av在线观看| 亚洲国产色片| 亚洲中文日韩欧美视频| 国产成人福利小说| 亚洲精品456在线播放app | 噜噜噜噜噜久久久久久91| 日韩大尺度精品在线看网址| 91字幕亚洲| 观看美女的网站| 天堂影院成人在线观看| 国产熟女xx| 久久亚洲精品不卡| 亚洲精品久久国产高清桃花| 麻豆一二三区av精品| 大型黄色视频在线免费观看| 一夜夜www| 久久热在线av| 好男人电影高清在线观看| 俄罗斯特黄特色一大片| 午夜成年电影在线免费观看| 久久久久久人人人人人| 中文字幕熟女人妻在线| 精品99又大又爽又粗少妇毛片 | 国产黄a三级三级三级人| bbb黄色大片| 757午夜福利合集在线观看| 国产精品一及| 老司机午夜十八禁免费视频| 成年版毛片免费区| 十八禁网站免费在线| 操出白浆在线播放| 成人高潮视频无遮挡免费网站| 欧美日韩乱码在线| 级片在线观看| 国产野战对白在线观看| 18禁裸乳无遮挡免费网站照片| 在线播放国产精品三级| 日韩免费av在线播放| 国产精品女同一区二区软件 | 国产精品国产高清国产av| 久久人人精品亚洲av| 亚洲无线在线观看| 在线观看日韩欧美| 婷婷亚洲欧美| 色综合站精品国产| 一边摸一边抽搐一进一小说| 日日摸夜夜添夜夜添小说| 精品国产三级普通话版| av女优亚洲男人天堂 | 免费一级毛片在线播放高清视频| 成人av一区二区三区在线看| 偷拍熟女少妇极品色| 国产 一区 欧美 日韩| 久久久久久九九精品二区国产| 天堂动漫精品| 亚洲一区高清亚洲精品| 麻豆久久精品国产亚洲av| av视频在线观看入口| 中文字幕最新亚洲高清| 神马国产精品三级电影在线观看| 天天躁日日操中文字幕| 国产成人精品久久二区二区免费| 久久这里只有精品中国| 成人午夜高清在线视频| 国产av麻豆久久久久久久| 国产一级毛片七仙女欲春2| 国产精品av久久久久免费| 午夜激情福利司机影院| 天天一区二区日本电影三级| 黄片小视频在线播放| 在线免费观看不下载黄p国产 | 久久亚洲精品不卡| 亚洲色图 男人天堂 中文字幕| 欧美乱码精品一区二区三区| 97超视频在线观看视频| 欧美极品一区二区三区四区| 一进一出抽搐gif免费好疼| 久久中文字幕人妻熟女| 丰满人妻熟妇乱又伦精品不卡| 亚洲 国产 在线| 亚洲精品久久国产高清桃花| 91在线观看av| 18禁国产床啪视频网站| 国产亚洲欧美98| 又紧又爽又黄一区二区| 90打野战视频偷拍视频| 久久人妻av系列| 一级a爱片免费观看的视频| 日本一二三区视频观看| 亚洲最大成人中文| 色综合站精品国产| 最近最新中文字幕大全电影3| 免费看a级黄色片| 最近最新免费中文字幕在线| 成年女人毛片免费观看观看9| 好男人电影高清在线观看| 久久天躁狠狠躁夜夜2o2o| 亚洲自偷自拍图片 自拍| 男女午夜视频在线观看| 天天躁日日操中文字幕| 国产1区2区3区精品| 国产乱人视频| 色精品久久人妻99蜜桃| 在线国产一区二区在线| 一进一出好大好爽视频| 成熟少妇高潮喷水视频| 午夜福利在线在线| 久久久色成人| 欧美日韩国产亚洲二区| 国产野战对白在线观看| av女优亚洲男人天堂 | 成人精品一区二区免费| 999精品在线视频| 狠狠狠狠99中文字幕| 午夜日韩欧美国产| 99在线人妻在线中文字幕| 国内久久婷婷六月综合欲色啪| 日韩欧美精品v在线| 别揉我奶头~嗯~啊~动态视频| 中文亚洲av片在线观看爽| 精品国产三级普通话版| 亚洲精品456在线播放app | 免费在线观看影片大全网站| xxxwww97欧美| 亚洲一区二区三区色噜噜| 国产麻豆成人av免费视频| 亚洲真实伦在线观看| 狂野欧美白嫩少妇大欣赏| 国产一区二区三区视频了| 中国美女看黄片| 国产亚洲精品综合一区在线观看| 精品乱码久久久久久99久播| svipshipincom国产片| 美女免费视频网站| www国产在线视频色| 精品人妻1区二区| 久99久视频精品免费| 亚洲中文字幕一区二区三区有码在线看 | 少妇人妻一区二区三区视频| 两个人视频免费观看高清| 99久久99久久久精品蜜桃| 亚洲中文字幕日韩| 母亲3免费完整高清在线观看| 中文资源天堂在线| 国产午夜精品久久久久久| 国产精品 欧美亚洲| 日本免费a在线| 成人欧美大片| 成在线人永久免费视频| 757午夜福利合集在线观看| 又大又爽又粗| 久久久久久九九精品二区国产| 国产亚洲欧美98| 亚洲国产看品久久| 久久久精品大字幕| av在线天堂中文字幕| 成年版毛片免费区| 香蕉国产在线看| 欧美+亚洲+日韩+国产| 欧美xxxx黑人xx丫x性爽| 夜夜看夜夜爽夜夜摸| 久久这里只有精品19| 麻豆成人av在线观看| 精品一区二区三区视频在线 | 免费在线观看影片大全网站| 夜夜爽天天搞| 两性午夜刺激爽爽歪歪视频在线观看| 久久久精品欧美日韩精品| 成人一区二区视频在线观看| 午夜免费成人在线视频| 成年女人毛片免费观看观看9| 欧美+亚洲+日韩+国产| 国产成人aa在线观看| 九九在线视频观看精品| 国产精品久久视频播放| 国产精品99久久久久久久久| 国内毛片毛片毛片毛片毛片| 亚洲 欧美一区二区三区| 2021天堂中文幕一二区在线观| 一级作爱视频免费观看| 午夜影院日韩av| 欧美zozozo另类| 淫妇啪啪啪对白视频| 精品国产乱码久久久久久男人| 中国美女看黄片| 国产极品精品免费视频能看的| 午夜精品在线福利| 国产精品综合久久久久久久免费| 免费观看的影片在线观看| 人人妻人人澡欧美一区二区| 制服人妻中文乱码| 精品久久久久久久毛片微露脸| 国产精品99久久99久久久不卡| 麻豆成人午夜福利视频| 老熟妇仑乱视频hdxx| 亚洲五月婷婷丁香| 国产极品精品免费视频能看的| 欧美黑人欧美精品刺激| 亚洲专区字幕在线| 我要搜黄色片| 亚洲午夜理论影院| 日本 欧美在线| 欧美日韩精品网址| 18美女黄网站色大片免费观看| 日韩免费av在线播放| 一进一出抽搐动态| 嫩草影视91久久| 国产精品永久免费网站| 日本 av在线| 色尼玛亚洲综合影院| 麻豆国产97在线/欧美| 免费观看精品视频网站| 一二三四社区在线视频社区8| 亚洲欧美精品综合久久99| 美女午夜性视频免费| 午夜福利成人在线免费观看| 午夜福利高清视频| 精品99又大又爽又粗少妇毛片 | 国产精品自产拍在线观看55亚洲| 窝窝影院91人妻| 99久久无色码亚洲精品果冻| 黄色日韩在线| 在线观看午夜福利视频| 精品人妻1区二区| 桃色一区二区三区在线观看| 亚洲自拍偷在线| 男女视频在线观看网站免费| 久久精品人妻少妇| 国产成人av教育| 999久久久国产精品视频| 久久伊人香网站| 国产欧美日韩一区二区三| 最新在线观看一区二区三区| 丰满人妻一区二区三区视频av | 国产亚洲精品久久久com| 精品人妻1区二区| xxx96com| 性色av乱码一区二区三区2| 欧美日韩黄片免| 一卡2卡三卡四卡精品乱码亚洲| 亚洲精品粉嫩美女一区| 日韩大尺度精品在线看网址| cao死你这个sao货| 国产精品99久久99久久久不卡| 国内精品久久久久久久电影| 日本与韩国留学比较| 美女cb高潮喷水在线观看 | 母亲3免费完整高清在线观看| 婷婷亚洲欧美| 亚洲黑人精品在线| 精品人妻1区二区| 99久国产av精品| 午夜福利欧美成人| 日韩欧美国产在线观看| 国产野战对白在线观看| 国产97色在线日韩免费| 久久午夜亚洲精品久久| 亚洲精品在线美女| 精品日产1卡2卡| 午夜福利免费观看在线| 日本黄色视频三级网站网址| 少妇的丰满在线观看| 麻豆国产97在线/欧美| 久久国产精品影院| 欧美中文日本在线观看视频| 19禁男女啪啪无遮挡网站| 一个人看的www免费观看视频| 97超级碰碰碰精品色视频在线观看| 亚洲18禁久久av| 网址你懂的国产日韩在线| 桃色一区二区三区在线观看| 长腿黑丝高跟| 在线观看66精品国产| 国产视频内射| 欧美黑人欧美精品刺激| 日日摸夜夜添夜夜添小说| 亚洲中文av在线| 久久国产精品影院| 欧美成狂野欧美在线观看| 无遮挡黄片免费观看| 亚洲精品乱码久久久v下载方式 | 久久久水蜜桃国产精品网| 人人妻人人看人人澡| 一本一本综合久久| 国产亚洲精品久久久久久毛片| 免费在线观看成人毛片|