How can I get each turbine AEP in a farm ?

[mufannao]

How can I obtain the Annual Energy Production (AEP) for each individual turbine in a wind farm, rather than the total AEP of the entire wind farm?

Thank you for all attention.

[Ooussamaa]

Hello @mufannao,
with the command "fi.get_turbine_powers()"

Try this script it should give you the average roduction of each turbine.

powers = fi.get_turbine_powers() / 1000.0
N_TURBINES = fi.floris.farm.n_turbines
print()
print("Turbine powers for the first given wind speed ()")
for i in range(3):
    print(f"Wind direction {i+1}")
    for j in range(N_TURBINES):
        print(f"  Turbine {j} - {powers[i, 0, j]:7,.2f} kW")
    print()

Best regards

[mufannao]

@Ooussamaa ,Thank you for your reply. Perhaps I did not explain my problem clearly, I originally wanted to calculate the annual power generation of a single wind turbine in a wind farm through the wind rose diagram, and now I should have solved this problem. Here is the complete code.

 Copyright 2022 NREL

# Licensed under the Apache License, Version 2.0 (the "License"); you may not
# use this file except in compliance with the License. You may obtain a copy of
# the License at http://www.apache.org/licenses/LICENSE-2.0

# Unless required by applicable law or agreed to in writing, software
# distributed under the License is distributed on an "AS IS" BASIS, WITHOUT
# WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. See the
# License for the specific language governing permissions and limitations under
# the License.

# See https://floris.readthedocs.io for documentation


import numpy as np
import pandas as pd
from scipy.interpolate import NearestNDInterpolator
from floris.tools import FlorisInterface

"""
This example demonstrates how to calculate the Annual Energy Production (AEP)
of a wind farm using wind rose information stored in a .csv file.

The wind rose information is first loaded, after which we initialize our Floris
Interface. A 3 turbine farm is generated, and then the turbine wakes and powers
are calculated across all the wind directions. Finally, the farm power is
converted to AEP and reported out.
"""

# Read the windrose information file and display
df_wr = pd.read_csv("inputs/wind_rose.csv")
print("The wind rose dataframe looks as follows: \n\n {} \n".format(df_wr))

# Derive the wind directions and speeds we need to evaluate in FLORIS
wd_array = np.array(df_wr["wd"].unique(), dtype=float)
ws_array = np.array(df_wr["ws"].unique(), dtype=float)

# Format the frequency array into the conventional FLORIS v3 format, which is
# an np.array with shape (n_wind_directions, n_wind_speeds). To avoid having
# to manually derive how the variables are sorted and how to reshape the
# one-dimensional frequency array, we use a nearest neighbor interpolant. This
# ensures the frequency values are mapped appropriately to the new 2D array.
wd_grid, ws_grid = np.meshgrid(wd_array, ws_array, indexing="ij")
freq_interp = NearestNDInterpolator(df_wr[["wd", "ws"]], df_wr["freq_val"])
freq = freq_interp(wd_grid, ws_grid)

# Normalize the frequency array to sum to exactly 1.0
freq = freq / np.sum(freq)
######################################################################################
# Load the FLORIS object
fi = FlorisInterface("inputs/gch.yaml") # GCH model


# Assume a three-turbine wind farm with 5D spacing. We reinitialize the
# floris object and assign the layout, wind speed and wind direction arrays.
D = fi.floris.farm.rotor_diameters[0] # Rotor diameter for the NREL 5 MW
fi.reinitialize(
    layout_x=[0.0, 5 * D, 10 * D],
    layout_y=[0.0, 0.0, 0.0],
    wind_directions=wd_array,
    wind_speeds=ws_array,
)
fi.calculate_wake()
powers = fi.get_turbine_powers()
power_use = np.nan_to_num(powers[:,:,:])

N_TURBINES = fi.floris.farm.n_turbines
all = []
for i in range(N_TURBINES):
    farm_power = power_use[:,:,i]
    aep = np.sum(np.multiply(freq, farm_power) * 365 * 24)
    print(aep/ 1.0e9)#   The annual electricity generation is converted to GWH
    all.append(aep/ 1.0e9)#A list of the annual power generation of each wind turbine in the wind farm

Best regards