Summary
Characterization of wormhole growth and propagation dynamics has been made possible with the semi-analytical models developed in this work by integrating both rate transient analysis (RTA) and a pressure-gradient-based (PGB) sand failure criterion in the presence of foamy oil behavior. The nonlinearity caused by the foamy oil properties is linearized using pseudofunctions, while the source function and finite difference methods are applied to obtain the solutions for the linearized models in the matrix and wormhole subsystems, respectively. A sequential method is adopted to solve the coupling fluid-solid problem, where new wormhole segments can be generated once the PGB sand failure criterion has been achieved, while their updated petrophysical properties can be sent back to obtain the solutions for the fluid flow in the next time interval. Furthermore, the influence of sand failure/fluidizing and foamy oil properties on the dynamic wormhole network can be investigated with the generated type curves. Both wormhole growth and foamy oil flow dictate an upward fluid production at the early times, while the production-induced pressure depletion dominates the declining production at the late times. The fractal wormhole networks can be dynamically characterized by history matching the field fluid and sand production profiles. Both the sand failure degree and foamy oil properties dictate the effective wormhole coverage and intensity. Sand production is found to be influenced by both the breakdown pressure gradient and wormhole conductivities, while oil production rate can be dominated by the wormhole coverage and intensity. Foamy oil flow can increase the oil production rate and increase the wormhole coverage compared with the conventional oil. The newly developed method in this work has been validated and then applied at field scale to dynamically characterize the wormhole growth and propagation by considering both the sand failure phenomenon and foamy oil properties within a unified, consistent, and accurate framework. Compared to the conventional numerical simulations, the novelty of such a semi-analytical method can not only allow us to characterize the growth and propagation dynamics of the fractal-like wormholes without the need for grid refinement near the wormhole grids but also take the complex sand failure and foamy oil behavior into account. In addition, it can minimize the impact of grid orientation effects in numerical simulations, enabling the growth and propagation of the wormholes in arbitrary directions rather than orthogonal directions conditioned to the largest pressure gradient.