CN104408534B - 基于联立方程模型的塔式太阳能热电系统运行优化方法 - Google Patents
基于联立方程模型的塔式太阳能热电系统运行优化方法 Download PDFInfo
- Publication number
- CN104408534B CN104408534B CN201410734256.5A CN201410734256A CN104408534B CN 104408534 B CN104408534 B CN 104408534B CN 201410734256 A CN201410734256 A CN 201410734256A CN 104408534 B CN104408534 B CN 104408534B
- Authority
- CN
- China
- Prior art keywords
- heat
- model
- power station
- working medium
- storage tank
- 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.)
- Active
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
- G06Q50/06—Electricity, gas or water supply
-
- 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
- Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
- Y02E40/70—Smart grids as climate change mitigation technology in the energy generation sector
-
- 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
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
Abstract
本发明公开了一种基于联立方程模型的塔式太阳能热电系统运行优化方法。包括：对热电系统各子系统进行分析，建立各子系统的动态数学模型；采用联立方程法建立塔式太阳能热电系统整体模型；以一天内电站收益最高为目标，热电系统参数的上下限为约束条件，构造优化问题；利用控制向量参数化方法来求解优化问题，得到最优的热电系统运行调度策略。基于联立方程模型的塔式太阳能热电系统运行优化方法，能够保证电站在市场电价下获得最高的收益，并且发电功率、热能存储、工质流量以及工质温度都能满足电站正常安全运行的要求，有利于保护电站设备、提高系统稳定性，为塔式太阳能热电站的运行提供参考依据。
Description
技术领域
本发明涉及塔式太阳能热电系统领域，特别是涉及一种基于联立方程模型的塔式太阳能热电系统运行优化方法。
背景技术
太阳能是一种可再生的新能源，储量丰富，并且具有清洁、环保、持续、长久的优势，越来越受到人们的广泛关注。美国能源部研究表明，目前太阳能聚光类发电中，塔式发电系统聚光倍数高，可得到较高的光热转换效率，更可能取代常规能源、实现大功率发电。太阳能热发电的主要原理是：定日镜场反射聚焦太阳光至集热器，加热集热器中的工作介质，工作介质进入蒸汽发生器换热，加热给水使其产生高温蒸汽推动汽轮机的运作，从而发电。
太阳能间歇性这一特性的存在，要求太阳能热电站加入储热系统，从而使得发电产出变得可控。这就要求对太阳能电站的操作运营包括发电量等进行优化调度。加入储热系统的意义在于，通过镜场集热器收集得到的太阳能，不再只能输送向发电装置部分，取而代之的，可以使其全部或者部分储进储热罐中，并在合适的时候再次取出来进行发电。
塔式太阳能热电系统一般包括集热器、储热罐、管道和蒸汽发生器。正常运行时一般要保证系统有稳定的功率输出，但经济效益上不是最优的。如果考虑电价的预测，通过储热系统能在高价格时间段释放更多热量，太阳能电站能够收到更大的收益。一种优秀的优化调度方法能够应对天气状况、环境温度以及电价水平等方面的变化，以最佳的策略，尽可能快的对电站的操作运营作出诸如调整发电量等合理的调度，从而使得电站取得更大的收益。
现有的针对塔式太阳能热电系统的优化调度策略，有的是针对一年的长时间段的根据不同季节进行调整的太阳能电站优化操作策略，有的是根据自由市场电价和子系统的输入输出关系进行一天内子系统间的静态能量调度，也有调度策略使得电站产出量尽可能大但未与市场电价进行关联。这些优化方法并不是基于电站整体机理模型，也没有进行根据整体模型进行实时动态调度，因此没有考虑运行中各参数的变化情况，无法对系统内存在的各种影响系统稳定安全运行的参数变化进行实施调控。
发明内容
本发明提供了一种基于联立方程模型的塔式太阳能热电系统运行优化方法，使电站在市场电价下获得最高的收益，并且发电功率、热能存储、工质流量以及工质温度都能满足电站正常安全运行的要求，达到保护电站设备、提高系统稳定性的目的。
本发明采用的技术方案如下：
(1)对热电系统各子系统进行分析，建立各子系统的动态数学模型。
塔式太阳能电站的热电系统包括集热器、储热系统、蒸汽发生器和管路系统。在实际的生产过程中，由于太阳辐射强度DNI和天气状况因素的变化，系统内的各变量也会随着时间变化，因此需要对各个模块建立动态数学模型，也就是系统在动态过程中各种参数随着时间变化的数学描述。
建立动态模型可以采用机理建模方法，根据各子系统的工作原理和换热方式，用质量平衡、能量平衡和流体传热平衡方程，以及传热工质的物性参数方程进行描述，分别建立接收器、管道、储热罐和蒸汽发生器子系统的动态数学模型。
接收器模型表示为：
其中，Amirror为定日镜总面积，I为太阳辐射强度，ξ为定日镜场的效率；ε为黑度，σ为黑体辐射常数；Am为吸热器外部面积，Dm为吸热器外部直径，Tf为吸热器内工质温度，Ta为环境温度，Af为吸热器内部面积，m1为流过吸热器内的工质流量，hm、hf为体积对流传热系数。
管道模型表示为：
其中，τ1是与流体密度、比热、管径相关的常数，τ12、τ2是与管壁金属密度、比热、管径相关的常数，Ta为管道外环境温度，Tf为管内流体温度，Tw为管壁温度。
储热罐模型表示为：
其中，A为填充床的横截面积，m为流过储热罐的空气流量，Tf为储热罐内的工质温度，U为流体与环境之间的传热系数，Ta为环境温度，D为储热罐的直径。
蒸汽发生器模型表示为：
式中，ρ为流体密度，D为单位面积的质量流速，h为流体比焓，Q为单位长度的热流率，A为流道的横截面积，P为流体压力。
(2)根据各子系统的动态数学模型，采用联立方程法建立塔式太阳能热电系统整体模型。
对整体模型机理建模常用的方法有序贯模块法和联立方程法。就求解速度而言，特别是对于大规模流程工业，联立方程法所占的优势相对突出。塔式太阳能热电系统具有模型复杂且方程、变量数多的特点，具体表现在建模时涉及到大量的系统模型方程组，包括物性计算方程、能量平衡方程、流体换热方程，并且其中有大量的微分方程。故而，传统的序贯模块法求解这种变量多、稀疏性强、高度非线性的流程模拟方程时速度较慢；联立方程法则有较快的求解速度，更适合对塔式太阳能热电系统整体模型的求解。
模拟的平台选用计算功能强大的过程模拟软件gPROMS。根据建立的接收器、管道、储热罐和蒸汽发生器等主要子系统的动态数学模型，在gPROMS平台下建立起热电系统的各个子模型，然后按照输入输出关系通过模型连接方程进行联立。考虑系统在不同工况下的运行情况，按照不同运行模式时的流程结构分别进行连接，接着将DNI和传热工质流量作为模型输入，将蒸汽发生器功率和接收器出口传热工质温度作为模型输出，就可以得到塔式太阳能热电系统整体模型。
(3)考虑实时电价，以一天内电站收益最高为目标，热电系统参数的上下限为约束条件，构造优化问题。
正常运行时要保证系统有稳定的功率输出，但经济效益上不是最优的。考虑电价的预测，储热系统能把低需求时间段的热量转换到高需求时间段的热量，这对于储存的热量来说，这是一个升值的过程，使得太阳能电站能够收到更大的收益。在给定当天每小时DNI数据以及市场电价的情况下，计划安排太阳能热电站每小时的发电量，使得当天电站所得收益最大化。
研究塔式太阳能热电系统正常运行的优化策略，满足热电系统安全运行要求的同时尽可能使得能量输出达到最大值，优化目标为一天内电站的收益最高。假定市场电价每小时变动一次。因此第i小时的电站收益为该时段内的市场电价乘以该时段内的发电量：
η(i)＝Price(i)·Qi
其中，η(i)为电站第i小时的收益，Price(i)为第i小时的市场电价，Qi是电站第i小时的发电量。
当天总收益为各个小时收益的累加。目标函数可以表示为：
电站运行时需要满足整个系统的安全运行要求，需要考虑以下约束条件：
1)储热罐有固定的容量，当储热罐中储存的热量达到上限则无法继续储存热量，相反，当储热罐中储存的热量过少时也不能提取热量；2)发电装置有其最大发电功率，当发电装置达到满功率运作时，即使蒸汽发生器提供更多的高温蒸汽去驱动汽轮机，也无法汽轮机达到更大的发电功率；3)阀门开度有一定的范围，同时工质流量需要保证系统正常运行，因此各管路中的工质流量有上下限；4)考虑传热工质在不同温度下的稳定性、腐蚀性等特性，以及接收器的使用寿命等因素，传热工质的温度也有上下限。约束条件可以表示为：
储热罐储能约束:Qminch≤Qich≤Qmaxch
发电功率约束:0≤Pi≤Pmax
传热工质流量约束:mjmin≤mj≤mjmax
传热工质温度约束：Tfmin≤Tjfa≤Tfmax
其中，Qich为第i小时储热罐储存的能量，Qmaxch、Qminch为储热罐最大、最小储热量；Pi为第i小时的发电功率，Pmax为发电装置的最大发电功率；mj为第j根管道内传热工质的流量，mjmax、mjmin为第j根管道内传热工质的流量的上下限；Tjfa为第j根管道内传热工质的温度，Tfmax、Tfmin为传热工质温度的上下限。
因此，整个优化问题可以描述为：
(4)利用控制向量参数化方法来求解优化问题，得到最优的热电系统运行调度策略。
基于联立方程的塔式太阳能热电系统整体模型是复杂的非线性模型，存在各种约束条件，一般的优化方法难以解决这类模型的动态优化问题。控制向量参数化方法是一种解决过程动态优化问题的有效方法.该方法首先将控制变量离散化为有限时间段；然后在优化时域内对系统模型积分,将优化问题中的系统模型约束完全消去，减少了非线性规划中的变量个数；最后用NLP方法求解，得到优化控制参数。
CVP_SS方法基于控制向量参数化方法，并采用了打靶策略，该方法的思想是将两点边界值问题转化为初值问题进行计算，通过优化算法确定每个控制区间的控制变量值以及持续时间。对于NLP非线性规划问题，借助gPROMS提供的基于SQP算法的SRQPD求解器和混合整数规划OAERAP求解器进行求解。
CVP_SS算法求解动态优化问题的实施流程如下：
步骤一：设置控制变量、控制间隔的初始值，确定时间节点数、约束变量上下限和优化时域；
步骤二：通过控制参数将控制变量离散化，将优化时域[0,tf]分割为N个时间单元，在各时间单元内用分段常数函数逼近连续控制量,即u(t)＝ui，其中ti-1≤t<ti,i＝1,2……N；
步骤三：在整个优化时域内确定各变量的动态过程；
步骤四：计算优化目标函数和约束变量，判断目标函数是否最优、约束条件是否满足。此时如果满足终止条件则算法终止转到步骤六，否则转到步骤五；
步骤五：计算目标函数关于控制参数的梯度，调整控制变量和控制间隔，转到步骤二；
步骤六：输出控制变量、约束变量的动态曲线以及优化目标值，得到最优的运行调度策略。
当算法终止时，得到最优的热电系统运行策略，确定了系统各主要变量在运行时间内的动态变化过程。同时算法结果也给出了主要变量的变化曲线、最高收益值等信息。
本发明的有益效果是：按照本发明提出基于联立方程模型的塔式太阳能热电系统运行优化方法，得到的最优运行优化策略，能够保证电站在市场电价下获得最高的收益，并且发电功率、热能存储、工质流量以及工质温度都能满足电站正常安全运行的要求，有利于保护电站设备、提高系统稳定性，为塔式太阳能热电站的运行提供参考依据。
附图说明
图1是基于联立方程模型的塔式太阳能热电系统运行优化方法的流程图；
图2是实施例中的塔式太阳能电站原理图；
图3是实施例中某日8小时的实时电价；
图4是实施例中优化后的蒸汽发生器功率曲线图；
图5是实施例中优化后的储热罐中储存能量曲线图。
具体实施方式
如图1所示，一种基于联立方程模型的塔式太阳能热电系统运行优化方法，实施步骤如下：
(1)对热电系统各子系统进行分析，建立各子系统的动态数学模型。
塔式太阳能电站的热电系统包括集热器、储热系统、蒸汽发生器和管路系统。在实际的生产过程中，由于太阳辐射强度DNI和天气状况因素的变化，系统内的各个主要变量也会随着时间变化，因此需要对各个模块建立动态数学模型，也就是系统在动态过程中各种参数随着时间变化的数学描述。
建立动态模型采用机理建模方法，根据各子系统的工作原理和换热方式，用质量平衡、能量平衡和流体传热平衡方程，以及传热工质的物性参数方程进行描述，分别建立接收器、管道、储热罐和蒸汽发生器等主要子系统的动态数学模型。
在建模时做以下假设：接收器表面接收到的能流密度相同，只需要考虑DNI作为能量输入而不用考虑能流密度分布；传热工质为非压缩流体，可以认为每个子模型的入口工质流量等于出口工质流量。
接收器模型表示为：
其中，Amirror为定日镜总面积，I为太阳辐射强度，ξ为定日镜场的效率；ε为黑度；σ为黑体辐射常数；Am为吸热器外部面积，Dm为吸热器外部直径，Tf为吸热器内工质温度，Ta为环境温度，Af为吸热器内部面积，m1为流过吸热器内的工质流量，hm、hf为体积对流传热系数。
管道模型表示为：
其中，τ1是与流体密度、比热、管径相关的常数，τ12、τ2是与管壁金属密度、比热、管径相关的常数，Ta为管道外环境温度，Tf为管内流体温度，Tw为管壁温度。
储热罐模型表示为：
其中，A为填充床的横截面积，m为流过储热罐的空气流量，Tf为储热罐内的工质温度，U为流体与环境之间的传热系数，Ta为环境温度，D为储热罐的直径。
蒸汽发生器模型表示为：
式中，ρ为流体密度，D为单位面积的质量流速，h为流体比焓，Q为单位长度的热流率，A为流道的横截面积，P为流体压力。
(2)根据各子系统的动态数学模型，采用联立方程法建立塔式太阳能热电系统整体模型。
对整体模型机理建模常用的方法有序贯模块法和联立方程法。在序贯法中，系统可以看作由一个个子模块组成，由于各个子模块都有其自身特定的解法以及收敛准则，求解算法与模块不可分离。序贯模块法与实际过程联系紧密，直观性强，但是由于模块中信息流动方向固定，且拓扑结构具有多层次特点，因此在处理再循环流股和优化计算时都需要在多个不同层次上循环迭代，耗时较大，计算效率不高。与序贯法相对的是联立方程法，联立方程法建模时，打破以往的模块化结构，将所有方程联立起来，将物性计算模块统一在模型内部，这样模拟与优化计算过程中就不需要多层迭代计算，从而可以大大提高求解效率。由于其收敛速度快，特别适用于多回路和交互作用比较强的情况。就求解速度而言，特别是对于大规模流程工业，联立方程法所占的优势就相对突出。
塔式太阳能热电系统具有模型复杂且方程、变量数多的特点，具体表现在建模时涉及到大量的系统模型方程组，包括物性计算方程、能量平衡方程、流体换热方程，并且其中有大量的微分方程。故而，传统的序贯模块法求解这种变量多、稀疏性强、高度非线性的流程模拟方程时速度较慢；联立方程法则有较快的求解速度，更适合对塔式太阳能热电系统整体模型的求解。
模拟的平台选用过程模拟软件gPROMS。它是目前在过程建模、仿真和优化应用中最先进和最通用的软件包之一。gPROMS采用联立方程法求解，并且直接面向具有大规模、非线性、稀疏性强的化工数学模型，计算功能强大，求解速度快。
根据建立的接收器、管道、储热罐和蒸汽发生器等主要子系统的动态数学模型，在gPROMS平台下建立起热电系统的各个子模型，然后按照输入输出关系通过模型连接方程进行联立。考虑系统在不同工况下的运行情况，按照不同运行模式时的流程结构分别进行连接，接着将DNI和传热工质流量作为模型输入，将蒸汽发生器功率和接收器出口传热工质温度作为模型输出，就可以得到塔式太阳能热电系统整体模型。
(3)考虑实时电价，以一天内电站收益最高为目标，热电系统参数的上下限为约束条件，构造优化问题。
电站正常运行时要保证系统有稳定的功率输出，但经济效益上不是最优的。考虑电价的预测，储热系统能把低需求时间段的热量转换到高需求时间段的热量，两个时间段的电价有着很大的差别，而对于储存的热量来说，这是一个升值的过程，这使得太阳能电站能够收到更大的收益。在给定当天每小时DNI数据以及市场电价的情况下，计划安排太阳能热电站每小时的发电量，使得当天电站所得收益最大化。
研究塔式太阳能热电系统正常运行的优化策略，满足热电系统安全运行要求的同时尽可能使得能量输出达到最大值，优化目标为一天内电站的收益最高。对市场电价与电价预测进行简化，在优化控制策略中使用已统计公布的历史电价数据作为预测得到的数据，而各个时段中电站产出的电量都能以当前时段的自由市场的电价卖出，并且市场电价每小时变动一次。因此第i小时的电站收益为该时段内的市场电价乘以该时段内的发电量：
η(i)＝Price(i)·Qi
其中，η(i)为电站第i小时的收益，Price(i)为第i小时的市场电价，Qi是电站第i小时的发电量。
当天总收益为各个小时收益的累加。目标函数可以表示为：
电站运行时需要满足整个系统的安全运行要求，需要考虑以下约束条件：
1)储热罐有固定的容量，当储热罐中储存的热量达到上限则无法继续储存热量，相反，当储热罐中储存的热量过少时也不能提取热量；2)发电装置有其最大发电功率，当发电装置达到满功率运作时，即使蒸汽发生器提供更多的高温蒸汽去驱动汽轮机，也无法汽轮机达到更大的发电功率；3)阀门开度有一定的范围，同时工质流量需要保证系统正常运行，因此各管路中的工质流量有上下限；4)考虑传热工质在不同温度下的稳定性、腐蚀性等特性，以及接收器的使用寿命等因素，传热工质的温度也有上下限。约束条件表示为：
储热罐储能约束:Qminch≤Qich≤Qmaxch
发电功率约束:0≤Pi≤Pmax
传热工质流量约束:mjmin≤mj≤mjmax
传热工质温度约束：Tfmin≤Tjfa≤Tfmax
其中，Qich为第i小时储热罐储存的能量，Qmaxch、Qminch为储热罐最大、最小储热量；Pi为第i小时的发电功率，Pmax为发电装置的最大发电功率；mj为第j根管道内传热工质的流量，mjmax、mjmin为第j根管道内传热工质的流量的上下限；Tjfa为第j根管道内传热工质的温度，Tfmax、Tfmin为传热工质温度的上下限。
因此，整个优化问题可以描述为：
(4)利用控制向量参数化方法来求解优化问题，得到最优的热电系统运行调度策略。
基于联立方程的塔式太阳能热电系统整体模型是复杂的非线性模型，存在各种约束条件，一般用间接方法和迭代动态规划方法难以解决这类模型的动态优化问题。控制向量参数化方法(CVP)是一种解决过程动态优化问题的有效方法。该方法首先将控制变量离散化为有限时间段；然后在优化时域内对系统模型积分，将优化问题中的系统模型约束完全消去，减少了非线性规划(NLP)中的变量个数；最后用NLP方法求解,以得到优化控制参数。
用于动态优化的CVP_SS方法基于CVP方法，并采用了打靶策略,该方法的思想是将两点边界值问题转化为初值问题进行计算，通过优化算法确定每个控制区间的控制变量值以及持续时间。
CVP_SS的优化计算部分描述如下：
设定决策变量(即控制变量)的初始值，其中，将连续型变量表示为x,而离散型的表示为y。
Step0：初始化，设置目标函数的当前可行值Φbest＝+∞，最优解对应的目标函数的值ΦLB＝+∞。
Step1：求解全松弛问题
①求解连续的NLP优化问题，将所有离散变量视为连续变量(即在其上下边界内取一定值)，以确定决策变量(包括连续型的x和离散型的y)的最优值和目标函数的对应值ΦFR；
②如果上述问题不可行，终止，即原问题不可行；
③如果在上述问题的解处，所有的优化决策变量有离散值，则结束(此时,原问题的优化解即为xFR、yFR，而目标函数的值为ΦFR)。
Step2：求解主问题
①构造混合整数线性规划问题，首先需要解决连续优化问题求解时目标函数和约束的线性化，并且排列所有离散型决策变量的组合；
②求解混整线性规划问题以确定连续和离散型决策变量的值(xMP，yMP)和对应的目标函数值ΦMP；
③如果上述问题不可行，或者目标函数的最优值Φbest和当前值ΦMP的差的绝对值小于εmax(1；|Φbest|)，则终止(如果Φbest＝+∞则原问题也不可行)，否则(xbest；ybest)为最优值，而Φbest为对应的目标函数值。
④由于混整线性规划问题改变了最优解的边界，因此需要更新目标值ΦLB:＝ΦMP。
Step3：求解原优化问题
①将所有需优化的离散决策变量固定于它们的当前值；
②再次求解连续优化问题(NLP)，以确定新的目标函数的值ΦPR和决策变量的优化值xPR；
③如果上述NLP可行并且ΦPR<Φbest，则更新Φbest＝ΦPR；xbest＝xPR；ybest＝yPR；
Step 4：进行迭代
①给定离散的决策变量的下一个解集，令yPR＝yMP；
②回Step 2；
对于NLP非线性规划问题,方法，借助gPROMS提供的基于SQP算法的SRQPD求解器和混合整数规划OAERAP求解器进行求解。CVP_SS算法终止时得到控制参数的最优解。
CVP_SS算法求解动态优化问题的实施流程如下：
步骤一：设置控制变量、控制间隔的初始值，确定时间节点数、约束变量上下限和优化时域；
步骤二：通过控制参数将控制变量离散化，将优化时域[0,tf]分割为N个时间单元，在各时间单元内用分段常数函数逼近连续控制量,即u(t)＝ui，其中ti-1≤t<ti,i＝1,2……N；
步骤三：在整个优化时域内确定各变量的动态过程；
步骤四：计算优化目标函数和约束变量，判断目标函数是否最优、约束条件是否满足。此时如果满足终止条件则算法终止转到步骤六，否则转到步骤五；
步骤五：计算目标函数关于控制参数的梯度，调整控制变量和控制间隔，转到步骤二；
步骤六：输出控制变量、约束变量的动态曲线以及优化目标值，得到最优的运行调度策略。
当算法终止时，得到最优的热电系统运行策略，确定了系统各主要变量在运行时间内的动态变化过程。此策略保证在市场电价下电站获得最高的收益，并使得发电功率、热能存储、工质流量以及工质温度都满足电站正常安全运行的要求，达到保护电站设备、提高系统稳定性的目的。同时算法结果也给出了主要变量的变化曲线、最高收益值等信息。
本发明实施例应用于西班牙的一个塔式太阳能热发电站，该电站以常压空气作为传热工质，其工作原理如图2所示。电站采用容积式开路空气接收器、单罐式储热系统以及螺旋管式直流蒸汽发生器，发电功率1MW，接收器出口温度设定值为700℃。模拟的时间为9:00-17:00，时长8小时。假设某日的市场电价如图3所示，在gPROMS平台下对电站运行进行模拟。采用稳定功率输出时，电站得到的总收益为230EUR；实施本发明提出的基于联立方程模型的塔式太阳能热电系统运行优化方法后，电站得到的总收益为270EUR，优化后的蒸汽发生器功率和储热罐能量的变化曲线如图4和图5所示。两者对比结果表明，经过运行优化后，电站能在市场电价下获得最高的收益，同时主要运行参数都能满足电站正常安全运行的要求。优化算法实施过程中，共进行了8次迭代。
Claims (1)
1.一种基于联立方程模型的塔式太阳能热电系统运行优化方法，其特征在于它的步骤如下：
(1)对热电系统各子系统进行分析，建立各子系统的动态数学模型；
所述的步骤(1)具体为：
塔式太阳能电站的热电系统包括的子系统有接收器、管道、储热罐和蒸汽发生器，在生产过程中，由于太阳辐射强度DNI因素和天气状况的变化，系统内的各变量也会随着时间变化，因此需要对各个子系统建立动态数学模型；分析各子系统的工作原理和换热方式，采用机理建模方法，用质量平衡、能量平衡和流体传热平衡方程，以及传热工质的物性参数方程进行描述，分别建立接收器、管道、储热罐和蒸汽发生器子系统的动态数学模型,
接收器模型表示为：
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&rho;</mi>
<mi>m</mi>
</msub>
<msub>
<mi>c</mi>
<mi>m</mi>
</msub>
<msub>
<mi>A</mi>
<mi>m</mi>
</msub>
<mi>&Delta;</mi>
<mi>x</mi>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>T</mi>
<mi>m</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>t</mi>
</mrow>
</mfrac>
<mo>=</mo>
<msub>
<mi>&xi;IA</mi>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>r</mi>
<mi>r</mi>
<mi>o</mi>
<mi>r</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>h</mi>
<mi>f</mi>
</msub>
<msub>
<mi>A</mi>
<mi>f</mi>
</msub>
<mi>&Delta;</mi>
<mi>x</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>T</mi>
<mi>m</mi>
</msub>
<mo>-</mo>
<msub>
<mi>T</mi>
<mi>f</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>&epsiv;&sigma;A</mi>
<mi>m</mi>
</msub>
<msubsup>
<mi>T</mi>
<mi>m</mi>
<mn>4</mn>
</msubsup>
<mo>-</mo>
<msub>
<mi>D</mi>
<mi>m</mi>
</msub>
<msub>
<mi>&pi;h</mi>
<mi>m</mi>
</msub>
<mi>&Delta;</mi>
<mi>x</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>T</mi>
<mi>m</mi>
</msub>
<mo>-</mo>
<msub>
<mi>T</mi>
<mi>a</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>A</mi>
<mi>f</mi>
</msub>
<msub>
<mi>h</mi>
<mi>f</mi>
</msub>
<mi>&Delta;</mi>
<mi>x</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>T</mi>
<mi>m</mi>
</msub>
<mo>-</mo>
<msub>
<mi>T</mi>
<mi>f</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>&rho;</mi>
<mi>f</mi>
</msub>
<msub>
<mi>c</mi>
<mi>f</mi>
</msub>
<msub>
<mi>A</mi>
<mi>f</mi>
</msub>
<mi>&Delta;</mi>
<mi>x</mi>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>T</mi>
<mi>f</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>t</mi>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<msub>
<mi>m</mi>
<mn>1</mn>
</msub>
<msub>
<mi>c</mi>
<mi>f</mi>
</msub>
<mi>&Delta;</mi>
<mi>x</mi>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>T</mi>
<mi>f</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
其中，Amirror为定日镜总面积，I为太阳辐射强度，ξ为定日镜场的效率，ε为黑度，σ为黑体辐射常数，Am为吸热器外部面积，Dm为吸热器外部直径，Tf为吸热器内工质温度，Ta为环境温度，Af为吸热器内部面积，m1为流过吸热器内的工质流量，hm、hf为体积对流传热系数；
管道模型表示为：
<mrow>
<mo>-</mo>
<mi>v</mi>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>T</mi>
<mi>f</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>&tau;</mi>
<mn>1</mn>
</msub>
</mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>T</mi>
<mi>w</mi>
</msub>
<mo>-</mo>
<msub>
<mi>T</mi>
<mi>f</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>T</mi>
<mi>f</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>t</mi>
</mrow>
</mfrac>
</mrow>
<mrow>
<mfrac>
<mn>1</mn>
<msub>
<mi>&tau;</mi>
<mn>12</mn>
</msub>
</mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>T</mi>
<mi>f</mi>
</msub>
<mo>-</mo>
<msub>
<mi>T</mi>
<mi>W</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>&tau;</mi>
<mn>2</mn>
</msub>
</mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>T</mi>
<mi>W</mi>
</msub>
<mo>-</mo>
<msub>
<mi>T</mi>
<mi>a</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>T</mi>
<mi>W</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>t</mi>
</mrow>
</mfrac>
</mrow>
其中，τ1是与流体密度、比热、管径相关的常数，τ12、τ2是与管壁金属密度、比热、管径相关的常数，Ta为管道外环境温度，Tf为管内流体温度，Tw为管壁温度；
储热罐模型表示为：
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&rho;</mi>
<mi>f</mi>
</msub>
<msub>
<mi>c</mi>
<mi>f</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>A</mi>
<mi>&Delta;</mi>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mi>&epsiv;</mi>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>T</mi>
<mi>f</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>t</mi>
</mrow>
</mfrac>
<mo>=</mo>
<msub>
<mi>mc</mi>
<mi>f</mi>
</msub>
<msub>
<mi>T</mi>
<mi>f</mi>
</msub>
<mo>-</mo>
<msub>
<mi>mc</mi>
<mi>f</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>T</mi>
<mi>f</mi>
</msub>
<mo>+</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>T</mi>
<mi>f</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<mi>&Delta;</mi>
<mi>x</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<msub>
<mi>h</mi>
<mi>v</mi>
</msub>
<mi>A</mi>
<mi>&Delta;</mi>
<mi>x</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>T</mi>
<mi>b</mi>
</msub>
<mo>-</mo>
<msub>
<mi>T</mi>
<mi>f</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>U</mi>
<mi>D</mi>
<mi>&pi;</mi>
<mi>&Delta;</mi>
<mi>x</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>T</mi>
<mi>a</mi>
</msub>
<mo>-</mo>
<msub>
<mi>T</mi>
<mi>f</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mrow>
<msub>
<mi>&rho;</mi>
<mi>b</mi>
</msub>
<msub>
<mi>c</mi>
<mi>b</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>A</mi>
<mi>&Delta;</mi>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mi>&epsiv;</mi>
<mo>)</mo>
</mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>T</mi>
<mi>b</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>t</mi>
</mrow>
</mfrac>
<mo>=</mo>
<msub>
<mi>h</mi>
<mi>v</mi>
</msub>
<mi>A</mi>
<mi>&Delta;</mi>
<mi>x</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>T</mi>
<mi>f</mi>
</msub>
<mo>-</mo>
<msub>
<mi>T</mi>
<mi>b</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
其中，A为填充床的横截面积，m为流过储热罐的空气流量，Tf为储热罐内的工质温度，U为流体与环境之间的传热系数，Ta为环境温度，D为储热罐的直径；
蒸汽发生器模型表示为：
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>D</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>t</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<mo>&part;</mo>
<mrow>
<mo>&part;</mo>
<mi>l</mi>
</mrow>
</mfrac>
<mrow>
<mo>(</mo>
<mfrac>
<msup>
<mi>D</mi>
<mn>2</mn>
</msup>
<mi>&rho;</mi>
</mfrac>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>P</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>t</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mfrac>
<msup>
<mi>D</mi>
<mn>2</mn>
</msup>
<mrow>
<mn>2</mn>
<mi>&rho;</mi>
</mrow>
</mfrac>
<msub>
<mi>f</mi>
<mi>f</mi>
</msub>
<mfrac>
<mi>S</mi>
<mi>A</mi>
</mfrac>
<mo>-</mo>
<mi>&rho;</mi>
<mi>g</mi>
</mrow>
式中，ρ为流体密度，D为单位面积的质量流速，h为流体比焓，Q为单位长度的热流率，A为流道的横截面积，P为流体压力；
(2)根据各子系统的动态数学模型，采用联立方程法建立塔式太阳能热电系统整体模型；
所述的步骤(2)具体为：
针对塔式太阳能热电系统具有的模型复杂且方程、变量数多的特点，采用联立方程法对系统整体模型进行求解；
模拟平台选用过程模拟软件gPROMS，根据建立的接收器、管道、储热罐和蒸汽发生器子系统的动态数学模型，在gPROMS平台下建立热电系统的各个子模型，然后按照输入输出关系通过模型连接方程进行联立，考虑系统在不同工况下的运行情况，按照不同运行模式时的流程结构分别进行连接，接着将DNI和传热工质流量作为模型输入，将蒸汽发生器功率和接收器出口传热工质温度作为模型输出，就得到塔式太阳能热电系统整体模型；
(3)考虑实时电价，以一天内电站收益最高为目标，热电系统参数的上下限为约束条件，构造优化问题；
所述的步骤(3)具体为：
考虑对电价的预测，储热罐把低需求时间段的热量转换到高需求时间段的热量，使储存的热量升值，也使太阳能电站能够得到更大的收益；在给定当天每小时DNI数据以及市场电价的情况下，计划安排太阳能热电站每小时的发电量，使得当天电站所得收益最大化；
研究塔式太阳能热电系统正常运行的优化策略，满足热电系统安全运行要求的同时尽可能使得能量输出达到最大值，优化目标为一天内电站的收益最高；在优化控制策略中使用已统计公布的历史电价数据作为预测得到的数据，假设各个时段中电站产出的电量都能以当前时段的自由市场的电价卖出，并且市场电价每小时变动一次，因此第i小时的电站收益为该时段内的市场电价乘以该时段内的发电量：
η(i)＝Price(i)·Qi
其中，η(i)为电站第i小时的收益，Price(i)为第i小时的市场电价，Qi是电站第i小时的发电量；
当天总收益为各个小时收益的累加，目标函数可以表示为：
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
<mi>&eta;</mi>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<mi>Pr</mi>
<mi>i</mi>
<mi>c</mi>
<mi>e</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<msub>
<mi>Q</mi>
<mi>i</mi>
</msub>
</mrow>
电站运行时需要满足整个系统的安全运行要求，需要考虑以下约束条件：
1)储热罐有固定的容量，当储热罐中储存的热量达到上限则无法继续储存热量，相反，当储热罐中储存的热量过少时也不能提取热量；2)发电装置有其最大发电功率，当发电装置达到满功率运作时，即使蒸汽发生器提供更多的高温蒸汽去驱动汽轮机，也无法汽轮机达到更大的发电功率；3)阀门开度有一定的范围，同时工质流量需要保证系统正常运行，因此各管路中的工质流量有上下限；4)考虑传热工质在不同温度下的稳定性、腐蚀性，以及接收器的使用寿命，因此传热工质的温度也有上下限；约束条件表示为：
储热罐储能约束:Qminch≤Qich≤Qmaxch
发电功率约束:0≤Pi≤Pmax
传热工质流量约束:mjmin≤mj≤mjmax
传热工质温度约束：Tfmin≤Tjfa≤Tfmax
其中，Qich为第i小时储热罐储存的能量，Qmaxch、Qminch为储热罐最大、最小储热量，Pi为第i小时的发电功率，Pmax为发电装置的最大发电功率，mj为第j根管道内传热工质的流量，mjmax、mjmin为第j根管道内传热工质的流量的上下限，Tjfa为第j根管道内传热工质的温度，Tfmax、Tfmin为传热工质温度的上下限；
因此，整个优化问题描述为：
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>max</mi>
<mi>&eta;</mi>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<mi>Pr</mi>
<mi>i</mi>
<mi>c</mi>
<mi>e</mi>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<msub>
<mi>Q</mi>
<mi>i</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mtable>
<mtr>
<mtd>
<mrow>
<mi>s</mi>
<mo>.</mo>
<mi>t</mi>
<mo>.</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>Q</mi>
<mrow>
<mi>min</mi>
<mi>c</mi>
<mi>h</mi>
</mrow>
</msub>
<mo>&le;</mo>
<msub>
<mi>Q</mi>
<mrow>
<mi>i</mi>
<mi>c</mi>
<mi>h</mi>
</mrow>
</msub>
<mo>&le;</mo>
<msub>
<mi>Q</mi>
<mrow>
<mi>max</mi>
<mi>c</mi>
<mi>h</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow>
<mn>0</mn>
<mo>&le;</mo>
<msub>
<mi>P</mi>
<mi>i</mi>
</msub>
<mo>&le;</mo>
<msub>
<mi>P</mi>
<mi>max</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>m</mi>
<mrow>
<mi>j</mi>
<mi>min</mi>
</mrow>
</msub>
<mo>&le;</mo>
<msub>
<mi>m</mi>
<mi>j</mi>
</msub>
<mo>&le;</mo>
<msub>
<mi>m</mi>
<mrow>
<mi>j</mi>
<mi>max</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>T</mi>
<mrow>
<mi>f</mi>
<mi>min</mi>
</mrow>
</msub>
<mo>&le;</mo>
<msub>
<mi>T</mi>
<mrow>
<mi>j</mi>
<mi>f</mi>
<mi>a</mi>
</mrow>
</msub>
<mo>&le;</mo>
<msub>
<mi>T</mi>
<mrow>
<mi>f</mi>
<mi>max</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>;</mo>
</mrow>
(4)利用控制向量参数化方法求解优化问题，得到最优的热电系统运行调度策略；
所述的步骤(4)具体为：
针对过程动态优化问题，采用控制向量参数化方法进行解决，该方法将控制变量离散化为有限时间段，从而把原问题转化为NLP问题进行求解，而CVP_SS方法基于控制向量参数化方法，并采用了打靶策略；应用CVP_SS方法求解塔式太阳能热电系统运行优化问题，首先根据实际需要设置控制变量、控制间隔的初始值，确定时间节点数、约束变量上下限和优化时域，然后将两点边界值问题转化为初值问题进行计算，通过优化算法确定每个控制区间的控制变量值以及持续时间，对于NLP非线性规划问题，借助gPROMS提供的基于SQP算法的SRQPD求解器和混合整数规划OAERAP求解器进行求解；
当算法终止时，得到最优的热电系统运行策略，确定了系统各主要变量在运行时间内的动态变化过程，此策略保证在市场电价下电站获得最高的收益，并满足电站正常安全运行的要求，同时算法结果也给出了主要变量的变化曲线、最高收益值信息；
CVP_SS算法求解动态优化问题的实施流程如下：
步骤一：设置控制变量、控制间隔的初始值，确定时间节点数、约束变量上下限和优化时域；
步骤二：通过控制参数将控制变量离散化，将优化时域[0,tf]分割为N个时间单元，在各时间单元内用分段常数函数逼近连续控制量,即u(t)＝ui，其中ti-1≤t<ti,i＝1,2……N；
步骤三：在整个优化时域内确定各变量的动态过程；
步骤四：计算优化目标函数和约束变量，判断目标函数是否最优、约束条件是否满足；此时如果满足终止条件则算法终止转到步骤六，否则转到步骤五；
步骤五：计算目标函数关于控制参数的梯度，调整控制变量和控制间隔，转到步骤二；
步骤六：输出控制变量、约束变量的动态曲线以及优化目标值，得到最优的运行调度策略。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410734256.5A CN104408534B (zh) | 2014-12-05 | 2014-12-05 | 基于联立方程模型的塔式太阳能热电系统运行优化方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410734256.5A CN104408534B (zh) | 2014-12-05 | 2014-12-05 | 基于联立方程模型的塔式太阳能热电系统运行优化方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104408534A CN104408534A (zh) | 2015-03-11 |
| CN104408534B true CN104408534B (zh) | 2018-04-20 |
Family
ID=52646163
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410734256.5A Active CN104408534B (zh) | 2014-12-05 | 2014-12-05 | 基于联立方程模型的塔式太阳能热电系统运行优化方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104408534B (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US11777441B2 (en) | 2021-04-01 | 2023-10-03 | Hamilton Sundstrand Corporation | Thermoelectric power generation using radiant and conductive heat dissipation |
Families Citing this family (15)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105406767B (zh) * | 2015-10-12 | 2017-05-17 | 武汉理工大学 | 一种汽车尾气热电转化装置电气拓扑结构的优化方法 |
| CN105205562B (zh) * | 2015-11-19 | 2019-03-29 | 浙江大学 | 塔式太阳能电站接收器的运行优化方法 |
| EP3184914A1 (en) * | 2015-12-21 | 2017-06-28 | E.ON Sverige AB | A thermal server plant and a method for controlling the same |
| CN106096085B (zh) * | 2016-05-30 | 2019-07-30 | 浙江大学 | 一种塔式太阳能热发电系统的辐射能密度模拟方法 |
| CN107247843A (zh) * | 2017-06-10 | 2017-10-13 | 中国科学院电工研究所 | 太阳能热发电站油水型蒸汽发生器模型的建模方法 |
| CN107632962B (zh) * | 2017-08-04 | 2023-12-19 | 中国电力科学研究院 | 一种塔式光热电站的仿真分析方法及其系统 |
| CN108984875B (zh) * | 2018-07-02 | 2022-11-08 | 中国电力工程顾问集团西北电力设计院有限公司 | 一种计算光热机组参与电力平衡容量的方法 |
| CN108983839A (zh) * | 2018-07-25 | 2018-12-11 | 浙江工业大学 | 一种带uwb定位的调温系统及方法 |
| CN109871607B (zh) * | 2019-02-14 | 2023-04-18 | 浙江可胜技术股份有限公司 | 一种汽轮机定功率输出方法 |
| CN109829654B (zh) * | 2019-02-14 | 2020-12-22 | 浙江中控太阳能技术有限公司 | 一种匹配吸热器运行的镜场调度方法 |
| CN110334878B (zh) * | 2019-07-11 | 2022-05-03 | 国网甘肃省电力公司经济技术研究院 | 一种基于典型静态模型的光热储能电站发电量优化方法 |
| CN112528542B (zh) * | 2020-12-15 | 2022-09-20 | 河海大学 | 一种考虑集热器工作模式的光热发电系统动态建模方法及控制方法 |
| CN113420428B (zh) * | 2021-06-09 | 2022-04-22 | 中国科学院理化技术研究所 | 固相蓄冷介质筛选方法及液态空气储能电站的设计方法 |
| CN113946935B (zh) * | 2021-08-17 | 2023-04-14 | 中国能源建设集团新疆电力设计院有限公司 | 一种基于并行计算的槽式聚光集热系统性能算法 |
| CN113983702B (zh) * | 2021-10-28 | 2022-11-25 | 东北大学 | 一种光-储-热供电机组及其控制方法 |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103093017A (zh) * | 2011-11-04 | 2013-05-08 | 新奥科技发展有限公司 | 分布式能源系统设计方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20130151022A1 (en) * | 2011-12-07 | 2013-06-13 | General Electric Company | Systems and Methods for Assessing Future Power Plant Capabilities |
-
2014
- 2014-12-05 CN CN201410734256.5A patent/CN104408534B/zh active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103093017A (zh) * | 2011-11-04 | 2013-05-08 | 新奥科技发展有限公司 | 分布式能源系统设计方法 |
Non-Patent Citations (4)
| Title |
|---|
| 具有切换结构的非线性系统最优控制方法研究;赵瑞艳;《中国博士学位论文全文数据库基础科学辑》;20111015(第10期);第A002-89页 * |
| 基于综合性能指标的天然气冷热电联供系统优化运行与配置;王庆华;《中国优秀硕士学位论文全文数据库工程科技II辑》;20121215(第12期);第C042-78页 * |
| 基于联立方程的塔式太阳能热电系统模拟;任涛 等;《计算机与应用化学》;20141128;第31卷(第11期);第1-6页 * |
| 风电-储能联合日运行优化调度模型;程世军等;《电网与清洁能源》;20140425;第30卷(第4期);第68-80页 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US11777441B2 (en) | 2021-04-01 | 2023-10-03 | Hamilton Sundstrand Corporation | Thermoelectric power generation using radiant and conductive heat dissipation |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104408534A (zh) | 2015-03-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104408534B (zh) | 基于联立方程模型的塔式太阳能热电系统运行优化方法 | |
| Chen et al. | Low-carbon economic dispatch of integrated energy system containing electric hydrogen production based on VMD-GRU short-term wind power prediction | |
| CN111144668B (zh) | 考虑场景模拟的综合能源系统随机优化模型的建立方法 | |
| Ammar et al. | Artificial neural network based control for PV/T panel to track optimum thermal and electrical power | |
| Sharma | A particle swarm optimization algorithm for optimization of thermal performance of a smooth flat plate solar air heater | |
| Zheng et al. | Bio-inspired optimization of sustainable energy systems: a review | |
| Osorio et al. | Forecasting solar-thermal systems performance under transient operation using a data-driven machine learning approach based on the deep operator network architecture | |
| Qi et al. | Optimal configuration of concentrating solar power in multienergy power systems with an improved variational autoencoder | |
| Wu et al. | Two-stage robust optimization model for park integrated energy system based on dynamic programming | |
| Sharma et al. | Optimization of a smooth flat plate solar air heater using stochastic iterative perturbation technique | |
| Kuang et al. | Stochastic dynamic solution for off-design operation optimization of combined cooling, heating, and power systems with energy storage | |
| He et al. | The many-objective optimal design of renewable energy cogeneration system | |
| Guo et al. | Multi-objective optimization of multi-energy heating systems based on solar, natural gas, and air-energy | |
| CN112347607A (zh) | 一种基于凸松弛的热电联合调度方法 | |
| CN107832873B (zh) | 基于双层母线式结构的综合能源系统优化规划方法及装置 | |
| Vasilj et al. | Robust nonlinear economic MPC based management of a multi energy microgrid | |
| CN113363976B (zh) | 一种基于场景图的风光水互补发电系统中期优化调度方法 | |
| CN116432824A (zh) | 基于多目标粒子群的综合能源系统优化方法及系统 | |
| CN112886571B (zh) | 一种基于边界变量可行域的电热综合能源系统分解协调优化运行方法及装置 | |
| CN112084680B (zh) | 一种基于dqn算法的能源互联网优化策略方法 | |
| CN112214883B (zh) | 一种微能源网综合能源系统运行优化方法和系统 | |
| CN105205562A (zh) | 塔式太阳能电站接收器的运行优化方法 | |
| CN116207796B (zh) | 计及可变流量调节的电-热综合能源系统调度方法及系统 | |
| de Oliveira et al. | Dynamic online optimization of a house heating system in a fluctuating energy price scenario | |
| Wang et al. | Rolling optimization based on holism for the operation strategy of solar tower power plant |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |