亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        近水平煤層開采地表移動邊界的預(yù)測研究

        2023-10-07 13:07:50孫學(xué)陽袁官思龍苗霖田杜榮軍
        煤礦安全 2023年9期
        關(guān)鍵詞:等值線圖曲率邊界

        孫學(xué)陽 ,袁官思龍 ,苗霖田 ,杜榮軍 ,程 正

        (1.西安科技大學(xué) 地質(zhì)與環(huán)境學(xué)院,陜西 西安 710054;2.陜西省煤炭綠色開發(fā)地質(zhì)保障重點實驗室,陜西 西安 710054;3.煤炭資源勘查與綜合利用重點實驗室,陜西 西安 710054)

        獲取采煤沉陷特征值是研究煤礦開采地表移動變形規(guī)律和揭示其形成機(jī)理的基礎(chǔ),現(xiàn)行規(guī)范和規(guī)程[1-2]根據(jù)采煤沉陷特征值對開采沉陷盆地邊界和損毀程度做了界定。理論分析方、相似材料模擬實驗、數(shù)值模擬是目前地表移動和變形預(yù)測的3 種方法。理論分析受到相關(guān)變量的限制,相似材料模擬實驗受到實驗條件的限制,20 世紀(jì)60年代R W Clough 首次將有限元引入土石壩的穩(wěn)定性分析中,從此數(shù)值模擬技術(shù)成為了巖土工程領(lǐng)域主流的研究和設(shè)計方法[3]。因此對于采煤沉陷特征值的獲取和地表移動邊界的預(yù)測研究中,數(shù)值模擬具有優(yōu)勢。

        地表損毀邊界的主要影響因子為下沉、傾斜、水平變形。胡振琪等[4]以理論分析結(jié)合現(xiàn)場調(diào)查,提出了以耕地?fù)p毀臨界下沉值、臨界水平變形值和臨界傾斜值中的最小值圈定耕地?fù)p毀范圍;鄧偉男[5]利用FISH 語言代碼,提取模型主斷面上的下沉、水平移動、傾斜、水平變形及曲率值來繪制曲線圖,為數(shù)值模擬中數(shù)據(jù)的處理提供了參考;文獻(xiàn)[6-8]采用FLAC3D軟件進(jìn)行模擬開采時,通過在走向和傾向主斷面上布設(shè)監(jiān)測點的方式獲得其采煤沉陷特征值;LI 等[9]從應(yīng)力場和位移場的角度出發(fā),結(jié)合地表影響范圍的放大比,分析了地表沉降影響區(qū)隨采空區(qū)不斷擴(kuò)大而發(fā)生的動態(tài)變化;WANG 等[10]研究了在FLAC3D建立含復(fù)雜地層邊界的三維數(shù)值計算模型的方法,并通過MATLAB 軟件導(dǎo)出地表沉降和水平移動數(shù)據(jù),進(jìn)而由計算公式轉(zhuǎn)化出傾斜、水平變形和曲率的等值線圖來直觀反映地表移動變形的影響范圍;蔡音飛等[11]利用Mathematica 軟件對FLAC3D模擬結(jié)果進(jìn)行再次分析得到地表移動變形量,并通過列表繪圖函數(shù)ListPlot 繪制各次試驗的變形圖對地表受損情況進(jìn)行評價。

        基于此,從沉陷盆地全面積的角度,運(yùn)用FISH 語言提取計算采煤沉陷特征值并進(jìn)行運(yùn)用,圈定了霍寶干河煤礦多工作面開采后地表移動盆地最外邊界、地表危險移動邊界、地表移動裂縫邊界;研究結(jié)果不僅為開采沉陷數(shù)值模擬數(shù)據(jù)后處理提供了新思路,還為生態(tài)恢復(fù)治理規(guī)劃提供參考。

        1 地表移動變形數(shù)據(jù)的提取

        地下巖土體的移動變形是采煤沉陷問題的核心,是地面沉陷的原因[12]。沉陷就是移動變形的總體描述,地表沉陷狀態(tài)按移動變形性質(zhì)、方向、方式的不同進(jìn)行描述,可分為垂直方向的移動變形和水平方向的移動變形。工程中通常采用下沉、傾斜、曲率,水平移動、水平變形等5 個采煤沉陷特征值進(jìn)行定量評價[13]。無論是實測數(shù)據(jù)分析還是形變預(yù)計,采煤沉陷特征值都是基于下沉和水平位移2 個移動分量[14]。因此先將水平移動分量和垂直分量提取出來,再進(jìn)行相應(yīng)計算就能得到傾斜、水平變形和曲率。

        在FLAC3D軟件中,力作用在單元中心坐標(biāo)上,位移產(chǎn)生在網(wǎng)格節(jié)點上。針對這一特征,利用FISH 語言對六面體網(wǎng)格主導(dǎo)的計算模型中各個網(wǎng)格單元的特定節(jié)點進(jìn)行遍歷,進(jìn)行數(shù)據(jù)的提取并保存為Surfer 能夠識別的數(shù)據(jù)格式。

        除此之外,F(xiàn)LAC3D中提供了豐富的內(nèi)置函數(shù),用戶可以運(yùn)用它們進(jìn)行簡單的程序設(shè)計,并得到相應(yīng)結(jié)果。針對采煤沉陷問題,有必要先對FLAC3D中模型網(wǎng)格的節(jié)點編號和單元編號進(jìn)行一下梳理,同時也是提取數(shù)據(jù)的代碼的編寫依據(jù)。

        為了便于闡述,只建立了3× 3× 1 的模型,利用內(nèi)置函數(shù)gp.pos(gp)、data.label.create(v)、data.label.text(lp)可以將單元編號和每個單元的節(jié)點編號在FLAC3D中顯示出來。FLAC3D中單元和節(jié)點編號示意圖如圖1,圖中:紅色數(shù)字為每個單元的編號;藍(lán)色數(shù)字為1 號單元的節(jié)點編號。規(guī)律如下:單元編號在走向上(x方向)依次排列,在傾向上(y方向)按等差數(shù)列排列;節(jié)點編號在單個單元上遵循右手螺旋法則。

        圖1 FLAC3D 中單元和節(jié)點編號示意圖Fig.1 Schematic diagrams of unit and node number in FLAC3D

        在FLAC3D中,對特定范圍的指定有2 種方式,自定義分組和自動分面;遍歷的對象也有2 種,一種是對單元進(jìn)行遍歷,另外一種是對節(jié)點進(jìn)行遍歷。

        對于下沉值和水平移動值僅對節(jié)點進(jìn)行遍歷即可;傾斜值和水平變形值先對單元進(jìn)行遍歷,然后在具體編號的單元內(nèi)通過節(jié)點指針找到相應(yīng)編號節(jié)點的坐標(biāo)和位移,曲率需要同時在2 個單元內(nèi)通過節(jié)點指針找到相應(yīng)編號節(jié)點的坐標(biāo)和位移。最后,通過FISH 語言編寫出傾斜、水平變形和曲率的函數(shù)來提取數(shù)據(jù)并進(jìn)行文本文件的保存,為繪制等值線提供數(shù)據(jù)基礎(chǔ)。

        1.1 下沉值數(shù)據(jù)提取及保存

        采用自動分面命令,模型表面將自動被命名為“Top i”,這就是下沉值的遍歷范圍。垂直移動分量為負(fù),由下沉值正負(fù)號規(guī)定可知:正值表示測點下沉,負(fù)值表示測點上升,所以要在變量前面添加負(fù)號。另外,F(xiàn)LAC3D的默認(rèn)單位為m,下沉值的單位需要轉(zhuǎn)換成mm,具體代碼如下:

        1.2 水平移動值數(shù)據(jù)提取及保存

        水平移動值遍歷范圍同上,它的正負(fù)號規(guī)定:指向礦層走向斷面走向方向的移動為正,反之為負(fù);指向傾斜斷面上山方向的移動為正,下山方向的為負(fù)。水平移動值涉及傾向和走向兩個方向的分量,在FISH 變量中分別用dx 和dy 表示。單位為mm,具體代碼如下:

        1.3 傾斜值數(shù)據(jù)提取及保存

        傾斜值的計算涉及同一單元的相鄰兩節(jié)點,根據(jù)每個單元的節(jié)點編號規(guī)則,對提取傾斜值的代碼設(shè)計如下:走向上傾斜值為每個單元的8 號節(jié)點和6 號節(jié)點(或7 號節(jié)點和4 號節(jié)點)的下沉值之差與它們之間沿x軸水平距離的比值;傾向上傾斜值為8 號節(jié)點和7 號節(jié)點(或6 號節(jié)點和4號節(jié)點)的下沉值之差與它們之間沿y軸水平距離的比值[15]。

        傾斜值反映了地表移動盆地沿某一方向的坡度,通常以i 表示,下列代碼中ix 代表走向上的傾斜值,iy 代表傾向上的傾斜值。它的正負(fù)號規(guī)定:指向礦層走向斷面走向方向的傾斜為正,反之為負(fù);指向傾斜斷面上山方向的傾斜為正,下山方向的為負(fù)。單位需轉(zhuǎn)換為mm/m,具體代碼如下(以走向為例):

        1.4 水平變形值數(shù)據(jù)提取及保存

        水平變形簡言之就是兩點間單位長度上的水平移動差,用相鄰兩點的水平移動差與兩點間水平距離的比值來表示。如圖1(b),走向上水平變形同樣用8 號和6 號節(jié)點(或7 號和4 號節(jié)點)的x軸位移分量之差與它們之間沿x軸水平距離的比值表示;傾向上水平變形用6 號和4 號節(jié)點(或8號和7 號節(jié)點)的y軸位移分量之差與它們之間沿y軸方向水平距離的比值表示。豎直方向的變化發(fā)生傾斜,水平方向的變化發(fā)生水平變形。所以參考提取傾斜值的代碼,將內(nèi)置函數(shù)gp.disp.z()在走向上改成gp.disp.x(),在傾向上改成gp.disp.y(),節(jié)點編號和坐標(biāo)函數(shù)對應(yīng)修改即可。它的正負(fù)號規(guī)定:單位長度線段的拉伸變形為正,壓縮變形為負(fù)[16],單位為mm/m。

        1.5 曲率值數(shù)據(jù)提取及保存

        曲率反映了地表垂直斷面上觀測線的彎曲程度,它可以用相鄰線段的傾斜差與兩線段中點間的水平距離之比來表示。

        如圖1,走向上的曲率為相鄰單元8 號節(jié)點和6 號節(jié)點連成線段的傾斜值之差與兩線段x軸中點間水平距離的比值;傾向上的曲率為相鄰單元的6 號節(jié)點和4 號節(jié)點連成線段的傾斜值之差與兩線段y軸中點間水平距離的比值。它的正負(fù)號規(guī)定:上凸表示地表下沉曲線為正,下凹則為負(fù)。單位需轉(zhuǎn)換成mm/m2,具體代碼如下(以傾向為例):

        其中a、b、c、d為網(wǎng)格單元編號,e為沿走向方向的網(wǎng)格單元個數(shù),g=a-e。

        2 近水平煤層開采沉陷模擬

        霍寶干河煤礦為近水平煤層,第四系黃土平均厚100 m,屬于厚松散層,煤層平均厚度為4.5m,開采深度為460~500 m。

        模型尺寸為2 000 m ×2 050 m ×600 m,共分為8 層,平均采深為480 m,采厚4.5 m,煤層傾角為0°,共計116 850 個單元,132 396 個節(jié)點;數(shù)值計算模型分為3 個工作面進(jìn)行開采模擬,分別是1101 首采工作面,1102 工作面和1103 工作面。3 個工作面尺寸均為600 m ×210 m,工作面之間均留設(shè)10 m 寬的護(hù)巷煤柱。模型剖切圖如圖2,巖層物理力學(xué)參數(shù)見表1。

        表1 巖層物理力學(xué)參數(shù)Table 1 Physical and mechanical parameters of rock strata

        由于地表移動觀測站只設(shè)置在1101 首采工作面,為了保證3 個工作面全采后數(shù)值模擬結(jié)果的正確性以及用FISH 語言提取特征值的合理性,有必要先將該工作面下沉盆地主斷面的特征值進(jìn)行提取,1101 工作面主斷面特征值如圖3~圖5。

        圖4 1101 工作面主斷面傾斜值和水平變形值Fig.4 Inclination values and horizontal deformation values of main section of 1101 working face

        圖5 1101 工作面主斷面曲率值Fig.5 Curvature values of main section of 1101 working face

        將圖3~圖5 中各特征值曲線的極值進(jìn)行統(tǒng)計,1101 工作面地表移動變形量極值見表2。

        表2 1101 工作面地表移動變形量極值Table 2 Extreme values of surface movement and deformation in 1101 working face

        從表2 中可以看出用FISH 語言提取的各特征值與現(xiàn)場實測資料相吻合,驗證了該數(shù)值模擬的正確性及用FISH 語言提取地表移動變形數(shù)據(jù)的合理性。

        3 應(yīng)用實例

        3.1 Surfer 繪圖流程

        Surfer 軟件的主要優(yōu)勢是利用其插值功能進(jìn)行等值線圖的繪制。用它來繪制等值線圖的前提是數(shù)據(jù)文件格式要轉(zhuǎn)換成“.grd”文件格式[17]。

        由于Surfer 的插值功能,即使數(shù)據(jù)不是等間距的,它依然能夠作圖。如果不進(jìn)行數(shù)據(jù)的插值,只是繪制原始數(shù)據(jù)的圖,可以在“.text”轉(zhuǎn)“.grd”的對話框中選擇距離平方反比法或克里金插值法,本實例所有等值線圖的繪制過程皆采用克里金插值法。

        運(yùn)行提取采煤沉陷特征值的各段FISH 代碼會得到相應(yīng)數(shù)據(jù)文件,按上述步驟操作之后新建等值線圖即可,按照項目的要求進(jìn)行加密、著色等設(shè)置。除此之外,Surfer 中還可以添加邊界文件,可用于工作面輪廓線的繪制,為測量地表各邊界與采空區(qū)邊界線間的距離提供了便利。

        3.2 地表移動邊界成圖

        地表必然受到煤炭開采的影響,提前預(yù)判這一影響的1 個重要環(huán)節(jié)就是對變形后的地表進(jìn)行等值線的繪制[18]。實測法獲取地表移動盆地角值參數(shù)存在困難,建立地表移動觀測站成本高,觀測數(shù)據(jù)匱乏;類比法獲取角值參數(shù)又受到地質(zhì)采礦條件的制約;統(tǒng)計模型法通過統(tǒng)計經(jīng)驗公式獲取角值參數(shù)的方法又受到實測數(shù)據(jù)和地質(zhì)采礦條件的雙重制約。

        綜上所述,隨著1101 工作面、1102 工作面和1103 工作面的開采,將整個模型通過FISH 代碼提取采煤沉陷數(shù)據(jù)并繪制等值線圖,從數(shù)值模擬的角度預(yù)測近水平煤層多工作面開采后地表移動邊界。地表下沉等值線圖如圖6,地表移動邊界等值線圖如圖7,圖中:藍(lán)色線框代表的是工作面輪廓,紅色和紫色等值線為臨界值等值線。

        圖6 地表下沉等值線圖Fig.6 Contour map of surface subsidence

        圖7 地表移動邊界等值線圖Fig.7 Contour diagrams of surface movement boundary

        3.2.1 地表移動盆地最外邊界

        地表移動盆地最外邊界以下沉10 mm 確定,該邊界在走向方向上距離采空區(qū)邊界518 m,在傾向方向上距離采空區(qū)邊界502 m。該邊界以內(nèi)的區(qū)域都屬于開采沉陷影響區(qū),這個邊界以外的區(qū)域一般認(rèn)為不受開采沉陷的影響。

        3.2.2 地表危險移動邊界

        地表危險移動邊界按照現(xiàn)行規(guī)范[1]中對磚混結(jié)構(gòu)建筑物臨界變形值的規(guī)定:傾斜變形值i=±3 mm/m,水平變形值ε=±2 mm/m,曲率K=±0.2 mm/m2來確定。

        圖7(a)和圖7(b)中的傾斜臨界值等值線為紅色橢圓線,圖7(c)和圖7(d)走向上曲率最大值為0.045 6,傾向上曲率最大值為0.036 4 均未達(dá)到曲率臨界值,因此未能繪制出曲率的臨界等值線,圖7(e)和圖7(f)中紫色等值線為水平變形臨界值等值線;由此可以看出走向上地表危險移動邊界與采空區(qū)邊界左側(cè)相距299 m,右側(cè)相距291 m,傾向上的地表危險移動邊界與采空區(qū)邊界上側(cè)相距308 m,下側(cè)相距316 m。上下左右間距不相等的原因是由于工作面推進(jìn)方向和工作面開采順序不同造成的。地表破壞的程度和范圍在等值線圖上的直觀反映就是傾斜和曲率大的區(qū)域[19]。

        3.2.3 地表移動裂縫邊界

        開采沉陷引起的土體拉伸破壞臨界值與土體變形值的定量關(guān)系如式(1)[20]:

        式中: εxc為 土體產(chǎn)生裂縫的水平變形臨界值,mm/m;h為 土體埋深,m; ρ為平均密度,t/m3;φ為摩擦角,(°); μ為泊松比;E為彈性模量,MPa;c為黏聚力,MPa。

        εxc隨土體的埋深增加而增大,土體中任意一點的水平拉伸變形值εx≥εxc時,該處土體產(chǎn)生拉伸破壞,否則僅產(chǎn)生連續(xù)的沉陷變形。地表(h=0)處土體破壞極限狀態(tài)的水平變形臨界值εh=0為最小,可通過式(2)計算。

        將黃土的物理力學(xué)參數(shù)E=10 MPa、μ=0.3、c=0.02 MPa、φ=15°代入式(2)得到地表產(chǎn)生拉伸變形的水平變形臨界值為εh=0=2.15 mm/m。

        將該臨界值的等值線繪制于圖7(e)和圖7(f)中,以紅色等值線標(biāo)示出,可以看出,走向上左側(cè)裂縫邊界距離采空區(qū)邊界498 m,右側(cè)裂縫邊界距離采空區(qū)邊界484 m;傾向上上側(cè)裂縫邊界距離采空區(qū)邊界476 m,下側(cè)裂縫邊界距離采空區(qū)邊界493 m。同樣,上下左右間距不相等的原因是由于工作面推進(jìn)方向和工作面開采順序不同造成的。紅色等值線不僅確定了裂縫邊界,而且圈定了產(chǎn)生裂縫區(qū)域的范圍,如圖7(e)和圖7(f)。

        4 結(jié) 語

        1)依據(jù)單元和節(jié)點的編號規(guī)則,繪制曲率等值線圖的數(shù)據(jù)需要多個循環(huán)控制語句進(jìn)行提取,F(xiàn)ISH 代碼能夠在FLAC3D中實現(xiàn)采煤沉陷特征值的數(shù)據(jù)提取。

        2)對霍寶干河煤礦進(jìn)了數(shù)值模擬,提取1101工作面的采煤沉陷特征值曲線,結(jié)果與現(xiàn)場實測值吻合較好,驗證了數(shù)值模擬結(jié)果的可靠性。并提取了該工作面走向上的曲率值+0.007 mm/m2和-0.014 mm/m2,傾向上的曲率值+0.01 mm/m2和-0.03 mm/m2。

        3)工作面推進(jìn)方向和開采順序的不同會導(dǎo)致地表移動盆地邊界距離采空區(qū)邊界左右間距不相等。4)通過FISH 語言編程提取3 個工作面開采后地表移動變形數(shù)據(jù),繪制的等值線圖預(yù)測了地表移動盆地最外邊界、地表危險移動邊界、地表移動裂縫邊界和裂縫區(qū)域,從數(shù)值模擬的角度為預(yù)測地表移動盆地邊界提供了新思路。

        猜你喜歡
        等值線圖曲率邊界
        湖南省似大地水準(zhǔn)面模型對比研究
        大曲率沉管安裝關(guān)鍵技術(shù)研究
        一類雙曲平均曲率流的對稱與整體解
        拓展閱讀的邊界
        半正迷向曲率的四維Shrinking Gradient Ricci Solitons
        論中立的幫助行為之可罰邊界
        如何來解決等值線問題
        利用ACAD線型編輯功能實現(xiàn)SURFER繪制等值線圖的修改初探
        江蘇水利(2015年9期)2015-12-12 08:17:18
        “偽翻譯”:“翻譯”之邊界行走者
        Esn+1中具有至多兩個不同主曲率的2-調(diào)和超曲面
        在线播放偷拍一区二区| 300部国产真实乱| 人妻体内射精一区二区三区| 朝鲜女人大白屁股ass| 国产精品美女| 女同性恋亚洲一区二区| 亚洲国产精品久久无人区| 国产成人午夜无码电影在线观看| 日本老熟欧美老熟妇| 音影先锋色天堂av电影妓女久久| 久久久精品少妇—二区| 国产精品永久久久久久久久久| 国产精品麻豆成人av电影艾秋| 亚洲国产精品久久久性色av| 精品人妻日韩中文字幕| 最新欧美精品一区二区三区| 无码人妻黑人中文字幕| 一区二区三区放荡人妻| 中文字幕日韩精品永久在线| 久久精品国产亚洲av香蕉| 永久黄网站色视频免费| 国产喷白浆精品一区二区| 91精品国产综合久久久密臀九色 | 久久久精品国产三级精品| 亚洲妇熟xxxx妇色黄| 欧美白人最猛性xxxxx| 精品亚洲不卡一区二区| 亚洲丝袜美腿在线视频| 欧美一区二区三区红桃小说| 人妻丰满av无码中文字幕| 亚洲发给我的在线视频| 国产a在亚洲线播放| 99久热re在线精品99 6热视频| 精品亚洲不卡一区二区| 免费一级淫片日本高清 | av永远在线免费观看| 久久精品中文字幕女同免费| 国产亚洲av综合人人澡精品| 亚洲伦理一区二区三区| 中文字幕一区二区三区在线看一区| 97人妻人人揉人人躁九色|