A mining planning method, system, storage medium and device

Abstract

The present application provides a mining planning method, including: obtaining the 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 drag coefficient of the ore particle group at different flow rates, and generating a mathematical relationship between the drag 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; and using the minimum lifting velocity to guide mining planning. The present application accurately calculates the drag coefficient and settling velocity of the particle group through numerical simulation, and then determines the minimum lifting velocity, ensuring that the ore particles can be smoothly lifted, and providing an important reference basis for the design and optimization of the deep-sea mining hydraulic lifting system. The present application also provides a mining planning system, a computer-readable storage medium, and an electronic device, which have the above beneficial effects.

Description

Technical Field

[0001] This application relates to the field of deep - sea mining, and particularly 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 sea surface through a hydraulic lifting system. To ensure that ore particles can be smoothly lifted, it is necessary to determine their minimum lifting speed. Traditional methods usually focus on the drag coefficient and settling velocity of a single particle, 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 needs to be urgently solved by those skilled in the art. 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. By numerically simulating to calculate the drag coefficient of a group of ore particles and combining with the settling velocity of the ore particle group, the minimum lifting speed is finally calculated, providing a reliable basis for setting hydraulic lifting parameters in the deep - sea mining process.

[0004] To solve the above - mentioned 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] Perform numerical simulation on the ore particle group according to the physical characteristics and the initial parameters, calculate the drag coefficient of the ore particle group at different flow velocities, and generate a mathematical relationship between the drag coefficient and the flow velocity change;

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

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

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

[0010] 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;

[0011] Establish a coupled model of the fluid domain and the solid domain based on the initial geometric model.

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

[0013] Use the Navier-Stokes equations to describe the motion of ore particles in the fluid domain;

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

[0015] Establish a coupling model based on the motion of the fluid domain and the motion of the solid domain.

[0016] Optionally, establishing a coupling model based on the motion of the fluid domain and the motion of the solid domain includes:

[0017] Define Lagrangian points and Eulerian grids;

[0018] Use an interpolation function to correlate the Lagrangian points and the Eulerian grids; wherein, the interpolation function is used to convert the force acting on the solid boundary into the force acting on the fluid grid points, and to transfer the reaction force of the fluid on the solid boundary from the fluid grid points to the Lagrangian points.

[0019] Optionally, using the interpolation function to correlate the Lagrangian points and the Eulerian grids includes:

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

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

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

[0023] Step D: Use the interpolation function to transfer the force on the solid boundary to the fluid grid points and calculate the reaction force of the fluid on the solid boundary;

[0024] Step E: 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 equations of the solid domain;

[0025] Iterate steps A to E, perform iterative calculations of the time step until the iteration duration meets the predetermined simulation time or meets the convergence condition and stop the iteration.

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

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

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

[0029] Determine the settling velocity of the ore particle group according to the overall drag coefficient and average diameter of the ore particle group.

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

[0031] Calculate the minimum lifting velocity of the ore particles according to the settling velocity and the design parameters of the hydraulic lifting system.

[0032] This application also provides a mining planning system, including:

[0033] A data acquisition module, configured to acquire the physical properties and initial parameters of the ore particle group;

[0034] A numerical simulation module, configured to perform numerical simulation on the ore particle group according to the physical properties and the initial parameters, calculate the drag coefficient of the ore particle group at different flow velocities, and generate a mathematical relationship between the drag coefficient and the flow velocity change;

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

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

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

[0038] This application also provides an electronic device, including a memory and a processor, where 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] This application provides a mining planning method, including: acquiring the physical properties and initial parameters of the ore particle group; performing numerical simulation on the ore particle group according to the physical properties and the initial parameters, calculating the drag coefficient of the ore particle group at different flow velocities, and generating a mathematical relationship between the drag coefficient and the flow velocity 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 accurately calculates the drag coefficient and settling velocity of the particle group through numerical simulation, and then determines the minimum lifting velocity to ensure that the ore particles can be smoothly lifted, providing an important reference basis for the design and optimization of the deep-sea mining hydraulic lifting system.

[0041] The present application also provides a mining planning system, a computer-readable storage medium, and an electronic device, which have the above beneficial effects and will not be elaborated herein. BRIEF DESCRIPTION OF THE DRAWINGS

[0042] In order to more clearly illustrate the technical solutions in the embodiments of the present application or the prior art, the following will briefly introduce the drawings required for use in the description of the embodiments or the prior art. Obviously, the drawings in the following description are only the embodiments of the present application. For those of ordinary skill in the art, without creative efforts, other drawings can also be obtained based on the provided drawings.

[0043] Figure 1 It is a flowchart of a mining planning method provided by an embodiment of the present application;

[0044] Figure 2 It is a flowchart of numerical simulation using the immersed boundary method provided by an embodiment of the present application;

[0045] Figure 3 It is a schematic diagram of an initial set model of an ore particle group provided by an embodiment of the present application;

[0046] Figure 4 It is a flowchart of another mining planning method provided by an embodiment of the present application;

[0047] Figure 5 It is a schematic diagram of the structure of a mining planning system provided by an embodiment of the present application;

[0048] Figure 6 It is a structural diagram of an electronic device provided by an embodiment of the present application. DETAILED DESCRIPTION OF THE EMBODIMENTS

[0049] To make the objectives, technical solutions, and advantages of the embodiments of the present application clearer, the following will clearly and completely describe the technical solutions in the embodiments of the present application with reference to the accompanying drawings in the embodiments of the present application. Obviously, the described embodiments are some, but not all, of the embodiments of the present application. Based on the embodiments in the present application, all other embodiments obtained by those of ordinary skill in the art without creative efforts fall within the scope of protection of the present application.

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

[0051] SeeFigure 1 , Figure 1 is a flowchart of a mining planning method provided by an embodiment of the present application. The method includes:

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

[0053] S102: Perform numerical simulation on the ore particle group according to the physical properties and the initial parameters, calculate the drag coefficient of the ore particle group at different flow velocities, and generate a mathematical relationship between the drag coefficient and the flow velocity change;

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

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

[0056] First, determine the physical properties and initial parameters of the ore particle group. 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] Use the immersed boundary method to perform direct numerical simulation of the particle group. The specific process of numerical simulation by the immersed boundary method is as follows:

[0058] The first step: Construct the initial geometric model of the particle group:

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

[0060] The second step: 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 the key part. This process involves transferring the force on 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, including the momentum equation and the continuity equation.

[0063] The momentum equation is expressed as follows:

[0064]

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

[0066] The continuity equation (mass conservation) 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, mainly considering the force balance and expressed as:

[0069]

[0070] where 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), and F b is the interaction force between solids (N).

[0071] The immersed boundary method realizes the coupling of the fluid domain and the solid domain through the interaction between Lagrangian points and Eulerian grids. The specific process is as follows:

[0072] First, define the Lagrangian points and Eulerian grids: Lagrangian points are used to represent the solid boundary, and Eulerian grids are used to describe the fluid domain. The two are connected through an interpolation function. In the following text, X s represents the coordinates of the Lagrangian points on the solid boundary, and X i represents the coordinates of the fluid grid points.

[0073] After that, define the interpolation function. To realize the interaction between Lagrangian points and Eulerian grids, define the interpolation function to represent

[0074]

[0075] Through the interpolation function, the force F s exerted on the solid boundary is applied to the fluid domain and converted into the force f i at the fluid grid points, expressed as:

[0076]

[0077] Similarly, through the interpolation function, the reaction force of the fluid on the solid boundary is transferred from the fluid grid points to the Lagrangian points, expressed as:

[0078]

[0079] Subsequently, coupling calculations are performed. Associating the Lagrangian points and the Eulerian grid using interpolation functions 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 Eulerian grid of the fluid domain and the Lagrangian points of the solid boundary;

[0081] Step B: Solve the velocity field and pressure field of the fluid domain using the Navier-Stokes equations;

[0082] Step C: Calculate the acting forces on the solid boundary, apply the acting forces to the fluid domain, and update the positions and velocities of the Lagrangian points;

[0083] Step D: Transfer the acting forces on the solid boundary to the fluid grid points using the interpolation functions, and calculate the reaction forces of the fluid on the solid boundary;

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

[0085] Iterate steps A to E, perform iterative calculations of the time step, and stop the iteration until the iteration duration meets the predetermined simulation time or the convergence condition is satisfied.

[0086] In step A, the physical parameters and initial positions of the fluid domain and solid particles are initialized. The Eulerian grid 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 field and pressure field of the fluid domain.

[0088] In step C, the acting force F on the solid boundary is calculated s , and it is applied to the fluid domain. Subsequently, the positions and velocities of the Lagrangian points are updated.

[0089] In step D, the acting forces on the solid boundary are transferred to the fluid grid points using the interpolation functions, and the reaction forces of the fluid on the solid boundary are calculated.

[0090] See Figure 2 , Figure 2 , which is the flowchart of numerical simulation using the immersed boundary method provided by the embodiments of the present application. Figure 2 Shows the complete process of numerical simulation using the immersed boundary method, including the following process:

[0091] First, construct an initial geometric model. The initial geometric model of the ore particle group includes multiple ore particles of different sizes and shapes.

[0092] Secondly, a coupled model of the fluid domain and the solid domain is established: In the fluid domain, the Navier-Stokes equation is used to describe the fluid motion. In the solid domain, Newton's second law is used to describe the motion of the particles. The coupling of the fluid domain and the solid domain is achieved through the interaction between Lagrangian points and Eulerian grids.

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

[0094] When updating the fluid and solid domains, the velocity field of the fluid domain and the positions and velocities of the solid particles are updated.

[0095] Iterative calculations are performed until convergence or the final time step is reached.

[0096] Repeat the above steps for iterative calculations of the time step until the iterative duration meets the predetermined final time step or the convergence condition, i.e., perform iterative calculations of the time step until the predetermined simulation time or the convergence condition is reached.

[0097] An exemplary calculation process is as follows:

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

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

[0100] Particle diameter d = 0.01 m;

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

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

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

[0104] Calculate the initial fluid velocity field 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 computational fluid on the solid boundary is as follows:

[0111]

[0112] By iterating the above steps repeatedly, 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 numerically simulating the ore particle swarm, the drag coefficient of the particle swarm is calculated. Here, 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 (N) on the particle swarm, ρ f represents the fluid density (kg / m 3 ), u represents the fluid velocity (m / s), A represents the frontal area of the particle swarm (m 2 ), and Vs represents the particle velocity.

[0118] By numerically simulating the drag coefficients at different flow velocities, the mathematical relationship between the drag coefficient and the flow velocity change can be determined. For example, a curve graph of the drag coefficient varying with the flow velocity can be plotted for subsequent calculation of the settling velocity of the particle swarm.

[0119] Calculating the settling velocity of the particle swarm is one of the key steps in determining the minimum lifting velocity. The settling velocity reflects the movement velocity of the particle swarm after the force balance in a static or flowing fluid. The process of calculating the settling velocity is described in detail below, and the relevant equations and physical models are listed. Specifically, it can include the following steps:

[0120] First step, generate the force model of the ore particles; the force model includes the interaction of gravity, buoyancy, and fluid drag;

[0121] Second step, calculate the settling velocity of the ore particles according to the gravity, the buoyancy, and the fluid drag;

[0122] Step 3: Determine the sedimentation velocity of the ore particle group according to the overall resistance coefficient and average diameter of the ore particle group.

[0123] In a fluid, a particle group is subjected to the actions of 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] where F g is gravity, F b is buoyancy, and F d is fluid resistance.

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

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

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

[0129]

[0130] where 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] The fluid resistance F d is usually related to the velocity of the fluid and the shape of the particle, and can be expressed by the resistance coefficient C d :

[0132]

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

[0134]

[0135] The resistance coefficient C d is obtained through numerical simulation and represents the resistance coefficient of the particle group at different flow velocities.

[0136] The sedimentation velocity V tIt is the velocity under the condition of force balance. By balancing gravity, buoyancy, and resistance, the calculation formula for the settling velocity can be obtained:

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

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

[0139]

[0140] After arrangement, we get:

[0141]

[0142] The final settling velocity formula is:

[0143]

[0144] For a particle swarm, the settling velocity needs to consider the interaction between particles and the fluid flow characteristics. The settling velocity of the particle swarm can be calculated through the overall drag coefficient C d of the particle swarm and the average diameter :

[0145]

[0146] In the above formula, is the average diameter of the particle swarm, which is obtained through statistical analysis.

[0147] Through numerical simulation by the immersed boundary method, the drag coefficient C d of the particle swarm at different flow velocities is obtained, and combined with the physical characteristics of the particle swarm, the settling velocity is calculated. The specific steps are as follows:

[0148] Step 1: Determine the initial parameters of the particle swarm: including the density ρ s of the particles, the average diameter , and the fluid density ρ f .

[0149] Step 2: Numerically simulate to obtain the drag coefficient: Through numerical simulation by the immersed boundary method, the drag coefficient C d of the particle swarm at different flow velocities is obtained.

[0150] Step 3: Calculate the settling velocity: Using the above settling velocity formula and combining with C d obtained from the numerical simulation, the settling velocity of the particle swarm can be calculated.

[0151] To better understand the above calculation process, assume the parameters of a certain ore particle swarm are as follows:

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

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

[0154]

[0155] C obtained from numerical simulation d = 0.05;

[0156] Substitute into the sedimentation velocity formula:

[0157]

[0158] Calculated as:

[0159] V t ≈ 0.9 m / s;

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

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

[0162] In the hydraulic lifting process of deep - sea mining, determining 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 velocity is usually set to three times the sedimentation velocity (Settling Velocity, V t ) of the particle swarm. This empirical formula is derived from a comprehensive consideration of multiple factors, including hydrodynamics, particle - swarm interactions, and the safety margin of actual operation.

[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 particle movement. Setting the minimum lifting velocity to three times the sedimentation velocity helps to ensure that all particles can be stably lifted under turbulent conditions and will not settle due to instantaneous velocity fluctuations.

[0164] The particles in the particle swarm are not only affected by gravity and fluid resistance but also have collision and agglomeration phenomena. These interactions increase the complexity of particle lifting. By increasing the minimum lifting velocity, the adverse effects of particle - particle interactions on the overall lifting effect can be reduced, ensuring the overall upward movement of the particle swarm.

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

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

[0167] See Figure 2 , Figure 2 which is a schematic diagram of the initial set model of the ore particle group provided by the embodiment of the present application. The appendix Figure 1 shows a schematic diagram of the initial geometric model of the ore particle group, demonstrating the distribution of particles of different sizes and shapes in the fluid. Through Figure 2 it can be seen the mechanical actions exerted on the ore particle group in the fluid and its sedimentation and lifting processes.

[0168] The appendix Figure 3 shows a schematic diagram of the initial geometric model of the ore particle group. 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 of d A , and the shape factor

[0170] For particle B: It is ellipsoidal with a major axis of a B and a minor axis of b B , and the shape factor

[0171] For particle C: It is polyhedral with a characteristic side length of l C , and the shape factor

[0172] Figure 2 In the appendix, the particles are randomly distributed in the fluid. Considering the complexity of the particle group in actual situations, a random distribution method is adopted. There is a certain distance between the particles to avoid collisions and interactions in the initial state.

[0173] Figure 2 In the appendix, a three-dimensional coordinate system is used for representation, 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 (x i , y i , z i ) of its center point.

[0174] The fluid boundary condition is assumed to be a no-slip boundary condition, that is, the velocity of the fluid at the particle surface is the same as the particle velocity. The initial flow direction and velocity distribution are set in the fluid domain to simulate the motion of the particle swarm in the fluid.

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

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

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

[0178] When defining the fluid domain, a rectangular region is adopted for the fluid domain, and the boundary conditions and initial flow velocity distribution are set. The initial geometric model and fluid domain parameters are imported through numerical simulation software.

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

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

[0181] Determine the physical properties and initial parameters of the ore particle swarm. 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] Construct the initial geometric model. Construct the geometric model of the ore particle swarm, including particles of different sizes and shapes, placed in the fluid domain.

[0183] Establish a coupled model of the fluid domain and the solid domain. The Navier - Stokes equation is used to describe the fluid flow in the fluid domain. Newton's second law is used to describe the particle motion in the solid domain. Through the interaction between Lagrangian points and Eulerian grids, the coupling of the fluid and the solid is realized.

[0184] Numerically calculate the flow field and the forces acting on it. Solve the velocity field and pressure field in the fluid domain through numerical simulation. Calculate the forces acting on the solid boundary and apply them to the fluid domain to update the positions and velocities of the particles.

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

[0186] Iteratively calculate until convergence or reach the final time step. Repeat the above steps and perform iterative calculations through time steps until the predetermined simulation time or convergence condition is reached.

[0187] Calculate the drag coefficient of the particle swarm. Calculate the drag coefficient of the particle swarm (the ratio of the force acting on the particle swarm in the fluid to the fluid flow velocity) according to the flow velocity obtained from the numerical simulation and the forces acting on the particles.

[0188] Calculate the settling velocity of the particle swarm. Calculate the settling velocity of the particle swarm (the movement velocity of the particle swarm after force balance in a stationary or flowing fluid) through the drag coefficient of the particle swarm and the fluid properties.

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

[0190] Verify and optimize the minimum lifting velocity. Verify the accuracy of the calculation results through actual operation and experiments. According to 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, accurate simulation of the dynamic characteristics of the ore particle swarm during the hydraulic lifting process can be achieved, thereby calculating the minimum lifting velocity. Based on the minimum lifting velocity, it is used to guide the mining plan and provide a reliable basis for setting the hydraulic lifting parameters in the deep - sea mining process.

[0192] See Figure 5 , Figure 5 which is the schematic structural diagram of the mining planning system provided by the embodiment of the present application. The present application also provides a mining planning system, including:

[0193] A data acquisition module for acquiring the physical properties and initial parameters of the ore particle swarm;

[0194] A numerical simulation module for performing numerical simulation on the ore particle swarm according to the physical properties and the initial parameters, calculating the drag coefficient of the ore particle swarm at different flow velocities, and generating the mathematical relationship between the drag coefficient and the flow velocity change;

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

[0196] A mining planning module, configured to calculate the minimum lifting velocity of the ore particles according to the settlement velocity; the minimum lifting velocity is used to guide mining planning.

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

[0198] A first construction sub-module, configured to construct an initial geometric model of the ore particle group according to the physical properties and the initial parameters; the initial geometric model includes ore particles of different sizes and shapes;

[0199] A second construction sub-module, configured to establish a coupled model of the fluid domain and the solid domain based on the initial geometric model.

[0200] Based on the above embodiments, in a feasible embodiment, the second construction sub-module includes:

[0201] A first description unit, configured to describe the motion of the fluid domain of the ore particles by using the Navier-Stokes equation;

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

[0203] A construction unit, configured to establish a coupled model according to the motion of the fluid domain and the motion of the solid domain.

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

[0205] A definition sub-unit, configured to define Lagrangian points and Eulerian grids;

[0206] An association sub-unit, configured to associate the Lagrangian points and the Eulerian grids by using an interpolation function; wherein, the interpolation function is used to convert the acting force on the solid boundary into the acting force on the fluid grid points, and to transfer the reaction force of the fluid on the solid boundary from the fluid grid points to the Lagrangian points.

[0207] Based on the above embodiments, in a feasible embodiment, the association sub-unit includes:

[0208] An execution subunit, configured to initialize the physical parameters and initial positions of the fluid domain and solid particles, construct the Eulerian grid of the fluid domain and the Lagrangian points of the solid boundary; solve the velocity field and pressure field of the fluid domain using the Navier-Stokes equations; calculate the acting forces on the solid boundary, apply the acting forces to the fluid domain, and update the positions and velocities of the Lagrangian points; transfer the acting forces of the solid boundary to the fluid grid points using the interpolation function, and calculate the reaction forces of the fluid on the solid boundary; update the velocity field of the fluid domain and the positions and velocities of the solid particles through the flow field in the fluid domain and the motion equations of the solid domain.

[0209] A loop subunit, configured to iterate the execution subunit to perform iterative calculations of time steps until the iteration duration meets the predetermined simulation time or the convergence condition, and then stop the iteration.

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

[0211] A generation sub-module, configured to generate a force model for ore particles; the force model includes the interaction of gravity, buoyancy, and fluid resistance.

[0212] A first calculation sub-module, configured to calculate the settlement velocity of the ore particles according to the gravity, the buoyancy, and the fluid resistance.

[0213] A determination sub-module, configured to determine the settlement velocity of the ore particle group according to the overall resistance coefficient and average diameter of the ore particle group.

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

[0215] A second calculation sub-module, configured to calculate the minimum lifting velocity of the ore particles according to the settlement velocity 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. When the computer program is executed, the steps provided in the above embodiments can be implemented. The storage medium may include various media that can store program codes, such as USB flash drives, external hard drives, read-only memories (ROMs), random access memories (RAMs), magnetic disks, or optical discs.

[0217] The present application also provides an electronic device. Refer to Figure 6 , the structural diagram of an electronic device provided in the embodiments of the present application, as shown in Figure 6 , which may include a processor 1410 and a memory 1420.

[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 used to process data in the wake state, also known as the CPU (Central Processing Unit); the coprocessor is a low-power processor used to process data in the standby state. In some embodiments, the processor 1410 may be integrated with a GPU (Graphics Processing Unit), and the GPU is responsible for rendering and drawing the content to be displayed on the display screen. In some embodiments, the processor 1410 may further include an AI (Artificial Intelligence) processor, and the AI processor is used to process computational operations related to machine learning.

[0219] The memory 1420 may include one or more computer-readable storage media, and the computer-readable storage media may be non-transitory. The memory 1420 may further include high-speed random access memory and non-volatile memory, such as one or more disk storage devices and flash storage devices. In this embodiment, the memory 1420 is at least used to store the following computer program 1421. After the computer program is loaded and executed by the processor 1410, it can implement the relevant steps in the method executed by the electronic device side disclosed in any of the foregoing embodiments. In addition, the resources stored in the memory 1420 may further 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 supply 1470, and a communication bus 1480.

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

[0222] The various embodiments in the specification are described in a progressive manner. Each embodiment focuses on the differences from other embodiments. For the same or similar parts among the embodiments, reference can be made 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. For the relevant parts, reference can be made to the description in the method section.

[0223] Specific examples are used in this article to elaborate on the principles and implementation manners of the present application. The descriptions of the above embodiments are only used to help understand the method of the present application and its core idea. It should be noted that for those of ordinary skill in the art in this technical field, without departing from the principle of the present application, several improvements and modifications can be made to the present application, and these improvements and modifications also fall within the protection scope of the claims of the present application.

[0224] It should also be noted that in this specification, relational terms such as "first" and "second" are only used to distinguish one entity or operation from another entity or operation, and do not necessarily require or imply any actual relationship or order between these entities or operations. Moreover, the term "comprising", "including" or any other variant thereof is intended to cover non-exclusive inclusion, so that a process, method, article or device including a series of elements not only includes those elements, but also includes other elements not expressly listed, or further includes elements inherent to such process, method, article or device. Without further limitation, an element defined by the statement "including a..." does not exclude the existence of additional identical elements in the process, method, article or device including the said 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; Determining the settling velocity of the ore particle group according to the overall resistance coefficient and the average diameter of the ore particle group; Wherein, the force model is expressed as: F g -F b -F d =0; Among them, F g is the gravity, F b is the buoyancy, F d is the fluid resistance.

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