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

        ?

        火災(zāi)工況下隧道通風(fēng)系統(tǒng)控制策略編制技術(shù)研究

        2024-03-06 00:00:00代東林
        交通科技與管理 2024年7期

        摘要 文章圍繞公路隧道通風(fēng)系統(tǒng)控制策略編制所需的斷面風(fēng)速測(cè)試和臨界風(fēng)速計(jì)算兩方面關(guān)鍵技術(shù),詳細(xì)介紹其基本原理和現(xiàn)有技術(shù)方案,在此基礎(chǔ)上提出一種多斷面風(fēng)速監(jiān)測(cè)新思路并深入分析風(fēng)速傳感器選型及其無(wú)線(xiàn)組網(wǎng)方式。此外,提出一種火災(zāi)工況下臨界風(fēng)速快速預(yù)測(cè)新思路,重點(diǎn)闡述基于FDS模擬試驗(yàn)的臨界風(fēng)速樣本獲取要點(diǎn),以及基于機(jī)器學(xué)習(xí)算法的多因素,即臨界風(fēng)速預(yù)測(cè)模型構(gòu)建要點(diǎn)。

        關(guān)鍵詞 隧道通風(fēng);火災(zāi)工況;風(fēng)速監(jiān)測(cè);臨界風(fēng)速;FDS

        中圖分類(lèi)號(hào) U453.5 文獻(xiàn)標(biāo)識(shí)碼 A 文章編號(hào) 2096-8949(2024)07-0001-04

        0 引言

        近年來(lái),公路隧道建設(shè)蓬勃發(fā)展,隧道數(shù)量和里程不斷增加,單座隧道長(zhǎng)度不斷突破,在給公路交通創(chuàng)造有利條件的同時(shí),隧道內(nèi)各類(lèi)事故隱患也不斷增多,其中隧道火災(zāi)事故更是時(shí)有發(fā)生。隧道內(nèi)一旦發(fā)生火災(zāi),隧道內(nèi)的溫度和煙氣對(duì)人員安全產(chǎn)生嚴(yán)重威脅,稍有不慎將導(dǎo)致嚴(yán)重的人員傷亡和財(cái)產(chǎn)損失。

        2020—2021年全國(guó)公路隧道養(yǎng)護(hù)運(yùn)行管理現(xiàn)狀調(diào)研期間,關(guān)于通風(fēng)系統(tǒng)應(yīng)急控制策略的調(diào)研,共收到495份問(wèn)卷,其中245家管理單位明確表示沒(méi)有制定突發(fā)事件的通風(fēng)控制策略,占比接近50%。在隧道火災(zāi)突發(fā)時(shí),通風(fēng)系統(tǒng)不能盡快發(fā)揮控?zé)?、排煙作用,有時(shí)甚至?xí)鸱醋饔谩?/p>

        事先通過(guò)現(xiàn)場(chǎng)試驗(yàn)探明隧道通風(fēng)資源特征及利用價(jià)值,是編制通風(fēng)系統(tǒng)控制策略的必要前提條件。對(duì)于典型的公路隧道,分別在單組風(fēng)機(jī)獨(dú)立運(yùn)行、多組風(fēng)機(jī)不同組合運(yùn)行條件下,分析得出:斷面風(fēng)速隨時(shí)間及縱向位置變化曲線(xiàn),是揭示通風(fēng)資源特征的一種重要途徑。因此,通風(fēng)系統(tǒng)控制策略編制中,斷面風(fēng)速的可靠監(jiān)測(cè)是通風(fēng)供應(yīng)側(cè)需要解決的關(guān)鍵技術(shù)問(wèn)題。

        從需求側(cè)來(lái)看,預(yù)測(cè)隧道火災(zāi)臨界風(fēng)速需求是編制通風(fēng)系統(tǒng)控制策略的重要客觀(guān)依據(jù)。分別分析隧道海拔高度、火源功率、斷面阻塞比以及不同區(qū)段路面坡度等因素對(duì)臨界風(fēng)速的影響程度,從而提出:綜合考慮多因素影響、普遍適用且及時(shí)性好的臨界風(fēng)速快速預(yù)測(cè)算法,是通風(fēng)需求側(cè)應(yīng)當(dāng)解決的關(guān)鍵技術(shù)問(wèn)題。

        該文針對(duì)通風(fēng)系統(tǒng)控制策略編制所需的斷面風(fēng)速測(cè)試、臨界風(fēng)速計(jì)算兩方面關(guān)鍵技術(shù),詳細(xì)介紹其基本原理和現(xiàn)有方案。在此基礎(chǔ)上,分別提出多斷面風(fēng)速監(jiān)測(cè)和臨界風(fēng)速快速預(yù)測(cè)的新技術(shù)總體思路,并對(duì)其中的傳感器選型與組網(wǎng)、臨界風(fēng)速樣本獲取與預(yù)測(cè)模型構(gòu)建進(jìn)行重點(diǎn)分析。

        1 斷面風(fēng)速測(cè)試原理與多斷面監(jiān)測(cè)新技術(shù)

        1.1 常用測(cè)試方案介紹

        關(guān)于隧道斷面風(fēng)速測(cè)試,多采用傳統(tǒng)風(fēng)速儀在斷面內(nèi)逐點(diǎn)進(jìn)行測(cè)量。交通運(yùn)輸部公路科學(xué)研究院李俊峰等將公路隧道斷面劃分為16個(gè)測(cè)量區(qū)域[1],開(kāi)發(fā)了斷面風(fēng)速多通道檢測(cè)系統(tǒng),從而實(shí)現(xiàn)各區(qū)域中心點(diǎn)風(fēng)速同步測(cè)量;公安部天津消防研究所研發(fā)的隧道防排煙多點(diǎn)風(fēng)速測(cè)量裝置[2],將隧道全斷面細(xì)分成32個(gè)區(qū)域,基于多通道熱線(xiàn)風(fēng)速傳感技術(shù)來(lái)測(cè)量斷面風(fēng)速,且測(cè)量裝置適用于不同尺寸、不同形狀的各類(lèi)隧道。

        這些斷面風(fēng)速測(cè)量方法以隧道全斷面為對(duì)象,近似等面積劃分區(qū)域測(cè)量風(fēng)速,對(duì)于隧道火災(zāi)工況下排煙來(lái)說(shuō),隧道斷面上部或者上中部煙氣集中位置才是通風(fēng)關(guān)鍵區(qū)域,因此有必要在斷面內(nèi)進(jìn)一步劃分關(guān)鍵區(qū)域并監(jiān)測(cè)其平均風(fēng)速。另外,鐵道部科學(xué)研究院進(jìn)行了鐵路隧道斷面平均風(fēng)速位置搜索,獲得了斷面內(nèi)平均風(fēng)速特征點(diǎn)分布曲線(xiàn),對(duì)研究公路隧道斷面關(guān)鍵區(qū)域平均風(fēng)速特征位置分布特點(diǎn)具有一定的參考意義。

        1.2 風(fēng)速傳感器分析

        在隧道斷面風(fēng)速多通道同步測(cè)量方面,風(fēng)速傳感器共有葉片式、風(fēng)杯式、超聲波式和熱敏式四種類(lèi)型,其測(cè)量范圍都能覆蓋1~25 m/s且示值誤差在5%以?xún)?nèi)。從起動(dòng)風(fēng)速大小、現(xiàn)場(chǎng)安裝難易程度以及是否有利于組網(wǎng)出發(fā),分析各類(lèi)風(fēng)速傳感器特點(diǎn)并匯總?cè)绫?所示??傮w來(lái)說(shuō),熱敏式風(fēng)速傳感器對(duì)于隧道斷面風(fēng)速測(cè)量需求來(lái)說(shuō)具有明顯優(yōu)勢(shì)。

        1.3 無(wú)線(xiàn)組網(wǎng)方式分析

        針對(duì)多熱敏風(fēng)速傳感器無(wú)線(xiàn)組網(wǎng),重點(diǎn)考慮數(shù)據(jù)傳輸實(shí)時(shí)性、通信距離長(zhǎng)短和外部服務(wù)依賴(lài)度等因素,常用無(wú)線(xiàn)組網(wǎng)技術(shù)優(yōu)缺點(diǎn)如表2所示。

        LoRa基于Sub-GHz的頻段使其更易以較低功耗遠(yuǎn)距離通信。雖然數(shù)據(jù)速率較低,但網(wǎng)絡(luò)終端容量大、通信距離遠(yuǎn)、信號(hào)穿透力強(qiáng)且功耗低。既可以小范密集部署終端,也可以長(zhǎng)距離分散通信部署終端,與斷面風(fēng)速監(jiān)測(cè)中數(shù)據(jù)通信量低、通信終端多、使用電池供電以及終端既要密集部署又要分散部署等需求最為匹配。

        1.4 多斷面監(jiān)測(cè)新技術(shù)總體思路

        基于以上分析,提出一種隧道風(fēng)速多斷面監(jiān)測(cè)新技術(shù),其總體技術(shù)路線(xiàn)如下:

        (1)分析隧道火災(zāi)工況下煙氣高度對(duì)人員逃生的影響,以此確定斷面通風(fēng)關(guān)鍵區(qū)域。

        (2)在斷面關(guān)鍵區(qū)域使用重型伸縮牛頭架和龍門(mén)架等設(shè)置數(shù)根等間距立柱,每根立柱上等間距布設(shè)若干風(fēng)速計(jì),如圖1所示。

        (3)在各風(fēng)速計(jì)旁設(shè)置LoRa無(wú)線(xiàn)傳輸節(jié)點(diǎn),用于將風(fēng)速數(shù)據(jù)傳輸至LoRa網(wǎng)關(guān),并由LoRa網(wǎng)關(guān)以有線(xiàn)局域網(wǎng)傳輸方式將風(fēng)速數(shù)據(jù)發(fā)送到風(fēng)速監(jiān)測(cè)軟件平臺(tái),如圖2所示。

        (4)在各斷面關(guān)鍵區(qū)域內(nèi)分別確定代表該關(guān)鍵區(qū)域平均風(fēng)速的特征位置。

        (5)對(duì)全隧道多斷面特征位置風(fēng)速同步測(cè)量,擬在隧道內(nèi)每間隔100 m布設(shè)一個(gè)三腳架伸縮桿,伸縮桿頂端布設(shè)風(fēng)速傳感器和LoRa無(wú)線(xiàn)傳輸節(jié)點(diǎn);風(fēng)速數(shù)據(jù)首先通過(guò)LoRa無(wú)線(xiàn)傳輸至附近數(shù)百米遠(yuǎn)的LoRa網(wǎng)關(guān),再經(jīng)無(wú)線(xiàn)網(wǎng)橋連接的局域網(wǎng)傳輸至風(fēng)速監(jiān)測(cè)軟件平臺(tái)。

        2 臨界風(fēng)速計(jì)算原理及快速預(yù)測(cè)新技術(shù)

        2.1 現(xiàn)有計(jì)算原理分析

        關(guān)于隧道火災(zāi)臨界風(fēng)速研究,國(guó)內(nèi)外學(xué)者開(kāi)展了大量工作。臨界風(fēng)速是指隧道內(nèi)發(fā)生火災(zāi)時(shí),能有效地將煙氣控制于火源下風(fēng)方向而不發(fā)生逆流的最小縱向通風(fēng)風(fēng)速,最早由英國(guó)的Thomas于1968年提出。根據(jù)不同火源功率及火源至出口距離對(duì)臨界風(fēng)速的影響研究,結(jié)果顯示:臨界風(fēng)速隨著火源功率的增大而增大,與火源功率的1/3次方成正比,且隨著火源與隧道出口距離的增大呈現(xiàn)出線(xiàn)性增長(zhǎng)趨勢(shì);部分學(xué)者研究了隧道坡度對(duì)射流風(fēng)機(jī)臨界風(fēng)速的影響,結(jié)果顯示上坡時(shí)坡度越大臨界風(fēng)速越小,下坡時(shí)坡度(絕對(duì)值)越大臨界風(fēng)速越大,并對(duì)現(xiàn)有Kennedy臨界風(fēng)速模型中的坡度修正系數(shù)進(jìn)行了優(yōu)化[3];武漢科技大學(xué)等單位研究了縱向通風(fēng)下隧道長(zhǎng)度對(duì)臨界風(fēng)速的影響,結(jié)果顯示:臨界風(fēng)速隨著隧道長(zhǎng)度的增加呈上升趨勢(shì)但增速逐步變緩趨于0,并采用量綱分析法構(gòu)建出臨界風(fēng)速與隧道長(zhǎng)度關(guān)系公式[4];西南交通大學(xué)研究了隧道海拔高度對(duì)臨界風(fēng)速的影響,結(jié)果顯示:臨界風(fēng)速隨著海拔高度的增加而增大,當(dāng)火災(zāi)熱釋放率大于30 MW時(shí)海拔高度對(duì)臨界風(fēng)速影響較小,火災(zāi)熱釋放率小于30 MW時(shí)海拔高度對(duì)臨界風(fēng)速的影響顯著增強(qiáng),且隨著熱釋放率的減小影響不斷增大[5]。

        長(zhǎng)安大學(xué)王永東等人系統(tǒng)總結(jié)了國(guó)內(nèi)外關(guān)于火源熱釋放率、燃料類(lèi)型、火源尺寸、火源位置、縱向風(fēng)溫度、隧道阻塞比、縱向坡度、斷面形狀、線(xiàn)形曲率、海拔高度、排煙模式等不同因素對(duì)臨界風(fēng)速影響的研究成果[6],這些成果大部分都是基于FDS模擬試驗(yàn)得到。當(dāng)采用FDS來(lái)模擬多因素變化下臨界風(fēng)速變化規(guī)律時(shí),發(fā)現(xiàn)明顯存在時(shí)間成本過(guò)高的問(wèn)題,如在構(gòu)建好的50 m長(zhǎng)隧道內(nèi)完成單次5 s燃燒試驗(yàn)的運(yùn)行時(shí)間長(zhǎng)達(dá)2 h,且單次運(yùn)行時(shí)間會(huì)隨著隧道長(zhǎng)度、模擬燃燒時(shí)間以及影響因素個(gè)數(shù)的增長(zhǎng)顯著增加至數(shù)天。對(duì)此,有學(xué)者利用FDS軟件就其中五個(gè)影響因素分別進(jìn)行了研究,得到臨界風(fēng)速隨各影響因素的變化規(guī)律,在此基礎(chǔ)上引入BP神經(jīng)網(wǎng)絡(luò)技術(shù)研究多種影響因素下臨界風(fēng)速的預(yù)測(cè)建模,為發(fā)展預(yù)測(cè)臨界風(fēng)速的工程計(jì)算模型提供了一種新思路[7]。值得注意的是,該預(yù)測(cè)模型選取的因素中火源寬度、火源橫向位置對(duì)臨界風(fēng)速的影響小到可以忽略不計(jì),而其他研究證明對(duì)臨界風(fēng)速影響較大的海拔高度、斷面尺寸(緊急停車(chē)帶)等因素未在模型中得以體現(xiàn)。因此,研究一種適用于火災(zāi)工況下隧道通風(fēng)系統(tǒng)控制策略、無(wú)需FSD建模分析、充分考慮主要影響因素的臨界風(fēng)速快速預(yù)測(cè)技術(shù)十分必要。

        2.2 快速預(yù)測(cè)新技術(shù)總體思路

        該文提出一種臨界風(fēng)速快速預(yù)測(cè)新思路,其總體技術(shù)路線(xiàn)如圖3所示。

        主要步驟:

        (1)在FDS中構(gòu)建典型隧道模型,將一組火源功率、縱向坡度、隧道阻塞比等臨界風(fēng)速影響因素量化后作為試驗(yàn)參數(shù)進(jìn)行配置,并在模型中給定一個(gè)通風(fēng)初始風(fēng)速(約0.5 m/s),如圖4(a)所示。

        (2)運(yùn)行FDS模擬試驗(yàn),重點(diǎn)觀(guān)察試驗(yàn)過(guò)程中火源附近是否存在煙氣逆流現(xiàn)象,如果有煙氣逆流則增大通風(fēng)風(fēng)速后再次運(yùn)行FDS模擬試驗(yàn),直至風(fēng)速增大到恰好不存在煙氣逆流為止,如圖4(b)所示。

        (3)保存此時(shí)的各影響因素量化值和對(duì)應(yīng)的通風(fēng)風(fēng)速數(shù)據(jù),生成一個(gè)多因素—臨界風(fēng)速樣本。

        (4)如果還需生成更多的樣本數(shù)據(jù),則修改影響因素量化值及試驗(yàn)參數(shù),繼續(xù)執(zhí)行步驟(3)和(4)。

        (5)在生成足夠多數(shù)量的多因素—臨界風(fēng)速樣本后,歸集形成樣本庫(kù)。

        (6)采用機(jī)器學(xué)習(xí)算法來(lái)訓(xùn)練樣本庫(kù),建立基于多特征分析的火災(zāi)臨界風(fēng)速快速預(yù)測(cè)模型。

        (7)在臨界風(fēng)速預(yù)測(cè)時(shí),將火源功率、縱向坡度等各影響因素值輸入預(yù)測(cè)模型,即可獲得臨界風(fēng)速預(yù)測(cè)輸出結(jié)果。

        2.3 基于機(jī)器學(xué)習(xí)的模型構(gòu)建

        快速預(yù)測(cè)技術(shù)采用神經(jīng)網(wǎng)絡(luò)算法建立多影響因素作用的臨界風(fēng)速映射模型,以火源功率、縱向坡度、隧道阻塞比等作為神經(jīng)網(wǎng)絡(luò)的輸入計(jì)算及預(yù)測(cè)臨界風(fēng)速值。BP神經(jīng)網(wǎng)絡(luò)作為一種比較典型的人工神經(jīng)網(wǎng)絡(luò),是基于誤差反向傳播算法(Back-Propagation)的多層前向神經(jīng)網(wǎng)絡(luò),其主要結(jié)構(gòu)由一個(gè)輸入層、一個(gè)或多個(gè)隱含層和一個(gè)輸出層組成,各層由若干個(gè)神經(jīng)元構(gòu)成,每一個(gè)節(jié)點(diǎn)的輸出值由輸入值、激發(fā)函數(shù)和閾值決定。

        用于臨界風(fēng)速預(yù)測(cè)建模的三層BP神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)如圖5所示,其輸入層設(shè)計(jì)為8個(gè)節(jié)點(diǎn),分別對(duì)應(yīng)于樣本中的各輸入指標(biāo);網(wǎng)絡(luò)的隱含層沒(méi)有明確的物理意義,其節(jié)點(diǎn)個(gè)數(shù)根據(jù)經(jīng)驗(yàn)選取為輸入節(jié)點(diǎn)數(shù)的2倍左右,激發(fā)函數(shù)采用雙曲函數(shù);網(wǎng)絡(luò)的輸出層僅包含1個(gè)節(jié)點(diǎn),對(duì)應(yīng)于樣本中的臨界風(fēng)速值,由于該輸出不能取負(fù)值,其激發(fā)函數(shù)選取非負(fù)的Sigmoid函數(shù)。同時(shí),選取性能指標(biāo)函數(shù)并按照梯度下降學(xué)習(xí)法逐步調(diào)整BP神經(jīng)網(wǎng)絡(luò)的內(nèi)部參數(shù),獲得滿(mǎn)足精度要求的臨界風(fēng)速預(yù)測(cè)模型。

        在以上預(yù)測(cè)模型中,神經(jīng)網(wǎng)絡(luò)求解臨界風(fēng)速y的步驟如下:

        由輸入層第i個(gè)節(jié)點(diǎn)的輸出ai、輸入層第i個(gè)節(jié)點(diǎn)到隱含層第j個(gè)節(jié)點(diǎn)的權(quán)值wij、隱含層第j個(gè)節(jié)點(diǎn)的閾值θj計(jì)算隱含層第j個(gè)節(jié)點(diǎn)的輸入:

        由隱含層激發(fā)函數(shù)bj計(jì)算隱含層第j個(gè)節(jié)點(diǎn)的輸出:

        由隱含層第j個(gè)節(jié)點(diǎn)的輸出bi、隱含層第j個(gè)節(jié)點(diǎn)到輸出層節(jié)點(diǎn)的權(quán)值wj、輸出層節(jié)點(diǎn)的閾值θ計(jì)算輸出層節(jié)點(diǎn)的輸入net:

        由輸出層激發(fā)函數(shù)f(·)計(jì)算輸出層節(jié)點(diǎn)的輸出及臨界風(fēng)速值:

        3 結(jié)語(yǔ)

        該文針對(duì)火災(zāi)工況下公路隧道通風(fēng)系統(tǒng)的控制策略編制,一方面,在介紹斷面風(fēng)速測(cè)試現(xiàn)有方案后,分析了各類(lèi)風(fēng)速傳感器用于斷面風(fēng)速測(cè)量的優(yōu)缺點(diǎn)以及多傳感器常用的無(wú)線(xiàn)組網(wǎng)技術(shù),進(jìn)而提出適用于長(zhǎng)大隧道多斷面關(guān)鍵區(qū)域風(fēng)速監(jiān)測(cè)的新技術(shù)。另一方面,介紹了火災(zāi)工況下臨界風(fēng)速計(jì)算的現(xiàn)有方法及其存在的不足,提出臨界風(fēng)速快速預(yù)測(cè)新技術(shù),并對(duì)臨界風(fēng)速樣本獲取要點(diǎn)和多因素作用的臨界風(fēng)速預(yù)測(cè)模型構(gòu)建要點(diǎn)進(jìn)行了重點(diǎn)闡述,為通風(fēng)系統(tǒng)控制策略編制提供了強(qiáng)有力的理論支撐。圍繞火災(zāi)工況下通風(fēng)系統(tǒng)控制策略編制這一主題,開(kāi)展其他內(nèi)涵關(guān)鍵技術(shù)的深入研究,是后續(xù)工作的重點(diǎn)。

        參考文獻(xiàn)

        [1]李俊峰, 朱立偉. 公路隧道斷面風(fēng)速多通道檢測(cè)系統(tǒng)研究[J]. 公路交通科技(應(yīng)用技術(shù)版), 2012(4): 164-166.

        [2]王同喜, 李樹(shù)超, 齊楸, 等. 隧道風(fēng)速測(cè)試方法研究[J]. 消防科學(xué)與技術(shù), 2016(8): 1069-1072.

        [3]李俊梅, 許鵬, 李炎鋒, 等. 大坡度隧道臨界風(fēng)速的數(shù)值模擬和實(shí)驗(yàn)研究[J]. 北京工業(yè)大學(xué)學(xué)報(bào), 2017(11): 1706-1712.

        [4]姜學(xué)鵬, 朱珠, 喬意峰, 等. 縱向通風(fēng)下隧道長(zhǎng)度對(duì)臨界風(fēng)速的影響分析[J]. 地下空間與工程學(xué)報(bào), 2019(4): 1267-1272.

        [5]王峰, 張路華, 袁松, 等. 高海拔公路隧道火災(zāi)煙氣控制臨界風(fēng)速研究[J]. 中國(guó)公路學(xué)報(bào), 2022(5): 153-160.

        [6]王永東, 張如, 張一龍. 國(guó)內(nèi)外隧道火災(zāi)臨界風(fēng)速研究綜述[J]. 現(xiàn)代隧道技術(shù), 2018(4): 14-24.

        [7]梁華剛, 龐麗琴, 高冬梅. 基于ANN的公路隧道火災(zāi)臨界風(fēng)速預(yù)測(cè)研究[J]. 自動(dòng)化與儀表, 2019(1): 93-98.

        精品国产精品三级在线专区| 国产激情视频在线观看首页| 淫妇日韩中文字幕在线| 亚洲国产国语对白在线观看| 波多野结衣av一区二区全免费观看| 午夜成人无码福利免费视频| 亚洲国产精品嫩草影院久久| 中文字幕一区二区三在线| 穿着白丝啪啪的av网站| 久久人妻少妇嫩草av无码专区| 久久久精品免费观看国产| 91亚洲精品久久久蜜桃| 国产精品熟女一区二区三区| 中文字幕无码毛片免费看| 福利一区视频| 亚洲av一区二区三区网站| 美女脱了内裤露出奶头的视频| 99久久精品日本一区二区免费| 中文字幕一区二区三区乱码不卡| 亚洲女同av一区二区在线观看| 中文字幕无码乱人伦| 熟女人妻在线视频| 亚洲中文欧美日韩在线| 久久精见国产亚洲av高清热| 久久狠狠爱亚洲综合影院| 亚洲av永久无码精品秋霞电影影院| 国产丰满乱子伦无码专| 国产精品人伦一区二区三| 无码国产69精品久久久久孕妇| 国产中文制服丝袜另类| 国产精品亚洲av无人区一区蜜桃| 国产精品偷窥熟女精品视频| 欧美gv在线观看| 美女视频很黄很a免费国产| 午夜视频一区二区三区四区| 免费人妻无码不卡中文字幕系| 国产亚洲亚洲精品777| 一本久道视频无线视频试看 | 欧美性高清另类videosex| 亚洲国产综合精品 在线 一区| 国产精品不卡无码AV在线播放 |