Evaluating Reservoir Risks and Their Influencing Factors during CO 2 Injection into Multilayered Reservoirs

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 increasemore 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.


Introduction
Carbon capture and storage (CCS) is widely recognized as an effective approach for greatly reducing CO 2 levels in the atmosphere [1][2][3][4][5].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 largescale CO 2 injection, multiple reservoirs are employed for simultaneous injection to achieve a preset injection target amount of CO 2 (see [13][14][15][16]).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 CO 2 injection through deep wellbores requires an in-depth understanding of the reservoir performance and safety and their influencing factors in terms of CO 2 injection.
The injection of CO 2 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 2 Geofluids 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 CO 2 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 CO 2 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.

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 threedimensional 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 CO 2 pressure can be determined when the injection pressure at the wellhead is known; then, the fluid pressure  +1 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/s 2 ),  0 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 CO 2 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 CO 2 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 CO 2 -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 CO 2 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  0 is the initial formation pressure of the reservoir (Pa);   is the injection pressure on the wellface of the reservoir (Pa);  0 is the maximum influence radius of the flow in the reservoir (m);   and  max are the radii of the CO 2 plume at the bottom and top of the reservoir, respectively (m);  is the absolute permeability of the reservoir (m 2 );  is the thickness of the reservoir (m);  is the density of the CO 2 (kg/m 3 ); and   and   are the mobility of CO 2 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 CO 2 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),  0 , and   , and  max , are used to describe the distributions of CO 2 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- (m 3 );  is the reservoir porosity;   is the pore compressibility (m 2 /N);   is the compressibility of brine (m 2 /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 CO 2 injections.The wellbore constraints include the pressure and flowrate 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 as where   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).

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, 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 CO 2 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 CO 2 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 CO 2 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 CO 2 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.

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 CO 2 , 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.

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 1A to 1D), injection rates (Cases 1D to 1F), injection temperatures (Cases 1G to 1J), and combinations of the first two factors (Cases 1C, 1E, and 1G) on the percentages of the CO 2 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 CO 2 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 CO 2 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 CO 2 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 CO 2 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 CO 2 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 CO 2 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 1D to Case 1F, 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 CO 2 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 CO 2 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 1G, 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.

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 CO 2 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.

Influence of the Formation-Fracturing Pressure.
As presented above, the concept of the CR depends on the   /[  ] ratio.[  ] is typically obtained from the formationfracturing 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 2B and 2C, the formation-fracturing pressures of Reservoirs 2 and 3, respectively, are obtained by slightly decreasing the corresponding values in Case 2A.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 2A to those of 2C.In Cases 2D and 2E, 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 CO 2 pressure or injectivity.Therefore, modifying the formation-fracturing pressures does not alter the proportion of CO 2 that flows into each reservoir.As shown in Figure 5, the majority of CO 2 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 2A, 2B, and 2C 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' formationfracturing pressures should decrease the risks of the reservoirs without altering the position of the CR.

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 CO 2 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 formationfracturing pressures are presented, and the settings of the formation-fracturing pressures of S-1 and S-2 are the same as those of Cases 3A and 3C, 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 CO 2 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 CO 2 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 CO 2 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 6-8 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 CO 2 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 CO 2 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 4B and 4C 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 4B and 4C 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 4B and 4C 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.

Conclusions
This paper presented a systematic study on the key issues that arise from multiple CO 2 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 CO 2 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 CO 2 .Therefore, a reservoir could still be the CR even when the flow rate of CO 2 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 CO 2 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 CO 2 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 flowrate distribution.The position of the CR and injection flow rate may vary in actual CO 2 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.

Figure 3 :
Figure 3: Distributions of the flow rates that correspond to the different injection parameters at the wellhead of the Shenhua CCS project.(a), (b), and (c) are the results for different injection pressures, flow rates, and temperatures, respectively.(d) shows the investigation of the influences of simultaneous changes of the injection pressures and flow rate.

Figure 4 :
Figure 4: Risk factors of each reservoir which correspond to the different injection parameters at the wellhead of the Shenhua CCS project.(a), (b), and (c) are the results for the different injection pressures, flow rates, and temperatures, respectively.(d) Investigation of the influences of simultaneous changes of the injection pressures and flow rates.

Figure 6 :
Figure 6: Variations in the risk factors of the reservoirs at different formation-fracturing pressures for the conceptual engineering model.

Figure 7 :
Figure 7: Proportions of the inflows (a) and risk factors of each reservoir under different initial formation pressures for the two scenarios of the conceptual engineering model.(b) and (c) are the risk-factor results for Scenarios S-1 and S-2, respectively.

Figure 8 :
Figure 8: Proportions of the inflows (a) and risk factors of each reservoir for different reservoir thicknesses for the two scenarios of the conceptual engineering model.(b) and (c) are the risk-factor results for Scenarios S-1 and S-2, respectively.

Figure 9 :
Figure 9: Proportions of the inflows (a) and risk factors of each reservoir for different reservoir permeabilities for the two scenarios of the conceptual engineering model.(b) and (c) are the risk-factor results for Scenarios S-1 and S-2, respectively.

Figure 10 :
Figure 10: Proportions of the inflows (a) and risk factors of each reservoir for different reservoir porosities for the two scenarios of the conceptual engineering model.(b) and (c) are the risk-factor results for Scenarios S-1 and S-2, respectively.

Table 1 :
Parameters of the reservoir-caprock units of the Shenhua CCS project.
ki ] Figure 2: Risk factors of each reservoir for the base case of the Shenhua CCS project.

Table 2 :
Injection parameters for the described cases of the Shenhua CCS project.

Table 3 :
Characteristic parameters of the reservoirs in the conceptual engineering model.
Figure 5: Flow-rate distribution of a conceptual model with the corresponding parameters of the formations in Table3; the injection pressure, flow, and temperature are 10 MPa, 4 kg/s, and −5 ∘ C, respectively.

Table 5 :
Five cases with different initial formation pressures (MPa).

Table 6 :
Five cases with different reservoir thicknesses (m).