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

Grouting pressure analogy method and device based on the equivalent grouting pressure vector of subregion Download PDF

Info

Publication number
CN110532609A
CN110532609A CN201910672411.8A CN201910672411A CN110532609A CN 110532609 A CN110532609 A CN 110532609A CN 201910672411 A CN201910672411 A CN 201910672411A CN 110532609 A CN110532609 A CN 110532609A
Authority
CN
China
Prior art keywords
grouting
region
subregion
equivalent
grouting 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.)
Granted
Application number
CN201910672411.8A
Other languages
Chinese (zh)
Other versions
CN110532609B (en
Inventor
林鹏
魏鹏程
樊启祥
汪志林
黄纪村
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tsinghua University
China Huaneng Group Co Ltd
China Three Gorges Projects Development Co Ltd CTG
Original Assignee
Tsinghua University
China Huaneng Group Co Ltd
China Three Gorges Projects Development Co Ltd CTG
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Tsinghua University, China Huaneng Group Co Ltd, China Three Gorges Projects Development Co Ltd CTG filed Critical Tsinghua University
Priority to CN201910672411.8A priority Critical patent/CN110532609B/en
Publication of CN110532609A publication Critical patent/CN110532609A/en
Application granted granted Critical
Publication of CN110532609B publication Critical patent/CN110532609B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Consolidation Of Soil By Introduction Of Solidifying Substances Into Soil (AREA)

Abstract

The invention discloses a kind of grouting pressure analogy methods and device based on the equivalent grouting pressure vector of subregion.Wherein, this method comprises: dividing multiple grouting regions according to the tilt angle of elevation where grout compartment and grout hole;The sensitivity coefficient in the grouting region is set according to the lithology in the grouting region and grouting construction arrangement condition;According to the sensitivity coefficient in the grouting region and the grouting region, the equivalent grouting pressure vector model of subregion is established, wherein the expression formula of the equivalent grouting pressure vector of subregion are as follows:The equivalent grouting pressure vector model of subregion is used for complex distributions of the simulated grouting pressure in Fracture Networks.The present invention solves the technical issues of can not instructing engineering construction by the practical grouting pressure of numerical method real simulation.

Description

Grouting pressure analogy method and device based on the equivalent grouting pressure vector of subregion
Technical field
The present invention relates to grouting construction fields, are based on the equivalent grouting pressure vector of subregion in particular to a kind of Grouting pressure analogy method and device.
Background technique
Grouting construction process has very strong concealment, complexity and uncertainty, and the quality and effect of construction determine The success or failure of engineering.Grouting pressure is the main source of grouting energy, is control and a key factor for improving grouting quality. The simulation of the practical grouting pressure in construction site, and then the choosing of guiding construction scene grouting pressure are carried out by the method for numerical simulation Arrangement of selecting and be in the milk is of great significance.However currently without a kind of method that can more actually simulate practical grouting pressure, lead The result of numerical simulation is caused to substantially reduce the directive significance of engineering.
Technical problem present in routine techniques is mainly as follows: (1) previous grouting pressure analogy method seldom considers not With the different elevations where the grout compartment of grout hole, often approximation takes an identical elevation;(2) previous grouting pressure simulation Method seldom considers the tilt angle of grout hole, is usually approximately considered grout hole and is each perpendicular to grouting surface;(3) previous grouting Pressure simulation method seldom considers the sequence of grouting construction, i.e. I, II and III sequence grout hole and suitable with the grouting of sequence grout hole Sequence;(4) previous grouting pressure analogy method is difficult to consider the specific geological conditions in grouting region.
For above-mentioned problem, currently no effective solution has been proposed.
Summary of the invention
The embodiment of the invention provides a kind of grouting pressure analogy methods and dress based on the equivalent grouting pressure vector of subregion It sets, at least to solve the technical issues of can not instructing engineering construction by the practical grouting pressure of numerical method real simulation.
According to an aspect of an embodiment of the present invention, a kind of grouting pressure based on the equivalent grouting pressure vector of subregion is provided Power analogy method, comprising: multiple grouting regions are divided according to elevation where grout compartment and grout pipe tilt angle;According to the filling The sensitivity coefficient in the grouting region is arranged in the lithology and grouting construction arrangement condition for starching region;According to the grouting region and The sensitivity coefficient in the grouting region, establishes the equivalent grouting pressure vector model of subregion, wherein the equivalent grouting pressure vector of subregion Expression formula are as follows:Wherein, N is the quantity of the multiple grouting subregion, i For i-th of grouting region in the multiple grouting region;SicFor the sensitivity coefficient in i-th of region of being in the milk;Si1, Si2And Si3Table Show the area of I, the II and III sequence grout hole in i-th of grouting region;ni1, ni2And ni3Indicate the I, II and III of i-th of grouting region The number of sequence grout hole;Pi1, Pi2And Pi3Indicate the grouting pressure of I, the II and III sequence grout hole in i-th of grouting region;rijIt is The grouting radius of influence of jth sequence grout hole in i grouting region;Z is unit vector straight up;FiFor filling process generation The equivalent grouting pressure vector to the above basement rock in grouting region or dam body with lifting effect, the equivalent grouting pressure of subregion Vector model is used for complex distributions of the simulated grouting pressure in Fracture Networks.
Further, according to the sensitivity coefficient in the grouting region and the grouting region, the equivalent grouting of subregion is established After pressure vector model, the method also includes: obtain the boundary condition and material parameter in the grouting region;According to described Boundary condition, the material parameter and the equivalent grouting pressure vector F of the subregioni, it is calculated by method for numerical simulation Grouting lifts displacement Hz' (X), wherein X is the X monitoring point in the grouting region;Establish inverting optimization object function G =∑ r2(X)=∑ (Δ Hz(X)-ΔHz′(X))2, wherein Δ Hz(X) displacement is lifted for practical grouting;Adjust the grout compartment The sensitivity coefficient S in domainicWith it is described grouting region subregion so that the value of the inverting optimization object function G reaches preset range; When the inverting optimization object function G reaches the preset range, according to the sensitivity coefficient in the grouting region after tune and grouting Region calculates the equivalent grouting pressure vector F of subregioni
Further, the sensitivity coefficient S in the grouting region is adjustedicWith the subregion in the grouting region so that described anti- The value for drilling optimization object function G reaches the sensitivity coefficient that preset range includes: the adjustment grouting region, so that the inverting is excellent The value for changing objective function G reaches preset range;If the value of the inverting optimization object function G is unable to reach preset range, The subregion in the grouting region is adjusted so that the value of the inverting optimization object function G reaches preset range.
Further, according to the sensitivity coefficient in the grouting region and the grouting region, the equivalent grouting of subregion is established After pressure vector model, the method also includes: it is determined according to the lithology in the grouting region and grouting construction arrangement condition The value range of the sensitivity coefficient in the grouting region;It is calculated according to the value range of the sensitivity coefficient in the grouting region The value range of the equivalent grouting pressure vector of subregion.
Other side according to an embodiment of the present invention additionally provides a kind of filling based on the equivalent grouting pressure vector of subregion Starch pressure simulation device, comprising: division module, for dividing multiple fillings according to elevation where grout compartment and grout pipe tilt angle Starch region;Setup module, for the grouting region to be arranged according to the lithology and grouting construction arrangement condition in the grouting region Sensitivity coefficient;Module is established, for the sensitivity coefficient according to the grouting region and the grouting region, it is equivalent to establish subregion Grouting pressure vector model, wherein the expression formula of the equivalent grouting pressure vector of subregion are as follows: Wherein, N is the quantity of the multiple grouting subregion, and i is i-th of grouting region in the multiple grouting region;SicIt is i-th The sensitivity coefficient in grouting region;Si1, Si2And Si3Indicate the area of I, II and III the sequence grout hole in i-th of grouting region;ni1, ni2And ni3Indicate the number of I, II and III the sequence grout hole in i-th of grouting region;Pi1, Pi2And Pi3Indicate i-th of grouting region I, II and III sequence grout hole grouting pressure;rijFor the grouting radius of influence of jth sequence grout hole in i-th of grouting region;z For unit vector straight up;FiThere is lifting effect to the above basement rock in grouting region or dam body for what filling process generated Equivalent grouting pressure vector, the equivalent grouting pressure vector model of subregion are used for simulated grouting pressure answering in Fracture Networks Miscellaneous distribution.
Further, described device further include: module is obtained, for obtaining the boundary condition and material in the grouting region Parameter;First computing module is used for according to the boundary condition, the material parameter with the equivalent grouting pressure vector of the subregion Fi, grouting is calculated by method for numerical simulation and lifts displacement Hz' (X), wherein X is the X prison in the grouting region Measuring point;Function module, for establishing inverting optimization object function G=∑ r2(X)=∑ (Δ Hz(X)-ΔHz′(X))2, wherein ΔHz(X) displacement is lifted for practical grouting;Module is adjusted, for adjusting the sensitivity coefficient s in the grouting regionicWith the grouting The subregion in region, so that the value of the inverting optimization object function G reaches preset range;Second computing module, for when described When inverting optimization object function G reaches the preset range, according to the sensitivity coefficient in the grouting region after tune and grouting region meter The equivalent grouting pressure vector F in point counting areai
Further, the adjustment module includes: the first adjustment unit, for adjusting the sensitivity coefficient in the grouting region Sic, so that the value of the inverting optimization object function G reaches preset range;Second adjustment unit, it is excellent for working as inverting described in fruit When the value of change objective function G is unable to reach preset range, the subregion in the grouting region is adjusted so that the inverting optimization aim The value of function G reaches preset range.
Further, described device further include: value module, for the lithology and grouting construction according to the grouting region Arrangement condition determines the value range of the sensitivity coefficient in the grouting region;Third computing module, for according to the grout compartment The value range of the equivalent grouting pressure vector of the subregion is calculated in the value range of the sensitivity coefficient in domain.
In embodiments of the present invention, multiple grout compartments are divided using according to elevation where grout compartment and grout pipe tilt angle Domain;The sensitivity coefficient in the grouting region is set according to the lithology in the grouting region and grouting construction arrangement condition;According to institute The sensitivity coefficient for stating grouting region and the grouting region, establishes the equivalent grouting pressure vector model of subregion, wherein subregion etc. Imitate the expression formula of grouting pressure vector are as follows:Wherein, N is the multiple filling The quantity of subregion is starched, i is i-th of grouting region in the multiple grouting region;SicFor the sensitive system in i-th of region of being in the milk Number;Si1, Si2And Si3Indicate the area of I, II and III the sequence grout hole in i-th of grouting region;ni1, ni2And ni3Indicate i-th of filling Starch the number of I, II and III the sequence grout hole in region;Pi1, Pi2And Pi3Indicate I, II and III the sequence grouting in i-th of grouting region The grouting pressure in hole;rijFor the grouting radius of influence of jth sequence grout hole in i-th of grouting region;Z be unit straight up to Amount;FiIt is filling process generation to the equivalent grouting pressure vector of the above basement rock in region or dam body with lifting effect that is in the milk Mode is used for complex distributions of the simulated grouting pressure in Fracture Networks by the equivalent grouting pressure vector model of subregion, reaches The purpose that the simulation of the practical grouting pressure in construction site is carried out using the method for numerical simulation, to realize true simulation The technical effect of practical grouting pressure, and then the technical issues of solve.
Detailed description of the invention
The drawings described herein are used to provide a further understanding of the present invention, constitutes part of this application, this hair Bright illustrative embodiments and their description are used to explain the present invention, and are not constituted improper limitations of the present invention.In the accompanying drawings:
Fig. 1 is a kind of grouting pressure analogy method based on the equivalent grouting pressure vector of subregion according to an embodiment of the present invention Flow chart;
Fig. 2 is that single equivalent grouting pressure in area according to an embodiment of the present invention assumes schematic diagram;
Fig. 3 is the equivalent grouting pressure schematic diagram of subregion according to an embodiment of the present invention;
Fig. 4 is a kind of equivalent grouting pressure vector applicating flow chart of subregion according to an embodiment of the present invention.
Specific embodiment
In order to enable those skilled in the art to better understand the solution of the present invention, below in conjunction in the embodiment of the present invention Attached drawing, technical scheme in the embodiment of the invention is clearly and completely described, it is clear that described embodiment is only The embodiment of a part of the invention, instead of all the embodiments.Based on the embodiments of the present invention, ordinary skill people The model that the present invention protects all should belong in member's every other embodiment obtained without making creative work It encloses.
It should be noted that description and claims of this specification and term " first " in above-mentioned attached drawing, " Two " etc. be to be used to distinguish similar objects, without being used to describe a particular order or precedence order.It should be understood that using in this way Data be interchangeable under appropriate circumstances, so as to the embodiment of the present invention described herein can in addition to illustrating herein or Sequence other than those of description is implemented.In addition, term " includes " and " having " and their any deformation, it is intended that cover Cover it is non-exclusive include, for example, the process, method, system, product or equipment for containing a series of steps or units are not necessarily limited to Step or unit those of is clearly listed, but may include be not clearly listed or for these process, methods, product Or other step or units that equipment is intrinsic.
According to embodiments of the present invention, a kind of side of grouting pressure simulation based on the equivalent grouting pressure vector of subregion is provided Method embodiment, it should be noted that step shown in the flowchart of the accompanying drawings can be in such as a group of computer-executable instructions Computer system in execute, although also, logical order is shown in flow charts, in some cases, can be with Shown or described step is executed different from sequence herein.
Fig. 1 is a kind of grouting pressure simulation side based on the equivalent grouting pressure vector of subregion according to an embodiment of the present invention Method, as shown in Figure 1, this method comprises the following steps:
Step S102 divides multiple grouting regions according to elevation where grout compartment and grout pipe tilt angle;
Step S104, according to the sensitivity coefficient of the lithology in grouting region and grouting construction arrangement condition setting grouting region;
Step S106 establishes the equivalent grouting pressure vector of subregion according to the sensitivity coefficient in grouting region and region of being in the milk Model, wherein the expression formula of the equivalent grouting pressure vector of subregion are as follows:Its In, N is the quantity of multiple grouting subregions, and i is i-th of grouting region in multiple grouting regions;SicFor i-th of region of being in the milk Sensitivity coefficient;Si1, Si2And Si3Indicate the area of I, II and III the sequence grout hole in i-th of grouting region;ni1, ni2And ni3It indicates The number of I, II and III the sequence grout hole in i-th of grouting region;Pi1, Pi2And Pi3Indicate I, II and the III in i-th of grouting region The grouting pressure of sequence grout hole;rijFor the grouting radius of influence of jth sequence grout hole in i-th of grouting region;Z is straight up Unit vector;FiThe equivalent grouting pressure to the above basement rock in grouting region or dam body with lifting effect generated for filling process Vector, the equivalent grouting pressure vector model of subregion are used for complex distributions of the simulated grouting pressure in Fracture Networks.
The definition and expression of the equivalent grouting pressure of subregion of this invention once are described first here.Referring to fig. 2, Fig. 3.
Optionally, the foundation of the equivalent grouting pressure vector of subregion is based on following hypothesis:
Assuming that 1: ignoring by the weight of filling body, only consider grouting pressure, the A of grout hole, the upper lift in the face B are tilted in Fig. 2 It is equal in magnitude with pressing force, take A, the face B (is directed to vertical direction, i.e., perpendicular to building base as shown in formula (1) and Fig. 2 up and down The direction in face);
Wherein L1For the length of any line of apsides in the face A, E is the length of semicircle, and a, b are any pair of arch of opposite (180 °) Line, β are grout hole and vertical direction angle.
Assuming that 2: being contacted between grout pipe and basement rock well, be not in opposite sliding and grout leaking;
Assuming that 3: there is no the big Geological Defects similar to solution cavity in grout compartment;
Assuming that 4: lifting may be the region more than this crack because caused by slurries in grout compartment in any crack To lift region, the top and bottom vertical direction pressure for lifting any one crack in region can be with self-balancing, i.e. formula (2) and figure Shown in 2:
Wherein x is the crack sum for lifting region, ciIt is numbered for crack,WithFor crack ciIt erects top and bottom Histogram is to joint efforts.
The definition and expression of the equivalent grouting pressure of subregion:
Equivalent grouting pressure acting surface is assumed to obtain grout compartment by equilibrium condition in the bottommost of grout compartment Equation (3).Equation (4) and (3) are of equal value, are the expression formula of the equivalent grouting pressure vector of subregion.Equation (5) is the equivalent grouting of subregion The scalar expression of pressure,The referred to as equivalent grouting pressure vector F of subregioniAverage value.
Fi=[f1, f2, f3..., fN] (4)
Optionally, according to the sensitivity coefficient in grouting region and region of being in the milk, the equivalent grouting pressure vector mould of subregion is established After type, method further include: obtain the boundary condition and material parameter in grouting region;According to boundary condition, material parameter and The equivalent grouting pressure vector F of subregioni, grouting is calculated by method for numerical simulation and lifts displacement Hz' (X), wherein X is to fill Starch the X monitoring point in region;Establish inverting optimization object functionIts In, Δ Hz(X) displacement is lifted for practical grouting;The sensitivity coefficient s in adjustment grouting regionicWith grouting region subregion so that instead The value for drilling optimization object function G reaches preset range;When inverting optimization object function G reaches preset range, after tune The sensitivity coefficient in grouting region and grouting region calculate the equivalent grouting pressure vector F of subregioni.Optionally, it is in the milk the boundary in region Condition is the Normal Constraint in the grouting region, three-dimensional constraint etc.;Material parameter is by the material parameter of filling body (such as rock);Grouting Region includes multiple monitoring points.
Optionally, the sensitivity coefficient S in adjustment grouting regionicWith grouting region subregion so that inverting optimization object function G Value reach preset range include: adjustment grouting region sensitivity coefficient so that the value of inverting optimization object function G reaches default Range;If the value of inverting optimization object function G is unable to reach preset range, the subregion in grouting region is adjusted so that inverting is excellent The value for changing objective function G reaches preset range.
Optionally, according to the sensitivity coefficient in grouting region and region of being in the milk, the equivalent grouting pressure vector mould of subregion is established After type, method further include: the sensitivity coefficient in grouting region is determined according to the lithology in grouting region and grouting construction arrangement condition Value range;The value of the equivalent grouting pressure vector of subregion is calculated according to the value range of the sensitivity coefficient in grouting region Range.
As shown in figure 4, in an alternative embodiment,
1) rough grouting subregion is carried out according to construction geology condition, grout compartment height and sequence of construction.
2) an equivalent grouting pressure vector F is choseni, range are as follows: A≤Fi≤ B (i=1,2,3,4 ..., n), wherein A and B is the bound of equivalent grouting pressure, generally according to geological prospecting, sonic test, rock physical property actual measurement grouting pressure etc. Data estimates ScBound so that obtain FiBound.
3) by boundary condition, material parameter and FiIt is brought into model through method for numerical simulation (such as: finite element side Method) it is calculated, FiPlane where being applied to grouting hole section bottom (see Fig. 2 and Fig. 3).
4) by the component Δ H for lifting displacement of field monitoringz(X) with numerical simulation calculation result Δ HzIt is poor that ' (X) makees, actual measurement Shown in the difference for the lifting displacement such as formula (4) obtained with numerical simulation calculation, X represents the X monitoring point for participating in making difference.
R (X)=Δ Hz(X)-ΔHz(X) (4)
The objective function of inverting optimization is established accordingly:
G=∑ r2(X)=∑ (Δ Hz(X)-ΔHz′(X))2 (5)
By adjusting SicConstantly reduce optimization object function value G with numerical model, finally obtains suitable for specific arch dam The equivalent grouting pressure vector of subregion of dam foundation grouting, by the S after optimizationicBring into model carry out grouting opportunity, subregion grouting and Simulation of the grouting pressure size on the safety problems such as lifting, cracking and influence, can be predicted different grouting construction sequences, force combination etc. Under the conditions of arch dam lift value.
Optionally, if by adjusting SicMajorized function G can not be made to reach precise requirements, then adjust grouting subregion.
Other side according to an embodiment of the present invention additionally provides a kind of filling based on the equivalent grouting pressure vector of subregion Starch pressure simulation device, comprising: division module, for dividing multiple fillings according to elevation where grout compartment and grout pipe tilt angle Starch region;Setup module, for being according to the lithology and the sensitive of grouting construction arrangement condition setting grouting region in grouting region Number;Module is established, for the sensitivity coefficient according to grouting region and region of being in the milk, establishes the equivalent grouting pressure vector mould of subregion Type, wherein the expression formula of the equivalent grouting pressure vector of subregion are as follows:Its In, N is the quantity of multiple grouting subregions, and i is i-th of grouting region in multiple grouting regions;SicFor i-th of region of being in the milk Sensitivity coefficient;Si1, Si2And Si3Indicate the area of I, II and III the sequence grout hole in i-th of grouting region;ni1, ni2And ni3It indicates The number of I, II and III the sequence grout hole in i-th of grouting region;Pi1, Pi2And Pi3Indicate I, II and the III in i-th of grouting region The grouting pressure of sequence grout hole;rijFor the grouting radius of influence of jth sequence grout hole in i-th of grouting region;Z is straight up Unit vector;FiThe equivalent grouting pressure to the above basement rock in grouting region or dam body with lifting effect generated for filling process Vector, the equivalent grouting pressure vector model of subregion are used for complex distributions of the simulated grouting pressure in Fracture Networks.
Optionally, device further include: module is obtained, for obtaining the boundary condition and material parameter in grouting region;First Computing module is used for according to boundary condition, material parameter with the equivalent grouting pressure vector F of subregioni, pass through method for numerical simulation meter Calculation obtains grouting and lifts displacement Hz' (X), wherein X is the X monitoring point in grouting region;Function module, it is anti-for establishing Drill optimization object function G=∑ r2(X)=∑ (Δ Hz(X)-ΔHz′(X))2, wherein Δ Hz(X) displacement is lifted for practical grouting; Module is adjusted, for adjusting the sensitivity coefficient s in grouting regionicWith the subregion in grouting region so that inverting optimization object function G Value reaches preset range;Second computing module, for when inverting optimization object function G reaches preset range, after tune The sensitivity coefficient in grouting region and grouting region calculate the equivalent grouting pressure vector F of subregioni
Optionally, adjustment module includes: the first adjustment unit, for adjusting the sensitivity coefficient s in grouting regionic, so that instead The value for drilling optimization object function G reaches preset range;Second adjustment unit, for when fruit inverting optimization object function G value without When method reaches preset range, the subregion in adjustment grouting region is so that the value of inverting optimization object function G reaches preset range.
Optionally, device further include: value module, for true according to the lithology and grouting construction arrangement condition in grouting region Surely the value range of the sensitivity coefficient in grouting region;Third computing module, for according to grouting region sensitivity coefficient value Range computation obtains the value range of the equivalent grouting pressure vector of subregion.
The serial number of the above embodiments of the invention is only for description, does not represent the advantages or disadvantages of the embodiments.
In the above embodiment of the invention, it all emphasizes particularly on different fields to the description of each embodiment, does not have in some embodiment The part of detailed description, reference can be made to the related descriptions of other embodiments.
In several embodiments provided herein, it should be understood that disclosed technology contents can pass through others Mode is realized.Wherein, the apparatus embodiments described above are merely exemplary, such as the division of the unit, Ke Yiwei A kind of logical function partition, there may be another division manner in actual implementation, for example, multiple units or components can combine or Person is desirably integrated into another system, or some features can be ignored or not executed.Another point, shown or discussed is mutual Between coupling, direct-coupling or communication connection can be through some interfaces, the INDIRECT COUPLING or communication link of unit or module It connects, can be electrical or other forms.
The unit as illustrated by the separation member may or may not be physically separated, aobvious as unit The component shown may or may not be physical unit, it can and it is in one place, or may be distributed over multiple On unit.It can some or all of the units may be selected to achieve the purpose of the solution of this embodiment according to the actual needs.
It, can also be in addition, the functional units in various embodiments of the present invention may be integrated into one processing unit It is that each unit physically exists alone, can also be integrated in one unit with two or more units.Above-mentioned integrated list Member both can take the form of hardware realization, can also realize in the form of software functional units.
If the integrated unit is realized in the form of SFU software functional unit and sells or use as independent product When, it can store in a computer readable storage medium.Based on this understanding, technical solution of the present invention is substantially The all or part of the part that contributes to existing technology or the technical solution can be in the form of software products in other words It embodies, which is stored in a storage medium, including some instructions are used so that a computer Equipment (can for personal computer, server or network equipment etc.) execute each embodiment the method for the present invention whole or Part steps.And storage medium above-mentioned includes: that USB flash disk, read-only memory (ROM, Read-Only Memory), arbitrary access are deposited Reservoir (RAM, Random Access Memory), mobile hard disk, magnetic or disk etc. be various to can store program code Medium.
The above is only a preferred embodiment of the present invention, it is noted that for the ordinary skill people of the art For member, various improvements and modifications may be made without departing from the principle of the present invention, these improvements and modifications are also answered It is considered as protection scope of the present invention.

Claims (8)

1. a kind of grouting pressure analogy method based on the equivalent grouting pressure vector of subregion characterized by comprising
Multiple grouting regions are divided according to elevation where grout compartment and grout pipe tilt angle;
The sensitivity coefficient in the grouting region is set according to the lithology in the grouting region and grouting construction arrangement condition;
According to the sensitivity coefficient in the grouting region and the grouting region, the equivalent grouting pressure vector model of subregion is established, Wherein, the expression formula of the equivalent grouting pressure vector of subregion are as follows:Wherein, N For the quantity of the multiple grouting subregion, i is i-th of grouting region in the multiple grouting region;SicFor i-th of grout compartment The sensitivity coefficient in domain;Si1, Si2And Si3Indicate the area of I, II and III the sequence grout hole in i-th of grouting region;ni1, ni2And ni3 Indicate the number of I, II and III the sequence grout hole in i-th of grouting region;Pi1, Pi2And Pi3Indicate the I, II in i-th of grouting region With the grouting pressure of III sequence grout hole;rijFor the grouting radius of influence of jth sequence grout hole in i-th of grouting region;Z is vertical Upward unit vector;FiThe equivalent filling to the above basement rock in grouting region or dam body with lifting effect generated for filling process Pressure vector is starched, the equivalent grouting pressure vector model of subregion is for complexity of the simulated grouting pressure in Fracture Networks point Cloth.
2. the grouting pressure analogy method according to claim 1 based on the equivalent grouting pressure vector of subregion, feature exist In, according to it is described grouting region and it is described grouting region sensitivity coefficient, establish the equivalent grouting pressure vector model of subregion it Afterwards, the method also includes:
Obtain the boundary condition and material parameter in the grouting region;
According to the boundary condition, the material parameter and the equivalent grouting pressure vector F of the subregioni, by numerical simulation side Method is calculated grouting and lifts displacement Hz' (X), wherein X is the X monitoring point in the grouting region;
Establish inverting optimization object function G=∑ r2(X)=∑ (Δ Hz(X)-ΔHz′(X))2, wherein Δ Hz(X) it is filled to be practical Slurry lifts displacement;
Adjust the sensitivity coefficient S in the grouting regionicWith it is described grouting region subregion so that the inverting optimization object function The value of G reaches preset range;
When the inverting optimization object function G reaches the preset range, according to after tune grouting region sensitivity coefficient and Be in the milk the region calculating equivalent grouting pressure vector F of subregioni
3. the grouting pressure analogy method according to claim 2 based on the equivalent grouting pressure vector of subregion, feature exist In the sensitivity coefficient S in the adjustment grouting regionicWith it is described grouting region subregion so that the inverting optimization object function G Value reach preset range and include:
The sensitivity coefficient for adjusting the grouting region, so that the value of the inverting optimization object function G reaches preset range;
If the value of the inverting optimization object function G is unable to reach preset range, adjust it is described grouting region subregion with The value of the inverting optimization object function G is set to reach preset range.
4. the grouting pressure analogy method according to claim 1 based on the equivalent grouting pressure vector of subregion, feature exist In, according to it is described grouting region and it is described grouting region sensitivity coefficient, establish the equivalent grouting pressure vector model of subregion it Afterwards, the method also includes:
The value of the sensitivity coefficient in the grouting region is determined according to the lithology in the grouting region and grouting construction arrangement condition Range;
Taking for the equivalent grouting pressure vector of the subregion is calculated according to the value range of the sensitivity coefficient in the grouting region It is worth range.
5. a kind of grouting pressure simulator based on the equivalent grouting pressure vector of subregion characterized by comprising
Division module, for dividing multiple grouting regions according to elevation where grout compartment and grout pipe tilt angle;
Setup module, for the quick of the grouting region to be arranged according to the lithology and grouting construction arrangement condition in the grouting region Feel coefficient;
Module is established, for the sensitivity coefficient according to the grouting region and the grouting region, establishes the equivalent grouting of subregion Pressure vector model, wherein the expression formula of the equivalent grouting pressure vector of subregion are as follows: Wherein, N is the quantity of the multiple grouting subregion, and i is i-th of grouting region in the multiple grouting region;SicIt is i-th The sensitivity coefficient in grouting region;Si1, Si2And Si3Indicate the area of I, II and III the sequence grout hole in i-th of grouting region;ni1, ni2And ni3Indicate the number of I, II and III the sequence grout hole in i-th of grouting region;Pi1, Pi2And Pi3Indicate i-th of grouting region I, II and III sequence grout hole grouting pressure;rijFor the grouting radius of influence of jth sequence grout hole in i-th of grouting region;z For unit vector straight up;FiThere is lifting effect to the above basement rock in grouting region or dam body for what filling process generated Equivalent grouting pressure vector, the equivalent grouting pressure vector model of subregion are used for simulated grouting pressure answering in Fracture Networks Miscellaneous distribution.
6. the grouting pressure simulator according to claim 5 based on the equivalent grouting pressure vector of subregion, feature exist In described device further include:
Module is obtained, for obtaining the boundary condition and material parameter in the grouting region;
First computing module is used for according to the boundary condition, the material parameter with the equivalent grouting pressure vector of the subregion Fi, grouting is calculated by method for numerical simulation and lifts displacement Hz' (X), wherein X is the X prison in the grouting region Measuring point;
Function module, for establishing inverting optimization object function G=∑ r2(X)=∑ (Δ Hz(X)-ΔHz′(X))2, wherein Δ Hz(X) displacement is lifted for practical grouting;
Module is adjusted, for adjusting the sensitivity coefficient S in the grouting regionicWith the subregion in the grouting region so that described anti- The value for drilling optimization object function G reaches preset range;
Second computing module, for when the inverting optimization object function G reaches the preset range, according to the grouting after tune The sensitivity coefficient in region and grouting region calculate the equivalent grouting pressure vector F of subregioni
7. the grouting pressure simulator according to claim 6 based on the equivalent grouting pressure vector of subregion, feature exist In the adjustment module includes:
The first adjustment unit, for adjusting the sensitivity coefficient S in the grouting regionic, so that the inverting optimization object function G Value reaches preset range;
Second adjustment unit, when the value for the inverting optimization object function G described in the fruit is unable to reach preset range, described in adjustment The subregion in grouting region is so that the value of the inverting optimization object function G reaches preset range.
8. the grouting pressure simulator according to claim 5 based on the equivalent grouting pressure vector of subregion, feature exist In described device further include:
Value module, for determining the quick of the grouting region according to the lithology and grouting construction arrangement condition in the grouting region Feel the value range of coefficient;
Third computing module, for that the subregion is calculated to be equivalent according to the value range of the sensitivity coefficient in the grouting region The value range of grouting pressure vector.
CN201910672411.8A 2019-07-24 2019-07-24 Grouting pressure simulation method and device based on partition equivalent grouting pressure vector Active CN110532609B (en)

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 true CN110532609A (en) 2019-12-03
CN110532609B 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 (6)

* Cited by examiner, † Cited by third party
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
CN102704497A (en) * 2012-05-31 2012-10-03 中冶集团武汉勘察研究院有限公司 RPW (ranked position weight) process for treating high fill slope
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
CN108316347A (en) * 2018-03-09 2018-07-24 中冶京诚工程技术有限公司 Partition waterproof assembly type underground comprehensive pipe gallery and waterproof method thereof

Patent Citations (6)

* Cited by examiner, † Cited by third party
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
CN102704497A (en) * 2012-05-31 2012-10-03 中冶集团武汉勘察研究院有限公司 RPW (ranked position weight) process for treating high fill slope
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
CN108316347A (en) * 2018-03-09 2018-07-24 中冶京诚工程技术有限公司 Partition waterproof assembly type underground comprehensive pipe gallery and waterproof method thereof

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
PENG LIN等: "Study on Optimal Grouting Timing for Controlling Uplift Deformation of a Super High Arch Dam", 《ROCK MECH ROCK ENG》 *
周正军等: "覆盖层地基注浆扩散数值模拟方法及应用", 《地下空间与工程学报》 *
朱晓旭: "特高拱坝开裂分析方法与应用", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 *

Also Published As

Publication number Publication date
CN110532609B (en) 2021-02-19

Similar Documents

Publication Publication Date Title
US20180293789A1 (en) Three-dimensional geomechanical modeling of casing deformation for hydraulic fracturing treatment design
US10571601B2 (en) Geologic fault representation for flow simulation
CN103472483B (en) Velocity modeling method based on true earth's surface or floating datum
CN108089227B (en) Novel stratum pore pressure prediction method based on three-dimensional seismic data
MXPA02001002A (en) Method for predicting quantitative values of a rock or fluid property in a reservoir using seismic data.
CN110838175B (en) Geological model building method for gas injection development oil reservoir
CN107219159B (en) Method and device for determining reservoir heterogeneity
AU2012375233B2 (en) System and method for automatic local grid refinement in reservoir simulation systems
US10684392B2 (en) Method to generate the in-situ state of stress in a domain Ω in a geological structure
CN108345719A (en) A kind of high angle hole/horizontal well is with boring inversion method after electromagnetic wave logging data is bored
CN106338765B (en) Method and system for determining the composite elastic parameter oil-gas recognition factor
WO2017127610A1 (en) Generating dynamically calibrated geo-models in green fields
CN107607996A (en) Based on phased sequential co-simulation Geological Modeling
CN111980662A (en) Method for rapidly processing anisotropic stratum array lateral logging data of inclined shaft
NO20190509A1 (en) Global surface paleo-temperature modeling tool
CN107609265A (en) A kind of formation stress field Finite Element Method and system based on ant tracking
Eyike et al. Lithospheric structure of the West and Central African Rift System from regional three-dimensional gravity modelling
Petersen et al. The subduction dichotomy of strong plates and weak slabs
CN105447319B (en) Arch Dam Abutment sliding block difference geology level meters calculate determination method for parameter
CN110532609A (en) Grouting pressure analogy method and device based on the equivalent grouting pressure vector of subregion
CN113591328A (en) Deep shale hydraulic fracture long-term support opening change calculation method
CA3019487C (en) Systems and methods for unfaulting point clouds
CN106321082A (en) Method and system for obtaining maximum external casing load of salt gypsum rock stratum
Manchanda et al. Creation of multi-stranded fractures in homogeneous rocks
WO2021047752A1 (en) Methods and systems for well-to-cell coupling in reservoir simulation

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