ホーム > 特許ランキング > 日鉄ステンレス株式会社
知財求人 - 知財ポータルサイト「IP Force」
特許7607746鋼材の腐食量の予測方法、鋼材の腐食量予測システム、鋼材の腐食量予測プログラム及び鋼材の提案方法
(19)【発行国】日本国特許庁(JP)
(12)【公報種別】特許公報(B2)
(11)【特許番号】
(24)【登録日】2024-12-19
(45)【発行日】2024-12-27
(54)【発明の名称】鋼材の腐食量の予測方法、鋼材の腐食量予測システム、鋼材の腐食量予測プログラム及び鋼材の提案方法
(51)【国際特許分類】
G01N 17/00 20060101AFI20241220BHJP
【ＦＩ】
G01N17/00
【請求項の数】
19
(21)【出願番号】P 2023511051
(86)(22)【出願日】2022-03-23
(86)【国際出願番号】 JP2022013419
(87)【国際公開番号】W WO2022210151
(87)【国際公開日】2022-10-06
【審査請求日】2023-09-26
(31)【優先権主張番号】P 2021059457
(32)【優先日】2021-03-31
(33)【優先権主張国・地域又は機関】JP
(73)【特許権者】
【識別番号】503378420
【氏名又は名称】日鉄ステンレス株式会社
(74)【代理人】
【識別番号】100149548
【弁理士】
【氏名又は名称】松沼 泰史
(74)【代理人】
【識別番号】100140774
【弁理士】
【氏名又は名称】大浪 一徳
(74)【代理人】
【識別番号】100134359
【弁理士】
【氏名又は名称】勝俣 智夫
(74)【代理人】
【識別番号】100188592
【弁理士】
【氏名又は名称】山口 洋
(74)【代理人】
【識別番号】100217249
【弁理士】
【氏名又は名称】堀田 耕一郎
(74)【代理人】
【識別番号】100221279
【弁理士】
【氏名又は名称】山口 健吾
(74)【代理人】
【識別番号】100207686
【弁理士】
【氏名又は名称】飯田 恭宏
(74)【代理人】
【識別番号】100224812
【弁理士】
【氏名又は名称】井口 翔太
(72)【発明者】
【氏名】吉見 敏彦
(72)【発明者】
【氏名】平出 信彦
【審査官】外川 敬之
(56)【参考文献】
【文献】国際公開第０３／００６９５７（ＷＯ，Ａ１）
【文献】国際公開第２０２０／２３４９２１（ＷＯ，Ａ１）
【文献】特開２０１８－０２５４９７（ＪＰ，Ａ）
【文献】特開２０００－１３１３６３（ＪＰ，Ａ）
【文献】米国特許出願公開第２０１６／０２９０９７４（ＵＳ，Ａ１）
【文献】石川雄一，尾崎敏範，極地統計解析による局部腐食損傷の予測，日本金属学会会報，日本，日本金属学会，1982年03月20日，21(3)，133-140
(58)【調査した分野】(Int.Cl.，ＤＢ名)
Ｇ０１Ｎ １７／００
(57)【特許請求の範囲】
【請求項1】
鋼材を屋外に曝した場合の腐食量を評価する指標として、観測地点Ｘにおける観測時刻ｔＸｉ（ただしｉは１～ｎの整数）毎の風速、風向及び降雨量の気象観測値ＷＸから求まる鋼材の腐食性指標Ｑ（ｔＸｉ，ＷＸ）を用いる、鋼材の腐食量を予測する方法であって、前記観測地点Ｘが鋼材の利用予定地点Ｍである気象観測値ＷＭに基づいて、過去の観測時刻ｔＭｉ毎に鋼材の腐食性指標Ｑ（ｔＭｉ，ＷＭ）を求め、観測時刻ｔＭ１～ｔＭｎ間の腐食性指標Ｑ（ｔＭｉ，ＷＭ）の時系列データに対して極値統計解析を行うことにより、将来の評価期間ｔＦ1～ｔＦｎ内における腐食性指標の極大値Ｑｍａｘ（ｔＦｉ）を推測する極値推測ステップと、
試験地点Ｚにおける鋼材の腐食試験から得た実測腐食量と、前記観測地点Ｘが前記試験地点Ｚである気象観測値ＷＺに基づいて前記腐食試験の試験期間中の観測時刻ｔＺｉ毎に得た鋼材の腐食性指標Ｑ（ｔＺｉ，ＷＺ）の最大値との関係式を予め求めておき、この関係式に、前記極値推測ステップにおいて推測した腐食性指標の極大値Ｑｍａｘ（ｔＦｉ）を導入して、前記評価期間における鋼材の腐食量の推測値を得る腐食量推測ステップと、を備えることを特徴とする、鋼材の腐食量の予測方法。
【請求項2】
前記極値推測ステップ及び前記腐食量推測ステップにはそれぞれ、前記腐食性指標Ｑ（ｔＸｉ，ＷＸ）を求める計算ステップが備えられており、前記計算ステップは、
前記観測地点Ｘにおける前記観測時刻ｔＸｉ毎の風速ｕ（ｔＸｉ）、風向θｗ（ｔＸｉ）及び降雨量を少なくとも取得する第１ステップと、
前記観測時刻ｔＸｉ毎に、一次腐食性指標Ｑ１（ｔＸｉ，ＷＸ）を下記式（１）により計算する第２ステップと、
前記観測時刻ｔＸｉ毎に、二次腐食性指標Ｑ２（ｔＸｉ，ＷＸ）を下記式（２）により計算する第３ステップと、
前記観測時刻ｔＸｉ毎に、腐食性指標Ｑ（ｔＸｉ，ＷＸ）を下記式（３Ａ）により計算する第４ステップと、を備えている、請求項１に記載の、鋼材の腐食量の予測方法。
Ｑ１（ｔＸｉ，ＷＸ）＝（ｄ＋１）－０．６｛ｕ（ｔＸｉ）・ｃｏｓ｜θｓ－θｗ（ｔＸｉ）｜｝２ …（１）
Ｑ２（ｔＸｉ，ＷＸ）＝｛Ｑ２（ｔＸｉ－１，ＷＸ）＋Ｑ１（ｔＸｉ，ＷＸ）・ｓ（ｔＸｉ）・（ｔＸｉ－ｔＸｉ－１）｝・ｃ（ｔＸｉ） …（２）
Ｑ（ｔＸｉ，ＷＸ）＝Ｑ２（ｔＸｉ，ＷＸ）・ｐ（ｔＸｉ） …（３Ａ）
ただし、式（１）におけるｄは、前記観測地点Ｘに最も近い海岸と前記観測地点Ｘとの距離（ｍ）であり、ｕ（ｔＸｉ）は、観測時刻ｔＸｉの風速（ｍ／ｓ）であり、θｗ（ｔＸｉ）は北を０°とした場合の観測時刻ｔＸｉでの風向（°）であり、θｓは、前記観測地点Ｘに最も近い海岸からの前記観測地点Ｘの方位であって北を０°とした場合の方位（°）である。
式（２）におけるｃ（ｔＸｉ）は腐食性消失係数であって、観測時刻ｔＸｉ－１～ｔＸｉ間の降雨量が閾値超の場合にｃ（ｔＸｉ）＝０、観測時刻ｔＸｉ－１～ｔＸｉ間の降雨量が閾値以下の場合にｃ（ｔＸｉ）＝１とする。
式（２）におけるｓ（ｔＸｉ）は腐食性蓄積係数であって、観測時刻ｔＸｉ－１～ｔＸｉ間が昼間（日の出～日没までの期間）の場合にｓ（ｔＸｉ）＝０、観測時刻ｔＸｉ－１～ｔＸｉ間が夜間（日没～日の出までの期間）の場合にｓ（ｔＸｉ）＝１とする。
式（３Ａ）におけるｐ（ｔＸｉ）は腐食性発現係数であって、観測時刻ｔＸｉ－１～ｔＸｉ間が夜間（日没～日の出までの期間）の場合にｐ（ｔＸｉ）＝０、観測時刻ｔＸｉ－１～ｔＸｉ間が昼間（日の出～日没までの期間）の場合にｐ（ｔＸｉ）＝１、とする。
【請求項3】
前記極値推測ステップ及び前記腐食量推測ステップにはそれぞれ、前記腐食性指標Ｑ（ｔＸｉ，ＷＸ）を求める計算ステップが備えられており、前記計算ステップは、
前記観測地点Ｘにおける前記観測時刻ｔＸｉ毎の風速ｕ（ｔＸｉ）、風向θｗ（ｔＸｉ）及び降雨量を少なくとも取得する第１ステップと、
前記観測時刻ｔＸｉ毎に、一次腐食性指標Ｑ１（ｔＸｉ，ＷＸ）を下記式（１）により計算する第２ステップと、
前記観測時刻ｔＸｉ毎に、二次腐食性指標Ｑ２（ｔＸｉ，ＷＸ）を下記式（２）により計算する第３ステップと、
前記観測時刻ｔＸｉ毎に、腐食性指標Ｑ（ｔＸｉ，ＷＸ）を下記式（３Ｂ）により計算する第４ステップと、を備えている、請求項１に記載の、鋼材の腐食量の予測方法。
Ｑ１（ｔＸｉ，ＷＸ）＝（ｄ＋１）－０．６｛ｕ（ｔＸｉ）・ｃｏｓ｜θｓ－θｗ（ｔＸｉ）｜｝２ …（１）
Ｑ２（ｔＸｉ，ＷＸ）＝｛Ｑ２（ｔＸｉ－１，ＷＸ）＋Ｑ１（ｔＸｉ，ＷＸ）・ｓ（ｔＸｉ）・（ｔＸｉ－ｔＸｉ－１）｝・ｃ（ｔＸｉ） …（２）
Ｑ（ｔＸｉ，ＷＸ）＝１００．０４Ｔｂ・Ｑ２（ｔＸｉ，ＷＸ）・ｐ（ｔＸｉ） …（３Ｂ）
ただし、式（１）におけるｄは、前記観測地点Ｘに最も近い海岸と前記観測地点Ｘとの距離（ｍ）であり、ｕ（ｔＸｉ）は、観測時刻ｔＸｉの風速（ｍ／ｓ）であり、θｗ（ｔＸｉ）は北を０°とした場合の観測時刻ｔＸｉでの風向（°）であり、θｓは、前記観測地点Ｘに最も近い海岸からの前記観測地点Ｘの方位であって北を０°とした場合の方位（°）である。
式（２）におけるｃ（ｔＸｉ）は腐食性消失係数であって、観測時刻ｔＸｉ－１～ｔＸｉ間の降雨量が閾値超の場合にｃ（ｔＸｉ）＝０、観測時刻ｔＸｉ－１～ｔＸｉ間の降雨量が閾値以下の場合にｃ（ｔＸｉ）＝１とする。
式（２）におけるｓ（ｔＸｉ）は腐食性蓄積係数であって、観測時刻ｔＸｉ－１～ｔＸｉ間が昼間（日の出～日没までの期間）の場合にｓ（ｔＸｉ）＝０、観測時刻ｔＸｉ－１～ｔＸｉ間が夜間（日没～日の出までの期間）の場合にｓ（ｔＸｉ）＝１とする。
式（３Ｂ）におけるＴｂは、日の出時の気温Ｔｅｍｐ（℃）とした場合にＴｅｍｐ（℃）～（Ｔｅｍｐ＋２）℃の範囲の何れかの温度である。
式（３Ｂ）におけるｐ（ｔＸｉ）は腐食性発現係数であって、観測時刻ｔＸｉ－１～ｔＸｉ間が夜間（日没～日の出までの期間）の場合にｐ（ｔＸｉ）＝０、観測時刻ｔＸｉ－１～ｔＸｉ間が昼間（日の出～日没までの期間）の場合にｐ（ｔＸｉ）＝１、とする。
【請求項4】
前記極値推測ステップ及び前記腐食量推測ステップにはそれぞれ、前記腐食性指標Ｑ（ｔＸｉ，ＷＸ）を求める計算ステップが備えられており、前記計算ステップは、
前記観測地点Ｘにおける前記観測時刻ｔＸｉ毎の風速ｕ（ｔＸｉ）、風向θｗ（ｔＸｉ）及び降雨量を少なくとも取得する第１ステップと、
前記観測時刻ｔＸｉ毎に、一次腐食性指標Ｑ１（ｔＸｉ，ＷＸ）を下記式（１）により計算する第２ステップと、
前記観測時刻ｔＸｉ毎に、二次腐食性指標Ｑ２（ｔＸｉ，ＷＸ）を下記式（２）により計算する第３ステップと、
前記観測時刻ｔＸｉ毎に、腐食性指標Ｑ（ｔＸｉ，ＷＸ）を下記式（３Ｃ）により計算する第４ステップと、を備えている、請求項１に記載の、鋼材の腐食量の予測方法。
Ｑ１（ｔＸｉ，ＷＸ）＝（ｄ＋１）－０．６｛ｕ（ｔＸｉ）・ｃｏｓ｜θｓ－θｗ（ｔＸｉ）｜｝２ …（１）
Ｑ２（ｔＸｉ，ＷＸ）＝｛Ｑ２（ｔＸｉ－１，ＷＸ）＋Ｑ１（ｔＸｉ，ＷＸ）・ｓ（ｔＸｉ）・（ｔＸｉ－ｔＸｉ－１）｝・ｃ（ｔＸｉ） …（２）
Ｑ（ｔＸｉ，ＷＸ）＝１００．０４Ｔｃ・Ｑ２（ｔＸｉ，ＷＸ）・ｐ（ｔＸｉ） …（３Ｃ）
ただし、式（１）におけるｄは、前記観測地点Ｘに最も近い海岸と前記観測地点Ｘとの距離（ｍ）であり、ｕ（ｔＸｉ）は、観測時刻ｔＸｉの風速（ｍ／ｓ）であり、θｗ（ｔＸｉ）は北を０°とした場合の観測時刻ｔＸｉでの風向（°）であり、θｓは、前記観測地点Ｘに最も近い海岸からの前記観測地点Ｘの方位であって北を０°とした場合の方位（°）である。
式（２）におけるｃ（ｔＸｉ）は腐食性消失係数であって、観測時刻ｔＸｉ－１～ｔＸｉ間の降雨量が閾値超の場合にｃ（ｔＸｉ）＝０、観測時刻ｔＸｉ－１～ｔＸｉ間の降雨量が閾値以下の場合にｃ（ｔＸｉ）＝１とする。
式（２）におけるｓ（ｔＸｉ）は腐食性蓄積係数であって、観測時刻ｔＸｉにおける鋼材表面の水分付着量ｗｓ（ｔＸｉ）が０の場合にｓ（ｔＸｉ）＝０とし、前記ｗｓ（ｔＸｉ）が０超の場合にｓ（ｔＸｉ）＝１とする。前記ｗｓ（ｔＸｉ）は、前記観測時刻ｔＸｉにおける降雨および結露生成による水分付着量から水の蒸発量を差し引いた推測値であって、鋼材の表面温度、気温（℃）、気圧（ｈＰａ）、鋼材表面における風速（ｍ／ｓ）および相対湿度（％）から求まる値である。
式（３Ｃ）におけるＴｃは、前記観測時刻ｔＸｉより前の時点であって鋼材表面の水分付着量ｗｓ（ｔＸｉ）が０になる時点の気温（℃）である。
式（３Ｃ）におけるｐ（ｔＸｉ）は腐食性発現係数であって、前記ｗｓ（ｔＸｉ）が０を超える場合にｐ（ｔＸｉ）＝０とし、前記ｗｓ（ｔＸｉ）が０以下の場合にｐ（ｔＸｉ）＝１、とする。
【請求項5】
前記鋼材が、ステンレス鋼からなる鋼材であることを特徴とする、請求項１乃至請求項４の何れか一項に記載の、鋼材の腐食量の予測方法。
【請求項6】
鋼材を屋外に曝した場合の腐食量を評価する指標として、観測地点Ｘにおける観測時刻ｔＸｉ（ただしｉは１～ｎの整数）毎の風速、風向及び降雨量の気象観測値ＷＸから求まる鋼材の腐食性指標Ｑ（ｔＸｉ，ＷＸ）を用いる、電子計算機による鋼材の腐食量予測システムであって、前記観測地点Ｘが鋼材の利用予定地点Ｍである気象観測値ＷＭに基づいて、過去の観測時刻ｔＭｉ毎に鋼材の腐食性指標Ｑ（ｔＭｉ，ＷＭ）を求め、観測時刻ｔＭ１～ｔＭｎ間の腐食性指標Ｑ（ｔＭｉ，ＷＭ）の時系列データに対して極値統計解析を行うことにより、将来の評価期間ｔＦ1～ｔＦｎ内における腐食性指標の極大値Ｑｍａｘ（ｔＦｉ）を推測する極値推測部と、
試験地点Ｚにおける鋼材の腐食試験から得た実測腐食量と、前記観測地点Ｘが前記試験地点Ｚである気象観測値ＷＺに基づいて前記腐食試験の試験期間中の観測時刻ｔＺｉ毎に得た鋼材の腐食性指標Ｑ（ｔＺｉ，ＷＺ）の最大値との関係式を予め求めておき、この関係式に、前記極値推測部において推測した腐食性指標の極大値Ｑｍａｘ（ｔＦｉ）を導入して、前記評価期間における鋼材の腐食量の推測値を得る腐食量推測部と、を備えることを特徴とする、鋼材の腐食量予測システム。
【請求項7】
前記極値推測部及び前記腐食量推測部にはそれぞれ、前記腐食性指標Ｑ（ｔＸｉ，ＷＸ）を求める計算部が備えられており、前記計算部は、
前記観測地点Ｘにおける前記観測時刻ｔＸｉ毎の風速ｕ（ｔＸｉ）、風向θｗ（ｔＸｉ）及び降雨量を少なくとも取得する第１処理部と、
前記観測時刻ｔＸｉ毎に、一次腐食性指標Ｑ１（ｔＸｉ，ＷＸ）を下記式（１）により計算する第２処理部と、
前記観測時刻ｔＸｉ毎に、二次腐食性指標Ｑ２（ｔＸｉ，ＷＸ）を下記式（２）により計算する第３処理部と、
前記観測時刻ｔＸｉ毎に、腐食性指標Ｑ（ｔＸｉ，ＷＸ）を下記式（３Ａ）により計算する第４処理部と、を備えている、請求項６に記載の、鋼材の腐食量予測システム。
Ｑ１（ｔＸｉ，ＷＸ）＝（ｄ＋１）－０．６｛ｕ（ｔＸｉ）・ｃｏｓ｜θｓ－θｗ（ｔＸｉ）｜｝２ …（１）
Ｑ２（ｔＸｉ，ＷＸ）＝｛Ｑ２（ｔＸｉ－１，ＷＸ）＋Ｑ１（ｔＸｉ，ＷＸ）・ｓ（ｔＸｉ）・（ｔＸｉ－ｔＸｉ－１）｝・ｃ（ｔＸｉ） …（２）
Ｑ（ｔＸｉ，ＷＸ）＝Ｑ２（ｔＸｉ，ＷＸ）・ｐ（ｔＸｉ） …（３Ａ）
ただし、式（１）におけるｄは、前記観測地点Ｘに最も近い海岸と前記観測地点Ｘとの距離（ｍ）であり、ｕ（ｔＸｉ）は、観測時刻ｔＸｉの風速（ｍ／ｓ）であり、θｗ（ｔＸｉ）は北を０°とした場合の観測時刻ｔＸｉでの風向（°）であり、θｓは、前記観測地点Ｘに最も近い海岸からの前記観測地点Ｘの方位であって北を０°とした場合の方位（°）である。
式（２）におけるｃ（ｔＸｉ）は腐食性消失係数であって、観測時刻ｔＸｉ－１～ｔＸｉ間の降雨量が閾値超の場合にｃ（ｔＸｉ）＝０、観測時刻ｔＸｉ－１～ｔＸｉ間の降雨量が閾値以下の場合にｃ（ｔＸｉ）＝１とする。
式（２）におけるｓ（ｔＸｉ）は腐食性蓄積係数であって、観測時刻ｔＸｉ－１～ｔＸｉ間が昼間（日の出～日没までの期間）の場合にｓ（ｔＸｉ）＝０、観測時刻ｔＸｉ－１～ｔＸｉ間が夜間（日没～日の出までの期間）の場合にｓ（ｔＸｉ）＝１とする。
式（３Ａ）におけるｐ（ｔＸｉ）は腐食性発現係数であって、観測時刻ｔＸｉ－１～ｔＸｉ間が夜間（日没～日の出までの期間）の場合にｐ（ｔＸｉ）＝０、観測時刻ｔＸｉ－１～ｔＸｉ間が昼間（日の出～日没までの期間）の場合にｐ（ｔＸｉ）＝１、とする。
【請求項8】
前記極値推測部及び前記腐食量推測部にはそれぞれ、前記腐食性指標Ｑ（ｔＸｉ，ＷＸ）を求める計算部が備えられており、前記計算部は、
前記観測地点Ｘにおける前記観測時刻ｔＸｉ毎の風速ｕ（ｔＸｉ）、風向θｗ（ｔＸｉ）及び降雨量を少なくとも取得する第１処理部と、
前記観測時刻ｔＸｉ毎に、一次腐食性指標Ｑ１（ｔＸｉ，ＷＸ）を下記式（１）により計算する第２処理部と、
前記観測時刻ｔＸｉ毎に、二次腐食性指標Ｑ２（ｔＸｉ，ＷＸ）を下記式（２）により計算する第３処理部と、
前記観測時刻ｔＸｉ毎に、腐食性指標Ｑ（ｔＸｉ，ＷＸ）を下記式（３Ｂ）により計算する第４処理部と、を備えている、請求項６に記載の、鋼材の腐食量予測システム。
Ｑ１（ｔＸｉ，ＷＸ）＝（ｄ＋１）－０．６｛ｕ（ｔＸｉ）・ｃｏｓ｜θｓ－θｗ（ｔＸｉ）｜｝２ …（１）
Ｑ２（ｔＸｉ，ＷＸ）＝｛Ｑ２（ｔＸｉ－１，ＷＸ）＋Ｑ１（ｔＸｉ，ＷＸ）・ｓ（ｔＸｉ）・（ｔＸｉ－ｔＸｉ－１）｝・ｃ（ｔＸｉ） …（２）
Ｑ（ｔＸｉ，ＷＸ）＝１００．０４Ｔｂ・Ｑ２（ｔＸｉ，ＷＸ）・ｐ（ｔＸｉ） …（３Ｂ）
ただし、式（１）におけるｄは、前記観測地点Ｘに最も近い海岸と前記観測地点Ｘとの距離（ｍ）であり、ｕ（ｔＸｉ）は、観測時刻ｔＸｉの風速（ｍ／ｓ）であり、θｗ（ｔＸｉ）は北を０°とした場合の観測時刻ｔＸｉでの風向（°）であり、θｓは、前記観測地点Ｘに最も近い海岸からの前記観測地点Ｘの方位であって北を０°とした場合の方位（°）である。
式（２）におけるｃ（ｔＸｉ）は腐食性消失係数であって、観測時刻ｔＸｉ－１～ｔＸｉ間の降雨量が閾値超の場合にｃ（ｔＸｉ）＝０、観測時刻ｔＸｉ－１～ｔＸｉ間の降雨量が閾値以下の場合にｃ（ｔＸｉ）＝１とする。
式（２）におけるｓ（ｔＸｉ）は腐食性蓄積係数であって、観測時刻ｔＸｉ－１～ｔＸｉ間が昼間（日の出～日没までの期間）の場合にｓ（ｔＸｉ）＝０、観測時刻ｔＸｉ－１～ｔＸｉ間が夜間（日没～日の出までの期間）の場合にｓ（ｔＸｉ）＝１とする。
式（３Ｂ）におけるＴｂは、日の出時の気温Ｔｅｍｐ（℃）とした場合にＴｅｍｐ（℃）～（Ｔｅｍｐ＋２）℃の範囲の何れかの温度である。
式（３Ｂ）におけるｐ（ｔＸｉ）は腐食性発現係数であって、観測時刻ｔＸｉ－１～ｔＸｉ間が夜間（日没～日の出までの期間）の場合にｐ（ｔＸｉ）＝０、観測時刻ｔＸｉ－１～ｔＸｉ間が昼間（日の出～日没までの期間）の場合にｐ（ｔＸｉ）＝１、とする。
【請求項9】
前記極値推測部及び前記腐食量推測部にはそれぞれ、前記腐食性指標Ｑ（ｔＸｉ，ＷＸ）を求める計算部が備えられており、前記計算部は、
前記観測地点Ｘにおける前記観測時刻ｔＸｉ毎の風速ｕ（ｔＸｉ）、風向θｗ（ｔＸｉ）及び降雨量を少なくとも取得する第１処理部と、
前記観測時刻ｔＸｉ毎に、一次腐食性指標Ｑ１（ｔＸｉ，ＷＸ）を下記式（１）により計算する第２処理部と、
前記観測時刻ｔＸｉ毎に、二次腐食性指標Ｑ２（ｔＸｉ，ＷＸ）を下記式（２）により計算する第３処理部と、
前記観測時刻ｔＸｉ毎に、腐食性指標Ｑ（ｔＸｉ，ＷＸ）を下記式（３Ｃ）により計算する第４処理部と、を備えている、請求項６に記載の、鋼材の腐食量予測システム。
Ｑ１（ｔＸｉ，ＷＸ）＝（ｄ＋１）－０．６｛ｕ（ｔＸｉ）・ｃｏｓ｜θｓ－θｗ（ｔＸｉ）｜｝２ …（１）
Ｑ２（ｔＸｉ，ＷＸ）＝｛Ｑ２（ｔＸｉ－１，ＷＸ）＋Ｑ１（ｔＸｉ，ＷＸ）・ｓ（ｔＸｉ）・（ｔＸｉ－ｔＸｉ－１）｝・ｃ（ｔＸｉ） …（２）
Ｑ（ｔＸｉ，ＷＸ）＝１００．０４Ｔｃ・Ｑ２（ｔＸｉ，ＷＸ）・ｐ（ｔＸｉ） …（３Ｃ）
ただし、式（１）におけるｄは、前記観測地点Ｘに最も近い海岸と前記観測地点Ｘとの距離（ｍ）であり、ｕ（ｔＸｉ）は、観測時刻ｔＸｉの風速（ｍ／ｓ）であり、θｗ（ｔＸｉ）は北を０°とした場合の観測時刻ｔＸｉでの風向（°）であり、θｓは、前記観測地点Ｘに最も近い海岸からの前記観測地点Ｘの方位であって北を０°とした場合の方位（°）である。
式（２）におけるｃ（ｔＸｉ）は腐食性消失係数であって、観測時刻ｔＸｉ－１～ｔＸｉ間の降雨量が閾値超の場合にｃ（ｔＸｉ）＝０、観測時刻ｔＸｉ－１～ｔＸｉ間の降雨量が閾値以下の場合にｃ（ｔＸｉ）＝１とする。
式（２）におけるｓ（ｔＸｉ）は腐食性蓄積係数であって、観測時刻ｔＸｉにおける鋼材表面の水分付着量ｗｓ（ｔＸｉ）が０の場合にｓ（ｔＸｉ）＝０とし、前記ｗｓ（ｔＸｉ）が０超の場合にｓ（ｔＸｉ）＝１とする。前記ｗｓ（ｔＸｉ）は、前記観測時刻ｔＸｉにおける降雨および結露生成による水分付着量から水の蒸発量を差し引いた推測値であって、鋼材の表面温度、気温（℃）、気圧（ｈＰａ）、鋼材表面における風速（ｍ／ｓ）および相対湿度（％）から求まる値である。
式（３Ｃ）におけるＴｃは、前記観測時刻ｔＸｉより前の時点であって鋼材表面の水分付着量ｗｓ（ｔＸｉ）が０になる時点の気温（℃）である。
式（３Ｃ）におけるｐ（ｔＸｉ）は腐食性発現係数であって、前記ｗｓ（ｔＸｉ）が０を超える場合にｐ（ｔＸｉ）＝０とし、前記ｗｓ（ｔＸｉ）が０以下の場合にｐ（ｔＸｉ）＝１、とする。
【請求項10】
前記鋼材が、ステンレス鋼からなる鋼材であることを特徴とする、請求項６乃至請求項９の何れか一項に記載の、鋼材の腐食量予測システム。
【請求項11】
鋼材を屋外に曝した場合の腐食量を評価する指標として、観測地点Ｘにおける観測時刻ｔＸｉ（ただしｉは１～ｎの整数）毎の風速、風向及び降雨量の気象観測値ＷＸから求まる鋼材の腐食性指標Ｑ（ｔＸｉ，ＷＸ）を用いる、電子計算機に利用される鋼材の腐食量予測プログラムであって、前記観測地点Ｘが鋼材の利用予定地点Ｍである気象観測値ＷＭに基づいて、過去の観測時刻ｔＭｉ毎に鋼材の腐食性指標Ｑ（ｔＭｉ，ＷＭ）を求め、観測時刻ｔＭ１～ｔＭｎ間の腐食性指標Ｑ（ｔＭｉ，ＷＭ）の時系列データに対して極値統計解析を行うことにより、将来の評価期間ｔＦ1～ｔＦｎ内における腐食性指標の極大値Ｑｍａｘ（ｔＦｉ）を推測する極値推測ステップと、
試験地点Ｚにおける鋼材の腐食試験から得た実測腐食量と、前記観測地点Ｘが前記試験地点Ｚである気象観測値ＷＺに基づいて前記腐食試験の試験期間中の観測時刻ｔＺｉ毎に得た鋼材の腐食性指標Ｑ（ｔＺｉ，ＷＺ）の最大値との関係式を予め求めておき、この関係式に、前記極値推測ステップにおいて推測した腐食性指標の極大値Ｑｍａｘ（ｔＦｉ）を導入して、前記評価期間における鋼材の腐食量の推測値を得る腐食量推測ステップと、を備えることを特徴とする、鋼材の腐食量予測プログラム。
【請求項12】
前記極値推測ステップ及び前記腐食量推測ステップにはそれぞれ、前記腐食性指標Ｑ（ｔＸｉ，ＷＸ）を求める計算ステップが備えられており、前記計算ステップは、
前記観測地点Ｘにおける前記観測時刻ｔＸｉ毎の風速ｕ（ｔＸｉ）、風向θｗ（ｔＸｉ）及び降雨量を少なくとも取得する第１ステップと、
前記観測時刻ｔＸｉ毎に、一次腐食性指標Ｑ１（ｔＸｉ，ＷＸ）を下記式（１）により計算する第２ステップと、
前記観測時刻ｔＸｉ毎に、二次腐食性指標Ｑ２（ｔＸｉ，ＷＸ）を下記式（２）により計算する第３ステップと、
前記観測時刻ｔＸｉ毎に、腐食性指標Ｑ（ｔＸｉ，ＷＸ）を下記式（３Ａ）により計算する第４ステップと、を備えている請求項１１に記載の、鋼材の腐食量予測プログラム。
Ｑ１（ｔＸｉ，ＷＸ）＝（ｄ＋１）－０．６｛ｕ（ｔＸｉ）・ｃｏｓ｜θｓ－θｗ（ｔＸｉ）｜｝２ …（１）
Ｑ２（ｔＸｉ，ＷＸ）＝｛Ｑ２（ｔＸｉ－１，ＷＸ）＋Ｑ１（ｔＸｉ，ＷＸ）・ｓ（ｔＸｉ）・（ｔＸｉ－ｔＸｉＳ－１）｝・ｃ（ｔＸｉ） …（２）
Ｑ（ｔＸｉ，ＷＸ）＝Ｑ２（ｔＸｉ，ＷＸ）・ｐ（ｔＸｉ） …（３Ａ）
ただし、式（１）におけるｄは、前記観測地点Ｘに最も近い海岸と前記観測地点Ｘとの距離（ｍ）であり、ｕ（ｔＸｉ）は、観測時刻ｔＸｉの風速（ｍ／ｓ）であり、θｗ（ｔＸｉ）は北を０°とした場合の観測時刻ｔＸｉでの風向（°）であり、θｓは、前記観測地点Ｘに最も近い海岸からの前記観測地点Ｘの方位であって北を０°とした場合の方位（°）である。
式（２）におけるｃ（ｔＸｉ）は腐食性消失係数であって、観測時刻ｔＸｉ－１～ｔＸｉ間の降雨量が閾値超の場合にｃ（ｔＸｉ）＝０、観測時刻ｔＸｉ－１～ｔＸｉ間の降雨量が閾値以下の場合にｃ（ｔＸｉ）＝１とする。
式（２）におけるｓ（ｔＸｉ）は腐食性蓄積係数であって、観測時刻ｔＸｉ－１～ｔＸｉ間が昼間（日の出～日没までの期間）の場合にｓ（ｔＸｉ）＝０、観測時刻ｔＸｉ－１～ｔＸｉ間が夜間（日没～日の出までの期間）の場合にｓ（ｔＸｉ）＝１とする。
式（３Ａ）におけるｐ（ｔＸｉ）は腐食性発現係数であって、観測時刻ｔＸｉ－１～ｔＸｉ間が夜間（日没～日の出までの期間）の場合にｐ（ｔＸｉ）＝０、観測時刻ｔＸｉ－１～ｔＸｉ間が昼間（日の出～日没までの期間）の場合にｐ（ｔＸｉ）＝１、とする。
【請求項13】
前記極値推測ステップ及び前記腐食量推測ステップにはそれぞれ、前記腐食性指標Ｑ（ｔＸｉ，ＷＸ）を求める計算ステップが備えられており、前記計算ステップは、
前記観測地点Ｘにおける前記観測時刻ｔＸｉ毎の風速ｕ（ｔＸｉ）、風向θｗ（ｔＸｉ）及び降雨量を少なくとも取得する第１ステップと、
前記観測時刻ｔＸｉ毎に、一次腐食性指標Ｑ１（ｔＸｉ，ＷＸ）を下記式（１）により計算する第２ステップと、
前記観測時刻ｔＸｉ毎に、二次腐食性指標Ｑ２（ｔＸｉ，ＷＸ）を下記式（２）により計算する第３ステップと、
前記観測時刻ｔＸｉ毎に、腐食性指標Ｑ（ｔＸｉ，ＷＸ）を下記式（３Ｂ）により計算する第４ステップと、を備えている請求項１１に記載の、鋼材の腐食量予測プログラム。
Ｑ１（ｔＸｉ，ＷＸ）＝（ｄ＋１）－０．６｛ｕ（ｔＸｉ）・ｃｏｓ｜θｓ－θｗ（ｔＸｉ）｜｝２ …（１）
Ｑ２（ｔＸｉ，ＷＸ）＝｛Ｑ２（ｔＸｉ－１，ＷＸ）＋Ｑ１（ｔＸｉ，ＷＸ）・ｓ（ｔＸｉ）・（ｔＸｉ－ｔＸｉＳ－１）｝・ｃ（ｔＸｉ） …（２）
Ｑ（ｔＸｉ，ＷＸ）＝１００．０４Ｔｂ・Ｑ２（ｔＸｉ，ＷＸ）・ｐ（ｔＸｉ） …（３Ｂ）
ただし、式（１）におけるｄは、前記観測地点Ｘに最も近い海岸と前記観測地点Ｘとの距離（ｍ）であり、ｕ（ｔＸｉ）は、観測時刻ｔＸｉの風速（ｍ／ｓ）であり、θｗ（ｔＸｉ）は北を０°とした場合の観測時刻ｔＸｉでの風向（°）であり、θｓは、前記観測地点Ｘに最も近い海岸からの前記観測地点Ｘの方位であって北を０°とした場合の方位（°）である。
式（２）におけるｃ（ｔＸｉ）は腐食性消失係数であって、観測時刻ｔＸｉ－１～ｔＸｉ間の降雨量が閾値超の場合にｃ（ｔＸｉ）＝０、観測時刻ｔＸｉ－１～ｔＸｉ間の降雨量が閾値以下の場合にｃ（ｔＸｉ）＝１とする。
式（２）におけるｓ（ｔＸｉ）は腐食性蓄積係数であって、観測時刻ｔＸｉ－１～ｔＸｉ間が昼間（日の出～日没までの期間）の場合にｓ（ｔＸｉ）＝０、観測時刻ｔＸｉ－１～ｔＸｉ間が夜間（日没～日の出までの期間）の場合にｓ（ｔＸｉ）＝１とする。
式（３Ｂ）におけるＴｂは、日の出時の気温Ｔｅｍｐ（℃）とした場合にＴｅｍｐ（℃）～（Ｔｅｍｐ＋２）℃の範囲の何れかの温度である。
式（３Ｂ）におけるｐ（ｔＸｉ）は腐食性発現係数であって、観測時刻ｔＸｉ－１～ｔＸｉ間が夜間（日没～日の出までの期間）の場合にｐ（ｔＸｉ）＝０、観測時刻ｔＸｉ－１～ｔＸｉ間が昼間（日の出～日没までの期間）の場合にｐ（ｔＸｉ）＝１、とする。
【請求項14】
前記極値推測ステップ及び前記腐食量推測ステップにはそれぞれ、前記腐食性指標Ｑ（ｔＸｉ，ＷＸ）を求める計算ステップが備えられており、前記計算ステップは、
前記観測地点Ｘにおける前記観測時刻ｔＸｉ毎の風速ｕ（ｔＸｉ）、風向θｗ（ｔＸｉ）及び降雨量を少なくとも取得する第１ステップと、
前記観測時刻ｔＸｉ毎に、一次腐食性指標Ｑ１（ｔＸｉ，ＷＸ）を下記式（１）により計算する第２ステップと、
前記観測時刻ｔＸｉ毎に、二次腐食性指標Ｑ２（ｔＸｉ，ＷＸ）を下記式（２）により計算する第３ステップと、
前記観測時刻ｔＸｉ毎に、腐食性指標Ｑ（ｔＸｉ，ＷＸ）を下記式（３Ｃ）により計算する第４ステップと、を備えている請求項１１に記載の、鋼材の腐食量予測プログラム。
Ｑ１（ｔＸｉ，ＷＸ）＝（ｄ＋１）－０．６｛ｕ（ｔＸｉ）・ｃｏｓ｜θｓ－θｗ（ｔＸｉ）｜｝２ …（１）
Ｑ２（ｔＸｉ，ＷＸ）＝｛Ｑ２（ｔＸｉ－１，ＷＸ）＋Ｑ１（ｔＸｉ，ＷＸ）・ｓ（ｔＸｉ）・（ｔＸｉ－ｔＸｉＳ－１）｝・ｃ（ｔＸｉ） …（２）
Ｑ（ｔＸｉ，ＷＸ）＝１００．０４Ｔｃ・Ｑ２（ｔＸｉ，ＷＸ）・ｐ（ｔＸｉ） …（３Ｃ）
ただし、式（１）におけるｄは、前記観測地点Ｘに最も近い海岸と前記観測地点Ｘとの距離（ｍ）であり、ｕ（ｔＸｉ）は、観測時刻ｔＸｉの風速（ｍ／ｓ）であり、θｗ（ｔＸｉ）は北を０°とした場合の観測時刻ｔＸｉでの風向（°）であり、θｓは、前記観測地点Ｘに最も近い海岸からの前記観測地点Ｘの方位であって北を０°とした場合の方位（°）である。
式（２）におけるｃ（ｔＸｉ）は腐食性消失係数であって、観測時刻ｔＸｉ－１～ｔＸｉ間の降雨量が閾値超の場合にｃ（ｔＸｉ）＝０、観測時刻ｔＸｉ－１～ｔＸｉ間の降雨量が閾値以下の場合にｃ（ｔＸｉ）＝１とする。
式（２）におけるｓ（ｔＸｉ）は腐食性蓄積係数であって、観測時刻ｔＸｉにおける鋼材表面の水分付着量ｗｓ（ｔＸｉ）が０の場合にｓ（ｔＸｉ）＝０とし、前記ｗｓ（ｔＸｉ）が０超の場合にｓ（ｔＸｉ）＝１とする。前記ｗｓ（ｔＸｉ）は、前記観測時刻ｔＸｉにおける降雨および結露生成による水分付着量から水の蒸発量を差し引いた推測値であって、鋼材の表面温度、気温（℃）、気圧（ｈＰａ）、鋼材表面における風速（ｍ／ｓ）および相対湿度（％）から求まる値である。
式（３Ｃ）におけるＴｃは、前記観測時刻ｔＸｉより前の時点であって鋼材表面の水分付着量ｗｓ（ｔＸｉ）が０になる時点の気温（℃）である。
式（３Ｃ）におけるｐ（ｔＸｉ）は腐食性発現係数であって、前記ｗｓ（ｔＸｉ）が０を超える場合にｐ（ｔＸｉ）＝０とし、前記ｗｓ（ｔＸｉ）が０以下の場合にｐ（ｔＸｉ）＝１、とする。
【請求項15】
前記鋼材が、ステンレス鋼からなる鋼材であることを特徴とする、請求項１１乃至請求項１４の何れか一項に記載の、鋼材の腐食量予測プログラム。
【請求項16】
請求項１乃至請求項５の何れか一項に記載の鋼材の腐食量の予測方法により、利用予定地点Ｍにおける将来にわたる鋼材の腐食量を予測する予測工程と、前記予測工程によって得られた腐食量の予測値または腐食量の予測値に対応する鋼材表面の画像を、営業担当者が顧客に提示する工程と、を備えた、鋼材の提案方法。
【請求項17】
請求項１乃至請求項５の何れか一項に記載の鋼材の腐食量の予測方法により、利用予定地点Ｍにおける将来にわたる鋼材の腐食量を、複数種類の鋼材についてそれぞれ予測する予測工程と、前記予測工程によって得られた鋼材毎の腐食量の予測値または腐食量の予測値に対応する鋼材表面の画像を、営業担当者が顧客に提示する工程と、を備えた、鋼材の提案方法。
【請求項18】
前記腐食量の予測値が、少なくとも小数点第一位までのレイティングナンバーである、請求項１６または請求項１７に記載の鋼材の提案方法。
【請求項19】
前記腐食量の予測値に対応する鋼材表面の画像が、少なくとも小数点第一位までのレイティングナンバーに対応する画像である、請求項１６または請求項１７に記載の鋼材の提案方法。
【発明の詳細な説明】
【技術分野】
【0001】
本発明は、鋼材の腐食量の予測方法、鋼材の腐食量予測システム、鋼材の腐食量予測プログラム及び鋼材の提案方法に関する。
本願は、２０２１年３月３１日に、日本に出願された特願２０２１－５９４５７号に基づき優先権を主張し、その内容をここに援用する。
本願は、２０２１年３月３１日に、日本に出願された特願２０２１－５９４５７号に基づき優先権を主張し、その内容をここに援用する。
【背景技術】
【0002】
建築用鋼材は、建物外壁、柱、梁などの素材として用いられる。建築用鋼材の購入を検討若しくは購入を決めた需要者は、建築用鋼材を長期間に渡り屋外環境下で使用した場合の将来の腐食量を予測したいという要望を持っている。
【0003】
特許文献１には、電子計算機によって、耐候性鋼を使用する予定の使用予定位置における年間ぬれ時間、年平均風速、年平均気温、飛来塩分量及び硫黄酸化物量を含む外因性の腐食情報並びに耐候性鋼の成分に関する内因性の腐食情報を用いて、耐候性鋼の予測腐食量を計算する工程を有する耐候性鋼の腐食予測方法が記載されている。
【0004】
ところで、鋼材の腐食は、突発的な気象変動に大きく影響される場合がある。例えば、海岸に近い地点において、台風接近に伴い風速が通常の気象条件に比べて著しく増大するような事象が発生したり、気圧配置の影響により短時間の間に風速が著しく増大する事象が発生すると、比較的短期間のうちに建物用鋼材に多量の塩分が付着する。これにより、鋼材の腐食量が著しく増大する場合がある。
【0005】
特許文献１では、将来の予測腐食量を求める際に、年平均風速や年平均気温などの年間平均値に基づいた予測を行っている。しかし、年間平均値に基づいた予測では、データの構造上、上述のような鋼材の腐食に大きな影響を与え得る短期間の突発的な気象変動が考慮されない。よって、従来の予測方法よりも腐食量の予測精度を向上させる余地がある。
【先行技術文献】
【特許文献】
【0006】
【文献】国際公開第２００３／００６９５７号
【発明の概要】
【発明が解決しようとする課題】
【0007】
本発明は、鋼材の腐食量の予測精度を向上させることが可能な、鋼材の腐食量の予測方法、鋼材の腐食量予測システム及び鋼材の腐食量予測プログラムを提供することを課題とする。また、本発明は、鋼材の腐食量の予測方法によって予測された鋼材の腐食量を利用した、鋼材の提案方法を提供することを課題とする。
【課題を解決するための手段】
【0008】
上記課題を解決するため、本発明は以下の構成を採用する。
［１］ 鋼材を屋外に曝した場合の腐食量を評価する指標として、観測地点Ｘにおける観測時刻ｔＸｉ（ただしｉは１～ｎの整数）毎の風速、風向及び降雨量の気象観測値ＷＸから求まる鋼材の腐食性指標Ｑ（ｔＸｉ，ＷＸ）を用いる、鋼材の腐食量を予測する方法であって、
前記観測地点Ｘが鋼材の利用予定地点Ｍである気象観測値ＷＭに基づいて、過去の観測時刻ｔＭｉ毎に鋼材の腐食性指標Ｑ（ｔＭｉ，ＷＭ）を求め、観測時刻ｔＭ１～ｔＭｎ間の腐食性指標Ｑ（ｔＭｉ，ＷＭ）の時系列データに対して極値統計解析を行うことにより、将来の評価期間ｔＦ1～ｔＦｎ内における腐食性指標の極大値Ｑｍａｘ（ｔＦｉ）を推測する極値推測ステップと、
試験地点Ｚにおける鋼材の腐食試験から得た実測腐食量と、前記観測地点Ｘが前記試験地点Ｚである気象観測値ＷＺに基づいて前記腐食試験の試験期間中の観測時刻ｔＺｉ毎に得た鋼材の腐食性指標Ｑ（ｔＺｉ，ＷＺ）の最大値との関係式を予め求めておき、この関係式に、前記極値推測ステップにおいて推測した腐食性指標の極大値Ｑｍａｘ（ｔＦｉ）を導入して、前記評価期間における鋼材の腐食量の推測値を得る腐食量推測ステップと、を備えることを特徴とする、鋼材の腐食量の予測方法。
なお、上記関係式は、鋼材の腐食試験から得た実測腐食量と、鋼材の腐食性指標Ｑ（ｔＺｉ，ＷＺ）の最大値の対数との関係式としてもよい。
［２］ 前記極値推測ステップ及び前記腐食量推測ステップにはそれぞれ、前記腐食性指標Ｑ（ｔＸｉ，ＷＸ）を求める計算ステップが備えられており、
前記計算ステップは、
前記観測地点Ｘにおける前記観測時刻ｔＸｉ毎の風速ｕ（ｔＸｉ）、風向θｗ（ｔＸｉ）及び降雨量を少なくとも取得する第１ステップと、
前記観測時刻ｔＸｉ毎に、一次腐食性指標Ｑ１（ｔＸｉ，ＷＸ）を下記式（１）により計算する第２ステップと、
前記観測時刻ｔＸｉ毎に、二次腐食性指標Ｑ２（ｔＸｉ，ＷＸ）を下記式（２）により計算する第３ステップと、
前記観測時刻ｔＸｉ毎に、腐食性指標Ｑ（ｔＸｉ，ＷＸ）を下記式（３Ａ）により計算する第４ステップと、を備えている、［１］に記載の、鋼材の腐食量の予測方法。
Ｑ１（ｔＸｉ，ＷＸ）＝（ｄ＋１）－０．６｛ｕ（ｔＸｉ）・ｃｏｓ｜θｓ－θｗ（ｔＸｉ）｜｝２ …（１）
Ｑ２（ｔＸｉ，ＷＸ）＝｛Ｑ２（ｔＸｉ－１，ＷＸ）＋Ｑ１（ｔＸｉ，ＷＸ）・ｓ（ｔＸｉ）・（ｔＸｉ－ｔＸｉ－１）｝・ｃ（ｔＸｉ） …（２）
Ｑ（ｔＸｉ，ＷＸ）＝Ｑ２（ｔＸｉ，ＷＸ）・ｐ（ｔＸｉ） …（３Ａ）
ただし、式（１）におけるｄは、前記観測地点Ｘに最も近い海岸と前記観測地点Ｘとの距離（ｍ）であり、ｕ（ｔＸｉ）は、観測時刻ｔＸｉの風速（ｍ／ｓ）であり、θｗ（ｔＸｉ）は北を０°とした場合の観測時刻ｔＸｉでの風向（°）であり、θｓは、前記観測地点Ｘに最も近い海岸からの前記観測地点Ｘの方位であって北を０°とした場合の方位（°）である。
式（２）におけるｃ（ｔＸｉ）は腐食性消失係数であって、観測時刻ｔＸｉ－１～ｔＸｉ間の降雨量が閾値超の場合にｃ（ｔＸｉ）＝０、観測時刻ｔＸｉ－１～ｔＸｉ間の降雨量が閾値以下の場合にｃ（ｔＸｉ）＝１とする。
式（２）におけるｓ（ｔＸｉ）は腐食性蓄積係数であって、観測時刻ｔＸｉ－１～ｔＸｉ間が昼間（日の出～日没までの期間）の場合にｓ（ｔＸｉ）＝０、観測時刻ｔＸｉ－１～ｔＸｉ間が夜間（日没～日の出までの期間）の場合にｓ（ｔＸｉ）＝１とする。
式（３Ａ）におけるｐ（ｔＸｉ）は腐食性発現係数であって、観測時刻ｔＸｉ－１～ｔＸｉ間が夜間（日没～日の出までの期間）の場合にｐ（ｔＸｉ）＝０、観測時刻ｔＸｉ－１～ｔＸｉ間が昼間（日の出～日没までの期間）の場合にｐ（ｔＸｉ）＝１、とする。
［３］ 前記極値推測ステップ及び前記腐食量推測ステップにはそれぞれ、前記腐食性指標Ｑ（ｔＸｉ，ＷＸ）を求める計算ステップが備えられており、
前記計算ステップは、
前記観測地点Ｘにおける前記観測時刻ｔＸｉ毎の風速ｕ（ｔＸｉ）、風向θｗ（ｔＸｉ）及び降雨量を少なくとも取得する第１ステップと、
前記観測時刻ｔＸｉ毎に、一次腐食性指標Ｑ１（ｔＸｉ，ＷＸ）を下記式（１）により計算する第２ステップと、
前記観測時刻ｔＸｉ毎に、二次腐食性指標Ｑ２（ｔＸｉ，ＷＸ）を下記式（２）により計算する第３ステップと、
前記観測時刻ｔＸｉ毎に、腐食性指標Ｑ（ｔＸｉ，ＷＸ）を下記式（３Ｂ）により計算する第４ステップと、を備えている、［１］に記載の、鋼材の腐食量の予測方法。
Ｑ１（ｔＸｉ，ＷＸ）＝（ｄ＋１）－０．６｛ｕ（ｔＸｉ）・ｃｏｓ｜θｓ－θｗ（ｔＸｉ）｜｝２ …（１）
Ｑ２（ｔＸｉ，ＷＸ）＝｛Ｑ２（ｔＸｉ－１，ＷＸ）＋Ｑ１（ｔＸｉ，ＷＸ）・ｓ（ｔＸｉ）・（ｔＸｉ－ｔＸｉ－１）｝・ｃ（ｔＸｉ） …（２）
Ｑ（ｔＸｉ，ＷＸ）＝１００．０４Ｔｂ・Ｑ２（ｔＸｉ，ＷＸ）・ｐ（ｔＸｉ） …（３Ｂ）
ただし、式（１）におけるｄは、前記観測地点Ｘに最も近い海岸と前記観測地点Ｘとの距離（ｍ）であり、ｕ（ｔＸｉ）は、観測時刻ｔＸｉの風速（ｍ／ｓ）であり、θｗ（ｔＸｉ）は北を０°とした場合の観測時刻ｔＸｉでの風向（°）であり、θｓは、前記観測地点Ｘに最も近い海岸からの前記観測地点Ｘの方位であって北を０°とした場合の方位（°）である。
式（２）におけるｃ（ｔＸｉ）は腐食性消失係数であって、観測時刻ｔＸｉ－１～ｔＸｉ間の降雨量が閾値超の場合にｃ（ｔＸｉ）＝０、観測時刻ｔＸｉ－１～ｔＸｉ間の降雨量が閾値以下の場合にｃ（ｔＸｉ）＝１とする。
式（２）におけるｓ（ｔＸｉ）は腐食性蓄積係数であって、観測時刻ｔＸｉ－１～ｔＸｉ間が昼間（日の出～日没までの期間）の場合にｓ（ｔＸｉ）＝０、観測時刻ｔＸｉ－１～ｔＸｉ間が夜間（日没～日の出までの期間）の場合にｓ（ｔＸｉ）＝１とする。
式（３Ｂ）におけるＴｂは、日の出時の気温Ｔｅｍｐ（℃）とした場合にＴｅｍｐ（℃）～（Ｔｅｍｐ＋２）℃の範囲の何れかの温度である。
式（３Ｂ）におけるｐ（ｔＸｉ）は腐食性発現係数であって、観測時刻ｔＸｉ－１～ｔＸｉ間が夜間（日没～日の出までの期間）の場合にｐ（ｔＸｉ）＝０、観測時刻ｔＸｉ－１～ｔＸｉ間が昼間（日の出～日没までの期間）の場合にｐ（ｔＸｉ）＝１、とする。
［４］ 前記極値推測ステップ及び前記腐食量推測ステップにはそれぞれ、前記腐食性指標Ｑ（ｔＸｉ，ＷＸ）を求める計算ステップが備えられており、
前記計算ステップは、
前記観測地点Ｘにおける前記観測時刻ｔＸｉ毎の風速ｕ（ｔＸｉ）、風向θｗ（ｔＸｉ）及び降雨量を少なくとも取得する第１ステップと、
前記観測時刻ｔＸｉ毎に、一次腐食性指標Ｑ１（ｔＸｉ，ＷＸ）を下記式（１）により計算する第２ステップと、
前記観測時刻ｔＸｉ毎に、二次腐食性指標Ｑ２（ｔＸｉ，ＷＸ）を下記式（２）により計算する第３ステップと、
前記観測時刻ｔＸｉ毎に、腐食性指標Ｑ（ｔＸｉ，ＷＸ）を下記式（３Ｃ）により計算する第４ステップと、を備えている、［１］に記載の、鋼材の腐食量の予測方法。
Ｑ１（ｔＸｉ，ＷＸ）＝（ｄ＋１）－０．６｛ｕ（ｔＸｉ）・ｃｏｓ｜θｓ－θｗ（ｔＸｉ）｜｝２ …（１）
Ｑ２（ｔＸｉ，ＷＸ）＝｛Ｑ２（ｔＸｉ－１，ＷＸ）＋Ｑ１（ｔＸｉ，ＷＸ）・ｓ（ｔＸｉ）・（ｔＸｉ－ｔＸｉ－１）｝・ｃ（ｔＸｉ） …（２）
Ｑ（ｔＸｉ，ＷＸ）＝１００．０４Ｔｃ・Ｑ２（ｔＸｉ，ＷＸ）・ｐ（ｔＸｉ） …（３Ｃ）
ただし、式（１）におけるｄは、前記観測地点Ｘに最も近い海岸と前記観測地点Ｘとの距離（ｍ）であり、ｕ（ｔＸｉ）は、観測時刻ｔＸｉの風速（ｍ／ｓ）であり、θｗ（ｔＸｉ）は北を０°とした場合の観測時刻ｔＸｉでの風向（°）であり、θｓは、前記観測地点Ｘに最も近い海岸からの前記観測地点Ｘの方位であって北を０°とした場合の方位（°）である。
式（２）におけるｃ（ｔＸｉ）は腐食性消失係数であって、観測時刻ｔＸｉ－１～ｔＸｉ間の降雨量が閾値超の場合にｃ（ｔＸｉ）＝０、観測時刻ｔＸｉ－１～ｔＸｉ間の降雨量が閾値以下の場合にｃ（ｔＸｉ）＝１とする。
式（２）におけるｓ（ｔＸｉ）は腐食性蓄積係数であって、観測時刻ｔＸｉにおける鋼材表面の水分付着量ｗｓ（ｔＸｉ）が０の場合にｓ（ｔＸｉ）＝０とし、前記ｗｓ（ｔＸｉ）が０超の場合にｓ（ｔＸｉ）＝１とする。前記ｗｓ（ｔＸｉ）は、前記観測時刻ｔＸｉにおける降雨および結露生成による水分付着量から水の蒸発量を差し引いた推測値であって、鋼材の表面温度、気温（℃）、気圧（ｈＰａ）、鋼材表面における風速（ｍ／ｓ）および相対湿度（％）から求まる値である。
式（３Ｃ）におけるＴｃは、前記観測時刻ｔＸｉより前の時点であって鋼材表面の水分付着量ｗｓ（ｔＸｉ）が０になる時点の気温（℃）である。
式（３Ｃ）におけるｐ（ｔＸｉ）は腐食性発現係数であって、前記ｗｓ（ｔＸｉ）が０を超える場合にｐ（ｔＸｉ）＝０とし、前記ｗｓ（ｔＸｉ）が０以下の場合にｐ（ｔＸｉ）＝１、とする。
［５］ 前記鋼材が、ステンレス鋼からなる鋼材であることを特徴とする、［１］乃至［４］の何れか一項に記載の、鋼材の腐食量の予測方法。
［６］ 鋼材を屋外に曝した場合の腐食量を評価する指標として、観測地点Ｘにおける観測時刻ｔＸｉ（ただしｉは１～ｎの整数）毎の風速、風向及び降雨量の気象観測値ＷＸから求まる鋼材の腐食性指標Ｑ（ｔＸｉ，ＷＸ）を用いる、電子計算機による鋼材の腐食量予測システムであって、
前記観測地点Ｘが鋼材の利用予定地点Ｍである気象観測値ＷＭに基づいて、過去の観測時刻ｔＭｉ毎に鋼材の腐食性指標Ｑ（ｔＭｉ，ＷＭ）を求め、観測時刻ｔＭ１～ｔＭｎ間の腐食性指標Ｑ（ｔＭｉ，ＷＭ）の時系列データに対して極値統計解析を行うことにより、将来の評価期間ｔＦ1～ｔＦｎ内における腐食性指標の極大値Ｑｍａｘ（ｔＦｉ）を推測する極値推測部と、
試験地点Ｚにおける鋼材の腐食試験から得た実測腐食量と、前記観測地点Ｘが前記試験地点Ｚである気象観測値ＷＺに基づいて前記腐食試験の試験期間中の観測時刻ｔＺｉ毎に得た鋼材の腐食性指標Ｑ（ｔＺｉ，ＷＺ）の最大値との関係式を予め求めておき、この関係式に、前記極値推測部において推測した腐食性指標の極大値Ｑｍａｘ（ｔＦｉ）を導入して、前記評価期間における鋼材の腐食量の推測値を得る腐食量推測部と、を備えることを特徴とする、鋼材の腐食量予測システム。
なお、上記関係式は、鋼材の腐食試験から得た実測腐食量と、鋼材の腐食性指標Ｑ（ｔＺｉ，ＷＺ）の最大値の対数との関係式としてもよい。
［７］ 前記極値推測部及び前記腐食量推測部にはそれぞれ、前記腐食性指標Ｑ（ｔＸｉ，ＷＸ）を求める計算部が備えられており、
前記計算部は、
前記観測地点Ｘにおける前記観測時刻ｔＸｉ毎の風速ｕ（ｔＸｉ）、風向θｗ（ｔＸｉ）及び降雨量を少なくとも取得する第１処理部と、
前記観測時刻ｔＸｉ毎に、一次腐食性指標Ｑ１（ｔＸｉ，ＷＸ）を下記式（１）により計算する第２処理部と、
前記観測時刻ｔＸｉ毎に、二次腐食性指標Ｑ２（ｔＸｉ，ＷＸ）を下記式（２）により計算する第３処理部と、
前記観測時刻ｔＸｉ毎に、腐食性指標Ｑ（ｔＸｉ，ＷＸ）を下記式（３Ａ）により計算する第４処理部と、を備えている、［６］に記載の、鋼材の腐食量予測システム。
Ｑ１（ｔＸｉ，ＷＸ）＝（ｄ＋１）－０．６｛ｕ（ｔＸｉ）・ｃｏｓ｜θｓ－θｗ（ｔＸｉ）｜｝２ …（１）
Ｑ２（ｔＸｉ，ＷＸ）＝｛Ｑ２（ｔＸｉ－１，ＷＸ）＋Ｑ１（ｔＸｉ，ＷＸ）・ｓ（ｔＸｉ）・（ｔＸｉ－ｔＸｉ－１）｝・ｃ（ｔＸｉ） …（２）
Ｑ（ｔＸｉ，ＷＸ）＝Ｑ２（ｔＸｉ，ＷＸ）・ｐ（ｔＸｉ） …（３Ａ）
ただし、式（１）におけるｄは、前記観測地点Ｘに最も近い海岸と前記観測地点Ｘとの距離（ｍ）であり、ｕ（ｔＸｉ）は、観測時刻ｔＸｉの風速（ｍ／ｓ）であり、θｗ（ｔＸｉ）は北を０°とした場合の観測時刻ｔＸｉでの風向（°）であり、θｓは、前記観測地点Ｘに最も近い海岸からの前記観測地点Ｘの方位であって北を０°とした場合の方位（°）である。
式（２）におけるｃ（ｔＸｉ）は腐食性消失係数であって、観測時刻ｔＸｉ－１～ｔＸｉ間の降雨量が閾値超の場合にｃ（ｔＸｉ）＝０、観測時刻ｔＸｉ－１～ｔＸｉ間の降雨量が閾値以下の場合にｃ（ｔＸｉ）＝１とする。
式（２）におけるｓ（ｔＸｉ）は腐食性蓄積係数であって、観測時刻ｔＸｉ－１～ｔＸｉ間が昼間（日の出～日没までの期間）の場合にｓ（ｔＸｉ）＝０、観測時刻ｔＸｉ－１～ｔＸｉ間が夜間（日没～日の出までの期間）の場合にｓ（ｔＸｉ）＝１とする。
式（３Ａ）におけるｐ（ｔＸｉ）は腐食性発現係数であって、観測時刻ｔＸｉ－１～ｔＸｉ間が夜間（日没～日の出までの期間）の場合にｐ（ｔＸｉ）＝０、観測時刻ｔＸｉ－１～ｔＸｉ間が昼間（日の出～日没までの期間）の場合にｐ（ｔＸｉ）＝１、とする。
［８］ 前記極値推測部及び前記腐食量推測部にはそれぞれ、前記腐食性指標Ｑ（ｔＸｉ，ＷＸ）を求める計算部が備えられており、
前記計算部は、
前記観測地点Ｘにおける前記観測時刻ｔＸｉ毎の風速ｕ（ｔＸｉ）、風向θｗ（ｔＸｉ）及び降雨量を少なくとも取得する第１処理部と、
前記観測時刻ｔＸｉ毎に、一次腐食性指標Ｑ１（ｔＸｉ，ＷＸ）を下記式（１）により計算する第２処理部と、
前記観測時刻ｔＸｉ毎に、二次腐食性指標Ｑ２（ｔＸｉ，ＷＸ）を下記式（２）により計算する第３処理部と、
前記観測時刻ｔＸｉ毎に、腐食性指標Ｑ（ｔＸｉ，ＷＸ）を下記式（３Ｂ）により計算する第４処理部と、を備えている、［６］に記載の、鋼材の腐食量予測システム。
Ｑ１（ｔＸｉ，ＷＸ）＝（ｄ＋１）－０．６｛ｕ（ｔＸｉ）・ｃｏｓ｜θｓ－θｗ（ｔＸｉ）｜｝２ …（１）
Ｑ２（ｔＸｉ，ＷＸ）＝｛Ｑ２（ｔＸｉ－１，ＷＸ）＋Ｑ１（ｔＸｉ，ＷＸ）・ｓ（ｔＸｉ）・（ｔＸｉ－ｔＸｉ－１）｝・ｃ（ｔＸｉ） …（２）
Ｑ（ｔＸｉ，ＷＸ）＝１００．０４Ｔｂ・Ｑ２（ｔＸｉ，ＷＸ）・ｐ（ｔＸｉ） …（３Ｂ）
ただし、式（１）におけるｄは、前記観測地点Ｘに最も近い海岸と前記観測地点Ｘとの距離（ｍ）であり、ｕ（ｔＸｉ）は、観測時刻ｔＸｉの風速（ｍ／ｓ）であり、θｗ（ｔＸｉ）は北を０°とした場合の観測時刻ｔＸｉでの風向（°）であり、θｓは、前記観測地点Ｘに最も近い海岸からの前記観測地点Ｘの方位であって北を０°とした場合の方位（°）である。
式（２）におけるｃ（ｔＸｉ）は腐食性消失係数であって、観測時刻ｔＸｉ－１～ｔＸｉ間の降雨量が閾値超の場合にｃ（ｔＸｉ）＝０、観測時刻ｔＸｉ－１～ｔＸｉ間の降雨量が閾値以下の場合にｃ（ｔＸｉ）＝１とする。
式（２）におけるｓ（ｔＸｉ）は腐食性蓄積係数であって、観測時刻ｔＸｉ－１～ｔＸｉ間が昼間（日の出～日没までの期間）の場合にｓ（ｔＸｉ）＝０、観測時刻ｔＸｉ－１～ｔＸｉ間が夜間（日没～日の出までの期間）の場合にｓ（ｔＸｉ）＝１とする。
式（３Ｂ）におけるＴｂは、日の出時の気温Ｔｅｍｐ（℃）とした場合にＴｅｍｐ（℃）～（Ｔｅｍｐ＋２）℃の範囲の何れかの温度である。
式（３Ｂ）におけるｐ（ｔＸｉ）は腐食性発現係数であって、観測時刻ｔＸｉ－１～ｔＸｉ間が夜間（日没～日の出までの期間）の場合にｐ（ｔＸｉ）＝０、観測時刻ｔＸｉ－１～ｔＸｉ間が昼間（日の出～日没までの期間）の場合にｐ（ｔＸｉ）＝１、とする。
［９］ 前記極値推測部及び前記腐食量推測部にはそれぞれ、前記腐食性指標Ｑ（ｔＸｉ，ＷＸ）を求める計算部が備えられており、
前記計算部は、
前記観測地点Ｘにおける前記観測時刻ｔＸｉ毎の風速ｕ（ｔＸｉ）、風向θｗ（ｔＸｉ）及び降雨量を少なくとも取得する第１処理部と、
前記観測時刻ｔＸｉ毎に、一次腐食性指標Ｑ１（ｔＸｉ，ＷＸ）を下記式（１）により計算する第２処理部と、
前記観測時刻ｔＸｉ毎に、二次腐食性指標Ｑ２（ｔＸｉ，ＷＸ）を下記式（２）により計算する第３処理部と、
前記観測時刻ｔＸｉ毎に、腐食性指標Ｑ（ｔＸｉ，ＷＸ）を下記式（３Ｃ）により計算する第４処理部と、を備えている、［６］に記載の、鋼材の腐食量予測システム。
Ｑ１（ｔＸｉ，ＷＸ）＝（ｄ＋１）－０．６｛ｕ（ｔＸｉ）・ｃｏｓ｜θｓ－θｗ（ｔＸｉ）｜｝２ …（１）
Ｑ２（ｔＸｉ，ＷＸ）＝｛Ｑ２（ｔＸｉ－１，ＷＸ）＋Ｑ１（ｔＸｉ，ＷＸ）・ｓ（ｔＸｉ）・（ｔＸｉ－ｔＸｉ－１）｝・ｃ（ｔＸｉ） …（２）
Ｑ（ｔＸｉ，ＷＸ）＝１００．０４Ｔｃ・Ｑ２（ｔＸｉ，ＷＸ）・ｐ（ｔＸｉ） …（３Ｃ）
ただし、式（１）におけるｄは、前記観測地点Ｘに最も近い海岸と前記観測地点Ｘとの距離（ｍ）であり、ｕ（ｔＸｉ）は、観測時刻ｔＸｉの風速（ｍ／ｓ）であり、θｗ（ｔＸｉ）は北を０°とした場合の観測時刻ｔＸｉでの風向（°）であり、θｓは、前記観測地点Ｘに最も近い海岸からの前記観測地点Ｘの方位であって北を０°とした場合の方位（°）である。
式（２）におけるｃ（ｔＸｉ）は腐食性消失係数であって、観測時刻ｔＸｉ－１～ｔＸｉ間の降雨量が閾値超の場合にｃ（ｔＸｉ）＝０、観測時刻ｔＸｉ－１～ｔＸｉ間の降雨量が閾値以下の場合にｃ（ｔＸｉ）＝１とする。
式（２）におけるｓ（ｔＸｉ）は腐食性蓄積係数であって、観測時刻ｔＸｉにおける鋼材表面の水分付着量ｗｓ（ｔＸｉ）が０の場合にｓ（ｔＸｉ）＝０とし、前記ｗｓ（ｔＸｉ）が０超の場合にｓ（ｔＸｉ）＝１とする。前記ｗｓ（ｔＸｉ）は、前記観測時刻ｔＸｉにおける降雨および結露生成による水分付着量から水の蒸発量を差し引いた推測値であって、鋼材の表面温度、気温（℃）、気圧（ｈＰａ）、鋼材表面における風速（ｍ／ｓ）および相対湿度（％）から求まる値である。
式（３Ｃ）におけるＴｃは、前記観測時刻ｔＸｉより前の時点であって鋼材表面の水分付着量ｗｓ（ｔＸｉ）が０になる時点の気温（℃）である。
式（３Ｃ）におけるｐ（ｔＸｉ）は腐食性発現係数であって、前記ｗｓ（ｔＸｉ）が０を超える場合にｐ（ｔＸｉ）＝０とし、前記ｗｓ（ｔＸｉ）が０以下の場合にｐ（ｔＸｉ）＝１、とする。
［１０］ 前記鋼材が、ステンレス鋼からなる鋼材であることを特徴とする、［６］乃至［９］の何れか一項に記載の、鋼材の腐食量予測システム。
［１１］ 鋼材を屋外に曝した場合の腐食量を評価する指標として、観測地点Ｘにおける観測時刻ｔＸｉ（ただしｉは１～ｎの整数）毎の風速、風向及び降雨量の気象観測値ＷＸから求まる鋼材の腐食性指標Ｑ（ｔＸｉ，ＷＸ）を用いる、電子計算機に利用される鋼材の腐食量予測プログラムであって、
前記観測地点Ｘが鋼材の利用予定地点Ｍである気象観測値ＷＭに基づいて、過去の観測時刻ｔＭｉ毎に鋼材の腐食性指標Ｑ（ｔＭｉ，ＷＭ）を求め、観測時刻ｔＭ１～ｔＭｎ間の腐食性指標Ｑ（ｔＭｉ，ＷＭ）の時系列データに対して極値統計解析を行うことにより、将来の評価期間ｔＦ1～ｔＦｎ内における腐食性指標の極大値Ｑｍａｘ（ｔＦｉ）を推測する極値推測ステップと、
試験地点Ｚにおける鋼材の腐食試験から得た実測腐食量と、前記観測地点Ｘが前記試験地点Ｚである気象観測値ＷＺに基づいて前記腐食試験の試験期間中の観測時刻ｔＺｉ毎に得た鋼材の腐食性指標Ｑ（ｔＺｉ，ＷＺ）の最大値との関係式を予め求めておき、この関係式に、前記極値推測ステップにおいて推測した腐食性指標の極大値Ｑｍａｘ（ｔＦｉ）を導入して、前記評価期間における鋼材の腐食量の推測値を得る腐食量推測ステップと、を備えることを特徴とする、鋼材の腐食量予測プログラム。
なお、上記関係式は、鋼材の腐食試験から得た実測腐食量と、鋼材の腐食性指標Ｑ（ｔＺｉ，ＷＺ）の最大値の対数との関係式としてもよい。
［１２］ 前記極値推測ステップ及び前記腐食量推測ステップにはそれぞれ、前記腐食性指標Ｑ（ｔＸｉ，ＷＸ）を求める計算ステップが備えられており、
前記計算ステップは、
前記観測地点Ｘにおける前記観測時刻ｔＸｉ毎の風速ｕ（ｔＸｉ）、風向θｗ（ｔＸｉ）及び降雨量を少なくとも取得する第１ステップと、
前記観測時刻ｔＸｉ毎に、一次腐食性指標Ｑ１（ｔＸｉ，ＷＸ）を下記式（１）により計算する第２ステップと、
前記観測時刻ｔＸｉ毎に、二次腐食性指標Ｑ２（ｔＸｉ，ＷＸ）を下記式（２）により計算する第３ステップと、
前記観測時刻ｔＸｉ毎に、腐食性指標Ｑ（ｔＸｉ，ＷＸ）を下記式（３Ａ）により計算する第４ステップと、を備えている［１１］に記載の、鋼材の腐食量予測プログラム。
Ｑ１（ｔＸｉ，ＷＸ）＝（ｄ＋１）－０．６｛ｕ（ｔＸｉ）・ｃｏｓ｜θｓ－θｗ（ｔＸｉ）｜｝２ …（１）
Ｑ２（ｔＸｉ，ＷＸ）＝｛Ｑ２（ｔＸｉ－１，ＷＸ）＋Ｑ１（ｔＸｉ，ＷＸ）・ｓ（ｔＸｉ）・（ｔＸｉ－ｔＸｉＳ－１）｝・ｃ（ｔＸｉ） …（２）
Ｑ（ｔＸｉ，ＷＸ）＝Ｑ２（ｔＸｉ，ＷＸ）・ｐ（ｔＸｉ） …（３Ａ）
ただし、式（１）におけるｄは、前記観測地点Ｘに最も近い海岸と前記観測地点Ｘとの距離（ｍ）であり、ｕ（ｔＸｉ）は、観測時刻ｔＸｉの風速（ｍ／ｓ）であり、θｗ（ｔＸｉ）は北を０°とした場合の観測時刻ｔＸｉでの風向（°）であり、θｓは、前記観測地点Ｘに最も近い海岸からの前記観測地点Ｘの方位であって北を０°とした場合の方位（°）である。
式（２）におけるｃ（ｔＸｉ）は腐食性消失係数であって、観測時刻ｔＸｉ－１～ｔＸｉ間の降雨量が閾値超の場合にｃ（ｔＸｉ）＝０、観測時刻ｔＸｉ－１～ｔＸｉ間の降雨量が閾値以下の場合にｃ（ｔＸｉ）＝１とする。
式（２）におけるｓ（ｔＸｉ）は腐食性蓄積係数であって、観測時刻ｔＸｉ－１～ｔＸｉ間が昼間（日の出～日没までの期間）の場合にｓ（ｔＸｉ）＝０、観測時刻ｔＸｉ－１～ｔＸｉ間が夜間（日没～日の出までの期間）の場合にｓ（ｔＸｉ）＝１とする。
式（３Ａ）におけるｐ（ｔＸｉ）は腐食性発現係数であって、観測時刻ｔＸｉ－１～ｔＸｉ間が夜間（日没～日の出までの期間）の場合にｐ（ｔＸｉ）＝０、観測時刻ｔＸｉ－１～ｔＸｉ間が昼間（日の出～日没までの期間）の場合にｐ（ｔＸｉ）＝１、とする。
［１３］ 前記極値推測ステップ及び前記腐食量推測ステップにはそれぞれ、前記腐食性指標Ｑ（ｔＸｉ，ＷＸ）を求める計算ステップが備えられており、
前記計算ステップは、
前記観測地点Ｘにおける前記観測時刻ｔＸｉ毎の風速ｕ（ｔＸｉ）、風向θｗ（ｔＸｉ）及び降雨量を少なくとも取得する第１ステップと、
前記観測時刻ｔＸｉ毎に、一次腐食性指標Ｑ１（ｔＸｉ，ＷＸ）を下記式（１）により計算する第２ステップと、
前記観測時刻ｔＸｉ毎に、二次腐食性指標Ｑ２（ｔＸｉ，ＷＸ）を下記式（２）により計算する第３ステップと、
前記観測時刻ｔＸｉ毎に、腐食性指標Ｑ（ｔＸｉ，ＷＸ）を下記式（３Ｂ）により計算する第４ステップと、を備えている［１１］に記載の、鋼材の腐食量予測プログラム。
Ｑ１（ｔＸｉ，ＷＸ）＝（ｄ＋１）－０．６｛ｕ（ｔＸｉ）・ｃｏｓ｜θｓ－θｗ（ｔＸｉ）｜｝２ …（１）
Ｑ２（ｔＸｉ，ＷＸ）＝｛Ｑ２（ｔＸｉ－１，ＷＸ）＋Ｑ１（ｔＸｉ，ＷＸ）・ｓ（ｔＸｉ）・（ｔＸｉ－ｔＸｉＳ－１）｝・ｃ（ｔＸｉ） …（２）
Ｑ（ｔＸｉ，ＷＸ）＝１００．０４Ｔｂ・Ｑ２（ｔＸｉ，ＷＸ）・ｐ（ｔＸｉ） …（３Ｂ）
ただし、式（１）におけるｄは、前記観測地点Ｘに最も近い海岸と前記観測地点Ｘとの距離（ｍ）であり、ｕ（ｔＸｉ）は、観測時刻ｔＸｉの風速（ｍ／ｓ）であり、θｗ（ｔＸｉ）は北を０°とした場合の観測時刻ｔＸｉでの風向（°）であり、θｓは、前記観測地点Ｘに最も近い海岸からの前記観測地点Ｘの方位であって北を０°とした場合の方位（°）である。
式（２）におけるｃ（ｔＸｉ）は腐食性消失係数であって、観測時刻ｔＸｉ－１～ｔＸｉ間の降雨量が閾値超の場合にｃ（ｔＸｉ）＝０、観測時刻ｔＸｉ－１～ｔＸｉ間の降雨量が閾値以下の場合にｃ（ｔＸｉ）＝１とする。
式（２）におけるｓ（ｔＸｉ）は腐食性蓄積係数であって、観測時刻ｔＸｉ－１～ｔＸｉ間が昼間（日の出～日没までの期間）の場合にｓ（ｔＸｉ）＝０、観測時刻ｔＸｉ－１～ｔＸｉ間が夜間（日没～日の出までの期間）の場合にｓ（ｔＸｉ）＝１とする。
式（３Ｂ）におけるＴｂは、日の出時の気温Ｔｅｍｐ（℃）とした場合にＴｅｍｐ（℃）～（Ｔｅｍｐ＋２）℃の範囲の何れかの温度である。
式（３Ｂ）におけるｐ（ｔＸｉ）は腐食性発現係数であって、観測時刻ｔＸｉ－１～ｔＸｉ間が夜間（日没～日の出までの期間）の場合にｐ（ｔＸｉ）＝０、観測時刻ｔＸｉ－１～ｔＸｉ間が昼間（日の出～日没までの期間）の場合にｐ（ｔＸｉ）＝１、とする。
［１４］ 前記極値推測ステップ及び前記腐食量推測ステップにはそれぞれ、前記腐食性指標Ｑ（ｔＸｉ，ＷＸ）を求める計算ステップが備えられており、
前記計算ステップは、
前記観測地点Ｘにおける前記観測時刻ｔＸｉ毎の風速ｕ（ｔＸｉ）、風向θｗ（ｔＸｉ）及び降雨量を少なくとも取得する第１ステップと、
前記観測時刻ｔＸｉ毎に、一次腐食性指標Ｑ１（ｔＸｉ，ＷＸ）を下記式（１）により計算する第２ステップと、
前記観測時刻ｔＸｉ毎に、二次腐食性指標Ｑ２（ｔＸｉ，ＷＸ）を下記式（２）により計算する第３ステップと、
前記観測時刻ｔＸｉ毎に、腐食性指標Ｑ（ｔＸｉ，ＷＸ）を下記式（３Ｃ）により計算する第４ステップと、を備えている［１１］に記載の、鋼材の腐食量予測プログラム。
Ｑ１（ｔＸｉ，ＷＸ）＝（ｄ＋１）－０．６｛ｕ（ｔＸｉ）・ｃｏｓ｜θｓ－θｗ（ｔＸｉ）｜｝２ …（１）
Ｑ２（ｔＸｉ，ＷＸ）＝｛Ｑ２（ｔＸｉ－１，ＷＸ）＋Ｑ１（ｔＸｉ，ＷＸ）・ｓ（ｔＸｉ）・（ｔＸｉ－ｔＸｉＳ－１）｝・ｃ（ｔＸｉ） …（２）
Ｑ（ｔＸｉ，ＷＸ）＝１００．０４Ｔｃ・Ｑ２（ｔＸｉ，ＷＸ）・ｐ（ｔＸｉ） …（３Ｃ）
ただし、式（１）におけるｄは、前記観測地点Ｘに最も近い海岸と前記観測地点Ｘとの距離（ｍ）であり、ｕ（ｔＸｉ）は、観測時刻ｔＸｉの風速（ｍ／ｓ）であり、θｗ（ｔＸｉ）は北を０°とした場合の観測時刻ｔＸｉでの風向（°）であり、θｓは、前記観測地点Ｘに最も近い海岸からの前記観測地点Ｘの方位であって北を０°とした場合の方位（°）である。
式（２）におけるｃ（ｔＸｉ）は腐食性消失係数であって、観測時刻ｔＸｉ－１～ｔＸｉ間の降雨量が閾値超の場合にｃ（ｔＸｉ）＝０、観測時刻ｔＸｉ－１～ｔＸｉ間の降雨量が閾値以下の場合にｃ（ｔＸｉ）＝１とする。
式（２）におけるｓ（ｔＸｉ）は腐食性蓄積係数であって、観測時刻ｔＸｉにおける鋼材表面の水分付着量ｗｓ（ｔＸｉ）が０の場合にｓ（ｔＸｉ）＝０とし、前記ｗｓ（ｔＸｉ）が０超の場合にｓ（ｔＸｉ）＝１とする。前記ｗｓ（ｔＸｉ）は、前記観測時刻ｔＸｉにおける降雨および結露生成による水分付着量から水の蒸発量を差し引いた推測値であって、鋼材の表面温度、気温（℃）、気圧（ｈＰａ）、鋼材表面における風速（ｍ／ｓ）および相対湿度（％）から求まる値である。
式（３Ｃ）におけるＴｃは、前記観測時刻ｔＸｉより前の時点であって鋼材表面の水分付着量ｗｓ（ｔＸｉ）が０になる時点の気温（℃）である。
式（３Ｃ）におけるｐ（ｔＸｉ）は腐食性発現係数であって、前記ｗｓ（ｔＸｉ）が０を超える場合にｐ（ｔＸｉ）＝０とし、前記ｗｓ（ｔＸｉ）が０以下の場合にｐ（ｔＸｉ）＝１、とする。
［１５］ 前記鋼材が、ステンレス鋼からなる鋼材であることを特徴とする、［１１］乃至［１４］の何れか一項に記載の、鋼材の腐食量予測プログラム。
［１６］ ［１］乃至［５］の何れか一項に記載の鋼材の腐食量の予測方法により、利用予定地点Ｍにおける将来にわたる鋼材の腐食量を予測する予測工程と、
前記予測工程によって得られた腐食量の予測値または腐食量の予測値に対応する鋼材表面の画像を、営業担当者が顧客に提示する工程と、を備えた、鋼材の提案方法。
［１７］ ［１］乃至［５］の何れか一項に記載の鋼材の腐食量の予測方法により、利用予定地点Ｍにおける将来にわたる鋼材の腐食量を、複数種類の鋼材についてそれぞれ予測する予測工程と、
前記予測工程によって得られた鋼材毎の腐食量の予測値または腐食量の予測値に対応する鋼材表面の画像を、営業担当者が顧客に提示する工程と、を備えた、鋼材の提案方法。
［１８］ 前記腐食量の予測値が、少なくとも小数点第一位までのレイティングナンバーである、［１６］または［１７］に記載の鋼材の提案方法。
［１９］ 前記腐食量の予測値に対応する鋼材表面の画像が、少なくとも小数点第一位までのレイティングナンバーに対応する画像である、［１６］または［１７］に記載の鋼材の提案方法。
［１］ 鋼材を屋外に曝した場合の腐食量を評価する指標として、観測地点Ｘにおける観測時刻ｔＸｉ（ただしｉは１～ｎの整数）毎の風速、風向及び降雨量の気象観測値ＷＸから求まる鋼材の腐食性指標Ｑ（ｔＸｉ，ＷＸ）を用いる、鋼材の腐食量を予測する方法であって、
前記観測地点Ｘが鋼材の利用予定地点Ｍである気象観測値ＷＭに基づいて、過去の観測時刻ｔＭｉ毎に鋼材の腐食性指標Ｑ（ｔＭｉ，ＷＭ）を求め、観測時刻ｔＭ１～ｔＭｎ間の腐食性指標Ｑ（ｔＭｉ，ＷＭ）の時系列データに対して極値統計解析を行うことにより、将来の評価期間ｔＦ1～ｔＦｎ内における腐食性指標の極大値Ｑｍａｘ（ｔＦｉ）を推測する極値推測ステップと、
試験地点Ｚにおける鋼材の腐食試験から得た実測腐食量と、前記観測地点Ｘが前記試験地点Ｚである気象観測値ＷＺに基づいて前記腐食試験の試験期間中の観測時刻ｔＺｉ毎に得た鋼材の腐食性指標Ｑ（ｔＺｉ，ＷＺ）の最大値との関係式を予め求めておき、この関係式に、前記極値推測ステップにおいて推測した腐食性指標の極大値Ｑｍａｘ（ｔＦｉ）を導入して、前記評価期間における鋼材の腐食量の推測値を得る腐食量推測ステップと、を備えることを特徴とする、鋼材の腐食量の予測方法。
なお、上記関係式は、鋼材の腐食試験から得た実測腐食量と、鋼材の腐食性指標Ｑ（ｔＺｉ，ＷＺ）の最大値の対数との関係式としてもよい。
［２］ 前記極値推測ステップ及び前記腐食量推測ステップにはそれぞれ、前記腐食性指標Ｑ（ｔＸｉ，ＷＸ）を求める計算ステップが備えられており、
前記計算ステップは、
前記観測地点Ｘにおける前記観測時刻ｔＸｉ毎の風速ｕ（ｔＸｉ）、風向θｗ（ｔＸｉ）及び降雨量を少なくとも取得する第１ステップと、
前記観測時刻ｔＸｉ毎に、一次腐食性指標Ｑ１（ｔＸｉ，ＷＸ）を下記式（１）により計算する第２ステップと、
前記観測時刻ｔＸｉ毎に、二次腐食性指標Ｑ２（ｔＸｉ，ＷＸ）を下記式（２）により計算する第３ステップと、
前記観測時刻ｔＸｉ毎に、腐食性指標Ｑ（ｔＸｉ，ＷＸ）を下記式（３Ａ）により計算する第４ステップと、を備えている、［１］に記載の、鋼材の腐食量の予測方法。
Ｑ１（ｔＸｉ，ＷＸ）＝（ｄ＋１）－０．６｛ｕ（ｔＸｉ）・ｃｏｓ｜θｓ－θｗ（ｔＸｉ）｜｝２ …（１）
Ｑ２（ｔＸｉ，ＷＸ）＝｛Ｑ２（ｔＸｉ－１，ＷＸ）＋Ｑ１（ｔＸｉ，ＷＸ）・ｓ（ｔＸｉ）・（ｔＸｉ－ｔＸｉ－１）｝・ｃ（ｔＸｉ） …（２）
Ｑ（ｔＸｉ，ＷＸ）＝Ｑ２（ｔＸｉ，ＷＸ）・ｐ（ｔＸｉ） …（３Ａ）
ただし、式（１）におけるｄは、前記観測地点Ｘに最も近い海岸と前記観測地点Ｘとの距離（ｍ）であり、ｕ（ｔＸｉ）は、観測時刻ｔＸｉの風速（ｍ／ｓ）であり、θｗ（ｔＸｉ）は北を０°とした場合の観測時刻ｔＸｉでの風向（°）であり、θｓは、前記観測地点Ｘに最も近い海岸からの前記観測地点Ｘの方位であって北を０°とした場合の方位（°）である。
式（２）におけるｃ（ｔＸｉ）は腐食性消失係数であって、観測時刻ｔＸｉ－１～ｔＸｉ間の降雨量が閾値超の場合にｃ（ｔＸｉ）＝０、観測時刻ｔＸｉ－１～ｔＸｉ間の降雨量が閾値以下の場合にｃ（ｔＸｉ）＝１とする。
式（２）におけるｓ（ｔＸｉ）は腐食性蓄積係数であって、観測時刻ｔＸｉ－１～ｔＸｉ間が昼間（日の出～日没までの期間）の場合にｓ（ｔＸｉ）＝０、観測時刻ｔＸｉ－１～ｔＸｉ間が夜間（日没～日の出までの期間）の場合にｓ（ｔＸｉ）＝１とする。
式（３Ａ）におけるｐ（ｔＸｉ）は腐食性発現係数であって、観測時刻ｔＸｉ－１～ｔＸｉ間が夜間（日没～日の出までの期間）の場合にｐ（ｔＸｉ）＝０、観測時刻ｔＸｉ－１～ｔＸｉ間が昼間（日の出～日没までの期間）の場合にｐ（ｔＸｉ）＝１、とする。
［３］ 前記極値推測ステップ及び前記腐食量推測ステップにはそれぞれ、前記腐食性指標Ｑ（ｔＸｉ，ＷＸ）を求める計算ステップが備えられており、
前記計算ステップは、
前記観測地点Ｘにおける前記観測時刻ｔＸｉ毎の風速ｕ（ｔＸｉ）、風向θｗ（ｔＸｉ）及び降雨量を少なくとも取得する第１ステップと、
前記観測時刻ｔＸｉ毎に、一次腐食性指標Ｑ１（ｔＸｉ，ＷＸ）を下記式（１）により計算する第２ステップと、
前記観測時刻ｔＸｉ毎に、二次腐食性指標Ｑ２（ｔＸｉ，ＷＸ）を下記式（２）により計算する第３ステップと、
前記観測時刻ｔＸｉ毎に、腐食性指標Ｑ（ｔＸｉ，ＷＸ）を下記式（３Ｂ）により計算する第４ステップと、を備えている、［１］に記載の、鋼材の腐食量の予測方法。
Ｑ１（ｔＸｉ，ＷＸ）＝（ｄ＋１）－０．６｛ｕ（ｔＸｉ）・ｃｏｓ｜θｓ－θｗ（ｔＸｉ）｜｝２ …（１）
Ｑ２（ｔＸｉ，ＷＸ）＝｛Ｑ２（ｔＸｉ－１，ＷＸ）＋Ｑ１（ｔＸｉ，ＷＸ）・ｓ（ｔＸｉ）・（ｔＸｉ－ｔＸｉ－１）｝・ｃ（ｔＸｉ） …（２）
Ｑ（ｔＸｉ，ＷＸ）＝１００．０４Ｔｂ・Ｑ２（ｔＸｉ，ＷＸ）・ｐ（ｔＸｉ） …（３Ｂ）
ただし、式（１）におけるｄは、前記観測地点Ｘに最も近い海岸と前記観測地点Ｘとの距離（ｍ）であり、ｕ（ｔＸｉ）は、観測時刻ｔＸｉの風速（ｍ／ｓ）であり、θｗ（ｔＸｉ）は北を０°とした場合の観測時刻ｔＸｉでの風向（°）であり、θｓは、前記観測地点Ｘに最も近い海岸からの前記観測地点Ｘの方位であって北を０°とした場合の方位（°）である。
式（２）におけるｃ（ｔＸｉ）は腐食性消失係数であって、観測時刻ｔＸｉ－１～ｔＸｉ間の降雨量が閾値超の場合にｃ（ｔＸｉ）＝０、観測時刻ｔＸｉ－１～ｔＸｉ間の降雨量が閾値以下の場合にｃ（ｔＸｉ）＝１とする。
式（２）におけるｓ（ｔＸｉ）は腐食性蓄積係数であって、観測時刻ｔＸｉ－１～ｔＸｉ間が昼間（日の出～日没までの期間）の場合にｓ（ｔＸｉ）＝０、観測時刻ｔＸｉ－１～ｔＸｉ間が夜間（日没～日の出までの期間）の場合にｓ（ｔＸｉ）＝１とする。
式（３Ｂ）におけるＴｂは、日の出時の気温Ｔｅｍｐ（℃）とした場合にＴｅｍｐ（℃）～（Ｔｅｍｐ＋２）℃の範囲の何れかの温度である。
式（３Ｂ）におけるｐ（ｔＸｉ）は腐食性発現係数であって、観測時刻ｔＸｉ－１～ｔＸｉ間が夜間（日没～日の出までの期間）の場合にｐ（ｔＸｉ）＝０、観測時刻ｔＸｉ－１～ｔＸｉ間が昼間（日の出～日没までの期間）の場合にｐ（ｔＸｉ）＝１、とする。
［４］ 前記極値推測ステップ及び前記腐食量推測ステップにはそれぞれ、前記腐食性指標Ｑ（ｔＸｉ，ＷＸ）を求める計算ステップが備えられており、
前記計算ステップは、
前記観測地点Ｘにおける前記観測時刻ｔＸｉ毎の風速ｕ（ｔＸｉ）、風向θｗ（ｔＸｉ）及び降雨量を少なくとも取得する第１ステップと、
前記観測時刻ｔＸｉ毎に、一次腐食性指標Ｑ１（ｔＸｉ，ＷＸ）を下記式（１）により計算する第２ステップと、
前記観測時刻ｔＸｉ毎に、二次腐食性指標Ｑ２（ｔＸｉ，ＷＸ）を下記式（２）により計算する第３ステップと、
前記観測時刻ｔＸｉ毎に、腐食性指標Ｑ（ｔＸｉ，ＷＸ）を下記式（３Ｃ）により計算する第４ステップと、を備えている、［１］に記載の、鋼材の腐食量の予測方法。
Ｑ１（ｔＸｉ，ＷＸ）＝（ｄ＋１）－０．６｛ｕ（ｔＸｉ）・ｃｏｓ｜θｓ－θｗ（ｔＸｉ）｜｝２ …（１）
Ｑ２（ｔＸｉ，ＷＸ）＝｛Ｑ２（ｔＸｉ－１，ＷＸ）＋Ｑ１（ｔＸｉ，ＷＸ）・ｓ（ｔＸｉ）・（ｔＸｉ－ｔＸｉ－１）｝・ｃ（ｔＸｉ） …（２）
Ｑ（ｔＸｉ，ＷＸ）＝１００．０４Ｔｃ・Ｑ２（ｔＸｉ，ＷＸ）・ｐ（ｔＸｉ） …（３Ｃ）
ただし、式（１）におけるｄは、前記観測地点Ｘに最も近い海岸と前記観測地点Ｘとの距離（ｍ）であり、ｕ（ｔＸｉ）は、観測時刻ｔＸｉの風速（ｍ／ｓ）であり、θｗ（ｔＸｉ）は北を０°とした場合の観測時刻ｔＸｉでの風向（°）であり、θｓは、前記観測地点Ｘに最も近い海岸からの前記観測地点Ｘの方位であって北を０°とした場合の方位（°）である。
式（２）におけるｃ（ｔＸｉ）は腐食性消失係数であって、観測時刻ｔＸｉ－１～ｔＸｉ間の降雨量が閾値超の場合にｃ（ｔＸｉ）＝０、観測時刻ｔＸｉ－１～ｔＸｉ間の降雨量が閾値以下の場合にｃ（ｔＸｉ）＝１とする。
式（２）におけるｓ（ｔＸｉ）は腐食性蓄積係数であって、観測時刻ｔＸｉにおける鋼材表面の水分付着量ｗｓ（ｔＸｉ）が０の場合にｓ（ｔＸｉ）＝０とし、前記ｗｓ（ｔＸｉ）が０超の場合にｓ（ｔＸｉ）＝１とする。前記ｗｓ（ｔＸｉ）は、前記観測時刻ｔＸｉにおける降雨および結露生成による水分付着量から水の蒸発量を差し引いた推測値であって、鋼材の表面温度、気温（℃）、気圧（ｈＰａ）、鋼材表面における風速（ｍ／ｓ）および相対湿度（％）から求まる値である。
式（３Ｃ）におけるＴｃは、前記観測時刻ｔＸｉより前の時点であって鋼材表面の水分付着量ｗｓ（ｔＸｉ）が０になる時点の気温（℃）である。
式（３Ｃ）におけるｐ（ｔＸｉ）は腐食性発現係数であって、前記ｗｓ（ｔＸｉ）が０を超える場合にｐ（ｔＸｉ）＝０とし、前記ｗｓ（ｔＸｉ）が０以下の場合にｐ（ｔＸｉ）＝１、とする。
［５］ 前記鋼材が、ステンレス鋼からなる鋼材であることを特徴とする、［１］乃至［４］の何れか一項に記載の、鋼材の腐食量の予測方法。
［６］ 鋼材を屋外に曝した場合の腐食量を評価する指標として、観測地点Ｘにおける観測時刻ｔＸｉ（ただしｉは１～ｎの整数）毎の風速、風向及び降雨量の気象観測値ＷＸから求まる鋼材の腐食性指標Ｑ（ｔＸｉ，ＷＸ）を用いる、電子計算機による鋼材の腐食量予測システムであって、
前記観測地点Ｘが鋼材の利用予定地点Ｍである気象観測値ＷＭに基づいて、過去の観測時刻ｔＭｉ毎に鋼材の腐食性指標Ｑ（ｔＭｉ，ＷＭ）を求め、観測時刻ｔＭ１～ｔＭｎ間の腐食性指標Ｑ（ｔＭｉ，ＷＭ）の時系列データに対して極値統計解析を行うことにより、将来の評価期間ｔＦ1～ｔＦｎ内における腐食性指標の極大値Ｑｍａｘ（ｔＦｉ）を推測する極値推測部と、
試験地点Ｚにおける鋼材の腐食試験から得た実測腐食量と、前記観測地点Ｘが前記試験地点Ｚである気象観測値ＷＺに基づいて前記腐食試験の試験期間中の観測時刻ｔＺｉ毎に得た鋼材の腐食性指標Ｑ（ｔＺｉ，ＷＺ）の最大値との関係式を予め求めておき、この関係式に、前記極値推測部において推測した腐食性指標の極大値Ｑｍａｘ（ｔＦｉ）を導入して、前記評価期間における鋼材の腐食量の推測値を得る腐食量推測部と、を備えることを特徴とする、鋼材の腐食量予測システム。
なお、上記関係式は、鋼材の腐食試験から得た実測腐食量と、鋼材の腐食性指標Ｑ（ｔＺｉ，ＷＺ）の最大値の対数との関係式としてもよい。
［７］ 前記極値推測部及び前記腐食量推測部にはそれぞれ、前記腐食性指標Ｑ（ｔＸｉ，ＷＸ）を求める計算部が備えられており、
前記計算部は、
前記観測地点Ｘにおける前記観測時刻ｔＸｉ毎の風速ｕ（ｔＸｉ）、風向θｗ（ｔＸｉ）及び降雨量を少なくとも取得する第１処理部と、
前記観測時刻ｔＸｉ毎に、一次腐食性指標Ｑ１（ｔＸｉ，ＷＸ）を下記式（１）により計算する第２処理部と、
前記観測時刻ｔＸｉ毎に、二次腐食性指標Ｑ２（ｔＸｉ，ＷＸ）を下記式（２）により計算する第３処理部と、
前記観測時刻ｔＸｉ毎に、腐食性指標Ｑ（ｔＸｉ，ＷＸ）を下記式（３Ａ）により計算する第４処理部と、を備えている、［６］に記載の、鋼材の腐食量予測システム。
Ｑ１（ｔＸｉ，ＷＸ）＝（ｄ＋１）－０．６｛ｕ（ｔＸｉ）・ｃｏｓ｜θｓ－θｗ（ｔＸｉ）｜｝２ …（１）
Ｑ２（ｔＸｉ，ＷＸ）＝｛Ｑ２（ｔＸｉ－１，ＷＸ）＋Ｑ１（ｔＸｉ，ＷＸ）・ｓ（ｔＸｉ）・（ｔＸｉ－ｔＸｉ－１）｝・ｃ（ｔＸｉ） …（２）
Ｑ（ｔＸｉ，ＷＸ）＝Ｑ２（ｔＸｉ，ＷＸ）・ｐ（ｔＸｉ） …（３Ａ）
ただし、式（１）におけるｄは、前記観測地点Ｘに最も近い海岸と前記観測地点Ｘとの距離（ｍ）であり、ｕ（ｔＸｉ）は、観測時刻ｔＸｉの風速（ｍ／ｓ）であり、θｗ（ｔＸｉ）は北を０°とした場合の観測時刻ｔＸｉでの風向（°）であり、θｓは、前記観測地点Ｘに最も近い海岸からの前記観測地点Ｘの方位であって北を０°とした場合の方位（°）である。
式（２）におけるｃ（ｔＸｉ）は腐食性消失係数であって、観測時刻ｔＸｉ－１～ｔＸｉ間の降雨量が閾値超の場合にｃ（ｔＸｉ）＝０、観測時刻ｔＸｉ－１～ｔＸｉ間の降雨量が閾値以下の場合にｃ（ｔＸｉ）＝１とする。
式（２）におけるｓ（ｔＸｉ）は腐食性蓄積係数であって、観測時刻ｔＸｉ－１～ｔＸｉ間が昼間（日の出～日没までの期間）の場合にｓ（ｔＸｉ）＝０、観測時刻ｔＸｉ－１～ｔＸｉ間が夜間（日没～日の出までの期間）の場合にｓ（ｔＸｉ）＝１とする。
式（３Ａ）におけるｐ（ｔＸｉ）は腐食性発現係数であって、観測時刻ｔＸｉ－１～ｔＸｉ間が夜間（日没～日の出までの期間）の場合にｐ（ｔＸｉ）＝０、観測時刻ｔＸｉ－１～ｔＸｉ間が昼間（日の出～日没までの期間）の場合にｐ（ｔＸｉ）＝１、とする。
［８］ 前記極値推測部及び前記腐食量推測部にはそれぞれ、前記腐食性指標Ｑ（ｔＸｉ，ＷＸ）を求める計算部が備えられており、
前記計算部は、
前記観測地点Ｘにおける前記観測時刻ｔＸｉ毎の風速ｕ（ｔＸｉ）、風向θｗ（ｔＸｉ）及び降雨量を少なくとも取得する第１処理部と、
前記観測時刻ｔＸｉ毎に、一次腐食性指標Ｑ１（ｔＸｉ，ＷＸ）を下記式（１）により計算する第２処理部と、
前記観測時刻ｔＸｉ毎に、二次腐食性指標Ｑ２（ｔＸｉ，ＷＸ）を下記式（２）により計算する第３処理部と、
前記観測時刻ｔＸｉ毎に、腐食性指標Ｑ（ｔＸｉ，ＷＸ）を下記式（３Ｂ）により計算する第４処理部と、を備えている、［６］に記載の、鋼材の腐食量予測システム。
Ｑ１（ｔＸｉ，ＷＸ）＝（ｄ＋１）－０．６｛ｕ（ｔＸｉ）・ｃｏｓ｜θｓ－θｗ（ｔＸｉ）｜｝２ …（１）
Ｑ２（ｔＸｉ，ＷＸ）＝｛Ｑ２（ｔＸｉ－１，ＷＸ）＋Ｑ１（ｔＸｉ，ＷＸ）・ｓ（ｔＸｉ）・（ｔＸｉ－ｔＸｉ－１）｝・ｃ（ｔＸｉ） …（２）
Ｑ（ｔＸｉ，ＷＸ）＝１００．０４Ｔｂ・Ｑ２（ｔＸｉ，ＷＸ）・ｐ（ｔＸｉ） …（３Ｂ）
ただし、式（１）におけるｄは、前記観測地点Ｘに最も近い海岸と前記観測地点Ｘとの距離（ｍ）であり、ｕ（ｔＸｉ）は、観測時刻ｔＸｉの風速（ｍ／ｓ）であり、θｗ（ｔＸｉ）は北を０°とした場合の観測時刻ｔＸｉでの風向（°）であり、θｓは、前記観測地点Ｘに最も近い海岸からの前記観測地点Ｘの方位であって北を０°とした場合の方位（°）である。
式（２）におけるｃ（ｔＸｉ）は腐食性消失係数であって、観測時刻ｔＸｉ－１～ｔＸｉ間の降雨量が閾値超の場合にｃ（ｔＸｉ）＝０、観測時刻ｔＸｉ－１～ｔＸｉ間の降雨量が閾値以下の場合にｃ（ｔＸｉ）＝１とする。
式（２）におけるｓ（ｔＸｉ）は腐食性蓄積係数であって、観測時刻ｔＸｉ－１～ｔＸｉ間が昼間（日の出～日没までの期間）の場合にｓ（ｔＸｉ）＝０、観測時刻ｔＸｉ－１～ｔＸｉ間が夜間（日没～日の出までの期間）の場合にｓ（ｔＸｉ）＝１とする。
式（３Ｂ）におけるＴｂは、日の出時の気温Ｔｅｍｐ（℃）とした場合にＴｅｍｐ（℃）～（Ｔｅｍｐ＋２）℃の範囲の何れかの温度である。
式（３Ｂ）におけるｐ（ｔＸｉ）は腐食性発現係数であって、観測時刻ｔＸｉ－１～ｔＸｉ間が夜間（日没～日の出までの期間）の場合にｐ（ｔＸｉ）＝０、観測時刻ｔＸｉ－１～ｔＸｉ間が昼間（日の出～日没までの期間）の場合にｐ（ｔＸｉ）＝１、とする。
［９］ 前記極値推測部及び前記腐食量推測部にはそれぞれ、前記腐食性指標Ｑ（ｔＸｉ，ＷＸ）を求める計算部が備えられており、
前記計算部は、
前記観測地点Ｘにおける前記観測時刻ｔＸｉ毎の風速ｕ（ｔＸｉ）、風向θｗ（ｔＸｉ）及び降雨量を少なくとも取得する第１処理部と、
前記観測時刻ｔＸｉ毎に、一次腐食性指標Ｑ１（ｔＸｉ，ＷＸ）を下記式（１）により計算する第２処理部と、
前記観測時刻ｔＸｉ毎に、二次腐食性指標Ｑ２（ｔＸｉ，ＷＸ）を下記式（２）により計算する第３処理部と、
前記観測時刻ｔＸｉ毎に、腐食性指標Ｑ（ｔＸｉ，ＷＸ）を下記式（３Ｃ）により計算する第４処理部と、を備えている、［６］に記載の、鋼材の腐食量予測システム。
Ｑ１（ｔＸｉ，ＷＸ）＝（ｄ＋１）－０．６｛ｕ（ｔＸｉ）・ｃｏｓ｜θｓ－θｗ（ｔＸｉ）｜｝２ …（１）
Ｑ２（ｔＸｉ，ＷＸ）＝｛Ｑ２（ｔＸｉ－１，ＷＸ）＋Ｑ１（ｔＸｉ，ＷＸ）・ｓ（ｔＸｉ）・（ｔＸｉ－ｔＸｉ－１）｝・ｃ（ｔＸｉ） …（２）
Ｑ（ｔＸｉ，ＷＸ）＝１００．０４Ｔｃ・Ｑ２（ｔＸｉ，ＷＸ）・ｐ（ｔＸｉ） …（３Ｃ）
ただし、式（１）におけるｄは、前記観測地点Ｘに最も近い海岸と前記観測地点Ｘとの距離（ｍ）であり、ｕ（ｔＸｉ）は、観測時刻ｔＸｉの風速（ｍ／ｓ）であり、θｗ（ｔＸｉ）は北を０°とした場合の観測時刻ｔＸｉでの風向（°）であり、θｓは、前記観測地点Ｘに最も近い海岸からの前記観測地点Ｘの方位であって北を０°とした場合の方位（°）である。
式（２）におけるｃ（ｔＸｉ）は腐食性消失係数であって、観測時刻ｔＸｉ－１～ｔＸｉ間の降雨量が閾値超の場合にｃ（ｔＸｉ）＝０、観測時刻ｔＸｉ－１～ｔＸｉ間の降雨量が閾値以下の場合にｃ（ｔＸｉ）＝１とする。
式（２）におけるｓ（ｔＸｉ）は腐食性蓄積係数であって、観測時刻ｔＸｉにおける鋼材表面の水分付着量ｗｓ（ｔＸｉ）が０の場合にｓ（ｔＸｉ）＝０とし、前記ｗｓ（ｔＸｉ）が０超の場合にｓ（ｔＸｉ）＝１とする。前記ｗｓ（ｔＸｉ）は、前記観測時刻ｔＸｉにおける降雨および結露生成による水分付着量から水の蒸発量を差し引いた推測値であって、鋼材の表面温度、気温（℃）、気圧（ｈＰａ）、鋼材表面における風速（ｍ／ｓ）および相対湿度（％）から求まる値である。
式（３Ｃ）におけるＴｃは、前記観測時刻ｔＸｉより前の時点であって鋼材表面の水分付着量ｗｓ（ｔＸｉ）が０になる時点の気温（℃）である。
式（３Ｃ）におけるｐ（ｔＸｉ）は腐食性発現係数であって、前記ｗｓ（ｔＸｉ）が０を超える場合にｐ（ｔＸｉ）＝０とし、前記ｗｓ（ｔＸｉ）が０以下の場合にｐ（ｔＸｉ）＝１、とする。
［１０］ 前記鋼材が、ステンレス鋼からなる鋼材であることを特徴とする、［６］乃至［９］の何れか一項に記載の、鋼材の腐食量予測システム。
［１１］ 鋼材を屋外に曝した場合の腐食量を評価する指標として、観測地点Ｘにおける観測時刻ｔＸｉ（ただしｉは１～ｎの整数）毎の風速、風向及び降雨量の気象観測値ＷＸから求まる鋼材の腐食性指標Ｑ（ｔＸｉ，ＷＸ）を用いる、電子計算機に利用される鋼材の腐食量予測プログラムであって、
前記観測地点Ｘが鋼材の利用予定地点Ｍである気象観測値ＷＭに基づいて、過去の観測時刻ｔＭｉ毎に鋼材の腐食性指標Ｑ（ｔＭｉ，ＷＭ）を求め、観測時刻ｔＭ１～ｔＭｎ間の腐食性指標Ｑ（ｔＭｉ，ＷＭ）の時系列データに対して極値統計解析を行うことにより、将来の評価期間ｔＦ1～ｔＦｎ内における腐食性指標の極大値Ｑｍａｘ（ｔＦｉ）を推測する極値推測ステップと、
試験地点Ｚにおける鋼材の腐食試験から得た実測腐食量と、前記観測地点Ｘが前記試験地点Ｚである気象観測値ＷＺに基づいて前記腐食試験の試験期間中の観測時刻ｔＺｉ毎に得た鋼材の腐食性指標Ｑ（ｔＺｉ，ＷＺ）の最大値との関係式を予め求めておき、この関係式に、前記極値推測ステップにおいて推測した腐食性指標の極大値Ｑｍａｘ（ｔＦｉ）を導入して、前記評価期間における鋼材の腐食量の推測値を得る腐食量推測ステップと、を備えることを特徴とする、鋼材の腐食量予測プログラム。
なお、上記関係式は、鋼材の腐食試験から得た実測腐食量と、鋼材の腐食性指標Ｑ（ｔＺｉ，ＷＺ）の最大値の対数との関係式としてもよい。
［１２］ 前記極値推測ステップ及び前記腐食量推測ステップにはそれぞれ、前記腐食性指標Ｑ（ｔＸｉ，ＷＸ）を求める計算ステップが備えられており、
前記計算ステップは、
前記観測地点Ｘにおける前記観測時刻ｔＸｉ毎の風速ｕ（ｔＸｉ）、風向θｗ（ｔＸｉ）及び降雨量を少なくとも取得する第１ステップと、
前記観測時刻ｔＸｉ毎に、一次腐食性指標Ｑ１（ｔＸｉ，ＷＸ）を下記式（１）により計算する第２ステップと、
前記観測時刻ｔＸｉ毎に、二次腐食性指標Ｑ２（ｔＸｉ，ＷＸ）を下記式（２）により計算する第３ステップと、
前記観測時刻ｔＸｉ毎に、腐食性指標Ｑ（ｔＸｉ，ＷＸ）を下記式（３Ａ）により計算する第４ステップと、を備えている［１１］に記載の、鋼材の腐食量予測プログラム。
Ｑ１（ｔＸｉ，ＷＸ）＝（ｄ＋１）－０．６｛ｕ（ｔＸｉ）・ｃｏｓ｜θｓ－θｗ（ｔＸｉ）｜｝２ …（１）
Ｑ２（ｔＸｉ，ＷＸ）＝｛Ｑ２（ｔＸｉ－１，ＷＸ）＋Ｑ１（ｔＸｉ，ＷＸ）・ｓ（ｔＸｉ）・（ｔＸｉ－ｔＸｉＳ－１）｝・ｃ（ｔＸｉ） …（２）
Ｑ（ｔＸｉ，ＷＸ）＝Ｑ２（ｔＸｉ，ＷＸ）・ｐ（ｔＸｉ） …（３Ａ）
ただし、式（１）におけるｄは、前記観測地点Ｘに最も近い海岸と前記観測地点Ｘとの距離（ｍ）であり、ｕ（ｔＸｉ）は、観測時刻ｔＸｉの風速（ｍ／ｓ）であり、θｗ（ｔＸｉ）は北を０°とした場合の観測時刻ｔＸｉでの風向（°）であり、θｓは、前記観測地点Ｘに最も近い海岸からの前記観測地点Ｘの方位であって北を０°とした場合の方位（°）である。
式（２）におけるｃ（ｔＸｉ）は腐食性消失係数であって、観測時刻ｔＸｉ－１～ｔＸｉ間の降雨量が閾値超の場合にｃ（ｔＸｉ）＝０、観測時刻ｔＸｉ－１～ｔＸｉ間の降雨量が閾値以下の場合にｃ（ｔＸｉ）＝１とする。
式（２）におけるｓ（ｔＸｉ）は腐食性蓄積係数であって、観測時刻ｔＸｉ－１～ｔＸｉ間が昼間（日の出～日没までの期間）の場合にｓ（ｔＸｉ）＝０、観測時刻ｔＸｉ－１～ｔＸｉ間が夜間（日没～日の出までの期間）の場合にｓ（ｔＸｉ）＝１とする。
式（３Ａ）におけるｐ（ｔＸｉ）は腐食性発現係数であって、観測時刻ｔＸｉ－１～ｔＸｉ間が夜間（日没～日の出までの期間）の場合にｐ（ｔＸｉ）＝０、観測時刻ｔＸｉ－１～ｔＸｉ間が昼間（日の出～日没までの期間）の場合にｐ（ｔＸｉ）＝１、とする。
［１３］ 前記極値推測ステップ及び前記腐食量推測ステップにはそれぞれ、前記腐食性指標Ｑ（ｔＸｉ，ＷＸ）を求める計算ステップが備えられており、
前記計算ステップは、
前記観測地点Ｘにおける前記観測時刻ｔＸｉ毎の風速ｕ（ｔＸｉ）、風向θｗ（ｔＸｉ）及び降雨量を少なくとも取得する第１ステップと、
前記観測時刻ｔＸｉ毎に、一次腐食性指標Ｑ１（ｔＸｉ，ＷＸ）を下記式（１）により計算する第２ステップと、
前記観測時刻ｔＸｉ毎に、二次腐食性指標Ｑ２（ｔＸｉ，ＷＸ）を下記式（２）により計算する第３ステップと、
前記観測時刻ｔＸｉ毎に、腐食性指標Ｑ（ｔＸｉ，ＷＸ）を下記式（３Ｂ）により計算する第４ステップと、を備えている［１１］に記載の、鋼材の腐食量予測プログラム。
Ｑ１（ｔＸｉ，ＷＸ）＝（ｄ＋１）－０．６｛ｕ（ｔＸｉ）・ｃｏｓ｜θｓ－θｗ（ｔＸｉ）｜｝２ …（１）
Ｑ２（ｔＸｉ，ＷＸ）＝｛Ｑ２（ｔＸｉ－１，ＷＸ）＋Ｑ１（ｔＸｉ，ＷＸ）・ｓ（ｔＸｉ）・（ｔＸｉ－ｔＸｉＳ－１）｝・ｃ（ｔＸｉ） …（２）
Ｑ（ｔＸｉ，ＷＸ）＝１００．０４Ｔｂ・Ｑ２（ｔＸｉ，ＷＸ）・ｐ（ｔＸｉ） …（３Ｂ）
ただし、式（１）におけるｄは、前記観測地点Ｘに最も近い海岸と前記観測地点Ｘとの距離（ｍ）であり、ｕ（ｔＸｉ）は、観測時刻ｔＸｉの風速（ｍ／ｓ）であり、θｗ（ｔＸｉ）は北を０°とした場合の観測時刻ｔＸｉでの風向（°）であり、θｓは、前記観測地点Ｘに最も近い海岸からの前記観測地点Ｘの方位であって北を０°とした場合の方位（°）である。
式（２）におけるｃ（ｔＸｉ）は腐食性消失係数であって、観測時刻ｔＸｉ－１～ｔＸｉ間の降雨量が閾値超の場合にｃ（ｔＸｉ）＝０、観測時刻ｔＸｉ－１～ｔＸｉ間の降雨量が閾値以下の場合にｃ（ｔＸｉ）＝１とする。
式（２）におけるｓ（ｔＸｉ）は腐食性蓄積係数であって、観測時刻ｔＸｉ－１～ｔＸｉ間が昼間（日の出～日没までの期間）の場合にｓ（ｔＸｉ）＝０、観測時刻ｔＸｉ－１～ｔＸｉ間が夜間（日没～日の出までの期間）の場合にｓ（ｔＸｉ）＝１とする。
式（３Ｂ）におけるＴｂは、日の出時の気温Ｔｅｍｐ（℃）とした場合にＴｅｍｐ（℃）～（Ｔｅｍｐ＋２）℃の範囲の何れかの温度である。
式（３Ｂ）におけるｐ（ｔＸｉ）は腐食性発現係数であって、観測時刻ｔＸｉ－１～ｔＸｉ間が夜間（日没～日の出までの期間）の場合にｐ（ｔＸｉ）＝０、観測時刻ｔＸｉ－１～ｔＸｉ間が昼間（日の出～日没までの期間）の場合にｐ（ｔＸｉ）＝１、とする。
［１４］ 前記極値推測ステップ及び前記腐食量推測ステップにはそれぞれ、前記腐食性指標Ｑ（ｔＸｉ，ＷＸ）を求める計算ステップが備えられており、
前記計算ステップは、
前記観測地点Ｘにおける前記観測時刻ｔＸｉ毎の風速ｕ（ｔＸｉ）、風向θｗ（ｔＸｉ）及び降雨量を少なくとも取得する第１ステップと、
前記観測時刻ｔＸｉ毎に、一次腐食性指標Ｑ１（ｔＸｉ，ＷＸ）を下記式（１）により計算する第２ステップと、
前記観測時刻ｔＸｉ毎に、二次腐食性指標Ｑ２（ｔＸｉ，ＷＸ）を下記式（２）により計算する第３ステップと、
前記観測時刻ｔＸｉ毎に、腐食性指標Ｑ（ｔＸｉ，ＷＸ）を下記式（３Ｃ）により計算する第４ステップと、を備えている［１１］に記載の、鋼材の腐食量予測プログラム。
Ｑ１（ｔＸｉ，ＷＸ）＝（ｄ＋１）－０．６｛ｕ（ｔＸｉ）・ｃｏｓ｜θｓ－θｗ（ｔＸｉ）｜｝２ …（１）
Ｑ２（ｔＸｉ，ＷＸ）＝｛Ｑ２（ｔＸｉ－１，ＷＸ）＋Ｑ１（ｔＸｉ，ＷＸ）・ｓ（ｔＸｉ）・（ｔＸｉ－ｔＸｉＳ－１）｝・ｃ（ｔＸｉ） …（２）
Ｑ（ｔＸｉ，ＷＸ）＝１００．０４Ｔｃ・Ｑ２（ｔＸｉ，ＷＸ）・ｐ（ｔＸｉ） …（３Ｃ）
ただし、式（１）におけるｄは、前記観測地点Ｘに最も近い海岸と前記観測地点Ｘとの距離（ｍ）であり、ｕ（ｔＸｉ）は、観測時刻ｔＸｉの風速（ｍ／ｓ）であり、θｗ（ｔＸｉ）は北を０°とした場合の観測時刻ｔＸｉでの風向（°）であり、θｓは、前記観測地点Ｘに最も近い海岸からの前記観測地点Ｘの方位であって北を０°とした場合の方位（°）である。
式（２）におけるｃ（ｔＸｉ）は腐食性消失係数であって、観測時刻ｔＸｉ－１～ｔＸｉ間の降雨量が閾値超の場合にｃ（ｔＸｉ）＝０、観測時刻ｔＸｉ－１～ｔＸｉ間の降雨量が閾値以下の場合にｃ（ｔＸｉ）＝１とする。
式（２）におけるｓ（ｔＸｉ）は腐食性蓄積係数であって、観測時刻ｔＸｉにおける鋼材表面の水分付着量ｗｓ（ｔＸｉ）が０の場合にｓ（ｔＸｉ）＝０とし、前記ｗｓ（ｔＸｉ）が０超の場合にｓ（ｔＸｉ）＝１とする。前記ｗｓ（ｔＸｉ）は、前記観測時刻ｔＸｉにおける降雨および結露生成による水分付着量から水の蒸発量を差し引いた推測値であって、鋼材の表面温度、気温（℃）、気圧（ｈＰａ）、鋼材表面における風速（ｍ／ｓ）および相対湿度（％）から求まる値である。
式（３Ｃ）におけるＴｃは、前記観測時刻ｔＸｉより前の時点であって鋼材表面の水分付着量ｗｓ（ｔＸｉ）が０になる時点の気温（℃）である。
式（３Ｃ）におけるｐ（ｔＸｉ）は腐食性発現係数であって、前記ｗｓ（ｔＸｉ）が０を超える場合にｐ（ｔＸｉ）＝０とし、前記ｗｓ（ｔＸｉ）が０以下の場合にｐ（ｔＸｉ）＝１、とする。
［１５］ 前記鋼材が、ステンレス鋼からなる鋼材であることを特徴とする、［１１］乃至［１４］の何れか一項に記載の、鋼材の腐食量予測プログラム。
［１６］ ［１］乃至［５］の何れか一項に記載の鋼材の腐食量の予測方法により、利用予定地点Ｍにおける将来にわたる鋼材の腐食量を予測する予測工程と、
前記予測工程によって得られた腐食量の予測値または腐食量の予測値に対応する鋼材表面の画像を、営業担当者が顧客に提示する工程と、を備えた、鋼材の提案方法。
［１７］ ［１］乃至［５］の何れか一項に記載の鋼材の腐食量の予測方法により、利用予定地点Ｍにおける将来にわたる鋼材の腐食量を、複数種類の鋼材についてそれぞれ予測する予測工程と、
前記予測工程によって得られた鋼材毎の腐食量の予測値または腐食量の予測値に対応する鋼材表面の画像を、営業担当者が顧客に提示する工程と、を備えた、鋼材の提案方法。
［１８］ 前記腐食量の予測値が、少なくとも小数点第一位までのレイティングナンバーである、［１６］または［１７］に記載の鋼材の提案方法。
［１９］ 前記腐食量の予測値に対応する鋼材表面の画像が、少なくとも小数点第一位までのレイティングナンバーに対応する画像である、［１６］または［１７］に記載の鋼材の提案方法。
【発明の効果】
【0009】
本発明の鋼材の腐食量の予測方法、鋼材の腐食量予測システム及び鋼材の腐食量予測プログラムでは、鋼材の利用予定地点Ｍにおける過去の観測時刻ｔＭｉ毎の腐食性指標ＱＭ（ｔＭｉ，ＷＭ）の時系列データに対して極値統計解析を行うことで、将来の評価期間ｔＦ1～ｔＦｎ内における腐食性指標の極大値Ｑｍａｘ（ｔＦｉ）を推測する。次いで、腐食試験の結果から予め求めておいた鋼材の実測腐食量と鋼材の腐食性指標Ｑ（ｔＺｉ，ＷＺ）との関係式に、腐食性指標の極大値Ｑｍａｘ（ｔＦｉ）を導入する。本発明は、このようにして将来の鋼材の腐食量を推測するので、従来の年平均値に基づいて腐食量を予測する方法に比べて、腐食量の予測精度を高めることができる。
また、本発明の鋼材の腐食量の予測方法、鋼材の腐食量予測システム及び鋼材の腐食量予測プログラムは、試験地点Ｚにおける鋼材の腐食試験の結果に基づき、利用予定地点Ｍにおける鋼材の将来の腐食量を推測するものである。本発明では、試験地点Ｚにおける鋼材の実測腐食量と腐食性指標Ｑ（ｔＺｉ，ＷＺ）との対応関係と、利用予定地点Ｍの気象観測値ＷＭから得た腐食性指標Ｑ（ｔＭｉ，ＷＭ）の時系列データとが得られていればよい。よって、本発明は、汎用性に優れる。
また、本発明の鋼材の腐食量の予測方法、鋼材の腐食量予測システム及び鋼材の腐食量予測プログラムは、試験地点Ｚにおける鋼材の腐食試験の結果に基づき、利用予定地点Ｍにおける鋼材の将来の腐食量を推測するものである。本発明では、試験地点Ｚにおける鋼材の実測腐食量と腐食性指標Ｑ（ｔＺｉ，ＷＺ）との対応関係と、利用予定地点Ｍの気象観測値ＷＭから得た腐食性指標Ｑ（ｔＭｉ，ＷＭ）の時系列データとが得られていればよい。よって、本発明は、汎用性に優れる。
【0010】
また、本発明の鋼材の腐食量の予測方法、鋼材の腐食量予測システム及び鋼材の腐食量予測プログラムによれば、一次腐食性指標Ｑ１（ｔＸｉ，ＷＸ）及び二次腐食性指標Ｑ２（ｔＸｉ，ＷＸ）を順次求め、更に、二次腐食性指標Ｑ２（ｔＸｉ，ＷＸ）から腐食性指標Ｑ（ｔＸｉ，ＷＸ）を求める。ここで、一次腐食性指標Ｑ１（ｔＸｉ，ＷＸ）は、鋼材への塩分付着量の指標となるものであり、風速ｕ（ｔＸｉ）及び風向θｓ、θｗ（ｔＸｉ）に基づいて算出される。また、二次腐食性指標Ｑ２（ｔＸｉ，ＷＸ）は、観測時刻ｔＸｉの直前の観測時刻である観測時刻ｔＸｉ－１における二次腐食性指標Ｑ２（ｔＸｉ－１，ＷＸ）に、観測時刻ｔＸｉにおける一次腐食性指標Ｑ１（ｔＸｉ，ＷＸ）を加算し、更に腐食の主原因である塩分が降雨により洗い流されることを腐食性消失係数ｃ（ｔＸｉ）として導入するとともに、鋼材表面に結露水が発生しやすい夜間に塩分の付着が増大することを腐食性蓄積係数ｓ（ｔＸｉ）として導入することにより算出される。そして、腐食性指標Ｑ（ｔＸｉ，ＷＸ）は、夜間に付着した塩分を含む結露水が昼間の日照により乾燥されることで腐食が進展する事象を腐食性発現係数ｐ（ｔＸｉ）により数値化し、この腐食性発現係数ｐ（ｔＸｉ）を二次腐食性指標Ｑ２（ｔＸｉ，ＷＸ）に乗じることで算出される。これにより、鋼材の塩分付着量に影響を与える気象状況が反映された腐食性指標Ｑ（ｔＸｉ，ＷＸ）を得ることができ、これを用いることで、腐食量の予測精度を高めることができる。
【0011】
本発明の鋼材の提案方法は、本発明の鋼材の腐食量の予測方法により、利用予定地点Ｍにおける将来にわたる鋼材の腐食量を予測し、得られた腐食量の予測値を、営業担当者が顧客に提示する。これにより本発明は、鋼材を長期間に渡り屋外環境下で使用した場合の将来の腐食量の予測値を知りたいという顧客の要望に応えることができる。また、本発明は、精度の高い腐食量の予測値を提供できるので、適切な鋼材を顧客に勧めることができる。
【0012】
また、本発明の鋼材の提案方法によれば、本発明の鋼材の腐食量の予測方法により、利用予定地点Ｍにおける将来にわたる鋼材の腐食量を、複数種類の鋼材についてそれぞれ予測し、得られた鋼材毎の腐食量の予測値を、営業担当者が顧客に提示する。これにより本発明は、鋼材を長期間に渡り屋外環境下で使用した場合の将来の腐食量の予測値を知りたいという顧客の要望に応えることができる。また、本発明は、精度の高い腐食量の予測値を鋼種毎に提供できるので、適切な鋼材選択して顧客に勧めることができる。
【図面の簡単な説明】
【0013】
【図1】本発明の実施形態である鋼材の腐食量予測システムを説明する模式図。
【図2】本発明の実施形態である鋼材の腐食量予測プログラムを説明する模式図。
【図3】腐食性指標Ｑ（ｔＭｉ，ＷＭ）の経時変化の一例を示すグラフ。
【図4A】観測時刻ｔＭ１～ｔＭｎ間の腐食性指標Ｑ（ｔＭｉ，ＷＭ）の時系列データの一例を示すグラフ。
【図4B】将来の評価期間ｔＦ1～ｔＦｎ内における腐食性指標の極大値Ｑｍａｘ（ｔＦｉ）を示すグラフ。
【図5】ＳＵＳ３０４についての最大値ＱＲｍａｘの対数と実測腐食量（ＲＮ）との関係を示すグラフ。
【図6】本実施形態の腐食量の推測方法のアウトプットの一例を示すグラフ。
【図7】ＳＵＳ３０４についての最大値ＱＲｍａｘの対数と実測腐食量（ＲＮ）との関係を示すグラフ。
【図8】ＳＵＳ３０４についての温度補正後の最大値ＱＲｍａｘの対数と実測腐食量（ＲＮ）との関係を示すグラフ。
【発明を実施するための形態】
【0014】
屋外環境下に置かれた鋼材は、風の影響などにより、例えば塩分のような腐食を促進させる物質が鋼材表面に付着することで、徐々に腐食が進行する。ビル等の一般的な建築物や、橋梁、塔などの土木構造物は、屋外環境下において長期間にわたる耐久性が要求される。従って、これら建築物や構造物に使用される鋼材は、将来の腐食量の進行度を予測できることが好ましい。従来技術として、年間ぬれ時間、年平均風速、年平均気温等に基づいて、腐食量を予測する方法が知られていた。
【0015】
しかし、鋼材の腐食は、塩分が付着することで急激に腐食が進展する場合があり、突発的な気象変動に大きく影響される場合がある。そのため、年間平均値に基づいた従来の予測では、データの構造上、鋼材の腐食に大きな影響を与え得る短期間の突発的な気象変動が考慮されない。よって、従来の予測方法よりも腐食量の予測精度を向上させる余地があった。
【0016】
すなわち、屋外環境下に置かれた鋼材の腐食量は、鋼材に付着する塩分量に影響される。塩分は、主に、海岸方向から吹く風に運ばれて鋼材に付着するので、従来は、風速の月間平均値や年平均値並びに風向の推移から、腐食量をある程度予測可能と考えられていた。しかし、実際に屋外環境下に鋼材を放置して腐食試験を行うと、風速の月間平均値、年平均値および風向を考慮するだけでは、腐食量を説明できない場合があった。
【0017】
また、ステンレス鋼からなる鋼材は、普通鋼などに比べて塩分の付着による腐食がより急激に進む場合がある。更に、ステンレス鋼は、意匠性に優れることから、未塗装のままで建物の外壁等に利用されることが多く、普通鋼などに比べて塩分の影響を受けやすい。
【0018】
よって、鋼材の腐食量を精度よく予測できる手法の開発が望まれている。
【0019】
そこで、発明者らが鋭意検討したところ、鋼材を屋外に放置した際の腐食が進展しやすい気象条件として、夜間から昼間に至る間で降雨が起きない気象条件が見出された。即ち、夜間の気温の低下に伴い鋼材表面に結露水が付着すると、そこに海岸方向から吹く海風によって運ばれた塩分が結露水中に溶け込む。そして、この塩分を含む結露水が昼間の日照によって蒸発することで塩分が鋼材表面に濃縮された結果、腐食が発生しやすくなることを見出した。特に、台風、ハリケーン、サイクロンなどの接近や気圧配置の影響によって、最寄りの海岸の方角から強風が吹いた場合は、鋼材表面に多くの塩分が付着しやすいことが判明した。また、降雨がない限り、夜間における鋼材表面への塩分の蓄積が続き、これにより、より腐食が起きやすくなることが判明した。一方で、降雨があった場合は、雨水によって鋼材表面から塩分が洗い流されるため、腐食が進行しにくくなることが判明した。そして、このような気象条件を考慮し、観測時刻ｔＸｉ毎の風速、風向及び降雨量の気象観測値ＷＸから求まる腐食性指標Ｑ（ｔＸｉ，ＷＸ）を用いることで、将来における鋼材の腐食量を高精度で推測できる方法を見出した。
【0020】
なお、日照の有無に関しては、本発明では、観測時刻ｔＸｉが昼間か夜間かで区別することとする。日照の有無を判定するために日照時間を気象観測値ＷＸに含めてもよい。また、日照の有無を判定するために、観測地点Ｘの緯度、経度及び標高から日の出時刻及び日の入り時刻を計算で求め、観測時刻ｔＸｉが昼間または夜間のいずれに含まれるかを判定してもよい。詳細は後述する。
【0021】
また、鋼材の腐食量を正確に予測するためには、鋼材の腐食試験の結果を参照する必要がある。従来の予測方法では、腐食試験を行った地点での腐食量の予測に限定されていた。本発明の腐食量の予測方法では、腐食試験の試験地点Ｚと鋼材の利用予定地点Ｍとが異なる場所であったとしても、本発明者らが見出した腐食性指標Ｑ（ｔＸｉ，ＷＸ）を利用することで、試験地点Ｚでの腐食試験の結果に基づき、利用予定地点Ｍでの将来の腐食量を高精度に予測できる。
【0022】
以下、本発明の実施形態について説明する。
【0023】
本実施形態において腐食量の予測対象となる鋼材は特に制限はなく、普通鋼（炭素鋼）、合金鋼、特殊鋼等よりなる鋼材であってもよい。また、これらの鋼材にめっきを施しためっき鋼材でもよい。めっきの種類に特に制限はなく、一般に鋼に用いられるめっきであればよい。また、鋼材は、ステンレス鋼よりなる鋼材であってもよい。
【0024】
また、本実施形態に係る鋼材は、建材用に限定されるものではなく、屋外において使用される鋼材の腐食量の予測に好適に用いることができる。
【0025】
（鋼材の腐食量予測システム及び鋼材の腐食量予測プログラム）
まず、本実施形態の腐食量予測システム及び鋼材の腐食量予測プログラムについて説明する。
まず、本実施形態の腐食量予測システム及び鋼材の腐食量予測プログラムについて説明する。
【0026】
本実施形態の腐食量予測システムは、図１に示すように、極値推測部と、腐食量推測部とを備える。これら極値推測部及び腐食量推測部は、例えば、電子計算機の中央演算装置に備えられた機能として実現される。
【0027】
極値推測部は、観測地点Ｘが鋼材の利用予定地点Ｍである気象観測値ＷＭに基づいて、過去の観測時刻ｔＭｉ毎に鋼材の腐食性指標Ｑ（ｔＭｉ，ＷＭ）を求め、観測時刻ｔＭ１～ｔＭｎ間の腐食性指標Ｑ（ｔＭｉ，ＷＭ）の時系列データに対して極値統計解析を行うことにより、将来の評価期間ｔＦ1～ｔＦｎ内における腐食性指標の極大値Ｑｍａｘ（ｔＦｉ）を推測する機能を備える。
【0028】
また、腐食量推測部は、試験地点Ｚにおける鋼材の腐食試験から得た実測腐食量と、観測地点Ｘが試験地点Ｚである気象観測値ＷＺに基づいて腐食試験の試験期間中の観測時刻ｔＺｉ毎に得た鋼材の腐食性指標Ｑ（ｔＺｉ，ＷＺ）の最大値との関係式に、極値推測部において推測した腐食性指標の極大値Ｑｍａｘ（ｔＦｉ）を導入して、推定期間における鋼材の腐食量の推測値を得る機能を備える。腐食量推測部には、予め求めておいた上記関係式が導入されている。
なお、上記関係式は、鋼材の腐食試験から得た実測腐食量と、鋼材の腐食性指標Ｑ（ｔＺｉ，ＷＺ）の最大値の対数との関係式としてもよい。
なお、上記関係式は、鋼材の腐食試験から得た実測腐食量と、鋼材の腐食性指標Ｑ（ｔＺｉ，ＷＺ）の最大値の対数との関係式としてもよい。
【0029】
また、極値推測部及び腐食量推測部には、それぞれ、腐食性指標Ｑ（ｔＭｉ，ＷＭ）、Ｑ（ｔＺｉ，ＷＺ）を計算するための計算部が備えられていてもよい。
【0030】
計算部には、観測地点Ｘにおける観測時刻ｔＸｉ毎の風速ｕ（ｔＸｉ）、風向θｗ（ｔＸｉ）及び降雨量を少なくとも取得する第１処理部と、観測時刻ｔＸｉ毎に、一次腐食性指標Ｑ１（ｔＸｉ，ＷＸ）を計算する第２処理部と、観測時刻ｔＸｉ毎に、二次腐食性指標Ｑ２（ｔＸｉ，ＷＸ）を計算する第３処理部と、観測時刻ｔＸｉ毎に、腐食性指標Ｑ（ｔＸｉ，ＷＸ）を計算する第４処理部と、が備えられていてもよい。
【0031】
なお、極値推測部において腐食性指標Ｑ（ｔＭｉ，ＷＭ）を求める場合の計算部では、地点Ｘを鋼材の利用予定地点Ｍとし、ＷＸをＷＭとし、ｔＸｉをｔＭｉとし、ｕ（ｔＸｉ）をｕ（ｔＭｉ）とし、θｗ（ｔＸｉ）をθｗ（ｔＭｉ）とする。また、ｃ（ｔＸｉ）をｃ（ｔＭｉ）とし、ｓ（ｔＸｉ）をｓ（ｔＭｉ）とし、ｐ（ｔＸｉ）をｐ（ｔＭｉ）とする。
【0032】
また、腐食量推測部において腐食性指標Ｑ（ｔＺｉ，ＷＺ）を求める場合の計算部では、地点Ｘを鋼材の試験地点Ｚとし、ＷＸをＷＺとし、ｔＸｉをｔＺｉとし、ｕ（ｔＸｉ）をｕ（ｔＺｉ）とし、θｗ（ｔＸｉ）をθｗ（ｔＺｉ）とする。また、ｃ（ｔＸｉ）をｃ（ｔＺｉ）とし、ｓ（ｔＸｉ）をｓ（ｔＺｉ）とし、ｐ（ｔＸｉ）をｐ（ｔＺｉ）とする。
【0033】
本実施形態の腐食量予測システムは、極値推測部及び腐食量推測部により、鋼材の腐食量の推測値を生成し、これを出力する。
【0034】
次に、本実施形態の腐食量予測プログラムは、図２に示すように、極値推測ステップと、腐食量推測ステップとを備える。本実施形態の腐食量予測プログラムは、電子計算機において実行される。腐食量予測プログラムを備えた電子計算機は、極値推測ステップ及び腐食量推測ステップを順次実行することで、鋼材の腐食量の推測値を出力する。また、腐食量予測プログラムは、上記の腐食量予測システムのプログラムとして利用することもできる。
【0035】
腐食量予測プログラムの極値推測ステップは、観測地点Ｘが鋼材の利用予定地点Ｍである気象観測値ＷＭに基づいて、過去の観測時刻ｔＭｉ毎に鋼材の腐食性指標Ｑ（ｔＭｉ，ＷＭ）を求める。次いで、観測時刻ｔＭ１～ｔＭｎ間の腐食性指標Ｑ（ｔＭｉ，ＷＭ）の時系列データに対して極値統計解析を行うことにより、将来の評価期間ｔＦ1～ｔＦｎ内における腐食性指標の極大値Ｑｍａｘ（ｔＦｉ）を推測する。極値推測ステップは、図１に示す極値推測部において行ってもよい。
【0036】
また、腐食量予測プログラムの腐食量推測ステップは、試験地点Ｚにおける鋼材の腐食試験から得た実測腐食量と、観測地点Ｘが試験地点Ｚである気象観測値ＷＺに基づいて腐食試験の試験期間中の観測時刻ｔＺｉ毎に得た鋼材の腐食性指標Ｑ（ｔＺｉ，ＷＺ）の最大値との関係式に、極値推測ステップにおいて推測した腐食性指標の極大値Ｑｍａｘ（ｔＦｉ）を導入して、推定期間における鋼材の腐食量の推測値を得る。上記関係式は、予め求めておいたものを用いる。上記関係式は、鋼材の腐食試験から得た実測腐食量と、鋼材の腐食性指標Ｑ（ｔＺｉ，ＷＺ）の最大値の対数との関係式としてもよい。腐食量推測ステップは、図１に示す腐食量推測部において行ってもよい。
【0037】
また、極値推測ステップ及び腐食量推測ステップには、それぞれ、腐食性指標Ｑ（ｔＭｉ，ＷＭ）、Ｑ（ｔＺｉ，ＷＺ）を計算するための計算ステップが備えられていてもよい。
【0038】
計算ステップには、第１ステップ、第２ステップ、第３ステップおよび第４ステップが備えられていてもよい。第１ステップは、観測地点Ｘにおける観測時刻ｔＸｉ毎の風速ｕ（ｔＸｉ）、風向θｗ（ｔＸｉ）及び降雨量を少なくとも取得する。第２ステップは、観測時刻ｔＸｉ毎に、一次腐食性指標Ｑ１（ｔＸｉ，ＷＸ）を計算する。第３ステップは、観測時刻ｔＸｉ毎に、二次腐食性指標Ｑ２（ｔＸｉ，ＷＸ）を計算する。第４ステップは、観測時刻ｔＸｉ毎に、腐食性指標Ｑ（ｔＸｉ，ＷＸ）を計算する。
【0039】
なお、極値推測ステップにおいて腐食性指標Ｑ（ｔＭｉ，ＷＭ）を求める場合の計算ステップでは、地点Ｘを鋼材の利用予定地点Ｍとし、ＷＸをＷＭとし、ｔＸｉをｔＭｉとし、ｕ（ｔＸｉ）をｕ（ｔＭｉ）とし、θｗ（ｔＸｉ）をθｗ（ｔＭｉ）とする。また、ｃ（ｔＸｉ）をｃ（ｔＭｉ）とし、ｓ（ｔＸｉ）をｓ（ｔＭｉ）とし、ｐ（ｔＸｉ）をｐ（ｔＭｉ）とする。
【0040】
また、腐食量推測ステップにおいて腐食性指標Ｑ（ｔＺｉ，ＷＺ）を求める場合の計算ステップでは、地点Ｘを鋼材の試験地点Ｚとし、ＷＸをＷＺとし、ｔＸｉをｔＺｉとし、ｕ（ｔＸｉ）をｕ（ｔＺｉ）とし、θｗ（ｔＸｉ）をθｗ（ｔＺｉ）とする。また、ｃ（ｔＸｉ）をｃ（ｔＺｉ）とし、ｓ（ｔＸｉ）をｓ（ｔＺｉ）とし、ｐ（ｔＸｉ）をｐ（ｔＺｉ）とする。
【0041】
本実施形態の腐食量予測システムの極値推測部及び腐食量推測部並びに計算部の動作、並びに、腐食量予測プログラムの極値推測ステップ及び腐食量推測ステップ並びに計算ステップの動作についての詳細な説明は、鋼材の腐食量の予測方法の説明において述べる。
【0042】
なお、本実施形態の腐食量予測システム及び腐食量予測プログラムにおける入力値は、少なくとも、腐食試験の試験地点Ｚの気象観測値ＷＺ、鋼材の利用予定地点Ｍの気象観測値ＷＭ、腐食試験によって得られた実測腐食量を含む。また、出力値は、鋼材の利用予定地点Ｍにおいて鋼材を屋外に暴露させた際の将来の腐食量の予測値である。腐食量としては、腐食減量、腐食深さ、腐食の外観評点であるＲＮ値（ＲＮ：レイティングナンバー）などが利用可能である。ＲＮ値の評価方法として、ステンレス鋼ではＪＩＳ Ｇ ０５９５：２００４（ステンレス鋼の表面さび発生程度評価方法）などがある。
【0043】
（鋼材の腐食量の予測方法）
次に、本実施形態の鋼材の腐食量の予測方法について説明する。
次に、本実施形態の鋼材の腐食量の予測方法について説明する。
【0044】
本実施形態の鋼材の腐食量の予測方法は、鋼材を屋外に曝した場合の腐食量を評価する指標として、観測地点Ｘにおける観測時刻ｔＸｉ（ただしｉは１～ｎの整数）毎の風速、風向及び降雨量の気象観測値ＷＸから求まる鋼材の腐食性指標Ｑ（ｔＸｉ，ＷＸ）を用いる。
【0045】
本実施形態の鋼材の腐食量の予測方法は、腐食性指標の極大値Ｑｍａｘ（ｔＦｉ）を推測する極値推測ステップと、利用予定地点Ｍにおける、推定期間経過後の鋼材の腐食量を推測する腐食量推測ステップと、を備える。
【0046】
上記予測方法の極値推測ステップは、観測地点Ｘが鋼材の利用予定地点Ｍである気象観測値ＷＭに基づいて、過去の観測時刻ｔＭｉ毎に鋼材の腐食性指標Ｑ（ｔＭｉ，ＷＭ）を求め、観測時刻ｔＭ１～ｔＭｎ間の腐食性指標Ｑ（ｔＭｉ，ＷＭ）の時系列データに対して極値統計解析を行うことにより、将来の評価期間ｔＦ1～ｔＦｎ内における腐食性指標の極大値Ｑｍａｘ（ｔＦｉ）を推測する。
【0047】
上記予測方法の腐食量推測ステップは、鋼材の実測腐食量と、鋼材の腐食性指標Ｑ（ｔＺｉ，ＷＺ）の最大値との関係式に、極値推測ステップにおいて推測した腐食性指標の極大値Ｑｍａｘ（ｔＦｉ）を導入して、推定期間経過後の鋼材の腐食量を推測する。
【0048】
鋼材の実測腐食量と、鋼材の腐食性指標Ｑ（ｔＺｉ，ＷＺ）の最大値との関係式は、予め求めておくことが好ましい。この関係式は、試験地点Ｚにおける鋼材の腐食試験から得た実測腐食量と、観測地点Ｘが試験地点Ｚである気象観測値ＷＺに基づいて腐食試験の試験期間中の観測時刻ｔＺｉ毎に得た鋼材の腐食性指標Ｑ（ｔＺｉ，ＷＺ）の最大値との関係式である。上記関係式は、鋼材の腐食試験から得た実測腐食量と、鋼材の腐食性指標Ｑ（ｔＺｉ，ＷＺ）の最大値の対数との関係式としてもよい。
【0049】
また、極値推測ステップ及び腐食量推測ステップには、それぞれ、腐食性指標Ｑ（ｔＭｉ，ＷＭ）、Ｑ（ｔＺｉ，ＷＺ）を計算するための計算ステップが備えられていてもよい。
【0050】
まず、本実施形態の鋼材の腐食量の評価方法を説明するにあたり、腐食性指標Ｑ（ｔＸｉ，ＷＸ）について説明する。
【0051】
本実施形態の鋼材の腐食量の予測方法では、鋼材を屋外に曝した場合の腐食量を評価する指標として、観測地点Ｘにおける観測時刻ｔＸｉ（ただしｉは１～ｎの整数）毎の風速、風向及び降雨量の気象データＷＸから求まる鋼材の腐食性指標Ｑ（ｔＸｉ，ＷＸ）を用いる。
【0052】
先に述べたように、鋼材の腐食が発生しやすい気象条件として、本発明者らは、夜間から昼間に至る間で、風が観測され、かつ降雨が起きない気象条件を見出した。この気象条件では、日没から日の出までの夜間に、鋼材表面の結露水が生成するとともに風に運ばれた塩分が結露水中に溶け込み、その後の昼間の日照によって結露水が蒸発するとともに塩分が濃縮され、その結果として腐食が進展する。塩分は降雨がないまま日数が重なることで蓄積される。その一方で、降雨があると、鋼材表面から塩分が消失する。本実施形態に係る腐食性指標Ｑ（ｔＸｉ，ＷＸ）は、このような気象変動を考慮した鋼材表面における塩分量の指標となる値である。
【0053】
また、腐食性指標Ｑ（ｔＸｉ，ＷＸ）の算出に用いる腐食性蓄積係数ｓ（ｔｉ）及び腐食性発現係数ｐ（ｔＸｉ）について説明する。ｓ（ｔｉ）及びｐ（ｔＸｉ）は、鋼材表面の濡れ状態に関する値であり、０または１の値になる。値が０か１になるかは、観測時刻ｔＸｉが昼間または夜間のいずれかであるかによって決まる。
【0054】
すなわち、夜間において、結露水が付着することで鋼材表面が濡れた状態になる。一方、昼間において、夜間に付着した結露水が蒸発することで鋼材表面が濡れた状態から脱することになる。よって本実施形態では、昼間か夜間かで、鋼材表面が濡れているか否かを間接的に判別し、その判別結果に基づき、ｓ（ｔｉ）及びｐ（ｔＸｉ）の値を定めるものとする。
【0055】
なお、本実施形態では、日没時刻から日の出時刻までの期間を夜間とし、日の出時刻から日没時刻までの期間を昼間とする。
【0056】
昼間か夜間かの判別を気象観測値ＷＸから求める場合は、気象観測値ＷＸとして、観測時刻ｔＸｉ－１～ｔＸｉの間の日照時間を取得する。観測時刻ｔＸｉが夜間である場合は日照時間の観測値が存在しないため、観測値の有無により、観測時刻ｔＸｉが昼間か夜間かを区別することができる。これにより、気象観測値ＷＸから腐食性蓄積係数ｓ（ｔｉ）と腐食性発現係数ｐ（ｔＸｉ）を決定することが出来る。
【0057】
また、日の出時刻と日没時刻は、日付（年、月、日）と地点Ｘの経度、緯度及び標高から計算できるので、計算により求めた日の出時刻と日没時刻から昼間か夜間かを判定し、これにより、腐食性蓄積係数ｓ（ｔｉ）と腐食性発現係数ｐ（ｔＸｉ）を決定してもよい。
【0058】
以下、腐食性指標Ｑ（ｔＸｉ，ＷＸ）の算出法を説明する。本実施形態では、以下の第１ステップ～第４ステップからなる計算ステップにより、一次腐食性指標Ｑ１（ｔＸｉ，ＷＸ）及び二次腐食性指標Ｑ２（ｔＸｉ，ＷＸ）を順次求め、更に、二次腐食性指標Ｑ２（ｔＸｉ，ＷＸ）から本実施形態に係る腐食性指標Ｑ（ｔＸｉ，ＷＸ）を求める。
【0059】
なお、本実施形態の説明において使用する各種のパラメータ記号における添え字Ｘは、観測地点Ｘを意味するものであり、利用可能地点であるＭまたは試験地点であるＺに置換されることが予定されている。また、観測時刻ｔＸｉは、観測測時刻ｔＸ１、ｔＸ２、ｔＸ3、…ｔＸｎを一般化した表記であり、任意の観測時刻を観測時刻ｔＸｉと表記する場合がある。また、観測時刻ｔＸｉ～ｔＸｎを要素とする集合を、観測期間ｔＸｉ～ｔＸｎと表記とする場合がある。
【0060】
まず、第１ステップでは、観測地点Ｘにおける観測時刻ｔＸｉ毎の風速ｕ（ｔＸｉ）、風向θｗ（ｔＸｉ）及び降雨量を少なくとも取得する。また、観測時刻ｔＸｉ毎に日照時間を取得してもよい。観測地点Ｘにおける風速ｕ（ｔＸｉ）、風向θｗ（ｔＸｉ）、降雨量及び日照時間を総称して、気象観測値ＷＸという場合がある。
【0061】
第１ステップにおける観測地点Ｘ及び気象観測値ＷＸは、予測方法の段階に応じて適宜選択すればよい。すなわち、極値推測ステップでは、観測地点Ｘを鋼材の利用予定地点Ｍとし、気象観測値ＷＸを利用予定地点Ｍにおける気象観測値ＷＭとする。また、腐食量推測ステップでは、観測地点Ｘを試験地点Ｚとし、気象観測値ＷＸを試験地点Ｚにおける気象観測値ＷＺとする。
【0062】
気象観測値ＷＸは、例えば、公的な気象観測機関が管理する気象観測所において測定され、公表されているデータを用いればよい。例えば、日本国の場合は、気象庁が管理する気象観測所において測定され、公表されているデータを用いてもよい。すなわち、気象観測値ＷＸは、観測地点Ｘに近い気象観測所で測定されたデータとすればよい。また、気象観測値ＷＸは、公的な気象観測機関が管理する気象観測所で測定されたものに限らず、鋼材の利用予定地点Ｍや試験地点Ｚの近傍に設置された臨時の気象観測施設において測定された測定データを用いてもよい。
【0063】
観測地点Ｘにおける気象の観測期間ｔＸｉ～ｔＸｎの長さは、特に制限されないが、以下に説明する期間長さであると好ましい。
【0064】
利用予定地点Ｍにおける過去の観測期間ｔＭｉ～ｔＭｎの長さは、例えば、１年以上、５０年以下のうちの何れかの期間とすればよく、好ましくは５年以上、５０年以下、より好ましくは１０年以上、５０年以下とすればよい。観測期間ｔＭｉ～ｔＭｎが長いほど、極値統計解析時に利用される時系列データに多数のデータが含まれることになり、極値統計解析の信頼性が向上する。そのため、観測期間ｔＭｉ～ｔＭｎは、なるべく長期間であるとよい。本実施形態の予測方法では、塩分の飛散状況から将来の腐食量を極値統計解析により予測するが、塩分の飛散量は、年単位で変動し、例えば、海塩が多く飛散する年と少なく飛散する年とがある。そして、観測期間ｔＭｉ～ｔＭｎが短く、かつ、その期間内に海塩が多く飛散する年が含まれると、将来の予測腐食量が高めになる可能性がある。このような予測値の偏りを防止するためには、観測期間ｔＭｉ～ｔＭｎは長いほうがよい。
【0065】
また、利用地点Ｍにおける将来の評価期間ｔＦ1～ｔＦｎ長さは、例えば、１ヶ月以上、１００年以下のうちの何れかの期間とすればよく、好ましくは２年以上、１００年以下のうちの何れかの期間とすればよい。
【0066】
更に、腐食量推測ステップでは、試験地点Ｚにおける観測期間ｔＺｉ～ｔＺｎの長さとして例えば、１日以上、３０年以下のうちの何れかの期間とすればよく、好ましくは１ヶ月以上、１０年以下のうちの何れかの期間とすればよく、１ヶ月以上５年以下でもよい。
各観測時刻ｔＸ１、…ｔＸｎの間隔は例えば、１０分毎、１５分毎または３０分毎の観測時刻とすればよい。試験地点Ｚにおける観測期間ｔＺｉ～ｔＺｎは、腐食量推測ステップにおいて腐食性指標Ｑ（ｔＺｉ，ＷＺ）の最大値と実測腐食量との関係式を求める際に、Ｑ（ｔＺｉ，ＷＺ）の最大値の範囲が一定の範囲になればよい。Ｑ（ｔＺｉ，ＷＺ）の最大値は、対数に変換したものでもよい。
各観測時刻ｔＸ１、…ｔＸｎの間隔は例えば、１０分毎、１５分毎または３０分毎の観測時刻とすればよい。試験地点Ｚにおける観測期間ｔＺｉ～ｔＺｎは、腐食量推測ステップにおいて腐食性指標Ｑ（ｔＺｉ，ＷＺ）の最大値と実測腐食量との関係式を求める際に、Ｑ（ｔＺｉ，ＷＺ）の最大値の範囲が一定の範囲になればよい。Ｑ（ｔＺｉ，ＷＺ）の最大値は、対数に変換したものでもよい。
【0067】
風速ｕ（ｔＸｉ）は、観測時刻ｔＸｉにおける瞬間風速（ｍ／ｓ）でもよく、観測時刻ｔＸｉ－１～ｔＸｉの間における平均風速（ｍ／ｓ）でもよい。風向θｗ（ｔＸｉ）は、北の方角を０°とした場合の観測時刻ｔＸｉにおける風向の方位とする。観測時刻ｔＸｉにおける降雨量及び日照時間はそれぞれ、観測時刻ｔＸｉ－１から観測時刻ｔＸｉまでの間における降雨量の積算降雨量及び日照時間の積算時間とする。例えば、観測時刻ｔＸｉが１０分毎である場合は、観測時刻ｔＸｉ－１～ｔＸｉの１０分間の積算降雨量及び積算日照時間とする。
【0068】
次に第２ステップでは、観測時刻ｔＸｉ毎に、一次腐食性指標Ｑ１（ｔＸｉ，ＷＸ）を下記式（１）により計算する。
【0069】
Ｑ１（ｔＸｉ，ＷＸ）＝（ｄ＋１）－０．６｛ｕ（ｔＸｉ）・ｃｏｓ｜θｓ－θｗ（ｔＸｉ）｜｝２ …（１）
【0070】
式（１）におけるｕ（ｔＸｉ）は、観測時刻ｔＸｉの風速（ｍ／ｓ）であり、θｗ（ｔＸｉ）は観測時刻ｔＸｉの北を０°とした場合の風向（°）であり、θｓは、観測地点Ｘに最も近い海岸からの観測地点Ｘの方位であって北を０°とした場合の方位（°）である。
【0071】
また、式（１）におけるｄは、観測地点Ｘに最も近い海岸と観測地点Ｘとの距離（ｍ）である。
【0072】
一次腐食性指標Ｑ１（ｔＸｉ，ＷＸ）は、鋼材への塩分付着量の指標となるものであり、風速ｕ（ｔＸｉ）及びθｓ、θｗ（ｔＸｉ）に基づいて算出される。式（１）における「ｃｏｓ｜θｓ－θｗ（ｔＸｉ）｜」の項は、観測地点Ｘにおける風速ｕ（ｔＸｉ）のうち、海岸方向から吹く風の風速成分を求めるための係数である。また、鋼材への塩分の付着量は、風速ｕ（ｔＸｉ）の２乗に比例する経験則がある。更に、（ｄ＋１）－０．６は、海岸からの距離ｄによって飛来してくる海塩が減衰することを考慮した減衰係数である。よって、鋼材への塩分付着量の指標である一次腐食性指標Ｑ１（ｔＸｉ，ＷＸ）を上記式（１）の通りとする。
【0073】
次に第３ステップでは、観測時刻ｔＸｉ毎に、二次腐食性指標Ｑ２（ｔＸｉ，ＷＸ）を下記式（２）により計算する。
【0074】
Ｑ２（ｔＸｉ，ＷＸ）＝｛Ｑ２（ｔＸｉ－１，ＷＸ）＋Ｑ１（ｔＸｉ，ＷＸ）・ｓ（ｔＸｉ）・（ｔＸｉ－ｔＸｉ－１）｝・ｃ（ｔＸｉ） …（２）
【0075】
ただし、式（２）におけるｃ（ｔＸｉ）は腐食性消失係数であって、観測時刻ｔＸｉ－１～ｔＸｉ間の降雨量が閾値超の場合にｃ（ｔＸｉ）＝０、観測時刻ｔＸｉ－１～ｔＸｉ間の降雨量が閾値以下の場合にｃ（ｔＸｉ）＝１とする。閾値は予め試験を行い、鋼材表面に付着させた塩分が消失する降雨量を実験的に定めることができる。閾値は、例えば０～５ｍｍ／ｈｏｕｒの範囲で決定すればよい。
【0076】
また、式（２）におけるｓ（ｔＸｉ）は腐食性蓄積係数であって、観測時刻ｔＸｉ－１～ｔＸｉ間が昼間（日の出～日没までの期間）の場合にｓ（ｔＸｉ）＝０、観測時刻ｔＸｉ－１～ｔＸｉ間が夜間（日没～日の出までの期間）の場合にｓ（ｔＸｉ）＝１とする。
腐食性蓄積係数ｓ（ｔＸｉ）の決定方法は、先に述べたとおりであり、気象データＷＸに基づき決定してもよく、地点Ｘの緯度、経度、標高に基づき計算された日の出時刻及び日の入り時刻に基づき決定してもよい。
腐食性蓄積係数ｓ（ｔＸｉ）の決定方法は、先に述べたとおりであり、気象データＷＸに基づき決定してもよく、地点Ｘの緯度、経度、標高に基づき計算された日の出時刻及び日の入り時刻に基づき決定してもよい。
【0077】
式（２）では、観測時刻ｔＸｉの直前の観測時刻である観測時刻ｔＸｉ－１における二次腐食性指標Ｑ２（ｔＸｉ－１，ＷＸ）に、観測時刻ｔｉにおける一次腐食性指標Ｑ１（ｔＸｉ，ＷＸ）が加算される。また、一次腐食性指標Ｑ１（ｔＸｉ，ＷＸ）には、腐食性蓄積係数ｓ（ｔｉ）及び（ｔＸｉ－ｔＸｉ－１）が乗算される。従って、式（２）の一部をなすＱ１（ｔＸｉ，ＷＸ）・ｓ（ｔＸｉ）・（ｔＸｉ－ｔＸｉ－１）の多項式は、観測時刻ｔＸｉ－１～ｔＸｉ間が昼間（日の出～日没までの期間）の場合には０となり、観測時刻ｔＸｉ－１～ｔＸｉ間が夜間（日没～日の出までの期間）の場合には０超の値になる。
【0078】
また、式（２）では、｛Ｑ２（ｔＸｉ－１，ＷＸ）＋Ｑ１（ｔＸｉ，ＷＸ）・ｓ（ｔＸｉ）・（ｔＸｉ－ｔＸｉ－１）｝の多項式に、ｃ（ｔＸｉ）が乗算される。従って、Ｑ２（ｔＸｉ，ＷＸ）は、観測時刻ｔＸｉ－１～ｔＸｉ間の降雨量が閾値超の場合に０となり、観測時刻ｔＸｉ－１～ｔＸｉ間の降雨量が閾値以下の場合には０超の値になる。
【0079】
よって、二次腐食性指標Ｑ２（ｔＸｉ，ＷＸ）は、日没から日の出までの夜間に、鋼材表面の結露水が生成するとともに、風に運ばれた塩分が結露水中に溶け込むことによって蓄積された塩分量の指標となる。
【0080】
次に第４ステップでは、観測時刻ｔＸｉ毎に、腐食性指標Ｑ（ｔＸｉ，ＷＸ）を下記式（３Ａ）により計算する。
【0081】
Ｑ（ｔＸｉ，ＷＸ）＝Ｑ２（ｔＸｉ，ＷＸ）・ｐ（ｔＸｉ） …（３Ａ）
【0082】
式（３Ａ）におけるｐ（ｔＸｉ）は腐食性発現係数であって、観測時刻ｔＸｉ－１～ｔＸｉ間が夜間（日没～日の出までの期間）の場合にｐ（ｔＸｉ）＝０、観測時刻ｔＸｉ－１～ｔＸｉ間が昼間（日の出～日没までの期間）の場合にｐ（ｔＸｉ）＝１、とする。腐食性発現係数ｐ（ｔＸｉ）の決定方法は、先に述べたとおりであり、気象データＷＸに基づき決定してもよく、地点Ｘの緯度、経度、標高に基づき計算された日の出時刻及び日の入り時刻に基づき決定してもよい。
【0083】
式（３Ａ）では、Ｑ２（ｔＸｉ，ＷＸ）にｐ（ｔＸｉ）が乗算される。従って、Ｑ（ｔＸｉ，ＷＸ）は、観測時刻ｔＸｉ－１～ｔＸｉ間が夜間（日没～日の出までの期間）の場合０となり、観測時刻ｔＸｉ－１～ｔＸｉ間が昼間（日の出～日没までの期間）の場合に０超の値になる。
【0084】
よって、腐食性指標Ｑ（ｔＸｉ，ＷＸ）は、日没から日の出までの夜間に、鋼材表面の結露水が生成するとともに風に運ばれた塩分が結露水中に溶け込み、その後の昼間の日照によって結露水が蒸発するとともに塩分が濃縮された場合の塩分量の指標となる。
一般に、塩分の付着量が大きいほど、鋼材が腐食しやすくなるので、腐食性指標Ｑ（ｔＸｉ，ＷＸ）が大きいほど、鋼材の腐食量が増大する傾向にあると言える。
一般に、塩分の付着量が大きいほど、鋼材が腐食しやすくなるので、腐食性指標Ｑ（ｔＸｉ，ＷＸ）が大きいほど、鋼材の腐食量が増大する傾向にあると言える。
【0085】
上記の計算ステップを実行することで、時刻ｔＸｉ毎に腐食性指標Ｑ（ｔＸｉ，ＷＸ）が得られる。
極値推測ステップの計算ステップでは、地点Ｘを鋼材の利用予定地点Ｍとし、ＷＸをＷＭとし、ｔＸｉをｔＭｉとし、ｕ（ｔＸｉ）をｕ（ｔＭｉ）とし、θｗ（ｔＸｉ）をθｗ（ｔＭｉ）とし、ｃ（ｔＸｉ）をｃ（ｔＭｉ）とし、ｓ（ｔＸｉ）をｓ（ｔＭｉ）とし、ｐ（ｔＸｉ）をｐ（ｔＭｉ）とすることで、腐食性指標Ｑ（ｔＭｉ，ＷＭ）が得られる。
腐食量推測ステップの計算ステップでは、地点Ｘを鋼材の試験地点Ｚとし、ＷＸをＷＺとし、ｔＸｉをｔＺｉとし、ｕ（ｔＸｉ）をｕ（ｔＺｉ）とし、θｗ（ｔＸｉ）をθｗ（ｔＺｉ）とし、ｃ（ｔＸｉ）をｃ（ｔＺｉ）とし、ｓ（ｔＸｉ）をｓ（ｔＺｉ）とし、ｐ（ｔＸｉ）をｐ（ｔＺｉ）とすることで、腐食性指標Ｑ（ｔＺｉ，ＷＺ）が得られる。
極値推測ステップの計算ステップでは、地点Ｘを鋼材の利用予定地点Ｍとし、ＷＸをＷＭとし、ｔＸｉをｔＭｉとし、ｕ（ｔＸｉ）をｕ（ｔＭｉ）とし、θｗ（ｔＸｉ）をθｗ（ｔＭｉ）とし、ｃ（ｔＸｉ）をｃ（ｔＭｉ）とし、ｓ（ｔＸｉ）をｓ（ｔＭｉ）とし、ｐ（ｔＸｉ）をｐ（ｔＭｉ）とすることで、腐食性指標Ｑ（ｔＭｉ，ＷＭ）が得られる。
腐食量推測ステップの計算ステップでは、地点Ｘを鋼材の試験地点Ｚとし、ＷＸをＷＺとし、ｔＸｉをｔＺｉとし、ｕ（ｔＸｉ）をｕ（ｔＺｉ）とし、θｗ（ｔＸｉ）をθｗ（ｔＺｉ）とし、ｃ（ｔＸｉ）をｃ（ｔＺｉ）とし、ｓ（ｔＸｉ）をｓ（ｔＺｉ）とし、ｐ（ｔＸｉ）をｐ（ｔＺｉ）とすることで、腐食性指標Ｑ（ｔＺｉ，ＷＺ）が得られる。
【0086】
なお、上記の計算ステップは、上述の腐食量予測システムの極値推測部または腐食量推測部において実行されてもよい。また、第１ステップ～第４ステップは、それぞれ、極値推測部または腐食量推測部に備えられた第１計算部～第４計算部において実行されてもよい。
【0087】
次に、本実施形態の予測方法の極値推測ステップ及び腐食量推測ステップについて説明する。なお、極値推測ステップ及び腐食量推測ステップは、上述の腐食量予測システムの極値推測部及び腐食量推測部において実行されてもよい。
【0088】
（極値推測ステップ）
極値推測ステップでは、まず、鋼材の利用予定地点Ｍの気象観測値ＷＭに基づいて、過去の観測時刻ｔＭｉ毎に鋼材の腐食性指標Ｑ（ｔＭｉ，ＷＭ）を求める。腐食性指標Ｑ（ｔＭｉ，ＷＭ）の算出は、上述した計算ステップによって行う。
極値推測ステップでは、まず、鋼材の利用予定地点Ｍの気象観測値ＷＭに基づいて、過去の観測時刻ｔＭｉ毎に鋼材の腐食性指標Ｑ（ｔＭｉ，ＷＭ）を求める。腐食性指標Ｑ（ｔＭｉ，ＷＭ）の算出は、上述した計算ステップによって行う。
【0089】
これにより、鋼材の利用予定地点Ｍにおける過去の期間ｔＭ１～ｔＭｎに含まれる観測時刻ｔＭ１、ｔＭ２、…ｔＭｎ毎に、腐食性指標Ｑ（ｔＭ１，ＷＭ）、Ｑ（ｔＭ２，ＷＭ）、…Ｑ（ｔＭｎ，ＷＭ）が求まる。このようにして、観測時刻ｔＭ１～ｔＭｎ間の腐食性指標Ｑ（ｔＭｉ，ＷＭ）の時系列データが得られる。
【0090】
図３に、観測時刻ｔＭ１～ｔＭｎ間の腐食性指標Ｑ（ｔＭｉ，ＷＭ）の時系列データをグラフ化した一例を示す。図３は、観測時刻ｔＭｉを横軸とし、腐食性指標Ｑ（ｔＭｉ，ＷＭ）を縦軸とするグラフである。なお、図３は、日本の瀬戸内地方の海岸から約７００ｍ離れた地点を利用予定地点Ｍとした場合の腐食性指標Ｑ（ｔＭｉ，ＷＭ）の時系列データの一例である。
【0091】
次に、極値推測ステップでは、先に得られた観測時刻ｔＭ１～ｔＭｎ間の腐食性指標Ｑ（ｔＭｉ，ＷＭ）の時系列データに対して極値統計解析を行うことにより、将来の評価期間ｔＦ1～ｔＦｎ内における腐食性指標の極大値Ｑｍａｘ（ｔＦｉ）を推測する。
【0092】
極値統計解析は、与えられた時間的、空間的なデータから、どの様な大きな値がどれくらいの確率で出現するのかを推測するための統計手法である。すなわち、観測時刻ｔＭ１～ｔＭｎ間の腐食性指標Ｑ（ｔＭｉ，ＷＭ）をデータにして極値統計解析を行うことによって、過去の気象観測値から、将来、どのような大きさの腐食性指標Ｑ（ｔＭｉ，ＷＭ）が出現するかを推測できる。具体的には、腐食性指標Ｑ（ｔＭｉ，ＷＭ）を大きい順に並べたデータとし、このデータ分布がＧｕｍｂｅｌ分布に従うとしてＧｕｍｂｅｌ分布確率紙を作成する。そして、Ｇｕｍｂｅｌ分布確率紙を利用して、将来の評価期間ｔＦ1～ｔＦｎ内における腐食性指標の極大値Ｑｍａｘ（ｔＦｉ）を推測する。
【0093】
図４Ａ及び図４Ｂには、極値統計解析の解析結果の一例を示す。極値統計解析を行うことで、図４Ａに示す観測時刻ｔＭ１～ｔＭｎ間の腐食性指標Ｑ（ｔＭｉ，ＷＭ）の時系列データから、図４Ｂに示す将来の評価期間ｔＦ1～ｔＦｎ内における腐食性指標の極大値のグラフが得られる。図４Ｂのグラフを参照することで、例えば１０年後における利用地点Ｍにおける腐食性指標の極大値Ｑｍａｘ（ｔＦｉ）を知ることができる。なお、図４Ａ及び図４Ｂは、日本の沖縄県の海岸から約７５ｍ離れた地点を利用予定地点Ｍとした場合の解析結果の一例である。
【0094】
（腐食量推測ステップ）
腐食量推測ステップでは、予め求めておいた、腐食試験結果の実測腐食量と腐食性指標Ｑ（ｔＺｉ，ＷＺ）の最大値との関係式を用いる。以下、この関係式の導出手順について説明する。
腐食量推測ステップでは、予め求めておいた、腐食試験結果の実測腐食量と腐食性指標Ｑ（ｔＺｉ，ＷＺ）の最大値との関係式を用いる。以下、この関係式の導出手順について説明する。
【0095】
関係式を求めるためには、試験地点Ｚにおける鋼材の腐食試験から得た実測腐食量と、観測地点Ｘが試験地点Ｚである気象観測値ＷＺに基づいて腐食試験の試験期間中の観測時刻ｔＺｉ毎に得た鋼材の腐食性指標Ｑ（ｔＺｉ，ＷＺ）とを求める。
【0096】
鋼材の実測腐食量は、試験地点Ｚにおいて実際に鋼材を屋外に配置し、その腐食量を観察することにより行う。腐食量の評価は、例えば、腐食減量や孔食深さでもよく、ＲＮによる評価でもよい。
【0097】
鋼材の腐食量は鋼種によって異なるため、利用予定地点Ｍにおいて利用予定の鋼材と同じ鋼種について腐食試験を行うとよい。
【0098】
また、試験地点Ｚは、鋼材の利用予定地点Ｍと同じ地点であってもよく、異なる地点であってもよい。更に、試験地点Ｚは１箇所でもよく、複数箇所でもよい。予測精度を高めるために好ましくは、試験地点Ｚは複数箇所がよい。すなわち、試験地点Ｚは一つの地点Ｚαでもよいが、好ましくは複数の異なる地点Ｚα、Ｚβ、Ｚγ…がよい。試験地点Ｚを複数箇所とする場合は、それら複数の試験地点には、海からの距離が異なる地点が含まれることがより好ましい。また、複数の試験地点には、気象条件の異なる地点が含まれることがより好ましい。複数の試験地点をこのようにすることで、腐食試験で得られる腐食性指標Ｑ（ｔＺｉ，ＷＺ）の最大値ＱＲｍａｘの範囲が広くなり、予測の精度を向上できる。また、各試験地点Ｚα、Ｚβ、Ｚγ…における試験期間はそれぞれ、同じ長さの期間であってもよく、異なる長さの期間であってもよい。更に、各試験地点Ｚα、Ｚβ、Ｚγ…における試験期間に含まれる観測時刻ｔＺ１、ｔＺ２、…ｔＺｎ、具体的にはｔＺα１～ｔＺαｎ、ｔＺβ１～ｔＺβｎ、ｔＺγ１～ｔＺγｎ等はそれぞれ、同じ時刻であってもよく、異なる時刻であってもよい。
【0099】
また、腐食性指標Ｑ（ｔＺｉ，ＷＺ）については、試験地点Ｚにおける腐食試験の試験期間ｔＺ１～ｔＺｎの気象観測値ＷＺに基づき、観測時刻ｔＺｉ毎に鋼材の腐食性指標Ｑ（ｔＺｉ，ＷＺ）を求める。腐食性指標Ｑ（ｔＺｉ，ＷＺ）の算出は、上述した計算ステップによって行う。これにより、試験地点Ｚにおける試験期間ｔＺ１～ｔＺｎに含まれる観測時刻ｔＺ１、ｔＺ２、…ｔＺｎ毎に、腐食性指標Ｑ（ｔＺ１，ＷＺ）、Ｑ（ｔＺ２，ＷＺ）、…Ｑ（ｔＺｎ，ＷＺ）が求まる。
試験地点Ｚが複数の場合は、試験地点毎に腐食性指標Ｑ（ｔＺ１，ＷＺ）を求める。この場合の気象観測値ＷＺは、試験地点Ｚα、Ｚβ、Ｚγ…におけるそれぞれの試験期間に含まれる観測時刻ｔＺ１、ｔＺ２、…ｔＺｎ毎、具体的にはｔＺα１～ｔＺαｎ、ｔＺβ１～ｔＺβｎ、ｔＺγ１～ｔＺγｎ等の気象観測値を用いる。
試験地点Ｚが複数の場合は、試験地点毎に腐食性指標Ｑ（ｔＺ１，ＷＺ）を求める。この場合の気象観測値ＷＺは、試験地点Ｚα、Ｚβ、Ｚγ…におけるそれぞれの試験期間に含まれる観測時刻ｔＺ１、ｔＺ２、…ｔＺｎ毎、具体的にはｔＺα１～ｔＺαｎ、ｔＺβ１～ｔＺβｎ、ｔＺγ１～ｔＺγｎ等の気象観測値を用いる。
【0100】
次に、腐食試験毎に、得られた腐食性指標Ｑ（ｔＺ１，ＷＺ）、Ｑ（ｔＺ２，ＷＺ）、…Ｑ（ｔＺｎ，ＷＺ）の中から最大値ＱＲｍａｘを抽出し、抽出した最大値ＱＲｍａｘを各腐食試験毎の試験結果である実測腐食量に紐付ける。そして、紐付けられた最大値ＱＲｍａｘの対数と実測腐食量との関係式を求める。図５には、ＳＵＳ３０４についての最大値ＱＲｍａｘの対数と実測腐食量（ＲＮ）とのプロット及び関係式のグラフの一例を示す。図５における点線が関係式の線である。
【0101】
関係式は、最大値ＱＲｍａｘの対数と実測腐食量（ＲＮ）とのプロットに対して最小二乗法などによって近似された近似線の関数式とすればよい。例えば、関係式として、最大値ＱＲｍａｘの対数を独立変数Ｘとし、実測腐食量（ＲＮ）を従属変数ＹとするＹ＝αＸ＋βの一次関数式としてもよい。α、βはそれぞれ定数である。この関係式は、試験地点Ｚの地点数を複数にすることで、その精度の向上を図ることができる。なお、最大値ＱＲｍａｘを対数に変換しない場合は、Ｙ＝αＸ＋βの一次関数式ではなく、近似が容易な別の関数式を用いて、回帰計算等により係数を求めてもよい。
【0102】
図５に示すように、ＳＵＳ３０４は、最大値ＱＲｍａｘが増加するに伴いＲＮが低下し、鋼材の外観が悪化する傾向になる。なお、図５は、日本の瀬戸内地方および沖縄県の海岸からの距離１０～１０００ｍ離れた複数の地点を試験地点Ｚとした場合の、最大値ＱＲｍａｘの対数と実測腐食量との関係を示すグラフである。
【0103】
そして、腐食量推測ステップでは、予め求めておいた最大値ＱＲｍａｘの対数と実測腐食量との関係式に、極値推測ステップにおいて推測した腐食性指標の極大値Ｑｍａｘ（ｔＦｉ）を導入する。これにより、利用予定地点Ｍでの将来の推定期間における鋼材の腐食量の推定値を得ることができる。
【0104】
図６には、本実施形態の推測方法のアウトプットの一例を示す。例えば、各種の鋼（Ｓｔｅｅｌ Ａ～Ｓｔｅｅｌ Ｅ）について、１年先から数十年先のＲＮ値の推測値を知ることが可能になる。
【0105】
以上説明したように、本実施形態の鋼材の腐食量の予測方法、腐食量予測システム及び腐食量予測プログラムでは、鋼材の利用予定地点Ｍにおける過去の観測時刻ｔＭｉ毎の腐食性指標ＱＭ（ｔＭｉ，ＷＭ）の時系列データに対して極値統計解析を行うことで、将来の評価期間ｔＦ1～ｔＦｎ内における腐食性指標の極大値Ｑｍａｘ（ｔＦｉ）を推測する。次いで、腐食試験の結果から予め求めておいた鋼材の実測腐食量と鋼材の腐食性指標Ｑ（ｔＺｉ，ＷＺ）との関係式に、腐食性指標の極大値Ｑｍａｘ（ｔＦｉ）を導入する。本発明は、このようにして、将来の鋼材の腐食量を推測するので、従来の年平均値に基づいて腐食量を予測する方法に比べて、腐食量の予測精度を高めることができる。
また、本実施形態の鋼材の腐食量の予測方法、腐食量予測システム及び腐食量予測プログラムは、試験地点Ｚにおける鋼材の腐食試験の結果に基づき、利用予定地点Ｍにおける鋼材の将来の腐食量を推測するものである。本実施形態では、試験地点Ｚにおける鋼材の実測腐食量と腐食性指標Ｑ（ｔＺｉ，ＷＺ）との対応関係と、利用予定地点Ｍの気象観測値ＷＭから得た腐食性指標Ｑ（ｔＭｉ，ＷＭ）の時系列データとが得られていればよい。よって、本発明は汎用性に優れる。
また、本実施形態の鋼材の腐食量の予測方法、腐食量予測システム及び腐食量予測プログラムは、試験地点Ｚにおける鋼材の腐食試験の結果に基づき、利用予定地点Ｍにおける鋼材の将来の腐食量を推測するものである。本実施形態では、試験地点Ｚにおける鋼材の実測腐食量と腐食性指標Ｑ（ｔＺｉ，ＷＺ）との対応関係と、利用予定地点Ｍの気象観測値ＷＭから得た腐食性指標Ｑ（ｔＭｉ，ＷＭ）の時系列データとが得られていればよい。よって、本発明は汎用性に優れる。
【0106】
また、本実施形態の鋼材の腐食量の予測方法、腐食量予測システム及び腐食量予測プログラムによれば、一次腐食性指標Ｑ１（ｔＸｉ，ＷＸ）及び二次腐食性指標Ｑ２（ｔＸｉ，ＷＸ）を順次求め、更に、二次腐食性指標Ｑ２（ｔＸｉ，ＷＸ）から腐食性指標Ｑ（ｔＸｉ，ＷＸ）を求める。ここで、一次腐食性指標Ｑ１（ｔＸｉ，ＷＸ）は、鋼材への塩分付着量の指標となるものであり、風速ｕ（ｔＸｉ）及び風向θｓ、θｗ（ｔＸｉ）に基づいて算出される。また、二次腐食性指標Ｑ２（ｔＸｉ，ＷＸ）は、観測時刻ｔＸｉの直前の観測時刻である観測時刻ｔＸｉ－１における二次腐食性指標Ｑ２（ｔＸｉ－１，ＷＸ）に、観測時刻ｔＸｉにおける一次腐食性指標Ｑ１（ｔＸｉ，ＷＸ）を加算し、更に腐食の主原因である塩分が降雨により洗い流されることを腐食性消失係数ｃ（ｔＸｉ）として導入するとともに、鋼材表面に結露水が発生しやすい夜間に塩分の付着が増大することを腐食性蓄積係数ｓ（ｔＸｉ）として導入することにより算出される。そして、腐食性指標Ｑ（ｔＸｉ，ＷＸ）は、夜間に付着した塩分を含む結露水が昼間の日照により乾燥されることで腐食が進展する事象を腐食性発現係数ｐ（ｔＸｉ）で数値化して、腐食性発現係数ｐ（ｔＸｉ）を二次腐食性指標Ｑ２（ｔＸｉ，ＷＸ）に乗じることで算出される。これにより、鋼材の塩分付着量に影響を与える気象状況が反映された腐食性指標Ｑ（ｔＸｉ，ＷＸ）を得ることができ、これを用いることで、腐食量の予測精度を高めることができる。
【0107】
次に、本実施形態の鋼材の腐食量の予測方法、腐食量予測システム及び腐食量予測プログラムにおける腐食性指標Ｑ（ｔＸｉ，ＷＸ）の変形例について説明する。第１の変形例は、腐食性指標Ｑ（ｔＸｉ，ＷＸ）の温度補正に関する変形例であり、第２の変形例は、鋼材表面の濡れ状態の判定結果に基づいて腐食性指標Ｑ（ｔＸｉ，ＷＸ）を補正する変形例である。これらの変形例を採用することで、腐食量の予測精度をより一層高めることができるようになる。
【0108】
（第１の変形例）
第１の変形例では、腐食量の予測方法または腐食量予測プログラムの第４ステップにおいて、また、腐食量予測システムの第４処理部において、上記の式（３Ａ）に代えて、下記式（３Ｂ）によって、観測時刻ｔＸｉ毎に、腐食性指標Ｑ（ｔＸｉ，ＷＸ）を計算する。
第１の変形例では、腐食量の予測方法または腐食量予測プログラムの第４ステップにおいて、また、腐食量予測システムの第４処理部において、上記の式（３Ａ）に代えて、下記式（３Ｂ）によって、観測時刻ｔＸｉ毎に、腐食性指標Ｑ（ｔＸｉ，ＷＸ）を計算する。
【0109】
Ｑ（ｔＸｉ，ＷＸ）＝１００．０４Ｔｂ・Ｑ２（ｔＸｉ，ＷＸ）・ｐ（ｔＸｉ） …（３Ｂ）
【0110】
ただし、式（３Ｂ）におけるＴｂは、日の出時の気温Ｔｅｍｐ（℃）とした場合にＴｅｍｐ（℃）～（Ｔｅｍｐ＋２）℃の範囲の何れかの温度である。好ましくは、Ｔｂは日の出時の気温Ｔｅｍｐ（℃）としてもよい。
【0111】
鋼材の腐食は、鋼材表面に付着する塩分量のほか、温度の影響も受ける。このため、観測地点Ｘである鋼材の利用予定地点Ｍおよび試験地点Ｚの気温差が大きくなると、腐食量の予測精度が低下する場合がある。そこで、腐食量の予測精度を高めるために、観測地点Ｘの気温の補正項をＱ（ｔＸｉ，ＷＸ）の計算式に加入することが望ましい。
【0112】
本発明者らが検討したところ、補正項を１０の累乗の形とし、また、累乗の指数は気温Ｔｂ（℃）と係数との積とすることが好ましいことがわかった。そして、鋼材の実測腐食量および鋼材の腐食性指標Ｑ（ｔＺｉ，ＷＺ）から関係式を求める際に、補正項における指数を回帰計算によって求めたところ、０．０４の値を得た。よって、温度の補正項における指数の係数は、０．０４とする。
【0113】
温度の補正項の指数に含まれる気温Ｔｂ（℃）は、日の出時の気温をＴｅｍｐとした場合にＴｅｍｐ（℃）～（Ｔｅｍｐ＋２）℃の範囲の何れかの温度とする。鋼材は、その表面が濡れた状態にあるときに表面の水分に塩分が溶け込み、その後、水分が蒸発して塩分濃度が高まった際に腐食が発生しやすく、また、水分が蒸発したときの気温が腐食量に影響する。このような水分の蒸発は、日の出前後の気温上昇時に起きる。そこで、補正項に加入する気温Ｔｂを、日の出時の気温をＴｅｍｐ（℃）とした場合にＴｅｍｐ（℃）～（Ｔｅｍｐ＋２）℃の範囲の何れかの温度とする。好ましくは気温Ｔｂは日の出時の気温Ｔｅｍｐ（℃）とする。気温Ｔｅｍｐは、気象観測値ＷＸに含まれる気温データとすればよい。
【0114】
以上のことから、上記式（３Ｂ）に組み込むべき温度の補正項を、１００．０４Ｔｂとする。
【0115】
図７は、ＳＵＳ３０４についての最大値ＱＲｍａｘの対数と実測腐食量（ＲＮ）との関係を示すグラフであって、温度補正前のグラフである。また、図８は、図７のグラフに対して、温度補正を行った後のグラフである。図７および図８は、図５と同様に、日本の瀬戸内地方および沖縄県の海岸からの距離１０～１０００ｍ離れた複数の地点を試験地点Ｚとした場合の、最大値ＱＲｍａｘの対数と実測腐食量との関係を示すグラフであるが、図５に示されたプロットのほかに、新たなプロットを追加している。図７および図８における点線はそれぞれ、最大値ＱＲｍａｘの対数と実測腐食量（ＲＮ）との関係式を示す線である。関係式は、最大値ＱＲｍａｘの対数を独立変数Ｘとし、実測腐食量（ＲＮ）を従属変数ＹとするＹ＝αＸ＋βの一次関数式としている。α、βはそれぞれ定数である。なお、図５についての説明と同様に、最大値ＱＲｍａｘは必ずしも対数に変換しなくてもよい。
【0116】
図７における最大値ＱＲｍａｘの対数と実測腐食量（ＲＮ）との相関係数の２乗値（Ｒ２値）は約０．８１であり、図８におけるＲ２値は約０．８６である。図７におけるＲ２値は０．８０を超えており、温度補正を行わない場合でも関係式の信頼性は十分に高いが、温度補正を行うことによって関係式の信頼性がより高くなることがわかる。よって、温度補正によって得られる関係式に、極値推測ステップにおいて推測した腐食性指標の極大値Ｑｍａｘ（ｔＦｉ）を導入することで、利用予定地点Ｍでの将来の推定期間における鋼材の腐食量の推定値を得ることができるとともに、その信頼性を高めることができる。
【0117】
（第２の変形例）
第２の変形例では、上記式（２）におけるｓ（ｔＸｉ）および上記式（３Ａ）におけるｐ（ｔＸｉ）をそれぞれ以下のように規定した上で、それらを式（２）または式（３Ｃ）に導入し、式（２）によって二次腐食性指標Ｑ２（ｔＸｉ，ＷＸ）を求め、更に、式（３Ｃ）によって腐食性指標Ｑ（ｔＸｉ，ＷＸ）を求める。
第２の変形例では、上記式（２）におけるｓ（ｔＸｉ）および上記式（３Ａ）におけるｐ（ｔＸｉ）をそれぞれ以下のように規定した上で、それらを式（２）または式（３Ｃ）に導入し、式（２）によって二次腐食性指標Ｑ２（ｔＸｉ，ＷＸ）を求め、更に、式（３Ｃ）によって腐食性指標Ｑ（ｔＸｉ，ＷＸ）を求める。
【0118】
式（２）におけるｓ（ｔＸｉ）は腐食性蓄積係数であって、観測時刻ｔＸｉにおける鋼材表面の水分付着量ｗｓ（ｔＸｉ）が０の場合にｓ（ｔＸｉ）＝０とし、ｗｓ（ｔＸｉ）が０超の場合にｓ（ｔＸｉ）＝１とする。
式（３Ｃ）におけるｐ（ｔＸｉ）は腐食性発現係数であって、ｗｓ（ｔＸｉ）が０を超える場合にｐ（ｔＸｉ）＝０とし、ｗｓ（ｔＸｉ）が０以下の場合にｐ（ｔＸｉ）＝１、とする。
式（３Ｃ）におけるｐ（ｔＸｉ）は腐食性発現係数であって、ｗｓ（ｔＸｉ）が０を超える場合にｐ（ｔＸｉ）＝０とし、ｗｓ（ｔＸｉ）が０以下の場合にｐ（ｔＸｉ）＝１、とする。
【0119】
なお、ｗｓ（ｔＸｉ）は、観測時刻ｔＸｉにおける降雨および結露生成による水分付着量から水の蒸発量を差し引いた推測値であって、鋼材の表面温度、気温（℃）、気圧（ｈＰａ）、鋼材表面における風速（ｍ／ｓ）および相対湿度（％）から求まる値である。水分付着量ｗｓ（ｔＸｉ）は、後述する計算方法によって得られた値を用いる。鋼材表面における風速（ｍ／ｓ）は、気象観測データから得られる風速値を代用できる。
【0120】
本発明に係る腐食性指標Ｑ（ｔＸｉ，ＷＸ）は、気象変動を考慮した鋼材表面における塩分量の指標である。鋼材においては、夜間において鋼材表面に結露水が生成するとともに、塩分が結露水中に溶け込み、昼間において気温の上昇に伴って塩分を含む結露水が蒸発することで鋼板表面に塩分が濃縮されて腐食が進行する。このように、鋼材の腐食は、鋼材表面の水分の付着量が大きく影響する。よって、上記の実施形態では、鋼材表面の水分量が夜間に結露水が生成することで増加し、昼間に温度が上昇して水分が蒸発することにより減少することを前提として、観測時刻ｔＸｉ－１～ｔＸｉ間が昼間か夜間かによって、ｓ（ｔＸｉ）およびｐ（ｔＸｉ）を０または１のいずれかに決定する。
【0121】
一方、本変形例は、観測時刻ｔＸｉにおける鋼材表面の水分付着量ｗｓ（ｔＸｉ）を推測し、推測したｗｓ（ｔＸｉ）に応じて、ｓ（ｔＸｉ）およびｐ（ｔＸｉ）を０または１のいずれかに決定する。これにより、昼間か夜間かによって鋼材表面の水分の状態を推定する場合に比べて、腐食量の予測をより高い精度で行うことができる。
【0122】
また、第２の変形例では、腐食量の予測方法または腐食量予測プログラムの第４ステップにおいて、また、腐食量予測システムの第４処理部において、上記の式（３Ａ）に代えて、下記式（３Ｃ）によって、観測時刻ｔＸｉ毎に、腐食性指標Ｑ（ｔＸｉ，ＷＸ）を計算する。
【0123】
Ｑ（ｔＸｉ，ＷＸ）＝１００．０４Ｔｃ・Ｑ２（ｔＸｉ，ＷＸ）・ｐ（ｔＸｉ） …（３Ｃ）
【0124】
ただし、式（３Ｃ）におけるＴｃは、観測時刻ｔＸｉより前の時点であって鋼材表面の水分付着量ｗｓ（ｔＸｉ）が０になる時点の気温（℃）である。
【0125】
鋼材の腐食は、鋼材表面に付着する塩分量のほか、温度の影響も受ける。このため、観測地点Ｘである鋼材の利用予定地点Ｍおよび試験地点Ｚの気温差が大きくなると、腐食量の予測精度が低下する場合がある。そこで、腐食量の予測精度を高めるために、観測地点Ｘの気温の補正項をＱ（ｔＸｉ，ＷＸ）の計算式に加入することが望ましい。
【0126】
本発明者らが検討したところ、補正項を１０の累乗の形とし、また、累乗の指数は気温Ｔｃ（℃）と係数との積とすることが好ましいことがわかった。そして、鋼材の実測腐食量および鋼材の腐食性指標Ｑ（ｔＺｉ，ＷＺ）から関係式を求める際に、補正項における指数を回帰計算によって求めたところ、０．０４の値を得た。よって、温度の補正項における指数の係数は、０．０４とする。
【0127】
補正項における気温Ｔｃ（℃）は、観測時刻ｔＸｉより前の時点であって鋼材表面の水分付着量ｗｓ（ｔＸｉ）が０になる時点の気温（℃）とする。鋼材は、その表面が濡れた状態にあるときに表面の水分に塩分が溶け込み、その後、水分が蒸発して塩分濃度が高まった際に腐食が発生しやすく、また、水分が蒸発したときの気温が腐食量に影響する。そこで、補正項に加入する気温Ｔｃを、観測時刻ｔＸｉより前の時点であって鋼材表面の水分付着量ｗｓ（ｔＸｉ）が０になる時点の気温（℃）とする。この気温Ｔｃは、気象観測値ＷＸに含まれる気温データとすればよい。また、鋼材表面の水分付着量ｗｓ（ｔＸｉ）は、後述する計算方法によって得られた値を用いる。
【0128】
本変形例では、観測時刻ｔＸｉにおける鋼材表面の水分付着量ｗｓ（ｔＸｉ）を推測し、水分付着量ｗｓ（ｔＸｉ）が０になる気温Ｔｃを求める。そして、Ｔｃを用いてＱ（ｔＸｉ，ＷＸ）を温度補正する。これにより、昼間か夜間かによって鋼材表面の水分の状態を推定する場合に比べて、腐食量の予測をより高い精度で行うことができる。
【0129】
（鋼材表面の水分付着量ｗｓ（ｔＸｉ））
次に、上記第１の変形例および第２の変形例において利用する鋼材表面の水分付着量ｗｓ（ｔＸｉ）について説明する。水分付着量ｗｓ（ｔＸｉ）は、下記式（Ａ）によって導かれる。
次に、上記第１の変形例および第２の変形例において利用する鋼材表面の水分付着量ｗｓ（ｔＸｉ）について説明する。水分付着量ｗｓ（ｔＸｉ）は、下記式（Ａ）によって導かれる。
【0130】
ｗｓ（ｔＸｉ）＝（ＷＲ＋ΔＷ－Ｖａ）・Δｔ …（Ａ）
【0131】
式（Ａ）において、ＷＲは、降雨により鋼材表面に付着した水分量（ｋｇ／ｍ２／ｓ）であり、ΔＷは、結露水の生成速度（ｋｇ／ｍ２／ｓ）であり、Ｖａは水の蒸発量（ｋｇ／ｍ２／ｓ）であり、ΔｔはΔｔ＝ｔＸｉ－１－ｔＸｉ（ｓ）である。
【0132】
ΔＷは、下記式（Ｂ）によって導かれる。
【0133】
ΔＷ＝α’・（ＶＨ－ＶＨｓ） …（Ｂ）
【0134】
α’は鋼材を取り巻く空気と鋼材表面との間の湿気伝達率（ｋｇ／ｍ２／ｓ）であり、ＶＨは観測時刻ｔＸｉにおける空気中の絶対湿度（ｋｇ／ｋｇ）であり、ＶＨｓは鋼材表面の温度で決まる表面飽和絶対湿度（ｋｇ／ｋｇ）である。
【0135】
Ｖａは、下記式（Ｃ）によって導かれる。
【0136】
Ｖａ＝Ｓｈ・Ｄ・（ｃ１－ｃ２）／Ｌ …（Ｃ）
【0137】
Ｓｈはシャーウッド数であってシュミット数Ｓｃ及び鋼材表面におけるレイノルズ数Ｒｅから下記式（Ｃ－１）により導かれる。Ｄは水蒸気の拡散係数であって、観測時刻ｔＸｉにおける気温Ｔ（℃）、気圧ｐ（ｈＰａ）から下記式（Ｃ－２）により導かれる。ｃ１は観測時刻ｔＸｉの気温Ｔ（℃）での飽和水蒸気量（ｋｇ／ｍ３）であり、ｃ２は観測時刻ｔＸｉの気温Ｔ（℃）における空気中の水蒸気量（ｋｇ／ｍ３）であり、Ｌは代表長さである。レイノルズ数Ｒｅは下記式（Ｃ－３）より導かれる。ρは空気密度（ｋｇ／ｍ３）であり、ｖは鋼材表面における風速（ｍ／ｓ）であり、Ｌは特性長さであり、μは空気の粘性係数（Ｐａ／ｓ）である。
【0138】
Ｓｈ＝０．３３２・Ｒｅ１／２・Ｓｃ１／３ …（Ｃ－１）
Ｄ＝０．２４１×１０－４・｛（Ｔ＋２７３．１５）／２８８｝１．７５・１０１３．２５／（ｐ・ｔ） …（Ｃ－２）
Ｒｅ＝ρｖ２／（μｖ／Ｌ） …（Ｃ－３）
Ｄ＝０．２４１×１０－４・｛（Ｔ＋２７３．１５）／２８８｝１．７５・１０１３．２５／（ｐ・ｔ） …（Ｃ－２）
Ｒｅ＝ρｖ２／（μｖ／Ｌ） …（Ｃ－３）
【0139】
このようにして求めた水分付着量ｗｓ（ｔＸｉ）を用いて腐食性指標Ｑ（ｔＸｉ，ＷＸ）を求めることで、より高い精度での腐食量の予測を行うことができる。
【0140】
鋼材表面の水分付着量ｗｓ（ｔＸｉ）の上限は、鋼材の表面状態、鋼材が設置された状態における表面の傾斜角度などの影響を受けて変化する場合があるので、上限値を設定する場合は、鋼材に種類、傾斜状態などに応じて適宜設定してもよい。
【0141】
なお、上述の第１の変形例および第２の変形例では、時刻は地点Ｘの時刻ｔＸｉとして、気象観測値は地点Ｘにおける気象観測値ＷＸとしてそれぞれ説明したが、地点Ｘを鋼材の利用地点Ｍとする場合は、時刻ｔＸｉを利用地点Ｍでの観測時刻ｔＭｉに置き換えればよく、気象観測値ＷＸを利用地点Ｍでの気象観測値ＷＭに置き換えればよい。また、地点Ｘを試験地点Ｚとする場合は、時刻ｔＸｉを試験期間中の時刻ｔＺｉに置き換えればよく、気象観測値ＷＸを試験地点Ｚの気象観測値ＷＺに置き換えればよい。
【0142】
次に、本実施形態の鋼材の提案方法について説明する。
本実施形態の鋼材の提案方法は、先に記載の鋼材の腐食量の予測方法によって利用予定地点Ｍにおける将来にわたる鋼材の腐食量を予測する予測工程と、予測工程によって得られた腐食量の予測値を、営業担当者が顧客に提示する工程と、を備えた、鋼材の提案方法である。
また、本実施形態の鋼材の提案方法は、先に記載の鋼材の腐食量の予測方法によって利用予定地点Ｍにおける将来にわたる鋼材の腐食量を、複数種類の鋼材についてそれぞれ予測する予測工程と、予測工程によって得られた鋼材毎の腐食量の予測値を、営業担当者が顧客に提示する工程と、を備えた、鋼材の提案方法である。
腐食量の予測値は、例えば、レイティングナンバーであってもよい。レイティングナンバーは、ＪＩＳ Ｇ ０５９５：２００４（ステンレス鋼の表面さび発生程度評価方法）において規定されるものであってもよい。
本実施形態の鋼材の提案方法は、先に記載の鋼材の腐食量の予測方法によって利用予定地点Ｍにおける将来にわたる鋼材の腐食量を予測する予測工程と、予測工程によって得られた腐食量の予測値を、営業担当者が顧客に提示する工程と、を備えた、鋼材の提案方法である。
また、本実施形態の鋼材の提案方法は、先に記載の鋼材の腐食量の予測方法によって利用予定地点Ｍにおける将来にわたる鋼材の腐食量を、複数種類の鋼材についてそれぞれ予測する予測工程と、予測工程によって得られた鋼材毎の腐食量の予測値を、営業担当者が顧客に提示する工程と、を備えた、鋼材の提案方法である。
腐食量の予測値は、例えば、レイティングナンバーであってもよい。レイティングナンバーは、ＪＩＳ Ｇ ０５９５：２００４（ステンレス鋼の表面さび発生程度評価方法）において規定されるものであってもよい。
【0143】
顧客に提示する工程においては、腐食量の予測値を提示することに代えて、または予測値を提示するとともに、腐食量の予測値に対応する鋼材表面の画像を提示してもよい。具体的には、例えば、腐食量の予測値に対応する鋼材表面の画像を、鋼材の腐食量予測装置の表示部に表示させてもよく、または、書類の形で顧客に提示してもよい。
【0144】
また、顧客に提示する腐食量の予測値は、少なくとも小数点第一位までのレイティングナンバーであってもよい。同様に、腐食量の予測値に対応する鋼材表面の画像は、少なくとも小数点第一位までのレイティングナンバーに対応する画像であってもよい。
【0145】
先に記載の鋼材の腐食量の予測方法によれば、例えば、図６に示したように、各種の鋼からなる鋼材について、例えば５０年先のＲＮ値の推測値を知ることが可能である。従って、本実施形態の鋼材の提案方法によれば、利用予定地点Ｍにおける将来にわたる鋼材の腐食量を予測し、得られた腐食量の予測値を、営業担当者が顧客に提示することで、鋼材を長期間に渡り屋外環境下で使用した場合の将来の腐食量の予測値を知りたいという顧客の要望に応えることができる。また、精度の高い腐食量の予測値を提供できるので、適切な鋼材を顧客に勧めることができる。
【0146】
また、先に記載の鋼材の腐食量の予測方法によれば、例えば、図６に示したように、複数種類の鋼材について、例えば５０年先のＲＮ値の推測値を鋼種毎に知ることが可能である。従って、本実施形態の鋼材の提案方法によれば、利用予定地点Ｍにおける将来にわたる鋼材の腐食量を、複数種類の鋼材についてそれぞれ予測し、得られた鋼材毎の腐食量の予測値を、営業担当者が顧客に提示することで、精度の高い腐食量の予測値を鋼種毎に提供することができ、適切な鋼材を選択することが可能になる。これにより、適切な鋼材を顧客に勧めることができる。
【0147】
また、腐食量の予測値に対応する鋼材表面の画像を提示することで、腐食の程度を直感的に顧客に理解してもらうことができる。また、腐食量の予測値として、少なくとも小数点第一位までのレイティングナンバー、またはそれに対応する鋼材表面の画像を示すことで、高い精度の予測結果を示すことができる。
【産業上の利用可能性】
【0148】
本発明は、鋼材の腐食量の予測精度を向上させることが可能な、鋼材の腐食量の予測方法、鋼材の腐食量予測システム及び鋼材の腐食量予測プログラムを提供できる。また、本発明は、鋼材の腐食量の予測方法によって予測された鋼材の腐食量を利用した、鋼材の提案方法を提供できる。これらは、産業上の利用可能性がある。
【図1】
【図2】
【図3】
【図4A】
【図4B】
【図5】
【図6】
【図7】
【図8】