3DEC煤层巷道开挖钻孔泄压， 地层分布模型， 地应力平衡， 巷道开挖， 分层钻孔泄压， 显示位移，裂隙形成过程，
巷道开挖遇到高应力煤层是真头疼，钻孔泄压这招用好了能救命。今天咱们用3DEC撸个完整流程，从建模到裂隙发育全过程走一遍。别急着跑模型，先把地层结构理清楚——上方15米砂岩，中间2米煤层，底下8米泥岩，这个组合够典型了吧。
建模这块直接上代码最实在：
block create brick -50 50 -30 30 0 25 block cut joint-set dip 90 dip-direction 0 origin 0,0,15 spacing 5 block cut joint-set dip 90 dip-direction 90 origin 0,0,15 spacing 3砂岩层节理间距5米，到了煤层直接加密到3米，这才是符合实际的各向异性。注意z=15的位置是煤层顶板，切割的时候别手滑。
地应力平衡容易翻车，关键在边界条件设置：
block face apply stress -12e6 0 0 range pos-x -50,50 pos-y -30,30 pos-z 24.9,25.1 block face apply stress 0 -8e6 0 range pos-y -30.1,-29.9 block property mat=1 density 2650 young 18e9 poisson 0.25 block property mat=2 density 1450 young 3.5e9 poisson 0.31水平应力给12MPa，垂直8MPa，这个比值刚好能触发煤岩差异变形。材料参数里煤层杨氏模量只有砂岩的1/5，软硬夹层效应就出来了。
开挖模拟要分步走，别一次性挖穿：
block tunnel face 0,0,10 direction 1,0,0 height 3.5 width 4.5 block relax range group 'tunnel' block solve每掘进5米暂停计算，这时候位移场还没稳定呢。见过新手直接挖20米，结果模型直接发散，煤壁位移飙到半米——现实里早塌方了。
钻孔泄压的骚操作来了：
block drill hole origin -5,0,13 direction 1,0,0 length 10 diameter 0.1 block drill hole origin -5,0,14 direction 1,0,0 length 10 diameter 0.1 block property mat=3 range group 'drill_holes'上下两排孔径10cm，间距1米。重点是把钻孔材料设为低刚度介质，用property mat=3 young 1e8模拟卸压效果。注意钻孔要在开挖面前方5米提前施工，这叫应力超前释放。
运行到3000时步的时候，调用这个看裂隙：
block plot contour displacement magnitude block plot crack这时候会看到煤帮表面先出现龟裂，随后裂隙沿着节理面向深部扩展。有个细节——卸压孔周围的裂纹呈现放射状，而两孔之间容易形成贯通裂纹，这个和现场监测结果对得上。
最后来个数据对比：没卸压时顶板下沉42cm，卸压后降到17cm。但别光看位移量，注意看应力重分布——原始集中应力系数2.1，卸压后降到1.3，这才是保命的关键。新手常犯的错是只盯着位移看，忽略了应力场演化。
模拟到5000时步记得保存状态，因为这时候可能突然出现煤炮现象——应力曲线会有一个陡降，对应现场的能量释放。模型里虽然没法完全复现动力效应，但通过声发射计数模块还是能捕捉到前兆信号的。
