Abstract
We have studied the hydrocarbon production from oil shale reservoirs filled with diverse initial saturations of fluid phases by implementing numerical simulations of various thermal in-situ upgrading processes. We use our in-house fully functional, fully implicit, and non-isothermal simulator, which describes the in-situ upgrading processes and hydrocarbon recovery by multiphase-multicomponent systems. We have conducted two sets of simulation cases—five-spot well pattern problems and Shell In-situ Conversion Process (ICP) problems. In the five-spot well pattern problems, we have analyzed the effects of initial fluid phase that fills the single-phase reservoir and thermal processes by four cases—electrical heating in the single-phase-aqueous reservoir, electrical heating in the single-phase-gaseous reservoir, hot water injection in the single-phase-aqueous reservoir, and hot CO2 injection in the single-phase-gaseous reservoir. In the ICP problems, we have analyzed the effects of initial saturations of fluid phases that fill two-phase-aqueous-and-gaseous reservoir by three cases—initial aqueous phase saturations of 0.16, 0.44, and 0.72. Through the simulation cases, system response and production behavior including temperature profile, kerogen fraction profile, evolution of effective porosity and absolute permeability, phase production, and product selectivity are analyzed. In the five-spot well pattern problems, it is found that the hot water injection in the aqueous phase reservoir shows the highest total hydrocarbon production, but also shows the highest water-oil-mass-ratio. Productions of phases and components show very different behavior in the cases of electrical heating in the aqueous phase reservoir and the gaseous phase reservoir. In the ICP problems, it is found that the speed of kerogen decomposition is almost identical in the cases, but the production behavior of phases and components is very different. It is found that more liquid organic phase has been produced in the case with the higher initial saturation of aqueous phase by the less production of gaseous phase.
Keywords
Introduction
Increasing energy demand necessitates the stable energy supply by technically and economically feasible production of energy sources. Oil shales are organic-rich sedimentary rocks, where solid organic matter—kerogen—occupies the pores. Kerogen is decomposed into multicomponent materials of fluids and solids in the environments of high temperature. Fluids include a number of species including liquid and gaseous hydrocarbons. Among the unconventional resources that have been actively researched, oil shales are believed to have tremendous reserves in the United States, approximately 75% of the world’s oil shales (Knaus et al., 2010).
To heat the kerogen in oil shales to its pyrolysis temperature, a few technologies can be applied such as mining followed by surface retorting with direct/indirect heating, in-situ upgrading process, and in-situ combustion approach. Among them, the in-situ upgrading process introduces heat into kerogen while it is located in the subsurface as in its natural state. In-situ upgrading process has several advantages over surface retorting and in-situ combustion; it can be applied to deep and thick deposits with less water requirement, emissions of pollutants, and environmental impacts (Crawford and Killen, 2010). There have been plenty of studies on the diverse methods of in-situ upgrading of oil shales including Shell In-situ Conversion Process (ICP), ExxonMobil Electrofrac, EGL Technology of American Shale Oil, LLC (AMSO), and Texas A&M Steamfrac (Crawford and Killen, 2010; Symington, 2006; Thoram and Ehlig-Economides, 2011; Vinegar, 2006).
In-situ upgrading process has a few challenges as well; far longer heating time is required than surface retorting, but excess heat should be avoided because it can reduce oil quality and energy efficiency and cause environmental impacts; shale is not a good conductor of heat, and natural or induced fracture systems are necessary to distribute the heat effectively (Crawford and Killen, 2010); and the cost of in-situ upgrading is estimated to be high by significant operating cost in the case of hot fluid injection process and significant fixed cost in the case of electrical heating, respectively (Lee, 2014). In this regard, Shell ICP is an effective method to evenly heat the oil shales using densely located multiple heaters in hexagonal patterns. Shell has been developing oil shale technologies for more than 50 years including experiments and field production test and produced over 1000 bbls of oil by applying ICP in the Mahogany Demonstration Project-South (MDP-S) (Fowler and Vinegar, 2009). The delineation of general ICP method is shown in Figure 1. Converted fluids by the conductive heating process from the hexagonally installed heaters are produced at the producing wells.
Delineation of general Shell ICP method. Correlations of porosity and permeability.
In the implementation of ICP, the initial saturation of fluid phases in oil shales, especially the amount of water, has the significant impact on the project efficiency by affecting flow of converted fluids. There have been several simulation studies on the ICP, but they did not consider the effect of differing water content in the system on the fluid flow in production process. Shen (2009) has conducted the simulation of pilot production by ICP in the Mahogany Field Experiment; heating and production processes after the 30 days of dewatering process have been simulated. Fan et al. (2010) have simulated the application of ICP using diverse heating temperatures and heater patterns; they have omitted the water component and aqueous phase in the simulations, and the effect of water in the system has been ignored.
Comparing with the ICP method of electrical heating, hot fluid injection can be considered to heat the oil shales by heat convection as well as conduction, especially in the systems with pre-existing fracture networks. Hot gas such as CO2 and hot water are favorable to be injected to gas-dominant reservoirs and water-dominant reservoirs, respectively. Lee et al. (2016a) have conducted the comparison study of diverse in-situ upgrading processes including ICP and ExxonMobil Electrofrac of electrical heating, and Texas A&M Steamfrac of hot water injection to analyze the effect of heating methods on the productivity and product selectivity in the oil shale system initially filled with 100% aqueous phase. Youtsos et al. (2013) have conducted the numerical simulation of thermal and reaction fronts during the in-situ upgrading by hot CO2 injection in the oil shale system initially filled with 100% gaseous phase.
From the previous simulation studies on the in-situ upgrading process, we have found that the effect of differing initial saturations of fluid phases, especially the water content, has not been considered thoroughly, even though it should be fully accounted for the realistic simulations. Also, the comparison of heating methods has not been done in the same geometry of reservoir system. In this regard, we conduct two sets of numerical simulations in this study.
The first set of simulations involves four cases of a five-spot well pattern problem to analyze the effects of heating methods and the initial saturations of fluid phases in the same reservoir geometry. In this set of simulation cases, we use the five-spot well pattern of simple geometry to clearly estimate and compare the system responses including the distributions of temperature, kerogen volume fraction, and effective porosity, and the production behaviors of fluid phases and components.
The second set of simulations involves three cases of ICP method as the application of actually proposed in-situ upgrading technology. We use the well patterns introduced in the pilot production of MDP-S in the Green River Formation of the United States (Fowler and Vinegar, 2009). From this set of simulation cases, we analyze the effect of water content on the implementation of ICP method in the two-phase-aqueous-and-gaseous reservoirs. In each simulation case, the system response, production behavior of fluids, product selectivity, and productivity will be analyzed.
Simulation model description
In the numerical simulation of in-situ upgrading process, heating and production procedures should be described thoroughly to predict the likelihood of feasible production of hydrocarbons from oil shales. However, the physical and chemical phenomena occurring in the thermal processes applied to oil shales are highly complex to be described, because they involve a number of subsequent/concurrent chemical reactions, continuous evolution of porosity and permeability, and changes of fluid composition and phase properties.
In this study, we use a fully functional simulator developed for the in-situ upgrading process, which describes the thermophysical and chemical phenomena in the oil shale reservoir systems (Lee et al., 2016b). The simulator is capable of describing the coupled flow of multiphase-multicomponent fluid and heat in porous and fractured media involving chemical reactions, where the porosity and permeability continuously alter by the changing amount of solid phase during the thermal process.
Chemical reactions of in-situ upgrading process.
Source: Wellington et al. (2005).
Materials in the reservoir system present as in either single-phase or multiphase-equilibrium states by the thermophysical conditions. Phase transition by fluid flow, change of reservoir pressure or temperature, chemical reactions, vaporization, condensation, and gas evolution is considered. The simulator computes the thermodynamic properties of phases by diverse equations of state and mixing rules. The properties of bulk oil shale rock are computed by considering the composition of phases in its pores and fractures.
Many oil shales are composed of dolomite and calcine, and natural fracture systems occur in such systems (Dyni, 2006). The simulator describes the naturally fractured media using the concept of Multiple Interacting Continua (MINC) (Pruess, 1985). MINC simulates the fractured rock by dividing the domain into homocentric grids for the fracture at the external position and the rock matrix at the internal position. The fracture domain mainly provides the flowing path of fluid and heat. In our simulations, the reservoir contains a well-connected natural fracture system with the dense fracture spacing of 3 m and the high fracture permeability of 1.480 × 10−13 m2 (= 150 md), and it satisfies the assumption of MINC method, which considers the naturally fractured reservoir as a continuum. We use the volumetric fraction of 0.01 of bulk oil shale rock for the fracture network.
Continuous change of porosity and permeability of rock matrix is computed at every moment by accounting kerogen decomposition, subsequent generation of fluids and solids, fluid expansion by the increase of pressure and temperature during the heating process, and pressure decrease by production. Two concepts of porosities are accounted in the simulations—media porosity and effective porosity. Media porosity involves whole phases of fluids and solid, while effective porosity involves only fluid phases. The permeability of rock is computed every moment according to the evolution of effective porosity.
In our two sets of simulation case study, we examine the effects of initial saturation of fluid phases in the pores and fractures and the diverse heating methods. In the first set of five-spot well pattern problems, we conduct four cases of simulations in the reservoir filled with single-phase fluid initially—electrical heating in the single-phase-aqueous reservoir, electrical heating in the single-phase-gaseous reservoir, hot water injection in the single-phase-aqueous reservoir, and hot CO2 injection in the single-phase-gaseous reservoir. In the second set of ICP problems, we conduct three cases of simulations in the reservoir filled with two-phase-aqueous-and-gaseous fluids with different water contents to each other—SiA = 0.16, 0.44, and 0.72 (SiG = 0.84, 0.56, and 0.28).
As well as the different initial saturation of fluid phases and the diverse heating methods, we use the different irreducible saturations of phases, different correlation models of porosity and permeability, and different initial effective porosity between the two sets of problems. Their conditions of the five-spot well pattern problems have been obtained from the validated model of previous study (Lee et al., 2016b). First, in the set of five-spot well pattern problems, we use the irreducible saturations of SiA = 0.30, SiO = 0.20, and SiG = 0.05 at the matrix domain, and SiA = 0.03, SiO = 0.02, and SiG = 0.005 at the fracture domain, respectively; in the set of ICP problems, we use the smaller irreducible saturations of SiA = 0.16, SiO = 0.12, and SiG = 0.04 at the matrix domain, and SiA = 0.016, SiO = 0.012, and SiG = 0.004 at the fracture domain, respectively, to counterbalance the flow of two-phase fluid presenting initially. Second, we use different correlation models of porosity and permeability in the two sets of problems as shown in Figure 2, even both models are based on the core data of Green River Formation (Kibodeaux, 2014). In the set of five-spot well pattern problems, we use the correlation of Nimblett and Ruppel (2003), which was used in the validated model of previous study (Lee et al., 2016b); in the set of ICP problems, we use the modified Kozeny–Carman equation for the better representation of core data (Krauss and Mays, 2014). Third, in the set of five-spot well pattern problems, we use the initial effective porosity of 0.04 and corresponding absolute permeability of 1.228 × 10−16 m2 (= 0.124 md) at the rock matrix; in the set of ICP problems, we use the initial effective porosity of 0.015 and 3.092 × 10−19 m2 (=3.133 × 10−4 md) at the rock matrix, for more realistic description of low-permeable raw shales.
Simulation of five-spot well pattern problems
Simulation model and initial reservoir properties
A quarter of a five-spot well pattern system in a well-connected naturally fractured reservoir is simulated by using a 2D model with 1 m thickness (Figure 3). The heater (or the injection well) is located at the corner of the model, and the production well is located at the other corner on a diagonal line. The reservoir properties and the initial conditions have been obtained from the Green River Oil Shale, as provided in Table 2 (Lee et al., 2016b). Two initial conditions of fluid phases have been considered—single phase aqueous and single phase gaseous. The initial aqueous phase and gaseous phase are occupied by 100% of water component and 100% of CO2, respectively. For the computation of relative permeability and capillary pressure in the matrix domain, Van Genuchten model and Parker’s model have been used for the two-phase flow and the three-phase flow, respectively.
Reservoir geometry and simulation model of five-spot well pattern cases. Reservoir properties and initial conditions in the simulation cases of five-spot well pattern problems.
Four cases of thermal and production processes for 1000 days are simulated—electrical heating in the aqueous phase reservoir, electrical heating in the gaseous phase reservoir, hot water injection in the aqueous phase reservoir, and hot CO2 injection in the gaseous phase reservoir. In the cases of electrical heating by an electric heater, the heater grid block is continuously heated during the whole process of heating and production, and the constant heater temperature of 618.15 K is used. In the cases of hot fluid injection, the heated fluids are injected by the constant rate of 5 × 10−4 kg/s in the simulation model (= 2 × 10−3 kg/s in one five-spot well pattern system), and by the variable enthalpy that makes the injection fluid temperature = 618.15 K. The heater/injection well temperature for each case during the process is provided in Figure 4. It shows that the heater temperature in the electrical heating cases was kept constant from t = 0 days to t = 1000 days by the constant temperature at the heater well grid, while the injection well temperature in the cases of hot fluid injection was gradually increasing from the initial reservoir temperature to the objective temperature. The injection well temperature approached the objective temperature of 618.15 K at t = 400 days in the injection cases.
Temperature of the heater/injection well of five-spot well pattern cases.
From t = 0 days to 50 days, only thermal process is applied in the cases of electrical heating, while the same rate of fluid to the injection rate is produced in the cases of hot fluid injection. Then the converted fluids are produced from t = 50 days to 1000 days in each case using the gradually decreasing wellbore flowing pressure.
Simulation results and discussion
System behavior
The distributions of temperature on the diagonal line AB in Figure 3 during the heating and production processes are provided in Figure 5. The case of hot water injection shows the fastest temperature increase among the cases. This has been caused by that the flowing hot water has allowed the heat to be transferred more vigorously by convection, than the cases of electrical heating. Comparing with the case of hot CO2 injection, hot water has the higher heat capacity than hot CO2 of same temperature; thus the hot water injection case shows even faster temperature increase in the reservoir than the hot CO2 injection case.
Reservoir profiles of five-spot well pattern cases—distributions of temperature on the line AB.
Figure 6 shows the distributions of kerogen volume fraction in solid phase. By the fast increase of temperature, the case of hot water injection shows the active kerogen decomposition in the wide area from the early time of 250 days. The distributions of effective porosity of the matrix domain are provided in Figure 7. The case of hot water injection has the fastest increase of effective porosity in the wide area among the cases, by the fast increase of temperature and the active decomposition of kerogen. At the end of the process at t = 1000 days, the effective porosity in each case have approached 0.2.
Reservoir profiles of five-spot well pattern cases—distributions of kerogen volume fraction in solid phase on the line AB. Reservoir profiles of five-spot well pattern cases—distributions of effective porosity on the line AB.
The increase of effective porosity in each case during the processes of heating and production has been resulted by the kerogen decomposition and fluid expansion, but the final effective porosity is still smaller than the initial media porosity of 0.316 because of the generation of solid products. In the cases of electrical heating and hot CO2 injection in the gaseous phase reservoir, it is found that the effective porosity of the matrix domain at distance = 0.75 m at t = 750 days is higher than the final effective porosity. This has been caused by that the increase of effective porosity by the gaseous phase expansion in the matrix domain at t = 750 days is higher than that at t = 1000 days in those cases, by the higher reservoir pressure at t = 750 days.
Kerogen mass in place and productions of phases
Figure 8 provides the kerogen mass in place and the production rates of phases in one five-spot well pattern system in the simulation cases. In Figure 8(a), it is found that the kerogen mass in place rapidly decreases from t = 400 days, 370 days, 200 days, and 370 days, until it reaches to zero at t = 750 days, 710 days, 470 days, and 710 days, in the cases of electrical heating in the aqueous phase reservoir, electrical heating in the gaseous phase reservoir, hot water injection in the aqueous phase reservoir, and hot CO2 injection in the gaseous phase reservoir, respectively.
Results of five-spot well pattern cases—kerogen mass in place and production rates of phases.
Figure 8(b) and (c) show that aqueous phase and liquid organic phase are not produced in the cases of electrical heating and hot CO2 injection in the gaseous phase reservoir during the process, and they have remained at the matrix domain after being generated from the kerogen decomposition due to their far lower mobility than that of gaseous phase in such system where gaseous phase is dominant.
In Figure 8(b), it is found that the peak of the production rate of aqueous phase is 8.103 × 10−5 kg/s (= 0.0587 STB/day) at t = 410 days in the case of electrical heating in the aqueous phase reservoir. In the case of hot water injection in the aqueous phase reservoir, the aqueous phase production rate remains almost same as the injection rate after t = 600 days.
In Figure 8(c), it is found that the peaks of the production rate of liquid organic phase are 7.267 × 10−5 kg/s (= 0.0702 STB/day) and 3.229 × 10−4 kg/s (= 0.312 STB/day) in the cases of electrical heating in the aqueous phase reservoir and hot water injection in the aqueous phase reservoir, respectively.
Cumulative productions of phases in the five-spot well pattern problems.
GOR: gaseous phase–liquid organic phase; WOR: aqueous phase–liquid organic phase.
Productions of components
The production rates of fluid components are provided in Figure 9, and the cumulative productions of components are listed in Table 4. From Figure 9(a), (b), and (f), it is found that the heavy oil, light oil, and water components are produced in the cases of electrical heating and hot CO2 injection in the gaseous phase reservoir, even though liquid organic phase and aqueous phase have not been produced in those cases as shown in Figure 8. In these cases, those components are produced in the form of vapor in gaseous phase.
Results of five-spot well pattern cases—production rates of components. Cumulative productions of components in the five-spot well pattern problems.
Comparing Figure 9(d) and (f) with the Figure 8(b) and (d), it is found that the aqueous phase and gaseous phase are composed of almost 100% water and 100% CO2 in the cases of hot water injection in the aqueous phase reservoir and hot CO2 injection in the gaseous phase reservoir, respectively. From Table 4, it is found that the case of hot water injection in the aqueous phase reservoir shows the highest cumulative production of every fluid component except for CO2, which shows the highest value of cumulative production in the case of hot CO2 injection.
Simulation of ICP
Simulation model and initial reservoir properties
The simulation examples include the three cases having the diverse water contents in the reservoir systems. The three reservoirs are naturally fractured and initially filled with two phases of fluid—aqueous and gaseous. The initial volume fraction of the aqueous phase in the void spaces of pores and fractures are 0.16, 0.44, and 0.72 in case 1, case 2, and case 3, respectively. The initial aqueous phase and gaseous phase are occupied by 100% of water component and 100% of hydrocarbon gas, respectively. The effects of initial saturation of aqueous phase on the productivity and product selectivity are evaluated through the simulation cases.
In the simulation cases, electrical heating by multiple electric heaters of ICP has been applied, with the specific well patterns proposed by MDP-S pilot production in Green River Formation (Fowler and Vinegar, 2009). Multiple electric heaters and producers are located in hexagonal patterns in the reservoir (Figure 10(a) and (b)). A quarter of one hexagon has been simulated using a 2D model as shown in Figure 10(c), which has 34.4 m thickness as in the pilot production implementation. The closed boundary condition for the flows of heat and fluid has been applied at the boundaries. In the simulations, heating and production are continued until t = 300 days. The electric heaters are heated up to T = 618.15 K (=345℃). The converted fluids by the heating process are produced using the gradually decreasing wellbore flowing pressure.
Reservoir configuration and simulation model of ICP cases.
Reservoir properties and initial conditions in the simulation cases of ICP.
ICP: in-situ conversion process.
Simulation results and discussion
System behavior
The distribution profiles of temperature and kerogen mass in place during the process in case 1 are shown in Figure 11. The formation temperature increases gradually from region around the heaters and approaches evenly high temperature at t = 230 days. The kerogen in the system gradually decomposes from the region around the heaters, and whole decomposes at t = 230 days. Cases 2 and 3 have showed the similar distribution profiles to case 1.
Reservoir profiles of ICP case 1.
The system responses of each case are shown in Figure 12. Oil shale formation of 34.44 m interval has been considered. Figure 12(a) shows that kerogen in the system actively decomposes from t = 100 days in each case. The curves of kerogen mass in place show insignificant difference among the cases. This is because, the initial volume fraction of void space of pores and fractures was 0.025, and the difference of thermal conductivity of bulk shale is less than 2% between the cases considering the following composite thermal conductivity of bulk oil shale rock.
System behavior of ICP cases—kerogen mass in place in the system, and effective porosity and absolute permeability of rock matrix at the monitoring point (x, y) = (2.56 m, 2.56 m).
Figure 12(b) and (c) indicate the evolution of effective porosity and corresponding matrix permeability by the kerogen decomposition at the monitoring point (x, y) = (2.56 m, 2.56 m). The values of final effective porosity are 0.1748, 0.1742, and 0.1739 in the cases, respectively. They are smaller than the initial media porosity of 0.291 because of the solid products remaining in the pores. The effective porosity and matrix permeability slightly decrease from t = 220 days because the system pressure decreases due to the production, while temperature has stopped increasing.
Production of phases
Production is continued from t = 15 days to t = 300 days in each case. The production rates and cumulative productions of fluid phases are provided in Figure 13. In each case, the production rates of liquid organic phase and gaseous phase increase as kerogen actively decomposes, and gradually decreases after their peaks as the supply of hydrocarbons from kerogen decomposition decreases. The production rates of aqueous phase and liquid organic phase approach zero around t = 150 days and t = 250 days, respectively. After t = 250 days, only a small amount of gaseous phase is produced.
Results of ICP cases—production rates and cumulative productions of fluid phases.
Production rate of liquid organic phase shows the highest peaks in each case, and gaseous phase and aqueous phase follow it. The peaks of the production rate of liquid organic phase are 1.100 × 10−2 kg/s, 1.216 × 10−2 kg/s, and 1.234 × 10−2 kg/s; the peaks of the production rates of gaseous phase are 1.021 × 10−2 kg/s, 9.252 × 10−3 kg/s, and 9.204 × 10−3 kg/s; the peaks of the production rates of aqueous phase are 2.566 × 10−3 kg/s, 2.484 × 10−3 kg/s, and 4.863 × 10−3 kg/s in cases 1, 2, and 3, respectively. It is noticeable that case 2 with SiA = 0.44 has the lower peak of the production rate of aqueous phase than case 1 with SiA = 0.16, by the higher production rate of liquid organic phase.
Cumulative productions of phases in the ICP cases.
ICP: in-situ conversion process.
The saturation curves of fluid phases at the monitoring point (x, y) = (2.56 m, 2.56 m) are provided in Figure 14. In the fracture domain, gaseous phase replaces the liquid organic phase and aqueous phase as production proceeds, while above 60% of liquid organic phase still has remained in the matrix domain at the end of the process in each case. The curves of relative permeability from Parker’s model and the mobility of liquid organic phase and gaseous phase in the system with SA = 10% are provided in Figure 15. Despite the higher relative permeability of liquid organic phase above So = 0.61, the far higher mobility of gaseous phase has caused the flow of only gaseous phase from the matrix domain to the fracture domain and the trapping of liquid organic phase in the matrix domain at the end of the process.
Results of ICP cases—saturations of fluid phases in the fracture domain and matrix domain at the monitoring point (x, y) = (2.56 m, 2.56 m). Relative permeability and mobility curves of liquid organic phase and gaseous phase in the matrix domain, where SA = 0.1.
Production of components
The production rates and cumulative productions of fluid components during the process are provided in Figure 16. Hydrocarbon gas, light oil, and water components have taken the highest portion of the produced fluids in the case 1, 2, and 3, respectively. The cumulative productions of the components are summarized in Table 7. Cumulative production of every component except for hydrocarbon gas increases with increasing initial saturation of aqueous phase. More carbon dioxide and hydrogen have been produced in the case with higher initial saturation of aqueous phase; they have been produced as dissolved in aqueous phase and liquid organic phase, as well as in the form of free gas. In every case, the higher cumulative production of light oil than heavy oil has been caused by the cracking reactions of heavy oil into lighter hydrocarbons, as well as by the higher mobility of light oil.
Results of ICP cases—production rates and cumulative productions of components. Cumulative productions of components in the ICP cases. ICP: in-situ conversion process.
Conclusion
In this study, we have investigated the effects of initial saturations of fluid phases in the reservoir and diverse thermal methods on production behavior by conducting two sets of simulations—five-spot well pattern problems and ICP problems. From the first set of five-spot well pattern problems, we have the following conclusions.
The increase of system temperature and the decomposition of kerogen in the system have been fastest in the case of hot water injection in the aqueous phase reservoir, by the vigorous heat transfer by conduction and convection. The final effective porosity and absolute permeability of rock matrix in each case has approached 0.2 and 3.111 md, respectively. The case of hot water injection in the aqueous phase reservoir has showed the largest cumulative productions of hydrocarbons (heavy oil, light oil, and hydrocarbon gas) among the cases, even also has showed the largest WOR. Enormous energy loss and huge operating cost are expected in the case of hot water injection. The envisioned work will include the analysis of energy efficiency and economic feasibility of each case. Even the cases have different initial saturations of fluid phases in the system, the difference of speed of kerogen decomposition is insignificant by the low initial effective porosity of the system. The final effective porosity and absolute permeability of rock matrix in the cases have approached 0.17 and 0.642 md, respectively. In every case, gaseous phase has showed the highest mass of cumulative production among the phases. More liquid organic phase has been produced with the higher SiA, by the less production of gaseous phase. At the end of the heating and production process, only gaseous phase has flowed through the fracture domain, while above 60% of liquid organic phase has been trapped in the matrix domain due to its far lower mobility than the gaseous phase.
From the second set of ICP problems, we have the following conclusions.
Footnotes
Acknowledgments
The authors would like to appreciate the Crisman Institute for Petroleum Research of Texas A&M University for supporting this study.
Declaration of conflicting interests
The author(s) declared no potential conflict of interest with respect to the research, authorship, and/or publication of this article.
Funding
The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This study was financially supported by the Crisman Institute for Petroleum Research of Texas A&M University (Crisman project number: 1.7.05b).