《高压交流架空送电线路对埋地金属油气管道的干扰确定方法.pdf》由会员分享，可在线阅读，更多相关《高压交流架空送电线路对埋地金属油气管道的干扰确定方法.pdf（9页珍藏版）》请在专利查询网上搜索。
1、(10)申请公布号 CN 102435877 A (43)申请公布日 2012.05.02 CN 102435877 A *CN102435877A* (21)申请号 201110275055.X (22)申请日 2011.09.16 G01R 31/00(2006.01) (71)申请人 国网电力科学研究院 地址 430074 湖北省武汉市洪山区珞喻路 143 号 (72)发明人 万保权 谢辉春 张建功 刘兴发 (74)专利代理机构 武汉开元知识产权代理有限 公司 42104 代理人 潘杰 (54) 发明名称 高压交流架空送电线路对埋地金属油气管道 的干扰确定方法 (57) 摘要 本发明涉及。
2、一种高压交流架空送电线路对埋 地金属管道干扰计算方法， 所采用的方法是 ： 通 过对送电线路与管道理想以及非理想平行时的 感性耦合进行简化计算推导， 将以往计算模型需 求解的复杂微分方程模型简化为多元一次方程组 模型。本发明优势在于使得复杂微分方程简化 为多元一次方程组， 编程简单， 且大大提高了计 算速度及对计算内存资源的占用， 利用 Matlab、 Mathematica 等数学计算软件很容易求解该方程 组， 通过求解简单的多元一次方程组可以得到管 道上的对地电压、 电流分布。普通技术人员均可 使用该项发明， 具有较好的实用性。需要补充的 是， 使用该简化方法时， 由于采用的是等效方法， 。
3、因此管道细分段越多， 越接近平行段情况， 计算越 准确。 (51)Int.Cl. (19)中华人民共和国国家知识产权局 (12)发明专利申请 权利要求书 2 页 说明书 4 页 附图 2 页 CN 102435883 A1/2 页 2 1. 一种高压交流架空送电线路对埋地金属油气管道的干扰确定方法， 所采用的方法 是 ： 通过对送电线路与管道理想以及非理想平行时的感性耦合进行简化计算推导， 将以往 计算模型需求解的复杂微分方程模型简化为多元一次方程组模型， 其特征在于 ： 该方法的 步骤是 ： 第一步骤 ： 获得高压交流架空送电线路导线参数、 运行电压及电流参数、 埋地金属管道 特性参数 ； 。
4、第二步骤 ： 获得高压交流架空送电线路、 埋地金属管道三维坐标， 包括两者的走向、 两 者相对位置关系 ； 第三步骤 ： 利用计算公式计算高压交流架空送电线路对埋地金属管道干扰程度。 2. 如权利要求 1 所述的高压交流架空送电线路对埋地金属油气管道的干扰确定方法， 其特征在于 ： 第三步骤中送电线路与管道理想以及非理想平行时的感性耦合进行简化计算 推导， 将以往计算模型需求解的复杂微分方程模型简化为多元一次方程组模型， 从而获得 高压交流架空送电线路对埋地金属管道干扰程度。 3. 如权利要求 1 所述的高压交流架空送电线路对埋地金属油气管道的干扰确定方法， 其特征在于 ： 该方法包括 1)、。
5、 送电线路与管道理想平行时感性耦合的确定方法 ： 当送电线路正常运行时， 对管道的干扰主要是相线交变电流产生的磁场在两者互感 M 的作用下于金属管道上产生纵向电动势 E， 计算时考虑交流送电线路与管道并行， L 为并 行长度， 管道所经过区域土壤电阻率保持不变， 其本身特性参数也不变， 考察管道基本微元 段， 依据传输线理论列写电压、 电流回路方程如下 ： Edx -U+IZ dx+(U+dU) IZ dx+dU (1) 其中， E 为交流送电线路在管道基本微元段感性耦合产生的纵向电动势， I 为管道基本 微元段感性耦合产生的电流， 其余参数为管道特性参数， 对于 (1) 式中的 E， 按照 。
6、E -jMIphase (3) 式中， 负号表示方向， 与送电线路电流方向相反 ； M送电线路与管道间的互感系数， 在理想并行段可看作常数 ； Iphase送电线路运行电流 ； 其中 式中 ： a 为架空影响线和地下被影响线之间的距离 ； b 为架空影响线的架设高度 ； c 为 地下被影响线的埋深 ； 为空气导磁率 ； 为大地电导率， 将式 (4) 做如下简化 ： 权 利 要 求 书 CN 102435877 A CN 102435883 A2/2 页 3 通过求解上述方程 (1)-(5) 即可获得送电线路与管道理想平行时感性耦合干扰 ； 2)、 送电线路与管道非理想平行时感性耦合简化计算方法。
7、 ： 对于斜接近段 P1-P2 的管道， 不同位置的感应纵电动势 E 与送电线路和管道间的互感 有关， 因此不同位置的 E 不同， 即使用平均间距代替进行互感系 数的求解， 其中 a1、 a2 为 P1、 P2 点与送电线路的最短间距， 此时将斜接近段等效为平行段， 根据式 (3)， 得到 E 近似为常数， 进而因此在斜接近段 P1-P2 管道上的电压、 电流分 布有如下形式解， Up1-2 -ZA1exp(x)-B1exp(-x) (6) Up1-2P1P2段电压分布， x 为从 P1到待求点的距离 ； Ip1-2P1P2段电流分布。A1、 B1P1P2段待定系数 ； 将连续曲折的管道细分，。
8、 利用等效方法将斜接近段、 交叉段转化为平行段， 结合土壤电 阻率及平行间距查找诺模互感系数图、 采用简化计算公式得到互感系数， 进而求得曲折段 的近似等效单位感应纵电动势， 接着利用单回传输线方程的近似解， 得到每个细分段的电 压、 电流表达式 ； 对于同一端点而言， 其即为前一细分段的末端点， 也为后一细分段的起始 点， 因此再根据同一点电压电流值唯一得到一个规则关于待定系数 Ai、 Bi的多元一次方程 组， 联合管道首末端边界条件， 求解该多元一次方程组， 求出各细分段管道对地电压、 电流 分布， 从而获得送电线路与管道非理想平行时感性耦合干扰。 权 利 要 求 书 CN 1024358。
9、77 A CN 102435883 A1/4 页 4 高压交流架空送电线路对埋地金属油气管道的干扰确定方 法 技术领域 0001 本发明专利属于属于交流输变电工程电磁环境领域， 具体地讲是一种高压交流架 空送电线路对埋地金属油气管道的干扰确定方法。 背景技术 0002 为了满足我国经济社会可持续发展的用电需求， 建设以特高压电网为核心的坚强 电网已成为国家电力建设的战略目标， 与此同时由于能源需求的飞速增长， 输油输气管道 建设也以前所未有的速度在全国各地展开。由于我国土地资源有限， 高压送电线路和输油 输气管道经常面临同走廊并行架设和交叉跨越的状况。高压送电线路与埋地金属管道如 果长距离接近。
10、， 可能产生交流电磁影响主要是长期存在着的线路运行电流通过空间磁场耦 合， 在邻近的金属管道上感应出纵向电动势， 可能影响管道维护人员的安全作业， 并在管道 防腐层引起泄漏电流从而加速管道表面腐蚀。因此， 如何准确快速评估高压送电线路与沿 线邻近的埋地输油输气管道间电磁干扰的是迫切需要解决的问题。 0003 中国专利 “200810197006.7” 采用电磁场理论计算并构造特高压交流送电线路和 埋地金属管道所组成的多导体传输线系统等效电路模型， 通过广义双侧消去法或矩量法计 算送电线路通过感性耦合对管道的干扰影响， 分析线路通过感性耦合及阻性耦合对管道的 干扰影响， 最后综合上述技术提出有效。
11、的防护措施。该专利重点在于对干扰的提供了一种 评估方法， 而其中对于感性耦合的计算技术， 采用的广义双侧消去法或矩量法， 程序编制复 杂， 对非专业人员而言， 存在较大的难度， 可操作性不够。 发明内容 0004 本发明的目的是通过对送电线路与管道理想以及非理想平行时的感性耦合进行 简化计算推导， 将以往计算模型需求解的复杂微分方程模型简化为多元一次方程组模型， 而提供一种高压交流架空送电线路对埋地金属油气管道的干扰确定方法。 0005 为了实现上述目的， 本发明所采用的技术方案是 ： 0006 一种高压交流架空送电线路对埋地金属油气管道的干扰确定方法， 其方法是 ： 通 过对送电线路与管道理。
12、想以及非理想平行时的感性耦合进行简化计算推导， 将以往计算模 型需求解的复杂微分方程模型简化为多元一次方程组模型， 具体步骤为 ： 0007 第一步骤 ： 获得高压交流架空送电线路导线参数、 运行电压及电流参数、 埋地金属 管道特性参数 ； 0008 第二步骤 ： 获得高压交流架空送电线路、 埋地金属管道三维坐标， 包括两者的走 向、 两者相对位置关系 ； 0009 第三步骤 ： 利用计算技术计算高压交流架空送电线路、 对埋地金属管道干扰程度。 0010 本发明优势在于使得复杂微分方程简化为多元一次方程组， 编程简单， 且大大提 高了计算速度及对计算内存资源的占用， 利用 Matlab、 Ma。
13、thematica 等数学计算软件很容 说 明 书 CN 102435877 A CN 102435883 A2/4 页 5 易求解该方程组， 通过求解简单的多元一次方程组可以得到管道上的对地电压、 电流分布。 普通技术人员均可使用该项发明， 具有较好的实用性。需要补充的是， 使用该简化方法时， 由于采用的是等效方法， 因此管道细分段越多， 越接近平行段情况， 计算越准确。 附图说明 0011 图 1 为本发明送电线路与管道理想平行时的感性耦合等效电路图。 0012 图 2 为本发明受干扰的 “管道土壤” 回路基本微元段图。 0013 图 3 为本发明实际的送电线路与管道并行接近情况图。 00。
14、14 图 4 为本发明管道与送电线路斜接近图。 0015 图 5 为本发明实施例管道与送电线路相对位置关系图。 具体实施方式 0016 下面结合附图和实施例对本发明作进一步详细的说明。 0017 本发明包括送电线路与管道理想平行时感性耦合干扰的确定方法以及送电线路 与管道非理想平行时感性耦合干扰的确定方法。 0018 1、 送电线路与管道理想平行时感性耦合的确定方法 ： 0019 当送电线路正常运行时， 对管道的干扰主要是相线交变电流产生的磁场在两者互 感 M 的作用下于金属管道上产生纵向电动势 E。对于管道而言， 在进行干扰研究时可用如 下的电路等效 ( 图 1)， 其中 R 为管道分布电阻。
15、， L 为分布电感， G 为分布电导， C 为分布 电容。这些参数的获得可参照我国石油行业标准 SY/T0032-2000 的条文说明 A.3.3 计算得 到。 0020 计算方法基于传输线理论。计算时考虑交流送电线路 ( 单相 ) 与管道并行， L 为 并行长度， 管道所经过区域土壤电阻率保持不变， 其本身特性参数也不变。 考察管道基本微 元段， 如附图 2 ： 0021 据传输线理论列写电压、 电流回路方程如下 ： 0022 Edx -U+IZ dx+(U+dU) IZ dx+dU (1) 0023 0024 其中， E 为交流送电线路在管道基本微元段感性耦合产生的纵向电动势， I 为管道。
16、 基本微元段感性耦合产生的电流， Z 、 Y 为管道特性参数， 由 L 、 G 、 C 决定。 0025 对于 (1) 式中的 E， 可按照 0026 E -jMIphase (3) 0027 式中， 负号表示方向， 与送电线路电流方向相反 ； 0028 M送电线路与管道间的互感系数， 在理想并行段可看作常数 ； 0029 Iphase送电线路运行电流 ； 0030 角频率 ； 0031 其中关键参数为M， 计算比较复杂。 国际电话电报咨询委员会CCITT在防护导则中 给出架空影响线和地下被影响线之间互感系数基本表达式如下 ： 说 明 书 CN 102435877 A CN 102435883。
17、 A3/4 页 6 0032 0033 式中 ： a为架空影响线和地下被影响线之间的距离 ； b为架空影响线的架设高度 ； c 为地下被影响线的埋深 ； 0为空气导磁率 ； 为大地电导率。 0034 由于上述无穷积分的数值计算比较困难， 做如下简化 ： 0035 0036 式中 k 1.7811， 通过求解上述方程 (1)-(5) 即可获得送电线路与管道理想平行 时感性耦合干扰。 0037 2、 送电线路与管道非理想平行时感性耦合简化计算方法 ： 0038 实际上管道与送电线路不可能理想平行的， 在有些区域发生斜接近， 有些区域甚 至交叉， 如图 3。对于斜接近段 P1P2 的管道 ( 如附图。
18、 4)， 不同位置的感应纵电动势 E 与送 电线路和管道间的互感有关， 因此不同位置的 E 不同， 即此时求解上述传输线 方程是非常复杂的。 0039 图中 a1、 a2 为 P1、 P2 点与送电线路的最短间距。简化计算时， 可以使用平均 间距代替进行互感系数的求解。此时将斜接近段等效为平行段， 根据 E -jMIphase， 得到 E 近似为常数， 进而因此在斜接近段 P1-P2 管道上的电压、 电流 分布有如下形式解， 如下 0040 Up1-2 -ZA1exp(x)-B1exp(-x) (6) 0041 0042 Up1-2P1P2段电压分布， x 为从 P1到待求点的距离。 0043。
19、 Ip1-2P1P2段电流分布。A1、 B1P1P2段待定系数。 0044 将附图 3 的连续曲折的管道细分， 利用等效方法将斜接近段、 交叉段转化为平行 段， 结合土壤电阻率及平行间距查找诺模互感系数图、 采用简化计算公式得到互感系数， 进 而求得曲折段的近似等效单位感应纵电动势， 接着利用单回传输线方程的近似解， 得到每 个细分段的电压、 电流表达式。对于同一端点而言， 其即为前一细分段的末端点， 也为后一 细分段的起始点， 因此再根据同一点电压电流值唯一得到一个规则关于待定系数 Ai、 Bi的 多元一次方程组， 联合管道首末端边界条件， 求解该多元一次方程组， 可以求出各细分段管 道对地。
20、电压、 电流分布， 从而获得送电线路与管道非理想平行时感性耦合干扰。 0045 实施例 ： 0046 假设对如下情况进行分析 ： 管道与送电线路相对位置关系如附图 5， 管道为三 层 PE 防腐层， 单回送电线路电流 I 1000A， 送电线路高 20m， 管道埋深 2m， 土壤电阻率 50m。管道参数如表 1。 说 明 书 CN 102435877 A CN 102435883 A4/4 页 7 0047 表 1 管道参数 0048 0049 运用本发明的技术方法， 计算获得的管道不同点的干扰电压与国外知名软件 CDEGS 计算结果对比如下表 2， 对比显示两者误差较小， 本发明有较好的准确性。 0050 表 2 本发明计算方法与国外知名软件计算结果对比 0051 0052 0053 本说明书中未作详细描述的内容属于本领域专业技术人员公知的现有技术。 说 明 书 CN 102435877 A CN 102435883 A1/2 页 8 图 1 图 2 图 3 说 明 书 附 图 CN 102435877 A CN 102435883 A2/2 页 9 图 4 图 5 说 明 书 附 图 CN 102435877 A 。