An underground iron mine in China has been used as a case-study to research the subsidence due to ore extraction and backfilling during open stoping operations. A 3-D discontinuum numerical model was built incorporati...An underground iron mine in China has been used as a case-study to research the subsidence due to ore extraction and backfilling during open stoping operations. A 3-D discontinuum numerical model was built incorporating geologic complexities including faults and interfaces between different lithologies,and the stoping and backfilling sequence adopted from the mine plans. The stoping was carried out in two vertically stacked horizontal layers, with a total of 16 stopes. Large displacements of up to 50 cm were observed along the roof of the stopes, and a maximum surface subsidence of 22.5 cm was observed.Backfilling was found to eliminate subsequent displacements and subsidence. The extraction of the upper orebody was found to influence displacements in the lower orebody. Finally, a subsidence profile was constructed to show the subsidence at all locations along the length of the surface and region of influence on the surface.展开更多
爆破开采是露天矿采选总成本控制的首要环节,数值模拟是进行露天矿爆破开采优化设计及爆破效果分析的有效手段。利用连续-非连续单元方法(continuum-discontinuum element method,CDEM)对露天矿的三维台阶爆破过程进行了模拟,通过朗道...爆破开采是露天矿采选总成本控制的首要环节,数值模拟是进行露天矿爆破开采优化设计及爆破效果分析的有效手段。利用连续-非连续单元方法(continuum-discontinuum element method,CDEM)对露天矿的三维台阶爆破过程进行了模拟,通过朗道爆炸模型实现了爆炸作用力的精确计算,通过弹性-损伤-断裂本构实现了岩体损伤破裂过程的描述,通过半弹簧-目标面及半棱-目标棱的联合接触算法实现了破碎岩块碰撞、飞散及堆积过程的高效模拟。开展了小尺度单自由面爆破过程的数值模拟,计算给出的块度分布曲线、爆破漏斗体积等参数与文献中模型实验的结果基本一致,证明了CDEM及本文所述各类模型在模拟爆炸破岩方面的精确性。以鞍千矿南采区的露天铁矿爆破开采为研究对象,建立了3排21炮孔的三维台阶爆破概化模型,模拟了从炸药起爆、岩体损伤破裂到最后爆堆形成的全过程;计算结果表明,除后缘拉裂槽外,数值计算给出的爆堆形态、顶部鼓起高度等与现场的测试结果基本一致,证明了利用CDEM开展三维露天台阶爆破全过程模拟的可行性。展开更多
The greatest challenges of rigorously modeling coupled hydro-mechanical(HM)processes in fractured geological media at different scales are associated with computational geometry.These challenges include dynamic sheari...The greatest challenges of rigorously modeling coupled hydro-mechanical(HM)processes in fractured geological media at different scales are associated with computational geometry.These challenges include dynamic shearing and opening of intersecting fractures at discrete fracture scales as a result of coupled processes,and contact alteration along rough fracture surfaces that triggers structural and physical changes of fractures at micro-asperity scale.In this paper,these challenges are tackled by developing a comprehensive modeling approach for coupled processes in fractured geological media based on numerical manifold method(NMM)at multiple scales.Based on their distinct geometric features,fractures are categorized into three different scales:dominant fracture,discrete fracture,and discontinuum asperity scales.Here the scale is relative,that of the fracture relative to that of the research interest or domain.Different geometric representations of fractures at different scales are used,and different governing equations and constitutive relationships are applied.For dominant fractures,a finite thickness zone model is developed to treat a fracture as a porous nonlinear domain.Nonlinear fracture mechanical behavior is accurately modeled with an implicit approach based on strain energy.For discrete fractures,a zero-dimensional model was developed for analyzing fluid flow and mechanics in fractures that are geometrically treated as boundaries of the rock matrix.With the zero-dimensional model,these fractures can be modeled with arbitrary orientations and intersections.They can be fluid conduits or seals,and can be open,bonded or sliding.For the discontinuum asperity scale,the geometry of rough fracture surfaces is explicitly represented and contacts involving dynamic alteration of contacts among asperities are rigorously calculated.Using this approach,fracture alteration caused by deformation,re-arrangement and sliding of rough surfaces can be captured.Our comprehensive model is able to handle the computational challenges with accurate representation of intersections and shearing of fractures at the discrete fracture scale and rigorously treats contacts along rough fracture surfaces at the discontinuum asperity scale.With future development of three-dimensional(3D)geometric representation of discrete fracture networks in porous rock and contacts among multi-body systems,this model is promising as a basis of 3D fully coupled analysis of fractures at multiple scales,for advancing understanding and optimizing energy recovery and storage in fractured geological media.展开更多
This paper focuses on the progress in geomechanical modeling associated with carbon dioxide(CO2)geological storage.The detailed review of some geomechanical aspects,including numerical methods,stress analysis,ground d...This paper focuses on the progress in geomechanical modeling associated with carbon dioxide(CO2)geological storage.The detailed review of some geomechanical aspects,including numerical methods,stress analysis,ground deformation,fault reactivation,induced seismicity and crack propagation,is presented.It is indicated that although all the processes involved are not fully understood,integration of all available data,such as ground survey,geological conditions,microseismicity and ground level deformation,has led to many new insights into the rock mechanical response to CO2injection.The review also shows that in geomechanical modeling,continuum modeling methods are predominant compared with discontinuum methods.It is recommended to develop continuum-discontinuum numerical methods since they are more convenient for geomechanical modeling of CO2geological storage,especially for fracture propagation simulation.The Mohr-Coulomb criterion is widely used in prediction of rock mass mechanical behavior.It would be better to use a criterion considering the effect of the intermediate principal stress on rock mechanical behavior,especially for the stability analysis of deeply seated rock engineering.Some challenges related to geomechanical modeling of CO2geological storage are also discussed.展开更多
基金the support from the National Natural Science Foundation of China (No. 51404024)partially funded by the NIOSH of the Centers for Disease Control and Prevention, USA through Contract No.200-201139886 given to the second author
文摘An underground iron mine in China has been used as a case-study to research the subsidence due to ore extraction and backfilling during open stoping operations. A 3-D discontinuum numerical model was built incorporating geologic complexities including faults and interfaces between different lithologies,and the stoping and backfilling sequence adopted from the mine plans. The stoping was carried out in two vertically stacked horizontal layers, with a total of 16 stopes. Large displacements of up to 50 cm were observed along the roof of the stopes, and a maximum surface subsidence of 22.5 cm was observed.Backfilling was found to eliminate subsequent displacements and subsidence. The extraction of the upper orebody was found to influence displacements in the lower orebody. Finally, a subsidence profile was constructed to show the subsidence at all locations along the length of the surface and region of influence on the surface.
文摘爆破开采是露天矿采选总成本控制的首要环节,数值模拟是进行露天矿爆破开采优化设计及爆破效果分析的有效手段。利用连续-非连续单元方法(continuum-discontinuum element method,CDEM)对露天矿的三维台阶爆破过程进行了模拟,通过朗道爆炸模型实现了爆炸作用力的精确计算,通过弹性-损伤-断裂本构实现了岩体损伤破裂过程的描述,通过半弹簧-目标面及半棱-目标棱的联合接触算法实现了破碎岩块碰撞、飞散及堆积过程的高效模拟。开展了小尺度单自由面爆破过程的数值模拟,计算给出的块度分布曲线、爆破漏斗体积等参数与文献中模型实验的结果基本一致,证明了CDEM及本文所述各类模型在模拟爆炸破岩方面的精确性。以鞍千矿南采区的露天铁矿爆破开采为研究对象,建立了3排21炮孔的三维台阶爆破概化模型,模拟了从炸药起爆、岩体损伤破裂到最后爆堆形成的全过程;计算结果表明,除后缘拉裂槽外,数值计算给出的爆堆形态、顶部鼓起高度等与现场的测试结果基本一致,证明了利用CDEM开展三维露天台阶爆破全过程模拟的可行性。
基金supported by Laboratory Directed Research and Development(LDRD)funding from Berkeley Labsupported by Open Fund of the State Key Laboratory of Geomechanics and Geotechnical Engineering,Institute of Rock and Soil Mechanics,Chinese Academy of Sciences(Grant No.Z017004)。
文摘The greatest challenges of rigorously modeling coupled hydro-mechanical(HM)processes in fractured geological media at different scales are associated with computational geometry.These challenges include dynamic shearing and opening of intersecting fractures at discrete fracture scales as a result of coupled processes,and contact alteration along rough fracture surfaces that triggers structural and physical changes of fractures at micro-asperity scale.In this paper,these challenges are tackled by developing a comprehensive modeling approach for coupled processes in fractured geological media based on numerical manifold method(NMM)at multiple scales.Based on their distinct geometric features,fractures are categorized into three different scales:dominant fracture,discrete fracture,and discontinuum asperity scales.Here the scale is relative,that of the fracture relative to that of the research interest or domain.Different geometric representations of fractures at different scales are used,and different governing equations and constitutive relationships are applied.For dominant fractures,a finite thickness zone model is developed to treat a fracture as a porous nonlinear domain.Nonlinear fracture mechanical behavior is accurately modeled with an implicit approach based on strain energy.For discrete fractures,a zero-dimensional model was developed for analyzing fluid flow and mechanics in fractures that are geometrically treated as boundaries of the rock matrix.With the zero-dimensional model,these fractures can be modeled with arbitrary orientations and intersections.They can be fluid conduits or seals,and can be open,bonded or sliding.For the discontinuum asperity scale,the geometry of rough fracture surfaces is explicitly represented and contacts involving dynamic alteration of contacts among asperities are rigorously calculated.Using this approach,fracture alteration caused by deformation,re-arrangement and sliding of rough surfaces can be captured.Our comprehensive model is able to handle the computational challenges with accurate representation of intersections and shearing of fractures at the discrete fracture scale and rigorously treats contacts along rough fracture surfaces at the discontinuum asperity scale.With future development of three-dimensional(3D)geometric representation of discrete fracture networks in porous rock and contacts among multi-body systems,this model is promising as a basis of 3D fully coupled analysis of fractures at multiple scales,for advancing understanding and optimizing energy recovery and storage in fractured geological media.
基金finically supported by the National Natural Science Foundation of China(Grant Nos.41272349 and 51322906)Key Research Program of Frontier Sciences,Chinese Academy of Sciences(Grant No.QYZDB-SSW-DQC029)
文摘This paper focuses on the progress in geomechanical modeling associated with carbon dioxide(CO2)geological storage.The detailed review of some geomechanical aspects,including numerical methods,stress analysis,ground deformation,fault reactivation,induced seismicity and crack propagation,is presented.It is indicated that although all the processes involved are not fully understood,integration of all available data,such as ground survey,geological conditions,microseismicity and ground level deformation,has led to many new insights into the rock mechanical response to CO2injection.The review also shows that in geomechanical modeling,continuum modeling methods are predominant compared with discontinuum methods.It is recommended to develop continuum-discontinuum numerical methods since they are more convenient for geomechanical modeling of CO2geological storage,especially for fracture propagation simulation.The Mohr-Coulomb criterion is widely used in prediction of rock mass mechanical behavior.It would be better to use a criterion considering the effect of the intermediate principal stress on rock mechanical behavior,especially for the stability analysis of deeply seated rock engineering.Some challenges related to geomechanical modeling of CO2geological storage are also discussed.
