Transient Pressure Analysis of Inclined Well in Continuous Triple-Porosity Reservoirs With Dual-Permeability Behavior

Abstract

Inclined wells has recently been adopted to develop fractured-vuggy carbonate hydrocarbon reservoirs (FVCHRs) composed by matrix, fracture and vug system. Therefore, it is significant for us to describe pressure transient of inclined well for FVCHRs. In this paper, it is assumed that vug and fracture system connect with wellbore, and inter-porosity flow from vug to fracture system appears. Therefore, the pressure transient responses model of inclined well with triple-porosity dual-permeability behavior was built. The model is solved by employing Laplace integral transform and finite cosine Fourier transform. Real-domain solution of the model is obtained by Stehfest inversion algorithm. On the basis model of the published paper, the solution of the simplified model of this paper was validated with horizontal and vertical well of FVCHRs with triple-porosity dual-permeability behavior, and results reach a good agreement. Type curve, according to pressure derivative curve characteristic, can be divided into eight flow regimes, which includes wellbore storage, skin reflect, early vertical radial flow, top and bottom boundary reflection, linear flow, inter-porosity flow from vug to fracture, inter-porosity from matrix to vug and fracture and pseudo-radial flow regimes. The influence of some vital parameters (inclination angle, inter-porosity flow coefficient, well length and permeability ratio etc.) on dimensionless pressure and its derivative curves are discussed in detail. The presented model can be used to understand pressure transient response characteristic of inclined wells in FVCHRs with dual-permeability behavior.

Introduction

Since the 1960s, numerous scholars mainly studied the flow characteristic of the fluid in naturally fractured reservoirs by the method of the pressure transient analysis (Warren and Root, 1963; Pruess and Narasimhan, 1985; Kuchuk and Biryukov, 2014; Kui et al., 2017; Tan et al., 2018; Yin, 2018; Song et al., 2019; Pang et al., 2020; Tontiwachwuthikul et al., 2020; Wang et al., 2020). Since a large number of the fractured-vuggy carbonate hydrocarbon reservoirs (FVCHRs) were found around the world, most of scholars have focused on studies of FVCHRs to understand flow characteristic of the fluid in the FVCHRs (Abadassah and Ershaghi, 1986; Liu et al., 2003; Sun et al., 2018). Compared with naturally fractured reservoirs, the complex structures of the FVCHRs posed significant research challenges towards scholars and scientists about pressure transient response of the FVCHRs. Therefore, the mathematic model descripting flow process of the fluid in FVCHRs is established to evaluate the properties and pressure transient characteristics of this type reservoirs.

During the last 20 years, pressure transient analysis is used widely to obtain some vital parameters of reservoirs. With the help of these parameters, reservoir engineers can evaluate oil and gas reservoirs effectively and develop a reasonable development plan for the FVCHRs. The source function method plays a vital role in the process of well test model solution. Gringarten and Ramey (1973) first put forward the real-domain point source analytical solution under the conditions of the different outer boundary. On the basis of the real-domain point source solution, Ozkan and Raghavan (1991; 2004) applied Green’s function theory to expand the point source to the Laplace domain and the corresponding calculation method is given.

Most of the world’s giant oil ans gas produce from naturally fractured and vuggy carbonate reserviors with complex pore system. Vug porosity is common pore system in FVCHRs, and the size of the vug have an imoprtant influence on pressure and production performance. The FVCHRs can be subdivided into triple-continuum medium model and discontinuous medium according to the size of the vug. For the research of the pressure and production of triple-continuum medium model, Abadassah and Ershaghi (1986) first presented triple-porosity single-permeability well-test mathematical model for vertical well, and semi-log and log-log curves were drawn to analysis flow characteristic of the triple-continuum medium. Camacho-Velázquez et al. (2005) established and solved mathematical model of triple-porosity dual-permeability of the vertical well, which assumed that fluid of the natural fracture and vuggy porosity flow into wellbore directly, and the characteristic of pressure transient and decline curve are analyzed. Wu et al. (2004; 2011) solved a triple-continuum pressure transient model of the multi-phase fluid by employing analytical and numerical method. Corbett et al. (2010) studied on fractured carbonate rocks by adopting the numerical well test, and summarized that numerical well testing has its limitations. Guo et al. (2012) established a triple-porosity single-permeability flow mathematical model of the horizontal well in a triple-continuum FVCHRs under different boundary condition, and standard pressure and its derivative log–log curves were drawn. Nie et al. (2012) presented pressure response mathematical model of the horizontal well in coal seam. Pressure and its derivative response curves are drawn in log-log plot, and influence of some vital parameters on type curves were analyzed. Jia et al. (2013) studied on flow issues of a porous–vuggy carbonate reservoir without fracture system. It is assumed that vug system is spherical and pressure distribute symmetrically. Corresponding mathematical model was established and solved by employing Laplace transform. Wang et al. (2018) developed a pressure transient mathematical model of acid fracturing wells in FVCHRs. Therefore, a composite system is considered, the inner region is made up with finite number of artificial fractures and wormholes, and the outer region is made up triple-porosity medium. Wang et al. (2014) derived well test mathematical model of hydraulic fractured well in FVCHRs by point sink integral method. This model assumed that fluid of matrix and natural fractures system flow into wormholes directly. The existence of large-scale cavities will have a noticeable influence on the transient pressure behavior. For the research of the pressure and production of discontinuous medium model, Gao et al. (2016) and Xing et al. (2018) presented a flow mathematical model, it is assumed that well is drilled through large-scale cavities and fluid flow into a filled cavity. Wellbore storage and skin factor were considered in their model. Du et al. (2019) presented an analytical well test model of fractured vuggy carbonate reservoirs. Pressure difference between the wellbore and large vug was considered in their model.

Because the reservoir geological condition becomes more and more complex, the technology of inclined wells has recently been adopted to solve the problem of onshore oil and gas reservoirs extraction. To estimate pressure transient behavior of inclined well, many scholars established pressure response mathematical model of inclined well. Cinco-Ley et al. (1975; 2007) presented an analytical solution of inclined well in real domain and shown the existence of the pseudo-radial flow period at large time values. Yildiz and Ozkan (Yildiz and Ozkan, 1999; Ozkan et al., 1999; Ozkan and Raghavan, 2000) studied transient pressure behavior and rate distributions of slanted well and perforated wells. Corresponding calculation method of slanted and perforated wells was shown to improve the computational efficiency and increase the accuracy. Feng and Liu (2015) established a well test mathematical model of inclined well with an impermeable fault by employing point source theory. Dong et al. (2018) presented a skin factor model of partially penetrated directionally-drilled wells in anisotropic reservoirs to study the stimulation effect of slanted.

This paper summarizes the published literature about well test model of the FVCHRs and slanted well, absorbs their highlights and studies on pressure transient behavior of the inclined well in FVCHRs. Therefore, the objective of this paper is to establish and solve well test mathematical model of inclined well in FVCHRs with triple-porosity dual-permeability behavior by employing Laplace transform and finite cosine Fourier transform. The fluid of natural fractures and vug system can flow into wellbore directly. Inter-porosity flows from matrix system to both natural fracture and vug system, and inter-porosity flow from vug system to a natural fracture system were also considered in this paper. Pressure and its derivative type curves under condition of constant production by adopting the Stehfest numerical inversion are drawn. Influence of some vital parameters on type curves are discussed in detail.

Physical Modeling and Assumptions

FIGURE 1

FIGURE 2

FIGURE 3

Mathematical Modeling

Based on the scheme of an inclined well as shown in Figure 2, the mathematical model of the inclined well is established in FVCHRs. In order to help reader to understand coordinate system as shown in Figure 2, three different views of an inclined well in top and bottom closed boundary are given in Figure 4.

FIGURE 4

For convenience of derivation, the relevant dimensionless variables are defined in Table 1: Laplace transformation with respect to is then adopted to invert the above seepage model into Laplace domain (Van Everdingen and Hurst, 1949).where: is dimensionless pressure in Laplace domain;

According to the dimensionless definition, the dimensionless forms of the differential equation and boundary conditions can be given by:where: is permeability ratio between vug and fracture, dimensionless; is permeability ratio, dimensionless; is inter-porosity flow coefficient from matrix to fracture, dimensionless; is inter-porosity flow coefficient from matrix to vuggy, dimensionless; is inter-porosity flow coefficient from vuggy to fracture, dimensionless; is fracture storativity ratio, dimensionless; is matrix storativity ratio, dimensionless; is vuggy storativity ratio, dimensionless; is Laplace variables; subscripts is dimensionless variables; superscripts – is variables in the Laplace domain;

By employing the finite cosine Fourier transform with respect to in Eq. 10, and expressing dimensionless pressures of natural fracture and vug systems by dimensionless matrix system pressure in Laplace domain after the finite cosine Fourier transform, Eq. 10 can be described as:where: ; ; subscripted variable.

By solving Eq. 10, the general point source solution for the top and bottom closed boundary, constant production and infinite lateral boundary FVCHRs with triple-porosity dual-permeability behavior can be obtained through Laplace transform and Fourier transform methods (Liu et al., 2019; Wang et al., 2019; Xue et al., 2019):where: , ;

It is useful to express the Eq. 12 in polar coordinates, let the polar coordinates of the points and be given by and respectively according to Raghavan and Ozkan (2017). We can do integration to the point source solution along inclined well in polar coordinates system as shown in Figure 4C. Therefore, the solution of the inclined well is obtain as following.where: is the second kind modified Bessel function; ; .

Although Gringarten et al. (1974) pointed out that the wellbore pressure of horizontal well and infinite conductivity hydraulic fractured well can be calculate by an equivalent-pressure point , there is difference between horizontal well and inclined well. Fortunately, Cinco-Ley et al. (1975) given the equivalent-pressure point of inclined well in the polar coordinates. As following:

In order to enhance calculation speed of the inclined well, Ozkan and Raghavan (2000) provides an efficient algorithm to compute transient pressure responses of inclined wells in Laplace domain. The principle of convolution can be used to obtain solution with the effects of wellbore storage and skin factor, which could be given bywhere: is dimensionless wellbore pressure with wellbore storage and skin; is dimensionless wellbore pressure without wellbore storage and skin; is dimensionless wellbore storage; is skin factor.

The solution of inclined well in real domain for FVCHRs with triple-porosity dual-permeability behavior is obtained by employing Stehfest (1970) numerical inversion algorithm. Dimensionless pressure and its derivative type curves are drawn analyzed.

Model Verification

Compared with vertical wells and horizontal wells, wellbore pressure solution of the inclined well is more universal, so our work expands analytical solution of vertical well and horizontal well in FVCHRs with triple-porosity dual-permeability behavior. However, in order to verify the correctness of our model and solution, the simplified model of this paper is compared with analytical solution of horizontal well for FVCHRs with triple-porosity dual-permeability behavior (Guo et al., 2012) and vertical well for FVCHRs with triple-porosity dual-permeability behavior (Camacho-Velázquez et al., 2005).

Guo et al. (2012) established the mathematical model of the horizontal well under condition of different outer boundary in FVCHRs with triple-porosity dual-permeability behavior. What is different with our model is that they assumed that fluid of the naturally fractured and the matrix system are directly flow into the wellbore, inter-porosity flow from matrix to fracture and inter-porosity flow from vug to matrix and fracture is considered. In order to make verification result reliable, if parameters of the Guo’ model are replaced by parameters of this paper correspondingly and corresponding value of these parameters is same. On the basis of replacement, we can prove that our model is same with Guo’ model. If is 90° in our paper, and of our paper, inclined well of this paper will become horizontal well, which is similar to the Guo et al.’ model with horizontal well length is 500. Other common parameters are same and their values are: , , , , , , , , , . As shown in Figure 5A, result reaches a good agreement.

FIGURE 5

Camacho-Velázquez et al. (2005) established the mathematical model of the vertical well in FVCHRs with triple-porosity dual-permeability behavior. It is assumed that fluid of the fracture and the vug system are directly flow into the wellbore. Inter-porosity flow from vug to fracture and inter-porosity flow from matrix to vug and fracture is considered, which is same with our assumption. If inclination angle of wellbore is less than 5° and well is penetrated fully, we can think that our model is equal to vertical well. Other common parameters are same and their values are: , , . As shown in Figure 5A, result reaches a good agreement.

By comparing with horizontal and vertical well in FVCHRs with triple-porosity dual-permeability behavior, the results reached a good agreement, which verifies the correctness of our model and solution.

Type Curve and Sensitivity of Parameters

Type Curve

Figure 6 shows type curves of dimensionless pressure and its derivative of inclined well with dual-permeability behavior. According to the characteristic of the dimensionless pressure derivative curve, the type curve of inclined well can be divided into eight flow regimes for FVCHRs with triple-porosity dual-permeability behavior, which mainly includes well storage regime, skin effect regime, early vertical radial flow regime (Figure 7A), top and bottom boundary reflection (Figure 7B), linear flow regime (Figure 7C), inter-porosity flow and pseudo-radial flow regime (Figure 7D).

FIGURE 6

FIGURE 7

Regime①: Wellbore storage regime. During the pure wellbore storage, pressure and its derivative curves coincide and is a unit-slope straight line.

Regime②: Skin effect regime. Due to consideration of skin, pressure derivative curve looks like a “hump.”

Regime③: Early vertical radial flow regime. Pressure derivative curve is a constant value horizontal line in log-log plot (Figure 7A).

Regime④: Top and bottom boundary reflection regime. Pressure derivative curve take on a small step in log-log plot. However, this stage may disappear because of change of the length of wellbore and inclination angle of wellbore (Figure 7B).

Regime⑤: Linear flow regime. Pressure derivative curve is 1/2-slope straight line in log-log plot. Shorter well length lead to disappearance of 1/2-slope straight line (Figure 7C).

Regime⑥: Inter-porosity flow from vug to fracture. Because vug permeability is better than matrix permeability, inter-porosity flow from vug to fracture appear firstly. The curve of pressure derivative is V-shaped.

Regime⑦: Inter-porosity flow from matrix to vug and fracture. The curve of pressure derivative is V-shaped, which is controlled by parameters ω_m and ω_v.

Regime⑧: Pseudo-radial flow regime. Inter-porosity flows have already finished and the curve of pressure derivative is a 0.5-value horizontal line in log-log plot (Figure 7D).

Sensitivity of Parameters

Figure 8 shows type curve affected by length of inclined well in FVCHRs with triple-porosity dual-permeability behavior. Because short inclined well length can decrease contract area between reservoirs and wellbore, the short inclined well length can increase pressure loss that fluid flows into wellbore. Therefore, short inclined well length can lead to high position of pressure curve in the whole flow regime in log-log plot and make top and bottom boundary reflection regime disappear and spherical flow characteristic more obvious in pressure derivative curve.

FIGURE 8

Figure 9 shows type curve affected by inclination angle in FVCHRs with triple-porosity dual-permeability behavior. It is assumed that length of inclined well keep constant when inclination angle changes. The large inclination angle change can only lead to disappearance of the top and bottom boundary reflection regime of pressure derivative curve without consideration of anisotropy (Figure 9A). However, if the inclined well is penetrated fully, inclination angle has a more effect on pressure and its derivative curves than inclined well penetrated partly. With the increasing of inclination angle length of inclined well become length, which leads to low position of pressure curve in whole regime and low position of pressure derivative curve during early vertical radial flow regime (Figure 9B).

FIGURE 9

Figure 10 shows type curve affected by vug storativity ratio in FVCHRs with triple-porosity dual-permeability behavior. According the flow model shown Figure 3A, fluid of the vug system do not flow into fracture, but also flow into wellbore directly. Therefore, as shown in Figure 10, small vug storativity ratio can make the second V-shaped of the pressure derivative curve become more and more deep while the first V-shaped become more and more shallow.

FIGURE 10

Figure 11 shows type curve affected by fracture storativity ratio in FVCHRs with triple-porosity dual-permeability behavior. According the flow model shown Figure 3A, fluid of fracture system only flow into wellbore directly. Therefore, as is shown in Figure 11, fracture storativity ratio has an influence on the first V-shaped. The first V-shaped become more and more deeply with the decrease of fracture storativity ratio. The position of pressure curve become higher during vertical radial flow regime and top and bottom boundary reflection regime.

FIGURE 11

Figure 12 shows type curve affected by inter-porosity flow coefficient in FVCHRs with triple-porosity dual-permeability behavior. According the flow model shown Figure 3A, fluid of fracture system only flow into wellbore directly. Therefore, as shown in Figure 12, inter-porosity flow coefficient has an influence on the first V-shaped. As increase, starting time of inter-porosity flow from vug system to fracture system become earlier, which makes the first V-shaped of pressure derivative move towards left in log-log plot.

FIGURE 12

Figure 13 shows type curve affected by inter-porosity flow coefficient in FVCHRs with triple-porosity dual-permeability behavior. As shown in Figure 13, inter-porosity flow coefficient λ_mv has an influence on the second V-shaped. As increase, starting time of inter-porosity flow from matrix system to vug system become earlier, which makes the second V-shaped of pressure derivative move towards left in log-log plot. Of course, with the flow coefficient keep increasing, first V-shaped may be covered and first V-shaped become more deeply.

FIGURE 13

Figure 14 shows type curve affected by inter-porosity flow coefficient in FVCHRs with triple-porosity dual-permeability behavior. As shown in Figure 14, inter-porosity flow coefficient has an influence on the second V-shaped. As increase, starting time of inter-porosity flow from matrix system to fracture system become earlier, which makes the second V-shaped of pressure derivative move towards left in log-log plot. Of course, because some of fluid of fracture system come from supplement of vug system, but also others come from supplement of matrix system. Therefore, when inter-porosity flow coefficient , the change of inter-porosity flow coefficient has an influence on pressure derivative curve slightly.

FIGURE 14

Figure 15 shows type curve affected by permeability ratio in FVCHRs with triple-porosity dual-permeability behavior. The large permeability ratio indicates that permeability of fracture system is higher than vug system, and vug storativity capacity is lower than fracture storativity capacity. Therefore, the first V-shaped of pressure derivative curve becomes more and more deeply with the increase of permeability ratio.

FIGURE 15

Figure 16 shows type curve affected by permeability ratio between horizontal permeability and vertical permeability in FVCHRs with triple-porosity dual-permeability behavior. The large indicates that permeability of horizontal permeability is higher than vertical permeability. Because inclination angle of wellbore is considered in vertical direction, difference between vertical and horizontal flow is shown obviously during early vertical radial flow regime. The greater the value of , the greater the loss of wellbore flow pressure in the vertical direction. Therefore, large will lead to higher position of pressure and its derivative curve in early flow regime.

FIGURE 16

Conclusion

Funding

This article was supported by the National Natural Science Foundation of China (51704246) and Sichuan Science and Technology Project (18SYXHZ0090).

References

• 1

  AbadassahD.ErshaghiI. (1986). Triple-porosity systems for representing naturally fractured reservoirs. SPE Form. Eval.1 (2), 113–127. 10.2118/13409-pa

  • CrossRef
  • Google Scholar

• 2

  Camacho-VelázquezR.Vásquez-CruzMCastrejón-AivarR.Arana-OrtizV. (2005). Pressure-transient and decline-curve behavior in naturally fractured vuggy carbonate reservoirs. SPE Reservoir Eval. Eng.8 (2), 95–111. 10.2118/77689-pa

  • CrossRef
  • Google Scholar

• 3

  Cinco-LeyH.MillerF. G.RameyH. J. (1975). Unsteady-State pressure distribution created by a directionally drilled well. J. Petrol. Technol.27, 1392–1400. 10.2118/5131-PA

  • CrossRef
  • Google Scholar

• 4

  Cinco-LeyH.RameyH. J.MillerF. (2007). “Pseudo-skin factors for partially-penetrating directionally-drilled wells,” in Fall meeting of the society of petroleum engineers of AIME, Dallas, TX, 28 September–1 October. 10.2523/5589-ms

  • CrossRef
  • Google Scholar

• 5

  CorbettP. W.M.GeigerS.BorgesL.GarayevM.GonzalezJ. G.ValdezC. (2010). “Limitations in numerical well rest modelling of fractured carbonate rocks.” in 72nd European association of geoscientists and engineers conference and exhibition 2010: a new Spring for geoscience, Barcelona, Spain, 14–17 June, 1135–1145. 10.2523/130252-ms

  • CrossRef
  • Google Scholar

• 6

  DongW.WangX.WangJ. (2018). A new skin factor model for partially penetrated directionally-drilled wells in anisotropic reservoirs. J. Petrol. Sci. Eng.161, 334–348. 10.1016/j.petrol.2017.11.062

  • CrossRef
  • Google Scholar

• 7

  DuX.LuZLiDXuYLiPLuD (2019). A novel analytical well test model for fractured vuggy carbonate reservoirs considering the coupling between oil flow and wave propagation. J. Petrol. Sci. Eng.173, 447–461. 10.1016/j.petrol.2018.09.077

  • CrossRef
  • Google Scholar

• 8

  FengG. Q.LiuQ.-G. (2015). Pressure transient behavior of a slanted well with an impermeable fault. J. Hydrodyn.26 (6), 980–985. 10.1016/S1001-6058(14)60108-6

  • CrossRef
  • Google Scholar

• 9

  GaoB.HuangX.YaoJ.LvX.WuY. (2016). Pressure transient analysis of a well penetrating a filled cavity in naturally fractured carbonate reservoirs. J. Petrol. Sci. Eng.145, 392–403. 10.1016/j.petrol.2016.05.037

  • CrossRef
  • Google Scholar

• 10

  GringartenA. C.RameyH. J. (1973). Use of source and green’S functions in solving unsteady-flow problems in reservoirs. Soc. Pet. Eng. AIME J.13 (5), 285–296. 10.2118/3818-PA

  • CrossRef
  • Google Scholar

• 11

  GringartenA. C.RameyH. J.RaghavanR. (1974). Unsteady-state pressure distributions created by a well with a single infinite-conductivity vertical fracture. Soc. Pet. Eng. AIME J.14 (4), 347–360. 10.2118/4051-pa

  • CrossRef
  • Google Scholar

• 12

  GuoJ. C.NieR. S.JiaY. L. (2012). Dual permeability flow behavior for modeling horizontal well production in fractured-vuggy carbonate reservoirs. J. Hydrol.464–465, 281–293. 10.1016/j.jhydrol.2012.07.021

  • CrossRef
  • Google Scholar

• 13

  JiaY. L.FanX. Y.NieR. S.HuangQ. H.JiaY. L. (2013). Flow modeling of well test analysis for porous-vuggy carbonate reservoirs. Transp. Porous Media97 (2), 253–279. 10.1007/s11242-012-0121-y

  • CrossRef
  • Google Scholar

• 14

  KuchukF.BiryukovD. (2014). Pressure-transient behavior of continuously and discretely fractured reservoirs. SPE Reservoir Eval. Eng.17 (1), 82–97. 10.2118/158096-PA

  • CrossRef
  • Google Scholar

• 15

  KuiM. Q.TanX. H.LiX. P.CaiJ. (2017). Analysis on unsteady flow for power-law fluids in dual fractal media. J. Porous Media20 (12), 1071–1086. 10.1615/JPorMedia.v20.i12.20

  • CrossRef
  • Google Scholar

• 16

  LiuJ.BodvarssonG. S.WuY. S. (2003). Analysis of flow behavior in fractured lithophysal reservoirs. J. Contam. Hydrol.62–63, 189–211. 10.1016/S0169-7722(02)00169-9

  • CrossRef
  • Google Scholar

• 17

  LiuQ.XuY.PengX.LiuY.QiS. (2019). Pressure transient analysis for multi-wing fractured wells in dual-permeability hydrocarbon reservoirs. J. Petrol. Sci. Eng.180, 278–288. 10.1016/j.petrol.2019.05.013

  • CrossRef
  • Google Scholar

• 18

  NieR. S.MengY. F.GuoJ. C.JiaY. L. (2012). Modeling transient flow behavior of a horizontal well in a coal seam. Int. J. Coal Geol.92, 54–68. 10.1016/j.coal.2011.12.005

  • CrossRef
  • Google Scholar

• 19

  OzkanE.RaghavanR. (1991). New solutions for well-test-analysis problems. Part 2. Computational considerations and applications. SPE Form. Eval.6 (3), 369–378. 10.2118/18616-PA

  • CrossRef
  • Google Scholar

• 20

  OzkanE.RaghavanR. (2000). A computationally efficient, transient-pressure solution for inclined wells. SPE Reservoir Eval. Eng.3 (5), 414–425. 10.2118/66206-PA

  • CrossRef
  • Google Scholar

• 21

  OzkanE.RaghavanR. (2004). New solutions for well-test-analysis problems: Part 1-Analytical Considerations (includes associated papers 28666 and 29213). SPE Form. Evaluation.6 (03), 359–368. 10.2118/18615-PA

  • CrossRef
  • Google Scholar

• 22

  OzkanE.YildizT.RaghavanR. (1999). “Pressure-transient analysis of perforated slant and horizontal wells,” in SPE Annual Technical Conference and Exhibition, Houston, TX, October 3–6, 1999.

  • Google Scholar

• 23

  PangY.HuXWangSChenS.SolimanM. Y.DengH. (2020). Characterization of adsorption isotherm and density profile in cylindrical nanopores: modeling and measurement. Chem. Eng. J.396, 125212. 10.1016/j.cej.2020.125212

  • CrossRef
  • Google Scholar

• 24

  PruessK.NarasimhanT. N. (1985). Practical method for modeling fluid and heat flow in fractured porous media. Soc. Petrol. Eng. J.25 (1), 14–26. 10.2118/10509-pa

  • CrossRef
  • Google Scholar

• 25

  RaghavanROzkanE. (2017). A method for computing unsteady flows in porous media. Oxon, UK: Routledge. 10.1201/9781315136387

  • CrossRef
  • Google Scholar

• 26

  SongF.BoL.ZhangS.SunY. (2019). Nonlinear flow in low permeability reservoirs: modelling and experimental verification. Adv. Geo-Energy Res.3, 76–81. 10.26804/ager.2019.01.06

  • CrossRef
  • Google Scholar

• 27

  StehfestH. (1970). Algorithm 368: numerical inversion of Laplace transforms [D5]. Commun. ACM13 (1), 47–49.

  • Google Scholar

• 28

  SunQ.ZhangN.MohamedF.WangY. (2018). Structural regeneration of fracture-vug network in naturally fractured vuggy reservoirs. J. Petrol. Sci. Eng.165, 28–41. 10.1016/j.petrol.2017.11.030

  • CrossRef
  • Google Scholar

• 29

  TanX. H.JiangL.LiX. P.ZhangB. J.LiX. C. (2018). Flow model of a multi-stage hydraulic fractured horizontal well based on tree-shaped fractal fracture networks. J. Petrol. Sci. Eng.169, 494–503. 10.1016/j.petrol.2018.06.008

  • CrossRef
  • Google Scholar

• 30

  TontiwachwuthikulP.ChanC. W.ZengF.LiangZ.SemaT.ChaoM. (Forthcoming2020). Recent progress and new developments of applications of artificial intelligence (AI), knowledge-based systems (KBS), and machine learning (ML) in the petroleum industry. Petroleum. 10.1016/j.petlm.2020.08.001

  • CrossRef
  • Google Scholar

• 31

  Van EverdingenA. F.HurstW. (1949). The application of the Laplace transformation to flow problems in reservoirs. J. Petrol. Technol.1 (12), 305–324. 10.2118/949305-g

  • CrossRef
  • Google Scholar

• 32

  WangL.ChengD.XiangL.ChenX.XiaZ. (2019). Pressure transient analysis for asymmetrically fractured wells in dual-permeability organic compound reservoir of hydrogen and carbon. Int. J. Hydrog. Energy44 (11), 5254–5261. 10.1016/j.ijhydene.2018.08.082

  • CrossRef
  • Google Scholar

• 33

  WangL.WangX.LuoE.WangJ. (2014). Analytical modeling of flow behavior for wormholes in naturally fractured–vuggy porous media. Transp. Porous Media105 (3), 539–558. 10.1007/s11242-014-0383-7

  • CrossRef
  • Google Scholar

• 34

  WangM.FanZ.DongX.SongH.ZhaoW.XuG (2018). Analysis of flow behavior for acid fracturing wells in fractured-vuggy carbonate reservoirs. Math. Probl. Eng.2018 (PT. 3), 6431910. 10.1155/2018/6431910

  • CrossRef
  • Google Scholar

• 35

  WangQ.YangS.GloverP. W. J.LorincziP.QianK.WangL. (2020). Effect of pore-throat microstructures on formation damage during miscible CO2 flooding of tight sandstone reservoirs. Energy Fuels34 (4), 4338–4352. 10.1021/acs.energyfuels.0c00158.

  • CrossRef
  • Google Scholar

• 36

  WarrenJ.E.RootP. J. (1963). The behavior of naturally fractured reservoirs. Soc. Petrol. Eng. J.3 (3), 245–255. 10.2118/426-pa

  • CrossRef
  • Google Scholar

• 37

  WuY. S.LiuH. H.BodvarssonG. S. (2004). A triple-continuum approach for modeling flow and transport processes in fractured rock. J. Contam. Hydrol.73 (1–4), 145–179. 10.1016/j.jconhyd.2004.01.002

  • CrossRef
  • Google Scholar

• 38

  WuY. S.YuanD.KangZ.FakcharoenpholP. (2011). A multiple-continuum model for simulating single-phase and multiphase flow in naturally fractured vuggy reservoirs. J. Petrol. Sci. Eng.78 (1), 13–22. 10.1016/j.petrol.2011.05.004

  • CrossRef
  • Google Scholar

• 39

  XingC.YinH.LiuK.LiXFuJ. (2018). Well test analysis for fractured and vuggy carbonate reservoirs of well drilling in large scale cave. Energies11 (1), 1–15. 10.3390/en11010080

  • CrossRef
  • Google Scholar

• 40

  XueL.ChenX.WangL. (2019). Pressure transient analysis for fluid flow through horizontal fractures in shallow organic compound reservoir of hydrogen and carbon. Int. J. Hydrog. Energy44 (11), 5245–5253. 10.1016/j.ijhydene.2018.08.084

  • CrossRef
  • Google Scholar

• 41

  YildizT.OzkanE. (1999). Pressure-transient analysis for perforated wells. SPE J.4. 10.2118/56878-PA

  • CrossRef
  • Google Scholar

• 42

  YinC. (2018). Test and analysis on the permeability of induced fractures in shale reservoirs. Nat. Gas. Ind. B5 (5), 513–522. 10.1016/j.ngib.2018.03.006

  • CrossRef
  • Google Scholar