CN116306173A - Optimization method for seepage calculation of tailing pond based on upstream boundary conditions - Google Patents
Optimization method for seepage calculation of tailing pond based on upstream boundary conditions Download PDFInfo
- Publication number
- CN116306173A CN116306173A CN202310536901.1A CN202310536901A CN116306173A CN 116306173 A CN116306173 A CN 116306173A CN 202310536901 A CN202310536901 A CN 202310536901A CN 116306173 A CN116306173 A CN 116306173A
- Authority
- CN
- China
- Prior art keywords
- seepage
- pond
- calculation
- boundary condition
- tailing pond
- 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.)
- Granted
Links
- 238000004364 calculation method Methods 0.000 title claims abstract description 105
- 238000011144 upstream manufacturing Methods 0.000 title claims abstract description 80
- 238000000034 method Methods 0.000 title claims abstract description 41
- 238000005457 optimization Methods 0.000 title claims abstract description 17
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 122
- 238000009826 distribution Methods 0.000 claims abstract description 18
- 238000013461 design Methods 0.000 claims description 14
- 238000012876 topography Methods 0.000 claims description 6
- 238000013519 translation Methods 0.000 claims description 6
- 230000035699 permeability Effects 0.000 claims description 4
- 239000004035 construction material Substances 0.000 claims description 3
- 230000008595 infiltration Effects 0.000 abstract description 8
- 238000001764 infiltration Methods 0.000 abstract description 8
- 238000005259 measurement Methods 0.000 abstract description 5
- 238000004458 analytical method Methods 0.000 description 7
- 230000008901 benefit Effects 0.000 description 4
- 238000007654 immersion Methods 0.000 description 4
- 238000004519 manufacturing process Methods 0.000 description 4
- 238000009933 burial Methods 0.000 description 3
- 238000010276 construction Methods 0.000 description 3
- 238000012795 verification Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 238000007599 discharging Methods 0.000 description 2
- 229910052500 inorganic mineral Inorganic materials 0.000 description 2
- 239000010410 layer Substances 0.000 description 2
- 239000011707 mineral Substances 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 239000013049 sediment Substances 0.000 description 2
- 239000002689 soil Substances 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 239000004575 stone Substances 0.000 description 2
- 238000009825 accumulation Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000000151 deposition Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000003912 environmental pollution Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 229910052755 nonmetal Inorganic materials 0.000 description 1
- 230000003204 osmotic effect Effects 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000002791 soaking Methods 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
- 239000002344 surface layer Substances 0.000 description 1
- 239000010878 waste rock Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Computing Systems (AREA)
- Fluid Mechanics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses an optimization method for tailings pond seepage calculation based on an upstream boundary condition, which comprises the following steps: s1, obtaining basic data of a tailing pond, and establishing a seepage calculation model based on the basic data; s2, determining a seepage calculation upstream boundary condition range of the tailing pond based on the fixed water head boundary line; s3, adding the seepage calculation upstream boundary condition according to the seepage calculation upstream boundary condition range; s4, calculating seepage of the tailing pond, and analyzing and calculating to obtain the spatial distribution of the seepage surface of the tailing pond. The invention considers the influence of the ore drawing water on the position of the infiltration surface, more reasonably determines the infiltration boundary of the tailing pond, and enables the calculation result to be more approximate to the actual measurement value.
Description
Technical Field
The invention relates to the technical field of mine construction, in particular to an optimization method for seepage calculation of a tailing pond based on upstream boundary conditions.
Background
The tailing pond is a special industrial building, which is one of three basic projects of mines. So far, various tailing ponds in use in the world are over 2 tens of thousands, and nearly ten thousands exist in China. The operation condition of the tailing pond not only relates to the production safety and economic benefit of the mine, but also is closely related to life, property and surrounding environment of residents at the downstream of the storage area. The geographical location of many tailings ponds is very important, some are located upstream of great rivers, great lakes and important water sources, some are located upstream of important public transportation facilities, and even some are located upstream of dense residential areas, so that serious casualties and property loss are caused once a dam collapses, and serious environmental pollution is generated. Thus, the tailing pond is an important production facility of mines and is an important hazard source. Cases of serious damage are frequent due to engineering failure of tailing ponds. With the development of industry, the number of tailing ponds is increased, and the accumulation of the dam body is also increased, so that the safety problem of the tailing dam is more and more important.
It is well known that seepage conditions in tailings ponds are one of the important factors affecting the stability of tailings pond dams. Meanwhile, in the safety and stability analysis of the tailing pond, most of dam break accidents of the tailing pond are directly related to the too high position of the seepage surface of the dam body. Therefore, the characteristics of the seepage field of the tailing pond are accurately analyzed, the seepage infiltration surface (line) position is accurately simulated and analyzed, the method is a precondition for analyzing the dam slope stability and the dam static power of the tailing pond, and meanwhile, the method is also very important for the design, construction and later operation management of the tailing pond.
In the seepage calculation of the tailings pond, whether the simulated space position of the seepage surface of the tailings pond is close to the actual space position or not, wherein one key influence factor is the setting of the seepage calculation boundary condition, particularly the determination of the upstream boundary condition. At present, the method for determining the upstream boundary condition in the seepage calculation adopts the distribution of the sedimentary beach surface in the tailing pond as the upstream boundary condition of the seepage calculation, and the method can lead to the lower position of the determined seepage surface in the seepage calculation than the actual position because of the influence of the mineral water, because the tailing pond is provided with a mineral placement pipe at the top of a dam, the tailing pond continuously discharges tailing pulp into the pond, and part of water infiltrates into the dam body at the dry beach position in the process of the pulp flowing to the tailing pond, thereby influencing the spatial distribution of the position of the seepage surface in the pond. Therefore, a more reasonable method for determining the boundary conditions of seepage calculation is put forward to be studied, and is very important for accurately simulating and analyzing the spatial distribution of the seepage surface.
On the basis of the existing method, the invention provides an optimization method for tailings pond seepage calculation based on an upstream boundary condition. Meanwhile, through the comparison and verification of a plurality of engineering examples, the space distribution of the infiltration surface simulated by the novel method is more consistent with the actual situation.
Disclosure of Invention
Aiming at the problems in the related art, the invention provides an optimization method for seepage calculation of a tailing pond based on an upstream boundary condition, so as to overcome the technical problems in the prior art.
For this purpose, the invention adopts the following specific technical scheme:
an optimization method for tailings pond seepage calculation based on upstream boundary conditions, which comprises the following steps:
s1, obtaining basic data of a tailing pond, and establishing a seepage calculation model based on the basic data;
s2, determining a seepage calculation upstream boundary condition range of the tailing pond based on the fixed water head boundary line;
s3, adding the seepage calculation upstream boundary condition according to the seepage calculation upstream boundary condition range;
s4, calculating seepage of the tailing pond, and analyzing and calculating to obtain the spatial distribution of the seepage surface of the tailing pond.
Further, the step of obtaining basic data of the tailing pond and establishing a seepage calculation model based on the basic data comprises the following steps:
s11, acquiring basic design data of a tailing pond;
s12, establishing a three-dimensional numerical calculation model of the tailing pond by utilizing basic design data in finite element analysis software;
the three-dimensional numerical calculation model is established for simulating a tailing pond when the stacking elevation reaches the actual operation elevation.
Further, the basic design data of the tailing pond comprise the dam height, the inner side slope and outer side slope ratio, the platform width, the designed normal water level value, the actual running water level line position, the actual water level elevation value, the topography condition of the tailing pond, the dam construction material of the initial dam, the volume weight, the cohesive force and the physical and mechanical parameters and the permeability coefficient of the tailing sand in the pond.
Further, the range of the upstream boundary condition of the seepage calculation is defined by a fixed water head boundary line I and a fixed water head boundary line II.
Further, the determining the seepage calculation upstream boundary condition range of the tailing pond based on the constant head boundary line comprises the following steps:
s21, determining a first fixed water head boundary line according to an actual running water level elevation value of the tailing pond;
s22, determining a fixed water head boundary line II according to the actual running water level elevation value and the actual running water level line position;
s23, the range of the upstream boundary condition of the seepage calculation is defined by utilizing the first fixed water head boundary line and the second fixed water head boundary line.
Further, the step of determining the boundary line I of the constant head according to the elevation value of the actual running water level of the tailing pond comprises the following steps:
s211, acquiring an actual water level elevation value of the tailing pond according to the basic data of the tailing pond;
s212, making a horizontal plane with the same elevation value as the actual water level elevation value to obtain a horizontal plane I;
s213, analyzing to obtain an intersecting line of the horizontal plane I and the mountain topography of the tail part of the tailing pond, wherein the intersecting line is a constant head boundary line I of the tailing pond.
Further, the step of determining the fixed water head boundary line II according to the actual running water level elevation value and the actual running water level line position comprises the following steps:
s221, acquiring the actual running water line position of the tailing pond according to the basic data of the tailing pond;
s222, obtaining the distance L between the actual running water line position and the current beach top of the tailing pond according to the actual running water line position;
s223, reducing the length of the dry beach, and translating the actual running water line to a distance D along the current beach top direction of the tailing pond to obtain a fixed water head boundary line II, wherein the calculation formula of the translation distance D is as follows:
further, the range of the upstream boundary condition calculated by the seepage is defined by the first fixed water head boundary line and the second fixed water head boundary line, and the range comprises the following steps:
and a closed area is defined by the fixed water head boundary line I and the fixed water head boundary line II on the horizontal plane I, and the closed area is the range of the upstream boundary condition calculated for the seepage of the tailing pond to be determined.
Further, the adding the seepage calculation upstream boundary condition according to the seepage calculation upstream boundary condition range comprises the following steps:
and adding an upstream water head in the range of the upstream boundary condition calculated by the seepage, wherein the water head value is an actual running water level elevation value.
Further, the steps of calculating the seepage of the tailing pond and analyzing and calculating to obtain the spatial distribution of the seepage surface of the tailing pond include the following steps:
s41, adding an upstream boundary condition of seepage calculation in the three-dimensional numerical calculation model;
s42, performing seepage calculation on the tailings pond by using finite element analysis software, and analyzing and calculating to obtain the spatial distribution of the seepage surface of the tailings pond.
The beneficial effects of the invention are as follows: the invention considers the influence of the ore drawing water on the position of the infiltration surface, so that the calculation result is more approximate to the actual measurement value. By carrying out seepage calculation on a large number of actual tailing pond projects and combining with the contrast verification of the engineering field actual measurement seepage line data, the seepage calculation boundary conditions determined by the method can more truly simulate the influence of the ore drawing water on the seepage surface, the space position of the seepage surface obtained by analysis calculation is more accurate and is more approximate to the actual measured value, and especially the critical areas (dam tops and dam slopes) considering the stability of the tailing pond are almost coincident with the actual measured value.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings that are needed in the embodiments will be briefly described below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and other drawings may be obtained according to these drawings without inventive effort for a person skilled in the art.
FIG. 1 is a schematic flow diagram of an optimization method for tailings pond seepage calculation based on upstream boundary conditions in accordance with an embodiment of the present invention;
fig. 2 is a schematic diagram of an optimization method of tailings pond seepage calculation based on an upstream boundary condition according to an embodiment of the present invention.
In the figure:
1. actually running a water line; 2. actually running a beach top line; 3. a fixed water head boundary line I; 4. a fixed water head boundary line II; 5. a first horizontal plane; D. translation distance.
Detailed Description
For the purpose of further illustrating the various embodiments, the present invention provides the accompanying drawings, which are a part of the disclosure of the present invention, and which are mainly used to illustrate the embodiments and, together with the description, serve to explain the principles of the embodiments, and with reference to these descriptions, one skilled in the art will recognize other possible implementations and advantages of the present invention, wherein elements are not drawn to scale, and like reference numerals are generally used to designate like elements.
For a better understanding of the technical solutions of the present invention, the following terms are explained with reference to the present invention as follows:
tailing pond: the device is used for storing metal and nonmetal mines, and discharging tailings after ore sorting.
Initial dam: is constructed by soil, stone materials and the like and is used as a seepage-discharging or supporting body of the tailing accumulating dam.
Tailing stacking dam: a dam formed by piling up tailings in the production process.
Infiltration line: the free surface of the seepage water in the dam body is a curve in the cross section.
And (3) soaking the surface: the location of the free surface of the seepage water in the dam.
Normal water level: the water level meeting the production backwater and discharge requirements in the tailing pond can be realized.
Length of dry beach: horizontal distance from the water edge line in the warehouse to the beach top.
Deposit beach: the surface layer of a sediment formed by the hydraulic alluvial tailings is divided into an upper water part and a lower water part according to the water surface of a water collecting area in the reservoir.
Beach top: and depositing the intersection line of the beach surface and the slope surface outside the sub-dam.
According to the embodiment of the invention, an optimization method for tailings pond seepage calculation based on an upstream boundary condition is provided.
The invention will be further described with reference to the accompanying drawings and detailed description, as shown in fig. 1-2, an optimization method for tailings pond seepage calculation based on upstream boundary conditions according to an embodiment of the invention, the optimization method comprising the following steps:
s1, obtaining basic data of a tailing pond, and establishing a seepage calculation model based on the basic data;
the method for obtaining the basic data of the tailing pond and establishing the seepage calculation model based on the basic data comprises the following steps of:
s11, acquiring basic design data of a tailing pond;
specifically, the basic data are basic design data of the tailing pond, and the basic design data of the tailing pond comprise physical and mechanical parameters such as dam height of an initial dam and a heap dam, inner and outer slope ratios, platform width, design normal water level values, actual running water level line positions, actual water level elevation values, topography conditions of the tailing pond, initial dam construction materials, volume weight of tailings in the pond, cohesive force, internal friction angle and the like, permeability coefficients and the like.
S12, establishing a three-dimensional numerical calculation model of the tailing pond by utilizing basic design data in finite element analysis software; the three-dimensional numerical calculation model is established for simulating a tailing pond when the stacking elevation reaches the actual operation elevation.
S2, determining a seepage calculation upstream boundary condition range of the tailing pond based on the fixed water head boundary line;
the range of the upstream boundary condition of the seepage calculation is composed of two parts, namely a fixed water head boundary line I3 and a fixed water head boundary line II 4, which are positioned on the same horizontal plane and jointly encircle the range of the upstream boundary condition of the seepage calculation.
Specifically, the determining the seepage calculation upstream boundary condition range of the tailing pond based on the constant head boundary line comprises the following steps:
s21, determining a fixed water head boundary line I3 according to the actual running water level elevation value of the tailing pond;
specifically, the step of determining the fixed water head boundary line 3 according to the actual running water level elevation value of the tailing pond comprises the following steps:
s211, acquiring an actual water level elevation value of the tailing pond according to the basic data of the tailing pond;
s212, making a horizontal plane with the same elevation value as the actual water level elevation value to obtain a horizontal plane I5;
s213, analyzing to obtain an intersecting line of the horizontal plane I5 and the mountain topography at the tail part of the tailing pond, wherein the intersecting line is a constant head boundary line I3 of the tailing pond.
S22, determining a fixed water head boundary line II 4 according to the actual running water level elevation value and the actual running water level line position;
specifically, the step of determining the fixed water head boundary line two 4 according to the actual running water level elevation value and the actual running water level line position includes the following steps:
s221, acquiring the actual running water line position of the tailing pond according to the basic data of the tailing pond;
s222, obtaining the distance L between the actual running water line position and the current beach top of the tailing pond according to the actual running water line position;
s223, translating the actual running water level boundary line to a distance D along the current beach top direction of the tailing pond to obtain a fixed water head boundary line II 4, wherein the calculation formula of the translation distance D is as follows:
s23, the range of the upstream boundary condition of the seepage calculation is defined by the aid of the fixed water head boundary line I3 and the fixed water head boundary line II 4.
Specifically, the range for calculating the upstream boundary condition by utilizing the constant head boundary line I3 and the constant head boundary line II 4 together to define the seepage flow comprises the following steps:
and a closed area is defined on the horizontal plane I5 by the fixed water head boundary line I3 and the fixed water head boundary line II 4, and the closed area is the range of the upstream boundary condition calculated for the seepage of the tailing pond to be determined.
S3, adding the seepage calculation upstream boundary condition according to the seepage calculation upstream boundary condition range;
wherein the adding the seepage calculation upstream boundary condition according to the seepage calculation upstream boundary condition range comprises the following steps:
and adding an upstream water head in the range of the upstream boundary condition calculated by the seepage, wherein the water head value is an actual running water level elevation value.
S4, calculating seepage of the tailing pond, and analyzing and calculating to obtain the spatial distribution of the seepage surface of the tailing pond.
The method for calculating the seepage of the tailing pond and analyzing and calculating the spatial distribution of the seepage surface of the tailing pond comprises the following steps:
s41, adding an upstream boundary condition of seepage calculation in the three-dimensional numerical calculation model;
s42, performing seepage calculation on the tailings pond by using finite element analysis software, and analyzing and calculating to obtain the spatial distribution of the seepage surface of the tailings pond.
In order to facilitate understanding of the above technical solutions of the present invention, the following details are about specific embodiments of the present invention in practical processes.
DETAILED DESCRIPTION OF EMBODIMENT (S) OF INVENTION
S1, obtaining basic data of a tailing pond, and establishing a three-dimensional numerical calculation model based on the basic data:
a certain tailing pond in Hebei province is located in a valley between mountains with V-shaped cross sections, the initial dam of the tailing pond is a waste rock pile dam body, the elevation of the dam top is 654.0m, the elevation of the dam bottom is 637.0m, the dam length is about 74m, the outer dam slope ratio of the initial dam is 1:1.75, the inner dam slope ratio is 1:1.70, and the outer dam slope adopts dry masonry stone for slope protection. The average outer dam slope ratio of the dam is about 1:3.0 when the beach top standard height reaches 745.0 m. The corresponding actual running water level value is 740.54m when the beach top standard height 745.0 m. The average slope ratio of the deposit beach is about 1:60.3. As the tailing pond is an existing pond, the permeability index of each rock and soil layer is obtained through a geological survey report, and the details are shown in Table 1.
TABLE 1 osmotic coefficient reference values for each petrography layer in reservoir area
According to the basic data, a three-dimensional numerical calculation model of a tailing pond area and a dam body is established in finite element analysis software, seepage calculation under actual operation working conditions is carried out on the tailing pond when the stacking elevation reaches 745.0m, and the space distribution form of the tailing stacking dam body and the infiltration surface in a sediment beach is simulated.
S2, determining a seepage calculation upstream boundary condition range of the tailing pond based on the fixed water head boundary line;
1) Determining a boundary line of a constant head
Firstly, according to the basic data, the actual running water level elevation of the tailing pond is 740.54m, and a horizontal plane I with the elevation of 740.54m is made. And the horizontal plane I is intersected with the mountain terrain at the tail part in the warehouse to obtain a fixed water head boundary line I.
2) Determining a fixed water head boundary line II
Firstly, according to the position of the actual running water line of the tailing pond, the distance L between the actual running water line of the tailing pond and the current beach top of the tailing pond is 268.94m. Substituting the L into a calculation formula of the translation distance D to obtain the translation distance D of 185.53m.
And further, translating the actual running water line to 185.53m in the current beach top direction of the tailing pond, and obtaining a fixed water head boundary line II.
3) Determining the range of boundary conditions upstream of a seepage calculation
The first fixed water head boundary line and the second fixed water head boundary line jointly encircle a closed area on the first horizontal plane, so that the range of the upstream boundary condition calculated by the seepage of the tailing pond is obtained.
S3, adding the seepage calculation upstream boundary condition according to the seepage calculation upstream boundary condition range;
the upstream boundary conditions are added according to the boundary condition range. An upstream constant head is added within the range of the upstream boundary condition, and the head value is 740.54m.
S4, calculating seepage of the tailing pond, and analyzing and calculating to obtain the spatial distribution of the seepage surface of the tailing pond;
further, according to the determined upstream boundary range and the water head value, an upstream boundary condition is added to the three-dimensional numerical calculation model of the tailing pond. And adding other calculation conditions and attributes to the three-dimensional numerical calculation model of the tailing pond according to the basic data. Further, seepage calculation is carried out on the tailing pond by utilizing finite element analysis software. Through analysis and calculation, the spatial distribution form of the infiltration surface in the tailing pond can be simulated.
In order to better demonstrate the advantages of the method provided by the invention, the calculation result of the embodiment is compared with the calculation result and the actual observation result of the original method. To better show the difference between the two calculated results and the measured values, the immersion line burial depths at different positions of the main section of the tailing pond are compared, as shown in table 2.
TABLE 2 comparison of immersion line burial depths
The method has the advantages that the burial depth of the immersion line calculated by the method is closer to the actual observed value, and the space distribution of the immersion surface obtained by analysis and simulation is closer to the actual value, so that the boundary condition determined by the method is more reasonable. The method for calculating seepage of the tailing pond can better provide technical support for dam slope stability and dam static and dynamic analysis of the tailing pond, and simultaneously provide more reliable analysis data for design, construction and later operation management of the tailing pond.
In summary, by means of the above technical solution of the present invention, the present invention considers the influence of the ore drawing water on the position of the wetted surface, so that the calculation result is closer to the actual measurement value. By carrying out seepage calculation on a large number of actual tailing pond projects and combining with the contrast verification of the engineering field actual measurement seepage line data, the seepage calculation boundary conditions determined by the method can consider the influence of the ore drawing water on the seepage surface in a simple and convenient way, the space position of the seepage surface obtained by analysis calculation is more accurate and is more approximate to the actual measured value, and especially the critical areas (dam tops and dam slopes) considering the stability of the tailing pond are almost coincident with the actual measured value.
The foregoing description of the preferred embodiments of the invention is not intended to be limiting, but rather is intended to cover all modifications, equivalents, alternatives, and improvements that fall within the spirit and scope of the invention.
Claims (10)
1. An optimization method for tailings pond seepage calculation based on upstream boundary conditions is characterized by comprising the following steps:
s1, obtaining basic data of a tailing pond, and establishing a seepage calculation model based on the basic data;
s2, determining a seepage calculation upstream boundary condition range of the tailing pond based on the fixed water head boundary line;
s3, adding the seepage calculation upstream boundary condition according to the seepage calculation upstream boundary condition range;
s4, calculating seepage of the tailing pond, and analyzing and calculating to obtain the spatial distribution of the seepage surface of the tailing pond.
2. The method for optimizing seepage calculation of a tailings pond based on an upstream boundary condition according to claim 1, wherein the steps of obtaining basic data of the tailings pond and establishing a seepage calculation model based on the basic data comprise the following steps:
s11, acquiring basic design data of a tailing pond;
s12, establishing a three-dimensional numerical calculation model of the tailing pond by utilizing basic design data in finite element analysis software;
the three-dimensional numerical calculation model is established for simulating a tailing pond when the stacking elevation reaches the actual operation elevation.
3. The optimization method of tailings pond seepage calculation based on the upstream boundary condition according to claim 1, wherein the basic design data of the tailings pond comprises the dam height, the inner and outer slope ratios, the platform width, the design normal water level value, the actual running water level line position, the actual water level elevation value, the topography condition of the tailings pond, the dam construction material of the initial dam, the volume weight, the cohesive force, the physical and mechanical parameters of the internal friction angle and the permeability coefficient of the tailings in the pond.
4. The method for optimizing tailings pond seepage calculation based on the upstream boundary condition according to claim 1, wherein the seepage calculation upstream boundary condition range is defined by a constant head boundary line I and a constant head boundary line II.
5. The method for optimizing tailings pond seepage calculation based on the upstream boundary condition of claim 4 wherein the determining the range of the seepage calculation upstream boundary condition of the tailings pond based on the constant head boundary line comprises the steps of:
s21, determining a first fixed water head boundary line according to an actual running water level elevation value of the tailing pond;
s22, determining a fixed water head boundary line II according to the actual running water level elevation value and the actual running water level line position;
s23, the range of the upstream boundary condition of the seepage calculation is defined by utilizing the first fixed water head boundary line and the second fixed water head boundary line.
6. The optimization method for tailings pond seepage calculation based on the upstream boundary condition according to claim 5, wherein the step of determining the constant head boundary line according to the actual running water level elevation value of the tailings pond comprises the following steps:
s211, acquiring an actual water level elevation value of the tailing pond according to the basic data of the tailing pond;
s212, making a horizontal plane with the same elevation value as the actual water level elevation value to obtain a horizontal plane I;
s213, analyzing to obtain an intersecting line of the horizontal plane I and the mountain topography of the tail part of the tailing pond, wherein the intersecting line is a constant head boundary line I of the tailing pond.
7. The optimization method for tailings pond seepage calculation based on the upstream boundary condition according to claim 6, wherein the step of determining the constant head boundary line II according to the actual running water level elevation value and the actual running water level line position comprises the following steps:
s221, acquiring the actual running water line position of the tailing pond according to the basic data of the tailing pond;
s222, obtaining the distance L between the actual running water line position and the current beach top of the tailing pond according to the actual running water line position;
s223, reducing the length of the dry beach, and translating the actual running water line to a distance D along the current beach top direction of the tailing pond to obtain a fixed water head boundary line II, wherein the calculation formula of the translation distance D is as follows:
8. the method for optimizing seepage calculation of tailings pond based on the upstream boundary condition according to claim 5, wherein the range for obtaining the seepage calculation upstream boundary condition by the joint delineation of the fixed water head boundary line I and the fixed water head boundary line II comprises the following steps:
and a closed area is defined by the fixed water head boundary line I and the fixed water head boundary line II on the horizontal plane I, and the closed area is the range of the upstream boundary condition calculated for the seepage of the tailing pond to be determined.
9. The method for optimizing tailings pond seepage calculation based on the upstream boundary condition according to claim 1, wherein the adding the seepage calculation upstream boundary condition according to the seepage calculation upstream boundary condition range comprises the steps of:
and adding an upstream water head in the range of the upstream boundary condition calculated by the seepage, wherein the water head value is an actual running water level elevation value.
10. The optimization method for tailings pond seepage calculation based on the upstream boundary condition according to claim 2, wherein the steps of calculating the seepage of the tailings pond and analyzing and calculating to obtain the spatial distribution of the seepage surface of the tailings pond comprise the following steps:
s41, adding an upstream boundary condition of seepage calculation in the three-dimensional numerical calculation model;
s42, performing seepage calculation on the tailings pond by using finite element analysis software, and analyzing and calculating to obtain the spatial distribution of the seepage surface of the tailings pond.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310536901.1A CN116306173B (en) | 2023-05-13 | 2023-05-13 | Optimization method for seepage calculation of tailing pond based on upstream boundary conditions |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310536901.1A CN116306173B (en) | 2023-05-13 | 2023-05-13 | Optimization method for seepage calculation of tailing pond based on upstream boundary conditions |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN116306173A true CN116306173A (en) | 2023-06-23 |
| CN116306173B CN116306173B (en) | 2023-08-04 |
Family
ID=86789028
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310536901.1A Active CN116306173B (en) | 2023-05-13 | 2023-05-13 | Optimization method for seepage calculation of tailing pond based on upstream boundary conditions |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116306173B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117251992A (en) * | 2023-08-28 | 2023-12-19 | 华南理工大学 | Calculation method of infiltration line of anchorage circular foundation pit |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109800500A (en) * | 2019-01-16 | 2019-05-24 | 中国恩菲工程技术有限公司 | Storage medium, the calculation method of Tailings Dam analysis of seepage flow stability and its device |
| WO2021003690A1 (en) * | 2019-07-10 | 2021-01-14 | 中国科学院地质与地球物理研究所 | Tailings pond dam burst disaster simulation system and method |
| CN114707222A (en) * | 2022-04-19 | 2022-07-05 | 南京信息工程大学 | Seepage-considered tailing dam body stability real-time online analysis and prediction method |
| CN115099092A (en) * | 2022-06-24 | 2022-09-23 | 西安理工大学 | Tailing pond seepage calculation method based on three-dimensional modeling |
-
2023
- 2023-05-13 CN CN202310536901.1A patent/CN116306173B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109800500A (en) * | 2019-01-16 | 2019-05-24 | 中国恩菲工程技术有限公司 | Storage medium, the calculation method of Tailings Dam analysis of seepage flow stability and its device |
| WO2021003690A1 (en) * | 2019-07-10 | 2021-01-14 | 中国科学院地质与地球物理研究所 | Tailings pond dam burst disaster simulation system and method |
| CN114707222A (en) * | 2022-04-19 | 2022-07-05 | 南京信息工程大学 | Seepage-considered tailing dam body stability real-time online analysis and prediction method |
| CN115099092A (en) * | 2022-06-24 | 2022-09-23 | 西安理工大学 | Tailing pond seepage calculation method based on three-dimensional modeling |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117251992A (en) * | 2023-08-28 | 2023-12-19 | 华南理工大学 | Calculation method of infiltration line of anchorage circular foundation pit |
| CN117251992B (en) * | 2023-08-28 | 2024-05-24 | 华南理工大学 | Calculation method of infiltration line of anchorage circular foundation pit |
Also Published As
| Publication number | Publication date |
|---|---|
| CN116306173B (en) | 2023-08-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Essaid | A multilayered sharp interface model of coupled freshwater and saltwater flow in coastal systems: Model development and application | |
| CN108647473A (en) | A kind of shield method builds seabed tunnel Rational Depth computational methods | |
| CN116306173B (en) | Optimization method for seepage calculation of tailing pond based on upstream boundary conditions | |
| CN105810075A (en) | Water-pumping triggered karst collapse process experimental device | |
| CN110929390B (en) | Numerical simulation detection method based on groundwater hydrogeology test | |
| CN112836270B (en) | Method for predicting influence of diving precipitation on building settlement | |
| CN108446489A (en) | Measurement method and its processing unit of the Coastline Changes to wetland influence on groundwater | |
| CN103437318B (en) | The settlement after construction computational prediction method of floated sea wall | |
| CN109284567A (en) | A kind of buried fault acts on the research method of lower deeply covered layer cut-pff wall stress state | |
| CN103614984A (en) | Ecological treatment method oriented to geological environment disasters caused by riverway goafs | |
| Toumi et al. | Evaluation of geology and hydrogeology of the water leakage in Hammam-Grouz Dam, Algeria | |
| Jaswal et al. | Delineation of phreatic surface in soil type slope—a comparative study using physical and numerical modeling | |
| Polomčić et al. | Application of Groundwater modeling for designing a dewatering system: Case study of the Buvač Open Cast Mine, Bosnia and Herzegovina | |
| Agarkov et al. | Optimization of the drainage system of overburden dumps using geofiltration modeling | |
| CN109147047B (en) | Seepage channel type riverside water source and place numerical model construction method | |
| Yu et al. | A groundwater flow model for excavation of coastal earthen sites: Case study Jintoushan site, China | |
| Ogilvy et al. | The possibilities of geophysical methods applied for investigating the impact of man on the geological medium | |
| Kohfahl et al. | Modelling of artificial water oscillations in a flooded lignite mine | |
| Bahgat | Optimum use of dredged materials for sustainable shoreline management in Nile Delta | |
| Ingles | Geomorphology as the prime determinant of current engineering problems in the Tamar River basin at Launceston, Tasmania | |
| Zhang et al. | Research on the Effective Control of Ground Settlement during Double‐Layered Foundation Pit Dewatering Based on Seepage Control‐Recharge Coupling Model | |
| Saretzky | Hydrological characterization of a sulphide waste rock dump | |
| Yi et al. | Two-dimensional simulation of underground seepage in a dangerous piping zone of the Jingjiang Great Levee, the middle reach of the Yangtze River | |
| Xiao et al. | Dewatering Design of a Foundation Pit | |
| Fraile del Río | Hydrogeological Model of Aitik Mine Using Leapfrog Geo Software |
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 |