Wellbore and site safety must be ensured during CO2 injection into multiple reservoirs during carbon capture and storage projects. This study focuses on multireservoir injection and investigates the characteristics of the flow-rate distribution and reservoir-risk evaluation as well as their unique influences on multireservoir injection. The results show that more CO2 enters the upper layers than the lower layers. With the increase in injection pressure, the risks of the upper reservoirs increase more dramatically than those of the low reservoirs, which can cause the critical reservoir (CR) to shift. The CO2 injection temperature has a similar effect on the injection flow rate but no effect on the CR’s location. Despite having no effect on the flow-rate distribution, the formation-fracturing pressures in the reservoirs determine which layer becomes the CR. As the thickness or permeability of a layer increases, the inflows exhibit upward and downward trends in this layer and the lower layers, respectively, whereas the inflows of the upper layers remain unchanged; meanwhile, the risks of the lower layer and those of the others decrease and remain constant, respectively. Compared to other parameters, the reservoir porosities have a negligible effect on the reservoir risks and flow-rate distributions.

1. Introduction

Carbon capture and storage (CCS) is widely recognized as an effective approach for greatly reducing CO2 levels in the atmosphere [15]. Many CCS projects have been conducted worldwide, including the Sleipner project [6, 7] in Norway, the Weyburn project in Canada [8], the Otway Pilot project in Australia [9], the In Salah project in Algeria [6], and the Shenhua CCS demonstration project in China [10]. Wellbore and site safety must be ensured in all fluid injection projects [11, 12]. In many CCS projects, particularly those with large-scale CO2 injection, multiple reservoirs are employed for simultaneous injection to achieve a preset injection target amount of CO2 (see [1316]). Layers of caprocks and reservoirs are sequentially spaced to form multiple suits of cap rock-reservoir combinations, which significantly increase the complexity of fluid migrations [14]. A mature design methodology for safe and effective CO2 injection through deep wellbores requires an in-depth understanding of the reservoir performance and safety and their influencing factors in terms of CO2 injection.

The injection of CO2 with multiple layers differs considerably from single-layer injection in several respects. First, safe injection requires that the pressure in each reservoir does not exceed the maximum allowable value. Second, the flow rate of each reservoir is unknown in advance of the initial injection. Third, the most dangerous reservoir, based on which wellbore working parameters should be used, is also not known in advance. Therefore, the analysis and evaluation of the risks involved with multiple reservoir injection are difficult.

Bai et al. [14] derived an analytical solution for two-phase flows based on the work of Nordbotten et al. [17], and a solution was derived to designate the wellhead injection pressure. These authors characterized the risks of a reservoir using the ratio of the actual pressure of the fluid to the maximum allowable pressure. A higher ratio indicates that the corresponding reservoir is more dangerous. The reservoir that has the largest ratio is the weakest reservoir, which is also defined as the critical reservoir (CR). As the shortest slab in terms of the Cask principle, the CR actually constrains the maximum allowable wellhead injection pressure. Although the actual injection pressures are typically less than the maximum allowable value, the injection flow rates can vary, which may lead to variation of CR’s location. On such occasions, one cannot determine which layer is the CR because the downhole pressures and flow rates of all the reservoirs are initially unknown. Therefore, investigations of the effects of the injection parameters on the allocations of the total flow rates among multiple layers are extremely valuable.

This paper investigates the determination of CRs and their influencing factors during multireservoir CO2 injection. For this purpose, the allocation percentage of the total flow rates among the multiple reservoir layers and the risk evaluation of the reservoirs must be simultaneously investigated. Such an investigation is expected to enrich the analysis and design methodologies for CO2 injection operations. First, we outline the basic theory of the wellbore pressure and temperature. Then, a base case for determining the CR is presented. Next, we analyze the influences of the injection parameters and reservoir properties on the allocation percentage of the total flow rate, the reservoir risk, and the position of the CR. Finally, we summarize our main findings and conclude the study.

2. Fundamental Theory

To achieve the goal of this study, the coupling calculations for the flow and heat transfer of the wellbore and formation must be performed. Although various methods, including three-dimensional numerical simulations, could be employed, a semianalytical method developed by Bai et al. [14] and based on a fast, explicit numerical method and an analytical solution is used to calculate the wellbore flow and heat exchange between the wellbore and formations in this paper. The method is fast and can be applied successfully. In this method, the wellbore is discretized into a series of one-dimensional elements, and the model should include the following assumptions:(1)One-dimensional flow with homogenous fluid is assumed in the vertical wellbore, and all the state variables and properties are assumed to be uniform at the same section.(2)Only radial heat transfer is considered, and the temperature at each point is updated using Ramy’s solution [18].(3)The influence of the phase change on the fluid properties is not considered.(4)All of the fluid properties should be constant within an element.(5)The portion of the wellbore in the reservoir is simplified as an element node; that is, the variation in the state variables of the wellbore fluid is held constant in a reservoir.

The CO2 pressure can be determined when the injection pressure at the wellhead is known; then, the fluid pressure at the end point of the -th well segment can be acquired from the following equation when is already known:where is the mass flow velocity of the wellbore cross section (kg·m−2·s−1), is the discretized segment (m), is the acceleration of gravity (m/s2), is the inner radius of the tubing (m), is the gas molar mass (kg/mol), is the universal gas constant, is the coefficient of friction, is the compression factor obtained by solving the Peng-Robinson equation [19], and is the thermodynamic temperature (K). The subscript is used to number the discretized segments of the wellbore in the finite difference method.

As noted above, Ramy’s analytical solution is used to obtain the CO2 temperature of the wellbore. The details of the derivations are provided in the studies of Liu et al. [15] and Wu et al. [16], and the heat transfer between the wellbore and the surrounding earth is detailed in the study by Streit and Hillis [13].

The steady two-phase flow for the CO2 flooding in the reservoir can be characterized by an equation proposed by Nordbotten et al. [17], with the assumption that a reservoir can be divided into CO2-saturated and brine-saturated zones with a sharp interface. Obtaining an analytic solution from Nordbotten’s equation is almost impossible. According to the derivation, by introducing two-phase mobility into the Darcy formula of a single-phase flow and assuming that the CO2 plume is radially symmetric, as suggested by Wu et al. [20, 21] and Bai et al. [12], the mass flow rate , which enters the -th reservoir from the wellbore, can be expressed as follows:where is the initial formation pressure of the reservoir (Pa); is the injection pressure on the wellface of the reservoir (Pa); is the maximum influence radius of the flow in the reservoir (m); and are the radii of the CO2 plume at the bottom and top of the reservoir, respectively (m); is the absolute permeability of the reservoir (m2); is the thickness of the reservoir (m); is the density of the CO2 (kg/m3); and and are the mobility of CO2 and the mobility of brine (m·s/kg), which are defined as the ratios of their relative permeability to their fluid viscosity, that is, and . Since the saturations of brine in the CO2 and brine domains of a reservoir are zero and one, respectively, the relative permeabilities of and are both equal to one.

Moreover, three of the parameters in (2), , and , and , are used to describe the distributions of CO2 and brine in a reservoir. These parameters are time-dependent and can be calculated as follows:where is the injection time (s); is the total flow into a reservoir during 0– (m3); is the reservoir porosity; is the pore compressibility (m2/N); is the compressibility of brine (m2/N); and is the viscosity of brine (kg/m/s).

According to the flow equilibrium conditions, the inflow of the section of the wellbore in the -th reservoir is equal to the summation of the outflow into the next well segment and the corresponding reservoir [12]. Therefore, the following equation is obtained:

When (4) is used for the bottom reservoir, the second term on the right-hand side disappears.

Typically, the site stability and the flow-rate target must also be satisfied and form the constraints during CO2 injections. The wellbore constraints include the pressure and flow-rate constraints. The lower limit of the latter can be determined during the project-feasibility stage, whereas the upper limit requires each branch flow to meet the corresponding constraint conditions. Therefore, the complete wellbore constraints can be expressed aswhere represents the capillary pressure (Pa); is the maximum allowable pressure of the i-th reservoir (Pa), which is the formation-fracturing pressure multiplied by a synthetic design coefficient ; and is the available reservoir capacity (kg/s).

3. Example of the Shenhua CCS Demonstration Project

The Shenhua CCS demonstration project is the first fully implemented CCS project in China; the storage target of this project was set to 100,000 tons/year. A single injection well, which is referred to as ZSZ1, was drilled to a depth of 2,450 m, penetrating 21 reservoir-cap rock pairs, which were then combined and reduced to eight reservoir-cap rock units for analysis by Bai et al. [14]. From top to bottom, the geological formations include the Zhifang group, Heshanggou group, Liujiagou group, Shiqianfeng group, Shihezi group, Shanxi group, Taiyuan group, Benxi group, and Majiagou group. Table 1 lists the computational parameters of the reservoir-cap rock units from top to bottom. The other parameters, that is, the inner radius of the tubing , pore compressibility  , brine compressibility , CO2 viscosity , and brine viscosity , were set to 31 mm,  m2/N,  m2/N, 88 kg/m/s, and 552 kg/m/s, respectively.

As noted above, a forward analysis is employed in this study. Therefore, the parameters of the wellhead, such as the wellhead injection pressure and injection flow rate, must be prepared before the aforementioned method can be used to obtain the pressure and CO2 injection flow rates of the reservoirs. Based on these values, the ratio can be obtained for the -th reservoir. The synthetic design coefficient for obtained is 0.8. The boundary conditions, that is, the injection pressure and injection rate at the wellhead, are 5 MPa and 4.63 kg/s, respectively, in the base case, which is consistent with the actual operational values at a certain injection stage. In addition, the calculated injection time is 3 years.

According to the distribution of the flow rates in the eight reservoirs plotted in Figure 1, the CO2 flows are greater in the upper reservoirs under the given injection parameters. The flow rates of the upper four reservoirs account for 59.2% of the CO2 flow, although the sixth reservoir has the largest flow rate, accounting for 23.2% of the total. Reservoirs 5, 7, and 8 comprise considerably lower percentages of the total flow rate.

The risk factors of all eight reservoirs, which are defined as the ratios of the actual pressures of the fluid to the maximum allowable pressures at the entrances of the reservoirs, are shown in Figure 2. Reservoir 8 has the highest risk factor and is thus the CR. Bai et al. [14] found that a different reservoir was the CR when using different injection parameters. Hence, the CR is closely related to the injection parameters. Because all the risk factors are less than 1.0, the given combination of injection parameters will not induce reservoir failure; the CR has a safety reserve of nearly 18%. In addition, Reservoir 8, which exhibits the minimum proportions of CO2 entry, is the most dangerous reservoir, which indicates that the pressure might build up sharply at the wellface because of the poor injectivity under the given injection conditions.

4. Analysis of the Influencing Factors

Because the injection parameters are typically not constant during the actual injection operations, it is valuable to know how the risks and proportions of the flow rates in each reservoir depend on the injection conditions. Knowing whether the CR will shift because of the variations in the injection parameters is more important. In addition to these injection parameters, the characteristic parameters of the reservoirs are believed to be important influencing factors. The typical characteristic parameters of reservoirs include the reservoir thickness, porosity, permeability, initial formation pressure, and formation-fracturing pressure [2]. The two classes of factors, that is, the injection and characteristic factors, will be investigated in a subsequent study. The injection fluid is assumed to be pure CO2, so the effects of impurities are not considered. Although the saturation of the fluid can significantly affect the distribution of the wellbore pressure and the wellhead pressure in particular, this factor will not change the CR [12]. Therefore, the effects of the characteristic physical parameters of the fluids will not be further investigated in this study.

4.1. Injection Parameters

In this subsection, we utilize the wellbore model and analysis method presented in the previous section. The main injection parameters of interest in this section are the injection pressure, injection flow rate and injection temperature. Ten injection cases are described in Table 2 based on combinations of these three parameters. The settings of the parameter values, which are chosen by considering the base case in the previous section and the actual injection history, are likely to be used in actual injection practices. Therefore, the results should provide valuable guidance for any subsequent injection operations during the Shenhua CCS project. The cases shown in Table 2 were designed with four targets, that is, determining the effect of variations in the injection pressures (Cases A to D), injection rates (Cases D to F), injection temperatures (Cases G to J), and combinations of the first two factors (Cases C, E, and G) on the percentages of the CO2 inflows and the risks of each reservoir. All the designed parameter values ensure the safety of each reservoir because the target of this study is not the failure properties of the formations but the sensitivities of the CR to the parameters. The responses of the percentages of the CO2 flow rates and the risk factors for each reservoir are illustrated in Figures 3 and 4, respectively.

Figure 3(a) shows that the proportion of CO2 that flows into each reservoir varies with changes in the injection pressure under the given flow rates and temperatures at the wellhead. With greater injection pressures, the inflows into Reservoirs 1–5 increase monotonically, the inflows into Reservoir 6 initially increase and then decrease, and the inflows into Reservoirs 7-8 gradually decrease. These variation trends indicate that a greater injection pressure causes more CO2 to enter the upper reservoirs because of the increase in the fluid pressure at the wellbores in the reservoirs and the constant total flux at the wellheads. Figure 4(a) indicates that the reservoir risk also increases due to the increased injection pressure. More importantly, the risk-increasing trends of the upper reservoirs are more dramatic than those of the lower reservoirs. When the injection pressure reaches 7.0 MPa, the CR transfers from Reservoir 8 to Reservoir 1, indicating that the injection pressure can control the overall risks of the reservoirs and that concentrating CO2 in the upper reservoirs causes the CR to shift to the top reservoir because of the increased injection pressure.

Figure 3(b) shows the variations of the distributions of the CO2 flow rates among the reservoirs with increases in the injection flow rates when the injection pressure and injection temperature are known. In contrast to Figure 3(a), a higher total injection flow rate increases the proportion of CO2 that flows into the lower reservoirs. This is mainly because the pressure attenuation increases with the injection flow rate such that the flows into the upper reservoirs are reduced. Moreover, as shown in Figure 4(b), the risks of all the reservoirs gradually decrease from Case D to Case F, and the CR shifts back to Reservoir 8, which differs from the results presented in Figure 4(a). However, the risk reductions of these reservoirs are relatively small, particularly for the lower reservoirs, even when their flow rates are increased. Thus, the injection flow rate has a smaller impact on the risk of a reservoir than the injection pressure does, and both parameters have larger impacts on the risks of the upper reservoirs than on the risks of the lower reservoirs.

As the injection temperature of CO2 increases, the variation trends of the flow rates in each layer become opposite to those induced by the increases in the injection pressure, despite agreeing well with those shown in Figure 3(b) for the gradually increasing concentrations of CO2 in the upper reservoirs, as shown in Figure 3(c). Similarly, Figure 4(c) indicates that the changes in the reservoir risks are approximately opposite to those in Figure 4(a). This decrease in the reservoir risk is likely induced by a drop in the viscosity, which can accelerate the transport of fluids and lead to a partial dissipation of pressure because of the increasing injection temperature. Furthermore, Reservoir 1 remains the CR, even though the reservoir risk decreases with the increasing injection temperature.

Given the synchronous variations in injection pressures and flow rates, the injection proportion of each reservoir remains almost invariant, except for that of Reservoir 6, which exhibits a 7.5% decrease in Case G, as shown in Figure 3(d). However, the risks of all the reservoirs in Figure 4(d) increase significantly, but Reservoir 1 remains the CR. In combination with the results shown in Figures 4(a) and 4(b), the injection pressure is shown to predominantly control the overall risks of the reservoirs, and the injection flow rate is shown to be able to affect the position of the CR to a certain extent.

4.2. Characteristic Parameters of the Reservoirs
4.2.1. Parameter Settings

The injectivity of the reservoir and the downhole limitations of the injection pressure depend on the reservoir’s characteristic parameters; variations in these parameters can affect the intake percentages and risks of each reservoir. For a CCS project, these parameters are typically considered to be constant. The objective of this section is to determine the effects of these parameters; therefore, the assumed parameters rather than the experimental parameters of the Shenhua CCS project are employed to conduct a parameter sensitivity analysis. A conceptual engineering model that contains four reservoir-cap rock units with a single injection well that penetrates 2,000 m underground is established. The formations in the conceptual model, from top to bottom, are as follows: Cap rock 1, Reservoir 1, Cap rock 2, Reservoir 2, Cap rock 3, Reservoir 3, Cap rock 4, and Reservoir 4. The depths of the bottoms of these formations are 1390 m, 1400 m, 1590 m, 1600 m, 1790 m, 1800 m, 1990 m, and 2000 m, respectively. The basic characteristic parameters of the reservoirs are listed in Table 3, and the parameters of the wellbores are the same as those for the Shenhua CCS project. For the base case, the wellhead injection pressure, injection flow rate, and injection temperature of CO2 are set to 10 MPa, 4 kg/s, and −5°C, respectively. Five cases with the same initial states as the base case (Table 3) are designed to study the influence of each reservoir’s parameters. The injection time is set to ten years.

4.2.2. Influence of the Formation-Fracturing Pressure

As presented above, the concept of the CR depends on the ratio. is typically obtained from the formation-fracturing pressure multiplied by the synthetic design factor, which is smaller than 1.0. In this conceptual model, the synthetic design factor is set to 0.6. Therefore, represents the intrinsic ability of a formation to resist damage and is a decisive factor in determining the risk of a reservoir and the position of the CR, which will be investigated first. Five groups of fracturing pressures are designed and listed in Table 4. In Cases B and C, the formation-fracturing pressures of Reservoirs 2 and 3, respectively, are obtained by slightly decreasing the corresponding values in Case A. Therefore, the effect of the decrease in the formation-fracturing pressure of a single reservoir on the intake percentage and the risk of each reservoir can be investigated by comparing the results from Cases A to those of 2C. In Cases D and E, the fracturing pressures of all the reservoirs increased by 0.5 and 1.0 MPa, respectively. Using these results, the effects of a synchronous increase in the formation-fracturing pressures can be studied. The percentages of the inflows and the risks of each reservoir in these five cases are illustrated in Figures 5 and 6, respectively.

The formation-fracturing pressure does not influence the CO2 pressure or injectivity. Therefore, modifying the formation-fracturing pressures does not alter the proportion of CO2 that flows into each reservoir. As shown in Figure 5, the majority of CO2 enters the upper two reservoirs. In contrast, variations of the formation-fracturing pressure considerably affect the risks of the reservoirs (Figure 6) and can cause the CR to shift. A comparison of the results of Cases A, B, and C indicates that when the fracturing pressure of only one reservoir formation is reduced, the risk of this reservoir immediately increases, such that this reservoir will become the CR when the formation-fracturing pressure decreases to a certain value. Furthermore, an identical increment in the formation-fracturing pressures of all the reservoirs should improve the safety of the reservoirs, and the reservoir with the smallest fracturing pressure will exhibit the largest decrease in risk. Similarly, increasing all the reservoirs’ formation-fracturing pressures should decrease the risks of the reservoirs without altering the position of the CR.

4.2.3. Influence of the Initial Formation Pressure

In engineering practices, the initial formation pressure is used to determine the lower limit of the injection pressure of the reservoir. Because the injection pressure is typically considerably higher than the initial formation pressure, variations in the initial formation pressure will not affect the reservoir risk directly but will affect the intake percentage of each reservoir. As discussed above, the uppermost reservoir is the CR for the original set of parameters used in this study. Because the pressure at this reservoir can be calculated directly from the wellhead conditions and is not related to the proportion of the CO2 inflow, the uppermost reservoir will remain the CR with the use of the original set of parameters. To further investigate these characteristics, another base case in which the CR is not the uppermost reservoir must be constructed, and the initial case can then be used as a reference to investigate the effects of varying the initial formation pressure.

Because the formation-fracturing pressure can affect the position of the CR considerably, directly changing its value allows the construction of a reference case. Therefore, two scenarios (i.e., S-1 and S-2) with two sets of formation-fracturing pressures are presented, and the settings of the formation-fracturing pressures of S-1 and S-2 are the same as those of Cases A and C, respectively, as shown in Table 4. Five different initial formation pressures for each scenario are provided in Table 5.

As noted above, the proportion of CO2 that flows into each reservoir is not related to the formation-fracturing pressure. The distributions of the flow rates among the reservoirs in S-1 and S-2 shown in Figure 7(a) are identical. Overall, when the initial formation pressure of a layer increases, the injection proportion of this layer decreases, whereas the flow rate that enters the layer beneath also increases. The inflow proportion for the uppermost reservoir changes only when its initial formation pressure changes. However, the initial formation pressures of the first and second reservoirs affect the inflow proportion of the second reservoir. The remaining inflow proportions can be deduced via an analogy; that is, the percentage of CO2 that enters a reservoir is related to the variation of the initial formation pressure of the reservoir itself and its upper reservoirs.

As the initial formation pressure increases, the risks of the reservoirs below the uppermost reservoir also increase, whereas the risk of the uppermost reservoir remains constant. Moreover, the initial formation pressure has a greater effect on the lower reservoirs than on the upper reservoirs, as shown in Figures 7(b) and 7(c). A comparison between the five cases of the two scenarios indicates that changes in the initial formation pressures have greater effects on the risks of the upper reservoirs than on the risks of the lower reservoirs. Reservoirs 1 and 3 are the CRs in S-1 and S-2, respectively. Thus, changes in the initial formation pressures also influence the position of the CR. In S-1, when the initial formation pressure of Reservoir 1 continues to increase, the risks of the lower reservoirs increase, and both Reservoirs 2 and 3 have the potential to become the CR. In S-2, when the initial formation pressure of Reservoir 1 or Reservoir 2 decreases to a certain value, then the risk of Reservoir 3 decreases and Reservoir 1 becomes the CR. Thus, when the first layer is not the CR, the risk of the CR can be reduced by decreasing the initial formation pressure of the upper reservoirs, which can help improve the safety and injectivity of the target site.

4.2.4. Influences of the Thickness, Permeability, and Porosity of the Reservoir

The thickness, permeability, and porosity of the reservoirs affect the capacity and injectivity of the reservoirs, which then affect the inflow and CO2 pressure in the wellbore, which is similar to the role of the initial formation pressure. Therefore, in this section, the aforementioned scenarios, S-1 and S-2, are used as the working models. The cases in Tables 68 were designed to investigate the effects of these three parameters, respectively. Changing the parameters of the last reservoir alone does not affect any of the results for the entire wellbore; hence, only the parameter sensitivities of the first three reservoirs are investigated.

Figure 8(a) shows the proportion of CO2 injected into each reservoir, considering the different reservoir thicknesses. The results in Figures 7(a) and 8(a) illustrate that as the thickness of a reservoir increases, the CO2 inflows into its upper reservoirs remain unchanged, whereas those of the reservoir itself and its lower reservoirs increase and decrease, respectively. As shown in Figures 8(b) and 8(c), an increase in the thickness of a single reservoir can effectively reduce the risks of its lower reservoirs without affecting those of its upper reservoirs, which is similar to the decreasing effect of the initial formation pressure. Therefore, the risk of the first reservoir is always fixed, regardless of whether the thickness of the reservoir varies. Moreover, increasing the thicknesses of the upper reservoirs to reduce the reservoir risk is significantly more effective than increasing the thicknesses of the lower reservoirs. In Cases B and C of S-2, the CR shifts from Reservoir 3 to Reservoir 1, which can effectively improve the safety of the wellbore and verify the conclusion that the CR might shift with changes in the formation-fracturing pressure. Additionally, the results for Cases B and C are the same in each of the three subpanels of Figure 8, which does not mean that doubling the thickness of Reservoir 1 produces the same effect as scaling up the thicknesses of all the reservoirs at the same time. This phenomenon occurs because the flow rate into Reservoir 1 accounts for 94% of the total flow after scaling up its thickness by one under the given injection conditions, and the remaining flow rate does not satisfy the capacity of the lower reservoir. In other words, when the total injection flow rate is adequate, the differences between the results of Cases B and C will clearly manifest themselves.

The scheme for generating the cases used to examine the reservoir permeability is the same as that used for the reservoir thicknesses. As shown in Figure 9, the influence of the permeability of the reservoir on the distribution of the flow rates and the risks of the reservoirs is nearly the same as that of the thickness of the reservoir, as shown in Figure 8, because the inflow of a reservoir, as expressed by (2), is proportional to the thickness and permeability of the reservoir. Therefore, we do not describe Figure 9 in further detail. However, subtle differences do exist in the flow-rate distributions and risks of the results shown in Figures 8 and 9 for the two scenarios. In short, the permeability has less significant effects than the thickness of the reservoirs.

Figure 10 shows the results when considering different reservoir porosities. The distributions of the flow rates and risks of the reservoirs for each case agree well with each other. Thus, the effect of the variations in porosity can be neglected.

5. Conclusions

This paper presented a systematic study on the key issues that arise from multiple CO2 injections during CCS projects, namely, the distributions of the flow rates among the reservoirs, the CRs, and their influencing factors. Understanding these issues is essential for both pressure design and field operations at wellheads. The calculation methods and programs were based on those in our previous work. Some specific conclusions of the study are given as follows.

(1) Assuming that all the reservoirs had identical properties, including thickness, porosity, and permeability, the injected fluid was mainly distributed in the upper reservoirs, and less CO2 entered the lower reservoirs. This trend mainly occurred because deeper reservoirs have greater formation pressure; therefore, high formation pressures require higher injection pressures to maintain injectivity.

(2) The risk of the reservoir was mainly controlled by the injection pressure and the formation-fracturing pressure and was minimally related to the inflow of CO2. Therefore, a reservoir could still be the CR even when the flow rate of CO2 into that reservoir is close to zero.

(3) The injection parameters at the wellhead considerably affected the flow-rate distributions and risks of the reservoirs. More CO2 entered the upper reservoirs than the lower reservoirs with increasing injection pressure (or decreasing injection rate). Moreover, the risks of all the reservoirs increased, and the risk-increasing trend of the upper reservoirs was notably greater than that of the lower reservoirs. Because the extent of the variation in each reservoir was different, the CR could shift. This result demonstrates that the injection pressure controls the overall risk level of the reservoir and the flow-rate distribution among the reservoirs can control the position of the CR. The CO2 injection temperature had a similar effect as the injection flow rate, but the fluid injection temperature had no effect on the position of the CR.

(4) Although the formation-fracturing pressure does not affect the flow-rate distribution, it is one of the key factors for determining the reservoir risk and CR. Decreasing the formation-fracturing pressure increased the reservoir risk and cause the CR to shift. When the formation-fracturing pressures of all the reservoirs were decreased by the same factor, the risk of each reservoir increased equally without changing the CR. Under a constant flow rate, variations in the reservoirs’ porosities had a minimal effect on the flow-rate distribution and the reservoir risk. The thickness and permeability had similar effects on the flow-rate distributions and reservoir risks. For a specific reservoir, the flow rate increased with increasing thickness or permeability, and the flow rates of the lower reservoirs simultaneously decreased, whereas the flow rates of the upper reservoirs remained unchanged. Furthermore, only the risks of the lower reservoirs decreased with increases in these two parameters; therefore, the CR can shift to another reservoir. The initial formation pressure had the opposite effect as the reservoir thickness. Therefore, certain reservoirs can be selectively reformed according to their risks to decrease these risks and adjust the position of the CR, which can effectively improve the injection safety and utilization efficiency of the reservoirs.

(5) The injection flow rates notably influenced the flow-rate distribution. The position of the CR and injection flow rate may vary in actual CO2 injection processes. The design of the maximum allowable wellhead injection pressure should consider multiple flow-rate scenarios to reduce uncertainties and to improve the reliability of the project.

(6) The results obtained from our investigations of the influence of the injection parameters on the flow-rate distributions and reservoir risks in the case study of the Shenhua CCS demonstration project can be used to guide subsequent injection operations. This research was based on a vertical wellbore that penetrated a number of reservoirs. Although the results for horizontal and inclined wells may be similar, further research on these topics is still required.

Conflicts of Interest

The authors declare that they have no conflicts of interest.


This work was sponsored by the International Science & Technology Cooperation Program of China (S2016G9005) under the cooperation framework of the US-China Clean Energy Research Center (CERC).