榆神府地區(qū)煤炭開采對地下水資源的影響

柳寧 趙曉光 解海軍 李瑜

摘 要：通過對覆巖含水性和地下水流場的研究，探討了榆神府地區(qū)煤炭開采對地下水資源造成的影響。通過公式估算法、瞬變電磁法和GMS地下水?dāng)?shù)值模擬法，探討了榆神府地區(qū)煤炭開采對地下水資源量、上覆巖層含水性和地下水流場的影響。研究發(fā)現(xiàn)煤炭開采會引起地下水水位下降，覆巖含水性降低，引起含水滲漏現(xiàn)象，最終改變原始地下水流暢，形成地下水降落漏斗。榆神府地區(qū)在產(chǎn)煤礦礦井水排水量（2014年）約為50×106 m3/a，其中有46.6%的礦井水會作為生產(chǎn)用水回用，隨著在建和規(guī)劃煤礦逐步建成，地下水外排量會增加140×106 m3/a，其中生產(chǎn)回用水量占礦井總排水量的50%以上;上覆巖層含水性在煤炭開采的不同開采階段表現(xiàn)出不一樣的特征，開采前覆巖含水性良好，開采中覆巖含水性逐漸降低，開采結(jié)束后，覆巖含水性會逐漸恢復(fù);GMS模擬結(jié)果顯示，煤炭開采初期地下水水位會以20 m/a的速度快速下降，到2018年已經(jīng)形成地下水降落漏斗，2028年漏斗范圍持續(xù)發(fā)展，到2048年，降落漏斗范圍基本穩(wěn)定，并形成以井田北部為中心的新地下水流動場。

關(guān)鍵詞：地質(zhì)資源與地質(zhì)工程;地下水;瞬變電磁;地下水?dāng)?shù)值模擬;覆巖含水性;

地下水流場

中圖分類號：P 641?文獻標(biāo)志碼：ADOI：10.13800/j.cnki.xakjdxxb.2019.0111文章編號：1672-9315（2019）01-0071-08

Influence of coal mining on groundwater

resources in Yushenfu AreaLIU Ning，ZHAO Xiao?guang，XIE Haijun，LI Yu

（College of Geology and Environment，Xi’an University of Science and Technology，Xi’an 710054，China）Abstract：The influence of coal mining on groundwater in Yushenfu is discussed by studying the aquosity of overburden rock and groundwater flow field.The influence of coal mining on the amount of groundwater resources，the aquosity of overburden rock and the groundwater flow field in Yushenfu is investigated through formula estimation，transient electromagnetic method and GMS groundwater numerical simulation method.It is found that coal mining will cause the decline of groundwater level，reduce the aquosity of overburden rock，cause water leakage，and eventually change the flow of the original ground water to form a groundwater cone of depression.The displacement of mine water in Yushenfu（2014）is about

50×106m3/a，of which 47.6% of the mine water will be used as the production water.With the construction and planning of the coal mine，the external displacement of the groundwater will increase by 140×106m3/a，of which 50% of the total amount of displacement of mine water will be used as production water.The aquosity of overburden rock shows different characteristics in the different stages of coal mining.The aquosity of overburden rock before mining is good，then it will gradually reduce，and after the mining it will be gradually restored.The GMS simulation results show that the groundwater level in the early stage of coal mining will fall sharply with the speed of 20 m/a，and the groundwater cone of depression has been formed by 2018.In 2028，the range of the cone of depression will continue to expand.By 2048，the range of the cone of depression is basically stable，and the new groundwater cone of depression is formed in the north of the mine field.

Key words：geological resources and geological engineering;groundwater;time domain electromagnetic methods;groundwater model system simulation;groundwater flow field

0?引?言

煤炭開采產(chǎn)生較大范圍的巖體結(jié)構(gòu)改變，將會對地下水流場以及覆巖含水性造成影響[1]。位于陜西北部的榆神府礦區(qū)是西北地區(qū)的主要產(chǎn)煤地，煤礦開采時當(dāng)?shù)氐闹饕a(chǎn)業(yè)，是重要的經(jīng)濟地位。該地區(qū)屬于干旱、半干旱大陸性季風(fēng)氣候，水資源短缺，而煤炭開采影響地下水賦存狀態(tài)，進一步加重了當(dāng)?shù)厮Y源匱乏的現(xiàn)象[2-4]。近幾年隨著陜北榆神府煤炭資源的大規(guī)模開發(fā)，所引起的地下水問題也日益凸顯，諸如地下水水位下降、地下水流場發(fā)生變化、地下水水質(zhì)惡化[5-7]。因此，研究采煤對地下水的影響是非常有意義的。國內(nèi)外研究者圍繞煤炭開發(fā)對地下水的影響展開了大量的研究[8-10]，冀瑞君等以窟野河中下游地區(qū)為例，通過現(xiàn)場調(diào)查和水化學(xué)特征分析，總結(jié)了神府礦區(qū)地下水循環(huán)規(guī)律，提出了采煤對地下水的影響機制[11];劉海濤等運用地下水?dāng)?shù)值模擬技術(shù)，深入研究分析了太原西山煤礦采前、采中地下水特征，陳述了地下水流場以及采煤對地下水均衡的影響。結(jié)果表明采煤會造成地下水位顯著下降，地下水動態(tài)平衡受到嚴重干擾[12]。范立民等以陜北地區(qū)的第四系薩拉烏蘇組含水層為例，結(jié)合該地區(qū)含水層的特征及分布，闡述了采煤與地下含水層之間的聯(lián)系，對該地區(qū)采煤過程中地下水保護和預(yù)防措施提出了中肯的建議[13]。張發(fā)旺等結(jié)合干旱區(qū)地下水的特征，針對性的提出“含水層再造”的新理念，闡述了如何協(xié)調(diào)煤層頂板和地板含水層再造與地下水系統(tǒng)保護之間的關(guān)系，使二者共存[14]。董震雨等以杭來灣煤礦為例通過冒裂帶高度、地下水靜儲量、地下水動儲量、采空區(qū)積水的計算和實地采集水樣的檢測，定性及定量分析了煤礦開采對榆溪河流域地下水資源的破壞影響，發(fā)現(xiàn)煤炭開采會影響地下水儲量，但對地下水水質(zhì)影響較小[15]。

從上述研究成果來看，目前大多數(shù)研究都是采用單一的模擬或是實測方法來研究煤炭開采對地下水問題，文中采用模擬

與實測相結(jié)合的方法深入的探討研究采煤對上覆巖層含水性、地下水水量和地下水流場的影響，以期為進一步認識煤炭開采對地下水的影響的基本特征和規(guī)律提供參考。

1?礦井水資源現(xiàn)狀

榆神府地區(qū)現(xiàn)有煤礦477座，其中9個煤礦生產(chǎn)規(guī)模為1.0×106 t/a，0.50×106 ～1.0×106 t/a的煤礦22座;小于0.5×106 t/a的42座;小于0.1×106 t/a的404座。由于整合及停產(chǎn)，煤礦數(shù)量下降了40%，但產(chǎn)能卻不斷增加。榆神府地區(qū)煤礦眾多，排水量很大，生產(chǎn)中的煤礦年排水量為41.7×106 m3/a，其中榆神礦區(qū)礦井水排放量為

27.3×106 m3/a，占該地區(qū)排放量的65%;神府礦區(qū)礦井水排放量為10.7×106 m3/a，占總排放量的26%，榆橫礦區(qū)大部分煤礦已閉礦，外排水量非常小，約為

37.5×106 m3.因此本次研究以礦井水外排量較大的榆神礦區(qū)（錦界煤礦）、神府礦區(qū)（大柳塔煤礦）為重點研究對象，研究煤炭開采對地下水資源的影響。

2?煤炭開采對地下水水量的影響

煤礦排水對地下水動態(tài)影響是第一因素，對地下水水量的計算是研究井下采煤對地下水影響的第一步[16-17]，根據(jù)表1中對榆神府地區(qū)3大煤礦區(qū)主要煤礦礦井水排放量現(xiàn)狀調(diào)查數(shù)據(jù)，根據(jù)下式計算噸煤排水系數(shù)，然后估算煤炭開采年排水量。

α=M/Q式中?α為噸煤排水系數(shù)，m3/t;M為原煤產(chǎn)量，t;Q為年排水量，m3/a.

通過對榆神府地區(qū)3大礦區(qū)進行實地調(diào)查得排水量數(shù)據(jù)，結(jié)合上述公式可計算出噸煤排水系數(shù)，見表1，根據(jù)表1可以估算神府榆地區(qū)3大煤礦區(qū)因采煤而流失的地下水總水量，見表2.

由表2可知，榆神府地區(qū)生產(chǎn)煤礦礦井水現(xiàn)狀排水量約為50×106 m3/a，隨著在建煤礦及規(guī)劃煤礦的陸續(xù)建成投產(chǎn)，礦井排水量會逐步增加，在建煤礦建成投產(chǎn)增加排水量為42.6×106 m3/a，規(guī)劃煤礦建成投產(chǎn)可再增加排水量為101×106 m3/a;礦井疏干水只是部分外排，有很大一部分會作為生產(chǎn)用水回用，查閱榆林市礦井資料顯示，已投產(chǎn)煤礦礦井生產(chǎn)回用水占總排水量的47.6%，在建和規(guī)劃中的礦井建成后預(yù)計礦井生產(chǎn)回用水占總礦排水量的50%以上。

3?煤炭開采對上覆巖層含水性的影響

選用瞬變電磁法對上覆巖層含水性進行研究，該方法屬于電磁感應(yīng)類探測方法[18]，它遵循電磁感應(yīng)原理，其機理就是利用不接地回線或接地線源向地下發(fā)射一次脈沖磁場，在一次脈沖磁場間歇期間利用線圈或接地電極觀測地下介質(zhì)中引起的二次感應(yīng)渦流場[19-20]，從而探測介質(zhì)電阻率的一種方法。

3.1?探測儀器及探測點布設(shè)

瞬變電磁法勘探使用澳大利亞Alpha公司生產(chǎn)的TerraTEM中功率瞬變電磁系統(tǒng)[21]。選取榆神府地區(qū)具有代表性的紅柳林煤礦作為研究對象，布點時采用中心回線法，發(fā)射線圈為10 m×10 m，兩匝。接收線圈為中心探頭，其等效面積為1 000 m2，測線總長1 000 m，每20 m布設(shè)一個測點，共布設(shè)51個測點，布點過程嚴格按照相關(guān)操作規(guī)范進行。

3.2?探測結(jié)果分析

圖1是3-1煤層開采前巖層富水性瞬變電磁探測等視電阻率擬斷面圖，圖中可以看出，3-1煤層所處位置顯示為相對高阻，巖層的成層性較好，中、上部大部分區(qū)域顯示相對低阻，可能是表層基巖風(fēng)化后裂隙相對富水和中部砂巖、粉砂巖層相對富水后電阻率降低造成。

圖2是3-1煤層開采后巖層富水性瞬變電磁探測等視電阻率擬斷面圖，可清晰看出，煤層所在巖層的電性發(fā)生很大變化。其中80 m處附近原3-1煤層上部位置顯示為相對低阻，這與開采前的相對高阻正好相反，推測是由于3-1煤層開采后造成上覆巖層裂隙增大，隔水層發(fā)生破裂，使上部的含水層中的水導(dǎo)入采空區(qū)造成電阻率降低。200～280，440和720 m處附近顯示相對高阻區(qū)域，推測是由于煤層開采后采空區(qū)未充水導(dǎo)致電阻率較高。而560 m附近處和840 m后的區(qū)域都顯示低阻，推測原因同樣是煤層開采后，下部承載力發(fā)生變換，造成上覆巖層的裂隙增大，進而導(dǎo)通了上部的含水層，其中840 m后的大片低阻區(qū)域可能是含水層中的水匯入采空區(qū)后導(dǎo)致。

4?煤炭開采對地下水流場的影響

選用GMS地下水流場模擬軟件對煤炭開采對地下水流場進行模擬研究，GMS（Groundwater Modeling System）是由Brigham Yong University的環(huán)境模型實驗室和美國軍工水道實驗站合作開發(fā)。GMS是一種先進的地下水建模軟件，可以進行復(fù)雜2D，3D地下水動力場流動模擬[22-23]。GMS軟件核心計算模塊為MODFLOW，其數(shù)學(xué)模型的計算采用有限差分法[24]，見下式[25]

式中?h0（x，y，z）為初始水頭，L;Гi為邊界。

4.1?數(shù)值模擬構(gòu)建

以錦界煤礦2011年10月～2013年10月地下水動態(tài)觀測數(shù)據(jù)為依據(jù)建立模型，預(yù)測煤炭開采對地下水流場的影響。該區(qū)主要含水層是第四系松散層潛水，受采煤影響顯著，因此，將該含水層作為本次模擬的目標(biāo)含水層。

4.1.1?基礎(chǔ)數(shù)據(jù)導(dǎo)入

鑒于數(shù)值模擬需要，將錦界礦區(qū)的MAPGIS底圖、含水層頂?shù)装宓臄?shù)據(jù)、區(qū)內(nèi)各觀測井觀測水位標(biāo)高（2012年10月）等數(shù)據(jù)導(dǎo)入GMS軟件進行模擬。

4.1.2?邊界條件設(shè)立[26]

根據(jù)錦界井田水文地質(zhì)單元特征，井田東部以大樹粱南分水嶺為界，西部以大曼梁地表分水嶺為界，南部以青草界溝域為界，北部以爾林兔與宮泊溝之間的地表分水嶺為界。

4.1.3?關(guān)鍵參數(shù)設(shè)定

本次數(shù)值模擬所需的關(guān)鍵水文地質(zhì)參數(shù)是巖土體滲透系數(shù)和目標(biāo)含水層的貯水率。依據(jù)相關(guān)數(shù)據(jù)資料，并參考《陜西省水資源調(diào)查評價》，采用巖土體滲透系數(shù)經(jīng)驗值如下，見表3.

4.2?模型的辨別與校核

4.2.1?模型的辨別

模型的辨別是數(shù)值模擬過程中一個重要過程，目的是使建立的模型能夠精確的模擬出地下將參數(shù)區(qū)劃法與試點法相結(jié)合進行模型的辨別，使模型可以得到連續(xù)、精確的計算結(jié)果[27-28]。選取水文觀測資料比較完整的時間段（2011年10月～2012年10月）來作為本次模型的辨別期，初始流場圖由2011年10月的觀測數(shù)據(jù)輸出，如圖3所示。

經(jīng)過反復(fù)調(diào)整模型參數(shù)，使得辨識期末的地下水是水位標(biāo)高的觀測值和計算值的擬合度達到98.9%，擬合效果非常好，如圖6所示。辨別期末的地下水流場和地下水流場擬合圖，如圖4和圖5所示。

4.2.2?模型的校核

為了后期達到準(zhǔn)確的模擬預(yù)測效果，對辨識后的模型進行校核。以圖4作為校核的初始流場圖，以2013年10月的地下水觀測數(shù)據(jù)作為依據(jù)，進行模型校核，如圖7～圖9所示。

選取2013年10月份的地下水監(jiān)測數(shù)據(jù)（圖7～圖9），作為模型校核的依據(jù)，將辨別后的2012年10月的計算流場作為初始流場，進行模擬。

通過對模型進行校核，校核期末的地下水視為計算值與觀測值的擬合度達到了99.3%，擬合效果非常好，地下水流場的運動規(guī)律基本一致。因此，該模型水文參數(shù)設(shè)定相對合理，能很好的模擬榆神府地區(qū)地下水動態(tài)變化過程。

4.3?榆神府地下水流動態(tài)預(yù)測

在假設(shè)該礦區(qū)開采條件和自然環(huán)境相對穩(wěn)定的前提下，將研究區(qū)錦界礦區(qū)2013年的地下水觀測數(shù)據(jù)，輸入GMS模擬軟件，來預(yù)測該礦區(qū)未來5年、15年及25年（即2018年、2028及2048）的地下水流場動態(tài)變化情況。預(yù)測結(jié)果如圖10～圖12所示。

分析模擬結(jié)果發(fā)現(xiàn)，2013—2048年的地下水流場圖，由于井下采煤，地下水水位持續(xù)下降，水位年降幅超過20 m.2018年時，在礦區(qū)北部，已初步形成地下水降落漏斗;到2028年時，地下水水位繼續(xù)下降，漏斗范圍繼續(xù)擴大;到2048年，地下水水位達到新的動態(tài)平衡，漏斗范圍基本穩(wěn)定，在礦區(qū)北部形成新的地下水流場。

由于井下采動引起圍巖移動，導(dǎo)致含水層滲漏，使得地下水水位下降，形成以采區(qū)為中心的降落漏斗，并且降落漏斗的范圍和深度會隨著采煤活動的進行持續(xù)增大，直到達到相對平衡狀態(tài)。閉礦后，地下水通過自我恢復(fù)能力，漏斗范圍和深度逐漸縮小，直到達到新的動態(tài)平衡，形成新的地下水流場，這與趙琳在研究漏斗區(qū)地下水運移規(guī)律時的結(jié)論相一致[29]。

錦界礦區(qū)的地下水水位會隨著采煤活動的持續(xù)及時間的法杖逐漸降低，嚴重時會出現(xiàn)地下水降落漏斗，因此需要做好預(yù)防措施。

5?結(jié)?論

1）對地下水水資源量的影響。榆神府地區(qū)生產(chǎn)煤礦礦井水現(xiàn)狀排水量約為50×106 m3/a，其中46.6%作為生產(chǎn)水回用，隨著在建煤礦及規(guī)劃煤礦的陸續(xù)建成投產(chǎn)，礦井排水量會逐步增加，增加排水量大約為140×106 m3/a，50%以上的排水用于生產(chǎn)，合理處理外排礦井水可在一定程度上緩解當(dāng)?shù)厮Y源緊缺這一難題;

2）對上覆巖層含水性的影響。開采前，煤系覆巖中地下含水層含水性較好，煤層基本干燥不含水;開采中，地下含水區(qū)因采動影響，發(fā)生含水層滲漏，含水性降低;開采后，地表層和地下松散含水層低阻特性明顯，含水性水平達到采前，但基巖段總體呈現(xiàn)為煤層頂板砂巖層含水性較弱;

3）對地下水流場的影響。模擬結(jié)果表明：在自然條件及開采速度保持穩(wěn)定的情況下，隨著煤炭的開采，研究區(qū)地下水位逐年遞減，平均減幅約為20 m/a.2018年，在礦區(qū)北部已經(jīng)明顯形成地下水降落漏斗，并且降落漏斗隨著煤炭進一步開采會逐漸擴大;2028年，漏斗范圍和深度持續(xù)增大;到2048年，地下水水位和漏斗范圍基本穩(wěn)定，并形成以井田北部為中心的新地下水流場，直到達到新的動態(tài)平衡狀態(tài)。

參考文獻（References）：

[1] 陳佩佩，許延春，王恩志，等.玉舍河下煤炭開采安全性分析[J].煤田地質(zhì)與勘探，2007，35（4）：46-50.CHEN Pei?pei，XU Yan?chun，WANG En?zhi，et al.Analyses of upper limit of coal mining under Yushe river[J].Coal Geology & Exploration，2007，35（4）：46-50.

[2]張茂省，董?英，杜榮軍，等.陜北能源化工基地采煤對地下水資源的影響及對策[J].地學(xué)前緣，2010，17（6）：235-246.ZHANG Mao?sheng，DONG Ying，DU Rong?jun，et al.The strategy and influence of coal mining on the groundwater resouces at energy and chemical base in the North Shaanxi[J].Earth Science Frontiers，2010，17（6）：235-246.

[3]陶?虹，李?成，柴小兵，等.陜西神府煤田環(huán)境地質(zhì)問題及成因[J].地質(zhì)與資源，2010，19（3）：249-251.TAO Hong，LI Cheng，CHAI Xiao?bing et al.The environmental geological problems and causes of the Shenfu Coalfield in Shaanxi Province[J].Geology and Resources，2010，19（3）：249-251.

[4]常金源，李文平，李?濤，等.神南礦區(qū)煤炭開采水資源滲漏失量評價分析[J].煤田地質(zhì)與勘探，2011，39（5）：42-45.CHANG Jin?yuan，LI Wen?ping，LI Tao，et al.Zonation of water resources leakage due to coal mining in Shennan mining area[J].Coal Geology & Exploration，2011，39（5）：42-45.

[5]曾慶銘，施龍青.山東省煤炭開采對水資源的影響分析及對策研究[J].山東科技大學(xué)學(xué)報（自然科學(xué)版），2009，28（2）：42-46.ZENG Qing?min，SHI Long?qing.Analysis of the influence of coal mining on water resources in Shandong Province and its counter measures[J].Shandong University of Science and Technology（Natural Science），2009，28（2）：42-46.

[6]Wu Q，Wang M，Wu X.Investigations of groundwater bursting into coal mine seam floors from fault zones[J].International Journal of Rock Mechanics and Mining Sciences，2004，41（2）：557-571.

[7]SHI Long?qing，Slngh R N.Study of mine water in rush from floor strata though faults[J].Mine Water and the Environment，2000（3）：140-147.

[8]韓寶平，鄭世書.煤礦開采誘發(fā)的水文地質(zhì)效應(yīng)研究[J].中國礦業(yè)大學(xué)學(xué)報，1994，23（3）：70-77.HAN Bao?ping，ZHENG Shi?shu.An approach to the hydrogeologic effect induced by underground coal mining[J].Journal of China University of Mining & Technology，1994，23（3）：70-77.

[9]Dong S，Liu B，Shi X，et al.The spatial distribution and hydrogeological conteols of fluoride in the confined and unconfined groundwater of Tuoketuo County，Hohhot，Inner Mongolia，China[J].Environmental Earth Sciences， 2015， 74（1）：1-11.

[10]馮海波.內(nèi)蒙古呼倫貝爾草原露天煤礦區(qū)地下水系統(tǒng)演化硏究[D].呼和浩特：內(nèi)蒙古大學(xué)，2017.

FENG Hai?bo.Groundwater sysytem evolution of the opencast coalmine area in Hulun Buir prairie Inner Mongolia[D].Huhehaote：Mongolian University，2017.

[11]冀瑞君，彭蘇萍，范立民，等.神府礦區(qū)采煤對地下水循環(huán)的影響——以窟野河中下游流域為例[J].煤炭學(xué)報，2015，40（4）：938-943.

JI Rui?jun，PENG Su?ping，F(xiàn)AN Li?min，et al.Effect of coal exploitation on groundwater circulation in the Shenfu mine area：an example from middle and lower reaches of the Kuye River basin[J].Journal of China Coal Society，2015，40（4）：938-943.

[12]劉海濤.太原西山礦區(qū)煤炭開采對地下水流場影響的數(shù)值模擬[D].太原：太原理工大學(xué)，2005.

LIU Hai?tao.A numerrical simulation on the influence of underground water flow regime caused by coal minging in Taiyuan Xishan coal area[D].Taiyuan：Taiyuan University of Technology，2005.

[13]范立民，王國柱.薩拉烏蘇組地下水及采煤影響與保護[J].采礦技術(shù)，2006，6（3）：422-425.FAN Li?min，WANG Guo?zhu.Influence and protection of underground water and coal mining in Sara Wusu Formation[J].Mining Technology，2006，6（3）：422-425.

[14]張發(fā)旺，周駿業(yè)，申保宏，等.干旱地區(qū)采煤條件下煤層頂板含水層再造與地下水資源保護[M].北京：地質(zhì)出版社，2006.

ZHANG Fa?wang，ZHOU Jun?ye，SHEN Bao?hong，et al.Reconstruction of coal seam roof aquifer and protection of groundwater resources under coal mining conditions in arid areas[M].Beijing：Geological Press，2006.

[15]董震雨，王雙明.采煤對陜北榆溪河流域地下水資源的影響分析——以杭來灣煤礦開采區(qū)為例[J].干旱區(qū)資源與環(huán)境，2017，31（3）：185-190.DONG，Zhen?yu，WANG Shuang?ming.Influence of coal exploitation on groundwater resources in Yuxi river valley of northern Shaanxi[J].Journal of Arid Land Resources and Environment，2017，31（3）：185-190.

[16]Scheibe T，Yabusakl S.Scaling of flow and transport behavior in heterogeneous groundwater system[J].Advances in Water Resources，1998，22（3）：223-238.

[17]王江莉.煤礦開采地下水環(huán)境影響評價研究——以山西古交煤焦集團金之中煤礦為例[D].太原：太原理工大學(xué)，2011.WANG Jiang?li.Study on groundwater environmental impact assessment of coal mining：take Shanxi Gujiao

char group Jinzhizhong Coalmine as an example[J].Taiyuan：Taiyuan University of Technology，2011.

[18]崔?杰，劉?洋，王振榮.錦界煤礦水文地質(zhì)條件分析及探查技術(shù)研究[J].煤炭工程，2009（9）：72-74.CUI Jie，LIU Yang，WANG Zhen?rong.Hydrogeological condition analysis and exploration technology study of Jinjie Mine[J].Coal Engineering，2009（9）：72-74.

[19]Mazzia A，Putti M.Mixed?finite element and finite volume discretization for heavy brine simulations in groundwater[J].Journal of Computational and Applied Mathematics，2002，147（1）： 191-213.

[20]Hertrich M，Yaramanci U.Joint inversion of surface nuclear magnetic resonance and vertical electrical sounding[J].Journal of Applied Geophysics，2002，50（1）：179-191.

[21]解海軍，孟小紅，王信文，等.煤礦積水采空區(qū)瞬變電磁法探測的附加效應(yīng)[J].煤田地質(zhì)與勘探，2009，37（2）：71-74.XIE Hai?jun，MENG Xiao?hong，WANG Xin?wen，et al.The additional effect of TEM detection and its application in coal mined?out area with accumulated water[J].Coalfield Geology and Exploration，2009，37（2）：71-74.

[22]顏文珠.煤礦開采對地下水影響的數(shù)值模擬研究[D].青島：山東科技大學(xué)，2011.YAN Wen?zhu.Study on effect of mining in groundwater with numerical simulation[D].Qingdao：Shandong University of Science and Technology，2011.

[23]王婷婷.基于GMS的某礦區(qū)地下水?dāng)?shù)值模擬研究[D].太原：中國輻射防護研究院，2015.

WANG Ting?ting.Groundwater numerical simulation research on one mining area based on GMS[D].Taiyuan：China Radiation Protection Research Institute，2015.

[24]侯恩科，童仁劍，馮?潔，等.燒變巖富水特征與采動水量損失預(yù)計[J].煤炭學(xué)報，2017，42（1）：175-182.HOU En?ke，TONG Ren?jian，F(xiàn)ENG Jie，Water enrichment characteristics of burnt rock and prediction on water loss cause by coal mining[J].Journal of Coal，2017，42（1）：175-182.

[25]楊東東，李文平，李小琴，等.深部礦井水煤共采水文地質(zhì)數(shù)值分析[J].地質(zhì)論評，2016，61（4）：835-842.YANG Dong?dong，LI Wen?ping，LI Xiao?qin，et al.Hydrogeological numerical analysis of coal mining in deep mine[J].Geological Review，2016，61（4）：835-842.

[26]盧文喜.地下水運動數(shù)值模擬過程中邊界條件問題探討[J].水利學(xué)報，2003（3）：33-36.LU Wen?xi.Discussion on boundary conditions in numerical simulation of groundwater movement[J].Journal of Water Conservancy，2003（3）：33-36.

[27]武?強，徐?華.地下水模擬的可視化設(shè)計環(huán)境[J].計算機工程，2003，29（6）：69-70.WU Qiang，XU Hua.Visual environment in 3D for groundwater modeling[J].Computer Engineering，2003，29（6）：69-70.

[28]田?廣，王曉波.東灘煤礦奧灰放水試驗數(shù)值模擬[J].煤田地質(zhì)與勘探，2014，42（3）：61-64.

TIAN Guang，WANG Xiao?bo.Numerical simulation of dewatering test of Ordovician limestone in Dongtan mine[J].Coal Geology and Exploration，2014，42（3）：61-64.

[29]師本強.陜北淺埋煤層礦區(qū)保水開采影響因素研究[D].西安：西安科技大學(xué)，2012.

SHI Ben?qiang.Study on influencing factors of water conservation mining in

shallow seam mine area in the north Shaanxi Province[D].Xi’an：Xi’an University of Science and Technology，2012.