劉東海，鄧念東，姚 婷

(1.西安科技大學(xué) 地質(zhì)與環(huán)境學(xué)院，陜西 西安 710054； 2.重慶市地質(zhì)礦產(chǎn)勘查開發(fā)局107地質(zhì)隊(duì)，重慶 401120)

潞安礦區(qū)是我國(guó)重要的煤炭生產(chǎn)基地，位于山西省上黨盆地沁水煤田東南部，礦區(qū)內(nèi)煤炭資源儲(chǔ)量豐富，分布有五陽(yáng)、常村、王莊、漳村、司馬、屯留等大型煤礦[1]。長(zhǎng)期的煤炭開采活動(dòng)導(dǎo)致地表沉陷面積達(dá)64.23 km2，造成了土地資源、建筑物、水體等的破壞，截至2013年，潞安礦區(qū)因采空區(qū)塌陷造成的經(jīng)濟(jì)損失高達(dá)24 854.33 萬(wàn)元[2]，開采沉陷問題嚴(yán)重制約了潞安礦區(qū)可持續(xù)發(fā)展的進(jìn)程。研究表明，煤炭開采沉陷的產(chǎn)生和發(fā)展與地質(zhì)、采礦條件有著密切的聯(lián)系[3]，不同的地質(zhì)、采礦條件下，煤礦開采引起的巖層和地表移動(dòng)存在較大的差異，地表沉陷破壞類型也各不相同?；谏鲜鲈蛩碌牡乇沓料荩嚓P(guān)科技工作者已做過許多深入的研究，并取得了較為成熟的研究成果。崔鋒等根據(jù)峰峰礦區(qū)現(xiàn)有的地質(zhì)條件，采用數(shù)值模擬手段研究了斷層活化造成地表異常沉陷的影響因素，并建立了斷層活化力學(xué)模型[4]；于秋鴿等通過引入地表下沉偏態(tài)程度，分析了各種因素對(duì)地表偏態(tài)下沉的影響，并利用數(shù)值模擬手段分析了松散層厚度對(duì)地表偏態(tài)下沉的影響[5]；高玉兵等采用現(xiàn)場(chǎng)實(shí)測(cè)和數(shù)值模擬等方法對(duì)不同基巖厚度、不同巖性地層進(jìn)行模擬分析，研究了不同采動(dòng)情況及不同地質(zhì)條件下的地層變形特征[6]；宋世杰、趙曉光等采用灰色關(guān)聯(lián)分析、主成分分析等方法篩選出影響陜北榆神府礦區(qū)開采沉陷的關(guān)鍵地礦因素[7-8]；羅海燕等利用正交試驗(yàn)系統(tǒng)分析了條帶工作面采寬、留寬和煤體強(qiáng)度對(duì)采場(chǎng)沉陷的綜合影響規(guī)律[9]；許國(guó)勝等通過多元線性回歸分析，指出影響地表最大下沉速度滯后角正切值的主要因素是煤層上覆巖層的巖性強(qiáng)弱[10]；趙兵朝等通過數(shù)值模擬及正交試驗(yàn)分析得出厚松散層淺埋深煤層開采條件下地表最大下沉值的影響因素及其影響程度的排序[11-12]。

有關(guān)潞安礦區(qū)煤炭開采沉陷影響因素的研究較少。筆者從地質(zhì)賦存條件和開采設(shè)計(jì)工藝2個(gè)方面初選出9個(gè)影響因子，并利用收集的潞安礦區(qū)10個(gè)工作面的地質(zhì)采礦資料及觀測(cè)資料對(duì)其進(jìn)行灰色關(guān)聯(lián)度分析和回歸預(yù)計(jì)，進(jìn)而確定出潞安礦區(qū)開采沉陷的主要控制因素，為潞安礦區(qū)的開采沉陷預(yù)計(jì)和防治提供理論依據(jù)，研究分析對(duì)礦區(qū)生態(tài)環(huán)境保護(hù)和推動(dòng)綠色礦山建設(shè)具有重要的現(xiàn)實(shí)意義。

1 潞安礦區(qū)開采沉陷機(jī)理

潞安礦區(qū)地處長(zhǎng)治盆地西部，廣為厚濕陷性黃土覆蓋。主采二疊系山西組3#煤層，煤層傾角為3°～8°，屬近水平煤層，平均厚度為6.42 m，采高為5.61～7.23 m，煤層開采厚度大；煤層埋深一般在300 m左右?；鶐r以砂巖和砂質(zhì)泥巖為主，覆巖類型屬中硬，松散層厚度為0～150 m。采用綜采放頂煤開采、全部陷落法控制頂板[13]。在煤炭開采過程中，隨著工作面的不斷推進(jìn)，上覆巖層移動(dòng)經(jīng)歷破裂→平衡→再破裂→再平衡的動(dòng)態(tài)發(fā)展過程。隨著開采范圍的不斷擴(kuò)大，巖層的移動(dòng)變形范圍逐漸波及到地表，地表點(diǎn)經(jīng)歷開始移動(dòng)→劇烈移動(dòng)→移動(dòng)停止的全過程，工作面煤炭開采停止后，經(jīng)過長(zhǎng)期緩慢的移動(dòng)變形，地表會(huì)呈現(xiàn)為下沉盆地的形態(tài)[14]。當(dāng)深厚比小于30時(shí)，地表沉陷為非連續(xù)變形；當(dāng)深厚比大于30時(shí)，地表沉陷為連續(xù)變形[15]。潞安礦區(qū)的深厚比普遍在30以上，其地表沉陷一般為連續(xù)變形。潞安礦區(qū)煤層上覆基巖在開采后易發(fā)生冒落，頂板隨采隨垮，從而形成3種不同的采動(dòng)影響帶，即垮落帶、斷裂帶和彎曲帶。工作面采用綜放開采方式一次采全高，采高及基巖厚度均較大，由于工作面推進(jìn)速度快，開采強(qiáng)度大，頂板卸壓速率快，導(dǎo)致基巖內(nèi)部產(chǎn)生近似連續(xù)變形，離層、裂隙未得到充分發(fā)育便成塊垮落，使得上覆巖層破壞異常劇烈，致使活躍期內(nèi)地表移動(dòng)變形更為劇烈且集中，從而加劇了地表的沉陷變形；由于松散層抗彎能力很差，相當(dāng)于荷載上覆于基巖，加劇了巖層的移動(dòng)破壞。與一般地質(zhì)、采礦條件相比，上述因素導(dǎo)致地表下沉量增大，地表沉陷變形異?；钴S；對(duì)采動(dòng)影響范圍廣，其破壞性大，對(duì)地面建筑物的影響明顯。

2 煤炭開采沉陷影響因素的灰色關(guān)聯(lián)分析

2.1 灰色關(guān)聯(lián)分析原理

灰色關(guān)聯(lián)分析是通過灰色關(guān)聯(lián)度來分析和確定系統(tǒng)因素間的影響程度或因素對(duì)系統(tǒng)主行為的貢獻(xiàn)測(cè)度的一種方法，其計(jì)算步驟主要為：①確定母因素與子因素。將反映系統(tǒng)行為特征的某個(gè)因素作為母因素，將影響系統(tǒng)行為的因素作為子因素。②構(gòu)建原始數(shù)據(jù)矩陣。通過搜集多組在不同時(shí)間不同空間內(nèi)觀測(cè)得到的母因素和子因素的相關(guān)數(shù)據(jù)，建立原始數(shù)據(jù)矩陣。③數(shù)據(jù)無(wú)量綱化處理。常用的方法有初值化和均值化。④絕對(duì)差值計(jì)算。計(jì)算各個(gè)子因素與母因素之間的絕對(duì)差值，并構(gòu)建絕對(duì)差值矩陣。⑤計(jì)算關(guān)聯(lián)系數(shù)矩陣。⑥計(jì)算關(guān)聯(lián)度。計(jì)算關(guān)聯(lián)系數(shù)各序列的平均值，即關(guān)聯(lián)度，從而對(duì)子因素和母因素之間的相關(guān)程度進(jìn)行總體比較。

2.2 影響因素選取

開采沉陷的影響因素較多，其中地質(zhì)因素和采礦因素起主要影響作用。地質(zhì)因素對(duì)開采沉陷起主要控制作用，影響開采沉陷的地質(zhì)因素主要有覆巖結(jié)構(gòu)、覆巖的物理力學(xué)性質(zhì)、煤層傾角、松散層性質(zhì)、構(gòu)造環(huán)境等。采礦因素對(duì)開采沉陷起誘發(fā)作用，采空區(qū)的形成是地表沉陷的決定因素，煤炭資源的賦存狀態(tài)和煤炭開采方式直接影響著地表沉陷的類型和特征；采礦因素主要有采高、采深、開采尺寸、開采工藝和頂板管理方法等[16]。筆者依據(jù)前人對(duì)潞安礦區(qū)煤炭開采沉陷機(jī)理的分析[17-18]，充分考慮潞安礦區(qū)的自身特點(diǎn)，結(jié)合實(shí)際的地質(zhì)、采礦條件，遵循科學(xué)性、代表性和易操作性的原則，初選出 9個(gè)影響因素，如表1所示。

表1 開采沉陷影響因素

2.3 原始數(shù)據(jù)采集

筆者收集了潞安礦區(qū)10個(gè)典型工作面的地質(zhì)和采礦條件等資料、觀測(cè)數(shù)據(jù)[1]，將其整理分析后，以表征開采沉陷強(qiáng)度的地表下沉系數(shù)作為母因素X0，將選取出的9個(gè)影響因素作為子因素，得到灰色關(guān)聯(lián)分析的原始數(shù)據(jù)矩陣，如表2所示。

表2 潞安礦區(qū)煤炭開采沉陷影響因素原始數(shù)據(jù)

2.4 灰色關(guān)聯(lián)度計(jì)算結(jié)果

經(jīng)過計(jì)算，可得到各影響因素與地表下沉系數(shù)的灰色關(guān)聯(lián)度。關(guān)聯(lián)度用G(0,i)表示，其中i代表第i個(gè)影響因素。選取的9個(gè)子因素與母因素的灰色關(guān)聯(lián)度大小如表3所示。

表3 灰色關(guān)聯(lián)度計(jì)算結(jié)果

由表3可以看出，將灰色關(guān)聯(lián)度按照由大到小的順序排列，得到灰色關(guān)聯(lián)度排序如下：G(0,5)>G(0,7)>G(0,4)>G(0,9)>G(0,6)>G(0,8)>G(0,3)>G(0,1)>G(0,2)。即各因素對(duì)煤炭開采沉陷影響程度的大小依次為：采高M(jìn)>工作面傾向長(zhǎng)度D1>覆巖綜合硬度Q>擾動(dòng)系數(shù)n>采深H0>工作面走向長(zhǎng)度D3>基巖厚度H巖>煤層傾角α>松散層厚度H土。

由灰色關(guān)聯(lián)度計(jì)算結(jié)果可以看出，采高M(jìn)、工作面傾向長(zhǎng)度D1和覆巖綜合硬度Q的灰色關(guān)聯(lián)度均在0.800以上，其中采高M(jìn)與地表下沉系數(shù)的關(guān)聯(lián)度最大為0.850，工作面傾向長(zhǎng)度D1和覆巖綜合硬度Q次之，說明上覆巖層的力學(xué)性質(zhì)和煤炭開采強(qiáng)度在很大程度上決定著煤炭開采沉陷損害的程度。上覆巖層力學(xué)性質(zhì)體現(xiàn)了礦山地質(zhì)環(huán)境抗損害的能力，而煤炭開采強(qiáng)度則體現(xiàn)了人類活動(dòng)對(duì)礦山地質(zhì)環(huán)境的破壞能力，二者之間的抗衡則通過地表沉陷程度得以表現(xiàn)[19]。擾動(dòng)系數(shù)n、采深H0、工作面走向長(zhǎng)度D3及基巖厚度H巖與地表下沉系數(shù)的關(guān)聯(lián)度介于0.700～0.800，煤層傾角α與地表下沉系數(shù)的關(guān)聯(lián)度超過0.650，表明開采沉陷損害的程度與覆巖的結(jié)構(gòu)、煤層自身傾角大小及工作面的尺寸有著密切的聯(lián)系。松散層厚度H土與地表下沉系數(shù)的關(guān)聯(lián)度最小為0.531。

3 煤炭開采沉陷影響因素的逐步回歸分析

3.1 回歸模型建立

在自變量較多的情況下可通過逐步回歸分析確定主要影響因素。逐步回歸分析的基本思路是將變量逐個(gè)引入回歸方程，每引入一個(gè)新的變量，便對(duì)引入的所有變量進(jìn)行檢驗(yàn)，當(dāng)已引入的變量由于新變量的引入而變得不再顯著時(shí)，便將其剔除，以保證引入的變量都是顯著的，如此反復(fù)進(jìn)行，直到不再引入新變量為止。

根據(jù)以上灰色關(guān)聯(lián)度分析結(jié)果，將采高、工作面傾向長(zhǎng)度、覆巖綜合硬度、擾動(dòng)系數(shù)、采深、工作面走向長(zhǎng)度、基巖厚度、煤層傾角、松散層厚度9個(gè)影響因素作為自變量，地表最大下沉系數(shù)作為因變量，選取表2中第1～8組數(shù)據(jù)進(jìn)行逐步回歸分析。選取不同的閾值時(shí)，引入的自變量及回歸方程的顯著性檢驗(yàn)結(jié)果如表4所示。

表4 逐步回歸分析結(jié)果

由表4可知，當(dāng)閾值為7.2時(shí)，引入方程的變量只有采高和覆巖綜合硬度；將閾值減小至1.7時(shí)，同時(shí)引入了基巖厚度和松散層厚度；當(dāng)閾值為1.0時(shí)，擾動(dòng)系數(shù)、工作面傾向長(zhǎng)度依次被引入，此時(shí)回歸方程的復(fù)相關(guān)系數(shù)最大為1.000，標(biāo)準(zhǔn)差最小為0.004，回歸方程最顯著，對(duì)應(yīng)的地表最大下沉系數(shù)與煤炭開采沉陷影響因素之間的線性回歸方程為：

Y=-0.537+0.036X5+0.219X4-0.001X3+0.001X2-0.097X9+0.002X7

由偏回歸系數(shù)可知，采高、覆巖綜合硬度、松散層厚度，以及工作面傾向長(zhǎng)度與地表下沉系數(shù)呈正相關(guān)關(guān)系；基巖厚度、擾動(dòng)系數(shù)與地表下沉系數(shù)呈負(fù)相關(guān)關(guān)系。由偏回歸系數(shù)可知各變量對(duì)地表下沉系數(shù)的貢獻(xiàn)率大小依次為：覆巖綜合硬度>擾動(dòng)系數(shù)>采高>工作面傾向長(zhǎng)度>基巖厚度、松散層厚度。

3.2 模型驗(yàn)證

運(yùn)用回歸模型，對(duì)潞安礦區(qū)10個(gè)工作面煤炭開采的下沉系數(shù)進(jìn)行預(yù)計(jì)，并與實(shí)測(cè)值進(jìn)行對(duì)比，實(shí)測(cè)值和回歸預(yù)計(jì)值對(duì)比如表5所示。

表5 地表下沉系數(shù)實(shí)測(cè)值和回歸預(yù)計(jì)值對(duì)比

由表5可知，最小偏差為0.003，最大偏差為0.057，平均偏差為0.032。說明該預(yù)計(jì)模型對(duì)潞安礦區(qū)典型工作面煤炭開采下沉系數(shù)預(yù)計(jì)效果較好，精度能夠滿足實(shí)際需求。

3.3 主要控制因子篩選

根據(jù)灰色關(guān)聯(lián)度計(jì)算與逐步回歸分析結(jié)果，將各因素對(duì)煤炭開采沉陷的影響程度進(jìn)行對(duì)比，最終篩選出主要控制因子。各影響因素貢獻(xiàn)率大小排序如表6所示。

表6 影響因素貢獻(xiàn)率大小排序

由表6可知，灰色關(guān)聯(lián)分析和逐步回歸分析的結(jié)果均表明：采高、覆巖綜合硬度、工作面傾向長(zhǎng)度、擾動(dòng)系數(shù)對(duì)煤炭開采沉陷的影響程度較大，因此，最終選擇覆巖綜合硬度、采高、工作面傾向長(zhǎng)度、擾動(dòng)系數(shù)作為潞安礦區(qū)煤炭開采沉陷的主要控制因素。

4 結(jié)論

1)從地質(zhì)賦存條件和開采設(shè)計(jì)工藝兩個(gè)方面初選出9個(gè)影響因素，并對(duì)其進(jìn)行灰色關(guān)聯(lián)度分析，得到各因素對(duì)煤炭開采沉陷影響程度的排序：采高>工作面傾向長(zhǎng)度>覆巖綜合硬度>擾動(dòng)系數(shù)>采深>工作面走向長(zhǎng)度>基巖厚度>煤層傾角>松散層厚度。

2)以地表最大下沉系數(shù)為因變量，各影響因素為自變量，建立煤炭開采沉陷回歸模型，得到了地表最大下沉系數(shù)與煤炭開采沉陷影響因素之間的線性回歸方程，其復(fù)相關(guān)系數(shù)最大為1.000，標(biāo)準(zhǔn)差最小為0.004，顯著性較高，且據(jù)其得到的回歸地表最大下沉系數(shù)與實(shí)測(cè)值的平均偏差為0.032，能較準(zhǔn)確地反映煤炭開采沉陷與各影響因子之間的線性關(guān)系。

3)通過灰色關(guān)聯(lián)度計(jì)算和逐步回歸分析，最終確定潞安礦區(qū)采煤沉陷的主要控制因素為覆巖綜合硬度、采高、工作面傾向長(zhǎng)度及擾動(dòng)系數(shù)。