Abstract:In this study, a fast multipole method (FMM) is used to decrease the computational time of a fully-coupled poroelastic hydraulic fracture model with a controllable effect on its accuracy. The hydraulic fracture model is based on the poroelastic formulation of the displacement discontinuity method (DDM) which is a special formulation of the boundary element method (BEM). DDM is a powerful and efficient method for problems involving fractures. However, this method becomes slow as the number of temporal, or spatial elements increases, or necessary details such as poroelasticity, that makes the solution history-dependent, are added to the model. FMM is a technique to expedite matrix-vector multiplications within a controllable error without forming the matrix explicitly. Fully-coupled poroelastic formulation of DDM involves the multiplication of a dense matrix with a vector in several places. A crucial modification to DDM is suggested in two places in the algorithm to leverage the speed efficiency of FMM for carrying out these multiplications. The first modification is in the time-marching scheme, which accounts for the solution of previous time steps to compute the current time step. The second modification is in the generalized minimal residual method (GMRES) to iteratively solve for the problem unknowns. Several examples are provided to show the efficiency of the proposed approach in problems with large degrees of freedom (in time and space). Examples include hydraulic fracturing of a horizontal well and randomly distributed pressurized fractures at different orientations with respect to horizontal stresses. The results are compared to the conventional DDM in terms of computational processing time and accuracy. Accordingly, the proposed algorithm may be used for fracture propagation studies while substantially reducing the processing time with a controllable error.

What problem does this paper attempt to address?

The main problem that this paper attempts to solve is to improve the computational efficiency of the fully - coupled poroelastic hydraulic fracturing model. Specifically, the researchers use the Fast Multipole Method (FMM) to reduce the computation time while maintaining the accuracy of the model. The hydraulic fracturing model is based on the poroelastic formula of the Displacement Discontinuity Method (DDM), which is a special Boundary Element Method (BEM) formula and is especially suitable for problems involving fractures. However, as the number of time and space elements increases, or when details such as poroelasticity are added to the model, the DDM method becomes slow. FMM is a technique that can accelerate matrix - vector multiplication within a controllable error without explicitly forming a matrix. By introducing important modifications to DDM at two key positions in the algorithm to take advantage of the speed advantage of FMM for these multiplications, the computational efficiency is significantly improved. The paper demonstrates the effectiveness of the proposed method in dealing with problems with a large number of degrees of freedom (time and space) through several examples, such as the hydraulic fracturing of horizontal wells and pressure fractures with different directional distributions relative to the horizontal stress. The results are compared with those of the traditional DDM in terms of computation processing time and accuracy. The research shows that for a hydraulic fracturing example with 20,000 constant space elements and one time element, FMM can reduce the computation time by up to 70 times with a relative error of less than 4%; while for the randomly distributed pressure fracture case, the relative error is less than 0.5%. In addition, the crack - tip displacements calculated using the two methods are further used to compare the Stress Intensity Factors (SIF) in Mode I and Mode II, and the results show that the error in calculating SIF using the proposed modified method is almost negligible and thus does not affect the estimation of the crack propagation direction. Therefore, the proposed algorithm can be used for crack propagation research while significantly reducing the processing time.

Application of the Fast Multipole Fully Coupled Poroelastic Displacement Discontinuity Method to Hydraulic Fracturing Problems

Coupled Hydro-Mechanical XFEM Analysis for Multi-Fracturing Through an Excavation Driven by an Underlying Aquifer: a Forensic Case Study

Numerical Simulations of the Hydraulic Fracture Propagation in Poroelastic Media Using the Coupled Hydro-Mechanical Field-Enriched Finite Element Method

A 2D, Fully-Coupled, Hydro-Mechanical, FDEM Formulation for Modelling Fracturing Processes in Discontinuous, Porous Rock Masses

A Fully Coupled Hydro-Mechanical Approach for Multi-Fracture Propagation Simulations

Simulation of dynamic pulsing fracking in poroelastic media by a hydro-damage-mechanical coupled cohesive phase field model

A hydro-mechanical-damage fully coupled cohesive phase field model for complicated fracking simulations in poroelastic media

A Hybrid Embedded Discrete Fracture Model and Dual-Porosity, Dual-Permeability Workflow for Hierarchical Treatment of Fractures in Practical Field Studies

Coupled hydro-mechanical two-phase flow model in fractured porous medium with the combined finite-discrete element method

Coupled Geomechanics and Fluid Flow Modeling for Petroleum Reservoirs Accounting for Multi-Scale Fractures

Multiscale fractures integrated equivalent porous media method for simulating flow and solute transport in fracture-matrix system

Embedded Fracture Model for Coupled Flow and Geomechanics

A computational dual‐porosity approach for the coupled hydro‐mechanical analysis of fractured porous media

Numerical Study of Simultaneous Growth of Multiple Hydraulic Fractures from a Horizontal Wellbore Combining Dual Boundary Element Method and Finite Volume Method

Fully coupled dynamic hydraulic fracturing of saturated porous media based on the numerical manifold method

Superposition‐based Concurrent Multiscale Approaches for Porodynamics

A Multi-Field and Fluid–solid Coupling Method for Porous Media Based on DEM-PNM

Matrix‐based implementation and GPU acceleration of hybrid FEM and peridynamic model for hydro‐mechanical coupled problems

A Hybrid Finite Volume and Extended Finite Element Method for Hydraulic Fracturing with Cohesive Crack Propagation in Quasi-Brittle Materials

Simulation of coupled multiphase flow and geomechanics in porous media with embedded discrete fractures

Coupling mixed hybrid and extended finite element methods for the simulation of hydro-mechanical processes in fractured porous media