A method and system for dynamic monitoring of ion-adsorption rare earth ore leaching process

Abstract

This invention discloses a method and system for dynamic monitoring of the leaching process of ion-adsorption rare earth ores, relating to the field of leaching monitoring technology. The dynamic monitoring system for the leaching process of ion-adsorption rare earth ores includes a leaching monitoring module and a leaching adjustment module. This invention comprehensively analyzes multi-dimensional data such as ore pore ion liquid strength, ore pore strain, and ion exchange reactions to monitor the ore body state in real time. It adjusts the leaching injection rate according to dynamic changes, optimizing the leaching process to achieve higher leaching efficiency and reach the optimal leaching effect in a shorter time. Through a dynamic game evolution model, it automatically adjusts the injection threshold based on the ore body stability parameters and leaching state characteristics, ensuring efficient leaching while maintaining the ore body structure stability. This avoids damage to the ore body from excessive injection and effectively improves the long-term stability of the ore body.

Description

Technical Field

[0001] This invention relates to the field of leaching monitoring technology, and in particular to a method and system for dynamic monitoring of the leaching process of ion-type rare earth minerals. Background Technology

[0002] Traditional rare earth ore leaching processes typically rely on manual monitoring and experience-based control, estimating leaching efficiency and progress by measuring data such as ore pore ion liquid strength and ore pore strain. However, existing technologies often suffer from slow response speed, low monitoring accuracy, and difficulty in finding the optimal injection rate while ensuring ore body safety. As the demands for ore leaching efficiency and safety increase, the limitations of traditional methods are becoming increasingly apparent, necessitating an efficient, automated, and dynamically adjustable leaching process monitoring and optimization scheme. Summary of the Invention

[0003] The present invention aims to provide a dynamic monitoring method and system for the leaching process of ion-type rare earth minerals. Through dynamic monitoring and intelligent adjustment, the optimal injection parameters can be found to improve leaching efficiency and ensure the safety of the ore body at all times.

[0004] A method for dynamic monitoring of the leaching process of ion-adsorption rare earth ores includes the following steps: Within the rare earth ore leaching area, real-time monitoring is conducted to acquire ore pore ion liquid strength data and ore pore strain data. Based on the ore pore ion liquid strength data and ore pore strain data, an interlayer mutual exclusion mapping function is constructed. The interlayer expansion pressure data is determined according to the interlayer mutual exclusion mapping function. By comparing the spatiotemporal evolution data of the interlayer expansion pressure data and the ore pore strain data, the stability parameter characteristics of the ore body are obtained. The complex impedance spectrum distribution inside the ore body is scanned to obtain the ore body pore phase angle data; the ion exchange reaction state inside the ore body is identified based on the ore body pore phase angle data to obtain the ore body ion exchange characteristics; the leaching process is judged based on the ore body ion exchange characteristics to obtain the ore body leaching state characteristics. Based on the stability parameters and leaching state characteristics of the ore body, a game-theoretic evolution is performed to obtain the leaching injection threshold matrix corresponding to the current ore body. The ion liquid injection parameters are adjusted according to the leaching injection threshold matrix to obtain new ore pore ion liquid intensity data.

[0005] As a preferred embodiment of the present invention, the specific steps for comparing the spatiotemporal evolution data of interlayer expansion pressure data and ore pore strain data include: Based on the ore pore ion liquid strength data, a double electric layer osmotic pressure function is constructed to represent the repulsive interaction between like-electric particles within the ore body; the effective stress of the ore body is obtained; based on the ore pore strain data and the effective stress of the ore body, a structural stiffness calibration coefficient for the ore body is constructed; the double electric layer osmotic pressure function is corrected according to the structural stiffness calibration coefficient of the ore body to obtain the mutual repulsion mapping function between leaching layers; The ore pore ion liquid intensity data is input into the interlayer mutual repulsion mapping function for calculation to obtain the particle repulsion increment of the ore body; based on the particle repulsion increment of the ore body, the interlayer expansion pressure data is obtained; the interlayer expansion pressure data and ore pore strain data are subjected to trend correlation analysis in the same time domain coordinate system to obtain the internal strain components of the ore body. A three-dimensional space inside the ore body is constructed; the strain components inside the ore body are calculated in the three-dimensional space to obtain the local gradient difference value of the ore body; based on the local gradient difference value of the ore body, the shear slip potential zone of the ore body caused by non-uniform chemical expansion is identified to obtain the deformation region of the ore body; the deformation region of the ore body is calculated to obtain the stability parameter characteristics of the ore body.

[0006] As a preferred embodiment of the present invention, the specific steps for identifying the ion exchange reaction state inside the ore body based on the ore body pore phase angle data include: Extract the corresponding phase angle peak and corresponding phase angle characteristic frequency from the ore body pore phase angle data; map the real-time coordinates of the corresponding phase angle peak and corresponding phase angle characteristic frequency into the ion exchange standard feature space; The response state of ion exchange is determined within the standard characteristic space of ion exchange to obtain the ion exchange phase state of the ore body; The distribution of ion exchange phases in the ore body at different monitoring times was statistically analyzed to obtain the evolution gradient of ion exchange saturation corresponding to different levels within the ore body; the ion exchange characteristics of the ore body were generated by combining the three-dimensional space within the ore body and the evolution gradient of ion exchange saturation.

[0007] As a preferred embodiment of the present invention, the specific steps for judging the leaching process based on the ion exchange characteristics of the ore body include: Spatial interpolation algorithms are used to connect the ion exchange characteristics of the ore body at different monitoring times to obtain the reaction envelope diagram of the leaching process. By comparing the spatial displacement vectors at different monitoring times in the reaction envelope diagram of the leaching process with the preset theoretical liquid flow propagation vector, the uniformity of leaching seepage is obtained. The leaching state is identified based on the uniformity of leaching seepage, and the leaching state characteristics of the ore body are obtained.

[0008] As a preferred technical solution of the present invention, the specific steps for game evolution based on the stability parameter characteristics and leaching state characteristics of the ore body include: A game theory model for leaching evolution is constructed, in which the first participant is the entity responsible for leaching efficiency recovery, and the second participant is the entity responsible for leaching safety constraints. The leaching safety penalty factor is calculated based on the stability parameter characteristics of the ore body; the leaching retention loss factor is calculated based on the retention loss component of the ore body leaching state characteristics; the leaching safety penalty factor and the leaching retention loss factor are input into the leaching evolution game model for evolution to obtain the leaching injection threshold matrix. The specific steps of evolution are as follows: A displacement-concentration adjustment operator is introduced. Within each evolution cycle, the fitness function is maximized by maximizing the uniformity of the three-dimensional space within the ore body. Within the evolution cycle, the displacement-concentration adjustment operator is adjusted. When the change in the displacement-concentration adjustment operator causes the growth rate of the leaching safety penalty factor to exceed the growth rate of the leaching retention loss factor, a first adjustment strategy is implemented for the displacement-concentration adjustment operator. When the change in the displacement-concentration adjustment operator causes the growth rate of the leaching retention loss factor to exceed the growth rate of the leaching safety penalty factor, a second adjustment strategy is implemented for the displacement-concentration adjustment operator.

[0009] As a preferred embodiment of the present invention, the specific steps for adjusting the ion liquid injection parameters according to the leaching injection threshold matrix include: The threshold matrix of leaching injection is analyzed to extract the ion liquid injection parameters. The ion liquid injection parameters include the ion liquid injection pressure boundary, chemical potential energy gradient coefficient and shear strain buffer frequency corresponding to different layers inside the ore body. During the adjustment of the ion liquid injection parameters, new ion liquid intensity data for ore pores were obtained.

[0010] A dynamic monitoring system for the leaching process of ion-adsorption rare earth ores includes: The leaching monitoring module includes a data monitoring unit and a leaching feature identification unit. The data monitoring unit is used to monitor and acquire real-time data on the ore pore ion liquid intensity and ore pore strain within the rare earth ore leaching area. Based on the ore pore ion liquid intensity and ore pore strain data, it constructs a mutual exclusion mapping function between leaching layers. It determines the interlayer expansion pressure data based on the interlayer mutual exclusion mapping function. By comparing the spatiotemporal evolution data of the interlayer expansion pressure data and the ore pore strain data, it obtains the stability parameter characteristics of the ore body. The leaching feature identification unit is used to scan the complex impedance spectrum distribution within the ore body to obtain the ore body pore phase angle data. Based on the ore body pore phase angle data, it identifies the ion exchange reaction state within the ore body to obtain the ore body ion exchange characteristics. Based on the ore body ion exchange characteristics, it judges the leaching process to obtain the ore body leaching state characteristics. The leaching adjustment module includes a leaching dynamic adjustment unit. The leaching dynamic adjustment unit is used to perform game evolution based on the stability parameter characteristics and leaching state characteristics of the ore body to obtain the leaching injection threshold matrix corresponding to the current ore body. The ion liquid injection parameters are adjusted according to the leaching injection threshold matrix to obtain new ore pore ion liquid intensity data.

[0011] The present invention has the following advantages: 1. This invention comprehensively analyzes multi-dimensional data such as ore pore ion liquid strength, ore pore strain, and ion exchange reaction to monitor the ore body status in real time. It also adjusts the leaching injection rate according to dynamic changes to optimize the leaching process, achieve higher leaching efficiency, and reach the best leaching effect in a shorter time. Through a dynamic game evolution model, it automatically adjusts the injection threshold based on the ore body stability parameters and leaching status characteristics to ensure efficient leaching while maintaining the stability of the ore body structure. This avoids damage to the ore body caused by excessive injection and effectively improves the long-term stability of the ore body.

[0012] 2. This invention accurately identifies the ion exchange reaction state within the ore body, evaluates the evolution gradient of ion exchange saturation in real time, and combines it with the three-dimensional spatial information of the ore body to optimize the leaching process more precisely, which helps to improve the overall leaching effect of the ore body and reduce resource waste. By intelligently adjusting and optimizing the injection speed, it reduces excessive injection and the use of chemical reagents, thereby significantly reducing environmental pollution and resource waste, and ensuring the safety and stability of the ore body. Attached Figure Description

[0013] Figure 1 This is a schematic diagram of the structure of a dynamic monitoring system for the ion-type rare earth ore leaching process used in an embodiment of the present invention. Detailed Implementation

[0014] To enable those skilled in the art to better understand the technical solutions of this invention, the technical solutions of this invention will be clearly and completely described below with reference to the accompanying drawings in the embodiments of this invention.

[0015] Example 1: A method for dynamic monitoring of the leaching process of ion-adsorption rare earth ores, comprising the following steps: Within the rare earth ore leaching area, real-time monitoring is conducted to acquire ore pore ion liquid strength data and ore pore strain data. Based on the ore pore ion liquid strength data and ore pore strain data, an interlayer mutual exclusion mapping function is constructed. The interlayer expansion pressure data is determined according to the interlayer mutual exclusion mapping function. By comparing the spatiotemporal evolution data of the interlayer expansion pressure data and the ore pore strain data, the stability parameter characteristics of the ore body are obtained. In the leaching area of ​​ion-adsorption rare earth minerals, multi-layered, multi-directional integrated sensor arrays of multiple physical quantities are pre-embedded inside the ore body to collect real-time data on the ore pore ion liquid intensity and ore body pore strain. The ore pore ion liquid intensity data refers to the dynamic changes in the concentration of the electrolyte solution in the ore body pores and its corresponding conductivity. It is usually collected using an embedded four-electrode conductivity sensor to reflect the penetration degree and chemical potential distribution of the leaching solution in the ore micropores. The ore pore strain data refers to the microscopic size changes or displacements of the ore body structure under the combined action of hydraulic pressure and chemical expansion. It is generally acquired using fiber optic strain sensors or high-precision flexible displacement sensors, which can capture the physical signals of mineral lattice expansion or structural slippage caused by ion exchange.

[0016] The specific steps for comparing the spatiotemporal evolution data of interlayer expansion pressure data and ore pore strain data include: Based on the ore pore ion liquid strength data, a double electric layer osmotic pressure function is constructed to represent the repulsive interaction between like-electric particles within the ore body; the effective stress of the ore body is obtained; based on the ore pore strain data and the effective stress of the ore body, a structural stiffness calibration coefficient for the ore body is constructed; the double electric layer osmotic pressure function is corrected according to the structural stiffness calibration coefficient of the ore body to obtain the mutual repulsion mapping function between leaching layers; By utilizing real-time acquired ore pore ion liquid intensity data and combining it with Graham's double-layer theory, a double-layer osmotic pressure function describing the interaction between the surface charge of mineral particles and pore liquid ions is constructed to simulate the repulsion phenomenon caused by changes in ion concentration. Simultaneously, the current effective stress of the ore body is obtained by monitoring the overlying pressure and pore water pressure, and the real-time measured ore pore strain data is correlated with this effective stress to calculate a structural stiffness calibration coefficient that reflects the ore skeleton's resistance to deformation under leaching conditions. The initial double-layer osmotic pressure function is nonlinearly corrected using this ore body structural stiffness calibration coefficient to balance chemical repulsion and physical skeleton constraints, thus synthesizing the final interlayer mutual repulsion mapping function.

[0017] When constructing the interlayer mutual repulsion mapping function, the interlayer mutual repulsion mapping function is a mathematical model specifically used to quantitatively describe the coupling relationship between the microscopic repulsive force between particles and the macroscopic ore body strain caused by chemical ion exchange during the leaching process of ion-adsorption rare earth minerals. From a physical perspective, since ion-adsorption rare earth minerals are mostly negatively charged clay minerals, when high-concentration ammonium salt plasma liquids enter the ore pores and displace rare earth ions, it will cause drastic changes in the thickness and potential of the double electric layer on the particle surface, thereby generating a strong electrochemical repulsive force between particles. The interlayer mutual repulsion mapping function transforms the microscopic charge repulsion effect, which is difficult to measure directly, into a calculable mechanical index, thus accurately explaining why significant expansion or slippage occurs inside the ore body when the injection pressure remains unchanged.

[0018] The ore pore ion liquid intensity data is input into the interlayer mutual repulsion mapping function for calculation to obtain the particle repulsion increment of the ore body; based on the particle repulsion increment of the ore body, the interlayer expansion pressure data is obtained; the interlayer expansion pressure data and ore pore strain data are subjected to trend correlation analysis in the same time domain coordinate system to obtain the internal strain components of the ore body. The measured values ​​of the ion liquid intensity in the ore pores were substituted into the interlayer mutual repulsion mapping function. The incremental repulsive force between mineral particles due to the overlap of charge layers was calculated through the function and further converted into interlayer expansion pressure data. The obtained interlayer expansion pressure data and the ore pore strain data fed back by the sensor were placed in the same time-domain coordinate system. Trend correlation analysis was performed by the sliding window cross-correlation algorithm to identify the part of the strain signal that is highly synchronized with the expansion pressure fluctuation, thereby accurately extracting the internal strain component induced by chemical expansion inside the ore body.

[0019] Construct a three-dimensional space inside the ore body; calculate the strain components inside the ore body within the three-dimensional space to obtain the local gradient difference value of the ore body; identify the shear slip potential zone of the ore body caused by non-uniform chemical expansion based on the local gradient difference value of the ore body to obtain the deformation region of the ore body; calculate the deformation region of the ore body to obtain the stability parameter characteristics of the ore body. Based on the geological survey data of the mine, a three-dimensional digital model of the ore body's interior was constructed, including the coordinates of all sensor points. The strain components inside the ore body were mapped to the corresponding data nodes in this three-dimensional space. By performing spatial vector operations on the strain components in the three-dimensional space, the local gradient difference values ​​of the ore body between different spatial points were calculated. Non-uniform chemical expansion regions with drastic changes in strain gradient and a high risk of shear fracture were identified, namely, the shear slip potential zone of the ore body. By statistically analyzing the spatial volume, growth rate, and deformation degree of the potential zone, the ore body stability parameter characteristics representing the overall safety margin of the ore body were calculated, providing a key mechanical decision-making basis for subsequent adjustment of injection intensity.

[0020] The complex impedance spectrum distribution inside the ore body is scanned to obtain the ore body pore phase angle data; the ion exchange reaction state inside the ore body is identified based on the ore body pore phase angle data to obtain the ore body ion exchange characteristics; the leaching process is judged based on the ore body ion exchange characteristics to obtain the ore body leaching state characteristics. In the initial stage of identifying the ion exchange reaction state inside the ore body, the complex impedance spectrum distribution inside the ore body is detected in real time using a high-frequency sweep frequency device, and the ore body pore phase angle data that changes dynamically over time is extracted from it. The role of the ore body pore phase angle data is that when ion-adsorption rare earth minerals undergo ion exchange reactions, the double-layer capacitance effect on the surface of mineral particles will change significantly, and the fluctuation of microscopic electrical properties will be reflected in the phase shift of the AC signal. The specific steps for identifying the ion exchange reaction state inside the ore body based on ore body pore phase angle data include: The corresponding phase angle peak and corresponding phase angle characteristic frequency are extracted from the pore phase angle data of the ore body; the real-time coordinates of the corresponding phase angle peak and corresponding phase angle characteristic frequency are mapped into the ion exchange standard feature space; the response state of ion exchange is determined in the ion exchange standard feature space to obtain the ion exchange phase state of the ore body; The peak height of the phase angle curve and its corresponding characteristic frequency points are accurately extracted from the complex impedance data stream. The extracted phase angle peak and characteristic frequency are converted into corresponding spatial coordinates in real time and mapped to a pre-constructed ion exchange standard characteristic space coordinate system. The ion exchange standard characteristic space is a benchmark model established based on a large number of indoor leaching experiments and historical engineering data. It contains characteristic regions representing different reaction maturity levels. By comparing the geometric map of the current monitoring point coordinates in this space, the ion exchange response state that the ore body is undergoing can be determined. For example, it can be determined whether the region is in the active period of large-scale ion exchange, the saturation period that is close to completion, or the latent period when the reagent has not yet arrived, thereby determining the specific ion exchange phase of the ore body.

[0021] The distribution of ion exchange phases in the ore body at different monitoring times was statistically analyzed to obtain the evolution gradient of ion exchange saturation corresponding to different layers within the ore body; the ion exchange characteristics of the ore body were generated by combining the three-dimensional space within the ore body and the evolution gradient of ion exchange saturation. In the process of statistically analyzing the ion exchange phase distribution of the ore body, the phase change data of monitoring points at different depths and horizontal orientations of the ore body under continuous time series are summarized, thereby capturing the nonlinear migration law of the reaction state from active to saturated. By analyzing the time difference of each monitoring layer reaching a specific exchange stage, the migration rate of the ion exchange front in the vertical and horizontal directions is calculated, and the evolution gradient of ion exchange saturation corresponding to different physical layers inside the ore body is derived. The evolution gradient of ion exchange saturation quantifies the spatial decay and increment relationship of chemical reaction, which not only reflects the chemical potential energy loss of the reagent between ore layers, but also clearly reveals whether there is a non-uniform reaction phenomenon caused by permeability differences in the ore body.

[0022] When generating ion exchange characteristics of ore bodies, the calculated ion exchange saturation evolution gradient is deeply fused with the pre-constructed three-dimensional spatial digital twin model inside the ore body, transforming the discrete gradient values ​​into a dynamic attribute field with continuous spatial properties. Through this coupling of multi-dimensional data, the diffusion profile of chemical reaction waves inside the ore body is displayed in a three-dimensional view, forming ion exchange characteristics of the ore body covering the entire field.

[0023] The specific steps for judging the leaching process based on the ion exchange characteristics of the ore body include: Spatial interpolation algorithms are used to connect the ion exchange characteristics of the ore body at different monitoring times to obtain the reaction envelope diagram of the leaching process; by comparing the spatial displacement vectors at different monitoring times in the reaction envelope diagram of the leaching process with the preset theoretical liquid flow propagation vector, the uniformity of leaching seepage is obtained; the leaching state is identified based on the uniformity of leaching seepage, and the leaching state characteristics of the ore body are obtained. In the process of constructing the reaction envelope diagram of the leaching process, a high-precision spatial interpolation algorithm is used to logically associate and spatially smooth the ion exchange characteristics of the ore body at different monitoring layers and different time nodes distributed within the ore body, thereby generating a continuous and dynamically changing equivalent reaction interface in the three-dimensional digital model. It should be noted that the preset theoretical fluid propulsion vector is a theoretical velocity and direction vector of the leaching fluid in the ore pores, calculated by a hydraulic model under ideal conditions based on the physical properties of the ore body and the injection process parameters. It represents the spatial trajectory reference value of the leaching fluid, driven by gravity and pressure, moving uniformly from the injection well to the extraction area under the assumption that the ore body is completely homogeneous and ideally distributed. The preset theoretical fluid propulsion vector not only includes the instantaneous velocity of the fluid movement, but also the direction of motion defined by the formation dip angle, permeability, and injection pressure gradient.

[0024] When assessing the uniformity of leaching flow, geometric analysis of the morphological evolution of the reaction envelope diagram over a continuous time series is performed to accurately extract the actual spatial displacement vectors of the reaction interfaces in each direction. These vectors are then compared in multiple dimensions with the theoretical flow propagation vectors under a pre-defined ideal Darcy velocity. By calculating the Euclidean distance and angular offset between the actual and theoretical displacement vectors, the flow deviations caused by ore heterogeneity or differences in pore structure are quantified. Based on the uniformity of leaching flow, combined with the expansion rate and saturation coverage of the reaction envelope diagram, the current leaching conditions are intelligently classified and identified. For example, if the monitored flow uniformity is extremely high and the reaction front advances smoothly, it is identified as a high-efficiency leaching state; if significant displacement vector offsets or local reaction lags are found, it is determined to be a short-circuit flow or leaching dead zone state. By integrating spatiotemporal evolution patterns, the ore body leaching state characteristics are generated, indicating the resource extraction efficiency and process adaptability of the current injection scheme at the chemical recovery level, providing direct decision support for subsequent execution of game-theoretic optimization strategies and adjustment of injection thresholds.

[0025] Based on the stability parameter characteristics and leaching state characteristics of the ore body, a game evolution is performed to obtain the leaching injection threshold matrix corresponding to the current ore body; the ion liquid injection parameters are adjusted according to the leaching injection threshold matrix to obtain new ore pore ion liquid intensity data. The specific steps for game evolution based on the stability parameters and leaching state characteristics of the ore body include: A game theory model for leaching evolution is constructed. The first participant is the leaching efficiency recovery entity, and the second participant is the leaching safety constraint entity. The leaching efficiency recovery entity tends to accelerate rare earth ion replacement by increasing injection concentration and pressure to shorten the production cycle and increase recovery rate. The leaching safety constraint entity aims to control the chemical expansion force and interlayer repulsion caused by ion exchange, preventing slippage or collapse of the ore body due to structural damage. These two participants counteract each other through the game theory model and seek a Nash equilibrium point, enabling the leaching process to overcome the limitations of a single objective and, while ensuring environmental and structural safety, seek an efficiency range for injection parameters.

[0026] The leaching safety penalty factor is calculated based on the stability parameter characteristics of the ore body; the leaching retention loss factor is calculated based on the retention loss component of the ore body leaching state characteristics; the leaching safety penalty factor and the leaching retention loss factor are input into the leaching evolution game model for evolution to obtain the leaching injection threshold matrix. The specific steps of evolution are as follows: A displacement-concentration adjustment operator is introduced. Within each evolution cycle, the fitness function is maximized by maximizing the uniformity of the three-dimensional space within the ore body. During the evolution cycle, the displacement-concentration adjustment operator is adjusted. When the change in the displacement-concentration adjustment operator causes the growth rate of the leaching safety penalty factor to exceed the growth rate of the leaching retention loss factor, a first adjustment strategy is implemented for the displacement-concentration adjustment operator. When the change in the displacement-concentration adjustment operator causes the growth rate of the leaching retention loss factor to exceed the growth rate of the leaching safety penalty factor, a second adjustment strategy is implemented for the displacement-concentration adjustment operator. Based on the strain increment and displacement trend in the stability parameters of the ore body, a leaching safety penalty factor is calculated to quantify the potential risks and costs of illegal operations. Simultaneously, combined with the acquired ore body leaching state characteristics, a leaching retention loss factor is defined by calculating the retention loss component. This factor represents the resource recovery loss caused by uneven injection, the existence of blind zones, or seepage short circuits. These two factors are input as core dynamic variables into an evolutionary game model to dynamically adjust their weights in the payoff function. This iterative calculation yields a set of equilibrium solutions that balance safety risks and leaching efficiency, resulting in the final leaching injection threshold matrix.

[0027] In the specific game evolution process, a displacement-concentration adjustment operator is introduced as the core control variable, and the global fitness function is to maximize the uniformity of the three-dimensional space inside the ore body in each evolution cycle. When the change of the displacement-concentration adjustment operator causes the growth rate of the leaching safety penalty factor to exceed the growth rate of the retention loss factor, it means that the pursuit of output is inducing unacceptable engineering risks. The first adjustment strategy will be implemented to prioritize the stability of the ore body by forcibly lowering the upper limit of concentration or reducing the injection pressure gradient. Conversely, if the growth rate of the leaching retention loss factor is dominant, it means that the current conservative strategy has led to serious resource waste. The second adjustment strategy will be implemented to moderately relax the injection parameter constraints to enhance the penetration of the seepage field. Through this dynamic game evolution, the corresponding leaching injection threshold is accurately output at each grid point in the three-dimensional space.

[0028] In actual leaching monitoring scenarios, the displacement-concentration adjustment operator is a composite control parameter used to coordinate the injection pump pressure (driving the liquid flow to generate displacement) and the initial concentration of the leaching agent (driving the ion exchange chemical reaction). It is used to establish a dynamic balance between the physical flow velocity and the chemical reaction rate. The displacement-concentration adjustment operator is essentially a coupling coefficient that specifies how the corresponding ion liquid concentration should compensate in the same or opposite direction when the injection pressure changes, causing the liquid flow displacement vector to change. This ensures that the chemical potential energy of the leaching agent can be adapted to the rare earth content in the pore structure in real time during the migration of the ore body. It should be noted that the first adjustment strategy, which implements the deceleration and efficiency enhancement logic by reducing the operator weight, specifically involves lowering the frequency of the injection pump to slow down fluid scouring and simultaneously finely adjusting the drug concentration to avoid excessive expansion of mineral particles induced by excessive chemical potential energy. Conversely, the second adjustment strategy implements the pressure boosting and penetration logic by increasing the operator weight. Specifically, it involves appropriately increasing the injection pressure within a safe threshold to enhance the long-range spatial directionality of the liquid flow and increasing the ionic liquid concentration to strengthen the penetration force in areas with poor permeability, thereby making the penetration more uniform across the entire field.

[0029] The specific steps for adjusting the ion liquid injection parameters based on the leaching injection threshold matrix include: The threshold matrix of leaching injection is analyzed to extract the ion liquid injection parameters. The ion liquid injection parameters include the ion liquid injection pressure boundary, chemical potential energy gradient coefficient and shear strain buffer frequency corresponding to different layers inside the ore body. During the adjustment of the ion liquid injection parameters, new ion liquid intensity data for ore pores were obtained.

[0030] In the process of analyzing the leaching injection threshold matrix, a set of ion liquid injection parameters that can guide the execution of production site operations is extracted using a multi-dimensional matrix inversion algorithm. These parameters are specifically refined into ion liquid injection pressure boundaries for different formation depths and horizontal layers within the ore body, used to mitigate the risk of ore body fracturing due to pressure overload. Simultaneously, the chemical potential energy gradient coefficient is extracted to define the step distribution requirements of the leaching agent concentration; and the shear strain buffer frequency serves as a reference indicator for adjusting the injection pump pulsation cycle to mitigate physical shear damage. During the adjustment of ion liquid injection parameters, the automated control system, based on the analyzed pressure boundaries and potential energy gradients, implements... By adjusting the output frequency of the pump station and the solute ratio of the liquid preparation system, the liquid injected into the ore pores is made to approach the optimal values ​​predicted by the game theory in terms of both physical flow rate and chemical characteristics. With the implementation of the new injection scheme, the diffusion and adsorption states of the liquid in the micro-pore structure of the ore undergo significant changes. These subtle changes are captured in real time by integrated sensors, and new ore pore ionic liquid intensity data are calculated by combining the solute transport model. This reflects the filling saturation of the liquid in the macro-pores in real time, and more accurately quantifies the instantaneous energy distribution of ion exchange reactions at the micro level. This provides basic data support for verifying the effectiveness of the game theory strategy and initiating the next round of evolution cycle.

[0031] In this embodiment, when monitoring and controlling the injection of ion-adsorption rare earth in a southern ore mine, a high-density multilayer potential sensor network is first deployed inside the ore pile to acquire resistivity and pore liquid titration data in real time. The system uses a spatial interpolation algorithm to connect the collected data and generate a three-dimensional leaching process reaction envelope diagram, which visually shows the irregular expansion of the liquid front towards the southeast corner of the ore body. By comparing the actual displacement vector extracted from the envelope diagram with the preset theoretical liquid flow propagation vector, the system identifies that there is an obvious seepage channel in the area, resulting in a decrease in uniformity, and determines it to be a state of dominant flow deviation. Subsequently, the system starts an evolutionary game model, introduces a displacement-concentration adjustment operator, and iteratively calculates the safety penalty factor (preventing slope instability) and the retention loss factor (reducing the loss of mineral resources). When the operator adjustment detects an increase in safety risk, the first adjustment strategy is triggered, which lowers the local injection pressure boundary and corrects the chemical potential gradient coefficient. Finally, the variable frequency pump parameters are adjusted according to the generated threshold matrix to make the ion liquid intensity data in the pores reach a new equilibrium, effectively suppressing the expansion of the blind zone.

[0032] Example 2: A dynamic monitoring system for the leaching process of ion-adsorption rare earth ore, see [link to example]. Figure 1 As shown, it includes: The leaching monitoring module includes a data monitoring unit and a leaching feature identification unit. The data monitoring unit is used to monitor and acquire real-time data on the ore pore ion liquid intensity and ore pore strain within the rare earth ore leaching area. Based on the ore pore ion liquid intensity and ore pore strain data, it constructs a mutual exclusion mapping function between leaching layers. It determines the interlayer expansion pressure data based on the interlayer mutual exclusion mapping function. By comparing the spatiotemporal evolution data of the interlayer expansion pressure data and the ore pore strain data, it obtains the stability parameter characteristics of the ore body. The leaching feature identification unit is used to scan the complex impedance spectrum distribution within the ore body to obtain the ore body pore phase angle data. Based on the ore body pore phase angle data, it identifies the ion exchange reaction state within the ore body to obtain the ore body ion exchange characteristics. Based on the ore body ion exchange characteristics, it judges the leaching process to obtain the ore body leaching state characteristics. The leaching adjustment module includes a leaching dynamic adjustment unit. The leaching dynamic adjustment unit is used to perform game evolution based on the stability parameter characteristics and leaching state characteristics of the ore body to obtain the leaching injection threshold matrix corresponding to the current ore body. The ion liquid injection parameters are adjusted according to the leaching injection threshold matrix to obtain new ore pore ion liquid intensity data.

[0033] It should be understood that those skilled in the art can make improvements or modifications based on the above description, and all such improvements and modifications should fall within the protection scope of the appended claims. Parts not described in detail in this specification are prior art known to those skilled in the art.

Claims

1. A method for dynamic monitoring of the leaching process of ion-adsorption rare earth ores, characterized in that, Includes the following steps: Within the rare earth ore leaching area, real-time monitoring is conducted to acquire ore pore ion liquid strength data and ore pore strain data. Based on the ore pore ion liquid strength data and ore pore strain data, an interlayer mutual exclusion mapping function is constructed. The interlayer expansion pressure data is determined according to the interlayer mutual exclusion mapping function. By comparing the spatiotemporal evolution data of the interlayer expansion pressure data and the ore pore strain data, the stability parameter characteristics of the ore body are obtained. The complex impedance spectrum distribution inside the ore body is scanned to obtain the ore body pore phase angle data; the ion exchange reaction state inside the ore body is identified based on the ore body pore phase angle data to obtain the ore body ion exchange characteristics; the leaching process is judged based on the ore body ion exchange characteristics to obtain the ore body leaching state characteristics. Based on the stability parameters and leaching state characteristics of the ore body, a game-theoretic evolution is performed to obtain the leaching injection threshold matrix corresponding to the current ore body. The ion liquid injection parameters are adjusted according to the leaching injection threshold matrix to obtain new ore pore ion liquid intensity data.

2. The method for dynamic monitoring of ion-adsorption rare earth ore leaching process according to claim 1, characterized in that, The specific steps for comparing the spatiotemporal evolution data of interlayer expansion pressure data and ore pore strain data include: Based on the ore pore ion liquid strength data, a double electric layer osmotic pressure function is constructed to represent the repulsive interaction between like-electric particles within the ore body; the effective stress of the ore body is obtained; based on the ore pore strain data and the effective stress of the ore body, a structural stiffness calibration coefficient for the ore body is constructed; the double electric layer osmotic pressure function is corrected according to the structural stiffness calibration coefficient of the ore body to obtain the mutual repulsion mapping function between leaching layers; The ore pore ion liquid intensity data is input into the interlayer mutual repulsion mapping function for calculation to obtain the particle repulsion increment of the ore body; based on the particle repulsion increment of the ore body, the interlayer expansion pressure data is obtained; the interlayer expansion pressure data and ore pore strain data are subjected to trend correlation analysis in the same time domain coordinate system to obtain the internal strain components of the ore body. A three-dimensional space inside the ore body is constructed; the strain components inside the ore body are calculated in the three-dimensional space to obtain the local gradient difference value of the ore body; based on the local gradient difference value of the ore body, the shear slip potential zone of the ore body caused by non-uniform chemical expansion is identified to obtain the deformation region of the ore body; the deformation region of the ore body is calculated to obtain the stability parameter characteristics of the ore body.

3. The method for dynamic monitoring of ion-adsorption rare earth ore leaching process according to claim 2, characterized in that, The specific steps for identifying the ion exchange reaction state inside the ore body based on ore body pore phase angle data include: Extract the corresponding phase angle peak and corresponding phase angle characteristic frequency from the ore body pore phase angle data; map the real-time coordinates of the corresponding phase angle peak and corresponding phase angle characteristic frequency into the ion exchange standard feature space; The response state of ion exchange is determined within the standard characteristic space of ion exchange to obtain the ion exchange phase state of the ore body; The distribution of ion exchange phases in the ore body at different monitoring times was statistically analyzed to obtain the evolution gradient of ion exchange saturation corresponding to different levels within the ore body; the ion exchange characteristics of the ore body were generated by combining the three-dimensional space within the ore body and the evolution gradient of ion exchange saturation.

4. The method for dynamic monitoring of ion-adsorption rare earth ore leaching process according to claim 3, characterized in that, The specific steps for judging the leaching process based on the ion exchange characteristics of the ore body include: Spatial interpolation algorithms are used to connect the ion exchange characteristics of the ore body at different monitoring times to obtain the reaction envelope diagram of the leaching process. By comparing the spatial displacement vectors at different monitoring times in the reaction envelope diagram of the leaching process with the preset theoretical liquid flow propagation vector, the uniformity of leaching seepage is obtained. The leaching state is identified based on the uniformity of leaching seepage, and the leaching state characteristics of the ore body are obtained.

5. The method for dynamic monitoring of ion-adsorption rare earth ore leaching process according to claim 4, characterized in that, The specific steps for game evolution based on the stability parameters and leaching state characteristics of the ore body include: A game theory model for leaching evolution is constructed, in which the first participant is the entity responsible for leaching efficiency recovery, and the second participant is the entity responsible for leaching safety constraints. The leaching safety penalty factor is calculated based on the stability parameter characteristics of the ore body; the leaching retention loss factor is calculated based on the retention loss component of the ore body leaching state characteristics; the leaching safety penalty factor and the leaching retention loss factor are input into the leaching evolution game model for evolution to obtain the leaching injection threshold matrix. The specific steps of evolution are as follows: A displacement-concentration adjustment operator is introduced. Within each evolution cycle, the fitness function is maximized by maximizing the uniformity of the three-dimensional space within the ore body. Within the evolution cycle, the displacement-concentration adjustment operator is adjusted. When the change in the displacement-concentration adjustment operator causes the growth rate of the leaching safety penalty factor to exceed the growth rate of the leaching retention loss factor, a first adjustment strategy is implemented for the displacement-concentration adjustment operator. When the change in the displacement-concentration adjustment operator causes the growth rate of the leaching retention loss factor to exceed the growth rate of the leaching safety penalty factor, a second adjustment strategy is implemented for the displacement-concentration adjustment operator.

6. The method for dynamic monitoring of ion-adsorption rare earth ore leaching process according to claim 5, characterized in that, The specific steps for adjusting the ion liquid injection parameters based on the leaching injection threshold matrix include: The threshold matrix of leaching injection is analyzed to extract the ion liquid injection parameters. The ion liquid injection parameters include the ion liquid injection pressure boundary, chemical potential energy gradient coefficient and shear strain buffer frequency corresponding to different layers inside the ore body. During the adjustment of the ion liquid injection parameters, new ion liquid intensity data for ore pores were obtained.

7. A dynamic monitoring system for the leaching process of ion-adsorption rare earth ores, characterized in that, The system employs a dynamic monitoring method for the leaching process of ion-adsorption rare earth ores as described in any one of claims 1-6, comprising: The leaching monitoring module includes a data monitoring unit and a leaching feature identification unit. The data monitoring unit is used to monitor and acquire real-time data on the ore pore ion liquid intensity and ore pore strain within the rare earth ore leaching area. Based on the ore pore ion liquid intensity and ore pore strain data, it constructs a mutual exclusion mapping function between leaching layers. It determines the interlayer expansion pressure data based on the interlayer mutual exclusion mapping function. By comparing the spatiotemporal evolution data of the interlayer expansion pressure data and the ore pore strain data, it obtains the stability parameter characteristics of the ore body. The leaching feature identification unit is used to scan the complex impedance spectrum distribution within the ore body to obtain the ore body pore phase angle data. Based on the ore body pore phase angle data, it identifies the ion exchange reaction state within the ore body to obtain the ore body ion exchange characteristics. Based on the ore body ion exchange characteristics, it judges the leaching process to obtain the ore body leaching state characteristics. The leaching adjustment module includes a leaching dynamic adjustment unit. The leaching dynamic adjustment unit is used to perform game evolution based on the stability parameter characteristics and leaching state characteristics of the ore body to obtain the leaching injection threshold matrix corresponding to the current ore body. The ion liquid injection parameters are adjusted according to the leaching injection threshold matrix to obtain new ore pore ion liquid intensity data.

Images