CN116306332A - Particle impact rock breaking process simulation optimization method and system - Google Patents
Particle impact rock breaking process simulation optimization method and system Download PDFInfo
- Publication number
- CN116306332A CN116306332A CN202211548562.0A CN202211548562A CN116306332A CN 116306332 A CN116306332 A CN 116306332A CN 202211548562 A CN202211548562 A CN 202211548562A CN 116306332 A CN116306332 A CN 116306332A
- Authority
- CN
- China
- Prior art keywords
- rock
- particle
- impact
- model
- particles
- 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.)
- Pending
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/25—Design optimisation, verification or simulation using particle-based methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/14—Pipes
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- Fluid Mechanics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Physics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computing Systems (AREA)
- Algebra (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention provides a particle impact rock breaking process simulation optimization method and system, which adopts computational fluid dynamics to simulate the flow of drilling fluid in a three-dimensional pipeline and adopts a digital elevation model to simulate the flow of impact particles in the drilling fluid and the process of spraying the impact particles on a rock mass; the rock mass breaking degree is characterized based on the breaking of parallel bonding bonds in the rock stratum model; representing the wear range of the pipeline based on the stress distribution condition in the spray pipe; the blocking condition in the spray pipe is represented based on the particle accumulation position in the spray pipe, and the damage effect of the particle impact rock breaking process on the rock mass under the current simulation parameters is obtained; and determining a particle impact rock breaking process optimization scheme which is suitable for actual operation conditions based on the analysis result. The invention can overcome the defect of the existing particle impact rock breaking simulation method in parameter characterization, optimize the particle impact rock breaking process and improve the production efficiency.
Description
Technical Field
The invention belongs to the technical field of high-speed particle impact rock breaking simulation, and relates to a particle impact rock breaking process simulation optimization method and system.
Background
The statements in this section merely provide background information related to the present disclosure and may not necessarily constitute prior art.
In the early exploitation process of petroleum, natural gas and other energy sources, the deep rock body is hard, so that the exploitation range is limited in the shallow rock layer, with the continuous development of the industrial field, the oil and gas resources of the shallow earth surface gradually cannot meet the industrial requirements, and the exploitation of the oil and gas resources of the deep stratum becomes an important work in the oil and gas exploration field. But during the production process, the hardness of the rock stratum and the production difficulty increase along with the increase of the depth of the produced stratum. Thus, how to increase the rate of rock breaking in hard formations is a current challenge.
The existing rock breaking process can be divided into three types of hydraulic rock breaking, mechanical rock breaking and particle impact rock breaking according to mechanisms, wherein the hydraulic rock breaking refers to a technology for breaking rock by applying high-pressure water jet flow exceeding 20MPa to break the rock so as to break the rock; mechanical rock breaking refers to the tunneling in the stratum by using a special cutting bit; the particle impact rock breaking means that particles are injected into high-speed liquid through a particle injection system after drilling fluid is pumped out, the particles are conveyed into a drill bit at the bottom of a well by a high-pressure pipeline, the particles are finally accelerated to the speed of the drilling fluid, then are sprayed out from a special PID drill bit nozzle, impact the stratum to achieve the purpose of breaking the rock, and finally, the intact particles are separated from the drilling fluid returned from a wellhead to achieve recycling of the particles. The drilling speed of the particle impact rock breaking is generally 3-4 times of that of the conventional rock breaking, dust is not generated, centralized control is easy to realize, the drilling period is greatly shortened, and the particle impact rock breaking has great potential in the application of deep hard stratum.
With the progress of computer technology and the perfection of computing theory, numerical simulation technology has become one of the mainstream means for researching engineering problems. The traditional particle impact rock breaking simulation is mostly based on a finite element method, a finite difference method and the like, and in actual engineering, a complex two-phase flow system is formed by high-speed particles and high-speed drilling fluid, interaction between the drilling fluid and impact particles cannot be simulated by a traditional single simulation calculation method, and theoretical guidance on rock breaking process parameters is difficult to be provided.
Disclosure of Invention
In order to solve the problems, the invention provides a particle impact rock breaking process simulation optimization method and system, which can overcome the defects of the existing particle impact rock breaking simulation method in the aspects of impact particle size distribution, drilling fluid flow rate, pipeline pressure, particle volume fraction and other parameter characterization, optimize the particle impact rock breaking process and improve the production efficiency.
According to some embodiments, the present invention employs the following technical solutions:
a particle impact rock breaking process simulation optimization method comprises the following steps:
acquiring space data of a rock stratum, acquiring structural surface information of the rock stratum, and establishing a model for the rock stratum;
acquiring macroscopic mechanical parameters of a rock sample at the bottom of a target well;
filling particles into an established rock stratum model according to a particle grading curve, transmitting force and moment between the particles through a parallel bonding model, and carrying out parameter calibration according to macroscopic mechanical parameters to ensure that the strength index of the generated rock stratum model is consistent with that of an obtained rock sample;
according to a rock breaking spray pipe used in a working place, a three-dimensional pipeline model is established;
simulating the flow of drilling fluid in a three-dimensional pipeline by adopting computational fluid dynamics, and simulating the flow of impact particles in the drilling fluid and the process of spraying the impact particles on a rock mass by adopting a digital elevation model;
the rock mass breaking degree is characterized based on the breaking of parallel bonding bonds in the rock stratum model; representing the wear range of the pipeline based on the stress distribution condition in the spray pipe; the blocking condition in the spray pipe is represented based on the particle accumulation position in the spray pipe, and the damage effect of the particle impact rock breaking process on the rock mass under the current simulation parameters is obtained;
changing rock breaking parameters, carrying out multiple simulation processes, calculating response relations between various rock breaking parameters and simulation results, and determining a particle impact rock breaking process optimization scheme suitable for actual operation conditions based on analysis results.
As an alternative embodiment, the macro-mechanical parameters include several of particle size distribution curve, uniaxial compressive strength, tensile strength, modulus of elasticity, poisson's ratio, friction angle, cohesion and coefficient of friction.
As an alternative implementation mode, the concrete process of establishing the three-dimensional pipeline model comprises the steps of establishing the three-dimensional model according to the same size ratio according to a spray pipe used for on-site drilling operation, and setting the length from an inlet to an outlet of the spray pipe;
and filling the wall body units into the three-dimensional model, and setting the distance from the outlet of the matched spray pipe to the rock mass according to the length from the inlet to the outlet of the spray pipe.
As an alternative embodiment, the process of simulating the flow of drilling fluid in a three-dimensional pipeline by adopting computational fluid dynamics and simulating the flow of impact particles in the drilling fluid and spraying the impact particles on a rock mass by adopting a digital elevation model specifically comprises the following steps:
setting parameters and boundary conditions for calculating the coupling of the fluid dynamics-digital elevation model;
dividing proper fluid calculation grids at the same spatial position of the pipeline according to the spatial size of the established three-dimensional pipeline model, and simulating the flow of drilling fluid in the pipeline;
generating impact particle particles at the inlet of the pipeline according to the volume fraction according to the set parameters of the coupling of the computational fluid dynamics-digital elevation model, setting the generation positions of the particles in the x direction, setting contact models between the impact particles and the impact particles, between the impact particles and the pipeline wall, and between the impact particles and the rock stratum to be linear contact models, and simulating the physical process of the impact particles in the rock breaking process;
setting a fluid calculation domain time step and a digital elevation model calculation domain time step, and calculating according to the set step until the total simulation time is met.
As a further example, the parameters include impact particle diameter, particle density, particle poisson's ratio, particle static friction system, particle rolling friction coefficient, particle generation rate, particle recovery coefficient, particle volume fraction, and conduit wall friction coefficient.
Boundary conditions include setting the inlet as a velocity boundary, the initial velocity of the particles and drilling fluid, the nozzle outlet as a pressure boundary, and the outlet pressure value.
In an alternative embodiment, the specific process of analyzing the effect of the particle impact rock breaking simulation comprises the following steps: recording the fracture condition of the parallel bond, and representing the fracture degree of the rock mass by using the fracture condition of the parallel bond;
the stress distribution conditions of the pipeline and the nozzles in the x direction, the y direction and the z direction in the coupling process are obtained, and the abrasion degree of the pipeline is represented by the stress in the pipeline;
and obtaining the damage effect of the rock stratum under the current simulation parameters and the abrasion degree of the three-dimensional high-pressure spray pipe.
As an alternative implementation mode, the rock breaking parameters are changed, the specific process of carrying out the multiple simulation processes comprises changing drilling fluid speed, particle volume fraction, pipeline pressure or/and particle size, carrying out multiple computational fluid dynamics-digital elevation model coupling simulation calculation, obtaining simulation results of rock breaking degree, pipeline abrasion range or/and nozzle breakage range under different rock breaking parameters, and obtaining the response relation between the rock breaking parameters and the simulation results.
A particle impact rock breaking process simulation optimization system, comprising:
the rock stratum model building module is configured to acquire space data of a rock stratum, acquire structural plane information of the rock stratum and build a model for the rock stratum;
a macro mechanical parameter acquisition module configured to acquire macro mechanical parameters of a rock sample at a bottom of a target well;
the parameter calibration module is configured to fill particles into the established rock stratum model according to a particle grading curve, transfer force and moment between the particles through the parallel bonding model, and perform parameter calibration according to macroscopic mechanical parameters so that the strength index of the generated rock stratum model is consistent with that of the acquired rock sample;
the three-dimensional pipeline model building module is configured to build a three-dimensional pipeline model according to the rock breaking spray pipe used in the workplace;
an impact simulation module configured to simulate the flow of drilling fluid in a three-dimensional pipeline using computational fluid dynamics, and simulate the flow of impact particles in the drilling fluid and the process of jetting onto a rock mass using a digital elevation model;
a computing module configured to implement characterization of a degree of rock mass fragmentation based on fracture of parallel bond bonds in the formation model; representing the wear range of the pipeline based on the stress distribution condition in the spray pipe; the blocking condition in the spray pipe is represented based on the particle accumulation position in the spray pipe, and the damage effect of the particle impact rock breaking process on the rock mass under the current simulation parameters is obtained;
the scheme confirming module is configured to change rock breaking parameters, develop a plurality of simulation processes, calculate response relations between various rock breaking parameters and simulation results, and determine a particle impact rock breaking process optimization scheme suitable for actual operation conditions based on analysis results.
A computer readable storage medium having stored therein a plurality of instructions adapted to be loaded by a processor of a terminal device and to perform the steps in the method.
A terminal device comprising a processor and a computer readable storage medium, the processor configured to implement instructions; the computer readable storage medium is for storing a plurality of instructions adapted to be loaded by a processor and to perform the steps in the method.
Compared with the prior art, the invention has the beneficial effects that:
(1) The traditional rock stratum modeling can not accurately obtain structural plane information such as joints, cracks and the like in the rock stratum, the space data of the working rock stratum is obtained through a three-dimensional laser scanning technology and a digital photogrammetry technology, the real state of the rock stratum can be accurately obtained, and the high-precision modeling of the rock stratum is realized.
(2) In the particle jet impact rock breaking process, high-speed particles and high-speed fluid form a complex two-phase flow, the traditional particle impact rock breaking simulation is based on a finite element method, a finite difference method and the like, and a single simulation calculation method cannot simulate the interaction of drilling fluid and impact particles.
(3) The rock mass breaking degree is represented by adopting the parallel bond breaking condition, the abrasion ranges of the pipeline and the nozzle are represented by adopting the stress distribution condition, and the CFD-DEM coupling simulation calculation is carried out for a plurality of times based on the rock breaking parameters such as the drilling fluid speed, the particle volume fraction, the pipeline pressure, the particle size and the like, so that the corresponding relation between the rock breaking parameters and the rock breaking effect can be obtained, the rock breaking process is optimized, and theoretical guidance is provided for the site.
Drawings
The accompanying drawings, which are included to provide a further understanding of the invention and are incorporated in and constitute a part of this specification, illustrate embodiments of the invention and together with the description serve to explain the invention.
FIG. 1 is a flow chart of a CFD-DEM simulation process for breaking rock by high-speed particle impact;
fig. 2 is a schematic diagram of a calculation model of the present invention, wherein 1 is a jet pipe model, 2 is high-speed impact particles, 3 is a DEM rock mass model, and 4 is broken rock mass DEM particles.
Detailed Description
The invention will be further described with reference to the drawings and examples.
It should be noted that the following detailed description is illustrative and is intended to provide further explanation of the invention. Unless defined otherwise, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention belongs.
It is noted that the terminology used herein is for the purpose of describing particular embodiments only and is not intended to be limiting of exemplary embodiments according to the present invention. As used herein, the singular is also intended to include the plural unless the context clearly indicates otherwise, and furthermore, it is to be understood that the terms "comprises" and/or "comprising" when used in this specification are taken to specify the presence of stated features, steps, operations, devices, components, and/or combinations thereof.
As shown in fig. 1, in one or more embodiments, a CFD-DEM coupling-based particle impact rock breaking process simulation optimization simulation method is disclosed, which specifically includes the following steps:
step 1: high-precision modeling of a working rock stratum is realized based on a three-dimensional laser scanning technology and a digital photogrammetry technology, and the method specifically comprises the following steps:
step 1.1: the three-dimensional laser scanner and the digital photographic measuring instrument are transported to the bottom of a working well, and three-dimensional scanning is carried out on the working surface of a rock stratum at the bottom of the well to obtain structural surface information such as a joint and a crack of the rock stratum;
step 1.2: importing the acquired working rock stratum space data into solidworks software, and establishing a high-precision three-dimensional model of the rock stratum;
step 1.3: the built three-dimensional stratum model is imported into DEM discrete meta-software, and then particles are filled;
step 2: acquiring a rock sample at the bottom of a working well, carrying out an indoor rock mechanics experiment, and acquiring macroscopic mechanical parameters of the working rock stratum, wherein the method specifically comprises the following steps:
step 2.1: obtaining 7-9 groups of rock samples by using a drilling machine, and conveying the rock samples to a laboratory to carry out an indoor rock mechanics experiment;
step 2.2: screening the obtained rock sample in a laboratory to obtain a particle grading curve of the rock sample, preparing the rock sample into a standard sample, and performing a uniaxial compression experiment, a Brazil split experiment and a triaxial compression experiment to obtain the Poisson epsilon, the uniaxial compressive strength UCS, the elastic modulus E, the tensile strength ten and the friction angle of the rock sample
Macroscopic mechanical parameters such as cohesive force theta, friction coefficient upsilon and the like;
step 3: filling particles into a rock three-dimensional model in DEM software, and carrying out DEM discrete element parameter calibration, wherein the method specifically comprises the following steps:
step 3.1: filling particles into a rock three-dimensional model according to the particle grading curve of the rock sample obtained in the step 2, transmitting force and moment between the particles through a parallel bonding model, and carrying out DEM parameter calibration according to macroscopic mechanical physical mechanical parameters to ensure that the strength index of the generated rock model is consistent with the obtained rock sample;
step 3.2: poisson's ratio epsilon, uniaxial compressive strength UCS, young's modulus E, tensile strength ten and friction angle obtained according to step 2
Carrying out DEM parameter calibration on macroscopic mechanical parameters such as cohesive force theta, friction coefficient upsilon and the like to obtain parameters of a parallel bonding model, including elastic modulus E p Normal bond strength delta p Tangential bond strength τ p Normal stiffness k n Tangential stiffness k s Normal phase stiffness and tangential stiffness ratio k n / s The strength index of the DEM rock model generated according to the parallel bonding model is consistent with that of the scene;
step 4: establishing a spray pipe three-dimensional model in sol idworks, importing the model into DEM software, and filling wall units into the model, wherein the method specifically comprises the following steps of:
step 4.1: according to the spray pipe used in the field drilling operation, a three-dimensional model is built in the solidworks according to the same size ratio, and the length from the inlet to the outlet of the spray pipe is L R ；
Step 4.2: introducing the three-dimensional spray pipe model into DEM software, filling the wall units into the model, and measuring the distance L from the outlet of the spray pipe to the rock mass p Is L R /10；
Step 5: performing CFD-DEM coupling simulation, simulating the flow of drilling fluid in a three-dimensional pipeline by adopting the CFD, and simulating the flow of impact particles in the drilling fluid and the process of jetting the impact particles on a rock mass by adopting the DEM, wherein the method specifically comprises the following steps of:
step 5.1: setting parameters and boundary conditions of CFD-DEM coupling, wherein initial values comprise particle diameter d of impact particles and particle density ρ p Poisson ratio gamma of particles p Particle static friction system mu r Coefficient of rolling friction mu of particles s Rate of particle formation v s Coefficient of particle recovery mu j Volume fraction omega of particles p Coefficient of friction mu for pipe wall t The boundary condition setting procedure is as follows: setting the inlet as a speed boundary, and the initial speeds of the particles and the drilling fluid are v 1 The outlet of the spray pipe is set as a pressure boundary, and the outlet pressure is set as P 1 。
Step 5.2: generating a fluid field
And dividing proper fluid calculation grids at the same spatial position of the pipeline according to the spatial dimension of the established three-dimensional pipeline model, and simulating the flow of drilling fluid in the pipeline.
Step 5.3: generation of impact particles
In the DEM software, the particle parameters set according to the step 5.1 comprise the particle diameter d and the particle density ρ p Particle poisson ratio gamma p Particle static friction system mu r Coefficient of rolling friction mu of particles s Particle recovery coefficient mu j Particle generation rate v s Particle volume fraction omega p Etc. generating impact particle particles at the inlet of the pipeline according to the volume fraction, and generating positions of the particles in the x directionSet to 0<x<And 10d, setting a contact model between the impact particles and the impact particles, between the impact particles and the pipeline wall body and between the impact particles and the rock layer as a linear contact model for simulating the physical process of the impact particles in the rock breaking process.
Step 5.4: CFD-DEM coupling calculation
Fluid calculation domain time step is set to Δt 1 Setting the DEM calculation domain time step to be deltat 2 And DEM time step Deltat 2 <Δt 1 Setting the total simulation time as T, carrying out CFD-DEM simulation calculation, and when T>And at T, stopping calculation.
Step 6: the particle impact rock breaking simulation effect is analyzed, and the method is specifically as follows:
step 6.1: in the DEM software, the self-contained language of the DEM software is utilized to record the fracture condition of the parallel bond, and the fracture condition of the parallel bond is used for representing the fracture degree of the rock mass.
Step 6.2: in DEM software, a post-processing function is utilized to obtain stress distribution conditions of a pipeline and a nozzle in the x direction, the y direction and the z direction in the CFD-DEM coupling process, and the abrasion degree of the pipeline is represented by the stress in the pipeline;
step 6.3: and obtaining the damage effect of the rock stratum under the current simulation parameters by using post-processing software, and the abrasion degree of the three-dimensional high-pressure spray pipe.
Step 7: and (3) changing rock breaking parameters such as drilling fluid speed, particle volume fraction, pipeline pressure, particle size and the like, carrying out CFD-DEM coupling simulation calculation for a plurality of times, obtaining simulation results such as rock breaking degree, pipeline abrasion range, nozzle breakage range and the like under different rock breaking parameters, and obtaining a response relation between the rock breaking parameters and the simulation results.
Step 8: based on the analysis result, a particle impact rock breaking process optimization scheme suitable for actual operation conditions is obtained, and guidance is provided for an actual site.
It will be appreciated by those skilled in the art that embodiments of the present invention may be provided as a method, system, or computer program product. Accordingly, the present invention may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present invention may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The present invention is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems) and computer program products according to embodiments of the invention. It will be understood that each flow and/or block of the flowchart illustrations and/or block diagrams, and combinations of flows and/or blocks in the flowchart illustrations and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
The above description is only of the preferred embodiments of the present invention and is not intended to limit the present invention, but various modifications and variations can be made to the present invention by those skilled in the art. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present invention should be included in the protection scope of the present invention.
While the foregoing description of the embodiments of the present invention has been presented in conjunction with the drawings, it should be understood that it is not intended to limit the scope of the invention, but rather, it is intended to cover all modifications or variations within the scope of the invention as defined by the claims of the present invention.
Claims (10)
1. A particle impact rock breaking process simulation optimization method is characterized by comprising the following steps:
acquiring space data of a rock stratum, acquiring structural surface information of the rock stratum, and establishing a model for the rock stratum;
acquiring macroscopic mechanical parameters of a rock sample at the bottom of a target well;
filling particles into an established rock stratum model according to a particle grading curve, transmitting force and moment between the particles through a parallel bonding model, and carrying out parameter calibration according to macroscopic mechanical parameters to ensure that the strength index of the generated rock stratum model is consistent with that of an obtained rock sample;
according to a rock breaking spray pipe used in a working place, a three-dimensional pipeline model is established;
simulating the flow of drilling fluid in a three-dimensional pipeline by adopting computational fluid dynamics, and simulating the flow of impact particles in the drilling fluid and the process of spraying the impact particles on a rock mass by adopting a digital elevation model;
the rock mass breaking degree is characterized based on the breaking of parallel bonding bonds in the rock stratum model; representing the wear range of the pipeline based on the stress distribution condition in the spray pipe; the blocking condition in the spray pipe is represented based on the particle accumulation position in the spray pipe, and the damage effect of the particle impact rock breaking process on the rock mass under the current simulation parameters is obtained;
changing rock breaking parameters, carrying out multiple simulation processes, calculating response relations between various rock breaking parameters and simulation results, and determining a particle impact rock breaking process optimization scheme suitable for actual operation conditions based on analysis results.
2. The method of claim 1, wherein the macro-mechanical parameters include parameters selected from the group consisting of particle size distribution curve, uniaxial compressive strength, tensile strength, modulus of elasticity, poisson's ratio, angle of friction, cohesion and coefficient of friction.
3. The particle impact rock breaking process simulation optimization method according to claim 1, wherein the specific process of establishing a three-dimensional pipeline model comprises the steps of establishing a three-dimensional model according to the same size ratio according to a spray pipe used for on-site drilling operation, and setting the length from an inlet to an outlet of the spray pipe;
and filling the wall body units into the three-dimensional model, and setting the distance from the outlet of the matched spray pipe to the rock mass according to the length from the inlet to the outlet of the spray pipe.
4. The method for simulating and optimizing a particle impact rock breaking process according to claim 1, wherein the process of simulating the flow of drilling fluid in a three-dimensional pipeline by adopting computational fluid dynamics and simulating the flow of impact particles in the drilling fluid and spraying the impact particles on a rock mass by adopting a digital elevation model comprises the following steps:
setting parameters and boundary conditions for calculating the coupling of the fluid dynamics-digital elevation model;
dividing proper fluid calculation grids at the same spatial position of the pipeline according to the spatial size of the established three-dimensional pipeline model, and simulating the flow of drilling fluid in the pipeline;
generating impact particle particles at the inlet of the pipeline according to the volume fraction according to the set parameters of the coupling of the computational fluid dynamics-digital elevation model, setting the generation positions of the particles in the x direction, setting contact models between the impact particles and the impact particles, between the impact particles and the pipeline wall, and between the impact particles and the rock stratum to be linear contact models, and simulating the physical process of the impact particles in the rock breaking process;
setting a fluid calculation domain time step and a digital elevation model calculation domain time step, and calculating according to the set step until the total simulation time is met.
5. The method of claim 4, wherein the parameters include impact particle diameter, particle density, particle poisson's ratio, particle static friction system, particle rolling friction coefficient, particle generation rate, particle recovery coefficient, particle volume fraction, and coefficient of friction of the pipe wall;
boundary conditions include setting the inlet as a velocity boundary, the initial velocity of the particles and drilling fluid, the nozzle outlet as a pressure boundary, and the outlet pressure value.
6. The method for optimizing the simulation of a particle impact rock breaking process according to claim 1, wherein the specific process of analyzing the effect of the particle impact rock breaking process comprises the following steps: recording the fracture condition of the parallel bond, and representing the fracture degree of the rock mass by using the fracture condition of the parallel bond;
the stress distribution conditions of the pipeline and the nozzles in the x direction, the y direction and the z direction in the coupling process are obtained, and the abrasion degree of the pipeline is represented by the stress in the pipeline;
and obtaining the damage effect of the rock stratum under the current simulation parameters and the abrasion degree of the three-dimensional high-pressure spray pipe.
7. The method for optimizing the simulation of the particle impact rock breaking process according to claim 1, wherein the specific process of changing the rock breaking parameters and carrying out the multiple simulation processes comprises changing drilling fluid speed, particle volume fraction, pipeline pressure or/and particle size, carrying out multiple computational fluid dynamics-digital elevation model coupling simulation calculation, obtaining simulation results of rock breaking degree, pipeline abrasion range or/and nozzle breakage range under different rock breaking parameters, and obtaining the response relation between the rock breaking parameters and the simulation results.
8. A particle impact rock breaking process simulation optimization system, comprising:
the rock stratum model building module is configured to acquire space data of a rock stratum, acquire structural plane information of the rock stratum and build a model for the rock stratum;
a macro mechanical parameter acquisition module configured to acquire macro mechanical parameters of a rock sample at a bottom of a target well;
the parameter calibration module is configured to fill particles into the established rock stratum model according to a particle grading curve, transfer force and moment between the particles through the parallel bonding model, and perform parameter calibration according to macroscopic mechanical parameters so that the strength index of the generated rock stratum model is consistent with that of the acquired rock sample;
the three-dimensional pipeline model building module is configured to build a three-dimensional pipeline model according to the rock breaking spray pipe used in the workplace;
an impact simulation module configured to simulate the flow of drilling fluid in a three-dimensional pipeline using computational fluid dynamics, and simulate the flow of impact particles in the drilling fluid and the process of jetting onto a rock mass using a digital elevation model;
a computing module configured to implement characterization of a degree of rock mass fragmentation based on fracture of parallel bond bonds in the formation model; representing the wear range of the pipeline based on the stress distribution condition in the spray pipe; the blocking condition in the spray pipe is represented based on the particle accumulation position in the spray pipe, and the damage effect of the particle impact rock breaking process on the rock mass under the current simulation parameters is obtained;
the scheme confirming module is configured to change rock breaking parameters, develop a plurality of simulation processes, calculate response relations between various rock breaking parameters and simulation results, and determine a particle impact rock breaking process optimization scheme suitable for actual operation conditions based on analysis results.
9. A computer readable storage medium, characterized in that a plurality of instructions are stored, said instructions being adapted to be loaded by a processor of a terminal device and to perform the steps of the method of any of claims 1-7.
10. A terminal device, comprising a processor and a computer readable storage medium, the processor configured to implement instructions; a computer readable storage medium for storing a plurality of instructions adapted to be loaded by a processor and to perform the steps of the method of any of claims 1-7.
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211548562.0A CN116306332A (en) | 2022-12-05 | 2022-12-05 | Particle impact rock breaking process simulation optimization method and system |
| PCT/CN2023/095865 WO2024119726A1 (en) | 2022-12-05 | 2023-05-23 | Particle impacting rock breaking process simulation and optimization method and system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211548562.0A CN116306332A (en) | 2022-12-05 | 2022-12-05 | Particle impact rock breaking process simulation optimization method and system |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN116306332A true CN116306332A (en) | 2023-06-23 |
Family
ID=86822787
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202211548562.0A Pending CN116306332A (en) | 2022-12-05 | 2022-12-05 | Particle impact rock breaking process simulation optimization method and system |
Country Status (2)
| Country | Link |
|---|---|
| CN (1) | CN116306332A (en) |
| WO (1) | WO2024119726A1 (en) |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109977505B (en) * | 2019-03-13 | 2024-02-09 | 南京大学 | Discrete element pore system quick searching method based on GPU matrix calculation |
| WO2021022042A1 (en) * | 2019-07-30 | 2021-02-04 | William Marsh Rice University | System and method for design of rapid excavating and wear-resistant drill bits |
| CN110473597B (en) * | 2019-09-05 | 2022-05-03 | 中国石油大学(北京) | Conglomerate mechanical property evaluation and analysis method and system |
| CN111024588B (en) * | 2019-12-31 | 2021-04-13 | 山东大学 | DEM contact model construction method for reflecting weakening of seepage to rock and soil mass strength |
-
2022
- 2022-12-05 CN CN202211548562.0A patent/CN116306332A/en active Pending
-
2023
- 2023-05-23 WO PCT/CN2023/095865 patent/WO2024119726A1/en unknown
Also Published As
| Publication number | Publication date |
|---|---|
| WO2024119726A1 (en) | 2024-06-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109063257B (en) | A kind of coal and rock subregion water filling seepage flow-damage-stress coupling method for numerical simulation | |
| Manchanda et al. | Integrating reservoir geomechanics with multiple fracture propagation and proppant placement | |
| Kong et al. | Fluid-driven micro-cracking behaviour of crystalline rock using a coupled hydro-grain-based discrete element method | |
| Jaime | Numerical modeling of rock cutting and its associated fragmentation process using the finite element method | |
| CN105787206A (en) | Design method of anchor bolt support for highway tunnel | |
| CA2919553C (en) | Injection treatment simulation using condensation | |
| CN103913772B (en) | Micro-seismic event the Forward Modeling based on reservoir geology mechanics parameter | |
| CN104963657A (en) | Oil reservoir numerical simulation method and apparatus based on high water-cut period monitoring data constraints | |
| Lin et al. | A finite-discrete element based appoach for modelling the hydraulic fracturing of rocks with irregular inclusions | |
| CN105089595A (en) | Oil reservoir numerical simulation method and device under the action of horizontal fracturing fracture diversion | |
| Pettitt et al. | Fracture network engineering for hydraulic fracturing | |
| Sanchez et al. | Numerical simulation of three-dimensional fracture interaction | |
| Shahverdi et al. | Numerical simulation of abrasive water jet cutting process using the SPH and ALE methods | |
| Zhang et al. | Visualization and investigation of the erosion process for natural gas hydrate using water jet through experiments and simulation | |
| Settgast et al. | Simulation of hydraulic fracture networks in three dimensions utilizing massively parallel computing resources | |
| Shen | Modeling rock fracturing processes with FRACOD | |
| Chen et al. | Hydraulic fracturing analysis in fluid‐saturated porous medium | |
| Kong et al. | Experimental and numerical investigations on crack propagation characteristics of rock-like specimens with preexisting flaws subjected to combined actions of internal hydraulic pressure and shear force | |
| Basson et al. | Simulating hydraulic fracturing preconditioning in mines with the material point method | |
| CN116306332A (en) | Particle impact rock breaking process simulation optimization method and system | |
| CN115292975A (en) | Three-dimensional continuous-discontinuous coupling numerical simulation method for shield type TBM excavation card machine | |
| Rasouli et al. | A true triaxial stress cell to simulate deep downhole drilling conditions | |
| Lisjak et al. | Development of a fully-coupled, hydro-mechanical model for finite-discrete element simulations of fluid-driven fracturing | |
| Sidorová et al. | Plastic Flow Modeling in Rock Fracture. | |
| Qin et al. | Solid-gas coupling law during methane seepage from a coal mass in the advanced pressure relief area of a mining seam |
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 |