CN116258094A - Oil reservoir air injection flooding feasibility judging method based on oxidation kinetic parameters - Google Patents
Oil reservoir air injection flooding feasibility judging method based on oxidation kinetic parameters Download PDFInfo
- Publication number
- CN116258094A CN116258094A CN202310118126.8A CN202310118126A CN116258094A CN 116258094 A CN116258094 A CN 116258094A CN 202310118126 A CN202310118126 A CN 202310118126A CN 116258094 A CN116258094 A CN 116258094A
- Authority
- CN
- China
- Prior art keywords
- oxidation
- crude oil
- stage
- temperature
- conversion rate
- 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
Links
- 238000007254 oxidation reaction Methods 0.000 title claims abstract description 222
- 230000003647 oxidation Effects 0.000 title claims abstract description 199
- 238000000034 method Methods 0.000 title claims abstract description 65
- 238000002347 injection Methods 0.000 title claims abstract description 26
- 239000007924 injection Substances 0.000 title claims abstract description 26
- 239000010779 crude oil Substances 0.000 claims abstract description 123
- 238000006243 chemical reaction Methods 0.000 claims abstract description 101
- 230000004913 activation Effects 0.000 claims abstract description 58
- 238000010438 heat treatment Methods 0.000 claims abstract description 26
- 239000003921 oil Substances 0.000 claims abstract description 26
- 230000008021 deposition Effects 0.000 claims description 34
- 239000000446 fuel Substances 0.000 claims description 34
- 238000002411 thermogravimetry Methods 0.000 claims description 14
- 230000007246 mechanism Effects 0.000 claims description 12
- 230000035484 reaction time Effects 0.000 claims description 10
- 238000004364 calculation method Methods 0.000 claims description 9
- 230000000630 rising effect Effects 0.000 claims description 9
- 238000002474 experimental method Methods 0.000 claims description 8
- 230000008569 process Effects 0.000 claims description 6
- 238000012935 Averaging Methods 0.000 claims description 4
- 230000006870 function Effects 0.000 description 9
- 239000000295 fuel oil Substances 0.000 description 5
- 239000007789 gas Substances 0.000 description 5
- 230000008901 benefit Effects 0.000 description 3
- 238000011161 development Methods 0.000 description 3
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 description 2
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000010354 integration Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 239000001301 oxygen Substances 0.000 description 2
- 229910052760 oxygen Inorganic materials 0.000 description 2
- BASFCYQUMIYNBI-UHFFFAOYSA-N platinum Chemical compound [Pt] BASFCYQUMIYNBI-UHFFFAOYSA-N 0.000 description 2
- 238000010926 purge Methods 0.000 description 2
- 230000007704 transition Effects 0.000 description 2
- UGFAIRIUMAVXCW-UHFFFAOYSA-N Carbon monoxide Chemical compound [O+]#[C-] UGFAIRIUMAVXCW-UHFFFAOYSA-N 0.000 description 1
- 238000005452 bending Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000002485 combustion reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000003546 flue gas Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 229910052757 nitrogen Inorganic materials 0.000 description 1
- 229910052697 platinum Inorganic materials 0.000 description 1
- 239000000243 solution Substances 0.000 description 1
- 230000002269 spontaneous effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Images
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B43/00—Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
- E21B43/16—Enhanced recovery methods for obtaining hydrocarbons
- E21B43/166—Injecting a gaseous medium; Injecting a gaseous medium and a liquid medium
- E21B43/168—Injecting a gaseous medium
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Theoretical Computer Science (AREA)
- Fluid Mechanics (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- General Physics & Mathematics (AREA)
- Computing Systems (AREA)
- Pure & Applied Mathematics (AREA)
- Geochemistry & Mineralogy (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Algebra (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Environmental & Geological Engineering (AREA)
- Production Of Liquid Hydrocarbon Mixture For Refining Petroleum (AREA)
Abstract
The application discloses a method for judging feasibility of oil reservoir air injection driving based on oxidation kinetic parameters, which comprises the steps of dividing three oxidation reaction stages according to distribution characteristics of oxidation characteristic curves of crude oil at a plurality of different heating rates, respectively obtaining respective different conversion rate values of the three oxidation reaction stages, then obtaining activation energy corresponding to the three oxidation stages according to the oxidation characteristic curves by combining expression of conversion rate methods such as OFW and the like, simultaneously obtaining pre-pointing factors of the crude oil at the different oxidation stages according to an Arrhenius equation, and comparing the obtained activation energy with corresponding preset data respectively to judge the feasibility of oil reservoir air injection driving. Therefore, on the basis of optimizing an oxidation dynamics parameter model, the method solves the accurate activation energy and factor parameters before finger, and thereby provides a reliable reference index for predicting the air injection flooding feasibility of the oil reservoir.
Description
Technical Field
The application relates to the technical field of energy exploitation, in particular to a reservoir air injection flooding feasibility judging method based on oxidation kinetic parameters.
Background
Air flooding technology has gradually attracted wide interest in industry in recent years due to the unique advantages of sufficient air source, high economy, strong injection capability, flue gas flooding mechanism, capability of forming low-temperature oxidation or high-temperature oxidation heat release under formation conditions, improvement of crude oil properties, and the like. Light oil reservoir air flooding and heavy oil reservoir fire flooding gradually become leading edge technologies for field test attack. The difference of the air flooding of the light oil reservoir and the heavy oil reservoir is mainly reflected in the difference of contact relation between oxygen and crude oil of stratum and different oxidation mechanisms. The light crude oil under the reservoir condition can generate a middle-low temperature spontaneous oxidation reaction, and the heavy crude oil under the high temperature condition of artificial ignition can generate a combustion reaction. The oxidation kinetic parameters such as activation energy, pre-factor and the like are important indexes for representing the difficulty of the air injection flooding oxidation effect, and become key parameters for predicting the feasibility of air injection flooding of the oil reservoir.
The currently commonly used oxidation kinetics calculation method mainly comprises a Coats-Redfren integration method, an OFW conversion rate method and an Arrhenius equation, but each method has respective advantages and disadvantages. The method comprises the steps of obtaining kinetic parameters such as activation energy, pre-finger factors and the like through a single heating rate and multiple conversion functions by using a Coats-Redfren integration method. However, the method obtains oxidation kinetic parameters only by fitting a transformation rate gradient relation of a TG/DTG curve, and the result is not accurate enough. The Arrhenius equation is a classical formula of oxidation kinetic parameters, and the method obtains kinetic parameters such as activation energy, pre-finger factors and the like through a crude oil oxidation characteristic curve at a single heating rate. However, from a quantitative point of view, this method is not entirely accurate. Because the assumptions in the formula regarding the independence of activation energy from temperature are not strict, the resulting straight lines from the mapping of oxidation kinetics parameters often exhibit varying degrees of bending. The conversion rate methods such as Ozawa-Flynn-Wall (OFW) and the like do not need to select reaction mechanism functions in advance, a thermogravimetric curve is not required to be recorded accurately, parameters such as activation energy and the like are not influenced by the temperature rising rate, errors caused by the parameters are avoided, but the method simplifies the reaction mechanism functions and brings certain difficulty to calculation of factors before finger.
In summary, because the oxidation reaction mechanism of the air-injected crude oil is complex, oxidation kinetic parameters obtained by calculation of numerous oxidation kinetic models at present are different, and the matched reaction mechanism model function may not be completely suitable for the actual oxidation reaction process of the air-injected crude oil. Therefore, the accurate parameters such as activation energy, pre-finger factors and the like must be solved on the basis of optimizing an oxidation dynamics parameter model, and a reference index is provided for predicting the air injection flooding feasibility of the oil reservoir.
Disclosure of Invention
The application aims to provide an oil reservoir air injection flooding feasibility judging method based on oxidation kinetic parameters, which is used for solving the defects in the prior art.
In order to achieve the above purpose, the present application provides a method for determining the feasibility of reservoir air flooding based on oxidation kinetic parameters, comprising the steps of:
s10: dividing the crude oil oxidation reaction into three oxidation stages according to the distribution characteristics of oxidation characteristic curves of the crude oil at a plurality of different heating rates, wherein the three oxidation stages are respectively a low-temperature oxidation stage, a fuel deposition stage and a high-temperature oxidation stage;
s20: obtaining a plurality of characteristic values of the conversion rate of crude oil in three oxidation stages, determining the corresponding characteristic temperature of the characteristic conversion rate of a single oxidation stage in other temperature rise rate experiments, calculating the reaction rate value in the crude oil oxidation reaction process, and making a relation graph of ln (beta) changing along with 1/T, wherein T represents absolute temperature and is expressed as K; beta represents the temperature rising rate, and the unit is ℃/min;
s30: according to the oxidation characteristic curve, the activation energy corresponding to the three oxidation stages is obtained by combining the expression of the conversion rate method such as OFW and the like, and meanwhile, the pre-finger factors of crude oil in different oxidation stages are obtained according to an Arrhenius equation;
s40: and comparing the obtained activation energy and the pre-finger factor with corresponding preset data respectively to judge the feasibility of air injection driving of the oil reservoir.
In one possible implementation manner, the step S10 includes the steps of:
s100: the crude oil oxidation reaction thermogravimetry experiment is repeated for a plurality of times, crude oil oxidation thermogravimetry TG and DTG experimental data under each heating rate are derived, and a crude oil TG curve graph and a crude oil DTG curve graph are respectively manufactured;
s110: fitting all the crude oil TG curve graphs and the crude oil DTG curve graphs to obtain an oxidation characteristic curve, and dividing the low-temperature oxidation stage, the fuel deposition stage and the high-temperature oxidation stage according to the turning points in the oxidation characteristic curve.
In one possible embodiment, three of the oxidation stages are divided by the oxidation profile of the crude oil at least four different ramp rates.
In one possible implementation manner, the step S20 includes the steps of:
s200: dividing intervals of the conversion rate corresponding to a crude oil TG curve graph and a crude oil DTG curve of crude oil in the low-temperature oxidation stage, the fuel deposition stage and the high-temperature oxidation stage respectively according to the oxidation characteristic curve;
s210: and selecting the mass corresponding to the crude oil with the start-stop temperatures of the three oxidation stages in the full temperature range under one heating rate, so as to calculate the conversion rate gamma value at the reaction time t of the low-temperature oxidation stage, the fuel deposition stage and the high-temperature oxidation stage.
In one possible implementation manner, the step S210 includes the steps of:
s211: selecting the quality Z of crude oil with different start and stop temperatures in different oxidation stages in the full temperature range at one heating rate, wherein the quality Z is the quality of the crude oil in the original state in the low-temperature oxidation stage Li Quality Z of crude oil in ending state of low-temperature oxidation stage Lw Quality Z of raw state of crude oil fuel deposition stage Fi Quality Z of crude oil fuel deposition end state Fw Quality Z of crude oil in original state of high-temperature oxidation Hi Quality Z of crude oil in high temperature oxidation end state Hw ;
S212: will Z Li 、Z Lw 、Z Fi 、Z Fw 、Z Hi And Z Hw Substituting the conversion rate gamma of the low-temperature oxidation stage, the fuel deposition stage and the high-temperature oxidation stage at the reaction time t L Numerical value, gamma F Numerical value and gamma H In the numerical calculation formula:
for determining 20 characteristic values of the conversion gamma of crude oil in three oxidation stages from 0.1 to 1 in a gradient of 0.05, 60 characteristic values in total for the three oxidation stages, wherein Z Lt 、Z Ft And Z Ht And the crude oil quality corresponding to the reaction time t of the three oxidation stages is respectively obtained.
In a possible implementation manner, the step S210 further includes the steps of:
s213: when the activation energy of the low-temperature oxidation stage is calculated, the conversion rate gamma is 0-0.9 for eighteen data points; when the activation energy of the fuel deposition stage is calculated, the conversion rate gamma is 0.2-0.9 for sixteen data points; when the activation energy of the high-temperature oxidation stage is calculated, the conversion rate gamma is 0.1-1 eighteen data points altogether.
In one possible embodiment, the reaction rate value during the oxidation reaction of the crude oil in step S20 is dγ/dT.
In one possible implementation, the step S30 includes the steps of:
s300: according to the oxidation characteristic curve, fitting the slope of a straight line, substituting the slope into a solving formula of conversion rate methods such as OFW and the like, and solving the activation energy corresponding to each characteristic conversion rate;
s310: and carrying out activation energy averaging treatment on the activation energy corresponding to each characteristic conversion rate obtained in each oxidation stage, thereby obtaining the activation energy corresponding to the low-temperature oxidation stage, the fuel deposition stage and the high-temperature oxidation stage.
In one possible embodiment, the expression of the conversion method such as OFW is:
wherein G (gamma) represents a reaction mechanism integral function; beta represents the temperature rising rate, and the unit is ℃/min; a represents a factor before finger, in min -1 The method comprises the steps of carrying out a first treatment on the surface of the E represents reaction activation energy, and the unit is KJ/mol; r represents a molar gas constant of 8.314J/(mol.K); t represents absolute temperature in K.
In one possible embodiment, the determining the pre-finger factor of the crude oil at different oxidation stages according to the arrhenius equation comprises:
based on crude oil thermogravimetric experiment data at a single heating rate, fittingAnd obtaining a fitting straight line of each reaction stage by a 1/T relation, and determining pre-finger factors of different oxidation stages by combining the Arrhenius equation through the intercept of a linear curve and a Y axis.
Compared with the prior art, the beneficial effect of this application:
according to the oil reservoir air injection driving feasibility judging method based on the oxidation kinetic parameters, the activation energy is obtained through calculation according to the corresponding relation of the same conversion rate under a plurality of heating rates by a conventional OFW conversion rate method, and meanwhile, the distribution characteristics of different conversion rates and activation energy can be obtained. Therefore, the selection of the conversion rate in the conversion rate method such as OFW is important for the accuracy of the activation energy parameter. However, the conventional OFW conversion method uses the whole reaction zone of crude oil oxidation in the full temperature range as the selection range of the conversion rate gamma, and the calculated activation energy has low precision and poor referenceability. The method for selecting the conversion rate of the traditional OFW and other conversion rate methods is improved, namely three oxidation reaction stages are divided according to the distribution characteristics of oxidation characteristic curves of crude oil at a plurality of different heating rates, respectively obtaining different conversion rate values of the three oxidation reaction stages, then according to the oxidation characteristic curves, and according to the expression of the OFW and other conversion rate methods, the activation energy corresponding to the three oxidation stages is obtained, and meanwhile, the pre-finger factors of the crude oil in the different oxidation stages are obtained according to an Arrhenius equation. Therefore, on the basis of optimizing an oxidation kinetic parameter model, the method solves the accurate activation energy and factor parameters before finger, and provides a reliable reference index for predicting the air injection flooding feasibility of the oil reservoir.
Additional features and advantages of the present application will be set forth in the detailed description which follows.
Drawings
The accompanying drawings, which are included to provide a further understanding of the application and are incorporated in and constitute a part of this specification, illustrate only some embodiments of the application and are therefore not to be considered limiting of its scope, for the purpose of providing additional related drawings from which the invention may be practiced by those of ordinary skill in the art without the exercise of inventive faculty. In the drawings:
FIG. 1 shows a graph of crude TG at 4 ramp rates;
FIG. 2 shows crude oil DTG graphs at 4 ramp rates;
FIG. 3 shows oxidation profile at a 10 ℃/min heating rate of crude oil;
FIG. 4 shows the oxidation profile at a 10℃/min heating rate of crude oil for improved conversion selection criteria;
fig. 5 shows a graph of the pre-finger factor calculated by the arrhenius equation.
Detailed Description
The following describes in detail the implementation of the embodiments of the present application with reference to the accompanying drawings. It should be understood that the detailed description is presented herein by way of illustration and explanation of the present application examples, and is not intended to limit the present application examples.
It should be noted that, in the case of no conflict, the embodiments and features in the embodiments may be combined with each other.
In the embodiments of the present application, it should be understood that the terms "center", "longitudinal", "lateral", "length", "width", "thickness", "upper", "lower", "front", "rear", "left", "right", "vertical", "horizontal", "top", "bottom", "inner", "outer", "clockwise", "counterclockwise", "axial", "radial", "circumferential", etc. indicate orientations or positional relationships based on the orientations or positional relationships shown in the drawings are merely for convenience in describing the present application and simplifying the description, and do not indicate or imply that the device or element being referred to must have a specific orientation, be configured and operated in a specific orientation, and therefore should not be construed as limiting the present application.
Furthermore, the terms "first," "second," and the like, are used for descriptive purposes only and are not to be construed as indicating or implying a relative importance or implicitly indicating the number of technical features indicated. Thus, a feature defining "a first" or "a second" may explicitly or implicitly include one or more such feature. In the description of the present application, the meaning of "a plurality" is two or more, unless explicitly defined otherwise.
In this application, unless specifically stated and limited otherwise, the terms "mounted," "connected," "secured," and the like are to be construed broadly, and may be, for example, fixedly connected, detachably connected, or integrally formed; can be mechanically or electrically connected; can be directly connected or indirectly connected through an intermediate medium, and can be communicated with the inside of two elements or the interaction relationship of the two elements. The specific meaning of the terms in this application will be understood by those of ordinary skill in the art as the case may be.
The present application will be described in detail below with reference to the attached drawings in conjunction with exemplary embodiments.
Example 1
The embodiment provides a method for judging the feasibility of oil reservoir air injection flooding based on oxidation kinetic parameters, which is hereinafter collectively called as a judging method, and comprises the following steps:
s10: dividing the crude oil oxidation reaction into three oxidation stages according to the distribution characteristics of oxidation characteristic curves of the crude oil at a plurality of different heating rates, wherein the three oxidation stages are respectively a low-temperature oxidation stage, a fuel deposition stage and a high-temperature oxidation stage;
s20: obtaining a plurality of characteristic values of the conversion rate of crude oil in three oxidation stages, determining the corresponding characteristic temperature of the characteristic conversion rate of a single oxidation stage in other temperature rise rate experiments, calculating the reaction rate value in the crude oil oxidation reaction process, and making a relation graph of ln (beta) changing along with 1/T, wherein T represents absolute temperature and is expressed as K; beta represents the temperature rising rate, and the unit is ℃/min; wherein, the reaction rate value in the crude oil oxidation reaction process is dgamma/dT value;
s30: according to the oxidation characteristic curve, the activation energy corresponding to the three oxidation stages is obtained by combining the expression of the conversion rate method such as OFW and the like, and meanwhile, the pre-finger factors of crude oil in different oxidation stages are obtained according to an Arrhenius equation;
s40: and comparing the obtained activation energy and the pre-finger factor with corresponding preset data respectively to judge the feasibility of air injection driving of the oil reservoir.
It can be understood that, comparing the oxidation kinetic parameters of the oil injection gas flooding development successfully developed, the light oil is compared with the kinetic parameters of the low-temperature oxidation stage, and the heavy oil is compared with the kinetic parameters of the high-temperature oxidation stage, if the activation energy parameter is smaller than the oxidation kinetic parameters of the oil deposit successfully developed by air injection, the oil deposit is indicated to be more suitable for air injection development.
It should be noted that the crude oils mentioned in this example are all crude oil samples, and three of the oxidation stages are divided by the oxidation profile of the crude oil at least four different heating rates.
Further, the step S10 further includes the following steps:
s100: the crude oil oxidation reaction thermogravimetry experiment is repeated for a plurality of times, crude oil oxidation thermogravimetry TG and DTG experimental data under each heating rate are derived, and a crude oil TG curve graph and a crude oil DTG curve graph are respectively manufactured;
s110: fitting all the crude oil TG curve graphs and the crude oil DTG curve graphs to obtain an oxidation characteristic curve, and dividing the low-temperature oxidation stage, the fuel deposition stage and the high-temperature oxidation stage according to the turning points in the oxidation characteristic curve.
The step S20 includes the following steps:
s200: dividing intervals of the conversion rate corresponding to a crude oil TG curve chart and a crude oil DTG curve chart of crude oil in the low-temperature oxidation stage, the fuel deposition stage and the high-temperature oxidation stage respectively according to the oxidation characteristic curves, so that the intervals of the conversion rate are more accurate;
s210: and selecting the mass corresponding to the crude oil with the start-stop temperature of the three oxidation stages in the full temperature range under one heating rate, and calculating the conversion rate gamma value at the moment of the reaction time t (unit s) of the low-temperature oxidation stage, the fuel deposition stage and the high-temperature oxidation stage.
Specifically, the step S210 includes:
s211: selecting the quality Z of crude oil with different start and stop temperatures in different oxidation stages in the full temperature range at one heating rate, wherein the quality Z is the quality of the crude oil in the original state in the low-temperature oxidation stage Li Quality Z of crude oil in ending state of low-temperature oxidation stage Lw Quality Z of raw state of crude oil fuel deposition stage Fi Quality Z of crude oil fuel deposition end state Fw Quality Z of crude oil in original state of high-temperature oxidation Hi Quality Z of crude oil in high temperature oxidation end state Hw ;
S212: will Z Li 、Z Lw 、Z Fi 、Z Fw 、Z Hi And Z Hw Substituting the conversion rate gamma of the low-temperature oxidation stage, the fuel deposition stage and the high-temperature oxidation stage at the reaction time t L Numerical value, gamma F Numerical value and gamma H In the numerical calculation formula:
for determining 20 characteristic values of the conversion gamma of crude oil in three oxidation stages from 0.1 to 1 in a gradient of 0.05, 60 characteristic values in total for the three oxidation stages, wherein Z Lt 、Z Ft And Z Ht And the crude oil quality corresponding to the reaction time t of the three oxidation stages is respectively obtained.
Further, in order to exclude the influence of the transition zone on the accuracy of the kinetic parameters, the following step S213 is further performed to exclude the influence of the transition zone when calculating the activation energy of the low-temperature oxidation stage, the activation energy of the fuel deposition stage, and the activation energy of the high-temperature oxidation stage. Wherein:
s213: when the activation energy of the low-temperature oxidation stage is calculated, the conversion rate gamma is 0-0.9 for eighteen data points; when the activation energy of the fuel deposition stage is calculated, the conversion rate gamma is 0.2-0.9 for sixteen data points; when the activation energy of the high-temperature oxidation stage is calculated, the conversion rate gamma is 0.1-1 eighteen data points altogether.
The step S30 includes the following steps:
s300: according to the oxidation characteristic curve, fitting the slope of a straight line, substituting the slope into a solving formula of conversion rate methods such as OFW and the like, and solving the activation energy corresponding to each characteristic conversion rate;
s310: and carrying out activation energy averaging treatment on the activation energy corresponding to each characteristic conversion rate obtained in each oxidation stage, thereby obtaining the activation energy corresponding to the low-temperature oxidation stage, the fuel deposition stage and the high-temperature oxidation stage.
Wherein the reaction rate in the oxidation reaction process of the air-injected crude oil is
Wherein k represents a reaction rate constant in mg/s; f (gamma) is a function of the reaction mechanism.
The reaction mechanism function f (gamma) is simplified into an n-order exponential form, and the expression is that
f(γ)=(1-γ) n #(5)
In the formula, n represents a reaction progression, and the value thereof is usually in the range of 0 to 1. The method improves the conversion rate method of OFW and the like, and the conversion rate method of OFW and the like does not need to select the reaction mechanism function f (gamma) in advance, so that the time cost is greatly saved, and the calculation step is simplified.
The reaction rate constant K is expressed as
Wherein A represents a pre-finger factor, and the unit is min-1; e represents reaction activation energy, and the unit is KJ/mol; r represents a molar gas constant of 8.314J/(mol.K); t represents absolute temperature in K.
The integral of the term transfer of the formula (6) is taken, and the expression of the conversion method such as OFW is deduced as follows:
wherein G (gamma) represents a reaction mechanism integral function; beta represents the temperature rising rate, and the unit is ℃/min; a represents a factor before finger, in min -1 The method comprises the steps of carrying out a first treatment on the surface of the E represents reaction activation energy, and the unit is KJ/mol; r represents a molar gas constant of 8.314J/(mol.K); t represents absolute temperature in K.
Therefore, according to the oxidation characteristic curve, the slope of the fitting straight line is substituted into a solving formula of the conversion rate method such as OFW (the slope of the fitting straight line corresponds to the slope in the conversion rate method expression such as OFW) so as to obtain the activation energy corresponding to each characteristic conversion rate.
Further, the calculating the pre-finger factors of the crude oil in different oxidation stages according to the Arrhenius equation comprises:
based on crude oil thermogravimetric experiment data at a single heating rate, fittingAnd obtaining a fitting straight line of each reaction stage by a 1/T relation, and determining pre-finger factors of different oxidation stages by combining the Arrhenius equation through the intercept of a linear curve and a Y axis.
The final expression of the Arrhenius equation is shown in formula (9).
Wherein t represents reaction time, and the unit is s; a represents a factor before finger, in min -1 The method comprises the steps of carrying out a first treatment on the surface of the E represents reaction activation energy, and the unit is KJ/mol; r represents a molar gas constant of 8.314J/(mol.K); t represents absolute temperature, and the unit is K; beta represents the rate of temperature rise in units of ℃/min.
Example two
In order to describe the technical scheme of the application more clearly, the embodiment provides an illustration of a reservoir air injection flooding feasibility judging method based on oxidation kinetic parameters, which is specifically as follows:
carrying out crude oil oxidation reaction thermogravimetric experiments: about 20mg of crude oil samples are weighed in each group of experiments and placed in a platinum tray, oxygen purging and nitrogen purging are set to be 50ml/min under the atmospheric pressure, 4 heating rates of 5, 10, 15 and 20 ℃/min are adopted respectively, the temperature setting range is 25-750 ℃, and the mass loss condition of the crude oil is recorded in real time. The experimental results were recorded accurately with a TG 50 thermogravimetry and analysis software was introduced. All thermogravimetric experiments were repeated 3 times to ensure experimental data accuracy. And sequentially leading out the crude oil oxidation thermal weight TG and DTG experimental data at 4 heating rates of 5, 10, 15 and 20 ℃/min, and making a crude oil TG curve graph (shown in figure 1) and a crude oil DTG curve graph (shown in figure 2).
As shown in fig. 3, the oxidation characteristic curve is obtained by fitting the crude oil TG graph and the crude oil DTG graph at a temperature rising rate of 10 ℃/min. The temperature ranges of three oxidation stages can be divided into a low-temperature oxidation stage temperature interval (25-402 ℃), a fuel deposition stage temperature interval (402-505 ℃) and a high-temperature oxidation stage temperature interval (505-616 ℃) at the inflection points of the crude oil oxidation thermogravimetry TG and DTG curve graphs.
The quality of crude oil samples at 25℃was recorded as Z Li (quality of crude oil in its original state by low-temperature oxidation), 402 ℃ quality as Z Lw (quality of crude oil at the end of low-temperature oxidation) and Z Fi (quality of crude oil fuel deposition raw state); mass at 505℃as Z Fw (quality of crude oil fuel deposition end state), Z Hi (quality of crude oil in the original state of high-temperature oxidation); quality at 616 ℃ as Z Hw (quality of the high-temperature oxidation end state of crude oil). Substituting the formulas (1), (2) and (3) to obtain 20 characteristic values of conversion rate gamma of the oxidation sample in 0.05 to 1 in a gradient manner, wherein the total of 60 characteristic values of the three oxidation stages.
Wherein, when the activation energy calculation of the low-temperature oxidation stage is carried out, the conversion rate gamma is 0-0.9 and eighteen data points are selected. Sixteen data points are sequentially added to the fuel deposition stage (alpha is 0.2-0.9), eighteen data points are added to the high-temperature oxidation stage (gamma is 0.1-1), as shown in fig. 4 (only the characteristic gamma value of the conversion rate is marked from 0.1 to 1 in a gradient of 0.1).
And determining the corresponding characteristic temperatures of the characteristic conversion rates of different oxidation stages in other temperature rise rate experiments (3 temperature rise rates of 5, 15 and 20 ℃/min), calculating dgamma/dT values, and making a relation graph of ln (beta) changing along with 1/T. On the basis, the slope of the straight line is fitted and then substituted into the formula (7), so that the activation energy corresponding to each characteristic conversion rate can be obtained.
Then respectively obtaining eighteen data points selected from the low-temperature oxidation conversion rate gamma of 0-0.9, sixteen data points are obtained for fuel deposition (alpha is 0.2-0.9), eighteen data points are obtained for high-temperature oxidation (gamma is 0.1-1), and the activation energy averaging treatment is carried out according to different oxidation stages, thus obtaining the low-temperature oxidation reaction (LTO),The activation energy of the fuel deposition reaction (FD) and the high temperature oxidation reaction (HTO) was 41.82KJ/mol, 143.7KJ/mol, and 81.3KJ/mol, respectively. In addition, based on the crude oil thermogravimetric experiment data at the temperature rising rate of 10 ℃/mol, the method comprises the following steps of fittingAnd a 1/T relation (formula 9) to obtain a fitting straight line of each reaction stage, and further determining that the pre-finger factors of different oxidation stages are 18.52/min, 1081357/min and 43052661/min respectively, as shown in figure 4.
The experimental crude oil sample is taken from a light oil reservoir, and compared with oxidation kinetic parameters which have been developed successfully by air injection flooding of crude oil, the experimental crude oil sample has an activation energy parameter of low-temperature oxidation of less than 63.3KJ/mol and a pre-finger factor of more than 12.09/min (the successfully implemented oxidation kinetic parameters of the air injection flooding oil reservoir), so that the experimental oil reservoir is indicated to be more suitable for air injection development. If the experimental crude oil sample is taken from a heavy oil reservoir, comparing oxidation kinetic parameters of the heavy oil reservoir which is successfully implemented and developed by injecting air and burning.
Compared with the prior art, the method for solving the parameters such as the accurate activation energy, the pre-finger factors and the like on the basis of optimizing the oxidation dynamics parameter model is provided, and the reference index can be provided for predicting the air injection flooding feasibility of the oil reservoir according to the parameters.
The foregoing details of the optional implementation manner of the embodiment of the present application have been described in detail with reference to the accompanying drawings, but the embodiment of the present application is not limited to the specific details of the foregoing implementation manner, and various simple modifications may be made to the technical solution of the embodiment of the present application within the scope of the technical concept of the embodiment of the present application, and these simple modifications all belong to the protection scope of the embodiment of the present application.
In addition, the specific features described in the foregoing embodiments may be combined in any suitable manner, and in order to avoid unnecessary repetition, various possible combinations are not described in detail in this application.
Moreover, any combination of the various embodiments of the present application may be made, so long as it does not deviate from the idea of the embodiment of the present application, and it should also be regarded as the disclosure of the embodiment of the present application.
Claims (10)
1. The method for judging the feasibility of the oil reservoir air injection flooding based on the oxidation kinetic parameters is characterized by comprising the following steps:
s10: dividing the crude oil oxidation reaction into three oxidation stages according to the distribution characteristics of oxidation characteristic curves of the crude oil at a plurality of different heating rates, wherein the three oxidation stages are respectively a low-temperature oxidation stage, a fuel deposition stage and a high-temperature oxidation stage;
s20: obtaining a plurality of characteristic values of the conversion rate of crude oil in three oxidation stages, determining the corresponding characteristic temperature of the characteristic conversion rate of a single oxidation stage in other temperature rise rate experiments, calculating the reaction rate value in the crude oil oxidation reaction process, and making a relation graph of ln (beta) changing along with 1/T, wherein T represents absolute temperature and is expressed as K; beta represents the temperature rising rate, and the unit is ℃/min;
s30: according to the oxidation characteristic curve, the activation energy corresponding to the three oxidation stages is obtained by combining the expression of the conversion rate method such as OFW and the like, and meanwhile, the pre-finger factors of crude oil in different oxidation stages are obtained according to an Arrhenius equation;
s40: and comparing the obtained activation energy and the pre-finger factor with corresponding preset data respectively to judge the feasibility of air injection driving of the oil reservoir.
2. The method for determining the feasibility of injecting air into a reservoir based on oxidation kinetic parameters according to claim 1, wherein said step S10 comprises the steps of:
s100: the crude oil oxidation reaction thermogravimetry experiment is repeated for a plurality of times, crude oil oxidation thermogravimetry TG and DTG experimental data under each heating rate are derived, and a crude oil TG curve graph and a crude oil DTG curve graph are respectively manufactured;
s110: fitting all the crude oil TG curve graphs and the crude oil DTG curve graphs to obtain an oxidation characteristic curve, and dividing the low-temperature oxidation stage, the fuel deposition stage and the high-temperature oxidation stage according to the turning points in the oxidation characteristic curve.
3. The method for discriminating the feasibility of injecting air into a reservoir based on oxidation kinetics according to claim 1 or 2 wherein three of said oxidation stages are divided by said oxidation profile of crude oil at least four different heating rates.
4. The method for determining the feasibility of injecting air into a reservoir based on oxidation kinetic parameters according to claim 1, wherein said step S20 comprises the steps of:
s200: dividing intervals of the conversion rate corresponding to a crude oil TG curve graph and a crude oil DTG curve of crude oil in the low-temperature oxidation stage, the fuel deposition stage and the high-temperature oxidation stage respectively according to the oxidation characteristic curve;
s210: and selecting the mass corresponding to the crude oil with the start-stop temperatures of the three oxidation stages in the full temperature range under one heating rate, so as to calculate the conversion rate gamma value at the reaction time t of the low-temperature oxidation stage, the fuel deposition stage and the high-temperature oxidation stage.
5. The method for determining the feasibility of injecting air into a reservoir based on oxidation kinetic parameters of claim 4, wherein said step S210 comprises the steps of:
s211: selecting the quality Z of crude oil with different start and stop temperatures in different oxidation stages in the full temperature range at one heating rate, wherein the quality Z is the quality of the crude oil in the original state in the low-temperature oxidation stage Li Quality Z of crude oil in ending state of low-temperature oxidation stage Lw Quality Z of raw state of crude oil fuel deposition stage Fi Quality Z of crude oil fuel deposition end state Fw Quality Z of crude oil in original state of high-temperature oxidation Hi Quality Z of crude oil in high temperature oxidation end state Hw ;
S212: will Z Li 、Z Lw 、Z Fi 、Z Fw 、Z Hi And Z Hw Substituting the conversion rate gamma of the low-temperature oxidation stage, the fuel deposition stage and the high-temperature oxidation stage at the reaction time t L Numerical value, gamma F Numerical value and gamma H In the numerical calculation formula:
for determining 20 characteristic values of the conversion gamma of crude oil in three oxidation stages from 0.1 to 1 in a gradient of 0.05, 60 characteristic values in total for the three oxidation stages, wherein Z Lt 、Z Ft And Z Ht And the crude oil quality corresponding to the reaction time t of the three oxidation stages is respectively obtained.
6. The method for determining the feasibility of injecting air into a reservoir based on oxidation kinetic parameters according to claim 5, wherein said step S210 further comprises the steps of:
s213: when the activation energy of the low-temperature oxidation stage is calculated, the conversion rate gamma is 0-0.9 for eighteen data points; when the activation energy of the fuel deposition stage is calculated, the conversion rate gamma is 0.2-0.9 for sixteen data points; when the activation energy of the high-temperature oxidation stage is calculated, the conversion rate gamma is 0.1-1 eighteen data points altogether.
7. The method for determining the feasibility of injecting air into an oil reservoir based on oxidation kinetic parameters according to claim 1, wherein the reaction rate value in the oxidation reaction of crude oil in step S20 is dγ/dT.
8. The method for determining the feasibility of injecting air into a reservoir based on oxidation kinetic parameters according to claim 1, wherein said step S30 comprises the steps of:
s300: according to the oxidation characteristic curve, fitting the slope of a straight line, substituting the slope into a solving formula of conversion rate methods such as OFW and the like, and solving the activation energy corresponding to each characteristic conversion rate;
s310: and carrying out activation energy averaging treatment on the activation energy corresponding to each characteristic conversion rate obtained in each oxidation stage, thereby obtaining the activation energy corresponding to the low-temperature oxidation stage, the fuel deposition stage and the high-temperature oxidation stage.
9. The method for judging the feasibility of injecting air into an oil reservoir based on oxidation kinetic parameters according to claim 1, wherein the expression of the conversion rate method such as OFW is as follows:
wherein G (gamma) represents a reaction mechanism integral function; beta represents the temperature rising rate, and the unit is ℃/min; a represents a factor before finger, in min -1 The method comprises the steps of carrying out a first treatment on the surface of the E represents reaction activation energy, and the unit is KJ/mol; r represents a molar gas constant of 8.314J/(mol.K); t represents absolute temperature in K.
10. The method for determining the feasibility of injecting air into a reservoir based on oxidation kinetic parameters according to claim 1, wherein the determining the pre-finger factors of the crude oil in different oxidation stages according to the arrhenius equation comprises:
based on crude oil thermogravimetric experiment data at a single heating rate, fittingWith 1/T relationship to obtainFitting straight lines to each reaction stage, and determining pre-finger factors of different oxidation stages by combining the Arrhenius equation through the intercept of a linear curve and a Y axis. />
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310118126.8A CN116258094B (en) | 2023-02-02 | 2023-02-02 | Oil reservoir air injection flooding feasibility judging method based on oxidation kinetic parameters |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310118126.8A CN116258094B (en) | 2023-02-02 | 2023-02-02 | Oil reservoir air injection flooding feasibility judging method based on oxidation kinetic parameters |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116258094A true CN116258094A (en) | 2023-06-13 |
CN116258094B CN116258094B (en) | 2024-02-06 |
Family
ID=86687548
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310118126.8A Active CN116258094B (en) | 2023-02-02 | 2023-02-02 | Oil reservoir air injection flooding feasibility judging method based on oxidation kinetic parameters |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116258094B (en) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5275719A (en) * | 1992-06-08 | 1994-01-04 | Mobil Oil Corporation | Production of high viscosity index lubricants |
CA2645101A1 (en) * | 1998-07-10 | 2000-01-20 | Stephen L. Buchwald | Ligands for metals and metal-catalyzed processes |
CN111610225A (en) * | 2019-02-25 | 2020-09-01 | 中国石油天然气股份有限公司 | Method for measuring oxidation exothermic property of crude oil |
CN113012763A (en) * | 2021-02-24 | 2021-06-22 | 西南石油大学 | Crude oil oxidation reaction kinetic model building method based on four-group components |
CN114970181A (en) * | 2022-06-06 | 2022-08-30 | 西南石油大学 | Shale oil reservoir air injection oxidation reaction kinetic model construction method |
-
2023
- 2023-02-02 CN CN202310118126.8A patent/CN116258094B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5275719A (en) * | 1992-06-08 | 1994-01-04 | Mobil Oil Corporation | Production of high viscosity index lubricants |
CA2645101A1 (en) * | 1998-07-10 | 2000-01-20 | Stephen L. Buchwald | Ligands for metals and metal-catalyzed processes |
CN111610225A (en) * | 2019-02-25 | 2020-09-01 | 中国石油天然气股份有限公司 | Method for measuring oxidation exothermic property of crude oil |
CN113012763A (en) * | 2021-02-24 | 2021-06-22 | 西南石油大学 | Crude oil oxidation reaction kinetic model building method based on four-group components |
CN114970181A (en) * | 2022-06-06 | 2022-08-30 | 西南石油大学 | Shale oil reservoir air injection oxidation reaction kinetic model construction method |
Also Published As
Publication number | Publication date |
---|---|
CN116258094B (en) | 2024-02-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Yang et al. | Effect of Temperature, Pressure and Equivalence Ratio on Ignition Delay in Ignition Quality Tester (IQT) Diesel, n-Heptane, and iso-Octane Fuels under Low Temperature Conditions | |
KR101430436B1 (en) | Method for determining fuel property and device for determining fuel property | |
CN100565197C (en) | A kind of measuring method of heat productivity of heating substance for calorimeter | |
CN116258094B (en) | Oil reservoir air injection flooding feasibility judging method based on oxidation kinetic parameters | |
WO2007001384A9 (en) | Flammability tester | |
CN110728039B (en) | Carbonized combustible unimodal pyrolysis kinetic parameter calculation method based on genetic algorithm | |
CN107421987A (en) | A kind of apparatus and method for measuring low-temperature oxidation of coal heat generation rate | |
JP2014005535A (en) | Method for evaluating solid fuel | |
EP2095105A1 (en) | A method and system for measuring ron and mon values for light distillates | |
Demuynck | A fuel independent heat transfer correlation for premixed spark ignition engines | |
CN110598185A (en) | Coke oven coking unit consumption influence factor sequencing method and system | |
CN113012763B (en) | Crude oil oxidation reaction kinetic model building method based on four-group components | |
CN108519301B (en) | Method for evaluating coal and biomass coke reactivity by using thermogravimetric analyzer | |
CN106569518A (en) | Temperature control method for volatile component measurement of coal of muffle furnace | |
CN112196514A (en) | Method for measuring deposition amount of fuel for thick oil air injection development by using thermogravimetric analyzer | |
CN207163948U (en) | A kind of device for measuring low-temperature oxidation of coal heat generation rate | |
EP2685258B1 (en) | Fuel physical property determination method and device | |
CN113670980A (en) | Method for determining shortest natural ignition period of coal | |
Bussmann et al. | Model predictions of temperature and velocity profiles in turbulent diffusion buoyant flames | |
Al-Hamamre et al. | Thermal partial oxidation of diesel in porous reactors for synthesis gas production | |
TWI818869B (en) | Method and device for evaluating an optimum transverse furnace temperature of a coke oven | |
CN109253943A (en) | The monitoring method of oil nature is mixed in a kind of coal-based needle coke production process | |
Hartl | Flamelet/progress variable modelling and flame structure analysis of partially premixed flames | |
KR102291910B1 (en) | Test equipment and test methd | |
JP2021156258A (en) | Model construction apparatus, model construction method, engine control parameter determination method, and engine |
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 |