Mining planning method and system, storage medium and equipment

Abstract

The invention provides a mining planning method. The mining planning method comprises the following steps: acquiring physical characteristics and initial parameters of an ore particle swarm; carrying out numerical simulation on the ore particle swarms according to the physical characteristics and the initial parameters, calculating resistance coefficients of the ore particle swarms at different flow velocities, and generating a mathematical relationship between the resistance coefficients and flow velocity changes; calculating the settling velocity of the ore particle swarm according to the mathematical relationship; calculating the minimum lifting speed of the ore particles according to the settling speed; the minimum lifting speed is used for guiding mining planning. According to the method, the resistance coefficient and the settling velocity of the particle swarm are accurately calculated through numerical simulation, then the minimum lifting speed is determined, it is ensured that ore particles can be smoothly lifted, and an important reference basis is provided for design and optimization of a deep-sea mining hydraulic lifting system. The invention also provides a mining planning system, a computer readable storage medium and electronic equipment, which have the above beneficial effects.

Description

Technical Field

[0001] The present application relates to the field of deep-sea mining, and in particular to a mining planning method, system, storage medium and device. Background Art

[0002] Deep-sea mining involves the process of transporting seabed ore to the surface through a hydraulic lifting system. In order to ensure that the ore particles can be lifted smoothly, their minimum lifting speed must be determined. Traditional methods usually focus on the drag coefficient and settling velocity of individual particles, but this method may not be accurate enough when dealing with a group of particles. Therefore, how to more accurately determine the lifting speed is a technical problem that technicians in this field need to solve urgently. Summary of the invention

[0003] The purpose of this application is to provide a mining planning method, system, computer-readable storage medium and electronic device, which calculates the resistance coefficient of a particle group through numerical simulation, and combines the settling velocity of the ore particle group to finally calculate the minimum lifting velocity, thereby providing a reliable basis for the setting of hydraulic lifting parameters in deep-sea mining.

[0004] In order to solve the above technical problems, this application provides a mining planning method, and the specific technical solution is as follows:

[0005] Obtain the physical characteristics and initial parameters of the ore particle group;

[0006] Performing numerical simulation on the ore particle group according to the physical characteristics and the initial parameters, calculating the resistance coefficient of the ore particle group at different flow rates, and generating a mathematical relationship between the resistance coefficient and the flow rate change;

[0007] Calculating the settling velocity of the ore particle group according to the mathematical relationship;

[0008] The minimum lifting speed of the ore particles is calculated according to the settling speed; the minimum lifting speed is used to guide mining planning.

[0009] Optionally, performing numerical simulation on the ore particle group according to the physical properties and the initial parameters includes:

[0010] Constructing an initial geometric model of the ore particle group according to the physical characteristics and the initial parameters; the initial geometric model includes ore particles of different sizes and shapes;

[0011] A coupling model of the fluid domain and the solid domain is established based on the initial geometric model.

[0012] Optionally, establishing a coupling model of a fluid domain and a solid domain based on the initial geometric model includes:

[0013] The Navier-Stokes equations are used to describe the movement of ore particles in the fluid domain;

[0014] Using Newton's second law to describe the solid domain motion of the ore particles;

[0015] A coupling model is established according to the fluid domain motion and the solid domain motion.

[0016] Optionally, establishing a coupling model according to the fluid domain motion and the solid domain motion includes:

[0017] Define Lagrange points and Eulerian grids;

[0018] The Lagrangian point and the Euler grid are associated by using an interpolation function, wherein the interpolation function is used to convert the force on the solid boundary into the force on the fluid grid point, and to transfer the reaction force of the fluid on the solid boundary from the fluid grid point to the Lagrangian point.

[0019] Optionally, the associating the Lagrange point and the Euler grid by using an interpolation function comprises:

[0020] Step A: Initialize the physical parameters and initial positions of the fluid domain and solid particles, and construct the Euler grid of the fluid domain and the Lagrange points of the solid boundary;

[0021] Step B: Use the Navier-Stokes equations to solve the velocity and pressure fields in the fluid domain;

[0022] Step C: Calculate the force on the solid boundary, apply the force to the fluid domain, and update the position and velocity of the Lagrangian point;

[0023] Step D: using the interpolation function to transfer the force of the solid boundary to the fluid grid point, and calculating the reaction force of the fluid on the solid boundary;

[0024] Step E: Update the velocity field in the fluid domain and the position and velocity of the solid particles through the flow field in the fluid domain and the motion equation of the solid domain;

[0025] Iteration step A to step E, perform iterative calculation of the time step, and stop the iteration when the iteration time meets the predetermined simulation time or meets the convergence condition.

[0026] Optionally, calculating the settling velocity of the ore particle group according to the mathematical relationship includes:

[0027] Generate a force model of the ore particles; the force model includes the interaction of gravity, buoyancy and fluid resistance;

[0028] Calculating the settling velocity of the ore particles according to the gravity, the buoyancy and the fluid resistance;

[0029] The settling velocity of the ore particle group is determined according to the overall resistance coefficient and the average diameter of the ore particle group.

[0030] Optionally, calculating the minimum lifting speed of ore particles according to the settling speed includes:

[0031] The minimum lifting velocity of the ore particles is calculated based on the settling velocity and the design parameters of the hydraulic lifting system.

[0032] The present application also provides a mining planning system, comprising:

[0033] A data acquisition module, used to obtain the physical characteristics and initial parameters of the ore particle group;

[0034] A numerical simulation module, used for performing numerical simulation on the ore particle group according to the physical characteristics and the initial parameters, calculating the resistance coefficient of the ore particle group at different flow rates, and generating a mathematical relationship between the resistance coefficient and the flow rate change;

[0035] A settling velocity calculation module, used to calculate the settling velocity of the ore particle group according to the mathematical relationship;

[0036] A mining planning module is used to calculate the minimum lifting speed of ore particles according to the settling speed; the minimum lifting speed is used to guide mining planning.

[0037] The present application also provides a computer-readable storage medium having a computer program stored thereon, and when the computer program is executed by a processor, the steps of the mining planning method described above are implemented.

[0038] The present application also provides an electronic device, including a memory and a processor, wherein a computer program is stored in the memory, and when the processor calls the computer program in the memory, the steps of the mining planning method described above are implemented.

[0039] The present application provides a mining planning method, comprising: obtaining physical properties and initial parameters of an ore particle group; performing numerical simulation on the ore particle group according to the physical properties and the initial parameters, calculating the resistance coefficient of the ore particle group at different flow rates, and generating a mathematical relationship between the resistance coefficient and the flow rate change; calculating the settling velocity of the ore particle group according to the mathematical relationship; calculating the minimum lifting velocity of the ore particles according to the settling velocity; the minimum lifting velocity is used to guide mining planning.

[0040] This application uses numerical simulation to accurately calculate the resistance coefficient and settling velocity of the particle group, and then determines the minimum lifting speed to ensure that the ore particles can be lifted smoothly, providing an important reference for the design and optimization of deep-sea mining hydraulic lifting systems.

[0041] The present application also provides a mining planning system, a computer-readable storage medium and an electronic device, which have the above-mentioned beneficial effects and are not described in detail here. BRIEF DESCRIPTION OF THE DRAWINGS

[0042] In order to more clearly illustrate the embodiments of the present application or the technical solutions in the prior art, the drawings required for use in the embodiments or the description of the prior art will be briefly introduced below. Obviously, the drawings described below are only embodiments of the present application. For ordinary technicians in this field, other drawings can be obtained based on the provided drawings without paying any creative work.

[0043] Figure 1 A flowchart of a mining planning method provided in an embodiment of the present application;

[0044] Figure 2 A flowchart of numerical simulation using the immersed boundary method provided in an embodiment of the present application;

[0045] Figure 3 A schematic diagram of an initial collection model of a group of ore particles provided in an embodiment of the present application;

[0046] Figure 4 A flowchart of another mining planning method provided in an embodiment of the present application;

[0047] Figure 5 A schematic diagram of the mining planning system structure provided in an embodiment of the present application;

[0048] Figure 6 A structural diagram of an electronic device provided in an embodiment of the present application. DETAILED DESCRIPTION

[0049] In order to make the purpose, technical solution and advantages of the embodiments of the present application clearer, the technical solution in the embodiments of the present application will be clearly and completely described below in conjunction with the drawings in the embodiments of the present application. Obviously, the described embodiments are part of the embodiments of the present application, not all of the embodiments. Based on the embodiments in the present application, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of this application.

[0050] The object information involved in this application, including but not limited to the object device information, object personal information, etc., and data, including but not limited to data used for analysis, stored data, displayed data, etc., are all information and data authorized by the object or fully authorized by all parties, and the collection, use and processing of relevant data must comply with the laws, regulations and standards of relevant countries and regions.

[0051] See also Figure 1 , Figure 1 A flowchart of a mining planning method provided in an embodiment of the present application, the method comprising:

[0052] S101: obtaining physical properties and initial parameters of the ore particle group;

[0053] S102: performing numerical simulation on the ore particle group according to the physical characteristics and the initial parameters, calculating the resistance coefficient of the ore particle group at different flow rates, and generating a mathematical relationship between the resistance coefficient and the flow rate change;

[0054] S103: Calculating the settling velocity of the ore particle group according to the mathematical relationship;

[0055] S104: Calculate the minimum lifting speed of ore particles according to the settling speed, and the minimum lifting speed is used to guide mining planning.

[0056] First, the physical properties and initial parameters of the ore particle group are determined. The physical properties of the ore particle group include the density, shape, size and initial distribution of the particles. These parameters are used to construct the initial geometric model of the particle group.

[0057] Direct numerical simulation of particle swarms is performed using the immersed boundary method. The specific process is as follows:

[0058] The first step is to build the initial geometric model of the particle group:

[0059] The model should include a variety of ore particles of different sizes and shapes to accurately reflect the actual situation.

[0060] Step 2: Establish a coupling model between the fluid domain and the solid domain:

[0061] This step requires coupling calculation of the fluid domain and the solid domain. In the calculation process based on the immersed boundary method, the coupling calculation of the fluid domain and the solid domain is a key part. This process involves transferring the force of the solid particle boundary to the fluid domain and affecting the movement of the solid particles through the reaction force of the fluid domain. The following is a detailed description and formula.

[0062] The motion of the fluid domain is described by the Navier-Stokes equations, which include the momentum equation and the continuity equation.

[0063] The momentum equation is expressed as follows:

[0064]

[0065] Where ρ represents the particle density (kg / m 3), d is the particle diameter (m), φ is the particle shape factor. u is the fluid velocity vector (m / s), t is the time (s), ρ is the fluid density (kg / m 3 ), p is the pressure (Pa), ν is the fluid dynamic viscosity (m 2 / s), f is the body force (N / m 3 ).

[0066] The continuity equation (conservation of mass) is expressed as follows:

[0067] ▽·u=0.

[0068] In the solid domain equation, the motion of solid particles is described by Newton's second law, which mainly considers the force balance expression:

[0069]

[0070] Among them, m s is the particle mass (kg), X s is the position vector of the solid particle, F s is the force exerted by the fluid on the solid (N), F b is the interaction force between solids (N).

[0071] The immersed boundary method realizes the coupling of fluid domain and solid domain through the interaction of Lagrange points and Euler grid. The specific process is as follows:

[0072] First, define Lagrange points and Euler grids: Lagrange points are used to represent solid boundaries, and Euler grids are used to describe fluid domains. The two are connected by interpolation functions. s represents the coordinates of the Lagrange point on the solid boundary, X i Represents the coordinates of the fluid grid points.

[0073] After that, the interpolation function is defined. In order to realize the interaction between Lagrange points and Euler grids, the interpolation function is defined to represent

[0074]

[0075] Through the interpolation function, the force F on the solid boundary is s Applied to the fluid domain, converted into a force f on the fluid grid point i express:

[0076]

[0077] Similarly, the reaction force of the fluid on the solid boundary is transferred from the fluid grid point to the Lagrange point through the interpolation function:

[0078]

[0079] Thereafter, coupling calculation is performed to associate the Lagrange point and the Euler grid using an interpolation function, which may include the following steps:

[0080] Step A: Initialize the physical parameters and initial positions of the fluid domain and solid particles, and construct the Euler grid of the fluid domain and the Lagrange points of the solid boundary;

[0081] Step B: Use the Navier-Stokes equations to solve the velocity and pressure fields in the fluid domain;

[0082] Step C: Calculate the force on the solid boundary, apply the force to the fluid domain, and update the position and velocity of the Lagrangian point;

[0083] Step D: using the interpolation function to transfer the force of the solid boundary to the fluid grid point, and calculating the reaction force of the fluid on the solid boundary;

[0084] Step E: Update the velocity field in the fluid domain and the position and velocity of the solid particles through the flow field in the fluid domain and the motion equation of the solid domain;

[0085] Iteration step A to step E, perform iterative calculation of the time step, and stop the iteration when the iteration time meets the predetermined simulation time or meets the convergence condition.

[0086] In step A, the physical parameters and initial positions of the fluid domain and solid particles are initialized. The Eulerian mesh of the fluid domain and the Lagrangian points of the solid boundary are constructed.

[0087] In step B, the Navier-Stokes equations are used to solve the velocity and pressure fields in the fluid domain.

[0088] In step C, the force F acting on the solid boundary is calculated. s , and applies it to the fluid domain. After that, the position and velocity of the Lagrangian point are updated.

[0089] In step D, the interpolation function is used to transfer the force of the solid boundary to the fluid grid points, and the reaction force of the fluid on the solid boundary is calculated.

[0090] See also Figure 2 , Figure 2 A flow chart of numerical simulation using the immersed boundary method provided in an embodiment of the present application. Figure 2 The complete process of numerical simulation using the immersed boundary method is demonstrated, including the following processes:

[0091] Firstly, an initial geometric model is constructed. The initial geometric model of the ore particle group includes a plurality of ore particles of different sizes and shapes.

[0092] Secondly, a coupling model of the fluid domain and the solid domain is established: the Navier-Stokes equations are used in the fluid domain to describe the fluid motion. Newton's second law is used in the solid domain to describe the motion of particles. The coupling of the fluid domain and the solid domain is achieved through the interaction of Lagrange points and Euler grids.

[0093] Then, numerical calculations are applied to calculate the flow field and the forces. The velocity field and pressure field of the fluid domain are solved by numerical simulation. The forces on the solid boundary are calculated and applied to the fluid domain.

[0094] When you update the fluid and solid domains, you update the velocity field of the fluid domain and the positions and velocities of the solid particles.

[0095] The calculation is iterated until convergence or the final time step is reached.

[0096] Repeat the above steps to perform iterative calculation of the time step until the iteration time meets the predetermined final time step or meets the convergence condition, that is, perform iterative calculation of the time step until the predetermined simulation time or convergence condition is reached.

[0097] An exemplary calculation process is as follows:

[0098] Assume that the initial conditions of the solid particles are:

[0099] Particle density ρ s =2500kg / m 3 ;

[0100] Particle diameter d = 0.01m;

[0101] Fluid density ρ f =1000kg / m 3 ;

[0102] Dynamic viscosity ν=1×10 -6 m 2 / s;

[0103] Initialize the parameters of the fluid domain and solid particles, and construct the Eulerian mesh and Lagrange points.

[0104] The initial fluid velocity field is calculated as:

[0105] u(t=0)=0;

[0106] Calculate the initial solid boundary force as:

[0107]

[0108] Apply the force and update the fluid velocity field as:

[0109]

[0110] The reaction force of the fluid on the solid boundary is calculated as:

[0111]

[0112] By repeatedly iterating the above steps, the coupled calculation of the fluid domain and the solid domain is realized until the convergence condition or the predetermined simulation time is reached.

[0113] u(t+Δt)=u(t)+Δt[-(u·▽)u+ν▽ 2 u+f];

[0114] X s (t+Δt)=X s (t)+V s Δt;

[0115] After numerical simulation of the ore particle group, the drag coefficient of the particle group is calculated. There is no limitation on how to calculate the drag coefficient. A feasible drag coefficient calculation formula is as follows:

[0116] Cd=Fd / (0.5*ρ*u 2 *A);

[0117] Where Cd represents the drag coefficient, Fd represents the drag force on the particle group (N), and ρ f Indicates fluid density (kg / m 3 ), u represents the fluid velocity (m / s), and A represents the frontal area of ​​the particle group (m 2 ), Vs represents the particle velocity.

[0118] By calculating the resistance coefficient at different flow rates through numerical simulation, the mathematical relationship between the resistance coefficient and the flow rate change can be determined. For example, a curve chart showing the resistance coefficient changing with the flow rate can be drawn to facilitate the subsequent calculation of the settling velocity of the particle group.

[0119] Calculating the sedimentation velocity of a particle group is one of the key steps in determining the minimum lifting velocity. The sedimentation velocity reflects the movement speed of a particle group after the force balance in a static or flowing fluid. The following describes in detail the process of calculating the sedimentation velocity and lists the relevant equations and physical models, which can specifically include the following steps:

[0120] The first step is to generate a force model of the ore particles; the force model includes the interaction of gravity, buoyancy and fluid resistance;

[0121] Step 2: calculating the settling velocity of the ore particles according to the gravity, the buoyancy and the fluid resistance;

[0122] The third step is to determine the settling velocity of the ore particle group according to the overall resistance coefficient and the average diameter of the ore particle group.

[0123] In a fluid, a group of particles is subject to gravity, buoyancy, and fluid resistance. For a single particle, the force balance equation can be expressed as

[0124] F g -F b -F d =0;

[0125] Among them, F g is gravity, F b is the buoyancy, F d is the fluid resistance.

[0126] For a group of particles, the force balance can be expanded to the sum of all particles:

[0127] Σ(F g -F b -F d )=0;

[0128] For the gravity and buoyancy in the above formula, the gravity F of a single particle g and buoyancy F b They are:

[0129]

[0130] Among them, m s is the particle mass, m f is the mass displaced by the fluid, d is the particle diameter, ρ s is the particle density, ρ f is the fluid density and g is the acceleration due to gravity.

[0131] Fluid resistance F d It is usually related to the velocity of the fluid and the shape of the particles, which can be expressed by the drag coefficient C d To express:

[0132]

[0133] For a single spherical particle, the frontal area A is:

[0134]

[0135] Drag coefficient C d It is obtained through numerical simulation and represents the resistance coefficient of particle groups at different flow rates.

[0136] The settling velocity V of the particle group tIt is the speed under the condition of force balance. By balancing gravity, buoyancy and resistance, we can get the calculation formula of sedimentation speed:

[0137] F g -F b =F d ;

[0138] Substituting into the above formula, we can get:

[0139]

[0140] After sorting, we get:

[0141]

[0142] The final settling velocity formula is:

[0143]

[0144] For a particle group, the sedimentation velocity needs to take into account the interaction between particles and the fluid flow characteristics. The sedimentation velocity of the particle group can be calculated by the overall resistance coefficient C of the particle group. d and average diameter To calculate:

[0145]

[0146] In the above formula, is the average diameter of a population of particles, obtained through statistical analysis.

[0147] Through the numerical simulation of the immersed boundary method, the resistance coefficient C of the particle group at different flow rates is obtained. d , and combined with the physical characteristics of the particle group, calculate the sedimentation velocity. The specific steps are as follows:

[0148] The first step is to determine the initial parameters of the particle group: including the density of the particles ρ s , average diameter and fluid density ρ f .

[0149] Step 2: Numerical simulation to obtain the resistance coefficient: Through the numerical simulation of the immersed boundary method, the resistance coefficient C of the particle group at different flow rates is obtained. d .

[0150] Step 3: Calculate the sedimentation velocity: Use the above sedimentation velocity formula combined with the C obtained by numerical simulation d , the settling velocity of the particle group can be calculated.

[0151] In order to better understand the above calculation process, assume that the parameters of a certain ore particle group are as follows:

[0152] ρs =2500kg / m 3 ;

[0153] ρ f =1000kg / m 3 ;

[0154]

[0155] The numerical simulation results show that the C d =0.05;

[0156] Substitute into the settling velocity formula:

[0157]

[0158] The calculation results are:

[0159] V t ≈0.9m / s;

[0160] Determine the sedimentation behavior of particle groups under different conditions to provide basic data for calculating the minimum lifting velocity.

[0161] Calculate the minimum lifting speed. First, analyze the basis for the minimum lifting speed being three times the settling speed.

[0162] In the hydraulic lifting process of deep-sea mining, it is necessary to determine the minimum lifting velocity (Minimum Lifting Velocity, V min ) is a key parameter to ensure that ore particles can be smoothly lifted from the deep sea bottom to the sea surface. In actual engineering, the minimum lifting speed is usually set to the particle group settling velocity (Settling Velocity, V t This empirical formula is derived from a combination of factors, including fluid dynamics, particle group interactions, and safety margins for practical operations.

[0163] In the lifting system, the flow state of the fluid has a significant impact on the movement of particles. The fluid flow in the lifting pipeline is usually in a turbulent state, that is, the velocity and pressure of the fluid fluctuate greatly, resulting in unstable movement of particles. Setting the minimum lifting speed to three times the settling speed helps to ensure that all particles can be stably lifted under turbulent conditions without settling due to instantaneous velocity fluctuations.

[0164] The particles in the particle group are not only affected by gravity and fluid resistance, but also collide and agglomerate with each other. These interactions will increase the complexity of particle lifting. By increasing the minimum lifting speed, the adverse effects of particle interactions on the overall lifting effect can be reduced, ensuring the overall rise of the particle group.

[0165] In addition, in engineering practice, a certain safety margin is usually considered to cope with unforeseen factors and errors in actual operation. Setting the minimum lifting speed to three times the settling speed provides a large safety margin to ensure the reliability and stability of the system under different operating conditions.

[0166] A number of practical operations and experimental studies have shown that setting the minimum lifting speed to three times the settling speed can effectively improve the efficiency and reliability of the lifting system. This empirical formula has been verified by long-term engineering practice and has become an important reference for design and operation.

[0167] See also Figure 2 , Figure 2 Schematic diagram of the initial collection model of the ore particle group provided in the embodiment of the present application, Figure 1 A schematic diagram of the initial geometry of a cluster of ore particles is shown, illustrating the distribution of particles of different sizes and shapes in the fluid. Figure 2 It can be seen that the ore particle group is subjected to mechanical action in the fluid and its sedimentation and lifting process.

[0168] Attached Figure 3 A schematic diagram of the initial geometry of a swarm of ore particles is shown in Figure 1. The figure contains particles of different sizes and shapes, and the distribution of these particles in the fluid is as follows:

[0169] For particle A: it is spherical, with a diameter d A , shape factor

[0170] For particle B: its ellipsoidal major axis a B Minor axis, b B Form Factor,

[0171] For particle C: it is a polyhedron with a characteristic side length of l C , shape factor

[0172] Figure 2 In the simulation, particles are randomly distributed in the fluid. Considering the complexity of the particle group in the actual situation, a random distribution method is adopted. A certain distance is maintained between the particles to avoid collision and interaction in the initial state.

[0173] Figure 2 The figure shows a three-dimensional coordinate system, where the X, Y, and Z axes represent the horizontal and vertical directions, respectively. The initial position of each particle is determined by the coordinates of its center point (x i ,y i ,z i )Sure.

[0174] The fluid boundary condition is assumed to be a no-slip boundary condition, that is, the velocity of the fluid on the particle surface is consistent with the particle velocity. The initial flow direction and velocity distribution are set in the fluid domain to simulate the movement of the particle group in the fluid.

[0175] The construction process of the initial geometric model is as follows:

[0176] First, use a mathematical model to generate particles of different sizes and shapes. For spherical particles, such as particle A, a standard sphere model is used. For ellipsoidal particles, such as particle B, an ellipsoid model defined by the major and minor axes is used. For polyhedral particles, such as particle C, a polyhedral model defined by the characteristic side length and vertex coordinates is used.

[0177] The shape factor is used to describe the effect of the particle geometry on its motion behavior. The Monte Carlo method is used to generate the initial distribution of particles. By setting the minimum spacing parameter, the initial distance between particles is ensured to be far enough to avoid collisions in the initial state.

[0178] When defining the fluid domain, the fluid domain uses a rectangular area, sets boundary conditions and initial velocity distribution, and imports the initial geometric model and fluid domain parameters through numerical simulation software.

[0179] Figure 3 The initial state of the ore particle group is shown, which helps to understand the force and movement behavior of the particle group in the fluid. By describing the shape, size and distribution of the particles in detail, numerical simulation and analysis can be better carried out, and the calculation accuracy and reliability can be improved.

[0180] See also Figure 4 , Figure 4 A flowchart of another mining planning method provided in an embodiment of the present application is shown in FIG. Figure 4 As shown, it may include the following steps:

[0181] Determine the physical characteristics and initial parameters of the ore particle population. Determine the density, shape, size and initial distribution of the ore particles, as well as the initial conditions and flow parameters of the fluid.

[0182] Build the initial geometry. Build the geometry of the ore particle cluster, including particles of different sizes and shapes, and place them in the fluid domain.

[0183] Establish a coupled model between the fluid domain and the solid domain. Use the Navier-Stokes equations to describe fluid flow in the fluid domain. Use Newton's second law to describe particle motion in the solid domain. The coupling of fluid and solid is achieved through the interaction between Lagrange points and Euler grids.

[0184] Numerical calculation of flow fields and forces. Solve the velocity and pressure fields in the fluid domain through numerical simulation. Calculate the forces on the solid boundary and apply them to the fluid domain to update the position and velocity of the particles.

[0185] Update fluid and solid domains. Update the velocity field of the fluid domain and the positions and velocities of the solid particles based on the numerical calculation results.

[0186] Iterate until convergence or the final time step is reached. Repeat the above steps and iterate through the time step until the predetermined simulation time or convergence condition is reached.

[0187] Calculate the drag coefficient of the particle group. According to the flow velocity and particle force obtained by numerical simulation, calculate the drag coefficient of the particle group (the ratio of the force on the particle group in the fluid to the fluid flow velocity).

[0188] Calculate the settling velocity of a particle group. Calculate the settling velocity of a particle group (the velocity of a particle group after force balance in a static or flowing fluid) based on the drag coefficient of the particle group and the fluid properties.

[0189] Calculate the minimum lifting speed. Calculate the minimum lifting speed based on the settling speed of the particle group and the design parameters of the hydraulic lifting system. It is usually set to three times the settling speed to ensure that the particle group can be lifted smoothly in the complex deep-sea mining environment.

[0190] Verify and optimize the minimum lifting speed. Verify the accuracy of the calculation results through actual operation and experiments. Based on the verification results, further optimize the design parameters of the lifting system to ensure the reliability and stability of the system under different operating conditions.

[0191] Through these steps, the dynamic characteristics of ore particle groups during hydraulic lifting can be accurately simulated, so that the minimum lifting speed can be calculated. The minimum lifting speed can be used to guide mining planning and provide a reliable basis for setting hydraulic lifting parameters in deep-sea mining.

[0192] See also Figure 5 , Figure 5 This is a schematic diagram of the structure of a mining planning system provided in an embodiment of the present application. The present application also provides a mining planning system, including:

[0193] A data acquisition module, used to obtain the physical characteristics and initial parameters of the ore particle group;

[0194] A numerical simulation module, used for performing numerical simulation on the ore particle group according to the physical characteristics and the initial parameters, calculating the resistance coefficient of the ore particle group at different flow rates, and generating a mathematical relationship between the resistance coefficient and the flow rate change;

[0195] A settling velocity calculation module, used to calculate the settling velocity of the ore particle group according to the mathematical relationship;

[0196] A mining planning module is used to calculate the minimum lifting speed of ore particles according to the settling speed; the minimum lifting speed is used to guide mining planning.

[0197] Based on the above embodiment, in a feasible embodiment, the numerical simulation module includes:

[0198] A first construction submodule is used to construct an initial geometric model of the ore particle group according to the physical characteristics and the initial parameters; the initial geometric model includes ore particles of different sizes and shapes;

[0199] The second construction submodule is used to establish a coupling model of the fluid domain and the solid domain based on the initial geometric model.

[0200] Based on the above embodiment, in a feasible embodiment, the second construction submodule includes:

[0201] A first description unit is used to describe the fluid domain movement of ore particles using the Navier-Stokes equations;

[0202] A first description unit is used to describe the solid domain motion of the ore particles using Newton's second law;

[0203] A construction unit is used to establish a coupling model according to the motion of the fluid domain and the motion of the solid domain.

[0204] Based on the above embodiment, in a feasible embodiment, the construction unit includes:

[0205] Define subcells to define Lagrange points and Eulerian grids;

[0206] An association subunit is used to associate the Lagrangian point and the Euler grid using an interpolation function; wherein the interpolation function is used to convert the force on the solid boundary into the force on the fluid grid point, and to transfer the reaction force of the fluid on the solid boundary from the fluid grid point to the Lagrangian point.

[0207] Based on the above embodiment, in a feasible embodiment, the association subunit includes:

[0208] The execution subunit is used to initialize the physical parameters and initial positions of the fluid domain and solid particles, construct the Euler grid of the fluid domain and the Lagrange points of the solid boundary; use the Navier-Stokes equations to solve the velocity field and pressure field of the fluid domain; calculate the force on the solid boundary, apply the force to the fluid domain, and update the position and velocity of the Lagrange point; use the interpolation function to transfer the force of the solid boundary to the fluid grid point, and calculate the reaction force of the fluid on the solid boundary; update the velocity field of the fluid domain and the position and velocity of the solid particles through the flow field in the fluid domain and the motion equation of the solid domain;

[0209] The loop subunit is used to iteratively execute the subunit and perform iterative calculation of the time step until the iteration time meets the predetermined simulation time or the convergence condition is met and the iteration is stopped.

[0210] Based on the above embodiment, in a feasible embodiment, the sedimentation velocity calculation module includes:

[0211] A generating submodule, used to generate a force model of ore particles; the force model includes the interaction of gravity, buoyancy and fluid resistance;

[0212] A first calculation submodule, configured to calculate the settling velocity of the ore particles according to the gravity, the buoyancy and the fluid resistance;

[0213] The determination submodule is used to determine the settling velocity of the ore particle group according to the overall resistance coefficient and the average diameter of the ore particle group.

[0214] Based on the above embodiment, in a feasible embodiment, the mining planning module includes:

[0215] The second calculation submodule is used to calculate the minimum lifting speed of the ore particles according to the settling speed and the design parameters of the hydraulic lifting system.

[0216] The present application also provides a computer-readable storage medium on which a computer program is stored, and when the computer program is executed, the steps provided in the above embodiment can be implemented. The storage medium may include various media that can store program codes, such as a USB flash drive, a mobile hard disk, a read-only memory (ROM), a random access memory (RAM), a magnetic disk or an optical disk.

[0217] The present application also provides an electronic device, see Figure 6 , a structural diagram of an electronic device provided in an embodiment of the present application, such as Figure 6 As shown, a processor 1410 and a memory 1420 may be included.

[0218] Among them, the processor 1410 may include one or more processing cores, such as a 4-core processor, an 8-core processor, etc. The processor 1410 may be implemented in at least one hardware form of DSP (Digital Signal Processing), FPGA (Field-Programmable Gate Array), and PLA (Programmable Logic Array). The processor 1410 may also include a main processor and a coprocessor. The main processor is a processor for processing data in the awake state, also known as a CPU (Central Processing Unit); the coprocessor is a low-power processor for processing data in the standby state. In some embodiments, the processor 1410 may be integrated with a GPU (Graphics Processing Unit), which is responsible for rendering and drawing the content to be displayed on the display screen. In some embodiments, the processor 1410 may also include an AI (Artificial Intelligence) processor, which is used to process computing operations related to machine learning.

[0219] The memory 1420 may include one or more computer-readable storage media, which may be non-transitory. The memory 1420 may also include a high-speed random access memory, and a non-volatile memory, such as one or more disk storage devices, flash memory storage devices. In this embodiment, the memory 1420 is at least used to store the following computer program 1421, wherein, after the computer program is loaded and executed by the processor 1410, it can implement the relevant steps in the method performed by the electronic device side disclosed in any of the aforementioned embodiments. In addition, the resources stored in the memory 1420 may also include an operating system 1422 and data 1423, etc., and the storage method may be temporary storage or permanent storage. Among them, the operating system 1422 may include Windows, Linux, Android, etc.

[0220] In some embodiments, the electronic device may further include a display screen 1430 , an input / output interface 1440 , a communication interface 1450 , a sensor 1460 , a power source 1470 , and a communication bus 1480 .

[0221] certainly, Figure 6 The structure of the electronic device shown does not constitute a limitation on the electronic device in the embodiments of the present application. In actual applications, the electronic device may include Figure 6 More or fewer components than shown, or combinations of certain components.

[0222] The various embodiments in the specification are described in a progressive manner, and each embodiment focuses on the differences from other embodiments. The same or similar parts between the various embodiments can be referred to each other. For the system provided in the embodiment, since it corresponds to the method provided in the embodiment, the description is relatively simple, and the relevant parts can be referred to the method part description.

[0223] Specific examples are used herein to illustrate the principles and implementation methods of the present application, and the description of the above embodiments is only used to help understand the method and core ideas of the present application. It should be pointed out that for ordinary technicians in this technical field, without departing from the principles of the present application, several improvements and modifications can be made to the present application, and these improvements and modifications also fall within the scope of protection of the claims of the present application.

[0224] It should also be noted that, in this specification, relational terms such as first and second, etc. are only used to distinguish one entity or operation from another entity or operation, and do not necessarily require or imply any such actual relationship or order between these entities or operations. Moreover, the terms "comprise", "include" or any other variants thereof are intended to cover non-exclusive inclusion, so that a process, method, article or device including a series of elements includes not only those elements, but also other elements not explicitly listed, or also includes elements inherent to such process, method, article or device. In the absence of further restrictions, an element defined by the statement "comprises a ..." does not exclude the presence of other identical elements in the process, method, article or device including the element.

Claims

1. A mining planning method, characterized in that: include: Obtain the physical characteristics and initial parameters of the ore particle group; Performing numerical simulation on the ore particle group according to the physical characteristics and the initial parameters, calculating the resistance coefficient of the ore particle group at different flow rates, and generating a mathematical relationship between the resistance coefficient and the flow rate change; Calculating the settling velocity of the ore particle group according to the mathematical relationship; The minimum lifting speed of the ore particles is calculated according to the settling speed; the minimum lifting speed is used to guide mining planning.

2. The mining planning method according to claim 1, characterized in that: The numerical simulation of the ore particle group according to the physical characteristics and the initial parameters comprises: Constructing an initial geometric model of the ore particle group according to the physical characteristics and the initial parameters; the initial geometric model includes ore particles of different sizes and shapes; A coupling model of the fluid domain and the solid domain is established based on the initial geometric model.

3. The mining planning method according to claim 2, characterized in that: The step of establishing a coupling model of a fluid domain and a solid domain based on the initial geometric model comprises: The Navier-Stokes equations are used to describe the movement of ore particles in the fluid domain; Using Newton's second law to describe the solid domain motion of the ore particles; A coupling model is established according to the fluid domain motion and the solid domain motion.

4. The mining planning method according to claim 3, characterized in that: Establishing a coupling model according to the fluid domain motion and the solid domain motion includes: Define Lagrange points and Eulerian grids; The Lagrangian point and the Euler grid are associated by using an interpolation function, wherein the interpolation function is used to convert the force on the solid boundary into the force on the fluid grid point, and to transfer the reaction force of the fluid on the solid boundary from the fluid grid point to the Lagrangian point.

5. The mining planning method according to claim 4, characterized in that: The associating the Lagrange point and the Euler grid by using an interpolation function comprises: Step A: Initialize the physical parameters and initial positions of the fluid domain and solid particles, and construct the Euler grid of the fluid domain and the Lagrange points of the solid boundary; Step B: Use the Navier-Stokes equations to solve the velocity and pressure fields in the fluid domain; Step C: Calculate the force on the solid boundary, apply the force to the fluid domain, and update the position and velocity of the Lagrangian point; Step D: using the interpolation function to transfer the force of the solid boundary to the fluid grid point, and calculating the reaction force of the fluid on the solid boundary; Step E: Update the velocity field in the fluid domain and the position and velocity of the solid particles through the flow field in the fluid domain and the motion equation of the solid domain; Iteration step A to step E is to perform iterative calculation of the time step, and the iteration is stopped when the iteration time meets the predetermined simulation time or meets the convergence condition.

6. The mining planning method according to claim 1, characterized in that: Calculating the settling velocity of the ore particle group according to the mathematical relationship includes: Generate a force model of the ore particles; the force model includes the interaction of gravity, buoyancy and fluid resistance; Calculating the settling velocity of the ore particles according to the gravity, the buoyancy and the fluid resistance; The settling velocity of the ore particle group is determined based on the overall resistance coefficient and the average diameter of the ore particle group.

7. The mining planning method according to claim 1, characterized in that: Calculation of the minimum lifting velocity of ore particles based on the settling velocity includes: The minimum lifting velocity of the ore particles is calculated based on the settling velocity and the design parameters of the hydraulic lifting system.

8. A mining planning system, characterized in that: include: A data acquisition module, used to obtain the physical characteristics and initial parameters of the ore particle group; A numerical simulation module, used for performing numerical simulation on the ore particle group according to the physical characteristics and the initial parameters, calculating the resistance coefficient of the ore particle group at different flow rates, and generating a mathematical relationship between the resistance coefficient and the flow rate change; A settling velocity calculation module, used to calculate the settling velocity of the ore particle group according to the mathematical relationship; a mining planning module for calculating a minimum lifting velocity of ore particles based on the settling velocity; The minimum lifting speed is used to guide mining planning.

9. An electronic device, characterized in that: Including representation Memory for storing computer programs; A processor, configured to implement the steps of the method according to any one of claims 1 to 7 when executing the computer program.

10. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores a computer program, which implements the steps of the method according to any one of claims 1 to 7 when executed.

Images