CN108062440A - 先进绝热压缩空气储能系统全动态模型的生成方法及装置 - Google Patents
先进绝热压缩空气储能系统全动态模型的生成方法及装置 Download PDFInfo
- Publication number
- CN108062440A CN108062440A CN201711318890.0A CN201711318890A CN108062440A CN 108062440 A CN108062440 A CN 108062440A CN 201711318890 A CN201711318890 A CN 201711318890A CN 108062440 A CN108062440 A CN 108062440A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- operating mode
- msubsup
- mfrac
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
Abstract
本发明提出了一种先进绝热压缩空气储能系统全动态模型的生成方法及装置,其中方法包括:建立空气压缩热能和压力势能的解耦生产单元和耦合释能单元的压缩机和膨胀机宽工况动态模型;建立压缩热能传输收集单元和压缩热能传输释放单元的换热系统宽工况动态模型;建立压缩热能和压力势能存储单元的储能系统宽工况动态模型;建立先进绝热压缩空气储能与电网耦合接口的联系,并且构建从系统组件动态特性到系统整体性能间的关联关系,以为智能电网不同场景应用及能效提升提供模型支撑。该方法可以充分计及相关组件宽工况运行特性,并建立全动态模型,从而可以适当简化并建立面向智能电网辅助服务等应用场景的备用模型等,具有很强的应用价值。
Description
技术领域
本发明涉及大规模物理储能技术领域,尤其涉及一种先进绝热压缩空气储能系统全动态模型的生成方法及装置。
背景技术
AA-CAES(Advanced Adiabatic Compressed Air Energy Storage,先进绝热压缩空气储能系统)是一种摒弃燃料补燃的清洁大规模物理储能技术,被视为继抽蓄储能之后最具吸引力的大规模物理储能技术。AA-CAES主要应用于智能电网削峰填谷、频率调节、旋转备用、无功支撑、黑启动等场景。AA-CAES电站利用弃风(光)、低谷电等低品位电能驱动压缩机,通过绝热压缩回收压缩热,解耦存储空气压力势能(储气室)和压缩热能(储热系统),再通过膨胀利用压缩热,实现空气压力势能和压缩热能的耦合释能,具有效率高、成本低等特点,其理论电-电转换效率可达70%左右。AA-CAES系统包含压缩机、换热器、膨胀机、储热系统、储气室等组件,各组件功-能转换高度耦合。AA-CAES动态模型是研究和衡量各组件性能对系统电-电转换效率等整体指标影响程度的关键,如何构建AA-CAES动态模型极为必要。
一般采用基于固定效率模型来设计、描述上述组件的功-能转换关系。然而,智能电网削峰填谷、旋转备用、黑启动、无功支撑等应用场景要求AA-CAES具有宽工况运行能力(20%-110%压缩/发电工况运行),在宽工况运行条件下,压缩机(压力势能和压缩热能解耦生产单元)、膨胀机(压力势能和压缩热能耦合释能单元)效率变化明显,如透平在额定工况下可达85%-90%效率,而在50%工况下效率将至65%-75%。换热及储热系统的存在是AA-CAES不同于传统CAES的重要特点,相应地,换热器(压缩热能传输收集及回馈单元)的传热系数受到热工质流速(压缩、发电负载情况)的影响,不同流速下传热系数变化明显,影响压缩热能的收集和回馈,进而影响系统整体电-电转换性能。因此,当前采用的固定效率模型难以计及压缩机、膨胀机、换热系统等组件宽工况运行对系统整体性能的影响,进而在AA-CAES的智能电网应用场景中受限。
发明内容
本发明旨在至少在一定程度上解决相关技术中的技术问题之一。
为此,本发明的一个目的在于提出一种先进绝热压缩空气储能系统全动态模型的生成方法,该方法可以适当简化并建立面向智能电网辅助服务等应用场景的备用模型等,高度契合AA-CAES工程运行特性,具有很强的应用价值。
本发明的另一个目的在于提出一种先进绝热压缩空气储能系统全动态模型的生成装置。
为达到上述目的,本发明一方面实施例提出了一种先进绝热压缩空气储能系统全动态模型的生成方法,包括以下步骤:根据压缩机和膨胀机的宽工况运行特性,建立空气压缩热能和压力势能的解耦生产单元和耦合释能单元的压缩机和膨胀机宽工况动态模型;根据压缩侧换热系统和膨胀机侧换热系统的宽工况运行特性,建立压缩热能传输收集单元和压缩热能传输释放单元的换热系统宽工况动态模型;根据储热罐温度流量和储气罐压力温度的动态方程,建立压缩热能和压力势能存储单元的储能系统宽工况动态模型;以及基于电动机及发电机动态模型,建立先进绝热压缩空气储能与电网耦合接口的联系,并且构建从系统组件动态特性到系统整体性能间的关联关系,以为智能电网不同场景应用及能效提升提供模型支撑。
本发明实施例的先进绝热压缩空气储能系统全动态模型的生成方法,可以充分计及压缩机、膨胀机、换热系统、储热系统、储气系统等组件宽工况运行特性,并建立AA-CAES全动态模型,进而刻画组件动态性能到系统整体电电转换效率等指标间的联系,并且可以衡量外部热源等对系统性能的影响,进而为AA-CAES能效提升提供分析方法,从而可以适当简化并建立面向智能电网辅助服务等应用场景的备用模型等,高度契合AA-CAES工程运行特性,具有很强的应用价值。
另外,根据本发明上述实施例的先进绝热压缩空气储能系统全动态模型的生成方法还可以具有以下附加的技术特征:
进一步地,在本发明的一个实施例中,所述压缩机的透平宽工况运行特性描述为:
其中,βc,ηc为压缩机的压缩比与等熵效率,Gc,nc分别表示宽工况运行时的压缩机质量流量与工质流速;c相应系数,a1,a2,a3为常数;
并且,类似压缩机的透平宽工况运行特性为:
其中,ηt为透平等熵效率,Gt,nt分别表示宽工况运行时的透平质量流量与工质流速。
进一步地,在本发明的一个实施例中,其中,
并且,
其中,Tc,in为实际运行工况对应的压缩机进口温度,Pc,in为实际运行工况对应的压缩机进口压力、mc为实际运行工况对应的压缩机进口质量流量,为额定设计工况下的压缩机进口温度、为额定设计工况下的压缩机进口压力,为额定设计工况下的压缩机进口质量流量,Tt,in为实际运行工况对应的透平进口温度,Pt,in为实际运行工况对应的透平进口压力,mt为实际运行工况对应的透平进口质量流量,额定设计工况下的为额定设计工况下的透平进口温度、为额定设计工况下的透平进口压力,为额定设计工况下的透平进口质量流量,b1,b2分别为相应系数。
进一步地,在本发明的一个实施例中,所述换热系统宽工况动态模型为:
其中,ΦHEX为换热器实际换热量,Cmin为最小热容,代表换热器两侧冷热流体入口温度,ε为由换热器换热系数、换热面积等界定的系数,受运行工况影响。
进一步地,在本发明的一个实施例中,所述动态方程为:
其中,Tw,hw,Aw分别为环境温度、换热系数及接触面积,mr、Tr为储气库高压空气质量及温度,cp,a为高压空气定压比热容,k为比热系数。
本发明另一方面实施例提出了一种先进绝热压缩空气储能系统全动态模型的生成装置,包括:第一构建模块,用于根据压缩机和膨胀机的宽工况运行特性,构建空气压缩热能和压力势能的解耦生产单元和耦合释能单元的压缩机和膨胀机宽工况动态模型;第二构建模块,用于根据压缩侧换热系统和膨胀机侧换热系统的宽工况运行特性,构建压缩热能传输收集单元和压缩热能传输释放单元的换热系统宽工况动态模型;第三构建模块,用于根据储热罐温度流量和储气罐压力温度的动态方程,构建压缩热能和压力势能存储单元的储能系统宽工况动态模型;以及生成模块,用于基于电动机及发电机动态模型,建立先进绝热压缩空气储能与电网耦合接口的联系,并且构建从系统组件动态特性到系统整体性能间的关联关系,以为智能电网不同场景应用及能效提升提供模型支撑。
根据本发明实施例的先进绝热压缩空气储能系统全动态模型的生成装置,可以充分计及压缩机、膨胀机、换热系统、储热系统、储气系统等组件宽工况运行特性,并建立AA-CAES全动态模型,进而刻画组件动态性能到系统整体电电转换效率等指标间的联系,并且可以衡量外部热源等对系统性能的影响,进而为AA-CAES能效提升提供分析方法,从而可以适当简化并建立面向智能电网辅助服务等应用场景的备用模型等,高度契合AA-CAES工程运行特性,具有很强的应用价值。
另外,根据本发明上述实施例的先进绝热压缩空气储能系统全动态模型的生成装置还可以具有以下附加的技术特征:
进一步地,在本发明的一个实施例中,所述压缩机的透平宽工况运行特性描述为:
其中,βc,ηc为压缩机的压缩比与等熵效率,Gc,nc分别表示宽工况运行时的压缩机质量流量与工质流速,c相应系数,a1,a2,a3为常数;
并且,类似压缩机的透平宽工况运行特性为:
其中,ηt为透平等熵效率,Gt,nt分别表示宽工况运行时的透平质量流量与工质流速。
进一步地,在本发明的一个实施例中,其中,
并且,
其中,Tc,in为实际运行工况对应的压缩机进口温度,Pc,in为实际运行工况对应的压缩机进口压力、mc为实际运行工况对应的压缩机进口质量流量,为额定设计工况下的压缩机进口温度、为额定设计工况下的压缩机进口压力,为额定设计工况下的压缩机进口质量流量,Tt,in为实际运行工况对应的透平进口温度,Pt,in为实际运行工况对应的透平进口压力,mt为实际运行工况对应的透平进口质量流量,额定设计工况下的为额定设计工况下的透平进口温度、为额定设计工况下的透平进口压力,为额定设计工况下的透平进口质量流量,b1,b2分别为相应系数,a0为。
进一步地,在本发明的一个实施例中,所述换热系统宽工况动态模型为:
其中,ΦHEX为换热器实际换热量,Cmin为最小热容,代表换热器两侧冷热流体入口温度,ε为由换热器换热系数、换热面积等界定的系数,受运行工况影响。
进一步地,在本发明的一个实施例中,所述动态方程为:
其中,Tw,hw,Aw分别为环境温度、换热系数及接触面积,mr、Tr为储气库高压空气质量及温度,cp,a为高压空气定压比热容,k为比热系数。
本发明附加的方面和优点将在下面的描述中部分给出,部分将从下面的描述中变得明显,或通过本发明的实践了解到。
附图说明
本发明上述的和/或附加的方面和优点从下面结合附图对实施例的描述中将变得明显和容易理解,其中:
图1是根据本发明实施例的先进绝热压缩空气储能系统全动态模型的生成方法的流程图;
图2是根据本发明实施例的先进绝热压缩空气储能系统全动态模型的生成装置的结构示意图。
具体实施方式
下面详细描述本发明的实施例,所述实施例的示例在附图中示出,其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的,旨在用于解释本发明,而不能理解为对本发明的限制。
下面参照附图描述根据本发明实施例提出的先进绝热压缩空气储能系统全动态模型的生成方法及装置,首先将参照附图描述根据本发明实施例提出的先进绝热压缩空气储能系统全动态模型的生成方法。
图1是根据本发明实施例的先进绝热压缩空气储能系统全动态模型的生成方法的流程图。
如图1所示,该先进绝热压缩空气储能系统全动态模型的生成方法包括以下步骤:
在步骤S101中,根据压缩机和膨胀机的宽工况运行特性,建立空气压缩热能和压力势能的解耦生产单元和耦合释能单元的压缩机和膨胀机宽工况动态模型。
可以理解的是,首先,计及压缩机、膨胀机宽工况运行导致的效率变化问题,建立压缩机、膨胀机等空气压缩热能和压力势能解耦生产单元和耦合释能单元的动态模型。
进一步地,在本发明的一个实施例中,压缩机的透平宽工况运行特性描述为:
其中,βc,ηc为压缩机的压缩比与等熵效率,Gc,nc分别表示宽工况运行时的压缩机质量流量与工质流速,c相应系数,a1,a2,a3为常数;
并且,类似压缩机的透平宽工况运行特性为:
其中,ηt为透平等熵效率,Gt,nt分别表示宽工况运行时的透平质量流量与工质流速。
可选地,在本发明的一个实施例中,其中,
并且,
其中,Tc,in为实际运行工况对应的压缩机进口温度,Pc,in为实际运行工况对应的压缩机进口压力、mc为实际运行工况对应的压缩机进口质量流量,为额定设计工况下的压缩机进口温度、为额定设计工况下的压缩机进口压力,为额定设计工况下的压缩机进口质量流量,Tt,in为实际运行工况对应的透平进口温度,Pt,in为实际运行工况对应的透平进口压力,mt为实际运行工况对应的透平进口质量流量,额定设计工况下的为额定设计工况下的透平进口温度、为额定设计工况下的透平进口压力,为额定设计工况下的透平进口质量流量,b1,b2分别为相应系数。
具体而言,压缩机、透平宽工况运行特性
其中,βc,ηc为压缩机的压缩比与等熵效率;Gc,nc分别表示宽工况运行时的压缩机质量流量与工质流速,可由实际运行工况对应的压缩机进口温度Tc,in、进口压力Pc,in、进口质量流量mc与额定设计工况下的及确定;b1,b2,c分别为相应系数,由下式确定,其中a1,a2,a3为常数:
类似压缩机,透平宽工况运行特性可描述为:
其中,ηt为透平等熵效率,Gt,nt分别表示宽工况运行时的透平质量流量与工质流速,可由实际运行工况对应的透平进口温度Tt,in、进口压力Pt,in、进口质量流量mt与额定设计工况下的及确定,具体如下式所示:
在步骤S102中,根据压缩侧换热系统和膨胀机侧换热系统的宽工况运行特性,建立压缩热能传输收集单元和压缩热能传输释放单元的换热系统宽工况动态模型。
可以理解的是,本发明实施例可以计及压缩侧换热系统和膨胀机侧换热系统由于宽工况运行引起的传热系数变化问题,构建压缩热能传输收集单元和压缩热能传输释放单元的动态模型。
进一步地,在本发明的一个实施例中,换热系统宽工况动态模型为:
其中,ΦHEX为换热器实际换热量,Cmin为最小热容,代表换热器两侧冷热流体入口温度,ε为由换热器换热系数、换热面积等界定的系数,受运行工况影响。
可以理解的是,换热系统宽工况动态模型为:
其中,ΦHEX为换热器实际换热量,Cmin为最小热容,代表换热器两侧冷热流体入口温度,ε为由换热器换热系数、换热面积等界定的系数,受运行工况影响。
在步骤S103中,根据储热罐温度流量和储气罐压力温度的动态方程,建立压缩热能和压力势能存储单元的储能系统宽工况动态模型。
可以理解的是,本发明实施例可以计及储热罐温度流量、储气罐压力温度动态方程,构建压缩热能和压力势能存储单元的动态模型。
进一步地,在本发明的一个实施例中,动态方程为:
其中,Tw,hw,Aw分别为环境温度、换热系数及接触面积,mr、Tr为储气库高压空气质量及温度,cp,a为高压空气定压比热容,k为比热系数。
可以理解的是,计及储气系统与周围环境换热特性,则储气库动态模型可由能量平衡和质量平衡动态方程描述,即:
其中,Tw,hw,Aw分别为环境温度、换热系数及接触面积,mr、Tr为储气库高压空气质量及温度,cp,a为高压空气定压比热容,k为比热系数。
在步骤S104中,基于电动机及发电机动态模型,建立先进绝热压缩空气储能与电网耦合接口的联系,并且构建从系统组件动态特性到系统整体性能间的关联关系,以为智能电网不同场景应用及能效提升提供模型支撑。
可以理解的是,本发明实施例可以计及电动机及发电机动态模型,建立先进绝热压缩空气储能与电网耦合接口的联系,构建从系统组件动态特性到系统整体性能间的关联关系,为其智能电网不同场景应用及能效提升提供模型支撑。
综上,本发明实施例在经典常效率模型的基础上计及压缩机、膨胀机、换热器等组件因宽工况运行引起的效率变化特性,进而构建从AA-CAES组件性能衡量系统整体电-电转换效率的全动态模型,并可以在动态组件性能和系统整体电电转换效率等指标间建立联系,进而可以衡量引入外部热源(光热系统)等对系统整体性能提升的效益,建立起面向智能电网辅助服务等应用场景的备用模型等,同时也可适当简化,建立起面向智能电网辅助服务等应用场景的备用模型等,符合工程实际,具有很强的应用价值。并且,本发明实施例的方法契合先进绝热压缩空气储能系统工程实际宽工况运行要求,具有广泛的应用价值。
根据本发明实施例的先进绝热压缩空气储能系统全动态模型的生成方法,可以充分计及压缩机、膨胀机、换热系统、储热系统、储气系统等组件宽工况运行特性,并建立AA-CAES全动态模型,进而刻画组件动态性能到系统整体电电转换效率等指标间的联系,并且可以衡量外部热源等对系统性能的影响,进而为AA-CAES能效提升提供分析方法,从而可以适当简化并建立面向智能电网辅助服务等应用场景的备用模型等,高度契合AA-CAES工程运行特性,具有很强的应用价值。
其次参照附图描述根据本发明实施例提出的先进绝热压缩空气储能系统全动态模型的生成装置。
图2是本发明一个实施例的先进绝热压缩空气储能系统全动态模型的生成装置的结构示意图。
如图2所示,该先进绝热压缩空气储能系统全动态模型的生成装置10包括:第一构建模块100、第二构建模块200、第三构建模块300和生成模块400。
其中,第一构建模块100用于根据压缩机和膨胀机的宽工况运行特性,构建空气压缩热能和压力势能的解耦生产单元和耦合释能单元的压缩机和膨胀机宽工况动态模型。第二构建模块200用于根据压缩侧换热系统和膨胀机侧换热系统的宽工况运行特性,构建压缩热能传输收集单元和压缩热能传输释放单元的换热系统宽工况动态模型。第三构建模块300用于根据储热罐温度流量和储气罐压力温度的动态方程,构建压缩热能和压力势能存储单元的储能系统宽工况动态模型。生成模块400用于基于电动机及发电机动态模型,建立先进绝热压缩空气储能与电网耦合接口的联系,并且构建从系统组件动态特性到系统整体性能间的关联关系,以为智能电网不同场景应用及能效提升提供模型支撑。本发明实施例的装置10可以充分计及压缩机、膨胀机、换热系统、储热系统、储气系统等组件宽工况运行特性,并建立AA-CAES全动态模型,从而可以适当简化并建立面向智能电网辅助服务等应用场景的备用模型等,高度契合AA-CAES工程运行特性,具有很强的应用价值。
进一步地,在本发明的一个实施例中,压缩机的透平宽工况运行特性描述为:
其中,βc,ηc为压缩机的压缩比与等熵效率,Gc,nc分别表示宽工况运行时的压缩机质量流量与工质流速,c相应系数,a1,a2,a3为常数;
并且,类似压缩机的透平宽工况运行特性为:
其中,ηt为透平等熵效率,Gt,nt分别表示宽工况运行时的透平质量流量与工质流速。
进一步地,在本发明的一个实施例中,其中,
并且,
其中,Tc,in为实际运行工况对应的压缩机进口温度,Pc,in为实际运行工况对应的压缩机进口压力、mc为实际运行工况对应的压缩机进口质量流量,为额定设计工况下的压缩机进口温度、为额定设计工况下的压缩机进口压力,为额定设计工况下的压缩机进口质量流量,Tt,in为实际运行工况对应的透平进口温度,Pt,in为实际运行工况对应的透平进口压力,mt为实际运行工况对应的透平进口质量流量,额定设计工况下的为额定设计工况下的透平进口温度、为额定设计工况下的透平进口压力,为额定设计工况下的透平进口质量流量,b1,b2分别为相应系数。
进一步地,在本发明的一个实施例中,换热系统宽工况动态模型为:
其中,ΦHEX为换热器实际换热量,Cmin为最小热容,代表换热器两侧冷热流体入口温度,ε为由换热器换热系数、换热面积等界定的系数,受运行工况影响。
进一步地,在本发明的一个实施例中,动态方程为:
其中,Tw,hw,Aw分别为环境温度、换热系数及接触面积,mr、Tr为储气库高压空气质量及温度,cp,a为高压空气定压比热容,k为比热系数。
需要说明的是,前述对先进绝热压缩空气储能系统全动态模型的生成方法实施例的解释说明也适用于该先进绝热压缩空气储能系统全动态模型的生成装置,此处不再赘述。
根据本发明实施例的一种先进绝热压缩空气储能系统全动态模型的生成装置,可以充分计及压缩机、膨胀机、换热系统、储热系统、储气系统等组件宽工况运行特性,并建立AA-CAES全动态模型,进而刻画组件动态性能到系统整体电电转换效率等指标间的联系,并且可以衡量外部热源等对系统性能的影响,进而为AA-CAES能效提升提供分析方法,从而可以适当简化并建立面向智能电网辅助服务等应用场景的备用模型等,高度契合AA-CAES工程运行特性,具有很强的应用价值。
在本发明的描述中,需要理解的是,术语“中心”、“纵向”、“横向”、“长度”、“宽度”、“厚度”、“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”“内”、“外”、“顺时针”、“逆时针”、“轴向”、“径向”、“周向”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本发明的限制。
此外,术语“第一”、“第二”仅用于描述目的,而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量。由此,限定有“第一”、“第二”的特征可以明示或者隐含地包括至少一个该特征。在本发明的描述中,“多个”的含义是至少两个,例如两个,三个等,除非另有明确具体的限定。
在本发明中,除非另有明确的规定和限定,术语“安装”、“相连”、“连接”、“固定”等术语应做广义理解,例如,可以是固定连接,也可以是可拆卸连接,或成一体;可以是机械连接,也可以是电连接;可以是直接相连,也可以通过中间媒介间接相连,可以是两个元件内部的连通或两个元件的相互作用关系,除非另有明确的限定。对于本领域的普通技术人员而言,可以根据具体情况理解上述术语在本发明中的具体含义。
在本发明中,除非另有明确的规定和限定,第一特征在第二特征“上”或“下”可以是第一和第二特征直接接触,或第一和第二特征通过中间媒介间接接触。而且,第一特征在第二特征“之上”、“上方”和“上面”可是第一特征在第二特征正上方或斜上方,或仅仅表示第一特征水平高度高于第二特征。第一特征在第二特征“之下”、“下方”和“下面”可以是第一特征在第二特征正下方或斜下方,或仅仅表示第一特征水平高度小于第二特征。
在本说明书的描述中,参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中,对上述术语的示意性表述不必须针对的是相同的实施例或示例。而且,描述的具体特征、结构、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外,在不相互矛盾的情况下,本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。
尽管上面已经示出和描述了本发明的实施例,可以理解的是,上述实施例是示例性的,不能理解为对本发明的限制,本领域的普通技术人员在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。
Claims (10)
1.一种先进绝热压缩空气储能系统全动态模型的生成方法,其特征在于,包括以下步骤:
根据压缩机和膨胀机的宽工况运行特性,建立空气压缩热能和压力势能的解耦生产单元和耦合释能单元的压缩机和膨胀机宽工况动态模型;
根据压缩侧换热系统和膨胀机侧换热系统的宽工况运行特性,建立压缩热能传输收集单元和压缩热能传输释放单元的换热系统宽工况动态模型;
根据储热罐温度流量和储气罐压力温度的动态方程,建立压缩热能和压力势能存储单元的储能系统宽工况动态模型;以及
基于电动机及发电机动态模型,建立先进绝热压缩空气储能与电网耦合接口的联系,并且构建从系统组件动态特性到系统整体性能间的关联关系,以为智能电网不同场景应用及能效提升提供模型支撑。
2.根据权利要求1所述的先进绝热压缩空气储能系统全动态模型的生成方法,其特征在于,所述压缩机的透平宽工况运行特性描述为:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&beta;</mi>
<mi>c</mi>
</msub>
<mo>=</mo>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>G</mi>
<mi>c</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>a</mi>
<mn>2</mn>
</msub>
<msub>
<mi>G</mi>
<mi>c</mi>
</msub>
<mo>+</mo>
<msub>
<mi>a</mi>
<mn>3</mn>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&eta;</mi>
<mi>c</mi>
</msub>
<mo>=</mo>
<mo>&lsqb;</mo>
<mn>1</mn>
<mo>-</mo>
<mi>c</mi>
<msup>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>n</mi>
<mi>c</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>&rsqb;</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>n</mi>
<mi>c</mi>
</msub>
<mo>/</mo>
<msub>
<mi>G</mi>
<mi>c</mi>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>-</mo>
<mo>(</mo>
<mrow>
<msub>
<mi>n</mi>
<mi>c</mi>
</msub>
<mo>/</mo>
<msub>
<mi>G</mi>
<mi>c</mi>
</msub>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>,</mo>
</mrow>
其中,βc,ηc为压缩机的压缩比与等熵效率,Gc,nc分别表示宽工况运行时的压缩机质量流量与工质流速,c相应系数,a1,a2,a3为常数;
并且,类似压缩机的透平宽工况运行特性为:
<mrow>
<msub>
<mi>&eta;</mi>
<mi>t</mi>
</msub>
<mo>/</mo>
<msubsup>
<mi>&eta;</mi>
<mi>t</mi>
<mn>0</mn>
</msubsup>
<mo>=</mo>
<mo>[</mo>
<mn>1</mn>
<mo>-</mo>
<mi>t</mi>
<msup>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>n</mi>
<mi>t</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<msub>
<mi>n</mi>
<mi>t</mi>
</msub>
<mo>/</mo>
<msub>
<mi>G</mi>
<mi>t</mi>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>-</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>n</mi>
<mi>t</mi>
</msub>
<mo>/</mo>
<msub>
<mi>G</mi>
<mi>t</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
其中,ηt为透平等熵效率,Gt,nt分别表示宽工况运行时的透平质量流量与工质流速。
3.根据权利要求2所述的先进绝热压缩空气储能系统全动态模型的生成方法,其特征在于,其中,
<mrow>
<msub>
<mi>G</mi>
<mi>c</mi>
</msub>
<mo>=</mo>
<mo>&lsqb;</mo>
<msub>
<mi>m</mi>
<mi>c</mi>
</msub>
<mfrac>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>T</mi>
<mrow>
<mi>c</mi>
<mo>,</mo>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mn>0.5</mn>
</msup>
<msub>
<mi>P</mi>
<mrow>
<mi>c</mi>
<mo>,</mo>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
</mfrac>
<mo>&rsqb;</mo>
<mo>/</mo>
<mo>&lsqb;</mo>
<msubsup>
<mi>m</mi>
<mi>c</mi>
<mn>0</mn>
</msubsup>
<mfrac>
<msup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>T</mi>
<mrow>
<mi>c</mi>
<mo>,</mo>
<mi>i</mi>
<mi>n</mi>
</mrow>
<mn>0</mn>
</msubsup>
<mo>)</mo>
</mrow>
<mn>0.5</mn>
</msup>
<msubsup>
<mi>P</mi>
<mrow>
<mi>c</mi>
<mo>,</mo>
<mi>i</mi>
<mi>n</mi>
</mrow>
<mn>0</mn>
</msubsup>
</mfrac>
<mo>&rsqb;</mo>
<mo>,</mo>
</mrow>
<mrow>
<msub>
<mi>&eta;</mi>
<mi>c</mi>
</msub>
<mo>=</mo>
<mo>&lsqb;</mo>
<msub>
<mi>n</mi>
<mi>c</mi>
</msub>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>T</mi>
<mrow>
<mi>c</mi>
<mo>,</mo>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>0.5</mn>
</mrow>
</msup>
<mo>&rsqb;</mo>
<mo>/</mo>
<mo>&lsqb;</mo>
<msubsup>
<mi>n</mi>
<mi>c</mi>
<mn>0</mn>
</msubsup>
<msup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>T</mi>
<mrow>
<mi>c</mi>
<mo>,</mo>
<mi>i</mi>
<mi>n</mi>
</mrow>
<mn>0</mn>
</msubsup>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>0.5</mn>
</mrow>
</msup>
<mo>&rsqb;</mo>
<mo>,</mo>
</mrow>
<mrow>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mfrac>
<msub>
<mi>n</mi>
<mi>c</mi>
</msub>
<msub>
<mi>a</mi>
<mn>0</mn>
</msub>
</mfrac>
<mo>,</mo>
<msub>
<mi>a</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>b</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<mn>2</mn>
<msub>
<mi>b</mi>
<mn>2</mn>
</msub>
<msubsup>
<mi>n</mi>
<mi>c</mi>
<mn>2</mn>
</msubsup>
<mo>)</mo>
</mrow>
<msub>
<mi>a</mi>
<mn>0</mn>
</msub>
</mfrac>
<mo>,</mo>
<msub>
<mi>a</mi>
<mn>3</mn>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msup>
<msub>
<mi>b</mi>
<mn>2</mn>
</msub>
<mn>2</mn>
</msup>
<msubsup>
<mi>n</mi>
<mi>c</mi>
<mn>3</mn>
</msubsup>
<mo>-</mo>
<msub>
<mi>b</mi>
<mn>1</mn>
</msub>
<msub>
<mi>b</mi>
<mn>2</mn>
</msub>
<msub>
<mi>n</mi>
<mi>c</mi>
</msub>
</mrow>
<msub>
<mi>a</mi>
<mn>0</mn>
</msub>
</mfrac>
<mo>,</mo>
</mrow>
<mrow>
<msub>
<mi>a</mi>
<mn>0</mn>
</msub>
<mo>=</mo>
<mo>&lsqb;</mo>
<msub>
<mi>b</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mfrac>
<msub>
<mi>b</mi>
<mn>2</mn>
</msub>
<msub>
<mi>n</mi>
<mi>c</mi>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>n</mi>
<mi>c</mi>
</msub>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>n</mi>
<mi>c</mi>
</msub>
<mo>-</mo>
<msub>
<mi>b</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>&rsqb;</mo>
<mo>,</mo>
</mrow>
并且,
<mrow>
<msub>
<mi>G</mi>
<mi>t</mi>
</msub>
<mo>=</mo>
<mo>&lsqb;</mo>
<msub>
<mi>m</mi>
<mi>t</mi>
</msub>
<mfrac>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>T</mi>
<mrow>
<mi>t</mi>
<mo>,</mo>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mn>0.5</mn>
</msup>
<msub>
<mi>P</mi>
<mrow>
<mi>t</mi>
<mo>,</mo>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
</mfrac>
<mo>&rsqb;</mo>
<mo>/</mo>
<mo>&lsqb;</mo>
<msubsup>
<mi>m</mi>
<mi>t</mi>
<mn>0</mn>
</msubsup>
<mfrac>
<msup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>T</mi>
<mrow>
<mi>t</mi>
<mo>,</mo>
<mi>i</mi>
<mi>n</mi>
</mrow>
<mn>0</mn>
</msubsup>
<mo>)</mo>
</mrow>
<mn>0.5</mn>
</msup>
<msubsup>
<mi>P</mi>
<mrow>
<mi>t</mi>
<mo>,</mo>
<mi>i</mi>
<mi>n</mi>
</mrow>
<mn>0</mn>
</msubsup>
</mfrac>
<mo>&rsqb;</mo>
<mo>,</mo>
</mrow>
其中,Tc,in为实际运行工况对应的压缩机进口温度,Pc,in为实际运行工况对应的压缩机进口压力、m c为实际运行工况对应的压缩机进口质量流量,为额定设计工况下的压缩机进口温度、为额定设计工况下的压缩机进口压力,为额定设计工况下的压缩机进口质量流量,Tt,in为实际运行工况对应的透平进口温度,Pt,in为实际运行工况对应的透平进口压力,mt为实际运行工况对应的透平进口质量流量,额定设计工况下的为额定设计工况下的透平进口温度、为额定设计工况下的透平进口压力,为额定设计工况下的透平进口质量流量,b1,b2分别为相应系数。
4.根据权利要求2或3所述的先进绝热压缩空气储能系统全动态模型的生成方法,其特征在于,所述换热系统宽工况动态模型为:
<mrow>
<msub>
<mi>&eta;</mi>
<mi>t</mi>
</msub>
<mo>=</mo>
<mo>&lsqb;</mo>
<msub>
<mi>n</mi>
<mi>t</mi>
</msub>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>T</mi>
<mrow>
<mi>t</mi>
<mo>,</mo>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>0.5</mn>
</mrow>
</msup>
<mo>&rsqb;</mo>
<mo>/</mo>
<mo>&lsqb;</mo>
<msubsup>
<mi>n</mi>
<mi>t</mi>
<mn>0</mn>
</msubsup>
<msup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>T</mi>
<mrow>
<mi>t</mi>
<mo>,</mo>
<mi>i</mi>
<mi>n</mi>
</mrow>
<mn>0</mn>
</msubsup>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>0.5</mn>
</mrow>
</msup>
<mo>&rsqb;</mo>
<mo>,</mo>
</mrow>
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&Phi;</mi>
<mrow>
<mi>H</mi>
<mi>E</mi>
<mi>X</mi>
</mrow>
</msub>
<mo>=</mo>
<mi>&epsiv;</mi>
<mo>&CenterDot;</mo>
<msub>
<mi>C</mi>
<mi>min</mi>
</msub>
<mrow>
<mo>(</mo>
<msubsup>
<mi>T</mi>
<mrow>
<mi>i</mi>
<mi>n</mi>
</mrow>
<mi>h</mi>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>T</mi>
<mrow>
<mi>i</mi>
<mi>n</mi>
</mrow>
<mi>c</mi>
</msubsup>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>&epsiv;</mi>
<mo>=</mo>
<mfrac>
<mrow>
<mn>1</mn>
<mo>-</mo>
<mi>exp</mi>
<mo>&lsqb;</mo>
<mo>-</mo>
<mi>N</mi>
<mi>T</mi>
<mi>U</mi>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mi>&chi;</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<mn>1</mn>
<mo>-</mo>
<mi>&chi;</mi>
<mi>exp</mi>
<mo>&lsqb;</mo>
<mo>-</mo>
<mi>N</mi>
<mi>T</mi>
<mi>U</mi>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mi>&chi;</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>,</mo>
</mrow>
其中,ΦHEX为换热器实际换热量,Cmin为最小热容,代表换热器两侧冷热流体入口温度,ε为由换热器换热系数、换热面积等界定的系数,受运行工况影响。
5.根据权利要求1所述的先进绝热压缩空气储能系统全动态模型的生成方法,其特征在于,所述动态方程为:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<msub>
<mi>dT</mi>
<mi>r</mi>
</msub>
</mrow>
<mrow>
<mi>d</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>m</mi>
<mi>r</mi>
</msub>
</mfrac>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mi>k</mi>
</mfrac>
<mo>(</mo>
<msub>
<mi>m</mi>
<mrow>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
<msub>
<mi>T</mi>
<mrow>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>m</mi>
<mrow>
<mi>o</mi>
<mi>u</mi>
<mi>t</mi>
</mrow>
</msub>
<msub>
<mi>T</mi>
<mi>r</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mfrac>
<mrow>
<msub>
<mi>h</mi>
<mi>w</mi>
</msub>
<msub>
<mi>A</mi>
<mi>w</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>T</mi>
<mi>w</mi>
</msub>
<mo>-</mo>
<msub>
<mi>T</mi>
<mi>r</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<msub>
<mi>c</mi>
<mrow>
<mi>p</mi>
<mo>,</mo>
<mi>a</mi>
</mrow>
</msub>
</mfrac>
<mo>)</mo>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<msub>
<mi>dm</mi>
<mi>r</mi>
</msub>
</mrow>
<mrow>
<mi>d</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>=</mo>
<msub>
<mi>m</mi>
<mrow>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>m</mi>
<mrow>
<mi>o</mi>
<mi>u</mi>
<mi>t</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>,</mo>
</mrow>
其中,Tw,hw,Aw分别为环境温度、换热系数及接触面积,mr、Tr为储气库高压空气质量及温度,cp,a为高压空气定压比热容,k为比热系数。
6.一种先进绝热压缩空气储能系统全动态模型的生成装置,其特征在于,包括:
第一构建模块,用于根据压缩机和膨胀机的宽工况运行特性,构建空气压缩热能和压力势能的解耦生产单元和耦合释能单元的压缩机和膨胀机宽工况动态模型;
第二构建模块,用于根据压缩侧换热系统和膨胀机侧换热系统的宽工况运行特性,构建压缩热能传输收集单元和压缩热能传输释放单元的换热系统宽工况动态模型;
第三构建模块,用于根据储热罐温度流量和储气罐压力温度的动态方程,构建压缩热能和压力势能存储单元的储能系统宽工况动态模型;以及
生成模块,用于基于电动机及发电机动态模型,建立先进绝热压缩空气储能与电网耦合接口的联系,并且构建从系统组件动态特性到系统整体性能间的关联关系,以为智能电网不同场景应用及能效提升提供模型支撑。
7.根据权利要求6所述的先进绝热压缩空气储能系统全动态模型的生成装置,其特征在于,所述压缩机的透平宽工况运行特性描述为:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&beta;</mi>
<mi>c</mi>
</msub>
<mo>=</mo>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>G</mi>
<mi>c</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>a</mi>
<mn>2</mn>
</msub>
<msub>
<mi>G</mi>
<mi>c</mi>
</msub>
<mo>+</mo>
<msub>
<mi>a</mi>
<mn>3</mn>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&eta;</mi>
<mi>c</mi>
</msub>
<mo>=</mo>
<mo>&lsqb;</mo>
<mn>1</mn>
<mo>-</mo>
<mi>c</mi>
<msup>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>n</mi>
<mi>c</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>&rsqb;</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>n</mi>
<mi>c</mi>
</msub>
<mo>/</mo>
<msub>
<mi>G</mi>
<mi>c</mi>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>-</mo>
<mo>(</mo>
<mrow>
<msub>
<mi>n</mi>
<mi>c</mi>
</msub>
<mo>/</mo>
<msub>
<mi>G</mi>
<mi>c</mi>
</msub>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>,</mo>
</mrow>
其中,βc,ηc为压缩机的压缩比与等熵效率,Gc,nc分别表示宽工况运行时的压缩机质量流量与工质流速,c相应系数,a1,a2,a3为常数;
并且,类似压缩机的透平宽工况运行特性为:
<mrow>
<msub>
<mi>&eta;</mi>
<mi>t</mi>
</msub>
<mo>/</mo>
<msubsup>
<mi>&eta;</mi>
<mi>t</mi>
<mn>0</mn>
</msubsup>
<mo>=</mo>
<mo>&lsqb;</mo>
<mn>1</mn>
<mo>-</mo>
<mi>t</mi>
<msup>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>n</mi>
<mi>t</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>&rsqb;</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>n</mi>
<mi>t</mi>
</msub>
<mo>/</mo>
<msub>
<mi>G</mi>
<mi>t</mi>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>-</mo>
<mo>(</mo>
<mrow>
<msub>
<mi>n</mi>
<mi>t</mi>
</msub>
<mo>/</mo>
<msub>
<mi>G</mi>
<mi>t</mi>
</msub>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
其中,ηt为透平等熵效率,Gt,nt分别表示宽工况运行时的透平质量流量与工质流速。
8.根据权利要求7所述的先进绝热压缩空气储能系统全动态模型的生成装置,其特征在于,其中,
<mrow>
<msub>
<mi>G</mi>
<mi>c</mi>
</msub>
<mo>=</mo>
<mo>&lsqb;</mo>
<msub>
<mi>m</mi>
<mi>c</mi>
</msub>
<mfrac>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>T</mi>
<mrow>
<mi>c</mi>
<mo>,</mo>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mn>0.5</mn>
</msup>
<msub>
<mi>P</mi>
<mrow>
<mi>c</mi>
<mo>,</mo>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
</mfrac>
<mo>&rsqb;</mo>
<mo>/</mo>
<mo>&lsqb;</mo>
<msubsup>
<mi>m</mi>
<mi>c</mi>
<mn>0</mn>
</msubsup>
<mfrac>
<msup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>T</mi>
<mrow>
<mi>c</mi>
<mo>,</mo>
<mi>i</mi>
<mi>n</mi>
</mrow>
<mn>0</mn>
</msubsup>
<mo>)</mo>
</mrow>
<mn>0.5</mn>
</msup>
<msubsup>
<mi>P</mi>
<mrow>
<mi>c</mi>
<mo>,</mo>
<mi>i</mi>
<mi>n</mi>
</mrow>
<mn>0</mn>
</msubsup>
</mfrac>
<mo>&rsqb;</mo>
<mo>,</mo>
</mrow>
<mrow>
<msub>
<mi>&eta;</mi>
<mi>c</mi>
</msub>
<mo>=</mo>
<mo>&lsqb;</mo>
<msub>
<mi>n</mi>
<mi>c</mi>
</msub>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>T</mi>
<mrow>
<mi>c</mi>
<mo>,</mo>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>0.5</mn>
</mrow>
</msup>
<mo>&rsqb;</mo>
<mo>/</mo>
<mo>&lsqb;</mo>
<msubsup>
<mi>n</mi>
<mi>c</mi>
<mn>0</mn>
</msubsup>
<msup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>T</mi>
<mrow>
<mi>c</mi>
<mo>,</mo>
<mi>i</mi>
<mi>n</mi>
</mrow>
<mn>0</mn>
</msubsup>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>0.5</mn>
</mrow>
</msup>
<mo>&rsqb;</mo>
<mo>,</mo>
</mrow>
<mrow>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mfrac>
<msub>
<mi>n</mi>
<mi>c</mi>
</msub>
<msub>
<mi>a</mi>
<mn>0</mn>
</msub>
</mfrac>
<mo>,</mo>
<msub>
<mi>a</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>b</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<mn>2</mn>
<msub>
<mi>b</mi>
<mn>2</mn>
</msub>
<msubsup>
<mi>n</mi>
<mi>c</mi>
<mn>2</mn>
</msubsup>
<mo>)</mo>
</mrow>
<msub>
<mi>a</mi>
<mn>0</mn>
</msub>
</mfrac>
<mo>,</mo>
<msub>
<mi>a</mi>
<mn>3</mn>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msup>
<msub>
<mi>b</mi>
<mn>2</mn>
</msub>
<mn>2</mn>
</msup>
<msubsup>
<mi>n</mi>
<mi>c</mi>
<mn>3</mn>
</msubsup>
<mo>-</mo>
<msub>
<mi>b</mi>
<mn>1</mn>
</msub>
<msub>
<mi>b</mi>
<mn>2</mn>
</msub>
<msub>
<mi>n</mi>
<mi>c</mi>
</msub>
</mrow>
<msub>
<mi>a</mi>
<mn>0</mn>
</msub>
</mfrac>
<mo>,</mo>
</mrow>
<mrow>
<msub>
<mi>a</mi>
<mn>0</mn>
</msub>
<mo>=</mo>
<mo>&lsqb;</mo>
<msub>
<mi>b</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mfrac>
<msub>
<mi>b</mi>
<mn>2</mn>
</msub>
<msub>
<mi>n</mi>
<mi>c</mi>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>n</mi>
<mi>c</mi>
</msub>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>n</mi>
<mi>c</mi>
</msub>
<mo>-</mo>
<msub>
<mi>b</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>&rsqb;</mo>
<mo>,</mo>
</mrow>
并且,
<mrow>
<msub>
<mi>G</mi>
<mi>t</mi>
</msub>
<mo>=</mo>
<mo>&lsqb;</mo>
<msub>
<mi>m</mi>
<mi>t</mi>
</msub>
<mfrac>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>T</mi>
<mrow>
<mi>t</mi>
<mo>,</mo>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mn>0.5</mn>
</msup>
<msub>
<mi>P</mi>
<mrow>
<mi>t</mi>
<mo>,</mo>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
</mfrac>
<mo>&rsqb;</mo>
<mo>/</mo>
<mo>&lsqb;</mo>
<msubsup>
<mi>m</mi>
<mi>t</mi>
<mn>0</mn>
</msubsup>
<mfrac>
<msup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>T</mi>
<mrow>
<mi>t</mi>
<mo>,</mo>
<mi>i</mi>
<mi>n</mi>
</mrow>
<mn>0</mn>
</msubsup>
<mo>)</mo>
</mrow>
<mn>0.5</mn>
</msup>
<msubsup>
<mi>P</mi>
<mrow>
<mi>t</mi>
<mo>,</mo>
<mi>i</mi>
<mi>n</mi>
</mrow>
<mn>0</mn>
</msubsup>
</mfrac>
<mo>&rsqb;</mo>
<mo>,</mo>
</mrow>
其中,Tc,in为实际运行工况对应的压缩机进口温度,Pc,in为实际运行工况对应的压缩机进口压力、mc为实际运行工况对应的压缩机进口质量流量,为额定设计工况下的压缩机进口温度、为额定设计工况下的压缩机进口压力,为额定设计工况下的压缩机进口质量流量,Tt,in为实际运行工况对应的透平进口温度,Pt,in为实际运行工况对应的透平进口压力,mt为实际运行工况对应的透平进口质量流量,额定设计工况下的为额定设计工况下的透平进口温度、为额定设计工况下的透平进口压力,为额定设计工况下的透平进口质量流量,b1,b2分别为相应系数。
9.根据权利要求7或8所述的先进绝热压缩空气储能系统全动态模型的生成装置,其特征在于,所述换热系统宽工况动态模型为:
<mrow>
<msub>
<mi>&eta;</mi>
<mi>t</mi>
</msub>
<mo>=</mo>
<mo>&lsqb;</mo>
<msub>
<mi>n</mi>
<mi>t</mi>
</msub>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>T</mi>
<mrow>
<mi>t</mi>
<mo>,</mo>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>0.5</mn>
</mrow>
</msup>
<mo>&rsqb;</mo>
<mo>/</mo>
<mo>&lsqb;</mo>
<msubsup>
<mi>n</mi>
<mi>t</mi>
<mn>0</mn>
</msubsup>
<msup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>T</mi>
<mrow>
<mi>t</mi>
<mo>,</mo>
<mi>i</mi>
<mi>n</mi>
</mrow>
<mn>0</mn>
</msubsup>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>0.5</mn>
</mrow>
</msup>
<mo>&rsqb;</mo>
<mo>,</mo>
</mrow>
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&Phi;</mi>
<mrow>
<mi>H</mi>
<mi>E</mi>
<mi>X</mi>
</mrow>
</msub>
<mo>=</mo>
<mi>&epsiv;</mi>
<mo>&CenterDot;</mo>
<msub>
<mi>C</mi>
<mi>min</mi>
</msub>
<mrow>
<mo>(</mo>
<msubsup>
<mi>T</mi>
<mrow>
<mi>i</mi>
<mi>n</mi>
</mrow>
<mi>h</mi>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>T</mi>
<mrow>
<mi>i</mi>
<mi>n</mi>
</mrow>
<mi>c</mi>
</msubsup>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>&epsiv;</mi>
<mo>=</mo>
<mfrac>
<mrow>
<mn>1</mn>
<mo>-</mo>
<mi>exp</mi>
<mo>&lsqb;</mo>
<mo>-</mo>
<mi>N</mi>
<mi>T</mi>
<mi>U</mi>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mi>&chi;</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<mn>1</mn>
<mo>-</mo>
<mi>&chi;</mi>
<mi>exp</mi>
<mo>&lsqb;</mo>
<mo>-</mo>
<mi>N</mi>
<mi>T</mi>
<mi>U</mi>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mi>&chi;</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>,</mo>
</mrow>
其中,ΦHEX为换热器实际换热量,Cmin为最小热容,代表换热器两侧冷热流体入口温度,ε为由换热器换热系数、换热面积等界定的系数,受运行工况影响。
10.根据权利要求6所述的先进绝热压缩空气储能系统全动态模型的生成装置,其特征在于,所述动态方程为:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<msub>
<mi>dT</mi>
<mi>r</mi>
</msub>
</mrow>
<mrow>
<mi>d</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>m</mi>
<mi>r</mi>
</msub>
</mfrac>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mi>k</mi>
</mfrac>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<msub>
<mi>m</mi>
<mrow>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
<msub>
<mi>T</mi>
<mrow>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>m</mi>
<mrow>
<mi>o</mi>
<mi>u</mi>
<mi>t</mi>
</mrow>
</msub>
<msub>
<mi>T</mi>
<mi>r</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mfrac>
<mrow>
<msub>
<mi>h</mi>
<mi>w</mi>
</msub>
<msub>
<mi>A</mi>
<mi>w</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>T</mi>
<mi>w</mi>
</msub>
<mo>-</mo>
<msub>
<mi>T</mi>
<mi>r</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<msub>
<mi>c</mi>
<mrow>
<mi>p</mi>
<mo>,</mo>
<mi>a</mi>
</mrow>
</msub>
</mfrac>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<msub>
<mi>dm</mi>
<mi>r</mi>
</msub>
</mrow>
<mrow>
<mi>d</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>=</mo>
<msub>
<mi>m</mi>
<mrow>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>m</mi>
<mrow>
<mi>o</mi>
<mi>u</mi>
<mi>t</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>,</mo>
</mrow>
其中,Tw,hw,Aw分别为环境温度、换热系数及接触面积,mr、Tr为储气库高压空气质量及温度,cp,a为高压空气定压比热容,k为比热系数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711318890.0A CN108062440A (zh) | 2017-12-12 | 2017-12-12 | 先进绝热压缩空气储能系统全动态模型的生成方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711318890.0A CN108062440A (zh) | 2017-12-12 | 2017-12-12 | 先进绝热压缩空气储能系统全动态模型的生成方法及装置 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN108062440A true CN108062440A (zh) | 2018-05-22 |
Family
ID=62138366
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711318890.0A Pending CN108062440A (zh) | 2017-12-12 | 2017-12-12 | 先进绝热压缩空气储能系统全动态模型的生成方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108062440A (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109405349A (zh) * | 2018-07-30 | 2019-03-01 | 东南大学 | 一种使用纯工质的高温热泵系统的控制方法 |
CN109583012A (zh) * | 2018-10-18 | 2019-04-05 | 国网安徽省电力有限公司 | 先进绝热压缩空气储能与风电协同运行调度方法与装置 |
CN109611320A (zh) * | 2018-10-17 | 2019-04-12 | 国网安徽省电力有限公司 | 先进绝热压缩空气储能系统的风储协调控制方法及装置 |
CN110489703A (zh) * | 2019-07-05 | 2019-11-22 | 清华大学 | 先进绝热压缩空气储能热电联供运行可行域分析方法 |
CN111946409A (zh) * | 2019-05-17 | 2020-11-17 | 中国中原对外工程有限公司 | 匹配钠堆的超临界co2布雷顿循环设计方法 |
CN113067353A (zh) * | 2021-03-05 | 2021-07-02 | 华中科技大学 | 一种考虑csp电站耦合aa-caes电站的系统优化调度方法 |
CN115223436A (zh) * | 2022-07-07 | 2022-10-21 | 西安热工研究院有限公司 | 一种模拟重力压缩空气储能系统运行的实验模型装置 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2447501A2 (en) * | 2010-10-29 | 2012-05-02 | Nuovo Pignone S.p.A. | Systems and methods for pre-heating compressed air in advanced adiabatic compressed air energy storage systems |
WO2016166095A1 (fr) * | 2015-04-13 | 2016-10-20 | IFP Energies Nouvelles | Systeme et procede de stockage et de recuperation d'energie par air comprime avec chauffage a volume constant |
CN106194587A (zh) * | 2016-07-27 | 2016-12-07 | 山东大学 | 一种应用于微型风电系统的压缩空气混合储能系统及方法 |
CN106246269A (zh) * | 2016-07-27 | 2016-12-21 | 山东大学 | 一种可重构压缩空气储能系统及其优化控制方法 |
-
2017
- 2017-12-12 CN CN201711318890.0A patent/CN108062440A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2447501A2 (en) * | 2010-10-29 | 2012-05-02 | Nuovo Pignone S.p.A. | Systems and methods for pre-heating compressed air in advanced adiabatic compressed air energy storage systems |
WO2016166095A1 (fr) * | 2015-04-13 | 2016-10-20 | IFP Energies Nouvelles | Systeme et procede de stockage et de recuperation d'energie par air comprime avec chauffage a volume constant |
CN106194587A (zh) * | 2016-07-27 | 2016-12-07 | 山东大学 | 一种应用于微型风电系统的压缩空气混合储能系统及方法 |
CN106246269A (zh) * | 2016-07-27 | 2016-12-21 | 山东大学 | 一种可重构压缩空气储能系统及其优化控制方法 |
Non-Patent Citations (1)
Title |
---|
ADRIANO SCIACOVELLI 等: "Dynamic simulation of Adiabatic Compressed Air Energy Storage(A-CAES) plant with integrated thermal storage–Link between components performance and plant performance", 《APPLIED ENERGY》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109405349A (zh) * | 2018-07-30 | 2019-03-01 | 东南大学 | 一种使用纯工质的高温热泵系统的控制方法 |
CN109611320A (zh) * | 2018-10-17 | 2019-04-12 | 国网安徽省电力有限公司 | 先进绝热压缩空气储能系统的风储协调控制方法及装置 |
CN109611320B (zh) * | 2018-10-17 | 2020-04-14 | 国网安徽省电力有限公司 | 绝热压缩空气储能系统的风储协调控制方法及装置 |
CN109583012A (zh) * | 2018-10-18 | 2019-04-05 | 国网安徽省电力有限公司 | 先进绝热压缩空气储能与风电协同运行调度方法与装置 |
CN109583012B (zh) * | 2018-10-18 | 2023-03-24 | 国网安徽省电力有限公司 | 先进绝热压缩空气储能与风电协同运行调度方法与装置 |
CN111946409A (zh) * | 2019-05-17 | 2020-11-17 | 中国中原对外工程有限公司 | 匹配钠堆的超临界co2布雷顿循环设计方法 |
CN110489703A (zh) * | 2019-07-05 | 2019-11-22 | 清华大学 | 先进绝热压缩空气储能热电联供运行可行域分析方法 |
CN113067353A (zh) * | 2021-03-05 | 2021-07-02 | 华中科技大学 | 一种考虑csp电站耦合aa-caes电站的系统优化调度方法 |
CN115223436A (zh) * | 2022-07-07 | 2022-10-21 | 西安热工研究院有限公司 | 一种模拟重力压缩空气储能系统运行的实验模型装置 |
CN115223436B (zh) * | 2022-07-07 | 2023-11-17 | 西安热工研究院有限公司 | 一种模拟重力压缩空气储能系统运行的实验模型装置 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108062440A (zh) | 先进绝热压缩空气储能系统全动态模型的生成方法及装置 | |
EP3575576A1 (en) | Compressed-air-storing power generation method and compressed-air-storing power generation device | |
CN109026241A (zh) | 一种热泵压缩空气储能系统 | |
CN103814199A (zh) | 使用涡轮机的压缩气体能量储存系统 | |
CN105043147A (zh) | 一种采用液态蓄冷工质的液化压缩空气储能系统 | |
WO2013037658A1 (en) | Thermal energy storage and recovery arrangement | |
CN108331618B (zh) | 一种10mw压缩空气储能方法和系统 | |
Li et al. | Exploring thermodynamic potential of multiple phase change thermal energy storage for adiabatic compressed air energy storage system | |
CN104775856A (zh) | 一种联合循环发电机组汽轮机快速启动控制方法和装置 | |
Garvey et al. | Analysis of a wind turbine power transmission system with intrinsic energy storage capability | |
Zhang et al. | Modelling and experimental validation of advanced adiabatic compressed air energy storage with off‐design heat exchanger | |
CN102352777B (zh) | 可逆单螺杆压缩膨胀机储能发电系统和储能发电方法 | |
CN106050344A (zh) | 一种利用低温气态蓄冷工质的深冷液化空气储能系统 | |
CN105507971A (zh) | 太阳能热动力储能系统 | |
Vongmanee | The renewable energy applications for uninterruptible power supply based on compressed air energy storage system | |
CN106410855A (zh) | 一种独立微能源网及其优化配置方法 | |
CN205503202U (zh) | 太阳能热动力储能系统 | |
CN202560330U (zh) | 利用两级螺杆膨胀机蒸汽余压的动力系统 | |
CN110112781A (zh) | 一种综合能源电力系统分布互联建模分析方法 | |
Chen et al. | Performance study of salt cavern air storage based non-supplementary fired compressed air energy storage system | |
CN105890285B (zh) | 一种利用混合工质的深冷液化空气储能方法和系统 | |
CN204457895U (zh) | 用太阳能光热/工业余热发电的汽轮发电机组 | |
CN108869213A (zh) | 光子增强热离子发射与二氧化碳循环联合发电装置及方法 | |
CN114626587A (zh) | 一种考虑天然气管道掺氢的气电联合系统运行优化方法 | |
CN112065518A (zh) | 一种用于燃煤机组液态压缩空气储能的汽电双驱发电系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20180522 |
|
RJ01 | Rejection of invention patent application after publication |