1.一种城市电网变压器中性点直流计算方法，其特征在于，包括以下步骤：步骤1：获取城市地铁车站地理信息、牵引系统及回流系统相关电气参数、列车位置及牵引电流，获取城市电网中变电站地理信息、变电站主接线图、变压器及线路相关电气参数；

步骤2：建立地铁与城市电网直流电阻网络模型；

步骤3：模型参数赋值，计算城市电网直流电阻网络模型各节点电压；

步骤4：建立城市电网变压器母线及中性点节点阻抗矩阵，计算主变中性点直流。

2.根据权利要求1所述的一种城市电网变压器中性点直流计算方法，其特征在于，所述步骤1具体为：

S11：记录地铁线路车站总量N，包含牵引变电站的地铁车站数量为n1，地铁车辆段及停车场总数量为n2，地铁主所数量为A，选取地铁最长线路末端车站为原点，建立坐标系，记录各地铁车站坐标W＝[w1(x1,y1)，w2(x2,y2)，…，wN(xN,yN)]，记录地铁主所坐标B＝[b1(x1,y1)，b2(x2,y2)，…，bA(xA,yA)]，记录地铁相邻车站间距L＝[L1，L2，…，LN‑1]，牵引变电站内阻r，地铁直流供电电压等级为Vr，单根钢轨每千米直流电阻r1，每千米钢轨对地电阻r2，接触网每千米直流电阻，地铁贯通地线每千米直流电阻，35kV电缆铠装每千米直流电阻，35kV电电缆长度，车站及场段接地电阻rg1，场段与进入场段的车站之间线路距离为La＝[l1，l2，…，ln2]；车站钢轨限位装置等效电阻rov，单向导通装置等效电阻rd；地铁线路中行驶列车数量为n，记录列车在线路中的位置，列车牵引电流It＝[I1，I2，…，In]；

S12：记录城市电网变电站总数M，中性点接地自耦变压器数量为m1，中性点接地三绕组变压器数量为m2，各变电站坐标D＝[d1(x1,y1)，d2(x2,y2)，…，dM(xM,yM)]，变压器高压绕组直流电阻，变压器中压绕组电阻，三相输电线路每千米直流电阻r6，线路回数，架空地线每千米直流电阻，三相输电电缆缆芯每千米直流电阻，输电电缆回数、电缆铠装每千米直流电阻，变电站接地电阻，各变电站之间输电线路H＝[h1，h2，…，hp]，各城市电网变电站之间输电电缆距离为[g1，g2，…，gq]；

S13：记录为地铁牵引主所110kV电缆长度，三相电缆铠装层每千米纵向直流电阻r10。

3.根据权利要求2所述的一种城市电网变压器中性点直流计算方法，其特征在于，所述步骤2具体为：

S21：根据地铁线路拓扑，将地铁线路等效为集中参数的直流电阻网络模型，模型包含上/下行接触网、钢轨、贯通地线、参考地结构，模型节点总数Nm＝2×n1+A+2×N+2×n+2×n2，模型中，上/下行接触网结构模型包含牵引变电站节点及列车节点，钢轨结构模型包含车站节点、列车节点、车辆段及停车场节点；贯通地线结构模型包含车站节点、车辆段及停车场节点，接触网、钢轨及贯通地线结构模型中纵向相邻节点之间直流电阻＝接触网、钢轨及贯通地线每千米直流电阻×相邻节点距离；

S22：地铁牵引变电站等效为直流电阻并联直流电流源结构，连接在接触网牵引变电站节点与钢轨牵引变电站节点之间，直流电阻值为r，并联电流源为Vr/r；地铁线路中钢轨车站节点与贯通地线节点之间的直流过渡电阻R2＝r2/([L1，L1，L2，…，LN‑1]/2)+r2/([L1，L2，…，LN‑1，LN‑1]/2)，车辆段及停车场钢轨节点与贯通地线节点之间过渡电阻为0.35，车站及场段钢轨限位装置等效为直流电阻rov，连接在钢轨车站节点与贯通地线节点之间，与R2并联；场段钢轨节点与车站钢轨节点之间电阻为rd+r1×La/2，地铁主所贯通地线节点与牵引变电站节点之间通关电缆铠装层连接，节点之间电阻为，车站、场段及主所接地系统等效为电压源串连直流电阻，直流电阻为rg1，电压源为[U1，U2，…，UN+n2+A]，连接车站及场段贯通地线节点与参考地节点之间，列车等效为直流电流源，电流源电流为列车牵引电流It，连接牵引网列车节点与钢轨列车节点；

S23：根据城市电网拓扑结构，将城市电网等效为集中参数的直流电阻网络模型，模型包含变压器母线节点、变电站接地节点，节点总数为Nn＝M+2×m1+m2，各变电站母线节点间电阻＝三相输电线路或电缆每千米纵向直流电阻×输电线路或电缆长度/输电线路或电缆回数，各变电站接地节点之间电阻＝架空地线或电缆铠装每千米纵向直流电阻×输电线路或电缆长度/输电线路或电缆回数；自耦变压器母线节点间电阻为自耦变压器高压绕组直流电阻，自耦变压器母线节点与变电站接地节点之间的电阻为变压器中压绕组直流电阻，三绕组变压器母线节点与接地节点之间电阻为变压器高压绕组与中压绕组直流电阻之和；

变电站接地节点与大地节点通过直流电阻与电压源连接，直流电阻为rg2，电压源为[U1，U2，…，UM]；

S24：基于地铁主所及城市电网变电站接地网之间电缆铠装连接关系，耦合地铁直流电阻网络模型及城市电网直流电阻网络模型，模型节点总数为NA＝Nm+Nn，地铁直流电阻网络模型中的主所接地节点与城市电网直流电阻网络模型变电站接地节点间直流电阻＝地铁牵引主所110kV电缆长度×r10。

4.根据权利要求3所述的一种城市电网变压器中性点直流计算方法，其特征在于，所述步骤3具体为：

S31：根据模型节点之间直流电阻，确定模型节点导纳矩阵Y，Y矩阵大小为NA×NA；确定模型电流激励I，I矩阵大小为NA×NA，矩阵中非列车节点及牵引变电站节点电流激励均为0；

确定接地节点与模型节点的关联矩阵Hg，Hg矩阵大小为NA×(N+n2+A+M)，模型节点与接地节点为同一节点时数值为1，其余为0；确定模型接地节点接地电阻矩阵Gg，Gg矩阵大小为N+n2+A+M，数值为节点节点与参考地节点之间的直流电阻；接地节点之间互阻抗矩阵为Mg，Mg矩阵大小为(N+n2+A+M)×(N+n2+A+M)；接地节点串联电压源矩阵U，U矩阵大小为N+n2+A+M，数值为各接地节点串联的电压U＝[U1，U2，…，UN+n2+A+M]；变电站接地节点电压VM，VM矩阵大小为M；

S32：基于上述矩阵，计算模型接地节点等效电流矩阵Jg＝GgU，模型节点电流矩阵J＝I+T

HgJg；节点点电压V与节点电流J满足YV＝J；接地节点电压为Vg＝(Hg) V，Vg与接地节点入地T

电流Ig满足关系Ig＝Gg(Vg‑U)；节点电压V与模型电流激励I满足关系：V＝((E+HgGgMgHg)Y‑T T

HgGgMgGgHg)(E+HgGgMgHg)I；

5 5

S33：默认地铁钢轨限位装置及单向导通装置不导通，rov＝10 Ω，rd＝10Ω，根据节点电压与模型电流激励之间的关系，计算节点电压；判断钢轨车站节点电压是否超过120V，如‑5

果超过则该节点rov＝10 Ω；计算场段钢轨节点电压‑车站钢轨节点电压，若该值大于0.8V‑5

或者小于‑8V，则rd＝10 Ω；重新计算模型节点电压。

5.根据权利要求4所述的一种城市电网变压器中性点直流计算方法，其特征在于，所述步骤4具体为：

S41：变压器母线节点与变电站接地节点电阻矩阵为Φ，Φ矩阵大小为m1+m2；变压器母线节点电压为VΦ，变压器母线节点电压VΦ及变电站接地节点电压VM与主变中性点直流In之间满足：In＝Φ(VΦ‑Vg)；

S42：基于节点电压计算结果，基于In＝Φ(VΦ‑Vg)计算主变中性点直流电流。