CN111538108A - Original geostress field-based mine rock joint structure simulation method - Google Patents
Original geostress field-based mine rock joint structure simulation method Download PDFInfo
- Publication number
- CN111538108A CN111538108A CN202010534086.1A CN202010534086A CN111538108A CN 111538108 A CN111538108 A CN 111538108A CN 202010534086 A CN202010534086 A CN 202010534086A CN 111538108 A CN111538108 A CN 111538108A
- Authority
- CN
- China
- Prior art keywords
- joint structure
- field
- rock
- ground stress
- simulation
- 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.)
- Pending
Links
- 239000011435 rock Substances 0.000 title claims abstract description 47
- 238000000034 method Methods 0.000 title claims abstract description 26
- 238000004088 simulation Methods 0.000 title claims abstract description 25
- 238000005065 mining Methods 0.000 claims abstract description 22
- 238000011835 investigation Methods 0.000 claims abstract description 10
- 238000000547 structure data Methods 0.000 claims abstract description 7
- 238000009826 distribution Methods 0.000 claims description 8
- 238000005553 drilling Methods 0.000 claims description 4
- 230000003746 surface roughness Effects 0.000 claims description 4
- 238000009933 burial Methods 0.000 claims description 3
- 238000011156 evaluation Methods 0.000 abstract description 8
- 229910052500 inorganic mineral Inorganic materials 0.000 abstract description 5
- 239000011707 mineral Substances 0.000 abstract description 5
- 238000011161 development Methods 0.000 abstract description 4
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical compound [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 description 8
- 229910052742 iron Inorganic materials 0.000 description 4
- RYGMFSIKBFXOCR-UHFFFAOYSA-N Copper Chemical compound [Cu] RYGMFSIKBFXOCR-UHFFFAOYSA-N 0.000 description 3
- 230000006835 compression Effects 0.000 description 3
- 238000007906 compression Methods 0.000 description 3
- 229910052802 copper Inorganic materials 0.000 description 3
- 239000010949 copper Substances 0.000 description 3
- 238000002474 experimental method Methods 0.000 description 2
- 239000010438 granite Substances 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000000691 measurement method Methods 0.000 description 2
- 238000013461 design Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- SZVJSHCCFOBDDC-UHFFFAOYSA-N iron(II,III) oxide Inorganic materials O=[Fe]O[Fe]O[Fe]=O SZVJSHCCFOBDDC-UHFFFAOYSA-N 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 239000011044 quartzite Substances 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V20/00—Geomodelling in general
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/64—Geostructures, e.g. in 3D data cubes
- G01V2210/646—Fractures
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)
Abstract
The invention belongs to the technical field of mine geology, and particularly relates to a method for simulating a rock and ore joint structure based on an original geostress field, which is characterized by comprising the following steps of: (1) investigating a field joint structure; (2) establishing a three-dimensional model; (3) and (4) performing inversion numerical simulation on the ground stress structure. The method utilizes mechanical software to restore the mineral rock stress environment during the mining period, namely a rock model is built in the software, and stresses in different directions and different magnitudes are applied, when the development condition of the joint crack is matched with the existing data, the mechanical environment is considered to be close to the stress condition during the mining period, and the stress condition is loaded in other position areas, so that the joint structure condition of an unknown area is obtained by inversion by the method, a large amount of reliable structure data is obtained, and the accuracy of the collapse evaluation can be effectively improved; a large amount of simulation data can be obtained through limited data, and the problems of high cost, high difficulty and long period of joint structure field investigation are thoroughly solved.
Description
Technical Field
The invention belongs to the technical field of mine geology, and particularly relates to a method for simulating a rock and ore joint structure based on an original geostress field.
Background
The natural caving mining method is a safe and efficient mining method capable of realizing large-scale continuous ore removal, and is mainly suitable for mining thick and large ore bodies with relatively-developed joint cracks and steep dip angles. At present, more than 50 mines mined by a natural caving method in the world are mainly distributed in the countries such as the United states, Australia, Canada, Chile, south Africa, China and the like, and domestic application mines mainly comprise copper mine, valley copper mine and plain copper mine. The ore rock collapsibility evaluation is an important content of natural caving method mine exploitation design and research, and has decisive influence on mining stoping sequence, bottom-drawing direction, bottom-drawing area, ore removal mode, ore drawing control, cut wall presplitting, safe production, technical and economic indexes and the like. In the past, the ore rock cavibility evaluation mainly depends on the joint structure condition of field investigation, the physical and mechanical parameters of the ore rock mass, the regional ground stress field condition and the like, and because the field investigation work is difficult to cover the whole ore region and only a small amount of limited data can be utilized to carry out the ore rock cavibility evaluation, the ore rock cavibility evaluation result of the ore region is usually local and inaccurate. The conventional joint structure data survey is usually carried out in the disclosed underground rock mass engineering or through exploration drilling, so that the difficulty is high, the cost is high, the period is long, and in addition, the obtained data is limited and the structural characteristics of the rock mass of the whole mining area are difficult to reflect. In order to more accurately evaluate the collapsibility of the ore rocks in the whole mining area and guide the mine to reasonably determine the first mining position, the bottom mining direction and the like, the collapsibility evaluation of the ore rocks needs to be carried out by relying on a more comprehensive mining area joint structural model.
Disclosure of Invention
The invention aims to provide a method for simulating a rock joint structure based on an original geostress field, which is suitable for evaluating the caving property of rock ores during mining by adopting a natural caving method, mainly solves the problem of insufficient geological data of field engineering, and can obtain more related data by utilizing the method.
The purpose of the invention is realized by the following technical scheme:
the invention discloses a rock joint structure simulation method based on an original geostress field, which is characterized by comprising the following steps of:
(1) on-site joint structure survey
Firstly, carrying out on-site geological survey work through existing underground engineering and exploration drilling engineering of a mine to obtain distribution conditions of rock mass joint structures in different areas; the dominant joint group occurrence, the structural surface line density and the structural surface roughness characteristic parameters of the regional ore rock mass are known;
(2) building a three-dimensional model
Establishing a three-dimensional model of lithological structure of the mining area, and preparing for subsequent ground stress simulation;
(3) inversion numerical simulation of ground stress structure
Carrying out structure inversion simulation under different stresses on a local area in which joint structure data are obtained, namely applying a three-dimensional ground stress field to a three-dimensional model of a rock mass in a mining area, and referring to the ground stress field measured on the spot of the current mine by the initial value of the ground stress field; randomly changing the magnitude and direction of the applied ground stress field; and extracting failure surface attitude information by simulating failure modes of the three-dimensional model under different stress field conditions, recording the ground stress parameters when the characteristics of the rock mass structural surface appearing in the simulation are close to or reach characteristic parameters of field investigation, and simulating rock blocks in an unknown region in the ranges of similar burial depths and rock types by using the parameters so as to obtain the joint structure distribution conditions of other regions.
The invention has the advantages that:
the invention relates to a method for simulating a mineral rock joint structure based on an original geostress field, which is characterized in that a mechanical software is used for reducing a mineral rock stress environment during a mining period, namely a rock model is established in the software, and stresses in different directions and different magnitudes are applied, when the joint crack development condition is matched with the existing data, the mechanical environment is considered to be close to the stress condition during the mining period, and the stress condition is loaded in other position areas, so that the method is used for inverting to obtain the joint structure condition of an unknown area, thereby obtaining a large amount of reliable structure data and effectively increasing the accuracy of collapsibility evaluation; a large amount of simulation data can be obtained through limited data, and the problems of high cost, high difficulty and long period of joint structure field investigation are thoroughly solved.
Drawings
FIG. 1 is a perspective view of a fault and joint fracture.
FIG. 2 is a diagram of a stratigraphic configuration model of the present invention.
FIG. 3 is a plan view of the segment joint density of example 203 of the present invention.
Detailed Description
The following further describes the embodiments of the present invention with reference to the drawings.
As shown in fig. 1-3, the method for simulating the joint structure of the ore rock based on the original geostress field of the invention can accurately acquire the distribution condition of the joint structure of an unknown region by adopting a software simulation method, and is characterized by comprising the following steps:
(1) on-site joint structure survey
Firstly, carrying out on-site geological survey work through existing underground engineering and exploration drilling engineering of a mine to obtain distribution conditions of rock mass joint structures in different areas; the dominant joint group occurrence, the structural surface line density and the structural surface roughness characteristic parameters of the regional ore rock mass are known; obtaining a certain amount of joint structure distribution conditions of different areas through field geological exploration work, wherein the data have certain representativeness, such as the data are respectively distributed in areas with different depths and different lithologies;
(2) building a three-dimensional model
Establishing a three-dimensional model of lithological structure of the mining area, and preparing for subsequent ground stress simulation;
(3) inversion numerical simulation of ground stress structure
Carrying out structure inversion simulation under different stresses on a local area in which joint structure data are obtained, namely applying a three-dimensional ground stress field to a three-dimensional model of a rock mass in a mining area, and referring to the ground stress field measured on the spot of the current mine by the initial value of the ground stress field; randomly changing the magnitude and direction of an applied ground stress field, wherein the change increment of the magnitude of the ground stress is 2 MPa, and the change increment of the direction of the ground stress is 2 degrees; and extracting failure surface attitude information (tendency, inclination angle, roughness and the like) by simulating failure modes of the three-dimensional model under different stress field conditions, recording the ground stress parameters when the rock mass structural surface characteristics appearing in the simulation are close to or reach characteristic parameters of field investigation, and simulating rock blocks in an unknown region in the ranges of similar burial depths (+ -20 m) and rock types by using the parameters so as to obtain the joint structure distribution conditions of other regions.
Examples
(1) Investigation of anterior hill joint fissure
The method for investigating the fault structure of the iron ore in the front of the eye is large-scale geological measurement and gallery recording, the joint structure investigation uses a line measurement method and a three-dimensional laser scanning method, wherein the line measurement method is used for investigating and analyzing 28 roadways in the horizontal middle section of the iron ore in the front of the eye, namely 159, 177, 195, 213 and 303, and the investigating and analyzing roadway 2788 remaining meters, and the obtained field data specifically comprise the following contents: structural surface occurrence (including dip angle, tendency, structural surface grouping and the like), structural surface spacing, structural surface roughness, structural surface openness, structural surface continuity and the like. In the three-dimensional laser scanning method roadway survey, a maptekI-Site8200 three-dimensional laser scanner is used for scanning rock mass structural planes of 16 roadways at the middle sections of iron ores-159, -177, -195 and-213 in the front of eyes, and the total length of the scanning roadway is more than 3000 meters. The attitude projection of the fault and the joint fissure is shown in figure 1;
(2) building a three-dimensional model
The method for establishing the stratum structure model is based on a geological exploration profile, and utilizes datamine software to spatially connect the stratum and the structure of the exploration profile to generate a solid model, so that preparation is made for subsequent ground stress simulation. Lithology in the earth model of the iron ore of the mountain before the eye is mainly divided into three types, phyllite, magnetite quartzite and taigu granite. The stratigraphic structure model is shown in figure 2;
(3) inversion simulation of ground stress structure
Exporting the three-dimensional established by the datamine into a file with a dxf format, and importing the dxf file into finite element software for numerical simulation. Rock mechanical Properties according to the results of previous experiments, the uniaxial compression strength of magnetoquartzite was defined as 120MPa, the modulus of elasticity was defined as 22GPa, the Poisson's ratio was defined as 0.2, the uniaxial compression strength of phyllite was defined as 50MPa, the modulus of elasticity was defined as 18GPa, the Poisson's ratio was defined as 0.13, the uniaxial compression strength of ancient granite was defined as 30MPa, the modulus of elasticity was defined as 8.35GPa, and the Poisson's ratio was defined as 0.24. The method comprises the steps of performing fracturing simulation on rock in a combined change mode in three directions of XYZ, taking 20 degrees change in each direction as an experiment, simulating 5832 times under the action of a fixed force, then re-simulating by changing the magnitude of the force, comparing a simulation result with a known joint structure, and displaying that the joint structure development is closest to an actual measurement result when a stress of 10MPa is applied in the vertical direction and 180MPa is applied in the direction of 80-160 degrees, writing a simulated joint crack into a model to perform subsequent collapsibility and block size prediction, wherein a joint density plan in the middle section of-203 is shown in figure 3.
The invention relates to a method for simulating a mineral rock joint structure based on an original geostress field, which is characterized in that a mechanical software is used for reducing a mineral rock stress environment during a mining period, namely a rock model is established in the software, and stresses in different directions and different magnitudes are applied, when the joint crack development condition is matched with the existing data, the mechanical environment is considered to be close to the stress condition during the mining period, and the stress condition is loaded in other position areas, so that the method is used for inverting to obtain the joint structure condition of an unknown area, thereby obtaining a large amount of reliable structure data and effectively increasing the accuracy of collapsibility evaluation; a large amount of simulation data can be obtained through limited data, and the problems of high cost, high difficulty and long period of joint structure field investigation are thoroughly solved.
Claims (1)
1. A rock joint structure simulation method based on an original geostress field is characterized by comprising the following steps:
(1) on-site joint structure survey
Firstly, carrying out on-site geological survey work through existing underground engineering and exploration drilling engineering of a mine to obtain distribution conditions of rock mass joint structures in different areas; the dominant joint group occurrence, the structural surface line density and the structural surface roughness characteristic parameters of the regional ore rock mass are known;
(2) building a three-dimensional model
Establishing a three-dimensional model of lithological structure of the mining area, and preparing for subsequent ground stress simulation;
(3) inversion numerical simulation of ground stress structure
Carrying out structure inversion simulation under different stresses on a local area in which joint structure data are obtained, namely applying a three-dimensional ground stress field to a three-dimensional model of a rock mass in a mining area, and referring to the ground stress field measured on the spot of the current mine by the initial value of the ground stress field; randomly changing the magnitude and direction of the applied ground stress field; and extracting failure surface attitude information by simulating failure modes of the three-dimensional model under different stress field conditions, recording the ground stress parameters when the characteristics of the rock mass structural surface appearing in the simulation are close to or reach characteristic parameters of field investigation, and simulating rock blocks in an unknown region in the ranges of similar burial depths and rock types by using the parameters so as to obtain the joint structure distribution conditions of other regions.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010534086.1A CN111538108A (en) | 2020-06-12 | 2020-06-12 | Original geostress field-based mine rock joint structure simulation method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010534086.1A CN111538108A (en) | 2020-06-12 | 2020-06-12 | Original geostress field-based mine rock joint structure simulation method |
Publications (1)
Publication Number | Publication Date |
---|---|
CN111538108A true CN111538108A (en) | 2020-08-14 |
Family
ID=71980878
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010534086.1A Pending CN111538108A (en) | 2020-06-12 | 2020-06-12 | Original geostress field-based mine rock joint structure simulation method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111538108A (en) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104252009A (en) * | 2013-06-27 | 2014-12-31 | 中国石油化工股份有限公司 | Method for restoring brittleness and plasticity evolution history of petroleum and natural gas cap rocks |
CN104656124A (en) * | 2015-02-06 | 2015-05-27 | 山东大学 | Multi-parameter comprehensive rock burst predicting method based on geophysical exploration method |
CN105672971A (en) * | 2016-01-05 | 2016-06-15 | 中国石油大学(华东) | Method for forecasting opening pressure, opening sequence and water injection pressure of reservoir fissure |
CN106485015A (en) * | 2016-10-20 | 2017-03-08 | 辽宁工程技术大学 | A kind of determination method of mine tomography coverage |
CN108693572A (en) * | 2018-03-12 | 2018-10-23 | 太原理工大学 | A kind of ground stress field inversion method based on three-dimensional modeling |
-
2020
- 2020-06-12 CN CN202010534086.1A patent/CN111538108A/en active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104252009A (en) * | 2013-06-27 | 2014-12-31 | 中国石油化工股份有限公司 | Method for restoring brittleness and plasticity evolution history of petroleum and natural gas cap rocks |
CN104656124A (en) * | 2015-02-06 | 2015-05-27 | 山东大学 | Multi-parameter comprehensive rock burst predicting method based on geophysical exploration method |
CN105672971A (en) * | 2016-01-05 | 2016-06-15 | 中国石油大学(华东) | Method for forecasting opening pressure, opening sequence and water injection pressure of reservoir fissure |
CN106485015A (en) * | 2016-10-20 | 2017-03-08 | 辽宁工程技术大学 | A kind of determination method of mine tomography coverage |
CN108693572A (en) * | 2018-03-12 | 2018-10-23 | 太原理工大学 | A kind of ground stress field inversion method based on three-dimensional modeling |
Non-Patent Citations (5)
Title |
---|
刘健等: "深部花岗岩中天然开启节理剪切性能试验研究", 《地下空间与工程学报》 * |
刘宁等: "双江口水电站初始地应力场反演回归分析", 《山东大学学报(工学版)》 * |
周春梅,章泽军,李先福: "古构造应力场数值模拟方法探讨—以山西省晋城成庄煤矿区为例", 《地质力学学报》 * |
周春梅等: "古构造应力场数值模拟方法探讨――以山西省晋城成庄煤矿区为例", 《地质力学学报》 * |
张国强等: "基于神经网络结构分解与FLAC~(3D)的初始地应力场反演", 《大连大学学报》 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Stead et al. | A critical review of rock slope failure mechanisms: The importance of structural geology | |
Stead et al. | Understanding the mechanics of large landslides | |
Bai et al. | A general review on longwall mining‐induced fractures in near‐face regions | |
CN106485015B (en) | A kind of determination method of mine tomography coverage | |
Palmström et al. | The deformation modulus of rock masses—comparisons between in situ tests and indirect estimates | |
Dezayes et al. | Structure of the low permeable naturally fractured geothermal reservoir at Soultz | |
Kun et al. | Influence of the fault zone in shallow tunneling: A case study of Izmir Metro Tunnel | |
CN113820750B (en) | Method for quantitatively predicting mudstone structural cracks based on elastoplastic mechanics | |
CN110162882B (en) | Accurate control method for surrounding rock of coal mine tunnel crossing fault area | |
CN106372297A (en) | Method for determining safe vertical distance between shield and karst cave in sand karst stratum | |
Wang et al. | Development characteristics and finite element simulation of fractures in tight oil sandstone reservoirs of Yanchang Formation in western Ordos Basin | |
CN112965136A (en) | Multi-stage advanced detection method for water-rich karst tunnel | |
Peacock et al. | Evaluation of the controls on fracturing in reservoir rocks | |
Wakabayashi et al. | Geological characterization of melanges for practitioners | |
Zhang et al. | Simulation of tectonic stress field and prediction of tectonic fracture distribution in Longmaxi Formation in Lintanchang area of eastern Sichuan Basin | |
Kushwaha et al. | Stability evaluation of old and unapproachable underground mine workings below surface structures | |
Chang et al. | Three-dimensional analysis of a deep-seated landslide in central Taiwan | |
Zhang | In situ stress inversion and distribution characteristics of tunnel based on numerical simulation and neural network technology | |
Meng et al. | Three‐Dimensional Numerical Modeling and Roof Deformation Analysis of Yuanjue Cave Based on Point Cloud Data | |
Wilson et al. | Fracture model of the Upper Freeport coal: Marshall County West Virginia pilot ECBMR and CO2 sequestration site | |
Umar et al. | Rock mass characterization and conceptual modeling of the printzsköld orebody of the malmberget mine, Sweden | |
CN111538108A (en) | Original geostress field-based mine rock joint structure simulation method | |
Xing et al. | Rock Mass Stability Around Underground Excavations in a Mine: A Case Study | |
Pierce et al. | Synthetic rock mass applications in mass mining | |
Danielsen et al. | The applicability of geoelectrical methods in pre-investigation for construction in rock |
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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20200814 |