CN106155980B - 基于不均匀激励电场的架空线路场线耦合建模方法 - Google Patents
基于不均匀激励电场的架空线路场线耦合建模方法 Download PDFInfo
- Publication number
- CN106155980B CN106155980B CN201510152893.6A CN201510152893A CN106155980B CN 106155980 B CN106155980 B CN 106155980B CN 201510152893 A CN201510152893 A CN 201510152893A CN 106155980 B CN106155980 B CN 106155980B
- Authority
- CN
- China
- Prior art keywords
- transmission line
- overhead transmission
- field
- electric field
- method based
- 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
- 230000005540 biological transmission Effects 0.000 title claims abstract description 37
- 230000005684 electric field Effects 0.000 title claims abstract description 32
- 238000000034 method Methods 0.000 title claims abstract description 25
- 230000008878 coupling Effects 0.000 title claims abstract description 15
- 238000010168 coupling process Methods 0.000 title claims abstract description 15
- 238000005859 coupling reaction Methods 0.000 title claims abstract description 15
- 230000005284 excitation Effects 0.000 title claims abstract description 14
- 238000005259 measurement Methods 0.000 claims abstract description 8
- 239000013598 vector Substances 0.000 claims description 6
- 238000004422 calculation algorithm Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 239000004020 conductor Substances 0.000 description 2
- 230000005611 electricity Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000005855 radiation Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 230000003321 amplification Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000005474 detonation Methods 0.000 description 1
- 238000004880 explosion Methods 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 239000004615 ingredient Substances 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Landscapes
- Measurement Of Resistance Or Impedance (AREA)
Abstract
本发明涉及基于不均匀激励电场的架空线路场线耦合建模方法，所述方法包括，建立坐标；根据线路参数构建架空线路模型，获得架空线路的电报方程；将测量得到的非均匀水平电场的场强幅值进行拟合，生成场强幅值与距离x的函数表达式；代入架空线路的电报方程，获得在非均匀波入射下的线缆响应；该方法不但提高了计算精度，也使计算效率大大提升。
Description
技术领域
本发明涉及一种电力系统领域的建模方法，具体涉及基于不均匀激励电场的架空线路场线耦合建模方法。
背景技术
核爆炸产生的高空核爆电磁脉冲(HEMP)包含的E1成分频谱覆盖范围包括中频、高频、高频甚至一些超高频波段的信号，具有辐射范围广、峰值场强高和频谱宽的特点。此外，HEMP在输电线路等多导体传输线上，可以感应出严重威胁电力系统的幅值较高的感应电压和电流。因此对架空输电线路进行建模，进而对其HEMP响应进行仿真计算一直是研究的热点。
经典的场线耦合计算模型一般针对均匀场入射的情况，而在实际情况中，线缆上的激励场大多为非均匀场，这对线缆响应计算造成了一定程度的困难。若仍旧按照经典的场线耦合计算模型对线缆进行建模，势必会对计算结果带来较大误差。
针对激励电场分布不均的情况，现有技术中一般采用时域FDTD算法解决非均匀平面波入射下的架空线路电磁脉冲响应仿真计算；该算法分别在时间和空间上离散，空间上的波形分布不均匀，只需利用离散化空间上的电场值来进行近似计算即可，虽然该方法具有较好的适用性，但仍存在以下不足：
1.由于FDTD算法是基于离散化的时间和空间进行计算，若入射波在空间上的分布越不均匀，则计算精度越低；
2.FDTD方法是一种时域的建模方法，而在频域模型中，大地阻抗一般为频变参数，利用时域模型处理频变参数较为困难，极大地影响计算效率。
发明内容
为了弥补现有技术的不足，本发明提供一种基于不均匀激励电场的场线耦合的建模方法；通过构建架空线路模型，获得架空线路的电报方程，并将测量得到的非均匀水平电场的场强幅值进行拟合，获得在非均匀波入射下的线缆响应；由此克服了在非均匀波入射情况下，计算精度低和计算效率差的缺陷。
本发明的目的是采用下述技术方案实现的：
基于不均匀激励电场的架空线路场线耦合建模方法，所述方法包括下述步骤：
(1)建立坐标；
(2)根据线路参数构建架空线路模型，获得架空线路的电报方程；
(3)拟合测量得到的非均匀水平电场的场强幅值，生成场强幅值与距离x的函数表达式；
(4)代入架空线路的电报方程，获得非均匀波入射下的线缆响应。
优选的，所述步骤(1)包括：在x轴上选取x1～xk的k个节点：所述x1～xk的k个节点范围内存在无线缆时，对应的场强幅值为E1～Ek。
优选的，所述步骤(2)包括，地面设置长为L、高为h、首端负载Z1、末端负载为ZL的架空导线，根据Agrawal频域电报方程得到如下式所示的架空线路电报方程：
其中，x为线缆上距离首端的距离；
边界条件为:
进一步地，解联立式(1)和式(2)得通解，其通解形式由方程的激励源Ex决定；其中，所述Ex为沿线分布的非均匀水平电场，Vs为沿线的散射电压，I为沿线电流，Z’和Y’为分布阻抗，Z1和ZL分别为线缆首末端的负载，z为垂直方向坐标，Ez inc为外界电场的垂直分量。
优选的，所述步骤(3)中，用一组正交基函数拟合不同位置处的Ex(x)场强幅值按照傅里叶级数形式，其表达式为：
式(3)中，w和n皆为拟合得到的参数，分别代表三角函数的基波和谐波次数；x为不同位置；
设
拟合后为：
其中，a0～an，b1～bn，C-n～Cn皆为待定系数；j指虚部常数，k为谐波次数n的上限。
优选的，所述步骤(4)中，依次将n个向量作为激励源带入架空线路电报方程，得到每个向量的时域响应，将n个响应进行叠加。
与最接近的现有技术比，本发明提供的技术方案的有益效果是：
本发明提供的技术方案根据Agrawal频域电报方程获得架空线路的电报方程，使得在处理频变参数时尤为方便，提升了计算效率。
利用正交基向量对沿线的电场强度幅值进行拟合，获得沿线电场强度随距离变化的函数表达式，并将其带入架空线路的电报方程；这种利用函数曲线对离散的实测电场进行拟合的方法，极大提高了计算精度。
附图说明
图1为本发明提供的基于不均匀激励电场的架空线路场线耦合建模方法流程图；
图2为本发明提供的地面上的导线结构示意图；
图3为本发明实施例的水平场分布测量示意图；
图4为本发明实施例的水平场幅值水平分布拟合示意图。
具体实施方式
下面结合附图对本发明的具体实施方式作进一步的详细说明。
如图1所示，基于不均匀激励电场的架空线路场线耦合建模方法，所述方法包括下述步骤：
(1)建立坐标；所述步骤(1)包括：在x轴上选取x1～xk的k个节点：所述x1～xk的k个节点范围内存在无线缆时，对应的场强幅值为E1～Ek。
(2)根据线路参数构建架空线路模型，获得架空线路的电报方程；如图2所示，所述步骤(2)包括，地面设置长为L、高为h、首端负载Z1、末端负载为ZL的架空导线，根据Agrawal频域电报方程得到如下式所示的架空线路电报方程：
其中，x为线缆上距离首端的距离；
边界条件为:
解联立式(1)和式(2)得通解，其通解形式由方程的激励源Ex决定；其中，所述Ex为沿线分布的非均匀水平电场，Vs为沿线的散射电压，I为沿线电流，Z’和Y’为分布阻抗，Z1和ZL分别为线缆首末端的负载，z为垂直方向坐标，Ez inc为外界电场的垂直分量。
(3)拟合测量得到的非均匀水平电场的场强幅值，生成场强幅值与距离x的函数表达式；所述步骤(3)中，用一组正交基函数拟合不同位置处的Ex(x)场强幅值按照傅里叶级数形式，其表达式为：
式(3)中，w和n皆为拟合得到的参数，分别代表三角函数的基波和谐波次数；x为不同位置；
设
拟合后为：
其中，a0～an，b1～bn，C-n～Cn皆为待定系数；j指虚部常数，k为谐波次数n的上限。
(4)代入架空线路的电报方程，获得非均匀波入射下的线缆响应。
所述步骤(4)中，依次将n个向量作为激励源带入架空线路电报方程，得到每个向量的时域响应，将n个响应进行叠加。
实施例：选取一有界波模拟器作为辐射源，首先选取了其内部每隔1m处的几点作为电场测点，测量了其水平电场分布，测量示意图如图3所示：
按照前述拟合方法对几个测点处的电场幅值进行拟合，拟合结果为：
|Ex|＝2.83e(-0.99x)+0.84e(0.21x) (6)
拟合结果如图4所示，拟合曲线较好的反应了真实测量结果。
最后应当说明的是：以上实施例仅用以说明本发明的技术方案而非对其限制，尽管参照上述实施例对本发明进行了详细的说明，所属领域的普通技术人员应当理解：依然可以对本发明的具体实施方式进行修改或者等同替换，而未脱离本发明精神和范围的任何修改或者等同替换，其均应涵盖在本发明的权利要求范围当中。
Claims (5)
1.基于不均匀激励电场的架空线路场线耦合建模方法，其特征在于，所述方法包括下述步骤：
(1)建立坐标；
(2)根据线路参数构建架空线路模型，获得架空线路的电报方程；
(3)拟合测量得到的非均匀水平电场的场强幅值，生成场强幅值与距离x的函数表达式；
(4)代入架空线路的电报方程，获得非均匀波入射下的线缆响应。
2.如权利要求1所述的基于不均匀激励电场的架空线路场线耦合建模方法，其特征在于，所述步骤(1)包括：在x轴上选取x1～xk的k个节点：所述x1～xk的k个节点范围内存在无线缆时，对应的场强幅值为E1～Ek。
3.如权利要求1所述的基于不均匀激励电场的架空线路场线耦合建模方法，其特征在于，所述步骤(2)包括，地面设置长为L、高为h、首端负载Z1、末端负载为ZL的架空导线，根据Agrawal频域电报方程得到如下式所示的架空线路电报方程：
其中，x为线缆上距离首端的距离；
边界条件为:
解联立式(1)和式(2)得通解，其通解形式由方程的激励源Ex决定；其中，所述Ex为沿线分布的非均匀水平电场，Vs为沿线的散射电压，I为沿线电流，Z’和Y’为分布阻抗，Z1和ZL分别为线缆首末端的负载，z为垂直方向坐标，Ez inc为外界电场的垂直分量。
4.如权利要求1所述的基于不均匀激励电场的架空线路场线耦合建模方法，其特征在于，所述步骤(3)中，用一组正交基函数拟合不同位置处的Ex(x)场强幅值按照傅里叶级数形式，其表达式为：
式(3)中，w和n皆为拟合得到的参数，分别代表三角函数的基波和谐波次数；x为不同位置；
设
拟合后为：
其中，a0～an，b1～bn，C-n～Cn皆为待定系数；j指虚部常数，k为谐波次数n的上限。
5.如权利要求4所述的基于不均匀激励电场的架空线路场线耦合建模方法，其特征在于，所述步骤(4)中，依次将n个向量作为激励源带入架空线路电报方程，得到每个向量的时域响应，将n个响应进行叠加。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510152893.6A CN106155980B (zh) | 2015-04-01 | 2015-04-01 | 基于不均匀激励电场的架空线路场线耦合建模方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510152893.6A CN106155980B (zh) | 2015-04-01 | 2015-04-01 | 基于不均匀激励电场的架空线路场线耦合建模方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106155980A CN106155980A (zh) | 2016-11-23 |
| CN106155980B true CN106155980B (zh) | 2019-03-22 |
Family
ID=57337888
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510152893.6A Active CN106155980B (zh) | 2015-04-01 | 2015-04-01 | 基于不均匀激励电场的架空线路场线耦合建模方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106155980B (zh) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107608931B (zh) * | 2017-08-07 | 2020-09-22 | 三峡大学 | 一种基于Chebyshev拟谱-二级边界值耦合方法的传输线瞬态响应快速计算方法 |
| CN107895078A (zh) * | 2017-11-09 | 2018-04-10 | 西安交通大学 | 强电磁脉冲对长电力输电线路耦合响应的宏模型构建方法 |
| CN111797541B (zh) * | 2020-07-20 | 2022-08-26 | 吉林大学 | 一种场线耦合不确定性量化及全局灵敏度计算方法 |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103714239A (zh) * | 2013-11-15 | 2014-04-09 | 贵州电力试验研究院 | 雷击大地时低压线路绝缘子雷电感应电压计算方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101981774B (zh) * | 2008-04-03 | 2014-06-25 | 西门子公司 | 用于产生故障信号的方法、装置和现场设备 |
-
2015
- 2015-04-01 CN CN201510152893.6A patent/CN106155980B/zh active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103714239A (zh) * | 2013-11-15 | 2014-04-09 | 贵州电力试验研究院 | 雷击大地时低压线路绝缘子雷电感应电压计算方法 |
Non-Patent Citations (4)
| Title |
|---|
| 《Lightning-Induced Voltages on Overhead Lines-Application of the Extended Rusck Model》;Alexandre Piantini 等;《IEEE Transactions on Electromagnetic Compatibility》;20090831;第51卷(第3期);全文 |
| 《Parameters on Lightning-Induced Voltages on Overhead Lines》;Albert De Conti 等;《Proc.of the 2014 International Symposium on Electromagnetic Compatibility》;20141231;全文 |
| 《架空多导体传输线HEMP响应建模及分析》;李妮 等;《现代应用物理》;20141231;第5卷(第4期);全文 |
| 《架空配电线路雷电感应过电压计算研究》;边凯 等;《中国电机工程学报》;20121105;第32卷(第31期);全文 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106155980A (zh) | 2016-11-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103869173B (zh) | 一种测量地表到地深数十公里土壤电阻率分布的方法 | |
| CN105953762B (zh) | 一种用于净空距离测量的输电线路建模方法和测距装置 | |
| CN106155980B (zh) | 基于不均匀激励电场的架空线路场线耦合建模方法 | |
| CN103412242B (zh) | 一种基于快速独立分量分析和互信息的谐波源定位方法 | |
| CN103792433A (zh) | 以火花系数修正杆塔接地装置的低幅值冲击电阻测量方法 | |
| US9746498B2 (en) | System and method for monitoring a power line | |
| CN108959768A (zh) | 基于分区域模拟电荷法的复杂场景工频电场计算方法 | |
| CN106018970B (zh) | 一种双层土壤中接地网接地电阻的测量方法 | |
| Gatta et al. | Generalized pi-circuit tower grounding model for direct lightning response simulation | |
| Zhang et al. | An improved approach for modeling lightning transients of wind turbines | |
| Lima et al. | Transient grounding impedance and transient resistivity measurements using a very short current lead | |
| CN109458988A (zh) | 基于uwb无线测距技术的抱杆倾角测量方法 | |
| US9664710B2 (en) | Sensory assembly system and method | |
| CN102590686A (zh) | 一种双极直流输电线路无线电干扰确定方法 | |
| CN105510733A (zh) | 一种高压、特高压直流输电线路的参数测量方法 | |
| Xie et al. | Altitude correction of radio interference of HVDC transmission lines part I: converting method of measured data | |
| Wang et al. | Calculation of current distribution in the lightning protective system of a residential house | |
| CN105260583B (zh) | 一种特高压工频电磁场对人体生物效应的计算方法及系统 | |
| CN105486961B (zh) | 基于电晕笼的直流输电线路无线电干扰激发函数测试方法 | |
| Araneo et al. | Hybrid and pi-circuit approaches for grounding system lightning response | |
| CN110927484A (zh) | 一种交流输电线路无线电干扰三维分布计算方法及系统 | |
| CN104778286B (zh) | 掠海飞行器电磁散射特性快速仿真方法 | |
| CN113297526B (zh) | 一种基于Wenner四极和大地电磁数据的水平分层土壤结构联合反演方法 | |
| De Silva et al. | Modeling and measuring of antenna array S-parameters and radiation efficiency | |
| CN203164295U (zh) | 一种用于测量输电线路铁塔输入阻抗和谐振频率的装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| TA01 | Transfer of patent application right | ||
| TA01 | Transfer of patent application right |
Effective date of registration: 20170427 Address after: 100192 Beijing city Haidian District Qinghe small Camp Road No. 15 Applicant after: China Electric Power Research Institute Applicant after: State Grid Corporation of China Applicant after: Xi'an Jiaotong University Applicant after: State Grid Zhejiang Electric Power Company Address before: 100031 Xicheng District West Chang'an Avenue, No. 86, Beijing Applicant before: State Grid Corporation of China Applicant before: China Electric Power Research Institute Applicant before: Xi'an Jiaotong University |
|
| GR01 | Patent grant | ||
| GR01 | Patent grant |