一种预测立筒井筒冻结情况的方法
【专利摘要】一种预测立筒井筒冻结情况的方法，包括以下步骤:(1)测量或查看冻结管的各参数；(2)确定每天送冷量和推出单根冻结管第n天冻土发展半径公式；(3)推出计算冻结交圈时间公式；(4)推出冻土发展到井帮所需要的冻结天数N计算公式；(5)推出受单圈孔影响的井帮温度计算公式；(6)推出受两个或以上冻结孔影响的井帮温度的计算公式；(7)确定距离影响系数ε的数值；(8)将测量或查看冻结管的各参数和确定距离影响系数ε的数值代入第(5)步中的计算公式中；其中步骤(1)步骤(3)可以放在步骤(8)后。该步骤方法避免了经验法和工程类比法的苛刻条件，以科学的方式提前获得数据。
【专利说明】
一种预测立筒井筒冻结情况的方法
技术领域
[0001] 本发明涉及立井井筒领域，尤其涉及一种预测立筒井筒冻结情况的方法。
【背景技术】
[0002] 冻结技术己广泛应用到立井和隧道等地下工程中，冻结技术就是对未知和复杂的 表土层进仔预先冻结处理。立井工程中表土段采用冻结技术施工通过。近20年来，采用冻结 法施工的立井达到几百个。比如，当冻结深度超过400米，各国冻结孔的设计难度较大，现有 技术是全冻结段上下同时送冷或停冻，这样容易造成井帮温度上下不均匀的情况，从而给 井筒掘砌施工带来巨大的困难和安全隐患，造成成本浪费和工期延误。
[0003] 冻结法施工立井井筒，冻结壁温度场发展的预测是非常重要的，预测的准确与否 直接影响到井筒掘砌的速度和下一阶段的施工安排，随着冻结深度的不断增大，冻结施工 中很难做到精细控制和准备把握，在通常情况下是超冻或者过冻，浪费了大量的电力和材 料消耗，同时也给井筒安全掘砌带来很多不确定的安全隐患，造成冻结立井施工时常发生 断管、透水等严重工程事故。在总结血的教训的同时，广大技术科技人员都在努力研究更加 准确可靠的立井冻结温度场计算与预测理论，期望能够更加准确的掌握冻结壁发展状况。
[0004] 如果在工程设计方案之前就做好预测工作，可以更好的指导冻结方案设计工作， 从技术源头抓好掘砌和冻结的协调关系。保证立井掘砌安全最根本、也是最有效的方法就 是制定一个好的冻结方案，也是避免安全问题的最好的"预防针"。因此，在冻结方案设计完 成后应根据该方案首先进行预测，再根据预测的结果进一步优化方案，才能最大程度地保 证掘砌的安全、连续进行。
[0005] 虽然目前施工单位和科研单位都在积极研究冻结计算和模拟技术，其中就包括 了冻结预测技术，但大多都处在完善和实验阶段，有的还一直在使用经验法，虽然这些技术 对冻结的预测有很大帮助，但有时误差太大，仍然不能作为可靠的井筒施工依据。因此，急 需一种准确地获得冻结情况的预测的方法，其中确定冻结壁交圈时间和井帮温度的提前预 测是获得冻结情况的两个重要的因素，为掘砌施工提供重要指导，及时采取最为合适的措 施，仓il造最优的施工条件。
【发明内容】
[0006] 本发明的所要解决的技术问题在于不能用科学的方法提供一种预测立筒井筒的 冻结情况的方法，特别涉及不能准确预测冻结壁交圈时间和井帮温度，本发明采用以下技 术方案解决上述技术问题的：
[0007] -种预测立筒井筒冻结情况的方法，其特征在于，包括以下步骤：
[0008] (1)测量或查看冻结管的各参数；
[0009] (2)确定每天送冷量和推出单根冻结管第η天冻土发展半径公式；
[0010] (3)推出计算冻结交圈时间公式；
[0011] (4)推出冻土发展到井帮所需要的冻结天数Ν计算公式；
[0012] (5)推出受单圈孔影响的井帮温度计算公式；
[0013] (6)推出受两个或以上冻结孔影响的井帮温度的计算公式；
[0014] (7)确定距离影响系数ε的数值；
[0015] (8)将测量或查看冻结管的各参数和确定距离影响系数ε的数值代入第(5)步中的 计算公式中；
[0016] 详细地说，所述的步骤(1)需要测量或查看的参数分别为:查看盐水状态和冻结管 最大散热期的盐水温度，从而确定冻结管最大散热系数Kt、测量冻结管的外径de、立井处的 土的种类，确定相应的土比热容m和密度P、该土形成的岩层的温度扩散系数α和测量岩层 的原始温度to、冻结管距井帮的距离R b，统计冻结站前期盐水温度没达到设计值的准备时间 Τι〇
[0017] 详细地说，所述的步骤(2)中，每天的送冷量一定且为
[0018] 243idehekt,
[0019] 推出单根冻结管第η天冻土发展半径，计算公式如下：
[0020]
[0021] 其中he为岩层厚度，m;a为岩层的温度扩散系数;de为冻结管的外径，m;Kt为冻结 管最大散热系数，J/(m2 · h) ;m为土的比热容，J/kg · °C ;P为土的密度，kg/m3;t〇为岩层的原 始温度，°C; η为冻结天数。
[0022] 详细地说，所述的步骤(3)中，计算冻结交圈第二时间Τ2的计算公式如下：
[0023]
[0024] 计算冻结交圈时间Τ的计算公式如下：
[0025]
[0026] 其中Ti为统计冻结站前期盐水温度没达到设计值的准备时间；LC为每个层位的最 大成孔间距，m;a为岩层的温度扩散系数;de为冻结管的外径，m;Kt为冻结管最大散热系数， J/(m2 · h);m为土的比热容，J/kg · °C;P为土的密度，kg/m3;t〇为岩层的原始温度，°C。
[0027] 详细地说，所述的步骤(4)中，冻土发展到井帮所需要的冻结天数N的计算公式为：
[0028]
[0029]其中Rb为冻结管距井帮的距离，m;a为岩层的温度扩散系数;de为冻结管的外径， m;Kt为冻结管最大散热系数，J/(m2 · h) ;m为土的比热容，J/kg · °C ;P为土的密度，kg/m3;t〇 为岩层的原始温度，°C。
[0030]详细地说，所述的步骤(5)中，受单圈孔影响的井帮温度的计算公式为：
[0031]
[0032] 其中Rb为冻结管距井帮的距离，m;a为岩层的温度扩散系数;de为冻结管的外径， m;Kt为冻结管最大散热系数，J/(m2 · h) ;m为土的比热容，J/kg · °C ;P为土的密度，kg/m3;t〇 为岩层的原始温度，°C;n为冻结天数，该公式中的Rb>500mm，该公式适用于受单圈孔影响的 冻结壁任意点的温度。
[0033] 详细地说，所述的步骤(6)中，计算出受两个或以上冻结孔影响的井帮温度，计算 公式如下：
[0034]
[0035]其中：Rb为冻结管距井帮的距离，m;a为岩层的温度扩散系数;de为冻结管的外径， m;Kt为冻结管最大散热系数，J/(m2 · h) ;m为土的比热容，J/kg · °C ;P为土的密度，kg/m3;t〇 为岩层的原始温度，°C ;n为冻结天数，ε为距离影响系数，该公式中的Rb>500mm。
[0036] 详细地说，Rn是第η天的冻结壁发展半径，Rbmin、Rbmax是测点到最近、最远冻结管的 距离。单排冻结管rn<Rbmin时，ε取1，双排孔时ε最小值取0.4,最大值为1。
[0037]详细地说，所述的步骤(2)之前还包括以下步骤：
[0038] (1)假设原岩温度从to降到0°C为一个循环，盐水和岩层热能交换可以用如下公式 表不：
[0039]
[0040] 其中he为岩层厚度，m;de为冻结管的外径，m;Kt为冻结管最大散热系数，J/(m2 · h) ;m为土的比热容，J/kg · °C ;P为土的密度，kg/m3;t〇为岩层的原始温度，°C ;
[0041] (2)从上述公式推导出权利要求1中的单根冻结管第一天冻土发展半径的公式：
[0042]
[0043] (3)假设第一天发展半径，第二天发展到R2· · ·第η天发展到Rn;
[0044] 冻结管每天送冷量近似相等，故有：
[0045]
[0046] 式中Q表示冻结管每天的送冷量，At = t〇-〇表示每天发展的冻结壁的温度差为常 数，根据上式可得到：
[0047] ?]
1 V -L I ^ v j Z v a ? -1
[0048]由上式可得单根冻结管第η天冻土发展半径：
[0049]
[0050] 详细地说，所述的步骤(4)之前还包括以下步骤：
[0051] (1)确定冻结范围内某点下降温度，该点第η天温度下降了
[0052]
[0053] (2)确定冻结第η天后的井帮温度tb，公式为
[0054]
[0055]本发明的优点在于：
[0056] (1)本发明所得出的冻结壁温度场计算公式中只有一两个主观参数，其他全部为 冻土试验所得到的客观参数，大大减少人为主观所造成的误差，计算得到的结果与实际结 果非常接近。
[0057] (2)本发明可以满足立井冻结施工所需的绝大多数的工程数据。
[0058] (3)本发明可以优化冻结方案，在冻结方案讨论过程中，提供可靠的预测数据，进 而根据预测数据发现方案的缺陷和不足，从而提前改进冻结方案。
【附图说明】 图1为本发明的立井冻结温度场发展与分布计算模型图； 图2为本发明的冻结冻结孔布置图。
【具体实施方式】
[0059]冻结是经过电能将氨或者二氧化碳等材料液化冷却，和盐水循环系统进行热能交 换，低温盐水再和岩层进行热能交换，从而达到降低岩土温度的目的，最后形成完整的冻结 壁。
[0060]如附图1-2所示，图中，1为冻结管；2区间范围为第一天冻结壁发展位置;3代表的 区间范围冻结影响范围;4的区间范围冻结壁计算半径Rn。
[0061 ]预测立井井筒冻结情况的步骤如下：
[0062]第一步:测量或查看冻结管的各参数：
[0063] 查看盐水状态和冻结管最大散热期的盐水温度，从而确定冻结管最大散热系数 Kt、测量冻结管的外径de，查看立井处的土的种类，确定相应的该土比热容Μ和密度P，该土 形成的岩层的温度扩散系数α和测量岩层的原始温度to,冻结管距井帮的距离R b，统计冻结 站前期盐水温度没达到设计值的准备时间Tu
[0064] 第二步:获得冻结管每天的送冷量和推出单根冻结管第η天冻土发展半径公式；
[0065] (1)假设原岩温度从to降到0°C为一个循环，盐水和岩层热能交换可以用如下公式 表不：
[0066]
( 1 )
[0067]其中he为岩层厚度，m;a为岩层的温度扩散系数;de为冻结管的外径，m;Kt为冻结 管最大散热系数，J/(m2 · h) ;m为土的比热容，J/kg · °C ;P为土的密度，kg/m3;t〇为岩层的原 始温度，°C; η为冻结天数。
[0068] (2)由公式（1)可知冻结管每天送冷量一定且为
[0069] 243idehekt (2)
[0070] (3)从公式(2)可推导出单根冻结管第一天冻土发展半径的公式：
[0071]
(3)
[0072] (4)推导出单根冻结管第η天冻土发展半径的公式
[0073] 假设第一天发展半径Ri，第二天发展到R2· · ·第η天发展到Rn;
[0074] 冻结管每天送冷量近似相等，故有：
[0075](4) ^ H-1
[0076]式中Q表示冻结管每天的送冷量，At = t〇-〇表示每天发展的冻结壁的温度差为常 数，根据上式可得到：
[0077]
[0078] 由（5)式可以推出：
[0079]
[0080]由上式可得单根冻结管第η天冻土发展半径： _
(7)
[0082] 冻结壁发展是从冻结管外壁开始，逐渐向四周发展，岩土温度也是从冻结管外壁 开始慢慢降温，但降温的幅度与离冻结管距离的远近有关，离冻结管越远的地方，温度下 降的越慢，这是由冷量在冻结管周围分部不均造成的，绝大多数冷量被冻结管附近的冻土 吸收了，只有少量冷量扩散到较远的地方。通过模拟可知，在冻结扩散半径以内的冷量占总 送冷量的 60%左右，这和土的导热能力密切相关(导热能力越大，温度扩散范围越大，反之 亦然），也有一小部分冷量扩散到影响范围以外的地方，还有一部分是地热影响，被消耗掉 了（这就是为何在冻结立井开机前，要弄清楚周围地下水流情况和水源井使用情况的原 因），所以冻土扩散半径计算公式应做相应调整：
[0083]
[0084] 其中:α为温度扩散系数，η为冻结天数。
[0085] 将公式(3)带入公式(8)中，即单根冻结管第η天冻土发展半径的公式为：
[0086]
(9)
[0087] 不同岩层的温度扩散速度不一样，且能量损失大小不一样，故不同岩层温度扩散 系数见下表：
[0088]
[0089]表1温度扩散系数取值表
[0090]第三步:推出计算冻结交圈时间的公式；
[0091 ]根据以上分析可知，η为冻结天数，单个冻结孔冻土发展速度为Rn，冻结层位的冻 结孔成孔最大间距假设为Lc，显然冻结壁交圈受此最大间距控制，相邻冻结孔的冻土相交 后即为冻结交圈，故：
[0092]
[0093] 根据公式(7)推导出冻结交圈第二时间T2的计算公式：
[0094]
:(11)[0095]取等于号就是交圈时间。即：
(12)
[0096]
[0097] 由于计算式（3)中kt基于盐水温度取值是-20度，同时公式（5)也是在假设冻结站 送冷量稳定的基础上，所以根据以上公式计算出的冻结交圈时间不包括冻结站前期盐水温 度没达到设计值的时间。通常冻结站盐水温度从20度降到-20度需要一段时间，盐水降温阶 段时间为Τι;
[0098] 即冻结交圈时间Τ的计算公式为：
(13)
[0099]
[0100] 第四步:计算冻土发展到井帮所需要的冻结天数Ν
[0101] 通过模拟可知，冻结范围内的某点温度随时间变化曲线近似是一个冻结天数η的 1/2次幂的函数，第η天温度下降了 ：
[0102] （14)
[0103] 其中η代表冻结天数，Ν代表动土发展到井帮所需要的冻结天数。
[0104] 故冻结η天后的井帮温度：
[0105]
(15)
[0106] (16)
[0107] 第五步，推出受单圈孔影响的井帮温度计算公式：
[0108]
(17.)
[0109] 其中：Rb为冻结管距井帮的距离，m;a为岩层的温度扩散系数;de为冻结管的外径， m;Kt为冻结管最大散热系数，J/(m2 · h) ;m为土的比热容，J/kg · °C ;P为土的密度，kg/m3;
[0110] 公式（17)随着Rb取值不同既可以计算井帮温度也可以计算冻结壁内任意一点的 温度，但不适用于Rb<500mm，因为在冻结管500mm左右的范围内，冻土已不再继续降温了， 可以认为冻结管附近的冻土处于恒温状态，它吸收的热量全部扩散到恒温带以外的地方。 所以公式(17)中要求R b>500mm，并且该公式适用于受单圈孔影响的冻结壁任意点的温度。
[0111] 第六步:推出受两个或以上冻结孔影响的井帮温度的计算公式
[0112] 通过现场实测数据对比分析，公式（17)在计算受单圈孔影响的冻结壁任意一点温 度时是准确的，而计算两排冻结孔之间的温度有较大差异，如主排孔和辅助孔之间的测点。 原因是双排孔之间，测点同时受两个或以上冻结孔影响，影响大小受单冻结管冻土发展半 径与测点到最近冻结孔的距离比控制，定义为距离影响系数ε，公式(17)相应修改为：
[0113]
(18)
[0114] 第七步:确定距离影响系数ε的数值
[0115]
[0117] 表2:距离影响系数表
[0118] 表2中匕是第η天的冻结壁发展半径，Rbmin、Rbmax是测点到最近、最远冻结管的距离。 单排冻结管KRbmin时，ε取1，双排孔时ε最小值取0.4,最大值为1。
[0119] 第八步:将测量或查看冻结管的各参数和确定的ε的值代入公式（18)中的计算公 式中。
[0120] 具体实施例如下：
[0121] -、计算冻结交圈时间：
[0122] 3个水文孔分别为S1孔、S2孔、S3孔，S1孔、S2孔、S3孔深度分别为92m，153m、338m， 所对应的层位的最大钻孔间距分别为1.816m，2.026m，2.269m，岩性分别是粘质粉砂、粗砂、 砾石层，统计出S1孔、S2孔、S3孔盐水降温阶段时间Ti均为9天，。
[0123] 第一步，根据温度扩散系数取值表和S1孔、S2孔、S3孔的岩性特质，α取值分别为 Ο·53、0·58、0·63·
[0124] 测得S1孔盐水状态为层流，冻结管最大散热期的盐水温度零下26/°C;
[0125] 测得S2孔盐水状态为层流，冻结管最大散热期的盐水温度零下28/°C;
[0126] 测得S3孔盐水状态为层流，冻结管最大散热期的盐水温度零下30/°C;
[0127] 根据表3，Kt取值范围均为kt= (882000~1080000)J/(m2 · h);
[0128]
[0129] 表3冻结管最大散热系数(单位热流量)Kt参考值
[0130] 测得si层位其他参数分别为：t〇= 19°C，m= 1200J/kg · °C，P = 2300kg/m3，de = 0.159m,Lc= 1.816m；
[0131] 测得52层位其他参数分别为：1：() = 2〇.4<€，111=12〇〇]71^.<€，0 = 22〇〇1^/1113，(16 = 0.159m,Lc = 2.026m；
[0132] 测得S3层位其他参数分别为：t〇 = 25.8°C，m=1200J/kg · °C，P = 2600kg/m3,、de = 0· 159m，Lc = 2.269m;
[0133] 根据公式（13)可得：
[0134] SI孔的交圈时间：T = 37天；
[0135] S2孔的交圈时间：Τ = 39天；
[0136] S3孔的交圈时间:Τ = 51天。
[0137] 实际的S1孔、S2孔、S3孔的实际交圈时间为35、38、48天。
[0138] 由上式计算出来的交圈时间和水文孔实际冒水时间相差很小。有误差时间是由于 在测试过程中有参数的范围的选择，该误差也在正常范围内。
[0139] 二、计算井帮温度
[0140] 冻结孔到井帮距离Rb在200m以上为1.5m，200m以下2.1m。根据冻结孔偏斜实测数 据并结合下表2可得相应冻结层位的距离影响系数ε的取值。
[0141] 根据公式(18)综合上述数据和表4中参数确定最后一列井帮温度。
[0142]
[0143] 表4预测井帮温度对应表
[0144]
[0145] 表5实测井帮温度应对表
[0146] 根据表4和表5的对比，预测的井帮温度和实际测量的井帮温度误差小于1°C，产生 误差原因是由于测点无法准确定位和环境温度的影响所致。
[0147] 所述的步骤(1)和\或步骤(3)可以放在步骤(8)后。
[0148] 以上所述仅为本发明创造的较佳实施例而已，并不用以限制本发明创造，凡在本 发明创造的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明创造 的保护范围之内。
【主权项】
1. 一种预测立筒井筒冻结情况的方法，其特征在于，包括W下步骤： (1) 测量或查看冻结管的各参数； (2) 确定每天送冷量和推出单根冻结管第η天冻±发展半径公式； (3) 推出计算冻结交圈时间公式； (4) 推出冻±发展到井帮所需要的冻结天数的十算公式； (5) 推出受单圈孔影响的井帮溫度计算公式； (6) 推出受两个或W上冻结孔影响的井帮溫度的计算公式； (7) 确定距离影响系数ε的数值； (8) 将测量或查看冻结管的各参数和确定距离影响系数ε的数值代入第(5)步中的计算 公式中。2. 根据权利要求1所述的一种预测立筒井筒冻结情况的方法，其特征在于，所述的步骤 (1) 需要测量或查看的参数分别为:查看盐水状态和冻结管最大散热期的盐水溫度，从而确 定冻结管最大散热系数Kt，测量冻结管的外径de、立井处的±的种类，确定相应的±比热容 m和密度0、该±形成的岩层的溫度扩散系数α和测量岩层的原始溫度to、冻结管距井帮的距 离化，统计冻结站前期盐水溫度没达到设计值的准备时间Τι。3. 根据权利要求1所述的一种预测立筒井筒冻结情况的方法，其特征在于，所述的步骤 (2) 中，每天的送冷量一定且为 24JT^dehekt ? 推出单根冻结管第η天冻±发展半径，计算公式如下：其中he为岩层厚度，m;a为岩层的溫度扩散系数;de为冻结管的外径，m;Kt为冻结管最 大散热系数，J/(m2 · }1);111为±的比热容，J/kg · °(：;0为±的密度，kg/m3;t日为岩层的原始溫 度，°C;n为冻结天数。4. 根据权利要求1所述的一种预测立筒井筒冻结情况的方法，其特征在于，所述的步骤 (3) 中，计算冻结交圈第二时间T2的计算公式如下：计算冻结交圈时间T的计算公式如下：其中Τι为统计冻结站前期盐水溫度没达到设计值的准备时间；Lc为每个层位的最大成 孔间距，m;a为岩层的溫度扩散系数;de为冻结管的外径，m;Kt为冻结管最大散热系数，J/ (1112'}1);111为±的比热容，J/kg.。0为±的密度，kg/m3;t日为岩层的原始溫度，°C。5. 根据权利要求1所述的一种预测立筒井筒冻结情况的方法，其特征在于，所述的步骤 (4) 中，冻±发展到井帮所需要的冻结天数N的计算公式为：其中化为冻结管距井帮的距离，m;a为岩层的溫度扩散系数;de为冻结管的外径，m;Kt为 冻结管最大散热系数，J/(m2 · }1);111为±的比热容，J/kg · °(：;0为±的密度，kg/m3;t〇为岩层 的原始溫度，C。6. 根据权利要求1所述的一种预测立筒井筒冻结情况的方法，其特征在于，所述的步骤 (5) 中，受单圈孔影响的井帮溫度的计算公式为：其中化为冻结管距井帮的距离，m;a为岩层的溫度扩散系数;de为冻结管的外径，m;Kt为 冻结管最大散热系数，J/(m2 · }1);111为±的比热容，J/kg · °(：;0为±的密度，kg/m3;t〇为岩层 的原始溫度，°C;n为冻结天数，该公式中的化巧00mm，该公式适用于受单圈孔影响的冻结壁 任意点的溫度。7. 根据权利要求1所述的一种预测立筒井筒冻结情况的方法，其特征在于，所述的步骤 (6) 中，计算出受两个或W上冻结孔影响的井帮溫度，计算公式如下：其中：Rb为冻结管距井帮的距离，m;a为岩层的溫度扩散系数;de为冻结管的外径，m;Kt 为冻结管最大散热系数，J/(m2 · }1);111为±的比热容，J/kg · °(：;0为±的密度，kg/m3;t日为岩 层的原始溫度，°C ;η为冻结天数，ε为距离影响系数，该公式中的化巧00mm。8. 根据权利要求7所述的一种预测立筒井筒冻结情况的方法，其特征在于，Rn是第η天的 冻结壁发展半径，Rbmin、Rbmax是测点到最近、最远冻结管的距离。单排冻结管rn《Rbmin时，ε取 1，双排孔时ε最小值取0.4,最大值为1。9. 根据权利要求3所述的一种预测立筒井筒冻结情况的方法，其特征在于，所述的步骤 (2)之前还包括W下步骤： (1) 假设原岩溫度从to降到(TC为一个循环，盐水和岩层热能交换可W用如下公式表示：其中he为岩层厚度，m;de为冻结管的外径，m;Kt为冻结管最大散热系数，J/(m2 · h);m为 ±的比热容，J/kg · °(：;0为±的密度，kg/m3;t〇为岩层的原始溫度，°C; (2) 从上述公式推导出权利要求1中的单根冻结管第一天冻±发展半径的公式：(3) 假设第一天发展半径化，第二天发展到R2· ··第η天发展到Rn; 冻结管每天送冷量近似相等，故有：式中Q表示冻结管每天的送冷量，At = t〇-〇表示每天发展的冻结壁的溫度差为常数，根 据上式可得到：由上式可得单根冻结管第η天冻±发展半径：10.根据权利要求5所述的一种预测立筒井筒冻结情况的方法，其特征在于，所述的步 骤(4)之前还包括W下步骤： (1) 确定冻结范围内某点下降溫度，该点第η天溫度下降了(2) 确定冻结第η天后的井帮溫度tb，公式为
【文档编号】G06Q10/04GK105976060SQ201610305516
【公开日】2016年9月28日
【申请日】2016年5月9日
【发明人】张海骄, 姚直书, 汤江明, 王成博, 常军, 魏玮, 杨海若, 李世辉, 谢玉兵
【申请人】淮南矿业(集团)有限责任公司