Abstract
Aimed at the development of shale gas reservoirs with large reservoir thickness and multiple layers, this paper carried out a numerical simulation study on the optimization of three different well types: horizontal well, deviated well, and vertical well. To make the model more in line with the characteristics of shale gas reservoirs, a two-phase gas−water seepage mathematical model of shale gas reservoirs was established, considering the adsorption and desorption of shale gas, Knudsen diffusion effect, and stress sensitivity effect. The embedded discrete fracture model was used to describe hydraulic fracture and natural fracture. Based on Fortran language, a numerical simulator for multi-layer development of shale gas reservoirs was compiled, and the calculation results were compared with the actual production data of Barnett shale gas reservoirs to verify the reliability of the numerical simulator. The spread range of hydraulic fractures in the reservoir with different natural fracture densities is calculated by the simulation to determine well spacing and fracture spacing. The orthogonal experimental design method is then used to optimize the best combination of well spacing and fracture spacing for different well types. The results show that the well productivity of the high-density (0.012 m/m2) natural fractures reservoir > the well productivity of the medium-density (0.006 m/m2) natural fractures reservoir > the well productivity of the low-density (0.001 m/m2) natural fractures reservoir. According to the design of the orthogonal test, it can be seen that the most significant factor affecting the productivity of horizontal wells is the fracture spacing in the Y direction. For deviated wells and vertical wells, the X-direction well spacing has the greatest impact on its productivity.