矿山隧道掘进围岩稳定性动态监测研究.pdf

矿山隧道掘进围岩稳定性动态监测研究 ① 关祥宏 （中铁第四勘察设计院集团有限公司，湖北 武汉 ４３００６３） 摘 要 在考虑矿区隧道围岩级别的同时，以矿区地质力学特征为基础，采用对数螺线法对矿山隧道掘进围岩稳定性进行分析，以 隧道围岩沉降、位移与埋深的关系为研究对象，进行了矿山隧道掘进围岩稳定性动态监测验证。 结果表明侧壁、底部沉降与埋深 呈线性关系；侧壁、底部位移与埋深呈线性关系；矿区隧道掘进过程中可通过降低沉降或减少位移来增加围岩稳定性。 关键词 隧道； 掘进； 围岩； 稳定性； 监测； 对数螺线 中图分类号 Ｕ４５１．２；Ｕ４５５．６文献标识码 Ａｄｏｉ１０．３９６９／ ｊ．ｉｓｓｎ．０２５３－６０９９．２０２０．０３．００９ 文章编号 ０２５３－６０９９（２０２０）０３－００３４－０５ Ｄｙｎａmｉｃ Ｍｏｎｉｔｏｒｉｎｇ ｏｆ ｔｈｅ Ｓｕｒｒｏｕｎｄｉｎｇ Ｒｏｃｋ Ｓｔａｂｉｌｉｔｙ Ｄｕｒｉｎｇ ｔｈｅ Ａｄｖａｎｃｅmｅｎｔ ｏｆ Ｍｉｎｅ Ｔｕｎｎｅｌｓ ＧＵＡＮ Ｘｉａｎｇ-ｈｏｎｇ （Ｃｈｉｎａ Ｒａｉｌｗａｙ Ｓｉｙｕａｎ Ｓｕｒｖｅｙ ａｎｄ Ｄｅｓｉｇｎ Ｇｒｏｕｐ Ｃｏ Ｌｔｄ， Ｗｕｈａｎ ４３００６３， Ｈｕｂｅｉ， Ｃｈｉｎａ） Ａｂｓｔｒａｃｔ Ｃｏｎｓｉｄｅｒｉｎｇ ｔｈｅ ｇｒａｄｅ ｏｆ ｓｕｒｒｏｕｎｄｉｎｇ ｒｏｃｋ ｉｎ ｔｈｅ ｍｉｎｅ ｔｕｎｎｅｌｓ， ａ ｌｏｇａｒｉｔｈｍ-ｓｐｉｒａｌ ｍｅｔｈｏｄ ｗａｓ ｕｓｅｄ ｔｏ ａｎａｌｙｚｅ ｔｈｅ ｓｔａｂｉｌｉｔｙ ｏｆ ｔｈｅ ｓｕｒｒｏｕｎｄｉｎｇ ｒｏｃｋ ｉｎ ｔｈｅ ｅｘｃａｖａｔｉｏｎ ｏｆ ｍｉｎｉｎｇ ｔｕｎｎｅｌｓ ｂａｓｅｄ ｏｎ ｔｈｅ ｇｅｏｌｏｇｉｃａｌ ａｎｄ ｍｅｃｈａｎｉｃａｌ ｃｈａｒａｃｔｅｒｉｓｔｉｃｓ ｏｆ ｔｈｅ ｍｉｎｉｎｇ ａｒｅａ． Ｂａｓｅｄ ｏｎ ｔｈｅ ｓｔｕｄｉｅｓ ｏｆ ｔｈｅ ｒｅｌａｔｉｏｎｓｈｉｐ ａｍｏｎｇ ｔｈｅ ｓｅｔｔｌｅｍｅｎｔ， ｄｉｓｐｌａｃｅｍｅｎｔ ａｎｄ ｂｕｒｉａｌ ｄｅｐｔｈ ｏｆ ｔｈｅ ｓｕｒｒｏｕｎｄｉｎｇ ｒｏｃｋ ｄｕｒｉｎｇ ｔｈｅ ｔｕｎｎｅｌ ｅｘｃａｖａｔｉｏｎ， ｔｈｅ ｄｙｎａｍｉｃ ｍｏｎｉｔｏｒｉｎｇ ｏｆ ｔｈｅ ｓｕｒｒｏｕｎｄｉｎｇ ｒｏｃｋ ｓｔａｂｉｌｉｔｙ ｗａｓ ｖｅｒｉｆｉｅｄ． Ｒｅｓｕｌｔｓ ｓｈｏｗｅｄ ｔｈａｔ ｔｈｅ ｓｕｂｓｉｄｅｎｃｅ ｏｆ ｓｉｄｅ ｗａｌｌ ａｎｄ ｆｏｏｔ ｗａｌｌ ｅｘｈｉｂｉｔｓ ａ ｌｉｎｅａｒ ｒｅｌａｔｉｏｎｓｈｉｐ ｗｉｔｈ ｔｈｅ ｂｕｒｉａｌ ｄｅｐｔｈ， ａｎｄ ｔｈｅ ｄｉｓｐｌａｃｅｍｅｎｔ ｏｆ ｓｉｄｅ ｗａｌｌ ａｎｄ ｆｏｏｔｗａｌｌ ａｌｓｏ ｅｘｈｉｂｉｔｓ ａ ｌｉｎｅａｒ ｒｅｌａｔｉｏｎｓｈｉｐ ｗｉｔｈ ｔｈｅ ｂｕｒｉａｌ ｄｅｐｔｈ． Ｉｔ ｉｓ ｃｏｎｃｌｕｄｅｄ ｔｈａｔ ｔｈｅ ｓｕｒｒｏｕｎｄｉｎｇ ｒｏｃｋ ｓｔａｂｉｌｉｔｙ ｃａｎ ｂｅ ｉｍｐｒｏｖｅｄ ｂｙ ｒｅｄｕｃｉｎｇ ｔｈｅ ｓｕｂｓｉｄｅｎｃｅ ｏｒ ｄｉｓｐｌａｃｅｍｅｎｔ ｄｕｒｉｎｇ ｔｈｅ ｔｕｎｎｅｌ ａｄｖａｎｃｅｍｅｎｔ． Ｋｅｙ ｗｏｒｄｓ ｔｕｎｎｅｌ； ａｄｖａｎｃｅｍｅｎｔ； ｓｕｒｒｏｕｎｄｉｎｇ ｒｏｃｋ； ｓｔａｂｉｌｉｔｙ； ｍｏｎｉｔｏｒｉｎｇ； ｌｏｇａｒｉｔｈｍｉｃ-ｓｐｉｒａｌ 矿山隧道掘进包括隧道内的施工测量防线、洞身 开挖、弃渣运输等工序［１］。 其中洞身开挖使用的机械 设备会对矿山隧道围岩稳定性造成影响，若事先未进 行岩石验证，会对后期矿山隧道带来一定的安全隐 患［２］。 因此，对矿山隧道掘进围岩稳定性进行动态监 测具有重要意义。 对矿山隧道掘进围岩稳定性进行动态监测，既可 以确定矿山开采方式，又可以为矿山安全开采提供可 靠参考。 近年来，我国学者针对矿山隧道掘进技术进 行了优化研究［３－７］。 本文在考虑矿区隧道围岩级别的 同时，以矿区地质力学特征为基础，采用对数螺线法对 矿山隧道掘进围岩稳定性进行分析，以隧道围岩沉降、 位移与埋深的关系为研究对象，进行了矿山隧道掘进 围岩稳定性动态监测验证。 １ 工况概述 此次研究对象为湖北省地下深部煤矿矿山隧道，该 煤矿井田中含有一条地下 ８００ ｍ 轨道开拓巷道，巷道的 主要作用是行人通行、通风以及资源运输。 巷道所在围 岩为粉细砂岩，在隧道不断掘进过程中，由于巷道较深， 受到明显的构造应力作用，导致隧道围岩严重变形形成 破碎区，因此对该矿山隧道掘进围岩稳定性进行分析。 ２ 矿山隧道掘进围岩稳定性分析 ２．１ 隧道围岩关键块力学模型 隧道围岩在掘进过程中，掘进机产生的振动作用 在周围围岩上，导致周围围岩形成裂痕严重时会出现 破裂。 而相比近距离的掘进，较远处的围岩却处于弹 ①收稿日期 ２０２０－０１－０５ 作者简介 关祥宏（１９６４－），男，河南辉县人，高级工程师，主要从事工程测量与工程监测工作。 第 ４０ 卷第 ３ 期 ２０２０ 年 ０６ 月 矿矿 冶冶 工工 程程 ＭＩＮＩＮＧ ＡＮＤ ＭＥＴＡＬＬＵＲＧＩＣＡＬ ＥＮＧＩＮＥＥＲＩＮＧ Ｖｏｌ．４０ №３ Ｊｕｎｅ ２０２０ 性介质范围内［８］。 因此，一般矿区开发的隧道围岩应 力场能够使用经典的弹性力学模型，即假设围岩岩体 处于介质连续均匀状态，符合弹性范围，则矿区隧道围 岩关键块受力情况如图 １ 所示。 图中 ｒ 为隧道中心 Ｏ 到围岩裂痕的距离，ｍ；ａ 为圆形隧道半径，ｍ。  θ O a r σxσx σZ σZ 图 １ 圆形隧道围岩应力状态示意 ２．１．１ 关键块侧面平均正应力 隧道围岩静态应力场中不同角度力的计算表达式 分别为 σｒ＝ σｚ ＋ σ ｘ ２ １ － ａ２ ｒ２ － σｚ ＋ σ ｘ ２ １ ＋ ａ４ ｒ４ － ４ａ２ ｒ２ ｃｏｓ２θ （１） σθ＝ σｚ ＋ σ ｘ ２ １ ＋ ａ２ ｒ２ ＋ σｚ － σ ｘ ２ １ ＋ ３ａ４ ｒ４ ｃｏｓ２θ （２） τｒθ＝ σｚ － σ ｘ ２ １ ＋ ２ａ２ ｒ２ － ３ａ４ ｒ４ ｓｉｎ２θ（３） 式中 σｘ ＝ｋ ０σｚ，σｚ为抗压强度，ＭＰａ。 在隧道顶部围岩处，θ＝９０，当 ｒ＝ａ 时，σθ＝－σｚ ＋３σ ｘ； 当 ｒ＝ ２ａ 时，σθ＝ σｚ＋３９σｘ ／３２。 当围岩的有效水平 应力和垂直应力比值 ｋ０不同时，隧道掘进过程中围岩 顶部的环向正应力 σθ可能为拉应力。 依据有关调查 研究分析可知，在埋深 ４００ ｍ 左右处进行掘进时，不同 种类的围岩具备不同的侧压系数，其中岩浆岩层侧压 系数最大，变质岩层的侧压系数最小，沉积岩则处于此 二者之间［９］。 一般情况下，设置隧道底部的围岩测压次数在 ０．８～６ 之间，变质岩在 ０．７～４ 之间，沉积岩在 ０．５～４ 之 间。 可见，ｋ０＞１／３ 时，则－σｚ＋３σｘ＞０，当隧道顶部围岩 处于压应力状态，通过积分运算，得到关键围岩块体的 侧面平均正应力 σ 计算式为 σ ＝ １ ｈ∫ ａ＋ｈ ａ σθｄｒ ＝ σ ｘ ＋ ａ（σｚ ＋ σ ｘ） ａ（ａ ＋ ｈ） ＋ （σｚ ＋ σ ｘ） ａ ４ － ａ（ａ ＋ ｈ）３[] ２ｈ（ａ ＋ ｈ）３ （４） 式中 ｈ 为关键块的高度，ｍ。 由此可以看出，围岩关键 块体侧向平均正应力的值与埋深、洞径及关键块体高 度相关。 ２．１．２ 关键块体临界荷载 当围岩关键块是长方体时，设置其高度为 ｈ，底面 长为 ｂ、宽为 ｄ，则在掘进过程中，当出现扰动荷载时， 关键块受力如图 ２ 所示。 图中 Ｇ 为块体自身重力；ｆ１、 ｆ２为从不同面受到的冲击力。 Fd Fr h b d G f1 f2 图 ２ 围岩顶部关键块体受力图 关键块体的自重 Ｇ ＝ γｂｄｈ（５） 式中 γ 为块体数量。 在弹性力学中，接触面正应力相等的条件下，能够 获取关键块体上表面静态压力为 Ｆｒ＝ ηｂｄＧ（６） 式中 η 为接触面正应力，ＭＰａ。 此时忽略掘进机自身对块体侧向产生的振动作 用，依据掘进速度导致的振动作用，得到关键块体上的 荷载与竖向速度关系式为 Ｆｄ＝ ρＣｐＶＦｒ（７） 式中 ρ 为围岩密度，ｋＮ／ ｍ３；Ｃｐ为块体结构面的粘聚 力，ＭＰａ；Ｖ 为竖向振动速度，ｍ／ ｓ。 关键体侧面阻滑力分为两部分，以 ｘ 轴法向侧面 阻滑力表达式为 ｆ１＝ ２βｈＦｄ（σｔａｎφ ＋ Ｃ）（８） 式中 φ 为结构面壁岩的内摩擦角，（）；Ｃ 为结构面壁 岩的粘聚力，ＭＰａ；β 为结构面抗剪力系数，与结构面 中凸起体大小和数量有关。 假定矿区隧道掘进过程中两个水平方向的压力相 等σｙ ＝σ ｘ ［１０］，以 ｙ 轴法向侧面阻滑力的表达式为 ｆ２＝ ２βｈＦｒ（σｘｔａｎφ ＋ Ｃ）（９） ５３第 ３ 期关祥宏 矿山隧道掘进围岩稳定性动态监测研究 那么联立两个方向的侧向阻力为 ｆ ＝ ２βｈσｄ ＋ σｘｂ ｔａｎφ ＋ｂ ＋ ｄ Ｃ []（１０） 依据刚体平衡原理，掘进过程中关键块遇到扰动 荷载时，达到临界平衡条件为 ｆ临＝ ｆ１ ＋ ｆ ２ ｆ （１１） ２．２ 围岩剪切支护阻力矩分析 由弹性力学理论可以得到圆形隧道围岩的剪切应 力解，取最大切应力区位于 ４５处的块体进行受力分 析，如图 ３ 所示。 Fd Fr h b d G f f 图 ３ 围岩最大切应力区关键块体受力图 地应力场产生的剪应力在块体两侧呈径向平衡， 导致块体发生剪切旋转［１１］。 因此，依据推导结果和隧 道的径向刚体平衡原理，能够得到剪应力区块体振动 速度的计算公式为 Ｖｍａｘ＝ ２βｈσｄ ＋ σｘｂ ｔａｎφ ＋ｂ ＋ ｄ Ｃ [] ρＣｐｂｄ － η ρＣｍ （１２） 其中 σ ＝ １ ｈ∫ ａ＋ｈ ａ σθｄｒ ＝ σｘ＋ ａ σｚ ＋ σ ｘ ２ ａ ＋ ｈ 在隧道掘进过程中围岩的掌子面存在不稳定因 素，主要包括前方产生的滑动土块的自重 Ｗｆ及滑动土 块上部松动区域的松动土压 Ｑ，将以上力对 Ｏ 点取力 矩平衡，就可以得到开挖面的稳定系数 ｆ′ ＝ ＭＺ ＋ Ｍ ｔ Ｍｆ ＋ Ｍ Ｑ ＝ ＭＺ ＋ Ｍ ｔ Ｗｆｌｆ＋ Ｑｌａ （１３） 式中 Ｍｚ为约束措施产生的支护阻力矩，Ｎｍ；Ｍｔ为土 体抗力矩，Ｎｍ；Ｍｆ为滑动土体自重力矩，Ｎｍ；ＭＱ为 松动土力矩，Ｎｍ。 在不加约束时，若要预测掌子面的稳 定性，则掘进面的稳定系数为１，其力矩应满足平衡条件 Ｍｔ≥ Ｍｆ ＋ Ｍ Ｑ （１４） 对数螺线的方程为 Ｒ ＝ Ｂ ２ ２ ＋ Ｆ ２ Ｆ ｃｏｓ π ４ － φ ２ ｅｘｐ π ４ － φ ２ ｔａｎφ ｆ′ ｅθｔａｎφ （１５） Ｍｆ计算公式为 Ｍｆ ＝ Ｍ ＯＡＣ － Ｍ ＯＡＢ － Ｍ ＯＢＣ （１６） 计算松动土压力矩 ＭＱ ＭＱ ＝ ｌ ｂ{ γＢ′ － ｃ Ｋ０ｔａｎφ １ － ｅｘｐ － ＨＫ０ｔａｎφ Ｂ′ ＋ ｑｅｘｐ － ＨＫ０ｔａｎφ Ｂ′ } ｌ ａ ＋ ｌｂ ２ （１７） 式中 ｑ 为荷载，ｋＮ／ ｍ２；ｑ０为地面荷载，ｋＮ／ ｍ２；Ｋ０为有 效水平应力和垂直应力比值；Ｈ 为隧道掘进的垂直影 响范围。 围岩体抗剪强度为 τｔ＝ ｃ ＋ ｔａｎφ（１８） 计算围岩体抗力矩 Ｍｔ Ｍｔ＝ Ｒ０２ｃ ｅｘｐ ２ θ１ － θ ２ ｔａｎφ[] － １{} ｃｏｓφ ２ｔａｎφ （１９） 当 Ｍｔ＜Ｍｆ ＋Ｍ Ｑ时，显然此时隧道掌子面处于危险 状态，一般需要给予加固约束措施才能保证掘进过程 中掌子面的稳定性。 支护阻力矩的计算公式为 Ｍｚ ＝ ｌ ｂσｔ Ｒ０ｓｉｎφ ＋ Ｍｔ ２ （２０） 以上判定条件皆是基于掌子面无任何人工支护或 支护阻力的情况下围岩的稳定性分析，而实际上在矿 区隧道掘进过程中，一般要做到超前支护措施，而 ＡＤＥＣＯ-ＲＳ 施工时更有超前预约束措施，其范围达到 隧道开挖宽度的 ２～４ 倍。 ３ 稳定性监测验证 矿区隧道掘进围岩稳定性的动态监测主要通过监 测掘进过程中隧道围岩产生的振动位移来研究稳定 性。 对此，在掘进之前埋设监测点，在掘进的同时采集 位移值，其监测点布置如图 ４ 所示。 每间隔 ５ ｍｉｎ 监 测 １ 次，持续监测 ３ ｄ。 1 2 3 4 5 6 7 8 9 图 ４ 测点分布图 实验建模过程中，需要考虑掘进对隧道围岩产生 的边界效应，确定隧道模型的横断面两侧及下部分为 ６３矿 冶 工 程第 ４０ 卷 ３ 倍洞径，上部为实际地表。 模型边界四周约束水平 方向，底部约束竖直方向，上部为自由面，围岩模型如 图 ５ 所示。 图 ５ 数值模拟模型 通过数值模拟得到矿区隧道掘进围岩位移的理想 模型和实际得到位移模型如图 ６ 所示。 图 ６ 数值模拟模型 （ａ） 理想模型； （ｂ） 实际位移模型 由图 ６ 可以看出，在掘进过程中矿区隧道围岩发 生破坏主要集中在侧面及底板处。 此外，隧道的拱腰 处和拱肩处均产生了破损情况。 由以上分析可以看 出，在掘进过程中应在侧壁及底板处加强支护，提高锚 喷强度以提高围岩稳定性。 经 ３ ｄ 持续监测的数据可直接将围岩级别、埋深 位移值进行提取，绘制出矿区隧道围岩相对变形在不 同级别及埋深下的关系图如图 ７ 所示。 从图 ７ 可以看 出，共同规律是相对位移都随埋深增加而增加，随着围 岩级别降低而减小。 且各围岩界别之间并没有明显的 界限，也印证了各围岩级别的位移是在一定范围内。 由此可知，隧道埋深在 １５０ ｍ下，位移取 １００ ～ ２００ ｍ 埋深段的中间值，则拱脚数量位移Ⅲ级为 ０．０９７ ６％， Ⅳ级为 ０．１６９ ６％，Ⅴ级为 ０．２３５ ２％，此数据在此只作 为定性论证。 5m 1.8 1.5 1.2 0.9 0.6 0.3 0.0 20004006008001000 1 5m  1.5 1.0 0.5 0.0 20004006008001000 1 0 0 0   0 0 0   图 ７ 沉降与埋深关系散点图 不同围岩级别下的周边相对位移与隧道埋深的关 系图如图 ８ 所示。 5m  1.6 1.2 0.8 0.4 0.0 20004006008001000 ,AD 0 0 0 0     5m  1.6 1.2 0.8 0.4 0.0 20004006008001000 ,AD 0 0 0 0    图 ８ 位移与埋深关系图 从图 ８ 可以看出，位移的变化规律符合不同围岩 ７３第 ３ 期关祥宏 矿山隧道掘进围岩稳定性动态监测研究 级别与埋深的基本规律，侧壁位移及底部位移均随着 埋深增加而增加，其中Ⅲ级围岩位移变化较为显著，且 呈现规律性变化。 显然，各级围岩的周边位移量与埋 深的相关性很强，可以根据回归的方式寻求这其中的 关系。 ４ 结 论 利用对数螺线法可以有效地确定围岩的稳定性， 同时可知侧壁、底部沉降与埋深呈线性关系；侧壁、底 部位移与埋深呈线性关系；矿区隧道掘进过程可通过 降低沉降或减少位移来增加围岩稳定性。 参考文献 ［１］ 张 洋，刘陕南，吴 俊，等． 盾构隧道掘进时地层参数变化对地 表沉降的敏感性研究［Ｊ］． 现代隧道技术， ２０１９，５６（４）１２７－１３４． ［２］ 刘小鸣，陈士海． 隧道掘进中掏槽孔爆破引起的地表振动波形预 测［Ｊ］． 岩土工程学报， ２０１９，４１（９）１７３１－１７３７． ［３］ 周 磊，朱哲明，刘 邦． 裂纹对直墙拱形隧道围岩损伤破坏模式 的影响规律研究［Ｊ］． 岩土力学， ２０１７（１２）３６８８－３６９７． ［４］ 汪学清，石瑞明，曹新玲，等． 月亮山隧道围岩变形规律及稳定性［Ｊ］． 科学技术与工程， ２０１７（３１）３３９－３４５． ［５］ 刘建刚，郁钟铭． 近距离薄煤层联采围岩稳定性及合理错距研究［Ｊ］． 煤炭技术， ２０１７，３６（３）４９－５１． ［６］ 关 玲． 不同的开挖方式对公路隧道围岩稳定性影响分析 以 重庆公路隧道为例［Ｊ］． 公路工程， ２０１９（１）１５６－１６０． ［７］ 魏秀泉． 采空区群对１＃竖井稳定性影响数值分析及实时监测［Ｊ］． 矿冶工程， ２０１７，３７（６）２６－２９． ［８］ 李梓源，王海亮． 基于钻爆法下的爆破振动对围岩稳定性的数值 模拟研究［Ｊ］． 公路， ２０１７（１２）２８８－２９１． ［９］ 惠 鑫，马凤山，郭 捷，等． 金属矿山深部巷道围岩分区破裂现 场监测研究［Ｊ］． 工程地质学报， ２０１８，２６（Ｓ１）７１０－７１５． ［１０］ 张 杰，路增祥，杨宇江． 地下选厂硐室开挖过程围岩应力动态 分布特征分析［Ｊ］． 矿冶工程， ２０１８，３８（４）２７－３１． ［１１］ 崔建斌，姬安召，熊贵明． 基于 Ｓｃｈｗａｒｚ Ｃｈｒｉｓｔｏｆｆｅｌ 变换的圆形隧 道围岩应力分布特征研究［Ｊ］． 应用数学和力学， ２０１９（１０） １０８９－１０９８． 引用本文 关祥宏． 矿山隧道掘进围岩稳定性动态监测研究［Ｊ］． 矿冶 工程， ２０２０，４０（３）３４－３８． �������������������������������������������������������������������������������������������������� 矿冶工程杂志 ２０２１ 年征订启事 矿冶工程（双月刊）由中国金属学会、长沙矿冶研究院有限责任公司主办，面向国内外公开发行。 本刊是 中国期刊方阵“双效期刊”、全国中文核心期刊、中国科学引文数据库（ＣＳＣＤ）及中国学术期刊综合评价数据 库来源期刊、中国核心学术期刊（ＲＣＣＳＥ），是集学术性和技术性于一体的综合性刊物，已被中国知网（ＣＮＫＩ）、 万方数据库、重庆维普资讯、台湾华艺数据库等全文收录，是国外多家知名检索机构的检索对象。 矿冶工程读者对象是采矿、选矿、冶金、材料、地质、化工等系统的有关生产人员、院校师生和管理人员。 主要栏目为采矿、选矿、冶金、材料、矿冶行业企业管理等，内容新颖，是开拓、激发创造力的良师益友。 矿冶工程编辑部承接彩色、黑白及文字广告业务，欢迎各企事业单位来电来函联络。 矿冶工程真诚欢迎新、老订户向全国各地邮局订阅本刊，也可直接向编辑部订阅。 邮发代号４２－５８，大 １６ 开，定价 ２０ 元，全年 １２０ 元。 地 址 湖南省长沙市麓山南路 ９６６ 号联 系 人 黄小芳 邮 编 ４１００１２开户名称 矿冶工程杂志（长沙）有限公司 电 话 （０７３１）８８６５７０７０／８８６５７１７６／８８６５７１７３开户银行 工商银行长沙左家垅支行 传 真 （０７３１）８８６５７１８６帐号 １９０１０１３００９２０１０９５５０２ Ｅ－ｍａｉＩ ｋｕａｎｇｙｅｇｏｎｇｃｈｅｎｇｚｚ＠ １６３．ｃｏｍ网址 ｗｗｗ．ｋｙｇｃ．ｎｅｔ ８３矿 冶 工 程第 ４０ 卷