Optimal Variable Renewable Energy Generation Schedules Considering Market Prices and System Operational Constraints

: The maximization of output from variable renewable energy (VRE) sources considering system operational constraints (SOCs) is a traditional method for maximizing VRE generators’ proﬁts. However, in wholesale electricity markets, VRE participation tends to reduce marginal prices (MP) because of its low marginal costs. This circumstance, called the “merit-order effect” (MOE), reduces the generators’ proﬁts. Thus, the traditional method is possibly no longer the best and only method to maximize the generators’ proﬁts. Moreover, the VRE support schemes also affect MP, making MOE more severe. VRE curtailment can relieve MOE, but VRE output must be decreased, thereby reducing the generators’ proﬁts. This paper proposes a method to ﬁnd the optimal VRE generation schedules that maximize VRE generators’ proﬁts while considering the trade-off among the VRE output, MP, and SOCs. The method combines the merit-order model and the unit-commitment model solved by the optimization tools in MATLAB. Thailand’s electrical system was the test system. The result shows that VRE generators’ proﬁts from the proposed method are signiﬁcantly higher than from the traditional method when the system has high wind penetration, and the generators have no support scheme. Curtailing approximately 7–10% of wind output can increase the average MP by 23.6–30%.


Introduction
Renewable energies (RE) continue to grow in importance for electrical systems because of a rapid decline in investment costs, their use of free energy from nature and environmental friendliness [1]. RE generation resources include both dispatchable RE sources, such as hydro, and non-dispatchable RE sources known as variable renewable energies (VRE), such as solar and wind. VRE forms a substantial proportion of RE targets in many countries, including Thailand [2,3]. In wholesale electricity markets, VRE is generally prioritized for supplying electricity because of its low marginal costs (MC). A traditional method to leverage low-cost resources, i.e., VRE, is to maximize its generated electricity (output). The maximization basically considers system operational constraints (SOCs) consisting of electrical system constraints and generation characteristic constraints. However, when there is VRE proportion in the markets, the marginal price (MP) at that time is inevitably diminished. The greater the VRE output, the greater the drop in MP. Moreover, VRE generators often supply electricity uncorrelated with electricity demand because of their non-dispatchable characteristics. The output could highly exceed the demand during windy or sunny hours, contributing to very low MP (possibly zero or even negative). Therefore, VRE output affects MP during the day in both quantity and distribution form. This circumstance is called the "merit-order effect" (MOE) [4][5][6][7][8][9][10][11][12][13][14][15][16].
Any generators in wholesale electricity markets gain their revenue based on MP; thus, MOE contributes to reductions in the generators' revenue. Additionally, many countries around the world provide support schemes for RE generation to reach their RE targets. The common support schemes are feed-in tariffs (FIT), contract for difference feed-in tariffs

Merit-Order Effect
In wholesale electricity markets with the current market design, the MP at a specific time is the MC of the last power plant needed to meet the electricity demands at that time. Among all generators offering their energy to the markets, VRE generators are prioritized because of their low MC. When VRE generators supply their output into the markets, the most expensive generators on the markets are driven out, and the MP is diminished. The greater the VRE output, the greater the drop in MP. Moreover, VRE generators often supply electricity with no regard for demand because of their non-dispatchable characteristics. The output could highly exceed the demand during windy or sunny hours, contributing to very low MP (possibly zero or even negative). That means VRE output affects the MP quantity and distribution during the day. This circumstance is called the "merit-order effect" (MOE). Any generators in wholesale electricity markets gain their revenue based on MP; thus, MOE contributes to reductions in generators' revenue [4][5][6][7][8][9][10][11][12][13][14][15][16].
Many studies confirmed that MP is declined by VRE penetration and its output. Research [27] found that the MP was decreased around 0.63 $/MWh in Germany and 0.95 $/MWh in Spain per additional percent of wind infeed. In the case of Italy, [28] proved that 1 GWh from solar and wind reduces the average MP by 2.73 $/MWh and 4.99 $/MWh, respectively. In the United States, research [11] found that the average loadweighed MP for each additional percentage of VRE penetration declined by 0.2-0.9 $/MWh (CAISO, NYISO, SPP, and ERCOT), similarly to [12], in which case it was 0.1-0.8 $/MWh (CAISO and ISO-NE). Research [29] indicated that European wholesale electricity prices had dropped by nearly two thirds since their all-time high in 2008. The largest factor depressing the prices was the expansion of VRE. VRE is expected to become competitive in energy markets. However, the competitiveness might not be enough to ensure profitability on wholesale electricity markets if the MP falls too low [30][31][32][33].
The methods to illustrate MOE are based on two main approaches: first, the development of electricity market models which simulate the operation of a wholesale electricity market and calculate the resulted MP for various scenarios; second, the regression analysis approach, which uses historical prices and generation data to quantify the actual reductions in MP for a given period. Both the approaches were combined in some studies [15].

Renewable Energy Support Schemes
Many countries provide support schemes for RE generation to reach their RE targets. The schemes are not only provided to VRE generators, but also other RE generators, such as biomass and hydro plants, but this paper focuses on the support provided to VRE generators. The schemes help covering the cost disadvantages faced on liberalized electricity markets [34]. However, the schemes involve prices VRE generators have offered to the markets [18] that unintentionally affect MP and probably make MOE more severe. The common support schemes are feed-in tariff (FIT), contract for difference feed-in tariff (CFD-FIT), and feed-in premiums (FIP). These support schemes have been applied in 23 out of 27 EU countries [17], and many other countries around the world. The approaches of the schemes are as follows: • Feed-in Tariff (FIT): Generators receive a fixed price per kWh for each unit of electricity generated, differing according to the generation sources (wind, solar, etc.) [34]. The fixed prices, which are independent from the MP, are mostly determined by the government. This means that generators do not receive any revenue directly from the markets [35]. • Contract for difference feed-in tariff (CFD-FIT): Generators receive a fixed price per kWh for each unit of electricity generated. The price called the "strike" price or "reference" price is established by the government through bidding. At a specific time, generators sell their energy at the MP that can be above, below, or the same as the strike price. If the MP is equal to the strike price, then there is no further action. If the MP is below the strike price, generators will get payment on top of the MP to reach Feed-in premiums (FIP): Generators receive the MP from the market and an additional fixed payment per kWh on top of the MP. The fixed payment could vary according to the associated risk sharing between the generators and the public [35].
The differences in support schemes are shown in Figure 1. FIT and CFD-FIT-supported generators receive fixed revenue regardless of the MP. FIP-supported generators' revenue depends on the MP at a specific time. In wholesale electricity markets, generators need to ensure that they will be committed to selling energy; thus, they will offer the lowest prices they can accept without loss. Supported VRE generators will offer negative prices equal the support prices they receive, and VRE generators without support will offer their MC [18].

•
Contract for difference feed-in tariff (CFD-FIT): Generators receive a fixed price per kWh for each unit of electricity generated. The price called the "strike" price or "reference" price is established by the government through bidding. At a specific time, generators sell their energy at the MP that can be above, below, or the same as the strike price. If the MP is equal to the strike price, then there is no further action. If the MP is below the strike price, generators will get payment on top of the MP to reach the strike price. If the MP is above the strike price, generators have to pay back the difference [35,36].
• Feed-in premiums (FIP): Generators receive the MP from the market and an additional fixed payment per kWh on top of the MP. The fixed payment could vary according to the associated risk sharing between the generators and the public [35].
The differences in support schemes are shown in Figure 1. FIT and CFD-FIT-supported generators receive fixed revenue regardless of the MP. FIP-supported generators' revenue depends on the MP at a specific time. In wholesale electricity markets, generators need to ensure that they will be committed to selling energy; thus, they will offer the lowest prices they can accept without loss. Supported VRE generators will offer negative prices equal the support prices they receive, and VRE generators without support will offer their MC [18]. There are other RE support schemes, such as green certificates, investment support, low interest loans, and tax exemptions [35]. The details of them depend on the policies of the countries implementing them. However, they are not directly relevant to energy selling. Such schemes are beyond the interest of this paper.

Variable Renewable Energy Curtailment
VRE curtailment is a reduction in the output of wind or solar generation from the output possible with the available wind or sunlight [18][19][20][21][22][23]. There is no standard method to measure curtailment. However, the common metric to measure it is as a percent of the output that the generation could have produced [21]. Many studies stated that VRE curtailment levels grow with VRE penetration [24].
There are both technical and economic reasons to curtail VRE. For technical reasons, the most common ones are to avoid insufficient transmission, local congestion, and excessive supply during low-load periods (oversupply). It is possible that these different reasons correlate in time. The curtailment called technical curtailment is done by system operators; the curtailed generators could gain compensation or remuneration for their cur- There are other RE support schemes, such as green certificates, investment support, low interest loans, and tax exemptions [35]. The details of them depend on the policies of the countries implementing them. However, they are not directly relevant to energy selling. Such schemes are beyond the interest of this paper.

Variable Renewable Energy Curtailment
VRE curtailment is a reduction in the output of wind or solar generation from the output possible with the available wind or sunlight [18][19][20][21][22][23]. There is no standard method to measure curtailment. However, the common metric to measure it is as a percent of the output that the generation could have produced [21]. Many studies stated that VRE curtailment levels grow with VRE penetration [24].
There are both technical and economic reasons to curtail VRE. For technical reasons, the most common ones are to avoid insufficient transmission, local congestion, and excessive supply during low-load periods (oversupply). It is possible that these different reasons correlate in time. The curtailment called technical curtailment is done by system operators; the curtailed generators could gain compensation or remuneration for their curtailed energy based on the regulations of each system. For economic reasons, VRE curtailment contributes to significant savings in both grid and storage extension investments. Avoiding the curtailment would require investing in transmission lines and storage, which would be very costly if it were only used for a few hours per year. The curtailment called economic curtailment is done by both system operators to minimize system costs and generators to maximize their profits [18,[21][22][23][24]. Moreover, MP is decreased by VRE output, as mentioned in the previous section. In markets, a small amount of capacity at the steepest part of the merit-order curve makes a significant difference in MP. One of generators' strategies in bidding is withholding electricity generation with low marginal costs from the market to increase the MP when the MP is low. VRE generators might curtail their output following their strategic bidding to maximize profits from the markets [18,25,26]. On the other hand, curtailment can be problematic, since it decreases the capacity factor of the generators. When the electricity generation is lower than intended by the design, it can be considered as a reduction in generators' revenue from selling less electricity than their capability [4, [7][8][9][10]. Curtailment also decreases generators' ability to recover their capital costs because of the reductions in the revenue [21]. Compensation to generators for revenue loss from curtailment varies greatly across the U.S. and Europe [23]. For technical curtailment, the costs terms of lost generation are discussed based on MP and support levels including the rationale for compensating the curtailed energy. For economic curtailment, it is allowed without compensation [18].

Proposed Method
This paper proposes a method to find the VRE generation schedules that maximize the profits of VRE generators while considering the trade-off among the amount of VRE output, the MP, and the SOCs. Moreover, the VRE support schemes involving the prices of VRE offered to the market were considered. The method in this paper is the combination of the merit-order model (Section 5.1) and the unit-commitment model (Section 5.2). The first model is for optimizing VRE output, and the second one is for satisfying the SOCs. The traditional method to maximize VRE generators' profits, which is the maximization of VRE output, was also demonstrated to compare the VRE profit with the one from the proposed method (Section 5.3).

The Merit-Order Model
The merit-order model simulates the operation of wholesale electricity markets for a variety of cases. In wholesale electricity markets with the theoretically perfect competition, generators offer two parameters to the markets at a specific time: first, their capability to produce energy; second, the price they would like to sell their energy at. Generators will offer the lowest price they can accept without loss to make sure that they can be committed to selling the energy [18]. The merit-order curve over time is determined by those offers. The MP at that time is then set by the intersection of the merit-order curve and the electricity demand at the time. There is no market participant that is able to affect the MP [37]. However, in real-life situations, generators could use strategies to drive up MP to gain more profit from the markets. First, they could curtail their output (offered less energy). Second, they could offer to sell their energy at high prices. These two different approaches lead to the same results: higher MP, higher profits, and withheld output [25,26]. This paper applied the first approach into the merit-order model to illustrate the relationship between VRE output and MP. Additionally, the VRE support schemes involving the prices VRE offered to the market were included in the model. Figure 2 demonstrates the concept of the merit-order model. In Figure 2, the support schemes are classified into two types: FIP schemes where VRE generators receive the fix support price on top of the MP; and FIT and CFD-FIT schemes where VRE generators receive only the fixed price. Both types of supported VRE generators will offer the negative prices equal to the support prices they receive. VRE In Figure 2, the support schemes are classified into two types: FIP schemes where VRE generators receive the fix support price on top of the MP; and FIT and CFD-FIT schemes where VRE generators receive only the fixed price. Both types of supported VRE generators will offer the negative prices equal to the support prices they receive. VRE generators without support will offer their MC. This strategy guarantees that they can be committed to selling the energy. Moreover, even if they are the last power plant committed to supplying energy (marginal unit), the MP will be at least equal to their support prices (if they are supported) or MC (if they are not supported). Therefore, their revenue gained from both the MP and support schemes is at least zero. That means no negative revenue from selling energy is possible.
In a specific time, the MP depends on the electricity demands at that time, the energy offered by the VRE (E VRE ), the energy offered by other generators, and the prices offered by all other generators. Note that this paper focuses on the relationship of VRE output and MP. Thus, other parameters involved in MP, such as all thermal and hydro generators' energy offers and price, were assumed to be fixed. Their collective offered energy was assumed to be the maximum energy they could provide, and their offered prices were assumed to be their MC at their maximum capability. Lastly, consumers were assumed to not react to the MP. In Figure 2, if VRE-offered energy is increased from E VRE (A) to E VRE (B), the merit-order curve will be shifted to the right. The MP will decline from MP(A) to MP(B). E VRE always affects the MP regardless of the support schemes the VRE generators receive because they shift the merit-order curve. Thus, a greater E VRE , contributes to a greater drop in MP (MOE), whereas a low E VRE means generators sell less electricity. As a result, if generators offer the optimal E VRE into the markets, they will gain the maximum profits.
In the merit-order model, the objective function is the maximization of daily VRE profit. The optimal VRE generation schedules of a considered day are determined at a resolution of one hour. The total profits of all VRE generators are maximized, rather than the profits of each individual generator, to avoid sub-optimal results. All VRE generators in the system are classified into two groups based on their resources, i.e., solar and wind. Thus, any parameters relevant to VRE in this paper refer to the total values of all solar or wind generators in the system.
The VRE generators' daily profits are calculated by summing the generators' hourly revenues (Revenue(t)), and subtracting the generators' hourly variable costs (VC(t)), and their capital costs per day (CC), as shown in Equation (1); t is a specific time.
As shown in Figure 1, VRE generators that receive no support scheme will gain their revenue only from the MP. FIP-supported VRE generators will gain their revenue from the MP and the FIP support price (SP FIP ). FIT and CFD-FIT supported VRE generators will gain their revenue only from the FIT support price (SP FIT ). Both SP FIP and SP FIT are constants. The Revenue(t) is determined by summation of the MP and the support schemes (if any) multiplied by the VRE output (E(t)). The Revenue(t) calculations differentiated by VRE support schemes and resources are shown in Equations (2) and (3). The MP(t) is the function of the merit-order curve and electricity demand (D(t)) at a specific time, as shown in Equation (4).
; Solar generators without support (MP(t) +SP FIP,Solar ) × E Solar (t) ; Solar generators with FIP support SP FIT,Solar × E Solar (t) ; Solar generators with FIT support ; Wind generators with FIT support The VC(t) is the generators' hourly variable costs calculated by multiplying their marginal costs (MC) and their E(t), as shown in Equations (5) and (6). CC is the generators' capital costs per day calculated, as shown in Equations (7) and (8). C is the generators' capital costs per installed capacity. ICAP is the generators' total installed capacity. Note that MC, C, and ICAP are constants.
The objective function is optimized, subject to the VRE resource constraints. The E(t) have to be less than or equal to the VRE generation capability at the time (Profile(t)), which is determined by the available solar irradiance and wind speed, multiplied by the installed capacity (ICAP), as shown in Equations (9) and (10).
The outputs from the merit-order model are E Solar (t) and E Wind (t), which are the VRE outputs offered to the markets that provide the maximum profit to VRE generators. The time series of the VRE output during the day is herein referred to as the VRE strategic schedule. The merit-order model consisting of nonlinear multivariable functions was solved by the optimization tool "Fmincon" in MATLAB. Fmincon has an interior-point algorithm that can handle various types of nonlinear problems. Moreover, the algorithm uses little memory and can solve large problems quickly [38].

The Unit-Commitment Model
After getting the VRE strategic schedules from the merit-order model, the unitcommitment model was then used to find whether SOCs can be satisfied when VRE supplies energy follows the VRE strategic schedules. If the VRE strategic schedules contribute to the unsatisfiable SOCs in some period during the day, VRE strategic schedules will be modified by curtailing the output at the time.
The objective function of the unit-commitment model is based on the unit commitment problem (UCP) with a resolution of one hour. The UPC minimizes the daily VC of all thermal hydropower plants incurred from supply energy to support residual demand (RD), as shown in Equations (11) and (12); n is a given power plant, nthermal is the total number of thermal power plants in the system, and nhydro is the total number of hydropower plants in the system.
The VC(t) is calculated from the generators' MC multiplied by their E(t). Equation (13)  indicate the cost of producing one more MW of power from the plant. This paper treats the curve as a piecewise linear function.MC 1 n,thermal , MC 2 n,thermal , MC 3 n,thermal are the piecewise costs from the incremental cost curve of the thermal power plant n. P n,thermal (t) is the output power of thermal power plant n at time t; P min n,thermal is minimum output power of the thermal power plant n; and P R1 n,thermal , P R2 n,thermal , P R3 n,thermal are ranges of piecewise power derived from the incremental cost curve of the thermal power plant n. VC n,hydro (t) = MC hydro × E n,hydro (t) (13) VC n,thermal (t) =      MC 1 n,thermal × E n,thermal (t) ; P min n,thermal ≤ P n,thermal (t) ≤ P R1 n,thermal MC 2 n,thermal × E n,thermal (t) ; P R1 n,thermal < P n,thermal (t) ≤ P R2 n,thermal MC 3 n,thermal × E n,thermal (t) ; P R2 n,thermal < P n,thermal (t) ≤ P R3 n,thermal The objective function is optimized, subject to the SOCs: firstly, electrical system constraints such as serving electricity demand, committing must-run units, and providing operating reserves requirement that cover demand and VRE forecast errors, along with spinning reserves requirement for contingency events fixed by the N-1 approach; second, generation characteristic constraints, i.e., minimum/maximum generation, ramp capability, minimum up/downtime, and the limitations of hydro units which depend on the amount of water reserved on the considered day. The UCP is defined to be mixed-integer programming (MIP) because it can address issues with non-convexity related to the SOCs [39]. The UCP, which is mixed-integer linear function was solved by the mixed-integer linear programming optimization tool "Intlinprog" in MATLAB.
The outputs from the combination of the merit-order model and the unit-commitment model are optimal VRE generation schedules that provides the maximum VRE profits while considering the trade-off among the amount of VRE output, the MP, and the SOCs. The MP and VRE generators' profits are then determined according to the schedules. The outputs of the proposed method are the optimal VRE generation schedules, the MP, and the VRE generators' profits that are consistent.

The Traditional Method
This paper illustrates the VRE generators' profits when using the proposed method and when using the traditional method, which is the maximization of VRE output. The UCP was also used to find VRE generators' profits and generation schedules based on the traditional method. The objective function is shown in the Equation (15). The generators' (i.e., VRE, thermal, and hydro; D(t)) total VC in the system was minimized; thus, the output from low MC generators, i.e., VRE, was automatically maximized.

Min
The constraints of the optimization are also the electrical system and the generators' characteristic constraints, such as the constraints of the unit-commitment model mentioned in Section 5.2. The UCP of traditional method was also solved by the optimization tool "Intlinprog" in MATLAB. The outputs are the VRE generation schedules, the VRE generators' profits, and the MP based on the traditional method.
Later, the outputs from the proposed method and the traditional method are compared. Moreover, we performed sensitivity analysis on VRE penetration and proportion, and the characteristics of electricity demand (workday and holiday), to find the effects of these parameters on the optimal VRE generation schedules and the maximized profits. Figure 3 shows the flow chart of the optimization and sensitivity analysis.
tool "Intlinprog" in MATLAB. The outputs are the VRE generation schedules, the VRE generators' profits, and the MP based on the traditional method.
Later, the outputs from the proposed method and the traditional method are compared. Moreover, we performed sensitivity analysis on VRE penetration and proportion, and the characteristics of electricity demand (workday and holiday), to find the effects of these parameters on the optimal VRE generation schedules and the maximized profits. Figure 3 shows the flow chart of the optimization and sensitivity analysis.

Data
This paper used Thailand's electrical system as the sample system. Although the system has a vertically integrated structure, this paper assumed it as a liberalized structure. The load (electricity demand) profiles of both workdays and holidays, and VRE generation profiles, are shown in Figure 4. The forecast error of solar was around 12-16.6%, and that of wind was 6.7-12.4% (mean absolute percentage error, MAPE). The system has 18 hydropower plants, one thermal power plant consisting of 18 combined cycle gas turbines (CCGT) power plants, and 8 coal power plants. The generation characteristics of each power plant depend on the individual configuration. Table 1 shows a summary of the characteristics. The generation capital costs and operation costs, i.e., variable costs, startup costs, and load-following costs, are shown in Table 2. The VRE support schemes prices of many countries are collected in Figure 5.

Data
This paper used Thailand's electrical system as the sample system. Although the system has a vertically integrated structure, this paper assumed it as a liberalized structure. The load (electricity demand) profiles of both workdays and holidays, and VRE generation profiles, are shown in Figure 4. The forecast error of solar was around 12-16.6%, and that of wind was 6.7-12.4% (mean absolute percentage error, MAPE). The system has 18 hydropower plants, one thermal power plant consisting of 18 combined cycle gas turbines (CCGT) power plants, and 8 coal power plants. The generation characteristics of each power plant depend on the individual configuration. Table 1 shows a summary of the characteristics. The generation capital costs and operation costs, i.e., variable costs, startup costs, and load-following costs, are shown in Table 2. The VRE support schemes prices of many countries are collected in Figure 5.      [41] (Thailand is a Non-OECD country; if the analysis is used for OECD country, the discount rate would be 7% [42]). 2 Variable costs and startup costs (assumed to be all hot start) of CCGT and coal were provided by [40]. The exchange rate THB/USD was 33/1 (on 4 February 2021). The other data were from [41]. 3 Load following costs were taken from [43]. As shown in Figure 5, some countries provide different prices depending on the installed capacity of individual generators, and some proportions are substantially higher or lower than the others; thus, we calculated the medians of the data and used them as the VRE support schemes' prices for the calculation. The medians of the data: For sensitivity analysis, seven cases of VRE penetration and proportion are shown in Table 3. VRE penetration was varied and the VRE proportion in the system was divided into three types, i.e., wind based, solar based, and mixed, to investigate effects of system As shown in Figure 5, some countries provide different prices depending on the installed capacity of individual generators, and some proportions are substantially higher or lower than the others; thus, we calculated the medians of the data and used them as For sensitivity analysis, seven cases of VRE penetration and proportion are shown in Table 3. VRE penetration was varied and the VRE proportion in the system was divided into three types, i.e., wind based, solar based, and mixed, to investigate effects of system configuration on the results. More VRE penetration than 15 GW is impossible because the SOCs cannot be satisfied by the existing system's configuration. All the cases were sensitivity analyzed by adjusting types of VRE support schemes and the characteristics of electricity demand, i.e., workdays and holidays.

Result and Discussion
The method was applied to the test system to find the optimal VRE generation schedules that maximized VRE generators' profits while considering MP and system reliability. The traditional method is also illustrated to compare the results. The seven cases of the VRE penetration and proportion were sensitivity analyzed in more than 80 simulations, adjusting types of VRE support schemes and the characteristics of electricity demand, i.e., workdays and holidays. The VRE generators' daily profits and the VRE outputs from both methods during a workday and a holiday are shown in Tables 4 and 5, respectively. In the tables, the comparisons of the profit and output are presented in the "difference" columns.   From Tables 4 and 5, the overview of the results is that generators' profits depend on the quantity of output and the revenue they gain from the MP and support schemes. FIP-supported generators made the most profits in most scenarios because they received revenue from both MP and support schemes. However, both FIP and FIT-supported generators made considerably greater profits than the generators without support. Wind generators' profits were generally higher than solar generators' profits with the same installed capacity and support schemes, although the capital costs and the variable costs of wind generation are higher than solar. That is because the wind generation profile is more distributable than solar, contributing to a higher capacity factor and more profits.
To illustrate the benefit of the proposed method, the "difference" columns in Tables 4 and 5 show the differences between profit and output from the proposed method and the traditional method. If there is no difference (0), it means the results from the proposed method and the traditional method were the same; maximizing the VRE output is still the method that provides the maximized profit in that situation. The results show that the VRE output of the proposed method was diverse, depending on the support schemes that involve MP and revenue. The VRE output of traditional method was the same regardless of the support schemes because the method maximized VRE output in any case.
Moreover, the results prove that VRE generators' profits from the proposed method were higher than from the traditional method in the cases of the system having high VRE penetration and low electricity demands-i.e., cases 6 and 7 for workdays, and cases 4-7 for the holidays. For example, in case 6 for workdays, the profits from the proposed method were higher than for the traditional method by $1.267 million, though the output from the proposed method was less than that from the traditional method by 15.503 GWh.
That means selling less electricity to gain high MP (the proposed method) provided more profits than continually selling maximized electricity at low MP (the traditional method). The reason is that in cases where the system has high VRE penetration, yet the electricity demand is low, the MP is greatly diminished, contributing to severe MOE. The marginal units in these situations tend to be hydro or coal power plants that have low MC; thus, it is worth curtailing some VRE output to change the marginal unit to be CCGT power plants because the MP will be significantly driven up. Therefore, the proposed method provides more VRE profits than the traditional method on holidays than workdays, because the electricity demand on holidays is lower.
Additionally, the proposed method provided more profits than the traditional method in the scenarios where VRE generators received no support and VRE generators received FIP support schemes, because their generators' profits depend on the MP. However, the proposed method provided significantly more profits than the traditional method when VRE generators received no support because MP was the only factor involved in their profits. FIT-supported generators are always maximizing their outputs (same as traditional method) because they gain a fixed price for every MWh they produce, and their profits are independent of the MP.
Furthermore, selling electricity at low MP decreases the generators' ability to recover their capital costs and leads to negative profits. The proposed method avoided negative profits in cases 5-7 of the holidays. For instance, in case 5, the profits of solar generators from traditional method were −$0.310 million because the daily revenue did not cover the daily capital costs. The profits from the proposed method were $0.134 million. That means curtailing 0.165 GWh avoided the negative profits and provided more than $0.444 million in profits.
The results also showed that the differences in wind generators' profits and outputs between the proposed method and the traditional method were higher than for solar. This is because, from Figure 4, wind generation generates a high output when the demand is low, i.e., during 12 a.m. to 8 a.m.; thus, curtailing wind at the time to drive up the MP is worth it. Solar power generates electricity during daytime, i.e., 8 a.m. to 6 p.m. The electricity demand during that time is high. Thus, the MP at the time is high. As a result, selling maximal electricity provides maximized profits for solar generators, especially on working days. However, during 6 p.m. to 11 p.m., there is no curtailment of either wind or solar because the demand is high, but VRE output is low, contributing to high MP. Thus, maximized electricity production also provides maximized profits for VRE generators during that time.
The VRE optimal generation schedules and the MP of every case are presented in Tables 4 and 5. This paper presents a case where the proposed method provided greater profits than the traditional method, which is case 6, wherein VRE received no support. Figure 6 shows the optimal VRE generation schedules and the MP of that case. Figure 6a,c shows that the wind output from the proposed method was lower than from the traditional method, contributing to a higher MP, as shown in Figure 6b,d. The outputs from the proposed method and the traditional method were higher on the holiday. The proposed method could increase the MP from the traditional method by 26.81 $/MWh at 12 a.m. on the workday, and by 20.80 $/MWh at 6 a.m. on the holiday. The average MP during the day was driven up to around 8.8 $/MWh (23.6%) on the workday, and 7.9 $/MWh (30%) on the holiday. Table 6 shows the optimal VRE output levels during a day in the unit of percent of the maximum output that generators could provide. The average VRE outputs from the proposed method and traditional method were 90% and 97% on the workday and 78% and 88% on the holiday. That means, in the case where there was high VRE penetration (15 GW) and where wind generation was the main VRE proportion, i.e., case 6-and the wind generators received no support-curtailing approximately 7% of the wind output on the workday and 10% of the wind output on the holiday could provide the maximized profits of wind generators. Note that there was no solar penetration in case 6; thus, there was no solar generation in Table 6. or solar because the demand is high, but VRE output is low, contributing to high MP. Thus, maximized electricity production also provides maximized profits for VRE generators during that time.
The VRE optimal generation schedules and the MP of every case are presented in Tables 4 and 5. This paper presents a case where the proposed method provided greater profits than the traditional method, which is case 6, wherein VRE received no support. Figure 6 shows the optimal VRE generation schedules and the MP of that case.  In all cases, VRE were not the marginal units because the VRE penetration was not enough to serve the electricity demands. More VRE penetration was impossible because of the unsatisfiable SOCs. However, if the system's flexibility is improved and more VRE can be integrated, the MP would be very low if VRE were to be the marginal unit. The MP could be equal to the VRE MC, if VRE generators are not supported, or there could be negative support prices, if VRE generators are supported. In these cases, the proposed method will show more significant benefits than the traditional method. Moreover, VRE generators can gain more revenue by using energy storage. The curtailed VRE output from the proposed method can be stored and sold back to the system during the time when electricity demand is high. However, the costs incurred from using energy storage, such as installation costs, and costs from lost energy due to the efficiency of the energy storage, must be less than the revenue to avoid negative profits.
Acknowledgments: This paper was supported by The 100th Anniversary Chulalongkorn University Fund for Doctoral Scholarship. The data used in this paper was supported by the Electricity Generating Authority of Thailand (EGAT).

Conflicts of Interest:
The authors declare no conflict of interest.

CCGT
Combined Cycle Gas Turbines CFD-FIT Contract for Difference Feed-in Tariff The electricity demand at specific time t Piecewise costs from the incremental cost curve of thermal plant n ($/MWh) P R1 n,thermal , P R2 n,thermal , P R3 n,thermal Piecewise power from the incremental cost curve of thermal plant n (MW) Profile(t) Generation profile of VRE (%of installed capacity)