Boundary condition conversion method and device for drain hole group

Abstract

The invention relates to the technical field of water conservancy and hydropower engineering, in particular to a drain hole group boundary condition conversion method and device.The method comprises the steps that a drain hole group is divided to obtain all drain downward holes in the drain hole group; based on each drainage pitching hole, obtaining an initial boundary as a water head boundary; obtaining a seepage field simulation result based on the initial boundary, and obtaining a water head value and / or a flow value of the hole wall of the drainage downward hole according to the seepage field simulation result; judging whether the flow value meets the drainage condition of the water head boundary or not; if the flow value does not meet the drainage condition of the water head boundary, the water head boundary is converted until iterative convergence is achieved, and the actual boundary of the drainage hole group meeting the drainage condition and the seepage field calculation result are obtained. Wherein the connotation of the water head-potential overflow mixed boundary is creatively expanded particularly. Therefore, the problems that the essence of the drainage hole cannot be reflected, the simulation result cannot meet the engineering precision requirement, the error rate is high, the numerical stability is poor, convergence is not easy, and the performance evaluation and optimization design requirements of an engineering anti-seepage drainage system cannot be met in the related technology are solved.

Description

Technical Field

[0001] This application relates to the field of water conservancy and hydropower engineering technology, and in particular to a method and apparatus for transforming the boundary conditions of a drainage hole group. Background Technology

[0002] In numerous engineering fields such as water conservancy, hydropower, civil engineering, mining, and agriculture, drainage is the most important and cost-effective engineering measure for controlling groundwater movement. Especially in large-scale water conservancy and hydropower projects, drainage plays a crucial role in reducing uplift pressure at dam foundations, decreasing pore water pressure in tunnel surrounding rock, and improving the seepage stability of rock masses. Detailed simulation of drainage well systems is a prerequisite for performance evaluation and optimized design; however, due to the characteristics of drainage well systems in terms of geometry, quantity, and spatial arrangement, detailed simulation is extremely difficult.

[0003] In related technologies, numerical simulation methods for drainage hole groups can be divided into four categories: The first category treats the function of drainage holes as a sink term in the seepage model, with its drainage volume proportional to the difference between the aquifer head and the drainage hole elevation; the second category treats the flow in the drainage holes as pipe flow and uses an equivalent permeability coefficient to simulate the function of the drainage holes; the third category combines a steady flow model of groundwater flowing to a well and uses a semi-analytical method to simplify the calculation of drainage volume of drainage holes or the distribution of head near them; the fourth category simulates drainage holes by setting reasonable head or flow boundary conditions and iterating repeatedly.

[0004] However, among the related technologies, the first three types of methods cannot reflect the essence of drainage holes, the simulation results are difficult to meet the engineering accuracy requirements, and the error rate is high. The fourth type of method has poor numerical stability and is not easy to converge, which cannot meet the performance evaluation and optimization design requirements of engineering seepage prevention and drainage systems, and urgently needs to be improved. Summary of the Invention

[0005] This application provides a method and apparatus for transforming boundary conditions of a drainage hole group to solve the problems in related technologies, such as the inability to reflect the essence of drainage holes, the difficulty in meeting engineering accuracy requirements in simulation results, the high error rate, poor numerical stability, and difficulty in convergence, which cannot meet the performance evaluation and optimization design needs of engineering seepage prevention and drainage systems.

[0006] The first aspect of this application provides a method for transforming the boundary conditions of a drainage hole group, comprising the following steps: dividing the drainage hole group to obtain all drainage downholes in the drainage hole group; obtaining an initial boundary of any drainage downhole as a head boundary based on any one drainage downhole; obtaining a seepage field simulation result of the any one drainage downhole based on the initial boundary, and obtaining the head value and / or flow rate value of the hole wall of the any one drainage downhole according to the seepage field simulation result; determining whether the flow rate value satisfies the drainage conditions of the head boundary; if the flow rate value does not satisfy the drainage conditions of the head boundary, transforming the head boundary until iterative convergence is obtained to obtain the actual boundary of the drainage hole group that satisfies the drainage conditions, so as to obtain the seepage field calculation result of the drainage hole group according to the actual boundary, wherein the seepage field calculation result includes at least one of head distribution, seepage flow rate, and hydraulic gradient.

[0007] Optionally, in one embodiment of this application, obtaining the initial boundary of any one drainage hole as a head boundary based on any one drainage hole includes: obtaining the orifice elevation or the corridor floor elevation of any one drainage hole; determining whether the head value of any one drainage hole is equal to the orifice elevation or the corridor floor elevation; if the head value is equal to the orifice elevation or the corridor floor elevation, then the initial boundary of any one drainage hole is obtained, and the initial boundary is used as the head boundary.

[0008] Optionally, in one embodiment of this application, the step of converting the head boundary if the flow rate value does not meet the drainage conditions of the head boundary includes: if the actual unit area normal flow rate of any one of the drainage holes does not meet the drainage conditions of the preset unit area normal flow rate of the head boundary, then the head boundary is converted into a water-proof boundary, and it is determined whether the actual boundary of any one of the drainage holes is a water-proof boundary to obtain a first determination result; if the actual total flow rate of any one of the drainage holes does not meet the drainage conditions of the preset total flow rate of the head boundary, then the head boundary is converted into a head-potential overflow mixing boundary, and it is determined whether the actual boundary of any one of the drainage holes is a head-potential overflow mixing boundary to obtain a second determination result.

[0009] Optionally, in one embodiment of this application, the step of converting the head boundary if the flow rate value does not meet the drainage conditions of the head boundary further includes: obtaining the actual head of any one of the drainage holes based on the first determination result; determining that the actual head of any one of the drainage holes is greater than or equal to the target head; and if the actual head of any one of the drainage holes is greater than or equal to the target head, converting the actual boundary of any one of the drainage holes from a water-blocking boundary to the head boundary.

[0010] Optionally, in one embodiment of this application, the step of converting the head boundary if the flow rate value does not meet the preset condition of the head boundary further includes: obtaining the actual head of any one of the drainage holes based on the second judgment result; determining whether the actual head of any one of the drainage holes is greater than or equal to the target head; if the actual head of any one of the drainage holes is greater than or equal to the target head, then converting the actual boundary of any one of the drainage holes from the head-potential overflow mixing boundary to the head boundary.

[0011] Optionally, in one embodiment of this application, the head-potential overflow mixing boundary is a flow self-balancing state within a single drainage orifice where there is no overflow at the orifice, i.e., the inflow rate at the upper part of the sidewall of the single drainage orifice is equal to the outflow rate at the lower part. The relationship between the head and flow rate on the sidewall of the drainage orifice can be, but is not limited to, the following:

[0012]

[0013] Where h is the water head on the wall of the drainage hole, z is the elevation of the node on the wall of the drainage hole, and q n Г represents the normal flow rate per unit area on the wall surface of the drainage hole. upper The boundary of the wall above the water level inside the drainage hole, Г lower The wall boundary below the water level inside the drainage hole. Let be the head at the head boundary, and Q be the total flow rate of the drainage orifice. Flow rate is negative for outflow and positive for inflow.

[0014] A second aspect of this application provides a boundary condition transformation device for a drainage hole group, comprising: an acquisition module for dividing the drainage hole group to acquire all drainage holes in the drainage hole group; a first generation module for obtaining an initial boundary of any drainage hole as a head boundary based on any drainage hole; a second generation module for obtaining a seepage field simulation result of any drainage hole based on the initial boundary, and obtaining the head value and / or flow rate value of the hole wall of any drainage hole according to the seepage field simulation result; a judgment module for judging whether the flow rate value satisfies the drainage condition of the head boundary; and a transformation module for transforming the head boundary when the flow rate value does not satisfy the drainage condition of the head boundary, until iterative convergence is obtained to obtain the actual boundary of the drainage hole group that satisfies the drainage condition, so as to obtain the seepage field calculation result of the drainage hole group according to the actual boundary, wherein the seepage field calculation result includes at least one of head distribution, seepage flow rate, and hydraulic gradient.

[0015] Optionally, in one embodiment of this application, the first generation module includes: a first acquisition unit, configured to acquire the orifice elevation or the corridor floor elevation of any one of the drainage holes; a first judgment unit, configured to judge whether the water head value of any one of the drainage holes is equal to the orifice elevation or the corridor floor elevation; and a generation unit, configured to obtain the initial boundary of any one of the drainage holes when the water head value is equal to the orifice elevation or the corridor floor elevation connected thereto, and to use the initial boundary as the water head boundary.

[0016] Optionally, in one embodiment of this application, the conversion module includes: a first conversion unit, configured to convert the head boundary into a water-proof boundary when the actual unit area normal flow rate of any one of the drainage holes does not meet the drainage condition of the preset unit area normal flow rate of the head boundary, and determine whether the actual boundary of any one of the drainage holes is a water-proof boundary, thereby obtaining a first determination result; and a second conversion unit, configured to convert the head boundary into a head-potential overflow mixing boundary when the actual total flow rate of any one of the drainage holes does not meet the drainage condition of the preset total flow rate of the head boundary, and determine whether the actual boundary of any one of the drainage holes is a head-potential overflow mixing boundary, thereby obtaining a second determination result.

[0017] Optionally, in one embodiment of this application, the conversion module further includes: a second acquisition unit, configured to acquire the actual water head of any one of the drainage holes based on the first judgment result; a second judgment unit, configured to determine that the actual water head of any one of the drainage holes is greater than or equal to the target water head; and a third conversion unit, configured to convert the actual boundary of any one of the drainage holes from a water-proof boundary to the water head boundary when the actual water head of any one of the drainage holes is greater than or equal to the target water head.

[0018] Optionally, in one embodiment of this application, the conversion module further includes: a third acquisition unit, configured to acquire the actual head of any one of the drainage holes based on the second judgment result; a third judgment unit, configured to determine whether the actual head of any one of the drainage holes is greater than or equal to the target head; and a fourth conversion unit, configured to convert the actual boundary of any one of the drainage holes from a head-potential overflow mixing boundary to a head boundary when the actual head of any one of the drainage holes is greater than or equal to the target head.

[0019] Optionally, in one embodiment of this application, the head-potential overflow mixing boundary is a flow self-balancing state within a single drainage orifice where there is no overflow at the orifice, i.e., the inflow rate at the upper part of the sidewall of the single drainage orifice is equal to the outflow rate at the lower part. The relationship between the head and flow rate on the sidewall of the drainage orifice can be, but is not limited to, the following:

[0020]

[0021] Where h is the water head on the wall of the drainage hole, z is the elevation of the node on the wall of the drainage hole, and q n Г represents the normal flow rate per unit area on the wall surface of the drainage hole. upper The boundary of the wall above the water level inside the drainage hole, Г lower The wall boundary below the water level inside the drainage hole. Let Q be the head at the head boundary, and let Q be the total flow rate of the drainage orifice, where the flow rate is negative for outflow and positive for inflow.

[0022] A third aspect of this 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 boundary condition transformation method for a group of drainage holes as described in the above embodiments.

[0023] A fourth aspect of this application provides a computer-readable storage medium storing a computer program that, when executed by a processor, implements the boundary condition transformation method for the drainage hole group described above.

[0024] A fifth aspect of this application provides a computer program product, including a computer program that, when executed, implements the boundary condition transformation method for the drainage hole group as described above.

[0025] This application embodiment can obtain the head value and / or flow rate of any drainage hole wall using seepage field simulation results when the initial boundary of the drainage hole is the head boundary. It then determines the relationship between the flow rate and the drainage conditions of the head boundary, transforms the head boundary, and iteratively converges to obtain the actual boundary of the drainage hole group that satisfies the drainage conditions. Based on the actual boundary, the seepage field calculation results of the drainage hole group are obtained. This achieves the application of reasonable boundary conditions to drainage holes under complex conditions based on a strict boundary condition transformation algorithm, thereby avoiding the functional errors in drainage hole simulation caused by related technologies. It can efficiently and accurately simulate the groundwater seepage field containing a three-dimensional drainage hole group. Therefore, it solves the problems in related technologies, such as the inability to reflect the essence of drainage holes, the difficulty in meeting engineering accuracy requirements, high error rate, poor numerical stability, difficulty in convergence, and inability to meet the performance evaluation and optimization design needs of engineering seepage prevention and drainage systems.

[0026] Additional aspects and advantages of this application will be set forth in part in the description which follows, and in part will be obvious from the description, or may be learned by practice of this application. Attached Figure Description

[0027] The above and / or additional aspects and advantages of this application will become apparent and readily understood from the following description of the embodiments taken in conjunction with the accompanying drawings, wherein:

[0028] Figure 1 This is a flowchart illustrating a method for transforming boundary conditions for a group of drainage holes according to an embodiment of this application;

[0029] Figure 2 This is a schematic diagram of a model provided according to an embodiment of this application;

[0030] Figure 3 This is a comparative schematic diagram showing the groundwater flow characteristics when the drainage borehole provided according to an embodiment of this application is configured with different boundary conditions;

[0031] Figure 4 A block diagram illustrating the transformation relationship of the boundary conditions of a drainage hole according to an embodiment of this application;

[0032] Figure 5 This is a flowchart illustrating the working principle of a boundary condition transformation method for a drainage hole group according to an embodiment of this application;

[0033] Figure 6 This is a block diagram of a boundary condition transformation device for a group of drainage holes provided according to an embodiment of this application;

[0034] Figure 7 This is a schematic diagram of the structure of an electronic device provided according to an embodiment of this application. Detailed Implementation

[0035] The embodiments of this application are described in detail below. Examples of these embodiments are shown in the accompanying drawings, wherein the same or similar reference numerals denote the same or similar elements or elements having the same or similar functions throughout. The embodiments described below with reference to the accompanying drawings are exemplary and intended to explain this application, and should not be construed as limiting this application.

[0036] The following describes a method and apparatus for transforming boundary conditions of a drainage hole group according to embodiments of this application, with reference to the accompanying drawings. Addressing the problems mentioned in the background art, such as the inability to reflect the essence of drainage holes, simulation results failing to meet engineering accuracy requirements, high error rates, poor numerical stability, and difficulty in convergence, thus failing to meet the performance evaluation and optimization design needs of engineering seepage prevention and drainage systems, this application provides a method for transforming boundary conditions of a drainage hole group. In this method, the head value and / or flow rate value of any drainage hole wall can be obtained using seepage field simulation results, given that the initial boundary of the drainage hole is a head boundary. The relationship between the flow rate value and the drainage conditions of the head boundary is determined, and the head boundary is transformed until iterative convergence is achieved, obtaining the actual boundary of the drainage hole group that meets the drainage conditions. The seepage field calculation results of the drainage hole group are then obtained based on the actual boundary. This achieves the application of reasonable boundary conditions to drainage holes under complex conditions based on a strict boundary condition transformation algorithm, thereby avoiding functional errors in drainage hole simulation caused by related technologies and enabling efficient and accurate simulation of groundwater seepage fields containing three-dimensional drainage hole groups. This solves the problems in related technologies, such as the inability to reflect the essence of drainage holes, the simulation results failing to meet engineering accuracy requirements, the high error rate, poor numerical stability, difficulty in convergence, and the inability to meet the performance evaluation and optimization design needs of engineering seepage prevention and drainage systems.

[0037] Specifically, Figure 1 This is a flowchart of a method for transforming boundary conditions of a drainage hole group according to an embodiment of this application.

[0038] like Figure 1 As shown, the boundary condition transformation method for this drainage hole group includes the following steps:

[0039] In step S101, the drainage hole group is divided to obtain all the drainage holes in the drainage hole group.

[0040] It is understood that, based on the arrangement of the drainage holes, the drainage hole group can be divided into two categories: drainage holes drilled upwards and drainage holes drilled downwards. Drainage holes drilled upwards can be understood as drainage holes drilled upwards from the top slab or sidewall of the drainage gallery or drainage holes drilled between two drainage galleries at different elevations; drainage holes drilled downwards can be understood as drainage holes drilled downwards into the bedrock from the bottom slab of the drainage gallery.

[0041] Furthermore, in this embodiment of the application, a potential overflow boundary condition is set for each drainage hole in the drainage hole group, and this condition remains unchanged in subsequent steps.

[0042] As one possible implementation method, embodiments of this application can divide the drainage hole group and obtain all the drainage holes in the drainage hole group.

[0043] In step S102, based on any one drainage hole, the initial boundary of any one drainage hole is obtained as the head boundary.

[0044] In actual implementation, for any one of the drainage holes in the drainage hole group, the initial boundary of each drainage hole is set as a head boundary.

[0045] Optionally, in one embodiment of this application, obtaining the initial boundary of any drainage hole as the head boundary based on any drainage hole includes: obtaining the orifice elevation or the floor elevation of any drainage hole; determining whether the head value of any drainage hole is equal to the orifice elevation or the floor elevation of the corridor; if the head value is equal to the orifice elevation or the floor elevation of the corridor, then the initial boundary of any drainage hole is obtained, and the initial boundary is used as the head boundary.

[0046] It is understandable that the water head boundary can be understood as the water head value of the drainage hole being equal to the orifice elevation of the corresponding drainage hole or the floor elevation of the corridor connected to it.

[0047] As one possible implementation method, embodiments of this application can obtain the orifice elevation or the floor elevation of any drainage hole, and then determine whether the water head value of any drainage hole is equal to the orifice elevation or the floor elevation of the corridor. If the water head value is equal to the orifice elevation or the floor elevation of the corridor connected to it, the water head boundary is used as the initial boundary.

[0048] For example, the model diagram used in the embodiments of this application is as follows: Figure 2 As shown, the model in this embodiment measures 100×27×140m, with a top elevation of 130m and a bottom elevation of -10m. Two drainage corridors, horizontally spaced 5m apart and at elevations of 100m and 55m respectively, are constructed within the stratum. Each corridor has a cross-sectional area of ​​3×3m². The bottom of both corridors has drainage holes 40m long, spaced 3m apart, and with a diameter of 15cm. The permeability coefficient of the stratum rock mass is 1×10⁻⁶. -5 cm / s, the model has a water head boundary of 120m on the left and right sides, the front, back and bottom surfaces of the model are assumed to be water-proof boundaries, the walls of the drainage corridor are all potential overflow boundaries, the lower drainage hole is set to a water head boundary of 55m (i.e. the bottom plate elevation of the lower corridor), and the upper drainage hole is set to a water head boundary of 100m (i.e. the bottom plate elevation of the upper corridor).

[0049] In step S103, the seepage field simulation results of any drainage hole are obtained based on the initial boundary, and the head value and / or flow rate value of the hole wall of any drainage hole are obtained according to the seepage field simulation results.

[0050] In some embodiments, the present application embodiments can simulate the seepage field based on the initial boundary set by the drainage hole group, obtain the seepage field simulation results, and obtain the head value and / or flow rate value of each drainage hole wall according to the seepage field simulation results.

[0051] Furthermore, in the embodiments of this application, the head value of each drainage hole wall may include, but is not limited to, the node head and head distribution on each drainage hole wall, and this application does not impose specific limitations.

[0052] The flow rate value of each drainage hole wall in the embodiments of this application may include, but is not limited to, the node normal flow rate, the total flow rate, and the infiltration and outflow rates of each drainage hole wall. This application does not impose specific limitations.

[0053] For example, in this embodiment of the application, the initial boundary of the upper drainage hole of the model is set, and then the seepage field simulation results are output. Figure 3 This diagram illustrates the groundwater flow characteristics when the drainage boreholes are transformed into different boundaries. The left side shows the hydraulic head distribution within the formation and the flow distribution on the borehole wall, while the right side shows the hydraulic head distribution in the upper drainage borehole curtain section. When the initial boundaries of the upper drainage boreholes are all set to a 100m hydraulic head boundary, such as... Figure 3 As shown in (a), the infiltration flow rate of each drainage hole is Q. + =20.3cm 3 / s, the outflow rate is Q - = -3.0cm 3 / s, the total flow rate is Q = 17.3cm 3 / s.

[0054] In step S104, it is determined whether the flow rate value meets the drainage conditions of the head boundary.

[0055] As one possible implementation method, embodiments of this application can determine whether each drainage hole meets the drainage conditions of the head boundary based on the flow rate value of any drainage hole wall.

[0056] In step S105, if the flow rate does not meet the drainage conditions of the head boundary, the head boundary is transformed until the iteration converges to obtain the actual boundary of the drainage hole group that meets the drainage conditions, so as to obtain the seepage field calculation results of the drainage hole group based on the actual boundary. The seepage field calculation results include at least one of head distribution, seepage flow rate and hydraulic gradient.

[0057] In practical implementation, embodiments of this application can transform the head boundary for drainage conditions where the flow rate does not meet the head boundary, until iterative convergence, to obtain the actual boundary of the drainage hole group that meets the drainage conditions, and obtain the seepage field calculation results of the drainage hole group based on the actual boundary. The drainage conditions can be set by those skilled in the art according to actual conditions, and this application does not impose specific limitations.

[0058] Furthermore, in the embodiments of this application, the calculation results of the seepage field may include, but are not limited to, the following: when the boundary conditions of the drainage hole and the seepage field converge iteratively, the calculation results of the head distribution, seepage flow rate and hydraulic gradient of the seepage field are output as needed. This application does not impose specific limitations.

[0059] For example, embodiments of this application can yield the following results: Figure 3 (c) shows the calculated results of the head distribution and seepage flow rate of the seepage field, where, by Figure 3 (c) It can be seen that when the upper drainage hole is a water head-potential overflow mixed boundary, the calculated groundwater level is much lower than the water head boundary and the impermeable boundary under this boundary condition.

[0060] Optionally, in one embodiment of this application, if the flow rate does not meet the drainage conditions of the head boundary, the head boundary is transformed, including: if the actual unit area normal flow rate of any drainage orifice does not meet the preset unit area normal flow rate drainage conditions of the head boundary, the head boundary is transformed into a water-proof boundary, and it is determined whether the actual boundary of any drainage orifice is a water-proof boundary to obtain a first determination result; if the actual total flow rate of any drainage orifice does not meet the preset total flow rate drainage conditions of the head boundary, the head boundary is transformed into a head-potential overflow mixed boundary, and it is determined whether the actual boundary of any drainage orifice is a head-potential overflow mixed boundary to obtain a second determination result.

[0061] It is understood that, in the embodiments of this application, when all drainage holes are installed in the horizontal gallery with the lowest elevation (i.e., seepage can only overflow from the drainage hole openings), if the initial boundary of the drainage hole is set as the head boundary, the normal flow rate per unit area of ​​the nodes on the drainage hole wall is less than zero (wherein, in the embodiments of this application, it can be specified that outflow is negative and inflow is positive, and the specific settings can be made by those skilled in the art according to the actual situation, and this application does not impose specific restrictions).

[0062] In some embodiments, when the actual normal flow rate per unit area of ​​any drainage orifice does not meet the drainage condition of a certain normal flow rate per unit area of ​​the head boundary, such as when the actual normal flow rate per unit area of ​​a node on the wall of any drainage orifice is greater than or equal to zero, the head boundary is transformed into a water-proof boundary, and it is determined whether the actual boundary of any drainage orifice is a water-proof boundary, thus obtaining a first determination result. The drainage condition of a certain normal flow rate per unit area can be set by those skilled in the art according to actual conditions, and this application does not impose specific limitations.

[0063] For example, in the embodiments of this application, when all drainage holes are installed in the lowest-elevation horizontal corridor, i.e., seepage can only overflow from the drainage hole openings, if the drainage holes are set as head boundaries, the normal flow rate q per unit area on the drainage hole wall should be satisfied. n <0 (wherein, in this embodiment of the application, outflow is defined as negative and inflow as positive, but this application does not impose specific limitations). If there is a node with normal flow rate q on the orifice wall. n If the value is ≥0, then the water head boundary should be transformed into a watertight boundary.

[0064] Furthermore, it should be noted that in the embodiments of this application, when the drainage holes are installed in an inclined stepped corridor or a multi-level three-dimensional corridor (i.e., not only the hole opening is the seepage overflow outlet, but other drainage structures near the bottom of the hole may also constitute its potential discharge outlet), if the initial boundary of any drainage hole is set as the head boundary, then the total flow rate of the corresponding drainage hole should be less than zero.

[0065] In some embodiments, when the actual total flow rate of any drainage orifice does not meet the drainage condition of a certain total flow rate at the head boundary (e.g., the actual total flow rate of the drainage orifice is greater than or equal to zero), the head boundary is transformed into a head-potential overflow mixing boundary, and it is determined whether the actual boundary of any drainage orifice is a head-potential overflow mixing boundary, thus obtaining a second determination result. The drainage condition of a certain total flow rate can be set by those skilled in the art according to actual conditions, and this application does not impose specific limitations.

[0066] For example, in the embodiments of this application, when the drainage hole is installed in an inclined stepped corridor or a multi-level three-dimensional corridor, that is, not only is the hole opening a seepage overflow outlet, but other drainage structures near the bottom of the hole may also constitute its potential discharge outlet, if the drainage hole is set as a head boundary, the total flow rate of the drainage hole should be satisfied. Among them, Г dra This represents the boundary of the drainage orifice wall. If the total flow rate Q on the orifice wall is greater than or equal to 0, then the head boundary should be transformed into a head-potential overflow mixing boundary.

[0067] Optionally, in one embodiment of this application, if the flow rate does not meet the drainage conditions of the head boundary, the head boundary is transformed, which further includes: obtaining the actual head of any drainage hole based on the first judgment result; determining that the actual head of any drainage hole is greater than or equal to the target head; if the actual head of any drainage hole is greater than or equal to the target head, then the actual boundary of any drainage hole is transformed from a water-blocking boundary to a head boundary.

[0068] As one possible implementation method, embodiments of this application can obtain the actual water head of any drainage hole based on the first judgment result, and if the actual water head of any drainage hole is greater than or equal to the target water head, the actual boundary of any drainage hole can be transformed from a water-proof boundary to a water head boundary.

[0069] In this embodiment, the watertight boundary can be understood as follows: when all drainage holes are installed in the lowest-elevation horizontal corridor (i.e., seepage can only overflow from the drainage hole openings), the nodal water head on the drainage hole wall is less than the target water head (wherein, the target water head can be understood as the elevation of the drainage hole opening or the elevation of the corridor floor slab connected to it, which can be specifically set by those skilled in the art according to the actual situation, and this application does not impose specific limitations), which can be expressed as the actual water head h on the drainage hole wall. <z top , where z top This refers to the elevation of the drainage borehole opening or the elevation of the floor slab of the corridor connected to it.

[0070] Furthermore, when the node water head on the wall of the drainage hole in this embodiment is greater than or equal to the elevation of the drainage hole opening or the elevation of the corridor floor slab connected to it, the water-proof boundary should be transformed into a water head boundary.

[0071] For example, in the embodiments of this application, when all drainage holes are installed in the lowest-elevation horizontal corridor, i.e., seepage can only overflow from the drainage hole openings, if the drainage holes are set as a water-proof boundary, then h should be satisfied. <z top If h≥z top Therefore, in this embodiment of the application, the water-proof boundary is transformed into a water head boundary.

[0072] Specifically, in the embodiments of this application, when the boundary conditions of the drainage hole are transformed into a watertight boundary, such as Figure 3 (b) As shown in the drainage hole head distribution diagram on the right, the head on the wall surface near the drainage hole opening is higher than 100m (i.e., the elevation of the bottom slab of the upper corridor connected to it), which does not meet the head requirement of the waterproof boundary. Furthermore, embodiments of this application can be based on... Figure 4 The boundary condition transformation relationship of the drainage downhole shown is that the boundary condition of the upper drainage downhole should be transformed from a watertight boundary to a head boundary (and then automatically transformed into a head-potential overflow mixed boundary).

[0073] Optionally, in one embodiment of this application, if the flow rate does not meet the drainage conditions of the head boundary, the head boundary is transformed, which further includes: obtaining the actual head of any drainage orifice based on the second judgment result; determining whether the actual head of any drainage orifice is greater than or equal to the target head; if the actual head of any drainage orifice is greater than or equal to the target head, then the actual boundary of any drainage orifice is transformed from the head-potential overflow mixing boundary into the head boundary.

[0074] Those skilled in the art will understand that, based on the second judgment result, the embodiments of this application can obtain the actual head of any drainage hole, and if the actual head of any drainage hole is greater than or equal to the target head, the actual boundary of any drainage hole can be transformed from the head-potential overflow mixing boundary into the head boundary.

[0075] In this embodiment, the head-potential overflow mixing boundary can be understood as follows: when the drainage borehole is installed in an inclined stepped gallery or a multi-level three-dimensional gallery (i.e., not only the borehole opening is a seepage overflow outlet, but other drainage structures near the bottom of the borehole may also constitute its potential discharge outlet), the nodal head on the borehole wall is less than the elevation of the borehole opening or the elevation of the gallery floor plate connected to it. For example, it can be expressed as the actual head h on the borehole wall. <z top , where z top This refers to the elevation of the drainage borehole opening or the elevation of the floor slab of the corridor connected to it.

[0076] Furthermore, when the node water head on the wall of the drainage hole in this embodiment is greater than or equal to the elevation of the drainage hole orifice or the elevation of the floor slab of the corridor connected to it, the water head-potential overflow mixed boundary should be transformed into a water head boundary.

[0077] For example, in the embodiments of this application, when the drainage hole is installed in an inclined stepped corridor or a multi-level three-dimensional corridor, that is, not only the hole opening is a seepage overflow outlet, but other drainage structures near the bottom of the hole may also constitute its potential discharge outlet, if the drainage hole is set as a water head-potential overflow mixing boundary, it should meet the following requirements: <z top If h≥z top Therefore, in this embodiment of the application, the head-potential overflow mixing boundary is transformed into a head boundary.

[0078] Specifically, in the embodiments of this application, when the boundary conditions of the drainage orifice are transformed into a head-potential overflow mixed boundary, such as... Figure 3 As shown in (c), in this embodiment of the application, the infiltration flow rate of the drainage hole is strictly equal to the outflow flow rate, i.e., Q + =-Q - =8.7cm 3 / s, the water head distribution on the wall of the drainage hole also meets the requirements of the water head-potential overflow mixing boundary, therefore the boundary conditions of the drainage hole are set reasonably and correctly.

[0079] Optionally, in one embodiment of this application, the head-potential overflow mixing boundary is a flow self-balancing state within a single drainage orifice where there is no overflow at the orifice, i.e., the infiltration flow rate in the upper part of the sidewall of the single drainage orifice is equal to the outflow flow rate in the lower part. The relationship between the head and flow rate on the sidewall of the drainage orifice can be, but is not limited to, the following:

[0080]

[0081] Where h is the water head on the wall of the drainage hole, z is the elevation of the node on the wall of the drainage hole, and q n Г represents the normal flow rate per unit area on the wall surface of the drainage hole. upper The boundary of the wall above the water level inside the drainage hole, Г lower The wall boundary below the water level inside the drainage hole. Let Q be the head at the head boundary, and let Q be the total flow rate of the drainage orifice, where the flow rate is negative for outflow and positive for inflow.

[0082] As one possible implementation, the head-potential overflow mixing boundary in this application embodiment differs from the mixing boundary on a conventional drainage well with a known pumping flow rate. Instead, it corresponds to a flow self-balancing state within a single drainage orifice where there is no overflow at the orifice. That is, the inflow rate at the upper part of the sidewall of the single drainage orifice is equal to the outflow rate at the lower part. The relationship between the head and flow rate on the sidewall of the drainage orifice can be, but is not limited to, the following:

[0083]

[0084] Where h is the water head on the wall of the drainage hole, z is the elevation of the node on the wall of the drainage hole, and q n Г represents the normal flow rate per unit area on the wall surface of the drainage hole. upper The boundary of the wall above the water level inside the drainage hole, Г lower The wall boundary below the water level inside the drainage hole. Let be the head at the head boundary, and Q be the total flow rate of the drainage orifice. Here, flow rate is negative for outflow and positive for inflow.

[0085] The working principle of the boundary condition transformation method for the drainage hole group proposed in this application will be described in detail below with reference to a specific embodiment.

[0086] in, Figure 5 This is a flowchart illustrating the working principle of a boundary condition transformation method for a drainage hole group according to an embodiment of this application.

[0087] Step S501: Setting initial boundary conditions for drainage borehole.

[0088] In other words, in this embodiment of the application, the initial boundary of all drainage holes can be set as the head boundary, wherein the head boundary can be understood as the head of the drainage hole being equal to the elevation of its orifice or the elevation of the floor slab of the corridor connected to it.

[0089] Specifically, the model diagram used in the embodiments of this application is as follows: Figure 2 As shown, the model in this embodiment measures 100×27×140m, with a top elevation of 130m and a bottom elevation of -10m. Two drainage corridors, horizontally spaced 5m apart and at elevations of 100m and 55m respectively, are constructed within the stratum. Each corridor has a cross-sectional area of ​​3×3m². The bottom of both corridors has drainage holes 40m long, spaced 3m apart, and with a diameter of 15cm. The permeability coefficient of the stratum rock mass is 1×10⁻⁶. -5 cm / s, the model has a water head boundary of 120m on the left and right sides, the front, back and bottom surfaces of the model are assumed to be water-proof boundaries, the walls of the drainage corridor are all potential overflow boundaries, the lower drainage hole is set to a water head boundary of 55m (i.e. the bottom plate elevation of the lower corridor), and the upper drainage hole is set to a water head boundary of 100m (i.e. the bottom plate elevation of the upper corridor).

[0090] Step S502: Simulation of seepage field and calculation of seepage flow rate through drainage holes.

[0091] In other words, the embodiments of this application can set the initial boundary of the drainage hole, and then output the seepage field simulation results, so as to obtain the head distribution of each drainage hole and its infiltration and outflow rates.

[0092] Specifically, such as Figure 3 As shown in the embodiments of this application, the groundwater flow characteristics when the drainage borehole is transformed into different boundaries are illustrated. The left side shows the hydraulic head distribution within the formation and the flow distribution on the drainage borehole wall, while the right side shows the hydraulic head distribution in the upper drainage borehole curtain section. When the initial boundaries of the upper drainage boreholes are all set to a 100m hydraulic head boundary, as shown in the embodiments of this application, the groundwater flow characteristics when the drainage boreholes are transformed into different boundaries are illustrated in the embodiments of this application. Figure 3 As shown in (a), the infiltration flow rate of each drainage hole is Q. + =20.3cm 3 / s, the outflow rate is Q - = -3.0cm 3 / s, the total flow rate is Q = 17.3cm 3 / s.

[0093] Step S503: Determine whether any drainage orifice meets the head or flow rate of its set boundary conditions.

[0094] Specifically, when the boundary conditions of a drainage orifice are transformed into head boundaries, it can be determined whether any drainage orifice meets the head or flow rate of its set boundary conditions.

[0095] When all drainage holes are installed in the lowest-elevation horizontal gallery, if there is a node on the hole wall with a normal flow rate q n If the value is ≥0, then the water head boundary should be transformed into a watertight boundary.

[0096] When a drainage hole is installed in an inclined stepped gallery or a multi-level three-dimensional gallery, if the total flow rate Q on the hole wall is greater than or equal to 0, the water head boundary should be transformed into a water head-potential overflow mixed boundary.

[0097] Step S504: Adjust the boundary conditions of the drainage hole according to the transformation conditions.

[0098] In other words, the embodiments of this application can be based on Figure 4 The boundary condition transformation relationship of the drainage hole shown transforms the boundary of the drainage hole.

[0099] Specifically, when the upper drainage hole does not meet the flow rate requirement of the head boundary, but simultaneously meets the flow rate requirement towards the water-proof boundary (q) n The critical condition for the transformation of the head-potential overflow mixing boundary (Q≥0) or the head-potential overflow mixing boundary (Q≥0) will be discussed separately below:

[0100] When the boundary conditions of the drainage hole are transformed into a watertight boundary, such as Figure 3 (b) As shown in the drainage hole head distribution diagram on the right, the water head on the wall near the drainage hole opening is higher than 100m (i.e., the bottom elevation of the upper corridor connected to it), which does not meet the water head requirements of the water-proof boundary. Therefore, the boundary conditions of the upper drainage hole should be changed from the water-proof boundary to the water head boundary.

[0101] When the boundary conditions of the drainage orifice are transformed into a head-potential overflow mixed boundary, such as Figure 3 As shown in (c), in this embodiment of the application, the infiltration flow rate of the drainage hole is strictly equal to the outflow flow rate, i.e., Q + =-Q - =8.7cm 3 / s, the water head distribution on the wall of the drainage hole also meets the requirements of the water head-potential overflow mixing boundary, therefore the boundary conditions of the drainage hole are set reasonably and correctly.

[0102] Step S505: Output the seepage field results, such as head distribution, seepage flow rate, and seepage gradient.

[0103] Specifically, the embodiments of this application can achieve the following: Figure 3 (c) shows the calculated results of the head distribution and seepage flow rate of the seepage field, where, by Figure 3(c) It can be seen that when the upper drainage hole is a water head-potential overflow mixed boundary, the calculated groundwater level is much lower than the water head boundary and the impermeable boundary under this boundary condition.

[0104] The boundary condition transformation method for drainage hole groups proposed in this application can obtain the head value and / or flow rate value of any drainage hole wall node using seepage field simulation results when the initial boundary of the drainage hole is a head boundary. The relationship between the flow rate value and the drainage conditions of the head boundary is determined, and the head boundary is transformed until iterative convergence is achieved, resulting in the actual boundary of the drainage hole group that satisfies the drainage conditions. The seepage field calculation results of the drainage hole group are then obtained based on the actual boundary. This method applies reasonable boundary conditions to drainage holes under complex conditions based on a strict boundary condition transformation algorithm, thus avoiding functional errors in drainage hole simulation caused by related technologies. It can efficiently and accurately simulate the groundwater seepage field containing a three-dimensional drainage hole group. Therefore, it solves the problems in related technologies, such as the inability to reflect the essence of drainage holes, simulation results that fail to meet engineering accuracy requirements, high error rates, poor numerical stability, difficulty in convergence, and inability to meet the performance evaluation and optimization design needs of engineering seepage prevention and drainage systems.

[0105] Next, the boundary condition transformation device for the drainage hole group proposed according to the embodiments of this application is described with reference to the accompanying drawings.

[0106] Figure 6 This is a block diagram of a boundary condition transformation device for a group of drainage holes provided according to an embodiment of this application.

[0107] like Figure 6 As shown, the boundary condition conversion device 10 for the drainage hole group includes: an acquisition module 100, a first generation module 200, a second generation module 300, a judgment module 400, and a conversion module 500.

[0108] The acquisition module 100 is used to divide the drainage hole group into sections to acquire all the drainage holes in the drainage hole group.

[0109] The first generation module 200 is used to obtain the initial boundary of any drainage hole as the head boundary based on any drainage hole.

[0110] The second generation module 300 is used to obtain the seepage field simulation results of any drainage hole based on the initial boundary, and to obtain the head value and / or flow rate value of the hole wall of any drainage hole based on the seepage field simulation results.

[0111] The judgment module 400 is used to determine whether the flow rate value meets the drainage conditions of the head boundary.

[0112] The transformation module 500 is used to transform the head boundary when the flow rate does not meet the drainage conditions of the head boundary until the iteration converges to obtain the actual boundary of the drainage hole group that meets the drainage conditions, so as to obtain the seepage field calculation results of the drainage hole group based on the actual boundary. The seepage field calculation results include at least one of head distribution, seepage flow rate and hydraulic gradient.

[0113] Optionally, in one embodiment of this application, the first generation module 200 includes: a first acquisition unit, a first judgment unit, and a generation unit.

[0114] The first acquisition unit is used to acquire the elevation of the orifice of any drainage hole or the elevation of the floor slab of the corridor.

[0115] The first judgment unit is used to determine whether the water head value of any drainage hole is equal to the hole elevation or the corridor floor elevation.

[0116] The generation unit is used to obtain the initial boundary of any drainage hole when the head value is equal to the orifice elevation or the floor elevation of the corridor connected to it, and to use the initial boundary as the head boundary.

[0117] Optionally, in one embodiment of this application, the conversion module 500 includes: a first conversion unit and a second conversion unit.

[0118] The first conversion unit is used to convert the head boundary into a water-proof boundary when the actual unit area normal flow of any drainage hole does not meet the drainage condition of the preset unit area normal flow of the head boundary, and to determine whether the actual boundary of any drainage hole is a water-proof boundary, thereby obtaining the first judgment result.

[0119] The second conversion unit is used to convert the head boundary into a head-potential overflow mixing boundary when the actual total flow of any drainage hole does not meet the drainage condition of the preset total flow of the head boundary, and to determine whether the actual boundary of any drainage hole is a head-potential overflow mixing boundary, thereby obtaining a second judgment result.

[0120] Optionally, in one embodiment of this application, the conversion module 500 further includes: a second acquisition unit, a second judgment unit, and a third conversion unit.

[0121] The second acquisition unit is used to acquire the actual water head of any drainage hole based on the first judgment result.

[0122] The second judgment unit is used to determine whether the actual water head of any drainage hole is greater than or equal to the target water head.

[0123] The third conversion unit is used to convert the actual boundary of any drainage hole from a water-proof boundary to a head boundary when the actual head of any drainage hole is greater than or equal to the target head.

[0124] Optionally, in one embodiment of this application, the conversion module 500 further includes: a third acquisition unit, a third judgment unit, and a fourth conversion unit.

[0125] The third acquisition unit is used to acquire the actual water head of any drainage hole based on the second judgment result.

[0126] The third judgment unit is used to determine whether the actual water head of any drainage hole is greater than or equal to the target water head.

[0127] The fourth conversion unit is used to convert the actual boundary of any drainage hole from a head-potential overflow mixed boundary to a head boundary when the actual head of any drainage hole is greater than or equal to the target head.

[0128] Optionally, in one embodiment of this application, the head-potential overflow mixing boundary is a flow self-balancing state within a single drainage orifice where there is no overflow at the orifice, i.e., the infiltration flow rate in the upper part of the sidewall of the single drainage orifice is equal to the outflow flow rate in the lower part. The relationship between the head and flow rate on the sidewall of the drainage orifice can be, but is not limited to, the following:

[0129]

[0130] Where h is the water head on the wall of the drainage hole, z is the elevation of the node on the wall of the drainage hole, and q n Г represents the normal flow rate per unit area on the wall surface of the drainage hole. upper The boundary of the wall above the water level inside the drainage hole, Г lower The wall boundary below the water level inside the drainage hole. Let Q be the head at the head boundary, and let Q be the total flow rate of the drainage orifice, where the flow rate is negative for outflow and positive for inflow.

[0131] It should be noted that the explanation of the aforementioned embodiment of the boundary condition transformation method for the drainage hole group also applies to the boundary condition transformation device for the drainage hole group in this embodiment, and will not be repeated here.

[0132] The boundary condition transformation device for drainage hole groups proposed in this application can obtain the head value and / or flow rate value of any drainage hole wall node using seepage field simulation results when the initial boundary of the drainage hole is a head boundary. It then determines the relationship between the flow rate value and the drainage conditions of the head boundary, transforms the head boundary, and iteratively converges to obtain the actual boundary of the drainage hole group that satisfies the drainage conditions. Based on the actual boundary, the seepage field calculation results of the drainage hole group are obtained. This achieves the application of reasonable boundary conditions to drainage holes under complex conditions based on a strict boundary condition transformation algorithm, thereby avoiding functional errors in the simulation of drainage holes caused by related technologies. It can efficiently and accurately simulate the groundwater seepage field containing a three-dimensional drainage hole group. Therefore, it solves the problems in related technologies, such as the inability to reflect the essence of drainage holes, the difficulty in meeting engineering accuracy requirements, high error rate, poor numerical stability, difficulty in convergence, and inability to meet the performance evaluation and optimization design needs of engineering seepage prevention and drainage systems.

[0133] Figure 7 This is a schematic diagram of the structure of an electronic device provided according to an embodiment of this application. The electronic device may include:

[0134] The memory 701, the processor 702, and the computer program stored on the memory 701 and executable on the processor 702.

[0135] When the processor 702 executes the program, it implements the boundary condition transformation method for the drainage hole group provided in the above embodiments.

[0136] Furthermore, electronic devices also include:

[0137] Communication interface 703 is used for communication between memory 701 and processor 702.

[0138] The memory 701 is used to store computer programs that can run on the processor 702.

[0139] The memory 701 may include high-speed RAM memory, and may also include non-volatile memory, such as at least one disk storage device.

[0140] If the memory 701, processor 702, and communication interface 703 are implemented independently, then the communication interface 703, memory 701, and processor 702 can be interconnected via a bus to complete communication between them. The bus can be an Industry Standard Architecture (ISA) bus, a Peripheral Component Interconnect (PCI) bus, or an Extended Industry Standard Architecture (EISA) bus, etc. Buses can be categorized as address buses, data buses, control buses, etc. For ease of representation, Figure 7 The bus is represented by a single thick line, but this does not mean that there is only one bus or one type of bus.

[0141] Optionally, in a specific implementation, if the memory 701, processor 702, and communication interface 703 are integrated on a single chip, then the memory 701, processor 702, and communication interface 703 can communicate with each other through an internal interface.

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

[0143] This application also provides a computer-readable storage medium storing a computer program that, when executed by a processor, implements the above-described method for transforming the boundary conditions of a group of drainage holes.

[0144] This application also provides a computer program product, including a computer program that, when executed, implements the boundary condition transformation method for the drainage hole group as described above.

[0145] In the description of this specification, the references to terms such as "one embodiment," "some embodiments," "example," "specific example," or "some examples," etc., indicate that a specific feature, structure, material, or characteristic described in connection with that embodiment or example is included in at least one embodiment or example of this application. In this specification, the illustrative expressions of the above terms do not necessarily refer to the same embodiment or example. Furthermore, the specific features, structures, materials, or characteristics described may be combined in any suitable manner in one or more embodiments or examples. Moreover, without contradiction, those skilled in the art can combine and integrate the different embodiments or examples described in this specification, as well as the features of different embodiments or examples.

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

[0147] Any process or method described in the flowchart or otherwise herein can be understood as representing a module, segment, or portion of code comprising one or N executable instructions for implementing custom logic functions or processes, and the scope of the preferred embodiments of this application includes additional implementations in which functions may be performed not in the order shown or discussed, including substantially simultaneously or in reverse order depending on the functions involved, as should be understood by those skilled in the art to which embodiments of this application pertain.

[0148] The logic and / or steps represented in the flowchart or otherwise described herein, for example, can be considered as a sequenced list of executable instructions for implementing logical functions, and can be embodied in any computer-readable medium for use by, or in conjunction with, an instruction execution system, apparatus, or device (such as a computer-based system, a processor-included system, or other system that can fetch and execute instructions from, an instruction execution system, apparatus, or device). For the purposes of this specification, "computer-readable medium" can be any means that can contain, store, communicate, propagate, or transmit programs for use by, or in conjunction with, an instruction execution system, apparatus, or device. More specific examples (a non-exhaustive list) of computer-readable media include: an electrical connection having one or more wires (electronic device), a portable computer disk drive (magnetic device), random access memory (RAM), read-only memory (ROM), erasable and editable read-only memory (EPROM or flash memory), fiber optic devices, and portable optical disc read-only memory (CDROM). Alternatively, the computer-readable medium may be paper or other suitable media on which the program can be printed, since the program can be obtained electronically by optically scanning the paper or other medium, followed by editing, interpreting, or otherwise processing as necessary, and then stored in a computer memory.

[0149] It should be understood that the various parts of this application can be implemented using hardware, software, firmware, or a combination thereof. In the above embodiments, the N steps or methods can be implemented using software or firmware stored in memory and executed by a suitable instruction execution system. If implemented in hardware, as in another embodiment, it can be implemented using any one or more of the following techniques known in the art: discrete logic circuits having logic gates for implementing logical functions on data signals, application-specific integrated circuits (ASICs) having suitable combinational logic gates, programmable gate arrays (PGAs), field-programmable gate arrays (FPGAs), etc.

[0150] Those skilled in the art will understand that all or part of the steps of the methods described in the above embodiments can be implemented by a program instructing related hardware. The program can be stored in a computer-readable storage medium, and when executed, it includes one or a combination of the steps of the method embodiments.

[0151] Furthermore, the functional units in the various embodiments of this application can be integrated into a processing module, or each unit can exist physically separately, or two or more units can be integrated into a module. The integrated module can be implemented in hardware or as a software functional module. If the integrated module is implemented as a software functional module and sold or used as an independent product, it can also be stored in a computer-readable storage medium.

[0152] The storage medium mentioned above can be a read-only memory, a disk, or an optical disk, etc. Although embodiments of this application have been shown and described above, it is understood that the above embodiments are exemplary and should not be construed as limiting this application. Those skilled in the art can make changes, modifications, substitutions, and variations to the above embodiments within the scope of this application.

Claims

1. A method for transforming boundary conditions of a group of drainage holes, characterized in that, Includes the following steps: The drainage hole group is divided to obtain all the drainage holes in the drainage hole group; Based on any one drainage hole, the initial boundary of the arbitrary drainage hole is obtained as the head boundary; Based on the initial boundary, the seepage field simulation results of any one of the drainage holes are obtained, and the head value and / or flow rate value of the hole wall of any one of the drainage holes are obtained according to the seepage field simulation results. Determine whether the flow rate value meets the drainage conditions of the head boundary; If the flow rate value does not meet the drainage conditions of the head boundary, the head boundary is transformed until iterative convergence is obtained to obtain the actual boundary of the drainage hole group that meets the drainage conditions. The seepage field calculation result of the drainage hole group is obtained based on the actual boundary, wherein the seepage field calculation result includes at least one of head distribution, seepage flow rate and hydraulic gradient.

2. The boundary condition transformation method for the drainage hole group according to claim 1, characterized in that, The method of obtaining the initial boundary of any one drainage hole as the head boundary includes: Obtain the elevation of the orifice of any one of the drainage holes or the elevation of the corridor floor slab; Determine whether the head value of any one of the drainage holes is equal to the elevation of the hole or the elevation of the corridor floor. If the water head value is equal to the orifice elevation or the corridor floor elevation, then the initial boundary of any one of the drainage holes is obtained, and the initial boundary is taken as the water head boundary.

3. The boundary condition transformation method for the drainage hole group according to claim 1, characterized in that, If the flow rate value does not meet the drainage conditions of the head boundary, then the head boundary is transformed, including: If the actual unit area normal flow rate of any of the drainage holes does not meet the drainage condition of the preset unit area normal flow rate of the head boundary, then the head boundary is converted into a water-proof boundary, and it is determined whether the actual boundary of any of the drainage holes is a water-proof boundary, and a first judgment result is obtained. If the actual total flow rate of any one of the drainage holes does not meet the drainage condition of the preset total flow rate of the head boundary, the head boundary is transformed into a head-potential overflow mixing boundary, and it is determined whether the actual boundary of any one of the drainage holes is a head-potential overflow mixing boundary, thus obtaining a second determination result.

4. The boundary condition transformation method for the drainage hole group according to claim 3, characterized in that, The step of transforming the head boundary if the flow rate value does not meet the drainage conditions of the head boundary further includes: Based on the first judgment result, the actual water head of any one of the drainage holes is obtained; Determine that the actual head of any one of the drainage holes is greater than or equal to the target head; If the actual water head of any one of the drainage holes is greater than or equal to the target water head, then the actual boundary of any one of the drainage holes is transformed from the water-blocking boundary to the water head boundary.

5. The boundary condition transformation method for the drainage hole group according to claim 3, characterized in that, The step of transforming the head boundary if the flow rate value does not meet the drainage conditions of the head boundary further includes: Based on the second judgment result, the actual water head of any one of the drainage holes is obtained; Determine whether the actual water head of any one of the drainage holes is greater than or equal to the target water head; If the actual head of any one of the drainage holes is greater than or equal to the target head, then the actual boundary of any one of the drainage holes will be transformed from the head-potential overflow mixed boundary to the head boundary.

6. The boundary condition transformation method for the drainage hole group according to claim 1, characterized in that, The head-potential overflow mixing boundary is a self-balancing flow state within a single drainage orifice where there is no overflow at the orifice. Specifically, the inflow rate at the upper part of the orifice wall is equal to the outflow rate at the lower part. The relationship between the head and flow rate on the orifice wall is as follows: Where h is the water head on the wall of the drainage hole, z is the elevation of the node on the wall of the drainage hole, and q n Г represents the normal flow rate per unit area on the wall surface of the drainage hole. upper The boundary of the wall above the water level inside the drainage hole, Г lower The wall boundary below the water level inside the drainage hole. Let Q be the head at the head boundary, and let Q be the total flow rate of the drainage orifice, where the flow rate is negative for outflow and positive for inflow.

7. A boundary condition transformation device for a group of drainage holes, characterized in that, include: The acquisition module is used to divide the drainage hole group to acquire all the drainage holes in the drainage hole group; The first generation module is used to obtain the initial boundary of any one drainage hole as the head boundary. The second generation module is used to obtain the seepage field simulation result of any one of the drainage holes based on the initial boundary, and to obtain the head value and / or flow rate value of the hole wall of any one of the drainage holes based on the seepage field simulation result. The judgment module is used to determine whether the flow rate value meets the drainage conditions of the head boundary; The transformation module is used to transform the head boundary when the flow rate does not meet the drainage conditions of the head boundary until the iteration converges to obtain the actual boundary of the drainage hole group that meets the drainage conditions, so as to obtain the seepage field calculation result of the drainage hole group based on the actual boundary, wherein the seepage field calculation result includes at least one of head distribution, seepage flow rate and hydraulic gradient.

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, the processor executing the program to implement the boundary condition transformation method for a group of drainage holes as described in any one of claims 1-6.

9. A computer-readable storage medium having a computer program stored thereon, characterized in that, The program is executed by the processor to implement the boundary condition transformation method for the drainage hole group as described in any one of claims 1-6.

10. A computer program product, characterized in that, Includes a computer program, which, when executed, is used to implement the boundary condition transformation method for the drainage hole group as described in any one of claims 1-6.

Images