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

    多采空區(qū)下堅(jiān)硬厚層破斷頂板群結(jié)構(gòu)的失穩(wěn)規(guī)律

    2014-06-07 05:55:00劉長(zhǎng)友楊敬軒楊培舉
    煤炭學(xué)報(bào) 2014年3期
    關(guān)鍵詞:厚層礦壓塊體

    劉長(zhǎng)友,楊敬軒,于 斌,楊培舉

    (1.中國(guó)礦業(yè)大學(xué)礦業(yè)工程學(xué)院,江蘇徐州 221116;2.中國(guó)礦業(yè)大學(xué)煤炭資源與安全開(kāi)采國(guó)家重點(diǎn)實(shí)驗(yàn)室,江蘇徐州 221116;3.大同煤礦集團(tuán)有限責(zé)任公司,山西大同 037003)

    多采空區(qū)下堅(jiān)硬厚層破斷頂板群結(jié)構(gòu)的失穩(wěn)規(guī)律

    劉長(zhǎng)友1,2,楊敬軒1,2,于 斌1,3,楊培舉1,2

    (1.中國(guó)礦業(yè)大學(xué)礦業(yè)工程學(xué)院,江蘇徐州 221116;2.中國(guó)礦業(yè)大學(xué)煤炭資源與安全開(kāi)采國(guó)家重點(diǎn)實(shí)驗(yàn)室,江蘇徐州 221116;3.大同煤礦集團(tuán)有限責(zé)任公司,山西大同 037003)

    為弄清多采空區(qū)下堅(jiān)硬頂板群結(jié)構(gòu)的破斷失穩(wěn)對(duì)工作面礦壓的影響,根據(jù)大同礦區(qū)多采空區(qū)堅(jiān)硬厚層破斷頂板群的賦存條件,采用理論分析、相似模擬實(shí)驗(yàn)和現(xiàn)場(chǎng)實(shí)測(cè)分析相結(jié)合的研究方法,對(duì)多采空區(qū)破斷頂板群結(jié)構(gòu)的失穩(wěn)規(guī)律及其對(duì)工作面來(lái)壓的影響進(jìn)行了研究探討。研究表明,多采空區(qū)下堅(jiān)硬厚層破斷頂板群結(jié)構(gòu)的失穩(wěn)具有一定的概率特征;采用威布爾多參數(shù)分布函數(shù)對(duì)破斷頂板群結(jié)構(gòu)的失穩(wěn)形式具有很好的描述,并建立了近距離煤層群多采空區(qū)頂板群結(jié)構(gòu)的失穩(wěn)模型,確定了大同礦區(qū)堅(jiān)硬頂板群結(jié)構(gòu)的失穩(wěn)參數(shù),得到了侏羅系煤層群頂板結(jié)構(gòu)的失穩(wěn)率;通過(guò)對(duì)大同礦區(qū)永定莊煤礦15號(hào)煤層端頭及中部位置礦壓和支架阻力的實(shí)測(cè)分析表明,工作面上方破斷頂板群結(jié)構(gòu)的失穩(wěn)率與工作面支架阻力大小具有相關(guān)性,驗(yàn)證了工作面堅(jiān)硬厚層破斷頂板的失穩(wěn)規(guī)律。

    多采空區(qū);破斷頂板群;結(jié)構(gòu)失穩(wěn);礦壓特征;威布爾分布

    大同礦區(qū)侏羅系煤層埋藏深度較淺、煤層數(shù)量多、層間距小,煤層傾角3°~7°,可采煤層厚度0.8~9.0 m,在下部煤層開(kāi)采時(shí)采空區(qū)頂板冒落易與上部采空區(qū)溝通。因此,對(duì)于大同礦區(qū)堅(jiān)硬厚層頂板近距離煤層群開(kāi)采,當(dāng)上部近距離多煤層已開(kāi)采完畢時(shí),下部煤層在上覆多重煤柱壓力作用以及上部煤層開(kāi)采的多次采動(dòng)影響下,下部煤層工作面開(kāi)采時(shí),覆巖應(yīng)力較為復(fù)雜,頂板也較為破碎,這給下部煤層工作面的回采與支護(hù)帶來(lái)了新的技術(shù)難題[1]。例如大同礦區(qū)侏羅紀(jì)下組煤主要可采煤層有11,12,14,15號(hào)煤層,煤層間距6~15 m,分岔與合并現(xiàn)象十分頻繁,煤層頂板巖性較為堅(jiān)硬且基本由粉砂巖、細(xì)砂巖、中粗粒砂巖互層組成。上部煤層開(kāi)采后形成的多重堅(jiān)硬厚層破斷頂板群失穩(wěn)垮落后經(jīng)過(guò)一定時(shí)間的運(yùn)動(dòng)調(diào)整,塊體間有可能再次形成結(jié)構(gòu)而具有一定的承載能力。當(dāng)下部煤層開(kāi)采時(shí),已經(jīng)破斷的堅(jiān)硬頂板將再次失穩(wěn)形成承載結(jié)構(gòu),并由下而上形成多層堅(jiān)硬破斷頂板群結(jié)構(gòu),而頂板群結(jié)構(gòu)的破斷失穩(wěn)又會(huì)對(duì)下部煤層工作面的礦壓顯現(xiàn)產(chǎn)生影響。目前對(duì)于堅(jiān)硬頂板破斷機(jī)理與控制問(wèn)題的研究[2-9]多側(cè)重于單層頂板結(jié)構(gòu)的破斷與控制以及現(xiàn)場(chǎng)實(shí)測(cè)方面的分析,而對(duì)于多采空區(qū)條件下多層堅(jiān)硬頂板破斷與控制問(wèn)題的探討[10-13]則相對(duì)較少,尤其對(duì)于多采空區(qū)下多層堅(jiān)硬厚層破斷頂板群結(jié)構(gòu)的失穩(wěn)及其與工作面來(lái)壓特征間關(guān)系的探討則更為鮮見(jiàn),因此,針對(duì)大同礦區(qū)侏羅系煤層群的賦存及開(kāi)采條件,分析研究近距離煤層群多采空區(qū)條件下,下部煤層開(kāi)采覆巖離散結(jié)構(gòu)承載失穩(wěn)特征及其與礦壓顯現(xiàn)的相關(guān)性影響規(guī)律,可為大同礦區(qū)多采空區(qū)下開(kāi)采的圍巖控制提供依據(jù)。

    1 破斷頂板群結(jié)構(gòu)的失穩(wěn)特征

    堅(jiān)硬厚層頂板條件下,隨著工作面的推進(jìn),煤層頂板周期性有規(guī)則斷裂失穩(wěn),基于頂板巖性堅(jiān)硬以及厚度較大的原因,破斷塊體形狀保持一定的幾何相似性。在整層煤層開(kāi)采完畢后,厚層頂板破斷塊體在上覆巖層靜載作用下經(jīng)過(guò)一定時(shí)間的運(yùn)動(dòng)調(diào)整,塊體間有可能再次相互咬合,形成具有一定承載能力的砌體結(jié)構(gòu)。隨著下部工作面的推進(jìn),采空區(qū)空間范圍的增大,頂板破斷塊體砌體結(jié)構(gòu)易發(fā)生單層失穩(wěn)或多層同步失穩(wěn)。以三層采空區(qū)堅(jiān)硬厚層頂板條件下采煤為例,頂板破斷塊體砌體結(jié)構(gòu)自下而上依次為Ⅰ,Ⅱ與Ⅲ號(hào),下部煤層開(kāi)采過(guò)程中,多層堅(jiān)硬厚層頂板砌體結(jié)構(gòu)運(yùn)動(dòng)特征如圖1所示。

    圖1 堅(jiān)硬厚層頂板群結(jié)構(gòu)Fig.1 Group structure of hard thick roof

    堅(jiān)硬厚層頂板破斷塊體形成的砌體結(jié)構(gòu)具有一定的承載能力,多采空區(qū)條件下工作面的來(lái)壓特征,主要取決于工作面上方堅(jiān)硬厚層頂板群破斷塊體間的組合承載作用[14-16]。以三層采空區(qū)及三層堅(jiān)硬厚層頂板條件下的煤層開(kāi)采為例,破斷塊體組合承載結(jié)構(gòu)形式如圖2所示。

    圖2 厚層頂板群破斷塊體組合結(jié)構(gòu)Fig.2 Fault block combination structure of hard thick roof

    多采空區(qū)堅(jiān)硬厚層頂板群結(jié)構(gòu)的失穩(wěn)垮斷隨下部煤層的開(kāi)采具有一定的概率特征,即頂板結(jié)構(gòu)在其下位頂板破斷塊體失穩(wěn)后發(fā)生單層失穩(wěn)或多層同步失穩(wěn)的現(xiàn)象具有一定的隨機(jī)性。隨著距離下部工作面高度的不同,不同位置頂板的垮斷失穩(wěn)會(huì)給工作面礦壓帶來(lái)不同影響,此時(shí)工作面來(lái)壓將對(duì)應(yīng)一定的特征強(qiáng)度。因此下部工作面來(lái)壓特征通常與其上部頂板結(jié)構(gòu)距離煤層高度、頂板破斷塊體厚度以及頂板發(fā)生單層失穩(wěn)還是多層同步失穩(wěn)有關(guān)。

    為弄清多采空區(qū)堅(jiān)硬厚層頂板群結(jié)構(gòu)下工作面來(lái)壓特征與頂板破斷塊體結(jié)構(gòu)失穩(wěn)率間的關(guān)系[17],設(shè)每分層頂板破斷塊體結(jié)構(gòu)失穩(wěn)事件為Ai,則工作面開(kāi)采影響范圍內(nèi),第i層頂板破斷塊體結(jié)構(gòu)的失穩(wěn)率為

    式中,i為頂板層位數(shù),自下往上依次為1,2,……;j為頂板破斷塊體序號(hào),自右向左依次為1,2,…,Ni,其中Ni為工作面開(kāi)采影響范圍內(nèi)第i巖層破斷塊體總數(shù);σij為第i層頂板第j塊體承載量;Fij(σij)為第i層頂板第j塊體失穩(wěn)率。

    由式(1)得到破斷頂板群結(jié)構(gòu)多組塊體結(jié)構(gòu)同時(shí)失穩(wěn)的概率為

    式中,n為堅(jiān)硬頂板失穩(wěn)總層數(shù)。

    頂板賦存及巖性基本相同條件下,近水平煤層頂板承載基本一致,工作面推進(jìn)過(guò)程中,頂板來(lái)壓步距趨于常量,頂板破斷塊體間具有一定的幾何相似性。若將每層頂板破斷塊體單元的失穩(wěn)率視為不變,即滿足關(guān)系

    式中,σi為工作面開(kāi)采影響范圍內(nèi)第i層頂板破斷塊體承載量;Fi(σi)為單層頂板破斷塊體失穩(wěn)率,抽樣計(jì)算中其樣本值為同一頂板巖層失穩(wěn)塊體數(shù)與工作面采動(dòng)影響范圍內(nèi)塊體總數(shù)比值。

    聯(lián)立式(1)與式(3),計(jì)算得到工作面開(kāi)采影響范圍內(nèi)第i層頂板塊體結(jié)構(gòu)的失穩(wěn)率為

    同理,計(jì)算得到工作面頂板群結(jié)構(gòu)多組塊體的同步失穩(wěn)率為

    隨著工作面的推進(jìn)與采空區(qū)空間范圍的增加,靠近工作面的頂板巖層首先趨于垮斷失穩(wěn),而上部頂板破斷塊體結(jié)構(gòu)的失穩(wěn)幾率一般小于其下部頂板塊體失穩(wěn)率。上部頂板破斷塊體在其下部頂板結(jié)構(gòu)失穩(wěn)后出現(xiàn)失穩(wěn)的概率由式(5)計(jì)算得到。

    鑒于多采空區(qū)堅(jiān)硬頂板下的工作面采場(chǎng)來(lái)壓影響因素較多,條件較為復(fù)雜,而威布爾多參數(shù)函數(shù)分布又是適用范圍較廣的一種概率分布形式[18-19]。當(dāng)參數(shù)取值不同時(shí),威布爾函數(shù)可以蛻化為指數(shù)分布、瑞利分布和正態(tài)分布形式,且大量實(shí)踐證明,凡是因?yàn)槟骋痪植渴Щ蚬收纤鸬南到y(tǒng)失效問(wèn)題均服從威布爾分布。因此,設(shè)工作面開(kāi)采影響范圍內(nèi)第i層頂板破斷塊體的失穩(wěn)密度f(wàn)(σi)服從威布爾函數(shù)分布,即

    式中,mi為第i巖層形狀參數(shù);ηi為第i巖層尺度參數(shù);σi0為第i巖層位置參數(shù)。

    由此得到工作面開(kāi)采影響范圍內(nèi)第i層頂板單個(gè)塊體的失穩(wěn)率為

    聯(lián)立式(4)與式(9),得到工作面開(kāi)采影響范圍內(nèi)第i層頂板破斷塊體結(jié)構(gòu)的失穩(wěn)率為

    為確定第i層頂板破斷塊體結(jié)構(gòu)失穩(wěn)時(shí)的分布位置參數(shù),對(duì)式(9)進(jìn)行簡(jiǎn)化,兩邊取對(duì)數(shù)得

    式中,xi為第i層頂板承載狀態(tài)換算參數(shù);yi為第i層頂板巖層失穩(wěn)率換算參數(shù)。

    式(11)簡(jiǎn)化為

    定義第i分層頂板的換算參數(shù)xi與yi在第j塊體位置處失穩(wěn)的樣本參數(shù)為(xij,yij),則換算樣本數(shù)(xij,yij)的相關(guān)系數(shù)γi為

    式中,sixixi,siyiyi與sixiyi分別為第i分層關(guān)于換算樣本數(shù)的方差及協(xié)方差。

    為了使堅(jiān)硬厚層頂板破斷塊體在相應(yīng)載荷作用下的失穩(wěn)換算樣本數(shù)具有良好的線性相關(guān)性,對(duì)相應(yīng)載荷下的失穩(wěn)換算樣本(xij,yij)相關(guān)系數(shù)取極值得

    根據(jù)極值條件式(14),可確定第i巖層的位置參數(shù)σi0。

    埋深相差不大的近距離煤層群開(kāi)采,在頂板巖性基本相同條件下,堅(jiān)硬厚層頂板的破斷形式基本相同,若將頂板破斷塊體失穩(wěn)參數(shù)采用相同參數(shù),即滿足如下關(guān)系:式中,m,η,σ0分別為巖性相近的頂板巖層的平均形狀參數(shù)、平均尺度參數(shù)與平均位置參數(shù)。

    將式(9)代入式(6),聯(lián)立式(15),計(jì)算得到近距離煤層頂板群多組破斷塊體結(jié)構(gòu)的同步失穩(wěn)率為

    式中,σ為近距離煤層群頂板破斷塊體平均承載。

    聯(lián)立式(7)與式(16),計(jì)算得到近距離煤層群相似頂板條件下破斷塊體間的失穩(wěn)率為

    由式(10)與式(17)可以看出,近距離煤層群相似頂板結(jié)構(gòu)的條件失穩(wěn)率等于單層頂板破斷塊體的失穩(wěn)率,即

    由式(17)與式(18)可知,在賦存條件基本相同,巖性相近的近距離頂板群結(jié)構(gòu)中,上部巖層破斷塊體在其下部頂板破斷塊體失穩(wěn)后出現(xiàn)失穩(wěn)現(xiàn)象的概率取決于工作面影響范圍內(nèi)本層破斷塊體數(shù)目。

    近距離煤層群開(kāi)采,多采空區(qū)頂板巖性相近條件下,同分層頂板破斷塊體失穩(wěn)造成的工作面礦壓顯現(xiàn)特征基本相似。為探討多采空區(qū)巖性相似頂板下的采場(chǎng)來(lái)壓特征,定義頂板結(jié)構(gòu)失穩(wěn)的加權(quán)平均失穩(wěn)率,其計(jì)算式為

    2 頂板群結(jié)構(gòu)失穩(wěn)的實(shí)驗(yàn)研究

    2.1 煤層賦存條件

    大同礦區(qū)侏羅統(tǒng)大同組煤系煤田面積約772 km2,含煤地層總厚75~264 m,平均210 m,可采煤層15層。煤層結(jié)構(gòu)相對(duì)較為致密,裂隙不發(fā)育,單向抗壓強(qiáng)度較高,普氏系數(shù)f>3,除2號(hào)煤層頂板為云崗組礫巖、砂礫巖組成外,其他巖層均為粉砂巖、細(xì)砂巖、中粗粒砂巖互層組成,節(jié)理裂隙不發(fā)育,連續(xù)性強(qiáng),結(jié)構(gòu)完整,抗壓強(qiáng)度在55.20~65.63 MPa。

    大同礦區(qū)上組煤層基本采完,主要開(kāi)采下組11, 12,14,15號(hào)煤層,煤層層間距離6~15 m,分叉合并頻繁,在開(kāi)采過(guò)程中下部煤層頂板受上部煤層開(kāi)采的影響很大。大同礦區(qū)11~15號(hào)煤層柱狀圖如圖3所示。

    2.2 厚層頂板群結(jié)構(gòu)的相似模擬實(shí)驗(yàn)

    為了分析侏羅系煤層群多采空區(qū)條件下,堅(jiān)硬厚層頂板群結(jié)構(gòu)的破斷失穩(wěn)規(guī)律及其時(shí)空影響關(guān)系,采用二維平面應(yīng)力實(shí)驗(yàn)臺(tái)進(jìn)行模擬分析。實(shí)驗(yàn)臺(tái)的尺寸為:長(zhǎng)×寬×高=5 000 mm×300 mm×3 000 mm,實(shí)驗(yàn)?zāi)P蛯?shí)際鋪設(shè)高度為2 000 mm。采用砂子、碳酸鈣與石膏和砂子、石灰與石膏相似材料,模擬不同巖性和強(qiáng)度的頂板巖層。模型幾何相似比例為1∶150;容重相似之比為1∶1.5。實(shí)驗(yàn)?zāi)P弯佋O(shè)與侏羅系煤層群開(kāi)挖情況如圖4所示。

    由圖4可知,15號(hào)煤層的開(kāi)采是在上部11,12及14號(hào)煤層采空區(qū)下進(jìn)行。實(shí)驗(yàn)觀測(cè)發(fā)現(xiàn)頂板有效失穩(wěn)塊體數(shù)在1~3塊間變化,而工作面開(kāi)挖采動(dòng)影響范圍內(nèi)的塊體數(shù)一般為1~6塊??梢?jiàn),工作面推進(jìn)過(guò)程中頂板失穩(wěn)變化率一般在16.7%~100%。

    實(shí)驗(yàn)中巖層應(yīng)力與位移的觀測(cè)由埋設(shè)于模型巖層內(nèi)部的壓力盒與位移計(jì)通過(guò)計(jì)算機(jī)自動(dòng)采集。實(shí)驗(yàn)采集設(shè)備埋設(shè)與連線、數(shù)據(jù)采集以及15號(hào)煤層開(kāi)挖前12,14,15號(hào)煤層頂板應(yīng)力分析結(jié)果如圖5所示。

    圖3 11~15號(hào)煤層柱狀圖Fig.3 Seam histogram of No.11-15 coal seam

    圖4 相似模型鋪設(shè)與侏羅系煤層群開(kāi)挖Fig.4 Similar model laying and Jurassic coal seam group excavation

    圖5 實(shí)驗(yàn)測(cè)點(diǎn)布置與多層頂板應(yīng)力分析結(jié)果Fig.5 Experimental measuring points and multi-layered roof stress analysis results

    由圖5(c)可知,15號(hào)煤層開(kāi)挖前,多采空區(qū)堅(jiān)硬厚層破斷頂板群結(jié)構(gòu)全長(zhǎng)范圍內(nèi)頂板應(yīng)力分布較為復(fù)雜,但相同位置處多層破斷頂板應(yīng)力大小與分布卻具有一定的相似性,近距離多分層頂板不同位置處的應(yīng)力基本趨于一穩(wěn)定值,且其范圍一般在0~14 MPa,從實(shí)驗(yàn)開(kāi)挖全程看,具有代表性的頂板應(yīng)力為7.5 MPa。

    對(duì)15號(hào)煤層開(kāi)挖過(guò)程中的頂板承載破斷特征進(jìn)行樣本取樣,得到堅(jiān)硬頂板失穩(wěn)條件下的樣本統(tǒng)計(jì)數(shù)據(jù),見(jiàn)表1。

    表1 15號(hào)煤層頂板失穩(wěn)時(shí)的樣本統(tǒng)計(jì)數(shù)據(jù)Table 1 No.15 coal seam roof instability sample data

    為了使實(shí)測(cè)數(shù)據(jù)具有較好的相關(guān)性,威布爾分布函數(shù)位置參數(shù)需滿足式(14)。根據(jù)實(shí)測(cè)數(shù)據(jù)(表1),聯(lián)立式(13),得到大同礦區(qū)砂巖頂板群結(jié)構(gòu)失穩(wěn)的線性相關(guān)系數(shù)與位置參數(shù)σ0的關(guān)系曲線(圖6)。

    圖6 失穩(wěn)相關(guān)系數(shù)與位置參數(shù)關(guān)系Fig.6 The relationship between instability correlation coefficient and positional parameters

    由圖6可知,在位置參數(shù)σ0=0時(shí),頂板群結(jié)構(gòu)的失穩(wěn)相關(guān)系數(shù)達(dá)到最大,此時(shí)大同礦區(qū)砂巖頂板群結(jié)構(gòu)的失穩(wěn)具有較好的線性關(guān)系,故位置參數(shù)σ0=0可視為頂板破斷失穩(wěn)的應(yīng)力閾值,其物理意義:只有在頂板破斷塊體承載量大于0的條件下,堅(jiān)硬厚層頂板群才將出現(xiàn)失穩(wěn)趨勢(shì)。

    根據(jù)大同礦區(qū)近距離砂巖頂板群結(jié)構(gòu)的承載與失穩(wěn)特征,對(duì)實(shí)驗(yàn)中堅(jiān)硬厚層頂板的受力及破斷失穩(wěn)率進(jìn)行觀測(cè),計(jì)算分析得到線性回歸方程系數(shù)分別為

    其中,a與b為近距離相似頂板條件下線性回歸方程系數(shù);sxy與sxx為根據(jù)表1所得樣本數(shù)據(jù)計(jì)算得到的換算樣本參數(shù)的協(xié)方差及方差,計(jì)算方法參照式(13);與為實(shí)驗(yàn)數(shù)據(jù)折算平均值。實(shí)驗(yàn)數(shù)據(jù)的折算平均值為

    考慮工作面正常推進(jìn)過(guò)程中,影響工作面回采的有效破斷塊體數(shù)目Ni一般為1~3[15-16],故堅(jiān)硬厚層頂板的單層破斷失穩(wěn)在威布爾分布條件下,根據(jù)表1中數(shù)據(jù),聯(lián)立式(20)與式(21),計(jì)算得到大同礦區(qū)近距離砂巖頂板群結(jié)構(gòu)失穩(wěn)的線性回歸方程系數(shù)分別為

    為檢驗(yàn)實(shí)驗(yàn)所得數(shù)據(jù)是否服從威布爾函數(shù)分布,即實(shí)測(cè)數(shù)據(jù)雙對(duì)數(shù)折算參數(shù)是否線性相關(guān),對(duì)實(shí)驗(yàn)采集數(shù)據(jù)(表1)進(jìn)行折算擬合(xi,yi),如圖7所示。

    由圖7可知,大同礦區(qū)砂巖堅(jiān)硬厚層頂板群結(jié)構(gòu)相似模擬實(shí)驗(yàn)數(shù)據(jù)的折算結(jié)果(xi,yi)具有較好的線性特征,說(shuō)明多采空區(qū)下頂板群結(jié)構(gòu)的失穩(wěn)形式基本服從威布爾分布。

    根據(jù)式(22),計(jì)算得到大同礦區(qū)砂巖堅(jiān)硬厚層頂板群結(jié)構(gòu)失穩(wěn)服從威布爾函數(shù)分布條件下的頂板巖層平均形狀參數(shù)及平均尺度參數(shù)分別為

    圖7 實(shí)測(cè)換算數(shù)據(jù)線性相關(guān)性Fig.7 Linear correlation of measured data

    至此,得到大同砂巖堅(jiān)硬厚層頂板群結(jié)構(gòu)的失穩(wěn)率表達(dá)式為

    根據(jù)式(24)分析得到頂板群結(jié)構(gòu)的失穩(wěn)率同頂板破斷塊體承載及下部煤層開(kāi)采影響范圍內(nèi)頂板破斷塊體數(shù)間的關(guān)系,如圖8所示。

    圖8 頂板群結(jié)構(gòu)失穩(wěn)曲線Fig.8 The roof group instability curves

    由圖8可知,相同破斷塊體數(shù)條件下,頂板群結(jié)構(gòu)的失穩(wěn)率隨頂板承載量的增加而增大;頂板相同承載條件下,頂板群結(jié)構(gòu)的失穩(wěn)率隨頂板破斷塊體數(shù)目的增加而增加。故大同礦區(qū)砂巖頂板群結(jié)構(gòu)的穩(wěn)定性與工作面來(lái)壓最終取決于頂板相應(yīng)承載條件下的塊體失穩(wěn)數(shù)。在頂板破斷塊體承載條件已知的情況下,判定堅(jiān)硬厚層頂板群結(jié)構(gòu)的失穩(wěn)率,只要知道當(dāng)層頂板在工作面開(kāi)采影響范圍內(nèi)的塊體失穩(wěn)數(shù),即可由圖8直接找到對(duì)應(yīng)應(yīng)力條件下的頂板失穩(wěn)率。

    以大同礦區(qū)12,14與15號(hào)近距離煤層開(kāi)采為例,當(dāng)工作面堅(jiān)硬厚層頂板層數(shù)n=3時(shí),依次將頂板破斷塊體結(jié)構(gòu)自上而下命名為Ⅲ,Ⅱ與Ⅰ號(hào)結(jié)構(gòu)。根據(jù)上述理論分析與相似模擬實(shí)測(cè)結(jié)果,計(jì)算得到15號(hào)煤層工作面開(kāi)采過(guò)程中頂板塊體結(jié)構(gòu)的失穩(wěn)特征,如圖9所示。

    圖9 15號(hào)煤層開(kāi)采時(shí)頂板塊體結(jié)構(gòu)的失穩(wěn)曲線Fig.9 The roof group instability curves in No.15 coal seam

    由圖9可知,15號(hào)煤層正常開(kāi)采過(guò)程中,相同應(yīng)力場(chǎng)條件下,靠近下部的堅(jiān)硬厚層頂板破斷塊體結(jié)構(gòu)的失穩(wěn)率要比其上部頂板塊體結(jié)構(gòu)的失穩(wěn)率高,即在相同頂板應(yīng)力條件下,堅(jiān)硬厚層頂板破斷塊體結(jié)構(gòu)的失穩(wěn)率滿足關(guān)系:P(Ⅰ號(hào))>P(Ⅱ號(hào))>P(Ⅲ號(hào))。還可以看出,破斷頂板群塊體結(jié)構(gòu)間同步失穩(wěn)率的大小隨著同步失穩(wěn)頂板群結(jié)構(gòu)數(shù)目的增加而有所降低,且單層頂板破斷塊體結(jié)構(gòu)的失穩(wěn)率要高于其組合頂板條件下的同步失穩(wěn)率。如大同礦區(qū)12,14與15號(hào)煤層堅(jiān)硬厚層頂板條件下,Ⅰ,Ⅱ與Ⅲ號(hào)塊體結(jié)構(gòu)間同步失穩(wěn)的幾率有:P(Ⅰ或Ⅱ號(hào))>P(Ⅰ,Ⅱ號(hào))>P(Ⅰ,Ⅱ和Ⅲ號(hào))。

    大同礦區(qū)侏羅系煤層頂板群破斷失穩(wěn)的相似模擬實(shí)驗(yàn)表明,近距離煤層破斷頂板群承載相似條件下(表1),15號(hào)煤層堅(jiān)硬厚層頂板Ⅰ號(hào)結(jié)構(gòu)隨工作面推進(jìn)發(fā)生單層失穩(wěn)的幾率約64%;14號(hào)煤層頂板Ⅱ號(hào)結(jié)構(gòu)單層失穩(wěn)率約50%;12號(hào)煤層頂板Ⅲ號(hào)結(jié)構(gòu)單層失穩(wěn)率約28%;單層頂板加權(quán)平均失穩(wěn)率約53%;Ⅰ與Ⅱ號(hào)頂板結(jié)構(gòu)發(fā)生同步失穩(wěn)的幾率約30%;Ⅰ,Ⅱ與Ⅲ號(hào)頂板結(jié)構(gòu)同步失穩(wěn)率僅8.9%左右;多層頂板加權(quán)平均失穩(wěn)率約17.3%。

    3 工作面礦壓現(xiàn)場(chǎng)實(shí)測(cè)分析

    大同礦區(qū)永定莊煤礦15號(hào)煤層厚度為2.2~4.6 m,煤層傾角1°~2°,距14號(hào)層間距為13~15 m,工作面裝備102架ZZS6 000/1.7/3.7型支撐掩護(hù)式液壓支架。支架額定初撐力5 105 kN,額定工作阻力6 000 kN。工作面長(zhǎng)150 m,沿工作面布置10個(gè)壓力分機(jī),采用KBJ60III型煤礦在線連續(xù)頂板動(dòng)態(tài)監(jiān)測(cè)系統(tǒng)進(jìn)行工作面支架阻力采集。

    鑒于近距離煤層多采空區(qū)堅(jiān)硬頂板破斷塊體結(jié)構(gòu)發(fā)生單層失穩(wěn)或多分層同步失穩(wěn)時(shí)采場(chǎng)礦壓顯現(xiàn)程度的不同,對(duì)工作面礦壓數(shù)據(jù)分析時(shí)需對(duì)實(shí)測(cè)數(shù)據(jù)進(jìn)行分區(qū)間處理。區(qū)間邊界以支架初撐力、支架額定工作阻力為上下限,對(duì)于區(qū)間內(nèi)數(shù)值接近而分布又較為集中的實(shí)測(cè)數(shù)據(jù)進(jìn)一步劃分為子區(qū)間。根據(jù)大同礦區(qū)永定莊煤礦15號(hào)煤層工作面礦壓實(shí)測(cè)數(shù)據(jù),對(duì)工作面端部和中部的兩組支架阻力觀測(cè)結(jié)果進(jìn)行分析,得到多采空區(qū)堅(jiān)硬厚層頂板條件下的工作面礦壓顯現(xiàn)特征,如圖10所示。

    圖10 工作面實(shí)測(cè)支架阻力區(qū)間分布Fig.10 Support resistance interval distribution in working face

    由圖10可知,永定莊煤礦15號(hào)煤層工作面來(lái)壓期間,頂板不同部位礦壓顯現(xiàn)特征:

    (1)端頭5號(hào)支架工作阻力位于3.9~5.7 MN區(qū)域的比率約占60%,最接近于15號(hào)煤層頂板發(fā)生單層失穩(wěn)率64%,由此可見(jiàn)工作面端頭礦壓顯現(xiàn)主要由上方煤層頂板的單層失穩(wěn)形式引起,原因在于端頭邊界效應(yīng)的影響導(dǎo)致多層上位頂板仍然受到邊界實(shí)體煤部分的支撐,從而引起上位采空區(qū)頂板的失穩(wěn)相對(duì)滯后;支架工作阻力位于5.7~6.0 MN區(qū)域的比率約占22%,接近于12號(hào)煤層頂板Ⅲ號(hào)結(jié)構(gòu)單層失穩(wěn)率28%,可認(rèn)為上位Ⅲ號(hào)頂板結(jié)構(gòu)受到擾動(dòng)處于失穩(wěn)狀態(tài)對(duì)下位工作面采場(chǎng)帶來(lái)一定影響,但此時(shí)多層頂板結(jié)構(gòu)保持單層逐次失穩(wěn)狀態(tài);工作面端頭支架工作阻力位于6.0~6.6 MN區(qū)域的比率約占18%,最接近多層頂板的加權(quán)平均失穩(wěn)率17.3%,由此可認(rèn)為造成工作面端頭來(lái)壓較高的主要影響因素來(lái)自多層頂板的同步失穩(wěn)。

    (2)中部55號(hào)支架工作阻力位于3.9~5.1 MN區(qū)域的比率約占55%,最接近煤層上覆頂板的單層加權(quán)平均失穩(wěn)幾率53%,由此可見(jiàn)工作面中部礦壓顯現(xiàn)主要由煤層頂板群結(jié)構(gòu)的單層失穩(wěn)形式引起;支架工作阻力位于5.1~6.0 MN區(qū)域的比率占18%,最接近多層頂板加權(quán)平均失穩(wěn)率17.3%,可認(rèn)為造成工作面中部來(lái)壓的主要影響因素為多層頂板的同步失穩(wěn);工作阻力位于6.0~6.6 MN區(qū)域的比率約占21%,接近12號(hào)煤層頂板單層失穩(wěn)率28%,可認(rèn)為上位Ⅲ號(hào)頂板結(jié)構(gòu)受到擾動(dòng)處于失穩(wěn)狀態(tài),導(dǎo)致工作面采場(chǎng)壓力變化。

    由此可見(jiàn),15號(hào)煤層頂板來(lái)壓與上覆多采空區(qū)破斷頂板群結(jié)構(gòu)的失穩(wěn)具有一定關(guān)聯(lián),上覆頂板破斷塊體的不同組合失穩(wěn)形式導(dǎo)致工作面回采過(guò)程中采場(chǎng)來(lái)壓的不同特征,而各種塊體的組合失穩(wěn)又具有一定的概率特征;實(shí)測(cè)支架工作阻力與推進(jìn)距離的關(guān)系曲線表明了多采空區(qū)下工作面頂板來(lái)壓的頻繁性與復(fù)雜性,說(shuō)明多采空區(qū)多分層堅(jiān)硬頂板條件下的采場(chǎng)來(lái)壓具有多樣性特征。

    4 結(jié) 論

    (1)多采空區(qū)下堅(jiān)硬厚層頂板群結(jié)構(gòu)的失穩(wěn)呈現(xiàn)出一定的概率特征,可采用威布爾多參數(shù)分布函數(shù)對(duì)頂板破斷群結(jié)構(gòu)失穩(wěn)進(jìn)行分析。通過(guò)相似模擬實(shí)驗(yàn)結(jié)果對(duì)分布函數(shù)參數(shù)的確定,得到了堅(jiān)硬多分層頂板的破斷失穩(wěn)率,可較好地分析多采空區(qū)破斷頂板群結(jié)構(gòu)的穩(wěn)定性及下部工作面采場(chǎng)礦壓顯現(xiàn)特征。

    (2)近距離煤層頂板群結(jié)構(gòu)的失穩(wěn)率取決于工作面采動(dòng)影響范圍內(nèi)的頂板破斷塊體數(shù)。通常,下部煤層開(kāi)采影響范圍內(nèi)的單層頂板塊體數(shù)目越多,頂板發(fā)生單層失穩(wěn)的幾率也越高,而多層頂板發(fā)生同步失穩(wěn)的幾率一般小于其單層失穩(wěn)率。

    (3)工作面來(lái)壓期間支架工作阻力的頻率分布與多層破斷頂板的失穩(wěn)率具有較好的一致性,表明了多采空區(qū)堅(jiān)硬破斷頂板群結(jié)構(gòu)的失穩(wěn)將造成工作面采場(chǎng)來(lái)壓的多樣性,而且單層頂板失穩(wěn)或多層頂板同步失穩(wěn)形式將影響工作面礦壓顯現(xiàn)強(qiáng)度。

    [1] 于 斌.大同礦區(qū)綜采40 a開(kāi)采技術(shù)研究[J].煤炭學(xué)報(bào),2010, 35(11):1772-1777.

    Yu Bin.Study on fully mechanized coa1 mining technology in passed 40 years in Datong mining area[J].Journal of China Coal Society, 2010,35(11):1772-1777.

    [2] 黃晉兵.王臺(tái)礦井綜采工作面處理堅(jiān)硬頂板技術(shù)[J].煤炭科學(xué)技術(shù),2008,36(9):4-7.

    Huang Jinbing.Hard roof treatment technology for fully mechanized longwall coal mining face in Wangtai mine[J].Coal Science and Technology,2008,36(9):4-7.

    [3] 王金安,尚新春,劉 紅,等.采空區(qū)堅(jiān)硬頂板破斷機(jī)理與災(zāi)變塌陷研究[J].煤炭學(xué)報(bào),2008,33(8):850-855.

    Wang Jin’an,Shang Xinchun,Liu Hong,et al.Study on fracture mechanism and catastrophic collapse of strong roof strata above the mined area[J].Journal of China Coal Society,2008,33(8):850-855.

    [4] 于 斌.“兩硬”條件下4.5~5 m大采高綜采技術(shù)[J].煤炭科學(xué)技術(shù),2004,32(8):35-37.

    Yu Bin.4.5~5 m high cutting fully mechanized mining technology in seam with hard seam and hard roof[J].Coal Science and Technology,2004,32(8):35-37.

    [5] 朱德仁,錢鳴高,徐林生.堅(jiān)硬頂板來(lái)壓控制的探討[J].煤炭學(xué)報(bào),1991,16(2):11-20.

    Zhu Deren,Qian Minggao,Xu Linsheng.Discussion on control of hard roof weighting[J].Journal of China Coal Society,1991,16 (2):11-20.

    [6] 朱衛(wèi)兵.淺埋近距離煤層重復(fù)采動(dòng)關(guān)鍵層結(jié)構(gòu)失穩(wěn)機(jī)理研究[D].徐州:中國(guó)礦業(yè)大學(xué),2010.

    [7] 靳鐘銘,徐林生.煤礦堅(jiān)硬頂板控制[M].北京:煤炭工業(yè)出版社,1994.

    [8] 韓貴雷,韓立軍,王延寧,等.厚層復(fù)合石灰?guī)r頂板垮落特性試驗(yàn)研究[J].采礦與安全工程學(xué)報(bào),2007,24(1):60-64.

    Han Guilei,Han Lijun,Wang Yanning,et al.Experimental study of collapsing characteristic of thick complex lime-stone roof[J].Journal of Mining&Safety Engineering,2007,24(1):60-64.

    [9] 楊 文,王玉民.大同礦區(qū)采空區(qū)覆巖破壞特點(diǎn)分析[J].煤炭技術(shù),2005,24(12):107-108.

    Yang Wen,Wang Yumin.Analysis on the empty area cover breakage characteristics of Datong Coal Mine[J].Coal Technology,2005,24 (12):107-108.

    [10] 張春山.永定莊煤礦多采空區(qū)下煤層頂板礦壓實(shí)測(cè)研究[J].同煤科技,2011,129(3):15-18.

    Zhang Chunshan.Roof mine pressure research of coal seam below multi gob in Yongdingzhuang Coal Mine[J].Datong Coal Mining and Technology,2011,129(3):15-18.

    [11] 趙延林,王衛(wèi)軍,趙伏軍,等.多層采空區(qū)隔離頂板安全性強(qiáng)度折減法[J].煤炭學(xué)報(bào),2010,35(8):1257-1262.

    Zhao Yanlin,Wang Weijun,Zhao Fujun,et al.Strength reduction method to study safety of multilayer goafs isolation roof[J].Journal of China Coal Society,2010,35(8):1257-1262.

    [12] 姚邦華,周海峰,陳 龍.重復(fù)采動(dòng)下覆巖裂隙發(fā)育規(guī)律模擬研究[J].采礦與安全工程學(xué)報(bào),2010,27(3):443-446.

    Yao Banghua,Zhou Haifeng,Chen Long.Numerical simulation about fracture development in overlying rocks under repeated mining[J].Journal of Mining&Safety Engineering,2010,27(3): 443-446.

    [13] 張耀平,曹 平,袁海平,等.復(fù)雜采空區(qū)穩(wěn)定性數(shù)值模擬分析[J].采礦與安全工程學(xué)報(bào),2010,27(2):233-238.

    Zhang Yaoping,Cao Ping,Yuan Haiping,et al.Numerical simulation on stability of complicated goaf[J].Journal of Mining&Safety Engineering,2010,27(2):233-238.

    [14] 錢鳴高,繆協(xié)興,何富連.采場(chǎng)“砌體梁”結(jié)構(gòu)的關(guān)鍵塊分析[J].煤炭學(xué)報(bào),1994,19(6):557-563.

    Qian Minggao,Miao Xiexing,He Fulian.Analysis of key blocks in the structure of voussoir in longwall mining[J].Journal of China Coal Society,1994,19(6):557-563.

    [15] 錢鳴高,繆協(xié)興,許家林,等.巖層控制的關(guān)鍵層理論[M].徐州中國(guó)礦業(yè)大學(xué)出版社,2000.

    [16] 錢鳴高,石平五,鄒喜正,等.礦山壓力與巖層控制[M].徐州:中國(guó)礦業(yè)大學(xué)出版社,2003.

    [17] 安偉光,蔡蔭林,陳衛(wèi)東.隨機(jī)結(jié)構(gòu)系統(tǒng)可靠性分析與優(yōu)化[M].哈爾濱:哈爾濱工程大學(xué)出版社,2007.

    [18] 徐以鋒,翁永剛,王明星,等.威布爾分布三參數(shù)的擬合特性[J].鄭州大學(xué)學(xué)報(bào)(理學(xué)版),2003,35(3):30-35.

    Xu Yifeng,Weng Yonggang,Wang Mingxing,et al.Characteristics of three fitted parameters for Weibull distribution[J].Journal of Zhengzhou University(Natural Science Edition),2003,35(3):30-35.

    [19] 高 峰,謝和平,趙 鵬.Weibull模量和巖石強(qiáng)度的分形性質(zhì)[J].科學(xué)通報(bào),1993,38(15):1435-1438.

    Gao Feng,Xie Heping,Zhao Peng.The fractal nature of the Weibull modulus and rock strength[J].Chines Science Bulletin,1993,38 (15):1435-1438.

    Destabilization regularity of hard thick roof group under the multi gob

    LIU Chang-you1,2,YANG Jing-xuan1,2,YU Bin1,3,YANG Pei-ju1,2

    (1.School of Mines,China University of Mining&Technology,Xuzhou 221116,China;2.State Key Laboratory of Coal Resources and Safe Mining,China University of Mining&Technology,Xuzhou 221116,China;3.Datong Coal Mine Group Company,Datong 037003,China)

    In order to understand the coal face pressure effects from the hard roof group structure breaking and instability,the relation was proved between the hard roof group structure and strata behavior.It investigated the instability characteristics of roof group structure combining the research methods of the theoretical analysis,similar simulation experiment and field measured according to the conditions of hard thick roof group structure under multi gob in Datong coal mining area.It shows that:the roof group structure instability has some certain probability characteristic;multiple parameters Weibull function describes well on instability forms in breaking roof group structure,and establishes the similar roof structure instability rate under the condition of the close distance coal seam group,then gives the roof instability parameters,and analyzes on the goaf roof instability rate in Jurassic coal seam group.It is shown that the roof group structure instability rate have some influence on the face support resistance.Through the mine pressure and support resistance measurement results on the end and central position of the coal face,verify the hard thick roof instability regularity in the No.15 face of Yongding Village Coal Mine in Datong mining area.

    multi gob;breaking roof group;structure instability;ground pressure behavior;Weibull distribution

    TD821;TD322

    A

    0253-9993(2014)03-0395-09

    劉長(zhǎng)友,楊敬軒,于 斌,等.多采空區(qū)下堅(jiān)硬厚層破斷頂板群結(jié)構(gòu)的失穩(wěn)規(guī)律[J].煤炭學(xué)報(bào),2014,39(3):395-403.

    10.13225/j.cnki.jccs.2013.0091

    Liu Changyou,Yang Jingxuan,Yu Bin,et al.Destabilization regularity of hard thick roof group under the multi gob[J].Journal of China Coal Society,2014,39(3):395-403.doi:10.13225/j.cnki.jccs.2013.0091

    2013-01-21 責(zé)任編輯:張曉寧

    國(guó)家自然科學(xué)基金資助項(xiàng)目(51174192);江蘇省“333”培養(yǎng)基金資助項(xiàng)目(BRA2010024);江蘇省研究生培養(yǎng)創(chuàng)新工程資助項(xiàng)目(CXLX12_0964)

    劉長(zhǎng)友(1965—),男,山東東營(yíng)人,教授,博士生導(dǎo)師,博士。E-mail:cyliu@cumt.edu.cn

    猜你喜歡
    厚層礦壓塊體
    一種新型單層人工塊體Crablock 的工程應(yīng)用
    溫州地區(qū)厚層軟土UU三軸試驗(yàn)與FLAC3D數(shù)值模擬
    Winter Walking in Wales威爾士的冬季徒步旅行
    采煤工作面礦壓顯現(xiàn)規(guī)律影響因素分析
    厚層塊狀特低滲礫巖油藏水平井壓裂參數(shù)優(yōu)化
    綜放工作面礦壓顯現(xiàn)規(guī)律研究
    山西煤炭(2015年4期)2015-12-20 11:36:18
    一種Zr 基塊體金屬玻璃的納米壓入蠕變行為研究
    上海金屬(2015年3期)2015-11-29 01:09:58
    塊體非晶合金及其應(yīng)用
    哮天龍水庫(kù)大壩厚層碾壓混凝土技術(shù)研究
    波浪作用下斜坡上護(hù)面塊體斷裂破壞的數(shù)值模擬
    水道港口(2015年1期)2015-02-06 01:25:35
    亚洲国产成人一精品久久久| 免费观看a级毛片全部| av卡一久久| 久久99精品国语久久久| 久久精品国产a三级三级三级| 搡老乐熟女国产| 色网站视频免费| 国产伦理片在线播放av一区| 中文字幕最新亚洲高清| 十八禁高潮呻吟视频| 成人三级做爰电影| 久久精品人人爽人人爽视色| 欧美精品亚洲一区二区| 欧美日韩亚洲国产一区二区在线观看 | 午夜福利,免费看| 国产老妇伦熟女老妇高清| 国产 一区精品| 亚洲av成人不卡在线观看播放网 | 久久人妻熟女aⅴ| 国产精品三级大全| 成人亚洲欧美一区二区av| 亚洲精品中文字幕在线视频| 欧美另类一区| 亚洲av电影在线进入| 国产精品香港三级国产av潘金莲 | 1024视频免费在线观看| 天天躁狠狠躁夜夜躁狠狠躁| av福利片在线| 男人添女人高潮全过程视频| 老司机亚洲免费影院| 街头女战士在线观看网站| 国产日韩欧美在线精品| 国产精品久久久久成人av| 婷婷色麻豆天堂久久| 亚洲国产成人一精品久久久| 亚洲欧美一区二区三区国产| 大码成人一级视频| 看免费成人av毛片| 国产99久久九九免费精品| 在线天堂最新版资源| 王馨瑶露胸无遮挡在线观看| 国产 一区精品| 捣出白浆h1v1| 精品少妇一区二区三区视频日本电影 | 男人操女人黄网站| 久久精品aⅴ一区二区三区四区| 各种免费的搞黄视频| 国产av码专区亚洲av| 丝袜喷水一区| 在线免费观看不下载黄p国产| 国产视频首页在线观看| 汤姆久久久久久久影院中文字幕| 国产精品国产av在线观看| 色视频在线一区二区三区| 美女脱内裤让男人舔精品视频| 欧美人与性动交α欧美软件| 久久久久久久久久久久大奶| 国产一级毛片在线| 最新在线观看一区二区三区 | 视频区图区小说| 中文字幕人妻丝袜一区二区 | 国产免费视频播放在线视频| 精品国产乱码久久久久久小说| 汤姆久久久久久久影院中文字幕| 久久99精品国语久久久| 日韩人妻精品一区2区三区| 国产 一区精品| 午夜福利网站1000一区二区三区| 满18在线观看网站| 咕卡用的链子| 纵有疾风起免费观看全集完整版| 国产精品一国产av| 欧美精品av麻豆av| 中国国产av一级| 下体分泌物呈黄色| 嫩草影视91久久| 久久国产精品男人的天堂亚洲| 99re6热这里在线精品视频| 国产成人精品久久久久久| 在线观看免费日韩欧美大片| 亚洲免费av在线视频| av一本久久久久| 精品国产国语对白av| 少妇 在线观看| 国产亚洲av片在线观看秒播厂| 熟女av电影| 国产av码专区亚洲av| 国产精品无大码| 午夜激情久久久久久久| 亚洲精品av麻豆狂野| 国产片特级美女逼逼视频| 日韩一区二区三区影片| 亚洲美女搞黄在线观看| 日韩电影二区| 好男人视频免费观看在线| 最近手机中文字幕大全| 午夜免费鲁丝| 国产成人免费无遮挡视频| 好男人视频免费观看在线| 亚洲一级一片aⅴ在线观看| 又粗又硬又长又爽又黄的视频| 2018国产大陆天天弄谢| 久久久久人妻精品一区果冻| 黄色毛片三级朝国网站| 高清欧美精品videossex| 久久久久精品性色| 精品国产超薄肉色丝袜足j| 老熟女久久久| 久久国产精品男人的天堂亚洲| 女人高潮潮喷娇喘18禁视频| 国产色婷婷99| 国产又色又爽无遮挡免| kizo精华| 国产一区二区激情短视频 | 精品免费久久久久久久清纯 | 男女下面插进去视频免费观看| 国产精品麻豆人妻色哟哟久久| 午夜免费鲁丝| 精品福利永久在线观看| 岛国毛片在线播放| 亚洲成人免费av在线播放| 国产精品99久久99久久久不卡 | 黑人欧美特级aaaaaa片| 国产国语露脸激情在线看| 制服丝袜香蕉在线| 在线亚洲精品国产二区图片欧美| 日日爽夜夜爽网站| 只有这里有精品99| 男女午夜视频在线观看| 韩国高清视频一区二区三区| 免费人妻精品一区二区三区视频| 免费少妇av软件| www日本在线高清视频| 国产片内射在线| 久久久亚洲精品成人影院| 日韩,欧美,国产一区二区三区| 成年动漫av网址| 亚洲国产精品成人久久小说| 国产精品女同一区二区软件| www.精华液| 中国三级夫妇交换| 亚洲精品自拍成人| 99国产精品免费福利视频| 欧美人与性动交α欧美软件| 成人黄色视频免费在线看| 亚洲国产精品999| 国产亚洲av高清不卡| 国产精品欧美亚洲77777| 欧美日韩福利视频一区二区| 伦理电影大哥的女人| 韩国av在线不卡| 国产激情久久老熟女| 母亲3免费完整高清在线观看| 丝袜喷水一区| 在线观看www视频免费| 最近2019中文字幕mv第一页| 18禁动态无遮挡网站| 久久青草综合色| 少妇的丰满在线观看| 色婷婷av一区二区三区视频| 视频区图区小说| 伊人久久大香线蕉亚洲五| 中文字幕色久视频| 亚洲第一区二区三区不卡| 午夜福利视频在线观看免费| 久久久久精品性色| 国产视频首页在线观看| e午夜精品久久久久久久| 99久久精品国产亚洲精品| 别揉我奶头~嗯~啊~动态视频 | av卡一久久| 精品国产国语对白av| 国产精品久久久久久精品电影小说| 亚洲精品第二区| 亚洲一区中文字幕在线| 日本欧美国产在线视频| 亚洲国产欧美日韩在线播放| 日韩精品免费视频一区二区三区| 久久久久久人人人人人| 久久精品久久久久久噜噜老黄| 伦理电影大哥的女人| 亚洲第一av免费看| 一边摸一边做爽爽视频免费| 国产av国产精品国产| 香蕉国产在线看| 男的添女的下面高潮视频| 深夜精品福利| 国产xxxxx性猛交| 欧美国产精品va在线观看不卡| 久久影院123| 亚洲欧美成人综合另类久久久| 欧美日韩综合久久久久久| 嫩草影视91久久| 免费黄色在线免费观看| 国产97色在线日韩免费| 久久精品亚洲av国产电影网| 久久性视频一级片| 纵有疾风起免费观看全集完整版| 亚洲成色77777| 男女免费视频国产| 久久久久久久国产电影| 久久久久久久精品精品| 十分钟在线观看高清视频www| 欧美精品一区二区免费开放| 热99久久久久精品小说推荐| 欧美av亚洲av综合av国产av | 热re99久久国产66热| 国产成人免费无遮挡视频| 青青草视频在线视频观看| 黄色怎么调成土黄色| 天美传媒精品一区二区| 精品人妻熟女毛片av久久网站| 精品国产一区二区三区四区第35| 建设人人有责人人尽责人人享有的| 看免费av毛片| 精品午夜福利在线看| 日本色播在线视频| 男女免费视频国产| 久久精品aⅴ一区二区三区四区| 女的被弄到高潮叫床怎么办| 王馨瑶露胸无遮挡在线观看| 久久韩国三级中文字幕| 欧美少妇被猛烈插入视频| 国产成人免费观看mmmm| 国产97色在线日韩免费| 街头女战士在线观看网站| 男的添女的下面高潮视频| 久久女婷五月综合色啪小说| 欧美在线黄色| 亚洲美女搞黄在线观看| 国产乱来视频区| 亚洲五月色婷婷综合| 又黄又粗又硬又大视频| 91精品伊人久久大香线蕉| 亚洲情色 制服丝袜| 国产男女内射视频| 日本猛色少妇xxxxx猛交久久| 桃花免费在线播放| 一区二区av电影网| 看非洲黑人一级黄片| 看非洲黑人一级黄片| 国产欧美日韩综合在线一区二区| 亚洲av中文av极速乱| 青草久久国产| 国产成人精品久久二区二区91 | 国产在视频线精品| 好男人视频免费观看在线| 老司机深夜福利视频在线观看 | 久久免费观看电影| 国产一区二区三区综合在线观看| 少妇精品久久久久久久| 国产又色又爽无遮挡免| videos熟女内射| 色视频在线一区二区三区| 在线观看人妻少妇| 18禁观看日本| 99精品久久久久人妻精品| 天天躁夜夜躁狠狠久久av| 大话2 男鬼变身卡| 美女高潮到喷水免费观看| 免费女性裸体啪啪无遮挡网站| 国产熟女欧美一区二区| 欧美精品一区二区大全| 午夜av观看不卡| 色吧在线观看| 国产精品欧美亚洲77777| 午夜福利视频精品| 丝袜美足系列| 不卡av一区二区三区| 亚洲欧美精品综合一区二区三区| 色婷婷久久久亚洲欧美| 亚洲熟女毛片儿| 久久毛片免费看一区二区三区| 香蕉丝袜av| 一本—道久久a久久精品蜜桃钙片| 99九九在线精品视频| 久久久久久人妻| 午夜福利一区二区在线看| 久久精品国产亚洲av高清一级| 亚洲av福利一区| 无限看片的www在线观看| 午夜福利视频在线观看免费| 免费女性裸体啪啪无遮挡网站| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲av国产av综合av卡| 国产成人免费观看mmmm| 99久久99久久久精品蜜桃| 中文字幕精品免费在线观看视频| 高清在线视频一区二区三区| 亚洲av电影在线进入| 伊人亚洲综合成人网| 各种免费的搞黄视频| 免费观看性生交大片5| 国产精品熟女久久久久浪| 免费黄频网站在线观看国产| 捣出白浆h1v1| 伊人亚洲综合成人网| 久久综合国产亚洲精品| 日韩制服丝袜自拍偷拍| 国产成人精品无人区| 中文字幕人妻丝袜制服| 国产成人啪精品午夜网站| 黑人欧美特级aaaaaa片| 乱人伦中国视频| 国产有黄有色有爽视频| 亚洲伊人久久精品综合| 人人澡人人妻人| 丰满饥渴人妻一区二区三| 国产成人精品在线电影| 国产国语露脸激情在线看| 99久久99久久久精品蜜桃| 久久久国产欧美日韩av| 狠狠婷婷综合久久久久久88av| 国产精品人妻久久久影院| 国产成人av激情在线播放| 欧美日韩亚洲高清精品| 国产免费福利视频在线观看| 久久午夜综合久久蜜桃| 麻豆乱淫一区二区| 波野结衣二区三区在线| 亚洲熟女精品中文字幕| 国产av国产精品国产| 极品少妇高潮喷水抽搐| 五月开心婷婷网| 精品国产乱码久久久久久小说| 欧美日韩视频高清一区二区三区二| 免费黄色在线免费观看| 久久婷婷青草| 国产熟女午夜一区二区三区| 日韩成人av中文字幕在线观看| 丝瓜视频免费看黄片| 在线天堂中文资源库| 日本一区二区免费在线视频| 国产1区2区3区精品| 老司机靠b影院| 日韩熟女老妇一区二区性免费视频| 国产精品亚洲av一区麻豆 | 欧美精品一区二区大全| 青春草亚洲视频在线观看| 丝袜在线中文字幕| 成年人免费黄色播放视频| 国产国语露脸激情在线看| 色婷婷久久久亚洲欧美| 午夜福利,免费看| 色婷婷av一区二区三区视频| 下体分泌物呈黄色| av视频免费观看在线观看| 欧美人与性动交α欧美精品济南到| √禁漫天堂资源中文www| 国产精品99久久99久久久不卡 | 色婷婷av一区二区三区视频| 又大又黄又爽视频免费| 丰满迷人的少妇在线观看| 悠悠久久av| 欧美最新免费一区二区三区| av.在线天堂| 男人操女人黄网站| 精品国产国语对白av| 免费少妇av软件| 久久免费观看电影| 精品少妇久久久久久888优播| 丁香六月天网| 又大又爽又粗| 美女福利国产在线| 久久精品熟女亚洲av麻豆精品| 久久精品亚洲av国产电影网| 亚洲精品国产av成人精品| 亚洲视频免费观看视频| 日韩不卡一区二区三区视频在线| a级毛片黄视频| 亚洲av日韩精品久久久久久密 | 国产女主播在线喷水免费视频网站| 欧美久久黑人一区二区| 十八禁高潮呻吟视频| 免费高清在线观看视频在线观看| 欧美日韩国产mv在线观看视频| 一区二区日韩欧美中文字幕| 久久久久精品久久久久真实原创| 国产精品嫩草影院av在线观看| 亚洲成av片中文字幕在线观看| 久久久久国产一级毛片高清牌| 热re99久久精品国产66热6| 黑丝袜美女国产一区| 欧美另类一区| 伊人久久大香线蕉亚洲五| 午夜日本视频在线| 亚洲成人手机| 少妇被粗大的猛进出69影院| 午夜老司机福利片| 一本久久精品| 亚洲精品久久成人aⅴ小说| 91aial.com中文字幕在线观看| 亚洲国产最新在线播放| avwww免费| 日本一区二区免费在线视频| 最近2019中文字幕mv第一页| 性高湖久久久久久久久免费观看| av网站在线播放免费| 国产成人欧美| 亚洲av欧美aⅴ国产| 黄色怎么调成土黄色| 欧美精品高潮呻吟av久久| 亚洲人成电影观看| 国产免费一区二区三区四区乱码| 久久天堂一区二区三区四区| 国产熟女午夜一区二区三区| 曰老女人黄片| 国产成人精品福利久久| 18禁裸乳无遮挡动漫免费视频| 韩国av在线不卡| 最近中文字幕2019免费版| 国产高清国产精品国产三级| 黄色 视频免费看| 国产成人精品无人区| 在线观看免费高清a一片| 亚洲精华国产精华液的使用体验| 精品国产露脸久久av麻豆| 在线看a的网站| 国产精品免费大片| 久久精品久久久久久久性| 免费观看av网站的网址| 亚洲,一卡二卡三卡| 一区福利在线观看| 日韩制服骚丝袜av| 国产精品99久久99久久久不卡 | 一边亲一边摸免费视频| 男人舔女人的私密视频| 女的被弄到高潮叫床怎么办| 亚洲欧美激情在线| 免费在线观看完整版高清| 少妇人妻久久综合中文| 啦啦啦啦在线视频资源| 久久久久久久久久久久大奶| 日本爱情动作片www.在线观看| 极品少妇高潮喷水抽搐| 午夜老司机福利片| 9热在线视频观看99| 国产精品二区激情视频| 成年女人毛片免费观看观看9 | 午夜日本视频在线| 美国免费a级毛片| 欧美黑人欧美精品刺激| 熟女少妇亚洲综合色aaa.| 亚洲国产精品国产精品| 欧美日韩成人在线一区二区| 亚洲精品国产av成人精品| 桃花免费在线播放| 一本—道久久a久久精品蜜桃钙片| 欧美日韩精品网址| 欧美最新免费一区二区三区| 国产成人av激情在线播放| 精品亚洲成国产av| 大香蕉久久成人网| 国产免费视频播放在线视频| 18禁国产床啪视频网站| 色综合欧美亚洲国产小说| 999久久久国产精品视频| 国产成人精品久久二区二区91 | 婷婷色综合大香蕉| 侵犯人妻中文字幕一二三四区| 精品酒店卫生间| 国产成人精品无人区| 咕卡用的链子| 女性生殖器流出的白浆| xxxhd国产人妻xxx| 国产不卡av网站在线观看| 欧美国产精品一级二级三级| 国产xxxxx性猛交| 中文字幕色久视频| 高清在线视频一区二区三区| 久久精品久久久久久噜噜老黄| 久久久久久免费高清国产稀缺| 伊人亚洲综合成人网| 秋霞在线观看毛片| 免费女性裸体啪啪无遮挡网站| 日韩欧美精品免费久久| 国产精品av久久久久免费| 精品少妇内射三级| 成人免费观看视频高清| 色播在线永久视频| 亚洲中文av在线| 亚洲av电影在线观看一区二区三区| 午夜福利乱码中文字幕| 日韩不卡一区二区三区视频在线| 一区二区三区激情视频| 亚洲人成电影观看| 麻豆精品久久久久久蜜桃| 天堂中文最新版在线下载| 日日摸夜夜添夜夜爱| 大香蕉久久网| 成人国语在线视频| 一边摸一边抽搐一进一出视频| 叶爱在线成人免费视频播放| av在线播放精品| 大香蕉久久网| 成人18禁高潮啪啪吃奶动态图| 尾随美女入室| 国产一区二区 视频在线| 午夜免费观看性视频| 少妇的丰满在线观看| 99久久人妻综合| 国产极品粉嫩免费观看在线| 日日啪夜夜爽| 日日撸夜夜添| 18禁观看日本| 日韩精品有码人妻一区| 亚洲av在线观看美女高潮| 美女主播在线视频| 老汉色∧v一级毛片| 麻豆av在线久日| 黑人欧美特级aaaaaa片| 美女视频免费永久观看网站| 肉色欧美久久久久久久蜜桃| 性少妇av在线| 如何舔出高潮| 伊人亚洲综合成人网| 亚洲综合精品二区| 国产成人一区二区在线| 大片电影免费在线观看免费| 午夜精品国产一区二区电影| 婷婷色麻豆天堂久久| 老司机影院成人| 久久精品国产亚洲av涩爱| 亚洲成国产人片在线观看| 久久久久久免费高清国产稀缺| 男女下面插进去视频免费观看| 伦理电影免费视频| 亚洲欧美色中文字幕在线| 亚洲一区二区三区欧美精品| 王馨瑶露胸无遮挡在线观看| 日韩视频在线欧美| 国产成人91sexporn| 黄网站色视频无遮挡免费观看| 国产精品久久久久久人妻精品电影 | 美女视频免费永久观看网站| 午夜av观看不卡| 如何舔出高潮| 久久人妻熟女aⅴ| 天堂8中文在线网| 丝瓜视频免费看黄片| 国产国语露脸激情在线看| 在线免费观看不下载黄p国产| 国产女主播在线喷水免费视频网站| 精品久久久精品久久久| 国产成人午夜福利电影在线观看| 精品一品国产午夜福利视频| 成人漫画全彩无遮挡| 人人妻人人爽人人添夜夜欢视频| 可以免费在线观看a视频的电影网站 | 国产在线一区二区三区精| 日韩精品有码人妻一区| 18禁观看日本| 日韩一卡2卡3卡4卡2021年| 久久人人97超碰香蕉20202| 亚洲av男天堂| 男人添女人高潮全过程视频| 黑人猛操日本美女一级片| 人人妻人人爽人人添夜夜欢视频| 综合色丁香网| 久久久久精品人妻al黑| 91老司机精品| 日本黄色日本黄色录像| 亚洲,一卡二卡三卡| 天天添夜夜摸| 免费观看性生交大片5| 日韩熟女老妇一区二区性免费视频| 亚洲国产欧美网| 免费女性裸体啪啪无遮挡网站| 无限看片的www在线观看| 国产精品久久久久久人妻精品电影 | 久久久精品免费免费高清| 一区二区av电影网| 国产精品蜜桃在线观看| 久久这里只有精品19| 久久精品国产a三级三级三级| 最近最新中文字幕免费大全7| 国产一区有黄有色的免费视频| 2018国产大陆天天弄谢| 国产毛片在线视频| 久久97久久精品| 欧美人与善性xxx| 午夜日韩欧美国产| 这个男人来自地球电影免费观看 | 新久久久久国产一级毛片| 制服诱惑二区| 亚洲欧美一区二区三区国产| 高清黄色对白视频在线免费看| 精品少妇一区二区三区视频日本电影 | 香蕉国产在线看| 国产乱来视频区| 国产精品成人在线| 久久久久精品性色| 电影成人av| 蜜桃国产av成人99| 最新在线观看一区二区三区 | 超碰97精品在线观看| 人人妻人人澡人人看| 涩涩av久久男人的天堂| 亚洲国产欧美在线一区| 国产午夜精品一二区理论片| 精品一区二区三区四区五区乱码 | 国产亚洲欧美精品永久| 国产亚洲av片在线观看秒播厂| 啦啦啦在线观看免费高清www| 涩涩av久久男人的天堂| 久久久久视频综合| 午夜福利一区二区在线看| 成人影院久久| 亚洲,欧美精品.| 国产成人a∨麻豆精品| 成人亚洲欧美一区二区av| 在线 av 中文字幕| 波多野结衣一区麻豆| 97在线人人人人妻| 日韩成人av中文字幕在线观看|