#### Productivity Analysis of Volume Fractured Vertical Well Model in Tight Oil Reservoirs

Productivity Analysis of Volume Fractured Vertical Well Model in Tight Oil Reservoirs
Jiahang Wang,1 Xiaodong Wang,1 Wenli Xu,1 Cheng Lu,1,2 Wenxiu Dong,1 and Cong Zhang1,2
1China University of Geosciences, Beijing 100083, China
2The Key Laboratory of Unconventional Petroleum Geology, CGS, Beijing 100029, China
Correspondence should be addressed to Cheng Lu; moc.qq@010071244
Received 18 October 2016; Accepted 27 December 2016; Published 26 January 2017
Academic Editor: Jian G. Zhou
Copyright © 2017 Jiahang Wang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Abstract
This paper presents a semianalytical model to simulate the productivity of a volume fractured vertical well in tight oil reservoirs. In the proposed model, the reservoir is a composite system which contains two regions. The inner region is described as formation with finite conductivity hydraulic fracture network and the flow in fracture is assumed to be linear, while the outer region is simulated by the classical Warren-Root model where radial flow is applied. The transient rate is calculated, and flow patterns and characteristic flowing periods caused by volume fractured vertical well are analyzed. Combining the calculated results with actual production data at the decline stage shows a good fitting performance. Finally, the effects of some sensitive parameters on the type curves are also analyzed extensively. The results demonstrate that the effect of fracture length is more obvious than that of fracture conductivity on improving production in tight oil reservoirs. When the length and conductivity of main fracture are constant, the contribution of stimulated reservoir volume (SRV) to the cumulative oil production is not obvious. When the SRV is constant, the length of fracture should also be increased so as to improve the fracture penetration and well production.
1. Introduction
Reservoir stimulation technologies have been widely applied to develop tight oil reservoirs. Volume fracturing technique is one of these methods which has been widely used to improve the productivity of low permeability and tight reservoirs. After repeated acid fracturing treatment to the fractured brittle reservoirs, hydraulic fracture, natural fracture, and shear cracks are mutually staggered and form a certain stimulated zone of joint fracture network near the wellbore, which changes the flow pattern, reduces the flow resistance, and improves production of a single well [1–4]. Testing and evaluating fracture network reconstruction along the well and productivity behavior are essential for improving the performance of production well in tight reservoirs after volume fracturing. Rate decline analyzing is one of the key methodologies to estimate reservoir parameters, such as permeability, porosity, length, width, and skin factor.
Due to the advanced techniques in fracturing, the behavior of rate decline curves in a fractured well has attracted increasing attention recently. In terms of numerical simulation, Al-Salem et al. established the model by using the vertical and horizontal orthogonal crack network to approximately substitute volume reconstruction [5], and this model has been widely used since then [6–9]. Combining the micro seismic exploration results, Arvind approximately simulated the volume and the degree of the reconstruction region around wells [10]. However, Du et al. described the transformation region volume by using the Kazemi dual medium model [11]. In terms of the analytic model, Liu et al. [12, 13] and Lei and Gang [14] described the fracture distribution of volume transformation region of vertical wells by using the fractal theory, and the production of cold and heavy oil with carrying sand is studied based on their model. Recently, a composite reservoir model with fractal permeability was applied to evaluate the productivity of tight oil reservoirs [15–17].
Compared with the analytical method, the numerical simulation method is able to deal with the complicated seepage problem to a large extent by the grid division, while the procedure is complicated and causes large amount of computation resources. The fractal theory can describe the spatial distribution of fracture better, but it is not suitable for the pressure transmission behavior and the artificial fracture parameter optimization research. Liao and Chen described the pressure transient analysis of volume fracturing well without considering the wellbore storage effect and skin effect, and 5 flow regimes were recognized in the transient pressure type curves [17]. So far, there has not been an overall well test interpretation model for the vertical wells with stimulated volume in fractured tight reservoirs.
In this work, we propose a semianalytical model to simulate the productivity of volume fractured vertical well in a tight oil reservoir by Laplace transform and Stehfest numerical inversion [18]. By modeling the stimulated volume as the dual medium and coupling the stimulated formation with the discretized artificial fracture, a semianalytical solution is obtained. Finally the effects of some sensitive parameters, such as storage capacity ratio, crossflow coefficient, fracture conductivity, fracture length, stimulated volume, and skin factor on the type curve,s are also analyzed extensively.
2. Mathematical Model
A stimulated volume with joint network is formed near the wellbore in brittle tight reservoirs after repeated acid fracturing treatments. The stimulated fracture network can be normally subdivided into two parts, the inner artificial main fracture and the outer natural fractured zone, respectively. The fluid flow in the main fracture is linear and follows Darcy’s law and the classical Warren-Root model [19] is used to describe the fracture distribution and seepage flow in the outer area where there are no artificial fractures and the permeability is very low (<0.1 mD) because of low connectivity. Affected by the extension of the artificial fracture and the brittle shear of reservoir rock in outer area, the artificial fracture and natural fractures are arranged in a crisscross pattern and change the flow pattern mainly in fractures. In this work, the fluid supplied to the stimulated volume region by natural fractures is neglected since it is much smaller compared with that from the artificial main fracture [20, 21]. In summary, the model assumptions are listed below:(1)The model is homogeneous and isotropic along the radial direction.(2)The pressure is constant and both fluid and rock are slightly compressible.(3)Fractures are the main flow channel, and the seepage flow is laminar and isothermal.(4)The conductivity of vertical artificial fracture is finite and the fracture fully penetrated the formation with the height equal to the thickness of the reservoir.Let us consider a volume fractured well in a circular closed reservoir. The main vertical fracture with finite conductivity has a half length , a width , and a permeability , and penetrates the formation vertically. The classical Warren-Root model is used to simulate the microcracks produced by volume fracturing in the reservoir formation. The reservoir is composed of a fracture network and matrix blocks. The fracture network possesses a bulk fracture porosity and total compressibility . The matrix blocks are slabs of thickness , permeability , porosity , and total compressibility .
The fluid in the reservoir is slightly compressible and its viscosity is . The flow process in the system can be studied by breaking up the medium into three parts and taking the interaction among different parts into account. These regions are hydraulic fracture in the inner area and reservoir fracture network and reservoir matrix in the outer area (see in Figure 1).
Figure 1: Schematic diagram: (a) stimulated reservoir volume model. (b) Physical modeling scheme of artificial main fracture. (c) The system of matrix. (d) Classic Warren-Root dual medium model.
2.1. Region I
The flow within the hydraulic fracture is considered to be linear because the fracture width, , is much smaller than fracture length and fracture height. It is assumed that the fluid flow into the wellbore takes place only through the hydraulic fracture and the fluid flow from the reservoir into the hydraulic fracture occurs only through the reservoir fracture network because is much larger than . Figure 1 illustrates the characteristic of this model, and here is the flow rate going to the fracture per unit length.
Cinco et al. have demonstrated that the compressibility of the hydraulic fracture can be neglected for practical purposes because the fracture volume is very limited [22]. Hence the flow within the fracture can be considered to be incompressible. Under these conditions, the transient flow in the hydraulic fracture can be described by the following equations in terms of dimensionless variables:Inner boundary conditionOuter boundary conditionThe flow correlation formula for surface of the fracture isWe can obtain the solution in Laplace domain through combining (1)–(4); that is,where
2.2. Region II
As mentioned above, the stimulated formation is full of microcracks. The reservoir is represented by a fracture network and matrix blocks. It is assumed that the characteristics of both the fracture network and matrix blocks remain constant. The flow from the reservoir into the hydraulic fracture only occurs through the fracture network, as generally considered in the literatures for dual medium reservoirs.
The transient flow in the formation can be described bybecauseWe can get the equation below:Equation (7) can be further simplified:Inner boundary conditionOuter boundary conditionwhereCombining Equations (7)–(12) can obtain the point source solution of formation transient flow:whereFurther, the solution of plane source is obtained by integrating point source in terms of Bessel functions. The pressure distribution of this system is then given belowCombining Equations (5) and (16), we can obtainConsidering the fracture symmetry,Equation (17) becomesEquation (19) gives the transient solution of the finite conductivity vertical fracture.
2.3. Solution
Assuming the fracture can be divided into segments, the first part on the right side of (19) can be written asand the second part can be expressed asBy using variable substitution method, the integral form of and can be expressed as follows:where(1)If ,(2)If ,(3)If ,(4) And the second-order integral of (19) can be expressed asAssumeAnd then
:::.
In addition to above expressions, due to steady flow, we also haveCombining (30)–(33), we can get the linear equations as follows:where is coefficient of .
The wellbore pressure of a hydraulic fractured well under constant production in naturally fractured reservoirs is obtained after solving (34) by using a Gaussian elimination approach.
In Laplace domain, we can easily consider the effect of skin factorAccording to the Duhamel principle [23], the pressure solution under constant rate and the rate solution under constant pressure have the following correlation in Laplace domain:Combining (35) and (36), we can obtain the correlation between dimensionless rate and dimensionless time for any given parameters.
3. Results and Discussion3.1. Validation of Solution
Riley [24] gave an analytical solution for elliptical finite conductivity fractures without volume fracturing. To validate the solution presented in this paper, we compared our solution with Riley’s result. In our model, storage capacity ratio , crossflow coefficient , and skin factor are considered to be equal to zero. Figure 2 shows the comparison of the two solutions under different fracture conductivity . The good agreement validates the solution obtained in this work.
Figure 2: The comparison for the results of this paper and Riley [24].
Combined with fracturing design, micro seismic detection and well test data, the parameters of a volume fractured well in tight oil reservoir are obtained, shown in Table 1. In Figure 3 it shows that the fractured well has produced for more than 200 days. After a short increasing period and stable period, the production began to decline. The calculation of established model is used to fit the actual production data at the decline stage. The calculated results are in good agreement with the actual production.
Table 1: Parameters’ value for influential factors analysis.
Figure 3: Actual and theoretical production data fitted curves.
3.2. Typical Curves Analysis
The pressure and its derivative curves are presented in Figure 4, which shows basic flow characteristics for a volume fractured well in tight oil reservoir under different parameters by using Stehfest numerical inversion. The parameters are given as , , and ; , , and ; (infinite boundary), , and ; , , and ; , , and . As shown in Figure 4, the flow can be divided into 7 stages:(1)Stage A: Early Bilinear Flow (Artificial Fracture and Fractured Reservoirs Near the Wellbore). In this stage, the segment has a straight line with slope equal to 1/4, demonstrating the bilinear flow region (see Figure 4). In this region, fluid flows through fracture to wellbore and from reservoirs to fracture at the same time. This region could be identified only if the fracture conductivity is relatively small.(2)Stage B: Early Coupled Boundary Flow. In this stage, when the effect of volume fracturing is good, the artificial main fracture conductivity is much larger than that of the double medium fracture system. The fluid in the artificial main fracture reaches wellbore rapidly. However, the double medium fracture system cannot provide adequate fluid supply. Both pressure and pressure derivate curves increase, similar to the transient pressure response of weak energy supply or closed boundary reservoirs.(3)Stage C: Early Linear Flow (Fractured Reservoir Near the Wellbore). The pressure and pressure derivative curves are both straight lines with slope equal to 1/2, which clearly demonstrates early linear flow. In this region, linear fluid flow from formation to the artificial fracture reduces the seepage resistance. The early and flow regions do not necessarily occur for each fracturing treatment, which depends on the conductivity of artificial fracture flow. As shown in Figure 4, the bilinear flow is more obvious in the case of smaller fracture conductivity.(4)Stage D: Mid Radial Flow (Microcracks in the Fractured Formation). The segment has a straight line with 0.5 constant.(5)Stage E: Mid Bilinear Flow or Pseudo Steady-State Flow. The larger the crossflow factor is, the earlier the crossflow happens. Before fluid crossflow between microcracks and formation occurs, it exhibits pseudo steady flow or mid linear flow briefly which is affected by the crossflow coefficient . Specifically, if is small it will exhibit pseudo steady flow when the pressure reaches the boundary or linear flow when the pressure wave disturbance does not reach the boundary.(6)Stage F: Mid Crossflow (Matrix and Fracture). Due to the low permeability of matrix and the slow pressure drop, crossflow will occur between the matrix and fracture. And the pressure derivative curve is concave. Compared to the conventional dual media the crossflow will happen earlier. At the same time, due to the dimensionless setting, the crossflow coefficient is 2 to 4 orders of magnitude larger than that of the conventional dual medium.(7)Stage G: Late Pseudo Steady Flow. For infinite outer boundary, the pressure derivative curve is a horizontal line. However, for closed outer boundary, the slopes of pressure and pressure derivative cures are 1. In some cases, affected by the storage capacity ratio and crossflow coefficient the medium segment has different flow characteristics: D + E or D + F.
Figure 4: Flow stage division.
3.3. Blasingame Type Curves and Discussions
The general definitions of the base decline for type curve variables can be given byThe definitions of (37) were presented by Fetkovich [26]. However, to introduce material balance time, we use (35) as a constant rate solution [27] and then use result as . To eliminate multiple solutions and errors, integral average method of rate was created by Blasingame et al. [28]. The auxiliary variables in this method typically used for decline type curve analysis are given by(1)rate integral function: (2)rate integral derivative function: Now incorporating (37)–(39) with the flow model assumed previously, we can consequently establish new Blasingame type curves for volume fractured vertical wells in fractured tight oil reservoirs, which are shown in Figures 5–17.
Figure 5: Comparison of vertical well productivity between conventional fracturing and volume fracturing under the same pressure drop.
Figure 6: Comparison of vertical well productivity between conventional fracturing and volume fracturing at different time.
Figure 7: The effect of storage capacity ratio factor on type curves.
Figure 8: The effect of crossflow coefficient factor on type curves.
Figure 9: The effect of fracture conductivity on type curves.
Figure 10: The effect of fracture length on type curves.
Figure 11: The effect of drainage radius for volume fracturing on type curves.
Figure 12: The effect of drainage radius for conventional fracturing on type curves [25].
Figure 13: Effect of stimulated reservoir volume on vertical well productivity.
Figure 14: Effect of fracture conductivity on cumulative oil production.
Figure 15: The effect of skin factor on type curves.
Figure 16: Flow distribution along fracture surface.
Figure 17: Pressure distribution along fracture surface.
3.3.1. Comparison between Volume Fracturing and Conventional Fracturing
The comparison of vertical well productivity with the same main fracture length between conventional fracturing and volume fracturing is shown in Figures 5 and 6. From Figure 5 we can see that the productivity of volume fracturing is obviously higher than that of conventional fracturing and it has nonlinear characteristics. Due to tightness of formation, complexity of pore structure, and low permeability of matrix, the vertical well productivity is still very low after conventional fracturing. However, volume fracturing can form a fracture network in formation which greatly improves the permeability of the whole reservoir, reduces the seepage resistance from the matrix to the main fracture, increases the contact area between main fracture and formation matrix, and shortens the flow distance from the matrix to the main fracture. Therefore, in order to get higher productivity and achieve economic demand, volume fracturing must be carried out.
3.3.2. Storage Capacity Ratio and Crossflow Coefficient
The effects of storage capacity ratios () and crossflow coefficients () on the rate decline curves are given in Figures 7 and 8, respectively. From these two figures, we can see that storage capacity ratio and crossflow coefficient affect the time and degree of crossflow, respectively, similar to the common dual media. The storage capacity ratio has an influence on the production of transitional flow. The smaller the storage capacity ratio is, the more obvious the crossflow is. At the intermediate time, the curve of production is sunken. The values of both and decrease with the increase of storage capacity ratio . In pseudo steady-state flow, curves of both groups normalize, respectively. The larger the storage capacity ratio , the smaller the peak value of production. The crossflow coefficient ratio has an influence on the production of matrix-fracture crossflow. The greater the crossflow coefficient is, the earlier the crossflow happens. Figure 8 shows that if the crossflow coefficient is too small, when the pressure wave does not touch the boundary, the medium linear flow occurs and after the pressure wave reaches the boundary, pseudo steady flow will occur. Similarly, under the condition of same storage capacity ratio , the values of both and increase with the increase of crossflow coefficient in transient flow and normalize, respectively, in pseudo steady-state flow.
3.3.3. Fracture Conductivity and Fracture Half Length
The effect of artificial fracturing is to leave a high permeability channel near the well formation, which is convenient for fluid to flow from the far well zone to the bottom hole. Figures 9 and 10 show the effect of fracture conductivity and length on dimensionless rate integral and rate integral derivative for the same , , , and . From Figure 9, we find that the effect of on and is only in the early stage. The values of both and increase with the increase of fracture conductivity. However, to a certain extent, the effect of fracture length is more obvious than that of fracture conductivity on improving production effect. By combining Figures 9 and 10, we can see that, with the increase of fracture conductivity, the increase rate of production gradually reduces. When the fracture conductivity is more than 300, the effect of on and is little and can be ignored. However, with the increase of fracture length, the increase rate of production is basically changeless. The longer the fracture length is, the higher the production and the longer the stable period will be. The essence of fracture conductivity is the amount of fluid from reservoir to fracture per unit pressure gradient. In tight oil reservoirs, the fluid supply of the reservoir matrix to the main fracture is limited; that is, the conductivity is low. As the length of fracture increases, more fracture network formed by the volume fracturing is effectively connected with the main fracture, forming a larger supply area and improving the overall yield. In conventional fractured reservoirs, the main fracture conductivity is dominant, whereas in volume fractured tight oil reservoirs the effect of fracture length on productivity is dominant which is contrary to normal expectation. Therefore, in order to achieve the desired effect of volume fracturing design, a certain length of the main fracture should also be ensured.
3.3.4. Stimulated Reservoir Volume Area
Figure 11 shows the effect of drainage radius on dimensionless production, rate integral , and rate integral derivative for the same , , , and , respectively. It can be seen from Figure 11 that when dimensionless decline time is relatively smaller , the production increases with the decrease of , which is contrary to the conventional fracturing [28] (Figure 12). This is because the early flow mainly occurs between the main fracture and the microcracks. In the case of the same length of the main fracture, the smaller the drainage radius is, the higher the fracture penetration ratio and the higher the production in the early and middle period will be. When the values of time become larger , the flow is mainly between the reservoir matrix and microcracks, and the smaller the volume of fracture network, the lower the fluid supply capacity, and the faster the production decline. can also affect the time and degree of crossflow. The greater the is, the earlier the crossflow happens, and the more obvious the crossflow is.
From Figure 13 we can see that, with the drainage radius of stimulated reservoir volume and length of main fracture increasing, the area of stimulated volume, well-controlled drainage area, and well productivity also increase, but the increase rate reduces.
From Figure 14 we can see that when the stimulated reservoir volume is constant, the greater the fracture conductivity and the higher the cumulative oil production. When the length and conductivity of main fracture are constant, the contribution of stimulated reservoir volume to the cumulative oil production is not obvious. This is mainly because the single well control reserve is very limited, and there is no sense in increasing the stimulated reservoir volume without limitation. At the same time, increasing the volume of reservoir reconstruction will also increase the scale of fracturing and construction difficulty. Therefore, when the stimulated reservoir volume is constant, the length of fracture should be increased so as to improve the fracture penetration, conductivity, and well production.
3.3.5. Skin Factor
Figure 15 shows the effect of skin factor on dimensionless production, rate integral , and rate integral derivative for the same , , , and , respectively. From Figure 15 we can conclude that the production, rate integral , and rate integral derivative decrease with the increase of skin factor. When the volume fracturing effect is good, the value of skin factor will be small and even be negative. Skin factor does not affect the time and degree of crossflow.
3.3.6. Flow and Pressure Distribution of Fracture Surface
By using (34) we can obtain the flow and pressure distribution along the fracture surface. Basic parameters of the artificial main fracture and formation after volume fracturing are given in Table 2.
Table 2: Basic parameters of fracture and formation.
The flow and pressure distribution along the fracture surface are given in Figures 16 and 17. From Figures 16 and 17 we can see that the yield on both ends of the fracture is higher and the center of the fracture is lower. This is because the drainage area on both ends of the fracture is much larger than the center. The pressure distribution on the fracture is proportional to the flow distribution. The better the volume modification is, the longer the stable production period will be.
4. Conclusions
In this study, we have investigated the productivity characteristics of a volume fractured vertical well in tight oil reservoir. The specific conclusions are as follows:(1)Using Laplace transform and Stehfest numerical inversion method, a semianalytical solution for the vertical well with reconstructed fractured network is established in fractured tight oil reservoirs. The effects of an artificial main fracture near wellbore were also taken into account in this model and it can simply reflect the flow characteristics of the production wells in each stage after fracturing and acidizing treatment. Combining the calculated results with actual production data shows a good fitting performance.(2)Based on the established model, new type curves are established to analyze the flow characteristics, which can be divided into seven stages: (a) linear flow in artificial main fracture; (b) coupled boundary flow; (c) early linear flow in fractured formation; (d) mid radial flow in the microcracks of the formation; (e) mid radial flow or pseudo steady flow; (f) mid crossflow; (g) closed boundary flow.(3)Effects of some sensitive parameters such as storage capacity ratio, crossflow coefficient, fracture conductivity, fracture length, and skin factor are also analyzed in detail. Storage capacity ratios and crossflow coefficients affect the time and degree of crossflow, respectively. Artificial fracturing can leave a high permeability channel in the near well formation, which is convenient for the fluid to flow from the well zone to the bottom hole. The production increases with the increase of artificial main fracture conductivity. To a certain extent, the effect of fracture length is more obvious than that of fracture conductivity on improving production effect. When the length and conductivity of the main fracture are constant, the contribution of stimulated reservoir volume to the cumulative oil production is not obvious. And there is no sense in increasing the stimulated reservoir volume without limitation. When the stimulated reservoir volume is constant, the length of fracture should also be increased so as to improve the fracture penetration and well production.
NomenclatureDimensionless Variables: Real Domain:Dimensionless pseudo steady-state constant:Dimensionless artificial main fracture conductivity:Dimensionless well bottom pressure:Dimensionless micro fracture pressure in volume modification region:Dimensionless artificial main fracture pressure:Dimensionless decline rate:Dimensionless decline rate integral:Dimensionless decline rate integral derivative:Dimensionless time:Dimensionless decline time:Midpoint of the segment:Euler constant, 05771.Dimensionless Variables: Laplace Domain:The pressure in Laplace domain:The pressure in Laplace domain:Dimensionless micro fracture pressure in Laplace domain:Artificial main fracture pressure in volume modification region in Laplace domain:Artificial main fracture pressure in Laplace domain:The fracture rate in Laplace domain:Time variable in Laplace domain, dimensionless.Field Variables:Reservoir drainage area, m2:Compressibility for matrix, 1/Mpa:Compressibility for micro fracture, 1/Mpa:Porosity for matrix, fraction:Porosity for micro fracture:Permeability of artificial main fracture, mD:Permeability of micro fracture, mD:Permeability of matrix, mD:Formation pressure, Mpa:Initial formation pressure, Mpa:Reservoir radius, m:Equivalent drainage radius, m:Time variable, days:Fracture half length, m:Fracture width, m:Elastic storativity ratio, fraction:Crossflow coefficient, fraction.Special Functions:Modified Bessel function (2nd kind, zero order):Modified Bessel function (2nd kind, first order):Modified Bessel function (1st kind, zero order):Modified Bessel function (1st kind, first order).Special Subscripts:Dimensionless decline variable:Integral function (or initial value):Integral derivative function:Pseudo steady-state.Competing Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
This paper was supported by the Fundamental Research Funds for the Central Universities (Grant no. 53200859545), the Ministry of Land and Resources Special Geological Survey: Upper Paleozoic Marine Shale Gas Geological Survey in Yunnan, Guizhou, Guangxi Region (Grant no. DD20160178), and the Key Laboratory of Unconventional Petroleum Geology of Geological Survey of China Open Fund and Major National R&D Projects: Study on the Test Method for Shale Structure and Composition at Different Scales (Grant no. 2016ZX05034-003-006).
References C. Jia, C. Zou, J. Li, D. Li, and M. Zheng, “Assessment criteria, main types, basic features and resource prospects of the tight oil in China,” Acta Petrolei Sinica, vol. 33, no. 3, pp. 343–350, 2012. View at Google Scholar · View at ScopusC. Zou, R. Zhu, S. Wu et al., “Types, characteristics, genesis and prospects of conventional and unconventional hydrocarbon accumulations: taking tight oil and tight gas in China as an instance,” Acta Petrolei Sinica, vol. 33, no. 2, pp. 173–187, 2012. View at Google Scholar · View at ScopusJ. Nan, S. Wang, W. Yao et al., “Micro-fractures in extra-low permeability reservoir of Yanchang Formation in Ordos Basin,” Lithologic Reservoir, vol. 19, no. 4, pp. 40–44, 2007. View at Google ScholarY. Guo, J. Liu, H. Yang et al., “Hydrocarbon accumulation mechanism of low permeable tight lithologic oil reservoirs in the Yanchang Formation, Ordos Basin, China,” Petroleum Exploration and Development, vol. 39, no. 4, pp. 417–425, 2012. View at Google Scholar · View at ScopusK. M. Al-Salem, M. A. Ali, and C. Lin, “Tight oil reservoir development feasibility study using finite difference simulation and streamlines,” in Proceedings of the SPE Saudi Arabia Section Technical Symposium, Al-Khobar, Saudi Arabia, May 2009. View at Publisher · View at Google ScholarC. L. Cipolla, C. Carbo, E. P. Lolon et al., “Reservoir modeling and production evaluation in shale-gas reservoirs,” in Proceedings of the International Petroleum Technology Conference, Paper IPTC 13185, Doha, Qatar, 2009. R. Barry, “Accurate simulation of non-darcy flow in stimulated fractured shale reservoirs,” in Proceedings of the SPE Western Regional Meeting, Paper SPE 132093, Anaheim, Calif, USA, May 2010. W. Wang, G. Zhao, Y. Su et al., “Application of network fracturing technology to tight oil reservoir,” Xinjiang Petroleum Geology, vol. 34, no. 3, pp. 345–348, 2013. View at Google ScholarW.-D. Wang, Y.-L. Su, L.-J. Mu, M.-R. Tang, and L. Gao, “Influencing factors of stimulated reservoir volume of vertical wells in tight oil reservoirs,” Journal of China University of Petroleum, vol. 37, no. 3, pp. 93–97, 2013. View at Publisher · View at Google Scholar · View at ScopusH. Arvind, D. Franz, C. Martin et al., “Volumetric Fracture Modeling Approach (VFMA): incorporating microseismic data in the simulation of shale gas reservoirs,” in Proceedings of the SPE Annual Technical Conference and Exhibition, Paper SPE 134683, Florence, Italy, September 2010. C. M. Du, X. Zhang, L. Zhan et al., “Modeling hydraulic fracturing induced fracture networks in shale gas reservoirs as a dual porosity system,” in Proceedings of the International Oil and Gas Conference and Exhibition in China, Beijing, China, June 2010. View at Publisher · View at Google ScholarX. Liu and G. Zhao, “A fractal wormhole model for cold heavy oil production,” Journal of Canadian Petroleum Technology, vol. 44, no. 9, pp. 31–36, 2005. View at Google Scholar · View at ScopusX. Liu, G. Zhao, and Y. C. Jin, “Coupled reservoir/wormholes model for cold heavy oil production wells,” Journal of Petroleum Science and Engineering, vol. 50, no. 3-4, pp. 258–268, 2006. View at Publisher · View at Google Scholar · View at ScopusX. Lei and Z. Gang, “A novel approach for determining wormhole coverage in CHOPS wells,” in Proceedings of the SPE Heavy Oil Conference, Paper SPE 157935, Calgary, Canada, June 2012. X. Liu, The Research of Deliverability Evaluated for Volume Reconstruction in Tight oil Reservoir, China University of Geosciences, Beijing, China, 2013. X. Liu, C. B. Tian, and L. Y. Jiang, “Steady state deliverability evaluation model of fracture network reconstruction vertical well in tight oil reservoir,” Journal of Northeast Petroleum University, vol. 38, no. 1, pp. 91–97, 2014. View at Google ScholarX. Liao and X. Chen, “Pressure transient analysis of volume fracturing well in low permeability oil reservoir,” Science & Technology Review, vol. 34, no. 7, pp. 117–122, 2016. View at Google ScholarH. Stehfest, “Algorithm 368: numerical inversion of Laplace transforms [D5],” Communications of the ACM, vol. 13, no. 1, pp. 47–49, 1970. View at Publisher · View at Google Scholar · View at ScopusJ. E. Warren and P. J. Root, “The behavior of naturally fractured reservoirs,” Society of Petroleum Engineers Journal, vol. 3, no. 3, 1963. View at Publisher · View at Google ScholarJ. Ge, Z. Ning, Y. Liu et al., Modern Reservoir Seepage Rule, Petroleum Industry Press, 2003. F. Medeiros, B. Kurtoglu, E. Ozkan et al., “Pressure-transient performance of hydraulically horizontal well in locally and globally naturally fractured formations,” in Proceedings of the International Petroleum Technology Conference, Paper IPTC 11781-MS, Dubai, UAE, December 2007. L. Cinco, V. Samaniego, and A. Dominguez, “Transient pressure behavior for a well with a finite-conductivity vertical fracture,” Society of Petroleum Engineers Journal, vol. 18, no. 4, pp. 253–264, 1978. View at Publisher · View at Google ScholarX. D. Wang, Mechanic Basis of Fluids Flowing in Porous Media, Petroleum Industry, Beijing, China, 2006. M. F. Riley, “Analytical solutions for elliptical finite-conductivity fractures,” in Proceedings of the 66th Annual Technical Conference and E14, Paper SPE 22656, pp. 31–36, Dallas, Tex, USA, 1991. View at Publisher · View at Google ScholarL. Wang, X.-D. Wang, X.-M. Ding, L. Zhang, and C. Li, “Rate decline curves analysis of a vertical fractured well with fracture face damage,” Journal of Energy Resources Technology, vol. 134, no. 3, Article ID 032803, 9 pages, 2012. View at Publisher · View at Google Scholar · View at ScopusM. Fetkovich, “Decline curve analysis using type curves,” Journal of Petroleum Technology, vol. 32, no. 6, pp. 1065–1077, 2013. View at Publisher · View at Google ScholarT. A. Blasingame, J. L. Johnston, and W. J. Lee, “Type-curve analysis using the pressure integral method,” in Proceedings of the SPE California Regional Meeting, Bakersfield, Calif, USA, April 1989. View at ScopusT. A. Blasingame, T. L. McCray, and W. J. Lee, “Decline curve analysis for variable pressure drop/variable flowrate systems,” in Proceedings of the SPE Gas Technology Symposium, Paper no. SPE 21513, Houston, Tex, USA, January 1991. View at Publisher · View at Google Scholar