化学回热循环与水系统仿真研究
|
摘要
化学回热循环以其优越的热力学性能(包括输出比功率大、热效率高和NO X排放量小等),成为解决现阶段能源危机和环境污染问题的重要方法之一。基于我国先进循环装置(STIG循环、HAT循环和联合循环等)的研究思路,本文以某型船用三轴简单循环燃机为母型机,采用面向对象模块化建模方法,对CRGT循环进行对比仿真研究。本文主要包括以下工作内容:
1)采用VC语言编制了工质热力性质(包括水和水蒸气、干空气和湿燃气等)的热力计算程序;采用MEX接口建立了VC软件和Matlab软件之间函数调用的动态链接库;
2)采用Matlab语言编制了先进循环部件(包括燃机部件、化学回热器和蒸汽发生器部件)的热力计算程序,编制了错流式翅片管换热器的结构设计程序;采用BP神经网络法拟合了部件特性曲线,编制了部件特性的查询程序;
3)采用S函数建立了包括所有程序的Simulink模块库。基于该模块库,分别进行了先进循环的性能优化和动态仿真研究。通过仿真计算,得出先进循环(包括CRGT循环和STIG循环)中水系统(蒸汽发生器)的性能优化规律、变工况性能网络图(即输出功率、热效率与喷油量、蒸汽回注量关系)和变工况控制策略(即快速、安全地变工况运行的供油-供汽规律)。
With the excellent thermodynamic properties,including the large specific power output,the high thermal efficiency and the low emissions of NOx and so on,the CRGT cycle (Chemically recuperated gas turbine) has become one of the most important methods to solve the energy crisis and environmental problems at this stage.Base on our research ideas of the advanced cycle gas turbine plants (including STIG cycle, HAT cycle and combined cycle,etc),this thesis use the modular object-oriented modeling method to study the CRGT performance contrastively ,with the basis of one certain marine-turbine type.In this thesis,the following elements are completed:
1)Use the VC Language to compile the thermodynamic calculation programs of the cycle medium’s thermodynamic properties(including water and water vapor, dry air and wet gas);use the MEX interface to establish the dynamic link library between the VC software and Matlab software.
2)Use the Matlab Language to compile the thermodynamic calculation programs of the advanced cycle components(including the gas turbine parts,chemi -cal recuperator and steam generator,etc)and the structural design process of the cross-flow type finned tube heat exchanger;use the BP neural network method to fit the components characteristic curves and compile its inquiry programs.
3)Use the S function to build the Simulink block library of all programs. Based on this module library,carry out the advanced cycle’s performance optimiza -tion and the dynamic simulation study.
Through the simulation,the optimization rules of the advanced cycle water system (including CRGT cycle and STIG cycle),the variable load performance network diagrams (including the correlation of the output power and thermal efficiency ,with the quality of fuel injection and steam injection),and the control strategy of variable operating conditions (namely the fuel and steam supply rate’s operate strategy,which made the turbine system operating quickly and safely between the off-design conditions) have been reached.
1)采用VC语言编制了工质热力性质(包括水和水蒸气、干空气和湿燃气等)的热力计算程序;采用MEX接口建立了VC软件和Matlab软件之间函数调用的动态链接库;
2)采用Matlab语言编制了先进循环部件(包括燃机部件、化学回热器和蒸汽发生器部件)的热力计算程序,编制了错流式翅片管换热器的结构设计程序;采用BP神经网络法拟合了部件特性曲线,编制了部件特性的查询程序;
3)采用S函数建立了包括所有程序的Simulink模块库。基于该模块库,分别进行了先进循环的性能优化和动态仿真研究。通过仿真计算,得出先进循环(包括CRGT循环和STIG循环)中水系统(蒸汽发生器)的性能优化规律、变工况性能网络图(即输出功率、热效率与喷油量、蒸汽回注量关系)和变工况控制策略(即快速、安全地变工况运行的供油-供汽规律)。
With the excellent thermodynamic properties,including the large specific power output,the high thermal efficiency and the low emissions of NOx and so on,the CRGT cycle (Chemically recuperated gas turbine) has become one of the most important methods to solve the energy crisis and environmental problems at this stage.Base on our research ideas of the advanced cycle gas turbine plants (including STIG cycle, HAT cycle and combined cycle,etc),this thesis use the modular object-oriented modeling method to study the CRGT performance contrastively ,with the basis of one certain marine-turbine type.In this thesis,the following elements are completed:
1)Use the VC Language to compile the thermodynamic calculation programs of the cycle medium’s thermodynamic properties(including water and water vapor, dry air and wet gas);use the MEX interface to establish the dynamic link library between the VC software and Matlab software.
2)Use the Matlab Language to compile the thermodynamic calculation programs of the advanced cycle components(including the gas turbine parts,chemi -cal recuperator and steam generator,etc)and the structural design process of the cross-flow type finned tube heat exchanger;use the BP neural network method to fit the components characteristic curves and compile its inquiry programs.
3)Use the S function to build the Simulink block library of all programs. Based on this module library,carry out the advanced cycle’s performance optimiza -tion and the dynamic simulation study.
Through the simulation,the optimization rules of the advanced cycle water system (including CRGT cycle and STIG cycle),the variable load performance network diagrams (including the correlation of the output power and thermal efficiency ,with the quality of fuel injection and steam injection),and the control strategy of variable operating conditions (namely the fuel and steam supply rate’s operate strategy,which made the turbine system operating quickly and safely between the off-design conditions) have been reached.
引文
[1] K.F.Kesser,M.A.Hoffman.Analysis of a Basic Chemically Eecuper -ated Gas Turbine Power Plant. Journal of Engineering for Gas Turbines and Power,1994,116:277-284.
[2] S.T.Adelman,M.A.Hoffman.A Methane-Steam Reformer for a Basic Chemically Recuperated Gas Turbine.Transactions of the ASME, 1995,17:16-23
[3] H.Abdallah,S.Harvey.Thermodynamic analysis of chemically recu -perated gas turbines.Int.J.Therm.Sci,2001,40:372-384
[4] Y.H.Chin,R.Dagle.Steam reforming of methanol over highly active Pd/ZnO catalst.Catalysis Today,2002,77:79-88
[5] T.Nakagaki,T.Ogawa.Development of Chemically Recuperated Micro Gas Turbine.Journal of Engineering for Gas Turbines and Power 2003,125:391-397
[6] G.Kolb,R.Zapf.Propane steam reforming in micro-channels-resul ts from catalyst screening and optimisation.Applied Catalysis A:General,2004,277:155-166
[7] R.Carapellucci,A.Milazzo.Thermodynamic optimization of a reh -eat chemically recuperated gas turbine.Energy Conversion and Management,2005,46:2936-2953
[8] D.Cocco,V.Tola.Performance ecaluation of chemically recuperat -ed gas turbine(CRGT) power plants fuelled by di-methyl-ether (DME).Energy,2006,31:1446-1458
[9]肖云汉,蔡睿贤,林汝谋.化学回热循环的热力学优化.工程热物理学报,1996,17(增刊):1-4
[10]李志敏,张密林,刘顺隆.CRGT循环的热力学性能分析.化学工程师,2001,82(1):28-30
[11]李政,王德慧.以甲醇为燃料的燃气轮机技术及研究进展综述.热力透平,2004,33(3):179-183
[12]徐威.应用天然气重整技术的新型动力系统研究.中国科学院研究生院硕士学位论文.2005:70-72
[13]曹文,郑丹星.天然气CO2转化化学回热热力循环.东南大学学报,2006,36(5):790-794
[14]苟晨华,蔡睿贤,洪慧.新型甲醇燃料化学回热燃气轮机循环.中国电机工程学报,2006,26(7):1:5
[15]张筱松,金红光等.基于甲醇化学链燃烧的新型燃气轮机间冷循环.工程热物理学报,2007,28(5):721-724
[16]张志.甲烷裂解高效催化剂制备.中国工程物理研究院硕士学位论文.2007:54-56
[17] E.N.Houstis,A.C.Catlin.GasTurbnLab:a multidisciplinary proble -m solving environment for gas turbine engine design on a network of nonhomogeneous machines.Journal of Computational and Appli -ed Mathematics,2002,149:83-100
[18] A.E.Ruano,P.J.Fleming.Nonlinear identification of aircraft gas -turbine dynamics.Neurocomputing,2003,55:551-579
[19] A.M.Bassily.Performance improvements of the intercooled reheat recuperated gas-turbine cycle using absorption inlet-cooling and evaporative after-cooling.Applied Energy,2004,77:249- 272
[20] M.Jonsson,J.Yan.Humidified gas turbines-a review of proposed and implemented cycles.Energy,2005,30:1013-1078
[21] T.S.Kim,S.H.Hwang.Part load performance analysis of recuperated gas turbines considering engine configuration and operationstrategy.Energy,2006,31:260-277
[22] J.A.Gonzalez,J.M.Sala.Modelling and dynamic simulation of processes with‘MATLAB’.An application of a natural gas installation in a power plant.Energy,2007:1271-1282
[23] A.L.Polyzakis,C.Koroneos,G.Xydis.Optimum gas turbine cycle for combined cycle power plant.Energy Conversion and anagement, 2008,49:551-563
[24]谢志武,陈德来,翁史烈.面向对象的燃气轮机仿真建模:综述与展望.热能动力工程,1998,13(76):243-246
[25]敖晨阳,张宁,陈华清.基于MATLAB的三轴燃气轮机动态仿真模型研究.热能动力工程,2001,16(95):523-526
[26]高阳.基于联合循环的燃气轮机动态模型与仿真.重庆大学硕士学位论文.2002:44-51
[27]赵景峰,叶春,秦春申.燃气轮机及其控制系统的综合建模研究.华东电力,2005,33(4):13-16
[28]刘喜超.单轴燃气轮机动态仿真.重庆大学硕士学位论文.2006:31-45
[29]周文祥.航空发动机及控制系统建模与面向对象的仿真研究.南京航空航天大学博士学位论文.2006:81-82
[30]肖琳.气垫船用三轴燃气轮机推进系统动态性能研究.上海交通大学硕士学位论文.2007:50-63
[31]王敏,李扬,李淑英.基于Matlab/Simulink的单轴燃气轮机动态仿真.燃气轮机技术,2006,19(4):47-50
[32]张泰岩.基于微型燃气轮机的冷热电联产系统仿真.华北电力大学硕士学位论文.2006:2-10,43
[33]贾省伟.舰船双轴燃气轮机性能仿真.华中科技大学硕士学位论文.2006:20-53
[34]李海英.生物质能利用BI/STIG系统性能的研究.华北电力大学硕士学位论文.2006:15-50
[35]张士杰,肖云汉.基于微燃机的HAT循环变工况性能分析.工程热物理学报,2006,27(4):545-548
[36]余又红,王义等.舰用中冷回热燃气轮机装置动态性能分析.海军工程大学学报,2007,19(1):81-85
[37]李一兴.湿空气透平(HAT)循环饱和器若干问题的研究.上海交通大学博士学位论文.2007:121-124
[38]刘春琪.HAT循环饱和器性能研究.上海交通大学硕士学位论文.2008:73-80
[39]孙士恩.焦炉煤气与气化煤气基焦炭、化工与动力多联产系统及集成机理.中国科学院研究生院博士学位论文.2009:119-121
[40]杨锦成.微型燃气轮机冷热电联供系统集成与性能仿真研究.上海交通大学博士学位论文.2009:31-63
[41]姚寿广.船舶热力系统分析.北京:科学出版社,2003:12-55
[42]刘志刚.工质热物理性质计算程序的编制及应用.北京:科学出版社1992:131-136
[43] W.瓦格纳,A.克鲁泽[德].水和蒸汽的性质.北京:高等教育出版社2003:5-26
[44]杨文滨,苏明.工质加湿对燃气透平性能的影响分析.中国动力工程学报,2005,25(2):185-187
[45]陈安斌,尚德敏.适用于燃煤气的STIG循环中湿燃气的状态方程.热能动力工程,2000,15(87):290-293
[46]陈安斌.湿燃气热力性质的计算方法.哈尔滨工业大学学报,2001, 33(6):799-801
[47]杨文滨,苏明.高湿空气在高温高压区域内的热物性计算及分析.动力工程,2004,24(3):383-387
[48]李传军.C语言与MATLAB接口编程与实例.北京:北京邮电大学出版社,2004:8-26
[49]严家禄,尚德敏.湿空气和烃燃气热力性质图表.北京:高等教育出版社,1989:20-27
[50]翁史烈.燃气轮机性能分析.上海:上海交通大学出版社,1985:191-196
[51]刘光宇.船舶燃气轮机装置原理与设计.哈尔滨:哈尔滨船舶工程学院出版社,1992:27-67
[52]倪维斗,徐向东,李政等.热动力系统建模与控制的若干问题.北京:科技出版社,1996:33-70
[53]钱颂文主编.换热器设计手册.北京:化学工业出版社,2002:305-325
[54]张娜,林汝谋,蔡睿贤.压气机特性通用数学表达式.工程热物理学报,1996,17(1):21-24
[55]胡松岩,黄忠湖,王月奇.燃气涡轮通用理论特性计算方法.航空发动机,1994(3):34-47
[56]张娜,蔡睿贤,单轴恒速燃气轮机及其功热并供装置的变工况显式解析解.工程热物理学报,1998,19(2):141-144
[57]李涛,王德慧,李政.船用分轴燃气轮机主要部件变工况性能的初步研究.汽轮机技术,2004,46(5):365-368
[58]彭靖波,谢寿生.基于RBF神经网络的压气机特性仿真.推进技术2006,27(1):30-32
[59]冯凤珍,孙允标,余又红.基于神经网络的压气机特性计算新方法.计算机与数字工程,2004,32(6):45-47
[60] K.Ghorbanian,M.Gholamrezaei.An artificial neural networkappro -ach to compressor performance prediction.Applied Energy, 2009,86:1210-1221
[61] Youhong Yu,Lingen Chen.Neural-network based analysis and predic -tion of a compressor’s characteristic performance map. Applied Energy,2007,84:48-55
[62] O.Cortes,G.Urquiza,J.A.Hernandez.Optimzation of operating co -nditions for compressor performance by means of neural network inverse.Applied Energy,2009,86:2487-2493
[63]刘喜超,唐胜利.基于偏最小二乘法的压气机特性曲线的拟合.汽轮机技术,2006,48(5):327-329
[64]杨欣毅,沈伟,刘海峰.一种应用滑动最小二乘求取压气机特性的方法.航空动力学报,2009,24(8):1741-1746
[65] R.Tirnovan,S.Giurgea,A.Miraoui.Modelling the characteristics of fuel systems using hybrid method based on moving least squa -res.Applied Energy,2009,86:1283-1289
[66] R.Tirnovan,S.Giurgea,A.Miraoui,M.Cirrincione.Surrogate modell -ing of compressor characteristics for fuel-cell applications. Applied Energy,2008,85:394-403
[67]周开利,康耀红.神经网络模型及其MATLAB仿真程序设计.北京:清华大学出版社,2005:70-101
[68]张德丰.MATLAB程序设计与典型应用.北京:电子工业出版社2009:207-222
[69]张志涌编著.精通MATLAB6.5版.北京:北京航空航天大学出版社,2007:360-380
[70]黄永安,马路.MATLAB7.0/Simulink6.0建模仿真开发与高级工程应用.北京:清华大学出版社,2006:76-183
[2] S.T.Adelman,M.A.Hoffman.A Methane-Steam Reformer for a Basic Chemically Recuperated Gas Turbine.Transactions of the ASME, 1995,17:16-23
[3] H.Abdallah,S.Harvey.Thermodynamic analysis of chemically recu -perated gas turbines.Int.J.Therm.Sci,2001,40:372-384
[4] Y.H.Chin,R.Dagle.Steam reforming of methanol over highly active Pd/ZnO catalst.Catalysis Today,2002,77:79-88
[5] T.Nakagaki,T.Ogawa.Development of Chemically Recuperated Micro Gas Turbine.Journal of Engineering for Gas Turbines and Power 2003,125:391-397
[6] G.Kolb,R.Zapf.Propane steam reforming in micro-channels-resul ts from catalyst screening and optimisation.Applied Catalysis A:General,2004,277:155-166
[7] R.Carapellucci,A.Milazzo.Thermodynamic optimization of a reh -eat chemically recuperated gas turbine.Energy Conversion and Management,2005,46:2936-2953
[8] D.Cocco,V.Tola.Performance ecaluation of chemically recuperat -ed gas turbine(CRGT) power plants fuelled by di-methyl-ether (DME).Energy,2006,31:1446-1458
[9]肖云汉,蔡睿贤,林汝谋.化学回热循环的热力学优化.工程热物理学报,1996,17(增刊):1-4
[10]李志敏,张密林,刘顺隆.CRGT循环的热力学性能分析.化学工程师,2001,82(1):28-30
[11]李政,王德慧.以甲醇为燃料的燃气轮机技术及研究进展综述.热力透平,2004,33(3):179-183
[12]徐威.应用天然气重整技术的新型动力系统研究.中国科学院研究生院硕士学位论文.2005:70-72
[13]曹文,郑丹星.天然气CO2转化化学回热热力循环.东南大学学报,2006,36(5):790-794
[14]苟晨华,蔡睿贤,洪慧.新型甲醇燃料化学回热燃气轮机循环.中国电机工程学报,2006,26(7):1:5
[15]张筱松,金红光等.基于甲醇化学链燃烧的新型燃气轮机间冷循环.工程热物理学报,2007,28(5):721-724
[16]张志.甲烷裂解高效催化剂制备.中国工程物理研究院硕士学位论文.2007:54-56
[17] E.N.Houstis,A.C.Catlin.GasTurbnLab:a multidisciplinary proble -m solving environment for gas turbine engine design on a network of nonhomogeneous machines.Journal of Computational and Appli -ed Mathematics,2002,149:83-100
[18] A.E.Ruano,P.J.Fleming.Nonlinear identification of aircraft gas -turbine dynamics.Neurocomputing,2003,55:551-579
[19] A.M.Bassily.Performance improvements of the intercooled reheat recuperated gas-turbine cycle using absorption inlet-cooling and evaporative after-cooling.Applied Energy,2004,77:249- 272
[20] M.Jonsson,J.Yan.Humidified gas turbines-a review of proposed and implemented cycles.Energy,2005,30:1013-1078
[21] T.S.Kim,S.H.Hwang.Part load performance analysis of recuperated gas turbines considering engine configuration and operationstrategy.Energy,2006,31:260-277
[22] J.A.Gonzalez,J.M.Sala.Modelling and dynamic simulation of processes with‘MATLAB’.An application of a natural gas installation in a power plant.Energy,2007:1271-1282
[23] A.L.Polyzakis,C.Koroneos,G.Xydis.Optimum gas turbine cycle for combined cycle power plant.Energy Conversion and anagement, 2008,49:551-563
[24]谢志武,陈德来,翁史烈.面向对象的燃气轮机仿真建模:综述与展望.热能动力工程,1998,13(76):243-246
[25]敖晨阳,张宁,陈华清.基于MATLAB的三轴燃气轮机动态仿真模型研究.热能动力工程,2001,16(95):523-526
[26]高阳.基于联合循环的燃气轮机动态模型与仿真.重庆大学硕士学位论文.2002:44-51
[27]赵景峰,叶春,秦春申.燃气轮机及其控制系统的综合建模研究.华东电力,2005,33(4):13-16
[28]刘喜超.单轴燃气轮机动态仿真.重庆大学硕士学位论文.2006:31-45
[29]周文祥.航空发动机及控制系统建模与面向对象的仿真研究.南京航空航天大学博士学位论文.2006:81-82
[30]肖琳.气垫船用三轴燃气轮机推进系统动态性能研究.上海交通大学硕士学位论文.2007:50-63
[31]王敏,李扬,李淑英.基于Matlab/Simulink的单轴燃气轮机动态仿真.燃气轮机技术,2006,19(4):47-50
[32]张泰岩.基于微型燃气轮机的冷热电联产系统仿真.华北电力大学硕士学位论文.2006:2-10,43
[33]贾省伟.舰船双轴燃气轮机性能仿真.华中科技大学硕士学位论文.2006:20-53
[34]李海英.生物质能利用BI/STIG系统性能的研究.华北电力大学硕士学位论文.2006:15-50
[35]张士杰,肖云汉.基于微燃机的HAT循环变工况性能分析.工程热物理学报,2006,27(4):545-548
[36]余又红,王义等.舰用中冷回热燃气轮机装置动态性能分析.海军工程大学学报,2007,19(1):81-85
[37]李一兴.湿空气透平(HAT)循环饱和器若干问题的研究.上海交通大学博士学位论文.2007:121-124
[38]刘春琪.HAT循环饱和器性能研究.上海交通大学硕士学位论文.2008:73-80
[39]孙士恩.焦炉煤气与气化煤气基焦炭、化工与动力多联产系统及集成机理.中国科学院研究生院博士学位论文.2009:119-121
[40]杨锦成.微型燃气轮机冷热电联供系统集成与性能仿真研究.上海交通大学博士学位论文.2009:31-63
[41]姚寿广.船舶热力系统分析.北京:科学出版社,2003:12-55
[42]刘志刚.工质热物理性质计算程序的编制及应用.北京:科学出版社1992:131-136
[43] W.瓦格纳,A.克鲁泽[德].水和蒸汽的性质.北京:高等教育出版社2003:5-26
[44]杨文滨,苏明.工质加湿对燃气透平性能的影响分析.中国动力工程学报,2005,25(2):185-187
[45]陈安斌,尚德敏.适用于燃煤气的STIG循环中湿燃气的状态方程.热能动力工程,2000,15(87):290-293
[46]陈安斌.湿燃气热力性质的计算方法.哈尔滨工业大学学报,2001, 33(6):799-801
[47]杨文滨,苏明.高湿空气在高温高压区域内的热物性计算及分析.动力工程,2004,24(3):383-387
[48]李传军.C语言与MATLAB接口编程与实例.北京:北京邮电大学出版社,2004:8-26
[49]严家禄,尚德敏.湿空气和烃燃气热力性质图表.北京:高等教育出版社,1989:20-27
[50]翁史烈.燃气轮机性能分析.上海:上海交通大学出版社,1985:191-196
[51]刘光宇.船舶燃气轮机装置原理与设计.哈尔滨:哈尔滨船舶工程学院出版社,1992:27-67
[52]倪维斗,徐向东,李政等.热动力系统建模与控制的若干问题.北京:科技出版社,1996:33-70
[53]钱颂文主编.换热器设计手册.北京:化学工业出版社,2002:305-325
[54]张娜,林汝谋,蔡睿贤.压气机特性通用数学表达式.工程热物理学报,1996,17(1):21-24
[55]胡松岩,黄忠湖,王月奇.燃气涡轮通用理论特性计算方法.航空发动机,1994(3):34-47
[56]张娜,蔡睿贤,单轴恒速燃气轮机及其功热并供装置的变工况显式解析解.工程热物理学报,1998,19(2):141-144
[57]李涛,王德慧,李政.船用分轴燃气轮机主要部件变工况性能的初步研究.汽轮机技术,2004,46(5):365-368
[58]彭靖波,谢寿生.基于RBF神经网络的压气机特性仿真.推进技术2006,27(1):30-32
[59]冯凤珍,孙允标,余又红.基于神经网络的压气机特性计算新方法.计算机与数字工程,2004,32(6):45-47
[60] K.Ghorbanian,M.Gholamrezaei.An artificial neural networkappro -ach to compressor performance prediction.Applied Energy, 2009,86:1210-1221
[61] Youhong Yu,Lingen Chen.Neural-network based analysis and predic -tion of a compressor’s characteristic performance map. Applied Energy,2007,84:48-55
[62] O.Cortes,G.Urquiza,J.A.Hernandez.Optimzation of operating co -nditions for compressor performance by means of neural network inverse.Applied Energy,2009,86:2487-2493
[63]刘喜超,唐胜利.基于偏最小二乘法的压气机特性曲线的拟合.汽轮机技术,2006,48(5):327-329
[64]杨欣毅,沈伟,刘海峰.一种应用滑动最小二乘求取压气机特性的方法.航空动力学报,2009,24(8):1741-1746
[65] R.Tirnovan,S.Giurgea,A.Miraoui.Modelling the characteristics of fuel systems using hybrid method based on moving least squa -res.Applied Energy,2009,86:1283-1289
[66] R.Tirnovan,S.Giurgea,A.Miraoui,M.Cirrincione.Surrogate modell -ing of compressor characteristics for fuel-cell applications. Applied Energy,2008,85:394-403
[67]周开利,康耀红.神经网络模型及其MATLAB仿真程序设计.北京:清华大学出版社,2005:70-101
[68]张德丰.MATLAB程序设计与典型应用.北京:电子工业出版社2009:207-222
[69]张志涌编著.精通MATLAB6.5版.北京:北京航空航天大学出版社,2007:360-380
[70]黄永安,马路.MATLAB7.0/Simulink6.0建模仿真开发与高级工程应用.北京:清华大学出版社,2006:76-183
© 2004-2018 中国地质图书馆版权所有 京ICP备05064691号 京公网安备11010802017129号 地址：北京市海淀区学院路29号 邮编：100083 电话：办公室：(+86 10)66554848；文献借阅、咨询服务、科技查新：66554700 |