• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Integration of a fused silica capillary and in-situ Raman spectroscopy for investigating CO2 solubility in n-dodecane at near-critical and supercritical conditions of CO2

    2022-02-17 17:52:26JunLingWngZiHoSongLinJunLiLiLiYngQunYunWngMingChouZhiYnPn
    Petroleum Science 2022年6期

    Jun-Ling Wng , Zi-Ho Song , Lin-Jun Li , Li-Li Yng , Qun-Yun Wng ,I-Ming Chou , Zhi-Yn Pn ,*

    a College of Environment, Zhejiang University of Technology, Hangzhou, Zhejiang, 310032, China

    b Hangzhou Academy of Environmental Protection, Hangzhou, Zhejiang, 310032, China

    c Laboratory for Experimental Study Under Deep-sea Extreme Conditions, Institute of Deep-sea Science and Engineering, Chinese Academy of Sciences,Sanya, Hainan, 572000, China

    A B S T R A C T To determine the solubility of CO2 in n-dodecane at T = 303.15-353.15 K, P ≤11.00 MPa, an integrated fused silica capillary and in-situ Raman spectroscopy system was built. The Raman peak intensity ratio(ICO2/IC-H)between the upper band of CO2 Fermi diad(ICO2)and the C-H stretching band of n-dodecane(IC-H)was employed to determine the solubility of CO2 in n-dodecane based on the calibrated correlation equation between the known CO2 molality in n-dodecane and the ICO2/IC-H ratio with R2 = 0.9998. The results indicated that the solubility of CO2 decreased with increasing temperature and increased with increasing pressure. The maximum CO2 molality (30.7314 mol/kg) was obtained at 303.15 K and 7.00 MPa. Finally, a solubility prediction model (lnS = (P - A)/B) based on the relationship with temperature (T in K) and pressure (P in MPa) was developed, where S is CO2 molality, A = - 8× 10-6T2+0.0354T- 8.1605,and B=0.0405T - 10.756.The results indicated that the solubilities of CO2 derived from this model were in good agreement with the experimental data.

    Keywords:CO2 solubility n-dodecane Raman spectroscopy Fused silica capillary In-situ

    1. Introduction

    Recently,mitigating climate change by reducing greenhouse gas emissions has become a global focus,and the ability to reduce CO2emissions accordingly has become an urgent task. Carbon capture and storage(CCS)has been widely recognized as a viable method to mitigate carbon emissions (Cormos, 2012; Peng et al., 2013; You et al., 2014; Li et al., 2016; Sun et al., 2017). Among all the storage methods, CO2geological storage has a direct effect on emission reduction and is considered to be the most potential storage method(Leung et al.,2014).Geological storage sites include oil and gas fields, saline aquifers, and deep coal seams (Michael et al.,2010). Injection of CO2in oil and gas fields can not only achieve CO2emission reductions but also enhance oil recovery (EOR) (Liu et al., 2016). Currently, many countries have conducted field tests of CO2flooding in oil fields using mature CO2-EOR technology(Awan et al.,2008;Hill et al.,2013;Lacy et al.,2013;Lv et al.,2015).

    The dissolution of CO2in oil can result in oil swelling and viscosity reduction, thereby enhancing oil recovery. CO2solubility data are one of the major parameters to determine the performance of CO2-EOR process (Mosavat et al., 2014). The solubility of CO2in crude or simulated oil (of which the primary components are alkanes) systems at different temperatures and pressures have been investigated previously.Chung et al.(1988)and Henni et al.(1996)employed gas chromatography to analyze the solubility of CO2in heavy oil at temperatures of 297,333 and 391 K and pressures of up to 34.5 MPa, and in n-dodecane at 313, 353, and 393 K and 0-9.6 MPa.Kavousi et al.(2014)studied CO2solubility in heavy oils with different viscosities using the pressure drop method at pressures of 1.73-4.48 MPa and temperatures of 295-305 K.

    Previous methods for investigating the solubility of CO2in crude/simulated oil systems at high pressure conditions have been performed primarily with the PVT apparatus (Lay et al., 2006;Nourozieh et al.,2013),equilibrium liquid sampling analysis(Yang et al., 2013), and chromatography (Forte et al., 2011). Although these methods offer distinct advantages,there continue to be some shortcomings.For example,the temperature and pressure ranges of the entire system of the PVT apparatus are limited, and the equilibrium liquid sampling analysis easily changes the original temperature and pressure conditions during the sampling process and upsets the system balance. Analysis of CO2quantity at the equilibrium point is the key to determine the solubility of CO2in crude oil and simulated crude oil. The equilibrium points of these methods generally are judged according to the pressure or the bubble point. It takes plenty of time to reach phase equilibrium because of the higher viscosity of the crude or simulated oil,which leads to an inaccurate judgment of the equilibrium condition in subjective or experimental environments(Mosavat et al.,2014;Han et al., 2015; Gui et al., 2017). In addition, the temperatures and pressures investigated in previous studies frequently have been confined to narrow ranges. Raman spectroscopy is an effective optical analysis method that is noninvasive,sensitive,and fast.(Liu et al.,2012;Belgodere et al.,2015).Therefore,Raman spectroscopy has been used in numerous studies of solubility measurements(e.g., Caumon et al., 2017; Wang et al.,2018).

    The change in Raman peak intensity(peak height or peak area)ratio between two species in a fluid reflects variations in their relative quantities. For example, Wang et al. (2019) and Guo et al.(2014) demonstrated that the ratio between the upper band of CO2Fermi diad and the O-H stretching band of water in homogenous solutions could be used to determine the solubility of CO2in water at various pressure-temperature conditions. Similarly, the ratio between the upper band of CO2Fermi diad and the C-H stretching band of n-dodecane enables quantification of CO2solubility in n-dodecane.

    In this study,to determine the solubility of CO2in n-dodecane at temperatures of 303.15-335.50 K and pressures up to 11.0 MPa,an integrated fused silica capillary and in-situ Raman spectroscopy system was built (Pan et al., 2013a, 2013b; Wang et al., 2017a,2017b). n-dodecane, one of the primary components of crude oil,was used as a simulated oil in this experiment. The physical properties of n-dodecane (ρ = 753 kg/m3,μ = 1.36 cp at 298.15 K,and atmospheric pressure)are similar to those of light crude oil(Bei et al.,2018).The Raman spectra of the system with CO2dissolved in n-dodecane are presented in Fig.1. It shows that the main Raman peaks (2800-3000 cm-1) of n-dodecane do not coincide with the CO2Fermi diad upper band (1385 cm-1) and its minor peaks at

    Fig.1. Raman spectra of CO2, n-dodecane, and fused silica capillary (FSC).

    1250-1500 cm-1do not block the upper Fermi diad band of CO2at 1385 cm-1.Therefore,any influence of the peaks of the n-dodecane at 1250-1500 cm-1to the intensity of the upper Fermi diad band of CO2at 1385 cm-1can be ignored. In this study, the ratio between the upper band of CO2Fermi diad and the C-H stretching band of n-dodecane(ICO2/IC-H)was used to determine the CO2solubility in n-dodecane,based on the calibrated correlation equation between the known CO2molality in n-dodecane and theICO2/IC-Hratio.Finally, a CO2solubility model based on its relationship with temperature and pressure was developed.

    2. Materials and methods

    2.1. Materials

    Technologies, LLC (Phoenix, AZ, USA) with an outer diameter of 665 μm and inner diameter of 300 μm. All of the valves and highpressure stainless-steel tubes were purchased from the Nantong Huaxing Petroleum Instrument Co., Ltd. (Nantong,Jiangsu, China).

    2.2. Apparatus

    An experimental setup was built to measure the solubility of CO2in n-dodecane at 303.15-353.15 K and up to 11.00 MPa(Fig.2).The apparatus primarily consisted of a fused silica capillary (FSC) and the balance kettle observation window, combined with a heatingcooling stage (CAP500, Linkam Scientific Instruments, Tadworth,UK), a confocal Raman spectrometer (Horiba JobinYvon, HR800,HORIBA FRANCE,Palaiseau,France),a horizontal phase equilibrium kettle with a magnetic stirrer, a manual pressure pump (JY-80,Jiangsu,China),a 70.00 MPa pressure sensor,a magnetic circulating pump,and a quantitative pump(30.00 MPa full scale).A circulation line was formed by the FSC connected to a circulating pump and a phase equilibrium kettle with an electric heating jacket, which could adjust the temperature of the fluid in the kettle. Controlled the temperature within the FSC utilizing a heating-cooling stage in conjunction with a digital temperature controller (accurate to 0.1 K). The pressure of the phase equilibrium kettle was adjusted using a manual pressure pump connected to a pressure transducer(70.00 MPa full scale,accurate to±0.25%FS)and an electric heater connected to a digital temperature controller (accurate to 0.1 K).Measurements of pressures in the quantitative pump were achieved by using a pressure transducer(30.00 MPa full scale,accurate to±0.25% FS).

    Fig. 2. Schematic diagram of the CO2 + n-dodecane solubility measurement system.1 to 13 are high-pressure valves. FSC = fused silica capillary.

    Fig. 3. Variation in the Raman peak intensity ratio of the CO2 + n-dodecane system with time (T = 293.15 K, P = 6.67 MPa, xCO2 = 0.1071 mol/kg).

    Fig. 4. The change in the CO2 + n-dodecane Raman peak intensity ratio with time under different CO2 concentrations at 303.15 K.

    Fig. 5. Changes in the CO2-saturated solution during heating and depressurization (a) →(d), and cooling and pressurization (e) →(h) in a FSC.

    Fig. 6. Three positions for Raman spectroscopic measurements in the CO2-saturated solution in a FSC.

    2.3. Experimental procedures

    The procedure for measuring the solubility of CO2in n-dodecane at different temperatures and pressures followed two steps: (1)determination of the relationship between the known CO2molality and the Raman peak intensity ratio of the upper band of CO2and ndodecane in a homogeneous CO2/n-dodecane system; and (2)measurement of the Raman peak intensity of the upper band of CO2and n-dodecane in a CO2-saturated CO2/n-dodecane system.

    The procedures for the first step included the following:(1)Valve 11 was opened,and a certain amount of n-dodecane was fed into the phase equilibrium kettle using the liquid sampler at roomtemperature(about 293.15 K).Then,the valve was closed;(2)Valves 2, 4,and 8 were opened, and the quantitative CO2, which was estimated in advance,was charged into the quantitative pump,and then the three valves were closed.The CO2was pressed into the cell using the quantitative pump, and valves 9 and 10 were successively opened and closed.During this step,a pressure transducer was used to measure the initial and final pressures,and a digital temperature controller was used to measure the initial and final temperatures.The void volume of the quantitative pump and the pipe gap between the quantitative pump and the phase equilibrium kettle was known.Thus, using the Peng-Robinson equation to calculate the amount of CO2which was injected into the system(accurate to±0.13%);(3)The mixture was slowly compressed to 30.00 MPa using a manual pressure pump that was connected to the phase equilibrium kettle.The stirring device was used for several hours to accelerate the CO2dissolution; (4) The circulating pump was turned on for approximately 30 min to ensure the homogeneity of the solution in the circulation line,and then it was turned off.The Raman spectrum of a homogeneous CO2/n-dodecane solution was collected;(5)Repeated step (4) until the deviation of the Raman peak intensity ratio between the upper band of CO2Fermi diad and the C-H stretching band of the n-dodecane from previous measurement could be neglected, proving that the system had reached dissolution equilibrium; (6) The temperature or pressure was changed using the same procedures to ensure that the system reached dissolution

    equilibrium,and the Raman spectra of the solution were collected at different pressures and temperatures. (7) The CO2molality was changed, and repeated steps (2)-(6); (8) The experiment was repeated to obtain the Raman spectra of the CO2-containing homogeneous solution under different pressures and at temperatures of 303.15-353.15 K.The CO2molality was plotted versus the Raman peak intensity ratio of the CO2/n-dodecane solution, and the equation of the calibrated line was obtained.

    The procedures for measuring the Raman peak intensity ratio of CO2and n-dodecane and judging the phase equilibrium in the saturated CO2/n-dodecane system are the same as those for the homogeneous CO2/n-dodecane system. After reaching phase equilibrium, some bubbles appeared as the temperature of the system was increased or as the pressure was decreased in the saturated CO2/n-dodecane system. The Raman peak intensities of the upper band of CO2Fermi diad and the C-H stretching band of the n-dodecane under different temperature-pressure conditions(303.15-353.15 K,0-15.00 MPa)were determined for the solution near the CO2bubbles(approximately 100 μm from the edge of the CO2bubbles),and the Raman peak intensity ratios were calculated.

    In the experiment,a JY/Horiba LabRam HR800 system equipped with a frequency-doubled Nd:YAG 531.95 nm laser with 20 mW of output laser power to collect the Raman spectra was employed.In addition,a charge coupled device(CCD)detector(multichannel,air cooled) was utilized to analyze the CO2+ n-dodecane system in situ. Under various temperature-pressure conditions, the Raman spectra of the system in the range of 1100-4000 cm-1were collected to obtain the peak intensity ratio of the CO2Fermi diad upper band to the C-H stretching band of the n-dodecane. The acquisition time was 20 s with two accumulations.The Raman peak intensity ratio of each group of temperature-pressure-composition(T-P-x) conditions was measured five times and the average value was invoked as the experimental datum.

    The details of how to determine the molality of CO2(x) have been reported in our previous study (Wang et al., 2019). Raman peak intensity (height) ratio (λ) was calculated using Eq. (1).

    whereICO2is the Raman peak intensity (height) of the CO2Fermi diad upper band,IC-His the Raman peak intensity (height) of the C-H stretching band of the n-dodecane, and the baseline is subtracted during the calculation.

    3. Results and discussion

    3.1. Determination of the phase equilibrium of the CO2 + ndodecane system

    The prerequisite for determining the solubility of CO2in ndodecane is the phase equilibrium of the CO2+n-dodecane system.Confocal Raman spectroscopy can be utilized to quantitatively monitor the concentration of gases in aqueous solutions because the band intensity of an active Raman species is proportional to its concentration in a fluid (Han et al., 2015; Wang et al., 2017a).Therefore, in this study, the peak intensity ratio of CO2Fermi diad upper band to the C-H stretching band of n-dodecane was employed to verify whether the system reached thermodynamic equilibrium under a particularT-P-xcondition.

    3.1.1. Phase equilibrium determination of the CO2+n-dodecane system during homogenization

    After feeding quantitative CO2into cells at room temperature,Raman spectroscopy was used to track the Raman peak intensity of the CO2and the n-dodecane in the FSC in real time and calculated the peak intensity ratio to understand the mass transfer of CO2in the system. Variation in the Raman peak intensity ratio of the CO2+ n-dodecane with time at 293 K and 6.67 MPa is shown in Fig. 3. It can be seen that the Raman peak intensity ratio was relatively high during the preliminary stage, but it decreased rapidly after the stirring device was turned on. During continuous stirring of the fluid,the CO2gradually dissolved in the n-dodecane,and the Raman peak intensity ratio increased and finally stabilized around 68 h, indicating the system had nearly reached phase equilibrium.

    After the system reached the phase equilibrium for this sample fluid with knownxCO2,the temperature of the system was adjusted and the Raman peak intensity of the system was measured at a newT-Pcondition. While the deviation of the obtained Raman peak intensity ratio from the previous measurement was negligible, it can be considered that the system has reached phase equilibrium at this time. Fig. 4 shows that the time required to reach the phase equilibrium was approximately 7 h at 303.15 K for CO2concentrations of 0.1071-0.2469 mol/kg.

    3.1.2. Phase equilibrium determination of a CO2-saturated solution

    While increasing the system temperature or decreasing the system pressure from the CO2-saturated equilibriumT-Pcondition,bubbles could be formed in the CO2+ n-dodecane system due to the amount of available CO2is more than the solubility value.Fig.5 shows the formation and growth of CO2bubble(s) during heating and decompression (a→d), and the shrinkage of a CO2bubble during cooling and pressurization (e→h). By adjusting the system to a specific temperature and pressure,the size of the bubbles in the CO2+ n-dodecane system can be observed under a microscope.When the sizes of bubbles did not change within 1 h, in order to further ensure the phase equilibrium of the system,Raman spectra at positions 1,2,and 3 were collected (Fig.6) and the Raman peak intensity ratios of the upper band of CO2Fermi diad and the C-H stretching band of n-dodecane were calculated. Results, listed in Table 2 for threeT-Pconditions,show that the system had reached phase equilibrium.

    Table 1 CAS registry number, source,mass fraction purity and molecular weight of the chemicals used in this work.

    Table 2 Raman peak intensity ratios of the CO2 + n-dodecane system at three positions in the FSC under three different temperature-pressure conditions.

    Table 3 The Raman peak intensity ratio of the CO2 + n-dodecane system for the same CO2 molality and different temperature-pressure conditions.a.

    To confirm the feasibility of the method,the peak intensities at position 2 were continuously tracked every 20 min for 2 h after phase equilibrium was reached to determine whether the phase balance was reached in this state. As shown in Fig. 7, the Raman spectra of the CO2+n-dodecane system were acquired at position 2 at 333.15 K and 5.55 MPa from 0 to 120 min, and the peak intensity ratios were found to be 0.2247, 0.2278, 0.2309,0.2313,and 0.2315. The variations in the peak intensity ratio became increasingly smaller as time went on,which were negligible during these times.

    Fig. 7. Raman spectra of the CO2 + n-dodecane system at 333.15 K and 5.55 MPa from 0 to 120 min after phase equilibrium was reached.

    Fig. 8. The relationship between CO2 molality and the Raman peak intensity ratio between the upper band of CO2 Fermi diad and n-dodecane derived from spectra collected from homogeneous CO2/n-dodecane solutions.

    3.2. The relationship between the CO2 concentration and the Raman peak intensity ratio in the CO2 + n-dodecane system

    After reaching phase equilibrium in a homogeneous CO2system,Raman spectroscopy was used to measure the peak intensity ratio of the system at different pressures and temperatures.As shown in Table 3,for the same CO2molality and temperatures,the changes in the Raman peak intensity ratio of the system with increasing pressure were very slight, with standard deviations of less than 0.0015, which were negligible. Thus, pressure had little effect on the peak intensity ratio. According to Table 4, the Raman peak intensity ratio of the system slightly decreased with increasing temperature at a constant CO2concentration and constant pressure, but the standard deviation of the peak intensity ratio from

    303.15 to 353.15 K was less than 0.0007 for eachxCO2, which indicated that the temperature also had little effect on the peak intensity ratio. As a result, the peak intensity ratio of the homogeneous solution did not depend on temperature and pressure within the range of these experimental conditions.Therefore,the linear relationship between the CO2molality and the CO2+ndodecane Raman intensity ratio does not need to consider the effect of temperature and pressure.

    In this experiment, the average value of the Raman peak intensity ratio of the CO2+ n-dodecane system at the same CO2molality and under design temperature conditions was used to establish a correlation equation for CO2molality.As shown in Fig.8(the data see Table S1), there was a positive correlation(R2= 0.9998) between the CO2molality and the Raman peak intensity ratio of the CO2+ n-dodecane system. The correlation equation is as following:

    Table 4 Raman peak intensity ratio of a homogeneous solution at different temperatures and CO2 molality.

    Table 5 Raman peak intensity ratio and solubility of the CO2 in n-dodecane at 303.15-353.15 K and 1.00-11.00 MPaa.

    Fig. 9. Raman peak intensity ratios of the CO2 + n-dodecane system at 303.15-353.15 K and 1.00-11.00 MPa.

    where λ is the peak intensity ratio of the CO2+n-dodecane system;andxis the molality of CO2(mol/kg).

    Fig.10. Solubility of CO2 in n-dodecane at 303.15-353.15 K and 1.00-11.00 MPa.

    Fig.11. Comparison of the CO2 solubilities in n-dodecane obtained in this study with those reported previously.

    3.3. Determination of the solubility of CO2 in n-dodecane

    The solubility of the system based on the peak intensity ratio of the system under different temperature and pressure conditions in the saturated solution could be calculated from the relationship between the CO2molality and the peak intensity ratio of the CO2+ n-dodecane (Eq. (2)). At temperatures of 303.15-353.15 K and pressures of 1.00-11.00 MPa,the Raman spectra of CO2and ndodecane in CO2-saturated solution were collected, and the peak intensity ratios were calculated (shown in Table 5 and plotted in Fig. 9). The reproducibility of peak intensity ratio is characterized by relative standard deviation,which ranges from 0.23%to 0.90%,it proved that samples were appropriate to be analyzed by this method with satisfactory results in our works.

    As shown in Fig.10,the solubility of the CO2-saturated solution decreased with increasing temperature and increased with increasing pressure.The maximum solubility obtained in this study is 30.7314 mol/kg at 303.15 K and 7.00 MPa.

    According to Fig.10,the CO2solubility in n-dodecane gradually increased with increasing pressure and decreased with increasing temperature. At lower pressure, CO2had a weaker solubility in ndodecane, and the solubility increased nearly linearly with increasing pressure. As the pressure increased, the distance between the CO2and the n-dodecane molecule decreased,and the interaction force increased and gradually approached the force between the n-dodecane molecules, causing the CO2solubility to increase (Yang et al., 2013). In particular, the CO2solubility in ndodecane-rich liquid phase changes significantly for pressures and temperatures close to the critical pressure and temperature of pure CO2. For example, the pressure increased from 1.08 MPa to 5.03 MPa at 303.15 K, whereas the solubility increased only from 0.6331 mol/kg to 7.9775 mol/kg.In addition,the pressure increased from 5.03 MPa to 7.00 MPa, whereas the solubility increased dramatically from 7.9775 mol/kg to 30.7314 mol/kg.This indicated that the CO2+n-dodecane system reached the miscible phase with increasing pressure, thereby making the solubility value change significantly with increasing pressure(Lashkarbolooki et al.,2017).In contrast, at lower pressure, the CO2solubility in n-dodecane decreased with increasing temperature.The temperature increased from 303.15 to 353.15 K at 5.03 MPa,resulting in a decrease in the solubility of CO2in n-dodecane from 7.9775 mol/kg to 1.6270 mol/kg. Studies have shown that increasing temperature can increase the distance between the CO2and alkane molecules, reduce their interactions, enhance the Brownian motion of the CO2molecules,and allow the CO2molecules to escape from the liquid phase,which reduces the CO2solubility (Yang et al., 2012).

    The CO2solubilities in n-dodecane obtained in this study were compared with those reported previously (Fig.11). Our CO2solubility data were basically consistent with the experimental data of Henni et al. (1996) and Brain Bufkin (1986) at pressures below 4.00 MPa.However,at higher pressures,our data show higher CO2solubility than those published in the literature,except at 353.15 K.The reasons of lower solubility may be that the system did not reach phase equilibrium while analyzing the results during their experiments.Another reason may be that previous studies need to extract sample fluids, which may disrupt the original equilibrium condition of the system,lead to the escape of CO2. In addition,the CO2solubility becomes more sensitive to the changes of pressure and temperature when theP-Tconditions are near the critical point of CO2.

    Fig.12. Comparison of the experimental values of CO2 solubility in n-dodecane with those calculated from a model.

    3.4. Solubility prediction model for a CO2 + n-dodecane system

    In CO2-EOR technology, CO2solubility is one of the most important factors required to improve oil recovery. Because of a lack of experimental data for CO2solubility in a wider range of temperatures and pressures, semi-empirical models often have been adopted in previous studies(Tsuji et al.,2004)to predict CO2solubility in organic liquids. Jou and Mather (2005) fitted a CO2solubility model based on experimental data, as shown in Eq. (3).

    wherePis the pressure(MPa);Sis the CO2solubility(mol/kg);andAandBare functions of the temperatureT.

    The model is relatively simple and suitable for predicting the solubility of CO2in different solvents.Fornari et al.(2009),Paninho et al. (2013), and Howlader et al. (2017) used this model to determine the solubility of CO2in organic liquids, such as ethyl lactate and triglycerides. The prediction results were in good agreement with the experimental results.Therefore,this model was also used in this work to predict the solubility of CO2in n-dodecane.

    In this study, the model developed by Jou and Mather (2005)was improved, as shown in Eq. (4). First, the experimental CO2solubilities were fitted with pressure,P,under isothermal conditions.Then,based on the improved model,the relationships amongA,B,andTwere obtained,as shown in Eqs.(5)and(6).Thus,the CO2solubility correlation model obtained by correlating the temperature and pressure was determined to be as following:

    The CO2solubility calculated using the prediction model(Fig.12)were consistent with the experimental data,and the mean average percentage error between the experimental and calculated data was within 4.0%. Therefore, the model can be expanded in future work to predict the CO2solubility in a wider range of temperatures and pressures.

    4. Conclusions

    In this study, to determine the solubility of CO2in n-dodecane under near-critical and supercritical conditions(T= 303.15-353.15 K,P≤11.00 MPa), an integrated fused silica capillary and in-situ Raman spectroscopy system was built. The calculated Raman peak intensity ratio was used to qualitatively and quantitatively analyze the CO2and its concentrations. The relationship between the known CO2molality in the n-dodecane and the Raman intensity ratio,ICO2/IC-H, was used to establish the correlation equation.Subsequently,the CO2solubility of the system was calculated by incorporating the measuredICO2/IC-Hratio into the correlation equation. The maximum CO2solubility(30.7314 mol/kg)was obtained at 303.15 K and 7.00 MPa.Finally,a solubility correlation model (lnS= (P-A)/B) based on the relationship between temperature and pressure was developed. The CO2solubility data can be extended using the model equation to provide an empirical correlation and experimental data for the CO2-EOR technology. This is a promising method for the determination of CO2solubility in organic liquids under supercritical conditions of CO2,provided the occasionally encountered interference of fluorescence on Raman signals can be effectively reduced.

    Declaration of competing interest

    The authors declare no competing financial interests.

    Acknowledgment

    This work was supported by the National Key Research and Development Program of China(2019YFE0117200)and the Natural Science Foundation of China (41977304). We thank the editors of Petroleum Science and the reviewers, whose constructive comments improve the paper.

    Appendix A. Supplementary data

    Supplementary data to this article can be found online at https://doi.org/10.1016/j.petsci.2022.06.014.

    久久久国产欧美日韩av| 久久香蕉国产精品| 久久久久久人人人人人| 97人妻天天添夜夜摸| xxxhd国产人妻xxx| 日韩有码中文字幕| 久久久久久久精品吃奶| 成人亚洲精品一区在线观看| 国产蜜桃级精品一区二区三区 | 精品国产一区二区三区久久久樱花| 老司机影院毛片| 国产三级黄色录像| 中文字幕人妻丝袜一区二区| 黄片小视频在线播放| 老熟女久久久| 精品亚洲成国产av| 日本欧美视频一区| 新久久久久国产一级毛片| www日本在线高清视频| 正在播放国产对白刺激| tocl精华| 国产精品久久久久久人妻精品电影| 大香蕉久久网| 欧美午夜高清在线| 婷婷精品国产亚洲av在线 | 成人免费观看视频高清| 免费少妇av软件| 乱人伦中国视频| 老汉色∧v一级毛片| 亚洲专区国产一区二区| 丁香六月欧美| 国产99久久九九免费精品| 中文字幕色久视频| e午夜精品久久久久久久| 亚洲五月婷婷丁香| 男女下面插进去视频免费观看| 极品教师在线免费播放| 91精品三级在线观看| 日本精品一区二区三区蜜桃| 欧美黑人精品巨大| 国产精品国产高清国产av | 亚洲 国产 在线| 丰满饥渴人妻一区二区三| 天堂俺去俺来也www色官网| 国产精品亚洲一级av第二区| 中文字幕制服av| 超碰97精品在线观看| 亚洲中文字幕日韩| 日本黄色视频三级网站网址 | 51午夜福利影视在线观看| 大香蕉久久网| 国产一区在线观看成人免费| 日韩视频一区二区在线观看| 天堂动漫精品| 久久精品国产清高在天天线| av视频免费观看在线观看| 黄频高清免费视频| 他把我摸到了高潮在线观看| 午夜福利欧美成人| 91成人精品电影| 久久久国产欧美日韩av| 在线视频色国产色| 久久国产乱子伦精品免费另类| av线在线观看网站| 少妇粗大呻吟视频| 国产精品免费视频内射| 好看av亚洲va欧美ⅴa在| 丁香六月欧美| av国产精品久久久久影院| 国产精品一区二区在线不卡| 亚洲精品国产一区二区精华液| 亚洲成国产人片在线观看| 亚洲欧美激情综合另类| 香蕉国产在线看| 欧美亚洲 丝袜 人妻 在线| 麻豆av在线久日| 亚洲黑人精品在线| 久久99一区二区三区| 性色av乱码一区二区三区2| 久久精品人人爽人人爽视色| 啪啪无遮挡十八禁网站| 亚洲人成电影观看| 自拍欧美九色日韩亚洲蝌蚪91| xxxhd国产人妻xxx| 激情在线观看视频在线高清 | 欧美黑人精品巨大| 国产精品九九99| 欧美激情极品国产一区二区三区| 国产精品.久久久| 在线看a的网站| 一边摸一边抽搐一进一小说 | 九色亚洲精品在线播放| 少妇 在线观看| 男女高潮啪啪啪动态图| 欧美黄色淫秽网站| 日韩大码丰满熟妇| 中文欧美无线码| av福利片在线| 18禁黄网站禁片午夜丰满| 在线观看日韩欧美| 在线观看免费午夜福利视频| 99热国产这里只有精品6| 成人黄色视频免费在线看| 两性夫妻黄色片| 欧美乱妇无乱码| 欧美亚洲 丝袜 人妻 在线| 大香蕉久久成人网| 精品午夜福利视频在线观看一区| 制服人妻中文乱码| 国产免费男女视频| 自线自在国产av| 18禁国产床啪视频网站| 99国产综合亚洲精品| 高清毛片免费观看视频网站 | 亚洲男人天堂网一区| 国产亚洲精品久久久久5区| cao死你这个sao货| 国产熟女午夜一区二区三区| 欧美色视频一区免费| 又大又爽又粗| 香蕉国产在线看| 欧美丝袜亚洲另类 | 久久香蕉国产精品| 80岁老熟妇乱子伦牲交| 深夜精品福利| 亚洲九九香蕉| 啪啪无遮挡十八禁网站| 久久香蕉激情| 国产单亲对白刺激| 淫妇啪啪啪对白视频| 丝袜美足系列| 国产精品98久久久久久宅男小说| 999精品在线视频| 国产精品久久久人人做人人爽| 免费高清在线观看日韩| 99精国产麻豆久久婷婷| 精品福利永久在线观看| 很黄的视频免费| 久久人人爽av亚洲精品天堂| 国产99久久九九免费精品| 人人澡人人妻人| 欧美中文综合在线视频| 欧美成人免费av一区二区三区 | 亚洲色图综合在线观看| 色综合婷婷激情| 老司机深夜福利视频在线观看| 一本大道久久a久久精品| 免费日韩欧美在线观看| 9热在线视频观看99| 欧美 亚洲 国产 日韩一| 不卡av一区二区三区| a级毛片黄视频| 亚洲专区中文字幕在线| 国产精品免费大片| √禁漫天堂资源中文www| 成熟少妇高潮喷水视频| 最新美女视频免费是黄的| xxxhd国产人妻xxx| 亚洲中文日韩欧美视频| 美女福利国产在线| 女人爽到高潮嗷嗷叫在线视频| 高清欧美精品videossex| 一本综合久久免费| 亚洲国产精品合色在线| 午夜两性在线视频| 动漫黄色视频在线观看| 丝袜在线中文字幕| 老熟妇仑乱视频hdxx| 好男人电影高清在线观看| 久久精品国产a三级三级三级| 国产日韩欧美亚洲二区| 巨乳人妻的诱惑在线观看| 一区福利在线观看| 午夜福利欧美成人| av视频免费观看在线观看| 国产一卡二卡三卡精品| 天堂√8在线中文| 国产三级黄色录像| 国产精品自产拍在线观看55亚洲 | 黄色怎么调成土黄色| 成人精品一区二区免费| 午夜福利影视在线免费观看| 久久久久国产精品人妻aⅴ院 | 黄频高清免费视频| 亚洲黑人精品在线| 如日韩欧美国产精品一区二区三区| 国产成人免费观看mmmm| 777米奇影视久久| 日韩欧美在线二视频 | 亚洲第一av免费看| 在线永久观看黄色视频| 美女高潮喷水抽搐中文字幕| 五月开心婷婷网| 欧美激情高清一区二区三区| 久久亚洲真实| 欧美精品一区二区免费开放| 99香蕉大伊视频| 操美女的视频在线观看| av一本久久久久| 午夜免费鲁丝| 午夜91福利影院| 变态另类成人亚洲欧美熟女 | 99riav亚洲国产免费| 一区二区三区精品91| 热99re8久久精品国产| 成年女人毛片免费观看观看9 | 亚洲午夜精品一区,二区,三区| 动漫黄色视频在线观看| 日韩 欧美 亚洲 中文字幕| 少妇被粗大的猛进出69影院| 黄色视频不卡| 在线看a的网站| 青草久久国产| 成人国语在线视频| 国产成人啪精品午夜网站| 9色porny在线观看| 在线观看舔阴道视频| 国产精品1区2区在线观看. | 国产精品久久久人人做人人爽| 91精品三级在线观看| 亚洲欧美一区二区三区黑人| 又大又爽又粗| 精品久久久久久电影网| 多毛熟女@视频| 免费在线观看日本一区| 亚洲自偷自拍图片 自拍| 建设人人有责人人尽责人人享有的| 伊人久久大香线蕉亚洲五| av在线播放免费不卡| 日韩三级视频一区二区三区| 亚洲七黄色美女视频| 两性夫妻黄色片| 欧美一级毛片孕妇| 国产激情欧美一区二区| 水蜜桃什么品种好| 国产亚洲精品久久久久久毛片 | 女人精品久久久久毛片| 免费一级毛片在线播放高清视频 | 91av网站免费观看| 久久久国产成人免费| 伦理电影免费视频| 99国产精品一区二区三区| 黑人巨大精品欧美一区二区蜜桃| 久热这里只有精品99| 久久人妻福利社区极品人妻图片| 亚洲精华国产精华精| 夜夜夜夜夜久久久久| 黄色a级毛片大全视频| 国产精品99久久99久久久不卡| 成人国产一区最新在线观看| 国产精品国产av在线观看| 亚洲av熟女| 亚洲人成电影免费在线| www.精华液| 亚洲中文日韩欧美视频| a级毛片黄视频| 久久精品国产a三级三级三级| 亚洲,欧美精品.| 亚洲av第一区精品v没综合| 久久午夜亚洲精品久久| 久久久久国产一级毛片高清牌| 看免费av毛片| 欧美日韩精品网址| 国产又爽黄色视频| 午夜福利免费观看在线| 久久久久国内视频| 国产97色在线日韩免费| 国产精品一区二区免费欧美| 国产精品乱码一区二三区的特点 | 青草久久国产| 成熟少妇高潮喷水视频| 99国产精品99久久久久| 精品国产超薄肉色丝袜足j| av欧美777| 欧美日韩亚洲综合一区二区三区_| 如日韩欧美国产精品一区二区三区| 老司机福利观看| 久久国产亚洲av麻豆专区| 午夜激情av网站| 色在线成人网| 国产午夜精品久久久久久| 国产一区二区三区综合在线观看| 久久午夜亚洲精品久久| 嫁个100分男人电影在线观看| 国产欧美日韩一区二区精品| 欧美久久黑人一区二区| 欧美日韩亚洲高清精品| 捣出白浆h1v1| 99热国产这里只有精品6| 成人精品一区二区免费| 无遮挡黄片免费观看| а√天堂www在线а√下载 | 国产在线一区二区三区精| 极品教师在线免费播放| 精品久久蜜臀av无| 水蜜桃什么品种好| 午夜福利,免费看| 三级毛片av免费| 99久久99久久久精品蜜桃| 国产亚洲av高清不卡| 在线免费观看的www视频| 亚洲欧美一区二区三区黑人| 亚洲第一欧美日韩一区二区三区| 99精品在免费线老司机午夜| 麻豆乱淫一区二区| 一进一出好大好爽视频| 亚洲七黄色美女视频| 中出人妻视频一区二区| 黄色毛片三级朝国网站| 天天躁夜夜躁狠狠躁躁| 免费在线观看完整版高清| 国产伦人伦偷精品视频| 精品久久蜜臀av无| 欧美乱码精品一区二区三区| 在线视频色国产色| 久久精品人人爽人人爽视色| 夜夜夜夜夜久久久久| 9191精品国产免费久久| 天天躁夜夜躁狠狠躁躁| 91av网站免费观看| 色婷婷av一区二区三区视频| 99久久99久久久精品蜜桃| 日韩制服丝袜自拍偷拍| 国产欧美亚洲国产| 国产欧美日韩一区二区精品| 国产精品免费大片| 成人特级黄色片久久久久久久| 十八禁网站免费在线| 久99久视频精品免费| 99久久99久久久精品蜜桃| 久久久国产欧美日韩av| 日本a在线网址| 欧美另类亚洲清纯唯美| 成年版毛片免费区| 少妇被粗大的猛进出69影院| 青草久久国产| 国产成人免费无遮挡视频| 新久久久久国产一级毛片| 91av网站免费观看| 美女扒开内裤让男人捅视频| 侵犯人妻中文字幕一二三四区| 亚洲国产毛片av蜜桃av| 午夜亚洲福利在线播放| 狂野欧美激情性xxxx| av不卡在线播放| 亚洲国产精品sss在线观看 | 亚洲成人免费av在线播放| 中文字幕另类日韩欧美亚洲嫩草| 中文字幕色久视频| 欧美国产精品va在线观看不卡| 午夜福利视频在线观看免费| 国产精品成人在线| 脱女人内裤的视频| 午夜福利乱码中文字幕| 欧美人与性动交α欧美精品济南到| 午夜福利影视在线免费观看| 亚洲全国av大片| 19禁男女啪啪无遮挡网站| 99香蕉大伊视频| 国产精品久久久久成人av| 美女高潮喷水抽搐中文字幕| 国产成人av教育| 黑人操中国人逼视频| 久久精品国产a三级三级三级| 国产精品一区二区在线观看99| 夫妻午夜视频| 美女国产高潮福利片在线看| 国产精品久久久人人做人人爽| 亚洲一区二区三区欧美精品| 国产片内射在线| 91字幕亚洲| 国内毛片毛片毛片毛片毛片| av网站在线播放免费| 午夜激情av网站| 1024香蕉在线观看| 香蕉国产在线看| 成人精品一区二区免费| 少妇被粗大的猛进出69影院| 国产精华一区二区三区| 99精品在免费线老司机午夜| 久久久久久久国产电影| 色94色欧美一区二区| 日本撒尿小便嘘嘘汇集6| 男女下面插进去视频免费观看| 亚洲综合色网址| 九色亚洲精品在线播放| 日韩欧美三级三区| 久久香蕉激情| 亚洲国产欧美网| 在线视频色国产色| 久久久久久久精品吃奶| 搡老乐熟女国产| 999久久久精品免费观看国产| 亚洲精品美女久久久久99蜜臀| 中文字幕人妻熟女乱码| 精品国产国语对白av| 日韩成人在线观看一区二区三区| 亚洲一码二码三码区别大吗| 国产精品免费大片| 黄网站色视频无遮挡免费观看| 91精品三级在线观看| 国产精品美女特级片免费视频播放器 | 日韩欧美国产一区二区入口| 精品久久久久久久毛片微露脸| 亚洲免费av在线视频| 免费日韩欧美在线观看| 久热这里只有精品99| 久久午夜亚洲精品久久| 真人做人爱边吃奶动态| 亚洲成人免费av在线播放| 国产成人精品无人区| 高清在线国产一区| 999精品在线视频| 亚洲性夜色夜夜综合| 成熟少妇高潮喷水视频| 麻豆乱淫一区二区| 女人久久www免费人成看片| 首页视频小说图片口味搜索| 亚洲色图av天堂| 99热只有精品国产| 一进一出抽搐动态| 又黄又粗又硬又大视频| 日韩欧美一区视频在线观看| 十八禁人妻一区二区| 99精品在免费线老司机午夜| 女人高潮潮喷娇喘18禁视频| cao死你这个sao货| 91大片在线观看| 男人舔女人的私密视频| 天堂中文最新版在线下载| 久久ye,这里只有精品| 午夜福利视频在线观看免费| 亚洲熟女精品中文字幕| 一级a爱视频在线免费观看| 露出奶头的视频| 久久久久精品国产欧美久久久| 一级a爱片免费观看的视频| 成人手机av| 色婷婷久久久亚洲欧美| 大香蕉久久网| 在线播放国产精品三级| 伊人久久大香线蕉亚洲五| 国产成人影院久久av| 精品高清国产在线一区| 国产免费av片在线观看野外av| 十分钟在线观看高清视频www| 国产亚洲av高清不卡| videosex国产| 亚洲精品国产色婷婷电影| 国产在线一区二区三区精| 国产精华一区二区三区| av不卡在线播放| 亚洲精品粉嫩美女一区| 国产精品免费大片| av线在线观看网站| 桃红色精品国产亚洲av| 亚洲一区二区三区不卡视频| 人人妻人人爽人人添夜夜欢视频| 国产亚洲欧美98| 国产视频一区二区在线看| 欧美激情极品国产一区二区三区| 大片电影免费在线观看免费| 国产亚洲av高清不卡| 亚洲欧美一区二区三区黑人| 日日摸夜夜添夜夜添小说| 黄频高清免费视频| 999久久久精品免费观看国产| 日韩视频一区二区在线观看| 两个人看的免费小视频| 久久精品熟女亚洲av麻豆精品| 国产成人av教育| 超碰97精品在线观看| 丝袜人妻中文字幕| 亚洲午夜精品一区,二区,三区| 老司机影院毛片| 丰满迷人的少妇在线观看| 久久精品aⅴ一区二区三区四区| av片东京热男人的天堂| 午夜激情av网站| 精品久久久久久久毛片微露脸| xxxhd国产人妻xxx| 精品视频人人做人人爽| 黄色片一级片一级黄色片| 丰满饥渴人妻一区二区三| 18禁美女被吸乳视频| 亚洲av成人一区二区三| 捣出白浆h1v1| 国产在线精品亚洲第一网站| 黄色女人牲交| 多毛熟女@视频| 国产一区有黄有色的免费视频| 亚洲欧美激情综合另类| 日韩制服丝袜自拍偷拍| 午夜精品在线福利| 久久香蕉精品热| 51午夜福利影视在线观看| 高潮久久久久久久久久久不卡| 亚洲美女黄片视频| 国产精华一区二区三区| 欧美乱色亚洲激情| 十八禁网站免费在线| 国产精品久久电影中文字幕 | 国产又爽黄色视频| 国产成人免费观看mmmm| 亚洲精品国产区一区二| 一级,二级,三级黄色视频| 亚洲专区字幕在线| 99精国产麻豆久久婷婷| 可以免费在线观看a视频的电影网站| 少妇粗大呻吟视频| 亚洲五月婷婷丁香| 久热爱精品视频在线9| 国产亚洲精品久久久久5区| 色综合婷婷激情| 国产男女超爽视频在线观看| 亚洲情色 制服丝袜| 日韩欧美一区二区三区在线观看 | svipshipincom国产片| 天天躁狠狠躁夜夜躁狠狠躁| 国产成人av教育| 身体一侧抽搐| 一级片免费观看大全| 香蕉久久夜色| 日本黄色日本黄色录像| 高清在线国产一区| 亚洲国产精品sss在线观看 | 999久久久精品免费观看国产| 亚洲人成77777在线视频| 免费久久久久久久精品成人欧美视频| 国产欧美日韩一区二区三区在线| av超薄肉色丝袜交足视频| 中文字幕av电影在线播放| 欧美老熟妇乱子伦牲交| 少妇的丰满在线观看| 国产精品久久久久久精品古装| 天天操日日干夜夜撸| 99精品在免费线老司机午夜| 国产精品免费一区二区三区在线 | 久久久国产成人精品二区 | 亚洲情色 制服丝袜| 免费看a级黄色片| 欧美日韩视频精品一区| 欧美日韩瑟瑟在线播放| 成人免费观看视频高清| 搡老乐熟女国产| 久久国产乱子伦精品免费另类| 欧美色视频一区免费| 手机成人av网站| 美女高潮到喷水免费观看| 久久久国产一区二区| 婷婷成人精品国产| 老鸭窝网址在线观看| 黑人操中国人逼视频| 亚洲精品美女久久久久99蜜臀| 成人影院久久| 亚洲第一av免费看| 亚洲欧美激情综合另类| 欧美亚洲日本最大视频资源| 亚洲aⅴ乱码一区二区在线播放 | 国产成人啪精品午夜网站| 国产一区在线观看成人免费| 成熟少妇高潮喷水视频| 国产男靠女视频免费网站| 99久久国产精品久久久| 精品电影一区二区在线| 久久久国产成人免费| 亚洲精品中文字幕在线视频| 国产精华一区二区三区| 黑人巨大精品欧美一区二区mp4| 国产成人精品无人区| 欧美精品啪啪一区二区三区| 精品国产乱码久久久久久男人| 在线观看舔阴道视频| 人人妻人人澡人人爽人人夜夜| 一级毛片精品| 国产一区二区激情短视频| 亚洲 欧美一区二区三区| 一二三四社区在线视频社区8| xxx96com| 人人妻人人爽人人添夜夜欢视频| 中文字幕最新亚洲高清| 亚洲av美国av| 黄色怎么调成土黄色| 日本黄色日本黄色录像| 亚洲精品av麻豆狂野| 国产精品九九99| 日日摸夜夜添夜夜添小说| 成年人免费黄色播放视频| xxxhd国产人妻xxx| 国产欧美日韩精品亚洲av| 淫妇啪啪啪对白视频| 久久久久精品人妻al黑| 国精品久久久久久国模美| 美女高潮到喷水免费观看| 一区二区三区国产精品乱码| 亚洲av熟女| 岛国在线观看网站| 母亲3免费完整高清在线观看| 一二三四在线观看免费中文在| 久久久久国内视频| 热re99久久精品国产66热6| 如日韩欧美国产精品一区二区三区| 女人被躁到高潮嗷嗷叫费观| 无人区码免费观看不卡| 午夜免费成人在线视频| 狂野欧美激情性xxxx| 国产淫语在线视频| 国产精品永久免费网站| 国产欧美日韩一区二区三| 黄片播放在线免费| 欧美日韩中文字幕国产精品一区二区三区 | 久久精品国产99精品国产亚洲性色 | 国产欧美日韩一区二区三区在线| 精品熟女少妇八av免费久了| 99久久人妻综合|