陳陸望，任星星，張 杰，陳逸飛，鄭 忻

(合肥工業(yè)大學(xué) 資源與環(huán)境工程學(xué)院，安徽 合肥 230009)

隨著淺部煤炭資源的逐漸枯竭，我國(guó)華北型隱伏煤田的煤炭開(kāi)采正快速地向深部進(jìn)行[1]。深部煤層開(kāi)采造成地下空間應(yīng)力的變化和巖層結(jié)構(gòu)的破壞，改變了深部含水層地下水天然的循環(huán)狀態(tài)，破壞了原有的平衡，導(dǎo)致水文地質(zhì)條件復(fù)雜化，煤礦突水事故頻繁發(fā)生[2-3]。因此，許多學(xué)者對(duì)煤田深部地下水循環(huán)與水文地球化學(xué)特征以及水文地球化學(xué)過(guò)程等進(jìn)行了研究。HUANG等[4]分析了深部含水層的補(bǔ)給來(lái)源與補(bǔ)給機(jī)制；YANG等[5]揭示了深部含水層水質(zhì)演化機(jī)制；QIAN等[6]分析了深部含水層水力聯(lián)系機(jī)制和水文地球化學(xué)演化規(guī)律；LIU等[7]分析了深部含水層水質(zhì)演變的水文地球化學(xué)過(guò)程；QIAO等[8]和LIN等[9]探討了深部含水層水文地球化學(xué)的采動(dòng)響應(yīng)；ZHANG等[10]分析了煤礦關(guān)閉后深部含水層水質(zhì)和水文地球化學(xué)過(guò)程；ZHANG等[11]分析了深部含水層水文地球化學(xué)演化的構(gòu)造控制機(jī)制。以上研究成果從一定程度上揭示了煤田深部含水層水循環(huán)與水文地球化學(xué)演化機(jī)制及其影響因素等。

淮北煤田是我國(guó)華北型隱伏煤田的重要組成部分，隨著深部煤炭資源開(kāi)采及區(qū)域地下水資源開(kāi)發(fā)利用，淮北煤田深部含水層地下水化學(xué)組成、演化及含水層間水力聯(lián)系受到密切關(guān)注[12]。文獻(xiàn)[13-15]通過(guò)多元統(tǒng)計(jì)方法、同位素示蹤、水文地球化學(xué)模擬等手段對(duì)淮北煤田深部含水層水力聯(lián)系、突水機(jī)制以及水文地球化學(xué)演化進(jìn)行了較為系統(tǒng)分析；殷曉曦等[16-17]從時(shí)空角度對(duì)淮北煤田境內(nèi)臨渙礦區(qū)采動(dòng)影響下深部含水層水文地球化學(xué)演化的本質(zhì)進(jìn)行了較為深入探討。近年來(lái)，眾多學(xué)者通過(guò)常規(guī)水化學(xué)、多元統(tǒng)計(jì)、非線(xiàn)性分析及水文地球化學(xué)模擬等方法對(duì)典型華北型隱伏煤田深部含水層開(kāi)展了突水水源識(shí)別[18]、地下水化學(xué)演化[19]、水文地球化學(xué)過(guò)程[20-21]及地下水水質(zhì)評(píng)價(jià)[22]等一系列研究。然而，在這些研究中，多數(shù)研究主要利用煤田范圍內(nèi)單一礦井或某一礦區(qū)地下含水層水化學(xué)組分所呈現(xiàn)的時(shí)空分布規(guī)律進(jìn)行探討，缺乏從整個(gè)煤田區(qū)域尺度上對(duì)深部含水層地下水中主要組分的演化及水文地球化學(xué)過(guò)程進(jìn)行系統(tǒng)分析。

淮北煤田太灰含水層埋藏深，與開(kāi)采煤層近，受采動(dòng)影響強(qiáng)烈，并與下伏奧陶系灰?guī)r巖溶裂隙含水層存在不同程度的水力聯(lián)系，直接關(guān)系到深部煤炭資源開(kāi)采的安全。因此，筆者以淮北煤田太灰含水層為研究對(duì)象，運(yùn)用離子比例分析法和主成分分析法對(duì)太灰含水層水化學(xué)常規(guī)組分歷年數(shù)據(jù)進(jìn)行分析，探討采動(dòng)影響下深部太灰含水層水文地球化學(xué)形成作用與時(shí)空演化規(guī)律。并以淮北煤田境內(nèi)臨渙礦區(qū)為例，應(yīng)用PHREEQC軟件進(jìn)行水文地球化學(xué)模擬，定量描述太灰含水層地下水中主要組分的變化，以此驗(yàn)證水文地球化學(xué)形成作用及其空間差異性。研究成果為華北型隱伏煤田中的煤礦水害防治及區(qū)域地下水資源開(kāi)發(fā)與利用奠定理論基礎(chǔ)。

1 研究區(qū)概況

淮北煤田位于安徽省淮北平原北部，在地貌單元上屬華北大平原的一部分。東起郯廬斷裂帶，西至固始斷裂，南接板橋—固鎮(zhèn)斷層，北臨豐沛隆起。淮北煤田四周大的斷裂構(gòu)造控制了該區(qū)地下水的循環(huán)條件，使其基本上形成一個(gè)封閉-半封閉的網(wǎng)格狀水文地質(zhì)單元?；幢泵禾镏胁堪l(fā)育有宿北斷裂，為區(qū)內(nèi)規(guī)模最大的一條東西向斷裂。宿北斷裂延伸使南、北地質(zhì)差異明顯，對(duì)深部地下水化學(xué)與水動(dòng)力起到控制作用。以宿北斷裂為界可將淮北煤田劃分為2個(gè)一級(jí)水文地質(zhì)單元(南部和北部)。南部以南坪斷層、豐渦斷裂為界劃分為宿縣礦區(qū)和臨渙礦區(qū)(屬二級(jí)水文地質(zhì)單元)，北部以蕭西向斜為界劃分為濉蕭礦區(qū)(屬二級(jí)水文地質(zhì)單元)(圖1)[12]。

根據(jù)區(qū)域地層巖性以及含水賦存空間分布，淮北煤田可劃分為新生界松散層孔隙含水層、二疊系煤系砂巖裂隙含水層、石炭系太原組灰?guī)r巖溶裂隙(簡(jiǎn)稱(chēng)“太灰”)含水層及奧陶系灰?guī)r巖溶裂隙含水層。其中，太灰含水層由灰或淺灰色石灰?guī)r、深灰色泥巖、砂質(zhì)泥巖夾數(shù)層薄煤層組成，以石灰?guī)r為主。石灰?guī)r有12～14層，占地層厚度40%～60%，具有細(xì)晶-粗晶結(jié)構(gòu)，普遍充填方解石脈，并見(jiàn)黃鐵礦晶體[12]。

淮北煤田太灰?guī)r溶裂隙發(fā)育不均一，淺部巖溶裂隙發(fā)育，富水性較強(qiáng)，向深部逐漸減弱?；幢泵禾锉辈块l河向斜附近太原組灰?guī)r出露，接受大氣降水的直接補(bǔ)給，區(qū)內(nèi)太灰含水層與奧陶系灰?guī)r巖溶裂隙含水層以及煤系砂巖裂隙含水層地下水通過(guò)導(dǎo)水?dāng)鄬?、巖溶陷落柱以及采動(dòng)裂隙等溝通，相互間存在不同程度的水力聯(lián)系?；幢泵禾锬喜坑捎诒缓袼缮痈采w，太灰含水層在基巖露頭帶與奧陶系灰?guī)r巖溶裂隙含水層、新生界松散層底部孔隙含水層貫通互補(bǔ)，深部太灰含水層與奧陶系灰?guī)r巖溶裂隙含水層和煤系砂巖裂隙含水層同樣通過(guò)導(dǎo)水?dāng)鄬?、巖溶陷落柱以及采動(dòng)裂隙等溝通形成不同程度的水力聯(lián)系[17]。

2 數(shù)據(jù)收集與整理

考慮到測(cè)試過(guò)程中主客觀誤差的存在，首先對(duì)數(shù)據(jù)進(jìn)行陰陽(yáng)離子平衡檢驗(yàn)，檢驗(yàn)公式[23]為

(1)

式中，E為相對(duì)誤差，%；mc，ma分別為陽(yáng)離子與陰離子毫克當(dāng)量濃度，meq/L。

若E值在±5%范圍內(nèi)，則認(rèn)為數(shù)據(jù)可行。經(jīng)過(guò)式(1)檢驗(yàn)，收集的153個(gè)太灰含水層地下水樣與補(bǔ)充測(cè)試的23個(gè)太灰含水層地下水樣均滿(mǎn)足要求。

3 結(jié)果與分析

3.1 水化學(xué)類(lèi)型分析

將太灰含水層地下水樣常規(guī)水化學(xué)數(shù)據(jù)按不同礦區(qū)分為5個(gè)時(shí)間段分別繪制Piper三線(xiàn)圖，如圖2所示。

圖2 淮北煤田各礦區(qū)太灰含水層地下水樣水化學(xué)Piper三線(xiàn)圖

3.2 水文地球化學(xué)形成機(jī)制

3.2.1離子比例分析

離子比例分析可用來(lái)探究地下水中主要化學(xué)組分的來(lái)源及可能存在的水文地球化學(xué)過(guò)程，如圖3所示。

圖3 淮北煤田各礦區(qū)內(nèi)煤礦太灰含水層地下水樣離子比例關(guān)系

由圖3(a)可看出，當(dāng)?shù)叵滤蠳a+和Cl-僅來(lái)自巖鹽溶濾溶解時(shí)，ρ(Cl-)/ρ(Na+)應(yīng)為1。除宿縣礦區(qū)內(nèi)的祁南礦與臨渙礦區(qū)內(nèi)的五溝礦、任樓礦太灰含水層地下水樣點(diǎn)位于ρ(Cl-)軸與ρ(Na+)軸組成坐標(biāo)系的1∶1線(xiàn)上部，其他大部分水樣點(diǎn)在1∶1線(xiàn)下部，可知巖鹽溶濾溶解不是太灰含水層地下水中的Na+的惟一來(lái)源，富余的Na+可能來(lái)源于陽(yáng)離子交替吸附。

3.2.2主成分分析

將淮北煤田太灰含水層地下水樣歷年常規(guī)水化學(xué)數(shù)據(jù)利用SPSS軟件進(jìn)行標(biāo)準(zhǔn)化，得到各離子之間的相關(guān)系數(shù)矩陣，然后計(jì)算特征值和特征向量，最后計(jì)算主成分貢獻(xiàn)率和累積貢獻(xiàn)率。并采用方差最大旋轉(zhuǎn)法對(duì)主成分軸進(jìn)行旋轉(zhuǎn)，使每個(gè)主成分具有最高因子荷載的變量數(shù)最少，從而簡(jiǎn)化對(duì)主成分的解釋?zhuān)逦亟沂舅瘜W(xué)數(shù)據(jù)的內(nèi)在信息[24-25]。

圖4 淮北煤田太灰含水層水化學(xué)主成分因子荷載

(2)

(3)

(4)

(5)

(6)

3.3 水文地球化學(xué)形成時(shí)空演化規(guī)律分析

3.3.1隨時(shí)間的演化規(guī)律

為了揭示采動(dòng)后淮北煤田太灰含水層地下水主要水文地球化學(xué)形成作用隨時(shí)間演化規(guī)律，按時(shí)間梯度繪制了主成分1與主成分2荷載得分散點(diǎn)圖(圖5)。

由圖5可知，宿縣礦區(qū)太灰含水層地下水樣點(diǎn)在開(kāi)采初期主要位于第1象限，主成分1荷載得分與主成分2荷載得分均較大，太灰含水層地下水以溶濾溶解并伴隨不同程度的陽(yáng)離子交替吸附作用。隨時(shí)間的延續(xù)，宿縣礦區(qū)太灰含水層地下水樣點(diǎn)總體上向左移動(dòng)，主成分1荷載得分減小，主成分2荷載得分增加，受pH值、CO2分壓等影響，太灰含水層方解石和白云石溶濾溶解作用受限，后期以陽(yáng)離子交替吸附作用為主。

臨渙礦區(qū)早期太灰含水層地下水樣點(diǎn)主要位于第1象限，且主成分1荷載得分大于主成分2荷載得分，以溶濾溶解作用為主。隨時(shí)間的延續(xù)，水樣點(diǎn)逐漸向主成分1荷載得分變小的方向移動(dòng)，后期水樣主成分2荷載得分大于主成分1荷載得分，以陽(yáng)離子交替吸附作用為主。

濉蕭礦區(qū)早期太灰含水層地下水樣點(diǎn)位于第3象限，主成分1荷載得分與主成分2荷載得分均較小，水文地球化學(xué)形成作用不明顯。隨著時(shí)間的延續(xù)，濉蕭礦區(qū)太灰含水層地下水樣點(diǎn)逐漸向主成分1荷載得分變大方向移動(dòng)，由于基巖埋藏淺，采動(dòng)后使地下水?dāng)y帶充足的CO2與O2，溶濾溶解作用逐漸增強(qiáng)，而陽(yáng)離子交替吸附作用不明顯。

3.3.2隨空間的演化規(guī)律

繪制淮北煤田各礦區(qū)太灰含水層地下水樣點(diǎn)主成分1荷載得分(SF1)、主成分2荷載得分(SF2)等值線(xiàn)與地質(zhì)背景關(guān)系如圖6～8所示。

由圖6可知，宿縣礦區(qū)SF1值中部區(qū)大，桃園礦附近SF1出現(xiàn)最大值；SF1值西部區(qū)小，鄒莊礦附近SF1值出現(xiàn)最小值。宿縣礦區(qū)SF2值西部區(qū)大，鄒莊礦附近SF2出現(xiàn)最大值；SF2值中部區(qū)小，桃園礦附近SF2出現(xiàn)最小值。桃園礦位于宿南向斜附近，基巖露頭帶太灰?guī)r溶裂隙發(fā)育，含水豐富且水動(dòng)力條件良好，使碳酸鹽、硫酸鹽溶濾溶解作用增強(qiáng)，所以SF1值大。鄒莊礦位于宿南背斜西南翼，太灰含水層埋深逐漸加大并處于相對(duì)封閉的環(huán)境中，陽(yáng)離子交替吸附作用增強(qiáng)，所以SF2值大；朱仙莊礦和蘆嶺礦位于宿東向斜內(nèi)，西坡寺斷層、宿北斷裂、東三鋪斷層及固鎮(zhèn)斷裂使得含水層相對(duì)封閉，水流滯緩，水文地球化學(xué)形成作用較為穩(wěn)定，所以SF1值和SF2值均小。

由圖7可知，臨渙礦區(qū)SF1值西部區(qū)大，青東礦附近出現(xiàn)最大值；SF1值東部區(qū)小，任樓礦附近出現(xiàn)最小值。臨渙礦區(qū)SF2值東部區(qū)大，任樓礦附近出現(xiàn)最大值；SF2值西部區(qū)小，袁店二礦附近出現(xiàn)最小值。青東礦位于臨渙向斜附近，向斜構(gòu)造有利于大氣降水補(bǔ)給，加快了碳酸鹽、硫酸鹽溶濾溶解和黃鐵礦氧化，所以SF1值大；任樓礦位于童亭背斜附近東南翼，太灰含水層埋深逐漸加大，陽(yáng)離子交替吸附作用增強(qiáng)，所以SF2值大。五溝礦、楊柳礦、袁店二礦、孫疃礦被大劉家斷層、楊柳斷層、界溝斷層所圍限，太灰含水層較為封閉，水文地球化學(xué)形成作用穩(wěn)定，所以SF1值和SF2值均小。

由圖8可知，濉蕭礦區(qū)SF1值與SF2值均在西部區(qū)大，劉橋二礦附近SF1與SF2均出現(xiàn)最大值；SF1值與SF2值均是東部區(qū)小，石臺(tái)礦附近SF1出現(xiàn)最小值，雙龍礦附近SF2出現(xiàn)最小值。濉蕭礦區(qū)SF1值與SF2值總體表現(xiàn)出西高東低，溶濾溶解和陽(yáng)離子交替吸附作用由東向西逐漸增強(qiáng)，與地下水徑流路徑一致，由此可知濉蕭礦區(qū)太灰含水層水文地球化學(xué)形成作用受構(gòu)造影響較小，主要與地下水徑流有關(guān)。

圖8 濉蕭礦區(qū)太灰含水層地下水樣主成分荷載得分等值線(xiàn)

3.4 反向水文地球化學(xué)模擬

臨渙礦區(qū)位于淮北煤田的中南部，與宿縣礦區(qū)、濉蕭礦區(qū)相鄰且均屬于二級(jí)水文地質(zhì)單元。筆者選取臨渙礦區(qū)為例，應(yīng)用PHREEQC軟件開(kāi)展反向水文地球化學(xué)模擬，定量描述太灰含水層地下水中主要組分的變化，用以驗(yàn)證上述水文地球化學(xué)形成作用及其演化規(guī)律[25]。

PHREEQC軟件模擬，選擇合適模擬路徑和“可能礦物相”是模型建立的關(guān)鍵，從而確定主要礦物飽和指數(shù)(SI)和模擬路徑上各礦物相的轉(zhuǎn)化量，以此分析水文地球化學(xué)形成作用。沿地下水徑流方向，隨著水文地球化學(xué)形成作用的不斷進(jìn)行，地下水中TDS值呈增加趨勢(shì)[22]。因此利用TDS空間分布特征和區(qū)域水文地質(zhì)條件選定模擬路徑，每條模擬路徑由2個(gè)水樣點(diǎn)構(gòu)成[26]，最終確定了5條模擬路徑。繪制臨渙礦區(qū)太灰含水層地下水TDS等值線(xiàn)及據(jù)此確定的5條模擬路徑如圖9所示。

圖9 臨渙礦區(qū)太灰含水層地下水TDS等值線(xiàn)與模擬路徑

根據(jù)臨渙礦區(qū)太灰含水層巖性與水文地球化學(xué)形成機(jī)制，確定“可能礦物相”為方解石、白云石、石膏、巖鹽、NaX、CaX2、H2S、O2、CO2。其中溶濾溶解模型加入方解石、白云石、石膏、巖鹽；碳酸鹽的溶濾溶解加入CO2，黃鐵礦氧化加入H2S和O2，陽(yáng)離子交替吸附加入NaX和CaX2[2]。礦物的SI表示礦物相對(duì)于地下水的飽和狀態(tài)。反向水文地球化學(xué)模擬出各礦井太灰含水層地下水樣的主要礦物SI見(jiàn)表1。由表1中可知，臨渙礦區(qū)各礦井石膏和巖鹽的SI均小于0，表明其在地下水中始終處于溶解狀態(tài)；除海孜礦之外，方解石和白云石的SI均大于0，表明有沉淀析出的趨勢(shì)。

表1 臨渙礦區(qū)各礦井太灰含水層地下水礦物飽和指數(shù)(SI)

表2 臨渙礦區(qū)太灰含水層不同路徑反向水文地球化學(xué)模擬結(jié)果

4 結(jié) 論

(1)宿縣礦區(qū)太灰含水層地下水化學(xué)類(lèi)型以Cl-Na,HCO3-Na,Cl·SO4-Na·Ca型為主，臨渙礦區(qū)太灰含水層地下水化學(xué)類(lèi)型主要有Cl-Na,Cl·SO4-Ca·Na,SO4-Ca·Mg型，濉蕭礦區(qū)太灰含水層地下水化學(xué)類(lèi)型由前期的HCO3-Na,HCO3·Cl-Na型向SO4·Cl-Ca·Na,SO4-Ca·Mg型轉(zhuǎn)變。

(2)通過(guò)離子比例分析和主成分分析揭示了淮北煤田太灰含水層主要水文地球化學(xué)形成作用類(lèi)型為溶濾溶解和陽(yáng)離子交替吸附作用。隨時(shí)間的延續(xù)，采動(dòng)后宿縣礦區(qū)與臨渙礦區(qū)太灰含水層溶濾溶解作用逐漸減弱，陽(yáng)離子交替吸附作用逐漸增強(qiáng)；濉蕭礦區(qū)太灰含水層埋藏淺，地下水CO2含量充足，溶濾溶解作用逐漸增強(qiáng)。

(3)宿縣礦區(qū)與臨渙礦區(qū)中部太灰含水層受向斜構(gòu)造影響SF1值較大，西部太灰含水層受背斜構(gòu)造影響SF2值較大，東部太灰含水層被斷層圍限，含水層封閉，SF1,SF2值均較小。濉蕭礦區(qū)太灰含水層受構(gòu)造影響不明顯，主要與地下水徑流有關(guān)，地下水由東向西徑流，SF1和SF2逐漸增大。

(4)應(yīng)用PHREEQC軟件對(duì)臨渙礦區(qū)開(kāi)展了反向水文地球化學(xué)模擬，模擬結(jié)果表明沿徑流路徑太灰含水層主要發(fā)生白云石、石膏和巖鹽的溶濾溶解，方解石沉淀，黃鐵礦氧化與陽(yáng)離子交替吸附等，模擬結(jié)果定量驗(yàn)證了離子比例分析和主成分分析定性揭示的水文地球化學(xué)形成作用。