湯靜雅 楊志強(qiáng) 蘆家欣

摘要：村莊下壓煤多采取房屋搬遷方式來實(shí)現(xiàn)安全開采。煤礦村落選址區(qū)的地理適宜度受諸多自然和人文社會(huì)因素的影響，目前缺乏普遍適用的評(píng)價(jià)方法。文中針對(duì)黃土礦區(qū)復(fù)雜的地理環(huán)境和采煤沉陷特征，選取黃土溝壑區(qū)地面坡度、坡向、地表水和地下水源等自然地理指標(biāo)以及交通便捷性、耕作半徑等社會(huì)經(jīng)濟(jì)指標(biāo)作為主要影響因子，將各指標(biāo)劃分為3個(gè)等級(jí)，確定相應(yīng)的評(píng)價(jià)標(biāo)準(zhǔn)，采用專家打分及層次分析法計(jì)算各指標(biāo)相應(yīng)權(quán)重，構(gòu)建黃土高原礦區(qū)村落遷址的地理適宜性評(píng)價(jià)模型。為驗(yàn)證模型的適用性，以陜西某礦塌陷區(qū)的村落遷址為例，基于ArcGIS的疊加分析功能，提取各指標(biāo)等級(jí)對(duì)應(yīng)區(qū)域，并按評(píng)價(jià)模型的對(duì)應(yīng)指標(biāo)權(quán)重進(jìn)行疊加，得到礦區(qū)地理適宜度的綜合評(píng)價(jià)分布圖;對(duì)采煤塌陷區(qū)內(nèi)的建筑面積進(jìn)行統(tǒng)計(jì)，結(jié)合評(píng)價(jià)圖分析并確定村落搬遷的最佳選址位置。結(jié)果表明，該模型在黃土高原煤礦塌陷區(qū)村落遷址中具有較好的適用性，可為煤礦塌陷區(qū)村落遷址規(guī)劃提供技術(shù)支持。

關(guān)鍵詞：黃土高原;采煤塌陷區(qū);村落遷址;地理適宜度

中圖分類號(hào)：TD 325文獻(xiàn)標(biāo)志碼：A

DOI：10.13800/j.cnki.xakjdxxb.2019.0221文章編號(hào)：1672-9315（2019）02-0334-07

0引言

西部黃土高原及其過渡地帶是我國(guó)主要的煤炭生產(chǎn)基地之一，地表溝壑縱橫，地貌復(fù)雜多變，水資源匱乏，人居生態(tài)環(huán)境本身較為脆弱，在大規(guī)模的地下采煤導(dǎo)致的沉陷變形影響下，地表建（構(gòu)）筑物不可避免地會(huì)產(chǎn)生破壞，嚴(yán)重危及到采煤沉陷區(qū)內(nèi)建筑物的安全和居民的正常生活[1-4]。目前，在村（鎮(zhèn)）居民區(qū)下實(shí)現(xiàn)安全采煤的途徑主要包括充填開采、條帶開采、建筑加固等方式[5-6]。由于黃土高原自然村落多零星散布在礦區(qū)各處，且多為建筑強(qiáng)度較低的窯洞或磚瓦房屋，上述安全采煤方式會(huì)使生產(chǎn)成本、作業(yè)效率和資源利用率均大幅降低。因此，當(dāng)前解決村落區(qū)煤炭開采問題大多采用村莊搬遷的方式，即在開采之前將村落搬遷到規(guī)劃的區(qū)域集中安置，并在地下煤層中留設(shè)一定范圍的保護(hù)煤柱，以確保新選址區(qū)的安全。

對(duì)于黃土高原采煤沉陷區(qū)村落搬遷的選址評(píng)價(jià)問題，涉及采礦、地理、地質(zhì)、人文、規(guī)劃等多學(xué)科理論，目前還沒有普適性的評(píng)價(jià)模型與方法。文中以陜西某煤礦開采沉陷區(qū)村落搬遷為實(shí)例模型，根據(jù)煤礦工作面布設(shè)計(jì)劃和開采沉陷預(yù)計(jì)理論，確定村落搬遷區(qū)范圍和遷址區(qū)保護(hù)煤柱的大小，采用運(yùn)籌學(xué)的層次分析法、模糊綜合評(píng)價(jià)法，結(jié)合特爾菲法，確定礦區(qū)村落地貌條件、水文狀況、交通條件、勞作范圍等多個(gè)影響因子的權(quán)重模型，形成綜合評(píng)價(jià)指標(biāo)體系[7-9]。基于 ArcGIS軟件平臺(tái)，定量地分析出沉陷區(qū)村落搬遷選址區(qū)的地理環(huán)境適宜度[10-12]，為西部黃土高原礦區(qū)村莊下壓煤的安全開采提供技術(shù)支持。

1采煤沉陷區(qū)村落遷址評(píng)價(jià)方法

西部礦井多采用走向長(zhǎng)壁工作面頂板全部垮落法開采。按照礦井生產(chǎn)設(shè)計(jì)要求，每年布設(shè)數(shù)個(gè)長(zhǎng)壁工作面有計(jì)劃地進(jìn)行回采，在地表逐漸形成由小到大的沉陷區(qū)域，而沉陷區(qū)內(nèi)的村落建筑也會(huì)不斷破壞。因此，在地下工作面計(jì)劃回采之前，必須先將地面沉陷區(qū)域內(nèi)的村落實(shí)施搬遷，新的搬遷選址區(qū)評(píng)價(jià)工作流程包括：制定開采工作面計(jì)劃、預(yù)計(jì)地表沉陷范圍、確定搬遷村落面積、確定選址區(qū)大小、選址區(qū)地理環(huán)境評(píng)價(jià)建模、選址區(qū)適宜度分析、適宜遷址區(qū)的確定。

1.1開采沉陷范圍預(yù)計(jì)與遷址區(qū)面積確定

每一個(gè)長(zhǎng)壁工作面開采后地表一般會(huì)形成沉陷盆地，地表沉陷變形的大小可依據(jù)開采沉陷預(yù)計(jì)理論進(jìn)行預(yù)計(jì)分析?，F(xiàn)有開采沉陷預(yù)計(jì)方法可分為影響函數(shù)法，剖面函數(shù)法，數(shù)值分析法和力學(xué)分析法等[13-19]。在我國(guó)，基于隨機(jī)介質(zhì)理論的概率積分法應(yīng)用最為普遍[20-24]，該法較精確合理地描述了地表沉陷盆地的分布形態(tài)，具有理論嚴(yán)密及便于編程計(jì)算等優(yōu)點(diǎn)。針對(duì)黃土高原礦區(qū)的特點(diǎn)，文獻(xiàn)[25]利用隨機(jī)介質(zhì)理論原理，將地表沉陷視為基巖面不均勻沉降在土層中傳播影響至地表的結(jié)果，導(dǎo)出概率積分法分層預(yù)計(jì)模型，并編制了專門的預(yù)計(jì)程序。實(shí)踐表明，只要按礦區(qū)具體地質(zhì)采礦條件選用合適的預(yù)計(jì)參數(shù)，上述模型完全能滿足黃土高原礦區(qū)開采沉陷預(yù)計(jì)和村落搬遷工作的要求。

在多個(gè)相鄰工作面開采后，工作面之間的間隔煤柱由于寬度僅為10～20 m左右，在覆巖移動(dòng)過程中會(huì)失去煤柱支撐作用，使地表形成整體性的沉陷盆地。按照上述概率積分法分層預(yù)計(jì)模型對(duì)各個(gè)計(jì)劃的回采工作面開采沉陷變形進(jìn)行預(yù)計(jì)。計(jì)算工作在專門開發(fā)的開采沉陷預(yù)計(jì)軟件平臺(tái)上完成。在ArcGIS中對(duì)各個(gè)工作面開采后的地表下沉和變形值，分別予以疊加并進(jìn)行可視化表達(dá)[26-28]。根據(jù)開采沉陷理論，按照地表下沉量小于10 mm來劃定地表沉陷區(qū)范圍。顯然，沉陷區(qū)范圍是按年度不斷增加的，沉陷區(qū)內(nèi)的村落面積也會(huì)逐年增加，根據(jù)計(jì)劃年限可以確定搬遷村落的總面積。在此基礎(chǔ)上，根據(jù)塌陷區(qū)村落建筑面積、相關(guān)搬遷安置政策和村鎮(zhèn)規(guī)劃原理來確定遷址區(qū)的規(guī)劃面積。

1.2遷址區(qū)地理適宜度評(píng)價(jià)模型

針對(duì)黃土高原村落的地理環(huán)境復(fù)雜多變且影響因子眾多的特點(diǎn)，文中選擇運(yùn)籌學(xué)中的AHP法，來構(gòu)建黃土礦區(qū)村落遷址區(qū)的地理適宜度評(píng)價(jià)模型。利用AHP法將定量分析與定性分析相結(jié)合，進(jìn)行黃土礦區(qū)村落搬遷選址評(píng)價(jià)，用專家經(jīng)驗(yàn)給出自然地理指標(biāo)和社會(huì)經(jīng)濟(jì)指標(biāo)在某一標(biāo)準(zhǔn)體系下的相對(duì)重要程度，并有效地計(jì)算出各個(gè)影響因子對(duì)應(yīng)的權(quán)重值，構(gòu)建地理適宜度評(píng)價(jià)模型[29-30]。

在黃土高原開采沉陷區(qū)村落遷址的評(píng)價(jià)中，利用AHP法進(jìn)行黃土高原沉陷區(qū)村落遷址的地理環(huán)境適宜性評(píng)價(jià)，將決策目標(biāo)、考慮因素和決策對(duì)象按其相互關(guān)系分為3個(gè)層次，其結(jié)構(gòu)模型如圖1所示，包括目標(biāo)層、準(zhǔn)則層、方案層。文中利用圖1所示的模型進(jìn)行黃土礦區(qū)村落遷址的綜合評(píng)價(jià)。

1.3地理適宜度評(píng)價(jià)因子及評(píng)分依據(jù)

針對(duì)黃土高原采煤沉陷區(qū)村落遷址問題，評(píng)價(jià)目標(biāo)為本礦區(qū)覆蓋范圍，首先要避開預(yù)計(jì)的開采沉陷區(qū)，以及不適宜居住建房的灘涂地、濕陷性黃土區(qū)、滑坡風(fēng)險(xiǎn)區(qū)等地質(zhì)條件惡劣和災(zāi)害易發(fā)區(qū)。選址區(qū)地理適宜度評(píng)價(jià)的影響因子可分為2類：自然地理指標(biāo)和社會(huì)經(jīng)濟(jì)指標(biāo)，細(xì)分為6個(gè)主要影響因子[31]，見表1.

地面坡度很大程度上影響居民點(diǎn)布局、土建費(fèi)用和土地利用效率等，是影響選址區(qū)地理適宜度的重要因子;坡向影響選址區(qū)日照時(shí)長(zhǎng)、太陽(yáng)輻射強(qiáng)度、建筑朝向和通風(fēng)條件等宜居環(huán)境;水源是影響農(nóng)村居民生活用水和農(nóng)業(yè)灌溉的重要因素，其中地下水位高低對(duì)村民打井用水的便利性至關(guān)重要;地表水系對(duì)農(nóng)業(yè)生產(chǎn)有重要影響，但黃土高原礦區(qū)水系的流域面積一般較小;遷址區(qū)距離交通干道的遠(yuǎn)近直接影響搬遷配套費(fèi)用和村民日常生活便利程度;沉陷區(qū)村落實(shí)施搬遷后，其農(nóng)田仍然位于原址，因此選址區(qū)距離搬遷區(qū)的距離（即村民勞作半徑）是影響村民生產(chǎn)、生活便利性的重要因素。對(duì)上述二級(jí)影響因子使用專家打分法和現(xiàn)場(chǎng)調(diào)查法進(jìn)行定級(jí)和打分。設(shè)置3個(gè)等級(jí)，對(duì)應(yīng)分值為20分、15分，10分。各二級(jí)影響因子評(píng)分依據(jù)及分值見表2.

計(jì)算矩陣A的最大特征值λmax及其對(duì)應(yīng)的特征向量B=[b1，b2，…，bn]，并對(duì)向量B做歸一化處理，即求得各因子的權(quán)重向量W=[W1，W2，…，Wn].

針對(duì)黃土高原采煤沉陷區(qū)村落選址的實(shí)際情況，確定自然地理一級(jí)指標(biāo)的相對(duì)權(quán)重為0.60，社會(huì)經(jīng)濟(jì)一級(jí)指標(biāo)的相對(duì)權(quán)重0.40.評(píng)價(jià)模型中自然地理及社會(huì)經(jīng)濟(jì)指標(biāo)下的各二級(jí)指標(biāo)間的相對(duì)重要關(guān)系見表3及表4.

indicators of physical geography影響因子坡度坡向距平均水位距地表水坡度1412坡向1/411/41/2距地下水位1413距地表水系1/221/31注：一致性指標(biāo)C.I.=0.006 9;平均隨機(jī)一致性指標(biāo)R.I.=0.890 0;隨機(jī)一致性比率C.R.=0.007 7.

由上述分析可知，各指標(biāo)間保持較好的相對(duì)一致性。采用上述方法計(jì)算出自然地理和社會(huì)經(jīng)濟(jì)中各二級(jí)指標(biāo)（因子）的權(quán)重向量及對(duì)應(yīng)的權(quán)重值，見表5.

針對(duì)評(píng)價(jià)區(qū)任意一個(gè)地理單元，根據(jù)表2的評(píng)分標(biāo)準(zhǔn)及表5所對(duì)應(yīng)的影響因子權(quán)重值，加權(quán)計(jì)算出各單元的綜合評(píng)價(jià)分值，通過GIS的疊加分析生成煤礦遷址區(qū)地理適宜度評(píng)價(jià)圖。將選址區(qū)的規(guī)劃面積作為搜索窗口，在評(píng)價(jià)圖中確定區(qū)域分值（取窗口平均值）最高的區(qū)域作為新的遷址區(qū)。

2.1塌陷搬遷區(qū)范圍確定

評(píng)價(jià)區(qū)為黃土高原的陜西某煤礦，開采煤層包含4上煤層和4下煤層。4上煤層開采厚度為3.8 m，4下煤層開采厚度為9 m，煤層平均傾角2°，開采深度350～450 m.礦區(qū)地表為黃土臺(tái)塬及溝壑地貌，中部和東部為溝壑區(qū)，其余部分較為平坦。地面高程最大1 270 m，最小750 m，地下水位平均高程約700 m.按照礦井生產(chǎn)規(guī)劃，在2020—2030年間計(jì)劃回采16個(gè)相鄰的工作面，地面將逐漸形成大范圍的塌陷盆地，必須提前實(shí)施村落搬遷。為此，根據(jù)文中提出的方法進(jìn)行村落遷址區(qū)的地理適宜度評(píng)價(jià)。本井田的井上下對(duì)照情況如圖2所示。

圖2中評(píng)價(jià)區(qū)范圍即井田邊界圈定的總面積約88 km2，2030年之前計(jì)劃回采區(qū)面積12 km2，地表塌陷范圍由概率積分法開采沉陷預(yù)計(jì)模型來確定。根據(jù)礦區(qū)地質(zhì)采礦條件采用類比法確定預(yù)計(jì)參數(shù)分別為：下沉系數(shù)q=0.7，主要影響角正切TANB=2.2，拐點(diǎn)偏距S0=0.05H，水平移動(dòng)系數(shù)b=0.3.其中平均采深H=400 m，煤層開采厚度m1=3.8 m，m2=9.0 m.按照上述預(yù)計(jì)參數(shù)用專門的開采沉陷預(yù)計(jì)軟件確定地表沉陷區(qū)范圍為14.6 km2，如圖2所示。沉陷區(qū)共有農(nóng)戶224戶，建筑總面積0.033×106 m2.按照新農(nóng)村規(guī)劃要求確定新的遷址區(qū)占地面積為0.1×106 m2.

2.2遷址區(qū)地理適宜度綜合評(píng)價(jià)結(jié)果

將圖2導(dǎo)入ArcGIS軟件中，將整個(gè)井田范圍（沉陷區(qū)除外）按照10 m分辨率進(jìn)行影響因子提取和評(píng)分，其中地下水位值用地面格網(wǎng)高程值減去區(qū)域平均水位而得到。將表5中每個(gè)因子的權(quán)重值賦予ArcGIS的各個(gè)圖層中，對(duì)6個(gè)二級(jí)影響因子進(jìn)行緩沖區(qū)分析以及疊加分析，通過聯(lián)合功能計(jì)算出每個(gè)地理單元的綜合分值。在ArcGIS中生成地理適宜度綜合分值的分級(jí)圖，如圖3所示。

由于文中對(duì)各影響因子的評(píng)價(jià)分值為10，15，20等3個(gè)等級(jí)，當(dāng)各因子分值均為10分時(shí)，則綜合評(píng)分值最低應(yīng)為50分，當(dāng)各因子分值均為20分時(shí)，綜合評(píng)分值最高為100分。圖3中綜合評(píng)價(jià)分值在52-97分之間，各分值區(qū)間面積占比見表6.其中綜合評(píng)價(jià)分值在90分以上的區(qū)域面積為034×106 m2，分值在85～90分的區(qū)域面積為627×106 m2.以遷址區(qū)規(guī)劃面積0.1×106 m2為搜索窗口，基于ArcGIS疊加得到最高分值區(qū)，如圖3中的A區(qū)，綜合評(píng)價(jià)平均分值為92分。

但是，該區(qū)域與待搬遷村落不屬于同一個(gè)縣級(jí)行政區(qū)，在拆遷實(shí)施和后續(xù)行政管理方面存在諸多不便。按照上述方法可搜索得到次高分值區(qū)，如圖3中的B區(qū)，綜合評(píng)價(jià)平均分值為87分。該區(qū)域與待搬遷村落屬于同一縣級(jí)行政區(qū)域，亦屬于地理適宜度較佳的遷址區(qū)域。將上述評(píng)價(jià)結(jié)果交由決策部門參考。

3結(jié)論

1）黃土高原煤礦區(qū)受開采影響村落的搬遷選址，主要考慮礦區(qū)地理環(huán)境、水文資源、交通便捷性、塌陷區(qū)相對(duì)位置等影響因素及其相對(duì)重要程度，形成綜合評(píng)價(jià)體系;

2）文中利用層次分析法結(jié)合專家打分法構(gòu)建了黃土溝壑區(qū)村落遷址區(qū)地理適宜度的綜合評(píng)價(jià)模型。應(yīng)用實(shí)例表明，基于該模型通過ArcGIS疊加分析可便捷地進(jìn)行村落遷址區(qū)地理環(huán)境的量化評(píng)價(jià)，其結(jié)果具有科學(xué)性與合理性;

3）黃土高原煤礦區(qū)村落搬遷選址的科學(xué)評(píng)價(jià)提供了有效的技術(shù)方法，具有實(shí)際應(yīng)用價(jià)值。

參考文獻(xiàn)（References）：

[1]范立民，向茂西，彭捷，等.西部生態(tài)脆弱礦區(qū)地下水對(duì)高強(qiáng)度采煤的響應(yīng)[J].煤炭學(xué)報(bào)，2016，41（11）：2672-2678.

FAN Limin，XIANG Maoxi，PENG Jie，et al.Groundwater response to intensive mining in ecologically fragile area[J].Journal of China Coal Society，2016，41（11）：2672-2678.

[2]GUO Qingbiao，GUO Guangli，ZHA Jianfeng，et al.Research on the surface movement in a mountain mining area： a case study of Sujiagou Mountain，China[J].Environmental Earth Sciences，2016，75（6）：472-490.

[3]ZHOU Dawei，WU Kan，CHENG Gonglin，et al.Mechanism of mining subsidence in coal mining area with thick alluvium soil in China[J].Arabian Journal of Geosciences，2015，8（4）：1855-1867.

[4]湯伏全，白峰.厚黃土層礦區(qū)采動(dòng)失水引起的地表沉降計(jì)算方法[J].西安科技大學(xué)學(xué)報(bào)，2011，31（4）：448-452.

TANG Fuquan，BAI Feng.Calculation method of surface subsidence caused by water loss in thick loess mining area[J].Journal of Xi’an University of Science and Technology，2011，31（4）：448-452.

[5]范立民，馬雄德，冀瑞君.西部生態(tài)脆弱礦區(qū)保水采煤研究與實(shí)踐進(jìn)展[J].煤炭學(xué)報(bào)，2015，40（8）：1711-1717.

FAN Limin，MA Xiongde，JI Ruijun.Progress in engineering practice of waterpreserved coal mining in western ecoenvironment frangible area[J].Journal of China Coal Society，2015，40（8）： 1711-1717.

[6]戴華陽(yáng)，郭俊廷，閻躍觀，等.“采-充-留”協(xié)調(diào)開采技術(shù)原理與應(yīng)用[J].煤炭學(xué)報(bào)，2014，39（8）：1602-1610.

DAI Huayang，GUO Junting，YAN Yueguan，et al.Principe and application of subsidence control technology of mining coordinately mixed with backfilling and keeping[J].Journal of China Coal Society，2014，39（8）：1602-1610.

[7]丁紹剛，朱嫣然.基于層次分析法與模糊綜合評(píng)價(jià)法的醫(yī)院戶外環(huán)境綜合評(píng)價(jià)體系構(gòu)建[J].浙江農(nóng)林大學(xué)學(xué)報(bào)，2017，34（6）：1104-1112.

DING Shaogang，ZHU Yanran.Constructing a synthetic evaluation system of hospital outdoor environment based on analytic hierarchy process and fuzzy synthetic evaluation[J].Journal of Zhejiang A&F University，2017，34（6）：1104-1112.

[8]連海波，趙法鎖，王雁林，等.陜南移民搬遷安置區(qū)選址適宜性評(píng)價(jià)模型及實(shí)例應(yīng)用[J].災(zāi)害學(xué)，2017，32（4）：120-125.

LIAN Haibo，ZHAO Fasuo，WANG Yanlin，et al.Immigrant relocation area in southern Shaanxi site suitability evaluation model and application instance[J].Journal of Catastrophology，2017，32（4）：120-125.

[9]孫穎，邢文杰，游?？?基于模糊層次分析法的公路梁式橋地震震害預(yù)測(cè)研究[J].福州大學(xué)學(xué)報(bào)（自然科學(xué)版），2013，41（4）：776-781.

SUN Yin，XING Wenjie，YOU Fuke.The research on the earthquake damage prediction of highway girder bridges based on fuzzy analytical hierarchy approach[J].Journal of Fuzhou University（Natural Science Edition），2013，41（4）：776-781.

[10]WU Quanyuan，PANG Jiewu，QI Shanzhong，et al.Impacts of coal mining subsidence on the surface landscape in Longkou city，Shandong Province of China[J].Environmental Earth Sciences，2009，59（4）：783-791.

[11]Yenilmez Firdes，Kuter Nazan，Emil MustafaKemal，et al.Evaluation of pollution levels at an abandoned coal mine site in Turkey with the aid of GIS[J].International Journal of Coal Geology，2011，86（1）：12-19.

[12]尚慧，倪萬(wàn)魁，羅東海.寧南山區(qū)泥石流災(zāi)害危險(xiǎn)性量化評(píng)價(jià)研究[J].西安科技大學(xué)學(xué)報(bào)，2017，37（1）：63-70.

SHANG Hui，NI Wankui，LUO Donghai，et al.Quantitative assessment of debris flow hazards in the mountainareas of southern Ningxia[J].Journal of Xi’an University of Science and Technology，2017，37（1）：63-70.

[13]ZHANG Jixiong，ZHANG Qiang，Sun Qiang，et al.Surface subsidence control theory and application to backfill coal mining technology[J].Environmental Earth Sciences，2015，74（2）：1439-1448.

[14]Howladar M Farhad，Hasan Kamrul.A study on the development of subsidence due to the extraction of 1203 slice with its associated factors around Barapukuria underground coal mining industrial area，Dinajpur，Bangladesh[J].Environmental Earth Sciences，2014，72（9）：3699-3713.

[15]lvarezFernndez M I，GonzlezNicieza C，MenéndezDíaz? A，et al.Generalization of the nk influence function to predict mining subsidence[J].Engineering Geology，2005，80：1-36.

[16]ZHU Xiaojun，GUO Guangli，F(xiàn)ANG Qi.Coupled discrete element finite difference method for analyzing subsidence control in fully mechanized solid backfilling mining[J].Environmental Earth Sciences，2016，75（8）：683-694.

[17]YANG Dejun，BIAN Zhenfu，LEI Shaogang.Impact on soil physical qualities by the subsidence of coal mining： a case study in Western China[J].Environmental Earth Sciences，2016，75（8）：652-665.

[18]王小華，胡海峰，廉旭剛.基于Weibull時(shí)間序列函數(shù)與負(fù)指數(shù)法的動(dòng)態(tài)沉陷預(yù)計(jì)[J].金屬礦山，2015，44（8）：158-162.

WANG Xiaohua，HU Haifeng，LIAN Xugang.Dynamic subsidence prediction method based on negative exponential method and Weibull time sequence function[J].Metal Mine，2015，44（8）： 158-162.

[19]WU Xiong，JIANG Xiaowei，CHEN Yufu，et al.The influences of mining subsidence on the ecological environment and public infrastructure： a case study at the Haolaigou Iron Ore Mine in Baotou，China[J].Environmental Earth Sciences，2009，59（4）：803-810.

[20]郭麒麟，喬世范，劉寶琛.開采影響下的巖土體移動(dòng)與變形規(guī)律[J].采礦與安全工程學(xué)報(bào)，2011，28（1）：109-114.

GUO Qilin，QIAO Shifan，LIU Baochen.The movement and deformation of soil and rock mass resulted from mining activity[J].Journal of Mining and Safety Engineering，2011，28（1）：109-114.

[21]劉寶琛，戴華陽(yáng).概率積分法的由來與研究進(jìn)展[J].煤礦開采，2016，21（2）：1-3.

Liu Baochen，DAI Huayang.Research development and origin of probability integral method[J].Coal Mining Technology，2016，21（2）：1-3.

[22]龔云，湯伏全.西部黃土山區(qū)開采沉陷變形數(shù)值模擬研究[J].西安科技大學(xué)學(xué)報(bào)，2012，32（4）：490-494.

GONG Yun，TANG Fuquan.Numerical simulation of mining subsidence in the western loess mountain area[J].Journal of Xi’an University of Science and Technology，2012，32（4）：490-494.

[23]趙曉東，陳陽(yáng)，蔣建.任意形狀工作面沉陷預(yù)測(cè)的概率積分法及其應(yīng)用[J].巖土力學(xué)，2016（12）：3387-3392，3400.

ZHAO Xiaodong，CHEN Yang，JIANG Jian.Probability integral method of arbitrary shape face to predict mining subsidence and its application[J].Rock and Soil Mechanics，2016（12）：3387-3392，3400.

[24]李培現(xiàn)，譚志祥，閆麗麗，等.基于支持向量機(jī)的概率積分法參數(shù)計(jì)算方法[J].煤炭學(xué)報(bào)，2010，35（8）：1247-1251.

LI Peixian，TAN Zhixiang，YAN Lili，et al.Calculation method of probability integration method parameters based on support vector machine[J].Journal of China Coal Society，2010，35（8）：1247-1251.

[25]湯伏全.西部厚黃土層礦區(qū)開采沉陷預(yù)計(jì)模型[J].煤炭學(xué)報(bào)，2011，36（增刊1）：74-78.

TANG Fuquan.Mining subsidence prediction model in western thick loess layer mining areas[J].Journal of China Coal Society，2011，36（Supp.1）：74-78.

[26]程鋼，張富杰，楊杰，等.GIS支持下的城市小區(qū)的居住適宜性分析研究[J].測(cè)繪通報(bào)，2016（1）：49-52.

CHENG Gang，ZHANG Fujie，YANG Jie，et al.Evaluation for livability for urban residential areas based on GIS[J].Bulletin of Surveying and Mapping，2016（1）：49-52.

[27]Gupta D K.Comprehensive suitability evaluation of tea crops using GIS and a modified land ecological suitability evaluation model[J].Pedosphere，2012，22（1）：122-130.

[28]金貴，王占岐，胡學(xué)東，等.基于模糊證據(jù)權(quán)模型的青藏高原區(qū)土地適宜性評(píng)價(jià)[J].農(nóng)業(yè)工程學(xué)報(bào)，2013，29（18）：241-250.

JIN Gui，WANG Zhanqi，HU Xuedong，et al.Land suitability evaluation in QinghaiTibet Plateau based on fuzzy weight of evidence model[J].Transactions of the Chinese Society of Agricultural Engineering，2013，29（18）： 241-250.

[29]章正高.采用層次分析法對(duì)地質(zhì)災(zāi)害進(jìn)行危險(xiǎn)性評(píng)價(jià)——以漢陰縣地質(zhì)災(zāi)害為例[J].西安科技大學(xué)學(xué)報(bào)，2014，34（1）：87-92.

ZHANG Zhenggao.Risk assessment of geologic hazards by analytic hierarchy process：a case study of the geologic hazard of Hanyin county[J].Journal of Xi’an University of Science and Technology，2014，34（1）：87-92.

[30]李寧，王李管，賈明濤.基于層次分析法的礦井六大系統(tǒng)模糊綜合評(píng)價(jià)[J].中南大學(xué)學(xué)報(bào)（自然科學(xué)版），2015，46（2）：631-637.

LI Ning，WANG Liguan，JIA Mingtao.An analytic hierarchy process based fuzzy evaluation of underground mine sixsystem[J].Journal of Central South University （Science and Technology），2015，46（2）：631-637.

[31]Giulia Dapueto，F(xiàn)rancesco Massa，Sara Costa，et al.A spatial multicriteria evaluation for site selection of offshore marine fish farm in the Ligurian Sea，Italy[J].Ocean and Coasta Management，2015，11（6）：64-77.