1、辽宁工业大学电力系统自动化课程设计（论文）题目机组自动发电控制系统设计（1）院（系）电气工程学院专业班级电气113学号110303089学生姓名指导教师起止时间201412120131214本科生课程设计（论文）课程设计（论文）任务及评语院（系）电气工程学院教研室电气工程及其自动化注成绩平时20论文质量60答辩20以百分制计算学号学生姓名专业班级课程设计题目机组自动发电控制系统设计（1）课程设计（论文）任务基本参数1发电厂有三机组，一号机组功率为50MW，二号机组功率为100MW，三号机组功率为70MW，功率因数均为09。2有功功率调差系数为006。3负荷频率调节效应系数（有功功率）为15。4
2、各发电机均以80负荷运行。设计要求1阐述发电机有功频率调节的基本原理。2系统负荷一定和变化时，分析机组有功功率调差系数和负荷频率调节效应系数对功率分配的影响，推导有功功率分配的公式。3确定合适的方案，使负荷变化时，机组的有功功率按各自的额定功率合理分配。4采用单片机，设计某台机组的自动发电控制系统，实现机组功率的合理分配。5对设计进行分析总结。进度计划1、布置任务，查阅资料，理解掌握系统设计要求。（1天）2、阐述发电机有功频率调节和自动发电的基本原理和实现方法。（1天）3、推导有功功率分配的公式，确定机组有功功率分配方案（2天）4、自动发电控制系统总体方案设计。（1天）5、系统硬件电路设计。（
3、2天）6、系统软件设计（2天）7、撰写、打印设计说明书（1天）指导教师评语及成绩平时论文质量答辩总成绩指导教师签字年月日本科生课程设计（论文）摘要自动发电控制（AGC）就是通过监视电厂出力和负荷系统之间的差异，来控制调频机组的出力，以满足不断变化的用户电力需要，达到电能的发供平衡，并且使整个系统处于经济的运行状态。本次课程设计对三台机组并联运行时各自有功率的分配进行设计当某台发电机组负荷改变导致频率改变时，由SCADA发送改变信号通过A/D转换器给单片机89C51，再由单片机89C51通过程序识别信号，再由输出口通过D/A转换器向AGC发送执行指令。最后AGC向发电机组发送改变运行状态指令。关
4、键词AGC；89C51；SCADA；经济运行状态本科生课程设计（论文）目录第1章绪论111自动发电控制概述112本文主要内容1第2章机组并联运行有功功率分配计算321机组有功功率频率控制及自动发电的基本原理322单台机组有功控制的基本方法623负荷变化时的功率分配计算6第3章自动发电系统硬件设计931自动发电系统功能932自动发电总体设计方案933单片机最小系统设计1033189C51单片机引脚功能10332复位电路设计11333时钟电路设计12334直流稳压电源设计12335单片机最小系统1334自动发电控制原理图设计14第4章自动发电系统软件设计1541软件功能1542软件流程图15第5章
5、课程设计总结16参考文献17本科生课程设计（论文）1第1章绪论11自动发电控制概述自动发电控制是并网发电厂提供的有偿辅助服务之一，发电机组在规定的出力调整范围内，跟踪电力调度交易机构下发的指令，按照一定调节速率实时调整发电出力，以满足电力系统频率和联络线功率控制要求的服务。或者说,自动发电控制AGC对电网部分机组出力进行二次调整,以满足控制目标要求其基本功能为负荷频率控制（LFC），经济调度控制（EDC），备用容量监视（RM），AGC性能监视（AGCPM），联络线偏差控制（TBC）等；以达到其基本的目标保证发电出力与负荷平衡，保证系统频率为额定值，使净区域联络线潮流与计划相等，最小区域化运行成
6、本。自动发电控制着重解决电力系统在运行中的频率调节和负荷分配问题，以及与相邻电力系统间按计划进行功率交换。电力系统的供电频率是系统正常运行的主要参数之一。系统电源的总输出功率与包括电力负荷在内的功率消耗相平衡时，供电频率保持恒定；若总输出功率与总功率消耗之间失去平衡时,频率就发生波动,严重时会出现频率崩溃。电力系统的负荷是不断变化的，这种变化有时会引起系统功率不平衡，导致频率波动。要保证电能的质量，就必须对电力系统频率进行监视和调整。当频率偏离额定值后，调节发电机的出力以使电力系统的有功功率达到新的平衡，从而使频率能维持在允许范围之内。所以，自动发电控制是通过对供电频率的监测、调整实现的。一个
7、大电力系统是由几个区域电力系统通过联络线互联构成。各区域电力系统按预定计划进行功率交换。每一个区域电力系统的负荷、线路损耗与联络线净交换功率之和必须与该地区的发电出力相等。12本文主要内容电力系统自动发电控制的实施是一项复杂的系统工程，在一个大型电力系统中，自动发电控制的实施涉及到几十个、乃至上百个单位的协作；在专业技术的跨度上，则涉及调度、运行方式、通信、自动化、热工、汽机、锅炉、水工等专业；即使是在技术条件基本具备的情况下，自动发电控制的实施从起步到比较完本科生课程设计（论文）2善，也要经过几年时间的努力才能完成。因此，一个大型电力系统（特别是互联电力系统）自动发电控制的实施，应该制定实施
8、规划，明确总体目标和实施步骤，研究技术原则和技术经济政策，对发电资源的利用进行优化组合，对发电机组、信息系统等的技术改造作出安排。电力系统自动发电控制实施规划是开展自动发电控制工作的纲领性文件，通过制定规划和实施规划，有助于在互联电力系统中对这项工作的开展统一思想、统一政策、统一步调，经济、高效地实现预定目标。华东电力系统的成功实践充分证明了自动发电控制实施规划的作用。自动发电控制着重解决电力系统在运行中的频率调节和负荷分配问题，以及与相邻电力系统间按计划进行功率交换。因此本次课程设计主要针对发电厂机组自动发电控制系统进行设计。首先对机组的有功功率频率特性进行系统研究，阐述自动发电控制的基本原
9、理、功能与实现方法。最后，采用单片机，设计某台机组的自动发电控制系统，实现机组功率的合理分配。本科生课程设计（论文）3第2章机组并联运行有功功率分配计算21机组有功功率频率控制及自动发电的基本原理频率是衡量电能质量的重要指标，频率质量的下降不仅影响用户的用电质量，同时对电力系统本身影响也很大，严重时可造成系统崩溃。为保证用户和电厂的正常运行和安全，我们需要进行频率调整与控制，使得系统的频率波动不超过02HZ。根据符合变动的不同特点，可将调整划归为一次调节和二次调节及互联系统的频率调节。在这里着重分析前两种调整的原理和方法。电力系统频率的一次调节是指利用系统固有的负荷频率特性，以及发电机的调速器
10、作用，来阻止系统频率偏离标准的调节方式。（1）其原理如下当电力系统中原动机功率或负荷功率发生变化时，必然引起电力系统频率的变化，此时，存储在系统负荷的电磁场和旋转质量（如电动机、照明镇流器等）中的能量会发生变化，以阻止系统频率的变化，即当系统频率下降时，系统负荷减少；当系统频率上升时，系统负荷会增加。这称为系统负荷的惯性作用，它用负荷的频率调节效应系数（又称为系统负荷阻尼常数）D来表示/HZMWFPD（21）系统负荷阻尼常数D常用标幺值来表示，其典型值为12。D2意味着1的频率变化会引起系统负荷2的变化。（2）其作用当电力系统频率发生变化时，系统中所有的发电机转速即发生变化，如转速的变化超出发
11、电机组规定的不灵敏去，该发电机的调速器就会动作，改变其原动机的阀门位置，调整原动机的功率，以求改善原动机功率或负荷功率的不平衡状况，即当系统频率下降时，发电机的蒸汽阀门或进水阀门的开度就会增大，增加原动机的功率；当系统频率上升时，发电机的蒸汽阀门或进水阀门的开度就减少，减少原动机的功率。发电机调速器的这种特性称为机组的调差特性，它用调差率R来表示/RONNNR（22）式中ON表示无载静态转速（主阀在无载位置）N表示满载静态转速（主阀全开）RN表示额定转速本科生课程设计（论文）4调差率R的实际含义是，如R5，则系统频率变化5，将引起主阀位置变化100。（3）特点第一，响应速度快，一般在在10S；
12、第二，作用时间短暂，一般在05到2分钟不等；第三，一次调节是有差调节，所有机组的调整只与一个变量有关，机组之间互相影响小，但无法实现无差调整。电力系统的二次调节由于发电机组一次调节实行的是频率有差调节，因此，早期的频率二次调节，是通过控制调速系统的同步电机，改变发电机组的调差特性曲线的位置，实现频率的无差调整。但未实现对火力发电机组的燃烧系统的控制，为使原动机的功率与负荷功率保持平衡，需要依靠人工调整原动机功率的基准值，达到改变原动机功率的目的。随着科学技术的进步，火力发电机组普遍采用了协调的控制系统，由自动控制来代替人工进行此类操作。在现代化电力系统中，各控制区则采用集中的计算机控制。这就是
13、电力系统频率的二次调节，即自动发电控制（AGC）。（1）其作用第一，响应时间较慢，能有效地调整分钟级及更长周期的负荷波动；第二，实现频率的无差调整。（2）特点第一，采用的调整方式对系统频率是无差的；第二，响应比较慢，一般需要12分钟；第三，对机组管理往往是比例分配，是发电机组偏离经济运行点。综上所述保证输电系统的频率对于单个用户和发电厂都有重要的意义。电力系统的稳定和安全需要频率的一次调整和二次调整。一次调整靠调速器完成。二次调整由系统中的调频机组实现，通过调频器控制。既可为有差调节，也可以做到无差调节。选择的调频厂有主调频厂和辅助调频厂之分并且满足第一，具有足够的容量。第二，具有较快的调整速
14、度。第三，调整范围内的经济性好。对于具有多个联络点和发电机组的实际电力系统，则AGC将变为包含许多并联发电机组控制回路的形式，如图21所示。其内部控制回路和外部控制回路的基本结构并未改变。G1、G2、G3为发电机组；AGC称为区域内控制误差，用来根据系统频率偏差以及输电线路功率偏差来确定控制信号；负荷分配器根据输入的控制信号大小并且根据等微增率准则或其他原则来控制各台发电机输出功率的大小。自动发电控制系统包括两大部分（1）负荷分配器。根据电力系统频率和其他有关测量信号，按照一定的调节控制准则确定各发电机组的最佳设定输出功率。（2）发电机组控制器。根据负荷分配器所确定的各发电机组最佳输出功率，控
15、制调速器的调节特性，使发电机组在电力系统额定频率下所发出的实际功率与本科生课程设计（论文）5设定的输出功率相一致。图21具有多台发电机的AGC系统自动发电控制系统中的负荷分配配齐使根据所测量的发电机实际输出功率和频率偏差等信号按照一定的准则分配各台发电机组输出功率。决定各台发电机组设定的功率CIP的负荷分配器，目前广泛采用以“基点经济功率BIP”和“分配系数I”来表示每台发电机组的输出功率的方法，即每台发电机组的设定调整功率按以下公式分配11NINIBIGIIBICIPAGCPPP（23）式中CIP各台发电机组的设定调整功率BIP各台发电机的基点经济功率；GIP每台发电机的实际输出功率；I分配
16、系数。也就是说，系统各台发电机组的设定功率，取决于系统发电机组总的实际输出GIP和每日台发电机组的基点经济功率BIP，以及系统频率偏差和功率偏差AEC。偏差越大，各大电机组的设定调整功率的变动就越大。当频率偏差和功率偏差趋于零时，AGC系统发电机组总的设定调整功率就与发电机总的实际输出功率相等。分配到每台发电机组的设定功率值则有分配系数I来决定。这种方法把自动调频N1SN2SN3S电力系统负荷分配器区域控制误差ACEKSG1G2G3PFPG3PG1PG2PC3PC2PC1本科生课程设计（论文）6与经济功率分配联系起来了。其中BIP和I的值可以在每次经济分配计算时加以修正。22单台机组有功控制的
17、基本方法最简单的AGC系统的结构如图22所示，它是具有一台发电机组和联络的AGC系统。图中ZDP为输电线路功率的整定值；ZDF为系统频率整定值；P为触电线路功率的实际值；F为系统频率的实际值；FB为频率修正系数；KS为外部控制回路，用来根据电力系统频率偏差和输电线路上的功率偏差来确定输出控制信号；CP为系统要求调整的控制信号功率；NS为内部控制回路，用来控制调整调速器阀门开度，以达到所需要的输出功率。图22单台发电机组的AGC系统23负荷变化时的功率分配计算课程设计基本参考设计参数如下（1）发电厂有三机组，一号机组功率为50MW，二号机组功率为100MW，三号机组功率为70MW，功率因数均为0
18、9。（2）有功功率调差系数为006。（3）负荷频率调节效应系数（有功功率）为15。机组调速器输电线路NSKSBF电力系统FZDPZDAGC系统PPGPC本科生课程设计（论文）7（4）各发电机均以80负荷运行。系统负荷增加时，经过频率的一次调整，频率由EF降为LF，由发电机组的静态调节公式0FPG，（24）因此可得出1221PP（25）表明并列运行的发电机组之间的功率分配与调差系数成发比关系（标幺值），与单位调节功率成正比。实际并联运行的发电机组的调速器均为有差调节，由其共同承担负荷的波动。假设系统由N台机组运行，则连理调解方程式（1）有铭值,2,10NIFPGIINIINIGIGFPP111（
19、26）按照调差系数的定义，可定义等效发电机组的调差系数为NIIGGPF11127等效发电机组的单位调节功率为NIGINIGGGKK111128这样系统的三台发电机组同样可以当做单独一台发电机对待，满足调节公式。同时，考虑负荷调节效应后，全系统负荷的变化量为FKKFKPPPLGGLGL29系统的调节功率LGLSKKFPK210物理意义表示全系统的频率没变化1HZ时，其调节的负荷有功功率大小。（2）标幺值由于标幺值的基准是不一致的，故采用标幺值要涉及基准值的转换。本科生课程设计（论文）8发电机组以自身的额定容量为GEP负荷以额定负荷容量为基准LEP全系统以负荷额定容量为基准LEP等效发电机组的单位
20、调节功率NIGEGEIGIGGPPKFPK1211系统单位调节功率LRGLSKKKFPK212其中LEGERPPK称为系统热备用系数，表示系统的备用的发电容量的多少。这样采用系统的单位调节功率，可以算出发电出力以及负荷调节效应后的实际功率。根据任务书给的参数可进行以下计算PP1P2P35010070220（MW）213将标幺值化为有名值，等效发电机的单位调节功率为/37350060/220HZMWFPKEG214假设系统的额定负荷为500MW荷的调节效应系数为/2050/2500/HZMWFKPKELLEL215当负荷增加200MW时，频率变化量为22222070/200/HZKKPFLG21
21、6系统频率为77847222250HZFFFE217等效发电机组发出的功率增加量为54155702222MWKFPGG218系统的实际负荷为5465554155500MWPPPGLEL219本科生课程设计（论文）9第3章自动发电系统硬件设计31自动发电系统功能自动发电控制（AGC）就是通过监视电厂出力和系统负荷之间的差异，来控制调频机组的出力，以满足不断变化的用户电力需要，达到电能的发供平衡，并且使整个系统处于经济的运行状态。在联合电力系统中，AGC是以区域系统为单位。AGC能实现机组出力的自动调节，是电力调度EMS系统中最重要的控制功能。在正常的系统运行状态下，AGC的基本功能是（1）使发电
22、机组自动跟踪电力系统负荷变化；（2）响应负荷和发电的随机变化，维持电力系统频率为额定值（50HZ）；（3）在各区域间分配系统发电功率，维持区域间净交换功率为计划值；（4）对周期性的负荷变化按发电计划调整发电功率；监视和调整备用容量，满足电力系统安全要求。32自动发电总体设计方案图31如图31，自动发电控制系统是建立在SCADA系统之上，且与网络分析等应用软件有着密切关系的独立的控制系统。为实现完整的自动发电控制系统，需要建立一套与之相适应的数据库管理系统来实现其各种功能。自动发电控制系统的数据库是能量管理级应用软件的共用数据库。它包括实时发电控制、机组发电计划、机组组合、联络线交易计划、燃料计
23、划、检修计划等。当SCADA传来负荷改变的实时信息时，经过A/D转换器，将信号传给单片发电机组AT89C51AGC时钟电路复位电路SCADA实时信息A/DD/A本科生课程设计（论文）10机89C51。然后单片机89C51经过程序，将信号经过D/A转换器传给AGC，最后通过AGC区域控制模块，将区域控制设定功率与区域当前的发电功率相比较，按照AGC程序设定的控制参数，重新分配各电厂控制器的目标功率。33单片机最小系统设计33189C51单片机引脚功能（1）XTAL1接外部晶体和微调电容的一端。在89C51片内，它是振荡电路反向放大器的输入端及内部时钟发生器的输入端，振荡电路的频率就是晶体的固有频
24、率。当采用外部振荡器是次引脚输入外部时钟脉冲。（2）XTAL2解外部晶体和微调电容的另一端。在89C51片内，它是振荡电路反向放大器的输出端。在采用外部振荡器时此引脚应悬浮。通过示波器查看XTAL2端是否有脉冲信号输出，可以确认89C51的振荡电路是否正常工作。（3）RST复位信号输入端，高电平有效。当振荡器工作时，在此引脚上出现两个机器周期一上的高电平，就可以使单片机复位。（4）ALE/PROG地址锁存允许信号。ALE锁存P0口传送的低8位地址信号，实现低8位地址与数据的分离。（5）PESN外部程序存储器的读选通信号。当89C51由外部程序存储器取指令（或常数）时，每个机器周期内PESN两次
25、有效输出。（6）EA/VPP内，外ROM选择端。当EA端接高电平时，CPU访问并执行内部ROM的指令；但当PC值超过4KB时，将自动转去执行外部ROM中的程序。但端接低电平时，CPU只访问外部ROM中的指令。（7）P0口双向8位三态I/O口，在访问外部存储器时，可分时用做低8位地址线和8位数据线。无上拉电阻，能驱动8个LSTTL门电路。P1口8位双向I/O口，用做普通I/O口。有上拉电阻，能驱动4个LSTTL门电路。P2口8位双向I/O口，做高8位地址线。有上拉电阻，能驱动4个LSTTL门电路。P3口8位双向I/O口，具有第二功能。有上拉电阻，能驱动4个LSTTL门电路。89C51单片机的引脚
26、图（40脚双列直插封装），如图32。EA本科生课程设计（论文）11图3289C51单片机的引脚图332复位电路设计复位电路是使单片机的CPU或系统中的其他部件处于某一确定的初始状态，并从这上状态开始工作，除了进入系统的正常初始化之外，当由于程序运行出现错误或操作错误使系统处于死锁状态时，为了摆脱困境，也需按复位电路以重新启动。复位电路包括上电复位，按键电平复位，按键脉冲复位。本设计中采用按键电平复位。按键电平复位是通过是复位端经过电阻与VCC电源接通而实现的，如下图S1SWPBC122UFR21KR1200VCCVCCRESETVSS图33按键电平复位电路复位信号及其产生RST引脚是复位信号的
27、输入端，复位信号时高电平有效，其有效时间持续24个振荡脉冲周期（即2个机器周期）以上，若使用频率为12MHZ本科生课程设计（论文）12的晶振，则复位信号持续时间应超过2US才能完成复位操作。整个复位电路包括片内外两部分，外部电路产生的复位电路送施密特触发器，再由片内复位电路在每个机器周期的S5P2时刻对施密特触发器的输出进行采样。然后才得到内部复位操作所需要的信号。333时钟电路设计时钟电路由一个晶体振荡器12MHZ和两个33PF的瓷片电容组成。时钟电路产生单片机工作所需要的时钟信号，而时序所研究的是指令执行中各信号之间的相互关系。单片机本身就如一个复杂的同步时序电路，为了保证同步工作方式的实
28、现，电路应在唯一的时钟信号控制下严格工作。单片机芯片内部有一个高增益反相放大器，其输入端为芯片引脚XTAL1，输出端为引脚XTAL2，在芯片外部通过两个引脚跨接晶体振荡器和微调电容，形成反馈电路，就构成了一个稳定的自激振荡器，如下图C133UFC233UF晶振12MHZXTAL2XTAL1图34时钟电路334直流稳压电源设计T11234D1C13300UFC2100UFC3100UFC4100UFW7805VINVOUT123GND图35直流稳压电源本科生课程设计（论文）13电源电路设计的要求是将220V交流电转化为5V的直流电，如下图所示变压器将220V交流电变换成12V交流电，经过整流桥又
29、变换为12V的直流电，三断集成稳压器W7805的作用是将12V的直流电转化成5V的直流电，图中的4个电容是滤波电容。集成串联型稳压电路有三个脚，分别为输入端、输出端和公共端，因此称为三端稳压器。按功能分为固定式和可调节式稳压电路；前者输出地电压不能进行调节，为固定值；后者可通过外接元件使输出电压得到很宽的调节范围。本设计中使用的W7805三端稳压器为固定式稳压器，型号最后两位数为输出电压值，所以W7805的输出电压为5V。335单片机最小系统P101P112P123P134P145P156P167RESET9P3010P3212P3313P3414P3515P3616P3717XTAL218X
30、TAL119VSS20P2021P2122P2223P2324P2425P2526P2627P2728PSEN29ALE/PROG30EA/VPP31P0732P0633P0534P0435P0336P0237P0138P0039VCC40P178P3111AT89C51C1C2XTAL2XTAL1R1200R21KC22UFRESETVCCVCCRESETVSS图36单片机最小系统本科生课程设计（论文）1434自动发电控制原理图设计C233UFC333UFY112MHZSWPBC122UFR1200R21KVCCGNDGNDRDWRP26P0032P1237P1138P1039ALE30EA
31、3189C51VCC40RESET9VSS20XTAL218XTAL119AGCSCADA89C51指令1指令2指令3图37输入输处接口图如图37所示,为自动发电控制原理原理图。SCADA实时监测发电机组的运行状态，不断地向单片机89C51发送信号，当发电机组运行状态改变时，SCADA发送改变信号，单片机接受信号并且辨别信号的内容，经过程序辨别后向AGC发出执行指令类别信号。再由AGC向发电机组发送改变运行状态的指令。本科生课程设计（论文）15第4章自动发电系统软件设计41软件功能当发电机组1的负荷改变时，SCADA经过A/D向单片机AT89C51发一个信号，AT89C51经过程序然后通过P1
32、0通道经过D/A向AGC发送执行指令，AGC再向发电机组发送执行指令。当发电机组2的负荷改变时，SCADA经过A/D向单片机AT89C51发一个信号，单片机通过P11向AGC发执行指令，AGC在向发电机组发出调节的指令。同理，在发电机组3的负荷改变时SCADA向单片机发信号，单片机经过P12向AGC发送信号，AGC向发电机组调节指令。42软件流程图当单片机接受到频率信号时，通过程序先识别频率是否改变，如果频率没有改变，说明发电机组中的负荷没有发生改变。如果频率发生了改变，程序再继续识别是否是发电机组1负荷发生了改变，如果是，则执行指令1改变发电机组1的运行状态，如果不是发电机组1，则在继续识别
33、。发电机组2、发电机组3与发电机组1一样。主程序流程图如图41开始初始化频率检测频率是否改变指令1指令2指令3结束NY图41软件流程图发电机1发电机2发电机3YY本科生课程设计（论文）16第5章课程设计总结本次设计阐述了发电机有功功率的基本原理，分析了系统负荷一定和变化时机组有功功率调差系数和负荷频率调节效应系数对功率分配的影响，并推倒相关有功功率的分配公式。设计的同时，分析了AGC控制系统控制原理和AT89C51单片机原理，在SCADA平台上对时钟电路和复位电路进行设计，阐述了AT89C51各引脚功能并对AT89C51单片机最小系统进行讨论和分析设计。通过对复位电路、时钟电路、AT89C51
34、、AGC等小系统的原理进行分析和讨论，最终选出最优方案，完成本次设计。本次设计由于时间和条件限制，还有很多地方需要改进，比如AGC和SCADA系统的理解不足，原理运用不到位等。但是总体完全满足设计要求。同时，本次设计还得到了几位老师的指导，比如单片机的软件设计，控制功能的完善等。在此特别感谢本科生课程设计（论文）17参考文献1方富淇等配电网自动化北京中国电力出版社，20002赵晶PROTE199高级应用人民邮电出版社，20003何仰赞等电力系统分析华中科技大学出版社，200234于海生微型计算机控制技术清华大学出版社，200345王士政电网调度自动化与配网自动化技术中国水利水电出版社，200736梅丽凤等单片机原理及接口技术清华大学出版社，200977李广弟单片机基础。北京航空航天大学出版社，19948郭培源等电力系统自动控制技术北京科学出版社，20019方富淇等配电网自动化北京中国电力出版社，200010商国才电力系统自动化天津大学出版社，2000
