Method and device for analyzing injection rate interval of curtain grouting of upper reservoir dam toe of pumping and storage project

Abstract

The invention relates to an injection rate interval analysis method and device for curtain grouting of an upper reservoir dam toe of a pumping and storage project, and the method comprises the steps: adding grouting data into a historical data set; the type of the injection rate duration process line of each injection section is detected; the fluctuation section unit average injection rate of the standard injection section is calculated within the fluctuation section range; calculating grouting similar indexes of the standard grouting section; the standard grouting sections meeting the preset same grouting similar mode conditions are clustered; outliers in the parameter array are removed; and according to the normal distribution, calculating a confidence interval of the unit average injection rate of the fluctuation section when the confidence coefficient of the similar grouting section set meets a preset percentage condition, so as to determine an injection rate interval when the upper reservoir dam toe curtain grouting of the pumping and storage project is in a grouting similar mode. Therefore, the problem that the quality of a finally formed curtain is affected due to the fact that the subjectivity and the empirical property of curtain grouting construction operation are too high because historical data are difficult to feed back construction in the related technology is solved.

Description

Technical Field

[0001] The present application relates to the technical field of grouting in water conservancy projects, and in particular to an injection rate interval analysis method and device for curtain grouting at the toe of a reservoir dam in a pumped storage project. Background Art

[0002] The construction of pumped storage power stations is an indispensable part of my country's strategy to achieve energy revolution, build a new energy system and ensure energy security. It has many components, among which the upper reservoir is the main water storage part of the project. It is located in a high-altitude position and is used to store water pumped up from the lower reservoir to maintain sufficient water volume for power generation. In order to enhance the foundation stability and bearing capacity of the upper reservoir of the pumped storage project and improve the safety of the project; at the same time, in order to reduce leakage, maintain the water volume of the reservoir, and prevent the formation of leakage channels that reduce the economic efficiency of the project, it is very necessary to perform curtain grouting on the bedrock of the upper reservoir.

[0003] Curtain grouting is a common construction technology in water conservancy projects, mainly used to reinforce the foundation, prevent seepage, and plug leaks to form underground anti-seepage walls. The construction process of curtain grouting projects mainly includes drilling, flushing, pressure testing, segmented grouting, and sealing. In order to ensure that the grouting material fully fills the holes and removes air and moisture in the holes, the grouting speed and pressure need to be controlled during the grouting process. In related technologies, the construction parameters of curtain grouting can be monitored in real time through an automated grouting data recorder.

[0004] However, the grouting construction party in the relevant technology can only confirm the reasonable range of injection rate based on the on-site water pressure test before the formal grouting and rely on the personal experience of technical personnel, resulting in rough pressure regulation during the grouting process. The grouting construction is empirical and vague, making it difficult to guarantee the grouting quality, which urgently needs to be improved. Summary of the invention

[0005] The present application provides an injection rate interval analysis method and device for curtain grouting at the toe of a reservoir dam in a pumped storage project, so as to solve the problem that the grouting construction party in the related technology can only confirm the reasonable injection rate interval based on the on-site water pressure test before the formal grouting and rely on the personal experience of the technicians, resulting in rough pressure regulation during the grouting process, empirical and vague grouting construction, and difficulty in ensuring the grouting quality.

[0006] The first aspect of the present application provides an injection rate interval analysis of the curtain grouting of the toe of the reservoir dam in a pumped storage project, comprising the following steps: acquiring the grouting data of the grouting site, and adding the grouting data to a historical data set; based on the injection rate time data of each grouting section in the historical data set, detecting the type of the injection rate time process line of each grouting section, and screening out the grouting sections that meet the preset normal conditions from the types of the injection rate time process lines, and obtaining a standard grouting section based on the grouting sections that meet the preset normal conditions; based on the injection rate inflection point time of the standard grouting section, determining the fluctuation section range of the injection rate time process line of the standard grouting section, and calculating the fluctuation section unit average injection rate of the standard grouting section within the fluctuation section range; based on the grouting parameters of the standard grouting section, calculating the grouting similarity index of the standard grouting section; obtaining the grouting similarity index of the standard grouting section according to the grouting similarity index. The method comprises the following steps: determining a grouting similarity mode, clustering the standard grouting sections that meet the preset same grouting similarity mode conditions, and obtaining at least one similar grouting section set; transmitting the unit average injection rate of the fluctuation section in the at least one similar grouting section set to the parameter array of the corresponding similar grouting section set, and removing the outliers in the parameter array to obtain the parameter array after removing the outliers; performing a normality test on the parameter array after removing the outliers to obtain a normal test result, and judging whether the normal test result meets the preset qualified conditions; when the normal test result meets the preset qualified conditions, calculating the confidence interval of the unit average injection rate of the fluctuation section of the similar grouting section set when the confidence level meets the preset percentage condition according to the normal distribution, and determining the injection rate interval when the toe curtain grouting of the upper reservoir dam of the pumped storage project is in the grouting similarity mode according to the confidence interval.

[0007] Optionally, in one embodiment of the present application, determining the type of the injection rate time process line of each irrigation section, and selecting the irrigation sections that meet the preset normal conditions from the types of the injection rate time process lines, includes: determining the occurrence position of the injection rate inflection point of each irrigation section based on the continuous change rule between the injection rate and time of each irrigation section; determining whether the injection rate meets the preset specification conditions according to the occurrence position; and when the injection rate meets the preset specification conditions, classifying the line shape of the injection rate time process line to obtain the irrigation sections that meet the preset normal conditions.

[0008] Optionally, in one embodiment of the present application, the calculation formula for the unit average injection rate of the fluctuation section of the standard injection section is:

[0009]

[0010] Wherein, L is the length of the standard injection section; Q is the injection rate time series array; i is the injection rate inflection point index; j is the index; Q[j] represents the injection rate of the j+1th time period, and its value range is from 0 to i.

[0011] Optionally, in one embodiment of the present application, the grouting similarity index of the standard irrigation section is calculated based on the grouting parameters of the standard irrigation section, including: calculating the GIN value of the irrigation section of the standard irrigation section that meets preset conditions according to the injection volume duration data of the standard irrigation section, and obtaining the energy index of the standard irrigation section according to the GIN value of the irrigation section that meets the preset conditions; obtaining the irrigation sequence index of the standard irrigation section according to the irrigation sequence of the standard irrigation section; determining the average slurry density of the fluctuation section of the standard irrigation section according to the slurry density duration data of the standard irrigation section, and judging whether the average slurry density of the fluctuation section is in the slurry ratio interval; when the average slurry density of the fluctuation section is in the slurry ratio interval, obtaining the fracture index of the standard irrigation section; obtaining the grouting similarity index of the standard irrigation section according to the energy index of the standard irrigation section, the irrigation sequence index and the fracture index.

[0012] Optionally, in one embodiment of the present application, after determining whether the normal test result meets the preset qualification conditions, it also includes: if the normal test result does not meet the preset qualification conditions, converting the non-normal data in the parameter array into normal data so that the normal data conforms to the normal distribution.

[0013] Optionally, in one embodiment of the present application, the standard deviation of the unit average injection rate of the fluctuation section is:

[0014]

[0015] Where n is the number of irrigation segments in a certain similar irrigation segment set; i is a counter, representing the average injection rate per unit of the fluctuation segment of the i-th irrigation segment; is the average value;

[0016] The upper and lower limits of the confidence interval are:

[0017]

[0018] Wherein, n is the number of irrigation segments in the set of similar irrigation segments; s is the standard deviation; z is the z score of the standard normal distribution; is the average value;

[0019] The injection rate interval is:

[0020]

[0021] Among them, Q max is the upper limit of injection rate, Q minis the lower limit of injection rate.

[0022] The second aspect of the present application provides an injection rate interval analysis device for curtain grouting at the toe of a reservoir dam in a pumped storage project, comprising: an acquisition module for acquiring grouting data of a grouting site and adding the grouting data to a historical data set; a detection module for detecting the type of the injection rate time process line of each irrigation section based on the injection rate time data of each irrigation section in the historical data set, and screening out irrigation sections that meet preset normal conditions from the types of the injection rate time process lines, and obtaining a standard irrigation section based on the irrigation sections that meet the preset normal conditions; a determination module for determining the fluctuation range of the injection rate time process line of the standard irrigation section based on the injection rate inflection point time of the standard irrigation section, and calculating the fluctuation unit average injection rate of the standard irrigation section within the fluctuation range; a calculation module for calculating the grouting similarity index of the standard irrigation section based on the grouting parameters of the standard irrigation section; a clustering module for obtaining the standard irrigation section based on the grouting similarity index. The grouting similarity mode of the standard grouting section is determined, and the standard grouting sections that meet the preset same grouting similarity mode conditions are clustered to obtain at least one similar grouting section set; a removal module is used to transmit the unit average injection rate of the fluctuation section in the at least one similar grouting section set to the parameter array of the corresponding similar grouting section set, and remove the outliers in the parameter array to obtain the parameter array after removing the outliers; a judgment module is used to perform a normality test on the parameter array after removing the outliers to obtain a normal test result, and judge whether the normal test result meets the preset qualified conditions; an analysis module is used to calculate the confidence interval of the unit average injection rate of the fluctuation section of the similar grouting section set when the confidence level meets the preset percentage condition according to the normal distribution when the normal test result meets the preset qualified conditions, and determine the injection rate interval when the toe curtain grouting of the upper reservoir dam of the pumped storage project is in the grouting similarity mode according to the confidence interval.

[0023] Optionally, in one embodiment of the present application, the detection module includes: a first determination unit, used to determine the occurrence position of the injection rate inflection point of each irrigation section based on the continuous change law between the injection rate and time of each irrigation section; a judgment unit, used to judge whether the injection rate meets the preset specification conditions according to the occurrence position; and a classification unit, used to classify the line shape of the injection rate time process line when the injection rate meets the preset specification conditions, so as to obtain the irrigation section that meets the preset normal conditions.

[0024] Optionally, in one embodiment of the present application, the calculation formula for the unit average injection rate of the fluctuation section of the standard injection section is:

[0025]

[0026] Wherein, L is the length of the standard injection section; Q is the injection rate time series array; i is the injection rate inflection point index; j is the index; Q[j] represents the injection rate of the j+1th time period.

[0027] Optionally, in one embodiment of the present application, the calculation module includes: a calculation unit, used to calculate the GIN value of the standard irrigation section that meets the preset conditions according to the injection volume duration data of the standard irrigation section, and obtain the energy index of the standard irrigation section according to the GIN value of the irrigation section that meets the preset conditions; a first acquisition unit, used to obtain the irrigation sequence index of the standard irrigation section according to the irrigation sequence of the standard irrigation section; a second determination unit, used to determine the average slurry density of the fluctuation section of the standard irrigation section according to the slurry density duration data of the standard irrigation section, and determine whether the average slurry density of the fluctuation section is in the slurry ratio interval; a second acquisition unit, used to obtain the fracture index of the standard irrigation section when the average slurry density of the fluctuation section is in the slurry ratio interval; a third acquisition unit, used to obtain the grouting similarity index of the standard irrigation section according to the energy index of the standard irrigation section, the irrigation sequence index and the fracture index.

[0028] Optionally, in one embodiment of the present application, it also includes: a conversion module, which is used to convert the non-normal data in the parameter array into normal data after determining whether the normal test result meets the preset qualification conditions, if the normal test result does not meet the preset qualification conditions, so that the normal data conforms to the normal distribution.

[0029] Optionally, in one embodiment of the present application, the standard deviation of the unit average injection rate of the fluctuation section is:

[0030]

[0031] Where n is the number of irrigation segments in a certain similar irrigation segment set; i is a counter, representing the average injection rate per unit of the fluctuation segment of the i-th irrigation segment; is the average value;

[0032] The upper and lower limits of the confidence interval are:

[0033]

[0034] Wherein, n is the number of irrigation segments in the set of similar irrigation segments; s is the standard deviation; z is the z score of the standard normal distribution; is the average value;

[0035] The injection rate interval is:

[0036]

[0037] in, is the upper limit of injection rate, is the lower limit of injection rate.

[0038] The third aspect of the present application provides an electronic device, including: a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the program to implement the injection rate interval analysis method for curtain grouting at the toe of a reservoir dam in a pumped storage project as described in the above embodiment.

[0039] The fourth aspect of the present application provides a computer-readable storage medium, which stores a computer program. When the program is executed by a processor, it implements the injection rate interval analysis method for the curtain grouting at the toe of the reservoir dam in the pumped storage project as mentioned above.

[0040] The fifth aspect of the present application provides a computer program product, which stores a computer program that, when executed by a processor, implements the above-mentioned injection rate interval analysis method for curtain grouting at the toe of a reservoir dam in a pumped storage project.

[0041] The embodiment of the present application can use data mining technology, starting from historical grouting data, analyze the linear characteristics of the injection rate of each grouting section, use the linear classification algorithm to select the standard grouting section with normal linear shape, calculate the unit average injection rate of the fluctuation section, and then calculate the grouting similarity index of multiple categories to subdivide and cluster the standard grouting section to obtain a similar grouting section set, use a variety of mathematical statistical methods to remove outliers and normalize the unit average injection rate data of the fluctuation section in the set, and finally calculate the confidence interval of the unit average injection rate of the fluctuation section under each set when the confidence level is 95%, and obtain the injection rate interval of the curtain grouting at the toe of the reservoir dam in the pumped storage project under each grouting similarity mode, so as to realize the feedback of the curtain grouting construction based on the historical grouting data. Thus, the problem that the grouting construction party in the related technology can only confirm the reasonable injection rate interval based on the on-site water pressure test before the formal grouting and the personal experience of the technicians is solved, resulting in the pressure adjustment in the grouting process being relatively extensive, the grouting construction being empirical and fuzzy, and the grouting quality being difficult to guarantee.

[0042] Additional aspects and advantages of the present application will be given in part in the description below, and in part will become apparent from the description below, or will be learned through the practice of the present application. BRIEF DESCRIPTION OF THE DRAWINGS

[0043] The above and / or additional aspects and advantages of the present application will become apparent and easily understood from the following description of the embodiments in conjunction with the accompanying drawings, in which:

[0044] Figure 1 A flow chart of an injection rate interval analysis method for curtain grouting at the toe of a dam in a pumped storage project provided according to an embodiment of the present application;

[0045] Figure 2 A flowchart of calculating an injection rate trend array in a method for determining a line shape of an injection rate time course line according to an embodiment of the present application;

[0046] Figure 3 A flow chart for determining the linear shape in a method for determining the linear shape of an injection rate time course line according to an embodiment of the present application;

[0047] Figure 4 It is a normal type schematic diagram of the injection rate time course of the curtain grouting injection rate of the dam toe of a pumped storage project according to one embodiment of the present application;

[0048] Figure 5 A flowchart of clustering similar grouting sections based on grouting similarity indicators according to one embodiment of the present application;

[0049] Figure 6 A flow chart for calculating a confidence interval of a unit average injection rate of a fluctuation section of each set of similar injection sections according to an embodiment of the present application;

[0050] Figure 7 It is a structural schematic diagram of an injection rate interval analysis device for curtain grouting of the upper reservoir dam of a pumped storage project provided according to an embodiment of the present application;

[0051] Figure 8 It is a schematic diagram of the structure of an electronic device provided according to an embodiment of the present application. DETAILED DESCRIPTION

[0052] Embodiments of the present application are described in detail below, and examples of the embodiments are shown in the accompanying drawings, wherein the same or similar reference numerals throughout represent the same or similar elements or elements having the same or similar functions. The embodiments described below with reference to the accompanying drawings are exemplary and are intended to be used to explain the present application, and should not be construed as limiting the present application.

[0053] The following describes the injection rate interval analysis method and device for curtain grouting at the toe of a reservoir dam in a pumped storage project with reference to the accompanying drawings. In view of the problem that the grouting construction party in the related technology mentioned in the above background technology can only confirm the reasonable injection rate interval based on the on-site water pressure test before the formal grouting and rely on the personal experience of the technicians, resulting in rough pressure regulation during the grouting process, and the empiricism and fuzziness of the grouting construction, which makes it difficult to ensure the grouting quality, the present application provides a method for analyzing the injection rate interval of curtain grouting at the toe of a reservoir dam in a pumped storage project. In this method, data mining technology can be used to analyze the linear characteristics of the injection rate process line of each grouting section based on historical grouting data, and linear classification can be used. The algorithm selects the standard grouting section with normal linear shape, calculates the unit average injection rate of the fluctuation section, and then subdivides and clusters the standard grouting section to obtain a set of similar grouting sections by calculating multiple categories of grouting similarity indicators. A variety of mathematical statistical methods are used to remove outliers and normalize the unit average injection rate data of the fluctuation section in the set. Finally, by calculating the confidence interval of the unit average injection rate of the fluctuation section under each set with a confidence level of 95%, the injection rate interval of the curtain grouting at the toe of the reservoir dam in the pumped storage project under each grouting similarity mode is obtained, and feedback on the curtain grouting construction based on historical grouting data is realized. In this way, the grouting construction party in the relevant technology can only confirm the reasonable injection rate interval based on the on-site water pressure test before the formal grouting and the personal experience of the technicians, resulting in the pressure adjustment in the grouting process being relatively extensive, the grouting construction being empirical and fuzzy, and the grouting quality being difficult to guarantee.

[0054] Specifically, Figure 1 A flow chart of an injection rate interval analysis method for curtain grouting at the toe of a dam in a pumped storage project provided in an embodiment of the present application.

[0055] like Figure 1 As shown in FIG. 1 , the injection rate interval analysis method for the curtain grouting of the upper reservoir dam toe of the pumped storage project includes the following steps:

[0056] In step S101, the grouting data of the grouting site is obtained, and the grouting data is added to the historical data set.

[0057] It can be understood that the grouting site in the embodiment of the present application can be a grouting construction site.

[0058] During the actual execution process, the embodiment of the present application can request the grouting data of the fixed structure from the grouting construction site server through the interface, that is, when the relevant interface detects that there is data returned, the grouting data of the fixed structure can be collected from the grouting construction site at a regular interval according to the interface protocol, and the original grouting data structure can be organized into a processable grouting data structure using a processing algorithm, and then added to the historical data set, and the subsequent calculation process is triggered, thereby realizing real-time monitoring of the grouting data, and facilitating the use of historical data for feedback on the construction.

[0059] In step S102, based on the injection rate time data of each irrigation section in the historical data set, the type of the injection rate time process line of each irrigation section is detected, and the irrigation sections that meet the preset normal conditions are screened out from the types of the injection rate time process lines, and the standard irrigation sections are obtained based on the irrigation sections that meet the preset normal conditions.

[0060] It can be understood that, in the embodiment of the present application, the irrigation section that meets the preset normal condition may be an irrigation section whose injection rate time process line has a normal line shape.

[0061] As a possible implementation method, the embodiment of the present application can calculate the change trend of the injection rate of each irrigation section in each time period based on the injection rate history data of each irrigation section in the historical data set, represent each change trend with a number, and then put the trend numbers into an array in time sequence, which is called a trend array. Table 1 is a judgment basis table of injection rate change trend, where, as shown in Table 1:

[0062] Table 1

[0063]

[0064] like Figure 3 As shown, the embodiment of the present application can use the sliding window technology to perform linear judgment on the injection rate time process line of each irrigation section, and obtain the injection rate inflection point and the position where the injection rate enters the end segment node. First, initialize a window of size 2, and the window will contain two consecutive elements in the trend array at any time point. Then, starting from the first element of the trend array, slide the window across the trend array, and each time it moves, the window will move forward by a trend number. After each window slide, the embodiment of the present application can execute two monitoring goals:

[0065] Objective 1: Check whether both numbers in the window are -2. If yes, the embodiment of the present application can record the index of the first number in the window in the array, that is, the position of the injection rate inflection point, but the window continues to slide;

[0066] Objective 2: Check whether the number behind the window is 0. If the number behind the window is 0, the embodiment of the present application can record the position of this 0 and stop the sliding of the window.

[0067] If both results are returned when the window stops sliding, and the index value of the first target record is greater than 1, the line shape of the injection rate duration process line of the irrigation section is judged to be a normal line shape, and the irrigation section is judged to be a standard irrigation section. The line shape of the normal injection rate duration process line is as follows: Figure 4 shown.

[0068] It should be noted that the preset normal conditions can be set by those skilled in the art according to actual conditions and are not specifically limited here.

[0069] Optionally, in one embodiment of the present application, the type of the injection rate time process line of each irrigation section is detected, and the irrigation sections that meet the preset normal conditions are screened out from the types of the injection rate time process lines, including: determining the occurrence position of the injection rate inflection point of each irrigation section based on the continuous change law between the injection rate and time of each irrigation section; judging whether the injection rate meets the preset specification conditions according to the occurrence position; and when the injection rate meets the preset specification conditions, classifying the line shape of the injection rate time process line to obtain the irrigation sections that meet the preset normal conditions.

[0070] It can be understood that the determination of whether the injection rate meets the preset specification condition in the embodiment of the present application may be a condition for determining whether the injection rate drops below the specification value.

[0071] In the actual implementation process, the embodiment of the present application can use the continuous change law of the injection rate over time to determine the location of the injection rate inflection point, and judge whether the injection rate drops below the standard value by the location of the inflection point, so as to classify the line shape of the injection rate process line. Among them, the basis for judging the occurrence of the injection rate inflection point is: let a certain time node be t, the interval between each time node be T, and the injection rate be represented by Q, if The node t is considered to be the injection rate inflection point, and the corresponding time is the injection rate inflection point appearance time. If in the injection rate time course line of a certain irrigation section, there is a time node where the injection rate is less than 1L / min, and this node appears after the injection rate inflection point, then the injection rate time course line is judged to be normal, and the irrigation section is a standard irrigation section.

[0072] It should be noted that the preset standard conditions can be set by those skilled in the art according to actual conditions and are not specifically limited here.

[0073] In step S103, based on the injection rate inflection point time of the standard irrigation section, the fluctuation range of the injection rate duration process line of the standard irrigation section is determined, and the fluctuation unit average injection rate of the standard irrigation section is calculated within the fluctuation range.

[0074] It can be understood that the grouting injection rate represents the amount of slurry actually injected into the rock mass per minute. The normal curtain grouting injection rate process line has three stages: injection rate fluctuation section, injection rate sharp drop section and grouting end section. The injection rate in the fluctuation section changes less with time, and the injection rate in the sharp drop section will drop rapidly with time. The intersection of the injection rate fluctuation section and the injection rate sharp drop section is the injection rate inflection point, and the corresponding time is the injection rate inflection point occurrence time. Therefore, the key to determining the fluctuation section time interval is to determine the injection rate inflection point time. The average injection rate of the fluctuation section can be obtained by averaging the injection rate in the fluctuation section.

[0075] In the actual implementation process, the embodiment of the present application has obtained the injection rate inflection point position of each standard irrigation section in the above step S102, and the part of the injection rate process line before the injection rate inflection point is a fluctuation section (including the injection rate inflection point itself), and the irrigation section is determined, that is, the length of the irrigation section is L.

[0076] It should be noted that the embodiment of the present application proposes the concept of a unit average injection rate in a fluctuation section. In the injection rate time course line, the portion before the injection rate inflection point is the fluctuation section. The injection rate of this portion does not change much over time. The average injection rate of this portion can be obtained by averaging the injection rate of this portion, and the unit is L / min. Because the length of a single grouting section is determined according to the total depth of the entire grouting hole, there are slight differences in the length of the grouting sections between different grouting holes. The present application assumes that the internal cracks of a single grouting section are evenly distributed along the depth direction, and the average injection rate of the fluctuation section per unit length is taken as the unit average injection rate of the fluctuation section, and the unit is L / (min*times*m).

[0077] In one embodiment of the present application, the calculation formula for the unit average injection rate of the fluctuation section of the standard injection section is:

[0078]

[0079] Where L is the length of the standard injection section; Q is the injection rate time series array; i is the injection rate inflection point index; j is the index; Q[j] represents the injection rate of the j+1th time period, and its value range is from 0 to i.

[0080] The embodiment of the present application can determine the fluctuation range of the injection rate duration process line of the standard irrigation section based on the injection rate inflection point time of the standard irrigation section, and use this formula for each standard irrigation section to obtain the unit average injection rate of the fluctuation section of each standard irrigation section.

[0081] In step S104, based on the grouting parameters of the standard grouting section, the grouting similarity index of the standard grouting section is calculated.

[0082] In the actual implementation process, the embodiment of the present application can calculate or select the grouting similarity index of the standard grouting section based on the grouting parameters of the standard grouting section, thereby providing support for subsequently obtaining the grouting similarity mode of the standard grouting section.

[0083] Optionally, in one embodiment of the present application, based on the grouting parameters of the standard grouting section, the grouting similarity index of the standard grouting section is calculated, including: calculating the GIN value of the grouting section that meets the preset conditions according to the injection volume duration data of the standard grouting section, and obtaining the energy index of the standard grouting section according to the GIN value of the grouting section that meets the preset conditions; obtaining the grouting sequence index of the standard grouting section according to the grouting sequence of the standard grouting section; determining the average grouting density of the fluctuation section of the standard grouting section according to the grouting density duration data of the standard grouting section, and judging whether the average grouting density of the fluctuation section is in the grouting ratio interval; when the average grouting density of the fluctuation section is in the grouting ratio interval, obtaining the fracture index of the standard grouting section; obtaining the grouting similarity index of the standard grouting section according to the energy index, grouting sequence index and fracture index of the standard grouting section.

[0084] It can be understood that the GIN value of the irrigation section that meets the preset conditions in the embodiment of the present application can be the maximum GIN value of each standard irrigation section.

[0085] In the actual implementation process, the embodiment of the present application can determine the maximum injection volume and the corresponding gauge pressure value of each standard irrigation section according to the injection volume duration data of each standard irrigation section, and then calculate the maximum GIN value of each standard irrigation section, and determine the interval in which it is located, and the determination result is used as the energy index. Specifically, the maximum injection volume V in the injection volume duration array of each standard irrigation section is calculated. max = max(V), record the index of the maximum injection volume in the injection volume duration array, and then take out the table pressure value P of the corresponding index in the table pressure duration array max , and then calculate the maximum GIN value of each standard irrigation section: GIN = P max ×V max , where V max is the maximum injection volume per irrigation section, in L / m, P max It is the surface pressure when the grouting section reaches the maximum injection volume, in MPa. This application sets 6 energy indicators according to the recommended GIN grouting limit envelope, which are very high, high, medium, low, very low and almost no energy input. Table 2 is the energy indicator judgment basis table, as shown in Table 2:

[0086] Table 2

[0087]

[0088] The maximum GIN value, maximum injection volume, and maximum injection volume correspond to the table pressure of each standard irrigation section, and the interval is judged according to the energy index judgment table, and finally the energy index of each standard irrigation section is obtained.

[0089] Furthermore, the embodiment of the present application can take the grouting sequence of each standard grouting section as the grouting sequence indicator. Specifically, the grouting sequence is a concept related to the grouting holes in the curtain grouting, which has been determined at the beginning of the design, and each grouting section belongs to a certain grouting hole, so the grouting sequence is included in the returned data structure. The embodiment of the present application only needs to select it from the data set when needed.

[0090] Furthermore, the embodiment of the present application can calculate the average density of the grouting liquid between the start time of grouting and the time of the inflection point of the injection rate, that is, determine the average grouting density of the fluctuation section of the standard grouting section based on the grouting density duration data of the standard grouting section, and judge whether the average grouting density of the fluctuation section is in the grouting ratio interval, obtain the judgment result, and use the judgment result as the crack index of the standard grouting section, and the crack index is determined by the average grouting density of the fluctuation section of the grouting section. Specifically, the embodiment of the present application has obtained the injection rate inflection point position of each standard grouting section in the above step S102, and the injection rate duration process line before the injection rate inflection point is the fluctuation section, so the average grouting density of the fluctuation section is the average density of the grout injected during the construction duration of the fluctuation section, which can be calculated based on the grouting density duration data of each standard grouting section and the injection rate inflection point position:

[0091]

[0092] Among them, ρ is the time series array of slurry density; i is the injection rate inflection point index; j is the index; ρ[j] represents the average slurry density in the j+1th time period, and its value range is from 0 to i.

[0093] The crack index in the embodiment of the present application is defined by the average density of the slurry. According to the actual recommended construction slurry replacement mix ratio, 5 crack indexes are set, namely very fine, fine, medium, coarse, and very coarse. The grouting section is further subdivided according to the grouting similarity mode, and a more targeted recommended injection rate interval can be obtained. Table 3 is a judgment basis table for crack indexes, wherein, as shown in Table 3:

[0094] Table 3

[0095]

[0096] It should be noted that the preset conditions can be set by those skilled in the art according to actual conditions and are not specifically limited here.

[0097] In step S105, the grouting similarity pattern of the standard grouting section is obtained according to the grouting similarity index, and the standard grouting sections that meet the preset same grouting similarity pattern conditions are clustered to obtain at least one similar grouting section set.

[0098] As a possible implementation method, the embodiment of the present application can combine the three similarity indicators obtained in each grouting section according to the above step S104, and the combination is called a grouting similarity pattern. The standard grouting sections with the same grouting similarity pattern are clustered to obtain multiple similar grouting section sets;

[0099] Specifically, the specific indicators obtained according to the above step S104 in the embodiment of the present application include:

[0100] (1) Energy index: very high, high, medium, low, very low, almost no energy input;

[0101] (2) Crack index: very fine, fine, medium, coarse, and very coarse;

[0102] (3) Grouting sequence index (determined according to the construction process): 0 sequence, 1 sequence, 2 sequence, 3 sequence, 4 sequence, 5 sequence.

[0103] In the embodiment of the present application, each grouting section can be combined according to specific indicators in the three types of similar indicators. Each combination is called a grouting similarity mode: grouting similarity mode = {energy indicator, crack indicator, grouting sequence indicator}, and the standard grouting sections with the same grouting similarity mode are clustered to obtain multiple similar grouting section sets: similar grouting section set = {X|X is a grouting section and X has a grouting similarity mode {mode 1, mode 2, mode 3...}}.

[0104] It should be noted that the preset same grouting similar mode conditions can be set by those skilled in the art according to actual conditions and are not specifically limited here.

[0105] In step S106, the unit average injection rate of the fluctuation segment in at least one similar irrigation segment set is transmitted to the parameter array of the corresponding similar irrigation segment set, and the outliers in the parameter array are removed to obtain the parameter array after the outliers are removed.

[0106] In the actual implementation process, the embodiment of the present application can put the average injection rate of the fluctuation section unit of all the irrigation sections in each similar irrigation section set into the parameter array of the corresponding similar irrigation section set, use the DBSCAN algorithm to identify the outliers in the array of the average injection rate of the fluctuation section unit of each similar irrigation section set, and remove the outliers in the parameter array to obtain the parameter array after removing the outliers. Specifically, Figure 5 As shown, in the embodiment of the present application, the average injection rate per unit of the fluctuation section of all irrigation sections in each similar irrigation section set can be placed in the parameter array of the corresponding similar irrigation section set: for the similar irrigation section set Y = {set Y1, set Y2, set Y3 ...}, the parameter array Then, for each set Y, the outliers in the parameter array Z are removed using the DBSCAN algorithm:

[0107] New parameter array In practical applications, the distance threshold eps=0.8 and the minimum number of samples min_samples=5 are taken.

[0108] In step S107, a normality test is performed on the parameter array after removing outliers to obtain a normal test result, and it is determined whether the normal test result meets the preset qualification conditions.

[0109] It can be understood that the preset qualified condition in the embodiment of the present application can be a condition for passing the normality test.

[0110] In the actual implementation process, the embodiment of the present application can use the Shapiro-Wilk method to perform a normality test on the parameter array of each similar segment set after removing outliers, obtain a normal test result, and determine whether the normal test result meets certain qualification conditions. Specifically, if Figure 5 As shown, the Shapiro-Wilk method is used to perform a normality test on the parameter set in the parameter array of each similar injection segment set. The algorithm steps include: (1) data sorting (2) calculating the order statistic (3) constructing the test statistic W (4) calculating the p-value. The value of the test statistic W is between 0 and 1. When W is close to 1, it means that the unit fluctuation segment injection rate data in the parameter set is very close to the normal distribution; when W is much less than 1, it means that the unit fluctuation segment injection rate data in the parameter set is significantly different from the normal distribution. If the p-value is less than the significance level of 0.05, the null hypothesis is rejected, and it is considered that the unit fluctuation segment injection rate data in the parameter set does not conform to the normal distribution.

[0111] It should be noted that the preset qualification conditions can be set by those skilled in the art according to actual conditions and are not specifically limited here.

[0112] Optionally, in one embodiment of the present application, after determining whether the normal test result meets the preset qualification conditions, it also includes: if the normal test result does not meet the preset qualification conditions, converting the non-normal data in the parameter array into normal data so that the normal data conforms to the normal distribution.

[0113] In the actual implementation process, the embodiment of the present application can use Box-Cox transformation to convert the non-normal data in the parameter array into normal data when the normal test result does not meet certain qualification conditions, that is, when the normality test believes that the unit fluctuation segment injection rate data in the parameter set does not conform to the normal distribution, so that the normal data conforms to the normal distribution. The parameter λ is estimated by the moment estimation method: λ = (second-order central moment - first-order central distance^2) / first-order central distance, where the second-order central moment is the variance of the parameter set and the first-order central distance is the mean of the parameter set.

[0114] It should be noted that the preset qualification conditions can be set by those skilled in the art according to actual conditions and are not specifically limited here.

[0115] In step S108, when the normal test results meet the preset qualification conditions, the confidence interval of the unit average injection rate of the fluctuation section of the similar irrigation section set is calculated according to the normal distribution when the confidence level meets the preset percentage conditions, and the injection rate interval when the curtain grouting at the toe of the upper reservoir dam of the pumped storage project is in the grouting similarity mode is determined according to the confidence interval.

[0116] It is understandable that during the curtain grouting construction process, insufficient injection rate may lead to under-grouting, and excessive injection rate may cause splitting, both of which will affect the quality of the curtain. Only by maintaining a suitable injection rate in the fluctuation section can the quality of the final curtain be up to standard. Therefore, the interval of the injection rate in the fluctuation section is calculated, and the injection rate is kept within the interval as much as possible, which is conducive to the final formation of a good quality waterproof curtain. Taking into account the coupling relationship between pressure and injection rate, the interval of the injection rate in the fluctuation section can also provide conditional constraints for pressure control. At the same time, based on the interval of the injection rate in the fluctuation section, the construction quality of the historical grouting section can also be reasonably evaluated. The preset percentage condition in the embodiment of the present application can be a 95% condition.

[0117] In the actual implementation process, the embodiment of the present application can use the standard normal distribution to calculate the confidence interval of the unit average injection rate of the fluctuation section for each similar irrigation section set when the confidence level is 95% when the normal test results meet certain qualification conditions, that is, when the parameter array of each similar irrigation section set conforms to the normal distribution. These confidence intervals are the injection rate intervals of the curtain grouting at the toe of the reservoir dam in the pumped storage project under each grouting similarity mode.

[0118] It should be noted that the preset qualification condition and the preset percentage condition can be set by those skilled in the art according to actual conditions and are not specifically limited here.

[0119] Among them, in one embodiment of the present application, Figure 6 As shown, for each set of similar injection sections, the standard deviation of the unit average injection rate of the fluctuation section is calculated when the confidence level is 95%:

[0120]

[0121] Where n is the number of irrigation segments in a certain similar irrigation segment set; i is a counter, representing the average injection rate of the fluctuation segment unit of the i-th irrigation segment, and its value range is 1 to n; is the average value.

[0122] The upper and lower limits of the confidence interval are:

[0123]

[0124] Where n is the number of irrigation sections in a certain set of similar irrigation sections; s is the standard deviation; z is the z score of the standard normal distribution, and z is 1.96; is the average value.

[0125] When the confidence intervals of all similar grouting sections are calculated, the injection rate intervals (grouting similarity modes) of the upper reservoir dam toe grouting in each grouting similarity mode are obtained.

[0126]

[0127] Among them, Q max is the upper limit of injection rate, Q min is the lower limit of injection rate.

[0128] The embodiment of the present application can calculate the injection rate range when curtain grouting is carried out at the toe of the reservoir dam in the pumped storage project. If the injection rate can be kept within this range, it will help increase the probability of finally forming a good quality waterproof curtain. This range also provides constraints on the pressure control during the grouting process. The timing of pressurization or depressurization and the size of the adjustment range will no longer rely solely on the experience of the technicians. This lays a strong technical foundation for truly realizing the automated construction of curtain grouting at the toe of the reservoir dam in the pumped storage project.

[0129] Specifically, it can be combined Figure 2 As shown, the working principle of the injection rate interval analysis method for curtain grouting at the toe of the upper reservoir dam of the pumped storage project in the embodiment of the present application is described in detail with a specific embodiment.

[0130] like Figure 2 As shown, the embodiment of the present application may include the following steps:

[0131] Step S201: Enter the injection rate trend calculation module.

[0132] Step S202: Read the injection rate array (1 to N) to define the state array, and define i=1.

[0133] Step S203: Determine whether array (i)>1 holds, if so, execute step S204; if not, execute step S205.

[0134] Step S204: Array (i+1) = Array (i).

[0135] Step S205: state array (i)=0, i=i+1.

[0136] Step S206: (array (i+1)-array (i)) / array (i)<-0.15 is established. If so, execute step S208; if not, execute step S209.

[0137] Step S207: state array (i)=1, i=i+1.

[0138] Step S208: state array (i) = -2, i = i + 1.

[0139] Step S209: state array (i) = -1, i = i+1.

[0140] Step S210: i=N.

[0141] Step S211: exit the injection rate trend calculation module.

[0142] Furthermore, if Figure 3 As shown, the embodiment of the present application may include the following steps:

[0143] Step S301: Enter the line shape determination module.

[0144] Step S302: Read the status array to define the window size as 2.

[0145] Step S303: Window sliding.

[0146] Step S304: Determine whether the number of windows (2) = 0 is true, if so, execute step S305; if not, execute step S306.

[0147] Step S305: Record the position B of the window number (2) in the state array.

[0148] Step S306: Determine whether the number of windows (1) = -2 and the number of windows (2) = -2 are true. If so, execute step S307; if not, execute step S308.

[0149] Step S307: Record the position A of the window number (1) in the state array.

[0150] Step S308: Slide to the end.

[0151] Step S309: Determine whether both A and B exist. If so, execute step S310; if not, execute step S311.

[0152] Step S310: The line shape is judged to be normal.

[0153] Step S311: exit the flag configuration module.

[0154] In the embodiment of the present application, the normal injection rate duration process line has a line shape as follows: Figure 4 shown.

[0155] Furthermore, if Figure 5 As shown, the embodiment of the present application may include the following steps:

[0156] Step S501: Enter the similar irrigation segment clustering module.

[0157] Step S502: Circularly read the standard filling sections.

[0158] Step S503: read three similarity indicators of the current standard irrigation section.

[0159] Step S504: grouting similarity mode J=energy index X+grouting sequence index Y+crack index Z.

[0160] Step S505: Determine whether the grouting similarity pattern J already exists. If so, execute step S506; if not, execute step S507.

[0161] Step S506: Add grouting similarity mode J to the grouting section.

[0162] Step S507: create a new grouting similarity mode J, create a new corresponding grouting section set, and create a new corresponding parameter array.

[0163] Step S508: The unit average injection rate of the fluctuation section of the grouting section is added to the grouting similarity model J.

[0164] Step S509: Determine whether there are still standard filling sections. If so, execute step S502; if not, execute step S510.

[0165] Step S510: exit the similar segment clustering module.

[0166] Furthermore, if Figure 6 As shown, the embodiment of the present application may include the following steps:

[0167] Step S601: Enter the confidence interval calculation module.

[0168] Step S602: Read a set of similar segments.

[0169] Step S603: Read the average injection rate array of the fluctuation segment unit of the current similar injection segment set.

[0170] Step S604: Use DBSCAN algorithm to remove outliers.

[0171] Step S605: Determine whether the SW method test array p>0.05 is established. If so, execute step S606; if not, execute step S607.

[0172] Step S606: 95% confidence interval calculation.

[0173] Step S607: Box-Cox transformation.

[0174] Step S608: Determine whether there is still a similar segment set, if so, execute step S602; if not, execute step S609.

[0175] Step S609: exit the confidence interval calculation module and obtain the recommended injection rate interval.

[0176] According to the injection rate interval analysis method of the curtain grouting at the toe of the reservoir dam in the pumped storage project proposed in the embodiment of the present application, data mining technology can be used to analyze the linear characteristics of the injection rate of each grouting section from the historical grouting data, and the standard grouting section with normal linear shape can be selected by using the linear classification algorithm, and the unit average injection rate of the fluctuation section is calculated. Then, by calculating the grouting similarity index of multiple categories, the standard grouting section is subdivided and clustered to obtain a similar grouting section set, and a variety of mathematical statistical methods are used to remove outliers and normalize the unit average injection rate data of the fluctuation section in the set. Finally, by calculating the confidence interval of the unit average injection rate of the fluctuation section under each set when the confidence level is 95%, the injection rate interval of the curtain grouting at the toe of the reservoir dam in the pumped storage project under each grouting similarity mode is obtained, and feedback on the curtain grouting construction based on the historical grouting data is realized. Thus, the problem that the related technology has realized real-time monitoring of grouting data, but still cannot use historical data to provide feedback on the construction, resulting in the excessive subjectivity and experience of the curtain grouting construction operation, thereby affecting the quality of the final curtain.

[0177] Next, the injection rate interval analysis device for curtain grouting at the toe of the upper reservoir dam of a pumped storage project proposed in accordance with an embodiment of the present application will be described with reference to the accompanying drawings.

[0178] Figure 7 It is a structural schematic diagram of an injection rate interval analysis device for curtain grouting at the toe of a dam in a pumped storage project according to an embodiment of the present application.

[0179] like Figure 7 As shown, the injection rate interval analysis device 10 for the curtain grouting at the toe of the upper reservoir dam of the pumped storage project includes: an acquisition module 100, a detection module 200, a determination module 300, a calculation module 400, a clustering module 500, a removal module 600, a judgment module 700 and an analysis module 800.

[0180] Specifically, the acquisition module 100 is used to acquire the grouting data of the grouting site and add the grouting data to the historical data set.

[0181] The detection module 200 is used to detect the type of the injection rate time process line of each irrigation section based on the injection rate time data of each irrigation section in the historical data set, and to filter out the irrigation sections that meet the preset normal conditions from the types of the injection rate time process lines, and to obtain the standard irrigation sections based on the irrigation sections that meet the preset normal conditions.

[0182] The determination module 300 is used to determine the fluctuation range of the injection rate duration process line of the standard irrigation section based on the injection rate inflection point time of the standard irrigation section, and calculate the fluctuation unit average injection rate of the standard irrigation section within the fluctuation range.

[0183] The calculation module 400 is used to calculate the grouting similarity index of the standard grouting section based on the grouting parameters of the standard grouting section.

[0184] The clustering module 500 is used to obtain the grouting similarity pattern of the standard grouting section according to the grouting similarity index, and cluster the standard grouting sections that meet the preset same grouting similarity pattern conditions to obtain at least one similar grouting section set.

[0185] The elimination module 600 is used to transmit the unit average injection rate of the fluctuation segment in at least one similar irrigation segment set to the parameter array corresponding to the similar irrigation segment set, and eliminate the outliers in the parameter array to obtain the parameter array after eliminating the outliers.

[0186] The judgment module 700 is used to perform a normality test on the parameter array after removing outliers, obtain a normal test result, and judge whether the normal test result meets a preset qualification condition.

[0187] Analysis module 800 is used to calculate the confidence interval of the unit average injection rate of the fluctuation section of the similar irrigation section set according to the normal distribution when the confidence level meets the preset percentage condition when the normal test results meet the preset qualification conditions, and determine the injection rate interval when the curtain grouting at the toe of the upper reservoir dam of the pumped storage project is in the grouting similarity mode according to the confidence interval.

[0188] Optionally, in one embodiment of the present application, the detection module 200 includes: a first determination unit, a judgment unit and a classification unit.

[0189] The first determination unit is used to determine the position of the injection rate inflection point of each irrigation section based on the continuous variation rule between the injection rate of each irrigation section and time.

[0190] The judging unit is used to judge whether the injection rate meets the preset specification conditions according to the occurrence position.

[0191] The classification unit is used to classify the line shape of the injection rate process line when the injection rate meets the preset standard conditions, so as to obtain the irrigation section that meets the preset normal conditions.

[0192] Optionally, in one embodiment of the present application, the calculation formula for the unit average injection rate of the fluctuation section of the standard injection section is:

[0193]

[0194] Wherein, L is the length of the standard injection section; Q is the injection rate time series array; i is the injection rate inflection point index; j is the index; Q[j] represents the injection rate of the j+1th time period.

[0195] Optionally, in one embodiment of the present application, the calculation module 400 includes: a calculation unit, a first acquisition unit, a second determination unit, a second acquisition unit and a third acquisition unit.

[0196] The calculation unit is used to calculate the GIN value of the standard irrigation section that meets the preset conditions according to the injection volume duration data of the standard irrigation section, and obtain the energy index of the standard irrigation section according to the GIN value of the irrigation section that meets the preset conditions.

[0197] The first acquisition unit is used to obtain the filling sequence index of the standard filling section according to the filling sequence of the standard filling section.

[0198] The second determination unit is used to determine the average slurry density of the fluctuation section of the standard filling section according to the slurry density history data of the standard filling section, and judge whether the average slurry density of the fluctuation section is in the slurry ratio range.

[0199] The second acquisition unit is used to acquire the crack index of the standard irrigation section when the average slurry density of the fluctuation section is within the slurry ratio range.

[0200] The third acquisition unit is used to obtain the grouting similarity index of the standard grouting section according to the energy index, grouting sequence index and crack index of the standard grouting section.

[0201] Optionally, in one embodiment of the present application, the injection rate interval analysis device 10 for curtain grouting of the reservoir dam toe in a pumped storage project further includes: a conversion module.

[0202] Among them, the conversion module is used to convert the non-normal data in the parameter array into normal data after judging whether the normal test result meets the preset qualification conditions, if the normal test result does not meet the preset qualification conditions, so that the normal data conforms to the normal distribution.

[0203] Optionally, in one embodiment of the present application, the standard deviation of the unit average injection rate of the fluctuation section is:

[0204]

[0205] Where n is the number of irrigation segments in a certain similar irrigation segment set; i is a counter, representing the average injection rate per unit of the fluctuation segment of the i-th irrigation segment; is the average value.

[0206] The upper and lower limits of the confidence interval are:

[0207]

[0208] Where n is the number of irrigation segments in a set of similar irrigation segments; s is the standard deviation; z is the z score of the standard normal distribution; is the average value.

[0209] The injection rate range is:

[0210]

[0211] in, is the upper limit of injection rate, is the lower limit of injection rate.

[0212] It should be noted that the aforementioned explanation of the embodiment of the injection rate interval analysis method for curtain grouting at the toe of a reservoir dam in a pumped storage project is also applicable to the injection rate interval analysis device for curtain grouting at the toe of a reservoir dam in a pumped storage project in this embodiment, and will not be repeated here.

[0213] According to the injection rate interval analysis device for curtain grouting at the toe of the reservoir dam in the pumped storage project proposed in the embodiment of the present application, data mining technology can be used to analyze the linear characteristics of the injection rate of each grouting section from the historical grouting data, and the standard grouting section with normal linear shape can be selected by using the linear classification algorithm, and the unit average injection rate of the fluctuation section is calculated. Then, by calculating the grouting similarity index of multiple categories, the standard grouting section is subdivided and clustered to obtain a similar grouting section set, and a variety of mathematical statistical methods are used to remove outliers and normalize the unit average injection rate data of the fluctuation section in the set. Finally, by calculating the confidence interval of the unit average injection rate of the fluctuation section under each set when the confidence level is 95%, the injection rate interval of the curtain grouting at the toe of the reservoir dam in the pumped storage project under each grouting similarity mode is obtained, and feedback on the curtain grouting construction based on the historical grouting data is realized. Thus, the problem that the related technology has realized real-time monitoring of grouting data, but still cannot use historical data to provide feedback on the construction, resulting in the excessive subjectivity and experience of the curtain grouting construction operation, thereby affecting the quality of the final curtain.

[0214] Figure 8 A schematic diagram of the structure of an electronic device provided in an embodiment of the present application. The electronic device may include:

[0215] A memory 801 , a processor 802 , and a computer program stored in the memory 801 and executable on the processor 802 .

[0216] When the processor 802 executes the program, the injection rate interval analysis method for the curtain grouting of the upper reservoir dam toe of the pumped storage project provided in the above embodiment is implemented.

[0217] Furthermore, the electronic device further comprises:

[0218] The communication interface 803 is used for communication between the memory 801 and the processor 802 .

[0219] The memory 801 is used to store computer programs that can be executed on the processor 802 .

[0220] The memory 801 may include a high-speed RAM memory, and may also include a non-volatile memory (non-volatile memory), such as at least one disk memory.

[0221] If the memory 801, the processor 802 and the communication interface 803 are implemented independently, the communication interface 803, the memory 801 and the processor 802 can be connected to each other through a bus and communicate with each other. The bus can be an Industry Standard Architecture (ISA) bus, a Peripheral Component (PCI) bus or an Extended Industry Standard Architecture (EISA) bus. The bus can be divided into an address bus, a data bus, a control bus, etc. For ease of representation, Figure 8 Only one thick line is used in the diagram, but this does not mean that there is only one bus or only one type of bus.

[0222] Optionally, in a specific implementation, if the memory 801, the processor 802 and the communication interface 803 are integrated on a chip, the memory 801, the processor 802 and the communication interface 803 can communicate with each other through an internal interface.

[0223] The processor 802 may be a central processing unit (CPU), or an application specific integrated circuit (ASIC), or one or more integrated circuits configured to implement the embodiments of the present application.

[0224] This embodiment also provides a computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the above injection rate interval analysis method for curtain grouting at the toe of a reservoir dam in a pumped storage project.

[0225] An embodiment of the present application also provides a computer program product on which a computer program is stored, and when the program is executed by a processor, the injection rate interval analysis method for curtain grouting at the toe of a reservoir dam in a pumped storage project as described above is implemented.

[0226] In the description of this specification, the description with reference to the terms "one embodiment", "some embodiments", "example", "specific example", or "some examples" etc. means that the specific features, structures, materials or characteristics described in conjunction with the embodiment or example are included in at least one embodiment or example of the present application. In this specification, the schematic representations of the above terms do not necessarily refer to the same embodiment or example. Moreover, the specific features, structures, materials or characteristics described may be combined in any one or N embodiments or examples in a suitable manner. In addition, those skilled in the art may combine and combine the different embodiments or examples described in this specification and the features of the different embodiments or examples, without contradiction.

[0227] In addition, the terms "first" and "second" are used for descriptive purposes only and should not be understood as indicating or implying relative importance or implicitly indicating the number of technical features indicated. Therefore, a feature defined as "first" or "second" may explicitly or implicitly include at least one of the features. In the description of this application, "N" means at least two, such as two, three, etc., unless otherwise clearly and specifically defined.

[0228] Any process or method description in a flowchart or otherwise described herein may be understood to represent a module, fragment or portion of code comprising one or N executable instructions for implementing the steps of a custom logical function or process, and the scope of the preferred embodiments of the present application includes alternative implementations in which functions may not be performed in the order shown or discussed, including performing functions in a substantially simultaneous manner or in reverse order depending on the functions involved, which should be understood by technicians in the technical field to which the embodiments of the present application belong.

[0229] The logic and / or steps represented in the flowchart or otherwise described herein, for example, can be considered as an ordered list of executable instructions for implementing logical functions, and can be embodied in any computer-readable medium for use by an instruction execution system, device or apparatus (such as a computer-based system, a system including a processor, or other system that can fetch instructions from an instruction execution system, device or apparatus and execute instructions), or in combination with these instruction execution systems, devices or apparatuses. For the purpose of this specification, "computer-readable medium" can be any device that can contain, store, communicate, propagate or transmit a program for use by an instruction execution system, device or apparatus, or in combination with these instruction execution systems, devices or apparatuses. More specific examples of computer-readable media (a non-exhaustive list) include the following: an electrical connection with one or N wirings (electronic devices), a portable computer disk box (magnetic device), a random access memory (RAM), a read-only memory (ROM), an erasable and programmable read-only memory (EPROM or flash memory), a fiber optic device, and a portable compact disk read-only memory (CDROM). In addition, the computer-readable medium may even be paper or other suitable medium on which the program is printed, since the program may be obtained electronically by optically scanning the paper or other medium and then editing, interpreting or processing in other suitable ways as necessary and then storing it in a computer memory.

[0230] It should be understood that the various parts of the present application can be implemented by hardware, software, firmware or a combination thereof. In the above-mentioned embodiment, the N steps or methods can be implemented by software or firmware stored in a memory and executed by a suitable instruction execution system. For example, if implemented by hardware, as in another embodiment, it can be implemented by any one of the following technologies known in the art or their combination: a discrete logic circuit having a logic gate circuit for implementing a logic function for a data signal, a dedicated integrated circuit having a suitable combination of logic gate circuits, a programmable gate array (PGA), a field programmable gate array (FPGA), etc.

[0231] A person skilled in the art may understand that all or part of the steps in the method for implementing the above-mentioned embodiment may be completed by instructing related hardware through a program, and the program may be stored in a computer-readable storage medium, which, when executed, includes one or a combination of the steps of the method embodiment.

[0232] In addition, each functional unit in each embodiment of the present application may be integrated into a processing module, or each unit may exist physically separately, or two or more units may be integrated into one module. The above-mentioned integrated module may be implemented in the form of hardware or in the form of a software functional module. If the integrated module is implemented in the form of a software functional module and sold or used as an independent product, it may also be stored in a computer-readable storage medium.

[0233] The storage medium mentioned above may be a read-only memory, a magnetic disk or an optical disk, etc. Although the embodiments of the present application have been shown and described above, it can be understood that the above embodiments are exemplary and cannot be understood as limiting the present application. A person of ordinary skill in the art may change, modify, replace and modify the above embodiments within the scope of the present application.

Claims

1. A method for analyzing the injection rate interval of curtain grouting at the toe of a pumped storage project, characterized in that: The following steps are involved: Acquire grouting data of the grouting site, and add the grouting data to the historical data set; Based on the injection rate duration data of each irrigation section in the historical data set, the type of the injection rate duration process line of each irrigation section is detected, and the irrigation sections that meet the preset normal conditions are screened out from the types of the injection rate duration process lines, and the standard irrigation sections are obtained according to the irrigation sections that meet the preset normal conditions; Based on the injection rate inflection point time of the standard irrigation section, determining the fluctuation section range of the injection rate duration process line of the standard irrigation section, and calculating the fluctuation section unit average injection rate of the standard irrigation section within the fluctuation section range; Calculating a grouting similarity index of the standard grouting section based on the grouting parameters of the standard grouting section; Obtaining a grouting similarity pattern of the standard grouting section according to the grouting similarity index, and clustering the standard grouting sections that meet the preset same grouting similarity pattern conditions to obtain at least one similar grouting section set; The unit average injection rate of the fluctuation section in the at least one similar irrigation section set is transmitted to the parameter array of the corresponding similar irrigation section set, and outliers in the parameter array are removed to obtain the parameter array after the outliers are removed; Performing a normality test on the parameter array after removing outliers to obtain a normal test result, and determining whether the normal test result meets a preset qualification condition; When the normal test result satisfies the preset qualification condition, the confidence interval of the unit average injection rate of the fluctuation section of the similar irrigation section set is calculated according to the normal distribution when the confidence level meets the preset percentage condition, and the injection rate interval of the curtain grouting at the toe of the reservoir dam of the pumped storage project when it is in the grouting similarity mode is determined according to the confidence interval.

2. The method according to claim 1, characterized in that The determining of the type of the injection rate time process line of each irrigation section, and selecting the irrigation section that meets the preset normal conditions from the types of the injection rate time process lines, includes: Based on the continuous variation law between the injection rate and time of each irrigation section, determining the occurrence position of the injection rate inflection point of each irrigation section; Determining whether the injection rate meets a preset specification condition according to the occurrence position; When the injection rate meets the preset standard conditions, the line shape of the injection rate time process line is classified to obtain the irrigation section that meets the preset normal conditions.

3. The method according to claim 1, characterized in that The calculation formula for the unit average injection rate of the fluctuation section of the standard irrigation section is: Wherein, L is the length of the standard injection section; Q is the injection rate time series array; i is the injection rate inflection point index; j is the index; Q[j] represents the injection rate of the j+1th time period.

4. The method according to claim 1, characterized in that: The calculating of the grouting similarity index of the standard grouting section based on the grouting parameters of the standard grouting section includes: Calculating the GIN value of the standard irrigation section that meets the preset conditions according to the injection volume duration data of the standard irrigation section, and obtaining the energy index of the standard irrigation section according to the GIN value of the irrigation section that meets the preset conditions; Obtaining an irrigation sequence index of the standard irrigation section according to the irrigation sequence of the standard irrigation section; Determine the average slurry density of the fluctuation section of the standard filling section according to the slurry density history data of the standard filling section, and judge whether the average slurry density of the fluctuation section is within the slurry ratio interval; When the average slurry density of the fluctuation section is within the slurry ratio range, obtaining the crack index of the standard irrigation section; The grouting similarity index of the standard grouting section is obtained according to the energy index of the standard grouting section, the grouting sequence index and the crack index.

5. The method according to claim 1, characterized in that After determining whether the normality test result meets the preset qualification condition, the method further includes: If the normal test result does not meet the preset qualification condition, the non-normal data in the parameter array is converted into normal data so that the normal data conforms to the normal distribution.

6. The method according to claim 1, characterized in that The standard deviation of the unit average injection rate of the fluctuation section is: Where n is the number of irrigation segments in a certain similar irrigation segment set; i is a counter, representing the average injection rate per unit of the fluctuation segment of the i-th irrigation segment; is the average value; The upper and lower limits of the confidence interval are: Wherein, n is the number of irrigation segments in the set of similar irrigation segments; s is the standard deviation; z is the z score of the standard normal distribution; is the average value; The injection rate interval is: Among them, Q max is the upper limit of injection rate, Q min is the lower limit of injection rate.

7. An injection rate interval analysis device for curtain grouting at the toe of a pumped storage project, characterized in that: include: An acquisition module, used for acquiring grouting data of a grouting site and adding the grouting data to a historical data set; A detection module, for detecting the type of the injection rate duration process line of each irrigation section based on the injection rate duration data of each irrigation section in the historical data set, and screening out irrigation sections that meet preset normal conditions from the types of the injection rate duration process lines, and obtaining standard irrigation sections according to the irrigation sections that meet the preset normal conditions; A determination module, for determining the fluctuation range of the injection rate duration process line of the standard irrigation section based on the injection rate inflection point time of the standard irrigation section, and calculating the fluctuation unit average injection rate of the standard irrigation section within the fluctuation range; A calculation module, used for calculating a grouting similarity index of the standard grouting section based on the grouting parameters of the standard grouting section; A clustering module, used for obtaining the grouting similarity mode of the standard grouting section according to the grouting similarity index, and clustering the standard grouting sections that meet the preset same grouting similarity mode conditions to obtain at least one similar grouting section set; A removal module, used for transmitting the unit average injection rate of the fluctuation section in the at least one similar irrigation section set to the parameter array of the corresponding similar irrigation section set, and removing outliers in the parameter array to obtain a parameter array after removing outliers; A judgment module is used to perform a normality test on the parameter array after removing outliers, obtain a normal test result, and judge whether the normal test result meets a preset qualification condition; An analysis module is used to calculate, when the normal test result satisfies the preset qualification condition, the confidence interval of the unit average injection rate of the fluctuation section when the confidence level of the similar irrigation section set meets the preset percentage condition according to the normal distribution, and determine the injection rate interval when the curtain grouting at the toe of the upper reservoir dam of the pumped storage project is in the grouting similarity mode according to the confidence interval.

8. An electronic device, characterized in that: include: A memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the program to implement the injection rate interval analysis method for curtain grouting at the toe of a reservoir dam in a pumped storage project as described in any one of claims 1 to 6.

9. A computer-readable storage medium having a computer program stored thereon, characterized in that: The program is executed by a processor to implement the injection rate interval analysis method for curtain grouting of the upper reservoir dam toe of a pumped storage project as described in any one of claims 1 to 6.

10. A computer program product, comprising a computer program, characterized in that The computer program is executed to implement the injection rate interval analysis method for curtain grouting of the upper reservoir dam toe of a pumped storage project as described in any one of claims 1 to 6.

Images