CN107219560B - Mine worked-out section deflection and stability assessment method based on gravity anomaly - Google Patents
Mine worked-out section deflection and stability assessment method based on gravity anomaly Download PDFInfo
- Publication number
- CN107219560B CN107219560B CN201710388355.6A CN201710388355A CN107219560B CN 107219560 B CN107219560 B CN 107219560B CN 201710388355 A CN201710388355 A CN 201710388355A CN 107219560 B CN107219560 B CN 107219560B
- Authority
- CN
- China
- Prior art keywords
- out section
- spheroid
- country rock
- mine worked
- barycenter
- 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.)
- Expired - Fee Related
Links
- 230000005484 gravity Effects 0.000 title claims abstract description 116
- 238000000034 method Methods 0.000 title claims abstract description 64
- 238000013097 stability assessment Methods 0.000 title claims abstract description 19
- 239000011435 rock Substances 0.000 claims abstract description 195
- 238000005065 mining Methods 0.000 claims abstract description 35
- 238000005070 sampling Methods 0.000 claims description 50
- PCTMTFRHKVHKIS-BMFZQQSSSA-N (1s,3r,4e,6e,8e,10e,12e,14e,16e,18s,19r,20r,21s,25r,27r,30r,31r,33s,35r,37s,38r)-3-[(2r,3s,4s,5s,6r)-4-amino-3,5-dihydroxy-6-methyloxan-2-yl]oxy-19,25,27,30,31,33,35,37-octahydroxy-18,20,21-trimethyl-23-oxo-22,39-dioxabicyclo[33.3.1]nonatriaconta-4,6,8,10 Chemical compound C1C=C2C[C@@H](OS(O)(=O)=O)CC[C@]2(C)[C@@H]2[C@@H]1[C@@H]1CC[C@H]([C@H](C)CCCC(C)C)[C@@]1(C)CC2.O[C@H]1[C@@H](N)[C@H](O)[C@@H](C)O[C@H]1O[C@H]1/C=C/C=C/C=C/C=C/C=C/C=C/C=C/[C@H](C)[C@@H](O)[C@@H](C)[C@H](C)OC(=O)C[C@H](O)C[C@H](O)CC[C@@H](O)[C@H](O)C[C@H](O)C[C@](O)(C[C@H](O)[C@H]2C(O)=O)O[C@H]2C1 PCTMTFRHKVHKIS-BMFZQQSSSA-N 0.000 claims description 21
- 238000004364 calculation method Methods 0.000 claims description 3
- 238000013507 mapping Methods 0.000 claims description 3
- 238000012544 monitoring process Methods 0.000 abstract description 10
- 238000013461 design Methods 0.000 abstract description 5
- 238000001514 detection method Methods 0.000 description 6
- 238000004458 analytical method Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 239000003245 coal Substances 0.000 description 3
- 230000006378 damage Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000009412 basement excavation Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 229910052500 inorganic mineral Inorganic materials 0.000 description 2
- 230000010354 integration Effects 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 239000011707 mineral Substances 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000011084 recovery Methods 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000011439 discrete element method Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 239000011229 interlayer Substances 0.000 description 1
- 239000010410 layer Substances 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 230000036244 malformation Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000000717 retained effect Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
- 238000012876 topography Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V7/00—Measuring gravitational fields or waves; Gravimetric prospecting or detecting
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
- Testing Or Calibration Of Command Recording Devices (AREA)
Abstract
The invention discloses a kind of mine worked-out section deflection and stability assessment method based on gravity anomaly, comprise the following steps：Step 1: the arrangement of absolute gravimeter observation station and the erection of absolute gravimeter；Step 2: the acquisition of absolute gravity observation value；3rd, rectangular coordinate system in space is established；4th, the center-of-mass coordinate of mine worked-out section country rock spheroid calculates；5th, the drafting of the barycenter change curve of mine worked-out section country rock spheroid and the acquisition of steady-state deformation amount.Rapid simple, the reasonable in design and cost of the inventive method is low, using easy to operate, the deflection in mining area can be monitored exactly using gravity anomaly, monitoring result is accurate, and the stability in mining area is assessed according to the deflection in minery, and the steady-state deformation amount in mining area is obtained, it is practical.
Description
Technical field
The invention belongs to mining engineering and technology field, is deformed more particularly, to a kind of mine worked-out section based on gravity anomaly
Amount and stability assessment method.
Background technology
Mine worked-out section is a series of underground interlayer hole region for referring to that underground mineral is left after man-made recoveries.
Natural rock mass is in natural original stress equilibrium state originally, but because the exploitation in mine, exploitation work will necessarily be in rock mass
The miscellaneous space of inner excavation, such as tunnel, stope, so necessarily destroy mining area subsurface rock layer stress equilibrium and
Produce secondary stress-strain effect.Just because of this secondary stress-strain effect so that the country rock of mine worked-out section
Body structure is deformed, or even stress equilibrium structural damage, movement, and until region, Lithospheric Structure stress equilibrium is adjusted to new
Balance untill.
The State-owned mining enterprise of eighties of last century, due to the limitation of mining technique and equipment, largely using open-stope method, avalanche
The mining methods such as method, shrinkage method, although recovery process passes through regular design, goaf all information, because completing filling
Goaf cost is higher, and many mines also fail in time all goafs be completed to fill, and have more than 40% goaf still
Each mine is retained in, the exploitation of more leisure opinion Mineral Resources in China was also subject to the serious stage out of control, and the people adopt, wildness of mining massively,
Grow thickly mine worked-out section.The emphasis colliery in China, average excavation collapse area account for the 1/10 of mining area area containing coal.At present, according to
Incomplete statistics, in more than 20 individual province of China, autonomous region, large-scale goaf geological disaster more than 200 occurs altogether, it is accumulative to collapse
Fall into area more than 700,000 hectares, caused by direct economic loss be even more exceeded 50,000,000,000 yuan.As can be seen that due to goaf
Buckling deformation and the social harm of Mine Geological Disasters that triggers is very big, the lighter causes the local stopping production of production, and severe one causes
Casualties and huge economic losses, seriously affect country, the assets and people's property safety in area.
In addition, the problem of mine worked-out section stability is one extremely complex, it not only with exploitation ore body it is buried, incline
Angle, thickness, the lithology of overlying rock, occurrence status, thickness, physico-mechanical properties, geological structure, site topography landforms, water, text
The geological mining condition such as geological conditions and coal mining area, roof controlling method, the exploitation frequency, extraction technique is relevant, and
It is and also close related to the type, size, position of upper load etc..
The deformation measurement on outer mine worked-out section mainly has with stability assessment analysis method at present：
Firstth, Method for Numerical：Mine worked-out section numerical simulation can use boundary element method, Finite Element, limited
Calculus of finite differences and discrete element method etc..Wherein, finite element method is adapted to widest in area, and development is more ripe.
Secondth, analytic method：Analytic method is that mine worked-out section chamber is simplified, and establishes corresponding geological model, then press
A preferable mathematics physics model is abstracted as according to certain principle or assumed condition, is then asked according to numerical method
Solution.
3rd, prediction method：Prediction method mainly by calculate mine worked-out section top plate bearing capacity, residual cavity stability,
Remaining earth's surface deflection and destruction of surface scope are carried out.Coal Mining Subsidence estimated theoretical and method are more, extensive in China
The probability integration process used,
4th, half prediction semi analytical method：This method is the combination of predicted method and analytic method, as B.Dzezi l professors exist
Fouier two-dimensional integrations are introduced in Bduryky-noteh theoretical foundations and convert the method to be formed.The major defect of current method
It is：These researchs are all built upon on specific experiment basis, and research method relies on mainly conventional with research conclusion
Experience, and the emphasis of these methods is generally the evaluation in terms of geological disaster caused by the deformation of mine worked-out section,
It has ignored to mine worked-out section from mining phase to the observation for the structural strain's process, changing rule for adopting rear vacant deformation stage.
Nowadays the method for not occurring also that the stability of mine worked-out section can be assessed, it is impossible to fully and effectively right
The deformation of mine worked-out section is detected, while the stability of mine worked-out section is evaluated, and it is mined out to embody mine
The changing rule and feature in area.As shown in the above, using traditional deformation monitoring method, lacked due to itself can not be evaded
Fall into so that the monitoring of the deflection of mine worked-out section has deviateed with Practical Project, so as to cause the deflection in minery and stably
Property assess in error, it is impossible to meet minery engineering accuracy demand.
The content of the invention
In view of the above-mentioned deficiencies in the prior art, the technical problem to be solved by the present invention is that providing one kind is based on gravity
Abnormal mine worked-out section deflection and stability assessment method, its method and step is simple, reasonable in design and cost is low, uses behaviour
To make easy, the deflection in mining area can be monitored exactly using gravity anomaly, monitoring result is accurate, and according to adopting
The deflection in mining area is assessed the stability in mining area, and obtains the steady-state deformation amount in mining area, practicality
By force.
In order to solve the above technical problems, the technical solution adopted by the present invention is：A kind of mine based on gravity anomaly is mined out
Area's deflection and stability assessment method, it is characterised in that this method comprises the following steps：
Step 1: the arrangement of absolute gravimeter observation station and the erection of absolute gravimeter：First in rectangle mine worked-out section
Study in place, multiple absolute gravimeter observation stations are arranged at intervals along the length direction of mine worked-out section ground surface；Wherein, it is multiple
The absolute gravimeter observation station is denoted as absolute gravity observation point A respectively1, absolute gravity observation point A2..., absolute gravity observation
Point An, n is the quantity of absolute gravity observation point, and the mine worked-out section ground surface is rectangle mine worked-out section plane ABCD；
Then, respectively in the absolute gravity observation point A1, absolute gravity observation point A2..., absolute gravity observation point An
Place sets up absolute gravimeter；
Step 2: the acquisition of absolute gravity observation value：Multiple absolute gravimeters are according to the sampling time set in advance
The absolute gravity observation value of multiple absolute gravity observation points in step 1 is detected, and it is multiple absolute by what is detected
Observed gravimetric data transmits to data processor, data processor and obtains the exhausted of each multiple absolute gravimeter observation stations of sampling instant
To observed gravimetric data；
Step 3: establish rectangular coordinate system in space：With the plane of rectangle mine worked-out section described in step 1 ABCD summit A
For origin of coordinates O, the straight line for crossing rectangle mine worked-out section plane ABCD long side AB is Y-axis, crosses rectangle mine worked-out section plane
ABCD broadside AD straight line is X-axis, crosses rectangle mine worked-out section plane ABCD summit A and perpendicular to rectangle mine worked-out section
The straight line of plane is Z axis where plane ABCD；Wherein, the summit C of the rectangle mine worked-out section plane ABCD be located at by X-axis and
In the first quartile for the plane right-angle coordinate XOY that Y-axis is formed；
Step 4: the center-of-mass coordinate of mine worked-out section country rock spheroid calculates：The data processor is first according to the sampling time
Multiple absolute gravity observation values that order is gathered to each sampling instant obtained in step 2 afterwards are respectively processed, and are obtained
The center-of-mass coordinate of mine worked-out section country rock spheroid, and the processing of the multiple absolute gravity observation values gathered to each sampling instant
Method all same, wherein, the mine worked-out section country rock spheroid is the least radius ball of country rock body residing for rectangle mine worked-out section
Body, the multiple absolute gravity observation values gathered to the sampling instant of any of step 2 are handled, and process is as follows：
The Mass Calculation of step 401, mine worked-out section country rock spheroid：Using the data processor according to formulaThe mass M of mine worked-out section country rock spheroid is obtained, wherein, G is universal gravitational constant, and Δ g is
Absolute gravity observation value；
Step 402, mine worked-out section country rock spheroid barycenter setting coordinate：Using what is connected with the data processor
Parameter input module is by the coordinate P (x of the barycenter of mine worked-out section country rock spheroidi,yi,zi) abscissa xi, ordinate yiWith it is perpendicular
Coordinate ziInput in the data processor；Wherein, abscissa xi>0, ordinate yi>0, ordinate zi>0；
Step 403, mine worked-out section country rock spheroid barycenter mapping：Using the data processor according to formulaWithBy the matter of mine worked-out section country rock spheroid
Coordinate P (the x of the hearti,yi,zi) be mapped in the first quartile for the plane right-angle coordinate XOY being made up of X-axis and Y-axis, thrown
Shadow point P ' (xi′,yi'), then obtain mine worked-out section country rock spheroid barycenter coordinate P (xi,yi,zi) abscissa xiWith vertical seat
Mark yi；Wherein, xi=xi', yi=yi′；
Step 404, mine worked-out section country rock spheroid gravitation line drafting：Using the data processor in step 2
The multiple absolute gravity observation values obtained are respectively processed, and obtain a plurality of gravitation line of mine worked-out section country rock spheroid, then more
The intersection point of gravitation line described in bar is the barycenter of mine worked-out section country rock spheroid, and obtain the barycenter of mine worked-out section country rock spheroid
Coordinate, wherein, handled to obtain multiple gravitation lines of mine worked-out section country rock spheroid for multiple absolute gravity observation values
Method all same, then any one absolute gravity observation value is handled to obtain the gravitation line of mine worked-out section country rock spheroid, had
Body comprises the following steps：
Step 4041, using the data processor and according to formulaObtain normal gravity corresponding to absolute gravity observation point
Value γ0, wherein,Represent the latitude residing for absolute gravity observation point；
Step 4042, using the data processor and according to formulaAnd lead to
Cross normal gravity γ corresponding to absolute gravity observation point described in the parameter input module input step 40410, obtain ore deposit
Angle theta between projection and Z axis of the gravitation line of mountain goaf country rock spheroid on plane right-angle coordinate XOZxoz,z, wherein, r
Represent any point Q (x on the surface of mine worked-out section country rock spheroid*, y*, z*) arrive absolute gravity observation point space length, A
Represent any point Q (x on the surface of mine worked-out section country rock spheroid*, y*, z*) relative to the azimuth of absolute gravity observation point；
Step 4043, using the data processor and according to formulaAnd lead to
Cross normal gravity γ corresponding to absolute gravity observation point described in the parameter input module input step 40410, obtain ore deposit
Angle theta between projection and Z axis of the gravitation line of mountain goaf country rock spheroid on plane right-angle coordinate YOZyoz,z；
Step 4044, the drawing the mine worked-out section country rock spheroid obtained in step 4042 using the parameter input module
Angle theta between projection and Z axis of the line of force on plane right-angle coordinate XOZxoz,zIt is mined out with the mine that is obtained in step 4043
Angle theta between projection and Z axis of the gravitation line of area's country rock spheroid on plane right-angle coordinate YOZyoz,zInput the data
In processor, then gravitation line drafting module called using data processor, draw out the gravitation line of mine worked-out section country rock spheroid；
Step 4045, step 4041 is repeated several times to step 4045, obtains a plurality of gravitation of mine worked-out section country rock spheroid
Line, then the intersection of a plurality of gravitation line is the barycenter of mine worked-out section country rock spheroid；
Step 405, determine mine worked-out section country rock spheroid barycenter coordinate P (xi,yi,zi) ordinate：Using described
Data processor then obtains mine worked-out section country rock ball according to the barycenter of the mine worked-out section country rock spheroid obtained in step 4045
Coordinate P (the x of the barycenter of bodyi,yi,zi) ordinate zi；
Step 406, step 401 is repeated several times to step 405, completes the multiple definitely weights gathered to each sampling time
Power observation is handled, and obtains the coordinate of the barycenter of the mine worked-out section country rock spheroid in each sampling time, and described in use
Data processor is arranged according to time order and function order the coordinate of the barycenter of mine worked-out section country rock spheroid, and is denoted as respectively
P(x1,y1,z1)、P(x2,y2,z2)、P(x3,y3,z3)、...、P(xn,yn,zn)；Wherein, n represents absolute gravimeter sampled point
Sum；
Step 5: the drafting of the barycenter change curve of mine worked-out section country rock spheroid and the acquisition of steady-state deformation amount：Using
The data processor calls Drawing of Curve module, draws the barycenter change curve of mine worked-out section country rock spheroid, and according to institute
The barycenter change curve for stating mine worked-out section country rock spheroid obtains the deflection in mining area.
The above-mentioned mine worked-out section deflection and stability assessment method based on gravity anomaly, it is characterised in that：Step
The barycenter change curve of mine worked-out section country rock spheroid is drawn in five, and is become according to the barycenter of the mine worked-out section country rock spheroid
Change the deflection in curve acquisition mining area, comprise the following steps：
Step 501, using the data processor call Drawing of Curve module, draw out mine worked-out section country rock ball respectively
The abscissa x of the barycenter of bodyiOver time t change curve, mine worked-out section country rock spheroid barycenter ordinate yiWith when
Between t change curve and mine worked-out section country rock spheroid barycenter ordinate ziThe curve of t changes over time, obtains mine
The ordinate change curve of the abscissa change curve of the barycenter of goaf country rock spheroid, the barycenter of mine worked-out section country rock spheroid
With the ordinate change curve of the barycenter of mine worked-out section country rock spheroid；
Step 502, using the data processor respectively to the matter of the mine worked-out section country rock spheroid obtained in step 501
The ordinate change curve and mine worked-out section country rock ball of the abscissa change curve of the heart, the barycenter of mine worked-out section country rock spheroid
The ordinate change curve of the barycenter of body is judged, when the abscissa rate of change of two neighboring sampling instant is less than or equal to horizontal seat
During mark change setting value, illustrate that the variable quantity of the abscissa of the barycenter of mine worked-out section country rock spheroid has been stablized, and by adjacent two
The abscissa value of later moment in time is recorded as mine worked-out section country rock spheroid after the steady-state deformation of mine worked-out section in individual sampling instant
The abscissa of barycenter；
When the ordinate rate of change of two neighboring sampling instant is less than or equal to ordinate change setting value, illustrate that mine is adopted
The variable quantity of the ordinate of the barycenter of dead zone country rock spheroid has been stablized, and by the vertical seat of later moment in time in two neighboring sampling instant
Scale value is recorded as the ordinate of the barycenter of mine worked-out section country rock spheroid after the steady-state deformation of mine worked-out section；
When the ordinate rate of change of two neighboring sampling instant is less than or equal to ordinate change setting value, illustrate that mine is mined out
The variable quantity of the ordinate of the barycenter of area's country rock spheroid has been stablized, and by the ordinate of later moment in time in two neighboring sampling instant
Value is recorded as the ordinate of the barycenter of mine worked-out section country rock spheroid after the steady-state deformation of mine worked-out section；
It is step 503, respectively that mine after mine worked-out section steady-state deformation in step 502 is mined out using the data processor
The abscissa x of the abscissa of the barycenter of area's country rock spheroid and the barycenter of mine worked-out section country rock spheroid1, mine worked-out section is stable becomes
The ordinate y of the ordinate of the barycenter of mine worked-out section country rock spheroid and the barycenter of mine worked-out section country rock spheroid after shape1And ore deposit
The barycenter of the ordinate of the barycenter of mine worked-out section country rock spheroid and mine worked-out section country rock spheroid after the steady-state deformation of mountain goaf
Ordinate z1Mathematic interpolation is carried out, obtains stabilization of the barycenter in X-direction, Y-direction and Z-direction of mine worked-out section country rock spheroid
Deflection.
The above-mentioned mine worked-out section deflection and stability assessment method based on gravity anomaly, it is characterised in that：Step
The span in sampling time set in advance described in two is 12h~24h.
The above-mentioned mine worked-out section deflection and stability assessment method based on gravity anomaly, it is characterised in that：It is described
Data processor is computer.
The above-mentioned mine worked-out section deflection and stability assessment method based on gravity anomaly, it is characterised in that：Step
The level interval of two neighboring absolute gravimeter observation station is 0.2km~1km in one.
The above-mentioned mine worked-out section deflection and stability assessment method based on gravity anomaly, it is characterised in that：Step
The value model of the change of abscissa described in 502 setting value, ordinate change setting value and ordinate change setting value
Enclose is 1%~5%.
The present invention has advantages below compared with prior art：
1st, the inventive method step is simple and realizes conveniently, and monitoring efficiency and the degree of accuracy are high, and using effect is good.
2nd, deformation quantity measuring method of the invention is easy, fast and accurately, passes through the length in mine worked-out section ground surface
Direction is arranged at intervals multiple absolute gravimeter observation stations, and the absolute gravity observation value of the multiple gravimetric observations in mining area point is entered
Row detection, by handling multiple absolute gravity observation values, obtains the center-of-mass coordinate of mine worked-out section country rock spheroid, passes through
Monitor the change of the center-of-mass coordinate of mine worked-out section country rock spheroid and obtain the detection of the deflection in mining area, so as to realize
The continuous monitoring of mine worked-out section malformation amount.
3rd, when multiple absolute gravity observation values that the present invention is gathered to any one sampling instant are handled, use first
Data processor utilizes absolute gravity observation value, the quality of mine worked-out section country rock spheroid is obtained, then using data processor
The abscissa and ordinate of mine worked-out section country rock spheroid are obtained, mine worked-out section country rock ball is then obtained using data processor
The a plurality of gravitation line of body, then the intersection point of a plurality of gravitation line is the barycenter of mine worked-out section country rock spheroid；Finally use data processing
Device obtains the ordinate of the barycenter of mine worked-out section country rock spheroid, passes through the multiple absolute gravities gathered to each sampling time
Observation is handled, and is obtained the center-of-mass coordinate of the mine worked-out section country rock spheroid in each sampling time, is obtained mine worked-out section
The barycenter change curve of country rock spheroid, the detection of the mine worked-out section deflection based on gravity anomaly is realized, and according to mine
Goaf deflection carries out the assessment of stability, and the monitoring of mine worked-out section deflection and stability assessment are in mine worked-out section
Research is carried out in place, so can be detected in real time to the deflection in mining area exactly, monitoring result is accurate.
4th, the present invention is mined out on pertaining only to mine to the appraisal procedure of mine worked-out section TERM DEFORMATION state and stability
The absolute gravity observation of area's ground surface and the center-of-mass coordinate of mine worked-out section country rock body, pass through the matter of mine worked-out section country rock spheroid
Heart coordinate changes with time, fluctuation tendency is exactly mine worked-out section deforms the integrative feedback of process with the time, enters
And eliminate a series of physical during mining, geometry influence factor on goaf may caused by influence, Jin Erji
The earth simplifies former problematic complexity.
5th, the present invention is seen by absolute gravimeter according to the sampling time set in advance to the absolute gravity of gravity observation station
The long-term Continuous Observation of measured value, so as to obtain the barycenter change curve of mine worked-out section country rock spheroid, and enclosed according to mine worked-out section
The barycenter change curve of orbicule body judges the steady-state deformation amount of mine worked-out section, realize to the deformation of mine worked-out section and
The dynamic evaluation of stability, and the stability in mining area is assessed according to the deflection in minery, and obtain mine
The steady-state deformation amount in minery, it is practical.
In summary, the inventive method step is simple, reasonable in design and cost is low, and use is easy to operate, different using gravity
Often the deflection in mining area can be monitored exactly, monitoring result is accurate, and according to the deflection in minery to ore deposit
The stability in mountain minery is assessed, and obtains the steady-state deformation amount in mining area, practical.
Below by drawings and examples, technical scheme is described in further detail.
Brief description of the drawings
Fig. 1 is the method flow block diagram of the present invention.
Fig. 2 is the schematic block circuit diagram of the present invention.
Fig. 3 is the rectangular coordinate system in space schematic diagram that the present invention establishes.
Fig. 4 is the barycenter schematic diagram of country rock spheroid in mine worked-out section of the present invention.
Description of reference numerals:
1-gravitation line；2-rectangle mine worked-out section；
3-mine worked-out section country rock spheroid；The barycenter of 4-mine worked-out section country rock spheroid；
5-absolute gravimeter, 6-data processor；
7-parameter input module.
Embodiment
As shown in Figure 1, Figure 2, Figure 3 and Figure 4, a kind of mine worked-out section deflection based on gravity anomaly of the present invention and stably
Property appraisal procedure, this method comprise the following steps：
Step 1: the arrangement of absolute gravimeter observation station and the erection of absolute gravimeter：First in rectangle mine worked-out section 2
Study in place, multiple absolute gravimeter observation stations are arranged at intervals along the length direction of mine worked-out section ground surface；Wherein, it is multiple
The absolute gravimeter observation station is denoted as absolute gravity observation point A respectively1, absolute gravity observation point A2..., absolute gravity observation
Point An, n is the quantity of absolute gravity observation point, and the mine worked-out section ground surface is rectangle mine worked-out section plane ABCD；
Then, respectively in the absolute gravity observation point A1, absolute gravity observation point A2..., absolute gravity observation point An
Place sets up absolute gravimeter 5；
Step 2: the acquisition of absolute gravity observation value：Multiple absolute gravimeters 5 are according to the sampling time set in advance
The absolute gravity observation value of multiple absolute gravity observation points in step 1 is detected, and it is multiple absolute by what is detected
Observed gravimetric data is transmitted to data processor 6, and data processor 6 obtains the multiple absolute gravimeter observation stations of each sampling instant
Absolute gravity observation value；
Step 3: establish rectangular coordinate system in space：With the plane of rectangle mine worked-out section described in step 1 ABCD summit A
For origin of coordinates O, the straight line for crossing rectangle mine worked-out section plane ABCD long side AB is Y-axis, crosses rectangle mine worked-out section plane
ABCD broadside AD straight line is X-axis, crosses rectangle mine worked-out section plane ABCD summit A and perpendicular to rectangle mine worked-out section
The straight line of plane is Z axis where plane ABCD；Wherein, the summit C of the rectangle mine worked-out section plane ABCD be located at by X-axis and
In the first quartile for the plane right-angle coordinate XOY that Y-axis is formed；
Step 4: the center-of-mass coordinate of mine worked-out section country rock spheroid calculates：The data processor 6 is first according to the sampling time
Multiple absolute gravity observation values that order is gathered to each sampling instant obtained in step 2 afterwards are respectively processed, and are obtained
The coordinate of barycenter 4 of mine worked-out section country rock spheroid, and the place of the multiple absolute gravity observation values gathered to each sampling instant
Reason method all same, wherein, the mine worked-out section country rock spheroid 3 is the most smaller part of country rock body residing for rectangle mine worked-out section 2
Footpath spheroid, the multiple absolute gravity observation values gathered to the sampling instant of any of step 2 are handled, and process is as follows：
The Mass Calculation of step 401, mine worked-out section country rock spheroid：Using the data processor 6 according to formulaThe mass M of mine worked-out section country rock spheroid 3 is obtained, wherein, G is universal gravitational constant, Δ g
For absolute gravity observation value；
Step 402, mine worked-out section country rock spheroid barycenter setting coordinate：Connect using with the data processor 6
Parameter input module 7 by the coordinate P (x of the barycenter 4 of mine worked-out section country rock spheroidi,yi,zi) abscissa xi, ordinate yi
With ordinate ziInput in the data processor 6；Wherein, abscissa xi>0, ordinate yi>0, ordinate zi>0；
Step 403, mine worked-out section country rock spheroid barycenter mapping：Using the data processor 6 according to formulaWithBy the matter of mine worked-out section country rock spheroid
Coordinate P (the x of the heart 4i,yi,zi) be mapped in the first quartile for the plane right-angle coordinate XOY being made up of X-axis and Y-axis, thrown
Shadow point P ' (xi′,yi'), then obtain mine worked-out section country rock spheroid barycenter 4 coordinate P (xi,yi,zi) abscissa xiWith it is vertical
Coordinate yi；Wherein, xi=xi', yi=yi′；
Step 404, mine worked-out section country rock spheroid gravitation line drafting：Using the data processor 6 to step 2
Multiple absolute gravity observation values of middle acquisition are respectively processed, and obtain a plurality of gravitation line 1 of mine worked-out section country rock spheroid 3,
Then the intersection point of a plurality of gravitation line 1 is the barycenter 4 of mine worked-out section country rock spheroid, and obtains mine worked-out section country rock spheroid
The coordinate of barycenter 4, wherein, handled to obtain the multiple of mine worked-out section country rock spheroid 3 for multiple absolute gravity observation values
The method all same of gravitation line 1, then any one absolute gravity observation value is handled to obtain mine worked-out section country rock spheroid 3
Gravitation line 1, specifically includes following steps：
Step 4041, using the data processor 6 and according to formulaObtain normal gravity corresponding to absolute gravity observation point
Value γ0, wherein,Represent the latitude residing for absolute gravity observation point；
Step 4042, using the data processor 6 and according to formulaAnd lead to
Cross normal gravity γ corresponding to absolute gravity observation point described in the input step 4041 of parameter input module 70, obtain ore deposit
Angle theta between projection and Z axis of the gravitation line 1 of mountain goaf country rock spheroid 3 on plane right-angle coordinate XOZxoz,z, its
In, r represents any point Q (x on the surface of mine worked-out section country rock spheroid 3*, y*, z*) to absolute gravity observation point space away from
From A represents any point Q (x on the surface of mine worked-out section country rock spheroid 3*, y*, z*) relative to the side of absolute gravity observation point
Parallactic angle；
Step 4043, using the data processor 6 and according to formulaAnd lead to
Cross normal gravity γ corresponding to absolute gravity observation point described in the input step 4041 of parameter input module 70, obtain ore deposit
Angle theta between projection and Z axis of the gravitation line 1 of mountain goaf country rock spheroid 3 on plane right-angle coordinate YOZyoz,z；
Step 4044, using the parameter input module 7 by the mine worked-out section country rock spheroid 3 obtained in step 4042
Angle theta between projection and Z axis of the gravitation line 1 on plane right-angle coordinate XOZxoz,zAdopted with the mine obtained in step 4043
Angle theta between projection and Z axis of the gravitation line 1 of dead zone country rock spheroid 3 on plane right-angle coordinate YOZyoz,zDescribed in input
In data processor 6, then using the calling gravitation line drafting module of data processor 6, draw out mine worked-out section country rock spheroid 3
Gravitation line 1；
Step 4045, step 4041 is repeated several times to step 4045, obtains a plurality of gravitation of mine worked-out section country rock spheroid 3
Line 1, then the intersection of a plurality of gravitation line 1 is the barycenter 4 of mine worked-out section country rock spheroid；
Step 405, determine mine worked-out section country rock spheroid barycenter 4 coordinate P (xi,yi,zi) ordinate：Using institute
Barycenter 4 of the data processor 6 according to the mine worked-out section country rock spheroid obtained in step 4045 is stated, then must be enclosed to mine worked-out section
Coordinate P (the x of the barycenter 4 of orbicule bodyi,yi,zi) ordinate zi；
Step 406, step 401 is repeated several times to step 405, completes the multiple definitely weights gathered to each sampling time
Power observation is handled, and obtains the coordinate of the barycenter 4 of the mine worked-out section country rock spheroid in each sampling time, and described in use
Data processor 6 is arranged according to time order and function order the coordinate of the barycenter 4 of mine worked-out section country rock spheroid, and is remembered respectively
Make P (x1,y1,z1)、P(x2,y2,z2)、P(x3,y3,z3)、...、P(xn,yn,zn)；Wherein, n represents the sampled point of absolute gravimeter 5
Sum；
Step 5: the drafting of the barycenter change curve of mine worked-out section country rock spheroid and the acquisition of steady-state deformation amount：Using
The data processor calls Drawing of Curve module, draws the barycenter change curve of mine worked-out section country rock spheroid, and according to institute
The barycenter change curve for stating mine worked-out section country rock spheroid obtains the deflection in mining area.
In the present embodiment, the barycenter change curve of mine worked-out section country rock spheroid is drawn in step 5, and according to the ore deposit
The barycenter change curve of mountain goaf country rock spheroid obtains the deflection in mining area, comprises the following steps：
Step 501, using the data processor 6 call Drawing of Curve module, draw out mine worked-out section country rock respectively
The abscissa x of the barycenter 4 of spheroidiOver time t change curve, mine worked-out section country rock spheroid barycenter 4 ordinate yiWith
The ordinate z of the curve of time t changes and the barycenter 4 of mine worked-out section country rock spheroidiThe curve of t changes over time, is obtained
Take the abscissa change curve of the barycenter 4 of mine worked-out section country rock spheroid, mine worked-out section country rock spheroid barycenter 4 ordinate
The ordinate change curve of the barycenter 4 of change curve and mine worked-out section country rock spheroid；
Step 502, using the data processor 6 respectively to the mine worked-out section country rock spheroid that is obtained in step 501
The abscissa change curve of barycenter 4, the ordinate change curve of the barycenter 4 of mine worked-out section country rock spheroid and mine worked-out section are enclosed
The ordinate change curve of the barycenter 4 of orbicule body judged, when the abscissa rate of change of two neighboring sampling instant is less than etc.
When abscissa changes setting value, illustrate that the variable quantity of the abscissa of the barycenter 4 of mine worked-out section country rock spheroid has been stablized, and will
The abscissa value of later moment in time is recorded as mine worked-out section country rock after the steady-state deformation of mine worked-out section in two neighboring sampling instant
The abscissa of the barycenter 4 of spheroid；
When the ordinate rate of change of two neighboring sampling instant is less than or equal to ordinate change setting value, illustrate that mine is adopted
The variable quantity of the ordinate of the barycenter 4 of dead zone country rock spheroid has been stablized, and by two neighboring sampling instant later moment in time it is vertical
Coordinate value is recorded as the ordinate of the barycenter 4 of mine worked-out section country rock spheroid after the steady-state deformation of mine worked-out section；
When the ordinate rate of change of two neighboring sampling instant is less than or equal to ordinate change setting value, illustrate that mine is mined out
The variable quantity of the ordinate of the barycenter 4 of area's country rock spheroid has been stablized, and by the perpendicular seat of later moment in time in two neighboring sampling instant
Scale value is recorded as the ordinate of the barycenter 4 of mine worked-out section country rock spheroid after the steady-state deformation of mine worked-out section；
Step 503, using the data processor 6 mine after mine worked-out section steady-state deformation in step 502 is adopted respectively
The abscissa x of the abscissa of the barycenter 4 of dead zone country rock spheroid and the barycenter 4 of mine worked-out section country rock spheroid1, mine worked-out section it is steady
The ordinate of the ordinate and the barycenter 4 of mine worked-out section country rock spheroid of the barycenter 4 of mine worked-out section country rock spheroid after fixed deformation
y1With the ordinate of the barycenter 4 of mine worked-out section country rock spheroid after the steady-state deformation of mine worked-out section and mine worked-out section country rock spheroid
Barycenter 4 ordinate z1Mathematic interpolation is carried out, obtains the barycenter 4 of mine worked-out section country rock spheroid in X-direction, Y-direction and Z side
To steady-state deformation amount.
In the present embodiment, the span in sampling time set in advance described in step 2 is 12h~24h.
In the present embodiment, the data processor 6 is computer.
In the present embodiment, the level interval of two neighboring absolute gravimeter observation station is 0.2km~1km in step 1.
In the present embodiment, abscissa described in step 502 changes setting value, the ordinate changes setting value and described perpendicular
The span of changes in coordinates setting value is 1%~5%.
In the present embodiment, considered according to detection demand and financial cost, three definitely weights are set in actual mine worked-out section
Power instrument observation station, absolute gravity observation point A is denoted as respectively1, absolute gravity observation point A2With absolute gravity observation point A3, actually make
During, the absolute gravimeter observation station of other quantity can be set according to mine worked-out section detection demand, meet detection demand
.
As shown in figure 4, in the present embodiment, three according to acquired in the three of setting absolute gravimeter observation stations are definitely heavy
Power observation is respectively processed, and obtains three gravitation lines 1 of mine worked-out section country rock spheroid 3, then three gravitation lines 1
Intersection point is the barycenter 4 of mine worked-out section country rock spheroid.
In the present embodiment, changed with time by the coordinate of the barycenter 4 of mine worked-out section country rock spheroid, fluctuation tendency just
Be mine worked-out section deforms the integrative feedback of process with the time, and then eliminates a series of during mining
Physics, the influence caused by goaf is possible of geometry influence factor, and then greatly simplifie former problematic complexity.
As shown in the above, the inventive method step is simple, reasonable in design and cost is low, and use is easy to operate, utilizes
Gravity anomaly can be monitored to the deflection in mining area exactly, and monitoring result is accurate, and according to the deformation in minery
Measure and the stability in mining area is assessed, and obtain the steady-state deformation amount in mining area, it is practical.
It is described above, only it is presently preferred embodiments of the present invention, not the present invention is imposed any restrictions, it is every according to the present invention
Any simple modification, change and the equivalent structure change that technical spirit is made to above example, still fall within skill of the present invention
In the protection domain of art scheme.
Claims (6)
1. a kind of mine worked-out section deflection and stability assessment method based on gravity anomaly, it is characterised in that this method bag
Include following steps：
Step 1: the arrangement of absolute gravimeter observation station and the erection of absolute gravimeter：Ground first in rectangle mine worked-out section (2)
Study carefully in place, multiple absolute gravimeter observation stations are arranged at intervals along the length direction of mine worked-out section ground surface；Wherein, Duo Gesuo
State absolute gravimeter observation station and be denoted as absolute gravity observation point A respectively1, absolute gravity observation point A2..., absolute gravity observation point
An, n is the quantity of absolute gravity observation point, and the mine worked-out section ground surface is rectangle mine worked-out section plane ABCD；
Then, respectively in the absolute gravity observation point A1, absolute gravity observation point A2..., absolute gravity observation point AnLocate frame
If absolute gravimeter (5)；
Step 2: the acquisition of absolute gravity observation value：Multiple absolute gravimeters (5) are according to the sampling time pair set in advance
The absolute gravity observation value of multiple absolute gravity observation points is detected in step 1, and the multiple definitely weights that will be detected
Power observation is transmitted to data processor (6), data processor (6) and obtains each multiple absolute gravimeter observation stations of sampling instant
Absolute gravity observation value；
Step 3: establish rectangular coordinate system in space：Using the plane of rectangle mine worked-out section described in step 1 ABCD summit A as seat
Origin O is marked, the straight line for crossing rectangle mine worked-out section plane ABCD long side AB is Y-axis, crosses rectangle mine worked-out section plane ABCD
Broadside AD straight line be X-axis, cross rectangle mine worked-out section plane ABCD summit A and perpendicular to rectangle mine worked-out section plane
The straight line of plane is Z axis where ABCD；Wherein, the summit C of the rectangle mine worked-out section plane ABCD is located at by X-axis and Y-axis
In the plane right-angle coordinate XOY of composition first quartile；
Step 4: the center-of-mass coordinate of mine worked-out section country rock spheroid calculates：The data processor (6) is according to sampling time priority
Multiple absolute gravity observation values that order is gathered to each sampling instant obtained in step 2 are respectively processed, and obtain ore deposit
Barycenter (4) coordinate of mountain goaf country rock spheroid, and the place of the multiple absolute gravity observation values gathered to each sampling instant
Reason method all same, wherein, the mine worked-out section country rock spheroid (3) be country rock body residing for rectangle mine worked-out section (2) most
Minor radius spheroid, the multiple absolute gravity observation values gathered to the sampling instant of any of step 2 are handled, and process is such as
Under：
The Mass Calculation of step 401, mine worked-out section country rock spheroid：Using the data processor (6) according to formulaThe mass M of mine worked-out section country rock spheroid (3) is obtained, wherein, G is universal gravitational constant, Δ
G is absolute gravity observation value；
Step 402, mine worked-out section country rock spheroid barycenter setting coordinate：Using what is connected with the data processor (6)
Parameter input module (7) is by the coordinate P (x of the barycenter (4) of mine worked-out section country rock spheroidi,yi,zi) abscissa xi, ordinate
yiWith ordinate ziInput in the data processor (6)；Wherein, abscissa xi>0, ordinate yi>0, ordinate zi>0；
Step 403, mine worked-out section country rock spheroid barycenter mapping：Using the data processor (6) according to formulaWithBy the matter of mine worked-out section country rock spheroid
Coordinate P (the x of the heart (4)i,yi,zi) be mapped in the first quartile for the plane right-angle coordinate XOY being made up of X-axis and Y-axis, obtain
Subpoint P ' (xi′,yi'), then obtain mine worked-out section country rock spheroid barycenter (4) coordinate P (xi,yi,zi) abscissa xi
With ordinate yi；Wherein, xi=x 'i, yi=y 'i；
Step 404, mine worked-out section country rock spheroid gravitation line drafting：Using the data processor (6) in step 2
The multiple absolute gravity observation values obtained are respectively processed, and obtain a plurality of gravitation line of mine worked-out section country rock spheroid (3)
(1), then the intersection point of a plurality of gravitation line (1) is the barycenter (4) of mine worked-out section country rock spheroid, and obtains mine worked-out section and enclose
The coordinate of the barycenter (4) of orbicule body, wherein, handled to obtain mine worked-out section country rock ball for multiple absolute gravity observation values
The method all same of multiple gravitation lines (1) of body (3), then any one absolute gravity observation value is handled to obtain mine mined out
The gravitation line (1) of area's country rock spheroid (3), specifically includes following steps：
Step 4041, using the data processor (6) and according to formulaObtain normal gravity corresponding to absolute gravity observation point
Value γ0, wherein,Represent the latitude residing for absolute gravity observation point；
Step 4042, using the data processor (6) and according to formulaAnd pass through
Normal gravity γ corresponding to absolute gravity observation point described in parameter input module (7) input step 40410, obtain ore deposit
Angle theta between projection and Z axis of the gravitation line (1) of mountain goaf country rock spheroid (3) on plane right-angle coordinate XOZxoz,z,
Wherein, r represents any point Q (x on the surface of mine worked-out section country rock spheroid (3)*, y*, z*) arrive absolute gravity observation point sky
Between distance, A represent mine worked-out section country rock spheroid (3) surface on any point Q (x*, y*, z*) relative to absolute gravity observation
The azimuth of point；
Step 4043, using the data processor (6) and according to formulaAnd pass through
Normal gravity γ corresponding to absolute gravity observation point described in parameter input module (7) input step 40410, obtain ore deposit
Angle theta between projection and Z axis of the gravitation line (1) of mountain goaf country rock spheroid (3) on plane right-angle coordinate YOZyoz,z；
Step 4044, using the parameter input module (7) by the mine worked-out section country rock spheroid (3) obtained in step 4042
Angle theta between projection and Z axis of the gravitation line (1) on plane right-angle coordinate XOZxoz,zWith the mine obtained in step 4043
Angle theta between projection and Z axis of the gravitation line (1) of goaf country rock spheroid (3) on plane right-angle coordinate YOZyoz,zIt is defeated
Enter in the data processor (6), then gravitation line drafting module is called using data processor (6), draw out mine worked-out section
The gravitation line (1) of country rock spheroid (3)；
Step 4045, step 4041 is repeated several times to step 4045, obtains a plurality of gravitation line of mine worked-out section country rock spheroid (3)
(1), then the intersection of a plurality of gravitation line (1) is the barycenter (4) of mine worked-out section country rock spheroid；
Step 405, determine mine worked-out section country rock spheroid barycenter (4) coordinate P (xi,yi,zi) ordinate：Using described
Data processor (6) then obtains mine worked-out section according to the barycenter (4) of the mine worked-out section country rock spheroid obtained in step 4045
Coordinate P (the x of the barycenter (4) of country rock spheroidi,yi,zi) ordinate zi；
Step 406, step 401 is repeated several times to step 405, the multiple absolute gravities for completing to gather each sampling time are seen
Measured value is handled, and obtains the coordinate of the barycenter (4) of the mine worked-out section country rock spheroid in each sampling time, and use the number
The coordinate of the barycenter (4) of mine worked-out section country rock spheroid is arranged according to time order and function order according to processor (6), and respectively
It is denoted as P (x1,y1,z1)、P(x2,y2,z2)、P(x3,y3,z3)、...、P(xn,yn,zn)；Wherein, n represents that absolute gravimeter (5) is adopted
The sum of sampling point；
Step 5: the drafting of the barycenter change curve of mine worked-out section country rock spheroid and the acquisition of steady-state deformation amount：Using described
Data processor calls Drawing of Curve module, draws the barycenter change curve of mine worked-out section country rock spheroid, and according to the ore deposit
The barycenter change curve of mountain goaf country rock spheroid obtains the deflection in mining area.
2. according to the mine worked-out section deflection and stability assessment method based on gravity anomaly described in claim 1, it is special
Sign is：The barycenter change curve of mine worked-out section country rock spheroid is drawn in step 5, and according to the mine worked-out section country rock
The barycenter change curve of spheroid obtains the deflection in mining area, comprises the following steps：
Step 501, using the data processor (6) call Drawing of Curve module, draw out mine worked-out section country rock ball respectively
The abscissa x of the barycenter (4) of bodyiOver time t change curve, mine worked-out section country rock spheroid barycenter (4) ordinate yi
The ordinate z of the barycenter (4) of the curve of t changes and mine worked-out section country rock spheroid over timeiThe song of t changes over time
Line, obtain abscissa change curve, the barycenter of mine worked-out section country rock spheroid of the barycenter (4) of mine worked-out section country rock spheroid
(4) ordinate change curve and the ordinate change curve of the barycenter (4) of mine worked-out section country rock spheroid；
Step 502, using the data processor (6) respectively to the matter of the mine worked-out section country rock spheroid obtained in step 501
The abscissa change curve of the heart (4), the ordinate change curve of the barycenter (4) of mine worked-out section country rock spheroid and mine worked-out section
The ordinate change curve of the barycenter (4) of country rock spheroid is judged, when the abscissa rate of change of two neighboring sampling instant is small
When setting value is changed equal to abscissa, illustrate that the variable quantity of the abscissa of the barycenter (4) of mine worked-out section country rock spheroid is steady
Calmly, mine is mined out after and the abscissa value of later moment in time in two neighboring sampling instant being recorded as into mine worked-out section steady-state deformation
The abscissa of the barycenter (4) of area's country rock spheroid；
When the ordinate rate of change of two neighboring sampling instant is less than or equal to ordinate change setting value, illustrate mine worked-out section
The variable quantity of the ordinate of the barycenter (4) of country rock spheroid has been stablized, and by the vertical seat of later moment in time in two neighboring sampling instant
Scale value is recorded as the ordinate of the barycenter (4) of mine worked-out section country rock spheroid after the steady-state deformation of mine worked-out section；
When the ordinate rate of change of two neighboring sampling instant is less than or equal to ordinate change setting value, illustrate that mine worked-out section is enclosed
The variable quantity of the ordinate of the barycenter (4) of orbicule body has been stablized, and by the ordinate of later moment in time in two neighboring sampling instant
Value is recorded as the ordinate of the barycenter (4) of mine worked-out section country rock spheroid after the steady-state deformation of mine worked-out section；
It is step 503, respectively that mine after mine worked-out section steady-state deformation in step 502 is mined out using the data processor (6)
The abscissa x of the abscissa of the barycenter (4) of area's country rock spheroid and the barycenter (4) of mine worked-out section country rock spheroid1, mine worked-out section
The ordinate of the barycenter (4) of mine worked-out section country rock spheroid and the barycenter (4) of mine worked-out section country rock spheroid after steady-state deformation
Ordinate y1With the ordinate of the barycenter (4) of mine worked-out section country rock spheroid and mine worked-out section after the steady-state deformation of mine worked-out section
The ordinate z of the barycenter (4) of country rock spheroid1Mathematic interpolation is carried out, obtains the barycenter (4) of mine worked-out section country rock spheroid in X side
To, Y-direction and the steady-state deformation amount of Z-direction.
3. according to the mine worked-out section deflection and stability assessment method based on gravity anomaly described in claim 1 or 2, its
It is characterised by：The span in sampling time set in advance described in step 2 is 12h~24h.
4. according to the mine worked-out section deflection and stability assessment method based on gravity anomaly described in claim 1 or 2, its
It is characterised by：The data processor (6) is computer.
5. according to the mine worked-out section deflection and stability assessment method based on gravity anomaly described in claim 1 or 2, its
It is characterised by：The level interval of two neighboring absolute gravimeter observation station is 0.2km~1km in step 1.
6. according to the mine worked-out section deflection and stability assessment method based on gravity anomaly described in claim 2, it is special
Sign is：The change of abscissa described in step 502 setting value, ordinate change setting value and ordinate change setting
The span of value is 1%~5%.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710388355.6A CN107219560B (en) | 2017-05-27 | 2017-05-27 | Mine worked-out section deflection and stability assessment method based on gravity anomaly |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710388355.6A CN107219560B (en) | 2017-05-27 | 2017-05-27 | Mine worked-out section deflection and stability assessment method based on gravity anomaly |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107219560A CN107219560A (en) | 2017-09-29 |
| CN107219560B true CN107219560B (en) | 2018-03-16 |
Family
ID=59947962
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710388355.6A Expired - Fee Related CN107219560B (en) | 2017-05-27 | 2017-05-27 | Mine worked-out section deflection and stability assessment method based on gravity anomaly |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107219560B (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115199336B (en) * | 2022-07-15 | 2023-06-06 | 中钢集团马鞍山矿山研究总院股份有限公司 | Real-time monitoring system and modeling method for mine goaf morphology |
| CN116777075B (en) * | 2023-06-30 | 2024-03-29 | 江苏中煤地质工程研究院有限公司 | Goaf influence range determination method and goaf influence range determination system |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103189760A (en) * | 2010-09-03 | 2013-07-03 | Geosoft公司 | Method and system for modeling anomalous density zones in geophysical exploration |
| CN104166170A (en) * | 2014-08-13 | 2014-11-26 | 昆明理工大学 | Method for positioning and detecting low-density concealed ore bodies in gallery gravity total spatial domain |
| CN104536057A (en) * | 2014-12-30 | 2015-04-22 | 山东科技大学 | Method for monitoring relative gravity of surface subsidence in coal mining process |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| GB201008993D0 (en) * | 2010-05-28 | 2010-07-14 | Arkex Ltd | Processing geophysical data |
-
2017
- 2017-05-27 CN CN201710388355.6A patent/CN107219560B/en not_active Expired - Fee Related
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103189760A (en) * | 2010-09-03 | 2013-07-03 | Geosoft公司 | Method and system for modeling anomalous density zones in geophysical exploration |
| CN104166170A (en) * | 2014-08-13 | 2014-11-26 | 昆明理工大学 | Method for positioning and detecting low-density concealed ore bodies in gallery gravity total spatial domain |
| CN104536057A (en) * | 2014-12-30 | 2015-04-22 | 山东科技大学 | Method for monitoring relative gravity of surface subsidence in coal mining process |
Non-Patent Citations (1)
| Title |
|---|
| 微重力法在采空区勘查中的应用;王延涛 等;《物探与化探》;20121015;第36卷;第61-64页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107219560A (en) | 2017-09-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Yazdani et al. | Displacement-based numerical back analysis for estimation of rock mass parameters in Siah Bisheh powerhouse cavern using continuum and discontinuum approach | |
| CN103233777B (en) | A kind of roof safety ore pressure monitor method and device | |
| Ren et al. | Model experimental research on deformation and subsidence characteristics of ground and wall rock due to mining under thick overlying terrane | |
| CN104133996A (en) | Ground settlement risk grade evaluation method based on cloud model and data field | |
| Liu et al. | Mechanism of formation of sliding ground fissure in loess hilly areas caused by underground mining | |
| CN103791884A (en) | Research method of integration rule for coal mining overlying strata and ground surface movement deformation | |
| CN108693572A (en) | A kind of ground stress field inversion method based on three-dimensional modeling | |
| CN103606019B (en) | Mine goaf overlying stratum sedimentation dynamic prediction method based on time-space relationship | |
| CN107219560B (en) | Mine worked-out section deflection and stability assessment method based on gravity anomaly | |
| CN104794315A (en) | Coal seam floor damaged zone depth determination method based on improved rupture mechanical model | |
| CN103778480A (en) | Fissure-zone height prediction method based on sensitivity analysis | |
| Zhao et al. | Monitoring and analysis of the mining-induced ground movement in the Longshou Mine, China | |
| CN111815471A (en) | Coal mine resource mining area geological disaster monitoring and early warning management system based on big data | |
| Pei et al. | Optimized back analysis method for stress determination based on identification of local stress measurements and its application | |
| CN112488455A (en) | Slope stability analysis and evaluation method, system, medium, equipment and terminal | |
| CN114692456A (en) | Determination method for inducing ground stress critical criterion of mine strong dynamic disaster | |
| Li et al. | Prediction of mining subsidence in shallow coal seam | |
| CN104866682B (en) | Ground deformation-based method for inversion of tectonic stress field in shale gas exploration area | |
| CN107036575B (en) | Mine worked-out section deformation and Detection of Stability method based on bouguer gravity anomaly | |
| Shi et al. | Study on the mechanisms and prevention of water inrush events in a deeply buried high-pressure coal seam—a case study of the Chensilou coal mine in China | |
| Jia et al. | Stability analysis of shallow goaf based on field monitoring and numerical simulation: A case study at an Open-Pit Iron Mine, China | |
| Lu et al. | Monitoring and analysis of ground subsidence and backfill stress distribution in Jinchuan Mine, China | |
| Yuanhui et al. | Analysis of Mechanical Characteristics and Instability Law of Fault under the Influence of Mining | |
| Lu et al. | Floor Disturbance and Failure Characteristics of Super‐Large Mining Height Working Face | |
| CN114048582A (en) | Method and device for predicting mine water disaster spreading process, electronic equipment and storage medium |
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 | ||
| CB03 | Change of inventor or designer information |
Inventor after: Wu Xiaolong Inventor after: Gong Yun Inventor before: Wu Xiaolong Inventor before: Gong Yun Inventor before: Tang Fuquan |
|
| CB03 | Change of inventor or designer information | ||
| CB03 | Change of inventor or designer information |
Inventor after: Wu Xiaolong Inventor after: Tang Jingya Inventor before: Wu Xiaolong Inventor before: Gong Yun |
|
| CB03 | Change of inventor or designer information | ||
| TA01 | Transfer of patent application right |
Effective date of registration: 20180108 Address after: Yanta Road 710054 Shaanxi city of Xi'an province No. 58 Applicant after: Xi'an University of Science and Technology Applicant after: Changan Univ. Address before: Yanta Road 710054 Shaanxi city of Xi'an province No. 58 Applicant before: Xi'an University of Science and Technology |
|
| TA01 | Transfer of patent application right | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180316 Termination date: 20190527 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |