CN107526920A - Method for analyzing corresponding relation between new surface subsidence dynamic process and overlying strata abscission layer development - Google Patents

Method for analyzing corresponding relation between new surface subsidence dynamic process and overlying strata abscission layer development Download PDF

Info

Publication number
CN107526920A
CN107526920A CN201710590471.6A CN201710590471A CN107526920A CN 107526920 A CN107526920 A CN 107526920A CN 201710590471 A CN201710590471 A CN 201710590471A CN 107526920 A CN107526920 A CN 107526920A
Authority
CN
China
Prior art keywords
subsidence
surface subsidence
corresponding relation
development
dynamic
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
Application number
CN201710590471.6A
Other languages
Chinese (zh)
Other versions
CN107526920B (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.)
Qingdao University of Technology
Original Assignee
Qingdao University of Technology
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 Qingdao University of Technology filed Critical Qingdao University of Technology
Priority to CN201710590471.6A priority Critical patent/CN107526920B/en
Publication of CN107526920A publication Critical patent/CN107526920A/en
Application granted granted Critical
Publication of CN107526920B publication Critical patent/CN107526920B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16ZINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
    • G16Z99/00Subject matter not provided for in other main groups of this subclass

Landscapes

  • Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)

Abstract

The invention discloses a method for analyzing the corresponding relation between a new surface subsidence dynamic process and overburden abscission layer development, which comprises the steps of measuring a sticking measuring point on a model overburden section to obtain the time process of abscission layer from cracking to expanding to closing and quantitative development indexes of the length and the width of the abscission layer; arranging a subsidence measuring point on the model ground surface, and measuring the whole process data of the ground surface subsidence gradually developing from zero to the maximum value at the same time with the separation layer; and obtaining a time schedule by a model test, displaying the dynamic growth data of the surface subsidence and the dynamic development data of the overburden bed separation, statistically analyzing the corresponding relation of the dynamic growth data and the dynamic development data of the overburden bed separation, finding the corresponding relation of the initial period, the active period and the decline period of the surface subsidence and the initiation, development and closing time points of the bed separation, and newly proposing and dividing bed separation time nodes and each development stage. The characteristic of breaking the overburden separation layer can be more clearly reflected by indirectly obtaining data when the physical model and the digital model are used.

Description

A kind of new surface subsidence dynamic process and overlying strata separation layer development corresponding relation analysis method
Technical field
The present invention relates to surface subsidence dynamic process and overlying strata separation layer technical field, specially a kind of new surface subsidence dynamic Process and overlying strata separation layer development corresponding relation analysis method.
Background technology
During coal mine underground is exploited, with the passage of underground mining working face, the spontaneous combustion stress of overlying strata body is destroyed Poised state, due to the effect of the secondary stress field of redistribution, the stress inside rock stratum also divides goaf overlying rock again Cloth is to reach new balance.Fracture, inbreak, bending and sinking will occur for overlying rock, and this motion change is from bottom to top progressively Development, forms subsidence, overlying rock apart from the difference of mined out zone position, thickness etc., and is caused due to lithology to earth's surface The motion of rock stratum is not that Complete Synchronization is consistent, and so as to displacement difference of the overlying strata motion on vertical occur, it is empty to form absciss layer Between.
The Production conditions of overlying strata separation layer:The absciss layer of rock stratum refers to generate Normal Displacement between two adjacent rock stratum, and it is produced Condition depends primarily on, the bending rigidity of lower rock stratum, when the bending rigidity of upper rock stratum is more than the next rock stratum, then can produce Absciss layer, and the bending rigidity of rock stratum depends on the bullet mould and thickness of rock stratum.
The content of the invention
It is an object of the invention to provide a kind of new surface subsidence dynamic process and overlying strata separation layer development corresponding relation analysis Method.
To achieve the above object, the present invention provides following technical scheme:The measurement that measuring point is pasted in model overlying strata section obtains Absciss layer by crack initiation, expand to the time course of closure and its quantization development indicators of absciss layer length and width;
Sinking measuring point is installed in model earth's surface, surface subsidence is measured in the same time with absciss layer and progressively develops into maximum from zero Overall process data;
Time schedule is drawn by model test, exhibition row surface subsidence dynamic increases data and overlying strata separation layer active development number According to, the corresponding relation of both statistical analyses, initial period, active period and the decline phase of surface subsidence and the crack initiation of absciss layer, hair are found Exhibition and the corresponding relation of closing time point, propose from new and divide absciss layer timing node and develop each stage.
A kind of above-mentioned new surface subsidence dynamic process and overlying strata separation layer development corresponding relation analysis method, are cutd open with model overlying strata The measurement that measuring point is pasted in face is analyzed, and overlying strata main key stratum is when it is 30-40m that working face, which adopts hole depth, surface subsidence speed Most fast, its subsidence velocity reaches 7.5mm/d, it can be seen that, earth's surface maximum subsidence velocity is synchronous with key stratum maximum subsidence velocity, The first breaking of key stratum result in the quick sinking of earth's surface,
Calculation formula is:
DW (t)/dt=c [W (t)/(1-1-e-ct)-Wk(1-1-e-ct)]
W (t) is a certain moment Dynamic Subsidence value
WkFor earth's surface final subsidence value
1-e-ctTo obtain variable quantity in certain time
DW (t)/dt is surface subsidence speed, as W (t)/(1-1-e-ct)-Wk(1-1-e-ct) value be more than 0 when surface subsidence Speed can be accelerated, as W (t)/(1-1-e-ct)-Wk(1-1-e-ct) value be less than 0 when earth's surface decrease speed can reduce, but due to by The condition limitation of plagioclase and length into model so that the numerical value and actual numerical value drawn has certain deviation.
Preferably, a kind of new surface subsidence dynamic process and overlying strata separation layer development corresponding relation analysis method, pass through The parameter setting at joint between key stratum block is obtained very big, make not producing motion between block, broken with this to simulate key stratum rock mass Moved before disconnected, when adopting width and reaching key stratum broken mechanics, reduce block joint parameters, moved after simulation key stratum is disrumpent feelings, overlying strata Key stratum obtains the disrumpent feelings surface subsidence that will influence and obtains dynamic process,
Calculation formula is:
W (t)=W0ω(1-ae-kt)
W (t) is a certain moment Dynamic Subsidence value
W0The largest amount of subsidence of earth's surface at the end of for the decline phase
ΔωTerminate the remaining sinking of rear earth's surface for the decline phase
Due to W0For solid value, the amount that W (t) sinks is relevant with a, k and t, works as 1-ae-ktDuring more than 0, W (t) sinkings can add Greatly, 1-ae is worked as-ktDuring less than 0, W (t) sinkings can slow down, and the disrumpent feelings of key strata of covering rock will cause the obvious of surface subsidence speed Increase is in not cyclically-varying.
Preferably, a kind of new surface subsidence dynamic process and overlying strata separation layer development corresponding relation analysis method, are carried out Simulation is excavated in coal seam according to irregular depth, after each excavate model calculation to starting again in next step after steady Excavate, when depth is advanced into 150m, key stratum occurs first breaking, first time periodic failure and second of cycle and broken respectively Disconnected, before key stratum is disrumpent feelings, surface subsidence increment is smaller, and after key stratum is disrumpent feelings, surface subsidence increment can increase, and earth's surface Moving influence angle and moving influence border are not unalterable,
Calculation formula is:
W (x, y, t)=Wm(1-x2/a2-y2/b2)(1-e-ct)k
W (x, y, t) is subsidence basin any point deflection at any time
1-e-ctTo obtain variable quantity in certain time
(1-x2/a2-y2/b2)(1-e-ct)kWhen product is positive number, W (x, y, t) deflection can increase, (1-x2/a2-y2/b2) (1-e-ct)kWhen product is negative, W (x, y, t) deflection be able to can slow down, and with a and b increasing, surface subsidence influences border will It is obvious to expand.
Compared with prior art, beneficial effects of the present invention are:Disrumpent feelings due to key strata of covering rock will cause surface subsidence Speed significantly increase, therefore can will ensure the principle that be designed as coal mining underbuildings of the not disrumpent feelings unstability of key strata of covering rock, be Guarantee coal mining underbuildings both has a preferable economic benefit, while also ensures that above ground structure does not suffer damage, and research is closed Key is to determine reasonably to subtract heavy production technique and parameter according to sand coated iron mold under actual conditions and key stratum feature.
Brief description of the drawings
Fig. 1 develops settling velocity under corresponding relation analysis method for a kind of new surface subsidence dynamic process of the present invention with overlying strata separation layer Degree and ground distance meander line structure schematic diagram;
Fig. 2 develops corresponding relation analysis method rock stratum power for a kind of new surface subsidence dynamic process of the present invention and overlying strata separation layer Learn parameter table schematic diagram;
Fig. 3 develops in corresponding relation analysis method for a kind of new surface subsidence dynamic process of the present invention with overlying strata separation layer adopts width Angle form schematic diagram is influenceed with surface movement.
Embodiment
The present invention provides a kind of technical scheme:The measurement that measuring point is pasted in model overlying strata section obtains absciss layer by crack initiation, expansion Open up the time course of closure and its quantization development indicators of absciss layer length and width;
Sinking measuring point is installed in model earth's surface, surface subsidence is measured in the same time with absciss layer and progressively develops into maximum from zero Overall process data;
Time schedule is drawn by model test, exhibition row surface subsidence dynamic increases data and overlying strata separation layer active development number According to, the corresponding relation of both statistical analyses, initial period, active period and the decline phase of surface subsidence and the crack initiation of absciss layer, hair are found Exhibition and the corresponding relation of closing time point, propose from new and divide absciss layer timing node and develop each stage.
A kind of new surface subsidence dynamic process and overlying strata separation layer development corresponding relation analysis method, are sticked with model overlying strata section The measurement of patch measuring point is analyzed, and for overlying strata main key stratum when it is 30-40m that working face, which adopts hole depth, surface subsidence is fastest, Its subsidence velocity reaches 7.5mm/d, it can be seen that, earth's surface maximum subsidence velocity is synchronous with key stratum maximum subsidence velocity, crucial The first breaking of layer result in the quick sinking of earth's surface,
Calculation formula is:
DW (t)/dt=c [W (t)/(1-1-e-ct)-Wk(1-1-e-ct)]
W (t) is a certain moment Dynamic Subsidence value
WkFor earth's surface final subsidence value
1-e-ctTo obtain variable quantity in certain time
DW (t)/dt is surface subsidence speed, as W (t)/(1-1-e-ct)-Wk(1-1-e-ct) value be more than 0 when surface subsidence Speed can be accelerated, as W (t)/(1-1-e-ct)-Wk(1-1-e-ct) value be less than 0 when earth's surface decrease speed can reduce, but due to by The condition limitation of plagioclase and length into model so that the numerical value and actual numerical value drawn has certain deviation.
A kind of new surface subsidence dynamic process and overlying strata separation layer development corresponding relation analysis method, by by key stratum block Between the parameter setting at joint obtain very big, make not producing motion between block, with this come simulate key stratum rock mass it is disrumpent feelings before move, when adopting When width reaches key stratum broken mechanics, reduce block joint parameters, simulation key stratum it is disrumpent feelings after move, key strata of covering rock obtain it is disrumpent feelings will Influence surface subsidence and obtain dynamic process,
Calculation formula is:
W (t)=W0ω(1-ae-kt)
W (t) is a certain moment Dynamic Subsidence value
W0The largest amount of subsidence of earth's surface at the end of for the decline phase
ΔωTerminate the remaining sinking of rear earth's surface for the decline phase
Due to W0For solid value, the amount that W (t) sinks is relevant with a, k and t, works as 1-ae-ktDuring more than 0, W (t) sinkings can add Greatly, 1-ae is worked as-ktDuring less than 0, W (t) sinkings can slow down, and the disrumpent feelings of key strata of covering rock will cause the obvious of surface subsidence speed Increase is in not cyclically-varying.
A kind of new surface subsidence dynamic process and overlying strata separation layer development corresponding relation analysis method, carry out in simulation coal seam by Excavated according to irregular depth, after each excavate model calculation to starting to excavate in next step again after steady, work as depth When being advanced into 150m, it is disrumpent feelings that first breaking, first time periodic failure and second of periodic failure, key stratum occur respectively for key stratum Before, surface subsidence increment is smaller, and after key stratum is disrumpent feelings, and surface subsidence increment can increase, and surface movement influence angle and Moving influence border is not unalterable,
Calculation formula is:
W (x, y, t)=Wm(1-x2/a2-y2/b2)(1-e-ct)k
W (x, y, t) is subsidence basin any point deflection at any time
1-e-ctTo obtain variable quantity in certain time
(1-x2/a2-y2/b2)(1-e-ct)kWhen product is positive number, W (x, y, t) deflection can increase, (1-x2/a2-y2/b2) (1-e-ct)kWhen product is negative, W (x, y, t) deflection be able to can slow down, and with a and b increasing, surface subsidence influences border will It is obvious to expand.
On analog simulation test model, the measurement that measuring point is pasted by model overlying strata section obtains absciss layer by crack initiation, expansion Open up the time course of closure and its quantization development indicators of absciss layer length and width;Meanwhile sinking measuring point is installed in model earth's surface, Surface subsidence is measured in the same time with absciss layer from the zero overall process data for progressively developing into maximum.Row ground is opened up according to time schedule Table, which sinks, dynamically increases data and overlying strata separation layer active development data, the corresponding relation of both statistical analyses, finds surface subsidence Initial period, active period and decline phase and the crack initiation of absciss layer, the corresponding relation of development and closing time point, count each period absciss layer To the contribution rate of surface subsidence, propose from new and divide absciss layer timing node and develop each stage;According to experiment similarity criterion Surface subsidence on the spot is converted into developing corresponding relation with overlying strata separation layer, provides overlying strata separation layer grouting control surface subsidence most Good time standard.
The preferred embodiments of the present invention are the foregoing is only, are not intended to limit the invention, although with reference to foregoing reality Apply example the present invention is described in detail, for those skilled in the art, it still can be to foregoing each implementation Technical scheme described in example is modified, or carries out equivalent substitution to which part technical characteristic.All essences in the present invention God any modification, equivalent substitution and improvements made etc., should be included in the scope of the protection with principle.

Claims (4)

1. a kind of new surface subsidence dynamic process and overlying strata separation layer development corresponding relation analysis method, it is characterised in that:
Model overlying strata section paste measuring point measurement obtain absciss layer by crack initiation, expand to the time course and its absciss layer of closure The quantization development indicators of length and width;
Sinking measuring point is installed in model earth's surface, surface subsidence is measured in the same time with absciss layer and progressively develops into the complete of maximum from zero Process data;
Time schedule is drawn by model test, exhibition row surface subsidence dynamic increases data and overlying strata separation layer active development data, The corresponding relation of both statistical analyses, find initial period, active period and the decline phase of surface subsidence and the crack initiation of absciss layer, development and The corresponding relation of closing time point, propose from new and divide absciss layer timing node and develop each stage.
2. a kind of new surface subsidence dynamic process according to claim 1 and overlying strata separation layer development corresponding relation analysis side Method, it is characterised in that:
The measurement that measuring point is pasted with model overlying strata section is analyzed, and it is 30-40m that overlying strata main key stratum adopts hole depth in working face When, surface subsidence is fastest, and its subsidence velocity reaches 7.5mm/d, it can be seen that, earth's surface maximum subsidence velocity and key stratum are most Rapid sink rate is synchronous, and the first breaking of key stratum result in the quick sinking of earth's surface,
Calculation formula is:
DW (t)/dt=c [W (t)/(1-1-e-ct)-Wk(1-1-e-ct)]
W (t) is a certain moment Dynamic Subsidence value
WkFor earth's surface final subsidence value
1-e-ctTo obtain variable quantity in certain time
DW (t)/dt is surface subsidence speed, as W (t)/(1-1-e-ct)-Wk(1-1-e-ct) value be more than 0 when surface subsidence speed meeting Accelerate, as W (t)/(1-1-e-ct)-Wk(1-1-e-ct) value when being less than 0 earth's surface decrease speed can reduce, but due to by model In plagioclase and length condition limitation so that the numerical value and actual numerical value drawn has certain deviation.
3. a kind of new surface subsidence dynamic process according to claim 1 and overlying strata separation layer development corresponding relation analysis side Method, it is characterised in that:
It is very big by the way that the parameter setting at joint between key stratum block is obtained, make not produce motion between block, key is simulated with this Moved before layer rock mass is disrumpent feelings, when adopting width and reaching key stratum broken mechanics, reduce block joint parameters, transported after simulation key stratum is disrumpent feelings Dynamic, key strata of covering rock obtains the disrumpent feelings surface subsidence that will influence and obtains dynamic process,
Calculation formula is:
W (t)=W0ω(1-ae-kt)
W (t) is a certain moment Dynamic Subsidence value
W0The largest amount of subsidence of earth's surface at the end of for the decline phase
ΔωTerminate the remaining sinking of rear earth's surface for the decline phase
Due to W0For solid value, the amount that W (t) sinks is relevant with a, k and t, works as 1-ae-ktDuring more than 0, W (t) sinkings can increase, when 1-ae-ktDuring less than 0, W (t) sinkings can slow down, and the disrumpent feelings of key strata of covering rock will cause significantly increasing simultaneously for surface subsidence speed Non- is in cyclically-varying.
4. a kind of new surface subsidence dynamic process according to claim 1 and overlying strata separation layer development corresponding relation analysis side Method, it is characterised in that:
Excavated in simulation coal seam according to irregular depth, after each excavate model calculation to opening again after steady Begin to excavate in next step, when depth is advanced into 150m, first breaking, first time periodic failure and second occur respectively for key stratum Periodic failure, before key stratum is disrumpent feelings, surface subsidence increment is smaller, and after key stratum is disrumpent feelings, surface subsidence increment can increase, and And surface movement influences angle and moving influence border is not unalterable,
Calculation formula is:
W (x, y, t)=Wm(1-x2/a2-y2/b2)(1-e-ct)k
W (x, y, t) is subsidence basin any point deflection at any time
1-e-ctTo obtain variable quantity in certain time
(1-x2/a2-y2/b2)(1-e-ct)kWhen product is positive number, W (x, y, t) deflection can increase, (1-x2/a2-y2/b2)(1-e-ct)kWhen product is negative, W (x, y, t) deflection be able to can slow down, and with a and b increasing, surface subsidence influences border will be obvious Expand.
CN201710590471.6A 2017-07-19 2017-07-19 Method for analyzing corresponding relation between new surface subsidence dynamic process and overlying strata abscission layer development Active CN107526920B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710590471.6A CN107526920B (en) 2017-07-19 2017-07-19 Method for analyzing corresponding relation between new surface subsidence dynamic process and overlying strata abscission layer development

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710590471.6A CN107526920B (en) 2017-07-19 2017-07-19 Method for analyzing corresponding relation between new surface subsidence dynamic process and overlying strata abscission layer development

Publications (2)

Publication Number Publication Date
CN107526920A true CN107526920A (en) 2017-12-29
CN107526920B CN107526920B (en) 2022-01-14

Family

ID=60748967

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710590471.6A Active CN107526920B (en) 2017-07-19 2017-07-19 Method for analyzing corresponding relation between new surface subsidence dynamic process and overlying strata abscission layer development

Country Status (1)

Country Link
CN (1) CN107526920B (en)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103216236A (en) * 2013-03-13 2013-07-24 中国矿业大学(北京) Method for reducing subsidence of stope overburden bed separation by grouting
CN103558646A (en) * 2013-10-31 2014-02-05 中国矿业大学 Method of detecting dynamic developmental situation of overburden mining-induce fissure in radon gas earth surface
CN105067294A (en) * 2015-07-17 2015-11-18 青岛理工大学 Mining overburden bed separation different-speed growth test monitoring system and analysis method
CN106772678A (en) * 2016-12-16 2017-05-31 安徽理工大学 A kind of wellhole many reference amounts method for surveying of rock deformation destructive characteristics
CN106884657A (en) * 2017-03-31 2017-06-23 青岛理工大学 Method for displaying surface subsidence dynamic process based on overburden bed separation dynamic development

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103216236A (en) * 2013-03-13 2013-07-24 中国矿业大学(北京) Method for reducing subsidence of stope overburden bed separation by grouting
CN103558646A (en) * 2013-10-31 2014-02-05 中国矿业大学 Method of detecting dynamic developmental situation of overburden mining-induce fissure in radon gas earth surface
CN105067294A (en) * 2015-07-17 2015-11-18 青岛理工大学 Mining overburden bed separation different-speed growth test monitoring system and analysis method
CN106772678A (en) * 2016-12-16 2017-05-31 安徽理工大学 A kind of wellhole many reference amounts method for surveying of rock deformation destructive characteristics
CN106884657A (en) * 2017-03-31 2017-06-23 青岛理工大学 Method for displaying surface subsidence dynamic process based on overburden bed separation dynamic development

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
胡晓阳: "基于异速生长理论的采动覆岩离层时空分布规律与沉陷模型研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅰ辑》 *

Also Published As

Publication number Publication date
CN107526920B (en) 2022-01-14

Similar Documents

Publication Publication Date Title
CN102678118B (en) Method for determining height of caving zone of stope
CN106884677B (en) Tight roof super high seam exploits strong mine pressure prediction pre-control method
Deng et al. Strata behavior in extra-thick coal seam mining with upward slicing backfilling technology
Liu et al. Voussoir beam model for lower strong roof strata movement in longwall mining–Case study
CN106354972A (en) Construction technology optimization method for tunnel in nearly horizontal soft and hard interbedded surrounding rock
Toraño et al. FEM modeling of roadways driven in a fractured rock mass under a longwall influence
Liu et al. Characteristics analysis of roof overburden fracture in thick coal seam in deep mining and engineering application of super high water material in backfill mining
CN102733812A (en) Small interval tunnel side wall expanded excavation tunnel construction method for turnout section of river bottom of city
He et al. A method for predicting the water-flowing fractured zone height based on an improved key stratum theory
Liu et al. A prediction model and numerical simulation of the location of the longwall face during the highest possible failure period of gob gas ventholes
CN115903078A (en) Hydraulic fracturing layer position, determination method thereof and determination method of hydraulic fracturing position
Zhang et al. Strata movement and stress evolution when mining two overlapping panels affected by hard stratum
CN110259443A (en) A kind of coal seam borehole wall stability prediction method based on 3DEC discrete element
CN107526920A (en) Method for analyzing corresponding relation between new surface subsidence dynamic process and overlying strata abscission layer development
Gao et al. Understanding the evolution of mining-induced fractures using physical and numerical modeling
Fuenzalida et al. REBOP–FLAC3D hybrid approach to cave modelling
CN108108519A (en) A kind of filling stope mine pressing reduction characterizing method
Ren et al. Strata movement and fracture evolution characteristics in adjacent seam mining
Liu et al. Numerical simulation on instability mechanism of section coal pillar in shallow buried coal seams
CN116822008B (en) Method and system for constructing three-dimensional perspective stope of coal mine
Guayacán-Carrillo et al. Poroelastic analysis of tunnel face advance in low permeability anisotropic ground
Shlash et al. Stress Analysis Around Tunnels During Construction Stages Using Finite Elements Method
CN104360013B (en) Deep mining floor strata lesion ruptures develops and the analysis method of failure law
Lv et al. Breaking of key layers and surface subsidence in the loess mountainous areas
Beer et al. Numerical simulation of tunnel construction in volcanic rocks

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