CN110532609B - Grouting pressure simulation method and device based on partition equivalent grouting pressure vector - Google Patents
Grouting pressure simulation method and device based on partition equivalent grouting pressure vector Download PDFInfo
- Publication number
- CN110532609B CN110532609B CN201910672411.8A CN201910672411A CN110532609B CN 110532609 B CN110532609 B CN 110532609B CN 201910672411 A CN201910672411 A CN 201910672411A CN 110532609 B CN110532609 B CN 110532609B
- Authority
- CN
- China
- Prior art keywords
- grouting
- area
- equivalent
- grouting pressure
- pressure
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Landscapes
- Consolidation Of Soil By Introduction Of Solidifying Substances Into Soil (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a grouting pressure simulation method and device based on a subarea equivalent grouting pressure vector. Wherein, the method comprises the following steps: dividing a plurality of grouting areas according to the elevation of the grouting section and the inclination angle of the grouting hole; setting the sensitivity coefficient of the grouting area according to the lithology of the grouting area and the grouting construction arrangement condition; establishing a subarea equivalent grouting pressure vector model according to the grouting area and the sensitivity coefficient of the grouting area, wherein the expression of the subarea equivalent grouting pressure vector is as follows:
the subarea equivalent grouting pressure vector model is used for simulating the complex distribution of grouting pressure in the fracture network. The invention solves the technical problem that the actual grouting pressure can not be truly simulated by a numerical method to guide the engineering construction.
Description
Technical Field
The invention relates to the field of grouting construction, in particular to a grouting pressure simulation method and device based on a subarea equivalent grouting pressure vector.
Background
The grouting construction process has strong concealment, complexity and uncertainty, and the construction quality and effect determine the success or failure of the project. Grouting pressure is a major source of grouting energy and is an important factor in controlling and improving grouting quality. The simulation of the actual grouting pressure of the construction site is carried out by a numerical simulation method, and the selection of the grouting pressure and the grouting arrangement of the construction site are guided. However, at present, no method for truly simulating actual grouting pressure exists, so that the guiding significance of the numerical simulation result on the engineering is greatly reduced.
The technical problems existing in the conventional techniques are mainly as follows: (1) the existing grouting pressure simulation method rarely considers the different elevations of the grouting sections of different grouting holes, and the same elevation is often approximately taken; (2) the traditional grouting pressure simulation method rarely considers the inclination angle of a grouting hole, and generally approximately considers that the grouting holes are all vertical to the grouting surface; (3) the conventional grouting pressure simulation method rarely considers the grouting construction sequence, namely the grouting sequence of I, II and III grouting holes and the grouting sequence of the same grouting hole; (4) the concrete geological conditions of the grouting area are difficult to consider by the traditional grouting pressure simulation method.
In view of the above problems, no effective solution has been proposed.
Disclosure of Invention
The embodiment of the invention provides a grouting pressure simulation method and device based on a subarea equivalent grouting pressure vector, which at least solve the technical problem that the actual grouting pressure cannot be truly simulated through a numerical method to guide engineering construction.
According to an aspect of an embodiment of the present invention, there is provided a grouting pressure simulation method based on a partitioned equivalent grouting pressure vector, including: marking according to the elevation of the grouting section and the inclination angle of the grouting pipeDividing a plurality of grouting areas; setting the sensitivity coefficient of the grouting area according to the lithology of the grouting area and the grouting construction arrangement condition; establishing a subarea equivalent grouting pressure vector model according to the grouting area and the sensitivity coefficient of the grouting area, wherein the expression of the subarea equivalent grouting pressure vector is as follows:
wherein N is the number of the plurality of grouting sections, and i is the ith grouting area in the plurality of grouting areas; sicThe susceptibility of the ith grouting area; si1，Si2And Si3Showing the areas of I, II and III grouting holes of the ith grouting area; n isi1，ni2And ni3The number of I, II and III grouting holes of the ith grouting area is shown; pi1，Pi2And Pi3The grouting pressure of I, II and III grouting holes of the ith grouting area is shown; r isijThe grouting influence radius of the jth grouting hole in the ith grouting area is set; z is a unit vector vertically upward; fiAnd the subarea equivalent grouting pressure vector model is used for simulating the complex distribution of the grouting pressure in the fracture network.
Further, after a zonal equivalent grouting pressure vector model is established according to the grouting area and the sensitivity coefficient of the grouting area, the method further comprises the following steps: acquiring boundary conditions and material parameters of the grouting area; according to the boundary condition, the material parameter and the partition equivalent grouting pressure vector FiCalculating grouting uplift displacement delta H by a numerical simulation methodz' (X), wherein X is the Xth monitoring point of the grout region; establishing an inversion optimization objective function G ═ Σ r2(X)＝∑(ΔHz(X)-ΔHz′(X))2Wherein, Δ Hz(X) actual grouting lifting displacement; adjusting the coefficient of sensitivity S of the grouted areaicAnd the zoning of the grout zone to optimize the objective function for the inversionThe value of G reaches a preset range; when the inversion optimization objective function G reaches the preset range, calculating a subarea equivalent grouting pressure vector F according to the adjusted sensitive coefficient of the grouting area and the grouting areai。
Further, adjusting the sensitivity coefficient S of the grouting areaicAnd the partition of the grouting area, so that the value of the inversion optimization objective function G reaches a preset range comprises the following steps: adjusting the sensitivity coefficient of the grouting area to enable the value of the inversion optimization target function G to reach a preset range; and if the value of the inversion optimization target function G cannot reach the preset range, adjusting the partition of the grouting area so as to enable the value of the inversion optimization target function G to reach the preset range.
Further, after a zonal equivalent grouting pressure vector model is established according to the grouting area and the sensitivity coefficient of the grouting area, the method further comprises the following steps: determining the value range of the sensitivity coefficient of the grouting area according to the lithology of the grouting area and the grouting construction arrangement condition; and calculating to obtain the value range of the partition equivalent grouting pressure vector according to the value range of the sensitive coefficient of the grouting area.
According to another aspect of the embodiments of the present invention, there is also provided a grouting pressure simulation apparatus based on a partitioned equivalent grouting pressure vector, including: the dividing module is used for dividing a plurality of grouting areas according to the elevation of the grouting section and the inclination angle of the grouting pipe; the setting module is used for setting the sensitivity coefficient of the grouting area according to the lithology of the grouting area and the grouting construction arrangement condition; the establishing module is used for establishing a subarea equivalent grouting pressure vector model according to the grouting area and the sensitivity coefficient of the grouting area, wherein the expression of the subarea equivalent grouting pressure vector is as follows:
wherein N is the number of the plurality of grouting sections, and i is the ith grouting area in the plurality of grouting areas; sicThe susceptibility of the ith grouting area; si1，Si2And Si3Is shown asThe areas of the I, II and III grouting holes of the I grouting areas; n isi1，ni2And ni3Representing the number of I, II and III grouting holes of the ith grouting area; pi1，Pi2And Pi3Indicating grouting pressure of I, II and III grouting holes of the ith grouting area; r isijThe grouting influence radius of the jth grouting hole in the ith grouting area is set; z is a unit vector vertically upward; fiAnd the subarea equivalent grouting pressure vector model is used for simulating the complex distribution of the grouting pressure in the fracture network.
Further, the apparatus further comprises: the acquisition module is used for acquiring boundary conditions and material parameters of the grouting area; a first calculation module for calculating the equivalent grouting pressure vector F of the subarea according to the boundary condition and the material parameteriCalculating grouting uplift displacement delta H by a numerical simulation methodz' (X), wherein X is the Xth monitoring point of the grout region; a function module for establishing an inversion optimization target function G ═ Σ r2(X)＝∑(ΔHz(X)-ΔHz′(X))2Wherein, Δ Hz(X) actual grouting lifting displacement; an adjusting module for adjusting the sensitivity coefficient s of the grouting areaicAnd the partition of the grouting area, so that the value of the inversion optimization target function G reaches a preset range; a second calculation module, configured to calculate a zonal equivalent grouting pressure vector F according to the adjusted sensitivity coefficient of the grouting area and the grouting area when the inversion optimization objective function G reaches the preset rangei。
Further, the adjustment module includes: a first adjusting unit for adjusting the sensitivity coefficient S of the grouting areaicSo that the value of the inversion optimization target function G reaches a preset range; and the second adjusting unit is used for adjusting the subareas of the grouting area to enable the value of the inversion optimization target function G to reach the preset range when the value of the inversion optimization target function G cannot reach the preset range.
Further, the apparatus further comprises: the value taking module is used for determining the value taking range of the sensitivity coefficient of the grouting area according to the lithology of the grouting area and the grouting construction arrangement condition; and the third calculation module is used for calculating the value range of the partition equivalent grouting pressure vector according to the value range of the sensitive coefficient of the grouting area.
In the embodiment of the invention, a plurality of grouting areas are divided according to the elevation of a grouting section and the inclination angle of a grouting pipe; setting the sensitivity coefficient of the grouting area according to the lithology of the grouting area and the grouting construction arrangement condition; establishing a subarea equivalent grouting pressure vector model according to the grouting area and the sensitivity coefficient of the grouting area, wherein the expression of the subarea equivalent grouting pressure vector is as follows:
wherein N is the number of the plurality of grouting sections, and i is the ith grouting area in the plurality of grouting areas; sicThe susceptibility of the ith grouting area; si1，Si2And Si3Indicating the areas of I, II and III grouting holes of the ith grouting area; n isi1，ni2And ni3Representing the number of I, II and III grouting holes of the ith grouting area; pi1，Pi2And Pi3Indicating grouting pressure of I, II and III grouting holes of the ith grouting area; r isijThe grouting influence radius of the jth grouting hole in the ith grouting area is set; z is a unit vector vertically upward; fiThe method is characterized in that a partition equivalent grouting pressure vector model is used for simulating the complex distribution of grouting pressure in a fracture network in a mode of generating an equivalent grouting pressure vector which has a lifting effect on bedrock or a dam body above a grouting area in the grouting process, so that the aim of simulating the actual grouting pressure on a construction site by using a numerical simulation method is fulfilled, the technical effect of truly simulating the actual grouting pressure is achieved, and the technical problem is solved.
Drawings
The accompanying drawings, which are included to provide a further understanding of the invention and are incorporated in and constitute a part of this application, illustrate embodiment(s) of the invention and together with the description serve to explain the invention without limiting the invention. In the drawings:
fig. 1 is a flowchart of a grouting pressure simulation method based on a zonal equivalent grouting pressure vector according to an embodiment of the present invention;
FIG. 2 is a schematic diagram of a single zone equivalent grouting pressure assumption according to an embodiment of the present invention;
FIG. 3 is a schematic illustration of zonal equivalent grouting pressures according to an embodiment of the present invention;
fig. 4 is a flow chart of a partitioned equivalent grouting pressure vector application according to an embodiment of the present invention.
Detailed Description
In order to make the technical solutions of the present invention better understood, 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.
It should be noted that the terms "first," "second," and the like in the description and claims of the present invention and in the drawings described above are used for distinguishing between similar elements and not necessarily for describing a particular sequential or chronological order. It is to be understood that the data so used is interchangeable under appropriate circumstances such that the embodiments of the invention described herein are capable of operation in sequences other than those illustrated or described herein. Furthermore, the terms "comprises," "comprising," and "having," and any variations thereof, are intended to cover a non-exclusive inclusion, such that a process, method, system, article, or apparatus that comprises a list of steps or elements is not necessarily limited to those steps or elements expressly listed, but may include other steps or elements not expressly listed or inherent to such process, method, article, or apparatus.
In accordance with an embodiment of the present invention, there is provided a method embodiment for grouting pressure simulation based on partitioned equivalent grouting pressure vectors, it is noted that the steps illustrated in the flowchart of the drawings may be performed in a computer system such as a set of computer executable instructions, and that while a logical order is illustrated in the flowchart, in some cases the steps illustrated or described may be performed in an order different than here.
Fig. 1 is a grouting pressure simulation method based on a zonal equivalent grouting pressure vector according to an embodiment of the present invention, as shown in fig. 1, the method includes the following steps:
step S102, dividing a plurality of grouting areas according to the elevation of a grouting section and the inclination angle of a grouting pipe;
step S104, setting a sensitivity coefficient of a grouting area according to the lithology of the grouting area and grouting construction arrangement conditions;
step S106, establishing a subarea equivalent grouting pressure vector model according to the grouting area and the sensitivity coefficient of the grouting area, wherein the expression of the subarea equivalent grouting pressure vector is as follows:
wherein N is the number of the grouting subareas, and i is the ith grouting area in the grouting areas; sicThe susceptibility of the ith grouting area; si1，Si2And Si3Indicating the areas of I, II and III grouting holes of the ith grouting area; n isi1，ni2And ni3Representing the number of I, II and III grouting holes of the ith grouting area; pi1，Pi2And Pi3Indicating grouting pressure of I, II and III grouting holes of the ith grouting area; r isijThe grouting influence radius of the jth grouting hole in the ith grouting area is set; z is a unit vector vertically upward; fiAnd the subarea equivalent grouting pressure vector model is used for simulating the complex distribution of the grouting pressure in the fracture network.
First, the definition and representation of the zonal equivalent grouting pressure according to the present invention will be described. See fig. 2, 3.
Alternatively, the establishment of the zonal equivalent grouting pressure vector is based on several assumptions:
assume that 1: neglecting the weight of the grouted body, only considering the grouting pressure, in fig. 2, the lifting and pressing acting forces of the surfaces A and B of the inclined grouting hole are equal, and the surfaces A and B are taken as shown in the formula (1) and fig. 2 (the upper and lower surfaces are directed to the vertical direction, namely the direction perpendicular to the building base surface);
wherein L is1Is the length of any arch line in the A surface, E is the length of a semicircle, a and b are any pair of opposite (180 degrees) arch lines, and beta is the included angle between a grouting hole and the vertical direction.
Assume 2: the contact between the grouting pipe and the bedrock is good, and relative sliding and slurry mixing cannot occur;
assume that 3: large geological defects similar to karst caves do not exist in the grouting section;
assume 4: the lifting may be caused by slurry in any crack in the grouting section, the region above the crack is a lifting region, and the vertical direction pressure of the upper surface and the lower surface of any crack in the lifting region can be self-balanced, that is, formula (2) and fig. 2 show that:
wherein x is the total number of cracks in the lifting area, ciThe number of the crack is the number of the crack,
and
is a crack ciThe upper and lower vertical directions resultant force.
And (3) defining and representing the equivalent grouting pressure of the subareas:
assuming an equivalent grouting pressure acting surface at the very bottom of the grouting section, equation (3) can be derived from the equilibrium conditions for the grouting zone. Equations (4) and (3) are equivalent and are an expression of the partition equivalent grouting pressure vector. Equation (5) is a scalar expression of the zonal equivalent grouting pressure,
referred to as zonal equivalent grouting pressure vector FiAverage value of (a).
Fi＝[f1，f2，f3，...，fN] (4)
Optionally, after the partitioned equivalent grouting pressure vector model is established according to the grouting area and the sensitivity coefficient of the grouting area, the method further includes: acquiring boundary conditions and material parameters of a grouting area; according to boundary conditions, material parameters and partition equivalent grouting pressure vector FiCalculating grouting uplift displacement delta H by a numerical simulation methodz' (X), wherein X is the Xth monitoring point of the grout region; establishing an inversion optimization objective function
Wherein, Δ Hz(X) actual grouting lifting displacement; adjusting the susceptibility s of the grouted areaicPartitioning the grouting area to enable the value of the inversion optimization target function G to reach a preset range; when the inversion optimization objective function G reaches a preset range, calculating a subarea equivalent grouting pressure vector F according to the adjusted sensitive coefficient of the grouting area and the grouting areai. Optionally, the boundary condition of the grouting area is normal constraint, three-way constraint and the like of the grouting area;material parameters the material parameters of the body (e.g., rock, etc.) being irrigated; the grout region includes a plurality of monitoring points.
Optionally, adjusting the susceptibility S of the grouted areaicAnd partitioning the grouting area so that the value of the inversion optimization objective function G reaches a preset range comprises the following steps: adjusting the sensitivity coefficient of the grouting area to enable the value of the inversion optimization target function G to reach a preset range; and if the value of the inversion optimization target function G cannot reach the preset range, adjusting the partition of the grouting area so as to enable the value of the inversion optimization target function G to reach the preset range.
Optionally, after the partitioned equivalent grouting pressure vector model is established according to the grouting area and the sensitivity coefficient of the grouting area, the method further includes: determining the value range of the sensitivity coefficient of the grouting area according to the lithology of the grouting area and the grouting construction arrangement condition; and calculating to obtain the value range of the partition equivalent grouting pressure vector according to the value range of the sensitive coefficient of the grouting area.
As shown in fig. 4, in an alternative embodiment,
1) and performing rough grouting partition according to the construction geological conditions, the height of the grouting section and the construction sequence.
2) Selecting an equivalent grouting pressure vector FiThe range is as follows: a is less than or equal to FiB (i ═ 1, 2, 3, 4.., n), where a and B are the upper and lower limits of equivalent grouting pressure, S is usually estimated from geological exploration, sonic testing, actual grouting pressure measurement of petrophysical properties, and other datacTo obtain FiThe upper and lower limits of (2).
3) Combining boundary conditions, material parameters and FiIs brought into the model and calculated by numerical simulation methods (e.g. finite element method), FiApplied to the plane of the bottom of the section of grout hole (see figures 2 and 3).
4) Component Δ H of lift displacement to be monitored in situz(X) and the result of numerical simulation calculation Δ Hz'X' is used for making difference, the difference value of the lifting displacement obtained by actual measurement and numerical simulation calculation is shown as a formula (4), and X represents the X-th monitoring point participating in making difference.
r(X)＝ΔHz(X)-ΔHz(X) (4)
And establishing an inversion optimization objective function according to the following steps:
G＝∑r2(X)＝∑(ΔHz(X)-ΔHz′(X))2 (5)
by adjusting SicAnd continuously reducing the optimized objective function value G by the numerical model to finally obtain the partition equivalent grouting pressure vector suitable for concrete arch dam foundation grouting, and optimizing the optimized SicThe model is brought in to simulate the influence of grouting opportunity, partition grouting and grouting pressure on the lifting, cracking and other safety problems, and the lifting value of the arch dam under different grouting construction sequences, pressure combinations and other conditions can be predicted.
Optionally, if by adjusting SicAnd if the optimization function G can not meet the accuracy requirement, adjusting the grouting subareas.
According to another aspect of the embodiments of the present invention, there is also provided a grouting pressure simulation apparatus based on a partitioned equivalent grouting pressure vector, including: the dividing module is used for dividing a plurality of grouting areas according to the elevation of the grouting section and the inclination angle of the grouting pipe; the setting module is used for setting the sensitivity coefficient of the grouting area according to the lithology of the grouting area and the grouting construction arrangement condition; the establishing module is used for establishing a subarea equivalent grouting pressure vector model according to the grouting area and the sensitivity coefficient of the grouting area, wherein the expression of the subarea equivalent grouting pressure vector is as follows:
wherein N is the number of the grouting subareas, and i is the ith grouting area in the grouting areas; sicThe susceptibility of the ith grouting area; si1，Si2And Si3Indicating the areas of I, II and III grouting holes of the ith grouting area; n isi1，ni2And ni3Representing the number of I, II and III grouting holes of the ith grouting area; pi1，Pi2And Pi3Indicating grouting pressure of I, II and III grouting holes of the ith grouting area; r isijGrouting for j-th grouting hole in ith grouting areaThe radius of influence; z is a unit vector vertically upward; fiAnd the subarea equivalent grouting pressure vector model is used for simulating the complex distribution of the grouting pressure in the fracture network.
Optionally, the apparatus further comprises: the acquisition module is used for acquiring boundary conditions and material parameters of a grouting area; a first calculation module for calculating a zonal equivalent grouting pressure vector F according to the boundary conditions and the material parametersiCalculating grouting uplift displacement delta H by a numerical simulation methodz' (X), wherein X is the Xth monitoring point of the grout region; a function module for establishing an inversion optimization target function G ═ Σ r2(X)＝∑(ΔHz(X)-ΔHz′(X))2Wherein, Δ Hz(X) actual grouting lifting displacement; an adjusting module for adjusting the sensitivity coefficient s of the grouting areaicPartitioning the grouting area to enable the value of the inversion optimization target function G to reach a preset range; a second calculation module, configured to calculate a zonal equivalent grouting pressure vector F according to the adjusted sensitivity coefficient of the grouting area and the grouting area when the inversion optimization objective function G reaches a preset rangei。
Optionally, the adjusting module comprises: a first adjusting unit for adjusting the sensitivity coefficient s of the grouting areaicSo that the value of the inversion optimization target function G reaches a preset range; and the second adjusting unit is used for adjusting the partition of the grouting area to enable the value of the inversion optimization target function G to reach the preset range when the value of the inversion optimization target function G cannot reach the preset range.
Optionally, the apparatus further comprises: the value taking module is used for determining the value taking range of the sensitivity coefficient of the grouting area according to the lithology of the grouting area and the grouting construction arrangement condition; and the third calculation module is used for calculating the value range of the partition equivalent grouting pressure vector according to the value range of the sensitive coefficient of the grouting area.
The above-mentioned serial numbers of the embodiments of the present invention are merely for description and do not represent the merits of the embodiments.
In the above embodiments of the present invention, the descriptions of the respective embodiments have respective emphasis, and for parts that are not described in detail in a certain embodiment, reference may be made to related descriptions of other embodiments.
In the embodiments provided in the present application, it should be understood that the disclosed technology can be implemented in other ways. The above-described embodiments of the apparatus are merely illustrative, and for example, the division of the units may be a logical division, and in actual implementation, there may be another division, for example, multiple units or components may be combined or integrated into another system, or some features may be omitted, or not executed. In addition, the shown or discussed mutual coupling or direct coupling or communication connection may be an indirect coupling or communication connection through some interfaces, units or modules, and may be in an electrical or other form.
The units described as separate parts may or may not be physically separate, and parts displayed as units may or may not be physical units, may be located in one place, or may be distributed on a plurality of units. Some or all of the units can be selected according to actual needs to achieve the purpose of the solution of the embodiment.
In addition, functional units in the embodiments of the present invention may be integrated into one processing unit, or each unit may exist alone physically, or two or more units are integrated into one unit. The integrated unit can be realized in a form of hardware, and can also be realized in a form of a software functional unit.
The integrated unit, if implemented in the form of a software functional unit and sold or used as a stand-alone product, may be stored in a computer readable storage medium. Based on such understanding, the technical solution of the present invention may be embodied in the form of a software product, which is stored in a storage medium and includes instructions for causing a computer device (which may be a personal computer, a server, or a network device) to execute all or part of the steps of the method according to the embodiments of the present invention. And the aforementioned storage medium includes: a U-disk, a Read-Only Memory (ROM), a Random Access Memory (RAM), a removable hard disk, a magnetic or optical disk, and other various media capable of storing program codes.
The foregoing is only a preferred embodiment of the present invention, and it should be noted that, for those skilled in the art, various modifications and decorations can be made without departing from the principle of the present invention, and these modifications and decorations should also be regarded as the protection scope of the present invention.
Claims (4)
1. A grouting pressure simulation method based on a subarea equivalent grouting pressure vector is characterized by comprising the following steps:
dividing a plurality of grouting areas according to the elevation of the grouting section and the inclination angle of the grouting pipe;
setting the sensitivity coefficient of the grouting area according to the lithology of the grouting area and the grouting construction arrangement condition;
establishing a subarea equivalent grouting pressure vector model according to the grouting area and the sensitivity coefficient of the grouting area, wherein the expression of the subarea equivalent grouting pressure vector is as follows:
wherein N is the number of the plurality of grouting sections, and i is the ith grouting area in the plurality of grouting areas; sicThe susceptibility of the ith grouting area; si1，Si2And Si3Indicating the areas of I, II and III grouting holes of the ith grouting area; n isi1，ni2And ni3Representing the number of I, II and III grouting holes of the ith grouting area; pi1，Pi2And Pi3Indicating grouting pressure of I, II and III grouting holes of the ith grouting area; r isijThe grouting influence radius of the jth grouting hole in the ith grouting area is set; z is a unit vector vertically upward; fiFor grouting processThe generated equivalent grouting pressure vector which has a lifting effect on bedrock or a dam body above a grouting area is used for simulating the complex distribution of grouting pressure in a fracture network;
after a partitioned equivalent grouting pressure vector model is established according to the grouting area and the sensitivity coefficient of the grouting area, the method further comprises the following steps:
acquiring boundary conditions and material parameters of the grouting area;
according to the boundary condition, the material parameter and the partition equivalent grouting pressure vector FiCalculating grouting uplift displacement delta H by a numerical simulation methodz' (X), wherein X is the Xth monitoring point of the grout region;
establishing an inversion optimization objective function G ═ Σ r2(X)＝∑(ΔHz(X)-ΔHz′(X))2Wherein, Δ Hz(X) actual grouting lifting displacement;
adjusting the coefficient of sensitivity S of the grouted areaicAnd the partition of the grouting area, so that the value of the inversion optimization target function G reaches a preset range;
when the inversion optimization target function G reaches the preset range, calculating a subarea equivalent grouting pressure vector F according to the adjusted sensitive coefficient of the grouting area and the grouting areai；
Adjusting the coefficient of sensitivity S of the grouted areaicAnd the partition of the grouting area, so that the value of the inversion optimization objective function G reaches a preset range comprises the following steps:
adjusting the sensitivity coefficient of the grouting area to enable the value of the inversion optimization target function G to reach a preset range;
and if the value of the inversion optimization target function G cannot reach the preset range, adjusting the partition of the grouting area so as to enable the value of the inversion optimization target function G to reach the preset range.
2. The method for simulating grouting pressure based on partitioned equivalent grouting pressure vectors according to claim 1, wherein after establishing a partitioned equivalent grouting pressure vector model according to the grouting area and the sensitivity coefficient of the grouting area, the method further comprises:
determining the value range of the sensitivity coefficient of the grouting area according to the lithology of the grouting area and the grouting construction arrangement condition;
and calculating to obtain the value range of the partition equivalent grouting pressure vector according to the value range of the sensitive coefficient of the grouting area.
3. The utility model provides a grouting pressure analogue means based on subregion equivalence grouting pressure vector which characterized in that includes:
the dividing module is used for dividing a plurality of grouting areas according to the elevation of the grouting section and the inclination angle of the grouting pipe;
the setting module is used for setting the sensitivity coefficient of the grouting area according to the lithology of the grouting area and the grouting construction arrangement condition;
the establishing module is used for establishing a subarea equivalent grouting pressure vector model according to the grouting area and the sensitivity coefficient of the grouting area, wherein the expression of the subarea equivalent grouting pressure vector is as follows:
wherein N is the number of the plurality of grouting sections, and i is the ith grouting area in the plurality of grouting areas; sicThe susceptibility of the ith grouting area; si1，Si2And Si3Indicating the areas of I, II and III grouting holes of the ith grouting area; n isi1，ni2And ni3Representing the number of I, II and III grouting holes of the ith grouting area; pi1，Pi2And Pi3Indicating grouting pressure of I, II and III grouting holes of the ith grouting area; r isijThe grouting influence radius of the jth grouting hole in the ith grouting area is set; z is a unit vector vertically upward; fiEquivalent grouting pressure vectors generated for the grouting process and having a lifting effect on bedrock or dam bodies above the grouting area, the subareasThe equivalent grouting pressure vector model is used for simulating the complex distribution of grouting pressure in the fracture network;
the device further comprises:
the acquisition module is used for acquiring boundary conditions and material parameters of the grouting area;
a first calculation module for calculating the equivalent grouting pressure vector F according to the boundary condition, the material parameter and the subareaiCalculating grouting uplift displacement delta H by a numerical simulation methodz' (X), wherein X is the Xth monitoring point of the grout region;
a function module for establishing an inversion optimization target function G ═ Σ r2(X)＝∑(ΔHz(X)-ΔHz′(X))2Wherein, Δ Hz(X) actual grouting lifting displacement;
an adjusting module for adjusting the sensitivity coefficient S of the grouting areaicAnd the partition of the grouting area, so that the value of the inversion optimization target function G reaches a preset range;
a second calculation module, configured to calculate a zonal equivalent grouting pressure vector F according to the adjusted sensitivity coefficient of the grouting area and the grouting area when the inversion optimization objective function G reaches the preset rangei；
The adjustment module includes:
a first adjusting unit for adjusting the sensitivity coefficient S of the grouting areaicSo that the value of the inversion optimization target function G reaches a preset range;
and the second adjusting unit is used for adjusting the subareas of the grouting area to enable the value of the inversion optimization target function G to reach the preset range when the value of the inversion optimization target function G cannot reach the preset range.
4. A grouting pressure simulation device based on a zoned equivalent grouting pressure vector according to claim 3, wherein the device further comprises:
the value taking module is used for determining the value taking range of the sensitivity coefficient of the grouting area according to the lithology of the grouting area and the grouting construction arrangement condition;
and the third calculation module is used for calculating the value range of the partition equivalent grouting pressure vector according to the value range of the sensitive coefficient of the grouting area.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910672411.8A CN110532609B (en) | 2019-07-24 | 2019-07-24 | Grouting pressure simulation method and device based on partition equivalent grouting pressure vector |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910672411.8A CN110532609B (en) | 2019-07-24 | 2019-07-24 | Grouting pressure simulation method and device based on partition equivalent grouting pressure vector |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110532609A CN110532609A (en) | 2019-12-03 |
| CN110532609B true CN110532609B (en) | 2021-02-19 |
Family
ID=68661303
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910672411.8A Active CN110532609B (en) | 2019-07-24 | 2019-07-24 | Grouting pressure simulation method and device based on partition equivalent grouting pressure vector |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110532609B (en) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102629107A (en) * | 2012-04-18 | 2012-08-08 | 长沙理工大学 | Pressure closed loop control method for large-circulation grouting monitor system |
| CN103018152A (en) * | 2012-12-20 | 2013-04-03 | 湖南科技大学 | Indoor grouting testing device under simulation of complex stress effect and testing method thereof |
| CN106284211A (en) * | 2016-07-15 | 2017-01-04 | 葛洲坝集团试验检测有限公司 | A kind of method improving distress in concrete inclined hole forensic chemistry grouting quality |
| CN107506557A (en) * | 2017-09-01 | 2017-12-22 | 宁波市交通规划设计研究院有限公司 | A kind of non-linear inversion analysis method of tunnel weak surrounding rock mechanics parameter |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102704497B (en) * | 2012-05-31 | 2015-05-27 | 中冶集团武汉勘察研究院有限公司 | RPW (ranked position weight) process for treating high fill slope |
| CN108316347A (en) * | 2018-03-09 | 2018-07-24 | 中冶京诚工程技术有限公司 | The assembled underground pipe gallery and its waterproof method of partitioning waterproofing |
-
2019
- 2019-07-24 CN CN201910672411.8A patent/CN110532609B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102629107A (en) * | 2012-04-18 | 2012-08-08 | 长沙理工大学 | Pressure closed loop control method for large-circulation grouting monitor system |
| CN103018152A (en) * | 2012-12-20 | 2013-04-03 | 湖南科技大学 | Indoor grouting testing device under simulation of complex stress effect and testing method thereof |
| CN106284211A (en) * | 2016-07-15 | 2017-01-04 | 葛洲坝集团试验检测有限公司 | A kind of method improving distress in concrete inclined hole forensic chemistry grouting quality |
| CN107506557A (en) * | 2017-09-01 | 2017-12-22 | 宁波市交通规划设计研究院有限公司 | A kind of non-linear inversion analysis method of tunnel weak surrounding rock mechanics parameter |
Non-Patent Citations (3)
| Title |
|---|
| Study on Optimal Grouting Timing for Controlling Uplift Deformation of a Super High Arch Dam;Peng Lin等;《Rock Mech Rock Eng》;20150318;第115-142页 * |
| 特高拱坝开裂分析方法与应用;朱晓旭;《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》;20130115(第1期);第C037-56页 * |
| 覆盖层地基注浆扩散数值模拟方法及应用;周正军等;《地下空间与工程学报》;20171130;第12卷;第646-650页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110532609A (en) | 2019-12-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US20180293789A1 (en) | Three-dimensional geomechanical modeling of casing deformation for hydraulic fracturing treatment design | |
| CN112012712B (en) | Numerical simulation method and device for water injection growth seam of embedded discrete seam | |
| US9157319B2 (en) | Stratigraphic modeling using production data density profiles | |
| Xu et al. | Response of framed buildings on raft foundations to tunneling | |
| CN103605900A (en) | Cross-scale complicated geologic body ground stress field identification method and device | |
| AU2013399193B2 (en) | Static earth model grid cell scaling and property re-sampling methods and systems | |
| CN115077473A (en) | Monitoring and analyzing method and system for whole forced landing and inclination correction construction process of building | |
| CN110532609B (en) | Grouting pressure simulation method and device based on partition equivalent grouting pressure vector | |
| CN105447319B (en) | Arch Dam Abutment sliding block difference geology level meters calculate determination method for parameter | |
| CN110705168A (en) | Simulation method of structural stress field | |
| CN108446251B (en) | Method for calculating active rock pressure of rock slope | |
| Sabzi et al. | The performance of buildings adjacent to excavation supported by inclined struts | |
| CN107451383B (en) | Calibration method for initial bed sand gradation of plane two-dimensional water sand mathematical model | |
| RU2745152C1 (en) | Method for combining a model of geological well passing with operational petrophysical interpretation of gis data in real time and a system implementing the method | |
| CN110320562B (en) | Method for correcting acoustic time difference in horizontal well of shale gas reservoir | |
| CN106777556A (en) | A kind of spacial analytical method for assessing slope excavating phase stable state | |
| Vahdati et al. | Inverse hardening soil parameter identification of an earth and rockfill dam by genetic algorithm optimization | |
| RU2631445C2 (en) | Method for determining number of workings while conducting engineering and geological surveys | |
| CN112861207B (en) | Method, apparatus and computer storage medium for subsidence prediction of a composite formation | |
| CN111980662B (en) | Rapid processing method for lateral logging data of inclined shaft anisotropic stratum array | |
| CN112558179A (en) | Intelligent inversion method of far field stress state | |
| EP3143597A1 (en) | Geomechanical modeling using dynamic boundary conditions from time-lapse data | |
| CN115510651B (en) | Ground surface transmission damage research simulation method for bedrock and aeolian sand in mining area | |
| CN111722276A (en) | Seismic guiding method and system for rock drilling | |
| Calvello | Calibration of soil constitutive laws by in-verse analysis |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |