METHOD AND SYSTEM FOR QUANTITATIVELY EVALUATING KEROGEN SWELLING OIL IN SHALE

Abstract

A method and system for quantitatively evaluating kerogen swelling oil in shale is provided. The method includes: establishing different types of kerogen molecular models, and loading each of the kerogen molecular models into a graphene slit-type pore composed of a lamellar structure; performing energy minimization (EM), relaxation and annealing to obtain a kerogen slit-type pore; loading a shale oil molecule into the kerogen slit-type pore to obtain an initial model of swelling and adsorption of shale oil in kerogen; assigning a value to a force field of the shale oil molecule and the kerogen molecule in the swelling and adsorption model to obtain a density of the kerogen and the shale oil; plotting a density curve of the kerogen and the shale oil; calculating kerogen swelling oil mass; determining per-unit kerogen swelling oil mass; and determining the kerogen swelling oil mass in different evolution stages.

Claims

1. A method for quantitatively evaluating kerogen swelling oil in shale, comprising: establishing different types of kerogen molecular models, and loading each of the kerogen molecular models into a graphene slit-type pore composed of a lamellar structure to obtain an initial model; subjecting each initial model to energy minimization (EM) and relaxation to obtain a compacted kerogen aggregate model; simulating an annealing process of the compacted kerogen aggregate model to obtain a kerogen slit-type pore; loading a shale oil molecule into the kerogen slit-type pore to obtain an initial model of swelling and adsorption of shale oil in kerogen; assigning a value to a force field of the shale oil molecule and the kerogen molecule in the initial model of swelling and adsorption of the shale oil in the kerogen to obtain a density of the kerogen and the shale oil; plotting a density curve of the kerogen and the shale oil according to the density of the kerogen and the shale oil; calculating kerogen swelling oil mass according to the density curve of the kerogen and the shale oil; determining per-unit kerogen swelling oil mass according to the kerogen swelling oil mass; and determining the kerogen swelling oil mass in different evolution stages according to the per-unit kerogen swelling oil mass.

2. The method for quantitatively evaluating kerogen swelling oil in shale according to claim 1, wherein the subjecting each initial model to EM and relaxation to obtain a compacted kerogen aggregate model specifically comprises: subjecting the initial model to EM and 200 ps relaxation at 75° C. under 20 MPa by using Gromacs software to obtain a compacted kerogen aggregate model.

3. The method for quantitatively evaluating kerogen swelling oil in shale according to claim 1, wherein the simulating an annealing process of the compacted kerogen aggregate model to obtain a kerogen slit-type pore specifically comprises: subjecting the compacted kerogen aggregate model to 200 ps relaxation for warming; and subjecting the kerogen aggregate model after the relaxation warming to 2 ns simulation, cooling and pressurization by using an isothermal-isobaric ensemble (NPT ensemble) at 800° C. under normal pressure to obtain a kerogen slit-type pore.

4. The method for quantitatively evaluating kerogen swelling oil in shale according to claim 1, wherein the calculating kerogen swelling oil mass according to the density curve of the kerogen and the shale oil specifically comprises: calculating the kerogen swelling oil mass by Q.sub.oil=∫.sub.L.sub.o1.sup.L.sup.o2S.sub.model.Math.ρ.sub.oildL according to the density curve of the kerogen and the shale oil; wherein, Q.sub.oil is the kerogen swelling oil mass; L.sub.o1 is a start position of an intersection between the kerogen density curve and the shale oil density curve; L.sub.o2 is a stop position of the intersection between the kerogen density curve and the shale oil density curve; S.sub.model is a cross-sectional area of the swelling and adsorption model of the shale oil in the kerogen; ρ.sub.oil is the shale oil density curve.

5. The method for quantitatively evaluating kerogen swelling oil in shale according to claim 1, wherein the determining per-unit kerogen swelling oil mass according to the kerogen swelling oil mass specifically comprises: deriving the mass of the kerogen; and determining the per-unit kerogen swelling oil mass according to the kerogen adsorbed oil mass and the kerogen mass.

6. The method for quantitatively evaluating kerogen swelling oil in shale according to claim 1, wherein the determining the kerogen swelling oil mass in different evolution stages according to the per-unit kerogen swelling oil mass specifically comprises: deriving the mass of the kerogen in different evolution stages; calculating the kerogen swelling oil mass in different evolution stages by Q.sub.s=Q.sub.w.Math.m.sub.k.Math.f.sub.s according to the per-unit kerogen swelling oil mass and the kerogen mass in different evolution stages; wherein, Q.sub.w is the kerogen swelling oil mass; m.sub.k is the kerogen mass in different evolution stages; f.sub.s is a swelling ratio reduction coefficient in different evolution stages.

7. A system for quantitatively evaluating kerogen swelling oil in shale, comprising: a first initial model establishment module, for establishing different types of kerogen molecular models, and loading each of the kerogen molecular models into a graphene slit-type pore composed of a lamellar structure to obtain an initial model; a compaction model, for subjecting each initial model to EM and relaxation to obtain a compacted kerogen aggregate model; an annealing model, for simulating an annealing process of the compacted kerogen aggregate model to obtain a kerogen slit-type pore; a second initial model establishment module, for loading a shale oil molecule into the kerogen slit-type pore to obtain an initial model of swelling and adsorption of shale oil in kerogen; an assigning module, for assigning a value to a force field of the shale oil molecule and the kerogen molecule in the initial model of swelling and adsorption of the shale oil in the kerogen to obtain a density of the kerogen and the shale oil; a density curve plotting module, for plotting a density curve of the kerogen and the shale oil according to the density of the kerogen and the shale oil; a kerogen swelling oil mass determination module, for deriving the kerogen swelling oil mass according to the density curve of the kerogen and the shale oil; a per-unit kerogen swelling oil mass determination module, for determining per-unit kerogen swelling oil mass according to the kerogen swelling oil mass; and a different-evolution-stage kerogen swelling oil mass determination module, for determining the kerogen swelling oil mass in different evolution stages according to the per-unit kerogen swelling oil mass.

8. The system for quantitatively evaluating kerogen swelling oil in shale according to claim 7, wherein the compaction module specifically comprises: a compaction unit, for subjecting the initial model to EM and 200 ps relaxation at 75° C. under 20 MPa by using Gromacs software to obtain a compacted kerogen aggregate model.

9. The system for quantitatively evaluating kerogen swelling oil in shale according to claim 7, wherein the annealing module specifically comprises: a warming unit, for subjecting the compacted kerogen aggregate model to 200 ps relaxation for warming; and an annealing unit, for subjecting the kerogen aggregate model after the relaxation warming to 2 ns simulation, cooling and pressurization by using a NPT ensemble at 800° C. under normal pressure to obtain a kerogen slit-type pore.

10. The system for quantitatively evaluating kerogen swelling oil in shale according to claim 7, wherein the kerogen swelling oil mass determination module specifically comprises: a kerogen swelling oil mass determination unit, for deriving the kerogen swelling oil mass by Q.sub.oil=∫.sub.L.sub.o1.sup.L.sup.o2S.sub.model.Math.ρ.sub.oildL according to the density curve of the kerogen and the shale oil; wherein, Q.sub.oil is the kerogen swelling oil mass; L.sub.o1 is a start position of an intersection between the kerogen density curve and the shale oil density curve; L.sub.o2 is a stop position of the intersection between the kerogen density curve and the shale oil density curve; S.sub.model is a cross-sectional area of the swelling and adsorption model of the shale oil in the kerogen; ρ.sub.oil is the shale oil density curve.

Description

BRIEF DESCRIPTION OF DRAWINGS

[0052] In order to describe the technical solutions in the examples of the present invention or in the prior art more clearly, the accompanying drawings required for the examples are briefly described below. Apparently, the accompanying drawings in the following description show merely some examples of the present invention, and a person of ordinary skill in the art may still derive other accompanying drawings from these accompanying drawings without creative efforts.

[0053] FIG. 1 is a flowchart of a method for quantitatively evaluating kerogen swelling oil in shale according to the present invention.

[0054] FIG. 2 shows a process of establishing a kerogen wall model according to the present invention.

[0055] FIG. 3 shows a simulation of an annealing process of a kerogen aggregate model according to the present invention.

[0056] FIG. 4 shows an initial model of swelling and adsorption of shale oil in different types of kerogen according to the present invention.

[0057] FIG. 5 shows a molecular dynamics simulation of swelling and adsorption in type II kerogen according to the present invention.

[0058] FIG. 6 shows a parameter necessary for a quantitative calculation process of kerogen swelling oil and a calculation result according to the present invention.

[0059] FIG. 7 is a structural diagram of a system for quantitatively evaluating kerogen swelling oil in shale according to the present invention.

DETAILED DESCRIPTION

[0060] The technical solutions in the examples of the present invention are clearly and completely described with reference to the accompanying drawings in the examples of the present invention. Apparently, the described examples are merely a part rather than all of the examples of the present invention. All other examples obtained by a person of ordinary skill in the art based on the examples of the present invention without creative efforts should fall within the protection scope of the present invention.

[0061] An objective of the present invention is to provide a method and system for quantitatively evaluating kerogen swelling oil in shale. The present invention improves the accuracy of the quantitative evaluation of the kerogen swelling oil.

[0062] In order to make the above objectives, features, and advantages of the present invention clearer and more comprehensible, the present invention is described in further detail below with reference to the accompanying drawings and specific implementations.

[0063] FIG. 1 is a flowchart of a method for quantitatively evaluating kerogen swelling oil in shale according to the present invention. As shown FIG. 1, a method for quantitatively evaluating kerogen swelling oil in shale includes:

[0064] Step 101: establish different types of kerogen molecular models, and load each of the kerogen molecular models into a graphene slit-type pore composed of a lamellar structure to obtain an initial model.

[0065] Avogadro software was used to establish type I, II and III kerogen molecular models. As a preferred example, only the type II kerogen molecular model is used for illustration. Specifically, Packmol software was used to load 100 type II kerogen molecules into the graphene slit-type pore composed of a lamellar structure (approximately 7.38 nm×7.67 nm×0.85 nm). FIG. 2a shows an initial kerogen wall model.

[0066] Step 102: subject each initial model to energy minimization (EM) and relaxation to obtain a compacted kerogen aggregate model, specifically as follows:

[0067] Subject the initial model to EM and 200 ps relaxation at 75° C. under 20 MPa by using an isothermal-isobaric ensemble (NPT ensemble) of Gromacs software to obtain a compacted kerogen aggregate model.

[0068] FIG. 2 shows a process of establishing a kerogen wall model according to the present invention, where a shows an initial model, b shows a model after 50 ps relaxation, c shows a model after 100 ps relaxation, d shows a model after 150 ps relaxation, and e shows a model after 200 ps relaxation.

[0069] Step 103: simulate an annealing process of the compacted kerogen aggregate model to obtain a kerogen slit-type pore, specifically as follows:

[0070] Subject the compacted kerogen aggregate model to 200 ps relaxation for warming.

[0071] Subject the kerogen aggregate model after the relaxation warming to 2 ns simulation, cooling and pressurization by using the NPT ensemble at 800° C. under normal pressure to obtain a kerogen slit-type pore.

[0072] FIG. 3 shows a simulation of an annealing process of a kerogen aggregate model according to the present invention, where, a shows a result after 2 ns simulation at 800° C. under normal pressure, and b shows a result after 2 ns simulation at 75° C. under 20 MPa. In this step, the compacted kerogen aggregate model was warmed by 200 ps relaxation, and the system was reduced to normal pressure. Then the entire system was simulated by using the NPT ensemble at 800° C. under normal pressure for 2 ns. FIG. 3a shows a final frame of the simulation. Then the system was cooled and pressurized. The entire system was simulated by using the NPT ensemble at 75° C. under 20 MPa for 2 ns. FIG. 3b shows a final frame of the simulation.

[0073] Step 104: load a shale oil molecule into the kerogen slit-type pore to obtain an initial model of swelling and adsorption of shale oil in kerogen.

[0074] Packmol software was used to load the shale oil molecule into the kerogen slit-type pore obtained in step 103. FIG. 4 shows an initial model after loading, where the model is composed of a kerogen wall model (after the annealing process simulated) on double sides and a shale oil model in the middle.

[0075] FIG. 4 shows an initial model of swelling and adsorption of shale oil in different types of kerogen. In the figure, a shows a swelling and adsorption model of shale oil in type I kerogen, b shows a swelling and adsorption model of shale oil in type II kerogen, and c shows a swelling and adsorption model of shale oil in type III kerogen.

[0076] Step 105: assign a value to a force field of the shale oil molecule and the kerogen molecule in the initial model of swelling and adsorption of the shale oil in the kerogen to obtain a density of the kerogen and the shale oil.

[0077] FIG. 5 shows a molecular dynamics simulation of swelling and adsorption in type II kerogen, where a shows a final frame of the simulation, and b shows a density curve of the kerogen and the shale oil. The Chemistry at Harvard Macromolecular Mechanics (CHARMM) 36/CHARMM general (CGenFF) force field was used to assign a value to the force field of the shale oil molecule and the kerogen molecule in the initial model of swelling and adsorption of the shale oil in the kerogen obtained in step 104. The Lorentz-Berthelot (LB) mixing rule was used to calculate an interaction force between the shale oil molecule and the kerogen molecule. The Particle-Mesh-Ewald (PME) model was used as an electrostatic force model. The van der Waals radius was 1.4 nm. The model after the force field assignment was simulated by using the NPT ensemble of Gromacs software at 75° C. under 20 MPa for 200 ns. Taking type II kerogen as an example, the simulation result is shown in FIG. 5a. The density curves of the kerogen and the shale oil were plotted (step 106), as shown in FIG. 5b.

[0078] Step 106: plot a density curve of the kerogen and the shale oil according to the density of the kerogen and the shale oil.

[0079] Step 107: derive the kerogen swelling oil mass according to the density curve of the kerogen and the shale oil, specifically as follows:

[0080] calculate the kerogen swelling oil mass by Q.sub.oil=∫.sub.L.sub.o1.sup.L.sup.o2S.sub.model.Math.ρ.sub.oildL according to the density curve of the kerogen and the shale oil.

[0081] In the formula, Q.sub.oil is the kerogen swelling oil mass; L.sub.o1 is a start position of an intersection between the kerogen density curve and the shale oil density curve; L.sub.o2 is a stop position of the intersection between the kerogen density curve and the shale oil density curve; S.sub.model is a cross-sectional area of the swelling and adsorption model of the shale oil in the kerogen; ρ.sub.oil is the shale oil density curve.

[0082] Step 108: determine per-unit kerogen swelling oil mass according to the kerogen swelling oil mass, specifically as follows:

[0083] derive the mass of the kerogen; and

[0084] determine the per-unit kerogen swelling oil mass according to the kerogen adsorbed oil mass and the kerogen mass.

[0085] The swelling oil mass in the type I kerogen, the swelling oil mass in the type II kerogen and the swelling oil mass in the type III kerogen were 161.04 mg/g TOC, 104.96 mg/g TOC and 70.29 mg/g TOC, respectively.

[0086] Step 109: determine the kerogen swelling oil mass in different evolution stages according to the per-unit kerogen swelling oil mass, specifically as follows:

[0087] derive the mass of the kerogen in different evolution stages; and

[0088] calculate the kerogen swelling oil mass in different evolution stages by Q.sub.s=Q.sub.w.Math.m.sub.k.Math.f.sub.s according to the per-unit kerogen swelling oil mass and the kerogen mass in different evolution stages.

[0089] In the formula, Q.sub.w is the kerogen swelling oil mass; m.sub.k is the kerogen mass in different evolution stages; f.sub.s is a swelling ratio reduction coefficient in different evolution stages.

[0090] FIG. 6 shows a parameter necessary for a quantitative calculation process of kerogen swelling oil and a calculation result according to the present invention. Taking the kerogen swelling oil mass of the Qing-1 Member in the Songliao Basin as an example, as the shale in the Qing-1 Member is dominated by the type I kerogen, a swelling oil mass in the type I kerogen derived by a molecular dynamics simulation was used as an initial swelling mass. The ability of the oil to swell the kerogen decreases continuously with the evolution stage. 1 g of original organic carbon was used to normalize the kerogen swelling oil mass. The initial swelling mass multiplied by the mass of the type I kerogen in different evolution stages and by the swelling ratio reduction coefficient (FIG. 6a) to obtain the swelling oil mass in the type I kerogen in different evolution stages.

[0091] Swelling mass in the type I kerogen in different evolution stages:

Q.sub.s=Q.sub.w.Math.m.sub.k.Math.f.sub.s

[0092] In the formula, Q.sub.w is the swelling oil mass in the type I kerogen, 161.04 mg/g; m.sub.k is the kerogen mass in different evolution stages; f.sub.s is the swelling ratio reduction coefficient in different evolution stages, dimensionless.

[0093] Kerogen mass correspondingto 1 g of original organic carbon:

m.sub.k=m.sub.f.Math.F.sub.t+m.sub.s=(HI.sup.0/1000).Math.F.sub.t+(1−HI.sup.0.Math.0.083/100)  (3)

[0094] In the formula, m.sub.per-unit is the mass of a convertible part of the kerogen, g; m.sub.s is the mass of a non-convertible part of the kerogen, g; HI.sup.0 is an original hydrogen index, mg/g TOC; 0.083 is a carbon conversion coefficient of the hydrogen index, dimensionless; F.sub.t is a conversion rate, dimensionless. The HI.sup.0 of the type I kerogen in the shale of the Qing-1 Member in the north of the Songliao Basin was determined as 750 mg/g TOC according to the geochemical data thereof.

[0095] Taking the shale of the Qing-1 Member in the north of the Songliao Basin as an example, a chemical kinetics parameter of the primary cracking of the kerogen was calibrated based on a pyrolysis gas chromatography (PY-GC) test result of an immature shale sample of the Qing-1 Member in the Du 402 well in the Taikang uplift in the north of the Songliao Basin, as well as the principles of chemical kinetics. Based on the burial history and thermal history of the Songliao Basin, the conversion rates corresponding to different R.sub.o were calculated (FIG. 6b), and the kerogen mass and kerogen swelling oil mass in different evolution stages were calculated according to Formulas 2 and 3 (FIGS. 6c and 6d).

[0096] Steps 101-106 of the present invention used a real kerogen model, which overcame the problems caused by the conventional method of using simple graphene instead of kerogen to study the kerogen-shale oil interaction. As graphene is a two-dimensional carbon material with a smooth surface and a surface structure different from kerogen molecules, the shale oil molecules cannot pass through the graphene and enter its lamellar structure to make the graphene swell like the real kerogen structure. The kerogen aggregate model was preprocessed to prevent incomplete compaction of a kerogen aggregate. If the kerogen aggregate was incompletely compacted, it would result in the existence of “large pores” inside, making the density of the kerogen aggregate model lower than that of the kerogen sample.

[0097] In steps 107-109 of the present invention, the processing of the per-unit swelling oil mass overcame the problem that a molecular dynamics simulation system is too small to be applied to the swelling shale oil. A shale oil-kerogen system used for molecular dynamics simulation is usually less than 20 nm, while the pore size of most shale reservoirs is greater than 20 nm. The present invention calculated the per-unit kerogen swelling oil mass, and calculated the actual kerogen swelling oil mass based on actual geological parameters such as swelling coefficient reduction and kerogen mass change, which greatly improved the accuracy of the results.

[0098] The present invention further provides a system for quantitatively evaluating kerogen swelling oil in shale. FIG. 7 is a structural diagram of a system for quantitatively evaluating kerogen swelling oil in shale according to the present invention. As shown FIG. 7, a system for quantitatively evaluating kerogen swelling oil in shale includes:

[0099] a first initial model establishment module 201, for establishing different types of kerogen molecular models, and loading each of the kerogen molecular models into a graphene slit-type pore composed of a lamellar structure to obtain an initial model;

[0100] a compaction model 202, for subjecting each initial model to EM and relaxation to obtain a compacted kerogen aggregate model;

[0101] an annealing model 203, for simulating an annealing process of the compacted kerogen aggregate model to obtain a kerogen slit-type pore;

[0102] a second initial model establishment module 204, for loading a shale oil molecule into the kerogen slit-type pore to obtain an initial model of swelling and adsorption of shale oil in kerogen;

[0103] an assigning module 205, for assigning a value to a force field of the shale oil molecule and the kerogen molecule in the initial model of swelling and adsorption of the shale oil in the kerogen to obtain a density of the kerogen and the shale oil;

[0104] a density curve plotting module 206, for plotting a density curve of the kerogen and the shale oil according to the density of the kerogen and the shale oil;

[0105] a kerogen swelling oil mass determination module 207, for deriving the kerogen swelling oil mass according to the density curve of the kerogen and the shale oil;

[0106] a per-unit kerogen swelling oil mass determination module 208, for determining per-unit kerogen swelling oil mass according to the kerogen swelling oil mass; and

[0107] a different-evolution-stage kerogen swelling oil mass determination module 209, for determining the kerogen swelling oil mass in different evolution stages according to the per-unit kerogen swelling oil mass.

[0108] The compaction module 202 specifically includes:

[0109] a compaction unit, for subjecting the initial model to EM and 200 ps relaxation at 75° C. under 20 MPa by using Gromacs software to obtain a compacted kerogen aggregate model.

[0110] The annealing module 203 specifically includes:

[0111] a warming unit, for subjecting the compacted kerogen aggregate model to 200 ps relaxation for warming; and

[0112] an annealing unit, for subjecting the kerogen aggregate model after the relaxation warming to 2 ns simulation, cooling and pressurization by using a NPT ensemble at 800° C. under normal pressure to obtain a kerogen slit-type pore.

[0113] The kerogen swelling oil mass determination module 207 specifically includes:

[0114] a kerogen swelling oil mass determination unit, for deriving the kerogen swelling oil mass by Q.sub.oil=∫.sub.L.sub.o1.sup.L.sup.o2S.sub.model.Math.ρ.sub.oildL according to the density curve of the kerogen and the shale oil.

[0115] In the formula, Q.sub.oil is the kerogen swelling oil mass; L.sub.o1 is a start position of an intersection between the kerogen density curve and the shale oil density curve; L.sub.o2 is a stop position of the intersection between the kerogen density curve and the shale oil density curve; S.sub.model is a cross-sectional area of the swelling and adsorption model of the shale oil in the kerogen; ρ.sub.oil is the shale oil density curve.

[0116] Each example of the present specification is described in a progressive manner, each example focuses on the difference from other examples, and the same and similar parts between the examples may refer to each other. For a system disclosed in the examples, since it corresponds to the method disclosed in the examples, the description is relatively simple, and reference can be made to the method description.

[0117] In this paper, several examples are used for illustration of the principles and implementations of the present invention. The description of the foregoing examples is used to help illustrate the method of the present invention and the core principles thereof. In addition, those of ordinary skill in the art can make various modifications in terms of specific implementations and scope of application in accordance with the teachings of the present invention. In conclusion, the content of this specification should not be construed as a limitation to the present invention.