Merge pull request #99 from ESMAP-World-Bank-Group/feature/96-phase-out-max-power-output

Feature/96 phase out max power output

Author: lucas-vivier

.claude/commands/diagnose-model.md

@@ -27,6 +27,9 @@
 *.grib
 *.zip
 
+# Claude working files
+.claude/work/
+
 # IDE / Editor configs
 .vscode/
 !.vscode/launch.json.example
@@ -215,3 +218,4 @@ celerybeat.pid
 .dmypy.json
 dmypy.json
 .pyre/
+.claude/

.gitignore

@@ -203,6 +203,7 @@ Parameters
    fEnableCapacityExpansion
    fApplyRampConstraint           'Enable ramping constraints'
    fApplyFuelConstraint           'Enable fuel availability limits'
+   fApplyGenerationPhaseout       'Enable max generation by fuel constraint (phase-out)'
    fApplyCapitalConstraint        'Enable total capital constraint'
    fEnableCSP                     'Enable CSP features'
    fEnableStorage                 'Enable storage operations'
@@ -418,6 +419,7 @@ Equations
 * ------------------------------
    eFuel(z,f,y)                   'Fuel consumption accounting'
    eFuelLimit(c,f,y)              'Fuel availability constraint'
+   eMaxGenerationByFuel(z,tech,f,y) 'Maximum annual generation by zone-tech-fuel [GWh]'
 
 * ------------------------------
 * Ramp and reserve constraints
@@ -778,6 +780,10 @@ eFuel(zfmap(z,f),y)..
 eFuelLimit(c,f,y)$(fApplyFuelConstraint and pMaxFuelLimit(c,f,y) > 0)..
    sum((zcmap(z,c),zfmap(z,f)), vFuel(z,f,y)) =l= pMaxFuelLimit(c,f,y)*1e6;
 
+* Phase-out constraint: limits annual generation by zone-technology-fuel combination
+eMaxGenerationByFuel(z,tech,f,y)$(fApplyGenerationPhaseout and pMaxGenerationByFuel(z,tech,f,y))..
+   sum((gzmap(g,z),gtechmap(g,tech),gfmap(g,f),q,d,t), vPwrOut(g,f,q,d,t,y)*pHours(q,d,t))/1e3 =l= pMaxGenerationByFuel(z,tech,f,y);
+
 * Applies the minimum output requirement per capacity share.
 eMinGen(g,q,d,t,y)$((fApplyMinGenShareAllHours and pGenData(g,"MinGenShareAllHours") > 0) and FD(q,d,t))..
     sum(gfmap(g,f), vPwrOut(g,f,q,d,t,y)) =g= vCap(g,y)*pGenData(g,"MinGenShareAllHours");
@@ -1215,6 +1221,7 @@ Model PA /
 
    eVREProfile
    eFuelLimit
+   eMaxGenerationByFuel
    eCapitalConstraint
    eZonalEmissions
    eEmissionsCountry

epm/base.gms

@@ -33,6 +33,7 @@ System-wide spinning reserve requirement (Unit: MW),pSpinningReserveReqSystem,re
 System-level annual emissions cap (Unit: tCO2),pEmissionsTotal,constraint/pEmissionsTotal.csv
 Country-level annual emissions caps (Unit: tCO2),pEmissionsCountry,constraint/pEmissionsCountry.csv
 Maximum fuel consumption limits by zone and fuel (Unit: MMBtu),pMaxFuellimit,constraint/pMaxFuellimit.csv
+Maximum generation by zone-tech-fuel for phase-out scenarios (Unit: GWh),pMaxGenerationByFuel,constraint/pMaxGenerationByFuel.csv
 Carbon pricing trajectory (Unit: USD/tCO2),pCarbonPrice,constraint/pCarbonPrice.csv
 ,TRADE,
 External trading zones (Unit: unitless),zext,trade/zext.csv

epm/input/data_test/config.csv

@@ -0,0 +1,3 @@
+z,tech,f,2025,2040
+Angola,CCGT,Gas,100,0.001
+Angola,OCGT,Gas,100,0.001

epm/input/data_test/constraint/pMaxGenerationByFuel.csv

@@ -55,6 +55,7 @@ Apply system CO2 constraints,fApplySystemCo2Constraint,0
 Minimum share of RE  (%),sMinRenewableSharePct,0
 RE share target year,sRenewableTargetYear,
 Apply fuel constraints,fApplyFuelConstraint,0
+Apply generation phase-out constraint,fApplyGenerationPhaseout,0
 Apply max capital constraint,fApplyCapitalConstraint,0
 "Total maximum capital investments, $ billion",sMaxCapitalInvestment,
 ,,

epm/input/data_test/pSettings.csv

@@ -103,6 +103,7 @@ $if not set pCarbonPrice $set pCarbonPrice %FOLDER_INPUT%/constraint/pCarbonPric
 $if not set pEmissionsCountry $set pEmissionsCountry %FOLDER_INPUT%/constraint/pEmissionsCountry.csv
 $if not set pEmissionsTotal $set pEmissionsTotal %FOLDER_INPUT%/constraint/pEmissionsTotal.csv
 $if not set pMaxFuellimit $set pMaxFuellimit %FOLDER_INPUT%/constraint/pMaxFuellimit.csv
+$if not set pMaxGenerationByFuel $set pMaxGenerationByFuel %FOLDER_INPUT%/constraint/pMaxGenerationByFuel.csv
 
 * H2 RELATED
 $if not set pH2DataExcel $set pH2DataExcel %FOLDER_INPUT%/h2/pH2DataExcel.csv
@@ -610,6 +611,15 @@ $onEmbeddedCode Connect:
     header: [1]
     type: par
 
+- CSVReader:
+    trace: %TRACE%
+    file: %pMaxGenerationByFuel%
+    name: pMaxGenerationByFuel
+    indexSubstitutions: {.nan: ""}
+    indexColumns: [1, 2, 3]
+    header: [1]
+    type: par
+
 
 # H2 RELATED
 

epm/input_readers.gms

@@ -51,8 +51,9 @@
     'pDemandForecast',
     'pCapexTrajectories',
     'pExtTransferLimit',
-    'pTransferLimit', 
-    'pFuelPrice'
+    'pTransferLimit',
+    'pFuelPrice',
+    'pMaxGenerationByFuel'
 ]
 
 ZONE_RESTRICTED_PARAMS = {

epm/input_treatment.py

@@ -185,6 +185,7 @@ Parameter
 * Fuel data
    pFuelCarbonContent(f)                                 'Carbon content by fuel (tCO2/MMBtu)'
    pMaxFuellimit(c,f,y)                                  'Fuel limit in MMBTU*1e6 (million) by country'
+   pMaxGenerationByFuel(z,tech,f,y)                      'Max annual generation by zone-tech-fuel [GWh]'
    pFuelPrice(c,f,y)                                     'Fuel price forecasts'
 
 * Storage and transmission
@@ -266,7 +267,7 @@ $load pDemandData pDemandForecast pDemandProfile pEnergyEfficiencyFactor sReleva
 $load pFuelCarbonContent pCarbonPrice pEmissionsCountry pEmissionsTotal pFuelPrice
 
 * Load constraints and technical data
-$load pMaxFuellimit pTransferLimit pLossFactorInternal pVREProfile pVREgenProfile pAvailabilityInput pEvolutionAvailability
+$load pMaxFuellimit pMaxGenerationByFuel pTransferLimit pLossFactorInternal pVREProfile pVREgenProfile pAvailabilityInput pEvolutionAvailability
 * Use $loadM to merge storage units into set g (first dimension of pStorageDataInput)
 $loadM g<pStorageDataInput.Dim1
 $load pStorageDataInput pStorageDataInputDefault pStorageDataInputGeneric pCSPData pCapexTrajectories pSpinningReserveReqCountry pSpinningReserveReqSystem 
@@ -468,6 +469,7 @@ fCountIntercoForReserves           = pSettings("fCountIntercoForReserves");
 * --- Settings: Policy and operational switches
 fApplyMinGenShareAllHours      = pSettings("fApplyMinGenShareAllHours");
 fApplyFuelConstraint               = pSettings("fApplyFuelConstraint");
+fApplyGenerationPhaseout           = pSettings("fApplyGenerationPhaseout");
 fApplyCapitalConstraint            = pSettings("fApplyCapitalConstraint");
 fEnableCSP                         = pSettings("fEnableCSP");
 fEnableCapacityExpansion           = pSettings("fEnableCapacityExpansion");

epm/main.gms

@@ -37,6 +37,7 @@ fApplySystemCo2Constraint
 sMinRenewableSharePct
 sRenewableTargetYear
 fApplyFuelConstraint
+fApplyGenerationPhaseout
 fApplyCapitalConstraint
 sMaxCapitalInvestment
 fApplyMinGenShareAllHours