CN113570681B - 接地金属材料用土壤腐蚀等级分布图的绘制方法及系统 - Google Patents
接地金属材料用土壤腐蚀等级分布图的绘制方法及系统 Download PDFInfo
- Publication number
- CN113570681B CN113570681B CN202110805733.2A CN202110805733A CN113570681B CN 113570681 B CN113570681 B CN 113570681B CN 202110805733 A CN202110805733 A CN 202110805733A CN 113570681 B CN113570681 B CN 113570681B
- Authority
- CN
- China
- Prior art keywords
- soil
- corrosion
- index
- point
- grade
- 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
- 230000007797 corrosion Effects 0.000 title claims abstract description 212
- 238000005260 corrosion Methods 0.000 title claims abstract description 212
- 239000002689 soil Substances 0.000 title claims abstract description 182
- 238000010586 diagram Methods 0.000 title claims abstract description 48
- 238000000034 method Methods 0.000 title claims abstract description 41
- 238000009826 distribution Methods 0.000 title claims abstract description 39
- 239000007769 metal material Substances 0.000 title claims abstract description 33
- 238000012544 monitoring process Methods 0.000 claims abstract description 39
- 238000001514 detection method Methods 0.000 claims abstract description 25
- 238000013210 evaluation model Methods 0.000 claims abstract description 18
- 239000003086 colorant Substances 0.000 claims abstract description 9
- 230000005540 biological transmission Effects 0.000 claims abstract description 8
- 238000004040 coloring Methods 0.000 claims abstract description 8
- 230000001502 supplementing effect Effects 0.000 claims abstract description 8
- 238000012423 maintenance Methods 0.000 claims description 18
- 230000009466 transformation Effects 0.000 claims description 16
- 238000012937 correction Methods 0.000 claims description 10
- 238000004364 calculation method Methods 0.000 claims description 9
- 239000000463 material Substances 0.000 claims description 9
- 229910001335 Galvanized steel Inorganic materials 0.000 claims description 8
- 239000008397 galvanized steel Substances 0.000 claims description 8
- 238000005259 measurement Methods 0.000 claims description 8
- 238000012360 testing method Methods 0.000 claims description 8
- 238000004458 analytical method Methods 0.000 claims description 7
- 229910000975 Carbon steel Inorganic materials 0.000 claims description 4
- 239000010962 carbon steel Substances 0.000 claims description 4
- 238000006243 chemical reaction Methods 0.000 claims description 4
- 238000011156 evaluation Methods 0.000 claims description 4
- 238000012284 sample analysis method Methods 0.000 claims description 4
- 239000002184 metal Substances 0.000 claims description 3
- 238000013507 mapping Methods 0.000 claims description 2
- 238000005536 corrosion prevention Methods 0.000 description 4
- 238000009412 basement excavation Methods 0.000 description 3
- 238000007689 inspection Methods 0.000 description 3
- 230000002950 deficient Effects 0.000 description 2
- 230000009469 supplementation Effects 0.000 description 2
- 238000010276 construction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000004870 electrical engineering Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
- 238000004162 soil erosion Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/20—Drawing from basic elements, e.g. lines or circles
- G06T11/206—Drawing of charts or graphs
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/40—Filling a planar surface by adding surface attributes, e.g. colour or texture
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T3/00—Geometric image transformations in the plane of the image
- G06T3/40—Scaling of whole images or parts thereof, e.g. expanding or contracting
- G06T3/4023—Scaling of whole images or parts thereof, e.g. expanding or contracting based on decimating pixels or lines of pixels; based on inserting pixels or lines of pixels
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Testing Resistance To Weather, Investigating Materials By Mechanical Methods (AREA)
Abstract
本发明公开了一种接地金属材料用土壤腐蚀等级分布图的绘制方法及系统。本发明的绘制方法包括：构建相应的评估模型，采集准确的土壤腐蚀速率；利用监测点已获得的土壤腐蚀速率，通过空间插值算法进行插值，得到非监测区域的土壤腐蚀速率；根据腐蚀速率进行土壤腐蚀性等级划分；结合收集的实际运行信息修正和补充腐蚀等级数据，当实际运行信息和通过检测得到的腐蚀等级两者不一致时，按运行经验信息确定腐蚀等级；在素色地理图上根据划定的不同等级范围分别用不同的鲜明颜色予以着色，同时标识变电站、换流站、发电厂、输电线路和各监测点符号及名称。本发明绘制的土壤腐蚀等级分布图准确度高，能简单、清晰地展示区域内的土壤腐蚀性等级。
Description
技术领域
本发明属于电网接地金属材料防腐蚀技术领域，具体地说是一种电力工程接地金属材料用土壤腐蚀等级分布图的绘制方法及系统。
背景技术
电力接地金属设备的自然环境腐蚀规律与特征具有明显的地域差异性，接地材料长期深埋于土壤环境中，不同土壤环境的金属接地材料腐蚀速率及使用寿命往往差别很大。
现阶段，电力工程所处土壤环境腐蚀性及相关接地金属材料腐蚀数据十分欠缺，部分地区腐蚀环境数据缺失严重，导致电网工程建设中无法根据土壤腐蚀环境进行合理选址、选材和个性化防腐设计，接地金属材料腐蚀寿命的不足给电网安全造成很大的风险，且大幅增加运行后续维护成本。大型接地网工程应是免维护的，提前根据土壤腐蚀性等级特点进行针对性选址、选材，是降低接地工程维护成本的关键举措。
目前，电力企业生产运维单位中缺乏防腐专业人员及技术支撑力量，难以通过测量及查询零散资料的方式获得当地土壤腐蚀性等级，也就无法选择接地工程的个性化防腐运维方案，故而防腐效率低下。
显而易见，目前还没有一种比较合理的绘制电力工程接地金属材料用土壤腐蚀等级分布图的方法。
因此，开发出一种准确度高、合理的绘制电力工程接地金属材料用土壤腐蚀等级分布图的方法，将以往分散的土壤腐蚀性数据以图表直观化显示。供给接地工程设计人员、运维管理人员查询使用，将大幅提升工作效率，具有明显的工程实际价值。
发明内容
本发明所要解决的技术问题是提供一种准确度高、合理可行的绘制接地金属材料用土壤腐蚀等级分布图的方法及系统，以克服现有技术存在的缺陷和满足实际情况的需求。
为此，本发明采用的一种技术方案为：接地金属材料用土壤腐蚀等级分布图的绘制方法，其包括以下步骤：
基于土壤环境对金属材料腐蚀的检测方法，针对土样分析检测法构建灰色星座聚类评估模型，采集准确的土壤腐蚀速率；
利用监测点已获得的土壤腐蚀速率，通过空间插值算法进行插值，得到非监测区域的土壤腐蚀速率；
根据腐蚀速率进行土壤腐蚀性等级划分；
结合收集的实际运行信息修正和补充腐蚀等级数据，当实际运行信息和通过检测得到的腐蚀等级两者不一致时，按实际运行信息确定腐蚀等级；
在素色地理图上根据划定的不同等级范围分别用不同的鲜明颜色予以着色，同时标识变电站、换流站、发电厂、输电线路和各监测点符号及名称，进而绘制电力工程接地金属材料用土壤腐蚀等级分布图。
进一步的，所述的针对土样分析检测法的灰色星座聚类评估模型，由于土壤分布不均匀的情况，在某一区域范围内，在对角线方向均匀采集多个土样，按照DL/T 1554-2016《接地网土壤腐蚀性评价导则》中的三指标法或八指标法分析土壤腐蚀特性，然后利用灰色星座聚类模型，获得土样的聚类中心，进而得到能够代表该区域的土壤腐蚀指标值，再根据指标值对应的腐蚀速率范围内任选一腐蚀速率值参与空间插值算法。
进一步地，所述的灰色星座聚类评估模型，将每个土样的指标测量值，经过原始数据的极差变换等手段转化为无纲量的坐标系，一个土样点用一颗“星点”表示，计算星座图中距离其余各点距离和最小的一点作为聚类中心，用该点的测试数据代表该区域的土壤腐蚀数据；
对土样的每个指标值进行极差变换，并使变换后的数值均落在[0°，180°]的闭区间内，计算公式为：
式中：aij为变换后的数据，以角度表示；Xij为原始数据；Xjmax为第j个指标的最大值；Xjmin第j个指标的最小值；i为样本下标，j为指标下标；
对每个指标，根据其对系统变化的影响程度，分别给一个权重Wj,使
式中：Wj为第j个指标的权数，0<Wj<1；
对各个指标进行直角坐标计算：利用极坐标与直角坐标的变换关系，先求出每个点各项指标的Xi，Yi值，然后将各点每项指标的Xi，Yi值相加，即为各土样点的坐标值，其变换公式为：
式中：Xi为第i个样点的横坐标；Yi为第i个样点的纵坐标；
绘制一个半径为1的上半圆，以圆中心为坐标原点，以上半圆底为横坐标X轴，并作出过原点的Y轴，根据Xi，Yi的值确定每一个样点在星座图内的位置，将性质相似和接近的样点聚在一起，形成一个“星座”；
计算星座图中距离其余各点距离和最小的一点
式中：Xi为第i个样点的横坐标；Yi为第i个样点的纵坐标；X,Y为星座图中一点；
将上式求得的样点作为聚类中心，用该点的腐蚀测量数据代表该区域的土壤腐蚀速率并参与空间插值算法。
进一步地，所述的空间插值算法包括将监测点检测所得的土壤腐蚀速率数据及监测点地理位置信息代入公式：
式中，Z是插值点估计值；Zi是第i个样本点观测值；λi是第i个样本点观测值的修正系数，其值为1/(1+Hi)；Hi为海拔高度，单位为km；Di是插值点与第i个样本点之间的距离；n是用于估算插值点的样点数；p值为2；
其中R为地球半径，插值点的纬度角为β1，经度角为/>观测点的纬度角为β2，经度角为/>
进一步地，所述的土壤腐蚀性等级根据DL/T 1554-2016《接地网土壤腐蚀性评价导则》中分级标准划分，以碳钢材料的腐蚀速率rcorr进行土壤腐蚀性等级划分，分为：微，rcorr＜1；弱，1≤rcorr＜5；中，5≤rcorr≤7；强，rcorr＞7；rcorr的单位为g/(dm2*a)。
进一步地，所述的结合收集的实际运行信息修正和补充腐蚀等级数据，包括：对新建镀锌钢接地装置，在5年内开挖检查即发现腐蚀或在5年以上且10年内开挖检查即发现明显腐蚀的地区，判断为强腐蚀环境；对新建镀锌钢接地装置，在10年以上且15年内开挖检查即发现明显腐蚀的地区，判断为中腐蚀环境；统计接地网维修记录，维修周期小于15年的地区，判断为强腐蚀环境；维修周期大于等于15年且小于20年的地区，判断为中腐蚀环境；近三年曾发生接地事故的电力工程，其所在地的土壤腐蚀性认定为强；近三年测试发现接地电阻和跨步电压有明显上升趋势的地区，其土壤腐蚀性认定为强。
进一步地，当实际运行信息得到的腐蚀等级与插值法得到的腐蚀等级只差别一个等级时，按照实际运行信息直接进行等级修正；如果两者相差在2个等级以上时，进行腐蚀等级校核。
本发明采用的另一种技术方案为：接地金属材料用土壤腐蚀等级分布图的绘制系统，其包括：
土壤腐蚀速率采集单元，基于土壤环境对金属材料腐蚀的检测方法，针对土样分析检测法构建灰色星座聚类评估模型，采集准确的土壤腐蚀速率；
非监测区域土壤腐蚀速率的获取单元，利用监测点已获得的土壤腐蚀速率，通过空间插值算法进行插值，得到非监测区域的土壤腐蚀速率；
土壤腐蚀性等级划分单元，根据腐蚀速率进行土壤腐蚀性等级划分；
腐蚀等级确定单元，结合收集的实际运行信息修正和补充腐蚀等级数据，当实际运行信息和通过检测得到的腐蚀等级两者不一致时，按实际运行信息确定腐蚀等级；
土壤腐蚀等级分布图绘制单元，在素色地理图上根据划定的不同等级范围分别用不同的鲜明颜色予以着色，同时标识变电站、换流站、发电厂、输电线路和各监测点符号及名称，进而绘制电力工程接地金属材料用土壤腐蚀等级分布图。
本发明的绘制方法及系统合理且准确度高，将以往分散的土壤腐蚀性数据以图表直观化显示土壤腐蚀等级，使用起来简单清晰；供给接地工程设计人员、运维管理人员查询使用，将大幅提升防腐效率，避免大量电网腐蚀故障的产生、降低设备运维成本，具有明显的工程实际价值。
附图说明
图1是本发明实施例1中某地区土壤腐蚀等级分布图的绘制流程图；
图2是本发明实施例1中所得某地区土壤腐蚀等级分布图；
图3为本发明实施例2中土壤腐蚀等级分布图绘制系统的结构图。
具体实施方式
下面结合附图并通过实施例对本发明作进一步的详细说明，以下实施例是对本发明的解释而本发明并不局限于以下实施例。
实施例1
参照图1所示，为本发明提供的一种电力工程接地金属材料用土壤腐蚀等级分布图的绘制方法的一个实施例的流程图，包括：构建相应的评估模型，采集准确的土壤腐蚀速率；利用监测点已获得的土壤腐蚀速率，通过空间插值算法进行插值，得到非监测区域的土壤腐蚀速率；根据腐蚀速率进行土壤腐蚀性等级划分；结合收集的实际运行信息修正和补充腐蚀等级数据，当实际运行信息和通过检测得到的腐蚀等级两者不一致时，按运行经验信息确定腐蚀等级；在素色地理图上根据划定的不同等级范围分别用不同的鲜明颜色予以着色，同时标识变电站、换流站、发电厂、输电线路和各监测点符号及名称。本发明绘制的土壤腐蚀等级分布图准确度高，能简单、清晰地展示区域内的土壤腐蚀性等级。
在某地区布设28个监测点，其中27个监测点采用埋片法检测土壤腐蚀特性，1个监测点采用土样分析法检测土壤腐蚀特性。
表1 27组监测点的腐蚀数据及腐蚀等级统计表
针对土样分析检测法构建灰色聚类评估模型，由于土壤分布不均匀的情况，在0.01km2范围内，在对角线方向均匀采集9个土样，按照八指标法分析土壤腐蚀特性，然后利用灰色聚类模型，获得每个指标的聚类中心，进而得到能够代表该区域的土壤腐蚀参数；
表2不同样本各个指标的原始数据
对土样的每个指标值进行极差变换，并使变换后的数值均落在[0°，180°]的闭区间内。
表3不同样本各个指标的极差变换值
对各个指标进行直角坐标计算。利用极坐标与直角坐标的变换关系，先求出每个点各项指标的Xi，Yi值，然后将各点每项指标的Xi，Yi值相加，即为各土样点的坐标值。
表4直角坐标转换
计算星座图中距离其余各点距离和最小的一点，作为聚类中心，由表5可以得知第8个样本为该组数据的聚类中心。八指标法得到的指标数值为3，对应腐蚀速率范围为3≤rcorr＜5，取中间值rcorr为4作为本监测点的土壤腐蚀速率并参与空间插值算法。
表5各样本点与其余各点距离之和
利用28个监测点已获得的土壤腐蚀速率，通过空间插值算法进行插值，得到非监测区域的土壤腐蚀速率。将收集、检测所得土壤腐蚀速率数据及监测点地理位置信息代入公式：
式中，Z是插值点估计值；Zi是第i个样本点观测值；λi是第i个样本点观测值的修正系数，其值为1/(1+Hi)；Hi为海拔高度，单位为km；Di是插值点与第i个样本点之间的距离；n是用于估算插值点的样点数；p值通常被默认为2；其中R为地球半径，插值点的纬度角为β1，经度角为/>观测点的纬度角为β2，经度角为/>
然后，根据碳钢材料的腐蚀速率rcorr(g/(dm2*a))进行土壤腐蚀性等级划分，分为微(rcorr＜1)、弱(1≤rcorr＜5)、中(5≤rcorr≤7)、强(rcorr＞7)。
最后，结合收集的实际运行信息修正、补充腐蚀等级数据，包括对新建镀锌钢接地装置，在5年内开挖检查即发现腐蚀或在5年以上且10年内开挖检查即发现呈明显腐蚀的地区，可判断为强腐蚀环境；对新建镀锌钢接地装置，在10年以上且15年内开挖检查即发现呈明显腐蚀的地区，可判断为中腐蚀环境；统计接地网维修记录，维修周期小于15年的地区，可判断为强腐蚀环境；维修周期大于等于15年且小于20年的地区，可判断为中腐蚀环境；近三年曾发生接地事故的电力工程，其所在地的土壤腐蚀性应认定为强；近三年测试发现接地电阻、跨步电压等电学参数有明显上升趋势的地区，其土壤腐蚀性应认定为强。当实际运行信息得到的腐蚀等级与插值法得到的腐蚀等级只差别一个等级时，可按照实际运行信息直接进行等级修正，如果两者相差在2个等级以上时，应进行该点腐蚀等级校核。
在素色地理图上根据划定的不同等级范围分别用不同的鲜明颜色予以着色，同时标识变电站、换流站、发电厂、输电线路和各监测点符号及名称，进而绘制该地区电力工程接地金属材料用土壤腐蚀等级分布图。本实施例所用腐蚀等级对应颜色如表6所示，绘制的电力工程接地金属材料用土壤腐蚀等级分布图如附图2所示。
表6土壤腐蚀等级对应颜色
实施例2
一种电力工程接地金属材料用土壤腐蚀等级分布图的绘制系统，其包括：
土壤腐蚀速率采集单元，基于土壤环境对金属材料腐蚀的检测方法，针对土样分析检测法构建灰色星座聚类评估模型，采集准确的土壤腐蚀速率；
非监测区域土壤腐蚀速率的获取单元，利用监测点已获得的土壤腐蚀速率，通过空间插值算法进行插值，得到非监测区域的土壤腐蚀速率；
土壤腐蚀性等级划分单元，根据腐蚀速率进行土壤腐蚀性等级划分；
腐蚀等级确定单元，结合收集的实际运行信息修正和补充腐蚀等级数据，当实际运行信息和通过检测得到的腐蚀等级两者不一致时，按实际运行信息确定腐蚀等级；
土壤腐蚀等级分布图绘制单元，在素色地理图上根据划定的不同等级范围分别用不同的鲜明颜色予以着色，同时标识变电站、换流站、发电厂、输电线路和各监测点符号及名称，进而绘制电力工程接地金属材料用土壤腐蚀等级分布图。
所述的针对土样分析检测法的灰色星座聚类评估模型，由于土壤分布不均匀的情况，在某一区域范围内，在对角线方向均匀采集多个土样，按照DL/T 1554-2016《接地网土壤腐蚀性评价导则》中的三指标法或八指标法分析土壤腐蚀特性，然后利用灰色星座聚类模型，获得土样的聚类中心，进而得到能够代表该区域的土壤腐蚀指标值，再根据指标值对应的腐蚀速率范围内任选一腐蚀速率值参与空间插值算法。
所述的灰色星座聚类评估模型，将每个土样的指标测量值，经过原始数据的极差变换等手段转化为无纲量的坐标系，一个土样点用一颗“星点”表示，计算星座图中距离其余各点距离和最小的一点作为聚类中心，用该点的测试数据代表该区域的土壤腐蚀数据；
对土样的每个指标值进行极差变换，并使变换后的数值均落在[0°，180°]的闭区间内，计算公式为：
式中：aij为变换后的数据，以角度表示；Xij为原始数据；Xjmax为第j个指标的最大值；Xjmin第j个指标的最小值；i为样本下标，j为指标下标；
对每个指标，根据其对系统变化的影响程度，分别给一个权重Wj,使
式中：Wj为第j个指标的权数，0<Wj<1；
对各个指标进行直角坐标计算：利用极坐标与直角坐标的变换关系，先求出每个点各项指标的Xi，Yi值，然后将各点每项指标的Xi，Yi值相加，即为各土样点的坐标值，其变换公式为：
式中：Xi为第i个样点的横坐标；Yi为第i个样点的纵坐标；
绘制一个半径为1的上半圆，以圆中心为坐标原点，以上半圆底为横坐标X轴，并作出过原点的Y轴，根据Xi，Yi的值确定每一个样点在星座图内的位置，将性质相似和接近的样点聚在一起，形成一个“星座”；
计算星座图中距离其余各点距离和最小的一点
式中：Xi为第i个样点的横坐标；Yi为第i个样点的纵坐标；X,Y为星座图中一点；
将上式求得的样点作为聚类中心，用该点的腐蚀测量数据代表该区域的土壤腐蚀速率并参与空间插值算法。
所述的空间插值算法包括将监测点检测所得的土壤腐蚀速率数据及监测点地理位置信息代入公式：
式中，Z是插值点估计值；Zi是第i个样本点观测值；λi是第i个样本点观测值的修正系数，其值为1/(1+Hi)；Hi为海拔高度，单位为km；Di是插值点与第i个样本点之间的距离；n是用于估算插值点的样点数；p值为2；
其中R为地球半径，插值点的纬度角为β1，经度角为/>观测点的纬度角为β2，经度角为/>
所述的土壤腐蚀性等级划分包括根据碳钢材料的腐蚀速率rcorr进行土壤腐蚀性等级划分，分为：微，rcorr＜1；弱，1≤rcorr＜5；中，5≤rcorr≤7；强，rcorr＞7；rcorr的单位为g/(dm2*a)。
所述的结合收集的实际运行信息修正和补充腐蚀等级数据，包括：对新建镀锌钢接地装置，在5年内开挖检查即发现腐蚀或在5年以上且10年内开挖检查即发现明显腐蚀的地区，判断为强腐蚀环境；对新建镀锌钢接地装置，在10年以上且15年内开挖检查即发现明显腐蚀的地区，判断为中腐蚀环境；统计接地网维修记录，维修周期小于15年的地区，判断为强腐蚀环境；维修周期大于等于15年且小于20年的地区，判断为中腐蚀环境；近三年曾发生接地事故的电力工程，其所在地的土壤腐蚀性认定为强；近三年测试发现接地电阻和跨步电压有明显上升趋势的地区，其土壤腐蚀性认定为强。
当实际运行信息得到的腐蚀等级与插值法得到的腐蚀等级只差别一个等级时，按照实际运行信息直接进行等级修正；如果两者相差在2个等级以上时，进行腐蚀等级校核。
以上所述的仅为本发明的具体实施例，并不以此限定本发明的保护范围，应当指出，对于本领域的技术人员来说，在不违反本发明构思的基础上所作的任何替换与改进，均属本发明的保护范围。
Claims (8)
1.接地金属材料用土壤腐蚀等级分布图的绘制方法，其特征在于，包括以下步骤：
基于土壤环境对金属材料腐蚀的检测方法，针对土样分析检测法构建灰色星座聚类评估模型，采集准确的土壤腐蚀速率；
利用监测点已获得的土壤腐蚀速率，通过空间插值算法进行插值，得到非监测区域的土壤腐蚀速率；
根据腐蚀速率进行土壤腐蚀性等级划分；
结合收集的实际运行信息修正和补充腐蚀等级数据，当实际运行信息和通过检测得到的腐蚀等级两者不一致时，按实际运行信息确定腐蚀等级；
在素色地理图上根据划定的不同等级范围分别用不同的鲜明颜色予以着色，同时标识变电站、换流站、发电厂、输电线路和各监测点符号及名称，进而绘制电力工程接地金属材料用土壤腐蚀等级分布图；
所述的灰色星座聚类评估模型，将每个土样的指标测量值，经过原始数据的极差变换手段转化为无纲量的坐标系，一个土样点用一颗“星点”表示，计算星座图中距离其余各点距离和最小的一点作为聚类中心，用该点的测试数据代表该区域的土壤腐蚀数据；
对土样的每个指标值进行极差变换，并使变换后的数值均落在[0°，180°]的闭区间内，计算公式为：
式中：aij为变换后的数据，以角度表示；Xij为原始数据；Xjmax为第j个指标的最大值；Xjmin第j个指标的最小值；i为样本下标，j为指标下标；
对每个指标，根据其对系统变化的影响程度，分别给一个权重Wj,使
式中：Wj为第j个指标的权数，0<Wj<1；
对各个指标进行直角坐标计算：利用极坐标与直角坐标的变换关系，先求出每个点各项指标的Xi，Yi值，然后将各点每项指标的Xi，Yi值相加，即为各土样点的坐标值，其变换公式为：
式中：Xi为第i个样点的横坐标；Yi为第i个样点的纵坐标；
绘制一个半径为1的上半圆，以圆中心为坐标原点，以上半圆底为横坐标X轴，并作出过原点的Y轴，根据Xi，Yi的值确定每一个样点在星座图内的位置，将性质相似和接近的样点聚在一起，形成一个“星座”；
计算星座图中距离其余各点距离和最小的一点
式中：Xi为第i个样点的横坐标；Yi为第i个样点的纵坐标；X,Y为星座图中一点；
将上式求得的样点作为聚类中心，用该点的腐蚀测量数据代表该区域的土壤腐蚀速率并参与空间插值算法。
2.根据权利要求1所述的接地金属材料用土壤腐蚀等级分布图的绘制方法，其特征在于，针对土样分析检测法，由于土壤分布不均匀的情况，在某一区域范围内，在对角线方向均匀采集多个土样，按照DL/T 1554-2016《接地网土壤腐蚀性评价导则》中的三指标法或八指标法分析土壤腐蚀特性，然后利用灰色星座聚类评估模型，获得土样的聚类中心，进而得到能够代表该区域的土壤腐蚀指标值，再根据土壤腐蚀指标值对应的腐蚀速率范围内任选一腐蚀速率值参与空间插值算法。
3.根据权利要求1或2所述的接地金属材料用土壤腐蚀等级分布图的绘制方法，其特征在于，所述的空间插值算法包括将监测点检测所得的土壤腐蚀速率数据及监测点地理位置信息代入公式：
式中，Z是插值点估计值；Zi是第i个样本点观测值；λi是第i个样本点观测值的修正系数，其值为1/(1+Hi)；Hi为海拔高度，单位为km；Di是插值点与第i个样本点之间的距离；n是用于估算插值点的样点数；p值为2；
其中R为地球半径，插值点的纬度角为β1，经度角为/>观测点的纬度角为β2，经度角为/>
4.根据权利要求1或2所述的接地金属材料用土壤腐蚀等级分布图的绘制方法，其特征在于，所述的土壤腐蚀性等级根据DL/T 1554-2016《接地网土壤腐蚀性评价导则》中分级标准划分，以碳钢材料的腐蚀速率rcorr进行土壤腐蚀性等级划分，分为：微，rcorr＜1；弱，1≤rcorr＜5；中，5≤rcorr≤7；强，rcorr＞7；rcorr的单位为g/(dm2*a)。
5.根据权利要求1或2所述的接地金属材料用土壤腐蚀等级分布图的绘制方法，其特征在于，所述的结合收集的实际运行信息修正和补充腐蚀等级数据，包括：对新建镀锌钢接地装置，在5年内开挖检查即发现腐蚀或在5年以上且10年内开挖检查即发现明显腐蚀的地区，判断为强腐蚀环境；对新建镀锌钢接地装置，在10年以上且15年内开挖检查即发现明显腐蚀的地区，判断为中腐蚀环境；统计接地网维修记录，维修周期小于15年的地区，判断为强腐蚀环境；维修周期大于等于15年且小于20年的地区，判断为中腐蚀环境；近三年曾发生接地事故的电力工程，其所在地的土壤腐蚀性认定为强；近三年测试发现接地电阻和跨步电压有明显上升趋势的地区，其土壤腐蚀性认定为强。
6.根据权利要求1或2所述的接地金属材料用土壤腐蚀等级分布图的绘制方法，其特征在于，当实际运行信息得到的腐蚀等级与插值法得到的腐蚀等级只差别一个等级时，按照实际运行信息直接进行等级修正；如果两者相差在2个等级以上时，进行腐蚀等级校核。
7.接地金属材料用土壤腐蚀等级分布图的绘制系统，其特征在于，包括以下步骤：
步骤1：土壤腐蚀速率采集单元，基于土壤环境对金属材料腐蚀的检测方法，针对土样分析检测法构建灰色星座聚类评估模型，采集准确的土壤腐蚀速率；
步骤2：非监测区域土壤腐蚀速率的获取单元，利用监测点已获得的土壤腐蚀速率，通过空间插值算法进行插值，得到非监测区域的土壤腐蚀速率；
步骤3：土壤腐蚀性等级划分单元，根据腐蚀速率进行土壤腐蚀性等级划分；
步骤4：腐蚀等级确定单元，结合收集的实际运行信息修正和补充腐蚀等级数据，当实际运行信息和通过检测得到的腐蚀等级两者不一致时，按实际运行信息确定腐蚀等级；
步骤5：土壤腐蚀等级分布图绘制单元，在素色地理图上根据划定的不同等级范围分别用不同的鲜明颜色予以着色，同时标识变电站、换流站、发电厂、输电线路和各监测点符号及名称，进而绘制电力工程接地金属材料用土壤腐蚀等级分布图；
所述的灰色星座聚类评估模型，将每个土样的指标测量值，经过原始数据的极差变换手段转化为无纲量的坐标系，一个土样点用一颗“星点”表示，计算星座图中距离其余各点距离和最小的一点作为聚类中心，用该点的测试数据代表该区域的土壤腐蚀数据；
对土样的每个指标值进行极差变换，并使变换后的数值均落在[0°，180°]的闭区间内，计算公式为：
式中：aij为变换后的数据，以角度表示；Xij为原始数据；Xjmax为第j个指标的最大值；Xjmin第j个指标的最小值；i为样本下标，j为指标下标；
对每个指标，根据其对系统变化的影响程度，分别给一个权重Wj,使
式中：Wj为第j个指标的权数，0<Wj<1；
对各个指标进行直角坐标计算：利用极坐标与直角坐标的变换关系，先求出每个点各项指标的xi，Yi值，然后将各点每项指标的Xi，Yi值相加，即为各土样点的坐标值，其变换公式为：
式中：Xi为第i个样点的横坐标；Yi为第i个样点的纵坐标；
绘制一个半径为1的上半圆，以圆中心为坐标原点，以上半圆底为横坐标X轴，并作出过原点的Y轴，根据Xi，Yi的值确定每一个样点在星座图内的位置，将性质相似和接近的样点聚在一起，形成一个“星座”；
计算星座图中距离其余各点距离和最小的一点
式中：Xi为第i个样点的横坐标；Yi为第i个样点的纵坐标；X,Y为星座图中一点；
将上式求得的样点作为聚类中心，用该点的腐蚀测量数据代表该区域的土壤腐蚀速率并参与空间插值算法。
8.根据权利要求7所述的接地金属材料用土壤腐蚀等级分布图的绘制系统，其特征在于，所述的针对土样分析检测法的灰色星座聚类评估模型，由于土壤分布不均匀的情况，在某一区域范围内，在对角线方向均匀采集多个土样，按照DL/T 1554-2016《接地网土壤腐蚀性评价导则》中的三指标法或八指标法分析土壤腐蚀特性，然后利用灰色星座聚类评估模型，获得土样的聚类中心，进而得到能够代表该区域的土壤腐蚀指标值，再根据指标值对应的腐蚀速率范围内任选一腐蚀速率值参与空间插值算法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110805733.2A CN113570681B (zh) | 2021-07-16 | 2021-07-16 | 接地金属材料用土壤腐蚀等级分布图的绘制方法及系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110805733.2A CN113570681B (zh) | 2021-07-16 | 2021-07-16 | 接地金属材料用土壤腐蚀等级分布图的绘制方法及系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113570681A CN113570681A (zh) | 2021-10-29 |
| CN113570681B true CN113570681B (zh) | 2023-09-12 |
Family
ID=78165134
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110805733.2A Active CN113570681B (zh) | 2021-07-16 | 2021-07-16 | 接地金属材料用土壤腐蚀等级分布图的绘制方法及系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113570681B (zh) |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102156089A (zh) * | 2011-01-18 | 2011-08-17 | 中国石油天然气股份有限公司 | 一种埋地管道内腐蚀评价方法 |
| CN103164619A (zh) * | 2013-03-06 | 2013-06-19 | 陕西电力科学研究院 | 基于最佳鉴别准则的变电站接地网腐蚀率模式识别方法 |
| CN107704669A (zh) * | 2017-09-20 | 2018-02-16 | 中国电力科学研究院 | 基于电化学极化曲线测试的接地极动态腐蚀过程仿真方法 |
| DE102018000418A1 (de) * | 2018-01-20 | 2019-07-25 | Bürkle Consulting Gmbh | Mechanochemisches Verfahren zur Herstellung von von persistenten organischen Schadstoffen und anderen Organohalogenverbindungen freien Wertprodukten aus Abfällen von Kunststoffen und Kunststofflaminaten |
| CN111341396A (zh) * | 2020-04-15 | 2020-06-26 | 国网浙江省电力有限公司电力科学研究院 | 一种大气环境对材料腐蚀安全评估方法及系统 |
| CN112251756A (zh) * | 2020-09-28 | 2021-01-22 | 北京科技大学 | 一种埋地金属管道动态直流腐蚀风险的评判系统和方法 |
| CN112633553A (zh) * | 2020-11-27 | 2021-04-09 | 合肥泽众城市智能科技有限公司 | 一种用于燃气管线-危化企业耦合隐患辨识与风险评估方法及系统 |
-
2021
- 2021-07-16 CN CN202110805733.2A patent/CN113570681B/zh active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102156089A (zh) * | 2011-01-18 | 2011-08-17 | 中国石油天然气股份有限公司 | 一种埋地管道内腐蚀评价方法 |
| CN103164619A (zh) * | 2013-03-06 | 2013-06-19 | 陕西电力科学研究院 | 基于最佳鉴别准则的变电站接地网腐蚀率模式识别方法 |
| CN107704669A (zh) * | 2017-09-20 | 2018-02-16 | 中国电力科学研究院 | 基于电化学极化曲线测试的接地极动态腐蚀过程仿真方法 |
| DE102018000418A1 (de) * | 2018-01-20 | 2019-07-25 | Bürkle Consulting Gmbh | Mechanochemisches Verfahren zur Herstellung von von persistenten organischen Schadstoffen und anderen Organohalogenverbindungen freien Wertprodukten aus Abfällen von Kunststoffen und Kunststofflaminaten |
| CN111341396A (zh) * | 2020-04-15 | 2020-06-26 | 国网浙江省电力有限公司电力科学研究院 | 一种大气环境对材料腐蚀安全评估方法及系统 |
| CN112251756A (zh) * | 2020-09-28 | 2021-01-22 | 北京科技大学 | 一种埋地金属管道动态直流腐蚀风险的评判系统和方法 |
| CN112633553A (zh) * | 2020-11-27 | 2021-04-09 | 合肥泽众城市智能科技有限公司 | 一种用于燃气管线-危化企业耦合隐患辨识与风险评估方法及系统 |
Non-Patent Citations (1)
| Title |
|---|
| 接地网腐蚀性评价方法与腐蚀速率预测;陈敬友等;《腐蚀与防护》;第42卷(第3期);第64-67页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113570681A (zh) | 2021-10-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Rhoades et al. | Adjustment of temperature and rainfall records for site changes | |
| Gugliani et al. | New methods to assess wind resources in terms of wind speed, load, power and direction | |
| CN105095589B (zh) | 一种山区电网风区分布图绘制方法 | |
| CN107391951B (zh) | 一种基于环形邻域梯度排序的空气污染溯源方法 | |
| CN114444259B (zh) | 一种雨污管网溯源追踪系统及方法 | |
| CN112685965B (zh) | 一种台风中输电杆塔风险监测方法及系统 | |
| CN111428942B (zh) | 基于变网格技术提取微地形因子的线路覆冰厚度预测方法 | |
| CN107677614A (zh) | 一种水中重金属污染风险在线预警系统及方法 | |
| CN114693281B (zh) | 一种基于云平台的工程勘察信息管理系统 | |
| CN116910041B (zh) | 一种基于尺度分析的遥感降水产品的逐日订正方法 | |
| CN113282576B (zh) | 一种气象数据质量控制方法 | |
| CN114637806A (zh) | 热带气旋预报检验指标数据的可视化分析方法 | |
| CN113935533B (zh) | 一种黄渤海海区大风推算方法 | |
| CN113570681B (zh) | 接地金属材料用土壤腐蚀等级分布图的绘制方法及系统 | |
| CN114398423A (zh) | 一种基于多源数据的河涌水质时空预测方法及系统 | |
| CN112257329A (zh) | 一种判定台风对线路影响的方法 | |
| CN113821895B (zh) | 输电线路覆冰厚度预测模型的构建方法、装置及存储介质 | |
| CN113610436B (zh) | 一种承灾体动态脆弱性评估方法及系统 | |
| CN113742929B (zh) | 一种针对格点气象实况的数据质量评价方法 | |
| CN115693666A (zh) | 基于卫星反演的海上风电场发电量确定方法及系统 | |
| CN115292425A (zh) | 一种考虑微环境影响的大气腐蚀等级地图绘制方法 | |
| CN111797545B (zh) | 一种基于实测数据的风电机组偏航折减系数计算方法 | |
| CN110751398B (zh) | 一种区域生态质量评价方法及装置 | |
| CN113361114A (zh) | 一种基于径流路径的多尺度面源污染物入河系数测算方法 | |
| CN113283100B (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 | ||
| CB03 | Change of inventor or designer information |
Inventor after: Liu Sen Inventor after: Hu Jiayuan Inventor after: Li Yanwei Inventor after: Qian Zhouhai Inventor before: Liu Sen Inventor before: Zheng Zhiming Inventor before: Hu Jiayuan Inventor before: Li Yanwei Inventor before: Qian Zhouhai Inventor before: Li Zhiguo Inventor before: Yu Zhiyong Inventor before: Zhu Xiaofeng Inventor before: Peng Yongliang Inventor before: Hu Wei |
|
| CB03 | Change of inventor or designer information | ||
| TA01 | Transfer of patent application right |
Effective date of registration: 20220107 Address after: The eight district of Hangzhou city in Zhejiang province 310014 Huadian Zhaohui under No. 1 Lane Applicant after: STATE GRID ZHEJIANG ELECTRIC POWER COMPANY LIMITED ELECTRIC POWER Research Institute Address before: The eight district of Hangzhou city in Zhejiang province 310014 Huadian Zhaohui under No. 1 Lane Applicant before: STATE GRID ZHEJIANG ELECTRIC POWER COMPANY LIMITED ELECTRIC POWER Research Institute Applicant before: HANGZHOU E ENERGY ELECTRIC POWER TECHNOLOGY Co.,Ltd. Applicant before: SHANGHAI JIUREN INFORMATION TECHNOLOGY CO.,LTD. |
|
| TA01 | Transfer of patent application right | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |