Regular Article
Estimating the chance of success of information acquisition for the Norne benchmark case
^{1}
University of Campinas – Center for Petroleum Studies, PO Box 6052, 13083970
Campinas, SP, Brazil
^{2}
SOLPE – Solutions for Petroleum Science and Engineering, Av. Alan Turing, 345, Room 12, 13083898
Campinas, SP, Brazil
^{*} Corresponding author: botechia@cepetro.unicamp.br
Received:
26
April
2018
Accepted:
27
August
2018
A key decision in field management is whether or not to acquire information to either improve project economics or reduce uncertainties. A widely spread technique to quantify the gain of information acquisition is Value of Information (VoI). However, estimating the possible outcomes of future information without the data is a complex task. While traditional VoI estimates are based on a single average value, the Chance of Success (CoS) methodology works as a diagnostic tool, estimating a range of possible outcomes that vary because of reservoir uncertainties. The objective of this work is to estimate the CoS of a 4D seismic before having the data, applied to a complex real case (Norne field). The objective is to assist the decision of whether, or not, to acquire further data. The methodology comprises the following steps: uncertainty quantification, selection of Representative Models (RMs), estimation of the acquisition period, production strategy optimization and, finally, quantification of the CoS. The estimates use numerical reservoir simulation, economic analysis, and uncertainty evaluation. We performed analyses considering perfect and imperfect information. We aim to verify the increment in economic return when the 4D data identifies the closesttoreality reservoir model. While the traditional expected VoI calculation provides only an average value, this methodology has the advantage of considering the increase in the economic return due to reservoir uncertainties, characterized by different RMs. Our results showed that decreased reliability of information affected the decision of which production strategy to select. In our case, information reliability less than 70% is insufficient to change the perception of the uncertain reservoir and consequently decisions. Furthermore, when the reliability reached around 50%, the information lost value, as the economic return became similar to that of the case without information acquisition.
© V. Botechia et al., published by IFP Energies nouvelles, 2018
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Nomenclature
EVII: Expected Value of Imperfect Information
EVoI: Expected Value of Information
EVPI: Expected Value of Perfect Information
Max: Maximum value of determined indicator
Min: Minimum value of determined indicator
NPVwi: Net Present Value with information
NPVwoi: Net Present Value without information
N_{SCNr} : Number of scenarios having the lowest distance value (d) with respect to the RM_{r}
P_{r} : Reliability probability
P_{RMr} : Probability of the representative model r
Sw_{i}^{base} : Water saturation of the base model for each grid block i
Sw_{i}^{scenario} : Water saturation of each scenario for each grid block i
T_{SCN} : Total number of scenarios
VOPI: Value Of Perfect Information
ΔNPV: Difference between NPVwi and NPVwoi
ε_{norm} : Normalized water saturation error
1 Introduction
An efficient way to evaluate strategies is fundamental to maximize economic performance of field operation and infrastructure decisions. However, this can be very complex since many decisions are made in scenarios involving high investments and multiple uncertainties. The decision to acquire information can aid the selection of a better strategy, but the costs and benefits must be justified.
A widely used tool to assist decision makers in how to reduce uncertainty before selecting a course of action is Value of Information (VoI) (Schlaifer, 1959; Grayson, 1960; Howard, 1966). A classical definition of VoI is the difference between the expected value with information and the expected value without information (Coopersmith and Cunningham, 2002). However, Bratvold et al. (2009) state that this is a misconception of VoI, while the concept is true in many cases, for example, under risk neutrality, it is not in general.
The value of gathering new information depends on several factors: degree of uncertainty (the greater the uncertainties, the greater the VoI), representativeness of the information, decision flexibility, the magnitude of future monetary impact on changing decisions, and the reliability of the information (Gerhardt and Haldorsen, 1989; Coopersmith and Cunningham, 2002).
Eidsvik et al. (2015) presented applications of VoI in earth sciences, integrating this topic with spatial modeling. Abellan and Noetinger (2010) proposed an approach that allows quantifying the information contained in measurements in a geoscience context. Recently, the VoI concept is being used in history matching contexts, involving data gathered at different times (Barros et al., 2015, 2016; He et al., 2016; Chen et al., 2017; Hong et al., 2018). Hong et al. (2018) discussed the variability of VoI over ensembles, including Bayesian inference with ensemble Kalman filter (EnKF) (Evensen, 1994).
An important source of information that can add value to an oil field project is 4D seismic data because it provides a better understanding of the variations of dynamic properties in the reservoir (Davolio and Schiozer, 2018). This type of information has been widely used in modelbased approaches, such as reservoir management, production optimization, and history matching techniques (Le Ravalec et al., 2012a, 2012b; Roggero et al., 2012; Tillier et al., 2012).
The presence of uncertainties and the complex interrelationships involved in large development decisions require probabilistic approaches to calculate VoI (Begg et al., 2002). Decision trees and Bayesian analyses (statistical method to revise probability estimates given new information – Waggoner, 2002) are widely used in the literature to assess VoI (Raiffa, 1968; Waggoner, 2002; Galli et al., 2004; Ballin et al., 2005; Bailey et al., 2011; Santos et al., 2017).
If the information provides perfect knowledge of the state of the world, then it is called perfect information, and VoI is referred to as Value Of Perfect Information (VOPI) or Expected Value of Perfect Information (EVPI) (Bratvold et al., 2009). However, in practice, information is rarely perfect because of uncertainty. This is imperfect information. In fact, perfect information works as an upper limit for imperfect information (Coopersmith and Cunningham, 2002), and the VoI, in this case, can be also referred to as Expected Value of Imperfect Information (EVII).
Usually, VoI determines the Expected Monetary Value (EMV) of acquiring and analyzing data to reduce uncertainty (Begg et al., 2002). Estimating VoI before data acquisition is complex because of the many possible outcomes and must consider the uncertainties of the reservoir models (Dunn, 1992). In this context, Ferreira and Schiozer (2014) use the term Expected Value of Information (EVoI) when the valuation is made before the data acquisition, and this is the term we use in the text from now on.
Ferreira et al. (2015) claimed that being a weighted average, the EVoI does not reflect the variation in the increase of the expected revenue due to information acquisition. In order to complement the traditional approach, the authors proposed the Chance of Success (CoS) methodology, a diagnostic tool that gives ranges of possible outcomes for the increase in the difference of the Net Present Value (NPV) with (NPVwi) and without information (NPVwoi) considering different scenarios. This allows the estimation of the probability of increased expected revenue to surpass the costs of information acquisition, considering perfect information and a simple synthetic reservoir model.
The idea of this approach is to decide whether the information should be acquired. If the information is acquired, it will inform which the mostlikely RM is; hence, we can decide which production strategy should be selected.
2 Objectives and assumptions
This work estimates the CoS of a 4D seismic project before acquiring 4D seismic data (based on the works of Ferreira and Schiozer, 2014 and Ferreira et al., 2015), applied to the Norne benchmark model, considering perfect and imperfect information. The objective is to identify which Representative Model (RM) is closest to the real model and assess the increment on the economic return.
3 Methodology
The proposed CoS methodology is based on the works of Ferreira and Schiozer (2014) and Ferreira et al. (2015), comprising the six steps shown in Figure 1.
Fig. 1. Methodology to estimate the chance of success of 4D seismic acquisition (adapted from Ferreira and Schiozer, 2014). 
3.1 Uncertainty and risk analyses
The first step defines the uncertain parameters of the reservoir model. Statistical procedures combine all uncertainties to generate many scenarios (simulation models). A base case is chosen from these models (usually that closest to the percentile P50).
We optimize a production strategy for the base case (base production strategy). This strategy is applied to all the other models, the NPV for each model is calculated, and a risk curve (descendent or complementary cumulative distribution function) is plotted.
3.2 Selection of representative models
Optimizing the hundreds of models generated in the first step is computationally impractical, so in this step, RMs are selected to reduce the number of models to be optimized (Schiozer et al., 2004; Costa et al., 2008; Marques et al., 2013; Sarma et al., 2013; Shirangi and Durlofsky, 2015; Meira et al., 2016). The RMs represent the variability of the uncertainties (input and output) in a few models. These are then optimized to have different project options, each representing the best decision for different scenarios. We base the selection of representative models on the procedure by Meira et al. (2017), which evaluates the variability of the models on the following parameters: NPV, oil Recovery Factor (RF), cumulative oil (Np), and cumulative water (Wp).
Each RM may have a specific probability of occurrence, which we determine according to Ferreira et al. (2015) – equation (1):(1)where P_{RMr} is the probability of the representative model r, N_{SCNr} is the number of scenarios with the lowest distance value (d) with respect to the representative model RMr and T_{SCN} is the total number of scenarios.
The measure referred to as “distance” (d) is a comparison between the normalized production and economic results of each representative model and is calculated using equation (2):(2)where d is the distance of the normalized production and economic results from the scenario j to the representative model (RMr) under analysis. Max and Min refer to the maximum and minimum values, respectively, obtained for each indicator among all scenarios.
3.3 Estimation of the acquisition period
The time for 4D seismic acquisition is best when: (1) only 4D seismic data can identify that the base model is nonrepresentative of the true earth model, and (2) there is time to make or change decisions (Ferreira et al., 2011).
Thus, the value of 4D seismic data is affected by the period of acquisition. We evaluate the information from the 4D seismic data using water saturation maps to consider two points: (1) ensure that the base reservoir model does not represent the true earth model, and (2) anticipate the water breakthrough for producers.
To find the best time to acquire information, we use the methods by Ferreira et al. (2011) and Ferreira and Schiozer (2013), which calculate an error function that represents the mismatch between the simulation scenario and the water saturation maps for the base model (Eq. (3)). The normalized error function is given by equation (4).(3) (4)where is the error, Ng is the number of grid cells, Sw_{i}^{scenario} is the water saturation of a specific scenario for each grid block i and Sw_{i}^{base} is the water saturation of the base model for each grid block.
3.4 Production strategy optimization
In this step, the production strategy for each RM is optimized for the time of data acquisition and processing. This enables the quantification of the economic impact of the 4D seismic data acquisition since these data identify the true representative model, and so the respective strategy should be implemented. The optimization aims to maximize the NPV.
This work uses a field model with a predeveloped strategy and 4 years of production history (i.e. the field is in the management phase), we consider the following variables in the optimization process: production and injection rates, optimal time to shutin producers and injectors, implementation of sidetracks in existing wells and allocation of new wells.
3.5 Chance of Success (CoS) analysis
We estimate the EVoI using the results obtained in the previous step. The probability of the EVoI must be high enough to cover the cost of data acquisition.
The first analysis considers perfect information. We compare the traditional EVoI calculation with the CoS analysis. The EVoI gives an average value for the information value and is calculated as the difference between the EMV of the project with information and the EMV without information (Eq. (5)). The EMV is the sum of the NPV of each scenario, weighted by its respective probability of occurrence.(5)
For the CoS analysis, we consider the cumulative probability curve of the ΔNPV for each RM. The ΔNPV is the difference between the NPV with information (NPV_{wi}) and the NPV without information (NPV_{woi}) (Eq. (6)) for each RM. NPV_{woi} refers to the economic return of the production strategy of each RM chosen if 4D seismic is not acquired, while NPV_{wi} refers to the economic return of the production strategies optimized for each RM.(6)
Tables 1 and 2 exemplify the ΔNPV calculation, as well as the difference between EVoI and ΔNPV. The number of the strategy in Table 1 refers to the number of the RM it was optimized for. For instance, S1 is the strategy optimized for RM1. This strategy is simulated for all RMs, so that the EMV of each strategy can be calculated. To simplify this example, we consider the RMs having the same probability of occurrence prior to information acquisition.
Example of NPV for three RMs (RM1 to RM3, equiprobable prior to information acquisition) simulated for all production strategies (S1–S3) and respective EMV. S3 is selected as the best strategy without further information.
Example of NPV without information (NPV_{woi}), NPV with information (NPV_{wi}) and ΔNPV for the three RMs.
Note that, in this example, S3 presents the highest EMV (EMV_{without information} = USD 3092 million). Thus, this strategy would be the chosen if we do not acquire information, and the ΔNPV is calculated as the difference between the NPV of the optimized strategy (with information) and the NPV of strategy S3 applied to each RM (without information).
The EMV with information is the average NPV of the optimized strategies (EMV_{with information} =] (3022 + 3172 + 3204)/3 = USD 3133 Million). The EVoI, in this example is (3133 − 3092) = USD 41 million.
Figure 2 shows the CoS analysis for this example. The economic chance of success depends on the cost of the information acquisition. Suppose that a seismic cost is [USD 30 million, thus the chance of success is around 60% (chance of the outputs be higher than this value).
Fig. 2. Example of CoS analysis – ΔNPV × cumulative probability. 
Then, we verify the EVoI when its reliability decreases (imperfect information). To deal with imperfect information, we update the probabilities of occurrence of each RM based on the Bayesian theorem (Clemen, 1996; Santos et al., 2017) shown in equations (7) and (8). This updating procedure is dependent on the level of reliability. A reliability of 100% reflects perfect information, and as the reliability decreases, the EVoI information decreases. (7) (8)where P_{u} is updated probability, P_{r} is reliability probability, P_{c} is the current probability and P_{P} is predicted probability, i is the index or number of the RM, and n is the total number of RM.
That is, P_{P} is the probability of the seismic data to predict which RM is the closest to the true reservoir model. The purpose of P_{u} is to answer the following question: if the seismic data indicate a specific RM to be the closest to the true reservoir model, what is the probability of this information to be correct (i.e. is the indicated RM truly the closest to the real reservoir model)? Finally, P_{Ci} is the initial probability, calculated as shown in the previous subsections.
3.6 Decisionmaker evaluation
Based on the results from previous analyses, the decision maker can end the process (already deciding whether the information must be acquired) or reevaluate by restarting the process with the following possibilities: (1) improving the accuracy of analyses by choosing more RMs, (2) evaluating a different acquisition period, or (3) improving the production strategy optimization process.
4 Application
This study is applied to the Norne benchmark case. The Norne field is located in the Norwegian continental shelf (200 km west of the midNorway coast), extending for 10 km, 2 km wide, and at a water depth of 380 m (Adlam, 1995). The Norne field comprises two separate oil compartments, the Norne Main structure (Norne C, D, and E segments) and the Northeast Segment (Norne GSegment) (Fig. 3a). It also consists of five formations, from top to base: Garn, Not, Ile, Tofte, and Tilje (Lind, 2004).
Fig. 3. (a) Norne segments (Lind, 2004) and (b) 3D view of porosity map for the simulation model. 
The simulation model is a cornerpoint grid with 46 × 112 × 22 blocks (44 431 active cells). Figure 3b shows a 3D view of a porosity map of the simulation model. This case has four years of production data (from 1997 until 2001), with a production strategy comprising 12 producers and 8 injectors.
4.1 Economic parameters
Table 3 shows the economic parameters for the optimization process. Since this is the Norne 2001 case, we considered the optimization for that time, and hence these were the values in force in 2001.
Economic parameters
5 Results
5.1 Risk analysis and selection of representative models
The uncertainties considered in this work were: rock compressibility, oilwater contact, gasoil contact, horizontal and vertical permeability multipliers, porosity multipliers, relative permeability tables, and fault transmissibility. Combining all these uncertainties, 2000 models were generated and a history matching process was performed, filtering 89 models which honored production data. Figure 4 shows the cumulative oil and water production for all scenarios.
Fig. 4. Cumulative (a) oil and (b) water production for all scenarios. 
We selected nine representative models based on risk curves and crossplots for the following indicators: NPV, Np, Wp, and RF. Figure 5 shows the risk curves for NPV and Np, while Figure 6 exemplifies crossplots analyses. In these graphs, each point represents a simulation model; the purple dots represent all scenarios, while the black dots are the selected RMs. Note that the RMs represent the variability of the uncertain outputs.
Fig. 5. Risk curves for (a) NPV and (b) Np. Purple dots represent all scenarios; black dots, the RMs. 
Fig. 6. Crossplots: (a) NPV × Np, (b) NPV × Wp, (c) NPV × RF, and (d) Np × Wp. Purple dots represent all scenarios; black dots, the RMs. 
Table 4 shows the NPV, cumulative oil and water production, and the probability of occurrence for each RM (calculated using Eqs. (1) and (2)).
NPV, cumulative oil and water production, and the probability of occurrence of each RM.
5.2 Acquisition period estimation
Figure 7 shows the curves of water saturation errors, the sum of the difference of water saturation over every grid cell. There are many scenarios with fluid distributions very different from the base model (S_{w} error higher than 50%) from the beginning of the analysis (year 4). Note that the first 4 years relate to past production of the production strategy already implemented, thus our analysis starts at this point.
Fig. 7. (a) Water saturation error curves and (b) normalized water saturation error curves. 
Figure 8 shows the histograms of production time for a normalized water saturation error of 60% (Fig. 8a) and 50% (Fig. 8b), confirming the high level of differences in water distribution between models.
Fig. 8. Histograms – number of scenarios that reached an Sw error of (a) 0.60 and (b) 0.50 in each year of production. 
To estimate the best time to acquire 4D seismic data, we also consider the water breakthrough. If the seismic data are acquired before breakthrough, the water flow path can be identified and so the well control optimized to improve project economics. Figure 9 shows the number of wells (considering all scenarios) which reached the breakthrough for each production time (Fig. 9a) and the percentage of the number of wells which reached the breakthrough for each production time (Fig. 9b). Note that nearly 50% of the wells reached breakthrough in the first year of the analysis (year 4).
Fig. 9. (a) Histogram of number of wells (considering all scenarios) according to time of breakthrough (years), and (b) rate of total number of wells which reached breakthrough in each year of production. 
Considering the water saturation error curves and the breakthrough time, we concluded that the best acquisition period is at the beginning of the management phase, i.e., 4 years after the start of production.
5.3 Production strategy optimization
After estimating the best time for 4D seismic data acquisition, we then optimize the production strategy. As this case is in the management phase (i.e., the production strategy is already implemented), there is little flexibility for changes. The variables we considered are: production and injection rates, optimal shutin time of producers and injectors, implementation of sidetracks in existing wells, and allocation of new wells.
The NPV_{woi} is the economic return for the selected strategy without information acquisition (i.e. the strategy that maximizes EMV) applied to each RM. The economic return after each RM is optimized is NPV_{wi}. Figure 10 shows the difference between NPV_{wi} and NPV_{woi}. Table 5 shows the results for ΔNPV calculation. As S3 is the strategy that maximizes EMV (Tab. 6) and, hence, it would be chosen if no information is acquired, the ΔNPV of this strategy is 0.
Fig. 10. ΔNPV for each strategy. 
Net Present Value without information (NPV_{woi}), Net Present Value with information (NPV_{wi}), ΔNPV and probability of occurrence (Prob.) for each strategy.
NPV of each RM (RM1 to RM9) simulated for all production strategies (S1–S9) and respective EMV. The bold values represent the NPV_{wi} of each strategy. The bold value in EMV line represents the maximum EMV among the strategies.
5.4 Chance of Success (CoS) analysis
5.4.1 Perfect information analysis
First, we considered perfect information and calculated the EVoI using the traditional decision tree, for comparison with the chance of success methodology.
Figure 11 shows the decision tree used to calculate the EVoI. The green squares are decision nodes, the brown circles are chance nodes. The branch labeled “Do not acquire 4DS” represents the possible outcomes from the decision not to acquire information. We simulated the set of production strategies (previously optimized deterministically for each RM) in all possible RMs. We named the strategies according to the RM they were optimized for. Thus, S1 is the strategy optimized for RM1, S2 is the strategy optimized for RM2, and so on. Then, we calculated the EMV of each strategy. If seismic data is not acquired, we chose the production strategy S3, which maximized EMV (Tab. 6).
Fig. 11. Compact decision tree for the calculation of the EVoI. 
The chance node labeled “Acquire 4DS” is in a compact form, since we treated the information as perfect. Thus, as we expect seismic data to provide enough information to identify the closesttoreality RM, the EMV of the 4DS project is the NPV_{wi} of each RM multiplied by its probability of occurrence.
The EVoI – the difference between the EMV when seismic data are acquired and are not acquired – reached the value of USD 47 million.
Being based on a single value, the EVoI cannot express the variability of the increase in the NPV caused by reservoir uncertainties. In the estimation of the chance of success of the 4DS project, we consider the cumulative probability curve of ΔNPV for each RM (Fig. 12). Moreover, the chance of success also depends on the cost of information acquisition, i.e. the increased economic return (ΔNPV) must be higher than the cost to acquire information. In this study, for a chance of success higher than 50%, the seismic data must cost less than USD 32 Million.
Fig. 12. CoS evaluation for perfect information: cumulative probability curve of the ΔNPV. 
Note the difference between the response of EVoI and CoS: while the former estimated an average economic gain of USD 47 Million, the latter estimated a 45% chance of the 4DS project yielding a higher return than that (Fig. 12). As these analyses assume perfect information, they provide an upper limit for value information estimate.
Ferreira et al. (2015) showed that the number of selected RMs impacts the analysis. Higher number of RMs leads to more accurate results, but at a certain point the results stabilize. However, choosing many RMs would be too much time consuming, hence it is necessary to balance the accuracy of the results and the time consumption.
In this sense, we performed a sensitivity analysis on the number of RMs (Fig. 13). Using only 3 RMs, the results are very different from the previously obtained. However, with 6 RMs, the results are very similar to the ones obtained with 9 RMs.
Fig. 13. Sensitivity analysis – CoS with different number of RMs. 
5.4.2 Imperfect information analysis
As information is usually imperfect, some unreliability is expected. In this case, there is a probability that incorrect information could be provided, diminishing the economic return.
Table 7 shows an example of P _{r} (reliability probability) for an information reliability of 0.9. That is, if the true model is RM1, there is a probability of 90% of the seismic data predicting RM1 to be the true model. The other 10% is divided between the other RMs, weighted by their current probability of occurrence. The same applies to the other RMs.
Example of P _{r} (reliability probability) for an information reliability of 0.9.
Table 8 shows the values of P _{p} (predicting probabilities), and Table 9 shows an example of P _{u} (updated probabilities) values, both for an information reliability of 0.9. In this example, if the seismic data predicts RM1 to be the true model, then there is a probability of 89.4% of RM1 being the true model.
Example of P _{p} values for an information reliability of 0.9.
Example of updated probabilities (P _{u}) for an information reliability of 0.9.
Figure 14 exemplifies a decision tree to calculate EVoI considering a reliability of 0.9. The example presents only the branch “Acquire 4DS” since the branch to calculate the EMV without information is the same as that shown in Figure 11. Also, we exemplify with one strategy for the first RM, but the branches should be expanded to all nine strategies for all RMs.
Fig. 14. Simplified decision tree for EVoI calculation considering imperfect information with a reliability of 0.9. 
We performed a sensitivity analysis on the EMV of the 4DS project according to the reliability of the information (Fig. 15). When information reliability nears 50%, it is better not to acquire the information. Similarly, we calculated the EVII, according to each level of reliability, which also showed to be closer to zero as the reliability tends to 50% (Fig. 16).
Fig. 15. Variation of the EMV of the 4DS project with information reliability. 
Fig. 16. Variation of EVII according to each level of reliability. 
The level of reliability also affected the strategy to be chosen when the seismic data indicates an RM to be closest to reality. Table 10 shows the best strategy (with maximum EMV) for each 4DS outcome according to each reliability level. Note that, depending on the reliability of the information, a strategy different to that optimized for the RM itself could be chosen (strategies in bold). For instance, when the seismic data suggests RM4 with a reliability of 70%, strategy S3 is selected (for RM3, the best on average without further information) instead of strategy S4 (optimized for RM4). Most of these changes occur when reliability is between 70% and 50%. Thus, highly imperfect information (less than 70% reliable in our case study) is insufficient to alter the perception of the uncertain reservoir and thus inadequate to change decisions.
Best strategy (maximum EMV) for each 4DS outcome (represented by the representative models) according to each level of information reliability. Strategies in bold represent changes in the decision when compared to the case with perfect information.
6 Conclusion
This work presented the practical application of the CoS estimation in a real case, providing decisionmakers with a tool to decide whether or not to acquire information (in the case of this work, the information would be provided by 4D seismic). This is a modelbased study, where information identifies the representative model closest to the true reservoir model, and decisions correspond to implementing the production strategy optimized for that model.
While the traditional expected value of information calculation provides only an average value, the chance of success considers the variability in the economic return due to reservoir uncertainties. We assessed perfect information to obtain the upper limit for the VoI.
We also showed that, when information is imperfect, the decision of which production strategy to implement changes with information reliability. Moreover, when the reliability reaches a level of around 50%, the information no longer has value, since the economic return, in this case, is similar to that of the case without information. Furthermore, highly imperfect information (less than 70% reliable in our case study) is insufficient to change the perception of the uncertain reservoir and consequently decisions.
Acknowledgments
This work was carried out in association with the ongoing Project registered as BG07 “Reduction of uncertainties through the incorporation of 4D seismic data in the modeling of the reservoir” (UNICAMP/Shell Brazil/ANP) funded by Shell Brazil under the ANP R&D levy as “Investment Commitment to Research and Development”. The authors also thank Energi Simulation, CEPETRO, UNICAMP, DEFEM for supporting this work and CMG for software licenses.
References
 Abellan A., Noetinger B. (2010) Optimizing subsurface field data acquisition using information theory, Math. Geosci. 42, 603–630. https://doi.org/10.1007/s1100401092856. [CrossRef] [Google Scholar]
 Adlam J. (1995) The Norne Field Development Overview. Offshore Technology Conference, Houston, 1–4 May, OTC7925MS. https://doi.org/10.4043/7925MS [Google Scholar]
 Bailey W.J., Couët B., Prange M. (2011) Forecast optimization and value of information under uncertainty, in: Ma Y.Z., La Pointe P.R. (eds), Uncertainty analysis and reservoir modeling: AAPG Memoir 96, 217–233. doi: 10.1306/13301416M963484. [Google Scholar]
 Ballin P.R., Ward G.S., Whorlow C.V., Khan T. (2005) Value of information for a 4D seismic acquisition project. SPE 94918MS, Annual Technical Conference and Exhibition, Denver, Colorado, USA. https://doi.org/10.2118/94918MS. [Google Scholar]
 Barros E.G.D., Jansen J.D., van den Hof P.M.J. (2015) Value of information in parameter identification and optimization of hydrocarbon reservoirs, IFACPapers OnLine 48, 6, 229–235. https://doi.org/10.1016/j.ifacol.2015.08.036. [CrossRef] [Google Scholar]
 Barros E.G.D., van den Hof P.M.J., Jansen J.D. (2016) Value of information in closedloop reservoir management, Comp. Geosci. 20, 3, 737–749. https://doi.org/10.1007/s105960159509. [CrossRef] [Google Scholar]
 Begg S., Bratvold R., Campbell J. (2002) The value of flexibility in managing uncertainty in oil and gas investments. SPE 77586MS, Annual Technical Conference and Exhibition, San Antonio, Texas, USA. https://doi.org/10.2118/77586MS. [Google Scholar]
 Bratvold R.B., Bickel J.E., Lohne H.P. (2009) Value of information in the oil and gas industry: Past, present and future. SPE 110378PA, SPE Reserv. Eval. Eng. 12, 4, 630–638. https://doi.org/10.2118/110378PA. [Google Scholar]
 Chen B., He J., Wen X.H., Chen W., Reynolds A.C. (2017) Uncertainty quantification and value of information assessment using proxies and Markov chain Monte Carlo method for a pilot project, J. Petrol. Sci. Eng. 157, 328–339. https://doi.org/10.1016/j.petrol.2017.07.039. [CrossRef] [Google Scholar]
 Clemen R.T. (1996) Making hard decisions: An introduction to decision analysis, Second Edition, Duxbury Press, USA. [Google Scholar]
 Coopersmith E.M., Cunningham P.C. (2002) A practical approach to evaluating the value of information and real option decisions in the upstream petroleum industry. SPE 77582MS, Annual Technical Conference and Exhibition, San Antonio, Texas, USA. https://doi.org/10.2118/77582MS. [Google Scholar]
 Costa A.P.A., Schiozer D.J., Moczydlower P., Bedrikovetsky P. (2008) Use of representative models to improve the decision making process of chemical flooding in a mature field. SPE 115442, SPE Russian Oil & Gas Technical Conference and Exhibition, Moscow, Russia. https://doi.org/10.2118/115442MS. [Google Scholar]
 Davolio A., Schiozer D.J. (2018) Probabilistic seismic history matching using binary images, J. Geophys. 15, 2018, 261–274. https://doi.org/10.1088/17422140/aa99f4. [Google Scholar]
 Dunn M.D. (1992) A method to estimate the value of well log information. SPE 24672MS, Annual Technical Conference and Exhibition, Washington, D.C., USA. https://doi.org/10.2118/24672MS. [Google Scholar]
 Eidsvik J., Mukerji T., Bhattacharjya D. (2015) Value of information in the Earth sciences: Integrating spatial modeling and decision analysis, 1st edn., Cambridge University Press, Cambridge, United Kingdom. [CrossRef] [Google Scholar]
 Evensen G. (1994) Sequential data assimilation with a nonlinear quasigeostrophic model using Monte Carlo methods to forecast error statistics, J. Geophys. 99, C5, 10143–10162. https://doi.org/10.1029/94JC00572. [CrossRef] [Google Scholar]
 Ferreira C.J., Schiozer D.J. (2013) Timelapse seismic and engineering data integration to estimate best time for seismic acquisition data. SPE 164914MS, EAGE Annual Conference & Exhibition Incorporating Europec, London, United Kingdom. https://doi.org/10.2118/164914MS. [Google Scholar]
 Ferreira C.J., Schiozer D.J. (2014) Methodology to estimate the chance of success of a timelapse seismic project using reservoir simulation. SPE 170674MS, Annual Technical Conference and Exhibition, Amsterdam, The Netherlands. https://doi.org/10.2118/170674MS. [Google Scholar]
 Ferreira C.J., Davolio A., Maschio C., Schiozer D.J. (2011) Evaluation of the seismic time_lapse acquisition period in reservoir monitoring. OTC22740, Offshore Technology Conference Brasil, Rio de Janeiro, Brazil. https://doi.org/10.4043/22740MS. [Google Scholar]
 Ferreira C.J., Davolio A., Schiozer D.J. (2015) Improving the estimation of the chance of success of a 4D seismic project based on representative models. SPE 174386MS, SPE Europec, Madrid, Spain. https://doi.org/10.2118/174386MS. [Google Scholar]
 Galli A., Armstrong M., Dias M.A.G. (2004) The value of information: A Bayesian real option approach. SPE 90418MS, Annual Technical Conference and Exhibition, Houston, Texas, USA. https://doi.org/10.2118/90418MS. [Google Scholar]
 Gerhardt J.H., Haldorsen H.H. (1989) On the Value of Information. SPE 19291MS, Offshore Europe, Aberdeen, United Kingdom. https://doi.org/10.2118/19291MS. [Google Scholar]
 Grayson C.J. Jr. (1960) Decisions under uncertainty: Drilling decisions by oil and gas operators, Harvard University Press, Boston, Massachusetts, USA. [Google Scholar]
 He J., Xie J., Sarma P., Wen X.H., Chen W.H., Kamath J. (2016) Proxybased Work flow for a priori evaluation of dataacquisition programs, SPE J. 21, 4, SPE173229PA. https://doi.org/10.2118/173229PA. [Google Scholar]
 Hong A.J., Bratvold R.B., Thomas P., Hanea R.G. (2018) Valueofinformation for model parameter updating through history matching, J. Petrol. Sci. Eng. 165, 253–268. https://doi.org/10.1016/j.petrol.2017.08.028. [CrossRef] [Google Scholar]
 Howard R.A. (1966) Information value theory, IEEE T. Syst. Man. Cyb. 2, 1, 22–26, http://dx.doi.org/10.1109/TSSC.1966.300074. [Google Scholar]
 Le Ravalec M., Tillier E., Da Veiga S., Enchéry G., Gervais V. (2012a) Advanced integrated workflows for incorporating both production and 4D seismicrelated data into reservoir models, Oil Gas Sci. Technol.  Rev. IFP Energies nouvelles 67, 2, 207–220. https://doi.org/10.2516/ogst/2011159. [CrossRef] [Google Scholar]
 Le Ravalec M., Da Veiga S., Derfoul R., Enchéry G., Gervais V., Roggero F. (2012b) Integrating Data of Different Types and Different Supports into Reservoir Models, Oil Gas Sci. Technol.  Rev. IFP Energies nouvelles 67, 5, 823–839. https://doi.org/10.2516/ogst/2012024. [CrossRef] [Google Scholar]
 Lind O. (2004) Annual Reservoir Development Plan – Norne Field, Statoil Internal Report. [Google Scholar]
 Marques M.D., Gaspar A.T., Schiozer D.J. (2013) Use of oil reservoir simulation to estimate value of flexibility. SPE 164878, EAGE Annual Conference & Exhibition incorporating Europec, London, United Kingdom. https://doi.org/10.2118/164878MS [Google Scholar]
 Meira L.A.A., Coelho G.P., Santos A.A.S., Schiozer D.J. (2016) Selection of representative models for decision analysis under uncertainty, Comput. Geosci. 88, 67–82. https://doi.org/10.1016/j.cageo.2015.11.012 [CrossRef] [Google Scholar]
 Meira L.A.A., Coelho G.P., Silva C.G., Santos A.A.S., Schiozer D.J. (2017) RMFinder 2.0: An improved interactive multicriteria scenario reduction methodology. SPE 185502, SPE Latin America and Caribbean Petroleum Engineering Conference. https://doi.org/10.1016/j.cageo.2015.11.012 [Google Scholar]
 Raiffa H. (1968) Decision analysis: Introductory lectures on choices under uncertainty, AddisonWesley, Reading, Massachusetts, USA. [Google Scholar]
 Roggero F., Lerat O., Ding D.Y., Berthet P., Bordenave C., Lefeuvre F., Perfetti P. (2012) History matching of production and 4D seismic data: Application to the Girassol Field, Offshore Angola, Oil Gas Sci. Technol.  Rev. IFP Energies nouvelles 67, 2, 237–262. https://doi.org/10.2516/ogst/2011148 [CrossRef] [Google Scholar]
 Santos S.M.G., Gaspar A.T.F.S., Schiozer D.J. (2017) Value of information in reservoir development projects: Technical indicators to prioritize uncertainties and information sources, J. Petrol. Sci. Eng. 157, 1179–1191. https://doi.org/10.1016/j.petrol.2017.08.028 [CrossRef] [Google Scholar]
 Sarma P., Chen W.H., Xie J. (2013) Selecting representative models from a large set of models. SPE 163671, SPE Reservoir Simulation Symposium, The Woodlands, Texas, USA. https://doi.org/10.2118/163671MS. [Google Scholar]
 Schiozer D.J., Ligero E.L., Suslick S.B., Costa A.P.A., Santos J.A.M. (2004) Use of representative models in the integration of risk analysis and production strategy definition, J. Petrol. Sci. Eng. 44, 1–2, 131–141. https://doi.org/10.1016/j.petrol.2004.02.010 [CrossRef] [Google Scholar]
 Schlaifer R. (1959) Probability and statistics for business decisions: an introduction to managerial economics under uncertainty, McGrawHill Book Company, New York City, USA. [Google Scholar]
 Shirangi M.G., Durlofsky L.J. (2015) Closedloop field development under uncertainty by use of optimization with sample validation. SPE 173219PA, SPE J. 20, 5, 908–922. https://doi.org/10.2118/173219PA [CrossRef] [Google Scholar]
 Tillier E., Le Ravalec M., Da Veiga S. (2012) Simultaneous inversion of production data and seismic attributes: application to a synthetic SAGD produced field case, Oil Gas Sci. Technol.  Rev. IFP Energies nouvelles 67, 2, 289–301. https://doi.org/10.2516/ogst/2012004 [CrossRef] [Google Scholar]
 Waggoner J.R. (2002) Quantifying the Economic Impact of 4D Seismic Projects. SPE 77969PA, SPE Reserv. Eval. Eng. 5, 2, 111–115. https://doi.org/10.2118/77969PA [Google Scholar]
All Tables
Example of NPV for three RMs (RM1 to RM3, equiprobable prior to information acquisition) simulated for all production strategies (S1–S3) and respective EMV. S3 is selected as the best strategy without further information.
Example of NPV without information (NPV_{woi}), NPV with information (NPV_{wi}) and ΔNPV for the three RMs.
NPV, cumulative oil and water production, and the probability of occurrence of each RM.
Net Present Value without information (NPV_{woi}), Net Present Value with information (NPV_{wi}), ΔNPV and probability of occurrence (Prob.) for each strategy.
NPV of each RM (RM1 to RM9) simulated for all production strategies (S1–S9) and respective EMV. The bold values represent the NPV_{wi} of each strategy. The bold value in EMV line represents the maximum EMV among the strategies.
Example of P _{r} (reliability probability) for an information reliability of 0.9.
Example of updated probabilities (P _{u}) for an information reliability of 0.9.
Best strategy (maximum EMV) for each 4DS outcome (represented by the representative models) according to each level of information reliability. Strategies in bold represent changes in the decision when compared to the case with perfect information.
All Figures
Fig. 1. Methodology to estimate the chance of success of 4D seismic acquisition (adapted from Ferreira and Schiozer, 2014). 

In the text 
Fig. 2. Example of CoS analysis – ΔNPV × cumulative probability. 

In the text 
Fig. 3. (a) Norne segments (Lind, 2004) and (b) 3D view of porosity map for the simulation model. 

In the text 
Fig. 4. Cumulative (a) oil and (b) water production for all scenarios. 

In the text 
Fig. 5. Risk curves for (a) NPV and (b) Np. Purple dots represent all scenarios; black dots, the RMs. 

In the text 
Fig. 6. Crossplots: (a) NPV × Np, (b) NPV × Wp, (c) NPV × RF, and (d) Np × Wp. Purple dots represent all scenarios; black dots, the RMs. 

In the text 
Fig. 7. (a) Water saturation error curves and (b) normalized water saturation error curves. 

In the text 
Fig. 8. Histograms – number of scenarios that reached an Sw error of (a) 0.60 and (b) 0.50 in each year of production. 

In the text 
Fig. 9. (a) Histogram of number of wells (considering all scenarios) according to time of breakthrough (years), and (b) rate of total number of wells which reached breakthrough in each year of production. 

In the text 
Fig. 10. ΔNPV for each strategy. 

In the text 
Fig. 11. Compact decision tree for the calculation of the EVoI. 

In the text 
Fig. 12. CoS evaluation for perfect information: cumulative probability curve of the ΔNPV. 

In the text 
Fig. 13. Sensitivity analysis – CoS with different number of RMs. 

In the text 
Fig. 14. Simplified decision tree for EVoI calculation considering imperfect information with a reliability of 0.9. 

In the text 
Fig. 15. Variation of the EMV of the 4DS project with information reliability. 

In the text 
Fig. 16. Variation of EVII according to each level of reliability. 

In the text 