陜西榆林市杭來灣煤礦三維地下水數(shù)值模擬

陜西榆林市杭來灣煤礦三維地下水數(shù)值模擬

1海陽大學井段預測頂板流量是一項相對復雜、困難的工作。國內(nèi)外生產(chǎn)實踐中積累了幾種方法來提高礦山流量的精度和準確性，但需要改進的精度和準確性。杭來灣煤礦位于陜西省榆林市北部，距榆林城北約20km。該礦可采儲量7.06億t，生產(chǎn)能力800萬t/a，服務年限68a。礦井采用斜井多水平開拓，分層大采高綜采。礦井按煤組設(shè)2個水平，其中，一水平開采3、3-1、4-1和5號煤層，水平標高+1010m，大巷沿3號煤層布置；二水平開采7、8號煤層，水平標高+860m，大巷沿8號煤層布置。礦井首采盤區(qū)為301盤區(qū)，首采煤層為3號煤層。2頂板冒落帶與第四系水害溝通根據(jù)井田水文地質(zhì)勘探資料顯示，井田地下水類型分為第四系松散層孔隙潛水和碎屑巖類裂隙水2類。通過實測裂采比計算3號煤層頂板冒落帶導水裂隙帶高度，當采用分層開采時，上分層累計采厚為4.80m，經(jīng)計算得到冒落帶高度為26.88m，導水裂隙帶高度為100.51m，而3號煤層上覆侏羅紀中統(tǒng)基巖厚度在67.51～211.13m，因此采空區(qū)冒裂隙帶一般不會與第四系潛水溝通；當按煤層采全高計算3號煤層頂板冒落帶導水裂隙帶高度，礦井部分地段采空區(qū)冒裂隙帶將與第四系潛水溝通，存在第四系水害影響。杭來灣煤礦井下涌水量主要來源于3號煤層頂板砂巖含水層。據(jù)調(diào)查，在3號煤層回采工作面、各個聯(lián)絡巷、3處斜井及掘進工作面等均有不同程度的淋水現(xiàn)象。依據(jù)井下涌水量實測數(shù)據(jù)，礦井月平均涌水量總體呈逐月上升的趨勢。3預測開采井的水量數(shù)值法預測礦井涌水量是采用三維地下水模擬軟件Modflow，該軟件是一款模塊化地下水模擬軟件，可進行三維地下水流模擬。3.1地質(zhì)條件的概念模型和數(shù)學描述3.1.1第四系潛水含水層由于導水裂隙帶之上的松散層潛水越流進入坑道，杭來灣煤礦先期開采地段3號煤層的礦坑涌水量預測模型區(qū)在垂向上應該包含3號煤層及其之上的巖土層（含水層和弱透水層），在平面上盡可能包含井田區(qū)松散層潛水較為完整的水文地質(zhì)單元。杭來灣井田分屬榆溪河流域，由于第四系潛水與地表水的分水嶺基本一致，所以地表水與松散層潛水均是按溝域形成、分布與運移的。在區(qū)域上，第四系潛水含水層的主體為上更新統(tǒng)薩拉烏蘇組砂層，其次為中上更新統(tǒng)黃土層。第四系潛水含水層的富水性具有西北部地區(qū)較強、東南部地區(qū)較弱的特點，而杭來灣井田位于中等富水性區(qū)域。綜合考慮區(qū)域水文地質(zhì)特征和溝域水文地質(zhì)單元的相對完整性以及井田邊界與模型邊界的適度距離（不宜太近與太遠）等因素，確定將水源取水影響范圍作為數(shù)值模型有效區(qū)的范圍。3.1.2第四系黃土地下水流結(jié)構(gòu)根據(jù)對評價區(qū)水文地質(zhì)條件、煤田開采范圍以及周圍水環(huán)境特征等綜合分析，確定本次地下水資源評價范圍：在水平方向上，北以公里網(wǎng)4277000為界，南至4240000，西以19410000為界，東至19380000，平面面積1110km評價區(qū)內(nèi)大部分地區(qū)被風積沙所覆蓋，部分地段有第四系黃土。其自上而下為第四系上更新統(tǒng)薩拉烏蘇組及黃土、侏羅系泥砂巖，組成2～3層結(jié)構(gòu)。為了便于地下水流數(shù)值模型計算，評價區(qū)垂向上統(tǒng)一為5層結(jié)構(gòu)，即第四系松散含水層、導水裂隙帶之上原基巖、導水冒裂帶、3號煤層采空區(qū)；水平方向包含井田所在的榆溪河流域區(qū)，從而構(gòu)建起模型區(qū)的5層結(jié)構(gòu)：第一結(jié)構(gòu)層，以上更新統(tǒng)薩拉烏蘇組為主體的孔隙含水層，厚度在0～64.92m，成片分布于煤礦現(xiàn)代灘地，含水層主要由松散的粉細砂夾中粗砂組成，厚度較大，一般在20～50m，是區(qū)內(nèi)第四系潛水主要含水巖組；第二結(jié)構(gòu)層，在基巖承壓水與第四系上更新統(tǒng)潛水之間連續(xù)分布的離石組黃土，其富水性極弱，是本區(qū)潛水與承壓水間主要隔水層，厚度在9.80～120.49m，一般在40～80m；第三結(jié)構(gòu)層，分布于全區(qū)的侏羅系砂泥巖完整巖層，為本區(qū)含煤地層，厚度在15.04～69.10m，總體由北向南，地層厚度增大，且規(guī)律較明顯；第四結(jié)構(gòu)層，分布于全區(qū)的侏羅系砂泥巖層，是透水性較差的孔隙裂隙承壓水含水層，在煤田開采區(qū)，是受開采影響的導水冒裂帶，厚度在15～30m，冒裂帶由于受巖層大裂隙影響，完全導水；第五結(jié)構(gòu)層，概化的3號煤層采空區(qū)，厚度為4.8m。3.1.3概化類邊界處理模型區(qū)第四系潛水含水層（1層）頂部邊界為潛水面；潛水含水層四周的西南邊角概化為第一類定水頭邊界，西南邊界南段榆溪河河段概化為變流量邊界（即第三類變流量邊界，按泉集河處理），其余按零流量邊界處理。模型區(qū)第四系弱透水、弱含水的碎屑巖層（2、3、4、5層）周界及5層底界均概化為零流量邊界。3.1.4評價區(qū)域內(nèi)地下水流數(shù)學模型評價區(qū)地下水流數(shù)學模型如下：3.2預測開采井的溢流量3.2.1分層開采涌水量預測依據(jù)區(qū)內(nèi)水文地質(zhì)特征以及杭來灣煤礦先期開采地段（301盤區(qū)及302盤區(qū)），本次3號煤層開采涌水量預測共設(shè)置2種方案：第1預測方案：3號煤層首采區(qū)（301盤區(qū)及302盤區(qū)）分層開采、平均開采厚度4.8m工況下的礦坑涌水量預測。第2預測方案：3號煤層首采區(qū)（301盤區(qū)及302盤區(qū)）不分層開采、平均開采厚度8.0m工況下的礦坑涌水量預測。3.2.2計算區(qū)有限差分剖分網(wǎng)將各預測方案的大氣降水入滲補給量按面源加入模型，將地下水流模型識別得到的末流場作為方案預測的初始流場，計算區(qū)有限差分剖分網(wǎng)格如圖1所示。依據(jù)所建立的數(shù)值模型，計算得到各預測方案的末流場、末降深場、降深曲線、泉集河流量以及補排均衡量成果如圖2～5、表1、2所示。4煤層開采對礦井總涌水量的影響(1）根據(jù)數(shù)值模型預測結(jié)果，杭來灣礦井采用分層開采煤層時礦坑涌水量為30805m(2）杭來灣煤礦開采煤層對地下水位的影響主要集中在開采區(qū)及其附近，煤層開采導致第四系松散層潛水位降深在0.5～2.5m、侏羅系碎屑巖層承壓水頭降深在3～25m