1108 化工学报第 57 卷温磷蒸汽及其聚合物对燃烧塔换热器等生产设备具有极强的腐蚀性, 且这种腐蚀性随温度升高而加剧. 工艺上为缓和这种腐蚀破坏, 通常都使用以夹套结构燃磷塔为核心装备的二步法工艺流程, 具体参见图 1 图. 夹套进口水温约为 30~35, 出口水温约为 50~55, 这样低

Size: px

Start display at page:

Download "1108 化工学报第 57 卷温磷蒸汽及其聚合物对燃烧塔换热器等生产设备具有极强的腐蚀性, 且这种腐蚀性随温度升高而加剧. 工艺上为缓和这种腐蚀破坏, 通常都使用以夹套结构燃磷塔为核心装备的二步法工艺流程, 具体参见图 1 图. 夹套进口水温约为 30~35, 出口水温约为 50~55, 这样低"

青壤宥马
7 years ago
Views:

1 第 57 卷第 5 期化工学报 Vol.57 No 年 5 月 Journal of Chemical Industry and Engineering China) May 006 檭檭檭檭檭檭檭檭檭檭檭檭檭檭研究论文膜式水冷壁结构燃磷塔的开发与传质传热的数值模拟 Ⅰ) 模型建立刘宝庆, 蒋家羚浙江大学化学工程与生物工程学系, 浙江杭州 31007) 摘要 : 针对热法磷酸生产中一方面要求外部供热, 另一方面又浪费自身燃烧热的现状, 提出能回收余热的新工艺, 开发成功膜式水冷壁结构新型燃磷塔. 建立了考虑新型燃磷塔内流动燃烧传热过程的三维数学模型, 通过数值模拟得到塔内湍流场温度场和组分浓度场的分布. 模拟结果显示, 燃磷塔内存在明显的涡流射流结构, 磷的燃烧主要在喷磷枪雾化角范围内, 塔内最高燃烧温度约为 1790K, 出现在喷磷枪轴线上偏离燃磷塔中心的一侧. 烟气出口温度约为 970K, 与测量值符合较好. 关键词 : 燃磷塔 ; 节能 ; 数学模型 ; 数值模拟中图分类号 :TQ16.3;TK11 5 文献标识码 :A 文章编号 : ) 犛犻犿狌犾犪狋犻狅狀狅犳犿犪狊犪狀犱犺犲犪狋狋犪狀狊犳犲犻狀狀犲狑狆犺狅狊狆犺狅犻犮犪犮犻犱狋狅狑犲狑犻狋犺犮犻犮狌犾犪狑犪狋犲狑犪犾 Ⅰ) Modeldevelopment 犔犐犝犅犪狅狇犻狀犵, 犑犐犃犖犌犑犻犪犾犻狀犵犇犲狆犪狋犿犲狀狋狅犳犆犺犲犿犻犮犪犾犪狀犱犅犻狅犮犺犲犿犻犮犪犾犈狀犵犻狀犲犲犻狀犵, 犣犺犲犼犻犪狀犵犝狀犻狏犲狊犻狋狔, 犎犪狀犵狕犺狅狌 31007, 犣犺犲犼犻犪狀犵, 犆犺犻狀犪 ) 犃犫狊狋犪犮狋 :Toavoidtheenergywasteintheproductionoffurnacegradephosphoricacid,anewphosphoric acidtowerwithcircularwaterwalwasdesigned,whichcouldrecoverwasteheattoproducesuperheated steam.athreedimensionalmathematicalmodelofphosphoruscombustioninthenew phosphoricacid towerwasdevelopedbyincorporatingthe 犽 εturbulence model,eddydissipationturbulentcombustion modelandp1radiationheattransfermodel.thestreamlinefielddistribution,temperaturedistributionand speciesconcentration distribution were predicted by numericalsimulation. The results showed that phosphoruscombustionmainlytookplaceintheinjectionzoneofatomizedphosphorus,andthemaximum combustiontemperatureandexittemperaturewere1790kand970krespectivelywhichconformedwelto practicaldata. 犓犲狔狑狅犱狊 :phosphoricacidtower;energyconservation;mathematicalmodel;numericalsimulation 引言热法磷酸的生产主要包括矿石冶炼黄磷燃烧及五氧化二磷的水化三个独立工艺流程. 其中燃烧流程中, 固态黄磷首先在熔磷槽中被 0.~0.5 MPa 的水蒸气加热熔化至液态, 然后熔融磷经喷嘴雾化, 与由空气吸入管进入的空气充分混合, 在燃烧室剧烈燃烧生成 P O 5 烟气. 黄磷燃烧可产生大量热能, 在无任何冷却的情况下, 即使过剩空气系数达 1.5, 混合烟气温度仍高达 1800 以上. 高收到初稿, 收到修改稿. 联系人 : 蒋家羚. 第一作者 : 刘宝庆 1978 ), 男, 博士, 讲师. 犚犲犮犲犻狏犲犱犱犪狋犲 : 犆狅犲狊狆狅狀犱犻狀犵犪狌狋犺狅 :LIUBaoqing. 犈 - 犿犪犻犾 :baoqingliu@16.com

2 1108 化工学报第 57 卷温磷蒸汽及其聚合物对燃烧塔换热器等生产设备具有极强的腐蚀性, 且这种腐蚀性随温度升高而加剧. 工艺上为缓和这种腐蚀破坏, 通常都使用以夹套结构燃磷塔为核心装备的二步法工艺流程, 具体参见图 1 图. 夹套进口水温约为 30~35, 出口水温约为 50~55, 这样低温的余热在工业上很难再利用. 同时为保证黄磷的流动及雾化燃烧, 必须使其处于熔融状态, 需配置专用燃煤锅炉来供给蒸汽. 为冷却循环水, 还必须配备风机泵换热器等一系列装置. 这种一方面要求外部供热, 另一方面又浪费自身燃烧热的现象是目前热法磷酸生产中普遍存在的问题. 因此从节能与环保角度考虑, 有必要对现有工艺进行改进, 改进的关键在于开发具有余热回收功效的新型燃磷塔. Fig.1 Twostepprocessoffurnace gradephosphoricacid 1 steamboiler; meltingphosphorustrough; 3 aircompressor;4 phosphoricacidtower; 5 watercoolingsystem;6 hydrator 1 改进的热法磷酸生产工艺与装备基本目标 : 回收磷的燃烧热用于固态黄磷的熔化, 实现能源综合利用. 11 改进工艺改进的热法磷酸燃烧水合二步法工艺流程, 如 [4] 图 3 所示. 固态单质磷在熔磷槽 1 中被高温水蒸气熔化至液态, 然后, 经压缩机 5 送入的空气雾化后通过喷磷枪 4 进入燃烧塔 ; 单质磷在燃烧塔中迅速燃烧生成五氧化二磷, 同时释放出大量热能. 这些热能被布置于四周的膜式水冷壁吸收, 生成大量蒸汽. 其中一部分蒸汽被送入熔磷槽 1, 用来提供固态磷熔化所需的热量, 大部分则通过外供蒸汽管道输出, 满足其他生产生活需要 ; 最后经吸热降温后的五氧化二磷混合烟气输出送入水化塔 3, 被水吸收后生成高浓度的磷酸. 1 新型燃磷塔的结构与工作原理燃磷塔采用了燃烧室和余热回收装置合二为一的整体结构设计, 具体参见图 4 [4]. 其燃烧室周壁布置有相应数量的抗磷酸腐蚀的不锈钢管, 这些管束与连于各管间的强化换热翅片围成环形膜式水冷壁. 水冷壁两端分别布有上下集箱, 外侧设有保温层, 共同构成封闭壳体, 完成磷的燃烧与热量回收. 同常规燃磷塔比较, 新型燃磷塔具有以下特点 : 1) 新塔的燃烧室与老塔形状相同, 体积相近 ; 燃磷塔结构形状及燃烧器与燃烧室的相对安装位置不变. 因此, 新塔能保证完成磷的燃烧, 并产生足量的 P O 5 烟气 ; ) 将水冷夹套换成膜式水冷壁后, 实现了对热法磷酸生产过程的余热回收 ; 3) 整个燃磷塔可由裙座支撑或用钢架悬吊. 工作 Fig. Structureofconventionalphosphoricacidtower numbersindicatecorresponding relationofcircularwater) 有关这方面的研究国外起步较早, 日本美国和德国都针对热法磷酸生产的余热回收进行了研究开发, 并取得一些成果 [13]. 作为最大的热法磷酸生产和需求国, 开发具有独立知识产权符合我国热法磷酸生产实际的余热利用技术, 十分必要. Fig.3 Improvedprocessoffurnacegradephosphoricacid 1 meltingphosphorustrough; phosphoricacidtower; 3 hydrator;4 phosphorusfeedsystem; 5 aircompressor;6 drum;7 phosphoricacidexit

3 第 5 期刘宝庆等 : 膜式水冷壁结构燃磷塔的开发与传质传热的数值模拟 Ⅰ) 1109 时炉体的热膨胀可向上或向下自由伸展 ; 集汽管两端与汽包及上集箱的连接, 下降管两端与汽包及下集箱的连接都设计有弯头, 可作为热补偿器, 自动补偿各部件热膨胀量不一致的差别. 行了内部传质传热的三维数值模拟. 1 物理模型几何形状和结构尺寸如图 5 所示. 喷磷枪位于塔高.1m 处, 和燃磷塔轴线夹角为 38. 烟气出口截面中心线在距塔底 4.95m 的位置, 与燃磷塔轴线成 30 角. 塔径取为 1.9 m, 考虑了燃烧产物在塔壁的结膜. Fig.4 Structureofnewphosphoricacidtower 1 drum; uppermanifold;3 burnernozzle; 4 lowermanifold;5 tubeofconcentratingsteam; 6 smogexit;7 downwardtube;8 coolingwaterwal 工作过程中高温烟气以辐射方式与膜式水冷壁内循环水进行换热. 上升管内的水因吸热汽化, 形成汽水混合的两相流体, 它的密度要比位于同一水平面的下降管内的水密度小, 这种密度差形成了水在汽包下降管下集箱上升管上集箱导汽管汽包之间自然循环的驱动力. 燃磷塔回收余热生产的蒸汽最后在汽包内经汽水分离后排出. 13 改造效果分析燃磷量 50kg h -1 的新型燃磷塔于 001 年投入运行, 目前, 设备状况良好, 无腐蚀无泄漏无爆管等现象发生. 它不仅能保证磷酸的原有生产能力和品质, 而且可副产 0.8 MPa 的饱和蒸汽, 余热回收率达 60% 以上. 由于采用新型燃磷塔的改进工艺省掉了专用燃煤锅炉, 一方面降低了设备投资, 另一方面每燃烧 1t 混合磷可节省标准煤 560kg 和循环冷却水 400t, 相当于每生产 1t 浓度为 85% 的磷酸可降低生产成本元, 经济效益和环保效益明显. 模型描述新型燃磷塔是进行磷氧化反应和余热回收的场所, 其内部流动燃烧传热过程同时共存, 理论解析十分困难. 同时, 由于炉内的高温强腐蚀环境, 使得实验测试很难进行. 因此在分析现有工艺流程装备结构和操作参数的基础上, 对燃磷塔进 Fig.5 Simplifiedphysicalmodelof 数学模型 phosphoricacidtower 实际磷的燃烧要涉及喷枪内磷的雾化混合以及液态黄磷的蒸发, 是一个伴有蒸发燃烧的湍流两相流动过程. 计算中假定磷在喷枪中已完全汽化, 并和空气混合均匀, 即认为塔内只存在单相多组分气体. 对燃磷塔内的混合气体湍流流动采用犽 ε 双方程模型 [5], 磷的湍流反应采用涡耗散模 [6] [7] 型, 高温烟气的辐射传热采用 P1 模型. 采用上述封闭模型后, 燃磷塔内稳态湍流流动燃烧传热的控制方程组由连续性方程动量方程犽 ε 方程组分方程和能量方程等组成, 在柱坐标下可表示成以下形式 [810]. 连续性方程狓 ρ 狌 ) ρυ ) θ ρ 狑 ) =0 1) 动量方程狓 ρ 狌 φ ) ρυφ ) θ ρ 狑 φ ) = φ 狓 μ e 狓 φ μ e φ 犛 ) θ μe θ φ 犛狌 =- 狆狌狓狓 μ e υ 狓 μ e 狑狓 θ μe 狓犛 υ =- 狆狌狓 μ e υ μ e [ ] 狑狑 - θ μe - μ e 狑 θ υ ρ 犵狓 3) ρ 狑 ρ 犵 r 4)

4 1110 化工学报第 57 卷 [ ] 犛 w =- 狆 θ 狌狓 μ e υ θ μ e θ - 狑 [ ] 狑 θ μe θ υ 犽 ε 方程 μe υ θ 狑狑 - - ρ υ 狑 ρ 犵 θ 狓 ρ 狌犽 ) ρυ 犽 ) θ ρ 狑犽 ) = 狓 μ e 犽 σ 犽 θ μ e σ 犽狓 ρ 狌 ε) ρυε ) θ ρ 狑 ε) = 狓 μ e ε σε θ 犌犽 =μ t 狌狓狌 υ 狓 μ e ε σε μ e σ 犽犽狓 5) 犽犌犽 -ρε 6) θ μ e ε σε 狓 ε 犆 1 犌犽 - 犆 θ 犽 ρ ε )7) υ ) 狑 {[ θ υ ) ] 气相组分方程狑狌狓 θ υ θ 狑狑 - } 8) 狓 ρ 狌犢 S ) ρυ 犢 S ) θ ρ 狑犢 S ) = Γ 犢 S 狓狓 3 能量方程犢犻犻犻 Γ 犢 S θ Γ 犢 S θ 犚 S 9) 犚 S1 =αs 犕 S 犃 ρmin α 犻犕 10) 犚 S =αs 犕 S 犃犅 ρ ε 犢犼犽 / β 犼犕犼 11) 犼犚 S = min 犚 S1, 犚 S ) 1) 狓 ρ 狌犆犘犜 ) ρυ 犆犘犜 ) θ ρ 狑犆犘犜 ) = 狓 θ μ e 犆犘犜 σ 狓 ) μ e 犆犘犜 σ 边界条件及初始参数 μ e 犆犘犜 σ θ 狇 ra 狇 re 13) 运行实践表明, 燃磷塔水冷壁上会结成约 10 mm 的酸膜, 并在 600 附近保持动态平衡. 数值模拟中燃磷塔的壁面实际是酸膜表面, 因此壁面温度假定为 600. 流体在壁面采用无滑移的边界条件, 即轴向径向和切向速度为 0, 利用标准壁面函数进行处理. 液态磷喷入燃磷塔之前, 首先在喷磷枪中和空气混合, 其中雾化用一次空气流量为 16 m 3 h -1, 温度为 0, 采用压缩空气 ; 二次空气流量为 981m 3 h -1, 温度为 5. 计算中采用的入口条件及磷燃烧热由表 1 给出. 犜犪犫犾犲 1 犐狀犻狋犻犪犾狆狅狆犲狋犻犲狊狅犳狊狆犲犮犻犲狊犻狀犻狀犾犲狋犵犪狊 Temperature/K Velocity/m s -1 Pressure/Pa Concentration ofp4 Concentration ofo Combustionheat ofp4/j kg 网格划分及数值求解方法对燃磷塔的网格划分采用分块方法, 即把喷磷枪及与其相连的部分圆柱腔划为一个区, 烟气出口及与其相连的部分锥形封头划分为一个区, 其余被分隔开来的部分各自成区. 每个区采用不同的网格划分方法, 以提高网格的质量. 外形较不规则的喷磷枪区和烟气出口区采用四面体网格, 其他规则区域采用六面体网格, 并在喷磷枪出口区域进行网格加密. 整个计算域共划分混合网格约 50 万. 燃磷塔内的物理场采用三维稳态模拟, 隐式分离求解的方法, 即在全区域按照一定次序分别对控制方程中的每个方程独立求解. 其中动量方程的离散采用二阶迎风格式, 用 SIMPLE 算法处理压力与速度的耦合. 为保证数值精度, 其他控制方程也采用二阶迎风格式离散. 各控制方程均非线性且彼此相互耦合, 因此在获得收敛解之前要进行多次循环迭代求解. 3 物理场模拟结果及讨论 31 流场燃磷塔内烟气流动情况是影响燃烧传热质量的关键因素, 其内部流场模拟结果如图 6 图 7 所示. 图 6 是塔高犣 = m 处横截面的流线图, 图 7 是燃磷塔烟气出口及喷磷枪所在纵剖面的流线图. 由图中可以看出 : 塔内气流是一强烈旋转的气流, 这对燃料与空气的混合有利, 便于燃烧充分进行. 其中在喷磷枪轴线上有一明显的射流, 该射流在冲至相对水冷壁附近受阻, 分别向上向下流动形成两个大涡. 结合横截面上的流线图可知, 涡内流体并没有封闭不动, 而是呈渐缩螺旋形向上流动. 流体速度在左侧壁面附近较高, 靠近喷磷枪的右侧壁面速度次之, 中心区域流动速度最低. 横截面上的流线呈微不对称性, 这是由非结构网格的不对称引起的. 3 温度场温度场是新型燃磷塔数值模拟的关键问题之

5 第５期刘宝庆等膜式水冷壁结构燃磷塔的开发与传质传热的数值模拟 Ⅰ １１１１６Ｓａｍｌｎｄｄａａｍｏｄｎｃｏｃ一图８给出了燃磷塔喷磷枪入口及烟气出口所在纵截面的温度等高线由图中可以看出塔内大部分区域温度在１０００Ｋ以上其中最高温度约为１８００Ｋ出现在喷磷枪轴线上偏离燃磷塔中心的一侧烟气出口温度约为９７０Ｋ与测量值１０２０１０３８Ｋ符合较好相对误差约为６燃磷塔内下部温度较高随着高度增加烟气温度逐渐降低这是由于燃磷塔膜式水冷壁沿途不断吸热的缘故但单位高度上温度下降的幅度却越来越小这说明达到一定临界高度后高度对烟气温度分布的影响已经很小因此确定一个合理的燃磷塔高度将在保证较大余热回收效率的前提下可以达到更大的容积热强度７Ｓａｍｌｎｄｄａａｍｏｌｕｄｎａｌｃ３３组分浓度场Ｐ４Ｏ２和Ｐ４Ｏ１０的质量分数是反映燃磷塔内燃烧情况的重要参考依据图９图１１给出了喷磷枪入口和烟气出口所在纵剖面这三种组分质量份额分布的三维立体图９Ｃｃｎａｄｂｕｏｐｈｏｐｈｏｕ图９是纵剖面上燃料磷的质量份额图由图中８Ｔｍｐａｕｃｏｕｏ可以看出磷被喷入后即充分迅速燃烧在塔高ｌｕｄｎａｌｃ２ｍ以上的大部分区域其质量份额基本为零这说

6 111 化工学报第 57 卷 Fig.10 Concentrationdistributionofoxygen 支撑系统组成, 并通过下降管导汽管与汽包相连, 依靠自然循环方式吸收高温混合烟气的辐射放热. 根据新型燃磷塔内湍流流动燃烧传热过程的特点, 建立了考虑三者相互影响的综合数学模型, 并运用二阶迎风格式和 SIMPLE 算法对控制方程组进行了离散求解, 得到塔内湍流场温度场和组分场的分布. 模拟结果显示, 燃磷塔内存在明显的涡流射流结构, 磷的燃烧主要在喷磷枪雾化角范围内, 在塔高.1m 位置磷的燃烧已基本完成. 塔内最高燃烧温度为 1790 K, 出现在喷磷枪轴线上偏离燃磷塔中心的一侧. 烟气出口温度约为 970 K, 与测量值 100~1038K) 符合较好. 符号说明 Fig.11 ConcentrationdistributionofP4O10 明磷的燃烧几乎都在燃磷塔的下部靠近喷磷枪的轴线附近完成. 图 10 是中心纵剖面上氧气的质量份额分布图. 其初始质量分数约为 0%, 全区域分布趋势与黄磷相似, 只是由于空气过剩, 在燃烧完全后, 质量分数仍保持在 6% 左右. 图 11 是生成物 P 4O 10 在纵剖面的质量份额图. 由图中可以看出,P 4O 10 的质量份额在喷磷枪附近迅速突增达到一定数值后就一直保持水平不变, 这也说明燃料磷在新型燃磷塔下部即已燃烧完全, 进而证明新型燃磷塔能满足正常的热法磷酸生产工艺要求, 结构设计合理. 4 结论针对国内热法磷酸生产存在的能源浪费的不合理现象, 对现有工艺进行了改进. 改进工艺可回收余热副产蒸汽, 不仅能解决自身所需热量, 又可为其他化工生产提供蒸汽, 达到节能降耗增利的效果. 具有余热回收功效的新型燃磷塔由膜式水冷壁上下集箱上下封头喷磷枪及底部犃, 犅涡耗散模型常数犮 1, 犮, 犮 μ 犽 ε 模型常数 -1 犇 s 组分 S 的扩散系数,m s - 犽湍动能,m s 狇 re, 狇 ra 分别为反应热热辐射通量,W m -3 犚 S 化学反应速率狌犡向速度分量,m s -1 狑 θ 周向速度分量,m s -1 犢 S 组分 S 的质量分数 α, β 反应物和生成物的当量系数 -3 ε 湍动能耗散率,m s μ 分子黏性系数,Pa s -1 μ 等效黏性系数,Pa s e -1 μ 湍流黏性系数,Pa s t -1 ρ 混合物密度,kg m -3 σ 犽湍动能方程模型常数 σε 湍动能耗散率方程模型常数 υ 犚向速度分量,m s -1 犚犲犳犲犲狀犮犲狊 [1] HudsonRB,MoK.Methodsofpreparingultraphosphoric acid:us,patent [] KobayashiToshitukasa.Processunitsofphosphoricacid. JP, [3] Shute J F, Rosenhouse H. Heat recovery in the manufactureofphosphoricacid:eu, [4] LiuBaoqing 刘宝庆 ).Developmentoftwokindsofnew specialreactorsandresearchontheirmassandheattransfer models [D].Hangzhou:ZhejiangUniversity,005 [5] Launder B E,Spalding D B.Lecturesin Mathematical

7 第 5 期刘宝庆等 : 膜式水冷壁结构燃磷塔的开发与传质传热的数值模拟 Ⅰ) 1113 ModelsofTurbulence.London:AcademicPress,197 [6] MagnussenBF,HjertagerBH.Onmathematicalmodelsof turbulent combustion with special emphasis onsoot formation and combustion//proc. 16th Int. Symp. on Combustion. London:The Combustion Institute,1976: [7] SiegelR,HowelJR.ThermalRadiation HeatTransfer. Washington:Hemisphere,199 [8] John D Anderson. Computational Fluid Dynamics the Basics with Applications. Beijing: Tsinghua University Press,00 [9] TaoWenquan 陶文铨 ).NumericalHeatTransfersecond edition).xi an:xi anjiaotonguniversitypress,00 [10] ZhaGuorong 查国荣 ),ChenJiayong 陈家镛 ).Principle and Application of Transmission Process. Beijing: MetalurgicalIndustryPress, 1997 櫒櫒櫒櫒櫒櫒櫒櫒櫒櫒櫒櫒櫒櫒櫒櫒櫒櫒櫒櫒櫒櫒櫒櫒櫒櫒櫒櫒櫒櫒櫒櫒櫒櫒櫒櫒櫒櫒櫒櫒櫒櫒櫒櫒櫒櫒櫒櫒櫒檭檭檭檭檭檭檭檭檭檭檭檭檭檭檭檭檭檭信息与交流第三届全国化学工程与生物化工年会 3rdChineseNationalChemicalandBiochemicalEngineeringAnnualMeeting 第一轮通知第一届全国化学工程与生物化工年会于 004 年 11 月 5~9 日在南京成功召开,18 位两院院士到会,738 位代表参会, 盛况空前 ; 第二届全国化学工程与生物化工年会于 005 年 11 月 4~8 日在北京召开, 参会代表千余名, 包括中国科学院和中国工程院院士 17 名荷兰皇家科学院院士和澳大利亚科学院院士各 1 名会议常设执委会确定第三届全国化学工程与生物化工年会于 006 年 11 月 3~7 日在南宁召开, 会议由广西大学广西科学院共同承办会议主题为化工 / 生物化工科学技术与生物质能源本次会议将有多位两院院士和国际著名学者出席并作大会报告, 会议期间还拟举行重大前沿领域的专题研讨会, 部分跨国公司国内大公司人力资源部门负责人与参会研究生见面会和相关业务推介活动热诚欢迎各高校科研院所 ) 以及企事业单位的科研工程技术和管理人员积极参加会议现将会议有关事项通知如下 : 一会议组织机构 : 1. 执委会常设 ) 主席 : 欧阳平凯 ; 副主席 : 汪文川, 徐南平, 谭天伟委员 : 王静康, 白骅, 刘汉兴, 李伯耿, 李静海, 芮新生, 沈寅初, 张嗣良, 何志敏, 欧阳藩, 袁渭康, 曹竹安, 戴猷元. 名誉学术委员会曹湘洪, 陈丙珍, 陈冠荣, 陈家镛, 陈俊武, 费维扬, 高从, 高金吉, 关兴亚, 郭慕孙, 何鸣元, 侯芙生, 胡英, 胡永康, 金涌, 李灿, 李大东, 李洪钟, 李静海, 李俊贤, 李正名, 林励吾, 陆婉珍, 闵恩泽, 欧阳平凯, 桑凤亭, 沈寅初, 时铭显, 舒兴田, 汪家鼎, 汪燮卿, 王静康, 王文兴, 魏可镁, 谢克昌, 徐南平, 杨锦宗, 杨胜利, 余国琮, 袁渭康, 袁晴棠, 袁权, 张懿 3. 执委会秘书长 : 宋宗文, 屈一新下转 113 页 )

犐犆犛犆 57 犌犅犣 GBZ 犛狆犲犮犻犳犻犮犪狋犻狅狀狊犳狅狉犻狀犱犻狏犻犱狌犪犾犿狅狀犻狋狅狉犻狀犵狅犳狅犮犮狌狆犪狋犻狅狀犪犾犲狓狋犲狉狀犪犾犲狓狆狅狊狌狉犲 !"

犐犆犛犆 57 犌犅犣 GBZ 犛狆犲犮犻犳犻犮犪狋犻狅狀狊犳狅狉犻狀犱犻狏犻犱狌犪犾犿狅狀犻狋狅狉犻狀犵狅犳狅犮犮狌狆犪狋犻狅狀犪犾犲狓狋犲狉狀犪犾犲狓狆狅狊狌狉犲 ! 犐犆犛 13.100 犆 57 GBZ128 2002 犛狆犲犮犻犳犻犮犪狋犻狅狀狊犳狅狉犻狀犱犻狏犻犱狌犪犾犿狅狀犻狋狅狉犻狀犵狅犳狅犮犮狌狆犪狋犻狅狀犪犾犲狓狋犲狉狀犪犾犲狓狆狅狊狌狉犲 2016 06 28 2016 11 01!" Ⅲ 1 1 2 1 3 1 4! 2 5 "#$%&' 3 6 $% 4 7 %() 5 8 *+,-./ 6 0+A(1 0+) "#*+.2 8 Ⅰ 4 5,!