• <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
    免费看日本二区| 男人狂女人下面高潮的视频| 热99在线观看视频| 在线天堂最新版资源| 真实男女啪啪啪动态图| 少妇人妻精品综合一区二区| 爱豆传媒免费全集在线观看| 街头女战士在线观看网站| 内射极品少妇av片p| 99热这里只有精品一区| 亚洲美女搞黄在线观看| 色尼玛亚洲综合影院| 精品99又大又爽又粗少妇毛片| 白带黄色成豆腐渣| 春色校园在线视频观看| 麻豆国产97在线/欧美| 少妇丰满av| 中文天堂在线官网| 欧美激情久久久久久爽电影| 80岁老熟妇乱子伦牲交| 嫩草影院精品99| 亚洲第一区二区三区不卡| 国产日韩欧美在线精品| 人妻系列 视频| 日日摸夜夜添夜夜爱| 有码 亚洲区| 日韩中字成人| 成人亚洲精品一区在线观看 | 狂野欧美白嫩少妇大欣赏| 青青草视频在线视频观看| 2018国产大陆天天弄谢| av播播在线观看一区| 国产精品一区二区三区四区免费观看| 亚洲av中文字字幕乱码综合| 久久久久久久大尺度免费视频| 五月天丁香电影| 日本爱情动作片www.在线观看| 久久亚洲国产成人精品v| 九九爱精品视频在线观看| 熟女人妻精品中文字幕| 深爱激情五月婷婷| 免费电影在线观看免费观看| 国产 一区精品| 永久免费av网站大全| 亚洲经典国产精华液单| 欧美潮喷喷水| 99久久人妻综合| 男人爽女人下面视频在线观看| 欧美一区二区亚洲| 欧美日韩视频高清一区二区三区二| 亚洲av电影在线观看一区二区三区 | 日韩成人av中文字幕在线观看| 男女边吃奶边做爰视频| 久久精品国产自在天天线| 国产一区二区三区av在线| 看非洲黑人一级黄片| 日韩制服骚丝袜av| 精品久久久久久久久av| av在线蜜桃| 九九在线视频观看精品| 国产成人精品福利久久| 大陆偷拍与自拍| 免费观看在线日韩| 国产白丝娇喘喷水9色精品| 欧美人与善性xxx| 国产成人一区二区在线| 午夜免费观看性视频| 成人美女网站在线观看视频| 午夜福利成人在线免费观看| 国产麻豆成人av免费视频| 狂野欧美白嫩少妇大欣赏| 久久久色成人| 亚洲av免费高清在线观看| 欧美日韩国产mv在线观看视频 | 精品国内亚洲2022精品成人| 午夜福利成人在线免费观看| 久久久久久久国产电影| 18禁动态无遮挡网站| 午夜激情福利司机影院| 五月天丁香电影| 国产精品一区www在线观看| 美女高潮的动态| 女人被狂操c到高潮| 极品教师在线视频| 国产成人aa在线观看| 成年人午夜在线观看视频 | 国产 一区 欧美 日韩| 99久久人妻综合| 噜噜噜噜噜久久久久久91| 两个人的视频大全免费| 精品一区二区三区人妻视频| 全区人妻精品视频| 日韩在线高清观看一区二区三区| 久久人人爽人人片av| 成人二区视频| 极品教师在线视频| 国产69精品久久久久777片| 久久久久久久久久久丰满| 国产久久久一区二区三区| 亚洲av不卡在线观看| 国产精品一二三区在线看| 欧美极品一区二区三区四区| 一级av片app| 久热久热在线精品观看| 国产乱人偷精品视频| 十八禁国产超污无遮挡网站| 免费看不卡的av| 亚洲国产日韩欧美精品在线观看| 免费看光身美女| 人人妻人人澡人人爽人人夜夜 | 国产大屁股一区二区在线视频| 午夜久久久久精精品| 99热这里只有是精品在线观看| 嫩草影院入口| 日本免费a在线| 国产精品精品国产色婷婷| 成人漫画全彩无遮挡| 国产精品女同一区二区软件| 久久久久久久午夜电影| 激情五月婷婷亚洲| 国产精品女同一区二区软件| 欧美一区二区亚洲| 久久久久久久久大av| 亚洲真实伦在线观看| 国产爱豆传媒在线观看| 国产成人精品久久久久久| 国产精品一区二区在线观看99 | 欧美+日韩+精品| 尤物成人国产欧美一区二区三区| 国产精品伦人一区二区| 99热这里只有是精品在线观看| 国产视频内射| 国产麻豆成人av免费视频| 国产黄色视频一区二区在线观看| 国产人妻一区二区三区在| 高清欧美精品videossex| 午夜老司机福利剧场| 欧美zozozo另类| 久久综合国产亚洲精品| 少妇裸体淫交视频免费看高清| 国产精品爽爽va在线观看网站| 成人漫画全彩无遮挡| 一级毛片久久久久久久久女| 亚洲av电影在线观看一区二区三区 | 国产日韩欧美在线精品| 欧美成人精品欧美一级黄| 中文字幕免费在线视频6| videos熟女内射| 亚洲精品一二三| 大陆偷拍与自拍| 久久精品人妻少妇| 精品久久久久久久久亚洲| 精品久久久久久电影网| 麻豆国产97在线/欧美| 亚洲国产欧美在线一区| 内地一区二区视频在线| 国产一区有黄有色的免费视频 | 亚洲精品一二三| 久久久久久久久久成人| 成年版毛片免费区| av天堂中文字幕网| 欧美三级亚洲精品| 中文字幕av成人在线电影| 亚洲精品视频女| 亚洲欧美成人综合另类久久久| 色吧在线观看| 激情五月婷婷亚洲| 少妇猛男粗大的猛烈进出视频 | 哪个播放器可以免费观看大片| 男人舔女人下体高潮全视频| 男人舔女人下体高潮全视频| 日韩一本色道免费dvd| av网站免费在线观看视频 | 久久久久久伊人网av| 成人国产麻豆网| 高清日韩中文字幕在线| 国产精品国产三级国产av玫瑰| 男人舔奶头视频| 国产精品一区二区三区四区免费观看| 精品酒店卫生间| 亚洲电影在线观看av| 在线观看av片永久免费下载| 亚洲精品aⅴ在线观看| 亚洲av中文av极速乱| 免费高清在线观看视频在线观看| 久久99热这里只有精品18| 久久精品国产自在天天线| 久久久成人免费电影| 国产精品三级大全| 一二三四中文在线观看免费高清| 好男人在线观看高清免费视频| 亚洲精品国产av蜜桃| 成人亚洲精品av一区二区| 天堂影院成人在线观看| 九九爱精品视频在线观看| 欧美97在线视频| 人体艺术视频欧美日本| 亚洲国产精品sss在线观看| 2018国产大陆天天弄谢| 亚洲高清免费不卡视频| 亚洲四区av| 久久人人爽人人片av| 在线观看人妻少妇| 欧美激情国产日韩精品一区| 边亲边吃奶的免费视频| 国产三级在线视频| 久久久久久久久大av| 视频中文字幕在线观看| 91精品伊人久久大香线蕉| 成人漫画全彩无遮挡| 欧美性猛交╳xxx乱大交人| 国产精品久久久久久精品电影| 三级男女做爰猛烈吃奶摸视频| 亚洲av成人精品一二三区| 91久久精品国产一区二区三区| 18禁在线播放成人免费| 国产av国产精品国产| 国产人妻一区二区三区在| 99久国产av精品| 国产伦在线观看视频一区| 国产 一区精品| 男女国产视频网站| 日本一二三区视频观看| 乱系列少妇在线播放| 亚洲精品乱码久久久v下载方式| 成人二区视频| 岛国毛片在线播放| 亚洲精品亚洲一区二区| 少妇丰满av| 国产成人91sexporn| 国产在视频线精品| 晚上一个人看的免费电影| 亚洲精华国产精华液的使用体验| 美女大奶头视频| 久久精品久久精品一区二区三区| 51国产日韩欧美| 亚洲av在线观看美女高潮| 小蜜桃在线观看免费完整版高清| 亚洲在久久综合| 岛国毛片在线播放| 男女边摸边吃奶| 国产成人精品婷婷| 日日干狠狠操夜夜爽| 丝袜美腿在线中文| 99热这里只有是精品在线观看| 欧美+日韩+精品| 国产精品国产三级专区第一集| 麻豆国产97在线/欧美| av卡一久久| 最近2019中文字幕mv第一页| 日韩电影二区| 男女啪啪激烈高潮av片| 国产成人精品婷婷| 久久午夜福利片| 亚洲四区av| 国产精品久久久久久久久免| 国产成人一区二区在线| 亚洲欧美成人综合另类久久久| 69人妻影院| 韩国高清视频一区二区三区| 亚洲av在线观看美女高潮| 国产精品人妻久久久久久| 欧美不卡视频在线免费观看| 人人妻人人澡欧美一区二区| 我要看日韩黄色一级片| 非洲黑人性xxxx精品又粗又长| 午夜视频国产福利| 日本与韩国留学比较| 免费大片18禁| 女人被狂操c到高潮| 精品少妇黑人巨大在线播放| 午夜福利在线观看免费完整高清在| 日韩av在线免费看完整版不卡| 少妇人妻一区二区三区视频| 婷婷六月久久综合丁香| 夫妻午夜视频| 国产精品嫩草影院av在线观看| 禁无遮挡网站| 国产精品蜜桃在线观看| 国产在线男女| 精品国产露脸久久av麻豆 | 亚洲精品视频女| 国产精品爽爽va在线观看网站| 国产精品伦人一区二区| 国产欧美日韩精品一区二区| 国产免费一级a男人的天堂| 神马国产精品三级电影在线观看| ponron亚洲| 国产一区亚洲一区在线观看| 青青草视频在线视频观看| 麻豆乱淫一区二区| 亚洲性久久影院| 少妇人妻一区二区三区视频| 嫩草影院入口| 亚洲三级黄色毛片| 亚洲精品视频女| 国产av国产精品国产| 天堂影院成人在线观看| 中文字幕av成人在线电影| 日本免费a在线| 久久精品国产亚洲av天美| 亚洲电影在线观看av| 91久久精品国产一区二区三区| 国产v大片淫在线免费观看| 成人高潮视频无遮挡免费网站| 99久久中文字幕三级久久日本| 狂野欧美白嫩少妇大欣赏| 国产一区二区三区综合在线观看 | 在线观看av片永久免费下载| 精品一区二区免费观看| av卡一久久| 极品教师在线视频| 亚洲国产日韩欧美精品在线观看| 男女啪啪激烈高潮av片| 色吧在线观看| 国产亚洲av嫩草精品影院| 18禁动态无遮挡网站| 麻豆成人av视频| 国产极品天堂在线| 日韩三级伦理在线观看| 国产在视频线精品| 中文字幕免费在线视频6| 亚洲,欧美,日韩| 久久热精品热| 国产黄a三级三级三级人| 大片免费播放器 马上看| 国产精品伦人一区二区| 国产在线男女| 亚洲精品aⅴ在线观看| 97超视频在线观看视频| 六月丁香七月| 蜜桃久久精品国产亚洲av| 日韩国内少妇激情av| 国产三级在线视频| or卡值多少钱| 亚洲精品成人av观看孕妇| 午夜视频国产福利| 在线观看av片永久免费下载| 中文精品一卡2卡3卡4更新| 能在线免费看毛片的网站| 熟妇人妻不卡中文字幕| 亚洲欧美日韩东京热| 黄片无遮挡物在线观看| 亚洲av电影在线观看一区二区三区 | 天堂中文最新版在线下载 | av专区在线播放| 天堂中文最新版在线下载 | 女人久久www免费人成看片| 亚洲乱码一区二区免费版| 免费看美女性在线毛片视频| 如何舔出高潮| 亚洲在久久综合| 少妇人妻一区二区三区视频| 能在线免费看毛片的网站| 又黄又爽又刺激的免费视频.| 赤兔流量卡办理| 高清午夜精品一区二区三区| 日韩三级伦理在线观看| 亚洲av免费高清在线观看| 99热这里只有是精品50| 免费观看在线日韩| 丰满人妻一区二区三区视频av| av国产久精品久网站免费入址| 一级爰片在线观看| 亚洲欧美一区二区三区黑人 | 国产探花在线观看一区二区| 国产午夜福利久久久久久| 亚洲国产欧美在线一区| 亚洲最大成人手机在线| 久久精品熟女亚洲av麻豆精品 | 少妇的逼好多水| 人妻少妇偷人精品九色| 午夜免费观看性视频| 在线观看av片永久免费下载| 日本猛色少妇xxxxx猛交久久| 精品久久久精品久久久| 极品少妇高潮喷水抽搐| 99视频精品全部免费 在线| 国产成人免费观看mmmm| 亚洲精品国产av成人精品| 尤物成人国产欧美一区二区三区| 亚洲av二区三区四区| 久久久久精品久久久久真实原创| 亚洲aⅴ乱码一区二区在线播放| 欧美成人午夜免费资源| 国产伦在线观看视频一区| 2021天堂中文幕一二区在线观| 国产精品精品国产色婷婷| 欧美丝袜亚洲另类| kizo精华| 看非洲黑人一级黄片| 亚洲婷婷狠狠爱综合网| 国产 一区 欧美 日韩| 国产成人精品久久久久久| 麻豆国产97在线/欧美| 麻豆乱淫一区二区| 人人妻人人澡人人爽人人夜夜 | 欧美一区二区亚洲| 99久国产av精品| 免费av不卡在线播放| 国产精品一区二区三区四区久久| 国产淫片久久久久久久久| 午夜精品一区二区三区免费看| 波野结衣二区三区在线| 国产伦精品一区二区三区视频9| 少妇猛男粗大的猛烈进出视频 | 精品亚洲乱码少妇综合久久| 夜夜看夜夜爽夜夜摸| 午夜激情久久久久久久| 午夜爱爱视频在线播放| 亚洲人成网站在线观看播放| 国产 一区精品| videossex国产| 免费观看精品视频网站| 久久久亚洲精品成人影院| 偷拍熟女少妇极品色| 精品99又大又爽又粗少妇毛片| 欧美xxⅹ黑人| 淫秽高清视频在线观看| 中文字幕制服av| 亚洲国产高清在线一区二区三| 综合色av麻豆| 熟妇人妻不卡中文字幕| 国产女主播在线喷水免费视频网站 | 69人妻影院| 中文在线观看免费www的网站| 91久久精品电影网| 久久精品久久精品一区二区三区| 丝袜喷水一区| 91午夜精品亚洲一区二区三区| 日韩av在线大香蕉| 国产 一区 欧美 日韩| 亚洲精品久久久久久婷婷小说| 欧美成人a在线观看| 免费观看精品视频网站| 亚洲欧美日韩无卡精品| 一区二区三区乱码不卡18| 一级毛片aaaaaa免费看小| 国产三级在线视频| 亚洲欧美清纯卡通| 看黄色毛片网站| 欧美日韩在线观看h| av国产久精品久网站免费入址| 一级毛片aaaaaa免费看小| 深夜a级毛片| 18禁裸乳无遮挡免费网站照片| 少妇的逼水好多| 日日撸夜夜添| 久久久久性生活片| 免费无遮挡裸体视频| 2018国产大陆天天弄谢| 麻豆av噜噜一区二区三区| 国产av码专区亚洲av| 欧美性感艳星| 国产在线一区二区三区精| 免费观看无遮挡的男女| 亚洲在线自拍视频| 亚洲精品,欧美精品| 午夜福利视频1000在线观看| 免费av观看视频| 99久久精品热视频| 日韩 亚洲 欧美在线| 国产白丝娇喘喷水9色精品| 欧美另类一区| 高清在线视频一区二区三区| 男人舔奶头视频| 99久国产av精品国产电影| 最近中文字幕2019免费版| 女人十人毛片免费观看3o分钟| 亚洲人成网站在线观看播放| 亚洲美女视频黄频| 亚洲精品第二区| 欧美最新免费一区二区三区| 精品一区在线观看国产| 中文字幕制服av| 男女啪啪激烈高潮av片| 成人无遮挡网站| 中文资源天堂在线| 亚洲av成人av| 男女那种视频在线观看| 亚洲精品日韩av片在线观看| 成年人午夜在线观看视频 | 久久久亚洲精品成人影院| 欧美成人一区二区免费高清观看| 国产精品人妻久久久久久| 有码 亚洲区| 久久精品国产亚洲av天美| 日日啪夜夜撸| 美女cb高潮喷水在线观看| 99re6热这里在线精品视频| 日韩av免费高清视频| 国产在视频线在精品| 91aial.com中文字幕在线观看| 男人爽女人下面视频在线观看| 欧美激情久久久久久爽电影| 97在线视频观看| 国产大屁股一区二区在线视频| 3wmmmm亚洲av在线观看| 国产精品美女特级片免费视频播放器| 天堂中文最新版在线下载 | freevideosex欧美| 春色校园在线视频观看| 免费大片黄手机在线观看| 97在线视频观看| 国产高清三级在线| 国产成人免费观看mmmm| 男女那种视频在线观看| 婷婷六月久久综合丁香| 麻豆av噜噜一区二区三区| 人妻夜夜爽99麻豆av| 嫩草影院入口| 人妻一区二区av| 精品久久国产蜜桃| 极品少妇高潮喷水抽搐| 最新中文字幕久久久久| 免费黄色在线免费观看| 精品不卡国产一区二区三区| 日本免费a在线| a级一级毛片免费在线观看| 精品国产三级普通话版| 女的被弄到高潮叫床怎么办| 边亲边吃奶的免费视频| 最后的刺客免费高清国语| 伊人久久精品亚洲午夜| 男女边吃奶边做爰视频| 十八禁国产超污无遮挡网站| 免费播放大片免费观看视频在线观看| 18禁动态无遮挡网站| 亚洲精品亚洲一区二区| 国产人妻一区二区三区在| 国产熟女欧美一区二区| 久久97久久精品| 久久久久免费精品人妻一区二区| 又粗又硬又长又爽又黄的视频| 日本一本二区三区精品| 国产精品福利在线免费观看| 国产片特级美女逼逼视频| 免费av不卡在线播放| 搡老妇女老女人老熟妇| 免费人成在线观看视频色| 日韩欧美三级三区| 精品久久久久久久久久久久久| 亚洲成人中文字幕在线播放| 色视频www国产| 国产爱豆传媒在线观看| 精品久久久久久久末码| 亚洲欧美日韩东京热| av在线天堂中文字幕| 中国美白少妇内射xxxbb| 成人特级av手机在线观看| 国产免费又黄又爽又色| 精品一区二区免费观看| 亚洲成人一二三区av| 国产亚洲5aaaaa淫片| 久久久久性生活片| 男女啪啪激烈高潮av片| 国产淫片久久久久久久久| 天堂中文最新版在线下载 | 日日啪夜夜撸| 国产精品一区www在线观看| 国产免费一级a男人的天堂| av女优亚洲男人天堂| 午夜福利在线观看免费完整高清在| 18禁在线无遮挡免费观看视频| 黄色配什么色好看| 亚洲av一区综合| 国产探花极品一区二区| 日本免费a在线| 91精品一卡2卡3卡4卡| 最近视频中文字幕2019在线8| 免费av观看视频| 一级毛片电影观看| 联通29元200g的流量卡| 一级片'在线观看视频| av在线老鸭窝| 国产成人免费观看mmmm| 乱码一卡2卡4卡精品| 一级毛片 在线播放| 中文在线观看免费www的网站| 国产一区亚洲一区在线观看| 老女人水多毛片| 在线免费观看的www视频| 日本一本二区三区精品| 黑人高潮一二区| 真实男女啪啪啪动态图| 亚洲av免费在线观看| 国产日韩欧美在线精品| 天堂√8在线中文| 日韩一区二区视频免费看| 国产高清国产精品国产三级 | 国国产精品蜜臀av免费| 日韩欧美一区视频在线观看 | 亚洲欧洲日产国产| 小蜜桃在线观看免费完整版高清| 成人特级av手机在线观看| 欧美高清成人免费视频www| 淫秽高清视频在线观看| 91精品一卡2卡3卡4卡| 国产黄片视频在线免费观看| 天天一区二区日本电影三级| 久久精品久久精品一区二区三区| 黄色配什么色好看| 午夜福利视频精品| 在线a可以看的网站| 日韩欧美精品免费久久| 两个人的视频大全免费| 51国产日韩欧美| 欧美日韩在线观看h| 女人被狂操c到高潮| 看非洲黑人一级黄片| 免费黄网站久久成人精品| 精品久久久久久电影网| 久久久a久久爽久久v久久| 天堂√8在线中文| 久久久精品免费免费高清|