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

    基于灰色理論對(duì)臨渙煤礦六采區(qū)煤層瓦斯壓力的預(yù)測優(yōu)化

    2020-02-03 02:37:14宣嘯羅文柯李賀黃妍陸俊翔王成龍
    礦業(yè)工程研究 2020年4期
    關(guān)鍵詞:采區(qū)殘差灰色

    宣嘯,羅文柯,李賀,黃妍,陸俊翔,王成龍

    (1.淮北礦業(yè)集團(tuán) 渦北煤礦,安徽 亳州 233611;2.湖南科技大學(xué) 資源環(huán)境與安全工程學(xué)院,湖南 湘潭 411201)

    煤礦開采突出煤層是指在礦井井田范圍內(nèi)發(fā)生過突出或經(jīng)鑒定、認(rèn)定有突出危險(xiǎn)性的煤層[1],對(duì)于突出煤層的瓦斯有效治理方法是在煤的整個(gè)采掘過程中嚴(yán)格實(shí)行“兩個(gè)四位一體”[2],無論是“區(qū)域四位一體”還是“局部四位一體”,突出危險(xiǎn)性預(yù)測是其首選任務(wù),也是提出防治突出危險(xiǎn)性措施和效果檢驗(yàn)的基礎(chǔ)[3].在區(qū)域突出危險(xiǎn)性預(yù)測指標(biāo)中,最為常用的方法是實(shí)測各突出煤層的瓦斯壓力(含量)[4].煤層瓦斯壓力是指瓦斯在煤層中所呈現(xiàn)的氣體壓力,是煤層孔隙和裂隙中的游離瓦斯自由熱運(yùn)動(dòng)對(duì)孔隙和裂隙空間壁面所產(chǎn)生的作用力[5],測定煤層瓦斯壓力工作程序復(fù)雜、工作量繁重,且測定結(jié)果誤差較大.

    灰色系統(tǒng)理論是中國學(xué)者鄧聚龍教授于1982年創(chuàng)立的新理論[6],即以“部分信息已知,部分信息未知”的 “小樣本” “貧信息”或“不確定性”系統(tǒng)為研究對(duì)象,通過對(duì)“部分”已知信息的生成、開發(fā),提出有價(jià)值的信息,實(shí)現(xiàn)對(duì)系統(tǒng)行為、演化規(guī)律的正確描述和有效監(jiān)控[7].經(jīng)過近30多年的改進(jìn)與發(fā)展,該理論以其橫斷面廣、內(nèi)涵豐富、滲透力強(qiáng)、方法實(shí)用、理論與實(shí)踐并重而不斷應(yīng)用到農(nóng)業(yè)、現(xiàn)代化工業(yè)、礦山、地質(zhì)、氣象、交通運(yùn)輸和經(jīng)濟(jì)決策等諸多領(lǐng)域[8].在灰色理論中,從簡單的灰色GM(1,1)出發(fā),經(jīng)不斷研究與改進(jìn),形成了單一類灰色預(yù)測模型和組合類灰色預(yù)測模型兩大類型[9].這些灰色系列預(yù)測模型,因其建模簡單、預(yù)測精度高、數(shù)據(jù)可靠而被廣泛應(yīng)用于包括災(zāi)害預(yù)測在內(nèi)的各工程領(lǐng)域的預(yù)測實(shí)踐中[10].

    本文將非線性科學(xué)技術(shù)引入到煤礦井下煤層瓦斯壓力測定的工作中,特以淮北礦業(yè)集團(tuán)臨渙煤礦六采區(qū)準(zhǔn)備巷道的7煤、9(8)煤和10煤在傾斜方向的實(shí)測煤層瓦斯壓力原始數(shù)據(jù)作為灰色預(yù)測模型的原始數(shù)據(jù)序列,利用灰色GM(1,1)[11]、新陳代謝GM(1,1)模型[12]、GM(1,1)殘差模型[13]和“灰色-馬爾科夫”組合模型[14]進(jìn)行預(yù)測,對(duì)其結(jié)果進(jìn)行精度等級(jí)對(duì)比分析,以獲取某種預(yù)測模型用于預(yù)測煤層瓦斯壓力為最優(yōu).

    1 工程背景

    淮北礦業(yè)集團(tuán)臨渙煤礦井田位于淮北平原,地勢平坦,海拔標(biāo)高一般為+27 m左右,西北高,東南低.礦井東西長13 km,南北寬4~5 km,面積約為50 km2.六采區(qū)位于礦井西翼(含一、二水平),二疊系含煤地層,主采7煤、9(8)煤、10煤.東及東北以DF232斷層上盤斷煤交線為界,南以F3和大吳家斷層下盤斷煤交線為界,西南以各主采煤層露頭為界,西、西北、北以騎路周斷層上盤斷煤交線為界.采區(qū)走向長3 145 m,傾向長1 180 m,面積3.71 km2;7煤、9(8)煤、10煤層回采上限標(biāo)高為-275 m,回采下限標(biāo)高分別為-560,-575,-650 m.采區(qū)平面投影如圖1所示.

    圖1 淮北礦業(yè)集團(tuán)臨渙煤礦六采區(qū)平面投影

    2 灰色預(yù)測模型預(yù)測精度等級(jí)

    預(yù)測目的是借助系統(tǒng)的歷史數(shù)據(jù)及狀態(tài)來推測和了解系統(tǒng)未來的發(fā)展趨勢,以掌握和發(fā)現(xiàn)系統(tǒng)發(fā)展規(guī)律,對(duì)系統(tǒng)的未來狀態(tài)作出科學(xué)的定量預(yù)測.任何預(yù)測模型能否可用于實(shí)際生產(chǎn)指導(dǎo),關(guān)鍵在于經(jīng)過多種精度檢驗(yàn)方法檢驗(yàn)后,來判斷預(yù)測值是否合理和有效.常用預(yù)測模型精度檢驗(yàn)方法有殘差合格檢驗(yàn)法、關(guān)聯(lián)度合格檢驗(yàn)法、均方差比合格檢驗(yàn)法和小誤差概率合格檢驗(yàn)法,如表1所示[15].

    表1 灰色預(yù)測模型精度等級(jí)

    3 灰色預(yù)測模型建立

    3.1 GM(1,1)模型建模

    GM(1,1)模型的原始形式為x(0)(k)+ax(1)(k)=b.

    (1)

    GM(1,1)模型的基本形式為x(0)(k)+az(1)(k)=b.

    (2)

    (3)

    (4)

    通過式(3)和式(4)求解得傳統(tǒng)GM(1,1)模型的時(shí)間響應(yīng)序列和其還原值公式分別為

    (5)

    (6)

    3.2 殘差GM(1,1)模型建模

    當(dāng)采用GM(1,1)模型的各種形式進(jìn)行模擬精度均達(dá)不到要求時(shí),可以考慮對(duì)殘差序列建立GM(1,1)模型,對(duì)原來的模型進(jìn)行修正,以提高模擬精度.

    設(shè)X(0)為原始序列,X(1)為X(0)的1-AGO序列,GM(1,1)模型的時(shí)間響應(yīng)式為

    (7)

    殘差序列為

    (8)

    若存在k0,滿足n-k0≥4,則稱(|ε(0)(k0)|,|ε(0)(k0+1)|,…,|ε(0)(n)|)為可建模殘差尾段,仍記為ε(0)=(ε(0)(k0),ε(0)(k0+1),…,ε(0)(n)).

    殘差修正也就是取k0=2,對(duì)GM(1,1)模型預(yù)測得到的殘差序列的絕對(duì)值重新按照GM(1,1)模型的計(jì)算步驟得到殘差序列的發(fā)展系數(shù)a和灰色作用量b,結(jié)合GM(1,1)模型的時(shí)間響應(yīng)式可以得到相應(yīng)的殘差修正時(shí)間響應(yīng)式:

    (9)

    3.3 新陳代謝GM(1,1)模型建模

    設(shè)原始數(shù)據(jù)序列為X(0)=(x(0)(1),x(0)(2),…,x(0)(n)),置入最新信息x(0)(n+1),去掉原信息x(0)(1),利用X(0)=(x(0)(2),x(0)(3),…,x(0)(n),x(0)(n+1))建立的模型為新陳代謝GM(1,1)模型.

    3.4 灰色狀態(tài)馬爾夫科模型建模

    3.4.1 灰色轉(zhuǎn)移馬爾科夫模型

    設(shè){Xn,n∈T}為隨機(jī)過程,若對(duì)于任意整數(shù)n∈T和任意狀態(tài)i0,i1,…,in+1∈I,若條件概率滿足:

    P(Xn+1=in+1|X0=i0,X1=i1,…,Xn=in)=P(Xn+1=in+1|Xn=in).

    (10)

    則稱{Xn,n∈T}為馬爾柯夫鏈.表示系統(tǒng)未來(t=n+1)所處的狀態(tài)僅與其現(xiàn)在(t=n)所處的狀態(tài)有關(guān),而與其過去(t≤n-1)所處的狀態(tài)無關(guān).若對(duì)任意整數(shù)n∈T和狀態(tài)i,j∈I,則稱

    pij(n)=P(Xn+1=j|Xn=i).

    (11)

    為馬爾柯夫鏈的轉(zhuǎn)移概率.設(shè)pij為轉(zhuǎn)移概率,則稱

    (12)

    為系統(tǒng)狀態(tài)轉(zhuǎn)移概率矩陣.

    3.4.2 灰色狀態(tài)馬爾科夫模型

    (13)

    (14)

    為狀態(tài)轉(zhuǎn)移概率.在實(shí)際中,一般只要考察一步轉(zhuǎn)移概率矩陣P,設(shè)預(yù)測對(duì)象處于⊕k狀態(tài),則考察P中第k行,若max{pk1,pk2,…,pkj}=pkl,則可認(rèn)為下一時(shí)刻系統(tǒng)最有可能由⊕k狀態(tài)轉(zhuǎn)向⊕l狀態(tài).若遇到矩陣P中第k行有2個(gè)或2個(gè)以上概率相同或相近時(shí),則狀態(tài)的未來轉(zhuǎn)向難確定.此時(shí),需要考察2步或s步轉(zhuǎn)移概率矩陣P(2)或P(s),其中s≥3.

    4 六采區(qū)7煤、9(8)煤和10煤瓦斯壓力灰色模型預(yù)測

    為體現(xiàn)預(yù)測模型的科學(xué)合理性,選取六采區(qū)7煤、9(8)煤和10煤在傾斜方向的實(shí)測瓦斯壓力前10個(gè)數(shù)據(jù),詳見表2所示.

    表2 臨渙六采區(qū)7煤、9(8)煤、10煤瓦斯壓力實(shí)測前10個(gè)觀測數(shù)據(jù) MPa

    4.1 GM(1,1)模型預(yù)測

    選取六采區(qū)7煤傾斜方向的實(shí)測瓦斯壓力前10個(gè)數(shù)據(jù)進(jìn)行模擬.

    1)初始化建模原始序列 :{0.400,0.420,0.350,0.320,0.300,0.400,0.300,0.480,0.460,0.420}.

    2)原始序列的1-AGO生成:

    {0.400 0,0.820 0,1.170 0,1.490 0,1.790 0,2.190 0,2.490 0,2.970 0,3.430 0,3.850 0}.

    3)1-AGO生成序列的緊鄰均值生成:

    {0.610 0,0.995 0,1.330 0,1.640 0,1.990 0,2.340 0,2.730 0,3.200 0,3.640 0}.

    5)由此得出7煤瓦斯壓力GM(1,1)模型預(yù)測結(jié)果.同上步驟可得9(8)煤和10煤瓦斯壓力GM(1,1)模型預(yù)測結(jié)果,如表3所示.

    表3 六采區(qū)7煤、9(8)煤和10煤瓦斯壓力傾斜方向GM(1,1)預(yù)測結(jié)果

    6)由此得到,對(duì)六采區(qū)主采的7煤瓦斯壓力預(yù)測結(jié)果平均相對(duì)誤差為13.168%,由誤差檢驗(yàn)計(jì)算得c=0.780,p=0.640,預(yù)測精度等級(jí)為Ⅳ級(jí).

    7)同理,對(duì)六采區(qū)主采的9(8)煤瓦斯壓力預(yù)測結(jié)果的平均相對(duì)誤差為10.591%,由誤差檢驗(yàn)計(jì)算得c=0.714,p=0.616,預(yù)測精度等級(jí)為Ⅳ級(jí).

    8)同理,對(duì)六采區(qū)主采的10煤瓦斯壓力預(yù)測結(jié)果的平均相對(duì)誤差為13.168%.由誤差檢驗(yàn)計(jì)算得c=0.887,p=0.691,預(yù)測精度為Ⅳ級(jí).

    4.2 新陳代謝GM(1,1)模型預(yù)測

    選取六采區(qū)7煤傾斜方向的實(shí)測瓦斯壓力前10個(gè)數(shù)據(jù)進(jìn)行模擬.

    1)根據(jù)新陳代謝GM(1,1)模型灰色建模步驟,先去掉一個(gè)舊序號(hào)1的信息0.400,再置入一個(gè)新實(shí)測數(shù)據(jù)0.400,分別建立原始數(shù)據(jù)序列.

    2)用灰色理論的建模原理進(jìn)行預(yù)測,依此類推,得到7煤瓦斯壓力新陳代謝GM(1,1)模擬預(yù)測結(jié)果如表4所示.

    同上步驟,亦可分別獲得9(8)煤和10煤瓦斯壓力新陳代謝GM(1,1)模型預(yù)測結(jié)果,如表4所示.

    表4 六采區(qū)7煤、9(8)煤和10煤傾斜方向瓦斯壓力新陳代謝GM(1,1)模型預(yù)測結(jié)果

    3)由此得到,對(duì)六采區(qū)主采的7煤瓦斯壓力預(yù)測結(jié)果的平均相對(duì)誤差為0.048 5%,由誤差檢驗(yàn)計(jì)算得c=0.483,p=0.827,預(yù)測精度等級(jí)為Ⅱ級(jí).

    4)同理,對(duì)六采區(qū)主采的9(8)煤瓦斯壓力預(yù)測結(jié)果的平均相對(duì)誤差為0.051 9%,由誤差檢驗(yàn)計(jì)算得c=0.573,p=0.700,預(yù)測精度等級(jí)為Ⅲ級(jí).

    5)同理,對(duì)六采區(qū)主采的10煤瓦斯壓力預(yù)測結(jié)果的平均相對(duì)誤差為0.047 1%,由誤差檢驗(yàn)計(jì)算得c=0.413,p=0.600,預(yù)測精度等級(jí)為Ⅲ級(jí).

    4.3 GM(1,1)殘差模型預(yù)測

    通過GM(1,1)模型預(yù)測得到的六采區(qū)7煤傾斜方向的實(shí)測瓦斯壓力前10個(gè)數(shù)據(jù)的殘差序列,將其絕對(duì)值重新按照GM(1,1)模型的步驟計(jì)算,得到殘差序列的發(fā)展系數(shù)a和灰色作用量b.

    1)初始化建模原始序列:{0.081,0.001,0.040,0.071,0.018,0.094,0.074,0.042,0.011}.

    2)原始序列的1-AGO生成:{0.081 0,0.082 0,0.122 0,0.193 0,0.211 0,0.305 0,0.379 0,0.421 0,0.432 0}.

    3)1-AGO生成序列的緊鄰均值生成:{0.081 5,0.102 0,0.157 5,0.202 0,0.258 0,0.342 0,0.400 0,0.426 5}.

    5)利用GM(1,1)殘差模型預(yù)測六采區(qū)7煤傾斜方向的瓦斯壓力結(jié)果如表5所示.

    同上步驟,亦可分別獲得9(8)煤和10煤瓦斯壓力GM(1,1)殘差模型預(yù)測結(jié)果,如表5所示.

    表5 六采區(qū)7煤、9(8)煤和10煤傾斜方向瓦斯壓力GM(1,1)殘差模型預(yù)測結(jié)果

    6)由此得到,對(duì)六采區(qū)主采的7煤瓦斯壓力預(yù)測結(jié)果的平均模擬相對(duì)誤差為0.071%,由誤差檢驗(yàn)計(jì)算得c=0.327,p=0.750,此模型預(yù)測精度為Ⅲ級(jí).

    7)同理,對(duì)六采區(qū)主采的9(8)煤瓦斯壓力預(yù)測結(jié)果的平均模擬相對(duì)誤差為0.054%,由誤差檢驗(yàn)計(jì)算得c=0.423,p=0.800,此模型預(yù)測精度為Ⅲ級(jí).

    8)同理,對(duì)六采區(qū)主采的10煤瓦斯壓力預(yù)測結(jié)果的平均模擬相對(duì)誤差為0.026%,由誤差檢驗(yàn)計(jì)算得c=0.537,p=0.800,此模型預(yù)測精度等級(jí)為Ⅱ級(jí).

    4.4 GM(1,1)-Markov模型預(yù)測

    選取六采區(qū)7煤傾斜方向的實(shí)測瓦斯壓力前10個(gè)數(shù)據(jù)進(jìn)行模擬.

    1)7煤傾斜方向煤層瓦斯壓力的灰色狀態(tài)劃分如表6所示.

    表6 六采區(qū)7煤傾斜方向煤層瓦斯壓力的灰色狀態(tài)劃分

    2)按照灰色GM(1,1)模型7煤的實(shí)際值和預(yù)測值,根據(jù)其相對(duì)值大小進(jìn)行煤層瓦斯壓力的狀態(tài)預(yù)測,如表7所示.

    表7 六采區(qū)7煤傾斜方向瓦斯壓力狀態(tài)預(yù)測

    經(jīng)計(jì)算得到狀態(tài)轉(zhuǎn)移概率矩陣:

    3)為了得到7煤層瓦斯壓力預(yù)測值,選取鄰近的4個(gè)序號(hào)(6,7,8,9),轉(zhuǎn)移步數(shù)分別定為1,2,3,4,在轉(zhuǎn)移步數(shù)對(duì)應(yīng)的轉(zhuǎn)移矩陣中,取起始狀態(tài)對(duì)應(yīng)的行向量組成新的概率矩陣,再對(duì)新的概率矩陣列向量求和,其和最大的列所對(duì)應(yīng)的狀態(tài)即為序號(hào)10的煤層瓦斯壓力所處的狀態(tài).計(jì)算結(jié)果如表8所示.

    表8 六采區(qū)7煤序號(hào)10的煤層瓦斯壓力所處狀態(tài)預(yù)測

    4)由列向量求和結(jié)果,10號(hào)預(yù)測值所處狀態(tài)為E3,即位于1.000~1.120,預(yù)測值為

    依此,利用GM(1,1)-Markov模型預(yù)測六采區(qū)7煤傾斜方向瓦斯壓力,如表9所示.

    同理利用GM(1,1)-Markov模型預(yù)測六采區(qū)9(8)煤和10煤傾斜方向瓦斯壓力,如表9所示.

    表9 六采區(qū)7煤、9(8)煤和10煤傾斜方向瓦斯壓力GM(1,1)-Markov模型預(yù)測結(jié)果

    5)由此得到,對(duì)六采區(qū)主采的7煤瓦斯壓力預(yù)測結(jié)果的平均模擬相對(duì)誤差為0.047%,由誤差檢驗(yàn)計(jì)算得c=0.493,p=1.000,此模型預(yù)測精度為Ⅱ級(jí).

    6)同理,對(duì)六采區(qū)主采的9(8)煤瓦斯壓力預(yù)測結(jié)果的平均模擬相對(duì)誤差為0.079%,由誤差檢驗(yàn)計(jì)算得c=0.437,p=0.800,此模型預(yù)測精度為Ⅱ級(jí).

    7)同理,對(duì)六采區(qū)主采的10煤瓦斯壓力預(yù)測結(jié)果的平均模擬相對(duì)誤差為0.022%,由誤差檢驗(yàn)計(jì)算得c=0.324,p=0.950,此模型預(yù)測精度為Ⅰ級(jí).

    5 結(jié)論

    1)GM(1,1)模型的預(yù)測精度均為Ⅳ級(jí),預(yù)測數(shù)據(jù)可信度較低,不宜用于井下煤層瓦斯壓力的預(yù)測.

    2)新陳代謝GM(1,1)模型和GM(1,1)殘差模型的預(yù)測精度為Ⅱ級(jí)、Ⅲ級(jí),預(yù)測數(shù)據(jù)可信度較高,可用于井下煤層瓦斯壓力的預(yù)測.

    3)GM(1,1)-Markov組合模型的預(yù)測精度為Ⅰ級(jí)、Ⅱ級(jí),預(yù)測數(shù)據(jù)可信度高,可推薦運(yùn)用到煤礦井下煤層瓦斯壓力的非線性預(yù)測中.

    猜你喜歡
    采區(qū)殘差灰色
    稠油熱采區(qū)塊冷采降粘技術(shù)應(yīng)用
    基于雙向GRU與殘差擬合的車輛跟馳建模
    基于殘差學(xué)習(xí)的自適應(yīng)無人機(jī)目標(biāo)跟蹤算法
    淺灰色的小豬
    基于遞歸殘差網(wǎng)絡(luò)的圖像超分辨率重建
    灰色時(shí)代
    Coco薇(2017年2期)2017-04-25 17:59:38
    她、它的灰色時(shí)髦觀
    Coco薇(2017年2期)2017-04-25 17:57:49
    基于Delphi-TOPSIS法的露天礦采區(qū)接續(xù)方案優(yōu)選
    中國煤炭(2016年9期)2016-06-15 20:29:54
    感覺
    平穩(wěn)自相關(guān)過程的殘差累積和控制圖
    河南科技(2015年8期)2015-03-11 16:23:52
    寂寞人妻少妇视频99o| 日本免费一区二区三区高清不卡| 中文欧美无线码| 欧美丝袜亚洲另类| 亚洲国产精品成人久久小说| 中文字幕精品亚洲无线码一区| 欧美3d第一页| av黄色大香蕉| 建设人人有责人人尽责人人享有的 | 一边亲一边摸免费视频| 麻豆国产97在线/欧美| 亚洲av.av天堂| 日本午夜av视频| 国产精品一区二区性色av| 日韩av在线免费看完整版不卡| 人妻夜夜爽99麻豆av| 亚洲自拍偷在线| 国产精品国产三级国产av玫瑰| 婷婷色综合大香蕉| 国产精品国产三级国产av玫瑰| 亚洲综合色惰| 婷婷色麻豆天堂久久 | 舔av片在线| 国产精品人妻久久久影院| 国产视频首页在线观看| 成人无遮挡网站| 免费黄色在线免费观看| 成人亚洲精品av一区二区| 啦啦啦观看免费观看视频高清| 成人欧美大片| 99久久精品一区二区三区| 国产黄片视频在线免费观看| 久久久久久久久久久免费av| 国产精品麻豆人妻色哟哟久久 | 国产成人免费观看mmmm| 91精品国产九色| 久久99热6这里只有精品| 精品无人区乱码1区二区| 国产单亲对白刺激| 亚洲人成网站在线播| 一级av片app| 中文字幕亚洲精品专区| 内地一区二区视频在线| 欧美极品一区二区三区四区| 99久久成人亚洲精品观看| 老女人水多毛片| 亚洲人成网站高清观看| 成人欧美大片| 内射极品少妇av片p| 国产精品一区二区三区四区免费观看| 色网站视频免费| 成人高潮视频无遮挡免费网站| 激情 狠狠 欧美| 亚洲欧美清纯卡通| 国产真实伦视频高清在线观看| 国产爱豆传媒在线观看| 熟妇人妻久久中文字幕3abv| 国产人妻一区二区三区在| 男女下面进入的视频免费午夜| 欧美高清成人免费视频www| 国产亚洲精品av在线| 亚洲色图av天堂| 身体一侧抽搐| 国产 一区精品| 最近最新中文字幕免费大全7| 色综合站精品国产| 欧美激情在线99| 国产日韩欧美在线精品| 亚洲电影在线观看av| 国产精品久久视频播放| 国产精品1区2区在线观看.| 久久久国产成人免费| 村上凉子中文字幕在线| 欧美人与善性xxx| 天堂网av新在线| 日本一本二区三区精品| 精品久久久噜噜| 久久精品国产亚洲av涩爱| 久久久久国产网址| 尾随美女入室| 国产激情偷乱视频一区二区| 女的被弄到高潮叫床怎么办| av卡一久久| 我要搜黄色片| 午夜福利视频1000在线观看| 久久99热6这里只有精品| av在线老鸭窝| 特级一级黄色大片| 三级国产精品欧美在线观看| av国产免费在线观看| 变态另类丝袜制服| 成人午夜高清在线视频| 在线免费十八禁| 国产一级毛片在线| 国产麻豆成人av免费视频| 一个人免费在线观看电影| 日韩国内少妇激情av| 婷婷色av中文字幕| 非洲黑人性xxxx精品又粗又长| 国产在线一区二区三区精 | 一级黄色大片毛片| 午夜福利在线在线| 日韩三级伦理在线观看| 久久99热这里只频精品6学生 | 国产激情偷乱视频一区二区| 夫妻性生交免费视频一级片| 国产精品国产三级专区第一集| 男的添女的下面高潮视频| 国产av码专区亚洲av| 日本黄色视频三级网站网址| 国产一区二区三区av在线| 国产精品99久久久久久久久| 一级黄片播放器| 赤兔流量卡办理| av天堂中文字幕网| 日韩精品有码人妻一区| 亚洲最大成人av| 国产精品嫩草影院av在线观看| 亚洲精品456在线播放app| 国产一区有黄有色的免费视频 | 成人av在线播放网站| 国产久久久一区二区三区| 免费看光身美女| 深夜a级毛片| 搡老妇女老女人老熟妇| 少妇熟女欧美另类| 国产精品电影一区二区三区| 在线观看66精品国产| 亚洲精品久久久久久婷婷小说 | 干丝袜人妻中文字幕| 波多野结衣巨乳人妻| 全区人妻精品视频| 村上凉子中文字幕在线| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产成人a区在线观看| 直男gayav资源| 老女人水多毛片| 午夜精品一区二区三区免费看| 国产成人a区在线观看| 黄色欧美视频在线观看| 亚洲欧美日韩东京热| 亚洲最大成人中文| 亚洲欧美日韩高清专用| 中文天堂在线官网| 日本黄大片高清| 亚洲精品自拍成人| 国产精品.久久久| 日本免费在线观看一区| 久久精品国产自在天天线| 亚洲自偷自拍三级| 一二三四中文在线观看免费高清| 狂野欧美白嫩少妇大欣赏| 国产 一区 欧美 日韩| 国产精品乱码一区二三区的特点| 日韩人妻高清精品专区| 色吧在线观看| 一二三四中文在线观看免费高清| 国产乱来视频区| 成人综合一区亚洲| 三级男女做爰猛烈吃奶摸视频| av播播在线观看一区| 变态另类丝袜制服| 亚洲欧美清纯卡通| 内地一区二区视频在线| 麻豆成人av视频| 日韩av在线免费看完整版不卡| 久久久久久久久久黄片| 欧美成人a在线观看| 桃色一区二区三区在线观看| 99热精品在线国产| 国产熟女欧美一区二区| 亚洲欧洲国产日韩| 少妇裸体淫交视频免费看高清| av女优亚洲男人天堂| 国产 一区 欧美 日韩| 久久久久久久久久久免费av| 只有这里有精品99| 丰满乱子伦码专区| 久久久久网色| 午夜久久久久精精品| 色综合色国产| 日韩欧美精品v在线| 久久久久久大精品| 在线免费十八禁| 国产成人a区在线观看| 国产老妇女一区| 精品人妻视频免费看| 国产精品爽爽va在线观看网站| 1024手机看黄色片| 观看免费一级毛片| 色网站视频免费| 亚洲美女视频黄频| av在线播放精品| 亚洲一区高清亚洲精品| 日韩成人伦理影院| 久久热精品热| 日韩三级伦理在线观看| 日韩高清综合在线| 国产精品久久电影中文字幕| 天天一区二区日本电影三级| av播播在线观看一区| 中文字幕熟女人妻在线| 三级毛片av免费| 日韩一区二区视频免费看| 精品一区二区三区人妻视频| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 99热这里只有是精品在线观看| 免费一级毛片在线播放高清视频| 亚洲第一区二区三区不卡| 99久久精品热视频| 色播亚洲综合网| 嘟嘟电影网在线观看| 欧美人与善性xxx| 2021天堂中文幕一二区在线观| 最近视频中文字幕2019在线8| 精品免费久久久久久久清纯| 午夜日本视频在线| 精华霜和精华液先用哪个| 婷婷色麻豆天堂久久 | 能在线免费看毛片的网站| 搡老妇女老女人老熟妇| 色综合站精品国产| 男女视频在线观看网站免费| 一本久久精品| 韩国高清视频一区二区三区| 久久久色成人| 在线免费十八禁| 日本欧美国产在线视频| 一个人免费在线观看电影| 精品久久久久久久久av| 国产真实伦视频高清在线观看| 在线a可以看的网站| 日韩,欧美,国产一区二区三区 | 哪个播放器可以免费观看大片| 亚洲精品自拍成人| 一级黄色大片毛片| 亚洲最大成人中文| 身体一侧抽搐| 亚洲18禁久久av| 国产片特级美女逼逼视频| 中文字幕人妻熟人妻熟丝袜美| 91精品一卡2卡3卡4卡| 久久久久久久久久久免费av| 国产高潮美女av| 欧美一区二区国产精品久久精品| 欧美激情在线99| 亚洲最大成人手机在线| 国产在线一区二区三区精 | 在线观看av片永久免费下载| 深夜a级毛片| 淫秽高清视频在线观看| 国产免费视频播放在线视频 | 午夜老司机福利剧场| 美女内射精品一级片tv| 国产私拍福利视频在线观看| 99久久中文字幕三级久久日本| 国产一区亚洲一区在线观看| av又黄又爽大尺度在线免费看 | 又粗又硬又长又爽又黄的视频| 中文字幕熟女人妻在线| 国产人妻一区二区三区在| 国产日韩欧美在线精品| 午夜福利在线观看免费完整高清在| 国产国拍精品亚洲av在线观看| 三级毛片av免费| 欧美激情国产日韩精品一区| 欧美激情久久久久久爽电影| 欧美潮喷喷水| 日本一本二区三区精品| 精品久久久久久久久久久久久| 美女高潮的动态| 欧美成人精品欧美一级黄| 久久精品国产亚洲网站| 蜜桃亚洲精品一区二区三区| 女人久久www免费人成看片 | 啦啦啦啦在线视频资源| 国产伦一二天堂av在线观看| 黄片无遮挡物在线观看| 国产男人的电影天堂91| 精品免费久久久久久久清纯| 免费在线观看成人毛片| 啦啦啦观看免费观看视频高清| av福利片在线观看| 日产精品乱码卡一卡2卡三| 国产私拍福利视频在线观看| 亚洲乱码一区二区免费版| 日韩一本色道免费dvd| 黄色日韩在线| 国内精品一区二区在线观看| 日韩在线高清观看一区二区三区| 亚洲精品aⅴ在线观看| 淫秽高清视频在线观看| 欧美成人免费av一区二区三区| 国产探花极品一区二区| 少妇丰满av| 男插女下体视频免费在线播放| 99久久精品一区二区三区| 精品一区二区三区人妻视频| 久久久久精品久久久久真实原创| 精华霜和精华液先用哪个| 99热网站在线观看| 床上黄色一级片| 亚洲av免费在线观看| 欧美精品国产亚洲| 亚洲国产欧美人成| 久久久精品欧美日韩精品| 亚洲av电影在线观看一区二区三区 | 精品免费久久久久久久清纯| 秋霞在线观看毛片| 内地一区二区视频在线| 老师上课跳d突然被开到最大视频| 精品99又大又爽又粗少妇毛片| 直男gayav资源| 亚洲国产精品合色在线| 久久婷婷人人爽人人干人人爱| 人体艺术视频欧美日本| 国产午夜精品久久久久久一区二区三区| 精品免费久久久久久久清纯| 天堂影院成人在线观看| 22中文网久久字幕| 免费看美女性在线毛片视频| 高清毛片免费看| 国产免费福利视频在线观看| 91久久精品国产一区二区三区| 亚洲精品456在线播放app| 亚洲成人av在线免费| 国产一区二区在线av高清观看| 麻豆av噜噜一区二区三区| 一个人看的www免费观看视频| 国产男人的电影天堂91| 韩国高清视频一区二区三区| 欧美最新免费一区二区三区| 九色成人免费人妻av| 干丝袜人妻中文字幕| 精品国产一区二区三区久久久樱花 | 亚洲欧美成人综合另类久久久 | 在现免费观看毛片| 国产高潮美女av| 国产高清有码在线观看视频| 中文乱码字字幕精品一区二区三区 | 欧美+日韩+精品| 91午夜精品亚洲一区二区三区| 久久久午夜欧美精品| 午夜a级毛片| 免费观看的影片在线观看| 成人特级av手机在线观看| 啦啦啦啦在线视频资源| av专区在线播放| 欧美丝袜亚洲另类| 午夜久久久久精精品| h日本视频在线播放| 日本wwww免费看| av国产久精品久网站免费入址| 久久婷婷人人爽人人干人人爱| 亚洲中文字幕日韩| 男女下面进入的视频免费午夜| 久久韩国三级中文字幕| 成人亚洲精品av一区二区| 国产日韩欧美在线精品| 永久免费av网站大全| 看非洲黑人一级黄片| 丰满人妻一区二区三区视频av| 噜噜噜噜噜久久久久久91| 亚洲av熟女| 国产人妻一区二区三区在| 青春草国产在线视频| 婷婷色av中文字幕| 老司机福利观看| 夜夜爽夜夜爽视频| 看片在线看免费视频| www.av在线官网国产| 国产私拍福利视频在线观看| 久久人人爽人人片av| 在线天堂最新版资源| 听说在线观看完整版免费高清| videos熟女内射| 日韩一本色道免费dvd| 一级毛片久久久久久久久女| 99热这里只有精品一区| 亚洲中文字幕一区二区三区有码在线看| 国产一区二区亚洲精品在线观看| 三级毛片av免费| 两性午夜刺激爽爽歪歪视频在线观看| 成人国产麻豆网| 美女内射精品一级片tv| 亚洲国产精品国产精品| 精品国产一区二区三区久久久樱花 | 人人妻人人看人人澡| 黄色配什么色好看| 亚洲成人精品中文字幕电影| 18禁在线无遮挡免费观看视频| 国产乱来视频区| 中文乱码字字幕精品一区二区三区 | 亚洲18禁久久av| 日日撸夜夜添| 午夜免费激情av| 老女人水多毛片| 国产乱来视频区| 亚洲精华国产精华液的使用体验| 国产熟女欧美一区二区| 国产在线男女| 免费观看在线日韩| 我的老师免费观看完整版| 日韩精品有码人妻一区| 大香蕉97超碰在线| 欧美潮喷喷水| 日日啪夜夜撸| 91久久精品国产一区二区三区| 国产精品福利在线免费观看| 淫秽高清视频在线观看| 成人av在线播放网站| 成人美女网站在线观看视频| 黄色配什么色好看| 99在线人妻在线中文字幕| 午夜爱爱视频在线播放| 欧美不卡视频在线免费观看| 在线观看av片永久免费下载| 69av精品久久久久久| 欧美丝袜亚洲另类| 国产在线男女| 在线观看av片永久免费下载| 五月玫瑰六月丁香| 国产精华一区二区三区| 国产视频内射| 黄色欧美视频在线观看| 国产午夜福利久久久久久| 免费不卡的大黄色大毛片视频在线观看 | 91狼人影院| 最近最新中文字幕免费大全7| 日韩亚洲欧美综合| 久久人妻av系列| 永久免费av网站大全| 国产亚洲91精品色在线| 午夜免费激情av| 乱码一卡2卡4卡精品| 女人久久www免费人成看片 | 一边亲一边摸免费视频| 性插视频无遮挡在线免费观看| 麻豆乱淫一区二区| 色尼玛亚洲综合影院| 欧美不卡视频在线免费观看| 亚洲电影在线观看av| 久久精品国产99精品国产亚洲性色| 久久这里有精品视频免费| 亚洲国产精品合色在线| 亚洲真实伦在线观看| 日本午夜av视频| 九九久久精品国产亚洲av麻豆| 国产精品美女特级片免费视频播放器| 97超视频在线观看视频| 麻豆乱淫一区二区| 少妇高潮的动态图| 久久久久久久午夜电影| 级片在线观看| 边亲边吃奶的免费视频| АⅤ资源中文在线天堂| 久久精品国产鲁丝片午夜精品| 亚洲av日韩在线播放| 国产精品日韩av在线免费观看| 又黄又爽又刺激的免费视频.| 久久婷婷人人爽人人干人人爱| 国产精品一二三区在线看| 日本黄色视频三级网站网址| 免费搜索国产男女视频| 国产免费男女视频| 亚洲精品国产成人久久av| 老师上课跳d突然被开到最大视频| 久久国产乱子免费精品| 久久精品国产亚洲网站| 天天躁夜夜躁狠狠久久av| 国产精品一区www在线观看| 日本黄色视频三级网站网址| 18禁在线无遮挡免费观看视频| 伦理电影大哥的女人| 精品免费久久久久久久清纯| 中文资源天堂在线| 久久久久久大精品| 久久这里有精品视频免费| 九色成人免费人妻av| 亚洲在线观看片| 黑人高潮一二区| 中文在线观看免费www的网站| 亚洲精品日韩在线中文字幕| 国产真实伦视频高清在线观看| 亚洲欧美一区二区三区国产| 国产免费福利视频在线观看| 国产亚洲精品久久久com| 三级经典国产精品| 成人午夜高清在线视频| 欧美激情久久久久久爽电影| 一区二区三区四区激情视频| 色综合色国产| 美女内射精品一级片tv| 夜夜看夜夜爽夜夜摸| 99在线视频只有这里精品首页| 两个人的视频大全免费| 少妇猛男粗大的猛烈进出视频 | 在线免费观看的www视频| 久久精品国产鲁丝片午夜精品| 最新中文字幕久久久久| 又爽又黄无遮挡网站| 免费观看a级毛片全部| 国产精品蜜桃在线观看| 日韩中字成人| 亚洲欧美精品综合久久99| 村上凉子中文字幕在线| 建设人人有责人人尽责人人享有的 | 观看免费一级毛片| 最近最新中文字幕免费大全7| 午夜久久久久精精品| 99久久中文字幕三级久久日本| 国产色爽女视频免费观看| 在线免费观看的www视频| 十八禁国产超污无遮挡网站| 99久久中文字幕三级久久日本| 久久久a久久爽久久v久久| 色吧在线观看| 久久人妻av系列| 五月玫瑰六月丁香| 国产日韩欧美在线精品| 五月玫瑰六月丁香| 欧美另类亚洲清纯唯美| 国产精品综合久久久久久久免费| 麻豆成人午夜福利视频| 成人毛片a级毛片在线播放| 国产精华一区二区三区| 成年女人看的毛片在线观看| 久久久精品欧美日韩精品| 国产精品久久电影中文字幕| 中文欧美无线码| 午夜激情福利司机影院| 人人妻人人澡人人爽人人夜夜 | 国产一区二区在线观看日韩| 欧美极品一区二区三区四区| 久久国产乱子免费精品| 18禁动态无遮挡网站| 久久久久精品久久久久真实原创| 蜜桃久久精品国产亚洲av| 精品一区二区免费观看| 国产又黄又爽又无遮挡在线| 爱豆传媒免费全集在线观看| 国产成人a区在线观看| 国产成人精品婷婷| 精品久久久久久久人妻蜜臀av| 噜噜噜噜噜久久久久久91| 中国美白少妇内射xxxbb| 校园人妻丝袜中文字幕| 大香蕉久久网| 国产成人a∨麻豆精品| 插逼视频在线观看| 精品欧美国产一区二区三| 日本免费一区二区三区高清不卡| 最近2019中文字幕mv第一页| 男人舔奶头视频| 免费看日本二区| 色综合亚洲欧美另类图片| 观看免费一级毛片| 长腿黑丝高跟| 国产午夜精品一二区理论片| 亚洲最大成人av| 国产精品久久久久久av不卡| 亚洲欧美成人精品一区二区| 久久久精品欧美日韩精品| 日日摸夜夜添夜夜添av毛片| 久久久精品94久久精品| 看黄色毛片网站| 99在线人妻在线中文字幕| 中文在线观看免费www的网站| 免费电影在线观看免费观看| 男女那种视频在线观看| 亚洲熟妇中文字幕五十中出| 国产精品av视频在线免费观看| 秋霞在线观看毛片| 国产免费福利视频在线观看| 欧美xxxx性猛交bbbb| 精品人妻偷拍中文字幕| 秋霞在线观看毛片| 中文欧美无线码| 亚洲av不卡在线观看| 99热网站在线观看| 亚洲精品乱码久久久v下载方式| 久久亚洲精品不卡| av线在线观看网站| 99久久中文字幕三级久久日本| 在线播放国产精品三级| 亚洲国产精品专区欧美| 老师上课跳d突然被开到最大视频| 久久久久久伊人网av| videos熟女内射| 美女国产视频在线观看| 老师上课跳d突然被开到最大视频| 日韩制服骚丝袜av| 老师上课跳d突然被开到最大视频| 高清日韩中文字幕在线| 又爽又黄a免费视频| av又黄又爽大尺度在线免费看 | 日本熟妇午夜| 美女内射精品一级片tv| 寂寞人妻少妇视频99o| 51国产日韩欧美| 成人三级黄色视频| 99久久成人亚洲精品观看| 国产成人福利小说| 搞女人的毛片| 一区二区三区免费毛片| 亚洲精品aⅴ在线观看| 麻豆精品久久久久久蜜桃| 26uuu在线亚洲综合色| 国产高清视频在线观看网站| 日韩欧美精品免费久久| 精品人妻偷拍中文字幕| 秋霞在线观看毛片| 最近中文字幕2019免费版| 在线免费观看不下载黄p国产| 亚洲精品成人久久久久久|