孫學陽 李玲華 李成 苗霖田

摘 要：為了研究工作面的布置方式對采煤沉陷的影響，以安山煤礦為例，采用相似材料模擬實驗、數值模擬等研究手段，在3-1號煤中布置面間煤柱寬度分別為20 m和40 m的2個工作面，在3-1號煤的下方3號煤中布置一個工作面，且使該工作面中心與上層煤留設的煤柱中心重合，構建了雙煤層工作面的2種布置方式。研究結果表明，留設不同寬度的煤柱，采用數值模擬和相似材料模擬得到的覆巖破壞變形特征基本一致。煤層重復采動工作面布置方式使工作面上覆巖層下沉形態(tài)呈現“W”型。工作面留設的煤柱越寬，煤柱上方覆巖下沉范圍顯著減小，對于覆巖的支撐效果越好，越能抵抗覆巖的移動變形，采煤沉陷越不明顯。通過對結果的分析，在布置工作面時，適當的增加煤柱寬度，可以有效減緩采煤沉陷，為科學預防采煤沉陷災害的發(fā)生提供了參考。關鍵詞：雙煤層;相似材料模擬;數值模擬;覆巖破壞;采煤沉陷

中圖分類號：TD 323??????????? 文獻標志碼：A

文章編號：1672-9315（2022）01-0099-08

DOI：10.13800/j.cnki.xakjdxxb.2022.0114開放科學（資源服務）標識碼（OSID）：

The influence of mining layout of double coal

seam faces on coal mining subsidence

SUN Xueyang1，2，3，LI Linghua1，2，3，LI Cheng1，3，MIAO Lintian4

（1.College of Geology and Environment，Xi’an University of Science and Technology，Xi’an 710054，China;

2.Shaanxi Provincial Key Laboratory of Geological Support for Coal Green Exploitation，Xi’an 710054，China;

3.Key Laboratory of Mine Geological Hazards Mechanism and Control，Xi’an 710054，China;

4.Key Laboratory of Coal Resources Exploration and Comprehensive Utilization，Xi’an 710021，China）

Abstract：

Anshan Coal Mine in northern Shaanxi is taken as the research object to explore

the influence of working face layout on coal mining subsidence.

Similar material simulation experiments and numerical simulation and other research methods

are carried out to arrange two working faces with the width of coal pillars

20 m and 40 m respectively.One working face is arranged in the No.3 coal below the No.3-1? coal，and the center of the working face coincides with the center of the coal pillar set in the upper coal to form a two layout methods of double coal seam working face.The research results show that the failure and deformation characteristics of the overburden rock obtained by numerical simulation and similar material simulation are basically the same after leaving coal pillars of different widths.The layout of the repeated mining face of the coal seam causes the overlying strata of the face to present the? shape of “W”.The wider the coal pillars left in the working face，the lower the subsidence range of the overlying rock above the coal pillars，the better the supporting effect of the overlying rock，the more resistant to the movement and deformation of the overlying rock，and the less obvious coal mining subsidence.

An analysis of the results shows that a proper increase in

the width of the coal pillars can effectively slow down the coal mining subsidence，a reference for scientific prevention of the occurrence of coal mining subsidence disasters.Key words：double coal seams;simulation of similar materials;numerical simulation;overburden failure;coal mining subsidence

0 引 言

采煤沉陷災害嚴重威脅礦區(qū)的安全與協(xié)調發(fā)展，因此科學預防采煤沉陷的災害，對煤礦安全開采以及有效保護礦區(qū)生態(tài)環(huán)境具有重要意義[1-4]。國內外相關專家采用物理模擬、數值模擬及理論分析等方法，研究了煤層開采覆巖移動破壞過程，“三帶”發(fā)育規(guī)律及破碎巖體裂隙分形特征，對有效緩解采煤沉陷災害的發(fā)生奠定了理論基礎[5-14]。中國多數的礦井賦存有2層及以上的可采煤層，多煤層與單煤層采煤沉陷特征顯著不同。實踐證明，多煤層開采對覆巖造成多次復合破壞和擾動[15-19]。潘瑞凱等采用物理相似模擬與PFC 2D軟件研究了淺埋近距離雙厚煤層覆巖裂隙發(fā)育規(guī)律[20];黃遠對雙煤層開采條件下，覆巖裂縫發(fā)育規(guī)律進行了深入研究[21];李闊運用相似材料模擬和數值模擬手段，研究了雙煤層在不同采寬、不同間距條件下開采時覆巖的地表沉陷規(guī)律[22];姜福興等提出雙煤層開采，采高增加會使采空區(qū)上方頂板“活化”，導致原鉸接平衡巖失穩(wěn)[23];李樹清等再現了重復開采條件下，覆巖裂隙發(fā)育發(fā)展的過程[24];李樹剛等提出煤層重復開采時覆巖裂隙發(fā)育呈橢拋帶動態(tài)形態(tài)，并建立了裂隙變化模型[25];張志祥等研究發(fā)現在雙層煤采動過程中，當下層煤工作面推進距離不斷增加，巖體裂隙分形維數經歷了從小到大再到小并穩(wěn)定的變化過程[26];侯恩科等采用PFC 2D軟件研究了淺埋雙煤層開采地表裂縫發(fā)育規(guī)律，認為在雙煤層開采中采高是影響覆巖破壞程度和地面裂縫發(fā)育程度的重要因素之一[27]。相關學者對雙煤層采動的研究大多集中在煤層采動后覆巖裂隙發(fā)育規(guī)律和覆巖移動特征方面，雙煤層工作面開采布置方式對采煤沉陷的影響仍需做進一步的研究?；诖?，筆者以安山煤礦4103工作面雙煤層開采為工程背景，分別建立了雙煤層開采物理相似材料模型和FLAC3D數值模型，分析雙煤層在不同工作面布置條件下開采對采煤沉陷的影響，以期為科學預防采煤沉陷災害的發(fā)生提供參考。

1 相似材料模擬實驗安山煤礦大地構造位于華北地臺鄂爾多斯盆地東緣河東斷褶帶西側。井田內地形切割強烈，大部分緩坡、山梁被新近系紅土和第四系黃土所覆蓋，溝谷、陡坡為基巖出露區(qū)，可采煤層為3-1號煤，厚度3.85～6.90 m，平均5.0 m，是近水平煤層;其次為3號煤，厚度0.87～4.32 m，平均2.5 m

是近水平煤層。3-1號煤為上層煤，3號煤為下層煤。

1.1 確定實驗參數以安山煤礦4103工作面3號煤和3-1號煤為研究對象，進行相似材料模擬實驗。通過在3-1號煤中布置2個工作面，兩工作面中分別留設20 m和40 m的煤柱，在3號煤中布置一個工作面，該工作面中心與上層煤留設的煤柱中心重合，構建了雙煤層工作面的2種布置方式，

研究雙煤層工作面開采布置方式對采煤沉陷的影響。

實驗所用臺架為一種鋼架結構，規(guī)格為2 m×2 m×0.2 m，模擬地層所用材料包括：河沙、石膏、大白粉和白云母片等，結合煤系賦存地區(qū)，地層結構的一般特點，對不同硬度的巖層進行模擬。實驗建立的模型厚度和巖石力學參數見表1。

1.2 實驗過程模型幾何相似比1∶250。本次實驗分別建立模型一和模型二，兩模型均在3-1號煤上布置第1和第2工作面，在3號煤中布置第3工作面，先采上層3-1號煤，再采下層3號煤，在3號煤與3-1號煤正中間布置第一排測點，在3-1號煤以上25 m處布置第二排測點，在地表布置第三排測點。各排測點間距為25 m。其中模型一留設煤柱寬度為20 m，模型二留設煤柱寬度為40 m（圖1）。

煤層在開挖過程中，上覆巖層發(fā)生垮落，且隨著開挖距離的增大，垮落范圍也在不斷擴大。不同工作面發(fā)生初次垮落的時間不同，當模型一3-1號煤第1，2個工作面分別推進至62.5和75 m時，頂板發(fā)生初次垮落，3號煤工作面頂板初次垮落發(fā)生在工作面推進至75 m時（圖2（a））。模型二3-1號煤在2個工作面均推進至75 m時，工作面頂板發(fā)生初次垮落，3號煤工作面頂板初次垮落發(fā)生在其工作面推進至87.5 m時（圖2（b））。模型一3號煤在工作面推進至150 m時，裂隙帶高度為61.25 m，垮落高度為26.25 m（圖2（c））。模型二3號煤工作面推進至145 m時，裂隙帶高度為53.25 m，垮落高度為24.5 m（圖2（d））。

1.3 實驗結果

1.3.1 覆巖破壞特征對比由于工作面間疊置關系和留設煤柱寬度不同，3-1號煤對3號煤所產生的影響程度不同，模型一3號煤對比模型二3號煤首次發(fā)生離層以及初次垮落的時間均有所提前，垮落高度也較大。模型開采過程中巖層破壞對比見表2。

1.3.2 覆巖下沉規(guī)律根據實驗所得數據，繪制不同模型開采后的地表下沉曲線如圖3所示。

分析地表位移變化曲線可以得出，模型一中地表下沉位移值在4號、6號監(jiān)測點處較大，對應工作面中心處，5號點位于煤柱上方，此處下沉值最小，并分別以4號、6號監(jiān)測點為中心，下沉值向兩側逐漸減小。模型二分別在13號、18號處，地表下沉值最大，由于留設煤柱的影響，位于煤柱上方的16號測點下沉值最小，且分別以13號、18號監(jiān)測點為中心，下沉值向兩側逐漸減小（圖3）。模型一與模型二地表下沉曲線形態(tài)都呈“W”型，但對比發(fā)現其形態(tài)又有所差異，模型一留設煤柱寬度較小，其地表下沉曲線顯示覆巖中心未下沉的區(qū)域較窄，而模型二留設煤柱寬度相較于模型一較大，其地表下沉曲線所顯示覆巖中心處未下沉的區(qū)域較寬。由此可知，在煤層開采過程中留設煤柱寬度大小對地表下沉量有一定影響，具體表現為在開采方式、程度等均相同的條件下，工作面留設煤柱寬度越大，其抵抗覆巖變形能力越強，可以承受的覆巖壓力范圍越大，地表下沉值越小。

2 雙層煤開采數值模擬煤層在開采過程中，由于圍巖失去了原有的力學平衡，在自重作用下，逐漸出現彎曲、張裂與冒落等現象。采用FLAC3D軟件，模擬不同工作面煤層開采過程中圍巖的應力應變特征，與物理相似材料模擬實驗結果進行對比論證，分析其開采過程中擾動破壞機理。

2.1 模型建立根據安山煤礦4103工作面的實際地質情況，建立有限元三維數值模型，模型上覆巖層物理力學參數、工作面布置情況與其相似材料模擬一致。

2.2 數值模擬結果

2.2.1 塑性區(qū)發(fā)育高度從圖4（a）中可以看出，兩模型在開采3-1號煤過程中，采空區(qū)上方頂板均首先出現剪切破壞，塑性區(qū)整體破壞程度不大，模型二相比于模型一其沿x方向所產生的剪切破壞范圍變大，而沿z方向的剪切破壞范圍有所減小。隨著3號煤工作面開始推進，兩模型采空區(qū)兩端邊緣處開始出現拉張破壞，且采空區(qū)上方的剪切破壞與拉張破壞均顯著發(fā)育。由于留設煤柱寬度增大，模型二相較于模型一其覆巖下沉量變小且各應力破壞范圍減小。將煤層開采結束后的最大主應力圖與塑性區(qū)形態(tài)圖進行對比得出，模型一裂隙帶高度為65.5 m，垮落帶高度為29.7 m;模型二裂隙帶高度為55.25 m，垮落帶高度為28.6 m。

2.2.2 覆巖下沉規(guī)律在煤層開采過程中兩模型頂板上覆巖層下沉量均在不斷增大，但模型二相比于模型一下沉影響范圍較小;當煤層開挖至靠近留設煤柱時，兩模型覆巖下沉值大小幾乎相同，隨著工作面繼續(xù)向前推進，由于留設煤柱寬度的不同，模型覆巖下沉曲線開始出現明顯的不同，采空區(qū)上方巖層下沉值都有不同程度的增大，在左右兩側煤柱以及中間留設煤柱附近，覆巖下沉值增量相對較小，可得出覆巖寬度對覆巖穩(wěn)定性起重要影響。

2.2.3 主應力分布特征煤層開采過程中，采空區(qū)覆巖主應力以拉應力為主，形態(tài)呈現出中間高而兩側低的特征（圖5紅色區(qū)域）。隨著煤層采動繼續(xù)進行，覆巖垮落、卸荷所產生的應力影響逐漸增大，上覆巖層的拉張應力逐漸向煤柱處集中，模型一3-1號煤開采完成時，拉張應力峰值出現在煤柱中心上方，大小為5.25 MPa，模型二應力峰值也出現于留設煤柱上覆巖層，應力峰值達5.16 MPa。對比兩模型3-1號煤開采完成時的應力分布現象可知，留設20 m煤柱開采的模型一，相較于留設煤柱40 m開采的模型二應力集中現象更明顯，在圖中表現為煤柱上方應力集中區(qū)內，呈紅色的范圍更大（圖5）。

3號煤開采完成后，3-1號煤采空區(qū)之上的應力集中現象基本消失，最大主應力集中于3號煤上方，此時模型一與模型二的應力峰值區(qū)均出現在3號煤采空區(qū)中心上方覆巖中（圖6紅色區(qū)域），模型一應力峰值躍升至9.88 MPa，模型二的應力峰值躍升至9.00 MPa。此時對比兩模型的應力分布現象可得，模型一相較于模型二，應力作用現象更明顯，則模型一中3號煤采空區(qū)上方的應力作用相對于模型二較強，最大主應力集中范圍更大，圖中紅色區(qū)域更大，即留設20 m煤柱的開采方式下，煤層覆巖的應力集中現象更明顯。

3 對比分析由數值模擬各巖層應力變化情況可知，留煤柱開采會影響工作面覆巖拉應力的分布，3-1號煤開采完成時，由于煤層中心留設煤柱的支撐作用，工作面覆巖拉應力分布呈現“W”型（圖5），對應相似材料實驗中煤層頂板覆巖下沉曲線形態(tài)呈現“W”型（圖3）。模型一、模型二中的3號煤、3-1號煤全部開采完成后，兩者應力集中現象皆基本趨于穩(wěn)定，對應相似材料模擬實驗開采完成時煤層覆巖的垮落現象，即采空區(qū)內部空間已經被基本壓實。其中，模型一應力作用相比于模型二較大，在數值模擬中表現為采空區(qū)上方應力集中區(qū)內紅色區(qū)域范圍的擴大。相似材料模擬試驗結果顯示，在留設20 m煤柱時，覆巖裂隙帶發(fā)育高度為61.25 m，垮落帶高度為26.25 m;留設40 m煤柱時，裂隙帶高度為53.25 m，垮落帶高度為24.5 m。數值模擬結果顯示，在留設煤柱為20 m時，裂隙帶高度為65.5 m，煤層覆巖垮落帶高度為29.7 m;留設40 m煤柱時，裂隙帶高度為55.25 m，垮落帶高度為28.6 m。通過兩種方法得到的垮落帶高度與裂隙帶高度結果基本一致。雙煤層在留煤柱開采時，煤柱留設寬度越小，兩煤層的疊置區(qū)域越大，壓力差越大，煤層開采后覆巖的破壞范圍越大，覆巖下沉值就越大。因此，在布置工作面時，適當增大留設煤柱的寬度，可以減緩覆巖移動破壞范圍和破壞程度。

4 結 論

1）在雙煤層開采過程中，隨著留設煤柱寬度的增大，煤層頂板上覆巖層下沉量減小，且下沉曲線形態(tài)呈現出“W”型，垮落帶高度與裂隙帶高度也隨之減小。

2）數值模擬過程表明，在煤層開采過程中，工作面上覆巖出現應力變化，主要以拉張應力為主，使得工作面上覆巖出現移動變形，隨著留設煤柱寬度的增加，煤柱上方的應力持續(xù)減小，采空區(qū)正上方應力值，隨著煤柱寬度的增加而減小。

參考文獻（References）：

[1]

湯伏全，黃韓，孫學陽，等.黃土溝壑區(qū)開采沉陷對地形因子的影響研究[J].干旱區(qū)資源與環(huán)境，2016，30（5）：124-128.

TANG Fuquan，HUANG Han，SUN Xueyang，et al.The influence of mining subsidence on terrain factors in gully region of loess[J].Journal of Arid Land Resources and Environment，2016，30（5）：124-128.

[2]蘆家欣，湯伏全，趙軍儀，等.黃土礦區(qū)開采沉陷與地表損害研究述評[J].西安科技大學學報，2019，39（5）：859-866.

LU Jiaxin，TANG Fuquan，ZHAO Junyi，et al.A review of research on mining subsidence and surface damage in loess mining area[J].Journal of Xi’an University of Science and Technology，2019，39（5）：859-866.

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

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

[4]湯伏全，原濤，汪桂生.渭北礦區(qū)厚黃土層采動變形數值模擬研究[J].西安科技大學學報，2011，31（1）：53-58.

TANG Fuquan，YUAN Tao，WANG Guisheng. Numerical simulation study on mining deformation of thick loess layer in Weibei mining area[J].Journal of Xi’an University of Science and Technology，2011，31（1）：53-58.

[5]劉秀英，張永波.采空區(qū)覆巖移動規(guī)律的相似模擬實驗研究[J].太原理工大學學報.2004，35（1）：29-31.

LIU Xiuying，ZHANG Yongbo.Simulaing experiment research on the law of movement of the mined rock[J].Journal of Taiyuan University of Technology，2004，35（1）：29-31.

[6]王創(chuàng)業(yè)，司建鋒，杜小雅，等.基于相似模擬實驗覆巖移動規(guī)律研究[J].煤炭技術.2017，36（1）：61-63.

WANG Chuangye，SI Jianfeng，DU Xiaoya，et al.Research on strata movement law based on similar simulation test[J].Coal Technology，2017，36（1）：61-63.

[7]王永國，王明，許蓬.巴彥高勒煤礦3-1煤層頂板垮落裂縫帶發(fā)育特征[J].煤田地質與勘探，2019，47（增1）：37-42.

WANG Yongguo，WANG Ming，XU Peng.Characteristics of collapsed fractured zone development of No.3-1seam roof in Bayangaoler coal mine[J].Coal Geology & Exploration，2019，47（S1）：37-42.

[8]孫學陽，楊旭，李鵬強，等.煤系地層與滑坡傾向相同條件下開采對滑坡擾動的影響[J].煤礦安全.2016，47（10）：36-39.

SUN Xueyang，YANG Xu，LI Pengqiang，et al.Influence of mining on landslide disturbance under the same inclination of coal measure strata and landslide[J].Safety in Coal Mines，2016，47（10）：36-39.

[9]張培河，張齊，孫學陽，等.煤炭開采覆巖移動導水裂隙帶發(fā)育高度相似材料模擬實驗研究[J].中國煤炭地質，2019，31（10）：49-52.

ZHANG Peihe，ZHANG Qi，SUN Xueyang，et al.Simulation experiment research on the development of highly similar materials in the moving water-conducting fissure zone of coal mining overlying strata[J].China Coal Geology，2019，31（10）：49-52.

[10]湯伏全.試驗模型變形測量數碼照相裝置研究[J].測繪通報，2012（3）：36-39.

TANG Fuquan.Study of digital photographic device to observe deformation for test model[J].Bulletin of Surveying and Mapping，2012（3）：36-39.

[11]錢鳴高，石平五.礦山壓力與巖層控制[M].徐州：中國礦業(yè)大學出版社，2010.

[12]錢鳴高，繆協(xié)興，許家林，等.巖層控制的關鍵層理論[M].徐州：中國礦業(yè)大學出版社，2003.

[13]郭文彬，劉俊.富水厚砂層下二分層綜采導水裂隙帶探測研究[J].煤炭工程，2014，46（12）：60-62.

GUO Wenbin，LIU Jun. Study on detection of water flow cracking zone in fully mechanized coal mining face in second slicing seam under watery thick sand stratum[J].Coal Engineering，2014，46（12）：60-62.

[14]余學義，王昭舜，楊云.大采深綜放開采地表移動變形規(guī)律[J].西安科技大學學報，2019，39（4）：555-563.

YU Xueyi，WANG Zhaoshun，YANG Yun.The law of surface movement and deformation in large mining depth fully mechanized caving mining[J].Journal of Xi’an University of Science and Technology，2019，39（4）：555-563.

[15]曹曉毅，劉小平，田延哲.煤層重復采動對水渠穩(wěn)定性及滲漏性影響評價[J].煤田地質與勘探，2018，46（4）：93-98.

CAO Xiaoyi，LIU Xiaoping，TIAN Yanzhe. Evaluation on influence of repeated coal mining on the stability and leakage of irrigation canal[J].Coal Geology & Exploration，2018，46（4）：93-98.

[16]武謀達.彬長礦區(qū)復合煤層聯(lián)合開采區(qū)涌水特征[J].煤田地質與勘探，2019，47（

1）：133-137.

WU Mouda.Analysis on water burst characteristics during united mining of multiple coal seams in Binchang mining area[J].Coal Geology & Exploration，2019，47（

1）：133-137.

[17]孫學陽.簡單地質-開采條件下的采煤沉陷預測程序設計[D].西安：西安科技大學，2005.

SUN Xueyang. Programming on prediction of coal-mining-induced subsidence under simple geological-mining conditions[D].Xi’an：Xi’an University of Science and Technology，2005.

[18]李智學，申小龍，李明培，等.榆神礦區(qū)最上可采煤層賦存規(guī)律及開采危害程度[J].煤田地質與勘探，2019，47（3）：130-139.

LI Zhixue，SHEN Xiaolong，LI Mingpei，et al. Occurrence regularity of uppermost minable coal seams and their harmful level of mining in Yushen mining area[J].Coal Geology & Exploration，2019，47（3）：130-139.

[19]陳鑫.杭來灣煤礦導水裂隙帶發(fā)育高度數值模擬研究[J].中國煤炭，2019，45（3）：55-59，74.

CHEN Xin. Study on numerical simulation of development height of water-conducting fracture zone in Hanglaiwan mine[J].China Coal，2019，45（3）：55-59，74.

[20]潘瑞凱，曹樹剛，李勇，等.淺埋近距離雙厚煤層開采覆巖裂隙發(fā)育規(guī)律[J].煤炭學報，2018，43（8）：2261-2268.

PAN Ruikai，CAO Shugang，LI Yong，et al.Development of overburden fractures for shallow double thick seams mining[J].Journal of China Coal Society，2018，43（8）：2261-2268.

[21]黃遠.雙煤層重復采動下覆巖變形破壞特征及導高發(fā)育研究[D].北京：中國礦業(yè)大學，2014.

HUANG Yuan. Research on the characteristics of deformation and failure of overlying strata and the development of conductive height under repeated mining of double coal seams[D].Beijing：China University of Mining and Technology，2014.

[22]李闊.雙層混合開采條件下覆巖及地表移動變形規(guī)律研究[D].太原：太原理工大學，2011.

LI Kuo. Research on the movement and deformation laws of overlying rock and ground surface under the condition of double-layer mixed mining[D].Taiyuan：Taiyuan University of Technology，2011.

[23]姜福興，姚順利，魏全德，等.重復采動引發(fā)礦震的機理探討及災害控制[J].采礦與安全工程學報，2015，32（3）：349-355.

JIANG Fuxing，YAO Shunli，WEI Quande，et al.Tremor mechanism and disaster control during repeated mining[J].Journal of Mining & Safety Engineering，2015，32（3）：349-355.

[24]李樹清，何學秋，李紹泉，等.煤層群雙重卸壓開采覆巖移動及裂隙動態(tài)演化的實驗研究[J].煤炭學報，2013，38（12）：2146-2152.

LI Shuqing，HE Xueqiu，LI Shaoquan，et al. Experimental research on strata movement and fracture dynamic evolution of double pressure-relief mining in coal seams group[J].Journal of China Coal Society，2013，38（12）：2146-2152.

[25]李樹剛，丁洋，安朝峰，等.近距離煤層重復采動覆巖裂隙形態(tài)及其演化規(guī)律實驗研究[J].采礦與安全工程學報，2016，33（5）：904-910.

LI Shugang，DING Yang，AN Zhaofeng，et al. Experimental research on the shape and dynamic evolution of repeated mining-induced fractures in short-distance coal seams[J].Journal of Mining & Safety Engineering，2016，33（5）：904-910.

[26]張志祥，張永波，趙雪花，等.雙煤層采動巖體裂隙分形特征實驗研究[J].太原理工大學學報.2014，45（3）：403-407.

ZHANG Zhixiang，ZHANG Yongbo，ZHAO Xuefu，et al.Experimental study on fractal characteristics of fractures in two-coal seam mining rock mass[J].Journal of Taiyuan University of Technology.2014，45（3）：403-407.

[27]侯恩科，從通，謝曉深，等.基于顆粒流的淺埋雙煤層斜交開采地表裂縫發(fā)育特征[J].采礦與巖層控制工程學報，2020，2（1）：20-28.

HOU Enke，CONG Tong，XIE Xiaoshen，et al.Ground surface fracture development characteristics of shallow double coal seam staggered mining based on particle flow[J].Journal of Mining and Strata Control Engineering，2020，2（1）：20-28.

1636501186302