JP5592247B2 - Recovery curve creation system, recovery curve creation method and program - Google Patents
Recovery curve creation system, recovery curve creation method and program Download PDFInfo
- Publication number
- JP5592247B2 JP5592247B2 JP2010288012A JP2010288012A JP5592247B2 JP 5592247 B2 JP5592247 B2 JP 5592247B2 JP 2010288012 A JP2010288012 A JP 2010288012A JP 2010288012 A JP2010288012 A JP 2010288012A JP 5592247 B2 JP5592247 B2 JP 5592247B2
- Authority
- JP
- Japan
- Prior art keywords
- damage probability
- damage
- recovery
- component
- calculation unit
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- 238000011084 recovery Methods 0.000 title claims description 311
- 238000000034 method Methods 0.000 title claims description 41
- 230000006378 damage Effects 0.000 claims description 446
- 238000004364 calculation method Methods 0.000 claims description 255
- 238000004519 manufacturing process Methods 0.000 claims description 97
- 230000008569 process Effects 0.000 claims description 30
- 230000004308 accommodation Effects 0.000 claims description 8
- 238000012545 processing Methods 0.000 claims description 7
- 239000000470 constituent Substances 0.000 claims 5
- 230000001133 acceleration Effects 0.000 description 37
- 230000008859 change Effects 0.000 description 16
- 238000007726 management method Methods 0.000 description 9
- 230000006870 function Effects 0.000 description 7
- 238000001514 detection method Methods 0.000 description 6
- 238000012502 risk assessment Methods 0.000 description 6
- 239000000284 extract Substances 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 238000006073 displacement reaction Methods 0.000 description 3
- 230000004044 response Effects 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- 238000004891 communication Methods 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 230000003247 decreasing effect Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 239000011229 interlayer Substances 0.000 description 2
- 230000002093 peripheral effect Effects 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000001066 destructive effect Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- 238000011027 product recovery Methods 0.000 description 1
- 230000008439 repair process Effects 0.000 description 1
- 230000009528 severe injury Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Description
本発明は、一つの企業の複数の事業所などの地震被害を定量的に評価し、企業全体としての事業の継続性を予測するための復旧曲線作成システム、復旧曲線作成方法及びプログラムに関するものである。 The present invention relates to a restoration curve creation system, a restoration curve creation method, and a program for quantitatively evaluating earthquake damage at a plurality of establishments of one company and predicting the continuity of the business as a whole company. is there.
ビジネス影響度分析（ＢＩＡ：Business Impact Analysis）手法とは、事業継続（ＢＣ：Business Continuity）のＰＤＣＡ（Plan do check act cycle）サイクルでは計画段階に対応し、この段階で初期に行う最も重要な分析の一つで、不測の事態が発生した場合のビジネスへ影響を定量的あるいは定性的に評価する分析手法である。
このＢＩＡ手法を用いて、地震等の災害などによる事業への影響を少なくするために事業継続計画（ＢＣＰ：Business Continuity Plan）が策定されている。
The Business Impact Analysis (BIA) method is the most important analysis that is performed at the initial stage at this stage corresponding to the planning stage in the PDCA (Plan do check act cycle) cycle of Business Continuity (BC). One of these is an analysis method that quantitatively or qualitatively evaluates the business impact when an unexpected situation occurs.
Using this BIA method, a business continuity plan (BCP) has been formulated to reduce the impact on the business due to disasters such as earthquakes.
しかし、現在のところ、一般的なＢＩＡ手法が確立されていない。このため、コンサルティング会社の多くは、ＢＩＡを提供する際に、業種ごとにカスタマイズした質問項目を用いて、経営層やリスク管理者、施設管理者などにヒアリング調査を行っている。
このヒアリング調査の結果を、経験的な指標を用いて点数化することにより、その会社のビジネスへの影響度合いを定量化している。
However, at present, a general BIA method has not been established. For this reason, many consulting companies conduct interview surveys with management, risk managers, facility managers, etc., using question items customized for each industry when providing BIA.
The results of this interview survey are scored using empirical indicators to quantify the degree of impact on the company's business.
上述したように、コンサルティング会社が行うＢＩＡにおいては、定量的な計算式などが重要視されている訳ではなく、数値的な根拠が明確でない場合が多い。
このため、コンサルティング会社の判定において、ＢＩＡを依頼した会社の経営判断が強く働くことにより、経営方針が反映された重要業務、目標復旧時間、対策優先順位、復旧優先順位が決められる場合が多い。
As described above, in BIA conducted by a consulting company, quantitative calculation formulas are not regarded as important, and the numerical basis is often not clear.
For this reason, in the judgment of the consulting company, the important business, target recovery time, countermeasure priority, and recovery priority reflecting the management policy are often determined by strong management judgment of the company that requested BIA.
一方、建築・不動産業界においては、不動産の証券化や企業のリスクマネージメント強化の動きに伴い、建物のリスク分析が行われている。このリスク分析は、将来発生する可能性のある大災害の発生確率に基づいて、保有資産に対して与えられるであろう経済的損失の大きさを分析する。
この際、経済的損失を表す指標として、最大予想損失率（ＰＭＬ：Probable Maximum Loss）を用い、このＰＭＬにより、災害に対する対策効果の定量化、災害に対する最適な対策方法の選定、優先度の判定などが行える（例えば、特許文献１参照）。
On the other hand, in the construction and real estate industry, building risk analysis is being carried out in line with movements to securitize real estate and strengthen corporate risk management. This risk analysis analyzes the magnitude of economic loss that would be inflicted on assets held based on the probability of a catastrophe that may occur in the future.
At this time, the maximum expected loss rate (PML: Probable Maximum Loss) is used as an indicator of economic loss, and this PML is used to quantify the effectiveness of countermeasures against disasters, to select optimal countermeasure methods against disasters, and to determine priority levels. (For example, refer to Patent Document 1).
ビジネスに対するリスク分析という観点から、リスク分析の対象となる施設の復旧時間が、リスクの対応を画策する上で非常に重要なファクターとなる。
一方、特許文献１に示すＰＭＬのみを用いた分析では、施設の復旧を明確に求めることが、災害のビジネスへの影響度合いを推定することができない。
このため、災害及びテロなどの不測の事態が発生した際、この不測の事態による損害のビジネスへの影響を定量的、定性的に評価するために、復旧度を時間経過で示した復旧曲線が必要となる。
しかしながら、上述した経営方針が反映されて経験的に復旧曲線が作成されるため、定量的に精度の高い復旧曲線が作成されていない。
From the perspective of risk analysis for business, the recovery time of facilities subject to risk analysis is a very important factor in planning risk responses.
On the other hand, in the analysis using only the PML shown in
Therefore, in the event of an unforeseen event such as a disaster or terrorism, a recovery curve showing the degree of recovery over time is used to quantitatively and qualitatively evaluate the business impact of damage due to this unforeseen event. Necessary.
However, since the above-mentioned management policy is reflected and a recovery curve is created empirically, a quantitatively highly accurate recovery curve is not created.
この発明は上述した課題を解決するためになされたもので、リスク分析に用いる復旧曲線を定量的に生成する復旧曲線作成システム、復旧曲線作成方法及びプログラムを提供する。 The present invention has been made to solve the above-described problems, and provides a restoration curve creation system, a restoration curve creation method, and a program that quantitatively generate a restoration curve used for risk analysis.
本発明の復旧曲線作成システムは、ハザードによる施設に対する被害を定量的に評価し、施設における事業の継続性を予測する復旧曲線を生成する復旧曲線生成システムであり、
診断対象の施設における少なくとも建物および各収容設備を構成要素とし、当該構成要素がハザードの大きさにより損傷を受ける損傷確率を、当該構成要素毎に記憶する損傷確率記憶部と、前記構成要素の全壊における最大復旧時間に対し、前記損傷確率を乗じ、前記構成要素毎の最大復旧期待時間を求めるＰＭＴ算出部と、所定の周期毎に、前記構成要素毎の全壊における最大損失額に対し、前記損傷確率を乗じて最大予想損失額を求め、前記最大予想損失額を前記構成要素の再構成費で除算して予想損失率を求め、１から前記予想損失率を減算し、減算結果を操業度として出力する操業度算出部と、前記操業度を前記周期毎にプロットし、復旧曲線を生成する復旧曲線生成部とを有し、前記操業度算出部が、前記周期の経過毎に経過時間をカウントし、当該経過時間が前記最大復旧期待時間を超えた前記構成要素の損傷確率を０とし、前記予想損失率を再度計算することを特徴とする。
The recovery curve creation system of the present invention is a recovery curve generation system that quantitatively evaluates damage to facilities due to hazards and generates a recovery curve that predicts business continuity in the facility,
A damage probability storage unit that stores, for each component, the damage probability that the component is damaged by the size of the hazard, at least a building and each accommodation facility in the facility to be diagnosed, and complete destruction of the component The PMT calculation unit for obtaining the maximum expected recovery time for each component by multiplying the maximum recovery time at the maximum recovery time by the damage probability, and for the maximum loss amount at the total destruction for each component, the damage Multiply the probability to obtain the maximum expected loss amount, divide the maximum expected loss amount by the reconfiguration cost of the component to obtain the expected loss rate, subtract the expected loss rate from 1, and use the subtraction result as the operating rate An operation rate calculation unit for outputting, and a recovery curve generation unit for plotting the operation rate for each cycle and generating a recovery curve, wherein the operation level calculation unit is an elapsed time for each elapse of the cycle. Counting, damage probability of the component to which the elapsed time exceeds the maximum recovery expected time to 0, and calculates the expected loss rate again.
本発明の復旧曲線作成システムは、前記周期の経過毎に経過時間をカウントし、前記構成要素毎に、当該経過時間を前記最大復旧期待時間で除算した結果に対して当該損傷確率を乗算し、この乗算した結果を前記損傷確率から減算して前記経過時間に対応した新たな損傷確率を算出し、当該算出された新たな損傷確率を前記操業度算出部へ出力する損傷確率計算部をさらに有することを特徴とする。 The recovery curve creation system of the present invention counts the elapsed time for each elapse of the cycle, and multiplies the damage probability for the result of dividing the elapsed time by the maximum expected recovery time for each component. A damage probability calculation unit that subtracts the result of multiplication from the damage probability to calculate a new damage probability corresponding to the elapsed time, and outputs the calculated new damage probability to the operation degree calculation unit. It is characterized by that.
本発明の復旧曲線作成システムは、前記操業度算出部が、フォールトツリーにおいて前記操業度を求める際に、複数の構成要素を組み合わせて構成される上位構成要素の最大復旧期待時間を、組み合わせた構成要素の事象に対する上位構成要素の事象がＡＮＤ事象の場合及びＯＲ事象の場合のいずれにおいても、下位の前記構成要素のいずれか大きい方の最大復旧時間を用いることを特徴とする。 The recovery curve creation system of the present invention has a configuration in which the operation level calculation unit combines the maximum expected recovery time of the upper component configured by combining a plurality of components when obtaining the operation level in the fault tree. In both cases where the event of the upper component with respect to the event of the element is an AND event and an OR event, the largest maximum recovery time of the lower component is used.
本発明の復旧曲線作成システムは、前記周期毎の前記構成要素を用いて生産する商品の生産額に対し、各構成要素の前記損傷確率を前記生産額に対して乗じた結果を、前記構成要素の前記最大復旧時間を超えるまでの期間において積算し、売り上げの損失額を求める売上損失額算出部をさらに有することを特徴とする。 The restoration curve creation system of the present invention is configured by multiplying the production value by the damage probability of each component with respect to the production value of a product produced using the component for each period. And a sales loss amount calculating unit for calculating a loss amount of sales by accumulating in a period until the maximum recovery time is exceeded.
本発明の復旧曲線作成方法は、施設に対する地震被害を定量的に評価し、施設における事業の継続性を予測するとともに、地震対策を施した場合の対策効果の予測を行う事業継続予測システムの制御方法であり、診断対象の施設における少なくとも建物および各収容設備を構成要素とし、当該構成要素が所定の最大地動加速度で損傷を受ける損傷確率を、当該構成要素毎に記憶する損傷確率記憶部と、ＰＭＴ算出部が、診断対象の施設における少なくとも建物および各収容設備を構成要素とし、当該構成要素が所定の最大地動加速度で損傷を受ける損傷確率を、当該要素毎に記憶する損傷確率記憶部から、前記構成要素の前記損傷確率を読み出し、前記構成要素の全壊における最大復旧時間に対し、読み出した前記損傷確率を乗じ、前記構成要素毎の最大復旧期待時間を求めるＰＭＴ算出過程と、操業度算出部が、所定の周期毎に、前記損傷確率記憶部から構成要素の前記損傷確率を読み出し、前記構成要素毎の全壊における最大損失額に対し、読み出した前記損傷確率を乗じて最大予想損失額を求め、前記最大予想損失額を前記構成要素の再構成費で除算して予想損失率を求め、１から前記予想損失率を減算し、減算結果を操業度として出力する操業度算出過程と、
復旧曲線生成部が前記操業度を前記周期毎にプロットし、復旧曲線を生成する復旧曲線生成過程とを有し、前記操業度算出部が、前記周期の経過毎に経過時間をカウントし、当該経過時間が前記最大復旧期待時間を超えた前記構成要素の損傷確率を０とし、前記予想損失率を再度計算することを特徴とする。
The restoration curve creation method of the present invention quantitatively evaluates earthquake damage to a facility, predicts business continuity in the facility, and controls a business continuity prediction system that predicts the effect of countermeasures when earthquake countermeasures are taken. A damage probability storage unit that stores at least a building and each accommodation facility in a facility to be diagnosed as a component, and stores the damage probability that the component is damaged at a predetermined maximum ground motion acceleration for each component; From the damage probability storage unit that stores, as a component, at least a building and each accommodation facility in the facility to be diagnosed as a component, and the damage probability that the component is damaged at a predetermined maximum ground acceleration, for each component. Read the damage probability of the component, multiply the maximum recovery time in the total destruction of the component by the read damage probability, and A PMT calculation process for obtaining a maximum expected recovery time for each element, and an operation degree calculation unit reads the damage probability of the component element from the damage probability storage unit for each predetermined period, and the maximum loss in total destruction for each component element Multiply the read damage probability to obtain the maximum expected loss amount, divide the maximum expected loss amount by the reconfiguration cost of the component to obtain the expected loss rate, and subtract the expected loss rate from 1 An operation degree calculation process for outputting the subtraction result as an operation degree;
A recovery curve generation unit plots the operation rate for each cycle, and a recovery curve generation process for generating a recovery curve, and the operation level calculation unit counts elapsed time for each elapse of the cycle, The damage probability of the component whose elapsed time exceeds the maximum expected recovery time is set to 0, and the expected loss rate is calculated again.
本発明のプログラムは、施設に対する地震被害を定量的に評価し、施設における事業の継続性を予測するとともに、地震対策を施した場合の対策効果の予測を行う事業継続予測システムの動作をコンピュータに実行させるプログラムであり、診断対象の施設における少なくとも建物および各収容設備を構成要素とし、当該構成要素が所定の最大地動加速度で損傷を受ける損傷確率を、当該構成要素毎に記憶する損傷確率記憶部と、ＰＭＴ算出部が、診断対象の施設における少なくとも建物および各収容設備を構成要素とし、当該構成要素が所定の最大地動加速度で損傷を受ける損傷確率を、当該要素毎に記憶する損傷確率記憶部から、前記構成要素の前記損傷確率を読み出し、前記構成要素の全壊における最大復旧時間に対し、読み出した前記損傷確率を乗じ、前記構成要素毎の最大復旧期待時間を求めるＰＭＴ算出処理と、操業度算出部が、所定の周期毎に、前記損傷確率記憶部から構成要素の前記損傷確率を読み出し、前記構成要素毎の全壊における最大損失額に対し、読み出した前記損傷確率を乗じて最大予想損失額を求め、前記最大予想損失額を前記構成要素の再構成費で除算して予想損失率を求め、１から前記予想損失率を減算し、減算結果を操業度として出力する操業度算出処理と、復旧曲線生成部が前記操業度を前記周期毎にプロットし、復旧曲線を生成する復旧曲線生成処理と、前記操業度算出部が、前記周期の経過毎に経過時間をカウントし、当該経過時間が前記最大復旧期待時間を超えた前記構成要素の損傷確率を０とし、前記予想損失率を再度計算する処理とをコンピュータに実行させるプログラムである。 The program of the present invention quantitatively evaluates earthquake damage to a facility, predicts business continuity in the facility, and operates the operation of a business continuity prediction system that predicts the countermeasure effect when earthquake countermeasures are taken. A damage probability storage unit that is a program to be executed and has at least a building and each accommodation facility in a facility to be diagnosed as components, and stores a damage probability that the component is damaged at a predetermined maximum ground motion acceleration for each component And a PMT calculation unit having at least a building and each accommodation facility in the facility to be diagnosed as components, and a damage probability storage unit for storing the damage probability that the component is damaged at a predetermined maximum ground motion acceleration for each of the components The damage probability of the component is read from the maximum recovery time in the complete destruction of the component before reading A PMT calculation process for obtaining a maximum expected recovery time for each component by multiplying by the damage probability, and an operation degree calculation unit reads the damage probability of the component from the damage probability storage unit for each predetermined period, and The maximum expected loss amount is obtained by multiplying the maximum loss amount in total destruction for each element by the read probability of damage, and the expected loss rate is obtained by dividing the maximum expected loss amount by the reconfiguration cost of the component. Subtracting the expected loss rate from the output, output the subtraction result as the operation degree calculation process, recovery curve generation unit plots the operation degree for each period, recovery curve generation process to generate a recovery curve, A process in which the operation degree calculation unit counts the elapsed time every time the cycle elapses, sets the damage probability of the component whose elapsed time exceeds the maximum expected recovery time to 0, and recalculates the expected loss rate Is a program for executing a computer.
この発明によれば、最大損失額を構成要素の再構成費にて除算し、各構成要素のＰＭＬを所定の周期毎に算出し、このＰＭＬから操業度を算出し、この求めた操業度を周期的にプロットし、復旧曲線を求めるため、経営方針が人為的に反映されることなく、リスク分析に用いる復旧曲線を定量的に、かつ高い精度で生成することができる。 According to the present invention, the maximum loss amount is divided by the reconfiguration cost of the component, the PML of each component is calculated for each predetermined period, the operation rate is calculated from the PML, and the obtained operation rate is calculated. Since the recovery curve is obtained by periodically plotting, the recovery curve used for risk analysis can be generated quantitatively and with high accuracy without artificially reflecting the management policy.
以下、図面を参照して、本発明の実施の形態について説明する。図１は、この発明の一実施形態による復旧曲線作成システムの構成例を示す概略ブロック図である。
図１に示す復旧曲線作成システム１００は、コンピュータに周辺デバイスが接続されており、損傷確率算出部１１、ＰＭＴ算出部１２、操業度算出部１３、復旧曲線生成部１４、入力装置１５、表示装置１６、損傷確率データベース１７、基本情報データベース１８及びモデルデータベース１９を備えている。ここで、入力装置１５は、ユーザーが、復旧曲線作成システム１００に与える命令やデータを入力するための装置であって、具体的には、キーボードやマウス等である。表示装置１６は、復旧曲線作成システム１００の各部によって算出された算出結果等が表示される装置であって、具体的には、ＣＲＴ（Cathode Ray Tube；陰極線管）ディスプレイや液晶ディスプレイである。
Embodiments of the present invention will be described below with reference to the drawings. FIG. 1 is a schematic block diagram showing a configuration example of a recovery curve creation system according to an embodiment of the present invention.
A recovery curve creation system 100 shown in FIG. 1 has peripheral devices connected to a computer, and includes a damage
損傷確率算出部１１は、最大地動加速度と、事業構成要素が小破、中破、または大破・全壊の損傷を受ける確率である事業構成要素の損傷確率を算出する。また、この事業構成要素の損傷確率は、復旧曲線作成システム１００の利用者が入力装置１５から、復旧曲線の作成を行う際に逐次入力するようにしても良い。この場合、損傷確率算出部１１は、入力される損傷確率を事業構成要素の識別番号とともに損傷確率データベース１７に書き込み、記憶させる。ここで、事業構成要素とは、例えば、工場における建物、電源、部品を作成するライン（製造機械）等の製品を製造するために必要なインフラの各要素を示している。ここでは、損傷確率を地震に起因するものとして説明しているが、以下に述べる復旧曲線の生成アルゴリズムは、水害、台風の被害、テロによる破壊行動などのハザードによる損傷確率を求めるフラジリティ曲線を生成し、損傷確率、再構成価格、復旧時間などを定義することにより、いずれのハザードに対しても対応させることができる。
The damage
操業度算出部１３は、事業構成要素毎に、事業構成要素の最大損失額に対し、損傷確率を乗算することにより、最大予想損失額を算出する。
また、操業度算出部１３は、事業構成要素毎に、事業構成要素を新たに構成するための費用である再構成費（部品の再調達額、再工事費、営業損失など）により、すでに求めた最大予想損失額を除算し、除算結果として最大予想損失率を求め、この算出した最大予想損失率を用い、この最大予想損失率（％）を通常の操業率を１００％から減算し、すなわち以下の式により操業度（％）を算出する。
操業度（％）＝１００−最大予想損失率（％）
The operation
In addition, the operation
Operating rate (%) = 100-Maximum expected loss rate (%)
ＰＭＴ算出部１２は、最大復旧時間に対し、損傷確率を乗算し、乗算結果として最大復旧期待時間（ＰＭＴ）を求める。この最大復旧期待時間は、各事業構成要素（あるいはサブ事業構成要素）損傷確率に対応した損傷度合いである場合の復旧時間の最大値を示している。
損傷確率データベース１７は、工場にある事業構成要素（あるいはサブ事業構成要素）毎の損傷確率が、地震が発生した際の最大地動加速度毎に記憶されている（詳細は後述）。
基本情報データベース１８は、各工場の事業構成要素毎の操業度を算出する際に必要な損傷確率を求めるために用いられる、建物のフラジリティ曲線を作成するために用いられる基本データが格納されている（詳細は後述）。
モデルデータベース１９は、復旧曲線を生成する際に、損傷確率算出部１１が算出した損傷確率や、ＰＭＴ算出部１２が算出した最大復旧期待時間等を含むモデルテーブルが記憶されている。
The
The
The
The
＜損傷確率の算出＞
次に、損傷確率算出部１１における損傷確率の算出の説明を行う。
復旧曲線作成システムは工場における製造ライン及び電力設備を含めたトータルなリスクからの復旧の推定を対象としている。図２は、本実施形態における復旧曲線を推定する対象である工場の構成要素（事業構成要素）のフォールトツリー（ＦＴ）をモデル化した例である。すなわち、この工場は、製造ラインと、電力設備とが事業構成要素（コンポーネント）となっている。さらに、製造ラインは機械１と機械２との下位のサブ事業構成要素（サブコンポーネント）から構成され、電力設備は商用電源と自家発電源との下位のサブ事業構成要素から構成されている。
このフォールトツリーにおいて、製造ラインは機械１及び機械２のいずれか１つが損傷を受けると機能しなくなる。すなわち、機械１及び機械２は、製造ラインに対するＯＲ（オア：論理和）事象となっている。
一方、電力設備は商用電源及び自家発電源の双方ともに損傷を受けた場合に、電力設備としての機能を失う。すなわち、商用電源及び自家発電源は、電力設備に対するＡＮＤ（アンド：論理積）事象となっている。
<Calculation of damage probability>
Next, calculation of damage probability in the damage
The restoration curve creation system is intended for the estimation of restoration from the total risk including the production line and power equipment in the factory. FIG. 2 is an example of modeling a fault tree (FT) of a component (business component) of a factory that is a target for estimating a recovery curve in the present embodiment. That is, in this factory, a production line and electric power equipment are business components (components). Further, the production line is composed of sub-business components (sub-components) of the lower level of the
In this fault tree, the production line fails when one of
On the other hand, the power equipment loses its function as a power equipment when both the commercial power source and the self-generated power source are damaged. That is, the commercial power source and the self-generated power source are AND (AND: logical product) events for the power equipment.
次に、図３は損傷確率と最大地動加速度との関係を示したフラジリティ曲線のグラフである。なお、このグラフは、耐震性能を示すＩｓ値（構造耐震指標）が０．６の機械１のフラジリティ曲線の例である。この図３のフラジリティ曲線のデータは、損傷確率データベース１７に記憶されている。このグラフの横軸は最大地動加速度ＰＧＡ（Peak Ground Acceleration、単位はｃｍ／ｓ２）、縦軸は損傷確率である。例えば、Ｉｓ値＝０．６の場合、最大地動加速度ＰＧＡが６５０ｃｍ／ｓ２では、損傷の程度が大破・全壊に至る確率は０．３３９すなわち３３．９％である。また、損傷の程度が中破に至る確率は、中破以上となる損傷確率と大破・全壊に至る損傷確率との差であるため、０．２４すなわち２４％になる。上述では、機械１や機械２が１台であり、それぞれが損傷する確率として損傷確率として求める説明となっているが、例えば機械１がＮ台あり、このＮ台の内、Ｍ（≦Ｎ）台が損傷を受け、（Ｎ−Ｍ）台が稼動可能な状態の機械１の集合の損傷確率とする場合も含む。
Next, FIG. 3 is a graph of a fragility curve showing the relationship between the damage probability and the maximum ground motion acceleration. In addition, this graph is an example of the fragility curve of the
すなわち、損傷確率算出部１１は、入力装置から入力される最大地動加速度に対応した損傷確率を、損傷確率データベース１７に事業構成要素（あるいはサブ事業構成要素）毎に設定された図３のフラジリティ曲線から読み出し、事業構成要素毎の損傷確率を求める。損傷確率算出部１１は、サブ事業構成要素毎に要素識別情報（例えば、サブ事業構成要素名や番号など）を付加し、この要素識別情報とともに、この要素識別情報に対応するフラジリティ曲線のグラフを、損傷確率データベース１７へ書き込み、記憶させる。そして、損傷確率算出部１１は、損傷確率データベース１７から各事業構成要素のフラジリティ曲線のグラフを検索する際、要素識別情報により検索して、要素識別情報に対応するフラジリティ曲線のグラフを読み出す。
また、損傷確率は、上述したように、損傷確率データベース１７における各事業構成要素のフラジリティ曲線のグラフを用いずに、事業構成要素毎の構造調査結果に基づいて、復旧曲線を求めようとするＢＩＡ実施者が直接に入力装置１５から入力するようにしても良い。事業構成要素がサブ事業構成要素に分割されない場合、事業分割要素のフラジリティ曲線が損傷確率データベース１７に記憶されて、設定されている。
That is, the damage
Further, as described above, the damage probability is a BIA that seeks a recovery curve based on the structural investigation result for each business component without using the graph of the fragility curve of each business component in the
図４は、上記の図３に示した建物のフラジリティ曲線を作成するために用いられた基本データを示す表である。この基本データも、地震リスク診断プログラムに含まれている。この表は、建物の耐震性能を示すＩｓ値（構造耐震指標）別に、小破、中破、大破・全壊の損傷確率を算出するために用いた中央値と、対数標準偏差とを示している。 FIG. 4 is a table showing basic data used to create the fragility curve of the building shown in FIG. This basic data is also included in the earthquake risk diagnosis program. This table shows the median and logarithmic standard deviation used to calculate the damage probability for small, medium, large, and total damage by Is value (structural earthquake resistance index) that indicates the earthquake resistance of the building. .
この表において、Ｉｓ＝０．６の建物の損傷の程度（損傷度）は、小破、中破または大破・全壊の中央値（最大地動加速度）を、阪神・淡路大震災において建物が受けた損傷の程度と、最大推定地動加速度との関係に基づいて定めている。 In this table, the damage level (damage) of the building with Is = 0.6 is the median value (maximum ground motion acceleration) of small, medium, large or total damage, and the damage that the building suffered in the Great Hanshin-Awaji Earthquake. And the maximum estimated ground motion acceleration.
また、Ｉｓ＝０．６以外の建物の損傷の程度は、小破、中破または大破・全壊の中央値（最大地動加速度）Ａを、下記の式（１）により算出する。
Ａ＝Ａｏ×（Ｉｓ／０．６） …（１）
ただし、Ａｏは、Ｉｓ＝０．６の場合の、小破、中破または大破・全壊の中央値（最大地動加速度）、Ａは、任意のＩｓ値の建物の損傷の程度が、小破、中破または大破・全壊の中央値（最大地動加速度）である。なお、上記の式（１）による、任意のＩｓ値の建物の最大地動加速度Ａの算出は、損傷確率算出部１１が行う。
As for the degree of damage to buildings other than Is = 0.6, the median value (maximum ground motion acceleration) A of small breakage, medium breakage, or large breakage / total destruction is calculated by the following equation (1).
A = Ao × (Is / 0.6) (1)
However, Ao is the median value (maximum ground motion acceleration) of small breakage, medium breakage or large breakage / total destruction when Is = 0.6, and A is the degree of damage of the building of any Is value, This is the median value (maximum ground acceleration) for moderately damaged, severely damaged or completely destroyed. The damage
例えば、Ｉｓ＝０．３の、小破の中央値（最大地動加速度）Ａは、Ｉｓ＝０．６の建物の、小破の中央値（最大地動加速度）Ａｏが４００ｃｍ／ｓ２なので、Ａ＝４００×（０．３／０．６）＝２００ｃｍ／ｓ２のようにして算出される。
なお、本実施形態における復旧曲線作成システムにおいて、建物の損傷の程度は、Is値以外にも、建物の応答変位（免震建物応答変位）、層間変位、層間変形角などから算出する方法もあり、上述したＩｓ値による場合と同様に、損傷確率を算出する構成としても良い。
また、上述した損傷確率は地震の最大地動加速度との関係のフラジリティ曲線を例にとり説明したが、所定のハザードとして、横軸を洪水の場合の浸水水位、停電の場合の停電時間等、縦軸を損傷確率としたフラジリティ曲線を用い、浸水水位または停電時間から損傷確率を求める構成としても良い。
For example, the median value (maximum ground acceleration) A of Is = 0.3 is A, because the median value (maximum ground acceleration) Ao of the building of Is = 0.6 is 400 cm / s 2 . = 400 × (0.3 / 0.6) = 200 cm / s 2
In addition, in the restoration curve creation system in this embodiment, there is a method of calculating the degree of damage to the building from the response displacement of the building (base-isolated building response displacement), interlayer displacement, interlayer deformation angle, etc. in addition to the Is value. As in the case of using the Is value described above, the damage probability may be calculated.
In addition, the damage probability described above has been explained by taking the example of the fragility curve in relation to the maximum ground motion acceleration of the earthquake. It is also possible to use a fragility curve with the damage probability as the damage probability and obtain the damage probability from the flooded water level or power failure time.
＜操業度の算出＞
図５は、各最大地動加速度での損害額を示すグラフ（地震ロス関数）である。この損害額グラフのデータは、損傷確率データベース１７にサブ事業構成要素毎に記憶されている。このグラフの横軸は最大地動加速度、縦軸は最大損害額である。例えば、この損失額グラフによれば、最大地動加速度が８００ｃｍ／ｓ２での最大損害額は約２４億円である。この損失額グラフは、予め損傷確率データベース１７に、サブ事業構成要素毎に書き込まれ、記憶されている。
<Calculation of operation rate>
FIG. 5 is a graph (earthquake loss function) showing the amount of damage at each maximum ground motion acceleration. The data of the damage amount graph is stored in the
すなわち、損傷確率算出部１１は、入力装置１５から入力される最大地動加速度に対応した損失額を、損傷確率データベース１７に事業構成要素毎に設定された図５の損失グラフから読み出し、サブ事業構成要素毎の損害額とする。損傷確率算出部１１は、サブ事業構成要素毎に要素識別情報を付加し、この要素識別情報とともに、この要素識別情報に対応する損失グラフを、損傷確率データベース１７へ書き込み、記憶させる。そして、損傷確率算出部１１は、損傷確率データベース１７から各事業構成要素の損失グラフを検索する際、要素識別情報により検索して、要素識別情報に対応する損失グラフを読み出す。
また、損害額は、上述したように、損傷確率データベース１７における各事業構成要素の損失額グラフを用いずに、工場の関係者からヒアリングをした結果に基づいて、復旧曲線を求めようとするＢＩＡ実施者が直接に入力装置１５から入力するようにしても良い。事業構成要素がサブ事業構成要素に分割されない場合、事業分割要素の損失額グラフが損傷確率データベース１７に設定されている。
That is, the damage
In addition, as described above, the amount of damage is the BIA that seeks a recovery curve based on the results of interviews with the persons concerned at the factory without using the loss amount graph of each business component in the
そして、操業度算出部１３は、損傷確率算出部１１が事業構成要素毎に求めた損失額と損傷確率とを、それぞれ事業構成要素毎に乗算し（損失額（円）×損傷確率）、各事業構成要素の最大予想損失額（円）を算出する。
また、操業度算出部１３は、入力装置１５から入力される各事業構成要素の再構成費（円：設備の再調達額、設備の再工事額、営業損失など）を事業構成要素毎に、基本情報データベース１８に対して、図６に示す基本データテーブルのデータとして予め書き込み、記憶させておく。操業度算出部１３は、サブ事業構成要素毎に要素識別情報を付加し、この要素識別情報とともに、この要素識別情報に対応する基本データテーブルの再構成費のデータを、基本情報データベース１８へ書き込み、記憶させる。そして、操業度算出部１３は、基本情報データベース１８から各事業構成要素の再構成費のデータを検索する際、要素識別情報により検索して、要素識別情報に対応する再構成費のデータを読み出す。図６は、モデルとしての工場における事業構成要素（コンポーネント）とサブ事業構成要素（サブコンポーネント）との関係を示し、各サブ事業構成要素の再構成費や復旧時間が示された基本データテーブルの構成を示す図である。
Then, the operation
In addition, the operation
次に、操業度算出部１３は、再構成費により最大予想損失額を除算し、除算した結果を最大予想損失率（ＰＭＬ）として算出する。
そして、操業度算出部１３は、求めた最大予想損失率を用いて以下の操業度関係式から、すなわち操業が１００％行われている状態から最大予想損失率を減算し、減算結果を操業度（％）として算出する。
操業度（％）＝１００（％）−最大予想損失率（％）
操業度は、一般的に、商品を製造する機器の稼働率により算定されている。しかしながら、稼働率から算定される操業度は、事業継続という観点から考えると、売上げに対する寄与率や材料の仕入れ量、キャッシュフローなどに対するリスクを適格に寄与しているとはいえない。
そこで、本実施形態においては、経営的な地震を含むハザードを考慮した最大予想損失率を用いることにより、容易に通常時の操業度を表す上記操業度関係式を生成した。この操業度関係式により、通常時から、ハザードによる被害発生時における操業度までを適格に表現することができる。
Next, the operation
Then, the operation
Operating rate (%) = 100 (%)-Maximum expected loss rate (%)
The operation rate is generally calculated based on the operation rate of the equipment that manufactures the product. However, from the viewpoint of business continuity, the operating rate calculated from the occupancy rate does not qualify for the contribution to sales, the amount of materials purchased, and the risk to cash flow.
Therefore, in the present embodiment, the above operating degree relational expression representing the operating degree at the normal time is easily generated by using the maximum expected loss rate in consideration of the hazard including a management earthquake. By this operational degree relational expression, it is possible to properly express from the normal time to the operational degree at the time of occurrence of damage due to the hazard.
＜最大復旧期待時間（ＰＭＴ：Probable Maximum Time）の算出＞
ＰＭＴ算出部１２は、ＢＩＡ実施者の入力装置１５からの操作により、このＢＩＡ実施者が予め工場の関係者にヒアリングした全壊した際の復旧時間を、基本情報データベース１８に対して、事業構成要素毎に、図６に示す基本データテーブルのデータとして予め書き込み、記憶させておく。ＰＭＴ算出部１２は、サブ事業構成要素毎に要素識別情報を付加し、この要素識別情報とともに、この要素識別情報に対応する基本データテーブルの復旧時間のデータを、基本データベース１８へ書き込み、記憶させる。そして、ＰＭＴ算出部１２は、基本情報データベース１８から各事業構成要素の復旧時間のデータを検索する際、要素識別情報により検索して、要素識別情報に対応する復旧時間のデータを読み出す。
<Calculation of maximum expected recovery time (PMT: Probable Maximum Time)>
The
そして、ＰＭＴ算出部１２は、基本情報データベース１８の基本データテーブルから読み出した復旧時間（日）に対して、損傷確率を乗算し、乗算した結果を最大復旧期待時間（ＰＭＴ）とする。
最大の損傷を受けたときの復旧時間（建物の新築や破損した機器の交換など）を用い、この復旧時間に対し、損傷確率を乗算することにより、容易に被害の程度に対応して、復旧に要する時間として最大復旧期待時間を算出することができる。
Then, the
By using the recovery time (such as a new building or replacement of damaged equipment) when the greatest damage has been received, multiplying this recovery time by the damage probability makes it easy to respond to the degree of damage. The maximum expected recovery time can be calculated as the time required for.
＜復旧曲線の生成−リソース＞
［リソース（事業構成要素）に対する操業度の算出］
次に、図２のフォールトツリーにおいて、各フォールトの状態の損傷確率を記述した図７のフォールトツリーを用いて、復旧曲線の生成を説明する。
損傷確率算出部１１は、事業構成要素である機械１、機械２、商用電源及び自家発電源の各々の損傷確率を、入力された最大地動加速度に対応させ、損傷確率データベース１７に記憶されている機械１、機械２、商用電源及び自家発電源の各々のフラジリティ曲線から読み出す。
ここで、損傷確率算出部１１により読み出された損傷確率において、機械１の損傷確率が０．２４０２０であり、機械２の損傷確率が０．０７４４であり、商用電源の損傷確率が０．００７４０であり、自家発電源の損傷確率が０．２２８６である。
<Generation of recovery curve-resource>
[Calculation of operation rate for resources (business components)]
Next, the generation of the recovery curve will be described using the fault tree of FIG. 7 in which the damage probability of each fault state in the fault tree of FIG. 2 is described.
The damage
Here, in the damage probability read by the damage
次に、損傷確率算出部１１は、機械１及び機械２の各々の損傷確率により、製造ラインの損傷確率を求める。
すなわち、損傷確率算出部１１は、図７のフォールトツリーにおいて、製造ラインが機械１及び機械２のＯＲ事象として定義されているため、機械１の損傷確率を１から減算した数値（損傷しない確率）と、機械２の損傷確率を１から減算した数値（損傷しない確率）とを乗算し、いずれも損傷しない非損傷確率を算出する。
そして、損傷確率算出部１１は、１から求めた非損傷確率を減算し、機械１または機械２のいずれかが損傷する損傷確率（製造ラインの損傷確率）を算出する。
すなわち、以下の式を用いて、機械１の損傷確率Ｐ１及び機械２の損傷確率Ｐ２から製造ラインの損傷確率ＰＡを算出する。
ＰＡ＝１−（１−Ｐ１）・（１−Ｐ２）
＝１−（１−０．２４０２０）・（１−０．０７４４０）
＝１−０．７５９８・０．９２５６
＝０．２９６７３
Next, the damage
That is, the damage
Then, the damage
That is, the damage probability P A of the production line is calculated from the damage probability P 1 of the machine 1 and the damage probability P 2 of the machine 2 using the following formula.
P A = 1- (1-P 1 ) · (1-P 2 )
= 1- (1-0.24020). (1-0.07440)
= 1-0.7598.0.9256
= 0.29673
また、同様に、損傷確率算出部１１は、図７のフォールトツリーにおいて、電力設備がが商用電源及び自家発電源のＡＮＤ事象として定義されているため、商用電源の損傷確率と、自家発電源の損傷確率とを乗算し、商用電源及び自家発電源の双方が損傷する損傷確率を算出する。
すなわち、以下の式を用いて、商用電源の損傷確率Ｐ３及び自家発電源の損傷確率Ｐ４から電力設備の損傷確率ＰＢを算出する。
ＰＢ＝Ｐ３・Ｐ４
＝０．００７４０×０．２２８６０
＝０．００１６９１６
Similarly, the damage
That is, the damage probability P B of the power equipment is calculated from the damage probability P 3 of the commercial power source and the damage probability P 4 of the self-generated power source using the following formula.
P B = P 3・ P 4
= 0.00740 x 0.22860
= 0.0016916
次に、操業度算出部１３は、イベントツリーに基づいて、最大予想損失率（ＰＭＬ）を、以下の様に算出する。
図８は、工場における製造ライン及び電力設備の損傷有りまたは無しの組み合わせ毎の確率を示すイベントツリーのテーブルである。このテーブルにおいて、損傷有りは状態が「Ｙ」で示され、損傷無しは状態が「Ｎ」で示されている。
そして、操業度算出部１３は、このイベントツリーにおけるＮ（損傷無し）に、各事業構成要素が所定の最大地動加速度で損傷を受けない確率を当てはめ、Ｙ（損傷有り）に、各事業構成要素が所定の最大地動加速度で損傷を受ける損傷確率を当てはめる。各事業構成要素が損傷を受ける確率とは、すでに説明したフォールトツリーにおいて説明した、各構成要素の損傷確率であり、各構成要素が損傷を受けない確率とは、（１−各構成要素の損傷確率）である。
そして、演算処理装置２は、損傷状態の組み合わせ毎に、各事業構成要素、本実施形態においては、製造ライン及び電力設備の各々の損傷を受けない確率または損傷を受ける確率どうしの積をとり、所定の最大地動加速度による各事業構成要素の各確率により、工場が各損傷モードに至る確率を算出する。
Next, the operation
FIG. 8 is an event tree table showing the probabilities for each combination of whether or not the production line and power equipment in the factory are damaged or not. In this table, the state with damage is indicated by “Y”, and the state without damage is indicated by “N”.
Then, the operation
And for each combination of damage state, the
例えば、図８に示すように、操業度算出部１３は、製造ラインが損傷が無く、電力設備も損傷が無い場合、以下の確率算出式により、いずれにも損傷が無い確率ＰＮＮを算出する。
ＰＮＮ＝（１−ＰＡ）・（１−ＰＢ）
＝０．７０２０８
また、操業度算出部１３は、製造ラインが損傷が無く、電力設備も損傷がある場合、以下の確率算出式により、いずれにも電力設備のみに損傷がある確率ＰＮＹを算出する。
ＰＮＹ＝（１−ＰＡ）・ＰＢ
＝０．００１８９
また、操業度算出部１３は、製造ラインが損傷が有り、電力設備に損傷が無い場合、以下の確率算出式により、製造ラインにのみ損傷がある確率ＰＹＮを算出する。
ＰＹＮ＝ＰＡ・（１−ＰＢ）
＝０．２９６２２
また、操業度算出部１３は、製造ライン及び電力設備の双方に損傷がある場合、以下の確率算出式により、いずれにも損傷がある確率ＰＹＹを算出する。
ＰＹＹ＝ＰＡ・ＰＢ
＝０．０００５０
For example, as shown in FIG. 8, when the production line is not damaged and the power equipment is not damaged, the operation
PNN = (1-P A ) · (1-P B )
= 0.70208
Further, when the production line is not damaged and the power equipment is also damaged, the operation
P NY = (1-P A ) · P B
= 0.00189
Further, when the production line is damaged and the power equipment is not damaged, the operation
P YN = P A · (1 -P B)
= 0.29622
In addition, when both the production line and the power equipment are damaged, the operation
P YY = P A・ P B
= 0.00050
次に、操業度算出部１３は、損傷確率を算出した最大地動加速度による各事業構成要素の損失額を、基本情報データベース１８の基本データテーブルから読み出す。
本実施形態の場合、例えば、製造ラインの損失額が２０００万円であり、電力設備の損失額が１５００万円である。説明を簡単にするため、再構成費も製造ラインが２０００万円とし、電力設備が１５００万円とする。
そして、操業度算出部１３は、製造ラインの損失額と、電力設備の損失額とを加算し、製造ラインと電力設備との双方が全壊した際の総損失額を算出する。
Next, the operation
In the case of the present embodiment, for example, the loss amount of the production line is 20 million yen, and the loss amount of the power equipment is 15 million yen. In order to simplify the explanation, the reconfiguration cost is assumed to be 20 million yen for the production line and 15 million yen for the power equipment.
And the operation
次に、操業度算出部１３は、製造ライン及び電力設備の損傷の有無によるそれぞれの組み合わせに対して、組み合わせに対応する損失額を乗算し、全ての組み合わせの乗算結果を加算し、所定の最大地動加速度の地震が発生した際に工場が受ける総損失額の推定値を算出する。
すなわち、操業度算出部１３は、以下の式により、工場の受ける総損失額Ｑ（円、本実施形態においては万円）を算出する。
Ｑ＝０×ＰＮＮ（損傷がない場合、損失額は０円）
＋２０００×ＰＮＹ（製造ラインにのみ損傷がある場合、損失額は２０００万円）
＋１５００×ＰＹＮ（電力設備にのみ損傷がある場合、損失額は１５００万円）
＋３５００×ＰＹＹ（双方に損傷がある場合、損失額は３５００万円）
＝４４８
Next, the operation
That is, the operation
Q = 0 × P NN (If there is no damage, the loss is 0 yen)
+ 2000 × P NY (If the production line is damaged, the loss is 20 million yen)
+ 1500 × P YN (If there is damage only to the power equipment, the loss is 15 million yen)
+ 3500 × P YY (If both are damaged, the loss is 35 million yen)
= 448
そして、操業度算出部１３は、工場全体の再構成費、すなわち製造ライン及び電力設備の再構成費の合計額として総再構成費Ｒ（円、本実施形態においては万円）を、製造ラインの再構成費と電力設備の再構成費とを加算して求める。
Ｒ＝２０００＋１５００＝３５００
次に、操業度算出部１３は、上述のようにして求めた総再構成費Ｒにより、先に求めた総損失額Ｑを除算することにより、最大予想損失率（ＰＭＬ）を算出する。
ＰＭＬ＝Ｑ／Ｒ
＝４４８／３５００
＝０．１２８
Then, the operation
R = 2000 + 1500 = 3500
Next, the operation
PML = Q / R
= 448/3500
= 0.128
次に、操業度算出部１３は、操業度ＯＰ（Ｏｐｅｒａｔｅ ｐｅｒｃｅｎｔａｇｅ）を以下の式により、すなわち１から最大予想損失率を減算することにより算出する。
ＯＰ＝１−ＰＭＬ
＝０．８７２
上記式により、最大予想損失率を用い、操業度を定義している。この最大予想損失率は経営的な地震等のハザードによるリスクを考慮しており、経営的な操業度を表現することができる。
すなわち、本実施形態においては、例えば地震が発生した際に、事業構成要素がどの程度の損失を受けるかを示す損失額の期待値である最大予想損失率を、通常の状態である１００％から減算することにより、操業可能な割合として操業度が定義できるとしている。
Next, the operation
OP = 1-PML
= 0.872
The above formula defines the operating rate using the maximum expected loss rate. This maximum expected loss rate takes into account risks due to hazards such as management earthquakes, and can express the level of management operation.
That is, in this embodiment, for example, when an earthquake occurs, the maximum expected loss rate, which is an expected value of the amount of loss indicating how much the business component receives loss, from 100% that is a normal state. By subtracting, the operation rate can be defined as the operable ratio.
また、ＰＭＴ算出部１２は、最大復旧期待時間（ＰＭＴ）の算出を以下の式により、すなわち基本情報データベース１８の基本データテーブルから各事業構成要素の復旧時間を読み出す。本実施形態においては、製造ラインにおける機械１の復旧時間が１０日、機械２の復旧時間が８日であり、電力設備の商用電源が３日、自家発電源が２０日である。
そして、ＰＭＴ算出部１２は、損傷確率算出部１１の求めた各事業構成要素の損傷確率を、それぞれ対応する事業構成要素の復旧時間に乗算し、各事業構成要素の復旧時間を算出する。
ここで、損傷確率算出部１１により読み出された損傷確率において、機械１の損傷確率が０．２４０２０であり、機械２の損傷確率が０．０７４４であり、商用電源の損傷確率が０．００７４０であり、自家発電源の損傷確率が０．２２８６である。
The
Then, the
Here, in the damage probability read by the damage
すなわち、ＰＭＴ算出部１２は、以下の式により、機械１、機械２、商用電源及び自家発電源の最大復旧期待時間ＰＭＴ１（機械１の最大復旧期待時間）、ＰＭＴ２（機械２の最大復旧期待時間）、ＰＭＴ３（商用電源の最大復旧期待時間）、ＰＭＴ４（自家発電源の最大復旧期待時間）を算出する。
ＰＭＴ１＝１０×０．２４０２０＝２．４０２００
ＰＭＴ２＝８×０．０７４４＝０．５９５２０
ＰＭＴ３＝３×０．００７４０＝０．０２２２０
ＰＭＴ４＝２０×０．２２８６＝４０５７２００
That is, the
PMT 1 = 10 × 0.24020 = 2.40200
PMT 2 = 8 × 0.0744 = 0.59520
PMT 3 = 3 × 0.00740 = 0.02220
PMT 4 = 20 × 0.2286 = 4057200
次に、以下に示すように、復旧曲線生成部１４は、上述した最大復旧期待時間ＰＭＬと操業度ＯＰとにより、復旧曲線を生成する。
図９は、モデルデータベース１９に記憶されているモデルテーブルの構成を示す図である。このモデルテーブルは、復旧曲線生成部１４によりモデル設定時に予め生成され、事業構成要素の単位としてタイプが示され、モデルとしての工場がコンポーネントの製造ライン及び電力設備から構成され、コンポーネントの製造ラインがサブコンポーネントの機械１及び機械２から構成され、コンポーネントの電力設備がサブコンポーネントの商用電源及び自家発電源から構成されていることを示している。また、それぞれのタイプが分解可能、すなわちモデルが下位のコンポーネント、あるいはサブコンポーネントで構成されているか否かを示す事象が示されている。
Next, as shown below, the recovery
FIG. 9 is a diagram showing the configuration of the model table stored in the
また、ＰＭＴ算出部１２は、基本情報データベース１８から、事業構成要素毎の復旧時間を読み出し、モデルテーブルに対して事業構成要素毎に書き込み、記憶させる。操業度算出部１３は、事業構成要素毎の再構成価格を基本情報データベース１８から読み出し、モデルテーブルに対して事業構成要素毎に書き込み、記憶させる。損傷確率算出部１１は、算出したサブ事業構成要素（機械１、機械２、商用電源及び自家発電源）毎の損傷確率、事業構成要素（製造ライン及び電力設備）毎の損傷確率、工場の損傷確率（事業構成要素の損傷率の合計）を、モデルテーブルに書き込み、記憶させる。また、損傷確率算出部１１は、損傷確率の算出に用いたサブ事業構成要素の最大地動加速度の中央値、標準偏差（最大地動加速度の対数標準偏差）、最大地動加速度、及びフラジリティ曲線の形態を、サブ事業構成要素毎にモデルテーブルに書き込み、記憶させる。以下、損傷確率算出部１１、操業度算出部１３及び復旧曲線生成部１４の各部は、モデルテーブルから各事業構成要素の再構成費、損傷確率、ＰＭＴを読み出して用いる。
Further, the
図１０及び図１１を用いて、復旧曲線生成の処理を説明する。図１０は、本実施形態による復旧曲線作成システムの復旧曲線の生成の動作例を示すフローチャートである。図１１は、図１０のフローチャートにより作成される復旧曲線であり、縦軸が操業度ＯＰであり、横軸が経過日数（日）である。
また、シミュレーションが開始された際、復旧曲線生成部１４は、図１１に示す復旧曲線の画像を表示装置１６に表示する。ここで、復旧曲線生成部１４は、ハザードが発生する経過日数を「０」とし、ハザードが発生する経過日数を前日として「−１」とし、ハザードが発生する「−１」日から発生した当日の「０」日までの操業度ＯＰを１として、グラフに表示する。
ＰＭＴ算出部１２は、すでに説明したように、機械１の最大復旧期待時間ＰＭＴ１、機械２の最大復旧期待時間ＰＭＴ２、商用電源の最大復旧期待時間ＰＭＴ３及び自家発電源の最大復旧期待時間ＰＭＴ４を算出し（ステップＳ１）、算出した最大復旧期待時間ＰＭＴ１からＰＭＴ４を復旧曲線生成部１４へ、それぞれの要素識別情報とともに、この要素識別情報に対応付けて出力する。
The recovery curve generation process will be described with reference to FIGS. 10 and 11. FIG. 10 is a flowchart showing an operation example of generating a recovery curve of the recovery curve creating system according to the present embodiment. FIG. 11 is a recovery curve created according to the flowchart of FIG.
Further, when the simulation is started, the recovery
ここで、復旧曲線生成部１４は、自身の内部に設けられた記憶部に最大復旧期待時間ＰＭＴ１からＰＭＴ４の各々を、それぞれ要素識別情報に対応付けて書き込み、記憶する。このとき、復旧曲線生成部１４は、最大復旧期待時間ＰＭＴ１からＰＭＴ４の各々の小数点以下を切り上げ、整数値とする。例えば、最大復旧期待時間ＰＭＴ１は３日、最大復旧期待時間ＰＭＴ２は１日、最大復旧期待時間ＰＭＴ３は１日、最大復旧期待時間ＰＭＴ１は５日とされ、記憶されることになる。
また、このとき、復旧曲線生成部１４は、最大復旧期待時間ＰＭＴ１からＰＭＴ４の各々を比較し、最大の経過日数である最大経過日数を求める。
Here, the recovery
At this time, the recovery
次に、復旧曲線生成部１４は、内部記憶部にある経過日数レジスタを初期化、すなわち経過日数レジスタに記憶されている経過日数を「０」とする（ステップＳ２）。
Next, the recovery
そして、復旧曲線生成部１４は、内部レジスタに記憶されている経過日数と、最大復旧期待時間ＰＭＴ１からＰＭＴ４における最大値である最大経過日数とを比較する（ステップＳ３）
このとき、復旧曲線生成部１４は、経過日数が最大経過日数を超えている場合、復旧曲線の生成を終了し、経過日数が最大経過日数を超えていない場合、処理をステップＳ４へ進める。したがって、本実施形態の場合、最大経過日数が最大復旧期待時間ＰＭＴ４の５日であるため、経過日数が６日となると、完全に復旧１日後となるため処理を終了する。
Then, the recovery
At this time, the recovery
次に、復旧曲線生成部１４は、内部レジスタに記憶されている経過日数と、最大復旧期待時間ＰＭＴ１からＰＭＴ４の各々とを比較し、いずれかと一致するか否かの判定を行う（ステップＳ４）。
このとき、復旧曲線生成部１４は、経過日数が最大復旧期待時間ＰＭＴ１からＰＭＴ４のいずれかと一致した場合、処理をステップＳ５へ進め、経過日数が最大復旧期待時間ＰＭＴ１からＰＭＴ４のいずれとも一致しない場合、処理をステップＳ６へ進める。
Next, the recovery
At this time, when the elapsed days coincide with any of the maximum expected recovery times PMT 1 to PMT 4 , the recovery
次に、復旧曲線生成部１４は、現在の経過日数を超える最大復旧期待時間を有するサブ事業構成要素の要素識別情報を付加し、この要素識別情報に対応するサブ事業構成要素の損傷確率により、損傷確率ＰＡ及びＰＢを算出することを指示する制御信号を、損傷確率算出部１１に対して出力する。
このとき、復旧曲線生成部１４は、経過日数以下の最大復旧時間を有するサブ事業構成要素の損傷確率を「０」とする情報も、対応する要素識別情報とともに損傷確率算出部１１に出力する制御信号に付加する。
これにより、損傷確率算出部１１は、入力される制御信号に含まれる要素識別情報に基づき、いずれのサブ事業構成要素の損失確率を用い、またいずれのサブ事業構成要素の損失確率を「０」とするかの情報を用い、すでに説明した式により、損傷確率ＰＡ及びＰＢを算出する（ステップＳ５）。そして、損傷確率算出部１１は、算出した損傷確率ＰＡ及びＰＢを、操業度算出部１３へ出力する。
Next, the recovery
At this time, the recovery
Thereby, the damage
次に、操業度算出部１３は、損傷確率算出部１１から入力される損傷確率ＰＡ及びＰＢに基づき、すでに説明した確率算出式により、製造ライン及び電力設備の双方に損傷の無い確率ＰＮＮ、電力設備のみに損傷の有る確率ＰＮＹ、製造ラインのみに損傷の有る確率ＰＹＮ、製造ライン及び電力設備の双方に損傷の有る確率ＰＹＹを算出する。
そして、操業度算出部１３は、算出した確率ＰＮＮ、ＰＮＹ、ＰＹＮ及びＰＹＹと、製造ライン及び電力設備の各々の再構成費とにより、総損失額Ｑを算出する。
現在の経過日数における総損失額Ｑを算出すると、操業度算出部１３は、算出した総損失額Ｑを、すでに算出してある工場全体の事業構成要素に対する総再構成費Ｒにより除算し、現在の経過日数における最大予想損失率を算出する。
次に、操業度算出部１３は、算出した最大予想損失率を１から減算し、現在の経過日数における操業度ＯＰを算出し（ステップＳ６）、算出した操業度ＯＰを復旧曲線生成部１４に対して出力する。
Next, based on the damage probabilities P A and P B input from the damage
Then, the operation
When calculating the total loss Q for the current number of days, the operation
Next, the operation
操業度ＯＰが供給されると、復旧曲線生成部１４は、表示装置１６の表示画面に表示している図１１に示すグラフにおいて、現在の経過日数の座標点に、供給された操業度ＯＰをプロットして表示する（ステップＳ７）。
When the operation level OP is supplied, the recovery
次に、復旧曲線生成部１４は、自身内部の記憶部にある経過日数レジスタに記憶されている経過日数をインクリメント（１を加算）し（ステップＳ８）、インクリメントし経過日数を経過日数レジスタに書き込み、記憶させる。
そして、復旧曲線生成部１４は、処理をステップＳ３へ進める。
Next, the recovery
Then, the recovery
次に、図１１を用いて、以下に、復旧曲線の生成の処理の実例を説明する。
ＢＩＡ実施者により入力装置１５からシミュレーション開始の制御信号が供給されると、復旧曲線生成部１４は、経過日数「−１」から経過日数「０」までの操業度ＯＰを１としてプロットする。
次に、復旧曲線生成部１４は、ハザード、例えば地震が発生した際の経過日数「０」における操業度ＯＰの計算を開始する制御信号を、損傷確率算出部１１、ＰＭＴ算出部１２及び操業度算出部１３に対して出力する。こにれにより、損傷確率算出部１１、ＰＭＴ算出部１２及び操業度算出部１３は、指定されたモデルの情報をモデルデータベース１０のモデルテーブルに書き込む。
損傷確率算出部１１は、製造ライン、電力設備、機械１、機械２、商用電源及び自家発電源の各々の損傷確率（ＰＡ、ＰＢ、Ｐ１、Ｐ２、Ｐ３及びＰ４）を算出し、モデルテーブルに書き込み、記憶させる。
ＰＭＴ算出部１２は、機械１、機械２、商用電源及び自家発電源の各々の最大復旧期待時間（ＰＭＴ１、ＰＭＴ２、ＰＭＴ３、ＰＭＴ４）を、それぞれ損傷確率と復旧日数とを乗算することにより算出し、モデルテーブルに書き込み、記憶させる。
Next, an example of the process of generating a recovery curve will be described below with reference to FIG.
When a simulation start control signal is supplied from the
Next, the recovery
The damage
The
このとき、復旧曲線生成部１４は、モデルテーブルにある最大復旧期待時間の中から最大の日数（本実施形態においては最大復旧期待時間ＰＭＴ４の５日）を抽出し、内部記憶部に最大経過日数として書き込んで記憶させる。
ここで、機械１、機械２、商用電源及び自家発電源の各々の最大復旧期待時間は、すでに説明したように、それぞれ３日、１日、１日、５日と求められる。
また、操業度算出部１３は、損傷確率算出部１１が算出した損傷確率により、それぞれの操業度ＯＰを算出する。
この経過日数「０」日における操業度ＯＰは、すでに算出したように、「０．８７２」である。
このため、復旧曲線生成部１４は、得られた操業度ＯＰ＝０．８７２を経過日数「０」にプロットし、操業度ＯＰが１から０．８７２へ低下してことを表示画面のグラフに表示する。また、復旧曲線生成部１４は、経過日数をインクリメントし、経過日数「１」とする。
At this time, the recovery
Here, the maximum expected recovery time of each of the
Further, the operation
The operation level OP on this elapsed day “0” is “0.872” as already calculated.
For this reason, the recovery
次に、経過日数「１」において、復旧曲線生成部１４は、経過日数が最大経過日数を超えたか否かの判定を行い、超えていないことを検出すると、経過日数と一致する最大復旧期待時間の事業構成要素をモデルテーブルにおいて検索する。
ここで、復旧曲線生成部１４は、経過日数と一致する最大復旧期待時間を有する事業構成要素として機械２と商用電源とを検出する。
検出後、復旧曲線生成部１４は、機械１と自家発電源とが経過日数を超えた事業構成要素とし、機械２と商用電源とが損傷確率を「０」とする事業構成要素とし、それぞれの構成識別情報を含む制御信号を、損傷確率算出部１１へ出力する。
そして、損傷確率算出部１１は、上述した制御信号が供給されると、機械２と商用電源との損傷確率を「０」とし、以下の式により、製造ラインの損傷確率ＰＡ及び電力設備の損傷確率ＰＢの算出を行う。
ＰＡ＝１−（１−Ｐ１）・（１−Ｐ２）
＝１−（１−０．２４０２０）・（１−０）
＝１−０．７５９８
＝０．２４０２
ＰＢ＝Ｐ３・Ｐ４
＝０．２２８６０×０
＝０
Next, in the elapsed days “1”, the recovery
Here, the recovery
After the detection, the recovery
Then, when the above-described control signal is supplied, the damage
P A = 1- (1-P 1 ) · (1-P 2 )
= 1- (1-0.24020). (1-0)
= 1-0.7598
= 0.2402
P B = P 3・ P 4
= 0.22860 × 0
= 0
そして、損傷確率算出部１１は、算出した経過日数「１」における損傷確率ＰＡ及びＰＢを、操業度算出部１３へ出力する。
操業度算出部１３は、損傷確率算出部１１から供給される損傷確率ＰＡ及びＰＢにより、製造ラインと電力設備とにおける損傷の有無の組み合わせの確率ＰＮＮ（製造ライン及び電力設備の双方に損傷無し）、ＰＮＹ（電力設備のみに損傷有り）、ＰＹＮ（製造ラインのみに損傷有り）、ＰＹＹ（製造ライン及び電力設備の双方に損傷有り）の算出を以下の式により行う。
ＰＮＮ＝（１−ＰＡ）・（１−ＰＢ）
＝（１−０．２４０２）×（１−０）
＝０．７５９８
ＰＮＹ＝（１−ＰＡ）・ＰＢ
＝（１−０．２４０２）×０
＝０
ＰＹＮ＝ＰＡ・（１−ＰＢ）
＝０．２４０２×（１−０）
＝０．２４０２
ＰＹＹ＝ＰＡ・ＰＢ
＝０．２４０２×０
＝０
Then, the damage
The operation
PNN = (1-P A ) · (1-P B )
= (1-0.2402) x (1-0)
= 0.7598
P NY = (1-P A ) · P B
= (1-0.2402) x 0
= 0
P YN = P A · (1 -P B)
= 0.2402 × (1-0)
= 0.2402
P YY = P A・ P B
= 0.2402 × 0
= 0
次に、操業度算出部１３は、上述のように求めた確率ＰＮＮ、ＰＮＹ、ＰＹＮ、ＰＹＹの各々を、以下の式に示すように、対応する事業構成要素の損失額に乗算し、乗算結果を加算して、工場の受ける総損失額Ｑを算出する。
Ｑ＝０×ＰＮＮ（損傷がない場合、損失額は０円）
＋２０００×ＰＮＹ（製造ラインにのみ損傷がある場合、損失額は２０００万円）
＋１５００×ＰＹＮ（電力設備にのみ損傷がある場合、損失額は１５００万円）
＋３５００×ＰＹＹ（双方に損傷がある場合、損失額は３５００万円）
＝０×０．７５９８＋２０００×０＋１５００×０．２４０２＋３５００×０
＝３６０．３
≒３６０
Next, the operation
Q = 0 × P NN (If there is no damage, the loss is 0 yen)
+ 2000 × P NY (If the production line is damaged, the loss is 20 million yen)
+ 1500 × P YN (If there is damage only to the power equipment, the loss is 15 million yen)
+ 3500 × P YY (If both are damaged, the loss is 35 million yen)
= 0 × 0.7598 + 2000 × 0 + 1500 × 0.2402 + 3500 × 0
= 360.3
≒ 360
そして、操業度算出部１３は、以下の式に示すように、算出した経過日数「１」における総損失額Ｑ＝３６０を、予め計算してある総再構成費Ｒ＝３５００により除算し、最大予想損失率ＰＭＬを算出する。
ＰＭＬ＝Ｑ／Ｒ
＝３６０／３５００
＝０．１０２８５７
最大予想損失率を算出すると、操業度算出部１３は、以下の式に示すように、１から最大予想損失率ＰＭＬ＝０．１０２８５７を減算し、操業度ＯＰを算出する。
ＯＰ＝１−０．１０２８５７
＝０．８９７１４３
Then, as shown in the following formula, the operation
PML = Q / R
= 360/3500
= 0.102857
When the maximum expected loss rate is calculated, the operation
OP = 1-0.1002857
= 0.897143
次に、復旧曲線生成部１４は、得られた操業度ＯＰ＝０．８９７１４３を経過日数「１」にプロットし、機械２と商用電源とが損傷から可動状態となり、操業度ＯＰが０．８７２から０．８９７１４３へ上昇したことを、表示装置６の表示画面のグラフに表示する。また、復旧曲線生成部１４は、経過日数をインクリメントし、経過日数「２」とする。
Next, the recovery
次に、経過日数「２」において、復旧曲線生成部１４は、経過日数が最大経過日数を超えたか否かの判定を行い、超えていないことを検出すると、経過日数と一致する最大復旧期待時間の事業構成要素をモデルテーブルにおいて検索する。
しかしながら、復旧曲線生成部１４は、モデルテーブルに２日の最大復旧期待時間がないため、経過日数「２」と一致する最大復旧期待時間を有する事業構成要素を抽出することができない。
このため、復旧曲線生成部１４は、新たに損傷状態から稼動状態に遷移した事業構成要素がないので、前回算出した操業度ＯＰから変化がないものとし、前回の操業度ＯＰを表示装置１６の表示画面のグラフに表示する。また、復旧曲線生成部１４は、経過日数をインクリメントし、経過日数「３」とする。
Next, in the elapsed days “2”, the recovery
However, the recovery
For this reason, since there is no business component which newly changed from the damage state to the operation state, the recovery
次に、経過日数「３」において、復旧曲線生成部１４は、経過日数が最大経過日数を超えたか否かの判定を行い、超えていないことを検出すると、経過日数と一致する最大復旧期待時間の事業構成要素をモデルテーブルにおいて検索する。
ここで、復旧曲線生成部１４は、経過日数と一致する最大復旧期待時間を有する事業構成要素として機械１を検出する。
検出後、復旧曲線生成部１４は、自家発電源が経過日数を超えた事業構成要素とし、機械１、機械２及び商用電源とが損傷確率を「０」とする事業構成要素とし、それぞれの構成識別情報を含む制御信号を、損傷確率算出部１１へ出力する。
そして、損傷確率算出部１１は、上述した制御信号が供給されると、機械１、機械２及び商用電源との損傷確率を「０」とし、以下の式により、製造ラインの損傷確率ＰＡ及び電力設備の損傷確率ＰＢの算出を行う。
ＰＡ＝１−（１−Ｐ１）・（１−Ｐ２）
＝１−（１−０）・（１−０）
＝１−０
＝０
ＰＢ＝Ｐ３・Ｐ４
＝０．２２８６０×０
＝０
Next, in the elapsed days “3”, the recovery
Here, the recovery
After the detection, the recovery
Then, damage
P A = 1- (1-P 1 ) · (1-P 2 )
= 1- (1-0). (1-0)
= 1-0
= 0
P B = P 3・ P 4
= 0.22860 × 0
= 0
そして、損傷確率算出部１１は、算出した経過日数「３」における損傷確率ＰＡ及びＰＢを、操業度算出部１３へ出力する。
操業度算出部１３は、損傷確率算出部１１から供給される損傷確率ＰＡ及びＰＢにより、製造ラインと電力設備とにおける損傷の有無の組み合わせの確率ＰＮＮ（製造ライン及び電力設備の双方に損傷無し）、ＰＮＹ（電力設備のみに損傷有り）、ＰＹＮ（製造ラインのみに損傷有り）、ＰＹＹ（製造ライン及び電力設備の双方に損傷有り）の算出を以下の式により行う。
ＰＮＮ＝（１−ＰＡ）・（１−ＰＢ）
＝（１−０）×（１−０）
＝１
ＰＮＹ＝（１−ＰＡ）・ＰＢ
＝（１−０）×０
＝０
ＰＹＮ＝ＰＡ・（１−ＰＢ）
＝０×（１−０）
＝０
ＰＹＹ＝ＰＡ・ＰＢ
＝０×０
＝０
Then, the damage
The operation
PNN = (1-P A ) · (1-P B )
= (1-0) × (1-0)
= 1
P NY = (1-P A ) · P B
= (1-0) × 0
= 0
P YN = P A · (1 -P B)
= 0 × (1-0)
= 0
P YY = P A・ P B
= 0x0
= 0
次に、操業度算出部１３は、上述のように求めた確率ＰＮＮ、ＰＮＹ、ＰＹＮ、ＰＹＹの各々を、以下の式に示すように、対応する事業構成要素の損失額に乗算し、乗算結果を加算して、工場の受ける総損失額Ｑを算出する。
Ｑ＝０×ＰＮＮ（損傷がない場合、損失額は０円）
＋２０００×ＰＮＹ（製造ラインにのみ損傷がある場合、損失額は２０００万円）
＋１５００×ＰＹＮ（電力設備にのみ損傷がある場合、損失額は１５００万円）
＋３５００×ＰＹＹ（双方に損傷がある場合、損失額は３５００万円）
＝０×１＋２０００×０＋１５００×０＋３５００×０
＝０
Next, the operation
Q = 0 × P NN (If there is no damage, the loss is 0 yen)
+ 2000 × P NY (If the production line is damaged, the loss is 20 million yen)
+ 1500 × P YN (If there is damage only to the power equipment, the loss is 15 million yen)
+ 3500 × P YY (If both are damaged, the loss is 35 million yen)
= 0x1 + 2000x0 + 1500x0 + 3500x0
= 0
そして、操業度算出部１３は、以下の式に示すように、算出した経過日数「３」における総損失額Ｑ＝０を、予め計算してある総再構成費Ｒ＝３５００により除算し、最大予想損失率ＰＭＬを算出する。
ＰＭＬ＝Ｑ／Ｒ
＝０／３５００
＝０
最大予想損失率を算出すると、操業度算出部１３は、以下の式に示すように、１から最大予想損失率ＰＭＬ＝０を減算し、操業度ＯＰを算出する。
ＯＰ＝１−０
＝１
Then, as shown in the following equation, the operation
PML = Q / R
= 0/3500
= 0
When the maximum expected loss rate is calculated, the operation
OP = 1-0
= 1
次に、復旧曲線生成部１４は、得られた操業度ＯＰ＝１を経過日数「３」にプロットし、機械１、機械２及び商用電源が損傷状態から可動状態となり、操業度ＯＰが０．８９７１４３から１のハザード前の操業度に戻ったことを、表示装置１６の表示画面のグラフに表示する。また、復旧曲線生成部１４は、経過日数をインクリメントし、経過日数「４」とする。
Next, the recovery
次に、経過日数「４」において、復旧曲線生成部１４は、経過日数が最大経過日数を超えたか否かの判定を行い、超えていないことを検出すると、経過日数と一致する最大復旧期待時間の事業構成要素をモデルテーブルにおいて検索する。
しかしながら、復旧曲線生成部１４は、モデルテーブルに２日の最大復旧期待時間がないため、経過日数「４」と一致する最大復旧期待時間を有する事業構成要素を抽出することができない。
このため、復旧曲線生成部１４は、新たに損傷状態から稼動状態に遷移した事業構成要素がないので、前回算出した操業度ＯＰから変化がないものとし、前回の操業度ＯＰを表示装置１６の表示画面のグラフに表示する。また、復旧曲線生成部１４は、経過日数をインクリメントし、経過日数「５」とする。
Next, in the elapsed days “4”, the recovery
However, the recovery
For this reason, since there is no business component which newly changed from the damage state to the operation state, the recovery
次に、経過日数「５」において、復旧曲線生成部１４は、経過日数が最大経過日数を超えたか否かの判定を行い、超えていないことを検出すると、経過日数と一致する最大復旧期待時間の事業構成要素をモデルテーブルにおいて検索する。
そして、復旧曲線生成部１４は、モデルテーブルに５日の最大復旧期待時間を有する事業構成要素として自家発電源があるため、経過日数「５」と一致する最大復旧期待時間を有する事業構成要素として自家発電源を抽出する。
しかしながら、復旧曲線生成部１４は、すでに操業度ＯＰが「１」となり、通常に工場が稼動している状態にあるため、前回算出した操業度ＯＰから変化がないものとし、前回の操業度ＯＰを表示装置１６の表示画面のグラフに表示する。また、復旧曲線生成部１４は、経過日数をインクリメントし、経過日数「６」とする。
Next, in the elapsed days “5”, the recovery
Then, since the recovery
However, the recovery
次に、経過日数「６」において、復旧曲線生成部１４は、経過日数が最大経過日数を超えたか否かの判定を行い、超えてることを検出すると、復旧曲線の生成処理を終了する。
上述したように、復旧曲線生成部１４は、図１１における工場のリソースの復旧曲線を生成する。
Next, in the elapsed days “6”, the recovery
As described above, the recovery
＜復旧曲線の生成−サービス（商品の生産）＞
次に、本実施形態による復旧曲線生成システムが行う商品Ａ及び商品Ｂの生産に対する操業度の算出について説明する。
以下の説明において、例えば、商品Ａは１日の売上げが６００万円で、製造ライン及び電力設備の双方が動作していなければ、生産ができず、一方、商品Ｂは一日の売上げが３００万円で、電力設備が動作していれば、生産ができるとする。
商品Ａ及び商品Ｂにおいても、各事業構成要素の損失確率の計算は、すでに説明したリソースの場合と、使用する式も数値も同様である。また、ＰＭＴ算出部１２が行う各事業構成要素に対する最大復旧期待時間の計算も同様である。
しかしながら、操業度計算部１３における総損失額Ｑを求める際、確率ＰＮＮ、ＰＮＹ、ＰＹＮ、ＰＹＹの算出に用いる損傷確率ＰＡ及びＰＢにおいて、商品Ａの場合は製造ライン及び電力設備の双方が稼動状態でなければ生産できないため、損傷確率ＰＡ及びＰＢともにリソースの計算に用いた数値と同様である。一方、商品Ｂの場合、電力設備さえあれば生産が可能であるため、製造ラインの損傷確率ＰＡを「０」とし、電力設備の損傷確率ＰＢはリソースの計算に用いた電力設備の損傷確率ＰＢと同一の数値である。
<Generation of recovery curve-service (product production)>
Next, calculation of the operation rate for the production of the products A and B performed by the recovery curve generation system according to the present embodiment will be described.
In the following description, for example, the product A has a daily sales of 6 million yen and cannot be produced unless both the production line and the power facility are in operation, while the product B has a daily sales of 300 If the power equipment is operating for 10,000 yen, it can be produced.
Also in the product A and the product B, the calculation of the loss probability of each business component is the same in the formula and numerical value used as in the case of the resource already described. The calculation of the maximum expected recovery time for each business component performed by the
However, when determining the total loss Q in the operating
また、図１０に示すフローチャートに従い、リソースの場合と同様に、商品Ａ及び商品Ｂの各々の復旧曲線を算出する。
ここで、商品の復旧曲線がリソースの復旧曲線の生成と異なる点は、総損失額Ｑの算出に用いる損失額が１日の売上額となっていることである。
以下に、操業度算出部１３の商品Ａ及び商品Ｂの操業度の算出処理について説明する。このとき、損傷確率算出部１１は、操業度算出部１３に対し、算出した損傷確率ＰＡ（＝０．２９６７３）及びＰＢ（＝０．００１６９）を出力する。
また、操業度算出部１３は、商品Ａが製造ライン及び電源設備の双方が可動状態でないと生産できないため、損傷確率ＰＡ及びＰＢをそのまま用い、商品Ｂが電源設備のみで生産できるため、損傷率ＰＡを「０」とし、損傷率ＰＢをそのまま用いるとの制御情報を、復旧曲線生成部１４により予め設定されている。この制御情報は、ＢＩＡ実施者が入力装置１５から入力する、復旧曲線を生成するための設定データである。
Further, according to the flowchart shown in FIG. 10, the recovery curves of the products A and B are calculated in the same manner as in the case of resources.
Here, the product recovery curve differs from the generation of the resource recovery curve in that the loss used for calculating the total loss Q is the daily sales.
Below, the calculation process of the operation degree of the goods A and the goods B of the operation
Moreover, the operating
まず、操業度算出部１３は、以下の様に、商品Ａの操業度ＯＰＡを、リソースで用いた確率算出式を用いて算出する。
すでに述べたように、商品Ａは、リソースと同様の製造ラインの損傷確率ＰＡと電力設備の損傷確率ＰＢとを用いて、製造ラインが損傷の有無と電力設備の損傷の有無との組み合わせの確率ＰＮＮ、ＰＮＹ、ＰＹＮ、ＰＹＹを算出するため、得られる数値は以下に示すようにリソースの場合と同様である。
ＰＮＮ＝０．７０２０８（製造ライン及び電力設備の双方に損傷が無い場合）
ＰＮＹ＝０．００１８９（電力設備のみに損傷がある場合）
ＰＹＮ＝０．２９６２２（製造ラインにのみ損傷がある場合）
ＰＹＹ＝０．０００５０（製造ライン及び電力設備の双方に損傷がある場合）
First, the operation
As described above, the commodity A uses the same damage probability P A of the production line as the resource and the damage probability P B of the power equipment to combine the presence / absence of the damage on the production line and the presence / absence of the damage on the power equipment. In order to calculate the probabilities P NN , P NY , P YN , and P YY , the obtained numerical values are the same as in the case of resources as shown below.
P NN = 0.70208 (when there is no damage to both the production line and the power equipment)
P NY = 0.00189 (when only the power equipment is damaged)
P YN = 0.29622 (when only the production line is damaged)
P YY = 0.00050 (when both the production line and power equipment are damaged)
次に、操業度算出部１３は、予め基本情報データベース１８に設定されている、商品Ａの１日の売上げ金額を読み出す。ここで、商品Ａの１日の売上げ金額は、損失額として商品Ａを識別する商品識別情報とともに対応付けて記憶されている。操業度算出部１３は、商品Ａの損失額を、商品Ａの商品識別情報により、基本情報データベース１８の基本データテーブルから読み出す。
本実施形態の場合、例えば、商品Ａの損失額が６００万円（１日）である。したがって、再構成費も商品Ａの損失額と同一の６００万円となる。
Next, the operation
In the case of the present embodiment, for example, the loss amount of the product A is 6 million yen (1 day). Therefore, the reconstruction cost is 6 million yen, which is the same as the loss amount of the product A.
次に、操業度算出部１３は、製造ライン及び電力設備の損傷の有無によるそれぞれの組み合わせに対して、商品Ａの損失額を乗算し、全ての組み合わせの乗算結果を加算し、所定の最大地動加速度の地震が発生した際に工場が受ける商品Ａの総損失額の推定値を算出する。
すなわち、操業度算出部１３は、以下の式により、工場の受ける商品Ａの総損失額Ｑ（円、本実施形態においては万円）を算出する。
Ｑ＝０×ＰＮＮ（損傷がない場合、損失額は０円）
＋６００×ＰＮＹ（製造ラインにのみ損傷がある場合、損失額は６００万円）
＋６００×ＰＹＮ（電力設備にのみ損傷がある場合、損失額は６００万円）
＋６００×ＰＹＹ（双方に損傷がある場合、損失額は６００万円）
＝１７８．７５
Next, the operation
That is, the operation
Q = 0 × P NN (If there is no damage, the loss is 0 yen)
+ 600 × P NY (If the production line is damaged, the loss is 6 million yen)
+ 600 × P YN (If there is damage only to the power equipment, the loss is 6 million yen)
+ 600 × P YY (If both are damaged, the loss is 6 million yen)
= 178.75
そして、操業度算出部１３は、以下の式に示すように、総損失額Ｑ＝１７８．７５（万円）を、予め設定されている再構成費Ｒ＝６００（万円）により除算し、最大予想損失率ＰＭＬを算出する。
ＰＭＬ＝Ｑ／Ｒ
＝１７８．５／６００
＝０．２９７９１９
最大予想損失率を算出すると、操業度算出部１３は、以下の式に示すように、１から最大予想損失率ＰＭＬ＝０．２９７９１９を減算し、操業度ＯＰＡを算出する。
ＯＰＡ＝１−０．２９７９１９
＝０．７０２０８１
Then, as shown in the following formula, the operation
PML = Q / R
= 178.5 / 600
= 0.297919
When the maximum expected loss rate is calculated, the operation
OPA = 1-0.297919
= 0.702081
次に、操業度算出部１３は、以下の様に、商品Ｂの操業度ＯＰＢを、リソースで用いた確率算出式を用いて算出する。
すでに述べたように、商品Ａは、製造ラインの損傷確率ＰＡを「０」とし、リソースと同様の電力設備の損傷確率ＰＢ（＝０．００１６９）とを用いて、製造ラインが損傷の有無と電力設備の損傷の有無との組み合わせの確率ＰＮＮ、ＰＮＹ、ＰＹＮ、ＰＹＹを算出する。
操業度算出部１３は、製造ラインが損傷が無く、電力設備も損傷が無い場合、以下の確率算出式により、いずれにも損傷が無い確率ＰＮＮを算出する。
ＰＮＮ＝（１−ＰＡ）・（１−ＰＢ）
＝（１−０）×（１−０．００１６９）
＝０．９９８３１
また、操業度算出部１３は、製造ラインが損傷が無く、電力設備も損傷がある場合、以下の確率算出式により、いずれにも電力設備のみに損傷がある確率ＰＮＹを算出する。
ＰＮＹ＝（１−ＰＡ）・ＰＢ
＝（１−０）×０．００１６９
＝０．００１６９
また、操業度算出部１３は、製造ラインが損傷が有り、電力設備に損傷が無い場合、以下の確率算出式により、製造ラインにのみ損傷がある確率ＰＹＮを算出する。
ＰＹＮ＝ＰＡ・（１−ＰＢ）
＝０×（１−０．００１６９）
＝０
また、操業度算出部１３は、製造ライン及び電力設備の双方に損傷がある場合、以下の確率算出式により、いずれにも損傷がある確率ＰＹＹを算出する。
ＰＹＹ＝ＰＡ・ＰＢ
＝０×０．００１６９
＝０
Next, the operation
As already mentioned, product A is a damage probability P A of the production line is "0", using a damage probability P B of similar power equipment and resources (= 0.00169), the production line is damaged The probabilities P NN , P NY , P YN , and P YY of the combination of presence / absence and presence / absence of damage to the power equipment are calculated.
When the production line is not damaged and the power equipment is not damaged, the operation
PNN = (1-P A ) · (1-P B )
= (1-0) × (1-0.00169)
= 0.99831
Further, when the production line is not damaged and the power equipment is also damaged, the operation
P NY = (1-P A ) · P B
= (1-0) x 0.00169
= 0.00169
Further, when the production line is damaged and the power equipment is not damaged, the operation
P YN = P A · (1 -P B)
= 0x (1-0.00169)
= 0
In addition, when both the production line and the power equipment are damaged, the operation
P YY = P A・ P B
= 0x0.00169
= 0
次に、操業度算出部１３は、予め基本情報データベース１８に設定されている、商品Ｂの１日の売上げ金額を読み出す。ここで、商品Ｂの１日の売上げ金額は、損失額として商品Ｂを識別する商品識別情報とともに対応付けて記憶されている。操業度算出部１３は、商品Ｂの損失額を、商品Ｂの商品識別情報により、基本情報データベース１８の基本データテーブルから読み出す。
本実施形態の場合、例えば、商品Ｂの損失額が３００万円（１日）である。したがって、再構成費も商品Ｂの損失額と同一の３００万円となる。
Next, the operation
In the case of the present embodiment, for example, the loss amount of the product B is 3 million yen (1 day). Therefore, the reconstruction cost is 3 million yen which is the same as the loss amount of the product B.
次に、操業度算出部１３は、製造ライン及び電力設備の損傷の有無によるそれぞれの組み合わせに対して、商品Ｂの損失額を乗算し、全ての組み合わせの乗算結果を加算し、所定の最大地動加速度の地震が発生した際に工場が受ける商品Ｂの総損失額の推定値を算出する。
すなわち、操業度算出部１３は、以下の式により、工場の受ける商品Ｂの総損失額Ｑ（円、本実施形態においては万円）を算出する。
Ｑ＝０×ＰＮＮ（損傷がない場合、損失額は０円）
＋３００×ＰＮＹ（製造ラインにのみ損傷がある場合、損失額は３００万円）
＋３００×ＰＹＮ（電力設備にのみ損傷がある場合、損失額は３００万円）
＋３００×ＰＹＹ（双方に損傷がある場合、損失額は３００万円）
＝０．５１
Next, the operation
That is, the operation
Q = 0 × P NN (If there is no damage, the loss is 0 yen)
+ 300 × P NY (If only the production line is damaged, the loss is 3 million yen)
+ 300 × P YN (If there is damage only to the power equipment, the loss is 3 million yen)
+ 300 × P YY (If both are damaged, the loss is 3 million yen)
= 0.51
そして、操業度算出部１３は、以下の式に示すように、総損失額Ｑ＝０．５１（万円）を、予め設定されている再構成費Ｒ＝３００（万円）により除算し、最大予想損失率ＰＭＬを算出する。
ＰＭＬ＝Ｑ／Ｒ
＝０．５１／３００
＝０．００１６２
最大予想損失率を算出すると、操業度算出部１３は、以下の式に示すように、１から最大予想損失率ＰＭＬ＝０．００１６２を減算し、操業度ＯＰＢを算出する。
ＯＰＢ＝１−０．００１６２
＝０．９９８３０８
復旧曲線生成部１４は、リソースの場合と同様に、表示装置１６の表示画面において、図１１に示すグラフに供給される操業度ＯＰＡ及びＯＰＢを算出した経過日数に対応する位置にプロットする。
Then, as shown in the following formula, the operation
PML = Q / R
= 0.51 / 300
= 0.00162
When the maximum expected loss rate is calculated, the operation
OPB = 1-0.00162
= 0.998308
The recovery
図１１に示す商品Ａ及び商品Ｂの操業度の経過日数毎における算出処理は、すでに述べたリソースの場合と同様であるため、詳細な説明を省略する。
また、操業度算出部１３は、商品Ａ及び商品Ｂの各々の損失額を、経過日数毎に積算し、操業度ＯＰＡ、ＯＰＢが通常に戻るまでの積算値を求める。
これにより、復旧曲線生成部１４は、製品単位に、それぞれの製品の生産に必要な事業構成要素の損傷状態により、生産が行えるまでの売上げに対する損失を、上述した各商品の損失額の積算値として算出することができる。ここで、復旧曲線生成部１４は、表示装置１６の表示画面における図１１のグラフの経過日数の部分に、算出した積算値を表示する。
また、復旧曲線生成部１４は、リソースで求めた、ハザードが発生した当日、すなわち経過日数が「０」の際の損失を、商品の積算値に加算することにより、発生したハザードによるリソース自体（事業構成要素）の損失額と、事業構成要素の損傷による商品の生産が行えないことによる売上げの損失額との総計を容易に算出することができる。
The calculation process for each elapsed day of the operation degrees of the products A and B shown in FIG. 11 is the same as that of the resource described above, and thus detailed description is omitted.
Further, the operation
As a result, the recovery
In addition, the recovery
＜損傷確率のステップ変化及びリニア変化＞
上述してきた説明において、損傷確率の変化は、ハザードが発生した時点に設定された損傷確率が最大復旧期待時間経過後に「０」となるステップでの変化になっている。
しかしながら、本実施形態に対し、最大復旧期待時間を損失率が「０」にいきなり変化する唯一のトリガとするのではなく、損失確率が設定値から「０」に経過日数により段階的に変化する、すなわちリニアに変化させ、経過日数毎に操業度ＯＰを算出する構成としても良い。
<Step change and linear change of damage probability>
In the above description, the change in the damage probability is a change in a step in which the damage probability set at the time of occurrence of the hazard becomes “0” after the maximum expected recovery time has elapsed.
However, with respect to the present embodiment, the maximum recovery expected time is not the only trigger that suddenly changes the loss rate to “0”, but the loss probability changes stepwise from the set value to “0” depending on the number of days elapsed. That is, it is good also as a structure which changes linearly and calculates the operation degree OP for every elapsed days.
上述したように、損傷率をステップ変化させた場合、事業構成要素の損傷確率は経過日数毎に変化はなく、稼動状態となる経過日数の日に、予め設定されている損傷確率の数値から０に変化し、このタイミングにおいて操業度ＯＰの算出を行う。
一方、損傷率をリニア変化させた場合、事業構成要素の損傷確率を、ハザードが発生した日の最高の損傷確率を基本確率として、経過日数ごとに以下の式により変化させて、経過日数毎の操業度の算出を行う。
ここで、損傷確率（ステップ）のリニア変化における損傷確率（経過日数毎）は、損傷確率算出部１１が以下の式により求める。
損傷確率（経過日数毎）＝基本損傷確率（１−（経過日数／ＰＭＴ））
すなわち損傷を受けてから徐々に回復する場合、例えば１００台の装置があり、５０台が損傷して、日数が経過する毎に、この損傷した５０台において数台ずつ可動状態となり、徐々に損傷確率が変化することを想定している。
As described above, when the damage rate is changed in steps, the damage probability of the business component does not change every elapsed day, and it is 0 from the numerical value of the preset damage probability on the day of the elapsed day to be in an operating state. The operating degree OP is calculated at this timing.
On the other hand, when the damage rate is changed linearly, the damage probability of the business components is changed according to the following formula for each elapsed day using the highest probability of damage on the day when the hazard occurs as the basic probability. Calculate the operating rate.
Here, the damage probability (for each elapsed day) in the linear change of the damage probability (step) is obtained by the damage
Damage probability (Elapsed days) = Basic damage probability (1- (Elapsed days / PMT))
That is, when recovering gradually after being damaged, there are 100 devices, for example, 50 units are damaged, and every time the number of days elapses, several units of these 50 units become movable and gradually damaged. It is assumed that the probability changes.
次に、図１２は、本実施形態による復旧曲線作成システムの復旧曲線の生成の動作例を示すフローチャートである。この図１２は、損傷確率をステップ変化とするか、あるいはリニア変化とするかを、事業構成要素毎に設定する機能を、損傷確率算出部１１が有する場合の復旧曲線作成システムの動作例を示す。
予めＢＩＡ実施者は、入力装置１５から、事業構成要素毎に、操業度ＯＰを算出する際に、損傷確率をステップ変化させるかリニア変化させるかの設定を行う。
これにより、復旧曲線生成部１４は、モデル設定時に、図９に示すモデルテーブルにおける途中復旧率の欄に、損傷確率の変化をステップかリニアかの情報が記憶される。
Next, FIG. 12 is a flowchart illustrating an operation example of generating a recovery curve of the recovery curve creation system according to the present embodiment. FIG. 12 shows an operation example of the recovery curve creation system in the case where the damage
The BIA practitioner sets in advance whether the damage probability is changed stepwise or linearly when calculating the operation level OP for each business component from the
Thereby, the restoration
図１２のフローチャートにおいて、図１０のフローチャートと同様のステップに、同一のステップ番号を付与している。
ステップＳ２までは、図１０のフローチャートと同様の動作であるため、説明を省略する。そして、ステップＳ２において、復旧曲線生成部１４は、経過日数Ｄを「０」にリセットした後、処理をステップＳ１０へ進める。
復旧曲線生成部１４は、モデルデータベース１９におけるモデルテーブルにおいて、各事業構成要素の途中復旧率の欄のデータを読み込み、各事業構成要素が損傷確率をステップあるいはリニアかであるかの判定を行い、リニアに設定されている事業構成要素の有無の検出を行う（ステップＳ１０）。
そして、復旧曲線生成部１４は、復旧曲線を生成するモデルのモデルテーブルにおける全ての事業構成要素がステップである場合、すなわち損傷確率をリニア処理して用いると設定された事業構成要素が無い場合、処理をステップＳ３へ進める。ここで、ステップＳ３に進んだ場合、後の処理は図１０のフローチャートの動作と同様のため、これ以降のステップＳ４からステップＳ８までの処理についての説明を省略する。
In the flowchart of FIG. 12, the same step numbers are assigned to the same steps as those in the flowchart of FIG.
The operations up to step S2 are the same as those in the flowchart of FIG. In step S2, the recovery
The recovery
Then, when all the business components in the model table of the model that generates the recovery curve are steps, that is, when there is no business component that is set to use the damage probability after linear processing, The process proceeds to step S3. Here, when the process proceeds to step S3, the subsequent process is the same as the operation of the flowchart of FIG.
一方、ステップＳ１０において、復旧曲線生成部１４は、復旧曲線を生成するモデルのモデルテーブルにおいて、リニアとして設定された事業構成要素が存在する場合、すなわち損傷確率をリニア処理して用いると設定された事業構成要素がある場合、処理をステップＳ１１へ進める。
On the other hand, in Step S10, the recovery
次に、復旧曲線生成部１４は、損傷確率算出部１１に対し、リニアと設定された事業構成要素の要素識別情報を付加し、この要素識別情報に対応する事業構成要素の損傷確率をリニア計算することを指示するリニア計算指示信号を出力する。
そして、損傷確率算出部１１は、上述したリニア計算指示信号が供給されると、付加されている要素識別情報に対応する事業構成要素の損傷確率（経過日数毎）の算出を行う（ステップＳ１１）。
ここで、損傷確率算出部１１は、例えば、図９に示すように、事業構成要素として機械１がリニアと設定されている場合、機械１の損傷確率から損傷確率（経過日数毎）の算出を行う。損傷確率算出部１１は、機械１の最大復旧期待時間が３日のため、以下のように各経過日数の損傷確率（経過日数毎）を算出する。
経過日数Ｄ＝０：
損傷確率（経過日数毎）＝０．２４０２０×（１−（０／３））
＝０．２４０２０
経過日数Ｄ＝１：
損傷確率（経過日数毎）＝０．２４０２０×（１−（１／３））
＝０．１６０１３
経過日数Ｄ＝２：
損傷確率（経過日数毎）＝０．２４０２０×（１−（２／３））
＝０．０８００６
経過日数Ｄ＝３：
損傷確率（経過日数毎）＝０．２４０２０×（１−（３／３））
＝０
Next, the recovery
And the damage
Here, for example, as shown in FIG. 9, the damage
Elapsed days D = 0:
Damage probability (every elapsed days) = 0.24020 × (1− (0/3))
= 0.24020
Elapsed days D = 1:
Damage probability (every elapsed days) = 0.24020 × (1− (1/3))
= 0.16013
Elapsed days D = 2:
Damage probability (every elapsed days) = 0.24020 × (1− (2/3))
= 0.08006
Elapsed days D = 3:
Damage probability (every elapsed days) = 0.24020 × (1− (3/3))
= 0
そして、損傷確率算出部１１は、算出した経過日数毎の機械１の損傷確率を、モデルテーブルの機械１の要素識別情報に対応した損傷確率の欄に、経過日数毎の損傷確率（経過日数毎）を書き込み、記憶させるとともに、復旧曲線生成部１４に対してリニア処理のための損傷確率（経過日数毎）の算出が終了したことを通知する。
次に、復旧曲線生成部１４は、損傷確率算出部１１から損傷確率（経過日数毎）の算出が終了したことが通知されると、処理をステップＳ１３へ進める。
Then, the damage
Next, when the recovery
そして、復旧曲線生成部１４は、ステップＳ３と同様に、内部レジスタに記憶されている経過日数と、最大復旧期待時間ＰＭＴ１からＰＭＴ４における最大値である最大経過日数とを比較する（ステップＳ１３）。
このとき、復旧曲線生成部１４は、経過日数が最大経過日数を超えている場合、復旧曲線の生成を終了し、経過日数が最大経過日数を超えていない場合、処理をステップＳ４へ進める。したがって、本実施形態の場合、最大経過日数が最大復旧期待時間ＰＭＴ４の５日であるため、経過日数が６日となると、完全に復旧１日後となるため処理を終了する。
Then, similarly to step S3, the recovery
At this time, the recovery
次に、復旧曲線生成部１４は、現在の経過日数を超える最大復旧期待時間を有するサブ事業構成要素の要素識別情報と、送信するサブ事業構成要素の要素識別情報にリニアの損傷確率を使用することを示すフラグを付加し、この要素識別情報に対応するサブ事業構成要素の損傷確率により、損傷確率ＰＡ及びＰＢを算出することを指示する制御信号を、損傷確率算出部１１に対して出力する。
このとき、復旧曲線生成部１４は、経過日数以下の最大復旧時間を有するサブ事業構成要素の損傷確率を「０」とする情報も、対応する要素識別情報とともに損傷確率算出部１１に出力する制御信号に付加する。
Next, the recovery
At this time, the recovery
これにより、損傷確率算出部１１は、入力される制御信号に含まれる要素識別情報に基づき、いずれのサブ事業構成要素の損失確率を用い、またいずれのサブ事業構成要素の損失確率を「０」とするかの情報を用い、すでに説明した式により、損傷確率ＰＡ及びＰＢを算出する（ステップＳ１５）。このとき、損傷確率算出部１１は、リニアの損傷確率を使用することを示すフラグが付加された要素識別情報を検出すると、この要素識別番号に対応する事業構成要素の現在の経過日数に対応する損傷確率（経過日数毎）を、モデルテーブルからこの損傷確率（経過日数毎）を読み出す。そして、損傷確率算出部１１は、リニアと指定された事業構成要素の損傷確率として、モデルテーブルから読み出した損傷確率（経過日数毎）を用いて、損傷確率ＰＡ及びＰＢを算出する。
そして、損傷確率算出部１１は、算出した損傷確率ＰＡ及びＰＢを、操業度算出部１３へ出力する。
Thereby, the damage
Then, the damage
次に、操業度算出部１３は、損傷確率算出部１１から入力される損傷確率ＰＡ及びＰＢに基づき、すでに説明した確率算出式により、製造ライン及び電力設備の双方に損傷の無い確率ＰＮＮ、電力設備のみに損傷の有る確率ＰＮＹ、製造ラインのみに損傷の有る確率ＰＹＮ、製造ライン及び電力設備の双方に損傷の有る確率ＰＹＹを算出する。
そして、操業度算出部１３は、算出した確率ＰＮＮ、ＰＮＹ、ＰＹＮ及びＰＹＹと、製造ライン及び電力設備の各々の再構成費とにより、総損失額Ｑを算出する。
現在の経過日数における総損失額Ｑを算出すると、操業度算出部１３は、算出した総損失額Ｑを、すでに算出してある工場全体の事業構成要素に対する総再構成費Ｒにより除算し、現在の経過日数における最大予想損失率を算出する。
次に、操業度算出部１３は、算出した最大予想損失率を１から減算し、現在の経過日数における操業度ＯＰを算出し（ステップＳ１６）、算出した操業度ＯＰを復旧曲線生成部１４に対して出力する。
Next, based on the damage probabilities P A and P B input from the damage
Then, the operation
When calculating the total loss Q for the current number of days, the operation
Next, the operation
操業度ＯＰが供給されると、復旧曲線生成部１４は、表示装置１６の表示画面に表示している図１１に示すグラフにおいて、現在の経過日数の座標点に、供給された操業度ＯＰをプロットして表示する（ステップＳ１７）。
When the operation level OP is supplied, the recovery
次に、復旧曲線生成部１４は、内部記憶部にある経過日数レジスタに記憶されている経過日数をインクリメント（１を加算）し（ステップＳ１８）、シンクリメントし経過日数を経過日数レジスタに書き込み、記憶させる。
そして、復旧曲線生成部１４は、処理をステップＳ１３へ進める。
Next, the recovery
Then, the recovery
次に、図１３を用いて、以下に、損傷確率をリニア処理した損傷確率を用いて操業度を算出し、復旧曲線の生成処理を行う実例を説明する。ここで、図９のモデルテーブルに示すように、機械１の損傷確率をリニア処理するとして以下の説明を行う。
ＢＩＡ実施者により入力装置１５からシミュレーション開始の制御信号が供給されると、復旧曲線生成部１４は、経過日数「−１」から経過日数「０」までの操業度ＯＰを１としてプロットする。
次に、復旧曲線生成部１４は、ハザード、例えば地震が発生した際の経過日数「０」における操業度ＯＰの計算を開始する制御信号を、損傷確率算出部１１、ＰＭＴ算出部１２及び操業度算出部１３に対して出力する。こにれにより、損傷確率算出部１１、ＰＭＴ算出部１２及び操業度算出部１３は、指定されたモデルの情報をモデルデータベース１９のモデルテーブルに書き込む。
損傷確率算出部１１は、製造ライン、電力設備、機械１、機械２、商用電源及び自家発電源の各々の損傷確率（ＰＡ、ＰＢ、Ｐ１、Ｐ２、Ｐ３及びＰ４）を算出し、モデルテーブルに書き込み、記憶させる。
ＰＭＴ算出部１２は、機械１、機械２、商用電源及び自家発電源の各々の最大復旧期待時間（ＰＭＴ１、ＰＭＴ２、ＰＭＴ３、ＰＭＴ４）を、それぞれ損傷確率と復旧日数とを乗算することにより算出し、モデルテーブルに書き込み、記憶させる。
Next, with reference to FIG. 13, an example will be described below in which an operation degree is calculated using a damage probability obtained by linearly processing a damage probability, and a recovery curve generation process is performed. Here, as shown in the model table of FIG. 9, the following description will be made assuming that the damage probability of the
When a simulation start control signal is supplied from the
Next, the recovery
The damage
The
このとき、復旧曲線生成部１４は、モデルテーブルにある最大復旧期待時間の中から最大の日数（本実施形態においては最大復旧期待時間ＰＭＴ４の５日）を抽出し、内部記憶部に最大経過日数として書き込んで記憶させる。
ここで、機械１、機械２、商用電源及び自家発電源の各々の最大復旧期待時間は、すでに説明したように、それぞれ３日、１日、１日、５日と求められる。
At this time, the recovery
Here, the maximum expected recovery time of each of the
また、復旧曲線生成部１４は、モデルテーブルデータベース１９のモデルテーブルにおける各事業構成要素（あるいはサブ事業構成要素）の途中復旧率の欄のデータを読み込み、機械１が損傷確率としてリニアで損傷確率（経過日数毎）を用いることを検出すると、損傷確率算出部１１に対して、機械１の損傷確率と最大復旧期待時間とから、すでに説明したように、経過日数毎に損傷確率（経過日数毎）の算出を指示する。
これにより、損傷確率算出部１１は、経過日数毎の損傷確率（経過日数毎）を算出し、モデルテーブルデータベース１９のモデルテーブルの損傷確率の欄に書き込み、記憶させる。このとき、各経過日数Ｄにおける損傷確率（経過日数毎）は、Ｄ＝０：損傷確率（経過日数毎）＝０．２４０２０、Ｄ＝１：損傷確率（経過日数毎）＝０．１６０１３、Ｄ＝２：損傷確率（経過日数毎）＝１０．８００６、Ｄ＝３：損傷確率（経過日数毎）＝０
また、操業度算出部１３は、損傷確率算出部１１が算出した損傷確率により、それぞれの操業度ＯＰを算出する。
この経過日数「０」日における操業度ＯＰは、すでに算出したように、「０．８７２」である。
このため、復旧曲線生成部１４は、得られた操業度ＯＰ＝０．８７２を経過日数「０」にプロットし、操業度ＯＰが１から０．８７２へ低下してことを表示画面のグラフに表示する。また、復旧曲線生成部１４は、経過日数をインクリメントし、経過日数「１」とする。
In addition, the recovery
Thereby, the damage
Further, the operation
The operation level OP on this elapsed day “0” is “0.872” as already calculated.
For this reason, the recovery
次に、経過日数「１」において、復旧曲線生成部１４は、経過日数が最大経過日数を超えたか否かの判定を行い、超えていないことを検出すると、経過日数を超える最大復旧期待時間の事業構成要素をモデルテーブルから、経過日数と各事業構成要素の最大復旧期待時間とを比較することにより検出する。
検出後、復旧曲線生成部１４は、機械１と自家発電源とが経過日数を超えた事業構成要素とし、機械２と商用電源とが損傷確率を「０」とする事業構成要素とし、それぞれの構成識別情報を含む制御信号を、損傷確率算出部１１へ出力する。このとき、復旧曲線生成部１４は、損傷確率をリニアとするフラグを機械１の構成識別情報に付加して、損傷確率算出部１１へ出力する。
Next, in the elapsed days “1”, the recovery
After the detection, the recovery
そして、損傷確率算出部１１は、上述した制御信号が供給されると、機械２と商用電源との損傷確率を「０」とし、かつモデルテーブルから機械１の経過日数Ｄ＝１に対応する損傷確率（経過日数毎）＝０．１６０１３を読み出し、以下の式により、製造ラインの損傷確率ＰＡ及び電力設備の損傷確率ＰＢの算出を行う。
ＰＡ＝１−（１−Ｐ１）・（１−Ｐ２）
＝１−（１−０．１６０１３）・（１−０）
＝１−０．８３９８７
＝０．１６０１３
ＰＢ＝Ｐ３・Ｐ４
＝０．２２８６０×０
＝０
Then, when the above-described control signal is supplied, the damage
P A = 1- (1-P 1 ) · (1-P 2 )
= 1- (1-0.16013). (1-0)
= 1-0.83987
= 0.16013
P B = P 3・ P 4
= 0.22860 × 0
= 0
そして、損傷確率算出部１１は、算出した経過日数「１」における損傷確率ＰＡ及びＰＢを、操業度算出部１３へ出力する。
操業度算出部１３は、損傷確率算出部１１から供給される損傷確率ＰＡ及びＰＢにより、製造ラインと電力設備とにおける損傷の有無の組み合わせの確率ＰＮＮ（製造ライン及び電力設備の双方に損傷無し）、ＰＮＹ（電力設備のみに損傷有り）、ＰＹＮ（製造ラインのみに損傷有り）、ＰＹＹ（製造ライン及び電力設備の双方に損傷有り）の算出を以下の式により行う。
ＰＮＮ＝（１−ＰＡ）・（１−ＰＢ）
＝（１−１６０１３）×（１−０）
＝０．８３９８７
ＰＮＹ＝（１−ＰＡ）・ＰＢ
＝（１−０．１６０１３）×０
＝０
ＰＹＮ＝ＰＡ・（１−ＰＢ）
＝０．１６０１３×（１−０）
＝０．１６０１３
ＰＹＹ＝ＰＡ・ＰＢ
＝０．１６０１３×０
＝０
Then, the damage
The operation
PNN = (1-P A ) · (1-P B )
= (1-160313) x (1-0)
= 0.83987
P NY = (1-P A ) · P B
= (1-0.16013) x 0
= 0
P YN = P A · (1 -P B)
= 0.16013 x (1-0)
= 0.16013
P YY = P A・ P B
= 0.16013 × 0
= 0
次に、操業度算出部１３は、上述のように求めた確率ＰＮＮ、ＰＮＹ、ＰＹＮ、ＰＹＹの各々を、以下の式に示すように、対応する事業構成要素の損失額に乗算し、乗算結果を加算して、工場の受ける総損失額Ｑを算出する。
Ｑ＝０×ＰＮＮ（損傷がない場合、損失額は０円）
＋２０００×ＰＮＹ（製造ラインにのみ損傷がある場合、損失額は２０００万円）
＋１５００×ＰＹＮ（電力設備にのみ損傷がある場合、損失額は１５００万円）
＋３５００×ＰＹＹ（双方に損傷がある場合、損失額は３５００万円）
＝０×０．８３９８７＋２０００×０＋１５００×０．１６０１３＋３５００×０
＝２４０．１９５
≒２４０
Next, the operation
Q = 0 × P NN (If there is no damage, the loss is 0 yen)
+ 2000 × P NY (If the production line is damaged, the loss is 20 million yen)
+ 1500 × P YN (If there is damage only to the power equipment, the loss is 15 million yen)
+ 3500 × P YY (If both are damaged, the loss is 35 million yen)
= 0 × 0.83987 + 2000 × 0 + 1500 × 0.16013 + 3500 × 0
= 240.195
≒ 240
そして、操業度算出部１３は、以下の式に示すように、算出した経過日数「１」における総損失額Ｑ＝２４０を、予め計算してある総再構成費Ｒ＝３５００により除算し、最大予想損失率ＰＭＬを算出する。
ＰＭＬ＝Ｑ／Ｒ
＝２４０／３５００
＝０．０６８５７
最大予想損失率を算出すると、操業度算出部１３は、以下の式に示すように、１から最大予想損失率ＰＭＬ＝０．０６８５７を減算し、操業度ＯＰを算出する。
ＯＰ＝１−０．０６８５７
＝０．９３１４３
Then, as shown in the following formula, the operation
PML = Q / R
= 240/3500
= 0.06857
When the maximum expected loss rate is calculated, the operation
OP = 1-0.06857
= 0.93143
次に、復旧曲線生成部１４は、得られた操業度ＯＰ＝０．９３１４３を経過日数「１」にプロットし、機械２と商用電源とが損傷から可動状態となり、操業度ＯＰが０．８９７１４３から０．９３１４３へ上昇したことを、表示装置６の表示画面のグラフに表示する。また、復旧曲線生成部１４は、経過日数をインクリメントし、経過日数「２」とする。
Next, the recovery
次に、経過日数「２」において、復旧曲線生成部１４は、経過日数が最大経過日数を超えたか否かの判定を行い、超えていないことを検出すると、経過日数を超える最大復旧期待時間の事業構成要素をモデルテーブルから、経過日数と各事業構成要素の最大復旧期待時間とを比較することにより検出する。
検出後、復旧曲線生成部１４は、機械１と自家発電源とが経過日数を超えた事業構成要素とし、機械２と商用電源とが損傷確率を「０」とする事業構成要素とし、それぞれの構成識別情報を含む制御信号を、損傷確率算出部１１へ出力する。このとき、復旧曲線生成部１４は、損傷確率をリニアとするフラグを機械１の構成識別情報に付加して、損傷確率算出部１１へ出力する。
Next, in the elapsed days “2”, the recovery
After the detection, the recovery
そして、損傷確率算出部１１は、上述した制御信号が供給されると、機械２と商用電源との損傷確率を「０」とし、かつモデルテーブルから機械１の経過日数Ｄ＝２に対応する損傷確率（経過日数毎）＝０．０８００６を読み出し、以下の式により、製造ラインの損傷確率ＰＡ及び電力設備の損傷確率ＰＢの算出を行う。
ＰＡ＝１−（１−Ｐ１）・（１−Ｐ２）
＝１−（１−０．０８００６）・（１−０）
＝１−０．９１９９４
＝０．０８００６
ＰＢ＝Ｐ３・Ｐ４
＝０．２２８６０×０
＝０
Then, when the above-described control signal is supplied, the damage
P A = 1- (1-P 1 ) · (1-P 2 )
= 1- (1-0.08006). (1-0)
= 1-0.91994
= 0.08006
P B = P 3・ P 4
= 0.22860 × 0
= 0
そして、損傷確率算出部１１は、算出した経過日数「２１」における損傷確率ＰＡ及びＰＢを、操業度算出部１３へ出力する。
操業度算出部１３は、損傷確率算出部１１から供給される損傷確率ＰＡ及びＰＢにより、製造ラインと電力設備とにおける損傷の有無の組み合わせの確率ＰＮＮ（製造ライン及び電力設備の双方に損傷無し）、ＰＮＹ（電力設備のみに損傷有り）、ＰＹＮ（製造ラインのみに損傷有り）、ＰＹＹ（製造ライン及び電力設備の双方に損傷有り）の算出を以下の式により行う。
ＰＮＮ＝（１−ＰＡ）・（１−ＰＢ）
＝（１−０．０８００６）×（１−０）
＝０．９１９９４
ＰＮＹ＝（１−ＰＡ）・ＰＢ
＝（１−０．０８００６）×０
＝０
ＰＹＮ＝ＰＡ・（１−ＰＢ）
＝０．０８００６×（１−０）
＝０．０８００６
ＰＹＹ＝ＰＡ・ＰＢ
＝０．０８００６×０
＝０
Then, the damage
The operation
PNN = (1-P A ) · (1-P B )
= (1-0.08006) × (1-0)
= 0.91994
P NY = (1-P A ) · P B
= (1-0.08006) x 0
= 0
P YN = P A · (1 -P B)
= 0.08006 × (1-0)
= 0.08006
P YY = P A・ P B
= 0.08006 × 0
= 0
次に、操業度算出部１３は、上述のように求めた確率ＰＮＮ、ＰＮＹ、ＰＹＮ、ＰＹＹの各々を、以下の式に示すように、対応する事業構成要素の損失額に乗算し、乗算結果を加算して、工場の受ける総損失額Ｑを算出する。
Ｑ＝０×ＰＮＮ（損傷がない場合、損失額は０円）
＋２０００×ＰＮＹ（製造ラインにのみ損傷がある場合、損失額は２０００万円）
＋１５００×ＰＹＮ（電力設備にのみ損傷がある場合、損失額は１５００万円）
＋３５００×ＰＹＹ（双方に損傷がある場合、損失額は３５００万円）
＝０×０．９１９９４＋２０００×０＋１５００×０．０８００６＋３５００×０
＝１２０．０９
≒１２０
Next, the operation
Q = 0 × P NN (If there is no damage, the loss is 0 yen)
+ 2000 × P NY (If the production line is damaged, the loss is 20 million yen)
+ 1500 × P YN (If there is damage only to the power equipment, the loss is 15 million yen)
+ 3500 × P YY (If both are damaged, the loss is 35 million yen)
= 0 × 0.91994 + 2000 × 0 + 1500 × 0.08006 + 3500 × 0
= 120.09
≒ 120
そして、操業度算出部１３は、以下の式に示すように、算出した経過日数「１」における総損失額Ｑ＝２４０を、予め計算してある総再構成費Ｒ＝３５００により除算し、最大予想損失率ＰＭＬを算出する。
ＰＭＬ＝Ｑ／Ｒ
＝１２０／３５００
＝０．０３４２８
最大予想損失率を算出すると、操業度算出部１３は、以下の式に示すように、１から最大予想損失率ＰＭＬ＝０．０３４２８を減算し、操業度ＯＰを算出する。
ＯＰ＝１−０．０３４２８
＝０．９６５７１
Then, as shown in the following formula, the operation
PML = Q / R
= 120/3500
= 0.03428
When the maximum expected loss rate is calculated, the operation
OP = 1-0.03428
= 0.96571
次に、復旧曲線生成部１４は、得られた操業度ＯＰ＝０．９６５７１を経過日数「２」にプロットし、機械２と商用電源とが損傷から可動状態となり、操業度ＯＰが０．９３１４３から０．９６５７１へ上昇したことを、表示装置１６の表示画面のグラフに表示する。また、復旧曲線生成部１４は、経過日数をインクリメントし、経過日数「３」とする。
Next, the recovery
次に、復旧曲線生成部１４は、経過日数を超える最大復旧期待時間の事業構成要素をモデルテーブルから、経過日数と各事業構成要素の最大復旧期待時間とを比較することにより検出する。
検出後、復旧曲線生成部１４は、自家発電源とが経過日数を超えた事業構成要素とし、機械１と機械２と商用電源とが損傷確率を「０」とする事業構成要素とし、それぞれの構成識別情報を含む制御信号を、損傷確率算出部１１へ出力する。
そして、損傷確率算出部１１は、上述した制御信号が供給されると、機械１、機械２及び商用電源との損傷確率を「０」とし、以下の式により、製造ラインの損傷確率ＰＡ及び電力設備の損傷確率ＰＢの算出を行う。
ＰＡ＝１−（１−Ｐ１）・（１−Ｐ２）
＝１−（１−０）・（１−０）
＝１−０
＝０
ＰＢ＝Ｐ３・Ｐ４
＝０．２２８６０×０
＝０
Next, the recovery
After the detection, the recovery
Then, damage
P A = 1- (1-P 1 ) · (1-P 2 )
= 1- (1-0). (1-0)
= 1-0
= 0
P B = P 3・ P 4
= 0.22860 × 0
= 0
そして、損傷確率算出部１１は、算出した経過日数「３」における損傷確率ＰＡ及びＰＢを、操業度算出部１３へ出力する。
操業度算出部１３は、損傷確率算出部１１から供給される損傷確率ＰＡ及びＰＢにより、製造ラインと電力設備とにおける損傷の有無の組み合わせの確率ＰＮＮ（製造ライン及び電力設備の双方に損傷無し）、ＰＮＹ（電力設備のみに損傷有り）、ＰＹＮ（製造ラインのみに損傷有り）、ＰＹＹ（製造ライン及び電力設備の双方に損傷有り）の算出を以下の式により行う。
ＰＮＮ＝（１−ＰＡ）・（１−ＰＢ）
＝（１−０）×（１−０）
＝１
ＰＮＹ＝（１−ＰＡ）・ＰＢ
＝（１−０）×０
＝０
ＰＹＮ＝ＰＡ・（１−ＰＢ）
＝０×（１−０）
＝０
ＰＹＹ＝ＰＡ・ＰＢ
＝０×０
＝０
Then, the damage
The operation
PNN = (1-P A ) · (1-P B )
= (1-0) × (1-0)
= 1
P NY = (1-P A ) · P B
= (1-0) × 0
= 0
P YN = P A · (1 -P B)
= 0 × (1-0)
= 0
P YY = P A・ P B
= 0x0
= 0
次に、操業度算出部１３は、上述のように求めた確率ＰＮＮ、ＰＮＹ、ＰＹＮ、ＰＹＹの各々を、以下の式に示すように、対応する事業構成要素の損失額に乗算し、乗算結果を加算して、工場の受ける総損失額Ｑを算出する。
Ｑ＝０×ＰＮＮ（損傷がない場合、損失額は０円）
＋２０００×ＰＮＹ（製造ラインにのみ損傷がある場合、損失額は２０００万円）
＋１５００×ＰＹＮ（電力設備にのみ損傷がある場合、損失額は１５００万円）
＋３５００×ＰＹＹ（双方に損傷がある場合、損失額は３５００万円）
＝０×１＋２０００×０＋１５００×０＋３５００×０
＝０
Next, the operation
Q = 0 × P NN (If there is no damage, the loss is 0 yen)
+ 2000 × P NY (If the production line is damaged, the loss is 20 million yen)
+ 1500 × P YN (If there is damage only to the power equipment, the loss is 15 million yen)
+ 3500 × P YY (If both are damaged, the loss is 35 million yen)
= 0x1 + 2000x0 + 1500x0 + 3500x0
= 0
そして、操業度算出部１３は、以下の式に示すように、算出した経過日数「３」における総損失額Ｑ＝０を、予め計算してある総再構成費Ｒ＝３５００により除算し、最大予想損失率ＰＭＬを算出する。
ＰＭＬ＝Ｑ／Ｒ
＝０／３５００
＝０
最大予想損失率を算出すると、操業度算出部１３は、以下の式に示すように、１から最大予想損失率ＰＭＬ＝０を減算し、操業度ＯＰを算出する。
ＯＰ＝１−０
＝１
Then, as shown in the following equation, the operation
PML = Q / R
= 0/3500
= 0
When the maximum expected loss rate is calculated, the operation
OP = 1-0
= 1
次に、復旧曲線生成部１４は、得られた操業度ＯＰ＝１を経過日数「３」にプロットし、機械１、機械２及び商用電源が損傷状態から可動状態となり、操業度ＯＰが０．８９７１４３から１のハザード前の操業度に戻ったことを、表示装置１６の表示画面のグラフに表示する。また、復旧曲線生成部１４は、経過日数をインクリメントし、経過日数「４」とする。
Next, the recovery
次に、経過日数「４」において、復旧曲線生成部１４は、経過日数が最大経過日数を超えたか否かの判定を行い、超えていないことを検出すると、経過日数を超える最大復旧期待時間の事業構成要素をモデルテーブルから、経過日数と各事業構成要素の最大復旧期待時間とを比較することにより検出する。
検出後、復旧曲線生成部１４は、自家発電源とが経過日数を超えた事業構成要素とし、機械１と機械２と商用電源とが損傷確率を「０」とする事業構成要素として検出するが、自家発電源がリニアの損傷確率を使用するとして設定されておらず、自家発電源の損傷確率が前回と同一であることを検出する。
このため、復旧曲線生成部１４は、新たに損傷状態から稼動状態に遷移した事業構成要素がないので、前回算出した操業度ＯＰから変化がないものとし、前回の操業度ＯＰを表示装置１６の表示画面のグラフに表示する。また、復旧曲線生成部１４は、経過日数をインクリメントし、経過日数「５」とする。
Next, in the elapsed days “4”, the recovery
After the detection, the recovery
For this reason, since there is no business component which newly changed from the damage state to the operation state, the recovery
次に、経過日数「５」において、復旧曲線生成部１４は、経過日数が最大経過日数を超えたか否かの判定を行い、超えていないことを検出すると、経過日数を超える最大復旧期待時間の事業構成要素をモデルテーブルから、経過日数と各事業構成要素の最大復旧期待時間とを比較することにより検出する。
しかしながら、復旧曲線生成部１４は、経過日数を超える最大復旧期待時間を超える事業構成要素が存在せず、すでに操業度ＯＰが「１」となっていることから、通常に工場が稼動している状態にあるため、前回算出した操業度ＯＰから変化がないものとし、前回の操業度ＯＰを表示装置１６の表示画面のグラフに表示する。また、復旧曲線生成部１４は、経過日数をインクリメントし、経過日数「６」とする。
Next, in the elapsed days “5”, the recovery
However, the restoration
次に、経過日数「６」において、復旧曲線生成部１４は、経過日数が最大経過日数を超えたか否かの判定を行い、超えてることを検出すると、復旧曲線の生成処理を終了する。
上述したように、復旧曲線生成部１４は、図１３における工場のリソースの復旧曲線を生成する。
また、リニアな損傷確率を使用する復旧曲線の生成−サービス（商品の生産）については、リソースと同様の演算を行うため、説明を省略する。
Next, in the elapsed days “6”, the recovery
As described above, the recovery
Moreover, since the calculation similar to a resource is performed about the generation | occurrence | production-service (production of goods) of the recovery curve which uses a linear damage probability, description is abbreviate | omitted.
上述した本実施形態によれば、損傷確率の変化にステップではなくリニアを用いることにより、複数の機械などが損傷し、数台ずつ復旧して生産が可能になる場合などにおいて、より詳細な復旧曲線を生成し、ハザードに対して耐性を持たせたり、あるいは損傷を受けた後の修理の優先順位などの対処方法の抽出の精度を向上させることができる。 According to the above-described embodiment, by using a linear instead of a step for changing the damage probability, more detailed restoration is possible when a plurality of machines are damaged and production is possible by recovering several units at a time. A curve can be generated to make it resistant to hazards, or to improve the accuracy of extraction of countermeasures such as the priority of repair after being damaged.
また、図１における損傷確率算出部１１、ＰＭＴ算出部１２、操業度算出部１３、復旧曲線生成部１４の機能を実現するためのプログラムをコンピュータ読み取り可能な記録媒体に記録して、この記録媒体に記録されたプログラムをコンピュータシステムに読み込ませ、実行することにより復旧曲線の生成処理を行ってもよい。なお、ここでいう「コンピュータシステム」とは、ＯＳや周辺機器等のハードウェアを含むものとする。
Further, a program for realizing the functions of the damage
また、「コンピュータシステム」は、ＷＷＷシステムを利用している場合であれば、ホームページ提供環境（あるいは表示環境）も含むものとする。
また、「コンピュータ読み取り可能な記録媒体」とは、フレキシブルディスク、光磁気ディスク、ＲＯＭ、ＣＤ−ＲＯＭ等の可搬媒体、コンピュータシステムに内蔵されるハードディスク等の記憶装置のことをいう。さらに「コンピュータ読み取り可能な記録媒体」とは、インターネット等のネットワークや電話回線等の通信回線を介してプログラムを送信する場合の通信線のように、短時間の間、動的にプログラムを保持するもの、その場合のサーバやクライアントとなるコンピュータシステム内部の揮発性メモリのように、一定時間プログラムを保持しているものも含むものとする。また上記プログラムは、前述した機能の一部を実現するためのものであっても良く、さらに前述した機能をコンピュータシステムにすでに記録されているプログラムとの組み合わせで実現できるものであっても良い。
Further, the “computer system” includes a homepage providing environment (or display environment) if a WWW system is used.
The “computer-readable recording medium” refers to a storage device such as a flexible medium, a magneto-optical disk, a portable medium such as a ROM and a CD-ROM, and a hard disk incorporated in a computer system. Furthermore, the “computer-readable recording medium” dynamically holds a program for a short time like a communication line when transmitting a program via a network such as the Internet or a communication line such as a telephone line. In this case, a volatile memory in a computer system serving as a server or a client in that case, and a program that holds a program for a certain period of time are also included. The program may be a program for realizing a part of the functions described above, and may be a program capable of realizing the functions described above in combination with a program already recorded in a computer system.
以上、この発明の実施形態を図面を参照して詳述してきたが、具体的な構成はこの実施形態に限られるものではなく、この発明の要旨を逸脱しない範囲の設計等も含まれる。 The embodiment of the present invention has been described in detail with reference to the drawings. However, the specific configuration is not limited to this embodiment, and includes design and the like within a scope not departing from the gist of the present invention.
１１…損傷確率算出部
１２…ＰＭＴ算出部
１３…操業度算出部
１４…復旧曲線生成部
１５…入力装置
１６…表示装置
１７…損傷確率データベース
１８…基本情報データベース
１９…モデルデータベース
１００…復旧曲線作成システム
DESCRIPTION OF
Claims (6)
診断対象の施設における少なくとも建物および各収容設備を構成要素とし、当該構成要素がハザードの大きさにより損傷を受ける損傷確率を、当該構成要素毎に記憶する損傷確率記憶部と、
前記構成要素の全壊における最大復旧時間に対し、前記損傷確率を乗じ、前記構成要素毎の最大復旧期待時間を求めるＰＭＴ算出部と、
所定の周期毎に、前記構成要素毎の全壊における最大損失額に対し、前記損傷確率を乗じて最大予想損失額を求め、前記最大予想損失額を前記構成要素の再構成費で除算して予想損失率を求め、１から前記予想損失率を減算し、減算結果を操業度として出力する操業度算出部と、
前記操業度を前記周期毎にプロットし、復旧曲線を生成する復旧曲線生成部と
を有し、
前記操業度算出部が、前記周期の経過毎に経過時間をカウントし、当該経過時間が前記最大復旧期待時間を超えた前記構成要素の損傷確率を０とし、前記予想損失率を再度計算することを特徴とする復旧曲線作成システム。 A recovery curve generation system that quantitatively evaluates damage to facilities caused by hazards and generates a recovery curve that predicts business continuity in the facility.
A damage probability storage unit that stores at least the building and each accommodation facility in the facility to be diagnosed as a component, and stores the damage probability that the component is damaged due to the size of the hazard for each component,
A PMT calculation unit that obtains the maximum expected recovery time for each component by multiplying the maximum recovery time in the total destruction of the component by the damage probability;
For each predetermined period, a maximum expected loss amount is obtained by multiplying the maximum loss amount in total destruction for each component by the probability of damage, and the maximum expected loss amount is divided by the reconfiguration cost of the component. An operation rate calculation unit that obtains a loss rate, subtracts the expected loss rate from 1, and outputs the subtraction result as an operation rate;
A plot of the degree of operation for each cycle, and a recovery curve generator for generating a recovery curve,
The operation degree calculation unit counts the elapsed time every time the period elapses, sets the damage probability of the component whose elapsed time exceeds the maximum expected recovery time to 0, and calculates the expected loss rate again. A recovery curve creation system characterized by
さらに有することを特徴とする請求項１に記載の復旧曲線作成システム。 The elapsed time is counted for each elapse of the cycle, and the damage probability is multiplied by the result obtained by dividing the elapsed time by the maximum expected recovery time for each component, and the result of the multiplication is calculated from the damage probability. 2. The apparatus according to claim 1, further comprising a damage probability calculation unit that calculates a new damage probability corresponding to the elapsed time by subtraction and outputs the calculated new damage probability to the operation degree calculation unit. The recovery curve creation system described.
さらに有することを特徴とする請求項１から請求項３のいずれか一項に記載の復旧曲線作成システム。 The result of multiplying the production value of the product produced using the component for each cycle by the damage probability of each component until the maximum recovery time of the component is exceeded. The recovery curve creation system according to any one of claims 1 to 3, further comprising a sales loss amount calculation unit that accumulates in a period and obtains a sales loss amount.
ＰＭＴ算出部が、診断対象の施設における少なくとも建物および各収容設備を構成要素とし、当該構成要素がハザードの大きさにより損傷を受ける損傷確率を、当該構成要素毎に記憶する損傷確率記憶部から、前記構成要素の前記損傷確率を読み出し、前記構成要素の全壊における最大復旧時間に対し、読み出した前記損傷確率を乗じ、前記構成要素毎の最大復旧期待時間を求めるＰＭＴ算出過程と、
操業度算出部が、所定の周期毎に、前記損傷確率記憶部から構成要素の前記損傷確率を読み出し、前記構成要素毎の全壊における最大損失額に対し、読み出した前記損傷確率を乗じて最大予想損失額を求め、前記最大予想損失額を前記構成要素の再構成費で除算して予想損失率を求め、１から前記予想損失率を減算し、減算結果を操業度として出力する操業度算出過程と、
復旧曲線生成部が前記操業度を前記周期毎にプロットし、復旧曲線を生成する復旧曲線生成過程と
を有し、
前記操業度算出部が、前記周期の経過毎に経過時間をカウントし、当該経過時間が前記最大復旧期待時間を超えた前記構成要素の損傷確率を０とし、前記予想損失率を再度計算することを特徴とする復旧曲線作成方法。 It is a control method of a recovery curve generation system that quantitatively evaluates damage to facilities due to hazards and generates a recovery curve that predicts business continuity in the facility,
From the damage probability storage unit that stores, for each component, the damage probability that the PMT calculation unit has at least the building and each accommodation facility in the facility to be diagnosed as components, and the component is damaged due to the size of the hazard. A PMT calculation step of reading the damage probability of the component, multiplying the read recovery probability by the maximum recovery time in the total destruction of the component, and obtaining a maximum expected recovery time for each component;
The operation rate calculation unit reads out the damage probability of the component from the damage probability storage unit every predetermined period, and multiplies the maximum loss amount in the total destruction for each component by the read out damage probability. An operation degree calculation process for obtaining a loss amount, dividing the maximum expected loss amount by a reconfiguration cost of the component to obtain an expected loss rate, subtracting the expected loss rate from 1, and outputting the subtraction result as an operation degree When,
A recovery curve generation unit plots the operation rate for each cycle and generates a recovery curve;
The operation degree calculation unit counts the elapsed time every time the period elapses, sets the damage probability of the component whose elapsed time exceeds the maximum expected recovery time to 0, and calculates the expected loss rate again. A recovery curve creation method characterized by
ＰＭＴ算出部が、診断対象の施設における少なくとも建物および各収容設備を構成要素とし、当該構成要素がハザードの大きさにより損傷を受ける損傷確率を、当該構成要素毎に記憶する損傷確率記憶部から、前記構成要素の前記損傷確率を読み出し、前記構成要素の全壊における最大復旧時間に対し、読み出した前記損傷確率を乗じ、前記構成要素毎の最大復旧期待時間を求めるＰＭＴ算出処理と、
操業度算出部が、所定の周期毎に、前記損傷確率記憶部から構成要素の前記損傷確率を読み出し、前記構成要素毎の全壊における最大損失額に対し、読み出した前記損傷確率を乗じて最大予想損失額を求め、前記最大予想損失額を前記構成要素の再構成費で除算して予想損失率を求め、１から前記予想損失率を減算し、減算結果を操業度として出力する操業度算出処理と、
復旧曲線生成部が前記操業度を前記周期毎にプロットし、復旧曲線を生成する復旧曲線生成処理と、
前記操業度算出部が、前記周期の経過毎に経過時間をカウントし、当該経過時間が前記最大復旧期待時間を超えた前記構成要素の損傷確率を０とし、前記予想損失率を再度計算する処理と
をコンピュータに実行させるプログラム。 A program that causes a computer to execute the operation of a recovery curve generation system that quantitatively evaluates damage to facilities due to hazards and generates a recovery curve that predicts business continuity in the facility,
From the damage probability storage unit that stores, for each component, the damage probability that the PMT calculation unit has at least the building and each accommodation facility in the facility to be diagnosed as components, and the component is damaged due to the size of the hazard. PMT calculation processing for reading the damage probability of the component, multiplying the read recovery probability by the maximum recovery time in the total destruction of the component, and obtaining a maximum expected recovery time for each component;
The operation rate calculation unit reads out the damage probability of the component from the damage probability storage unit every predetermined period, and multiplies the maximum loss amount in the total destruction for each component by the read out damage probability. An operation degree calculation process for obtaining a loss amount, dividing the maximum expected loss amount by a reconfiguration cost of the component to obtain an expected loss rate, subtracting the expected loss rate from 1, and outputting the subtraction result as an operation rate When,
A recovery curve generation unit plots the operation rate for each cycle, and generates a recovery curve, a recovery curve generation process;
A process in which the operation degree calculation unit counts the elapsed time every time the cycle elapses, sets the damage probability of the component whose elapsed time exceeds the maximum expected recovery time to 0, and recalculates the expected loss rate A program that causes a computer to execute and.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2010288012A JP5592247B2 (en) | 2010-12-24 | 2010-12-24 | Recovery curve creation system, recovery curve creation method and program |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2010288012A JP5592247B2 (en) | 2010-12-24 | 2010-12-24 | Recovery curve creation system, recovery curve creation method and program |
Related Child Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2014155397A Division JP5970505B2 (en) | 2014-07-30 | 2014-07-30 | Recovery curve creation system, recovery curve creation method, and program |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2012137827A JP2012137827A (en) | 2012-07-19 |
| JP5592247B2 true JP5592247B2 (en) | 2014-09-17 |
Family
ID=46675218
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2010288012A Expired - Fee Related JP5592247B2 (en) | 2010-12-24 | 2010-12-24 | Recovery curve creation system, recovery curve creation method and program |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP5592247B2 (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP6074264B2 (en) * | 2012-12-28 | 2017-02-01 | 大和ハウス工業株式会社 | Method, apparatus and program for creating fragility curves for existing wooden houses |
| CN108241912B (en) * | 2016-12-23 | 2021-01-29 | 中国石油天然气股份有限公司 | Method and device for determining profit and loss balance period of oil pumping well |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2008186285A (en) * | 2007-01-30 | 2008-08-14 | Takenaka Komuten Co Ltd | Earthquake damage estimation device, method, and program |
| JP5186990B2 (en) * | 2008-04-22 | 2013-04-24 | 株式会社大林組 | Computer for recovery process evaluation and recovery process evaluation program |
-
2010
- 2010-12-24 JP JP2010288012A patent/JP5592247B2/en not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| JP2012137827A (en) | 2012-07-19 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Fan et al. | Reliability of the Box–Jenkins model for forecasting construction demand covering times of economic austerity | |
| Stojadinović et al. | Rapid earthquake loss assessment based on machine learning and representative sampling | |
| US20090083089A1 (en) | Systems and methods for analyzing failure modes according to cost | |
| Cauzzi et al. | Earthquake early warning and operational earthquake forecasting as real-time hazard information to mitigate seismic risk at nuclear facilities | |
| US8768742B2 (en) | Assessment and rationalization of resiliency of data center strategies | |
| Hariri-Ardebili et al. | Seismic risk prioritization of a large portfolio of dams: Revisited | |
| US20170178039A1 (en) | Computer system for catastrophic event management | |
| Zorn et al. | Quantifying directional dependencies from infrastructure restoration data | |
| DeBock et al. | A comparative evaluation of probabilistic regional seismic loss assessment methods using scenario case studies | |
| Meng et al. | The asymmetric impact of natural disasters on China's bilateral trade | |
| Sweya et al. | Tool development to measure the resilience of water supply systems in Tanzania: Economic dimension | |
| JP2013105221A (en) | Maintenance parts delivery support system, maintenance parts delivery support device and maintenance parts delivery support program | |
| JP5970505B2 (en) | Recovery curve creation system, recovery curve creation method, and program | |
| JP5592247B2 (en) | Recovery curve creation system, recovery curve creation method and program | |
| Azarbakht et al. | A decision-making approach for operational earthquake forecasting | |
| JP2004054954A (en) | Risk diagnostic system, method of generating risk map data, and program | |
| US20130282410A1 (en) | Hazard risk assessment | |
| Markhvida et al. | Modeling future economic costs and interdependent industry recovery after earthquakes | |
| CN111666191A (en) | Data quality monitoring method and device, electronic equipment and storage medium | |
| WO2013061324A2 (en) | A method for estimating the total cost of ownership (tcp) for a requirement | |
| US20160026951A1 (en) | Mature practice assessment framework | |
| JP6550254B2 (en) | Recovery curve creation system, recovery curve creation method, and program | |
| Atef et al. | Risk assessment of infrastructure facilities considering spatial and operational interdependencies: temporal simulation model | |
| Sheikhalishahi et al. | Combinatorial optimization of petrochemical plants by asset integrity management indicators | |
| Hu et al. | Framework for prioritizing geospatial data processing tasks during extreme weather events |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| RD02 | Notification of acceptance of power of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7422 Effective date: 20130612 |
|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20130919 |
|
| RD04 | Notification of resignation of power of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7424 Effective date: 20130912 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20140522 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20140701 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20140731 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 5592247 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| LAPS | Cancellation because of no payment of annual fees |