基于BPA法计算特高压直流输电线路合成电场_邹岸新.pdf

第4 1卷第5期2023年5月水 电 能 源 科 学W a t e r R e s o u r c e s a n d P o w e rV o l.4 1 N o.5M a y 2 0 2 3D O I:1 0.2 0 0 4 0/j.c n k i.1 0 0 0-7 7 0 9.2 0 2 3.2 0 2 2 1 4 3 1基于B P A法计算特高压直流输电线路合成电场邹岸新1,王守刚1,杨 滔1,李永明2,刘 岩2(1.国网重庆市电力公司超高压分公司,重庆 4 0 0 0 3 9;2.重庆大学输配电装备及系统安全与新技术国家重点实验室,重庆 4 0 0 0 4 4)摘要:为研究特高压直流输电线路地面附近的离子流场问题,应用B P A法求解双极直流输电线路离子流场,结合逐次镜像法计算的标称电场,进而求得合成电场,并利用该方法计算了8 0 0 k V直流输电线路的离子流场与合成电场问题,分析了导线对地高度、极间距、子导线半径对特高压直流输电线路合成电场的影响。此外,还对比了求解合成电场强度的电场线法。结果表明,B P A法求解离子流场准确有效,且计算效率大大提高;提高导线架设高度和增加分裂子导线半径均可改善地面附近电磁环境,而减小导线极间距能够降低地面附近的合成电场强度,但效果不明显。关键词:B P A法;特高压直流输电;合成电场;对地高度;极间距;分裂子导线半径中图分类号:TM 7 5 文献标志码:A 文章编号:1 0 0 0-7 7 0 9(2 0 2 3)0 5-0 1 9 4-0 5收稿日期:2 0 2 2-0 7-1 3,修回日期:2 0 2 2-0 8-1 0作者简介:邹岸新(1 9 8 5-),硕士、高级工程师,研究方向为电网电磁环境与噪声控制技术,E-m a i l:z o u a x i n 1 6 3.c o m1 引言随着电力系统的发展,为满足远距离、大容量输电的需要,特高压直流输电工程在国内得到了快速的发展和推广。在直流输电系统中,当导线表面的电场强度超过起晕场强时,周围空气会发生局部放电、产生电晕,空气中游离的离子会在电场的作用下产生离子流,增强地面附近的电场。因此特高压直流输电引发的电磁环境问题日益成为关注的焦点。直流合成电场强度和离子流密度是直流输电线路的电磁环境评估指标。高压电气设备周围合成电场强度和离子流密度的控制方程均为耦合偏微分方程组1,其解析解较难获得,需对其进行简化后再求解。直流输电线路离子流场的计算方法总体可分为两大类,即以D e u t s c h假设为基础的通量线法和采用网格剖分的数值解法2-3。李永明等4基于D e u t s c h对地面离子流场进 行 解 析 求 解;乔 骥 等5提 出 一 种 消 除D e u t s c h假设的高精度迭代特征线法;江沨等6采用高精度上流有限元法求解离子流场。邦纳维尔电业局(B P A)法以其求解效率高、稳定性好等优点多次在C I G R E被阐述,但该方法在国内却鲜有介绍。为此,本文在介绍B P A法理论的基础上,编制了适应性更强的求解程序,在验证算法的有效性后,将其应用于双极直流输电线路离子流场的求解与分析中,大大提高了计算效率。2 基于 BPA 的离子流场求解算法2.1 BPA 法简介B P A法是在沿通量线方向积分时,预估导线表面初始电荷密度,推导得到电压电流特性方程,引入与e相关的经验常数P,由于P的取值不同在实际计算中通常会出现不收敛的情况。为提高计算效率,本文应用导线表面电荷密度的预估初值公式来求解导线表明电荷密度,迭代求解后收敛速度更快、精度更高。2.2 BPA 法求解双极直流输电线路离子流场当输电线路电场强度超过电晕起始电场强度产生电晕现象时,周围空气分子电离所形成的电离层见图1。以负极导线电离过程为例,电离层内部空气电离生成的正、负离子在电场力的作用下,正极性离子运动至电极导线为电晕电流形成通路;负极性离子离开电离层向空间运动,或向异极性极导线大地运动,最终达到稳态。直流输电线路离子流场问题的求解关键是对空间电荷分布的求解。B P A法求解时认为电场线形成的通量管中离子流保持连续,满足电流连续性方程。基于以第4 1卷第5期邹岸新等:基于B P A法计算特高压直流输电线路合成电场图1 周围空间空气电离层F i g.1 P e r i p h e r a l s p a c e a i r i o n o s p h e r e下假设可获得计算所需的通量管(图2):不考虑电离导体周围电离层的厚度;导线为无限长直导线;离子迁移率、风速视为常数,离子流场视为二维稳态场;空间电荷仅影响电场强度的大小,不改变其方向。dldl图2 单根通量管示意图F i g.2 S c h e m a t i c d i a g r a m o f a s i n g l e f l u x t u b e图2中的单根通量管示意图,起始于导线表面,两边沿着电场线,终端截面垂直于电场强度的方向,阴影部分为通量管的上、下端面。图2中H为线路中心对地高度;R为导线半径;A2、E2分别为沿着电场线方向、距离导线表面l处的通量管界面边长和电场强度;A1、E1分别为l+dl处的通量管界面边长和电场强度;m为单位长度;V为dl范围内通量管微元的体积。针对通量管微元,结合高斯通量定理可得:l i mV01Vs?Ed?s=l i mdl0E2A2m+2E(l)A(l)dl+2E(l)dl m-E1A1mA(l)dl m(1)式中,?E为电场强度向量;d?s为通量管微元的面积分;E(l)为导线表面l处的电场强度;A(l)为导线表面l处的边长。对图1中的通量管微元V、电场强度E1、E2分别为正交于阴影部分所示的上、下端面;而电场强度E(l)平行于通量管微元的四侧面,即E(l)A(l)dl=E(l)dl m=0,由此可将式(1)化简为:l i mV01Vs?Ed?s=l i mdl0E2A2-E1A1A(l)dl(2)不考虑环境风速的离子流场控制方程:-()=/0(3)()=0(4)式中,0为空气介电常数,8.8 51 01 2F/m;为电位的梯度;为单位体积内的总电荷。当通量管微元较小时,假设在通量管内部离子流恒定不变,结合式(3)、(4)推导可得方程:E(l)A(l)dE(l)A(l)=j(0)A(0)k 0A(l)dl(5)式中,j为离子流密度;k为等效离子迁移率。将式(5)两端沿着电场线积分可得:E(l)=2j(0)k 0Ee(0)L01Ee(l)dl+E20Ee(0)21/2Ee(l)(6)式中,Ee(0)、Ee(l)分别为标称电场时沿着电场线方向,距离导线表面0、l处的电场强度;E0为导线表面的电晕起始电场强度,其值可由P e e k公式获得。引入计算因子KB=2j(0)k 0Ee(0),式(6)可化简为:E(l)=KBL01Ee(l)dl+E20Ee(0)21/2Ee(l)(7)式中,KB、E0、Ee(0)均为常 量;合成电场 强度E(l)为标称电场强度Ee(l)的积分函数,因此在标称电场强度已经求得的前提下,可通过积分法获得地面处的合成电场分布情况。在地面处合成电场强度已求得的前提下,可得电荷密度(l):(l)=j(0)A(0)k E(l)A(l)=KB0Ee(l)2E(l)(8)获得地面处的合成电场强度E(l)和电荷密度(l)后,可求解地面处的离子流密度。导线表面电荷密度0 n e w修正公式为:0 n e w=0 o l d1+fU-UiU()(9)式中,0 n e w、0 o l d分别为修正后、修正前导线表面电荷密度;f为修正函数;U、Ui分别为实际电位、计算电位。2.3 求解流程(1)计算导线表面起晕场强。应用P e e k公式计算直流输电线路导线起晕场强:Eo n=E 0m(1+k/r)(1 0)式中,Eo n为起晕场强;E 0、k均为两个经验常数;591m为导线表面的粗糙系数;为空气的相对密度;r为导线半径。(2)逐次镜像法求解输电线路标称电场强度。逐步镜像法用一个或多个模拟电荷来代替相应的输电线路分裂子导线。然后解电位系数矩阵得到每根子导线中包含的电荷量,经过反复的镜像满足精度要求后,将电荷量代入电场计算公式。计算流程见图3。在每一根导线中心设置一个模拟电荷开始设置相应的匹配点和校验点逐一计算每根子导线的电荷量逐步移动模拟电荷位置根据公式计算导线表面电场结束校验点的值是否满足精度要求是否图3 逐次镜像法计算流程F i g.3 C a l c u l a t i o n f l o w o f m i r r o r-b y-i m a g e m e t h o d(3)预估导线表面的电荷密度。导线表面电荷密度的预估初值公式7为:m=U0UdE2dU0UdE2d(1 1)式中,m为每条电场线上的平均电荷密度,对m进行适当的缩放即可作为导线表面电荷密度的预估值;U为导线电压;为无空间电荷时空间某点的电位值;为积分变量。(4)绘制电场线,形成通量管。图4描绘了双极直流输电线路附近存在电场线的分布。由图4(b)可看出,电场线从分裂子导线发出不同于从一个单根导线发出,标称电场和离子流场或多或少会因此而发生改变。B P A法全部求解流程见图5。(b)分裂子导线附近电场线(a)没有空间电荷存在时电场线分布图y/my/mx/mx/m-20-15-10-50510 15-11.5-11-10.5-10-9.512840-4-8-120.90.70.50.30.1-0.1-0.3-0.5-0.7图4 双极直流输电线路电场线F i g.4 E l e c t r i c f i e l d l i n e o f b i p o l a r D C t r a n s m i s s i o n l i n e3 BPA 法求解直流输电线路合成电场3.1 算例验证以某实际8 0 0 k V特高压直流输电线路为求解导线电场强度和电晕起始电场强度预估导线表面电荷密度初值求解合成电场强度求解地面处电荷密度求解地面处离子流密度修正导线表面节点电荷密度值求解计算因子BK三地面出发沿着电场线方向求解导线表面节点的计算电位均值Ui是否|-|U UUi图5 B P A法计算流程F i g.5 B P A m e t h o d c a l c u l a t i o n p r o c e s s例,线 路 参 数 见 表1,导 线 表 面 粗 糙 度 系 数 取0.4 5。采用B P A法求解输电 线路下距离 地面1.5 m处的合成电场强度,计算结果见图6。由图6可知,空间电荷的存在使得地面附近的合成电场强度显著增加。电场强度最大值由1 0.1 3 k V/m变为2 7.4 2 k V/m。表1 8 0 0 k V直流输电线路分裂导线参数T a b.1 8 0 0 k V D C t r a n s m i s s i o n l i n e s p l i t c o n d u c t o r p a r a m e t e r s线路参数参数值线路参数参数值极间距/m2 2线路对地高度/m2 2子导线半径/c m1.6 8避雷线半径/c m0.8 7 5分裂间距/c m4 5注:导线的分裂形式为6分裂。-60-50-40-30-20-1001020304050距离线路中心的距离/m-30-20-100102030标称电场合成电场电场强度/(kV.m)-1图6 标称电场与合成电场计算结果F i g.6 C a l c u l a t e d r e s u l t s o f n o m i n a l a n d s y n t h e t i c e l e c t r i c f i e l d s3.2 与电场线法计算结果比较采用电场线法求解合成电场强度时,与B P A法相同,需绘制相应的电场线,并提取出每步电场线节点处的电场强度,沿着电场线方向对每条电场线采用弦截迭代积分,获得导线节点电荷密度,进而求得地面处电场强度分布情况。针对8 0 0 k V特高压直流输电线路,应用电场线法和B P A法求解输电线路离子流密度与合成电场强度,结果见图7。由图7(b)可看出,B P A法和电场线法求得合成电场强度曲线整体趋势相同,数值略微存在偏差。就最大合成电场强度而言,两种方法求得的691水 电 能 源 科 学 2 0 2 3年 第4 1卷第5期邹岸新等:基于B P A法计算特高压直流输电线路合成电场-60-40-2002040距离线路中心的距离/m-50-40-30201001020304050电场线法计算结果BPA法计算结果-60-40-2002040-40-30-20-10010203040电场线法BPA法电场强度/(kV.m)-1距离线路中心的距离/m电场强度/(kV.m)-1(a)离子流密度分布图(b)合成电场强度分布图图7 B P A法与电场线法求解结果对比F i g.7 C o m p a r i s o n o f t h e r e s u l t s o f t h e B P A m e t h o d a n d t h e e l e c t r i c f i e l d l i n e m e t h o d最大出 现 位 置 相 同,且 前 者 求 得 的 最 大 值 为2 7.4 2 k V/m,后者求得最大值为2 8.4 4 k V/m。若以电场线法求得的最大值为基准,则采用B P A法所得地面合成电场强度最大值偏差为3.6%,地面离子流密度最大值偏差为3.8%。在同样的计算平台上,采用B P A法的计算时间约为电场线法计算时间的1/3,计算效率更高,再次验证了B P A法计算特高压直流输电线路合成电场的优越性。3.3 线路结构参数对合成电场强度的影响(1)导线对地高度对合成电场强度的影响。设置导线对地高度分别为1 7、1 9、2 1、2 5 m,其他线路参数保持不变,分析导线对地高度对合成电场强度的影响结果见图8。由图8可知,地面合成电场强度随导线对地高度的增加逐渐递减,这是因为提高导线对地高度,使得标称电场强度和离子流密度在同等条件下在地面附近的强度降低,从而合成电场强度大幅度减小。另外,不同高度的合成电场峰值会出现在距离线路中心1 6 m左右,在单极侧距离线路中心3 6 m左右会出现交汇的现象,在此之后各高度导线的地面合成电场减小的趋势缓慢,辐射范围也因此增大。当极间距较大,计算过程中考虑由于导线自重产生的弧垂时,同样位置下合成电场强度的峰值会增加。-60-50-40-30-20-1001020304050-50-40-30-20-1001020304050H=17 mH=19 mH=21 mH=25 m电场强度/(kV.m)-1距离线路中心的距离/m图8 高度对合成电场的影响F i g.8 E f f e c t o f h e i g h t o n t h e s y n t h e t i c e l e c t r i c f i e l d(2)极间距对合成电场强度的影响。设置导线极间距分别为2 0、2 2、2 4、2 6 m,其他线路参数保持不变,分析导线极间距对合成电场强度的影-60-50-40-30-20-1001020304050-25-20-15-10-50510152025P=20 mP=22 mP=24 mP=26 m电场强度/(kV.m)-1距离线路中心的距离/m图9 极间距对合成电场的影响F i g.9 E f f e c t o f p o l e s p a c i n g o n t h e s y n t h e t i c e l e c t r i c f i e l d响结果见图9。由图9可知,双极导线极间距为2 0、2 2、2 4、2 6 m时对应的地面合成电场峰值分别为2 0.9 7、2 2.1 9、2 3.3、2 4.2 9 k V/m,随着极间距的增加,输电线路下地面附近的合成电场强度增加,但幅度很小,峰值相差不大。极间距的变化不会影响空间电荷到达地面的距离。因此输电导线横向的移动没有纵向移动对地面合成电场的影响大,在线路架设时,导线极间距对周围电磁环境的影响可不必作为首要考虑因素,但要注意极导线之间的绝缘。(3)分裂子导线半径对合成电场强度的影响。设置分裂子导线半径分别为1 5、1 7.2、1 8.3、2 0.3 mm,其他线路参数保持不变,分析分裂子导线半径对合成电场强度的影响结果见图1 0。由图1 0可知,子导线横截面积增加,地面合成电场减小;按子导线横截面积由小到大排序,地面合成电场峰值分别减小了0.6 1、0.3 4、0.7 3 k V/m。根据规定,对于无人居住的区域,地面合成电场限制规定不超过3 0 k V/m,由此可见计算结果均满足限制要求,但由于子导线截面积对地面合成电场强度的影响不大,因此为改善电磁环境改用更粗导线的措施适用性不高。-60-50-40-30-20-1001020304050-25-20-15-10-50510152025子导线半径为15.0 mm子导线半径为.17 2 mm子导线半径为.18 3 mm子导线半径为.20 3 mm电场强度/(kV.m)-1距离线路中心的距离/m图1 0 子导线半径对合成电场的影响F i g.1 0 E f f e c t o f r a d i u s o f t h e s u b c o n d u c t o r o n t h e s y n t h e t i c e l e c t r i c f i e l d4 结论a.本文应用B P A法,结合实例求解特高压直流输电线路下地面附近合成电场强度,并对比791应用电场线法求解结果,验证了B P A法的准确性和高效性。b.提高导线的对地高度、增加分裂子导线半径和减小极间距可有效降低地面附近的电场强度,但极间距与子导线半径的影响效果很小。参考文献:1 李伟,张波,何金良.直流输电线路离子流场计算方法J.陕西电力,2 0 0 9,3 7(1 0):1-5.2 邹军,程启问,乔骥,等.高压直流输电线路离子流场计算研究综述J.南方电网技术,2 0 2 0,1 4(6):1-1 0.3 方正瑚,顾翔宇.直流离子流场的求解问题J.中国电机工程学报,1 9 9 0(增刊1):7 5-8 0.4 李永明,陈晓鑫,张淮清,等.计及地电位升高的HV D C输电线路离子流场的计算J.高压电器,2 0 1 2,4 8(1):1 7-2 4.5 乔骥,徐志威,邹军,等.一种消除D e u t s c h假设的高精度迭代特征线方法求解高压直流输电线路离子流场 J.电 工 技 术 学 报,2 0 1 8,3 3(1 9):4 4 1 9-4 4 2 5.6 汪沨,范竞敏,李敏,等.高精度上流有限元法在特高压直流输电线路离子流场计算中的应用J.高电压技术,2 0 1 6,4 2(4):1 0 6 1-1 0 6 7.7 S A RMA M P,J AN I S CHEWS KY J W.A n a l y s i s o f c o r o n a l o s s e s o n d c t r a n s m i s s i o n l i n e s:i-u n i p o l a r l i n e sJ.I E E E t r a n s a c t i o n s o n p o w e r a p p a r a t u s a n d s y s t e m s,1 9 6 9,P A S-8 8(5):7 1 8-7 3 1.C a l c u l a t i o n o f S y n t h e t i c E l e c t r i c F i e l d o f U H V D C T r a n s m i s s i o n L i n e B a s e d o n B P A M e t h o dZ OU A n-x i n1,WANG S h o u-g a n g1,YANG T a o1,L I Y o n g-m i n g2,L I U Y a n2(1.S t a t e G r i d C h o n g q i n g E l e c t r i c P o w e r C o m p a n y U l t r a H i g h V o l t a g e B r a n c h,C h o n g q i n g 4 0 0 0 3 9,C h i n a;2.S t a t e K e y L a b o r a t o r y o f P o w e r T r a n s m i s s i o n E q u i p m e n t&S y s t e m S e c u r i t y a n d N e w T e c h n o l o g y,C h o n g q i n g U n i v e r s i t y,C h o n g q i n g 4 0 0 0 4 4,C h i n a)A b s t r a c t:I n o r d e r t o s t u d y t h e i o n f l o w f i e l d p r o b l e m n e a r t h e g r o u n d o f UHV D C t r a n s m i s s i o n l i n e s,t h e B P A m e t h o d i s a p p l i e d t o s o l v e t h e i o n f l o w f i e l d o f b i p o l a r D C t r a n s m i s s i o n l i n e s.T h e s y n t h e t i c e l e c t r i c f i e l d i s s o l v e d c o m-b i n e d w i t h t h e n o m i n a l e l e c t r i c f i e l d c a l c u l a t e d b y t h e s u c c e s s i v e i m a g e m e t h o d.T h e i o n f l o w f i e l d a n d s y n t h e t i c e l e c t r i c f i e l d p r o b l e m s o f 8 0 0 k V D C t r a n s m i s s i o n l i n e s a r e c a l c u l a t e d u s i n g t h e B P A m e t h o d,a n d t h e e f f e c t s o f e r e c t i o n h e i g h t,i n t e r e l e c t r o d e d i s t a n c e,a n d s u b-c o n d u c t o r r a d i u s o n t h e s y n t h e t i c e l e c t r i c f i e l d o f UHV D C t r a n s m i s s i o n l i n e s a r e a n a-l y z e d.I n a d d i t i o n,t h e p o w e r-l i n e m e t h o d f o r s o l v i n g t h e s y n t h e t i c e l e c t r i c f i e l d s t r e n g t h i s c o m p a r e d.T h e r e s u l t s s h o w t h a t t h e B P A m e t h o d i s a c c u r a t e a n d e f f e c t i v e i n s o l v i n g t h e i o n f l o w f i e l d,a n d t h e c a l c u l a t i o n e f f i c i e n c y i s g r e a t l y i m-p r o v e d.I n c r e a s i n g t h e h e i g h t o f t h e w i r e e r e c t i o n a n d i n c r e a s i n g t h e r a d i u s o f t h e s p l i t s u b-c o n d u c t o r c a n i m p r o v e t h e e-l e c t r o m a g n e t i c e n v i r o n m e n t n e a r t h e g r o u n d,w h i l e r e d u c i n g t h e i n t e r e l e c t r o d e d i s t a n c e c a n r e d u c e t h e s y n t h e t i c e l e c t r i c f i e l d i n t e n s i t y n e a r t h e g r o u n d,b u t t h e e f f e c t i s n o t o b v i o u s.K e y w o r d s:B P A(B o n n e v i l l e P o w e r A d m i n i s t r a t i o n)m e t h o d;UHV D C t r a n s m i s s i o n;s y n t h e t i c e l e c t r i c f i e l d;h e i g h t t o g r o u n d;p o l e s p a c i n g;s p l i t s u b c o n d u c t o r r a d i u s(上接第1 2 0页)3 王冰清.基于回溯式搜索算法的随机神经网络优化及应用D.济南:山东大学,2 0 1 6.4 田文凯.回溯搜索优化算法的研究和应用D.西安:西安电子科技大学,2 0 1 4.5 蒋中明,徐卫亚.三维初始地应力场反分析的径向基函数法J.岩土力学,2 0 0 2,2 3(6):7 3 7-7 4 1.6 王红彬,石焱炯,沈德虎,等.白鹤滩水电站巨型地下厂房高边墙开挖支护变形研究J.水电能源科学,2 0 2 0,3 8(5):1 1 8-1 2 1.I n i t i a l G r o u n d S t r e s s I n v e r s i o n a n d S u r r o u n d i n g R o c k S t a b i l i t y A n a l y s i s o f U n d e r g r o u n d P o w e r h o u s e o f W u y u e P u m p e d S t o r a g e P o w e r S t a t i o nR E N B i n1,Z HANG F a n1,L I Y o n g-s h e n g1,Z HOU Y o n g-j i a n2,Z HANG Z h o u-r u i2,CHE N K a i1,L I U Y u n1(1.H e n a n X i n h u a W u y u e P u m p e d S t o r a g e P o w e r G e n e r a t i o n C o.,L t d.,X i n y a n g 4 6 5 4 5 0,C h i n a;2.S t a t e K e y L a b o r a t o r y o f G e o h a z a r d P r e v e n t i o n a n d G e o e n v i r o n m e n t P r o t e c t i o n,C h e n g d u U n i v e r s i t y o f T e c h n o l o g y,C h e n g d u 6 1 0 0 5 9,C h i n a)A b s t r a c t:T h e u n d e r g r o u n d p o w e r h o u s e o f H e n a n Wu y u e P u m p e d S t o r a g e P o w e r S t a t i o n h a s a c o m p l e x s t r u c t u r e o f c a v e r n s.T h e m a i n p o w e r h o u s e h a s a m a x i m u m s p a n o f 2 6.0 m a n d a h e i g h t o f 5 5.1 m,w h i c h i s a l a r g e u n d e r g r o u n d c h a m b e r.T h e r o c k m a s s s t r u c t u r a l p l a n e s a r e d e v e l o p e d a n d t h e s t r e n g t h i s l o w.T h e f a u l t F 1 r u n s t h r o u g h t h e t h r e e p o w e r h o u s e,a n d t h e s a f e t y p r o b l e m s a r e p r o m i n e n t.T h e s t a b i l i t y o f s u r r o u n d i n g r o c k a f t e r e x c a v a t i o n s h o u l d b e a n a-l y z e d.T h e r e f o r e,F L A C3 D p r o g r a m w a s u s e d t o e s t a b l i s h t h e g e o l o g i c a l m o d e l.T h e B S A-B P n e u r a l n e t w o r k w a s u s e d t o c a l c u l a t e t h e i n i t i a l i n-s i t u s t r e s s f i e l d o f t h e a r e a.A n d t h e n t h e e x c a v a t i o n o f u n d e r g r o u n d p o w e r h o u s e w a s s i m u l a t e d.T h e s t r e s s f i e l d,d i s p l a c e m e n t f i e l d a n d p l a s t i c z o n e d e v e l o p m e n t a n d d i s t r i b u t i o n l a w o f s u r r o u n d i n g r o c k w e r e a n a l y z e d.T h e r e s u l t s s h o w t h a t t h e s u r r o u n d i n g r o c k o f u n d e r g r o u n d p o w e r h o u s e c a v e r n s g r o u p i s b a s i c a l l y s t a b l e.H o w e v e r,l o c a l d e f o r m a t i o n a n d f a i l u r e m a y o c c u r i n t h e f a u l t d e v e l o p m e n t a r e a a n d a t b o t h e n d s o f t h e b u s b a r h o l e.T h e r e s e a r c h r e s u l t s c a n p r o v i d e g u i d a n c e f o r t h e d e s i g n a n d c o n s t r u c t i o n o f t h e p r o j e c t,a n d p r o v i d e s o m e r e f e r e n c e f o r s i m i l a r p r o j e c t s.K e y w o r d s:p u m p e d s t o r a g e p o w e r s t a t i o n;u n d e r g r o u n d p o w e r h o u s e c a v e r n s;F L A C3 D;g e o s t r e s s i n v e r s i o n;B S A-B P n e u r a l n e t w o r k;s u r r o u n d i n g r o c k s t a b i l i t y891水 电 能 源 科 学 2 0 2 3年