Characteristics and mechanisms of supercritical CO2 flooding under different factors in low-permeability reservoirs

Zheng Chen ,Yu-Liang Su ,*,Lei Li ,Fan-Kun Meng ,Xiao-Mei Zhou

a School of Petroleum Engineering,China University of Petroleum(East China),Qingdao,266580,Shandong,PR China

b Key Laboratory of Unconventional Oil ＆ Gas Development(China University of Petroleum(East China)),Ministry of Education,Qingdao,266580,Shandong,PR China

c School of Petroleum Engineering,Yangtze University,Wuhan,430100,Hubei,PR China

Keywords:

ABSTRACT

1.Introduction

As the total global demand for oil continues to increase,the output of proven conventional oil reservoirs cannot meet the needs of consumption growth,and the contradiction between supply and demand has become increasingly prominent.This will inevitably require people to continue to explore unconventional oil reservoirs(Chen et al.,2021;Dong et al.,2020;Pei et al.,2020).Reasonably exploiting this type of oil reservoir plays an important role in the development of the petroleum industry(Zhao et al.,2020a).Lowpermeability reservoirs have tight lithology and small particles and exhibit low porosity and low permeability.The most commonly used development method is water injection.However,the high starting pressure and high remaining oil content make it difficult to effectively implement water flooding in lowpermeability reservoirs(Fang et al.,2019b;Fan et al.,2015;Shi et al.,2017;Wang et al.,2014;Hu et al.,2018).Therefore,gas injection becomes an effective method to enhance oil recovery after water flooding.The density of supercritical carbon dioxide is close to that of liquid,and it has the characteristics of good injection performance and low viscosity(Wang et al.,2018;Li et al.,2006;Xiong et al.,2015).In the context of global warming,CO2flooding is recognized as one of the most potent methods for enhancing oil recovery(Fang et al.,2019a;Wei et al.,2021;Zhao et al.,2015)and mitigating the greenhouse environment problem at the same time.In recent years,it has been widely used in the development of major unconventional oil fields.Field tests of CO2flooding have shown that it is easier to inject CO2in low-permeability reservoirs than water.CO2flooding can increase oil recovery by 10% over water flooding.The reservoir also provides a good area for carbon dioxide storage(Cheng et al.,2017).

Core flooding experiments and numerical simulation methods are effective approaches to investigate the characteristics of CO2flooding(Sander et al.,2017;Sun et al.,2016).There are few studies of supercritical CO2using numerical simulation methods,which mainly focus on using commercial software and self-built mathematical models to simulate the flow process and variation characteristics of various parameters of supercritical CO2in the displacement process(Cui et al.,2011;Zhang et al.,2020;Zhu et al.,2016).At present,many scholars still take experimental means to understand the displacement characteristics of supercritical CO2.It is believed that the permeability of a low permeability reservoir can be improved by injecting supercritical CO2first by adjusting injection and production parameters and injection methods.Meanwhile,the influence of gravity effect in immiscible flooding cannot be ignored(Bikkina et al.,2016;Han et al.,2016;Wang et al.,2011).However,the oil samples used in their studies are single hydrocarbon compounds,and the rock samples used are generally short cores,and only one factor is considered to affect the supercritical CO2displacement behavior.Given the above problems,Li et al.(2018)studied the oil-increasing laws of various supercritical CO2displacement methods by using long cores and live oil.They believe that the miscible flooding has the minimum injection pressure difference and the maximum recovery factor,and the supercritical CO2alternating water-miscible flooding has a lower injection pressure and smaller gas and oil production.At the same time,they also believe that the increase in gas-oil ratio does not represent a breakthrough of supercritical CO2,and gas injection at high water cut stage has a lag in water content change(Li et al.,2018).In subsequent studies,the influences of pressure,supercritical CO2content,soaking time,and permeability on supercritical CO2displacement characteristics were studied successively.However,the rock samples and oil samples used were short core and degassed oil,which differed greatly from the actual reservoir conditions(Chen et al.,2018;Li et al.,2019b;Xu et al.,2020).

In view of the strong heterogeneity of low-permeability reservoirs,Park et al.(2017)studied the influence of core sedimentary heterogeneity and displacement pressure difference on CO2flooding(Park et al.,2017).They revealed that the capillary pressure difference caused by core heterogeneity is one of the main reasons for the CO2displacement characteristics.At the same time,many scholars also use the experimental method of assembling cores of different permeability into parallel or vertical combination samples to study the CO2flooding effect under different injectionproduction parameters(Al-Bayati et al.,2019;Wang et al.,2021;Wei et al.,2014).

The above studies discussed the characteristics of supercritical CO2flooding at the macro scale of the core.Some scholars evaluated water flooding and supercritical CO2flooding in core flooding experiments by combining CT scanning experiment and nuclear magnetic resonance experiment with macro core experiment.Thus,the saturation and fluid distribution changes of porous cores under typical reservoir temperature and pressure conditions are clarified(Qian et al.,2020;Mahabadi et al.,2020;Zhao et al.,2020b).At the same time,to study the details of pore-scale flow patterns in the porous media that cannot be resolved,some microsimulation methods are used to study and understand the dynamic displacement process of carbon dioxide injection in porous media,for example the Boltzmann method(Liu et al.,2014,2020;Yamabe et al.,2015).

All previous studies have shown that gas channeling and low oil recovery under supercritical CO2flooding of low permeability cores are inevitable(Shen et al.,2021;Wang et al.,2019;Zhang et al.,2012).It is necessary to study the characteristics and mechanisms of supercritical CO2flooding under different factors.However,as shown in Table 1,for most of the published studies,the experimental design is focused on one kind of influencing factor.The rock samples used are short cores,and the gas breakout and interflow occurred earlier during gas injection flooding,so it is difficult to accurately reflect the seepage parameters during the displacement process.The oil sample is a mixture of degassed oil or other decanes,and the influence of the original gas-oil ratio of formation oil and the change of degassed oil composition is not considered.In general,there are few reports on the study of supercritical CO2displacement characteristics by using long core and live oil.Therefore,it is necessary to carry out supercritical CO2flooding long core experiments to comprehensively analyze the displacement characteristics and influence mechanism of supercritical CO2on live oil under different injection and production parameters.

Table1 Summary of experimental research on CO2 flooding.

In this study,the live oil configuration and inspection method are independently designed to ensure the matching of the gas-oil ratio and the stability of live oil.At the same time,the supercritical CO2flooding experiment is carried out using natural outcrop long cores in low-permeability reservoirs.This paper studies and discusses various parameters that affect CO2displacement efficiency,gas channeling characteristics,and inject capacities,such as reservoir physical properties,formation conditions,and injectionproduction system.This paper uses long core samples and live oil,which is closer to reservoir conditions than previously reported.Based on these discussions,we try to clarify the influence mechanism of different factors on the characteristics of CO2flooding.

2.Experimental

2.1.Experimental materials

The core samples are outcrop cores,with their basic parameters shown in Table 2.Permeability was determined by a standard measuring method of 99.99% pure nitrogen,while porosity was determined by the liquid weighing method.The purity of CO2gas used in the experiment was 99.99%.The live oil sample used for core saturation was based on the composition of the gas sample provided by the oilfield,original gas-oil ratio,and crude oil viscosity.The oil properties are shown in Table 3.Distilled water was used in the experiment,while the influence of salinity on the core was ignored.

Table2 The basic parameters of core samples.

2.2.Procedures of testing

The schematic of the extended core displacement experiment device is shown in Fig.1.A Vindum pump provides driving power and a confining pressure pump provides confining pressure and back pressure.Confining pressure,back pressure,and inlet pressure were measured using pressure gauges.The volumetric cylinder,shown in Fig.1,was used to measure the liquid volume.The gas volume was measured by discharging sodium carbonate solution by gas.

The block studied in this research is block G of the Shengli Oilfield in China.The original reservoir pressure and temperature are 42.6 MPa and 126?C.The development in recent years has caused a decrease in reservoir pressure to the current value of 19.0 MPa.The experimental temperature was set at 50?C,higher than the critical temperature of carbon dioxide.The experimental back pressure was set at the current reservoir pressure of 20 MPa,ensuring that the injected CO2was in a supercritical state.The experimental steps mainly included core processing and displacement as follows.

(1)According to the crude oil characteristic data given by the oil field,live oil was configured by a sample mixing device.Based on the original gas-oil ratio data provided by the oilfield,the required gas and dead oil were calculated.Then,the dead oil,gas,and a small amount of kerosene were mixed in the sample mixing equipment to fully dissolve the gas in crude oil.

(2)A portion of live oil was extracted and measured separately in liquid and gas volumes.The relative error between actual and theoretical gas-oil ratio was calculated and controlled within 10% to check whether the gas in live oil firmly dissolves.The relative error between the actual gas-oil ratio and the live gas-oil ratio was 6.38%.

(3)A part of the live oil was then introduced into the viscometer to measure its viscosity before calculating its relative error with respect to the viscosity of the actual formation crude oil to control it within 10%(Table 3).Suppose the configured live oil properties are different from the actual,repeat steps(1)-(3)to reconfigure the oil.The viscosity of the live oil was 2.55 mPa?s and the relative error was 3.66%.

Table3 Live oil properties.

(4)The live oil that meets the standard was pressure-maintained and exported to the accumulator.A pressure tank was used for real-time pressure monitoring of the export process.

(5)The outcrop core was cleaned and dried before measuring the permeability of the core using the gas measurement method.The porosity of the cores was measured by saturating the cores with distilled water.The results are shown in Table 2.

(6)The prepared live oil was injected into the water-saturated cores with a backpressure of 13 MPa.The irreducible water saturation was calculated by flooding with a displacement pressure of 20 MPa until no water flowed out at the end of the core.

(7)Supercritical CO2flooding was conducted on oil-bearing cores under different factors.The experimental design parameters are shown in Table 4.The injection methods were divided into three types:continuous gas injection(CGI),synchronous injection(SI),and asynchronous injection-production(AIP).Section 3.2 describes the schemes of synchronous and AIP.

Table4 The operation parameters for different cases.

Fig.1.Experimental setup for supercritical CO2 displacement:(a)experimental device;(b)schematic diagram.

3.Results and discussion

The experimental methods in Section 2 analyze the oil displacement characteristics under different injection rates,injection-production modes,permeability,pressure levels,and water saturation.This section discusses the influence mechanisms of different factors of supercritical CO2flooding in low-permeability reservoirs by comparing the recovery factor,gas-oil ratio,and injection capacity under various factors.

3.1.Displacement characteristics at different injection rates

The recovery of each stage at different injection rates is shown in Fig.2.The oil recovery rate rises slowly before the gas breakthrough at three different injection rates.In this stage,CO2dissolves in crude oil,reducing the viscosity of crude oil and improving the fluidity ratio.When supercritical CO2breaks through at the end of the core,the dynamic balance among the original CO2,crude oil,and water is broken due to the decrease in pressure,forming the three-phase flow of oil,gas,and water.At this point,gas begins to appear at the end of the core.This stage depicts maximum oil production due to the carrying effect of supercritical CO2on crude oil during the breakthrough.Compared with low injection rate,the stable production stage of high injection rate becomes more extended,which leads to an increase in final recovery(Ajoma et al.,2021;He et al.,2019;Li et al.,2019a).However,the recovery efficiency before the breakthrough at the gas injection rate of 0.1 mL/min is more significant than that at 0.5 mL/min,proving the optimal value.Fig.3 shows that the gas-oil ratio rises sharply with the increase in injection volume,while the gas forms a controlled channel.The utilization degree of crude oil by the subsequent injection of CO2decreases sharply,and the recovery degree no longer increases.When the injection rate is 0.05 mL/min,the final oil recovery is 80%.Likewise,with the injection rate of 0.1 mL/min,the final oil recovery is 88.42%,and with the injection rate of 0.5 mL/min,the final oil recovery is 92.86%.The higher the core pressure level and increase in injection rates,the more fully miscible supercritical CO2and crude oil,resulting in more significant oil recovery.The increase in the injection rate can prolong gas channeling occurrence(Matkivskyi et al.,2021),as shown in Fig.3.The first occurrence of gas channeling is when the gas injection rate is 0.05 mL/min.When the gas injection rate is 0.05 mL/min,the gas channeling occurs before the gas injection rate is 0.1 mL/min,proving that gas channeling can be delayed by increasing the gas injection rate.Still,it cannot increase indefinitely because the increased range has an optimal value.

The injection index was introduced to evaluate the injection capacity to facilitate the analysis of the injection capacity during the supercritical CO2flooding process.The meaning of injection index is the volume of gas injected per unit time under unit pressure difference,which can be expressed as follows.

where I is the gas injection index,viis the injection rate,P2is the inlet pressure,and P1is the back pressure.

Fig.2.Relationship between oil recovery and injection volume(a)and gas injection rate(b).

Fig.3.Curves of gas-oil ratio versus injection volume(a),and injection volumes occurring gas breakthrough and gas channeling(b)at different gas injection rates.

The variation tendency of the injection index at different injection rates is shown in Fig.4.As the injection rate increases,the inlet pressure increases to a certain extent,but this is not a sign of poor injection capability.It can be understood that the higher the injection rate before the occurrence of gas channeling,the larger the injection index,proving a better injection capability.In the initial injection stage,the injection index was higher.This is because the initially injected supercritical CO2is easily dissolved in the crude oil,resulting in better injection capacity.After gas channeling occurs,the injection index increases sharply due to the formation of gas channels in the core.At this time,the supercritical CO2can be easily injected.

3.2.Displacement characteristics of different injection-production modes

In order to study the displacement characteristics of different injection-production modes,three injection-production modes were designed.In Mode 1,supercritical CO2was continuously injected into the cores at an injection rate of 0.01 mL/min.In Mode 2,the outlet was closed for initial injection and opened after 0.5 PV of supercritical CO2was injected.The injection rate was kept at 0.01 mL/min in Mode 2.In Mode 3,the outlet was initially opened,0.25 PV of supercritical CO2was injected at an injection rate of 0.01 mL/min,and then the outlet was closed.After injecting another 0.25 PV,the outlet was opened.

Fig.4.Comparison of injection index at different injection volumes and gas injection rates.

The recovery rates at different stages under different injection and production methods are shown in Fig.5.The oil displacement efficiency of Mode 1 increases rapidly.However,the increased rate of the flooding efficiency of Mode 2 and Mode 3 after the soaking process is faster than that of Mode 1.After the injection and production coupling adjustment,the recovery rate before gas breakthrough greatly increased,and the most considerable improvement was in Mode 2.The highest final oil recovery of 95.56% was achieved in Mode 2,followed by 90.82% in Mode 3.Mode 1 has the lowest final oil recovery of 74%.The reason is that the increase in supercritical CO2and soaking time can make the crude oil fully contact with supercritical CO2.The decrease in crude oil viscosity and the enhancement of dissolved gas flooding increase final recovery.Gas injection in advance also has an energy enhancement effect on the reservoir.The more gas injection,the more pronounced the energy enhancement effect,the larger the gas swept area,and more increased the original recovery factor(Fan et al.,2020;Gao et al.,2021;Zhao et al.,2015).

Similarly,in Modes 1,2,and 3,the amount of supercritical CO2consumed is 2.01,1.93,and 1.58 PV,respectively.This shows that adjusting the injection-production system for continuous gas flooding reduces the volume of the injected supercritical CO2,thereby maximizing the use of injected supercritical CO2.The variation curve of gas-oil ratio under different injectionproduction modes is shown in Fig.6.Mode 1 has the shortest gas channeling time,followed by Mode 3,then Mode 2.The long soaking time of Mode 2 leads to more supercritical CO2being dissolved in crude oil in the soaking stage,which can effectively delay the occurrence of gas channeling.

The variation tendency of the injection index under different injection-production modes is shown in Fig.7.The injection index of each mode in the early stage tends to be the same.Because gas channeling occurs earlier in Mode 1,the injection index rises earlier.The variation tendency of the injection capacity of each mode is consistent with changes in the gas-oil ratio.It is concluded that the injection capacity of supercritical CO2is higher in the synchronous injection-production mode.The formation pressure increases and the supercritical CO2injection capacity deteriorates with the increase in soaking time.

3.3.Displacement characteristics at different permeability

Fig.5.Relationship between oil recovery and injection volume(a)and injection-production modes(b).

Fig.6.Curves of gas-oil ratio versus injection volume(a),and gas breakthrough time and gas channeling time(b)in different injection-production modes.

Fig.7.Comparison of injection index at different injection volumes and injection-production modes.

The influence of core permeability on oil recovery at each stage is shown in Fig.8.It can be seen from Fig.8 that the final oil recovery of the rock samples with permeability of 5.1,2.5,and 0.63 mD are 80%,75%,and 56.12%,respectively.It can be seen that with the increase in permeability,the recovery factor before gas breakthrough increases,leading to the recovery factor increasing and reaching a specific value.The core with low permeability has a large proportion of tiny pores.In supercritical CO2displacement,oil in small pores is challenging to be displaced,resulting in a low recovery rate.Meanwhile,the mixing rate of supercritical CO2and crude oil increases with permeability.The higher the permeability is,the more likely it is that the supercritical CO2fingering will occur(Niu et al.,2020;Xiao et al.,2017).As shown in Fig.9,When the core permeability increases further,the gas channeling occurs earlier.The low-permeability core has poor percolation capacity and supercritical CO2has a low passing rate,so the gas-oil ratio rises slowly.The gas diffusion rate also changes with the change of permeability.When the permeability is below a specific value,the diffusion rate is minimal with little difference.When the permeability increases to a particular value,the gas diffusion rate rises rapidly.Therefore,both gas breakthrough time and gas channeling time shorten with the increase in core permeability,and the latter shortens more sharply with the increase in permeability(Amarasinghe et al.,2020).

The variation tendency of the injection index of the cores with different permeability is shown in Fig.10.The core with a permeability of 5.8 mD has a higher injection index and the best injection capacity.The core with a permeability of 0.63 mD tends to be stable.The rapid growth phase of the injection index occurred earlier due to the early occurrence of gas channeling in high-permeability cores.In cores with high permeability,the injection capacity is more significant due to quicker supercritical CO2seepage rate.

3.4.Displacement characteristics at different back pressures

The recovery rates at different stages under different back pressures are shown in Fig.11.By comparing the recovery rates at the early stage of gas injection under different back pressures,it is found that at 10 MPa,the recovery factor before gas channeling is the lowest.With the gas injection,the lower the back pressure,the earlier the gas cut,the earlier the gas channeling,and the lower the final recovery factor.The oil recovery is 81.9%,87%,and 91%,respectively,at a back pressure of 10,20,and 30 MPa.Due to the slow dissolution rate of CO2and crude oil,the degree of mutual dissolution of oil and gas is low,promoting gas channeling and lower final recovery.At high injection pressure,most of the injected supercritical CO2can be diffused to the original stagnant area to replace crude oil to increase production.At the same time,increasing pressure can enhance the miscibility of supercritical CO2and crude oil and change the wettability of rock to improve oil recovery,which has been mentioned in previous reports(Ding et al.,2017;Du et al.,2018;Lu et al.,2021;Ren et al.,2011).When the pressure reaches above the miscible pressure,the influence of the pressure on the gas channeling time is negligible.Therefore,a reasonable injection pressure should be selected to prolong gas channel time,improve swept volume and oil recovery,and ensure the maximum economic benefit during field test implementation(see Fig.12).

Fig.8.Relationship between oil recovery and injection volume(a)and permeability(b).

Fig.9.Curves of gas-oil ratio versus injection volume(a),and gas breakthrough time and gas channeling time(b)in different permeability cores.

Fig.10.Comparison of injection index at different injection volumes and permeability.

The variation tendency of the injection index under different back pressures is shown in Fig.13.The injection index at a back pressure of 10 MPa is larger in the initial gas injection stage.With a back pressure increase during CO2injection,the injection index can be reduced to a certain value.The injection capacity is generally better when the back pressure is 10 MPa.The smaller the back pressure,the earlier the gas channeling occurs,resulting in better injection capabilities.At the same time,compared with 10 MPa,the dissolved amount of supercritical CO2increases when the back pressure is 30 MPa,which delays the occurrence of gas channeling,resulting in a lower injection index and poor injection capacity.

3.5.Displacement characteristics at different water saturation

In this part,the dry core is saturated with water and oil displacement is followed to establish bound water saturation.Secondly,water flooding is conducted on the oil-bearing cores,and different water saturation is established according to the amount of oil produced.This method simulates gas displacement characteristics at different water flooding stages considering bound water saturation.

Oil recovery at different stages with different water saturation is shown in Fig.14.The lower the water saturation,the more oil will be displaced in the early stage,resulting in greater displacement efficiency.Gas displacement mainly depends on the carrying effect after gas breakthrough with the increase in water saturation.At a water saturation of 42%,60%,and 80%,the final oil recovery was 81.05%,73.44%,and 62.5%,respectively.At high water saturation,oil distribution is relatively dispersed,mostly in tiny pores,resulting in a poorer displacement effect of supercritical CO2.The presence of water reduces the contact area between oil and CO2to a certain extent,reducing the diffusion of CO2-oil in water-wet pores and making more oil become immovable(Fern?et al.,2015).This shows that when switching to gas injection,the displacement efficiency of supercritical CO2is different(Fern?et al.,2015;Torabi and Asghari,2010).Considering the economic cost of supercritical CO2,choosing the right time for gas injection is very important.Furthermore,when the water saturation increases to a certain value,the water content in the cores is high.At the same time,as shown in Fig.15,the solubility of supercritical CO2in water is small,therefore the gas drive front is easy to break through,causing a reduction in gas breakthrough and gas channeling time.Previous studies in the literature have also confirmed the reduction in gas breakthrough and gas channeling time(Li and Yu,2020),proving the reliability of the experimental study.

Fig.11.Relationship between oil recovery and injection volume(a)and back pressure(b).

Fig.12.Curves of gas-oil ratio versus injection volume(a),and gas breakthrough time and gas channeling time(b)at different back pressures.

Fig.13.Comparison of injection index at different injection volumes and back pressures.

The variation tendency of the injection index under different water saturations is shown in Fig.16.In the early stage of gas injection,the injection index is also different due to the difference in water saturation.In gas flooding with different water saturation,the injection index of high water saturation is larger than that of low water saturation,so the injection capacity of high water saturation is better.Due to the small amount of supercritical CO2dissolved in water,the supercritical CO2has a strong piston displacement effect.The liquid and gas at the outlet are observed earlier to maintain the injection capacity at a reasonable level.However,the gas-oil ratio rises faster for cores with low water saturation gas channeling,resulting in better supercritical CO2injection capacity and a more extensive injection index.In general,the supercritical CO2injection capacity under high water saturation in the early stage of gas injection is better.Still,the supercritical CO2injection capacity under low water saturation in the later gas injection stage is better.

4.Conclusions

(1)The gas breakthrough stage is where the oil production rate increases the fastest,however,the injection capacity at this stage is poor.The occurrence of gas channeling produces less oil,a sharp rise in the injection index,and better injection capacity.

(2)With an increase in the gas injection rate,the core pressure increases.At the same time,the dissolution and mixing of CO2with crude oil are enhanced.Moreover,the sweep efficiency of injected gas is improved,leading to a higher final recovery factor and a higher supercritical CO2consumption.The increase in injection rate delays the occurrence of gas channeling to an extent,but there is an optimal range of increase in injection rate.Adopting a high injection speed will cause the inlet pressure to rise,but this is not a sign of poor injection capacity.At this stage,it is necessary to evaluate the injection capacity quantitatively.

Fig.14.Relationship between oil recovery and injection volume(a)and water saturation(b).

Fig.15.Curves of gas-oil ratio versus injection volume at different saturation(a),and gas breakthrough time and gas channeling time(b)at different water saturation.

Fig.16.Comparison of injection index at different injection volumes and water saturation.

(3)CGI has low displacement efficiency,early gas channeling,more supercritical CO2consumption,and better injection capabilities.The AIP mode maximizes the use of supercritical CO2,improves the energy enhancement effect,expands the spread area of supercritical CO2,and improves oil recovery,but the injection-production capacity is poor.

(4)The high-permeability core has few tiny pores,and the diffusion rate of supercritical CO2is fast.Supercritical CO2produces most crude oil in large pores with high recovery and good injection capacity.Similarly,the solubility of supercritical CO2and crude oil is high under high pressure,changing the wettability of rock and improving crude oil recovery.The high-pressure gas injection contributes to late gas channeling time and poor injection capacity.

(5)The presence of water reduces the contact between supercritical CO2and crude oil under the condition of high water saturation.It reduces the diffusion of the CO2-crude oil mixture system,resulting in low supercritical CO2recovery,early gas invasion time,and good injection capacity.The gas-oil ratio with low water saturation rises sharply after gas channeling,but the injection capacity is adequately controlled and improved.Therefore,it is crucial to choose the appropriate gas injection occasion in the different stages of water flooding.

Declaration of competing interest

The authors declare no competing financial interest.

Acknowledgments

The authors would like to acknowledge financial support from the National Natural Science Foundation of China(No.51904324,51974348,U19B6003).The authors are grateful for Shengli Oilfield,providing rock and fluid data for our study.