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

        ?

        掘進(jìn)工作面瓦斯動(dòng)態(tài)涌出規(guī)律數(shù)值模擬研究

        2017-05-02 09:33:18趙東云
        中國(guó)煤炭 2017年4期
        關(guān)鍵詞:煤壁煤體滲透率

        楊 濤 趙東云

        (1.華北科技學(xué)院安全工程學(xué)院,河北省三河市,101601; 2. 煤炭資源與安全開(kāi)采國(guó)家重點(diǎn)實(shí)驗(yàn)室,北京市海淀區(qū),100083)

        ?

        ★ 煤礦安全 ★

        掘進(jìn)工作面瓦斯動(dòng)態(tài)涌出規(guī)律數(shù)值模擬研究

        楊 濤1,2趙東云1

        (1.華北科技學(xué)院安全工程學(xué)院,河北省三河市,101601; 2. 煤炭資源與安全開(kāi)采國(guó)家重點(diǎn)實(shí)驗(yàn)室,北京市海淀區(qū),100083)

        為研究掘進(jìn)過(guò)程瓦斯涌出規(guī)律,基于有限體積法的C-N格式離散方法,編制程序并對(duì)移動(dòng)掘進(jìn)工作面巷道周圍瓦斯壓力分布、煤體透氣性變化及瓦斯涌出速率變化情況等進(jìn)行了數(shù)值模擬。結(jié)果表明,隨著工作面向前掘進(jìn),掘進(jìn)工作面一定區(qū)域內(nèi)巷道圍巖煤體瓦斯壓力分布向前移動(dòng)。掘進(jìn)過(guò)程中巷幫瓦斯壓力存在轉(zhuǎn)折點(diǎn),轉(zhuǎn)折點(diǎn)之前的瓦斯壓力梯度減小速率大于轉(zhuǎn)折點(diǎn)之后的瓦斯壓力梯度。不同掘進(jìn)長(zhǎng)度下巷道沿程煤壁瓦斯涌出速率曲線基本一致,且瓦斯涌出速率較大區(qū)域主要集中在掘進(jìn)頭巷道部分。掘進(jìn)過(guò)程中瓦斯涌出速率呈鋸齒狀變化,每一個(gè)掘進(jìn)循環(huán)中掘進(jìn)頭巷道部分剛被掘開(kāi)時(shí)瓦斯涌出異常劇烈。隨著掘進(jìn)巷道增長(zhǎng),瓦斯涌出速率總體上逐漸增大。

        掘進(jìn)工作面 瓦斯壓力 瓦斯涌出速率 滲透率 數(shù)值模擬

        瓦斯災(zāi)害是常見(jiàn)的礦井災(zāi)害之一,而掘進(jìn)工作面瓦斯災(zāi)害尤為突出。煤巷掘進(jìn)過(guò)程中揭露煤層會(huì)伴有大量煤體瓦斯解吸,加之通風(fēng)難度較大,常引起局部瓦斯積聚,給采掘作業(yè)帶來(lái)諸多困難。掘進(jìn)工作面瓦斯涌出規(guī)律受煤層瓦斯初始?jí)毫?、圍巖應(yīng)力狀態(tài)及掘進(jìn)工藝等因素影響,研究這些因素影響下的煤體瓦斯涌出規(guī)律,對(duì)于控制掘進(jìn)工作面瓦斯災(zāi)害具有重要意義。

        國(guó)內(nèi)外學(xué)者圍繞煤體中瓦斯?jié)B流過(guò)程及機(jī)理開(kāi)展了大量研究工作,通過(guò)理論分析建立了相關(guān)理論模型,有學(xué)者通過(guò)編制程序計(jì)算得到了相應(yīng)的物理解,也有學(xué)者利用數(shù)值計(jì)算軟件(FLUENT、COMSOL等)建立多孔介質(zhì)滲流模型開(kāi)展了相關(guān)模擬計(jì)算。但現(xiàn)有的計(jì)算程序在運(yùn)算時(shí)間和誤差控制方面仍有不足,數(shù)值模擬軟件的計(jì)算又難以綜合考慮采掘過(guò)程圍巖三區(qū)應(yīng)力動(dòng)態(tài)變化對(duì)瓦斯流動(dòng)的影響,本文在前人研究的基礎(chǔ)上通過(guò)改進(jìn)模型、優(yōu)化算法對(duì)掘進(jìn)工作面的瓦斯涌出規(guī)律進(jìn)行了研究。

        1 掘進(jìn)巷道煤壁瓦斯流動(dòng)模型

        1.1 物理模型

        煤巷掘進(jìn)過(guò)程中,隨著新鮮煤壁的揭露和圍巖應(yīng)力的變化,煤體內(nèi)瓦斯不斷向巷道涌出,瓦斯在煤層厚度方向上的運(yùn)移居于次要地位,可將煤層內(nèi)瓦斯運(yùn)移看作二維平面瓦斯流動(dòng)問(wèn)題。煤層巷道掘進(jìn)物理模型示意圖如圖1所示。

        圖1 煤層巷道掘進(jìn)平面圖

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

        1.2.1 基本假設(shè)

        (1)煤層為多孔介質(zhì)結(jié)構(gòu),假定煤層各向同性且非均質(zhì),同時(shí)認(rèn)為近煤壁處煤層各處煤體滲透率與瓦斯壓力不相關(guān)。

        (2)揭露煤層時(shí)認(rèn)為煤層各處瓦斯解吸過(guò)程瞬間完成,且解吸熱不予考慮,吸附瓦斯量與瓦斯壓力始終滿足朗格繆爾等溫吸附定律:

        (1)

        式中:a,b——Langmuir等溫吸附常數(shù);

        P——煤層內(nèi)部瓦斯壓力平方,MPa2;

        Q——單位體積煤體內(nèi)瓦斯含量,m3/ m3;

        ρm——煤的密度,t/m3;

        W——煤層的水分含量,%;

        Ad——煤層的灰分含量,%;

        c——修正系數(shù)。

        (3)瓦斯在煤層內(nèi)運(yùn)移屬達(dá)西流,符合達(dá)西定律,考慮瓦斯壓力變化ΔP與流量關(guān)系時(shí),給定溫度下單位體積瓦斯體積流量q為:

        q=-λΔP

        (2)

        式中:λ——煤層透氣性系數(shù),m2/(MPa2· s);

        (4)煤層頂?shù)装鍨橹旅苄詭r層,不透氣且不含瓦斯。

        1.2.2 巷道圍巖滲透率分布

        煤巷掘進(jìn)過(guò)程中煤層圍巖的原巖應(yīng)力狀態(tài)被打破,煤壁前方形成應(yīng)力三區(qū),煤體受到破壞,在卸壓區(qū)內(nèi)形成大量連通孔、裂隙結(jié)構(gòu),煤體滲透率隨之增大;在彈性區(qū)內(nèi),煤體原生裂隙在高應(yīng)力作用下閉合,氣體運(yùn)移通道減少,煤體滲透率降低;在塑性區(qū)內(nèi),煤體受高應(yīng)力作用發(fā)生破壞,產(chǎn)生大量新的裂隙網(wǎng)絡(luò),煤體滲透率增大。掘進(jìn)巷道煤體滲透率控制方程可表示為:

        (3)

        式中:l——煤體距巷道煤壁的最小距離,m;

        k0——煤層的原始瓦斯壓力, m2/ (MPa2·d)。

        利用VB程序計(jì)算后得滲透率分布規(guī)律如圖2所示。

        圖2 掘進(jìn)巷道圍巖煤體滲透率分布

        1.2.3 控制方程及求解條件

        將瓦斯在巷道圍巖煤體中流動(dòng)認(rèn)為是二維非穩(wěn)態(tài)流動(dòng),控制方程如下:

        (4)

        模型計(jì)算初始條件:當(dāng)t=0時(shí),煤層各處的瓦斯壓力為煤層原始瓦斯壓力,即p(x,y,0) =p0,取值為1.26 MPa(取振興二礦2#煤層瓦斯原始?jí)毫?,揭露后煤壁處瓦斯壓力與巷道內(nèi)的大氣壓力相當(dāng),即為0.10325 MPa。

        模型計(jì)算邊界條件:新揭露煤壁處瓦斯壓力為巷道內(nèi)的大氣壓力,即為0.10325 MPa;計(jì)算模型邊界處瓦斯壓力梯度為0,即:

        (5)

        2 掘進(jìn)巷道煤壁瓦斯?jié)B流數(shù)值模擬

        2.1 網(wǎng)格劃分

        計(jì)算時(shí),考慮掘進(jìn)工作按周期規(guī)律有序進(jìn)行,掘進(jìn)后揭露煤層逐漸由計(jì)算區(qū)域轉(zhuǎn)為邊界區(qū)域而不再屬于計(jì)算范圍;工作面推進(jìn)過(guò)程,計(jì)算區(qū)域?qū)?huì)逐步向揭露煤層內(nèi)部擴(kuò)展,以保證計(jì)算區(qū)域內(nèi)煤體始終處于原巖應(yīng)力狀態(tài),此時(shí)可將煤層內(nèi)部瓦斯壓力認(rèn)為是煤層原始瓦斯壓力。計(jì)算區(qū)域采用結(jié)構(gòu)式矩形網(wǎng)格劃分。

        2.2 控制方程離散化

        由質(zhì)量守恒定律知,單位時(shí)間內(nèi),模型中任一控制體積內(nèi),瓦斯流入或涌出質(zhì)量之差等于該控制體積內(nèi)瓦斯質(zhì)量變化量,可對(duì)控制方程進(jìn)行離散化處理。

        計(jì)算區(qū)域內(nèi),依據(jù)網(wǎng)格劃分形式,對(duì)瓦斯?jié)B流方程進(jìn)行有限體積C-N格式離散,網(wǎng)格中每個(gè)節(jié)點(diǎn)的離散方程均考慮有10個(gè)瓦斯壓力值,包括任一節(jié)點(diǎn)n0當(dāng)前時(shí)刻t0時(shí)周邊4個(gè)節(jié)點(diǎn)瓦斯壓力值、上一時(shí)刻t0-1時(shí)該節(jié)點(diǎn)瓦斯壓力值及與之相鄰的4個(gè)節(jié)點(diǎn)的瓦斯壓力值。其中,t0-1時(shí)該瓦斯壓力值為已知,t0時(shí)刻瓦斯壓力值為未知。根據(jù)控制方程組系數(shù)矩陣及常數(shù)矩陣中所包含的待求未知量,采用迭代求解的方法計(jì)算t0時(shí)刻煤體瓦斯壓力分布。計(jì)算時(shí)設(shè)置計(jì)算誤差值在小于規(guī)定誤差值時(shí),才進(jìn)行下一時(shí)刻各節(jié)點(diǎn)瓦斯壓力值迭代計(jì)算。

        煤層瓦斯?jié)B流微分方程離散方程經(jīng)轉(zhuǎn)化后得到方程組形式見(jiàn)式(6)。

        (6)

        Δyi,j——節(jié)點(diǎn)(i,j)控制體長(zhǎng),m;

        Δxi,j——節(jié)點(diǎn)(i,j)控制體長(zhǎng),m;

        Δt——時(shí)間步長(zhǎng),s;

        λi+0.5,j、λi,j-0.5、λi,j+0.5、λi-0.5,j——節(jié)點(diǎn)(i,j)與(i+1,j)間、節(jié)點(diǎn)(i,j)與(i,j-1)間、節(jié)點(diǎn)(i,j)與(i,j+1)間、節(jié)點(diǎn)(i,j)與(i-1,j)間甲烷透氣性系數(shù)。

        2.3 解算程序編制

        依據(jù)前文所述離散方程以及模型邊界的初始和邊界條件,編制了基于VB語(yǔ)言的解算程序,程序結(jié)構(gòu)流程如圖3所示。

        圖3 程序流程圖

        3 數(shù)值計(jì)算結(jié)果及分析

        數(shù)值模擬以振興二礦2#煤層基本參數(shù)為參數(shù)值,最大吸附瓦斯量a為44.81 m3/t,吸附常數(shù)b為0.522 MPa-1,煤的密度ρm為1.25 t/m3,煤的水分含量W為0.93%,煤的灰分含量Ad為8.17%,掘進(jìn)速度v為3 m/d,煤層原始滲透率K0為6×10-4μm2,煤體孔隙率n為11%,煤層平均厚度h為3 m。

        3.1 掘進(jìn)巷道煤壁瓦斯分布

        經(jīng)計(jì)算,工作面掘進(jìn)10 d、20 d時(shí)巷道圍巖煤層瓦斯壓力分布如圖4所示。隨著工作面勻速推進(jìn),在掘進(jìn)頭前方形成了8 m左右穩(wěn)定的瓦斯低壓區(qū)域,且瓦斯壓力分布基本相同;在不同掘進(jìn)時(shí)間條件下,在工作面掘進(jìn)頭20 m區(qū)域內(nèi)巷幫圍巖煤體瓦斯壓力分布基本相同,驗(yàn)證了隨著工作面向前掘進(jìn),掘進(jìn)工作面一定區(qū)域內(nèi)巷道圍巖煤體瓦斯壓力分布向前移動(dòng),即工作面勻速推進(jìn)中,工作面掘進(jìn)頭一定區(qū)域內(nèi)巷道圍巖煤體瓦斯壓力分布不隨時(shí)間變化,可看作穩(wěn)態(tài)過(guò)程。

        圖4 掘進(jìn)10 d、20 d時(shí)煤層瓦斯壓力分布云圖

        工作面掘進(jìn)22 d時(shí)巷道沿程巷幫瓦斯壓力分布曲線如圖5所示。

        圖5 掘進(jìn)22 d時(shí)沿程巷幫煤體瓦斯壓力分布曲線

        圖5中顯示巷幫內(nèi)向6 m處煤體為瓦斯壓力曲線轉(zhuǎn)折點(diǎn),該點(diǎn)之后的瓦斯壓力梯度較該點(diǎn)之前的瓦斯壓力梯度大??紤]該點(diǎn)之前為卸壓區(qū)與應(yīng)力集中區(qū)內(nèi)的塑性變形區(qū),該區(qū)域內(nèi)煤體滲透率大,致使該區(qū)域內(nèi)的煤體瓦斯瞬間大量涌入掘進(jìn)巷道,但該點(diǎn)處煤體處于應(yīng)力集中區(qū),煤體滲透率較受擾動(dòng)前小,該點(diǎn)后的煤體瓦斯涌出受到制約,難以向揭露煤體方向流動(dòng),因此轉(zhuǎn)折點(diǎn)前的瓦斯壓力梯度減小速率大于轉(zhuǎn)折點(diǎn)之后的瓦斯壓力梯度,表明巷幫煤體滲透率分布特征會(huì)影響煤層內(nèi)瓦斯運(yùn)移過(guò)程。

        3.2 掘進(jìn)巷道瓦斯涌出規(guī)律

        煤巷不同掘進(jìn)長(zhǎng)度下巷道沿程煤壁瓦斯涌出速率分布如圖6所示。由圖6可知,不同掘進(jìn)長(zhǎng)度下巷道沿程煤壁瓦斯涌出速率曲線基本一致,說(shuō)明掘進(jìn)頭一定區(qū)域內(nèi)巷道沿程煤壁瓦斯涌出速率隨著工作面推進(jìn)向前移動(dòng),可認(rèn)為掘進(jìn)頭一定區(qū)域內(nèi)巷道沿程煤壁瓦斯涌出速率不隨時(shí)間變化,看作穩(wěn)態(tài)過(guò)程,且瓦斯涌出速率較大區(qū)域主要集中在掘進(jìn)頭巷道部分,極易造成掘進(jìn)頭瓦斯超限;從圖中可看出隨著掘進(jìn)巷道增長(zhǎng),瓦斯涌出速率總體上逐漸增大。

        圖6 掘進(jìn)煤巷沿程煤壁瓦斯涌出速率分布

        煤巷掘進(jìn)過(guò)程中巷道總體瓦斯涌出速率如圖7所示。

        圖7 掘進(jìn)煤巷瓦斯涌出總速率

        由圖7可以看出,瓦斯涌出速率呈鋸齒狀變化,這是由于瓦斯涌出速率較大區(qū)域主要集中在掘進(jìn)頭巷道部分,每一個(gè)掘進(jìn)循環(huán)中掘進(jìn)頭巷道部分剛被掘開(kāi)時(shí)瓦斯涌出異常劇烈,但隨著時(shí)間延長(zhǎng),瓦斯涌出速率迅速降低,巷道瓦斯涌出總量最大值達(dá)到最小值的3倍左右。

        4 結(jié)論

        (1)隨著工作面向前掘進(jìn),掘進(jìn)工作面一定區(qū)域內(nèi)巷道圍巖煤體瓦斯壓力分布向前移動(dòng),即工作面勻速推進(jìn)中,工作面掘進(jìn)頭一定區(qū)域內(nèi)巷道圍巖煤體瓦斯壓力分布不隨時(shí)間變化,可看作穩(wěn)態(tài)過(guò)程。掘進(jìn)過(guò)程中巷幫瓦斯壓力存在轉(zhuǎn)折點(diǎn),轉(zhuǎn)折點(diǎn)之前的瓦斯壓力梯度減小速率大于轉(zhuǎn)折點(diǎn)之后的瓦斯壓力梯度,說(shuō)明巷幫煤體滲透率分布對(duì)于煤層內(nèi)瓦斯運(yùn)移具有重大影響。

        (2)不同掘進(jìn)長(zhǎng)度下巷道沿程煤壁瓦斯涌出速率曲線基本一致,說(shuō)明掘進(jìn)頭一定區(qū)域內(nèi)巷道沿程煤壁瓦斯涌出速率隨著工作面推進(jìn)向前移動(dòng),且瓦斯涌出速率較大區(qū)域主要集中在掘進(jìn)頭巷道部分。

        (3)掘進(jìn)過(guò)程中瓦斯涌出速率呈鋸齒狀變化,這是由于瓦斯涌出速率較大區(qū)域主要集中在掘進(jìn)頭巷道部分,每一個(gè)掘進(jìn)循環(huán)中掘進(jìn)頭巷道部分剛被掘開(kāi)時(shí)瓦斯涌出異常劇烈。但隨著時(shí)間延長(zhǎng),瓦斯涌出速率迅速降低,巷道瓦斯涌出總量最大值達(dá)到最小值的3倍左右,且由于瓦斯涌出主要集中于掘進(jìn)頭部位,極易造成掘進(jìn)頭瓦斯超限。隨著掘進(jìn)巷道增長(zhǎng),瓦斯涌出速率總體上逐漸增大。

        [1] 候三中,劉德民,李金嶺. 掘進(jìn)工作面煤層滲透率變化對(duì)瓦斯涌出量的影響分析[J].礦業(yè)安全與環(huán)保,2010(5)

        [2] 李宵尖,李炎濤,王兆豐. 煤層巷道周圍瓦斯排放寬度研究[J].煤礦安全,2012(11)

        [3] 劉三鈞,林柏泉,郝志勇等. 鉆孔密封機(jī)理及新型煤層瓦斯壓力測(cè)定技術(shù)研究[J].中國(guó)煤炭,2009(10)

        [4] 許滿貴,林海飛,李樹(shù)剛等. 鉆孔預(yù)抽煤層瓦斯影響規(guī)律研究[J].礦業(yè)安全與環(huán)保,2010(5)

        [5] 方可,吳桂義,馬帥.鉆孔預(yù)抽煤層瓦斯 CFD模擬參數(shù)優(yōu)化研究[J].煤,2012(10)

        [6] 許雙泉.鉆孔預(yù)抽煤層瓦斯數(shù)值模擬及應(yīng)用[D].西安科技大學(xué),2010

        [7] 徐明智, 李希建.煤層瓦斯抽放半徑及其影響因素的數(shù)值模擬[J].礦業(yè)安全與環(huán)保,2012(12)

        [8] 梁冰,劉薊南,孫維吉等.掘進(jìn)工作面瓦斯流動(dòng)規(guī)律數(shù)值模擬分析[J].中國(guó)地質(zhì)災(zāi)害與防治學(xué)報(bào),2011(4)

        [9] 秦躍平,王翠霞,王健等.煤粒瓦斯放散數(shù)學(xué)模型及數(shù)值解算[J].煤炭學(xué)報(bào),2012(9)

        (責(zé)任編輯 張艷華)

        Numerical simulation study of gas dynamic emission laws of driving face

        Yang Tao1, 2, Zhao Dongyun1

        (1. North China Institute of Science & Technology, Sanhe, Hebei 101601, China; 2. State Key Laboratory of Coal Resources and Safe Mining, Haidian, Beijing 100083, China)

        In order to study gas emission laws during driving process, based on C-N discretization method of finite volume method, the authors programmed to numerically simulate gas pressure distribution, coal mass permeability and gas emission rate change around moving driving face. The simulation results indicated that with the heading of driving face, gas pressure distribution of roadway surrounding rock and coal in certain region was moving forward. Drift gas pressure had turning point during driving process and the decreasing rate of gas pressure gradient before turning point was higher than that after turning point. Roadway with different heading length shared the gas emission rate and areas with high gas emission rate were gathered at heading roadway. In the driving process, the gas emission rate was serrate-shaped. Gas emission was extremely intense when every heading roadway was freshly driven. Generally, gas emission rate increased gradually with the driving process.

        driving face, gas pressure, gas emission rate, permeability, numerical simulation

        國(guó)家自然科學(xué)基金項(xiàng)目(51604116),煤炭資源與安全開(kāi)采國(guó)家重點(diǎn)實(shí)驗(yàn)室開(kāi)放基金(SKLCRSM16KFB10),河北省高等學(xué)校科學(xué)技術(shù)研究項(xiàng)目(Z2015104),河北省自然科學(xué)基金資助項(xiàng)目(E2016508036)

        楊濤,趙東云.掘進(jìn)工作面瓦斯動(dòng)態(tài)涌出規(guī)律數(shù)值模擬研究[J].中國(guó)煤炭,2017,43(4):114-118. Yang Tao, Zhao Dongyun . Numerical simulation study of gas dynamic emission laws of driving face [J]. China Coal, 2017,43(4):114-118.

        TD712.52

        A

        楊濤(1983-),男,山東棗莊人,講師,博士,主要從事煤巖瓦斯動(dòng)力災(zāi)害防治方面的研究工作。

        猜你喜歡
        煤壁煤體滲透率
        煤體強(qiáng)度對(duì)煤壁穩(wěn)定性的影響研究
        煤炭工程(2020年8期)2020-08-27 07:27:00
        注熱井周圍煤體蠕變過(guò)程的滲透率變化規(guī)律模擬研究
        正交試驗(yàn)下煤體滲透性影響因素評(píng)價(jià)
        煤炭工程(2019年4期)2019-05-05 02:30:16
        中煤階煤層氣井排采階段劃分及滲透率變化
        不同滲透率巖芯孔徑分布與可動(dòng)流體研究
        SAGD井微壓裂儲(chǔ)層滲透率變化規(guī)律研究
        以“中央廚房”為突破口探索時(shí)政報(bào)道的融煤體之路——以浙江之聲為例
        站在煤壁跟前
        高滲透率風(fēng)電并網(wǎng)對(duì)電力系統(tǒng)失步振蕩的影響
        硬煤工作面煤壁破壞與防治機(jī)理
        永久黄网站色视频免费| 国产亚洲人成在线观看| 三年片免费观看影视大全视频| 亚洲18色成人网站www| 精品中文字幕制服中文| 好看的国内自拍三级网站| 亚洲最大成人综合网720p| 97久久人人超碰超碰窝窝| 欧美日韩国产在线观看免费| 一区二区三区人妻在线| 国产一区二区自拍刺激在线观看| 日韩亚洲欧美中文在线| 亚洲国产成人精品福利在线观看| 免费国产自拍视频在线观看| 麻豆亚洲一区| 成人片黄网站色大片免费观看cn| 日韩在线视精品在亚洲| 久久国产精品亚洲我射av大全| 天天躁夜夜躁av天天爽| 老妇肥熟凸凹丰满刺激| 国产人禽杂交18禁网站| 色综合悠悠88久久久亚洲| 国产成+人欧美+综合在线观看| 国产成人精品午夜福利在线| 国产喷白浆精品一区二区豆腐| 国产激情视频免费在线观看| 国产精品亚韩精品无码a在线| 国产尤物二区三区在线观看| 久久久久亚洲AV无码专区喷| 麻豆三级视频网站在线观看| 久久中文骚妇内射| 女同久久精品国产99国产精品| 成人特黄特色毛片免费看| 91日韩东京热中文字幕| 欧美极品jizzhd欧美| 国产在线欧美日韩一区二区| 伊人精品成人久久综合97| 国产a国产片国产| 欧美日韩国产专区| 久久久噜噜噜噜久久熟女m| 精品综合久久久久久888蜜芽|