Spontaneous imbibition (SI) is regarded as an effective method to improve the oil recovery in a tight sandstone reservoir, which leads to a significant change in fracturing design and flowback treatment. However, a longtime shut-in period would aggravate the retention of fracturing fluid, which is in contradiction with high production in the field. It is imperative to understand how SI works during shut-in time, so as to maximize the effect of imbibition in oil recovery enhancement. In this study, a series of experiments were conducted to simulate the status of residual oil saturation so that the inner mechanism of imbibition on oil recovery can be investigated. Low-field nuclear magnetic resonance (LF-NMR) was used to provide direct observation of phase changes in different pore sizes. The experimental results show a positive effect of imbibition on residual oil reduction. This phenomenon further elucidates the observations made during the well shut-in, soaking period, and low flowback efficiency. This study aims to understand the mechanism of SI behavior and help to improve the accuracy of production prediction.

Based on the property discontinuity in the radial direction, this paper develops a new composite model to simulate the productivity of volume fractured wells in naturally fractured reservoirs. The analytical solution of this model is derived in detail and its accuracy is verified by the same model’s numerical solution. Detailed analyses of the traditional transient and cumulative rate are provided for the composite model. Results show that volume fracturing mainly contributes to the early-flow period’s production rate. As interregional diffusivity ratio increases or interregional conductivity ratio decreases, the transient rate at the same wellbore pressure increases. Three characteristic decline stages may be observed on transient rate curves and the shape of traditional rate curves in the middle- and late-flow periods depends on naturally fractured medium and boundary condition. By introducing a new pseudo-steady constant, new Blasingame type curves are also established and their features are more salient than those of traditional rate curves. Five typical flow regimes can be observed on these new type curves. Sensitivity analysis indicates that new Blasingame type curves for varied interregional diffusivity ratio, interregional conductivity ratio, interporosity coefficient and dimensionless reservoir radius, except storativity ratio, will normalize in the late-flow period.

A fractured horizontal well is an effective technology to obtain hydrocarbons from tight reservoirs. In this study, a new semi-analytical model for a horizontal well intercepted by multiple finite-conductivity reorientation fractures was developed in an anisotropic rectangular tight reservoir. Firstly, to establish the flow equation of the reorientation fracture, all reorientation fractures were discretized by combining the nodal analysis technique and the fracture-wing method. Secondly, through coupling the reservoir solution and reorientation fracture solution, a semi-analytical solution for multiple reorientation fractures along a horizontal well was derived in the Laplace domain, and its accuracy was also verified. Thirdly, typical flow regimes were identified on the transient-pressure curves. Finally, dimensionless pressure and pressure derivative curves were obtained to analyze the effect of key parameters on the flow behavior, including fracture angle, permeability anisotropy, fracture conductivity, fracture spacing, fracture number, and fracture configuration. Results show that, for an anisotropic rectangular tight reservoir, horizontal wells should be deployed parallel to the direction of principal permeability and fracture reorientation should be controlled to extend along the direction of minimum permeability. Meanwhile, the optimal fracture number should be considered for economic production and the fracture spacing should be optimized to reduce the flow interferences between reorientation fractures.

Many oil wells in closed reservoirs continue to produce in the pseudo-steady-state flow regime for a long time. The principal objective of this work is to investigate the characteristics of two key pseudo-steady-state parameters—pseudo-steady-state constant (b_Dpss) and pseudo-skin factor (S)—for a well penetrated by a fracture with an azimuth angle (θ) in an anisotropic reservoir. Firstly, a general analytical pressure solution for a finite-conductivity fracture with or without an azimuth angle in an anisotropic rectangular reservoir was developed by using the point-source function and spatial integral method, and two typical cases were employed to verify this solution. Secondly, with the asymptotic analysis method, the expressions of pseudo-steady-state constant and pseudo-skin factor were obtained on the basis of their definitions, and the effects of permeability anisotropy, fracture azimuth angle, fracture conductivity and reservoir shape on them were discussed in detail. Results show that all the b_Dpss-θ and S-θ curves are symmetric around the vertical line, θ = 90° and form a hump or groove shape. The optimized fracture direction in an anisotropic reservoir is perpendicular to the principal permeability axis. Furthermore, a new formula to calculate pseudo-skin factor was successfully proposed based on these two parameters’ relationship. Finally, as an application of pseudo-steady-state constant, a set of Blasingame format rate decline curves for the proposed model were established.

Refracturing treatment is an economical way to improve the productivity of poorly or damaged fractured horizontal wells in tight reservoirs. Fracture reorientation and fracture face damage may occur during refracturing treatment. At present, there is still no report on the rate decline solution for refractured horizontal wells in tight reservoirs. In this work, by taking a semi-analytical method, traditional rate decline and Blasingame-type rate decline solutions were derived for a refractured horizontal well intercepted by multiple reorientation fractures with fracture face damage in an anisotropic tight reservoir. The accuracy and reliability of the traditional rate decline solution were verified and validated by comparing it with a classic case in the literature and a numerical simulation case. The effects of fracture reorientation and fracture face damage on the rate decline were investigated in depth. These investigations demonstrate that fracture face damage is not conducive to increasing well productivity during the early flow period and there is an optimal matching relationship between the principal fracture section angle and permeability anisotropy, particularly for the reservoirs with strong permeability anisotropy. The fracture length ratio and fracture spacing have a weak effect on the production rate and cumulative production while the fracture number shows a strong influence on the rate decline. Furthermore, multifactor sensitivity analysis indicates that fracture conductivity has a more sensitive effect on well productivity than fracture face damage, implying the importance of improving fracture conductivity. Finally, a series of Blasingame-type rate decline curves were presented, and type curve fitting and parameter estimations for a field case were conducted too. This work deepens our understanding of the production performance of refractured horizontal wells, which helps to identify reorientation fracture properties and evaluate post-fracturing performance.