CN108020501A - 一种大气腐蚀等级图绘制方法 - Google Patents
一种大气腐蚀等级图绘制方法 Download PDFInfo
- Publication number
- CN108020501A CN108020501A CN201711131439.8A CN201711131439A CN108020501A CN 108020501 A CN108020501 A CN 108020501A CN 201711131439 A CN201711131439 A CN 201711131439A CN 108020501 A CN108020501 A CN 108020501A
- Authority
- CN
- China
- Prior art keywords
- atmospheric corrosion
- grade
- atmospheric
- data
- mrow
- 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.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N17/00—Investigating resistance of materials to the weather, to corrosion, or to light
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01D—MEASURING NOT SPECIALLY ADAPTED FOR A SPECIFIC VARIABLE; ARRANGEMENTS FOR MEASURING TWO OR MORE VARIABLES NOT COVERED IN A SINGLE OTHER SUBCLASS; TARIFF METERING APPARATUS; MEASURING OR TESTING NOT OTHERWISE PROVIDED FOR
- G01D21/00—Measuring or testing not otherwise provided for
- G01D21/02—Measuring two or more variables by means not covered by a single other subclass
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Analytical Chemistry (AREA)
- Environmental Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- Health & Medical Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Ecology (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Biodiversity & Conservation Biology (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Testing Resistance To Weather, Investigating Materials By Mechanical Methods (AREA)
- Preventing Corrosion Or Incrustation Of Metals (AREA)
Abstract
本发明公开了一种大气腐蚀等级图绘制方法,所述方法首先通过现场暴露试验收集大气腐蚀速率数据,再收集大气腐蚀等级图绘图区域内各个气象站点的气象环境数据,然后通过收集的大气腐蚀速率数据和气象环境数据运用回归分析方法建立剂量响应函数模型进行大气腐蚀等级图的绘制操作,最后进行电力输电系统的数据植入。利用所述方法绘制的大气腐蚀等级图能用于南方电网地区大气腐蚀性的评估及材料的服役寿命预测,特别增加电力输电系统的地理信息数据,能够为输电线路的布设、检修、防腐蚀等工作提供便利指导。
Description
技术领域
本发明涉及大气腐蚀等级图绘制领域,具体涉及一种区域大气腐蚀等级图的计算与绘制方法。
背景技术
电力行业是支撑国民经济和社会发展的基础性产业和公用事业,随着我国国民经济的快速发展和人民生活水平的不断提高,对电力的依赖程度也越来越高。输电杆塔是电力传输的桥梁,而长期暴晒于自然环境中的输电杆塔钢构件的防腐蚀性能,是保证输电线路安全可靠运行的关键要素之一,在沿海和因重工业生产排放所带来的环境污染并伴随自然环境气候的恶化而形成的重腐蚀环境中,在役杆塔钢构件等重要设备的腐蚀日渐加快,因此,针对不同腐蚀等级的区域,在生产建设过程中,应该选取耐腐蚀不同程度的防护材料和维护技术,才能实现钢构件长寿命、少维护或免维护的防腐效果,从而满足输电线路安全运行的需要。
目前,腐蚀实验大多关注某一大气环境下材质的腐蚀行为,都是对实验时所测数据进行分析,而从宏观上覆盖整个区域的大气腐蚀行为及腐蚀等级划分的研究与应用整体很缺乏。在现有腐蚀实验技术条件下,借助地学模型的构建,进行区域大气腐蚀速率计算和腐蚀等级划分,能够为输电线路的布设、检修、防腐蚀等工作提供科学的指导,具有非常现实的意义。
发明内容
本发明的目的是针对现有技术的不足,为弥补现有南方电网地区电力行业中大气腐蚀等级图绘制方法的空白,提供了一种大气腐蚀等级图绘制方法,所述方法基于挂片点实测腐蚀数据、气象数据、污源点数据及海岸线影响数据,建立了大气腐蚀的剂量响应函数模型,从而完成大气腐蚀等级图的绘制操作,同时通过电力输电系统数据的植入,生成了适用于电力系统的大气腐蚀等级图。
本发明的目的可以通过如下技术方案实现:
一种大气腐蚀等级图绘制方法,所述方法包括以下步骤:
S1、大气腐蚀速率数据收集:通过现场暴露试验,包括输电杆塔现场挂片、挂片回收和清洗以及挂片腐蚀速率计算,从而获得各个现场暴露试验点的实测大气腐蚀速率数据;
S2、气象环境数据收集:收集大气腐蚀等级图绘图区域各个气象站点的气象环境数据,收集对象包括各个气象站点的坐标、温度、相对湿度以及二氧化硫、二氧化氮和臭氧的浓度,并计算温度、相对湿度以及二氧化硫、二氧化氮和臭氧浓度的年平均数值;
S3、剂量响应函数模型建立:运用回归分析方法,建立大气腐蚀速率数据和气象环境数据间的函数模型1;建立污染源影响大气腐蚀速率的函数模型2;建立海洋影响大气腐蚀速率的函数模型3;
S4、大气腐蚀等级图绘制操作:根据步骤S3的三个函数模型和GB/T19292.1-2003标准中的大气腐蚀等级划分,运用ArcGIS软件中的统计分析模块、地图叠加分析模块、重分类模块绘制大气腐蚀等级图;
S5、电力输电系统数据植入:对输电线路杆塔、变电站、电厂和行政分区的地理信息数据进行矢量化处理,将所有进行矢量化处理后的数据的坐标系统一,并进行图形符号的配置,生成适用于电力输电系统的大气腐蚀等级图。
进一步地,步骤S1中,所述现场暴露试验是根据GB/T16145-2008《金属和合金大气腐蚀试验-现场试验的一般要求》中的现场挂片标准进行现场挂片,挂片试样的尺寸为100mm×50mm,厚度为3mm,试样表面粗糙度为2μm,试样卡槽为橡胶材料,以隔离铝合金架子与金属试片;同时,试样架平面与水平面成30°角,现场暴露架被安装在离地面高5~6m处,方向朝南,一个暴露架上4个试片,历时一年回收暴露架;
挂片回收和清洗是根据GB/T16145《金属和合金的腐蚀-腐蚀试样上腐蚀产物的清除》中的清洗标准,首先用流水冲刷掉挂片表面疏松的锈层,然后按照化学配比500mL浓度1.19g/mL的盐酸+3.5g六次甲基四胺+蒸馏水配置成1000mL的清洗溶液,将经流水冲刷过的挂片放入装有清洗溶液的容器内,在温度20℃~25℃下,人工清洗10min或超声振荡不超过5min,将腐蚀产物去除后,再用水和酒精先后清洗挂片,充分干燥后用电子天平称重,精确到±0.1mg;
通过腐蚀前和腐蚀后的挂片的质量差,计算大气腐蚀速率,计算公式为:
其中,ML表示大气腐蚀速率,单位为g·m-2·y-1;W0表示腐蚀前的挂片的质量,单位为g,WT表示腐蚀后经过清洗后的挂片的质量,单位为g,S表示挂片的面积,单位为m2,T表示现场暴露试验的时间,单位为y。
进一步地,步骤S2中,所述温度、相对湿度以及二氧化硫、二氧化氮和臭氧浓度的年平均数值计算公式为:
其中,xij表示温度、相对湿度或者二氧化硫、二氧化氮、臭氧浓度在第i天第j时的数值,表示温度、相对湿度或者二氧化硫、二氧化氮、臭氧浓度的年平均数值。
进一步地,步骤S3中,所述大气腐蚀速率数据和气象环境数据间的函数模型1的公式为:
其中,ML表示大气腐蚀速率,单位为g·m-2·y-1;T表示温度,单位为℃;Rh表示相对湿度;SO2表示二氧化硫的浓度,单位为μg/m3;O3表示臭氧的浓度,单位为μg/m3;NO2表示二氧化氮的浓度,单位为μg/m3;
所述污染源影响大气腐蚀速率的函数模型2的公式为:
ML=-A·lnDw+B
其中,ML表示大气腐蚀速率,单位为g·m-2·y-1;A、B为污染源腐蚀系数;Dw表示距离污染源的距离,单位为km;当污染源为重度污染时,A=63.78,B=99.063;当污染源为中度污染时,A=32.64,B=55.71;当污染源为轻度污染时,Dw的取值范围为0~1.5,A=0,B=25;
所述海洋影响大气腐蚀速率的函数模型3的公式为:
ML=A·Ds -B
其中,ML表示大气腐蚀速率,单位为g·m-2·y-1;A、B为污染源腐蚀系数;Ds表示距离海岸线的距离,单位为m;污染源腐蚀系数A=1.119×106,B=0.984。
进一步地,所述步骤S4的具体过程为:
步骤1、将收集的大气腐蚀等级图绘图区域各个气象站点的原始气象环境数据矢量化,将原始气象环境数据根据函数模型1计算出大气腐蚀速率,并通过气象站点坐标生成矢量点文件,矢量点文件坐标系为WGS84;
步骤2、通过ArcGIS软件中的统计分析模块,采用克里金插值法对步骤1中根据函数模型1计算出的大气腐蚀速率数据结合通过现场暴露试验得到的大气腐蚀速率数据进行空间插值,得到大气腐蚀等级图的底图;
步骤3、对污染源的坐标根据函数模型2以反距离权重插值法进行空间插值,以50个大气腐蚀变化量为间距,根据不同程度的污染源进行缓冲生成矢量面文件,并转为矢量线文件,再由矢量线文件生成矢量点文件,根据不同距离的污染源分别赋值,运用反距离权重插值法对矢量点文件进行空间插值;
步骤4、运用地图叠加分析模块,将步骤3中的插值结果与步骤2中的插值结果数据进行空间叠加运算,再依据海洋影响大气腐蚀速率的函数模型3的递减规律,采用每50个腐蚀速率递减量的距离为一个缓冲区,得到各缓冲区的大气腐蚀等级值,并以多环区的腐蚀值进行空间插值计算,生成连续的海洋影响下大气腐蚀速率数据,并根据GB/T19292.1-2003《金属和合金的腐蚀大气腐蚀性分类》,将大气腐蚀从低至高划分为6个等级,分别为:C1-干净、C2-轻微、C3-中度、C4-严重、C5-非常严重、Cx-极其严重,根据上述分级标准进行腐蚀等级的划分,从而生成最终的大气腐蚀等级图。
本发明与现有技术相比,具有如下优点和有益效果:
1、本发明通过对不同大气环境下的腐蚀数据进行分析处理,构建了基于气象数据与污染源等数据的大气腐蚀剂量响应函数模型,在充分考虑污染源和海洋对大气腐蚀影响的基础上,绘制了大气腐蚀等级图,所述大气腐蚀等级图既能根据不同地区的腐蚀等级程度因地制宜地设计指导输电线路的设计,也能在输电线路进行建设或改造时,指导建设方根据不同的大气环境腐蚀等级而选择不同的耐腐蚀材料,从而有效地节约成本;
2、本发明的方法解决了从传统的单一环境大气腐蚀数据向区域大气腐蚀数据转化的途径,所绘制出的大气腐蚀等级图能用于地区大气腐蚀性的评估及材料的服役寿命预测,特别增加电力输电系统的地理信息数据,能够为输电线路的布设、检修、防腐蚀等工作提供便利指导。
附图说明
图1为本发明实施例现场暴露试验的流程图。
图2为本发明实施例一种大气腐蚀等级图绘制方法的流程图。
具体实施方式
下面结合实施例及附图对本发明作进一步详细的描述,但本发明的实施方式不限于此。
实施例:
本实施例提供了一种大气腐蚀等级图绘制方法,所述方法的流程图如图2所示,包括以下步骤:
S1、大气腐蚀速率数据收集:通过现场暴露试验,流程图如图1所示,包括输电杆塔现场挂片、挂片回收和清洗以及挂片腐蚀速率计算,从而获得各个现场暴露试验点的实测大气腐蚀速率数据;
所述现场暴露试验是根据GB/T16145-2008《金属和合金大气腐蚀试验-现场试验的一般要求》中的现场挂片标准进行现场挂片,挂片试样的尺寸为100mm×50mm,厚度为3mm,试样表面粗糙度为2μm,试样卡槽为橡胶材料,以隔离铝合金架子与金属试片;同时,试样架平面与水平面成30°角,现场暴露架被安装在离地面高5~6m处,方向朝南,一个暴露架上4个试片,历时一年回收暴露架;
挂片回收和清洗是根据GB/T16145《金属和合金的腐蚀-腐蚀试样上腐蚀产物的清除》中的清洗标准,首先用流水冲刷掉挂片表面疏松的锈层,然后按照化学配比500mL浓度1.19g/mL的盐酸+3.5g六次甲基四胺+蒸馏水配置成1000mL的清洗溶液,将经流水冲刷过的挂片放入装有清洗溶液的容器内,在温度20℃~25℃下,人工清洗10min或超声振荡不超过5min,将腐蚀产物去除后,再用水和酒精先后清洗挂片,充分干燥后用电子天平称重,精确到±0.1mg;
通过腐蚀前和腐蚀后的挂片的质量差,计算大气腐蚀速率,计算公式为:
其中,ML表示大气腐蚀速率,单位为g·m-2·y-1;W0表示腐蚀前的挂片的质量,单位为g,WT表示腐蚀后经过清洗后的挂片的质量,单位为g,S表示挂片的面积,单位为m2,T表示现场暴露试验的时间,单位为y。
S2、气象环境数据收集:收集大气腐蚀等级图绘图区域各个气象站点的气象环境数据,收集对象包括各个气象站点的坐标、温度、相对湿度以及二氧化硫、二氧化氮和臭氧的浓度,并计算温度、相对湿度以及二氧化硫、二氧化氮和臭氧浓度的年平均数值;
所述温度、相对湿度以及二氧化硫、二氧化氮和臭氧浓度的年平均数值计算公式为:
其中,xij表示温度、相对湿度或者二氧化硫、二氧化氮、臭氧浓度在第i天第j时的数值,x表示温度、相对湿度或者二氧化硫、二氧化氮、臭氧浓度的年平均数值。
S3、剂量响应函数模型建立:运用回归分析方法,建立大气腐蚀速率数据和气象环境数据间的函数模型1;建立污染源影响大气腐蚀速率的函数模型2;建立海洋影响大气腐蚀速率的函数模型3;
材料的腐蚀和退化是由多种空气污染物的化学反应和气象因子共同作用导致的,空气污染物主要包括二氧化硫、二氧化氮以及臭氧,气象因子主要有温度、相对湿度等,都会影响材料的腐蚀及退化程度。所述大气腐蚀速率数据和气象环境数据间的函数模型1的公式为:
ML=8.484×10-6·T0.485·Rh2.719·SO2 0.413·O3 0.460·NO2 0.261
其中,ML表示大气腐蚀速率,单位为g·m-2·y-1;T表示温度,单位为℃;Rh表示相对湿度;SO2表示二氧化硫的浓度,单位为μg/m3;O3表示臭氧的浓度,单位为μg/m3;NO2表示二氧化氮的浓度,单位为μg/m3;
由于污染源的污染会导致周围大气腐蚀加剧,需要对污染源周围区域大气腐蚀速率进行处理,根据污染源的污染程度,分为重度污染、中度污染、轻度污染三个等级,根据在污源点附近的挂片所测得的大气腐蚀数据,构建污染源对大气腐蚀影响的修正模型,具体步骤:
1)收集绘制区域内工业污染源的GIS坐标;
2)在GIS系统中统计并记录暴露试验点周围是否存在工业污染源(距离≤4km),并确定污染源离暴露试验点的距离;
3)用暴露试验同一线路上受工业污染影响的暴露试验点的腐蚀量减去同一线路上无工业污染影响的暴露试验点中的最小腐蚀量,即为对应的工业污染源所叠加的腐蚀量;
4)由于不同污染源的污染程度不同,根据3)所得数据将污染源等级分为重、中、轻三级;
5)对污染源等级分类后的数据分别进行函数拟合,得到污染源的距离与叠加腐蚀量的函数,用于大气腐蚀图的绘制。
所述污染源影响大气腐蚀速率的函数模型2的公式为:
ML=-A·lnDw+B
其中,ML表示大气腐蚀速率,单位为g·m-2·y-1;A、B为污染源腐蚀系数;Dw表示距离污染源的距离,单位为km;当污染源为重度污染时,A=63.78,B=99.063;当污染源为中度污染时,A=32.64,B=55.71;当污染源为轻度污染时,Dw的取值范围为0~1.5,A=0,B=25;
由于海洋的盐雾作用会导致周围大气腐蚀加剧,需要对海岸线周围区域的大气腐蚀速率进行处理,根据现场暴露试验挂片点距海岸线的距离及挂片点的腐蚀数据,构建海洋对大气腐蚀影响的修正模型,具体步骤:
1)在GIS系统中测得靠近海洋的挂片点离海岸线的距离(≤10000m),并记录下来;
2)用暴露试验同一线路上受海洋影响的暴露试验点的腐蚀量减去同一线路上无海洋影响的暴露试验点中的最小腐蚀量,即为对应海洋影响大气腐蚀所叠加的腐蚀量;
3)对所得数据进行函数拟合,得到距海岸线的距离与叠加腐蚀量的函数,用于大气腐蚀图的绘制。
所述海洋影响大气腐蚀速率的函数模型3的公式为:
ML=A·Ds -B
其中,ML表示大气腐蚀速率,单位为g·m-2·y-1;A、B为污染源腐蚀系数;Ds表示距离海岸线的距离,单位为m;污染源腐蚀系数A=1.119×106,B=0.984。
S4、大气腐蚀等级图绘制操作:根据步骤S3的三个函数模型和GB/T19292.1-2003标准中的大气腐蚀等级划分,运用ArcGIS软件中的统计分析模块、地图叠加分析模块、重分类模块绘制大气腐蚀等级图;具体过程为:
步骤1、将收集的大气腐蚀等级图绘图区域各个气象站点的原始气象环境数据矢量化,将原始气象环境数据根据函数模型1计算出大气腐蚀速率,并通过气象站点坐标生成矢量点文件,矢量点文件坐标系为WGS84;
步骤2、通过ArcGIS软件中的统计分析模块,采用克里金插值法对步骤1中根据函数模型1计算出的大气腐蚀速率数据结合通过现场暴露试验得到的大气腐蚀速率数据进行空间插值,得到大气腐蚀等级图的底图;
步骤3、对污染源的坐标根据函数模型2以反距离权重插值法进行空间插值,以50个大气腐蚀变化量为间距,根据不同程度的污染源进行缓冲生成矢量面文件,并转为矢量线文件,再由矢量线文件生成矢量点文件,根据不同距离的污染源分别赋值,运用反距离权重插值法对矢量点文件进行空间插值;
步骤4、运用地图叠加分析模块,将步骤3中的插值结果与步骤2中的插值结果数据进行空间叠加运算,再依据海洋影响大气腐蚀速率的函数模型3的递减规律,采用每50个腐蚀速率递减量的距离为一个缓冲区,得到各缓冲区的大气腐蚀等级值,并以多环区的腐蚀值进行空间插值计算,生成连续的海洋影响下大气腐蚀速率数据,并根据GB/T19292.1-2003《金属和合金的腐蚀大气腐蚀性分类》,将大气腐蚀从低至高划分为6个等级,分别为:C1-干净、C2-轻微、C3-中度、C4-严重、C5-非常严重、Cx-极其严重,根据上述分级标准进行腐蚀等级的划分,从而生成最终的大气腐蚀等级图。
S5、电力输电系统数据植入:对输电线路杆塔、变电站、电厂和行政分区的地理信息数据进行矢量化处理,将所有进行矢量化处理后的数据的坐标系统一,并进行图形符号的配置,生成适用于电力输电系统的大气腐蚀等级图。
以上所述,仅为本发明专利较佳的实施例,但本发明专利的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明专利所公开的范围内,根据本发明专利的技术方案及其发明专利构思加以等同替换或改变,都属于本发明专利的保护范围。
Claims (5)
1.一种大气腐蚀等级图绘制方法,其特征在于,所述方法包括以下步骤:
S1、大气腐蚀速率数据收集:通过现场暴露试验,包括输电杆塔现场挂片、挂片回收和清洗以及挂片腐蚀速率计算,从而获得各个现场暴露试验点的实测大气腐蚀速率数据;
S2、气象环境数据收集:收集大气腐蚀等级图绘图区域各个气象站点的气象环境数据,收集对象包括各个气象站点的坐标、温度、相对湿度以及二氧化硫、二氧化氮和臭氧的浓度,并计算温度、相对湿度以及二氧化硫、二氧化氮和臭氧浓度的年平均数值;
S3、剂量响应函数模型建立:运用回归分析方法,建立大气腐蚀速率数据和气象环境数据间的函数模型1;建立污染源影响大气腐蚀速率的函数模型2;建立海洋影响大气腐蚀速率的函数模型3;
S4、大气腐蚀等级图绘制操作:根据步骤S3的三个函数模型和GB/T19292.1-2003标准中的大气腐蚀等级划分,运用ArcGIS软件中的统计分析模块、地图叠加分析模块、重分类模块绘制大气腐蚀等级图;
S5、电力输电系统数据植入:对输电线路杆塔、变电站、电厂和行政分区的地理信息数据进行矢量化处理,将所有进行矢量化处理后的数据的坐标系统一,并进行图形符号的配置,生成适用于电力输电系统的大气腐蚀等级图。
2.根据权利要求1所述的一种大气腐蚀等级图绘制方法,其特征在于:步骤S1中,所述现场暴露试验是根据GB/T16145-2008《金属和合金大气腐蚀试验-现场试验的一般要求》中的现场挂片标准进行现场挂片,挂片试样的尺寸为100mm×50mm,厚度为3mm,试样表面粗糙度为2μm,试样卡槽为橡胶材料,以隔离铝合金架子与金属试片;同时,试样架平面与水平面成30°角,现场暴露架被安装在离地面高5~6m处,方向朝南,一个暴露架上4个试片,历时一年回收暴露架;
挂片回收和清洗是根据GB/T16145《金属和合金的腐蚀-腐蚀试样上腐蚀产物的清除》中的清洗标准,首先用流水冲刷掉挂片表面疏松的锈层,然后按照化学配比500mL浓度1.19g/mL的盐酸+3.5g六次甲基四胺+蒸馏水配置成1000mL的清洗溶液,将经流水冲刷过的挂片放入装有清洗溶液的容器内,在温度20℃~25℃下,人工清洗10min或超声振荡不超过5min,将腐蚀产物去除后,再用水和酒精先后清洗挂片,充分干燥后用电子天平称重,精确到±0.1mg;
通过腐蚀前和腐蚀后的挂片的质量差,计算大气腐蚀速率,计算公式为:
<mrow>
<mi>M</mi>
<mi>L</mi>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>W</mi>
<mn>0</mn>
</msub>
<mo>-</mo>
<msub>
<mi>W</mi>
<mi>T</mi>
</msub>
</mrow>
<mrow>
<mi>S</mi>
<mo>&CenterDot;</mo>
<mi>T</mi>
</mrow>
</mfrac>
</mrow>
其中,ML表示大气腐蚀速率,单位为g·m-2·y-1;W0表示腐蚀前的挂片的质量,单位为g,WT表示腐蚀后经过清洗后的挂片的质量,单位为g,S表示挂片的面积,单位为m2,T表示现场暴露试验的时间,单位为y。
3.根据权利要求1所述的一种大气腐蚀等级图绘制方法,其特征在于,步骤S2中,所述温度、相对湿度以及二氧化硫、二氧化氮和臭氧浓度的年平均数值计算公式为:
<mrow>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>365</mn>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mn>365</mn>
</munderover>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>/</mo>
<mn>24</mn>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mn>24</mn>
</munderover>
<msub>
<mi>x</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
其中,xij表示温度、相对湿度或者二氧化硫、二氧化氮、臭氧浓度在第i天第j时的数值,表示温度、相对湿度或者二氧化硫、二氧化氮、臭氧浓度的年平均数值。
4.根据权利要求1所述的一种大气腐蚀等级图绘制方法,其特征在于:步骤S3中,所述大气腐蚀速率数据和气象环境数据间的函数模型1的公式为:
ML=8.484×10-6·T0.485·Rh2.719·SO2 0.413·O3 0.460·NO2 0.261
其中,ML表示大气腐蚀速率,单位为g·m-2·y-1;T表示温度,单位为℃;Rh表示相对湿度;SO2表示二氧化硫的浓度,单位为μg/m3;O3表示臭氧的浓度,单位为μg/m3;NO2表示二氧化氮的浓度,单位为μg/m3;
所述污染源影响大气腐蚀速率的函数模型2的公式为:
ML=-A·lnDw+B
其中,ML表示大气腐蚀速率,单位为g·m-2·y-1;A、B为污染源腐蚀系数;Dw表示距离污染源的距离,单位为km;当污染源为重度污染时,A=63.78,B=99.063;当污染源为中度污染时,A=32.64,B=55.71;当污染源为轻度污染时,Dw的取值范围为0~1.5,A=0,B=25;
所述海洋影响大气腐蚀速率的函数模型3的公式为:
ML=A·Ds -B
其中,ML表示大气腐蚀速率,单位为g·m-2·y-1;A、B为污染源腐蚀系数;Ds表示距离海岸线的距离,单位为m;污染源腐蚀系数A=1.119×106,B=0.984。
5.根据权利要求4所述的一种大气腐蚀等级图绘制方法,其特征在于,所述步骤S4的具体过程为:
步骤1、将收集的大气腐蚀等级图绘图区域各个气象站点的原始气象环境数据矢量化,将原始气象环境数据根据函数模型1计算出大气腐蚀速率,并通过气象站点坐标生成矢量点文件,矢量点文件坐标系为WGS84;
步骤2、通过ArcGIS软件中的统计分析模块,采用克里金插值法对步骤1中根据函数模型1计算出的大气腐蚀速率数据结合通过现场暴露试验得到的大气腐蚀速率数据进行空间插值,得到大气腐蚀等级图的底图;
步骤3、对污染源的坐标根据函数模型2以反距离权重插值法进行空间插值,以50个大气腐蚀变化量为间距,根据不同程度的污染源进行缓冲生成矢量面文件,并转为矢量线文件,再由矢量线文件生成矢量点文件,根据不同距离的污染源分别赋值,运用反距离权重插值法对矢量点文件进行空间插值;
步骤4、运用地图叠加分析模块,将步骤3中的插值结果与步骤2中的插值结果数据进行空间叠加运算,再依据海洋影响大气腐蚀速率的函数模型3的递减规律,采用每50个腐蚀速率递减量的距离为一个缓冲区,得到各缓冲区的大气腐蚀等级值,并以多环区的腐蚀值进行空间插值计算,生成连续的海洋影响下大气腐蚀速率数据,并根据GB/T19292.1-2003《金属和合金的腐蚀大气腐蚀性分类》,将大气腐蚀从低至高划分为6个等级,分别为:C1-干净、C2-轻微、C3-中度、C4-严重、C5-非常严重、Cx-极其严重,根据上述分级标准进行腐蚀等级的划分,从而生成最终的大气腐蚀等级图。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711131439.8A CN108020501B (zh) | 2017-11-15 | 2017-11-15 | 一种大气腐蚀等级图绘制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711131439.8A CN108020501B (zh) | 2017-11-15 | 2017-11-15 | 一种大气腐蚀等级图绘制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108020501A true CN108020501A (zh) | 2018-05-11 |
CN108020501B CN108020501B (zh) | 2019-08-20 |
Family
ID=62080825
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711131439.8A Expired - Fee Related CN108020501B (zh) | 2017-11-15 | 2017-11-15 | 一种大气腐蚀等级图绘制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108020501B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109187380A (zh) * | 2018-07-25 | 2019-01-11 | 北京农业智能装备技术研究中心 | 一种水质探测器、水体污染排放源监测系统及方法 |
CN109580464A (zh) * | 2018-11-22 | 2019-04-05 | 广西电网有限责任公司电力科学研究院 | 一种检测评价电网设备涂层质量的方法 |
CN112033879A (zh) * | 2020-07-16 | 2020-12-04 | 国网山东省电力公司电力科学研究院 | 一种大气腐蚀性数据插值方法及系统 |
US11009448B2 (en) * | 2019-10-11 | 2021-05-18 | Matergenics, Inc. | Atmospheric corrosivity mapping method and apparatus |
CN113408131A (zh) * | 2021-06-23 | 2021-09-17 | 武汉钢铁有限公司 | 耐候桥梁钢腐蚀寿命预测方法 |
CN113781307A (zh) * | 2021-09-13 | 2021-12-10 | 国网山东省电力公司电力科学研究院 | 大气腐蚀多因素空间插值算法、计算机设备、存储介质 |
CN113970515A (zh) * | 2020-07-24 | 2022-01-25 | 宝山钢铁股份有限公司 | 一种准确测量金属材料腐蚀增重的方法及应用 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2008157647A (ja) * | 2006-12-21 | 2008-07-10 | Tokyo Electric Power Co Inc:The | Acmセンサによる構造物の腐食速度推定方法 |
CN105095589A (zh) * | 2015-08-10 | 2015-11-25 | 贵州电网有限责任公司电力科学研究院 | 一种山区电网风区分布图绘制方法 |
CN105606869A (zh) * | 2015-09-10 | 2016-05-25 | 国家电网公司 | 一种输电线路大气环境腐蚀严酷性快速评估方法 |
CN105956389A (zh) * | 2016-04-26 | 2016-09-21 | 南方电网科学研究院有限责任公司 | 环境因素数据的采集方法及金属大气腐蚀等级图绘制方法 |
CN105973788A (zh) * | 2016-04-25 | 2016-09-28 | 南方电网科学研究院有限责任公司 | 一种大气腐蚀挂片的耐腐蚀速率测定方法 |
CN106295653A (zh) * | 2016-07-29 | 2017-01-04 | 宁波大学 | 一种水质图像分类方法 |
CN106483060A (zh) * | 2016-09-22 | 2017-03-08 | 国网山东省电力公司电力科学研究院 | 一种动态大气腐蚀区域地图数据处理方法、装置及系统 |
-
2017
- 2017-11-15 CN CN201711131439.8A patent/CN108020501B/zh not_active Expired - Fee Related
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2008157647A (ja) * | 2006-12-21 | 2008-07-10 | Tokyo Electric Power Co Inc:The | Acmセンサによる構造物の腐食速度推定方法 |
CN105095589A (zh) * | 2015-08-10 | 2015-11-25 | 贵州电网有限责任公司电力科学研究院 | 一种山区电网风区分布图绘制方法 |
CN105606869A (zh) * | 2015-09-10 | 2016-05-25 | 国家电网公司 | 一种输电线路大气环境腐蚀严酷性快速评估方法 |
CN105973788A (zh) * | 2016-04-25 | 2016-09-28 | 南方电网科学研究院有限责任公司 | 一种大气腐蚀挂片的耐腐蚀速率测定方法 |
CN105956389A (zh) * | 2016-04-26 | 2016-09-21 | 南方电网科学研究院有限责任公司 | 环境因素数据的采集方法及金属大气腐蚀等级图绘制方法 |
CN106295653A (zh) * | 2016-07-29 | 2017-01-04 | 宁波大学 | 一种水质图像分类方法 |
CN106483060A (zh) * | 2016-09-22 | 2017-03-08 | 国网山东省电力公司电力科学研究院 | 一种动态大气腐蚀区域地图数据处理方法、装置及系统 |
Non-Patent Citations (6)
Title |
---|
I.S.COLE ET AL.: "An Australia0wide map of corrosivity: a GIS approach", 《DURABILITY OF BUILDING AN MATERIALS AND COMPONENTS 8》 * |
K.SLAMOVA ET AL.: "Measurement and GIS-base spatial modeling of copper corrosion in different environments in europe", 《MATERIALS AND CORROSION 》 * |
REISS D.ET AL.: "Mapping stock at risk and release of zinc and copper in swizerland-dose response functions for runoff rates derived from corrosion rate data", 《WATER AIR AND SOIL POLLUTION》 * |
吴超等: "基于数据分析的大气腐蚀等级细化研究", 《装备环境工程》 * |
唐其环等: "国内外典型大气腐蚀图及其绘制", 《装备环境工程》 * |
蒋武斌: "南网区域输电杆塔大气腐蚀等级评估及腐蚀行为研究_海南省和云南省", 《中国优秀硕士学位论文全文数据库》 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109187380A (zh) * | 2018-07-25 | 2019-01-11 | 北京农业智能装备技术研究中心 | 一种水质探测器、水体污染排放源监测系统及方法 |
CN109580464A (zh) * | 2018-11-22 | 2019-04-05 | 广西电网有限责任公司电力科学研究院 | 一种检测评价电网设备涂层质量的方法 |
US11009448B2 (en) * | 2019-10-11 | 2021-05-18 | Matergenics, Inc. | Atmospheric corrosivity mapping method and apparatus |
CN112033879A (zh) * | 2020-07-16 | 2020-12-04 | 国网山东省电力公司电力科学研究院 | 一种大气腐蚀性数据插值方法及系统 |
CN113970515A (zh) * | 2020-07-24 | 2022-01-25 | 宝山钢铁股份有限公司 | 一种准确测量金属材料腐蚀增重的方法及应用 |
CN113970515B (zh) * | 2020-07-24 | 2023-11-14 | 宝山钢铁股份有限公司 | 一种准确测量金属材料腐蚀增重的方法及应用 |
CN113408131A (zh) * | 2021-06-23 | 2021-09-17 | 武汉钢铁有限公司 | 耐候桥梁钢腐蚀寿命预测方法 |
CN113781307A (zh) * | 2021-09-13 | 2021-12-10 | 国网山东省电力公司电力科学研究院 | 大气腐蚀多因素空间插值算法、计算机设备、存储介质 |
CN113781307B (zh) * | 2021-09-13 | 2023-09-29 | 国网山东省电力公司电力科学研究院 | 大气腐蚀多因素空间插值算法、计算机设备、存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN108020501B (zh) | 2019-08-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108020501B (zh) | 一种大气腐蚀等级图绘制方法 | |
CN112905560B (zh) | 一种多源时空大数据深度融合的空气污染预测方法 | |
CN104809511B (zh) | 绝缘子污染预测方法及设备 | |
Paxian et al. | Present-day and future global bottom-up ship emission inventories including polar routes | |
Slamova et al. | Mapping atmospheric corrosion in coastal regions: methods and results | |
CN104008278B (zh) | 一种基于特征向量和最小二乘支持向量机的pm2.5浓度预测方法 | |
CN105631537A (zh) | 基于气象服务平台的空气质量预报业务系统 | |
CN101059418B (zh) | 斜拉桥拉索腐蚀状况的评价方法 | |
CN110646867A (zh) | 一种城市排水监测预警方法和系统 | |
CN113436045A (zh) | 大气污染多尺度融合溯源方法、装置和可读存储介质 | |
CN109034252A (zh) | 空气质量站点监测数据异常的自动化识别方法 | |
CN111428942A (zh) | 基于变网格技术提取微地形因子的线路覆冰厚度预测方法 | |
CN111103405A (zh) | 空气质量预报预警系统 | |
Wong et al. | Estimation of aerosol sources and aerosol transport pathways using AERONET clustering and backward trajectories: a case study of Hong Kong | |
CN110889196A (zh) | 水环境承载力评估方法、装置及存储介质 | |
Otache et al. | ARMA modelling of Benue River flow dynamics: comparative study of PAR model | |
Abdolmanafi et al. | Evaluation of global ensemble prediction models for forecasting medium to heavy precipitations | |
CN116703004B (zh) | 一种基于预训练模型的水系流域智慧巡护方法和装置 | |
Obata et al. | Numerical simulation of adhesion of sea-salt particles on bridge girders | |
CN106056477A (zh) | 一种基于用电大数据的行业产能利用率的计算方法 | |
CN109697630A (zh) | 一种基于稀疏回归的商家客流量多因素分析与预测方法 | |
CN113780880A (zh) | 一种基于空间代表性的pm2.5观测站点布局评价方法 | |
CN102128914A (zh) | 一种水质量预测、评价的等浓度线方法 | |
Man et al. | Investigating metrological and geographical effect in remote sensing retrival of PM2. 5 concentration in Yangtze River Delta | |
Nykolaychuk et al. | Structure and functioning of information systems of background monitoring of landscape elements of gorgany nature reserve |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190820 Termination date: 20201115 |
|
CF01 | Termination of patent right due to non-payment of annual fee |