CN108228952A - A kind of electrical anisotropy method for numerical simulation of electrical survey (-ing) goaf fissure zone - Google Patents
A kind of electrical anisotropy method for numerical simulation of electrical survey (-ing) goaf fissure zone Download PDFInfo
- Publication number
- CN108228952A CN108228952A CN201710283502.3A CN201710283502A CN108228952A CN 108228952 A CN108228952 A CN 108228952A CN 201710283502 A CN201710283502 A CN 201710283502A CN 108228952 A CN108228952 A CN 108228952A
- Authority
- CN
- China
- Prior art keywords
- numerical simulation
- electrical
- goaf
- earth
- fissure zone
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH DRILLING; MINING
- E21B—EARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
Abstract
A kind of electrical anisotropy method for numerical simulation of electrical survey (-ing) goaf fissure zone, belongs to mine electromagnetic surveying method for numerical simulation.Method and step:1st, the electrical method theory deduction based on goaf fissure zone electrical anisotropy feature;2nd, according to electrical method theory deduction, collect Numerical Simulation Program;3rd, according to the practical geologic information in goaf, the required earth-electricity model of numerical simulation is established;4th, numerical simulation is carried out, to each current electrode, using A as the center of circle, is rotated by radius of AB, one numerical value of every 15 degree of calculating calculates each electrodes of A B;5th, change the length of AB, repeat (4) and calculated, until the length of AB meets the needs of design;6th, it after result of calculation is exported, is plotted in polar coordinate system, the shape and the length of AB being distributed according to apparent resistivity tentatively judge stratum electrical anisotropy feature.Advantage:The apparent resistivity of different direction can be measured, prevents and administer to provide detailed geologic information for mine water.
Description
Technical field
The present invention relates to a kind of mine electromagnetic surveying method for numerical simulation, particularly a kind of electrical survey (-ing) goaf fissure zone
Electrical anisotropy method for numerical simulation.
Background technology
It takes precautions against and the containment great water inrush accident of coal mine is always all kinds of coal mining enterprises, Coalmine Safety Supervision department, coal life
Production administrative department, management organization of governments at all levels implement safety responsibility, carry out the most important thing of hidden troubles removing work.In recent years
Come, great water inrush accident type, goaf water leak accident accounts for all kinds of more than 70% coal mine flood accidents.The permeable key in goaf
Point is that fissure zone forms conduit pipe, and the ponding in goaf is caused to penetrate into mine and causes mine water disaster accident.
Yin Changchun has delivered the method that layered medium electrical anisotropy feature is calculated using DC electrical method in 1999;
Shen Jinsong on this basis, in 2008 delivered calculate oil fracture reservoir electrical anisotropy feature;
On the basis of existing technology, the electrical anisotropy of calculating water flowing fractured zone delivered in 2016 is special by inventor
Sign, theory deduction specific method step are as follows:
▽ × E=0, ▽ J=0 (1)
▽ × H=J, ▽ B=0 (2)
Wherein:
Wherein, resistivityThe electrical information of fissure zone can be described.
Magnetic field and current density can be expressed as:
Wherein, TJAnd PHIt is electric field and magnetic field spiral type and pole respectively to type scalar potential;THAnd PJIt is magnetic field and electric field respectively
Spiral type and pole are to type scalar potential;
By formula:
Therefore,
Similarly, the expression formula of E can be obtained;By Fourier transformation and boundary condition;The table of horizontal component of electric field and apparent resistivity
Up to formula:
Wherein, G is K factor;ρaFor apparent resistivity;R is between coordinate origin O to MN midpoints
Distance;ErFor MN midpoints electric field;ExFor MN midpoints x directions electric field, EyFor MN midpoints y directions electric field.
Invention content
The technical problems to be solved by the invention are to provide a kind of electrical anisotropy number of electrical survey (-ing) goaf fissure zone
It is worth analogy method, the trend or tendency of development characteristics, i.e. crack according to water flowing fractured zone crack provide further crack
Information provides specifically geology data for the slurry injection technique in the preventing and administer of goaf water leak accident.
The object of the present invention is achieved like this:Detect method for numerical simulation:Based on electrical survey (-ing) theory, goaf is considered
Water flowing fractured zone electrical anisotropy feature, using between technique study fissure zone development condition and the apparent resistance distribution of numerical simulation
Relationship;It is as follows:
(1) electrical method theory deductions of the based on goaf fissure zone electrical anisotropy feature;
(2) is according to electrical method theory deduction, and collect Numerical Simulation Program;Compilation Numerical Simulation Program is divided into four modules:①
Model parameter inputs;2. electrical anisotropy earth-electricity model judges;3. electrical anisotropy earth-electricity model response computation;4. result of calculation
Output;
(3) establishes the required earth-electricity model of numerical simulation, the mined out of reality according to the practical geologic information in goaf
Area is reduced to three floor, wherein, middle layer is that desired homogeneous is situated between for the fissure zone of electrical anisotropy and fissure zone the upper and lower
Matter;
(4) carries out numerical simulation for the earth-electricity model that step (3) is established;The numerical simulation is:To each confession
Electrode using A as the center of circle, is rotated by radius of AB, one numerical value of every 15 degree of calculating, to each electrodes of A B, always
24 apparent resistance values are calculated altogether;
(5) changes the length of AB, repeats step (4) and is calculated, until the needs that the length of AB meets investigation depth are
Only;
(6) the compiled Numerical Simulation Programs of carry out numerical simulation, after result of calculation output module exports result of calculation,
It is plotted in polar coordinate system, pole is in oval figure shape distribution to figure apparent resistivity, and long axis direction is directed toward the trend of fractue spacing, and ellipse
Round property is related with crack tendency;Fractue spacing using the oval feature parameter estimation subsurface formations of orientation apparent resistivity is joined
Number.
Advantageous effect, as a result of said program, the present invention is based on electromagnetic prospecting theories, consider goaf fissure zone
Electrical anisotropy is utilized apparent resistivity and describes goaf fissure zone calculation formulaBy analyzing in earth's surface
The apparent resistance of measurement is worth distribution, can obtain the information of goaf fissure zone, is the detection of electrical method goaf water leak accident
With prevent from providing fundamental basis.
Advantages of the present invention:
1st, the present invention considers goaf fissure zone electrical anisotropy feature, obtains electricity from basic Theory of Electromagnetic Field
Relationship between magnetic response and apparent resistivity.
2nd, method proposed by the present invention can utilize the distributional pattern feature estimation water flowing fractured zone of earth's surface apparent resistivity to split
The trend and tendency of seam prevent and administer to provide detailed geologic information for mine water.
3rd, the apparent resistivity of different direction is measured, may determine that goaf fissure zone is sent out according to the distribution of apparent resistivity
State is educated, the detection for goaf fissure zone provides theoretical foundation, significant.
Description of the drawings
Fig. 1 is the goaf geological model figure of the present invention.
Fig. 2 (a) is the goaf earth-electricity model simplification figure of the present invention.
Fig. 2 (b) is measurement method schematic diagram in goaf of the present invention.
The fracture strike that Fig. 3 (a) is the present invention is consistent with Y-axis, crack tendency:The goaf earth-electricity model figure of α=45 °.
The fracture strike that Fig. 3 (b) is the present invention is consistent with Y-axis, crack tendency:The goaf earth-electricity model figure of α=60 °.
The fracture strike that Fig. 3 (c) is the present invention is consistent with Y-axis, crack tendency:The goaf earth-electricity model figure of α=90 °.
Fig. 4 (a) is the earth-electricity model numerical simulation result figure of Fig. 3 (a).
Fig. 4 (b) is the earth-electricity model numerical simulation result figure of Fig. 3 (b).
Fig. 4 (c) is the earth-electricity model numerical simulation result figure of Fig. 3 (c).
Fig. 5 (a) is the fracture strike of the present invention:β=0 °, tendency:The goaf earth-electricity model figure of α=45 °.
Fig. 5 (b) is the fracture strike of the present invention:β=30 °, tendency:The goaf earth-electricity model figure of α=45 °.
Fig. 5 (c) is the fracture strike of the present invention:β=60 °, tendency:The goaf earth-electricity model figure of α=45 °.
Fig. 6 (a) is the earth-electricity model numerical simulation figure of Fig. 5 (a).
Fig. 6 (b) is the earth-electricity model numerical simulation figure of Fig. 5 (b).
Fig. 6 (c) is the earth-electricity model numerical simulation figure of Fig. 5 (c).
Specific embodiment
The detection method for numerical simulation of the present invention:Based on electrical survey (-ing) theory, consider goaf water producing fractures electrification it is each to
Different in nature feature utilizes the relationship between technique study fissure zone development condition and the apparent resistance distribution of numerical simulation;Specific steps
It is as follows:
(1) electrical method theory deductions of the based on goaf fissure zone electrical anisotropy feature;
(2) is according to electrical method theory deduction, and collect Numerical Simulation Program;Compilation Numerical Simulation Program is divided into four modules:①
Model parameter inputs;2. electrical anisotropy earth-electricity model judges;3. electrical anisotropy earth-electricity model response computation;4. result of calculation
Output;
(3) establishes the required earth-electricity model of numerical simulation, the mined out of reality according to the practical geologic information in goaf
Area is reduced to three floor, wherein, middle layer is that desired homogeneous is situated between for the fissure zone of electrical anisotropy and fissure zone the upper and lower
Matter;
(4) carries out numerical simulation for the earth-electricity model that step (3) is established;The numerical simulation is:To each confession
Electrode using A as the center of circle, is rotated by radius of AB, one numerical value of every 15 degree of calculating, to each electrodes of A B, always
24 apparent resistance values are calculated altogether;
(5) changes the length of AB, repeats step (4) and is calculated, until the needs that the length of AB meets investigation depth are
Only;
(6) the compiled Numerical Simulation Programs of carry out numerical simulation, after result of calculation output module exports result of calculation,
It is plotted in polar coordinate system, pole is in oval figure shape distribution to figure apparent resistivity, and long axis direction is directed toward the trend of fractue spacing, and ellipse
Round property is related with crack tendency;Fractue spacing using the oval feature parameter estimation subsurface formations of orientation apparent resistivity is joined
Number.
Embodiment 1:The scalar potential of the max equation of electrical anisotropy is expressed as spiral type and pole to type scalar potential.
By fast algorithm, current density and magnetic vector are taken to the two-dimensional Fourier transform in x and y directions, must be shown up point
The wave-number domain expression formula of amount.
Max fundamental equation is solved in wave-number domain, the elementary solution of electromagnetic field scalar can be obtained.
According to electromagnetic field in the bed boundary condition of continuity, corresponding relational expression is obtained.
Current source is inserted in earth's surface, considers the coupling condition derived from earth's surface.
Source coupling condition based on the bed boundary condition of continuity and earth's surface, construction calculate recurrence relation, obtain scalar potential
The magnetic field of current density and scalar potential.
It is obtained using two-dimensional inverse Fourier transform, obtains electric field and the magnetic field of spatial domain;Wherein, it is carried out using horizontal component of electric field
Express the apparent resistivity of anisotropic formation.
Wherein,E is electric field, and x and y are respectively the horizontal seat of measurement point
Mark and ordinate, G are K factor, G=π L2, L=AB/2, AB are electrode spacings, and I is electric current.
Numerical simulation studies the relationship between crack tendency and distribution features of apparent resistivity first, with selecting three goafs
Electric model, wherein, fracture strike is all along Y direction, if moving towards angle as angle β=0 ° with Y direction, crack dip is not
Three same goaf earth-electricity models, wherein (a) inclination angle is α=45 °, (b) inclination angle is α=60 °, and (c) inclination angle is β=90 °, such as
Shown in Fig. 4, numerical simulation result is as shown in figure 5, the morphological feature of apparent resistivity distribution changes with the change at inclination angle.
Earth-electricity model Fig. 4 and corresponding numerical simulation result Fig. 5 is compared, it can be with it is concluded that (1):The pole of apparent resistivity distribution to figure,
Its shape is related to inclination alpha.
The relationship between fracture strike and distribution features of apparent resistivity is studied again, selects three goaf earth-electricity models, such as
Fig. 6 wherein crack tendencies are all α=45 ° always, and fracture strike is distributed as, (a) trend:β=0 °;(b) it moves towards:β=30 °;
(c) it moves towards:β=60 °, numerical simulation result as shown in fig. 6, apparent resistivity distribution elliptical symmetry axis be directed toward with
Move towards related.Earth-electricity model Fig. 5 and numerical simulation result Fig. 6 is compared, is obtained with draw a conclusion (3):The ellipse of apparent resistivity distribution
The symmetry axis of shape is directed toward related to trend.
Technical scheme of the present invention is further described in detail with reference to the accompanying drawings and detailed description.
Fig. 2 (a) is goaf earth-electricity model simplification figure;Fig. 2 (b) is goaf measurement method schematic diagram;Using A as the center of circle,
Using AB as radius, a value is calculated every 15 degree, 24 values were obtained around one week.
It is the fracture strike crack (along Y direction, if moving towards angle as angle β=0 ° with Y direction) always in Fig. 3
Three different goaf earth-electricity models of inclination angle, wherein (a) inclination angle be α=45 °, (b) inclination angle be α=60 °, (c) inclination angle for β=
90°。
Fig. 3 (a):It is inclined in crack:α=45 °, trend are consistent with Y direction:β=0 °;(b):It is inclined in crack:α=60 ° are walked
To consistent with Y direction:β=0 °;(c):It is inclined in crack:α=90 °, trend are consistent with Y direction:β=0 °.
Fig. 4 is the earth-electricity model numerical simulation result of Fig. 3 descriptions, wherein (a):It is inclined in crack:α=45 °, trend:β=
0 °, earth-electricity model numerical simulation result;(b):It is inclined in crack:α=60 °, trend:β=0 °, earth-electricity model numerical simulation result;
(c):It is inclined in crack:α=90 °, trend:β=0 °, earth-electricity model numerical simulation result.
Fig. 4 is the earth-electricity model numerical simulation result of Fig. 3 descriptions, wherein (a):It is inclined in crack:α=45 °, trend:β=
0 °, earth-electricity model numerical simulation result;(b):It is inclined in crack:α=60 °, trend:β=0 °, earth-electricity model numerical simulation result;
(c):It is inclined in crack:α=90 °, trend:β=0 °, earth-electricity model numerical simulation result.
It is inclined to goaf earth-electricity model figure in Fig. 5 cracks;When:(a) it moves towards:β=0 °, tendency:α=45 °;(b) it moves towards:β=
30 °, tendency:α=45 °;(c) it moves towards:β=60 °, tendency:α=45 °.
Fig. 6 is inclined to always for crack, is all α=45 °.Wherein, (a) is moved towards:β=0 °;(b) it moves towards:β=30 °;(c) it walks
To:β=60 °.
Trend always (along Y direction, if moving towards angle as angle β=0 ° from Y direction) and crack dip is different
Three goaf earth-electricity models, wherein (a) inclination angle is α=45 °, (b) inclination angle is α=60 °, and (c) inclination angle is β=90 °.
Fig. 6 is the numerical simulation result of the earth-electricity model described in Fig. 5.Wherein, (a) is moved towards:β=0 °, tendency:α=
45 °, earth-electricity model numerical simulation result;(b) it moves towards:β=30 °, tendency:α=45 °, earth-electricity model numerical simulation result;(c)
Trend:β=60 °, tendency:α=45 °, earth-electricity model numerical simulation result.
Claims (1)
1. a kind of electrical anisotropy method for numerical simulation of electrical survey (-ing) goaf fissure zone, it is characterised in that:Detect Numerical-Mode
Plan method:Based on electrical survey (-ing) theory, consider goaf water flowing fractured zone electrical anisotropy feature, utilize the method for numerical simulation
Study the relationship between fissure zone development condition and apparent resistance distribution;It is as follows:
(1) electrical method theory deductions of the based on goaf fissure zone electrical anisotropy feature;
(2) is according to electrical method theory deduction, and collect Numerical Simulation Program;The compilation Numerical Simulation Program is divided into four modules:
1. model parameter inputs;2. electrical anisotropy earth-electricity model judges;3. electrical anisotropy earth-electricity model response computation;4. calculate knot
Fruit exports;
(3) establishes the required earth-electricity model of numerical simulation, practical goaf letter according to the practical geologic information in goaf
Three layers are turned to, wherein, middle layer is desired homogeneous medium for the fissure zone of electrical anisotropy and fissure zone the upper and lower;
(4) carries out numerical simulation for the earth-electricity model that step (3) is established;The numerical simulation is:To each power supply electricity
Pole using A as the center of circle, is rotated by radius of AB, one numerical value of every 15 degree of calculating, total to amount to each electrodes of A B
Calculate 24 apparent resistance values;
(5) changes the length of AB, repeats step (4) and is calculated, until the length of AB meets the needs of investigation depth;
(6) compilations Numerical Simulation Program carries out numerical simulation, after result of calculation output module exports result of calculation, is plotted in
In polar coordinate system, pole is in that oval figure shape is distributed to figure apparent resistivity, the trend of long axis direction direction fractue spacing, and elliptical property
Matter is related with crack tendency;Utilize the fractue spacing parameter of the oval feature parameter estimation subsurface formations of orientation apparent resistivity.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710283502.3A CN108228952B (en) | 2017-04-26 | 2017-04-26 | Electrical anisotropy numerical simulation method for detecting goaf fracture zone by electrical method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710283502.3A CN108228952B (en) | 2017-04-26 | 2017-04-26 | Electrical anisotropy numerical simulation method for detecting goaf fracture zone by electrical method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108228952A true CN108228952A (en) | 2018-06-29 |
CN108228952B CN108228952B (en) | 2020-08-14 |
Family
ID=62656667
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710283502.3A Active CN108228952B (en) | 2017-04-26 | 2017-04-26 | Electrical anisotropy numerical simulation method for detecting goaf fracture zone by electrical method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108228952B (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108982594A (en) * | 2018-08-06 | 2018-12-11 | 山东建筑大学 | The ERT pollution detection system and setting method of detection cycle adaptive setting |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2006135510A1 (en) * | 2005-06-09 | 2006-12-21 | Exxonmobil Upstream Research Company | Method for determining earth vertical electrical anisotropy in marine electromagnetic surveys |
CN102808644A (en) * | 2012-07-17 | 2012-12-05 | 西安科技大学 | Coal and gas co-mining method based on mining fissure elliptic paraboloid zone |
CN103513281B (en) * | 2013-10-23 | 2015-12-02 | 中国矿业大学 | Solid filling coal mining overlying strata water flowing fractured zone development height method for predicting |
CN105223337A (en) * | 2015-10-22 | 2016-01-06 | 重庆大学 | The three-dimensional simulation method of goaf cracked district form |
CN105868472A (en) * | 2016-03-29 | 2016-08-17 | 辽宁工程技术大学 | Gob air leakage flow field dynamic numerical simulation method based on deformation geometry |
-
2017
- 2017-04-26 CN CN201710283502.3A patent/CN108228952B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2006135510A1 (en) * | 2005-06-09 | 2006-12-21 | Exxonmobil Upstream Research Company | Method for determining earth vertical electrical anisotropy in marine electromagnetic surveys |
CN102808644A (en) * | 2012-07-17 | 2012-12-05 | 西安科技大学 | Coal and gas co-mining method based on mining fissure elliptic paraboloid zone |
CN103513281B (en) * | 2013-10-23 | 2015-12-02 | 中国矿业大学 | Solid filling coal mining overlying strata water flowing fractured zone development height method for predicting |
CN105223337A (en) * | 2015-10-22 | 2016-01-06 | 重庆大学 | The three-dimensional simulation method of goaf cracked district form |
CN105868472A (en) * | 2016-03-29 | 2016-08-17 | 辽宁工程技术大学 | Gob air leakage flow field dynamic numerical simulation method based on deformation geometry |
Non-Patent Citations (2)
Title |
---|
苏本玉 等: "采空区裂隙带电各向异性特征研究", 《2014年中国地球科学联合学术年会--专题22:煤炭资源与矿山安全勘查技术论文集》 * |
陈爽爽 等: "断裂带电各向异性直流电阻率法探测研究", 《中州煤炭》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108982594A (en) * | 2018-08-06 | 2018-12-11 | 山东建筑大学 | The ERT pollution detection system and setting method of detection cycle adaptive setting |
CN108982594B (en) * | 2018-08-06 | 2020-09-25 | 山东建筑大学 | ERT pollution detection system with detection period set in self-adaptive mode and setting method |
Also Published As
Publication number | Publication date |
---|---|
CN108228952B (en) | 2020-08-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110231659B (en) | Railway tunnel aviation electromagnetic detection method with large buried depth change | |
Su et al. | Research of the electrical anisotropic characteristics of water-conducting fractured zones in coal seams | |
CN104198539A (en) | Method for detecting and evaluating grouting reinforcement effect of seam floor | |
CN106919770A (en) | A kind of damage variable based on numerical simulation determines method | |
Wang et al. | Potential of hydraulic tomography in exploring the preferential flowpaths of water inrush in coal mine areas | |
CN106443189A (en) | Method and system for three-dimensional detection on earth electrode field and surrounding soil resistivity | |
CN108760825A (en) | It is adapted to across the hole resistivity CT imaging methods of three-dimensional of the arbitrary cloth hole of urban environment | |
Shi et al. | Application of three-dimensional high-density resistivity method in roof water advanced detection during working stope mining | |
CN102182437B (en) | Method for determining and eliminating hydraulic fracture stress boundary of coal mine underground drilling | |
CN107065019A (en) | Applied to road disaster and the 3 D electromagnetic imaging device and application method that collapse detection | |
Gao et al. | Water detection within the working face of an underground coal mine using 3D electric resistivity tomography (ERT) | |
CN108804728A (en) | Horizontal well stratum reservoir grading analysis method and computer readable storage medium | |
CN109611073A (en) | A kind of densification grease horizontal well mud stone can the preferred method and device in fractured layer position | |
CN108873073B (en) | A kind of across hole resistivity tomography method of three-dimensional based on electrical method of network concurrency | |
CN108228952A (en) | A kind of electrical anisotropy method for numerical simulation of electrical survey (-ing) goaf fissure zone | |
CN106772630A (en) | A kind of earth lead source transient electromagnetic axial detection method | |
CN104407386B (en) | A kind of electromagnetic survey data collection and processing method and system | |
Vargemezis | 3D geoelectrical model of geothermal spring mechanism derived from VLF measurements: A case study from Aggistro (Northern Greece) | |
CN104749645B (en) | A kind of method and apparatus of coal seam substrate interface undulation form detection in ash top difficult to understand | |
CN1828016A (en) | Bipolar-bipolar detection method for mine working face floor strata | |
Magnúsdóttir | Fracture Characterization in Geothermal Reservoirs Using Time-lapse Electric Potential Data | |
CN113447991A (en) | Method and device for rebuilding underground electrical anomaly | |
Zhang et al. | Discussion on monitoring and characterising group drilling pumping test within a massive thickness aquifer using the time-lapse transient electromagnetic method (TEM). | |
Hu et al. | A parallel electric field penetration technique for identifying hazardous water sources beneath coal seam mining face floors | |
Wang et al. | Application of comprehensive geophysical exploration technology in the detection of fault structure in geothermal field target area |
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 |