(57)【要約】  金属材料の腐食量マッピング方法は、金属材料の使用期間と、複数の環境パラメータと、マップ上における環境パラメータの位置座標と、マップの地勢データと、金属材料の腐食量とを含む腐食量データを用いて腐食量予測マップを作成する方法であり、複数の環境パラメータと、マップの地勢データと、環境パラメータの位置座標とから任意のメッシュの間隔で環境マップを作成する環境マップ作成ステップと、予測要求点入力ステップと、類似度算出ステップと、腐食量予測ステップと、メッシュの予測要求点における腐食量の予測結果をマップ上に色付けすることにより、腐食量予測マップを作成する腐食量予測マップ作成ステップと、を含む。
2021100341.tif 000002
【特許請求の範囲】
【請求項１】
金属材料の使用期間と、前記使用期間内の前記金属材料の使用環境を示し、かつ前記金属材料が使用される地域を示すマップ上で既知の複数の環境パラメータと、前記マップ上における前記環境パラメータの位置座標と、前記マップの地勢データと、前記使用期間における前記金属材料の腐食量と、を含む腐食量データを用いて金属材料の腐食量を予測し、腐食量予測マップを作成する金属材料の腐食量マッピング方法であって、
前記複数の環境パラメータと、前記マップの地勢データと、前記マップ上における前記環境パラメータの位置座標とから、任意のメッシュの間隔で、前記環境パラメータごとの環境マップを作成する環境マップ作成ステップと、
前記腐食量予測マップを作成する前記金属材料の使用期間と、前記腐食量データにおける複数の環境パラメータと、を含む予測要求点を入力する予測要求点入力ステップと、
前記腐食量データにおける複数の環境パラメータと、前記予測要求点における複数の環境パラメータとの類似度を算出する類似度算出ステップと、
前記類似度を考慮して、前記腐食量データにおける複数の環境パラメータを、潜在変数に次元圧縮する次元圧縮ステップと、
前記潜在変数および前記類似度を用いて構築した予測式に基づいて、前記メッシュの予測要求点における前記金属材料の腐食量を予測する腐食量予測ステップと、
前記メッシュの予測要求点における腐食量の予測結果を、マップ上に色付けすることにより、腐食量予測マップを作成する腐食量予測マップ作成ステップと、
を含むことを特徴とする金属材料の腐食量マッピング方法。
【請求項２】
前記環境マップ作成ステップは、
前記複数の環境パラメータを前記マップ上の位置座標に入力する第一ステップと、
前記マップの地勢データに基づいて、前記環境マップを海抜０ｍの環境マップへと補正する第二ステップと、
前記海抜０ｍの環境マップにおいて、前記複数の環境パラメータ間の環境パラメータを補完する第三ステップと、
前記マップの地勢データに基づいて、環境パラメータを補完した海抜０ｍの環境マップを、元の海抜の環境マップへと補正する第四ステップと、
を含むことを特徴とする請求項１に記載の金属材料の腐食量マッピング方法。
【請求項３】
前記複数の環境パラメータは、飛来塩分量を含み、
前記第三ステップは、海岸線に近いメッシュの飛来塩分量を補完する際に、予め設定した上限値を超えないように飛来塩分量を補完することを特徴とする請求項２に記載の金属材料の腐食量マッピング方法。
【請求項４】
前記腐食量予測ステップは、
前記金属材料の所定期間の腐食量を示す第一のパラメータを予測する初期腐食量予測ステップと、
前記金属材料の腐食速度の減衰を示す第二のパラメータを予測する減衰予測ステップと、
前記金属材料の使用期間と、前記第一のパラメータと、前記第二のパラメータと、に基づいて前記金属材料の前記所定期間よりも長い期間の腐食量を予測する長期腐食量予測ステップと、
を含むことを特徴とする請求項１から請求項３のいずれか一項に記載の金属材料の腐食量マッピング方法。
【請求項５】
前記初期腐食量予測ステップは、前記潜在変数および前記類似度を用いて構築した予測式に基づいて、前記予測要求点の環境パラメータ下における、前記金属材料の所定期間の腐食量を予測することを特徴とする請求項４に記載の金属材料の腐食量マッピング方法。
【請求項６】
前記減衰予測ステップは、前記腐食量予測マップを作成する前記金属材料の使用期間と、前記腐食量データにおける複数の環境パラメータと、前記類似度とに基づいて、前記第二のパラメータを予測することを特徴とする請求項４に記載の金属材料の腐食量マッピング方法。
【請求項７】
前記複数の環境パラメータは、温度と、相対湿度、絶対湿度、濡れ時間および降雨量のうちの少なくとも一つと、飛来塩分量、ＳＯＸ濃度およびＮＯＸ濃度のうちの少なくとも一つと、を含むことを特徴とする請求項１から請求項６のいずれか一項に記載の金属材料の腐食量マッピング方法。
【請求項８】
前記金属材料は、鉄鋼材料であることを特徴とする請求項１から請求項４のいずれか一項に記載の金属材料の腐食量マッピング方法。
【請求項９】
前記腐食量予測マップは、前記腐食量予測マップ作成ステップで予測された腐食量の値に応じて、前記メッシュごとに色付けされたマップであることを特徴とする請求項１から請求項８のいずれか一項に記載の金属材料の腐食量マッピング方法。
【請求項１０】
請求項１から請求項９のいずれか一項に記載の金属材料の腐食量マッピング方法を用いて、使用環境に応じた金属材料を選定することを特徴とする金属材料の選定方法。
【請求項１１】
金属材料の使用期間と、前記使用期間内の前記金属材料の使用環境を示し、かつ前記金属材料が使用される地域を示すマップ上で既知の複数の環境パラメータと、前記マップ上における前記環境パラメータの位置座標と、前記マップの地勢データと、前記使用期間における前記金属材料の腐食量と、を含む腐食量データが保存されるデータベースと、
前記複数の環境パラメータと、前記マップの地勢データと、前記マップ上における前記環境パラメータの位置座標とから、任意のメッシュの間隔で、前記環境パラメータごとの環境マップを作成する環境マップ作成部と、
腐食量予測マップを作成する前記金属材料の使用期間と、前記腐食量データにおける複数の環境パラメータと、を含む予測要求点が入力される入力部と、
前記腐食量データにおける複数の環境パラメータと、前記予測要求点における複数の環境パラメータとの類似度を算出する類似度算出部と、
前記類似度を考慮して、前記腐食量データにおける複数の環境パラメータを、潜在変数に次元圧縮する次元圧縮部と、
前記潜在変数および前記類似度を用いて構築した予測式に基づいて、前記メッシュの予測要求点における前記金属材料の腐食量を予測する腐食量予測部と、
前記メッシュの予測要求点における腐食量の予測結果を、マップ上に色付けすることにより、前記腐食量予測マップを作成する腐食量予測マップ作成部と、
を備えることを特徴とする金属材料の腐食量マッピング装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、金属材料の腐食量マッピング方法、金属材料の選定方法および金属材料の腐食量マッピング装置に関する。
【背景技術】
【０００２】
非特許文献１に示すように、従来から大気腐食環境における金属材料の腐食量は、経験式として以下の式（１）で表されることが知られている。【０００３】
【数１】
2021100341.tif 000003
【０００４】
ここで、上記式（１）において、Ｙは金属材料の腐食量、Ｘは金属材料の使用期間、Ａは金属材料の初期一年間の腐食量を示すパラメータ、Ｂは腐食により形成されるさび層の効果による腐食速度の減衰を示すパラメータ、である。これらのパラメータＡ、Ｂの値は、金属材料の種類や大気腐食環境により変化する。そのため、現在は、長期腐食量を予測する際に、対象となる大気腐食環境において金属材料を複数期間暴露し、腐食量の経時変化を上記式（１）で外挿する手法が多く用いられている。【０００５】
ところで、金属材料の腐食量は、金属が持つ耐食性能と大気腐食環境因子、例えば温度、相対湿度、濡れ時間、降雨量、飛来塩分量およびＳＯ２濃度等が複雑に作用することにより決まる。そこで、以下に示すように、上記の環境因子を用いて定式化することにより、金属材料の腐食量を予測する技術が提案されている。【０００６】
例えば非特許文献２では、炭素鋼、亜鉛、銅およびアルミニウムについて、腐食量の対数を、温度、相対湿度および飛来塩の対数で重回帰した項と、温度、相対湿度およびＳＯ２量の対数で重回帰した項との和により算出している。【０００７】
また、特許文献１では、年間濡れ時間、年平均風速、飛来塩分量、硫黄酸化物量、腐食反応の活性化エネルギーおよび温度をパラメータとして、腐食指標Ｚを下記式（２）のように表し、この腐食指標Ｚの二次関数を長期腐食量の算出に用いている。【０００８】
【数２】
2021100341.tif 000004
【０００９】
ここで、上記式（２）において、ＴＯＷは年間濡れ時間（ｈ）、Ｗは年平均風速（ｍ／ｓ）、Ｃは飛来塩分量（ｍｄｄ）、Ｓは硫黄酸化物量（ｍｄｄ）、Ｅαは腐食反応の活性化エネルギー（Ｊ／ｍｏｌ）、Ｒは気体定数（Ｊ／（Ｋ／ｍｏｌ））、Ｔは年平均気温（Ｋ）、α，κ，δ，εは定数、である。また、前記した「ｍｄｄ」は、単位日数および単位面積あたりに捕集されたＮａＣｌの量であり、「ｍｇ  ＮａＣｌ・ｄｍ−２・ｄａｙ−１」の略である。【００１０】
また、特許文献２では、飛来塩分量Ｃを予測する技術が提案されている。また、特許文献３では、経験式として知られる上記式（１）のパラメータＡを、温度、相対湿度、飛来塩分量および濡れ確率をパラメータとして下記式（３）により算出し、パラメータＢを実験室的試験によりパラメータＡの関数として算出する腐食予測技術が提案されている。【００１１】
【数３】
2021100341.tif 000005
【００１２】
ここで、上記式（３）において、Ｔは温度（℃）、Ｈは相対湿度（％）、Ｓａは飛来塩分量（ｍｇ／ｄｍ２／ｄａｙ（＝ｍｄｄ））、Ｐｗ（Ｔ，Ｈ）は濡れ確率、α，β，γは鋼種に応じて設定される係数、である。【００１３】
また、特許文献４では、屋外の大気腐食環境において、経験式として知られる上記式（１）のパラメータＡを、温度、濡れ時間および飛来塩分量をパラメータとして下記式（４）により算出し、パラメータＢを０．３〜０．６の範囲に設定する腐食予測技術が提案されている。【００１４】
【数４】
2021100341.tif 000006
【００１５】
ここで、上記式（４）において、Ｔは温度（℃）、ＴＯＷは濡れ時間（ｈ）、Ｓａは飛来塩分量（ｍｇ／ｄｍ２／ｄａｙ（＝ｍｄｄ））、α，β，γは係数、である。【００１６】
また、特許文献５では、大気環境における鋼材の板厚減少量を予測する際に、経験式として知られる上記式（１）のパラメータＡを下記式（５）により、パラメータＢを下記式（６）により算出する技術が提案されている。【００１７】
【数５】
2021100341.tif 000007
【数６】
2021100341.tif 000008
【００１８】
ここで、上記式（５）および上記式（６）において、ＣＲ０は、環境因子をパラメータとする鋼材製造直後の初期腐食速度を表す関数であり、ＣＲ１は、環境因子をパラメータとする鋼材製造から１年後の腐食速度を表す関数である。また、ここでの環境因子とは、年平均温度（℃）、年平均湿度（％）、年平均風速（ｍ／ｓｅｃ）、飛来塩分量（ｍｇ／ｄｍ２／ｄａｙ（＝ｍｄｄ））、硫黄酸化物量（ｍｇ／ｄｍ２／ｄａｙ（＝ｍｄｄ））、のことを示している。【００１９】
また、特許文献６では、金属材料の腐食速度を予測する際に、腐食速度を目的変数としその腐食速度に影響を与える環境因子と地形因子を説明変数とする重回帰分析を行っている。そして、この重回帰分析を行うにあたり、少なくとも説明変数の一つとして相対湿度０％〜１００％に応じて重み付けした仮想濡れ時間を含めている。また、この仮想濡れ時間を、変化する相対湿度に応じて異なる重み係数を変化する相対湿度に対応した時間に乗算して得られた乗算値を総和して求め、測定した金属材料の腐食速度に基づき重回帰分析法により腐食速度推定式を作成する手法が提案されている。【００２０】
また、特許文献７では、腐食量の予測結果をマッピング（マップ化）する技術として、ステップワイズ法を用いて最適化した重回帰分析によって金属材料の腐食速度を予測し、当該金属材料の腐食速度を、クラスタリングによってマッピングする技術が提案されている。
【先行技術文献】
【特許文献】
【００２１】
【特許文献１】特許第３９０９０５７号公報
【特許文献２】特許第４１４３０１８号公報
【特許文献３】特許第４７０６２５４号公報
【特許文献４】特許第５８９５５２２号公報
【特許文献５】特許第５０６６１６０号公報
【特許文献６】特許第５０６６９５５号公報
【特許文献７】特許第５６８４５５２号公報
【非特許文献】
【００２２】
【非特許文献１】「耐候性鋼材の橋梁への適用に関する共同研究報告書（ＸＶＩＩＩ）」、建設省土木研究所、（社）鋼材倶楽部、（社）日本橋梁建設協会、平成５年３月
【非特許文献２】ISO  9223：1992  “Corrosion  of  metals  and  alloys  -  Corrosivity  of  atmospheres  -  Classification、  determination  and  estimation”
【発明の概要】
【発明が解決しようとする課題】
【００２３】
ここで、特許文献１〜５では、腐食量および腐食速度と各環境パラメータとの関係を評価することにより、説明変数として採用する環境パラメータを選定して定式化している。しかしながら、腐食量、腐食速度および各環境パラメータは、複雑に相関関係を持っている。例えば、腐食量と温度との関係は非線形であり、飛来塩分量とＳＯ２濃度とは疑似相関を持っている。このような関係性の中で、特許文献１〜５のように定式化した場合、高い精度の予測は望めない。【００２４】
また、特許文献６は、濡れ時間に着目して重み付けをすることを特徴としている。このように重み付けをすることにより精度は上がるものの、多数ある環境パラメータの中で濡れ時間にのみ重みを付けても有効な精度向上は望めない。また、特許文献６の手法では、現在データを保有している期間でのみ腐食予測が可能であり、長期腐食予測をすることはできない。【００２５】
また、特許文献７は、クラスタリングによってデータを分類し、クラスごとに多数ある環境パラメータから複数の環境パラメータを選定し、重回帰を繰り返し行って最良の式とすることにより予測精度向上を図っている。しかしながら、腐食量および腐食速度の関係は、重回帰で得られる各環境パラメータとの１次式で簡単に表現することができない。従って、特許文献７の手法では大きな精度向上は望めない。【００２６】
また、特許文献７で開示された方法では、予測地点の環境パラメータ、例えば日本国内の場合は１ｋｍ四方のメッシュごとの環境パラメータを、公に提供されているデータ（例えば農研機構が提供する「農業気象データ」等）等から取得している。しかしながら、このようなデータでは、飛来塩分量やＳＯ２量等の全ての環境パラメータは網羅されていない。また、特許文献７で開示された方法では、不足する環境パラメータの補完方法についても明示されていない。そのため、特許文献７で開示された方法を用いて、既知の環境パラメータのみでクラスタリングおよび腐食量の予測を行ったとしても、精度の低い腐食量予測マップしか作成できないものと考えられる。【００２７】
本発明は、上記に鑑みてなされたものであって、大気腐食環境において、金属材料の腐食量を精度良くマッピングすることができる金属材料の腐食量マッピング方法、金属材料の選定方法および金属材料の腐食量マッピング装置を提供することを目的とする。
【課題を解決するための手段】
【００２８】
上述した課題を解決し、目的を達成するために、本発明に係る金属材料の腐食量マッピング方法は、金属材料の使用期間と、前記使用期間内の前記金属材料の使用環境を示し、かつ前記金属材料が使用される地域を示すマップ上で既知の複数の環境パラメータと、前記マップ上における前記環境パラメータの位置座標と、前記マップの地勢データと、前記使用期間における前記金属材料の腐食量と、を含む腐食量データを用いて金属材料の腐食量を予測し、腐食量予測マップを作成する金属材料の腐食量マッピング方法であって、前記複数の環境パラメータと、前記マップの地勢データと、前記マップ上における前記環境パラメータの位置座標とから、任意のメッシュの間隔で、前記環境パラメータごとの環境マップを作成する環境マップ作成ステップと、前記腐食量予測マップを作成する前記金属材料の使用期間と、前記腐食量データにおける複数の環境パラメータと、を含む予測要求点を入力する予測要求点入力ステップと、前記腐食量データにおける複数の環境パラメータと、前記予測要求点における複数の環境パラメータとの類似度を算出する類似度算出ステップと、前記類似度を考慮して、前記腐食量データにおける複数の環境パラメータを、潜在変数に次元圧縮する次元圧縮ステップと、前記潜在変数および前記類似度を用いて構築した予測式に基づいて、前記メッシュの予測要求点における前記金属材料の腐食量を予測する腐食量予測ステップと、前記メッシュの予測要求点における腐食量の予測結果を、マップ上に色付けすることにより、腐食量予測マップを作成する腐食量予測マップ作成ステップと、を含むことを特徴とする。【００２９】
また、本発明に係る金属材料の腐食量マッピング方法は、上記発明において、前記環境マップ作成ステップが、前記複数の環境パラメータを前記マップ上の位置座標に入力する第一ステップと、前記マップの地勢データに基づいて、前記環境マップを海抜０ｍの環境マップへと補正する第二ステップと、前記海抜０ｍの環境マップにおいて、前記複数の環境パラメータ間の環境パラメータを補完する第三ステップと、前記マップの地勢データに基づいて、環境パラメータを補完した海抜０ｍの環境マップを、元の海抜の環境マップへと補正する第四ステップと、を含むことを特徴とする。【００３０】
また、本発明に係る金属材料の腐食量マッピング方法は、上記発明において、前記複数の環境パラメータが、飛来塩分量を含み、前記第三ステップは、海岸線に近いメッシュの飛来塩分量を補完する際に、予め設定した上限値を超えないように飛来塩分量を補完することを特徴とする。【００３１】
また、本発明に係る金属材料の腐食量マッピング方法は、上記発明において、前記腐食量予測ステップが、前記金属材料の所定期間の腐食量を示す第一のパラメータを予測する初期腐食量予測ステップと、前記金属材料の腐食速度の減衰を示す第二のパラメータを予測する減衰予測ステップと、前記金属材料の使用期間と、前記第一のパラメータと、前記第二のパラメータと、に基づいて前記金属材料の前記所定期間よりも長い期間の腐食量を予測する長期腐食量予測ステップと、を含むことを特徴とする。【００３２】
また、本発明に係る金属材料の腐食量マッピング方法は、上記発明において、前記初期腐食量予測ステップが、前記潜在変数および前記類似度を用いて構築した予測式に基づいて、前記予測要求点の環境パラメータ下における、前記金属材料の所定期間の腐食量を予測することを特徴とする。【００３３】
また、本発明に係る金属材料の腐食量マッピング方法は、上記発明において、前記減衰予測ステップが、前記腐食量予測マップを作成する前記金属材料の使用期間と、前記腐食量データにおける複数の環境パラメータと、前記類似度とに基づいて、前記第二のパラメータを予測することを特徴とする。【００３４】
また、本発明に係る金属材料の腐食量マッピング方法は、上記発明において、前記複数の環境パラメータが、温度と、相対湿度、絶対湿度、濡れ時間および降雨量のうちの少なくとも一つと、飛来塩分量、ＳＯＸ濃度およびＮＯＸ濃度のうちの少なくとも一つと、を含むことを特徴とする。【００３５】
また、本発明に係る金属材料の腐食量マッピング方法は、上記発明において、前記金属材料が、鉄鋼材料であることを特徴とする。【００３６】
また、本発明に係る金属材料の腐食量マッピング方法は、上記発明において、前記腐食量予測マップが、前記腐食量予測マップ作成ステップで予測された腐食量の値に応じて、前記メッシュごとに色付けされたマップであることを特徴とする。【００３７】
上述した課題を解決し、目的を達成するために、本発明に係る金属材料の選定方法は、前記した金属材料の腐食量マッピング方法を用いて、使用環境に応じた金属材料を選定することを特徴とする。【００３８】
上述した課題を解決し、目的を達成するために、本発明に係る金属材料の腐食量マッピング装置は、金属材料の使用期間と、前記使用期間内の前記金属材料の使用環境を示し、かつ前記金属材料が使用される地域を示すマップ上で既知の複数の環境パラメータと、前記マップ上における前記環境パラメータの位置座標と、前記マップの地勢データと、前記使用期間における前記金属材料の腐食量と、を含む腐食量データが保存されるデータベースと、前記複数の環境パラメータと、前記マップの地勢データと、前記マップ上における前記環境パラメータの位置座標とから、任意のメッシュの間隔で、前記環境パラメータごとの環境マップを作成する環境マップ作成部と、腐食量予測マップを作成する前記金属材料の使用期間と、前記腐食量データにおける複数の環境パラメータと、を含む予測要求点が入力される入力部と、前記腐食量データにおける複数の環境パラメータと、前記予測要求点における複数の環境パラメータとの類似度を算出する類似度算出部と、前記類似度を考慮して、前記腐食量データにおける複数の環境パラメータを、潜在変数に次元圧縮する次元圧縮部と、前記潜在変数および前記類似度を用いて構築した予測式に基づいて、前記メッシュの予測要求点における前記金属材料の腐食量を予測する腐食量予測部と、前記メッシュの予測要求点における腐食量の予測結果を、マップ上に色付けすることにより、腐食量予測マップを作成する腐食量予測マップ作成部と、を備えることを特徴とする。【００３９】
本発明は以上の知見に基づいてなされたものであり、以下を要旨とするものである。
［１］
金属材料の使用期間と、前記使用期間内の前記金属材料の使用環境を示し、かつ前記金属材料が使用される地域を示すマップ上で既知の複数の環境パラメータと、前記マップ上における前記環境パラメータの位置座標と、前記マップの地勢データと、前記使用期間における前記金属材料の腐食量と、を含む腐食量データを用いて金属材料の腐食量を予測し、腐食量予測マップを作成する金属材料の腐食量マッピング方法であって、
前記複数の環境パラメータと、前記マップの地勢データと、前記マップ上における前記環境パラメータの位置座標とから、任意のメッシュの間隔で、前記環境パラメータごとの環境マップを作成する環境マップ作成ステップと、
前記腐食量予測マップを作成する前記金属材料の使用期間と、前記腐食量データにおける複数の環境パラメータと、を含む予測要求点を入力する予測要求点入力ステップと、
前記腐食量データにおける複数の環境パラメータと、前記予測要求点における複数の環境パラメータとの類似度を算出する類似度算出ステップと、
前記類似度を考慮して、前記腐食量データにおける複数の環境パラメータを、潜在変数に次元圧縮する次元圧縮ステップと、
前記潜在変数および前記類似度を用いて構築した予測式に基づいて、前記メッシュの予測要求点における前記金属材料の腐食量を予測する腐食量予測ステップと、
前記メッシュの予測要求点における腐食量の予測結果を、マップ上に色付けすることにより、腐食量予測マップを作成する腐食量予測マップ作成ステップと、
を含むことを特徴とする金属材料の腐食量マッピング方法。
［２］
前記環境マップ作成ステップは、
前記複数の環境パラメータを前記マップ上の位置座標に入力する第一ステップと、
前記マップの地勢データに基づいて、前記環境マップを海抜０ｍの環境マップへと補正する第二ステップと、
前記海抜０ｍの環境マップにおいて、前記複数の環境パラメータ間の環境パラメータを補完する第三ステップと、
前記マップの地勢データに基づいて、環境パラメータを補完した海抜０ｍの環境マップを、元の海抜の環境マップへと補正する第四ステップと、
を含むことを特徴とする［１］に記載の金属材料の腐食量マッピング方法。
［３］
前記複数の環境パラメータは、飛来塩分量を含み、
前記第三ステップは、海岸線に近いメッシュの飛来塩分量を補完する際に、予め設定した上限値を超えないように飛来塩分量を補完することを特徴とする［２］に記載の金属材料の腐食量マッピング方法。
［４］
前記腐食量予測ステップは、
前記金属材料の所定期間の腐食量を示す第一のパラメータを予測する初期腐食量予測ステップと、
前記金属材料の腐食速度の減衰を示す第二のパラメータを予測する減衰予測ステップと、
前記金属材料の使用期間と、前記第一のパラメータと、前記第二のパラメータと、に基づいて前記金属材料の前記所定期間よりも長い期間の腐食量を予測する長期腐食量予測ステップと、
を含むことを特徴とする［１］から［３］のいずれか一項に記載の金属材料の腐食量マッピング方法。
［５］
前記初期腐食量予測ステップは、前記潜在変数および前記類似度を用いて構築した予測式に基づいて、前記予測要求点の環境パラメータ下における、前記金属材料の所定期間の腐食量を予測することを特徴とする［４］に記載の金属材料の腐食量マッピング方法。
［６］
前記減衰予測ステップは、前記腐食量予測マップを作成する前記金属材料の使用期間と、前記腐食量データにおける複数の環境パラメータと、前記類似度とに基づいて、前記第二のパラメータを予測することを特徴とする［４］に記載の金属材料の腐食量マッピング方法。
［７］
前記複数の環境パラメータは、温度と、相対湿度、絶対湿度、濡れ時間および降雨量のうちの少なくとも一つと、飛来塩分量、ＳＯＸ濃度およびＮＯＸ濃度のうちの少なくとも一つと、を含むことを特徴とする［１］から［６］のいずれか一項に記載の金属材料の腐食量マッピング方法。
［８］
前記金属材料は、鉄鋼材料であることを特徴とする［１］から［４］のいずれか一項に記載の金属材料の腐食量マッピング方法。
［９］
前記腐食量予測マップは、前記腐食量予測マップ作成ステップで予測された腐食量の値に応じて、前記メッシュごとに色付けされたマップであることを特徴とする［１］から［８］のいずれか一項に記載の金属材料の腐食量マッピング方法。
［１０］
［１］から［９］のいずれか一項に記載の金属材料の腐食量マッピング方法を用いて、使用環境に応じた金属材料を選定することを特徴とする金属材料の選定方法。
［１１］
金属材料の使用期間と、前記使用期間内の前記金属材料の使用環境を示し、かつ前記金属材料が使用される地域を示すマップ上で既知の複数の環境パラメータと、前記マップ上における前記環境パラメータの位置座標と、前記マップの地勢データと、前記使用期間における前記金属材料の腐食量と、を含む腐食量データが保存されるデータベースと、
前記複数の環境パラメータと、前記マップの地勢データと、前記マップ上における前記環境パラメータの位置座標とから、任意のメッシュの間隔で、前記環境パラメータごとの環境マップを作成する環境マップ作成部と、
腐食量予測マップを作成する前記金属材料の使用期間と、前記腐食量データにおける複数の環境パラメータと、を含む予測要求点が入力される入力部と、
前記腐食量データにおける複数の環境パラメータと、前記予測要求点における複数の環境パラメータとの類似度を算出する類似度算出部と、
前記類似度を考慮して、前記腐食量データにおける複数の環境パラメータを、潜在変数に次元圧縮する次元圧縮部と、
前記潜在変数および前記類似度を用いて構築した予測式に基づいて、前記メッシュの予測要求点における前記金属材料の腐食量を予測する腐食量予測部と、
前記メッシュの予測要求点における腐食量の予測結果を、マップ上に色付けすることにより、前記腐食量予測マップを作成する腐食量予測マップ作成部と、
を備えることを特徴とする金属材料の腐食量マッピング装置。
【発明の効果】
【００４０】
本発明によれば、大気腐食環境において、金属材料の腐食量を精度良くマッピングすることができ、使用環境に応じた耐食性を有する最適な金属材料を選定することが可能となる。
【図面の簡単な説明】
【００４１】
【図１】図１は、大気腐食環境において、温度（年平均）と腐食量（年間）との関係を示すグラフである。
【図２】図２は、大気腐食環境において、ＳＯ２濃度（年平均）と飛来塩分量（年平均）との擬似相関関係を示すグラフである。
【図３】図３は、本発明の実施形態に係る金属材料の腐食量マッピング装置の構成を示すブロック図である。
【図４】図４は、本発明の実施形態に係る金属材料の腐食量マッピング方法の流れを示すフローチャートである。
【図５】図５は、従来技術に係る金属材料の腐食量マッピング方法によって得た、金属材料の一年後の腐食量予測値を示す腐食量予測マップである。
【図６】図６は、本発明の実施形態に係る金属材料の腐食量マッピング方法によって得た、金属材料の五十年後の腐食量予測値を示す腐食量予測マップである。
【発明を実施するための形態】
【００４２】
以下、本発明の実施形態に係る金属材料の腐食量マッピング方法、金属材料の選定方法および金属材料の腐食量マッピング装置について、図面を参照しながら説明する。なお、本発明は以下の実施形態に限定されるものではない。
【００４３】
（環境パラメータの補完）
大気中における金属材料と腐食量と各種の環境パラメータとの間には複雑な相関関係があり、前述の非特許文献２の国際規格ＩＳＯ９２２３の中でも、年平均の温度、相対湿度、飛来塩分量，ＳＯＸ濃度に基づいて年間の腐食量が定式化されている。しかしながら、これらの環境パラメータを、例えばマップ上で１ｋｍ、２ｋｍ程度の狭いメッシュ間隔で入手することは困難であるため、入手可能な環境パラメータの位置座標に基づいて、各環境パラメータの位置座標間の環境パラメータを補完する必要がある。
【００４４】
その際、例えば温度は、経度、緯度のみならず標高によっても変化する。また、飛来塩分の飛来の仕方も、山や丘等の障害物の影響によって変化する。そのため、このような地勢の影響を考慮して環境パラメータを補完する必要がある。
【００４５】
但し、地勢データの反映および環境パラメータの位置座標間の補完を一度の演算で行うことは困難である。そのため、本実施形態では、まず既知の複数の環境パラメータがプロットされた複数の環境マップを、海抜０ｍ換算の複数の環境マップへと補正する。続いて、環境パラメータの位置座標間の補完した海抜０ｍの環境マップをそれぞれ作成する。続いて、実際の地勢データに合わせて補完後の環境パラメータを補正することにより、元の海抜の環境マップをそれぞれ作成する。なお、前記した「環境マップ」とは、環境パラメータがプロットされたマップを示している。以下、環境パラメータごとの具体的な補完方法について説明する。
【００４６】
（温度の補完）
温度の補完では、まず入手可能な温度データを、マップ上にプロットする。続いて、「標高が１００ｍ下がると気温が０．６℃上がる」という法則（気温減率）に基づいて、地勢データを用いて、各プロット点の温度の値を、海抜０ｍに換算した値へと補正する。続いて、任意のメッシュ間隔でデータ間の温度を補完する。なお、メッシュ間隔は、マップ上で既知の環境パラメータの数や、計算時のマシンスペック等に応じて決定すればよい。
【００４７】
データの補完方法は、線形補完や多項式補完等の様々な手法があるが、環境を予測した際のフルクロスバリデーションにより、本実施形態のユーザが求める精度に達していれば、どのような補完方法を用いてもよい。これは、温度以外の環境パラメータの補完においても同様である。そして、「標高が１００ｍ上がると気温が０．６℃下がる」という法則に基づいて、地勢データを用いて、海抜０ｍの環境マップの各プロット点の温度の値を、元の海抜に換算した値へと補正することにより、元の海抜の温度の環境マップ（温度マップ）を作成する。
【００４８】
（絶対湿度の補完）
絶対湿度は、環境パラメータのうちの温度および相対湿度から算出することができる。ここで、大気腐食環境は、「温度が変化しても絶対湿度はほぼ一定」という特徴を有している。そのため、上記の特徴を利用して、標高等の地勢データを反映させることなく、温度および相対湿度から各メッシュ点の絶対湿度を算出し、マップ上の絶対湿度のデータ間を補完することにより、絶対湿度の環境マップ（絶対湿度マップ）を作成する。
【００４９】
（相対湿度の補完）
相対湿度の補完では、上記の「温度が変化しても絶対湿度はほぼ一定」という大気腐食環境の特徴を利用して、温度および絶対湿度の環境マップの各メッシュ点で相対湿度を算出することにより、相対湿度の環境マップ（相対湿度マップ）を作成する。
【００５０】
（濡れ時間の補完）
濡れ時間に関して、国際規格では「相対湿度が８０％以上の時間」と定められている。そのため、上記の手法を用いた時間ごと、日ごとの相対湿度の経時変化マップから、各メッシュ点における相対湿度８０％以上の時間を積算して年の濡れ時間を算出することにより、濡れ時間の環境マップ（濡れ時間マップ）を作成する。
【００５１】
（降雨量の補完）
降雨量は、地勢の影響を受けないため、標高等の地勢データを反映させることなく、マップ上における降雨量のデータ間を補完することにより、降雨量の環境マップ（降雨量マップ）を作成する。
【００５２】
（飛来塩分量、ＳＯＸ濃度およびＮＯＸ濃度の補完）
飛来塩分量、ＳＯＸ濃度およびＮＯＸ濃度の補完では、経度、緯度、標高のユークリッド距離に基づいてデータ間を補完することにより、飛来塩分量、ＳＯＸ濃度およびＮＯＸ濃度の環境マップ（飛来塩分量マップ、ＳＯＸ濃度マップおよびＮＯＸ濃度マップ）を作成する。
【００５３】
ここで、飛来塩分量は、具体的にはモデル式「ｙ＝ａｘ−ｂ」により算出する。このモデル式において、ｘ：離岸距離（ｋｍ）、ｙ：飛来塩分量（ｍｄｄ）、ａ，ｂ：係数である。また、離岸距離ｘは、海岸線の形状データの各点からの距離の最小値である。なお、飛来塩分量は、上記のモデル式以外にも、Ｃｏｌｅモデルやメソ気象モデル等の既存の腐食予測式を用いて算出してもよい。
【００５４】
また、飛来塩分量の補完に際しては、メッシュ点が海岸線と近い場合に飛来塩分量が異常に高くなることを避けるために、飛来塩分量の上限を設けることが好ましい。飛来塩分量の上限としては、例えば「１．０ｍｄｄ＝６２．３ｍｍｄ」等の値を設定することができる。
【００５５】
（腐食量の予測）
金属材料の腐食量の予測では、上記のように作成した環境マップのメッシュごとに腐食量を予測する。本発明者らは、各大気腐食環境における金属材料の初期一年間の腐食量を示すパラメータＡ（第一のパラメータ）と、さび層による腐食速度の減衰を示すパラメータＢ（第二のパラメータ）とをそれぞれ分けて予測し、かつ予測に際し、各環境パラメータに重み付けをしたデータから予測することにより、腐食量の予測精度が向上することを見出した。
【００５６】
金属材料の腐食速度は、一般的に経時的に減衰していく。これは、金属材料表面に生成する腐食生成物（例えばさび層）による保護効果によるものである。また、この保護効果は、周囲の環境や金属材料の種類により大きく異なる。このように、金属材料の腐食速度は、様々な環境や金属材料が持つ耐食性等の、非常に多くの要素が絡み合う。そのため、任意の環境および期間における金属材料の腐食量を、各環境パラメータと腐食量との関係の原理原則から精度良く予測することは非常に困難であり、例えばデータベースに蓄積された腐食量および各環境パラメータのデータ群から統計的に予測することが現実的であり、精度向上にも繋がる。
【００５７】
一方、一般的に、蓄積された腐食量と各環境パラメータのデータ群の中に、長期間のデータは多くない。例えば腐食量を予測する任意の期間が数十年という長期間となる場合、期間を変数として予測に取り込んで直接腐食量を予測すると、任意の長期間の腐食量を、期間の離れたデータから予測することになるため、精度が低下してしまう。そこで、本発明では、金属材料の初期一年間の腐食量を示すパラメータＡと、さび層による腐食速度の減衰を示すパラメータＢとをそれぞれ分けて予測することにより、精度の向上を図った。
【００５８】
ここで、任意の期間における金属材料の腐食量に関わる主な環境パラメータとしては、例えば温度、相対湿度、絶対湿度、濡れ時間、降雨量、飛来塩分量、ＳＯＸ濃度およびＮＯＸ濃度等が挙げられる。これらの環境パラメータの中には、例えば図１に示した温度と腐食量との関係のように、非線形の関係を持つものが存在する。また、例えば図２に示した飛来塩分量とＳＯ２濃度との関係のように、環境パラメータ同士が多重共線性を示すものも存在する。金属材料の腐食に影響を及ぼす環境パラメータが複数存在することに加えて、この二点が、任意の環境および期間における金属材料の腐食量を精度良く予測することを困難にしている要因となっている。
【００５９】
腐食量に対して環境パラメータが非線形の関係を持つものに関しては、予測をしたい任意の環境および期間との類似性により各サンプルに重み付けを行い、かつ局所的に重回帰分析を行うことで予測精度を改善することができる。なお、前記した「サンプル」とは、データベースに保存された腐食量および各環境パラメータのデータ群（後記する腐食量データ）のことを意味している。
【００６０】
また、環境パラメータ同士が多重共線性を持つことに関しては、各環境パラメータに対して独立性を有するパラメータに次元圧縮し、新たなパラメータを作り出すことにより解決することができる。そして、これらを同時に実現する手法の一つとして、以下の参考文献１に示す「局所重み付き部分最小二乗法（ＬＷ−ＰＬＳ：Locally  weighted  partial  least  squares）」がある。
【００６１】
参考文献１：金尚弘、岡島亮太、加納学、長谷部伸治、「高精度な局所PLSモデル構築のためのサンプル選択」、第54回自動制御連合会講演会、54(2011)、p.1594
【００６２】
本発明では、各大気腐食環境における金属材料の初期一年間の腐食量を示すパラメータＡと、さび層による腐食速度の減衰を示すパラメータＢとを、予測要求点との類似性（類似度）をサンプルごとに求め、その類似度を利用して重み付けを行い、局所回帰することにより予測する手法を用いる。また、各環境パラメータを次元圧縮することにより新たなパラメータを作り出し（潜在変数の導出）、局所回帰の説明変数とする。またその際、類似度を利用した重み付けを用いた潜在変数と目的変数の内積が最大となるように潜在変数を決定し、局所的な重回帰を行う。以下、本発明の具体的な実施形態について、図面を参照しながら説明する。
【００６３】
（腐食量マッピング装置）
本発明の実施形態に係る金属材料の腐食量マッピング装置の構成について、図３を参照しながら説明する。腐食量マッピング装置１は、入力部１０と、データベース２０と、演算部３０と、表示部４０と、を備えている。
【００６４】
入力部１０は、例えばキーボード、マウスポインタ、テンキー等の入力装置によって実現される。後記するように、演算部３０には、この入力部１０を介して、後記する予測要求点が入力される。
【００６５】
データベース２０には、金属材料の腐食量の実績値である腐食量データが保存されている。腐食量データには、金属材料（例えば鉄鋼材料）の使用期間と、当該使用期間における金属材料の腐食量と、当該使用期間内の金属材料の使用環境を示し、かつ金属材料が使用される地域を示すマップ上で既知の複数の環境パラメータと、マップ上における環境パラメータの位置座標と、マップの地勢データと、が含まれる。
【００６６】
前記した「複数の環境パラメータ」には、温度（気温）と、相対湿度、絶対湿度、濡れ時間および降雨量のうちの少なくとも一つと、飛来塩分量、ＳＯＸ濃度およびＮＯＸ濃度のうちの少なくとも一つと、が含まれる。また、これらの環境パラメータは、例えば年平均のデータである。また、データベース２０には、鋼種ごとの腐食量データが保存されている。
【００６７】
演算部３０は、具体的には、ＣＰＵ（Central  Processing  Unit）、ＤＳＰ（Digital  Signal  Processor）、ＦＰＧＡ（Field-Programmable  Gate  Array）等からなるプロセッサと、ＲＡＭ（Random  Access  Memory）、ＲＯＭ（Read  Only  Memory）等からなるメモリ（主記憶部）と、により実現される。
【００６８】
演算部３０は、例えば図示しない記憶部に格納されたプログラムを主記憶部の作業領域にロードして実行し、プログラムの実行を通じて各構成部等を制御することにより、所定の目的に合致した機能を実現する。演算部３０は、プログラムの実行を通じて、マップ作成部３１、類似度算出部３２、次元圧縮部３３および腐食量予測部３４として機能する。なお、マップ作成部３１、類似度算出部３２、次元圧縮部３３および腐食量予測部３４の詳細については後記する。
【００６９】
表示部４０は、例えばＬＣＤディスプレイ、ＣＲＴディスプレイ等の表示装置によって実現され、演算部３０から入力される表示信号をもとに、金属材料の腐食量の予測結果として腐食量予測マップを表示する。なお、前記した「腐食量予測マップ」とは、腐食量の予測値が例えば色によって表示されたマップを示している（図５および図６参照）。
【００７０】
（腐食量マッピング方法）
本発明の実施形態に係る金属材料の腐食量マッピング方法について、図４を参照しながら説明する。腐食量マッピング方法では、環境マップ作成ステップと、予測要求点入力ステップと、第一の類似度算出ステップと、第一の次元圧縮ステップと、初期腐食量予測ステップ（第一のパラメータ予測ステップ）と、第二の類似度算出ステップと、第二の次元圧縮ステップと、減衰予測ステップ（第二のパラメータ予測ステップ）と、長期腐食量予測ステップと、腐食量予測マップ作成ステップと、を行う。なお、本実施形態に係る腐食量マッピング方法は、金属材料の中でも特に耐候性鋼等の鋼材に適用することにより、腐食量をより精度良く予測してマッピングすることが可能である。
【００７１】
環境マップ作成ステップでは、マップ作成部３１によって、入手可能なマップ上における環境パラメータから、任意のメッシュ間隔でデータ間を補完して、環境パラメータごとの環境マップを作成する（ステップＳ１）。
【００７２】
環境マップ作成ステップでは、マップ上で既知の複数の環境パラメータと、マップの地勢データと、マップ上における環境パラメータの位置座標とから、任意のメッシュの間隔で、環境パラメータごとの環境マップを作成する。環境マップ作成ステップでは、より具体的には以下の四つのステップを行う。
【００７３】
まず、既知の複数の環境データを、前記したマップ上の位置座標に入力する（第一ステップ）。続いて、マップの地勢データに基づいて、環境マップを海抜０ｍの環境マップへと補正する（第二ステップ）。続いて、海抜０ｍの環境マップにおいて、既知の環境パラメータ間の環境パラメータを補完する（第三ステップ）。続いて、マップの地勢データに基づいて、環境パラメータを補完した海抜０ｍの環境マップを、元の海抜の環境マップへと補正する（第四ステップ）。
【００７４】
ここで、前記した第三ステップでは、前記したように、海岸線に近いメッシュ点の飛来塩分量を補完する際に、予め設定した上限値（例えば１．０ｍｄｄ）を超えないように飛来塩分量を補完することが好ましい。飛来塩分量に上限値を設けることにより、海岸線に近いメッシュ点の飛来塩分量が異常に高くなることを避けることができる。
【００７５】
予測要求点入力ステップでは、入力部１０を介して、メッシュ点の環境パラメータである予測要求点が演算部３０に入力される（ステップＳ２）。この予測要求点は、腐食量を予測したい金属材料の使用期間（腐食量予測マップを作成する金属材料の使用期間）と、この使用期間内の金属材料の使用環境を示す、年平均の複数の環境パラメータ（腐食量データにおける複数の環境パラメータ）と、を含んでいる。
【００７６】
続いて、第一の類似度算出ステップでは、類似度算出部３２によって、データベース２０に保存されている、使用期間が一年である金属材料の腐食量データにおける複数の環境パラメータと、予測要求点における複数の環境パラメータとの類似度を算出する（ステップＳ３）。本ステップにおいて、類似度算出部３２は、例えば後記する式（８）によって上記類似度を算出する。なお、本ステップの具体例については後記する実施例で説明する。
【００７７】
続いて、第一の次元圧縮ステップでは、次元圧縮部３３によって、第一の類似度算出ステップで算出された類似度を考慮して、腐食量データにおける複数の環境パラメータ（説明変数）を、潜在変数に次元圧縮する（ステップＳ４）。本ステップにおいて、次元圧縮部３３は、例えば後記する式（７）によって上記潜在変数を算出する。なお、本ステップの具体例については後記する実施例で説明する。
【００７８】
続いて、初期腐食量予測ステップでは、腐食量予測部３４によって、第一の次元圧縮ステップで算出された潜在変数、および第一の類似度算出ステップで算出された類似度を用いて構築した予測式に基づいて、予測要求点の環境パラメータ下における、金属材料の初期一年間の腐食量を予測する（ステップＳ５）。この金属材料の初期一年間の腐食量は、上記式（１）のパラメータＡ（第一のパラメータ）のことを意味している。本ステップにおいて、腐食量予測部３４は、例えば後記する式（１０）に示す予測式を構築し、この予測式に基づいて金属材料の初期一年間の腐食量を予測する。なお、本ステップの具体例については後記する実施例で説明する。
【００７９】
続いて、第二の類似度算出ステップでは、類似度算出部３２によって、腐食量を予測したい金属材料の使用期間を考慮して、金属材料の腐食量データにおける複数の環境パラメータと、予測要求点における複数の環境パラメータとの類似度を算出する（ステップＳ６）。本ステップにおいて、類似度算出部３２は、例えば後記する式（１３）によって上記類似度を算出する。なお、本ステップの具体例については後記する実施例で説明する。
【００８０】
続いて、第二の次元圧縮ステップでは、次元圧縮部３３によって、第二の類似度算出ステップで算出された類似度および腐食量を予測したい金属材料の使用期間を考慮して、腐食量データにおける複数の環境パラメータ（説明変数）を、潜在変数に次元圧縮する（ステップＳ７）。本ステップにおいて、次元圧縮部３３は、例えば後記する式（１２）によって上記潜在変数を算出する。なお、本ステップの具体例については後記する実施例で説明する。
【００８１】
続いて、減衰予測ステップでは、腐食量予測部３４によって、第二の次元圧縮ステップで算出された潜在変数、および第二の類似度算出ステップで算出された類似度を用いて構築した予測式に基づいて、金属材料の腐食速度の減衰を示すパラメータを予測する（ステップＳ８）。この金属材料の腐食速度の減衰を示すパラメータは、上記式（１）のパラメータＢ（第二のパラメータ）のことを意味している。本ステップにおいて、腐食量予測部３４は、例えば後記する式（１６）に示す予測式を構築し、この予測式に基づいて金属材料の腐食速度の減衰を示すパラメータを予測する。なお、本ステップの具体例については後記する実施例で説明する。
【００８２】
続いて、長期腐食量予測ステップでは、腐食量予測部３４によって、金属材料の使用期間と、初期腐食量予測ステップで算出されたパラメータＡと、減衰予測ステップで算出されたパラメータＢと、に基づいて金属材料の一年よりも長い期間の腐食量（長期腐食量）を予測する（ステップＳ９）。本ステップにおいて、腐食量予測部３４は、上記式（１）に基づいて上記長期腐食量を予測する。なお、本ステップの具体例については後記する実施例で説明する。
【００８３】
続いて、腐食量予測マップ作成ステップでは、マップ作成部３１によって、メッシュの予測要求点における腐食量の予測結果を、マップ上に色付けすることにより、腐食量予測マップを作成する（ステップＳ１０）。この腐食量予測マップは、腐食量予測マップ作成ステップで予測された腐食量の値に応じて、前記メッシュごとに色付けされたマップである（例えば後記する図５および図６参照）。
【００８４】
以上のように、本実施形態に係る腐食量マッピング装置１を用いた金属材料の腐食量マッピング方法によれば、大気腐食環境において、金属材料の長期腐食予測を精度高く行うことができる。また、金属材料の腐食量を精度良くマッピングすることができ、使用環境に応じた耐食性を有する最適な金属材料を選定することが可能となる。
【００８５】
また、従来技術、例えば特許文献１〜５のように説明変数として採用する環境パラメータを選定して定式化してしまうと、ある地域では予測精度が高いが、ある地域では予測精度が低い、というように腐食量を予測しようとする地域によって予測精度に差が生じてしまう。一方、本実施形態に係る金属材料の腐食予測方法によれば、データベース２０内に、腐食量を予測しようとする地域のものと類似する腐食量データがあれば予測可能であるため、どのような地域であっても金属材料の腐食量を精度高く予測することができる。
【実施例】
【００８６】
（実施例１）
以下、実施例を挙げて本発明をより詳細に説明する。本実施例では、ベトナムにおいて、軒下環境における鋼材の年間腐食量（初期一年間の腐食量）を予測し、腐食量予測マップを作成した。ここで、本実施例の内容は、前記した腐食量マッピング方法のうち、環境マップ作成ステップ、予測要求点入力ステップ、第一の類似度算出ステップ、第一の次元圧縮ステップ、初期腐食量予測ステップおよび腐食量予測マップ作成ステップに相当する。【００８７】
本実施例では、データベースに保存されている、世界各地の５５地域における鋼材の年間腐食量と年平均の各環境パラメータのデータ群（腐食量データ）を用いて、年間腐食量の予測を行った。本実施例では、環境パラメータとして、温度（℃）、相対湿度（％）、飛来塩分量（ｍｍｄ（Ｃｌ−換算））、ＳＯ２濃度（ｍｍｄ（ＳＯ２換算））、の四つを使用した。ここで、前記した「ｍｍｄ」は、単位日数および単位面積あたりに捕集されたＣｌ−またはＳＯ２の量であり、「ｍｇ・ｍ−２・ｄａｙ−１」の略である。【００８８】
環境パラメータごとの環境マップは、前記した方法により作成した。今回の実施例ではデータ間の補完方法として線形補完を利用して計算を行った（図４の「環境マップ作成ステップ」）。得られた環境マップのメッシュ点での腐食予測において、前記したように、説明変数同士が相関を持つ場合、多重共線性によって予測精度が悪化することが知られている。前記したＬＷ−ＰＬＳでは、説明変数として使用する各環境パラメータが相関を持っているため、次元圧縮によって相関を排除し、新たなパラメータ（潜在変数）を作成する。ここで、ＬＷ−ＰＬＳは、前記した参考文献１に示された手順に従って計算している。四つの環境パラメータを潜在変数（パラメータｔ）に圧縮するための式は、下記式（７）のように示すことができる。【００８９】
【数７】
2021100341.tif 000009
【００９０】
ここで、上記式（７）において、Ｔは温度（℃）、ＲＨは相対湿度（％）、Ｃｌは飛来塩分量（ｍｇ／ｍ２／ｄａｙ（＝ｍｍｄ）（Ｃｌ−換算））、ＳＯ２はＳＯ２濃度（ｍｇ／ｍ２／ｄａｙ（＝ｍｍｄ）（ＳＯ２換算））、ｗ１〜ｗ４は係数、である。なお、上記式（７）では、一例として四つの環境パラメータのみを用いているが、実際には、腐食量の予測をしようとする地域において腐食に関係することが予想される環境パラメータについては全て含めることが好ましい。【００９１】
本実施例では、まず腐食量を予測したい鋼材の使用期間と、この使用期間内の鋼材の使用環境を示す年平均の複数の環境パラメータと、を含む予測要求点を、腐食量マッピング装置の演算部に入力した（図４の「予測点入力ステップ」）。【００９２】
続いて、腐食量を求めたい予測要求点の環境パラメータと、腐食量の予測のために参照する腐食量データの環境パラメータｉの類似度ωｉを、下記式（８）に示すようなユークリッド距離を用いて算出した（図４の「第一の類似度算出ステップ」）。なお、各環境パラメータによってデータの広がりが異なるため、数値は正規化を実施した。【００９３】
【数８】
2021100341.tif 000010
【００９４】
ここで、ωｉは局所化パラメータ、右下の添え字がｑの環境パラメータは腐食量を予測したい予測要求点の環境パラメータ、右下の添え字がｉの環境パラメータはデータベースから参照した腐食量データの環境パラメータ、σは下記式（９）の標準偏差である。また、φは調整パラメータであり、例えばφ＝１を基準として適宜に調整して値は決定される。【００９５】
【数９】
2021100341.tif 000011
【００９６】
続いて、上記式（８）で算出した類似度ωｉと、腐食量データの環境パラメータおよび腐食量から、参考文献１（２．１章）に示された手順に従って、類似度ωｉを重みとして係った潜在変数と目的変数（腐食量）の内積が最大となるように、上記式（７）の係数ｗ１〜ｗ４を決定した。そして、決定した係数ｗ１〜ｗ４を用いて、上記式（７）によって各環境パラメータの潜在変数を算出した（図４の「第一の次元圧縮ステップ」）。【００９７】
続いて、下記式（１０）に示す腐食量の予測式を局所回帰によって構築し、下記式（１０）に基づいて、予測要求点の環境パラメータ下における鋼材の年間腐食量（初期一年間の腐食量）を予測した（図４の「初期腐食量予測ステップ」）。【００９８】
【数１０】
2021100341.tif 000012
【００９９】
ここで、Ｙは腐食量の予測値、αは係数（回帰係数）である。なお、上記式（１０）では記載を省略したが、上記式（１０）に定数項を含ませたり、あるいは複数の潜在変数を用いることも可能である。【０１００】
このように、本実施例では、腐食量を予測したい予測要求点が入力されるたびに、予測要求点と各腐食量データとの類似度の算出、潜在変数の係数の算出、予測式の構築を行う。得られたメッシュ点ごとの年間腐食量の予測値を、例えば図５に示すような腐食量予測マップとしてマップ化し（図４の「腐食量予測マップ作成ステップ」）、表示部４０によって表示する。【０１０１】
（実施例２）
本実施例では、ベトナムにおいて、軒下環境における鋼材の長期（五十年）腐食量を予測し、腐食量予測マップを作成した。ここで、本実施例の内容は、前記した腐食量マッピング方法の全ステップに相当する。【０１０２】
本実施例では、軒下環境を対象として、データベースに保存されている、世界５５地域の１年間、５５地域の３年間、３９地域の５年間、３８地域の７年間、３８地域の９年間の腐食量と年平均の各環境パラメータのデータ群（腐食量データ）を用いて、実施例１と同様の手法で、鋼材について、初期一年間の腐食量を示すパラメータＡ（上記式（１）参照）を算出した（図４の環境マップ作成ステップ、予測要求点入力ステップ、第一の類似度算出ステップ、第一の次元圧縮ステップおよび初期腐食量予測ステップ）。【０１０３】
続いて、７年間の長期腐食試験結果に重み付けをし、さび層による腐食速度の減衰を示すパラメータＢを算出し、上記式（１）により鋼材の五十年後の腐食量を予測した。五十年後の腐食量の予測では、まず上記式（１）を下記式（１１）に示すように変形し、パラメータＢおよび使用期間Ｘを用いて左辺の初期一年間の腐食量の対数から五十年後の腐食量の対数の差を算出した。【０１０４】
【数１１】
2021100341.tif 000013
【０１０５】
上記式（１１）の左辺の予測にはＬＷ−ＰＬＳを用いる。具体的には、四つの環境パラメータと使用期間Ｘの対数の積を取ることにより、時間因子を考慮した新たな環境パラメータを作成し、これらの新たな環境パラメータを、下記式（１２）に示すように潜在変数（パラメータｕ）に圧縮する。【０１０６】
【数１２】
2021100341.tif 000014
【０１０７】
ここで、上記式（１２）において、Ｔは温度（℃）、ＲＨは相対湿度（％）、Ｃｌは飛来塩分量（ｍｇ／ｍ２／ｄａｙ（＝ｍｍｄ）（Ｃｌ−換算））、ＳＯ２はＳＯ２濃度（ｍｇ／ｍ２／ｄａｙ（＝ｍｍｄ）（ＳＯ２換算））、ｖ１〜ｖ４は係数、である。【０１０８】
続いて、腐食量を求めたい予測要求点の環境パラメータと、腐食量の予測のために参照する腐食量データの環境パラメータｉの類似度ωｉを、下記式（１３）に示すようなユークリッド距離を用いて算出した（図４の「第二の類似度算出ステップ」）。なお、各環境パラメータによってデータの広がりが異なるため、数値は正規化を実施した。【０１０９】
【数１３】
2021100341.tif 000015
【０１１０】
ここで、ωｉは局所化パラメータ、右下の添え字がｑの環境パラメータは腐食量を予測したい予測要求点の環境パラメータ、右下の添え字がｉの環境パラメータはデータベースから参照した腐食量データの環境パラメータ、σは下記式（１４）の標準偏差である。また、φは調整パラメータであり、例えばφ＝１を基準として適宜に調整して値は決定される。【０１１１】
【数１４】
2021100341.tif 000016
【０１１２】
続いて、上記式（１３）で算出した類似度ωｉと、腐食量データの環境パラメータおよび腐食量から、参考文献１（２．１章）に示された手順に従って、類似度ωｉを重みとして係った潜在変数と目的変数（腐食量）の内積が最大となるように、上記式（１２）の係数ｖ１〜ｖ４を決定した。そして、決定した係数ｖ１〜ｖ４を用いて、上記式（７）によって各環境パラメータの潜在変数を算出した（図４の「第二の次元圧縮ステップ」）。【０１１３】
続いて、下記式（１５）に示す腐食量の予測式を局所回帰によって構築し、予測したい環境パラメータにおける腐食量の対数の差を算出した。【０１１４】
【数１５】
2021100341.tif 000017
【０１１５】
ここで、上記式（１５）において、βは係数（回帰係数）である。なお、上記式（１５）では記載を省略したが、上記式（１５）に定数項を含ませたり、あるいは複数の潜在変数を用いることも可能である。【０１１６】
このように、本発明例では、腐食量を予測したい予測要求点が入力されるたびに、予測要求点と各腐食量データとの類似度の算出、潜在変数の係数の算出、予測式の構築を行う。【０１１７】
続いて、上記式（１２）および上記式（１５）に基づいて、上記式（１１）のパラメータＢを、下記式（１６）によって算出した（図４の「減衰予測ステップ」）。【０１１８】
【数１６】
2021100341.tif 000018
【０１１９】
続いて、算出したパラメータＡおよびパラメータＢを用いて、上記式（１）により９年後の腐食量Ｙを算出した（図４の長期腐食量予測ステップ）。得られたメッシュ点ごとの長期腐食量の予測値を、例えば図６に示すような腐食量予測マップとしてマップ化し（図４の「腐食量予測マップ作成ステップ」）、表示部４０によって表示する。【０１２０】
以上、本発明に係る金属材料の腐食量マッピング方法、金属材料の選定方法および金属材料の腐食量マッピング装置について、発明を実施するための形態および実施例により具体的に説明したが、本発明の趣旨はこれらの記載に限定されるものではなく、請求の範囲の記載に基づいて広く解釈されなければならない。また、これらの記載に基づいて種々変更、改変等したものも本発明の趣旨に含まれることはいうまでもない。【０１２１】
ここで、前記した実施形態では、金属材料の初期一年間の腐食量（パラメータＡ）と、金属材料の腐食速度の減衰を示すパラメータ（パラメータＢ）とをそれぞれ分けて予測し、初期一年間の腐食量を基準として長期腐食量を予測していたが、長期腐食量を予測する際の基準は初期一年間の腐食量には限定されない。【０１２２】
すなわち、初期腐食量予測ステップにおいて、予め定める任意の所定期間（初期期間）における金属材料の腐食量を予測し、長期腐食量予測ステップにおいて、前記した所定期間の腐食量を基準として長期腐食量を予測してもよい。【０１２３】
例えば初期腐食量として１．５年の腐食量がＡ’として与えられた場合、そこからＸ年後の腐食量の予測式は、上記式（１）を拡張して、下記式（１７）のように記述可能と考えられる。【０１２４】
【数１７】
2021100341.tif 000019
【０１２５】
これを一般化すれば、ある初期期間Ｘ０年の腐食量をＡ’、Ｘ０年を基準とした減衰パラメータをＢ’として、下記式（１８）を得ることができる。この式（１８）を用いることにより、Ｘ０年間を基準とした腐食量として、Ｘ＞Ｘ０の期間の腐食量を算出することができる。【０１２６】
【数１８】
2021100341.tif 000020
【０１２７】
金属材料の任意の初期期間の腐食量Ａ’と、減衰パラメータＢ’とをそれぞれ分けて予測し、上記式（１８）に示すように、初期期間以降の経過年数Ｘ’を減衰パラメータＢ’で累乗することにより、初期期間以降の長期腐食量を予測することができる。但し、上記式（１）の初期腐食量Ａは、１年間の腐食量を前提としている。そのため、上記式（１８）のＸ０の期間は、１年間から大きくずれた場合は想定しておらず、半年間から２年間程度であることが現実的な実用範囲と考えられる。
【符号の説明】
【０１２８】
１  腐食量マッピング装置
１０  入力部
２０  データベース
３０  演算部
３１  マップ作成部
３２  類似度算出部
３３  次元圧縮部
３４  腐食量予測部
４０  表示部
【図１】
2021100341.tif 000021
【図２】
2021100341.tif 000022
【図３】
2021100341.tif 000023
【図４】
2021100341.tif 000024
【図５】
2021100341.tif 000025
【図６】
2021100341.tif 000026
【国際調査報告】
2021100341.tif 000027
2021100341.tif 000028
2021100341.tif 000029
2021100341.tif 000030