CN115758851B - Method for selecting multi-scale propping agent for natural fracture-containing stratum fracture - Google Patents

Method for selecting multi-scale propping agent for natural fracture-containing stratum fracture Download PDF

Info

Publication number
CN115758851B
CN115758851B CN202211503074.8A CN202211503074A CN115758851B CN 115758851 B CN115758851 B CN 115758851B CN 202211503074 A CN202211503074 A CN 202211503074A CN 115758851 B CN115758851 B CN 115758851B
Authority
CN
China
Prior art keywords
fracture
propping agent
particle
crack
scale
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
Application number
CN202211503074.8A
Other languages
Chinese (zh)
Other versions
CN115758851A (en
Inventor
徐长贵
郭书生
梁豪
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China National Offshore Oil Corp CNOOC
CNOOC China Ltd Hainan Branch
Original Assignee
China National Offshore Oil Corp CNOOC
CNOOC China Ltd Hainan Branch
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China National Offshore Oil Corp CNOOC, CNOOC China Ltd Hainan Branch filed Critical China National Offshore Oil Corp CNOOC
Priority to CN202211503074.8A priority Critical patent/CN115758851B/en
Publication of CN115758851A publication Critical patent/CN115758851A/en
Application granted granted Critical
Publication of CN115758851B publication Critical patent/CN115758851B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A10/00TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
    • Y02A10/40Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping

Landscapes

  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention discloses a method for selecting a multi-scale propping agent for a natural fracture-containing stratum fracture, which comprises the steps of carrying out fracture expansion numerical simulation under the natural fracture condition according to the geomechanical condition of a reservoir stratum to obtain the fracture expansion condition of the natural fracture-containing stratum after fracturing; evaluating and dividing the main crack, the first-stage branch crack and the second-stage branch crack to obtain a multi-scale crack form; extracting the geometric type of the main crack according to the multi-scale crack morphology; establishing a multi-scale fracture solid-liquid two-phase flow numerical model; carrying out conveying numerical simulation by adopting propping agents with different particle sizes respectively to obtain the distribution condition of the propping agents; counting propping agent duty ratios of the main fracture, the first-stage branch fracture and the second-stage branch fracture under the condition of multi-scale fracture width according to propping agent distribution conditions; and the highest-efficiency propping agent filling particle sizes of the cracks with different dimensions are optimized according to the propping agent ratio in the main cracks, the primary branch cracks and the secondary branch cracks. The method can accurately obtain the optimal propping agent particle size of the multi-scale fracture.

Description

Method for selecting multi-scale propping agent for natural fracture-containing stratum fracture
Technical Field
The invention relates to a method for selecting a multi-scale propping agent for a natural fracture-containing stratum fracture, and belongs to the technical field of oil and gas well yield increase transformation.
Background
The method is characterized in that the fracture conductivity is supported under different proppant combination conditions based on single fracture expansion simulation analysis, and finally parameters such as the type, the particle size and the addition amount of the propping agent capable of meeting the full-purpose requirement of the oil reservoir are obtained through simulation by using fracturing software. For example, the patent "a non-conventional reservoir fracturing propping agent optimizing method and application" of China oil and gas stock company describes that a fracture image processing method is used for extracting fracture parameters, namely, a morphological refinement algorithm is used for acquiring the central axis of a fracture image, so that the fracture width is acquired. However, the main problem in the prior art is that the diversion capability of the combined propping agent is studied only for single cracks, and the optimal proportion of the combined propping agent is difficult to reach for a natural fracture reservoir due to the limited width of the cracks, so that the capability of entering the cracks with different widths of different particle sizes is inconsistent.
Disclosure of Invention
In order to overcome the problems in the prior art, the invention provides a method for selecting a multi-scale propping agent for a natural fracture-containing stratum fracture.
The technical scheme provided by the invention for solving the technical problems is as follows: a method for selecting a multi-scale propping agent for a natural fracture-containing stratum fracture, comprising the following steps:
step one, carrying out fracture expansion numerical simulation under natural fracture conditions according to the geomechanical conditions of the reservoir to obtain the fracture expansion condition of the natural fracture reservoir after fracturing;
step two, evaluating and dividing the main fracture, the first-stage branch fracture and the second-stage branch fracture according to the fracture expansion condition of the fractured natural fracture reservoir to obtain a multi-scale fracture morphology;
extracting the geometric type of the main crack according to the multi-scale crack shape;
step four, establishing a multi-scale fracture solid-liquid two-phase flow numerical model according to the geometric type of the main fracture;
step five, carrying out conveying numerical simulation by adopting propping agents with different particle sizes respectively to obtain propping agent distribution conditions;
step six, counting propping agent proportion of the main fracture, the first-stage branch fracture and the second-stage branch fracture under the condition of multi-scale fracture width according to propping agent distribution conditions;
and step seven, optimizing the highest-efficiency propping agent filling particle sizes of the cracks with different dimensions according to the propping agent ratio in the main cracks, the first-stage branch cracks and the second-stage branch cracks.
The further technical scheme is that the solid-liquid two-phase flow numerical model in the fourth step comprises:
(1) Equation of fluid control
Wherein: ρ l Is of liquid phase density of kg/m 3 ;ε l Liquid phase volume fraction,%; u (u) l Is the liquid phase flow speed, m/s; t is time, s; p is flow field pressure, pa; g is gravity acceleration, m/s 2 ;M pl Momentum exchange source items for the liquid phase and the particle phase; k is the turbulence energy of the fluid phase, m 2 /s 2 The method comprises the steps of carrying out a first treatment on the surface of the ε is the turbulent dissipation ratio, m 2 /s 3 ;G k kg/mS as a generating term of turbulent energy 3 ;σ k To turbulent energy pairThe plannd number of the response, dimensionless, is taken to be sigma k =1.0;σ ε Taking sigma for planner number corresponding to turbulence dissipation rate ε =1.3;S K 、S ε For turbulent flow exchange between liquid phase and solid phase, kg/m.s 3 ;C 、C Taking C as empirical constant and dimensionless =1.44、C =1.92;
(2) Equation of particle motion control
Wherein: m is m i The mass of the particles i is kg; u (u) p,i The linear velocity of the particles i, m/s; f (F) pc,ij A contact force N generated for the contact of the particle i with other particles; f (F) lp,i Acting force of fluid on the particles i, N; i pc,ij Kg.m for the moment of inertia of the particles i 2 ;ω p,i Rad/s is the angular velocity of particle i; t (T) pc,ij The contact moment generated by the contact of the particle i and the particle j is N.m;
(3) Fluid-particle force
F v =C v V p ρ l (Du slip /Dt)
Wherein: f (F) F Welfare for propping agent, N; v (V) p For proppant volume, m 3 ;F d N, the drag force of the fluid on the particulate phase(s); c (C) d As drag coefficient, dimensionless; epsilon p Is the volume fraction of the proppant, dimensionless; ρ l Is fluid density, kg/m 3 ;u p Is the movement speed of the propping agent, m/s; u (u) l Is the movement speed of the liquid phase, m/s; d is the particle size of the propping agent, m; f (F) v N, the proppant virtual mass force; u (u) slip Slip velocity, m/s, of the proppant relative to the continuous phase; c (C) v Is a virtual mass force coefficient, C V =0.5;
(4) Particle-particle force model
F c,ij =F cn,ij +F ct,ij
F cn,ij =k n,ij δ n,ij n+γ n,ij u n,ij
F ct,ij =k t,ij δ t,ij t+γ t,ij u t,ij
Wherein: f (F) c,ij N is the contact resultant force of the particles i and the particles j; f (F) cn,ij N, the normal contact force of particle i and particle j; f (F) ct,ij N, the tangential contact force of particle i and particle j; k (k) n,ij Is the normal stiffness coefficient; delta n,ij M is the normal displacement between the impinging particles; n is a normal unit vector between contact particles; gamma ray n,ij Is the normal dissipation coefficient; u (u) n,ij Is the normal component of the relative velocity between the impinging particles; k (k) t,ij Is a tangential stiffness coefficient; delta t,ij M is the normal displacement between the impinging particles; t is a tangential unit vector between contact particles; gamma ray t,ij Is the tangential dissipation factor; u (u) t,ij Is the tangential component of the relative velocity between impinging particles.
According to a further technical scheme, proppants with different particle sizes in the fifth step comprise 30/50 meshes, 40/70 meshes, 70/140 meshes and 100/200 meshes.
The technical scheme is that the proppant migration and transportation rule under the actual working condition is simulated in the fifth step.
According to a further technical scheme, the actual working conditions comprise displacement and fluid viscosity.
In the seventh step, the propping agent particle size corresponding to the optimal propping agent ratio of each crack in the main crack, the first-stage branch crack and the second-stage branch crack is respectively selected to be the propping agent filling particle size of the crack;
the filling degree is low;
filling degree is moderate;
the filling degree is good;
wherein: w (W) f1 The width of the first-stage crack is mm; w (W) f2 The width of the secondary crack is mm; v (V) P-f1 Is the volume of the first-stage fracture propping agent, m 3 ;V P-f2 Is the volume of the secondary fracture propping agent, m 3 The method comprises the steps of carrying out a first treatment on the surface of the α, β are coefficients related to the angle of the fracture, and are obtained through experiments, and α=0.089 and β=0.318 are taken as an example of a 90 ° orthogonal fracture.
The invention has the following beneficial effects: the method can accurately obtain the optimal propping agent particle size of the multi-scale fracture.
Drawings
FIG. 1 is a flow chart of the present invention;
FIG. 2 is a graph of the results of a natural fracture-containing reservoir fracture propagation simulation;
FIG. 3 is a statistical plot of multi-scale slit width distribution;
fig. 4 is a graph of proppant particle flow rate.
Detailed Description
The following description of the embodiments of the present invention will be made apparent and fully in view of the accompanying drawings, in which some, but not all embodiments of the invention are shown. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention.
As shown in fig. 1, the method for selecting the multi-scale proppants for the natural fracture-containing stratum fracture comprises the following steps:
step one, performing fracture expansion numerical simulation under a natural fracture condition according to the geomechanical condition of the reservoir to obtain the fracture expansion condition (shown in figure 2) of the natural fracture reservoir after fracturing;
step two, evaluating and dividing the main fracture, the primary branch fracture and the secondary branch fracture according to the fracture expansion condition after fracturing the natural fracture reservoir to obtain a multi-scale fracture morphology and multi-scale fracture width distribution (shown in figure 3);
extracting the geometric type of the main crack according to the multi-scale crack shape;
step four, establishing a multi-scale fracture solid-liquid two-phase flow numerical model according to the geometric type of the main fracture;
the solid-liquid two-phase flow numerical model comprises:
(1) Equation of fluid control
Wherein: ρ l Is of liquid phase density of kg/m 3 ;ε l Is the volume fraction of the liquid phase,%;u l Is the liquid phase flow speed, m/s; t is time, s; p is flow field pressure, pa; g is gravity acceleration, m/s 2 ;M pl Momentum exchange source items for the liquid phase and the particle phase; k is the turbulence energy of the fluid phase, m 2 /s 2 The method comprises the steps of carrying out a first treatment on the surface of the ε is the turbulent dissipation ratio, m 2 /s 3 ;G k kg/mS as a generating term of turbulent energy 3 ;σ k Taking sigma as planner number corresponding to turbulent energy k =1.0;σ ε Taking sigma for planner number corresponding to turbulence dissipation rate ε =1.3;S K 、S ε For turbulent flow exchange between liquid phase and solid phase, kg/m.s 3 ;C 、C Taking C as empirical constant and dimensionless =1.44、C =1.92;
(2) Equation of particle motion control
Wherein: m is m i The mass of the particles i is kg; u (u) p,i The linear velocity of the particles i, m/s; f (F) pc,ij A contact force N generated for the contact of the particle i with other particles; f (F) lp,i Acting force of fluid on the particles i, N; i pc,ij Kg.m for the moment of inertia of the particles i 2 ;ω p,i Is the angular velocity of particle i, rsd/s; t (T) pc,ij The contact moment generated by the contact of the particle i and the particle j is N.m;
(3) Fluid-particle force
F v =C v V p ρ l (Du slip /Dt)
Wherein: f (F) F Welfare for propping agent, N; v (V) p For proppant volume, m 3 ;F d N, the drag force of the fluid on the particulate phase(s); c (C) d As drag coefficient, dimensionless; epsilon p Is the volume fraction of the proppant, dimensionless; ρ l Is fluid density, kg/m 3 ;u p Is the movement speed of the propping agent, m/s; u (u) l Is the movement speed of the liquid phase, m/s; d is the particle size of the propping agent, m; f (F) v N, the proppant virtual mass force; u (u) slip Slip velocity, m/s, of the proppant relative to the continuous phase; c (C) v Is a virtual mass force coefficient, C V =0.5;
(4) Particle-particle force model
F c,ij =F cn,ij +F ct,ij
F cn,ij =k n,ij δ n,ij n+γ n,ij u n,ij
F ct,ij =k t,ij δ t,ij t+γ t,ij u t,ij
Wherein: f (F) c,ij N is the contact resultant force of the particles i and the particles j; f (F) cn,ij N, the normal contact force of particle i and particle j; f (F) ct,ij N, the tangential contact force of particle i and particle j; k (k) n,ij Is the normal stiffness coefficient; delta n,ij M is the normal displacement between the impinging particles; n is a normal unit vector between contact particles; gamma ray n,ij Is the normal dissipation coefficient; u (u) n,ij Is the normal component of the relative velocity between the impinging particles; k (k) t,ij Is a tangential stiffness coefficient; delta t,ij M is the normal displacement between the impinging particles; t is a tangential unit vector between contact particles; gamma ray t,ij For cuttingTo the dissipation factor; u (u) t,ij Tangential component of the relative velocity between impinging particles;
step five, propping agents with different particle sizes (30/50 meshes, 40/70 meshes, 70/140 meshes and 100/200 meshes) are adopted to simulate the migration and transportation rule of the propping agents under the condition of actual working conditions (displacement and fluid viscosity) so as to obtain the distribution condition of the propping agents;
step six, counting propping agent proportion of the main fracture, the primary branch fracture and the secondary branch fracture under the condition of multi-scale fracture width according to propping agent distribution conditions (shown in figure 4);
step seven, optimizing the highest-efficiency propping agent filling particle sizes of the cracks with different dimensions according to the propping agent occupation ratio in the main cracks, the first-stage branch cracks and the second-stage branch cracks;
in the seventh step, propping agent particle sizes corresponding to the optimal propping agent duty ratio of each crack in the main crack, the first-stage branch crack and the second-stage branch crack are respectively selected to be propping agent filling particle sizes of the crack;
the filling degree is low;
filling degree is moderate;
the filling degree is good;
wherein: w (W) f1 The width of the first-stage crack is mm; w (W) f2 The width of the secondary crack is mm; v (V) P-f1 Is the volume of the first-stage fracture propping agent, m 3 ;V P-f2 Is the volume of the secondary fracture propping agent, m 3 The method comprises the steps of carrying out a first treatment on the surface of the α, β are coefficients related to the angle of the fracture, and are obtained through experiments, taking 90 ° orthogonal fracture as an example, α=0.089, β=0.318;
and the propping agent grain size corresponding to the optimal propping agent ratio of each crack in the main crack, the first-stage branch crack and the second-stage branch crack is respectively selected to be the propping agent filling grain size of the crack.
The present invention is not limited to the above-mentioned embodiments, but is not limited to the above-mentioned embodiments, and any person skilled in the art can make some changes or modifications to the equivalent embodiments without departing from the scope of the technical solution of the present invention, but any simple modification, equivalent changes and modifications to the above-mentioned embodiments according to the technical substance of the present invention are still within the scope of the technical solution of the present invention.

Claims (4)

1. The method for selecting the multi-scale propping agent for the natural fracture-containing stratum fracture is characterized by comprising the following steps of;
step one, carrying out fracture expansion numerical simulation under natural fracture conditions according to the geomechanical conditions of the reservoir to obtain the fracture expansion condition of the natural fracture reservoir after fracturing;
step two, evaluating and dividing the main fracture, the first-stage branch fracture and the second-stage branch fracture according to the fracture expansion condition of the fractured natural fracture reservoir to obtain a multi-scale fracture morphology;
extracting the geometric type of the main crack according to the multi-scale crack shape;
step four, establishing a multi-scale fracture solid-liquid two-phase flow numerical model according to the geometric type of the main fracture;
step five, carrying out conveying numerical simulation by adopting propping agents with different particle sizes respectively to obtain propping agent distribution conditions;
the numerical model in the fifth step comprises the following steps:
(1) Equation of fluid control
Wherein: ρ l Is of liquid phase density of kg/m 3 ;ε l Liquid phase volume fraction,%; u (u) l Is the liquid phase flow speed, m/s; t is time, s; p is flow field pressure, pa; g is gravity acceleration, m/s 2 ;M pl Momentum exchange source items for the liquid phase and the particle phase; k is the turbulence energy of the fluid phase, m 2 /s 2 The method comprises the steps of carrying out a first treatment on the surface of the ε is the turbulent dissipation ratio, m 2 /s 3 ;G k kg/mS as a generating term of turbulent energy 3 ;σ k Taking sigma as planner number corresponding to turbulent energy k =1.0;σ ε Taking sigma for planner number corresponding to turbulence dissipation rate ε =1.3;S K 、S ε For turbulent flow exchange between liquid phase and solid phase, kg/m.s 3 ;C 、C Taking C as empirical constant and dimensionless =1.44、C =1.92;
(2) Equation of particle motion control
Wherein: m is m p,i The mass of the particles i is kg; u (u) p,i The linear velocity of the particles i, m/s; f (F) pc,ij A contact force N generated for the contact of the particle i with other particles; f (F) lp,i Acting force of fluid on the particles i, N; i pc,ij Kg.m for the moment of inertia of the particles i 2 ;ω p,i Rad/s is the angular velocity of particle i; t (T) pc,ij The contact moment generated by the contact of the particle i and the particle j is N.m;
(3) Fluid-particle force
F v =C v V p ρ l (Du slip /Dt)
Wherein: f (F) F The propping agent is subjected to buoyancy, N; v (V) p For proppant volume, m 3 ;F d N, the drag force of the fluid on the particle phase; c (C) d As drag coefficient, dimensionless; epsilon p Is the volume fraction of the proppant, dimensionless; ρ l Is fluid density, kg/m 3 ;u p Is the movement speed of the propping agent, m/s; u (u) l Is the movement speed of the liquid phase, m/s; d is the particle size of the propping agent, m; f (F) v N, the proppant virtual mass force; u (u) slip Slip velocity, m/s, of the proppant relative to the continuous phase; c (C) v Is a virtual mass force coefficient, C V =0.5;
(4) Particle-particle force model
F c,ij =F cn,ij +F ct,ij
F cn,ij =k n,ij δ n,ij n+γ n,ij u n,ij
F ct,ij =k t,ij δ t,ij t+γ t,ij u t,ij
Wherein: f (F) c,ij N is the contact resultant force of the particles i and the particles j; f (F) cn,ij N, the normal contact force of particle i and particle j; f (F) ct,ij N, the tangential contact force of particle i and particle j; k (k) n,ij Is the normal stiffness coefficient; delta n,ij M is the normal displacement between the impinging particles; n is a normal unit vector between contact particles; gamma ray n,ij Is the normal dissipation coefficient; u (u) n,ij Is the normal component of the relative velocity between the impinging particles; k (k) t,ij Is a tangential stiffness coefficient; delta t,ij M is the normal displacement between the impinging particles; t is a tangential unit vector between contact particles; gamma ray t,ij Is the tangential dissipation factor; u (u) t,ij Tangential component of the relative velocity between impinging particles;
step six, counting propping agent proportion of the main fracture, the first-stage branch fracture and the second-stage branch fracture under the condition of multi-scale fracture width according to propping agent distribution conditions;
step seven, optimizing the highest-efficiency propping agent filling particle sizes of the cracks with different dimensions according to the propping agent occupation ratio in the main cracks, the first-stage branch cracks and the second-stage branch cracks;
in the seventh step, propping agent particle sizes corresponding to the optimal propping agent duty ratio of each crack in the main crack, the first-stage branch crack and the second-stage branch crack are respectively selected to be propping agent filling particle sizes of the crack;
the filling degree is low;
filling degree is moderate;
the filling degree is good;
wherein: w (W) f1 The width of the first-stage crack is mm; w (W) f2 For the width of the secondary crack, mm;V P-f1 Is the volume of the first-stage fracture propping agent, m 3 ;V P-f2 Is the volume of the secondary fracture propping agent, m 3 The method comprises the steps of carrying out a first treatment on the surface of the α, β are coefficients related to the angle of the fracture, and are obtained through experiments, and α=0.089 and β=0.318 are taken as an example of a 90 ° orthogonal fracture.
2. The method of claim 1, wherein the proppants of different particle sizes in the fifth step comprise 30/50 mesh, 40/70 mesh, 70/140 mesh and 100/200 mesh.
3. The method for selecting the multi-scale propping agent for the natural fracture-containing stratum fracture according to claim 1, wherein the proppant migration and transportation rule under the actual working condition is simulated in the fifth step.
4. A method of selecting a multi-scale proppant for a fracture of a subterranean formation containing a natural fracture according to claim 3, wherein the actual conditions include displacement and fluid viscosity.
CN202211503074.8A 2022-11-28 2022-11-28 Method for selecting multi-scale propping agent for natural fracture-containing stratum fracture Active CN115758851B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211503074.8A CN115758851B (en) 2022-11-28 2022-11-28 Method for selecting multi-scale propping agent for natural fracture-containing stratum fracture

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211503074.8A CN115758851B (en) 2022-11-28 2022-11-28 Method for selecting multi-scale propping agent for natural fracture-containing stratum fracture

Publications (2)

Publication Number Publication Date
CN115758851A CN115758851A (en) 2023-03-07
CN115758851B true CN115758851B (en) 2024-01-05

Family

ID=85339456

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211503074.8A Active CN115758851B (en) 2022-11-28 2022-11-28 Method for selecting multi-scale propping agent for natural fracture-containing stratum fracture

Country Status (1)

Country Link
CN (1) CN115758851B (en)

Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107642348A (en) * 2016-07-21 2018-01-30 中国石油化工股份有限公司 A kind of method and its application for transforming crevice volume
CN108983312A (en) * 2018-05-25 2018-12-11 中国海洋石油集团有限公司 A kind of crack elimination anisotropic rock modulus evaluation method
CN110219635A (en) * 2018-03-01 2019-09-10 中国石油化工股份有限公司 A kind of full size fracture support method
CN110965980A (en) * 2020-01-17 2020-04-07 西南石油大学 Proppant conveying experimental device and method capable of obtaining proppant particle size distribution
CN111894545A (en) * 2019-05-06 2020-11-06 中国石油化工股份有限公司 Method for determining proppant pumping scheme
CN112580271A (en) * 2020-12-09 2021-03-30 中国石油大学(华东) Modeling method for integrating staged multi-cluster fracturing fracture expansion and proppant migration of horizontal well
CN113033049A (en) * 2021-03-22 2021-06-25 西南石油大学 Proppant conveying numerical simulation method in rough crack under stratum scale
CN113283197A (en) * 2021-06-10 2021-08-20 西南石油大学 Sand feeding parameter design method based on complex fracture width distribution
CN113530513A (en) * 2020-04-22 2021-10-22 中国石油化工股份有限公司 Fracturing method for graded support of proppants with different particle sizes in multi-scale fracture
CN114239308A (en) * 2021-12-24 2022-03-25 西南石油大学 Multi-scale high-density fracturing parameter optimization method
CN114458280A (en) * 2022-02-16 2022-05-10 西南石油大学 CFD-DEM (computational fluid dynamics-dynamic effect model) -based coated proppant conveying simulation method considering particle bonding
CN114611259A (en) * 2020-12-03 2022-06-10 中国石油化工股份有限公司 Method for evaluating effective reconstruction volume of multi-scale artificial seam
CN114622888A (en) * 2020-12-10 2022-06-14 中国石油化工股份有限公司 Complex fracture proppant dynamic migration rule testing method and application thereof

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9869649B2 (en) * 2015-09-03 2018-01-16 Saudi Arabian Oil Company Nano-level evaluation of kerogen-rich reservoir rock

Patent Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107642348A (en) * 2016-07-21 2018-01-30 中国石油化工股份有限公司 A kind of method and its application for transforming crevice volume
CN110219635A (en) * 2018-03-01 2019-09-10 中国石油化工股份有限公司 A kind of full size fracture support method
CN108983312A (en) * 2018-05-25 2018-12-11 中国海洋石油集团有限公司 A kind of crack elimination anisotropic rock modulus evaluation method
CN111894545A (en) * 2019-05-06 2020-11-06 中国石油化工股份有限公司 Method for determining proppant pumping scheme
CN110965980A (en) * 2020-01-17 2020-04-07 西南石油大学 Proppant conveying experimental device and method capable of obtaining proppant particle size distribution
CN113530513A (en) * 2020-04-22 2021-10-22 中国石油化工股份有限公司 Fracturing method for graded support of proppants with different particle sizes in multi-scale fracture
CN114611259A (en) * 2020-12-03 2022-06-10 中国石油化工股份有限公司 Method for evaluating effective reconstruction volume of multi-scale artificial seam
CN112580271A (en) * 2020-12-09 2021-03-30 中国石油大学(华东) Modeling method for integrating staged multi-cluster fracturing fracture expansion and proppant migration of horizontal well
CN114622888A (en) * 2020-12-10 2022-06-14 中国石油化工股份有限公司 Complex fracture proppant dynamic migration rule testing method and application thereof
CN113033049A (en) * 2021-03-22 2021-06-25 西南石油大学 Proppant conveying numerical simulation method in rough crack under stratum scale
CN113283197A (en) * 2021-06-10 2021-08-20 西南石油大学 Sand feeding parameter design method based on complex fracture width distribution
CN114239308A (en) * 2021-12-24 2022-03-25 西南石油大学 Multi-scale high-density fracturing parameter optimization method
CN114458280A (en) * 2022-02-16 2022-05-10 西南石油大学 CFD-DEM (computational fluid dynamics-dynamic effect model) -based coated proppant conveying simulation method considering particle bonding

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
基于最优支撑剂指数法优化低渗气藏裂缝参数;郭建春等;《西南石油大学学报(自然科学版)》;第35卷(第1期);93-98 *
泥灰岩储层缝网裂缝导流能力优化研究;段贵府;《中国优秀硕士学位论文全文数据库》(第2期);B019-12 *

Also Published As

Publication number Publication date
CN115758851A (en) 2023-03-07

Similar Documents

Publication Publication Date Title
Gomes et al. Effect of particle size and sphericity on the pickup velocity in horizontal pneumatic conveying
Yin et al. Thermophoresis and Brownian motion effects on nanoparticle deposition inside a 90 square bend tube
Dawen et al. Numerical investigations on the wake structures of micro-ramp and micro-vanes
CN115758851B (en) Method for selecting multi-scale propping agent for natural fracture-containing stratum fracture
Mishra et al. Thermal analysis of radiated (aluminum oxide)/water through a magnet based geometry subject to Cattaneo-Christov and Corcione’s models
Shi et al. Clusters and coherent voids in particle-laden wake flow
Da Silveira et al. Evaluation of collection efficiency methods for icing analysis
Hayati et al. Numerical modeling of particle motion and deposition in turbulent wavy channel flows
Wang et al. Drag model of finite-sized particle in turbulent wall-bound flow over sediment bed
Ernst et al. Quantification of preferential concentration of colliding particles in a homogeneous isotropic turbulent flow
Uchiyama et al. Numerical simulation for the collision between a vortex ring and solid particles
Guan et al. Numerical simulations of complex water spray flows around aircraft landing tire
Zhang et al. Numerical simulation of proppant transportation in hydraulic fracture based on DDPM-KTGF Model
Vaziri et al. Simulation of dynamic stall using direct-forcing immersed boundary method at low Reynolds number
Gu et al. Large-eddy simulation of the wind field and plume dispersion within different obstacle arrays using a dynamic mixing length subgrid-scale model
Munoz et al. Solid-state sintering simulation: surface, volume and grain-boundary diffusions
Yang et al. Influence of wall friction on granular column
Camata et al. FEM Simulation of coupled flow and bed morphodynamic interactions due to sediment transport phenomena
Damjanović et al. Car design as a new conceptual solution and CFD analysis in purpose of improving aerodynamics
Shen et al. The influence of particle density and diameter on the interactions between the finite-size particles and the turbulent channel flow
Smuts et al. A coupled CFD-DEM model for simulating the rheology of particulate suspensions
Barker et al. SPH Simulation of Packed-beds and Columns Applied to Heap-leaching
Zhou et al. Modeling and analysis of micro proppant transport and deposition in rough fracture networks
Tavouktsoglou et al. CFD simulation of clear water scour at complex foundations
Dritselis Direct numerical simulation of particle-laden turbulent channel flows with two-and four-way coupling effects: models of terms in the Reynolds stress budgets

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