CN110983072B - Method for calculating injection time of ore leaching agent solution for in-situ ore leaching of ionic rare earth mine - Google Patents

Method for calculating injection time of ore leaching agent solution for in-situ ore leaching of ionic rare earth mine Download PDF

Info

Publication number
CN110983072B
CN110983072B CN201911186033.9A CN201911186033A CN110983072B CN 110983072 B CN110983072 B CN 110983072B CN 201911186033 A CN201911186033 A CN 201911186033A CN 110983072 B CN110983072 B CN 110983072B
Authority
CN
China
Prior art keywords
ore
coefficient
soil
mineral
liquid injection
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201911186033.9A
Other languages
Chinese (zh)
Other versions
CN110983072A (en
Inventor
王观石
洪本根
秦磊
彭陈亮
胡世丽
罗嗣海
张春雷
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Jiangxi University of Science and Technology
Original Assignee
Jiangxi University of Science and Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Jiangxi University of Science and Technology filed Critical Jiangxi University of Science and Technology
Priority to CN201911186033.9A priority Critical patent/CN110983072B/en
Publication of CN110983072A publication Critical patent/CN110983072A/en
Application granted granted Critical
Publication of CN110983072B publication Critical patent/CN110983072B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • C—CHEMISTRY; METALLURGY
    • C22—METALLURGY; FERROUS OR NON-FERROUS ALLOYS; TREATMENT OF ALLOYS OR NON-FERROUS METALS
    • C22B—PRODUCTION AND REFINING OF METALS; PRETREATMENT OF RAW MATERIALS
    • C22B59/00—Obtaining rare earth metals
    • C—CHEMISTRY; METALLURGY
    • C22—METALLURGY; FERROUS OR NON-FERROUS ALLOYS; TREATMENT OF ALLOYS OR NON-FERROUS METALS
    • C22B—PRODUCTION AND REFINING OF METALS; PRETREATMENT OF RAW MATERIALS
    • C22B3/00—Extraction of metal compounds from ores or concentrates by wet processes
    • C22B3/04—Extraction of metal compounds from ores or concentrates by wet processes by leaching
    • C—CHEMISTRY; METALLURGY
    • C22—METALLURGY; FERROUS OR NON-FERROUS ALLOYS; TREATMENT OF ALLOYS OR NON-FERROUS METALS
    • C22B—PRODUCTION AND REFINING OF METALS; PRETREATMENT OF RAW MATERIALS
    • C22B3/00—Extraction of metal compounds from ores or concentrates by wet processes
    • C22B3/04—Extraction of metal compounds from ores or concentrates by wet processes by leaching
    • C22B3/12—Extraction of metal compounds from ores or concentrates by wet processes by leaching in inorganic alkaline solutions
    • C22B3/14—Extraction of metal compounds from ores or concentrates by wet processes by leaching in inorganic alkaline solutions containing ammonia or ammonium salts
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
    • Y02P10/00—Technologies related to metal processing
    • Y02P10/20—Recycling

Abstract

The invention relates to a method for calculating the injection time of an ore leaching agent solution for in-situ ore leaching of an ionic rare earth mine, which is suitable for calculating the injection time of the ore leaching agent solution for each liquid injection zone during zone liquid injection. The invention provides a basis for reasonably determining the injection time of the mineral leaching agent solution of each injection subarea by systematically researching a calculation method of the dry density, the nonuniform coefficient, the curvature coefficient, the spatial variability of the effective grain diameter and the permeability coefficient of the mineral soil and based on the ion type rare earth mine subarea injection technology based on the volume of the mineral soil. On the basis of the known consumption of the liquid injection partition mineral leaching agent, the dynamic regulation and control of the injection time of each liquid injection partition mineral leaching agent solution can be carried out, the leaching rate is improved, the usage amount of the mineral leaching agent is accurately controlled, and the environmental pollution is reduced.

Description

Method for calculating injection time of ore leaching agent solution for in-situ ore leaching of ionic rare earth mine
Technical Field
The invention relates to a method for calculating the injection time of an ore leaching agent solution for in-situ ore leaching of an ionic rare earth mine, which is suitable for calculating the injection time of the ore leaching agent solution for each liquid injection subarea during subarea liquid injection.
Background
At present, the extraction of ionic rare earth resources adopts an in-situ ore leaching process, the liquid injection step of the process is carried out according to the three-first principle of 'top-bottom-first', 'thick-then-thin' and 'liquid-first-water-last', and the dosage of the ore leaching agent cannot be adjusted according to the spatial distribution rule of ore bodies, so that the problems in two aspects are often caused: firstly, the excessive using amount of the mineral leaching agent at the mountain top is easily caused, so that the ammonia nitrogen in the ore soil exceeds the standard; and secondly, the mineral leaching agent dosage at the injection liquid boundary is easy to be insufficient, so that the rare earth resource recovery is insufficient, and the resource loss is caused. Therefore, researchers have proposed an ionic rare earth partition injection method based on the volume of the ore.
The technology calculates the consumption of the liquid injection partition mineral leaching agent according to the cation exchange capacity of the unit volume of the mineral soil, and on the basis, the liquid injection time of the liquid injection partition mineral leaching agent solution is calculated according to the average permeability coefficient. Because the spatial variability of the permeability coefficient is not considered, the actual injection time of the mineral leaching agent solution of each injection partition cannot be accurately calculated, and the problem of insufficient mineral leaching agent consumption of partial areas still exists for ore blocks with obvious permeability coefficient variation. Therefore, how to determine the permeability coefficient of each injection zone is a key issue in ensuring that a sufficient amount of the mineral leaching agent solution is injected.
Disclosure of Invention
The invention aims to provide a method for calculating the injection time of an ore leaching agent solution for in-situ ore leaching of an ionic rare earth mine, and the method can be used for achieving the purpose of optimizing the injection time of each liquid injection partition ore leaching agent solution during in-situ ore leaching.
The technical scheme of the invention is as follows: a method for calculating the injection time of an ore leaching agent solution for in-situ ore leaching of an ionic rare earth mine comprises the following steps:
step one, calculating the permeability coefficient of the ore soil;
selecting a test point of the permeability coefficient of the mineral soil in the mineral block planned to be mined, recording the coordinate position of the test point, and measuring the permeability coefficient of the mineral soil at the position of the test point where the mineral is found to be deep, the dry density of the mineral soil and the particle size distribution characterization parameters of the mineral soil by adopting the prior art, wherein the characterization parameters comprise a non-uniform coefficient, a curvature coefficient and an effective particle size; substituting the permeability coefficient, the dry density, the nonuniform coefficient, the curvature coefficient and the effective grain diameter of the tested ore soil into a relational expression (1) to determine a fitting parameter A;
relation (1):
Figure BDA0002292419040000021
in relation (1): k is the permeability coefficient; rhodDry density of the mineral soil'; cuIs a non-uniformity coefficient; ccIs a curvature coefficient; d10Is an effective granuleDiameter; a is a fitting parameter;
step two, testing the dry density and the particle size distribution of the ore soil in the ore block;
arranging sampling points according to the sampling mesh, recording the coordinate positions of the sampling points, completely taking out all the mineral soil of all the sampling points within the mineral depth range, testing the dry density and the grain composition of the mineral soil, and selecting characterization parameters of the grain composition, such as uneven coefficients, curvature coefficients and effective grain diameters;
analyzing the spatial variability of the dry density and the particle grading of the ore soil in the ore block;
determining a variation function value of each parameter by a geostatistical method (the prior art) according to test data of the dry density, the nonuniform coefficient, the curvature coefficient and the effective grain size of the ore soil at all sampling points in the ore block, fitting a spherical variation function model, an exponential variation function model and a Gaussian variation function model according to the variation function values, and determining an optimal variation function model and a variation range of the dry density, the nonuniform coefficient, the curvature coefficient and the effective grain size of the ore soil in the ore block;
step four, calculating the weight of the dry density of the ore soil and the corresponding parameter of the particle grading characterization parameter to the target point;
selecting the maximum variation of the four parameters of the dry density, the nonuniform coefficient, the curvature coefficient and the effective grain size of the ore soil in the ore block determined in the third step as a calculation unit, setting each liquid injection hole for in-situ ore leaching as a target point, and calculating according to a relational expression (2) to obtain the weight of the dry density, the nonuniform coefficient, the curvature coefficient and the effective grain size of the ore soil at each sampling point in the calculation unit on the corresponding parameters of the target point;
relation (2):
Figure BDA0002292419040000031
in relation (2): gamma rayijThe function values of the variation from the ith to the jth (i, j is 1 to N); n is the number of sampling points used to estimate the target point; w is akIs a sampling point xkFor target point xpThe weight of (k) is 1 to N; gamma rayNpThe function value of the variation from the target point to any sampling point; subscript p is the target point number; phi is the Lagrangian multiplier;
step five, calculating the dry density and particle size distribution characterization parameters of the ore soil at each liquid injection hole;
according to the dry density, the nonuniform coefficient, the curvature coefficient and the weight of the effective grain diameter of the ore soil at all sampling points in the calculation unit in the fourth step to the corresponding parameters at the liquid injection hole, the dry density, the nonuniform coefficient, the curvature coefficient and the effective grain diameter of the ore soil at the liquid injection hole are calculated according to the relation (3):
relation (3):
Figure BDA0002292419040000032
in relation (3): f (x)p) The dry density, uneven coefficient, curvature coefficient or effective grain size of the ore soil at the liquid injection hole; f (x)k) Is a sampling point xk(k-1, 2, …, N) dry density, coefficient of non-uniformity, coefficient of curvature, or effective particle size; w is akIs a sampling point xkFor target point xpThe weight of (k) is 1 to N; n is the number of sampling points; subscript p is the target point number;
step six, calculating the permeability coefficient of the ore soil at each liquid injection hole;
calculating the permeability coefficient of the ore soil at each liquid injection hole according to the relational expression (1) according to the calculation method of the permeability coefficient of the ore soil in the step one and the dry density, the uneven coefficient, the curvature coefficient and the effective particle size of the ore soil at each liquid injection hole calculated in the step five;
step seven, calculating the injection duration of the mineral leaching agent solution in each injection hole;
under the premise of determining the permeability coefficient of the ore soil at each injection hole, calculating the single-hole injection strength under the mesh injection condition according to the relational expression (4); on the basis, according to a published patent ionic rare earth partition liquid injection method based on the ore volume (application number 2019105596678), the project predicted leaching rate is selected, the consumption of the leaching agent per unit volume of the ore soil is determined, and the liquid injection duration of each liquid injection hole is calculated according to a relational expression (5);
relation (4):
Figure BDA0002292419040000041
in relation (4): q. q ofmThe single-hole injection strength is the unit m under the condition of hole-mesh injection3/d；laThe hole distance of the liquid injection holes is expressed in m; alpha is a fitting parameter and takes a value of 1.61; k is a permeability coefficient and has a unit of m/d; r is0The radius of the liquid injection hole is expressed in m; h is the height of the liquid level in the liquid injection hole and is unit m.
Relation (5):
Figure BDA0002292419040000042
in relation (5): t is the injection duration of the mineral leaching agent solution, unit d; laThe hole distance of the liquid injection holes is expressed in m; l. thebThe row spacing of the liquid injection holes is m; h is the thickness of the ore soil at the liquid injection hole, and the unit is m; beta is the consumption of the mineral leaching agent per unit volume of the mineral soil, and the unit is kg/m3(ii) a c is the concentration of the mineral leaching agent solution, and the unit is kg/m3；qm(K) The single-hole injection strength is the unit m under the condition of hole-mesh injection3/d。
In the first step, 2-5 test points of permeability coefficient of the mineral soil are selected, and the position of the test points is observed to have the mineral depth of 0.5-3 m.
In the second step, sampling points are arranged according to a sampling mesh of 2-30 m multiplied by 2-30 m, and the mining depth of all the sampling points is 0.5-3 m.
In the third step, the optimal variation function models of the dry density, the nonuniform coefficient, the curvature coefficient and the effective grain diameter of the ore soil in the ore block are all exponential models, the function expression of the model is shown as the relational expression (6),
relation (6):
Figure BDA0002292419040000051
relation formula(6) The method comprises the following steps: gamma (h) is a variogram; c0Is the gold lump constant; c is the arch height; a is a variable range and is in a unit of m; h is the distance between two sample points in m.
Therefore, the permeability coefficient of the ore soil is determined according to characterization parameters (such as nonuniform coefficients, curvature coefficients and effective particle diameters) of the dry density and the particle composition of the ore soil, and the permeability coefficient of each injection partition is further determined on the basis of analyzing the spatial variability of the characterization parameters of the dry density and the particle composition of the ore soil. And then, calculating the injection duration of the mineral leaching agent solution of each injection partition by taking the volume of the ore soil of each injection partition as a reference according to the permeability coefficient of each injection partition.
The invention provides a method for calculating the injection time of the leaching agent solution for in-situ leaching of the ionic rare earth mine by systematically researching the calculation methods of the dry density, the nonuniform coefficient, the curvature coefficient, the spatial variability of the effective particle size and the permeability coefficient of the ore soil and based on the ionic rare earth mine partitioned liquid injection technology based on the volume of the ore soil, and provides a basis for reasonably determining the injection time of the leaching agent solution of each liquid injection partition. By using the calculation method, the injection duration of the mineral leaching agent solution of each injection partition can be dynamically regulated and controlled on the basis of the known consumption of the mineral leaching agent of the injection partition, the leaching rate is improved (the embodiment is improved by 6.22%), the usage amount of the mineral leaching agent is accurately controlled, and the environmental pollution is reduced.
Detailed Description
In consideration of the problems of large test workload, high test cost and the like existing in the actual measurement of permeability coefficients of all the subareas, the invention provides a method for calculating the injection duration of the mineral leaching agent solution for in-situ mineral leaching of an ionic rare earth mine by systematically researching a calculation method of dry density, an uneven coefficient, a curvature coefficient, spatial variability of effective particle size and a permeability coefficient of mineral soil and based on an ionic rare earth mine subarea injection technology taking the volume of the mineral soil as a basis.
The invention is applied to the undisclosed test, calculates the injection duration of the mineral leaching agent solution of each injection subarea when injecting the liquid into a certain full-coverage rare earth mine subarea in Fujian, and is specifically described as follows:
step one, calculating the permeability coefficient of ore soil;
2 test points were selected in the block planned for mining, with X and Y coordinates [179943.525, 202036.03] respectively]And [180016.525, 202036.03]Excavating to 1m of the ore bed, and testing the permeability coefficients of the ore soil of the two test points by adopting a single-ring method, wherein the test values are 0.88m/d and 0.58m/d respectively; the ore soil at the depth of 1m to 1.5m of the ore is completely taken out by a portable QTZ-8 soil-taking drilling machine, and the dry density, the nonuniform coefficient, the curvature coefficient and the effective grain diameter of the ore soil are tested by adopting an indoor experimental method, wherein coordinate points [179943.525, 202036.03]]The dry density, the nonuniform coefficient, the curvature coefficient and the effective grain diameter of the ore soil are respectively 1.41g/cm318.91, 1.46 and 0.024mm, coordinate points [180016.525, 202036.03]The dry density, the nonuniform coefficient, the curvature coefficient and the effective grain diameter of the mineral soil are respectively 1.49g/cm318.64, 1.22 and 0.025 mm; thus, the fitting parameter a was determined to be 0.057 according to relation (1);
relation (1):
Figure BDA0002292419040000061
in relation (1): k is the permeability coefficient; ρ is a unit of a gradientdDry density of the ore soil; cuIs a non-uniformity coefficient; ccIs the curvature coefficient; d10Is an effective particle size; a is a fitting parameter;
step two, testing the dry density and the particle size distribution of the ore soil in the ore block;
the area of the unit is set to be 8m multiplied by 8m, the ore block is divided into a plurality of units, and the invention is explained by taking a ridge line as an example. Numbering units corresponding to the ridge line from No. 1 to No. 9, arranging sampling points in each unit, and acquiring actual coordinate values of the sampling points, wherein partial data are shown in an attached table 1; taking out all the mineral soil in the depth range of 0.5-3 m from each sampling point by using a portable QTZ-8 soil-taking drilling machine, testing the dry density and the particle gradation of the mineral soil, wherein the characterization parameters of the particle gradation of the mineral soil comprise a non-uniform coefficient, a curvature coefficient and an effective particle size, and part of data are shown in an attached table 1;
thirdly, analyzing the space variability of the dry density and the particle grading of the ore soil in the ore block,
determining the variation function values of all parameters by a geostatistical method according to the test data of the dry density, the nonuniform coefficient, the curvature coefficient and the effective grain size of the ore soil at all sampling points in the ore block, fitting a spherical variation function model, an exponential variation function model and a Gaussian variation function model according to the variation function values, determining that the optimal variation function models of the dry density, the nonuniform coefficient, the curvature coefficient and the effective grain size of the ore soil in the ore block are all exponential models, and determining that the variation ranges of the dry density, the nonuniform coefficient, the curvature coefficient and the effective grain size of the ore soil in the ore block are respectively 28m, 21m, 35m and 34m by using a relational expression (6) as shown in the relational expression (6);
relation (6):
Figure BDA0002292419040000071
in relation (6): gamma (h) is a variation function; c0Is the lump constant; c is the arch height; a is a variable range and is in a unit of m; h is the distance between two sampling points in m;
step four, calculating the weight of the corresponding parameters of the target point by the dry density and the particle grading characterization parameters of the mineral soil,
taking a certain liquid injection hole on the ridge line as a target point, wherein X and Y coordinates of the liquid injection hole are [179987.525, 202036.03], determining a calculation unit according to the variation of the curvature coefficient of the mineral soil in the mineral block in the step three, calculating to obtain the weight of the dry density, the non-uniform coefficient, the curvature coefficient and the effective grain diameter of the mineral soil at each sampling point to corresponding parameters of the liquid injection hole, wherein the calculation formula is shown as a relational expression (2), and part of data of the calculation result of the dry density of the mineral soil at the target point is shown in an attached table 2;
relation (2):
Figure BDA0002292419040000081
in relation (2): gamma rayijIs a firstThe function value of variation from the point i to the point j (i, j is 1 to N); n is the number of sampling points used to estimate the target point; w is akIs a sampling point xkFor target point xpThe weight of (k) is 1 to N; gamma rayNpThe function value of the variation from the target point to any sampling point; subscript p is the target point number;
Figure BDA0002292419040000083
is a Lagrangian multiplier;
step five, calculating the dry density of the ore soil and the particle grading characterization parameters at each injection hole,
according to the dry density, the nonuniform coefficient, the curvature coefficient and the weight of the effective grain diameter to the corresponding parameters of the liquid injection hole of the ore soil at each sampling point in the calculation unit selected in the step four, the dry density of the ore soil at the liquid injection hole is calculated to be 1.42g/cm3The coefficient of non-uniformity was 18.71, the coefficient of curvature was 1.23 and the effective particle diameter was 0.026mm, as calculated by the equation (3),
relation (3):
Figure BDA0002292419040000082
in relation (3): f (x)p) The dry density, the non-uniformity coefficient, the curvature coefficient or the effective grain diameter of a target point; f (x)k) Is a sampling point xk(k ═ 1,2, …, N) dry density, coefficient of non-uniformity, coefficient of curvature, or effective particle size; w is akIs a sampling point xkFor target point xpThe weight of (k) is 1 to N; n is the number of sampling points; subscript p is the target point number;
step six, calculating the permeability coefficient of the ore soil at each liquid injection hole,
according to the calculation method of the permeability coefficient of the ore soil in the first step and the dry density, the uneven coefficient, the curvature coefficient and the effective grain diameter of the ore soil at the liquid injection hole calculated in the fifth step, the permeability coefficient of the ore soil at the liquid injection hole is calculated to be 0.65 m/d;
step seven, calculating the injection time of the mineral leaching agent solution in each injection hole,
ore blockThe hole distance of the liquid injection holes is 2m, the row distance of the liquid injection holes is 2m, the radius of the liquid injection holes is 0.09m, the liquid level height in the liquid injection holes is 0.6m, and the single-hole liquid injection strength under the hole-mesh liquid injection condition is calculated to be 1.15m according to the relational expression (4)3D; an estimated leaching rate of 90% is selected, according to a patent titled "ion type rare earth partition liquid injection method based on ore volume", and the consumption of the leaching agent per unit ore volume is 7kg/m under the condition of the leaching rate3(ii) a The concentration of the injected mineral leaching agent solution is 20kg/m3If the thickness of the ore soil at the liquid injection hole is 8m, calculating the liquid injection time of the mineral leaching agent solution at the liquid injection hole to be 9.74d according to the relational expression (5);
relation (4):
Figure BDA0002292419040000091
in relation (4): q. q ofmThe single-hole injection strength is the unit m under the condition of hole-mesh injection3/d；laThe hole distance of the liquid injection holes is expressed in m; alpha is a fitting parameter and takes a value of 1.61; k is the permeability coefficient and the unit m/d; r is0The radius of the liquid injection hole is expressed in m; h is the height of the liquid level in the liquid injection hole and is unit m.
Relation (5):
Figure BDA0002292419040000092
in relation (5): t is the injection duration of the mineral leaching agent solution, unit d; laThe hole distance of the liquid injection hole is unit m; l. thebThe row distance of the liquid injection holes is m; h is the thickness of the ore soil at the liquid injection hole, and the unit is m; beta is the consumption of the mineral leaching agent per unit volume of the mineral soil, and the unit is kg/m3(ii) a c is the concentration of the mineral leaching agent solution, and the unit is kg/m3；qm(K) The single-hole injection strength is the unit m under the condition of hole-mesh injection3/d。
The experimental effect is as follows: performing an undisclosed test on a certain full-coverage rare earth ore block in Fujian province to determine the ore leaching agent consumption of 7kg/m per unit ore soil volume3Concentration of ammonium sulfate solution during injectionThe ore block leaching rate is calculated to be 94.46% by adopting the method of the invention when the injection time of the ore block leaching agent solution is 20g/L, and the leaching rate is improved by 6.22% compared with the prior art under the condition that the unit consumption of the ore leaching agent is the same.
Attached Table 1:
function of variation X coordinate Y coordinate Dry density Coefficient of non-uniformity Coefficient of curvature Effective particle size
1 179943.525 202036.03 1.41g/cm3 18.91 1.46 0.024mm
2 179951.525 202036.03 1.51g/cm3 18.60 1.25 0.025mm
3 179967.525 202036.03 1.53g/cm3 18.41 1.02 0.030mm
4 179975.525 202036.03 1.28g/cm3 19.10 1.45 0.025mm
5 179983.525 202036.03 1.30g/cm3 19.25 1.68 0.021mm
6 179991.525 202036.03 1.54g/cm3 18.16 0.79 0.032mm
7 179999.525 202036.03 1.53g/cm3 18.35 1.00 0.027mm
8 180007.525 202036.03 1.29g/cm3 19.05 1.22 0.026mm
9 180016.525 202036.03 1.49g/cm3 18.64 1.22 0.025mm
。
Attached table 2: dry density sampling point to target point distribution coefficient calculation table
Figure BDA0002292419040000101
。

Claims (4)

1. A method for calculating the injection time of an ore leaching agent solution for in-situ ore leaching of an ionic rare earth mine is characterized by comprising the following steps of:
step one, calculating the permeability coefficient of ore soil;
selecting a test point of the permeability coefficient of the ore soil in the ore block planned to be mined, recording the coordinate position of the test point, and measuring characterization parameters of the permeability coefficient of the ore soil, the dry density of the ore soil and the grain composition of the ore soil at the position of the test point at the ore depth by adopting the prior art, wherein the characterization parameters comprise a non-uniform coefficient, a curvature coefficient and an effective grain diameter; substituting the permeability coefficient, the dry density, the nonuniform coefficient, the curvature coefficient and the effective grain diameter of the tested ore soil into a relational expression (1) to determine a fitting parameter A;
relation (1):
Figure FDA0002292419030000011
in relation (1): k is the permeability coefficient; ρ is a unit of a gradientdIs the dry density of the ore`；CuIs a non-uniformity coefficient; ccIs the curvature coefficient; d10Is an effective particle size; a is a fitting parameter;
step two, testing the dry density and the particle size distribution of the ore soil in the ore block;
arranging sampling points according to the sampling mesh degree, recording the coordinate positions of the sampling points, completely taking out all the mineral soil in the mineral depth range of all the sampling points, testing the dry density and the particle grading of the mineral soil, and selecting characterization parameters of the particle grading, wherein the characterization parameters are a non-uniform coefficient, a curvature coefficient and an effective particle size;
analyzing the spatial variability of the dry density and the particle grading of the ore soil in the ore block;
according to the test data of the dry density, the nonuniform coefficient, the curvature coefficient and the effective grain size of the ore soil at all sampling points in the ore block, determining the variation function value of each parameter by a geostatistical method, fitting a spherical variation function model, an exponential variation function model and a Gaussian variation function model according to the variation function value, and determining the optimal variation function model and the variation range of the dry density, the nonuniform coefficient, the curvature coefficient and the effective grain size of the ore soil in the ore block;
calculating the weight of the dry density of the mineral soil and the corresponding parameter weight of the particle grading characterization parameter to the target point;
selecting the maximum variation of the four parameters of the dry density, the nonuniform coefficient, the curvature coefficient and the effective particle size of the mineral soil in the ore block determined in the third step as a calculation unit, setting each liquid injection hole for in-situ ore leaching as a target point, and calculating according to a relational expression (2) to obtain the weight of the dry density, the nonuniform coefficient, the curvature coefficient and the effective particle size of the mineral soil at each sampling point in the calculation unit to the corresponding parameters of the target point;
relation (2):
Figure FDA0002292419030000021
in relation (2): gamma rayijThe function value of the variation from the ith point to the jth point, i and j are 1 to N; n is the number of sampling points used to estimate the target point; w is akIs a sampling point xkTo target point xpK is 1 to N; gamma rayNpThe function value of the variation from the target point to any sampling point; subscript p is the target point number; φ is the Lagrangian multiplier;
step five, calculating the dry density and particle size distribution characterization parameters of the ore soil at each liquid injection hole;
according to the dry density, the nonuniform coefficient, the curvature coefficient and the weight of the effective grain size of the ore soil at all the sampling points in the unit to the corresponding parameters at the liquid injection hole, the dry density, the nonuniform coefficient, the curvature coefficient and the effective grain size of the ore soil at the liquid injection hole are calculated according to the relation (3):
relation (3):
Figure FDA0002292419030000022
in relation (3): f (x)p) The dry density, uneven coefficient, curvature coefficient or effective grain size of the ore soil at the liquid injection hole; f (x)k) Is a sampling point xkDry density, non-uniformity coefficient, curvature coefficient or effective particle size of; w is akIs a sampling point xkTo target point xpThe weight of (c); n is the number of sampling points; subscript p is the target Point numberSubscript k ═ 1,2, …, N;
step six, calculating the permeability coefficient of the ore soil at each liquid injection hole;
calculating the permeability coefficient of the ore soil at each liquid injection hole according to the relational expression (1) according to the calculation method of the permeability coefficient of the ore soil in the step one and the dry density, the uneven coefficient, the curvature coefficient and the effective particle size of the ore soil at each liquid injection hole calculated in the step five;
step seven, calculating the injection duration of the mineral leaching agent solution in each injection hole;
under the premise of determining the permeability coefficient of the ore soil at each injection hole, calculating the single-hole injection strength under the mesh injection condition according to the relational expression (4); on the basis, according to a published patent ionic rare earth partitioned liquid injection method based on the ore volume, project predicted leaching rate is selected, the ore leaching agent consumption per unit ore soil volume is determined, and the liquid injection duration of each liquid injection hole is calculated according to a relational expression (5);
relation (4):
Figure FDA0002292419030000031
in relation (4): q. q.smThe single-hole injection strength is the unit m under the condition of hole-mesh injection3/d；laThe hole distance of the liquid injection holes is expressed in m; alpha is a fitting parameter and takes a value of 1.61; k is a permeability coefficient and has a unit of m/d; r is0The radius of the liquid injection hole is expressed in m; h is the height of the liquid level in the liquid injection hole and is unit m;
relation (5):
Figure FDA0002292419030000032
in relation (5): t is the injection duration of the mineral leaching agent solution, unit d; laThe hole distance of the liquid injection hole is unit m; l. thebThe row spacing of the liquid injection holes is m; h is the thickness of the ore soil at the liquid injection hole, and the unit is m; beta is the consumption of the mineral leaching agent per unit volume of the mineral soil, and the unit is kg/m3(ii) a c is the concentration of the mineral leaching agent solution in unitskg/m3；qm(K) The single-hole injection strength is the unit m under the condition of hole-mesh injection3/d。
2. The method for calculating the injection time of the mineral leaching agent solution for in-situ mineral leaching of the ionic rare earth mine according to claim 1, wherein the method comprises the following steps: in the first step, 2-5 test points of permeability coefficient of the mineral soil are selected, and the position of the test points is observed to have the mineral depth of 0.5-3 m.
3. The method for calculating the injection time of the mineral leaching agent solution for in-situ mineral leaching of the ionic rare earth mine according to claim 1, wherein the method comprises the following steps: and in the second step, sampling points are arranged according to a sampling mesh of 2-30 m multiplied by 2-30 m, and the depth of all the sampling points to the mine is 0.5-3 m.
4. The method for calculating the injection time of the mineral leaching agent solution for in-situ mineral leaching of the ionic rare earth mine according to claim 1, wherein the method comprises the following steps: in the third step, the optimal variation function models of the dry density, the nonuniform coefficient, the curvature coefficient and the effective grain diameter of the ore soil in the ore block are all exponential models, the function expression of the models is shown as the relational expression (6),
relation (6):
Figure FDA0002292419030000041
in relation (6): gamma (h) is a variogram; c0Is the lump constant; c is the arch height; a is a variable range and is in a unit of m; h is the distance between two sample points in m.
CN201911186033.9A 2019-11-28 2019-11-28 Method for calculating injection time of ore leaching agent solution for in-situ ore leaching of ionic rare earth mine Active CN110983072B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911186033.9A CN110983072B (en) 2019-11-28 2019-11-28 Method for calculating injection time of ore leaching agent solution for in-situ ore leaching of ionic rare earth mine

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911186033.9A CN110983072B (en) 2019-11-28 2019-11-28 Method for calculating injection time of ore leaching agent solution for in-situ ore leaching of ionic rare earth mine

Publications (2)

Publication Number Publication Date
CN110983072A CN110983072A (en) 2020-04-10
CN110983072B true CN110983072B (en) 2022-07-22

Family

ID=70087554

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911186033.9A Active CN110983072B (en) 2019-11-28 2019-11-28 Method for calculating injection time of ore leaching agent solution for in-situ ore leaching of ionic rare earth mine

Country Status (1)

Country Link
CN (1) CN110983072B (en)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102787236A (en) * 2012-08-21 2012-11-21 江西稀有金属钨业控股集团有限公司 Ion-absorbing type rare earth in-situ mineral leaching mother liquor recovery process
CN106939374A (en) * 2017-05-24 2017-07-11 江西理工大学 In-situ ionic rare earth soaks the computational methods of ore deposit ammonium sulfate fluid injection mass concentration
US10030286B1 (en) * 2013-11-13 2018-07-24 Ii-Vi Incorporated Method of direct solvent extraction of rare earth metals from an aqueous acid-leached ore slurry
CN110195158A (en) * 2019-06-26 2019-09-03 江西理工大学 Using ore volume as the ion type rareearth subregion electrolyte filling method of foundation

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103014359B (en) * 2012-11-27 2014-03-19 益阳鸿源稀土有限责任公司 Separation and recovery method of monazite slag

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102787236A (en) * 2012-08-21 2012-11-21 江西稀有金属钨业控股集团有限公司 Ion-absorbing type rare earth in-situ mineral leaching mother liquor recovery process
US10030286B1 (en) * 2013-11-13 2018-07-24 Ii-Vi Incorporated Method of direct solvent extraction of rare earth metals from an aqueous acid-leached ore slurry
CN106939374A (en) * 2017-05-24 2017-07-11 江西理工大学 In-situ ionic rare earth soaks the computational methods of ore deposit ammonium sulfate fluid injection mass concentration
CN110195158A (en) * 2019-06-26 2019-09-03 江西理工大学 Using ore volume as the ion type rareearth subregion electrolyte filling method of foundation

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Simulation of one-dimensional column leaching of;龙平;《Transactions of Nonferrous Metals Society of China》;20190315;625-633 *
渗透系数空间变异条件下临界注液强度的计算;洪本根;《稀土》;20180515;25-33 *

Also Published As

Publication number Publication date
CN110983072A (en) 2020-04-10

Similar Documents

Publication Publication Date Title
CN109685378B (en) TBM construction surrounding rock digchability grading method based on data mining
US7869954B2 (en) Using microseismic data to characterize hydraulic fractures
CN111690831B (en) Liquid injection process optimization method of ionic rare earth ore
EP2795527B1 (en) Systems and methods for selecting facies model realizations
CN105388513B (en) The method for building up and device of earthquake-capturing observation system
CN108731561B (en) Automatic calculation method for blast hole charging amount
WO2022078516A1 (en) Tunnel resistivity advanced detection optimization method and system based on depth resolution
CN103984788A (en) Automatic intelligent design and optimization system for anchor bolt support of coal tunnel
CN107366534A (en) It is roughened the determination method and apparatus of permeability
CN105093331B (en) The method for obtaining Rock Matrix bulk modulus
Samani et al. MODFLOW equipped with a new method for the accurate simulation of axisymmetric flow
CN103065051A (en) Method for performing grading and sectionalizing on rock mass automatically
CN110983072B (en) Method for calculating injection time of ore leaching agent solution for in-situ ore leaching of ionic rare earth mine
CN104181593A (en) Three-dimensional ray-free tracing contraflexure wave tomographic imaging method and device
CN116341201A (en) Survey and drawing execution supervision system based on geological survey engineering
CN108227033B (en) Method and device for predicting porosity of sandstone reservoir and computer storage medium
EP2678718A2 (en) Computerized method for the estimation of a value for at least a parameter of a hydrocarbon-producing region, for planning the operation and operating the region
CN112948924A (en) Near unconsolidated formation mining water flowing fractured zone height determination method based on overlying strata structure
CN110863814B (en) Method and device for determining single-well subsection specific fluid production index of ultra-thick type cuttings limestone oil reservoir
CN105137041A (en) Soil parameter spatial distribution monitoring method and system
CN110195158B (en) Ion type rare earth partition liquid injection method based on ore volume
CN114236624B (en) Method and system for estimating fracturing modification space volume based on electromagnetic method
CN115688307A (en) Tunnel face instability risk assessment method considering soil body strength uncertainty
CN110157905B (en) Ion type rare earth ore partition liquid injection method based on resource reserves
CN107767282A (en) A kind of evaluation method of uranium resource distribution

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant