淺埋煤層房柱式采空巷道地震物理模型設(shè)計(jì)及制作

王 盼，朱書階，賈 茜，菅本啟，張淼淼，董蕊靜

淺埋煤層房柱式采空巷道地震物理模型設(shè)計(jì)及制作

王 盼1,2，朱書階2，賈 茜2，菅本啟3，張淼淼2，董蕊靜2

(1. 西安科技大學(xué) 地質(zhì)與環(huán)境學(xué)院，陜西 西安 710056；2. 中煤科工集團(tuán)西安研究院有限公司，陜西 西安 710077；3. 中國(guó)石油大學(xué)(北京) 地球物理學(xué)院，北京 102249)

針對(duì)陜北淺埋煤層房柱式采空區(qū)探測(cè)精度不高的問題，制作地震物理模型進(jìn)行淺埋煤層房柱式采空區(qū)物理模擬研究。按照物理模擬幾何尺寸和波阻抗相似比原則設(shè)計(jì)地震物理模型；進(jìn)行大量相似材料配比試驗(yàn)，確定低速黃土層相似材料為環(huán)氧樹脂和硅橡膠(1∶1.2)，泥巖相似材料為環(huán)氧樹脂、硅橡膠和滑石粉(1∶0.2∶0.6)，煤層為環(huán)氧樹脂和硅橡膠(1∶0.4)，泥質(zhì)砂巖相似材料為環(huán)氧樹脂和滑石粉(1∶0.8)，砂巖為環(huán)氧樹脂和滑石粉(1∶1.2)；采用數(shù)控雕刻機(jī)雕刻煤層采空區(qū)及巷道保證高精度，利用澆筑法和粘接法相結(jié)合制作采空區(qū)夾心層及煤層倒扣澆筑的方法解決采空區(qū)巷道內(nèi)留存空氣的模型制作工藝難題，完成了房柱式采空巷道地震物理模型制作。經(jīng)測(cè)量，物理模型的單層形態(tài)測(cè)量精度為0.2 mm、速度相對(duì)誤差小于5 %、密度絕對(duì)誤差為±0.3 g/cm3，滿足模型設(shè)計(jì)要求。

房柱式采空區(qū)；相似材料配比；地震物理模型；模型制作

淺埋煤層一般是以埋深在150 m以內(nèi)、基巖小于50 m、上覆松散砂土為賦存條件特征的煤層，廣泛分布于我國(guó)西北部，其具有可采煤層多、煤層厚、結(jié)構(gòu)簡(jiǎn)單、穩(wěn)定性好等優(yōu)越的開采條件。我國(guó)在煤炭大規(guī)模開采之前多采用房柱式即開采時(shí)預(yù)留煤柱支撐頂板以防地表沉陷[1]，加之小煤窯多年來利用房柱式無序開采，遺留了許多錯(cuò)綜復(fù)雜的采空區(qū)，采空區(qū)沉陷后形成采煤沉陷積水區(qū)，導(dǎo)致耕地減少、建筑物損壞、水土流失等問題[2]，嚴(yán)重威脅周邊民眾的正常生活和大型煤礦安全生產(chǎn)。

多年來，我國(guó)在淺埋煤層房柱式采空區(qū)勘探方面開展了大量的探測(cè)實(shí)踐及數(shù)值模擬研究[3-7]，但淺埋煤層房柱式采空區(qū)尤其是采空巷道的探測(cè)仍然是難題。

從方法研究的角度看，地震物理模擬是一條途徑。地震物理模擬是進(jìn)行理論基礎(chǔ)研究最直接有效的方法之一，可以為復(fù)雜構(gòu)造的地震勘探提供標(biāo)準(zhǔn)、可靠的數(shù)據(jù)，從而對(duì)各種地震成像方法和軟件系統(tǒng)的成像效果和質(zhì)量進(jìn)行測(cè)試和評(píng)價(jià)[8-9]。進(jìn)行地震物理模擬首先要制作物理模型，為此，開展淺埋煤層房柱式采空巷道地震物理模型設(shè)計(jì)和制作研究，以便為開展物理模擬提供物質(zhì)基礎(chǔ)，進(jìn)而為深入認(rèn)識(shí)淺埋房采區(qū)波場(chǎng)響應(yīng)特征，為淺埋房采區(qū)的地震數(shù)據(jù)采集及處理解釋提供基礎(chǔ)服務(wù)。

1 房柱式采空巷道地震物理模型設(shè)計(jì)

根據(jù)陜北某煤礦的采掘工程平面圖、地質(zhì)資料及鉆孔資料設(shè)計(jì)一個(gè)雙煤層采空巷道地震物理模型。野外地質(zhì)體實(shí)際尺寸為1 300 m×1 000 m× 200 m，分為6層，頂層為疏松黃土層，厚度為60 m；泥巖為上組煤頂板，厚度42 m；上組煤為厚煤層，厚度為9 m，其中采空區(qū)高度為5 m；上組煤底板和下組煤頂板為泥質(zhì)砂巖，厚度為30 m；下組煤厚度為5 m，下組煤底板為砂巖，厚度為54 m。

首先，根據(jù)物理模型采集系統(tǒng)和超聲換能器主頻等因素確定模型比例因子，再根據(jù)采空巷道形態(tài)、分布情況設(shè)計(jì)采空區(qū)、巷道的形態(tài)和具體位置。

1.1 相似準(zhǔn)則下比例因子確定

模型設(shè)計(jì)一定要考慮模型與實(shí)際介質(zhì)的相似性，相似性的要求有很多方面，但是現(xiàn)實(shí)的模型材料尚達(dá)不到全部要求，所以在解決不同問題時(shí)，有不同的側(cè)重點(diǎn)。在采空區(qū)探測(cè)中主要用到反射地震勘探，主要研究縱波在不同構(gòu)造形態(tài)的界面和多層界面反射所形成的圖形，因此，對(duì)實(shí)際地質(zhì)構(gòu)造的模擬主要是幾何尺寸和波阻抗兩個(gè)方面。

幾何尺寸的相似常常用視波長(zhǎng)作標(biāo)度，不僅包含模型幾何尺寸的大小，還要考慮震源激發(fā)的脈沖主頻和介質(zhì)傳播中波的速度[10-11]。所以幾何尺寸包含時(shí)間比例、幾何尺寸比例和速度比例，其對(duì)應(yīng)關(guān)系如下：

波長(zhǎng)比例因子為：

即波在模型介質(zhì)中旅行時(shí)間為1 μs，相當(dāng)于在實(shí)際介質(zhì)中的旅行時(shí)間1 ms。

由長(zhǎng)度比例因子(L)和波長(zhǎng)比例因子相等，可計(jì)算模型的尺寸：

即模型尺寸長(zhǎng)度1 m，相當(dāng)于野外工作面長(zhǎng)1 000 m。

根據(jù)時(shí)間比例、幾何尺寸比例和速度比例，最終確定模型的模擬比例為1 000∶1，則物理模型尺寸為1.3 m×1 m×0.2 m。

1.2 模型設(shè)計(jì)

根據(jù)淺埋煤層房柱式采掘的實(shí)際情況設(shè)計(jì)的采空巷道地震物理模型包含巷道、斜交采空巷道、正交采空巷道3類采空區(qū)地質(zhì)體。為最大限度發(fā)揮模型的功效，可在同一模型上完成多項(xiàng)測(cè)試，設(shè)計(jì)上組煤包含斜交和正交采空巷道，可分別分析斜交和正交采空巷道的波場(chǎng)響應(yīng)；設(shè)計(jì)下組煤只包含正交采空巷道，且與上組煤的正交采空巷道空間上有部分重疊，可分析重疊和非重疊區(qū)域的波場(chǎng)特征響應(yīng)。為最大程度降低模型的邊界效應(yīng)，采空區(qū)巷道設(shè)計(jì)在煤層的中間位置。模型設(shè)計(jì)前必須考慮模型材料速度、密度以及模型制作工藝的可實(shí)現(xiàn)性，采集系統(tǒng)的發(fā)射能量、頻率以及換能器尺寸、頻率等參數(shù)能滿足要求。

根據(jù)以上設(shè)計(jì)思路，設(shè)計(jì)模型參數(shù)見表1。

表1 實(shí)際地層和物理模型參數(shù)對(duì)比

模型的剖面如圖1a所示，采空區(qū)和巷道的采高、寬度均為0.005 m，圖1b為上組煤采空區(qū)平面圖，左邊斜交采空巷道，斜交角度為 145°，右邊為正交采空巷道；圖1c下組煤采空區(qū)只設(shè)計(jì)了正交采空巷道。

2 房柱式采空巷道地震物理模型制作

模型制作的關(guān)鍵是材料配比，難點(diǎn)在于采空區(qū)及巷道的制作，因?yàn)槌Ｒ?guī)的模型澆筑方法無法保證采空區(qū)內(nèi)保留空氣，為此采用制作采空區(qū)夾心層的新工藝來解決這一問題。

2.1 模型材料的選擇及配置

材料的選擇是制作模型的關(guān)鍵，特別是在制作多層介質(zhì)模型時(shí)，在選擇材料的時(shí)候要考慮以下幾個(gè)方面[12-15]：

①可模擬性，能滿足模擬實(shí)際地層的彈性參數(shù)，特別在速度、密度和衰減的聲學(xué)特征方面；

②可制作性，在模型制作過程中所選材料要使得界面之間要能很好的粘合、無縫隙；

③制作工藝，要考慮模型的固化性能，便于后期打磨處理。

根據(jù)近幾十年來學(xué)者的研究，樹脂和橡膠類材料或與其他材料混合后，流動(dòng)性、收縮性、固化性都較好，且速度差異大，固化后易于處理。所以在進(jìn)行地震物理模型制作時(shí)?？紤]環(huán)氧樹脂、硅橡膠添加滑石粉等。模型制作之前通過大量的材料配比試驗(yàn)，制作試塊并測(cè)試速度、密度，配比出與設(shè)計(jì)的各層參數(shù)相近的材料。經(jīng)過多次調(diào)試配比得到了采空區(qū)各層材料配比試塊(圖2)，主料為環(huán)氧樹脂(記為A)，配料為硅橡膠和滑石粉分別記為B和C，具體配比比例見表2。

圖1 采空巷道設(shè)計(jì)

從這些試塊中測(cè)得縱波速度和密度的比例為1.5∶1，所以實(shí)際地層轉(zhuǎn)換為模型的比例參數(shù)時(shí)，速度1∶1的情況下，密度按照1∶1.5轉(zhuǎn)換，最終的轉(zhuǎn)換密度和模型設(shè)計(jì)基本一致。經(jīng)計(jì)算，黃土層速度與設(shè)計(jì)的相對(duì)誤差為4.4 %，泥巖層為2.3 %，上組煤為1.1 %，泥質(zhì)砂巖為 0.6%，下組煤為1.9 %，砂巖為4.8 %，速度相對(duì)誤差控制在5 %以內(nèi)；黃土層密度與設(shè)計(jì)的絕對(duì)誤差為0.3 g/cm3，泥巖層為–0.2 g/cm3，上組煤為0.3 g/cm3，泥質(zhì)砂巖為 0.01 g/cm3，下組煤為0.2 g/cm3，砂巖為0.1 g/cm3，密度絕對(duì)誤差控制在±0.3 g/cm3以內(nèi)。

圖2 材料配比試塊

表2 模型參數(shù)配比

在速度測(cè)試方面，由于在工業(yè)和復(fù)合材料中很難找到速度低于1 000 m/s或者高于3 000 m/s的材料，所以第一層黃土層和第六層砂巖誤差較大，其他各層基本吻合。

2.2 模型制作及工藝

整個(gè)模型制作的最大難點(diǎn)是采空區(qū)的制作，一是需要保證采空區(qū)內(nèi)填充為空氣，二是采空區(qū)巷道寬度和高度僅有5 mm。整個(gè)模型制作分為兩部分，第一部分采用澆筑法和粘接法混合制作采空區(qū)夾心層的新技術(shù)方法制作上、下組煤層采空區(qū)，第二部分利用逐層澆筑的方法整體澆筑模型。

2.2.1 采空區(qū)制作

1) 下組煤采空區(qū)

采空區(qū)平面面積為750 mm×350 mm，以下組煤的上下相鄰層作為采空區(qū)夾心層的頂部和底部，分別澆筑20 mm；凝固后底層上面澆筑厚度為5 mm的煤層，利用數(shù)控雕刻機(jī)雕刻煤層采空區(qū)及巷道，如圖3a所示；最后粘貼夾心層的頂部最關(guān)鍵，要把頂部放在下面，薄涂相同材料，煤層倒扣上，再用重物壓制，薄涂和煤層倒扣的方法可以達(dá)到防止材料入侵采空區(qū)和粘貼時(shí)產(chǎn)生氣泡，如圖3b所示。

2) 上組煤采空區(qū)

與下組煤采空區(qū)制作方法相同，煤層厚度9 mm，采空區(qū)高度5 mm。

圖3 上下組煤層采空區(qū)部分制作

2.2.2 模型整體制作

模型整體制作采用從底部逐層澆筑的方法，厚度較大層每次澆筑1～2 cm，否則由于材料混合放熱過快會(huì)出現(xiàn)氣泡、裂縫等現(xiàn)象。

砂巖層制作：澆筑至34 mm高度后放置下組煤采空區(qū)夾心層，繼續(xù)澆筑至砂巖層完成，如圖4所示；下組煤層制作：按照下組煤配比澆筑除采空區(qū)的煤層部分；泥質(zhì)砂巖層制作：作為下組煤的頂板和上組煤的底板，要注意高度計(jì)算準(zhǔn)確；上組煤層和剩余其他層按照相同的方法依次澆筑完成，最后去模具，對(duì)模型進(jìn)行測(cè)量和精細(xì)化打磨，完成模型制作，制作的模型如圖5所示。

模型制作過程中采用高精度三維測(cè)量機(jī)對(duì)每一層進(jìn)行精確測(cè)量，精確控制模型形態(tài)的幾何尺寸，誤差保證在0.2 mm以內(nèi)。

圖4 底板砂巖層和下組煤

圖5 模型實(shí)物

3 結(jié)論

a.經(jīng)過大量材料配比試驗(yàn)，最終采用環(huán)氧樹脂、硅橡膠添加滑石粉按照不同比例得到各層參數(shù)，經(jīng)過測(cè)量，速度相對(duì)誤差小于5 %，密度絕對(duì)誤差為±0.3 g/cm3。

b.采用澆筑法和粘接法制作采空區(qū)夾心層的新技術(shù)方法，解決采空區(qū)巷道及其他巷道內(nèi)留置空氣的模型制作工藝難題。

c. 利用數(shù)控雕刻技術(shù)保證采空區(qū)巷道、巷道的精度，使制作模擬的單層形態(tài)絕對(duì)誤差不大于0.2 mm，實(shí)現(xiàn)了模型制作的高精度。

[1] 劉清洲. 淺埋近距離房柱式采空區(qū)之上綜采圍巖結(jié)構(gòu)穩(wěn)定性研究[D]. 西安：西安科技大學(xué)，2020.

LIU Qingzhou. Research on structural stability of fully mechanized mining surrounding rock in shallow buried and close distance room and pillar goaf[D]. Xi’an：Xi’an University of Science and Technology，2020.

[2] 趙艷玲，李素萃，肖武. 我國(guó)采煤沉陷領(lǐng)域研究現(xiàn)狀可視化分析與展望[J]. 煤炭科學(xué)技術(shù)，2020，48(1)：202–210.

ZHAO Yanling，LI Sucui，XIAO Wu. Visualization analysis and prospect of the current research situation in the field of coal mining subsidence in China[J]. Coal Science and Technology，2020，48(1)：202–210.

[3] 程建遠(yuǎn)，孫洪星，趙慶彪，等. 老窯采空區(qū)的探測(cè)技術(shù)與實(shí)例研究[J]. 煤炭學(xué)報(bào)，2008，33(3)：251–255．

CHENG Jianyuan，SUN Hongxing，ZHAO Qingbiao，et al. The detection technology of excavated region in coal mine and case study[J]. Journal of China Coal Society，2008，33(3)：251–255.

[4] 覃思，程建遠(yuǎn)，胡繼武，等. 煤礦采空區(qū)及巷道的井地聯(lián)合地震超前勘探[J]. 煤炭學(xué)報(bào)，2015，40(3)：636–639.

QIN Si，CHENG Jianyuan，HU Jiwu，et al. Coal-seam-ground-seismic for advance detection of goaf and roadway[J]. Journal of China Coal Society，2015，40(3)：636–639.

[5] 王磊. 房柱式采空區(qū)穩(wěn)定性評(píng)價(jià)與數(shù)值模擬：以陜北府谷縣盛海煤礦為例[J]. 能源與節(jié)能，2020(10)：2–6.

WANG Lei. Stability evaluation and numerical simulation of room and pillar goaf：A case study of Shenghai coal mine in Fugu County，Northern Shaanxi[J]. Energy and Energy Conservation，2020(10)：2–6.

[6] 薛國(guó)強(qiáng)，潘冬明，于景邨. 煤礦采空區(qū)地球物理探測(cè)應(yīng)用綜述[J]. 地球物理學(xué)進(jìn)展，2018，33(5)：2187–2192.

XUE Guoqiang，PAN Dongming，YU Jingcun. Review the applications of geophysical methods for mapping coal-mine voids[J]. Progress in Geophysics，2018，33(5)：2187–2192.

[7] 張昭. 不同年限采空區(qū)下地震勘探效果實(shí)例研究[J/OL]. 煤田地質(zhì)與勘探. https：//kns.cnki.net/kcms/detail/61.1155.P.202 10715.1032.002.html

ZHANG Zhao. Case study on seismic exploration effect under goaf with different years[J/OL]. Coal Geology & Exploration. https：//kns.cnki.net/kcms/detail/61.1155.P.20210715.1032.002. html

[8] WU Mansheng，DI Bangrang，WEI Jianxin，et al. Large-scale complex physical modeling and precision analysis[J]. Applied Geophysics，2014，11(2)，245–251.

[9] 張福宏，黃平，黃開偉，等. 復(fù)雜裂縫地球物理模型制作及地震采集處理研究[J]. 物探與化探，2018，42(1)：87–95.

ZHANG Fuhong，HUANG Ping，HUANG Kaiwei，et al. The construction of complex fracture geophysical model and the gathering and processing of seismic data[J]. Geophysical and Geochemical Exploration，2018，42(1)：87–95.

[10] 趙鴻儒，唐文榜，郭鐵栓. 超聲地震模型試驗(yàn)技術(shù)及應(yīng)用[M].北京：石油工業(yè)出版社，1986.

ZHAO Hongru，TANG Wenbang，GUO Tieshuan. Technology and application of ultrasonic seismic modeling[M]. Beijing：Petroleum Industry Press，1986.

[11] DI Bangrang，XU Xiucang，WEI Jianxin. A seismic modeling analysis of wide and narrow 3D observation systems for channel sand bodies[J]. Applied Geophysics，2008，5(3)：294–300.

[12] 魏建新，狄?guī)妥? 地震物理模型中三維地質(zhì)模型材料特性研究[J]. 石油物探，2006，45(6)：586–590.

WEI Jianxin，DI Bangrang. Properties of materials forming the 3D geological model in seismic physical model[J]. Geophysical Prospecting for Petroleum，2006，45(6)：586–590.

[13] 李智宏，朱海龍，趙群，等. 地震物理模型材料研制與應(yīng)用研究[J]. 地球物理學(xué)進(jìn)展，2009，24(2)：408–417.

LI Zhihong，ZHU Hailong，ZHAO Qun，et al. Study and materialization of new seismic physical model building materials [J]. Progress in Geophysics，2009，24(2)：408–417.

[14] 裴宇翀，楊勤勇，趙群，等. 硅微粉改性新型三維地震物理模型材料特性研究[J]. 地球物理學(xué)進(jìn)展，2016，31(1)：455–460.

PEI Yuchong，YANG Qinyong，ZHAO Qun，et al. Research on silica micro-powders modified 3D seismic physical model materials[J]. Progress in Geophysics，2016，31(1)：455–460.

[15] 董金玉，楊繼紅，楊國(guó)香，等. 基于正交設(shè)計(jì)的模型試驗(yàn)相似材料的配比試驗(yàn)研究[J]. 煤炭學(xué)報(bào)，2012，37(1)：44–49.

DONG Jinyu，YANG Jihong，YANG Guoxiang，et al. Research on similar material proportioning test of model test based on orthogonal design[J]. Journal of China Coal Society，2012，37(1)：44–49.

Design and construction of a seismic physical model of room-pillar goafs in shallow coal seams

WANG Pan1,2, ZHU Shujie2, JIA Qian2, JIAN Benqi3, ZHANG Miaomiao2, DONG Ruijing2

(1. College of Geology and Environment, Xi’an University of Science and Technology, Xi’an 710054, China; 2. Xi’an Research Institute Co. Ltd., China Coal Technology and Engineering Group Corp., Xi’an 710077, China; 3. College of Geophysics, China University of Petroleum, Beijing102249, China)

In view of the problem of low detection accuracy of room-pillar goafs in shallow coal seams in northern Shaanxi, the seismic physical model is designed for physical simulation according to the principle of similarity ratio of geometry size and wave impedance. Through a large number of ratio tests of similar materials, the ratio of epoxy resin and silicone rubber is determined as 1 : 1.2 for similar materials of low velocity loess layer, 1 : 0.2 : 0.6 for epoxy resin, silicone rubber and talc for mudstone, 1 : 0.4 for epoxy resin and silicone rubber for coal seams, 1 : 0.8 for epoxy resin and talc for argillaceous sandstone, and 1 : 1.2 for epoxy resin and talc for sandstone.The numerical control engraving machine is used to carve the coal seam goaf and roadway to ensure high precision. By using the pouring method combined with the bonding method, the sandwich layer of goafs and coal seams is made to solve the problem of air model production in the goaf roadway, and the seismic physical model of the room-pillar goaf roadway is completed. The single-layer shape measurement accuracy of the physical model measures 0.2 mm, the speed relative error less than 5%, and the density absolute error ±0.3 g/cm3, which meets the design requirements of the model.

room-pillar goaf; ratio of similar materials; seismic physical model; model construction

語音講解

P315.8; P631

A

1001-1986(2021)06-0101-06

2021-09-14；

2021-10-27

國(guó)家自然科學(xué)基金項(xiàng)目(41974209)；中煤科工集團(tuán)西安研究院有限公司科技創(chuàng)新基金項(xiàng)目(2019XAYPT02，2020XAYDC02-03)

王盼，1987年生，女，陜西咸陽人，博士研究生，從事物探儀器、地震物理模擬研究工作. E-mail：wangpan@cctegxian.com

王盼，朱書階，賈茜，等. 淺埋煤層房柱式采空巷道地震物理模型設(shè)計(jì)及制作[J]. 煤田地質(zhì)與勘探，2021，49(6)：101–106. doi∶10.3969/j.issn.1001-1986.2021.06.012

WANG Pan，ZHU Shujie，JIA Qian，et al. Design and construction of a seismic physical model of room-pillar goafs in shallow coal seams[J]. Coal Geology & Exploration，2021，49(6)：101–106. doi∶10.3969/j.issn.1001-1986.2021.06.012

移動(dòng)閱讀

(責(zé)任編輯 聶愛蘭)