亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        多段進(jìn)水A/O生物接觸氧化工藝模擬

        2021-11-10 08:56:26朱繼濤冷峻彤王子捷
        關(guān)鍵詞:總氮氧池水流量

        吳 沛,朱繼濤,顧 鑫,張 凱,冷峻彤,王子捷,胡 博*

        (1. 長(zhǎng)安大學(xué) 建筑工程學(xué)院,陜西 西安 710061; 2. 長(zhǎng)安大學(xué) 住房和城鄉(xiāng)建設(shè)部給水排水重點(diǎn)實(shí)驗(yàn)室,陜西 西安 710061)

        0 引 言

        多段進(jìn)水缺氧/好氧(Anoxic/Oxic,A/O)生物脫氮工藝是一種基于傳統(tǒng)A/O工藝改進(jìn)的生物脫氮工藝[1]。多段進(jìn)水A/O生物脫氮工藝由多個(gè)缺氧和好氧反應(yīng)單元交替變化組成,進(jìn)水分多點(diǎn)進(jìn)入工藝的缺氧池, 為缺氧池發(fā)生的反硝化過(guò)程提供碳源[2]。工藝不同反應(yīng)單元內(nèi)的微生物濃度呈階梯分布[3]。進(jìn)水碳氮比(C/N)、進(jìn)水流量分配、缺氧池和好氧池的容積分配等都是影響多段進(jìn)水A/O生物脫氮工藝運(yùn)行性能的關(guān)鍵因素[2-3]。多段進(jìn)水A/O生物脫氮工藝可采用生物膜法和活性污泥法運(yùn)行。與活性污泥法相比,生物膜法由于無(wú)需污泥回流、不存在污泥膨脹等問(wèn)題,在運(yùn)行和管理上都更為方便[4]。在生物膜工藝中,底物的降解過(guò)程受生化反應(yīng)動(dòng)力學(xué)、液-固界面?zhèn)髻|(zhì)等過(guò)程影響,這些過(guò)程彼此相互聯(lián)系又相互制約[5],使得污染物在生物膜中的降解機(jī)制較為復(fù)雜,工藝設(shè)計(jì)只能依賴經(jīng)驗(yàn)公式或半經(jīng)驗(yàn)公式[6]。

        隨著國(guó)際水協(xié)會(huì)推出活性污泥系列數(shù)學(xué)模型,數(shù)學(xué)模型已成為一種重要的研究手段被廣泛應(yīng)用于污水生物處理工藝的數(shù)值模擬和優(yōu)化設(shè)計(jì)[7-8]。活性污泥1號(hào)數(shù)學(xué)模型(ASM1)于1987年提出,包含快速可生物利用基質(zhì)、慢速可生物利用基質(zhì)、自養(yǎng)菌、異養(yǎng)菌、氨氮、硝態(tài)氮等13個(gè)反應(yīng)組分,以及異養(yǎng)菌好氧生長(zhǎng)、異養(yǎng)菌缺氧生長(zhǎng)、自養(yǎng)菌好氧生長(zhǎng)、異養(yǎng)菌和自養(yǎng)菌衰減等8個(gè)生化反應(yīng)過(guò)程,可描述好氧和缺氧條件下的碳氧化過(guò)程以及自養(yǎng)硝化過(guò)程和異養(yǎng)反硝化過(guò)程[9]。生物膜數(shù)學(xué)模型的研究始于20世紀(jì)80年代,模型發(fā)展經(jīng)歷了從簡(jiǎn)單到復(fù)雜,從一維到多維,從穩(wěn)態(tài)到動(dòng)態(tài)的過(guò)程。目前,一維連續(xù)模型、DLA模型(Diffusion-limited Aggregation Model)、細(xì)胞自動(dòng)機(jī)模型等都已被應(yīng)用于生物膜工藝的數(shù)值模擬[10]。在生物膜工藝中,污染物需要先通過(guò)擴(kuò)散過(guò)程進(jìn)入生物膜內(nèi)部,與微生物接觸后,才能被利用和降解。因此,在生物膜系統(tǒng)中,除了要考慮基質(zhì)限制之外,還需要考慮擴(kuò)散過(guò)程對(duì)污染物降解過(guò)程的影響。對(duì)于建模而言,在滿足需求的前提下,模型應(yīng)做到盡可能地簡(jiǎn)單,以便于求解和應(yīng)用[11]。

        本研究基于ASM1模型,通過(guò)將描述物質(zhì)擴(kuò)散過(guò)程的Fick擴(kuò)散定律與ASM1模型結(jié)合,采用生物膜零級(jí)反應(yīng)動(dòng)力學(xué)的半級(jí)反應(yīng)解作為生物膜系統(tǒng)生化反應(yīng)的過(guò)程速率,構(gòu)建可用于模擬多段進(jìn)水A/O生物膜脫氮工藝的動(dòng)力學(xué)模型,并采用三段進(jìn)水A/O生物接觸氧化工藝對(duì)模型進(jìn)行校驗(yàn)和驗(yàn)證;最終,采用建立的動(dòng)力學(xué)模型對(duì)多段進(jìn)水A/O生物接觸氧化工藝的進(jìn)水流量分配進(jìn)行優(yōu)化,以期能夠提升工藝的運(yùn)行性能。

        1 材料與方法

        1.1 生物膜零級(jí)反應(yīng)動(dòng)力學(xué)模型的構(gòu)建

        假設(shè)生物膜為理想生物膜,并對(duì)生物膜工藝的反應(yīng)過(guò)程進(jìn)行以下簡(jiǎn)化:①忽略氨化過(guò)程,假設(shè)微生物衰減產(chǎn)生的有機(jī)氮直接轉(zhuǎn)化為氨氮;②化學(xué)需氧量(COD)對(duì)硝化過(guò)程有抑制作用;③好氧池與缺氧池的泥齡不同;④生物膜中發(fā)生的生化反應(yīng)過(guò)程條件及模型的邊界條件與ASM1模型一致。模型組分包括易生物降解基質(zhì)(SS)、慢速可生物降解基質(zhì)(XS)、異養(yǎng)菌(XB,H)、自養(yǎng)菌(XB,A)、微生物衰減的顆粒態(tài)產(chǎn)物(XP)、溶解氧(SO)、硝態(tài)氮(SNO)和氨氮(SNH)。模型的生化反應(yīng)過(guò)程包括慢速可生物降解有機(jī)物水解過(guò)程、異養(yǎng)菌好氧生長(zhǎng)過(guò)程、異養(yǎng)菌缺氧生長(zhǎng)過(guò)程、異養(yǎng)菌衰減過(guò)程、自養(yǎng)菌好氧生長(zhǎng)過(guò)程和自養(yǎng)菌衰減過(guò)程。由于溶解性不可生物降解性基質(zhì)(SI)和顆粒態(tài)不可生物降解基質(zhì)(XI)對(duì)生化過(guò)程無(wú)影響,為了簡(jiǎn)化模型,在構(gòu)建模型中不考慮這兩個(gè)組分。此外,由于進(jìn)水中堿度充足,模型組分中也不考慮堿度(SALK)對(duì)生化反應(yīng)的影響。此外,由于氨化速率難以獲取[12],為了簡(jiǎn)化模型中氮的轉(zhuǎn)化過(guò)程,故氨化過(guò)程也從模型中略去。

        在生物膜工藝中,液相中的基質(zhì)會(huì)通過(guò)擴(kuò)散作用進(jìn)入生物膜內(nèi),微生物利用擴(kuò)散到生物膜的基質(zhì)實(shí)現(xiàn)自身增殖。根據(jù)基質(zhì)進(jìn)入生物膜深度的不同,基質(zhì)可分為完全穿透生物膜基質(zhì)和部分穿透生物膜基質(zhì)。對(duì)于完全穿透生物膜基質(zhì),基質(zhì)可通過(guò)擴(kuò)散作用抵達(dá)載體表面;對(duì)于部分穿透生物膜基質(zhì),基質(zhì)僅能穿透一定厚度的生物膜,并不能完全穿透生物膜抵達(dá)載體表面[13]。對(duì)于部分穿透生物膜基質(zhì),生物膜中發(fā)生生化反應(yīng)速率的大小將受到傳質(zhì)阻力的影響和限制。穩(wěn)態(tài)條件下,生物膜內(nèi)距離膜表面z處的基質(zhì)濃度表達(dá)式[14]為

        (1)

        式中:z為生物膜內(nèi)距離膜表面的距離(m);S為生物膜內(nèi)距離膜表面z處的基質(zhì)濃度(g·m-3);Df為生物膜內(nèi)的基質(zhì)擴(kuò)散速率(m2·d-1);rf為反應(yīng)速率(g·m-3·d-1)。

        基于生物膜零級(jí)反應(yīng)動(dòng)力學(xué),以及穩(wěn)態(tài)條件下生物膜內(nèi)污染物的降解量與通過(guò)生物膜表面進(jìn)入生物膜的基質(zhì)相等,生物膜內(nèi)基質(zhì)降解速率計(jì)算公式為

        (2)

        根據(jù)式(2),假設(shè)生物膜中的微生物由異氧菌和自養(yǎng)菌組成,異養(yǎng)菌所占比例為fH,自養(yǎng)菌為1-fH??紤]電子受體濃度對(duì)反應(yīng)速率的影響,多基質(zhì)多種類微生物生物膜反應(yīng)動(dòng)力學(xué)模型的系數(shù)矩陣見(jiàn)表1,模型的反應(yīng)過(guò)程速率見(jiàn)表2。除擴(kuò)散系數(shù)外,模型中其他化學(xué)計(jì)量系數(shù)與動(dòng)力學(xué)參數(shù)的意義及取值均與ASM1模型[9]相同。

        表1 模型化學(xué)計(jì)量矩陣

        表2 模型的反應(yīng)過(guò)程速率

        1.2 試驗(yàn)裝置

        三段進(jìn)水A/O生物接觸氧化工藝的試驗(yàn)裝置采用有機(jī)玻璃制成[15](圖1)。采用組合填料作為微生物生長(zhǎng)的載體,總體積為44 L,將其等分為4部分,每部分體積為11 L(圖2)。最后一段用作沉淀池,其余3段作為反應(yīng)段。每一反應(yīng)段均由一個(gè)缺氧池(A池)和一個(gè)好氧池(O池)組成,缺氧池與好氧池的總體積比為2∶3,其中3個(gè)反應(yīng)段內(nèi)缺氧池與好氧池之間、各段之間均使用隔板分隔。進(jìn)水流量Q為0.916 7 L·h-1,水力停留時(shí)間為48 h,其中反應(yīng)段的水力停留時(shí)間為36 h,沉淀池的水力停留時(shí)間為12 h。各段進(jìn)水流量分配系數(shù)r1、r2和r3均為0.333 3,各段進(jìn)水流量分別以Q1、Q2和Q3進(jìn)入各段缺氧池。同時(shí),設(shè)置回流系統(tǒng),將經(jīng)沉淀池沉淀后的出水回流至第1段缺氧池,回流比(R)分別為0.50、0.75和1.00。

        圖1 試驗(yàn)裝置

        圖2 三段進(jìn)水A/O生物接觸氧化工藝圖

        1.3 試驗(yàn)水質(zhì)

        試驗(yàn)進(jìn)水為人工模擬的高濃度生活污水,進(jìn)水水質(zhì)見(jiàn)表3,配水組分見(jiàn)表4。有機(jī)碳、氮和磷分別由C6H12O6·H2O、NH4HCO3和KH2PO4提供。進(jìn)水中加入適量NaHCO3補(bǔ)充堿度,微量元素的添加與Cheong等的研究[16]相同。進(jìn)水COD的組成采用耗氧速率(OUR)間歇實(shí)驗(yàn)進(jìn)行測(cè)定,其中SS占COD的比例為95%,XS占COD的比例為5%。

        表3 試驗(yàn)進(jìn)水水質(zhì)

        表4 試驗(yàn)配水組分

        1.4 分析項(xiàng)目及分析方法

        1.5 模型的驗(yàn)證和校驗(yàn)

        經(jīng)過(guò)約3個(gè)月的運(yùn)行,三段進(jìn)水A/O生物膜工藝各池生物膜的厚度、生物量及出水水質(zhì)基本達(dá)到穩(wěn)定。由于好氧池中氣流和水流對(duì)生物膜的沖刷作用強(qiáng)于缺氧池,好氧池生物膜的脫附速率大于缺氧池生物膜,好氧池生物膜的泥齡短于缺氧池生物膜。采用回流比為0.5時(shí)反應(yīng)器達(dá)到穩(wěn)定處理效果的結(jié)果對(duì)模型進(jìn)行驗(yàn)證,并采用回流比為1.0時(shí)的反應(yīng)器運(yùn)行結(jié)果對(duì)模型進(jìn)行校驗(yàn)。模型特征參數(shù)取值見(jiàn)表5。

        表5 模型特征參數(shù)的取值

        1.6 進(jìn)水流量分配系數(shù)的優(yōu)化

        以試驗(yàn)水質(zhì)(COD為1 000 mg·L-1,硝態(tài)氮濃度為100 mg·L-1)為模型的進(jìn)水水質(zhì),采用校驗(yàn)后的模型對(duì)內(nèi)回流比為1.0時(shí),不同流量分配系數(shù)下工

        藝的運(yùn)行結(jié)果進(jìn)行模擬。模擬工況見(jiàn)表6。

        表6 模擬的三段進(jìn)水A/O生物接觸氧化工藝進(jìn)水流量分配系數(shù)

        2 結(jié)果分析與討論

        2.1 穩(wěn)態(tài)條件工藝運(yùn)行結(jié)果

        當(dāng)各段的進(jìn)水流量分配系數(shù)均為0.333 3,回流比分別為0.50、0.75和1.00時(shí),三段進(jìn)水A/O生物接觸氧化工藝對(duì)高濃度生活污水具有較好的處理效果(圖3)。回流比分別為0.50、0.75和1.00時(shí),出水總氮分別為18.72 mg·L-1、14.77 mg·L-1和13.58 mg·L-1,對(duì)應(yīng)的總氮去除效率分別為82.84%±0.31%、85.36%±0.49%和86.12%±1.83%。工藝的總氮去除效率隨回流比的增大而增大?;亓鞅鹊拇笮Q定了從第3段回流進(jìn)入第1段缺氧池進(jìn)行反硝化的硝態(tài)氮的量?;亓鞅仍酱?,進(jìn)入第1段缺氧池參與反硝化的硝態(tài)氮越多,工藝的總氮去除效率越高。但是過(guò)高的回流比會(huì)導(dǎo)致生化反應(yīng)有效反應(yīng)時(shí)間減小,反而會(huì)降低工藝的總氮去除效率[18]。此外,回流液中的溶解氧進(jìn)入缺氧池會(huì)破壞缺氧環(huán)境, 使反硝化細(xì)菌活性降低, 抑制了反硝化細(xì)菌的生長(zhǎng)繁殖,導(dǎo)致反硝化作用降低[19],實(shí)際工藝運(yùn)行過(guò)程中的回流比不宜過(guò)大。進(jìn)水進(jìn)入各段缺氧池會(huì)使缺氧池出水的COD和氨氮濃度升高,在經(jīng)過(guò)后續(xù)的好氧池降解和氧化后,工藝出水的COD和氨氮濃度維持在較低水平。

        1A、2A、3A分別表示缺氧池的3個(gè)反應(yīng)段;1O、2O、3O分別表示好氧池的3個(gè)反應(yīng)段

        當(dāng)三段進(jìn)水A/O生物接觸氧化工藝的回流比為0.75~1.00時(shí),工藝的總氮去除效率較高,大于85%。與其他進(jìn)水碳氮比為10的多段進(jìn)水A/O生物脫氮工藝[3]相比,本研究的總氮去除效率略低。其原因有以下兩點(diǎn):①本研究采用的進(jìn)水流量和容積分配并非最優(yōu)。最佳的流量分配系數(shù)條件下,進(jìn)入各段缺氧池的COD應(yīng)恰好將上一段好氧池產(chǎn)生的硝態(tài)氮完全反硝化,從而使得工藝可以達(dá)到最大的總氮去除效率。本研究進(jìn)水的碳氮比為10,碳源充足,工藝前段的進(jìn)水流量應(yīng)高于后段的進(jìn)水流量,工藝才可達(dá)到較高的總氮去除效率[20]。②從實(shí)驗(yàn)結(jié)果來(lái)看,由于進(jìn)水流量分配并非最優(yōu),缺氧池出水的總氮高于氨氮濃度,硝態(tài)氮去除不徹底,會(huì)累積至出水導(dǎo)致工藝總氮去除效率降低。通過(guò)對(duì)進(jìn)水流量和反應(yīng)段容積進(jìn)行優(yōu)化分配后,多段進(jìn)水A/O生物脫氮工藝的總氮去除效率可得到明顯提升[15]。

        2.2 模型的驗(yàn)證和校驗(yàn)

        采用回流比為0.5時(shí)的試驗(yàn)結(jié)果對(duì)建立模型進(jìn)行驗(yàn)證表明,模型取用的參數(shù)數(shù)值合適,模型的預(yù)測(cè)值和試驗(yàn)實(shí)測(cè)值吻合度較高,僅A段COD的預(yù)測(cè)值和模擬值出現(xiàn)了一定的偏差(圖4)。缺氧池屬于低溶解氧環(huán)境,當(dāng)有COD和硝態(tài)氮存在時(shí)會(huì)發(fā)生反硝化反應(yīng)。生物膜結(jié)構(gòu)獨(dú)特,外層形成缺氧環(huán)境,由于溶解氧無(wú)法擴(kuò)散至生物膜深層,生物膜內(nèi)部會(huì)形成厭氧環(huán)境[12]。COD可能會(huì)通過(guò)厭氧發(fā)酵途徑降解,而本研究建立的模型中并不包括COD的厭氧降解過(guò)程,因此,缺氧池的COD模擬值與實(shí)測(cè)值之間存在一定偏差。

        采用回流比為1.0時(shí)的試驗(yàn)數(shù)據(jù)對(duì)本研究建立的模型進(jìn)行校驗(yàn)表明,模型能較好地模擬三段進(jìn)水A/O生物接觸氧化工藝的運(yùn)行結(jié)果,可以較好地模擬各單元出水的氨氮、硝態(tài)氮濃度和COD(圖4)。在校驗(yàn)工況下,除第1段好氧池總氮外,模型預(yù)測(cè)的COD、氨氮濃度和總氮與實(shí)測(cè)值吻合度較高。生物膜的獨(dú)特結(jié)構(gòu)使生物膜內(nèi)部能夠在好氧條件下維持缺氧微環(huán)境,在此情況下生物膜內(nèi)部會(huì)發(fā)生同步硝化反硝化過(guò)程[21],使得試驗(yàn)得到的好氧池出水硝態(tài)氮濃度低于模擬值。本研究建立的模型并未考慮同步硝化反硝化過(guò)程,導(dǎo)致總氮的模擬結(jié)果和實(shí)測(cè)結(jié)果產(chǎn)生了一定偏差。此外,工藝出水COD實(shí)測(cè)值略高于模擬預(yù)測(cè)值。為簡(jiǎn)化模型,本研究建立的模型中有機(jī)物組分僅包含SS和XS,對(duì)于生化反應(yīng)過(guò)程沒(méi)有影響的SI和XI組分并未包含在模型中,然而在實(shí)際的出水中含不可生物降解的有機(jī)物[13],導(dǎo)致實(shí)際工藝出水COD比模型預(yù)測(cè)值略高。

        1A、2A、3A分別表示缺氧池的3個(gè)反應(yīng)段;1O、2O、3O分別表示好氧池的3個(gè)反應(yīng)段

        2.3 進(jìn)水流量分配系數(shù)的優(yōu)化

        采用校驗(yàn)后的模型對(duì)不同流量分配系數(shù)條件下三段進(jìn)水A/O生物接觸氧化工藝的模擬表明,當(dāng)進(jìn)水碳氮比為10,回流比為1,r1∶r2∶r3=0.5∶0.3∶0.2時(shí),工藝的總氮去除效率高于其他工況(圖5)。隨著第1段進(jìn)水流量分配系數(shù)的增加和第3段進(jìn)水流量分配系數(shù)的減小,雖然第1段缺氧池出水的COD、氨氮濃度和總氮增加,但最終工藝的出水總氮下降,工藝的脫氮效率提升。然而,進(jìn)入第1段的進(jìn)水流量不能過(guò)大,否則會(huì)使得進(jìn)入第2段和第3段的進(jìn)水中含有的COD不足以將上一段好氧池產(chǎn)生的硝態(tài)氮完全反硝化,最終積累至出水,導(dǎo)致工藝總氮去除效率降低,如表6中的工況1。同樣地,當(dāng)過(guò)量的污水被分配至第3段時(shí),由于好氧池產(chǎn)生的硝態(tài)氮僅有部分能回流至第1段缺氧池進(jìn)行部分去除,工藝的總氮去除效率也會(huì)下降。

        1A、2A、3A分別表示缺氧池的3個(gè)反應(yīng)段;1O、2O、3O分別表示好氧池的3個(gè)反應(yīng)段

        3 結(jié) 語(yǔ)

        (1)多段進(jìn)水A/O生物接觸氧化工藝對(duì)高濃度生活污水的COD、氨氮和總氮有很好的處理效果。工藝的內(nèi)回流比控制在0.75~1.00時(shí),工藝的總氮去除效率大于85%。

        (2)當(dāng)進(jìn)水碳氮比為10,回流比為1,三段進(jìn)水A/O生物接觸氧化工藝的流量分配系數(shù)分別為0.5、0.3和0.2時(shí),工藝可以取得最大的總氮去除效率。

        (3)基于活性污泥1號(hào)數(shù)學(xué)模型(ASM1),結(jié)合Fick擴(kuò)散定律建立的多段進(jìn)水A/O生物膜工藝數(shù)學(xué)模型能夠用于模擬多段進(jìn)水A/O生物接觸氧化工藝的運(yùn)行和對(duì)工藝進(jìn)水流量分配系數(shù)的優(yōu)化。

        猜你喜歡
        總氮氧池水流量
        M701F4燃?xì)廨啓C(jī)TCA系統(tǒng)冷卻水流量異常分析
        青海電力(2022年1期)2022-03-18 12:08:36
        焦化廢水硝化系統(tǒng)受沖擊后恢復(fù)系統(tǒng)穩(wěn)定參考實(shí)例
        蒙大公司污水裝置好氧池曝氣系統(tǒng)改造后節(jié)能降耗見(jiàn)成效
        國(guó)產(chǎn)納濾膜脫除硝酸鹽的試驗(yàn)研究
        基于重力方向影響的低壓渦輪葉片水流量測(cè)量數(shù)值計(jì)算
        甲醇污水站QWSTN法北池系統(tǒng)COD高的原因及應(yīng)對(duì)措施
        中氮肥(2019年2期)2019-04-22 10:46:02
        Enumeration and molecular characterization of circulating tumor cell using an in vivo capture system in squamous cell carcinoma of head and neck
        氣相分子吸收光譜法在線氧化消解測(cè)定水質(zhì)總氮研究
        SBBR進(jìn)水NH+4—N濃度對(duì)反硝化及總氮的影響
        綠色科技(2016年20期)2016-12-27 16:41:16
        松花湖水體中總磷總氮濃度變化特征的研究
        疯狂撞击丝袜人妻| 日韩在线一区二区三区中文字幕| 亚洲熟妇自偷自拍另类| 国精品午夜福利视频不卡| 乱码午夜-极国产极内射| 五月天无码| 亚洲综合网在线观看首页| 久天啪天天久久99久孕妇| 精品极品视频在线观看| 日本精品少妇一区二区三区| 久久精品人人做人人爽| 99热国产在线| 亚洲午夜久久久精品国产| 日韩女优图播一区二区| 青娱乐极品视觉盛宴国产视频 | 美腿丝袜一区二区三区| 国产一区二区三区在线大屁股| 国产免费一区二区三区免费视频| 亚洲欧美日韩在线观看一区二区三区| 狠狠久久精品中文字幕无码| 久久精品一区二区三区不卡牛牛| 亚洲人不卡另类日韩精品| 免费观看的a级毛片的网站| 成人xx免费无码| 久久亚洲国产高清av一级 | 成年女人免费v片| 粗大挺进孕妇人妻在线| 国产精品无码制服丝袜| 免费观看羞羞视频网站| 国产黑色丝袜在线观看下| 日韩狼人精品在线观看| av在线免费观看网站免费| 水蜜桃精品一二三| 亚洲欧美日韩中文字幕网址 | 亚洲熟妇AV一区二区三区宅男| 一区二区三区视频偷拍| 中文字幕免费在线观看动作大片| 比比资源先锋影音网| 女同成片av免费观看| 一道本久久综合久久鬼色| 白天躁晚上躁麻豆视频|