金属材料的腐蚀量测绘方法、金属材料的选定方法及金属材料的腐蚀量测绘装置与流程

1.本发明涉及金属材料的腐蚀量测绘方法、金属材料的选定方法及金属材料的腐蚀量测绘装置。


背景技术：

2.如非专利文献1所示，以往，关于大气腐蚀环境中的金属材料的腐蚀量，作为经验式已知使用以下的式(1)表示。
3.[数学式1]
[0004]
y＝axb…
(1)
[0005]
在此，在上述式(1)中，y是金属材料的腐蚀量，x是金属材料的使用期间，a是表示金属材料的自初始起一年期间的腐蚀量的参数，b是表示由腐蚀形成的锈层的效果引起的腐蚀速度的衰减的参数。这些参数a、b的值根据金属材料的种类、大气腐蚀环境而变化。因此，目前，在预测长期腐蚀量时，大多使用在成为对象的大气腐蚀环境下将金属材料暴露多个期间，用上述式(1)外推腐蚀量的经时变化的方法。
[0006]
然而，金属材料的腐蚀量由金属具有的耐腐蚀性能和大气腐蚀环境因素、例如温度、相对湿度、润湿时间、降雨量、空气含盐量(amount of airbornesea salt)及so2浓度等复杂地作用而决定。因此，如下所示，提出了通过使用上述环境因素进行公式化来预测金属材料的腐蚀量的技术。
[0007]
例如在非专利文献2中，对于碳钢、锌、铜及铝，通过以温度、相对湿度及飞来盐的对数进行多元回归而得的项与以温度、相对湿度及so2量的对数进行多元回归而得的项之和来计算腐蚀量的对数。
[0008]
此外，在专利文献1中，将年润湿时间、年平均风速、空气含盐量、硫氧化物量、腐蚀反应的活化能和温度作为参数，如下述式(2)那样表示腐蚀指标z，将该腐蚀指标z的二次函数用于长期腐蚀量的计算。
[0009]
[数学式2]
[0010][0011]
在此，在上述式(2)中，tow为年润湿时间(h)，w为年平均风速(m/s)，c为空气含盐量(mdd)，s为硫氧化物量(mdd)，e
α
为腐蚀反应的活化能(j/mol)，r为气体常量(j/(k/mol))，t为年平均气温(k)，α、κ、δ、ε为常量。此外，上述的[mdd]是每单位天数和单位面积收集的nacl的量，是[mg nacl
·
dm
－2
·
day
－1
]的缩写。
[0012]
此外，在专利文献2中，提出了预测空气含盐量c的技术。此外，在专利文献3中，提出了以温度、相对湿度、空气含盐量及润湿概率作为参数而通过下述式(3)算出将作为经验式而已知的上述式(1)的参数a，将参数b通过实验室试验作为参数a的函数而算出的腐蚀预测技术。
[0013]
[数学式3]
[0014]
a＝(α
·
t β)
·
pw(t，h)
·
(say)
…
(3)
[0015]
在此，上述式(3)中，t为温度(℃)，h为相对湿度(％)，sa为空气含盐量(mg/dm2/day(＝mdd))，pw(t，h)为润湿概率，α、β、γ为根据钢种设定的系数。
[0016]
此外，专利文献4中提出了在室外的大气腐蚀环境中，以温度、润湿时间和空气含盐量为参数，通过下述式(4)计算作为经验式而已知的上述式(1)的参数a，将参数b设定在0.3～0.6的范围的腐蚀预测技术。
[0017]
[数学式4]
[0018]
a＝kt
α
·
tow
β
·
say…
(4)
[0019]
这里，上述式(4)中，t为温度(℃)，tow为润湿时间(h)，sa为空气含盐量(mg/dm2/day(＝mdd))，α、β、γ为系数。
[0020]
此外，在专利文献5中，提出了在预测大气环境下的钢材的板厚减少量时，根据下述式(5)算出作为经验式而已知的上述式(1)的参数a，根据下述式(6)算出参数b的技术。
[0021]
[数学式5]
[0022]
a＝(cr0 cr1)
÷2…
(5)
[0023]
[数学式6]
[0024]
b＝2cr1
÷
(cr0 cr1)
…
(6)
[0025]
在此，上述式(5)及上述式(6)中，cr0是表示以环境因素为参数的钢材刚制造后的初始腐蚀速度的函数，cr1是表示以环境因素为参数的钢材制造后1年后的腐蚀速度的函数。此外，这里的环境因素表示年平均温度(℃)、年平均湿度(％)、年平均风速(m/sec)、空气含盐量(mg/dm2/day(＝mdd))、硫氧化物量(mg/dm2/day(＝mdd))。
[0026]
此外，在专利文献6中，在预测金属材料的腐蚀速度时，进行以腐蚀速度为目标变量、以对该腐蚀速度造成影响的环境因素和地形因素为解释变量的多元回归分析。而且，在进行该多元回归分析时，至少作为解释变量之一，包含根据相对湿度0％～100％进行了加权而得的假想润湿时间。此外，提出了如下方法：将根据变化的相对湿度而不同的加权系数乘以与变化的相对湿度对应的时间而得到的乘积值进行总和而求出该假想润湿时间，基于测定到的金属材料的腐蚀速度，通过多元回归分析法制作腐蚀速度推定式。
[0027]
此外，在专利文献7中提出了如下技术：作为测绘(地图化)腐蚀量的预测结果的技术，提出了通过使用逐级法来优化的多元回归分析来预测金属材料的腐蚀速度，通过聚类来测绘该金属材料的腐蚀速度。
[0028]
[现有技术文献]
[0029]
[专利文献]
[0030]
[专利文献1]日本专利第3909057号公报
[0031]
[专利文献2]日本专利第4143018号公报
[0032]
[专利文献3]日本专利第4706254号公报
[0033]
[专利文献4]日本专利第5895522号公报
[0034]
[专利文献5]日本专利第5066160号公报
[0035]
[专利文献6]日本专利第5066955号公报
[0036]
[专利文献7]日本专利第5684552号公报
[0037]
[非专利文献]
[0038]
[非专利文献1][关于耐候性钢材在桥梁上的应用的共同研究报告书(xviii)]，建设省土木研究所，(公司)钢材俱乐部，(公司)日本桥梁建设协会，平成5年3月
[0039]
[非专利文献2]iso9223：1992“corrosion of metals and alloys-corrosivity of atmospheres-classification，determination and estimation”

技术实现要素：

[0040]
[发明要解决的课题]
[0041]
在此，在专利文献1－5中，通过评价腐蚀量和腐蚀速度与各环境参数之间的关系，选定作为解释变量采用的环境参数并创建公式。但是，腐蚀量、腐蚀速度以及各环境参数具有复杂的相关关系。例如，腐蚀量与温度的关系为非线性，空气含盐量与so2浓度具有近似相关。在这样的关系中，在如专利文献1～5那样创建公式的情况下，不能期望高精度的预测。
[0042]
此外，专利文献6的特征在于，着眼于润湿时间进行加权。通过这样进行加权，虽然精度提高，但即使在多个环境参数中只对润湿时间进行加权，也不能期望有效的精度提高。此外，在专利文献6的方法中，仅在保有当前数据的期间能够进行腐蚀预测，无法进行长期腐蚀预测。
[0043]
此外，在专利文献7中，通过聚类对数据进行分类，针对每个类从多个环境参数中选定多个环境参数，反复进行多元回归而成为最佳的式子，由此实现预测精度提高。但是，腐蚀量及腐蚀速度的关系无法用与由多元回归得到的各环境参数的1次式简单地表现。因此，在专利文献7的方法中，不能期望较大的精度提高。
[0044]
此外，在专利文献7所公开的方法中，从公共提供的数据(例如农研机构提供的[农业气象数据]等)等中取得预测地点的环境参数，例如在日本国内的情况下取得1km见方的各个网格的环境参数。但是，在这样的数据中，空气含盐量、so2量等所有的环境参数没有被网罗。此外，在专利文献7所公开的方法中，也没有明示不足的环境参数的补充方法。因此，认为即使使用专利文献7中公开的方法，仅用已知的环境参数进行聚类和腐蚀量的预测，也只能制作精度低的腐蚀量预测地图。
[0045]
本发明是鉴于上述情况而完成的，其目的在于提供一种在大气腐蚀环境中能够高精度地测绘金属材料的腐蚀量的金属材料的腐蚀量测绘方法、金属材料的选定方法及金属材料的腐蚀量测绘装置。
[0046]
[用于解决课题的手段]
[0047]
为了解决上述课题、达成目的，本发明的金属材料的腐蚀量测绘方法，其使用腐蚀量数据来预测金属材料的腐蚀量，并制作腐蚀量预测地图，所述腐蚀量数据包含：金属材料的使用期间；表示所述使用期间内的所述金属材料的使用环境且在表示使用所述金属材料的地域的地图上已知的多个环境参数；所述地图上的所述环境参数的位置坐标；所述地图的地形数据；以及所述使用期间内的所述金属材料的腐蚀量，所述金属材料的腐蚀量测绘方法的特征在于，包括：环境地图制作步骤，根据所述多个环境参数、所述地图的地形数据、所述地图上的所述环境参数的位置坐标，以任意的网格间隔制作每个所述环境参数的环境地图；预测请求点输入步骤，输入预测请求点，该预测请求点包含制作所述腐蚀量预测地图
的所述金属材料的使用期间和所述腐蚀量数据中的多个环境参数；相似度算出步骤，算出所述腐蚀量数据中的多个环境参数与所述预测请求点中的多个环境参数的相似度；维数压缩步骤，考虑所述相似度来将所述腐蚀量数据中的多个环境参数维数压缩为潜在变量；腐蚀量预测步骤，基于使用所述潜在变量及所述相似度构建的预测公式，预测在所述网格的预测请求点处的所述金属材料的腐蚀量；以及腐蚀量预测地图制作步骤，通过将所述网格的预测请求点处的腐蚀量的预测结果在地图上进行涂色标注，由此制作腐蚀量预测地图。
[0048]
此外，本发明的金属材料的腐蚀量测绘方法的特征在于，在上述方案中，所述环境地图制作步骤包括：第一步骤，将所述多个环境参数输入到所述地图上的位置坐标；第二步骤，根据所述地图的地形数据，将所述环境地图修正为海拔0m的环境地图；第三步骤，在所述海拔0m的环境地图中，在所述多个环境参数之间补充环境参数；第四步骤，基于所述地图的地形数据，将补充了环境参数的海拔0m的环境地图修正为原来海拔的环境地图。
[0049]
此外，本发明的金属材料的腐蚀量测绘方法的特征在于，在上述方案中，所述多个环境参数包括空气含盐量，所述第三步骤中，在补充靠近海岸线的网格的空气含盐量时，以不超过预先设定的上限值的方式补充空气含盐量。
[0050]
此外，本发明的金属材料的腐蚀量测绘方法的特征在于，在上述方案中，所述腐蚀量预测步骤包括：初始腐蚀量预测步骤，预测表示所述金属材料的规定期间的腐蚀量的第一参数；衰减预测步骤，预测表示所述金属材料的腐蚀速度衰减的第二参数；以及长期腐蚀量预测步骤，基于所述金属材料的使用期间、所述第一参数和所述第二参数，预测所述金属材料的比所述规定期间长的期间的腐蚀量。
[0051]
此外，本发明的金属材料的腐蚀量测绘方法的特征在于，在上述方案中，所述初始腐蚀量预测步骤中，基于使用所述潜在变量及所述相似度构建的预测公式，预测所述预测请求点的环境参数下的所述金属材料的规定期间的腐蚀量。
[0052]
此外，本发明的金属材料的腐蚀量测绘方法的特征在于，在上述方案中，所述衰减预测步骤中，基于制作所述腐蚀量预测地图的所述金属材料的使用期间、所述腐蚀量数据中的多个环境参数和所述相似度，预测所述第二参数。
[0053]
此外，本发明的金属材料的腐蚀量测绘方法的特征在于，在上述方案中，所述多个环境参数包括：温度、相对湿度、绝对湿度、润湿时间及降雨量中的至少一个；和空气含盐量、so
x
浓度及no
x
浓度中的至少一个。
[0054]
此外，本发明的金属材料的腐蚀量测绘方法的特征在于，在上述方案中，所述金属材料为钢铁材料。
[0055]
此外，本发明的金属材料的腐蚀量测绘方法的特征在于，在上述方案中，所述腐蚀量预测地图为根据在所述腐蚀量预测地图制作步骤中预测出的腐蚀量的值对每个所述网格进行涂色标注而成的地图。
[0056]
为了解决上述课题、达成目的，本发明的金属材料的选定方法的特征在于，使用上述的金属材料的腐蚀量测绘方法，选定与使用环境相应的金属材料。
[0057]
为了解决上述课题、达成目的，本发明的金属材料的腐蚀量测绘装置，包括：数据库，其保存腐蚀量数据，该腐蚀量数据包括：金属材料的使用期间、表示所述使用期间内的所述金属材料的使用环境且在表示使用所述金属材料的地域的地图上已知的多个环境参数、所述地图上的所述环境参数的位置坐标、所述地图的地形数据、以及所述使用期间内的
所述金属材料的腐蚀量；环境地图制作部，其根据所述多个环境参数、所述地图的地形数据、所述地图上的所述环境参数的位置坐标，以任意的网格间隔，制作每个所述环境参数的环境地图；输入部，其输入预测请求点，该预测请求点包括：制作腐蚀量预测地图的所述金属材料的使用期间和所述腐蚀量数据中的多个环境参数；相似度算出部，其算出所述腐蚀量数据中的多个环境参数与所述预测请求点处的多个环境参数的相似度；维数压缩部，其考虑所述相似度来将所述腐蚀量数据中的多个环境参数维数压缩为潜在变量；腐蚀量预测部，其基于使用所述潜在变量及所述相似度构建的预测公式，预测所述网格的预测请求点处的所述金属材料的腐蚀量；以及腐蚀量预测地图制作部，其通过将所述网格的预测请求点处的腐蚀量的预测结果在地图上进行涂色标注，由此制作所述腐蚀量预测地图。
[0058]
本发明是基于以上见解而完成的，其主旨如下。
[0059]
[1]金属材料的腐蚀量测绘方法，其使用腐蚀量数据来预测金属材料的腐蚀量，并制作腐蚀量预测地图，所述腐蚀量数据包含：金属材料的使用期间；表示所述使用期间内的所述金属材料的使用环境且在表示使用所述金属材料的地域的地图上已知的多个环境参数；所述地图上的所述环境参数的位置坐标；所述地图的地形数据；以及所述使用期间内的所述金属材料的腐蚀量，所述金属材料的腐蚀量测绘方法的特征在于，包括：
[0060]
环境地图制作步骤，根据所述多个环境参数、所述地图的地形数据、所述地图上的所述环境参数的位置坐标，以任意的网格间隔制作每个所述环境参数的环境地图；
[0061]
预测请求点输入步骤，输入预测请求点，该预测请求点包含制作所述腐蚀量预测地图的所述金属材料的使用期间和所述腐蚀量数据中的多个环境参数；
[0062]
相似度算出步骤，算出所述腐蚀量数据中的多个环境参数与所述预测请求点中的多个环境参数的相似度；
[0063]
维数压缩步骤，考虑所述相似度来将所述腐蚀量数据中的多个环境参数维数压缩为潜在变量；
[0064]
腐蚀量预测步骤，基于使用所述潜在变量及所述相似度构建的预测公式，预测在所述网格的预测请求点处的所述金属材料的腐蚀量；以及
[0065]
腐蚀量预测地图制作步骤，通过将所述网格的预测请求点处的腐蚀量的预测结果在地图上进行涂色标注，由此制作腐蚀量预测地图。
[0066]
[2]根据[1]所述的金属材料的腐蚀量测绘方法，其中，
[0067]
所述环境地图制作步骤包括：
[0068]
第一步骤，将所述多个环境参数输入到所述地图上的位置坐标；
[0069]
第二步骤，根据所述地图的地形数据，将所述环境地图修正为海拔0m的环境地图；
[0070]
第三步骤，在所述海拔0m的环境地图中，在所述多个环境参数之间补充环境参数；
[0071]
第四步骤，基于所述地图的地形数据，将补充了环境参数的海拔0m的环境地图修正为原来海拔的环境地图。
[0072]
[3]根据[2]所述的金属材料的腐蚀量测绘方法，其中，
[0073]
所述多个环境参数包括空气含盐量，
[0074]
所述第三步骤中，在补充靠近海岸线的网格的空气含盐量时，以不超过预先设定的上限值的方式补充空气含盐量。
[0075]
[4]根据[1]～[3]中任一项所述的金属材料的腐蚀量测绘方法，其中，所述腐蚀量
预测步骤包括：
[0076]
初始腐蚀量预测步骤，预测表示所述金属材料的规定期间的腐蚀量的第一参数；
[0077]
衰减预测步骤，预测表示所述金属材料的腐蚀速度衰减的第二参数；以及
[0078]
长期腐蚀量预测步骤，基于所述金属材料的使用期间、所述第一参数和所述第二参数，预测所述金属材料的比所述规定期间长的期间的腐蚀量。
[0079]
[5]根据[4]所述的金属材料的腐蚀量测绘方法，其中，
[0080]
所述初始腐蚀量预测步骤中，基于使用所述潜在变量及所述相似度构建的预测公式，预测所述预测请求点的环境参数下的所述金属材料的规定期间的腐蚀量。
[0081]
[6]根据[4]所述的金属材料的腐蚀量测绘方法，其中，
[0082]
所述衰减预测步骤中，基于制作所述腐蚀量预测地图的所述金属材料的使用期间、所述腐蚀量数据中的多个环境参数和所述相似度，预测所述第二参数。
[0083]
[7]根据[1]～[6]中任一项所述的金属材料的腐蚀量测绘方法，其中，
[0084]
所述多个环境参数包括：温度、相对湿度、绝对湿度、润湿时间及降雨量中的至少一个；和空气含盐量、sox浓度及nox浓度中的至少一个。
[0085]
[8]根据[1]～[4]中任一项所述的金属材料的腐蚀量测绘方法，其中，
[0086]
所述金属材料为钢铁材料。
[0087]
[9]根据[1]～[8]中任一项所述的金属材料的腐蚀量测绘方法，其中，
[0088]
所述腐蚀量预测地图为根据在所述腐蚀量预测地图制作步骤中预测出的腐蚀量的值对每个所述网格进行涂色标注而成的图。
[0089]
[10]金属材料的选定方法，其特征在于，使用[1]～[8]中任一项所述的金属材料的腐蚀量测绘方法，选定与使用环境相应的金属材料。
[0090]
[11]金属材料的腐蚀量测绘装置，其包括：
[0091]
数据库，其保存腐蚀量数据，该腐蚀量数据包括：金属材料的使用期间、表示所述使用期间内的所述金属材料的使用环境且在表示使用所述金属材料的地域的地图上已知的多个环境参数、所述地图上的所述环境参数的位置坐标、所述地图的地形数据、以及所述使用期间内的所述金属材料的腐蚀量；
[0092]
环境地图制作部，其根据所述多个环境参数、所述地图的地形数据、所述地图上的所述环境参数的位置坐标，以任意的网格间隔，制作每个所述环境参数的环境地图；
[0093]
输入部，其输入预测请求点，该预测请求点包括：制作腐蚀量预测地图的所述金属材料的使用期间和所述腐蚀量数据中的多个环境参数；
[0094]
相似度算出部，其算出所述腐蚀量数据中的多个环境参数与所述预测请求点处的多个环境参数的相似度；
[0095]
维数压缩部，其考虑所述相似度来将所述腐蚀量数据中的多个环境参数维数压缩为潜在变量；
[0096]
腐蚀量预测部，其基于使用所述潜在变量及所述相似度构建的预测公式，预测所述网格的预测请求点处的所述金属材料的腐蚀量；以及
[0097]
腐蚀量预测地图制作部，其通过将所述网格的预测请求点处的腐蚀量的预测结果在地图上进行涂色标注，由此制作所述腐蚀量预测地图。
[0098]
[发明的效果]
[0099]
根据本发明，在大气腐蚀环境中，能够高精度地测绘金属材料的腐蚀量，能够选定具有与使用环境相应的耐腐蚀性的最佳的金属材料。
附图说明
[0100]
[图1]图1是表示在大气腐蚀环境中温度(年平均)与腐蚀量(年间)的关系的曲线图。
[0101]
[图2]图2是表示在大气腐蚀环境中so2浓度(年平均)与空气含盐量(年平均)的模拟相关关系的曲线图。
[0102]
[图3]图3是表示本发明的实施方式的金属材料的腐蚀量测绘装置的结构的框图。
[0103]
[图4]图4是表示本发明的实施方式的金属材料的腐蚀量测绘方法的流程的流程图。
[0104]
[图5]图5是通过现有技术的金属材料的腐蚀量测绘方法得到的表示金属材料的一年后的腐蚀量预测值的腐蚀量预测地图。
[0105]
[图6]图6是表示通过本发明的实施方式的金属材料的腐蚀量测绘方法得到的金属材料的五十年后的腐蚀量预测值的腐蚀量预测地图。
具体实施方式
[0106]
以下，参照附图对本发明的实施方式的金属材料的腐蚀量测绘方法、金属材料的选定方法及金属材料的腐蚀量测绘装置进行说明。需要说明的是，本发明并不限定于以下的实施方式。
[0107]
(环境参数的补充)
[0108]
大气中的金属材料、腐蚀量和各种环境参数之间存在复杂的相关关系，在上述非专利文献2的国际标准iso9223中，基于年平均的温度、相对湿度、空气含盐量、so
x
浓度，将年度的腐蚀量公式化。但是，例如在地图上以1km、2km左右的窄的网格间隔获得这些环境参数是困难的，因此，需要基于能够获得的环境参数的位置坐标来补充各环境参数的位置坐标间的环境参数。
[0109]
此时，例如温度不仅根据经度、纬度而变化，还根据标高而变化。此外，飞来盐分的飞来方式也因山或山丘等障碍物的影响而变化。因此，需要考虑这样的地形的影响来补充环境参数。
[0110]
但是，难以通过一次运算来进行地形数据的反映和在位置坐标之间的补充环境参数。因此，在本实施方式中，首先将标绘有已知的多个环境参数的多个环境地图修正为以海拔0m换算的多个环境地图。接着，分别制作在位置坐标之间的补充环境参数后的海拔0m的环境地图。接着，通过结合实际的地形数据修正补充后的环境参数，分别生成原来海拔的环境地图。需要说明的是，上述的“环境地图”表示标绘有环境参数的地图。以下，对每个环境参数的具体的补充方法进行说明。
[0111]
(温度的补充)
[0112]
在温度的补充中，首先将可获得的温度数据标绘在地图上。接着，基于“若标高下降100m则气温上升0.6℃”这一法则(气温减率)，使用地形数据，将各标绘点的温度的值修正为换算成海拔0m的值。接着，以任意的网格间隔补充数据间的温度。需要说明的是，网格
间隔可以根据地图上已知的环境参数的数量、计算时的机器规格等来决定。
[0113]
数据的补充方法有线性补充、多项式补充等各种方法，但只要通过预测环境时的全交叉验证达到本实施方式的用户所要求的精度，则可以使用任何补充方法。这在温度以外的环境参数的补充中也同样。然后，基于“若标高上升100m则气温下降0.6℃”这一法则，使用地形数据，将海拔0m的环境地图的各标绘点的温度的值修正为换算为原来海拔的值，由此制作原来海拔的温度的环境地图(温度地图)。
[0114]
(绝对湿度的补充)
[0115]
绝对湿度可以根据环境参数的温度和相对湿度来计算。在此，大气腐蚀环境具有“即使温度变化绝对湿度也大致恒定”的特征。因此，利用上述特征，无需反映标高等地形数据，而根据温度和相对湿度计算各网格点的绝对湿度，并对地图上的绝对湿度的数据之间进行补充，由此生成绝对湿度的环境地图(绝对湿度地图)。
[0116]
(相对湿度的补充)
[0117]
在相对湿度的补充中，利用上述的“即使温度变化绝对湿度也大致恒定”这样的大气腐蚀环境的特征，通过在温度和绝对湿度的环境地图的各网格点计算相对湿度，制作相对湿度的环境地图(相对湿度地图)。
[0118]
(润湿时间的补充)
[0119]
关于润湿时间，在国际标准中规定为“相对湿度为80％以上的时间”。因此，根据使用了上述方法的每小时、每天的相对湿度的经时变化地图，累计各网格点处的相对湿度80％以上的时间来计算年的润湿时间，由此制作润湿时间的环境地图(润湿时间地图)。
[0120]
(降雨量的补充)
[0121]
降雨量不受地形的影响，因此无需反映标高等地形数据，而是通过对地图上的降雨量的数据之间进行补充来制作降雨量的环境地图(降雨量地图)。
[0122]
(空气含盐量、so
x
浓度和no
x
浓度的补充)
[0123]
在空气含盐量、so
x
浓度和no
x
浓度的补充中，通过基于经度、纬度和海拔的欧几里得距离在数据之间进行补充，由此创建空气含盐量、so
x
浓度和no
x
浓度的环境地图(空气含盐量、so
x
浓度和no
x
浓度地图)。
[0124]
在此，空气含盐量具体通过模型公式“y＝ax
－b”计算。在该模型公式中，x：离岸距离(km)，y：空气含盐量(mdd)，a，b：系数。此外，离岸距离x是与海岸线的形状数据的各点相距的距离的最小值。需要说明的是，除了上述的模型公式以外，也可以使用cole模型、中气象模型等已有的腐蚀预测公式来计算空气含盐量。
[0125]
此外，在补充空气含盐量时，为了避免在网格点接近海岸线的情况下空气含盐量变得异常高，优选设置空气含盐量的上限。作为空气含盐量的上限，例如可以设定“1.0mdd＝62.3mmd”等值。
[0126]
(腐蚀量的预测)
[0127]
在金属材料的腐蚀量的预测中，按如上所述制作的环境地图的每个网格来预测腐蚀量。本技术发明人发现，通过分别分开预测表示各大气腐蚀环境下的金属材料的自初始起一年期间的腐蚀量的参数a(第一参数)和表示锈层导致的腐蚀速度的衰减的参数b(第二参数)，并且在预测时，根据对各环境参数加权而得的数据进行预测，腐蚀量的预测精度得以提高。
[0128]
金属材料的腐蚀速度一般随时间而衰减。这是由于在金属材料表面生成的腐蚀生成物(例如锈层)带来的保护效果。此外，该保护效果根据周围环境和金属材料的种类而有很大不同。这样，金属材料的腐蚀速度与各种环境、金属材料所具有的耐腐蚀性等非常多的要素相互关联。因此，从各环境参数与腐蚀量之间关系的原理原则来高精度地预测任意的环境及期间中的金属材料的腐蚀量是非常困难的，例如从数据库中存储的腐蚀量及各环境参数的数据组统计地预测是现实的，也关系到精度提高。
[0129]
另一方面，一般而言，在蓄积的腐蚀量和各环境参数的数据组中，长期间的数据不多。例如在预测腐蚀量的任意的期间为数十年这样的长期间的情况下，如果将期间作为变量取入到预测中而直接预测腐蚀量，则根据远离期间的数据预测任意的长期间的腐蚀量，所以精度降低。因此，在本发明中，通过分别分开预测表示金属材料的自初始起一年期间的腐蚀量的参数a和表示由锈层引起的腐蚀速度的衰减的参数b，实现精度的提高。
[0130]
在此，作为与任意期间的金属材料的腐蚀量相关的主要的环境参数，例如可以举出温度、相对湿度、绝对湿度、润湿时间、降雨量、空气含盐量、so
x
浓度及no
x
浓度等。在这些环境参数中，存在例如图1所示的温度与腐蚀量的关系那样具有非线性关系的环境参数。此外，也存在例如图2所示的空气含盐量与so2浓度的关系那样环境参数彼此显示多重共线性的情况。除了存在多个对金属材料的腐蚀造成影响的环境参数之外，这两点成为难以高精度地预测任意的环境及期间中的金属材料的腐蚀量的主要原因。
[0131]
关于环境参数相对于腐蚀量具有非线性的关系的情况，通过根据与想要预测的任意的环境以及期间的相似性对各样本进行加权，并且局部地进行多元回归分析，由此能够改善预测精度。需要说明的是，上述的“样本”是指保存在数据库中的腐蚀量和各环境参数的数据组(后述的腐蚀量数据)。
[0132]
此外，关于环境参数彼此具有多重共线性的情况，可以通过对各环境参数进行维数压缩为具有独立性的参数，生成新的参数来解决。而且，作为同时实现这些的方法之一，有以下的参考文献1所示的“局部加权偏最小二乘法(lw－pls：locally weighted partial least squares)”。
[0133]
参考文献1：金尚弘，冈岛亮太，加纳学，长谷部伸治，“用于构建高精度的局部pls模型的样本选择”，第54次自动控制联合会演讲会，54(2011)，p.1594
[0134]
在本发明中，对于表示各大气腐蚀环境下的金属材料的自初始起一年期间的腐蚀量的参数a、和表示由锈层引起的腐蚀速度的衰减的参数b，使用如下的预测方法：对每个样本求出与预测请求点的相似性(相似度)，利用该相似度进行加权，通过局部回归进行预测。此外，通过对各环境参数进行维数压缩来生成新的参数(潜在变量的导出)，作为局部回归的解释变量。此外，此时，以使用了利用相似度的加权的潜在变量与目标变量的内积成为最大的方式决定潜在变量，进行局部的多元回归。以下，参照附图对本发明的具体实施方式进行说明。
[0135]
(腐蚀量测绘装置)
[0136]
参照图3说明本发明的实施方式的金属材料的腐蚀量测绘装置的结构。腐蚀量测绘装置1具备输入部10、数据库20、运算部30和显示部40。
[0137]
输入部10例如通过键盘、鼠标指针、数字键等输入装置来实现。如后所述，经由该输入部10向运算部30输入后述的预测请求点。
[0138]
在数据库20中保存有作为金属材料的腐蚀量的实际值的腐蚀量数据。腐蚀量数据中包含金属材料(例如钢铁材料)的使用期间、该使用期间中的金属材料的腐蚀量、表示该使用期间内的金属材料的使用环境且在表示使用金属材料的地域的地图上的已知的多个环境参数、地图上的环境参数的位置坐标、地图的地形数据。
[0139]
上述的“多个环境参数”包括温度(气温)、相对湿度、绝对湿度、润湿时间及降雨量中的至少一个；空气含盐量、so
x
浓度及no
x
浓度中的至少一个。此外，这些环境参数例如是年平均的数据。此外，在数据库20中保存有每种钢种的腐蚀量数据。
[0140]
具体而言，运算部30通过由cpu(central processing unit)，dsp(digital signal processor)、fpga(field－programmable gate array)等构成的处理器和由ram(random access memory)、rom(read only memory)等构成的存储器(主存储部)来实现。
[0141]
运算部30例如将存储在未图示的存储部中的程序加载到主存储部的作业区域中并执行，通过程序的执行来控制各构成部等，由此实现与规定目的一致的功能。运算部30通过程序的执行，作为地图制作部31、相似度算出部32、维数压缩部33以及腐蚀量预测部34发挥功能。需要说明的是，关于地图制作部31、相似度算出部32、维数压缩部33以及腐蚀量预测部34的详细内容在后面叙述。
[0142]
显示部40例如由lcd显示器、crt显示器等显示装置实现，基于从运算部30输入的显示信号，显示腐蚀量预测地图作为金属材料的腐蚀量的预测结果。需要说明的是，上述的“腐蚀量预测地图”表示例如通过颜色显示腐蚀量的预测值的地图(参照图5及图6)。
[0143]
(腐蚀量测绘方法)
[0144]
参照图4对本发明的实施方式的金属材料的腐蚀量测绘方法进行说明。在腐蚀量地图方法中，进行环境地图制作步骤、预测请求点输入步骤、第一相似度算出步骤、第一维数压缩步骤、初始腐蚀量预测步骤(第一参数预测步骤)、第二相似度算出步骤、第二维数压缩步骤、衰减预测步骤(第二参数预测步骤)、长期腐蚀量预测步骤和腐蚀量预测地图制作步骤。需要说明的是，本实施方式的腐蚀量测绘方法通过应用于金属材料中尤其是耐气候性钢等钢材，能够更高精度地预测腐蚀量并进行测绘。
[0145]
在环境地图制作步骤中，由地图制作部31根据能够获得的地图上的环境参数，以任意的网格间隔对数据间进行补充，制作每个环境参数的环境地图(步骤s1)。
[0146]
在环境地图制作步骤中，根据地图上已知的多个环境参数、地图的地形数据、地图上的环境参数的位置坐标，以任意的网格间隔制作每个环境参数的环境地图。在环境地图制作步骤中，更具体而言，进行以下四个步骤。
[0147]
首先，将已知的多个环境数据输入到上述地图上的位置坐标(第一步骤)。接着，基于地图的地形数据，将环境地图修正为海拔0m的环境地图(第二步骤)。接着，在海拔0m的环境地图中，对已知的环境参数之间补充环境参数(第三步骤)。接着，基于地图的地形数据，将补充了环境参数的海拔0m的环境地图修正为原来海拔的环境地图(第四步骤)。
[0148]
在此，在上述第三步骤中，如上所述，在补充靠近海岸线的网格点的空气含盐量时，优选以不超过预先设定的上限值(例如1.0mdd)的方式补充空气含盐量。通过对空气含盐量设置上限值，能够避免接近海岸线的网格点的空气含盐量变得异常高。
[0149]
在预测请求点输入步骤中，经由输入部10向运算部30输入网格点的环境参数即预测请求点(步骤s2)。该预测请求点包括想要预测腐蚀量的金属材料的使用期间(制作腐蚀
量预测地图的金属材料的使用期间)、和表示该使用期间内的金属材料的使用环境的年平均的多个环境参数(腐蚀量数据中的多个环境参数)。
[0150]
接着，在第一相似度算出步骤中，由相似度算出部32算出保存于数据库20的使用期间为一年的金属材料的腐蚀量数据中的多个环境参数与预测请求点处的多个环境参数的相似度(步骤s3)。在本步骤中，相似度算出部32例如通过后述的式(8)算出上述相似度。需要说明的是，关于本步骤的具体例，在后述的实施例中进行说明。
[0151]
接着，在第一维数压缩步骤中，通过维数压缩部33，考虑在第一相似度算出步骤中算出的相似度，将腐蚀量数据中的多个环境参数(解释变量)维数压缩为潜在变量(步骤s4)。在本步骤中，维数压缩部33例如通过后述的式(7)来算出上述潜在变量。需要说明的是，关于本步骤的具体例，在后述的实施例中进行说明。
[0152]
接着，在初始腐蚀量预测步骤中，通过腐蚀量预测部34，基于使用在第一维数压缩步骤中算出的潜在变量以及在第一相似度算出步骤中算出的相似度而构建的预测公式，预测在预测请求点的环境参数下的金属材料的自初始起一年期间的腐蚀量(步骤s5)。该金属材料的自初始起一年期间的腐蚀量是指上述式(1)的参数a(第一参数)。在本步骤中，腐蚀量预测部34例如构建后述的式(10)所示的预测公式，基于该预测公式预测金属材料的自初始起一年期间的腐蚀量。需要说明的是，关于本步骤的具体例，在后述的实施例中进行说明。
[0153]
接着，在第二相似度算出步骤中，通过相似度算出部32，考虑想要预测腐蚀量的金属材料的使用期间，算出金属材料的腐蚀量数据中的多个环境参数与预测请求点处的多个环境参数的相似度(步骤s6)。在本步骤中，相似度算出部32例如通过后述的式(13)算出上述相似度。需要说明的是，关于本步骤的具体例，在后述的实施例中进行说明。
[0154]
接着，在第二维数压缩步骤中，通过维数压缩部33，考虑在第二相似度算出步骤中算出的相似度以及想要预测腐蚀量的金属材料的使用期间，将腐蚀量数据中的多个环境参数(解释变量)维数压缩为潜在变量(步骤s7)。在本步骤中，维数压缩部33例如通过后述的式(12)来算出上述潜在变量。需要说明的是，关于本步骤的具体例，在后述的实施例中进行说明。
[0155]
接着，在衰减预测步骤中，通过腐蚀量预测部34，基于使用在第二维数压缩步骤中算出的潜在变量及在第二相似度算出步骤中算出的相似度而构建的预测公式，预测表示金属材料的腐蚀速度的衰减的参数(步骤s8)。表示该金属材料的腐蚀速度的衰减的参数是指上述式(1)的参数b(第二参数)。在本步骤中，腐蚀量预测部34例如构建后述的式(16)所示的预测公式，基于该预测公式预测表示金属材料的腐蚀速度的衰减的参数。需要说明的是，关于本步骤的具体例，在后述的实施例中进行说明。
[0156]
接着，在长期腐蚀量预测步骤中，通过腐蚀量预测部34，基于金属材料的使用期间、在初始腐蚀量预测步骤中算出的参数a和在衰减预测步骤中算出的参数b，预测金属材料的比一年长的期间的腐蚀量(长期腐蚀量)(步骤s9)。在本步骤中，腐蚀量预测部34基于上述式(1)预测上述长期腐蚀量。需要说明的是，关于本步骤的具体例，在后述的实施例中进行说明。
[0157]
接着，在腐蚀量预测地图制作步骤中，通过地图制作部31将网格的预测请求点处的腐蚀量的预测结果在地图上涂色标注，由此制作腐蚀量预测地图(步骤s10)。该腐蚀量预
测地图是根据在腐蚀量预测地图制作步骤中预测出的腐蚀量的值，对每个所述网格进行涂色标注而成的地图(例如参照后述的图5及图6)。
[0158]
如上所述，根据使用了本实施方式的腐蚀量测绘装置1的金属材料的腐蚀量测绘方法，在大气腐蚀环境中，能够高精度地进行金属材料的长期腐蚀预测。此外，能够高精度地测绘金属材料的腐蚀量，能够选定具有与使用环境相应的耐腐蚀性的最佳的金属材料。
[0159]
此外，如果如现有技术、例如专利文献1～5那样选定作为解释变量采用的环境参数并创建公式，则会如在某个地域中预测精度较高，而在某个地域中预测精度较低那样，根据要预测腐蚀量的地域而预测精度会产生差异。另一方面，根据本实施方式的金属材料的腐蚀预测方法，由于只要在数据库20内存在与要预测腐蚀量的地域的腐蚀量数据相似的腐蚀量数据就能够预测，所以无论在怎样的地域都能够高精度地预测金属材料的腐蚀量。
[0160]
[实施例]
[0161]
(实施例1)
[0162]
以下，列举实施例更详细地说明本发明。在本实施例中，在越南，预测屋檐下环境下的钢材的年间腐蚀量(自初始起一年期间的腐蚀量)，制作腐蚀量预测地图。在此，本实施例的内容相当于上述的腐蚀量测绘方法中的环境测绘制作步骤、预测请求点输入步骤、第一相似度算出步骤、第一维数压缩步骤、始腐蚀量预测步骤及腐蚀量预测地图制作步骤。
[0163]
在本实施例中，使用保存在数据库中的世界各地的55个地区的钢材的年间腐蚀量和年平均的各环境参数的数据组(腐蚀量数据)，进行年间腐蚀量的预测。在本实施例中，作为环境参数，使用温度(℃)、相对湿度(％)、空气含盐量(mmd(cl
－
换算))、so2浓度(mmd(so2换算))这四个。在此，上述的“mmd”是每单位天数和单位面积所捕集的cl
－
或so2的量，是“mg
·m－2
·
day
－1”的缩写。
[0164]
每个环境参数的环境地图通过上述方法制作。在本次的实施例中，作为数据间的补充方法，利用线性补充进行了计算(图4的“环境地图制作步骤”)。在得到的环境地图的网格点处的腐蚀预测中，如上所述，在解释变量彼此具有相关性的情况下，已知由于多重共线性而预测精度变差。在上述的lw－pls中，由于作为解释变量使用的各环境参数具有相关性，所以通过维数压缩排除相关性，生成新的参数(潜在变量)。在此，lw－pls按照上述参考文献1所示的步骤进行计算。用于将四个环境参数压缩为潜在变量(参数t)的式子可以如下式(7)所示。
[0165]
[数学式7]
[0166]
t＝w1t w2rh w3cl w4so2
…
(7)
[0167]
在此，上述式(7)中，t为温度(℃)，rh为相对湿度(％)，cl为空气含盐量(mg/m2/day(＝mmd)(cl
－
换算))，so2为so2浓度(mg/m2/day(＝mmd)(so2换算))，w1～w4为系数。需要说明的是，在上述式(7)中，作为一个例子仅使用了四个环境参数，但实际上，优选包含所有在要预测腐蚀量的地域中预想与腐蚀有关的环境参数。
[0168]
在本实施例中，首先将预测请求点输入到腐蚀量测绘装置的运算部(图4的“预测点输入步骤”)，该预测请求点包含想要预测腐蚀量的钢材的使用期间和表示该使用期间内的钢材的使用环境的年平均的多个环境参数。
[0169]
接着，使用下式(8)所示的欧几里得距离算出想要求出腐蚀量的预测请求点的环境参数与为了预测腐蚀量而参照的腐蚀量数据的环境参数i的相似度ωi(图4的“第一相似
度算出步骤”)。需要说明的是，由于数据的范围因各环境参数而不同，因此对数值实施了标准化。
[0170]
[数学式8]
[0171][0172]
在此，ωi为局部化参数，右下的下标为q的环境参数为想要预测腐蚀量的预测请求点的环境参数，右下的下标为i的环境参数为从数据库参照的腐蚀量数据的环境参数，σ为下述式(9)的标准偏差。此外，φ是调整参数，例如以φ＝1为基准适当调整来决定值。
[0173]
[数学式9]
[0174][0175]
接着，根据由上述式(8)算出的相似度ωi和腐蚀量数据的环境参数及腐蚀量，按照参考文献1(2.1章)所示的步骤，以相似度ωi为权重，以关联的潜在变量和目标变量(腐蚀量)的内积为最大的方式，决定上述式(7)的系数w1～w4。然后，使用所决定的系数w1～w4，通过上述式(7)算出各环境参数的潜在变量(图4的“第一维数压缩步骤”)。
[0176]
接着，通过局部回归构建下述式(10)所示的腐蚀量的预测公式，基于下述式(10)，对预测请求点的环境参数下的钢材的年间腐蚀量(自初始起一年期间的腐蚀量)进行预测(图4的“初始腐蚀量预测步骤”)。
[0177]
[数学式10]
[0178]
y＝αt
…
(10)
[0179]
在此，y是腐蚀量的预测值，α是系数(回归系数)。需要说明的是，在上述式(10)中省略了记载，但也可以在上述式(10)中包含常数项，或者使用多个潜在变量。
[0180]
这样，在本实施例中，每当想要预测腐蚀量的预测请求点被输入时，进行预测请求点与各腐蚀量数据的相似度的算出、潜在变量的系数的算出、预测公式的构建。将得到的每个网格点的年间腐蚀量的预测值做成地图而作为例如图5所示的腐蚀量预测地图(图4的“腐蚀量预测地图制作步骤”)，并通过显示部40进行显示。
[0181]
(实施例2)
[0182]
在本实施例中，在越南，预测在屋檐下环境的钢材的长期(五十年)腐蚀量，制作腐蚀量预测地图。在此，本实施例的内容相当于上述的腐蚀量测绘方法的全部步骤。
[0183]
在本实施例中，以屋檐下环境为对象，使用保存在数据库中的世界55个地区的1年间、55个地区的3年间、39个地区的5年间、38个地区的7年间、38个地区的9年间的腐蚀量和年平均的各环境参数的数据组(腐蚀量数据)，以与实施例1同样的方法，对钢材算出了表示自初始起一年期间的腐蚀量的参数a(参照上述式(1)参照)(图4的环境地图制作步骤、预测请求点输入步骤、第一相似度算出步骤、第一维数压缩步骤和初始腐蚀量预测步骤)。
[0184]
接着，对7年间的长期腐蚀试验结果进行加权，算出表示由锈层引起的腐蚀速度的衰减的参数b，通过上述式(1)预测钢材的五十年后的腐蚀量。在五十年后的腐蚀量的预测中，首先将上述式(1)变形为下述式(11)所示，使用参数b及使用期间x算出左边的自初始起
一年期间的腐蚀量的对数与五十年后的腐蚀量的对数之差。
[0185]
[数学式11]
[0186]
logy-loga＝blogx
…
(11)
[0187]
在上述式(11)的左边的预测中使用lw－pls。具体而言，通过取四个环境参数与使用期间x的对数的积，制作考虑了时间因子的新的环境参数，将这些新的环境参数如下述式(12)所示压缩为潜在变量(参数u)。
[0188]
[数学式12]
[0189]
u＝(v1t v2rh v3cl v4so2)logx
…
(12)
[0190]
在此，上述式(12)中，t为温度(℃)，rh为相对湿度(％)，cl为空气含盐量(mg/m2/day(＝mmd)(cl
－
换算))，so2为so2浓度(mg/m2/day(＝mmd)(so2换算))，v1～v4为系数。
[0191]
接着，使用下式(13)所示的欧几里得距离算出想要求出腐蚀量的预测请求点的环境参数与为了预测腐蚀量而参照的腐蚀量数据的环境参数i的相似度ωi(图4的“第二相似度算出步骤”)。需要说明的是，由于数据的范围因各环境参数而不同，因此对数值实施了标准化。
[0192]
[数学式13]
[0193][0194]
在此，ωi为局部化参数，右下的下标为q的环境参数为想要预测腐蚀量的预测请求点的环境参数，右下的下标为i的环境参数为从数据库参照的腐蚀量数据的环境参数，σ为下述式(14)的标准偏差。此外，φ是调整参数，例如以φ＝1为基准适当调整来决定值。
[0195]
[数学式14]
[0196][0197]
接着，根据由上述式(13)算出的相似度ωi和腐蚀量数据的环境参数及腐蚀量，按照参考文献1(2.1章)所示的步骤，以相似度ωi为权重，以关联的潜在变量和目标变量(腐蚀量)的内积为最大的方式，决定上述式(12)的系数v1～v4。然后，使用所决定的系数v1～v4，通过上述式(7)计算出各环境参数的潜在变量(图4的“第二维数压缩步骤”)。
[0198]
接着，通过局部回归构建下述式(15)所示的腐蚀量的预测公式，算出想要预测的环境参数中的腐蚀量的对数之差。
[0199]
[数学式15]
[0200]
logy-loga＝βu
…
(15)
[0201]
在此，在上述式(15)中，β是系数(回归系数)。需要说明的是，在上述式(15)中省略了记载，但也可以在上述式(15)中包含常数项，或者使用多个潜在变量。
[0202]
这样，在本发明例中，每当想要预测腐蚀量的预测请求点被输入时，进行预测请求
点与各腐蚀量数据的相似度的算出、潜在变量的系数的算出和预测公式的构建。
[0203]
接着，基于上述式(12)和上述式(15)，通过下述式(16)计算出上述式(11)的参数b(图4的“衰减预测步骤”)。
[0204]
[数学式16]
[0205][0206]
接着，使用算出的参数a和参数b，通过上述式(1)计算出9年后的腐蚀量y(图4的长期腐蚀量预测步骤)。将得到的每个网格点的长期腐蚀量的预测值做成地图而作为例如图6所示的腐蚀量预测地图(图4的“腐蚀量预测地图制作步骤”)，并通过显示部40进行显示。
[0207]
以上，对于本发明的金属材料的腐蚀量测绘方法、金属材料的选定方法及金属材料的腐蚀量测绘装置，通过用于实施发明的方式及实施例具体地进行了说明，但本发明的主旨并不限定于上述记载，应基于权利要求书的记载而广泛地解释。此外，基于上述记载进行了各种变更、改变等的方式当然也包含在本发明的主旨中。
[0208]
在此，在上述的实施方式中，分别分开预测金属材料的自初始起一年期间的腐蚀量(参数a)和表示金属材料的腐蚀速度的衰减的参数(参数b)，以自初始起一年期间的腐蚀量为基准预测长期腐蚀量，但预测长期腐蚀量时的基准并不限定于自初始起一年期间的腐蚀量。
[0209]
即，也可以在初始腐蚀量预测步骤中，预测预定的任意的规定期间(初始期间)中的金属材料的腐蚀量，在长期腐蚀量预测步骤中，以所述的规定期间的腐蚀量为基准预测长期腐蚀量。
[0210]
例如，作为初始腐蚀量，1.5年的腐蚀量作为a'给出的情况下，认为由此开始x年后的腐蚀量的预测公式可以扩展上述式(1)，如下述式(17)那样记述。
[0211]
[数学式17]
[0212][0213]
如果将其一般化，则将某初始期间x0年的腐蚀量设为a'，将以x0年为基准的衰减参数设为b'，能够得到下述式(18)。通过使用该式(18)，作为以x0年期间为基准的腐蚀量，可以算出x》x0的期间的腐蚀量。
[0214]
[数学式18]
[0215]
y＝a
′
x
′b＝a
′
(x/x0)b′
…
(18)
[0216]
分别分开预测金属材料的任意初始期间的腐蚀量a'和衰减参数b'，如上述式(18)所示，通过用衰减参数b'对初始期间以后的经过年数x'求幂，可以预测初始期间以后的长期腐蚀量。但是，上述式(1)的初始腐蚀量a以1年期间的腐蚀量为前提。因此，对于上述式(18)的x0的期间，不设想大幅偏离1年的情况，认为从半年到2年左右是现实的实用范围。
[0217]
[附图标记的说明]
[0218]1ꢀꢀ
腐蚀量测绘装置
[0219]
10 输入部
[0220]
20 数据库
[0221]
30 运算部
[0222]
31 地图制作部
[0223]
32 相似度算出部
[0224]
33 维数压缩部
[0225]
34 腐蚀量预测部
[0226]
40 显示部