CN112908425A - Microbial oil recovery numerical simulation method based on reaction kinetics model - Google Patents

Microbial oil recovery numerical simulation method based on reaction kinetics model Download PDF

Info

Publication number
CN112908425A
CN112908425A CN202110124181.9A CN202110124181A CN112908425A CN 112908425 A CN112908425 A CN 112908425A CN 202110124181 A CN202110124181 A CN 202110124181A CN 112908425 A CN112908425 A CN 112908425A
Authority
CN
China
Prior art keywords
reaction
microbial
microorganism
oil
kinetic
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.)
Withdrawn
Application number
CN202110124181.9A
Other languages
Chinese (zh)
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 East China
Original Assignee
China University of Petroleum 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 East China filed Critical China University of Petroleum East China
Priority to CN202110124181.9A priority Critical patent/CN112908425A/en
Publication of CN112908425A publication Critical patent/CN112908425A/en
Withdrawn legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
    • G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
    • G16C20/10—Analysis or design of chemical reactions, syntheses or processes
    • 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/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/10—Numerical modelling

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Chemical & Material Sciences (AREA)
  • Computing Systems (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Crystallography & Structural Chemistry (AREA)
  • Chemical Kinetics & Catalysis (AREA)
  • Analytical Chemistry (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Preparation Of Compounds By Using Micro-Organisms (AREA)

Abstract

The invention relates to a microbial oil recovery numerical simulation method based on a reaction kinetics model, which specifically comprises the following steps: step S1: determining a microorganism reaction kinetic component; step S2: establishing a microbial reaction kinetic equation; step S3: determining a microbial reaction kinetic parameter; step S4: establishing a microorganism reaction kinetic model; and step S5, creating a geological conceptual model based on CMG software, and verifying the accuracy of the microbial reaction kinetic model. The invention utilizes a microorganism reaction kinetic equation to quantitatively describe the growth, reproduction, metabolism and death characteristics of microorganisms, combines mature oil reservoir numerical simulation software to carry out functional development on microorganism oil recovery numerical simulation, forms a set of practical microorganism oil recovery numerical simulation method, and provides an important reference basis for the optimized design of a microorganism oil-displacement scheme on an oil field site.

Description

Microbial oil recovery numerical simulation method based on reaction kinetics model
Technical Field
The invention relates to the field of microbial flooding for improving recovery efficiency, in particular to a microbial oil recovery numerical simulation method based on a reaction kinetics model.
Background
The microbe flooding technology for raising recovery ratio is characterized by that it utilizes the beneficial metabolic product and beneficial metabolic activity of microbe to change the interface property of crude oil-rock-water, and makes it react with residual oil in the oil reservoir to raise oil-water fluidity ratio and increase the flow capacity of crude oil in the pores of oil reservoir so as to attain the goal of reducing water and raising oil recovery ratio. Compared with other tertiary oil recovery technologies, the microbial enhanced oil recovery technology has the characteristics of low reservoir damage, no environmental pollution, simple field operation, strong oil reservoir adaptability and low raw material price.
At the end of the 20 th century, Islam and Chang established mathematical models of microbial oil recovery based on black oil models. Since then, the research of the numerical simulation theory of the microbial oil recovery is continuously developed at home and abroad, and the more classical models comprise an Islam model, a Chang model, a microbial oil recovery mathematical model established by Raynaud and a microbial oil recovery mathematical model established by Kjew. However, the description of the microbial oil recovery mechanism is mostly a theoretical description of the microbes and their metabolites so far, and a simulation method and a parameter determination method combined with mature reservoir numerical simulation software and an optimized design of the embodiment of the microbial oil recovery field are lacked.
Disclosure of Invention
The invention aims to provide a microbial oil recovery simulation method and a parameter determination method which are combined with mature oil reservoir numerical simulation software, so as to realize the optimization design of a microbial oil recovery field implementation scheme.
In order to achieve the purpose, the invention adopts the following technical scheme.
A microbial oil recovery numerical simulation method based on a reaction kinetics model specifically comprises the following steps:
step S1: determination of the kinetic components of the microbial reaction:
the microorganism reaction kinetic components comprise microorganisms, oxygen, nutrient substances, biosurfactant, biopolymer, methane and water;
determination of the microbial composition: experiments and microbial component element analysis prove that the composition of microbial cell elements for oil displacement is stable, and the chemical formula is represented as C4-7H7-10O1-3And N has a relative molecular mass of 85-156.
Determination of nutrient composition: the nutrient substance comprises carbon source and nitrogen source, and the carbon source comprises blend oil (C)18-60H30- 116O2-15) Glucose (C)6H12O6) Fructose (C)6H12O6) Sucrose (C)12H22O11) Maltose (C)12H22O11·H2O) and starch ((C)5H10O5)n) Composition, there is no fixed molecular formula for complex organic carbon sources. The nitrogen source is divided into inorganic nitrogen source and organic nitrogen source, and NH is used as inorganic nitrogen source3Instead, the organic nitrogen source consists of amino acids.
Determination of biosurfactant component: the surfactant for displacing oil mainly refers to rhamnolipid, sophoricoside and mannitol ester, and its molecular formula is C16-32H30-58O7-13The relative molecular mass was 334-650.
Determination of the biopolymer component: the biosurfactant for displacing oil is xanthan gum (monomer chemical formula C)67H99O56) Scleroglucan (monomer formula C)24H40O20) And welanGlue (monomer formula C)30H48O24) Assuming that the biopolymer is composed of a mixture of the above-mentioned substances, its chemical formula is C24-67H40-99O20-56The relative molecular mass is 648-1799.
Aiming at the field application of microbial oil recovery, an element analysis instrument can be used for analyzing the element composition of each component, so that the chemical formula and the relative molecular mass of the component can be accurately obtained.
Step S2: and establishing a microorganism reaction kinetic equation.
From the components of the microbial reaction kinetics determined in step S1, the following microbial reaction kinetics equation is established.
The microorganism growth metabolic reaction kinetic equation:
Figure BDA0002923353180000021
Figure BDA0002923353180000022
Figure BDA0002923353180000023
Figure BDA0002923353180000031
in the formula: ciHmOnNfIs a microorganism; cx1Hy1Oz1Is a carbon source; hy2Oz2Nj2Is a nitrogen source; cu1Hv1Ow1Is a biosurfactant; cu2Hv2Ow2Is a biopolymer; subscripts i, m, n, f, x1、y1、z1、y2、z2、j2、u1、v1、w1、u2、v2、w2Is one by oneAtomic number of the element; y isx/cCell productivity coefficients based on carbon; y isp/c' is the coefficient of biosurfactant yield based on carbon element; y isp/c"is the biopolymer yield coefficient based on carbon element; a. b and c are coefficients to be solved.
Sixthly, a kinetic equation of polymer degradation reaction:
Figure BDA0002923353180000032
in the formula: cu2Hv2Ow2Is a biopolymer; u. of2、v2、w2Is the atomic number of C, H, O elements.
Seventhly, a kinetic equation of the microorganism death reaction:
Figure BDA0002923353180000033
in the formula: ciHmOnNfIs a microorganism; i. m, n and f are the atomic number of C, H, O, N elements.
The reaction kinetics equation of the microbial degradation of the crude oil is as follows:
CiHmOnNf+aX1＝CiHmOnNf+bX2
here, because the crude oil components are complex and variable, specific chemical elements are not used to characterize the elemental compositions, and the reaction equation is balanced according to mass conservation assuming that no change occurs before and after the crude oil is degraded by microorganisms.
In the formula: ciHmOnNfIs a microorganism; x1Is a high-viscosity oil phase; x2Is a low-viscosity oil phase; a. b is the reaction coefficient of high viscosity oil phase and low viscosity oil phase.
Step S3: determining a microbial reaction kinetic parameter.
The microorganism reaction kinetic parameters comprise reaction order n and reaction frequency factorSub k0And activation energy Ea。
Determining the kinetic parameters of the growth and metabolic reaction of the microorganisms:
Figure BDA0002923353180000041
according to the Moro equation, when the concentration of the nutrient is low, the growth, reproduction and metabolism of the microorganism are restricted by the concentration of the nutrient and follow the first-order reaction rule, so that the reaction order n is 1.
Reaction frequency factor k0The arrhenius equation was used to determine:
Figure BDA0002923353180000042
Figure BDA0002923353180000043
in the formula: k is a reaction rate constant; k is a radical of0Is a reaction frequency factor; eakJ. mol for activation energy-1(ii) a R is a molar gas constant of 8.314 J.mol-1·K-1(ii) a T is absolute temperature, K; t is t1/2Is the nutrient concentration half-life, d.
During the growth and metabolism of the microorganism, the influence of temperature change on the reaction rate is not considered, so the activation energy Ea＝0。
Sixthly, determining the kinetic parameters of the polymer degradation reaction:
Figure BDA0002923353180000044
the degradation reaction of the biological polymer is a first-order reaction, the reaction order n is 1, and the reaction frequency factor k0K, activation energy Ea＝0。
Figure BDA0002923353180000045
In the formula: k is a reaction rate constant, d-1；t1/2Is the polymer half-life, d.
Determining kinetic parameters of the death reaction of the microorganisms:
Figure BDA0002923353180000046
the microbial death reaction is a first-order reaction, the reaction order n is 1, and the reaction frequency factor k is0K, activation energy Ea＝0。
Figure BDA0002923353180000051
In the formula: k is a reaction rate constant, d-1；t1/2Is the microbial half-life, d.
Determining reaction kinetic parameters of the microbial degradation crude oil:
CiHmOnNf+aX1＝CiHmOnNf+bX2
the reaction of the microbial degradation of the crude oil conforms to a first-order reaction rule, the influence of temperature change on the reaction rate is not considered, and the relationship between the reaction rate and the concentration of the high-viscosity crude oil is as follows:
Figure BDA0002923353180000052
the time integration is given by:
Figure BDA0002923353180000053
in the formula: c is the concentration of high viscous crude oil in the culture solution at any moment, g/L; c0The concentration of high viscous crude oil in the culture solution at the initial moment is g/L; k is a reaction rate constant, d-1。
Here, the reaction kinetics parameters for the microbial degradation of crude oil were determined as: reaction order n is 1, reaction frequency factor k0K, activation energy Ea＝0。
Step S4: and establishing a microorganism reaction kinetic model.
In CMG numerical simulation software, microorganisms and metabolites thereof obtained by a microorganism autocatalytic reaction equation are respectively defined as components dissolved in water, and the migration of the microorganisms and the metabolites thereof can be simulated and calculated by considering characteristics such as growth, migration, diffusion, adsorption and the like (a CMG default value or an experimental measurement value is adopted in a model).
Equation of microbial migration:
Figure BDA0002923353180000054
nutrient transport equation:
Figure BDA0002923353180000061
metabolite transport equation:
Figure BDA0002923353180000062
and step S5, creating a geological conceptual model based on CMG software, and verifying the accuracy of the microbial reaction kinetic model.
And (3) creating a geological concept model through a Builder module in the CMG, creating a microbial reaction kinetic model by utilizing a thermal recovery and chemical flooding simulator STARS, performing a microbial oil recovery numerical simulation experiment, and verifying the accuracy of the microbial reaction kinetic model.
The invention has the beneficial effects that: the invention utilizes a microorganism reaction kinetic equation to quantitatively describe the growth, reproduction, metabolism and death characteristics of microorganisms, combines mature oil reservoir numerical simulation software to carry out functional development on microorganism oil recovery numerical simulation, forms a set of practical microorganism oil recovery numerical simulation method, and provides an important reference basis for the optimized design of a microorganism oil-displacement scheme on an oil field site.
Drawings
FIG. 1 is a flow chart of a microbial oil recovery simulation method adopted by the invention.
FIG. 2 is a schematic representation of the Mornor equation employed in the practice of the present invention.
FIG. 3 is a graph of nutrient consumption used in a method embodying the present invention.
FIG. 4 is a schematic representation of the half-life of a polymer employed in an embodiment of the present invention.
FIG. 5 is a graph showing the growth and death of microorganisms used in a method of practicing the present invention.
FIG. 6 is a graph of the degradation curve for highly viscous crude oil used in an embodiment of the present invention.
FIG. 7 is a geological model map used in the practice of the invention.
FIG. 8 is a graph of relative permeability used in a method embodying the present invention.
FIG. 9 is a biopolymer viscosity concentration curve used in a method embodying the present invention.
FIG. 10 is a graph of the concentration of a biosurfactant versus interfacial tension used in a method embodying the present invention.
FIG. 11 is a graph comparing oil production in a method embodying the present invention.
FIG. 12 is a graph showing a comparison of water content in the method of the present invention.
Detailed description of the invention
The invention is further described with reference to the accompanying drawings and specific embodiments.
Examples
The microbial oil recovery numerical simulation method in the embodiment of the invention is realized by adopting the process shown in FIG. 1, and comprises the following specific steps:
and step S1, determining the components of the reaction kinetics of the microorganisms.
The microbial reaction kinetics components mainly comprise microorganisms, oxygen, nutrients, biosurfactants, biopolymers, methane and water.
Determination of the microbial composition: experiments and microbial component element analysis prove that the composition of microbial cell elements for oil displacement is stable, and the chemical formula can be expressed as C4-7H7-10O1-3N, relative molecular mass 85-156, in this example the microorganism chemical formula is C4H10O3N, relative molecular mass 120.
Determination of nutrient composition: the nutrient substance comprises carbon source and nitrogen source, and the carbon source comprises blend oil (C)18-60H30- 116O2-15) Glucose (C)6H12O6) Fructose (C)6H12O6) Sucrose (C)12H22O11) Maltose (C)12H22O11·H2O) and starch ((C)5H10O5)n) Composition, no fixed molecular formula for complex organic carbon source, in this example model blend oil was used as carbon source, molecular formula C18H34O2The relative molecular mass was 282. The nitrogen source is divided into inorganic nitrogen source and organic nitrogen source, and NH is used as inorganic nitrogen source3Instead, the organic nitrogen source consists of amino acids, in this example ammonia, of formula NH, is used as nitrogen source3·H2O, relative molecular mass 35.
Determination of biosurfactant component: the surfactant for displacing oil mainly refers to rhamnolipid, sophoricoside and mannitol ester, and its molecular formula is C16-32H30-58O7-13The relative molecular mass is 334-650, in this example, the chemical formula of the biosurfactant is C18H51O10The relative molecular mass was 427.
Determination of the biopolymer component: the biosurfactant for displacing oil is xanthan gum (monomer chemical formula C)67H99O56) Scleroglucan (monomer formula C)24H40O20) And welan gum (monomer formula C)30H48O24) Assuming that the biopolymer is composed of a mixture of the above-mentioned substances, its chemical formula is C24-67H40-99O20-56The relative molecular mass is 648-1799, in this example, the chemical formula of the biopolymer is C27H76O25The relative molecular mass was 800.
Aiming at the field application of microbial oil recovery, an element analysis instrument can be used for analyzing the element composition of each component, so that the chemical formula and the relative molecular mass of the component can be accurately obtained.
And step S2, establishing a microbial reaction kinetic equation.
From the components of the microbial reaction kinetics determined in step S1, the following microbial reaction kinetics equation is established.
The kinetic equation of the growth and metabolic reaction of the microorganisms is as follows:
Figure BDA0002923353180000081
a＝4.5Yx/c
b＝16.875Yx/c+10.25Yp/c′+13.67Yp/c″-10.5
c＝24.75Yx/c+10.5Yp/c′+10.67Yp/c″-19
in the formula: c4H10O3N is a microorganism; c18H34O2Is a carbon source; h5ON is a nitrogen source; c18H51O10Is a biosurfactant; c27H76O25Is a biopolymer; y isx/cCell productivity coefficients based on carbon; y isp/c' is the coefficient of biosurfactant yield based on carbon element; y isp/c"is the biopolymer yield coefficient based on carbon element.
The cell productivity coefficient Y was taken for the model of this examplex/c0.1, biosurfactant yield coefficient Yp/c' -0.3, biopolymer yield coefficient Yp/cAnd 0.5, the kinetic equation of the growth and reproduction reaction of the microorganisms is as follows:
C4H10O3N+C18H34O2+0.45H5ON+1.098O2+8.04H2O＝1.45C4H10O3N+0.3C18H51O10+0.333C27H76O25+18(1-Yx/c-Yp/c′-Yp/c″)CH4
kinetic equation of polymer degradation reaction:
C27H76O25＝44.444H2O
in the formula: c27H76O25Is a biopolymer.
③ kinetic equation of death reaction of microorganisms:
C4H10O3N＝6.667H2O
in the formula: c4H10O3N is a microorganism.
Fourthly, the reaction kinetic equation of the microbial degradation of the crude oil is as follows:
C4H10O3N+5X1＝C4H10O3N+6X2
here, because the crude oil composition is complicated and varied, no specific chemical element is used to characterize the elemental composition, and assuming no change before and after the microbial degradation of the crude oil, the reaction equation is balanced according to the mass conservation, i.e., the molecular weight X of the high-viscosity crude oil in the model of this example1600, low viscosity crude oil molecular weight X2Is 500.
In the formula: c4H10O3N is a microorganism; x1Is a high-viscosity oil phase; x2Is a low viscosity oil phase.
Step S3, determining the microbial reaction kinetic parameters.
The microorganism reaction kinetic parameters comprise reaction order n and reaction frequency factor k0And activation energy Ea。
Determination of the kinetic parameters of the growth and metabolic reaction of the microorganisms:
C4H10O3N+C18H34O2+0.45H5ON+1.098O2+8.04H2O＝1.45C4H10O3N+0.3C18H51O10+0.333C27H76O25+18(1-Yx/c-Yp/c′-Yp/c″)CH4
according to the morno equation, when the concentration of nutrients is low, the growth, reproduction and metabolism of microorganisms are restricted by the concentration of nutrients, and follow the first-order reaction law, so that the reaction order number n is 1, see fig. 2.
In this example model, the effect of temperature change on the reaction rate was not taken into account, so the activation energy Ea＝0。
Reaction frequency factor k0The arrhenius equation was used to determine:
Figure BDA0002923353180000091
Figure BDA0002923353180000092
in the formula: k is a reaction rate constant; k is a radical of0Is a reaction frequency factor; eakJ. mol for activation energy-1(ii) a R is a molar gas constant of 8.314 J.mol-1·K-1(ii) a T is absolute temperature, K; t is t1/2Is the nutrient concentration half-life, d.
Here, because E a0, so there is:
Figure BDA0002923353180000101
in this example, the nutrient consumption curve is shown in FIG. 3, the half-life t of nutrient concentration1/20.55d, reaction frequency factor k0＝1.26。
For the field application of microbial oil recovery, the nutrient consumption curve of a certain microorganism can be measured, so that the reaction frequency factor of the microorganism can be determined.
Determining kinetic parameters of polymer degradation reaction:
C27H76O25＝44.444H2O
the degradation reaction of the biological polymer is a first-order reaction, the reaction order n is 1, and the reaction frequency factor k0K, activation energy Ea＝0。
Figure BDA0002923353180000102
In the formula: k is a reaction rate constant, d-1；t1/2Is the polymer half-life, d.
Schematic of the biopolymer half-Life referring to FIG. 4, in this example model, the biopolymer half-life t1/290d, the frequency factor k is reflected0＝0.0077。
For microbial oil recovery field applications, the half-life of a certain biopolymer can be measured to determine its response frequency factor.
Determination of kinetic parameters of the death reaction of the microorganisms:
C4H10O3N＝6.667H2O
the microbial death reaction is a first-order reaction, the reaction order n is 1, and the reaction frequency factor k is0K, activation energy Ea＝0。
Figure BDA0002923353180000111
In the formula: k is a reaction rate constant, d-1；t1/2Is the microbial half-life, d.
Microorganism growth Curve schematic referring to FIG. 5, the half-life of the microorganism can be determined from the decay phase, in this example the half-life t of the microorganism1/260d, then the frequency factor k is reflected0＝0.0116。
For microbial oil recovery field applications, the half-life of a certain microorganism can be measured to determine its response frequency factor.
Fourthly, determining the reaction kinetic parameters of the microbial degradation of the crude oil:
C4H10O3N+5X1＝C4H10O3N+6X2
the reaction for degrading crude oil by microbe accords with the first-order reaction rule, the reaction order n is 1, the influence of temperature change on the reaction rate is not considered, and the activation energy EaWhen the reaction rate is 0, the relationship between the reaction rate and the concentration of the high-viscosity crude oil is as follows:
Figure BDA0002923353180000112
the time integration is given by:
Figure BDA0002923353180000113
in the formula: c is the concentration of high viscous crude oil in the culture solution at any moment, g/L; c0The concentration of high viscous crude oil in the culture solution at the initial moment is g/L; k is a reaction rate constant, d-1。
The curve of microbially degraded crude oil, see FIG. 6, can be used to determine the reaction frequency factor, k, according to the slope, in this example model0＝0.1623。
Aiming at the field application of microbial oil recovery, the reaction frequency factor can be determined by measuring the concentration change of high viscous crude oil in the culture solution.
And step S4, establishing a microorganism reaction kinetic model.
In CMG numerical simulation software, microorganisms and metabolites thereof obtained by a microorganism autocatalytic reaction equation are respectively defined as components dissolved in water, and the migration of the microorganisms and the metabolites thereof can be simulated and calculated by considering characteristics such as growth, migration, diffusion, adsorption and the like (a CMG default value or an experimental measurement value is adopted in a model).
Equation of microbial migration:
Figure BDA0002923353180000121
nutrient transport equation:
Figure BDA0002923353180000122
metabolite transport equation:
Figure BDA0002923353180000123
and step S5, creating a conceptual model based on CMG software, and optimally designing the microbial oil recovery scheme.
Firstly, a CMG-Builder module is utilized to establish a geological concept model by a one-injection-one-sampling five-point method, a grid system of 15 multiplied by 10 is adopted, a geological concept model diagram is shown in figure 7, relevant reservoir and fluid parameters are input, parameters adopted in the model are shown in table 1, and a relative permeability curve is shown in figure 8.
TABLE 1 reservoir and fluid parameters table
Figure BDA0002923353180000124
Figure BDA0002923353180000131
② introducing a biopolymer-surfactant oil displacement system by a Process Wizard Process guide, and the polymer viscosity-concentration relation curve and the surfactant interfacial tension-concentration relation curve adopted in the model are shown in figure 9 and figure 10.
And thirdly, introducing a kinetic equation of microorganism growth and metabolism, crude oil degradation, polymer decomposition and microorganism death Reaction through a CMG-Reaction module.
The model adopts a displacement mode of converting water drive to water drive from microbial drive with water content of 90 percent, and analyzes main factors influencing the effect of improving the recovery ratio by the microbial drive through an orthogonal design method to optimally design the microbial oil recovery scheme. The orthogonal design parameters used in the model of this example are shown in Table 2, the microbial displacement parameters are shown in Table 3, and the results of the orthogonal experiments are shown in Table 4.
TABLE 2 microbial flooding orthogonal design Table
Figure BDA0002923353180000132
TABLE 3 microbial displacement protocol parameter Table
Figure BDA0002923353180000133
Figure BDA0002923353180000141
TABLE 4 orthogonal test results Table
Injection scheme Cumulative oil production (m)3) Dimensionless oil increasing Recovery ratio (%) Enhanced recovery (%)
Water drive 24446.2 / 52.23 /
Scheme 1 24734.8 0.12 52.85 0.62
Scheme 2 25218.6 0.31 53.88 1.65
Scheme 3 25433.7 0.40 54.34 2.11
Scheme 4 25723.4 0.52 54.96 2.73
Scheme 5 24979.4 0.22 53.37 1.14
Scheme 6 25150.3 0.29 53.74 1.51
Scheme 7 26436 0.81 56.48 4.25
Scheme 8 26112.4 0.68 55.79 3.56
Scheme 9 24888.3 0.18 53.18 0.95
Scheme 10 26913.3 1.00 57.50 5.27
Scheme 11 26648.7 0.89 56.94 4.71
Scheme 12 26211.6 0.72 56.00 3.77
Scheme 13 25954.2 0.61 55.45 3.22
Scheme 14 26580 0.86 56.79 4.56
Scheme 15 26524.8 0.84 56.67 4.44
Scheme 16 25783.1 0.54 55.09 2.86
The optimal injection scheme 10 is taken as an example to be compared with a water flooding scheme, the cumulative oil production comparison graph is shown in figure 11, and the water content comparison graph is shown in figure 12.
Based on the optimal injection scheme 10, under the condition of ensuring that the total injection amount and the total gas-liquid ratio of nutrients are not changed, three injection modes of periodic continuous injection, half-period continuous injection and high-concentration short-time concentrated injection are simulated, and the influence of the three injection modes on the improvement of the recovery ratio is analyzed.
(ii) periodic continuous injection
On the premise of ensuring the total injection amount of the activator and the air, the activator and the air are continuously injected in a low-concentration mode.
Half cycle continuous injection
On the premise of ensuring the total injection amount of the activator and the air, the activator is alternately injected for 15 days and the water and the air are injected for 15 days in a higher concentration mode.
High concentration short time centralized injection
On the premise of ensuring the total injection amount of the activator and the air, the activator is alternately injected for 1 day and the water and the air are injected for 29 days in a high-concentration mode.
The different injection enhanced recovery results are shown in table 5.
TABLE 5 injection results Table
Injection mode Accumulated oil (m)3) Cumulative oil (m)3) Recovery ratio (%) Enhanced recovery (%)
Water drive 24446.2 / 52.23 /
Periodic continuous injection 26913.3 2467.1 57.50 5.27
Half cycle continuous injection 26906.4 2460.2 47.49 5.26
High concentration short time concentrated injection 26878.8 2432.6 57.43 5.20
It can be seen that the effect of improving the recovery ratio of the periodic continuous injection, the half-period continuous injection and the high-concentration short-time concentrated injection is not greatly different, the development effect of the high-concentration short-time concentrated injection is slightly lower than that of the periodic continuous injection and the half-period continuous injection, but the high-concentration short-time concentrated injection is short in time, convenient to operate and more reasonable from the economic perspective. Therefore, it is recommended to use a high concentration short time concentration implantation method.

Claims (1)

1. A microbial oil recovery numerical simulation method based on a reaction kinetics model is characterized in that: the method specifically comprises the following steps:
step S1: determination of the kinetic components of the microbial reaction: the microorganism reaction kinetic components comprise microorganisms, oxygen, nutrient substances, biosurfactant, biopolymer, methane and water; determination of the microbial composition: experiments and microbial component element analysis prove that the composition of microbial cell elements for oil displacement is stable, and the chemical formula is represented as C4-7H7-10O1-3N, the relative molecular mass is 85-156; determination of nutrient composition: the nutrient substance comprises carbon source and nitrogen source, and the carbon source comprises blend oil (C)18-60H30- 116O2-15) Glucose (C)6H12O6) Fructose, fructose(C6H12O6) Sucrose (C)12H22O11) Maltose (C)12H22O11·H2O) and starch ((C)5H10O5)n) Composition, there is no fixed molecular formula for complex organic carbon sources; the nitrogen source is divided into inorganic nitrogen source and organic nitrogen source, and NH is used as inorganic nitrogen source3Instead, the organic nitrogen source consists of amino acids; determination of biosurfactant component: the surfactant for displacing oil mainly refers to rhamnolipid, sophoricoside and mannitol ester, and its molecular formula is C16-32H30-58O7-13Relative molecular mass 334-650; determination of the biopolymer component: the biosurfactant for displacing oil is xanthan gum (monomer chemical formula C)67H99O56) Scleroglucan (monomer formula C)24H40O20) And welan gum (monomer formula C)30H48O24) Assuming that the biopolymer is composed of a mixture of the above-mentioned substances, its chemical formula is C24-67H40-99O20-56The relative molecular mass is 648-1799; aiming at the field application of microbial oil recovery, an element analysis instrument can be used for analyzing the element composition of each component, so that the chemical formula and the relative molecular mass of the component can be accurately obtained;
step S2: establishing a microbial reaction kinetic equation: establishing a following microorganism reaction kinetic equation according to the microorganism reaction kinetic components determined in the step S1;
the kinetic equation of the growth and metabolic reaction of the microorganisms is as follows:
Figure FDA0002923353170000011
Figure FDA0002923353170000012
Figure FDA0002923353170000013
Figure FDA0002923353170000021
in the formula: ciHmOnNfIs a microorganism; cx1Hy1Oz1Is a carbon source; hy2Oz2Nj2Is a nitrogen source; cu1Hv1Ow1Is a biosurfactant; cu2Hv2Ow2Is a biopolymer; subscripts i, m, n, f, x1、y1、z1、y2、z2、j2、u1、v1、w1、u2、v2、w2Is the atomic number of each element; y isx/cCell productivity coefficients based on carbon; y isp/c' is the coefficient of biosurfactant yield based on carbon element; y isp/c"is the biopolymer yield coefficient based on carbon element; a. b and c are coefficients to be solved;
kinetic equation of polymer degradation reaction:
Figure FDA0002923353170000022
in the formula: cu2Hv2Ow2Is a biopolymer; u. of2、v2、w2Is the atomic number of C, H, O element;
③ kinetic equation of death reaction of microorganisms:
Figure FDA0002923353170000023
in the formula: ciHmOnNfIs a microorganism; i. m, n and f are the atomic number of C, H, O, N elements;
fourthly, the reaction kinetic equation of the microbial degradation of the crude oil is as follows:
CiHmOnNf+aX1＝CiHmOnNf+bX2
here, because the crude oil components are complex and changeable, specific chemical elements are not used for representing the element compositions, and if no change occurs before and after the crude oil is degraded by microorganisms, the reaction equation is balanced according to mass conservation;
in the formula: ciHmOnNfIs a microorganism; x1Is a high-viscosity oil phase; x2Is a low-viscosity oil phase; a. b is the reaction coefficient of high viscosity oil phase and low viscosity oil phase;
step S3: determining the kinetic parameters of the microbial reaction: the microorganism reaction kinetic parameters comprise reaction order n and reaction frequency factor k0And activation energy Ea；
Determination of the kinetic parameters of the growth and metabolic reaction of the microorganisms:
Figure FDA0002923353170000024
according to the Moro equation, when the concentration of the nutrient is low, the growth, reproduction and metabolism of the microorganism are restricted by the concentration of the nutrient and follow the first-order reaction rule, so that the reaction order n is 1;
reaction frequency factor k0The arrhenius equation was used to determine:
Figure FDA0002923353170000031
Figure FDA0002923353170000032
in the formula: k is a reaction rate constant; k is a radical of0Is a reaction frequency factor; eakJ. mol for activation energy-1(ii) a R is a molar gas constant of 8.314 J.mol-1·K-1(ii) a T is absolute temperature, K; t is t1/2Half life of nutrient concentration, d;
during the growth and metabolism of the microorganism, the influence of temperature change on the reaction rate is not considered, so the activation energy Ea＝0；
Determining kinetic parameters of polymer degradation reaction:
Figure FDA0002923353170000033
the degradation reaction of the biological polymer is a first-order reaction, the reaction order n is 1, and the reaction frequency factor k0K, activation energy Ea＝0；
Figure FDA0002923353170000034
In the formula: k is a reaction rate constant, d-1；t1/2Is the polymer half-life, d;
determination of kinetic parameters of the death reaction of the microorganisms:
Figure FDA0002923353170000035
the microbial death reaction is a first-order reaction, the reaction order n is 1, and the reaction frequency factor k is0K, activation energy Ea＝0；
Figure FDA0002923353170000041
In the formula: k is a reaction rate constant, d-1；t1/2Is the microbial half-life, d;
fourthly, determining the reaction kinetic parameters of the microbial degradation of the crude oil:
CiHmOnNf+aX1＝CiHmOnNf+bX2
the reaction of the microbial degradation of the crude oil conforms to a first-order reaction rule, the influence of temperature change on the reaction rate is not considered, and the relationship between the reaction rate and the concentration of the high-viscosity crude oil is as follows:
Figure FDA0002923353170000042
the time integration is given by:
Figure FDA0002923353170000043
in the formula: c is the concentration of high viscous crude oil in the culture solution at any moment, g/L; c0The concentration of high viscous crude oil in the culture solution at the initial moment is g/L; k is a reaction rate constant, d-1；
Here, the reaction kinetics parameters for the microbial degradation of crude oil were determined as: reaction order n is 1, reaction frequency factor k0K, activation energy Ea＝0；
Step S4: establishing a microbial reaction kinetic model: in CMG numerical simulation software, respectively defining microorganisms and metabolites thereof obtained by a microorganism autocatalytic reaction equation as water-soluble components, and considering characteristics such as growth, migration, diffusion, adsorption and the like (a CMG default value or an experimental measurement value is adopted in a model), so that the simulation calculation of the migration of the microorganisms and the metabolites thereof can be realized;
equation of microbial migration:
Figure FDA0002923353170000044
nutrient transport equation:
Figure FDA0002923353170000051
metabolite transport equation:
Figure FDA0002923353170000052
step S5, creating a geological conceptual model based on CMG software, and verifying the accuracy of the microbial reaction kinetic model; and (3) creating a geological concept model through a Builder module in the CMG, creating a microbial reaction kinetic model by utilizing a thermal recovery and chemical flooding simulator STARS, performing a microbial oil recovery numerical simulation experiment, and verifying the accuracy of the microbial reaction kinetic model.
CN202110124181.9A 2021-01-29 2021-01-29 Microbial oil recovery numerical simulation method based on reaction kinetics model Withdrawn CN112908425A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110124181.9A CN112908425A (en) 2021-01-29 2021-01-29 Microbial oil recovery numerical simulation method based on reaction kinetics model

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110124181.9A CN112908425A (en) 2021-01-29 2021-01-29 Microbial oil recovery numerical simulation method based on reaction kinetics model

Publications (1)

Publication Number Publication Date
CN112908425A true CN112908425A (en) 2021-06-04

Family

ID=76120772

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110124181.9A Withdrawn CN112908425A (en) 2021-01-29 2021-01-29 Microbial oil recovery numerical simulation method based on reaction kinetics model

Country Status (1)

Country Link
CN (1) CN112908425A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117198385A (en) * 2023-11-07 2023-12-08 陕西延长石油（集团）有限责任公司 Oil reservoir microbial community competition evolution method

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103104238A (en) * 2013-01-16 2013-05-15 中国石油大学(华东) Microorganism oil displacement numerical simulation method
CN110593833A (en) * 2019-09-04 2019-12-20 中国石油化工股份有限公司 Simulation method for improving recovery ratio by microbial flooding
CN111222253A (en) * 2020-01-20 2020-06-02 中国科学院生态环境研究中心 CFD-ASM coupling simulation model, establishment method thereof and sewage treatment system optimization method
CN111985735A (en) * 2019-05-05 2020-11-24 中国石油天然气股份有限公司 Endogenous microbial oil displacement oil reservoir screening and evaluating method

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103104238A (en) * 2013-01-16 2013-05-15 中国石油大学(华东) Microorganism oil displacement numerical simulation method
CN111985735A (en) * 2019-05-05 2020-11-24 中国石油天然气股份有限公司 Endogenous microbial oil displacement oil reservoir screening and evaluating method
CN110593833A (en) * 2019-09-04 2019-12-20 中国石油化工股份有限公司 Simulation method for improving recovery ratio by microbial flooding
CN111222253A (en) * 2020-01-20 2020-06-02 中国科学院生态环境研究中心 CFD-ASM coupling simulation model, establishment method thereof and sewage treatment system optimization method

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
YAO CHUANJIN等: "Numerical Simulation On Oil Shale In-Situ Upgrading By Steam Injection", 《IFEDC》 *
王昌伟: "微生物驱油数学模型及数值模拟研究", 《硕士电子期刊》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117198385A (en) * 2023-11-07 2023-12-08 陕西延长石油（集团）有限责任公司 Oil reservoir microbial community competition evolution method
CN117198385B (en) * 2023-11-07 2024-03-19 陕西延长石油（集团）有限责任公司 Oil reservoir microbial community competition evolution method

Similar Documents

Publication Publication Date Title
Wang et al. Inhibitory effect of ethanol, acetic acid, propionic acid and butyric acid on fermentative hydrogen production
Chang et al. Evaluation of different pretreatment methods for preparing hydrogen-producing seed inocula from waste activated sludge
Ke et al. A pilot study on large-scale microbial enhanced oil recovery (MEOR) in Baolige Oilfield
Schink et al. Energetics of syntrophic methane formation and the influence of aggregation
Ke et al. Microbial enhanced oil recovery in Baolige Oilfield using an indigenous facultative anaerobic strain Luteimonas huabeiensis sp. nov
Xia et al. Conversion of petroleum to methane by the indigenous methanogenic consortia for oil recovery in heavy oil reservoir
Zhang et al. Potential of hydrolyzed polyacrylamide biodegradation to final products through regulating its own nitrogen transformation in different dissolved oxygen systems
Bauer et al. Reevaluation of the global warming impacts of algae-derived biofuels to account for possible contributions of nitrous oxide
CN104087534A (en) Activator for activating endogenous microbial oil displacement in polymer-displaced oil reservoir
CN1424484A (en) Method for driving crude oil out by microorgans in crude oil
CN102926728A (en) Indigenous microorganism activation and exogenous microorganism intensified oil production method in offshore oilfield
CN112908425A (en) Microbial oil recovery numerical simulation method based on reaction kinetics model
CN103103130A (en) Production method for lipid through mixed culture of microbes
US20220195466A1 (en) Process and device for the production of methane
Liu et al. Improved formation of biogenic methane by cultivable bacteria in highly volatile bituminous coals
Hellerschmied et al. Hydrogen storage and geo-methanation in a depleted underground hydrocarbon reservoir
Spirov et al. Modelling of microbial enhanced oil recovery application using anaerobic gas-producing bacteria
CN101575634B (en) Method for confirming oxygen consuming content and oxygen consumption rate of oil pool microorganisms
Ashley et al. Deuterium as a quantitative tracer of enhanced microbial methane production
US20160145978A1 (en) Method to enhance microbial gas production from unconventional reservoirs and kerogen deposits
CN108219765A (en) A kind of reservoir endogenous micro-organisms activator and its flooding method based on inorganic salts
Jeon et al. Enrichment of $ CO_2 $-Fixing Bacteria in Cylinder-Type Electrochemical Bioreactor with Built-In Anode Compartment
Liu et al. Effect of simple carbon source on degradation performance of longchain alkane degrading bacteria
Chen et al. Screening of microalgae strains for efficient biotransformation of small molecular organic acids from dark fermentation biohydrogen production wastewater
Zhu et al. Biogenic methane generation from lignite coal at different temperatures

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
WW01 Invention patent application withdrawn after publication
WW01 Invention patent application withdrawn after publication

Application publication date: 20210604