CN116628402B - Tailing pond flood regulating calculation method and device, electronic equipment and storage medium - Google Patents

Tailing pond flood regulating calculation method and device, electronic equipment and storage medium Download PDF

Info

Publication number
CN116628402B
CN116628402B CN202310890320.8A CN202310890320A CN116628402B CN 116628402 B CN116628402 B CN 116628402B CN 202310890320 A CN202310890320 A CN 202310890320A CN 116628402 B CN116628402 B CN 116628402B
Authority
CN
China
Prior art keywords
flood
water level
discharge flow
regulating reservoir
calculation
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
CN202310890320.8A
Other languages
Chinese (zh)
Other versions
CN116628402A (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.)
BGRIMM Technology Group Co Ltd
Original Assignee
BGRIMM Technology Group Co Ltd
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 BGRIMM Technology Group Co Ltd filed Critical BGRIMM Technology Group Co Ltd
Priority to CN202310890320.8A priority Critical patent/CN116628402B/en
Publication of CN116628402A publication Critical patent/CN116628402A/en
Application granted granted Critical
Publication of CN116628402B publication Critical patent/CN116628402B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/15—Correlation function computation including computation of convolution operations
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063—Operations research, analysis or management
    • G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/02—Agriculture; Fishing; Forestry; Mining
    • 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)
  • General Physics & Mathematics (AREA)
  • Business, Economics & Management (AREA)
  • Theoretical Computer Science (AREA)
  • Human Resources & Organizations (AREA)
  • Data Mining & Analysis (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Strategic Management (AREA)
  • Economics (AREA)
  • General Engineering & Computer Science (AREA)
  • Software Systems (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Tourism & Hospitality (AREA)
  • Marketing (AREA)
  • Operations Research (AREA)
  • Databases & Information Systems (AREA)
  • Algebra (AREA)
  • General Business, Economics & Management (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Probability & Statistics with Applications (AREA)
  • Educational Administration (AREA)
  • Development Economics (AREA)
  • Game Theory and Decision Science (AREA)
  • Evolutionary Biology (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Quality & Reliability (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Computing Systems (AREA)
  • Agronomy & Crop Science (AREA)
  • Animal Husbandry (AREA)
  • Marine Sciences & Fisheries (AREA)
  • Mining & Mineral Resources (AREA)
  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)

Abstract

The application provides a tailing pond flood regulating algorithm method, a tailing pond flood regulating algorithm device, electronic equipment and a storage medium, and relates to the technical field of data processing. The method comprises the following steps: determining a flood regulating reservoir capacity function according to the relation between the water level and the flood regulating reservoir capacity; determining a flood discharge flow function according to the relationship between the water level and the flood discharge flow; acquiring a water level step and a time step; calculating the variable quantity of the capacity of the flood regulating reservoir based on the capacity function of the flood regulating reservoir according to the step length of the water level; calculating the variation of the flood discharge flow based on the flood discharge flow function according to the water level step and the time step; determining the tolerance of iterative computation convergence of flood regulating calculation according to the variable quantity of the flood regulating reservoir capacity and the variable quantity of the flood discharging flow; and carrying out iterative computation of flood diversion calculation for a plurality of times until the absolute value of the computation result is smaller than the tolerance, stopping the iterative computation, and determining the currently obtained computation result as the calculation result of the flood diversion calculation of the tailing pond. The application can effectively solve the problem that flood regulating calculation in the prior art can not be converged well.

Description

Tailing pond flood regulating calculation method and device, electronic equipment and storage medium
Technical Field
The application relates to the technical field of data processing, in particular to a tailing pond flood regulating calculation method, a tailing pond flood regulating calculation device, electronic equipment and a storage medium.
Background
The tailing pond is an important facility in mine production, and the flood control capability of the tailing pond has special significance for the safe production of the tailing pond. According to the requirements of the safety regulations of tailing ponds (GB 39496-2020), the tailing ponds should carry out flood control calculation in the design stage and the operation stage so as to determine whether the flood control capability of tailings meets the safety requirement, especially, the flood control calculation is carried out by entrusted design units before each year according to the actual measured topography of the tailing ponds, the water level and the actual condition of the tailing deposit beach surface, the flood control capability of the tailing ponds is rechecked, and the operation water level, the dry beach length, the safety superelevation and other safety operation control parameters of the tailing ponds in the flood season are determined, so that the flood control calculation is a necessary safety work throughout the whole life cycle of the tailing ponds.
According to the requirements of the safety regulations of tailing ponds (GB 39496-2020) and the design specifications of tailing facilities (GB 50863-2013), the flood regulating algorithm of the tailing ponds should adopt a water balance method. The method for directly solving the water balance formula is easy to generate root increment, needs manual screening, and has extremely large manual screening workload for the calculation process of shorter unit time and longer flood duration. Therefore, the currently commonly adopted operation method is to set the water level of the unit step length for cyclic iterative calculation, set the tolerance at the same time, and consider the calculation result to be converged when the absolute value of the difference value of the flood bank capacity between the final moment and the initial moment of the unit time period is smaller than or equal to the tolerance.
For computer programming solution, a loop iteration method is more commonly used, but loop iteration also has problems: when machine calculation is involved, the tolerance is taken as a convergence criterion, but no clear method is used for determining the convergence criterion all the time, and the method can be manually set by experience. When the tolerance setting is unreasonable, if the tolerance is too small, the calculation result is scattered and cannot be solved; too large a tolerance would result in too early convergence of the calculation, resulting in a larger error in the calculation result.
Disclosure of Invention
In order to solve the technical problems, the method, the device, the electronic equipment and the storage medium for the flood control algorithm of the tailing pond are provided, and reasonable convergence criteria can be determined, so that the circulation iteration of the flood control algorithm of the tailing pond can be converged better.
Embodiments of the application may be implemented as follows:
in a first aspect, an embodiment of the present application provides a method for calculating flood diversion of a tailing pond, where the method includes:
determining a flood regulating reservoir capacity function according to the relation between the water level and the flood regulating reservoir capacity;
determining a flood discharge flow function according to the relationship between the water level and the flood discharge flow;
acquiring a water level step and a time step;
calculating the variable quantity of the flood regulating reservoir capacity based on the flood regulating reservoir capacity function according to the water level step length;
Calculating the variation of the flood discharge flow based on the flood discharge flow function according to the water level step and the time step;
determining the tolerance of iterative computation convergence of the flood regulating algorithm according to the variable quantity of the flood regulating reservoir capacity and the variable quantity of the flood discharging flow;
and carrying out iterative computation of flood diversion calculation for a plurality of times until the absolute value of the computation result is smaller than the tolerance, stopping the iterative computation, and determining the currently obtained computation result as the calculation result of the flood diversion calculation of the tailing pond.
In an embodiment, the performing the iterative computation of the flood diversion algorithm for a plurality of times until the absolute value of the computation result is smaller than the tolerance, and stopping the iterative computation includes:
the iterative calculation of the flood diversion calculation is performed according to the following formula,
wherein ,ΔTthe time period of the preset time period is indicated,Q s andQ z hong Liuliang indicating the initial and final moments of the preset time period respectively,q s andq z respectively representing the flood discharge flow at the initial moment and the final moment of the preset time period,V s andV z respectively representing the flood regulating reservoir capacity at the initial moment and the final moment of a preset time period;q s andq z determined according to the flood discharge flow function,V s andV z determining according to the flood regulating reservoir capacity function;
when (when)V s -V z And stopping the iterative calculation when the absolute value of (a) is smaller than the tolerance.
In one embodiment, the determining the flood diversion reservoir capacity function according to the water level and the flood diversion reservoir capacity relation comprises:
Obtaining a contour topography of a tailing pond;
acquiring the area of each contour line between the rising and adjusting water level and the beach top elevation of the contour line topographic map;
obtaining an average value of areas of adjacent contour lines, and multiplying the average value of the areas of the adjacent contour lines by a height difference of the adjacent contour lines to obtain a volume of the adjacent contour lines;
accumulating the volumes of all the adjacent contour lines between the lifting water level and the beach elevation to obtain the relation between the water level and the flood regulating reservoir capacity;
and fitting the water level and the flood regulating reservoir capacity relation to obtain the flood regulating reservoir capacity function.
In an embodiment, the determining the flood discharge flow function according to the relationship between the water level and the flood discharge flow comprises:
obtaining a flood drainage system arrangement form, wherein the flood drainage system arrangement form comprises at least one of a pit-type, a trench-type and a spillway;
and determining the flood discharge flow function according to at least one of a first preset formula corresponding to the well hole type, a second preset formula corresponding to the groove hole type and a third preset formula corresponding to the spillway and the arrangement form of the flood discharge system.
In an embodiment, the calculating the variation of the capacity of the flood diversion reservoir based on the capacity function of the flood diversion reservoir according to the water level step comprises:
The flood regulating reservoir capacity value is calculated according to the following formula,
representing the value of the flood regulating reservoir->Representing the flood regulating reservoir capacity function, wherein the definition domain of the flood regulating reservoir capacity function is [ [h 0 ，h 1 ]The value range of n is consistent with [0, ]>]Integer of deltahRepresents the step size of the water level,h 0 indicating the start-up water level to be adjusted,h 1 representing beach elevation;
calculating the difference value of the adjacent water level step-length flood regulating reservoir capacity according to the following formula to obtain the variation quantity of the flood regulating reservoir capacity,
indicating the variation of the capacity of the flood regulating reservoir, wherein the value range of n is in accordance with [0, ], and->]Is an integer of (a).
In an embodiment, the calculating the variation of the flood discharge flow based on the flood discharge flow function according to the water level step and the time step includes:
the discharge flow value is calculated according to the following formula,
represents the discharge flow value, +.>Representing the flood discharge flow function, the flood discharge flow function has a definition domain of [ [h 0 ，h 1 ]The value range of n is consistent with [0, ]>]Integer of deltahRepresents the water level step length, deltatThe time step is represented by a time step,h 0 indicating the start-up water level to be adjusted,h 1 representing beach elevation;
calculating the difference value of adjacent water level step-length flood discharge flow according to the following formula to obtain the variation of the flood discharge flow,
representing the variation of the flood discharge flow, wherein n has a value in the range of [0 ],/or ] >]Is an integer of (a).
In an embodiment, the determining the tolerance of convergence of the iterative computation of the flood diversion algorithm according to the variation of the capacity of the flood diversion reservoir and the variation of the flood discharge flow comprises:
and summing the maximum value of the variation of the capacity of the flood regulating reservoir and the maximum value of the variation of the flood discharge flow to obtain the tolerance.
In a second aspect, an embodiment of the present application provides a device for flood control calculation in a tailings pond, where the device includes:
the first determining module is used for determining a flood regulating reservoir capacity function according to the relation between the water level and the flood regulating reservoir capacity;
the second determining module is used for determining a flood discharge flow function according to the relation between the water level and the flood discharge flow;
the acquisition module is used for acquiring the water level step length and the time step length;
the first calculation module is used for calculating the variation of the flood diversion reservoir capacity based on the flood diversion reservoir capacity function according to the water level step length;
the second calculation module is used for calculating the variation of the flood discharge flow based on the flood discharge flow function according to the water level step and the time step;
the third determining module is used for determining the tolerance of iterative computation convergence of the flood regulating algorithm according to the variable quantity of the capacity of the flood regulating reservoir and the variable quantity of the flood discharging flow;
And the fourth determining module is used for carrying out iterative computation of multiple flood diversion operations until the absolute value of the computation result is smaller than the tolerance, stopping the iterative computation, and determining the currently obtained computation result as the operation result of the flood diversion operation of the tailing pond.
In a third aspect, an embodiment of the present application provides an electronic device, including a memory and a processor, where the memory stores a computer program, and the computer program executes the method for flood diversion calculation of a tailings pond according to the first aspect when the processor runs.
In a fourth aspect, an embodiment of the present application provides a computer readable storage medium storing a computer program, which when run on a processor performs the method of the first aspect for floodgate a tailings pond.
The beneficial effects of the embodiment of the application include, for example:
according to the tailing pond flood regulating calculation method provided by the embodiment of the application, the tolerance of the flood regulating calculation loop iteration is obtained based on the variable quantity of the capacity of the flood regulating pond and the variable quantity of the flood discharging flow, and the problem that the flood regulating calculation in the prior art cannot be converged well is effectively solved.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present application, the drawings that are needed in the embodiments will be briefly described below, it being understood that the following drawings only illustrate some embodiments of the present application and therefore should not be considered as limiting the scope, and other related drawings may be obtained according to these drawings without inventive effort for a person skilled in the art.
Fig. 1 is a schematic flow chart of a method for calculating flood control of a tailing pond according to an embodiment of the present application;
FIG. 2 is a diagram illustrating a fold line of the flood control algorithm with different tolerance according to the embodiment of the present application;
FIG. 3 is another schematic diagram of a broken line of a different tolerance flood control algorithm according to an embodiment of the present application;
fig. 4 is a schematic structural diagram of a device for flood control calculation of a tailings pond according to an embodiment of the present application.
Detailed Description
For the purpose of making the objects, technical solutions and advantages of the embodiments of the present application more apparent, the technical solutions of the embodiments of the present application will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present application, and it is apparent that the described embodiments are some embodiments of the present application, but not all embodiments of the present application. The components of the embodiments of the present application generally described and illustrated in the figures herein may be arranged and designed in a wide variety of different configurations.
Thus, the following detailed description of the embodiments of the application, as presented in the figures, is not intended to limit the scope of the application, as claimed, but is merely representative of selected embodiments of the application. All other embodiments, which can be made by those skilled in the art based on the embodiments of the application without making any inventive effort, are intended to be within the scope of the application.
It should be noted that: like reference numerals and letters denote like items in the following figures, and thus once an item is defined in one figure, no further definition or explanation thereof is necessary in the following figures.
In the description of the present invention, it should be noted that, if the terms "upper", "lower", "inner", "outer", and the like indicate an azimuth or a positional relationship based on the azimuth or the positional relationship shown in the drawings, or the azimuth or the positional relationship in which the inventive product is conventionally put in use, it is merely for convenience of describing the present invention and simplifying the description, and it is not indicated or implied that the apparatus or element referred to must have a specific azimuth, be configured and operated in a specific azimuth, and thus it should not be construed as limiting the present invention.
Furthermore, the terms "first," "second," and the like, if any, are used merely for distinguishing between descriptions and not for indicating or implying a relative importance.
It should be noted that the features of the embodiments of the present invention may be combined with each other without conflict.
Example 1
The embodiment of the disclosure provides a tailing pond flood regulating calculation method.
According to the requirements of the safety regulations of tailing ponds (GB 39496-2020) and the design specifications of tailing facilities (GB 50863-2013), the flood regulating algorithm of the tailing ponds should adopt a water balance method. Unit time of The water balance formula of (2) is:, in the formula /> and />Each is a unit time slot->Initial and final time of flood flow (unit m 3 /s）；/> and />Each is a unit time slot->Flood discharge flow (m) at initial and final moments 3 /s）；/> and />Each is a unit time slot->Flood control reservoir capacity (m) at final and initial moments 3 ). Wherein (1)> and />According to the flood discharge flow function->Determination of-> and />According to the stock capacity function->Determining that 2 functions are all corresponding to the water level of the reservoirzAnd (5) correlation. The method for directly solving the water balance formula is easy to generate root increment, needs manual screening, and is easy to add per unit time>The calculation process of shorter flood duration and longer flood duration is realized, and the manual screening workload is extremely high.
Therefore, the operation method commonly adopted at present is to set the water level of the unit stepPerforming loop iterative calculation while setting tolerance +.>When->When the calculation result is considered to be converged. The method is more standard, particularly for computer programming solution, is more common, and replaces manual calculation more and more in the time background of continuous development of the Internet of things and computer technology, so that the calculation efficiency and accuracy are improved.
However, for the current tailings pond flood algorithm using iterative methods, especially when machine calculation is involved, the rule of matching and convergence There is always no explicit method, typically set manually by experience, when the tolerance is + ->When the arrangement is not reasonable, such as +.>Too small, the calculation result diverges and cannot be solved; />Too large results in too early convergence of the calculation, resulting in a larger error in the calculation result.
Based on the above, the embodiment of the application provides a tailing pond flood regulating algorithm, which can better realize the judgment of iterative computation convergence criteria, can realize self-adaptive adjustment according to different pond capacities, different flood discharging system parameters and different iteration step sizes, does not need artificial experience judgment, improves the accuracy and the calculation efficiency, and is suitable for calculation of a tailing pond flood regulating algorithm program.
Specifically, referring to fig. 1, the tailing pond flood control algorithm method includes:
step S110, determining a flood diversion reservoir capacity function according to the relation between the water level and the flood diversion reservoir capacity;
taking flood season flood control calculation of a large tailing pond (second class pond) as an example, the calculation adopts iteration method flood control calculation, and the iteration tolerance determined by the embodiment of the application is compared with the manually set iteration tolerance calculation result. The flood process required for the flood control algorithm is shown in table 1.
Table 1:
in one embodiment, the determining the flood diversion reservoir capacity function according to the water level and the flood diversion reservoir capacity relation comprises:
obtaining a contour topography of a tailing pond; acquiring the area of each contour line between the rising and adjusting water level and the beach top elevation of the contour line topographic map; obtaining an average value of areas of adjacent contour lines, and multiplying the average value of the areas of the adjacent contour lines by a height difference of the adjacent contour lines to obtain a volume of the adjacent contour lines; accumulating the volumes of all the adjacent contour lines between the lifting water level and the beach elevation to obtain the relation between the water level and the flood regulating reservoir capacity; and fitting the water level and the flood regulating reservoir capacity relation to obtain the flood regulating reservoir capacity function.
Taking the specific embodiment as an example, a water level-flood regulating reservoir capacity curve is formed according to the water level and the flood regulating reservoir capacity relation:
firstly, obtaining an actual measurement topographic map of a contour line of a tailing pond, wherein in the embodiment, the actual measurement starting water level of the current flood control algorithm of the tailing pond is 1004.8m, and the beach top elevation is 1007.2m;
according to the actually measured topographic map, the area of each contour line from the lift water level 1004.8m to the beach top elevation 1007.2m is 0.2m, and the area enclosed by the contour lines of each elevation of the tailing pond is shown in table 2.
Table 2:
calculating the average value between the areas of adjacent contour lines;
multiplying the area of the adjacent contour lines by the height difference of the adjacent contour lines to obtain the volume between the adjacent contour lines;
and taking the elevation of the equal-altitude line of the water level 1004.8m as 0 point, and sequentially accumulating the volumes between adjacent equal-altitude lines to obtain the relation between the water level and the flood regulating reservoir capacity. The calculation process of the relation between the water level and the flood diversion reservoir capacity in the embodiment is shown in table 3.
Table 3:
fitting the relation between the water level and the flood regulating reservoir capacity to obtain a function expression of the water level-flood regulating reservoir capacity curve, wherein the embodiment adopts a polynomial of degree 3 for fitting, and the fitting function is as follows:
when the span of the contour line is too large, the contour line and the corresponding area thereof can be processed by adopting an interpolation method.
Step S120, determining a flood discharge flow function according to the relation between the water level and the flood discharge flow;
taking the foregoing specific embodiment as an example, the flood drainage system of this embodiment is in the form of "drainage well+tunnel", and the basic parameters of the flood drainage system are shown in table 4.
Table 4:
in an embodiment, the determining the flood discharge flow function according to the relationship between the water level and the flood discharge flow comprises:
obtaining a flood drainage system arrangement form, wherein the flood drainage system arrangement form comprises at least one of a pit-type, a trench-type and a spillway;
And determining the flood discharge flow function according to at least one of a first preset formula corresponding to the well hole type, a second preset formula corresponding to the groove hole type and a third preset formula corresponding to the spillway and the arrangement form of the flood discharge system.
Flood drainage system arrangements can be divided into the following types: first kind: well-tunnel type, including "drainage well-tunnel", "drainage well-culvert", wherein the framed drainage well is further divided into "frame type" and "window type" according to actual forms. Second kind: slot-tunnel, including "chute-tunnel", "chute-culvert". Third kind: spillway. Fourth kind: combinations of the three forms described above. The relation between the water level and the drainage capacity can be obtained by adopting an empirical formula according to the arrangement form of the actual drainage system of the tailing pond.
The working state of the well-hole drainage system is divided into self-flowing drainage, half-pressure flow and pressure flow according to the size of a drainage head. The well-hole drainage calculations can be found in the following formulas and table 5.
When the working state is frame type free drainage and the water level position is that the water level does not submerge the frame ring beam, the method can calculate according to the following formula:
；
when the water level position is the water level submerged frame ring beam, the water level position can be calculated according to the following formula:
；
When the water level position is the orifice drainage of the water level submerged wellhead, the water level position can be calculated according to the following formula:
；
when the working state is window type free drainage and the water level position is that the water level is between windows, the water level can be calculated according to the following formula:
；
when the working state is window type free drainage and the water level position is the water level at the window part (the window shape is square), the water level can be calculated according to the following formula:
；
when the working state is window type free drainage and the water level position is the water level at the window part (the window shape is Fang Yuanxing), the water level can be calculated according to the following formula:
；
when the operating state is a half pressure flow, it can be calculated according to the following formula:
；
when the operating state is pressure flow, it can be calculated according to the following formula:
；
the calculation parameters may be determined as follows:
H i calculating a water head, m, of a drainage flow representing an ith layer of total submerged working window; h 0 Calculating a water head, m, of a drainage flow representing an unsubmerged working window at the uppermost layer; h represents a calculated water head, and is the difference between the reservoir water level and the central elevation of the tunnel inlet section, and m; h z Representing calculated head as reservoir water level and tunnel downstreamThe difference between the elevation of the center of the outlet section and the elevation of the downstream water level is m when water exists in the downstream water; h y Represents the overflow weir drainage head, m; h j Representing a wellhead drainage head, m;representing the area of a drainage window, m; />Representing the reduced cross-sectional area of the water flow at the wellhead, wherein m is; />Representing the total water clearance area, m, between the frame upright posts and the ring beams; omega represents the total area of the window in the water depth range in the well, m; />Representing the cross-sectional area of a drainage well shaft, m; />Representing the reduced cross-sectional area of the tunnel (culvert) inlet flow, m,/or->；/>Represents the cross-sectional area of the tunnel (culvert) inlet, m < th >; />Representing the area of the section of the downstream outlet of the tunnel (culvert), m; />Representing the section area of a tunnel (culvert pipe) calculated hole section, m; zeta represents the local head loss coefficient along the tunnel, including corner, bifurcation, section change equalization, searching on page 331 of tailing facility design reference materials; />Representing coefficients, by tailings facilitiesPage 335 of the design reference; />The local head loss coefficient of the tunnel inlet is represented, the right angle inlet is 0.5, the round angle or oblique angle inlet is 0.2-0.25, and the bell mouth inlet is 0.1-0.2;the local head loss coefficient of the water flow steering in the drainage well is shown and checked by page 335 of tailing facility design reference materials; Representing the local head loss coefficient of the inlet of the drainage well, and searching by page 329 of tailing facility design reference materials; />The local head loss coefficient of the frame is represented as the sum of the local head loss coefficients of the upright post and the cross beam, namely: />The method comprises the steps of carrying out a first treatment on the surface of the Beta represents the shape coefficients of a beam and a column, the rectangular section is 2.42, and the circular section is 1.79; />Representing the effective section coefficients of beams and columns, the ratio of the clearance space to the center space can be ∈>The tailings facility design reference materials are searched on page 329; />Represents the lateral contraction coefficient, ++>；/>Representing the sudden reduction of area coefficient, and searching by page 335 of tailing facility design reference materials; />Represents the inner diameter of the drainage well, m; />Representing the inner diameter of the calculated section of the tunnel (culvert), m, for non-circular holes +.>；/>Representing the water depth, m, above the top of a tunnel (culvert pipe) in the drainage well; />Representing the length of a tunnel (culvert pipe) calculation tunnel section, m; />The hydraulic radius of a pipe section is calculated by a tunnel (culvert pipe), and m is represented; />The hydraulic radius of the section of the well bore of the drainage well is represented by m; />Representing a weir flow coefficient; />Calculated by thin wall weir, < >>；/>Calculated according to the practical weir,；/>representing the width of the weir crest, m; />Represents the width of one drain opening, m; / >Representing the number of water outlets on the same cross section; />Represents the water head loss coefficient along the drainage well, +.>；/>Represents the along-path head loss coefficient of the tunnel (culvert)>；/>Representing a thank you coefficient; n represents a hole wall roughness coefficient;
wherein ,；/>；/>；/>；/>；/>；/>；/>；/>。
the working state of the tank-hole drainage system is divided into self-flowing drainage, half-pressure drainage and pressure drainage according to the size of a drainage head. The slot-hole bleed calculation can be found in the following formulas and table 6.
When the working state is free drainage and the water level does not exceed the highest point of the upper edge of the arch plate, the water level can be calculated according to the following formula:
；
when the water level is the highest point of the upper edge of the arch plate, the water level position can be calculated according to the following formula:
；
wherein ,；
when the operating state is a half pressure flow, it can be calculated according to the following formula:
；
when the operating state is pressure flow, it can be calculated according to the following formula:
；
the calculation parameters may be determined as follows:
a free drainage head is expressed, and the lowest point of the water passing part of the side wall of the chute is counted by (m); />A free drainage head is represented, and (m) is calculated from the highest point of the upper edge of the arch plate; />Represents a half pressure flow bleed head (m); />Represents the pressure flow bleed head (m);brepresents a weir bottom width (m); />Representing the inclination angle of the chute; m is m 1 Representing a weir flow coefficient; m is m 2 Representing orifice flow coefficient; />Representing a flooding coefficient; />Representing the cross-sectional area (m) of the chute 2 ）；/>Represents the outlet cross-sectional area (m) of the tunnel (culvert pipe) 2 ）。
Step S130, acquiring a water level step and a time step;
in one embodiment, the water level step is divided into a range of 0.001-0.000001 mThe time step is any value between 3 and 60 minutes. In this embodiment, the water level step may be set to 0.0001m and the time step to 15min.
Step S140, calculating the variation of the flood diversion reservoir capacity based on the flood diversion reservoir capacity function according to the water level step;
in an embodiment, the calculating the variation of the capacity of the flood diversion reservoir based on the capacity function of the flood diversion reservoir according to the water level step comprises:
the flood regulating reservoir capacity value is calculated according to the following formula,
representing the value of the flood regulating reservoir->Representing the flood regulating reservoir capacity function, wherein the definition domain of the flood regulating reservoir capacity function is [ [h 0 ，h 1 ]The value range of n is consistent with [0, ]>]Integer of deltahRepresents the step size of the water level,h 0 indicating the start-up water level to be adjusted,h 1 representing beach elevation;
calculating the difference value of the adjacent water level step-length flood regulating reservoir capacity according to the following formula to obtain the variation quantity of the flood regulating reservoir capacity,
indicating the variation of the capacity of the flood regulating reservoir, wherein the value range of n is in accordance with [0, ], and- >]Is an integer of (a).
Step S150, calculating the variation of the flood discharge flow based on the flood discharge flow function according to the water level step and the time step;
in an embodiment, the calculating the variation of the flood discharge flow based on the flood discharge flow function according to the water level step and the time step includes:
the discharge flow value is calculated according to the following formula,
represents the discharge flow value, +.>Representing the flood discharge flow function, the flood discharge flow function has a definition domain of [ [h 0 ，h 1 ]The value range of n is consistent with [0, ]>]Integer of deltahRepresents the water level step length, deltatThe time step is represented by a time step,h 0 indicating the start-up water level to be adjusted,h 1 representing beach elevation;
calculating the difference value of adjacent water level step-length flood discharge flow according to the following formula to obtain the variation of the flood discharge flow,
representing the variation of the flood discharge flow, wherein n has a value in the range of [0 ],/or ]>]Is an integer of (a).
Step S160, determining the tolerance of iterative computation convergence of the flood regulating algorithm according to the variable quantity of the capacity of the flood regulating reservoir and the variable quantity of the flood discharging flow;
according to the specific embodiment, the variation of the capacity of the flood regulating reservoir and the variation of the flood discharging amount are calculated:
in the embodiment, the water level step is 0.0001, the time step is 15min, the water level is 1004.8m, and the beach height is 1007.2m;
And fitting to obtain a water level-flood regulating reservoir capacity curve. In this embodiment, the functional expression of the water level-flood regulating reservoir capacity curve is the formula:
and fitting the obtained water level-flood discharge flow curve by adopting an empirical formula or a hydraulic test result. In this embodiment, an empirical formula method is adopted, and according to the calculation, in the calculation interval from the water level 1004.8m to the beach top 1007.2m, the change of the drainage flow state is changed from free flow to orifice flow, so that the water level-flood discharge flow curve expression is divided into two sections, and the water level-flood discharge flow curve expression is as follows:
the values of the flood control reservoir capacity corresponding to all the water level step sizes on the definition domain are calculated, the initial water level of the embodiment is 1004.8m, the beach top is 1007.2m, and the unit step size is 0.0001m, and then the values of the flood control reservoir capacity corresponding to all the water level step sizes on the definition domain are shown in table 5.
Table 5:
and calculating the flood discharge amount corresponding to all the water level step sizes on the definition domain. In this embodiment, the initial water level is 1004.8m, the beach top is 1007.2m, the time step is set to 15min, the drainage volumes corresponding to all the water level steps in the defined area are calculated according to the water level-drainage flow curve expression, and the drainage volumes corresponding to all the water level steps in the defined area are obtained by multiplying the drainage volumes by the time, as shown in table 6.
Table 6:
and calculating the difference value of the adjacent water level step size flood regulating reservoir capacity. The difference between adjacent water level step sizes in the present embodiment is shown in table 7.
Table 7:
the difference between the adjacent water level step-size flood discharge amounts is calculated, and the difference between the adjacent water level step-size flood discharge amounts in this embodiment is shown in table 8.
Table 8:
in one embodiment, the tolerance is obtained by summing the maximum value of the variation of the capacity of the flood bank and the maximum value of the variation of the flood discharge flow.
Finally, determining the minimum tolerance of iterative computation convergence of the tailing pond flood control algorithm, in this embodiment, as can be seen from table 10, the maximum value of the pond tolerance between 1007.1999m and 1007.2m is 14.8551, the maximum value of the flood discharge difference between 1005.4916m and 1005.4917m is 5.1117, so that the minimum iterative tolerance under the conditions of water level step of 0.0001m and time step of 15min is。/>
And S170, performing iterative computation of multiple flood diversion algorithms until the absolute value of the computation result is smaller than the tolerance, stopping the iterative computation, and determining the currently obtained computation result as the computation result of the flood diversion algorithm of the tailing pond.
In an embodiment, the performing the iterative computation of the flood diversion algorithm for a plurality of times until the absolute value of the computation result is smaller than the tolerance, and stopping the iterative computation includes:
The iterative calculation of the flood diversion calculation is performed according to the following formula,
wherein ,ΔTthe time period of the preset time period is indicated,Q s andQ z flood-entering respectively representing initial time and final time of preset time periodThe flow rate of the liquid is controlled,q s andq z respectively representing the flood discharge flow at the initial moment and the final moment of the preset time period,V s andV z respectively representing the flood regulating reservoir capacity at the initial moment and the final moment of a preset time period;q s andq z determined according to the flood discharge flow function,V s andV z determining according to the flood regulating reservoir capacity function; when (when)V s -V z And stopping the iterative calculation when the absolute value of (a) is smaller than the tolerance.
In order to verify the implementation effect of the scheme, the tolerance 19.9668 calculated by the embodiment is compared with the flood regulating calculation results of other tolerance values, the flood regulating calculation adopts an iterative calculation method, and the results are shown in table 9:
table 9:
as can be seen from table 9, when the tolerance is set to 10, the highest water level and the total residual value cannot be obtained because the iterative calculation is not converged; when the residual values are different, the obtained highest water level and the total residual values are different, and when the residual value is 19.9668 in this embodiment, the total residual value is 2163.8503, which is the smallest in all the embodiments, so that it is proved that the water level 1005.59439m calculated from the tolerance is the closest to the actual one. The different tolerance flood control algorithm is shown in figure 2. Wherein 210 represents tolerance 19.9668;220 represents a tolerance of 50;230 denotes a tolerance of 100;240 denotes a tolerance 300;250 represents a tolerance 500.
According to the tailing pond flood regulating calculation method provided by the embodiment of the application, the tolerance of the flood regulating calculation loop iteration is obtained based on the variable quantity of the capacity of the flood regulating pond and the variable quantity of the flood discharging flow, and the problem that the flood regulating calculation in the prior art cannot be converged well is effectively solved. Meanwhile, the self-adaptive adjustment can be realized according to different reservoir capacities, different flood drainage system parameters and different iteration step sizes, artificial experience judgment is not needed, the accuracy is improved, and meanwhile, the calculation efficiency is improved, so that the method is suitable for programming calculation of the flood control algorithm of the tailing pond.
Example 2
In order to verify the self-adaptive characteristics of the method, the water level step and the time step in the embodiment 1 are adjusted, and the flood process, the reservoir capacity and the flood drainage system of the embodiment are the same as those in the embodiment 1. Therefore, the first 2 steps of this embodiment are the same as those of embodiment 1, and thus will not be described again.
Step S210: step S110 is the same as in example 1;
step S220: step S120 is the same as in example 1;
step S230: setting the water level step and the time step of iterative calculation. The value range is between 0.001 and 0.000001m and can be dividedThe time step is any value between 3 and 60 minutes. In this embodiment, the water level step may be set to 0.001m and the time step to 30min.
Step S240: calculating the variation of the capacity of the flood regulating reservoir and the variation of the flood discharging capacity:
in the embodiment, the water level step is 0.001, the time step is 30min, the water level is 1004.8m, and the beach height is 1007.2m;
and fitting to obtain a water level-flood regulating reservoir capacity curve. In this embodiment, the functional expression of the water level-flood regulating reservoir capacity curve is:
and fitting the obtained water level-flood discharge flow curve by adopting an empirical formula or a hydraulic test result. In this embodiment, an empirical formula method is adopted, and according to calculation, in a calculation interval from the water level 1004.8m to the beach top 1007.2m, the change of the drainage flow state is changed from free flow to orifice flow, so that the water level-flood discharge flow curve expression is divided into two sections, and the water level-flood discharge flow curve expression is:
calculating the values of the flood bank capacity corresponding to all the water level step sizes on the definition domain, wherein the initial water level of the embodiment is 1004.8m, the beach top is 1007.2m, and the unit step size is 0.001m, and the values of the flood bank capacity corresponding to all the water level step sizes on the definition domain are shown in table 10:
table 10:
and calculating the flood discharge amount corresponding to all the water level step sizes on the definition domain. In this embodiment, the initial water level is 1004.8m, the beach top is 1007.2m, the time step is set to 30min, the drainage volumes corresponding to all the water level steps on the definition domain are calculated according to the water level-drainage flow curve expression, and the drainage volumes corresponding to all the water level steps on the definition domain are obtained by multiplying the drainage volumes by the time, as shown in table 11.
Table 11:
and calculating the difference value of the adjacent water level step size flood regulating reservoir capacity. The difference between adjacent water level step sizes in the present embodiment is shown in table 12:
table 12:
the difference between the adjacent water level step-size flood discharge amounts is calculated, and in this embodiment, the difference between the adjacent water level step-size flood discharge amounts is shown in table 13:
table 13:
step S250: and determining the minimum tolerance of iterative computation convergence of the tailing pond flood diversion algorithm. In this example, as can be seen from tables 12 and 13, the reservoir capacity between 1007.1999m and 1007.2m is at maximum 148.5461, and the flood discharge difference between 1005.490m and 1005.491m is at maximum 102.0519, therefore, the minimum iterative tolerance of this embodiment under the conditions of a water level step of 0.001m and a time step of 30min is。
Effect verification of the embodiment:
in order to verify the implementation effect of the scheme, the tolerance 250.598 calculated by the embodiment is compared with the flood regulating calculation results of other tolerance values, the flood regulating calculation adopts an iterative calculation method, and the results are shown in table 14:
table 14:
as can be seen from table 14, when the tolerance is set to 150, the highest water level and total residual value cannot be obtained because the iterative calculation is not converged; when the residual values are different, the obtained highest water level and the total residual values are different, and when the residual value is 250.598 in this embodiment, the total residual value is 12224.8696, which is the smallest in all the embodiments, so that it is proved that the water level 1005.941m calculated from the tolerance is the closest to the actual one. The different tolerance flood control algorithm is shown in figure 3. Wherein 310 represents tolerance 250.598;320 represents a tolerance 500;330 denotes a tolerance 1000;340 represents a tolerance 3000;350 denotes a tolerance 3900.
It can be seen from embodiments 1 and 2 that the tolerance value determined by the tailing pond flood regulating algorithm provided by the embodiment of the application can be ensured to be relatively reasonable, and the total residual error is relatively small while the absolute convergence of the flood regulating algorithm can be ensured, so that the accuracy of the flood regulating algorithm is improved.
Example 3
In addition, the embodiment of the disclosure provides a tailing pond flood regulating calculation device.
Specifically, as shown in fig. 4, the tailing pond flood control calculation apparatus 400 includes:
a first determining module 410, configured to determine a flood diversion reservoir capacity function according to a relationship between a water level and a flood diversion reservoir capacity;
a second determining module 420 for determining a flood discharge flow function according to a relationship between the water level and the flood discharge flow;
an obtaining module 430, configured to obtain a water level step and a time step;
a first calculation module 440, configured to calculate a variation of the capacity of the flood diversion reservoir based on the capacity function of the flood diversion reservoir according to the water level step;
a second calculating module 450, configured to calculate a variation of the flood discharge flow based on the flood discharge flow function according to the water level step and the time step;
a third determining module 460, configured to determine a tolerance of iterative computation convergence of the flood diversion algorithm according to the variation of the capacity of the flood diversion reservoir and the variation of the flood drainage flow;
And a fourth determining module 470, configured to perform iterative computation of the flood diversion computation for multiple times until the absolute value of the computation result is smaller than the tolerance, stop the iterative computation, and determine the currently obtained computation result as the computation result of the flood diversion computation of the tailing pond.
The tailing pond flood control calculation device 400 provided in this embodiment can implement the tailing pond flood control calculation method provided in embodiment 1, and in order to avoid repetition, details are not repeated here.
According to the tailing pond flood regulating calculation device provided by the embodiment of the application, the tolerance of the flood regulating calculation loop iteration is obtained based on the variable quantity of the capacity of the flood regulating pond and the variable quantity of the flood discharging flow, and the problem that the flood regulating calculation in the prior art cannot be converged well is effectively solved. Meanwhile, the self-adaptive adjustment can be realized according to different reservoir capacities, different flood drainage system parameters and different iteration step sizes, artificial experience judgment is not needed, the accuracy is improved, and meanwhile, the calculation efficiency is improved, so that the method is suitable for programming calculation of the flood control algorithm of the tailing pond.
Example 4
Furthermore, an embodiment of the present disclosure provides an electronic device, including a memory and a processor, where the memory stores a computer program that, when run on the processor, performs the tailings pond flood control calculation method provided in embodiment 1.
The electronic device provided by the embodiment of the application can realize the tailing pond flood control calculation method provided by the embodiment 1, and in order to avoid repetition, the description is omitted.
The electronic equipment provided by the embodiment obtains the tolerance of the flood regulating calculation loop iteration based on the variable quantity of the flood regulating reservoir capacity and the variable quantity of the flood discharging flow, and effectively solves the problem that the flood regulating calculation in the prior art can not be converged well. Meanwhile, the self-adaptive adjustment can be realized according to different reservoir capacities, different flood drainage system parameters and different iteration step sizes, artificial experience judgment is not needed, the accuracy is improved, and meanwhile, the calculation efficiency is improved, so that the method is suitable for programming calculation of the flood control algorithm of the tailing pond.
Example 5
The application also provides a computer readable storage medium, wherein the computer readable storage medium stores a computer program, and the computer program realizes the tailing pond flood diversion calculation method provided by the embodiment 1 when being executed by a processor.
In the present embodiment, the computer readable storage medium may be a Read-Only Memory (ROM), a random access Memory (Random Access Memory, RAM), a magnetic disk, an optical disk, or the like.
The computer readable storage medium provided in this embodiment may implement the tailing pond flood control algorithm provided in embodiment 1, and in order to avoid repetition, details are not described here.
It should be noted that, in this document, the terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or terminal that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, article, or terminal. Without further limitation, an element defined by the phrase "comprising one … …" does not exclude the presence of other like elements in a process, method, article or terminal comprising the element.
From the above description of the embodiments, it will be clear to those skilled in the art that the above-described embodiment method may be implemented by means of software plus a necessary general hardware platform, but of course may also be implemented by means of hardware, but in many cases the former is a preferred embodiment. Based on such understanding, the technical solution of the present application may be embodied essentially or in a part contributing to the prior art in the form of a software product stored in a storage medium (e.g. ROM/RAM, magnetic disk, optical disk) comprising instructions for causing a terminal (which may be a mobile phone, a computer, a server, an air conditioner, or a network device, etc.) to perform the method according to the embodiments of the present application.
The foregoing is merely illustrative of the present invention, and the present invention is not limited thereto, and any changes or substitutions easily contemplated by those skilled in the art within the scope of the present invention should be included in the present invention. Therefore, the protection scope of the present invention shall be subject to the protection scope of the claims.

Claims (8)

1. A method for flood control algorithm in a tailings pond, the method comprising:
determining a flood regulating reservoir capacity function according to the relation between the water level and the flood regulating reservoir capacity;
determining a flood discharge flow function according to the relationship between the water level and the flood discharge flow;
acquiring a water level step and a time step;
calculating the variable quantity of the flood regulating reservoir capacity based on the flood regulating reservoir capacity function according to the water level step length;
calculating the variation of the flood discharge flow based on the flood discharge flow function according to the water level step and the time step;
determining the tolerance of iterative computation convergence of the flood regulating algorithm according to the variable quantity of the flood regulating reservoir capacity and the variable quantity of the flood discharging flow;
performing iterative computation of flood diversion calculation for a plurality of times until the absolute value of the computation result is smaller than the tolerance, stopping the iterative computation, and determining the currently obtained computation result as the calculation result of the flood diversion calculation of the tailing pond;
The calculating the variable quantity of the flood diversion reservoir capacity based on the flood diversion reservoir capacity function according to the water level step length comprises the following steps:
the flood regulating reservoir capacity value is calculated according to the following formula,
representing the value of the flood regulating reservoir->Representing the flood regulating reservoir capacity function, wherein the definition domain of the flood regulating reservoir capacity function is [ [h 0 ，h 1 ]The value range of n is consistent with [0, ]>]Integer of deltahRepresents the step size of the water level,h 0 indicating the start-up water level to be adjusted,h 1 representing beach elevation;
calculating the difference value of the adjacent water level step-length flood regulating reservoir capacity according to the following formula to obtain the variation quantity of the flood regulating reservoir capacity,
indicating the variation of the capacity of the flood regulating reservoir, wherein the value range of n is in accordance with [0, ], and->]Is an integer of (2);
the calculating the variable quantity of the flood discharge flow based on the flood discharge flow function according to the water level step length and the time step length comprises the following steps:
the discharge flow value is calculated according to the following formula,
represents the discharge flow value, +.>Representing the flood discharge flow function, the flood discharge flow function has a definition domain of [ [h 0 ，h 1 ]The value range of n is consistent with [0, ]>]Integer of deltahRepresents the water level step length, deltatThe time step is represented by a time step,h 0 indicating the start-up water level to be adjusted,h 1 representing beach elevation;
calculating the difference value of adjacent water level step-length flood discharge flow according to the following formula to obtain the variation of the flood discharge flow,
Representing the variation of the flood discharge flow, wherein n has a value in the range of [0 ],/or ]>]Is an integer of (a).
2. The method of claim 1, wherein performing the iterative computation of the plurality of flood diversion algorithms until the absolute value of the computation result is less than the tolerance, and stopping the iterative computation comprises:
the iterative calculation of the flood diversion calculation is performed according to the following formula,
wherein ,ΔTthe time period of the preset time period is indicated,Q s andQ z hong Liuliang indicating the initial and final moments of the preset time period respectively,q s andq z respectively representing the flood discharge flow at the initial moment and the final moment of the preset time period,V s andV z respectively representing the flood regulating reservoir capacity at the initial moment and the final moment of a preset time period;q s andq z determined according to the flood discharge flow function,V s andV z determining according to the flood regulating reservoir capacity function;
when (when)V s -V z And stopping the iterative calculation when the absolute value of (a) is smaller than the tolerance.
3. The method of claim 1, wherein determining a flood bank capacity function based on the water level and the flood bank capacity relationship comprises:
obtaining a contour topography of a tailing pond;
acquiring the area of each contour line between the rising and adjusting water level and the beach top elevation of the contour line topographic map;
obtaining an average value of areas of adjacent contour lines, and multiplying the average value of the areas of the adjacent contour lines by a height difference of the adjacent contour lines to obtain a volume of the adjacent contour lines;
Accumulating the volumes of all the adjacent contour lines between the lifting water level and the beach elevation to obtain the relation between the water level and the flood regulating reservoir capacity;
and fitting the water level and the flood regulating reservoir capacity relation to obtain the flood regulating reservoir capacity function.
4. The method of claim 1, wherein determining a flood discharge flow function based on a relationship between water level and flood discharge flow comprises:
obtaining a flood drainage system arrangement form, wherein the flood drainage system arrangement form comprises at least one of a pit-type, a trench-type and a spillway;
and determining the flood discharge flow function according to at least one of a first preset formula corresponding to the well hole type, a second preset formula corresponding to the groove hole type and a third preset formula corresponding to the spillway and the arrangement form of the flood discharge system.
5. The method of claim 1, wherein determining a tolerance for convergence of iterative calculation of a flood control algorithm based on the amount of change in the flood control reservoir capacity and the amount of change in the flood discharge flow comprises:
and summing the maximum value of the variation of the capacity of the flood regulating reservoir and the maximum value of the variation of the flood discharge flow to obtain the tolerance.
6. A tailings pond flood regulating algorithm apparatus, the apparatus comprising:
the first determining module is used for determining a flood regulating reservoir capacity function according to the relation between the water level and the flood regulating reservoir capacity;
the second determining module is used for determining a flood discharge flow function according to the relation between the water level and the flood discharge flow;
the acquisition module is used for acquiring the water level step length and the time step length;
the first calculation module is used for calculating the variation of the flood diversion reservoir capacity based on the flood diversion reservoir capacity function according to the water level step length;
the second calculation module is used for calculating the variation of the flood discharge flow based on the flood discharge flow function according to the water level step and the time step;
the third determining module is used for determining the tolerance of iterative computation convergence of the flood regulating algorithm according to the variable quantity of the capacity of the flood regulating reservoir and the variable quantity of the flood discharging flow;
the fourth determining module is used for carrying out iterative computation of flood diversion calculation for a plurality of times until the absolute value of the computation result is smaller than the tolerance, stopping the iterative computation, and determining the currently obtained computation result as the calculation result of the flood diversion calculation of the tailing pond;
the first computing module is further configured to:
the flood regulating reservoir capacity value is calculated according to the following formula,
Representing the value of the flood regulating reservoir->Representing the flood regulating reservoir capacity function, wherein the definition domain of the flood regulating reservoir capacity function is [ [h 0 ，h 1 ]The value range of n is consistent with [0, ]>]Integer of deltahRepresents the step size of the water level,h 0 indicating the start-up water level to be adjusted,h 1 representing beach elevation;
calculating the difference value of the adjacent water level step-length flood regulating reservoir capacity according to the following formula to obtain the variation quantity of the flood regulating reservoir capacity,
indicating the variation of the capacity of the flood regulating reservoir, wherein the value range of n is in accordance with [0, ], and->]Is an integer of (2);
the first computing module is further configured to:
the discharge flow value is calculated according to the following formula,
represents the discharge flow value, +.>Representing the flood discharge flow function, the flood discharge flow function has a definition domain of [ [h 0 ，h 1 ]The value range of n is consistent with [0, ]>]Integer of deltahRepresents the water level step length, deltatThe time step is represented by a time step,h 0 indicating the start-up water level to be adjusted,h 1 representing beach elevation;
calculating the difference value of adjacent water level step-length flood discharge flow according to the following formula to obtain the variation of the flood discharge flow,
representing the variation of the flood discharge flow, wherein n has a value in the range of [0 ],/or ]>]Is an integer of (a).
7. An electronic device comprising a memory and a processor, the memory storing a computer program that, when run by the processor, performs the tailings pond flood control algorithm of any one of claims 1 to 5.
8. A computer readable storage medium, characterized in that it stores a computer program which, when run on a processor, performs the tailings pond flood control algorithm of any one of claims 1 to 5.
CN202310890320.8A 2023-07-20 2023-07-20 Tailing pond flood regulating calculation method and device, electronic equipment and storage medium Active CN116628402B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310890320.8A CN116628402B (en) 2023-07-20 2023-07-20 Tailing pond flood regulating calculation method and device, electronic equipment and storage medium

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310890320.8A CN116628402B (en) 2023-07-20 2023-07-20 Tailing pond flood regulating calculation method and device, electronic equipment and storage medium

Publications (2)

Publication Number Publication Date
CN116628402A CN116628402A (en) 2023-08-22
CN116628402B true CN116628402B (en) 2023-09-29

Family

ID=87638544

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310890320.8A Active CN116628402B (en) 2023-07-20 2023-07-20 Tailing pond flood regulating calculation method and device, electronic equipment and storage medium

Country Status (1)

Country Link
CN (1) CN116628402B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117151309B (en) * 2023-10-30 2024-01-19 矿冶科技集团有限公司 Method and device for predicting flood regulating reservoir capacity of tailing reservoir, electronic equipment and medium

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106228242A (en) * 2016-07-28 2016-12-14 三峡大学 A kind of high accuracy flood routing method based on optimized algorithm
CN111275242A (en) * 2020-01-08 2020-06-12 昆明理工大学 Downstream protection tailing pond flood regulating method based on tabu search optimization algorithm
CN112101697A (en) * 2020-07-17 2020-12-18 中国安全生产科学研究院 Tailings pond flood regulating calculation system flow manager
CN113610693A (en) * 2021-10-09 2021-11-05 矿冶科技集团有限公司 Method and device for flood discharge of tailing pond

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20140278561A1 (en) * 2013-03-14 2014-09-18 Stoneriver National Flood Services, Inc. Computerized system and method for determining flood risk

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106228242A (en) * 2016-07-28 2016-12-14 三峡大学 A kind of high accuracy flood routing method based on optimized algorithm
CN111275242A (en) * 2020-01-08 2020-06-12 昆明理工大学 Downstream protection tailing pond flood regulating method based on tabu search optimization algorithm
CN112101697A (en) * 2020-07-17 2020-12-18 中国安全生产科学研究院 Tailings pond flood regulating calculation system flow manager
CN113610693A (en) * 2021-10-09 2021-11-05 矿冶科技集团有限公司 Method and device for flood discharge of tailing pond

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
尾矿库调洪演算的时变分析模型及应用;邓红卫;陈宜楷;雷涛;;中国安全科学学报(第10期);全文 *
尾矿库调洪过程分析与迭代演算;吴鹏程;郭修建;;现代矿业(第06期);全文 *

Also Published As

Publication number Publication date
CN116628402A (en) 2023-08-22

Similar Documents

Publication Publication Date Title
CN116628402B (en) Tailing pond flood regulating calculation method and device, electronic equipment and storage medium
CN110851765B (en) Method for acquiring river course roughness in gradual flow state
CN112116785B (en) Tailing pond disaster early warning method and device based on strong rainfall weather forecast
CN112270115B (en) Cellular automaton-based complex terrain flood inundation process simulation method
CN111046563B (en) Method for simulating continuous burst flood of cascade reservoir group
CN112464584A (en) Method for estimating water level and flow of free surface flow
CN105783994A (en) Online monitoring and rapid early warning method for tailing dam
CN113486556A (en) Improved efficient automatic history fitting method for oil and gas reservoir
CN110059443B (en) Rapid prediction method for water taking and discharging temperature of stratified reservoir
CN109469465B (en) Method and device for processing parameters of water injection well and storage medium
CN112328952B (en) Drainage well discharge flow calculation method and device and electronic equipment
Zhang et al. Discharge coefficient of a spillway with a riser perforated by rectangular orifices
CN111680834B (en) Method for calculating precipitation stabilization time of submerged well
CN115481353A (en) River channel comprehensive roughness calculation method based on beach trough partition
CN114896909A (en) Open channel flow calculation method based on water level height
CN115712800A (en) Reservoir water level fluctuation processing method
CN107015289B (en) Trap evaluation stock number determines method and device
CN104713920A (en) Method and system for uniformly and quantitatively expressing soil water and surface water
CN114842921A (en) Municipal sewage dephosphorization agent adding method based on linear regression and decision tree model
CN114970392A (en) Fracture fractal characterization method and system and readable storage medium
Munta et al. Study of the inception length of flow over stepped spillway models
Wang et al. Topographic evolution around twin piers in a tandem arrangement
CN113128039A (en) Vertical seam type fishway step-shaped water surface line calculating method and device and storage medium
CN111798108B (en) Urban drainage area cooperative scheduling method
Maturi et al. Experimental Study of the Effects of Flow Discharge, Diameter, and Depth on Shear Stress in a Rectangular Channel with Rigid Unsubmerged Vegetation

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