CN112749468B - Numerical simulation method for capacity of solid-phase suspended matters to pass through pores of oil-gas reservoir - Google Patents
Numerical simulation method for capacity of solid-phase suspended matters to pass through pores of oil-gas reservoir Download PDFInfo
- Publication number
- CN112749468B CN112749468B CN201911048410.2A CN201911048410A CN112749468B CN 112749468 B CN112749468 B CN 112749468B CN 201911048410 A CN201911048410 A CN 201911048410A CN 112749468 B CN112749468 B CN 112749468B
- Authority
- CN
- China
- Prior art keywords
- pores
- fluid
- particle
- particles
- oil
- 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.)
- Active
Links
Images
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention provides a numerical simulation method for the passing capacity of solid phase suspended matters in pores of an oil and gas reservoir, which comprises the following steps: step 1, determining a numerical simulation research object of the pore of the oil-gas reservoir according to research needs, and establishing boundary parameters; step 2, establishing a pore simulation model of the oil-gas reservoir according to the boundary parameters of the research object; step 3, determining fluid parameters in pores according to the basic model obtained in the step 2, and simulating the fluid in the oil-gas reservoir by adopting a computational fluid mechanics method; step 4, determining parameters of particles in pores, and calculating solid suspended matters in the oil-gas reservoir by adopting a discrete element method; and 5, coupling the step 3 and the step 4 to simulate the capacity of solid phase suspended matters to pass through pores of the oil and gas reservoir. The method for simulating the passing capacity of the solid-phase suspended matters in the pores of the oil-gas reservoir can intuitively obtain an experimental effect and can solve the technical problems that an experimental device is difficult to manufacture, experimental conditions are difficult to change and the like.
Description
Technical Field
The invention relates to the technical field of oilfield development, in particular to a numerical simulation method for the passing capacity of solid phase suspended matters in pores of an oil-gas reservoir stratum.
Background
The solid phase suspended matter injected into water blocks pores and throats in oil and gas reservoirs according to certain rules, and oil layer protection can be guided by carefully researching the conditions under which the blocking phenomena occur. For example, foreign scholars Barkman and dawdson set forth the well-known one-third to one-seventh laws in studying mud filter cakes. That is, when the particle size of the solid phase suspended matter is larger than one third of the diameter of the throat of the rock core, the solid phase suspended matter injected into water can not enter an oil layer, and only a filter cake can be formed on the outer surface of a stratum, so that a certain blocking effect is realized on the section of the stratum; when the particle size of the solid phase suspended matter injected into the water is smaller than one seventh of the diameter of the throat of the rock core, the solid phase suspended matter cannot block the stratum; when the particle size of the solid phase suspended matter in the injected water is larger than one third of the diameter of the core throat and smaller than one seventh of the diameter of the core throat, the solid phase suspended matter in the injected water can cause serious damage to the stratum, so that stratum pores are blocked, and the stratum permeability is reduced. After that, other scholars have obtained the compatibility principle of one third to one tenth and one third to one fourteen first-class solid suspension injected into water and the diameter of underground city throat by means of experiments.
However, the method for obtaining the matching relationship between the particle size of the solid-phase suspended matter in the injected water and the pore radius of the oil-gas reservoir by means of indoor experiments is complicated, and in some cases, the method does not have the necessary condition for carrying out the indoor experiments. Therefore, a novel numerical simulation method for the capacity of solid phase suspended matters passing through the pores of the oil-gas reservoir is invented, and the technical problems are solved.
Disclosure of Invention
The invention aims to provide a numerical simulation method for the capacity of solid-phase suspended matters to pass through pores of a hydrocarbon reservoir.
The object of the invention can be achieved by the following technical measures: a method for numerically simulating the ability of a solid phase suspension to pass through hydrocarbon reservoir pores, the method comprising: step 1, determining a numerical simulation research object of the pore of the oil-gas reservoir according to research needs, and establishing boundary parameters; step 2, establishing a pore simulation model of the oil-gas reservoir according to the boundary parameters of the research object; step 3, determining fluid parameters in pores according to the basic model obtained in the step 2, and simulating the fluid in the oil-gas reservoir by adopting a computational fluid mechanics method; step 4, determining parameters of particles in pores on the basis of the simulation model obtained in the step 2, and calculating solid suspended matters in the oil-gas reservoir by adopting a discrete element method; and 5, coupling and simulating the capacity of solid-phase suspended matters in the pores of the oil and gas reservoir by the coupling of the step 3 and the step 4.
The object of the invention can also be achieved by the following technical measures:
the step 1 also comprises the step of forming data information of different grain diameters and throat ratios according to geological conditions.
In step 2, the step size is calculated through discrete iterative optimization, and the throat part at the contraction part is reduced by the operation step size.
In step 3, determining the fluid calculation type through the calculation model obtained in step 2: transient calculation and a gravity acceleration value; selecting a fluid viscosity related value; determining fluid properties includes viscosity, density.
In step 4, the appropriate solid suspension parameters are selected: particle density, particle size, concentration, poisson's ratio, and interparticle interaction forces.
In step 4, because the pore throat is influenced by solid phase suspended matters, a corrected Navier-Stokes equation is adopted for calculation, and a continuity equation and a kinetic energy equation are respectively given:
where ρ is f Fluid density, kg/m 3 (ii) a t is time, s; alpha is the volume fraction of the fluid phase, and has no dimension; v, fluid velocity, m/s; μ, fluid viscosity, mPa · s;
is the fluid pressure drop, Pa; f. of b Is the volume force, N, experienced by the fluid.
In step 4, tracking calculation is carried out on the migration of the single particles through a discrete element method, and interaction between the particles is simulated; respectively tracking the speed change of each particle during the free motion in the system by integrating Newton's second motion law, and coupling the acting force exerted by the fluid on the particles; the contact force between the particles is calculated by:
in the formula, F i Is the contact force between the particles, N; k n Normal contact stiffness, N/m; k s Shear contact stiffness, N/m; delta n Is the normal action distance, m, between the particles; 6, delta s Is the normal action distance, m, between the particles; gamma ray n Is a critical normal damping ratio; gamma ray s Is the critical shear damping ratio; m is the mass of the particles, kg; v n The normal relative velocity of the particles is m/s; v s Is the relative velocity of particle shear, m/s; i is a particle label, i is 1 or 2;
the mechanical interaction between the particles and the contact shear will cause them to rotate, and the particle momentum is calculated by the following equation:
M i ＝K s δ s (x 0 -x i )
in the formula, M i The momentum at the point of particle contact, kg · m/s; x is a radical of a fluorine atom 0 Is the position of the particle contact point; x is the number of i Is the geometric center position of the particle i.
In step 3 and step 4, based on the basic model in step 2, the parameters of the fluid and the solid suspended matter in the pores need to be determined, and the determined parameters of the pore throat comprise: establishing a model according to the length of the pore, the diameter of the pore and the diameter of the throat at the contraction position of the pore; the solid suspension parameters determined include: particle size, density, roundness; and the experiments were numbered with different particle size-throat ratios.
In step 3 and step 4, based on the mathematical model obtained in step 2, a model boundary condition is established: ingress, egress and other boundary conditions and ingress velocity values.
In step 5, simulating the capacity of solid phase suspended solids to pass through pores of the oil and gas reservoir through pressure and speed coupling of the fluid and particle parameters calculated in the step 3 and the step 4; firstly, calculating the drag force of fluid to particles in each calculation unit, then calculating the interaction force among the particles and between the particles and the wall surface according to a force-displacement rule, and finally calculating the particle migration according to a Newton's second law; the porosity in the computing unit is changed due to the particle migration, so that the flow field is influenced; calculating and updating the flow field by calculating the fluid pressure and velocity in a fluid dynamics method between the force-displacement calculation and the motion calculation in the discrete element method;
the drag force exerted by the fluid on each particle is as follows:
wherein, the first and the second end of the pipe are connected with each other,
is the drag force exerted on each particle; d is the particle size;
for particles on each cellAverage volume force of (2).
In step 5, the best matching particle size-pore throat ratio is obtained by analyzing the results obtained from the coupled simulation of step 3 and step 4.
The numerical simulation method for the passing capacity of the solid-phase suspended matters in the pores of the oil-gas reservoir can simulate the distribution state and the filling state of different solid suspended matters under the conditions of different particle sizes and pore-throat ratios. The experimental effect can be intuitively obtained. The technical problems that an experimental device is difficult to manufacture, experimental conditions are difficult to change and the like can be solved.
Drawings
FIG. 1 is a flow chart of one embodiment of a method of numerical simulation of the ability of a solid phase suspension to pass through the pores of a hydrocarbon reservoir of the present invention;
FIG. 2 is a diagram of a pore throat model in an embodiment of the present invention;
FIG. 3 is a diagram of simulation results in an embodiment of the present invention.
Detailed Description
In order to make the aforementioned and other objects, features and advantages of the present invention comprehensible, preferred embodiments accompanied with figures are described in detail below.
As shown in FIG. 1, FIG. 1 is a flow chart of a method of numerical simulation of the ability of a solid phase suspension of the present invention to pass through the pores of a hydrocarbon reservoir.
102, establishing a hydrocarbon reservoir pore simulation model according to the boundary parameters of the research object; establishing a model boundary condition: ingress, egress and other boundary conditions and ingress velocity values. The model is divided by adopting an unstructured grid, the step length is calculated by discrete iterative optimization, and the operation step length is reduced at the throat part at the contraction part (as shown in figure 2).
103, determining parameters of the fluid in the pores on the basis of the basic model obtained in the step 102, and simulating the fluid in the oil-gas reservoir by adopting a Computational Fluid Dynamics (CFD) method;
from the calculation model obtained in step 102, the fluid calculation type is determined: transient calculation and a gravity acceleration value; selecting a fluid viscosity related value; determining fluid properties includes viscosity, density, and the like.
Because the pore throat is influenced by solid phase suspended matters, the modified Navier-Stokes equation is adopted for calculation, and a continuity equation and a kinetic energy equation are respectively given.
Where ρ is f Fluid density, kg/m 3 (ii) a t is time, s; alpha is the volume fraction of the fluid phase, and has no dimension; v, fluid velocity, m/s; μ, fluid viscosity, mPa · s;
is the fluid pressure drop, Pa; f. of b Is the volume force, N, experienced by the fluid.
The migration of single particles can be tracked and calculated by a discrete element method, and the interaction between the particles can be simulated. By integrating newton's second law of motion, the velocity change of each particle during free motion in the system (between successive collisions) can be tracked separately and the force exerted by the fluid on the particles coupled. The contact force between the particles (the elastic force and the damping force between the particles) can be calculated by the following equation.
In the formula, F i Is the contact force between the particles, N; k n Normal contact stiffness, N/m; k s Shear contact stiffness, N/m; delta n Is the normal action distance, m, between the particles; 6, delta s Is the normal action distance, m, between the particles; gamma ray n Is a critical normal damping ratio; gamma ray s Is the critical shear damping ratio; m is the mass of the particles, kg; v n The normal relative velocity of the particles is m/s; v s Is the relative velocity of particle shear, m/s; i is a particle number, i is 1 and 2.
The mechanical interaction between the particles and the contact shear will cause them to rotate, and the particle momentum can be calculated from the following equation.
M i ＝K s δ s (x 0 -x i )
In the formula, M i The momentum at the point of particle contact, kg · m/s; x is the number of 0 Is the position of the particle contact point; x is the number of i Is the geometric center position of the particle i.
Based on the mathematical model obtained in step 102, model boundary conditions are established in steps 103 and 104: ingress, egress and other boundary conditions and ingress velocity values.
In step 103 and step 104, based on the basic model in step 102, parameters of the fluid and the suspended solid in the pores are determined, and the determined parameters of the pore throat include: establishing a model according to the length of the pore, the diameter of the pore and the diameter of the throat at the contraction position of the pore; the solid suspension parameters determined include: particle size, density, roundness; and the experiments were numbered with different particle size-throat ratios.
The fluid and particle parameters calculated in step 103 and step 104 are coupled by pressure and velocity to simulate the ability of solid phase suspensions to pass through the pores of a hydrocarbon reservoir. The drag force of the fluid on the particles is calculated in each calculation unit, then the interaction force between the particles and the wall surface is calculated according to a force-displacement principle, and finally the particle migration is calculated according to Newton's second law. The particle migration causes the porosity in the computing unit to change, thereby affecting the flow field. And calculating and updating the flow field by calculating the fluid pressure and the fluid velocity in a fluid dynamics method between the force-displacement calculation and the motion calculation in the discrete element method.
The drag force exerted by the fluid on each particle is as follows.
Wherein the content of the first and second substances,
is the drag force exerted on each particle; d is the particle size;
is the average volumetric force experienced by the particle on each cell.
And analyzing the results obtained by coupling simulation in the step 103 and the step 104 to obtain the best matching particle size-pore throat ratio.
Therefore, simulation can be carried out on the matching relation between the particle size of solid phase suspended matters in the injected water and the pore radius of the oil and gas reservoir through a means of coupling a fluid dynamics model and a discrete element model.
Claims (11)
1. A method for numerically simulating the ability of a solid phase suspension to pass through the pores of a hydrocarbon reservoir, the method comprising:
step 1, determining a numerical simulation research object of the pore of the oil-gas reservoir according to research needs, and establishing boundary parameters;
step 2, establishing a pore simulation model of the oil-gas reservoir according to the boundary parameters of the research object;
step 3, determining fluid parameters in pores according to the basic model obtained in the step 2, and simulating the fluid in the oil-gas reservoir by adopting a computational fluid mechanics method;
step 4, determining parameters of particles in pores on the basis of the simulation model obtained in the step 2, and calculating solid suspended matters in the oil-gas reservoir by adopting a discrete element method;
and 5, coupling and simulating the capacity of solid-phase suspended matters in the pores of the oil and gas reservoir by the coupling of the step 3 and the step 4.
2. The method for numerical simulation of the ability of a solid suspension to pass through pores of a hydrocarbon reservoir of claim 1, wherein step 1 further comprises generating data information of different particle size-throat ratio according to geological conditions.
3. The method for numerical simulation of the transit ability of a solid phase suspension in hydrocarbon reservoir pores according to claim 1, wherein in step 2, the calculation step size is calculated by discrete iterative optimization, and the calculation step size is reduced at the throat part at the constriction.
4. The method for numerical simulation of the ability of a solid phase suspension to pass through the pores of a hydrocarbon reservoir as set forth in claim 1, wherein in step 3, the calculation type of fluid is determined from the calculation model obtained in step 2: transient calculation and a gravity acceleration value; selecting a fluid viscosity related value; determining fluid properties including viscosity, density.
5. The method of numerical simulation of the ability of a solid suspension to pass through hydrocarbon reservoir pores of claim 1, wherein in step 4, suitable solid suspension parameters are selected: particle density, particle size, concentration, poisson's ratio, and interparticle interaction forces.
6. The method for numerical simulation of the ability of solid phase suspensions to pass through the pores of an oil and gas reservoir according to claim 1, wherein in step 4, the modified Navier-Stokes equation is used for calculation due to the influence of the solid phase suspensions on the pore throats to respectively give a continuity equation and a kinetic energy equation:
wherein ρ f Fluid density, kg/m 3 (ii) a t is time, s; alpha is the volume fraction of the fluid phase, and has no dimension; v, fluid velocity, m/s; μ, fluid viscosity, mPa · s;
is the fluid pressure drop, Pa; f. of b Is the volume force, N, experienced by the fluid.
7. The numerical simulation method for the capacity of solid phase suspensions to pass through the pores of oil and gas reservoirs according to claim 6, characterized in that in step 4, the migration of single particles is tracked and calculated by a discrete element method, and the interaction between the particles is simulated; respectively tracking the speed change of each particle during the free motion in the system by integrating Newton second motion law, and coupling the acting force exerted by the fluid on the particles; the contact force between the particles is calculated by:
in the formula, F i Is the contact force between the particles, N; k is n Normal contact stiffness, N/m; k s Shear contact stiffness, N/m; delta n Is the normal action distance, m, between the particles; 6, delta s Is the normal action distance, m, between the particles; gamma ray n Is a critical normal damping ratio; gamma ray s Is the critical shear damping ratio; m is the mass of the particles, kg; v n The normal relative velocity of the particles is m/s; v s Is the relative velocity of particle shear, m/s; i is a particle label, i is 1 or 2;
the mechanical interaction between the particles and the contact shear will cause them to rotate, and the particle momentum is calculated by the following equation:
M i ＝K s δ s (x 0 -x i )
in the formula, M i The momentum at the point of particle contact, kg · m/s; x is the number of 0 Is the position of the particle contact point; x is the number of i Is the geometric center position of the particle i.
8. The method of numerical simulation of the ability of a solid suspension to pass through the pores of a hydrocarbon reservoir of claim 1 wherein in steps 3 and 4, based on the base model of step 2, parameters of the fluid and solid suspension in the pores are determined, the determined parameters of the pore throat comprising: establishing a model according to the length of the pore, the diameter of the pore and the diameter of the throat at the contraction position of the pore; the solid suspension parameters determined include: particle size, density, roundness; and the experiments were numbered with different particle size-throat ratios.
9. The method for numerically simulating the ability of a solid phase suspension to pass through the pores of a hydrocarbon reservoir as claimed in claim 1, wherein in step 3 and step 4, model boundary conditions are established based on the mathematical model obtained in step 2: ingress, egress and other boundary conditions and ingress velocity values.
10. The method for numerical simulation of the ability of a solid phase suspension to pass through hydrocarbon reservoir pores according to claim 1, wherein in step 5, the fluid and particle parameters calculated in steps 3 and 4 are used to simulate the ability of a solid phase suspension to pass through hydrocarbon reservoir pores by pressure-velocity coupling; firstly, calculating the drag force of fluid on particles in each calculation unit, then calculating the interaction force between the particles and the wall surface according to a force-displacement rule, and finally calculating the particle migration according to a Newton second law; the porosity in the computing unit is changed due to the particle migration, so that the flow field is influenced; calculating and updating the flow field by calculating the fluid pressure and velocity in a fluid dynamics method between the force-displacement calculation and the motion calculation in the discrete element method;
the drag force exerted by the fluid on each particle is as follows:
wherein the content of the first and second substances,
is the drag force exerted on each particle; d is the particle size;
is the average volumetric force experienced by the particle on each cell.
11. The method of numerical simulation of the ability of a solid suspension to pass through the pores of a hydrocarbon reservoir of claim 10 wherein in step 5 the best match particle size to pore throat ratio is determined by analyzing the results of the coupled simulation of step 3 and step 4.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911048410.2A CN112749468B (en) | 2019-10-29 | 2019-10-29 | Numerical simulation method for capacity of solid-phase suspended matters to pass through pores of oil-gas reservoir |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911048410.2A CN112749468B (en) | 2019-10-29 | 2019-10-29 | Numerical simulation method for capacity of solid-phase suspended matters to pass through pores of oil-gas reservoir |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112749468A CN112749468A (en) | 2021-05-04 |
| CN112749468B true CN112749468B (en) | 2022-09-02 |
Family
ID=75640998
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201911048410.2A Active CN112749468B (en) | 2019-10-29 | 2019-10-29 | Numerical simulation method for capacity of solid-phase suspended matters to pass through pores of oil-gas reservoir |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112749468B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113188976B (en) * | 2021-05-11 | 2022-08-16 | 中国石油大学(华东) | Method and system for determining anisotropic permeability of sandwich-shaped shale |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107291993A (en) * | 2017-05-27 | 2017-10-24 | 中国石油大学(华东) | The analogy method of pre-crosslinked gel suspension micro flow in a kind of porous media |
| CN109657322A (en) * | 2018-12-13 | 2019-04-19 | 中国科学院、水利部成都山地灾害与环境研究所 | A kind of solid-liquid multiphase is suitable for the Numerical Simulation On The Dynamics of mud-rock flow |
-
2019
- 2019-10-29 CN CN201911048410.2A patent/CN112749468B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107291993A (en) * | 2017-05-27 | 2017-10-24 | 中国石油大学(华东) | The analogy method of pre-crosslinked gel suspension micro flow in a kind of porous media |
| CN109657322A (en) * | 2018-12-13 | 2019-04-19 | 中国科学院、水利部成都山地灾害与环境研究所 | A kind of solid-liquid multiphase is suitable for the Numerical Simulation On The Dynamics of mud-rock flow |
Non-Patent Citations (1)
| Title |
|---|
| 沉积盆地地层孔隙动力学研究进展;刘震等;《石油学报》;20161025(第10期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112749468A (en) | 2021-05-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Ju et al. | Enhanced oil recovery by flooding with hydrophilic nanoparticles | |
| Choi et al. | Fluidization by lift of 300 circular particles in plane Poiseuille flow by direct numerical simulation | |
| Chang et al. | Physics of proppant transport through hydraulic fracture network | |
| CN104021277A (en) | Numerical analysis method for piping phenomenon | |
| Bonnecaze et al. | Patterns of sedimentation from polydispersed turbidity currents | |
| CN112199903B (en) | Multi-parameter-based discrete element nanoparticle plugging shale pore numerical simulation method | |
| CN106951641B (en) | A kind of method and system of fracture-pore reservoir numerical simulation | |
| Yang et al. | Nanoparticle plugging prediction of shale pores: A numerical and experimental study | |
| CN112749468B (en) | Numerical simulation method for capacity of solid-phase suspended matters to pass through pores of oil-gas reservoir | |
| CN115510778A (en) | Continental facies shale reservoir infinite stage fracturing process optimization method and system | |
| Xu et al. | Influence of reservoir water level variations on slope stability and evaluation of landslide tsunami | |
| Yin et al. | Clogging and unclogging of fine particles in porous media: Micromechanical insights from an analog pore system | |
| CN103206203B (en) | The analytical method that the single perforation of oil well shakes out | |
| CN111259564B (en) | Method for predicting sand production rule of hydrate pre-filled sand control screen pipe | |
| Kazidenov et al. | Coarse-graining methods for the modified jkr contact model on a triaxial compression test | |
| CN114638142A (en) | Prediction method for borehole abrasion and actual diameter in close-cut multi-cluster temporary plugging fracturing | |
| Yang et al. | Method for calculating pressure losses in the pipelines of slurry shield tunneling based on coupled simulation of computational fluid dynamics and discrete element method | |
| CN109505591B (en) | Method and system for determining permeability limit of unfilled karst cave of fracture-cavity oil reservoir | |
| Lo et al. | Channeling and stress during fluid and suspension flow in self-affine fractures | |
| CN107609257B (en) | Gravity flow slumping body deposition simulation method based on numerical simulation | |
| Xia et al. | A Bonded 3D LBM-DEM Approach for Sand Production Process | |
| Liu et al. | Particle migration behavior of fractured lost circulation control: Numerical simulation | |
| Dun et al. | Effects of shape and rotation of sand particles in saltation | |
| Li et al. | Flow behaviors of ellipsoidal suspended particles in porous reservoir rocks using CFD-DEM combined with multi-element particle model | |
| CN117669411A (en) | Gravel packing sand prevention simulation method based on CFD-DEM |
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 |