AQ1067-2008矿井风流势力状态预测方法.pdf

ICS D 备案号备案号 中华人民共和国安全生产行业标准 中华人民共和国安全生产行业标准 送审稿 发 布 A AQQ/ /T T 矿井风流热力状态预测方法 Mine Merry Heat Prediction Standard AQ/T10672008 20081119 发布 20090101 实施 国家安全生产监督管理总局国家安全生产监督管理总局 AQ/T10672008 2 目 录 前 言 Ⅱ 1 范 围 1 2 规范性引用文 件 1 3 术语和定 义 1 4 矿井风流热力状态预测的基础资 料 2 5 矿井风流热力状态预测内容、范围及精度要 求 . . 4 6 矿井风流热力状态预测方 法 . . 5 附录A规范性附录井底车场风流热状态预测方 法 6 附录B规范性附录 井巷风流热状态预测方 法 . . . . . . 8 附录C规范性附录 掘进工作面风流热状态预测方 法..12 附录D规范性附录 回采工作面风流热状态预测方 法..14 附录E资料性附录 大型机电设备放热量计 算 . . 1 6 附录F资料性附录 矿井微气候基本参数计 算 . . 1 8 AQ/T10672008 3 前 言 在深部采矿工程中，矿井热环境条件将是制约采矿生产、危及安全和健康的一个重要因素。 在新矿井、新水平和新采区的设计中，为了调查清楚投产后的矿井热环境状况，则必须对矿井风 流的热力状态进行预测，以为合理地节能、经济、安全和环保设计矿井开拓系统、生产系统和 生产工艺、通风系统与降温系统提供必要的热力参数，规范矿井风流热力状态预测方法，为此制 定本标准。 由于矿井热环境条件的复杂性和多变性，致使预测精度受到限制，因此，在选用矿井制冷降 温设备时，要有足够的备用量。进行矿井风流热力状态预测，除执行本规范的规定外，尚应符合 国家现行的有关标准、规范和规定。 本标准为推荐性标准。 本标准的附录 A、附录 B、附录 C、附录 D 为规范性附录，附录 E、附录 F 为资料性附录。 本标准由国家安全生产监督管理总局提出。 本标准由全国安全生产标准化技术委员会煤矿安全分技术委员会归口。 本标准负责起草单位中国煤炭工业劳动保护科学技术委员会。 本标准参加起草单位煤炭科学研究总院抚顺分院。 本标准主要起草人李红阳、杨德源、罗海珠、朱锦文 AQ/T10672008 4 矿井风流热力状态预测方法 1 范围 本标准规定了矿井风流热力状态预测的定义术语、预测的基础资料、内容、范围、精度要求 及预测方法。 本标准适用于全国井工开采的矿井，包括生产、新建和改、扩建矿井。 2 规范性引用文件 下列文件中的条款通过本标准的引用而成为本标准的条款。凡是注明日期的引用文件，其随 后所有的修改件不包括堪误的内容或修订版均不适用于本标准，然而，鼓励根据本标准达成协 议的各方研究后确定，是否可使用这些文件的最新版本。凡是不注日期的引用文件，其最新版本 适用于本标准。 GB50155-1992 采暖通风与空气调节术语标准 GB50215-2005 煤炭工业矿井设计规范 DZ/T0080-1993 煤田地球物理测井规范 DZ/T0215-2002 煤、泥炭地质勘查规范 煤矿安全规程2006 年版 3 术语和定义 GB50155-1992 中确立的术语和定义以及下列术语和定义适用于本标准。 3.1 矿井热源 mine heat source 处在井巷中的，并能对风流加热或吸热的载热体，称为矿井热源。 3.2 矿山地热 geothermal Mine 在矿区开采范围内的地热状况，称为矿山地热。 3.3 地温场 geothermal Field 地层中的温度分布，称为地温场。 3.4 地温异常区 temperature anomalies 地温梯度值，高于或低于正常值1.63.0/100m℃的地层区域称为地温异常区。 3.5 地面大气参数 ground atmospheric parameters 矿区地面空气的参数，包括温度、湿度、风速和气压等，可由矿区附近气象站台查出。 3.6 矿井风流热力状态预测 mine Merry heat prediction 在新矿井或新水平、新采区设计中，根据矿井的热环境条件，预先计算出基建和投产后矿内 风流的热力状态。 3.7 最热月平均气象参数 on average the hottest meteorological parameters 最热月空气的温度、相对湿度和大气压力的平均值。 AQ/T10672008 5 3.8 正向、逆向热力计算 positive, the reverse thermodynamic calculation 沿风流方向进行的热力计算称为正向热力计算，逆风流方向进行的热力计算称为逆向热力计 算。 3.9 矿内作业场所温度 mine temperature at work area 在矿井下从事掘进、回采、巷道维护以及设备操作等场所的风流温度即干球温度。 4 矿井风流热力状态预测的基础资料 4.1 矿区气象资料 调查矿区附近气象站台近 10 年的气温、相对湿度和大气压力历年的月平均值和年平 均值； 03.2m 之间不同深度的温度； 年平均日影响深度； 降水量； 气象站的地理坐标； 4.2 井田地质地温资料 提供井田的精查地质报告； 井田恒温带的深度，温度，地温梯度； 矿井不同开采深度的地温深钻孔或浅钻孔测温，等温线图，地温场分布； 井田地质构造，岩浆活动以及岩层中放射性元素含量； 井田地层和井巷围岩种类、分布及岩石的导热率、密度和比热； 4.3 矿井水文地质资料 地下水的来源、分布、通道，水量及温度； 矿井涌水量，水温及水质； 矿井排水动力及排水系统； 水沟断面规格； 顶板及两帮淋水情况，水温及水量； 4.4 地温、岩石热物理参数测试及深部地温预测 参照 DZ/T0080-1993 煤田地球物理测井规范、DZ/T0215-2002 煤、泥炭地质勘查规范； 4.5 矿井生产技术资料 4.5.1 矿井开拓系统 开拓系统布置； 矿井生产能力日产量及年产量； 矿井服务年限； 4.5.2 矿井通风系统 风量分配；矿井、采区及采面的通风方式； 矿井主要通风机的功率； AQ/T10672008 6 掘进通风方式及局部通风机的功率。 4.5.3 矿井运输系统 运输系统布置； 运输方式； 运输动力的容量； 运输量； 日运输作业时间； 4.5.4 矿井压缩空气供应系统； 4.5.5 井下冷、热水管道的分布； 4.5.6 矿井安全情况瓦斯，火灾，粉尘； 4.5.7 矿井供电系统 中央变电所、采区变电所及移动变电站的位置和设备容量； 主要机电设备的分布及容量； 4.6 矿井巷道状况 巷道编号、名称； 井巷的几何尺寸断面、长度及周长； 始、终端的标高； 围岩岩性； 围岩初始温度； 井巷通风时间； 断面形状及支护方式； 风量； 4.7 采掘工作面状况 4.7.1 回采工作面状况 回采工作面编号、名称； 采面上、下口的标高； 采面长度； 采高； 控顶距； 顶、底板岩性； 围岩初始温度； 日产量； 采面供风量，采面通风方式； 采面内及上、下平巷机电设备的容量kW； 采矿方法，回采工艺及顶板管理； 采面内的涌水量及水温； AQ/T10672008 7 采面内每斑的工作人数。 4.7.2 掘进工作面状况 掘进工作面编号、名称； 掘进工作面的几何尺寸断面，长度，周长及标高； 围岩岩性； 围岩初始温度； 通风方式，局部通风机功率，供风量，送风距离，风筒直径及材料； 日进尺； 每班工作人数； 掘进及转载设备的容量； 巷道支护方式； 涌水量及水温； 每斑火药耗量。 5 矿井风流热力状态预测内容、范围及精度要求 5.1 预测内容 预测采掘工作面最热月的风流热力状态参数。 5.1.1 预测矿井通风系统特征点井底车场风流热力状态参数； 5.1.2 预测矿井通风系统特征点井巷风流热力状态参数； 5.1.3预测回采工作面出口距回风道口15m处采面内的运输空间中的风流热力状态参数参考 煤 矿安全规程1992 年版执行说明； 5.1.4 进行逆向热力计算时，预测回采工作面进口距进风巷 3m 的采面内的风流热力状态参数； 5.1.5 预测掘进通风系统中特征点风机后、风筒出口及回风巷的风流热力状态参数； 5.1.6 预测掘进工作面距迎头 2m 处回风流风流的热力状态参数。 5.2 预测范围 对于新矿井、新水平、新采区设计，要分期基建期、投产期、达产期、中期和后期预测矿 井风流的热力状态参数。 5.3 预测精度 特征点井底车场、井巷、掘进工作面及回采工作面风流温度预测绝对误差≤1.5℃，相对误 差不得超过5。 5.4 对风流热力状态参数测试仪器仪表的精度要求 5.4.1 温度测量测量范围为-2550℃，最小分度值为 0.1℃； 5.4.2 大气压力测量测量范围为 80kPa～120kPa600～900mmHg，最小分度值为 50Pa； 5.4.3 风速测量测量范围低速为 0.2～5m/s，中速为 0.4～10m/s,高速为 0.8～30m/s；测量精度 要求低速为≤0.2m/s，中速为≤0.4m/s，高速为≤0.8m/s； 5.4.4 相对湿度测量测量范围为 50～100，测量精度要求为≤5。 AQ/T10672008 8 6 矿井风流热力状态预测方法 采用单线路预测法主干线路法， 从已知节点一般以入风井筒的地面井口进风点开始， 顺风 流方向，选择一条通往预测终点一般为回采工作面出口的主要风路，逐段预测，计算出相关特 征点如井底车场、一翼或采区进风点以及采掘工作面进风点等的风流热状态参数。 6.1 预测方法的选用 6.1.1 按照 GB50215-2005 标准 7.5.2 的规定，所采用的预测方法必须经过鉴定的、并经多个矿井 的预测验证，证明是可靠的； 6.1.2 必须调查清楚预测矿井的主要热、湿源和设计或实际通风状况； 6.1.3 在矿井通风网路中的风流汇合点，应考虑汇合风流的影响，计算混合风流的热、湿参数； 6.1.4 在通风网路解算中，应考虑风流热力状态的变化； 6.2 井底车场风流热状态预测方法按照规范性附录 A 进行。 6.3 井巷风流热状态预测方法按照规范性附录 B 进行。 6.4 掘进工作面风流热状态预测方法按照规范性附录 C 进行。 6.5 回采工作面风流热状态预测方法按照规范性附录 D 进行。 6.6 大型机电设备放热量计算方法按照资料性附录 E 进行。 6.7 矿井微气候基本参数计算按照资料性附录 F 进行。 6.8 预测结果处理 预测完成后，应编制矿井风流热力状态预测报告。 预测报告主要内容包括 预测矿井名称、预测时间、预测人员； 预测目的和要求； 当时矿井通风和生产情况； 列出矿井风流热力状态预测的基础资料； 矿井风流热力状态参数计算； 预测结果分析； 绘制预测矿内风流热状态曲线图； 提出矿井热害防治措施的建议方案。 AQ/T10672008 9 附录 A 规范性附录 井底车场风流热状态预测方法 A.1 井底车场风流热状态预测方法 A.1.1 井筒深度 H900m 时 井底车场风流温度 t2 mk cnk 2 1 mk lkHg10i mk cnk 4 1 t 22B p22B 22B 22B 3 1 2 22B p22B 2           ，℃..A-1 A.1.2 井筒深度 H≥900m 时 井底车场风流温度 t2 A 2 1 DA 4 1 t 2 2 ，℃...A-2 mk M HDk 2 1 cnk A 2B2 B p2B2      mk lkt 2 1 t M HDk M Q gH10i D 2B2 2B21gu BB M3 1        式中 t1地面井口入风温度，℃； t2井底车场风流温度，℃； i1地面井口入风焓值，kJ/kg； i2井底车场风流焓值，kJ/kg；  M Q井筒中各种局部热源放热量之和，kW； Qgu井筒围岩放热量，kW； 1 地面井口入风相对湿度，； 2 井底车场风流相对湿度，； kB1地面井口气压修正系数， 1 B1 B 08.101 k； kB2井底车场气压修正系数， 2 B2 B 08.101 k； B1地面井口大气压力，kPa； B2井底车场大气压力，kPa； tgu井筒围岩的平均初始温度，℃； MB通过井筒的风量，kg/s； kτ风流与围岩间的不稳定热交换系数，kW/m2℃； H井筒的深度，m； D井筒的直径，m； AQ/T10672008 10 cp空气的定压比热，cp1.005kJ/kg℃； γ水在 0℃时的汽化潜热，γ2.501kJ/g℃； m，n，l常系数，见表 A1。 表 A1 常系数 m，n，l 风温 ℃ 010 515 1020 1525 2030 2535 m 3.77 3.96 4.68 6.42 9.82 16.60 n 0.266 0.209 0.089 -0.114 -0.424 -0.940 l 0.0120 0.0158 0.0206 0.0264 0.0334 0.0432 AQ/T10672008 11 附录 B 规范性附录 井巷风流热状态预测方法 B.1 水平井巷风流热状态预测方法 B.1.1 通风时间小于 1 年的水平井巷 B.1.1.1 水平井巷终端温度 cM Q T e1 ett pB M 12        ，℃B-1   E cM LBk cM LUk cM ULk pB WW pB TT pB   F-t cM LBk t cM LUk t cM ULk T w pB WW T pB TT gu pB 12 ； m y PB N E  ； EF                0 E E1 E 1 ln 0 E1 1 1 1 B.1.1.2 水平井巷始端温度逆向计算 cM Q T 1e ett pB M 21        ，℃..B-2 B.1.2 通风时间大于 1 年的水平井巷 B.1.2.1 水平井巷终端温度 A 2 1 DA 4 1 t 2 2 ...B-3 lk k 2 1 cnk A 2B 1p2B    lk tk 2 1 mk M Q ki D 2B 112B B M 21     B WW B TT B 1 M LBk M LUk M ULk k  W B WW T B TT gu B 2 t M LBk t M LUk t M ULk k  B.1.2.2 水平井巷始端温度逆向计算 AQ/T10672008 12 A 2 1 DA 4 1 t 2 1 ，℃B-4 lk k 2 1 cnk A 1B 1p1B    lk tk 2 1 mk M Q ki D 1B 211B B M 22     B.2 倾斜井巷风流热状态预测方法 B.2.1 通风时间小于 1 年的倾斜井巷 B.2.1.1 倾斜井巷终端温度 cM Q YT e1 ett pB M 12        ，℃B-5   E cM LBk cM LUk cM ULk pB WW pB TT pB   F-t cM LBk t cM LUk t cM ULk T w pB WW T pB TT gu pB        sin 1 e1 1 L cM ULk Y - pB 12 ； m y PB N E  ；EF                0 E E1 E 1 ln 0 E1 1 1 1 B.2.1.2 倾斜井巷始端温度逆向计算 cM Q YT 1e ett pB M 21        ，℃B-6 B.2.2 通风时间大于 1 年的倾斜井巷 B.2.2.1 倾斜井巷终端温度 A 2 1 DA 4 1 t 2 2 ，℃....B-7 lk k 2 1 cnk A 2B2 1p2B2    AQ/T10672008 13 lk tk 2 1 mk M Q kHg10i D 2B2 112B2 B M 2 3 1      B WW B TT B 1 M LBk M LUk M ULk k  W B WW T B TT gu B 2 t M LBk t M LUk t M ULk k  B.1.2.2 井巷始端温度逆向计算 A 2 1 DA 4 1 t 2 1 ，℃B-8 lk k 2 1 cnk A 1B1 1p1B1    lk tk 2 1 mk M Q kHg10i D 1B1 211B1 B M 2 3 2       式中 t1巷道始端风流温度，℃； t2巷道终端风流温度，℃； tgu围岩原始温度，℃； MB风量，kg/s； kτ风流与围岩不稳定热交换系数，kW/m2℃； kT水管壁传热系数，kW/m2℃； kW水沟盖板传热系数，kW/m2℃； U巷道周长，m； UT水管周长，m； BW水沟的宽度，m； L巷道长度，m； tT水管内水的温度，℃ tw水沟中水的平均温度，℃ ∑QM各种局部热源放热量之和，W； σ地温梯度，℃/m； ψ巷道倾角，； φ1巷道始端风流的相对湿度，； φ2巷道终端风流的相对湿度，； kB巷道风流气压修正系数， B 08.101 kB； kB1巷道始端风流气压修正系数， 1 B1 B 08.101 k； AQ/T10672008 14 kB2巷道始端风流气压修正系数， 2 B2 B 08.101 k； B巷道大气压力，kPa； B1巷道始端大气压力，kPa； B2巷道终端大气压力，kPa； ε’，Pm，Ny常系数，见表 B1。 表 B1 空气温度 ℃ 110 1117 1723 2329 2935 3545 4550 ε’ 9.324 19.979 -3.770 -8.988 -14.288 -22.958 -27.685 井下 1016.20 1459.012108.053028.41 4281.276497.05 10193 Pm Pa 地面 734.16 1053.361522.082187.85 3105.554692.24 6974 Ny 95880 77805 223307306100 415160 608076 850589 AQ/T10672008 15 附录 C 规范性附录 掘进工作面风流热状态预测方法 C.1 掘进工作面风流温度预测图见图 E1 图 E1 掘进工作面风流预测系统 C.2 风机后风流温度特征点 2 B e t112 V N ktttt， ℃..C-1 C.3 风筒出口温度特征点 3 ]NRAAR[NN1 AR ]ARAA[N]ANNA[tR t 1ba2a1 c1baaa21 3   ， ℃C-2 M c nk 1R p 4B 1   ；M c nk 1R p 3B 2       p 34B M3gu c mk QZ2Mt NE 2 1 c nk 1A p 5B a    NE 2 1 c nk 1A p 4B b    p 45B PBm MB guc c mk cM Q EtA    PBm BB CM Fk E  ； PBm ff c2M Fk N； B2 B3 M M k； k1 2kM M Bm B3  ； PBmc 2kM k1 Z   MM 2 1 M B3B2Bm ； PBm 33 c4kM Fkk1 M    C.4 掘进迎头温度特征点 4 1 32 4 R tR t  ，℃..C-3 C.5 回风巷中温度特征点 5 AQ/T10672008 16 a c324b 5 A Att NtA t  ，℃..C-4 式中 t1特征点 1 风流温度，℃； t2特征点 2 风流温度，℃； t3特征点 3 风流温度，℃； t4特征点 4 风流温度，℃； t5特征点 5 风流温度，℃； Ne风机的额定功率，kW； VB通过风机的风量，m3/min； kt温度系数，可取 3034； ∑QM3掘进工作面近区各种热源放热量之和，kW； ∑QM掘进巷道中各种局部热源放热量之和，kW； Ff风筒的传热面积，m2； kf风筒的传热系数，kW/m2℃； F4掘进面近区传热面积，m2； kτ4掘进面近区不稳定热交换系数，kW/m2℃； FB掘进巷道的传热面积，m2； kτB掘进巷道围岩不稳定热交换系数，kW/m2℃； MB2特征点 2 风量，kg/s； MB3特征点 3 风量，kg/s； φ4特征点 4 的相对湿度，； φ5特征点 5 的相对湿度，； kB4特征点 4 气压修正系数， 4 B4 B 08.101 k kB5特征点 5 气压修正系数， 5 B5 B 08.101 k x4特征点 4 空气的含湿量，kg/kg干空气； x5特征点 5 空气的含湿量，kg/kg干空气； m，n常系数，其取值见表 C1 表 C1 空气温度 ℃ 16 612 12181421212626313137 3745 m 3.68 2.83 0.14 -1.70 -8.33 -18.80-35.00 -63.85 n 0.34 0.48 0.69 0.82 1.13 1.53 2.05 2.83 AQ/T10672008 17 附录 D 规范性附录 回采工作面风流热状态预测方法回采工作面风流热状态预测方法 D.1 回采工作面风流温度预测图图 F1。 图 F1 回采工作面风流预测图 D.2 回采工作面出口风流温度 cM Q T e1 ett pB M 12        ， ℃.. D-1   E cM ULk pB ； 12 ； tt 0024L. 0t fmk 8 . 0 k  m y PB N E  ；    F cM tcG cM ULtk T pB kkk pB gu ； EF                0 E E1 E 1 ln 0 E1 1 1 1 D.3 回采工作面进口风流温度逆向计算 cM Q T 1e tet pB M 21        ，℃. D-2 式中 t1特征点 1 风流温度，℃； t2特征点 2 风流温度，℃； Gk运煤量，kg/s； ck煤的比热容，kJ/kgk； U回采工作面周长，m； AQ/T10672008 18 L回采工作面长度，m； MB回采工作面风量，kg/s； kτ回采工作面的不稳定热交换系数，kW/m2℃； ∑QM采面内各种热源放热量之和，kW； tfm运输线路上风流的平均湿球温度，℃； tk回采工作面落煤温度，其取值见表 D1，℃； 表 D1 tgu ℃ ≤40 ≤50 ≤60 未降温 tgu-2 tgu-4 tgu-6 tk ℃ 降温 tgu-6 tgu-8 tgu-10 AQ/T10672008 19 附录 E 资料性附录 大型机电设备放热量计算方法 E.1 采掘机械运转时放热 cjcjcj N0.8kQ， kW..E-1 式中 cj k设备的时间利用系数，它等于每日实际工作时间时被 24 除； cj N采掘机械电机消耗的功率，kW； 采掘机械设备放热量而引起的风流温升按下式计算 Bp cj cj Mc N 15k. 0t...E-2 E.2 提升设备工作时放热量 ttt Nk1 Q，kW.....E-3 式中 t 提升机工作效率； t k提升机时间利用系数； t N设备功率，kW。 E.3 变压器工作时放热量 bbb NkQ，kW....E-4 式中 05. 0kb矿井下变压器的平均热损率； b N变压器功率，kW； E.4 水泵工作时放热量 shsh 2879N. 0Q，kW......E-5 或 sh shsh sh VH 2879. 0Q  ， kW...E-6 引起水的温升按下式计算 shshsh sh c H 2879. 0t  ，℃..E-7 式中 Hsh水泵扬程，kPa； Vsh流量，m3/s； sh 密度，kg/m3； sh c比热，kJ/kg℃； sh 水泵效率，一般取 0.60.85。 E.5 通风机工作时放热量 AQ/T10672008 20 BB 564N. 0Q，kW...E-8 通风机引起的风流温升为 B B B V N 3428t，℃...E-9 式中 B N通风机功率，kW； B V通风机的风量，m3/min。 E.6 电机车工作时放热量 ece kLA 1 Q  ， kW..E-10 式中 Τ每日运输时间，h； L运输距离，km； Ac运输量，t/d； ke吨公里能耗量。 E.7 电动机运转时放热量 tEiE k1 NQ， kW.E-11 式中 Ni电动机功率，kW； ηE电动机效率； kt时间利用系数。 E.8 蓄电池机车运行时的放热量 xtxx NkkQ，kW.E-12 式中 kx修正系数，可取 1.5； kt时间利用系数； Nx设备功率，kW。 AQ/T10672008 21 附录 F 资料性附录 矿井微气候基本参数计算 F.1 矿内空气温度 F.1.1 热力学温标T，K和摄氏温标t，℃ t15.273T，KF-1 F.1.2 干球温度tB和湿球温度tf，可用干、湿球温度计直接测出； F.1.3 露点温度te x275.17 3x.273 te  ，℃ ...F-2 式中 Pw水蒸气分压力，kPa； F.2 空气压力 F.2.1 矿井下大气压力B HkBB p0 或HgBB 0 ，kPa F-3 式中 B0矿井地面大气压力，kPa； kp压力梯度，kp0.01133夏季～0.01267冬季； H计算点水平距地表的深度，m； g重力加速度，m/s2； ρ密度，kg/m3。 F.2.2 水蒸气分压力Pw gw pBp，kPa.F-4 式中 Pg湿空气中干空气的分压力kPa。 F.2.3 饱和水蒸气分压力Pb Pb值可按空气温度值t在湿空气性质表中查出，也可以用下式计算 t3 .237 t275.17 exp6105. 0Pb  ，kPa.F-5 F.3 空气的密度 F.3.1 干空气的密度 T PB 00348. 0 w g  ，kg/m3 F-6 F.3.2 湿空气的密度 T P 00132. 0 T B 00348. 0 w B ，kg/m3 F-7 H.4 空气的湿度 H.4.1 相对湿度φ AQ/T10672008 22 b w P P ， ..F-8 H.4.2 含湿量x b b w w PB P 622. 0 PB P 622. 0 x     ，kg/kg干空气..F-9 H.5 空气的焓 i x84t. 12501005t. 1i，kJ/kg ....F-10 或ltntmktci 2 Bp ，kJ/kg.....F-11 式中 t空气温度，℃； x空气含湿量，kg/kg干空气 φ相对湿度，； cp空气的定压比热，cp1.005kJ/kg℃； γ水在 0℃时的汽化潜热，γ2.501kJ/g℃； kB气压修正系数； m，n，l常系数，取值见表 A1。