• <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二区三区四区| 欧美日本中文国产一区发布| 精品酒店卫生间| 亚洲电影在线观看av| 亚洲国产精品一区二区三区在线| 只有这里有精品99| 色婷婷av一区二区三区视频| 欧美精品国产亚洲| 国产高清国产精品国产三级| 免费观看无遮挡的男女| 国语对白做爰xxxⅹ性视频网站| 国产在视频线精品| 国产精品一区www在线观看| 少妇人妻精品综合一区二区| 成年av动漫网址| 狂野欧美白嫩少妇大欣赏| 欧美3d第一页| tube8黄色片| av专区在线播放| 国产精品99久久久久久久久| 亚洲第一区二区三区不卡| 亚洲精品视频女| 亚洲天堂av无毛| 亚洲伊人久久精品综合| √禁漫天堂资源中文www| 成人漫画全彩无遮挡| 成人国产av品久久久| 搡女人真爽免费视频火全软件| 精品国产乱码久久久久久小说| 九九久久精品国产亚洲av麻豆| 18禁裸乳无遮挡动漫免费视频| 国产亚洲91精品色在线| 色视频在线一区二区三区| 18禁在线无遮挡免费观看视频| 国产老妇伦熟女老妇高清| 麻豆乱淫一区二区| 亚洲国产色片| 亚洲美女搞黄在线观看| 国产欧美日韩精品一区二区| 夫妻性生交免费视频一级片| 日韩人妻高清精品专区| 免费播放大片免费观看视频在线观看| 91久久精品国产一区二区三区| 日本爱情动作片www.在线观看| av免费观看日本| 亚洲国产av新网站| 成人18禁高潮啪啪吃奶动态图 | 人妻一区二区av| 欧美老熟妇乱子伦牲交| 大话2 男鬼变身卡| 国产欧美日韩综合在线一区二区 | 国产老妇伦熟女老妇高清| 乱系列少妇在线播放| 另类精品久久| 十八禁高潮呻吟视频 | 爱豆传媒免费全集在线观看| 国产精品女同一区二区软件| 一区二区三区免费毛片| 久久精品久久久久久久性| 色视频www国产| 精品亚洲成a人片在线观看| 日韩欧美一区视频在线观看 | 国产片特级美女逼逼视频| 99精国产麻豆久久婷婷| 少妇高潮的动态图| av视频免费观看在线观看| 女性被躁到高潮视频| 水蜜桃什么品种好| 亚洲人与动物交配视频| 在线 av 中文字幕| 最近中文字幕高清免费大全6| 国产精品无大码| 大香蕉久久网| 伦精品一区二区三区| 蜜桃在线观看..| 国产欧美日韩精品一区二区| 蜜桃在线观看..| 亚洲精品一二三| 久久精品国产亚洲av涩爱| 久久这里有精品视频免费| 老女人水多毛片| 如何舔出高潮| 少妇 在线观看| 国产男女内射视频| 精品久久久久久电影网| 日韩伦理黄色片| 日韩制服骚丝袜av| 国产成人精品婷婷| 亚洲国产欧美在线一区| 视频中文字幕在线观看| 男男h啪啪无遮挡| 国产精品久久久久久av不卡| av女优亚洲男人天堂| 男人狂女人下面高潮的视频| av在线app专区| 国产精品国产av在线观看| 不卡视频在线观看欧美| 欧美激情极品国产一区二区三区 | 午夜精品国产一区二区电影| 亚洲,欧美,日韩| 人妻夜夜爽99麻豆av| av福利片在线| 91久久精品国产一区二区三区| 久久午夜福利片| 国产视频内射| 日韩一区二区视频免费看| 看十八女毛片水多多多| 永久网站在线| 国产极品粉嫩免费观看在线 | 22中文网久久字幕| 最近最新中文字幕免费大全7| 国产精品女同一区二区软件| 熟女人妻精品中文字幕| 国产淫语在线视频| 日本黄大片高清| www.色视频.com| 18禁在线播放成人免费| 久久精品久久久久久久性| 亚洲丝袜综合中文字幕| 三级国产精品欧美在线观看| 十分钟在线观看高清视频www | 国产精品女同一区二区软件| 色视频www国产| h视频一区二区三区| 美女脱内裤让男人舔精品视频| 日本与韩国留学比较| 中文精品一卡2卡3卡4更新| 久久午夜福利片| 极品教师在线视频| 国产欧美日韩精品一区二区| 久久久久久久久久久久大奶| 丰满乱子伦码专区| 高清欧美精品videossex| 九九爱精品视频在线观看| 国产伦理片在线播放av一区| 日韩av不卡免费在线播放| 亚洲一级一片aⅴ在线观看| 一级二级三级毛片免费看| 亚洲国产精品成人久久小说| 亚洲第一av免费看| 成年美女黄网站色视频大全免费 | 尾随美女入室| 中文天堂在线官网| 国产精品一区二区性色av| 少妇人妻一区二区三区视频| 9色porny在线观看| 在线 av 中文字幕| 国产成人一区二区在线| 婷婷色麻豆天堂久久| 一本色道久久久久久精品综合| 99九九线精品视频在线观看视频| 熟女人妻精品中文字幕| 人妻 亚洲 视频| 在线天堂最新版资源| 五月开心婷婷网| 天天躁夜夜躁狠狠久久av| 亚洲怡红院男人天堂| 天堂8中文在线网| 国产老妇伦熟女老妇高清| 久久久久久久国产电影| 国产伦在线观看视频一区| 中文资源天堂在线| 肉色欧美久久久久久久蜜桃| 美女国产视频在线观看| 一级片'在线观看视频| 国产成人a∨麻豆精品| 国产在视频线精品| 午夜激情福利司机影院| 亚洲三级黄色毛片| 色5月婷婷丁香| 国产av码专区亚洲av| 中文欧美无线码| 亚洲精品国产色婷婷电影| 成年人午夜在线观看视频| 国产精品久久久久久精品电影小说| 亚洲精品中文字幕在线视频 | 亚洲精品乱码久久久v下载方式| 又爽又黄a免费视频| 日本黄色日本黄色录像| freevideosex欧美| 大片电影免费在线观看免费| 黄色毛片三级朝国网站 | 男人爽女人下面视频在线观看| 中文字幕免费在线视频6| 男人狂女人下面高潮的视频| 精品一品国产午夜福利视频| 精品熟女少妇av免费看| 五月玫瑰六月丁香| 亚洲av免费高清在线观看| 人妻少妇偷人精品九色| 在线观看一区二区三区激情| 欧美日韩综合久久久久久| 日韩一本色道免费dvd| 国产精品三级大全| 国产精品一区www在线观看| 伊人久久精品亚洲午夜| 少妇的逼好多水| 久热久热在线精品观看| 一级,二级,三级黄色视频| 高清不卡的av网站| 97超碰精品成人国产| 亚洲av成人精品一二三区| 亚洲国产精品成人久久小说| 欧美日韩综合久久久久久| 搡老乐熟女国产| 国产日韩欧美亚洲二区| 精华霜和精华液先用哪个| 久久国内精品自在自线图片| 岛国毛片在线播放| a级毛片免费高清观看在线播放| 在线亚洲精品国产二区图片欧美 | 久久久久久伊人网av| 麻豆精品久久久久久蜜桃| 老司机影院毛片| 我要看日韩黄色一级片| 国产精品麻豆人妻色哟哟久久| 只有这里有精品99| 乱系列少妇在线播放| 国产白丝娇喘喷水9色精品| 97在线视频观看| 久久久久久人妻| 欧美精品人与动牲交sv欧美| 亚洲国产最新在线播放| 十八禁网站网址无遮挡 | 人人妻人人爽人人添夜夜欢视频 | 国产在线一区二区三区精| 亚洲av国产av综合av卡| 精品久久久精品久久久| 久久鲁丝午夜福利片| 精品久久国产蜜桃| 久久久国产欧美日韩av| 26uuu在线亚洲综合色| 曰老女人黄片| 美女脱内裤让男人舔精品视频| 9色porny在线观看| 国产av一区二区精品久久| 久热久热在线精品观看| 汤姆久久久久久久影院中文字幕| 下体分泌物呈黄色| 欧美xxⅹ黑人| 久久午夜福利片| 色5月婷婷丁香| 午夜激情久久久久久久| 伊人久久国产一区二区| 欧美日韩精品成人综合77777| 国产成人aa在线观看| 黑丝袜美女国产一区| 久久人人爽av亚洲精品天堂| 日韩一区二区视频免费看| www.色视频.com| 亚洲av综合色区一区| 久久人妻熟女aⅴ| 日韩成人av中文字幕在线观看| 亚洲av成人精品一区久久| 精华霜和精华液先用哪个| 日韩 亚洲 欧美在线| 久久女婷五月综合色啪小说| 欧美性感艳星| 亚洲国产毛片av蜜桃av| 边亲边吃奶的免费视频| 少妇丰满av| 国产免费福利视频在线观看| 大片电影免费在线观看免费| 五月天丁香电影| a级片在线免费高清观看视频| 欧美日韩精品成人综合77777| 日日爽夜夜爽网站| 美女视频免费永久观看网站| 99re6热这里在线精品视频| 一区在线观看完整版| 欧美日韩视频高清一区二区三区二| 免费看av在线观看网站| 国产黄片美女视频| 中国三级夫妇交换| 曰老女人黄片| 97超碰精品成人国产| 免费高清在线观看视频在线观看| 一级毛片 在线播放| 婷婷色综合www| 十八禁高潮呻吟视频 | 七月丁香在线播放| 亚洲欧美一区二区三区黑人 | 十八禁网站网址无遮挡 | 黄色配什么色好看| 我要看日韩黄色一级片| 免费看不卡的av| 免费不卡的大黄色大毛片视频在线观看| 久久久久久久国产电影| 亚洲国产日韩一区二区| 日韩人妻高清精品专区| 亚洲av成人精品一二三区| 王馨瑶露胸无遮挡在线观看| 日本91视频免费播放| 国产熟女欧美一区二区| 伦理电影免费视频| 乱系列少妇在线播放| 国产午夜精品一二区理论片| 欧美精品一区二区免费开放| 中国国产av一级| 建设人人有责人人尽责人人享有的| 中国美白少妇内射xxxbb| 日本午夜av视频| 一级毛片我不卡| a级一级毛片免费在线观看| 国产黄频视频在线观看| 精品少妇内射三级| 99热这里只有是精品50| 久久免费观看电影| 校园人妻丝袜中文字幕| 黑人巨大精品欧美一区二区蜜桃 | 美女视频免费永久观看网站| 国产有黄有色有爽视频| 国国产精品蜜臀av免费| 国产一级毛片在线| 一级毛片黄色毛片免费观看视频| 大片免费播放器 马上看| 人妻制服诱惑在线中文字幕| 天美传媒精品一区二区| 久久久久久久久久成人| 人人妻人人添人人爽欧美一区卜| 日日摸夜夜添夜夜添av毛片| 国产高清国产精品国产三级| 日韩成人伦理影院| 国产精品一区二区在线不卡| 久久毛片免费看一区二区三区| 日韩av免费高清视频| 在线观看人妻少妇| 日日撸夜夜添| 国产日韩欧美亚洲二区| 一本—道久久a久久精品蜜桃钙片| 99热这里只有精品一区| 26uuu在线亚洲综合色| 嘟嘟电影网在线观看| 久久久午夜欧美精品| 下体分泌物呈黄色| 久久免费观看电影| 久久人人爽人人片av| 欧美精品一区二区大全| 亚洲人与动物交配视频| 制服丝袜香蕉在线| 一本一本综合久久| 男女免费视频国产| 女性被躁到高潮视频| 成年美女黄网站色视频大全免费 | 国产成人精品一,二区| 亚洲综合色惰| 丝袜喷水一区| 午夜激情久久久久久久| 久久这里有精品视频免费| 一边亲一边摸免费视频| 国产成人精品一,二区| 亚洲精品色激情综合| 精品亚洲成a人片在线观看| 精品一区二区免费观看| 一个人免费看片子| 久久人人爽人人爽人人片va| 国产淫语在线视频| av一本久久久久| 久久婷婷青草| 亚洲精品国产色婷婷电影| 少妇被粗大的猛进出69影院 | 91午夜精品亚洲一区二区三区| 国产片特级美女逼逼视频| 两个人的视频大全免费| av女优亚洲男人天堂| 亚洲精品aⅴ在线观看| 国产精品嫩草影院av在线观看| 三级国产精品欧美在线观看| 精品人妻熟女av久视频| 伊人亚洲综合成人网| 男人爽女人下面视频在线观看| 纵有疾风起免费观看全集完整版| 97精品久久久久久久久久精品| 美女大奶头黄色视频| 日韩免费高清中文字幕av| 国产熟女欧美一区二区| 成年人午夜在线观看视频| 欧美精品一区二区大全| 欧美日韩av久久| 纵有疾风起免费观看全集完整版| 简卡轻食公司| 亚洲精品视频女| 韩国av在线不卡| freevideosex欧美| 另类精品久久| 国产一区二区在线观看日韩| 久久久久久久大尺度免费视频| 26uuu在线亚洲综合色| 韩国高清视频一区二区三区| 国产精品成人在线| 亚洲人成网站在线播| 五月伊人婷婷丁香| 精品亚洲乱码少妇综合久久| 丰满人妻一区二区三区视频av| 99热这里只有是精品50| 国产精品欧美亚洲77777| 国产伦精品一区二区三区四那| 国产亚洲午夜精品一区二区久久| 欧美高清成人免费视频www| 99九九线精品视频在线观看视频| 岛国毛片在线播放| 人妻人人澡人人爽人人| 日韩欧美 国产精品| 丝袜脚勾引网站| 国内少妇人妻偷人精品xxx网站| 少妇 在线观看| 国产精品偷伦视频观看了| 91精品国产九色| 国产欧美日韩综合在线一区二区 | 中文欧美无线码| 欧美 日韩 精品 国产| 日韩不卡一区二区三区视频在线| 国产日韩欧美亚洲二区| 国产精品免费大片| 国产成人一区二区在线| 毛片一级片免费看久久久久| 天堂中文最新版在线下载| 一本一本综合久久| 国产黄色免费在线视频| 亚洲av综合色区一区| 国产极品粉嫩免费观看在线 | 少妇被粗大的猛进出69影院 | 日本午夜av视频| 九九爱精品视频在线观看| 女人精品久久久久毛片| 久久狼人影院| av天堂久久9| 国产探花极品一区二区| 91精品伊人久久大香线蕉| 国产男女超爽视频在线观看| 久久久久网色| 久久久午夜欧美精品| 亚洲精品456在线播放app| 亚洲av成人精品一二三区| 久久久a久久爽久久v久久| 日本猛色少妇xxxxx猛交久久| 赤兔流量卡办理| 在线观看美女被高潮喷水网站| 女人久久www免费人成看片| 精品国产国语对白av| av.在线天堂| 成年人午夜在线观看视频| 交换朋友夫妻互换小说| 国产亚洲精品久久久com| 亚洲人成网站在线观看播放| 国产色婷婷99| 国产午夜精品久久久久久一区二区三区| 免费观看的影片在线观看| 大香蕉久久网| 免费观看性生交大片5| 激情五月婷婷亚洲| 中文乱码字字幕精品一区二区三区| 人妻制服诱惑在线中文字幕| 亚洲精品国产色婷婷电影| 精品视频人人做人人爽| 爱豆传媒免费全集在线观看| 老熟女久久久| 久久久久精品久久久久真实原创| 婷婷色av中文字幕| 亚洲av国产av综合av卡| 下体分泌物呈黄色| 亚洲欧美一区二区三区黑人 | 日韩三级伦理在线观看| 老司机亚洲免费影院| .国产精品久久| 在现免费观看毛片| 国产极品天堂在线| 黑人巨大精品欧美一区二区蜜桃 | 国内揄拍国产精品人妻在线| 成人综合一区亚洲| 人人妻人人添人人爽欧美一区卜| 美女脱内裤让男人舔精品视频| 欧美精品亚洲一区二区| 少妇人妻一区二区三区视频| 丁香六月天网| 97在线视频观看| 欧美日本中文国产一区发布| 大香蕉97超碰在线| 久久99蜜桃精品久久| 蜜桃久久精品国产亚洲av| 18禁裸乳无遮挡动漫免费视频| 少妇被粗大的猛进出69影院 | 人人妻人人澡人人看| 久久国产精品男人的天堂亚洲 | 日韩av免费高清视频| 伦精品一区二区三区| 亚洲国产最新在线播放| 亚洲精品视频女| 久久久久久久国产电影| 久久人人爽人人爽人人片va| 一个人看视频在线观看www免费| 美女内射精品一级片tv| 男女啪啪激烈高潮av片| 亚洲美女搞黄在线观看| 夫妻午夜视频| 人妻制服诱惑在线中文字幕| 女性被躁到高潮视频| 91aial.com中文字幕在线观看| 亚洲图色成人| 国产男女内射视频| 岛国毛片在线播放| 国产精品久久久久久精品电影小说| 如日韩欧美国产精品一区二区三区 | 综合色丁香网| 肉色欧美久久久久久久蜜桃| 性色av一级| 精品一区二区免费观看| 国产又色又爽无遮挡免| 精品国产国语对白av| 一级,二级,三级黄色视频| 中文字幕免费在线视频6| 国产亚洲欧美精品永久| 国产一区二区三区综合在线观看 | 六月丁香七月| 日韩av在线免费看完整版不卡| 99热这里只有精品一区| √禁漫天堂资源中文www| 国产精品一区二区三区四区免费观看| 日本欧美国产在线视频| 哪个播放器可以免费观看大片| 99热这里只有是精品50| 日本av手机在线免费观看| 国产av码专区亚洲av| 国产精品国产三级国产专区5o| 精品一区二区免费观看| kizo精华| 精品久久久精品久久久| 一区二区三区精品91| 久久久久久久久久久丰满| 男女国产视频网站| 三级经典国产精品| 国产老妇伦熟女老妇高清| 国产熟女午夜一区二区三区 | 免费观看在线日韩| 国产精品久久久久久精品古装| 我的女老师完整版在线观看| 国产成人一区二区在线| 亚洲四区av| 亚洲精品日韩av片在线观看| 美女中出高潮动态图| 美女xxoo啪啪120秒动态图| 国产有黄有色有爽视频| 妹子高潮喷水视频| 精品午夜福利在线看| 欧美激情国产日韩精品一区| 国产精品嫩草影院av在线观看| 街头女战士在线观看网站| 新久久久久国产一级毛片| 亚洲国产精品一区三区| 插逼视频在线观看| 国产亚洲一区二区精品| av又黄又爽大尺度在线免费看| 人妻一区二区av| av不卡在线播放| 国产精品无大码| 夜夜爽夜夜爽视频| 国产高清有码在线观看视频| 国产免费一级a男人的天堂| 国产欧美日韩一区二区三区在线 | 亚洲欧美日韩另类电影网站| 亚洲av男天堂| 内射极品少妇av片p| 国产片特级美女逼逼视频| 国内揄拍国产精品人妻在线| 国产精品国产av在线观看| 最后的刺客免费高清国语| 中文欧美无线码| 中文字幕久久专区| 毛片一级片免费看久久久久| 2022亚洲国产成人精品| 中文在线观看免费www的网站| 视频中文字幕在线观看| 日韩欧美 国产精品| av线在线观看网站| 一级二级三级毛片免费看| 制服丝袜香蕉在线| 成年av动漫网址| 丰满乱子伦码专区| 欧美日韩一区二区视频在线观看视频在线| 色视频www国产| 又粗又硬又长又爽又黄的视频| 国产在线视频一区二区| 欧美精品亚洲一区二区| 女人精品久久久久毛片| 国产在线视频一区二区| 51国产日韩欧美| 国产成人91sexporn| 亚洲久久久国产精品| 老司机影院成人| 少妇丰满av| 亚洲精品一二三| 五月天丁香电影| 九九久久精品国产亚洲av麻豆| 人妻一区二区av| av女优亚洲男人天堂| 少妇人妻久久综合中文| 欧美激情国产日韩精品一区| 人妻少妇偷人精品九色| av免费在线看不卡| 人妻制服诱惑在线中文字幕| 在线免费观看不下载黄p国产| 桃花免费在线播放| 成年人免费黄色播放视频 | av.在线天堂| 精品99又大又爽又粗少妇毛片| 激情五月婷婷亚洲| 亚洲第一区二区三区不卡| 蜜桃久久精品国产亚洲av| 偷拍熟女少妇极品色| 美女大奶头黄色视频|