矩形复合接地网基础接地电阻的分析计算

1、矩形复合接地网基础接地电阻的分析计算 平帅1引言接地涉及防雷系统的安全，是防雷系统极为关注的问题之一，GB50057建筑物防雷设计规范要求尽量使用建筑物基础接地作为建筑物防雷的接地，但是建筑物基础接地的工频接地电阻的计算一直存在疑意，为有关部门所关注。由于土壤的电阻率并非常数，对竖直方向不均匀的土壤，常用二层土壤模型代替，再进行分析5。用建筑物基础中的钢筋作为接地体，具有耐腐蚀、可节约钢材和节省占地面积等优点，均匀土壤中水平或竖直对称地网基础接地电阻的计算可见文献14。建筑物的地基常是多层水平埋设与竖直埋设的钢筋复合体，其接地电阻无法用解析法计算。本文对二层土壤中这种复合地网的基础接地电阻，用

2、矩量法与静电比拟法推导出计算式。外部防雷的环形接地体尚宜按以下方法敷设：1)当土壤电阻率小于或等于500m时，对环形接地体所包围的面积的等效圆半径大于或等于5 m的情况，环形接地体不需补加接地体；对等效圆半径小于5 m的情况，每一引下线处应补加水平接地体或垂直接地体。当补加水平接地体时，其长度应按下式确定： (4.2.4-1)式中 补加水平接地体的长度(m)； A环形接地体所包围的面积(m2)。当补加垂直接地体时，其长度应按下式确定： (4.2.4-2)式中 补加垂直接地体的长度(m)。2)当土壤电阻率为500m至3000m时，对环形接地体所包围的面积的等效圆半径大于或等于 m的情况，环形接地

3、体不需补加接地体；对等效圆半径小于 m的情况，每一引下线处应补加水平接地体或垂直接地体。当补加水平接地体时，其总长度应按下式确定： (4.2.4-3)当补加垂直接地体时，其总长度应按下式确定： (4.2.4-4)注：按本款方法敷设接地体时，每根引下线的冲击接地电阻可不作规定。共用接地装置的接地电阻按50 Hz电气装置以人身安全所要求的阻值确定。共用接地装置的接地电阻应按50 Hz电气装置以人身安全所要求的阻值确定。在土壤电阻率小于或等于3000m的条件下，外部防雷装置的接地体当符合下列规定之一时可不计及冲击接地电阻；但当每根专设引下线的冲击接地电阻不大于10时，可不按本条1、2款敷设接地体。1

4、 当土壤电阻率小于或等于800m时，对环形接地体所包围的面积的等效圆半径大于或等于5 m的情况，环形接地体不需补加接地体；对等效圆半径小于5 m的情况，每一引下线处应补加水平接地体或垂直接地体。当补加水平接地体时，其长度应按本章(4.2.4-1)式确定；当补加垂直接地体时，其长度应按本章(4.2.4-2)式确定。 当土壤电阻率为800m至3000m时，对环形接地体所包围的面积的等效圆半径大于或等于 m的情况，环形接地体不需补加接地体；对等效圆半径小于 m的情况，每一引下线处应补加水平接地体或垂直接地体。当补加水平接地体时，其总长度应按下式确定： (4.3.6-1)当补加垂直接地体时，其总长度应

5、按下式确定： (4.3.6-2)第五款，关于共同接地：由于防雷装置直接装在建、构筑物上，要保持防雷装置与各种金属物体之间的安全距离已成为不可能。此时，只能将屋内各种金属物体及进出建筑物的各种金属管线，进行严格的接地，而且所有接地装置都必须共用，并进行多处连接，使防雷装置和邻近的金属物体电位相等或降低其间的电位差，以防反击危险。一般说来，接地电阻越低，防雷得到的改善越多。但是，不能由于要达到某一很低的接地电阻而花费过大。出现反击危险可以从基本计算公式来评价，IR项对于建筑物内某一小范围中互相连接在一起的金属物（包括防雷装置）说来都是一样的，它们之间的电位差与防雷装置的接地电阻无关。此外，考虑到已

6、采取严格的各种金属物与防雷装置之间的连接和均压措施，故不必要求很低的接地电阻。从防雷观点出发，较好是设共用接地装置，它适合供所有接地之用（例如防雷、低压电力系统、电讯系统）。第六款，为了将雷电流流散入大地而不会产生危险的过电压，接地装置的布置和尺寸比接地电阻的特定值更重要。然而，通常建议有低的接地电阻。本款的规定完全采用IEC1024l防雷标准2.3.3.2的规定（接地体的B型布置）。图32系根据该标准的图2换成本规范的防雷建筑物类别的图。该标准对接地体B型布置的规定是：对于环形接地体（或基础接地体），其所包围的面积的平均几何半径r应不小于l1，即rl1，l1示于图 2（相当于本说明的图3.2

7、）；当l1大于r时，则必须增加附加的水平放射形或垂直（或斜形）导体，其长度lr（水平）为lr l1r或其长度lv（垂直）为。环形接地体（或基础接地体），其所包围的面积A的平均几何半径r为：r2 A，所以。根据图32，对于第一类防雷建筑物，当500m时l1为m，因此，导出本款第1项的规定；当5003000m时，l1与的关系是一根斜线，从该斜线上找出方便的任两点的坐标，则可求出l1与的关系式为，所以，导出本款第2项的规定。 由于接地体通常靠近墙、基础敷设，所以补加的水平接地体一般都是从引下线与环形接地体的连接点向外延伸，可为一根，也可为多根。 由于本条采用了若干等电位措施，本款的接地电阻值不是起主

8、要作用，因此，没有提出接地电阻值的具体要求。本款所要求的环形接地体的工频接地电阻R，在其半径r等于l1 的场合下，当5003000m时，大约处于3313；500m时，R0.067（）。 环形接地体的工频接地电阻的计算式为（），（m）。式中为土壤电阻率（mf），A为环形接地体所包围的面积（m2）。当500m、d10m时，（）。当500m、时，（）。2 计算原理设地中有一尺寸为abH的矩形混凝土立方体(图1)，内有一复合钢筋网，其水平网有K层，形状相同，外尺寸为L1L2。第一层距地面为Z0，相邻两层距离为OZ，钢筋半径为r，竖直钢筋长为l，水平与竖直钢筋的半径可相等或不等，混凝土及二层土壤的电阻率

9、分别为c、1和2，上层土壤厚度为HUP，设HUPH。由于是计算工频接地电阻，电流频率低，可按恒定电流场考虑。为了消去空气的影响 ，利用恒定电流场的镜象法将地面以下的部分对称地补充到地面以上来(如图2)，它以地表面为对称面上下对称。用图2计算时，原钢筋和镜象钢筋必须流有相同的电流，以保证边界条件与图1相同，即地表面与电流线重合。根据唯一性定理，图1、图2下半空间场的解答相同 ，可用图2计算图1的电流场。根据静电比拟原理，将图2的恒定电流场转化为图2的静电场处理，为了保证c/c1/12/2(其中c、1和2为混凝土及土壤的导电系数)，c1/c，11/1，21/2，故图2中的c、1和2应分别取c1/c

10、，11/1，21/2，这样，当给定电极电位0时，若能算出导体上的自由电荷Q，则可得出系统的对地电容C，由静电比拟法，可得出图2系统以及图1系统的接地电阻。图2中，存在3种电介质c、1和2，c与1的分界面为六面立方体，1与2的分界面为无穷大平面，c电介质内有导体，难以用解析法计算，宜用数值计算法求解。电介质在电场作用下被极化，在介质分界面上出现束缚电荷。为了便于后面用矩量法计算,设想将介质去掉，在介质c、1的分界面上和介质1、2的分界面上放置束缚电荷，在介质c内的接地电极上放置等效电荷，此等效电荷由束缚电荷及电极上的自由电荷所组成。当以束缚电荷代替电介质作用后，便可将图2的静电场看成置于真空0中

11、，场域内任一点M(x，y，z)的电位可表示为 (1) 式中ss1s2，s1为导体表面积，对应的为导体表面的等效面电荷密度，s2为介质分界面，对应的为介质分界面上的束缚电荷面密度；ds为面s上的面积元；rMP为源点P至场点M的距离。设导体电位为0，当将场点M取在导体上时，导体钢筋上的电位为 (2) 在介质c与1的分界面上，应满足边界条件,式中n为法线。用电荷密度来表示此式，利用式(1)分别算出此式中左右二项，在分界面P0点上，经推导可得1，2 (3)式中s0为边界面上以点P0为球心、半径为无穷小的半球面；nP0为边界面上点P0处的法线矢量；P0为P0处的束缚电荷面密度。 同理，对1与2分界面上的

12、边界条件，经类似计算后可得 (4)式(2)(3)(4)为积分方程，积分号内电荷密度(x，y，z)为待求量，用矩量法对这3式进行离散化，并求其数值解。以式(3)为例，水平网有K层，将每层剖分为m0个单元，水平钢筋上剖分数为mHKm0，将竖直钢筋剖分为mV个单元。钢筋上总的剖分数为mmHmV。将介质c与1的分界面剖分为m1个单元，介质1与2的分界面为无穷大平面，取有限大面积AB近似代替(图2)，将其剖分为Nmm1个单元，N是系统总剖分数。选脉冲函数为基函数，意味着每个剖分单元面积sj上的面电荷密度j为常数，则电荷密度可表示为其中fj取1或0，当点(x,y,z)在sj内时，则取1，否则为0。再选狄拉

13、克函数为权函数Wi，第i个剖分单元坐标定义为其面积中心，设点M的坐标为(x，y，z)，点i的坐标为(xi，yi，zi)，则有（5）对式(3)应用矩量法，注意到P0是在分界面上的任意一点，且设其在分界面上的第i个单元上，则式(3)左边为 （6）由于P0在第i个单元上，ds上的电荷面密度为，点P在ds上，则式(3)右边为（7）令，ni为边界面上第i个单元的面法线矢量。式（3）的离散结果为 （8）因系统上下对称，对称点上电荷密度相等，未知量仅为上半空间的电荷密度，只需对上半空间进行剖分，这样便可使方程数减少一半，节省了的内存，加快了运算速度。用0分别表示源点j以及j在下半空间的对称点至场点i的距离。

14、考虑到复合地网，则式(8)变为 （9）对式（2）（4）进行类似运算后，在钢筋上有 式中第一项对应水平钢筋，第二项对应竖直钢筋，它表示当ij时，水平钢筋网上第i个剖分单元对其本身的电位系数，当ij时，Laii0；Lbii的形式与Laii相同，它表示当ij时，竖直钢筋上第i个剖分单元对其本身的电位系数 ，此时r表示竖直钢筋半径，当ij时，Lbii0；j为第j个剖分单元上的面电荷密度；sj为第j个剖分单元面积；L为钢筋剖分单元长度；0为钢筋电位；N为系统剖分单元总数。在介质1与2之分界面上有 （11）式(9)(10)和(11)分别是式(3)(2)和(4)用矩量法离散后得出的线性代数方程，共N个。据此

15、可进行编程，给定0后，用LU分解法解方程组，求出各单元上的电荷，算出钢筋上的总等效电荷，钢筋上的总自由电荷Q=(c/0)Qeq,钢筋对地电容CQ/0，前面已取c1/c，11/1，21/2，如计算时再设置01，01，则接地电阻R1C。3 程序验证均匀土壤和二层土壤接地电阻的计算以及均匀土壤基础接地电阻的计算均可看为本文方法的特例。以下用文献78943的计算结果与本文程序计算结果进行比较。(1)复合地网验证文献7对均匀土壤中一复合地网进行了计算，水平钢筋半径r0007 m，L1L24572 m4572 m，N1N26，Z00457 m，在每个节点上放置长为1524 m、半径为0007 94 m的竖

16、直钢筋，共36根(图3)，土壤电阻率为100m，得接地电阻为0.774 。文献8算得接地电阻为0746。本文程序中取12c100m，算得接地电阻为0755 3。若将地网放入尺寸为abH50 m50 m18 m的混凝土中央，设土壤为二层，HUP20 m，本文程序可算得：当c99m，1100m，2101m时，接地电阻为0756 7；当c101m，1100m，299m时，接地电阻为0.754 0。(2)双层土壤验证文献9给出双层土壤中方形水平地网接地电阻计算公式为 式中 HUP!endif 0；hHUP；Cf为网格形状系数(正方形网格为09）；Lx和Ly分别为网格在a与b方向的长度，m；L为网格导体

17、总长，m；d0为导体直径，m；h为埋深，m；A0为地网总面积L1L2，m2；HUP为上层土壤厚度，m；1和2分别为上、下层土壤的电阻率，m。取hZ02 m，1100m，2300m，HUP5 m，r0008 m，L1L219 m39 m，N1N220，由式(12)得接地电阻R2692 7。本文程序取c1100m ，得接地电阻为2710 7。若地网放入尺寸为abH20 m40 m3 m的混凝土中，取c099199m时，本文程序算得接地电阻为2709 5；取c101m时，本文程序算得接地电阻为2712 0。(3)有混凝土时的验证均匀土壤中单层水平地网的基础接地电阻的计算见文献4，如：abH40 m2

18、0 m2 m，L1L239 m19 m，N120，N220，Z01 m，r0008 m，c500m，土壤电阻率为100m，则文献4算得接地电阻为1925。本文程序取12100m，则算得接地电阻为1917。而均匀土壤中多层水平地网基础接地电阻的验证可见文献3，本文程序计算结果与之相同，不再重复。以上数据表明，本文程序计算结果与文献78943的计算结果很接近，表明了本文计算方法和所编程序的正确性。4 算例取图1的地网参数如表1，竖直钢筋位于水平网的各个联接点上(图4)，顶端与地面相齐 ，将钢筋剖分为1 m、混凝土表面剖分为1 m1 m，AB面积剖分为2 m2 m的单元，取不同的A/a和B/b，可得

19、表2。随A/a和B/b的增大，电阻值越来越接近，A/a、B/b取10与取15的结果，其误差已在15%以内。计算时取A/a、B/b为10便可满足工程计算的精度要求。应用本文程序时，输入如表1的各项数据和c、1及2，选取适当的剖分单元尺寸及A/a、B/b的大小(一般为10)，便可算出接地电阻。限于篇幅，其它计算结果未列出。5 结论经文献78943的计算结果比较，表明本文的分析方法、计算式及所编程序是正确的。本文程序可方便地计算均匀土壤、二层土壤中常规接地和基础接地的接地电阻，钢筋可为水平网、竖直网以及由单层或多层水平网与竖直网所构成的复合体。参考文献1胡登宇应用矩量法计算混凝土中钢筋接地电阻J广西电力技术，1998，(1)：20222胡登宇矩形混凝土中水平地网接地电阻的计算J长沙电力学院学报，1999，(4)：3343363胡登宇，陈彩屏多层水平地网基础接地电阻的计算C99全国电磁兼容学术会议(青岛)论文集，1999：3543574Kostic M B，Popovic B D，Jov