文摘

在浅埋和近距离多个煤层的开采,周围的弱胶结软岩的流变失败下的低煤层巷道的集中压力是非常罕见的。然而,上部煤层的应力在屋顶上传播通过剩余支柱,导致这种情况。以Gaojialiang煤矿位于中国西部的矿区为研究对象,巷道顶板的破坏机理残余煤柱下研究了浅埋和近距离多个缝结合现场监测和数值模拟。此外,建议在这种地质条件提出了巷道支架。结果表明,残余煤柱的工作面煤层低逐渐崩溃在浅埋和近距离多个煤层的开采。煤柱的集中应力转移的进一步增加,这使得低煤层巷道的顶板压力不断增加。此外,巷道顶板的应力也增加进一步由于采空区上方的破碎岩石的旋转,和peek的压力转移到nongoaf区域。结合重大集中应力和弱胶结财产,浅埋围岩显示流变行为和失败。因此,我们必须更加注意的蠕变破坏巷道顶板作用下的残余煤柱甚至在浅埋煤层。

1。介绍

浅埋在中国西部侏罗纪煤田比比皆是,拥有丰富的储备。他们是丰富的可开采的煤层、关闭层间间距(1,2]。矿近距离煤层时,工作面临的上部和下部煤层通常安排在平行(3]。然而,由于断层等地质条件和coal-free区域,交错安排不可避免地出现(4- - - - - -6]。当停止下煤层巷道煤柱下安排部分留下的上部煤层,煤柱载荷可以计算基于假设的煤柱熊装运上覆地层(7,8)在其影响区域根据支流区域理论,和集中压力会导致失败的道路9- - - - - -12]。

为了找出采空区围岩的破坏机理,残余煤柱,或固体煤,一些学者已经做了很多研究矿山压力定律当挖掘浅埋近距离煤层通过现场试验、理论分析和数值模拟。首先,物理和机械性能和组件的研究了浅埋弱胶结的岩石。发现弱胶结岩石的岩性主要由砂质泥岩和细砂岩13,14]。,泥岩包含一个高比例的粘土矿物,这很容易导致解体和高通货膨胀在一定压力(15]。在此基础上,一些学者已经使用传统的压缩试验(16),剪切试验(17),而声发射测试(18),发现故障演化和弱胶结围岩的应力分布有明显的时间效应(19- - - - - -21]。此外,较轻的压力下的岩石稳步[22]。

随后,残余煤柱的集中荷载由上部煤层将向下通过地板层,形成一个stress-concentrated区域在地板上。因此,残留煤柱的集中应力转移规则应该确定(23,24]。使用支流区域理论,煤柱的稳定性分析方法(25)和威尔逊和卡尔的方法(26)开发。通过建立减压开采层的应力增量模型(27,28),它可以获得减压后的煤层开采的地板可分为四个领域方面的地板断裂区(弧),深区域(抛物线)stress-concentrated区域(椭圆形),和原岩应力区。只有当较低的煤巷深地区安排的煤柱底板巷道的稳定性可以更容易得到保证。

最后,有必要分析地板的稳定性在近距离煤层开采的巷道。原因是地板的煤柱巷道顶板的煤层低。基于弹塑性力学和断裂力学在卸荷条件下的动态分布和演变规律在不同应力骨折区域的地板可以实现矿业的影响下(29日- - - - - -31日]。同时,变形的上部煤层开采造成的地面决定(32]。的故障特点和深度间隔不同开采条件下岩石通常是由理论分析、现场监测和数值模拟。数值模拟,因为其简单的常用的长壁开采研究采动岩体响应(33]。

它可以发现成效显著的弱胶结岩石的强度特性的研究和残余煤柱的应力传播。然而,由于在中国西部的浅埋煤层,集中应力下的弱胶结的岩石的流变现象并没有吸引了人们的注意力。流变失败的结果,问题集中荷载的长期作用下巷道煤柱出现(34在生产过程中,带来了安全风险。本文以中国西部的浅埋近距离煤层为研究对象。结合现场监测和数值模拟,研究了巷道顶板的破坏机制下残留煤柱在浅埋近距离多个缝和提供类似地质条件下的巷道支架的基础(35,36]。

2。弱胶结软岩的力学性能

有6层煤Gaojialiang可以开采的煤矿在中国西部。目前,2 - 2和2 - 2煤层上部主要开采。2 - 2上的深度和平均厚度煤层是165和1.71,分别,而中间2 - 2煤层的170和2.73,分别。两个煤层之间的距离是6.65∼11.9 m,算作一个近距离煤层。两者之间的弱胶结的岩石,煤层,由砂质泥岩和粉砂岩。它属于典型的弱胶结软岩层。

为了找出变形和破坏的力学机制下的巷道围岩的残余煤柱和提供的基础类似条件下的巷道支架的设计,组件分析和地层的蠕变试验Gaojialiang层之间进行。

2.1。地层的分析组件

的组件和参数区间测量地层的物理力学和地层的组件如表所示1。间隔层主要由砂质泥岩、砂质泥岩是主要由高岭石和伊利石。都是常见的低硬度的粘土矿物。砂质泥岩的抗拉和抗压强度是11.5 MPa和3.5 MPa,分别,其弹性模量是13.7绩点。可以看出,其强度一般较低,因此巷道的围岩属于典型的弱胶结软岩。

2.2。对弱胶结岩石蠕变试验

正如上面提到的,一些学者进行了浅埋弱胶结的岩石流变测试,发现它会表现出低应力下的流变行为(37- - - - - -39]。因此,获得的砂质泥岩Gaojialiang矿山蠕变试验的选择。一个核心是准备一个标准的样本50 mm×100毫米,graded-loading蠕变测试执行。的故障强度试样在1.5 MPa围岩蠕变压缩14.5 MPa。时程曲线和1.5 MPa围压下的蠕变特性曲线如图所示1。

根据蠕变加载曲线和考虑到轴向,横向,和体积变形特性的岩石,1.5 MPa围压下的应力-应变等时曲线的集群,如图2(一个)- - - - - -2 (c)。三方蠕变等时曲线的集群中,有两个拐点其中一点一个蠕变变形和的符号吗B后强非线性塑性变形。压力点一个可以被视为破坏阈值的生成和B长期强度指数σ_∞岩石的样本。标本的长期强度在1.5 MPa下围压为11.7 MPa。长期强度略低于流变强度。

3所示。地质概况和巷道的破坏

为了进一步验证失败的流变特性和机制弱胶结围岩集中负荷下,我们在Gaojialiang煤矿进行了现场监测。

3.1。地质概况

目前,2 - 2煤层上部已经完成,2 - 2煤层时Gaojialiang煤矿的开采。coal-free影响区域,大部分道路中间的2 - 2煤层开采在这个阶段被安排在2 - 2的采空区上部煤层或残留煤柱。图3介绍了道路和采空区的布局。20314年辅助运输巷道位于中间2 - 2煤层,在交叉安排与多个采空区上覆煤层。下,巷道位于上部煤层采空区或残留煤柱。因此,20314年辅助运输巷道复杂的屋顶压力。特别是残余煤柱下,巷道顶板的下沉等问题,地板的凸起,肋骨剥落,打破螺栓和锚索经常发生。

3.2。现场监测方法

为了观察表面和深部位移下的巷道围岩的残留煤柱在上部煤层,变形和破坏的机理可以显示和围岩的基础支持措施可以提供。站是安排在20314年的辅助运输巷道确定巷道位移的变化规律和支持身体的压力。基于交叉口的特点的工作面临2 - 2和2 - 2煤层上部,对下煤层巷道监测计划制定。1 #安排的结下煤柱采空区,2 #,3 #,5 #煤柱下放置,和4 #安排下固体煤。两个测量点5米间距设置在每一个车站。每个站测试主要包括以下设备:锚(电缆)测力计,井下压力计,钻孔钻电视,等图4描述了监测仪器的安排。

3.3。巷道变形和破坏

3.3.1。巷道的围岩宏观失败

P20314工作面时搬到900米,巷道的围岩44米北四的辅助运输巷道巷道变形和破坏。问题,如巷道顶板的下沉,地板的凸起,肋骨剥落。图5揭示了变形和破坏范围的围岩(54 m)。再次巷道稳定时,巷道顶板在6∼18米右边弧形下沉,而在18∼36米的范围在左边显示了步进式下沉0.6∼0.8米的高度。由于更大的水槽的屋顶,一些螺栓(电缆)破碎和下降和一些单一液压道具被压。图6显示了围岩的变形和破坏形式的巷道。

3.3.2。巷道位移变化

(1)巷道表面位移。巷道表面位移包括巷道顶板的下沉收敛和双方的位移收敛。表面位移观测4 #和5 #车站随时间,如图7(一)和7 (b)。它可以知道沉降和位移收敛先增加,然后随着时间的推移,稳定。在第一个12天,表面位移几乎没有变化。然而,围岩表面位移的突然增加了13天,稳定后再19天。作为测量分在5 #站位于残余煤柱下的围岩巷道严重变形,因为它是一个高应力集中。此外,5 # 1是最糟糕的部分车站的最大沉降和位移收敛达到328毫米和282毫米,分别。相比之下,屋顶的沉降收敛和收敛位移的双方在4 #站稍有变化,分别显示27毫米和91毫米。原因是测量两点的4 #站位于下固体煤,和巷道的围岩的应力小于煤柱的下部。

(2)深埋巷道围岩的位移。数据8(一个)和8 (b)显示深埋巷道围岩的位移曲线屋顶在4 #和5 #站。类似于巷道表面位移的变化规律,在深水中的每个测点的位移也稳定下来,然后最后增加和稳定。浅基础的point-depth测点测量2.5 m2分在4 #站位于下固体煤4 # 1点远的边界固体煤和深顶的位移是14毫米。4 # 2点安排在边界,这大大地受到集中应力比4 # 1。因此,深顶的位移是16毫米。,5 #站的两个计量点残余煤柱下排列,也极大地受到集中应力,位移是16毫米和13毫米,分别。的深基础point-depth测点4.5张深顶的位移测量两点的4 #站稍有变化,而点的位移5 # 1点18毫米5 # 2是12毫米。

3.3.3。巷道支架的应力变化

数据9(一个)- - - - - -9 (c)螺栓的应力变化(电缆)4 #和5 #站。随着时间的移动,螺栓的应力(电缆)在屋顶和两侧,然后稳定增加。5 # 1测点位于残余煤柱下完成的20314工作面开采;这里的压力是最集中,最支持身体变化的力量。具体来说,锚索的力量和螺栓屋顶增加40 kN, 32 kN,分别和螺栓的力量双方增加18 kN。相应地,4 # 1测点安排下固体煤,这是影响较小的停止P20314工作面。因此,在这一点上的力略有变化,部队的锚索和螺栓在屋顶增加20 kN和12 kN,分别。,在双方的前螺栓增加7 kN。

3.3.4。钻孔电视监测结果

一些水井的监测结果如图所示10。根据钻孔窥视的结果,它可以获得有许多断裂区域在0 - 2.72在4 # 2,在0 - 2.60在5 # 1。屋顶的岩石在这些地区严重变形和表达岩石裂缝。分离的现象很少出现。当偷窥钻井(深度)的范围在2.72∼2.72点(4 # 2和2.60∼2.60点5 # 1,屋顶摇滚有点变形,整体的完整性。此外,当偷窥钻井深度超过4.7米,围岩屋顶基本上是完整的,没有小层分离。应该注意,不能观察到的失败的步进式沉降引起的围岩顶板巷道的一边,因为电视钻井是安排在中间道路的早期阶段。

4所示。数值模拟在近距离煤层巷道的流变损伤

4.1。数值模型

为了进一步研究的变形和破坏规律的中间2 - 2煤层巷道工作面的停止,建立的数值计算模型是通过使用ANSYS软件。模型的大小是400×300×119,如图11。网格元素的20314辅助运输巷道(变形和破坏的主要领域,观察研究巷道的围岩)强化,最小长度是0.5米。此外,模型的水平和垂直方向是有限的。垂直载荷为5.0 MPa是应用在模型来代替上述未建模岩层。

应用FLAC3D软件用于数值模拟和莫尔-库仑本构模型是整体计算模型中使用。研究围岩的变形和破坏特征下煤层巷道的集中荷载,采用CVISC流变本构模型在巷道的围岩残余煤柱的有效面积40]。表2描述了参数的数值计算。在数值模拟过程中,巷道和工作面2 - 2煤层上部第一次挖掘形成采空区煤柱和部分,这些中间2 - 2煤层开挖顺序。

4.2。数值模拟结果

在2 - 2上煤层的开采,一个残余煤柱宽度的23之间形成20120采空区和20113年主要运输巷道。屋顶的垂直应力残留下的辅助运输巷道煤柱大于12 MPa(图12(一个))。自2 - 2煤层中间没有被发掘,岩石残余煤柱下的位移变化,从0.05米到0.15米(图(13日))。

如图12 (b)期间,中间2 - 2煤层的开采,残余煤柱的应力和较低的岩石释放由于20314辅助运输巷道的开挖。与此同时,岩石的垂直应力显著降低,是5∼9 MPa。没有明显的巷道的围岩位移变化。P20314工作面进一步发掘。这时,屋顶的应力集中的辅助运输巷道再次出现,和垂直应力达到9∼11 MPa(图12 (c))。此外,20314年的屋顶的位移(图辅助运输巷道会显著上升13 (c)),达到0.15∼0.3米。

此外,进化定律间隔下岩石的塑性区域残留煤柱在两个煤层的开采进行了分析。可以看出,开挖后的2 - 2煤层,塑性区出现在地板上(20314年的辅助运输巷道的屋顶)的残留煤柱(图(14日)),但变形范围很小。辅助运输巷道开挖后的20314年,塑料的岩石区域残余煤柱下稍微扩展(图14 (b))。当P20314工作面移动的活动区域残余煤柱,采空区的2 - 2煤层上部和中部煤层连接在一起。破碎的屋顶向上发展,增加屋顶荷载通过残余煤柱进一步向下转移。因此,围岩的塑性区中间的2 - 2煤层巷道的大幅增加,以及整个厚度的间隔在一些地区岩石变形(图14 (c))。

基于上述分析,屋顶的垂直应力变化的20314辅助运输巷道经历了三个阶段。在第一阶段,其余煤柱形成由于开挖2 - 2上的每个巷道和工作面煤层,导致巷道顶板的垂直应力急剧增加,应力集中和屋顶。在第二阶段,垂直应力释放的屋顶部分是由于20314年的辅助运输巷道的开挖,垂直应力降低,巷道顶板不变形。在第三阶段,煤柱的应力集中点向20113年集中主要运输巷道的停止是由于P20314工作面。这时,巷道顶板的垂直压力增加,造成的变形和破坏屋顶和巷道的围岩。

5。分析变形和破坏机理

通过数值模拟部分4压力的屋顶上,20314后挡板后不同的开挖工作面被重复,和压力的演化过程,过程中位移和塑性区。节5失效机理进行了分析,在此基础上进化过程。煤层上方的地层抵抗屈服和容易面临相当大的长度,因此屋顶悬臂梁可以形成采空区上方边缘工作面后撤退。根据现场数据和理论分析,可以看出,20314年的屋顶压力变化辅助运输巷道经历以下三个阶段。

5.1。阶段我

如图(15日),由于形成的悬臂梁的停止2012年上部煤层工作面,残余煤柱下的地层的压力迅速增加,煤柱边缘的显然是坏了。与此同时,“a”是弹性和塑性区交界处应力集中是显而易见的。结合数值模拟结果,可以得出结论,屋顶的塑性区辅助运输巷道的小完成后2 - 2煤层上部(图(16日)),“a”是压力峰面积。这时,顶板岩层的应力和位移变化的数据所示(17日)和(18日),分别。

5.2。第二阶段

矿业中间2 - 2煤层时,2 - 2上煤层的顶板运动稳定。自20314年辅助运输巷道的大小相对较小,扩大范围和位移收敛的塑性区巷道顶板也小当挖掘(数字16 (b)和18 (b)辅助运输巷道。压力的峰值位置基本不变。

5.3。第三阶段

如图15 (b)的一部分P20314工作面坐落在20120采空区。的停止P20314工作面,20120采空区向下继续下降与P20314工作面采空区的联系。停止前的20120工作面,屋顶的上覆地层的一部分已经坏了。然而,屋顶的地层完整,和承载力足够大的残余煤柱高应力。停止的工作面,上覆采空区破碎地层的进一步发展。悬臂,由于弯曲和断裂的地层沉降,增加残余煤柱上的压力。采空区破碎和塑料领域的增加,并且应力峰值点移动到煤柱的中心。“b”此外,位置略高于20314辅助运输巷道的屋顶,与数值模拟结果相一致的应力集中后的残余煤柱下巷道顶板的停止P20314工作面(数字16 (c)和17 (c))。

在这种结构下,运动的基本顶断块确定工作面巷道的稳定性。骨折跨度一个骨折块的长度有关b工作面和周期性断裂l:

工作面长度b1611.1米,是周期性断裂跨越吗l是12米,一个是14.1米。重量直接转移到煤层(41,42] 在哪里上覆地层的平均体积密度,25 kN / m³,h岩层的厚度,3 m,然后呢问528.8 kN / m。传送到煤层的力量 在哪里的最大应力= 15.9 MPa的断块煤层,=问罪棕褐色/H;H是断块的垂直距离煤层,43 m;是压力影响的角块,60°。用数据,我们可以发现= 18.4 MPa。60°角影响时,x小于H,所以可以看出= 15.9 MPa,大于岩石的长期强度层11.7 MPa,所以流变失败在这个条件下会发生。

基于现场监测,可以知道巷道的围岩变形与破坏的20314辅助运输巷道发生在很长一段时间。因此,结合理论分析和数值模拟结果,得到,顶板岩层的应力超过其承载力,然后导致蠕变,从而带来巷道的围岩变形与破坏的。此外,从图的数值模拟结果19可以看到,它的力量在顶部角落20314辅助运输巷道和采空区顶板明显比旁边的P20313工作面。这符合步进式沉降的一侧屋顶P20314工作面采空区。

结合实验室实验结果,集中应力的20314采空区上隅角巷道开挖后的屋顶P20314工作面11∼12 MPa。数据是类似于顶板岩层的长期强度,导致缓慢的失败20314辅助运输巷道残余煤柱下屋顶,导致变形和不稳定围岩的巷道。

6。讨论

屋顶螺栓20314辅助运输巷道的消除了巷道顶板的岩层之间的相互错位,使屋顶结合成一个更厚的屋顶结构的复合岩石梁。与此同时,由于顶板岩石梁厚度的增加,屋顶的整体稳定性提高,巷道压力降低。此外,屈服拱的跨度相应地减少。结合相关研究,复合悬浮岩石梁的力学模型建立,如图20.。在模型中,压力作用于屋顶的结构岩石梁被定义为问屋顶的支持压力锚索,螺栓,液压支柱被定义为问₁,被定义为两端的约束力量R_O和R_B。

根据力学平衡,可获得下列条件:

根据力平衡方程和力矩平衡可以获得如下:

屋顶上的任何部分弯矩岩石梁

屋顶的最大拉伸应力岩石梁

屋顶的条件没有失败岩石梁 在哪里l(米)代表的净宽度巷道,h₁(米)代表复合悬浮岩石梁的厚度,和螺栓顶部的有效长度,和问(Pa)代表的力量覆盖在岩石上时梁结构应用的压力。根据负载评估方法,它可以获得问=Kmγ在哪里米(米)是指的厚度和间隔岩石(包括上部煤层的厚度),直接的屋顶,屋顶和主要的上部煤层工作面,γ(N / m³)代表表土容重;K代表的动载荷系数,问₁是单体液压支柱的力,屋顶上的螺栓,锚索梁。此外,问₁简化为一个向上立统一的负载(Pa),可以计算吗 在这象征着螺栓的锚固力,P_年代电缆的锚固力,P_z是指额定工作阻力单一液压道具,年代_b代表锚之间的间距,米每一行是锚的数量。锚的间隔一行nS_b是不可或缺的多个行距离的螺栓n将一个正整数1、2、3…,σ_马克斯(Pa)象征着屋顶的最大拉伸应力岩石梁、和σ(Pa)的抗拉强度是屋顶锚固后地层。

辅助运输巷道的P20314工作面,净宽度l5.2米,有效长度吗h₁螺栓是2.05米,顶部和间距年代_b主持人之间的0.8米。图3被选中的米和n。的最大厚度总和(米间隔的岩石(包括上部煤层),直接的屋顶,屋顶和主要的工作面上部煤层是30 m。统一的负载问₁的锚索、螺栓和单体液压支柱在屋顶上岩石梁是356.73 kN / m²(和P_年代锚的锚固力80 kN, 164 kN,分别单液压支柱的额定工作阻力250 kN)。

6.1。上覆地层的压力下重量

用上面的参数方程(9),它可以获得最大拉应力在屋顶上岩石梁为1.9 MPa。根据相关研究,岩石的抗拉强度通常在锚定可以达到5 MPa 8 MPa。它可以知道屋顶岩石梁上的最大拉伸应力小于屋顶锚固后的抗拉强度。因此,支持系统可以控制围岩的作用只在考虑上覆岩层重量。

6.2。在压力下的应力集中

的停止P20314工作面,20314年辅助运输巷道坐落在残余煤柱边缘,导致两个工作面临的附加轴承压力(工作面临20113年和20120年的2 - 2煤层上部)。因此,由于应力集中系数3。用相关参数方程(9),它可以知道,屋顶上的最大拉应力岩石梁为9.14 MPa。根据相关研究,岩石的抗拉强度在锚定可以一般从5 MPa 8 MPa。可以看出,最大拉应力在屋顶上岩石梁大于锚定后顶板岩石的抗拉强度。结果,在考虑应力集中时,支持系统不能有效地控制围岩和辅助支持等措施支持必须采用。

因此,降低煤层的巷道应合理的安排,以避免其围岩的变形和破坏过程中多个和近距离煤层的开采43,44]。原因在于,降低煤层的巷道不同固体上部煤层的煤,残余煤柱和采空区。通过一个完美的安排,它可以避免的问题的残余煤柱巷道穿透煤层上部和造成的变形和破坏现象过度巷道的围岩压力。如果道路安排必须穿透残余煤柱上煤层由于一些原因如煤层发生故障,可以采取一些措施根据特定站点的生产地质条件。在这种情况下,支持高强度和high-pretension螺栓,积极支持锚索加固,被动的支持组合钢梁和单一液压支柱,与围岩注浆加固可以选择(41,45- - - - - -47]。

7所示。结论

摘要弱胶结软岩的流变破坏机制下巷道残留煤柱在浅埋近距离煤层的开采Gaojialiang煤矿通过现场监测进行了分析和研究,数值模拟和实验室实验。现场监测结果表明,每个监控装置的数据辅助运输巷道的同时开始改变,他们是稳定后再改变同时以较慢的速度。此外,受损区域的监测结果观察到井下钻孔电视是一致的与位移的屋顶,这证明监测数据的可靠性。

数值模拟的结果表明,巷道顶板的应力增加,然后最后增加和减少。同时,挖掘后P20314工作面,屋顶的应力和位移的残留下的辅助运输巷道煤柱集中和增加,分别。结果符合应力和位移的变化数据在5 # 1车站现场监测。此外,顶部的压力P20314工作面采空区更集中,这也符合步进式沉降的屋顶。

实验室试验表明,屋顶层属于弱胶结地层,和试样蠕变压力小。在此基础上,20314年辅助运输巷道的破坏机理进行了分析。(1)煤柱离开是因为上部煤层的开采。与上部煤层开挖应力的屋顶上20314年的辅助运输巷道煤柱下逐渐增加和集中。然而,20314年的挖掘带来的压力在屋顶部分释放,进一步降低。因此,巷道围岩的完整性。(2)的停止P20314低煤层工作面,重叠部分的采空区上部和下部煤层被连接起来,和上覆岩石破碎的转动使残余煤柱的应力增加。集中应力转移到辅助运输巷道顶部的屋顶通过残余煤柱和应力集中区域位于正上方的屋顶。但是,由于浅埋煤层,集中应力不够高直接导致巷道顶板的变形。相反,屋顶爬随着时间的推移,在这样的条件下,巷道的围岩变形。(3)通过分析,建议在工作系统的设计,gateroad在一个孤立的残余煤柱上煤层开采时应避免近距离煤层。如果无法避免,被动的支持加上积极措施应该用于支持道路,确保煤矿安全生产。

数据可用性

使用的数据来支持本研究的结果包括在本文中。

的利益冲突

作者宣称他们没有利益冲突有关的出版。

确认

作者欣然承认,这项工作得到了中国国家重点研发项目(2018 yfc0604703),山东省自然科学基金(ZR2017BEE013),国家重点实验室的研究基金煤炭资源安全开采,CUMT (SKLCRSM19KF016),中国国家自然科学基金(51574154和51574154),山东省自然科学基金的主要程序(ZR2018ZA0603),山东省重点研发项目(2018 gsf116003),和科学研究基金会的山东科技大学招聘人才(2017 rcjj009)。