CN109854236A - A kind of method for numerical simulation and device for shale oil flowing - Google Patents
A kind of method for numerical simulation and device for shale oil flowing Download PDFInfo
- Publication number
- CN109854236A CN109854236A CN201910137987.4A CN201910137987A CN109854236A CN 109854236 A CN109854236 A CN 109854236A CN 201910137987 A CN201910137987 A CN 201910137987A CN 109854236 A CN109854236 A CN 109854236A
- Authority
- CN
- China
- Prior art keywords
- shale
- amount
- oil
- unit volume
- reaction
- 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
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH DRILLING; MINING
- E21B—EARTH DRILLING, e.g. DEEP 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
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH DRILLING; MINING
- E21B—EARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH DRILLING; MINING
- E21B—EARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
- E21B49/08—Obtaining fluid samples or testing fluids, in boreholes or wells
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
Abstract
The invention discloses a kind of method for numerical simulation and device for shale oil flowing, belong to unconventional oil and gas field Development Engineering technical field.Method for numerical simulation and device provided by the invention describe flow characteristics in different existence state shale oil production process by establishing reaction model, including free state crude oil seepage flow, adsorb absorption of the mutual solvent crude oil in organic matter, desorption behavior, influence etc. of the capillary force to solution gas occurrence status, in conjunction with laboratory physical simulation experiment, determine chemical reaction equation, the parameters such as reaction rate and the order of reaction (parameters are demarcated), finally by numerical reservoir simulation method, establish the numerical simulator for considering shale oil Complex Flows mechanism, to realize the simulation to unconventional shale oil seepage flow characteristics.
Description
Technical field
The present invention relates to unconventional oil and gas field Development Engineering technical field, in particular to a kind of number for shale oil flowing
It is worth analogy method and device.
Background technique
Shale oil refers to the petroleum being enriched in rich organic matter black shale stratum, with sides such as free, absorption and solubilised states
Formula exists, and general oil is mainly preserved in nanoscale pore throat and Fracture System, compared with light, viscosity is lower mostly along sheet bed plane
Or microcrack distribution in parallel.The main feature of shale oil shows themselves in that rich in organic matter there is absorption solution in recovery process
Suction phenomenon；Porosity type complexity (including organic hole, inorganic hole, microcrack)；Fluid occurrence form multiplicity (including it is free
State, ADSORPTION STATE, solubilised state) etc..
The oil reservoir routinely recognized can usually be stored up since fluid flowing meets Darcy's law by logging method
Layer physical property distribution situation, and then can fluid flow progress numerical simulation.
Inventor discovery in the prior art the prior art has at least the following problems:
Due to being rich in organic matter in shale reservoir and there are a large amount of micro-nano holes, there are hairs during oil-gas two-phase flow is dynamic
The effect of pipe power, so that there is the multiple fluids migration mechanism such as desorption, diffusion, Darcy Flow and Non-Darcy's flow in shale oil, and can not
The flow mechanism of shale oil is described using existing method for numerical simulation.
Summary of the invention
The embodiment of the invention provides a kind of method for numerical simulation and device for shale oil flowing, this method and device
The Adsorption and desorption process and solution gas of organic matter and oil in shale are described in micro-nano hole by chemical reaction equation
Conversion process is then based on reservoir numerical simulation software, flows to shale oil and carries out numerical simulation.
Specifically, including technical solution below:
On the one hand, a kind of method for numerical simulation for shale oil flowing is provided, comprising the following steps:
The pressure when chemical reaction equation and the chemical reaction equation for determining description shale oil seepage flow characteristics react
Power condition；
Determine the amount of maximum adsorption crude petroleum materials in unit volume shale；
According to the amount, dissolved gas oil ratio and strata pressure item of maximum adsorption crude petroleum materials in the unit volume shale
Part, determine can be desorbed in unit volume shale oil substances amount, can stripping gas substance amount and the amounts of crude petroleum materials can not be desorbed；
According to the amounts of oil substances can be desorbed in the unit volume shale, can stripping gas substance in the unit volume shale
Amount and the unit volume shale in can not desorb the amounts of crude petroleum materials, determine amendment porosity；
Determine the corresponding order of reaction of the chemical reaction equation and reaction rate；
According to maximum adsorption crude petroleum materials in the chemical reaction equation, the pressure condition, the unit volume shale
Amount, the amount, amount for the substance that can vent one's spleen in the unit volume shale, described that oil substances can be desorbed in the unit volume shale
Amount, the amendment porosity, the order of reaction and the reaction rate of crude petroleum materials can not be desorbed in unit volume shale,
Based on reservoir numerical simulation software, numerical simulator is established.
In a kind of possible design, the chemical reaction equation includes:
Reaction 1: kerogen oil → kerogen+oil；
Reaction 2: kerogen+oil → kerogen oil；
Reaction 3: kerogen solution gas → kerogen+solution gas；
Reaction 4: kerogen+solution gas → kerogen solution gas；
Reaction 5: solution gas → dispersion gas；
Reaction 6: dispersion gas → continuous gas.
In a kind of possible design, the amount of maximum adsorption crude petroleum materials passes through following formula in the unit volume shale
It obtains:
Wherein, NoFor the amount of maximum adsorption crude petroleum materials in unit volume shale；TOC is shale organic carbon content；K1It is normal
Number；MoFor the molecular weight of crude oil；ρrFor shale density.
In a kind of possible design, the amounts of oil substances, the unit volume page can be desorbed in the unit volume shale
In rock can stripping gas substance amount and the unit volume shale in can not desorb the amounts of crude petroleum materials and obtained by following formula:
Nk·o=TOCK2(Pi-P)Co；
Nk·g=TOCK2(Pi-P)Cg；
Nk1=No-Nk·o-Nk·g；
Wherein, Nk·oFor the amount that can desorb oil substances in unit volume shale；Nk·gFor can stripping gas in unit volume shale
The amount of substance；Nk1For the amount that can not desorb crude petroleum materials in unit volume shale；K2For constant；PiFor oil reservoir initial pressure, P is
Produce the reservoir pressure in latter stage, CoFor the amount score of substance oily under original state, CgFor the amount point of the substance of gas under original state
Number.
In a kind of possible design, the amendment porosity is obtained by following formula:
Wherein, φ1To correct porosity；φ0For primary porosity；ρk.oFor rubbing for the kerogen that can desorb and oily conjugate
That density；ρk.gFor the molar density of the kerogen and gas conjugate that can desorb；ρk1It is the kerogen that can not desorb in conjunction with crude oil
The molar density of object.
In a kind of possible design, the reaction rate is obtained by following formula:
R=rk·(Ck.o+Ck.g)ek；
Wherein, r is reaction rate；rkFor reaction rate constant；Ek is the order of reaction；Ck.oFor the kerogen that can desorb with
The molar concentration of oily conjugate；Ck.gFor the molar concentration of the kerogen and gas conjugate that can desorb.
On the other hand, a kind of numerical simulation device for shale oil flowing is provided, described device includes:
First determining module, for determining the chemical reaction equation for describing shale oil seepage flow characteristics and the chemical reaction side
Pressure condition when journey reacts；
Second determining module, for determining the amount of maximum adsorption crude petroleum materials in unit volume shale；
Third determining module, for according to the amount of maximum adsorption crude petroleum materials, dissolution gas and oil in the unit volume shale
Than and condition of formation pressure, determine can be desorbed in unit volume shale oil substances amount, can stripping gas substance amount and can not
Desorb the amount of crude petroleum materials；
4th determining module, for according to amount, the unit volume that can desorb oil substances in the unit volume shale
In shale can stripping gas substance amount and the unit volume shale in can not desorb the amounts of crude petroleum materials, determine amendment hole
Degree；
5th determining module, for determining the corresponding order of reaction of the chemical reaction equation and reaction rate；
Module is established, for according to maximum in the chemical reaction equation, the pressure condition, the unit volume shale
It adsorbs the amounts of crude petroleum materials, the amount that oil substances can be desorbed in the unit volume shale, can vent one's spleen in the unit volume shale
Can not be desorbed in the amount of substance, the unit volume shale amounts of crude petroleum materials, the amendment porosity, the order of reaction and
The reaction rate is based on reservoir numerical simulation software, establishes numerical simulator.
Technical solution bring beneficial effect provided in an embodiment of the present invention includes at least:
Method for numerical simulation and device provided in an embodiment of the present invention describe different existence state by establishing reaction model
Flow characteristics in shale oil production process, including free state crude oil seepage flow, adsorb absorption of the mutual solvent crude oil in organic matter,
Desorption behavior, influence etc. of the capillary force to solution gas occurrence status determine chemical reaction side in conjunction with laboratory physical simulation experiment
The parameters such as journey, reaction rate and the order of reaction (demarcate parameters), finally by numerical reservoir simulation method, build
The vertical numerical simulator for considering shale oil Complex Flows mechanism, to realize the simulation to unconventional shale oil seepage flow characteristics.
Detailed description of the invention
To describe the technical solutions in the embodiments of the present invention more clearly, make required in being described below to embodiment
Attached drawing is briefly described, it should be apparent that, drawings in the following description are only some embodiments of the invention, for
For those of ordinary skill in the art, without creative efforts, it can also be obtained according to these attached drawings other
Attached drawing.
Fig. 1 is a kind of schematic diagram of method for numerical simulation for shale oil flowing provided in an embodiment of the present invention；
Fig. 2 is a kind of schematic diagram of numerical simulation device for shale oil flowing provided in an embodiment of the present invention.
Specific embodiment
To keep technical solution of the present invention and advantage clearer, below in conjunction with attached drawing to embodiment of the present invention make into
One step it is described in detail.Unless otherwise defined, all technical terms used in the embodiment of the present invention all have and art technology
The normally understood identical meaning of personnel.
On the one hand, the embodiment of the invention provides a kind of method for numerical simulation for shale oil flowing, as shown in Figure 1,
The method for numerical simulation the following steps are included:
Step 101: it is anti-to determine that the chemical reaction equation of description shale oil seepage flow characteristics and the chemical reaction equation occur
Seasonable pressure condition；
Step 102: determining the amount of maximum adsorption crude petroleum materials in unit volume shale；
Step 103: according to amount, dissolved gas oil ratio and the stratum of maximum adsorption crude petroleum materials in the unit volume shale
Pressure condition, determine can be desorbed in unit volume shale oil substances amount, can stripping gas substance amount and crude oil object can not be desorbed
The amount of matter；
Step 104: according to can desorb the amounts of oil substances in the unit volume shale, can solve in the unit volume shale
The amount of crude petroleum materials can not be desorbed in the amount of getter material and the unit volume shale, determine amendment porosity；
Step 105: determining the corresponding order of reaction of the chemical reaction equation and reaction rate；
Step 106: according to maximum adsorption in the chemical reaction equation, the pressure condition, the unit volume shale
The amounts of oil substances can be desorbed in the amount of crude petroleum materials, the unit volume shale, can vent one's spleen in the unit volume shale substance
Amount, amount, the amendment porosity, the order of reaction and described that crude petroleum materials can not be desorbed in the unit volume shale
Reaction rate is based on reservoir numerical simulation software, establishes numerical simulator.
Method for numerical simulation provided in an embodiment of the present invention describes different existence state shale oil by establishing reaction model
Flow characteristics in production process, including free state crude oil seepage flow adsorb Adsorption and desorption row of the mutual solvent crude oil in organic matter
For influence etc. of the capillary force to solution gas occurrence status determines chemical reaction equation, reaction in conjunction with laboratory physical simulation experiment
The parameters such as rate and the order of reaction (demarcating to parameters) are established and are considered finally by numerical reservoir simulation method
The numerical simulator of shale oil Complex Flows mechanism, to realize the simulation to unconventional shale oil seepage flow characteristics.
For step 101, the chemical reaction equation and chemical reaction equation hair of description shale oil seepage flow characteristics are determined
Pressure condition when raw reaction；
It should be noted that shale oil itself does not participate in chemical reaction, but shale oil has complexity in flow process
Occurrence status, " occurrence status " have a major impact the flow characteristics in shale oil production process, therefore " chemical reaction herein
Equation " is used to describe the different existence state of shale oil, that is to say seepage flow characteristics.
In view of kerogen to the suction-operated of crude oil be shale from tight sand the flow mechanism in recovery process is different
The main reason for, therefore define following chemical reaction:
Reaction 1: kerogen oil → kerogen+oil；
Reaction 2: kerogen+oil → kerogen oil；
Reaction 3: kerogen solution gas → kerogen+solution gas；
Reaction 4: kerogen+solution gas → kerogen solution gas.
When only oil in reservoir, with pressure reduction, reaction 1 starts to carry out, on the contrary, with pressure rise, reaction 2 start into
Row；When in reservoir containing oily and solution gas, with pressure reduction, reaction 1, reaction 3 start to carry out, on the contrary, with pressure rise, instead
2, reaction 4 is answered to start to carry out.
Wherein, the mass ratio of the material value of kerogen absorption crude oil does not influence total adsorbance and reaction rate, can simplify setting
For 1:1.Pressure condition when for specifically reacting can be determined by institute's actual conditions to be simulated.
In addition to this, preservation shape of the capillary force effect of the micro-nano hole of shale to solution gas during crude oil two-phase flow
State also has larger impact.When pressure is down to saturation pressure, since capillary force acts on, gas cannot form continuous phase immediately, and
It is to be present in micro-nano hole with dispersed, forms the oil displacement process of similar foam oil, when pressure continues that view is reduced to be saturated
After pressure, gas forms continuous phase, and gas-oil ratio increases, and yield reduces.Therefore, and below definition chemically react:
Reaction 5: solution gas → dispersion gas；
Reaction 6: dispersion gas → continuous gas.
After reservoir pressure is down to saturation pressure, reaction 5 starts to carry out, after reservoir pressure is down to view saturation pressure, instead
6 are answered to start to carry out.Wherein, saturation pressure can be by flashing equation calculation, can also be by determination of experimental method, can also be direct
It tables look-up acquisition；It can be calculated by experimental method correction formula depending on saturation pressure, determination of experimental method can also be passed through.
For step 102, the amount of maximum adsorption crude petroleum materials in unit volume shale is determined；
In view of TOC value (organic carbon content) is the important parameter of cheese radical content in determining shale, and determine shale
Therefore the important evidence of middle absorption solubilised state crude oil amount can determine maximum adsorption crude oil in unit volume shale by TOC value
The amount of substance.Wherein, in the higher shale reservoir of TOC value, the crude oil amount of ADSORPTION STATE is higher, and (this part crude oil is dropped with reservoir pressure
It is low and desorb into hole, to slow down the reduction of reservoir pressure, increase reservoir yield, be conducive to the holding of yield).And it is existing
Result of study is shown, under isothermy, is increased with pressure, and the ability of shale absorption crude oil is gradually increased, and is adsorbed crude oil and contained
Amount has preferable linear positive correlation with organic carbon content.
Further, in DIFFERENT DEPOSITIONAL ENVIRONMENTS, the kerogen formed by separate sources organic matter, nature difference compared with
Greatly, can mainly three types be divided into:
I type kerogen (sapropel type), with containing based on lipids, linear paraffin is more, polycyclic aromatic hydrocarbon and oxygen-containing function
Group is less, has high hydrogen low oxygen content, and oil generation potential is big；II type kerogen, hydrogen content is higher, but slightly lower compared with I type kerogen, is
HI SA highly saturated polycyclic carbon skeleton, linear paraffin containing moderate-length and cycloalkane are more, also contain polycyclic aromatic hydrocarbon and heteroatom functional
Group, derives from marine facies planktonic organism and microorganism, and oil generation potential is medium；There is type III kerogen (humus type) low hydrogen hyperoxia to contain
Amount, it is unfavorable to oil generation with containing based on polycyclic aromatic hydrocarbon and oxygen-containing functional group, saturated hydrocarbons is seldom, derives from land higher plant, but bury
When hiding enough depth, advantageous Oil-gas source can be become.
For different types of kerogen, the molecular weight M of kerogen component can be calculated by the following formulak(Basic Physical Properties
Parameter) and unit volume shale in maximum adsorption crude petroleum materials amount No:
In formula, MoFor the molecular weight of crude oil, KcFor carbon mass fraction, ρ in kerogenrFor shale density, TOC is page
Rock organic carbon content, K1For constant.
MoIt illustratively can be by assuming that kerogen and oil or gas be adsorbed by the amount of 1:1 substance for the molecular weight of crude oil
In the case of, conversion obtains.
TOC value can directly be measured by testing, and can also be obtained by relevant historical data；Its range is usually 0~10%
Between.
K1For constant, K1The measurement of the methods of depressurization desorption, the experiment of vacuum self-priming can be used, concretely unit mass shale
The ratio of maximum adsorption crude quality and TOC；Illustratively, shale TOC value is 2%=0.02, the absorption of every gram of shale most multipotency
0.01g oil, then K1=0.01g/0.02g=0.5.
KcThe methods of combustion experiment, chromatography can be used, can also be obtained by searching for data, it is specific real in combination with research area
Border situation determines.KcIt is distributed in 0.1-0.9.
For step 103: according to the amount, dissolved gas oil ratio of maximum adsorption crude petroleum materials in the unit volume shale and
Condition of formation pressure, determine can be desorbed in unit volume shale oil substances amount, can stripping gas substance amount and original can not be desorbed
The amount of oil substances；
For shale, with pressure reduction, shale adsorbance is gradually decreased, and pressure and adsorbance are in preferable linear relationship.
Relationship between adsorbance and pressure can be described with following formula:
Nk=K2p+C；
In formula, NkFor the amount for actually adsorbing crude petroleum materials in unit volume shale, p is reservoir pressure, K2It is constant with C, it can
It is determined by experiment or is directly determined according to Kerogen type according to the actual situation.
Further, according to NkIt can in the case where known oil reservoir initial pressure and production terminal pressure with TOC content
Calculate can be desorbed from kerogen oil substances amount and can stripping gas substance amount, that is to say decompression recovery process in, unit bodies
The kerogen of reaction and the amount of oily conjugate, kerogen and gas conjugate substance are participated in product shale.Formula is as follows:
Nk·o=TOCK2(Pi-P)Co；
Nk·g=TOCK2(Pi-P)Cg；
Nk1=No-Nk·o-Nk·g；
Wherein, Nk·oFor the amount that can desorb oil substances in unit volume shale；Nk·gFor can stripping gas in unit volume shale
The amount of substance；Nk1For the amount that can not desorb crude petroleum materials in unit volume shale；K2For constant；PiFor oil reservoir initial pressure, P is
Produce the reservoir pressure in latter stage, CoFor the amount score of substance oily under original state, CgFor the amount point of the substance of gas under original state
Number.
PiIt can be measured according to actual production conditions with P, CoAnd CgIt can be directly calculated according to known dissolved gas oil ratio；
K2It can be obtained by the relation curve between above-mentioned adsorbance and pressure.
For step 104: according to the amounts of oil substances can be desorbed in the unit volume shale, in the unit volume shale
Can stripping gas substance amount and the unit volume shale in can not desorb the amounts of crude petroleum materials, determine amendment porosity；
Since the conjugate of the conjugate of kerogen and kerogen and oil, kerogen and gas is present in shale in solid form
In hole, and with pressure reduction, oil is gentle to be desorbed from kerogen, solid-state kerogen and fluid oil, solution gas is generated, in order to more
The flowing of true simulation shale oil, needs to be modified rammell primary porosity.
Specifically, it can be calculated by the following formula to obtain:
Wherein, φ1To correct porosity；φ0For primary porosity；ρk.oFor rubbing for the kerogen that can desorb and oily conjugate
That density；ρk.gFor the molar density of the kerogen and gas conjugate that can desorb；ρk1It is the kerogen that can not desorb in conjunction with crude oil
The molar density of object.
φ0ρ can be obtained by well-log informationk.o、ρk.g、ρk1It can be calculated by " molar density=density/molecular weight "
It arrives.Herein, density is it is known that kerogenic molecular weight MkIt has been computed and has provided above, oil gentle molecular weight is it is known that can desorb
Molecular weight=kerogenic molecular weight+oil molecular weight of kerogen and oily conjugate, the kerogen that can be desorbed and gas conjugate
Molecular weight=kerogenic molecular weight+gas molecular weight, the molecular weight of the kerogen that can not be desorbed and crude oil conjugate=dry
The molecular weight of junket root+crude oil average molecular weight.
For step 105, the corresponding order of reaction of chemical reaction equation and reaction rate are determined；
Reaction rate can be obtained by following formula:
R=rk·(Ck.o+Ck.g)ek；
Wherein, r is reaction rate；rkFor reaction rate constant；Ek is the order of reaction；Ck.oFor the kerogen that can desorb with
The molar concentration of oily conjugate；Ck.gFor the molar concentration of the kerogen and gas conjugate that can desorb.
rk, ek can pass through depressurization desorption measuring, Ck.o、Ck.gN can be passed throughk·o、Nk·gWith primary porosity φ0It calculates
It arrives.
For step 106, according to maximum in the chemical reaction equation, the pressure condition, the unit volume shale
It adsorbs the amounts of crude petroleum materials, the amount that oil substances can be desorbed in the unit volume shale, can vent one's spleen in the unit volume shale
Can not be desorbed in the amount of substance, the unit volume shale amounts of crude petroleum materials, the amendment porosity, the order of reaction and
The reaction rate is based on reservoir numerical simulation software, establishes numerical simulator.
Reservoir numerical simulation software is a kind of for simulating the software of Reservoir behavior in reservoir, and the embodiment of the present invention is used
Software can be CMG reservoir numerical simulation software.
In application, by determining reactional equation, reaction condition, reaction rate, the order of reaction, the amount of reactant species, repairing
Parameters for numerical simulation and the Basic Physical Properties parameters such as positive porosity, are updated to reservoir numerical simulation software, can establish can describe
The numerical simulator of shale oil feature, and then carry out shale oil exploitation numerical simulation work.
On the other hand, the embodiment of the invention also provides a kind of numerical simulation devices for shale oil flowing, such as Fig. 2 institute
Show, which includes:
First determining module 201, the chemical reaction equation and the chemistry for determining description shale oil seepage flow characteristics are anti-
Pressure condition when equation being answered to react；
Second determining module 202, for determining the amount of maximum adsorption crude petroleum materials in unit volume shale；
Third determining module 203, for according to the amount of maximum adsorption crude petroleum materials, solution gas in the unit volume shale
Oil ratio and condition of formation pressure, determine can be desorbed in unit volume shale oil substances amount, can stripping gas substance amount and not
The amount of crude petroleum materials can be desorbed；
4th determining module 204, for according to amount, the unit bodies that can desorb oil substances in the unit volume shale
In product shale can stripping gas substance amount and the unit volume shale in can not desorb the amounts of crude petroleum materials, determine amendment hole
Degree；
5th determining module 205, for determining the corresponding order of reaction of the chemical reaction equation and reaction rate；
Module 206 is established, for according in the chemical reaction equation, the pressure condition, the unit volume shale
The amounts of oil substances can be desorbed in the amount of maximum adsorption crude petroleum materials, the unit volume shale, can in the unit volume shale
Vent one's spleen the amount of substance, amount, the amendment porosity, the order of reaction of crude petroleum materials can not be desorbed in the unit volume shale
The several and reaction rate is based on reservoir numerical simulation software, establishes numerical simulator.
The numerical simulation device of shale oil flowing provided in an embodiment of the present invention, describes different taxes by establishing reaction model
The flow characteristics in state shale oil production process, including free state crude oil seepage flow are deposited, adsorbs mutual solvent crude oil in organic matter
Adsorption and desorption behavior, influence etc. of the capillary force to solution gas occurrence status determine chemistry in conjunction with laboratory physical simulation experiment
The parameters such as reactional equation, reaction rate and the order of reaction (demarcate parameters), finally by reservoir numerical simulation
Method establishes the numerical simulator for considering shale oil Complex Flows mechanism, to realize to unconventional shale oil seepage flow characteristics
Simulation.
It should be noted that the numerical simulation device of shale oil flowing provided by the above embodiment is in simulation shale oil flowing
When, only the example of the division of the above functional modules, in practical application, it can according to need and divide above-mentioned function
With being completed by different functional modules, i.e., the internal structure of device is divided into different functional modules, to complete above description
All or part of function.In addition, the numerical simulation device and shale oil of shale oil flowing provided by the above embodiment flow
Method for numerical simulation embodiment belong to same design, specific implementation process is detailed in embodiment of the method, and which is not described herein again.
Those of ordinary skill in the art will appreciate that realizing that all or part of the steps of above-described embodiment can pass through hardware
It completes, relevant hardware can also be instructed to complete by program, the program can store in a kind of computer-readable
In storage medium, storage medium mentioned above can be read-only memory, disk or CD etc..
The above is merely for convenience of it will be understood by those skilled in the art that technical solution of the present invention, not to limit
The present invention.All within the spirits and principles of the present invention, any modification, equivalent replacement, improvement and so on should be included in this
Within the protection scope of invention.
Claims (7)
1. a kind of method for numerical simulation for shale oil flowing, which comprises the following steps:
The pressure strip when chemical reaction equation and the chemical reaction equation for determining description shale oil seepage flow characteristics react
Part；
Determine the amount of maximum adsorption crude petroleum materials in unit volume shale；
According to amount, dissolved gas oil ratio and the condition of formation pressure of maximum adsorption crude petroleum materials in the unit volume shale, really
Determine can to desorb in unit volume shale oil substances amount, can stripping gas substance amount and the amounts of crude petroleum materials can not be desorbed；
According to the amounts of oil substances can be desorbed in the unit volume shale, in the unit volume shale can stripping gas substance amount
With the amount that can not desorb crude petroleum materials in the unit volume shale, amendment porosity is determined；
Determine the corresponding order of reaction of the chemical reaction equation and reaction rate；
According to the amount of maximum adsorption crude petroleum materials in the chemical reaction equation, the pressure condition, the unit volume shale,
Amount, the unit of the amounts of oil substances, the substance that can vent one's spleen in the unit volume shale can be desorbed in the unit volume shale
Amount, the amendment porosity, the order of reaction and the reaction rate that crude petroleum materials can not be desorbed in volume shale, are based on
Reservoir numerical simulation software, establishes numerical simulator.
2. method for numerical simulation according to claim 1, which is characterized in that the chemical reaction equation includes:
Reaction 1: kerogen oil → kerogen+oil；
Reaction 2: kerogen+oil → kerogen oil；
Reaction 3: kerogen solution gas → kerogen+solution gas；
Reaction 4: kerogen+solution gas → kerogen solution gas；
Reaction 5: solution gas → dispersion gas；
Reaction 6: dispersion gas → continuous gas.
3. method for numerical simulation according to claim 1, which is characterized in that maximum adsorption is former in the unit volume shale
The amount of oil substances is obtained by following formula:
Wherein, NoFor the amount of maximum adsorption crude petroleum materials in unit volume shale；TOC is shale organic carbon content；K1For constant；
MoFor the molecular weight of crude oil；ρrFor shale density.
4. method for numerical simulation according to claim 3, which is characterized in that can denuded oil object in the unit volume shale
In the amount of matter, the unit volume shale can stripping gas substance amount and the unit volume shale in can not desorb crude petroleum materials
Amount obtained by following formula:
Nk.o=TOCK2(Pi-P)Co；
Nk.g=TOCK2(Pi-P)Cg；
Nk1=N0-Nk.o-Nk.g；
Wherein, Nk.oFor the amount that can desorb oil substances in unit volume shale；Nk.gFor can stripping gas substance in unit volume shale
Amount；Nk1For the amount that can not desorb crude petroleum materials in unit volume shale；K2For constant；PiFor oil reservoir initial pressure, P is production
The reservoir pressure in latter stage, CoFor the amount score of substance oily under original state, CgFor the amount score of the substance of gas under original state.
5. method for numerical simulation according to claim 4, which is characterized in that the amendment porosity is obtained by following formula
It arrives:
Wherein, φ1To correct porosity；φ0For primary porosity；ρk.oFor the mole close of the kerogen that can desorb and oily conjugate
Degree；ρk.gFor the molar density of the kerogen and gas conjugate that can desorb；ρk1For the kerogen and crude oil conjugate that can not be desorbed
Molar density.
6. method for numerical simulation according to claim 1, which is characterized in that the reaction rate is obtained by following formula
It arrives:
R=rk·(Ck.o+Ck.g)ek；
Wherein, r is reaction rate；rkFor reaction rate constant；Ek is the order of reaction；Ck.oKerogen for that can desorb is combined with oil
The molar concentration of object；Ck.gFor the molar concentration of the kerogen and gas conjugate that can desorb.
7. a kind of numerical simulation device for shale oil flowing, which is characterized in that described device includes:
First determining module, for determining the chemical reaction equation and chemical reaction equation hair of description shale oil seepage flow characteristics
Pressure condition when raw reaction；
Second determining module, for determining the amount of maximum adsorption crude petroleum materials in unit volume shale；
Third determining module, for according to the amount, dissolved gas oil ratio of maximum adsorption crude petroleum materials in the unit volume shale with
And condition of formation pressure, determine can be desorbed in unit volume shale oil substances amount, can stripping gas substance amount and can not desorb
The amount of crude petroleum materials；
4th determining module, for according to the amount that can desorb oil substances in the unit volume shale, the unit volume shale
In can stripping gas substance amount and the unit volume shale in can not desorb the amounts of crude petroleum materials, determine amendment porosity；
5th determining module, for determining the corresponding order of reaction of the chemical reaction equation and reaction rate；
Module is established, for according to maximum adsorption in the chemical reaction equation, the pressure condition, the unit volume shale
The amounts of oil substances can be desorbed in the amount of crude petroleum materials, the unit volume shale, can vent one's spleen in the unit volume shale substance
Amount, amount, the amendment porosity, the order of reaction and described that crude petroleum materials can not be desorbed in the unit volume shale
Reaction rate is based on reservoir numerical simulation software, establishes numerical simulator.
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910137987.4A CN109854236B (en) | 2019-02-25 | 2019-02-25 | Numerical simulation method and device for shale oil flow |
| PCT/CN2019/114056 WO2020173114A1 (en) | 2019-02-25 | 2019-10-29 | Method and device for simulating shale oil flow |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910137987.4A CN109854236B (en) | 2019-02-25 | 2019-02-25 | Numerical simulation method and device for shale oil flow |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109854236A true CN109854236A (en) | 2019-06-07 |
| CN109854236B CN109854236B (en) | 2022-06-14 |
Family
ID=66898834
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910137987.4A Active CN109854236B (en) | 2019-02-25 | 2019-02-25 | Numerical simulation method and device for shale oil flow |
Country Status (2)
| Country | Link |
|---|---|
| CN (1) | CN109854236B (en) |
| WO (1) | WO2020173114A1 (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2020173114A1 (en) * | 2019-02-25 | 2020-09-03 | 中国石油大学(华东) | Method and device for simulating shale oil flow |
| CN111912960A (en) * | 2020-08-24 | 2020-11-10 | 东北石油大学 | Quantitative evaluation method and system for kerogen free oil in shale |
| CN113670960A (en) * | 2021-07-23 | 2021-11-19 | 西南石油大学 | Real shale gas adsorption capacity prediction method based on molecular simulation |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114594531A (en) * | 2020-12-07 | 2022-06-07 | 中国石油化工股份有限公司 | Continental facies shale oil earthquake rock physical modeling method for movable oil |
| CN117269452B (en) * | 2023-08-18 | 2024-02-27 | 中国地质大学（北京） | Method, device and terminal for determining helium gas migration-dissolution ratio in tight sandstone reservoir |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103907114A (en) * | 2011-10-26 | 2014-07-02 | 兰德马克绘图国际公司 | Methods and systems of modeling hydrocarbon flow from kerogens in a hydrocarbon bearing formation |
| CN106547966A (en) * | 2016-11-01 | 2017-03-29 | 中国石油大学(华东) | A kind of shale oil adsorbance with can momentum evaluation model and its foundation, application process |
| WO2017062531A2 (en) * | 2015-10-09 | 2017-04-13 | Schlumberger Technology Corporation | Reservoir simulation using an adaptive deflated multiscale solver |
| CN107506948A (en) * | 2017-09-27 | 2017-12-22 | 科为联合(青岛)能源技术研究院有限公司 | A kind of shale oil and gas composite volume analysis method based on dynamic aerial drainage volume |
| CN108665086A (en) * | 2017-04-01 | 2018-10-16 | 中国石油化工股份有限公司 | A kind of shale gas well yield prediction computational methods |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104156593B (en) * | 2014-08-08 | 2017-04-19 | 中国石油大学(华东) | Shale oil and gas yield evaluation model building and parameter calibrating method applied under closed system |
| CN107818188B (en) * | 2016-09-14 | 2021-07-27 | 中国石油化工股份有限公司 | Simulation prediction method for in-situ exploitation of fluid heating oil shale |
| CN106501145A (en) * | 2016-09-18 | 2017-03-15 | 中国石油大学(北京) | The bearing calibration of shale gas reservoir numerical simulation |input paramete and device |
| CN107133411B (en) * | 2017-05-16 | 2019-06-25 | 中国石油大学(华东) | A kind of shale oil gas compositional modeling method and system considering oil reservoir mechanism |
| CN109854236B (en) * | 2019-02-25 | 2022-06-14 | 中国石油大学(华东) | Numerical simulation method and device for shale oil flow |
-
2019
- 2019-02-25 CN CN201910137987.4A patent/CN109854236B/en active Active
- 2019-10-29 WO PCT/CN2019/114056 patent/WO2020173114A1/en active Application Filing
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103907114A (en) * | 2011-10-26 | 2014-07-02 | 兰德马克绘图国际公司 | Methods and systems of modeling hydrocarbon flow from kerogens in a hydrocarbon bearing formation |
| WO2017062531A2 (en) * | 2015-10-09 | 2017-04-13 | Schlumberger Technology Corporation | Reservoir simulation using an adaptive deflated multiscale solver |
| CN106547966A (en) * | 2016-11-01 | 2017-03-29 | 中国石油大学(华东) | A kind of shale oil adsorbance with can momentum evaluation model and its foundation, application process |
| CN108665086A (en) * | 2017-04-01 | 2018-10-16 | 中国石油化工股份有限公司 | A kind of shale gas well yield prediction computational methods |
| CN107506948A (en) * | 2017-09-27 | 2017-12-22 | 科为联合(青岛)能源技术研究院有限公司 | A kind of shale oil and gas composite volume analysis method based on dynamic aerial drainage volume |
Non-Patent Citations (1)
| Title |
|---|
| 高诚 等: "油页岩原位开采数值模拟研究进展", 《石油钻采工艺》 * |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2020173114A1 (en) * | 2019-02-25 | 2020-09-03 | 中国石油大学(华东) | Method and device for simulating shale oil flow |
| CN111912960A (en) * | 2020-08-24 | 2020-11-10 | 东北石油大学 | Quantitative evaluation method and system for kerogen free oil in shale |
| CN111912960B (en) * | 2020-08-24 | 2022-11-08 | 东北石油大学 | Quantitative evaluation method and system for kerogen free oil in shale |
| CN113670960A (en) * | 2021-07-23 | 2021-11-19 | 西南石油大学 | Real shale gas adsorption capacity prediction method based on molecular simulation |
| CN113670960B (en) * | 2021-07-23 | 2021-12-28 | 西南石油大学 | Real shale gas adsorption capacity prediction method based on molecular simulation |
Also Published As
| Publication number | Publication date |
|---|---|
| WO2020173114A1 (en) | 2020-09-03 |
| CN109854236B (en) | 2022-06-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109854236A (en) | A kind of method for numerical simulation and device for shale oil flowing | |
| Shurr et al. | Unconventional shallow biogenic gas systems | |
| Tilley et al. | Gas isotope reversals in fractured gas reservoirs of the western Canadian Foothills: Mature shale gases in disguise | |
| Firoozabadi et al. | Prospects for subsurface CO2 sequestration | |
| Liu et al. | Study on competitive adsorption and displacing properties of CO2 enhanced shale gas recovery: advances and challenges | |
| Ma et al. | The formation mechanism of high-quality dolomite reservoir in the deep of Puguang Gas Field | |
| Yanze et al. | The role of diffusion for nonequilibrium gas injection into a fractured reservoir | |
| Bustin et al. | Learnings from a failed nitrogen enhanced coalbed methane pilot: Piceance Basin, Colorado | |
| Hauge et al. | Methane Production from Natural Gas Hydrates by CO2 Replacement–Review of Lab Experiments and Field Trial | |
| CN108194065B (en) | Method for improving shale oil huff and puff recovery ratio by using carbon dioxide/cosolvent system | |
| John et al. | Experimental and simulation of CO2-foam flooding in fractured chalk rock at reservoir conditions: Effect of mode of injection on oil recovery | |
| CN107387068A (en) | A kind of method and system for being used to determine shale gas reservoir free gas content | |
| Ancell et al. | Analysis of the coalbed degasification process at a seventeen well pattern in the Warrior Basin of Alabama | |
| He et al. | A field study on simulation of CO2 injection and ECBM production and prediction of CO2 storage capacity in unmineable coal seam | |
| Pashin | Coal as a petroleum source rock and reservoir rock | |
| Zhou et al. | Influence and sensitivity study of matrix shrinkage and swelling on enhanced coalbed methane production and CO2 sequestration with mixed gas injection | |
| WG et al. | An overview of underground hydrogen storage with prospects and challenges for the Australian context | |
| Luo et al. | A dual-porosity dual-permeability model for acid gas injection process evaluation in hydrogen-carbonate reservoirs | |
| Yang et al. | Simulation study of hydrogen sulfide removal in underground gas storage converted from the multilayered sour gas field | |
| Da Silva | Primary and enhanced recovery of Ekofisk Field: A single-and double-porosity numerical simulation study | |
| Maricic | Parametric and predictive analysis of horizontal well configurations for coalbed methane reservoirs in Appalachian Basin | |
| Zuta et al. | Oil recovery during CO2-foam injection in fractured chalk rock at reservoir conditions | |
| Parakh | Experimental investigation of enhanced coal bed methane recovery | |
| Kamruzzaman et al. | Enhanced Oil Recovery Experiments in Wolfcamp Outcrop Cores and Synthetic Cores to Assess Contribution of Pore-Scale Processes | |
| Tiyntayev et al. | Simulation-Based Evaluation of Concurrent CH4 Storage Potentials During CO2-ECBM in Karaganda Coal Basin |
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 |