Abstract:Accurate and efficient numerical simulation of unconventional reservoirs is challenging. Long periods of transient flow and steep potential gradients occur due to the extreme conductivity contrast between matrix and fracture. Detailed near-well/near-fracture models are necessary to provide sufficient resolution, but they are computationally impractical for field cases with multiple hydraulic-fracture stages.Previous works in the literature of unconventional simulations mainly focus on the gridding level that adapts to wells and fractures. Limited research has been conducted on nonlinear strategies that exploit locality across timesteps and nonlinear iterations. It was reported that an individual Newton update is typically sparse and nonlinear convergence is constrained by a small portion of the model. To perform localized computations, an <math>a-priori</math>strategy is essential to first determine the active subset of simulation cells for the subsequent iteration. The active set flags the cells that will be updated, and then the corresponding localized linear system is solved.The objective of this work is to develop localization methods that are readily applicable to complex fracture networks and flow physics in unconventional reservoirs. By utilizing the diffusive nature of pressure updates, an adaptive algorithm is proposed to make adequate estimates for the active domains. In addition, we further develop a localized solver based on nonlinear domain decomposition (DD). Compared to a standard DD method, domain partitions are dynamically constructed. The new solver provides effective partitioning that adapts to flow dynamics and Newton updates.We evaluate the developed methods using several complex problems with discrete fracture networks. The problems consider multi-phase and compositional fluid systems with phase changes. The results show that large degrees of solution locality present across timesteps and iterations. Compared to a standard Newton solver, the new solvers enable superior computational performance. Moreover, Newton convergence behavior is preserved, without any impact on solution accuracy.

A Pod-Based Fast Algorithm For The Nonlocal Unsteady Problems

Fast Model Order Reduction Via Nonlinear Optimization

A POD based reduced-order local RBF collocation approach for time-dependent nonlocal diffusion problems

Proper Orthogonal Decomposition Method to Nonlinear Filtering Problems in Medium-High Dimension

An efficient proximal-based approach for solving nonlocal Allen-Cahn equations

Fast Algorithms for Finite Element nonlinear Discrete Systems to Solve the PNP Equations

A Reduced-Order FEM Based on POD for Solving Non-Fourier Heat Conduction Problems under Laser Heating

A spectral collocation method for nonlocal diffusion equations

A local character based method for solving linear systems of radiation diffusion problems

A fast finite volume method for spatial fractional diffusion equations on nonuniform meshes

An Ensemble-Proper Orthogonal Decomposition Method for the Nonstationary Navier-Stokes Equations

Fast numerical algorithm for the reaction-diffusion equations using an interpolating method

Localized nonlinear solution strategies for efficient simulation of unconventional reservoirs

Applying Proper Orthogonal Decomposition to Parabolic Equations: A Reduced Order Numerical Approach

Accelerating the solution of linear systems appearing in two-phase reservoir simulation by the use of POD-based deflation methods

Stability and Error Analysis for a Second-Order Fast Approximation of the Local and Nonlocal Diffusion Equations on the Real Line

An accelerated novel meshless coupled algorithm for non-local nonlinear behavior in 2D/3D space-fractional GPEs

A Non-Nested Multilevel Method for Meshless Solution of the Poisson Equation in Heat Transfer and Fluid Flow

A new alternating iteration strategy based on the proper orthogonal decomposition for solving large-scaled transient nonlinear heat conduction problems

A PARALLEL SOLVING DESIGN FOR UNSTEADY NAVIER-STOKES EQUATION

An Augmented Subspace Based Adaptive Proper Orthogonal Decomposition Method for Time Dependent Partial Differential Equations