小工作面開采地表下沉系數探討

成 樞 隋冰冰 王 濤

(山東科技大學測繪科學與工程學院，山東 青島 266590)

小工作面開采地表下沉系數探討

成 樞 隋冰冰 王 濤

(山東科技大學測繪科學與工程學院，山東 青島 266590)

生產實踐中，受采煤技術、地質條件等因素的影響，所開采的工作面往往是煤層走向或傾向不能同時達到充分采動，因此面臨的問題是如何將小工作面開采條件下的概率積分法計算參數換算為充分采動或接近充分采動條件下的計算參數，對于下沉系數的修正又是小工作面概率積分法計算參數的關鍵所在。通過對比某礦地質條件相似的不同工作面下的計算參數后，表明現有的下沉系數修正規律適用于某礦區。

小尺寸工作面 充分采動 下沉系數 計算參數

小尺寸工作面開采地表沉陷預計一直是礦區工作者所關心的問題，開采沉陷的預計方法有概率積分法、典型曲線法、剖面函數法等方法，經過國內外學者幾十年的大量研究，概率積分法預計模型被我國學者所采用，廣泛地應用于礦區的開采沉陷研究當中。但在非充分采動或極不充分采動時,由于巖層結構對地表沉陷有一定的控制作用,偏離概率積分法模型的假設較遠,此時如不對其進行修正。預計結果往往不再適用[1]。下沉系數又是概率積分法模型中的關鍵參數，運用概率積分法求取小尺寸工作面的角值參數時，存在的問題是如何將下沉系數進行修正，通過對某煤礦地質條件相似工作面的下沉系數依據現有規律進行修正，進而求取工作面的概率積分法計算參數，對比計算參數進一步驗證了現有的下沉系數修正規律符合該礦區。

1 小尺寸工作面概況

7269工作面為西九采區首采工作面，工作面標高為-729～ -805 m，工作面走向長度864 m，傾向長度170 m，煤層傾角7～12°，平均9°。回采時間為2010年10月25日至2011年9月13日，平均每天開采推進距離為2.7 m。7267工作面標高-730～-775 m，走向長度930 m，傾向長度160 m，煤層傾角7～12°，平均9°。回采時間為2013年3月1日至2013年10月30日，平均每天開采推進距離為3.7 m。從2個工作面的尺寸可以發現均為小尺寸工作面的開采，7269、7267工作面相對位置以及工作面地表移動觀測站布置如圖1所示。

2 工作面開采地表移動實測情況及資料處理

2.1 工作面實測情況

綜合考慮7269、7267工作面開采的影響范圍以及各測線相應的控制測點的布設問題，根據測站設計，部分觀測線所在的測點個數及位置如下。

圖1 7269、7267工作面地表移動觀測站布置

第1條觀測線：位于鐵路上有A1-A3共3個控制點；位于村莊道路上的編號1～19共19個測點，測線長度484 m；位于村莊內的編號20～29共10個測點，共241 m。

第2條觀測線：位于農田的編號222～229、231～232共10個測點，測線長度298 m；位于村莊內的編號204～209、234～233、210～221共20個測點，測線長度645 m。

第3條觀測線：位于農田的編號115、156～177、179～185共30個測點，測線長度939 m。

第4條觀測線：位于村莊內的編號113～128共16個測點，測線長度481 m。

同時，7269和7267工作面均布設了GPS控制網和一級導線。其中，7269工作面分別在2010年11月、2011年4月、2012年3月進行了3次沉降觀測點平面坐標測量，水準觀測周期為2010年11月30日至2012年4月17日，共31期沉降數據資料。7267工作面分別在2013年1月、2013年9月和2014年5月進行了3次沉降觀測點平面坐標測量，水準觀測觀測周期為2013年2月8日至2014年4月22日，共25期沉降數據資料。

2.2 工作面實測資料數據處理

將各觀測站所針對的主要開采工作面的基本條件和實際觀測資料，利用參數識別的方法，求出各觀測站的地表移動計算參數，并分析出計算參數的可靠性。同時，對各非充分采動條件下的觀測站資料，利用求取的計算參數，模擬計算出本地質開采條件下充分采動或接近充分采動時的地表移動變形。

2.2.1 地表移動計算參數

概率積分法是因其所用的移動和變形預計公式中含有概率積分(或其導數)而得名。由于這種方法的基礎是隨機介質理論，所以又叫隨機介質理論法。

隨機介質理論首先由波蘭學者李特威尼申于上世紀50年代引入巖層移動研究，后由我國學者劉寶琛、廖國華等發展為概率積分法。經過我國開采沉陷工作者幾十年的研究，目前已成為我國較成熟的、應用最為廣泛的預計方法之一[6]。因此基于概率積分法的理論，對7269、7267工作面進行了開采沉陷預計，計算參數如表1所示。

表1 7267、7269工作面概率積分法計算參數

將概率積分法預計結果與實測結果進行對比，表明了概率積分法計算參數的可靠性。7269、7267工作面預計情況與實測情況對比如圖2、圖3所示，曲線1為實測曲線，曲線2為預計曲線。

圖2 7267工作面觀測線下沉曲線擬合

圖3 7269工作面觀測線下沉曲線擬合

2.2.2 概率積分法下沉系數修正

由于7269、7267工作面具有小尺寸采空區的特點，因此在求取工作面的角值參數時需要對下沉系數進行修正，文獻[2]中給出了針對小尺寸采空區下沉系數修正系數規律表，如表2所示。

表2 小尺寸采空區下沉系數修正系數K值

根據表2中提供的數值，可以發現這些離散點的走勢符合指數函數的規律，因此運用曲線擬合的基本理論，對表2進行了公式化表達形式，如公式(1)所示：

(1)

式中，D3為工作面的長度，r為主要影響半徑。

采用指數模型對下沉系數進行修正后，7267、7269工作面充分采動條件下的下沉系數為0.52和0.56,7267和7269工作面充分采動條件下的下沉系數為0.54，從這3個數值可以看出地質情況相似的工作面利用修正模型求取出來的充分采動下沉系數變化不大，進一步表明下沉系數修正系數規律適用于此礦區，并可將下沉系數修正系數指數模型應用于和工作面地質條件相似的其他礦區中。

3 結 論

運用概率積分法對7269、7267小尺寸工作面開采沉陷進行了預計，將概率積分法參數中的下沉系數根據現有修正系數規律換算成充分采動條件對應的下沉系數。對比地質條件相似工作面充分采動條件下的下沉系數后，表明現有下沉系數修正系數規律符合該地質條件下的礦區，同時修正系數指數模型為其他地質條件相似的礦區提供了參考依據。

[1] 查劍鋒,郭廣禮,趙海濤,等.概率積分法修正體系現狀及發展展望[J].金屬礦山,2008(1):15-18. Zha Jianfeng，Guo Guangli，Zhao Haitao，et al.Present situation and prospect of correction system for probability integral method[J].Metal Mine,2008(1):15-18.

[2] 戴仔強,郭廣禮,王卷樂.極小工作面預計參數問題的探討[J].江蘇煤炭,2001(3):15-16. Dai Ziqiang，Guo Guangli，Wang Juanle.Discussion on predicted parameter of minimal working face[J].Jiangsu Coal,2001(3):15-16.

[3] 陳俊杰,鄒友峰,郭文兵.厚松散層下下沉系數與采動程度關系研究[J].采礦與安全工程學報,2012(3):250-254. Chen Junjie，Zou Youfeng，Guo Wenbing.Study on the relationship between subsidence coefficient and mining degree under a thick alluvium stratum[J].Journal of Mining and Safety Engineering,2012(3):250-254.

[4] 蘇仲杰,黃厚旭,趙 松,等.基于數值模擬的充填開采地表下沉系數分析[J].中國地質災害與防治學報,2014(2):98-102. Su Zhongjie，Huang Houxu，Zhao Song，et al.Numerical simulation analysis of surface subsidence coefficient of filling mining[J].The Chinese Journal of Geological Hazard and Control,2014(2):98-102.

[5] 中華人民共和國煤炭工業部.建筑物、水體、鐵路及主要井巷煤柱留設與壓煤開采規程[M].北京:煤炭工業出版社,1986. The Ministry of Coal Industry of the People's Republic of China.Buildings，Water，Railway and Main Well Alley Coal Pillar Left Set Coal Mining with Pressure Regulation[M].Beijing:China Coal Industry Press,1986.

[6] 何國清,楊 倫,凌賡娣,等.礦山開采沉陷學[M].徐州:中國礦業大學出版社,1994. He Guoqing，Yang Lun，Ling Gengdi，et al.Mining Subsidence Science[M].Xuzhou：China University of Mining and Technology Press,1994.

[7] 趙海濱.MATLAB應用大全[M].北京:清華大學出版社，2012. Zhao Haibin.Complete Application of MATLAB[M].Beijing:Tsinghua University Press,2012.

[8] 欒 穎.MATLAB R2013a工程分析與仿真[M].北京:清華大學出版社，2014. Luan Ying.Analysis and Simulation of MATLAB R2013a Project[M].Beijing:Tsinghua University Press,2014.

(責任編輯 石海林)

Discussion on Surface Subsidence Coefficient of Small Mining Face

Cheng Shu Sui Bingbing Wang Tao

(CollegeofGeomatics，ShandongUniversityofScienceandTechnology，Qingdao266590，China)

In production practice，influenced by coal mining technology，geological conditions and other factors，the working face cannot achieve full mining at the same time in coal seam strike or dip.So，the issue faced is how to convert some calculation parameters of probability integral method under the condition of small face mining into parameters under full or close to full mining conditions.Also，correction in factor of subsidence is the key to calculate the parameters by small-face probability integration method.After contrast on calculating parameters of different working faces under similar geological conditions，it is indicated that the correction rules of existing subsidence coefficient is applicable to this mining area.

Small face，Full mining，Factor of subsidence，Calculation parameters

2015-03-03

成 樞(1963—)，男，教授，博士，碩士研究生導師。

TD823

A

1001-1250(2015)-04-071-03