多煤層開(kāi)采條帶錯(cuò)動(dòng)程度對(duì)巖層控制的影響

王 冰 郭廣禮 朱曉峻 王 炯

(1.中國(guó)礦業(yè)大學(xué)環(huán)境與測(cè)繪學(xué)院，江蘇 徐州 221008；2.國(guó)土環(huán)境與災(zāi)害監(jiān)測(cè)國(guó)家測(cè)繪地理信息局重點(diǎn)實(shí)驗(yàn)室，江蘇 徐州 221008；3.江蘇省資源環(huán)境信息工程重點(diǎn)實(shí)驗(yàn)室，江蘇 徐州 221008)

多煤層開(kāi)采條帶錯(cuò)動(dòng)程度對(duì)巖層控制的影響

王 冰1,2,3郭廣禮1,2,3朱曉峻1,2,3王 炯1,2,3

(1.中國(guó)礦業(yè)大學(xué)環(huán)境與測(cè)繪學(xué)院，江蘇 徐州 221008；2.國(guó)土環(huán)境與災(zāi)害監(jiān)測(cè)國(guó)家測(cè)繪地理信息局重點(diǎn)實(shí)驗(yàn)室，江蘇 徐州 221008；3.江蘇省資源環(huán)境信息工程重點(diǎn)實(shí)驗(yàn)室，江蘇 徐州 221008)

針對(duì)我國(guó)多煤層條帶開(kāi)采的實(shí)際情況，提出了煤層錯(cuò)動(dòng)系數(shù)的概念。利用FLAC3D數(shù)值模擬軟件，系統(tǒng)地研究了多煤層條帶開(kāi)采時(shí)，上下煤層錯(cuò)動(dòng)系數(shù)的大小對(duì)上下條帶保護(hù)煤柱應(yīng)力、中間巖體應(yīng)力以及地表移動(dòng)變形的影響規(guī)律，分別給出了相應(yīng)的函數(shù)關(guān)系式，并指出了在實(shí)際多煤層條帶開(kāi)采時(shí)，極限煤層錯(cuò)動(dòng)系數(shù)的求取方法。研究成果對(duì)于指導(dǎo)多煤層條帶開(kāi)采工作面設(shè)計(jì)具有十分重要的理論和現(xiàn)實(shí)意義。

FLAC3D數(shù)值模擬 多煤層條帶開(kāi)采 煤層錯(cuò)動(dòng)系數(shù)

條帶開(kāi)采技術(shù)是解決“三下”采煤?jiǎn)栴}的重要開(kāi)采技術(shù)之一，因其能有效地控制巖層及地表的下沉與變形，且生產(chǎn)管理簡(jiǎn)單，被我國(guó)煤礦區(qū)廣泛利用[1]。針對(duì)條帶開(kāi)采，近年來(lái)國(guó)內(nèi)外進(jìn)行了大量的研究并取得了一些成果[2-6]，但大多數(shù)僅僅是針對(duì)單煤層條帶開(kāi)采的研究，對(duì)于多煤層條帶開(kāi)采研究不多。鄧喀中教授[7-8]采用彈塑性有限單元模擬法得到了上下煤柱空間位置與巖層及地表移動(dòng)間的關(guān)系，并進(jìn)一步提出了煤柱穩(wěn)定性的評(píng)價(jià)方法；張俊英[9-10]采用相似材料模型實(shí)驗(yàn)與有限元數(shù)值模擬相結(jié)合的綜合研究方法，揭示了多煤層條帶過(guò)程中的開(kāi)采覆巖破壞和地表移動(dòng)變形規(guī)律；胡炳南[11]從地質(zhì)采礦因素對(duì)開(kāi)采沉陷的影響程度出發(fā),提出了煤層群條帶開(kāi)采中選擇開(kāi)采煤層的最小采動(dòng)影響原則、條帶布置的層間距和層間條帶煤柱的對(duì)齊原則以及條帶參數(shù)計(jì)算的“小變形”準(zhǔn)則等新內(nèi)容。

在研究上下煤層空間位置對(duì)巖層以及地表移動(dòng)變形的影響機(jī)理時(shí)，大部分文獻(xiàn)只是將上下條帶工作面的位置關(guān)系簡(jiǎn)單地分為“對(duì)齊”、“半對(duì)齊”以及“完全錯(cuò)開(kāi)”3種情況或是簡(jiǎn)單地分為錯(cuò)開(kāi)10、20 m等具體數(shù)值。這種對(duì)上下條帶工作面空間位置的分類(lèi)過(guò)于簡(jiǎn)單，往往達(dá)不到研究的目的，或者會(huì)因引入條帶開(kāi)采的采寬面干擾對(duì)工作面空間位置影響的研究。基于以上缺點(diǎn)，本研究考慮利用煤層錯(cuò)動(dòng)系數(shù)來(lái)說(shuō)明多煤層條帶開(kāi)采時(shí)上下條帶工作面的位置關(guān)系對(duì)巖層及地表下沉、變形的影響機(jī)理，煤層錯(cuò)動(dòng)系數(shù)的定義為上下煤層錯(cuò)動(dòng)距離與上煤層條帶采寬的比值。從定義中可以看出，煤層錯(cuò)動(dòng)系數(shù)很好地將條帶采寬這一因素排除在外，實(shí)現(xiàn)了研究影響因素的單一性。

FLAC3D是一種以巖石力學(xué)理論為基礎(chǔ)、以各巖層以及表土層的物理力學(xué)參數(shù)和地層的構(gòu)造特性為基本計(jì)算依據(jù)的數(shù)值模擬方法，它可以靈活有效地表現(xiàn)出各種不同煤層開(kāi)挖情況下引起的地表或巖體內(nèi)部的移動(dòng)變形情況。因此，為了得到煤層錯(cuò)動(dòng)系數(shù)與巖層及地表下沉變形的關(guān)系，本研究結(jié)合某礦地質(zhì)采礦條件背景，采用FLAC3D有限差分軟件建立了不同方案的數(shù)值模型，討論了煤層錯(cuò)動(dòng)系數(shù)變化時(shí)，地表移動(dòng)下沉、煤層間巖體的應(yīng)力以及上下條帶開(kāi)采煤柱的應(yīng)力變化規(guī)律。

1 上下煤層條帶開(kāi)采模擬方案

1.1 模型的建立以及邊界條件

某礦區(qū)可采煤層從上到下分別有5、8這2個(gè)煤層,平均厚度分別為3、4 m，層間距為32 m，底板厚30 m，礦區(qū)整體煤層平均傾角僅為3°，可按水平煤層考慮。設(shè)計(jì)上下煤層開(kāi)采寬度為60 m，留寬90 m，采出率為40%，上下均開(kāi)采3個(gè)條帶工作面。模擬巖層總厚度為501 m，走向長(zhǎng)度為1 000 m，傾向?qū)挾葹?00 m。模型共392 000個(gè)單元，482 679個(gè)節(jié)點(diǎn)。為方便研究，x方向和y方向均為每10 m 1個(gè)單元格，幾何模型如圖1所示。整體模型均采用摩爾-庫(kù)侖(Mohr-Coulomb)模型[12]，根據(jù)現(xiàn)場(chǎng)實(shí)際條件，選取的模擬參數(shù)見(jiàn)表1。

圖1 開(kāi)采煤層幾何模型Fig.1 Geometric model of coal seam mining

在模型左右邊界，固定x方向位移為零；在模型前后邊界，固定y方向位移為零；在模型下邊界，固定x、y方向位移均為零；模型上邊界為自由變形邊界。在分析計(jì)算過(guò)程中,不考慮構(gòu)造應(yīng)力對(duì)原巖應(yīng)力的影響,僅考慮巖體自重引起的應(yīng)力。巖體內(nèi)部初始應(yīng)力狀態(tài)取決于上覆巖層的重力和性質(zhì)。

表1 模擬礦區(qū)巖石力學(xué)參數(shù)Table 1 Rock mechanics parameters in simulated mining area

1.2 數(shù)值模擬方案

為研究煤層錯(cuò)動(dòng)系數(shù)這一單因素對(duì)巖層及地表下沉的影響，本研究在固定上下煤層采留寬為60和90 m即采出率均為40%、采深432 m、層間距32 m的基礎(chǔ)上，每次模擬時(shí)，上煤層條帶開(kāi)采位置保持不變，下煤層條帶開(kāi)采工作面距離上煤層條帶工作面錯(cuò)開(kāi)距離依次為10、20、30、40、50、60 m(此時(shí)為完全錯(cuò)開(kāi))，即算上工作面對(duì)齊共7個(gè)模擬方案。根據(jù)定義可算出各方案煤層錯(cuò)動(dòng)系數(shù)如表2所示。

表2 數(shù)值模擬方案Table 2 The numerical simulation program

注：其他影響因素為采留寬60 m和90 m、采深432 m、層間距32 m。

2 數(shù)值模擬結(jié)果分析

2.1 煤層錯(cuò)動(dòng)系數(shù)對(duì)上下條帶煤柱的影響

由條帶開(kāi)采的理論可知，中間留設(shè)的煤柱受力最大，因此本次研究以中間的條帶煤柱為例進(jìn)行分析。數(shù)值模擬結(jié)果表明：上、下煤層留設(shè)的條帶煤柱所受應(yīng)力兩端較大，向中間逐漸減小并趨于穩(wěn)定，呈馬鞍狀分布。中間條帶煤柱的最大壓應(yīng)力值與錯(cuò)動(dòng)系數(shù)呈二次函數(shù)關(guān)系(如圖2)，隨著煤層錯(cuò)動(dòng)系數(shù)的增加，兩煤層中間條帶煤柱壓應(yīng)力均增加，最后趨于某一穩(wěn)定值。當(dāng)壓應(yīng)力穩(wěn)定時(shí)，相對(duì)于煤柱對(duì)齊時(shí)，所受壓應(yīng)力增加較大。因此錯(cuò)動(dòng)系數(shù)的增大會(huì)進(jìn)一步減弱上下煤柱的穩(wěn)定性，嚴(yán)重影響煤礦的安全生產(chǎn)。

圖2 上下煤層壓應(yīng)力值與錯(cuò)動(dòng)系數(shù)關(guān)系Fig.2 The relationship of coal dislocation coefficient with pressure stress value of upper and lower coal seam

根據(jù)模擬結(jié)果，可建立上下煤柱壓應(yīng)力σ上,σ下與煤層錯(cuò)動(dòng)系數(shù)s的回歸關(guān)系式為

σ上=-21.476s2+40.116s+17.907,

(1)

σ下=-6.878s2+17.886s+16.434.

(2)

2.2 煤層錯(cuò)動(dòng)系數(shù)對(duì)中間巖體的影響

數(shù)值模擬結(jié)果表明，中間巖體的最大拉應(yīng)力值與最大壓應(yīng)力值均與錯(cuò)動(dòng)系數(shù)呈指數(shù)函數(shù)關(guān)系(如圖3)，隨著煤層錯(cuò)動(dòng)系數(shù)的增加，拉應(yīng)力與壓應(yīng)力均增加，且拉應(yīng)力增加較大，使中間巖體穩(wěn)定性下降。這是因?yàn)殡S著下煤層開(kāi)采工作面的錯(cuò)動(dòng)，使位于上煤層煤柱下面的中間巖體位于下煤層的煤柱之上，而是慢慢位于下煤層開(kāi)采工作面之上，中間巖體下部得不到支撐，致使上煤層煤柱對(duì)中間巖體的壓力以及下煤層開(kāi)采工作面對(duì)中間巖體的拉力都隨著煤層錯(cuò)動(dòng)系數(shù)的增大而增大。

根據(jù)模擬結(jié)果，可建立中間巖體拉應(yīng)力σ拉和壓應(yīng)力σ壓與煤層錯(cuò)動(dòng)系數(shù)s的回歸關(guān)系式為

(3)

(4)

2.3 煤層錯(cuò)動(dòng)系數(shù)對(duì)地表下沉的影響

從數(shù)值模擬結(jié)果可以看出，各方案地表下沉等值線均呈近似圓形分布，最大下沉值隨錯(cuò)動(dòng)系數(shù)的增大近似呈線性函數(shù)增大趨勢(shì)(如圖4)。由2.1節(jié)和2.2節(jié)可知，隨著錯(cuò)動(dòng)系數(shù)的增大，上下煤柱以及中間巖體的壓應(yīng)力逐漸增大，使得煤柱與中間巖體的壓縮量增大，因此地表下沉值也隨之增大。

圖3 中間巖體拉應(yīng)力和壓應(yīng)力值與錯(cuò)動(dòng)系數(shù)關(guān)系Fig.3 The relationship of coal dislocation coefficient with the tensile stress and compressive stress of intermediate rock

圖4 地表下沉值與錯(cuò)動(dòng)系數(shù)關(guān)系Fig.4 The relationship of coal dislocation coefficient with surface subsidence value

由圖4可以看出，上下條帶采空區(qū)完全對(duì)齊時(shí)，地表下沉值為174.13 mm，當(dāng)完全錯(cuò)開(kāi)即錯(cuò)動(dòng)系數(shù)為1時(shí)，地表最大下沉值為309.45 mm，最大下沉值幾乎增大了1倍，說(shuō)明多煤層條帶開(kāi)采上下工作面是否對(duì)齊對(duì)地表下沉有著極大的影響。因此，在布設(shè)條帶開(kāi)采工作面時(shí)，應(yīng)盡量使上下工作面在豎直方向上保持一致。根據(jù)模擬結(jié)果，可建立地表最大下沉值W與煤層錯(cuò)動(dòng)系數(shù)的回歸關(guān)系式為

W=137.35s+175.81.

(5)

2.4 煤層錯(cuò)動(dòng)系數(shù)對(duì)地表水平移動(dòng)的影響

根據(jù)數(shù)值模擬結(jié)果，得出地表水平移動(dòng)最大值隨著錯(cuò)動(dòng)系數(shù)也近似呈線性增大趨勢(shì)(如圖5)，因此在實(shí)際布設(shè)條帶開(kāi)采工作面時(shí)，使上下工作面在豎直方向上保持一致，不僅能夠有效地控制地表下沉，而且可以很好地控制水平移動(dòng)，這對(duì)于“三下”采煤尤其是建筑物下采煤有非常大的實(shí)際意義。

圖5 地表水平移動(dòng)值與錯(cuò)動(dòng)系數(shù)關(guān)系Fig.5 The relationship of coal dislocation coefficient with the surface horizontal movement value

根據(jù)模擬結(jié)果，可建立地表最大水平移動(dòng)值U與煤層錯(cuò)動(dòng)系數(shù)s的回歸關(guān)系式為

U=17.387s+29.292.

(6)

2.5 極限煤層錯(cuò)動(dòng)系數(shù)的求取

多煤層條帶開(kāi)采過(guò)程中，在條帶采留寬一定的情況下，開(kāi)采時(shí)必須保證上下煤層煤柱以及中間巖體不破壞。由上面的分析可以看出，煤層錯(cuò)動(dòng)系數(shù)越大，中間巖體以及煤柱越不穩(wěn)定，因此在實(shí)際布設(shè)工作面時(shí)要求出極限煤層錯(cuò)動(dòng)系數(shù)大小。

考慮下煤層煤柱的穩(wěn)定性時(shí)[8]，由于上下煤柱不對(duì)齊，導(dǎo)致承受的覆巖重力不一樣，因此要考慮煤層間受力的差異。左右煤柱受力分別為

(7)

(8)

式中，q為上層煤柱的載荷集度，t/m；l為上下煤層錯(cuò)動(dòng)距離，m；a2、b2分別為下煤層留寬和采寬，m；γ為覆巖平均容重，N/m3；h為上下煤層距離，m；[σ]為煤層單軸抗壓強(qiáng)度，MPa。

考慮中間巖體的穩(wěn)定性時(shí)，由梁理論可知:當(dāng)

巖層間最大應(yīng)力為

(9)

當(dāng)

巖層間最大應(yīng)力為

(10)

式中，

由于巖體的抗拉強(qiáng)度小于抗壓和抗剪強(qiáng)度，所以層間巖體的破壞主要為拉破壞，因此必須滿(mǎn)足式(11)才能保證中間巖體的穩(wěn)定性：

(11)

式中，[σ′]為層間巖體的平均抗拉強(qiáng)度。

在考慮上煤層煤柱的穩(wěn)定性時(shí)，由于下煤層煤柱以及中間巖體的穩(wěn)定性已得到保證，因此這時(shí)上煤層的煤柱穩(wěn)定性可按單一煤層條帶煤柱的穩(wěn)定性進(jìn)行考慮[13]，即煤柱在單向或者三向受力狀態(tài)下的實(shí)際載荷強(qiáng)度要小于極限載荷強(qiáng)度。因在前期工作面采留寬設(shè)計(jì)的時(shí)候已經(jīng)考慮該因素，所以這時(shí)上煤層煤柱的穩(wěn)定性不受煤層錯(cuò)動(dòng)系數(shù)的影響。

將式(7)～式(11)聯(lián)立，可求得在保證下煤層煤柱以及中間巖體均穩(wěn)定的前提下的最大煤層錯(cuò)動(dòng)距離 ，進(jìn)而根據(jù)定義可以求得極限煤層錯(cuò)動(dòng)系數(shù)。

通過(guò)上面分析，得出了多煤層條帶開(kāi)采條件下，地表移動(dòng)變形、中間巖體的應(yīng)力以及上下煤層的應(yīng)力均不同程度地隨煤層錯(cuò)動(dòng)系數(shù)增大而增大。因此在實(shí)際進(jìn)行多煤層條帶開(kāi)采時(shí)，應(yīng)盡量保持上下條帶工作面的一致性。

3 結(jié) 論

(1)中間巖體的拉應(yīng)力和壓應(yīng)力均與煤層錯(cuò)動(dòng)系數(shù)呈指數(shù)函數(shù)關(guān)系，隨著錯(cuò)動(dòng)系數(shù)的增加，中間巖體拉應(yīng)力及壓應(yīng)力均增大，但拉應(yīng)力增加較大，壓應(yīng)力增加較小。

(2)上、下條帶煤柱的壓應(yīng)力兩端大、中間小，呈馬鞍狀分布，且與煤層錯(cuò)動(dòng)系數(shù)均呈二次函數(shù)關(guān)系，隨著煤層錯(cuò)動(dòng)系數(shù)的增加，兩煤層壓應(yīng)力先增加，且增加量較大，最后趨于某一穩(wěn)定值不再變化。

(3)地表下沉值與煤層錯(cuò)動(dòng)系數(shù)呈線性函數(shù)關(guān)系，隨著錯(cuò)動(dòng)系數(shù)的增大，地表下沉值增大，且變化量較大，當(dāng)煤柱完全錯(cuò)開(kāi)時(shí)，地表下沉值相對(duì)于煤柱對(duì)齊時(shí)增大了近1倍。

(4)地表水平移動(dòng)值隨著煤層錯(cuò)動(dòng)系數(shù)也近似呈線性函數(shù)關(guān)系，隨著錯(cuò)動(dòng)系數(shù)的增大，地表水平移動(dòng)值增大。

[1] 郭文兵.深部大采寬條帶開(kāi)采地表移動(dòng)的預(yù)計(jì)[J].煤炭學(xué)報(bào),2008,33(4):368-372. Guo Wenbing.Surface movement predicting problems of deep strip pillar mining[J].Journal of China Coal Society,2008,33(4):368-372.

[2] 劉 貴,張華興,徐乃忠.深部厚煤層條帶開(kāi)采地表沉陷規(guī)律模擬研究[J].礦山測(cè)量,2008(1):53-55. Liu Gui,Zhang Huaxing,Xu Naizhong.Simulation research on deep and thick coal seam strip mining surface subsidence[J].Mine Surveying,2008(1):53-55.

[3] 郭增長(zhǎng),謝和平,王金莊.條帶開(kāi)采保留煤柱寬度和采出寬度與地表變形的關(guān)系[J].湘潭礦業(yè)學(xué)院學(xué)報(bào),2003,18(2):13-17. Guo Zengzhang,Xie Heping,Wang Jinzhuang.The relationship of surface deformation with reserved width and mining width of coal pillar in strip mining[J].Journal of Xiangtan Mining Institute,2003,18(2):13-17.

[4] 王金莊,吳立新.建(構(gòu))筑物下條帶開(kāi)采覆巖移動(dòng)機(jī)理與實(shí)踐[J].煤,1998(6):6-7. Wang Jinzhuang,Wu Lixin.Mechanism and practice of overlying rock in building(structure) building materials strip mining[J].Coal,1998(6):6-7.

[5] 郭文兵.條帶開(kāi)采的非線性理論模型與應(yīng)用研究[D].徐州:中國(guó)礦業(yè)大學(xué),2004. Guo Wenbing.Nonlinear Theory Model and Application Research of Strip Mining[D].Xuzhou:China University of Mining Technology,2004.

[6] 郭文兵,劉義新.深部條帶開(kāi)采下沉系數(shù)與采厚關(guān)系的數(shù)值模擬[J].河南理工大學(xué)學(xué)報(bào):自然科版,2007,26(3):254-258. Guo Wenbing,Liu Yixin.Deep relationship between strip mining subsidence coefficient and mining thickness of numerical simulation[J].Journal of Henan Polytechnic University:Natural Science，2007,26(3):254-258.

[7] 鄧喀中,馬偉民,何國(guó)清.多煤層開(kāi)采時(shí)條帶空間位置對(duì)巖層移動(dòng)的影響[J].中國(guó)礦業(yè)大學(xué)學(xué)報(bào),1991,20(2):75-81. Deng Kazhong,Ma Weimin,He Guoqing.Influence of stripe space position on strata movement in multiple coal seams mining[J].Journal of China University of Mining & Technology,1991,20(2):75-81.

[8] 鄧喀中,范洪冬,譚志祥,等.多煤層條帶開(kāi)采煤柱穩(wěn)定性評(píng)價(jià)方法研究[J].中國(guó)科技論文在線,2009,4(11):824-829. Deng Kazhong,Fan Hongdong,Tan Zhixiang,et al.Evaluation method for pillar stability during strip mining in multiple coal seams[J].Science Paper Online,2009,4(11):824-829.

[9] 張俊英.多煤層條帶開(kāi)采實(shí)驗(yàn)?zāi)M研究[J].礦山測(cè)量,2000(1):28-30. Zhang Junying.The experimental simulation studies in multiple coal seams strip mining[J].Mine Surveying,2000(1):28-30.

[10] 張俊英.多煤層條帶開(kāi)采模擬理論研究[J].煤炭學(xué)報(bào),2000,25(S1):67-70. Zhang Junying.Simulation theory research in multiple coal seams strip mining[J].Journal of China Coal Society,2000,25(S1):67-70.

[11] 胡炳南.煤層群條帶開(kāi)采優(yōu)化設(shè)計(jì)原則[J].煤礦設(shè)計(jì),2000(4):3-6. Hu Bingnan.Strip mining coal seam group optimization design principle[J].Design of Coal Mine,2000(4):3-6.

[12] 查劍鋒,郭廣禮,劉元旭,等.矸石變形非線性及其對(duì)巖層移動(dòng)的影響[J].煤炭學(xué)報(bào),2009,34(8):1071-1075. Zha Jianfeng,Guo Guangli,Liu Yuanxu,et al.The nonlinear features of waste deformation and its impaction on strata movement[J].Journal of China Coal Society,2009,34(8):1071-1075.

[13] 何國(guó)清,楊 倫,凌賡娣,等.礦山開(kāi)采沉陷學(xué)[M].徐州:中國(guó)礦業(yè)大學(xué)出版社,1990. He Guoqing,Yang Lun,Ling Gengdi,et al.Mining Subsidence in Mine[M].Xuzhou：China University of Mining and Technology Press,1994.

(責(zé)任編輯 徐志宏)

Influence of Strip Dislocation Degree on the Strata Control in Multi-coal Seam Mining

Wang Bing1,2,3Guo Guangli1,2,3Zhu Xiaojun1,2,3Wang Jiong1,2,3

(1.SchoolofEnvironmentScienceandSpatialInformatics，ChinaUniversityofMiningandTechnology，Xuzhou221008，China;2.KeyLaboratoryforLandEnvironmentandDisasterMonitoringofSBSM，Xuzhou221008，China;3.TheMainLaboratoryofResourceEnvironmentInformationofJiangsu，Xuzhou221008，China)

In view of the actual situation of multi-coal seam strip mining in China，the concept of coal dislocation coefficient is put forward.Based on the FLAC3Dnumerical simulation software，the influencing laws of different coal dislocation coefficients on the stress of upper and lower bands protection coal pillar，the stress of intermediate rock and surface movement deformation in multi-coal seam strip mining were studied systematically.The reason for the influencing law is analyzed，and the function relations were given.The calculating method of ultimate coal dislocation coefficient in multi-coal seam strip mining was pointed out.The research results have important theoretical and practical significance for the coal face design of multi-coal seam strip mining.

FLAC3D，Numerical simulation，Multi-coal seam strip mining，Coal dislocation coefficient

2015-09-14

“十二五”國(guó)家科技支撐計(jì)劃項(xiàng)目(編號(hào):2012BAB13B03),國(guó)家自然科學(xué)基金青年基金項(xiàng)目(編號(hào)：41104011),江蘇省高校優(yōu)勢(shì)學(xué)科建設(shè)工程項(xiàng)目(編號(hào):SZBF2011-6-B35),江蘇省資源環(huán)境信息工程重點(diǎn)實(shí)驗(yàn)室基金項(xiàng)目(編號(hào):JS201309)。

王 冰(1990—)，男，碩士研究生。通訊作者 郭廣禮(1965—)，男，教授，博士研究生導(dǎo)師。

TD325

A

1001-1250(2015)-11-153-05