9.3 Energy Resources
The ocean is both a major source of fossil fuels (offshore oil and gas supply ~30% of global production) and a vast reservoir of renewable energy. Tidal power, wave energy, ocean thermal energy conversion (OTEC), offshore wind, salinity gradient power, and ocean current energy collectively represent an enormous untapped potential that is critical for the global energy transition.
Offshore Oil & Gas
Offshore hydrocarbon production occurs on continental shelves (shallow water, <200 m), continental slopes (deep water, 200–1500 m), and ultra-deep water (>1500 m). Technology has advanced from fixed steel platforms to floating production, storage, and offloading (FPSO) vessels and tension-leg platforms (TLPs).
~30%
Of global oil production is offshore
~25%
Of global gas production is offshore
>3000 m
Ultra-deep water drilling depth record
Environmental risks include blowouts (e.g., Deepwater Horizon 2010, 4.9 million barrels spilled), chronic operational discharges, seismic survey impacts on marine mammals, and greenhouse gas emissions from production and consumption. Decommissioning of aging platforms is a growing challenge.
Tidal Power
Tidal energy is extracted either from the potential energy of the tidal range (barrage systems) or from the kinetic energy of tidal currents (tidal stream turbines). The power available from a tidal stream turbine follows the same cubic velocity law as wind turbines:
$$P = \frac{1}{2} \rho A v^3 C_p$$
$\rho \approx 1025$ kg/m³ (seawater), $A$ = rotor swept area,$v$ = current speed, $C_p$ = power coefficient ($\leq 16/27$ Betz limit)
For a tidal barrage, the energy per tidal cycle from a basin of area $A_b$ with tidal range $R$ is:
$$E = \frac{1}{2} \rho g A_b R^2$$
Energy scales with the square of the tidal range
La Rance (France)
World's first tidal barrage (1966). 240 MW capacity. Tidal range ~8 m. Still operational after >55 years.
MeyGen (Scotland)
World's largest tidal stream array. Pentland Firth. Up to 398 MW planned. Peak currents >3.5 m/s.
Sihwa Lake (South Korea)
World's largest tidal barrage by capacity: 254 MW. Tidal range ~5.6 m. Completed 2011.
Bay of Fundy (Canada)
Highest tides on Earth (~16 m range). Multiple tidal stream projects in testing. Enormous energy potential.
Wave Power
Ocean waves carry energy transferred from the wind. The power per unit crest width for deep-water irregular waves characterized by significant wave height $H_s$ and energy period $T_e$ is:
$$P_w = \frac{\rho g^2}{64\pi} H_s^2 T_e \approx 0.49 \, H_s^2 T_e \;\;\text{(kW/m)}$$
Power scales as $H_s^2 T_e$; typical values: 20–80 kW/m in temperate latitudes
The capture width ratio $\eta = P_{\text{captured}} / (P_w \cdot D)$ characterizes device efficiency, where $D$ is the device width. Wave energy converter (WEC) types include point absorbers, attenuators (Pelamis), oscillating water columns (OWC), and overtopping devices. The global theoretical wave power resource is estimated at ~2 TW.
Ocean Thermal Energy Conversion (OTEC)
OTEC exploits the temperature difference between warm surface water ($\sim 25\text{--}28°C$ in the tropics) and cold deep water ($\sim 4\text{--}5°C$ at 1000 m). The maximum thermodynamic efficiency is limited by the Carnot efficiency:
$$\eta_{\text{Carnot}} = \frac{T_h - T_c}{T_h} = 1 - \frac{T_c}{T_h}$$
For $T_h = 300$ K and $T_c = 278$ K: $\eta_{\text{Carnot}} \approx 7.3\%$; practical efficiency ~3–5%
Despite low efficiency, the resource is enormous (~3–5 TW potential). OTEC can also produce desalinated water and support deep-water aquaculture with nutrient-rich cold water.
Offshore Wind & Salinity Gradient Power
Offshore Wind
The fastest-growing ocean energy sector with >60 GW installed globally. Fixed-bottom turbines dominate in shallow water (<60 m), while floating platforms (spar, semi-submersible, TLP) unlock deep-water sites. Turbine capacity now exceeds 15 MW.
$P_{\text{wind}} = \frac{1}{2} \rho_{\text{air}} A v_{\text{wind}}^3 C_p$
Salinity Gradient Power
Extracts energy from the mixing of freshwater and seawater at river mouths. Two main technologies: Pressure Retarded Osmosis (PRO) and Reverse Electrodialysis (RED). The theoretical mixing energy is:
$\Delta G_{\text{mix}} = -RT \sum c_i \ln(a_i) \approx 2.5$ MJ/m³
Levelized Cost of Energy (LCOE) Comparison
$$\text{LCOE} = \frac{\sum_{t=1}^{n} \frac{I_t + M_t + F_t}{(1+r)^t}}{\sum_{t=1}^{n} \frac{E_t}{(1+r)^t}}$$
$I_t$ = capital cost, $M_t$ = O&M, $F_t$ = fuel,$E_t$ = energy produced, $r$ = discount rate
Ocean Current Energy & Environmental Considerations
Major ocean currents like the Gulf Stream, Kuroshio, and Agulhas carry enormous kinetic energy. The power density of a current flowing at velocity $v$ is:
$$P_{\text{current}} = \frac{1}{2} \rho v^3 \quad (\text{W/m}^2)$$
Gulf Stream: $v \approx 1.5$ m/s ⇒ $P \approx 1700$ W/m²; vast but technically challenging to extract
Environmental Risks of Offshore Energy
Oil spills (chronic and catastrophic), seismic survey impacts on marine mammals, electromagnetic fields from submarine cables, collision risk for marine animals with turbine blades, and seabed disturbance during installation.
Artificial Reef Effect
Offshore wind foundations and tidal turbine structures act as artificial reefs, increasing local biodiversity. Studies show 2–5x increases in fish abundance around offshore wind monopiles compared to surrounding seabed.
Derivation: Tidal Power from a Barrage
Step 1: Setup — Potential Energy of a Water Column
Consider a tidal basin of surface area $A_b$ impounded behind a barrage. At high tide the water surface inside is elevated by a tidal range $R$ above the low-tide level outside. As the basin empties through turbines, a thin horizontal slab of water at height $z$ above the exterior level has mass:
$$dm = \rho \, A_b \, dz$$
Step 2: Potential Energy of Each Slab
Each slab falls through height $z$ (its center of mass drops from $z$ to the exterior level). Its potential energy is:
$$dE = dm \cdot g \cdot z = \rho g A_b z \, dz$$
Step 3: Integrate Over the Tidal Range
Integrate from $z = 0$ (low-tide level) to $z = R$ (high-tide level):
$$E = \int_0^R \rho g A_b z \, dz = \rho g A_b \left[\frac{z^2}{2}\right]_0^R = \frac{1}{2}\rho g A_b R^2$$
Step 4: Average Power Over a Tidal Cycle
A semidiurnal tidal cycle has period $T \approx 12.42$ hours ($\approx 44712$ s). The energy is released once per ebb (or twice for ebb-flood generation). Average power for single-effect operation:
$$\bar{P} = \frac{E}{T/2} = \frac{\rho g A_b R^2}{T}$$
Step 5: Practical Power with Turbine Efficiency
Including turbine-generator efficiency $\eta_t$ and accounting for two ebb tides per day:
$$P_{\text{annual}} = \eta_t \cdot \rho g A_b \overline{R^2} \cdot \frac{n_{\text{cycles}}}{t_{\text{year}}}$$
Since the tidal range varies with spring-neap cycles, $\overline{R^2}$ must be averaged over the lunar month. Typically $\overline{R^2} \approx 0.5 R_{\text{spring}}^2$.
Step 6: Numerical Example (La Rance)
For La Rance: $A_b = 22 \times 10^6$ m$^2$, $R = 8$ m, $\rho = 1025$ kg/m$^3$:
$$E = \frac{1}{2}(1025)(9.81)(22 \times 10^6)(8)^2 \approx 7.1 \times 10^{12} \text{ J} = 7.1 \text{ TJ per cycle}$$
Derivation: Wave Energy Flux in Deep Water
Step 1: Energy Density of a Monochromatic Wave
For a sinusoidal deep-water wave of amplitude $a$ and wave number $k$, the total energy per unit horizontal area (kinetic + potential) is:
$$E = \frac{1}{2}\rho g a^2$$
Step 2: Group Velocity in Deep Water
In deep water the dispersion relation is $\omega^2 = gk$. The phase speed is $c = \omega/k = g/\omega$ and the group velocity (energy propagation speed) is half the phase speed:
$$c_g = \frac{d\omega}{dk} = \frac{1}{2}\sqrt{\frac{g}{k}} = \frac{c}{2} = \frac{gT}{4\pi}$$
Step 3: Power Flux for a Regular Wave
The energy flux (power per unit crest width) is energy density times group velocity. Substituting $a = H/2$ (wave height $H = 2a$):
$$P = E \cdot c_g = \frac{1}{2}\rho g \left(\frac{H}{2}\right)^2 \cdot \frac{gT}{4\pi} = \frac{\rho g^2 H^2 T}{32\pi}$$
Step 4: Extend to Irregular Seas
For a wave spectrum $S(f)$, the significant wave height is $H_s = 4\sqrt{m_0}$ where $m_0 = \int S(f)\,df$. The energy period is $T_e = m_{-1}/m_0$. The power flux generalizes to:
$$P_w = \rho g \int_0^\infty c_g(f) S(f) \, df = \frac{\rho g^2}{64\pi} H_s^2 T_e$$
Step 5: Numerical Approximation
Substituting $\rho = 1025$ kg/m$^3$ and $g = 9.81$ m/s$^2$:
$$P_w = \frac{1025 \times 9.81^2}{64\pi} H_s^2 T_e \approx 490 \, H_s^2 T_e \;\; \text{W/m} \approx 0.49 \, H_s^2 T_e \;\; \text{kW/m}$$
For the North Atlantic with $H_s = 3$ m and $T_e = 9$ s: $P_w \approx 0.49 \times 9 \times 9 \approx 40$ kW/m.
Step 6: Capture Width and Device Efficiency
A wave energy converter of characteristic width $D$ captures power $P_{\text{cap}}$. The capture width ratio (hydrodynamic efficiency) is:
$$\eta = \frac{P_{\text{cap}}}{P_w \cdot D}, \quad \text{CW} = \eta D = \frac{P_{\text{cap}}}{P_w}$$
Point absorbers can have capture widths exceeding their physical size ($\eta > 1$) near resonance, analogous to the effective cross-section of an antenna.
Derivation: OTEC Carnot Efficiency
Step 1: The Second Law and Carnot's Theorem
A heat engine operating between a hot reservoir at temperature $T_h$ and a cold reservoir at $T_c$ (both in Kelvin) receives heat $Q_h$ and rejects $Q_c$. By the second law, for a reversible engine the entropy change of the universe is zero:
$$\Delta S_{\text{universe}} = -\frac{Q_h}{T_h} + \frac{Q_c}{T_c} = 0 \quad \Longrightarrow \quad \frac{Q_c}{Q_h} = \frac{T_c}{T_h}$$
Step 2: Derive Maximum Efficiency
The work output is $W = Q_h - Q_c$. The thermal efficiency is:
$$\eta_{\text{Carnot}} = \frac{W}{Q_h} = 1 - \frac{Q_c}{Q_h} = 1 - \frac{T_c}{T_h} = \frac{T_h - T_c}{T_h}$$
Step 3: Apply to OTEC Conditions
Tropical surface water: $T_h \approx 28°\text{C} = 301.15$ K. Deep water at 1000 m: $T_c \approx 4°\text{C} = 277.15$ K. The Carnot limit is:
$$\eta_{\text{Carnot}} = 1 - \frac{277.15}{301.15} = \frac{24}{301.15} \approx 7.97\%$$
Step 4: Practical OTEC Efficiency
Real OTEC systems achieve only 40–60% of the Carnot limit due to finite heat-exchanger temperature differences ($\Delta T_{\text{pinch}}$), pump parasitic loads (~30% of gross power), and working-fluid irreversibilities:
$$\eta_{\text{practical}} \approx (0.4\text{--}0.6) \times \eta_{\text{Carnot}} \approx 3\text{--}5\%$$
Step 5: Gross Power Output
The thermal power input is set by the warm-water mass flow rate $\dot{m}$, specific heat $c_p$, and the temperature drop across the evaporator $\Delta T_e$:
$$P_{\text{gross}} = \eta_{\text{practical}} \cdot \dot{m} c_p \Delta T_e$$
Step 6: Net Power and Parasitic Losses
Pumping cold water from 1000 m depth requires significant power. The net output is:
$$P_{\text{net}} = P_{\text{gross}} - P_{\text{pump,warm}} - P_{\text{pump,cold}} - P_{\text{aux}}$$
The cold-water pipe pumping power scales as $P_{\text{pump}} \propto \dot{m}_c g H / \eta_p$, where $H \approx 1000$ m is the pipe depth. Typically $P_{\text{net}} \approx 0.5 P_{\text{gross}}$, giving a net efficiency of $\sim 1.5\text{--}2.5\%$.
Python: Tidal, Wave, and OTEC Analysis
Python: Tidal, Wave, and OTEC Analysis
Python!/usr/bin/env python3
Click Run to execute the Python code
Code will be executed with Python 3 on the server
Fortran: Tidal Stream Turbine Power Output Model
This program models the power output of a tidal stream turbine over multiple tidal cycles, accounting for the sinusoidal variation of current speed with M2 and S2 tidal constituents and applying a cut-in speed and rated power limit.
Fortran: Tidal Stream Turbine Power Output Model
FortranTidal stream turbine power output with M2/S2 tidal variation
Click Run to execute the Fortran code
Code will be compiled with gfortran and executed on the server