CN111081320A - Method for determining high-pressure physical property parameters of thickened oil-methane-carbon dioxide-propane system - Google Patents

Method for determining high-pressure physical property parameters of thickened oil-methane-carbon dioxide-propane system Download PDF

Info

Publication number
CN111081320A
CN111081320A CN201911228850.6A CN201911228850A CN111081320A CN 111081320 A CN111081320 A CN 111081320A CN 201911228850 A CN201911228850 A CN 201911228850A CN 111081320 A CN111081320 A CN 111081320A
Authority
CN
China
Prior art keywords
propane
methane
oil
formula
coefficient
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201911228850.6A
Other languages
Chinese (zh)
Other versions
CN111081320B (en
Inventor
孙晓飞
蔡林峰
宋兆尧
张艳玉
尉雪梅
李朋
李晓宇
龚航飞
王翔宇
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China University of Petroleum UPC East China
Original Assignee
China University of Petroleum UPC East China
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China University of Petroleum UPC East China filed Critical China University of Petroleum UPC East China
Priority to CN201911228850.6A priority Critical patent/CN111081320B/en
Publication of CN111081320A publication Critical patent/CN111081320A/en
Application granted granted Critical
Publication of CN111081320B publication Critical patent/CN111081320B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16CCOMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
    • G16C20/00Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
    • G16C20/10Analysis or design of chemical reactions, syntheses or processes
    • GPHYSICS
    • G06COMPUTING; CALCULATING; COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/18Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16CCOMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
    • G16C20/00Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
    • G16C20/30Prediction of properties of chemical compounds, compositions or mixtures
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16CCOMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
    • G16C60/00Computational materials science, i.e. ICT specially adapted for investigating the physical or chemical properties of materials or phenomena associated with their design, synthesis, processing, characterisation or utilisation

Abstract

The invention relates to a method for determining high-pressure physical property parameters of a thickened oil-methane-carbon dioxide-propane system, which is based on a PR state equation, adopts an optimized α function and a binary interaction coefficient suitable for the thickened oil-methane-carbon dioxide-propane system, and combines an optimized volume correction method to establish a high-pressure physical property prediction method.

Description

Method for determining high-pressure physical property parameters of thickened oil-methane-carbon dioxide-propane system
Technical Field
The invention relates to a method for determining high-pressure physical property parameters of a heavy oil-methane-carbon dioxide-propane system, belonging to the technical field of heavy oil reservoir development.
Background
At present, the reserves of conventional petroleum resources in China are continuously reduced, the external dependence of crude oil is continuously high, and the energy contradiction is very sharp. Under a new strategic situation, the center of gravity of oil exploration and development is shifted to complex oil reservoirs represented by heavy oil reservoirs and the like. And the thick oil reserves in China are rich, and the thick oil reserves are widely distributed in more than ten oil fields such as victory, Changqing, Xinjiang, Liaohe and the like, so that the thick oil has huge development potential.
China has wide distribution of thick oil resources and abundant reserves, but often has complex geological characteristics such as thin layers, edge and bottom water, low permeability and the like. The conventional heavy oil thermal recovery process has the disadvantages of serious heat loss, low efficiency and high cost. The thickened oil gas injection technology is a new technology developed in recent years, and injected gas can reduce the viscosity of thickened oil, improve the flowability of the thickened oil, promote the volume expansion of the thickened oil, delay the reduction of the formation pressure and has good application prospect.
After the thick oil is injected, high-pressure physical properties such as bubble point pressure, expansion coefficient and the like can be changed, and the high-pressure physical properties of the thick oil have important significance for researching the thick oil exploitation mechanism, optimizing injection and production parameters, designing a shaft and a ground oil pipeline and the like, so that the determination of the high-pressure physical properties of the thick oil is important for the thick oil injection technology. Heavy oil gas injection techniqueThe gas is a single gas or a mixed gas, and the common gas comprises CO2、CH4、C3H8And the like. Propane can effectively reduce viscosity of thickened oil, natural gas (CH)4、CO2Mainly) can inhibit propane liquefaction and is cheap and easy to obtain, so that the injection of methane-carbon dioxide-propane mixed gas comprehensively utilizes C3H8、CO2、CH4The advantages of the three gases are better than the advantages of injecting a single gas.
At present, with the development of indoor experiments and corresponding theoretical researches, experts and scholars at home and abroad obtain a plurality of achievements on the research on the phase behavior of thickened oil gas injection. However, some problems still exist, which are summarized as follows:
① lack the study on the change rules of high-pressure physical parameters such as bubble point pressure, expansion coefficient and the like of a thickened oil-methane-carbon dioxide-propane system under different gas injection quantities, ② most of the current thickened oil gas injection behavior studies still use single gas, the study on the thickened oil gas injection mixed gas is still in the development stage, the study on the thickened oil-mixed gas system phase behavior formed by the thickened oil gas injection of the methane-carbon dioxide-propane mixed gas is not related, ③ because the study object is thickened oil, the mixed gas is difficult to dissolve, and stirring is needed when a sample is prepared, so that the time consumption of the thickened oil gas injection high-pressure physical property experiment is very long, and the requirement that field workers can quickly master the high-pressure physical property parameters of the thickened oil-methane-carbon dioxide-propane system is difficult to meet only by indoor experiments.
Due to the defects, the high-pressure physical property parameters of the conventional thickened oil-methane-carbon dioxide-propane system are difficult to accurately calculate, the requirements of a thickened oil methane-carbon dioxide-propane mixed gas injection development technology cannot be met, and the production efficiency is reduced. Therefore, the invention forms a method for determining the high-pressure physical property parameters of the thickened oil-methane-carbon dioxide-propane system by deeply researching the phase behavior of the thickened oil after injecting the methane-carbon dioxide-propane mixed gas, realizes the rapid and accurate prediction of the high-pressure physical property parameters of the system, and lays a foundation for improving the recovery efficiency of the thickened oil reservoir injected with the mixed gas.
Disclosure of Invention
Aiming at the defects of the prior art, the invention provides a method for determining high-pressure physical property parameters of a thickened oil-methane-carbon dioxide-propane system.
The specific technical scheme of the invention is as follows:
a method for determining high-pressure physical property parameters of a thickened oil-methane-carbon dioxide-propane system comprises the following steps:
step 1, determining the oil reservoir temperature T, the heavy oil specific gravity gamma and the relative molecular mass M of each component of the systemi(i ═ 1,2,3,4, respectively representing propane, carbon dioxide, methane, heavy oil), crude oil volume coefficient BoAnd the mole fraction x of the methane-carbon dioxide-propane mixed gas in the system is givengMolar fraction N of each gas component in the gas mixturei(i is 1,2,3), and C in the system is calculated according to the formulas (1) to (4)3H8、CO2、CH4And the mole fraction x of the heavy oil componenti
x1=N1xg(1)
x2=N2xg(2)
x3=N3xg(3)
x4=1-xg(4)
Step 2, determining critical parameters including critical temperature T of each gas component by looking up a common gas main physical and chemical property tableciCritical pressure PciCritical molar volume Vci(i is 1,2,3), and calculating the critical parameters of the thickened oil according to the formulas (5) to (7);
in the formula, TbRIs the standard boiling point of the thickened oil, K; t isc4Is the critical temperature of the thick oil, K; pc4Critical pressure of thick oil, kPa;
step 3, determining the attraction coefficient a of each component in the thickened oil-methane-carbon dioxide-propane system according to the formulas (8) to (10)iAnd van der Waals molar volume bi
ai=aciα(Trii)(i=1,2,3,4) (8)
Wherein R is a general gas constant of 8.314kPa · m3/(K·kmol);TriThe conversion temperature of each component in the thickened oil-methane-carbon dioxide-propane system is dimensionless; omegaiThe eccentricity factors of all components in a thickened oil-methane-carbon dioxide-propane system can be obtained by looking up a pure material thermodynamic basic parameter table; pciCritical pressure, kPa, of each component in a thickened oil-methane-carbon dioxide-propane system; t isciIs the critical temperature of each component in the thickened oil-methane-carbon dioxide-propane system, K, α is a dimensionless parameter;
step 4, determining the system binary interaction coefficient delta according to the formulas (11) to (14)ijWherein δ12Is C3H8And CO2Binary coefficient of interaction, delta13Is C3H8And CH4Binary coefficient of interaction, delta23Is CO2And CH4Binary coefficient of interaction, delta14Is C3H8Coefficient of interaction, delta, with heavy oil24Is CO2Coefficient of interaction, delta, with heavy oil34Is CH4The coefficient of interaction with the thickened oil is binary;
δ14=-0.4560T/Tc4+0.1817 (12)
δ24=-0.5462T/Tc4-0.4596γ-0.0238ω4+0.7523 (13)
δ34=-0.8060T/Tc4-0.8550γ-0.0809ω4+1.1880 (14)
step 5, correcting the interaction between different components according to the formulas (15) and (16) to determine the attractive force coefficient a and the van der Waals molar volume b of the thickened oil-methane-carbon dioxide-propane system;
step 6, calculating the bubble point pressure P of the thickened oil-methane-carbon dioxide-propane system according to a PR state equation, namely a formula (17);
V=V1Bo/(1/x4) (19)
in the formula, P is the bubble point pressure of a thickened oil-methane-carbon dioxide-propane system, kPa; v is the molar volume of the thickened oil-methane-carbon dioxide-propane system, m3/kmol;V1In order to degas the molar volume of the thick oil, m3/kmol;
Step 7, calculating the expansion coefficient SF of the thickened oil-methane-carbon dioxide-propane system according to a formula (20);
wherein, SF is expansion coefficient and has no dimension; v2Corrected molar volume, m, for thickened oil-methane-carbon dioxide-propane system3/kmol。
According to the present invention, the dimensionless parameter α in step 3 is calculated by formula (21):
in the formula, TriThe calculation formula of (2) is as follows:
Tri=T/Tci(i=1,2,3,4) (22)
in the formula (21), the thick oil eccentricity factor omega4The calculation formula of (2) is as follows:
when converted to boiling point Tbr=TbR/Tc4When the content of the carbon dioxide is less than 0.8,
when converted to boiling point Tbr=TbR/Tc4When the content of the carbon dioxide is more than or equal to 0.8,
according to the present invention, the volume correction method for reducing the error of the expansion coefficient calculation in step 7 is shown in formula (25):
in the formula, CiFor correction factor, m3Per kmol; wherein the correction coefficient of the thickened oil is calculated according to the following formula:
in the formula ZRAFor the Rackett compression coefficient, the calculation formula is:
each gas Component (CH) in the formula (25)4、CO2And C3H8) The correction coefficient of (d) is calculated by the following formula:
Ci=si×bi(i=1,2,3) (28)
in the formula, siDimensionless parameters; e and d are constants with values of 0.182 and 2.258, respectively.
The invention has the advantages of
1. At present, most of thick oil gas-injection phase behavior researches still use single gas, the calculation method provided by the invention can accurately calculate the bubble point pressure and the expansion coefficient of a thick oil-methane-carbon dioxide-propane system, and the change characteristics of the thick oil gas-injection mixed gas phase behavior are determined, so that the method has important significance for the development of the thick oil gas-injection mixed gas.
2. The calculation method provided by the invention solves the problem that the high-pressure physical property parameters of the thickened oil-methane-carbon dioxide-propane system can only be determined by relying on indoor experiments, saves the manpower, capital and time cost, and meets the requirement of field workers for obtaining a large number of high-pressure physical property parameters of the system.
3. The calculation method adopted by the invention is simple and easy to implement, needs fewer basic physical property parameters, and is convenient for field workers to quickly obtain the requirements of the high-pressure physical property parameters of the thickened oil-methane-carbon dioxide-propane system in the initial development stage of injecting the mixed gas into the oil field.
4. The high-pressure physical property parameter of the thickened oil-methane-carbon dioxide-propane system obtained by the calculation method has important significance for researching the exploitation mechanism of the thickened oil injected mixed gas, optimizing the injection and production parameters, designing a shaft and a ground oil pipeline and the like, so the method has wide practicability in the field of the development of the thickened oil reservoir injected mixed gas.
Drawings
FIG. 1 is a flow chart for calculating high-pressure physical property parameters of a thickened oil-methane-carbon dioxide-propane system.
FIG. 2 is a graph comparing the bubble point pressure and expansion coefficient calculated for a thickened oil-methane-carbon dioxide-propane system with experimental values for different gas mole fractions in an example of the invention.
FIG. 3 is a graph of regression analysis of calculated bubble point pressure and experimental values in an example embodiment of the present invention.
FIG. 4 is a graph of regression analysis of calculated expansion coefficients and experimental values in an example embodiment of the present invention.
Detailed Description
The invention is explained in more detail below with reference to exemplary embodiments and the drawing of the description, without limiting the scope of protection.
In the examples, the experiments of methane-carbon dioxide-propane mixed gas phase behavior of heavy oil injection were carried out using degassed heavy oil samples produced from venezuela Orinoco reservoirs. The reliability of the method is verified by comparing the calculation result of the method with the high-pressure physical property experiment result of a thickened oil-methane-carbon dioxide-propane system.
Example 1
A method for determining high-pressure physical property parameters of a thickened oil-methane-carbon dioxide-propane system is disclosed, wherein a calculation flow chart is shown in figure 1, and calculation is carried out according to the following steps:
step 1, determining the oil reservoir temperature T, the heavy oil specific gravity gamma and the relative molecular mass M of each component of the systemi(i ═ 1,2,3,4, respectively representing propane, carbon dioxide, methane, heavy oil), crude oil volume coefficient BoAnd the mole fraction x of the methane-carbon dioxide-propane mixed gas in the system is givengMolar fraction N of each gas component in the gas mixturei(i is 1,2,3), and C in the system is calculated according to the formulas (1) to (4)3H8、CO2、CH4And the mole fraction x of the heavy oil componenti
x1=N1xg(1)
x2=N2xg(2)
x3=N3xg(3)
x4=1-xg(4)
In this embodiment, the reservoir temperature T is 327.35K; the specific gravity gamma of the thickened oil is 0.981; in the system C3H8Relative molecular mass M of1Is 44.10; CO 22Relative molecular mass M of244.01; CH (CH)4Relative molecular mass M of3Is 16.04; relative molecular mass M of thickened oil4557.81; volume coefficient of crude oil BoIs 1.17; the mole fraction x of the methane-carbon dioxide-propane mixed gas in the systemg10.29 percent; c in the mixed gas3H8Mole fraction N of1Is 28%; CO 22Mole fraction N of28 percent; CH (CH)4Mole fraction N of3The content was 64%.
C in the system is obtained by calculation according to the formulas (1) to (4)3H8Mole fraction x of12.88%; CO 22Mole fraction x of20.82%; CH (CH)4Mole fraction x of36.59 percent; mole fraction x of heavy oil4The content was 89.71%.
Step 2, determining critical parameters of each gas component by referring to a common gas main physical and chemical property table (oil layer physics, Leifen eds.), including critical temperature TciCritical pressure PciCritical molar volume Vci(i is 1,2,3), and calculating the critical parameters of the thickened oil according to the formulas (5) to (7);
in the formula, TbRIs the standard boiling point of the thickened oil, K; m4The relative molecular mass of the thickened oil is dimensionless; t isc4Is critical temperature of thick oil,K;Pc4The critical pressure of the thick oil is kPa.
In this embodiment, the critical parameters of the gas components are determined by looking up a table, the critical parameters of the heavy oil are calculated according to formulas (5) to (7), and the critical parameters of each component of the obtained heavy oil-methane-carbon dioxide-propane system are shown in table 1.
Table 1 critical properties table for each component in the thickened oil-methane-carbon dioxide-propane system
Step 3, determining the attraction coefficient a of each component in the thickened oil-methane-carbon dioxide-propane system according to the formulas (8) to (10)iAnd van der Waals molar volume bi
ai=aciα(Trii)(i=1,2,3,4) (8)
Wherein R is a general gas constant of 8.314kPa · m3/(K·kmol);TriThe conversion temperature of each component in the thickened oil-methane-carbon dioxide-propane system is dimensionless; omegaiThe eccentricity factors of the components in the thickened oil-methane-carbon dioxide-propane system and the eccentricity factors of the gas components can be obtained by consulting a pure material thermodynamic basic parameter table (Perry chemical engineering handbook, R.H. Perry), wherein omega is1Is 0.008 omega2Is 0.224 omega3Is 0.152; pciCritical pressure, kPa, of each component in a thickened oil-methane-carbon dioxide-propane system; t isciIs the critical temperature of each component in the thickened oil-methane-carbon dioxide-propane system, and K and α are dimensionless parameters.
The dimensionless parameter α in equation (8) is calculated by equation (21):
in the formula, TriThe calculation formula of (2) is as follows:
Tri=T/Tci(i=1,2,3,4) (22)
in the formula (21), the thick oil eccentricity factor omega4The calculation formula of (2) is as follows:
when converted to boiling point Tbr=TbR/Tc4When the content of the carbon dioxide is less than 0.8,
when converted to boiling point Tbr=TbR/Tc4When the content of the carbon dioxide is more than or equal to 0.8,
in the present example, the attraction coefficient a of each component in the thickened oil-methane-carbon dioxide-propane system was determined according to the equations (8) to (10)iAnd van der Waals molar volume biThe calculation results are shown in table 2:
TABLE 2 values of the components a and b of the system
Step 4, determining the system binary interaction coefficient delta according to the formulas (11) to (14)ijWherein δ12Is C3H8And CO2Binary coefficient of interaction, delta13Is C3H8And CH4Binary coefficient of interaction, delta23Is CO2And CH4Binary coefficient of interaction, delta14Is C3H8Coefficient of interaction, delta, with heavy oil24Is CO2Coefficient of interaction, delta, with heavy oil34Is CH4The coefficient of interaction with the thickened oil is binary;
δ14=-0.4560T/Tc4+0.1817 (12)
δ24=-0.5462T/Tc4-0.4596γ-0.0238ω4+0.7523 (13)
δ34=-0.8060T/Tc4-0.8550γ-0.0809ω4+1.1880 (14)
in the formula Vci、VcjCritical molar volumes of gas component i and component j, respectively, m3·kmol-1
In the present example, C is determined according to equations (11) to (14)3H8-CO212)、C3H8-CH413)、CO2-CH423)、C3H8Heavy oil (. delta.)14)、CO2Heavy oil (. delta.)24) And CH4Heavy oil (. delta.)34) Of the binary interaction coefficient deltaij(ii) a Calculating the binary interaction coefficient delta between the obtained system componentsijAs shown in table 3:
TABLE 3 binary interaction coefficient between components in thickened oil-methane-carbon dioxide-propane system
Step 5, correcting the interaction between different components according to Van der Waals mixing rules, namely formulas (15) and (16), and determining an attractive force coefficient a and a Van der Waals molar volume b of the thickened oil-methane-carbon dioxide-propane system;
the true bookIn the examples, the interaction between the various components was corrected according to the equations (15), (16) to give a value for the attraction coefficient a of 3.93X 10 for the thickened oil-methane-carbon dioxide-propane system4The van der waals molar volume b has a value of 0.551.
Step 6, calculating the bubble point pressure P of the thickened oil-methane-carbon dioxide-propane system according to a PR state equation, namely a formula (17);
V=V1Bo/(1/x4) (19)
in the formula, P is the bubble point pressure of a thickened oil-methane-carbon dioxide-propane system, kPa; v is the molar volume of the thickened oil-methane-carbon dioxide-propane system, m3/kmol;V1In order to degas the molar volume of the thick oil, m3/kmol。
In this example, the value of the bubble point pressure P of the thickened oil-methane-carbon dioxide-propane system was calculated to be 1.42MPa according to the PR equation of state, i.e., formula (17).
Step 7, calculating the expansion coefficient SF of the thickened oil-methane-carbon dioxide-propane system according to a formula (20);
wherein, SF is expansion coefficient and has no dimension; v2Corrected molar volume, m, for thickened oil-methane-carbon dioxide-propane system3/kmol。
The volume correction method to reduce the error of the expansion coefficient calculation is shown in equation (25):
in the formula, CiFor correction factor, m3Per kmol; wherein the correction coefficient of the thickened oil is calculated according to the following formula:
in the formula ZRAFor the Rackett compression coefficient, the calculation formula is:
each gas Component (CH) in the formula (25)4、CO2And C3H8) The correction coefficient of (d) is calculated by the following formula:
Ci=si×bi(i=1,2,3) (28)
in the formula, siDimensionless parameters; e and d are constants with values of 0.182 and 2.258, respectively.
In this example, the value of the expansion coefficient SF of the thick oil-methane-carbon dioxide-propane system was calculated to be 1.0077 according to the formula (20).
Example 2
Method for determining high-pressure physical property parameters of thickened oil-methane-carbon dioxide-propane system
The method for calculating the high-pressure physical property parameters of the thick oil-methane-carbon dioxide-propane system in example 1 is different from the method for calculating the high-pressure physical property parameters of the thick oil-methane-carbon dioxide-propane system in that the mole fraction of the methane-carbon dioxide-propane mixed gas in the thick oil is changed, as shown in table 4, the mole fraction of the methane-carbon dioxide-propane mixed gas in the system is respectively 4 different values, and the bubble point pressure and the expansion coefficient of the thick oil-methane-carbon dioxide-propane system are calculated according to the calculation steps in example 1.
Table 4 example 2 high pressure physical properties calculation results of thickened oil-methane-carbon dioxide-propane system
Examples of effects
In order to verify the accuracy of the calculation results obtained by the present invention, experimental values of bubble point pressure and expansion coefficient at 5 mole fractions of mixed gas in examples 1 and 2 were measured by a stepwise depressurization method using a high temperature and high pressure formation fluid analyzer (PVT apparatus). In the experimental process, a certain amount of mixed gas and thick oil are injected into a PVT cylinder, the pressure is increased to reach a balanced state, a thick oil-methane-carbon dioxide-propane system is formed, the volume of the PVT cylinder is increased step by step, the pressure is reduced to reach the balanced state, the pressure and the volume of the PVT cylinder are recorded, and finally a curve of the pressure along with the change of the volume is drawn, so that bubble point pressure and expansion coefficient experimental values under each gas mole fraction are obtained. The calculated values and experimental values are shown in table 5 and fig. 2. As can be seen from fig. 2, the calculated values of bubble point pressure and expansion coefficient are very small in difference from the experimental values.
TABLE 5 calculation and experiment results of high-pressure physical properties of thickened oil-methane-carbon dioxide-propane system
The accuracy of the calculation result of the high-pressure physical property parameters of the thickened oil-methane-carbon dioxide-propane system in the embodiment 1 and the embodiment 2 is quantitatively evaluated by adopting three methods of linear regression analysis, average error analysis and variance analysis.
From the results of the high pressure physical property parameter calculations of the examples, a regression analysis and an average error table of the calculated bubble point pressure and the experimental value were obtained as shown in table 6, and a regression analysis chart is shown in fig. 3. Bubble point pressure calculated values and an analysis of variance table of experimental values are shown in table 7.
TABLE 6 regression analysis and average error tables for calculated and experimental values of bubble point pressure
TABLE 7 analysis table of calculated bubble point pressure and variance of experimental values
As shown in table 6, the slope (a) of the linear relation (y ═ Ax + B) obtained by regression of the calculated bubble point pressure and the experimental values was approximately equal to 1, and the average error of the calculated bubble point pressure was 1.77%.
According to the analysis of variance method, the F value is calculated according to the following formula:
where MSA is the mean square between groups and MSE is the mean square within a group.
As shown in Table 7, the F value is 64604.8837, the significance level α is 0.05, and the critical value F is consulted according to the degrees of freedom between groups and in groupscThe distribution table can know the critical value Fc10.13 (applied statistics, main compilation for constant mega light), so the F value is much larger than the critical value Fc. Therefore, the predictive model can accurately calculate the bubble point pressure.
The results of the calculation of the high pressure physical property parameters according to the examples were used to obtain a regression analysis and an average error table of the calculated values of the expansion coefficient and the experimental values, as shown in Table 8, and a regression analysis chart is shown in FIG. 4. The calculated values of the expansion coefficient and the analysis of variance of the experimental values are shown in Table 9.
TABLE 8 calculated values of expansion coefficients and regression analysis and mean error tables of experimental values
TABLE 9 calculated values of expansion coefficient and analysis of variance of experimental values
As shown in table 8, the slope (a) of the linear relationship (y ═ Ax + B) obtained by regression of the calculated expansion coefficients and the experimental values was approximately equal to 1, and the average error of the calculated expansion coefficients was 0.07%.
According to the analysis of variance method, the F value is calculated according to the following formula:
where MSA is the mean square between groups and MSE is the mean square within a group.
As shown in Table 9, F is 1780, significance level α is 0.05, and the threshold value F is consulted according to the degree of freedom between and within groupscThe distribution table can know the critical value FcIs 10.13, so the F value is much larger than the critical value Fc. Therefore, the prediction model can accurately calculate the expansion coefficient.

Claims (3)

1. A method for determining high-pressure physical property parameters of a thickened oil-methane-carbon dioxide-propane system comprises the following steps:
step 1, determining the oil reservoir temperature T, the heavy oil specific gravity gamma and the relative molecular mass M of each component of the systemi(i ═ 1,2,3,4, respectively representing propane, carbon dioxide, methane, heavy oil), crude oil volume coefficient BoAnd the mole fraction x of the methane-carbon dioxide-propane mixed gas in the system is givengMolar fraction N of each gas component in the gas mixturei(i is 1,2,3), and C in the system is calculated according to the formulas (1) to (4)3H8、CO2、CH4And the mole fraction x of the heavy oil componenti
x1=N1xg(1)
x2=N2xg(2)
x3=N3xg(3)
x4=1-xg(4)
Step 2, determining critical parameters including critical temperature T of each gas component by looking up a common gas main physical and chemical property tableciCritical pressure PciCritical molar volume Vci(i is 1,2,3), and calculating the critical parameters of the thickened oil according to the formulas (5) to (7);
in the formula, TbRIs the standard boiling point of the thickened oil, K; t isc4Is the critical temperature of the thick oil, K; pc4Critical pressure of thick oil, kPa;
step 3, determining the attraction coefficient a of each component in the thickened oil-methane-carbon dioxide-propane system according to the formulas (8) to (10)iAnd van der Waals molar volume bi
ai=aciα(Trii) (i=1,2,3,4) (8)
Wherein R is a general gas constant of 8.314kPa · m3/(K·kmol);TriThe conversion temperature of each component in the thickened oil-methane-carbon dioxide-propane system is dimensionless; omegaiThe eccentricity factors of all components in a thickened oil-methane-carbon dioxide-propane system can be obtained by looking up a pure material thermodynamic basic parameter table; pciCritical pressure, kPa, of each component in a thickened oil-methane-carbon dioxide-propane system; t isciIs the critical temperature of each component in the thickened oil-methane-carbon dioxide-propane system, K, α is a dimensionless parameter;
step 4, determining system binary interaction according to formulas (11) to (14)By a factor deltaijWherein δ12Is C3H8And CO2Binary coefficient of interaction, delta13Is C3H8And CH4Binary coefficient of interaction, delta23Is CO2And CH4Binary coefficient of interaction, delta14Is C3H8Coefficient of interaction, delta, with heavy oil24Is CO2Coefficient of interaction, delta, with heavy oil34Is CH4The coefficient of interaction with the thickened oil is binary;
δ14=-0.4560T/Tc4+0.1817(12)
δ24=-0.5462T/Tc4-0.4596γ-0.0238ω4+0.7523(13)
δ34=-0.8060T/Tc4-0.8550γ-0.0809ω4+1.1880(14)
step 5, correcting the interaction between different components according to the formulas (15) and (16) to determine the attractive force coefficient a and the van der Waals molar volume b of the thickened oil-methane-carbon dioxide-propane system;
step 6, calculating the bubble point pressure P of the thickened oil-methane-carbon dioxide-propane system according to a PR state equation, namely a formula (17);
V=V1Bo/(1/x4)(19)
in the formula, P is the bubble point pressure of a thickened oil-methane-carbon dioxide-propane system, kPa; v is the molar volume of the thickened oil-methane-carbon dioxide-propane system, m3/kmol;V1In order to degas the molar volume of the thick oil, m3/kmol;
Step 7, calculating the expansion coefficient SF of the thickened oil-methane-carbon dioxide-propane system according to a formula (20);
wherein, SF is expansion coefficient and has no dimension; v2Corrected molar volume, m, for thickened oil-methane-carbon dioxide-propane system3/kmol。
2. The method for determining high pressure physical property parameters of thick oil-methane-carbon dioxide-propane system as claimed in claim 1, wherein the dimensionless parameter α in step 3 is calculated by equation (21):
in the formula, TriThe calculation formula of (2) is as follows:
Tri=T/Tci(i=1,2,3,4)(22)
in the formula (21), the thick oil eccentricity factor omega4The calculation formula of (2) is as follows:
when converted to boiling point Tbr=TbR/Tc4When the content of the carbon dioxide is less than 0.8,
when converted to boiling point Tbr=TbR/Tc4When the content of the carbon dioxide is more than or equal to 0.8,
3. the method for determining high pressure physical property parameters of thick oil-methane-carbon dioxide-propane system according to claim 1, wherein the volume correction method for reducing the error of expansion coefficient calculation in step 7 is shown in formula (25):
in the formula, CiFor correction factor, m3Per kmol; wherein the correction coefficient of the thickened oil is calculated according to the following formula:
in the formula ZRAFor the Rackett compression coefficient, the calculation formula is:
each gas Component (CH) in the formula (25)4、CO2And C3H8) The correction coefficient of (d) is calculated by the following formula:
Ci=si×bi(i=1,2,3)(28)
in the formula, siDimensionless parameters; e and d are constants with values of 0.182 and 2.258, respectively.
CN201911228850.6A 2019-12-04 2019-12-04 Method for determining high-pressure physical property parameters of thickened oil-methane-carbon dioxide-propane system Active CN111081320B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911228850.6A CN111081320B (en) 2019-12-04 2019-12-04 Method for determining high-pressure physical property parameters of thickened oil-methane-carbon dioxide-propane system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911228850.6A CN111081320B (en) 2019-12-04 2019-12-04 Method for determining high-pressure physical property parameters of thickened oil-methane-carbon dioxide-propane system

Publications (2)

Publication Number Publication Date
CN111081320A true CN111081320A (en) 2020-04-28
CN111081320B CN111081320B (en) 2021-06-15

Family

ID=70312949

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911228850.6A Active CN111081320B (en) 2019-12-04 2019-12-04 Method for determining high-pressure physical property parameters of thickened oil-methane-carbon dioxide-propane system

Country Status (1)

Country Link
CN (1) CN111081320B (en)

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101497805A (en) * 2009-03-13 2009-08-05 北京化工大学 Method for preparing biodiesel by catalytic cracking acidified solid catalysis
US20090250224A1 (en) * 2008-04-04 2009-10-08 Halliburton Energy Services, Inc. Phase Change Fluid Spring and Method for Use of Same
CN101555500A (en) * 2009-05-22 2009-10-14 天津大学 Method for preparing diglyceride by enzyme catalysis in supercritical fluids system
CN102448580A (en) * 2009-03-25 2012-05-09 氟石科技公司 Improved configurations and methods for high pressure acid gas removal
US20140190074A1 (en) * 2013-01-09 2014-07-10 C2O Technologies, Llc Process for treating coal to improve recovery of condensable coal derived liquids
CN105481211A (en) * 2014-09-20 2016-04-13 海安华达石油仪器有限公司 Ultrasonic-strengthened supercritical extraction method for oily sludge
CN106837271A (en) * 2017-04-19 2017-06-13 中国石油大学(华东) A kind of computational methods for determining foam oil physical data under high pressure
CN107652963A (en) * 2017-10-13 2018-02-02 中国石油大学(华东) A kind of gas foam foam stabilizer system and preparation method thereof

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090250224A1 (en) * 2008-04-04 2009-10-08 Halliburton Energy Services, Inc. Phase Change Fluid Spring and Method for Use of Same
CN101497805A (en) * 2009-03-13 2009-08-05 北京化工大学 Method for preparing biodiesel by catalytic cracking acidified solid catalysis
CN102448580A (en) * 2009-03-25 2012-05-09 氟石科技公司 Improved configurations and methods for high pressure acid gas removal
CN101555500A (en) * 2009-05-22 2009-10-14 天津大学 Method for preparing diglyceride by enzyme catalysis in supercritical fluids system
US20140190074A1 (en) * 2013-01-09 2014-07-10 C2O Technologies, Llc Process for treating coal to improve recovery of condensable coal derived liquids
CN105481211A (en) * 2014-09-20 2016-04-13 海安华达石油仪器有限公司 Ultrasonic-strengthened supercritical extraction method for oily sludge
CN106837271A (en) * 2017-04-19 2017-06-13 中国石油大学(华东) A kind of computational methods for determining foam oil physical data under high pressure
CN107652963A (en) * 2017-10-13 2018-02-02 中国石油大学(华东) A kind of gas foam foam stabilizer system and preparation method thereof

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
BEI LIU ET AL.: "Evaluation of different CH4-CO2 replacement processes in hydrate-bearing sediments by measuring p-wave velocity", 《ENERGICS》 *
张勇: "苏北低渗油藏二氧化碳混相驱油方案优化研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅰ辑》 *

Also Published As

Publication number Publication date
CN111081320B (en) 2021-06-15

Similar Documents

Publication Publication Date Title
CN108959767B (en) Method for simulating damage numerical values of condensate oil in different well types of narrow-river-channel gas reservoir
CN103114830B (en) Enriched-gas-drive water-altering-gas (WAG) injection method
Vulin et al. Slim-tube simulation model for CO2 injection EOR
CN103726815A (en) CO2 flooding extraction well shaft flow state determining and parameter optimizing method
Lei et al. Effects of boundary layer and stress sensitivity on the performance of low-velocity and one-phase flow in a shale oil reservoir: Experimental and numerical modeling approaches
CN111081320B (en) Method for determining high-pressure physical property parameters of thickened oil-methane-carbon dioxide-propane system
Cook et al. A beta-type reservoir simulator for approximating compositional effects during gas injection
CN109296363B (en) Ultra-low permeability reservoir carbon dioxide flooding initial-stage capacity prediction method
Hou et al. Phase behavior and physical parameters of natural gas mixture with
Sheng et al. An experimental study of steam-assisted gravity drainage
CN112630090A (en) Method for calculating diffusion coefficient of methane-carbon dioxide-propane mixed gas by considering volume expansion of thickened oil
Hongwei et al. Experimental and numerical simulation study on single layer injectivity for CO2 flooding in low permeability oil reservoir
CN109188557B (en) Method for calculating hydrocarbon discharge efficiency by using hydrocarbon generation pressurization model
Lins et al. Three-phase hydrocarbon thermodynamic liquid-liquid-vapour equilibrium in CO2 process
CN108694255B (en) Supercritical gas phase displacement medium-condensate oil gas equilibrium phase behavior evaluation method
Zhao et al. The influence of negative pressure during gas extraction in low-permeability coal seams on the effect of gas extraction
Xu et al. Experimental investigation and numerical simulation of dynamic characteristics for multithermal fluid-assisted SAGD in extraheavy oil reservoir
Shen et al. Optimization Strategy to Reduce Asphaltene Deposition-Associated Damage During CO $$ _2 $$2 Huff-n-Puff Injection in Shale
Patil et al. Phase behavior, solid organic precipitation, and mobility characterization studies in support of enhanced heavy oil recovery on the Alaska North Slope
Jafari et al. Experimental study and simulation of different EOR techniques in a non-fractured carbonate core from an Iranian offshore oil reservoir
CN110940796B (en) Continuous testing method for gas injection expansion of crude oil
CN111472764B (en) Method for calculating recovery rate of rich gas in rich gas flooding process
Fan et al. Study on the change law of dew point pressure after gas injection in fracture-cavity condensate gas reservoir
CN107975356B (en) Oilfield gas injection pressure prediction method
Sun et al. Molecular Dynamics Simulation of Nanoconfined Methane Flow Through Realistic Organic Shale Matrix

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant