A spatial oligopolistic electricity model under uncertain demands

This paper mainly investigates a single-period spatial oligopolistic model of electricity generators under uncertain demands. We maximize the β\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta $$\end{document}-optimistic value of each generator’s uncertain profit function in this game. Emphasis is first put on computing the β\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta $$\end{document}-optimistic value of the uncertain profit function through the inverse uncertainty distribution of the uncertain variable in the inverse demand function. We later focus on transforming the single-period spatial oligopolistic game, which is a generalized Nash equilibrium (short for GNE ) problem, into a variational inequality problem. Then the model is applied to simulate the transmission price of the electricity system of America, which reveals that our model is effective.


Introduction
The electricity available to consumers requires three basic steps: generation, transmission, and distribution. Those involve two markets: energy market and transmission service market. Both markets are operated under a monopoly of the past. To avoid the dependence on a single supply region, which brings the ineffectiveness of electricity market, the electricity systems currently are restructuring or going to restructure in several countries. This process has different paradigms, leading to different ways to develop electricity models for the power market. One way which attracts considerable attention is modeling through Nash equilibria. As far as we are concerned, there are mainly three approaches to model the spatial oligopolistic electricity market by Nash equilibria. One is the supply function equilibrium model, see references (Bolle 1992;Green and Newbery 1992;Klem-perer and Meyer 1989). Another is the Bertrand equilibrium model, please refer to Hobbs and Kelly (1992). The third is the Cournot equilibrium model, see (Borenstein and Bushnell 1996;Oren 1997;Schmalensee and Golub 1984).
Historically speaking, researchers studied either the longrun equilibria or the short-term equilibria of those models. For instance, Wei and Smeer (1999) looked at a long-term spatial oligopolistic electricity equilibrium model. Fehr and Harbord (1993) investigated a short-run supply function equilibrium model of the electricity in the UK. Smeers and Wei (1997) discussed different short-run equilibria on the electricity power market. All work we mentioned above only focused on a static situation. However, the electricity market involves some indeterminancies in reality. For instance, not only some unexpected weather condition like hurricanes, tornadoes, etc., but also some natural disasters like floods, seaquakes, earthquakes, etc., even some emergencies like train crashes and rumors will influence the power demands in a region. To cope with the indeterminancies, Pineau and Murto (2003) considered a dynamic stochastic oligopolistic model of multi-market segments (the demand growth rate of the supply was modeled as a stochastic variable) and applied it to the Finnish electricity market.
However, the indeterminancy lying the nature and society does not confine to the randomness. Beside randomness, there is another type of indeterminancy-uncertainties. Here "Uncertainties", as a jargon, come from the uncertainty theory, which is a branch of mathematics to deal with indeterminancies and concerned with the analysis of belief degree (Please refer to Liu 2019). The belief degree appeared here describes to what degree someone believes an event will happen under the case that there is no sample or very few samples to derive from what is studied. For example, it is difficult to know the exact oil reserves in a region. For an oil company, there is no choice but to invite some experts to estimate and give their belief degrees to the quantities of the oil underground. Then the oil company will decide whether to bid the oilfield according to the analysis of experts' belief degrees.
Inspired by the aforementioned studies, this paper mainly investigates a single-period spatial oligopolistic electricity model, regarding the electricity demand as an uncertain variable. The case we study only involves the generation and the transmission of electricity. All producers compete for supplying the electricity to consumers in the Cournot manner, and the transmission service is provided by a transmission company. All producers aim at maximizing their profits under some uncertain factors. Profit, in this situation, has no conceptual meaning. We maximize uncertain profits of producers in the sense of optimistic values. That is the bigger the optimistic value of a producer's uncertain profit is, the larger profit a producer earns.
The present paper is constructed as follows. Some basic concepts and lemmas in uncertainty theory are covered in Sect. 2. A single-period spatial oligopolistic electricity model with uncertain demands is given in Sect. 3. The description of the single-period spatial oligopolistic electricity model with uncertain demands is described in Sect. 3.1 while its formulation into a variational inequality is presented in Sect. 3.2. An illustrative example is constructed in Sect. 4. We conclude our work in Sect. 5. Appendix discusses the β−optimistic/pessimistic value of a crisp number, which is a special uncertain variable in uncertainty theory.

Preliminary
In this section, we recall some basics in uncertainty theory, which is a new branch of mathematics founded by Liu (2007) in 2007.
Let Γ be a nonempty set and let L be a σ −algebra on it. We call the ordered pair (Γ , L) a measurable space. The members of L are measurable sets in Γ , which are also renamed as events in uncertainty theory. Liu (2004) defined an uncertain measure on (Γ , L) to describe the belief degree with which a person believes an event Λ may happen. It is a set function M defined on the measurable space (Γ , L) which satisfies the following three axioms: Axiom 1 (Normality Axiom) M(Γ ) = 1 for the universal set Γ .
Axiom 2 (Duality Axiom) M{Λ} + M{Λ c } = 1 for any event Λ. Axiom 3 (Subadditivity Axiom) For every countable sequence of events Λ 1 , Λ 2 , . . . , we have The triplet (Γ , L, M) is called an uncertainty space. It is one vital concept in uncertainty theory. Another one is the uncertain variable used to represent the quantity of uncertainty. An uncertain variable is a function ξ from an uncertainty space (Γ , L, M) to the set of real numbers such that {ξ ∈ B} = {γ ∈ Γ | ξ(γ ) ∈ B} is an event for any Borel set B. Uncertain variables ξ 1 , ξ 2 , . . . , ξ n are said to be inde- The function Φ(t) is said to be regular if it is continuous and strictly increasing with respect to t at which 0 < Φ(t) < 1, and lim We now present what we need to use in the sequel.
Definition 1 (Liu 2014) An uncertain variable ξ is called linear if it has a linear uncertainty distribution denoted by L(a, b), where a and b are real numbers with a < b.
Definition 2 (Liu 2004) An uncertain variable ξ is called normal if it has a normal uncertainty distribution denoted by N (e, σ ), where e and σ are real numbers with σ > 0.
Remark 1 (1) An uncertain variable ξ with linear uncertainty distribution L(a, b) has an inverse uncertainty distribution (2) An uncertain variable ξ with normal uncertainty distribution N (e, σ ) has an inverse uncertainty distribution Definition 3 (Liu 2007) Let ξ be an uncertain variable, and β ∈ (0, 1]. Then The following conclusions can be drawn from definitions of linear uncertainty distribution, normal uncertainty distribution, β-optimistic/pessimistic values and Remark 1.
Furthermore, β−optimistic/pessimistic values share the properties as follows.
Lemma 4 (Liu 2014) Let ξ be an uncertain variable with uncertainty distribution Φ(t), and let g(y) be a strictly decreasing function with respect to y. Then g(ξ ) has an uncertainty distribution where g −1 (t) is the inverse function of g(y).

Problem description
We describe in this part a single-period spatial oligopolistic electricity model under uncertain demands. The generators are assumed to behave in the Cournot manner. The relevant notions are introduced first as a preliminary to define strategies and equilibria of this game. s i j : Power supplied by firm i to region j. s j = i∈I s i j : Total power supplied in region j. p j (s j , ξ j ): Uncertain inverse demand function of energy in region j, where the uncertain variable ξ j is an uncertain factor. q k : Electricity produced by plant k ∈ G i j . C k (q k ): Production cost function of plant k. q k : Generation capacity at plant k. f i, ( j, j ) : Power generated by firm i and sent from node j through transmission capacity ( j, j ) to node j .
Vector of prices for transmission services, where A is the set of electric transmission capacities in the network.
Thus, the strategy set for each producer can now be introduced. Assume that all producers play the game in a noncooperative manner. Both of the strategy set and the payoff function of a producer depend on the other produc-ers' strategies. Consequently, the equilibrium of this spatial oligopolistic model is a generalized Nash equilibrium.
For simplicity, we denote as producer i's decision vector. Additionally, let x −i = (x i : i = i) be the decision vector of all firms but i. The constraints that must be satisfied by x i are as follows: 1. The energy flow has a bound. That is the flow in each link ( j, j ) ∈ A is less than or equal to the available transmission capacity. 2. The power output of each plant is limited. That is to say, the generation level is less than or equal to the available generation capacity for each plant. 3. The power in any region keeps balance. That is the flow, the generation and the supply of the power are balanced at each node.
When fixed x −i , we denote K i (x −i ) as firm i's feasible set satisfied all above-mentioned constraints. That is ( j, j ) . Note that firm i's feasible set K i (x −i ) depends on other firms' decision variables through the following capacity constraint: Let x = (x i : i ∈ I ) be the vector of all firms' decision variables. Then firm i's uncertain profit function equals its uncertain revenue minus its generation cost and transmission cost: Note that θ i (x, ξ i ) depends on other firms' strategies only through the demand variables s i j with i = i in s j .
In this paper, we attempt to use the β−optimistic value criterion to study the equilibrium in the electricity power market. This leads to solve the following optimization problem: fixed x −i , which is composed of other firms' decision variables, and the transmission price vector ρ, solve the problem for all i ∈ I , where β is a given confidence level in (0, 1). For simplicity, we denote problem (10) The overall equilibrium of this model is to solve each β − O PT (x −i , ρ) for all i ∈ I , which is further denoted as problem β − O PT . When solving the uncertain profit maximization problem under the criterion of β−optimistic value, we treat other firms' decision variables and transmission prices as endogenous. Now we present the formal definition of the equilibrium of this spatial oligopolistic electricity market.

Definition 4
The equilibrium of the spatial oligopolistic electricity market with the Cournot generators under uncertain demands by the criterion of β−optimistic value is to find a point and a set of transmission prices In other words, the overall equilibrium problem of this spatial oligopolistic electricity market is to seek a vector Obviously, producer i's strategy is dependent on other producers' through constraint (8). It is this dependence that has the implication for the interpretation of the model in terms of game theory. As mentioned previously, it also implies that the equilibrium introduced in Definition 4 is a generalized Nash equilibrium.

Problem formulation
In this section, we formulate the model of β − O PT as a single variational inequality problem. To this end, we formulate the optimization problem β − O PT (x −i , ρ) first. Note that it is naturally a quasi-variational inequality problem because of the dependence of the set K i (x −i ) on x −i for any i. The special structure of the producers' strategy sets leads to reformulate the problem β−O PT (x −i , ρ) into a variational inequality problem. The whole process consists of two phases. One is recasting the objective function θ i (x, ξ i ) sup (β); the other is reformulating the feasible set K i (x −i ). We deal with them one by one.
We handle θ i (x, ξ i ) sup (β) first. Recall that the uncertain profit for firm i is ( j, j ) . From the point of uncertainty theory, we need the following assumptions on uncertain variables.

Assumption 1
The uncertain variable ξ j , j ∈ J in the uncertain inverse demand function p j (s j , ξ j ) has a regular uncertainty distribution Φ j and the function p j (s j , ξ j ) is decreasing on ξ j for each j ∈ J . Furthermore, all uncertain variables ξ j , j ∈ J are independent of each other.
Theorem 1 Suppose that ξ j and p j (s j , ξ j ) satisfy Assumption 1, then we have Proof Since the β-optimistic value of ξ is Φ −1 (1−β) when ξ has an inverse uncertainty distribution Φ, it is easy to obtain Theorem 1 with the help of Lemma 4.
2) We can also derive both the β−optimistic and β−pessimistic value of θ i (x, ξ i ) when p j (s j , ξ j ) is increasing on ξ j for each j ∈ J similarly.
For convenience, when given β ∈ (0, 1), we denote Up till now, we had formulated the objective function in the problem β − O PT (x −i , ρ) to be a function which only involves the decision variable x and the confidence level β. Before we push our discussion further, in conformity with the generally accepted economic behavior, we present some assumptions as follows: has the following characteristics: 1-1: where S j is a constant. 1-2:p j (s j , β) is strictly decreasing and continuously twice differentiable with respect to s i j . 2) For all j ∈ J , β ∈ (0, 1), the associated revenue function s i jp j (s j , β) is strictly concave on s i j , where s j = i∈I s i j .
Assumption 3 For all i ∈ I , j ∈ J , k ∈ G i j , each cost function C k (q k ) is convex and twice continuously differentiable with respect to q k . Murphy et al. (1982) had shown that Assumptions 2 and 3 imply the following result.
Lemma 5 When given β ∈ (0, 1), s i jp j (s i j +S, β) is concave on s i j for 0 S S j − s i j and 0 S .
As a result of all Assumptions and Lemma 5, when given β ∈ (0, 1), the objective functionθ i (x, β) is continuously differentiable and concave on x.
We next turn to formulate the feasible set K i (x −i ) of firm i. To this end, let Note thatK i is independent of x −i , that is,K i only involves the decision variable of firm i. Hence, Thus,K i (x −i ) is a mapping of x −i , which is the constraint imposed on producer i's strategy through other producers' strategies. Let where M R i j (s i j , β) = ∂θ i (x, β) ∂s i j denotes firm i's marginal return in region j with a given confidence level β ∈ (0, 1), MC k (q k ) = dC k q k dq k is its marginal cost at generation plant k. Here we recall that given a closed convex set K ∈ R n and a continuous function F : K → R n , the variational inequality defined by K and F (which is denote by V I (K , F)) means to find an x * ∈ K such that Now we go back to our problem.
We reveal the relation between the model β − O PT and the V I (K , F) in the theorem followed. And under what conditions can guaranteed the existence of a solution to the V I (K , F) are also presented.
(2) Suppose that (x * , ρ * ) is an equilibrium of the model β − O P T and let λ * i,( j, j ) be a dual vector of constraint Proof The proof of this theorem is essentially the same as its counterpart in Wei and Smeer'work (1999) except that there is a parameter β in the objective function, so it is omitted here.

Remark 3
The model β−O PT is equivalent to the V I (K , F) when the transmission price ρ( f * ) is such that at all equilibria, λ * i, ( j, j ) = 0 for all i, ( j, j ) from Theorem 2.

An application
To demonstrate what we discussed in the previous section, we present an application in this part.

Market assumptions
The market investigation here consists of two investor-owned electric utilities that supply electricity to two sectors: residential sector and commercial sector. Note that in order to correspond with the data we obtained, the electric utility will be used to take place of the power producer/firm in this subsection and the followed one. Figure 1 describes the network under investigation. The data 1 is from US Energy Information Administration mainly for the year 2015.
A linear uncertain inverse demand function is calibrated for each consumption sector. This can be done by computing the price elasticity at the reference prices and considering some indeterminant factors. The reference prices, in 2015, for the residential sector and commercial sector, respectively, are 0.1265 and 0.1064 USD/KWh. The total net generation by nuclear source and thermal source is 205,174,5 GWh. It is divided between the residential and commercial sectors with the ration 48/52 in this paper. Suppose that one electric utility has only the thermal power plants while the other owns only nuclear power plants. The average power plant operating expense for a thermal power plant is 25.71 while the average power plant operating expense for a nuclear power plant is 37.26 Mills/KWh, respectively. Note that 1Mill = 1 1000 USD. In this application, we assume that the total supply of the power is synonymous to the total production as well as the total flow. Denote the two electric utilities as U 1 and U 2 . The decision variable of U i is x i , which is further composed of x i1 and x i2 (i = 1, 2). Note that x i j represents the power supplied by the utility U i (i = 1, 2) to the sector j( j = 1, 2).

Computation and analysis
The linear uncertain inverse demand functions are computed from the price elasticity between the data of sales and the average price of electricity to ultimate customers for 2014-2015 together with uncertain variables. Here, the uncertain variable depicts some uncertain factors in the electricity power market. Those uncertainties are caused by a certain natural disasters such as hurricanes, tornadoes, seaquakes, earthquakes, etc., and the uncertain demands of consumption sectors.
The linear uncertain inverse demand function of the electricity power at sectors j 1 , j 2 are, respectively, represented by p 1 (x, ξ 1 ) = 713,079.9804 − 0.4178(x 11 + x 21 + ξ 1 ), is a linear uncertain variable with the uncertainty distribution L(−100,100), while ξ 2 is a normal uncertain variable with uncertainty distribution N (100, 10). Moreover, ξ 1 and ξ 2 are independent of each other. Denote ρ = (ρ 1 , ρ 2 ) T as the transmission price vector, where ρ 1 and ρ 2 are the unit transmission prices delivered to two sectors. For convenience, the operating expense of the production is scaled up as USD/GWh, the consumption of the electricity power is weighted as GWh. Given the decision variable x −1 = x 2 and the transmission price vector ρ, the uncertain profit function of utility U 1 is The corresponding constraint set is Then the problem β − O PT (x −1 , ρ) with β = 0.75 for utility U 1 is maximize θ 1 (x, ξ ) sup (0.75) subject to x 1 ∈ K 1 (x −1 ).
From the inverse uncertainty distributions of the linear and the normal uncertain variables, together with Theorem 1, it leads to the following result.
In other words, it is ∇θ 1 (x * , 0.75) T (x * 1 − x 1 ) + ∇θ 2 (x * , 0.75) T (x * 2 − x 2 ) 0, for all x ∈ K , is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/.
It means that for any r c, the expression (22) is true, so the supremum of r is c. That is c sup (β) = c.
Hence c sup (β) has no meaning for β ∈ (0, 1] in this case. Thus, we have Similarly, we also have