JP7084344B2 - Corrosion rate estimation device, corrosion rate estimation method and program - Google Patents

Corrosion rate estimation device, corrosion rate estimation method and program Download PDF

Info

Publication number
JP7084344B2
JP7084344B2 JP2019044872A JP2019044872A JP7084344B2 JP 7084344 B2 JP7084344 B2 JP 7084344B2 JP 2019044872 A JP2019044872 A JP 2019044872A JP 2019044872 A JP2019044872 A JP 2019044872A JP 7084344 B2 JP7084344 B2 JP 7084344B2
Authority
JP
Japan
Prior art keywords
corrosion rate
corrosion
dry
metal material
estimation
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
JP2019044872A
Other languages
Japanese (ja)
Other versions
JP2020148541A (en
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.)
Hokkaido Electric Power Co Inc
Original Assignee
Hokkaido Electric Power Co Inc
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 Hokkaido Electric Power Co Inc filed Critical Hokkaido Electric Power Co Inc
Priority to JP2019044872A priority Critical patent/JP7084344B2/en
Publication of JP2020148541A publication Critical patent/JP2020148541A/en
Application granted granted Critical
Publication of JP7084344B2 publication Critical patent/JP7084344B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y04INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04SSYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
    • Y04S10/00Systems supporting electrical power generation, transmission or distribution
    • Y04S10/50Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications

Landscapes

  • Testing Resistance To Weather, Investigating Materials By Mechanical Methods (AREA)

Description

本発明は、腐食速度推定装置、腐食速度推定方法及びプログラムに関する。 The present invention relates to a corrosion rate estimation device, a corrosion rate estimation method and a program.

鉄塔に用いられる亜鉛メッキ鋼材や送電線に用いられるアルミ系電線は、異種金属同士が接触する内部で腐食が発生し、目視による腐食点検が困難である。また、鉄塔や送電線は施設数が多く、かつ広域に設置されているため、老朽化状況を考慮し合理的な改修計画を策定する必要がある。このため、鉄塔や送電線への腐食の影響を高精度に評価する手法の研究が進められている。例えば、特許文献1には、亜鉛系めっき機材の腐食速度と銅片の光沢度又は重量増加量との相関図を作成することで、亜鉛系めっき機材の余寿命を評価する方法が開示されている。 Galvanic steel used for steel towers and aluminum-based electric wires used for power transmission lines are corroded inside where dissimilar metals come into contact with each other, making visual corrosion inspection difficult. In addition, since steel towers and transmission lines have a large number of facilities and are installed over a wide area, it is necessary to formulate a rational renovation plan in consideration of the aging situation. For this reason, research is underway on methods for highly accurate evaluation of the effects of corrosion on steel towers and transmission lines. For example, Patent Document 1 discloses a method for evaluating the remaining life of zinc-based plating equipment by creating a correlation diagram between the corrosion rate of zinc-based plating equipment and the glossiness or weight increase of copper pieces. There is.

特開2017-106785号公報JP-A-2017-106785

特許文献1の手法では、鉄塔や送電線の設置箇所によって異なる周囲環境が鉄塔や送電線の腐食に与える影響を考慮しておらず、推定される腐食速度の精度の点で改善の余地がある。そして、このような問題は、鉄塔や送電線の場合に限られず、例えば、屋外の腐食環境等に設置されたその他の金属材料の腐食速度を推定する場合にも存在している。 The method of Patent Document 1 does not consider the influence of the surrounding environment, which differs depending on the installation location of the tower or transmission line, on the corrosion of the tower or transmission line, and there is room for improvement in terms of the accuracy of the estimated corrosion rate. .. And, such a problem is not limited to the case of a steel tower or a transmission line, and also exists in the case of estimating the corrosion rate of other metal materials installed in an outdoor corrosive environment or the like, for example.

本発明は、このような背景に基づいてなされたものであり、周囲環境の影響を考慮して金属材料の腐食速度を定量化することが可能な腐食速度推定装置、腐食速度推定方法及びプログラムを提供することを目的とする。 The present invention has been made based on such a background, and provides a corrosion rate estimation device, a corrosion rate estimation method, and a program capable of quantifying the corrosion rate of a metal material in consideration of the influence of the surrounding environment. The purpose is to provide.

上記目的を達成するために、本発明の第1の観点に係る腐食速度推定装置は、
金属材料の腐食速度を目的変数とし、前記金属材料の腐食に影響を与える年間飛来海塩粒子量、濡れ時間及び乾湿サイクルを含む腐食環境因子を説明変数とする重回帰分析を行うことで、前記金属材料の腐食速度の推定式を作成する推定式作成手段と、
前記推定式に年間飛来海塩粒子量、濡れ時間及び乾湿サイクルを含む腐食環境因子を入力することで、前記金属材料の腐食速度を推定する腐食速度推定手段と、
を備え、
乾湿サイクルは、所定時間以上の湿潤状態と所定時間以上の乾燥状態とが交互に継続するサイクル数である。
In order to achieve the above object, the corrosion rate estimation device according to the first aspect of the present invention is
By performing multiple regression analysis with the corrosion rate of the metal material as the objective variable and the corrosive environment factors including the annual amount of sea salt particles flying in the year, the wetting time, and the dry-wet cycle that affect the corrosion of the metal material as the explanatory variables. An estimation formula creation method for creating an estimation formula for the corrosion rate of metal materials,
A corrosion rate estimation means for estimating the corrosion rate of the metal material by inputting a corrosion environment factor including an annual amount of sea salt particles, a wetting time, and a dry / wet cycle into the estimation formula.
Equipped with
The dry-wet cycle is the number of cycles in which a wet state for a predetermined time or longer and a dry state for a predetermined time or longer are alternately continued.

前記推定式は、以下の式に示すとおりであってもよい。
Y=α+β+β+β
ただし、Yは金属材料の腐食速度、Xは年間飛来海塩粒子量、Xは濡れ時間、Xは乾湿サイクル、αは定数、β~βは係数である。
The estimation formula may be as shown in the following formula.
Y = α + β 1 X 1 + β 2 X 2 + β 3 X 3
However, Y is the corrosion rate of the metal material, X 1 is the amount of sea salt particles flying annually, X 2 is the wetting time, X 3 is the dry / wet cycle, α is a constant, and β 1 to β 3 are coefficients.

前記推定式作成手段における年間飛来海塩粒子量、濡れ時間及び乾湿サイクルは、いずれも気象解析コードから得られたシミュレーション値であってもよい。 The annual amount of sea salt particles, the wetting time, and the dry / wet cycle in the estimation formula creating means may all be simulation values obtained from the meteorological analysis code.

前記金属材料は、鉄塔を構成する鋼鉄部材であり、
前記腐食速度は、前記鋼鉄部材に施された亜鉛めっき厚さの一年間当たりの減少量であってもよい。
The metal material is a steel member constituting a steel tower.
The corrosion rate may be a reduction in the galvanized thickness applied to the steel member per year.

前記金属材料は、送電線を構成するアルミ素線であり、
前記腐食速度は、前記アルミ素線の引張強度の一年間当たりの減少量であってもよい。
The metal material is an aluminum wire constituting a power transmission line.
The corrosion rate may be the amount of decrease in the tensile strength of the aluminum wire per year.

前記濡れ時間及び前記乾湿サイクルは、複数の定義のパターンを用意しておき、重回帰分析を行った場合に最も解析結果が良好な定義のパターンを選択してもよい。 For the wetting time and the dry-wet cycle, a plurality of defined patterns may be prepared, and the defined pattern with the best analysis result may be selected when the multiple regression analysis is performed.

推定された前記金属材料の腐食速度に基づいて、所定の地域における金属材料の腐食速度を示す腐食速度マップを作成する腐食速度マップ作成手段を備えてもよい。 A corrosion rate map creating means may be provided for creating a corrosion rate map showing the corrosion rate of the metal material in a predetermined area based on the estimated corrosion rate of the metal material.

上記目的を達成するために、本発明の第2の観点に係る腐食速度推定方法は、
金属材料の腐食速度を目的変数とし、前記金属材料の腐食に影響を与える年間飛来海塩粒子量、濡れ時間及び乾湿サイクルを含む腐食環境因子を説明変数とする重回帰分析を行うことで、前記金属材料の腐食速度の推定式を作成する推定式作成ステップと、
前記推定式に年間飛来海塩粒子量、濡れ時間及び乾湿サイクルを含む腐食環境因子を入力することで、前記金属材料の腐食速度を推定する腐食速度推定ステップと、
を含み、
乾湿サイクルは、所定時間以上の湿潤状態と所定時間以上の乾燥状態とが交互に継続するサイクル数である。
In order to achieve the above object, the corrosion rate estimation method according to the second aspect of the present invention is:
By performing multiple regression analysis with the corrosion rate of the metal material as the objective variable and the corrosive environment factors including the annual amount of sea salt particles flying in the year, the wetting time, and the dry-wet cycle that affect the corrosion of the metal material as the explanatory variables. An estimation formula creation step to create an estimation formula for the corrosion rate of metal materials, and
A corrosion rate estimation step for estimating the corrosion rate of the metal material by inputting a corrosion environment factor including an annual amount of sea salt particles, a wetting time, and a dry / wet cycle into the estimation formula.
Including
The dry-wet cycle is the number of cycles in which a wet state for a predetermined time or longer and a dry state for a predetermined time or longer are alternately continued.

上記目的を達成するために、本発明の第3の観点に係るプログラムは、
コンピュータを、
金属材料の腐食速度を目的変数とし、前記金属材料の腐食に影響を与える年間飛来海塩粒子量、濡れ時間及び乾湿サイクルを含む腐食環境因子を説明変数とする重回帰分析を行うことで、前記金属材料の腐食速度の推定式を作成する推定式作成手段、
前記推定式に年間飛来海塩粒子量、濡れ時間及び乾湿サイクルを含む腐食環境因子を入力することで、前記金属材料の腐食速度を推定する腐食速度推定手段、
として機能させるためのプログラムであって、
乾湿サイクルは、所定時間以上の湿潤状態と所定時間以上の乾燥状態とが交互に継続するサイクル数である。
In order to achieve the above object, the program according to the third aspect of the present invention is
Computer,
By performing multiple regression analysis with the corrosion rate of the metal material as the objective variable and the corrosive environment factors including the annual amount of sea salt particles flying in the year, the wetting time, and the dry-wet cycle that affect the corrosion of the metal material as the explanatory variables. An estimation formula creation method for creating an estimation formula for the corrosion rate of metal materials,
A corrosion rate estimation means for estimating the corrosion rate of the metal material by inputting the corrosion environment factors including the annual amount of sea salt particles, the wetting time and the dry / wet cycle into the estimation formula.
It is a program to function as
The dry-wet cycle is the number of cycles in which a wet state for a predetermined time or longer and a dry state for a predetermined time or longer are alternately continued.

本発明によれば、周囲環境の影響を考慮して金属材料の腐食速度を定量化することが可能な腐食速度推定装置、腐食速度推定方法及びプログラムを提供できる。 INDUSTRIAL APPLICABILITY According to the present invention, it is possible to provide a corrosion rate estimation device, a corrosion rate estimation method and a program capable of quantifying the corrosion rate of a metal material in consideration of the influence of the surrounding environment.

本発明の実施の形態に係る腐食速度推定装置のハードウェア構成を示すブロック図である。It is a block diagram which shows the hardware composition of the corrosion rate estimation apparatus which concerns on embodiment of this invention. (a)は、年間飛来海塩粒子量データテーブルの一例、(b)は、濡れ時間データテーブルの一例、(c)は、乾湿サイクルデータテーブルの一例、(d)は、腐食実績データテーブルの一例を示す図である。(A) is an example of an annual sea salt particle amount data table, (b) is an example of a wet time data table, (c) is an example of a dry / wet cycle data table, and (d) is an example of a corrosion record data table. It is a figure which shows an example. (a)~(c)は、乾湿サイクルの集計方法の具体例を説明するための図である。(A) to (c) are diagrams for explaining a specific example of the method of totaling the dry / wet cycle. 年間飛来海塩粒子量、濡れ時間及び乾湿サイクルの説明変数を説明するための図である。It is a figure for demonstrating the explanatory variables of the annual amount of sea salt particles, the wetting time, and the dry-wet cycle. 本発明の実施の形態に係る腐食速度推定処理の流れを示すフローチャートである。It is a flowchart which shows the flow of the corrosion rate estimation process which concerns on embodiment of this invention. 本発明の実施の形態に係る推定式作成処理の流れを示すフローチャートである。It is a flowchart which shows the flow of the estimation formula making process which concerns on embodiment of this invention. 本発明の実施の形態に係る外れ値除外処理の流れを示すフローチャートである。It is a flowchart which shows the flow of the outlier exclusion processing which concerns on embodiment of this invention. 実施例1における重回帰分析の結果を示す図である。It is a figure which shows the result of the multiple regression analysis in Example 1. FIG. 実施例1における鉄塔腐食速度マップを示す図である。It is a figure which shows the steel tower corrosion rate map in Example 1. FIG. 実施例2における重回帰分析の結果を示す図である。It is a figure which shows the result of the multiple regression analysis in Example 2. FIG. 実施例2における電線腐食速度マップを示す図である。It is a figure which shows the electric wire corrosion rate map in Example 2. FIG.

以下、本発明に係る腐食速度推定装置、腐食速度推定方法及びプログラムの実施の形態を、図面を参照しながら詳細に説明する。各図面においては、同一又は同等の部分に同一の符号を付す。以下、鉄塔の腐食速度を算出する場合を例に説明する。 Hereinafter, embodiments of the corrosion rate estimation device, the corrosion rate estimation method, and the program according to the present invention will be described in detail with reference to the drawings. In each drawing, the same or equivalent parts are designated by the same reference numerals. Hereinafter, a case of calculating the corrosion rate of a steel tower will be described as an example.

実施の形態に係る腐食速度推定方法では、鉄塔の腐食速度及び腐食環境因子のペアに基づいて、腐食環境因子から鉄塔の腐食速度を求める推定式を作成し、作成された推定式に各地点の位置の腐食環境因子の実測値を代入することで各地点の鉄塔の腐食速度を算出する。 In the corrosion rate estimation method according to the embodiment, an estimation formula for obtaining the corrosion rate of the steel tower from the corrosion environment factor is created based on the pair of the corrosion rate of the steel tower and the corrosion environment factor, and the created estimation formula is used for each point. The corrosion rate of the steel tower at each point is calculated by substituting the measured value of the corrosion environment factor at the position.

鉄塔の腐食速度を求める推定式の作成には、腐食速度を目的変数とし、腐食環境因子を説明変数とした重回帰分析を用いる。重回帰分析では、複数の説明変数X~Xを用いて目的変数Yを表す重回帰式を作成する。重回帰式は、αを定数、β~βを係数とすると、以下の式(1)で表される。
Y=α+β・X+β・X+…β・X …(1)
Multiple regression analysis is used with the corrosion rate as the objective variable and the corrosion environment factor as the explanatory variable to create the estimation formula for determining the corrosion rate of the tower. In the multiple regression analysis, a multiple regression equation representing the objective variable Y is created using a plurality of explanatory variables X 1 to X n . The multiple regression equation is expressed by the following equation (1), where α is a constant and β 1 to β n are coefficients.
Y = α + β 1・ X 1 + β 2・ X 2 +… β n・ X n … (1)

実施の形態に係る腐食速度推定方法では、目的変数Yである腐食速度と目的変数Yに対応する説明変数X~Xである腐食環境因子とを複数箇所で取得し、複数箇所で取得された目的変数Y及び説明変数X~Xに基づいて定数α及び係数β~βを決定することで、鉄塔の腐食速度を表す推定式(重回帰式)を作成する。 In the corrosion rate estimation method according to the embodiment, the corrosion rate, which is the objective variable Y, and the corrosive environment factors, which are the explanatory variables X1 to Xn corresponding to the objective variable Y, are acquired at a plurality of locations and acquired at a plurality of locations. By determining the constant α and the coefficients β 1 to β n based on the objective variable Y and the explanatory variables X 1 to X n , an estimation formula (multiple regression equation) expressing the corrosion rate of the steel tower is created.

腐食影響因子は、金属材料の腐食に影響を及ぼす環境因子であり、例えば、年間飛来海塩粒子量、濡れ時間、乾湿サイクルを含む。腐食環境因子は、気象解析コードから得られるシミュレーション値である。腐食環境因子は、シミュレーション値に限られず、例えば、観測機を用いて直接観測することで取得された実測値であってもよい。 Corrosion influencing factors are environmental factors that affect the corrosion of metallic materials and include, for example, the annual amount of sea salt particles, wetting time, and wet / dry cycle. The corrosive environmental factor is a simulation value obtained from the meteorological analysis code. The corrosive environmental factor is not limited to the simulated value, and may be, for example, an actually measured value obtained by direct observation using an observation aircraft.

腐食実績データは、屋外設置物からの採取時における部材の腐食の程度を示すデータである。例えば、屋外設置物が鉄塔の場合、腐食実績データは、鉄塔部材の採取時における亜鉛めっきの厚さ(μm)である。腐食実績データに基づいて、一年間に進行する対象部材の腐食の程度を示す指標である腐食速度を算出できる。鉄塔の部材の腐食速度(μm/年)は、亜鉛めっきの初期厚さから採取時点と亜鉛めっきの厚さとの差分を部材の設置年数で除算することで算出される。 Corrosion record data is data showing the degree of corrosion of members when collected from an outdoor installation. For example, when the outdoor installation is a steel tower, the corrosion record data is the thickness (μm) of the zinc plating at the time of collecting the steel tower member. Based on the actual corrosion data, the corrosion rate, which is an index showing the degree of corrosion of the target member progressing in one year, can be calculated. The corrosion rate (μm / year) of a member of a steel tower is calculated by dividing the difference between the time of collection and the thickness of zinc plating from the initial thickness of zinc plating by the number of years of installation of the member.

鉄塔の腐食速度Yの推定式(1)は、説明変数を年間飛来海塩粒子量X、濡れ時間X、乾湿サイクルXとすると、以下の式(2)で表される。
Y=α+β・X+β・X+β・X …(2)
The estimation formula (1) for the corrosion rate Y of the steel tower is expressed by the following formula (2), where the explanatory variables are the annual amount of sea salt particles X 1 , the wetting time X 2 , and the wet / dry cycle X 3 .
Y = α + β 1・ X 1 + β 2・ X 2 + β 3・ X 3 … (2)

次に、実施の形態に係る腐食速度推定装置1のハードウェア構成を説明する。図1は、実施の形態に係る腐食速度推定装置1のハードウェア構成を示すブロック図である。 Next, the hardware configuration of the corrosion rate estimation device 1 according to the embodiment will be described. FIG. 1 is a block diagram showing a hardware configuration of the corrosion rate estimation device 1 according to the embodiment.

腐食速度推定装置1は、例えば、汎用コンピュータである。腐食速度推定装置1は、操作部11と、表示部12と、通信部13と、記憶部14と、制御部15と、を備える。腐食速度推定装置1の各部は、内部バス等で相互に接続されている。 The corrosion rate estimation device 1 is, for example, a general-purpose computer. The corrosion rate estimation device 1 includes an operation unit 11, a display unit 12, a communication unit 13, a storage unit 14, and a control unit 15. Each part of the corrosion rate estimation device 1 is connected to each other by an internal bus or the like.

操作部11は、ユーザの指示を受け付け、受け付けた操作に対応する操作信号を制御部15に供給する。操作部11は、例えば、キーボード、マウスを備える。 The operation unit 11 receives a user's instruction and supplies an operation signal corresponding to the accepted operation to the control unit 15. The operation unit 11 includes, for example, a keyboard and a mouse.

表示部12は、制御部15から供給される画像データに基づいて、ユーザに向けて各種の画像を表示する。表示部12は、例えば、各地点の腐食速度の推定結果を表示する。 The display unit 12 displays various images to the user based on the image data supplied from the control unit 15. The display unit 12 displays, for example, the estimation result of the corrosion rate at each point.

通信部13は、アンテナ等を備え、インターネット等の通信ネットワークに接続することが可能なインターフェースである。通信部13は、外部の端末、サーバ、メモリ等とインターネット等の通信ネットワークを介して通信する。 The communication unit 13 is an interface provided with an antenna or the like and capable of connecting to a communication network such as the Internet. The communication unit 13 communicates with an external terminal, a server, a memory, or the like via a communication network such as the Internet.

記憶部14は、RAM(Random Access Memory)、ROM(Read Only Memory)、フラッシュメモリ、ハードディスクドライブ等を備える。記憶部14は、制御部15に実行されるプログラムや各種のデータを記憶する。また、記憶部14は、制御部15が処理を実行するためのワークメモリとして機能する。 The storage unit 14 includes a RAM (Random Access Memory), a ROM (Read Only Memory), a flash memory, a hard disk drive, and the like. The storage unit 14 stores programs and various data to be executed in the control unit 15. Further, the storage unit 14 functions as a work memory for the control unit 15 to execute the process.

さらに、記憶部14は、年間飛来海塩粒子量データテーブル141と、濡れ時間データテーブル142と、乾湿サイクルデータテーブル143と、腐食実績データテーブル144と、を備える。各データテーブル141~144は、それぞれ年間飛来海塩粒子量、濡れ時間、乾湿サイクルを位置情報に対応付けて記憶する。腐食実績データテーブル144は、腐食実績データを鉄塔の位置情報及びデータの取得日時に関する情報と対応付けて記憶する。 Further, the storage unit 14 includes an annual flying sea salt particle amount data table 141, a wet time data table 142, a dry / wet cycle data table 143, and a corrosion record data table 144. Each data table 141 to 144 stores the annual amount of sea salt particles flying in the sea, the wetting time, and the dry / wet cycle in association with the position information. The corrosion record data table 144 stores the corrosion record data in association with the position information of the tower and the information regarding the acquisition date and time of the data.

図2(a)は、年間飛来海塩粒子量データテーブル141の一例である。年間飛来海塩粒子量データテーブル141は、例えば、5km×5kmのメッシュで区切られた領域毎に年間飛来海塩粒子量データを記憶する。年間飛来海塩粒子量は、一年間に飛来した海塩の粒子量(μg/ms)の平均値である。年間飛来海塩粒子量は、例えば、一般財団法人電力中央研究所が開発した解析コードNuWiCC-STを用い、気流や塩分濃度等の計算条件を設定してシミュレーションを実行することで算出される。 FIG. 2A is an example of the annual sea salt particle amount data table 141. The annual sea salt particle amount data table 141 stores, for example, the annual sea salt particle amount data for each region divided by a mesh of 5 km × 5 km. The annual amount of sea salt particles flying is the average value of the amount of sea salt particles (μg / m 2 s) flying in one year. The annual amount of sea salt particles is calculated by, for example, using the analysis code NuWiCC-ST developed by the Central Research Institute of Electric Power Industry, setting calculation conditions such as airflow and salinity, and executing a simulation.

図2(b)は、濡れ時間データテーブル142の一例である。濡れ時間データテーブル142は、例えば、5km×5kmのメッシュで区切られた領域毎に濡れ時間データを記憶する。濡れ時間は、ISO(International Organization for Standardization)9223によると、気温0℃以上で相対湿度80%以上の時間(hour/年)である。濡れ時間は、例えば、一般財団法人電力中央研究所が開発した解析コードNuWFASに基づいて決定される。 FIG. 2B is an example of the wetting time data table 142. The wet time data table 142 stores, for example, wet time data for each region divided by a mesh of 5 km × 5 km. According to ISO (International Organization for Standardization) 9223, the wetting time is a time (hour / year) at a temperature of 0 ° C. or higher and a relative humidity of 80% or higher. The wetting time is determined, for example, based on the analysis code NuWFAS developed by the Central Research Institute of Electric Power Industry.

図2(c)は、乾湿サイクルデータテーブル143の一例である。乾湿サイクルデータテーブル143は、例えば、5km×5kmのメッシュで区切られた領域毎に乾湿サイクルデータを記憶する。乾湿サイクルは、一般的に、気温0℃以上かつ相対湿度80%以上の状態(湿潤状態)と湿潤状態以外の状態(乾燥状態)とが、それぞれ所定時間(例えば、12時間)以上、交互に継続した場合を1サイクルとする年間サイクル数(cycle/年)である。乾湿サイクルは、例えば、濡れ時間と同様に解析コードNuWFASに基づいて決定される。 FIG. 2C is an example of the dry / wet cycle data table 143. The wet / dry cycle data table 143 stores the dry / wet cycle data for each region divided by a mesh of, for example, 5 km × 5 km. In the dry-wet cycle, in general, a state where the temperature is 0 ° C. or higher and a relative humidity of 80% or more (wet state) and a state other than the wet state (dry state) are alternately performed for a predetermined time (for example, 12 hours) or more. It is the number of annual cycles (cycle / year) with the case of continuation as one cycle. The wet-wet cycle is determined, for example, based on the analysis code NuWFAS as well as the wetting time.

図2(d)は、腐食実績データテーブル144の一例である。腐食実績データテーブル144は、部材が設置されていた位置に関する情報及びデータ取得日時に関する情報と対応付けて、腐食実績データであるデータの取得日時における部材の亜鉛めっき厚さ(μm)を記憶する。 FIG. 2D is an example of the corrosion actual data table 144. The corrosion record data table 144 stores the zinc plating thickness (μm) of the member at the acquisition date and time of the data which is the corrosion record data in association with the information regarding the position where the member was installed and the information regarding the data acquisition date and time.

図3は、乾湿サイクルの集計方法の具体例を示す図である。図3(a)~(c)における横線から垂直に延びる複数の縦線は、1時間毎の時間の区切りを意味する。(a)の場合は、乾湿サイクルにおける湿潤状態が継続した時間(定義時間)が12時間以上であり、その後の乾湿サイクルにおける乾燥状態が継続した時間(非定義時間)が12時間以上であるため、乾湿サイクルは1サイクルである。他方、(b)及び(c)の場合は、定義時間又は非定義時間が12時間未満であるため、乾湿サイクルはゼロである。 FIG. 3 is a diagram showing a specific example of a method for tabulating the dry / wet cycle. A plurality of vertical lines extending vertically from the horizontal lines in FIGS. 3 (a) to 3 (c) mean an hourly time division. In the case of (a), the time during which the wet state continues in the dry-wet cycle (defined time) is 12 hours or more, and the time during which the dry state continues in the subsequent dry-wet cycle (undefined time) is 12 hours or more. , The dry / wet cycle is one cycle. On the other hand, in the case of (b) and (c), the dry-wet cycle is zero because the defined time or the undefined time is less than 12 hours.

図4は、年間飛来海塩粒子量、濡れ時間及び乾湿サイクルの説明変数を説明するための図である。年間飛来海塩粒子量は、粒子量という物理量であるため一義的に定義されるが、濡れ時間及び乾湿サイクルについては、様々な定義づけを行うことが可能であり、腐食速度推定を精度よく行うためには、一般的な定義に縛られずに適切な定義づけを行う必要がある。このため、濡れ時間及び乾湿サイクルについては、上述した一般的な定義を拡張して、それぞれ図4に示す9種類の濡れ時間及び乾湿サイクルを説明変数とする。 FIG. 4 is a diagram for explaining the explanatory variables of the annual amount of sea salt particles, the wetting time, and the dry / wet cycle. The annual amount of sea salt particles is uniquely defined because it is a physical quantity called the amount of particles. However, various definitions can be made for the wetting time and the dry-wet cycle, and the corrosion rate can be estimated accurately. In order to do so, it is necessary to make an appropriate definition without being bound by a general definition. Therefore, regarding the wetting time and the dry-wet cycle, the above-mentioned general definitions are extended to use the nine types of wet time and the dry-wet cycle shown in FIG. 4, respectively, as explanatory variables.

濡れ時間は、相対湿度だけではなく、降水量0mm以上の場合も加味され、乾湿サイクルは、湿潤状態と乾燥状態の繰り返しサイクル数だけではなく、降水量0mm以上の期間と降水量0mmの期間との繰り返しサイクル数も加味される。このため、以下の処理では、9通りの濡れ時間の説明変数と9通りの乾湿サイクルの説明変数との組み合わせ(パターン)、言い換えると計81通りの説明変数の組み合わせを用いて重回帰分析を実行することになる。なお、説明変数の組み合わせ数は計81種類に限られず、適切な種類の説明変数の組み合わせを作成できる。 The wetting time is taken into consideration not only in relative humidity but also in the case of precipitation of 0 mm or more, and the dry-wet cycle includes not only the number of repeated cycles of wet and dry states, but also the period of precipitation of 0 mm or more and the period of precipitation of 0 mm. The number of repeated cycles of is also taken into consideration. Therefore, in the following processing, multiple regression analysis is executed using a combination (pattern) of 9 explanatory variables of wet time and 9 explanatory variables of the dry-wet cycle, in other words, a total of 81 combinations of explanatory variables. Will be done. The number of combinations of explanatory variables is not limited to 81 in total, and appropriate combinations of explanatory variables can be created.

図1に戻り、制御部15は、CPU(Central Processing Unit)等を備え、腐食速度推定装置1の各部の制御を行う。制御部15は、記憶部14に記憶されているプログラムを実行することにより、図5の腐食速度推定処理、図6の推定式作成処理及び図7の外れ値除外処理を実行する。 Returning to FIG. 1, the control unit 15 includes a CPU (Central Processing Unit) and the like, and controls each unit of the corrosion rate estimation device 1. By executing the program stored in the storage unit 14, the control unit 15 executes the corrosion rate estimation process of FIG. 5, the estimation formula creation process of FIG. 6, and the outlier exclusion process of FIG. 7.

制御部15は、機能的には、推定式作成部151と、腐食速度推定部152と、腐食速度マップ作成部153と、を備える。 The control unit 15 functionally includes an estimation formula creation unit 151, a corrosion rate estimation unit 152, and a corrosion rate map creation unit 153.

推定式作成部151は、記憶部14に記憶された鉄塔部材の腐食実績データ及び複数種類の腐食環境因子のペアに基づいて重回帰分析を実行することで、複数種類の腐食環境因子から鉄塔の腐食速度を求める推定式を作成する。 The estimation formula creation unit 151 performs a multiple regression analysis based on the corrosion record data of the tower member stored in the storage unit 14 and a pair of a plurality of types of corrosion environment factors, so that the tower can be converted from a plurality of types of corrosion environment factors. Create an estimation formula to calculate the corrosion rate.

推定式作成部151は、計81通りの組み合わせの説明変数を用いてそれぞれ重回帰分析を実行し、計81通りの説明変数の組み合わせから最も適切な濡れ時間及び乾湿サイクルの説明変数を選択する。具体的には、計81通りの説明変数の組み合わせの中から、有意確率がゼロであって決定係数Rが最も高い濡れ時間及び乾湿サイクルの説明変数の組み合わせを選択する。 The estimation formula creation unit 151 performs multiple regression analysis using a total of 81 combinations of explanatory variables, and selects the most appropriate wet time and wet / wet cycle explanatory variables from the total of 81 combinations of explanatory variables. Specifically, from a total of 81 combinations of explanatory variables, the combination of the explanatory variables of the wet time and the wet-wet cycle with the highest significance probability of zero and the highest coefficient of determination R2 is selected.

決定係数Rは回帰モデルの当てはまりのよさを示す指標である。決定係数Rは、0≦R≦1の範囲内であり、1に近づくほど作成した重回帰式の精度が高いことを意味している。有意確率は、P値とも呼ばれ、決定係数の有意性検定に用いられている指標である。有意確率はゼロに近づくほど有意性ありと判断でき、一般的には、0.05よりも小さい場合に、有意性ありと判断できる。 The coefficient of determination R 2 is an index showing the goodness of the fit of the regression model. The coefficient of determination R 2 is within the range of 0 ≦ R 2 ≦ 1, and the closer it is to 1, the higher the accuracy of the created multiple regression equation. The significance probability, also called a P-value, is an index used for the significance test of the coefficient of determination. The significance probability can be judged to be significant as it approaches zero, and generally, it can be judged to be significant when it is smaller than 0.05.

推定式作成部151は、最も精度が良い説明変数の組み合わせを含む目的変数及び説明変数のペアから、外れ値の影響力評価に使用されているDfBetas、Cookの距離、てこ比の3つの指標を用いて外れ値を除外して重回帰分析を実行し、3つの指標のうち最も結果のよい指標を用いた重回帰分析の結果を選択する。 From the pair of objective variable and explanatory variable including the most accurate combination of explanatory variables, the estimation formula creation unit 151 obtains three indexes of DfBetas, Cook distance, and lever ratio used for evaluating the influence of outliers. The outliers are excluded and the multiple regression analysis is performed, and the result of the multiple regression analysis using the index with the best result among the three indicators is selected.

DfBetasは、ある1つのデータを取り除いたときの回帰係数の変化を表す。DfBetasの絶対値が2/√N(Nは総サンプル数)を超えるとき、当該データは外れ値とみなす。 DfBetas represents the change in the regression coefficient when one piece of data is removed. When the absolute value of DfBetas exceeds 2 / √N (N is the total number of samples), the data is considered outlier.

Cookの距離は、全てのデータを用いた場合と1つのデータを除いた後に求めた回帰式による予測値を用いた場合との食い違いに関する距離を表す。Cookの距離が大きいと、除かれた1つのデータが回帰式による予測値に大きく影響していることを意味する。 The Cook distance represents the distance related to the discrepancy between the case where all the data are used and the case where the predicted value by the regression equation obtained after excluding one data is used. If the distance of Cook is large, it means that one of the excluded data has a great influence on the predicted value by the regression equation.

てこ比は、回帰分析の観察点(サンプル)ごとに、説明変数のデータを変えずに目的変数Yの値を1だけ変えたときの予測値の変化量を表す。てこ比が大きいと、目的変数Yの値を1だけ変えられたデータが回帰式による予測値に大きく影響していることを意味する。 The lever ratio represents the amount of change in the predicted value when the value of the objective variable Y is changed by 1 without changing the data of the explanatory variables for each observation point (sample) of the regression analysis. When the leverage ratio is large, it means that the data in which the value of the objective variable Y is changed by 1 has a great influence on the predicted value by the regression equation.

腐食速度推定部152は、推定式作成部151で作成された推定式に各地点の腐食環境因子を代入することで、各地点の鉄塔部材の腐食速度を推定する。 The corrosion rate estimation unit 152 estimates the corrosion rate of the tower member at each point by substituting the corrosion environmental factor at each point into the estimation formula created by the estimation formula creation unit 151.

腐食速度マップ作成部153は、腐食速度推定部152により推定された各地点の腐食速度に基づいて、腐食速度マップを作成する。腐食速度マップは、例えば、図9及び図11に示すように、対象領域の各地点を腐食速度に応じた区分された色で表示したマップである。以上が、腐食速度推定装置1のハードウェア構成である。 The corrosion rate map creating unit 153 creates a corrosion rate map based on the corrosion rate at each point estimated by the corrosion rate estimation unit 152. As shown in FIGS. 9 and 11, for example, the corrosion rate map is a map in which each point in the target area is displayed in a color divided according to the corrosion rate. The above is the hardware configuration of the corrosion rate estimation device 1.

図5は、実施の形態に係る腐食速度推定方法の流れを示すフローチャートである。以下、複数の地点で取得された年間飛来海塩粒子量X、濡れ時間X、乾湿サイクルXに関するデータと、これらのデータが取得された各地点で取得された腐食速度Yに関するデータとは、いずれも既に記憶部14に記憶されているものとする。 FIG. 5 is a flowchart showing the flow of the corrosion rate estimation method according to the embodiment. Below, the data on the annual amount of sea salt particles X 1 , the wetting time X 2 , and the dry / wet cycle X 3 acquired at multiple points, and the data on the corrosion rate Y acquired at each point where these data were acquired. Are already stored in the storage unit 14.

まず、制御部15は、記憶部14から説明変数である年間飛来海塩粒子量X、濡れ時間X、乾湿サイクルXを読み出す(ステップS1)。濡れ時間X、乾湿サイクルXについては、それぞれ図4に示す9種類の説明変数を用いた場合のデータを取得する。 First, the control unit 15 reads out the annual flying sea salt particle amount X 1 , the wetting time X 2 , and the dry / wet cycle X 3 which are explanatory variables from the storage unit 14 (step S1). For the wet time X 2 and the dry / wet cycle X 3 , data are acquired when the nine types of explanatory variables shown in FIG. 4 are used.

次に、制御部15は、記憶部14から目的変数である鉄塔の腐食速度Yを読み出す(ステップS2)。鉄塔の腐食速度Yは、記憶部14に記憶された鉄塔部材の採取時の亜鉛めっき厚さと亜鉛めっきの初期厚さとを読み出し、亜鉛めっきの厚さの差分を鉄塔部材の設置年数で除算することで算出される。鉄塔部材の設置年数は、記憶部14に記憶された鉄塔部材の設置日時と当該鉄塔部材の採取日時との差分を算出することで得られる。 Next, the control unit 15 reads out the corrosion rate Y of the steel tower, which is the objective variable, from the storage unit 14 (step S2). For the corrosion rate Y of the iron tower, read out the zinc plating thickness at the time of collecting the iron tower member stored in the storage unit 14 and the initial thickness of the zinc plating, and divide the difference in the zinc plating thickness by the number of years of installation of the iron tower member. It is calculated by. The number of years of installation of the tower member is obtained by calculating the difference between the installation date and time of the tower member stored in the storage unit 14 and the collection date and time of the tower member.

次に、推定式作成部151は、ステップS1で取得された各地点の説明変数X~X、ステップS2で取得された各地点の目的変数Yに基づいて重回帰分析を実行し、式(2)の定数α及び係数β~βを算出することで、腐食速度Yの推定式を作成する(ステップS3)。以下、図6のフローチャートを参照して、ステップS3の推定式を作成する処理(推定式作成処理)を説明する。 Next, the estimation formula creation unit 151 executes multiple regression analysis based on the explanatory variables X1 to X3 of each point acquired in step S1 and the objective variable Y of each point acquired in step S2, and formulates the equation. By calculating the constant α and the variables β 1 to β 3 in (2), an estimation formula for the corrosion rate Y is created (step S3). Hereinafter, the process of creating the estimation formula in step S3 (estimation formula creation process) will be described with reference to the flowchart of FIG.

まず、推定式作成部151は、記憶部14から読み出された計81種類の説明変数の組み合わせのそれぞれについて重回帰分析を実行し、有意確率がゼロで決定係数Rが最も1に近づいた説明変数の組み合わせを選択する(ステップS31)。 First, the estimation formula creation unit 151 executed multiple regression analysis for each of a total of 81 combinations of explanatory variables read from the storage unit 14, and the significance probability was zero and the coefficient of determination R2 was closest to 1. Select a combination of explanatory variables (step S31).

次に、推定式作成部151は、腐食実績データから外れ値を除外する処理を実行する(ステップS32)。以下、図7のフローチャートを参照して、ステップS32の外れ値を除外する処理(外れ値除外処理)を説明する。 Next, the estimation formula creation unit 151 executes a process of excluding outliers from the actual corrosion data (step S32). Hereinafter, the process of excluding the outliers in step S32 (outlier exclusion process) will be described with reference to the flowchart of FIG. 7.

図7は、推定式作成部151が実行する腐食実績データの外れ値除外処理の流れを示すフローチャートである。まず、推定式作成部151は、記憶部14から読み出された全ての腐食実績データ(データ数:N)に基づいて重回帰分析を行う(ステップS321)。 FIG. 7 is a flowchart showing a flow of outlier exclusion processing of corrosion actual data executed by the estimation formula creation unit 151. First, the estimation formula creation unit 151 performs multiple regression analysis based on all the corrosion record data (number of data: N) read from the storage unit 14 (step S321).

次に、推定式作成部151は、腐食実績データのDfBetasを算出し、算出されたDfBetasによるα個の外れ値を除外し、DfBetasによる外れ値を除外した腐食実績データ(データ数:N-α)に基づいて重回帰分析を行う(ステップS322)。 Next, the estimation formula creation unit 151 calculates DfBetas of the actual corrosion data, excludes α outliers due to the calculated DfBetas, and excludes the outliers due to DfBetas (number of data: N-α). ) Is performed (step S322).

次に、推定式作成部151は、腐食実績データのCookの距離を算出し、算出されたCookの距離による外れ値を上位α個まで除外し、Cookの距離による外れ値を除外した腐食実績データ(データ数:N-α)に基づいて重回帰分析を行う(ステップS323)。 Next, the estimation formula creation unit 151 calculates the Cook distance of the corrosion actual data, excludes the calculated outliers due to the Cook distance to the top α, and excludes the outliers due to the Cook distance. Multiple regression analysis is performed based on (number of data: N—α) (step S323).

次に、推定式作成部151は、腐食実績データのてこ比を算出し、算出されたてこ比による外れ値を上位α個まで除外し、てこ比による外れ値を除外した腐食実績データ(データ数:N-α)に基づいて重回帰分析を行う(ステップS324)。 Next, the estimation formula creation unit 151 calculates the lever ratio of the corrosion actual data, excludes the calculated outliers due to the iron ratio to the top α, and excludes the outliers due to the iron ratio (number of data). : N—α) is used for multiple regression analysis (step S324).

次に、推定式作成部151は、ステップS322~ステップS324で実行された重回帰分析の結果として得られた各決定係数Rを比較し、決定係数Rが最も1に近づいた外れ値の指標を選択し、当該指標を用いた場合の外れ値を腐食実績データから除外する(ステップS325)。以上が、ステップS32の腐食実績データの外れ値を除外する処理の流れである。 Next, the estimation formula creation unit 151 compares each coefficient of determination R2 obtained as a result of the multiple regression analysis executed in steps S322 to S324, and the outliers in which the coefficient of determination R2 is closest to 1. An index is selected, and outliers when the index is used are excluded from the actual corrosion data (step S325). The above is the flow of the process for excluding the outliers of the corrosion record data in step S32.

図6に戻り、推定式作成部151は、ステップS32の外れ値除外処理で除外された腐食実績データと当該腐食実績データに対応する計81種類の中から選択された説明変数とを用いて再回帰分析を実行し、式(2)の定数α及び係数β~βを算出することで、腐食速度Yの推定式を作成する(ステップS33)。以上が、ステップS3の推定式作成処理である。 Returning to FIG. 6, the estimation formula creation unit 151 re-uses the corrosion actual data excluded by the outlier exclusion process in step S32 and the explanatory variables selected from a total of 81 types corresponding to the corrosion actual data. By executing regression analysis and calculating the constant α and the coefficients β 1 to β 3 of the equation (2), an estimation equation of the corrosion rate Y is created (step S33). The above is the estimation formula creation process in step S3.

図5に戻り、腐食速度推定部152は、ステップS3で作成された腐食速度の推定式に、各データテーブル141~143に記憶された各地点の説明変数X~Xを入力し、各地点における鉄塔部材の腐食速度Yを推定する(ステップS4)。 Returning to FIG. 5, the corrosion rate estimation unit 152 inputs the explanatory variables X1 to X3 of each point stored in the data tables 141 to 143 into the corrosion rate estimation formula created in step S3 , and inputs the explanatory variables X1 to X3 at each location. The corrosion rate Y of the steel tower member at the point is estimated (step S4).

次に、腐食速度マップ作成部153は、ステップS4で推定された各地点の腐食速度Yに基づいて、所定領域の腐食速度の分布を示す腐食速度マップを作成する(ステップS5)。ステップS5で作成された腐食速度マップは、腐食速度推定装置1の表示部12に表示され、ユーザによる鉄塔の補修計画の策定に利用される。以上が、腐食速度推定処理の一連の流れである。 Next, the corrosion rate map creating unit 153 creates a corrosion rate map showing the distribution of the corrosion rate in a predetermined region based on the corrosion rate Y at each point estimated in step S4 (step S5). The corrosion rate map created in step S5 is displayed on the display unit 12 of the corrosion rate estimation device 1 and is used by the user to formulate a repair plan for the steel tower. The above is a series of steps of the corrosion rate estimation process.

実施の形態に係る腐食速度推定装置1は、鉄塔の腐食速度を目的変数とし、鉄塔の腐食に影響を与える年間飛来海塩粒子量、濡れ時間、乾湿サイクルを説明変数とする重回帰分析を行うことで、鉄塔の腐食速度の推定式を作成するように構成されている。このため、周囲環境の影響を加味して鉄塔の腐食速度を精度よく推定することができる。 The corrosion rate estimation device 1 according to the embodiment performs multiple regression analysis with the corrosion rate of the steel tower as the objective variable and the annual amount of sea salt particles flying in the sea, the wetting time, and the dry / wet cycle that affect the corrosion of the steel tower as explanatory variables. Therefore, it is configured to create an estimation formula for the corrosion rate of the steel tower. Therefore, the corrosion rate of the steel tower can be estimated accurately in consideration of the influence of the surrounding environment.

以下、実施例を挙げて本発明を具体的に説明する。ただし、本発明はこれらの実施例に限定されるものではない。 Hereinafter, the present invention will be specifically described with reference to examples. However, the present invention is not limited to these examples.

(実施例1)
実施例1では、既設鉄塔から採取した部材から取得した腐食実績データに基づき重回帰分析を行うことで鉄塔の腐食速度マップの作成を作成した。重回帰分析には、例えば、統計解析ソフトウェア“Statistical Package for Social Science Ver.21”を用いた。
(Example 1)
In Example 1, a corrosion velocity map of the tower was created by performing multiple regression analysis based on the corrosion record data obtained from the members collected from the existing tower. For the multiple regression analysis, for example, statistical analysis software "Statistical Package for Social Science Ver.21" was used.

まず、腐食実績データである部材採取時の亜鉛めっき厚さの測定方法を説明する。複数地点の鉄塔より採取した部材を約50cmの長さにカットし、写真を撮影した後、部材から腐食生成物を除去し、磁気膜厚計により腐食性生物を含む亜鉛めっき厚さを測定した(150箇所/本)。カットされた部材上でXYZ駆動アームに支持された膜厚計を走査させ、部材の150箇所で膜厚を測定した。部材一本当たり150箇所の亜鉛めっき厚さを測定したのは、亜鉛めっき厚さの平均値を取得するためである。 First, a method of measuring the zinc plating thickness at the time of collecting members, which is the actual corrosion data, will be described. Members collected from steel towers at multiple points were cut to a length of about 50 cm, photographs were taken, corrosion products were removed from the members, and the thickness of zinc plating containing corrosive organisms was measured with a magnetic film thickness meter. (150 locations / book). A film thickness meter supported by the XYZ drive arm was scanned on the cut member, and the film thickness was measured at 150 points on the member. The reason why the zinc plating thickness at 150 points per member was measured is to obtain the average value of the zinc plating thickness.

次に、部材の亜鉛めっきの初期厚さについて検討した。鉄塔部材の亜鉛めっきの初期厚さは、少なくとも規定値以上であるが、実際の厚さについては不明である。このため、亜鉛めっきの初期厚さを部材の腐食の影響を受けないと推定される部分(ラップ部)における実測値と仮定し、推定式による腐食速度計算値と腐食実績データとの際が大きい場には、初期値を100μmと仮定した。 Next, the initial thickness of the zinc plating of the member was examined. The initial thickness of the zinc plating of the tower member is at least the specified value or more, but the actual thickness is unknown. Therefore, assuming that the initial thickness of zinc plating is the measured value in the part (wrapped part) that is estimated not to be affected by the corrosion of the member, there is a large difference between the corrosion rate calculation value by the estimation formula and the actual corrosion data. In the field, the initial value was assumed to be 100 μm.

図8は、鉄塔部材の腐食実績データに基づき重回帰分析を行った結果を示す。DfBetas、Cookの距離、てこ比のいずれの指標を採用した場合であっても、同一の外れ値を除外する結果となったため、当該外れ値のみを除外して腐食速度の推定式を作成した。この場合の決定係数Rは、0.403、有意確率は、0.000であった。 FIG. 8 shows the results of multiple regression analysis based on the actual corrosion data of the tower member. Regardless of which index of DfBetas, Cook distance, and lever ratio was adopted, the same outliers were excluded. Therefore, the estimation formula of the corrosion rate was created by excluding only the outliers. In this case, the coefficient of determination R 2 was 0.403 and the significance probability was 0.000.

鉄塔の腐食速度推定式は、以下の式(3)で表される。
鉄塔部材の腐食速度(μm/年)=-0.774545+1.211769×年間飛来海塩粒子量+0.000159×濡れ時間(ISO又は降水量0mm以上)+0.003600×乾湿サイクル(ISO又は0℃以上、降水量0mm以上、6h) …(3)
The corrosion rate estimation formula for the steel tower is expressed by the following formula (3).
Corrosion rate of steel tower members (μm / year) = -0.774545 + 1.211769 x Annual amount of sea salt particles +0.000159 x Wetting time (ISO or precipitation 0 mm or more) +0.003600 x Wet and dry cycle (ISO or 0 ° C or more) , Precipitation 0 mm or more, 6h) ... (3)

図9は、腐食速度を明確に区分するために腐食速度を速、中、遅の3段階に分類した鉄塔腐食速度マップを示す。図9の腐食速度の単位は、μm/年であり、メッシュの単位は、5km×5kmである。鉄塔腐食速度マップから、鉄塔の腐食速度は、日本海側で速い傾向にあり、太平洋側やオホーツク海側で遅い傾向があることが理解できる。また、道北地方や道南地方、道央地方の一部では、内陸部においても腐食速度が速い傾向にあることも理解できる。これらの結果は、腐食実績データと合致することも確認できた。 FIG. 9 shows a tower corrosion rate map in which the corrosion rate is classified into three stages of high, medium, and slow in order to clearly classify the corrosion rate. The unit of the corrosion rate in FIG. 9 is μm / year, and the unit of the mesh is 5 km × 5 km. From the tower corrosion rate map, it can be understood that the tower corrosion rate tends to be faster on the Sea of Japan side and slower on the Pacific Ocean side and the Sea of Okhotsk side. It can also be understood that in some parts of Northern Hokkaido, Southern Hokkaido, and Central Hokkaido, the corrosion rate tends to be high even in the inland areas. It was also confirmed that these results are in agreement with the actual corrosion data.

図9における腐食速度の分類は、亜鉛めっきの初期厚さの平均値、亜鉛めっきの寿命、鉄塔の電気事業固定資産耐用年数、亜鉛めっきの断面組織で純亜鉛層であるη層の厚さ等を考慮して、以下のように設定している。
腐食速度「速」 1.49μm/年以上
腐食速度「中」 0.66μm/年~1.49μm/年未満
腐食速度「遅」 0.66μm/年未満
The classification of the corrosion rate in FIG. 9 is the average value of the initial thickness of zinc plating, the life of zinc plating, the useful life of the electric business fixed assets of the steel tower, the thickness of the η layer which is a pure zinc layer in the cross-sectional structure of galvanization, etc. In consideration of, the settings are as follows.
Corrosion rate "fast" 1.49 μm / year or more Corrosion rate "medium" 0.66 μm / year to less than 1.49 μm / year Corrosion rate “slow” less than 0.66 μm / year

(実施例2)
実施例2では、既設送電線から採取した鋼心アルミ撚り線(電線)から取得した腐食実績データに基づき重回帰分析を行うことで電線腐食速度マップを作成した。重回帰分析には、実施例1の場合と同一の統計解析ソフトウェアを用いた。
(Example 2)
In Example 2, a wire corrosion rate map was created by performing multiple regression analysis based on the corrosion record data obtained from the steel core aluminum stranded wire (electric wire) collected from the existing transmission line. The same statistical analysis software as in Example 1 was used for the multiple regression analysis.

電線の腐食実績データは、鋼心アルミ撚り線から採取したアルミ素線の引張強さ(MPa)である。アルミ素線の引張強さを腐食実績データとしたのは、アルミ素線が鋼素線よりも経年や環境要因に起因して劣化しやすく、鋼心アルミ撚り線の余寿命をより適切に推定できるためである。なお、防食材が塗布された防食電線については、推定式の精度に影響を及ぼすことを避けるため、腐食実績データから除外している。 The actual corrosion data of the electric wire is the tensile strength (MPa) of the aluminum wire collected from the steel core aluminum stranded wire. The reason why the tensile strength of the aluminum wire is used as the corrosion record data is that the aluminum wire is more likely to deteriorate due to aging and environmental factors than the steel wire, and the remaining life of the steel core aluminum stranded wire is estimated more appropriately. Because it can be done. Corrosion-proof wires coated with anti-corrosion ingredients are excluded from the actual corrosion data in order to avoid affecting the accuracy of the estimation formula.

電線の腐食実績データは、安全側に評価するため、アルミ素線の引張試験を複数回実行し、測定された引張強さの最小値とした。アルミ素線の初期状態の引張強さと試験時の引張強さの最小値との差分を、電線設置後の経過年数で除算することで、アルミ素線の腐食速度(MPa/年)を算出した。なお、アルミ素線の初期状態の引張強さは、規格値の1.2倍とした。 In order to evaluate the corrosion record data of the electric wire on the safe side, the tensile test of the aluminum wire was executed multiple times, and the measured tensile strength was set to the minimum value. The corrosion rate (MPa / year) of the aluminum wire was calculated by dividing the difference between the initial tensile strength of the aluminum wire and the minimum value of the tensile strength at the time of the test by the number of years elapsed after the wire was installed. .. The tensile strength of the aluminum wire in the initial state was 1.2 times the standard value.

図10は、電線の腐食実績データに基づき重回帰分析を行った結果を示す。てこ比により6つの外れ値を除外した場合で最も良好な解析結果が得られため、この6つの外れ値のみを除外して腐食速度の推定式を作成した。この場合の決定係数Rは、0.228、有意確率は、0.001であった。 FIG. 10 shows the result of multiple regression analysis based on the actual corrosion data of the electric wire. Since the best analysis results were obtained when 6 outliers were excluded by the lever ratio, an estimation formula for the corrosion rate was created by excluding only these 6 outliers. In this case, the coefficient of determination R 2 was 0.228 and the significance probability was 0.001.

電線の腐食速度推定式は、以下の式(4)で表される。
電線の腐食速度(MPa/年)=-0.972796+5.169522×年間飛来海塩粒子量+0.000172×濡れ時間(温度0℃以上で湿度80%以上又は降水量0mm以上)-0.039616×乾湿サイクル(温度0℃以上で湿度80%以上又は降水量0mm以上、12h) …(4)
The corrosion rate estimation formula for electric wires is expressed by the following formula (4).
Corrosion rate of electric wire (MPa / year) = -0.972796 + 5.169522 x Annual amount of sea salt particles +0.000172 x Wetting time (humidity 80% or more or precipitation 0 mm or more at temperature 0 ° C or higher) -0.039616 x Wet and dry cycle (temperature 0 ° C or higher, humidity 80% or higher or precipitation 0 mm or higher, 12h) ... (4)

図11は、腐食速度を明確に区分するために腐食速度を速、中、遅の3段階に分類した電線の腐食速度マップを示す。図11の腐食速度の単位は、MPa/年である。電線の腐食速度は、年間飛来海塩粒子量に大きく依存しており、海岸付近で速く、内陸部で遅いことが理解できる。また、日本海側では、腐食速度が速く、太平洋側やオホーツク海側では、腐食速度が遅いことが理解できる。さらに、稚内地方では、内陸部まで海塩粒子が飛来するため、広範囲に亘って腐食速度が速いことが理解できる。これらの結果は、腐食実績データに合致することも確認できた。 FIG. 11 shows a corrosion rate map of an electric wire in which the corrosion rate is classified into three stages of high, medium, and slow in order to clearly classify the corrosion rate. The unit of corrosion rate in FIG. 11 is MPa / year. It can be understood that the corrosion rate of electric wires greatly depends on the amount of sea salt particles flying annually, and is fast near the coast and slow in the inland area. In addition, it can be understood that the corrosion rate is high on the Sea of Japan side and slow on the Pacific Ocean side and the Sea of Okhotsk side. Furthermore, in the Wakkanai region, it can be understood that the corrosion rate is high over a wide area because sea salt particles fly to the inland area. It was also confirmed that these results are in agreement with the actual corrosion data.

図11における腐食速度の分類は、鋼心アルミ撚り線ACSR160mmの初期値及び規格値、電線の電気事業固定資産耐用年数等を考慮して、以下のように設定している。
腐食速度「速」 1.69MPa/年以上
腐食速度「中」 0.845MPa/年~1.69MPa/年未満
腐食速度「遅」 0.845MPa/年未満
The classification of the corrosion rate in FIG. 11 is set as follows in consideration of the initial value and standard value of the steel core aluminum stranded wire ACSR160 mm 2 , the useful life of the electric business fixed assets of the electric wire, and the like.
Corrosion rate "fast" 1.69 MPa / year or more Corrosion rate "medium" 0.845 MPa / year to less than 1.69 MPa / year Corrosion rate "slow" less than 0.845 MPa / year

なお、図9及び図11では、腐食速度を速、中、遅の3つに区分しているが、あくまで例示であって、腐食速度を2つに区分してもよく、4つ以上に区分してもよい。 In addition, in FIG. 9 and FIG. 11, the corrosion rate is classified into three categories of fast, medium, and slow, but this is just an example, and the corrosion rate may be classified into two or more. You may.

本発明は上記の実施形態に限られず、以下に述べる変形も可能である。 The present invention is not limited to the above embodiment, and the modifications described below are also possible.

(変形例)
上記実施の形態では、説明変数として年間飛来海塩粒子量、濡れ時間、乾湿サイクルを用いていたが、本発明はこれに限られない。年間飛来海塩粒子量、濡れ時間、乾湿サイクルに加えて、他の腐食環境因子を説明変数として追加してもよい。
(Modification example)
In the above embodiment, the annual amount of sea salt particles flying in the sea, the wetting time, and the dry / wet cycle are used as explanatory variables, but the present invention is not limited to this. In addition to the annual amount of sea salt particles, wetting time, and dry-wet cycle, other corrosive environmental factors may be added as explanatory variables.

上記実施の形態では、鉄塔を構成する鋼鉄部材に施された亜鉛めっき厚さの一年間当たりの減少量を腐食速度としていたが、本発明はこれに限られない。例えば、送電線の腐食を評価する場合、アルミ素線の引張強度の一年間当たりの減少量を腐食速度としてもよい。 In the above embodiment, the amount of decrease in the galvanized thickness applied to the steel member constituting the steel tower per year is defined as the corrosion rate, but the present invention is not limited to this. For example, when evaluating the corrosion of a transmission line, the amount of decrease in the tensile strength of the aluminum wire per year may be used as the corrosion rate.

上記実施の形態では、81種類の説明変数の組み合わせから説明変数を選択し、選択された説明変数を用いて重回帰分析を行うことで外れ値の指標を選択していたが、本発明はこれに限られない。例えば、説明変数と外れ値を除外していない全ての腐食実績データ(仮の目的変数)とに基づいて多変量解析を実行し、多変量解析の結果からDfBatas、Cookの距離、てこ比のそれぞれを用いて腐食実績データの外れ値を除外することで、説明変数(81通り)及び目的変数の決め方(3通り)からなる全243通りの組み合わせから最適な目的変数を決定してもよい。 In the above embodiment, an explanatory variable is selected from a combination of 81 types of explanatory variables, and an index of outliers is selected by performing multiple regression analysis using the selected explanatory variable. Not limited to. For example, a multivariate analysis is performed based on the explanatory variables and all the actual corrosion data (provisional objective variables) that do not exclude the outliers, and the DfBatas, Cook distances, and the lever ratios are obtained from the results of the multivariate analysis. By excluding the outliers of the actual corrosion data using the above, the optimum objective variable may be determined from all 243 combinations consisting of the explanatory variables (81 ways) and the method of determining the objective variable (3 ways).

上記実施の形態では、腐食速度推定装置1を用いて鉄塔の腐食速度を推定していたが、本発明はこれに限られない。例えば、電線、架空地線、鋼管、柱上変電器、車体、各種筐体等の腐食速度の推定に腐食速度推定装置1を用いてもよい。 In the above embodiment, the corrosion rate of the steel tower is estimated by using the corrosion rate estimation device 1, but the present invention is not limited to this. For example, the corrosion rate estimation device 1 may be used to estimate the corrosion rate of electric wires, overhead ground wires, steel pipes, pillar transformers, vehicle bodies, various housings, and the like.

上記実施の形態では、腐食速度マップを腐食速度推定装置1の表示部12にて表示していたが、本発明はこれに限られない。例えば、腐食速度推定装置1の外部にある端末等に出力してもよい。 In the above embodiment, the corrosion rate map is displayed on the display unit 12 of the corrosion rate estimation device 1, but the present invention is not limited to this. For example, it may be output to a terminal or the like outside the corrosion rate estimation device 1.

上記実施の形態においては、年間飛来海塩粒子量、濡れ時間、乾湿サイクル等のデータは腐食速度推定装置1の記憶部14に記憶されていたが、本発明はこれに限定されない。例えば、各種データは、その全部又は一部が通信ネットワークを介して外部のサーバやコンピュータ等に記憶されていてもよい。 In the above embodiment, data such as the annual amount of sea salt particles flying, the wetting time, and the dry / wet cycle are stored in the storage unit 14 of the corrosion rate estimation device 1, but the present invention is not limited thereto. For example, all or part of the various data may be stored in an external server, computer, or the like via a communication network.

上記実施の形態では、重回帰分析を実行することで、腐食速度の推定式を導出していたが、本発明はこれに限られない。入力データと出力データとの関係を規定しており、入力層にデータを供給して出力層から出力データを得ることができる手法であれば、いかなる手法で腐食速度の推定式を導出してもよい。 In the above embodiment, the estimation formula of the corrosion rate has been derived by performing the multiple regression analysis, but the present invention is not limited to this. Any method that defines the relationship between input data and output data and can supply data to the input layer and obtain output data from the output layer can be used to derive the corrosion rate estimation formula. good.

例えば、入力層と、中間層と、出力層と、を備える階層型ニューラルネットワークを用いた予測モデルを用いて腐食速度の推定式を導出してもよい。より具体的には、入力用ニューロンと中間ニューロンとの間、中間ニューロンと出力用ニューロンとの間に、それぞれニューロン間の結合状態を設定する重み係数を設定し、誤差逆伝播法による学習演算処理を実行することで、入力データ及び出力データからなる学習データを用いて重み係数を最適化してもよい。 For example, an estimation formula of the corrosion rate may be derived using a prediction model using a hierarchical neural network including an input layer, an intermediate layer, and an output layer. More specifically, a weighting coefficient that sets the connection state between neurons is set between the input neuron and the intermediate neuron, and between the intermediate neuron and the output neuron, and learning calculation processing by the error back propagation method is performed. By executing the above, the weighting coefficient may be optimized using the training data consisting of the input data and the output data.

上記実施の形態では、通信ネットワークとしてインターネットを用いていたが、本発明はこれに限られない。例えば、通信ネットワークは、LAN(Local Area Network)や専用線等を用いて実現してもよい。 In the above embodiment, the Internet is used as the communication network, but the present invention is not limited to this. For example, the communication network may be realized by using a LAN (Local Area Network), a dedicated line, or the like.

上記実施の形態においては、腐食速度推定装置1は、記憶部14に記憶されたプログラムに基づいて動作していたが、本発明はこれに限定されない。例えば、プログラムにより実現された機能的な構成をハードウェアにより実現してもよい。 In the above embodiment, the corrosion rate estimation device 1 operates based on the program stored in the storage unit 14, but the present invention is not limited thereto. For example, the functional configuration realized by the program may be realized by hardware.

上記実施の形態では、腐食速度推定装置1は、例えば汎用コンピュータであったが、本発明はこれに限られない。例えば、腐食速度推定装置1は、専用のシステムで実現してもよく、クラウド上に設けられたコンピュータであってもよい。また、腐食速度推定装置1が実行する処理は、例えば、上述の物理的な構成を備える装置が、記憶部14に記憶されたプログラムを実行することによって実現されていたが、本発明は、プログラムとして実現されてもよく、そのプログラムが記録された記憶媒体として実現されてもよい。 In the above embodiment, the corrosion rate estimation device 1 is, for example, a general-purpose computer, but the present invention is not limited to this. For example, the corrosion rate estimation device 1 may be realized by a dedicated system or may be a computer provided on the cloud. Further, the process executed by the corrosion rate estimation device 1 is realized by, for example, an apparatus having the above-mentioned physical configuration executing a program stored in the storage unit 14, but the present invention has a program. It may be realized as a storage medium in which the program is recorded.

また、上述の処理動作を実行させるためのプログラムを、フレキシブルディスク、CD-ROM(Compact Disk Read-Only Memory)、DVD(Digital Versatile Disk)、MO(Magneto-Optical Disk)等のコンピュータにより読み取り可能な記録媒体に格納して配布し、そのプログラムをコンピュータにインストールすることにより、上述の処理動作を実行する装置を構成してもよい。 Further, the program for executing the above-mentioned processing operation can be read by a computer such as a flexible disk, a CD-ROM (Compact Disk Read-Only Memory), a DVD (Digital Versatile Disk), or an MO (Magneto-Optical Disk). A device that executes the above-mentioned processing operation may be configured by storing it in a recording medium, distributing it, and installing the program on a computer.

上記実施の形態は例示であり、本発明はこれらに限定されるものではなく、特許請求の範囲に記載した発明の趣旨を逸脱しない範囲でさまざまな実施の形態が可能である。各実施の形態や変形例で記載した構成要素は自由に組み合わせることが可能である。また、特許請求の範囲に記載した発明と均等な発明も本発明に含まれる。 The above embodiments are examples, and the present invention is not limited thereto, and various embodiments are possible without departing from the spirit of the invention described in the claims. The components described in each embodiment and modification can be freely combined. The present invention also includes inventions equivalent to those described in the claims.

1 腐食速度推定装置
11 操作部
12 表示部
13 通信部
14 記憶部
141 年間飛来海塩粒子量データテーブル
142 濡れ時間データテーブル
143 乾湿サイクルデータテーブル
144 腐食実績データテーブル
15 制御部
151 推定式作成部
152 腐食速度推定部
153 腐食速度マップ作成部

1 Corrosion rate estimation device 11 Operation unit 12 Display unit 13 Communication unit 14 Storage unit 141 Annual arrival sea salt particle amount data table 142 Wet time data table 143 Wet and dry cycle data table 144 Corrosion record data table 15 Control unit 151 Estimate formula creation unit 152 Corrosion rate estimation unit 153 Corrosion rate map creation unit

Claims (9)

金属材料の腐食速度を目的変数とし、前記金属材料の腐食に影響を与える年間飛来海塩粒子量、濡れ時間及び乾湿サイクルを含む腐食環境因子を説明変数とする重回帰分析を行うことで、前記金属材料の腐食速度の推定式を作成する推定式作成手段と、
前記推定式に年間飛来海塩粒子量、濡れ時間及び乾湿サイクルを含む腐食環境因子を入力することで、前記金属材料の腐食速度を推定する腐食速度推定手段と、
を備え、
乾湿サイクルは、所定時間以上の湿潤状態と所定時間以上の乾燥状態とが交互に継続するサイクル数である、
腐食速度推定装置。
By performing multiple regression analysis with the corrosion rate of the metal material as the objective variable and the corrosive environment factors including the annual amount of sea salt particles flying in the year, the wetting time, and the dry-wet cycle that affect the corrosion of the metal material as the explanatory variables. An estimation formula creation method for creating an estimation formula for the corrosion rate of metal materials,
A corrosion rate estimation means for estimating the corrosion rate of the metal material by inputting a corrosion environment factor including an annual amount of sea salt particles, a wetting time, and a dry / wet cycle into the estimation formula.
Equipped with
The dry-wet cycle is the number of cycles in which a wet state for a predetermined time or longer and a dry state for a predetermined time or longer are alternately continued.
Corrosion rate estimation device.
前記推定式は、以下の式に示すとおりである請求項1に記載の腐食速度推定装置。
Y=α+β+β+β
ただし、Yは金属材料の腐食速度、Xは年間飛来海塩粒子量、Xは濡れ時間、Xは乾湿サイクル、αは定数、β~βは係数である。
The corrosion rate estimation device according to claim 1, wherein the estimation formula is as shown in the following formula.
Y = α + β 1 X 1 + β 2 X 2 + β 3 X 3
However, Y is the corrosion rate of the metal material, X 1 is the amount of sea salt particles flying annually, X 2 is the wetting time, X 3 is the dry / wet cycle, α is a constant, and β 1 to β 3 are coefficients.
前記推定式作成手段における年間飛来海塩粒子量、濡れ時間及び乾湿サイクルは、いずれも気象解析コードから得られたシミュレーション値である、
請求項1又は2に記載の腐食速度推定装置。
The annual amount of sea salt particles, the wetting time, and the dry / wet cycle in the estimation formula creating means are all simulation values obtained from the meteorological analysis code.
The corrosion rate estimation device according to claim 1 or 2.
前記金属材料は、鉄塔を構成する鋼鉄部材であり、
前記腐食速度は、前記鋼鉄部材に施された亜鉛めっき厚さの一年間当たりの減少量である、
請求項1から3のいずれか1項に記載の腐食速度推定装置。
The metal material is a steel member constituting a steel tower.
The corrosion rate is the amount of decrease in the galvanized thickness applied to the steel member per year.
The corrosion rate estimation device according to any one of claims 1 to 3.
前記金属材料は、送電線を構成するアルミ素線であり、
前記腐食速度は、前記アルミ素線の引張強度の一年間当たりの減少量である、
請求項1から3のいずれか1項に記載の腐食速度推定装置。
The metal material is an aluminum wire constituting a power transmission line.
The corrosion rate is the amount of decrease in the tensile strength of the aluminum wire per year.
The corrosion rate estimation device according to any one of claims 1 to 3.
前記濡れ時間及び前記乾湿サイクルは、複数の定義のパターンを用意しておき、重回帰分析を行った場合に最も解析結果が良好な定義のパターンを選択する、
請求項1から5のいずれか1項に記載の腐食速度推定装置。
For the wetting time and the dry-wet cycle, a plurality of defined patterns are prepared, and the defined pattern with the best analysis result is selected when the multiple regression analysis is performed.
The corrosion rate estimation device according to any one of claims 1 to 5.
推定された前記金属材料の腐食速度に基づいて、所定の地域における金属材料の腐食速度を示す腐食速度マップを作成する腐食速度マップ作成手段を備える、
請求項1から6のいずれか1項に記載の腐食速度推定装置。
A means for creating a corrosion rate map is provided, which creates a corrosion rate map showing the corrosion rate of the metal material in a predetermined area based on the estimated corrosion rate of the metal material.
The corrosion rate estimation device according to any one of claims 1 to 6.
金属材料の腐食速度を目的変数とし、前記金属材料の腐食に影響を与える年間飛来海塩粒子量、濡れ時間及び乾湿サイクルを含む腐食環境因子を説明変数とする重回帰分析を行うことで、前記金属材料の腐食速度の推定式を作成する推定式作成ステップと、
前記推定式に年間飛来海塩粒子量、濡れ時間及び乾湿サイクルを含む腐食環境因子を入力することで、前記金属材料の腐食速度を推定する腐食速度推定ステップと、
を含み、
乾湿サイクルは、所定時間以上の湿潤状態と所定時間以上の乾燥状態とが交互に継続するサイクル数である、
腐食速度推定方法。
By performing multiple regression analysis with the corrosion rate of the metal material as the objective variable and the corrosive environment factors including the annual amount of sea salt particles flying in the year, the wetting time, and the dry-wet cycle that affect the corrosion of the metal material as the explanatory variables. An estimation formula creation step to create an estimation formula for the corrosion rate of metal materials, and
A corrosion rate estimation step for estimating the corrosion rate of the metal material by inputting a corrosion environment factor including an annual amount of sea salt particles, a wetting time, and a dry / wet cycle into the estimation formula.
Including
The dry-wet cycle is the number of cycles in which a wet state for a predetermined time or longer and a dry state for a predetermined time or longer are alternately continued.
Corrosion rate estimation method.
コンピュータを、
金属材料の腐食速度を目的変数とし、前記金属材料の腐食に影響を与える年間飛来海塩粒子量、濡れ時間及び乾湿サイクルを含む腐食環境因子を説明変数とする重回帰分析を行うことで、前記金属材料の腐食速度の推定式を作成する推定式作成手段、
前記推定式に年間飛来海塩粒子量、濡れ時間及び乾湿サイクルを含む腐食環境因子を入力することで、前記金属材料の腐食速度を推定する腐食速度推定手段、
として機能させるためのプログラムであって、
乾湿サイクルは、所定時間以上の湿潤状態と所定時間以上の乾燥状態とが交互に継続するサイクル数である、
プログラム。
Computer,
By performing multiple regression analysis with the corrosion rate of the metal material as the objective variable and the corrosive environment factors including the annual amount of sea salt particles flying in the year, the wetting time, and the dry-wet cycle that affect the corrosion of the metal material as the explanatory variables. An estimation formula creation method for creating an estimation formula for the corrosion rate of metal materials,
A corrosion rate estimation means for estimating the corrosion rate of the metal material by inputting the corrosion environment factors including the annual amount of sea salt particles, the wetting time and the dry / wet cycle into the estimation formula.
It is a program to function as
The dry-wet cycle is the number of cycles in which a wet state for a predetermined time or longer and a dry state for a predetermined time or longer are alternately continued.
program.
JP2019044872A 2019-03-12 2019-03-12 Corrosion rate estimation device, corrosion rate estimation method and program Active JP7084344B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2019044872A JP7084344B2 (en) 2019-03-12 2019-03-12 Corrosion rate estimation device, corrosion rate estimation method and program

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2019044872A JP7084344B2 (en) 2019-03-12 2019-03-12 Corrosion rate estimation device, corrosion rate estimation method and program

Publications (2)

Publication Number Publication Date
JP2020148541A JP2020148541A (en) 2020-09-17
JP7084344B2 true JP7084344B2 (en) 2022-06-14

Family

ID=72431903

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2019044872A Active JP7084344B2 (en) 2019-03-12 2019-03-12 Corrosion rate estimation device, corrosion rate estimation method and program

Country Status (1)

Country Link
JP (1) JP7084344B2 (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2771144C2 (en) * 2020-07-13 2022-04-27 Федеральное Государственное Бюджетное Учреждение Научно-Исследовательский Институт Проблем Хранения Федерального Агентства По Государственным Резервам (Фгбу Ниипх Росрезерва) Method for predicting shelf life of non-ferrous metals in closed rooms
WO2023233465A1 (en) * 2022-05-30 2023-12-07 三菱電機株式会社 Corrosive environment diagnostic system
KR102616805B1 (en) * 2023-05-08 2023-12-21 주식회사 이케이전력 IoT-based corrosion pre-detection switchgear preventive maintenance system

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3066162U (en) 1999-07-28 2000-02-18 スガ試験機株式会社 Spray test equipment
JP2005337838A (en) 2004-05-26 2005-12-08 Chugoku Electric Power Co Inc:The Corrosion deterioration evaluation device and method, computer program and recording medium
JP2006208346A (en) 2004-02-24 2006-08-10 Jfe Steel Kk Method and device for predicting life of steel product, and computer program
JP2008224405A (en) 2007-03-13 2008-09-25 Tokyo Electric Power Co Inc:The Corrosion rate evaluating method
JP2013011498A (en) 2011-06-29 2013-01-17 Kita Nippon Electric Cable Co Ltd Method for determining lifetime of overhead transmission lines

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3066162U (en) 1999-07-28 2000-02-18 スガ試験機株式会社 Spray test equipment
JP2006208346A (en) 2004-02-24 2006-08-10 Jfe Steel Kk Method and device for predicting life of steel product, and computer program
JP2005337838A (en) 2004-05-26 2005-12-08 Chugoku Electric Power Co Inc:The Corrosion deterioration evaluation device and method, computer program and recording medium
JP2008224405A (en) 2007-03-13 2008-09-25 Tokyo Electric Power Co Inc:The Corrosion rate evaluating method
JP2013011498A (en) 2011-06-29 2013-01-17 Kita Nippon Electric Cable Co Ltd Method for determining lifetime of overhead transmission lines

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
國武義高,架空送電設備の電線腐食推定におけるNuWiCC-STの適用について,エネルギア総研レビュー,No.33,日本,2013年,p.11~14,https://www.energia.co.jp/eneso/kankoubutsu/review/no33/index.html

Also Published As

Publication number Publication date
JP2020148541A (en) 2020-09-17

Similar Documents

Publication Publication Date Title
JP7084344B2 (en) Corrosion rate estimation device, corrosion rate estimation method and program
JP3909057B2 (en) Corrosion prediction method for weathering steel
Costantini et al. A hierarchical procedure for the combination of forecasts
JP5066955B2 (en) Corrosion rate evaluation method
CN109726865A (en) Customer charge probability density prediction technique, device and storage medium based on EMD-QRF
WO2019135361A1 (en) Method for predicting corrosion amount of metal material, method for selecting metal material, and device for predicting corrosion amount of metal material
Vera et al. On the prediction of atmospheric corrosion of metals and alloys in Chile using artificial neural networks
CN113138158B (en) Corrosion-resistant life prediction method and device for metal material and electronic equipment
JP2003329573A (en) Method for evaluating corrosion resistance of metallic material, method for predicting corrosion life thereof, metallic material, method for designing thereof, and method for manufacturing thereof
Wi et al. Customer baseline load (CBL) calculation using exponential smoothing model with weather adjustment
JP5034021B2 (en) Wave estimation method, system and program
JP2003215024A (en) Method for predicting amount of corrosion of metallic material due to galvanic corrosion, life predicting method, metallic material, structure designing method, and method for manufacturing metallic material
Jiang et al. Application of accelerated cyclic test with synthetic ocean water salt-deposition process to the evaluation on corrosion characteristics of weathering steel
CN112686693A (en) Method, system, equipment and storage medium for predicting marginal electricity price of electric power spot market
CN114729875A (en) Method for plotting corrosion amount of metal material, method for selecting metal material, and apparatus for plotting corrosion amount of metal material
JP2003169415A (en) Aged deterioration management system utilizing map information for plated structure
Panchenko et al. Evaluating the reliability of predictions of first-year corrosion losses of structural metals calculated using dose-response functions for territories with different categories of atmospheric corrosion aggressiveness
CN114564487B (en) Meteorological raster data updating method combining forecast prediction
Li et al. Evaluating the corrosion resistance of marine steels under different exposure environments via machine learning
CN110728466B (en) Method for determining target accessory demand of new product and computer equipment
CN114066004A (en) GDP (graphics device performance) current prediction method based on electric power big data
Hameed Smoothing techniques for time series forecasting
Park et al. An appraisal of composite forecasting methods
JP5352530B2 (en) Method for estimating the corrosion state of steel
Mwenda et al. Analyzing the impact of historical data length in non-seasonal ARIMA models forecasting

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20201112

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20211130

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20220121

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20220322

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20220418

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20220510

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20220602

R150 Certificate of patent or registration of utility model

Ref document number: 7084344

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150