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

        ?

        礦井高密度三維電法成像數(shù)值模擬及應(yīng)用

        2019-08-05 01:45:40萬(wàn)豪豪陳金昊李文喜
        煤礦安全 2019年6期
        關(guān)鍵詞:切片電阻率底板

        萬(wàn)豪豪,朱 魯,陳金昊,李文喜

        (山東科技大學(xué) 地球科學(xué)與工程學(xué)院,山東 青島 266590)

        直流高密度電阻率法在礦井中的應(yīng)用已漸為成熟[1-2],但隨著煤礦開(kāi)采深度的逐漸增加,地質(zhì)條件變得尤為復(fù)雜,其中巖體突水已成為危害礦井安全的重大災(zāi)害之一,為確保井下開(kāi)采工作的安全高效進(jìn)行,科研工作者進(jìn)行了深入研究。2003年,底青云分別運(yùn)用佐迪法和2D積分法進(jìn)行了高密度電阻率成像方面的研究[3]。2005年,為解決工作面采空區(qū)滯后突水問(wèn)題,翟培合開(kāi)發(fā)了底板突水動(dòng)態(tài)監(jiān)測(cè)系統(tǒng),實(shí)現(xiàn)了工作面底板突水監(jiān)測(cè)過(guò)程[4];2008年,施龍青等運(yùn)用三維高密度電法探測(cè)了煤礦巷道底板富水性[5];劉志新、徐志剛和王大慶等[6-7]介紹了電阻率成像技術(shù)的特殊問(wèn)題并對(duì)電阻率成像技術(shù)進(jìn)行了展望;2016年,高衛(wèi)富等使用ANSYS進(jìn)行全空間直流電法異常體正演模擬,精準(zhǔn)地探測(cè)出了異常體的位置[8]。

        1 地電場(chǎng)三維有限元理論

        有限元理論的數(shù)學(xué)基礎(chǔ)是變分原理和剖分插值,將研究區(qū)剖分成有限個(gè)小單元,在每個(gè)單元上,將函數(shù)簡(jiǎn)化成線(xiàn)性函數(shù)、二次函數(shù)或高次函數(shù)。它適用于物理性質(zhì)及參數(shù)復(fù)雜區(qū)域,其求解過(guò)程規(guī)范化,適于解決復(fù)雜數(shù)學(xué)模型問(wèn)題[9]。

        在直流電法數(shù)值模擬時(shí),由于點(diǎn)電源周?chē)妶?chǎng)會(huì)因存在諸多干擾因素使得無(wú)法精確計(jì)算,采用二次場(chǎng)方法解決數(shù)據(jù)的誤差問(wèn)題,使用正方形的網(wǎng)格對(duì)其進(jìn)行劃分,每個(gè)網(wǎng)格節(jié)點(diǎn)的潛在值被設(shè)置為1個(gè)未知數(shù),未知數(shù)被代入拉普拉斯方程進(jìn)行計(jì)算,即將模型的電導(dǎo)率σ細(xì)分為背景電導(dǎo)率σp和異常電導(dǎo)率σs,其分別產(chǎn)生一次場(chǎng)up和二次場(chǎng)us,則總場(chǎng)u為二者疊加:

        三維電場(chǎng)電位的邊值問(wèn)題與下列變分問(wèn)題等價(jià):

        式中:σ為地下介質(zhì)的電導(dǎo)率;r為點(diǎn)電源到邊界點(diǎn)的距離;Г∞為區(qū)域Ω的下邊界;Г為2種介質(zhì)的分界面;n為邊界外法線(xiàn)方向的向量坐標(biāo);u0為正常電位。

        2 模型建立與完善

        在礦井中進(jìn)行三維直流電法勘探時(shí),根據(jù)礦井環(huán)境的特點(diǎn),井下巷道之間大多數(shù)都是通過(guò)聯(lián)絡(luò)巷連接,但是每條聯(lián)絡(luò)巷并不都是在1條直線(xiàn)上,無(wú)法沿用地面勘測(cè)的傳統(tǒng)布線(xiàn)方法,且巷道與聯(lián)絡(luò)巷兩幫的圍巖均可導(dǎo)電,使得頂板與底板之間有著電力聯(lián)系,故并沒(méi)有采用傳統(tǒng)模擬實(shí)驗(yàn)中,將巷道、聯(lián)絡(luò)巷設(shè)定為類(lèi)空氣的超高阻體,而采用井下全空間地電模型[10],將巷道與聯(lián)絡(luò)巷相結(jié)合形成測(cè)量通道。

        巷道設(shè)置于地表深度600 m處,巷道長(zhǎng)度設(shè)置為300 m,位于x坐標(biāo)-100~200 m處,聯(lián)絡(luò)巷長(zhǎng)度設(shè)置為100 m,位于y坐標(biāo)-50~50 m處,總長(zhǎng)為600 m,模型所要測(cè)量的每條邊長(zhǎng)250 m,寬100 m,無(wú)窮遠(yuǎn)邊界長(zhǎng)2 400 m,寬2 400 m,深1 500 m,地質(zhì)模型位置平面圖如圖1。正方形GHKM表示無(wú)窮遠(yuǎn)邊界,ABCD表示模型所要測(cè)量的范圍,ABEF表示巷道和聯(lián)絡(luò)巷的范圍。

        對(duì)所建立的模型進(jìn)行網(wǎng)格剖分[11],沿著x軸剖分的平面數(shù)110個(gè),沿著y軸剖分的平面數(shù)為110個(gè),沿著z軸剖分的平面數(shù)為100個(gè),三維立體空間總共剖分區(qū)塊數(shù)約為1.2×106個(gè)。為了提高正演計(jì)算的精度,在網(wǎng)格剖分過(guò)程中對(duì)所研究的目標(biāo)區(qū)域進(jìn)行加密剖分,巷道,聯(lián)絡(luò)巷網(wǎng)格單元剖分大小為5 m,研究區(qū)域到無(wú)窮遠(yuǎn)邊界處網(wǎng)格剖分逐漸變大,網(wǎng)格剖分大小不均勻。

        圖1 地質(zhì)模型位置平面圖

        供電電極A與測(cè)量電極M分布在ABCD的巷道上,電極之間的距離選擇10 m,電極總數(shù)為61個(gè),電極編排的順序?yàn)?~61,電極測(cè)量的順序?yàn)镃BAD,采用固定斷面連續(xù)測(cè)量掃描方式,測(cè)量深度為A到M距離的一半,測(cè)量點(diǎn)O位于A到M距離的中點(diǎn),二級(jí)裝置AM[12]在三維空間中測(cè)量時(shí),測(cè)量點(diǎn)O的位置如圖2。二極裝置A-M三維測(cè)量散點(diǎn)圖如圖2。

        圖2 二極裝置A-M三維測(cè)量散點(diǎn)圖

        根據(jù)實(shí)際礦井突水隱患建立在層狀介質(zhì)中的地電模型,巷道和聯(lián)絡(luò)巷的電阻率為10 000 Ω·m,圍巖的電阻率為2 000 Ω·m。模型1設(shè)置塊狀低阻異常體尺寸為20 m×20 m×20 m,電阻率值為10 Ω·m,位于2條巷道之間的煤層下方,頂部距離煤層底板20 m;模型2設(shè)置低阻陷落柱異常體模型尺寸為80 m×80 m×20 m,電阻率值為 10 Ω·m,位于 2條巷道之間的煤層下方,頂部距離煤層底板30 m;模型3在模型2的基礎(chǔ)上改變陷落柱在底板的位置,其余設(shè)置同模型1;模型4設(shè)置正方體低阻異常體尺寸為 40 m×40 m×40 m,電阻率值為 1 Ω·m,位于 2條巷道之間的煤層下方,頂部距離煤層底板20 m。

        3 數(shù)值模擬結(jié)果與分析

        3.1 塊狀低阻異常體地電模型

        塊狀低阻異常體地電模型經(jīng)過(guò)有限元數(shù)值模擬計(jì)算后在經(jīng)過(guò)反演得到xz切面圖(圖3)。相比于正常地層明顯出現(xiàn)了低阻異常區(qū),視電阻率值在10~30 Ω·m,異常體的位置和大小已經(jīng)在圖中用正方形線(xiàn)圈標(biāo)出,對(duì)比地電模型中異常體的位置發(fā)生了偏移,大小形態(tài)發(fā)生了改變,但總體來(lái)說(shuō)二極裝置(A M)對(duì)于本次的塊狀低阻異常體的捕捉能力較強(qiáng),分辨低阻異常體的能力較高。

        圖3 塊狀低阻異常體地電模型xz切面圖(y=0)

        塊狀低阻異常體地電模型經(jīng)模擬計(jì)算后的三維切片圖如圖4。其中y=50 m和y=-50 m的xz切片在位于巷道底板處,有2條相互平行的條帶,且條帶所顯示的視電阻率較高,推測(cè)是由于巷道高阻區(qū)的影響所致,無(wú)異常體區(qū)域模擬結(jié)果良好,巷道底板下方綠色區(qū)域的視電阻率值在100 Ω·m左右,與地電模型相符。但異常體x的坐標(biāo)范圍與實(shí)際存在偏差,因此,二極裝置(AM)在對(duì)異常體測(cè)量過(guò)程中,在y方向上測(cè)量效果好,分辨率高。

        圖4 塊狀低阻異常體地電模型xz方向三維切片圖

        3.2 陷落柱低阻異常體地電模型

        通過(guò)陷落柱低阻異常體地電模型xz切面圖(圖5),二極裝置對(duì)本次分辨陷落柱低阻異常體的能力較強(qiáng),斷面圖可以明顯地反映低阻異常區(qū)的大小以及埋深,而且也能夠判定異常體的形態(tài),異常體的位置和大小已經(jīng)在圖中用長(zhǎng)方形線(xiàn)圈標(biāo)出,視電阻率值在10~30 Ω·m,其頂?shù)茁裆钆c實(shí)際模型基本吻合。

        圖5 陷落柱低阻異常體地電模型xz切面圖(y=0)

        由陷落柱低阻異常體地電模型xz方向三維切片圖(圖6)可知,異常體視電阻率值在100 Ω·m左右。從低阻區(qū)的視電阻率范圍可推測(cè)異常體的y坐標(biāo)范圍為0~15 m,與實(shí)際地電模型中y坐標(biāo)的范圍0~20 m基本一致,異常體x的坐標(biāo)范圍與實(shí)際模型一致。因此,二極裝置(A-M)在對(duì)異常體測(cè)量過(guò)程中,在y方向上測(cè)量效果好,分辨率高。中間切片的地層,除了相對(duì)低阻區(qū)外,與正常地層和塊狀異常體的視電阻率基本相同。

        圖6 陷落柱低阻異常體地電模型xz方向三維切片圖

        3.3 柱狀低阻異常體地電模型

        柱狀低阻異常區(qū)的大小和埋深與塊狀低阻異常體地電模型保持一致,異常體的長(zhǎng)度基本和地電模型中的一樣,其頂?shù)茁裆钆c實(shí)際模型基本吻合。柱狀低阻異常體地電模型xz切面圖如圖7。但低阻體的左邊有明顯的上移,視電阻率值在20~50 Ω·m,其相較于陷落柱低阻異常體視電阻率值偏大,出現(xiàn)這種現(xiàn)象的原因有可能是異常體的中心坐標(biāo)發(fā)生偏移,離巷道較近,受到巷道高阻體的影響。

        從柱狀低阻異常體地電模型xz方向三維切片圖(圖8)觀察可知,低阻區(qū)的視電阻率范圍可推測(cè)異常體的y坐標(biāo)范圍為-25~15 m,與實(shí)際地電模型中y坐標(biāo)的范圍0~20 m基本一致。

        圖7 柱狀低阻異常體地電模型xz切面圖(y=0)

        圖8 柱狀低阻異常體地電模型xz方向三維切片圖

        3.4 正方體低阻異常體地電模型

        正方體低阻異常體地電模型xz切面圖如圖9。正方體低阻體模擬結(jié)果可以反映低阻異常區(qū)的大小和埋深,視電阻率值在1~30 Ω·m,圖中的空白區(qū)域視電阻率值小于10 Ω·m。受低阻異常體的影響,主要異常區(qū)域集中在異常體的上部,地面下660 m處是異常區(qū)域的分界面,但660 m下方異常反應(yīng)在實(shí)際成像中并沒(méi)有出現(xiàn),主要是受探測(cè)深度的影響,探測(cè)的越深,數(shù)據(jù)量越少,分辨率越低,因此,在700 m以下的位置視電阻率值并不準(zhǔn)確。

        圖9 正方體低阻異常體地電模型xz切面圖(y=0)

        從正方體低阻異常體地電模型xz方向三維切片圖(圖10)可知,從低阻區(qū)的視電阻率范圍可推測(cè)異常體的y坐標(biāo)范圍為-25~25 m,與實(shí)際地電模型中y坐標(biāo)的范圍-20~20 m基本一致,在y=0方向切片圖中低阻異常區(qū)的范圍大于實(shí)際地電模型,這可能是由于低阻異常體的背景電阻率值非常小,電阻率值為1 Ω·m所致。

        圖10 正方體低阻異常體地電模型xz方向三維切片圖

        4 工程應(yīng)用

        某礦1302上工作面(3煤)位于山西組中下部,3煤層開(kāi)采過(guò)程中深受頂板砂巖水害的威脅,為保證1302上工作面開(kāi)采的順利、安全進(jìn)行,對(duì)該礦的1302上工作面進(jìn)行井下高密度三維電法勘探。要求頂板最大探測(cè)高度為頂板上150 m,構(gòu)建頂板砂巖地層電阻率三維數(shù)據(jù)體,應(yīng)用切片技術(shù)圈定頂板砂巖地層的富水區(qū)域,結(jié)合具體鉆探情況指導(dǎo)1302上工作面開(kāi)采過(guò)程中的防治水工作。

        本次井下高密度三維電法勘探數(shù)據(jù)采集使用的是WDJD-4型高密度電阻率系統(tǒng),采用A-M二極裝置形式的三維電法采集方法,采集參數(shù)中儀器所使用的電極道數(shù)是60道,測(cè)量的剖面數(shù)為30個(gè),采用滾動(dòng)的測(cè)量方式,2個(gè)相鄰電極之間的道間距為10 m,供電電壓選擇180 V的直流電壓。

        把RES3DINV反演軟件所得出的三維立體成像圖數(shù)據(jù)輸出,然后在用Slicer Dicer和Voxler對(duì)三維立體成像圖進(jìn)行多方向的三維切片,即為最終的成像成果圖。1302上工作面頂板砂巖地層電阻率三維數(shù)據(jù)體原點(diǎn)位于運(yùn)輸巷與1聯(lián)絡(luò)巷向交匯處,沿運(yùn)輸巷向里為y、沿1聯(lián)絡(luò)巷向下為正方向x,高度方向?yàn)榇怪钡貙酉蛏?,本次勘探建立?302上工作面外側(cè)部分頂板砂巖地層電阻率三維數(shù)據(jù)體長(zhǎng)度750 m、寬度140 m、高度280 m。

        1302上工作面頂板地層電阻率三維數(shù)據(jù)體如圖11,1302上工作面頂板地層電阻率三維數(shù)據(jù)體切片技術(shù)如圖12,1302上工作面頂板地層電阻率三維數(shù)據(jù)體垂層切片如圖13。

        根據(jù)數(shù)據(jù)整體(圖11)和數(shù)據(jù)切片(圖12)可知,低阻區(qū)域明顯,對(duì)比圖13,可圈定3個(gè)低阻區(qū)進(jìn)行解釋。第1區(qū)域位于橫向起始測(cè)點(diǎn)到280 m,頂板向上至高度180 m;第2區(qū)域位于橫向320~400 m之間,頂板向上至高度90 m處;第3區(qū)域位于橫向400~750 m之間,頂板向上至高度120 m處,以上3個(gè)區(qū)域,低阻情況明顯,推測(cè)為含水區(qū)域。

        圖11 1302上工作面頂板地層電阻率三維數(shù)據(jù)體

        圖12 1302上工作面頂板地層電阻率三維數(shù)據(jù)體切片技術(shù)

        圖13 1302上工作面頂板地層電阻率三維數(shù)據(jù)體垂層切片(平行切眼、50 m間距)

        5 結(jié)論

        1)使用Visual Fortran編程軟件,對(duì)一定深度范圍內(nèi)所建立的三維全空間地電模型進(jìn)行正演數(shù)值模擬,反演后的結(jié)果能夠準(zhǔn)確反應(yīng)出所建地電模型異常體區(qū)域。Surfer軟件繪成的ρs視電阻率斷面圖,在一定的深度范圍內(nèi),能大體上反映異常體的形態(tài)、位置和大小,分辨率高,基本上與地電模型吻合。

        2)三維直流電法在煤礦中的應(yīng)用實(shí)例,圈定了3煤層頂板砂巖地層富水區(qū),后期經(jīng)過(guò)鉆探驗(yàn)證了富水區(qū)的存在,針對(duì)性強(qiáng),極大地提高了鉆探的效率,從生產(chǎn)實(shí)際角度論證了高密度電法的可行性,為指導(dǎo)礦井水防治補(bǔ)充了新的方法。

        猜你喜歡
        切片電阻率底板
        基于SDN與NFV的網(wǎng)絡(luò)切片架構(gòu)
        三維電阻率成像與高聚物注漿在水閘加固中的應(yīng)用
        軟土地基上超長(zhǎng)底板節(jié)制閘沉降及應(yīng)力分析
        腎穿刺組織冷凍切片技術(shù)的改進(jìn)方法
        底板巷一巷兩用的回采工作面防火實(shí)踐
        底板隔水層破壞的力學(xué)模型及破壞判據(jù)
        隨鉆電阻率測(cè)井的固定探測(cè)深度合成方法
        冰凍切片、快速石蠟切片在中樞神經(jīng)系統(tǒng)腫瘤診斷中的應(yīng)用價(jià)值比較
        海洋可控源電磁場(chǎng)視電阻率計(jì)算方法
        粉煤灰摻量對(duì)水泥漿體電阻率與自收縮的影響
        乱色欧美激惰| 狠狠亚洲婷婷综合久久久 | 中文乱码人妻系列一区二区| 一本一道AⅤ无码中文字幕| 亚洲av天堂在线免费观看| 7777精品久久久大香线蕉| 国产一精品一aⅴ一免费| 91精品亚洲熟妇少妇| 亚洲成a人片在线观看无码专区| 丰满少妇大力进入av亚洲| 日本a在线免费观看| 亚州中文字幕乱码中文字幕| 国产av一区二区三区无码野战| 亚洲国产无套无码av电影| 中文字幕天天躁日日躁狠狠| 精品国产中文久久久免费| 国产欧美成人一区二区a片| 丰满爆乳无码一区二区三区| 久久久精品国产亚洲麻色欲| 风流熟女一区二区三区| 性无码专区无码| 欧美巨大精品欧美一区二区| 日本不卡的一区二区三区| 日韩女同视频在线网站| 18禁裸男晨勃露j毛免费观看 | 人妻少妇喷水意淫诱惑| 亚洲日本一区二区三区四区| 成人免费xxxxx在线观看| 亚洲综合中文字幕乱码在线| 亚洲一区二区三区新视频| 亚洲av成人片色在线观看| 久久午夜无码鲁丝片直播午夜精品| 亚洲AV无码日韩综合欧亚| 最新国产女主播在线观看| 久久超碰97人人做人人爱| 精品无码久久久久久久久粉色| 不卡av一区二区在线| 亚洲欧美中文字幕5发布| 亚洲成a人片在线| 久久网站在线免费观看| 亚洲国产精品综合久久网络|