Abstract

A number of different factors can affect flow performance in perforated completions, such as perforation density, perforation damage, and tunnel geometry. In perforation damage, any compaction at the perforation tunnels will lead to reduced permeability, more significant pressure drop, and lower productivity of the reservoir. The reduced permeability of the crushed zone around the perforation can be formulated as a crushed-zone skin factor. For reservoir flow, earlier research studies show how crushed (compacted) zones cause heightened resistance in radially converging vertical and horizontal flow entering perforations. However, the effects related to crushed zones on the total skin factor are still a moot point, especially for horizontal flows in perforations. Therefore, the present study will look into the varied effects occurring in the crushed zone in relation to the vertical and horizontal flows. The experimental test was carried out using a geotechnical radial flow set-up to measure the differential pressure in the perforation tunnel with a crushed zone. Computational fluid dynamics (CFD) software was used for simulating pressure gradient in a cylindrical perforation tunnel. The single-phase water was radially injected into the core sample with the same flow boundary conditions in the experimental and numerical procedures. In this work, two crushed zone configuration scenarios were applied in conjunction with different perforation parameters, perforation length, crushed zone radius, and crushed zone permeability. In the initial scenario, the crushed zone is assumed to be located at the perforation tunnel’s side only, while in the second scenario, the crushed zone is assumed to be located at a side and the tip of perforation (a tip-crushed zone). The simulated results indicate a good comparison with regard to the two scenarios’ pressure gradients. Furthermore, the simulations’ comparison reveals another pressure drop caused by the tip crushed zone related to the horizontal or plane flow in the perforations. The differences between the two simulations’ results show that currently available models for estimating the skin factor for vertical perforated completions need to be improved based on which of the two cases is closer to reality. This study has presented a better understanding of crushed zone characteristics by employing a different approach to the composition and shape of the crushed zone and permeability reduction levels for the crushed zone in the axial direction of the perforation.

1. Introduction

For cased hole completions, the well needs to be perforated in order for communication to be enabled between the formation and wellbore. This is because perforations enhance flow convergence at the near-wellbore region. In comparison with ideal open-hole wells, perforated wells can undergo added pressure gains and losses. Especially when a well has a short length of perforations with lower density, it may lead to a greater drop in pressure at the near-wellbore region, which will lower the overall well productivity. Moreover, when wells are perforated by deep penetrating tunnels, there may be a more expansive communication area between the formation and the perforated wells. When this occurs, the perforations might enhance well productivity due to less pressure drop.

Any added pressure gain or drop that is due to ideal perforation may be formulated as a perforation pseudoskin factor. In this case, this skin factor is related to parameters like wellbore diameter, perforation length, shot density, and phasing angle. Note that during the perforation process, rocks at the perforation tunnels’ sites are crushed. However, the perforation pseudoskin factor measures flow convergence at perforations in ideal perforated wells, neglecting the effect of the crushed (compacted) zone encircling the perforation tunnels. Experimental data indicate that the greatest contribution to total skin is the crushed zone, whereas formation damage skin and perforation pseudoskin make less contribution overall. Furthermore, based on studies (i.e., field and experimental) conducted over the past half-century, a major factor adversely affecting productivity is deficiencies in perforating procedures, perforator design, or both. For example, some studies on perforation damage showed that productivity declined in a gun-perforated hole formation when shots were fired into solid-containing fluid, such as drilling muds. There was also decreased productivity in cases where wellbore pressure exceeded formation pressure. In a typical process of perforation, a low-conductivity hole is created in the damaged formation. Postperforation, the crushed damage expands around the perforation radially from the center of the perforation tunnel. In this evolution, compacted and pulverized rock and other debris from a barrier block the formation’s natural pore spaces. Figure 1 illustrates this process [13].

The earlier studies mentioned above showed the importance of clearing perforating debris away from perforations to optimize the flow capacity. Debris removal, whether of some or most of the damaged material, enables the perforation tunnel to better perform its role as the wellbore’s fluid conduit. One highly efficient way to remove the debris is through underbalanced perforating, although the precise underbalance degree required for generating effective perforations is still a topic of debate. The underbalance is usually calculated by considering the matrix permeability, along with tunnel parameters and fluid, as described by Tariq [4], Walton [5], and Grove et al. [6]. More recent research has looked into various aspects of cleanup and productivity mechanisms in relation to perforation tunnels. Several studies have used computational analysis, due to a lack of experimental data, to explore a number of different challenges, such as isolating the competing effects of debris cleanup and dealing with the heterogeneity that often characterizes analog cores. One of the research teams who used the computational approach is Satti et al. [79]; they proposed a 3-D flow model to determine flow impediments based on characteristics of perforation damage, debris and blockages, and tunnel geometry. The authors used the latest micro-CT methods as well as simulations to determine the most viable cleanup strategies for optimal productivity.

Interactions between formation damage, flow convergence near perforation tunnels, and rock compaction need to be carefully formulated; rock compaction and formation damage can create a higher degree of dynamic interaction between plane flow and flow convergence at the perforations. The effects of rock compaction and formation damage should thus be included in any perforation model that is developed. However, these interactions can be extremely complex in nature, particularly when including parameters such as mechanical skin factor (formation damage), ideal perforation pseudoskin factor, and crushed-zone skin factor. Overall, perforation total skin factor () can be used to represent the combined impacts of reductions in compaction zone permeability near perforation tunnels, ideal perforations, and formation damage caused by drilling and other production operations near wellbores. In the Darcy flow case, any other pressure changes caused by these impacts in combination may be written as

A wide range of models (empirical, numerical, experimental, semianalytical, etc.) have been developed over the past several decades for predicting perforated total skin and for determining the effects of the crushed zone. For example, Klotz et al. [10] employed a 2-D finite element model to examine how perforations and formation damage at the crushed zone affect well productivity, while Locke [11] looked at a novel way (nomograph) to estimate skin for perforated completions. Nomographs can take skin factors into account by considering rock compaction and formation damage at perforation tunnels. The nomograph results can be integrated into models that measure well performance. A decade later, in a study conducted by Thomas et al. [12], a table look-up procedure was used for combining the results of Locke’s nomograph and software developed for well performance. McLeod [13] suggested that the combined effects of rock compaction, formation damage, and perforation pseudoskin at perforation tunnels could be captured in Equation (2), as written below:

The crushed-zone skin factor is expressed in

Using a finite element simulator, Karakas and Tariq [14] formulated several empirical equations for predicting total skin factor at fully perforated vertical wells. One example of their formulations was an equation that predicted perforation pseudoskin in cases where there was no drilling damage.

In this equation, the initial term indicates flow convergence in the horizontal plane, the second term denotes flow convergence in the vertical plane, and the third term describes the effect of the wellbore geometry. The researchers also employed a formulation somewhat like Equation (3) above to express the crushed-zone skin factor. This formulation is written as Equation (5) below:

In other related research, Bell et al. [15] formulated a skin factor equation to account for perforations that ended within the confines of drill-damaged zones. This formula was based on the previous Karakas and Tariq model [14] mentioned earlier and considered formation damage and crushed zone effects in combination. In 2006, Yildiz [16] proposed a novel way to assess total skin factor that accounted for a large number of other factors, including compaction zone skin, perforation pseudoskin, drill-damaged skin, and the impact of partial formation perforation. However, despite their achievements, the effects related to crushed zones on the total skin factor are still a moot point, especially for horizontal flows in perforations. Also, computational fluid dynamics (CFD) investigations have been conducted by Sun. et al. [17] for a perforated vertical well with different crushed zone scenarios. Their results showed an apparent discrepancy between the supposed crushed zone scenarios and available models.

All the preceding factors have derived several models based on the assumption that a crushed zone is uniform and cylindrical and has a homogeneous permeability and porosity. At the same time, recent studies using modern techniques computerized tomography (CT) and a scanning electron microscope (SEM)) have shown a different perception of the composition and shape of the crushed zone, especially in the axial direction of the perforation tunnel. In recent studies, the researchers have analyzed the results obtained based on the perforation process’s effect. The process of perforation involves shooting a perforation agent into the near-wellbore formation to make holes. However, during this process, the affected rock could experience grain bond breakages in addition to microfractures caused by shock waves [18]. Four levels of damage (i.e., four permeability and porosity change zones) were found in each of the simulated cases when permeability changes, porosity, and bond breakage in the axial direction and vicinity of the perforation tunnel. In Zone 1, there was a significant and rapid decrease in porosity after the minimum value was achieved. The force of the shock wave in this zone was sufficient to sever every bond between particles. The rocks typically had numerous microcracks, which became angular fragments during backward grain displacement. The result of this process was a major reduction in permeability due to the closing of pore throats. Zone 2 featured less variation in porosity compared to Zone 1. Even though most of the particle bonds in Zone 2 broke, a small portion of the grains did stay intact. Furthermore, after achieving a minimum value, both permeability and porosity saw slower increases in Zone 2 in comparison with Zone 1. In Zone 3, there were overall fewer decreases in permeability and porosity compared to Zones 1 and 2. Also, in Zone 3, the permeability and porosity generally stayed constant postperforation. The reason for these differences between Zone 3 and Zones 1 and 2 is that the force of the shock wave was attenuated prior to arriving at Zone 3. Moreover, even though many bonds broke in Zone 3, the grains stayed intact. From this, we can see that the porosity decrease in Zone 3 was likely caused by particle adjustment, where the particles stay constant due to forces from particles situated around them. In Zone 4, no bond breakage or porosity damage was detectable, as this zone is too far away for the shock wave to reach, leaving the rock matrix intact. Figure 2 illustrates the four mentioned crushed change zones described by Sarmadivaleh et al. [19], Pucknell and Behrmann [20] Nabipour et al. [21], and Craddock et al. [22].

In addition, Xue et al. [23] developed a mechanical model to describe influences of perforating damage, especially with regard to perforations caused by explosions affecting the permeability and porosity of sandstone at the compaction zone. Their study results indicated that there was a notable decrease in damage, particularly to the thickness of the sandstone, when the load and pressure decreased in the axial direction of perforation tunnel, and these results confirm Sarmadivaleh et al.’s [19] and others’ conclusions.

In addressing this research gap, the present study was conducted experimental and numerical investigations by employing a different approach to the composition and shape of the crushed zone and permeability reduction levels for the crushed zone in the axial direction of the perforation. The two crushed zone configuration scenarios were applied in conjunction with different perforation parameters, perforation length, crushed zone radius, and crushed zone permeability. In the initial scenario, the perforation zone’s tip is assumed to be too far away for the shock wave to reach, leaving the rock matrix intact, so the crushed zone is located at the perforation tunnel’s side. In the second scenario, the force of the shock wave in the perforation zone’s tip is considered to be sufficient to sever every bond between particles, so the crushed area is located at a side and tip of perforation. The two scenarios are applied in conjunction with different perforation parameters, perforation length, crushed zone radius, and a permeability ratio. In additional investigations, the effect of permeability anisotropy in crushed zones on the crushed skin factor has been studied considering these scenarios.

2. Methodology

2.1. Experimental Procedure

The experimental study was used to validate the numerical model for single-phase flow through the perforation tunnel that includes two different permeability zones. The first region surrounds the perforations and represents the crushed zone with low permeability; the second region exemplifies the formation region with high permeability. Statistical analysis was coupled with numerical simulation to expand the investigation of fluid flow in the near-wellbore region due to the limitations of the experimental setup, especially the small sample size. In the study, two crushed zone configuration scenarios are conducted in conjunction with different perforation parameters, perforation length, crushed zone, and permeability ratio, as shown in Figure 3.

In the present study, the experimental set-up initially designed and built by Ahammad et al. [24, 25] as a radial flow cell (RFC). The RFC, which was constructed at Memorial University of Newfoundland’s Drilling Technology Lab, was created for the purpose of carrying out experiments under wellbore conditions. Figure 4 shows the set-up, which features the three following main sections: flow lines extending from inlet to outlet, an inner chamber for holding samples with axial loads, and a data acquisition (DAQ) system.

In the laboratory experiments, a highly permeable synthetic sandstone sample was prepared. The cylindrical sample was made of sand particles measuring 0.18 to 1.18 mm. The synthetic sample has been created from two different sandstone grain sizes. The fine grain size used to create a crushed zone around the perforation and the coarse grain size used to create a formation reservoir zone are demonstrated in Figure 5. The dimensions and the index properties of the sample included permeability and porosity, as shown in Table 1.

Experiments carried out on perforation methods have primarily relied on rather simplistic assumptions, such as those presented by Rahman et al. [2630]. Moreover, as a result of laboratory constraints, most experimental investigations have neglected key reservoir characteristics, such as thermal effects, drawdown pressure, and actual reservoir pressure. In the experimental portion of our work, we injected a measured volume of water into our core sample. Also, we used a geotechnical radial flow test set-up to measure the differential pressure and single-phase flow rate of our synthetic sandstone sample, with water being radially injected into our core sample within Darcy flow and flow boundary conditions (L/min,  mPa·s). The outer side of the sample is considered an inlet while the perforation surface is an outlet. Furthermore, both inlet and outlet pressures were measured for our cylindrical samples using a specified water flow rate.

2.2. Numerical Procedure

In the present work, we used Ansys Fluent 18.1 for our computational fluid dynamics (CFD) model. Our aim was to present a single-phase fluid flow simulation for a reservoir described as three-dimensional, vertical, and cylindrically layered. We created a sample that is vertical with a single layer of uniform thickness (). Next, we assumed the well was centrally located and had a drilled and casted radius throughout the formation. Hence, we could apply the conditions and assumptions enumerated below in developing our model: (1)The formation is isotropic and porous, of uniform thickness, and is constantly permeable (i.e., features constant vertical permeability that is nonzero)(2)The flow through the reservoir can be described as single-phase water and either radial-vertical laminar or Darcy’s flow(3)The liquid is incompressible with a constant viscosity(4)Any flux proceeding into the well features uniform distribution across perforated intervals(5)Thermal effects were ignored(6)The crushed skin has been considered in the present research. Other skin factors and effects of perforations angle, formation permeability anisotropy, and wellbore radius were neglected

The crushed skin factor is affected by crushed zone parameters and the permeability anisotropy of the crushed zone. In contrast, the perforation skin factor is more affected by perforation angle, formation permeability anisotropy, and the wellbore radius. Therefore, additional CFD investigations analyzed the effect of permeability anisotropy in crushed zones on the crushed skin factor, considering the study’s two mentioned scenarios.

In the numerical work, we injected a measured volume of water into the cylindrical sample. The conservation equations for mass and momentum describing single-phase flow in a porous region could be expressed, respectively, as

The last term in Equation (7) represents the viscous and inertial loss imposed by the porous media on the fluid. The laminar flows in porous media generally feature a pressure drop proportional to permeability () and velocity (). By ignoring out an inertial loss term, we can reduce the porous media model to Darcy’s law, as expressed in

For the three coordinate directions (, , and ) of the porous region, the pressure drop calculated by Ansys Fluent can be expressed as in where indicates velocity components for , and directions and , , and denote the medium thickness for the, , and directions.

Uniform mesh and cut mesh methods (Figure 6) were used to generate high-quality mesh. This configuration helped to predict a good-quality and high-density mesh close to perforation borders.

2.3. Statistical Procedure

Various methods for examining how different parameters may affect experimental results are applied by using Design of Experiments (DoE) software. The initial step in DoE is identifying independent variables and/or factors that may affect the experimental outcomes. The next step involves identifying the dependent variables and/or factors [31]. The experiments usually run based on different factor values or levels, such that each run features a combination pertaining to the specific factor value(s) or level(s) being investigated. In the present work, we used Design Expert software with BBD in designing the runs needed for statistical analysis. The Box-Behnken design (BBD) is described as a response surface methodology (RSM) design and requires only three levels in its experimental or numerical runs [32, 33].

The model is then statistically validated through analysis of variance (ANOVA) [34]. Based on these initial processes, three dimensionless parameters were investigated, including the ratio of penetration space to perforation length , the ratio of crushed radius to perforation radius , and the crushed-zone damage permeability ratio:

Two boundary points were then selected, and one midpoint was determined by BBD for the intervals of the parameters, as presented in Table 2.

3. Results and Discussion

The comparison between the experimental and numerical results of the pressure buildup with the same flow boundary conditions is shown in Figure 7. The dimensions of the perforation geometry and the index properties are the same as those used in the experimental procedure (see Table 1). This research investigates a specific water flow rate (1 L/min) that was injected into the sample to determine the differential pressure. The experimental data and numerical results are in good agreement.

The validation of numerical results with experimental ones has given full confidence in using the numerical model to conduct huge investigations by creating a crushed zone with different crushed perforation parameters, perforation length, and crushed zone radius. The relative effect of three dimensionless parameters on the crushed skin factor was investigated before conducting statistical analysis. The numerical results showed that the ratio of penetration space to perforation length and the crushed-zone damage permeability ratio significantly affect the crushed skin factor. In contrast, the ratio of crushed radius to perforation radius has a moderate effect on the crushed skin factor, as shown in Table 3.

Twelve numerical runs were performed and analyzed to obtain a suitable statistical analysis using the ANOVA analysis with the BBD model (Table 4).

Therefore, crushed perforation parameters are analyzed by using statistical analysis coupled with the numerical simulation model. This study provided two correlations from the statistical analysis, based on the numerical results. These correlations were used to determine the relative impact of each factor for the two scenarios on the crushed skin factor.

The results clearly indicate the crushed perforation parameters effects on both the crushed skin factor and pressure gradient values. The results show that crushed skin increases the pressure drop and thus contributes to a reduction in the productivity index. As illustrated in Figure 8, the value of the crushed skin factor increases with the increase of the penetration ratio (), the ratio of the crushed zone to perforation radius (), and crushed-zone damage permeability ratio (). From this, we can surmise that the interaction among these parameters has a marked effect on the crushed skin factor value, whereas the perforations without a crushed zone represent an ideal case that may or may not be reproducible in practice. Furthermore, based on the results, we can see that the pressure gradient is more affected by a high ratio for the three parameters. The perforation length makes a more significant contribution to pressure drop decreasing, and the higher crushed-zone damage permeability ratio with short perforations leads to pressure drop increases around the perforations, as shown in Figure 9. Also, the results show the significant effect of a large thickness of crushed zone with a high crushed-zone damage permeability ratio on the pressure gradient due to the frame’s resistance around the perforations with a significantly high reduction in the permeability.

Moreover, the numerical results show a clear view of pressure distribution for the perforation with a crushed tip, without a crushed tip, and ideal perforations cases. For example, the pressure gradient for the three cases at dimensions’ parameters values (, , and =100) is shown in Figure 10.

In order to compare and discuss the accuracy of the common models, one model by Karakas and Tariq [14] was selected and used to calculate the crushed skin factor without formation damage. The model was compared with two novel correlations for two scenarios. The comparison results show a comprehensive realization of the effect of three dimensionless parameters (, , and ) on the crushed skin value and the deviation between the Karakas and Tariq [14] model and novel correlation for two scenarios.

The crushed skin factor results of the first correlation with a tip-crushed zone showed a slight deviation with the previous model of Karakas and Tariq [14] while the crushed skin results of the second correlation without a tip-crushed zone () demonstrated a large deviation for three dimensionless parameters, as shown in Figure 11. In general, the Karakas and Tariq [14] model applied to calculate the crushed skin factor shows a good convergence with the second scenario, so the model takes into consideration the effect of all frame resistance around the perforations. However, our CFD model gave almost identical results to the model if we assumed the length of the perforations was shorter than their length by a thickness of the crushed zone.

Also, the present study looks at the effect of permeability anisotropy in crushed zones on the skin factor with regard to the study’s two scenarios. In the CFD simulations, the interaction effect between permeability anisotropy in the crushed zone and the crushed-zone damage permeability ratio has been investigated. As illustrated in Figures 12 and 13, there is a seeming increase in crushed skin factor when permeability anisotropy in the crushed zone is assumed. It is plausible that this rise’s mechanism is related to nonradial flow near the perforation tunnel caused by permeability anisotropy at the crushed zone, particularly in cases of increased crushed zone thickness and perforation length. The results showed that the anisotropy in the crushed zone between has a significant effect on the crushed skin factor, and then, its effect decreases, due to the domination of the horizontal flow. These outcomes showed that a better understanding of crushed zone anisotropy is needed, using the improvement of techniques for determining which factors affect perforated completions and to what extent.

4. Conclusions

The study was conducted in order to further investigate the accuracy of Karakas and Tariq’s [14] model for crushed skin factor as well as to obtain new correlations which may be closer to reality, established on some assumptions. Based on the results of this investigations’ analysis, the following conclusions can be summarized: (1)The experimental data showed good agreement with the numerical model results used in this work to conduct more investigations(2)The study showed a clear view of the effect of the three dimensionless parameters (, , and ) on the crushed skin factor and pressure gradient(3)The comparison of the simulations reveals that there is a significant difference between each of the two tip-crushed zone scenarios(4)The numerical model gave almost identical results as for Karakas and Tariq’s model, if the length of the perforations was assumed to be shorter than their real length by a thickness of crushed zone(5)The differences between the two simulations’ results show that the currently available model [14] for estimating the skin factor for vertical perforated completions needs to be improved, based on which of the two cases is closer to reality(6)The study presented two novel correlations that give more than one option to calculate the crushed skin factor(7)The outcomes of this study underscore the need to include the crushed zone anisotropy effect through the improvement of available models for determining the crushed skin factor

Nomenclature

:Body force
:Inertial resistance factor
:Formation thickness
:Permeability of crushed zone
:Horizontal permeability of crushed zone
:Vertical permeability of crushed zone
:Permeability of damaged zone
:Formation permeability
:Permeability ratio
:Perforation length
:Number of perforations (number of shots per foot)
:Pressure
:Penetration ratio
:Flow rate
:Radius of crushed zone around perforation
:Reservoir radius
:Perforation radius
:Ratio of crushed zone to perforation radius
:Wellbore radius
:Skin due to rock crushed around perforations
:Skin due to formation damage
:Skin due to horizontal flow effect
:Skin due to ideal perforations
:Total perforation skin factor
:Skin due to vertical converging effect
:Skin due to wellbore effect
:Skin due to crushed zone without tip
:Skin due to crushed zone with tip
:Time
:Velocity
:Velocity components for , , and directions
:Porosity
:Porosity of formation zone
:Porosity of crushed zone
ϴ:Perforation angle
:Fluid density
:Stress tensor related to viscous flow
μ:Fluid viscosity
:Medium thickness
:Pressure drop due to total skin factor
:Pressure drop (ideal perforations without crushed zone)
:Pressure drop due to crushed zone without tip
:Pressure drop due to crushed zone with tip
Abbreviations
ANOV:Analysis of variance
BBD:Box-Behnken design
CFD:Computational fluid dynamics
DAQ:Data acquisition
DoE:Design of experiments
RFC:Radial flow cell
RSM:Response surface methodology.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Disclosure

The authors certify that this work is original, has not been published, and will not be submitted elsewhere for publication. Statements made herein are solely the responsibility of the authors.

Conflicts of Interest

The authors certify that they have no known competing financial interests that could have influenced the work reported in this manuscript.

Acknowledgments

This publication was also made possible by the grant NPRP10-0101-170091 from Qatar National Research Fund (a member of the Qatar Foundation).