Multi-Energy Systems
MES · Sector Coupling
Electricity · Heat · Gas · Hydrogen · Mobility
Multi-Energy Systems (MES) jointly optimize electricity, district heating, natural gas, hydrogen, and increasingly cooling and mobility as a single coupled system. Sector coupling via power-to-X technologies (electrolyzers, heat pumps, synthetic fuels) is the defining infrastructure shift of the 2020s energy transition. Mathematically, MES is a multi-commodity network flow LP/MIP with each commodity's dispatch coupled through conversion technologies. Mancarella's Energy (2014) review and Lund et al.'s smart-energy-systems programme are the foundational references.
The problem
Integrated optimization across energy carriers
Traditional energy-system OR models each carrier (electricity, gas, heat) as a silo. But modern decarbonization requires tight coupling. A heat pump: electricity in, heat out. An electrolyzer: electricity in, hydrogen out. A CHP: gas in, electricity + heat out. A methanation plant: hydrogen + CO2 in, methane out. A P2L (power-to-liquid) plant: electricity + CO2 in, diesel or kerosene out. These power-to-X (P2X) technologies are the physical bridges between carriers; optimizing them jointly with each carrier's network is the multi-energy systems problem.
The key OR concept is the energy hub (Geidl & Andersson, 2007): a node that can convert among multiple inputs and outputs via a conversion matrix $\mathbf{C}$. If inputs are $\mathbf{P}_{\mathrm{in}}$ (vector of commodity flows in) and outputs $\mathbf{L}_{\mathrm{out}}$, then $\mathbf{L}_{\mathrm{out}} = \mathbf{C} \mathbf{P}_{\mathrm{in}}$ where $C_{ij}$ is the coupling coefficient from input $j$ to output $i$. This hub abstraction reduces messy sector-coupling interactions to linear-algebra bookkeeping.
Mathematical formulation
Multi-commodity network LP with conversion technologies
Energy hub modelGeidl-Andersson 2007
Each hub $h$ has vector of input commodity flows $\mathbf{P}_h \in \mathbb{R}^{|\mathcal{C}|}$ and output demand vector $\mathbf{L}_h$. Conversion matrix $\mathbf{C}_h$ (rows: output commodities, columns: input commodities) represents installed technology mix:
Each entry $C_{ij,h}$ = dispatch fraction of input $j$ at hub $h$ converted to output $i$, times efficiency. Row sums $\le 1$ (efficiency ≤ 100%); column sums represent input allocation.
Full MES LP
Objective: minimize total operational cost across carriers:
Commodity balance at each hub/network node:
Technology capacity: $A_{k,t} \le K_k$. Network flow: gas, heat, electricity, hydrogen networks each with their own DC-OPF-like constraints (for electricity) or pipeline-flow constraints (for gas/H2/heat).
Storage: each carrier has its own storage with SOC dynamics $e_{s,t+1} = e_{s,t} + \eta_s^c c_{s,t} - d_{s,t}/\eta_s^d$.
P2X technology examples
Heat pump: electricity → heat, COP ≈ 3 (3 kWh heat per 1 kWh electricity).
Electrolyzer: electricity → H2, efficiency ≈ 70% (HHV).
CHP: gas → electricity + heat, efficiencies ≈ 40% elec + 50% heat = 90% total.
Methanation: H2 + CO2 → CH4, efficiency ≈ 80%.
P2L: electricity + CO2 → liquid fuel, efficiency ≈ 40% (multi-step).
Complexity
Typical sector-coupled LP with 50 nodes, 5 carriers, 8760 hours, 20 technology types ≈ 10–100 million variables. Solved via Benders decomposition, time-slice aggregation, or Gurobi's distributed LP. PyPSA-Eur-Sec runs continental-scale cases in 4–12 hours on a workstation.
Solution interpretation
MES optimization reveals which P2X technologies are optimal at what scale. In high-wind regions: heat pumps dominate heating (electrify heat → absorb wind surplus); electrolyzers come online once wind curtailment is economic at current capex. In gas-rich regions: CHP and gas heating persist longer; H2 takes over for industrial heat only when carbon prices are high.
Cross-carrier price linkages emerge as shadow prices: electricity LMP + heat-pump COP = implicit district-heating floor price; electrolyzer capacity puts an upper bound on electricity prices during surplus hours. These linkages are the value of sector coupling — visible only in the joint model.
A consistent finding across MES studies: sector coupling reduces total system cost by 10–25% at 95%+ renewable shares compared to electricity-only optimization. The primary mechanism: heat pumps, electrolyzers, and thermal storage absorb renewable surpluses that would otherwise be curtailed.
Extensions & variants
Stochastic multi-energy MPC
Uncertain weather affects multiple carriers simultaneously (cold day = electricity + heat + gas demand spike). Joint scenario tree.
Integrated planning (MES-GEP)
Co-plan investments across all carriers. Key question: build new H2 pipeline or electrolyzer + electric grid reinforcement?
District-heating network optimization
Thermal storage + heat-pump scheduling with network temperatures. Physics adds nonlinearity but approximations keep LP-tractable.
Natural-gas network coupling
Gas-grid pressure/flow constraints coupling with electricity system. Gas-fired generation becomes network-constrained.
Cross-carrier ancillary services
Heat pumps + electrolyzers as flexible electricity load providing frequency regulation or reserve.
100% RE smart-energy-systems
Aalborg-style national planning with fully integrated sectors. Denmark 2050 studies led this agenda.