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.

In recent years, CO₂ flooding has become an important technical measure for oil and gas field enterprises to further improve oil and gas recovery and achieve the goal of “dual carbon”. It is also one of the concrete application forms of CCUS. Numerical simulation based on CO₂-EOR plays an indispensable role in the study of the mechanism of CO₂ flooding and buried percolation, allowing for technical indicators to be selected and EOR/EGR prediction to be improved for reservoir engineers. This paper discusses the numerical simulation techniques related to CO₂ flooding and storage, including mathematical models and solving algorithms. A multiphase and multicomponent mathematical model is developed to describe the flow mechanism of hydrocarbon components–CO₂–water underground and to simulate the phase diagram of the components. The two-phase P-T flash calculation with SSI (+DEM) and the Newton method is adopted to obtain the gas–liquid phase equilibrium parameters. The extreme value judgment of the TPD function is used to form the phase stability test and miscibility identification model. A tailor-made multistage preconditioner is built to solve the linear equation of the strong-coupled, multiphase, multicomponent reservoir simulation, which includes the variables of pressure, saturation, and composition. The multistage preconditioner improves the computational efficiency significantly. A numerical simulation of CO₂ injection in a carbonate reservoir in the Middle East shows that it is effective for researching the recovery factor and storage quantity of CO₂ flooding based on the above numerical simulation techniques.

Predictive analysis of the reservoir surveillance data is crucial for the high-efficiency management of oil and gas reservoirs. Here we introduce a new approach to reservoir surveillance that uses the machine learning tree boosting method to forecast production data. In this method, the prediction target is the decline rate of oil production at a given time for one well in the low-permeability carbonate reservoir. The input data to train the model includes reservoir production data (e.g., oil rate, water cut, gas oil ratio (GOR)) and reservoir operation data (e.g., history of choke size and shut-down activity) of 91 producers in this reservoir for the last 20 years. The tree boosting algorithm aims to quantitatively uncover the complicated hidden patterns between the target prediction parameter and other monitored data of a high variety, through state-of-the-art automatic classification and multiple linear regression algorithms. We also introduce a segmentation technique that divides the multivariate time-series production and operation data into a sequence of discrete segments. This feature extraction technique can transfer key features, based on expert knowledge derived from the in-reservoir surveillance, into a data form that is suitable for the machine learning algorithm. Compared with traditional methods, the approach proposed in this article can handle surveillance data in a multivariate time-series form with different strengths of internal correlation. It also provides capabilities for data obtained in multiple wells, measured from multiple sources, as well as of multiple attributes. Our application results indicate that this approach is quite promising in capturing the complicated patterns between the target variable and several other explanatory variables, and thus in predicting the daily oil production rate.

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.