WO2017181444A1 - Method and device for measuring diffusion coefficient of carbon dioxide in crude oil - Google Patents
Method and device for measuring diffusion coefficient of carbon dioxide in crude oil Download PDFInfo
- Publication number
- WO2017181444A1 WO2017181444A1 PCT/CN2016/080898 CN2016080898W WO2017181444A1 WO 2017181444 A1 WO2017181444 A1 WO 2017181444A1 CN 2016080898 W CN2016080898 W CN 2016080898W WO 2017181444 A1 WO2017181444 A1 WO 2017181444A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- carbon dioxide
- crude oil
- pressure
- closed container
- container
- Prior art date
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N13/00—Investigating surface or boundary effects, e.g. wetting power; Investigating diffusion effects; Analysing materials by determining surface, boundary, or diffusion effects
- G01N13/04—Investigating osmotic effects
Definitions
- the present invention relates to the field of oil and gas field development, and more particularly to a method and apparatus for determining the diffusion coefficient of carbon dioxide in crude oil.
- the direct method is: sampling the fluid at different times and different diffusion distances, and then analyzing the samples to obtain the concentration data of carbon dioxide, and then deducing the diffusion coefficient of carbon dioxide in the fluid.
- the sampling process can interfere with the flow field in the system, causing experimental errors.
- the indirect method does not require sampling, and the NMR (nuclear magnetic resonance) and PVT (pressure-volume-temperature) methods are the two most widely used methods.
- the NMR method can directly measure the concentration of carbon dioxide in a closed system by nuclear magnetic resonance spectroscopy, but it is expensive and costly.
- the PVT method is a relatively conventional method.
- the pressure exhaustion method measures the pressure change of the gas-liquid diffusion system, and then calculates the diffusion coefficient by the curve of the pressure versus time.
- the direct or indirect measurement cannot directly measure the diffusion coefficient of carbon dioxide gas.
- these measurement methods are unable to obtain accurate diffusion coefficients of carbon dioxide gas, thereby failing to accurately confirm the equilibrium time of carbon dioxide gas in crude oil, which adversely affects oil recovery.
- the present invention proposes a method for determining the diffusion coefficient of carbon dioxide in crude oil.
- the diffusion coefficient of carbon dioxide in crude oil can be accurately determined to accurately predict the time when carbon dioxide reaches equilibrium in crude oil.
- the present invention also proposes a device for determining the diffusion coefficient of carbon dioxide in crude oil.
- a method for determining a diffusion coefficient of carbon dioxide in crude oil comprising the steps of
- Step 1 Filling a closed container of volume V C with V 0 volume of crude oil, where V 0 ⁇ V C ;
- Step two making the closed container after step one constant at the experimental temperature
- Step 3 Filling the sealed container after the second step with carbon dioxide gas, and measuring the pressure P(t) of the carbon dioxide gas every t time from the start time t 0 until the equilibrium pressure is reached, wherein at the initial time t 0 , The pressure of the carbon dioxide gas is P 0 (t);
- Step 4 According to the expansion coefficient ⁇ of the crude oil in which carbon dioxide is dissolved and using the relationship between the diffusion time t and the pressure P(t), the diffusion coefficient D AB of carbon dioxide in the crude oil is obtained.
- the factor of the volume expansion of the crude oil after the carbon dioxide is diffused into the crude oil is considered, so that the diffusion coefficient of carbon dioxide in the crude oil can be accurately determined.
- step four according to
- Z g is a compression factor of carbon dioxide, and Z g is approximately equal to 0.81774 in a pressure range of 5-7 MPa, and R is a general gas constant, that is, 8.314 cm 3 ⁇ MPa ⁇ mol -1 ⁇ K -1 , and T is an experimental temperature.
- x eq is the equilibrium interface concentration of carbon dioxide at the oil and gas interface, determined by experimental conditions, and A is the cross-sectional area of the core, determined by experimental conditions.
- V d is the volume of crude oil without carbon dioxide.
- step three the time interval t of the pressure P(t) of the carbon dioxide gas is measured to be 1 minute. This uniform time interval helps to automatically control the measurement process through the software, facilitating the operator's operation.
- step two the test temperature is 97.53 °C.
- the initial pressure P 0 (t) of the carbon dioxide gas is 7 MPa.
- step one the pressure in the closed vessel is maintained at a vacuum of 0.1 MPa prior to charging the closed vessel with crude oil. This lower pressure minimizes the effect of residual gas in the closed vessel on the measurement results.
- a device for determining the diffusion coefficient of carbon dioxide in crude oil comprising a closed container for testing, upstream of a closed container for testing and a closed container for testing Carbon dioxide a closed container, a crude oil sealed container connected to the closed container for testing downstream of the closed container for testing, wherein the carbon dioxide sealed container and the closed container are disposed in the incubator to thermostat the sealed container for testing at the experimental temperature .
- a vacuum pump that evacuates the closed container for testing is also included.
- An advantage of the present invention over the prior art is that the method of the present invention takes into account the effect of volume expansion of the crude oil after diffusion of carbon dioxide into the crude oil, thereby enabling accurate determination of the diffusion coefficient of carbon dioxide in the crude oil, and thus in carbon dioxide Oil production can accurately predict the balance time, which is conducive to the oil production work.
- the method of the present invention uses a uniform time interval to measure pressure, which facilitates automatic control and facilitates operator operation.
- Figure 1 is a schematic illustration of an experimental setup in accordance with the present invention.
- FIG. 2 is a schematic illustration of the steps in accordance with the present invention.
- Figure 3 is a graph showing the relationship between the volume expansion ratio and the pressure of the carbon dioxide-dissolved crude oil obtained according to the present invention.
- Figure 4 is a graph showing the relationship between the pressure correction term and the square root of the diffusion time in accordance with the present invention.
- Figure 5 is a graph of the relationship between the pressure term and the square root of the diffusion time in accordance with the prior art.
- FIG. 1 shows schematically an experimental setup 10 according to the invention.
- the experimental apparatus 10 includes a sealed container 15 for testing, a carbon dioxide sealed container 14 that communicates with the closed container 15 upstream of the sealed container 15, and a crude oil sealed container 29 that communicates with the closed container 15 downstream of the sealed container 15.
- the carbon dioxide sealed container 14 and the closed container 15 are provided in the incubator 16 so that the closed container 15 can be thermostated at the experimental temperature.
- a temperature sensor 18 is provided on the oven 16 to ensure a constant test temperature.
- a pressure sensor 17 is provided on the pipe 19 that connects the carbon dioxide container 14 and the hermetic container 15 to measure the pressure of carbon dioxide in the pipe 19. Since the conduit 19 is always in communication with the closed vessel 15, the pressure sensor 17 actually indicates the pressure of the carbon dioxide-crude oil system within the closed vessel 15. In one embodiment, in order to be able to evacuate the closed container 15, a vacuum pump 21 is also provided. As shown in Fig.
- the vacuum pump 21 is connected to the pipe 19 through a pipe 25, a valve 26 is provided between the junction of the pipe 25 and the pipe 19, and the outlet of the carbon dioxide container 14, and a valve 27 is provided downstream of the hermetic container 15, so that After the valves 26 and 27 are closed, the vacuum pump 21 is actually only in communication with the hermetic container 15, so that the hermetic container 15 can be evacuated.
- the pressure in the closed container 15 is maintained at a vacuum degree before the sealed container 15 is filled with the crude oil. 0.1 MPa.
- the experimental apparatus 10 further includes a carbon dioxide supply unit such as a pump 22, a crude oil supply unit such as a pump 23, and a data processing unit 24, which are well known to those skilled in the art and will not be further described herein for the sake of brevity.
- a carbon dioxide supply unit such as a pump 22
- a crude oil supply unit such as a pump 23
- a data processing unit 24 which are well known to those skilled in the art and will not be further described herein for the sake of brevity.
- FIG. 2 shows the steps of the method according to the invention
- Step 31 filling the sealed container 15 with a predetermined volume of V 0 crude oil
- Step 32 Constantly confining the closed container 15 to the experimental temperature
- Step 33 charging the sealed container 15 with carbon dioxide gas, and measuring the pressure P(t) of the carbon dioxide gas every t time from the start time t 0 until the equilibrium pressure is reached;
- Step 34 treating the measured pressure P(t) according to the expansion coefficient ⁇ of the crude oil in which the carbon dioxide is dissolved;
- Step 35 Obtain a diffusion coefficient D AB of carbon dioxide in the crude oil.
- starting time is defined as the time at which the carbon dioxide inflation ends.
- the volume V c of the hermetic container 15 was 174 ml.
- the experimental apparatus 10 was connected, it was washed with petroleum ether, dried with hot air, and then checked for airtightness of the entire apparatus to ensure airtightness.
- the valve 26 and the valve 27 were closed, and the vacuum pump 21 was started to evacuate the closed container 15 for 16 hours.
- the crude oil had a density of 0.87 g/cm 3 and a viscosity of 20.1 mPa ⁇ s.
- the closed vessel 15 and the carbon dioxide sealed vessel 14 were thermostated to an experimental temperature of 97.53 °C.
- the valve 26 is opened, and the sealed container 15 is filled with carbon dioxide gas by the pump 22 until the pressure sensor 17 indicates 7 MPa, and the valve 26 is closed.
- the pressure P(t) indicated by the pressure sensor 17 is recorded every 1 minute, where t is the time point at which the pressure is measured until the equilibrium pressure is reached. The measurement results are shown in Table 1.
- the expansion coefficient ⁇ of the crude oil in which carbon dioxide is dissolved is measured.
- the volume V L of the liquid phase is measured at the bubble point pressure of the crude oil in which different amounts of carbon dioxide are dissolved.
- the measurement was carried out under the conditions of a molar ratio of carbon dioxide to crude oil of 0.412, 1.813, 3.057, 3.806, respectively.
- V L /V d ⁇ P+1
- the pressure P is plotted on the abscissa
- V L /V d is plotted on the ordinate in the plane direct coordinate system, as shown in Fig. 3. Fitting It follows that the expansion coefficient ⁇ is equal to 0.0089.
- the parameters used were as follows: Z g carbon dioxide compression factor. Although those skilled in the art with a known amount g is a change in pressure change Z, the inventors found that in a pressure range of 5-8MPa the Z g may be considered to be approximately constant, and the inventors used this art
- the method well known to the skilled person was measured under the experimental conditions, and the result was 0.81774;
- R was a general gas constant, 8.314 cm 3 ⁇ MPa ⁇ mol -1 ⁇ K -1 ;
- T was an experimental temperature of 97.53 ° C;
- the diffusion coefficient D AB 0.000772233 cm 2 /s obtained by the method of the invention predicts that the time for the carbon dioxide to reach equilibrium in the crude oil is 450 days, and the method according to the prior art is used.
- the diffusion coefficient D' AB 0.001033 predicts that the time for carbon dioxide to reach equilibrium in crude oil is 352 days.
- Equation 1 used in the present invention.
- the boundary conditions are:
- V L V 0 [1+ ⁇ P(t)], (Equation 4)
- V L is the volume of the oil after absorption of carbon dioxide.
- Equation 2-5 is derived in parallel to obtain Equation 1 used in the present invention:
Abstract
A method and a device for measuring the diffusion coefficient of carbon dioxide in crude oil. The method comprises: step 1, filling a closed container of VC with crude oil of V0, where V0<VC; step 2, keeping the closed container at an constant experimental temperature; step 3, filling the closed container with carbon dioxide gas, measuring the pressure P(t) of the carbon dioxide gas once at every time interval t from the starting time t0 until a balanced pressure is reached, the initial pressure of the carbon dioxide gas at the starting time t0 being P0(t); and step 4, obtaining the diffusion coefficient DAB of carbon dioxide in crude oil according to the expansion coefficient β of crude oil with carbon dioxide dissolved therein and by using the relation between the diffusion time t and the pressure P(t). The device comprises: a closed container (15) for testing, a closed carbon dioxide container (14) provided upstream of the closed container (15) for testing, and a closed crude oil container (29) provided downstream of the closed container for testing.
Description
本发明涉及油气田的开发领域,更具体地涉及一种测定二氧化碳在原油中的扩散系数的方法和装置。The present invention relates to the field of oil and gas field development, and more particularly to a method and apparatus for determining the diffusion coefficient of carbon dioxide in crude oil.
在油气田开发中,采用二氧化碳驱油来可以提高原油的采收率,因此二氧化碳采油技术得到了广泛应用。使用二氧化碳采油的步骤如下:首先向油层中注入一定量的二氧化碳气体,然后关井使二氧化碳在原油中的扩散达到平衡,最后打开井进行采油。为了能预测二氧化碳在原油中达到平衡的时间,需要得知二氧化碳在原油中的扩散系数。如果扩散系数不准确,不能准确确认平衡时间,给采油带来不利影响。In the development of oil and gas fields, the use of carbon dioxide flooding can improve the recovery of crude oil, so carbon dioxide production technology has been widely used. The steps of using carbon dioxide to produce oil are as follows: first, a certain amount of carbon dioxide gas is injected into the oil layer, and then the well is closed to balance the diffusion of carbon dioxide in the crude oil, and finally the well is opened for oil recovery. In order to predict the time when carbon dioxide reaches equilibrium in crude oil, it is necessary to know the diffusion coefficient of carbon dioxide in crude oil. If the diffusion coefficient is not accurate, the equilibrium time cannot be accurately confirmed, which has an adverse effect on oil recovery.
在现有技术中,测量二氧化碳在液体中的扩散系数测定方法可分为两类:直接法和间接法。其中,直接法为:在不同的时间和不同的扩散距离对流体采样,然后把这些样本进行分析,得到二氧化碳的浓度数据,再推导出二氧化碳在流体中的扩散系数。但是,采样过程会干扰系统中的流场,从而引起实验的误差。间接法不需要取样,其中的NMR(核磁共振)和PVT(压力-体积-温度)方法是两种最广泛应用的方法。NMR法可以通过核磁共振光谱直接测量出密闭体系中二氧化碳的浓度,但是昂贵价格、成本较高。PVT法是比较常规的方法,以压力衰竭法为例,其测量气液扩散体系的压力变化,然后通过压力随时间的变化曲线计算得到扩散系数。但是,无论直接测量还是间接测量,都不能直接测量出二氧化碳气体的扩散系数。此外,这些测量方法都不能得到精确的二氧化碳气体的扩散系数由此不能准确确认二氧化碳气体在原油中的平衡时间,给采油带来不利影响。In the prior art, methods for measuring the diffusion coefficient of carbon dioxide in a liquid can be classified into two types: direct method and indirect method. Among them, the direct method is: sampling the fluid at different times and different diffusion distances, and then analyzing the samples to obtain the concentration data of carbon dioxide, and then deducing the diffusion coefficient of carbon dioxide in the fluid. However, the sampling process can interfere with the flow field in the system, causing experimental errors. The indirect method does not require sampling, and the NMR (nuclear magnetic resonance) and PVT (pressure-volume-temperature) methods are the two most widely used methods. The NMR method can directly measure the concentration of carbon dioxide in a closed system by nuclear magnetic resonance spectroscopy, but it is expensive and costly. The PVT method is a relatively conventional method. Taking the pressure exhaustion method as an example, it measures the pressure change of the gas-liquid diffusion system, and then calculates the diffusion coefficient by the curve of the pressure versus time. However, the direct or indirect measurement cannot directly measure the diffusion coefficient of carbon dioxide gas. In addition, these measurement methods are unable to obtain accurate diffusion coefficients of carbon dioxide gas, thereby failing to accurately confirm the equilibrium time of carbon dioxide gas in crude oil, which adversely affects oil recovery.
发明内容Summary of the invention
针对现有技术中所存在的上述技术问题,本发明提出了一种测定二氧化碳在原油中的扩散系数的方法。通过本发明的方法能够精确测定出二氧化碳在原油中扩散系数以能够准确预测二氧化碳在原油中达到平衡的时间。本发明还提出了一种测定二氧化碳在原油中的扩散系数的装置。In view of the above technical problems existing in the prior art, the present invention proposes a method for determining the diffusion coefficient of carbon dioxide in crude oil. By the method of the present invention, the diffusion coefficient of carbon dioxide in crude oil can be accurately determined to accurately predict the time when carbon dioxide reaches equilibrium in crude oil. The present invention also proposes a device for determining the diffusion coefficient of carbon dioxide in crude oil.
根据本发明的第一方面,提出了一种测定二氧化碳在原油中的扩散系数的方法,包括以下步骤,According to a first aspect of the invention, a method for determining a diffusion coefficient of carbon dioxide in crude oil is provided, comprising the steps of
步骤一:向容积为VC的密闭容器内充入V0体积的原油,其中V0<VC;
Step 1: Filling a closed container of volume V C with V 0 volume of crude oil, where V 0 <V C ;
步骤二:使步骤一之后的密闭容器恒定地处于实验温度下;Step two: making the closed container after step one constant at the experimental temperature;
步骤三:向步骤二之后的密闭容器内充入二氧化碳气体,从起始时间t0开始每隔t时间测量一次二氧化碳气体的压力P(t),直到达到平衡压力,其中在初始时间t0,二氧化碳气体的压力为P0(t);Step 3: Filling the sealed container after the second step with carbon dioxide gas, and measuring the pressure P(t) of the carbon dioxide gas every t time from the start time t 0 until the equilibrium pressure is reached, wherein at the initial time t 0 , The pressure of the carbon dioxide gas is P 0 (t);
步骤四:根据溶解有二氧化碳的原油的膨胀系数β并利用扩散时间t和压力P(t)之间的关系来得到二氧化碳在原油中的扩散系数DAB。Step 4: According to the expansion coefficient β of the crude oil in which carbon dioxide is dissolved and using the relationship between the diffusion time t and the pressure P(t), the diffusion coefficient D AB of carbon dioxide in the crude oil is obtained.
通过本发明的方法,考虑了二氧化碳扩散到原油中后原油体积膨胀的因素,因此能够精确地测定二氧化碳在原油中的扩散系数。By the method of the present invention, the factor of the volume expansion of the crude oil after the carbon dioxide is diffused into the crude oil is considered, so that the diffusion coefficient of carbon dioxide in the crude oil can be accurately determined.
在一个实施例中,在步骤四中,根据In one embodiment, in step four, according to
得到从而获得二氧化碳在原油中的扩散系数DAB, get Thereby obtaining the diffusion coefficient D AB of carbon dioxide in the crude oil,
其中,Zg为二氧化碳的压缩因子,在5-7MPa的压力范围内Zg近似等于0.81774,R为通用气体常数,即8.314cm3·MPa·mol-1·K-1,T为实验温度,xeq为油气界面处二氧化碳的平衡界面浓度,由实验条件确定,A为岩心截面积,由实验条件确定。Wherein, Z g is a compression factor of carbon dioxide, and Z g is approximately equal to 0.81774 in a pressure range of 5-7 MPa, and R is a general gas constant, that is, 8.314 cm 3 ·MPa·mol -1 ·K -1 , and T is an experimental temperature. x eq is the equilibrium interface concentration of carbon dioxide at the oil and gas interface, determined by experimental conditions, and A is the cross-sectional area of the core, determined by experimental conditions.
在一个实施例中,在步骤四中,在溶有不同量二氧化碳的原油在泡点压力P下测定液相体积VL,并通过VL/Vd=βP+1得到原油的膨胀系数β,其中,Vd是不含二氧化碳的原油的体积。通过本方法,能够精确测出原油因溶解二氧化碳而出产生的体积膨胀,以得到二氧化碳在原油中精确的扩散系数DAB。In one embodiment, in step four, different amount of carbon dioxide dissolved in the crude liquid volume V L measured at the bubble point pressure P, the via V L / V d = βP + 1 to give crude expansion coefficient β, Where V d is the volume of crude oil without carbon dioxide. Through the method, the volume expansion of the crude oil due to the dissolution of carbon dioxide can be accurately measured to obtain the precise diffusion coefficient D AB of carbon dioxide in the crude oil.
在一个实施例中,在步骤三中,测量二氧化碳气体的压力P(t)的时间间隔t为1分钟。这种均匀的时间间隔有助于通过软件来自动控制测定进程,方便了操作者的操作。In one embodiment, in step three, the time interval t of the pressure P(t) of the carbon dioxide gas is measured to be 1 minute. This uniform time interval helps to automatically control the measurement process through the software, facilitating the operator's operation.
在一个实施例中,在步骤二中,试验温度为97.53℃。在步骤三中,二氧化碳气体的初始压力P0(t)为7MPa。这些参数与实际采油工作的条件相一致,由此测定的二氧化碳在原油中的扩散系数能准确反映实际情况。In one embodiment, in step two, the test temperature is 97.53 °C. In the third step, the initial pressure P 0 (t) of the carbon dioxide gas is 7 MPa. These parameters are consistent with the conditions of the actual oil recovery work, and the measured diffusion coefficient of carbon dioxide in crude oil can accurately reflect the actual situation.
在一个实施例中,在步骤一中,在向密闭容器内充入原油之前,将密闭容器内的压力保持在真空度为0.1MPa。这种较低的压力能最大程度地降低密闭容器内残余气体对测定结果的影响。In one embodiment, in step one, the pressure in the closed vessel is maintained at a vacuum of 0.1 MPa prior to charging the closed vessel with crude oil. This lower pressure minimizes the effect of residual gas in the closed vessel on the measurement results.
根据本发明的第二方面,提出了一种测定二氧化碳在原油中的扩散系数的装置,包括用于测试的密闭容器,在用于测试的密闭容器的上游与用于测试的密闭容器相连通的二氧化碳
密闭容器,在用于测试的密闭容器的下游与用于测试的密闭容器相连通的原油密闭容器,其中二氧化碳密闭容器和密闭容器设置在恒温箱内以将用于测试的密闭容器恒温在实验温度。According to a second aspect of the invention, there is provided a device for determining the diffusion coefficient of carbon dioxide in crude oil, comprising a closed container for testing, upstream of a closed container for testing and a closed container for testing Carbon dioxide
a closed container, a crude oil sealed container connected to the closed container for testing downstream of the closed container for testing, wherein the carbon dioxide sealed container and the closed container are disposed in the incubator to thermostat the sealed container for testing at the experimental temperature .
在一个实施例中,还包括给用于测试的密闭容器抽真空的真空泵。In one embodiment, a vacuum pump that evacuates the closed container for testing is also included.
在本申请中,相同的参数使用了相同的代号。In this application, the same parameters are used for the same parameters.
与现有技术相比,本发明的优点在于,本发明的方法考虑了二氧化碳扩散到原油中后原油体积膨胀的影响,由此能够精确地测定二氧化碳在原油中的扩散系数,并由此在二氧化碳采油中能较为准确地预测平衡时间,有利于采油工作的进行。本发明的方法采用均匀的时间间隔来测量压力,这有助于实现自动控制,方便了操作者的操作。An advantage of the present invention over the prior art is that the method of the present invention takes into account the effect of volume expansion of the crude oil after diffusion of carbon dioxide into the crude oil, thereby enabling accurate determination of the diffusion coefficient of carbon dioxide in the crude oil, and thus in carbon dioxide Oil production can accurately predict the balance time, which is conducive to the oil production work. The method of the present invention uses a uniform time interval to measure pressure, which facilitates automatic control and facilitates operator operation.
在下文中将基于实施例并参考附图来对本发明进行更详细的描述。其中:The invention will be described in more detail hereinafter based on the embodiments and with reference to the accompanying drawings. among them:
图1是根据本发明的实验装置的示意图;Figure 1 is a schematic illustration of an experimental setup in accordance with the present invention;
图2是根据本发明的步骤示意图;Figure 2 is a schematic illustration of the steps in accordance with the present invention;
图3是根据本发明得到的溶解有二氧化碳的原油的体积膨胀倍数与压力的关系;Figure 3 is a graph showing the relationship between the volume expansion ratio and the pressure of the carbon dioxide-dissolved crude oil obtained according to the present invention;
图4是根据本发明的压力修正项与扩散时间平方根的关系图;Figure 4 is a graph showing the relationship between the pressure correction term and the square root of the diffusion time in accordance with the present invention;
图5是根据本现有技术的的压力项与扩散时间平方根的关系图。Figure 5 is a graph of the relationship between the pressure term and the square root of the diffusion time in accordance with the prior art.
附图并未按照实际的比例绘制。The drawings are not drawn to scale.
下面将结合附图对本发明做进一步说明。The invention will now be further described with reference to the accompanying drawings.
图1示意性地显示了根据本发明的实验装置10。实验装置10包括用于测试的密闭容器15,在密闭容器15的上游与密闭容器15相连通的二氧化碳密闭容器14,在密闭容器15的下游与密闭容器15相连通的原油密闭容器29。其中,二氧化碳密闭容器14和密闭容器15设置在恒温箱16内以能够将密闭容器15恒温在实验温度。在一个实施例中,在恒温箱16上设置有温度传感器18,以保证恒定的试验温度。Figure 1 shows schematically an experimental setup 10 according to the invention. The experimental apparatus 10 includes a sealed container 15 for testing, a carbon dioxide sealed container 14 that communicates with the closed container 15 upstream of the sealed container 15, and a crude oil sealed container 29 that communicates with the closed container 15 downstream of the sealed container 15. Among them, the carbon dioxide sealed container 14 and the closed container 15 are provided in the incubator 16 so that the closed container 15 can be thermostated at the experimental temperature. In one embodiment, a temperature sensor 18 is provided on the oven 16 to ensure a constant test temperature.
在连通二氧化碳容器14和密闭容器15的管道19上设置有压力传感器17以测定管道19内的二氧化碳的压力。由于管道19始终保持与密闭容器15连通,因此压力传感器17实际指示了密闭容器15内的二氧化碳-原油系统的压力。在一个实施方案中,为了能够给密闭容器15抽真空,还设置有真空泵21。如图1所示,真空泵21通过管道25连接于管道19,在管道25与管道19的结合点和二氧化碳容器14的出口之间设置有阀门26,在密闭容器15的下游设置阀门27,这样在关闭阀门26和27之后,真空泵21实际上仅与密闭容器15连通,从而能将密闭容器15抽真空。在一个实施例中,为了地降低密闭容器15内残余气体对测定结果的影响,在向密闭容器15内充入原油之前,将密闭容器15内的压力保持在真空度为
0.1MPa。此外,实验装置10还包括有二氧化碳供给装置,例如泵22;原油供给装置,例如泵23;数据处理装置24,这些装置是本领域的技术人员所熟知的,为了简单起见,这里不再赘述。A pressure sensor 17 is provided on the pipe 19 that connects the carbon dioxide container 14 and the hermetic container 15 to measure the pressure of carbon dioxide in the pipe 19. Since the conduit 19 is always in communication with the closed vessel 15, the pressure sensor 17 actually indicates the pressure of the carbon dioxide-crude oil system within the closed vessel 15. In one embodiment, in order to be able to evacuate the closed container 15, a vacuum pump 21 is also provided. As shown in Fig. 1, the vacuum pump 21 is connected to the pipe 19 through a pipe 25, a valve 26 is provided between the junction of the pipe 25 and the pipe 19, and the outlet of the carbon dioxide container 14, and a valve 27 is provided downstream of the hermetic container 15, so that After the valves 26 and 27 are closed, the vacuum pump 21 is actually only in communication with the hermetic container 15, so that the hermetic container 15 can be evacuated. In one embodiment, in order to reduce the influence of the residual gas in the closed container 15 on the measurement result, the pressure in the closed container 15 is maintained at a vacuum degree before the sealed container 15 is filled with the crude oil.
0.1 MPa. In addition, the experimental apparatus 10 further includes a carbon dioxide supply unit such as a pump 22, a crude oil supply unit such as a pump 23, and a data processing unit 24, which are well known to those skilled in the art and will not be further described herein for the sake of brevity.
图2显示了根据本发明的方法的步骤,Figure 2 shows the steps of the method according to the invention,
步骤31:向密闭容器15内充入预定体积V0的原油;Step 31: filling the sealed container 15 with a predetermined volume of V 0 crude oil;
步骤32:将密闭容器15恒定到实验温度;Step 32: Constantly confining the closed container 15 to the experimental temperature;
步骤33:向密闭容器15内充入二氧化碳气体,从起始时间t0开始每隔t时间测量一次二氧化碳气体的压力P(t),直到达到平衡压力;Step 33: charging the sealed container 15 with carbon dioxide gas, and measuring the pressure P(t) of the carbon dioxide gas every t time from the start time t 0 until the equilibrium pressure is reached;
步骤34:根据溶有溶解有二氧化碳的原油的膨胀系数β处理所测得的压力P(t);Step 34: treating the measured pressure P(t) according to the expansion coefficient β of the crude oil in which the carbon dioxide is dissolved;
步骤35:得到二氧化碳在原油中的扩散系数DAB。Step 35: Obtain a diffusion coefficient D AB of carbon dioxide in the crude oil.
在本申请中,用语“起始时间”规定为二氧化碳充气结束之时的时刻。In the present application, the term "starting time" is defined as the time at which the carbon dioxide inflation ends.
为了能够精确确定膨胀系数β,在一个实施例中,在溶有不同量二氧化碳的原油在泡点压力P下测定液相体积VL,并通过VL/Vd=βP+1得到原油的膨胀系数β,其中,Vd是不含二氧化碳的原油的体积。In order to be able to accurately determine the expansion coefficient β, in one embodiment, the liquid phase volume V L is measured at a bubble point pressure P in a crude oil in which different amounts of carbon dioxide are dissolved, and the expansion of the crude oil is obtained by V L /V d =βP+1 Coefficient β, where V d is the volume of crude oil free of carbon dioxide.
下面根据图1来详细描述本发明。The invention will be described in detail below with reference to FIG.
实施例1:Example 1:
密闭容器15的容积Vc为174毫升。将实验装置10连接之后,用石油醚清洗,热空气吹干,然后检查整个装置的气密性,确保不漏气。关闭阀门26和阀门27,启动真空泵21给密闭容器15抽真空16小时。抽气结束后,打开阀门27同时通过泵23向密闭容器15内供入V0=100毫升原油,关闭阀门27。原油的密度为0.87g/cm3,粘度为20.1mPa·s。将密闭容器15和二氧化碳密闭容器14恒温到实验温度97.53℃。打开阀门26,通过泵22向密闭容器15内充入二氧化碳气体,直到压力传感器17指示为7MPa,关闭阀门26。每隔1分钟记录一次压力传感器17指示的压力P(t),这里t是测定压力的时间点,直到达到平衡压力。测定结果如表1所示。The volume V c of the hermetic container 15 was 174 ml. After the experimental apparatus 10 was connected, it was washed with petroleum ether, dried with hot air, and then checked for airtightness of the entire apparatus to ensure airtightness. The valve 26 and the valve 27 were closed, and the vacuum pump 21 was started to evacuate the closed container 15 for 16 hours. After the end of the pumping, the valve 27 is opened while V 0 = 100 ml of crude oil is supplied into the hermetic container 15 through the pump 23, and the valve 27 is closed. The crude oil had a density of 0.87 g/cm 3 and a viscosity of 20.1 mPa·s. The closed vessel 15 and the carbon dioxide sealed vessel 14 were thermostated to an experimental temperature of 97.53 °C. The valve 26 is opened, and the sealed container 15 is filled with carbon dioxide gas by the pump 22 until the pressure sensor 17 indicates 7 MPa, and the valve 26 is closed. The pressure P(t) indicated by the pressure sensor 17 is recorded every 1 minute, where t is the time point at which the pressure is measured until the equilibrium pressure is reached. The measurement results are shown in Table 1.
表1Table 1
时间/minTime/min | 压力/MPaPressure / MPa | 时间/minTime/min | 压力/MPaPressure / MPa | 时间/minTime/min | 压力/MPaPressure / MPa | 时间/minTime/min |
压力/MPaPressure / |
00 | 7.001457.00145 | 11 | 7.001457.00145 | 22 | 6.986856.98685 | 33 | 6.990486.99048 |
44 | 6.980656.98065 | 55 | 6.971876.97187 | 66 | 6.972256.97225 | 77 | 6.965666.96566 |
88 | 6.96366.9636 | 99 | 6.959216.95921 | 1010 | 6.951196.95119 | 1111 | 6.942546.94254 |
1212 | 6.946426.94642 | 1313 | 6.942296.94229 | 1414 | 6.936596.93659 | 1515 | 6.927566.92756 |
1616 | 6.935956.93595 | 1717 | 6.920976.92097 | 1818 | 6.912836.91283 | 1919 | 6.904946.90494 |
2020 | 6.90656.9065 | 21twenty one | 6.908446.90844 | 22twenty two | 6.892036.89203 | 23twenty three | 6.883636.88363 |
24twenty four | 6.875496.87549 | 2525 | 6.879376.87937 | 2626 | 6.86896.8689 | 2727 | 6.858316.85831 |
2828 | 6.848486.84848 | 2929 | 6.850166.85016 | 3030 | 6.846156.84615 | 3131 | 6.846156.84615 |
3232 | 6.831176.83117 | 3333 | 6.816456.81645 | 3434 | 6.808686.80868 | 3535 | 6.812826.81282 |
3636 | 6.812566.81256 | 3737 | 6.808566.80856 | 3838 | 6.79026.7902 | 3939 | 6.787886.78788 |
4040 | 6.775096.77509 | 4141 | 6.768896.76889 | 4242 | 6.766696.76669 | 4343 | 6.773146.77314 |
4444 | 6.770576.77057 | 4545 | 6.768896.76889 | 4646 | 6.760616.76061 | 4747 | 6.754036.75403 |
4848 | 6.746146.74614 | 4949 | 6.733866.73386 | 5050 | 6.716946.71694 | 5151 | 6.712936.71293 |
5252 | 6.702466.70246 | 5353 | 6.702086.70208 | 5454 | 6.702086.70208 | 5555 | 6.700276.70027 |
5656 | 6.70666.7066 | 5757 | 6.700146.70014 | 5858 | 6.702216.70221 | 5959 | 6.698076.69807 |
6060 | 6.702216.70221 | 6161 | 6.69826.6982 | 6262 | 6.691616.69161 | 6363 | 6.687736.68773 |
6464 | 6.689936.68993 | 6565 | 6.685796.68579 | 6666 | 6.681536.68153 | 6767 | 6.674696.67469 |
6868 | 6.67746.6774 | 6969 | 6.675456.67545 | 7070 | 6.673256.67325 | 7171 | 6.66466.6646 |
7272 | 6.654016.65401 | 7373 | 6.653636.65363 | 7474 | 6.650136.65013 | 7575 | 6.645616.64561 |
7676 | 6.639916.63991 | 7777 | 6.637596.63759 | 7878 | 6.625446.62544 | 7979 | 6.621186.62118 |
8080 | 6.617176.61717 | 8181 | 6.614596.61459 | 8282 | 6.604256.60425 | 8383 | 6.604136.60413 |
8484 | 6.6046.604 | 8585 | 6.595856.59585 | 8686 | 6.591596.59159 | 8787 | 6.58596.5859 |
8888 | 6.579186.57918 | 8989 | 6.577116.57711 | 9090 | 6.581516.58151 | 9191 | 6.585266.58526 |
9292 | 6.591856.59185 | 9393 | 6.589656.58965 | 9494 | 6.587456.58745 | 9595 | 6.585526.58552 |
9696 | 6.580616.58061 | 9797 | 6.575186.57518 | 9898 | 6.570666.57066 | 9999 | 6.568726.56872 |
100100 | 6.564586.56458 | 101101 | 6.562136.56213 | 102102 | 6.559816.55981 | 103103 | 6.557996.55799 |
104104 | 6.556056.55605 | 105105 | 6.553996.55399 | 106106 | 6.551926.55192 | 107107 | 6.547536.54753 |
108108 | 6.545466.54546 | 109109 | 6.545466.54546 | 110110 | 6.537446.53744 | 111111 | 6.534746.53474 |
112112 | 6.532546.53254 | 113113 | 6.534616.53461 | 114114 | 6.530096.53009 | 115115 | 6.526346.52634 |
116116 | 6.520146.52014 | 117117 | 6.520266.52026 | 118118 | 6.514196.51419 | 119119 | 6.513556.51355 |
120120 | 6.507986.50798 | 121121 | 6.507736.50773 | 122122 | 6.503346.50334 | 123123 | 6.501536.50153 |
124124 | 6.499586.49958 | 125125 | 6.492626.49262 | 126126 | 6.490686.49068 | 127127 | 6.48686.4868 |
128128 | 6.478526.47852 | 129129 | 6.482276.48227 | 130130 | 6.472196.47219 | 131131 | 6.469876.46987 |
132132 | 6.465736.46573 | 133133 | 6.461726.46172 | 134134 | 6.459666.45966 | 135135 | 6.451386.45138 |
136136 | 6.451386.45138 | 137137 | 6.449196.44919 | 138138 | 6.445056.44505 | 139139 | 6.43696.4369 |
140140 | 6.434586.43458 | 141141 | 6.430576.43057 | 142142 | 6.428516.42851 | 143143 | 6.428256.42825 |
144144 | 6.421926.42192 | 145145 | 6.417916.41791 | 146146 | 6.415976.41597 | 147147 | 6.415976.41597 |
148148 | 6.411846.41184 | 149149 | 6.407066.40706 | 150150 | 6.407576.40757 | 151151 | 6.404866.40486 |
152152 | 6.399056.39905 | 153153 | 6.397116.39711 | 154154 | 6.394786.39478 | 155155 | 6.392466.39246 |
156156 | 6.388326.38832 | 157157 | 6.38476.3847 | 158158 | 6.37636.3763 | 159159 | 6.380566.38056 |
160160 | 6.376816.37681 | 161161 | 6.368036.36803 | 162162 | 6.364026.36402 | 163163 | 6.361956.36195 |
164164 | 6.357436.35743 | 165165 | 6.353176.35317 | 166166 | 6.345286.34528 | 167167 | 6.341276.34127 |
168168 | 6.334556.33455 | 169169 | 6.3336.333 | 170170 | 6.326546.32654 | 171171 | 6.320466.32046 |
172172 | 6.318266.31826 | 173173 | 6.314136.31413 | 174174 | 6.305736.30573 | 175175 | 6.305736.30573 |
176176 | 6.299276.29927 | 177177 | 6.297846.29784 | 178178 | 6.293066.29306 | 179179 | 6.282726.28272 |
180180 | 6.282726.28272 | 181181 | 6.276396.27639 | 182182 | 6.272646.27264 | 183183 | 6.265676.26567 |
184184 | 6.262046.26204 | 185185 | 6.258176.25817 | 186186 | 6.254546.25454 | 187187 | 6.24776.2477 |
188188 | 6.241626.24162 | 189189 | 6.239426.23942 | 190190 | 6.23366.2336 | 191191 | 6.233356.23335 |
192192 | 6.220816.22081 | 193193 | 6.221076.22107 | 194194 | 6.214616.21461 | 195195 | 6.212156.21215 |
196196 | 6.208536.20853 | 197197 | 6.201816.20181 | 198198 | 6.19386.1938 | 199199 | 6.195996.19599 |
200200 | 6.187596.18759 | 201201 | 6.185276.18527 | 202202 | 6.181266.18126 | 203203 | 6.17486.1748 |
204204 | 6.172736.17273 | 205205 | 6.166656.16665 | 206206 | 6.164466.16446 | 207207 | 6.160076.16007 |
208208 | 6.158136.15813 | 209209 | 6.150366.15036 | 210210 | 6.147926.14792 | 211211 | 6.141586.14158 |
212212 | 6.13996.1399 | 213213 | 6.12936.1293 | 214214 | 6.127236.12723 | 215215 | 6.125296.12529 |
216216 | 6.121036.12103 | 217217 | 6.114956.11495 | 218218 | 6.112376.11237 | 219219 | 6.108376.10837 |
220220 | 6.108116.10811 | 221221 | 6.098026.09802 | 222222 | 6.099846.09984 | 223223 | 6.09356.0935 |
224224 | 6.09356.0935 | 225225 | 6.087046.08704 | 226226 | 6.08736.0873 | 227227 | 6.079166.07916 |
228228 | 6.072826.07282 | 229229 | 6.072956.07295 | 230230 | 6.067126.06712 | 231231 | 6.06216.0621 |
232232 | 6.060036.06003 | 233233 | 6.05866.0586 | 234234 | 6.053826.05382 | 235235 | 6.05156.0515 |
236236 | 6.048136.04813 | 237237 | 6.043746.04374 | 238238 | 6.037416.03741 | 239239 | 6.033526.03352 |
240240 | 6.025256.02525 | 241241 | 6.029136.02913 | 242242 | 6.022936.02293 | 243243 | 6.016726.01672 |
244244 | 6.012726.01272 | 245245 | 6.010786.01078 | 246246 | 6.010656.01065 | 247247 | 6.004316.00431 |
248248 | 5.999795.99979 | 249249 | 5.995785.99578 | 250250 | 5.991525.99152 | 251251 | 5.991655.99165 |
252252 | 5.98575.9857 | 253253 | 5.983255.98325 | 254254 | 5.977565.97756 | 255255 | 5.977685.97768 |
256256 | 5.969285.96928 | 257257 | 5.971485.97148 | 258258 | 5.966955.96695 | 259259 | 5.958685.95868 |
260260 | 5.954935.95493 | 261261 | 5.950415.95041 | 262262 | 5.950415.95041 | 263263 | 5.952225.95222 |
264264 | 5.942395.94239 | 265265 | 5.939945.93994 | 266266 | 5.933865.93386 | 267267 | 5.936065.93606 |
268268 | 5.929475.92947 | 269269 | 5.92965.9296 | 270270 | 5.927665.92766 | 271271 | 5.921335.92133 |
272272 | 5.915125.91512 | 273273 | 5.913055.91305 | 274274 | 5.911375.91137 | 275275 | 5.910735.91073 |
276276 | 5.905045.90504 | 277277 | 5.902975.90297 | 278278 | 5.900395.90039 | 279279 | 5.898195.89819 |
280280 | 5.89035.8903 | 281281 | 5.887725.88772 | 282282 | 5.886045.88604 | 283283 | 5.890175.89017 |
284284 | 5.879585.87958 | 285285 | 5.879965.87996 | 286286 | 5.873385.87338 | 287287 | 5.873385.87338 |
288288 | 5.865485.86548 | 289289 | 5.869245.86924 | 290290 | 5.862785.86278 | 291291 | 5.861225.86122 |
292292 | 5.856835.85683 | 293293 | 5.855015.85501 | 294294 | 5.854885.85488 | 295295 | 5.848555.84855 |
296296 | 5.844425.84442 | 297297 | 5.842475.84247 | 298298 | 5.840025.84002 | 299299 | 5.836145.83614 |
300300 | 5.838215.83821 | 301301 | 5.829555.82955 | 302302 | 5.829815.82981 | 303303 | 5.825295.82529 |
304304 | 5.823735.82373 | 305305 | 5.821545.82154 | 306306 | 5.819215.81921 | 307307 | 5.817025.81702 |
308308 | 5.813015.81301 | 309309 | 5.80685.8068 | 310310 | 5.807065.80706 | 311311 | 5.804865.80486 |
312312 | 5.798275.79827 | 313313 | 5.798145.79814 | 314314 | 5.798785.79878 | 315315 | 5.796085.79608 |
316316 | 5.792455.79245 | 317317 | 5.794525.79452 | 318318 | 5.792075.79207 | 319319 | 5.783795.78379 |
320320 | 5.784055.78405 | 321321 | 5.781985.78198 | 322322 | 5.779665.77966 | 323323 | 5.771265.77126 |
324324 | 5.771645.77164 | 325325 | 5.7715.771 | 326326 | 5.76885.7688 | 327327 | 5.765185.76518 |
328328 | 5.763245.76324 | 329329 | 5.763375.76337 | 330330 | 5.756915.75691 | 331331 | 5.756655.75665 |
332332 | 5.754585.75458 | 333333 | 5.750965.75096 | 334334 | 5.75075.7507 | 335335 | 5.74365.7436 |
336336 | 5.746695.74669 | 337337 | 5.739855.73985 | 338338 | 5.741925.74192 | 339339 | 5.74015.7401 |
340340 | 5.733775.73377 | 341341 | 5.733775.73377 | 342342 | 5.732085.73208 | 343343 | 5.725245.72524 |
344344 | 5.727055.72705 | 345345 | 5.72555.7255 | 346346 | 5.719295.71929 | 347347 | 5.721235.72123 |
348348 | 5.719295.71929 | 349349 | 5.714775.71477 | 350350 | 5.715285.71528 | 351351 | 5.712965.71296 |
352352 | 5.713095.71309 | 353353 | 5.708955.70895 | 354354 | 5.704685.70468 | 355355 | 5.70085.7008 |
356356 | 5.700165.70016 | 357357 | 5.700555.70055 | 358358 | 5.697975.69797 | 359359 | 5.696155.69615 |
360360 | 5.690085.69008 | 361361 | 5.694085.69408 | 362362 | 5.694215.69421 | 363363 | 5.689955.68995 |
364364 | 5.685945.68594 | 365365 | 5.681675.68167 | 366366 | 5.681425.68142 | 367367 | 5.679225.67922 |
368368 | 5.677155.67715 | 369369 | 5.67565.6756 | 370370 | 5.673275.67327 | 371371 | 5.671335.67133 |
372372 | 5.673145.67314 | 373373 | 5.669145.66914 | 374374 | 5.66725.6672 | 375375 | 5.663325.66332 |
376376 | 5.665255.66525 | 377377 | 5.660995.66099 | 378378 | 5.660615.66061 | 379379 | 5.659185.65918 |
380380 | 5.656475.65647 | 381381 | 5.652845.65284 | 382382 | 5.654015.65401 | 383383 | 5.650525.65052 |
384384 | 5.648715.64871 | 385385 | 5.646645.64664 | 386386 | 5.648715.64871 | 387387 | 5.646385.64638 |
388388 | 5.644195.64419 | 389389 | 5.642125.64212 | 390390 | 5.640315.64031 | 391391 | 5.636045.63604 |
392392 | 5.635915.63591 | 393393 | 5.633975.63397 | 394394 | 5.635915.63591 | 395395 | 5.627775.62777 |
396396 | 5.629715.62971 | 397397 | 5.629585.62958 | 398398 | 5.625575.62557 | 399399 | 5.625195.62519 |
400400 | 5.625445.62544 |
测量溶有二氧化碳的原油的膨胀系数β。首先,在溶有不同量的二氧化碳的原油的泡点压力下测定液相的体积VL。这种测量方式是本领域的技术人员所熟知的,为了简单起见,这里不再赘述。然后,根据VL/Vd=βP+1得到膨胀系数β,这里Vd是不含二氧化碳的原油的体积。如表2所示,分别在二氧化碳与原油的摩尔比为0.412、1.813、3.057、3.806的条件进行测量。The expansion coefficient β of the crude oil in which carbon dioxide is dissolved is measured. First, the volume V L of the liquid phase is measured at the bubble point pressure of the crude oil in which different amounts of carbon dioxide are dissolved. Such measurement methods are well known to those skilled in the art and will not be described again for the sake of brevity. Then, the expansion coefficient β is obtained from V L /V d =βP+ 1, where V d is the volume of the crude oil containing no carbon dioxide. As shown in Table 2, the measurement was carried out under the conditions of a molar ratio of carbon dioxide to crude oil of 0.412, 1.813, 3.057, 3.806, respectively.
表2Table 2
摩尔比0.412Molar ratio 0.412 | 摩尔比1.813Molar ratio 1.813 | 摩尔比3.057Moore ratio 3.057 | 摩尔比3.806Moore ratio 3.806 | |
压力P/MPaPressure P/MPa | 4.1694.169 | 11.79011.790 | 18.50718.507 | 24.99524.995 |
Vd/mLV d /mL | 97.92997.929 | 96.65796.657 | 95.85295.852 | 95.20195.201 |
VL/mLV L /mL | 102.283102.283 | 105.199105.199 | 112.961112.961 | 116.243116.243 |
根据VL/Vd=βP+1,以压力P为横坐标,以VL/Vd为纵坐标在平面直接坐标系中作图,如图3所示。拟合得到由此得出膨胀系数β等于0.0089。According to V L /V d =βP+ 1, the pressure P is plotted on the abscissa and V L /V d is plotted on the ordinate in the plane direct coordinate system, as shown in Fig. 3. Fitting It follows that the expansion coefficient β is equal to 0.0089.
根据上文测得的膨胀系数β,对表1中的数据进行处理。以扩散时间平方根为横坐标,以压力修正项为纵坐标,在平面直角坐标系内作图,如图4所示。拟合得到y=0.8457x-26.512,其中y代表压力修正项,而x代表扩散时间平方根,由此得出斜率k等于0.8457。根据(以下称公式
1)得到,并因此得出二氧化碳在原油中的扩散系数DAB等于0.000772233cm2/s。The data in Table 1 was processed according to the expansion coefficient β measured above. Square root of diffusion time For the abscissa, the pressure correction term For the ordinate, plot in the plane Cartesian coordinate system, as shown in Figure 4. The fit yields y = 0.8457 x - 26.512, where y represents the pressure correction term and x represents the square root of the diffusion time, resulting in a slope k equal to 0.8457. according to (hereinafter referred to as formula 1), Therefore, it is concluded that the diffusion coefficient D AB of carbon dioxide in crude oil is equal to 0.000772233 cm 2 /s.
在本实施例中,所使用的参数如下:Zg为二氧化碳的压缩因子。尽管本领域的技术人员已知Zg是随着压力变化而变化的量,但是发明人发现在5-8MPa的压力范围内可以将Zg认为是近似恒定的,并且发明人使用了本领域的技术人员所熟知的方法在本实验条件下进行了测定,结果为0.81774;R为通用气体常数,8.314cm3·MPa·mol-1·K-1;T为实验温度97.53℃;xeq平衡状态下密闭容器15内P(t)压力下单位体积内二氧化碳的物质的量,其中A为岩心截面积,在本实施例中为4.71cm2;t是测定压力的时间点;t0则为开始测试的时间,在本实施例中,t0=0。In the present embodiment, the parameters used were as follows: Z g carbon dioxide compression factor. Although those skilled in the art with a known amount g is a change in pressure change Z, the inventors found that in a pressure range of 5-8MPa the Z g may be considered to be approximately constant, and the inventors used this art The method well known to the skilled person was measured under the experimental conditions, and the result was 0.81774; R was a general gas constant, 8.314 cm 3 · MPa·mol -1 · K -1 ; T was an experimental temperature of 97.53 ° C; x eq equilibrium state The amount of carbon dioxide in a unit volume under the pressure of P(t) in the closed vessel 15 , wherein A is the core cross-sectional area, which is 4.71 cm 2 in this embodiment; t is the time point at which the pressure is measured; t 0 is the time at which the test is started, and in the present embodiment, t 0 =0.
对比例1:Comparative example 1:
实验装置和实验过程与实施例1相同,区别仅在于步骤34,使用现有技术中的方法处理表1中的数据。在现有技术中,所依据的公式为:The experimental apparatus and experimental procedure were the same as in Example 1, except that in step 34, the data in Table 1 was processed using the methods of the prior art. In the prior art, the formula on which it is based is:
其中,P(t)、P0(t)、t、t0、Zg、R、T、xeq、A的含义与实施例1相同,这里不再赘述。The meanings of P(t), P 0 (t), t, t 0 , Z g , R, T, x eq , and A are the same as those in Embodiment 1, and are not described herein again.
以扩散时间的平方根为横坐标,以压力差P(t)-P0(t)为纵坐标,在平面直角坐标系内作图,如图5所示。拟合得到y=0.0124x-0.3833,其中y代表压力项,而x代表扩散时间平方根,斜率k′等于0.0124。根据得出二氧化碳在原油中的扩散系数D′AB等于0.001033cm2/s。Square root of diffusion time For the abscissa, the pressure difference P(t)-P 0 (t) is plotted on the ordinate and plotted in a plane Cartesian coordinate system, as shown in Figure 5. The fit yields y = 0.0124x - 0.3833, where y represents the pressure term and x represents the square root of the diffusion time and the slope k' is equal to 0.0124. according to It is concluded that the diffusion coefficient D' AB of carbon dioxide in crude oil is equal to 0.001033 cm 2 /s.
在实际采油过程中,使用本发明的方法得出的扩散系数DAB=0.000772233cm2/s,预测出二氧化碳在原油中达到平衡的时间为450天,而使用根据现有技术中的方法得出的扩散系数D′AB=0.001033,预测出二氧化碳在原油中达到平衡的时间为352天。而经实际验证,二氧化碳在原油中达到平衡的时间为465天,可看出通过本发明的方法得出的扩散系数更为准确。In the actual oil recovery process, the diffusion coefficient D AB = 0.000772233 cm 2 /s obtained by the method of the invention predicts that the time for the carbon dioxide to reach equilibrium in the crude oil is 450 days, and the method according to the prior art is used. The diffusion coefficient D' AB = 0.001033 predicts that the time for carbon dioxide to reach equilibrium in crude oil is 352 days. However, it has been verified that the time for carbon dioxide to reach equilibrium in crude oil is 465 days, and it can be seen that the diffusion coefficient obtained by the method of the present invention is more accurate.
下面来描述本发明所使用的公式1的得到过程:
The process of obtaining Equation 1 used in the present invention will be described below:
根据菲克第二定律According to Fick's second law
其中,x:二氧化碳在原油中的浓度;DAB:二氧化碳在原油中的扩散系数;t:扩散时间,z:扩散距离。在实验刚开始时,没有二氧化碳进入原油中,故初始条件为:0≤z≤z0,t=0,x=0Where x: concentration of carbon dioxide in crude oil; D AB : diffusion coefficient of carbon dioxide in crude oil; t: diffusion time, z: diffusion distance. At the beginning of the experiment, no carbon dioxide entered the crude oil, so the initial conditions were: 0 ≤ z ≤ z 0 , t = 0, x = 0
边界条件为:The boundary conditions are:
z=z0,t>0,x=xeq,z=z 0 , t>0, x=x eq ,
其中,xeq平衡状态下单位体积内二氧化碳的物质的量。Wherein, the amount of carbon dioxide per unit volume in the x eq equilibrium state.
求解菲克第二定律,得到Solve Fick's second law and get
这里仅考虑二氧化碳在原油中的扩散距离,因此使用y代替z。Only the diffusion distance of carbon dioxide in crude oil is considered here, so y is used instead of z.
根据理想气体的状态方程PV=nRT,在气体扩散时,随着时间t流逝气相减少的物质的量为:根据菲克第二定律,经过气液界面进入液相的气相的物质的量为:
According to the state equation PV=nRT of the ideal gas, the amount of the substance reduced by the gas phase as the time t elapses during gas diffusion is: According to Fick's second law, the amount of material in the gas phase that enters the liquid phase through the gas-liquid interface is:
根据气相减少的物质的量等于通过气液界面的进入液相的气相的物质的量这一公知原理,得出:According to the well-known principle that the amount of the substance reduced in the gas phase is equal to the amount of the substance in the gas phase entering the liquid phase through the gas-liquid interface, it is concluded that:
这里V是气相减少的体积。因为二氧化碳扩散进入原油之后会引起液相体积的膨胀,同时原油中的轻质组分也会被二氧化碳抽提。根据膨胀系数β的公式VL=Vd[1+βP(t)],由于原油的原始体积为V0,因此可将据膨胀系数β的公式改写为:Here V is the volume of the gas phase reduction. Because carbon dioxide diffuses into the crude oil, it causes expansion of the liquid phase volume, while the light components in the crude oil are also extracted by carbon dioxide. According to the formula V L =V d [1+βP(t)] of the expansion coefficient β, since the original volume of the crude oil is V 0 , the formula according to the expansion coefficient β can be rewritten as:
VL=V0[1+βP(t)],(公式4)V L =V 0 [1+βP(t)], (Equation 4)
这里VL是吸收了二氧化碳之后的原油的体积。在容器体积为VC时,Where V L is the volume of the oil after absorption of carbon dioxide. When the container volume is V C ,
V=VC-VL=VC-V0-βV0P(t),(公式5)
V=V C -V L =V C -V 0 -βV 0 P(t), (Equation 5)
将公式2-5联立推导,得到本发明所使用的公式1:Equation 2-5 is derived in parallel to obtain Equation 1 used in the present invention:
应注意地是,膨胀系数β也可以进行单独测定。It should be noted that the expansion coefficient β can also be measured separately.
虽然已经参考优选实施例对本发明进行了描述,但在不脱离本发明的范围的情况下,可以对其进行各种改进并且可以用等效物替换其中的部件。尤其是,只要不存在结构冲突,各个实施例中所提到的各项技术特征均可以任意方式组合起来。本发明并不局限于文中公开的特定实施例,而是包括落入权利要求的范围内的所有技术方案。
Although the present invention has been described with reference to the preferred embodiments thereof, various modifications may be made without departing from the scope of the invention. In particular, the technical features mentioned in the various embodiments can be combined in any manner as long as there is no structural conflict. The present invention is not limited to the specific embodiments disclosed herein, but includes all technical solutions falling within the scope of the claims.
Claims (3)
- 一种测定二氧化碳在原油中的扩散系数的方法,包括以下步骤,A method for determining a diffusion coefficient of carbon dioxide in crude oil, comprising the following steps,步骤一:向容积为VC的密闭容器内充入V0体积的原油,其中V0<VC;Step 1: Filling a closed container of volume V C with V 0 volume of crude oil, where V 0 <V C ;步骤二:使步骤一之后的密闭容器恒定地处于实验温度下;Step two: making the closed container after step one constant at the experimental temperature;步骤三:向步骤二之后的密闭容器内充入二氧化碳气体,从起始时间t0开始每隔t时间测量一次二氧化碳气体的压力P(t),直到达到平衡压力,其中在初始时间t0,二氧化碳气体的初始压力为P0(t);Step 3: Filling the sealed container after the second step with carbon dioxide gas, and measuring the pressure P(t) of the carbon dioxide gas every t time from the start time t 0 until the equilibrium pressure is reached, wherein at the initial time t 0 , The initial pressure of carbon dioxide gas is P 0 (t);步骤四:根据溶解有二氧化碳的原油的膨胀系数β并利用扩散时间t和压力P(t)之间的关系来得到二氧化碳在原油中的扩散系数DAB,在所述步骤四中,根据得到从而获得二氧化碳在原油中的扩散系数DAB,Step 4: According to the expansion coefficient β of the crude oil in which carbon dioxide is dissolved and using the relationship between the diffusion time t and the pressure P(t), the diffusion coefficient D AB of the carbon dioxide in the crude oil is obtained, in the step 4, according to get Thereby obtaining the diffusion coefficient D AB of carbon dioxide in the crude oil,其中,Zg为二氧化碳的压缩因子,在5-7MPa的压力范围内Zg近似等于0.81774,R为通用气体常数,即8.314cm3·MPa·mol-1·K-1,T为实验温度,xeq为油气界面处二氧化碳的平衡界面浓度,由实验条件确定,A为岩心截面积,由实验条件确定,在所述步骤四中,在溶有不同量二氧化碳的原油在泡点压力P下测定液相体积VL,并通过VL/Vd=βP+1得到原油的膨胀系数β,其中,Vd是不含二氧化碳的原油的体积,在所述步骤三中,测量二氧化碳气体的压力P(t)的时间间隔t为1分钟,,在所述步骤二中,实验温度为97.53℃,在所述步骤三中,所述二氧化碳气体的初始压力P0(t)为7MPa,在所述步骤一中,在向所述密闭容器内充入原油之前,将所述密闭容器内的压力保持在真空度为0.1MPa。Wherein, Z g is a compression factor of carbon dioxide, and Z g is approximately equal to 0.81774 in a pressure range of 5-7 MPa, and R is a general gas constant, that is, 8.314 cm 3 ·MPa·mol -1 ·K -1 , and T is an experimental temperature. x eq is the equilibrium interface concentration of carbon dioxide at the oil and gas interface, determined by the experimental conditions, A is the core cross-sectional area, determined by the experimental conditions. In the fourth step, the crude oil with different amounts of carbon dioxide is determined at the bubble point pressure P. The liquid phase volume V L , and the expansion coefficient β of the crude oil is obtained by V L /V d =βP+1, wherein V d is the volume of the crude oil without carbon dioxide, and in the third step, the pressure P of the carbon dioxide gas is measured. The time interval t of (t) is 1 minute, and in the second step, the experimental temperature is 97.53 ° C. In the third step, the initial pressure P 0 (t) of the carbon dioxide gas is 7 MPa, In the first step, the pressure in the sealed container is maintained at a degree of vacuum of 0.1 MPa before the crude oil is filled into the sealed container.
- 使用权利要求1方法的一种测定二氧化碳在原油中的扩散系数的装置,包括,用于测试的密闭容器,在所述用于测试的密闭容器的上游与所述用于测试的密闭容器相连通的二氧化碳密闭容器,在用于测试的密闭容器的下游与用于测试的密闭容器相连通的原油密闭容器,其中所述二氧化碳密闭容器和所述密闭容器设置在恒温箱内以将所述用于测试的密闭容器恒温在实验温度。A device for determining the diffusion coefficient of carbon dioxide in crude oil using the method of claim 1, comprising a closed container for testing, communicating with said closed container for testing upstream of said closed container for testing a carbon dioxide sealed container, a crude oil sealed container in communication with the closed container for testing downstream of the closed container for testing, wherein the carbon dioxide sealed container and the closed container are disposed in an incubator to be used for The sealed container tested was thermostated at the experimental temperature.
- 根据权利要求2所述的装置,其特征在于,还包括给所述用于测试的密闭容器抽真空的真空泵。 The apparatus of claim 2 further comprising a vacuum pump for evacuating said sealed container for testing.
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNPCT/CN2016/079900 | 2016-04-21 | ||
PCT/CN2016/079900 WO2017181392A1 (en) | 2016-04-21 | 2016-04-21 | Method for measuring diffusion coefficient of carbon dioxide in crude oil |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2017181444A1 true WO2017181444A1 (en) | 2017-10-26 |
Family
ID=60115461
Family Applications (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/CN2016/079900 WO2017181392A1 (en) | 2016-04-21 | 2016-04-21 | Method for measuring diffusion coefficient of carbon dioxide in crude oil |
PCT/CN2016/080898 WO2017181444A1 (en) | 2016-04-21 | 2016-05-03 | Method and device for measuring diffusion coefficient of carbon dioxide in crude oil |
Family Applications Before (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/CN2016/079900 WO2017181392A1 (en) | 2016-04-21 | 2016-04-21 | Method for measuring diffusion coefficient of carbon dioxide in crude oil |
Country Status (1)
Country | Link |
---|---|
WO (2) | WO2017181392A1 (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111024555A (en) * | 2018-10-10 | 2020-04-17 | 中国石油天然气股份有限公司 | Equipment and method for measuring fault characteristic gas diffusion speed of transformer |
CN113791101A (en) * | 2021-09-01 | 2021-12-14 | 大连理工大学 | Method for in-situ measurement of gas-liquid diffusion parameters in porous medium by nuclear magnetic resonance |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111855472B (en) * | 2019-04-26 | 2023-04-07 | 中国石油化工股份有限公司 | Evaluation device and method for extracting crude oil by carbon dioxide |
CN113295827B (en) * | 2021-05-24 | 2023-07-11 | 优百诺(成都)生物科技有限公司 | Detection system for oxygen release of oxygen-containing dressing |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2007205797A (en) * | 2006-01-31 | 2007-08-16 | Nissan Motor Co Ltd | Oxygen diffusion coefficient measurement method for porous body and its measurement instrument |
CN102809526A (en) * | 2012-08-27 | 2012-12-05 | 中国石油大学(华东) | Method for measuring diffusion coefficient of carbon dioxide in saturated oil core |
CN104502236A (en) * | 2014-12-15 | 2015-04-08 | 中国石油大学(华东) | Method for measuring diffusion coefficient and equilibrium concentration of CO2 in process of diffusion from water phase to oil phase |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6393897B1 (en) * | 2000-09-22 | 2002-05-28 | The Goodyear Tire & Rubber Company | Accelerated leakage testing |
FR2897434B1 (en) * | 2006-02-15 | 2014-07-11 | Commissariat Energie Atomique | METHOD AND DEVICE FOR PERMEATION MEASUREMENT |
CN101294949A (en) * | 2008-06-05 | 2008-10-29 | 江苏工业学院 | Evaluation system and method for oil product evaporating loss |
CN102314556A (en) * | 2011-08-29 | 2012-01-11 | 常州大学 | Method for evaluating oil gas recovery and system |
CN105300849B (en) * | 2015-11-27 | 2017-11-07 | 西南石油大学 | The test device and method of gas diffusivity in a kind of porous media |
-
2016
- 2016-04-21 WO PCT/CN2016/079900 patent/WO2017181392A1/en active Application Filing
- 2016-05-03 WO PCT/CN2016/080898 patent/WO2017181444A1/en active Application Filing
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2007205797A (en) * | 2006-01-31 | 2007-08-16 | Nissan Motor Co Ltd | Oxygen diffusion coefficient measurement method for porous body and its measurement instrument |
CN102809526A (en) * | 2012-08-27 | 2012-12-05 | 中国石油大学(华东) | Method for measuring diffusion coefficient of carbon dioxide in saturated oil core |
CN104502236A (en) * | 2014-12-15 | 2015-04-08 | 中国石油大学(华东) | Method for measuring diffusion coefficient and equilibrium concentration of CO2 in process of diffusion from water phase to oil phase |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111024555A (en) * | 2018-10-10 | 2020-04-17 | 中国石油天然气股份有限公司 | Equipment and method for measuring fault characteristic gas diffusion speed of transformer |
CN111024555B (en) * | 2018-10-10 | 2022-07-05 | 中国石油天然气股份有限公司 | Equipment and method for measuring fault characteristic gas diffusion speed of transformer |
CN113791101A (en) * | 2021-09-01 | 2021-12-14 | 大连理工大学 | Method for in-situ measurement of gas-liquid diffusion parameters in porous medium by nuclear magnetic resonance |
Also Published As
Publication number | Publication date |
---|---|
WO2017181392A1 (en) | 2017-10-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
WO2017181444A1 (en) | Method and device for measuring diffusion coefficient of carbon dioxide in crude oil | |
WO2017128479A1 (en) | Fully-automated system for testing gas permeability of rock and estimation method | |
CN108801860B (en) | Method for measuring diffusion coefficient of carbon dioxide in crude oil | |
US6595036B1 (en) | Method and apparatus for measuring amount of gas adsorption | |
WO2018236637A1 (en) | Method for correcting low permeability laboratory measurements for leaks | |
CN107655805B (en) | System and method for measuring permeability of low-permeability rock ore particles | |
CN109029619B (en) | Volume measuring device based on dynamic differential pressure attenuation | |
JPS63502450A (en) | Method and apparatus for measuring pore volume and permeability of close-core plugs | |
CN114412429B (en) | Method for testing relationship between crack size and nuclear magnetic T2 of Brazilian splitting method | |
JP2008209220A (en) | Leak inspection method and device | |
CN112505085B (en) | Method for measuring porosity effective stress coefficient based on nuclear magnetic resonance | |
CN105784054A (en) | System and method suitable for calibrating volume of gas cavity in irregular shape | |
JP4994494B2 (en) | Differential pressure measurement method and apparatus | |
CN109030783A (en) | A kind of rock-fluid reciprocation experimental system for simulating and method | |
CN110187011B (en) | Simulation experiment device and simulation method for dynamic adsorption and desorption research | |
CN112611675A (en) | Reaction kettle, device and method for testing solubility of natural gas under high-temperature and high-pressure conditions | |
CN113916748B (en) | Device and method for measuring shale matrix permeability and recovery ratio by light oil | |
CN208636318U (en) | A kind of rock-fluid reciprocation experimental system for simulating | |
CN110927359A (en) | Experimental test device and method for gas loss content in low-permeability porous medium coring process | |
CN107991215B (en) | Method for testing natural aperture and specific surface area of large-size low-permeability rock sample | |
CN110873674B (en) | Device for measuring saturated vapor pressure of solid matter | |
CN209783837U (en) | multi-medium leakage rate adjustable leakage rate calibration system | |
CN110455694B (en) | Rock sample permeability determination method | |
JP2017106738A (en) | Method of testing air tightness of tanks | |
CN111781228B (en) | Test device for measuring boiling point of viscous substance and application thereof |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
NENP | Non-entry into the national phase |
Ref country code: DE |
|
121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 16899031 Country of ref document: EP Kind code of ref document: A1 |
|
122 | Ep: pct application non-entry in european phase |
Ref document number: 16899031 Country of ref document: EP Kind code of ref document: A1 |