Seepage Model and Pressure Response Characteristics of Non-Orthogonal Multi-Fracture Vertical Well with Superimposed Sand Body in Tight Gas Reservoir

1. Introduction

There are many difficulties in the development process of tight sandstone gas reservoirs^[1-4]. The sand body configuration is complex and the discontinuous sand body is common. The current production dynamic analysis method lacks the relevant consideration of discontinuous overlapping sand bodies, as shown in Figure 1. The fracture characteristics after fracturing are more complex, showing non-orthogonal multi-fracture characteristics, which makes its seepage law more complicated ^[5-10]. In order to correctly understand the influence of the above factors on the production and production dynamics of gas wells, it is necessary to establish the corresponding seepage model.

At present, many scholars have studied the seepage model of fractured well in tight reservoir. Zerzar and Bettam[11] combined boundary element method and integral transformation method to provide analytical solutions for multi-stage fractured horizontal Wells in semi-infinite or anisotropic closed reservoirs. Further, Li [12] used point source function, Green function, Fourier and Laplace transform methods to give the formula of unstable pressure for multi-stage fractured horizontal wells with infinite conductivity, and took into account the influence of fracture wall pollution and wellbore storage on the unstable pressure.

Using superposition principle, Amini [13] innovatively used the volumetric source method for the first time to deduce the productivity of horizontal wells with multiple multi-dimensional fractures. Fan et al. [14] used mirror mapping and superposition principle to give the pressure solution of multi-stage fractured horizontal wells with different inclination angles. However, both of these models suffer from slower calculation speeds due to the inability to utilize the Newman product.

To improve the calculation speed, Yao et al. [15] proposed a trilinear seepage model for multi-stage fractured horizontal wells based on the characteristics of start-up fracturing gradient, and believed that a large start-up pressure gradient would lead to an upturn of the well test curve. However, due to the adoption of the trilinear seepage model and the consideration of the influence of start-up pressure gradient, the radial flow in the middle and middle stages of the well test curve was lacking. Thus, to address this issue, Wang et al. [16] proposed the approximate solution of multi-stage fractured horizontal well in rectangular reservoir by using the influence function of conductivity. However, the adoption of the influence function of conductivity led to errors in the early linear seepage stage, especially for low conductivity fractures.

By using the semi-analytical and Laplace transform methods, Meng et al [17] proposed a method for fracture evaluation and well test analysis of parameters of non-uniformly fractured multi-fracture horizontal Wells, which has the advantages of simple and rapid calculation. Liu et al. [18] used Saphir numerical simulation software to study the applicability of trilinear seepage model. They believed that trilinear seepage model was only applicable to multi-stage fractured horizontal wells with dimensionless conductivity greater than 1 of fractures.

Recently, there are some new proposed model for multi-stage fractured horizontal wells. For instance, Zhang et al. [19] proposes a composite model to investigate the pressure behavior and production performance of a multi-wing hydraulically fractured multiple fractured vertical well (MFVW) in a coalbed gas reservoir with a stimulated reservoir volume (SRV). Wang [20] presents a semi-analytical model to simulate the bottom-hole pressure (BHP) performances of a fractured well with multiple radial hydraulic fractures (MRHF) in a stress-sensitive coal seam gas reservoir. What’s more, Zou et al. [21] found that fractal geometry based on fractal theory can better describe reservoir porosity due to its self-similarity and scale invariance. The complexity and disorder of the medium are closer to the reality, so the fractal theory is applied to the well test model of multi-slot with limited diversion of coalbed methane reservoir. Using Fourier transform, perturbation technique, Laplace transform and other methods, Zhao [22] established a horizontal well model of double porosity complex tight gas reservoir with consideration of stress sensitivity effect, and explained the influence of stress sensitivity effect on pressure and pressure derivative curve.

At present, there is a lack of understanding of seepage flow model and law considering superimposed sand bodies and complex fractures comprehensively, necessitating the adoption of an analytical model. ^[23-30]. Therefore, it is necessary to establish corresponding models and methods based on actual gas reservoir status, analyze pressure response characteristics, and clarify the influence law of reservoir type on production and pressure.

2. Establishment and Solution of Seepage Model

Within a tight sandstone gas reservoir, Figure 2 illustrates a schematic diagram of a multi-fractured vertical well. Similarly, in order to facilitate the solution of a multi-fractured vertical well in a tight sandstone gas reservoir considering stress sensitivity and supply boundaries, this section assumes that the entire physical model consists of three zones: the inner zone (the sand body where the gas well is located), the semi-permeable strip zone (the sand body discontinuity zone described by the supply boundary), and the outer zone (the external sand body supply zone). The basic assumptions of the physical model are as follows ^[31-36]:

(1)

The tight sandstone gas reservoir presents a circular distribution on the whole. The initial pressure in the whole area is p_i, the thickness of the tight sandstone gas reservoir is h, and the porosity is φ;

(2)

The whole tight sandstone gas reservoir is homogeneous isotropy, in which the initial permeability of the inner zone is k_g1, the permeability of the semi-permeable zone is k_gm, and the initial permeability of the outer zone is k_g2;

(3)

It is assumed that any hydraulic fracture is distributed in the inner zone and each hydraulic fracture completely opens the reservoir vertically, and that the fracture is distributed in the inner zone laterally without passing through the semi-permeable zone. The half-length of the i fracture is marked as r_fi, and the Angle between the i fracture and the x-axis is θ_i;

(4)

The hydraulic fracture is a finite conductivity fracture, and the flow in each fracture is independent of each other. The tip of the fracture can be regarded as an impermeable boundary, so the flow of tight sandstone gas into the wellbore through the tip of the fracture is negligible;

(5)

Flow patterns throughout the reservoir and in fractures are Darcy's law. Constant production in vertical Wells (q_sc) is maintained.

Since the research object in the multi-fracture vertical well model is the fracture of arbitrary length, the concept of equivalent fracture half-length is used in the dimensionless definition, where the equivalent fracture half-length is:

$${\overline{r}}_f=\frac{{\displaystyle \sum _{i=1}^N{r}_{fi}}}{N}$$

(1)

The basic point source solution of the above physical model is:

$$\frac{{\tilde{\psi }}_{1D0}}{{\tilde{q}}_{Di}}=K_0\left(\sqrt{s}{r}_D\right)+X{I}_0\left(\sqrt{s}{r}_D\right)$$

(2)

where X is the mobility ratio between outer and inner zones.

In order to obtain the solution of any line source, the above equation can be deformed into:

$$\frac{{\tilde{\psi }}_{1D0}}{{\tilde{q}}_{Di}}={\displaystyle \underset{\Gamma }{\int }{K}_0\left(\sqrt{s}{r}_D\right)+X{I}_0\left(\sqrt{s}{r}_D\right)}dl$$

(3)

where dl is a micro source of arbitrary length along the fracture direction.

Therefore, for any line source, the above equation can be rewritten as:

$$\frac{{\tilde{\psi }}_{1D0}}{{\tilde{q}}_{Di}}={\displaystyle \underset{\Gamma }{\int }\left\{\begin{array}{l}{K}_0\left[\sqrt{s}\sqrt{{\left(x_D-x_{pD}\right)}^2+{\left(y_D-y_{pD}\right)}^2}\right]+\\ X{I}_0\left[\sqrt{s}\sqrt{{\left(x_D-x_{pD}\right)}^2+{\left(y_D-y_{pD}\right)}^2}\right]\end{array}\right\}}dl$$

(4)

where (x_D, y_D) represents any point in tight sandstone gas reservoir, and (x_pD, y_pD) represents the position of point source in tight sandstone gas reservoir. For any radial coordinate system, there are:

$$dl=\sqrt{d{x}^2+d{y}^2}$$

(5)

According to the superposition principle, for N cracks (the number of discrete segments in the k crack is n_k), its pseudo-pressure distribution can be written as:

$$\frac{{\tilde{\psi }}_{1D0}\left(i,j\right)}{{\tilde{q}}_{Di}}={\displaystyle \sum _{k=1}^N{\displaystyle \sum _{j=1}^{n_k}\frac{1}{\left|\cos {\theta }_i\right|}{\displaystyle \underset{x_{pD,i}-\frac{L_{fD,i}}{2}\cos {\theta }_i}{\overset{x_{pD,i}+\frac{L_{fD,i}}{2}\cos {\theta }_i}{\int }}\left\{\begin{array}{l}{K}_0\left[\sqrt{s}\sqrt{{\left(x_{D,j}-u\right)}^2+{\left(y_{D,j}-u\tan {\theta }_i\right)}^2}\right]+\\ X{I}_0\left[\sqrt{s}\sqrt{{\left(x_{D,j}-u\right)}^2+{\left(y_{D,j}-u\tan {\theta }_i\right)}^2}\right]\end{array}\right\}du}}}$$

(6)

Obviously, the analytical solution of Equation (6) cannot be determined in the pull space, so it needs to be solved by numerical integration. Considering the calculation speed and accuracy, the Gauss - Legendre numerical integration method is adopted in this paper.

Meanwhile, Equation (6) can be written in the following matrix form:

$${\overrightarrow{\mathrm{A}}}_{\mathrm{MFV}}\left(\begin{array}{c}{\tilde{q}}_{D,1}\\ ⋮\\ {\tilde{q}}_{D,MFV}\end{array}\right)-\left(\begin{array}{c}{\tilde{\psi }}_{1D0，1}\\ ⋮\\ {\tilde{\psi }}_{1D0,MFV}\end{array}\right)=0$$

(7)

In Equation (7),$MFV={\displaystyle \sum _{k=1}^N{\displaystyle \sum _{j=1}^{n_k}j}}$

The traditional fracture model solves the fracture system parallel to the X-axis. When there is an asymmetric distribution of cracks, the classical Cinco-Ley integral equation of finite conductivity is no longer applicable. In order to obtain the solution of the above model, the fracture wing equation proposed in previous studies is adopted in this paper.

$${\tilde{\psi }}_{1fD0}\left(r_D=0\right)-{\tilde{\psi }}_{1fD0}\left(r_D\right)=\frac{2\pi }{C_{fD}}\left[{\tilde{q}}_{fwD}{r}_D-{\displaystyle \underset{0}{\overset{x_D}{\int }}{\displaystyle \underset{0}{\overset{v}{\int }}{\tilde{q}}_{fD}\left(u,s\right)}dudv}\right]$$

(8)

Similarly, Equation (8) can be written in the following matrix form:

$$\left(\begin{array}{c}{\tilde{\psi }}_{1wD0}\\ ⋮\\ {\tilde{\psi }}_{1wD0}\end{array}\right)-\left(\begin{array}{c}{\tilde{\psi }}_{1fD0，1}\\ ⋮\\ {\tilde{\psi }}_{1fD0,SFV}\end{array}\right)={\overrightarrow{{\rm B}}}_{SFV}\left(\begin{array}{c}{\tilde{q}}_{fD,1}\\ ⋮\\ {\tilde{q}}_{fD,SFV}\end{array}\right)$$

(9)

In Equation (9),$SFV={\displaystyle \sum _{j=1}^{n_k}j}$

For a multi-fracture system, the above formula can be written as:

$$\left(\begin{array}{c}{\tilde{\psi }}_{1wD0}\\ ⋮\\ {\tilde{\psi }}_{1wD0}\end{array}\right)-\left(\begin{array}{c}{\tilde{\psi }}_{1fD0，1}\\ ⋮\\ {\tilde{\psi }}_{1fD0,MFV}\end{array}\right)={\overrightarrow{{\rm B}}}_{MFV}\left(\begin{array}{c}{\tilde{q}}_{fD,1}\\ ⋮\\ {\tilde{q}}_{fD,MFV}\end{array}\right)$$

(10)

On the fracture wall, the flow rate and pressure meet the continuity condition, namely:

$$\left(\begin{array}{c}{\tilde{\psi }}_{1D0，1}\\ ⋮\\ {\tilde{\psi }}_{1D0,MFV}\end{array}\right)=\left(\begin{array}{c}{\tilde{\psi }}_{1fD0，1}\\ ⋮\\ {\tilde{\psi }}_{1fD0,MFV}\end{array}\right)$$

(11)

$$\left(\begin{array}{c}{\tilde{q}}_{D,1}\\ ⋮\\ {\tilde{q}}_{D,MFV}\end{array}\right)=\left(\begin{array}{c}{\tilde{q}}_{fD，1}\\ ⋮\\ {\tilde{q}}_{fD,MFV}\end{array}\right)$$

(12)

Combining Equations (9) and (10), and substituting Equations (11) and (12), the following can be obtained:

$$\left({\overrightarrow{A}}_{MFV}+{\overrightarrow{{\rm B}}}_{MFV}\right)\left(\begin{array}{c}{\tilde{q}}_{fD,1}\\ ⋮\\ {\tilde{q}}_{fD,MFV}\end{array}\right)-\left(\begin{array}{c}{\tilde{\psi }}_{1wD0}\\ ⋮\\ {\tilde{\psi }}_{1wD0}\end{array}\right)=0$$

(13)

On the wellbore, the flow satisfies the normalization condition as follows:

$$\left(\begin{array}{ccc}\Delta l_{D1}& \cdots & \Delta l_{D,MFV}\end{array}\right)\left(\begin{array}{c}{\tilde{q}}_{fD,1}\\ ⋮\\ {\tilde{q}}_{fD,MFV}\end{array}\right)=\frac{1}{s}$$

(14)

By combining Equation (13) and Equation (14), the following matrix equation can be obtained:

$$\left(\begin{array}{cc}{\overrightarrow{A}}_{MFV}+{\overrightarrow{{\rm B}}}_{MFV}& -I_{MFV}\\ I_{MFV}^T& 0\end{array}\right)\left(\begin{array}{c}{\overrightarrow{\tilde{q}}}_{fD,MFV}\\ {\tilde{\psi }}_{1\mathrm{wD}0}\end{array}\right)=\left(\begin{array}{c}{0}_{MFV}\\ \frac{1}{\mathrm{s}}\end{array}\right)$$

(15)

In Equation (15):$\begin{array}{l}{\overrightarrow{\tilde{q}}}_{fD,MFV}={\left(\begin{array}{ccc}{\tilde{q}}_{fD,1}& \cdots & {\tilde{q}}_{fD,MFV}\end{array}\right)}^T,\hspace{1em}{I}_{MFV}={\left(\begin{array}{ccc}1& \cdots & 1\end{array}\right)}^T\\ I_{MFV}^T=\left(\begin{array}{ccc}1& \cdots & 1\end{array}\right),\hspace{1em}\hspace{1em}{0}_{MFV}={\left(\begin{array}{ccc}0& \cdots & 0\end{array}\right)}^T\end{array}$

The matrix in Equation (15) is a real matrix, which can be solved by Gaussian elimination method. By solving the above equation, the flow rate and bottom-hole pressure of any discrete line source in the pull space can be obtained. Then, the bottom-hole pressure and the flow rate of discrete fracture unit in real space can be obtained by Stehfest numerical inversion technique.

3. Model Verification and Analysis

For the mathematical model in this paper, there is no similar model to verify the accuracy of the mathematical model in this paper. Therefore, the approach adopted in this paper is to verify the accuracy of the new model by using the approximate model in the classical literature. Luo and Tang presented the unsteady pressure solution of multi-fractured vertical Wells in infinite formations. Therefore, this paper only needs to degenerate the multi-fractured model of tight sandstone gas reservoir with semi-permeable supply boundary to obtain the multi-fractured model of tight sandstone gas reservoir with finite boundary. The comparison between the solution of the degradation model in this paper and that in references is shown in Figure 3.

Figure 3 reveals that, with the exception of the early flow stage, the solution proposed in this paper exhibits relatively larger values; however, the difference between the two remains within an acceptable range. The error may be caused by the following two aspects: first, the number of discrete segments of different cracks; Second, different point source functions (the reference uses infinite point source function, this paper uses closed bounded point source function).

4. Analysis of Pressure Response Curve and Influencing Factors

4.1. Pressure Response Characteristics

In order to more accurately analyze the unsteady seepage characteristics of multi-fractured vertical Wells in tight sandstone gas reservoirs, the pseudo pressure response curve of multi-fractured vertical Wells considering pressure sensitive and semi-permeable supply boundary and corresponding flow section division are given in this paper, as shown in Figure 4. According to Figure 4, the pseudo-pressure derivative response curve of multi-fracture fractured vertical Wells can be roughly divided into the following flow segments: (1) bilinear flow; (2) interseam interference flow; (3) linear stratigraphic flow; (4) Radial flow; (5) Regional channeling; (6) Boundary control control flow. Compared with the pressure response curve of the single-fracture vertical well model in tight sandstone gas reservoir, the multi-fracture model has only one more fracture interference flow in the flow pattern. For the sake of space, the flow pattern of interseam interference flow is simply described in this section.

The interfracture interference flow mainly occurs in the early stage of flow, and its curve morphology is characterized by "convex" groove on the pseudo-pressure derivative curve, which usually occurs before linear flow in formation. In the phase of interfracture interference seepage, the flow on each fracture wall will rush to the dominant flow channel, so that the flow of fluid interferes with each other, resulting in additional pressure loss, resulting in the appearance of "bulge". The size of the "bulge" depends on the conductivity of the fracture and the Angle between adjacent fractures. Generally speaking: the greater the fracture conductivity, the "bulge" becomes smaller or even disappear; The larger the Angle between adjacent cracks is, the smaller the interference between cracks is, resulting in the smaller or even disappeared "bulge".

4.2. Influencing Factors of Pressure Response Curve of Multi-Fracture Vertical Well

(1)

Pressure sensitivity

The influence of pressure sensitivity on pseudo pressure and pseudo pressure derivative response curves of a multi-fractured vertical well in a tight sandstone gas reservoir, considering pressure sensitivity and supply boundary, is depicted in Figure 5.

As can be seen from Figure 5, the influence of pressure sensitivity on pseudo pressure and pseudo pressure derivative response curve of the multi-fracture vertical well model is mainly manifested in the middle and late flow stage, that is, the seepage stage after the pressure wave hits the boundary. With the increase of pressure sensitivity, the propagation speed of pressure waves in the reservoir becomes slower and slower, resulting in greater pressure loss in the reservoir. It is highlighted in the response curve as follows: with the increase of the pressure sensitive coefficient, the larger the pseudo-pressure derivative, the more obvious the upwardness.

(2)

Storage Capacity Ratio Of Inner And Outer Areas

Figure 6 shows the influence of reservoir capacity ratio on the response curve of multi-fractured vertical well model in tight sandstone gas reservoir. It can be seen that the influence of the variation of the storage capacity ratio in the inner and outer zones on the dimensionless pseudo pressure and pseudo pressure derivative curves is mainly manifested in the middle and late seepage stages, that is, the flow section after the middle radial seepage stage. With the increase of the storage capacity ratio in the inner and outer regions, the values of the dimensionless pseudo pressure and pseudo pressure derivative become smaller and smaller. The specific impact is as follows: In the channing phase, with the increase of the internal and external zone capacity ratio, "grooves" gradually appear and gradually become wider and larger. This is mainly because with the increase of the internal and external zone capacity ratio, the reservoir presents the flow characteristics of the two-pore medium reservoir (the flow characteristics of the strong flow ability of the fracture, and the flow characteristics of the poor flow ability of the matrix, Causing its curve to be similar to that of the two-pore medium); In the boundary control flow stage, with the increase of the storage capacity ratio in the inner and outer regions, the later the pressure wave arrives at the boundary, the later the pseudo pressure and pseudo pressure derivative curves coincide on the curve. This is mainly because with the increase of the storage capacity ratio in the outer region, the propagation speed of the pressure wave in the outer region is relatively slow, resulting in the pseudo pressure and pseudo pressure derivative curves coincide later.

(3)

Permeability Coefficient

The influence of inner and outer semi-permeability coefficients on the response curves of multi-fractured vertical Wells in tight sandstone gas reservoirs, considering pressure sensitive and supply boundaries, is revealed in Figure 7 and Figure 8, respectively.

Figure 7 focuses on the influence of the inner semi-permeability coefficient on the pseudo pressure and pseudo pressure derivative curves of the multi-fracture vertical well model. It can be seen from Figure 7 that the influence of the inner semi-permeability coefficient on the pseudo pressure and pseudo pressure derivative curve of the multi-fractured vertical well model of tight sandstone gas reservoir is mainly manifested in the middle and late seepage stage, that is, the flow stage after the middle radial seepage stage, which is similar to the influence of the storage volume ratio curve of the inner and outer zones. It can be seen from Figure 7 that in the middle and late percolation stages, pseudo-pressure and pseudo-pressure derivatives increase with the increase of the inner semi-permeability coefficient, indicating that the increase of the inner semi-permeability coefficient is not conducive to the rapid development of tight sandstone gas reservoirs. The specific effects are as follows: in the phase of channeling flow, with the increase of the inner zone semi-permeability coefficient, the duration of channeling flow is longer and longer. When the inner zone semi-permeability coefficient increases to a certain value, the "groove" begins to appear. At this time, the fluid flow pattern in the reservoir is similar to the fluid flow pattern in the two-hole medium model. In the boundary control flow stage, with the increase of the semi-permeability coefficient in the inner region, the later the boundary control flow appears, the later the pseudo pressure and pseudo pressure derivative curves coincide on the curve.

The influence of the outer semi-permeability coefficient on the curve is also manifested in the middle and late seepage stages, but its influence on the pseudo-pressure and derivative response curve is opposite to that of the inner semi-permeability coefficient on the curve as Figure 8 shown. The specific influence is shown as follows: in the middle and late percolation stage, the pseudo pressure and pseudo pressure derivative gradually decrease with the increase of the outer semi-permeability coefficient, which is a favorable factor for the rapid development of tight sandstone gas reservoirs.

Based on the influence characteristics of the inner and outer zone semi-permeability coefficients on the corresponding curves, it can be seen that the existence of the best inner and outer zone semi-permeability coefficients makes the development of tight sandstone gas reservoirs achieve the best production effect. In addition, during the channeling phase, with the increase of the outer zone semi-permeability coefficient, the "groove" is more obvious, showing a greater depth of "groove" on the curve. In the boundary control flow stage, with the increase of the outer semi-permeability coefficient, the pseudo pressure and pseudo pressure derivative curves coincide later.

(4)

Fracture Symmetry

In the actual fracturing process, the two sides of the hydraulic fracture are usually not strictly symmetrical about the wellbore. Therefore, it is of practical significance to study the effect of fracture symmetry on pseudo pressure response of multi-fracture vertical well in tight sandstone gas reservoir. In this section, the following symmetry factors are defined as follow^[33-35]:

$$\alpha =\frac{\left|L_{fl}-L_{fs}\right|}{\left|L_{fl}+L_{fs}\right|}$$

(16)

when α=0, it means that the two wings of the fracture are strictly symmetric about the wellbore. When α≠0, it means that the two wings of the fracture are asymmetrical with respect to the wellbore. Therefore, the value range of symmetry factor is usually [0,1). The more α value approaches 0, it means that the two wings are approximately symmetric about the wellbore, and the larger α value is, the stronger the asymmetry of the two wings. Figure 9 describes the influence of fracture symmetry factor on the corresponding pseudo-pressure curve of multi-fractured vertical well in tight sandstone gas reservoir.

It can be seen from Figure 9 that the influence of fracture symmetry on the pseudo-pressure response curve of multi-fracture vertical Wells in tight sandstone gas reservoirs is mainly reflected in the early transition stage, that is, the formation linear see-through stage and the excessive flow stage before radial flow. The specific effects are as follows: in the linear seepage stage of the formation, with the increase of fracture asymmetry (symmetry factor), the pseudo-pressure and pseudo-pressure derivatives become larger and larger, mainly because with the increase of asymmetry, there are differences in the utilization of the two wings of the fracture in tight sandstone gas reservoirs. In the ITD before the radial flow, the pseudo-pressure derivative tends to become slightly smaller with the increase of asymmetry, mainly due to the increase of fracture symmetry, which leads to the formation entering the intermediate radial flow earlier (as shown by the black line in Figure 9). According to the above curve characteristics, it can be considered that when the curve of tight sandstone gas reservoir reverses in formation linear flow and intermediate radial flow characteristics, it may be caused by asymmetrical hydraulic fractures. On the whole, in the process of hydraulic fracturing, the hydraulic fractures generated by volume fracturing should be symmetrical about both ends of the wellbore as far as possible to avoid the impact of asymmetric fractures on the initial productivity.

(5)

The Number Of Fractures

In hydraulic fracture operation, the relationship among total fracture length, fracture number and fracture initiation number is complex, and they affect each other, and the specific relationship is affected by rock properties, injection parameters, fracture transmission and other factors. Longer total fracture length may lead to changes in some characteristics of well test analysis curve, such as fluid transfer properties and permeability estimation. The increase of total fracture length may lead to more significant fluid transfer characteristics in the well test analysis curve, thus affecting the interpretation of reservoir properties. There is a relationship between the number of cracks and the pressure interference effect between cracks. Increasing the number of cracks may lead to increased pressure interference between cracks. As the number of cracks increases, fluid transfer and interaction between cracks can lead to changes in productivity and pressure distribution. This may affect the fracture production effect. The number of initiation bars will affect the distribution and morphology of cracks. Different fracture initiation strategies may lead to different fracture distribution, which may affect rock fracture coverage and fracture propagation path, which may affect the fluid transfer and production effect of the reservoir to some extent. In summary, the relationship between total fracture length, fracture number and fracture initiation number is complex, and the influence of these parameters may vary in different reservoirs and operating conditions.

The influence of fracture number on the curve of tight sandstone gas reservoir is mainly manifested in the early seepage stage, namely, bilinear seepage stage, fracture interference flow stage and linear formation seepage stage. Figure 10 reveals the influence of fracture number on pseudo pressure response curve of multi-fracture vertical well in tight sandstone gas reservoir. The specific impact is as follows: In the bilinear seepage stage, with the increase of the number of fractures, the dimensionless pseudo-pressure and pseudo-pressure derivatives become smaller and smaller, but the decreasing trend of pseudo-pressure derivatives gradually weakens when the number of fractures increases to a certain extent. This is mainly because when the number of fractures increases to a certain value, the interference between adjacent fractures is intensified, and the influence of increasing the number of fractures on improving the gas production of a single well is weakened. In the fracture interference flow stage, with the increase of the number of cracks, the "bulge" is more and more obvious and the amplitude of the "bulge" is larger and larger. This is mainly because with the increase of the number of cracks, the discharge area of each crack is gradually smaller. Although the Angle between cracks is consistent, the crack interference is more and more obvious, resulting in the pressure loss caused by the crack interference. The amplitude of "bulge" on the curve is increasing; In the linear seepage stage of formation, with the increase of the number of fractures, the pseudo-pressure and derivative of the pseudo-pressure become smaller and smaller, which shows that the curves of the pseudo-pressure and derivative of the pseudo-pressure get closer and even coincide in this seepage stage. In addition, with the increase of the number of fractures, the duration of the intermediate radial flow tends to be shorter and the arrival time of the radial flow becomes later and later.

In general, for multi-fracture vertical Wells, although the increase of the number of fractures intensifies the mutual interference between fractures to some extent, with the increase of the number of fractures, it is more and more conducive to production. Therefore, in the process of fracturing construction design, more and more fractures should be generated by hydraulic fracturing.

(6)

Fracture Network

Figure 11 mainly describes the influence of different fracture networks on response curves of multi-fracture vertical Wells in tight sandstone gas reservoirs. In this section, the following four fracture networks are mainly considered: one is symmetrical and equal length fracture network; Asymmetry isometric fracture network; Third, symmetrical unequal length fracture network (two long and one short); Fourth, symmetrical unequal length fracture network (two short and one long). In the above four sewing network situations: sewing network 1 and sewing network 2 form a contrast effect, sewing network 2 and sewing network 3 form a contrast effect, sewing network 3 and sewing network 4 form a contrast effect.

It can be seen from Figure 11 that the influence of different fracture network arrangements on the quasi-pressure and quasi-pressure derivative curves of the multi-fracture vertical well model in tight sandstone gas reservoirs is mainly concentrated in the linear see-through stage after the fracture interference flow. The specific effects are as follows: In the early linear flow stage, the pseudo-pressure and pseudo-pressure derivative values corresponding to the symmetric isometric fracture network distribution are the smallest, which indicates that the production is higher under the same pressure drop value, so the production of the symmetric isometric fracture network distribution well is higher in this flow stage. This is mainly because in the early seepage stage, under the condition of symmetrical distribution of fracture network, the fluid around each fracture can quickly and evenly flow into the wellbore through the fracture. In the linear seepage phase before the radial flow, the pseudo pressure and derivative values corresponding to the symmetric unequal length (two short and one long) fracture network are smaller, mainly because the existence of longer fractures leads to larger discharge area corresponding to longer fractures and relatively higher production. However, compared with the early case, the influence of different fracture network arrangement on pseudo pressure and pseudo pressure derivative is weak. Therefore, it is suggested that in the process of hydraulic fracturing, the fracture network should be equally distributed as far as possible, which is more conducive to the rapid production of tight sandstone gas reservoirs.

5. Discussion

Based on the method of point source function and line source function, the unsteady percolation law and gas-water two-phase percolation characteristics of tight sandstone gas reservoir considering the stress sensitive phenomenon and discontinuous reservoir boundary are studied in this paper. This is a preliminary exploration to solve the complex seepage problem of tight sandstone gas reservoir by using semi-analytical method, and the research in many aspects needs to be strengthened. Combined with the research results of this paper and my personal thinking on the complex seepage problem of tight sandstone gas reservoir, it is believed that relevant research can be carried out in the future in the following aspects:

(1)

For tight sandstone gas reservoirs, cross-fracture networks are formed in the fracturing process. However, the fracture network in this paper only considers non-coplanar fractures, and does not study the complex seepage of cross-fracture networks. Therefore, further studies are needed for the unstable seepage of such complex fracture networks.

(2)

In the development process of tight sandstone gas reservoir, the occurrence of gas-water two-phase seepage may lead to the temporary closure and opening of artificial fractures. However, this paper only considers the artificial fractures with constant conductivity and ignores the temporary closure and opening of fractures. Therefore, further studies on the temporary closure and opening of fractures in the gas-water two-phase seepage process of tight sandstone gas reservoir need to be strengthened.

6. Conclusions

(1)

The phenomenon of superimposed sand body deposition in tight gas reservoirs is common, and the well test curve warps up and then falls down. This phenomenon is often analyzed based on the constant pressure boundary or radial composite reservoir model, and the inversion results are prone to be misleading. In this paper, a mathematical model of unsteady seepage flow in fractured vertical Wells in tight sandstone gas reservoirs is established, which takes into account factors such as stress sensitivity, fracture density and fracture symmetry.

(2)

There are great differences in the seepage mechanism at different stages of complex fracture Wells with superposition sand bodies. In the early linear flow stage, interwing interference plays a dominant role, which lasts to the middle and late production stage. In the middle seepage stage, the interfracture interference plays a dominant role until the boundary control flow appears. Finally, the interference between the wings is reflected in the curve as a "bulge" in the pseudo pressure derivative, and the interference between the slots will affect the duration of the first radial flow.

(3)

For tight sandstone gas reservoirs, cross-fracture networks are formed in the fracturing process. However, the fracture network in this paper only considers non-coplanar fractures, and does not study the complex seepage of cross-fracture networks. Therefore, the unstable seepage of such complex fracture networks needs further study.