Three-dimensional coal seam modeling of fully mechanized working face based on transparent geology
-
摘要: 基于透明地质的三维煤层建模方法是间接解决煤岩识别难题的有效途径。现有三维煤层建模方法的研究大多集中于对空间三维实体的表达,对开采过程中煤层顶底板动态变化的过程缺乏研究,对于复杂地质条件的煤层顶底板高程的预测精度不高,难以满足采煤实际需求。针对上述问题,提出了一种基于透明地质的综采工作面三维煤层建模方法。基于进回风巷地质数据、钻孔测量数据、工作面切眼数据及利用三维地震再解释技术、槽波地震勘探技术与无线电磁波透视技术获得的煤层地质数据,应用离散平滑插值(DSI)算法预测煤层顶底板高程,构建综采工作面静态三维煤层模型。为提高工作面静态三维煤层模型精度,通过切眼开采新揭露的地质信息和DSI算法对其进行动态更新,获得更为精确的工作面动态三维煤层模型,基于更新后的三维煤层模型动态规划采煤机截割曲线,指导采煤机进行自动调高控制,从而实现自适应割煤。将该方法应用于黄陵一号煤矿810综采工作面,结果表明:DSI算法对煤层顶底板高程预测效果优于克里金插值算法和样条函数插值算法,插值平均绝对误差为0.015 5 m;每截割5 m对三维煤层模型更新1次,煤层顶底板高程预测误差≤6.3 cm,满足采煤机截割轨迹精确规划要求。Abstract: The three-dimensional coal seam modeling method based on transparent geology is an effective way to indirectly solve the problem of coal rock identification. Most of the existing three-dimensional coal seam modeling methods focus on the expression of spatial three-dimensional entities. There is a lack of research on the dynamic change of coal seam roof and floor in the mining process. And the prediction precision of coal seam roof and floor elevation under complex geological conditions is not high, which is difficult to meet the actual needs of coal mining. In order to solve the above problems, this paper proposes a three-dimensional coal seam modeling method of fully mechanized working face based on transparent geology. Based on the geological data of air inlet and return roadway, borehole measurement data, open-off cut data of working face and the coal seam geological data obtained by using three-dimensional seismic re-interpretation technology, in-seam seismic exploration technology and wireless electromagnetic wave perspective technology, the discrete smooth interpolation (DSI) algorithm is applied to predict the elevation of coal seam roof and floor. And the static three-dimensional coal seam model of fully mechanized working face is constructed. In order to improve the precision of the static three-dimensional coal seam model of the working face, the geological information newly revealed by open-off cut and DSI algorithm are used to dynamically update the model to obtain a more accurate dynamic three-dimensional coal seam model of the working face. Based on the updated three-dimensional coal seam model, the cutting curve of the shearer is dynamically planned to guide the shearer to automatically adjust height so as to achieve adaptive coal cutting. The method is applied to 810 fully mechanized working face of Huangling No.1 Coal Mine, the results show that the DSI algorithm is better than Kriging interpolation algorithm and spline function interpolation algorithm in the prediction of coal seam roof and floor elevation. The mean absolute error of interpolation is 0.015 5 m. The three-dimensional coal seam model is updated once every 5 m of cutting, and the elevation prediction error of coal seam roof and floor is ≤ 6.3 cm, which meets the requirements for precise planning of the cutting track of the shearer.
-
0. 引言
钢丝绳具有弯曲性能好、承载能力强、质量轻和结构多样性等优点,被广泛应用于矿井牵引和承载工作中[1]。钢丝绳作为牵引和承载设备的关键构件,在长期作业过程中会受到腐蚀、冲击等影响,再加上处于承力、弯曲等工作状态,会出现断丝、锈蚀等损伤[2-3]。断丝是在役钢丝绳常见的损伤类型,钢丝绳产生断丝后,剩余强度降低,削弱其承载能力,严重时可导致钢丝绳断裂,出现极大的安全隐患[4-6]。为确保钢丝绳安全服役,各行业制定了相应的钢丝绳报废定量标准,其中钢丝绳断丝数量是评判钢丝绳是否报废的主要指标[7-8]。漏磁检测因其原理简单、实施成本较低、检测效果好,成为目前钢丝绳无损检测最主要的方法之一[9-11]。现场漏磁检测时将钢丝绳置于检测仪器中,通过移动钢丝绳或检测仪器进行扫描检测[12],分析检测到的漏磁场,从而了解钢丝绳损伤状态,及时更换达到报废标准的钢丝绳,确保钢丝绳安全运行。
为建立断丝与漏磁场的定量关系,需通过理论计算方法分析钢丝绳断丝漏磁场分布规律。理论计算钢丝绳断丝漏磁场的方法分为解析法和数值法。解析法通过简化钢丝损伤模型建立等效磁偶极子模型,该模型最早用来模拟点状缺陷和裂纹的漏磁场[13-15]。K. Seleznyova等[16]建立了3种不同的磁偶极子模型,实现了磁偶极子模型的磁力线可视化。窦连城等[17]采用二维磁荷模型对钢丝绳断丝漏磁场进行分析,得出钢丝绳断丝轴向漏磁场计算公式与磁场分布形态。解析法能够精确地计算出漏磁场,但在计算钢丝绳绳股漏磁场时存在困难。数值法是利用麦克斯韦方程组分析漏磁场参数变化,且利用有限元分析钢丝绳缺陷漏磁场的数值法研究较为完善。李登蓬[18]使用有限元仿真软件建立了二维静态磁场仿真模型,对磁回路的磁化状态进行分析,验证了钢丝绳的磁化饱和状态。郭永亮[19]通过建立钢丝绳二维动态磁场仿真模型,研究钢丝绳绳股漏磁场,证明了用直线运动对钢丝绳绳股漏磁场进行仿真的方法是可行的。朱良[20]使用三维静态磁场仿真模型分析了断口参数变化对漏磁场的影响。目前钢丝绳断丝漏磁场分布规律的仿真研究主要使用有限元静态磁场仿真模型,该模型中钢丝绳与损伤检测仪器是相对静止的,而现场钢丝绳损伤检测过程中钢丝绳与损伤检测仪器之间有相对运动,因此静态磁场仿真漏磁场和实际检测获得的漏磁场存在偏差。
为准确仿真钢丝绳与损伤检测仪器相对运动状态下的漏磁场,本文使用Ansoft Maxwell电磁仿真软件,基于漏磁无损检测原理,建立钢丝绳三维动态磁场仿真模型,分析不同断口宽度、断丝数量和提离值下断丝轴向漏磁场分布规律。
1. 钢丝绳漏磁无损检测原理和漏磁场特征定义
钢丝绳漏磁无损检测原理如图1所示。使用永磁体将钢丝绳磁化至饱和,钢丝绳与气隙、永磁体、衔铁形成磁回路,钢丝绳出现断丝等损伤时,断口位置的磁阻会增大,有一部分磁通量泄漏至空气中形成漏磁场。使用霍尔元件等磁敏元件检测漏磁场,实现钢丝绳损伤状态的识别[21]。
钢丝绳轴向漏磁场磁感应强度如图2所示。中间位置为钢丝绳断丝轴向漏磁场,Pz为钢丝绳断丝轴向漏磁场峰峰值,两侧类正弦波为钢丝绳绳股漏磁场。
2. 基于三维动态磁场仿真模型的钢丝绳断丝漏磁场分析
Ansoft Maxwell软件可通过设置模型的尺寸、结构和材料电磁属性,仿真分析二维/三维模型的静态磁场或动态磁场,根据钢丝绳漏磁无损检测原理,提取待检测位置仿真漏磁场[22]。
Ansoft Maxwell依据麦克斯韦方程组进行仿真计算:
$$ \left\{\begin{array}{l} \nabla\cdot {\boldsymbol{D}}={\rho }_{0}\\ \nabla \times {\boldsymbol{E}}=-\dfrac{\partial {\boldsymbol{B}}}{\partial t}\\ \nabla\cdot {\boldsymbol{B}}=0\\ \nabla \times {\boldsymbol{H}}={{\boldsymbol{j}}}_{0}+\dfrac{\partial {\boldsymbol{D}}}{\partial t}\end{array}\right. $$ (1) 式中:$ \nabla $为矢量微分算子;D为电位移;$ {\rho }_{0} $为自由电荷密度;E为电场强度;B为磁感应强度;H为磁场强度;$ {{\boldsymbol{j}}}_{0} $为传导电流密度;$ t $为时间。
静态磁场仿真中模型是静止的,只能仿真模型单一时间节点的磁场状态;动态磁场仿真可以通过设置部分模型的运动速度和方向仿真模型运动时间内任意时间节点的瞬时磁场状态。动态磁场仿真相较于静态磁场仿真,其仿真结果含有断丝漏磁场和钢丝绳绳股漏磁场,更接近实际情况下钢丝绳漏磁信号,因此使用动态磁场仿真钢丝绳断丝漏磁场。
三维动态磁场仿真模型包括钢丝绳模型和励磁结构模型。根据6×37+FC的钢丝绳结构建立简化钢丝绳模型,模型直径为43 mm,断丝设置在钢丝绳中间位置表面,使用与钢丝绳绳股旋向相同的矩形凹槽表示断丝情况,不同断丝参数设置见表1。
表 1 不同断丝参数设置Table 1. Settings for different broken strands parameters断丝数量/根 断丝截面长度/mm 断丝深度/mm 1 2 2 2 4 2 3 6 2 4 8 2 励磁结构模型由永磁体和衔铁组成,如图3所示,采用1对径向充磁的环形永磁体进行励磁,永磁体材料为钕铁硼,磁铁之间的导磁部件为衔铁,材料为工业纯铁,其余部分设置为空气。
励磁结构模型轴线与钢丝绳模型轴线重合,励磁结构模型保持不动,钢丝绳模型沿轴线移动,从运动初始位置开始移动,经过励磁结构模型中心位置,最后移动到运动结束位置,总运动时间为50 s,初始位置和结束位置关于励磁结构模型中心对称,如图4所示。
基于漏磁无损检测原理,为使磁化钢丝绳断丝后的漏磁场易于检测与分辨,需要将钢丝绳磁化至饱和状态且磁化均匀。绘制钢丝绳磁通量密度云图(图5),判断钢丝绳的磁化状态和磁通分布。
从图5可看出,钢丝绳断丝周边已被均匀磁化,断丝断口端面处相对于断口周围位置的磁通量密度偏小,这是由于断丝断口断面处的空气磁阻远大于钢丝绳,阻碍了磁力线的通过。
以下分析不同因素对断丝漏磁场的影响。
1) 断口宽度对断丝漏磁场的影响分析。设置断丝断口宽度分别为3,5,7,10,13,15 mm,提取断丝上方5 mm处轴向漏磁场进行分析。通过三维动态磁场仿真模型得到的不同断口宽度下断丝轴向漏磁场磁感应强度及漏磁场峰峰值分别如图6和图7所示。可看出三维动态磁场仿真计算的漏磁场中包含断丝漏磁场和钢丝绳绳股漏磁场;在4种断丝数量情况下,随着断口宽度增加,断丝轴向漏磁场峰峰值均呈先增大后减小的趋势,且峰峰值在断口宽度为7,10 mm处开始减小。
2) 不同断丝数量对断丝漏磁场的影响分析。设置断丝数量分别为1,2,3,4根,提取断丝上方5 mm处轴向漏磁场进行分析。通过三维动态磁场仿真模型得到的不同断丝数量下断丝轴向漏磁场磁感应强度及漏磁场峰峰值分别如图8和图9所示。可看出在6种断口宽度情况下,断丝轴向漏磁场峰峰值均随着断丝数量的增加而增大。
3) 提离值对断丝漏磁场的影响分析。定义霍尔元件与钢丝绳表面之间的距离为提离值。设置钢丝绳断丝数量为3根、断口宽度为10 mm,提取断丝上方1~10 mm处轴向漏磁场进行分析。通过三维动态磁场仿真模型得到的不同提离值下断丝轴向漏磁场磁感应强度及漏磁场峰峰值分别如图10和图11所示。可看出断丝轴向漏磁场峰峰值随着提离值的增加而减小。
3. 基于三维磁偶极子模型的钢丝绳断丝漏磁场验证分析
为验证三维动态磁场仿真模型仿真结果的正确性,构建常用于钢丝绳断丝漏磁场分析的三维磁偶极子模型,如图12所示。
设定模型轴线中点o为原点,模型轴线方向为y轴,模型横截面内x轴和z轴正交,建立笛卡尔坐标系。用磁荷+Q和−Q产生的磁场代替钢丝绳两侧断面磁场,这2个磁场关于原点o对称。断面到原点o的距离为l,断口宽度为2l,磁荷+Q和−Q的坐标分别为($ {x}_{1},l,{z}_{1} $)和($ {x}_{1},-l,{z}_{1} $)。P为空间任一点,其坐标为($ {x}_{P} $,$ {y}_{P} $,$ {z}_{P} $),磁荷+Q和−Q到P点的磁场矢量分别为$ {{\boldsymbol{r}}_{1}} $和$ {{{\boldsymbol{r}}}_{2}} $,根据磁场库仑定律,+Q在P点处产生的磁场强度$ {{\boldsymbol{H}}}_{+Q} $和−Q在P点处产生的磁场强度$ {{\boldsymbol{H}}}_{-Q} $分别为
$$ {{\boldsymbol{H}}}_{+Q}=\dfrac{{q}_{1}{{\boldsymbol{ r}}}_{1}}{4{\text{π}} {\mu }_{0}{\left|{{\boldsymbol{ r}}}_{1}\right|}^{3}} $$ (2) $${{\boldsymbol{H}}}_{-Q}=\dfrac{{q}_{1}{{\boldsymbol{r}}}_{2}}{4{\text{π}} {\mu }_{0}{\left|{{\boldsymbol{r}}}_{2}\right|}^{3}} $$ (3) 式中:$ {q}_{1} $为两磁荷的磁荷量;$ {\mu }_{0} $为真空环境的磁导率。
根据磁场叠加原理,P点处磁场强度为
$$ {{\boldsymbol{H}}}_{{{P}}}={{\boldsymbol{H}}}_{+Q}+{{\boldsymbol{H}}}_{-Q}=\frac{{q}_{1}{ {\boldsymbol{r}}}_{1}}{4{\text{π}} {\mu }_{0}{\left|{{\boldsymbol{r}}}_{1}\right|}^{3}}+\frac{{q}_{1}{{\boldsymbol{r}}}_{2}}{4{\text{π}} {\mu }_{0}{\left|{{\boldsymbol{r}}}_{2}\right|}^{3}} $$ (4) 真空环境的磁导率定义为磁感应强度B与磁场强度H的比值,则
$$ {\boldsymbol{B}}={\mu }_{0}{\boldsymbol{H}} $$ (5) 结合式(4)和式(5),得到P点处磁感应强度:
$$ \boldsymbol{B}_{\mathrm{P}}=\boldsymbol{B}_{+Q}+\boldsymbol{B}_{-Q}=\frac{q_1\boldsymbol{r}_1}{4\text{π}\left|\boldsymbol{r}_1\right|^3}+\frac{q_1\boldsymbol{r}_2}{4\text{π}\left|\boldsymbol{r}_2\right|^3} $$ (6) 单根断丝时,P点处轴向磁感应强度分量为
$$ \begin{split} {B}_{P\mathrm{s}}=&\frac{{q}_{1}}{4{\text{π}} }\left\{\frac{{y}_{P}+l}{\left[\right({x}_{P}-{x}_{1}{)}^{2}+({y}_{P}+l{)}^{2}+({z}_{P}-{z}_{1}{)}^{2}{]}^{\tfrac{3}{2}}}\right.-\\ &\left.\frac{{y}_{h}-l}{\left[\right({x}_{P}-{x}_{1}{)}^{2}+({y}_{P}-l{)}^{2}+({z}_{P}-{z}_{1}{)}^{2}{]}^{\tfrac{3}{2}}}\right\} \end{split} $$ (7) 多根断丝时,P点处轴向磁感应强度分量可看作多对磁荷磁场的叠加,即
$$\begin{split} {B}_{P\mathrm{m}}=&\sum _{i=1}^{n}\frac{{q}_{i}}{4{\text{π}} }\left\{\frac{{y}_{P}+l}{\left[\right({x}_{P}-{x}_{i}{)}^{2}+({y}_{P}+l{)}^{2}+({z}_{P}-{z}_{i}{)}^{2}{]}^{\tfrac{3}{2}}}\right.-\\ &\left.\frac{{y}_{P}-l}{\left[\right({x}_{P}-{x}_{i}{)}^{2}+({y}_{P}-l{)}^{2}+({z}_{P}-{z}_{i}{)}^{2}{]}^{\tfrac{3}{2}}}\right\} \end{split} $$ (8) 式中:$ {x}_{i} $和$ {z}_{i} $为第i处断丝时磁荷在x轴和z轴的坐标;n为断丝数量。
通过三维磁偶极子模型得到的钢丝绳断丝轴向漏磁场磁感应强度如图13所示。
1) 不同断口宽度下断丝漏磁场验证分析。设置断丝断口宽度为1~20 mm,变化步长为1 mm,提取断丝上方5 mm处轴向漏磁场进行分析。通过三维磁偶极子模型得到的不同断口宽度下断丝轴向漏磁场峰峰值如图14所示。可看出在4种断丝情况下,随着断口宽度的增加,通过三维磁偶极子模型计算的断丝轴向漏磁场峰峰值均呈先增大后减小的趋势,且峰峰值在断口宽度为8,9 mm处开始减小。
2) 不同断丝数量下断丝漏磁场验证分析。设置断丝数量分别为1,2,3,4根,提取断丝上方5 mm处轴向漏磁场进行分析。通过三维磁偶极子模型得到的不同断丝数量下断丝轴向漏磁场峰峰值如图15所示。可看出在20种断口宽度情况下,通过三维磁偶极子模型计算的断丝轴向漏磁场峰峰值均随着断丝数量的增加而增大。
3) 不同提离值下断丝漏磁场验证分析。设置钢丝绳断丝数量为3根、断口宽度为10 mm,提取断丝上方1~10 mm处轴向漏磁场进行分析。通过三维磁偶极子模型得到的不同提离值下断丝轴向漏磁场峰峰值如图16所示。可看出通过三维磁偶极子模型计算的断丝轴向漏磁场峰峰值随着提离值的增加而减小。
由上述分析可知,通过三维磁偶极子模型得到的断丝轴向漏磁场峰峰值随断口宽度、断丝数量、提离值变化规律与通过三维动态磁场仿真模型得到的变化规律相同,验证了三维动态磁场仿真模型的正确性。
4. 结论
1) 三维动态磁场仿真模型可以模拟损伤检测仪器与钢丝绳间的相对运动状态,其计算的漏磁场含有断丝漏磁场和钢丝绳绳股漏磁场,更接近实际漏磁场。
2) 断丝轴向漏磁场峰峰值随着断口宽度的增加呈先增大后减小的趋势,随着断丝数量的增加而增大,随着提离值的增加而减小。
-
表 1 3种插值算法MAE的对比
Table 1 Comparison of mean absolute error (MAE) of three interpolation algorithms
m 克里金插值算法 DSI算法 样条函数插值算法 0.022 5 0.015 5 0.231 2 表 2 顶板高程预测误差
Table 2 Roof elevation prediction error
m 采样点 MAE 采样点 MAE 采样点 MAE 1 0.053840 16 0.910000 31 0.501333 2 0.091417 17 0.617000 32 0.562667 3 0.337333 18 0.114767 33 0.497667 4 0.468333 19 0.178380 34 0.431667 5 0.885000 20 0.420667 35 0.250667 6 1.010000 21 0.684667 36 0.108667 7 1.030000 22 0.742000 37 0.103667 8 1.050667 23 0.879333 38 0.155333 9 0.416667 24 1.036000 39 0.003667 10 0.051093 25 0.980000 40 0.004000 11 0.271667 26 0.890333 41 0.153333 12 0.507333 27 0.357183 42 0.286000 13 0.473333 28 0.140887 43 0.362667 14 0.619333 29 0.019613 44 0.552333 15 0.920000 30 0.287333 45 0.229667 -
[1] 王国法,刘峰,庞义辉,等. 煤矿智能化−煤炭工业高质量发展的核心技术支撑[J]. 煤炭学报,2019,44(2):349-357. WANG Guofa,LIU Feng,PANG Yihui,et al. Coal mine intellectualization:the core technology of high quality development[J]. Journal of China Coal Society,2019,44(2):349-357.
[2] 袁亮,张平松. 煤炭精准开采地质保障技术的发展现状及展望[J]. 煤炭学报,2019,44(8):2277-2284. YUAN Liang,ZHANG Pingsong. Development status and prospect of geological guarantee technology for precise coal mining[J]. Journal of China Coal Society,2019,44(8):2277-2284.
[3] 张小艳,朱圣凯,杨鑫磊. 采煤工作面煤层三维地质建模[J]. 科学技术与工程,2020,20(10):4049-4055. DOI: 10.3969/j.issn.1671-1815.2020.10.038 ZHANG Xiaoyan,ZHU Shengkai,YANG Xinlei. Three-dimensional geological modeling of coal seam in mining face[J]. Science Technology and Engineering,2020,20(10):4049-4055. DOI: 10.3969/j.issn.1671-1815.2020.10.038
[4] 贾庆仁,车德福,李佳徐,等. 动态精化的煤层三维建模方法[J]. 东北大学学报(自然科学版),2018,39(5):726-730. JIA Qingren,CHE Defu,LI Jiaxu,et al. Three-dimensional modeling method of coal seam with gradual refinement[J]. Journal of Northeastern University (Natural Science),2018,39(5):726-730.
[5] 刘万里,张学亮,王世博. 采煤工作面煤层三维模型构建及动态修正技术[J]. 煤炭学报,2020,45(6):1973-1983. LIU Wanli,ZHANG Xueliang,WANG Shibo. Modeling and dynamic correction technology of 3D coal seam model for coal-mining face[J]. Journal of China Coal Society,2020,45(6):1973-1983.
[6] 肖静. 基于伪点剔除与四域样条插值的三维煤层建模[J]. 科技通报,2014,30(12):166-168. DOI: 10.3969/j.issn.1001-7119.2014.12.056 XIAO Jing. 3D coal seam modeling based on pseudo point elimination and four fields spline interpolation[J]. Bulletin of Science and Technology,2014,30(12):166-168. DOI: 10.3969/j.issn.1001-7119.2014.12.056
[7] 朱德福,邢存恩,乔港介. 基于改进三棱柱体元构建三维煤层模型的实现[J]. 工矿自动化,2014,40(7):22-24. ZHU Defu,XING Cun'en,QIAO Gangjie. Construction of 3D coal seam model based on improved triangular prism element[J]. Industry and Mine Automation,2014,40(7):22-24.
[8] 修春华,车德福,贾国兵. 含复杂地质构造的三维煤层动态建模方法[J]. 矿山测量,2015(6):52-55,59. DOI: 10.3969/j.issn.1001-358X.2015.06.16 XIU Chunhua,CHE Defu,JIA Guobin. Dynamic modeling method of 3D coal seam containing complex geological structure[J]. Mine Surveying,2015(6):52-55,59. DOI: 10.3969/j.issn.1001-358X.2015.06.16
[9] 周为喜,陈玉华,杨永国,等. 基于角点网格的煤层三维建模与可视化研究[J]. 煤田地质与勘探,2016,44(5):53-57. DOI: 10.3969/j.issn.1001-1986.2016.05.010 ZHOU Weixi,CHEN Yuhua,YANG Yongguo,et al. 3D modeling and visualization of coal reservoir based on corner-point grid[J]. Coal Geology & Exploration,2016,44(5):53-57. DOI: 10.3969/j.issn.1001-1986.2016.05.010
[10] 李晓军,胡金虎,朱合华,等. 基于Kriging方法的煤层厚度估计及三维煤层建模[J]. 煤炭学报,2008,33(7):765-769. DOI: 10.3321/j.issn:0253-9993.2008.07.009 LI Xiaojun,HU Jinhu,ZHU Hehua,et al. The estimation of coal thickness based on Kriging technique and 3D coal seam modeling[J]. Journal of China Coal Society,2008,33(7):765-769. DOI: 10.3321/j.issn:0253-9993.2008.07.009
[11] 吴王文. 基于GIS的煤层厚度变化规律及三维建模研究[D]. 徐州: 中国矿业大学, 2016. WU Wangwen. Study on coal seam thickness variation and 3D modeling based on GIS[D]. Xuzhou: China University of Mining and Technology, 2016.
[12] 张龙正. 基于MicroStation的煤田地质三维建模研究[J]. 煤炭工程,2020,52(7):37-40. ZHANG Longzheng. Coalfield geology 3D modeling with MicroStation[J]. Coal Engineering,2020,52(7):37-40.
[13] 李章林,吴冲龙,张夏林,等. 煤炭三维地质建模信息系统的研制及关键技术[J]. 煤炭学报,2011,36(7):1117-1123. LI Zhanglin,WU Chonglong,ZHANG Xialin,et al. Key technologies and development of a 3D coal geological modeling information system[J]. Journal of China Coal Society,2011,36(7):1117-1123.
[14] 刘勇,崔洪庆. 基于单元划分的复杂岩层面三维建模方法研究[J]. 工矿自动化,2017,43(12):99-103. LIU Yong,CUI Hongqing. Research on 3D modeling method of complex rock strata surfaces based on unit division[J]. Industry and Mine Automation,2017,43(12):99-103.
[15] 荆永滨,杜学胜,张瑞林,等. 复杂地质构造煤层三维模型自动构建技术[J]. 辽宁工程技术大学学报(自然科学版),2016,35(3):243-247. JING Yongbin,DU Xuesheng,ZHANG Ruilin,et al. Techniques for automatic 3D modeling of coal seam with complicated geological structure[J]. Journal of Liaoning Technical University (Natural Science Edition),2016,35(3):243-247.