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 PDFInfo
- 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
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
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.
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)
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 |
-
2019
- 2019-07-24 CN CN201910672411.8A patent/CN110532609B/en active Active
Patent Citations (6)
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)
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 |