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

    ASM2d模型改良及在DE型氧化溝中的應(yīng)用

    2017-03-29 22:39:25郭彥雪李偉趙凱
    土木建筑與環(huán)境工程 2017年1期
    關(guān)鍵詞:數(shù)學(xué)模型污水處理

    郭彥雪++李偉++趙凱

    摘要:對國際水質(zhì)污染與控制協(xié)會(IAWQ)推出的活性污泥2d號模型(ASM2d)進(jìn)行優(yōu)化,使之能夠模擬含有菌膠團(tuán)的污水處理工藝。使用改良后的ASM2d模型模擬并分析了有機(jī)物、氮和磷在DE型雙溝式氧化溝中的降解過程,發(fā)現(xiàn)優(yōu)化后模型能夠較好地?cái)M合這些物質(zhì)的降解過程。依據(jù)所建立的模型計(jì)算出各組分在DE型氧化溝中1個周期4個階段的轉(zhuǎn)化速率,該結(jié)果有助于深入分析有機(jī)物、氮及磷不同形態(tài)之間的轉(zhuǎn)化規(guī)律。

    關(guān)鍵詞:DE型氧化溝;污水組分;數(shù)學(xué)模型;污水處理

    中圖分類號:X703.1文獻(xiàn)標(biāo)志碼:A文章編號:16744764(2017)01012507

    收稿日期:20160405

    基金項(xiàng)目:國家自然科學(xué)基金(51308437);西安建筑科技大學(xué)人才基金(RC1247)

    作者簡介:郭彥雪(1986),女,主要從事水質(zhì)模擬研究,(Email)guoyanxue1128@163.com。

    李偉(通信作者),男,副教授,碩士生導(dǎo)師,(Email) liwei@xauat.edu.cn。

    Received:20160405

    Foundation item:National Natural Science Foundation of China (NSFC)(No. 51308437);Science Foundation of Xi'an University of Architecture and Technology(No. RC1247)

    Author brief:Guo Yanxue (1986), main research interest: water quality simulation, (Email) guoyanxue1128@163.com.

    Li Wei (corresponding author), associate professor, master supervisor, (Email) liwei@xauat.edu.cn.Modified ASM2d model and application in the DEoxidation ditch

    Guo Yanxue1a, Li Wei1b,Zhao Kai2

    (1a. School of Management; 1b.School of Environmental and Municipal Engineering,Xi'an University

    of Architecture and Technology,Xi'an 710055, P.R. China; 2. Beishiqiao Wastewater Treatment

    Corporation, Xian Capital Environmental Protection Group Company Limited, Xi'an 710086, P.R. China)

    Abstract:Activated sludge model 2d (ASM2d) proposed by IAWQ was modified,and it was used to simulate wastewater treatment containing biomass flocs. Combined with the operation parameters of a doublechannel DEoxidation ditch, a calculated model was built to simulate carbon,nitrogen and phosphorus degradation processes. The results showed that it could fit the experiment data well. The removal rates of pollutant in the four phases in a cycle were also calculated. And the analyses provided an insight in the nitrification,denitrification and phosphorus removal mechanisms in such operation system,which was useful for the system diagnosis and optimization.

    Keywords:DEoxidation ditch;wastewater characterization;mathematical model;wastewater treatment

    由于污水組分和活性污泥法的復(fù)雜性,必須開展數(shù)學(xué)模型研究,才能對其進(jìn)行更加深入研究[1]。目前,常用的模型為國際水質(zhì)污染與控制協(xié)會(IAWQ)推出的活性污泥2d號模型(ASM2d),各國的專家與學(xué)者將該數(shù)學(xué)模型和實(shí)際污水處理工藝相結(jié)合建立了各種計(jì)算模型,用于預(yù)測污染物降解和指導(dǎo)污水廠的運(yùn)行[25]。根據(jù)報(bào)道,ASM2d模型能夠較好的模擬污染物在厭氧區(qū)/缺氧區(qū)/好氧區(qū)存在明顯邊界的生物反應(yīng)器(如A2/O工藝)中的降解過程,而對邊界不明顯的反應(yīng)器(如氧化溝)中降解過程的擬合結(jié)果往往與實(shí)測值相差較大。高揚(yáng)等[4]使用ASM2d模型模擬氮在改良型A2/O工藝中降解過程,發(fā)現(xiàn)出水模擬值與實(shí)測值偏差僅為054%~1.24%;然而,Xie等[6]使用相同模型對氮在Carrousel氧化溝中降解過程進(jìn)行模擬,發(fā)現(xiàn)模擬值與實(shí)測值的平均偏差達(dá)8%~28%。在實(shí)際污水處理工藝中(如DE氧化溝),細(xì)菌會與水中物質(zhì)結(jié)合形成菌膠團(tuán),溶解氧濃度在菌膠團(tuán)中由外而內(nèi)逐漸降低,形成溶解氧濃度梯度。當(dāng)外部溶解氧濃度相對較低時,菌膠團(tuán)內(nèi)部會形成缺氧區(qū)甚至是厭氧區(qū),從而會發(fā)生同步硝化反硝化(SND)脫氮現(xiàn)象,據(jù)報(bào)道氧化溝中50%以上的氮是依靠該過程去除的[7]。然而,ASM2d模型并未包含SND過程,因此,難以準(zhǔn)確的擬合污染物在存在SND工藝中的降解過程。

    本文依據(jù)菌膠團(tuán)的結(jié)構(gòu),提出“厭氧核”概念,并以此為依據(jù)優(yōu)化了ASM2d模型。使用優(yōu)化后的數(shù)學(xué)模型,結(jié)合某DE型雙溝式氧化溝的實(shí)際運(yùn)行特點(diǎn),建立了DE型雙溝式氧化溝中碳氧化、硝化和反硝化以及磷去除的計(jì)算模型。

    1材料與方法

    1.1污水廠簡介

    污水處理量為12萬m3/d。進(jìn)水經(jīng)過粗格柵、細(xì)格柵、曝氣撇油沉砂池、厭氧選擇池、DE型雙溝式氧化溝和終沉池后排放至皂河。

    污泥回流比為80%,回流污泥濃度約12 000 mg/L,污泥停留時間(SRT)為18 d,平均水力停留時間(HRT)為4.75 h。DE型雙溝式氧化溝的運(yùn)行方式是兩溝交替工作循環(huán)曝氣,1個周期歷時240 min,分4個階段進(jìn)行。DE型氧化溝1個周期4個階段的具體運(yùn)行方式見圖1。

    圖1DE型氧化溝4個階段運(yùn)行狀態(tài)圖

    Fig.1 Operation state in DE oxidation ditch1.2測量項(xiàng)目與方法

    由于DE型雙溝式氧化溝I溝第1、2階段與II溝3、4階段成鏡像關(guān)系,因此僅對I溝的4個階段進(jìn)行了測定。

    具體方式如下:在氧化溝I溝上固定取樣點(diǎn)(如圖1所示),以第1階段開端為始點(diǎn),10 min為一時間間隔進(jìn)行取樣(水力測試表明污水在氧化溝中循環(huán)1周的時間約為10 min)。取樣后立即測定混合液的溶解氧濃度(DO)和水溫(T); 然后用中速濾紙對混合液進(jìn)行過濾,在實(shí)驗(yàn)室測定濾液的以下項(xiàng)目:pH值、化學(xué)需氧量(COD)、總氮(TN)、氨氮(NH+4±N)、硝酸鹽氮(NO-3—N)、總磷(TP)和正磷酸鹽(PO3-4—P)。此外,實(shí)驗(yàn)還測定了SV30、MLSS和MLVSS以確定污泥的性質(zhì),以上各項(xiàng)目均按照國家標(biāo)準(zhǔn)方法進(jìn)行[8]。

    1.3污水特征分析

    具體氧化溝的污水進(jìn)水和出水水質(zhì)組分分析結(jié)果見表1。

    表1DE氧化溝系統(tǒng)進(jìn)出水水質(zhì)分析(其中進(jìn)水取自曝氣沉砂池出口,出水取自氧化溝出水)

    Table 1Analysis of influent and effluent sampled from DEoxidation ditch system符號名稱數(shù)值單位符號名稱數(shù)值單位SO2溶解氧0g O2·m-3CODtot,inf進(jìn)水總COD360g COD·m-3SF易生物降解有機(jī)物79.2g COD·m-3CODf,,tot進(jìn)水濾液中COD135.2g COD·m-3SA發(fā)酵產(chǎn)物 (醋酸)32.4g COD·m-3CODtot,eff出水濾液中COD25.2g COD·m-3SNH4氨氮38.8g N·m-3CODf,eff出水COD24.8g COD·m-3SNO3硝態(tài)氮2.6g N·m-3CODVFAsVFAs轉(zhuǎn)化的COD79.2g COD·m-3SPO4正磷酸鹽3.8g P·m-3SI不可降解有機(jī)物23.6g COD·m-3SALK堿度5mol HCO- 3·m-3續(xù)表1符號名稱數(shù)值單位符號名稱數(shù)值單位顆粒態(tài)組分SA=CODVFAsXI不可降解有機(jī)物23.6g COD·m-3Ss=CODf,tot- SIXS慢速可降解有機(jī)物156g COD·m-3SF=SS-SAXI=(1-α)CODsusp,infXH異養(yǎng)菌32.1g COD·m-3SA=CODVFAsXS=αCODsusp,infXPAO聚磷菌0g COD·m-3SI=0.95CODf,effα= [(BODtot/(1- YH)- SS)]/CODsusp,infXPP聚磷菌體內(nèi)正磷酸0g P·m-3Ss=CODf,tot- SIBODtot = BOD5/(1- e-5k)XPHA聚磷菌體內(nèi)有機(jī)物0g COD·m-3SA=CODVFAsCODsusp,inf = CODtot,inf- CODfilt,infXAUT硝化菌0g COD·m-3SI=0.95CODf,effXTSS顆粒物質(zhì)211.7g TSS·m-3

    2數(shù)學(xué)模型的建立

    2.1模型建立的假設(shè)

    在模型建立時做了如下假設(shè):1)進(jìn)水量不變;2)原污水中組分及含量恒定;3)每立方米混合液進(jìn)入氧化溝后直接與一份污水完全混合;4)空間平面上為理想推流式;5)系統(tǒng)運(yùn)行時溫度不變;6)pH值恒定;7)氧化溝內(nèi)無毒害微生物的物質(zhì);8)微生物對顆粒狀有機(jī)物的捕捉是瞬時進(jìn)行的;9)有機(jī)物和有機(jī)氮的水解是同時進(jìn)行的,且速率相等;10)污水在氧化溝內(nèi)的流動速度恒定。

    2.2模型建立的方法

    水力試驗(yàn)發(fā)現(xiàn)溝內(nèi)混合液平均流速約為0.5 m/s,因此在模擬中假定流速恒定為0.5 m/s。在氧化溝長度方向上以0.5 m為一個單位進(jìn)行分割,每一單位的污水看成一個處理單元,袁宏林等[9]通過測定發(fā)現(xiàn),從寬度和深度來看,DE氧化溝具有完全混合式反應(yīng)器的特點(diǎn),因此,每一單元的污水水質(zhì)在空間上可以近似認(rèn)為是均勻的。從整個系統(tǒng)長度方向上看,氧化溝系統(tǒng)又具有推流式反應(yīng)器的特點(diǎn),因此我們假定取出一單元混合液,讓其在氧化溝中循環(huán)(如圖2所示)。

    2.3建立模型中使用的物料平衡模型

    根據(jù)質(zhì)量守恒方程可以得到式(1),由水量平衡可以得到式(2),ASM2d給出的系統(tǒng)反應(yīng)量由化學(xué)劑量數(shù)vij和組分的過程速率ρj乘積之和表示,得到式(3)。d(V[t]·Ci[t])dt=Qin[t]·Cini[t]-

    Qout[t]·Ci[t]+ri·V[t](1)

    dV[t]dt=Qin[t]-Qout[t](2)

    ri=jvijρj(3)式中:V[t] 為 t時刻污水的體積;Ci[t]為t時刻i組分的濃度, g/ m3;Qin[t]為 t時刻進(jìn)水的流量;Cini[t]為t時刻進(jìn)水i組分的濃度,g/ m3;Qout[t]為 t時刻出水流量;ri為反應(yīng)速率,g/(0.1 min(m3);νij為化學(xué)計(jì)量系數(shù);ρj為過程速率式。

    圖2污水單元分割與反應(yīng)速率示意圖

    Fig.2Wastewater divided unit and reaction rate如圖2所示,上述方程在不同階段可進(jìn)行如下簡化。

    進(jìn)水段:氧化溝的進(jìn)水有兩種方式,一是進(jìn)水口進(jìn)水,二是在固定時刻由穿孔進(jìn)水。此時Qout[t] =0,若t時刻進(jìn)水量為Qin[t],則方程可以簡化為dCi[t]dt=Qin[t]·Cin[t]+V[t]·Ci[t]Qin[t]+V[t]+ri

    dV[t]dt=Qin[t](4)式中:C0為曝氣沉砂池出水口處污染物濃度;Ci,II為II溝進(jìn)水前1秒這份污水中污染物的濃度。

    反應(yīng)段:氧化溝在不進(jìn)水也不出水時認(rèn)為其處于反應(yīng)段,此時Qout[t]=0,Qin[t]=0,方程簡化為dCi[t]dt=ri[t]dV[t]dt=0(5)出水段:該系統(tǒng)出水有兩種方式:一是出水口出水,二是固定時刻通過穿孔出水,此時Qin[t]=0,方程簡化為dCi[t]dt=(V[t]-Qout[t])·Ci[t]V[t]+ri

    dV[t]dt=-Qout[t](6)2.4動力學(xué)參數(shù)的確定

    有研究表明測定動力學(xué)參數(shù)并不需要逐個確定,而是選擇一些對模型影響顯著的系數(shù),通過實(shí)驗(yàn)方法或者模擬值擬合的方法進(jìn)行確定[1012]。在模擬過程中,我們采用ASM2d模型中推薦的動力學(xué)參數(shù)。

    3ASM2d模型優(yōu)化

    使用ASM2d模型模擬有機(jī)物、氮和磷在氧化溝中降解過程時,NO-3—N和PO3-4—P的模擬值與實(shí)測值具有顯著差別(見圖3),這種差異可能是由于ASM2d模型沒有考慮SND作用的影響引起的[1315]。

    由于DO在菌膠團(tuán)外層和內(nèi)層的濃度不同,以此為依據(jù)提出“厭氧核”的概念。

    如圖4所示,假定存在一個薄層將菌膠團(tuán)分為兩個部分,分別為“厭氧核”和“好氧外層”,其所占菌膠團(tuán)體積比例分別為α和(1-α)。其中“好氧外層”中DO濃度與菌膠團(tuán)外部混合液中DO濃度相同,“厭氧核”中DO濃度采用式(7)進(jìn)行計(jì)算。SO2,f=SO2kO2SO2+kO2(7)式中:SO2,f為菌膠團(tuán)內(nèi)部DO濃度;SO2為混合液中DO濃度;kO2為溶解氧擴(kuò)散系數(shù)(本模型取0.10)。

    此時,總的反應(yīng)速率可以采用式(8)計(jì)算。rT=α×rSO2,f+(1-α)×rSO2(8)式中:rT為每種組分的總反應(yīng)速率;rSO2,f為菌膠團(tuán)內(nèi)反應(yīng)速率; rSO2為好氧外層和混合液中反應(yīng)速率。

    研究中,當(dāng)α取值為0.68時,DO的模擬值與實(shí)際測量值之間的平均偏差僅為0.3%;4種污染指標(biāo)COD、NH+4—N、NO-3—N和PO3-4—P的模擬值與實(shí)際測量值之間的平均偏差為2.3%、1.9%、5.1%和8.7%,表明優(yōu)化后的ASM2d模型能夠較好地?cái)M合DO和各污染物的濃度(見圖3)。4個階段各組分的變化速率見表2。圖3溶解氧與污染物在DE氧化溝

    中降解過程的實(shí)測值與模擬值變化

    Fig. 3The variation of measured value and simulation

    of DO and pollutants in a cycle of DE oxidant ditch圖4菌膠團(tuán)分層示意圖

    Fig. 4Proposed configuration of biomass flocs表2物質(zhì)SF, SA, NOx-N, NH+4—N和PO3-4—P在不同階段的平均降解速率

    Table 2Average reaction rates of SF, SA, NOx-N, NH+4—N and PO3-4—P in each phaseg/(h·m3)階段SFSASNH4SNO3SPO412341234123412341234141.2 364.5 301.1 218.2 0.4 3.6 3.0 2.2 2402.0 178.1 230.6 305.9 4.0 1.8 2.3 3.1 367.1 79.7 23.2 26.0 0.7 0.8 0.2 0.3 4-145.6 -412.2 -439.0 -648.4 -2.0 -5.7 -6.0 -8.9 -0.4 -1.0 -1.1 -1.6 5-30.4 -48.1 -7.0 -95.1 -1.3 -2.1 -0.3 -4.2 -0.4 -0.6 -0.1 -1.2 6-497.2 -178.6 -70.1 -169.1 -6.8 -2.5 -1.0 -2.3 -65.2 -23.4 -9.2 -22.2 -1.2 -0.4 -0.2 -0.4 7-106.5 -45.5 -1.1 -24.7 -4.7 -2.0 0.0 -1.1 -14.0 -6.0 -0.1 -3.2 -1.3 -0.6 0.0 -0.3 8-255.5 -108.2 -37.7 -91.4 255.5 108.2 37.7 91.4 7.7 3.2 1.1 2.7 2.6 1.1 0.4 0.9 910.6 10.6 10.6 10.6 3.3 3.3 3.3 3.3 10-379.2 -87.3 -32.5 -240.7 151.7 34.9 13.0 96.3 11-38.1 -143.7 -61.5 -111.4 12-9.3 -3.6 -2.5 -6.5 -133.1 -51.7 -35.1 -93.2 13-1.2 -9.9 -8.0 -6.4 -0.4 -2.8 -2.3 -1.8 14-4.4 -3.6 -4.4 -4.9 -13.1 -10.7 -13.2 -14.7 -1.2 -1.0 -1.3 -1.4 152.9 2.9 2.9 2.9 0.9 0.9 0.9 0.9 1685.9 86.1 86.3 86.2 175.7 6.0 3.6 3.8 18-8.9 -89.6 -66.3 -53.0 8.8 88.1 65.2 52.1 0.0 -0.4 -0.3 -0.3 190.4 0.4 0.4 0.4 0.1 0.1 0.1 0.1 Tot-388.0 -76.78.2-358.8-246.1-66.70.7-265.4-2.7 -92.0 -65.5 -58.7 -92.8 44.5 40.3 5.568.3 -75.9 2.1 -23.9 注:1.好氧水解過程; 2.缺氧水解過程; 3.厭氧水解過程; 4. 異養(yǎng)菌利用SF生長; 5.異養(yǎng)菌利用SA生長; 6.利用SF反硝化; 7.利用SA反硝化; 8.SF發(fā)酵; 9.異養(yǎng)菌衰亡; 10.聚磷菌儲存PHA; 11.聚磷菌好氧儲存正磷酸鹽; 12.聚磷菌缺氧儲存正磷酸鹽; 13. 聚磷菌好氧生長; 14.聚磷菌缺氧生長; 15.聚磷菌衰亡;16.正磷酸鹽釋放; 17.胞內(nèi)有機(jī)物釋放; 18.好氧菌生長; 19.好氧菌衰亡。4模擬結(jié)果與分析

    4.1第1階段

    此階段轉(zhuǎn)刷停止,混合液中DO濃度較低(0.2~0.3 mg/L,圖3(a)),是典型的缺氧狀態(tài)。在該條件下,硝化細(xì)菌生長速率緩慢,由表2可知,NH+4—N的平均去除效率僅為2.7 g/(h·m3),而進(jìn)水中帶入了較高濃度的NH+4—N,致使其在氧化溝內(nèi)積累,濃度不斷升高。進(jìn)水中含有的大量易生物降解的有機(jī)物SF和SA,為反硝化反應(yīng)提供了足夠能量,缺氧條件為反硝化細(xì)菌提供了合適的生存環(huán)境,因而反硝化速率高達(dá)92.8 g/(h·m3),其結(jié)果是混合液中NO-3—N濃度大幅下降。

    在較低DO條件下,聚磷菌(PAOs)開始分解體內(nèi)的PHB。由表2可知,本階段釋放PO3-4—P的速率高達(dá)68.3 g/(h·m3)。同時,氧化溝進(jìn)水中也帶來了較高濃度的PO3-4—P(見表2),致使本階段PO3-4—P濃度大幅升高。

    盡管進(jìn)水中帶入了大量有機(jī)物,但濾液中COD濃度基本保持不變,其原因是反硝化脫氮(過程6,見表2)和聚磷菌分解體內(nèi)PHA(過程10)需要消耗大量的溶解態(tài)有機(jī)物SF和SA(其去除速率分別達(dá)到497.2和379.2 g/(h·m3)),將進(jìn)水中帶入的和顆粒態(tài)有機(jī)物分解生成的溶解態(tài)有機(jī)物快速消耗,致使氧化溝混合液中溶解態(tài)COD基本保持不變。

    4.2第2階段

    本階段開始曝氣,DO濃度不斷升高。同時,由表2可知大量顆粒狀可降解有機(jī)物質(zhì)被水解為SF(好氧水解和缺氧水解速率達(dá)364.5和178.1 g/(h·m3)),為硝化細(xì)菌的生長提供了足夠的能量。在硝化細(xì)菌的快速生長過程中,NH+4—N被快速消耗,本階段平均硝化速率為92 g/(h·m3),高于其余3個階段,且沒有NH+4—N的額外補(bǔ)充,其結(jié)果是NH+4—N濃度直線下降。然而,本階段的NO-3—N濃度并沒有相應(yīng)的大幅升高。由表2可知,盡管混合液中DO濃度較高,但本階段反硝化速率依然可達(dá)23.4 g/(h·m3)(過程6)。這是由于混合液中的DO沒有穿透菌膠團(tuán),即在菌膠團(tuán)內(nèi)部依然處于厭氧或缺氧狀態(tài),出現(xiàn)了同步硝化反硝化脫氮現(xiàn)象[16]。

    在好氧條件下,PAOs會過量吸收正磷酸鹽,以胞內(nèi)聚合物的形式儲存。由表2可知,本階段好氧存儲PP和缺氧儲存PP的速率高達(dá)143.7和51.7 g/(h·m3),致使混合液中正磷酸鹽濃度快速下降。

    4.3第3階段

    DO處于較高水平,保持在1.5 mg/L左右,處于好氧狀態(tài)。本階段含有大量顆粒狀有機(jī)物、少量溶解態(tài)有機(jī)物的污水從II溝通過穿孔進(jìn)入I溝,因此,溶解態(tài)COD濃度基本沒有變化。由表2可知,本階段異養(yǎng)細(xì)菌快速生長所需的有機(jī)物來自顆粒狀有機(jī)物的好氧水解和缺氧水解,其速率分別為3011和230.6 g/(h·m3)。

    由II溝進(jìn)入的混合液中含有大量氨氮,硝化細(xì)菌在好氧、能量充足的條件下進(jìn)行快速硝化作用,將其帶入的NH+4—N去除,因此本階段NH+4—N濃度維持在較低水平。與第2階段相似,在同步硝化反硝化的作用下,NO-3—N濃度僅有小幅升高。

    由II溝流入I溝的混合液中含有較高的磷酸鹽濃度,但I(xiàn)溝處于好氧狀態(tài),PAOs在這種條件下將磷酸鹽合成PHA儲存[17],因此正磷酸鹽在第3階段維持在較低濃度。

    4.4第4階段

    含有高濃度有機(jī)物、NH+4—N和正磷酸鹽的污水重新進(jìn)入I溝,同時DO依然維持在較高的水平。由表2可知,SF和SA降解速率高于水解速率,然而進(jìn)水中含有較高濃度的SF和SA,因此溶解態(tài)COD濃度有小幅度升高。

    由表2可知,本階段硝化速率低于第2、3階段,不足以將進(jìn)水中高濃度的NH+4—N完全氧化,因此會有小幅度的積累。在同步硝化反硝化的作用下,硝酸鹽氮濃度基本維持不變。

    與第2、3階段類似,在好氧、有機(jī)物充足的條件下,PAOs依然能夠?qū)⒒旌弦褐姓姿猁}合成PHA儲存體內(nèi),因此正磷酸鹽濃度依然維持在較低水平。

    5結(jié)論

    根據(jù)菌膠團(tuán)的結(jié)構(gòu),提出“厭氧核”概念,并以此為基礎(chǔ)優(yōu)化了ASM2d模型,使之能夠應(yīng)用于模擬含有菌膠團(tuán)的污水處理工藝。將優(yōu)化的ASM2d模型結(jié)合某DE型雙溝式氧化溝的實(shí)際運(yùn)行特點(diǎn),建立了其計(jì)算模型,該模型可以很好的對氧化溝內(nèi)COD、NH+4—N、硝酸鹽氮和正磷酸鹽濃度的變化進(jìn)行模擬。根據(jù)計(jì)算結(jié)果深度分析并總結(jié)了DE氧化溝1個周期4個階段在處理污染物所起的作用,該結(jié)果對于氧化溝工藝改良具有一定意義。

    參考文獻(xiàn):

    [1] HENZE M. Activated sludge models ASM1, ASM2, ASM2d and ASM3. Vol. 9 [M].IWA publishing, 2000.

    [2] YANG S, GUO W, CHEN Y, et al. Simultaneous insitu sludge reduction and nutrient removal in an A2 MOM system: Performances, mechanisms, and modeling with an extended ASM2d model [J].Water Research, 2016, 88:524537.

    [3] MOHAMMED R N, LU X W. Modeling of nitritation and denitritation in a novel PITSFSEU process using an extension of ASM2d model [J]. Desalination and Water Treatment, 2015, 53(7):18621880.

    [4] 高揚(yáng),邵青. 基于ASM2D模型對龍王嘴污水處理廠的工藝模擬[J]. 給水排水,2014,40(2):117121.

    GAO Y, SHAO Q. Simulation of the process in Longwangzui wastewater treatment plant using ASM2D model [J]. Water & Wastewater Engineering, 2014, 40(2):117121. (in Chinese)

    [5] 陳文亮,呂錫武,姚重華,等. 基于 ASM2d 的數(shù)學(xué)模型在六箱一體化工藝分析中的應(yīng)用[J]. 環(huán)境科學(xué)研究,2014,27(8):929935.

    CHEN W L, LYU X W, YAO H H, et al. Application of mathematical model for analysis of sixtank integrative activated sludge process [J]. Research of Environmental Sciences, 2014, 27(8):929935. (in Chinese)

    [6] XIE W M, ZHANG R, LI W W, et al. Simulation and optimization of a fullscale carrousel oxidation ditch plant for municipal wastewater treatment [J]. Biochemical Engineering Journal, 2011, 56(1): 916.

    [7] LIU Y, SHI H, XIA L, et al. Study of operational conditions of simultaneous nitrification and denitrification in a carrousel oxidation ditch for domestic wastewater treatment [J]. Bioresource Technology, 2010, 101(3):901906.

    [8] 國家環(huán)境保護(hù)局編委會.水與廢水監(jiān)測分析方法[M].4版. 北京:中國環(huán)境科學(xué)出版社,2002.

    Environmental protection administration editorial board. Methods for the examination of water and wastewater[M].4th Edition.Beijing:China Environmental Science Press,2002.(in Chinese)

    [9] 袁宏林,王曉昌,李志華,等. DE型氧化溝系統(tǒng)處理城市污水的效果[J]. 給水排水,2000,26(4):68.

    YUAN H L, WANG X C, LI Z H, et al. DE oxidation ditch treating urban wastewater [J]. Water & Wastewater Engineering, 2000, 26(4):68. (in Chinese)

    [10] GIWA A, HASAN S W. Numerical modeling of an electrically enhanced membrane bioreactor (MBER) treating mediumstrength wastewater [J]. Journal of Environmental Management, 2015, 164:19.

    [11] BRUN R, KHNI M, SIEGRIST H, et al. Practical identifiability of ASM2d parameterssystematic selection and tuning of parameter subsets [J]. Water Research, 2002, 36(16): 41134127.

    [12] MIKOSZ J. Determination of permissible industrial pollution load at a municipal wastewater treatment plant [J]. International Journal of Environmental Science and Technology, 2015, 12(3): 827836.

    [13] ZHOU X, HAN Y, GUO X. Identification and evaluation of SND in a fullscale multichannel oxidation ditch system under different aeration modes [J]. Chemical Engineering Journal, 2015, 259:715723.

    [14] HE Q, PENG X, LI Z. The treatment of animal manure wastewater by coupled simultaneous methanogenesis and denitrification (SMD) and shortcut nitrificationdenitrification (SND) [J]. Journal of Chemical Technology and Biotechnology, 2014, 89(11):16971704.

    [15] LUO W, JIN X, YU Y, et al. Efficient nitrogen removal via simultaneous nitrification and denitrification in a penicillin wastewater biological treatment plant [J]. Environmental Technology, 2014, 35(22):28852893.

    [16] LI Y, ZOU J, ZHANG L, et al. Aerobic granular sludge for simultaneous accumulation of mineral phosphorus and removal of nitrogen via nitrite in wastewater [J]. Bioresource Technology, 2014, 154:178184.

    [17] MIELCZAREK A T, NGUYEN H T T, NIELSEN J L, et al. Population dynamics of bacteria involved in enhanced biological phosphorus removal in Danish wastewater treatment plants [J]. Water Research, 2013, 47(4):15291544.

    猜你喜歡
    數(shù)學(xué)模型污水處理
    AHP法短跑數(shù)學(xué)模型分析
    活用數(shù)學(xué)模型,理解排列組合
    農(nóng)村黑臭水體治理和污水處理淺探
    太原市61個村要建污水處理設(shè)施嚴(yán)禁直排入河
    淺談醫(yī)療污水處理要求
    基于電力機(jī)器人控制系統(tǒng)的數(shù)學(xué)模型簡述
    電子制作(2019年23期)2019-02-23 13:21:34
    涌浪機(jī)在污水處理中的應(yīng)用
    工業(yè)污水處理的再利用
    對一個數(shù)學(xué)模型的思考
    古塔形變的數(shù)學(xué)模型
    河南科技(2014年5期)2014-02-27 14:08:47
    av福利片在线| 999精品在线视频| 天堂俺去俺来也www色官网| 村上凉子中文字幕在线| 岛国视频午夜一区免费看| 男女下面插进去视频免费观看| 99精品久久久久人妻精品| 国产日韩一区二区三区精品不卡| 老司机亚洲免费影院| 老司机靠b影院| 久久精品aⅴ一区二区三区四区| 琪琪午夜伦伦电影理论片6080| 欧美乱妇无乱码| 国产精品 欧美亚洲| 一本大道久久a久久精品| 日韩精品青青久久久久久| 欧美黑人欧美精品刺激| 真人一进一出gif抽搐免费| 精品久久久久久成人av| 色婷婷久久久亚洲欧美| 精品人妻1区二区| 日本黄色日本黄色录像| 99精品久久久久人妻精品| 久久99一区二区三区| 国产麻豆69| 免费观看精品视频网站| 超碰97精品在线观看| 国产免费av片在线观看野外av| 天堂中文最新版在线下载| 日韩一卡2卡3卡4卡2021年| 国产成人欧美在线观看| 亚洲人成伊人成综合网2020| 久久这里只有精品19| 久久九九热精品免费| 欧美日韩乱码在线| 国产精品二区激情视频| 久久香蕉国产精品| 久久香蕉激情| 看片在线看免费视频| 国产精品98久久久久久宅男小说| 一级a爱片免费观看的视频| 成人免费观看视频高清| 亚洲人成伊人成综合网2020| 亚洲一区二区三区不卡视频| 久久国产乱子伦精品免费另类| 色综合站精品国产| 操美女的视频在线观看| 9191精品国产免费久久| 欧美日本亚洲视频在线播放| 亚洲av电影在线进入| 欧美人与性动交α欧美精品济南到| 天堂影院成人在线观看| 天堂俺去俺来也www色官网| 1024香蕉在线观看| 两性夫妻黄色片| 国产高清videossex| 一边摸一边抽搐一进一出视频| 成人三级黄色视频| 亚洲五月天丁香| 校园春色视频在线观看| 亚洲欧美精品综合一区二区三区| a级片在线免费高清观看视频| 午夜福利在线观看吧| 人妻久久中文字幕网| 国产成+人综合+亚洲专区| 好男人电影高清在线观看| a级片在线免费高清观看视频| 视频在线观看一区二区三区| 婷婷六月久久综合丁香| 欧美激情高清一区二区三区| 日韩精品中文字幕看吧| 黄片播放在线免费| 中国美女看黄片| 久久精品国产亚洲av香蕉五月| 国产精品免费视频内射| 久久久国产成人精品二区 | 丰满迷人的少妇在线观看| 制服诱惑二区| 日韩欧美一区二区三区在线观看| 国产视频一区二区在线看| 人人妻人人添人人爽欧美一区卜| 欧美在线一区亚洲| 国产主播在线观看一区二区| 中文字幕精品免费在线观看视频| 男女高潮啪啪啪动态图| 久久久国产一区二区| 制服诱惑二区| 在线观看免费午夜福利视频| 国产精品久久视频播放| 精品国产乱码久久久久久男人| 亚洲中文av在线| 人人妻人人爽人人添夜夜欢视频| 亚洲一区二区三区色噜噜 | 亚洲国产欧美网| 欧美黄色淫秽网站| 少妇裸体淫交视频免费看高清 | 成人18禁高潮啪啪吃奶动态图| a在线观看视频网站| 亚洲色图av天堂| 夜夜爽天天搞| 国产精品爽爽va在线观看网站 | 一级作爱视频免费观看| 日本a在线网址| 天天躁夜夜躁狠狠躁躁| 亚洲人成电影观看| 国产精品野战在线观看 | 午夜福利免费观看在线| 国产97色在线日韩免费| 亚洲 国产 在线| tocl精华| 亚洲中文日韩欧美视频| 日本 av在线| 极品人妻少妇av视频| 亚洲午夜精品一区,二区,三区| 国产精品香港三级国产av潘金莲| 午夜精品在线福利| 又黄又粗又硬又大视频| 成年人免费黄色播放视频| 丰满迷人的少妇在线观看| 涩涩av久久男人的天堂| 久久久久久久久久久久大奶| 亚洲一区中文字幕在线| 成人免费观看视频高清| 99国产精品一区二区蜜桃av| 最近最新中文字幕大全电影3 | 极品教师在线免费播放| 巨乳人妻的诱惑在线观看| 99久久99久久久精品蜜桃| 多毛熟女@视频| 成人影院久久| 久久影院123| 看免费av毛片| 91精品国产国语对白视频| 日日爽夜夜爽网站| av国产精品久久久久影院| 精品久久久久久久久久免费视频 | 国产精品野战在线观看 | 天堂√8在线中文| 18禁国产床啪视频网站| 日韩欧美一区二区三区在线观看| 精品乱码久久久久久99久播| 天堂动漫精品| 亚洲av日韩精品久久久久久密| 久久久久久大精品| videosex国产| 免费观看人在逋| 中文字幕人妻丝袜制服| 淫秽高清视频在线观看| 黄色视频不卡| 久久香蕉国产精品| 亚洲av日韩精品久久久久久密| 久久久精品欧美日韩精品| 国产一区二区在线av高清观看| 欧美激情 高清一区二区三区| 亚洲美女黄片视频| 中文亚洲av片在线观看爽| 一边摸一边做爽爽视频免费| 日韩精品青青久久久久久| 亚洲国产精品sss在线观看 | 国产一区在线观看成人免费| 88av欧美| 久久精品影院6| 日韩欧美一区二区三区在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 男人舔女人的私密视频| 国产无遮挡羞羞视频在线观看| 一级a爱视频在线免费观看| 老司机亚洲免费影院| 精品国产乱子伦一区二区三区| 久久影院123| 丰满迷人的少妇在线观看| 日本三级黄在线观看| 亚洲欧美激情综合另类| bbb黄色大片| 在线观看免费午夜福利视频| 69精品国产乱码久久久| 操出白浆在线播放| 精品熟女少妇八av免费久了| 久久九九热精品免费| 香蕉国产在线看| 男人操女人黄网站| 久久精品亚洲av国产电影网| 视频在线观看一区二区三区| 亚洲黑人精品在线| 满18在线观看网站| 男女之事视频高清在线观看| cao死你这个sao货| 久久人妻福利社区极品人妻图片| 亚洲国产精品一区二区三区在线| 国产黄色免费在线视频| 99在线视频只有这里精品首页| 99久久精品国产亚洲精品| 人妻久久中文字幕网| 黄色丝袜av网址大全| 怎么达到女性高潮| 超碰成人久久| 亚洲欧美日韩另类电影网站| 在线观看免费视频日本深夜| 亚洲伊人色综图| 国产又色又爽无遮挡免费看| 久久精品国产清高在天天线| 国产91精品成人一区二区三区| 久久久水蜜桃国产精品网| 99在线人妻在线中文字幕| 久久中文字幕人妻熟女| 日日干狠狠操夜夜爽| 免费搜索国产男女视频| 麻豆国产av国片精品| 操出白浆在线播放| 亚洲男人的天堂狠狠| 欧美日韩国产mv在线观看视频| 又紧又爽又黄一区二区| 亚洲成人久久性| 新久久久久国产一级毛片| 天天躁夜夜躁狠狠躁躁| 午夜福利免费观看在线| 中文字幕人妻丝袜一区二区| 19禁男女啪啪无遮挡网站| 99在线人妻在线中文字幕| 久久香蕉国产精品| 国产免费现黄频在线看| 久久中文字幕人妻熟女| а√天堂www在线а√下载| 狠狠狠狠99中文字幕| 欧美在线一区亚洲| 亚洲一区二区三区欧美精品| 嫁个100分男人电影在线观看| 国产av又大| 免费在线观看亚洲国产| 国产亚洲欧美精品永久| av天堂久久9| 夜夜躁狠狠躁天天躁| 久久精品国产综合久久久| 超碰成人久久| av免费在线观看网站| 国产成人欧美| 乱人伦中国视频| 99re在线观看精品视频| 久热这里只有精品99| 老汉色av国产亚洲站长工具| 国产一区二区三区综合在线观看| 国产三级黄色录像| 黄频高清免费视频| 这个男人来自地球电影免费观看| 成人黄色视频免费在线看| 最新在线观看一区二区三区| 精品午夜福利视频在线观看一区| 多毛熟女@视频| 一级毛片女人18水好多| 国产激情欧美一区二区| 真人做人爱边吃奶动态| 国产精品乱码一区二三区的特点 | 国产精品野战在线观看 | 亚洲国产看品久久| 在线永久观看黄色视频| 亚洲精品美女久久av网站| av网站免费在线观看视频| 99香蕉大伊视频| 成年女人毛片免费观看观看9| 成人影院久久| 91精品国产国语对白视频| 曰老女人黄片| 国产男靠女视频免费网站| 久久亚洲精品不卡| 国产成人精品在线电影| 一边摸一边抽搐一进一小说| 女人被狂操c到高潮| xxx96com| 亚洲午夜精品一区,二区,三区| 男人舔女人的私密视频| 在线观看免费视频网站a站| av有码第一页| 色精品久久人妻99蜜桃| 丰满人妻熟妇乱又伦精品不卡| 精品福利观看| 波多野结衣av一区二区av| 国产欧美日韩精品亚洲av| 亚洲精品中文字幕一二三四区| 亚洲一码二码三码区别大吗| 免费少妇av软件| 国产一区二区三区综合在线观看| 电影成人av| 国产成人欧美在线观看| 午夜免费鲁丝| av电影中文网址| 丰满饥渴人妻一区二区三| 夜夜看夜夜爽夜夜摸 | 91麻豆av在线| 国产精品 欧美亚洲| 欧美久久黑人一区二区| 久久精品国产综合久久久| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲中文日韩欧美视频| 热re99久久国产66热| 又黄又粗又硬又大视频| 两性午夜刺激爽爽歪歪视频在线观看 | 18禁黄网站禁片午夜丰满| 国产97色在线日韩免费| 亚洲国产中文字幕在线视频| 欧美日韩国产mv在线观看视频| 看片在线看免费视频| 亚洲欧美一区二区三区久久| 成年版毛片免费区| 亚洲欧美激情在线| 久久欧美精品欧美久久欧美| 久久婷婷成人综合色麻豆| 免费看a级黄色片| 欧美在线黄色| 久久天堂一区二区三区四区| 亚洲成av片中文字幕在线观看| 岛国视频午夜一区免费看| 日韩精品免费视频一区二区三区| 两个人免费观看高清视频| 长腿黑丝高跟| 极品教师在线免费播放| 久久99一区二区三区| 久久午夜亚洲精品久久| 久久久久久亚洲精品国产蜜桃av| 成人亚洲精品av一区二区 | a级毛片黄视频| 亚洲欧美精品综合一区二区三区| 亚洲色图av天堂| 999久久久精品免费观看国产| 18禁观看日本| 在线视频色国产色| 宅男免费午夜| 天天添夜夜摸| 一级,二级,三级黄色视频| 人人妻人人澡人人看| 麻豆国产av国片精品| aaaaa片日本免费| 色综合欧美亚洲国产小说| 91大片在线观看| 一区二区三区国产精品乱码| 国产精品1区2区在线观看.| 悠悠久久av| 日韩 欧美 亚洲 中文字幕| 岛国视频午夜一区免费看| 99国产综合亚洲精品| www.www免费av| 久久久精品国产亚洲av高清涩受| 精品无人区乱码1区二区| 高潮久久久久久久久久久不卡| 国产日韩一区二区三区精品不卡| av国产精品久久久久影院| 国产成+人综合+亚洲专区| 悠悠久久av| 久久狼人影院| 88av欧美| 亚洲精品国产精品久久久不卡| 免费人成视频x8x8入口观看| 成人黄色视频免费在线看| 自拍欧美九色日韩亚洲蝌蚪91| 久久午夜综合久久蜜桃| 精品卡一卡二卡四卡免费| 变态另类成人亚洲欧美熟女 | 国产片内射在线| 精品高清国产在线一区| 国产精品98久久久久久宅男小说| 午夜精品国产一区二区电影| 国产一区在线观看成人免费| 国产一区二区三区综合在线观看| 啪啪无遮挡十八禁网站| 久99久视频精品免费| 久久久精品国产亚洲av高清涩受| 高清在线国产一区| 韩国av一区二区三区四区| 丰满饥渴人妻一区二区三| 亚洲国产精品999在线| 99国产精品免费福利视频| 狠狠狠狠99中文字幕| 国产熟女午夜一区二区三区| 久久人人爽av亚洲精品天堂| 水蜜桃什么品种好| 亚洲欧美日韩无卡精品| 欧美日本中文国产一区发布| 亚洲av熟女| 国产成人av教育| 黑人欧美特级aaaaaa片| 岛国在线观看网站| av在线播放免费不卡| 极品教师在线免费播放| 日韩有码中文字幕| 国产精品 欧美亚洲| 免费日韩欧美在线观看| 日韩成人在线观看一区二区三区| 男女下面插进去视频免费观看| 精品日产1卡2卡| 久久青草综合色| 国产精品亚洲一级av第二区| 伊人久久大香线蕉亚洲五| a级片在线免费高清观看视频| 国产精品自产拍在线观看55亚洲| 午夜精品在线福利| 成人18禁高潮啪啪吃奶动态图| 在线十欧美十亚洲十日本专区| 日韩精品青青久久久久久| 国产精品二区激情视频| 交换朋友夫妻互换小说| 国产亚洲欧美98| 怎么达到女性高潮| 激情视频va一区二区三区| 精品国产一区二区三区四区第35| 精品久久久久久久毛片微露脸| av福利片在线| 精品国内亚洲2022精品成人| 丰满饥渴人妻一区二区三| 久久天堂一区二区三区四区| 亚洲黑人精品在线| 久久精品国产99精品国产亚洲性色 | 国产成人欧美在线观看| av国产精品久久久久影院| 国产精品98久久久久久宅男小说| 人人妻,人人澡人人爽秒播| 国产一卡二卡三卡精品| 91国产中文字幕| 久久人妻av系列| 亚洲五月色婷婷综合| 久久精品aⅴ一区二区三区四区| 亚洲av第一区精品v没综合| 精品久久久久久电影网| 脱女人内裤的视频| 精品国产超薄肉色丝袜足j| 国产91精品成人一区二区三区| 亚洲一码二码三码区别大吗| www国产在线视频色| 热re99久久国产66热| 亚洲专区中文字幕在线| 欧美激情极品国产一区二区三区| tocl精华| 看免费av毛片| 成人特级黄色片久久久久久久| 国产91精品成人一区二区三区| 视频在线观看一区二区三区| 欧美另类亚洲清纯唯美| 变态另类成人亚洲欧美熟女 | 亚洲欧美一区二区三区久久| 黑人巨大精品欧美一区二区蜜桃| 久久香蕉国产精品| 热99re8久久精品国产| 可以在线观看毛片的网站| 国产精品综合久久久久久久免费 | 免费观看精品视频网站| 精品国内亚洲2022精品成人| 天天躁狠狠躁夜夜躁狠狠躁| 丁香欧美五月| 亚洲久久久国产精品| 一级片免费观看大全| 夫妻午夜视频| 丰满人妻熟妇乱又伦精品不卡| 丝袜美足系列| 伊人久久大香线蕉亚洲五| 成人亚洲精品一区在线观看| 久久久国产成人免费| 99国产精品99久久久久| 好男人电影高清在线观看| 天堂√8在线中文| 午夜免费激情av| 亚洲情色 制服丝袜| 成人特级黄色片久久久久久久| 国产蜜桃级精品一区二区三区| 丝袜美腿诱惑在线| 美女 人体艺术 gogo| 成人影院久久| 欧美成人午夜精品| 女性生殖器流出的白浆| 三级毛片av免费| 黄网站色视频无遮挡免费观看| 最近最新中文字幕大全免费视频| av国产精品久久久久影院| av在线播放免费不卡| 久久久国产成人免费| 日韩免费高清中文字幕av| 一级a爱片免费观看的视频| 51午夜福利影视在线观看| 中出人妻视频一区二区| 热99国产精品久久久久久7| 精品一品国产午夜福利视频| 午夜亚洲福利在线播放| 在线观看免费午夜福利视频| 激情在线观看视频在线高清| 久久精品影院6| 国产麻豆69| 日韩国内少妇激情av| 国产又色又爽无遮挡免费看| 免费看十八禁软件| 夜夜躁狠狠躁天天躁| 村上凉子中文字幕在线| 无遮挡黄片免费观看| 午夜福利,免费看| 久久久国产成人免费| 女性被躁到高潮视频| 宅男免费午夜| 黄片大片在线免费观看| 我的亚洲天堂| 婷婷精品国产亚洲av在线| 18美女黄网站色大片免费观看| 丰满迷人的少妇在线观看| 天天躁夜夜躁狠狠躁躁| 精品免费久久久久久久清纯| 亚洲avbb在线观看| 超碰97精品在线观看| av有码第一页| 两人在一起打扑克的视频| 一边摸一边抽搐一进一出视频| 欧美黑人欧美精品刺激| 三级毛片av免费| 久久精品亚洲精品国产色婷小说| 久久这里只有精品19| 免费在线观看视频国产中文字幕亚洲| 一区福利在线观看| 一区二区三区精品91| 一本大道久久a久久精品| 日本五十路高清| 国产精品av久久久久免费| av网站在线播放免费| av免费在线观看网站| 国产精品一区二区三区四区久久 | 亚洲久久久国产精品| 69精品国产乱码久久久| 欧美乱色亚洲激情| 男女下面进入的视频免费午夜 | 又大又爽又粗| 欧美最黄视频在线播放免费 | 欧美在线黄色| 日韩大码丰满熟妇| 在线观看免费高清a一片| 国产欧美日韩一区二区三| 麻豆国产av国片精品| 精品福利观看| 国产精品爽爽va在线观看网站 | 午夜老司机福利片| 欧美激情久久久久久爽电影 | 亚洲全国av大片| 香蕉丝袜av| 伦理电影免费视频| 国产成人系列免费观看| 一区二区三区精品91| 一本大道久久a久久精品| 啦啦啦免费观看视频1| 99久久国产精品久久久| 午夜激情av网站| 日韩三级视频一区二区三区| 国产亚洲欧美精品永久| 新久久久久国产一级毛片| 1024视频免费在线观看| 亚洲色图综合在线观看| 久久狼人影院| 麻豆av在线久日| 欧美大码av| 久久婷婷成人综合色麻豆| 欧美日韩乱码在线| av天堂久久9| 波多野结衣高清无吗| a级毛片在线看网站| 国产伦一二天堂av在线观看| 老汉色av国产亚洲站长工具| 国产激情久久老熟女| 婷婷精品国产亚洲av在线| 一级,二级,三级黄色视频| 丝袜美足系列| 在线观看一区二区三区| 久9热在线精品视频| 欧美日韩av久久| 欧美日韩中文字幕国产精品一区二区三区 | x7x7x7水蜜桃| 麻豆成人av在线观看| 日韩视频一区二区在线观看| 亚洲精品国产色婷婷电影| 亚洲少妇的诱惑av| 电影成人av| 麻豆一二三区av精品| 国产真人三级小视频在线观看| 制服诱惑二区| 日韩精品中文字幕看吧| 啦啦啦 在线观看视频| 午夜久久久在线观看| 日韩欧美国产一区二区入口| 99riav亚洲国产免费| 欧美最黄视频在线播放免费 | 久久久久久亚洲精品国产蜜桃av| 91字幕亚洲| 99在线视频只有这里精品首页| 丰满饥渴人妻一区二区三| 1024视频免费在线观看| 咕卡用的链子| 国产av在哪里看| 亚洲黑人精品在线| 黄色a级毛片大全视频| 久久久水蜜桃国产精品网| 男女高潮啪啪啪动态图| 精品人妻1区二区| 后天国语完整版免费观看| 亚洲国产精品合色在线| 久久久久久大精品| 国产精品一区二区免费欧美| 精品久久久久久,| 制服诱惑二区| 亚洲欧美精品综合一区二区三区| 丝袜美腿诱惑在线| 中文字幕高清在线视频| 欧美黑人欧美精品刺激| 国产单亲对白刺激| 亚洲国产精品sss在线观看 | 夜夜爽天天搞| 曰老女人黄片| 国产精品自产拍在线观看55亚洲| 最近最新中文字幕大全电影3 | 国产伦一二天堂av在线观看| 亚洲精品美女久久av网站| 国产成人啪精品午夜网站| 国产成人欧美| 亚洲一区二区三区欧美精品| 老司机深夜福利视频在线观看| 国产三级黄色录像|