Reduced-Order Modeling of Coupled Flow and Quasistatic Geomechanics

SPE Journal ◽  
2019 ◽  
Vol 25 (01) ◽  
pp. 326-346 ◽  
Author(s):  
Zhaoyang Larry Jin ◽  
Timur Garipov ◽  
Oleg Volkov ◽  
Louis J. Durlofsky

Summary A reduced-order-modeling (ROM) framework is developed and applied to simulate coupled flow/quasistatic-geomechanics problems. The reduced-order model is constructed using POD-TPWL, in which proper orthogonal decomposition (POD), which enables representation of the solution unknowns in a low-dimensional subspace, is combined with trajectory-piecewise linearization (TPWL), where solutions with new sets of well controls are represented by means of linearization around previously simulated (training) solutions. The overdetermined system of equations is projected into the low-dimensional subspace using a least-squares Petrov-Galerkin (LSPG) procedure, which has been shown to maintain numerical stability in POD-TPWL models. The states and derivative matrices required by POD-TPWL, generated by an extended version of the Stanford University Automatic Differentiation General Purpose Research Simulator (AD-GPRS), are provided in an offline (preprocessing or training) step. Offline computational requirements correspond to the equivalent of five to eight full-order simulations, depending on the number of training runs used. Run-time (online) speedups of O(100) or more are typically achieved for new POD-TPWL test-case simulations. The POD-TPWL model is tested extensively for a 2D coupled problem involving oil/water flow and geomechanics. It is shown that POD-TPWL provides predictions of reasonable accuracy, relative to full-order simulations, for well-rate quantities, global pressure and saturation fields, global maximum- and minimum-principal-stress fields, and the Mohr-Coulomb rock-failure criterion, for the cases considered. A systematic study of POD-TPWL error is conducted using various training procedures for different levels of perturbation between test and training cases. The use of randomness in the well-bottomhole-pressure (BHP) profiles used in training is shown to be beneficial in terms of POD-TPWL solution accuracy. The procedure is also successfully applied to a prototype 3D example case.

SPE Journal ◽  
2014 ◽  
Vol 19 (05) ◽  
pp. 858-872 ◽  
Author(s):  
Jincong He ◽  
Louis J. Durlofsky

Summary Compositional simulation can be very demanding computationally as a result of the potentially large number of system unknowns and the intrinsic nonlinearity of typical problems. In this work, we develop a reduced-order modeling procedure for compositional simulation. The technique combines trajectory piecewise linearization (TPWL) and proper orthogonal decomposition (POD) to provide a highly efficient surrogate model. The compositional POD-TPWL method expresses new solutions in terms of linearizations around states generated (and saved) during previously simulated “training” runs. High-dimensional states are projected (optimally) into a low-dimensional subspace by use of POD. The compositional POD-TPWL model is based on a molar formulation that uses pressure and overall component mole fractions as the primary unknowns. Several new POD-TPWL treatments, including the use of a Petrov-Galerkin projection to reduce the number of equations (rather than the Galerkin projection, which was applied previously), and a new procedure for determining which saved state to use for linearization are incorporated into the method. Results are presented for heterogeneous 3D reservoir models containing oil and gas phases with up to six hydrocarbon components. Reasonably close agreement between full-order reference solutions and compositional POD-TPWL simulations is demonstrated for the cases considered. Construction of the POD-TPWL model requires preprocessing overhead computations equivalent to approximately three or four full-order runs. Runtime speedups by use of POD-TPWL are, however, very significant—up to a factor of 800 for the cases considered. The POD-TPWL model is thus well suited for use in computational optimization, in which many simulations must be performed, and we present an example demonstrating its application for such a problem.


2019 ◽  
Author(s):  
Zhaoyang Larry Jin ◽  
Timur Garipov ◽  
Oleg Volkov ◽  
Louis J. Durlofsky

Author(s):  
Jie Yuan ◽  
Fadi El-Haddad ◽  
Loic Salles ◽  
Chian Wong

This work presents an assessment of classical and state of the art reduced order modeling (ROM) techniques to enhance the computational efficiency for dynamic analysis of jointed structures with local contact nonlinearities. These ROM methods include classical free interface method (Rubin method, MacNeal method), fixed interface method Craig-Bampton (CB), Dual Craig-Bampton (DCB) method and also recently developed joint interface mode (JIM) and trial vector derivative (TVD) approaches. A finite element (FE) jointed beam model is considered as the test case taking into account two different setups: one with a linearized spring joint and the other with a nonlinear macroslip contact friction joint. Using these ROM techniques, the accuracy of dynamic behaviors and their computational expense are compared separately. We also studied the effect of excitation levels, joint region size, and number of modes on the performance of these ROM methods.


Author(s):  
Kazuto Hasegawa ◽  
Kai Fukami ◽  
Takaaki Murata ◽  
Koji Fukagata

Abstract We propose a reduced order model for predicting unsteady flows using a data-driven method. As preliminary tests, we use two-dimensional unsteady flow around bluff bodies with different shapes as the training datasets obtained by direct numerical simulation (DNS). Our machine-learned architecture consists of two parts: Convolutional Neural Network-based AutoEncoder (CNN-AE) and Long Short Term Memory (LSTM), respectively. First, CNN-AE is used to map into a low-dimensional space from the flow field data. Then, LSTM is employed to predict the temporal evolution of the low-dimensional data generated by CNN-AE. Proposed machine-learned reduced order model is applied to two-dimensional circular cylinder flows at various Reynolds numbers and flows around bluff bodies of various shapes. The flow fields reconstructed by the machine-learned architecture show reasonable agreement with the reference DNS data. Furthermore, it can be seen that our machine-learned reduced order model can successfully map the high-dimensional flow data into low-dimensional field and predict the flow fields against unknown Reynolds number fields and shapes of bluff body. As concluding remarks, we discuss the extension study of machine-learned reduced order modeling for various applications in experimental and computational fluid dynamics.


2014 ◽  
Author(s):  
Donald L. Brown ◽  
Jun Li ◽  
Victor M. Calo ◽  
Mehdi Ghommem ◽  
Yalchin Efendiev

Sign in / Sign up

Export Citation Format

Share Document