• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Optimal Operation of Integrated Heat and Electricity Systems:A Tightening McCormick Approach

    2022-01-26 06:39:14LirongDengHongbinSunBojuLiYongSunTinshuYngXunZhng
    Engineering 2021年8期

    Lirong Deng,Hongbin Sun,b,Boju Li,Yong Sun,Tinshu Yng,Xun Zhng,*

    a Tsinghua–Berkeley Shenzhen Institute,Tsinghua University,Shenzhen 518055,China

    b State Key Laboratory of Power Systems,Department of Electrical Engineering,Tsinghua University,Beijing 100084,China

    c State Grid Jilin Electric Power Supply Co.,Ltd.,Changchun 130000,China

    d Electrical College,Harbin Institute of Technology,Harbin 150006,China

    ABSTRACT

    Keywords:Integrated heat and electricity system Convex relaxation Operation McCormick envelopes

    1.Introduction

    1.1.Background and motivatio n

    The electric power system(EPS)is a crucial part of the national energy supply,and policymakers are gradually realizing that heat supply also plays a key role in the energy system.The International Energy Agency reports that over half of global energy use is for heating[1].A large proportion of the heat supply in high-density areas comes from district heating systems(DHSs),in comparison with the other heat alternatives—such as individual heat pumps,gas boilers,solar thermal heating,and electrical heating—that are often applied in low-density areas.Although different countries and areas use these heat supply methods in different proportions,the district heating approach has been proven to be extremely energy efficient[2].

    Electricity and heat can be produced simultaneously with centralized energy generation and district heating infrastructures.In general,those two large energy systems—EPS and DHS—are tightly connected by combined heat and power(CHP,also known as cogeneration)plants and power-to-heat facilities.By 2050,CHP will provide 26% of electricity in the European Union.In Denmark,the government aims to achieve 100%renewable heat and electricity production by 2035[3].It is expected that EPSs and DHSs will impact each other to a greater extent in the energy production and consumption process in the near future.Therefore,the idea of the integrated heat and electricity system has attracted interest from practitioners and researchers.

    Studies on different aspects of integrated heat and electricity systems are emerging,including modeling[4],state estimation[5],unit commitment[6],economic dispatch[7–9],market mechanisms[10,11],and planning[12,13].Among them,modeling plays a fundamental and substantial role in the commercialization of integrated heat and electricity systems.Although EPS modeling has been studied thoroughly during the past decades,a great deal of research on DHS modeling is still ongoing.In general,DHS has three regulation modes,as shown in Table 1.Quality regulation employs a constant mass flow rate and variable temperature strategy[7].Hence,the constraints in related optimization problems become linear and are thus easy to handle.However,a situation in which hydraulic conditions are predetermined may not lead to an economical solution.In contrast,quantity regulation maintains a constant supply temperature but regulates the quantity of mass flow rates.It has the advantages of flexibility and cost reduction because the mass flow rates will change according to the heat load,which will reduce the power consumption of the circulating water pumps.Apparently,greater economic efficiency and flexibility can be achieved by regulating both the temperature and the mass flow rates,which is referred to as quality–quantity regulation.So far,only a few works have explored the operation of integrated energy systems with a variable mass flow strategy[14–18].Furthermore,efficient and scalable algorithms to solve the nonconvex and nonlinear network flow constraints are still under investigation.To this end,this paper proposes a convex model and an efficient algorithm for integrated heat and electricity systems under the quality–quantity regulation,in which the algorithm is expected to find the global optimum or a near-global-optimal feasible solution with a relatively small computational burden.

    Table 1 Different regulation modes in DHS.

    1.2.Literature review

    The modeling of quality–quantity-regulated DHS is a type of pooling problem[19].The pooling problem is a network flow problem that aims to find the minimum-cost way to mix several inputs in intermediate pools such that the output meets the demand or certain requirements.Mixing inputs involves mixing the product of the flow quantity and the feature.As a result,the pooling problem becomes a bilinear program.In DHS modeling,the nonconvex and nonlinear network flow hinders the problem from being solved efficiently.Among the nonconvex terms,the bilinear term,which comes from the product of the mass flow rate and the nodal temperature,is one of the most difficult to deal with.

    Current achievements in dealing with the bilinearity of DHS optimization can be divided into four categories:nonlinear programming methods,generalized Benders decomposition,relaxation methods,and relaxation tightening methods.Nonlinear programming methods,such as the interior point method,sequential linear programming,and successive quadratic programming[14],are generally able to solve nonlinear programing with continuous variables,and are easy to implement with off-the-shelf solvers.However,they are only able to find local solutions and may have slow convergence or even convergence failure when the network becomes large.

    Generalized Benders decomposition can solve certain types of nonlinear programming and mixed-integer nonlinear program-ming.To deal with the bilinearity from the product of the mass flow rate and the nodal temperature,Ref.[15]presents an iterative algorithm:By fixing one set of variables and solving the other set,and then vice versa,the two subproblems are iteratively solved until convergence is reached.However,the resulting solution cannot be ensured to be a global optimum or local optimum,and the convergence property requires further investigation.

    In general,relaxation methods,such as conic relaxation[16,17]and polyhedral relaxation[18],enlarge the original nonconvex feasible set until it is convex.The relaxed problem becomes convex at the expense of sacrificing the feasibility of the solution to the original problem.The performance of relaxation methods largely depends on the relaxed boundaries,with tighter bounds leading to stronger relaxations.This gives rise to relaxation tightening methods,such as those shown in Refs.[18–20].The core principle of such methods is to provide tighter bounds to enhance the relaxation.Among the relaxation tightening methods,branch and bound has been successfully implemented in the latest version of Gurobi—that is,Gurobi 9.0[21]—to deal with bilinear programming.The bilinear solver in Gurobi 9.0 ensures a global optimum and can be used as a benchmark to evaluate the optimality of other methods.However,when dealing with large-scale problems,it may converge slowly.

    1.3.Contributions and the organization of this paper

    In this paper,we first reformulate the classic quality–quantityregulated DHS optimization model through an equivalent transformation and the first-order Taylor expansion.Compared with the original model,which has two bilinear terms in each nonconvex constraint,the reformulated model has fewer bilinear equalities;more specifically,it reduces the bilinear terms by approximately half.The reformulation not only ensures optimality,but also may have advantages in reducing the computational complexities of the original problem.Next,we perform McCormick envelopes to convexify the bilinear constraints and obtain an objective lower bound of the reformulated model.To improve the quality of the McCormick relaxation,we employ a piecewise McCormick technique to reduce the McCormick volume by deriving stronger lower and upper bounds of the bilinear terms.The piecewise McCormick technique involves partitioning the domain of one of the variables in the bilinear term into several disjoint regions and determining the optimal region such that the bounds of the selected variable have been tightened.Thus,a strengthened lower-bound solution of the original problem is obtained.Since the strengthened lower bound may not be feasible,a heuristic bound contraction algorithm is further established to constrict the strengthened bounds of the piecewise McCormick technique and obtain a nearby feasible solution in an iterative way.Compared with nonlinear optimization[14]and generalized Benders decomposition[15],the proposed method is based on relaxation and piecewise techniques,which are more scalable,and thus avoids being trapped into local infeasibility.In contrast to current relaxation tightening methods,such as the one implemented in the Gurobi bilinear solver[21],the proposed tightening McCormick algorithm can achieve a comparable solution quality with better computational efficiency.

    In summary,this paper makes three main contributions:

    (1)A classic integrated heat and electricity operation problem with a quality–quantity-regulated DHS model is reformulated through a variable substitution and equivalent transformation,which largely reduces the bilinear complexity of the classic model.

    (2)Based on the reformulated model,McCormick envelopes are applied to relax the remaining bilinear terms.To reduce the relaxation errors and strengthen the bounds used by the McCormick envelopes,a piecewise McCormick technique is proposed.

    (3)To improve the feasibility of the solution,a bound contraction algorithm is designed to constrict the upper and lower bounds of the piecewise McCormick envelopes through a perturbation of the latest optimal results.Case studies show that the tightening McCormick method,which combines the piecewise McCormick technique with a bound contraction algorithm,quickly solves the problem with an acceptable feasibility check and optimality.Given its convex property,the tightening McCormick method is promising for large-scale integrated heat and electricity optimization and can allow economic analysis using shadow prices.

    The remainder of this paper is organized as follows.Section 2 introduces the integrated operation problem with a DHS base model and a reformulated model.Section 3 details the convex relaxation and the tightening McCormick algorithm used to solve the integrated problem.Section 4 presents case studies,and Section 5 draws conclusions.

    2.Problem formulation

    An integrated heat and electricity system consists of DHSs and local EPSs.A local EPS is part of the whole multi-area power system,which has one or more interfaces to exchange power with other parts of the whole power system.For DHS modeling,we present a base model first,which is a nonlinear optimization model without any relaxation of the constraints.Then,we derive an equivalent reformulated model of the DHS through the firstorder Taylor expansion and variable substitution.The reformulated model turns out to reduce the nonlinearity of the base model;that is,the bilinear terms in the base model are reduced by nearly half.For EPS formulation,we adopt a state-independent linear power flow model that considers the reactive power and voltage magnitude.The coupling elements connecting the two systems are CHP units.The operation problem determines optimal solutions with multiple objectives such as the minimization of fuel costs,power transaction costs with other interconnected power systems,and network losses.For notational simplicity,we build single-time horizon models for the DHS and EPS,respectively.Then,we extend the whole operation optimization to the multi-time horizon case and add the notation of time indextto all the decision variables.

    2.1.District heating system modeling

    In the quality–quantity-regulated radial district heating network,we regulate the mass flow rates through the circulating pump.Mass flow rate has a magnitude and direction.In this paper,we assume the magnitude of the mass flow rate to be variable,while the direction is fixed.This assumption is reasonable because frequent changes of direction will result in supply instability.

    Base model:

    whereHGiis the heat generation at nodei,HLiis the heat load at nodei,cis the specific heat capacity of water,mjiis the mass flow rate of water transferred from nodejto nodeiin the pipeline of the heating network,τiis the outlet temperature of nodei,τjiis the outlet temperature of the pipeline from nodejto nodei,τais the ambient temperature,v is the heat transfer coefficient per unit length,Ljiis the length of the pipeline from nodejto nodei,is the lower/upper limit of τji,andis the lower/upper limit ofmji.In(i)is the set of indices of nodes flowing into nodeiin the heating network,and Out(i)is the set of indices of nodes flowing out of nodeiin the heating network.Inodeand Ipipeare the set of indices of nodes and pipelines in the heating network,respectively.Eq.(1a)defines the nodal heat balance,while Eq.(1b)defines the nodal flow balance.Eq.(1c)describes the process by which the water temperature drops along the pipeline,considering the heat-loss factors[22].To be specific,the outlet temperature of the pipeline τjirelies on the outlet temperature at the starting point of the pipeline,τj.If the pipeline length is longer,the heat transfer coefficient is larger,or the mass flow rate is smaller,the temperature drop(as well as the heat loss)will become more significant.Eq.(1d)are the minimum-and maximum-operating limits of the nodal outlet temperature and the pipeline outlet temperature.Eq.(1e)gives the minimum-and maximum-operating limits of the mass flow rates.Note that additional constraints that keep the supply nodal temperatures constant can be added into the base model;these would constitute the quantity-regulated model,whose computational complexity remains the same as the quality–quantity-regula ted model.

    2.2.DHS reformulation

    Nonconvexities of the DHS base model arise due to Eqs.(1a)and(1c).Eq.(1a)has the bilinear termmτ,and Eq.(1c)has the exponential function exp(–vL/(cm)).

    ·To deal with the bilinear term,we introduce an auxiliary variableH = mτ.With this variable substitution,it turns out that the bilinear terms are reduced.

    · To deal with the exponential function term,we use the firstorder Taylor expansion to approximate the constraints in Eq.(1c).

    By introducing auxiliary variables

    Eqs.(1a)and(1c)are equivalent to

    In practice,the total heat transfer coefficient of pipeline v is very small.According to the Chinese design code for a city heating network[23],the heat transfer coefficient of the thermal insulation material should be less than 0.08 W·(m·K)-1.Thus,we assume that vLji?cmji—an assumption that has been broadly studied and justified for heating networks[16,17,20].Under this assumption,we can use the first-order Taylor expansion to approximate Eq.(5).

    To transform the upper and lower limit constraints associated with temperature(i.e.,Eq.(1d))into constraints related to the heat powerH,the following constraints are used:

    Therefore,the district heating network formulation under quality–quantity regulation(the base model)can be reformulated as Eqs.(9a)–(9g),which represent the reformulated model.

    Reformulated model:

    It should be noted that Eq.(2)does not appear in the reformulated model.It has been eliminated becausecompletely represents τji.Introducingand eliminating variables τji(Eq.(2))will not influence the feasible region of the original problem.However,Eq.(9f)does not fall under the same case because τiis not only presented in Eq.(9f),but also appears in Eq.(9b).In the reformulated model(Eqs.(9a)–(9g)),the nonconvex quadratic constraints(Eq.(1a))are converted into linear constraints(Eq.(9a))and independent bilinear constraints(Eq.(9f));thus,all constraints are linear except for Eq.(9f).It can be argued that the reformulated model is equivalent to the base model,with negligible errors from the first-order Taylor approximation.

    2.3.Electric power system modeling

    A linear power flow[24]with accurate estimation of voltage magnitude is adopted to characterize the electric power network;this improves the direct current(DC)power flow results,since it accounts for the reactive power and voltage magnitude.

    2.4.Energy sources modeling

    There are three typical types of e nergy sources in the integrated system:heating boilers(HBs),CHP units,and non-CHP thermal units(TUs).The variable cost of an HB is typically expressed as a linear function with respect to the heat output,whereas the variable cost may include fuel prices and taxes[25].The generation cost of a CHP unit is usually formulated as a convex quadratic cost function of the power and heat output,including the product of the power and heat output[26].The cost of each non-CHP thermal unit is modeled as either a quadratic or piecewise linear function in MATPOWER[27],and the quadratic form has been selected in this paper.The objectives of the above energy sources are described as follows:

    2.5.The operation problem

    The optimal operation problem in the integrated heat and electricity system is to minimize the fuel costs,the power exchange costs with other interconnected power systems,the network losses,and some other appropriate objective functions.The constraints include Eqs.(1a)–(1e)for the base model or Eqs.(9a)–(9g)for the reformulated model,in addition to Eq.(10),Eq.(12),and three other constraints regarding the nodal electricity/heat production equalities:

    The problem(Eqs.(13a)–(13e)),whether with the DHS base model(Eqs.(1a)–(1e))or the reformulated model(Eqs.(9a)–(9g)),is a nonconvex program with quadratic constraints.It is non-deterministic polynomial hard.

    3.Convex relaxation and the solution algorithm

    In the DHS reformulated model,the bilinear term is the only difficulty.An intuitive idea is to remove it—that is,remove Eq.(9f)—directly in order to check whether the bilinear constraints have a great impact on the quality of the solution.However,simulations show that removing bilinear terms leads to solutions with large violation errors to the constraints Eq.(9f).A classic way to relax the bilinear terms is to use McCormick envelopes.To improve the relaxation and reduce the McCormick volume,a piecewise McCormick technique is presented,since the sum of the smaller volumes is smaller than the big volume in the standard McCormick approach.The optimal results from the piecewise McCormick technique can be regarded as the relaxed global optimum,but may not be feasible for the reformulated model.Hence,a bound contraction algorithm is further proposed in order to strengthen the boundaries of the McCormick envelopes and find a feasible result near the relaxed global optimum.

    3.1.McCormick convex relaxation

    We replace Eq.(9f)with the McCormick relaxation[31].

    McCormick Model:

    Under the McCormick relaxation,the reformulated model becomes a convex optimization problem,where the Karush–Kuh n–Tucker conditions are sufficient and necessary(under e.g.,Slater’s conditions).Therefore,a global minimum can be achieved in the relaxed McCormick model.This global minimum can be regarded as a lower bound to the reformulated model.However,the McCormick relaxation still induces relatively large errors in the bilinear constraints.

    3.2.The tightening McCormick algorithm

    We have obtained a solution through the convex McCormick model.Since the relaxed feasible region is not tight in the McCormick model,the solution may not be feasible due to violation of the bilinear constraints.To improve the quality of the relaxation,tighter lower and upper bounds are expected to help construct the McCormick envelopes.As such,we employ a piecewise McCormick technique and a bound contraction algorithm to improve the McCormick method.

    3.2.1.The piecewise McCormick technique

    In the piecewise McCormick technique,the domain of one of the variables in the bilinear term is partitioned into several disjoint regions and the optimal region is determined such that the bounds of the selected variable have been tightened.A typical partition pattern is the uniform partition[32],where the problem size increases linearly with the number of partitions.Other partition schemes with adaptive segment lengths or partition-dependent bounds can also be applied to improve the performance[20].

    The selected variables to be partitioned will influence the quality of the relaxation.In the district heating network flow optimization problem,the variable choice could be mass flow rates,nodal temperatures,or a combination of both.Finding an optimal set of variables for the best relaxation is beyond the scope of this paper.We will choose the nodal temperatures to partition.Letandrepresent the lower and upper bounds of variable τiin partitions.A binary variableyi,sis assigned to each disjunctions.yi,s=1 if the value of τidoes belong to this disjunction;otherwise,yi,s=0.The other variable of the bilinear term,mik,is also disaggregated asmik,s,?s∈S,whereSis the set of indices ofs.

    Piecewise McCormick:

    In the above formulation,if the binary variableyi,sis equal to 1,then all the variables in thesth partition—that is,τi,sandmik,s—play dominant roles in deciding the values of τiandmik.In comparison,the variables in all the other partitions are enforced to zero.Similarly,ifyi,s= 1,then all the constraints in thesth partition are enforced,while the constraints in all the other partitions are neglected.The increasing number of binaries will lead to stronger relaxations,but at the cost of increased computational effort to solve the resulting mixed-integer problem.In general,the algorithm performs well with three partitions[32].

    3.2.2.The bound contraction algorithm

    The piecewise McCormick technique presents tighter upper bounds and lower bounds of the nodal temperatures,resulting in a stronger lower-bound solution of the reformulated model.To further reduce errors of the bilinear constraints,a feasible solution near the lower-bound solution is expected.Hence,we propose a bound contraction algorithm to iteratively strengthen the variable bounds and approach near-optimal results that violate the bilinear constraints less.The main idea of the bound contraction algorithm is that at each iterationn,the upper and lower bounds of the decision variables—that is,the mass flow rates and nodal temperatures—are updated according to the results from the last iterationn–1 and a sequence of hyperparameters ε,where 0 <ε <1.The principle in choosing the sequence is to gradually reduce the value of ε by κ in order to tighten the bounds[33].To ensure the feasibility of the original problem,the updated bounds should be the intersection of the result-orientated bounds and the original bounds.Once the average relative error of the bilinear constraints(Eq.(9f))is reached to an acceptable level δ,the algorithm can be terminated.

    The abovementioned procedure is to search for a tighter lower bound.However,there is a need to obtain a feasible solution that can be regarded as an upper bound in order to assess the optimality gap ω,and that can act as another algorithm termination.The feasible solution can be recovered by fixing either the mass flow rates or the nodal temperatures obtained from the McCormick solution and re-optimizing the operation problem with the fixed value.

    a UB:upper bound;LB:lower bound.

    3.3.Overview of the proposed method and algorithm

    Fig.1 presents an intuitive illustration of the tightening McCormick algorithm,including the piecewise McCormick technique and the bound contraction algorithm.Fig.2 illustrates the main procedure of the proposed method and algorithm.After model reformulation and convex relaxation,a convex McCormick model is formulated and provides an objective lower bound.To tighten the McCormick relaxation,stronger upper and lower bounds of the nodal temperatures are derived by partitioning the variable domain and are solved by means of the piecewise McCormick technique.Meanwhile,a feasible solution is expected to provide an upper bound and to form the stopping criterion together with the lower bound.The feasible solution can be recovered by fixing the mass flow rates and solving the resulting convex reformulation model,Eqs.(13a)–(13e),with Eqs.(9a)–(9g).The upper and lower bounds of the mass flow rates and the temperatures are successively tightened based on the bound contraction algorithm.Finally,the algorithm converges to a local optimum near the objective lower bound.

    Fig.1.Illustration of the tightening McCormick algorithm.(a)The piecewise McCormick technique;(b)the bound contraction algorithm.

    Fig.2.The implementation flowchart for the proposed method and algorithm.

    4.Case studies

    In this section,we compare the performance of several models,as shown in Table 2.

    We test these models in two integrated heat and electricity systems.One is a small-scale system with a 6-bus EPS and an 8-node DHS,as shown in Fig.3.The 6-bus EPS is connected to the outer grid at bus 1.The 8-node DHS is modified based on a real DHS in Jilin Province,China[7].There are two TUs,three electric loads,a CHP unit that connects the EPS and DHS,an HB,and three heat loads.The total heat load and electric load are shown in Fig.4.The parameters are set to ε1= 0.02,κ = 0.01,and δ = 0.0001.The other is a large-scale system consisting of a modified 118-bus EPS and a 33-node DHS in Barry Island[4].The 118-bus EPS exchanges power with an outer grid at bus 69.The 33-node DHS disconnects branch number 25–28 to form a radial network.There are two CHP units and an HB in the DHS.Detailed data for the two integrated heat and electricity systems can be found in Ref.[34].

    Fig.4.The daily total heat and electricity profile.

    4.1.Optimality

    Table 3 and Table 4 show the optimal dispatch comparison in two cases,respectively.The base and reformulated models solved by the bilinear solver in Gurobi 9.0 provide the global optimal results.It is notable that there is almost no gap between those two models,which demonstrates that the reformulation preserves the same solution structure.As expected,the bilinearremoved model gives the lowest bound of the optimal solution,while the McCormick model provides a tighter lower bound.The base model solved by the local solver IPOPT presents an upper bound of the optimal solution.The tightening McCormick model has a small value gap and a fast solution time.When the network becomes large,the base model with the bilinear solver in Gurobi is time-consuming to solve.Since the problem is highly nonconvex,IPOPT even fails to converge in the large-scale case.However,with the convexified model,the tightening McCormick algorithm can be applied to the large-scale integrated system optimization with relatively small errors and is efficiently solved by off-the-shelf solvers.It is worth noting that,except for the preprocessed bound-strengthening process conducted by the piecewise McCormick technique,the remaining procedure in the tightening McCormick algorithm is actually applied on a convex model.In all,the tightening McCormick model performs well in terms of both solution accuracy and computational efficiency.Moreover,compared with using a global search,like branch and bound in Gurobi 9.0,the tightening McCormick model is a convex model,and it is easy to derive shadow prices for further economic analysis.

    From the perspective of the regulation methods,it can be observed that the regulation with variable mass flow rates has lower costs than that with constant mass flow rates,as shown in Table 3 and Table 4,since the variable mass flow rate case results in more flexibility in both the mass flow rates and the nodal temperatures.A detailed day-ahead operation is demonstrated in Fig.5,including the power output of each production unit,power selling to the grid,and the mass flow rates of typical pipelines.

    4.2.Feasibility

    Fig.6 shows violations of the constraints(Eq.(9f))in two cases.Violations are represented by the errors of the bilinear constraints;that is,Although Table 3 and Table 4 show that the bilinear-removed and McCormick models provide lower costs,it can be observed that their violations,as presented in Fig.6,are huge.This is unacceptable when applied to real-world operations.However,the tightening McCormick model gives relatively small violation errors.The maximum errors are no more than 0.358%,and the average errors are less than 0.133%in both cases.Detailed error data is provided in Table 5 and Table 6.

    Table 2 Model settings for comparison.

    Table 3 Optimal dispatch comparison in the small-scale case.

    Table 4 Optimal dispatch comparison in the large-scale case.

    Table 5 Errors of the bilinear constraints(Eq.(9f))during 24 h in the small-scale case.

    Table 6 Errors of the bilinear constraints(Eq.(9f))during 24 h in the large-scale case.

    4.3.Sensitivity analysis

    From Fig.7(a),it is obvious that a larger partition number presents tighter upper and lower bounds at the starting points for bound contraction iterations;thus,the error with caseS= 3 atiteration 1 is the smallest among all partition number cases.In the caseS= 1,the McCormick model is used directly without the piecewise process.The larger the partition number is,the more the binaries are required,and the more expensive the computation time will be.Choosing a partition number of 3 or 2 can present a satisfactory performance.Fig.7(b)demonstrates that a smaller initial value and a larger step size result in faster convergence.However,the initial value should not be too small and the step size should not be too large in case there is no feasible region.

    Fig.3.An integrated heat and electricity system with a 6-bus(i1–i6)EPS and an 8-node(n1–n8)DHS.Pd1–Pd3 refer to the electric demand 1–3;Hd1–Hd3 refer to the heating demand 1–3.

    Fig.5.The day-ahead operation of the small-scale case with different regulation methods.(a)Variable mass flow rates;(b)constant mass flow rates.

    Fig.6.Violations of constraints(Eq.(9f))at hour 1.(a)The small-scale case;(b)the large-scale case.

    5.Conclusions

    In this paper,we proposed a convex operation model of integrated heat and electricity systems.To reduce the bilinear terms coming from the quality–quantity-regulated heat flow,we first reformulated the model through a variable substitution and equivalent transformation.Next,we relaxed the remaining bilinear terms by means of McCormack envelopes.To further reduce the McCormick volume,a piecewise McCormick technique was presented to strengthen the bounds of the decision variables in bilinear terms.To ensure the feasibility of the results,a bound contraction algorithm was presented to improve the bounds of the McCormick model and obtain a nearby feasible solution.Case studies showed that the tightening McCormick method quickly solves the problem with an acceptable feasibility check and optimality.Meanwhile,with convexity,the tightening McCormick method is promising for large-scale integrated heat and electricity system optimization and can allow economic analysis using shadow prices.It is notable that the performance of the tightening McCormick algorithm greatly depends on the hyperparameter sequence chosen in the bound contraction algorithm.Further studies could focus on finding an efficient sequence and performing parameter analysis.

    In this paper,we applied the proposed model and algorithm to the energy production operation model at daily or weekly time scales with steady-state DHS.However,the proposed methods can also be extended to the intra-day or real-time economic dispatch of integrated systems,if the newly added thermal dynamics with transmission delay can be formulated as convex or bilinear constraints.The bilinear form of thermal dynamics has been realized in the so-called water mass method[15,17].Future work could consider the thermal dynamics and conduct real-time economic dispatch analysis with the proposed methods.

    Fig.7.The convergence performance in the small-scale case.(a)Partition number;(b)initial point and step size of the tightening ratio sequence. S refers to the partition number in the piecewise McCormick technique,and ε1 and κ are the initial point and step size of the tightening ratio sequence,respectively.Errors denote the mean value of the violation of the bilinear constraints in Eq.(9f).

    Acknowledgement

    This work was supported by the Science and Technology Program of State Grid Corporation of China(522300190008).

    Compliance with ethics guidelines

    Lirong Deng,Hongbin Sun,Baoju Li,Yong Sun,Tianshu Yang,and Xuan Zhang declare that they have no conflict of interest or financial conflicts to disclose.

    亚洲国产色片| 国产av不卡久久| 久久精品夜夜夜夜夜久久蜜豆| 国产女主播在线喷水免费视频网站 | 国产精品伦人一区二区| 麻豆成人午夜福利视频| 久久久色成人| 国产在视频线精品| 69av精品久久久久久| 午夜福利视频精品| 一夜夜www| 麻豆国产97在线/欧美| 黑人高潮一二区| 亚洲欧美成人综合另类久久久| 国产av在哪里看| 一区二区三区高清视频在线| 99视频精品全部免费 在线| 夫妻性生交免费视频一级片| 大香蕉久久网| 91精品一卡2卡3卡4卡| 国产伦精品一区二区三区四那| 大又大粗又爽又黄少妇毛片口| 亚洲图色成人| 亚洲av不卡在线观看| 免费看日本二区| 亚洲精品456在线播放app| 免费黄色在线免费观看| 国产黄色视频一区二区在线观看| 日本一二三区视频观看| 大片免费播放器 马上看| 久久6这里有精品| 亚洲欧美日韩卡通动漫| 校园人妻丝袜中文字幕| 激情 狠狠 欧美| 嫩草影院精品99| 淫秽高清视频在线观看| 一级片'在线观看视频| 国产欧美另类精品又又久久亚洲欧美| 九九久久精品国产亚洲av麻豆| 美女主播在线视频| 国产午夜精品久久久久久一区二区三区| 免费看a级黄色片| 日韩av在线免费看完整版不卡| 中文天堂在线官网| av网站免费在线观看视频 | 亚洲国产av新网站| 欧美bdsm另类| 国产一区二区在线观看日韩| 哪个播放器可以免费观看大片| 国产午夜精品久久久久久一区二区三区| 日产精品乱码卡一卡2卡三| 一区二区三区高清视频在线| 我的老师免费观看完整版| 国产免费福利视频在线观看| 寂寞人妻少妇视频99o| 全区人妻精品视频| av天堂中文字幕网| 精品国产一区二区三区久久久樱花 | 特级一级黄色大片| 亚洲欧美日韩东京热| 熟妇人妻久久中文字幕3abv| 91狼人影院| 国产成人一区二区在线| 久久久久免费精品人妻一区二区| 少妇的逼水好多| 国产亚洲精品av在线| 18+在线观看网站| 亚洲在久久综合| 特级一级黄色大片| 三级经典国产精品| 在线a可以看的网站| 大陆偷拍与自拍| 在线a可以看的网站| av专区在线播放| 尾随美女入室| 内射极品少妇av片p| 成人亚洲欧美一区二区av| 国产一级毛片七仙女欲春2| 熟女电影av网| 一级毛片黄色毛片免费观看视频| 国产有黄有色有爽视频| 99视频精品全部免费 在线| 身体一侧抽搐| av网站免费在线观看视频 | 又粗又硬又长又爽又黄的视频| 一边亲一边摸免费视频| av一本久久久久| av在线天堂中文字幕| 嫩草影院新地址| 女人十人毛片免费观看3o分钟| 嫩草影院入口| 最近手机中文字幕大全| 日本欧美国产在线视频| 午夜爱爱视频在线播放| 国产黄色视频一区二区在线观看| 久久精品国产亚洲网站| 成人二区视频| 欧美成人午夜免费资源| 国产亚洲av片在线观看秒播厂 | 亚洲欧美成人精品一区二区| 日韩在线高清观看一区二区三区| 欧美激情国产日韩精品一区| 国产白丝娇喘喷水9色精品| 精品国产一区二区三区久久久樱花 | 国产av码专区亚洲av| av在线蜜桃| 国产高潮美女av| 国产精品.久久久| 久久久久免费精品人妻一区二区| 国产在线一区二区三区精| 国产久久久一区二区三区| 精品人妻一区二区三区麻豆| 精品久久久噜噜| 少妇的逼水好多| 亚洲在久久综合| 看十八女毛片水多多多| 99视频精品全部免费 在线| 国产免费一级a男人的天堂| 日韩人妻高清精品专区| 亚洲乱码一区二区免费版| 久久精品夜夜夜夜夜久久蜜豆| 亚洲av福利一区| 亚洲最大成人av| 一个人观看的视频www高清免费观看| 性色avwww在线观看| 国产成人aa在线观看| 亚洲一区高清亚洲精品| 成年av动漫网址| 国产大屁股一区二区在线视频| 成年女人在线观看亚洲视频 | 天堂网av新在线| 中国美白少妇内射xxxbb| 国内精品一区二区在线观看| 亚洲成人av在线免费| 国产爱豆传媒在线观看| 91aial.com中文字幕在线观看| 久久久成人免费电影| 少妇的逼好多水| 99九九线精品视频在线观看视频| 久久6这里有精品| 一本久久精品| 秋霞在线观看毛片| 欧美高清性xxxxhd video| 国产一区有黄有色的免费视频 | 边亲边吃奶的免费视频| 97人妻精品一区二区三区麻豆| 亚洲av.av天堂| or卡值多少钱| 在线观看免费高清a一片| 毛片女人毛片| 久久人人爽人人片av| 伊人久久国产一区二区| 国产午夜精品一二区理论片| 九草在线视频观看| 插逼视频在线观看| 日韩精品青青久久久久久| 久久国产乱子免费精品| 久久久久久久久大av| 国产视频首页在线观看| 免费观看无遮挡的男女| 亚洲四区av| 国产v大片淫在线免费观看| 99热这里只有是精品50| 亚洲人成网站在线观看播放| 少妇熟女aⅴ在线视频| 国产免费一级a男人的天堂| 91av网一区二区| 亚洲自偷自拍三级| 18禁在线播放成人免费| 99久久精品国产国产毛片| 亚洲成人一二三区av| videos熟女内射| 波多野结衣巨乳人妻| 国产黄片美女视频| 日韩一本色道免费dvd| 五月伊人婷婷丁香| .国产精品久久| 日产精品乱码卡一卡2卡三| 免费av毛片视频| 亚洲欧美日韩卡通动漫| 成人av在线播放网站| 中文字幕免费在线视频6| 亚洲丝袜综合中文字幕| 国产精品一区二区三区四区免费观看| 日本黄大片高清| 亚洲成人中文字幕在线播放| 久久人人爽人人爽人人片va| 日韩成人av中文字幕在线观看| 亚洲精品国产av成人精品| 天堂俺去俺来也www色官网 | 国产成人精品福利久久| 日韩av不卡免费在线播放| 久久久成人免费电影| 一级毛片电影观看| 日本一二三区视频观看| 一级av片app| 欧美97在线视频| 少妇熟女aⅴ在线视频| 国产爱豆传媒在线观看| 一级毛片 在线播放| 国产高清国产精品国产三级 | h日本视频在线播放| 成人毛片a级毛片在线播放| 亚洲美女视频黄频| or卡值多少钱| 国产精品嫩草影院av在线观看| 99久国产av精品| 中文字幕av在线有码专区| 欧美日韩一区二区视频在线观看视频在线 | 亚洲四区av| 99热这里只有精品一区| 青春草视频在线免费观看| av天堂中文字幕网| 亚洲精品成人久久久久久| av一本久久久久| 亚洲精品国产av成人精品| 中文字幕亚洲精品专区| 亚洲人成网站高清观看| 亚洲在久久综合| 男女啪啪激烈高潮av片| 国产成年人精品一区二区| 一夜夜www| 三级国产精品片| 免费观看精品视频网站| 亚洲色图av天堂| 久久久午夜欧美精品| 美女国产视频在线观看| 欧美日韩视频高清一区二区三区二| 国产片特级美女逼逼视频| 国产精品蜜桃在线观看| 亚洲欧美成人精品一区二区| 亚洲第一区二区三区不卡| 精品国产三级普通话版| 欧美性猛交╳xxx乱大交人| 亚洲av免费高清在线观看| 蜜桃久久精品国产亚洲av| 日韩成人伦理影院| 国产成人精品久久久久久| 亚洲av成人精品一二三区| 国产精品爽爽va在线观看网站| 在现免费观看毛片| 蜜桃久久精品国产亚洲av| 日本熟妇午夜| 天堂中文最新版在线下载 | 国产精品.久久久| 中文欧美无线码| 国国产精品蜜臀av免费| 国产在视频线精品| 九九爱精品视频在线观看| 日日摸夜夜添夜夜添av毛片| 午夜福利网站1000一区二区三区| 免费av毛片视频| 69av精品久久久久久| 美女高潮的动态| av专区在线播放| 国产免费福利视频在线观看| 丰满人妻一区二区三区视频av| 少妇猛男粗大的猛烈进出视频 | 欧美3d第一页| 久久精品熟女亚洲av麻豆精品 | 久久这里只有精品中国| 日本爱情动作片www.在线观看| 午夜激情福利司机影院| 免费人成在线观看视频色| 夜夜看夜夜爽夜夜摸| 色综合亚洲欧美另类图片| 国精品久久久久久国模美| 啦啦啦啦在线视频资源| 亚洲精华国产精华液的使用体验| 日本黄大片高清| 97超视频在线观看视频| av国产久精品久网站免费入址| 成年av动漫网址| 国产成人精品福利久久| 国产一区二区三区av在线| 激情 狠狠 欧美| 三级毛片av免费| 国产在线男女| 久久99精品国语久久久| 偷拍熟女少妇极品色| 国产亚洲一区二区精品| 国产伦理片在线播放av一区| 亚洲真实伦在线观看| 亚洲最大成人av| 99热全是精品| 免费少妇av软件| 人体艺术视频欧美日本| 亚洲国产精品sss在线观看| 国产视频首页在线观看| 亚洲精品第二区| 亚洲人成网站在线观看播放| 亚洲欧美日韩无卡精品| 嫩草影院精品99| 三级毛片av免费| 亚洲欧美成人精品一区二区| 男女啪啪激烈高潮av片| 久久人人爽人人爽人人片va| 在线 av 中文字幕| 日韩电影二区| 国产一级毛片七仙女欲春2| 国产成年人精品一区二区| 亚洲最大成人手机在线| 国产大屁股一区二区在线视频| 熟女电影av网| 色综合亚洲欧美另类图片| 一区二区三区四区激情视频| 欧美日韩亚洲高清精品| 天天一区二区日本电影三级| 男女下面进入的视频免费午夜| 久久97久久精品| 91精品国产九色| 久久鲁丝午夜福利片| 精品人妻视频免费看| 少妇裸体淫交视频免费看高清| 亚洲成人一二三区av| 搞女人的毛片| 久久精品国产自在天天线| 六月丁香七月| 中文在线观看免费www的网站| 插逼视频在线观看| 丰满人妻一区二区三区视频av| 亚洲精品国产av成人精品| 国产 一区 欧美 日韩| 免费观看在线日韩| 丰满人妻一区二区三区视频av| 99久久精品国产国产毛片| 国产日韩欧美在线精品| 夜夜看夜夜爽夜夜摸| av又黄又爽大尺度在线免费看| 好男人在线观看高清免费视频| 夜夜爽夜夜爽视频| 我的老师免费观看完整版| 亚洲内射少妇av| 97精品久久久久久久久久精品| 亚洲精华国产精华液的使用体验| 亚洲真实伦在线观看| 成人漫画全彩无遮挡| 久久精品久久久久久久性| 国产一区二区三区综合在线观看 | av又黄又爽大尺度在线免费看| 国产精品嫩草影院av在线观看| 黄色欧美视频在线观看| 国产精品一区二区性色av| 国产精品三级大全| 午夜激情福利司机影院| 91狼人影院| 又大又黄又爽视频免费| 亚洲高清免费不卡视频| 亚洲精品国产av蜜桃| 久久久久久久久久久免费av| 国产伦精品一区二区三区视频9| 一级毛片我不卡| 亚洲人成网站高清观看| 国产爱豆传媒在线观看| 最近视频中文字幕2019在线8| 国产一级毛片在线| 亚洲经典国产精华液单| 免费av毛片视频| 亚洲四区av| 久久久精品免费免费高清| 国产亚洲5aaaaa淫片| 日日干狠狠操夜夜爽| 亚洲美女搞黄在线观看| 你懂的网址亚洲精品在线观看| xxx大片免费视频| 免费黄色在线免费观看| 精品少妇黑人巨大在线播放| 国内揄拍国产精品人妻在线| 欧美日韩国产mv在线观看视频 | 欧美精品国产亚洲| 狂野欧美白嫩少妇大欣赏| 特级一级黄色大片| 草草在线视频免费看| 日本一二三区视频观看| 成年女人在线观看亚洲视频 | 日韩av在线免费看完整版不卡| 国产黄色小视频在线观看| 尤物成人国产欧美一区二区三区| 日本黄大片高清| 国产高清国产精品国产三级 | 街头女战士在线观看网站| 男女视频在线观看网站免费| 国产亚洲av片在线观看秒播厂 | 又大又黄又爽视频免费| 男女边吃奶边做爰视频| 国产精品国产三级国产专区5o| 男女那种视频在线观看| 欧美三级亚洲精品| 亚洲欧洲日产国产| 又爽又黄无遮挡网站| 在线天堂最新版资源| 淫秽高清视频在线观看| 欧美日韩亚洲高清精品| 亚洲电影在线观看av| a级毛片免费高清观看在线播放| 久久99热6这里只有精品| 大香蕉久久网| 国产国拍精品亚洲av在线观看| 青春草国产在线视频| 国产高清不卡午夜福利| 国国产精品蜜臀av免费| 特大巨黑吊av在线直播| 在线天堂最新版资源| 激情五月婷婷亚洲| 午夜福利成人在线免费观看| 免费观看av网站的网址| 中国国产av一级| 国产久久久一区二区三区| 纵有疾风起免费观看全集完整版 | 午夜爱爱视频在线播放| 久久鲁丝午夜福利片| 国产一区亚洲一区在线观看| 久久人人爽人人片av| 成人亚洲精品一区在线观看 | 日韩av在线免费看完整版不卡| 中国国产av一级| 色尼玛亚洲综合影院| av在线蜜桃| 中国国产av一级| av黄色大香蕉| 日日啪夜夜撸| 美女内射精品一级片tv| 少妇丰满av| 老司机影院成人| 少妇熟女aⅴ在线视频| 最近最新中文字幕免费大全7| 免费黄色在线免费观看| 国产精品一区二区三区四区免费观看| 亚洲18禁久久av| 欧美3d第一页| 精品国产露脸久久av麻豆 | 一本一本综合久久| 久久亚洲国产成人精品v| 欧美日韩精品成人综合77777| 午夜福利视频1000在线观看| 国产精品99久久久久久久久| 亚洲av.av天堂| 青春草国产在线视频| 夜夜看夜夜爽夜夜摸| 噜噜噜噜噜久久久久久91| 日本三级黄在线观看| 日韩av在线免费看完整版不卡| 久久韩国三级中文字幕| 亚洲三级黄色毛片| 搡女人真爽免费视频火全软件| 人妻少妇偷人精品九色| 天天一区二区日本电影三级| 一边亲一边摸免费视频| 亚洲久久久久久中文字幕| 亚洲激情五月婷婷啪啪| 国产成人a∨麻豆精品| 国产亚洲精品av在线| 高清日韩中文字幕在线| 天堂俺去俺来也www色官网 | 草草在线视频免费看| 韩国av在线不卡| 精品国产露脸久久av麻豆 | 联通29元200g的流量卡| 亚洲真实伦在线观看| 美女cb高潮喷水在线观看| 伊人久久精品亚洲午夜| 一级片'在线观看视频| 日本三级黄在线观看| 婷婷色麻豆天堂久久| 国产在视频线精品| 一个人看视频在线观看www免费| 丰满人妻一区二区三区视频av| 成人二区视频| 边亲边吃奶的免费视频| 国产黄片视频在线免费观看| 大又大粗又爽又黄少妇毛片口| 日本午夜av视频| 午夜激情福利司机影院| 美女脱内裤让男人舔精品视频| 日韩成人av中文字幕在线观看| 国产免费又黄又爽又色| eeuss影院久久| 国产黄色视频一区二区在线观看| 国产男女超爽视频在线观看| 国产综合精华液| 久久这里有精品视频免费| 精品久久久久久电影网| 国产亚洲一区二区精品| 两个人的视频大全免费| 一个人观看的视频www高清免费观看| 国产精品人妻久久久久久| 日韩三级伦理在线观看| 午夜福利成人在线免费观看| 日韩精品青青久久久久久| 久久久久九九精品影院| 久久鲁丝午夜福利片| 一级毛片电影观看| 国产色爽女视频免费观看| 国产不卡一卡二| 少妇人妻精品综合一区二区| 国产熟女欧美一区二区| 女人十人毛片免费观看3o分钟| 午夜视频国产福利| 高清av免费在线| 乱码一卡2卡4卡精品| 免费观看a级毛片全部| 国产精品国产三级国产专区5o| 一级毛片黄色毛片免费观看视频| 91精品伊人久久大香线蕉| 国产精品不卡视频一区二区| 国产男人的电影天堂91| 丝瓜视频免费看黄片| 街头女战士在线观看网站| 亚洲欧美中文字幕日韩二区| 成人亚洲欧美一区二区av| 最近中文字幕高清免费大全6| 久久精品国产亚洲av天美| 99视频精品全部免费 在线| 啦啦啦中文免费视频观看日本| 中文字幕亚洲精品专区| 高清午夜精品一区二区三区| 国产成人a∨麻豆精品| 亚洲在久久综合| 国产亚洲最大av| 青春草视频在线免费观看| 七月丁香在线播放| 欧美极品一区二区三区四区| 日产精品乱码卡一卡2卡三| 亚洲精品成人av观看孕妇| 日韩欧美 国产精品| 能在线免费看毛片的网站| 亚洲色图av天堂| 搡老妇女老女人老熟妇| 免费av观看视频| 日韩av免费高清视频| 日韩视频在线欧美| 亚洲av成人精品一二三区| 亚洲精品456在线播放app| 国产黄频视频在线观看| 免费在线观看成人毛片| 两个人的视频大全免费| 亚洲欧美清纯卡通| 人体艺术视频欧美日本| 六月丁香七月| 国产白丝娇喘喷水9色精品| 亚洲av一区综合| 久久久欧美国产精品| 精品久久国产蜜桃| 国产高清三级在线| 国产精品一及| 乱码一卡2卡4卡精品| 精品久久久久久电影网| 亚洲不卡免费看| 97人妻精品一区二区三区麻豆| 国产精品一二三区在线看| 免费高清在线观看视频在线观看| 午夜精品在线福利| 夫妻性生交免费视频一级片| 久久韩国三级中文字幕| 久久鲁丝午夜福利片| 国产高清有码在线观看视频| 亚洲,欧美,日韩| 欧美日韩视频高清一区二区三区二| 日韩欧美一区视频在线观看 | 日本猛色少妇xxxxx猛交久久| 啦啦啦韩国在线观看视频| 亚洲一级一片aⅴ在线观看| 十八禁国产超污无遮挡网站| 丰满人妻一区二区三区视频av| 国产高清三级在线| 日韩欧美三级三区| 最近最新中文字幕大全电影3| 99久久中文字幕三级久久日本| 欧美成人午夜免费资源| 日韩不卡一区二区三区视频在线| 99热这里只有精品一区| 成人国产麻豆网| 久久精品国产亚洲av天美| 国内揄拍国产精品人妻在线| 熟妇人妻不卡中文字幕| 免费播放大片免费观看视频在线观看| 成年av动漫网址| 国产伦精品一区二区三区视频9| 亚洲图色成人| 色综合亚洲欧美另类图片| 又爽又黄a免费视频| 久久久久久伊人网av| videos熟女内射| 黑人高潮一二区| 亚洲三级黄色毛片| 精品久久久久久久久av| 成人午夜精彩视频在线观看| 一级毛片电影观看| kizo精华| 国产综合精华液| 国产成人freesex在线| 亚洲欧美日韩卡通动漫| 三级男女做爰猛烈吃奶摸视频| 性色avwww在线观看| av.在线天堂| 国产淫片久久久久久久久| 日韩人妻高清精品专区| 秋霞在线观看毛片| 黄片无遮挡物在线观看| 国产精品伦人一区二区| 免费播放大片免费观看视频在线观看| 亚洲久久久久久中文字幕| 久久这里只有精品中国| 观看免费一级毛片| 午夜激情久久久久久久| 80岁老熟妇乱子伦牲交| 中国国产av一级| 日韩国内少妇激情av| 久久午夜福利片| 肉色欧美久久久久久久蜜桃 | 一区二区三区乱码不卡18| 亚洲精品视频女| 亚洲人与动物交配视频| 国产av在哪里看| 国产精品蜜桃在线观看|