ResInsight/ApplicationLibCode/ProjectDataModel/RimSimWellInViewTools.cpp

/////////////////////////////////////////////////////////////////////////////////
//
//  Copyright (C) 2019-     Equinor ASA
//
//  ResInsight is free software: you can redistribute it and/or modify
//  it under the terms of the GNU General Public License as published by
//  the Free Software Foundation, either version 3 of the License, or
//  (at your option) any later version.
//
//  ResInsight is distributed in the hope that it will be useful, but WITHOUT ANY
//  WARRANTY; without even the implied warranty of MERCHANTABILITY or
//  FITNESS FOR A PARTICULAR PURPOSE.
//
//  See the GNU General Public License at <http://www.gnu.org/licenses/gpl.html>
//  for more details.
//
/////////////////////////////////////////////////////////////////////////////////

#include "RimSimWellInViewTools.h"

#include "RiaLogging.h"
#include "RiaSummaryTools.h"

#include "RifEclipseSummaryAddress.h"
#include "RifSummaryReaderInterface.h"

#include "RigSimWellData.h"
#include "RigWellResultPoint.h"

#include "Rim3dView.h"
#include "RimEclipseResultCase.h"
#include "RimGridSummaryCase.h"
#include "RimOilField.h"
#include "RimProject.h"
#include "RimSimWellInView.h"
#include "RimSimWellInViewCollection.h"
#include "RimSummaryCaseMainCollection.h"

#include "cvfAssert.h"

//--------------------------------------------------------------------------------------------------
///
//--------------------------------------------------------------------------------------------------
RimGridSummaryCase* RimSimWellInViewTools::gridSummaryCaseForWell( RimSimWellInView* well )
{
    RimProject* project = RimProject::current();
    if ( !project ) return nullptr;

    RimSummaryCaseMainCollection* sumCaseColl =
        project->activeOilField() ? project->activeOilField()->summaryCaseMainCollection() : nullptr;
    if ( !sumCaseColl ) return nullptr;

    RimEclipseResultCase* eclCase = nullptr;
    well->firstAncestorOrThisOfType( eclCase );
    if ( eclCase )
    {
        RimGridSummaryCase* gridSummaryCase =
            dynamic_cast<RimGridSummaryCase*>( sumCaseColl->findSummaryCaseFromEclipseResultCase( eclCase ) );
        if ( gridSummaryCase )
        {
            return gridSummaryCase;
        }
    }

    return nullptr;
}

//--------------------------------------------------------------------------------------------------
///
//--------------------------------------------------------------------------------------------------
std::vector<RimSummaryCase*> RimSimWellInViewTools::summaryCases()
{
    std::vector<RimSummaryCase*> cases;

    RimProject* project = RimProject::current();
    if ( project )
    {
        RimSummaryCaseMainCollection* sumCaseColl =
            project->activeOilField() ? project->activeOilField()->summaryCaseMainCollection() : nullptr;
        if ( sumCaseColl )
        {
            cases = sumCaseColl->allSummaryCases();
        }
    }

    return cases;
}

//--------------------------------------------------------------------------------------------------
///
//--------------------------------------------------------------------------------------------------
bool RimSimWellInViewTools::isInjector( RimSimWellInView* well )
{
    RigSimWellData* wRes = well->simWellData();
    if ( wRes )
    {
        Rim3dView* rimView = nullptr;
        well->firstAncestorOrThisOfTypeAsserted( rimView );

        int currentTimeStep = rimView->currentTimeStep();

        if ( wRes->hasWellResult( currentTimeStep ) )
        {
            const RigWellResultFrame* wrf = wRes->wellResultFrame( currentTimeStep );

            if ( RiaDefines::isInjector( wrf->m_productionType ) )
            {
                return true;
            }
        }
    }

    return false;
}

//--------------------------------------------------------------------------------------------------
///
//--------------------------------------------------------------------------------------------------
double RimSimWellInViewTools::extractValueForTimeStep( RifSummaryReaderInterface* summaryReader,
                                                       const QString&             wellName,
                                                       const std::string&         vectorName,
                                                       const QDateTime&           currentDate,
                                                       bool*                      isOk )

{
    CVF_ASSERT( summaryReader );

    if ( vectorName.empty() )
    {
        *isOk = true;
        return 0.0;
    }

    RifEclipseSummaryAddress addr( RifEclipseSummaryAddress::SUMMARY_WELL,
                                   vectorName,
                                   -1,
                                   -1,
                                   "",
                                   wellName.toStdString(),
                                   -1,
                                   "",
                                   -1,
                                   -1,
                                   -1,
                                   -1,
                                   false,
                                   -1 );

    if ( !summaryReader->hasAddress( addr ) )
    {
        QString message = "ERROR: no address found for well " + wellName + " " + QString::fromStdString( vectorName );
        RiaLogging::warning( message );

        *isOk = false;
        return 0.0;
    }

    std::vector<double> values;
    summaryReader->values( addr, &values );
    std::vector<time_t> timeSteps = summaryReader->timeSteps( addr );
    if ( values.empty() || timeSteps.empty() )
    {
        QString message = "ERROR: no data found for well " + wellName + " " + QString::fromStdString( vectorName );
        RiaLogging::warning( message );

        *isOk = false;
        return 0.0;
    }

    auto [resampledTimeSteps, resampledValues] =
        RiaSummaryTools::resampledValuesForPeriod( addr, timeSteps, values, RiaQDateTimeTools::DateTimePeriod::DAY );

    time_t currentTime_t = currentDate.toSecsSinceEpoch();

    for ( unsigned int i = 0; i < resampledTimeSteps.size(); i++ )
    {
        if ( resampledTimeSteps[i] > currentTime_t )
        {
            *isOk = true;
            return resampledValues[i];
        }
    }

    QString message = "ERROR: no resampled values found for well " + wellName + " " + QString::fromStdString( vectorName );
    RiaLogging::warning( message );

    *isOk = false;
    return -1;
}