• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    鉆孔瓦斯抽采效果影響因素的響應(yīng)面分析

    2022-12-01 08:53:20龔選平韓柏青年福田鄭春山成小雨
    煤礦安全 2022年11期
    關(guān)鍵詞:煤體曲面滲透率

    龔選平,薛 生,韓柏青,年福田,鄭春山,成小雨

    (1.中煤能源研究院有限責任公司,陜西 西安 710054;2.合肥綜合性國家科學中心能源研究院(安徽省能源實驗室),安徽 合肥 230031;3.安徽理工大學 安全科學與工程學院,安徽 淮南 232001;4.中煤新集能源股份有限公司,安徽 淮南 232000;5.中煤新集劉莊礦業(yè)有限公司,安徽 阜陽 236200)

    2021 年煤炭在我國能源消費總量中占比56.0%,與往年相比,煤炭的生產(chǎn)和消費比重雖有所降低,但其在我國能源結(jié)構(gòu)中仍然占據(jù)著主導地位,對我國未來的經(jīng)濟發(fā)展仍將發(fā)揮重要的作用[1-2]。隨著我國煤礦開采深度的增加,地質(zhì)條件日趨復雜,煤礦災害事故的發(fā)生幾率也在增加[3]。與其他煤礦事故相比,瓦斯事故一直是煤礦井下危險程度最大、死亡比例最高的事故類型之一[4],對于瓦斯事故的遏制不容輕視。我國一直堅持“可保盡保、應(yīng)抽盡抽、先抽后采、煤氣共采”的瓦斯治理原則,瓦斯抽采是預防瓦斯事故的重要方法之一,而瓦斯抽采效果受到眾多因素的影響。

    郭欣等[5]以瓦斯?jié)B流理論為基礎(chǔ),利用COMSOL Multiphysics 數(shù)值模擬軟件建立考慮煤巖變形控制方程的氣固耦合數(shù)學模型,研究相關(guān)因素對瓦斯有效抽采半徑大小的影響;徐剛等[6]以流固耦合原理為基礎(chǔ),構(gòu)建煤層瓦斯抽采流固耦合數(shù)學模型,借助FLUENT 數(shù)值模擬軟件對煤層瓦斯預抽鉆孔有效抽采半徑進行數(shù)值模擬研究,對建新煤礦4207 工作面煤層瓦斯預抽鉆孔布孔參數(shù)進行了優(yōu)化;陳月霞等[7]采用三維數(shù)值模擬方法,研究分析了單排不同數(shù)量鉆孔抽采條件下,煤層瓦斯壓力及有效抽采區(qū)域的空間分布特征;張?zhí)燔姷萚8]將相似模擬試驗和數(shù)值模擬相結(jié)合,分析了負壓變化對煤體瓦斯?jié)B流特性的影響;江明泉等[9]運用數(shù)值模擬方法,探討了無斷層影響下順層鉆孔有效抽采半徑變化規(guī)律;李守瑞等[10]提出了理論計算先做預測、數(shù)值模擬驗證理論計算結(jié)果并指導開展現(xiàn)場抽采半徑考察的綜合方法,解決艾維爾溝礦區(qū)瓦斯抽采鉆孔抽采半徑考察精度不高的問題;劉殿平等[11]提出了以抽采鉆孔瓦斯流量負指數(shù)衰減規(guī)律為基礎(chǔ)、以工作面抽采達標指標為判據(jù)的有效抽采半徑測定方法,推導出了以抽采率表示的鉆孔有效抽采半徑計算公式;Kong Xiangguo 等[12]通過構(gòu)建流固耦合滲透率數(shù)值模型,分析水力沖孔條件下多因素對鉆孔有效抽采半徑的影響;Fang Huihuang 等[13]模擬分析多種因素對鉆孔抽采效果以及有效抽采半徑的影響,并對鉆孔抽采條件下的多場多相耦合機理進行了探究。

    國內(nèi)外學者在瓦斯抽采效果的敏感性分析方面取得了豐碩的成果[14-16],但以往研究多偏向于單一因素對瓦斯抽采效果影響的定性分析,而多因素交互作用下的瓦斯抽采效果變化規(guī)律需要進一步的定量分析。為此,以中煤新集公司劉莊礦為研究依托礦井,采用套孔應(yīng)力解除法進行礦區(qū)地應(yīng)力測試,基于煤體基質(zhì)裂隙系統(tǒng)瓦斯解吸-擴散-滲流模型和COMSOLMultiphysics 數(shù)值模擬軟件,分析鉆孔瓦斯抽采過程,并利用Mintab 軟件設(shè)計響應(yīng)曲面試驗,研究多因素交互作用對鉆孔瓦斯抽采效果的影響,找出主次影響因素,確定鉆孔有效抽采半徑與多因素之間的關(guān)系式。研究結(jié)果對煤礦井下瓦斯抽采設(shè)計具有重要的參考價值。

    1 地應(yīng)力測試

    滲透率是影響煤層瓦斯抽采過程的重要因素,而滲透率與地應(yīng)力的大小密切相關(guān),在研究礦區(qū)開展地應(yīng)力實測試驗,為數(shù)值模型構(gòu)建提供地應(yīng)力數(shù)據(jù)。現(xiàn)場測試采用套孔應(yīng)力解除法,為保證測試數(shù)據(jù)準確可靠,在測點處布置了2 個地應(yīng)力測孔K1和K2,測孔均布置在巷道南幫,測孔距離底板約1.5 m,兩孔間距約1.2 m。

    由于地應(yīng)力測量過程中取得的巖心較為破碎,難以通過圍壓試驗獲取用于地應(yīng)力計算的巖石力學性質(zhì)參數(shù),因此將現(xiàn)場從K1 孔和K2 孔取得的能夠用于室內(nèi)巖石力學試驗的巖塊加工成標準巖樣(φ50 mm×100 mm),在中國科學院武漢巖土力學研究所研發(fā)的RMT-150C 巖石力學試驗機上進行單軸壓縮試驗,以獲取巖石的彈性模量、泊松比以及強度參數(shù),K1 孔和K2 孔取得巖心的單軸壓縮試驗結(jié)果見表1。

    表1 K1 孔和K2 孔取得巖心的單軸壓縮試驗結(jié)果Table 1 Uniaxial compression test results of cores obtained from boreholes K1 and K2

    基于套心應(yīng)力解除法的計算公式,以及測量得到的巖石力學參數(shù),計算得到的測點的三維地應(yīng)力測量結(jié)果如圖1。

    圖1 三維地應(yīng)力測量結(jié)果圖Fig.1 3D in-situ stress measurement results

    K1、K2 孔的套孔應(yīng)力解除法地應(yīng)力測試結(jié)果表明:現(xiàn)場測點的最大主應(yīng)力、中間主應(yīng)力的量值區(qū)間分別為18.92~19.80 MPa、14.52~14.57 MPa,最小主應(yīng)力的量值為11.75 MPa;方位角分別為265.9°~266.4°、74.6°~88.6°、176.2~176.5.6°;傾斜角分別為7.0°~11.9°、78.1°~84.6°、-0.43°~-4.44°。水平應(yīng)力大于自重應(yīng)力,最大水平主應(yīng)力與自重應(yīng)力比的范圍為1.38~1.41,說明測點處的應(yīng)力以構(gòu)造應(yīng)力為主。

    2 數(shù)值模型

    2.1 煤體瓦斯運移模型

    為了有效開展數(shù)值模擬研究,在理論上,首先構(gòu)建數(shù)值模型的瓦斯運移控制方程,煤層中的裂隙將煤體分割為一個個基質(zhì)單元體,煤基質(zhì)中的吸附瓦斯在解吸后擴散到裂隙中,然后以滲流的方式在裂隙中流動,即為雙重孔隙介質(zhì)瓦斯運移模型。運用此模型時,需將基質(zhì)滲透率定義為常數(shù),僅考慮裂隙滲透率的動態(tài)變化,將煤層中的瓦斯運移簡化為串聯(lián)過程。

    在煤層瓦斯抽采過程中,煤基質(zhì)單元體內(nèi)吸附瓦斯作為質(zhì)量源向外解吸,使得擴散和滲流持續(xù)進行,煤基質(zhì)與裂隙系統(tǒng)的質(zhì)量交換方程為[17-18]:

    式中:Qs為單位體積煤基質(zhì)同裂隙系統(tǒng)的質(zhì)量交換率,kg/(m3·s);a 為煤基質(zhì)形狀因子,m-2;D 為瓦斯擴散系數(shù),m2/s;cm為煤體基質(zhì)中的瓦斯含量,kg/m3;cf為煤體裂隙中的瓦斯含量,kg/m3;τ 為吸附時間,s。

    根據(jù)理想氣體狀態(tài)方程,煤體基質(zhì)與裂隙中瓦斯含量與壓力之間的關(guān)系如下:

    式中:Mg為氣體分子質(zhì)量,kg/mol;R 為普適氣體常數(shù),其值為8.314 J/(mol·K);T 為煤層溫度,K;pm為煤體基質(zhì)中的瓦斯壓力,MPa;pf為煤體裂隙中的瓦斯壓力,MPa。

    關(guān)于吸附時間τ 的基質(zhì)內(nèi)瓦斯向裂隙擴散方程可變?yōu)椋?/p>

    在煤層瓦斯抽采過程中,基質(zhì)系統(tǒng)是裂隙系統(tǒng)的正質(zhì)量源,裂隙系統(tǒng)是基質(zhì)系統(tǒng)的負質(zhì)量源,固由質(zhì)量守恒定律可知,煤基質(zhì)與裂隙系統(tǒng)的質(zhì)量交換率應(yīng)等于煤基質(zhì)內(nèi)瓦斯質(zhì)量隨時間的變化量,即:

    式中:mm為單位體積煤體基質(zhì)中總的瓦斯賦存質(zhì)量,kg;t 為時間,s。

    瓦斯的吸附作用只發(fā)生在煤基質(zhì)中,煤基質(zhì)中的瓦斯氣體質(zhì)量包含了吸附態(tài)瓦斯和游離態(tài)瓦斯,因此單位體積煤基質(zhì)內(nèi)的瓦斯賦存質(zhì)量可以表示為[19]:

    式中:φm為煤體基質(zhì)孔隙度,%;ρn為標準狀態(tài)下 的 瓦 斯 密 度,kg/m3;ρc為 煤 體 假 密 度,kg/m3;VL為Langmuir 體 積 常 數(shù),m3/t;pL為Langmuir 壓 力 常數(shù),MPa。

    其中瓦斯氣體在標準狀態(tài)時的密度可通過(7)式計算:

    式中:Vm為標準狀態(tài)下的理想氣體摩爾體積,m3/mol。

    整理可得煤基質(zhì)瓦斯壓力隨時間變化的控制方程為:

    由此單位體積煤體裂隙系統(tǒng)內(nèi)瓦斯質(zhì)量平衡方程為:

    式中:mf為單位體積煤體裂隙中的游離瓦斯質(zhì)量,kg;ρf為煤體裂隙系統(tǒng)中的瓦斯密度,kg/m3;q→為Darcy 定律 的速度矢量,m/s;φf為煤 體裂隙孔隙度,%。

    單位體積煤體裂隙中賦存的游離瓦斯質(zhì)量為:

    由于瓦斯氣體質(zhì)量很小,忽略重力對瓦斯在煤層中擴散和流動的影響,根據(jù)Darcy 定律可得到氣體的速度矢量為:

    式中:k 為煤體滲透率,m2;μ 為瓦斯的動力黏度系數(shù),Pa·s。

    最后可得裂隙系統(tǒng)中的瓦斯流動控制方程:

    該模型被用于控制瓦斯抽采數(shù)值模擬中的瓦斯流動過程,保證瓦斯抽采參數(shù)敏感性分析的準確性。

    2.2 幾何模型與邊界條件

    基于劉莊礦的實際地質(zhì)條件,建立數(shù)值模型,該模型的長度為40 m,高度為10 m,鉆孔之間的距離為10 m。邊界條件方面,在模型兩側(cè)設(shè)置輥支撐,模型的底部邊界設(shè)置固定約束,模型整體可以發(fā)生沉降。根據(jù)收集的礦井煤層參數(shù)資料,模型上表面施加均布載荷壓應(yīng)力19.8 MPa,煤層的初始瓦斯壓力為0.6 MPa,煤層的孔隙度為0.01,瓦斯密度為0.716 kg/m3,瓦斯動力黏度為1.8×10-5Pa·s。數(shù)值模型示意圖如圖2。

    圖2 數(shù)值模型示意圖Fig.2 Schematic diagram of numerical model

    3 數(shù)值模擬

    3.1 影響因素篩選

    劉莊礦8 煤層150804 工作面的煤層瓦斯壓力為0.6 MPa,不能以瓦斯壓力下降到0.74 MPa 以下作為有效抽采半徑的指標,根據(jù)相對壓力指標法,如果煤層預抽率為30%,煤層的瓦斯壓力下降量需達到51%[20]。當模擬煤層某點的瓦斯壓力降至0.294 MPa 時,即該點瓦斯壓力的下降量為51%,該點瓦斯抽采達標,則該點至鉆孔處的距離即為有效抽采半徑。

    在多變量正交試驗法中,通常結(jié)合極差分析探究正交試驗結(jié)果的影響,極差分析可直觀地依次得出各因素對試驗數(shù)據(jù)影響的強弱程度。采用正交試驗設(shè)計,開展鉆孔瓦斯抽采參數(shù)對抽采效果的影響程度分析,篩選出影響最不顯著的因素。以鉆孔有效抽采半徑R 為響應(yīng)指標,研究抽采時間、煤層初始滲透率、抽采負壓、孔徑等參數(shù)對鉆孔有效抽采半徑的影響。以上4 種因子均設(shè)置3 個水平,采用3 水平4 因子的正交水平設(shè)計表,正交方案為L9(34),共9 種試驗方案,正交試驗方案見表2。

    表2 正交試驗方案Table 2 Orthogonal test schemes

    對正交試驗的結(jié)果進行處理分析,計算各因子各個水平的均值和極差,并進行歸一化處理,再對各個因子進行排秩,了解各個因子對指標的影響程度。有效抽采半徑極差分析見表3。不同因素不同水平下的有效抽采半徑均值如圖3。

    圖3 不同因素不同水平下的有效抽采半徑均值Fig.3 The mean values of effective extraction radius under different factors and their different levels

    表3 有效抽采半徑極差分析Table 3 Range analysis of effective drainage radius

    由圖3 可以看出,煤層初始滲透率這一因素在3 個水平下的均值差距最大,這表明其影響程度最大;鉆孔孔徑這一因素在3 個水平下的均值差不多,這表明其影響程度很小。另外,結(jié)合表3 可以得出,4 個因素對鉆孔的有效抽采半徑有著不同程度的影響,各因素的影響程度按照從大到小的排序為:煤層初始滲透率、抽采時間、抽采負壓、鉆孔孔徑。其中,煤層初始滲透率這一因素的3 個水平中,鉆孔有效抽采半徑最大均值為2.618 m,最小均值為0.591 m,煤層初始滲透率的有效抽采半徑極差為2.027 m,說明其對鉆孔有效抽采半徑的影響程度最高;鉆孔孔徑這一因素的3 個水平中,鉆孔有效抽采半徑最大均值為1.645 m,最小均值為1.037 m,鉆孔孔徑的有效抽采半徑極差為0.608 m,說明其對鉆孔有效抽采半徑的影響程度最低。因此,在響應(yīng)曲面試驗中剔除鉆孔孔徑這一因素,著重分析煤層初始滲透率、抽采時間和抽采負壓對鉆孔瓦斯抽采效果的影響。

    3.2 響應(yīng)曲面模型

    響應(yīng)曲面法是利用數(shù)學和統(tǒng)計分析進行參數(shù)優(yōu)化的方法,通過擬合響應(yīng)函數(shù)與影響因素之間的一階或者二階模型作為真實響應(yīng)函數(shù)的1 個近似,在多元線性回歸的基礎(chǔ)上主動收集數(shù)據(jù),以獲得具有較好性質(zhì)的回歸方程。建立的復雜多維空間曲面較接近實際情況,所需要的試驗組數(shù)相對較少。響應(yīng)曲面法經(jīng)常使用的設(shè)計方法為中心復合設(shè)計和Box-Behnken 試驗設(shè)計。在因素相同時,Box-Behnken 試驗設(shè)計比中心復合設(shè)計試驗次數(shù)少,而且具有近似旋轉(zhuǎn)性、無序慣性,沒有試驗因素同時為高水平的試驗組合,因此采用Box-Behnken 試驗設(shè)計。Box-Behnken 試驗設(shè)計共15 種試驗方案,運用COMSOL軟件解算不同抽采參數(shù)條件下的鉆孔有效抽采半徑。響應(yīng)曲面試驗設(shè)計方案及結(jié)果見表4。

    根據(jù)結(jié)果可建立有效抽采半徑與多因素的耦合關(guān)系模型,協(xié)助指導抽采參數(shù)動態(tài)優(yōu)化調(diào)整。由于采用大于二階多項式的響應(yīng)曲面模型會使得高次項系數(shù)的數(shù)目增大,大大增加計算量,若采用二次多項式,較為靈活簡單,且擬合精度較高,應(yīng)用廣泛,因此采用二次多項式來表達。針對表4 中得到的模擬結(jié)果,采用響應(yīng)曲面試驗開展多因素回歸擬合分析,建立有效抽采半徑與多因素之間的多項式響應(yīng)曲面回歸方程,如下式:

    表4 響應(yīng)曲面試驗設(shè)計方案及結(jié)果Table 4 Response surface experimental design schemes and results

    式中:r 為有效抽采半徑;k 為煤層滲透率;t 為抽采時間;pQ為抽采負壓。

    3.3 有效抽采半徑響應(yīng)曲面

    二次多項式響應(yīng)面模型的近似函數(shù)可通過變量代換轉(zhuǎn)化為形式上的線性函數(shù),然后通過試驗樣本空間的參數(shù)矩陣得到基于響應(yīng)面模型的函數(shù)值,由此計算響應(yīng)值與實驗值之間的誤差,利用最小二乘法求解上式中二次多項式系數(shù)。方差分析見表5。表中P 值為各項顯著性的關(guān)鍵分析值,如果P 值越小,說明極端的假設(shè)情況發(fā)生的概率越小,則說明結(jié)果越顯著。

    由表5 可以看出,目標函數(shù)R 的響應(yīng)面模型P值遠遠小于0.01,即表示模型顯著性極好,所得多項式回歸方程可準確反映各因素對響應(yīng)值(有效抽采半徑)的影響。模型中的初始滲透率項、抽采時間項、初始滲透率的平方項以及初始滲透率與抽采時間交互項的P 值都小于0.05,說明這幾項在模型中都是顯著的,其他項的P 值都沒遠大于0.05,說明其它項的顯著性也可以接受。其決定系數(shù)為0.9957,表明99.57%以上的響應(yīng)值均可由這個模型解釋。

    表5 方差分析Table 5 Variance analysis

    根據(jù)二次多項回歸方程繪制出響應(yīng)曲面,煤層初始滲透率與抽采時間的響應(yīng)曲面如圖4,煤層初始滲透率與抽采負壓的響應(yīng)曲面如圖5,抽采時間與抽采負壓的響應(yīng)曲面如圖6。

    圖4 煤層初始滲透率與抽采時間的響應(yīng)曲面Fig.4 Response surface of initial permeability of coal seam and extraction time

    圖5 煤層初始滲透率與抽采負壓的響應(yīng)曲面Fig.5 Response surface of initial permeability of coal seam and negative drainage pressure

    圖6 抽采時間與抽采負壓的響應(yīng)曲面Fig.6 Response surface of drainage time and drainage negative pressure

    煤層初始滲透率、抽采時間與抽采負壓兩兩之間的響應(yīng)曲面圖表明,在煤層初始滲透率的坐標軸方向上,鉆孔的有效抽采半徑變化最快,而在抽采負壓的坐標軸方向上,鉆孔有效抽采半徑變化最慢。其中,圖4 的響應(yīng)曲面扭曲最大,說明煤層初始滲透率與抽采時間之間的交互作用顯著;圖5 的響應(yīng)曲面有一些扭曲,說明煤層初始滲透率與抽采負壓之間的交互作用有些許顯著;圖6 的響應(yīng)面沒有明顯的扭曲,說明兩者的交互作用不顯著。這也證實了煤層初始滲透率在鉆孔抽采瓦斯過程中,對鉆孔的有效抽采半徑影響顯著性極大;抽采時間對鉆孔的有效抽采半徑影響顯著性較大;抽采負壓對鉆孔的有效抽采半徑影響顯著性較小。因為影響鉆孔瓦斯抽采最為重要的因素是煤層的滲透率,當煤層滲透率高時,煤層中的裂隙發(fā)育好,透氣性高,瓦斯更容易在煤層中流動,也更容易被鉆孔抽采出來。綜上所述,對鉆孔有效抽采半徑的影響程度為:煤層初始滲透率>抽采時間>抽采負壓,這與正交試驗的極差分析結(jié)果相互吻合。故對于低透氣性煤層,在抽采前采取水力割縫、保護層開采、松動爆破等卸壓增透措施極為重要,可增加低透氣性煤層的滲透率,提高瓦斯抽采效率。

    4 結(jié) 語

    1)根據(jù)正交試驗結(jié)果,煤層初始滲透率對應(yīng)的有效抽采半徑極差最大,可達2.027 m;鉆孔孔徑對應(yīng)的有效抽采半徑極差最小,其值為0.608 m。煤層初始滲透率對鉆孔有效抽采半徑的影響程度最大,抽采時間、抽采負壓、鉆孔孔徑的影響程度依次減小。

    2)采用響應(yīng)面法得到煤層初始滲透率k、抽采時間t、抽采負壓pQ與有效抽采半徑r 的關(guān)系模型為r=-0.093+2.069k-0.004 20t+0.045 4pQ-2.204k2+0.000 020t2-0.000 985pQ2+0.018 02kt+0.006 89kpQ+0.000 008tpQ,該響應(yīng)面模型P 值小于0.05,表明其具有較好的顯著性,模型決定系數(shù)為0.995 7,該模型可以解釋99.57%以上的響應(yīng)值。

    3)多因素交互作用對鉆孔有效抽采半徑變化產(chǎn)生新影響,煤層初始滲透率的影響仍然最顯著,同時,煤層初始滲透率與抽采時間的響應(yīng)曲面扭曲程度最大,二者之間對抽采半徑的交互影響作用也最大;抽采時間與抽采負壓的響應(yīng)曲面扭曲程度最小,說明二者的交互影響作用不顯著。故提高煤體滲透率是瓦斯抽采提效的主要途徑之一。

    猜你喜歡
    煤體曲面滲透率
    注熱井周圍煤體蠕變過程的滲透率變化規(guī)律模擬研究
    正交試驗下煤體滲透性影響因素評價
    煤炭工程(2019年4期)2019-05-05 02:30:16
    相交移動超曲面的亞純映射的唯一性
    中煤階煤層氣井排采階段劃分及滲透率變化
    圓環(huán)上的覆蓋曲面不等式及其應(yīng)用
    不同滲透率巖芯孔徑分布與可動流體研究
    SAGD井微壓裂儲層滲透率變化規(guī)律研究
    以“中央廚房”為突破口探索時政報道的融煤體之路——以浙江之聲為例
    傳媒評論(2017年7期)2017-11-08 03:11:15
    基于曲面展開的自由曲面網(wǎng)格劃分
    高滲透率風電并網(wǎng)對電力系統(tǒng)失步振蕩的影響
    91字幕亚洲| 男女午夜视频在线观看| 天堂动漫精品| 麻豆国产av国片精品| 欧美成人免费av一区二区三区| 免费人成在线观看视频色| 亚洲色图av天堂| 国产精品久久久久久亚洲av鲁大| 亚洲熟妇中文字幕五十中出| 欧美成人一区二区免费高清观看| 欧美极品一区二区三区四区| 亚洲av中文字字幕乱码综合| 我的老师免费观看完整版| 精品国产三级普通话版| 亚洲黑人精品在线| www.熟女人妻精品国产| 亚洲精品一卡2卡三卡4卡5卡| 精华霜和精华液先用哪个| 欧美性猛交╳xxx乱大交人| 18禁黄网站禁片午夜丰满| 怎么达到女性高潮| 国产乱人视频| 99国产综合亚洲精品| 午夜激情欧美在线| netflix在线观看网站| 成人国产一区最新在线观看| 黄色视频,在线免费观看| 女警被强在线播放| 亚洲成av人片在线播放无| 在线观看av片永久免费下载| 国产三级黄色录像| 日本黄色片子视频| 老熟妇乱子伦视频在线观看| www.色视频.com| 久久久久久久精品吃奶| 久久久久精品国产欧美久久久| 国产毛片a区久久久久| 久久欧美精品欧美久久欧美| 欧美性感艳星| 最新在线观看一区二区三区| 婷婷亚洲欧美| 日韩欧美三级三区| 久99久视频精品免费| 可以在线观看毛片的网站| 人妻夜夜爽99麻豆av| 国产高清视频在线观看网站| 日韩精品中文字幕看吧| 美女免费视频网站| 一区二区三区免费毛片| 观看免费一级毛片| 一区二区三区高清视频在线| 亚洲av电影在线进入| 无人区码免费观看不卡| 亚洲成av人片在线播放无| 亚洲精品在线观看二区| 91在线精品国自产拍蜜月 | 欧美日本亚洲视频在线播放| 叶爱在线成人免费视频播放| 97超级碰碰碰精品色视频在线观看| 麻豆一二三区av精品| 麻豆久久精品国产亚洲av| 99精品欧美一区二区三区四区| 久久性视频一级片| 国产精品电影一区二区三区| 婷婷精品国产亚洲av在线| 国产精品爽爽va在线观看网站| 热99re8久久精品国产| 中文字幕久久专区| 免费看a级黄色片| 国产精品爽爽va在线观看网站| 亚洲欧美激情综合另类| 久久精品亚洲精品国产色婷小说| 国内揄拍国产精品人妻在线| 免费看十八禁软件| 国产激情欧美一区二区| 精品电影一区二区在线| 亚洲熟妇熟女久久| 美女高潮的动态| 国产探花极品一区二区| avwww免费| 夜夜爽天天搞| 国产色婷婷99| 人人妻人人澡欧美一区二区| 国产精品女同一区二区软件 | 香蕉丝袜av| 美女高潮的动态| 伊人久久精品亚洲午夜| 香蕉av资源在线| 久久精品91无色码中文字幕| 中文在线观看免费www的网站| 在线a可以看的网站| 一进一出好大好爽视频| 一进一出抽搐动态| 亚洲最大成人中文| 国产高清激情床上av| 日韩高清综合在线| 成人一区二区视频在线观看| 婷婷精品国产亚洲av在线| 少妇裸体淫交视频免费看高清| 老司机福利观看| 亚洲精品一卡2卡三卡4卡5卡| 久久精品国产清高在天天线| 成人av在线播放网站| 成人精品一区二区免费| 午夜a级毛片| 99riav亚洲国产免费| 国产成人av激情在线播放| 国内少妇人妻偷人精品xxx网站| 国产精品日韩av在线免费观看| 免费在线观看亚洲国产| 国产精品98久久久久久宅男小说| 哪里可以看免费的av片| 成人鲁丝片一二三区免费| 亚洲国产精品合色在线| 亚洲,欧美精品.| 好男人在线观看高清免费视频| 国产精品电影一区二区三区| 97人妻精品一区二区三区麻豆| 亚洲国产精品成人综合色| 国产精品亚洲av一区麻豆| 精品久久久久久久久久久久久| 最新美女视频免费是黄的| 国产伦在线观看视频一区| 国内少妇人妻偷人精品xxx网站| 五月玫瑰六月丁香| 欧美丝袜亚洲另类 | 免费人成在线观看视频色| ponron亚洲| 中文字幕熟女人妻在线| 老司机午夜福利在线观看视频| 老司机在亚洲福利影院| www.999成人在线观看| 亚洲最大成人手机在线| 午夜日韩欧美国产| 久久精品夜夜夜夜夜久久蜜豆| 一边摸一边抽搐一进一小说| 可以在线观看毛片的网站| 亚洲精品一区av在线观看| 欧美激情在线99| av黄色大香蕉| 精品电影一区二区在线| 中文字幕熟女人妻在线| 国产在视频线在精品| 1000部很黄的大片| 亚洲av第一区精品v没综合| 日本在线视频免费播放| 97碰自拍视频| 久久精品91无色码中文字幕| 成人精品一区二区免费| 搡女人真爽免费视频火全软件 | 成人国产综合亚洲| 国产成年人精品一区二区| 可以在线观看的亚洲视频| 亚洲一区二区三区不卡视频| 男女做爰动态图高潮gif福利片| 午夜福利视频1000在线观看| 亚洲一区二区三区不卡视频| 国产精品免费一区二区三区在线| 最新在线观看一区二区三区| 又紧又爽又黄一区二区| 国产国拍精品亚洲av在线观看 | 91久久精品电影网| 成人午夜高清在线视频| 成人欧美大片| 国产主播在线观看一区二区| bbb黄色大片| 日韩国内少妇激情av| 国产精品免费一区二区三区在线| 免费在线观看成人毛片| 免费av不卡在线播放| 亚洲av免费在线观看| 国产淫片久久久久久久久 | 黄色日韩在线| 久久久久久大精品| 成熟少妇高潮喷水视频| 国产免费一级a男人的天堂| 叶爱在线成人免费视频播放| 欧美一级a爱片免费观看看| av中文乱码字幕在线| 91久久精品国产一区二区成人 | 国产主播在线观看一区二区| 少妇人妻精品综合一区二区 | 国产一区在线观看成人免费| 免费看十八禁软件| 国产精品99久久久久久久久| 久久精品影院6| 欧美黄色淫秽网站| 九九在线视频观看精品| 免费观看人在逋| 级片在线观看| 国产男靠女视频免费网站| 亚洲精华国产精华精| 国产av麻豆久久久久久久| 国产成人系列免费观看| 性色avwww在线观看| 午夜日韩欧美国产| 国产91精品成人一区二区三区| 国产中年淑女户外野战色| 国产精品98久久久久久宅男小说| 国产精品女同一区二区软件 | 精品久久久久久成人av| 中文字幕av成人在线电影| 亚洲欧美日韩卡通动漫| 午夜免费男女啪啪视频观看 | 夜夜躁狠狠躁天天躁| 村上凉子中文字幕在线| 国产欧美日韩一区二区精品| 少妇熟女aⅴ在线视频| 国产高清视频在线观看网站| 欧洲精品卡2卡3卡4卡5卡区| 国产午夜精品久久久久久一区二区三区 | 亚洲av电影在线进入| 日本在线视频免费播放| 精品久久久久久,| 午夜福利视频1000在线观看| 亚洲色图av天堂| 国产三级黄色录像| 99热这里只有是精品50| 久久精品夜夜夜夜夜久久蜜豆| 在线观看免费午夜福利视频| 午夜两性在线视频| 亚洲av熟女| 国产精品野战在线观看| 国产免费av片在线观看野外av| 黄色片一级片一级黄色片| av福利片在线观看| 精品免费久久久久久久清纯| 一区二区三区国产精品乱码| 国产高清视频在线观看网站| 成人18禁在线播放| 男女午夜视频在线观看| 国产真实乱freesex| 激情在线观看视频在线高清| 欧美性猛交╳xxx乱大交人| 男女那种视频在线观看| 最新中文字幕久久久久| 午夜激情福利司机影院| 国产精品亚洲一级av第二区| 亚洲七黄色美女视频| 一级毛片女人18水好多| 国产高清激情床上av| 久久香蕉精品热| 91av网一区二区| 99在线视频只有这里精品首页| 手机成人av网站| 亚洲无线在线观看| 国产av一区在线观看免费| 国产精品影院久久| 最近在线观看免费完整版| 伊人久久精品亚洲午夜| 欧美性感艳星| 色综合婷婷激情| 宅男免费午夜| 欧美国产日韩亚洲一区| 3wmmmm亚洲av在线观看| 欧美色欧美亚洲另类二区| 国产老妇女一区| 男女之事视频高清在线观看| 韩国av一区二区三区四区| 97超视频在线观看视频| 尤物成人国产欧美一区二区三区| av国产免费在线观看| 99精品久久久久人妻精品| 国语自产精品视频在线第100页| 国产一区二区在线av高清观看| 蜜桃久久精品国产亚洲av| 男人舔奶头视频| 男插女下体视频免费在线播放| 欧美精品啪啪一区二区三区| 淫秽高清视频在线观看| 国产欧美日韩精品一区二区| 岛国视频午夜一区免费看| 又爽又黄无遮挡网站| 午夜福利高清视频| 精品欧美国产一区二区三| 少妇高潮的动态图| 免费看美女性在线毛片视频| 亚洲国产欧美网| 国产精品一区二区三区四区久久| 国产精品久久电影中文字幕| 亚洲精品亚洲一区二区| 成人性生交大片免费视频hd| 午夜两性在线视频| 亚洲精华国产精华精| 搡老熟女国产l中国老女人| 长腿黑丝高跟| 日韩大尺度精品在线看网址| 精品乱码久久久久久99久播| 亚洲不卡免费看| 精品久久久久久,| 国产成人aa在线观看| 亚洲国产高清在线一区二区三| 两个人视频免费观看高清| 欧美日本视频| 嫩草影院精品99| 免费观看精品视频网站| 他把我摸到了高潮在线观看| 人妻久久中文字幕网| 人人妻,人人澡人人爽秒播| 99在线视频只有这里精品首页| 国产老妇女一区| 岛国在线免费视频观看| 亚洲成av人片在线播放无| 床上黄色一级片| 男女那种视频在线观看| 97超视频在线观看视频| 观看免费一级毛片| 国产伦精品一区二区三区视频9 | 欧美又色又爽又黄视频| bbb黄色大片| av在线蜜桃| 日韩高清综合在线| 在线天堂最新版资源| 丝袜美腿在线中文| 国产伦在线观看视频一区| 亚洲无线在线观看| 在线国产一区二区在线| 久久久精品欧美日韩精品| 免费一级毛片在线播放高清视频| 九色成人免费人妻av| 国产精品一区二区免费欧美| 99精品欧美一区二区三区四区| 久久久久精品国产欧美久久久| 国产不卡一卡二| 久久草成人影院| 久久人人精品亚洲av| 亚洲成人久久爱视频| 精品国产超薄肉色丝袜足j| 午夜福利在线观看吧| 欧美不卡视频在线免费观看| 欧美激情在线99| 日本熟妇午夜| 丝袜美腿在线中文| 网址你懂的国产日韩在线| 岛国在线免费视频观看| 激情在线观看视频在线高清| 男女视频在线观看网站免费| 欧美日韩精品网址| 亚洲国产精品成人综合色| 少妇熟女aⅴ在线视频| 制服丝袜大香蕉在线| 亚洲成av人片在线播放无| 亚洲电影在线观看av| 91久久精品国产一区二区成人 | 90打野战视频偷拍视频| 在线播放无遮挡| 欧美日韩中文字幕国产精品一区二区三区| a级一级毛片免费在线观看| 他把我摸到了高潮在线观看| 制服人妻中文乱码| 久久中文看片网| 欧美大码av| 久久久久久久久大av| 高清在线国产一区| 热99在线观看视频| 丰满人妻熟妇乱又伦精品不卡| avwww免费| 欧美日韩亚洲国产一区二区在线观看| 国产主播在线观看一区二区| 99久久成人亚洲精品观看| 国产主播在线观看一区二区| 日本在线视频免费播放| 国产私拍福利视频在线观看| 亚洲av成人精品一区久久| 国产精品免费一区二区三区在线| 久久久精品大字幕| 88av欧美| 国产av麻豆久久久久久久| 97超视频在线观看视频| 久久精品国产99精品国产亚洲性色| 亚洲电影在线观看av| 国产精品日韩av在线免费观看| 亚洲18禁久久av| 久久精品国产亚洲av香蕉五月| 夜夜夜夜夜久久久久| 成人高潮视频无遮挡免费网站| 19禁男女啪啪无遮挡网站| 欧美国产日韩亚洲一区| 波多野结衣高清无吗| 日韩欧美精品免费久久 | 99热这里只有是精品50| 国产主播在线观看一区二区| 99久久综合精品五月天人人| 村上凉子中文字幕在线| 少妇人妻精品综合一区二区 | 亚洲成人久久爱视频| 久久精品91蜜桃| 最近最新免费中文字幕在线| 国产精品乱码一区二三区的特点| www.999成人在线观看| 蜜桃久久精品国产亚洲av| 亚洲国产精品sss在线观看| 亚洲国产色片| 一区二区三区激情视频| 俄罗斯特黄特色一大片| 高潮久久久久久久久久久不卡| 欧美又色又爽又黄视频| 99在线人妻在线中文字幕| 国产成人啪精品午夜网站| 麻豆成人av在线观看| 久久久久久久久大av| 欧美丝袜亚洲另类 | 国产精品嫩草影院av在线观看 | 舔av片在线| 内地一区二区视频在线| 一个人看的www免费观看视频| 国产精品一区二区三区四区久久| 亚洲精品456在线播放app | 精品久久久久久成人av| 久久久久久人人人人人| 免费观看人在逋| 国产欧美日韩精品一区二区| 久久精品国产亚洲av香蕉五月| 成人亚洲精品av一区二区| 国产伦精品一区二区三区视频9 | 国产免费av片在线观看野外av| 欧美大码av| 老汉色av国产亚洲站长工具| 天堂av国产一区二区熟女人妻| 国产精品自产拍在线观看55亚洲| 久久精品国产亚洲av涩爱 | 国产又黄又爽又无遮挡在线| 亚洲内射少妇av| 小蜜桃在线观看免费完整版高清| 搡老妇女老女人老熟妇| 日本成人三级电影网站| 一本一本综合久久| 午夜免费成人在线视频| 国产v大片淫在线免费观看| 热99在线观看视频| 亚洲欧美日韩高清在线视频| 精品久久久久久久毛片微露脸| 亚洲一区二区三区不卡视频| 欧美成人a在线观看| 欧美日韩综合久久久久久 | 国产精品,欧美在线| 亚洲在线自拍视频| 亚洲av成人av| 在线播放无遮挡| 狠狠狠狠99中文字幕| 欧美日韩精品网址| 日韩亚洲欧美综合| 国产精品影院久久| 色老头精品视频在线观看| 18禁美女被吸乳视频| 国产精品98久久久久久宅男小说| 久久久色成人| 亚洲一区二区三区色噜噜| 操出白浆在线播放| 日本五十路高清| 国产精品亚洲av一区麻豆| 久久精品夜夜夜夜夜久久蜜豆| 久久人妻av系列| 国产三级黄色录像| 国产精品乱码一区二三区的特点| 国产成人影院久久av| 老熟妇仑乱视频hdxx| 国产亚洲欧美98| 一a级毛片在线观看| 午夜久久久久精精品| 天堂√8在线中文| 在线观看舔阴道视频| 免费人成视频x8x8入口观看| 亚洲电影在线观看av| 无人区码免费观看不卡| 欧美在线黄色| 久久久久久久亚洲中文字幕 | 亚洲成人久久爱视频| 精品久久久久久,| 高清在线国产一区| 久久精品综合一区二区三区| 欧美激情久久久久久爽电影| 久久午夜亚洲精品久久| 亚洲自拍偷在线| 日本三级黄在线观看| 亚洲熟妇熟女久久| 日韩国内少妇激情av| a在线观看视频网站| 99久久久亚洲精品蜜臀av| 欧美3d第一页| 亚洲激情在线av| 怎么达到女性高潮| 性色av乱码一区二区三区2| 国产精品亚洲一级av第二区| 99热6这里只有精品| 最新中文字幕久久久久| 久久久精品大字幕| 深爱激情五月婷婷| 国产视频一区二区在线看| 狂野欧美白嫩少妇大欣赏| 18禁黄网站禁片免费观看直播| 天堂影院成人在线观看| 国产精品爽爽va在线观看网站| 成人国产综合亚洲| 别揉我奶头~嗯~啊~动态视频| 男女午夜视频在线观看| 中文字幕人妻熟人妻熟丝袜美 | 757午夜福利合集在线观看| 美女cb高潮喷水在线观看| 男人舔奶头视频| 欧美极品一区二区三区四区| 久久精品影院6| 99精品久久久久人妻精品| 少妇人妻一区二区三区视频| 久久国产精品人妻蜜桃| 最近视频中文字幕2019在线8| 国产黄片美女视频| 操出白浆在线播放| eeuss影院久久| 亚洲久久久久久中文字幕| 丰满人妻一区二区三区视频av | 天堂动漫精品| 韩国av一区二区三区四区| 成人无遮挡网站| 男人舔奶头视频| 国产精品 欧美亚洲| 久久精品影院6| 欧美乱色亚洲激情| 欧美乱妇无乱码| 最好的美女福利视频网| 观看免费一级毛片| 亚洲性夜色夜夜综合| 一级黄色大片毛片| 精品熟女少妇八av免费久了| 国产真实伦视频高清在线观看 | 国产精品久久电影中文字幕| 国产精品三级大全| 手机成人av网站| 麻豆成人午夜福利视频| 欧美在线黄色| 一级毛片女人18水好多| 欧美zozozo另类| 一区二区三区高清视频在线| 亚洲人成网站在线播| 变态另类丝袜制服| 国产乱人伦免费视频| 国产精品美女特级片免费视频播放器| eeuss影院久久| 国产激情欧美一区二区| 夜夜爽天天搞| 午夜福利在线在线| 观看免费一级毛片| 男女之事视频高清在线观看| 色哟哟哟哟哟哟| 狂野欧美激情性xxxx| 亚洲精品美女久久久久99蜜臀| 特级一级黄色大片| 麻豆国产av国片精品| 国产不卡一卡二| 亚洲av中文字字幕乱码综合| 一区二区三区免费毛片| 最好的美女福利视频网| 久久久精品大字幕| 久久6这里有精品| ponron亚洲| 丰满人妻熟妇乱又伦精品不卡| 国产一区二区三区在线臀色熟女| 在线天堂最新版资源| 国产单亲对白刺激| 变态另类丝袜制服| 黄色日韩在线| 夜夜夜夜夜久久久久| 88av欧美| 国产精品国产高清国产av| 18禁黄网站禁片午夜丰满| www.色视频.com| 好看av亚洲va欧美ⅴa在| 亚洲成av人片免费观看| www日本黄色视频网| 美女被艹到高潮喷水动态| 精品久久久久久久末码| 日本五十路高清| 亚洲在线观看片| 成人午夜高清在线视频| 国产三级中文精品| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 国内毛片毛片毛片毛片毛片| 日本黄大片高清| 99热这里只有精品一区| 午夜福利免费观看在线| 夜夜看夜夜爽夜夜摸| 成年女人毛片免费观看观看9| 国产av不卡久久| 黄色成人免费大全| 首页视频小说图片口味搜索| АⅤ资源中文在线天堂| 网址你懂的国产日韩在线| 国产日本99.免费观看| 日本精品一区二区三区蜜桃| 手机成人av网站| aaaaa片日本免费| 99热只有精品国产| 婷婷亚洲欧美| 国产视频内射| 亚洲熟妇中文字幕五十中出| 一个人看的www免费观看视频| e午夜精品久久久久久久| 天堂动漫精品| 黄色日韩在线| 国产成人欧美在线观看| 国产亚洲欧美在线一区二区| 中文亚洲av片在线观看爽| 国产精品久久久久久久久免 | 夜夜爽天天搞| 人人妻人人澡欧美一区二区| 脱女人内裤的视频| 亚洲精品久久国产高清桃花| 97人妻精品一区二区三区麻豆| 亚洲av中文字字幕乱码综合| 日本一二三区视频观看| 少妇丰满av| 国产亚洲精品av在线| av女优亚洲男人天堂| 日韩欧美一区二区三区在线观看| 搡老熟女国产l中国老女人|