Abstract
In the past few decades, scholars have made great breakthroughs in the study of well test analysis of carbonate rock. The previous studies are based on horizontal wells, straight wells, fractured wells, and inclined wells. With the development of fracturing technology, acid fracturing technology is considered to be the most effective measure to develop carbonate reservoirs. As the carbonate rock is easily dissolved in carbonic acid, multi-branched fractures will be produced near a vertical well. This article presented a semi-analytical model for multi-branched fractures in naturally fractured-vuggy reservoirs for the first time, which laid a theoretical foundation for solving well test analysis for finite conductivity multi-branched fractures. The model can quantify the wellbore flow pressure and applied to obtain more parameters reflecting comprehensive flow characteristics through using history matching procedure. The results were compared with numerical simulation and the existing analytical solutions of a single fracture model. Then in this paper, flow characteristics are recognized and there are five flow regimes found in the type curves, e.g. bi-linear flow region, linear flow region, inter-porosity flow region between vugs and fractures, inter-porosity region between matrix and fractures, and radial flow region. Finally, the influence factors analysis shows fracture number will mainly affect flow behavior of bi-linear flow and linear flow. The angle analysis showed that as the fractures were closer, their interaction became stronger. The conductivity would seriously affect the flow behavior in the early time. Linear flow cannot be observed when the conductivity is less than 1 and bi-linear flow cannot be observed when the conductivity is more than 20. And the effect of fracture length on flow behavior occurs in the early time. Bi-linear flow and linear flow characteristics cannot be observed when the fracture length is reduced.
Similar content being viewed by others
Introduction
Fractured-vuggy rock with complex pore structures has been widely studied by many countries1,2,3,4,5,6,7,8,9. Well test analysis in hydraulically fractured wells has been a topic of interest to the oil industry since early 1950s. Classic analytical and semi-analytical models have been developed for many years. In 1973, Gringarten and Ramey10 first introduced source functions for transient pressure analysis of uniform flux fractures and infinite conductivity fracture wells. Taking account for more realistic cases, Cinco-Ley et al.11 extended the Green’s functions to the transient pressure solutions. Cinco-Ley and Samaniego12 also proposed a bi-linear flow model for analyzing early-time pressure data and type curves for identifying all the flow regimes for wells intersected by finite-conductivity fractures. Furthermore, general solutions considering the effects of dual-porosity effect in the reservoir was presented by Cinco-Ley and Meng13. The above solutions become the basic solutions of production data analysis14,15,16,17 and some authors also considered fracture asymmetry based on the basic solutions18,19,20,21,22.
However, only a single fracture around a vertical well was considered in the above models10,11,12,13,18,19,20,21,22. Carbonate minerals are easy to be dissolved by carbonic acid23,24, which is often used to react with the rock to create a high channel as possible, namely wormhole (Fig. 1), and this process is called as carbonate acidizing. These wormholes are similar to multiple fractures systems from stimulated reservoir volume (SRV)25. Dora Patricia Restrepo in his ph.D’s dissertation25 investigated the pressure behavior of a system containing near wellbore and far field multiple fractures. He also extended the model into dual-porosity media. However, only infinite conductivity model was investigated so that some phenomenon in early time couldn’t be observed.
Fractured-vuggy reservoirs are generally composed of matrix, fractures and vugs systems with varying properties. Generally speaking, both porosity and permeability are different from those systems9,26. Both analytical and numerical methods can be used to simulate the flow behaviors in these systems27,28,29,30,31,32,33,34,35,36,37. However, pressure transient behavior for multi-branched fractures has been rarely reported in previous literatures.
In 2014, Wang et al.34 presented a multi-branched model with an infinite conductivity, however, the flow in the fractures was not considered in their model so that the regime in the early time was not observed. In this paper, a finite conductivity multi-branched fractures model for carbonate reservoirs was studied in details. We focused on a single phase transient flow behavior in fractured-vuggy rock, conceptualized as multiple-continuum medium, consisting permeable natural fractures, low-permeability rock matrix, vugs, and multi-branched fractures. It was assumed that vugs which provided storage space were dispersed throughout fractured carbonate reservoirs, and the fractures were responsible for global flow. Use of analytical techniques, a Laplace-transform model is established and the semi-analytical solutions are successfully derived. The solutions are presented and compared with numerical simulation and the existing analytical solutions of a single fracture model. The presented model can be used to analyze the well test data and production date through using history matching procedure in the future.
Mathematical Models
Basic assumptions
Most of the researchers consider the non-Darcy effect during the process of gas seepage near the wellbore because of the high velocity of gas flow38,39,40,41,42,43,44. However, they also think that although non-Darcy flow may occur, the assumption of Darcy (linear) flow for oil at near the wellbore is reasonable because of relatively low velocity of oil flow. Fractured-vuggy carbonate reservoirs with multi-branched fractures are naturally structured by matrix system, natural fractures system, vugs system. (See Fig. 1a,b). Figure 2 shows that the physical modeling scheme of fractured-vuggy carbonate medium. The assumptions of physical model are as follows:
-
(1)
Sugar cube fractures as seen in classical multiple porosity media models are used in this article4,5,37; shape factors of fracture-vug, fracture-matrix, and vug-matrix could be defined as αfv, αfm, and αvm5,29, respectively.
-
(2)
The reservoir permeability, porosity, compressibility, are constants.
-
(3)
Isothermal process and Darcy flow are assumed.
-
(4)
The total production of multi-branched fractures is equal to total production Q.
-
(5)
The well produces with constant-viscosity and slightly compressible fluid, but the flux of each fracture changes with the time t.
-
(6)
All fractures length may be unequal and the fractures are assumed to be finite conductivity.
-
(7)
The initial pressure is equal to pi.
Reservoir Model
Some reservoir models for vertical wells in carbonate reservoirs have been studied4,5,23,34. Figure 1c shows the flow process. Derivation and description of details for reservoir model can be found in Appendix B of Supplementary material. Through using the superposition principle and double Fourier transform, the dimensionless pressure of Laplace domain is given by
Where \({\tilde{p}}_{fD}\) is the dimensionless fracture pressure of Laplace domain; \({\tilde{q}}_{fD}\) is dimensionless fracture flow rate of Laplace domain; (xD, yD) is dimensionless coordinates. (xfDn, yfDn) are the coordinates of the starting point for sources of nth fracture; θfn is the angle between the nth fracture and x axis; uDn is point source position of nth fracture, integral variable; LfDn is dimensionless length of nth fracture; Nf is fracture number; s is dimensionless time variable of Laplace domain; ωf is fracture storage coefficient; ωv is vug storage coefficient; ωm is matrix storage coefficient; λfv is fracture-vug inter-porosity flow coefficient; λfm is fracture-matrix inter-porosity flow coefficient; λvm is vug-matrix inter-porosity flow coefficient; K0(x) is the modified Bessel function (2nd kind, 0 order).
Fracture Flow Model
The well is fixed in the endpoint. Linear flow occurs in the fractures (Fig. 3). Establishment and solution of fracture flow model can be found in Appendix C of Supplementary material. The pressure drop at kth segment on the nth fracture is given by
where
and
Where \({\tilde{p}}_{fDnk}\) is the dimensionless pressure of Laplace domain at kth segment on the nth fracture; \({\tilde{p}}_{fDn}(0,s)\) is dimensionless average pressure of Laplace domain at the nth fracture; \({\tilde{q}}_{fDni}\) is the dimensionless fracture flux of Laplace domain at ith segment on the nth fracture; ΔzDn is the dimensionless length of fracture segment on the nth fracture; ZfDnk is the dimensionless midpoint at kth segment on the nth fracture; ZwhfDn is the dimensionless well location on the nth fracture; CfDn is the fracture conductivity on the nth fracture.
Semi-analytical Solution
The reservoir model and fracture model are presented in details in Appendix B and C of Supplementary material. Now, we divide per fracture into Nsn segments to Eq. (1), thus, the jth pressure segment of nth fracture can be given by
where
The flow rate of each fracture should be equal to the sum of the flow flux of all fracture segment in this fracture, therefore, we have
The pressure should be equal for all fractures connected with horizontal wells, we have
In addition, total rate of all fractures connected with horizontal wells should be equal to flow rate of horizontal wells, we have
Coupling the multi-branched fractures and reservoir-flow models at the fracture faces, we obtain a system of Nts + Nf + Nf equations with Nts + Nf + Nf unknowns. Nts is the total number of fracture segments, and Nf is the total number of multi-branched fractures. Equations 7 to 11 constitute a system of linear equations that can be solved by Gauss elimination. Once the equations are solved, we can obtain the flux of each fracture segment. By substituting the flux of all fractures segment into equation 7, we can obtain bottom-hole pressure.
Results
Comparison with numerical simulation
Figure 4 shows the comparison between results from this work and those computed data from commercial software Eclipse with four fractures. All the coordinates of the starting point for sources (xfD1, yfD1), (xfD2, yfD2), (xfD3, yfD3) and (xfD4, yfD4) are set to (3,3). Initial fracture conductivity CfD1, CfD2, CfD3 and CfD4 for all the fractures is set to 1, 10, 100 for three cases to compare with numerical simulation. Fracture length LfD1, LfD2, LfD3 and LfD4 for all the fractures is set to 0.5. The fracture angles θf1, θf2, θf3 and θf4 are set to π/2, π, 3π/2, 2π, respectively. The g(s) in Eq. (8) is set to be 1. The software is used to simulate the oil, gas and water flow during the reservoir development. The code is written using finite difference method. The fracture number is set to 4 and fracture conductivity CfD is set to 1, 10, 100 in the software Eclipse. The comparison suggests that our model can be well validated by the software. However, for early time, the results from numerical simulation cannot well agree with the semi-analytical results and the accuracy of numerical simulation is controlled by the number of grids near the wellbore. Both pressure and pressure derivative data deviates from 1/4 slope straight line.
Comparison with existing an analytical solution for a single fracture
Figure 5 shows the comparison with existing analytical solutions presented by Cinco et al. (1988). The existing analytical solutions are special cases of our presented model. Therefore, the coordinates of the starting point for sources (xfD1, yfD1) and (xfD2, yfD2), are set to (3,3). Initial fracture conductivity CfD1 and CfD2 for all the fractures is set to 0.2π, π, 2π, 10π, and 1000π for five different cases. Fracture length LfD1 and LfD2 for both the fractures is set to 0.5. The fracture angles θf1 and θf2 are set to 0, π, respectively. The comparison shows that our model can be well validated by the solutions from literatures. However, in the literature, the solution is only considered into a single fracture as shown in Fig. 5. Therefore, our presented model is a general model, which could consider more complex branched fractures and inter-porosity flow effect.
Comparison with numerical simulation for a single fracture with non-Darcy flow
Figure 6 shows the graph of dimensionless pressure vs. dimensionless time for different values of non-Darcy Forchheimer number FND,F at the same fracture conductivity. We compared the presented semi-analytical solution with Darcy flow and numerical simulation with non-Darcy flow (different Forchheimer number values). The coordinates of the starting point for sources (xfD1, yfD1) and (xfD2, yfD2), are set to (3,3). Initial fracture conductivity CfD1 and CfD2 for all the fractures is set to 10. Fracture length LfD1 and LfD2 for both the fractures is set to 0.5. The fracture angles θf1 and θf2 are set to 0, π, respectively. The Forchheimer number values are set to 1, 3, 5, respectively. As expected, the Forchheimer number, FND,F, has a observable effect on transient pressure behavior. In the early time when tD < 1, A bigger Forchheimer number will lead to a stronger non-Darcy flow. We can see that in the early time the pressure becomes larger with the increase of Forchheimer number FND,F, which implies non-Darcy flow will lead to a bigger pressure depletion. However, in the late time when tD > 1, there is almost no error between Darcy flow and non Darcy flow.
Flow Characteristics analysis
The transient flow characteristics are analyzed by type curves, which can be used to analyze transient pressure so as to recognize the fluids flow characteristics in fractured-vuggy reservoirs21,22. Fracture number is set to 4. All the coordinates of the starting point for sources (xfD1, yfD1), (xfD2, yfD2), (xfD3, yfD3) and (xfD4, yfD4) are set to (3,3). Initial fracture conductivity CfD1, CfD2, CfD3 and CfD4 for all the fractures is set to 2π for three cases to compare with numerical simulation. Fracture length LfD1, LfD2, LfD3 and LfD4 for all the fractures is set to 0.5. The fracture angles θf1, θf2, θf3 and θf4 are set to π/2, π, 3π/2, 2π, respectively. Fracture-vug, fracture-matrix, and vug-matrix inter-porosity flow coefficients are set to 3, 0.01 and 0.01, respectively. Vugs, fracture, matrix storage coefficients are set to 0.04, 0.005 and 0.955, respectively. Five main flow regions are well shown in Fig. 7.
-
(1)
A 1/4 slope straight line on the pressure derivative curve represents bi-linear flow region.
-
(2)
A 1/2 slope straight line on the pressure derivative curve represents linear flow region.
-
(3)
The first V-shaped segment on the pressure derivative curve represents inter-porosity flow region from vug system to the fracture system shows. After bi-linear flow and linear flow, the fluid in the fracture gradually decreases, so the fluid in the vugs will flow towards the fracture 46.
-
(4)
The second V-shaped segment on the pressure derivative curve represents inter-porosity flow region from matrix system to fracture system and vug system to matrix system. During this period, the fluid in the fracture is further reduced, and the fluid in the matrix also begins to flow towards the fracture. Because of the large amount of fluid stored in the vugs, the part of fluid in the vugs also begins to flow towards the fracture and other part of fluids in vugs flow towards the matrix9.
-
(5)
The radial flow regime shows a zero slope and 0.5 constant straight line on the pressure derivative curve. At the late time of the flow, the fluid far from the wellbore began to flow. At this time, the fluid around the wellbore converge toward the wellbore, which is a radial flow. During this period, fluid exchange in matrix, fracture and vugs is in dynamic balance.
Parameter Influence
Figures 8–12 show the influence of parameters (fracture number N, fracture angle θ, conductivity CfD, and fracture length LfDn) on pressure and pressure derivative curves. The initial parameters are listed here. The initial fracture number is set to 4. All the coordinates of the starting point for sources (xfD1, yfD1), (xfD2, yfD2), (xfD3, yfD3) and (xfD4, yfD4) are set to (3,3). Initial fracture conductivity CfD1, CfD2, CfD3 and CfD4 for all the fractures is set to 2π for three cases to compare with numerical simulation. Fracture length LfD1, LfD2, LfD3 and LfD4 for all the fractures is set to 0.5. The fracture angles θf1, θf2, θf3 and θf4 are set to π/2, π, 3π/2, 2π, respectively. Fracture-vug, fracture-matrix, and vug-matrix inter-porosity flow coefficients are set to 3, 0.01 and 0.01, respectively. Vugs, fracture, matrix storage coefficients are set to 0.04, 0.005 and 0.955, respectively.
As shown in Fig. 8, the dimensionless pressure decreases with the increase of the fracture number, which implies that it is the more the number of fractures, the smaller the pressure depletion. Therefore, pressure depletion could be reduced and well production can be improved through producing more fractures. The pressure derivative curves show that fractures number is mainly affecting flow behavior of bi-linear and linear flow region.
Figure 9 shows that facture angle affects the pressure responses and fracture number is 4. It is shown that the smaller the fracture angle is, the greater the pressure loss is. Pressure derivative curves show that when the fracture angle is greater than 1°, the bi-linear flow and the linear flow can be affected importantly by fracture angle. When the angle of the fractures is greater than 50°, the flow characteristic in early time is nearly not affected by the angle.
Figure 10 shows the effect of the conductivity on type curves. Both the pressure and its derivative curves show this effect mainly occurs in the early time. The pressure curve shows the dimensionless pressure increases as the conductivity decreases, which means a low conductivity will cause big pressure depletion. Pressure derivative curve shows bi-linear flow will end early as the conductivity increases. Linear flow cannot be observed when the conductivity is less than 1 and bi-linear flow cannot be observed when the conductivity is more than 20.
If the conductivity for each fracture is different, how the conductivity affects the type curves? As shown in Fig. 11, the uniform conductivity for case 3 exhibits the complete flow characteristics, the starting time of the linear flow becomes late and the pressure depletion increases as the each fracture conductivity decreases. This point is coincided with the analysis above under the uniform conductivity. Figure 12 shows the effect of the fracture length on type curves. The curves exhibit the complete flow characteristics when all fracture length is set to be 0.5. However, the pressure depletion will increase as fracture length is reduced. The effect of fracture length on flow behavior occurs in the early time. Bi-linear flow and linear flow characteristics are not observed when the fracture length is reduced. In summary, effect of fracture number on type curves is big and influence of other parameters on them is small.
The model proposed in this paper can well reveal the formation bilinear flow and the linear flow pattern and effects of some important parameters on bilinear flow and linear low are presented. In well test, because of the short test time, late data cannot be obtained. Therefore, the bilinear flow and the linear flow patterns are very important for early interpretation of production data. In addition, in this paper, we analyzed the influence of fracture correlation parameters on type curves. Because the fracture related parameters are all affecting flow near wellbore, we can see that the parameters only affect bilinear flow and linear flow. However, the region 4 and region 5 can be affected by inter-porosity flow coefficient and storage coefficient, which have been analyzed in previous literatures9,34. However, flow pattern in early time is not revealed and production data in early time cannot be analyzed in previous literatures. We presented a comprehensive model considering all possible flow patterns in this paper and that we can analyze the production data of bilinear flow and linear flow is the important contribution of this paper.
Discussions
Advantages of the presented semi-analytical model
The semi-analytical method effectively eliminates the truncation error and has the same order of accuracy as the analytical method if the property varication is not considered. The semi-analytical method does not require discretization of reservoir at the first step of simulation and it also it does not need discretization of time. It can be done at any time. The disadvantage is that at present our semi-analytical model is now limited to single-phase flows. However, in the future, we are trying to study the semi-analytical solution of two-phase flow.
Application of a real field case
What is the most striking is that this semi-analytical model can be applied to obtain more parameters reflecting comprehensive flow characteristics through using history matching procedure. A carbonate reservoir can be characterized by many parameters, i.e. fracture conductivity, fracture number, fracture angle, fracture length. If those parameters are unknown, the determination of unknown parameters is in fact a subject of inverse problems. Determining those unknown parameters could use the present solution combining with the algorithm of auto history matching45.
In this section, a field case is selected to validate the presented solution. Well A located in the western oil field of China is selected. The well test data is shown in Fig. 13. The reservoir thickness is 25 m and the effective porosity (matrix + fractures + vugs) is 0.12. The oil viscosity is 0.00152 Pa▪s and the total compressibility is 0.0014 MPa−1. The production of well is 120 m3/d. According to well test data, we have drawn the pressure and pressure derivative curves in Fig. 13. Pressure derivative curves show 1/4 and 1/2 slopes straight lines, which implies that bi-linear flow and linear flow occur in the early time. However, the inter-porosity flow region is not found in the figure. Therefore, g(s) is set to 1 in our presented model. Furthermore, we can determine those unknown parameters use the present solution combining with the algorithm of auto history matching45. The final matching results are given: The permeability kf is 23.2Md; fracture half length is 430 m; fracture number is 4; fracture conductivity is 99.7. From fracture conductivity, we can see that fracturing effects are good. The matching results also implies that the presented model can be used to explain the well test data in the early time, which cannot be done in previous study34.
Conclusions
Based on our work, the following conclusions can be drawn.
-
(1)
In this paper, firstly a semi-analytical solution for finite conductivity multi-branched fractures is presented.
-
(2)
The presented semi-analytical solutions are validated by numerical simulation and existing analytical solutions.
-
(3)
The flow Characteristics are also proposed in this article, and flow in carbonate reservoirs can be divided into five regions, e.g. bi-linear flow region, linear flow region, inter-porosity flow region between vugs and fractures, inter-porosity region between matrix and fractures, and radial flow region.
-
(4)
The influence factors analysis shows fracture number will be the main factor influencing on flow. The present solution combining with the algorithm of auto history matching to obtain reservoir parameters, furthermore, production performance analysis can be done.
References
Kossack, C. A. & Gurpinar, O. A methodology for simulation of vuggy and fractured reservoirs. In SPE Reservoir Simulation Symposium, Houston, Texas. Society of Petroleum Engineers. https://doi.org/10.2118/66366-MS (2001, February 11–14).
Rivas-Gomez, S. & Gonzalez-Guevara, J. A. Numerical simulation of oil displacement by water in a vuggy fractured porous medium. In SPE Reservoir Simulation Symposium, Houston, Texas. Society of Petroleum Engineers. https://doi.org/10.2118/66386-MS (2001, February 11–14).
Camacho-Velázquez, R., Vásquez-Cruz, M. & Castrejón-Aivar, R. Pressure-transient and decline-curve behaviors in naturally fractured vuggy carbonate reservoirs. SPEREE. 8(2), 95–112 (2005).
Wu, Y. S., Qin, G. & Ewing, R. E. A multiple-continuum approach for modeling multiphase flow in naturally fractured vuggy petroleum reservoirs. In International Oil & Gas Conference and Exhibition in China, Beijing, China. Society of Petroleum Engineers. https://doi.org/10.2118/104173-MS (2006, December 5–7).
Wu, Y. S., Economides, C. E. & Qin, G. A triple-continuum pressure-transient model for a naturally fractured vuggy reservoir. In SPE Annual Technical Conference and Exhibition, Anaheim, California, U.S.A. Society of Petroleum Engineers. https://doi.org/10.2118/110044-MS (2007, November 11–14).
Jazayeri Noushabadi, M. R., Jourde, H. & Massonnat, G. Influence of the observation scale on permeability estimation at local and regional scales through well tests in a fractured and karstic aquifer. J. Hydrol. 403(3-4), 321–336 (2011).
Popov, P. & Qin, G. Multiphysics and multiscale methods for modeling fluid flow through naturally fractured carbonate karst reservoirs. SPEREE. 12(2), 218–231 (2009).
Rzonca, B. Carbonate aquifers with hydraulically non-active matrix: a case study from Poland. J. Hydrol. 355(1–4), 202–231 (2008).
Guo, J. C., Nie, R. S. & Jia, Y. L. Dual permeability flow behavior for modeling horizontal well production in fractured-vuggy carbonate reservoirs. J. Hydrol. 464-465, 281–293 (2012).
Gringarten, A. C. & Ramey, H. J. Jr. The use of source and green’s function in solving unsteady-flow problems in reservoirs. SPE J. 13(05), 285–296 (1973).
Cinco, L. H., Samaniego, V. F. & Dominguez, N. Transient pressure behavior for a well with a finite-conductivity vertical fracture. SPE J. 18(4), 253–264 (1978).
Cinco, L. H. & Samaniego, V. F. Transient pressure analysis for fractured wells. J. Pet. Technol. 33(9), 1749–1766 (1981).
Cinco, L. H. & Meng, H. Z. Pressure transient analysis of wells with finite conductivity vertical fractures in double porosity reservoirs. In SPE Annual Technical Conference and Exhibition, Houston, Texas. Society of Petroleum Engineers. https://doi.org/10.2118/18172-MS (1988, October 2–5).
Tiab, D. Analysis of pressure derivative without type-curve matching-skin and wellbore storage. J. Pet. Sci. Eng. 12(3), 171–181 (1995).
Agarwal, R. G., Gardner, D. C., Kleinsteiber, S. W. & Fussell, D. D. Analyzing well production data using combined type curve and decline curve concepts. SPE Form. Eval. 2(5), 478–486 (1998).
Liang, C., Cao, Y. C., Liu, K. Y., Wu, J. & Hao, F. Diagenetic variation at the lamina scale in lacustrine organic-rich shales: Implications for hydrocarbon migration and accumulation. GEOCHIM COSMOCHIM AC. 229, 112–128 (2018).
Tiab, D. Analysis of derivative data of hydraulically fractured wells by the Tiab’s direct synthesis technique. J. Pet. Sci. Eng. 49(1-2), 1–21 (2005).
Rodriguez, F., Cinco, L. H. & Samaniego, V. F. Evaluation of fracture asymmetry of finite-conductivity fractured wells. SPE Form. Eval. 8(2), 233–239 (1992).
Wang, L., Dai, C. & Xue, L. A semi-analytical model for pumping tests in finite heterogeneous confined aquifers with arbitrarily shaped boundary. Water Resour Res. https://doi.org/10.1002/2017WR022217 (2018).
Berumen, S., Tiab, D. & Rodriguez, F. Constant rate solutions for a fractured well with an asymmetric fracture. J. Pet. Sci. Eng. 25, 49–58 (2000).
Wang et al. Simulation of pressure transient behavior for asymmetrically finite-conductivity fractured wells in coal reservoirs. Transp. Porous. Med. 97(3), 353–372 (2013).
Wang, L. Wang X.D.2013b. Type curves analysis for asymmetrically fractured wells. J Eenergy Resource-ASME. 136, 023101 (2013).
Liu, M., Zhang, S. C. & Mou, J. Y. Wormhole propagation behavior under reservoir condition in carbonate acidizing. Trans. Porous. Med. 96(1), 203–220 (2013).
Fredd, C. N. Dynamic model of wormhole formation demonstrates conditions for effective skin reduction during carbonate matrix acidizing. In SPE Permian Basin Oil and Gas Recovery Conference, Midland, Texas. Society of Petroleum Engineers. https://doi.org/10.2118/59537-MS (2000, March 21–23).
Dora, P. R. Pressure behavior of a system containing multiple vertical fractures, PhD thesis, The University of Oklahoma (2008).
Mai, A. & Kantzas, A. Porosity distributions in carbonate reservoirs using low-field NMR. J. Canadian Petroleum Technol. 46(7), 30–36 (2007).
Akram, A. H., Gherryo, Y. S. & Ali, S. M. Dynamic behavior of a fissured dual-carbonate reservoir modeled with DFN. In North Africa Technical Conference and Exhibition, Cairo, Egypt. Society of Petroleum Engineers. https://doi.org/10.2118/127783-MS (2010, February 14-17).
Arana, V., Pena-Chaparro. O. & Cortes-Rubio, E. A practical numerical approach for modeling multi-porosity naturally fractured reservoirs. In Latin American and Caribbean Petroleum Engineering Conference, Cartagena de Indias, Colombia. Society of Petroleum Engineers. https://doi.org/10.2118/123006-MS (2009, May 31).
Cobett, P. W. M., et al Limitations in the numerical well test modeling of fractured carbonate rocks. In SPE EUROPEC/EAGE Annual Conference and Exhibition, Barcelona, Spain. Society of Petroleum Engineers. https://doi.org/10.2118/130252-MS (2010, June 14–17).
Djatmiko, W. & Hansamuit, V. Pressure buildup analysis in karstified carbonate reservoir. In SPE Asia Pacific Oil and Gas Conference and Exhibition, Brisbane, Queensland, Australia. Society of Petroleum Engineers. https://doi.org/10.2118/132399-MS (2010, October 18–20).
Gulbransen, A. F., Hauge, V. L. & Lie, K. A. A multiscale mixed finite element method for vuggy and naturally fractured reservoirs. SPE J. 15(2), 395–403 (2010).
Leveinen, J. Composite model with fractional flow dimensions for well test analysis in fractured rocks. J. Hydrol. 234(3-4), 116–141 (2000).
Yao, J., Huang, Z. Q., Li, Y. Z., Wang, C. C. & Lv, X. R. Discrete fracture-vug network model for modeling fluid flow in fractured vuggy porous media. In International Oil and Gas Conference and Exhibition in China, Beijing. Society of Petroleum Engineers. https://doi.org/10.2118/130287-MS (2010, June 8–10).
Wang, L., Wang, X. D., Luo, E. H. & Wang, J. L. Analytical modeling of flow behavior for wormholes in naturally fractured-vuggy porous media. Trans. Porous. Med. 105(3), 539–558 (2014).
Liu, R. C., Jiang, Y. Z., Huang, N. & Sugimoto, S. Hydraulic properties of 3D crossed rock fractures by considering anisotropic aperture distributions. Advances in Geo-Energy Research. 2(2), 113–121 (2018).
Huang, S. Y., Wu, Y. Y., Meng, X. B., Liu, L. W. & Ji, W. Recent advances on microscopic pore characteristics of low permeability sandstone reservoirs. Advances in Geo-Energy Research. 2(2), 122–134 (2018).
Warren, J. E. & Root, P. J. The behavior of naturally fractured reservoirs. SPE J. 228, 245–255 (1963).
Guppy, K. H., Cinco-Ley, H., Ramey, H. J. Jr. & Samaniego-V., F. Non-darcy flow in wells with finite-conductivity vertical fractures. SPE J. 22(5), 681–698 (1982).
Guppy, K. H., Cinco-Ley, H. & Ramey, H. J. Jr. Pressure buildup analysis of fractured wells producing at high flow rates. J. Pet. Technol. 34(11), 2656–2666 (1982).
Gidley, J. L. & John, L. A method for correcting dimensionless fracture conductivity for non-Darcy flow effects. SPE Prod. Eng. 6(4), 391–394 (1991).
Roldan, J., Camacho-V., R. & Samaniego, F. Well testing in solution gas-drive systems under non-laminar flow conditions. In International Petroleum Conference and Exhibition of Mexico, Mexico. Society of Petroleum Engineers. https://doi.org/10.2118/35344-MS (1996, March 5–7).
Umnuayponwiwat, S. & Ozkan, E. Water and gas coning toward finite-conductivity horizontal wells: cone buildup and breakthrough. In SPE Rocky Mountain Regional/Low-Permeability Reservoirs Symposium and Exhibition, Denver, Colorado. Society of Petroleum Engineers. https://doi.org/10.2118/60308-MS (2000, March 12–15).
Ramirez, B., Kazemi, H. & Ozkan, E. Non-Darcy flow effects in dual-porosity, dual-permeability naturally fractured gas condensate reservoirs. In SPE Annual Technical Conference and Exhibition, Anaheim, California, U.S.A. Society of Petroleum Engineers. https://doi.org/10.2118/109295-MS (2007, November 11–14).
Zeng, F. H. Modeling of non-Darcy flow in porous media and its application, PhD thesis, University of Regina (2008).
Nelder, J. A. & Mead, R. A simplex method for function minimization. Computer Journal. 7, 308–313 (1965).
Acknowledgements
This work is funded by China Postdoctoral Science Foundation (Grant No. 2017M620527) and the National Science and Technology Major Project of China (Grant No. 2016ZX05034001-007).
Author information
Authors and Affiliations
Contributions
As the first author, L. Wang established the novel mathematical model and wrote the main manuscript text. As the second author, X. Chen made substantial contributions to the conception/design of the work. As the third author, Z. Xia revised the manuscript. All the authors discussed the results and critically reviewed the manuscript.
Corresponding author
Ethics declarations
Competing Interests
The authors declare no competing interests.
Additional information
Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Electronic supplementary material
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Wang, L., Chen, X. & Xia, Z. A Novel Semi-Analytical Model for Multi-branched Fractures in Naturally Fractured-Vuggy Reservoirs. Sci Rep 8, 11586 (2018). https://doi.org/10.1038/s41598-018-30097-2
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-018-30097-2
This article is cited by
-
Pressure transient analysis for a fractured well in a stress-sensitive tight multi-medium oil reservoir
Frontiers of Earth Science (2021)
-
Velocity dependent up-winding scheme for node control volume finite element method for fluid flow in porous media
Scientific Reports (2020)
-
A new analytical model for flow in acidized fractured-vuggy porous media
Scientific Reports (2019)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.