• <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
    999久久久国产精品视频| 亚洲精品中文字幕在线视频| 免费在线观看视频国产中文字幕亚洲| 大片电影免费在线观看免费| 欧美人与性动交α欧美精品济南到| 亚洲av日韩在线播放| 夜夜爽天天搞| 啪啪无遮挡十八禁网站| 嫁个100分男人电影在线观看| 亚洲一区二区三区欧美精品| 国产一区二区三区在线臀色熟女 | 亚洲欧美精品综合一区二区三区| 一级黄色大片毛片| 国产无遮挡羞羞视频在线观看| 亚洲专区中文字幕在线| 老熟女久久久| www.自偷自拍.com| 天天躁狠狠躁夜夜躁狠狠躁| 人妻 亚洲 视频| 国产欧美日韩精品亚洲av| 国产成人精品在线电影| 欧美精品人与动牲交sv欧美| 国产熟女午夜一区二区三区| 免费在线观看亚洲国产| 色播在线永久视频| 日韩免费av在线播放| 中文字幕人妻丝袜一区二区| 久99久视频精品免费| 亚洲欧美日韩高清在线视频| 亚洲国产欧美日韩在线播放| 90打野战视频偷拍视频| 欧美丝袜亚洲另类 | 熟女少妇亚洲综合色aaa.| 精品久久蜜臀av无| 一区二区三区精品91| 国产av又大| 两性夫妻黄色片| 后天国语完整版免费观看| 日本wwww免费看| 国产精品 欧美亚洲| 国产在线精品亚洲第一网站| 美女 人体艺术 gogo| 欧美在线一区亚洲| 99在线人妻在线中文字幕 | 这个男人来自地球电影免费观看| 在线天堂中文资源库| 色在线成人网| 在线观看免费视频日本深夜| 一区二区三区国产精品乱码| 免费久久久久久久精品成人欧美视频| 国产色视频综合| 亚洲性夜色夜夜综合| 精品一品国产午夜福利视频| 妹子高潮喷水视频| 亚洲国产欧美网| 久久久久国产精品人妻aⅴ院 | 精品一区二区三卡| 久久久久久免费高清国产稀缺| 亚洲av片天天在线观看| 国产成人啪精品午夜网站| 在线观看日韩欧美| 亚洲精品中文字幕一二三四区| 深夜精品福利| 精品久久久久久,| 999久久久精品免费观看国产| 妹子高潮喷水视频| 超碰成人久久| 亚洲在线自拍视频| 中文字幕色久视频| 久久午夜亚洲精品久久| 亚洲五月天丁香| 亚洲人成伊人成综合网2020| 一级片免费观看大全| 黄色视频不卡| 久久亚洲精品不卡| 久久精品国产a三级三级三级| 久热爱精品视频在线9| 精品高清国产在线一区| 欧美精品av麻豆av| 精品久久久久久久毛片微露脸| 国产真人三级小视频在线观看| 久久久久久免费高清国产稀缺| 两性午夜刺激爽爽歪歪视频在线观看 | 日韩中文字幕欧美一区二区| 超碰97精品在线观看| 高清欧美精品videossex| 老司机午夜十八禁免费视频| 成人免费观看视频高清| 亚洲美女黄片视频| 免费久久久久久久精品成人欧美视频| 欧美日韩福利视频一区二区| 国产精品久久久久久人妻精品电影| 国产欧美亚洲国产| av视频免费观看在线观看| 黑人巨大精品欧美一区二区蜜桃| 亚洲精品一卡2卡三卡4卡5卡| 亚洲色图综合在线观看| 美女扒开内裤让男人捅视频| 人妻丰满熟妇av一区二区三区 | 我的亚洲天堂| 欧美成人免费av一区二区三区 | xxx96com| 99久久综合精品五月天人人| 国产极品粉嫩免费观看在线| 色精品久久人妻99蜜桃| 国产精品自产拍在线观看55亚洲 | 国产激情久久老熟女| 色婷婷av一区二区三区视频| 欧美精品一区二区免费开放| 日本a在线网址| 亚洲人成77777在线视频| 欧美老熟妇乱子伦牲交| 国产精品二区激情视频| 99re在线观看精品视频| 成人国产一区最新在线观看| 成人特级黄色片久久久久久久| 欧美乱色亚洲激情| 国产成人免费无遮挡视频| 亚洲视频免费观看视频| 亚洲久久久国产精品| 精品视频人人做人人爽| av中文乱码字幕在线| 亚洲国产欧美日韩在线播放| 国产精品一区二区精品视频观看| 亚洲自偷自拍图片 自拍| 一区二区三区国产精品乱码| 我的亚洲天堂| 国产片内射在线| 免费久久久久久久精品成人欧美视频| 两个人看的免费小视频| 老司机在亚洲福利影院| 在线观看免费视频日本深夜| 又黄又爽又免费观看的视频| 丝瓜视频免费看黄片| 久久国产精品影院| 免费日韩欧美在线观看| 精品卡一卡二卡四卡免费| 黄色怎么调成土黄色| 男女高潮啪啪啪动态图| 大片电影免费在线观看免费| 亚洲精品乱久久久久久| 日韩一卡2卡3卡4卡2021年| 国产精品美女特级片免费视频播放器 | 99国产精品99久久久久| 久久精品国产清高在天天线| 精品人妻1区二区| 超碰97精品在线观看| 美女国产高潮福利片在线看| 日日摸夜夜添夜夜添小说| svipshipincom国产片| 高清欧美精品videossex| 国产精品久久久久久精品古装| 精品国产一区二区三区四区第35| 999久久久精品免费观看国产| 免费日韩欧美在线观看| 久久国产亚洲av麻豆专区| 高清在线国产一区| 最新在线观看一区二区三区| 桃红色精品国产亚洲av| 69av精品久久久久久| 国产麻豆69| 一进一出抽搐动态| 国产精品秋霞免费鲁丝片| 免费在线观看亚洲国产| 18禁黄网站禁片午夜丰满| 午夜精品久久久久久毛片777| 91国产中文字幕| 一区二区三区精品91| 国产成人影院久久av| 亚洲一区二区三区欧美精品| 久久天躁狠狠躁夜夜2o2o| 满18在线观看网站| 性色av乱码一区二区三区2| 99精国产麻豆久久婷婷| 99国产综合亚洲精品| 精品国产一区二区三区久久久樱花| 亚洲人成77777在线视频| 亚洲人成77777在线视频| 日韩视频一区二区在线观看| 一个人免费在线观看的高清视频| 男人的好看免费观看在线视频 | 岛国在线观看网站| 一边摸一边抽搐一进一出视频| 欧美国产精品va在线观看不卡| 日韩中文字幕欧美一区二区| 亚洲精品在线美女| 美女视频免费永久观看网站| 久久久精品区二区三区| 婷婷成人精品国产| 三上悠亚av全集在线观看| 午夜精品国产一区二区电影| 日韩制服丝袜自拍偷拍| 嫩草影视91久久| 男女午夜视频在线观看| 国产成人系列免费观看| 亚洲片人在线观看| www.999成人在线观看| 99re6热这里在线精品视频| 捣出白浆h1v1| 国产亚洲欧美在线一区二区| 国产一区二区三区在线臀色熟女 | 视频区欧美日本亚洲| 亚洲av电影在线进入| 国内毛片毛片毛片毛片毛片| 狠狠狠狠99中文字幕| 日韩人妻精品一区2区三区| 一边摸一边抽搐一进一小说 | 手机成人av网站| 女同久久另类99精品国产91| 国产精品久久久久久精品古装| 精品电影一区二区在线| 国产男靠女视频免费网站| 国产免费现黄频在线看| 精品免费久久久久久久清纯 | 午夜福利影视在线免费观看| 欧美性长视频在线观看| av免费在线观看网站| 精品国产超薄肉色丝袜足j| 三上悠亚av全集在线观看| 婷婷精品国产亚洲av在线 | 日韩一卡2卡3卡4卡2021年| 可以免费在线观看a视频的电影网站| 熟女少妇亚洲综合色aaa.| 成人18禁高潮啪啪吃奶动态图| 久久国产精品人妻蜜桃| 久久人人爽av亚洲精品天堂| 啦啦啦在线免费观看视频4| 久久精品国产综合久久久| 国产免费男女视频| 在线观看一区二区三区激情| 国产不卡av网站在线观看| 啦啦啦 在线观看视频| 国产精品av久久久久免费| 欧美在线一区亚洲| 香蕉国产在线看| 日本wwww免费看| 欧美另类亚洲清纯唯美| 久久人人爽av亚洲精品天堂| 天天躁日日躁夜夜躁夜夜| 97人妻天天添夜夜摸| 亚洲精品在线美女| 嫩草影视91久久| 亚洲国产看品久久| 99精品欧美一区二区三区四区| 午夜福利免费观看在线| 高清黄色对白视频在线免费看| 久久精品亚洲av国产电影网| 亚洲人成伊人成综合网2020| 久久国产精品人妻蜜桃| 少妇 在线观看| ponron亚洲| 午夜福利欧美成人| 人人妻人人爽人人添夜夜欢视频| 一个人免费在线观看的高清视频| 免费人成视频x8x8入口观看| 午夜精品在线福利| 成人国语在线视频| 国产日韩一区二区三区精品不卡| 午夜日韩欧美国产| 欧美精品亚洲一区二区| 悠悠久久av| 搡老岳熟女国产| 人妻一区二区av| 91av网站免费观看| 一区二区三区精品91| aaaaa片日本免费| 夫妻午夜视频| 亚洲精品av麻豆狂野| 久久久久久人人人人人| 人妻 亚洲 视频| 美女高潮到喷水免费观看| 亚洲国产精品一区二区三区在线| 香蕉国产在线看| 久久久久国产一级毛片高清牌| 夜夜躁狠狠躁天天躁| 欧美 亚洲 国产 日韩一| 欧美日韩黄片免| 久久国产精品男人的天堂亚洲| 亚洲av成人av| 交换朋友夫妻互换小说| 国产欧美日韩综合在线一区二区| 久久久久国产一级毛片高清牌| 欧美人与性动交α欧美软件| 精品一品国产午夜福利视频| 这个男人来自地球电影免费观看| av国产精品久久久久影院| 岛国在线观看网站| 精品电影一区二区在线| 两个人免费观看高清视频| 免费高清在线观看日韩| 黄片大片在线免费观看| 欧美最黄视频在线播放免费 | 国产一区有黄有色的免费视频| 日本精品一区二区三区蜜桃| 午夜精品久久久久久毛片777| 一边摸一边抽搐一进一出视频| 欧美激情极品国产一区二区三区| 久久草成人影院| 精品少妇一区二区三区视频日本电影| 欧美激情极品国产一区二区三区| 最新在线观看一区二区三区| 不卡av一区二区三区| 丰满的人妻完整版| 午夜福利免费观看在线| 亚洲少妇的诱惑av| 中文字幕另类日韩欧美亚洲嫩草| 成人手机av| 亚洲avbb在线观看| 中文字幕高清在线视频| 亚洲五月色婷婷综合| 亚洲国产毛片av蜜桃av| 久久香蕉国产精品| 麻豆av在线久日| 女人被狂操c到高潮| 怎么达到女性高潮| 一级a爱视频在线免费观看| 精品亚洲成国产av| 日韩三级视频一区二区三区| 麻豆av在线久日| 久久草成人影院| 色播在线永久视频| 国产成人欧美| 亚洲成人国产一区在线观看| 最近最新中文字幕大全免费视频| 18禁黄网站禁片午夜丰满| 国产精品一区二区在线不卡| 嫩草影视91久久| 亚洲精品美女久久久久99蜜臀| 老鸭窝网址在线观看| 久久人人爽av亚洲精品天堂| 乱人伦中国视频| 日本a在线网址| 黄片播放在线免费| 国产精品1区2区在线观看. | 精品一区二区三卡| 久久中文字幕人妻熟女| av福利片在线| 91老司机精品| 精品少妇一区二区三区视频日本电影| 久久香蕉国产精品| 精品人妻熟女毛片av久久网站| 精品久久久久久久毛片微露脸| 麻豆国产av国片精品| 国产日韩一区二区三区精品不卡| 下体分泌物呈黄色| 搡老岳熟女国产| 久久人妻福利社区极品人妻图片| 日韩有码中文字幕| 黄片播放在线免费| 欧美成人午夜精品| 国产蜜桃级精品一区二区三区 | 一边摸一边抽搐一进一出视频| 国产日韩欧美亚洲二区| 欧美日韩中文字幕国产精品一区二区三区 | 国产亚洲精品久久久久5区| 久久久国产成人免费| 国产日韩一区二区三区精品不卡| 美女国产高潮福利片在线看| 亚洲国产精品sss在线观看 | 男人舔女人的私密视频| 国产欧美日韩精品亚洲av| 久久性视频一级片| 91九色精品人成在线观看| 黄色女人牲交| 后天国语完整版免费观看| 丰满迷人的少妇在线观看| 麻豆成人av在线观看| 国产av一区二区精品久久| 男女下面插进去视频免费观看| 中文字幕色久视频| 久久精品亚洲精品国产色婷小说| 女人高潮潮喷娇喘18禁视频| 女人久久www免费人成看片| 欧美精品一区二区免费开放| 黑人欧美特级aaaaaa片| 日韩中文字幕欧美一区二区| 18禁国产床啪视频网站| 欧美亚洲 丝袜 人妻 在线| 国产精品永久免费网站| 欧美精品亚洲一区二区| 亚洲片人在线观看| 电影成人av| 久久久久精品国产欧美久久久| 欧美日韩精品网址| 人妻 亚洲 视频| 久久精品91无色码中文字幕| 中文字幕高清在线视频| 两个人免费观看高清视频| 久久久久久久久久久久大奶| 亚洲国产中文字幕在线视频| 欧美日韩av久久| 国产av一区二区精品久久| 国产精品久久视频播放| www.熟女人妻精品国产| 久久香蕉精品热| 人人妻人人爽人人添夜夜欢视频| 成年人黄色毛片网站| 人人妻人人添人人爽欧美一区卜| 黄色 视频免费看| 女人久久www免费人成看片| 无人区码免费观看不卡| 9191精品国产免费久久| 新久久久久国产一级毛片| 国产亚洲精品久久久久5区| 亚洲第一欧美日韩一区二区三区| 欧美精品啪啪一区二区三区| 在线观看免费视频网站a站| 一边摸一边抽搐一进一小说 | 精品一品国产午夜福利视频| 欧美日韩视频精品一区| 国产亚洲一区二区精品| 一级,二级,三级黄色视频| 欧美大码av| 身体一侧抽搐| 久久精品人人爽人人爽视色| 99热网站在线观看| 欧美性长视频在线观看| 一区福利在线观看| 久久久久久久久久久久大奶| 搡老乐熟女国产| 亚洲国产欧美日韩在线播放| 视频区欧美日本亚洲| 国产97色在线日韩免费| 一边摸一边抽搐一进一出视频| 一进一出好大好爽视频| 国产av又大| 亚洲国产欧美网| 亚洲第一av免费看| 成人手机av| 精品国产国语对白av| 国产xxxxx性猛交| 久久久国产一区二区| 9热在线视频观看99| 欧美午夜高清在线| 69精品国产乱码久久久| 亚洲av片天天在线观看| 国产一区有黄有色的免费视频| 一级片免费观看大全| 91老司机精品| 国产精品影院久久| 亚洲av日韩在线播放| 搡老乐熟女国产| 三上悠亚av全集在线观看| 久久中文字幕一级| 多毛熟女@视频| 国产成人欧美| 国产有黄有色有爽视频| 多毛熟女@视频| 巨乳人妻的诱惑在线观看| 一区在线观看完整版| 韩国av一区二区三区四区| 日韩精品免费视频一区二区三区| 国产成人影院久久av| 一级毛片高清免费大全| 久久久国产成人精品二区 | 无限看片的www在线观看| 80岁老熟妇乱子伦牲交| 别揉我奶头~嗯~啊~动态视频| 色播在线永久视频| 国产主播在线观看一区二区| 国产av又大| 成人永久免费在线观看视频| 日韩欧美一区二区三区在线观看 | 国产精品一区二区在线不卡| 天天躁夜夜躁狠狠躁躁| 最近最新中文字幕大全免费视频| 久久国产精品人妻蜜桃| 日韩制服丝袜自拍偷拍| 欧美在线一区亚洲| av国产精品久久久久影院| 国产在线观看jvid| 99香蕉大伊视频| 欧美日韩黄片免| 免费观看人在逋| 精品第一国产精品| 涩涩av久久男人的天堂| 在线观看午夜福利视频| 亚洲成人手机| 在线永久观看黄色视频| 三级毛片av免费| 九色亚洲精品在线播放| 黄色视频不卡| 国产激情久久老熟女| ponron亚洲| 母亲3免费完整高清在线观看| 中文字幕av电影在线播放| 精品视频人人做人人爽| 别揉我奶头~嗯~啊~动态视频| 欧美国产精品va在线观看不卡| 深夜精品福利| 中出人妻视频一区二区| 国产成人欧美| 人成视频在线观看免费观看| 国精品久久久久久国模美| 精品一品国产午夜福利视频| 在线观看免费午夜福利视频| 成年版毛片免费区| 亚洲中文av在线| 高清在线国产一区| 久久精品熟女亚洲av麻豆精品| 亚洲免费av在线视频| 精品国产超薄肉色丝袜足j| 性色av乱码一区二区三区2| 高潮久久久久久久久久久不卡| 天天影视国产精品| 久久久国产成人免费| 成年版毛片免费区| 老司机影院毛片| 人人澡人人妻人| 国产一区在线观看成人免费| 操出白浆在线播放| 一级黄色大片毛片| 精品国产乱码久久久久久男人| 人人妻人人添人人爽欧美一区卜| 日本欧美视频一区| bbb黄色大片| 丰满饥渴人妻一区二区三| 亚洲国产精品一区二区三区在线| 中文字幕av电影在线播放| 国产精品久久久人人做人人爽| 国产亚洲欧美98| 每晚都被弄得嗷嗷叫到高潮| 亚洲五月天丁香| 国产在视频线精品| 久久香蕉国产精品| 91大片在线观看| 欧美在线黄色| aaaaa片日本免费| 伊人久久大香线蕉亚洲五| 热99re8久久精品国产| 超碰97精品在线观看| 电影成人av| 欧美国产精品va在线观看不卡| 成人免费观看视频高清| 欧美久久黑人一区二区| 在线观看免费视频网站a站| 99久久精品国产亚洲精品| 免费在线观看影片大全网站| 中国美女看黄片| 久久久久国产精品人妻aⅴ院 | 人妻丰满熟妇av一区二区三区 | 狠狠狠狠99中文字幕| 王馨瑶露胸无遮挡在线观看| 亚洲精品美女久久久久99蜜臀| 在线播放国产精品三级| 99re在线观看精品视频| 两人在一起打扑克的视频| 三上悠亚av全集在线观看| 高清毛片免费观看视频网站 | 久久久久视频综合| 乱人伦中国视频| 久99久视频精品免费| 妹子高潮喷水视频| 丝袜人妻中文字幕| 亚洲色图av天堂| 国产伦人伦偷精品视频| 美女国产高潮福利片在线看| 欧美精品人与动牲交sv欧美| 超碰成人久久| 久久中文字幕人妻熟女| 免费一级毛片在线播放高清视频 | 亚洲欧洲精品一区二区精品久久久| 国产一区有黄有色的免费视频| x7x7x7水蜜桃| 成人国产一区最新在线观看| 久9热在线精品视频| 动漫黄色视频在线观看| 又大又爽又粗| 国产精品1区2区在线观看. | 亚洲欧美激情综合另类| 亚洲专区中文字幕在线| 免费观看a级毛片全部| 热re99久久精品国产66热6| 免费看十八禁软件| 国产精品98久久久久久宅男小说| 久久99一区二区三区| 超碰成人久久| 免费女性裸体啪啪无遮挡网站| 搡老熟女国产l中国老女人| 午夜福利在线免费观看网站| 男男h啪啪无遮挡| 欧美日韩中文字幕国产精品一区二区三区 | 女警被强在线播放| 女人被躁到高潮嗷嗷叫费观| 如日韩欧美国产精品一区二区三区| 在线观看免费视频日本深夜| 亚洲精品自拍成人| 久热爱精品视频在线9| 久久精品91无色码中文字幕| 九色亚洲精品在线播放| 日韩人妻精品一区2区三区| 国产精品99久久99久久久不卡| 天天操日日干夜夜撸| 曰老女人黄片| 母亲3免费完整高清在线观看| 亚洲色图av天堂| av超薄肉色丝袜交足视频| 国产精品1区2区在线观看. | 成年女人毛片免费观看观看9 | 国产精品一区二区在线观看99| 亚洲成人手机| 俄罗斯特黄特色一大片| 女人爽到高潮嗷嗷叫在线视频| 老司机影院毛片| 午夜福利影视在线免费观看| 中文字幕人妻丝袜制服| 热99国产精品久久久久久7| 亚洲精品中文字幕一二三四区| 欧美国产精品va在线观看不卡| 狠狠狠狠99中文字幕| 亚洲熟女毛片儿| av片东京热男人的天堂| 亚洲精华国产精华精| 黑人巨大精品欧美一区二区mp4| 日韩大码丰满熟妇|