56 工程科学学报，第43卷，第1期 method and thermal cracking of rock.Chin Rock Mech Eng framework for subsurface flow and transport simulations.Procedia 2000,19(4):412 Comput Sci,2015,51:1098 (周克群，楚泽涵，张元中，等.岩石热开裂与检测方法研究岩 [38]Cao W J,Huang W B,Jiang F M.A thermal-hydraulic-mechanical 石力学与工程学报，2000,19(4)：412) fully coupled model for heat extraction in enhanced geothermal [23]Yang J,Kang Y L,Li C L,et al.Study on the time-space systems /Proceedings Word Geothermal Congress.Melbourne, distribution of liquid phase around wellbore in the fractured tight 2015:1 sandstone reservoir.Oil Drill Prod Technol,2010,32(1):57 [39]Vecchiarelli A.Application of the 3-D Hydro-Mechanical Model (杨建，康毅力，李朝林，等.裂缝性致密砂岩储层井周液相时空 GEOFRAC in Enhanced Geothermal Systems [Dissertation]. 分布规律.石油钻采工艺，2010,32(1)：57) Massachusetts Institute of Technology,2013 [24]Dong JJ,Hsu J Y,Wu W J.et al.Stress-dependence of the [40]Rutqvist J.Status of the TOUGH-FLAC simulator and recent permeability and porosity of sandstone and shale from TCDP applications related to coupled fluid flow and crustal deformations Hole-A.Int J Rock Mech Min Sci.2010,47(7):1141 Comput Geosci,2011,37(6):739 [25]Lysmer J,Kulemeyer R L.Finite dynamic models for infinite [41]Blair S C,Carlson S R,Lee K,et al.Coupled THM simulations of media.J Eng Mech Div,1969,95:759 the drift scale test at Yucca Mountain /5th American Rock [26]Jaswon M A,Poter A R.An integral equation solution of the Mechanics Symposium (NARMS)and 17th Tunneling Association torsion problem.Proc R Soc London Ser A,1963,273(1353):237 of Canada (TAC)Conference.Toronto,2002 [27]Cundall PA.A computer model for simulating progressive,large [42]Clauser C,Bartels J.Numerical Simulation of Reactive Flow in scale movements in rocky block systems.Proc Symp Int Soc Rock Hot Aquifers:SHEMAT and Processing SHEMAT.Springer Mech,1971,2:129 Science Business Media,2003 [28]Shi G H.Discontinuous Deformation Analysis for Jointed Rock [43]Xu T F,Sonnenthal E,Spycher N,et al.TOUGHREACT-a Masses and Other Dynamics of Block System [Dissertation]. simulation program for non-isothermal multiphase reactive Berkeley:University of California,1988 geochemical transport in variably saturated geologic media: [29]Yan B,Kouame K J,Lv W,et al.Influence of new hydrophobic applications to geothermal injectivity and CO2 geological agent on the mechanical properties of modified cemented paste sequestration.Comput Geosci,2006,32(2):145 backfill.J Mater Res Technol,2019,8(6):5716 [44]Settgast RR,Johnson S,Fu P,et al.Simulation of hydraulic [30]Liu Q S,Liu X W.Research on critical problem for fracture fracture networks in three dimensions utilizing massively parallel network propagation and evolution with multifield coupling of computing resources /Unconventional Resources Technology fractured rock mass.Rock Soil Mech,2014,35(2):305 Conference.Denver,Colorado,2014:1730 (刘泉声，刘学伟.多场耦合作用下岩体裂隙扩展演化关键问题 [45]Wong Z Y,Horne R,Voskov D.A geothermal reservoir simulator 研究.岩土力学，2014,35(2)：305) with AD-GPRS II Proceedings World Geothermal Congress [31]Steefel C I,Appelo CA,Arora B.et al.Reactive transport codes Melbourne,2015:1 for subsurface environmental simulation.Comput Geosci,2015, [46]Lie K A,Krogstad S,Ligaarde I S,et al.Open-source MATLAB 19(3):445 implementation of consistent discretisations on complex grids. [32]White M D,Phillips B R.Code comparison study fosters Comput Geosci,2012,16(2):297 confidence in the numerical simulation of enhanced geothermal [47]Diersch H J G.FEFLOW:Finite Eement Modeling of Flow,Mass systems ll Fortieth Workshop on Geothermal Reservoir and Heat Transport in Porous and fractured Media.Springer Engineering.Stanford,2015:SGP-TR-204 Science Business Media,2013 [33]Pimia P,Duhaime F,Ethier Y,et al.ICY:an interface between [48]Yang B,Xu T F,Li F Y,et al.Numerical simulation on impact of COMSOL multiphysics and discrete element code YADE for the water-rock interaction on reservoir permeability:a case study of modelling of porous media.Comput Geosci,2018,123:38 upper Paleozoic sandstone reservoirs in Northeastern Ordos Basin. [34]Xia Y D,Plummer M,Mattson E,et al.Design,modeling,and J Jilin Univ Earth Sci,2019,49(2):526 evaluation of a doublet heat extraction model in enhanced (杨冰，许天福，李凤显，等.水-岩作用对储层渗透性影响的数值 geothermal systems.Renewable Energy,2017,105:232 模拟研究一一以鄂尔多斯盆地东北部上古生界砂岩储层为例， [35]Pandey S N,Chaudhuri A,Kelkar S.A coupled thermo-hydro- 吉林大学学报（地球科学版），201949(2)：526) mechanical modeling of fracture aperture alteration and reservoir [49]Yan B,Guo Q,Ren F,et al.Modified Nishihara model and deformation during heat extraction from a geothermal reservoir. experimental verification of deep rock mass under the water-rock Geothermics,2017,65:17 interaction.Int J Rock Mech Min Sci,2020.128:104250 [36]Kolditz O,Shao H,Wang W Q,et al.Thermo-Hydro-Mechanical [50]Wang X.Deformation characteristics of deep roadway under the Chemical Processes in Fractured Porous Media:Modelling and influence of temperature and groundwater.Coal Sci Technol Mag Benchmarking.Berlin:Springer,2016 2018(1):22 [37]Scheibe T D,Yang X F,Chen X Y,et al.A hybrid multiscale (王晓.深部巷道在温度和地下水彩响下的变形特征.煤炭科技，method  and  thermal  cracking  of  rock. Chin J Rock Mech Eng, 2000, 19（4）: 412 （周克群, 楚泽涵, 张元中, 等. 岩石热开裂与检测方法研究. 岩 石力学与工程学报, 2000, 19（4）：412） Yang  J,  Kang  Y  L,  Li  C  L,  et  al.  Study  on  the  time-space distribution of liquid phase around wellbore in the fractured tight sandstone reservoir. Oil Drill Prod Technol, 2010, 32（1）: 57 （杨建, 康毅力, 李朝林, 等. 裂缝性致密砂岩储层井周液相时空 分布规律. 石油钻采工艺, 2010, 32（1）：57） [23] Dong  J  J,  Hsu  J  Y,  Wu  W  J,  et  al.  Stress-dependence  of  the permeability  and  porosity  of  sandstone  and  shale  from  TCDP Hole-A. Int J Rock Mech Min Sci, 2010, 47（7）: 1141 [24] Lysmer  J,  Kulemeyer  R  L.  Finite  dynamic  models  for  infinite media. J Eng Mech Div, 1969, 95: 759 [25] Jaswon  M  A,  Poter  A  R.  An  integral  equation  solution  of  the torsion problem. Proc R Soc London Ser A, 1963, 273（1353）: 237 [26] Cundall P A. A computer model for simulating progressive, large scale movements in rocky block systems. Proc Symp Int Soc Rock Mech, 1971, 2: 129 [27] Shi  G  H. Discontinuous Deformation Analysis for Jointed Rock Masses and Other Dynamics of Block System [Dissertation]. Berkeley: University of California, 1988 [28] Yan B, Kouame K J, Lv W, et al. Influence of new hydrophobic agent  on  the  mechanical  properties  of  modified  cemented  paste backfill. J Mater Res Technol, 2019, 8（6）: 5716 [29] Liu  Q  S,  Liu  X  W.  Research  on  critical  problem  for  fracture network  propagation  and  evolution  with  multifield  coupling  of fractured rock mass. Rock Soil Mech, 2014, 35（2）: 305 （刘泉声, 刘学伟. 多场耦合作用下岩体裂隙扩展演化关键问题 研究. 岩土力学, 2014, 35（2）：305） [30] Steefel C I, Appelo C A, Arora B, et al. Reactive transport codes for  subsurface  environmental  simulation. Comput Geosci,  2015, 19（3）: 445 [31] White  M  D,  Phillips  B  R.  Code  comparison  study  fosters confidence  in  the  numerical  simulation  of  enhanced  geothermal systems  // Fortieth Workshop on Geothermal Reservoir Engineering. Stanford, 2015: SGP-TR-204 [32] Pirnia  P,  Duhaime  F,  Ethier  Y,  et  al.  ICY:  an  interface  between COMSOL multiphysics and discrete element code YADE for the modelling of porous media. Comput Geosci, 2018, 123: 38 [33] Xia  Y  D,  Plummer  M,  Mattson  E,  et  al.  Design,  modeling,  and evaluation  of  a  doublet  heat  extraction  model  in  enhanced geothermal systems. Renewable Energy, 2017, 105: 232 [34] Pandey  S  N,  Chaudhuri  A,  Kelkar  S.  A  coupled  thermo-hydro￾mechanical  modeling  of  fracture  aperture  alteration  and  reservoir deformation  during  heat  extraction  from  a  geothermal  reservoir. Geothermics, 2017, 65: 17 [35] Kolditz O, Shao H, Wang W Q, et al. Thermo-Hydro-Mechanical Chemical Processes in Fractured Porous Media: Modelling and Benchmarking. Berlin: Springer, 2016 [36] [37] Scheibe  T  D,  Yang  X  F,  Chen  X  Y,  et  al.  A  hybrid  multiscale framework for subsurface flow and transport simulations. Procedia Comput Sci, 2015, 51: 1098 Cao W J, Huang W B, Jiang F M. A thermal-hydraulic-mechanical fully  coupled  model  for  heat  extraction  in  enhanced  geothermal systems  // Proceedings Word Geothermal Congress.  Melbourne, 2015: 1 [38] Vecchiarelli  A. Application of the 3-D Hydro-Mechanical Model GEOFRAC in Enhanced Geothermal Systems [Dissertation]. Massachusetts Institute of Technology, 2013 [39] Rutqvist  J.  Status  of  the  TOUGH-FLAC  simulator  and  recent applications related to coupled fluid flow and crustal deformations. Comput Geosci, 2011, 37（6）: 739 [40] Blair S C, Carlson S R, Lee K, et al. Coupled THM simulations of the  drift  scale  test  at  Yucca  Mountain  //  5th American Rock Mechanics Symposium (NARMS) and 17th Tunneling Association of Canada (TAC) Conference. Toronto, 2002 [41] Clauser  C,  Bartels  J. Numerical Simulation of Reactive Flow in Hot Aquifers: SHEMAT and Processing SHEMAT.  Springer Science & Business Media, 2003 [42] Xu  T  F,  Sonnenthal  E,  Spycher  N,  et  al.  TOUGHREACT —a simulation  program  for  non-isothermal  multiphase  reactive geochemical  transport  in  variably  saturated  geologic  media: applications  to  geothermal  injectivity  and  CO2 geological sequestration. Comput Geosci, 2006, 32（2）: 145 [43] Settgast  R  R,  Johnson  S,  Fu  P,  et  al.  Simulation  of  hydraulic fracture networks in three dimensions utilizing massively parallel computing  resources  // Unconventional Resources Technology Conference. Denver, Colorado, 2014: 1730 [44] Wong Z Y, Horne R, Voskov D. A geothermal reservoir simulator with  AD-GPRS  // Proceedings World Geothermal Congress. Melbourne, 2015: 1 [45] Lie K A, Krogstad S, Ligaarde I S, et al. Open-source MATLAB implementation  of  consistent  discretisations  on  complex  grids. Comput Geosci, 2012, 16（2）: 297 [46] Diersch H J G. FEFLOW: Finite Eement Modeling of Flow, Mass and Heat Transport in Porous and Fractured Media.  Springer Science & Business Media, 2013 [47] Yang B, Xu T F, Li F Y, et al. Numerical simulation on impact of water-rock  interaction  on  reservoir  permeability:  a  case  study  of upper Paleozoic sandstone reservoirs in Northeastern Ordos Basin. J Jilin Univ Earth Sci, 2019, 49（2）: 526 （杨冰, 许天福, 李凤昱, 等. 水-岩作用对储层渗透性影响的数值 模拟研究——以鄂尔多斯盆地东北部上古生界砂岩储层为例. 吉林大学学报(地球科学版), 2019, 49（2）：526） [48] Yan  B,  Guo  Q,  Ren  F,  et  al.  Modified  Nishihara  model  and experimental verification of deep rock mass under the water-rock interaction. Int J Rock Mech Min Sci, 2020, 128: 104250 [49] Wang  X.  Deformation  characteristics  of  deep  roadway  under  the influence of temperature and groundwater. Coal Sci Technol Mag, 2018（1）: 22 （王晓. 深部巷道在温度和地下水影响下的变形特征. 煤炭科技, [50] · 56 · 工程科学学报，第 43 卷，第 1 期