Assisted history matching workflow for unconventional reservoirs
The information of fractures geometry and reservoir properties can be retrieved from the production data, which is always available at no additional cost. However, in unconventional reservoirs, it is insufficient to obtain only one realization because the non-uniqueness of history matching and subsurface uncertainties cannot be captured. Therefore, the objective of this study is to obtain multiple realizations in shale reservoirs by adopting Assisted History Matching (AHM). We used multiple proxy-based Markov Chain Monte Carlo (MCMC) algorithm and Embedded Discrete Fracture Model (EDFM) to perform AHM. The reason is that MCMC has benefits of quantifying uncertainty without bias or being trapped in any local minima. Also, using MCMC with proxy model unlocks the limitation of an infeasible number of simulations required by a traditional MCMC algorithm. For fractures modeling, EDFM can mimic fractures flow behavior with a higher computational efficiency than a traditional local grid refinement (LGR) method and more accuracy than the continuum approach. We applied the AHM workflow to actual shale gas wells. We found that the algorithm can find multiple history matching solutions and quantify the fractures and reservoir properties posterior distributions. Then, we predicted the production probabilistically. Moreover, we investigated the performance of neural network (NN) and k-nearest neighbors (KNN) as a proxy model in the proxy-based MCMC algorithm. We found that NN performed better in term of accuracy than KNN but NN required twice running time of KNN. Lastly, we studied the effect of enhanced permeability area (EPA) and natural fractures existence on the history matching solutions and production forecast. We concluded that we would over-predict fracture geometries and properties and estimated ultimate recovery (EUR) if we assumed no EPA or no natural fractures even though they actually existed. The degree of over-prediction depends on fractures and reservoir properties, EPA and natural fractures properties, which can only be quantified after performing AHM. The benefits from this study are that we can characterize fractures geometry, reservoir properties, and natural fractures in a probabilistic manner. These multiple realizations can be further used for a probabilistic production forecast, future fracturing design improvement, and infill well placement decision.