A Transient Productivity Model of Fractured Wells in Shale Reservoirs Based on the Succession Pseudo-Steady State Method

After volume fracturing, shale reservoirs can be divided into nonlinear seepage areas controlled by micro- or nanoporous media and Darcy seepage areas controlled by complex fracture networks. In this paper, firstly, on the basis of calculating complex fracture network permeability in a stimulated zone, the steady-state productivity model is established by comprehensively considering the multi-scale flowing states, shale gas desorption and diffusion after shale fracturing coupling flows in matrix and stimulated region. Then, according to the principle of material balance, a transient productivity calculation model is established with the succession pseudo-steady state (SPSS) method, which considers the unstable propagation of pressure waves, and the factors affecting the transient productivity of fractured wells in shale gas areas are analyzed. The numerical model simulation results verify the reliability of the transient productivity model. The results show that: (1) the productivity prediction model based on the SPSS method provides a theoretical basis for the transient productivity calculation of shale fractured horizontal well, and it has the characteristics of simple solution process, fast computation speed and good agreement with numerical simulation results; (2) the pressure wave propagates from the bottom of the well to the outer boundary of the volume fracturing zone, and then propagates from the outer boundary of the fracturing zone to the reservoir boundary; (3) with the increase of fracturing zone radius, the initial average aperture of fractures, maximum fracture length, the productivity of shale gas increases, and the increase rate gradually decreases. When the fracturing zone radius is 150 m, the daily output is approximately twice as much as that of 75 m. If the initial average aperture of fractures is 50 μm, the daily output is about half of that when the initial average aperture is 100 μm. When the maximum fracture length increases from 50 m to 100 m, the daily output only increases about by 25%. (4) When the Langmuir volume is relatively large, the daily outputs of different Langmuir volumes are almost identical, and the effect of Langmuir volume on the desorption output can almost be ignored