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

        ?

        基于可調(diào)諧二極管激光吸收光譜技術的爐內(nèi)燃燒場參數(shù)在線監(jiān)測系統(tǒng)設計

        2022-05-11 05:21:52李源郭志成孟曉超陳科峰任利明毛睿岑可法
        發(fā)電技術 2022年2期
        關鍵詞:譜線鍋爐網(wǎng)格

        李源,郭志成,孟曉超,陳科峰,任利明,毛睿,岑可法

        (1.潤電能源科學技術有限公司,河南省 鄭州市 450052;2.能源清潔利用國家重點實驗室(浙江大學),浙江省 杭州市 310027)

        0 引言

        近年來,在煤炭價格高漲和環(huán)保日趨嚴格的背景下,電力企業(yè)為了提高自身競爭力,基于融合數(shù)字化技術和工業(yè)生產(chǎn)過程提出了智能電廠概念[1-2]。智能電廠一般分為4 個層次,從低至高依次為智能設備層、智能控制層、智能生產(chǎn)監(jiān)管層和智能管理層[3-5]。其中,智能設備層是智能電廠最為基礎的一層,智能設備實現(xiàn)了電廠運行過程中參數(shù)的測量。

        鍋爐是火力發(fā)電廠中的能量輸出設備,爐內(nèi)涉及復雜的多相湍流燃燒過程。爐內(nèi)燃燒場參數(shù)的測量對研究燃燒機理、提高燃燒效率和控制污染物排放等具有積極意義?;鹧鏈囟群徒M分濃度是爐內(nèi)最基礎的燃燒參數(shù),燃燒參數(shù)的測量技術主要涉及聲音、電磁、圖像、激光和光譜等技術[6-11]。

        燃燒溫度是表征燃燒強度和熱釋放率最直接的物理量;CO 和CO2是表征燃燒反應的關鍵指標,且CO 含量與壁面高溫腐蝕密切相關[12];NO是燃燒過程中主要的污染物,且NO 濃度在一定程度上也能反映出燃燒狀況。常規(guī)的熱電偶和煙氣分析儀等測量手段均為單點測量,但單點測量只能反映爐內(nèi)整體的燃燒狀況,而無法反映爐內(nèi)局部的燃燒狀況。單點測量很難監(jiān)測爐內(nèi)火焰中心是否偏離,無法有效監(jiān)測和治理鍋爐效率低、污染物排放濃度高、水冷壁高溫腐蝕和結渣積灰等問題。此外,常規(guī)的熱電偶和煙氣分析儀等接觸式測量手段具有設備壽命短、干擾流場、響應速度慢和靈敏度低等缺點,且這些設備無法長期適應高溫、多塵、多變的燃燒環(huán)境。針對目前鍋爐快速變負荷和工況多變的現(xiàn)狀,為探尋鍋爐在復雜工況時的最佳運行策略,迫切需要對爐內(nèi)三維燃燒溫度和多種氣體組分濃度進行在線測量。

        激光燃燒診斷技術是一種在線、非接觸式測量手段,可以實現(xiàn)燃燒參數(shù)的快速無干擾測量[13]??烧{(diào)諧二極管激光吸收光譜(tunable diode laser absorption spectroscopy,TDLAS)技術具有系統(tǒng)結構簡單、精度高、響應快和適應環(huán)境能力強等優(yōu)點[14]。傳統(tǒng)的TDLAS 技術是通過對待測區(qū)域內(nèi)某一組分沿光學路徑的積分吸收度進行測量,從而反演出光學路徑上的平均溫度和組分濃度[15]。有學者[16-17]提出將TDLAS 技術與計算機層析技術(computer tomography,CT)結合,用于重建火焰內(nèi)部的二維火焰溫度和組分濃度場。目前,基于TDLAS 技術測量溫度和組分濃度主要還局限于實驗室水平,尚未有在實際鍋爐上應用的報道。

        鍋爐燃燒優(yōu)化的目的是實現(xiàn)高效燃燒、低污染物排放及減小熱損失[18-21]。由于實際爐膛為三維結構,爐內(nèi)火焰溫度和組分濃度分布差異較大。但當前的優(yōu)化控制策略還局限于整體的燃燒狀態(tài),沒有對爐內(nèi)燃燒區(qū)域進行劃分,這樣的控制策略不能針對鍋爐進行精細化控制。為實現(xiàn)鍋爐的精細化優(yōu)化控制,有必要在優(yōu)化策略中引入燃燒場參數(shù),包括溫度場和多組分濃度場。此外,為減弱高溫腐蝕和結渣積灰,針對爐內(nèi)三維溫度和多組分濃度進行實時監(jiān)測也是十分必要的。本文基于TDLAS技術設計了一種爐內(nèi)三維溫度和多組分濃度在線監(jiān)測系統(tǒng)。

        1 基本原理

        1.1 TDLAS原理

        TDLAS 技術的理論基礎是Beer-Lambert 定律,其原理如圖1所示。

        圖1 Beer-Lambert原理Fig.1 Beer-Lambert principle

        假設一束頻率為v的激光束穿過待測區(qū)域,被均勻氣體介質(zhì)吸收,光路上的積分吸光度為

        式中:p為混合氣體總壓;X為介質(zhì)氣體組分濃度;L為激光束在介質(zhì)氣體中的傳輸距離;S為吸收線強度;T為溫度;av為氣體吸收系數(shù)。在熱平衡體系的條件下,譜線強度S隨溫度T的變化可以表示為

        式中:S(T0)為參考溫度T0時的譜線強度;E''為吸收線的低態(tài)能級值;h為Planck 常數(shù);k為Boltzmann 常數(shù);c為光速;v0為吸收譜線中心頻率;Q(T)為溫度T時的配分函數(shù)值,可用多項式表示為

        其中z0,z1,z2,z3為擬合多項式系數(shù)。

        在TDLAS雙線測溫法[22]中,挑選同種氣體分子的2條吸收譜線,將2個吸收線強值相除,消掉配分函數(shù),得到不同吸收譜線強度的比值R為

        則待測區(qū)域溫度T通過式(2)和式(4)可表示為

        獲取溫度T后,選擇吸收較大的譜線來計算介質(zhì)氣體組分濃度,可得

        為了測量煙氣溫度及CO、CO2和NO 濃度,在激光波段范圍內(nèi)選擇多條相對獨立的吸收譜線。選擇位于7 153.75 cm-1和7 154.35 cm-1處的2 條H2O 吸收譜線,利用雙線測溫法反演燃燒溫度;選擇位于6 338.59 cm-1處的CO 吸收譜線和6 337.99 cm-1處的CO2吸收譜線,同步反演CO 和CO2組分濃度;選擇位于1 900.07 cm-1處的NO 吸收譜線,反演NO 組分濃度。吸收譜線在特定條件下的吸光度可利用HITRAN數(shù)據(jù)庫模擬計算[22]。

        1.2 代數(shù)迭代算法

        代數(shù)迭代算法(algebraic reconstruction technique,ART)算法是一種典型的CT算法[23]。對待測區(qū)域進行二維重建時,將待測區(qū)域劃分成n×n個網(wǎng)格單元,并假定網(wǎng)格單元內(nèi)待測量分布均勻。圖2為將待測區(qū)域劃分成10×10的網(wǎng)格單元。

        圖2 網(wǎng)格化待測區(qū)域Fig.2 Girded measurement region

        單光路的積分吸收率可表示為光路經(jīng)過的所有網(wǎng)格對吸收的貢獻之和:

        式中:Lj為光路經(jīng)過第j個網(wǎng)格的光程;J為網(wǎng)格總數(shù);aij為第j個網(wǎng)格內(nèi)第i條射線的積分吸收系數(shù)。

        基于ART算法,網(wǎng)格單元的吸收系數(shù)的迭代計算式為

        式中:Ai為第i條射線的積分吸收率;β(k)為第k次迭代的松弛系數(shù);Lij為第i條射線穿過第j個網(wǎng)格的光程。

        利用修正型ART 算法[24]獲取譜線v1和v2的吸收系數(shù)aν1,i和aν2,i,并計算各網(wǎng)格內(nèi)吸收譜線強度的比值Ri:

        將式(9)代入式(5)便可以得到對應網(wǎng)格內(nèi)的溫度值Ti:

        再由式(2)得出相應吸收譜線強度Sv1,i(Ti)和Sv2,i(Ti),可計算得到對應網(wǎng)格內(nèi)氣體的組分濃度:

        2 系統(tǒng)設計

        2.1 一維燃燒參數(shù)監(jiān)測系統(tǒng)方案

        為了實現(xiàn)爐內(nèi)煙氣溫度和多組分濃度的在線監(jiān)測,設計一種基于TDLAS的單光纖爐內(nèi)燃燒參數(shù)監(jiān)測系統(tǒng)。圖3為一維燃燒參數(shù)監(jiān)測系統(tǒng)框圖,該系統(tǒng)基于TDLAS技術和光開關技術實現(xiàn)高溫煙氣溫度和多組分濃度的同時在線測量。該系統(tǒng)包括上位機、信號發(fā)生器、光開關、光纖準直器、燃燒場、凸透鏡、光電轉換器、鎖相放大器、采集卡和單模光纖,以及激光控制器A、B、C,分布式反饋(distributed feedback laser,DFB)激光器A、B、C。其中,光開關可以對傳輸激光進行物理切換或邏輯操作,實現(xiàn)多參數(shù)的同時測量,簡化系統(tǒng)結構,節(jié)約設備成本。

        圖3 一維燃燒參數(shù)監(jiān)測系統(tǒng)Fig.3 One-dimensional combustion parameter monitoring system

        根據(jù)待測燃燒場參數(shù)選擇DFB激光器的工作波段,DFB 激光器工作波段要涵蓋待測參數(shù)的吸收譜線范圍。DFB 激光器A 選擇中心波長為1 397.8 nm,一次掃描可獲取7 153.75 cm-1和7 154.35 cm-1處的2條H2O吸收譜線,利用雙線測溫法反演燃燒溫度;DFB 激光器B 選擇中心波長為1 579 nm,一次掃描可獲取6 338.59 cm-1處的CO 吸收譜線和6 337.99 cm-1處的CO2吸收譜線,同步反演CO 和CO2組分濃度;DFB 激光器C 選擇中心波長為5 263 nm,一次掃描可獲取1 900.07 cm-1處的NO 吸收譜線,反演NO 組分濃度。該監(jiān)測系統(tǒng)包括如下步驟:

        1)開啟儀器設備并調(diào)整光路;

        2)上位機控制信號發(fā)生器輸出電信號,3 個激光控制器接收電信號后分別驅動對應的DFB激光器發(fā)射激光信號;

        3)3個DFB激光器發(fā)射的激光信號通過光纖傳輸至光開關,光開關將3路激光信號耦合成1路激光信號輸出,由上位機分時復用控制光開關中三路激光信號的交替輸出;

        4)光纖準直器將光開關輸出的激光信號校正為一束平行激光;

        5)平行激光束穿過燃燒場,部分光譜被氣體吸收,剩余激光束發(fā)生微小偏折并由凸透鏡匯聚至光電轉換器;

        6)光電轉換器將入射光信號轉換為電信號,傳輸至鎖相放大器,鎖相放大器具有2 路輸入信號和1 路輸出信號,輸入信號包括光電轉換器輸出的電信號和信號發(fā)生器輸出的參考信號,輸入電信號經(jīng)調(diào)試后產(chǎn)生二次諧波輸出信號;

        7)采集卡接收二次諧波信號,并將模擬信號轉變?yōu)閿?shù)字信號傳輸至上位機;

        8)上位機對數(shù)字信號進行采集和處理,反演計算出燃燒溫度和多組分氣體濃度。

        該系統(tǒng)與火焰圖像系統(tǒng)[6]、聲學測溫系統(tǒng)[8]、常規(guī)火焰光譜系統(tǒng)[10-11]等相比,具有以下優(yōu)點:1)基于TDLAS技術實現(xiàn)了高溫煙氣溫度和多組分濃度的在線測量;2)基于光開關技術實現(xiàn)了多個燃燒參數(shù)的同時測量,簡化了系統(tǒng)結構,系統(tǒng)中僅需1臺光纖準直器、1臺光電轉換器和1臺鎖相放大器,節(jié)約了設備成本費用;3)光開關技術為分時測量技術,消除了多組分氣體測量時氣體間的相互干擾。

        2.2 多維燃燒參數(shù)監(jiān)測系統(tǒng)

        為了實現(xiàn)爐內(nèi)多維燃燒參數(shù)的在線監(jiān)測,設計一種基于TDLAS的爐內(nèi)三維燃燒參數(shù)監(jiān)測系統(tǒng)。三維燃燒參數(shù)測量包括3個步驟:1)將爐膛劃分為多層二維水平切片;2)基于ART算法重建二維燃燒參數(shù);3)利用二維燃燒參數(shù)數(shù)據(jù)重構三維燃燒參數(shù)。圖4為一種基于TDLAS的爐內(nèi)三維燃燒參數(shù)監(jiān)測系統(tǒng),包括TDLAS控制與處理單元、激光發(fā)射單元、測試單元和信號采集單元。基于分層切片法將測試單元中的爐膛劃分為N層切片,將三維燃燒參數(shù)測量轉化為N個二維燃燒參數(shù)測量。

        圖4 三維燃燒參數(shù)監(jiān)測系統(tǒng)Fig.4 Three-dimensional combustion parameter monitoring system

        將每一層切片分為n×n個網(wǎng)格單元,利用TDLAS技術和ART算法重建該層切片的二維燃燒參數(shù)。圖5為爐內(nèi)二維切片燃燒參數(shù)監(jiān)測系統(tǒng),包括上位機、信號發(fā)生器、光開關、單模光纖、分光器、光纖準直器、二維燃燒場、光電轉換器、采集卡、鎖相放大器,以及激光控制器A、B、C,DFB激光器A、B、C。其中,分光器可將單束激光均勻分成多束激光,既保證了多束激光的同步性又簡化了系統(tǒng)結構。與2.1節(jié)相同,DFB激光器的工作波段根據(jù)待測燃燒場參數(shù)進行選擇,DFB 激光器A 選擇中心波長為1 397.8 nm,DFB激光器B選擇中心波長為1 579 nm,DFB激光器C選擇中心波長為5 263 nm。

        圖5 二維燃燒參數(shù)監(jiān)測系統(tǒng)Fig.5 Two-dimensional combustion parameter monitoring system

        爐內(nèi)二維切片燃燒參數(shù)在線監(jiān)測系統(tǒng)包括如下步驟:

        1)打開所述儀器設備并調(diào)整光路;

        2)上位機控制信號發(fā)生器輸出模擬信號,3 個激光控制器接收模擬信號后分別驅動對應的DFB激光器發(fā)射激光信號;

        3)3個DFB激光器發(fā)射的激光信號通過光纖傳輸至光開關,光開關將3路激光信號耦合成1路激光信號輸出,由上位機分時復用控制光開關中3路激光信號的交替輸出;

        4)光開關輸出的激光信號通過單模光纖傳輸至分光器,分光器將輸入的1 路激光信號分為2n路激光信號;

        5)分光器輸出的2n路激光信號通過光纖分別傳輸至2n個光纖準直器,每個光纖準直器將對應的1路激光信號校正為1束平行激光;

        6)2n束平行激光穿過二維燃燒場,形成一個n×n的光路網(wǎng)格,2n束激光的部分光譜被氣體吸收后分別被對應的2n個光電轉換器接收;

        7)2n個光電轉換器分別將接收的激光信號轉換為電信號,并統(tǒng)一傳輸至采集卡;

        8)鎖相放大器具有2 路輸入信號和1 路輸出信號,輸入信號包括采集卡輸出的電信號和信號發(fā)生器輸出的參考信號,輸入電信號經(jīng)調(diào)試后產(chǎn)生二次諧波輸出信號,二次諧波信號傳輸至上位機;

        9)上位機對輸入信號進行采集和處理,計算反演出二維火焰溫度和組分濃度場。

        在獲取N層切片的二維火焰溫度和組分濃度場數(shù)據(jù)后,基于分層切片法重構出爐內(nèi)三維火焰溫度和組分濃度場。

        該系統(tǒng)與三維火焰圖像系統(tǒng)[7]、三維聲學測溫系統(tǒng)[9]等相比,具有以下優(yōu)點:1)基于TDLAS 技術和ART 算法實現(xiàn)了爐內(nèi)二維火焰溫度和組分濃度場的非接觸式在線測量,并在此基礎上基于分層切片法重構三維火焰溫度和組分濃度場;2)光開關技術作為一種分時測量技術實現(xiàn)了多個燃燒參數(shù)的同時測量,并消除了多組分氣體之間的相互干擾;3)分光器技術簡化了二維燃燒場參數(shù)測量系統(tǒng)的結構,節(jié)約了設備成本。

        3 結果與討論

        目前,由于現(xiàn)場條件及設備可靠性等限制,本系統(tǒng)尚未布置于實際電站鍋爐。為展示爐內(nèi)三維燃燒參數(shù)的測量效果,本文基于鍋爐燃燒仿真數(shù)據(jù)重建爐內(nèi)三維燃燒參數(shù)分布,并以某300 MW直流燃煤鍋爐為對象進行燃燒仿真研究。鍋爐型號為SG-1025/16.7-M313UP,采用中儲式制粉系統(tǒng)。燃燒方式為四角切圓燃燒,設有5層一次風、7層二次風、2層三次風,其中,二次風又分為分離燃盡風(separated overfire air,SOFA)、緊靠布置 燃 盡 風(closed-coupled overfire air,CCOFA)、周界風及其他二次風。表1為鍋爐燃用煤質(zhì)特性。表2 為鍋爐燃燒仿真邊界條件。利用ICEM 軟件進行全爐膛網(wǎng)格劃分。利用Fluent 軟件進行鍋爐燃燒數(shù)值計算,其中,湍流流動采用k-ε模型,輻射傳熱采用P1 模型,顆粒運動采用隨機軌道模型,氣相湍流燃燒采用混合分數(shù)/概率密度函數(shù)模型,煤粉揮發(fā)分析出采用單步反應模型,焦炭燃燒計算采用動力/擴散控制反應速率模型。

        表1 鍋爐燃用煤質(zhì)特性Tab.1 Properties of the burning coal

        表2 鍋爐燃燒仿真邊界條件Tab.2 Combustion simulation boundaries of the boiler

        根據(jù)鍋爐結構和燃燒仿真結果,將鍋爐劃分為7 層切片,每層切片劃分為10×10 的網(wǎng)格,如圖6 所示。根據(jù)鍋爐燃燒仿真結果,獲取每層切片的二維燃燒數(shù)據(jù)后,基于分層切片原理,重構獲取爐內(nèi)三維燃燒參數(shù)。

        圖6 鍋爐多層切片F(xiàn)ig.6 Boiler multilayer sections

        圖7 為基于仿真數(shù)據(jù)重建的爐內(nèi)三維溫度及CO、CO2、NO 濃度分布。通過三維燃燒場參數(shù)分布可以清晰地了解爐內(nèi)各區(qū)域的燃燒信息。從圖7(a)可以看出,在水平方向上,爐內(nèi)溫度表現(xiàn)為中心高、四周低;在高度方向上,爐內(nèi)溫度隨著標高的增加表現(xiàn)為先增加再減小的趨勢。爐內(nèi)火焰溫度的分布主要與火焰中心位置相關,火焰中心為爐內(nèi)燃燒最劇烈也是火焰溫度最高的地方。從圖7(b)可以看出,主燃燒器區(qū)域CO 濃度最大,這主要與鍋爐低氮運行時此處嚴重缺氧有關。從圖7(c)可以看出,爐底部和上部區(qū)域CO2濃度較高。爐底部CO2濃度較高主要是由于爐底部空氣量較少,大量反應生成的CO2積聚于爐底部。爐上部CO2濃度較高則主要是由于該區(qū)域煤粉已與空氣完全反應。從圖7(d)可以看出,主燃燒區(qū)域存在一個低NO 濃度區(qū)域,這主要是由于主燃燒區(qū)域嚴重缺氧,導致該區(qū)域生成大量CO、NHi及HCN 等產(chǎn)物[25],將形成的NO 還原成N2。

        圖7 爐內(nèi)三維燃燒參數(shù)重建Fig.7 Remonstration of three-dimensional combustion parameters in furnace

        利用爐內(nèi)三維溫度分布特征可以準確判斷爐內(nèi)火焰中心是否偏斜、火焰高度是否合理、熱負荷分布是否均勻等情況,有利于提高鍋爐效率、防治結渣積灰、消除汽溫/煙溫偏差。利用爐內(nèi)CO和CO2三維分布特征可以準確判斷爐內(nèi)各區(qū)域燃燒反應狀態(tài),有助于鍋爐精準配風/配煤、消除燃燒死區(qū)、提升燃燒效率。此外,由于高溫腐蝕與CO 分布存在一定的相關性[12],基于爐內(nèi)CO 三維分布情況可以有效防治水冷壁高溫腐蝕。利用爐內(nèi)NO 三維分布特征可以準確判斷爐內(nèi)各區(qū)域NO生成特征,有助于制定合理的配風、配煤、摻燒方式,以降低鍋爐出口NO 排放量。因此,爐內(nèi)三維燃燒參數(shù)在線監(jiān)測十分有利于鍋爐精細化燃燒調(diào)整。

        4 結論

        基于TDLAS技術設計了一種爐內(nèi)燃燒場參數(shù)在線監(jiān)測系統(tǒng),包括一維和多維燃燒參數(shù)在線監(jiān)測系統(tǒng),可監(jiān)測的燃燒參數(shù)包括溫度及CO、CO2和NO濃度等。得到以下結論:

        1)在一維燃燒參數(shù)在線監(jiān)測系統(tǒng)方案中,基于TDLAS技術和光開關技術實現(xiàn)爐內(nèi)多個燃燒參數(shù)的同時在線測量。其中,光開關技術主要用于同時測量多個燃燒參數(shù)。

        2)在二維燃燒參數(shù)在線監(jiān)測系統(tǒng)方案中,基于TDLAS 技術、分光器技術和ART 算法實現(xiàn)平面內(nèi)多個燃燒參數(shù)的同時在線測量。其中,分光器技術主要用于測量二維燃燒場參數(shù)?;诜謱忧衅韺⑷S燃燒參數(shù)測量問題轉化為多個二維燃燒參數(shù)測量問題。

        3)在三維燃燒參數(shù)在線監(jiān)測系統(tǒng)方案中,先將三維爐膛劃分為多個二維分層切片,再基于TDLAS技術和ART算法重建分層切片的二維燃燒參數(shù),最后利用多個二維燃燒參數(shù)數(shù)據(jù)重構三維燃燒場參數(shù)。

        猜你喜歡
        譜線鍋爐網(wǎng)格
        用全等三角形破解網(wǎng)格題
        基于HITRAN光譜數(shù)據(jù)庫的合并譜線測溫仿真研究
        對干熄焦余熱鍋爐運行爆管的幾點探討
        昆鋼科技(2020年6期)2020-03-29 06:39:50
        反射的橢圓隨機偏微分方程的網(wǎng)格逼近
        12CrlMoV鍋爐吊桿用鋼的開發(fā)生產(chǎn)實踐
        山東冶金(2018年5期)2018-11-22 05:12:06
        鐵合金光譜譜線分離實驗研究
        電子測試(2018年11期)2018-06-26 05:56:00
        重疊網(wǎng)格裝配中的一種改進ADT搜索方法
        關于鍋爐檢驗的探討
        鍶原子光鐘鐘躍遷譜線探測中的程序控制
        基于曲面展開的自由曲面網(wǎng)格劃分
        一区二区三区中文字幕在线播放| 国产精品亚洲A∨天堂| 天堂Av无码Av一区二区三区| 日本午夜a级理论片在线播放| 欧美精品色婷婷五月综合| 日韩激情无码免费毛片 | 亚洲中文字幕亚洲中文| 国产精品主播在线一区二区 | 无码精品人妻一区二区三区98| 在线视频观看一区二区| 国产免费无遮挡吸奶头视频| 国产suv精品一区二区| 啊v在线视频| 蜜桃精品视频一二三区| 中国老熟妇506070| 无遮无挡三级动态图| 精品视频在线观看一区二区有 | 2020亚洲国产| 97成人精品在线视频| 99在线精品视频在线观看| 少妇被粗大的猛进69视频| 国产一区二区三区免费精品| 中文字幕人妻互换av| 国产无吗一区二区三区在线欢| 成人片黄网站色大片免费观看app| 无码人妻丰满熟妇区免费| 日韩人妻久久中文字幕| 伊人久久大香线蕉av不卡| 欧美人妻日韩精品| 91中文字幕精品一区二区| 97se色综合一区二区二区| 越南女子杂交内射bbwxz| 天天澡天天揉揉AV无码人妻斩 | 午夜三级网| 国产av大片久久中文字幕| 人妻少妇被粗大爽.9797pw| 精品久久亚洲中文无码| 国产成人AⅤ| 亚洲色图三级在线观看| 在线综合亚洲欧洲综合网站| 亚洲V无码一区二区三区四区观看|