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

        ?

        特厚易自燃煤層回采工作面上隅角貧氧致因機理研究

        2022-11-09 02:29:22高瑞青劉紅威張小龍賀志宏劉振明王子邦
        煤礦安全 2022年10期

        高瑞青,王 飛,劉紅威,張小龍,賀志宏,劉振明,王子邦,宮 彪

        (1.太原理工大學(xué)安全與應(yīng)急管理工程學(xué)院,山西太原 030024;2.山西煤礦安全研究生教育創(chuàng)新中心,山西太原030024;3.西山煤電(集團)有限責(zé)任公司,山西太原 030000;4.晉能控股煤業(yè)集團有限公司,山西大同 037003)

        隨著煤炭開采技術(shù)水平的提升,我國煤礦綜采工作面采高不斷增加[1],在特厚自燃煤層放頂煤綜采過程中,采空區(qū)自然發(fā)火嚴重[2-3],工作面上隅角貧氧現(xiàn)象時有發(fā)生。根據(jù)GB 8958—88《缺氧危險作業(yè)安全規(guī)程》中定義在氧氣體積分數(shù)低于19.5%的工作環(huán)境作業(yè)時為缺氧作業(yè)。當氧氣體積分數(shù)低于18.5%工作能力下降,劇烈運動感到不適[4]??紤]煤礦井下環(huán)境復(fù)雜,人員工作強度大,以18.5%為界定義當工作面氧氣體積分數(shù)低于18.5%的現(xiàn)象稱為工作面貧氧現(xiàn)象[5]。為保證工人的安全生產(chǎn),對工作面上隅角貧氧致因機理的研究具有重要意義。

        為掌握特厚易自燃煤層回采工作面上隅角貧氧機理,研究人員進行了相關(guān)研究。郭永文[6]針對神東礦區(qū)工作面低氧現(xiàn)象提出采空區(qū)遺煤低溫氧化是低氧氣體產(chǎn)生的來源;潘榮錕等[7]采用CDEM 數(shù)值模擬以及SF6示蹤氣體漏風(fēng)測試等方法對淺埋近距離煤層群工作面上隅角貧氧致因進行了研究,指出工作面上隅角貧氧是多因素影響下的綜合結(jié)果,開采擾動產(chǎn)生的漏風(fēng)裂隙,是造成上隅角貧氧的主要致因;方保明[8]提出造成工作面低氧的主要原因是煤層賦存于CO2-N2帶,回采過程中高氮氣體涌出與采空區(qū)內(nèi)遺煤低溫氧化形成高氮低氧環(huán)境。現(xiàn)場研究還發(fā)現(xiàn)采空區(qū)二氧化碳涌出,外界大氣壓變化也是造成上隅角氧氣體積分數(shù)低的重要原因。上述研究表明工作面上隅角貧氧現(xiàn)象的出現(xiàn)不是單一因素造成的,是多因素耦合綜合作用的結(jié)果。

        目前大范圍采場氣體分布的研究主要采用數(shù)值模擬和現(xiàn)場測試相結(jié)合的手段[9-11]?;诖?,以馬道頭礦8105 工作面為工程背景,建立特厚易自燃煤層回采工作面流場數(shù)值計算模型,分析采空區(qū)耗氧速率、工作面進風(fēng)量、采空區(qū)埋管注氮(注氮量和注氮位置)對采場氧氣體積分數(shù)分布規(guī)律的影響效應(yīng),明晰工作面上隅角的貧氧致因機理。

        1 上隅角貧氧工程實例

        馬道頭礦井通風(fēng)方式為混合式,通風(fēng)方法為機械抽出式。8105 綜放工作面主要開采中5#、(3-5)#煤層,埋深300~460 m,平均傾角3°,走向長度2 288 m,傾斜長度220 m,總厚度為5.76~22.98 m,平均煤厚17.8 m,平均落底煤厚度1.17 m。8105 工作面采用一進一回的“U”型通風(fēng)方式,全負壓通風(fēng);所在煤層為易自燃煤層,自然發(fā)火期41 d;采用向采空區(qū)24 h 不間斷注氮為主的方式預(yù)防采空區(qū)遺煤自燃。在8105 工作面3 種生產(chǎn)狀態(tài)下對上隅角距煤壁1.5 m 處,分別在高度為2.5、1.5、0.8 m 位置用橡膠氣囊取低氧氣體氣樣,在實驗室用氣相色譜分析。8105工作面上隅角低氧氣體成分分析表見表1。由表1可知,生產(chǎn)過程中工作面附近存在嚴重的貧氧區(qū)域,對工作面的安全生產(chǎn)造成了極大的威脅。

        表1 8105 工作面上隅角低氧氣體成分分析表Table 1 Analysis table of low-oxygen gas composition in the upper corner of working surface

        2 計算模型

        2.1 模型建立與參數(shù)設(shè)定

        考慮到采空區(qū)中復(fù)雜的流場分布,提出合理假設(shè)[12]:視采場內(nèi)氣體為不可壓縮流體;采空區(qū)孔隙率及滲透率為采空區(qū)位置坐標的函數(shù),不隨時間發(fā)生改變;考慮重力對采空區(qū)流場的影響;忽略溫度及煤層傾角對采空區(qū)流場分布的影響。

        2.1.1 數(shù)學(xué)模型

        采空區(qū)內(nèi)煤巖受工作面開采影響,在采空區(qū)內(nèi)形成自由堆積的多孔介質(zhì)區(qū)域,流體在采空區(qū)內(nèi)的運動符合多孔介質(zhì)滲流理論。流體在采空區(qū)內(nèi)的運動滿足多孔介質(zhì)質(zhì)量及動量守恒方程,多孔介質(zhì)內(nèi)的質(zhì)量守恒方程為[13]:

        式中:γ 為多孔介質(zhì)空隙率,m-3;ρf為流體質(zhì)量密度,kg/m3;▽為哈密算子;v→為流體在多孔介質(zhì)中滲流速度,m/s;Sm為氣體質(zhì)量源項,kg/(m3·s);t 為時間,s。

        多孔介質(zhì)內(nèi)的動量守恒方程為:

        式中:p 為氣體壓力,Pa;τ=為黏性應(yīng)力張量;g→為重力加速度,m/s2;S 為多孔介質(zhì)中風(fēng)流運移過程中附加的動量損失源項,N/m3。

        對各向同性的多孔介質(zhì)動量損失源項為[14]:

        式中:μ 為流體的動力黏度系數(shù),kg/(m·s);α 為滲透率;vi為i 方向的速度;C2為慣性阻力系數(shù)。

        根據(jù)Blake-Kozeny 推導(dǎo)思想,在非線性流態(tài)條件下定義適用的經(jīng)驗公式為[15]:

        式中:Dm為采空區(qū)內(nèi)落巖的平均粒徑,m。

        氧氣擴散方程:

        式中:ρ 為流體密度,kg/m3;Yi為氧氣體積分數(shù),%;J→i為氧氣的擴散通量;Ri為氧氣化學(xué)反應(yīng)生成凈速率,kg/(m3·s);Si為氧氣生成率,kg/(m3·s)。

        J→i為氧氣的擴散通量,由Fick 定律給出,湍流模型中表達式為:

        式中:Di,m為氧氣的質(zhì)量擴散系數(shù);μt為湍動黏度,kg/(m·s);Sct為湍流施密特數(shù);DT,i為氧氣組分的熱擴散系數(shù);T 為溫度,K。

        2.1.2 幾何模型與邊界條件

        8105 進風(fēng)巷與工作面交界點為坐標原點,工作面傾向長度為y 軸方向,沿著采空區(qū)走向為x 軸方向,進風(fēng)巷、回風(fēng)巷斷面尺寸為5.1 m×3.9 m;工作面尺寸為6.864 m×3.9 m,傾向長度為220 m;采空區(qū)走向長度為300 m;高度為48 m。工作面采用“U”型通風(fēng),進風(fēng)量為2 100 m3/min,配合進風(fēng)巷雙埋管注氮距工作面30、60 m,24 h 不間斷注氮,總注氮量為1 745 m3/h,邊界條件和參數(shù)設(shè)置見表2。8105 工作面物理模型如圖1。采空區(qū)空隙率分布函數(shù)為[16]:

        表2 邊界條件和參數(shù)設(shè)置Table 2 Boundary conditions and parameter settings

        圖1 8105 工作面物理模型Fig.1 Physical model of 8105 working face

        式中:x 為采空區(qū)走向方向距離,m;y 為工作面傾向方向距離,m;z 為采空區(qū)高度,m;l 為工作面總長度,m。

        z=0 m 平面采空區(qū)空隙率分布如圖2。

        圖2 z=0 m 平面采空區(qū)空隙率分布Fig.2 Porosity distribution of goaf on z=0 m plane

        2.2 計算模型準確性校驗

        8105 模型工作面測點布置示意圖如圖3。利用后處理工具對工作面測點及束管監(jiān)測處的數(shù)據(jù)進行截取,通過模擬結(jié)果與實測數(shù)據(jù)對比分析。工作面氧氣體積分數(shù)實測與模擬結(jié)果對比如圖4。束管監(jiān)測氧氣體積分數(shù)與模擬結(jié)果對比如圖5。

        圖3 8105 模型工作面測點布置示意圖Fig.3 Schematic diagram of measuring point layout of 8105 model working face

        圖4 工作面氧氣體積分數(shù)實測與模擬結(jié)果對比Fig.4 Comparison of actual measurement and simulation results of oxygen volume fraction in working face

        圖5 束管監(jiān)測氧氣體積分數(shù)與模擬結(jié)果對比Fig.5 Comparison of tube bundle monitoring oxygen volume fraction and simulation results

        由圖4 可知:當工作面靠近進風(fēng)側(cè)y=0~20 m時,模擬結(jié)果與工作面實測分別為21%與20.9%;當y=20~50 m 時,受進風(fēng)側(cè)埋管注氮影響,氧氣體積分數(shù)小范圍急速下降;當y=50~115 m 時,氧氣體積分數(shù)逐步平緩下降;當y=115~210 m 時,采空區(qū)氣體涌入工作面,氧氣體積分數(shù)下降趨勢顯著;當y>210 m,由于工作面靠近回風(fēng)巷側(cè)壓力小,采空區(qū)氣體加速涌入工作面,氧氣體積分數(shù)急劇下降。工作面氧氣體積分數(shù)分布與現(xiàn)場實測結(jié)果相對一致。

        由圖5 可知:模擬顯示距離工作面0~10 m 時,氧氣體積分數(shù)從17.9%驟降至12.1%,距工作面10~80 m 范圍內(nèi),氧氣體積分數(shù)由12.1%平緩降6.4%;束管實測數(shù)據(jù)顯示在距工作面5~10 m 位置,氧氣體積分數(shù)急劇下降,在距工作面10~75 m 處氧氣體積分數(shù)緩步下降至6.7%,在75~80 m 驟然下降至6%;模擬結(jié)果與實測數(shù)據(jù)變化趨勢吻合較好。。

        3 回采工作面上隅角貧氧致因機理

        3.1 耗氧速率對采場氧氣體積分數(shù)分布的影響

        通過改變采空區(qū)耗氧速率,模擬分析工作面及上隅角區(qū)域氧氣體積分數(shù)分布的變化規(guī)律。不同耗氧速率下工作面測點氧氣體積分數(shù)分布如圖6。

        圖6 不同耗氧速率下工作面測點氧氣體積分數(shù)分布Fig.6 Oxygen volume fraction at measuring points in working faces under different oxygen consumption rates

        由圖6 可知:當耗氧速率為0 kg/(m3·s)時,工作面氧氣體積分數(shù)基本穩(wěn)定在21%,隨著耗氧速率的增加工作面測點0~200 m 測點氧氣體積分數(shù)逐步下降,耗氧速率與工作面氧氣體積分數(shù)呈現(xiàn)負相關(guān);當耗氧速率變化梯度較大時,在5×10-7~1.2×10-6kg/(m3·s),工作面氧氣體積分數(shù)急劇降低;當耗氧速率變化梯度最小,在1×10-8~1×10-7kg/(m3·s),工作面測點氧氣體積分數(shù)變化較?。缓难跛俾试酱?,工作面測點氧氣體積分數(shù)下降梯度越大。

        z=1.7 m 不同耗氧速率上隅角氧氣體積分數(shù)分布如圖7。不同耗氧速率上隅角貧氧效應(yīng)統(tǒng)計如圖8。

        圖7 z=1.7 m 不同耗氧速率上隅角氧氣體積分數(shù)分布Fig.7 Oxygen volume fraction distribution in upper corner(z=1.7 m)under different oxygen consumption rates

        圖8 不同耗氧速率上隅角貧氧效應(yīng)統(tǒng)計Fig.8 Statistics of oxygen depletion effect in upper corner of different oxygen consumption rates

        由圖7 與圖8 可知:當耗氧速率為0 kg/(m3·s)時,上隅角與下隅角氧氣體積分數(shù)基本一致,氧氣體積分數(shù)穩(wěn)定在20.8%以上;耗氧速率增至1.8×10-7kg/(m3·s)時,上隅角區(qū)域氧氣體積分數(shù)逐步下降,氧氣體積分數(shù)均大于18.5%,未出現(xiàn)貧氧現(xiàn)象;當耗氧速率大于1×10-7kg/(m3·s)后,上隅角貧氧現(xiàn)象出現(xiàn),隨著耗氧速率的增加,上隅角氧氣體積分數(shù)不斷降低,貧氧區(qū)域在上隅角附近逐步擴大。

        綜合分析可知采空區(qū)遺煤氧化是造成工作面上隅角貧氧現(xiàn)象的前提條件,采空區(qū)遺煤氧化加劇,耗氧速率增加,產(chǎn)生大量的低氧體積分數(shù)氣體,工作面負壓通風(fēng)方式幫助了采空區(qū)內(nèi)貧氧氣體運移到靠近工作面一側(cè),在上隅角聚集造成上隅角貧氧現(xiàn)象。

        3.2 進風(fēng)量對采場氧氣體積分數(shù)分布規(guī)律的影響

        通過改變巷道進風(fēng)量,模擬分析了工作面與上隅角區(qū)域氧氣體積分數(shù)分布的變化規(guī)律。不同進風(fēng)量下工作面測點氧氣體積分數(shù)分布如圖9。

        圖9 不同進風(fēng)量下工作面測點氧氣體積分數(shù)分布Fig.9 Oxygen volume fraction at measuring points in working faces under different air intake

        由圖9 可知:當進風(fēng)量大于2 100 m3/min 時,工作面測點0~200 m 氧氣體積分數(shù)穩(wěn)定在20.7%以上;進風(fēng)量增大,采空區(qū)涌入工作面的低氧氣體及時被風(fēng)流帶走,工作面氧氣體積分數(shù)顯著提升,隨著進風(fēng)量的增加工作面測點氧氣體積分數(shù)逐步提高;同時進風(fēng)量的增加,使進風(fēng)側(cè)壓力增大,采空區(qū)漏風(fēng)加劇,采空區(qū)與工作面回風(fēng)側(cè)壓差變大,采空區(qū)低氧體積分數(shù)氣體在工作面尾部區(qū)域(y>200 m)加速涌入,氧氣體積分數(shù)迅速下降;當工作面進風(fēng)量降至1 874 m3/min 以下,采空區(qū)氣體通過漏風(fēng)裂隙涌入工作面,由于工作面風(fēng)速較小不能及時稀釋低氧體積分數(shù)氣體造成工作面測點氧氣體積分數(shù)迅速降低。

        z=1.7 m 不同進風(fēng)量上隅角氧氣體積分數(shù)分布如圖10,不同進風(fēng)量上隅角貧氧效應(yīng)統(tǒng)計如圖11。

        由圖10 與圖11 可知:當工作面進風(fēng)量小于2 100 m3/min 時,工作面漏風(fēng)較小,采空區(qū)遺煤氧化減緩,低氧體積分數(shù)氣體減少,上隅角處氧氣體積分數(shù)均大于18.5%;隨著進風(fēng)量的增加上隅角處最低氧氣體積分數(shù)緩步下降,當工作面進風(fēng)量在1 648~2 100 m3/min 時,最低氧氣體積分數(shù)急劇下降,進風(fēng)量2 100~2 552 m3/min 之間時,上隅角出現(xiàn)貧氧現(xiàn)象,低氧范圍迅速擴大;當進風(fēng)量大于2 552 m3/min后,隨著進風(fēng)量的增加,低氧范圍趨于基本穩(wěn)定,上隅角氧氣體積分數(shù)緩步下降。

        圖10 z=1.7 m 不同進風(fēng)量上隅角氧氣體積分數(shù)分布Fig.10 Oxygen volume fraction distribution on upper corner(z=1.7 m)under different air intake

        圖11 不同進風(fēng)量上隅角貧氧效應(yīng)統(tǒng)計Fig.11 Statistical analysis of oxygen deficiency effect in upper corner of different air intake

        綜合分析可知隨著進風(fēng)量的增加,工作面漏風(fēng)增大,采空區(qū)內(nèi)遺煤氧化加劇產(chǎn)生大量的低氧體積分數(shù)氣體,同時采空區(qū)與工作面回風(fēng)側(cè)的壓差變大,采空區(qū)內(nèi)氣體在工作面尾部加速涌入,在上隅角渦流區(qū)域聚集加劇了上隅角貧氧現(xiàn)象。

        3.3 注氮對采場氧氣體積分數(shù)分布規(guī)律的影響

        3.3.1 工作面及上隅角氧氣體積分數(shù)分布

        當注氮位置為30、60 m 時,不同注氮量下工作面測點氧氣體積分數(shù)分布如圖12。

        圖12 不同注氮量下工作面測點氧氣體積分數(shù)分布Fig.12 Oxygen volume fraction at measuring points of working face under different nitrogen injection rates

        由圖12 可知:注氮量的變化對工作面測點氧氣體積分數(shù)分布有較大影響;當注氮流量在1 385~2 465 m3/h 時,采空區(qū)注氮口與工作面壓差均衡,增加注氮量對工作面氧氣體積分數(shù)分布影響較?。划斪⒌繌? 825 m3/h 增加至3 545 m3/h 時,注氮口附近區(qū)域與工作面壓差增大,氮氣通過漏風(fēng)通道涌入工作面,造成工作面測點氧氣體積分數(shù)迅速下降。

        z=1.7 m 不同注氮量工作面氧氣體積分數(shù)分布如圖13,不同注氮量上隅角貧氧效應(yīng)統(tǒng)計如圖14。

        圖13 z=1.7 m 不同注氮量工作面氧氣體積分數(shù)分布Fig.13 Oxygen volume fraction distribution on working face(z=1.7 m)under different nitrogen injection rates

        由圖13 與圖14 可知:隨著注氮量的增加,工作面測點氧氣體積分數(shù)逐漸下降,上隅角低氧范圍逐步擴大;當注氮量超過2 825 m3/h 后,注氮口與工作面壓力失衡,采空區(qū)低氧體積分數(shù)氣體大量涌入工作面,注氮口附近工作面氧氣體積分數(shù)降至18.5%以下造成工作面低氧。當注入氮量從0 m3/h 增至3 185 m3/h 上隅角氧氣體積分數(shù)逐步下降,上隅角低氧區(qū)域逐步擴大,當注氮量增加至3 545 m3/h 時,注氮口壓力升高,與工作面壓差增大,氮氣通過漏風(fēng)通道涌入工作面,工作面氧氣體積分數(shù)急劇下降,同時氮氣大量涌入工作面,在上隅角渦流區(qū)域聚集,造成上隅角低氧區(qū)域急劇擴大。

        圖14 不同注氮量上隅角貧氧效應(yīng)統(tǒng)計Fig.14 Statistics of oxygen-deficient effects in the upper corner of different nitrogen injection rates

        因此自燃礦井在采用注氮技術(shù)治理采空區(qū)遺煤自燃時,應(yīng)合理設(shè)置注氮量大小,減輕工作面低氧危險性,兼顧防治采空區(qū)遺煤自燃的發(fā)生,從而達到工作面低氧與煤自燃協(xié)調(diào)防治的目的。

        3.3.2 工作面及上隅角貧氧分析

        當注氮量為1 745 m3/h,分別模擬不同注氮位置時,工作面與上隅角區(qū)域氧氣體積分數(shù)分布的變化規(guī)律。不同注氮位置下工作面測點氧氣體積分數(shù)分布如圖15。不同注氮位置上隅角貧氧效應(yīng)統(tǒng)計如圖16。

        圖15 不同注氮位置下工作面測點氧氣體積分數(shù)分布Fig.15 Oxygen volume fraction at measuring points of working face under different nitrogen injection positions

        圖16 不同注氮位置上隅角貧氧效應(yīng)統(tǒng)計Fig.16 Statistics of oxygen depletion effect in upper corner of different nitrogen injection positions

        由圖15 可知:隨著注氮氣位置的不斷深入,工作面測點氧氣體積分數(shù)分布變化較小。綜合分析注氮位置的深入對工作面測點氧氣體積分數(shù)影響較小。

        由圖16 可知:當注氮量為1 745 m3/h 隨著注氮位置的變化,上隅角氧氣體積分數(shù)和低氧區(qū)域范圍均無顯著變化。即注氮位置的深入對工作面及上隅角氧氣體積分數(shù)分布有著微弱的影響效應(yīng)。

        4 結(jié) 論

        1)采空區(qū)耗氧速率為0 kg/(m3·s)時,工作面進風(fēng)側(cè)與回風(fēng)側(cè)氧氣體積分數(shù)基本一致。采空區(qū)耗氧速率大于0 kg/(m3·s)時,回風(fēng)側(cè)氧氣體積分數(shù)小于進風(fēng)側(cè)。采空區(qū)遺煤氧化是造成工作面上隅角貧氧現(xiàn)象的前提條件。

        2)隨著耗氧速率的升高,采空區(qū)遺煤氧化加劇產(chǎn)生大量的低氧體積分數(shù)氣體,通過漏風(fēng)裂隙涌入工作面氣體氧體積分數(shù)降低,工作面及上隅角氧氣體積分數(shù)逐漸下降,貧氧區(qū)域逐步擴大。

        3)隨著工作面進風(fēng)量的增加,工作面氧氣體積分數(shù)升高,同時采空區(qū)漏風(fēng)量增大,采空區(qū)內(nèi)遺煤氧化反應(yīng)加劇產(chǎn)生大量的低氧體積分數(shù)氣體,由于采空區(qū)內(nèi)與工作面回風(fēng)側(cè)壓差變大,采空區(qū)氣體在工作面尾部加速涌入,上隅角附近氧氣體積分數(shù)急劇降低,貧氧范圍迅速擴大。當進風(fēng)量大于2 778 m3/min時,上隅角貧氧區(qū)域趨于穩(wěn)定,但氧氣體積分數(shù)仍逐步下降。

        4)采空區(qū)埋管注氮在防滅火的同時為采空區(qū)內(nèi)提供了大量的氮氣,氮氣通過漏風(fēng)通道擴散到工作面中稀釋了風(fēng)流中的氧氣含量,隨著注氮量的增加,工作面及上隅角氧氣體積分數(shù)逐步降低。當注氮量大于2 825 m3/h 時,注氮口與工作面壓力失衡,大量氮氣涌入工作面造成工作面貧氧;當注氮量為1 745 m3/h 時,注氮位置的不斷深入,對工作面與上隅角氧氣體積分數(shù)分布影響較小。

        无码喷潮a片无码高潮| 扒开双腿操女人逼的免费视频| 九九久久精品一区二区三区av| 午夜人妻久久久久久久久| 无码av无码天堂资源网| 综合色久七七综合尤物| 亚洲女同恋中文一区二区| 国产精品国产三级国产密月| 无遮挡h肉动漫在线观看| 可以免费在线看黄的网站| 精品国产麻豆一区二区三区| 亚洲一区二区三区,日本| 日韩毛片免费无码无毒视频观看| 91在线在线啪永久地址| 精品专区一区二区三区| 一个色综合中文字幕人妻激情视频| 中文字幕一区二区人妻性色| 国产精品原创巨作av无遮| 日韩字幕无线乱码免费| 国产无套乱子伦精彩是白视频 | 精品国产亚洲av麻豆尤物| 北条麻妃在线中文字幕| 国产成人无码精品久久二区三区 | 一区二区在线观看日本免费| 狠狠综合亚洲综合亚洲色| 久久久久久久久久久国产| 精品国产午夜久久久久九九| 国产人妖视频一区二区| 亚洲av美国av产亚洲av图片| 久久AV中文一区二区三区 | 久久97久久97精品免视看| 亚洲av无码成人精品区天堂 | 97久久久久人妻精品专区| 人妻丝袜中文字幕久久| 亚洲一区二区免费在线观看视频| 18禁黄久久久aaa片| 亚洲国产精品500在线观看| 乳乱中文字幕熟女熟妇| 亚洲国产av玩弄放荡人妇| 久久噜噜噜| 在线观看二区视频网站二区|