Disclosure of Invention
The invention aims to provide a multi-objective optimization method for group hole dredging and reducing of a low-permeability aquifer, which aims at high-pressure-bearing water dredging and reducing engineering of a coal bed bottom plate, starts from the low-permeability characteristic of the aquifer and the requirements of economy, environmental protection and the like, constructs the minimum water drainage hole under the condition of meeting the safe production condition, achieves the optimal dredging and reducing effect by the minimum single-hole flow, simultaneously reduces the total water drainage amount and the drilling engineering amount to the maximum extent, realizes the optimal balance among safety, economy and environmental benefits, reduces the construction cost and achieves the aim of safe exploitation; the method comprises the steps of selecting positions of scattered drilling holes in a roadway according to underground positions, generalizing the roadway where the candidate positions of the drilling holes are located into a straight line segment, searching the best point position meeting conditions on the straight line segment, and obtaining accurate hole position coordinates. Because the aquifer to be dredged and degraded needs to simultaneously meet the constraint conditions that the flow of a single hole is minimum to achieve the optimal dredging and degrading effect and the number of holes is minimum, the method adopts a multi-target particle swarm (MOPSO) algorithm, has the advantages of high searching speed, high efficiency and simple algorithm, and can quickly obtain the optimal solution of the model; the characteristics of large water pressure, poor permeability and nonuniform water-rich property commonly exist in deep mining limestone aquifers, single-hole water drainage is difficult to meet the requirements, and the number of drilled holes and the flow rate of the drilled holes are optimized in group-hole water drainage, so that the coal mining cost and the damage degree to the ecological environment are effectively reduced while the dredging target is achieved.
The purpose of the invention can be realized by the following technical scheme:
a low-permeability aquifer group-hole sparse-reduction multi-objective optimization method comprises the following steps:
s1: and determining the water level control of the dewatering and dewatering area according to the mining engineering layout.
S2: and collecting geological and hydrogeological data of a research area, determining the permeability of the aquifer to be dredged and lowered, determining parameters such as safe water level and depth of a control point, and calculating the maximum depth of the water drainage hole, namely recording the maximum depth as Smax when the water level of the water drainage hole is lowered to the top plate of the aquifer.
S3: and constructing a drainage drilling hole candidate position.
S4: and (3) setting the degree of the drainage hole penetrating through the aquifer, setting the constraint conditions of the hole number n, the single-hole flow Q and the depth S of the water level control point, and definitely establishing an optimization model by taking the minimum of Q and n as an optimization function.
S5: and (4) solving the optimal Q, n value of the model by using a multi-objective optimization method.
S6: the optimization model is discussed and the production practice is guided.
Preferably, in S1, the water level control point of the drainage and precipitation area is determined according to the mining engineering layout, and in the actual drainage engineering, the boundary of the coverage mining working face range is used as the drainage and precipitation boundary, and the working face end point is selected as the water level control point in the drainage and precipitation boundary.
Preferably, in S2, geological and hydrogeological data of the research area are collected, the permeability of the aquifer to be drained and lowered is determined, parameters such as the safe water level and the depth of the control point are determined, and the maximum depth Smax of the drain hole is calculated.
The safe water head value is a key index for dredging, the safe water head is the water head pressure born by the water-resisting layer of the coal bed bottom plate, the smaller the safe water head value is, the larger the water head required to be reduced is, if the safe water head value is larger than the actual water head value, the direct mining under pressure can be carried out without water drainage, and the safe water head of the working face is generally determined by adopting a water bursting coefficient method in the actual engineering, as shown in a formula (1).
p=TsM (1)
Wherein M is the thickness of a water-resisting layer, the unit M is, Ts is a critical water bursting coefficient, 0.06MPa/M is taken in a section with structural damage and 0.1MPa/M is taken in a complete section without fracture according to the water prevention and control requirements of a coal mine.
In the actual water discharge project, the boundary of the range of the covered mining working face is used as a dredging and descending boundary, the end point of the working face is selected as a water level control point in the dredging and descending boundary to calculate the safe water level, the water-proof bottom plate elevation Hd (namely the water-bearing layer top plate elevation) and the initial water level elevation H0 of the coal seam at each water level control point are determined, the safe water pressure is converted into the safe water level by using the formula Hs (p × 100+ Hd), and the safe water level descending depth at each water level control point is obtained by using the formula Ss (H0-Hs).
Preferably, the candidate positions of the water drainage drill holes are constructed in the step S3, the water drainage roadway needing to be provided with the water drainage holes is generalized to be a straight line segment, the inclination angle is θ, coordinates of two end points A, B of the water drainage roadway are respectively (xp1, yp1) and (xp2, yp2), a straight line equation of the roadway is obtained by using x and y of any two points on the straight line segment through a formula y ═ ax + b, a coordinate range xp1 ≤ x ≤ xp2 where x of the straight line segment is located is marked, a in the formula is a straight line slope, b is a straight line intercept, (xp1, xp2) is a horizontal coordinate range of the straight line segment, and the total length of the roadway is L.
N water drainage drill holes are distributed in the water drainage roadway, wherein the n water drainage drill holes are respectively F1, F2... Fi... Fn and are arranged at equal intervals, the interval between adjacent holes is D, and the coordinate of any water drainage hole is (xci, yci); m water level control points are arranged in the working surface, wherein the m water level control points are C1, C2 … Cj and Cm respectively, the coordinate of any control point is (xj and yj), the distance from any water discharging hole to any water level control point is recorded as rij, the water discharging quantity of all the water discharging holes is equal to Q, and the water discharging holes and the water level control points are shown in the attached drawing 1.
Preferably, in S4, constraint conditions of the degree of penetration of the drainage hole through the aquifer, the number n of holes, the single-hole flow Q, and the depth reduction S of the water level control point are set, and an optimization model is established by definitely taking the minimum of Q and n as an optimization function; the optimization of the grouped hole dredging and landing engineering of the limestone aquifer comprises the optimization of flow and the number of water drainage holes, aiming at the limestone aquifer with high bearing pressure and low permeability, the single-hole flow and the number of the water drainage holes are generally used as target functions in an optimization model, namely, the minimum single-hole flow and the number of the water drainage holes are reached under the constraint conditions of meeting the water level depth reduction and the like, the underground water resource is protected to the maximum extent, and the target functions can be expressed as follows:
Z=opt{Q,n} (2)
when the water is discharged from the incomplete well, and the confined aquifer where the incomplete well is located meets the assumed conditions used for deriving the Theis formula, when the linear distance r between the observation point and the water discharge hole is less than or equal to 1.5M, the depth reduction equation is as follows:
in the formula, Q is the water discharge amount of the drill hole; r is the linear distance between the calculated point and the water drainage drill hole; s is water level depth reduction; k is the permeability coefficient of the aquifer; m is the total aquifer thickness; u is the elastic water release coefficient; w (u) is a function of the entire Thas well; u is a well function argument; ζ is the non-integrity add-on coefficient of resistance; d is the distance from the top plate of the aquifer to the top of the filter of the water drainage hole, l is the distance from the top plate of the aquifer to the bottom of the filter of the dredging drill hole, and z is the distance from the top plate of the aquifer to the bottom of the water level observation hole, namely the opening position; t is the time for draining water.
From equation (3), it can be seen that the incomplete well dip is composed of two parts, the former representing the corresponding complete well dip and the latter representing the additional dip caused by the bending of the flow line near the extraction well due to the imperfections of the extraction well.
When the drainage hole completely penetrates through the thickness of the whole aquifer, namely the drilled hole is a complete well, or the linear distance r between an observation point and the center of the drainage hole is more than 1.5M, the additional resistance coefficient can be ignored, the additional resistance coefficient is simplified into a corresponding complete well formula, and a depth reduction equation is converted into:
when water is discharged from the group of holes, the water amount is respectively Q1, Q2 and … Qn, and the water head depth reduction at any point can be calculated according to the seepage superposition principle, namely the sum of the water head depth reduction of each single water discharge hole is equal to the sum of the water head depth reduction of each single water discharge hole, as shown in the formula.
Setting the left lower corner of the roadway as the arrangement position of a first water drainage hole (x1, y1), wherein the water drainage hole cannot exceed the range of the water drainage roadway, and according to the coordinates of the end points of the roadway and the linear equation of the roadway, the constraint condition of the number n of the water drainage holes of the group of holes can be written as follows:
when the water level of the water drainage hole is reduced to the top plate of the aquifer (marked as Smax), the water drainage flow of the drilling hole is marked as Qmax, the flow is the maximum constraint condition of the single-hole flow in optimization, namely the upper limit value of the drainage flow of the drilling hole, the radius of the water drainage hole is set as rw, and the flow constraint condition is that the water drainage hole is mutually interfered: when the drainage holes are complete wells or incomplete wells but the distance D between the holes is more than or equal to 1.5M, the drainage holes for the ith have the following characteristics:
when the drainage hole is an incomplete well, if the distance D between holes is less than 1.5M, when the mutual interference of the incomplete additional resistance coefficient zeta is not considered, the drainage hole for the ith has the following characteristics:
the water level descending of all the control points is to achieve the purpose of safe mining, for any jth water level control point, the safe water level descending is Ssi, and the constraint conditions are as follows:
when the drainage hole is a complete well or an incomplete well but r is more than 1.5M, the following steps are carried out:
when the water discharge hole is an incomplete well and the linear distance r between the water level control point and the water discharge hole is less than 1.5M, the following steps are carried out:
the minimum optimization model of Q and n is as follows: min (Q) and min (n)
Preferably, in S5, a multi-objective optimization method is used to calculate an optimal Q, n value of the model, and for the problem of dredging and lowering a water-bearing stratum of limestone with high bearing pressure and low permeability, because the water-bearing stratum has a strong heterogeneity characteristic, safe mining in the whole working plane is generally realized in a form of encrypted hydrophobic drilling, but if the water-bearing stratum is too much drilled, the mining cost is inevitably increased, if the drilled holes are too many, the water discharge amount of a single hole is inevitably increased, but because of the low permeability characteristic, the water discharge amount of the single hole is not too large, otherwise, the dredging and lowering effect is poor, and the hidden danger of water inrush possibly caused by insufficient pressure lowering is caused.
The method adopts a multi-objective particle swarm optimization (MOPSO) algorithm for solving, the MOPSO algorithm is expanded into a multi-objective optimization design, and the calculation flow is shown in figure 2.
Preferably, the optimization model is discussed in S5, production practice is guided, and an optimal scheme (meeting the conditions of reducing single-hole flow, properly encrypting drainage holes, and the like) is formulated according to the optimization result to solve the problem of dredging and lowering the non-homogeneous and low-permeability limestone aquifer in the actual engineering.
The invention has the beneficial effects that:
1. aiming at the coal seam floor high-pressure-bearing water dredging and lowering project, starting from the low-permeability characteristic of a water-bearing stratum and the requirements of economy, environmental protection and the like, under the condition of meeting the safe production condition, the minimum drainage hole is constructed, the minimum single-hole flow reaches the optimal dredging and lowering effect, the total drainage amount and the drilling project amount are reduced to the maximum extent, the optimal balance among safety, economy and environmental benefits is realized, the construction cost is reduced, and the purpose of safe exploitation is also achieved;
2. the method is characterized in that the positions of the underground scattered drilling holes are selected in a roadway, the roadway where the drilling hole candidate positions are located is generalized into a straight line segment, the optimal point position meeting the conditions is searched on the straight line segment, and accurate hole position coordinates can be obtained. Because the aquifer to be dredged and degraded needs to simultaneously meet the constraint conditions that the flow of a single hole is minimum to achieve the optimal dredging and degrading effect and the number of holes is minimum, the method adopts the multi-target particle swarm algorithm, has the advantages of high searching speed, high efficiency and simple algorithm, and can quickly obtain the optimal solution of the model.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
In the embodiment, a certain mine is taken as an example, 10 coal seams are mainly mined, the distance between the 10 coal seams and the top of a taiyuan limestone confined aquifer is about 35-55m, 1-4 limestone aquifers at the upper section of the 10 coal seam mainly influencing the safe mining of the 10 coal are used, the mine strata run towards the south and the north, incline to the east, the inclination angle is about 10 degrees, the mine strata are of a monoclinic structure, the fault in a well field develops, but most of the water is not guided and does not contain water, and due to the water blocking effect of the fault, the taiyuan limestone group at the upper section of the taiyuan group is a closed weak aquifer, so that the feasibility is.
In the exploration and production stages, 22-23-1, 27-1, 07-3 and other drilling holes of a mine are sequentially used for carrying out a water pumping test on the ash aquifer at the upper section 1-4 of the Tai-gray, and the water-rich parameter results are shown in a table 1.
TABLE 1 mine Taogen group 1-4 limestone pumping test
The pumping results in table 1 show that the upper 1-4 ash aquifer of the taiyuan group has q equal to 0.0052-0.773L/sm, K equal to 0.026-3.93m/d, the shallow hidden outcrop has strong water-rich property and communication property, and gradually weakens to a weak permeable layer towards the deep.
The 1013 and 1011 working faces are located in the 101 mining area, the mining area is divided into several blocks by several large faults of F9, F10 and F11, wherein the 1013 and 1011 working faces are located between two broken layers of F9 and F10, the hydrogeological conditions are relatively independent, the water supply source is limited due to ash, and the dredging property is strong. The basement of the working face is a Taiyuan high pressure bearing limestone aquifer, wherein the thickness of the 1-4 limestone aquifers is 50m, 126 1-4 limestone water detection drill holes constructed in a mining area are drilled, and the water outlet is more than 5m3Perh is only 4 holes, most of which are dry holes, the water yield is less than 2%, and the characteristics of good water-dredging property and insufficient supply amount of bottom plate limestone are further proved。
Before stoping, the water level elevation of the Tai-grey water bit long observation hole 14-1 in the mining area is-50 m, the minimum thickness of the water-resisting layer of the coal bed and the aquifer is 45m, the maximum water pressure borne by the bottom of the water-resisting layer is 4.65MPa according to the calculation of the minimum elevation of the working surface which is-515 m, the water inrush coefficient is about 0.103, and the maximum water pressure is larger than the critical value in the water prevention and control requirements of the coal mine, so that the water inrush danger is caused.
In order to ensure the smooth extraction of two working faces, a special water discharge roadway is constructed in a mine, as shown in figure 3, the elevation of a bottom plate of the water discharge roadway is-541.9 to-540.7 m, the total length is 390m, according to a hydrological exploration report, the thickness of a water-containing layer of 1-4 ash in a mining area is 50m, the permeability coefficient K is 0.03m/d, the water storage coefficient is 1 × 10-6, the radius rw of a water discharge hole is 0.05m, and the maximum depth Smax of a single hole is 490 m.
And S1, in the actual water discharge project, taking the range boundary of the covered mining working face as a dredging and descending boundary, and selecting the end point of the working face as a water level control point in the dredging and descending boundary, as shown in figure 1.
S2: 1-4 ash has the characteristics of typical high pressure bearing and low permeability. According to the principle of selecting water level control points, the ends of the upper roadway and the lower roadway far away from the drainage drill holes and the working surface with large burial depth are used as control points, as shown in fig. 3, each control point is used as a water level observation point, the water pressure borne by the bottom of the water-resisting layer of the coal seam floor is observed, namely, z is 0, and meanwhile, d of the drainage drill holes is 0. For convenience of calculation, the lower left corner of the attached drawing 3 is taken as a coordinate origin to obtain relative coordinates of each control point, the relative coordinates are shown in the following table, according to the water prevention and control requirements of a coal mine, the interfacial water inrush coefficient of a working face is 0.06MPa/m, and the safe water pressure of the obtained control point is 2.7 MPa; the safe water head elevation of the control points and the safe water level depression of each control point can be obtained according to the water head formula and are shown in table 2.
TABLE 2 safe water level drawdown at water level control point
S3: the water drainage drill holes are planned to be arranged in an AB section of the roadway, and the roadway is generalized to be a straight line equation according to coordinates of two points A (2008.42, 582.83) and B (2146.8, 408.48): y is-1.23 x +3057.4(2008.42 ≦ x ≦ 2146.8).
S4: the established optimization model is utilized, a multi-objective particle swarm optimization algorithm is applied, the model is solved based on Matlab software, the calculation result is shown in table 3, the attached figure 4 is the solution which is given by taking the complete well as an example and meets the constraint condition, and the program gives the final result according to the optimization objective of less holes and small single-hole flow and the principle of minimum total flow.
TABLE 3 optimized solution results
S5: and (4) solving the optimal hole number and single-hole flow value of the model by using a multi-objective optimization method.
S6: as can be seen from Table 3, as l/M is reduced (the degree of penetration of the drill hole into the aquifer), n is gradually increased, the increasing amplitude is gradually increased, the flow rate of a single hole is gradually reduced, the difference of the total flow rate is not large, and the total drilling project amount shows a trend of reducing in a stepped manner; aiming at the low-permeability heterogeneous aquifer, a better dredging effect can be achieved through the encrypted drilling and the single-hole small flow rate; however, when the single-hole flow is small, the problems of more holes, high construction cost and the like exist; meanwhile, the aquifer is a composite aquifer consisting of 1-4 ashes, the 3-4 ashes are generally good in water-rich property, and if l/M is small, the problem that the 3-4 ashes are difficult to enter a water drainage hole exists; according to the optimization calculation result, when l/M is 0.6, the hole number and the drilling engineering quantity are moderate, the total water discharge quantity is also minimum, and the water discharge method is an ideal water discharge design scheme.
The mine refers to an optimization result in the actual water discharge engineering design, adopts an incomplete well scheme, simultaneously reduces the single-hole flow, properly encrypts water discharge holes and other measures to solve the dredging and descending problem of the heterogeneous and low-permeability limestone aquifer, and totally arranges 6 layer-crossing water discharge holes in two drilling sites of a water discharge roadway, as shown in figure 5; the final holes of the other 5 water discharging holes are all in the four ashes, except that the final hole of the F1-2 hole is in the three ashes; the water discharge amount of a single hole is concentrated at 52m3/d-100m3In addition, in order to better achieve the dredging effect, 4 approximately horizontal bedding water discharging holes are constructed in a No. 3 drill site of the water discharging lane along azimuth angles of 170 degrees, 196 degrees, 206 degrees and 222 degrees respectively, the water discharging holes are basically parallel to the trend of the four ashes, and the average length is 140 m.
Observing the water level change of the 14-view 1-hole, as shown in figure 6, the total water discharge amount of the over-ash in the 101 mining area is 50m3Per hour, cumulative limestone water drainage of about 49 ten thousand meters3. The water level of a 14-view 1 (outside the outburst of a mining area) of the ground space-ash observation hole is reduced to-170 m from-50 m, the accumulated water level is reduced by 120m, and the maximum daily reduction amplitude is 2.4 m; the water pressure of the underground pressure measuring hole on the working face is 0.2MPa, the corresponding water level is-320 m, the lowest point elevation of the extraction section on the working face is-401.1 m, the maximum water inrush coefficient of the working face is 0.026, and the aim of dry extraction is basically achieved.
In the description herein, references to the description of "one embodiment," "an example," "a specific example" or the like are intended to mean that a particular feature, structure, material, or characteristic described in connection with the embodiment or example is included in at least one embodiment or example of the invention. In this specification, the schematic representations of the terms used above do not necessarily refer to the same embodiment or example. Furthermore, the particular features, structures, materials, or characteristics described may be combined in any suitable manner in any one or more embodiments or examples.
The foregoing shows and describes the general principles, essential features, and advantages of the invention. It will be understood by those skilled in the art that the present invention is not limited to the embodiments described above, which are described in the specification and illustrated only to illustrate the principle of the present invention, but that various changes and modifications may be made therein without departing from the spirit and scope of the present invention, which fall within the scope of the invention as claimed.