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

    多層釩礦體開采應(yīng)力遷移規(guī)律和分布研究

    2023-10-19 08:32:42池秀文董浩然汪宗英王其洲范純超
    金屬礦山 2023年9期
    關(guān)鍵詞:礦柱礦層下層

    池秀文 董浩然 汪宗英 王其洲 趙 龍 范純超

    (1.礦物資源加工與環(huán)境湖北省重點實驗室,湖北 武漢 430070;2.武漢理工大學(xué)資源與環(huán)境工程學(xué)院,湖北 武漢 430070;3.山東黃金礦業(yè)股份有限公司,山東 濟南 250100)

    釩元素因其優(yōu)異的物理和化學(xué)特性,在提高金屬材料物理力學(xué)性能、改善非金屬材料的化學(xué)性質(zhì)以及作為蓄電池材料等方面得到了廣泛應(yīng)用。但釩礦資源復(fù)雜的賦存條件和針對這類賦存條件不成熟的開采技術(shù)導(dǎo)致我國釩礦生產(chǎn)能力小,安全性也較差。釩礦資源多賦存于含釩沉積頁巖中,常以緩傾斜—傾斜、薄—中厚礦體產(chǎn)出,一般為多層排列分布,是一種典型的多層礦體,并且各礦層的埋深、走向、厚度和品位變化大[1]。開采近距離的多層礦體時,由于礦體層間距離較小,相鄰礦層開采的相互影響顯著,采場圍巖結(jié)構(gòu)和多次擾動后的次生應(yīng)力場變化劇烈,導(dǎo)致相鄰待采礦體的開采環(huán)境趨向復(fù)雜,形成了多層礦體開采時特殊的地壓環(huán)境,增加了開采難度。

    針對上述問題,許多學(xué)者在開采順序、開采方向等方面開展了深入研究,如葉義成等[2]通過相似模擬試驗,分析了不同開采順序多層礦采空區(qū)圍巖的變形規(guī)律,得出多層礦適合采用前進式開采方式。談曉明等[3]經(jīng)過現(xiàn)場研究,對某多層礦體礦層回采順序、各礦塊的超前關(guān)系、礦房布置形式等進行了優(yōu)化分析。殷亮[4]根據(jù)某多層薄礦脈近似水平的分布特點,利用數(shù)值模擬得到合理的開采順序。胡倩等[5]根據(jù)不同的圍巖移動角、礦層傾角與礦層間距,確定相鄰多層礦同階段開采時的超前距離與回采順序。池秀文等[6-7]基于彈性梁和薄板模型,求解了多層礦頂板的應(yīng)力分量表達式,得到頂板的極限跨距和破斷距。在采用相似模擬、數(shù)值模擬和現(xiàn)場試驗等方法分析多層礦開采穩(wěn)定性時,一般利用模擬或試驗的應(yīng)力和應(yīng)變結(jié)果,判斷圍巖破壞情況,據(jù)此評價開采方案的合理性。此外,不少學(xué)者還針對多層礦體的開采方法開展研究并得到了不少有益結(jié)論,如林坤峰等[8]提出大盤區(qū)機械化分層充填沿傾向后退式開采和盤區(qū)房柱法兩種多層礦開采方法。古德生等[9]提出了深孔合采井下礦廢分離連續(xù)采礦法,避免了多層礦間柱回采帶來的一系列技術(shù)經(jīng)濟難題。孟祥凱等[10]基于某緩傾斜薄—中厚鐵礦床礦體儲量大、品位低的特點,提出兩步中深孔聯(lián)合采礦法。諸利一等[11]提出了單層淺孔條帶式充填采礦法、多層合采條帶式充填采礦法和中深孔條帶充填采礦法等多層礦開采方法。對于層間距小的多層礦體,解聯(lián)庫等[12]利用綜合技術(shù)經(jīng)濟比較分析方法,分別研究了合采和分采過程中的合理采礦方法。此外,陳慶發(fā)[13-14]提出的協(xié)同開采理論也為多層礦開采提供了新思路。綜上可知,已有的多層礦圍巖穩(wěn)定性研究,均是基于傳統(tǒng)單層礦開采條件做出的拓展,鮮有以多層礦體間開采擾動導(dǎo)致的特殊地壓環(huán)境為基礎(chǔ)條件開展研究。

    煤炭資源賦存特征與含釩頁巖礦體較為近似,均呈厚度不等的多層狀排布,但煤炭資源賦存條件相對簡單,已經(jīng)形成了系統(tǒng)性極強的開采理論和方法。如采用綜采放頂煤方式回采時,將開采工作面附近的圍巖按應(yīng)力分布特征分為應(yīng)力穩(wěn)定區(qū)、應(yīng)力增高區(qū)、應(yīng)力恢復(fù)區(qū)和原始應(yīng)力區(qū),4 個區(qū)域隨著開采工作面向前推進,根據(jù)其不同的應(yīng)力特征可采取針對性較強的地壓管理措施。為此,本研究借鑒煤礦地壓分析理論,在多層礦體開采過程對采場圍巖進行分區(qū),并探究不同區(qū)域的應(yīng)力遷移特征,為多層狀金屬礦體開采理論研究與實踐提供參考。

    1 單層礦應(yīng)力遷移特征

    1.1 采空區(qū)圍巖應(yīng)力分布規(guī)律

    在研究開采過程中的采空區(qū)穩(wěn)定性問題時,常將復(fù)雜的地質(zhì)條件進行簡化,以便建立數(shù)學(xué)模型并分析其應(yīng)力狀態(tài)。對于一個處于某種應(yīng)力狀態(tài)下的無限彈性各向同性平面,用U0(x,y)表示對應(yīng)這種應(yīng)力狀態(tài)的應(yīng)力函數(shù),并將無孔平面內(nèi)的應(yīng)力狀態(tài)稱為基本應(yīng)力狀態(tài)[15-16]。

    假如在所研究的平面內(nèi)作任意形狀的孔(圖1),那么孔的存在就會引起平面內(nèi)應(yīng)力重新分布。新的應(yīng)力狀態(tài)可以表示為

    圖1 簡化的含孔平面Fig.1 Simplified plane with holes

    對于基本應(yīng)力狀態(tài)和附加應(yīng)力分量,根據(jù)雙調(diào)和函數(shù)可以分別找出兩個復(fù)變函數(shù)的穆斯海里什維里函數(shù)φ0(z),ψ0(z)和φ?(z),ψ?(z)[17-18]。與新的應(yīng)力狀態(tài)(式(1))相對應(yīng)的函數(shù)φ1(z)和ψ1(z)為

    式中,z為復(fù)數(shù);φ1(z)和ψ1(z)為變化后的應(yīng)力函數(shù);φ0(z),ψ0(z)為初始應(yīng)力狀態(tài)對應(yīng)的穆斯海里什維里函數(shù);φ?(z),ψ?(z)為附加應(yīng)力狀態(tài)對應(yīng)的穆斯海里什維里函數(shù)。

    當(dāng)用與OX軸成α角的力p拉伸連續(xù)平板時,板的基本應(yīng)力狀態(tài)可表示為

    式中,σx,σy,τxy為基本應(yīng)力分量,MPa;p為作用于平板上的力,MPa;α為力與OX軸夾角,(°)。

    任一處的應(yīng)力函數(shù)可表示為

    式中,U0(x,y)為平板上任一處應(yīng)力函數(shù)表達式;x、y為該處對應(yīng)的坐標(biāo)值。

    假如在已知平板中作一孔使其中心和XOY坐標(biāo)系的原點相合,并且其周界不受外力,那么帶孔平板的應(yīng)力函數(shù)φ1(z)和ψ1(z)可表示為

    式中,i為虛數(shù)單位。

    當(dāng)所作孔的形狀為矩形、正方形、三角形、梯形等有限部分的平面時,如果能將已知區(qū)域映射為一單位圓的函數(shù)z=ω(ζ),那么問題就可以用泛函方程組確定邊界條件φ(ζ)和ψ(ζ)來解決[19]。z=ω(ζ)見下式:

    根據(jù)模擬需求,僅對α=π/2 的情況進行如下分析計算:

    式中,φ(ζ),ψ(ζ)為泛函方程組確定的邊界條件;R為矩形孔頂點處的曲率半徑。

    將式(8)代入式(9),分開實部和虛部,得到應(yīng)力分量σρ,σθ和τρθ的表達式:

    將該情況的矩形孔周界上的σθ被p除得的值列出,結(jié)果見表1。

    表1 長寬比為5 的條件下α=π/2 時的σθ/p 值Table 1 σθ/p values at α=π/2 for an aspect ratio of 5

    對表1 中的數(shù)據(jù)進行擬合,所得曲線如圖2所示。

    圖2 不同角度下的應(yīng)力變化率擬合曲線Fig.2 Fitting curve of stress change rates at different angles

    1.2 單采場應(yīng)力變化區(qū)域劃分

    根據(jù)某含釩頁巖多層礦的巖石力學(xué)參數(shù)建立FLAC3D數(shù)值模型(表2),對上述采空區(qū)圍巖應(yīng)力分布規(guī)律進行驗證。該礦體埋深185 m,采場寬度為10 m,以1 m/次的掘進速度推進,分別得到5、10、15、20、25 m采空區(qū)頂板圍巖的垂直應(yīng)力分布曲線(拉應(yīng)力為正值、壓應(yīng)力為負值),如圖3所示,0 m 為初始開采位置。采空區(qū)形成后,頂板產(chǎn)生變形,應(yīng)力得到釋放,頂板圍巖的垂直應(yīng)力迅速跌落至穩(wěn)定值。開采工作面附近的圍巖承受采空區(qū)上覆巖體向周圍巖體轉(zhuǎn)移的荷載,此時,垂直應(yīng)力向遠離工作面的方向迅速增加至峰值,應(yīng)力峰值隨著開采寬度增加而增高,一段距離后應(yīng)力逐漸降低至原巖應(yīng)力。根據(jù)圖3所示的應(yīng)力遷移規(guī)律,可將采場圍巖劃分為如圖4所示的應(yīng)力釋放區(qū)、應(yīng)力跌落區(qū)、應(yīng)力增高區(qū)和原始應(yīng)力區(qū),垂直應(yīng)力的變化率分布與圖2 擬合曲線基本符合。

    表2 巖石力學(xué)參數(shù)Table 2 Rock mechanical parameter

    圖3 不同開采進度礦體上部圍巖垂直應(yīng)力分布曲線Fig.3 Vertical stress distribution curves of the upper surrounding rock of orebody with different mining progress

    圖4 單采場應(yīng)力變化分區(qū)Fig.4 Stress variation zones of single stope

    應(yīng)力釋放區(qū)和應(yīng)力增高區(qū)對單采場的穩(wěn)定性影響較大,目前主要的研究集中于這兩個分區(qū),即采空區(qū)的頂板、側(cè)壁和礦柱。應(yīng)力釋放區(qū)對應(yīng)采場的頂板位置,該區(qū)域圍巖下部臨空面受拉應(yīng)力作用,極易產(chǎn)生冒落甚至垮塌。應(yīng)力增高區(qū)對應(yīng)采場側(cè)壁或礦柱位置,該區(qū)域承受高于原巖應(yīng)力的荷載,提高了圍巖破壞的風(fēng)險,若此時能量積聚較高且瞬時釋放,極易產(chǎn)生巖爆。

    1.3 單采場應(yīng)力增高區(qū)應(yīng)力分布特征

    針對單采場,根據(jù)式(10)計算不同開采過程中的圍巖垂直應(yīng)力變化率α:

    式中,α為圍巖垂直應(yīng)力變化率;σi為距開采起點im處的當(dāng)前垂直應(yīng)力,MPa;σ0為初始垂直應(yīng)力,MPa。

    如圖5所示,將α>0.1 的區(qū)域劃分為高應(yīng)力增加區(qū),0.01<α<0.1 的區(qū)域劃分為低應(yīng)力增加區(qū),并在圖中標(biāo)記出應(yīng)力峰值區(qū)。根據(jù)圖5 可知,開采5 m 后出現(xiàn)高應(yīng)力增加區(qū),隨著開采的推進,高應(yīng)力增加區(qū)最終穩(wěn)定在距離開采工作面5 m 內(nèi),低應(yīng)力增加區(qū)范圍穩(wěn)定在距離開采工作面5~19 m,應(yīng)力增高區(qū)峰值出現(xiàn)在距離開采工作面2 m 處。

    不同跨度采空區(qū)對應(yīng)的應(yīng)力增高區(qū)應(yīng)力峰值不同,采空區(qū)跨度越大,峰值越高。利用上述結(jié)果擬合得到的峰值應(yīng)力變化率表達式為

    式中,Lc為采空區(qū)跨度,m;R2為擬合優(yōu)度。

    選取不同跨度的采空區(qū),分析了應(yīng)力增高區(qū)的應(yīng)力變化率隨距離的變化特征如圖6所示。由圖6 可知:當(dāng)采空區(qū)寬度為5~24 m 時,在應(yīng)力峰值區(qū)以外的區(qū)域,應(yīng)力變化率滿足式(12),其中a1=0.129~0.217,a2=-0.240~-0.201。

    圖6 應(yīng)力增高區(qū)應(yīng)力變化率隨距離變化曲線Fig.6 Variation curves of stress change rate with distance in the area of increased stress area

    式中,L為與應(yīng)力峰值點距離,m;a1、a2為常數(shù)。

    1.4 單采場應(yīng)力跌落區(qū)應(yīng)力分布特征

    應(yīng)力跌落區(qū)是應(yīng)力增高區(qū)與應(yīng)力釋放區(qū)的過渡區(qū)域,不同寬度采空區(qū)的應(yīng)力變化率隨著應(yīng)力增高區(qū)峰值點距離的變化曲線如圖7所示。

    圖7 應(yīng)力跌落區(qū)應(yīng)力變化率隨距離變化曲線Fig.7 Variation curves of stress change rate with distance in attenuated stress area

    圖7所示曲線可用如下三次函數(shù)式進行表示:

    式中,b0,b1,b2,b3分別為擬合曲線三次函數(shù)式的常數(shù)系數(shù)。

    當(dāng)采空區(qū)尺寸為24 m 時,應(yīng)力跌落區(qū)的應(yīng)力變化率隨峰值點距離的關(guān)系可表示為

    1.5 多采場應(yīng)力增高區(qū)應(yīng)力變化率疊加特征

    煤層工作面開采時,放頂工作使應(yīng)力釋放區(qū)頂板圍巖下落充填采空區(qū),最終實現(xiàn)該區(qū)域內(nèi)的應(yīng)力穩(wěn)定。因此,煤礦開采中將此區(qū)域稱為應(yīng)力穩(wěn)定區(qū),使用崩落法開采層狀金屬或非金屬礦與之類似。使用充填法開采時,充填體代替崩落的廢石起到支承采空區(qū)頂板的作用,同樣能達到類似效果。在圍巖穩(wěn)固的情況下使用空場法開采,采空區(qū)頂板能夠依靠自身強度承受一定的荷載,但隨開采跨度增加,頂板中部所受拉力逐漸增大,考慮到巖體本身軟弱結(jié)構(gòu)面的存在,應(yīng)力釋放區(qū)的寬度不宜過大。除了利用充填等措施外,工程中會優(yōu)化采場參數(shù),控制采空區(qū)跨度,從而將應(yīng)力釋放區(qū)限制在合理寬度,具體方法是間隔一定距離留設(shè)礦柱。通過上文分析可知,采空區(qū)的應(yīng)力增高區(qū)必定位于礦柱中,兩個相鄰采空區(qū)的應(yīng)力增高區(qū)若在礦柱內(nèi)重疊,則進一步增加了礦柱荷載,增大了礦柱失穩(wěn)風(fēng)險(圖8)。

    圖8 不同寬度礦柱應(yīng)力疊加特征示意Fig.8 Schematic of stress superposition characteristics of pillars with different widths

    如圖9所示,設(shè)兩個跨度分別為20 m 和10 m 的采空區(qū)之間留設(shè)10 m 寬礦柱,模擬單獨開采10 m、單獨開采20 m 和兩者全部開采3 種工況,分別得到如圖10所示的采空區(qū)圍巖應(yīng)力變化率曲線。將兩個采空區(qū)單獨作用值相加,發(fā)現(xiàn)其與兩個采空區(qū)共同作用的模擬值一致,可以認為礦柱內(nèi)部應(yīng)力增高區(qū)的應(yīng)力變化規(guī)律由兩側(cè)采空區(qū)圍巖應(yīng)力代數(shù)疊加得到。當(dāng)采空區(qū)數(shù)目增加時,礦柱內(nèi)部應(yīng)力增高區(qū)的應(yīng)力變化規(guī)律同樣滿足上述規(guī)律。

    圖9 20 m 采空區(qū)-10 m 礦柱-10 m 采空區(qū)示意Fig.9 Schematic of 20 m goaf -10 m pillar -10 m goaf

    圖10 應(yīng)力增高區(qū)應(yīng)力變化率模擬值隨距離變化曲線Fig.10 Curves of simulated stress change rate with distance in increased stress area

    選取圖9所示的20 m 采空區(qū)-10 m 礦柱-10 m采空區(qū)的模型,利用式(12)求得3 種工況下的應(yīng)力變化率。如圖11所示,兩個采空區(qū)共同作用的應(yīng)力計算值與模擬值的變化規(guī)律基本一致。在對礦柱進行理論分析時,通常將其視為受均布荷載,然而通過上述分析可知,礦柱受到的荷載大小和規(guī)律應(yīng)由兩側(cè)采空區(qū)共同決定。

    圖11 應(yīng)力增高區(qū)應(yīng)力變化率計算值隨距離變化曲線Fig.11 Curves of calculated stress change rate with distance in increased stress area

    2 雙層礦體應(yīng)力遷移規(guī)律

    2.1 水平位置關(guān)系對礦巖應(yīng)力遷移特征的影響

    研究雙層礦體開采應(yīng)力遷移規(guī)律時,需要考慮兩層礦體的位置關(guān)系,按礦體的重疊關(guān)系與長度,可分為上覆型完全重疊、下懸型完全重疊、部分重疊和不重疊4 種類型。其物理含義可表述為:上覆型完全重疊的雙層礦,上層礦體長度大于下層礦體,且下層礦體完全與上層重疊;下懸型完全重疊的雙層礦,下層礦體長度大于上層礦體,且上層礦體完全與下層重疊;部分重疊的雙層礦,僅有部分區(qū)域重疊;不重疊的雙層礦,完全不存在重疊區(qū)域。由于礦體間距過大,層間應(yīng)力擾動不明顯,因此為便于觀察,建立若干組不同重疊關(guān)系的近距離礦體數(shù)值模型用于分析其應(yīng)力擾動特征,礦體間距設(shè)置為8 m。

    如圖12(a)所示,建立上覆型完全重疊的雙層礦數(shù)值模型,自上而下分別留設(shè)50、20 m 采空區(qū),先回采上層礦體,再回采下層,分析夾層應(yīng)力分布規(guī)律。分別單獨開采上下兩層礦體時,滿足應(yīng)力釋放區(qū)、應(yīng)力跌落區(qū)、應(yīng)力增高區(qū)和原始應(yīng)力區(qū)4 個應(yīng)力變化區(qū)域劃分(圖12(b))。兩層礦體全部開采后,應(yīng)力擾動疊加。下層采空區(qū)的4 個應(yīng)力變化區(qū)完全位于上層采空區(qū)的應(yīng)力釋放區(qū)下方,此時上層采空區(qū)的應(yīng)力釋放區(qū)可視為下層采空區(qū)的次生原始應(yīng)力區(qū),在此基礎(chǔ)上形成次生應(yīng)力釋放區(qū)、跌落區(qū)和增高區(qū)。根據(jù)圖12(c)最小主應(yīng)力曲線可知,夾層最大壓應(yīng)力在上層采空區(qū)的應(yīng)力增高區(qū)達到最大值,下層采空區(qū)的應(yīng)力增高區(qū)與上層采空區(qū)的應(yīng)力釋放區(qū)的重疊區(qū)域也出現(xiàn)較小的峰值。由圖12(d)最大主應(yīng)力分布曲線可知,拉應(yīng)力在兩層采空區(qū)的應(yīng)力釋放區(qū)重疊處最大。

    圖12 上覆型完全重疊雙層礦體Fig.12 Overlying type completely overlapped double-layer orebody

    本研究構(gòu)建的下懸型完全重疊的雙層礦數(shù)值模型如圖13(a)所示。模型自上而下開采分別留設(shè)的20、50 m 采空區(qū)。夾層的垂直應(yīng)力擾動主要受下層采空區(qū)控制,最大壓應(yīng)力出現(xiàn)在下層采空區(qū)的應(yīng)力增高區(qū),上層采空區(qū)的應(yīng)力增高區(qū)與下層采空區(qū)的應(yīng)力釋放區(qū)重疊處的拉應(yīng)力最大。

    圖13 下懸型完全重疊雙層礦體Fig.13 Overhanging type completely overlapped double-layer orebody

    如圖14所示,對于部分重疊的雙層礦,自上而下開采,左側(cè)非重疊區(qū)的應(yīng)力擾動受下層采空區(qū)控制,右側(cè)非重疊區(qū)的應(yīng)力擾動受上層采空區(qū)控制,重疊區(qū)域則因兩層采空區(qū)的擾動疊加產(chǎn)生變化。最大壓應(yīng)力出現(xiàn)在下層采空區(qū)非重疊區(qū)域的應(yīng)力增高區(qū),最大拉應(yīng)力出現(xiàn)在下層采空區(qū)的應(yīng)力釋放區(qū)與上層采空區(qū)的應(yīng)力增高區(qū)重疊處。

    圖14 部分重疊雙層礦體Fig.14 Partially overlapped double-layer orebody

    當(dāng)由上至下依次開采形成兩層不重疊的采空區(qū)時,若距離足夠遠,則應(yīng)力擾動區(qū)域互不影響;若距離小于任一采空區(qū)的應(yīng)力增高區(qū)擾動范圍,則會出現(xiàn)應(yīng)力增高區(qū)疊加。如圖15所示,最大垂直應(yīng)力和最大壓應(yīng)力均出現(xiàn)在兩層采空區(qū)的應(yīng)力增高區(qū)重疊部位,最大拉應(yīng)力處于下層采空區(qū)的應(yīng)力釋放區(qū)中部。

    圖15 不重疊雙層礦體Fig.15 Non-overlapped double-layer orebody

    上述4 種重疊關(guān)系中,下層采空區(qū)頂板受到的拉應(yīng)力較大,尤其當(dāng)下層采空區(qū)的應(yīng)力釋放區(qū)與上層采空區(qū)的應(yīng)力增高區(qū)重疊時,拉應(yīng)力最大;最大壓應(yīng)力位于較大跨距采空區(qū)的應(yīng)力增高區(qū),當(dāng)上下兩層采空區(qū)的應(yīng)力增高區(qū)重疊時,重疊處的壓應(yīng)力進一步增加。

    2.2 垂直距離對礦巖應(yīng)力遷移規(guī)律的影響

    雙層礦體的相互擾動程度和范圍受礦體間距影響,礦體間距足夠大時,相鄰礦層不會形成開采擾動。因此,有必要研究礦層間距對開采擾動的影響,確定不同礦層間距的應(yīng)力擾動特征。在分析不同采空區(qū)長度的垂直應(yīng)力變化率后,設(shè)置采空區(qū)長度為24 m,模擬得到采空區(qū)下方1,3,6,11,16,21 m 垂直應(yīng)力變化率分布曲線(圖16)。上層礦體的開采擾動隨著深度增加而逐漸降低,應(yīng)力增高區(qū)在采空區(qū)以下6 m 位置幾乎完全消失。應(yīng)力跌落區(qū)的應(yīng)力跌落速率逐漸放緩。應(yīng)力釋放區(qū)的影響范圍較大,采空區(qū)下方21 m處應(yīng)力擾動依然存在,但已經(jīng)回落至較小值。因此,可按應(yīng)力增高區(qū)和應(yīng)力釋放區(qū)的擾動范圍,將應(yīng)力擾動區(qū)域劃分為強擾動區(qū)、弱擾動區(qū)和無擾動區(qū),3 個區(qū)域分界線為強擾動距和弱擾動距。礦體間距小于強擾動距時,應(yīng)力增高區(qū)和應(yīng)力釋放區(qū)共同對相鄰礦層造成擾動,形成強擾動區(qū);礦體間距大于強擾動距且小于弱擾動距時,應(yīng)力釋放區(qū)對相鄰礦層造成擾動,應(yīng)力增高區(qū)不造成擾動,形成弱擾動區(qū);礦體間距大于弱擾動距時,應(yīng)力增高區(qū)和釋放區(qū)對相鄰礦體均不造成擾動,此處成為無擾動區(qū)。綜合上述數(shù)值模擬結(jié)果得到該模型中強擾動距為6 m,弱擾動距為22 m。

    圖16 單層采空區(qū)對下部巖體擾動分布Fig.16 Distribution of disturbance caused by single goaf to lower rock mass

    以不完全重疊的雙層礦體為例,根據(jù)24 m 采空區(qū)的應(yīng)力擾動區(qū)分布情況設(shè)置3 組數(shù)值模型,礦層間距為3,15,25 m。如圖17所示,礦層間距為3 m 時,下層采空區(qū)位于上層采空區(qū)的強擾動區(qū),受上層采空區(qū)應(yīng)力增高區(qū)和應(yīng)力釋放區(qū)共同影響。如圖17(c)和圖17(d)所示,上層采空區(qū)的應(yīng)力釋放區(qū)導(dǎo)致下層采空區(qū)右側(cè)壓應(yīng)力和拉應(yīng)力降低,上層采空區(qū)的應(yīng)力增高區(qū)導(dǎo)致下層采空區(qū)右側(cè)壁出現(xiàn)小型壓應(yīng)力峰值,中部拉應(yīng)力有較大幅度增加。

    如圖18所示,礦層間距為15 m 時,下層采空區(qū)位于上層采空區(qū)的弱擾動區(qū),受上層采空區(qū)應(yīng)力釋放區(qū)影響,表現(xiàn)為下層采空區(qū)右側(cè)壓應(yīng)力小幅度降低,拉應(yīng)力小幅度增高,但總體變化值較小(圖18(c)和圖18(d))。

    圖18 夾層為15 m 的應(yīng)力擾動情況Fig.18 Stress disturbance of a 15 m interlayer

    如圖19所示,礦層間距為25 m 時,下層采空區(qū)位于上層采空區(qū)的無擾動區(qū),單獨開采下層礦體與開采雙層礦體時的應(yīng)力變化幾乎相同,可認為上層采空區(qū)未對下層礦體開采造成影響。

    圖19 夾層為25 m 的應(yīng)力擾動情況Fig.19 Stress disturbance of a 25 m interlayer

    3 結(jié)論

    (1)根據(jù)含釩礦體開采后的垂直應(yīng)力變化規(guī)律,將采空區(qū)圍巖劃分為應(yīng)力釋放區(qū)、應(yīng)力跌落區(qū)、應(yīng)力增高區(qū)和原始應(yīng)力區(qū)等4 個區(qū)域。

    (2)單采場采空區(qū)跨度與其圍巖應(yīng)力峰值為正相關(guān)關(guān)系,即采空區(qū)跨度越大,其應(yīng)力峰值越高,并且應(yīng)力增高區(qū)以深圍巖內(nèi)的應(yīng)力變化率隨著與應(yīng)力峰值點距離的增大呈冪函數(shù)關(guān)系減小,應(yīng)力跌落區(qū)內(nèi)的應(yīng)力變化率隨著與應(yīng)力峰值點距離的增大而增大;相鄰采空區(qū)的應(yīng)力增高區(qū)會在兩者間的礦柱內(nèi)疊加,礦柱內(nèi)某點的應(yīng)力變化率為兩側(cè)采空區(qū)分別在該點引起的應(yīng)力變化率之和。

    (3)雙層礦雙采場開采后,當(dāng)下層采空區(qū)應(yīng)力釋放區(qū)與上層采空區(qū)應(yīng)力增高區(qū)重疊時導(dǎo)致圍巖內(nèi)拉應(yīng)力增高,當(dāng)上層采空區(qū)的應(yīng)力釋放區(qū)與下層采空區(qū)任一區(qū)域重疊時導(dǎo)致圍巖壓應(yīng)力降低,并且上下采空區(qū)重疊與非重疊區(qū)域的交匯處應(yīng)力變化顯著。

    (4)應(yīng)進一步分析多層礦體多采場采空區(qū)圍巖極限平衡狀態(tài),并根據(jù)其圍巖應(yīng)力遷移規(guī)律和失穩(wěn)特征合理確定回采參數(shù)與回采順序,以實現(xiàn)安全、高效、經(jīng)濟開采。

    猜你喜歡
    礦柱礦層下層
    某銅鈷礦柱浸試驗研究
    云南化工(2021年9期)2021-12-21 07:43:40
    新疆和田市沉積礦床礦體特征與巖層特征對比研究
    河南舞陽鐵礦田鐵山廟組礦層下盤地質(zhì)體成因
    開采殘余礦層存在的風(fēng)險及防治對策
    一類多個下層的雙層規(guī)劃問題
    積雪
    基于動態(tài)強度折減法礦柱漸進失穩(wěn)演變規(guī)律
    陜西橫山羅圪臺村元代壁畫墓發(fā)掘簡報
    考古與文物(2016年5期)2016-12-21 06:28:48
    有借有還
    貴州省岑鞏縣注溪釩礦含礦層特征與找礦前景
    河南科技(2014年3期)2014-02-27 14:06:02
    日本黄色片子视频| 99热全是精品| 综合色av麻豆| 免费看光身美女| 国产伦在线观看视频一区| 国产伦理片在线播放av一区| 在线观看人妻少妇| 亚洲精品中文字幕在线视频 | 熟女电影av网| 久久久a久久爽久久v久久| 嘟嘟电影网在线观看| 美女cb高潮喷水在线观看| 日本爱情动作片www.在线观看| 久久久精品94久久精品| 七月丁香在线播放| 99久久九九国产精品国产免费| 久久久欧美国产精品| 中文字幕亚洲精品专区| 成人av在线播放网站| 午夜福利视频1000在线观看| 午夜福利在线观看吧| 亚洲精品影视一区二区三区av| av在线播放精品| 国产精品三级大全| 99热网站在线观看| 又大又黄又爽视频免费| 少妇被粗大猛烈的视频| 亚洲成人中文字幕在线播放| 欧美精品一区二区大全| videos熟女内射| 国产精品一区www在线观看| 亚洲最大成人手机在线| 久久韩国三级中文字幕| av国产久精品久网站免费入址| 午夜福利视频精品| 欧美zozozo另类| 国产精品爽爽va在线观看网站| 成人特级av手机在线观看| 亚洲av福利一区| 日日撸夜夜添| 国产av不卡久久| 国产精品爽爽va在线观看网站| 99九九线精品视频在线观看视频| 日本av手机在线免费观看| 免费看美女性在线毛片视频| 亚洲国产欧美人成| 久久久国产一区二区| 国产精品福利在线免费观看| 国产成人精品久久久久久| 欧美三级亚洲精品| 在线观看av片永久免费下载| 精品熟女少妇av免费看| 欧美日韩综合久久久久久| 人体艺术视频欧美日本| 日本-黄色视频高清免费观看| 国产精品久久久久久久电影| 日本-黄色视频高清免费观看| 人体艺术视频欧美日本| 免费黄频网站在线观看国产| 一本一本综合久久| 国产淫片久久久久久久久| 激情五月婷婷亚洲| 精品人妻视频免费看| 黄片无遮挡物在线观看| 汤姆久久久久久久影院中文字幕 | 国产成人精品婷婷| 国产黄a三级三级三级人| 床上黄色一级片| 男人舔女人下体高潮全视频| 成人国产麻豆网| 精华霜和精华液先用哪个| 婷婷色综合大香蕉| 免费观看a级毛片全部| 国内揄拍国产精品人妻在线| 爱豆传媒免费全集在线观看| 色吧在线观看| 99久久精品热视频| 国产成人91sexporn| 99久久中文字幕三级久久日本| 国产成人精品婷婷| 亚洲最大成人av| 偷拍熟女少妇极品色| 久久99热这里只频精品6学生| 亚洲精品久久午夜乱码| 国产高清有码在线观看视频| 国产成人a∨麻豆精品| 插阴视频在线观看视频| av女优亚洲男人天堂| 麻豆乱淫一区二区| 女的被弄到高潮叫床怎么办| 欧美极品一区二区三区四区| 99久久人妻综合| 亚洲在线自拍视频| 国产伦精品一区二区三区四那| 欧美97在线视频| 国产成人一区二区在线| 亚洲av二区三区四区| 肉色欧美久久久久久久蜜桃 | 亚洲精品视频女| 美女主播在线视频| 国产男人的电影天堂91| 亚洲内射少妇av| 免费观看性生交大片5| 成人av在线播放网站| 成人鲁丝片一二三区免费| 欧美成人精品欧美一级黄| 免费黄网站久久成人精品| 老女人水多毛片| 人人妻人人看人人澡| 国产毛片a区久久久久| 我的老师免费观看完整版| 91aial.com中文字幕在线观看| 精品一区二区免费观看| 亚洲18禁久久av| 免费大片18禁| 91精品伊人久久大香线蕉| 国模一区二区三区四区视频| 男人和女人高潮做爰伦理| 国产成人aa在线观看| 欧美 日韩 精品 国产| 精品亚洲乱码少妇综合久久| 国产色婷婷99| 久久久久久久久中文| 特级一级黄色大片| 欧美日韩综合久久久久久| 国产爱豆传媒在线观看| av在线观看视频网站免费| 欧美bdsm另类| 99久久中文字幕三级久久日本| 男人狂女人下面高潮的视频| 激情 狠狠 欧美| 欧美xxxx性猛交bbbb| 久久久亚洲精品成人影院| 日韩av不卡免费在线播放| 色视频www国产| 如何舔出高潮| 国产大屁股一区二区在线视频| 久久久久久久亚洲中文字幕| 不卡视频在线观看欧美| 永久网站在线| 亚洲欧美日韩无卡精品| 国产精品av视频在线免费观看| 亚洲精品日韩在线中文字幕| 少妇丰满av| av又黄又爽大尺度在线免费看| 一级毛片久久久久久久久女| 能在线免费看毛片的网站| 非洲黑人性xxxx精品又粗又长| 中文字幕免费在线视频6| 国产成人精品一,二区| 亚洲在久久综合| 免费观看av网站的网址| 国产精品福利在线免费观看| 日本午夜av视频| 国产精品一区二区三区四区久久| 在线观看美女被高潮喷水网站| 好男人视频免费观看在线| 老司机影院成人| 超碰97精品在线观看| 午夜精品一区二区三区免费看| 熟妇人妻久久中文字幕3abv| 成年版毛片免费区| 亚洲精品自拍成人| 久久精品夜夜夜夜夜久久蜜豆| 欧美精品国产亚洲| 欧美日韩视频高清一区二区三区二| 国产精品久久视频播放| 久久99热6这里只有精品| 亚洲人成网站高清观看| 日韩国内少妇激情av| 国产高潮美女av| 少妇被粗大猛烈的视频| 一级a做视频免费观看| 免费大片18禁| 午夜爱爱视频在线播放| 97精品久久久久久久久久精品| 国产v大片淫在线免费观看| 久久久a久久爽久久v久久| 黄色配什么色好看| 国内精品美女久久久久久| 最新中文字幕久久久久| 97在线视频观看| 一本久久精品| 韩国高清视频一区二区三区| 夜夜爽夜夜爽视频| 国产成人freesex在线| 国产一级毛片七仙女欲春2| 91aial.com中文字幕在线观看| 免费电影在线观看免费观看| 欧美一区二区亚洲| 免费观看无遮挡的男女| 韩国av在线不卡| 国产高清三级在线| 十八禁网站网址无遮挡 | 日韩一本色道免费dvd| 亚洲av.av天堂| 激情 狠狠 欧美| 日韩亚洲欧美综合| 女的被弄到高潮叫床怎么办| 久久精品熟女亚洲av麻豆精品 | 国产午夜精品久久久久久一区二区三区| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 麻豆成人av视频| 国产伦一二天堂av在线观看| 国产精品久久久久久久电影| 亚洲精品日韩av片在线观看| 久久久精品94久久精品| 18禁裸乳无遮挡免费网站照片| 成人二区视频| 亚洲欧美成人精品一区二区| 麻豆久久精品国产亚洲av| 日韩伦理黄色片| 观看美女的网站| 99热6这里只有精品| 国产免费又黄又爽又色| 久久久久久久久久人人人人人人| 汤姆久久久久久久影院中文字幕 | 欧美日韩国产mv在线观看视频 | 日韩伦理黄色片| 观看美女的网站| 能在线免费观看的黄片| 久久久久久久久久久丰满| 99re6热这里在线精品视频| 男女啪啪激烈高潮av片| 亚州av有码| 精品久久久久久久末码| 欧美成人精品欧美一级黄| 亚洲丝袜综合中文字幕| 免费看不卡的av| 国产欧美日韩精品一区二区| 你懂的网址亚洲精品在线观看| 九九久久精品国产亚洲av麻豆| 日本免费a在线| 简卡轻食公司| 青春草视频在线免费观看| 寂寞人妻少妇视频99o| 国产欧美另类精品又又久久亚洲欧美| 欧美日韩国产mv在线观看视频 | 街头女战士在线观看网站| 亚洲av电影在线观看一区二区三区 | 老司机影院毛片| 亚洲精品,欧美精品| 色哟哟·www| 精品人妻熟女av久视频| 啦啦啦韩国在线观看视频| 你懂的网址亚洲精品在线观看| 久久99蜜桃精品久久| 我要看日韩黄色一级片| 99热6这里只有精品| 小蜜桃在线观看免费完整版高清| 身体一侧抽搐| xxx大片免费视频| 国产老妇伦熟女老妇高清| 国产综合懂色| 大陆偷拍与自拍| 国产成人91sexporn| 亚洲高清免费不卡视频| 蜜桃亚洲精品一区二区三区| 91狼人影院| 日韩人妻高清精品专区| 有码 亚洲区| 69av精品久久久久久| 国产精品伦人一区二区| 少妇的逼水好多| 国产极品天堂在线| 我要看日韩黄色一级片| 三级国产精品片| 欧美高清成人免费视频www| 日日啪夜夜撸| 深爱激情五月婷婷| 日本黄色片子视频| 久久久久免费精品人妻一区二区| 亚洲精品视频女| 国产男人的电影天堂91| 男人爽女人下面视频在线观看| 亚洲丝袜综合中文字幕| 国产免费福利视频在线观看| 国产黄a三级三级三级人| 在线免费观看不下载黄p国产| 亚洲丝袜综合中文字幕| 寂寞人妻少妇视频99o| 激情五月婷婷亚洲| 在线 av 中文字幕| 黄色欧美视频在线观看| 国产精品一区二区性色av| 国产一区亚洲一区在线观看| 五月天丁香电影| 大又大粗又爽又黄少妇毛片口| 精品人妻偷拍中文字幕| 日韩强制内射视频| 99热6这里只有精品| 精品久久久久久久人妻蜜臀av| 亚洲精品日韩在线中文字幕| 综合色av麻豆| 亚洲精品成人av观看孕妇| 午夜免费观看性视频| 午夜福利高清视频| 夫妻午夜视频| 国产麻豆成人av免费视频| 成人午夜精彩视频在线观看| 天天躁夜夜躁狠狠久久av| 一边亲一边摸免费视频| 亚洲四区av| 亚洲av中文字字幕乱码综合| 久久久亚洲精品成人影院| 欧美另类一区| 久久精品国产亚洲网站| 亚洲性久久影院| 免费观看a级毛片全部| 搡老妇女老女人老熟妇| 久久久久久久大尺度免费视频| 亚洲最大成人中文| 日本与韩国留学比较| 激情五月婷婷亚洲| 能在线免费观看的黄片| 夜夜骑夜夜射夜夜干| 欧美另类一区| 久久精品夜色国产| 国产xxxxx性猛交| 伦精品一区二区三区| 各种免费的搞黄视频| 久久久久久人人人人人| 欧美激情高清一区二区三区 | 午夜激情久久久久久久| 亚洲国产精品国产精品| 精品国产乱码久久久久久小说| 国产日韩欧美视频二区| 国产免费福利视频在线观看| 我要看黄色一级片免费的| 国产精品不卡视频一区二区| 伦理电影免费视频| 午夜福利一区二区在线看| 亚洲av男天堂| 丝瓜视频免费看黄片| 亚洲天堂av无毛| 麻豆乱淫一区二区| 午夜日本视频在线| 免费黄色在线免费观看| 亚洲精品aⅴ在线观看| 日韩欧美精品免费久久| 色网站视频免费| 中文字幕av电影在线播放| 国产成人午夜福利电影在线观看| 激情五月婷婷亚洲| 赤兔流量卡办理| av网站在线播放免费| 免费黄频网站在线观看国产| 中国三级夫妇交换| 欧美av亚洲av综合av国产av | av卡一久久| 亚洲精品美女久久久久99蜜臀 | 日韩中文字幕视频在线看片| 一个人免费看片子| 最近最新中文字幕免费大全7| 亚洲美女搞黄在线观看| 日本wwww免费看| 亚洲国产欧美日韩在线播放| 日本色播在线视频| 久久久久久人人人人人| 亚洲天堂av无毛| 婷婷色av中文字幕| 波多野结衣一区麻豆| 建设人人有责人人尽责人人享有的| 天天躁夜夜躁狠狠久久av| 久久久久网色| 成人毛片a级毛片在线播放| av不卡在线播放| 老司机影院毛片| 久久 成人 亚洲| 老汉色av国产亚洲站长工具| 在线观看美女被高潮喷水网站| 午夜免费观看性视频| 99香蕉大伊视频| 99久久精品国产国产毛片| 国产精品.久久久| 中文字幕最新亚洲高清| 涩涩av久久男人的天堂| 99热国产这里只有精品6| 欧美成人精品欧美一级黄| 久久精品亚洲av国产电影网| 激情视频va一区二区三区| 2022亚洲国产成人精品| 国产精品一二三区在线看| 亚洲av中文av极速乱| 狂野欧美激情性bbbbbb| 亚洲精品日韩在线中文字幕| 久久精品久久久久久久性| 免费人妻精品一区二区三区视频| 韩国高清视频一区二区三区| 中文欧美无线码| 麻豆精品久久久久久蜜桃| 18在线观看网站| 精品少妇一区二区三区视频日本电影 | 亚洲国产精品成人久久小说| 男女下面插进去视频免费观看| 9191精品国产免费久久| 国产熟女欧美一区二区| 欧美 亚洲 国产 日韩一| 精品国产乱码久久久久久男人| 丰满少妇做爰视频| 亚洲第一区二区三区不卡| 欧美人与性动交α欧美软件| 在线看a的网站| 中文乱码字字幕精品一区二区三区| 国产成人精品久久久久久| 亚洲四区av| 寂寞人妻少妇视频99o| 欧美国产精品va在线观看不卡| 免费女性裸体啪啪无遮挡网站| 亚洲av免费高清在线观看| 久久精品国产亚洲av高清一级| 2021少妇久久久久久久久久久| 狠狠婷婷综合久久久久久88av| 婷婷色综合www| 制服丝袜香蕉在线| 国产 精品1| 日韩av在线免费看完整版不卡| 久久久久国产精品人妻一区二区| 中文欧美无线码| 飞空精品影院首页| 少妇人妻久久综合中文| 少妇精品久久久久久久| 大片电影免费在线观看免费| 少妇被粗大的猛进出69影院| 国产日韩一区二区三区精品不卡| 日本黄色日本黄色录像| 国产成人一区二区在线| 亚洲,一卡二卡三卡| 老汉色av国产亚洲站长工具| 午夜日韩欧美国产| 日本-黄色视频高清免费观看| 老司机影院成人| 美女大奶头黄色视频| 一二三四在线观看免费中文在| 午夜日本视频在线| 飞空精品影院首页| 桃花免费在线播放| 久久狼人影院| 中国国产av一级| 深夜精品福利| 卡戴珊不雅视频在线播放| 亚洲经典国产精华液单| av不卡在线播放| 国产精品 国内视频| 精品少妇一区二区三区视频日本电影 | 亚洲在久久综合| 久久久久久久精品精品| 免费高清在线观看日韩| 亚洲精品日本国产第一区| 亚洲激情五月婷婷啪啪| 蜜桃国产av成人99| 精品一区二区免费观看| 国产又爽黄色视频| 国产男女内射视频| 秋霞伦理黄片| 久久综合国产亚洲精品| 亚洲欧洲日产国产| 一二三四在线观看免费中文在| 成年人免费黄色播放视频| 久久久精品国产亚洲av高清涩受| 免费在线观看视频国产中文字幕亚洲 | 成人18禁高潮啪啪吃奶动态图| 欧美亚洲日本最大视频资源| 亚洲美女视频黄频| 午夜福利视频在线观看免费| 亚洲欧洲日产国产| 国产在线视频一区二区| 国产在视频线精品| 日韩一本色道免费dvd| 丰满饥渴人妻一区二区三| 26uuu在线亚洲综合色| 精品人妻一区二区三区麻豆| 大片电影免费在线观看免费| 亚洲国产精品999| 啦啦啦在线免费观看视频4| 午夜福利在线免费观看网站| 老司机亚洲免费影院| 91午夜精品亚洲一区二区三区| 性色av一级| 国产国语露脸激情在线看| 国产亚洲av片在线观看秒播厂| 有码 亚洲区| 久久av网站| 美女福利国产在线| 大话2 男鬼变身卡| 性色av一级| 亚洲国产av影院在线观看| 久久国产精品男人的天堂亚洲| 精品一区在线观看国产| 岛国毛片在线播放| 免费在线观看视频国产中文字幕亚洲 | 丰满迷人的少妇在线观看| 在线天堂中文资源库| 一级毛片我不卡| 永久免费av网站大全| 亚洲精品在线美女| 一级毛片 在线播放| 久热这里只有精品99| 天堂8中文在线网| 午夜福利在线免费观看网站| 亚洲一区中文字幕在线| 香蕉国产在线看| 国产日韩欧美亚洲二区| √禁漫天堂资源中文www| 91精品三级在线观看| 人妻人人澡人人爽人人| 视频在线观看一区二区三区| 18禁观看日本| 在线观看三级黄色| 亚洲在久久综合| 夜夜骑夜夜射夜夜干| 18在线观看网站| 日韩制服丝袜自拍偷拍| 亚洲内射少妇av| 中文字幕人妻丝袜一区二区 | 91成人精品电影| 多毛熟女@视频| 看免费成人av毛片| 亚洲一区中文字幕在线| 美女主播在线视频| 亚洲三级黄色毛片| 国产成人av激情在线播放| 欧美日韩综合久久久久久| 黄色视频在线播放观看不卡| 欧美日韩精品成人综合77777| 丝袜喷水一区| 成人漫画全彩无遮挡| 国产精品一区二区在线观看99| 国产老妇伦熟女老妇高清| 日本av免费视频播放| 久久久a久久爽久久v久久| 国产人伦9x9x在线观看 | 丝袜脚勾引网站| av福利片在线| 秋霞伦理黄片| 嫩草影院入口| 国产极品天堂在线| 97在线人人人人妻| 国产精品女同一区二区软件| 精品国产露脸久久av麻豆| 日韩中文字幕欧美一区二区 | 这个男人来自地球电影免费观看 | 精品亚洲成国产av| 日韩中文字幕视频在线看片| 日韩视频在线欧美| 午夜激情久久久久久久| 国产淫语在线视频| 精品一区在线观看国产| 人妻一区二区av| 啦啦啦中文免费视频观看日本| 性高湖久久久久久久久免费观看| 热re99久久国产66热| 最近中文字幕2019免费版| 2021少妇久久久久久久久久久| 久久99一区二区三区| 在现免费观看毛片| 伊人亚洲综合成人网| 欧美日韩精品网址| 国产黄色免费在线视频| 蜜桃国产av成人99| 亚洲美女视频黄频| 亚洲国产毛片av蜜桃av| 欧美精品人与动牲交sv欧美| 日日摸夜夜添夜夜爱| 美女大奶头黄色视频| 亚洲图色成人| 制服人妻中文乱码| 丝袜美腿诱惑在线| 热re99久久精品国产66热6| 色吧在线观看| 久久这里有精品视频免费| 中文字幕最新亚洲高清| 久热这里只有精品99| 亚洲国产精品一区二区三区在线| 国产精品久久久久成人av| 美女大奶头黄色视频| 欧美成人精品欧美一级黄| 亚洲精品乱久久久久久| 国产精品不卡视频一区二区| 国产精品久久久久久av不卡| 免费高清在线观看日韩| 亚洲欧美精品综合一区二区三区 | 日本猛色少妇xxxxx猛交久久| 久久亚洲国产成人精品v| 黄色毛片三级朝国网站| 国产精品av久久久久免费| 中文乱码字字幕精品一区二区三区| 久久精品亚洲av国产电影网| 国产日韩欧美亚洲二区| 久久ye,这里只有精品| 女性生殖器流出的白浆| 人妻少妇偷人精品九色| 欧美精品一区二区免费开放| 三上悠亚av全集在线观看| 99热网站在线观看| 国产成人午夜福利电影在线观看| 美女福利国产在线| 成人毛片a级毛片在线播放| 亚洲成人一二三区av| 桃花免费在线播放| 国产又色又爽无遮挡免| 久久97久久精品| 永久网站在线| 男女下面插进去视频免费观看| 美国免费a级毛片| 中文天堂在线官网| 飞空精品影院首页| 免费不卡的大黄色大毛片视频在线观看| 精品亚洲成a人片在线观看| 精品午夜福利在线看| 国产成人午夜福利电影在线观看| 午夜免费鲁丝| 国产男女超爽视频在线观看| 亚洲国产精品成人久久小说| 青草久久国产|