基于Ansys、FLAC3D和Fortran的某抽水蓄能电站初始应力场反演回归分析.pdfVIP

基于Ansys 、FLAC3D 和Fortran 的某抽水蓄能电站初 始应力场反演回归分析 1 1 2 岑成汉 ，张燎军 ，左晓霞 1.河海大学水利水电工程学院，江苏南京（210098 ） 2.浙江省水利水电勘测设计院，浙江杭州（310002 ） E-mail ：lbcch@126.com 摘 要：本文结合某抽水蓄能电站工程地质条件以及实测地应力资料，分析了研究区域初始 地应力的影响因素。在各影响因素单独作用下，采用ansys 、FLAC3D 以及Fortran 三种计算 程序结合，对建立的地质概化计算模型进行计算；在实测点地应力值与计算得到的应力值之 间建立多元回归模型，通过多元回归分析，求出最优回归系数。通过实测点的计算应力值与 现场实测值的比较，两者在量值上很接近，表明经过回归得到的地应力场是合理的，从而获 得某抽水蓄能电站初始地应力场较为合理的分布规律，为电站开挖模拟及长期稳定性分析提 供了合理的三维初始地应力场。 关键词：FLAC3D ，地应力场模拟，多元回归，逐步回归分析，反演 1．引言 某抽水蓄能电站主副厂房洞室长158.50m，其中主机段长94.50m，副厂房长20.00m，安 装场长44.00m ；厂房顶拱高程130.80m，厂房总高度49.20m ；吊车梁以下开挖宽度为22.00m， 吊车梁以上宽度为23.60m 。地下厂房区受区域构造影响大，构造较发育，岩层产状扭曲变化 较大，上库进口至上水平段为N75~85°E NW ∠35~45°，下平段为背斜核部，下平段至下游岩 层产状为N70~80°E SE ∠15~35°，从上库进出水口至斜井段为一等斜褶皱段，其中上库进出 口至上平段为上库向斜的南翼，同时也是背斜的北翼，背斜轴与输水系统轴线近直交，背斜 轴部位于厂房顶部。 根据地勘资料，输水线路沿线揭露的断层共62条，主要位于尾水反坡段核上平段，与洞 线近直交，走向以NE方向为主，在此我们主要考虑对电站工程区影响较大的三个断层：F11 ， F521 ，F113 ，受断层核褶皱作用，输水系统部位挤压破碎带和节理近平行主干断层发育，且 分布疏密不均，在岩质较软弱的泥质砂岩夹层内多形成挤压破碎带，在断层之间多组成节理 密集带发育。 由于地质构造复杂，利用ansys 的强大建模功能建立合理的计算模型，利用FLAC3D和 Fortran程序，采用多元回归分析的方法对某抽水蓄能电站的三维初始地应力场进行了分析。 本文结合ansys，FLAC3D ，Fortran对该电站地下厂房区进行分析，得到合理的应力场分布， 使其能满足某抽水蓄能电站开挖模拟和长期稳定性分析的需要。 2 ．地应力实测结果 该电站的实测地应力值采用水压致裂法进行测量，在地应力场回归分析中，以计算坐标 系xyz 中坐标应力分量为基本对象，而地应力实测值所在的坐标系和我们模型的计算坐标系 不一致，需要将实测地应力进行正确的坐标变换以得到计算坐标系下的应力分量。由于计算 坐标系和实测值所用的坐标系中z轴是一致的，所以只需要转换两个坐标系的x ，y轴，于是 我们所需要的两个坐标系的之间的方向余弦为： L cosα ⎧ ⎨ （1） M cos β ⎩ -1- 式中L为计算坐标系中σi与x轴正方向之间的夹角，M为计算坐标系中σi与y轴正方向之 间的夹角 根据每组实测主应力量值及方位角、倾角，可由式（2 ）将主应力转换成计算坐标系下 的应力分量：