Robust Life-Cycle Production Optimization With a Support-Vector-Regression Proxy

We design a new and general work flow for efficient estimation of the optimal well controls for the robust production-optimization problem using support-vector regression (SVR), where the cost function is the net present value (NPV). Given a set of simulation results, an SVR model is built as a proxy to approximate a reservoir-simulation model, and then the estimated optimal controls are found by maximizing NPV using the SVR proxy as the forward model. The gradient of the SVR model can be computed analytically so the steepest-ascent algorithm can easily and efficiently be applied to maximize NPV. Then, the well-control optimization is performed using an SVR model as the forward model with a steepest-ascent algorithm. To the best of our knowledge, this is the first SVR application to the optimal well-control problem. We provide insight and information on proper training of the SVR proxy for life-cycle production optimization. In particular, we develop and implement a new iterative-sampling-refinement algorithm that is designed specifically to promote the accuracy of the SVR model for robust production optimization. One key observation that is important for reservoir optimization is that SVR produces a high-fidelity model near an optimal point, but at points far away, we only need SVR to produce reasonable approximations of the predicting output from the reservoir-simulation model. Because running an SVR model is computationally more efficient than running a full-scale reservoir-simulation model, the large computational cost spent on multiple forward-reservoir-simulation runs for robust optimization is significantly reduced by applying the proposed method. We compare the performance of the proposed method using the SVR runs with the popular stochastic simplex approximate gradient (StoSAG) and reservoir-simulations runs for three synthetic examples, including one field-scale example. We also compare the optimization performance of our proposed method with that obtained from a linear-response-surface model and multiple SVR proxies that are built for each of the geological models.