key: cord-0490314-398c1e1p authors: Chen, Xiaojun; Shi, Yun; Wang, Xiaozhou title: Equilibrium Oil Market Share under the COVID-19 Pandemic date: 2020-07-30 journal: nan DOI: nan sha: 92e3480b92e5f331bbc37f538c289b194fb590a4 doc_id: 490314 cord_uid: 398c1e1p Equilibrium models for energy markets under uncertain demand and supply have attracted considerable attentions. This paper focuses on modelling crude oil market share under the COVID-19 pandemic using two-stage stochastic equilibrium. We describe the uncertainties in the demand and supply by random variables and provide two types of production decisions (here-and-now and wait-and-see). The here-and-now decision in the first stage does not depend on the outcome of random events to be revealed in the future and the wait-and-see decision in the second stage is allowed to depend on the random events in the future and adjust the feasibility of the here-and-now decision in rare unexpected scenarios such as those observed during the COVID-19 pandemic. We develop a fast algorithm to find a solution of the two-stage stochastic equilibrium. We show the robustness of the two-stage stochastic equilibrium model for forecasting the oil market share using the real market data from January 2019 to May 2020. In this paper, we consider the two-stage stochastic equilibrium model under uncertain demand and supply for forecasting the market share in an oligopoly market. This model is developed based on noncooperative game theory which allows us to investigate how the market share of one product (oil, electricity, steel, etc.) depends on the strategies of a few agents whose decisions against and affect each other. A Nash-equilibrium for oligopolies is expected, in which costs and benefits are balanced so that no agent can gain by changing only their own strategies. The uncertainties in the future events are represented by random variables in the model. The COVID-19 pandemic has significant impact on the global energy industry and arouse great challenging in maintaining market stability. Whether we can expect a Nash-equilibrium for crude oil market share during the COVID-19 pandemic is an interesting question. In traditional oil industry, each producer has normally two strategies for maximizing its profits. One is to limit supply which leads to high oil price and this approach allows high-cost producers to remain profitable. The other is to drive up production and squeeze out high-cost producers, e.g., US shale oil, under expected low price market environment. It has been shown that under certain conditions, one strategy can be more appropriate than the other [2] . However, the change in production had inconsistent market responses due to the complex and uncertain nature of oil market [24, 34] . Furthermore, it has been observed that even during the periods of violently volatile in oil market, the relative market share remained stable. See Tables 3-4 for monthly oil market share from January 2019 to May 2020. The industry of oil is no stranger to volatilities in price, supply and demand [26, 45] . Economic growth can be largely affected by oil volatility, and certain countries may be more sensitive than the others [30, 33, 42] . Most of these changes, the so-called oil shocks, were caused by the occurrence of major events of global economics and dramatic trends shifting that affects supply or/and demand sides in the crude oil market, and the situation was no different in the latest oil shock during COVID-19 pandemic. Events of the crude oil market are capable of delivering significant impact on energy markets as well as non-energy markets [23, 27, 32] . It is believed, however, despite different causes of oil shock, their consequences for economy are very similar [22] . The most apparent observation on oil uncertainty is its price, and it is naturally assumed as of certain importance to market participants. However, being an oligopoly market implies that the price is largely manipulated by major producers. In particular, these producers may choose to sell greater quantity of oil at a lower price over the strategies of limiting supply to boost price, because the former can be more profitable in long run. There are many attempts on predicting crude oil price recently using tools of machine learning and artificial intelligence [1, 3, 43] . To analyze the oil market share under uncertain supply and demand, we model the production decisions as a solution of a two-stage stochastic game. The two-stage stochastic game is suitable to reflect the complexity of the trading process while market demand needs to be predicted/estimated at the second stage to guide the production at the first stage for a certain period of time in the future. In the two-stage stochastic game, we have two types of decisions: "here-and-now" and "wait-and-see" for each agent. Endowed with market data of production, supply, demand, cost and price of the commodity, each agent chooses its here-and-now decision and wait-and-see decision with and without observation of uncertainties in the future events, respectively. For each agent, the objective of the first stage is to maximize its expected utilities under capacity constraints, and the objective of the second stage is to maximize its utilities subject to constraints endowed with its here-and-now decision and recourse term for almost every realization. The structure of the two-stage stochastic game is a natural way to represent the equilibrium state under uncertain environment. Each agent must make a here-and-now decision for a production plan before knowing the future events. The here-and-now decision will affect future revenues, cost and feasibilities over the decision horizon. The wait-and-see decision is made by taking the uncertain demand, supply and price in the future market into account. Assuming that the utility functions at both stages are concave and continuously differentiable, we can derive the equivalent two-stage stochastic equilibrium model by the first-order optimality condition, which is a special case of the two-stage stochastic variational inequalities (SVI). Variational inequalities play a central role in operations research and optimization, which model equilibrium problems in engineering and economics [16] and represent the optimality conditions of the optimization problems [39] . Stochastic variational inequalities (also called one-stage SVI ) [4, 15, 17, 20, 21, 28, 29, 37] consider a here-and-now solution in an uncertain environment using expectation of random functions. Compared with the one-stage SVI, the two-stage SVI consider a pair of a here-and-now solution and a wait-and-see solution, which have inherently dynamic components that involve uncer-tain information and whose solutions depend on the outcome of random events to be revealed in the future. In the last few years, the two-stage SVI and multi-stage SVI have attracted considerable attention [5, 7, 6, 9, 10, 25, 35, 40, 41] . Chen, Pong and Wets [5] , and Rockafellar and Wets [40] introduced the the two-stage SVI and multi-stage SVI with examples including the first-order optimality conditions for stochastic programs, Walras equilibrium problems with an incomplete financial market and stochastic Wardrop traffic flow equilibrium problems. The existence and uniqueness of solutions of the two-stage SVI, and convergence of its sample average approximation have been investigated in [6, 7] . The Progressive Hedging Algorithm (PHA) [38] for scenarios and policy aggregation in optimization under uncertainty has been extended to solve monotone two-stage SVI in [41] . Our contributions in this paper are twofold. We develop new optimization theory and algorithms for the two-stage SVI arising from an oligopoly market. Furthermore, we apply the new optimization results to show the existence of a Nash-equilibrium for oil production game in an uncertain environment. In particular, we model the game of crude oil production with a few agents as the two-stage SVI, and show the existence and uniqueness of solutions of the two-stage SVI. Moreover, we develop a fast algorithm to find a solution of the two-stage SVI, which is used to analyze and forecast oil market share under the CODVI-19 pandemic. This paper is organized as follows. In section 2, we present the two-stage SVI arising from two-stage stochastic games. In section 3, we show the existence and uniqueness of the solution of the two-stage SVI and provide perturbation error bounds for the solution. In section 4, we develop a fast algorithm called Alternation Block Algorithm (ABA) to find a solution of the two-stage stochastic equilibrium. The computational cost of the ABA is much less than the popular PHA. We show outperformance of the ABA over the PHA using large-scale numerical examples. In section 5, we apply the two-stage SVI, our new theory and algorithm to study the impact and the production responses of oil producers during COVID-19 on oil market share. Our numerical results with real market data from January 2019 to May 2020 show the efficiency of our methods for modelling oil market share. We consider an oligopolistic market where J agents compete to supply a homogeneous product noncooperatively in the future. Each agent needs to make a decision on the quantity of production based on the anticipated future market demand and supply and other agents' decision in an uncertain environment. The uncertainties are represented by a random variable ξ : Ω → Ξ ⊂ R m defined in the probability space (Ω, F, P) with a support set Ξ and the space Y of measurable functions from Ξ to R J . We define the variables and functions for each agent i, i = 1, . . . , J. x i ∈ R, the production quantity θ i : R → R, the cost function of the production y i (ξ) ∈ R, the supply quantity ϕ i : R × Ξ → R, the cost function for supplying a quantity The uncertain market supply and demand are characterized by a random total supply η(ξ) = J i=1 y i (ξ) to the market and a random inverse demand function p(η(ξ), ξ) : R × Ξ → R. Here p(η(ξ), ξ) can be regarded as the spot price of trading in the future. Each agent aims to maximize its profit and make its decision in two stages. The first stage is to make optimal decision on production quantity based on the average of random demand, and the second stage is to make optimal decision on supply quantity based on the observation of the uncertainty in the future. The first stage decision is called "here-and-now" and the second stage decision is called "wait-and-see". For each realization of random variable ξ and a nonnegative vector x ∈ R J , agent i wants to find an optimal decision y i (ξ) by solving the following problem Moreover, agent i has to make an optimal decision before knowing the future events by solving Here x −i and y −i are decision variables of the agents other than the agent i. Problem (2) is the first stage, and problem (1) is the second stage of the two-stage stochastic games [6, 7, 25, 36] . We call (x * , y * ) ∈ R J × Y an optimal solution of the two-stage stochastic games, if for i = 1, . . . , J, (x * i , y * i (ξ)) is the optimal solution of the two-stage optimization problem In this paper, we consider the case where the objective functions in (1)-(2) are quadratic concave in the following forms where c i > 0, γ(ξ) > 0, h i (ξ) > 0, a i , r i , α(ξ) and β i (ξ) are real given numbers. In such setting, the function where s i (ξ) ≥ 0 are Lagrange multipliers for the constraints y i (ξ) ≤ x i . If x i = 0, then the optimal solution of the optimization problem (4) is y i (ξ) ≡ 0 for any ξ ∈ Ξ, and the subdifferential The Karush-Kuhn-Tucker (KKT) conditions of problems (1)-(2) with the functions defined in (5)-(7) derive the following two-stage stochastic linear complementarity prob- for almost every ξ ∈ Ξ, where and e ∈ R J is the vector with all elements being 1. When the matrix C + re T is positive definite, problems (1)-(2) with the functions defined in (5)-(7) are equivalent to problem (8) in the sense that if (x * , y * (·)) is a solution of problems (1)-(2), then there is s * (·) such that (x * , y * (·), s * (·)) is a solution of (8); conversely, if (x * , y * (·), s * (·)) is a solution of (8), then (x * , y * (·)) is a solution of problems (1)- (2) . A sufficient condition for the matrix C + re T being positive definite is Condition (9) implies that C + 1 2 (re T + er T ) is a symmetric diagonally dominate matrix with positive diagonally elements and thus a positive definite matrix. Let It is easy to see that for any where Note that we cannot apply [7, Proposition 2] for the existence of solutions of the two-stage stochastic LCP (8) , since the condition that there exists a positive continuous fails for any x = 0, v = (0, s) T , s ∈ R J and s = 0. In the next section, we will establish the existence of solutions of the two-stage stochastic LCP (8) with a finite number of realizations of the random variable. We first consider the following LCP(M (ξ), q(x, ξ)) for a fixed x ≥ 0 and ξ ∈ Ξ, where M (ξ) and v(ξ) are defined in (10) and Note that the matrix M (ξ) is not a P -matrix. Thus the existence and uniqueness of the solution of the LCP(M (ξ), q(x, ξ)) and the error bound in [11, 12, 14] cannot be guaranteed for any q(x, ξ). For example, if there is i such that x i < 0, then the LCP(M (ξ), q(x, ξ)) does not have a solution; if x = 0, then any v(ξ) with y(ξ) = 0 and where · is the Euclidean norm. The following lemma gives the form of the least norm solution of the LCP(M (ξ), q(x, ξ)) and shows that it is the unique solution of the LCP(M (ξ), q(x, ξ)) when x > 0. Lemma 3.1. For any x ≥ 0 and ξ ∈ Ξ, the LCP(M (ξ), q(x, ξ)) has the least norm solution v * (ξ) = (y * (ξ), s * (ξ)) T with where Π [0,x] (z) is the projection from z to the set [0, x]. Moreover, the least norm solution is the unique solution of the LCP(M (ξ), q(x, ξ)) if x > 0. Lemma 3.1 shows that for a fixed vector x ≥ 0, the least norm solution of the LCP(M (ξ), q(x, ξ)) is uniquely defined by H(ξ), ρ(ξ) and γ(ξ). In real applications, noise may be presented in the data set and we consider the perturbation bound of the solution regarding noise in H(ξ), ρ(ξ) and γ(ξ). Let us consider the following LCP(M ,q) with From Lemma 3.1, the LCP(M ,q) with x ≥ 0 has the least norm solutionv = (ū,t) T in the following form The following theorem provides the distance y * (ξ) −ū regarding the noise in data ρ(ξ), H(ξ), γ(ξ). Theorem 3.1. Suppose that x ≥ 0 and let Γ > 0 such that Γγ(ξ) ≥ 1. Then the least norm solution of the LCP(M (ξ), q(x, ξ)) is continuous with respect to H(ξ), ρ(ξ) and γ(ξ). Moreover, we have the following perturbation error bound In this section, we consider how to efficiently solve problem (8) with a finite number of realizations of the random variable. For Ξ ν = {ξ 1 , . . . , ξ ν }, the probability for each realization is 1/ν. For simplicity, we set B = (0, In such setting, problem (8) is the standard LCP(M,q) and the progressive hedging algorithm for the LCP(M,q) is as follows. Algorithm 1: Progressive Hedging Algorithm (PHA) [41] Step 0. Given an initial point k , and for = 1, . . . , ν, update Step 2. Set k := k + 1; go back to Step 1. The PHA has been widely used for solving stochastic optimization problems and stochastic variational inequalities. Due to the special structure of the M, we propose a new algorithm for solving the LCP(M,q), which is called Alternating Block Algorithm (ABA). At each iteration of the PHA, we need to solve ν linear complementarity problems in 3J dimension. At each iteration of the ABA, we only need to solve ν strongly convex quadratic programs with a simple constraint in J dimension and one linear complementarity problem in J dimension. Hence the computation cost of the ABA is much less than that of the PHA. Step 0. Given an initial point x 0 ∈ R J with x 0 ≥ 0. Set k = 0. Step 1. For = 1, . . . , ν, find y k that solves quadratic program Set Step 2. Find x k+1 that solves the LCP problem Set k := k + 1; go back to Step 1. Let C 2 ∈ R J×J be the symmetric positive definite matrix such that C T where M (ξ ) K is the nonsingular principal submatrix of M (ξ ) whose entries of M (ξ ) are indexed by the set K ⊆ {1, . . . , 2J}. Now we use randomly generated problems to compare the performance of the PHA and ABA. With uniform distribution, we randomly generate each element of vectors a, c,β on [0, 1],h on [2, 3] , and numbersγ on [0, 0.5] andᾱ on [5, 10] . Let r = 0.5e and C ii = 10 + c i + (r T e + (J − 2)r i ). We generate sample ξ 1 , ξ 2 , . . . , ξ ν from uniform distribution on [1, 2] and set Then, we let We terminate the two algorithms when one of the following three stop criteria is met: The number of iterations reaches 400, or v k − v k−1 ≤ 10 −6 , or Res := min(Mv + q, v) ≤ 10 −6 . The initial points for PHA and ABA are v 0 = ((−(C + re T ) −1 a) + , 0, . . . , 0) T and x 0 = −(C + re T ) −1 a) + , respectively. The step size t in the PHA is set to t = 1. We choose J = 5, 15, 20, and increase the sample size ν from 5 to 1000. generate 10 problems following the description above. The ABA and PHA methods are used to solve these 10 problems. The results reported in Table 1 are the average of iterations, cpu times(seconds), residuals and initial residuals. Table 1 shows that the ABA can solve the LCP(M, q) with all J and ν efficiently. Another advantage of the ABA is that the iteration numbers remain almost the same when J and ν increase, and the cpu time is roughly linearly increasing as ν increases. Moreover Table 1 shows that although the PHA can solve the problems with small J and ν, it fails to solve the problems with large J or ν within 400 iterations. In this section, we use the real data of oil price, demand and the market share of 14 major oil producers in the last 17 months to demonstrate the predicability of the twostage stochastic LCP model with adaptive parameters in the cost functions. Moreover, the model and the simulation results rationalize decisions made by major producers during COVID-19 pandemic. Since its first identification, COVID-19 has become a pandemic across all continents of the globe. To fight back the disease, large portions of social and commercial activities have been restricted if not suspended entirely. It seemed, during the worst of the outbreak, the whole world had come to a halt, so did the consumption of oil. This far reaching global event contributed to an already fragile market of oil. As a result of the pandemic, demand in transportation and factory output fell, which leads to decrease in overall demand for oil. Consequently, it caused oil prices to fall deeply. The last straw happened on 6 March 2020, when Russia rejected the demand of OPEC on further production cuts in response to the demand shrinkage. The oil price fell following the Russian announcement. We treat the oligopoly market of oil as a two-stage stochastic game where oil producers compete for profit by deciding the optimal production at the first stage of the game. In particular, we consider 15 oil producers in the following list as 15 agents in the game. Producer i makes a decision on its oil production quantity x i , based on its predicted further characteristics of the oil market at a later time, where the trading actually occurs. We suppose the trading occurs at the second stage where producer i supplies part if not all of its produced quantity in the first stage to generate revenue. The spot price of the trading is uncertain at the time of production decision, and is mainly depended on the demand supply relation at the second stage. To be more precise, supply quantity y i (ξ) of producer i is uncertain being further event at the time of production decision. To simulate the effect of excessive supply on price, we adapt a simple supply and demand relation and express the inverse demand function as: where p 0 (ξ) is the benchmark price and γ(ξ) > 0 represents the negative effect on price if there is excessive supply of oil with respect to the observed demand D(ξ) at the second stage. Traditionally, for the production of every barrel of crude oil, producers need to explore oil fields before building the extraction site, even after the oil is extracted refinement and shipment require both time and labour not to mention the cost that involves. The nature of the oil production has been evolving with the technological advance which enables, e.g., extraction from oil sands. There exist fundamental differences in energy infrastructures between traditional producers and oil sands producers, e.g., US, Canada [31] , and it is expected to be reflected at the first stage in choosing production strategies. In this paper, we assume that the costs in both stages are quadratic as expressed in (5) and (6) where the parameters c i , a i , h i (ξ), γ(ξ) are to be learned from market data. We also include an extra term (r i J j=1 x j )x i in the first stage to represent the strategic concern of producer i in response to global production quantity J j=1 x j . Note that we have no restriction on the sign of r i and when it is positive it represents the fact that producer i is willing to decrease its production when the global production is high. Typical of such agents are often price setters, since the action would give a boost to the oil price and may be more profitable despite the production cuts. In reality, Russia's refusal to production cuts agreement triggered huge volatility of the market during the COVID-19 pandemic. To express it with our model, it means that Russia adapted its value of r i to be negative given the market prediction back in March 2020. That is to say, Russia made a "squeeze" strategy and decided that if the price drop caused by over supplying is under control, the maintained market share is potentially more important for long term profitability. Other notable decisions of producers during the COVID-19 pandemic include that Russia refused to production cut on 6th March 2020; Saudi Arabia offered price discount on 8th March and increased production; U.S. demanded production cuts on 2nd April 2020; OPEC and Russia cut production on 9th April 2020. We will use our model to show that most decisions of producers were reasonable from the prospective of the producers, since they had different tolerance on low oil prices, and they all attempted the best action for their own benefits. The market data used in our study are obtained from the following sources. (i) Statistical Review of World Energy 1 , latest publish in June 2020 by bp. Inc; This data set provides average daily production quantities of major oil-producing countries, and the percentage of their production of the total production is regarded as their market shares respectively. Tables 3-4 . • a i : This parameter represents the linear production cost of producer i. It is widely agreed that traditional producers have very low unit cost of oil production. As for the oil sand producers, the unit cost is much higher. In our simulation, a i = c i , i = 1, 2, 4, 5, 7, . . . , 15, a 3 = 6c 3 , a 6 = 2c 6 . • r i : This parameter represents producer i's response to the total of production by all producers. Before the pandemic, it was widely agreed that supply should be kept in accordance to the demand but little preference was taken till Russia's refusal to further production cuts. In our simulation, r i = 0 for all producers in 2019, and r i were given in Table 2 for different producers in 2020. These are basic production cost parameters restricted by technological advance and complicated operations, and we do not expect them to change over short periods of time for all producers. For the purpose of forecasting current year production, these parameters are revised monthly taken based on the market share of the month before. The stochastic parameters are the risk-adjusted spot price ρ i (ξ) = α(ξ) − β i (ξ) and H(ξ) = diag(h i (ξ) + γ(ξ)), where α(ξ) is benchmark price, γ(ξ) is stochastic supply discount, α i and β i are the supply cost coefficients. For our experiments, we randomly choose ζ ∈ [0.05, 0.1], and let h i = β i = ζ × a i representing 5% to 10% of the unit production cost. The data (ii) gives the crude oil price change due to different factors of contributions, namely contribution of demand ∆D, supply ∆S and the residual contribution ∆R. Then, price change ∆price is computed as follows: ∆price = ∆D + ∆S + ∆R. These contributions ∆D, ∆S and ∆R over certain period of time are uncertain. We assume that it can be described by random variable ξ with unknown distribution, written as d(ξ), s(ξ) and r(ξ). For the purpose of our numerical tests, we formulate empirical distributions of historical data and use them as an approximation to the unknown distributions of different factors of contribution respectively. Recall that in our model the price is given by in which the demand D(ξ) is ignored as it would be a constant term and has no effect Table 4 : Average of daily market shares of producers in each month of January-May 2020 on solution. Then, for any realization of α k (ξ ) of k-th day, it corresponds to where α k 0 is the known price of prior day given in (iii), d k (ξ ) and r k (ξ ) are random scenarios taken from empirical distributions of d(ξ) and r(ξ), respectively. It follows that, we can generate a set of data of stochastic supply discount γ(ξ ) where absolute value | · | ensures that increase in quantity has a negative influence on price, ξ ∈ [0.99, 1.01] is uniformly distributed andη is the total supply obtained from data (i). We chose sample size ν = 800 of random variable ξ for both in-sample and out-of-sample in the numerical simulation. For long-term prediction (yearly market shares prediction), we refer interested readers to [25] for more details. Here, we focus on short-term prediction, namely the monthly in-sample and out-of-sample market shares. Table 3 gives average of daily market shares of producers in each month of 2019. Figures 1 and 2 display results for the recovered monthly market shares in 2019. For each month, the first column is the real market share, while the second and third column are the in sample and out sample recovered results, respectively. They show that our two-stage stochastic LCP model recovers and predicts the short-term real market shares from January 2019 to May 2020 very well. Although global oil demand has been hit hard by COVID-19 and oil price has fell to historically low, our results show that a Nash-equilibrium for the global oil market share during the COVID-19 pandemic can be expected. Of particular interests are what had happened in March and April among Russia, Saudi Arabia and U.S.A. In particular, the Brent spot price fell from around $70 per barrel to about $50 since the identification of the pandemic and was believed to decrease further. If no actions of change were taken by the major producers, what could have happened is that the high-cost producers would be forced to cut production because they are more vulnerable in low price environment. On one hand, if the production cut is significant, the price will boost and the non-limited producers can be more profitable. On the other hand, if the oil price continues to fall, the optimal production decision has little to change for the traditional producers with low unit cost, e.g., Russian, Saudi Arabia. In indeed, with our model of two-stage stochastic game, the estimated market share given no strategic changes can be seen in accordance with the description above. It is, given the estimated market share, the most rational decision for Russia is to refuse the production cut agreement. Saudi Arabia, who believed to have the lowest unit production cost, followed the strategy of Russia immediately by offering price discount and increased its production. For both countries, the low price environment has little effects in the sense of maintaining their market share respectively. The same cannot be said for U.S.A., who has high cost in unit oil production, and is estimated to loss its market share if not to change its production strategy, presented by choosing non-zero r in our model. It is most apparent from Table 2 , that all the major producers responded in choosing their strategies to increase production quantities. The reality was counter intuitive at first glance, but those were in fact all rational decisions and can be forecasted by our model In this paper, we model the oil market share using the two-stage stochastic LCP (8) via the two-stage stochastic games (1)-(2) with quadratic concave utility functions (5)- (7) . We show the existence and uniqueness of a Nash-equilibrium for the oil market share by the solution of the two-stage stochastic LCP (8) . Moreover, we propose the Alternating Block Algorithm (ABA) to find a solution of the two-stage stochastic LCP (8) with finite realizations. We derive the convergence theorems of the ABA and the PHA for solving (8) with finite realizations and show the ABA is much faster than the PHA by randomly generated problems. We apply the new theoretical results and the ABA to the analysis of the oil market share from January 2019 to May 2020. Simulation results show that our model is effective for forecasting the oil market share under an uncertain environment of COVID-19 pandemic. Proof. It is easy to see that v(ξ) = (x, |ρ(ξ)|) T ≥ 0 and M (ξ)v(ξ) + q(x, ξ) ≥ 0, that is, (x, |ρ(ξ)|) T is a feasible solution of the LCP(M (ξ), q(x, ξ)). Hence from the feasibility and the positive semi-definiteness of M (ξ), the LCP(M (ξ), q(x, ξ)) has at least one solution [14, Theorem 3.1.2] . The LCP(M (ξ), q(x, ξ)) is the first order optimality condition of the strongly convex quadratic program min 1 2 y(ξ) T (H(ξ) + γ(ξ)ee T )y(ξ) + ρ(ξ) T y(ξ) s.t. 0 ≤ y(ξ) ≤ x. (A.1) Because of the strong convexity, the first order optimality condition is necessary and sufficient for the unique optimal solution y * (ξ) of problem (A.1), which is a fixed point of the fixed point problem y(ξ) = Π [0,x] (y(ξ) − ρ(ξ) − (H(ξ) + γ(ξ)ee T )y(ξ)). Hence, y * (ξ) has the form of (13) . By the definition of the projection, we can easily obtain (ρ(ξ) + (H(ξ) + γ(ξ)ee T )y * (ξ)) i Hence, the multiplier s * (ξ) in (14) with y * (ξ) in (13) is a solution of the LCP(M (ξ), q(x, ξ)), where for the case y * i (ξ) < x i it is from s * i (ξ) = (ρ(ξ) + (H(ξ) + γ(ξ)ee T )y * (ξ)) i = 0 and s * i (ξ)(x i − y * i (ξ)) = 0, and for the case y * i (ξ) = x i , it is from y * i (ξ)(s * i (ξ) + (ρ(ξ) + (H(ξ) + γ(ξ)ee T )y * (ξ)) i ) = 0. If y * i (ξ) = x i > 0, s * i (ξ) in (14) is uniquely defined. If y * i (ξ) = x i = 0, then s * i (ξ) = max(0, −(ρ(ξ) + (H(ξ) + γ(ξ)ee T )y * (ξ)) i ) = 0. Hence v * (ξ) = (y * (ξ), s * (ξ)) T defined in (13)- (14) is the least norm solution of the LCP(M (ξ), q(x, ξ)) and the unique solution of LCP(M (ξ), q(x, ξ)) if x > 0. We complete the proof. Proof. We first prove (15) . Let w(ξ) = (H(ξ) + γ(ξ)ee T )y * (ξ) + ρ(ξ) andw = (H +γee T )ū +ρ. Then we have y * (ξ) = Π [0,x] (y * (ξ) − w(ξ)) andū = Π [0,x] (ū −w), which implies mid(y * (ξ), y * (ξ) − x, w(ξ)) = 0 and mid(ū,ū − x,w) = 0, where "mid" is the componentwise median operator. Following the proof of Lemma 2.1 in [8] , there is a diagonal matrix D =diag(d 1 , . . . , d J ) with 0 ≤ d i ≤ 1 such that 0 = (I − D)(y * (ξ) −ū) + D(w(ξ) −w). Hence we obtain (I − D + D(H(ξ) + γ(ξ)ee T ))(y * (ξ) −ū) = −D(ρ(ξ) −ρ + (H(ξ) −H + (γ(ξ) −γ)ee T )ū). (A.2) Since H(ξ) + γ(ξ)ee T is symmetric positive definite, by Theorem 2.7 in [12] , we have max d∈[0,1] J (I − D + D(H(ξ) + γ(ξ)ee T )) −1 D = (H(ξ) + γ(ξ)ee T )) −1 ≤ γ(ξ) −1 (I + ee T ) −1 . Therefore, using (I + ee T ) −1 = 1, ee T = J, γ(ξ) −1 ≤ Γ and 0 ≤ū ≤ x, we obtain (15) from (A.2). Next we show the continuity of the last J-components of the least norm solution of the LCP(M (ξ), q(x, ξ)). Without loss of generality, assume that σ = {i ∈ {1, . . . , J} | y * i (ξ) < x i } = ∅. From (15), for any 0 < < min i∈σ x i − y * i (ξ), there is δ > 0 such that if ρ(ξ) −ρ + H(ξ) −H + γ(ξ) −γ < δ, then y * (ξ) −ū < . This implies for i ∈ σ, Hence, we havet i = s * i (ξ) = 0. For i ∈ σ, that, is, y * i (ξ) = x i , we have s * i (ξ) = −(ρ(ξ) + (H(ξ) + γ(ξ)ee T )y * (ξ)) i ≥ 0. Ift i = 0, then −(ρ + (H +γee T )ū) i ≤ 0, otherwiset i = −(ρ + (H +γee T )ū) i ≥ 0. Hence, using (15) where L ≥ max{1, J x }. Hence, the solution of the LCP(M (ξ), q(x, ξ)) is continuous with respect to H(ξ), ρ(ξ) and γ(ξ). We complete the proof. Forecasting the term structure of crude oil futures prices with neural networks OPEC vs US shale oil: Analyzing the shift to a market-share strategy Evolutionary Neural Network model for West Texas Intermediate crude oil price prediction Expected residual minimization method for stochastic linear complementarity problems Two-stage stochastic variational inequalities: an ERM-solution procedure Convergence analysis of sample average approximation of two-stage stochastic generalized equations Discrete approximation of two-stage stochastic and distributionally robust linear complementarity problems Computational error bounds for differential linear variational inequality Stochastic Equilibrium and Variational Inequalities Stochastic variational inequalities: Residual minimization smoothing sample average approximations Computation of error bounds for P-matrix linearcomplementarity problem Perturbation bounds of P-matrix linear complementarity problems Newton iterations in implicit time-stepping scheme for differential linear complementarity systems The Linear Complementarity Problem Generation capacity expansion in a risky environment: a stochastic equilibrium analysis Engineering and economic applications of complementarity problems Sample-path solution of stochastic variational inequalities Generation capacity investments in electricity markets: perfect competition Approximations of Nash equilibria Dynamic oligopolistic games under uncertainty: a stochastic programming approach A Benders decomposition method for solving stochastic complementarity problems with an application in energy Causes and Consequences of the Oil Shock of 2007-08 How does oil price volatility affect non-energy commodity markets? Oil price volatility and oil-related events: An Internet concern study perspective Regularized two-stage stochastic variational inequalities for Cournot-Nash equilibrium under uncertainty Speculation in the oil market Worldwide energy use across global supply chains: Decoupled from economic growth? Strategic behavior in power markets under uncertainty Addressing supply-side risk in uncertain power markets: stochastic Nash models, scalable algorithms and error analysis Where is an oil shock? Energy infrastructure modeling for the oil sands industry: Current situation Gold and oil futures markets: Are markets efficient? Who's afraid of a big bad oil shock Preventive crude oil scheduling under demand uncertainty using structure adapted genetic algorithm Two-stage non-cooperative games with risk-averse players Convergence of stationary points of sample average two stage stochastic programs: a generalized equation approach On the characterization of solution sets of smooth and nonsmooth convex stochastic Nash games Scenarios and policy aggregation in optimization under uncertainty Variational Analysis Stochastic variational inequalities: single-stage to multistage Solving monotone stochastic variational inequalities and complementarity problems by progressive hedging Oil price volatility and economic growth: Evidence from advanced economies using more than a centurys data Modeling and computing two-settlement oligopolistic equilibrium in a congested electricity network How does investor attention affect international crude oil prices? Xiaojun Chen would like to thank Hong Kong Research Grant Council for grant PolyU153001/18P.