(57)【要約】【課題】リスク評価の対象である評価対象物に実際に損傷が発生する前であっても、実際の評価対象物により即した損傷確率に基づいてリスク評価を可能とする。【解決手段】特定の評価対象物の損傷確率を算出する損傷確率算出装置であって、損傷率と時間との関係を示す損傷率関数を記憶する記憶部１と、上記評価対象物の通常検査時に取得された検査データに基づいて上記損傷率関数を補正する損傷率関数補正部３と、該損傷率関数補正部３にて補正された上記損傷率関数に基づいて上記損傷確率を算出する損傷確率算出部４とを備える。【選択図】図１
2009288854.tif 000002
【特許請求の範囲】
【請求項１】
特定の評価対象物の損傷確率を算出する損傷確率算出装置であって、
損傷率と時間との関係を示す損傷率関数を記憶する記憶部と、
前記評価対象物の通常検査時に取得された検査データに基づいて前記損傷率関数を補正する損傷率関数補正部と、
該損傷率関数補正部にて補正された前記損傷率関数に基づいて前記損傷確率を算出する損傷確率算出部と
を備えることを特徴とする損傷確率算出装置。
【請求項２】
前記損傷率関数補正部は、前記損傷率関数が含む係数に対する尤度を前記検査データに基づいて算出し、ベイズの定理に基づいて前記尤度にて前記係数を補正することによって前記損傷率関数を補正することを特徴とする請求項１記載の損傷確率算出装置。
【請求項３】
前記損傷率関数補正部は、前記損傷率関数が複数の係数を含む場合に、値の変更が最も前記損傷率関数に対して大きな影響を与える前記係数のみを前記尤度にて補正することを特徴とする請求項２記載の損傷確率算出装置。
【請求項４】
過去の損傷データから前記損傷率関数を算出する損傷率関数算出部をさらに備えることを特徴する請求項１〜３いずれかに記載の損傷確率算出装置。
【請求項５】
特定の評価対象物の損傷確率を算出する損傷確率算出方法であって、
損傷率と時間との関係を示す損傷率関数を算出する損傷率関数算出工程と、
前記評価対象物の通常検査時に取得された検査データに基づいて前記損傷率関数を補正する損傷率関数補正工程と、
該損傷率関数補正工程にて補正された前記損傷率関数に基づいて前記損傷確率を算出する損傷確率算出工程と
を備えることを特徴とする損傷確率算出方法。
【請求項６】
前記損傷率関数補正工程にて、前記損傷率関数が含む係数に対する尤度を前記検査データに基づいて算出し、ベイズの定理に基づいて前記尤度にて前記係数を補正することによって前記損傷率関数を補正することを特徴とする請求項５記載の損傷確率算出方法。
【請求項７】
前記損傷率関数が複数の係数を含む場合に、前記損傷率関数補正工程にて、値の変更が最も前記損傷率関数に対して大きな影響を与える前記係数のみを前記尤度にて補正することを特徴とする請求項６記載の損傷確率算出方法。
【請求項８】
特定の評価対象物に対する損傷確率と、該特定の評価対象物が損傷することによる被害の影響度とに基づいて前記特定の評価対象物に対するリスクを評価し、該リスクに基づいて前記特定の評価対象物に対するメンテナンスを行うメンテナンス方法であって、
請求項５〜７いずれかに記載の損傷確率算出方法にて算出された損傷確率を用いて前記リスクを評価することを特徴とするメンテナンス方法。
【請求項９】
コンピュータを特定の評価対象物の損傷確率を算出する手段として機能させる損傷確率算出プログラムであって、
前記コンピュータを、
損傷率と時間との関係を示す損傷率関数を記憶する記憶手段と、
前記評価対象物の通常検査時に取得された検査データに基づいて前記損傷率関数を補正する損傷率関数補正手段と、
該損傷率関数補正手段にて補正された前記損傷率関数に基づいて前記損傷確率を算出する損傷確率算出手段と
して機能させることを特徴とする損傷確率算出プログラム。
【請求項１０】
前記コンピュータを、前記損傷率関数補正手段として、前記損傷率関数が含む係数に対する尤度を前記検査データに基づいて算出し、ベイズの定理に基づいて前記尤度にて前記係数を補正するように機能させることを特徴とする損傷確率算出プログラム。
【請求項１１】
前記コンピュータを、前記損傷率関数補正手段として、前記損傷率関数が複数の係数を含む場合に、値の変更が最も前記損傷率関数に対して大きな影響を与える前記係数のみを前記尤度にて補正するように機能させることを特徴とする請求項１０記載の損傷確率算出プログラム。
【請求項１２】
前記コンピュータを、過去の損傷データから前記損傷率関数を算出する損傷率関数算出手段として機能させることを特徴とする請求項９〜１１いずれかに記載の損傷確率算出プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、特定の評価対象物の損傷確率を算出する損傷確率算出装置及び損傷確率算出方法、該損傷確率算出方法にて算出された損傷確率を用いて特定の評価対象物のメンテナンスを行うメンテナンス方法、及びコンピュータを特定の評価対象物の損傷確率を算出する手段として機能させる損傷確率算出プログラムに関するものである。
【背景技術】
【０００２】
近年、合理的かつ経済的な保守管理として、複数の評価対象物をリスク評価してリスクの高い評価対象物の検査を優先的に実施するＲＢＭ（Risk-Based Maintenance）による保守管理が行われている。
このＲＢＭによる保守管理では、ある事象（損傷）が発生する確率とその影響度との組み合わせによって評価対象物のリスク評価を行っている。より詳細には、保守管理の対象となる評価対象物を、縦軸が損傷確率、横軸が被害の影響度とされたリスクマトリックスを用いて分類することによって、リスク評価が行われる。
【０００３】
ところで、リスク評価の際に用いられる損傷確率は、損傷率と時間との関係を示す関係式（損傷率関数）に基づいて算出される。そして、通常損傷率は、過去の複数の損傷データに基づいて設定される。つまり、リスク評価の際に用いられる損傷確率は、過去の損傷データに基づいて算出される。しかしながら、リスク評価の対象とされている評価対象物が置かれている環境は、必ずしも過去の評価対象物と完全に一致するわけではない。
【０００４】
このため、特許文献１には、リスク評価の対象とされている評価対象物に故障（損傷）が生じた場合に、この故障の発生を踏まえて故障確率（損傷確率）を改めて算出する方法が提案されている。
このような特許文献１に示される方法によれば、実際のリスク評価の対象となる評価対象物に生じた故障の発生に基づいて故障確率が改めて算出されるため、実際のリスク評価の対象となる評価対象物により即した故障確率を用いてリスク評価が可能となる。
【特許文献１】特開２００７−３２８５２２号公報
【発明の開示】
【発明が解決しようとする課題】
【０００５】
しかしながら、特許文献１の方法では、リスク評価の対象とされた評価対象物に実際に損傷が生じない限りは、新たな損傷確率を算出することができない。
実際の現場では、より合理的かつ経済的な保守管理を実現するために、実際に損傷が生じるより前に、実際のリスク評価の対象となる評価対象物により即した損傷確率を用いてリスク評価を行うことが望まれている。
【０００６】
本発明は、上述する問題点に鑑みてなされたもので、リスク評価の対象である評価対象物に実際に損傷が発生する前であっても、実際の評価対象物により即した損傷確率に基づいてリスク評価を可能とすることを目的とする。
【課題を解決するための手段】
【０００７】
本発明は、上記課題を解決するための手段として、以下の構成を採用する。
【０００８】
第１の発明は、特定の評価対象物の損傷確率を算出する損傷確率算出装置であって、損傷率と時間との関係を示す損傷率関数を記憶する記憶部と、上記評価対象物の通常検査時に取得された検査データに基づいて上記損傷率関数を補正する損傷率関数補正部と、該損傷率関数補正部にて補正された上記損傷率関数に基づいて上記損傷確率を算出する損傷確率算出部とを備えるという構成を採用する。
【０００９】
このような構成を採用する第１の発明によれば、評価対象物の通常検査時に取得された検査データに基づいて当該評価対象物に関する損傷率関数が補正され、補正された損傷率関数に基づいて損傷確率が算出される。
【００１０】
第２の発明は、上記第１の発明において、上記損傷率関数補正部が、上記損傷率関数が含む係数に対する尤度を上記検査データに基づいて算出し、ベイズの定理に基づいて上記尤度にて上記係数を補正することによって上記損傷率関数を補正するという構成を採用する。
【００１１】
第３の発明は、上記第２の発明において、上記損傷率関数補正部が、上記損傷率関数が複数の係数を含む場合に、値の変更が最も上記損傷率関数に対して大きな影響を与える上記係数のみを上記尤度にて補正するという構成を採用する。
【００１２】
第４の発明は、上記第１〜第３いずれかの発明において、過去の損傷データから上記損傷率関数を算出する損傷率関数算出部をさらに備えるという構成を採用する。
【００１３】
第５の発明は、特定の評価対象物の損傷確率を算出する損傷確率算出方法であって、損傷率と時間との関係を示す損傷率関数を算出する損傷率関数算出工程と、上記評価対象物の通常検査時に取得された検査データに基づいて上記損傷率関数を補正する損傷率関数補正工程と、該損傷率関数補正工程にて補正された上記損傷率関数に基づいて上記損傷確率を算出する損傷確率算出工程とを備えるという構成を採用する。
【００１４】
このような構成を採用する第５の発明によれば、評価対象物の通常検査時に取得された検査データに基づいて当該評価対象物に関する損傷率関数が補正され、補正された損傷率関数に基づいて損傷確率が算出される。
【００１５】
第６の発明は、上記第５の発明において、上記損傷率関数補正工程にて、上記損傷率関数が含む係数に対する尤度を上記検査データに基づいて算出し、ベイズの定理に基づいて上記尤度にて上記係数を補正することによって上記損傷率関数を補正するという構成を採用する。
【００１６】
第７の発明は、上記第６の発明において、上記損傷率関数が複数の係数を含む場合に、上記損傷率関数補正工程にて、値の変更が最も上記損傷率関数に対して大きな影響を与える上記係数のみを上記尤度にて補正するという構成を採用する。
【００１７】
第８の発明は、特定の評価対象物に対する損傷確率と、該特定の評価対象物が損傷することによる被害の影響度とに基づいて上記特定の評価対象物に対するリスクを評価し、該リスクに基づいて上記特定の評価対象物に対するメンテナンスを行うメンテナンス方法であって、上記第５〜第７いずれかの発明にて算出された損傷確率を用いて上記リスクを評価するという構成を採用する。
【００１８】
第９の発明は、コンピュータを特定の評価対象物の損傷確率を算出する手段として機能させる損傷確率算出プログラムであって、上記コンピュータを、損傷率と時間との関係を示す損傷率関数を記憶する記憶手段と、上記評価対象物の通常検査時に取得された検査データに基づいて上記損傷率関数を補正する損傷率関数補正手段と、該損傷率関数補正手段にて補正された上記損傷率関数に基づいて上記損傷確率を算出する損傷確率算出手段として機能させるという構成を採用する。
【００１９】
このような構成を採用する第９の発明によれば、評価対象物の通常検査時に取得された検査データに基づいて当該評価対象物に関する損傷率関数が補正され、補正された損傷率関数に基づいて損傷確率が算出される。
【００２０】
第１０の発明は、上記第９の発明において、上記コンピュータを、上記損傷率関数補正手段として、上記損傷率関数が含む係数に対する尤度を上記検査データに基づいて算出し、ベイズの定理に基づいて上記尤度にて上記係数を補正するように機能させるという構成を採用する。
【００２１】
第１１の発明は、上記第１０の発明において、上記コンピュータを、上記損傷率関数補正手段として、上記損傷率関数が複数の係数を含む場合に、値の変更が最も上記損傷率関数に対して大きな影響を与える上記係数のみを上記尤度にて補正するように機能させるという構成を採用する。
【００２２】
第１２の発明は、上記第９〜第１１のいずれかの発明において、上記コンピュータを、過去の損傷データから上記損傷率関数を算出する損傷率関数算出手段として機能させるという構成を採用する。
【発明の効果】
【００２３】
本発明によれば、評価対象物の通常検査時に取得された検査データに基づいて当該評価対象物に関する損傷率関数が補正され、補正された損傷率関数に基づいて損傷確率が算出される。
このため、通常検査にて評価対象物に損傷が発生していない場合であっても、当該評価対象物の損傷率関数が補正され、この補正された損傷率関数に基づいて損傷確率が算出される。
したがって、本発明によれば、リスク評価の対象である評価対象物に実際に損傷が発生する前であっても、実際の評価対象物により即した損傷確率に基づいてリスク評価を行うことが可能となる。
【発明を実施するための最良の形態】
【００２４】
以下、図面を参照して、本発明に係る損傷確率算出装置、損傷確率算出方法、メンテナンス方法及び損傷確率算出プログラムの一実施形態について説明する。
なお、以下の説明においては、本発明における評価対象物が、船の主機関が備えるシリンダカバーである例を挙げて説明する。
【００２５】
図１は、本実施形態の損傷確率算出装置Ｓ１の機能構成を示すブロック図である。この図に示すように、本実施形態の損傷確率算出装置Ｓ１は、記憶部１と、損傷率関数算出部２と、損傷率関数補正部３と、損傷確率算出部４と、入力装置５と、表示装置６とを備えている。なお、実際には、記憶部１は、メモリやハードディスク等によって構成されている。また、損傷率関数算出部２、損傷率関数補正部３及び損傷確率算出部４は、ＣＰＵ等の演算処理装置によって構成されている。
【００２６】
記憶部１は、本実施形態の損傷確率算出装置Ｓ１の動作を制御するためのプログラムや各種データを記憶するものである。
そして、本実施形態の損傷確率算出装置Ｓ１において記憶部１は、過去のシリンダカバーの損傷データを記憶している。なお、ここで言う、過去のシリンダカバーとは、本実施形態の損傷確率算出装置Ｓ１が損傷確率を算出する対象としているシリンダカバーに限られるものではなく、他の船に設置された過去のシリンダカバー等の損傷データを含むものである。
また、記憶部１は、損傷確率を算出するために用いる損傷率関数（損傷率と時間との関係式）及び演算式を記憶する。
【００２７】
なお、本実施形態においては、損傷率関数が損傷率関数算出部２にて算出されるため、記憶部１には損傷率関数算出部２にて算出された損傷率関数が記憶されるが、損傷率関数が損傷確率算出装置Ｓ１の外部にて算出される場合には、記憶部１には予め損傷率関数が記憶されている。
さらに、記憶部１に記憶される過去のシリンダカバーは、損傷率関数の算出のために用いられるデータである。このため、記憶部１に予め損傷率関数が記憶されている場合には、記憶部１が過去のシリンダカバーの損傷データを記憶する必要はない。
【００２８】
また、記憶部１は、本実施形態の損傷確率算出装置Ｓ１によって損傷確率を算出するシリンダカバーに対して行われる通常の検査時に取得された検査データが記憶される。
なお、ここで言う、通常の検査とは、予め定められた期間の経過後に行われる検査であり、シリンダカバーが損傷していない場合であっても行われる検査である。
また、記憶部１に記憶される検査データは、入力装置５を介して記憶部１に入力され、検査の時期、損傷の有無及び損傷の具合を含むデータである。
【００２９】
損傷率関数算出部２は、記憶部１に記憶された過去のシリンダカバーの損傷データに基づいて、損傷率関数を算出する。
具体的には、記憶部１に記憶された過去のシリンダカバーの損傷データには、損傷の発生タイミングが含まれるため、例えば、損傷率関数算出部２は、時間ごとの損傷の発生率（損傷率）を算出し、この発生率を時系列的に並べて、さらに損傷率曲線を求める。そして、損傷率曲線を数式として表し、この数式を損傷率関数とする。このような損傷率関数は、損傷率と時間との関係を示す関係式であり、確率変数である係数を少なくとも一つ含む関数として表される。
【００３０】
損傷率関数補正部３は、記憶部１に記憶された上記検査データに基づいて、同じく記憶部１に記憶された損傷率関数を補正する。
具体的には、損傷率関数補正部３は、損傷率関数が含む係数に対する尤度を検査データに基づいて算出し、ベイズの定理を用いて算出した尤度にて上記係数を補正することによって損傷率関数を補正する。
なお、ベイズの定理は、得られたある結果（検査データ）に基づいて事前確率（補正前の損傷率関数）から事後確率（補正後の損傷率関数）を求める定理である。なお、このベイズの定理は、周知の定理であるため、ここでの詳細な説明は省略する。
【００３１】
また、本実施形態の損傷確率算出装置Ｓ１において損傷率関数補正部３は、損傷率関数が複数の係数を含む場合に、値の変更が最も損傷率関数に対して大きな影響を与える係数（以下、最大寄与係数と称する）のみに対する尤度を算出し、この算出した尤度に基づいて最大寄与係数のみを補正する。
なお、本実施形態の損傷確率算出装置Ｓ１においては、計算作業の効率化を図るために、最大寄与係数のみを補正する構成を採用するが、迅速な計算が可能な環境であれば全ての係数を補正するようにしても良い。また、複数の係数うち、値の変更が損傷率関数に対して大きな影響を与える順に、予め決められた数の係数を補正するようにしても良い。
【００３２】
損傷確率算出部４は、損傷率関数（補正前の損傷率関数と補正後の損傷率関数の両方を含む）に基づいて、記憶部１に記憶された演算式を用いてシリンダカバーの損傷確率を算出するものである。
この損傷確率算出部４は、損傷率関数補正部３にて損傷率関数が補正される前には、補正前の損傷率関数に基づいてシリンダカバーの損傷確率を算出し、損傷率関数補正部３にて損傷率関数が補正された後には、補正後の損傷率関数に基づいてシリンダカバーの損傷確率を算出する。
【００３３】
入力装置５は、本実施形態の損傷確率算出装置Ｓ１に、必要なデータを入力する際等に、作業者によって操作されるものであり、キーボードやマウスによって構成されている。
そして、入力装置５を介して、シリンダカバーの検査時に取得された検査データ等が入力される。また、入力装置５は、上記検査データの他にも、本実施形態の損傷確率算出装置Ｓ１において必要とんされる各種データの入力に用いられる。
【００３４】
表示装置６は、損傷確率算出部４にて算出された損傷確率を視覚化して表示するものであり、ディスプレイやプリンタによって構成されている。
なお、表示装置６は、損傷確率の他にも、本実施形態の損傷確率算出装置Ｓ１において視覚化が必要とされる各種データを表示する。
【００３５】
次に、上述のように構成された本実施形態の損傷確率算出装置Ｓ１の動作（損傷確率算出方法）について、図２のフローチャートを参照して説明する。
【００３６】
図２に示すように、本実施形態の損傷確率算出方法は、損傷率関数算出工程（ステップＳ１）と、損傷率関数補正工程（ステップＳ２）と、損傷確率算出工程（ステップＳ３）とを有している。
【００３７】
損傷率関数算出工程（ステップＳ１）は、損傷率と時間との関係を示す上記損傷率関数を算出工程である。本損傷率関数算出工程（ステップＳ１）においては、損傷率関数算出部２が、記憶部１に記憶された過去のシリンダカバーの損傷データに基づいて、損傷率関数を算出する。
なお、実際に過去のシリンダカバーの損傷データに基づいて損傷率関数を算出すると、損傷率関数λ（損傷率（％）／時間（年））は、時間Ｔの関数として、下式（１）のように指数関数として表される。
したがって、本実施形態の損傷確率算出方法においても、式（１）が損傷率関数として算出されたものとする。
そして、本損傷率関数算出工程（ステップＳ１）にて算出された損傷率関数λは、記憶部１に記憶される。
【００３８】
【数１】
2009288854.tif 000003
【００３９】
続いて、損傷率関数補正工程（ステップＳ２）は、損傷確率の算出対象であるシリンダカバーの通常検査時に取得された検査データに基づいて、ステップＳ１にて算出された損傷率関数を補正する工程である。本損傷率関数補正工程（ステップＳ２）においては、損傷率関数補正部３が、記憶部１に記憶された上記検査データに基づいて、同じく記憶部１に記憶された損傷率関数を補正する。
【００４０】
ここで、本実施形態の損傷確率算出方法においては、式（１）に示される損傷率関数λがａ及びｂの２つの係数を含んでいる。これら２つの係数は、記憶部１に記憶された損傷データの分散を示す確率変数である。そして、この２つ係数のうち、値の変化が損傷率関数λに最も大きな影響を与える係数ａを最大寄与係数として確率変数として扱い、値の変化が損傷率関数λに与える影響の小さい係数ｂを固定値として扱う。
なお、係数ａは、ａ＞０の場合に、係数ａの平均値をａｍ、標準偏差をσａの正規分布をした確率変数とする。
【００４１】
なお、損傷率関数λ（Ｔ）とＴ年までに損傷する確率である損傷確率Ｆ（Ｔ）との関係は、下式（２）によって表すことができる。また、下式（２）から損傷確率Ｆ（Ｔ）は、下式（３）によって表すことができる。
【００４２】
【数２】
2009288854.tif 000004
【００４３】
【数３】
2009288854.tif 000005
【００４４】
そして、シリンダカバーは初期状態で健全性が確認されているため、Ｔ＞０においては、下式（４）となる。
【００４５】
【数４】
2009288854.tif 000006
【００４６】
この式（４）から、例えば、５年まで損傷が発見されない確率は下式（５）となり、１０年まで損傷が発見されず、１５年目で損傷が発見される確率は下式（６）となる。
【００４７】
【数５】
2009288854.tif 000007
【００４８】
【数６】
2009288854.tif 000008
【００４９】
よって、確率変数であるａを式（５）や式（６）を尤度［ｌ（ａ）］として扱うことができ、ベイズの定理により係数ａを補正することが可能となる。
【００５０】
ここで、ｉ回目の検査をＴｉ年目に通常検査を行った場合を考える。なお、ａの平均値をａｍ＿ｉ−１、標準偏差をσａ＿ｉ−１と表し、ａｍ＿０及びσａ＿０が初期値とする。
そして、ｉ−１回目の通常検査（Ｔｉ−１年目）まで損傷が発見されず、ｉ回目の通常検査（Ｔｉ年目）でも損傷が発見されなかった場合、ｉ回目の通常検査での尤度ｌｉ（ａ）は、下式（７）となる。
【００５１】
【数７】
2009288854.tif 000009
【００５２】
この場合、ｉ回目までの尤度は、下式（８）となる。
【００５３】
【数８】
2009288854.tif 000010
【００５４】
一方、ｉ−１界面の通常検査（Ｔｉ−１年目）まで損傷が発見されず、ｉ回目の通常検査（Ｔｉ年目）で損傷が発見された場合、ｉ回目の通常検査の尤度ｌｉ（ａ）は、下式（９）となる。
【００５５】
【数９】
2009288854.tif 000011
【００５６】
なお、ｉ−１回目までの尤度は、上式（４）で表されるため、ｉ回目までの尤度は、下式（１０）となる。
【００５７】
【数１０】
2009288854.tif 000012
【００５８】
このようにして、通常検査の検査データに基づいて、検査結果に応じた尤度を算出する。そして、算出された尤度に基づいて、係数ａが補正され、さらには損傷率関数が補正される。
【００５９】
本損傷率関数補正工程（ステップＳ２）において行われる損傷率関数の補正を図３に示すグラフを用いて説明する。なお、図３において、横軸は時間、縦軸は損傷率を示している。
この図３に示すように、係数ａが確率変数であり分布を有しており、補正前の損傷率関数は、その平均曲線が通常検査のタイミングＴｍにおいて補正前の係数ａの平均値を通過するグラフＡとして表される。そして、本実施形態の損傷確率算出方法においては、ベイズの定理を用いて、通常検査のタイミングＴｍにおいて算出した尤度ｌ（ａ）にて係数ａを補正する。この結果、補正後の係数ａである係数ａ’は、その平均値が補正前の係数ａに対して移動する。このため、損傷率関数は、その平均曲線が通常検査のタイミングＴｍにおいて補正後の係数ａ（係数ａ’）の平均値を通過するグラフＢとして表される。そして、このグラフＢが、補正された損傷率関数となる。
【００６０】
続いて、損傷確率算出工程（ステップＳ３）は、損傷率関数補正工程（ステップＳ２）にて補正された損傷率関数に基づいて損傷確率を算出する工程である。本損傷確率算出工程（ステップＳ３）においては、損傷確率算出部４が、損傷率関数に基づいて、記憶部１に記憶された演算式を用いてシリンダカバーの損傷確率を算出する。そして、損傷確率算出部４は、上式（２）及び上式（３）に基づいてシリンダカバーの損傷確率を算出する。
【００６１】
なお、損傷確率を算出する場合には、損傷確率を算出する期間までに、シリンダカバーに対してどのようなメンテナンスが行われたかを考慮する必要がある。つまり、損傷率関数の時間軸と、現実の時間軸とをずらす必要が生じる。
【００６２】
より詳細には、式（２）及び式（３）から分かるように、ある期間におけるシリンダカバーの損傷確率を算出するためには、ある期間に相当する期間で損傷率関数を積分することによって求められる。ここで、現実の時間におけるＴｉ年目からＴｊ年目の損傷確率を算出する場合について考える。
シリンダカバーの交換が行われていない場合には、損傷率関数の時間軸と現実の時間軸とが一致するため、上記積分の開始点（Ｔｍａｉｎｔ）は、下式（１１）に示すように、損傷率関数におけるＴｉとなる。このため、損傷確率の算出対象であるシリンダカバーの累積使用時間（Ｔｔｏｔａｌ）、つまり積分の終了点は、下式（１２）に示すように、Ｔｊとなる。
【００６３】
【数１１】
2009288854.tif 000013
【００６４】
【数１２】
2009288854.tif 000014
【００６５】
また、損傷確率を算出する期間までの間にシリンダカバーが既に使用されたことのあるシリンダカバーに置き換えられた場合には、損傷率関数の時間軸と現実の時間軸とが一致せず、上記積分の開始点（Ｔｍａｉｎｔ）は、下式（１３）に示すように、損傷率関数におけるＴｋとなる。このため、損傷確率の算出対象であるシリンダカバーの累積使用時間（Ｔｔｏｔａｌ）、つまり積分の終了点は、下式（１４）に示すように、Ｔｊ−Ｔｉ＋Ｔｋとなる。
【００６６】
【数１３】
2009288854.tif 000015
【００６７】
【数１４】
2009288854.tif 000016
【００６８】
また、損傷確率を算出する期間までの間にシリンダカバーが新品のシリンダカバーに交換された場合には、損傷率関数の時間軸と現実の時間軸とが一致せず、上記積分の開始点（Ｔｍａｉｎｔ）は、下式（１５）に示すように、損傷率関数におけるＴ０（始点）となる。このため、損傷確率の算出対象であるシリンダカバーの累積使用時間（Ｔｔｏｔａｌ）、つまり積分の終了点は、下式（１６）に示すように、Ｔｊ−Ｔｉとなる。
【００６９】
【数１５】
2009288854.tif 000017
【００７０】
【数１６】
2009288854.tif 000018
【００７１】
このように損傷確率を算出する期間までの間にシリンダカバーが置き換えあるいは交換された場合には、置き換えあるいは交換の時期及び置き換えられるシリンダカバーの使用済み期間に応じて、損傷率関数の時間軸と現実の時間軸とがずらされる。これによって、現実の時間におけるＴｉ年目からＴｊ年目までにおけるシリンダカバーの損傷確率を正確に算出することができる。
【００７２】
次に、尤度計算の一般化について説明する。
現在ｍｎ年目であり、ｍ１〜ｍｎ年目にメンテナンスを行ったとする。この場合のトータル尤度は、下式（１７）と表すことができる。
【００７３】
【数１７】
2009288854.tif 000019
【００７４】
そして、式（１７）において積算される右辺の各値は、下式（１８）のように一般化される。
【００７５】
【数１８】
2009288854.tif 000020
【００７６】
ここで、ｍｉ回目のメンテナンスがシリンダカバーの置き換え及び交換を伴わない場合は、下式（１９）となり、ｍｉ回目のメンテナンスが置き換えであった場合には、下式（２０）となり、ｍｉ回目のメンテナンスが交換であった場合には、下式（２１）となる。
【００７７】
【数１９】
2009288854.tif 000021
【００７８】
【数２０】
2009288854.tif 000022
【００７９】
【数２１】
2009288854.tif 000023
【００８０】
この場合、式（１８）における分母は、下式（２２）で表すことができる。
【００８１】
【数２２】
2009288854.tif 000024
【００８２】
一方、分子は、ｍｉ＋１回目の検査結果によって以下のように分類できる。
ｍｉ＋１回目の検査結果が問題なしの場合には下式（２３）となり、ｍｉ＋１回目の検査結果が３年後に損傷の場合には下式（２４）となり、ｍｉ＋１回目の検査結果が損傷ありの場合には下式（２５）となる。
【００８３】
【数２３】
2009288854.tif 000025
【００８４】
【数２４】
2009288854.tif 000026
【００８５】
【数２５】
2009288854.tif 000027
【００８６】
このように尤度計算は、一般化することができる。
【００８７】
以上のような本実施形態の損傷確率算出装置及び損傷確率算出方法によれば、評価対象物であるシリンダカバーの通常検査時に取得された検査データに基づいて当該シリンダカバーに関する損傷率関数が補正され、補正された損傷率関数に基づいて損傷確率が算出される。
このため、通常検査にてシリンダカバーに損傷が発生していない場合であっても、当該シリンダカバーの損傷率関数が補正され、この補正された損傷率関数に基づいて損傷確率が算出される。
したがって、本実施形態の損傷確率算出装置及び損傷確率算出方法によれば、リスク評価の対象であるシリンダカバーに実際に損傷が発生する前であっても、現在用いられているシリンダカバーにより即した損傷確率に基づいてリスク評価を行うことが可能となる。
【００８８】
また、シリンダカバーの損傷確率とシリンダカバーが損傷することによる被害の影響度とに基づいて、シリンダカバーに対するリスクを評価し、このリスクに基づいてシリンダカバーに対するメンテナンスを行う場合に、上述のような損傷確率算出方法を用いてシリンダカバーの損傷確率を算出することによって、シリンダカバーに対するメンテナンスの期間を延ばすことが可能な場合がある。
つまり、シリンダカバーに対する通常検査において、シリンダカバーに損傷が見つからない場合には、図３に示すように、損傷率関数の傾きが補正されることによって緩やかとなる。この結果、損傷確率が低下し、リスクが低減することとなる。よって、シリンダカバーに対するメンテナンスの期間が延ばされることとなる。
【００８９】
なお、上述の損傷確率算出装置がコンピュータである場合には、コンピュータを、上述の記憶部１（記憶手段）と、損傷率関数算出部２（損傷率算出手段）、損傷率関数補正部３（損傷率関数補正手段）、損傷確率算出部４（損傷確率算出手段）として機能させるプログラムをコンピュータに格納し、このプログラムを実行することによって、上記損傷確率算出装置を実現しても良い。
このようなプログラムは、記録メディアやネットワーク上の記憶装置に格納され、ドライブ装置やネットワークを介してコンピュータにインストールされる。
【００９０】
以上、添付図面を参照しながら本発明に係る損傷確率算出装置、損傷確率算出方法、メンテナンス方法及び損傷確率算出プログラムの好適な実施形態について説明したが、本発明は、上記実施形態に限定されないことは言うまでもない。上述した実施形態において示した各構成部材の諸形状や組み合わせ等は一例であって、本発明の主旨から逸脱しない範囲において設計要求等に基づき種々変更可能である。
【００９１】
例えば、上記実施形態においては、本発明の評価対象物が、船が備える機関のシリンダカバーである構成について説明した。
しかしながら、本発明はこれに限定されるものではなく、リスク評価すべき全ての装置や部品等の対象物に対する損傷確率を算出する際に用いることができる。
【００９２】
また、上記実施形態においては、検査データに基づいて損傷率関数をベイズの定理を用いて補正する方法について説明した。
しかしながら、本発明はこれに限定されるものではなく、検査データに基づいて他の方法を用いて損傷率関数を補正しても良い。
【図面の簡単な説明】
【００９３】
【図１】本発明の一実施形態における損傷確率算出装置の機能構成を示すブロック図である。
【図２】本発明の一実施形態における損傷確率算出方法を説明するためのフローチャートである。
【図３】本発明の一実施形態における損傷確率算出方法での損傷率関数の補正方法を説明するためのグラフである。
【符号の説明】
【００９４】
Ｓ１……損傷確率算出装置、１……記憶部、２……損傷率関数算出部、３……損傷率関数補正部、４……損傷確率算出部
【図１】
2009288854.tif 000028
【図２】
2009288854.tif 000029
【図３】
2009288854.tif 000030