A sensitivity analysis of factors affecting in geologic CO2 storage in the Ordos

Shi-xin Di, Yn-jio Dong, Feng Wng, Zhen-hn Xing, Pn Hu, Fu Yng

a School of Earth Sciences and Geospatial Information Engineering, Hunan University of Science and Technology, Xiangtan 411201, Chinab China Railway Times Architectural Design Institute Co.Ltd., Wuhu 241000, Chinac Key Laboratory of Coal Resources Exploration and Comprehensive Utilization, Ministry of Natural Resources, Xi’an 710021, China

Keywords:Geologic CO2 storage Influencing factors Sensitivity analysis Carbon neutrality Oil and gas exploration engineering Ordos Basin China

ABSTRACT To accelerate the achievement of China’s carbon neutrality goal and to study the factors affecting the geologic CO2 storage in the Ordos Basin, China’s National Key R&D Programs propose to select the Chang 6 oil reservoir of the Yanchang Formation in the Ordos Basin as the target reservoir to conduct the geologic carbon capture and storage (CCS) of 100000 t per year.By applying the basic theories of disciplines such as seepage mechanics, multiphase fluid mechanics, and computational fluid mechanics and quantifying the amounts of CO2 captured in gas and dissolved forms, this study investigated the effects of seven factors that influence the CO2 storage capacity of reservoirs, namely reservoir porosity,horizontal permeability, temperature, formation stress, the ratio of vertical to horizontal permeability,capillary pressure, and residual gas saturation.The results show that the sensitivity of the factors affecting the gas capture capacity of CO2 decreases in the order of formation stress, temperature, residual gas saturation, horizontal permeability, and porosity.Meanwhile, the sensitivity of the factors affecting the dissolution capture capacity of CO2 decreases in the order of formation stress, residual gas saturation,temperature, horizontal permeability, and porosity.The sensitivity of the influencing factors can serve as the basis for carrying out a reasonable assessment of sites for future CO2 storage areas and for optimizing the design of existing CO2 storage areas.The sensitivity analysis of the influencing factors will provide basic data and technical support for implementing geologic CO2 storage and will assist in improving geologic CO2 storage technologies to achieve China’s carbon neutralization goal.

1.Introduction

Carbon neutrality has become an international effort and a commitment of the Chinese government.However, China’s carbon neutrality goal is facing many challenges, such as high total carbon emissions, short time of carbon emission reduction, many challenges in economic transformation and upgrading, and difficult energy system transformation.As the world’s largest developing country, China needs to take tougher measures to reduce carbon emissions in wider areas in a shorter time to achieve carbon neutrality by 2060 (Mallapaty S, 2020; Salvia M, 2021; Xu T et al., 2022).Developing geologic CO2storage technologies is a strategy for China to reduce CO2emissions and ensure energy security in the future, and it is also an important means to build ecological civilization and to achieve sustainable development (Bachu S,2000, 2003; Zeng RS et al., 2004; Jiang HY et al., 2007; Li GJ and Zhang J, 2011; Goodarzi S et al., 2012; Chen B et al.,2018; Li Q et al., 2019; Ma X et al., 2021; Wen DG et al.,2014).To select and optimize the design of CO2storage areas to ensure that the injected CO2is stored underground safely,effectively, and permanently, it is necessary to assess the sensitivity of the influencing factors for CO2storage in advance.

Many scholars have studied the influencing factors of geologic CO2storage.(Wiese B et al., 2010) conducted a local sensitivity analysis of parameters affecting the CO2injection rate of a single well and concluded that permeability has the greatest impact, followed by the injection pressure,while vertical anisotropy has a lower impact.Ershadnia B et al.(2020) use a transition-probability based approach to simulate heterogeneous systems with binary facies distributions and the resulting petrophysical properties at the field scale.Results reveal that for a given volume of injected CO2, shorter injection times yield higher total amounts of trapped CO2.In the study of the CO2migration law on different scales in Mount Simon aquifers in the Illinois Basin,Zhou QL et al.(2010) found that relatively low permeability and high intake pressure have strong impacts on CO2storage.Ke YB et al.(2012) conducted a numerical simulation study of CO2storage and injection in saline aquifers in the Jianghan Basin.The sensitivity analysis of salinity in this study showed that the increasingly high salinity of reservoirs corresponds to increasingly strong pressure accumulation.Zhao RR et al.(2012) conducted a simulation study and parameter sensitivity analysis of fluid migration in deep saline aquifers in the Songliao Basin after CO2injection.This study showed that residual gas saturation and the anisotropy of saline aquifers have the greatest influence on CO2dissolution capture,followed by salinity, while residual liquid saturation and initial CO2concentration have the lowest influence.Zheng F et al.(2014) analyzed the effects of parameters(hydrogeological parameters and salinity) of saline aquifers in the Subei Basin on CO2storage, migration, and leakage.The results of these studies have strongly promoted the study of geologic CO2storage.However, existing literature tends to conduct comparative analyses of partial influencing factors and does not provide comprehensive analyses of influencing factors of geologic CO2storage.

To achieve the carbon neutrality goal, geologic CO2storage technologies are predicted to reduce CO2emissions by 0.6×109-1.4×109t in 2050 and by 1.0×109-1.8 ×109t in 2060 according to the current technology development.Taking the Upper Triassic Yanchang Formation in the central part of the Shaanbei Slope in the Ordos Basin as a case study, this study investigated the effects of seven factors affecting the CO2storage capacity of reservoirs, namely reservoir porosity,horizontal permeability, temperature, formation stress, the ratio of vertical to horizontal permeability (kv/kh), capillary pressure, and residual gas saturation.The purpose of this study is to provide data basis and technical support for the design and implementation of industrial CO2storage projects in the Ordos Basin in the future.Moreover, this study presented the actions China has taken in the Ordos Basin to achieve the great carbon neutrality goal.

2.Regional geology

This study focuses on the Huaziping oil area on both sides of Hezhuanggou District, western Huaziping Town, Ansai County, Shaanxi Province.The Huaziping area is a part of the Xingzichuan Oil Production Plant of the Yanchang Oilfield.In terms of regional tectonics, it is a monoclinal structure inclining from east to west lying in the middle section of the Shaanbei Slope in the Ordos Basin.It is mostly flat but has multiple small nose-shaped structures locally.The strata in the area incline westward, with a dip angle of about 0.6° (Lei XL et al., 2008).

Based on its evolutionary history and its current structural morphology, the Ordos Basin is divided into six major structural units, namely the western margin thrust belt, the Tianhuan Depression, the Yishan Slope, the Jinxi Flexure Belt, the Yimeng Uplift, and the Weibei Uplift (Fig.1).

Based on the drilling data of the Huaziping area (Fig.2),the strata in the study area are divided into the Quaternary loess layer, and the Luohe, Anding, Zhiluo, Yan’an, Fuxian,and Yanchang Formations.Table 1 shows brief lithological descriptions of these strata.The oil reservoir in the Chang 6 Member of the Yanchang Formation serves as both the primary block and stratum developing in the Huaziping area and the principal reservoir for geologic CO2storage (Fig.3).The stable Chang 6 Member is mainly composed of gray and dark gray medium-fine-grained feldspar quartz sandstones and feldspar lithic sandstones (Li B, 2014; Zhang ZM, 2015;Gao ZD, 2017).With a porosity of 7%-12% (average: 9.4%)and permeability of 0.1×10-3-5×10-3μm (average: about 0.94×10-3μm2), this member is a low-porosity and lowpermeability reservoir (Zheng YR et al., 2011; Zhao XS et al.,2019).The formation pressure in the Chang 6 Member is 4.68-9.46 MPa, with an average of 7.08 MPa (Fig.4).The average temperature and the geothermal gradient of the member ’s oil reservoir are 47.0°C and 2.7°C/100 m,respectively (Zhang TJ et al., 2016).The overlying Chang 4 and 5 members show the development of mudstones and thus are favorable cap rocks.Therefore, factors such as reservoir burial depth, lithological characteristics, the sealing of regional cap rocks, and geological trap conditions suggest that this series of strata are qualified to serve as a site for geologic CO2storage.

Table 1.Geological and lithological characteristics of the Huaziping area.

Fig.1.a-tectonic framework in China showing the spatial location of the Ordos Basin (modified from Groves DI et al., 2020); b-structureoutline map of the Ordos Basin (modified from Li SX et al., 2017).

The primary sedimentary facies in the study area include the subaqueous distributary channel microfacies and the interdistributary bay microfacies deposited in the delta front subfacies.The subaqueous distributary channel microfacies consists of large composite sand bodies that are distributed in bands in the NE-SE direction.The interdistributary bay microfacies is composed of thin mudstones, silty mudstones,and argillaceous siltstones.The sand body thickness profile of the study area (Fig.5) shows that the sand bodies parallel to the river channels take the shape of bands and have good connectivity, whereas those perpendicular to river channels are poorly connected and are flat in the upper part and convex in the lower part.

To determine the lithological characteristics of strata in the study area, five core samples were collected from production well H303-14 (Fig.6) and were sent for testing and analysis.These core samples were taken at depths of 1476.59-1478.37 m, 1450.77-1452.52 m, 1430.97-1432.72 m, 1421.43-1423.78 m, and 1418.08-1419.78 m,corresponding to the Chang 64, Chang 63, Chang 63, Chang 61-2, and Chang 61-1submembers of the Chang 6 oil reservoir,respectively.Mineral content analysis and conventional pore permeability tests were carried out on these core samples.

The types and contents of the minerals in these core samples are shown in Table 2.The main minerals include quartz, potassium feldspar, plagioclase, calcite, dolomite, and clay minerals.The core sample taken from the Chang 61-1Submember has a higher clay mineral content and a lower feldspar content compared to other four core samples,indicating that the Chang 61-1Submember has a higher argillaceous content.The core samples from other submembers have high feldspar content.Despite slightly different contents of various minerals, the overall lithology of the Chang 6 oil reservoir is dominated by sandstones.

The measured porosity and permeability of the core samples are shown in Table 3.Core samples Nos.2-5 are roughly similar to sandstones in the Chang 6 oil reservoir,while the porosity and permeability of core sample No.1 are more similar to those of the upper Chang 4 and Chang 5 members.

3.Simulation principle

The numerical simulation of geologic CO2storage in depleted low-porosity and low-permeability oil and gas reservoirs aims to study the gas-liquid-solid multiphase flow consisting of multi-component mixture H2O-NaCl-CO2.The fluid migration in depleted oil and gas reservoirs involves constant changes in temperature, salinity, pressure, and heat and is even accompanied by chemical and mechanical changes (Pruess K et al., 1999; Xu T, 2001, 2006; Ling LL et al., 2013; Zhang EY, 2012).

In depleted oil and gas reservoirs, supercritical CO2occurs in a single non-wetting phase and forms a multicomponent multiphase flow system together with the saline solution in the reservoirs.In this system, the flow of each component satisfies the law of conservation of mass.The system does not consider temperature changes.Its massconservation equation consists of a mass change item, a flow item, and a source-sink item.The change in the fluid mass in a volume is equal to the sum of the fluid mass flowing into the volume from the volume surface and the fluid mass of the source-sink item (Pruess K et al., 1999; Zhang MY et al.,2017; Zhang ZX et al., 2018).

Fig.2.Geological profile of the Huaziping area.

where,Mκrepresents the cumulative amount of mass per unit volume in an area;Fkis the mass flux, including convection flux and diffusion amount;qkis the mass source or energy per unit volume in an area;nis the internal unit normal vector; superscriptkrepresents the material component, andtrepresents the time.

The total mass of all components in a porous media can be expressed using the following equation:

where, φ is the porosity; β represents the phase state;Sβ,ρβ,andXβκrepresents the saturation, fluid density, and component mass fraction of each phase, respectively, and κ=w,i,gis the water, salt, and gas component.

In the process of CO2displacing saline water, it is believed that gas and water phases in the reservoirs obey the multiphase Darcy law (Pruess K et al., 1999):

where,krepresents the absolute permeability coefficient of a porous media;krβrepresents the relative permeability coefficient of phase β; ρβand µβrepresents the density and viscosity of phase β, respectively;Pβrepresents the pressure of phase β, andgrepresents the acceleration of gravity.

Fkin equation (1), which represents the mass flux of eachphase, can be deduced in the multiphase flow form of Darcy’s law:

Fig.3.The histogram of the Chang 6 oil reservoir.

4.Design and establishment of the geological model

Based on the typical lithological characteristics and physical parameters of the Chang 6 Member in the Huaziping oil area, the area to be modeled was generalized into a cylinder, and a geological model of a 140 m thick twodimensional homogeneous sandstone system was established(Fig.7).RZ2D was adopted for mesh generation of the model.This model was divided into seven simulated reservoirs with the same thickness (20 m) in the vertical direction.In the horizontal direction, the model was 100000 m long and was divided into 130 radial grids in total.The radius of the injection wellbore was set to 0.3 m.The hydrogeological parameters of the homogeneous model are listed in Table 4.For the sandstones in the simulated reservoirs, the average porosity was set to 12% and the average horizontal permeability was set to 3 mD.Given that the general ratio of horizontal to vertical permeability is 10∶1, the average vertical permeability was set to 0.3 mD.The vicinity of the injection well was the critical area for CO2diffusion and transport, and the outermost cells was set as boundaries with an infinite volume.CO2was injected into the reservoirs from the bottom of the injection well, with an injection depth of 20 m.It was assumed that the CO2injection was an isothermal process.The initial conditions of the model were calculated using the average formation pressure coefficient and the empirical formula of the geothermal gradient in the study area (Li C et al., 2011).

5.Similar simulation of geologic CO2 storage

To conduct a parallel comparison of the effects of various physical parameters of reservoirs on CO2capture in gas and dissolved forms, the standard, reduced, and magnified values of seven influencing factors (i.e., porosity, horizontalpermeability, temperature, formation stress, the ratio of vertical to horizontal permeability, capillary pressure, and residual gas saturation) were set (Table 5).In the simulation,the capillary pressure and residual gas saturation were set at empirical values.The CO2injection rate and injection time were set at 0.15 kg/s and 100 years, respectively, in order to keep the total CO2injection amount in the reservoirs consistent.

Table 2.Mineral types and contents of core samples.

Table 3.Sample rock porosity and permeability.

Table 4.Hydrogeological parameters in geological model.

Table 5.Parameters for sensitivity analysis of the model.

Fig.4.Profile of the oil reservoir of the Huaziping area (modified from Wang F et al., 2020).

Fig.5.Sand body thickness profile of the Huaziping area.

Numerical simulation was carried out according to the parameters selected in the table, taking porosity as an example, and the results are shown in Fig.8.The last six sets of simulation results were obtained by the same method.

5.1.Porosity

Porosity greatly affected the distribution of gas saturationSgand the concentration of dissolved CO2XCO2(aq) after 100 years of CO2storage.The spreading scope of the CO2diffusion halo decreased as the porosity increased.This is because a high porosity indicates a high water content in strata, which is conducive to the dissolution and capture of CO2.However, this limited the migration and distribution scope of the CO2diffusion halo, thereby reducing the chance of wide contact between CO2and water.

Fig.6.The core samples of well H303-14, Chang 6, Huaziping.

5.2.Horizontal permeability

Permeability is one of the influencing factors of geologic CO2storage.Its effects on CO2saturation were difficult to identify, but it had distinct effects on the concentration of CO2dissolved in water XCO2(aq).Theoretically, the CO2diffusion halo is distributed in an increasingly wide area as the horizontal permeability increases.Thus, CO2can contact water more widely.Therefore, the increase in the horizontal permeability is beneficial to the dissolution and capture of CO2, but it decreases the concentration of CO2dissolved in water.

5.3.Temperature

The distribution of the CO2saturation and the concentration of CO2dissolved in water after 100 years of CO2storage at 45°C, 47°C, and 70°C indicated that the migration scope of CO2increased as the temperature increased.This occurred because the CO2solubility decreased as the temperature increased, causing more CO2to continue to migrate in the form of gas at a higher temperature instead of dissolving in formation water.

5.4.Stratum stress

As revealed by the distribution of the CO2saturation and the concentration of CO2dissolved in water after 100 years of CO2storage under formation stress of 4.68 MPa, 7.08 MPa,and 9.46 MPa, the spreading scope of the CO2diffusion halo increased under formation stress of 4.68 MPa compared to that under formation stress of 9.46 MPa.This occurred because more CO2vertically migrated and then accumulated and diffused at the bottom of cap rocks under lower formation stress.

5.5.Ratio of vertical to horizontal permeability (kv/kh)

Reservoir heterogeneity is an important influencing factor in CO2migration, distribution, and storage, and the ratio of vertical to horizontal permeability (kv/kh) is an important parameter reflecting reservoir heterogeneity.This study investigated the effects of the ratio of vertical to horizontal permeability on CO2storage by changing vertical permeability.

According to simulated result, a decrease in the ratio of vertical to horizontal permeability corresponded to an increase in the vertical and horizontal migration capacity of CO2and a decrease in the direction scope of CO2gas.Whenkv/kh=10∶1, CO2rapidly migrated in the vertical direction in homogeneous reservoirs while migrating and diffusing in the horizontal direction.As a result, massive amounts of CO2migrated to cap rocks in a short time and laterally migrated under the cap rocks due to the blockage of the cap rocks.Since CO2migrated and diffused in a large scope, the contact area of CO2with fresh formation water grew.This increased contact was conducive to the conversion of free CO2into other stable CO2storage forms.Whenkv/kh= 1∶10, CO2migrated more rapidly in the vertical direction.In this case,CO2first migrated in the vertical direction along the injection well and then horizontally migrated along the bottom of cap rocks.As a result, the CO2gas had only a small influencing scope and a limited contact area with formation water.This was not conducive to the conversion of free CO2into other storage forms such as irreducible or dissolved states.

With a decrease in the ratio of vertical to horizontal permeability, the horizontal dissolution scope of CO2gradually increased, the scope of the dissolved gas pocket of CO2constantly shrank, and the transition zone with a decreasing CO2molar mass fraction decreased.Whenkv/kh=10∶1, highly dissolved CO2distributed in the form of an inverted triangle occurred in large areas, and more CO2

Fig.7.A diagram of the CO2 perfusion radial well flow model(modified after Wang F et al., 2020).

was stored in strata in a relatively stable state.This is conducive to the dissolution of CO2.In contrast, whenkv/kh=1∶10, highly dissolved CO2distributed in the form of an inverted triangle occurred in only small areas, and there was less CO2stored in strata in a relatively stable state.This is inconducive to the conversion of free CO2into a dissolved form.

5.6.Capillary pressure

Capillary pressure plays an important role in geologic CO2storage, especially in the storage of irreducible gas.Therefore,a comparative study was conducted to determine the effects of capillary pressure on geologic CO2storage under the presence of capillary and in the absence of capillary pressure (i.e.,minimum capillary pressure model).

The grids right above the top of the injection well were selected to observe the effects of capillary pressure on CO2dissolution in water and CO2gas saturation of the reservoirs.As shown in Fig.9a, capillary pressure has almost no effect on CO2dissolution in water.The reason is as follows.The sequestration amount of dissolved CO2is mainly affected by the salinity of formation water and temperature and pressure conditions.Therefore, as an important influencing factor in the sequestration of irreducible gas, capillary pressure has no direct influence on the sequestration of dissolved CO2.Fig.9b shows the CO2saturation curves of the reservoirs.After CO2injection stopped, the CO2saturation curve under the presence of capillary pressure became stable, while the CO2saturation curve of the minimum capillary force model slowly decreased until it was finally relatively stable.The reasons are as follows.After CO2injection stopped, CO2was sequestrated in stratum pores as irreducible gas under the action of capillary forces.In contrast, the minimum capillary pressure model had a relatively weak capacity to sequestrate irreducible gas, and thus the corresponding CO2saturation decreased.

5.7.Relative permeability

Studies (Wei L and Saaf F, 2009; Wang L, 2014) have revealed that the relative permeability of water has a small effect on the migration and distribution scopes of CO2and the amount of CO2storage.Therefore, this paper focuses on the effects of the relative permeability of gas on CO2storage.The relative permeability of gas is one of the important factors affecting CO2storage capacity, and a slight change in the morphology or value of the relative permeability curve may significantly affect the simulation results.However, there are no relevant data about the relative permeability curve and the capillary pressure curve of the study area since it is difficult to obtain capillary pressure and its relevant parameters on site.Given this, this study referred to the capillary pressure setting adopted in previous studies (Alkan H et al., 2010).

With a decrease in the residual gas saturation (Sgr), CO2stored in an irreducible state decreased, the permeability (Krg)of gas increased, the migration and distribution capacity ofCO2increased, and the influencing scope of CO2constantly expanded accordingly.As shown in simulated result, when the residual gas saturation (Sgr) was as small as 0.01, highly dissolved CO2distributed in the shape of an inverted triangle had a small area, and a small amount of CO2was stored in strata in a relatively stable state.In this case, CO2showed a high migration and diffusion capacity and can diffuse up to a maximum distance of 504 m.When the residual gas saturation(Sgr) was 0.05, the highly dissolved CO2had a larger area compared to the case of the residual gas saturation (Sgr) of 0.01, and CO2diffused up to a maximum distance of 446 m.When the residual gas saturation (Sgr) was 0.1, the highly dissolved CO2distributed in the form of an inverted triangle has the largest area, and a large amount of CO2was stored in strata in a relatively stable state.In this case, CO2showed the lowest migration and diffusion capacity and can diffuse up to a maximum of 395 m.In the model of this study, the total storage amount of irreducible CO2increased with an increase in the residual gas saturation.In the case of low residual gas saturation of CO2, the relative permeability of gas was high,and the CO2had high migration and diffusion capacity.As a result, more CO2was dissolved in the process of contact with formation water in a large scope.

6.Comparative analysis of the sensitivity of physical property parameters of CO2 reservoirs

6.1.Changes in capture amount under different capture mechanisms

For hydrodynamic capture, CO2exists in the form of gas in reservoir pores.Dissolution capture refers to the dissolution of CO2in formation water.Mineral capture means that CO2is sequestrated through precipitation as carbonate minerals in the process of CO2-water-rock interactions.The calculationformulas of capture amount under different capture mechanisms are as follows (Yang GD, 2015):

Fig.8.Distribution of saturation and the concentration of CO2 dissolved in water after 100 years of CO2 storage.

Fig.9.a-curve of dissolved CO2 in water; b-curve of CO2 gas saturation.

where,MCO2gas,MCO2aqueous, andMCO2mineralare the total amount of CO2captured in the form of CO2gas, dissolved CO2, carbonate minerals, kg;nmaxis the total number of grids;Vnis the volume of thenth grid, m3; φnis the porosity of thenth grid, dimensionless;SgnandSlnare gas saturation and liquid saturation, dimensionless;DgnandDlnare gas density and liquid density, kg/m3;Xco2is the mass fraction of dissolved CO2in liquid, dimensionless;SMco2is the amount of captured carbonate minerals, kg/m3.

6.2.Analysis and comparison of effects of various parameters on CO2 capture amount

The total gas and dissolution capture amounts of CO2after 100 years of CO2storage can be calculated using formulas 5 and 6, as shown in Tables 6 and 7.According to these tables,with a gradual increase in the temperature, formation stress,capillary pressure, and residual gas saturation, the total gas capture amount of CO2increased, while the total dissolution capture amount of CO2decreased.Moreover, with a gradual increase in porosity, horizontal permeability, and the ratio of vertical to horizontal permeability (kv/kh), the total gas capture amount of CO2decreased, while the total dissolution capture amount of CO2increased.

6.3.Sensitivity analysis

The degree of the influence of a parameter on the model can be expressed by the sensitivity of the model to this parameter.Sensitivity refers to the degree of influence of the change in a factor on other factors.To quantify the sensitivity of the capture amount of CO2to physical property parameters of reservoirs, the sensitivity of each physical property parameter of the reservoirs was calculated using the following formula (Xue YQ and Xie CH, 2007):

where, subscriptiis the time observation point;Hiis the CO2capture amountMCO2calculated at theith time observation point when thekth reservoir physical property parameter is taken as αk(kg); αkis the value of thekth physical property parameter; ∆ αkis the increment of αk; βi,kis the sensitivity ofHito thekth reservoir physical parameter at theith time observation point.

Because the physical parameters of reservoirs have different dimensions, their sensitivity cannot be compared directly.Therefore, the above formula was transformed into a dimensionless form:

Since sensitivity parameter βi,kchanges greatly with time,Hiof thekth reservoir physical parameter calculated at theith time observation point was taken as the average valueHof thekth reservoir physical parameter over the entire time interval.

The sensitivity of gas and dissolution capture amounts to various physical parameters can be calculated using formulas 9 and 10, as shown in Tables 5 and 6.Since these formulas do not apply to the sensitivity calculation of capillary pressure and the ratio of vertical to horizontal permeability (kv/kh),these two parameters were excluded from the parallel comparison.

Table 6.Effects of various parameters on gas capture amount of CO2.

Table 7.Effects of various parameters on dissolution capture amount of CO2.

As can be seen from Table 8, the standard sensitivity (i.e.,the sensitivity calculated using standard parameter values) of factors affecting the gas capture capacity of CO2decreased in the order of formation stress, temperature, residual gas saturation, horizontal permeability, and porosity (0.573 >0.072 > 0.067 > 0.052 > 0.024).Moreover, the gas capture amount of CO2had a positive correlation with the temperature, formation stress, and residual gas saturation but had a negative correlation with the porosity and horizontal permeability.According to Table 9, the sensitivity of factors affecting the dissolution capture capacity of CO2decreased in the order of formation stress, residual gas saturation,temperature, horizontal permeability, and porosity (0.763 >0.09 > 0.086 > 0.072 > 0.032).Furthermore, the dissolution capture capacity CO2was positively correlated with the porosity and horizontal permeability but was negatively correlated with the temperature, formation stress, and residual gas saturation.

7.Carbon neutrality-measures in the Ordos Basin

As the second largest sedimentary basin in China, the Ordos Basin enjoys extremely rich natural resources and is the most important energy production and supply base in China.The proven resources of coal, natural gas, and coalbed methane in the basin all rank first in China, and the proven reserves of petroleum in the basin rank fourth in the country.However, CO2emissions in the basin have sharply increased and account for a large proportion of China’s total CO2emissions due to the rapid development of the energy,chemical industry, and coal-fired power in the basin, as well as the CO2emissions from industries such as petroleum,natural gas, cement, and iron and steel.Therefore,decarbonization is crucial in the energy field, and further actions are urgently needed to reduce carbon emissions and mitigate the impacts of carbon emissions on climate change.

The sensitivity analysis of factors affecting geologic CO2storage conducted in this study is practically significant for achieving carbon neutrality, constructing ecological civilization, and obtaining sustainable development in the Ordos Basin.Based on the sensitivity analysis of factors affecting geologic CO2storage, the CO2storage areas in the Ordos Basin can be optimized.For example, the sensitivity analysis revealed that formation stress has more prominent effects in the process of CO2gas capture.Table 10 shows the changes in the gas capture amount of CO2induced by different formation stress in this study.

Table 8.Sensitivity of various parameters to gas capture of CO2.

Table 9.Sensitivity of various parameters of dissolution capture of CO2.

Table 10.Analysis of gas capture amount CO2 under different formation stress.

As can be seen, the gas capture amount of CO2increased by 40×106kg when the formation stress increased from 4.68 MPa to 7.08 MPa and increased by 64×106kg when formation stress increased from 7.08 MPa to 9.46 MPa.This indicates that the gas capture amount of CO2increased exponentially with an increase in the formation stress.In this case, CO2gas capture at the storage site can be further improved by increasing formation stress, which can be obtained by appropriately increasing CO2injection amount on the premise of no leakage of cap rocks, adequate storage capacity of the reservoirs, and given injection pressure or rate of the storage site.

Meanwhile, based on the analysis of the factors affecting geologic CO2storage, a reasonable assessment can be made from the following aspects for the CO2storage site in the Ordos Basin: (1) Whether the potential storage site can provide the required storage capacity; (2) Whether the potential storage site can provide given injection pressure or rate; (3) Whether the sequestrated CO2can be safely and permanently trapped at the potential storage site.

8.Conclusions

A study on geologic CO2storage was conducted in the Ordos Basin, China, aiming to achieve the goal of carbon neutrality.Based on the geological conditions of the lowpermeability oil reservoirs in the Chang 6 Member of the Upper Triassic Yanchang Formation in the central part of the northern Shaanxi Slope in the Ordos Basin, this study systematically explored the effects of seven influencing factors (i.e., reservoir porosity, horizontal permeability,temperature, formation stress, the ratio of vertical to horizontal permeability, capillary pressure, and residual gas saturation) on the main capture mechanisms (including gas capture and dissolution capture) after CO2injection into the reservoirs.The primary understandings obtained are as follows.

(i) The Ordos Basin is a large-scale regional energy and chemical industry base in China, where a large amount of CO2is intensively emitted.Based on available data and studies, the Ordos Basin has ideal reservoir - cap rockassemblage and potential strata for geologic CO2storage.Meanwhile, there are relatively complete well networks and sufficient basic data of this basin.The Huaziping area in the Ordos Basin is close to Yulin and Jingbian areas, where many large-scale energies and chemical enterprises lie nearby.The CO2emitted during the industrial production of these enterprises can be used as sufficient gas sources ofin-situCO2storage.Therefore, the Huaziping area is an excellent place for geologic CO2storage.

(ii) For the study area, the standard sensitivity of factors affecting the gas capture capacity of CO2decrease in the order of formation stress, temperature, residual gas saturation,horizontal permeability, and porosity (0.573>0.072>0.067>0.052>0.024).Among them, the standard sensitivity of the formation stress is as high as 0.57.The sensitivity of factors affecting the dissolution capture capacity of CO2decreased in the order of formation stress, residual gas saturation,temperature, horizontal permeability, and porosity (0.763>0.09>0.086>0.072>0.032).Among them, the standard sensitivity of the formation stress is up to 0.76.

(iii) As one of the effective means to achieve carbon neutralization, geologic CO2storage technology still needs to be further studied.In this study, only the physical parameters of the reservoirs were studied and analyzed as the influencing factors of geologic CO2storage.Subsequently, the fracture pressure of reservoirs and the sequestration safety of cap rocks will be further investigated.The systematical study of a series of scientific problems of geologic CO2storage will allow for enriching the theoretical basis of geologic CO2storage in low-permeability reservoirs, developing China’s geologic CO2storage technologies, and achieving net zero emission of greenhouse gases such as CO2for sustainable development.

CRediT authorship contribution statement

Shi-xin Dai conceived the presented idea, Shi-xin Dai,Feng Wang, and Yan-jiao Dong prepared the manuscript.Feng Wang, Yan-jiao Dong, Zhen-han Xing, Pan Hu, and Fu Yang drew all the figures.Feng Wang and Yan-jiao Dong supervised the findings of this work.All authors discussed the results and contributed to the final manuscript.

Declaration of competing interest

The authors declare no conflicts of interest.

Acknowledgment

This study was jointly supported by the National Key R&D Program of China (2018YFB0605503), the National Natural Science Foundation of China (51804112), the National Key R&D Program of China (2018YFC0807801),the Open Foundation of Key Laboratory of Coal Exploration and Comprehensive Utilization of Ministry of Natural Resources (KF2021-5), and the Natural Science Foundation of Hunan Province of China (2018JJ3169).