(57)【要約】【課題】未来の腐食速度を表す予測モデルを用いて土壌中の金属の腐食速度を推定する腐食速度推定装置を提供する。【解決手段】評価対象の金属を埋設した土壌に給水する給水サイクルを繰り返し、各々の給水サイクルのサイクル数と、サイクル内の一つ以上の時間と、当該時間に対応する上記の金属の腐食速度とを測定する腐食速度測定部１０と、サイクル数、サイクル内の時間、及び腐食速度を入力とし、機械学習アルゴリズムを用いて未来の腐食速度を表す予測モデルを求める学習部２０と、予測モデルに、推定したいサイクル数と当該サイクル内の時間を与え、該時間の上記の金属の腐食速度を推定する腐食速度推定部３０とを備える。【選択図】図１
2019203807.tif 000002
【特許請求の範囲】
【請求項１】
評価対象の金属を埋設した土壌に給水する給水サイクルを繰り返し、各々の給水サイクルのサイクル数と、サイクル内の一つ以上の時間と、当該時間に対応する前記金属の腐食速度とを測定する腐食速度測定部と、
前記サイクル数、前記時間、及び前記腐食速度を入力とし、機械学習アルゴリズムを用いて未来の前記腐食速度を表す予測モデルを求める学習部と、
前記予測モデルに、推定したい前記サイクル数と当該サイクル内の時間を与え、該時間の前記金属の腐食速度を推定する腐食速度推定部と
を備えることを特徴とする腐食速度推定装置。
【請求項２】
前記腐食速度測定部で測定した複数の腐食速度の変化を表す変化関数を計算し、前記学習部に与える変化関数計算部を
備え、
前記学習部は、前記変化関数を入力とし、機械学習アルゴリズムを用いて未来の前記腐食速度を表す予測モデルを求める
ことを特徴とする請求項１に記載の腐食速度推定装置。
【請求項３】
腐食速度推定装置が実行する腐食速度推定方法であって、
評価対象の金属を埋設した土壌に給水する給水サイクルを繰り返し、各々の給水サイクル内の前記金属の腐食速度を複数回測定し、
複数の前記給水サイクルのサイクル数と前記腐食速度を入力とし、機械学習アルゴリズムを用いて未来の前記腐食速度を表す予測モデルを求め、
前記予測モデルに、前記サイクル数と腐食速度を求めたいサイクル内の時間を与え、当該時間の前記金属の腐食速度を推定する
ことを特徴とする腐食速度推定方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、地中埋設された金属の腐食速度を推定する腐食速度推定装置とその方法に関する。
【背景技術】
【０００２】
インフラ設備の多くは高度経済成長期に集中的に整備されたため、今後一斉に老朽化する。その一方で、設備の維持管理に割ける費用及び保守に関わる人員は減少する傾向にある。そのため、将来的にはインフラ設備の維持管理を適切に行うことができず、安心安全の確保が困難になる恐れがある。このようなインフラ設備の老朽化問題に対処するため、近年では、設備の劣化進行を予測し、その結果に基づいて効率的な保守計画を立案する試みがなされている。
【０００３】
インフラ設備の劣化を予測する方法として、統計学的手法がある。統計学的手法は、多様な自然環境に晒され劣化したサンプルを観察・分析し、そのデータを元に回帰分析や極値統計といった統計解析を行うことで予測式を導く方法である。
【０００４】
また、実験的に予測式を求める方法もある。対象とする土壌を模擬した環境中に金属を埋設し、電気化学的手法を用いて土壌中の金属の腐食を計測する考えが提案されている（例えば非特許文献１）。
【先行技術文献】
【非特許文献】
【０００５】
【非特許文献１】宮田義一ほか１名、「電気化学的手法を中心とした土壌腐食計測（その２）」、材料と環境、46, 610〜619, 1997.
【発明の概要】
【発明が解決しようとする課題】
【０００６】
しかしながら、従来の統計学的手法では、統計解析に必要な情報がそもそも収集困難か、もしくは多大なコストと時間を掛ける必要があるという課題がある。また、実験的に土壌中の金属の腐食を計測して腐食速度を表す腐食速度関数を求めるには、実験に長期間を要する。つまり、実験的に腐食量を求める方法もコストと時間を要するという課題がある。
【０００７】
本発明は、この課題に鑑みてなされたものであり、土壌に埋設された金属の任意の経過時間における腐食速度を、簡便かつ短期間に予測することができる腐食速度推定装置とその方法を提供することを目的とする。
【課題を解決するための手段】
【０００８】
本実施形態の一態様に係る腐食速度定装置は、評価対象の金属を埋設した土壌に給水する給水サイクルを繰り返し、各々の給水サイクルのサイクル数と、サイクル内の一つ以上の時間と、当該時間に対応する前記金属の腐食速度とを測定する腐食速度測定部と、前記サイクル数、前記時間、及び前記腐食速度を入力とし、機械学習アルゴリズムを用いて未来の前記腐食速度を表す予測モデルを求める学習部と、前記予測モデルに、推定したい前記サイクル数と当該サイクル内の時間を与え、該時間の前記金属の腐食速度を推定する腐食速度推定部とを備えることを要旨とする。
【０００９】
本実施形態の一態様に係る腐食速度推定方法は、上記の腐食速度推定装置が実行する腐食速度推定方法であって、評価対象の金属を埋設した土壌に給水する給水サイクルを繰り返し、各々の給水サイクル内の前記金属の腐食速度を複数回測定し、複数の前記給水サイクルのサイクル数と前記腐食速度を入力とし、機械学習アルゴリズムを用いて未来の前記腐食速度を表す予測モデルを求め、前記予測モデルに、前記サイクル数と腐食速度を求めたいサイクル内の時間を与え、当該時間の前記金属の腐食速度を推定することを要旨とする。
【発明の効果】
【００１０】
本発明によれば、土壌に埋設された金属の任意の経過時間における腐食速度を、簡便かつ短期間に予測することができる。
【図面の簡単な説明】
【００１１】
【図１】本発明の第１実施形態に係る腐食速度推定装置の機能構成例を示す図である。
【図２】図１に示す腐食速度推定装置の動作フローを示す図である。
【図３】降雨と土壌含水率の関係を模式的に示す図である。
【図４】降雨と土壌３中の金属の腐食速度との関係を模式的に示す図である。
【図５】サイクル数を増やした場合の腐食速度の変化を模式的に示す図である。
【図６】ナイキスト線図を模式的に示す図である。
【図７】電荷移動抵抗を計算するのに仮定する等価回路の例を示す図である。
【図８】電荷移動抵抗を計算するのに仮定する等価回路の例を示す図である。
【図９】図１に示す腐食速度推定装置の腐食速度測定部で測定した腐食速度の例を模式的に示す図である。
【図１０】予測モデルに与える降雨データと、推定した腐食速度を模式的に示す図である。
【図１１】本発明の第２実施形態に係る腐食速度推定装置の機能構成例を示す図である。
【図１２】腐食速度の変化を模式的に示す図である。
【図１３】図１２に示す腐食速度の変化を直線近似によってカーブフィッティングした例を模式的に示す図である。
【図１４】図１２に示す腐食速度の変化を指数関数によってカーブフィッティングした例を模式的に示す図である。
【発明を実施するための形態】
【００１２】
以下、本発明の実施形態について図面を用いて説明する。複数の図面中同一のものに
は同じ参照符号を付し、説明は繰り返さない。
【００１３】
〔第１実施形態〕
図１は、本発明の第１実施形態に係る腐食速度推定装置の機能構成例を示す図である。図１に示す腐食速度推定装置１００は、土壌中の金属の腐食速度を、簡便かつ短期間に推定する装置である。
【００１４】
図１には、収容部２、土壌３、及びセンサ４も示す。センサ４は、評価対象の金属を含み、収容部２に収容された土壌３に埋設される。評価対象の金属は特に制限しない。例えば、鋼管柱、支持アンカー、及び地中鋼配管等に用いられる金属である。またこの他、ステンレス鋼や銅などの一般的な構造用金属でもよく、金属の種類を特に制限するものではない。
【００１５】
センサ４は、後述する3電極交流インピーダンス法の場合は評価対象の金属片を一つ、2電極交流インピーダンス法の場合は評価対象の金属片を二つ、それぞれの金属片が電気的絶縁を保持するように配置して構成される。
【００１６】
収容部２は、収容する土壌３の含有イオン成分を変化させない化学的に安定な素材で構成される。収容部２の材質は、例えばアクリルや塩ビなどのプラスチック、ガラスやセラミックスなどが好適である。
【００１７】
収容部２は、給水機能と排水機能を備える。給水機能と排水機能は、土壌３に対して十分な水を供給できることと、時間と共に土壌３中の水分量が減少できればよい。よって、給水は手動で行ってもよい。また、排水は収容部２に開口部（図示せず）が設けられていればよい。また、給水機能には、給水時間を設定可能なポンプを設けてもよい。また、排水のオン/オフを制御する開閉弁を設けてもよい。
【００１８】
土壌３は、対象とする土壌と分類上同種の土壌を用いる。現場の土壌を採取してもよいし、対象とする土壌を模擬した土壌を用いてもよい。土壌の分類には、例えば日本ペドロジー学会による日本土壌分類体系に基づく分類がある。
【００１９】
腐食速度推定装置１００は、腐食速度測定部１０、学習部２０、及び腐食速度推定部３０を備える。図２は、腐食速度推定装置１００の処理手順を示す動作フローである。図１と図２を参照してその動作を説明する。
【００２０】
腐食速度測定部１０は、評価対象の金属を含むセンサ４を埋設した土壌に給水する給水サイクルを繰り返し、各々の給水サイクルのサイクル数と、サイクル内の一つ以上の時間と、当該時間に対応するセンサ４中の金属の腐食速度とを測定する（ステップＳ１）。センサ４中の金属の腐食速度の測定方法について詳しくは後述する。以降において、「センサ４中の金属の」の文言は省略し、単に腐食速度と称する。
【００２１】
学習部２０は、腐食速度測定部１０で測定したサイクル数、時間、及び腐食速度を入力とし、機械学習アルゴリズムを用いて未来の腐食速度を表す予測モデルを求める（ステップＳ２）。機械学習アルゴリズムとは、例えばランダムフォレストや決定木等のアルゴリズムのことである。予測モデル等について詳しくは後述する。
【００２２】
腐食速度推定部３０は、学習部２０で求めた予測モデルに、腐食速度を求めたいサイクル数と、サイクル内の時間を与え、当該時間の腐食速度を推定する（ステップＳ３）。
【００２３】
以上述べた本実施形態に係る腐食速度推定装置１００によれば、給水サイクルのサイクル数、サイクル内の時間、及びサイクル内の時間の腐食速度から未来の腐食速度を表す予測モデルを学習し、当該予測モデルを用いて腐食速度を推定するので、土壌中の金属の腐食速度を、簡便かつ短期間に予測することができる。
【００２４】
土壌はSiやAl，Ti，Fe，Caなどの酸化物等からなる土壌粒子と、前記土壌粒子の間隙中に存在する気相および液相（水）から構成される3相共存環境である。土壌中の気相の割合と液相の割合の合計は一定と考えることができ、一方が高まればもう一方が低くなる相反関係にある。ここで、前記土壌中での金属の腐食反応には、基本的に水と酸素が必要であり、これらの状態に依存した腐食速度で腐食が進行する。従って、土壌中を占める水の割合を指す土壌含水率は腐食速度に寄与する主要な環境因子であり、土壌含水率と共に腐食速度は変化するといえる。
【００２５】
土壌含水率は、地中のよほど深い位置でない限り、常に一定に保持されているわけではない。土壌含水率は、例えば降雨などの自然現象に応じて変化する。
【００２６】
図３は、降雨と土壌含水率の関係を模式的に示す図である。図３の横軸は経過時間である。図３に示すように、土壌含水率の増減は降雨とよく連動しており，降雨時に急激に増加し，雨が止むと徐々に減少するというサイクルを繰り返す．従って，腐食速度の経時的な変化もまた，降雨を起点としたサイクルの繰り返しになると考えることができる。
【００２７】
図４は、降雨と土壌３中の金属の腐食速度との関係を模式的に示す図である。図４の横軸は経過時間である。ここで、１サイクルとは、降雨から次ぎの降雨までの期間を指しており、降雨間隔によって１サイクルの時間の長さは異なる。
【００２８】
１サイクルにおける腐食速度の時間変化は、土壌と金属によって異なる。また、同じ土壌と金属の組み合わせでも、何回目のサイクルであるかによって腐食速度は異なる。例えば鋼等の金属では、埋設してから10回目のサイクルにおける腐食速度の時間変化と、100回目の腐食速度とでは異なる。一般的に、サイクル数を重ねるほど腐食速度は低下する。
【００２９】
図５は、サイクル数を増やした場合の腐食速度の変化を模式的に示す図である。図５の横軸は経過時間である。
【００３０】
図５に示すように100回目のサイクルでは、10回目のサイクルよりも一般的に腐食速度は低下する。これは、腐食の進行に伴い土壌中の金属の表面に酸化被膜（錆）が形成され、この錆が腐食反応を抑制するためである。
【００３１】
このようにサイクル数を増やすごとに、金属表面を覆う錆の面積と厚みが増加するため腐食速度は低下する。ただし、金属表面に形成される錆の状態は、その金属が晒されて来た降雨パターンによって異なる。したがって、例えば一定期間における降雨頻度が多い場合と少ない場合とでは、同じサイクル数を経た金属でも腐食速度及び腐食量は異なる。
【００３２】
本実施形態に係る腐食速度推定装置１００が実行する腐食速度推定方法は、降雨から次ぎの降雨までの１サイクルを基準に、ｎ回目のサイクルにおける腐食速度の時間変化ｆｎ（ｔ）は、金属が晒された1回目〜ｎ−1回目までの降雨パターン（各サイクルの降雨間隔及び降雨量や連続降雨時間等の組み合わせ）、及びｎ回目の降雨間隔と降雨量や連続降雨時間等に依存するという考えに基づき、土壌中の金属の腐食速度を推定する方法である。
【００３３】
つまり、本実施形態に係る腐食速度推定装置１００は、評価対象の金属を埋設した土壌に給水する給水サイクルを繰り返し、各々の給水サイクルのサイクル数と、サイクル内の一つ以上の時間と、当該時間に対応す上記の金属の腐食速度とを測定する腐食速度測定部１０と、サイクル数、サイクル内の時間、及び腐食速度を入力とし、機械学習アルゴリズムを用いて未来の腐食速度を表す予測モデルを求める学習部２０と、予測モデルに、サイクル数と腐食速度を求めたいサイクル内の時間を与え、当該時間の金属の腐食速度を推定する。これによれば、腐食速度の推定に必要な複数サイクルにおける腐食速度の時間変化は、低コスト、短時間で取得可能であるため、任意の未来の腐食速度を簡便且つ短時間で予測することが可能になる。
【００３４】
次に、腐食速度推定装置１００の各機能構成部について詳しく説明する。
【００３５】
（腐食速度測定部）
腐食速度測定部１０は、交流インピーダンス測定機能を有する。交流インピーダンス測定は、センサ４中の金属を電極とし、電極間に微少電圧もしくは電流を交流で印加し、電気的な応答を測定する。印加する電圧もしくは電流は、センサ４中の金属の表面が変化しないように微少にするのがよい。例えば5mV程度である。周波数は、例えば0.1Hz〜500Hzの幅で変化させる。
【００３６】
交流インピーダンス測定を行うことでナイキスト線図を得ることができる。図６に、ナイキスト線図を模式的に示す。ナイキスト線図の横軸は実部、縦軸軸は虚部である。ナイキスト線図を元に、所定の等価回路に基づいてカーブフィッティングすることで電荷移動抵抗を導く。
【００３７】
図７と図８は、電荷移動抵抗を計算するのに仮定する等価回路の例である。何れの図も（ａ）は、3電極で交流インピーダンスを測定した場合の等価回路である。（ｂ）は、2電極で交流インピーダンスを測定した場合の等価回路である。
【００３８】
図中の電荷移動抵抗Ｒｃｔは、土壌３中のセンサ４中の金属の腐食反応の抵抗を表す。電気二重層Ｃｄｌは、センサ４中の金属と土壌３の界面に存在する容量である。抵抗成分Ｒｓ１，Ｒｓ２は、土壌３中及びその他の抵抗を表す。容量Ｃｓは、土壌３の容量成分である。ワールブルグインピーダンスＺｗは、拡散過程によるインピーダンスである。なお、カーブフィッティングする際は、電気二重層Ｃｄｌと容量ＣＳはＣＰＥ（Constant Phase Element）に置き代えてもよい。
【００３９】
図７と図８に示す等価回路によれば、図６に示すようにナイキスト線図上に理論上二つの円弧が描かれる。高周波数側の円弧は土壌３に由来する。低周波数側の円弧は腐食反応に起因するものである。
【００４０】
電荷移動抵抗は、ナイキスト線図の低周波数側の円弧が横軸（実部）と交差する幅で与えられる。なお、2電極で交流インピーダンスを測定した場合の電荷移動抵抗は、その幅の半分の値である。
【００４１】
腐食速度は、電荷移動抵抗の逆数に比例する。その比例係数Ｋは、アノード及びカソード分極曲線からＴａｆｅｌ勾配を導いて算出する（参考文献：「コンクリート中鋼材の腐食速度測定方法（CIPE法）の開発）」、さび、148号、p2-8,2015）。
【００４２】
比例係数Ｋを用いることで、電荷移動抵抗の逆数から、腐食電流密度（腐食速度）を算出することができる。腐食電流密度を用いて重量減肉速度及び体積減肉速度等を算出することもできる。
【００４３】
このように測定した1回のインピーダンス測定の結果から1つの腐食速度を得る。腐食速度測定部１０は、1回の給水サイクルにおいて複数回のインピーダンス測定を行う。そして、その給水サイクルを繰り返す。
【００４４】
図９は、腐食速度測定部１０で測定した腐食速度の例を模式的に示す図である。図８の横軸は経過時間、縦軸は腐食速度である。
【００４５】
このように、腐食速度測定部１０は、評価対象の金属を埋設した土壌３に給水する給水サイクルを繰り返し、各々の給水サイクルのサイクル数と、サイクル内の一つ以上の時間と、当該時間に対応する金属の腐食速度とを測定する。
【００４６】
（学習部）
学習部２０は、腐食速度測定部１０で測定したサイクル数、サイクル内の時間、及び当該時間の腐食速度を入力とし、機械学習アルゴリズムを用いて未来の腐食速度を表す予測モデルを求める。
【００４７】
予測モデルとは、予測に使う情報と予測したい情報を結ぶモデル（例えば式）である。予測したい情報は、腐食速度である。また、予測に使う情報は、腐食速度測定部１０で測定した各情報（サイクル数、サイクル内の時間、及び当該時間の腐食速度）である。
【００４８】
学習部２０の機械学習方法（アルゴリズム）は特に制限しない。例えば、ランダムフォレストや決定木などのアルゴリズムを用いることが考えられる。
【００４９】
上記のように土壌中の金属の腐食速度は、降雨を基点としたサイクルの繰り返しに伴い変化する。あるサイクルでの腐食速度の時間変化は、それ以前の全てのサイクルの変化に依存する。いわゆる自己相関の強い変化を示す。
【００５０】
したがって、予測モデルを求める説明変数に、例えば、埋設後からの経過時間ｔ、ｎ回目のサイクルの腐食速度の時間変化ｆｎ（ｔ）、ｎ回目の連続降雨時間Ｗｎ、ｎ回目のサイクルの時間間隔Ｔｎを用いることができる。また、説明変数には、これらの組み合わせを用いてもよい。
【００５１】
（腐食速度推定部）
腐食速度推定部３０は、学習部２０で求めた予測モデルに、推定したいサイクル数と当該サイクル内の時間を与え、該時間の腐食速度を推定する。
【００５２】
予測モデルに与える情報は、予測モデルを使って推定するのに必要な情報を含む。例えば、推定したい期間における降雨データである。
【００５３】
図１０は、予測モデルに与える降雨データと、推定した腐食速度を模式的に示す図である。このように腐食速度の時間変化を予測することができる。
【００５４】
なお、腐食速度を時間積分することにより、腐食量も容易に求めることができる。
【００５５】
〔第２実施形態〕
図１１は、本発明の第２実施形態に係る腐食速度推定装置の機能構成例を示す図である。図１１に示す腐食速度推定装置２００は、腐食速度推定装置１００（図１）に対して変化関数計算部４０を備える点で異なる。
【００５６】
変化関数計算部４０は、腐食速度測定部１０で測定した複数の腐食速度の変化を表す変化関数を計算し、学習部２０に与える。
【００５７】
変化関数計算部４０は、例えば図１２に示す様な腐食速度の変化から、一つ又は複数の変化関数を計算する。つまり、変化関数計算部４０は、経過時間に対する腐食速度の変化を表す変化関数を計算する。
【００５８】
変化関数は、例えば図１２に示す腐食速度の経過時間に対する変化を表す任意の関数Ｖ（ｔ）である。関数Ｖ（ｔ）は、例えば図１２に示す腐食速度の変化に、任意の関数をカーブフィッティングして求める。
【００５９】
図１３は、図１２に示す腐食速度の変化を直線近似によってカーブフィッティングした例を模式的に示す。給水開始から腐食速度が上昇し、ほぼ安定値になるまでの関数はＶ１（ｔ）＝ａｔ、安定値になってからの関数はＶ２（ｔ）＝ｂである。また、給水を停止してから腐食速度が低下する部分の関数はＶ３（ｔ）＝ｃｔである。
【００６０】
直線の他、二次関数、指数関数を用いてもよい。図１４は、Ｖ３（ｔ）の代わりにＶ４（ｔ）＝c×exp(dt)+bとした図である。なお、ａ，ｂ，ｃ，ｄは定数である。
【００６１】
このように求めた変化関数を学習部２０に入力するようにしてもよい。つまり、変化関数の予測モデルを求め、その予測モデルから未来の腐食速度を推定するようにしてもよい。
【００６２】
以上説明したように本実施形態の腐食速度推定装置１００，２００によれば、腐食速度の推定に必要な複数サイクルにおける腐食速度の時間変化は、低コスト、短時間で取得可能であるため、任意の未来の腐食速度を簡便且つ短時間で予測することが可能になる。
【００６３】
本発明は、上記の実施形態に限定されるものではなく、その要旨の範囲内で変形が可能である。例えば、腐食速度推定部３０は、腐食速度を時間積分して腐食量を推定するようにしてもよい。
【００６４】
なお、学習部２０のアルゴリズムとしてランダムフォレストと決定木を例に挙げたが、これ以外のアルゴリズムを用いても構わない。要するに、学習部２０は、予測に使う情報（説明変数）がどのようなもので有っても、未来の腐食速度（目標関数）を表す予測モデルが推定出来ればどのようなアルゴリズムを用いてもよい。
【００６５】
このように、本発明はここでは記載していない様々な実施形態等を含むことは勿論である。したがって、本発明の技術的範囲は上記の説明から妥当な特許請求の範囲に係る発明特定事項によってのみ定められるものである。
【符号の説明】
【００６６】
１００、２００：腐食速度推定装置
２：収容部
３：土壌
４：センサ（センサ中の金属）
１０：腐食速度測定部
２０：学習部
３０：腐食速度推定部
４０：変化関数計算部
【図１】
2019203807.tif 000003
【図２】
2019203807.tif 000004
【図３】
2019203807.tif 000005
【図４】
2019203807.tif 000006
【図５】
2019203807.tif 000007
【図６】
2019203807.tif 000008
【図７】
2019203807.tif 000009
【図８】
2019203807.tif 000010
【図９】
2019203807.tif 000011
【図１０】
2019203807.tif 000012
【図１１】
2019203807.tif 000013
【図１２】
2019203807.tif 000014
【図１３】
2019203807.tif 000015
【図１４】
2019203807.tif 000016