Figure 1.
Description of Random Forests model where, each of the branches from sample dataset correspond to a bootstrap sample and a random feature subset is chosen for each branch to construct a decision tree.
Figure 1.
Description of Random Forests model where, each of the branches from sample dataset correspond to a bootstrap sample and a random feature subset is chosen for each branch to construct a decision tree.
Figure 2.
(a) Heterogeneous channelized permeability field with one producer in red and one injector in white (b) Streamlines imposed on the permeability field.
Figure 2.
(a) Heterogeneous channelized permeability field with one producer in red and one injector in white (b) Streamlines imposed on the permeability field.
Figure 3.
(a) The total travel time of the particle injected at injector to any point in the reservoir (b) F — Φ plot to compute Lorenz coefficient as the area under the solid F — Φ curve and the dashed line F = Φ.
Figure 3.
(a) The total travel time of the particle injected at injector to any point in the reservoir (b) F — Φ plot to compute Lorenz coefficient as the area under the solid F — Φ curve and the dashed line F = Φ.
Figure 4.
(a) Homogeneous permeability reservoir with 1 injector (white) and 1 producer (red). Producer is the parameter and can change location anywhere in the reservoir (b) BHP profiles of injector and producer that is set constant for all well configurations considered.
Figure 4.
(a) Homogeneous permeability reservoir with 1 injector (white) and 1 producer (red). Producer is the parameter and can change location anywhere in the reservoir (b) BHP profiles of injector and producer that is set constant for all well configurations considered.
Figure 5.
Change in cumulative energy of the basis, as calculated by magnitude of eigenvalues for increasing number of parameters (well configurations) (a) Nζ = 10, (b) Nζ = 200, (c) Nζ = 400, in the state snapshot matrix.
Figure 5.
Change in cumulative energy of the basis, as calculated by magnitude of eigenvalues for increasing number of parameters (well configurations) (a) Nζ = 10, (b) Nζ = 200, (c) Nζ = 400, in the state snapshot matrix.
Figure 6.
(a) Training samples of the producer well locations shown by red crosses. (b) Test case 1 well location at (3,13) in a homogeneous reservoir.
Figure 6.
(a) Training samples of the producer well locations shown by red crosses. (b) Test case 1 well location at (3,13) in a homogeneous reservoir.
Figure 7.
Pressure solution comparison at (a) Time = 70 days and (b) Time = 360 days.
Figure 7.
Pressure solution comparison at (a) Time = 70 days and (b) Time = 360 days.
Figure 8.
Saturation solution comparison at (a) Time = 70 days and (b) Time = 360 days.
Figure 8.
Saturation solution comparison at (a) Time = 70 days and (b) Time = 360 days.
Figure 9.
Relative error in pressure at time = (a) 1 day, (b) 70 days and (c) 360 days.
Figure 9.
Relative error in pressure at time = (a) 1 day, (b) 70 days and (c) 360 days.
Figure 10.
Relative error in saturation at time = (a) 1 day, (b) 70 days and (c) 360 days.
Figure 10.
Relative error in saturation at time = (a) 1 day, (b) 70 days and (c) 360 days.
Figure 11.
Well block state solution comparison (a) Pressure and (b) Saturation.
Figure 11.
Well block state solution comparison (a) Pressure and (b) Saturation.
Figure 12.
Quantities of Interest comparison (a) Oil production rate and (b) Water cut.
Figure 12.
Quantities of Interest comparison (a) Oil production rate and (b) Water cut.
Figure 13.
Test case 2 with producer at (8,4) grid block and injector at (20,1) grid block in homogeneous permeability reservoir.
Figure 13.
Test case 2 with producer at (8,4) grid block and injector at (20,1) grid block in homogeneous permeability reservoir.
Figure 14.
Pressure solution comparison at (a) Time = 70 days and (b) Time = 360 days.
Figure 14.
Pressure solution comparison at (a) Time = 70 days and (b) Time = 360 days.
Figure 15.
Saturation solution comparison at (a) Time = 70 days and (b) Time = 360 days.
Figure 15.
Saturation solution comparison at (a) Time = 70 days and (b) Time = 360 days.
Figure 16.
Well block state solution comparison (a) Pressure and (b) Saturation.
Figure 16.
Well block state solution comparison (a) Pressure and (b) Saturation.
Figure 17.
Quantities of Interest comparison (a) Oil production rate and (b) Water cut.
Figure 17.
Quantities of Interest comparison (a) Oil production rate and (b) Water cut.
Figure 18.
True and ML predicted pressure basis coefficient comparison at time = (a) 70 days and (b) 360 days.
Figure 18.
True and ML predicted pressure basis coefficient comparison at time = (a) 70 days and (b) 360 days.
Figure 19.
True and ML predicted saturation basis coefficient comparison at time = (a) 70 days and (b) 360 days.
Figure 19.
True and ML predicted saturation basis coefficient comparison at time = (a) 70 days and (b) 360 days.
Figure 20.
Comparison of predicted orthogonal and true orthogonal solutions with true solution at time = 360 days for (a) Pressure and (b) Saturation.
Figure 20.
Comparison of predicted orthogonal and true orthogonal solutions with true solution at time = 360 days for (a) Pressure and (b) Saturation.
Figure 21.
Schematic of a feedforward Artificial Neural Network.
Figure 21.
Schematic of a feedforward Artificial Neural Network.
Figure 22.
Workflow of the ML-based predict and correct PMOR procedure to estimate well production rates for new well locations.
Figure 22.
Workflow of the ML-based predict and correct PMOR procedure to estimate well production rates for new well locations.
Figure 23.
(a) Test case with a producer well location at (2,6). (b) Comparison of oil production rate and (c) comparison of water cut, predicted using the global PMOR method alone using 100 samples (dotted green line) and after the implementation of the error correction model (red circled line) with the true solution (blue line).
Figure 23.
(a) Test case with a producer well location at (2,6). (b) Comparison of oil production rate and (c) comparison of water cut, predicted using the global PMOR method alone using 100 samples (dotted green line) and after the implementation of the error correction model (red circled line) with the true solution (blue line).
Figure 24.
(a) Test case with a producer well location at (3,16). (b) Comparison of the oil production rate and (c) comparison of water cut, predicted using the global PMOR method alone (dotted green line) and after the implementation of the error correction model (red circled line) with the true solution (blue line).
Figure 24.
(a) Test case with a producer well location at (3,16). (b) Comparison of the oil production rate and (c) comparison of water cut, predicted using the global PMOR method alone (dotted green line) and after the implementation of the error correction model (red circled line) with the true solution (blue line).
Figure 25.
(a) Test case with a producer well location at (16,19). (b) Comparison of the oil production rate and (c) comparison of water cut, predicted using global the PMOR method alone using 100 samples (dotted green line) and after the implementation of error correction model (red circled line) with the true solution (blue line).
Figure 25.
(a) Test case with a producer well location at (16,19). (b) Comparison of the oil production rate and (c) comparison of water cut, predicted using global the PMOR method alone using 100 samples (dotted green line) and after the implementation of error correction model (red circled line) with the true solution (blue line).
Figure 26.
Error in the prediction of (a) the oil production rate and (b) water cut for all the test cases before and after the error correction of the solutions.
Figure 26.
Error in the prediction of (a) the oil production rate and (b) water cut for all the test cases before and after the error correction of the solutions.
Figure 27.
(a) Test case with producer well location at (3,16). (b) aomparison of oil production rate and (c) comparison of water cut, predicted using the global PMOR method alone using 50 samples (dotted green line) and after the implementation of error correction model (red circled line) with the true solution (blue line).
Figure 27.
(a) Test case with producer well location at (3,16). (b) aomparison of oil production rate and (c) comparison of water cut, predicted using the global PMOR method alone using 50 samples (dotted green line) and after the implementation of error correction model (red circled line) with the true solution (blue line).
Figure 28.
(a) Test case with producer well location at (11,10). (b) Comparison of oil production rate and (c) comparison of water cut, predicted the using global PMOR method alone using 50 samples (dotted green line) and after the implementation of error correction model (red circled line) with the true solution (blue line).
Figure 28.
(a) Test case with producer well location at (11,10). (b) Comparison of oil production rate and (c) comparison of water cut, predicted the using global PMOR method alone using 50 samples (dotted green line) and after the implementation of error correction model (red circled line) with the true solution (blue line).
Figure 29.
Error in the prediction of (a) the oil production rate and (b) water cut for all the test cases before and after the error correction of the solutions.
Figure 29.
Error in the prediction of (a) the oil production rate and (b) water cut for all the test cases before and after the error correction of the solutions.
Figure 30.
Heterogeneous log permeability field (section of SPE10 model layer 50) with one producer and one injector.
Figure 30.
Heterogeneous log permeability field (section of SPE10 model layer 50) with one producer and one injector.
Figure 31.
Samples considered as potential well locations shown by red crosses. Injector location is fixed at location (10,1) (not shown here) and hence not sampled.
Figure 31.
Samples considered as potential well locations shown by red crosses. Injector location is fixed at location (10,1) (not shown here) and hence not sampled.
Figure 32.
(a) Training samples randomly chosen for basis coefficient prediction. (b) Cumulative energy of eigenvalues for pressure (red) and saturation basis (blue).
Figure 32.
(a) Training samples randomly chosen for basis coefficient prediction. (b) Cumulative energy of eigenvalues for pressure (red) and saturation basis (blue).
Figure 33.
Test case with producer well at grid block (28,50).
Figure 33.
Test case with producer well at grid block (28,50).
Figure 34.
Pressure solution comparison at (a) Time = 506 days and (b) Time = 1046 days.
Figure 34.
Pressure solution comparison at (a) Time = 506 days and (b) Time = 1046 days.
Figure 35.
Saturation solution comparison at (a) Time = 506 days and (b) Time = 1046 days.
Figure 35.
Saturation solution comparison at (a) Time = 506 days and (b) Time = 1046 days.
Figure 36.
Relative error in pressure at time = (a) 506 days (b) 1046 days.
Figure 36.
Relative error in pressure at time = (a) 506 days (b) 1046 days.
Figure 37.
True and ML predicted pressure basis coefficient comparison at time = (a) 506 days and (b) 1046 days.
Figure 37.
True and ML predicted pressure basis coefficient comparison at time = (a) 506 days and (b) 1046 days.
Figure 38.
True and ML predicted pressure basis coefficient comparison at time = (a) 506 days and (b) 1046 days.
Figure 38.
True and ML predicted pressure basis coefficient comparison at time = (a) 506 days and (b) 1046 days.
Figure 39.
Comparison of the relative error in pressure due to orthogonal projection (left) and ML-predicted solution (right).
Figure 39.
Comparison of the relative error in pressure due to orthogonal projection (left) and ML-predicted solution (right).
Figure 40.
(a) Test case 1 with producer well location at (36,26). (b) Comparison of oil production rate and (c) comparison of water cut, predicted using the global PMOR method alone (dotted green line) and after the implementation of error correction model (red circled line) with the true solution (blue line).
Figure 40.
(a) Test case 1 with producer well location at (36,26). (b) Comparison of oil production rate and (c) comparison of water cut, predicted using the global PMOR method alone (dotted green line) and after the implementation of error correction model (red circled line) with the true solution (blue line).
Figure 41.
(a) Test case 2 with the producer well location at (28,50). (b) Comparison of the oil production rate and (c) comparison of water cut, predicted using the global PMOR method alone (dotted green line) and after the implementation of error correction model (red circled line) with the true solution (blue line).
Figure 41.
(a) Test case 2 with the producer well location at (28,50). (b) Comparison of the oil production rate and (c) comparison of water cut, predicted using the global PMOR method alone (dotted green line) and after the implementation of error correction model (red circled line) with the true solution (blue line).
Figure 42.
(a) Test case 3 with producer well location at (14,12). (b) Comparison of oil production rate and (c) comparison of water cut, predicted using the global PMOR method alone (dotted green line) and after the implementation of error correction model (red circled line) with the true solution (blue line).
Figure 42.
(a) Test case 3 with producer well location at (14,12). (b) Comparison of oil production rate and (c) comparison of water cut, predicted using the global PMOR method alone (dotted green line) and after the implementation of error correction model (red circled line) with the true solution (blue line).
Figure 43.
Error in the prediction of (a) the oil production rate (b) water cut for all the test cases before and after the error correction of the solutions.
Figure 43.
Error in the prediction of (a) the oil production rate (b) water cut for all the test cases before and after the error correction of the solutions.
Figure 44.
Box plot of the accuracies using the predicted POD coefficients of (a) pressure for 95%, 97% and 99% energy of pressure singular values and (b) saturation for 80%, 85% and 95% energy of saturation singular values.
Figure 44.
Box plot of the accuracies using the predicted POD coefficients of (a) pressure for 95%, 97% and 99% energy of pressure singular values and (b) saturation for 80%, 85% and 95% energy of saturation singular values.
Figure 45.
(a) Test case 1 with the producer well location at (9,29). (b) Comparison of the oil production rate and (c) comparison of water cut, predicted using the global PMOR method alone (dotted green line) and after the implementation of error correction model (red circled line) with the true solution (blue line).
Figure 45.
(a) Test case 1 with the producer well location at (9,29). (b) Comparison of the oil production rate and (c) comparison of water cut, predicted using the global PMOR method alone (dotted green line) and after the implementation of error correction model (red circled line) with the true solution (blue line).
Figure 46.
(a) Test case 2 with the producer well location at (26,48). (b) Comparison of oil production rate and (c) comparison of water cut, predicted using the global PMOR method alone (dotted green line) and after the implementation of the error correction model (red circled line) with the true solution (blue line).
Figure 46.
(a) Test case 2 with the producer well location at (26,48). (b) Comparison of oil production rate and (c) comparison of water cut, predicted using the global PMOR method alone (dotted green line) and after the implementation of the error correction model (red circled line) with the true solution (blue line).
Figure 47.
Error in the prediction of (a) the oil production rate and (b) water cut for all the test cases before and after the error correction of the solutions.
Figure 47.
Error in the prediction of (a) the oil production rate and (b) water cut for all the test cases before and after the error correction of the solutions.
Figure 48.
Feature Rank as calculated by the Random Forest model.
Figure 48.
Feature Rank as calculated by the Random Forest model.
Figure 49.
Reservoir model with an injector at location (24,26) and red regions showing the regions where the producer 1 (top of the reservoir) and producer 2 (bottom of the reservoir) locations are sampled.
Figure 49.
Reservoir model with an injector at location (24,26) and red regions showing the regions where the producer 1 (top of the reservoir) and producer 2 (bottom of the reservoir) locations are sampled.
Figure 50.
(a) Test case 1 with producer 1 at (9,2) and producer 2 at (31,50). (b) Comparison of oil production rate for producer 1, (c) comparison of water cut for producer 1, (d) comparison of oil production rate for producer 2, and (e) comparison of water cut for producer 2, predicted using the global PMOR method alone (dotted green line) and after the implementation of the error correction model (red circled line) with the true solution (blue line).
Figure 50.
(a) Test case 1 with producer 1 at (9,2) and producer 2 at (31,50). (b) Comparison of oil production rate for producer 1, (c) comparison of water cut for producer 1, (d) comparison of oil production rate for producer 2, and (e) comparison of water cut for producer 2, predicted using the global PMOR method alone (dotted green line) and after the implementation of the error correction model (red circled line) with the true solution (blue line).
Figure 51.
(a) Test case 2 with producer 1 at (20,3) and producer 2 at (32,49). (b) Comparison of oil production rate for producer 1, (c) comparison of water cut for producer 1, (d) comparison of oil production rate for producer 2, (e) comparison of water cut for producer 2, predicted using the global PMOR method alone (dotted green line) and after the implementation of the error correction model (red circled line) with the true solution (blue line).
Figure 51.
(a) Test case 2 with producer 1 at (20,3) and producer 2 at (32,49). (b) Comparison of oil production rate for producer 1, (c) comparison of water cut for producer 1, (d) comparison of oil production rate for producer 2, (e) comparison of water cut for producer 2, predicted using the global PMOR method alone (dotted green line) and after the implementation of the error correction model (red circled line) with the true solution (blue line).
Figure 52.
Error in the prediction of (a) the oil production rate and (b) water cut for producer 1 and (c) oil production rate and (d) water cut for producer 2 for all the test cases before and after the error correction of the solutions.
Figure 52.
Error in the prediction of (a) the oil production rate and (b) water cut for producer 1 and (c) oil production rate and (d) water cut for producer 2 for all the test cases before and after the error correction of the solutions.
Table 1.
Geometric and physics-based features for ML model construction corresponding to the well configuration.
Table 1.
Geometric and physics-based features for ML model construction corresponding to the well configuration.
Feature Set for a Well Configuration |
---|
x—X coordinate of the well |
y—Y coordinate of the well |
r—Distance between well and injector |
—Angle between well and injector |
K—Permeability at the well location |
—Total time of flight at the producer well location |
—Lorenz coefficient for the well configuration |
—Well grid block number in the reservoir |
t—Time at which the POD coefficients are computed |
Table 2.
Hyperparameters chosen by 5-fold cross-validation for Random Forest Regressor using 100 training samples.
Table 2.
Hyperparameters chosen by 5-fold cross-validation for Random Forest Regressor using 100 training samples.
| RF Regression | Train Accuracy | Test Accuracy |
---|
Pressure | = 3, = 2 | 99.83 | 98.18 |
Saturation | = 3, = 3 | 98.89 | 93.21 |
Table 3.
Homogeneous reservoir case 1: the average accuracy of the oil production rate and water cut for all test samples.
Table 3.
Homogeneous reservoir case 1: the average accuracy of the oil production rate and water cut for all test samples.
| Average Accuracy (%) |
---|
Oil production rate | 95.28 |
Water cut | 98.4 |
Table 4.
Error NN model Hyperparameters found through grid search for homogeneous reservoir case.
Table 4.
Error NN model Hyperparameters found through grid search for homogeneous reservoir case.
Number of Nodes | L2 Regularization Parameter |
---|
8 | 0.1 |
Table 5.
Hyperparameters chosen by 5-fold cross-validation for Random Forest Regressor using 50 training samples.
Table 5.
Hyperparameters chosen by 5-fold cross-validation for Random Forest Regressor using 50 training samples.
| RF Regression | Train Accuracy | Test Accuracy |
---|
Pressure | = 2, = 2 | 99.55 | 97.55 |
Saturation | = 2, = 2 | 98.71 | 90.79 |
Table 6.
Homogeneous reservoir case 2: Average accuracy of the oil production rate and water cut for all test samples.
Table 6.
Homogeneous reservoir case 2: Average accuracy of the oil production rate and water cut for all test samples.
| Average Accuracy (%) |
---|
Oil production rate | 95.2 |
Water cut | 98.3 |
Table 7.
Error NN model hyperparameters found through grid search for the heterogeneous model with higher basis dimensions in the first step.
Table 7.
Error NN model hyperparameters found through grid search for the heterogeneous model with higher basis dimensions in the first step.
Number of Modes | L2 Regularization Parameter |
---|
10 | 0.4 |
Table 8.
Heterogeneous reservoir case 1: average accuracy of oil production rate and water cut for all test samples.
Table 8.
Heterogeneous reservoir case 1: average accuracy of oil production rate and water cut for all test samples.
| Average Accuracy (%) |
---|
Oil production rate | 92.6 |
Water cut | 95.04 |
Table 9.
Time (seconds) comparison for the two test cases between fine-scale simulation and reduced-order model with error correction.
Table 9.
Time (seconds) comparison for the two test cases between fine-scale simulation and reduced-order model with error correction.
| Fine-Scale Simulation | PMOR + Error Correction |
---|
Test Case 1 | 65 s | 1.3 s |
Test Case 2 | 52 s | 1 s |
Table 10.
Error NN model hyperparameters found through grid search for the heterogeneous model with lower basis dimensions in the first step.
Table 10.
Error NN model hyperparameters found through grid search for the heterogeneous model with lower basis dimensions in the first step.
Number of nodes | L2 Regularization Parameter |
---|
10 | 0.2 |
Table 11.
Heterogeneous reservoir case 2: average accuracy of the oil production rate and water cut for all test samples.
Table 11.
Heterogeneous reservoir case 2: average accuracy of the oil production rate and water cut for all test samples.
| Average Accuracy (%) |
---|
Oil production rate | 91.74 |
Water cut | 94.6 |
Table 12.
Time (seconds) comparison for the two test cases (lower-dimensional basis) between fine-scale simulation and reduced-order model with error correction.
Table 12.
Time (seconds) comparison for the two test cases (lower-dimensional basis) between fine-scale simulation and reduced-order model with error correction.
| Fine Scale Simulation | PMOR + Error Correction |
---|
Test Case 1 | 56 s | 0.7 s |
Test Case 2 | 46 s | 0.5 s |
Table 13.
Time (seconds) comparison for the two test cases between fine-scale simulation and reduced-order model with error correction.
Table 13.
Time (seconds) comparison for the two test cases between fine-scale simulation and reduced-order model with error correction.
| Fine Scale Simulation | PMOR + Error Correction |
---|
Test Case 1 | 102 s | 1 s |
Test Case 2 | 90 s | 0.8 s |