CN108104876B - Water damage real-time graded method for early warning and system based on mine electric method monitoring - Google Patents

Water damage real-time graded method for early warning and system based on mine electric method monitoring Download PDF

Info

Publication number
CN108104876B
CN108104876B CN201711265780.2A CN201711265780A CN108104876B CN 108104876 B CN108104876 B CN 108104876B CN 201711265780 A CN201711265780 A CN 201711265780A CN 108104876 B CN108104876 B CN 108104876B
Authority
CN
China
Prior art keywords
resistivity
early warning
local minimum
inverting
level
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201711265780.2A
Other languages
Chinese (zh)
Other versions
CN108104876A (en
Inventor
鲁晶津
宁殿艳
王冰纯
段建华
范涛
赵兆
岳辉
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xian Research Institute Co Ltd of CCTEG
Original Assignee
Xian Research Institute Co Ltd of CCTEG
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Xian Research Institute Co Ltd of CCTEG filed Critical Xian Research Institute Co Ltd of CCTEG
Priority to CN201711265780.2A priority Critical patent/CN108104876B/en
Publication of CN108104876A publication Critical patent/CN108104876A/en
Application granted granted Critical
Publication of CN108104876B publication Critical patent/CN108104876B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • EFIXED CONSTRUCTIONS
    • E21EARTH DRILLING; MINING
    • E21FSAFETY DEVICES, TRANSPORT, FILLING-UP, RESCUE, VENTILATION, OR DRAINING IN OR OF MINES OR TUNNELS
    • E21F17/00Methods or devices for use in mines or tunnels, not covered elsewhere
    • E21F17/18Special adaptations of signalling or alarm devices

Abstract

The present invention carries out automatic, processing in real time to monitoring data using 3 D resistivity inverting, according to the real-time change situation of coal measure strata resistivity obtained by inverting, automatically classification judgement is carried out to water gush face risk, automatically, grading forewarning system information is sent in real time, realize the automation and intelligence of working face water damage early warning, control the developmental state of stope conduit pipe in real time convenient for monitoring personnel, gushing water hidden danger is checked in time, provides sound assurance for the life security of miner and the property safety of mine side.

Description

Water damage real-time graded method for early warning and system based on mine electric method monitoring
Technical field
The present invention relates to a kind of water damage real time early warning method and systems, belong to pit mining technical field, are specifically related to A kind of water damage real-time graded method for early warning and system based on mine electric method monitoring.
Background technique
Working face water damage often by containing, water guide geological structure induce.Contain, water guide geological structure will not draw in a static condition Water inrush accident is sent out, but under the influence of mining work activities, since crustal stress distribution changes, the structural fracture of rock is caused to expand Or new rupture occurs, contain, water guide geological structure can further develop conduit pipe between coal seam and water-bearing layer, to lead Cause the generation of working face water damage.Contain, water guide geological structure is with respect to low-resistance is normally behaved as coal measure strata, mine electric can be used Method is detected.
Mine electric method detects identification abnormal area according to the electrical property difference between anomalous geologic body and country rock, with other physical prospectings Method is compared, and has unique advantage when detection shows the containing of low-resistance characteristic, permeable structure.Contain, the development of water guide geological structure During conduit pipe, the range of low resistance region can gradually change, and therefore, can answer during working face extraction It is monitored with growth course of the mine electric method to conduit pipe, provides reference frame for the prediction of analysis of possibility of water inrush.
However, mine electric method monitoring materials can only intuitively reflect the electrical situation of change of coal measure strata, how according to electricity Property situation of change judge that the developmental state of conduit pipe then needs the experience by data interpretation personnel to go to judge.Artificial experience is sentenced It is disconnected to have certain hysteresis quality in time, it, can not be right sometimes as electrically changing the faint generation for causing to fail to judge or judge by accident Water gush face risk carries out timely and effectively early warning.Currently, what is faced during working face extraction is mostly high artesian, water guide Once being formed, the generation of water damage is all extremely swift and violent in channel, has seriously threatened the life security of miner and the property safety of mine side, because This, the working face water damage based on mine electric method monitoring is automatic, in real time, the realization of grading forewarning system method very it is necessary to.
Summary of the invention
The present invention mainly solves water damage early warning present in the prior art, and there are hysteresis qualitys, and electrically variation is faint causes to leak The generation sentenced or judged by accident, can not to water gush face risk carry out timely and effective early warning the technical issues of, provide one kind and be based on The working face water damage real time early warning method and system of mine electric method monitoring.This method and system are according to the reality of coal measure strata resistivity When situation of change judged water gush face risk and carried out automatic, real-time, grading forewarning system, it is pre- to realize working face water damage Alert automation and intelligence, timeliness with higher also can effectively avoid artificial erroneous judgement.
Above-mentioned technical problem of the invention is mainly to be addressed by following technical proposals:
A kind of working face water damage real-time graded method for early warning based on mine electric method monitoring, comprising:
Initial modeling procedure establishes the initial 3 D resistivity mould of coal measure strata according to the geologic information of stope Type;
Inverting warning step is obtained using the initial 3 D resistivity model inversion each batch underground Real-time Monitoring Data Inversion of three dimensional resistivity models, when the inversion of three dimensional resistivity models of two neighboring batch are overlapped and in differential resistivity value sequence There are when negative value, being judged as that there are risk and early warning, otherwise continue inverting next batch underground Real-time Monitoring Data;
Wherein, the differential resistivity value sequence is resistivity value and the last time inverting three of this inversion of three dimensional resistivity models Tie up the difference of the resistivity value of resistivity models.
Optimization, a kind of above-mentioned working face water damage real-time graded method for early warning based on mine electric method monitoring, the number According to inverting warning step further include:
Three-level early warning sub-step finds out sequence of differences according to coordinate information for the point of small Yu 0 in differential resistivity value sequence Local minimum, save corresponding with its coordinate position of local minimum;Send includes the three of local minimum coordinate position Grade warning information.
Optimization, a kind of above-mentioned working face water damage real-time graded method for early warning based on mine electric method monitoring, the number According to inverting warning step further include:
Second level early warning sub-step judges the local minimum that this is calculated and the local minimum that last computation obtains With the presence or absence of position coincidence point, if then saving the coordinate position of coincidence, the second level early warning letter comprising being overlapped coordinate position is sent Breath.
Optimization, a kind of above-mentioned working face water damage real-time graded method for early warning based on mine electric method monitoring, the number According to inverting warning step further include:
Level-one early warning sub-step, for position coincidence point, primary calculate is obtained before being subtracted with the local minimum newly obtained Local minimum, calculate the difference of local minimum twice, judge whether there is point of the difference less than 0, and if it exists, then save The location information of the point, and send the level-one warning information comprising the coordinate position.
A kind of working face water damage real-time graded early warning system based on mine electric method monitoring, comprising:
Initial modeling module establishes the initial 3 D resistivity mould of coal measure strata according to the geologic information of stope Type;
Inverting warning module is obtained using the initial 3 D resistivity model inversion each batch underground Real-time Monitoring Data Inversion of three dimensional resistivity models, when the inversion of three dimensional resistivity models of two neighboring batch are overlapped and in differential resistivity value sequence There are when negative value, being judged as that there are risk and early warning, otherwise continue inverting next batch underground Real-time Monitoring Data;
Wherein, the differential resistivity value sequence is resistivity value and the last time inverting three of this inversion of three dimensional resistivity models Tie up the difference of the resistivity value of resistivity models.
Optimization, a kind of above-mentioned working face water damage real-time graded early warning system based on mine electric method monitoring, the number According to inverting warning module further include:
Three-level early warning submodule finds out sequence of differences according to coordinate information for the point of small Yu 0 in differential resistivity value sequence Local minimum, save corresponding with its coordinate position of local minimum;Send includes the three of local minimum coordinate position Grade warning information.
Optimization, a kind of above-mentioned working face water damage real-time graded early warning system based on mine electric method monitoring, the number According to inverting warning module further include:
Second level early warning submodule judges the local minimum that this is calculated and the local minimum that last computation obtains With the presence or absence of position coincidence point, if then saving the coordinate position of coincidence, the second level early warning letter comprising being overlapped coordinate position is sent Breath.
Optimization, a kind of above-mentioned working face water damage real-time graded early warning system based on mine electric method monitoring, the number According to inverting warning module further include:
Level-one early warning submodule, for position coincidence point, primary calculate is obtained before being subtracted with the local minimum newly obtained Local minimum, calculate the difference of local minimum twice, judge whether there is point of the difference less than 0, and if it exists, then save The location information of the point, and send the level-one warning information comprising the coordinate position.
Therefore, the present invention has the advantage that carrying out automatic, real-time place to monitoring data using 3 D resistivity inverting Reason carries out classification judgement to water gush face risk automatically according to the real-time change situation of coal measure strata resistivity obtained by inverting, Automatically, grading forewarning system information is sent in real time, realizes the automation and intelligence of working face water damage early warning, it is real convenient for monitoring personnel When control stope conduit pipe developmental state, gushing water hidden danger is checked in time, be miner life security and The property safety of mine side provides sound assurance.
Detailed description of the invention
Attached drawing 1 is a kind of flow chart of the invention.
Specific embodiment
Below with reference to the embodiments and with reference to the accompanying drawing the technical solutions of the present invention will be further described.
Embodiment:
The monitoring data of mine electric method monitoring system acquisition are generally current potential or potential difference, at 3 D resistivity inverting Reason can obtain the 3 D resistivity data that can more directly reflect coal measure strata Electrical distribution situation.Conduit pipe growth course In, the zone resistance rate being affected by it can gradually decrease.Pass through the real-time change for the coal measure strata resistivity that analysis inverting obtains Situation, can the developmental state to conduit pipe effectively judged, thus achieve the purpose that analysis of possibility of water inrush carry out early warning.
As shown in Figure 1, the embodiment of the present invention provide it is a kind of based on mine electric method monitoring working face water damage divide in real time automatically Grade method for early warning, the present invention realize the automation and intelligence of working face water damage early warning, are suitable for during working face extraction Analysis of possibility of water inrush early warning.
First according to the geologic information of stope, establish the 3 D resistivity model of coal measure strata, model be comprising The 3D data volume of resistivity value.System is monitored in stope arrangement mine electric method, starts downhole data acquisition, will monitor Data synchronous transfer is to ground-based server.Monitoring number is recorded with integer i, starts season i=0 in monitoring, completes primary monitoring I adds 1 automatically later, and the monitoring data of acquisition are stored with array d (i).Ground-based server is accessed automatically by process control, Obtain real-time monitoring data d (i).Using coal measure strata resistivity models as initial model m0, three-dimensional electricity is carried out to monitoring data Resistance rate inverting, the objective function of inverting are as follows:
M is the model for inverting iteration in formula, and d (m) is resulting to given conductivity model m progress forward modelling Data, d (i) are real-time monitoring data, and β is regularization parameter, and W is model regularization matrix.In order to obtain objective function Φ (m) minimum value obtains more new formula with above formula is to model parameter m derivation and to enable it be zero:
(JTJ+βWTW) δ m=- (JT(QA-1q-d(i))+βWTW(m-m0)) (2)
WhereinFor Jacobian matrix.Above formula is solved using quasi- Gauss-Newton method, obtains update Inverse model.The 3 D resistivity model that inverting obtains is stored to ground-based server, inverse model is deposited with array m (i) Storage.
When just starting monitoring, need first to carry out first monitoring data d (1) and second batch monitoring data d (2) anti- It drills, obtains corresponding inverse model m (1) and m (2), then carry out subsequent early warning judgement again.By adjacent batch monitoring data d (i-1) it is compared with the resulting 3 D resistivity model m (i-1) of d (i) inverting and m (i), judges two according to coordinate information Whether model has overlapping region.If overlapping region is not present in model m (i-1) and m (i), next group monitoring data are come back for Inverting, and by the inverse model newly obtained with last time obtain inverse model be compared, repeat the process until obtain Until being finished there are the two groups of adjacent inverse models or working face extraction of overlapping region.If there is coincidence in model m (i-1) and m (i) Region, then for overlapping region calculate two models resistivity differences δ m (i)=m (i)-m (i-1), and judge be in δ m (i) The no point existed less than 0.If meaning that gushing water hidden danger temporarily is not present there is no the point less than 0 in δ m (i), coming back for The inverse model that the inverse model newly obtained is obtained with the last time is compared for the inverting of next group monitoring data, and repeating should There are the two groups of adjacent inverse models or working face extraction that there is the point less than 0 in overlapping region and δ m (i) up to obtaining for process Until finishing.For the point in δ m (i) less than 0, the local minimum r (i) of δ m (i) is found out according to coordinate information, saves part Minimum value coordinate position corresponding with its sends the three-level warning information comprising local minimum coordinate position: " at so-and-so position There are gushing water hidden danger!".
It is primary before judging whether there is to calculate the local minimum r (i-1) obtained.Local minimum r (i- if it does not exist 1) inverting of next group observation data, is then come back for, and judges whether to need to send three-level warning information, it is straight to repeat the process Until obtaining two groups of front and backs adjacent local minimums or working face extraction and finishing.Local minimum r (i-1) if it exists, then Judge that the local minimum r (i) newly obtained and the preceding primary local minimum r (i-1) obtained that calculates are overlapped with the presence or absence of position Point.If the point that position is overlapped is not present in local minimum r (i) and r (i-1), the anti-of next group observation data is come back for It drills, and judges whether to need to send three-level warning information, repeat the process until obtaining two groups of phases of the point being overlapped there are position Until adjacent local minimum or working face extraction finish.If local minimum r (i) and r (i-1) there are the point that position is overlapped, Save the coordinate position being overlapped, send the second level warning information comprising being overlapped coordinate position: " analysis of possibility of water inrush increases at so-and-so position Greatly!".
For the point that position in two groups of local minimums is overlapped, difference s (i)=r (i)-r of two groups of local minimums is calculated (i-1), and judge in s (i) with the presence or absence of the point less than 0.If coming back for next group there is no the point less than 0 in s (i) The inverting of data is observed, and judges whether to need to send three-level warning information and second level warning information, repeats the process until obtaining Until the s (i) or working face extraction that must have the point less than 0 are finished.If there is the point less than 0 in s (i), the point less than 0 is saved Coordinate position, send include the coordinate position level-one warning information: " analysis of possibility of water inrush further increases at so-and-so position, It is recommended that investigation!".
Judge whether back production finishes working face, if then terminating this early warning, if otherwise coming back for next group observation number According to inverting, and judge whether to need to send warning information at different levels, repeat the process until working face extraction finishes.
Specific embodiment described herein is only an example for the spirit of the invention.The neck of technology belonging to the present invention The technical staff in domain can make various modifications or additions to the described embodiments or replace by a similar method In generation, however, it does not deviate from the spirit of the invention or beyond the scope of the appended claims.

Claims (2)

1. a kind of working face water damage real-time graded method for early warning based on mine electric method monitoring characterized by comprising
Initial modeling procedure establishes the initial 3 D resistivity model of coal measure strata according to the geologic information of stope;
Inverting warning step obtains inverting using the initial 3 D resistivity model inversion each batch underground Real-time Monitoring Data 3 D resistivity model, when the inversion of three dimensional resistivity models of two neighboring batch are overlapped and are existed in differential resistivity value sequence When negative value, it is judged as that there are risk and early warning, otherwise continues inverting next batch underground Real-time Monitoring Data;
Wherein, the differential resistivity value sequence is the resistivity value and last time inversion of three dimensional electricity of this inversion of three dimensional resistivity models The difference of the resistivity value of resistance rate model;
Wherein, the inverting warning step further include:
Three-level early warning sub-step finds out the office of sequence of differences according to coordinate information for the point of small Yu 0 in differential resistivity value sequence Portion's minimum value saves local minimum coordinate position corresponding with its;It is pre- to send the three-level comprising local minimum coordinate position Alert information;
Second level early warning sub-step judges whether is local minimum that local minimum that this is calculated and last computation obtain There are position coincidence points, if then saving the coordinate position of coincidence, send the second level warning information comprising being overlapped coordinate position;
Level-one early warning sub-step, it is primary before being subtracted with the local minimum newly obtained to calculate the office obtained for position coincidence point Portion's minimum value calculates the difference of local minimum twice, judges whether there is point of the difference less than 0, and if it exists, then save the point Location information, and send include the coordinate position level-one warning information.
2. a kind of working face water damage real-time graded early warning system based on mine electric method monitoring characterized by comprising
Initial modeling module establishes the initial 3 D resistivity model of coal measure strata according to the geologic information of stope;
Inverting warning module obtains inverting using the initial 3 D resistivity model inversion each batch underground Real-time Monitoring Data 3 D resistivity model, when the inversion of three dimensional resistivity models of two neighboring batch are overlapped and are existed in differential resistivity value sequence When negative value, it is judged as that there are risk and early warning, otherwise continues inverting next batch underground Real-time Monitoring Data;
Wherein, the differential resistivity value sequence is the resistivity value and last time inversion of three dimensional electricity of this inversion of three dimensional resistivity models The difference of the resistivity value of resistance rate model;
Wherein, the inverting warning module further include:
Three-level early warning submodule finds out the office of sequence of differences according to coordinate information for the point of small Yu 0 in differential resistivity value sequence Portion's minimum value saves local minimum coordinate position corresponding with its;It is pre- to send the three-level comprising local minimum coordinate position Alert information;
Second level early warning submodule judges whether is local minimum that local minimum that this is calculated and last computation obtain There are position coincidence points, if then saving the coordinate position of coincidence, send the second level warning information comprising being overlapped coordinate position;
Level-one early warning submodule, it is primary before being subtracted with the local minimum newly obtained to calculate the office obtained for position coincidence point Portion's minimum value calculates the difference of local minimum twice, judges whether there is point of the difference less than 0, and if it exists, then save the point Location information, and send include the coordinate position level-one warning information.
CN201711265780.2A 2017-12-05 2017-12-05 Water damage real-time graded method for early warning and system based on mine electric method monitoring Active CN108104876B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711265780.2A CN108104876B (en) 2017-12-05 2017-12-05 Water damage real-time graded method for early warning and system based on mine electric method monitoring

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711265780.2A CN108104876B (en) 2017-12-05 2017-12-05 Water damage real-time graded method for early warning and system based on mine electric method monitoring

Publications (2)

Publication Number Publication Date
CN108104876A CN108104876A (en) 2018-06-01
CN108104876B true CN108104876B (en) 2019-03-22

Family

ID=62209117

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711265780.2A Active CN108104876B (en) 2017-12-05 2017-12-05 Water damage real-time graded method for early warning and system based on mine electric method monitoring

Country Status (1)

Country Link
CN (1) CN108104876B (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112253236B (en) * 2020-10-13 2022-09-23 中煤科工集团西安研究院有限公司 Method for cleaning data of mine electrical method monitoring data by utilizing correlation analysis
CN113903149B (en) * 2021-10-08 2022-12-27 河北煤炭科学研究院有限公司 Mine water inrush early warning system and method and terminal equipment

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101382070A (en) * 2007-09-03 2009-03-11 中国石油天然气集团公司 Electromagnetical method for dynamically monitoring oil reservoir injection-production
CN102495428A (en) * 2011-12-12 2012-06-13 山东大学 Resistivity real-time imaging monitoring method and system for water-bursting geological disaster in construction period of underground engineering
CN105204071A (en) * 2015-07-31 2015-12-30 淮南矿业(集团)有限责任公司 Method for determining water-flowing fractured zone of overlaying rock stratum to be monitored
CN105223612A (en) * 2015-06-10 2016-01-06 中国矿业大学 A kind of coal mine flood prediction and evaluation method based on earthquake information
CN105277992A (en) * 2015-12-02 2016-01-27 山东科技大学 Method for dynamically monitoring stope bottom plate water regime through equator dipole device
JP2017037515A (en) * 2015-08-11 2017-02-16 能美防災株式会社 Evacuation guidance display device inside tunnel
CN106768736A (en) * 2016-11-24 2017-05-31 中国科学技术大学 A kind of subway tunnel monitoring of leakage and early warning system and its method
CN107144891A (en) * 2017-06-12 2017-09-08 山东大学 The monitoring of water burst precursor information dash forward with merging early warning system and method in tunnel

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10557346B2 (en) * 2016-04-11 2020-02-11 Schlumberger Technology Corporation Methods for determining saturation parameters of a formation

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101382070A (en) * 2007-09-03 2009-03-11 中国石油天然气集团公司 Electromagnetical method for dynamically monitoring oil reservoir injection-production
CN102495428A (en) * 2011-12-12 2012-06-13 山东大学 Resistivity real-time imaging monitoring method and system for water-bursting geological disaster in construction period of underground engineering
CN105223612A (en) * 2015-06-10 2016-01-06 中国矿业大学 A kind of coal mine flood prediction and evaluation method based on earthquake information
CN105204071A (en) * 2015-07-31 2015-12-30 淮南矿业(集团)有限责任公司 Method for determining water-flowing fractured zone of overlaying rock stratum to be monitored
JP2017037515A (en) * 2015-08-11 2017-02-16 能美防災株式会社 Evacuation guidance display device inside tunnel
CN105277992A (en) * 2015-12-02 2016-01-27 山东科技大学 Method for dynamically monitoring stope bottom plate water regime through equator dipole device
CN106768736A (en) * 2016-11-24 2017-05-31 中国科学技术大学 A kind of subway tunnel monitoring of leakage and early warning system and its method
CN107144891A (en) * 2017-06-12 2017-09-08 山东大学 The monitoring of water burst precursor information dash forward with merging early warning system and method in tunnel

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
《煤矿井下含/导水构造三维电阻率反演成像技术》;鲁晶津;《煤炭学报》;20160315;第41卷(第3期);第688页右栏至第690页左栏

Also Published As

Publication number Publication date
CN108104876A (en) 2018-06-01

Similar Documents

Publication Publication Date Title
CN106150485B (en) A kind of drilling well leakage forecasting system based on digital rock mass
CN105607146B (en) A kind of quantitatively characterizing method of meandering river sand body scale
CN107861157B (en) A kind of underground water seal cave depot operation phase micro seismic monitoring method
CN103615238B (en) A kind of cavernous formation dual laterolog physical simulating device of scaled down and experimental technique
CN108104876B (en) Water damage real-time graded method for early warning and system based on mine electric method monitoring
CN109236273B (en) Dynamic data processing method for oil field development and production
CN112330927B (en) Multi-parameter integrated karst collapse monitoring and early warning device and working method thereof
CN103064124A (en) Ratio method for correcting topographic influence in electromagnetic prospecting
CN110532507B (en) Method for improving drilling rate of compact oil reservoir of horizontal well
CN111563653A (en) Early warning construction method for water-rich broken stratum of underground engineering
CN104459800A (en) Sand body pinch-out prediction method and device
JP5557646B2 (en) Estimation method of formation parameters
CN205066964U (en) Working face bottom plate gushing water temperature distribution formula monitoring system
CN109991667A (en) A kind of quick judgment method of aqueous interlayer hydraulic connection
CN113433595B (en) Advanced prediction method based on natural electric field tunnel crack water
CN113586160B (en) Comprehensive-caving face water condition monitoring system based on artificial intelligence
CN108757036B (en) The system and method for liter state is led for monitoring stope bottom plate artesian water
CN108228952B (en) Electrical anisotropy numerical simulation method for detecting goaf fracture zone by electrical method
CN104749645B (en) A kind of method and apparatus of coal seam substrate interface undulation form detection in ash top difficult to understand
CN106405672B (en) A kind of TBM construction tunnels three-dimensional induced polarization advanced prediction disturbs minimizing technology
CN114087019A (en) Method for preventing and controlling water damage area of thick and non-homogeneous sandstone
CN104749641A (en) Method and device for recognizing mined-out areas
CN117726045B (en) Water inflow prediction method and system based on hydrological field inversion along tunnel
CN111123404A (en) Data fusion method for roadway advanced detection by earthquake and direct current method
KALENCHUK et al. 28 Downie Slide, British Columbia, Canada

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