CN110532609A - Grouting pressure analogy method and device based on the equivalent grouting pressure vector of subregion - Google Patents

Abstract

The invention discloses a method and device for simulating grouting pressure based on equivalent grouting pressure vectors of zones. Wherein, the method includes: 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 layout conditions; As well as the sensitivity coefficient of the grouting area, a partition equivalent grouting pressure vector model is established, wherein the expression of the partition equivalent grouting pressure vector is: The partition equivalent grouting pressure vector model is used to simulate the complex distribution of grouting pressure in the fracture network. The invention solves the technical problem that the actual grouting pressure cannot be truly simulated to guide engineering construction through a numerical method.

Description

Grouting pressure simulation method and device based on partition equivalent grouting pressure vector

technical field

The present invention relates to the field of grouting construction, in particular to a method and device for simulating grouting pressure based on equivalent grouting pressure vectors in zones.

Background technique

The grouting construction process has strong concealment, complexity and uncertainty, and its construction quality and effect determine the success or failure of the project. Grouting pressure is the main source of grouting energy and an important factor to control and improve grouting quality. It is of great significance to simulate the actual grouting pressure on the construction site by means of numerical simulation, and then guide the selection of grouting pressure and grouting layout on the construction site. However, there is currently no method that can simulate the actual grouting pressure more realistically, which leads to the greatly reduced guiding significance of the numerical simulation results for engineering.

The technical problems existing in the conventional technology are mainly as follows: (1) The previous grouting pressure simulation method seldom considered the different elevations of the grouting sections of different grouting holes, and often took approximately the same elevation; (2) The previous grouting pressure simulation method The inclination angle of the grouting hole is rarely considered, and it is generally considered that the grouting hole is perpendicular to the grouting surface; (3) The previous grouting pressure simulation method rarely considered the sequence of grouting construction, that is, grouting holes in sequence I, II and III and grouting in the same sequence (4) The previous grouting pressure simulation method is difficult to consider the specific geological conditions of the grouting area.

For the above problems, no effective solution has been proposed yet.

Contents of the invention

Embodiments of the present invention provide a method and device for simulating grouting pressure based on equivalent grouting pressure vectors in different regions, so as to at least solve the technical problem that the actual grouting pressure cannot be truly simulated to guide engineering construction through numerical methods.

According to an aspect of an embodiment of the present invention, a method for simulating grouting pressure based on the equivalent grouting pressure vector of a partition is provided, including: dividing multiple grouting areas according to the elevation of the grouting section and the inclination angle of the grouting pipe; The sensitivity coefficient of the grouting area is set according to the lithology and grouting construction layout conditions; according to the grouting area and the sensitivity coefficient of the grouting area, a partition equivalent grouting pressure vector model is established, wherein the expression of the partition equivalent grouting pressure vector for: Among them, N is the number of the multiple grouting areas, i is the i-th grouting area in the multiple grouting areas; S _ic is the sensitivity coefficient of the i-th grouting area; S _i1 , S _i2 and S _i3 represent The area of I, II and III sequence grouting holes in the i-th grouting area; n _i1 , _ni2 and n _i3 represent the number of I, II and III sequence grouting holes in the i-th grouting area; P _i1 , P _i2 and P _i3 represents the grouting pressure of the I, II and III sequence grouting holes in the i-th grouting area; r _ij is the grouting influence radius of the j-th sequence grouting hole in the i-th grouting area; z is the vertical upward unit vector; F _i is the equivalent grouting pressure vector generated during the grouting process that has a lifting effect on the bedrock or dam body above the grouting area, and the partition equivalent grouting pressure vector model is used to simulate the complex distribution of grouting pressure in the fracture network.

Further, after establishing the partition equivalent grouting pressure vector model according to the grouting area and the sensitivity coefficient of the grouting area, the method further includes: obtaining the boundary conditions and material parameters of the grouting area; , the material parameters and the equivalent grouting pressure vector F _i of the partition, the grouting displacement ΔH _z ′(X) is calculated by numerical simulation method, where X is the Xth monitoring point of the grouting area; Inversion optimization objective function G=∑r ² (X)=∑(ΔH _z (X)-ΔH _z ′(X)) ² , where ΔH _z (X) is the actual grouting displacement; adjust the grouting area Sensitivity coefficient S _ic and the division of the grouting area, so that the value of the inversion optimization objective function G reaches the preset range; when the inversion optimization objective function G reaches the preset range, according to the adjustment After the sensitivity coefficient of the grouting area and the grouting area calculate the partition equivalent grouting pressure vector F _i .

Further, adjusting the sensitivity coefficient S _ic of the grouting area and the division of the grouting area so that the value of the inversion optimization objective function G reaches a preset range includes: adjusting the sensitivity coefficient of the grouting area so that The value of the inversion optimization objective function G reaches a preset range; if the value of the inversion optimization objective function G cannot reach a preset range, then adjust the division of the grouting area so that the inversion optimization objective function G The value reaches the preset range.

Further, after establishing the zone equivalent grouting pressure vector model according to the grouting area and the sensitivity coefficient of the grouting area, the method further includes: determining the grouting area according to the lithology of the grouting area and the grouting construction layout conditions The value range of the sensitivity coefficient of the area; the value range of the equivalent grouting pressure vector of the partition is calculated according to the value range of the sensitivity coefficient of the grouting area.

According to another aspect of the embodiment of the present invention, there is also provided a grouting pressure simulation device based on the zone equivalent grouting pressure vector, including: a division module, which is used to divide multiple grouting areas according to the elevation of the grouting section and the inclination angle of the grouting pipe ; setting module, used to set the sensitivity coefficient of the grouting area according to the lithology of the grouting area and the grouting construction layout conditions; establishing module, used to establish partitions, etc. according to the grouting area and the sensitivity coefficient of the grouting area The effective grouting pressure vector model, where the expression of the partition equivalent grouting pressure vector is: Among them, N is the number of the multiple grouting areas, i is the i-th grouting area in the multiple grouting areas; S _ic is the sensitivity coefficient of the i-th grouting area; S _i1 , S _i2 and S _i3 represent The area of I, II and III order grouting holes in the i-th grouting area; n _i1 , n _i2 and n _i3 represent the number of I, II and III order grouting holes in the i-th grouting area; P _i1 , P _i2 and P _i3 represents the grouting pressure of the I, II and III sequence grouting holes in the i-th grouting area; r _ij is the grouting influence radius of the j-th sequence grouting hole in the i-th grouting area; z is the vertical unit vector; F _i is the equivalent grouting pressure vector generated during the grouting process that has a lifting effect on the bedrock or dam body above the grouting area, and the partition equivalent grouting pressure vector model is used to simulate the complex distribution of grouting pressure in the fracture network.

Further, the device further includes: an acquisition module, configured to acquire boundary conditions and material parameters of the grouting area; a first calculation module, configured to equivalently grout in the partition according to the boundary conditions and the material parameters The pressure vector F _i is calculated by the numerical simulation method to obtain the grouting displacement ΔH _z ′(X), where X is the Xth monitoring point of the grouting area; the function module is used to establish the inversion optimization objective function G= ∑r ² (X)=∑(ΔH _z (X)-ΔH _z ′(X)) ² , where ΔH _z (X) is the actual grouting displacement; the adjustment module is used to adjust the sensitivity of the grouting area The coefficient _sic and the division of the grouting area, so that the value of the inversion optimization objective function G reaches a preset range; the second calculation module is used for when the inversion optimization objective function G reaches the preset range , calculate the partition equivalent grouting pressure vector F _i according to the adjusted sensitivity coefficient of the grouting area and the grouting area.

Further, the adjustment module includes: a first adjustment unit, configured to adjust the sensitivity coefficient S _ic of the grouting area, so that the value of the inversion optimization objective function G reaches a preset range; a second adjustment unit, configured to If the value of the inversion optimization objective function G cannot reach a preset range, the division of the grouting area is adjusted so that the value of the inversion optimization objective function G reaches a preset range.

Further, the device also includes: a value module, used to determine the value range of the sensitivity coefficient of the grouting area according to the lithology of the grouting area and the grouting construction layout conditions; Calculate the value range of the sensitivity coefficient of the grouting area to obtain the value range of the equivalent grouting pressure vector of the partition.

In the embodiment of the present invention, a plurality of grouting areas are divided according to the elevation of the grouting section and the inclination angle of the grouting pipe; the sensitivity coefficient of the grouting area is set according to the lithology of the grouting area and the grouting construction layout conditions; according to the grouting area area and the sensitivity coefficient of the grouting area, and establish a partition equivalent grouting pressure vector model, where the expression of the partition equivalent grouting pressure vector is: Among them, N is the number of the multiple grouting areas, i is the i-th grouting area in the multiple grouting areas; S _ic is the sensitivity coefficient of the i-th grouting area; S _i1 , S _i2 and S _i3 represent The area of I, II and III order grouting holes in the i-th grouting area; n _i1 , n _i2 and n _i3 represent the number of I, II and III order grouting holes in the i-th grouting area; P _i1 , P _i2 and P _i3 represents the grouting pressure of the I, II and III sequence grouting holes in the i-th grouting area; r _ij is the grouting influence radius of the j-th sequence grouting hole in the i-th grouting area; z is the vertical unit vector; F _i The equivalent grouting pressure vector generated during the grouting process has a lifting effect on the bedrock above the grouting area or the dam body. The partition equivalent grouting pressure vector model is used to simulate the complex distribution of grouting pressure in the fracture network. The numerical simulation method is used to simulate the actual grouting pressure on the construction site, so as to realize the technical effect of simulating the actual grouting pressure and solve the technical problems.

Description of drawings

The accompanying drawings described here are used to provide a further understanding of the present invention and constitute a part of the application. The schematic embodiments of the present invention and their descriptions are used to explain the present invention and do not constitute improper limitations to the present invention. In the attached picture:

Fig. 1 is a flow chart of a method for simulating grouting pressure based on zone equivalent grouting pressure vectors according to an embodiment of the present invention;

Fig. 2 is a hypothetical schematic diagram of a single-zone equivalent grouting pressure according to an embodiment of the present invention;

Fig. 3 is a schematic diagram of the equivalent grouting pressure of a zone according to an embodiment of the present invention;

Fig. 4 is a flow chart of the application of zone equivalent grouting pressure vectors according to an embodiment of the present invention.

Detailed ways

In order to enable those skilled in the art to better understand the solutions of the present invention, the following will clearly and completely describe the technical solutions in the embodiments of the present invention in conjunction with the drawings in the embodiments of the present invention. Obviously, the described embodiments are only It is an embodiment of a part of the present invention, but not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by persons of ordinary skill in the art without making creative efforts shall fall within the protection scope of the present invention.

It should be noted that the terms "first" and "second" in the description and claims of the present invention and the above drawings are used to distinguish similar objects, but not necessarily used to describe a specific sequence or sequence. It is to be understood that the data so used are interchangeable under appropriate circumstances such that the embodiments of the invention described herein can be practiced in sequences other than those illustrated or described herein. Furthermore, the terms "comprising" and "having", as well as any variations thereof, are intended to cover a non-exclusive inclusion, for example, a process, method, system, product or device comprising a sequence of steps or elements is not necessarily limited to the expressly listed instead, may include other steps or elements not explicitly listed or inherent to the process, method, product or apparatus.

According to an embodiment of the present invention, a method embodiment of grouting pressure simulation based on the equivalent grouting pressure vector of a partition is provided. It should be noted that the steps shown in the flow chart of the accompanying drawings can be executed in a set of computer-executable instructions such as and, although a logical order is shown in the flowcharts, in some cases the steps shown or described may be performed in an order different from that shown or described herein.

Fig. 1 is a kind of grouting pressure simulation method based on partition 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 multiple grouting areas according to the elevation of the grouting section and the inclination angle of the grouting pipe;

Step S104, setting the sensitivity coefficient of the grouting area according to the lithology of the grouting area and the layout conditions of the grouting construction;

Step S106, according to the grouting area and the sensitivity coefficient of the grouting area, establish a partition equivalent grouting pressure vector model, where the expression of the partition equivalent grouting pressure vector is: Among them, N is the number of multiple grouting areas, i is the i-th grouting area in the multiple grouting areas; S _ic is the sensitivity coefficient of the i-th grouting area; S _i1 , S _i2 and S _i3 represent the i-th grouting area The area of grouting holes of sequence I, II and III in the area; n _i1 , n _i2 and n _i3 represent the number of grouting holes of sequence I, II and III in the i-th grouting region; P _i1 , P _i2 and P _i3 represent the number of grouting holes in the i-th grouting region The grouting pressure of the I, II and III sequence grouting holes in the i grouting area; r _ij is the grouting influence radius of the jth grouting hole in the i grouting area; z is the vertical upward unit vector; F _i is the grouting process The equivalent grouting pressure vector that has a lifting effect on the bedrock or dam body above the grouting area, the partition equivalent grouting pressure vector model is used to simulate the complex distribution of grouting pressure in the fracture network.

Firstly, the definition and expression of the partition equivalent grouting pressure in this invention will be described here. See Figure 2 and Figure 3.

Optionally, the establishment of the partition equivalent grouting pressure vector is based on the following assumptions:

Assumption 1: Neglect the weight of the grouted body and only consider the grouting pressure. In Figure 2, the upward and downward force of the A and B sides of the inclined grouting hole are equal, and the A and B sides are taken as formula (1) and figure 2 (up and down point to the vertical direction, that is, the direction perpendicular to the foundation plane);

Among them, L ₁ is the length of any arch line in surface A, E is the length of the semicircle, a, b are any pair of arch lines opposite (180°), and β is the angle between the grouting hole and the vertical direction.

Assumption 2: The contact between the grouting pipe and the bedrock is good, and there will be no relative sliding and grouting;

Hypothesis 3: There are no large geological defects similar to karst caves in the grouting section;

Hypothesis 4: The lifting may be caused by the grout in any crack in the grouting section, the area above this crack is the lifting area, and the vertical pressure of any crack in the lifting area can be self-balanced, that is Formula (2) and Figure 2 show:

where x is the total number of cracks in the lifted area, _ci is the number of cracks, and is the resultant force in the vertical direction above and below the crack _ci .

The definition and expression of the equivalent grouting pressure of the partition:

Assuming that the equivalent grouting pressure action surface is at the bottom of the grouting section, for the grouting area, equation (3) can be obtained through the equilibrium condition. Equation (4) is equivalent to (3), which is the expression of the equivalent grouting pressure vector of the partition. Equation (5) is the scalar expression of the equivalent grouting pressure of the partition, It is called the average value of the partition equivalent grouting pressure vector F _i .

F _i =[f ₁ , f ₂ , f ₃ , . . . , f _N ] (4)

Optionally, after establishing the zone equivalent grouting pressure vector model according to the grouting area and the sensitivity coefficient of the grouting area, the method further includes: obtaining the boundary conditions and material parameters of the grouting area; Vector F _i , the grouting lift displacement ΔH _z ′(X) is calculated by numerical simulation method, where X is the Xth monitoring point in the grouting area; the inversion optimization objective function is established Among them, ΔH _z (X) is the actual grouting displacement; adjust the sensitivity coefficient s _ic of the grouting area and the division of the grouting area so that the value of the inversion optimization objective function G reaches the preset range; when the inversion optimization objective function G When the preset range is reached, the partition equivalent grouting pressure vector F _i is calculated according to the adjusted sensitivity coefficient of the grouting area and the grouting area. Optionally, the boundary conditions of the grouting area are normal constraints, three-dimensional constraints, etc. of the grouting area; material parameters are material parameters of the body to be filled (such as rock); the grouting area includes multiple monitoring points.

Optionally, adjusting the sensitivity coefficient S _ic of the grouting area and the division of the grouting area so that the value of the inversion optimization objective function G reaches a preset range includes: adjusting the sensitivity coefficient of the grouting area so that the inversion optimization objective function G The value reaches the preset range; if the value of the inversion optimization objective function G cannot reach the preset range, adjust the partition of the grouting area so that the value of the inversion optimization objective function G reaches the preset range.

Optionally, after the zone 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 layout conditions; The value range of the equivalent grouting pressure vector of the partition is calculated according to the value range of the sensitivity coefficient of the grouting area.

As shown in Figure 4, in an optional embodiment,

1) According to the construction geological conditions, the height of the grouting section and the construction sequence, roughly divide the grouting.

2) Select an equivalent grouting pressure vector F _i within the range: A≤F _i ≤B (i=1, 2, 3, 4,..., n), where A and B are the upper and lower values of the equivalent grouting pressure The lower limit is usually based on geological exploration, acoustic wave testing, rock physical properties measured grouting pressure and other data to estimate the upper and lower limits of _Sc to obtain the upper and lower limits of F _i .

3) Boundary conditions, material parameters and F _i are brought into the model and calculated by numerical simulation methods (such as: finite element method), and F _i is applied to the plane where the bottom of the grouting hole section is located (see Figure 2 and Figure 3).

4) The component ΔH _z (X) of the lifting displacement monitored on site is compared with the numerical simulation calculation result ΔH _z ′(X), and the difference between the measured lifting displacement and the numerical simulation calculation is shown in formula (4) , X represents the Xth monitoring point participating in the differential operation.

r(X)=ΔH _z (X)-ΔH _z (X) (4)

Based on this, the objective function of inversion optimization is established:

G＝∑r ² (X)＝∑(ΔH _z (X)-ΔH _z ′(X)) ² (5)

By adjusting the S _ic and the numerical model to continuously reduce the value of the optimized objective function G, finally obtain the equivalent grouting pressure vector of the zone suitable for the specific arch dam foundation grouting, and bring the optimized S _ic into the model for grouting timing, zone grouting and The simulation of the impact of grouting pressure on safety issues such as lifting and cracking can predict the lifting value of the arch dam under different grouting construction sequences and pressure combinations.

Optionally, if the optimization function G cannot meet the accuracy requirement by adjusting S _ic , adjust the grouting partition.

According to another aspect of the embodiment of the present invention, there is also provided a grouting pressure simulation device based on the zone equivalent grouting pressure vector, including: a division module, which is used to divide multiple grouting areas according to the elevation of the grouting section and the inclination angle of the grouting pipe The setting module is used to set the sensitivity coefficient of the grouting area according to the lithology of the grouting area and the grouting construction layout conditions; the establishment module is used to establish a partition equivalent grouting pressure vector model according to the grouting area and the sensitivity coefficient of the grouting area, wherein, The expression of the partition equivalent grouting pressure vector is: Among them, N is the number of multiple grouting areas, i is the i-th grouting area in the multiple grouting areas; S _ic is the sensitivity coefficient of the i-th grouting area; S _i1 , S _i2 and S _i3 represent the i-th grouting area The area of grouting holes of sequence I, II and III in the area; n _i1 , n _i2 and n _i3 represent the number of grouting holes of sequence I, II and III in the i-th grouting region; P _i1 , P _i2 and P _i3 represent the number of grouting holes in the i-th grouting region The grouting pressure of the I, II and III sequence grouting holes in the i grouting area; r _ij is the grouting influence radius of the jth grouting hole in the i grouting area; z is the vertical upward unit vector; F _i is the grouting process The equivalent grouting pressure vector that has a lifting effect on the bedrock or dam body above the grouting area, the partition equivalent grouting pressure vector model is used to simulate the complex distribution of grouting pressure in the fracture network.

Optionally, the device further includes: an acquisition module, used to acquire the boundary conditions and material parameters of the grouting area; a first calculation module, used to partition the equivalent grouting pressure vector F _i according to the boundary conditions and material parameters, and use the numerical simulation method Calculate the grouting displacement ΔH _z ′(X), where X is the Xth monitoring point in the grouting area; the function module is used to establish the inversion optimization objective function G=∑r ² (X)=∑(ΔH _z (X)-ΔH _z ′(X)) ² , where ΔH _z (X) is the actual grouting displacement; the adjustment module is used to adjust the sensitivity coefficient s _ic of the grouting area and the division of the grouting area, so that the inversion The value of the optimization objective function G reaches the preset range; the second calculation module is used to calculate the partition equivalent grouting pressure vector according to the adjusted sensitivity coefficient of the grouting area and the grouting area when the inversion optimization objective function G reaches the preset range F _i .

Optionally, the adjustment module includes: a first adjustment unit, used to adjust the sensitivity coefficient s _ic of the grouting area, so that the value of the inversion optimization objective function G reaches a preset range; a second adjustment unit, used for Dangguo inversion When the value of the optimization objective function G cannot reach the preset range, adjust the partition of the grouting area so that the value of the inversion optimization objective function G reaches the preset range.

Optionally, the device also includes: a value module, used to determine the value range of the sensitivity coefficient of the grouting area according to the lithology of the grouting area and the layout conditions of the grouting construction; a third calculation module, used to determine the value range of the sensitivity coefficient of the grouting area according to The value range is calculated to obtain the value range of the partition equivalent grouting pressure vector.

The serial numbers of the above embodiments of the present invention are for description only, and do not represent the advantages and disadvantages of the embodiments.

In the above-mentioned embodiments of the present invention, the descriptions of each embodiment have their own emphases, and for parts not described in detail in a certain embodiment, reference may be made to relevant descriptions of other embodiments.

In the several embodiments provided in this application, it should be understood that the disclosed technical content can be realized in other ways. Wherein, the device embodiments described above are only illustrative. For example, the division of the units may be a logical function division. In actual implementation, there may be other division methods. For example, multiple units or components may be combined or may be Integrate into another system, or some features may be ignored, or not implemented. In another point, the mutual coupling or direct coupling or communication connection shown or discussed may be through some interfaces, and the indirect coupling or communication connection of units or modules may be in electrical or other forms.

The units described as separate components may or may not be physically separated, and the components displayed as units may or may not be physical units, that is, they may be located in one place, or may be distributed to multiple units. Part or all of the units can be selected according to actual needs to achieve the purpose of the solution of this embodiment.

In addition, each functional unit in each embodiment of the present invention may be integrated into one processing unit, each unit may exist separately physically, or two or more units may be integrated into one unit. The above-mentioned integrated units can be implemented in the form of hardware or in the form of software functional units.

If the integrated unit is realized in the form of a software function unit and sold or used as an independent product, it can be stored in a computer-readable storage medium. Based on this understanding, the essence of the technical solution of the present invention or the part that contributes to the prior art or all or part of the technical solution can be embodied in the form of a software product, and the computer software product is stored in a storage medium , including several instructions to make a computer device (which may be a personal computer, a server, or a network device, etc.) execute all or part of the steps of the methods described in various embodiments of the present invention. The aforementioned storage media include: U disk, read-only memory (ROM, Read-Only Memory), random access memory (RAM, Random Access Memory), mobile hard disk, magnetic disk or optical disk and other media that can store program codes. .

The above is only a preferred embodiment of the present invention, it should be pointed out that, for those of ordinary skill in the art, without departing from the principle of the present invention, some improvements and modifications can also be made, and these improvements and modifications can also be made. It should be regarded as the protection scope of the present invention.

Claims (8)

1. A grouting pressure simulation method based on partition equivalent grouting pressure vector, it is characterized in that, comprising: Divide multiple 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 arrangement conditions of the grouting construction; According to the grouting area and the sensitivity coefficient of the grouting area, a partition equivalent grouting pressure vector model is established, wherein the expression of the partition equivalent grouting pressure vector is: Among them, N is the number of the multiple grouting areas, i is the i-th grouting area in the multiple grouting areas; S _ic is the sensitivity coefficient of the i-th grouting area; S _i1 , S _i2 and S _i3 represent The area of I, II and III order grouting holes in the i-th grouting area; n _i1 , n _i2 and n _i3 represent the number of I, II and III order grouting holes in the i-th grouting area; P _i1 , P _i2 and P _i3 represents the grouting pressure of the I, II and III sequence grouting holes in the i-th grouting area; r _ij is the grouting influence radius of the j-th sequence grouting hole in the i-th grouting area; z is the vertical unit vector; F _i is the equivalent grouting pressure vector generated during the grouting process that has a lifting effect on the bedrock or dam body above the grouting area, and the partition equivalent grouting pressure vector model is used to simulate the complex distribution of grouting pressure in the fracture network. 2. the grouting pressure simulation method based on the zonal equivalent grouting pressure vector according to claim 1, is characterized in that, after setting up the zonal equivalent grouting pressure vector model according to the sensitivity coefficient of the grouting area and the grouting area, The method also includes: Obtain boundary conditions and material parameters of the grouting area; According to the boundary conditions, the material parameters and the equivalent grouting pressure vector F _i of the partition, the grouting displacement ΔH _z ′(X) is calculated by numerical simulation method, where X is the Xth grouting area of the grouting area monitoring points; Establish the inversion optimization objective function G=∑r ² (X)=∑(ΔH _z (X)-ΔH _z ′(X)) ² , where ΔH _z (X) is the actual grouting displacement; Adjusting the sensitivity coefficient S _ic of the grouting area and the division of the grouting area, so that the value of the inversion optimization objective function G reaches a preset range; When the inversion optimization objective function G reaches the preset range, the partition equivalent grouting pressure vector F _i is calculated according to the adjusted sensitivity coefficient of the grouting area and the grouting area. 3. the grouting pressure simulation method based on the partition equivalent grouting pressure vector according to claim 2, is characterized in that, adjust the sensitivity coefficient S _ic of the grouting area and the partition of the grouting area, so that the inversion Optimizing the value of the objective function G to reach the preset range includes: Adjusting the sensitivity coefficient of the grouting area so that the value of the inversion optimization objective function G reaches a preset range; If the value of the inversion optimization objective function G cannot reach a preset range, then adjust the division of the grouting area so that the value of the inversion optimization objective function G reaches a preset range. 4. the grouting pressure simulation method based on the zonal equivalent grouting pressure vector according to claim 1, is characterized in that, after setting up the zonal equivalent grouting pressure vector model according to the sensitivity coefficient of the grouting area and the grouting area, The method also includes: Determine the value range of the sensitivity coefficient of the grouting area according to the lithology of the grouting area and the grouting construction layout conditions; The value range of the equivalent grouting pressure vector of the partition is calculated according to the value range of the sensitivity coefficient of the grouting area. 5. A grouting pressure simulation device based on a partition equivalent grouting pressure vector, characterized in that it comprises: The division module is used to divide multiple grouting areas according to the elevation of the grouting section and the inclination angle of the grouting pipe; A setting module, configured to set the sensitivity coefficient of the grouting area according to the lithology of the grouting area and the arrangement conditions of the grouting construction; Establishing a module for establishing a partition equivalent grouting pressure vector model according to the grouting area and the sensitivity coefficient of the grouting area, wherein the expression of the partition equivalent grouting pressure vector is: Among them, N is the number of the multiple grouting areas, i is the i-th grouting area in the multiple grouting areas; S _ic is the sensitivity coefficient of the i-th grouting area; S _i1 , S _i2 and S _i3 represent The area of I, II and III order grouting holes in the i-th grouting area; n _i1 , n _i2 and n _i3 represent the number of I, II and III order grouting holes in the i-th grouting area; P _i1 , P _i2 and P _i3 represents the grouting pressure of the I, II and III sequence grouting holes in the i-th grouting area; r _ij is the grouting influence radius of the j-th sequence grouting hole in the i-th grouting area; z is the vertical unit vector; F _i is the equivalent grouting pressure vector generated during the grouting process that has a lifting effect on the bedrock or dam body above the grouting area, and the partition equivalent grouting pressure vector model is used to simulate the complex distribution of grouting pressure in the fracture network. 6. The grouting pressure simulation device based on the partition equivalent grouting pressure vector according to claim 5, wherein the device further comprises: An acquisition module, configured to acquire boundary conditions and material parameters of the grouting area; The first calculation module is used to calculate and obtain the grouting displacement ΔH _z '(X) by numerical simulation method according to the boundary conditions and the material parameters and the equivalent grouting pressure vector F _i of the partition, where X is The Xth monitoring point in the grouting area; Function module, used to establish inversion optimization objective function G=∑r ² (X)=∑(ΔH _z (X)-ΔH _z ′(X)) ² , where ΔH _z (X) is the actual grouting displacement ; An adjustment module, configured to adjust the sensitivity coefficient S _ic of the grouting area and the division of the grouting area, so that the value of the inversion optimization objective function G reaches a preset range; The second calculation module is used to calculate the partition equivalent grouting pressure vector F _i according to the adjusted sensitivity coefficient of the grouting area and the grouting area when the inversion optimization objective function G reaches the preset range. 7. The grouting pressure simulation device based on the partition equivalent grouting pressure vector according to claim 6, wherein the adjustment module comprises: A first adjustment unit, configured to adjust the sensitivity coefficient S _ic of the grouting area, so that the value of the inversion optimization objective function G reaches a preset range; The second adjustment unit is configured to adjust the division of the grouting area so that the value of the inversion optimization objective function G reaches a preset range when the value of the inversion optimization objective function G cannot reach a preset range. 8. The grouting pressure simulation device based on the partition equivalent grouting pressure vector according to claim 5, wherein the device further comprises: A value acquisition module, configured to determine the value range of the sensitivity coefficient of the grouting area according to the lithology of the grouting area and the arrangement conditions of the grouting construction; The third calculation module is configured to calculate the value range of the equivalent grouting pressure vector of the partition according to the value range of the sensitivity coefficient of the grouting area.