# Theoretical Study and Application of Rate Transient Analysis on Complex Fractured-Caved Carbonate Reservoirs

Authors: Qingyan Yu; Qi Wang; Pengcheng Liu; Jing Zhang; Qi Zhang; Xiaojuan Deng; Kai Feng
Journal: Geofluids
 (2021)

---

## Abstract

Carbonate reservoirs are mainly fractured-caved reservoirs with very well-developed dissolved pores, fractures, and caves. They have strong heterogeneity with various types of reservoir pore spaces. Using seismic inversion and reservoir static characterization, the result shows that the fractured-caved carbonate rocks in China are mainly caves with poor connectivity and complex oil-water distribution. Large-scale dissolved caves are mostly discrete and isolated, while the fractures are complex and various. The fracture features are observed either as a single large fractures or as a local fracture network. The characteristics of fluid flow in fracture-caved reservoirs vary as a result of the different combinations of fractures and caves. Currently, the static characterization technology of fractured-caved reservoirs is influenced by the limited resolution of seismic data, leading to large interpretation errors. In contrast, the dynamic method is a more reliable and effective method to determine reservoir parameters. However, traditional seepage equations cannot accurately characterize the flow pattern of fractured-caved carbonate reservoirs. In the case of a single large-scale dissolved fractured-caved reservoir, oil wells are usually connected to large caves through large fractures or directly drilled into large dissolved caves. In this study, the large-scale dissolved caved reservoir is simplified into two cases: (1) a single-cave and single-fracture series model composed of a single-cave and a single-fracture and (2) a composite model of dissolved caves and surrounding fracture networks. Note that the flow in a large cave is considered as free flow due to its large scale. The flow in a large fracture connected to the cave is considered as flow through porous media, and the flow in the reservoir surrounding the fracture network is considered as multiple-porosity model seepage flow. The corresponding seepage-free flow coupling mathematical model of different fractured-caved reservoirs has been established on this basis. We also obtained the rate transient analysis type curves of the oil well, conducted sensitivity analysis of each parameter, constructed the corresponding rate transient analysis curves, analyzed sensitivities of each parameter, and finally designed a dynamic evaluation method of well and reservoir parameters for different types of fractured-caved carbonate reservoirs. This study extensively applies this method in the Halahatang Oilfield of China and evaluates parameters such as reservoir reserves and physical properties to provide rational guidance for developing fractured-caved carbonate reservoirs.


---

## Body

## 1. Introduction

The fractured-caved reservoir is a vital type of carbonate reservoirs, characterized by strong heterogeneity, various types, complex connectivity, complex fluid flow mechanisms, and complicated gas-water contacts [1–3]. The reservoir is dominated by fractures, dissolved caves, and dissolved pores. Most of the dissolved caves are discretely distributed and isolated, and the fractures are highly heterogeneous, either with a single large fracture or a local fracture network [4]. Due to their very low porosity and permeability, matrix pores have minimal seepage and storage capacity. Therefore, they have little influence on the fluid flow of the entire reservoir. Generally, the influence of the matrix does not need to be taken into consideration in the fluid flow of fractured-caved reservoirs [5]. Fractures have significant permeability and are the main seepage channels. Nevertheless, their storage capacity is lower than that of large dissolved caves, so the fracture is generally considered to be a unit with low storage and high seepage and is the primary fluid flow channel [6–10]. Large-scale dissolved caves are units with high storage capacity and high permeability. However, due to the large scale of the caves, their fluid flow is considered to be free flow and no longer Darcy seepage flow [11–13]. The variable scale of fractured-caved reservoir pores determines the variable rates of fluid flow in the reservoir. The fluid flow of the reservoir presents different characteristics due to the different combinations of fractures and caves [14–16].For this type of fractured-caved reservoir, the mathematical characterization of continuum models such as dual-porosity model, triple-porosity model, multiple-porosity model, and their composite models proposed by previous researchers will no longer be applicable. The characterization method based on the discrete medium model can describe the details of the flow in different fractured caves, but this comes with modeling challenges and large computation limits when applied to well-testing and dynamic analysis theory [17, 18]. Other studies have used the pipe flow and percolation-pipe flow coupling model [19, 20] and conducted their research based on alternative theories for large fractured-caved reservoirs, but these studies were not comprehensive. Other researchers have conducted only qualitative research. They have not conducted quantitative research on reservoir parameters and oil well parameters of large fractured-caved reservoirs. As a result, there is no systematic dynamic evaluation method for reservoir parameter evaluation for the fractured-caved reservoirs [17, 18, 21]. Therefore, this study establishes a simplified model which mainly considers large dissolved caves, by simplifying an extensive dissolved caved reservoir. Specifically, there are two large models: a single cave connected to a large fracture and a composite model consisting of a single cave and surrounding reservoirs. A mathematical model of fluid flow, considering seepage-free flow couplings under different combinations, is established. The model is finally adjusted to obtain the corresponding typical curve after rate transient analysis of the oil well. Using this process and the fractured-caved parameter matching method, we can evaluate the dynamic parameter of production wells. What is more, the connectivity of fracture-caved reservoirs can be properly understood through dynamic inversion of fractured-caved reservoirs with large-scale dissolved caves and by obtaining the necessary parameters. These results are of great significance for the proper understanding and development of such fractured-caved reservoirs.


## 2. Simplified Discrete Mathematical Model and Rate Transient Analysis Type Curves

There are various types of connections between fractures and caves within fractured-caved carbonate reservoirs. Also, the relationship between dissolved caves and fractures is complex. However, most fractured-caved carbonate reservoirs are characterized by “bead-shaped” seismic reflection, generally indicating large caves (Figure1). The existence of these “bead-shaped” fractured-caved reservoirs suggests that this is one of the most common forms of connectivity [19, 20]. In cases where large-scale dissolved caves have been developed in carbonate reservoirs, the caves are connected to the wellbore through large fractures, or vertical wells are drilled directly on large caves; the capacity and permeability of the matrix cannot be taken into consideration in the reservoir and the reservoir can be described mathematically utilizing a simplified discrete model with seepage-free flow coupling. The following assumptions can be made for the simplified discrete large-scale fracture-cave model: (1) oil wells produce at a stable rate; (2) oil reservoir is a depletion drive reservoir; (3) the fluid is single phase and weakly compressible, and its compression coefficient and volume coefficient are constant; (4) large caves are weakly compressible, and their compressibility is constant; (5) large caves are unfilled or semifilled, as spheres or cylinders with radius R, and flow is considered as free flow; (6) the large fracture is flat and the fluid flow in the fracture follows Darcy’s law, and the fluid flow in the surrounding multimedium reservoir also follows Darcy’s law; and (7) gravity acceleration, wellbore storage effect, and skin effect are ignored.Figure 1
Seismic characteristics, inversion results, and development models of fractured-caved carbonate reservoirs: (a) seismic reflection characteristics of fractured-caved reservoirs; (b) seismic inversion results of fracture-cavity reservoirs; (c) fractured-caved reservoir development schematic diagram.
(a)(b)(c)


### 2.1. Single-Cave and Single-Fracture Series Model

When oil wells are drilled through a large fracture, and the fracture is connected to a large-scale cave (i.e., the large fracture and the cave form a single-cave and single-fracture series model), fluid production ultimately depends on the elastic energy of the large cave, the fracture, and the fluid. Large caves are the main reservoir space, and fractures are the effective secondary reservoir space and the main seepage channel. Large caves cannot supply fluid directly to the wellbore but can indirectly supply fluid through fractures.


#### 2.1.1. Physical Model

Figure2 shows the single-cave and single-fracture series physical model. Where the cave’s central point is the origin O, the intersection point of the cave, and the fracture is x1, and the connection point of the fracture and the wellbore is x2. It can be seen that the cave radius is R=x1, the fracture length is L1=x2−x1, the fracture height is L2, the fracture length is W, and the wellbore radius is rw. The oil well is directly connected to the fractures, and the dissolved caves produce fluid through the fractures. Figure 3 is a schematic diagram of the seepage flow of a single-cave and single-fracture series connection model.Figure 2
Schematic diagram of single-cave and single-fracture series physical model.Figure 3
Schematic diagram of single-cave and single-fracture series seepage flow model.


#### 2.1.2. Establishment and Solution of the Mathematical Model

Large caves are the main reservoir spaces, and fractures are the effective reservoir and seepage spaces. The fluid flow capacity in a large cave is robust, and it is considered as free flow, i.e., the fluid flowing in caves results in very little pressure drop. Then, the dimensionless partial differential equation describing fluid flow in the reservoir is:(1)∂2pFD∂xD2=ωF∂pFD∂tD,x1D≤xD≤x2D,∂pFD∂xDxD=x1D=ωVRD3L2DWDdpVDdtD,∂pFD∂xDxD=x2D=1L2DWD,pFDxD=x1D=pVD,pFDxD=x2D=pwD,pFDtD=0=pVDtD=0=pwDtD=0=0.In the above equations, dimensionless quantities are defined as follows:(2)xD=xrw,x1,2D=x1,2rw,pF,V,wD=86.4kFrwqμBpi−pF,V,w,ωF,V=ϕF,VCtF,tVϕVCtV+ϕFCtF,tD=3.6kFtμϕVCtV+ϕFCtFrw2,RD=π181/3Rrw,L2D=L2rw,WD=WrW.Employing the Laplace transform, the dimensionless Laplace solution can be obtained:(3)p~wD=a3a1+a3a4a2a4a2−a1,where a1=wFs−wVRD3/L2DWDs, a2=wFs+wVRD3/L2DWDs, a3=1/L2DWDsωFs, and a4=e−2wFsx1D+2wFsx2D.Using the Duhamel principle, the Laplace solution of the dimensionless bottomhole rate is(4)q~Ds=1s2p~wDs.


#### 2.1.3. Analysis of the Rate Transient Analysis Typical Curves and Influencing Factors

Numerical inversion is performed for the above-mentioned Laplace solutions. The logarithm of dimensionless productionqD, dimensionless rate integral qDi, and dimensionless rate integral derivative qDid are considered as the ordinate, and the logarithm of dimensionless time tD is taken as the abscissa to draw the rate transient analysis curves of the oil well, as shown in Figure 4.Figure 4
Vertical well rate transient analysis curve of single-cave and single-fracture model.It can be seen from Figure4 that the oil well rate transient analysis curves of the single-cave and single-fracture model can be divided into three stages. (I) The fracture linear flow stage, which is before the pressure wave reaches the cave, is mainly affected by the fracture’s characteristic parameters. The dimensionless rate integral curve is parallel to the rate integral derivative curve, and the slope is -1/2. (II) The cave’s response stage. This stage is mainly affected by the dissolved cave properties. The dimensionless integral curve tends to be horizontal, and the dimensionless integral derivative curve tends to “concave.” The concave show the fluid supplementation from dissolved caves to fractures. (III) The boundary control flow stage. The dimensionless rate integral curve and the rate integral derivative curve overlap, showing a straight line with a slope of -1.The sensitivity analysis of the factors affecting the vertical well rate transient analysis curves with the single-cave and single-fracture model is as follows. Figure5(a) shows the influence of dimensionless caves radius RD on the rate transient analysis curves of the oil well. With the larger RD, the “concave” will be deeper and wider, and the cave’s response stage will last longer. In addition, the boundary response stage will be later, i.e., the cave’s storage capacity will be more significant, and the stable production capacity of the well will be better. Figure 5(b) shows the influence of fracture length L1D on the rate transient analysis curves of the oil well. L1D affects the duration of the fracture linear flow phase. With the longer L1D, the parallel sections of the dimensionless rate integral and rate integral derivative curves are longer, the fracture linear flow stage lasts longer, and the cave response stage starts later. Figures 5(c) and 5(d) show the influence of fracture width WD and fracture height L2D on the well rate transient analysis curve. WD, L2D, and their products mainly affect fracture permeability and storage capacity. With a higher value, the fracture seepage performance is better, and the fluid will consume less energy. When the dimensionless curve in the fracture linear flow stage trends upward, the product of WD and L2D becomes more significant, which weakens the influence of caves, and then the corresponding concave is relatively shallower and narrower. Figure 5(e) shows the influence of fracture storage capacity ratio ωf on the rate transient analysis curve. With the increase of ωF, the fracture storage capacity will be better, and the rate transient analysis curve closer to upward. And as the corresponding ωV decreases, the concave becomes shallower and narrower.Figure 5
Vertical well rate transient analysis curves of single-cave and single-fracture model: (a) the influence ofRD on the rate transient analysis curve; (b) the influence of L1D on rate transient analysis curve; (c) the influence of WD on rate transient analysis curve; (d) the influence of L2D on rate transient analysis curve; (e) the influence of ωF on rate transient analysis curve.
(a)(b)(c)(d)(e)


#### 2.1.4. Rate Transient Analysis Typical Curve of Simplified Model

If the fracture lengthL1D=0, the single-cave and single-fracture model can be simplified to the single-cave model. In this case, the oil well is directly drilled in the large-scale cave, and the storage capacity and permeability of the surrounding matrix is negligible. The single-cave physical model is shown in Figure 6, and the single-cave seepage model is shown in Figure 7.Figure 6
Single-cave physical model.Figure 7
Single-cave seepage model.Figure8 is the vertical well rate transient analysis curve of the single-cave model. We found that the oil well rate transient analysis curve of the single-cave model can be divided into two stages. (I) The cave response stage. In the cave response stage, the dimensionless rate integral curve and the rate integral derivative curve overlap and form a line with a slope of -1. The wellbore is then directly drilled in the dissolved cave, which is equivalent to an enlarged wellbore and consistent with the wellbore storage effect. (II) The boundary control flow stage. For the free flow of large isolated caves, the pressure wave will quickly propagate to the boundary. The dimensionless rate integral curve coincides with the rate integral derivative curve and is also a straight line with a slope of -1.Figure 8
Vertical well rate transient analysis curve of single-cave model.


### 2.2. Single-Cave and Multiporosity Composite Model

When large-scale caves are developed in carbonate reservoirs and small-scale dissolved pores and fracture systems are evenly distributed around the caves, then the permeability and storage capacity of the matrix pores and fracture network around the caves must be taken into consideration. In this case, the simplified large-scale cave and multiple-porosity composite model can be used for reservoir characterization.


#### 2.2.1. Physical Model

Figure9 shows the composite physical model of single cave and multiple porosity. The oil well is directly drilled in large-scale dissolved caves. The area around the dissolved caves can be modeled as a single medium, double medium, or multiple medium model according to the characteristics of the small-scale pores, vugs and fractures. The fractured-caved reservoir is considered as a three-dimensional cylindrical composite model, in which the large-scale dissolved caves are internal cylinders. The cylinder height is h, the radius is RV, and the outside is multimedium and concentric with the cylinder. The fluid flow in large-scale dissolved cave is considered to be free flow, while that in the surrounded multiple-porosity model is considered to be seepage flow. Figure 10 is a fluid flow schematic. The surrounded improved matrix forms because of the development of small-scale caves.Figure 9
Single-cave and multimedium composite model.Figure 10
Single-cave and multimedium composite model.


#### 2.2.2. Establishment and Solution of the Mathematical Model

The fluid flow capacity in large caves is extraordinarily strong and therefore considered to be free flow, i.e., the fluid flow in the caves produces truly little pressure drawdown. The dimensionless partial differential equation describing the composite model composed of large-scale caves and surrounding formation is as follows (considering improved matrix that develops small-scale caves):(5)∂2pfD∂rD2+1rD∂pfD∂rD+λimpimD−pfD=ωf∂pfD∂tD,ωim∂pimD∂tD+λimpimD−pfD=0,1+RD∂pfD∂rDrD=RD=RD2ωV∂pfD∂tDrD=RD,pfDrD=RD=pwD=pvD,∂pfD∂rDrD=reD=0,pimDt=0=pfDt=0=pVDt=0=0.In the above equations, dimensionless quantities are defined as follows:(6)pf,im,VDrD,tD=kfh1.842×10−3Bqμpi−pf,im,Vr,t,rD=rrw,reD=rerw,RD=Rrw,tD=3.6kftϕfCtf+ϕimCtimμrw2,ωf,im,V=ϕf,im,VCtf,tim,tVϕfCtf+ϕimCtim+ϕVCtV,λim,v=αimkimkfrw2.The dimensionless Laplace space solution expression can be written as Equation (7) by combining the above equations and the Laplace transformation.
(7)p~wD=Y1I0RDσ+K0RDσsY3−Y1Y2,where
(8)Y1=K1reDσI1reDσ,Y2=RDσI1RDσ−RD2ωVsI0RDσ,Y3=RDσK1RDσ+RD2ωVsK0RDσ,σ=sfs=sωf+λimωimωims+λim.Using the Duhamel principle, the Laplace solution of the dimensionless bottomhole rate can be expressed as(9)q~Ds=1s2p~wDs.


#### 2.2.3. Analysis of Rate Transient Analysis Typical Curves and the Influencing Factors

Similarly, we perform a numerical inversion on the Laplace solution above. The dimensionless oil well rate transient analysis curves form, as shown in Figure11.Figure 11
Oil well rate transient analysis curves of single-cave and multiple-porosity composite model.As shown in Figure11, the characteristic curve of the composite model of a large-scale cave and multiple-porosity model can be divided into six stages. (I) The large-scale dissolved cave response stage is equivalent to an enlarged “wellbore,” showing the wellbore storage effect. The dimensionless rate integral and rate integral derivative curves overlap into a straight line with a slope of -1: (II) the fractured peripheral multiple-porosity model linear flow and radial flow response stages, mainly related to the fracture properties of the multiple-porosity model; (III) the fluid flow from improved matrix to fractures in the surrounding multiple media, the dimensionless rate integral derivative curve firstly appears “concave”; (IV) the radial flow stage of the model; and (V) the boundary control flow stage. The rate integral curve and the rate integral derivative curve overlap in a line with a slope of -1 for a closed boundary.The sensitivity analysis of the influencing factors on the oil well rate transient analysis curve of the single-cave and dual-porosity composite model is carried out as shown in Figure12.Figure 12
The influence of different parameters on the rate transient analysis curves: (a) the influence ofRD on the rate transient analysis curves; (b) the influence of ωV on the rate transient analysis curves; (c) the influence of reD on the rate transient analysis curves.
(a)(b)(c)As shown in Figure12(a), with the larger radius RD of the large central cave, the intersecting line of the dimensionless rate integral and the rate integral derivative curve is close to upward. The stage’s duration reflects the fact that the wellbore storage effect is longer. Figure 12(b) shows that the effect of ωV on the curve is similar to that of RD, but the effect of ωV is relatively small.IfωV=0 and ωf=0.8, the single-cave and dual-porosity composite model can be simplified to a single-cave and single-medium composite model. Figure 13 shows the oil well rate transient analysis curves with a single cave and single medium composite model.Figure 13
Oil well rate transient analysis curve with single-cave and single-medium composite model.


## 2.1. Single-Cave and Single-Fracture Series Model

When oil wells are drilled through a large fracture, and the fracture is connected to a large-scale cave (i.e., the large fracture and the cave form a single-cave and single-fracture series model), fluid production ultimately depends on the elastic energy of the large cave, the fracture, and the fluid. Large caves are the main reservoir space, and fractures are the effective secondary reservoir space and the main seepage channel. Large caves cannot supply fluid directly to the wellbore but can indirectly supply fluid through fractures.


### 2.1.1. Physical Model

Figure2 shows the single-cave and single-fracture series physical model. Where the cave’s central point is the origin O, the intersection point of the cave, and the fracture is x1, and the connection point of the fracture and the wellbore is x2. It can be seen that the cave radius is R=x1, the fracture length is L1=x2−x1, the fracture height is L2, the fracture length is W, and the wellbore radius is rw. The oil well is directly connected to the fractures, and the dissolved caves produce fluid through the fractures. Figure 3 is a schematic diagram of the seepage flow of a single-cave and single-fracture series connection model.Figure 2
Schematic diagram of single-cave and single-fracture series physical model.Figure 3
Schematic diagram of single-cave and single-fracture series seepage flow model.


### 2.1.2. Establishment and Solution of the Mathematical Model

Large caves are the main reservoir spaces, and fractures are the effective reservoir and seepage spaces. The fluid flow capacity in a large cave is robust, and it is considered as free flow, i.e., the fluid flowing in caves results in very little pressure drop. Then, the dimensionless partial differential equation describing fluid flow in the reservoir is:(1)∂2pFD∂xD2=ωF∂pFD∂tD,x1D≤xD≤x2D,∂pFD∂xDxD=x1D=ωVRD3L2DWDdpVDdtD,∂pFD∂xDxD=x2D=1L2DWD,pFDxD=x1D=pVD,pFDxD=x2D=pwD,pFDtD=0=pVDtD=0=pwDtD=0=0.In the above equations, dimensionless quantities are defined as follows:(2)xD=xrw,x1,2D=x1,2rw,pF,V,wD=86.4kFrwqμBpi−pF,V,w,ωF,V=ϕF,VCtF,tVϕVCtV+ϕFCtF,tD=3.6kFtμϕVCtV+ϕFCtFrw2,RD=π181/3Rrw,L2D=L2rw,WD=WrW.Employing the Laplace transform, the dimensionless Laplace solution can be obtained:(3)p~wD=a3a1+a3a4a2a4a2−a1,where a1=wFs−wVRD3/L2DWDs, a2=wFs+wVRD3/L2DWDs, a3=1/L2DWDsωFs, and a4=e−2wFsx1D+2wFsx2D.Using the Duhamel principle, the Laplace solution of the dimensionless bottomhole rate is(4)q~Ds=1s2p~wDs.


### 2.1.3. Analysis of the Rate Transient Analysis Typical Curves and Influencing Factors

Numerical inversion is performed for the above-mentioned Laplace solutions. The logarithm of dimensionless productionqD, dimensionless rate integral qDi, and dimensionless rate integral derivative qDid are considered as the ordinate, and the logarithm of dimensionless time tD is taken as the abscissa to draw the rate transient analysis curves of the oil well, as shown in Figure 4.Figure 4
Vertical well rate transient analysis curve of single-cave and single-fracture model.It can be seen from Figure4 that the oil well rate transient analysis curves of the single-cave and single-fracture model can be divided into three stages. (I) The fracture linear flow stage, which is before the pressure wave reaches the cave, is mainly affected by the fracture’s characteristic parameters. The dimensionless rate integral curve is parallel to the rate integral derivative curve, and the slope is -1/2. (II) The cave’s response stage. This stage is mainly affected by the dissolved cave properties. The dimensionless integral curve tends to be horizontal, and the dimensionless integral derivative curve tends to “concave.” The concave show the fluid supplementation from dissolved caves to fractures. (III) The boundary control flow stage. The dimensionless rate integral curve and the rate integral derivative curve overlap, showing a straight line with a slope of -1.The sensitivity analysis of the factors affecting the vertical well rate transient analysis curves with the single-cave and single-fracture model is as follows. Figure5(a) shows the influence of dimensionless caves radius RD on the rate transient analysis curves of the oil well. With the larger RD, the “concave” will be deeper and wider, and the cave’s response stage will last longer. In addition, the boundary response stage will be later, i.e., the cave’s storage capacity will be more significant, and the stable production capacity of the well will be better. Figure 5(b) shows the influence of fracture length L1D on the rate transient analysis curves of the oil well. L1D affects the duration of the fracture linear flow phase. With the longer L1D, the parallel sections of the dimensionless rate integral and rate integral derivative curves are longer, the fracture linear flow stage lasts longer, and the cave response stage starts later. Figures 5(c) and 5(d) show the influence of fracture width WD and fracture height L2D on the well rate transient analysis curve. WD, L2D, and their products mainly affect fracture permeability and storage capacity. With a higher value, the fracture seepage performance is better, and the fluid will consume less energy. When the dimensionless curve in the fracture linear flow stage trends upward, the product of WD and L2D becomes more significant, which weakens the influence of caves, and then the corresponding concave is relatively shallower and narrower. Figure 5(e) shows the influence of fracture storage capacity ratio ωf on the rate transient analysis curve. With the increase of ωF, the fracture storage capacity will be better, and the rate transient analysis curve closer to upward. And as the corresponding ωV decreases, the concave becomes shallower and narrower.Figure 5
Vertical well rate transient analysis curves of single-cave and single-fracture model: (a) the influence ofRD on the rate transient analysis curve; (b) the influence of L1D on rate transient analysis curve; (c) the influence of WD on rate transient analysis curve; (d) the influence of L2D on rate transient analysis curve; (e) the influence of ωF on rate transient analysis curve.
(a)(b)(c)(d)(e)


### 2.1.4. Rate Transient Analysis Typical Curve of Simplified Model

If the fracture lengthL1D=0, the single-cave and single-fracture model can be simplified to the single-cave model. In this case, the oil well is directly drilled in the large-scale cave, and the storage capacity and permeability of the surrounding matrix is negligible. The single-cave physical model is shown in Figure 6, and the single-cave seepage model is shown in Figure 7.Figure 6
Single-cave physical model.Figure 7
Single-cave seepage model.Figure8 is the vertical well rate transient analysis curve of the single-cave model. We found that the oil well rate transient analysis curve of the single-cave model can be divided into two stages. (I) The cave response stage. In the cave response stage, the dimensionless rate integral curve and the rate integral derivative curve overlap and form a line with a slope of -1. The wellbore is then directly drilled in the dissolved cave, which is equivalent to an enlarged wellbore and consistent with the wellbore storage effect. (II) The boundary control flow stage. For the free flow of large isolated caves, the pressure wave will quickly propagate to the boundary. The dimensionless rate integral curve coincides with the rate integral derivative curve and is also a straight line with a slope of -1.Figure 8
Vertical well rate transient analysis curve of single-cave model.


## 2.1.1. Physical Model

Figure2 shows the single-cave and single-fracture series physical model. Where the cave’s central point is the origin O, the intersection point of the cave, and the fracture is x1, and the connection point of the fracture and the wellbore is x2. It can be seen that the cave radius is R=x1, the fracture length is L1=x2−x1, the fracture height is L2, the fracture length is W, and the wellbore radius is rw. The oil well is directly connected to the fractures, and the dissolved caves produce fluid through the fractures. Figure 3 is a schematic diagram of the seepage flow of a single-cave and single-fracture series connection model.Figure 2
Schematic diagram of single-cave and single-fracture series physical model.Figure 3
Schematic diagram of single-cave and single-fracture series seepage flow model.


## 2.1.2. Establishment and Solution of the Mathematical Model

Large caves are the main reservoir spaces, and fractures are the effective reservoir and seepage spaces. The fluid flow capacity in a large cave is robust, and it is considered as free flow, i.e., the fluid flowing in caves results in very little pressure drop. Then, the dimensionless partial differential equation describing fluid flow in the reservoir is:(1)∂2pFD∂xD2=ωF∂pFD∂tD,x1D≤xD≤x2D,∂pFD∂xDxD=x1D=ωVRD3L2DWDdpVDdtD,∂pFD∂xDxD=x2D=1L2DWD,pFDxD=x1D=pVD,pFDxD=x2D=pwD,pFDtD=0=pVDtD=0=pwDtD=0=0.In the above equations, dimensionless quantities are defined as follows:(2)xD=xrw,x1,2D=x1,2rw,pF,V,wD=86.4kFrwqμBpi−pF,V,w,ωF,V=ϕF,VCtF,tVϕVCtV+ϕFCtF,tD=3.6kFtμϕVCtV+ϕFCtFrw2,RD=π181/3Rrw,L2D=L2rw,WD=WrW.Employing the Laplace transform, the dimensionless Laplace solution can be obtained:(3)p~wD=a3a1+a3a4a2a4a2−a1,where a1=wFs−wVRD3/L2DWDs, a2=wFs+wVRD3/L2DWDs, a3=1/L2DWDsωFs, and a4=e−2wFsx1D+2wFsx2D.Using the Duhamel principle, the Laplace solution of the dimensionless bottomhole rate is(4)q~Ds=1s2p~wDs.


## 2.1.3. Analysis of the Rate Transient Analysis Typical Curves and Influencing Factors

Numerical inversion is performed for the above-mentioned Laplace solutions. The logarithm of dimensionless productionqD, dimensionless rate integral qDi, and dimensionless rate integral derivative qDid are considered as the ordinate, and the logarithm of dimensionless time tD is taken as the abscissa to draw the rate transient analysis curves of the oil well, as shown in Figure 4.Figure 4
Vertical well rate transient analysis curve of single-cave and single-fracture model.It can be seen from Figure4 that the oil well rate transient analysis curves of the single-cave and single-fracture model can be divided into three stages. (I) The fracture linear flow stage, which is before the pressure wave reaches the cave, is mainly affected by the fracture’s characteristic parameters. The dimensionless rate integral curve is parallel to the rate integral derivative curve, and the slope is -1/2. (II) The cave’s response stage. This stage is mainly affected by the dissolved cave properties. The dimensionless integral curve tends to be horizontal, and the dimensionless integral derivative curve tends to “concave.” The concave show the fluid supplementation from dissolved caves to fractures. (III) The boundary control flow stage. The dimensionless rate integral curve and the rate integral derivative curve overlap, showing a straight line with a slope of -1.The sensitivity analysis of the factors affecting the vertical well rate transient analysis curves with the single-cave and single-fracture model is as follows. Figure5(a) shows the influence of dimensionless caves radius RD on the rate transient analysis curves of the oil well. With the larger RD, the “concave” will be deeper and wider, and the cave’s response stage will last longer. In addition, the boundary response stage will be later, i.e., the cave’s storage capacity will be more significant, and the stable production capacity of the well will be better. Figure 5(b) shows the influence of fracture length L1D on the rate transient analysis curves of the oil well. L1D affects the duration of the fracture linear flow phase. With the longer L1D, the parallel sections of the dimensionless rate integral and rate integral derivative curves are longer, the fracture linear flow stage lasts longer, and the cave response stage starts later. Figures 5(c) and 5(d) show the influence of fracture width WD and fracture height L2D on the well rate transient analysis curve. WD, L2D, and their products mainly affect fracture permeability and storage capacity. With a higher value, the fracture seepage performance is better, and the fluid will consume less energy. When the dimensionless curve in the fracture linear flow stage trends upward, the product of WD and L2D becomes more significant, which weakens the influence of caves, and then the corresponding concave is relatively shallower and narrower. Figure 5(e) shows the influence of fracture storage capacity ratio ωf on the rate transient analysis curve. With the increase of ωF, the fracture storage capacity will be better, and the rate transient analysis curve closer to upward. And as the corresponding ωV decreases, the concave becomes shallower and narrower.Figure 5
Vertical well rate transient analysis curves of single-cave and single-fracture model: (a) the influence ofRD on the rate transient analysis curve; (b) the influence of L1D on rate transient analysis curve; (c) the influence of WD on rate transient analysis curve; (d) the influence of L2D on rate transient analysis curve; (e) the influence of ωF on rate transient analysis curve.
(a)(b)(c)(d)(e)


## 2.1.4. Rate Transient Analysis Typical Curve of Simplified Model

If the fracture lengthL1D=0, the single-cave and single-fracture model can be simplified to the single-cave model. In this case, the oil well is directly drilled in the large-scale cave, and the storage capacity and permeability of the surrounding matrix is negligible. The single-cave physical model is shown in Figure 6, and the single-cave seepage model is shown in Figure 7.Figure 6
Single-cave physical model.Figure 7
Single-cave seepage model.Figure8 is the vertical well rate transient analysis curve of the single-cave model. We found that the oil well rate transient analysis curve of the single-cave model can be divided into two stages. (I) The cave response stage. In the cave response stage, the dimensionless rate integral curve and the rate integral derivative curve overlap and form a line with a slope of -1. The wellbore is then directly drilled in the dissolved cave, which is equivalent to an enlarged wellbore and consistent with the wellbore storage effect. (II) The boundary control flow stage. For the free flow of large isolated caves, the pressure wave will quickly propagate to the boundary. The dimensionless rate integral curve coincides with the rate integral derivative curve and is also a straight line with a slope of -1.Figure 8
Vertical well rate transient analysis curve of single-cave model.


## 2.2. Single-Cave and Multiporosity Composite Model

When large-scale caves are developed in carbonate reservoirs and small-scale dissolved pores and fracture systems are evenly distributed around the caves, then the permeability and storage capacity of the matrix pores and fracture network around the caves must be taken into consideration. In this case, the simplified large-scale cave and multiple-porosity composite model can be used for reservoir characterization.


### 2.2.1. Physical Model

Figure9 shows the composite physical model of single cave and multiple porosity. The oil well is directly drilled in large-scale dissolved caves. The area around the dissolved caves can be modeled as a single medium, double medium, or multiple medium model according to the characteristics of the small-scale pores, vugs and fractures. The fractured-caved reservoir is considered as a three-dimensional cylindrical composite model, in which the large-scale dissolved caves are internal cylinders. The cylinder height is h, the radius is RV, and the outside is multimedium and concentric with the cylinder. The fluid flow in large-scale dissolved cave is considered to be free flow, while that in the surrounded multiple-porosity model is considered to be seepage flow. Figure 10 is a fluid flow schematic. The surrounded improved matrix forms because of the development of small-scale caves.Figure 9
Single-cave and multimedium composite model.Figure 10
Single-cave and multimedium composite model.


### 2.2.2. Establishment and Solution of the Mathematical Model

The fluid flow capacity in large caves is extraordinarily strong and therefore considered to be free flow, i.e., the fluid flow in the caves produces truly little pressure drawdown. The dimensionless partial differential equation describing the composite model composed of large-scale caves and surrounding formation is as follows (considering improved matrix that develops small-scale caves):(5)∂2pfD∂rD2+1rD∂pfD∂rD+λimpimD−pfD=ωf∂pfD∂tD,ωim∂pimD∂tD+λimpimD−pfD=0,1+RD∂pfD∂rDrD=RD=RD2ωV∂pfD∂tDrD=RD,pfDrD=RD=pwD=pvD,∂pfD∂rDrD=reD=0,pimDt=0=pfDt=0=pVDt=0=0.In the above equations, dimensionless quantities are defined as follows:(6)pf,im,VDrD,tD=kfh1.842×10−3Bqμpi−pf,im,Vr,t,rD=rrw,reD=rerw,RD=Rrw,tD=3.6kftϕfCtf+ϕimCtimμrw2,ωf,im,V=ϕf,im,VCtf,tim,tVϕfCtf+ϕimCtim+ϕVCtV,λim,v=αimkimkfrw2.The dimensionless Laplace space solution expression can be written as Equation (7) by combining the above equations and the Laplace transformation.
(7)p~wD=Y1I0RDσ+K0RDσsY3−Y1Y2,where
(8)Y1=K1reDσI1reDσ,Y2=RDσI1RDσ−RD2ωVsI0RDσ,Y3=RDσK1RDσ+RD2ωVsK0RDσ,σ=sfs=sωf+λimωimωims+λim.Using the Duhamel principle, the Laplace solution of the dimensionless bottomhole rate can be expressed as(9)q~Ds=1s2p~wDs.


### 2.2.3. Analysis of Rate Transient Analysis Typical Curves and the Influencing Factors

Similarly, we perform a numerical inversion on the Laplace solution above. The dimensionless oil well rate transient analysis curves form, as shown in Figure11.Figure 11
Oil well rate transient analysis curves of single-cave and multiple-porosity composite model.As shown in Figure11, the characteristic curve of the composite model of a large-scale cave and multiple-porosity model can be divided into six stages. (I) The large-scale dissolved cave response stage is equivalent to an enlarged “wellbore,” showing the wellbore storage effect. The dimensionless rate integral and rate integral derivative curves overlap into a straight line with a slope of -1: (II) the fractured peripheral multiple-porosity model linear flow and radial flow response stages, mainly related to the fracture properties of the multiple-porosity model; (III) the fluid flow from improved matrix to fractures in the surrounding multiple media, the dimensionless rate integral derivative curve firstly appears “concave”; (IV) the radial flow stage of the model; and (V) the boundary control flow stage. The rate integral curve and the rate integral derivative curve overlap in a line with a slope of -1 for a closed boundary.The sensitivity analysis of the influencing factors on the oil well rate transient analysis curve of the single-cave and dual-porosity composite model is carried out as shown in Figure12.Figure 12
The influence of different parameters on the rate transient analysis curves: (a) the influence ofRD on the rate transient analysis curves; (b) the influence of ωV on the rate transient analysis curves; (c) the influence of reD on the rate transient analysis curves.
(a)(b)(c)As shown in Figure12(a), with the larger radius RD of the large central cave, the intersecting line of the dimensionless rate integral and the rate integral derivative curve is close to upward. The stage’s duration reflects the fact that the wellbore storage effect is longer. Figure 12(b) shows that the effect of ωV on the curve is similar to that of RD, but the effect of ωV is relatively small.IfωV=0 and ωf=0.8, the single-cave and dual-porosity composite model can be simplified to a single-cave and single-medium composite model. Figure 13 shows the oil well rate transient analysis curves with a single cave and single medium composite model.Figure 13
Oil well rate transient analysis curve with single-cave and single-medium composite model.


## 2.2.1. Physical Model

Figure9 shows the composite physical model of single cave and multiple porosity. The oil well is directly drilled in large-scale dissolved caves. The area around the dissolved caves can be modeled as a single medium, double medium, or multiple medium model according to the characteristics of the small-scale pores, vugs and fractures. The fractured-caved reservoir is considered as a three-dimensional cylindrical composite model, in which the large-scale dissolved caves are internal cylinders. The cylinder height is h, the radius is RV, and the outside is multimedium and concentric with the cylinder. The fluid flow in large-scale dissolved cave is considered to be free flow, while that in the surrounded multiple-porosity model is considered to be seepage flow. Figure 10 is a fluid flow schematic. The surrounded improved matrix forms because of the development of small-scale caves.Figure 9
Single-cave and multimedium composite model.Figure 10
Single-cave and multimedium composite model.


## 2.2.2. Establishment and Solution of the Mathematical Model

The fluid flow capacity in large caves is extraordinarily strong and therefore considered to be free flow, i.e., the fluid flow in the caves produces truly little pressure drawdown. The dimensionless partial differential equation describing the composite model composed of large-scale caves and surrounding formation is as follows (considering improved matrix that develops small-scale caves):(5)∂2pfD∂rD2+1rD∂pfD∂rD+λimpimD−pfD=ωf∂pfD∂tD,ωim∂pimD∂tD+λimpimD−pfD=0,1+RD∂pfD∂rDrD=RD=RD2ωV∂pfD∂tDrD=RD,pfDrD=RD=pwD=pvD,∂pfD∂rDrD=reD=0,pimDt=0=pfDt=0=pVDt=0=0.In the above equations, dimensionless quantities are defined as follows:(6)pf,im,VDrD,tD=kfh1.842×10−3Bqμpi−pf,im,Vr,t,rD=rrw,reD=rerw,RD=Rrw,tD=3.6kftϕfCtf+ϕimCtimμrw2,ωf,im,V=ϕf,im,VCtf,tim,tVϕfCtf+ϕimCtim+ϕVCtV,λim,v=αimkimkfrw2.The dimensionless Laplace space solution expression can be written as Equation (7) by combining the above equations and the Laplace transformation.
(7)p~wD=Y1I0RDσ+K0RDσsY3−Y1Y2,where
(8)Y1=K1reDσI1reDσ,Y2=RDσI1RDσ−RD2ωVsI0RDσ,Y3=RDσK1RDσ+RD2ωVsK0RDσ,σ=sfs=sωf+λimωimωims+λim.Using the Duhamel principle, the Laplace solution of the dimensionless bottomhole rate can be expressed as(9)q~Ds=1s2p~wDs.


## 2.2.3. Analysis of Rate Transient Analysis Typical Curves and the Influencing Factors

Similarly, we perform a numerical inversion on the Laplace solution above. The dimensionless oil well rate transient analysis curves form, as shown in Figure11.Figure 11
Oil well rate transient analysis curves of single-cave and multiple-porosity composite model.As shown in Figure11, the characteristic curve of the composite model of a large-scale cave and multiple-porosity model can be divided into six stages. (I) The large-scale dissolved cave response stage is equivalent to an enlarged “wellbore,” showing the wellbore storage effect. The dimensionless rate integral and rate integral derivative curves overlap into a straight line with a slope of -1: (II) the fractured peripheral multiple-porosity model linear flow and radial flow response stages, mainly related to the fracture properties of the multiple-porosity model; (III) the fluid flow from improved matrix to fractures in the surrounding multiple media, the dimensionless rate integral derivative curve firstly appears “concave”; (IV) the radial flow stage of the model; and (V) the boundary control flow stage. The rate integral curve and the rate integral derivative curve overlap in a line with a slope of -1 for a closed boundary.The sensitivity analysis of the influencing factors on the oil well rate transient analysis curve of the single-cave and dual-porosity composite model is carried out as shown in Figure12.Figure 12
The influence of different parameters on the rate transient analysis curves: (a) the influence ofRD on the rate transient analysis curves; (b) the influence of ωV on the rate transient analysis curves; (c) the influence of reD on the rate transient analysis curves.
(a)(b)(c)As shown in Figure12(a), with the larger radius RD of the large central cave, the intersecting line of the dimensionless rate integral and the rate integral derivative curve is close to upward. The stage’s duration reflects the fact that the wellbore storage effect is longer. Figure 12(b) shows that the effect of ωV on the curve is similar to that of RD, but the effect of ωV is relatively small.IfωV=0 and ωf=0.8, the single-cave and dual-porosity composite model can be simplified to a single-cave and single-medium composite model. Figure 13 shows the oil well rate transient analysis curves with a single cave and single medium composite model.Figure 13
Oil well rate transient analysis curve with single-cave and single-medium composite model.


## 3. Application of Type Curves and the Evaluation Processes of Dynamic Parameters

### 3.1. Dynamic Evaluation Method and Processes

For the above discretized simplified models of different fractured-caved reservoirs developed in a large-scale single cave, the specific evaluation method and dynamic matching process of oil well and reservoir parameters are as follows:(1)
Obtain the oil well actual production data and calculate the normalized production rate, rate integral, and rate integral derivative curves [22, 23]. The specific calculation method adopts Equations (10)–(13)The equation for calculating the material balance pseudo-time of the oil well:(10)td=Npq.The equation for calculating the normalized production rate of the oil well:(11)qd=qΔp=qpi−pwf.The equation for calculating the integral of normalized production rate for an oil well:(12)qdi=1td∫0tdqddτ.The equation for calculating the integral derivative of normalized oil well production rate:(13)qdid=tddqdidtd.(2)
In accordance with the initial static reservoir characterization model and the characteristics of the actual oil well production curve, a reasonable single-cave simplified model was initially selected(3)
In accordance with the initial static reservoir characterization model and the known dynamic and static description data, preliminary dimensionless parameters, such as large-scale dissolved caves and large fractures were assigned, and the dimensionless curves were initially matched with the field normalized curves(4)
We adjust the dimensionless parameters of large-scale dissolved caves, large fractures, or multiple porosity according to the initial matching and the analysis of the influence of various factors on the curve. We also recalculate the dimensionless curves. In addition, the normalized production rate curve is moved to fit the well actual curve and to record the given dimensionless cave parameters, fracture parameters, or multiple-porosity model parameters. We obtain the actual reservoir parameters through the following specific equations(14)kF=μB86.4rwqdqDM,rwa=3.6kFμϕFCtF+ϕVCtVtdtDM,LF=LDrwa,re=reDrwa,R=RDrwaπ/181/3.For the series model:(15)N=43πR3ϕV1−Swc+LFLF2WFϕF1−Swc.For the composite model:(16)N=1CttdtDdqdqDd+πR2hϕV1−Swc.Since we did not normalize the boundary control flow for the oil well rate transient analysis chart of the simplified discrete multiscale model (where seepage-free flow coupling is considered), we therefore have to comprehensively evaluate the initial static characterization model of the fractured-caved reservoir (i.e., basic oil well-setting data, seismic data, well logging data, and other dynamic and static data) to obtain more accurate reservoir parameters and oil well parameters in the final fitting.


### 3.2. Programming of Dynamic Evaluation Method

We use C# to compile a program to evaluate the corresponding reservoir parameters and oil well parameters in the discrete, simplified, and large-scale dissolved caved reservoir. The program includes four modules: crude oil and formation water physical parameter definition module, oil well production data input and analysis module, fitting model selection module, and specific matching module. The fitting model selection module can determine the number of large-scale fractures and caves and the contact relationship according to the known parameters. Then, the appropriate fitting model can be selected (Figure14).Figure 14
Flow chart of large-scale simplified model diagram matching.


## 3.1. Dynamic Evaluation Method and Processes

For the above discretized simplified models of different fractured-caved reservoirs developed in a large-scale single cave, the specific evaluation method and dynamic matching process of oil well and reservoir parameters are as follows:(1)
Obtain the oil well actual production data and calculate the normalized production rate, rate integral, and rate integral derivative curves [22, 23]. The specific calculation method adopts Equations (10)–(13)The equation for calculating the material balance pseudo-time of the oil well:(10)td=Npq.The equation for calculating the normalized production rate of the oil well:(11)qd=qΔp=qpi−pwf.The equation for calculating the integral of normalized production rate for an oil well:(12)qdi=1td∫0tdqddτ.The equation for calculating the integral derivative of normalized oil well production rate:(13)qdid=tddqdidtd.(2)
In accordance with the initial static reservoir characterization model and the characteristics of the actual oil well production curve, a reasonable single-cave simplified model was initially selected(3)
In accordance with the initial static reservoir characterization model and the known dynamic and static description data, preliminary dimensionless parameters, such as large-scale dissolved caves and large fractures were assigned, and the dimensionless curves were initially matched with the field normalized curves(4)
We adjust the dimensionless parameters of large-scale dissolved caves, large fractures, or multiple porosity according to the initial matching and the analysis of the influence of various factors on the curve. We also recalculate the dimensionless curves. In addition, the normalized production rate curve is moved to fit the well actual curve and to record the given dimensionless cave parameters, fracture parameters, or multiple-porosity model parameters. We obtain the actual reservoir parameters through the following specific equations(14)kF=μB86.4rwqdqDM,rwa=3.6kFμϕFCtF+ϕVCtVtdtDM,LF=LDrwa,re=reDrwa,R=RDrwaπ/181/3.For the series model:(15)N=43πR3ϕV1−Swc+LFLF2WFϕF1−Swc.For the composite model:(16)N=1CttdtDdqdqDd+πR2hϕV1−Swc.Since we did not normalize the boundary control flow for the oil well rate transient analysis chart of the simplified discrete multiscale model (where seepage-free flow coupling is considered), we therefore have to comprehensively evaluate the initial static characterization model of the fractured-caved reservoir (i.e., basic oil well-setting data, seismic data, well logging data, and other dynamic and static data) to obtain more accurate reservoir parameters and oil well parameters in the final fitting.


## 3.2. Programming of Dynamic Evaluation Method

We use C# to compile a program to evaluate the corresponding reservoir parameters and oil well parameters in the discrete, simplified, and large-scale dissolved caved reservoir. The program includes four modules: crude oil and formation water physical parameter definition module, oil well production data input and analysis module, fitting model selection module, and specific matching module. The fitting model selection module can determine the number of large-scale fractures and caves and the contact relationship according to the known parameters. Then, the appropriate fitting model can be selected (Figure14).Figure 14
Flow chart of large-scale simplified model diagram matching.


## 4. The Application of Conventional Type Curves in Simplified Large-Scale Discrete Fracture-Cave Model

We perform the dynamic evaluation of reservoir parameters using the W oil well of a fractured-caved reservoir in the Halahatang fractured-caved carbonate oilfield in Tarim Basin as a case study. Figure15 shows the production curve of well W with high oil pressure, stable pressure, and stable production.Figure 15
W well production dynamic curve: (a) daily oil production; (b) bottomhole pressure.
(a)(b)Based on the regular curve and static characterization model of well W, we can conclude that the well is indirectly connected to large caves through fractures. It is located on the west slope of a low uplift to the south of the uplift in the north of Tarim Basin. The productive zone is Ordovician and lies at a vertical depth of 6579 m. According to the current analysis, it can be simplified as a series model with a single cave and a single fracture, and this model can be used for fitting. Table1 shows the basic data of the W well. The final matching results are shown in Figure 16, and the parameters obtained by matching are shown in Table 2.Table 1
Basic data of W oil well.
Initial formation pressure (MPa)64Formation temperature (K)424.45Effective formation thickness (m)19Irreducible water saturation (%)20Porosity (%)40Well radius (m)0.086Rock compressibility (10−4 MPa−1)1.16Crude oil density (kg/m3)783.7Crude oil saturation pressure (MPa)14.6Crude oil viscosity (mPa·s)1.5345Crude oil volume factor (m3/m3)1.176Crude oil compression factor (10-4MPa-1)7.39Figure 16
Matching results.Table 2
Basic data of W oil well.
Fracture conductivity (mD·cm)1427.3Cave radius (m)90Oil well dynamic reserves (106 m3)1.12Fracture storativity ratio0.2Dissolved cave storativity ratio0.8Fracture length (m)100The dynamic controlled reserves of this well calculated by the material balance method are1.14×106m3, and the gap between this result and the result calculated by the method proposed in this article is relatively small. The type curve matching method can also accurately determine other caved parameters, large fracture parameters, and multiple surrounding porosity model parameters. If we consider the free flow of the cave, the result is more accurate. Through characterization of all wells in the Halahatang fractured-caved carbonate oilfield, we found that the single well dynamic reserves are mostly around 0.1×106m3~1×106m3, the permeability is mostly around 1 mD~100 mD (Figure 17), and the skin factor is mostly less than 0 (Figure 18). An accurate understanding of the scale of reserves and of reservoir parameters can further guide reservoir development.Figure 17
Histogram of single-well permeability evaluation results in Halahatang Oilfield.Figure 18
Distribution histogram of skin factor evaluation results of a single well in Halahatang Oilfield.


## 5. Conclusion

(1)
For fractured-caved reservoirs with large-scale single dissolved cave, we propose simplified models for two large-cave models—a beaded model composed of caves and large fractures and a composite model composed of caves and surrounding reservoirs(2)
The large-scale discrete fractured-caved reservoir is simplified. According to the static fractured-caved characterization, the scale of fractures and caves and the contact connection relationship can be obtained. Then, the reservoir can be simplified into different single-cave combinations. Free flow is considered in large-scale dissolved caves, and seepage is considered in large fractures and surrounding reservoirs. The corresponding simplified model is established, and the mathematical model considering the coupling of seepage and free flow is solved analytically. The rate transient analysis typical curves of the oil well are obtained. A systematic theory and method are formed for the dynamic evaluation of reservoir parameters and oil well parameters(3)
By matching the oil well actual dynamic production data in a fractured-caved reservoir with the typical rate transient analysis type curves, parameters such as equivalent radius, fracture permeability, and the multiple-porosity model storage capacity ratio of each cave can be calculated. Finally, the dynamic controlled reserves of the different single-cave reservoirs can be determined. These results provide parameters for rational development of the fractured-caved reservoir



---

*Source: 6611957-2021-01-25.xml*