com.farao_community.farao.cse.computation.BorderExchanges Maven / Gradle / Ivy

Go to download
Show more of this group Show more artifacts with this name
Show all versions of gridcapa-cse-computation Show documentation
Module that manages common computation for CSE processes
The newest version!
/*
 * Copyright (c) 2022, RTE (http://www.rte-france.com)
 * This Source Code Form is subject to the terms of the Mozilla Public
 * License, v. 2.0. If a copy of the MPL was not distributed with this
 * file, You can obtain one at http://mozilla.org/MPL/2.0/.
 */

package com.farao_community.farao.cse.computation;

import com.powsybl.balances_adjustment.util.CountryArea;
import com.powsybl.balances_adjustment.util.CountryAreaFactory;
import com.powsybl.iidm.network.Country;
import com.powsybl.iidm.network.Network;
import com.powsybl.loadflow.LoadFlow;
import com.powsybl.loadflow.LoadFlowParameters;
import com.powsybl.loadflow.LoadFlowResult;
import org.slf4j.Logger;
import org.slf4j.LoggerFactory;

import java.util.HashMap;
import java.util.Map;
import java.util.function.Function;
import java.util.stream.Collectors;
import java.util.stream.Stream;

/**
 * @author Joris Mancini {@literal }
 */
public final class BorderExchanges {
    private static final Logger LOGGER = LoggerFactory.getLogger(BorderExchanges.class);
    public static final String IT_CH = "IT-CH";
    public static final String IT_FR = "IT-FR";
    public static final String IT_AT = "IT-AT";
    public static final String IT_SI = "IT-SI";
    public static final String CH_FR = "CH-FR";
    public static final String FR_DE = "FR-DE";
    public static final String CH_DE = "CH-DE";

    private BorderExchanges() {
        // Should not be instantiated
    }

    public static double computeItalianImport(Network network) {
        runLoadFlow(network);
        CountryArea itArea = new CountryAreaFactory(Country.IT).create(network);
        return Stream.of(Country.FR, Country.AT, Country.CH, Country.SI)
            .map(country -> new CountryAreaFactory(country).create(network).getLeavingFlowToCountry(itArea))
            .reduce(0., Double::sum);
    }

    public static Map computeCseBordersExchanges(Network network) {
        return computeCseBordersExchanges(network, true);
    }

    public static Map computeCseBordersExchanges(Network network, boolean withLoadflow) {
        if (withLoadflow) {
            runLoadFlow(network);
        }
        Map borderExchanges = new HashMap<>();
        Map countryAreaPerCountry = Stream.of(Country.FR, Country.AT, Country.CH, Country.SI, Country.IT, Country.DE)
            .collect(Collectors.toMap(
                Function.identity(),
                country -> new CountryAreaFactory(country).create(network)
            ));
        borderExchanges.put(IT_CH, getBorderExchange(Country.IT, Country.CH, countryAreaPerCountry));
        borderExchanges.put(IT_FR, getBorderExchange(Country.IT, Country.FR, countryAreaPerCountry));
        borderExchanges.put(IT_AT, getBorderExchange(Country.IT, Country.AT, countryAreaPerCountry));
        borderExchanges.put(IT_SI, getBorderExchange(Country.IT, Country.SI, countryAreaPerCountry));
        borderExchanges.put(CH_FR, getBorderExchange(Country.CH, Country.FR, countryAreaPerCountry));
        borderExchanges.put(FR_DE, getBorderExchange(Country.FR, Country.DE, countryAreaPerCountry));
        borderExchanges.put(CH_DE, getBorderExchange(Country.CH, Country.DE, countryAreaPerCountry));
        return borderExchanges;
    }

    private static double getBorderExchange(Country fromCountry, Country toCountry, Map countryAreaPerCountry) {
        return countryAreaPerCountry.get(fromCountry).getLeavingFlowToCountry(countryAreaPerCountry.get(toCountry));
    }

    public static Map computeCseCountriesBalances(Network network) {
        runLoadFlow(network);
        Map countriesBalances = new HashMap<>();
        countriesBalances.put("AT", new CountryAreaFactory(Country.AT).create(network).getNetPosition());
        countriesBalances.put("CH", new CountryAreaFactory(Country.CH).create(network).getNetPosition());
        countriesBalances.put("FR", new CountryAreaFactory(Country.FR).create(network).getNetPosition());
        countriesBalances.put("SI", new CountryAreaFactory(Country.SI).create(network).getNetPosition());
        countriesBalances.put("IT", new CountryAreaFactory(Country.IT).create(network).getNetPosition());
        countriesBalances.put("DE", new CountryAreaFactory(Country.DE).create(network).getNetPosition());
        return countriesBalances;
    }

    private static void runLoadFlow(Network network) {
        LoadFlowResult result = LoadFlow.run(network, LoadFlowParameters.load());
        if (!result.isOk()) {
            LOGGER.error("Loadflow computation diverged on network '{}'", network.getId());
            throw new CseComputationException(String.format("Loadflow computation diverged on network %s", network.getId()));
        }
    }
}