CN113701689B - 一种耐张塔输电线路等值覆冰厚度测量方法及系统 - Google Patents

一种耐张塔输电线路等值覆冰厚度测量方法及系统 Download PDF

Info

Publication number
CN113701689B
CN113701689B CN202110847432.6A CN202110847432A CN113701689B CN 113701689 B CN113701689 B CN 113701689B CN 202110847432 A CN202110847432 A CN 202110847432A CN 113701689 B CN113701689 B CN 113701689B
Authority
CN
China
Prior art keywords
tension
transmission line
string
equivalent
power transmission
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
Application number
CN202110847432.6A
Other languages
English (en)
Other versions
CN113701689A (zh
Inventor
马晓红
牛唯
毛先胤
黄欢
杨旗
文屹
黄增浩
李�昊
何锦强
赵林杰
廖永力
张志强
朱登杰
王斌
魏发生
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China South Power Grid International Co ltd
Guizhou Power Grid Co Ltd
Original Assignee
China South Power Grid International Co ltd
Guizhou Power Grid Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China South Power Grid International Co ltd, Guizhou Power Grid Co Ltd filed Critical China South Power Grid International Co ltd
Priority to CN202110847432.6A priority Critical patent/CN113701689B/zh
Publication of CN113701689A publication Critical patent/CN113701689A/zh
Application granted granted Critical
Publication of CN113701689B publication Critical patent/CN113701689B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B21/00—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant
    • G01B21/02—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant for measuring length, width, or thickness
    • G01B21/08—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant for measuring length, width, or thickness for measuring thickness
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/15—Correlation function computation including computation of convolution operations

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Computational Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Algebra (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Computing Systems (AREA)
  • Insulators (AREA)
  • Length Measuring Devices With Unspecified Measuring Means (AREA)

Abstract

本发明公开了一种耐张塔输电线路等值覆冰厚度测量方法及系统，包括建立输电线路的有限元仿真模型，获得无覆冰工况下耐张段首尾端的第一耐张串轴向张力和第二耐张串轴向张力；向有限元仿真模型施加覆冰荷载，获得对应的第三耐张串轴向张力和第四耐张串轴向张力；当第一拉力监测值大于第一耐张串轴向张力，且第二拉力监测值大于第二耐张串轴向张力时，计算输电线路的第一等值覆冰厚度和第二等值覆冰厚度；基于第一等值覆冰厚度、第二等值覆冰厚度与覆冰厚度理论值的比较结果确定拉力传感器的布置位置，采用称重法计算最终等值覆冰厚度。本发明能够减小覆冰厚度计算结果的误差，从而提高输电线路覆冰监测准确度。

Description

一种耐张塔输电线路等值覆冰厚度测量方法及系统
技术领域
本发明涉及输电线路覆冰监测领域，尤其涉及一种耐张塔输电线路等值覆冰厚度测量方法及系统。
背景技术
输电线路覆冰严重影响电网的安全稳定运行，极端天气情况的出现使得输电线大面积覆冰事件频发，从而发生输电线路断线、导线舞动、倒塔、绝缘子闪络等事故，给电网部门造成巨大的经济财产损失。由于覆冰带来巨大危害，开展输电线路覆冰监测越来越迫切。称重法因其原理简单、技术成熟，在国内电网覆冰监测获得广泛应用，并逐步成为行业规范。而拉力传感器是称重法应用的关键器件，其安装位置的不规范导致覆冰计算结果与实际值的误差较大，从而影响覆冰监测结果的准确度。
发明内容
本发明提供一种耐张塔输电线路等值覆冰厚度测量方法及系统，以解决现有技术覆冰计算结果与实际值的误差较大的问题，本发明能够减小覆冰厚度计算结果的误差，从而提高输电线路覆冰监测准确度，避免覆冰监测系统预警误报。
本发明实施例提供了一种耐张塔输电线路等值覆冰厚度测量方法，包括：
建立输电线路的有限元仿真模型，获得无覆冰工况下耐张段的首端的第一耐张串轴向张力和所述耐张段的尾端的第二耐张串轴向张力；
向所述有限元仿真模型施加预设覆冰厚度理论值的覆冰荷载，获得与所述预设覆冰厚度理论值对应的所述首端的第三耐张串轴向张力和所述尾端的第四耐张串轴向张力；
当安装在所述首端的耐张塔上的第一拉力传感器监测到的第一拉力监测值大于所述第一耐张串轴向张力，且安装在所述尾端的耐张塔上的第二拉力传感器监测到的第二拉力监测值大于所述第二耐张串轴向张力时，根据所述第一耐张串轴向张力和所述第三耐张串轴向张力计算得到输电线路的第一等值覆冰厚度，根据所述第二耐张串轴向张力和所述第四耐张串轴向张力计算得到所述输电线路的第二等值覆冰厚度；
基于所述第一等值覆冰厚度、所述第二等值覆冰厚度与所述覆冰厚度理论值的比较结果确定拉力传感器的布置位置，采用称重法计算得到所述输电线路的最终等值覆冰厚度。
进一步地，在所述基于所述第一等值覆冰厚度、所述第二等值覆冰厚度与所述覆冰厚度理论值的比较结果确定拉力传感器的布置位置后，所述方法还包括：
根据所述布置位置，将第三拉力传感器布置在所述首端的耐张塔上或所述尾端的耐张塔上；
根据所述有限元仿真模型，获得无覆冰工况下的第五耐张串轴向张力；
向所述有限元仿真模型施加若干均匀覆冰或不均匀覆冰的覆冰荷载，获得相应的若干第二耐张串轴向张力，绘制覆冰厚度-耐张串轴向张力的散点图；
根据所述散点图，采用数据拟合方法，获取关于覆冰厚度与耐张串轴向张力的曲线拟合公式；
当第三拉力传感器监测到的第三拉力监测值大于所述第五耐张串轴向张力时，将所述第三拉力监测值代入所述曲线拟合公式，得到所述输电线路的最终等值覆冰厚度。
进一步地，在所述有限元仿真模型中，通过沿线向节点施加集中力等效模拟覆冰荷载，其中，根据下式计算得到所述集中力：
ω＝0.9πgb(b+D)×10-3
T＝ωL＝[0.9πgb(b+D)×10-3]L
其中，D为输电线外径，b为覆冰厚度，g为重力加速度，ω为输电线单位长度重力，L为有限元单元长度，T为施加到有限元仿真模型的集中力大小。
进一步地，所述根据所述第一耐张串轴向张力和所述第三耐张串轴向张力计算得到输电线路的第一等值覆冰厚度，包括：
根据所述第一耐张串轴向张力计算得到输电线路的第一平均等效长度：
式中，F01为所述第一耐张串轴向张力，n0为输电线分裂数，ω0为输电线单位长度重量，θ为绝缘子串夹角；
根据所述第一耐张串轴向张力计算得到所述输电线路的第一最大等效长度：
式中，F01为所述第一耐张串轴向张力，n0为输电线分裂数，ω0为输电线单位长度重量，θ为绝缘子串夹角；
根据所述第一平均等效长度、所述第一最大等效长度和所述第三耐张串轴向张力计算得到所述输电线路的第一等值覆冰厚度：
式中，F1为所述第三耐张串轴向张力，ω0为输电线单位长度重量，n0为输电线分裂数，D为输电线外径，lav1为所述第一平均等效长度，lmax1为所述第一最大等效长度。
进一步地，所述根据所述第二耐张串轴向张力和所述第四耐张串轴向张力计算得到所述输电线路的第二等值覆冰厚度，包括：
根据所述第二耐张串轴向张力计算得到所述输电线路的第二平均等效长度：
式中，F02为所述第二耐张串轴向张力，n0为输电线分裂数，ω0为输电线单位长度重量，θ为绝缘子串夹角；
根据所述第二耐张串轴向张力计算得到所述输电线路的第二最大等效长度：
式中，F02为所述第二耐张串轴向张力，n0为输电线分裂数，ω0为输电线单位长度重量，θ为绝缘子串夹角；
根据所述第二平均等效长度、所述第二最大等效长度和所述第四耐张串轴向张力计算得到所述输电线路的第二等值覆冰厚度：
式中，F2为所述第四耐张串轴向张力，ω0为输电线单位长度重量，n0为输电线分裂数，D为输电线外径，lav2为所述第二平均等效长度，lmax2为所述第二最大等效长度。
进一步地，所述绘制覆冰厚度-耐张串轴向张力的散点图，包括：
当所述覆冰荷载为不均匀覆冰时，所述散点图的覆冰厚度为所述耐张段的平均覆冰厚度；
其中，根据预设的第一公式计算所述平均覆冰厚度，或，根据预设的第二公式计算所述平均覆冰厚度；
所述第一公式为：
在所述第一公式中，bi为第i档的覆冰厚度，Si为第i档的输电线长度，为耐张段的平均覆冰厚度；
所述第二公式为：
在所述第二公式中，bi为第i档的覆冰厚度，li为第i档档距长度，为耐张段的平均覆冰厚度。
相应地，本发明实施例还提供一种耐张塔输电线路等值覆冰厚度测量系统，包括：
建立仿真模型模块，用于建立输电线路的有限元仿真模型，获得无覆冰工况下耐张段的首端的第一耐张串轴向张力和所述耐张段的尾端的第二耐张串轴向张力；
施加覆冰荷载模块，用于向所述有限元仿真模型施加预设覆冰厚度理论值的覆冰荷载，获得与所述预设覆冰厚度理论值对应的所述首端的第三耐张串轴向张力和所述尾端的第四耐张串轴向张力；
第一计算模块，用于当安装在所述首端的耐张塔上的第一拉力传感器监测到的第一拉力监测值大于所述第一耐张串轴向张力，且安装在所述尾端的耐张塔上的第二拉力传感器监测到的第二拉力监测值大于所述第二耐张串轴向张力时，根据所述第一耐张串轴向张力和所述第三耐张串轴向张力计算得到输电线路的第一等值覆冰厚度，根据所述第二耐张串轴向张力和所述第四耐张串轴向张力计算得到所述输电线路的第二等值覆冰厚度；
第二计算模块，用于基于所述第一等值覆冰厚度、所述第二等值覆冰厚度与所述覆冰厚度理论值的比较结果确定拉力传感器的布置位置，采用称重法计算得到所述输电线路的最终等值覆冰厚度。
与现有技术相比，本发明实施例所提供的耐张塔输电线路等值覆冰厚度测量方法及系统，通过建立有限元仿真模型，获得无覆冰工况下耐张段首尾两端的耐张串轴向张力，向有限元仿真模型施加预设覆冰厚度理论值的覆冰荷载，获得对应的耐张段首尾两端的耐张串轴向张力，根据上述耐张串轴向张力计算得到第一等值覆冰厚度和第二等值覆冰厚度，基于所述第一等值覆冰厚度、所述第二等值覆冰厚度与所述覆冰厚度理论值的比较结果确定拉力传感器的布置位置，采用称重法计算得到所述输电线路的最终等值覆冰厚度。本发明实施例通过耐张串轴向张力计算得到等值覆冰厚度，并与预设覆冰厚度理论值对比，确定拉力传感器的布置位置，减小覆冰厚度计算结果的误差，从而提高输电线路覆冰监测准确度，避免覆冰监测系统预警误报。
附图说明
图1是本发明实施例提供的一种耐张塔输电线路等值覆冰厚度测量方法的流程示意图；
图2是本发明实施例提供的一种耐张塔输电线路等值覆冰厚度测量系统的结构示意图；
图3是本发明实施例提供的一种耐张塔所在的耐张段为“三塔两档”的输电线路模型图；
图4是本发明实施例提供的覆冰厚度-耐张串轴向张力散点图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
参见图1，是本发明实施例提供的一种耐张塔输电线路等值覆冰厚度测量方法的流程示意图，所述测量方法包括：
S11、建立输电线路的有限元仿真模型，获得无覆冰工况下耐张段的首端的第一耐张串轴向张力和所述耐张段的尾端的第二耐张串轴向张力；
具体地，通过以下步骤建立输电线路的有限元仿真模型：
结合输电线路的结构特点，选择适当的有限元单元类型，并确定其实常数和材料属性；基于直接迭代法确定输电线的初始形态，求出满足预应力状态和边界条件的结构形状；对绝缘子串建模，建立输电线-绝缘子串耦合模型。其中，悬垂串上端位置固定不动，下端与输电线铰接，悬垂串与输电线连接点能够移动；根据仿真工况设定覆冰荷载；指定约束条件，定义分析类型、荷载数据和荷载步选项，打开大变形开关，进行有限元求解；有限元求解结果后处理部分，定义轴向张力单元表，列表显示单元轴向张力，根据单元编号提取耐张串轴向张力。可以理解的，对有限元仿真模型不施加静力，通过有限元求解，获得无覆冰工况下的耐张串轴向张力。
S12、向所述有限元仿真模型施加预设覆冰厚度理论值的覆冰荷载，获得与所述预设覆冰厚度理论值对应的所述首端的第三耐张串轴向张力和所述尾端的第四耐张串轴向张力；
可以理解的，将预设覆冰厚度理论值的覆冰荷载等效换算为集中力，并将该集中力施加到有限元仿真模型，通过有限元求解，获得与所述预设覆冰厚度理论值对应的所述首端的第三耐张串轴向张力和所述尾端的第四耐张串轴向张力。
S13、当安装在所述首端的耐张塔上的第一拉力传感器监测到的第一拉力监测值大于所述第一耐张串轴向张力，且安装在所述尾端的耐张塔上的第二拉力传感器监测到的第二拉力监测值大于所述第二耐张串轴向张力时，根据所述第一耐张串轴向张力和所述第三耐张串轴向张力计算得到输电线路的第一等值覆冰厚度，根据所述第二耐张串轴向张力和所述第四耐张串轴向张力计算得到所述输电线路的第二等值覆冰厚度；
可以理解的，当所述第一拉力监测值不大于所述第一耐张串轴向张力，或第二拉力监测值不大于所述第二耐张串轴向张力时，输电线路的等值覆冰厚度为0；
S14、基于所述第一等值覆冰厚度、所述第二等值覆冰厚度与所述覆冰厚度理论值的比较结果确定拉力传感器的布置位置，采用称重法计算得到所述输电线路的最终等值覆冰厚度。
具体地，基于所述第一等值覆冰厚度b1、所述第二等值覆冰厚度b2与所述覆冰厚度理论值b0的比较结果确定拉力传感器的布置位置：
当|b1-b0|＜|b2-b0|时，拉力传感器布置在耐张段首端的耐张塔上；
当|b1-b0|＞|b2-b0|时，拉力传感器布置在耐张段尾端的耐张塔上；
当|b1-b0|＝|b2-b0|时，拉力传感器布置在耐张段首端的耐张塔上或耐张段尾端的耐张塔上。
在另一可选实施例中，在所述基于所述第一等值覆冰厚度、所述第二等值覆冰厚度与所述覆冰厚度理论值的比较结果确定拉力传感器的布置位置后，所述方法还包括：
根据所述布置位置，将第三拉力传感器布置在所述首端的耐张塔上或所述尾端的耐张塔上；
根据所述有限元仿真模型，获得无覆冰工况下的第五耐张串轴向张力；
向所述有限元仿真模型施加若干均匀覆冰或不均匀覆冰的覆冰荷载，获得相应的若干第二耐张串轴向张力，绘制覆冰厚度-耐张串轴向张力的散点图；
根据所述散点图，采用数据拟合方法，获取关于覆冰厚度与耐张串轴向张力的曲线拟合公式；
具体地，分析数据点在平面直角坐标系的分布特征，把握覆冰厚度随耐张串轴向张力的变化趋势，计算覆冰厚度与耐张串轴向张力的相关系数，以此表征上述两个变量之间相关关系密切程度，相关系数的表达式为
其中，xi为第i个耐张串轴向张力数据，yi为第i个覆冰厚度数据，为耐张串轴向张力数据平均值，/>为覆冰厚度数据平均值，n为数据点总数。
根据相关系数，结合散点图中离散点的分布，选择合适的函数类型对数据点进行拟合，所选择的函数类型尽可能是初等函数，保证曲线拟合公式的简洁性；
当第三拉力传感器监测到的第三拉力监测值大于所述第五耐张串轴向张力时，将所述第三拉力监测值代入所述曲线拟合公式，得到所述输电线路的最终等值覆冰厚度。
具体地，当所述第三拉力监测值不大于所述第五耐张串轴向张力时，输电线路的等值覆冰厚度为0；当所述第三拉力监测值大于所述第五耐张串轴向张力时，将所述第三拉力监测值代入所述曲线拟合公式，得到所述输电线路的最终等值覆冰厚度
在本发明实施例中，通过有限元仿真得到与覆冰厚度对应的耐张串轴向张力，推导得到覆冰厚度与耐张串轴向张力的曲线拟合公式，实现耐张塔输电线路等值覆冰厚度，提高了输电线路覆冰计算方法的简洁性。
作为上述方案的改进，在所述有限元仿真模型中，通过沿线向节点施加集中力等效模拟覆冰荷载，其中，根据下式计算得到所述集中力：
ω＝0.9πgb(b+D)×10-3
T＝ωL＝[0.9πgb(b+D)×10-3]L
其中，D为输电线外径，b为覆冰厚度，g为重力加速度，ω为输电线单位长度重力，L为有限元单元长度，T为施加到有限元仿真模型的集中力大小。
作为上述方案的改进，所述根据所述第一耐张串轴向张力和所述第三耐张串轴向张力计算得到输电线路的第一等值覆冰厚度，包括：
根据所述第一耐张串轴向张力计算得到输电线路的第一平均等效长度：
式中，F01为所述第一耐张串轴向张力，n0为输电线分裂数，ω0为输电线单位长度重量，θ为绝缘子串夹角；
根据所述第一耐张串轴向张力计算得到所述输电线路的第一最大等效长度：
式中，F01为所述第一耐张串轴向张力，n0为输电线分裂数，ω0为输电线单位长度重量，θ为绝缘子串夹角；
根据所述第一平均等效长度、所述第一最大等效长度和所述第三耐张串轴向张力计算得到所述输电线路的第一等值覆冰厚度：
式中，F1为所述第三耐张串轴向张力，ω0为输电线单位长度重量，n0为输电线分裂数，D为输电线外径，lav1为所述第一平均等效长度，lmax1为所述第一最大等效长度。
作为上述方案的改进，所述根据所述第二耐张串轴向张力和所述第四耐张串轴向张力计算得到所述输电线路的第二等值覆冰厚度，包括：
根据所述第二耐张串轴向张力计算得到所述输电线路的第二平均等效长度：
式中，F02为所述第二耐张串轴向张力，n0为输电线分裂数，ω0为输电线单位长度重量，θ为绝缘子串夹角；
根据所述第二耐张串轴向张力计算得到所述输电线路的第二最大等效长度：
式中，F02为所述第二耐张串轴向张力，n0为输电线分裂数，ω0为输电线单位长度重量，θ为绝缘子串夹角；
根据所述第二平均等效长度、所述第二最大等效长度和所述第四耐张串轴向张力计算得到所述输电线路的第二等值覆冰厚度：
式中，F2为所述第四耐张串轴向张力，ω0为输电线单位长度重量，n0为输电线分裂数，D为输电线外径，lav2为所述第二平均等效长度，lmax2为所述第二最大等效长度。
作为上述方案的改进，所述绘制覆冰厚度-耐张串轴向张力的散点图，包括：
当所述覆冰荷载为不均匀覆冰时，所述散点图的覆冰厚度为所述耐张段的平均覆冰厚度；
其中，根据预设的第一公式计算所述平均覆冰厚度，或，根据预设的第二公式计算所述平均覆冰厚度；
所述第一公式为：
在所述第一公式中，bi为第i档的覆冰厚度，Si为第i档的输电线长度，为耐张段的平均覆冰厚度；
所述第二公式为：
在所述第二公式中，bi为第i档的覆冰厚度，li为第i档档距长度，为耐张段的平均覆冰厚度。
可以理解的，为简化起见，档距内输电线长度可用档距长度代替，因此可以根据预设的第二公式计算所述平均覆冰厚度。
需要说明的，在本发明实施例中，对于不均匀覆冰，通过有限元仿真的方法，获得不均匀覆冰下耐张串轴向张力，用耐张段覆冰量的平均覆冰厚度表征不均匀覆冰时耐张段覆冰程度，建立平均覆冰厚度与耐张串轴向张力的数据集；根据预设的第一公式或预设的第二公式计算得到平均覆冰厚度；对于均匀覆冰，通过改变覆冰厚度大小，获得不同覆冰厚度对应的耐张串张力大小，其覆冰厚度就是平均覆冰厚度；
以平均覆冰厚度为因变量，耐张串轴向张力为自变量，在平面直角坐标系中绘制平均覆冰厚度与耐张串轴向张力两个变量的散点图。
在一具体实施中，参见图3，上述方法作用于某110kV输电线路耐张该该该耐张段的设计冰厚为30mm，导线分裂数为1，导线型号为JLHA1/G1A-185/30-26/7，导线单位长度重量为0.73kg/m，导线外径为18.88mm，绝缘子串重量为150kg，绝缘子串型为I串，绝缘子串夹角为0，某时刻拉力传感器监测到的拉力监测值为1957kg。设档距长度分别为l1＝150m，l2＝196m，大号侧悬挂点与小号侧悬挂点的高度差分别为h1＝-40m，h2＝50m。
运用有限元法和全尺寸仿真模拟方法，建立输电线路的有限元仿真模型：根据输电线和绝缘子串的结构特点和连接方式，选择合适的单位类型进行仿真模拟，建立输电线-绝缘子串耦合模型，采用直接迭代法确定输电线初始形态以及运用集中力模拟方法等效模拟覆冰荷载。指定约束条件，定义分析类型、荷载数据和荷载步选项，然后进行有限元求解，提取耐张串轴向张力。
通过有限元仿真计算，获得无风且无覆冰工况下的耐张串轴向张力为263kg。
分别施加均匀覆冰和不均匀覆冰的覆冰荷载，通过有限元求解，获得对应的耐张串轴向张力，如下表所示：
由上表可知，对于均匀覆冰的情况，通过向大号侧档距与小号侧档距施加相同覆冰厚度的覆冰荷载模拟均匀覆冰情况，在大号侧档距和小号侧档距同时施加5mm、10mm、15mm、20mm、25mm或30mm五种不同覆冰厚度的覆冰荷载进行多次仿真计算，获得与其对应的耐张串轴向张力；对于不均匀覆冰的情况，通过向大号侧档距与小号侧杆塔档距不同覆冰厚度的覆冰荷载模拟不均匀覆冰情况。
根据有限元计算结果，绘制覆冰厚度-耐张串轴向张力散点图，如图4所示，其中，覆冰厚度和耐张串轴向张力呈正相关，且覆冰厚度和耐张串张力呈线性关系，计算相关系数R2＝0.9313。根据相关系数，采用线性拟合方式进行拟合，曲线拟合公式为b＝0.0185F-0.2089。
由于拉力监测值1957kg大于无风且无覆冰工况下的耐张串轴向张力263kg，故将拉力监测值1957kg代入曲线拟合公式b＝0.0185F-0.2089，求得输电线路的最终等值覆冰厚度b为36mm。
参见图2，是本发明实施例提供的一种耐张塔输电线路等值覆冰厚度测量系统结构示意图，包括：
建立仿真模型模块21，用于建立输电线路的有限元仿真模型，获得无覆冰工况下耐张段的首端的第一耐张串轴向张力和所述耐张段的尾端的第二耐张串轴向张力；
施加覆冰荷载模块22，用于向所述有限元仿真模型施加预设覆冰厚度理论值的覆冰荷载，获得与所述预设覆冰厚度理论值对应的所述首端的第三耐张串轴向张力和所述尾端的第四耐张串轴向张力；
第一计算模块23，用于当安装在所述首端的耐张塔上的第一拉力传感器监测到的第一拉力监测值大于所述第一耐张串轴向张力，且安装在所述尾端的耐张塔上的第二拉力传感器监测到的第二拉力监测值大于所述第二耐张串轴向张力时，根据所述第一耐张串轴向张力和所述第三耐张串轴向张力计算得到输电线路的第一等值覆冰厚度，根据所述第二耐张串轴向张力和所述第四耐张串轴向张力计算得到所述输电线路的第二等值覆冰厚度；
第二计算模块24，用于基于所述第一等值覆冰厚度、所述第二等值覆冰厚度与所述覆冰厚度理论值的比较结果确定拉力传感器的布置位置，采用称重法计算得到所述输电线路的最终等值覆冰厚度。
优选地，所述系统还包括：
布置模块，用于根据所述布置位置，将第三拉力传感器布置在所述首端的耐张塔上或所述尾端的耐张塔上；
第三计算模块，用于根据所述有限元仿真模型，获得无覆冰工况下的第五耐张串轴向张力；
绘制散点图模块，用于向所述有限元仿真模型施加若干均匀覆冰或不均匀覆冰的覆冰荷载，获得相应的若干第二耐张串轴向张力，绘制覆冰厚度-耐张串轴向张力的散点图；
曲线拟合模块，用于根据所述散点图，采用数据拟合方法，获取关于覆冰厚度与耐张串轴向张力的曲线拟合公式；
第四计算模块，用于当第三拉力传感器监测到的第三拉力监测值大于所述第五耐张串轴向张力时，将所述第三拉力监测值代入所述曲线拟合公式，得到所述输电线路的最终等值覆冰厚度。
优选地，在所述有限元仿真模型中，通过沿线向节点施加集中力等效模拟覆冰荷载，其中，根据下式计算得到所述集中力：
ω＝0.9πgb(b+D)×10-3
T＝ωL＝[0.9πgb(b+D)×10-3]L
其中，D为输电线外径，b为覆冰厚度，g为重力加速度，ω为输电线单位长度重力，L为有限元单元长度，T为施加到有限元仿真模型的集中力大小。
优选地，所述第一计算模块还用于：
根据所述第一耐张串轴向张力计算得到输电线路的第一平均等效长度：
式中，F01为所述第一耐张串轴向张力，n0为输电线分裂数，ω0为输电线单位长度重量，θ为绝缘子串夹角；
根据所述第一耐张串轴向张力计算得到所述输电线路的第一最大等效长度：
式中，F01为所述第一耐张串轴向张力，n0为输电线分裂数，ω0为输电线单位长度重量，θ为绝缘子串夹角；
根据所述第一平均等效长度、所述第一最大等效长度和所述第三耐张串轴向张力计算得到所述输电线路的第一等值覆冰厚度：
式中，F1为所述第三耐张串轴向张力，ω0为输电线单位长度重量，n0为输电线分裂数，D为输电线外径，lav1为所述第一平均等效长度，lmax1为所述第一最大等效长度。
优选地，所述第一计算模块还用于：
根据所述第二耐张串轴向张力计算得到所述输电线路的第二平均等效长度：
式中，F02为所述第二耐张串轴向张力，n0为输电线分裂数，ω0为输电线单位长度重量，θ为绝缘子串夹角；
根据所述第二耐张串轴向张力计算得到所述输电线路的第二最大等效长度：
式中，F02为所述第二耐张串轴向张力，n0为输电线分裂数，ω0为输电线单位长度重量，θ为绝缘子串夹角；
根据所述第二平均等效长度、所述第二最大等效长度和所述第四耐张串轴向张力计算得到所述输电线路的第二等值覆冰厚度：
式中，F2为所述第四耐张串轴向张力，ω0为输电线单位长度重量，n0为输电线分裂数，D为输电线外径，lav2为所述第二平均等效长度，lmax2为所述第二最大等效长度。
优选地，所述绘制散点图模块还用于：
当所述覆冰荷载为不均匀覆冰时，所述散点图的覆冰厚度为所述耐张段的平均覆冰厚度；
其中，根据预设的第一公式计算所述平均覆冰厚度，或，根据预设的第二公式计算所述平均覆冰厚度；
所述第一公式为：
在所述第一公式中，bi为第i档的覆冰厚度，Si为第i档的输电线长度，为耐张段的平均覆冰厚度；
所述第二公式为：
在所述第二公式中，bi为第i档的覆冰厚度，li为第i档档距长度，为耐张段的平均覆冰厚度。
需要说明的是，本发明实施例所提供的耐张塔输电线路等值覆冰厚度测量系统，能够实现上述任一实施例所述的耐张塔输电线路等值覆冰厚度测量方法的所有流程，系统中的各个单元、子单元的作用以及实现的技术效果分别与上述实施例所述的耐张塔输电线路等值覆冰厚度测量方法的作用以及实现的技术效果对应相同，这里不再赘述。
本发明实施例提供的一种耐张塔输电线路等值覆冰厚度测量方法及系统，通过建立有限元仿真模型，获得无覆冰工况下耐张段首尾两端的耐张串轴向张力，向有限元仿真模型施加预设覆冰厚度理论值的覆冰荷载，获得对应的耐张段首尾两端的耐张串轴向张力，根据上述耐张串轴向张力计算得到第一等值覆冰厚度和第二等值覆冰厚度，基于所述第一等值覆冰厚度、所述第二等值覆冰厚度与所述覆冰厚度理论值的比较结果确定拉力传感器的布置位置，采用称重法计算得到所述输电线路的最终等值覆冰厚度。本发明实施例通过耐张串轴向张力计算得到等值覆冰厚度，并与预设覆冰厚度理论值对比，确定拉力传感器的布置位置，减小覆冰厚度计算结果的误差，从而提高输电线路覆冰监测准确度，为输电线路覆冰监测预警提供重要技术支持。
以上所述是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也视为本发明的保护范围。

Claims (4)

1.一种耐张塔输电线路等值覆冰厚度测量方法，其特征在于，包括：
建立输电线路的有限元仿真模型，获得无覆冰工况下耐张段的首端的第一耐张串轴向张力和所述耐张段的尾端的第二耐张串轴向张力；
向所述有限元仿真模型施加预设覆冰厚度理论值的覆冰荷载，获得与所述预设覆冰厚度理论值对应的所述首端的第三耐张串轴向张力和所述尾端的第四耐张串轴向张力；
当安装在所述首端的耐张塔上的第一拉力传感器监测到的第一拉力监测值大于所述第一耐张串轴向张力，且安装在所述尾端的耐张塔上的第二拉力传感器监测到的第二拉力监测值大于所述第二耐张串轴向张力时，根据所述第一耐张串轴向张力和所述第三耐张串轴向张力计算得到输电线路的第一等值覆冰厚度，根据所述第二耐张串轴向张力和所述第四耐张串轴向张力计算得到所述输电线路的第二等值覆冰厚度；
基于所述第一等值覆冰厚度、所述第二等值覆冰厚度与所述覆冰厚度理论值的比较结果确定拉力传感器的布置位置，采用称重法计算得到所述输电线路的最终等值覆冰厚度；
在所述基于所述第一等值覆冰厚度、所述第二等值覆冰厚度与所述覆冰厚度理论值的比较结果确定拉力传感器的布置位置后，所述方法还包括：
根据所述布置位置，将第三拉力传感器布置在所述首端的耐张塔上或所述尾端的耐张塔上；
根据所述有限元仿真模型，获得无覆冰工况下的第五耐张串轴向张力；
向所述有限元仿真模型施加若干均匀覆冰或不均匀覆冰的覆冰荷载，获得相应的若干第二耐张串轴向张力，绘制覆冰厚度-耐张串轴向张力的散点图；
根据所述散点图，采用数据拟合方法，获取关于覆冰厚度与耐张串轴向张力的曲线拟合公式；
当第三拉力传感器监测到的第三拉力监测值大于所述第五耐张串轴向张力时，将所述第三拉力监测值代入所述曲线拟合公式，得到所述输电线路的最终等值覆冰厚度；
其中，所述根据所述第一耐张串轴向张力和所述第三耐张串轴向张力计算得到输电线路的第一等值覆冰厚度，包括：
根据所述第一耐张串轴向张力计算得到输电线路的第一平均等效长度：
式中，F01为所述第一耐张串轴向张力，n0为输电线分裂数，ω0为输电线单位长度重量，θ为绝缘子串夹角；
根据所述第一耐张串轴向张力计算得到所述输电线路的第一最大等效长度：
式中，F01为所述第一耐张串轴向张力，n0为输电线分裂数，ω0为输电线单位长度重量，θ为绝缘子串夹角；
根据所述第一平均等效长度、所述第一最大等效长度和所述第三耐张串轴向张力计算得到所述输电线路的第一等值覆冰厚度：
式中，F1为所述第三耐张串轴向张力，ω0为输电线单位长度重量，n0为输电线分裂数，D为输电线外径，lav1为所述第一平均等效长度，lmax1为所述第一最大等效长度；
所述根据所述第二耐张串轴向张力和所述第四耐张串轴向张力计算得到所述输电线路的第二等值覆冰厚度，包括：
根据所述第二耐张串轴向张力计算得到所述输电线路的第二平均等效长度：
式中，F02为所述第二耐张串轴向张力，n0为输电线分裂数，ω0为输电线单位长度重量，θ为绝缘子串夹角；
根据所述第二耐张串轴向张力计算得到所述输电线路的第二最大等效长度：
式中，F02为所述第二耐张串轴向张力，n0为输电线分裂数，ω0为输电线单位长度重量，θ为绝缘子串夹角；
根据所述第二平均等效长度、所述第二最大等效长度和所述第四耐张串轴向张力计算得到所述输电线路的第二等值覆冰厚度：
式中，F2为所述第四耐张串轴向张力，ω0为输电线单位长度重量，n0为输电线分裂数，D为输电线外径，lav2为所述第二平均等效长度，lmax2为所述第二最大等效长度；
所述基于所述第一等值覆冰厚度、所述第二等值覆冰厚度与所述覆冰厚度理论值的比较结果确定拉力传感器的布置位置包括：
当|b1-b0|＜|b2-b0|时，拉力传感器布置在耐张段首端的耐张塔上；
当|b1-b0|＞|b2-b0|时，拉力传感器布置在耐张段尾端的耐张塔上；
当|b1-b0|＝|b2-b0|时，拉力传感器布置在耐张段首端的耐张塔上或耐张段尾端的耐张塔上；
其中，b1为第一等值覆冰厚度，b2为第二等值覆冰厚度，b0为覆冰厚度理论值。
2.如权利要求1所述的耐张塔输电线路等值覆冰厚度测量方法，其特征在于，在所述有限元仿真模型中，通过沿线向节点施加集中力等效模拟覆冰荷载，其中，根据下式计算得到所述集中力：
ω＝0.9πgb(b+D)×10-3
T＝ωL＝[0.9πgb(b+D)×10-3]L
其中，D为输电线外径，b为覆冰厚度，g为重力加速度，ω为输电线单位长度重力，L为有限元单元长度，T为施加到有限元仿真模型的集中力大小。
3.如权利要求1所述的耐张塔输电线路等值覆冰厚度测量方法，其特征在于，所述绘制覆冰厚度-耐张串轴向张力的散点图，包括：
当所述覆冰荷载为不均匀覆冰时，所述散点图的覆冰厚度为所述耐张段的平均覆冰厚度；
其中，根据预设的第一公式计算所述平均覆冰厚度，或，根据预设的第二公式计算所述平均覆冰厚度；
所述第一公式为：
在所述第一公式中，bi为第i档的覆冰厚度，Si为第i档的输电线长度，为耐张段的平均覆冰厚度；
所述第二公式为：
在所述第二公式中，bi为第i档的覆冰厚度，li为第i档档距长度，为耐张段的平均覆冰厚度。
4.一种耐张塔输电线路等值覆冰厚度测量系统，其特征在于，包括：
建立仿真模型模块，用于建立输电线路的有限元仿真模型，获得无覆冰工况下耐张段的首端的第一耐张串轴向张力和所述耐张段的尾端的第二耐张串轴向张力；
施加覆冰荷载模块，用于向所述有限元仿真模型施加预设覆冰厚度理论值的覆冰荷载，获得与所述预设覆冰厚度理论值对应的所述首端的第三耐张串轴向张力和所述尾端的第四耐张串轴向张力；
第一计算模块，用于当安装在所述首端的耐张塔上的第一拉力传感器监测到的第一拉力监测值大于所述第一耐张串轴向张力，且安装在所述尾端的耐张塔上的第二拉力传感器监测到的第二拉力监测值大于所述第二耐张串轴向张力时，根据所述第一耐张串轴向张力和所述第三耐张串轴向张力计算得到输电线路的第一等值覆冰厚度，根据所述第二耐张串轴向张力和所述第四耐张串轴向张力计算得到所述输电线路的第二等值覆冰厚度；
第二计算模块，用于基于所述第一等值覆冰厚度、所述第二等值覆冰厚度与所述覆冰厚度理论值的比较结果确定拉力传感器的布置位置，采用称重法计算得到所述输电线路的最终等值覆冰厚度；
所述系统还包括：
布置模块，用于根据所述布置位置，将第三拉力传感器布置在所述首端的耐张塔上或所述尾端的耐张塔上；
第三计算模块，用于根据所述有限元仿真模型，获得无覆冰工况下的第五耐张串轴向张力；
绘制散点图模块，用于向所述有限元仿真模型施加若干均匀覆冰或不均匀覆冰的覆冰荷载，获得相应的若干第二耐张串轴向张力，绘制覆冰厚度-耐张串轴向张力的散点图；
曲线拟合模块，用于根据所述散点图，采用数据拟合方法，获取关于覆冰厚度与耐张串轴向张力的曲线拟合公式；
第四计算模块，用于当第三拉力传感器监测到的第三拉力监测值大于所述第五耐张串轴向张力时，将所述第三拉力监测值代入所述曲线拟合公式，得到所述输电线路的最终等值覆冰厚度；
所述第一计算模块还用于：
根据所述第一耐张串轴向张力计算得到输电线路的第一平均等效长度：
式中，F01为所述第一耐张串轴向张力，n0为输电线分裂数，ω0为输电线单位长度重量，θ为绝缘子串夹角；
根据所述第一耐张串轴向张力计算得到所述输电线路的第一最大等效长度：
式中，F01为所述第一耐张串轴向张力，n0为输电线分裂数，ω0为输电线单位长度重量，θ为绝缘子串夹角；
根据所述第一平均等效长度、所述第一最大等效长度和所述第三耐张串轴向张力计算得到所述输电线路的第一等值覆冰厚度：
式中，F1为所述第三耐张串轴向张力，ω0为输电线单位长度重量，n0为输电线分裂数，D为输电线外径，lav1为所述第一平均等效长度，lmax1为所述第一最大等效长度；
所述第一计算模块还用于：
根据所述第二耐张串轴向张力计算得到所述输电线路的第二平均等效长度：
式中，F02为所述第二耐张串轴向张力，n0为输电线分裂数，ω0为输电线单位长度重量，θ为绝缘子串夹角；
根据所述第二耐张串轴向张力计算得到所述输电线路的第二最大等效长度：
式中，F02为所述第二耐张串轴向张力，n0为输电线分裂数，ω0为输电线单位长度重量，θ为绝缘子串夹角；
根据所述第二平均等效长度、所述第二最大等效长度和所述第四耐张串轴向张力计算得到所述输电线路的第二等值覆冰厚度：
式中，F2为所述第四耐张串轴向张力，ω0为输电线单位长度重量，n0为输电线分裂数，D为输电线外径，lav2为所述第二平均等效长度，lmax2为所述第二最大等效长度；
所述基于所述第一等值覆冰厚度、所述第二等值覆冰厚度与所述覆冰厚度理论值的比较结果确定拉力传感器的布置位置包括：
当|b1-b0|＜|b2-b0|时，拉力传感器布置在耐张段首端的耐张塔上；
当|b1-b0|＞|b2-b0|时，拉力传感器布置在耐张段尾端的耐张塔上；
当|b1-b0|＝|b2-b0|时，拉力传感器布置在耐张段首端的耐张塔上或耐张段尾端的耐张塔上；
其中，b1为第一等值覆冰厚度，b2为第二等值覆冰厚度，b0为覆冰厚度理论值。
CN202110847432.6A 2021-07-26 2021-07-26 一种耐张塔输电线路等值覆冰厚度测量方法及系统 Active CN113701689B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110847432.6A CN113701689B (zh) 2021-07-26 2021-07-26 一种耐张塔输电线路等值覆冰厚度测量方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110847432.6A CN113701689B (zh) 2021-07-26 2021-07-26 一种耐张塔输电线路等值覆冰厚度测量方法及系统

Publications (2)

Publication Number Publication Date
CN113701689A CN113701689A (zh) 2021-11-26
CN113701689B true CN113701689B (zh) 2024-01-16

Family

ID=78650823

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110847432.6A Active CN113701689B (zh) 2021-07-26 2021-07-26 一种耐张塔输电线路等值覆冰厚度测量方法及系统

Country Status (1)

Country Link
CN (1) CN113701689B (zh)

Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2009049544A1 (fr) * 2007-10-12 2009-04-23 Lianggui Gong Dispositif de fusion de glace pour ligne de transmission à faisceaux et son procédé de réalisation
CN102507123A (zh) * 2011-09-29 2012-06-20 中国电力科学研究院 一种输电线路覆冰导线断线试验方法
CN102663215A (zh) * 2012-05-14 2012-09-12 重庆大学 架空输电线路塔-线耦合体系抗冰能力的评估方法
US20140067271A1 (en) * 2012-08-30 2014-03-06 International Business Machines Corporation Predicting ice coating status on transmission lines
CN104091046A (zh) * 2014-06-16 2014-10-08 中国南方电网有限责任公司超高压输电公司检修试验中心 一种基于称重法的耐张塔导线覆冰厚度计算方法和系统
CN104573192A (zh) * 2014-12-17 2015-04-29 华南理工大学 一种架空线路耐张塔等值覆冰厚度的在线监测方法
CN104655030A (zh) * 2015-02-16 2015-05-27 国网安徽省电力公司铜陵供电公司 一种输电线路覆冰检测与预警装置
CN110378066A (zh) * 2019-07-29 2019-10-25 贵州电网有限责任公司 一种耐张塔输电线路覆冰厚度计算方法
CN111442752A (zh) * 2020-03-26 2020-07-24 广州长川科技有限公司 一种输电线路等值覆冰厚度的监测方法
CN111998814A (zh) * 2020-08-26 2020-11-27 贵州电网有限责任公司 一种输电线路和高空输电线路覆冰厚度获取方法及计算仪
CN112329221A (zh) * 2020-11-02 2021-02-05 云南师范大学 输电线覆冰厚度确定方法

Patent Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2009049544A1 (fr) * 2007-10-12 2009-04-23 Lianggui Gong Dispositif de fusion de glace pour ligne de transmission à faisceaux et son procédé de réalisation
CN102507123A (zh) * 2011-09-29 2012-06-20 中国电力科学研究院 一种输电线路覆冰导线断线试验方法
CN102663215A (zh) * 2012-05-14 2012-09-12 重庆大学 架空输电线路塔-线耦合体系抗冰能力的评估方法
US20140067271A1 (en) * 2012-08-30 2014-03-06 International Business Machines Corporation Predicting ice coating status on transmission lines
CN104091046A (zh) * 2014-06-16 2014-10-08 中国南方电网有限责任公司超高压输电公司检修试验中心 一种基于称重法的耐张塔导线覆冰厚度计算方法和系统
CN104573192A (zh) * 2014-12-17 2015-04-29 华南理工大学 一种架空线路耐张塔等值覆冰厚度的在线监测方法
CN104655030A (zh) * 2015-02-16 2015-05-27 国网安徽省电力公司铜陵供电公司 一种输电线路覆冰检测与预警装置
CN110378066A (zh) * 2019-07-29 2019-10-25 贵州电网有限责任公司 一种耐张塔输电线路覆冰厚度计算方法
CN111442752A (zh) * 2020-03-26 2020-07-24 广州长川科技有限公司 一种输电线路等值覆冰厚度的监测方法
CN111998814A (zh) * 2020-08-26 2020-11-27 贵州电网有限责任公司 一种输电线路和高空输电线路覆冰厚度获取方法及计算仪
CN112329221A (zh) * 2020-11-02 2021-02-05 云南师范大学 输电线覆冰厚度确定方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
一种基于图像的覆冰厚度测量方法;邵前程等;《软件导刊》;第17卷(第4期);全文 *
输电线路覆冰监测新技术研究;朱传勇等;《地理空间信息》;第15卷(第6期);全文 *

Also Published As

Publication number Publication date
CN113701689A (zh) 2021-11-26

Similar Documents

Publication Publication Date Title
CN102507123B (zh) 一种输电线路覆冰导线断线试验方法
WO2015149411A1 (zh) 覆冰导线脱冰跳跃仿真测试方法
CN102692309B (zh) 一种应用于台风风场的输电线风洞试验系统及方法
CN104167076B (zh) 一种覆冰输电线路薄弱环节预警方法
CN104406777A (zh) 一种三维空间变角度加载装置
CN107907311A (zh) 一种输电线路导线覆冰不平衡张力试验方法及系统
CN104457594B (zh) 一种输电线路分布式覆冰监测方法
CN114896872A (zh) 一种高压输电线路覆冰状态综合评估方法
CN113701689B (zh) 一种耐张塔输电线路等值覆冰厚度测量方法及系统
CN101858945A (zh) 光纤复合架空地线（opgw）实验室模拟覆冰试验方法
Al Shboul et al. Intelligent approach for accurately predicting fatigue damage in overhead highway sign structures
Paluch et al. Experimental and numerical assessment of EPS wind action on long span transmission line conductors
CN102735468B (zh) 支座广义位移温度变化基于混合监测的受损索识别方法
CN108090263B (zh) 输电线脱冰振动缩尺试验气动阻尼的验证方法
CN113553747A (zh) 一种覆冰输电塔力学特性分析及状态评估方法
CN107358020B (zh) 一种弧垂过大配网导线紧线长度计算方法
CN103441450B (zh) 一种基于输电线路滑轮辗压铲刮的除冰方法
Solovev et al. Modeling of the mean wind loads on structures
CN116222466A (zh) 耐张塔导线覆冰厚度监测方法、装置、设备和介质
CN113704969B (zh) 一种直线塔地线等值覆冰厚度测量方法及系统
CN102735471B (zh) 支座广义位移温度变化基于应变监测的受损索识别方法
CN211824840U (zh) 超大钢结构合拢位置焊接残余应力监测系统
CN117494542A (zh) 一种基于激光点云数据和bp神经网络的输电线路直线塔线路覆冰厚度反演方法
CN110455511B (zh) 测定复合横担在断线工况下动态响应的试验装置及试验方法
CN113128125B (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