Extreme Conditions Modeling - Full Sea State Approach

Extreme conditions modeling consists of identifying the expected extreme (e.g. 100-year) response of some quantity of interest, such as WEC motions or mooring loads. Three different methods of estimating extreme conditions were adapted from WDRT: full sea state approach, contour approach, and MLER design wave. This noteboook presents the full sea state approach.

The full sea state approach consists of the following steps: 1. Take \(N\) samples to represent the sea state. Each sample represents a small area of the sea state and consists of a representative \((H_{s}, T_{e})_i\) pair and a weight \(W_i\) associated with the probability of that sea state area. 2. For each sample \((H_{s}, T_{e})_i\) calculate the short-term (e.g. 3-hours) extreme for the quantity of interest (e.g. WEC motions or mooring tension). 3. Integrate over the entire sea state to obtain the long-term extreme. This is a sum of the products of the weight of each sea state times the short-term extreme.

See more details and equations in

[1] Coe, Ryan G., Carlos A. Michelén Ströfer, Aubrey Eckert-Gallup, and Cédric Sallaberry. 2018. “Full Long-Term Design Response Analysis of a Wave Energy Converter.” Renewable Energy 116: 356–66.

NOTE: Prior to running this example it is recommended to become familiar with environmental_contours_example.ipynb and short_term_extremes_example.ipynb since some code blocks are adapted from those examples and used here without the additional description.

We start by importing the relevant modules, including waves.contours submodule which includes the sampling function, and loads.extreme which inlcudes the short-term extreme and full sea state integration functions.

from mhkit.wave import resource, contours, graphics
from mhkit.loads import extreme
import matplotlib.pyplot as plt
from mhkit.wave.io import ndbc
import pandas as pd
import numpy as np

Obtain and Process NDBC Buouy Data

The first step will be obtaining the environmental data and creating the contours. See environmental_contours_example.ipynb for more details and explanations of how this is being done in the following code block.

parameter = "swden"
buoy_number = "46022"
ndbc_available_data = ndbc.available_data(parameter, buoy_number)

years_of_interest = ndbc_available_data[ndbc_available_data.year < 2013]

filenames = years_of_interest["filename"]
ndbc_requested_data = ndbc.request_data(parameter, filenames)

ndbc_data = {}
for year in ndbc_requested_data:
    year_data = ndbc_requested_data[year]
    ndbc_data[year] = ndbc.to_datetime_index(parameter, year_data)

Hm0_list = []
Te_list = []

# Iterate over each year and save the result in the initalized dictionary
for year in ndbc_data:
    year_data = ndbc_data[year]
    Hm0_list.append(resource.significant_wave_height(year_data.T))
    Te_list.append(resource.energy_period(year_data.T))

# Concatenate list of Series into a single DataFrame
Te = pd.concat(Te_list, axis=0)
Hm0 = pd.concat(Hm0_list, axis=0)
Hm0_Te = pd.concat([Hm0, Te], axis=1)

# Drop any NaNs created from the calculation of Hm0 or Te
Hm0_Te.dropna(inplace=True)
# Sort the DateTime index
Hm0_Te.sort_index(inplace=True)

Hm0_Te_clean = Hm0_Te[Hm0_Te.Hm0 < 20]

Hm0 = Hm0_Te_clean.Hm0.values
Te = Hm0_Te_clean.Te.values

dt = (Hm0_Te_clean.index[2] - Hm0_Te_clean.index[1]).seconds

1. Sampling

The first step is sampling the sea state to get samples \((H_s, T_e)_i\) and associtated weights. For this we will use the waves.contours.samples_full_seastate function. We will sample 20 points between each return level, for 10 levels ranging from 0.001—100 years return periods. For more details on the sampling approach see

[1] Coe, Ryan G., Carlos A. Michelén Ströfer, Aubrey Eckert-Gallup, and Cédric Sallaberry. 2018. “Full Long-Term Design Response Analysis of a Wave Energy Converter.” Renewable Energy 116: 356–66.

[2] Eckert-Gallup, Aubrey C., Cédric J. Sallaberry, Ann R. Dallman, and Vincent S. Neary. 2016. “Application of Principal Component Analysis (PCA) and Improved Joint Probability Distributions to the Inverse First-Order Reliability Method (I-FORM) for Predicting Extreme Sea States.” Ocean Engineering 112 (January): 307–19.

# return levels
levels = np.array([0.001, 0.01, 0.05, 0.1, 0.5, 1, 5, 10, 50, 100])

# points per return level interval
npoints = 20

# Create samples
sample_hs, sample_te, sample_weights = contours.samples_full_seastate(
    Hm0, Te, npoints, levels, dt
)

We will now plot the samples alongside the contours. First we will create the different contours using contours.environmental_contours. See environmental_contours_example.ipynb for more details on using this function. There are 20 samples, randomly distributed, between each set of return levels.

# Create the contours
Te_contours = []
Hm0_contours = []

for period in levels:
    copula = contours.environmental_contours(
        Hm0, Te, dt, period, "PCA", return_PCA=True
    )
    Hm0_contours.append(copula["PCA_x1"])
    Te_contours.append(copula["PCA_x2"])

# plot
fig, ax = plt.subplots(figsize=(8, 4))

labels = [f"{period}-year Contour" for period in levels]

ax = graphics.plot_environmental_contour(
    sample_te,
    sample_hs,
    Te_contours,
    Hm0_contours,
    data_label="Samples",
    contour_label=labels,
    x_label="Energy Period, $Te$ [s]",
    y_label="Sig. wave height, $Hm0$ [m]",
    ax=ax,
)

2. Short-Term Extreme Distributions

Many different methods for short-term extremes were adapted from WDRT, and a summary and examples can be found in short_term_extremes_example.ipynb. The response quantity of interest is typically related to the WEC itself, e.g. maximum heave displacement, PTO extension, or load on the mooring lines. This requires running a simulation (e.g. WEC-Sim) for each of the 200 sampled sea states \((H_s, T_e)_i\). For the sake of example we will consider the wave elevation as the quantity of interest (can be thought as a proxy for heave motion in this example). Wave elevation time-series for a specific sea state can be created quickly without running any external software.

NOTE: The majority of the for loop below is simply creating the synthetic data (wave elevation time series). In a realistic case the variables time and data describing each time series would be obtained externally, e.g. through simulation software such as WEC-Sim or CFD. For this reason the details of creating the synthetic data are not presented here, instead assume for each sea state there is time-series data available.

The last lines of the for-loop create the short-term extreme distribution from the time-series using the loads.extreme.short_term_extreme function. The short-term period will be 3-hours and we will use 1-hour “simulations” and the Weibul-tail-fitting method to estimate the 3-hour short-term extreme distributions for each of the 200 samples.

For more details on short-term extreme distributions see short_term_extremes_example.ipynb and

[3] Michelén Ströfer, Carlos A., and Ryan Coe. 2015. “Comparison of Methods for Estimating Short-Term Extreme Response of Wave Energy Converters.” In OCEANS 2015 - MTS/IEEE Washington, 1–6. IEEE.

# create the short-term extreme distribution for each sample sea state
t_st = 3.0 * 60.0 * 60.0
gamma = 3.3
t_sim = 1.0 * 60.0 * 60.0

ste_all = []
i = 0
n = len(sample_hs)
for hs, te in zip(sample_hs, sample_te):
    tp = te / (0.8255 + 0.03852 * gamma - 0.005537 * gamma**2 + 0.0003154 * gamma**3)
    i += 1
    print(f"Sea state {i}/{n}. (Hs, Te) = ({hs} m, {te} s). Tp = {tp} s")
    # time & frequency arrays
    df = 1.0 / t_sim
    T_min = tp / 10.0  # s
    f_max = 1.0 / T_min
    Nf = int(f_max / df) + 1
    time = np.linspace(0, t_sim, 2 * Nf + 1)
    f = np.linspace(0.0, f_max, Nf)
    # spectrum
    S = resource.jonswap_spectrum(f, tp, hs, gamma)
    # 1-hour elevation time-series
    data = resource.surface_elevation(S, time).values.squeeze()
    # 3-hour extreme distribution
    ste = extreme.short_term_extreme(time, data, t_st, "peaks_weibull_tail_fit")
    ste_all.append(ste)

Sea state 1/200. (Hs, Te) = (2.851166377141435 m, 10.195736685047796 s). Tp = 11.282805679200564 s
Sea state 2/200. (Hs, Te) = (3.5862460438066304 m, 10.936985799409186 s). Tp = 12.103086741331571 s
Sea state 3/200. (Hs, Te) = (2.7751615384993915 m, 9.506986218875468 s). Tp = 10.520620668805236 s
Sea state 4/200. (Hs, Te) = (3.0820450802745007 m, 9.728717644182273 s). Tp = 10.765993088865644 s
Sea state 5/200. (Hs, Te) = (2.7162732683836195 m, 9.177186568411022 s). Tp = 10.155657794203384 s
Sea state 6/200. (Hs, Te) = (2.7916623630702304 m, 8.810734371427328 s). Tp = 9.75013448019444 s
Sea state 7/200. (Hs, Te) = (2.592821066032464 m, 8.374599835754573 s). Tp = 9.267499299629163 s
Sea state 8/200. (Hs, Te) = (2.4021986215575373 m, 8.348365892255037 s). Tp = 9.238468294234675 s
Sea state 9/200. (Hs, Te) = (2.327723379324664 m, 8.299620704291481 s). Tp = 9.184525896487639 s
Sea state 10/200. (Hs, Te) = (2.00764351315584 m, 7.416472871792525 s). Tp = 8.207216881170892 s
Sea state 11/200. (Hs, Te) = (1.951087921529932 m, 8.122015735790558 s). Tp = 8.987984694104963 s
Sea state 12/200. (Hs, Te) = (1.642063924746251 m, 7.997394534817154 s). Tp = 8.85007638619882 s
Sea state 13/200. (Hs, Te) = (1.5614901176199742 m, 8.156813011698352 s). Tp = 9.02649205403011 s
Sea state 14/200. (Hs, Te) = (1.3223442301926103 m, 8.373988391532727 s). Tp = 9.266822663251443 s
Sea state 15/200. (Hs, Te) = (2.310517824251443 m, 9.135556192896424 s). Tp = 10.109588789893749 s
Sea state 16/200. (Hs, Te) = (1.546312740170234 m, 9.978412702019654 s). Tp = 11.042310622719521 s
Sea state 17/200. (Hs, Te) = (1.9019087805024188 m, 10.44123868584584 s). Tp = 11.554483092459133 s
Sea state 18/200. (Hs, Te) = (2.3338262351351116 m, 9.930633820011005 s). Tp = 10.989437558425541 s
Sea state 19/200. (Hs, Te) = (2.4479576516446064 m, 9.760511705407653 s). Tp = 10.801177031491845 s
Sea state 20/200. (Hs, Te) = (2.651578075055819 m, 10.749317530563182 s). Tp = 11.895409289966368 s
Sea state 21/200. (Hs, Te) = (5.030844895631051 m, 14.282752593036674 s). Tp = 15.805579042430452 s
Sea state 22/200. (Hs, Te) = (4.463644749835009 m, 12.34746631098706 s). Tp = 13.663952622633797 s
Sea state 23/200. (Hs, Te) = (4.903918051881774 m, 11.272595420034712 s). Tp = 12.474479044855961 s
Sea state 24/200. (Hs, Te) = (4.37807740903193 m, 9.956362967463496 s). Tp = 11.017909946440787 s
Sea state 25/200. (Hs, Te) = (3.9498715792778443 m, 8.856270002839082 s). Tp = 9.800525118612159 s
Sea state 26/200. (Hs, Te) = (4.070115900693828 m, 8.307318836624022 s). Tp = 9.193044803348798 s
Sea state 27/200. (Hs, Te) = (2.736299341121807 m, 7.4381867053480155 s). Tp = 8.231245842699135 s
Sea state 28/200. (Hs, Te) = (2.4118215576759034 m, 6.537612647601825 s). Tp = 7.234652618770482 s
Sea state 29/200. (Hs, Te) = (1.9552243695320324 m, 6.321440801925461 s). Tp = 6.995432540474677 s
Sea state 30/200. (Hs, Te) = (1.801478730354882 m, 6.510603277357231 s). Tp = 7.204763510665695 s
Sea state 31/200. (Hs, Te) = (1.3679970074652714 m, 6.468135086563519 s). Tp = 7.157767363248965 s
Sea state 32/200. (Hs, Te) = (1.0743335950058275 m, 7.090922017059409 s). Tp = 7.846955808713215 s
Sea state 33/200. (Hs, Te) = (1.2357391850105406 m, 7.849764893906227 s). Tp = 8.686706479507244 s
Sea state 34/200. (Hs, Te) = (0.9213431748668424 m, 8.443600703410047 s). Tp = 9.343857036740465 s
Sea state 35/200. (Hs, Te) = (0.4814455431778253 m, 9.340126880987539 s). Tp = 10.33597079569597 s
Sea state 36/200. (Hs, Te) = (1.0933659724162896 m, 9.979135574932235 s). Tp = 11.043110568310052 s
Sea state 37/200. (Hs, Te) = (1.0123053762369256 m, 11.850902621257452 s). Tp = 13.114445334280386 s
Sea state 38/200. (Hs, Te) = (1.7075134960575675 m, 11.578239068080656 s). Tp = 12.812710405130463 s
Sea state 39/200. (Hs, Te) = (2.518274653405948 m, 12.53479278415886 s). Tp = 13.871251836085138 s
Sea state 40/200. (Hs, Te) = (3.529852197639782 m, 13.30774814454997 s). Tp = 14.72661966279441 s
Sea state 41/200. (Hs, Te) = (5.148223472776434 m, 15.586005564766051 s). Tp = 17.247784788330833 s
Sea state 42/200. (Hs, Te) = (5.556129509905438 m, 13.785131436466855 s). Tp = 15.254901540169124 s
Sea state 43/200. (Hs, Te) = (6.576702748036527 m, 13.803319512494747 s). Tp = 15.275028828058206 s
Sea state 44/200. (Hs, Te) = (6.326628174866466 m, 11.194121656530978 s). Tp = 12.387638412160275 s
Sea state 45/200. (Hs, Te) = (5.4962832147142935 m, 9.623068260891554 s). Tp = 10.649079373003927 s
Sea state 46/200. (Hs, Te) = (3.8010066164603415 m, 7.1969594306047 s). Tp = 7.964298926598074 s
Sea state 47/200. (Hs, Te) = (3.376615764959727 m, 7.194776215525951 s). Tp = 7.961882937224246 s
Sea state 48/200. (Hs, Te) = (2.6216876932404376 m, 6.227995589253622 s). Tp = 6.892024203363135 s
Sea state 49/200. (Hs, Te) = (1.6315497369764984 m, 5.243725558813215 s). Tp = 5.80281134583553 s
Sea state 50/200. (Hs, Te) = (1.3669420800943772 m, 5.356634348657379 s). Tp = 5.92775846585616 s
Sea state 51/200. (Hs, Te) = (0.8083067025753863 m, 5.521508218837552 s). Tp = 6.110211180778536 s
Sea state 52/200. (Hs, Te) = (0.6856869087958939 m, 5.632227381764594 s). Tp = 6.232735215957039 s
Sea state 53/200. (Hs, Te) = (0.40402171887623717 m, 6.179589350227688 s). Tp = 6.838456893274449 s
Sea state 54/200. (Hs, Te) = (0.18050092680749552 m, 7.8131070731262104 s). Tp = 8.646140203497934 s
Sea state 55/200. (Hs, Te) = (0.09436285952322665 m, 9.01492519221614 s). Tp = 9.976096117259395 s
Sea state 56/200. (Hs, Te) = (0.47347906604424894 m, 10.821566270980709 s). Tp = 11.975361187889881 s
Sea state 57/200. (Hs, Te) = (0.527437042017636 m, 13.504929376019787 s). Tp = 14.944824348437386 s
Sea state 58/200. (Hs, Te) = (1.7419771172383736 m, 14.110799728601373 s). Tp = 15.615292571198744 s
Sea state 59/200. (Hs, Te) = (2.4488493325170952 m, 14.752553677478932 s). Tp = 16.32547029770514 s
Sea state 60/200. (Hs, Te) = (4.087706649607579 m, 16.122641924784478 s). Tp = 17.841637293305862 s
Sea state 61/200. (Hs, Te) = (0.33176145271261115 m, 13.75930528117036 s). Tp = 15.22632180133784 s
Sea state 62/200. (Hs, Te) = (0.8855174001154404 m, 15.076157859092664 s). Tp = 16.683577142841592 s
Sea state 63/200. (Hs, Te) = (2.5427157275613133 m, 16.843426458663085 s). Tp = 18.63927184228866 s
Sea state 64/200. (Hs, Te) = (4.337377283702386 m, 17.077378007829093 s). Tp = 18.898167295273343 s
Sea state 65/200. (Hs, Te) = (4.836290692740213 m, 16.880304961749243 s). Tp = 18.680082329741825 s
Sea state 66/200. (Hs, Te) = (7.138210033996883 m, 15.883932531728718 s). Tp = 17.577476715326455 s
Sea state 67/200. (Hs, Te) = (7.04607811580354 m, 13.695319484323035 s). Tp = 15.15551384166231 s
Sea state 68/200. (Hs, Te) = (6.683515089328116 m, 11.949932510664025 s). Tp = 13.22403378611297 s
Sea state 69/200. (Hs, Te) = (6.127794684752878 m, 10.242888484924617 s). Tp = 11.334984801893572 s
Sea state 70/200. (Hs, Te) = (6.0852360496095255 m, 10.21048065086315 s). Tp = 11.299121645998666 s
Sea state 71/200. (Hs, Te) = (4.284337177258922 m, 7.5594798415804405 s). Tp = 8.365471247748896 s
Sea state 72/200. (Hs, Te) = (3.9700754639265234 m, 7.320993637781578 s). Tp = 8.101557655455082 s
Sea state 73/200. (Hs, Te) = (2.4674726025092806 m, 5.567536470063029 s). Tp = 6.161146962112716 s
Sea state 74/200. (Hs, Te) = (2.4411565850268797 m, 5.576064022637772 s). Tp = 6.170583721965596 s
Sea state 75/200. (Hs, Te) = (1.5027020575734593 m, 4.974528663747209 s). Tp = 5.504912689730755 s
Sea state 76/200. (Hs, Te) = (1.371060685104678 m, 4.955398811863336 s). Tp = 5.483743213885607 s
Sea state 77/200. (Hs, Te) = (0.9867265408067029 m, 5.015636588492398 s). Tp = 5.550403539593066 s
Sea state 78/200. (Hs, Te) = (0.4019223940784351 m, 5.424495697652933 s). Tp = 6.002855188878451 s
Sea state 79/200. (Hs, Te) = (0.3041167158650624 m, 6.175974380713595 s). Tp = 6.834456495870742 s
Sea state 80/200. (Hs, Te) = (0.04419855678833473 m, 7.035319486906778 s). Tp = 7.785424939256372 s
Sea state 81/200. (Hs, Te) = (1.1406640868177353 m, 16.99689554825793 s). Tp = 18.80910379942441 s
Sea state 82/200. (Hs, Te) = (2.088875590556972 m, 18.853554386717196 s). Tp = 20.86371952107473 s
Sea state 83/200. (Hs, Te) = (3.3670302274138626 m, 19.204362993751015 s). Tp = 21.251931326265645 s

Sea state 84/200. (Hs, Te) = (5.125297529462439 m, 18.672283211728097 s). Tp = 20.66312122143036 s
Sea state 85/200. (Hs, Te) = (7.5395641905378215 m, 17.827501744997704 s). Tp = 19.728269192102534 s
Sea state 86/200. (Hs, Te) = (7.077565789406247 m, 16.565463530565825 s). Tp = 18.331672519098774 s
Sea state 87/200. (Hs, Te) = (7.660797463748704 m, 14.340144866170789 s). Tp = 15.86909047718626 s
Sea state 88/200. (Hs, Te) = (8.267146386057341 m, 13.752524476289345 s). Tp = 15.218818027340491 s
Sea state 89/200. (Hs, Te) = (7.5585377867557515 m, 11.41873618959678 s). Tp = 12.636201336801353 s
Sea state 90/200. (Hs, Te) = (6.189724705513505 m, 9.943266929556025 s). Tp = 11.003417609551178 s
Sea state 91/200. (Hs, Te) = (4.985820233723228 m, 8.101640035422118 s). Tp = 8.965436537464955 s
Sea state 92/200. (Hs, Te) = (3.950449104550366 m, 7.032331054963992 s). Tp = 7.7821178808320965 s
Sea state 93/200. (Hs, Te) = (3.349526273174492 m, 6.2239567755115734 s). Tp = 6.887554771478645 s
Sea state 94/200. (Hs, Te) = (2.5305938620600323 m, 5.471413955273054 s). Tp = 6.054775869049686 s
Sea state 95/200. (Hs, Te) = (1.864284800127662 m, 4.692418833126264 s). Tp = 5.19272432145363 s
Sea state 96/200. (Hs, Te) = (1.3594583232962816 m, 4.600545483005879 s). Tp = 5.091055438809216 s
Sea state 97/200. (Hs, Te) = (1.0278889995696616 m, 4.580250801273515 s). Tp = 5.068596939008678 s
Sea state 98/200. (Hs, Te) = (0.7073568060165395 m, 4.792230763150695 s). Tp = 5.30317819504014 s
Sea state 99/200. (Hs, Te) = (0.508628205353712 m, 4.912581414426901 s). Tp = 5.43636062742936 s
Sea state 100/200. (Hs, Te) = (0.1887138929931016 m, 5.815845769563984 s). Tp = 6.435930988137664 s
Sea state 101/200. (Hs, Te) = (0.1871884572595377 m, 17.031553776651883 s). Tp = 18.847457286595947 s
Sea state 102/200. (Hs, Te) = (1.8255533228167624 m, 19.044007965006898 s). Tp = 21.07447925145326 s
Sea state 103/200. (Hs, Te) = (4.085461042668986 m, 19.625018339020407 s). Tp = 21.717436925831777 s
Sea state 104/200. (Hs, Te) = (5.26393200451075 m, 19.6875070497272 s). Tp = 21.78658818010872 s
Sea state 105/200. (Hs, Te) = (7.7652549552321455 m, 18.932000395608686 s). Tp = 20.950529439962647 s
Sea state 106/200. (Hs, Te) = (8.242537266166547 m, 17.662905419000868 s). Tp = 19.54612361311204 s
Sea state 107/200. (Hs, Te) = (9.27591707540523 m, 15.26409450321619 s). Tp = 16.89155158364453 s
Sea state 108/200. (Hs, Te) = (8.866267940605884 m, 13.700480521607233 s). Tp = 15.161225148513351 s
Sea state 109/200. (Hs, Te) = (7.32630515058675 m, 10.931360540476366 s). Tp = 12.096861717540277 s
Sea state 110/200. (Hs, Te) = (6.499223554458602 m, 9.489259188129742 s). Tp = 10.501003582825902 s
Sea state 111/200. (Hs, Te) = (4.677064608710481 m, 7.3713129233161805 s). Tp = 8.15724198098653 s
Sea state 112/200. (Hs, Te) = (3.760226861005709 m, 6.226071145112327 s). Tp = 6.889894575072661 s
Sea state 113/200. (Hs, Te) = (2.9941197640583335 m, 5.514672741901748 s). Tp = 6.102646905594336 s
Sea state 114/200. (Hs, Te) = (2.5717123084913123 m, 5.108974693146778 s). Tp = 5.653693348835068 s
Sea state 115/200. (Hs, Te) = (1.8783735128639163 m, 4.539471617179031 s). Tp = 5.023469880110703 s
Sea state 116/200. (Hs, Te) = (1.2752130858691124 m, 4.242007181700333 s). Tp = 4.694289799685399 s
Sea state 117/200. (Hs, Te) = (0.9793803265738974 m, 4.2470174592346845 s). Tp = 4.6998342727887366 s
Sea state 118/200. (Hs, Te) = (0.7124811539135159 m, 4.424785222214989 s). Tp = 4.89655562679098 s
Sea state 119/200. (Hs, Te) = (0.3588409389431 m, 4.705576661939675 s). Tp = 5.207285037392834 s
Sea state 120/200. (Hs, Te) = (0.18838257732596975 m, 4.92177034404479 s). Tp = 5.446529280316457 s
Sea state 121/200. (Hs, Te) = (1.1238727812408893 m, 19.89593130663322 s). Tp = 22.017234622062354 s
Sea state 122/200. (Hs, Te) = (3.256336956058525 m, 22.014360984983714 s). Tp = 24.361531178968573 s
Sea state 123/200. (Hs, Te) = (4.659638373417929 m, 20.55339693581869 s). Tp = 22.74479920753578 s
Sea state 124/200. (Hs, Te) = (8.409675325429806 m, 20.5323227793156 s). Tp = 22.721478125398956 s
Sea state 125/200. (Hs, Te) = (9.29875873352201 m, 19.828872472868778 s). Tp = 21.9430259784096 s
Sea state 126/200. (Hs, Te) = (9.628455421584839 m, 16.805889531372728 s). Tp = 18.597732729471783 s
Sea state 127/200. (Hs, Te) = (10.40162875212416 m, 15.193034921745538 s). Tp = 16.812915632742182 s
Sea state 128/200. (Hs, Te) = (9.514119089597632 m, 13.496075929567501 s). Tp = 14.935026947916162 s
Sea state 129/200. (Hs, Te) = (8.425241115023425 m, 11.273431425724304 s). Tp = 12.475404185435181 s
Sea state 130/200. (Hs, Te) = (7.181441955239233 m, 9.927469647992638 s). Tp = 10.985936022526603 s
Sea state 131/200. (Hs, Te) = (5.875841398493123 m, 8.165254125537952 s). Tp = 9.035833159053732 s
Sea state 132/200. (Hs, Te) = (4.792485900332849 m, 6.9990505762219986 s). Tp = 7.745289038919442 s
Sea state 133/200. (Hs, Te) = (3.1248928086261785 m, 5.482392746553494 s). Tp = 6.066925218570586 s
Sea state 134/200. (Hs, Te) = (2.6573548488075835 m, 5.066441779515841 s). Tp = 5.606625577835072 s
Sea state 135/200. (Hs, Te) = (1.7636422458226333 m, 4.351913064671034 s). Tp = 4.8159138430346085 s
Sea state 136/200. (Hs, Te) = (1.5014369679854003 m, 4.266948107278621 s). Tp = 4.721889925645097 s
Sea state 137/200. (Hs, Te) = (0.9816361473032175 m, 4.15497551695166 s). Tp = 4.597978822692763 s
Sea state 138/200. (Hs, Te) = (0.8127603228313202 m, 4.141264736044045 s). Tp = 4.582806198931542 s
Sea state 139/200. (Hs, Te) = (0.4006290776350021 m, 4.287345227405981 s). Tp = 4.744461785817772 s
Sea state 140/200. (Hs, Te) = (0.18188378189464394 m, 4.611357261437197 s). Tp = 5.103019968578412 s
Sea state 141/200. (Hs, Te) = (0.2805216324632074 m, 20.448028838143248 s). Tp = 22.62819676794919 s
Sea state 142/200. (Hs, Te) = (3.4488490911962097 m, 22.44173356304588 s). Tp = 24.834470202390577 s
Sea state 143/200. (Hs, Te) = (6.128032867857982 m, 22.368550725950417 s). Tp = 24.753484614442662 s
Sea state 144/200. (Hs, Te) = (8.0257433252404 m, 21.09485355189949 s). Tp = 23.343985904061235 s
Sea state 145/200. (Hs, Te) = (9.624239654696968 m, 20.397290291584074 s). Tp = 22.57204847980129 s
Sea state 146/200. (Hs, Te) = (10.880424233644302 m, 18.09569728837873 s). Tp = 20.025059732449993 s
Sea state 147/200. (Hs, Te) = (10.980286858713281 m, 15.648662372262576 s). Tp = 17.317122061869792 s
Sea state 148/200. (Hs, Te) = (9.626205932332264 m, 13.035460355521655 s). Tp = 14.425300561749852 s
Sea state 149/200. (Hs, Te) = (8.48145364001674 m, 10.969828551388684 s). Tp = 12.139431186073685 s
Sea state 150/200. (Hs, Te) = (7.496113729159182 m, 9.77133103334069 s). Tp = 10.813149915690298 s
Sea state 151/200. (Hs, Te) = (5.616718937684953 m, 7.69244938169326 s). Tp = 8.512617994344046 s
Sea state 152/200. (Hs, Te) = (4.8696012734631715 m, 6.889172832566847 s). Tp = 7.623696135098361 s
Sea state 153/200. (Hs, Te) = (3.1127279534654217 m, 5.146140131826839 s). Tp = 5.6948213649424595 s
Sea state 154/200. (Hs, Te) = (2.3851224952288668 m, 4.573133253276174 s). Tp = 5.060720518358844 s
Sea state 155/200. (Hs, Te) = (2.1405588081566735 m, 4.408339089583258 s). Tp = 4.87835600822587 s
Sea state 156/200. (Hs, Te) = (1.4376356173488731 m, 3.9651074597519202 s). Tp = 4.3878670416369 s
Sea state 157/200. (Hs, Te) = (1.0717947445845732 m, 3.797209276299676 s). Tp = 4.202067561295225 s
Sea state 158/200. (Hs, Te) = (0.7714659563414245 m, 3.8906591392370555 s). Tp = 4.305481044483412 s
Sea state 159/200. (Hs, Te) = (0.4675473016122803 m, 4.01891894221551 s). Tp = 4.447415901979359 s
Sea state 160/200. (Hs, Te) = (0.22912012469457133 m, 4.117183864360642 s). Tp = 4.556157825774942 s
Sea state 161/200. (Hs, Te) = (0.8868405022143149 m, 22.262946721786285 s). Tp = 24.636621116028003 s
Sea state 162/200. (Hs, Te) = (3.3435032701311083 m, 23.466210376684526 s). Tp = 25.968176688561687 s
Sea state 163/200. (Hs, Te) = (6.767304773171118 m, 23.562164689248302 s). Tp = 26.074361645684604 s
Sea state 164/200. (Hs, Te) = (10.473453427011934 m, 21.667662928997025 s). Tp = 23.977868191595526 s
Sea state 165/200. (Hs, Te) = (10.16869748664882 m, 20.88428094047033 s). Tp = 23.11096205012028 s
Sea state 166/200. (Hs, Te) = (12.592753558411193 m, 18.630634795246163 s). Tp = 20.61703225262625 s

Sea state 167/200. (Hs, Te) = (11.673289407430126 m, 16.196134650574198 s). Tp = 17.9229658102669 s
Sea state 168/200. (Hs, Te) = (10.377818028562164 m, 12.9398233046888 s). Tp = 14.319466692789565 s
Sea state 169/200. (Hs, Te) = (8.744961663025443 m, 10.648293741704611 s). Tp = 11.783614349210453 s
Sea state 170/200. (Hs, Te) = (6.859939603914331 m, 8.737652031185462 s). Tp = 9.669260104070206 s
Sea state 171/200. (Hs, Te) = (5.960600437241585 m, 8.015473102288261 s). Tp = 8.87008248973364 s
Sea state 172/200. (Hs, Te) = (4.693969291807187 m, 6.476150181500641 s). Tp = 7.166637027253579 s
Sea state 173/200. (Hs, Te) = (3.1828656644648836 m, 5.082363302518114 s). Tp = 5.624244652915248 s
Sea state 174/200. (Hs, Te) = (2.7310116880872317 m, 4.561878698618447 s). Tp = 5.048266003581575 s
Sea state 175/200. (Hs, Te) = (1.960761965106161 m, 3.9966122086956775 s). Tp = 4.422730825518815 s
Sea state 176/200. (Hs, Te) = (1.745663077932535 m, 3.8392198892428673 s). Tp = 4.248557343931261 s
Sea state 177/200. (Hs, Te) = (1.2054393707718805 m, 3.807490373589555 s). Tp = 4.213444828723166 s
Sea state 178/200. (Hs, Te) = (0.7359796530337633 m, 3.6471354786545067 s). Tp = 4.035992901986565 s
Sea state 179/200. (Hs, Te) = (0.6244103771848238 m, 3.805058507193115 s). Tp = 4.2107536768391585 s
Sea state 180/200. (Hs, Te) = (0.1947445590996868 m, 4.049025575814537 s). Tp = 4.48073250352036 s
Sea state 181/200. (Hs, Te) = (2.1148904507720747 m, 24.527579096070692 s). Tp = 27.142708493838484 s
Sea state 182/200. (Hs, Te) = (5.047569685078679 m, 24.918912509361533 s). Tp = 27.575765858336137 s
Sea state 183/200. (Hs, Te) = (6.420147613959171 m, 24.768168426047513 s). Tp = 27.40894944753028 s
Sea state 184/200. (Hs, Te) = (9.266139293031646 m, 23.421568304895782 s). Tp = 25.91877488105444 s
Sea state 185/200. (Hs, Te) = (12.41274872958178 m, 20.666780023774187 s). Tp = 22.870271195311382 s
Sea state 186/200. (Hs, Te) = (12.957910005068928 m, 18.230023846708807 s). Tp = 20.17370818248467 s
Sea state 187/200. (Hs, Te) = (12.2843220215601 m, 15.409966126505479 s). Tp = 17.052976032953453 s
Sea state 188/200. (Hs, Te) = (10.519991169667653 m, 12.528633985612458 s). Tp = 13.864436386710274 s
Sea state 189/200. (Hs, Te) = (9.450056877314221 m, 11.348957852297204 s). Tp = 12.5589832362669 s
Sea state 190/200. (Hs, Te) = (8.160919020315601 m, 9.81433166497037 s). Tp = 10.860735272761364 s
Sea state 191/200. (Hs, Te) = (5.587707137405241 m, 7.232570607296176 s). Tp = 8.00370696537244 s
Sea state 192/200. (Hs, Te) = (4.8308252936215474 m, 6.3910826892692425 s). Tp = 7.072499642765086 s
Sea state 193/200. (Hs, Te) = (3.9608488792601535 m, 5.641529187160236 s). Tp = 6.2430287794323185 s
Sea state 194/200. (Hs, Te) = (2.8150456982143615 m, 4.58212563012944 s). Tp = 5.070671662034252 s
Sea state 195/200. (Hs, Te) = (2.214607067395854 m, 4.085048800640139 s). Tp = 4.520596522982679 s
Sea state 196/200. (Hs, Te) = (1.649511190726026 m, 3.713958023766156 s). Tp = 4.1099400639008214 s
Sea state 197/200. (Hs, Te) = (0.9652307901744451 m, 3.456495828266831 s). Tp = 3.8250272605112143 s
Sea state 198/200. (Hs, Te) = (0.6126122310648106 m, 3.484511462758769 s). Tp = 3.8560299207128654 s
Sea state 199/200. (Hs, Te) = (0.468987733159132 m, 3.55586551038935 s). Tp = 3.9349917337440825 s
Sea state 200/200. (Hs, Te) = (0.10978529871712417 m, 3.772910684185013 s). Tp = 4.175178254364618 s

3. Long-Term Extreme Distribution

Finally we integrate the weighted short-term extreme distributions over the entire sea state space to obtain the extreme distribution, assuming a 3-hour sea state coherence. For this we use the loads.extreme.full_seastate_long_term_extreme function. The integral reduces to a sum over the 200 bins of the weighted short-term extreme distributions.

lte = extreme.full_seastate_long_term_extreme(ste_all, sample_weights)

Similar to the short-term extreme functions, the output of long-term extreme function is a probability distribution (scipy.stats.rv_continuous). This object provides common statistical functions (PDF, CDF, PPF, etc.) and metrics (expected value, median, etc). Here, we will look at the survival function and the 100-year return level.

The value of the survival function at a given return level (e.g. 100-years) may be calculated using the return_year_value function. This gives a 100-year wave of about 11.1 meters.

t_st_hr = t_st / (60.0 * 60.0)
t_return_yr = 100.0
x_t = extreme.return_year_value(lte.ppf, t_return_yr, t_st_hr)

print(f"100-year elevation: {x_t} m")

100-year elevation: 10.949955788933982 m

Finally we plot the survival function and show the 100-year return level (dashed grey lines). The 100-year value is about 11.1 m (where the grey line intersects the x-axis)

x = np.linspace(0, 20, 1000)

fig, ax = plt.subplots()

# plot survival function
ax.semilogy(x, lte.sf(x))

# format plot
plt.grid(True, which="major", linestyle=":")
ax.tick_params(axis="both", which="major", direction="in")
ax.xaxis.set_ticks_position("both")
ax.yaxis.set_ticks_position("both")
plt.minorticks_off()
ax.set_xticks([0, 5, 10, 15, 20])
ax.set_yticks(1.0 * 10.0 ** (-1 * np.arange(11)))
ax.set_xlabel("elevation [m]")
ax.set_ylabel("survival function (1-cdf)")
ax.set_xlim([0, x[-1]])
ylim = [1e-10, 1]
ax.set_ylim(ylim)

# 100-year return level
s_t = lte.sf(x_t)
ax.plot([0, x[-1]], [s_t, s_t], "--", color="0.5", linewidth=1)
ax.plot([x_t, x_t], ylim, "--", color="0.5", linewidth=1)

[<matplotlib.lines.Line2D at 0x27ff4207fd0>]