CN112765897B - 一种基于架空输电线路的二维磁感应强度控制方法 - Google Patents
一种基于架空输电线路的二维磁感应强度控制方法 Download PDFInfo
- Publication number
- CN112765897B CN112765897B CN202110131901.4A CN202110131901A CN112765897B CN 112765897 B CN112765897 B CN 112765897B CN 202110131901 A CN202110131901 A CN 202110131901A CN 112765897 B CN112765897 B CN 112765897B
- Authority
- CN
- China
- Prior art keywords
- magnetic induction
- current
- phase
- magnetic
- transmission line
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 230000006698 induction Effects 0.000 title claims abstract description 86
- 230000005540 biological transmission Effects 0.000 title claims abstract description 35
- 238000000034 method Methods 0.000 title claims abstract description 34
- 239000004020 conductor Substances 0.000 claims abstract description 47
- 238000004364 calculation method Methods 0.000 claims abstract description 42
- 239000002245 particle Substances 0.000 claims abstract description 9
- 238000004088 simulation Methods 0.000 claims description 19
- 238000005457 optimization Methods 0.000 claims description 4
- 230000035699 permeability Effects 0.000 claims description 3
- 238000004422 calculation algorithm Methods 0.000 abstract description 5
- 238000010276 construction Methods 0.000 description 3
- 230000007547 defect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000007792 addition Methods 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/27—Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/004—Artificial life, i.e. computing arrangements simulating life
- G06N3/006—Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/04—Power grid distribution networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/16—Cables, cable trees or wire harnesses
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- Software Systems (AREA)
- Artificial Intelligence (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- General Health & Medical Sciences (AREA)
- Mathematical Physics (AREA)
- Computational Linguistics (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- Biophysics (AREA)
- Biomedical Technology (AREA)
- Data Mining & Analysis (AREA)
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Medical Informatics (AREA)
- Computer Hardware Design (AREA)
- Geometry (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及一种基于架空输电线路的二维磁感应强度控制方法,通过粒子群算法,在等效导线横截面上,对构造出的适应度函数寻找模拟电流的个数和位置的最优值,然后再用得到的模拟电流最优个数和位置来对架空输电线路下方的磁感应强度进行求解。与现有技术相比,本发明具有提高磁感应强度分布计算结果的准确性、降低磁感应强度分布的计算量和计算成本、减少时间成本等优点。
Description
技术领域
本发明涉及电网建设技术领域,尤其是涉及一种基于架空输电线路的二维磁感应强度控制方法。
背景技术
随着超高压和特高压电网持续建设,现有的电网规模不断扩大,电网电能输送能力不断提高,输电线路的运行电压等级也随之提高。在满足用电负荷增加需求的同时,也使得架空输电导线下方的工频磁场环境有所恶化,导致居住在输电线路周围的居民对居住区域磁场干扰水平的担忧越来越严重。
目前常用的磁感应强度分布计算分析方法多为软件的仿真或理论的计算,仿真的软件有COMSOL,理论计算方法包括有限元法、模拟电流法等。虽然通过仿真能够比较准确的得出架空输电线路下方的磁感应强度分布情况,但仿真软件不仅对计算机的配置要求很高,操作起来还很复杂,尤其是建模十分繁琐,最后还存在计算时间过长的缺点。现有的有限元法也是随着计算机性能的提升而被使用在磁感应强度的计算上,是通过MAXWELL软件内置计算方法进行计算,为保证计算精度,通常会划分很密的网格,计算量十分庞大,耗时耗力。模拟电流法对磁感应强度的计算就较为简单,但是计算精度不够,误差较大。
发明内容
本发明的目的就是为了克服上述现有技术存在的磁感应强度分布情况的计算耗时较大、精度不足、误差较大的缺陷而提供一种基于架空输电线路的二维磁感应强度控制方法。
本发明的目的可以通过以下技术方案来实现:
一种基于架空输电线路的二维磁感应强度控制方法,具体包括以下步骤:
S1、获取架空输电导线的三相导线中分裂导线的参数,计算分裂导线的等效半径;
S2、获取三相导线中初始模拟电流的个数和位置,并在三相导线的表面设置匹配点和校验点;
S3、获取预设的零磁位参考点,通过构建方程组计算得到初始模拟电流的电流值和匹配点的磁位,根据所述初始模拟电流的电流值计算出校验点的磁位,根据匹配点的磁位和校验点的磁位计算得到平均磁位误差;
S4、根据初始模拟电流的电流值计算出架空输电导线下方的第一磁感应强度分布,并获取有限元法计算得到的第二磁感应强度分布,第一磁感应强度分布与第二磁感应强度分布对比得到平均磁感应强度分布误差;
S5、根据平均磁位误差和平均磁感应强度分布误差计算适应度函数的函数值,通过粒子群算法进行优化得到函数值最小时对应的初始模拟电流的个数和位置,作为最终模拟电流的个数和位置,并计算出每个最终模拟电流的数值;
S6、根据每个最终模拟电流的数值,计算出各个最终模拟电流产生的磁感应强度,并进行叠加,得到架空输电线路下方的最终磁感应强度分布。
所述零磁位参考点具体为三相导线中A相导线表面的匹配点。
所述三相导线中每相导线的初始模拟电流的位置分别位于相应的圆周上。
进一步地,所述圆周的半径小于三相导线的半径,所述圆周的圆心和三相导线的圆心重合。
所述等效半径的计算公式如下所示:
其中,req为等效半径,R为分裂导线的半径,n为分裂导线的根数,r为次导线的半径。
所述三相导线中每相导线上匹配点、校验点和初始模拟电流的数量相同,三相导线之间的匹配点、校验点和初始模拟电流的数量相同。
所述匹配点的磁位的计算公式如下所示:
进一步地,所述步骤S3中方程组根据磁位系数、磁感应强度和磁感应强度系数进行构建,磁位系数的计算公式如下所示:
磁感应强度的计算公式如下所示:
磁感应强度系数的计算公式如下所示:
所述步骤S3中构建的方程组具体如下所示:
其中,Ij为初始模拟电流的电流值,rij为第j个初始模拟电流到第i个匹配点的距离,r1j为第j个模拟电流到磁位参考点的距离,所述磁位参考点即为第一个匹配点,II、Iп和Iш分别为三相导线A相电流、B相电流和C相电流的大小,N为三相导线中单相导线上匹配点的数量,为B相导线表面匹配点的磁位,为C相导线表面匹配点的磁位。
所述适应度函数为粒子群算法的目标函数,适应度函数的计算公式具体如下所示:
fitness=0.5ε1+0.5ε2
其中,fitness为适应度函数的函数值,ε1为平均磁位误差,ε2为平均磁感应强度分布误差。
与现有技术相比,本发明具有以下有益效果:
本发明通过模拟电流的方程组计算得到平均磁位误差和平均磁感应强度分布误差,将包含平均磁位误差和平均磁感应强度分布误差的适应度函数作为粒子群算法的目标函数来寻优,根据最适应架空输电线路的最终模拟电流来计算架空输电线路下方的磁感应强度分布,提高了电网建设过程中磁感应强度分布计算结果的准确性,同时降低了磁感应强度分布的计算量和计算成本,缩短了计算用时,有效解决了普通模拟电流法中,模拟电流的个数和位置只能依靠经验来寻找的问题。
附图说明
图1为本发明分裂导线的结构示意图;
图2为本发明初始模拟电流、匹配点和校验点在导线截面分布示意图;
图3为本发明无限长直导线产生的磁位和磁感应强度示意图;
图4为本发明磁感应强度计算结果与传统模拟电流法计算结果的对比图
图5为本发明普通模拟电流法的流程示意图;
图6为本发明粒子群算法的流程示意图;
图7为本发明应用在架空输电线路下方的磁感应强度计算的流程示意图。
具体实施方式
下面结合附图和具体实施例对本发明进行详细说明。本实施例以本发明技术方案为前提进行实施,给出了详细的实施方式和具体的操作过程,但本发明的保护范围不限于下述的实施例。
实施例
如图7所示,一种基于架空输电线路的二维磁感应强度控制方法,具体包括以下步骤:
S1、获取架空输电导线的三相导线中分裂导线的参数,计算分裂导线的等效半径;
S2、获取三相导线中初始模拟电流的个数和位置,并在三相导线的表面设置匹配点和校验点;
S3、获取预设的零磁位参考点,通过构建方程组计算得到初始模拟电流的电流值和匹配点的磁位,根据初始模拟电流的电流值计算出校验点的磁位,根据匹配点的磁位和校验点的磁位计算得到平均磁位误差;
S4、根据初始模拟电流的电流值计算出架空输电导线下方的第一磁感应强度分布,并获取有限元法计算得到的第二磁感应强度分布,第一磁感应强度分布与第二磁感应强度分布对比得到平均磁感应强度分布误差;
S5、根据平均磁位误差和平均磁感应强度分布误差计算适应度函数的函数值,如图6所示,通过粒子群算法进行优化得到函数值最小时对应的初始模拟电流的个数和位置,作为最终模拟电流的个数和位置,按照如图5所示的方法计算出每个最终模拟电流的数值;
S6、根据每个最终模拟电流的数值,通过磁感应强度的计算公式计算出各个最终模拟电流产生的磁感应强度,并进行叠加,得到架空输电线路下方的最终磁感应强度分布,如图4所示,相较于普通的模拟电流法按照经验设置模拟电流的个数和位置,计算得到的磁感应强度分布与有限元计算的结果更加贴近。
零磁位参考点具体为三相导线中A相导线表面的匹配点。
三相导线中每相导线的初始模拟电流的位置分别位于相应的圆周上。
圆周的半径小于三相导线的半径,圆周的圆心和三相导线的圆心重合。
分裂导线如图1所示,等效半径的计算公式如下所示:
其中,req为等效半径,R为分裂导线的半径,n为分裂导线的根数,r为次导线的半径。
如图2所示,三相导线中每相导线上匹配点、校验点和初始模拟电流的数量相同,三相导线之间的匹配点、校验点和初始模拟电流的数量相同。
如图3所示,匹配点的磁位的计算公式如下所示:
步骤S3中方程组根据磁位系数、磁感应强度和磁感应强度系数进行构建,磁位系数的计算公式如下所示:
磁感应强度的计算公式如下所示:
磁感应强度系数的计算公式如下所示:
步骤S3中构建的方程组具体如下所示:
其中,Ij为初始模拟电流的电流值,rij为第j个初始模拟电流到第i个匹配点的距离,r1j为第j个模拟电流到磁位参考点的距离,磁位参考点即为第一个匹配点,II、Iп和Iш分别为三相导线A相电流、B相电流和C相电流的大小,N为三相导线中单相导线上匹配点的数量,为B相导线表面匹配点的磁位,为C相导线表面匹配点的磁位。
适应度函数为粒子群算法的目标函数,适应度函数的计算公式具体如下所示:
fitness=0.5ε1+0.5ε2
其中,fitness为适应度函数的函数值,ε1为平均磁位误差,ε2为平均磁感应强度分布误差。
本实施例中,等效导线中设置模拟电流的个数和位置,常规的做法是按照经验来设置,一般是在每根分解导线的等效导线的几何中心设置模拟电流,这样的计算虽然比较简单,但是计算出的磁感应强度误差较大。本发明使用粒子群算法,在等效导线横截面上,对构造出的适应度函数寻找模拟电流的个数和位置的最优值,然后再用得到的模拟电流最优个数和位置来对磁感应强度进行求解,提高了磁感应强度计算结果的准确性。
此外,需要说明的是,本说明书中所描述的具体实施例,所取名称可以不同,本说明书中所描述的以上内容仅仅是对本发明结构所做的举例说明。凡依据本发明构思的构造、特征及原理所做的等效变化或者简单变化,均包括于本发明的保护范围内。本发明所属技术领域的技术人员可以对所描述的具体实例做各种各样的修改或补充或采用类似的方法,只要不偏离本发明的结构或者超越本权利要求书所定义的范围,均应属于本发明的保护范围。
Claims (10)
1.一种基于架空输电线路的二维磁感应强度控制方法,其特征在于,具体包括以下步骤:
S1、获取架空输电导线的三相导线中分裂导线的参数,计算分裂导线的等效半径;
S2、获取三相导线中初始模拟电流的个数和位置,并在三相导线的表面设置匹配点和校验点;
S3、获取预设的零磁位参考点,通过构建方程组计算得到初始模拟电流的电流值和匹配点的磁位,根据所述初始模拟电流的电流值计算出校验点的磁位,根据匹配点的磁位和校验点的磁位计算得到平均磁位误差;
S4、根据初始模拟电流的电流值计算出架空输电导线下方的第一磁感应强度分布,并获取有限元法计算得到的第二磁感应强度分布,第一磁感应强度分布与第二磁感应强度分布对比得到平均磁感应强度分布误差;
S5、根据平均磁位误差和平均磁感应强度分布误差计算适应度函数的函数值,通过粒子群算法进行优化得到函数值最小时对应的初始模拟电流的个数和位置,作为最终模拟电流的个数和位置,并计算出每个最终模拟电流的数值;
S6、根据每个最终模拟电流的数值,计算出各个最终模拟电流产生的磁感应强度,并进行叠加,得到架空输电线路下方的最终磁感应强度分布。
2.根据权利要求1所述的一种基于架空输电线路的二维磁感应强度控制方法,其特征在于,所述零磁位参考点具体为三相导线中A相导线表面的匹配点。
3.根据权利要求1所述的一种基于架空输电线路的二维磁感应强度控制方法,其特征在于,所述三相导线中每相导线的初始模拟电流的位置分别位于相应的圆周上。
4.根据权利要求3所述的一种基于架空输电线路的二维磁感应强度控制方法,其特征在于,所述圆周的半径小于三相导线的半径,所述圆周的圆心和三相导线的圆心重合。
6.根据权利要求1所述的一种基于架空输电线路的二维磁感应强度控制方法,其特征在于,所述三相导线中每相导线上匹配点、校验点和初始模拟电流的数量相同,三相导线之间的匹配点、校验点和初始模拟电流的数量相同。
10.根据权利要求1所述的一种基于架空输电线路的二维磁感应强度控制方法,其特征在于,所述适应度函数的计算公式具体如下所示:
fitness=0.5ε1+0.5ε2
其中,fitness为适应度函数的函数值,ε1为平均磁位误差,ε2为平均磁感应强度分布误差。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110131901.4A CN112765897B (zh) | 2021-01-30 | 2021-01-30 | 一种基于架空输电线路的二维磁感应强度控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110131901.4A CN112765897B (zh) | 2021-01-30 | 2021-01-30 | 一种基于架空输电线路的二维磁感应强度控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112765897A CN112765897A (zh) | 2021-05-07 |
CN112765897B true CN112765897B (zh) | 2022-07-05 |
Family
ID=75704236
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110131901.4A Active CN112765897B (zh) | 2021-01-30 | 2021-01-30 | 一种基于架空输电线路的二维磁感应强度控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112765897B (zh) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103714197A (zh) * | 2013-11-18 | 2014-04-09 | 重庆大学 | 一种优化超/特高压输电线路电磁环境的结构设计方法 |
CN109902432A (zh) * | 2019-03-14 | 2019-06-18 | 国网上海市电力公司 | 一种基于高压输电线路周围的工频磁场获取方法 |
CN110738009A (zh) * | 2019-10-14 | 2020-01-31 | 山东科技大学 | 一种输电线路电场计算中导线内模拟电荷的设置方法 |
CN110927484A (zh) * | 2019-11-15 | 2020-03-27 | 国网江西省电力有限公司电力科学研究院 | 一种交流输电线路无线电干扰三维分布计算方法及系统 |
CN111709150A (zh) * | 2019-08-07 | 2020-09-25 | 电子科技大学 | 任意形状磁体磁场空间分布仿真模拟方法 |
-
2021
- 2021-01-30 CN CN202110131901.4A patent/CN112765897B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103714197A (zh) * | 2013-11-18 | 2014-04-09 | 重庆大学 | 一种优化超/特高压输电线路电磁环境的结构设计方法 |
CN109902432A (zh) * | 2019-03-14 | 2019-06-18 | 国网上海市电力公司 | 一种基于高压输电线路周围的工频磁场获取方法 |
CN111709150A (zh) * | 2019-08-07 | 2020-09-25 | 电子科技大学 | 任意形状磁体磁场空间分布仿真模拟方法 |
CN110738009A (zh) * | 2019-10-14 | 2020-01-31 | 山东科技大学 | 一种输电线路电场计算中导线内模拟电荷的设置方法 |
CN110927484A (zh) * | 2019-11-15 | 2020-03-27 | 国网江西省电力有限公司电力科学研究院 | 一种交流输电线路无线电干扰三维分布计算方法及系统 |
Non-Patent Citations (4)
Title |
---|
1000kV特高压交流输电...架建筑的感应电压、电流研究;徐静;《中国优秀硕士学位论文全文数据库》;20200915;全文 * |
Interference-rejecting current measurement method with tunnel magnetoresistive magnetic sensor array;Yafeng Chen, Qi Huang, Arsalan Habib Khawaja;《IET SCIENCE MEASUREMENT & TECHNOLOGY》;20180822;全文 * |
UHVAC同塔双回输电线路表面场强分析;马爱清;《电力系统及其自动化学报》;20190630;全文 * |
超高压输电线下电磁场分布状况分析及检测技术研究;姜志超;《中国优秀硕士学位论文全文数据库》;20111215;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN112765897A (zh) | 2021-05-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
KR101219545B1 (ko) | 전력계통에서의 최적화 기법을 적용한 파라미터 추정 방법 | |
CN112907075A (zh) | 一种电力系统综合负荷模型参数辨识方法 | |
CN109299564B (zh) | 变压器偏磁电流仿真计算过程中温度因素影响的修正方法 | |
CN109659934B (zh) | 基于短路电流裕度的新能源并网规划优选方法 | |
CN112765897B (zh) | 一种基于架空输电线路的二维磁感应强度控制方法 | |
CN110188480A (zh) | 一种直流偏磁条件下铁磁材料的磁滞特性模拟分析系统和方法 | |
CN108376208B (zh) | 一种接地网阴极保护的辅助阳极系统优化方法 | |
CN111313406B (zh) | 一种集电线路的确定方法、装置及电站 | |
CN111859714A (zh) | 工频电场强度计算方法、系统及工频电场屏蔽装置 | |
CN116742632A (zh) | 一种配电网最大承载力估算方法及系统 | |
CN116151068A (zh) | 一种多绕组高频变压器绝缘结构优化设计方法 | |
CN115714421A (zh) | 基于智能单粒子优化算法的分布式电源选址与定容方法 | |
CN105514981A (zh) | 一种基于建筑信息模型的光伏逆变器接入位置优化方法 | |
Shen et al. | Study on electromagnetic environment and conductor selection for 750 kv same-tower four-circuit lines | |
CN102254071A (zh) | 交流特高压单回紧凑型10分裂导线相间距离确定方法 | |
CN112347676A (zh) | 电机定子绕组损耗的快速计算方法、装置和系统 | |
CN112834022A (zh) | 一种交流输电线路交叉跨越区域可听噪声的计算方法 | |
CN117368648B (zh) | 配电网单相接地故障检测方法、系统、设备和存储介质 | |
CN111859713B (zh) | 一种户内变电站工频电场优化装置 | |
CN113486488A (zh) | 一种考虑天气条件的交流输电线路电磁环境计算平台 | |
CN116596199B (zh) | 一种配电网线损影响因子的叠加分析方法和系统 | |
Daniels et al. | Mathematical Modeling of a Electric Power Network Segment of Ijuí using the Pi Model | |
CN110427631B (zh) | 理论线损计算所需主网线路与变压器参数双重校核方法 | |
CN112134288B (zh) | 基于前/后向扫描谐波潮流的谐波污染配电网络重构方法 | |
CN116365555B (zh) | 一种计及局部三相不对称的电网状态估计方法和系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
TR01 | Transfer of patent right |
Effective date of registration: 20240709 Address after: Room 2307 and 2309, No.1, Fenghe Road, China (Shanghai) pilot Free Trade Zone, Pudong New Area, Shanghai, 200120 Patentee after: SHANGHAI HAIDA COMMUNICATION CO.,LTD. Country or region after: China Address before: 201620 No. 333, Longteng Road, Shanghai, Songjiang District Patentee before: SHANGHAI University OF ENGINEERING SCIENCE Country or region before: China |
|
TR01 | Transfer of patent right |