[仏教の最適化プログラムの続き。] 水害(洪水)のモデル化。…

jrf> [仏教の最適化プログラムの続き。]

水害(洪水)のモデル化。

水害の頻度は、月によって異なり、月に何度か襲うこともある。ただ、一日には一度にしても問題ないだろうから、一日の確率に対し 30 日分サイコロを振って、水害を発生させる感じにしよう。一日の確率は、6月・9月・10月は 4/30、7月・8月は 3/30、それ以外の月は、0.5/30 ぐらいにとりあえずしておくか。

水害の強さも、対防御レベルで、7 から 3 ぐらいまでで uniform なランダムでいいのではないか。とりあえずは。7 のほうが少なくなるよう傾斜を付けたほうがリアルかもしれないが、「災害」レベルになるかどうかでまず判断があり、「災害」レベルに決まれば、あとは運で uniform のほうが適当かな…と今のところは考えている。

で水害でどれぐらい被害が出るかだが、10000人の人口で、レベル 7 で何もしなければ 1000 人死ぬ。財産の被害は 4000人ぐらいに出る。農作物の被害は半数の土地で被害が出るが、1月-4月・11月・12月は、作物が農地の上にないので被害が出ないとしよう。倉庫の穀物の被害とかありうるが、それはとりあえず考えないことにする。そして、堤防にレベル 3 のダメージがある。…と。

レベル 5 だと、100人死に、レベル 1 だと一ケタしか死なない。…というのを scipy.optimze.curve_fit で求めようとして、まず、二次関数だとどうも妙な形になる…というので、a * np.exp(X - b) + c を fit させると、curve_fit の仕様のせいなのかどうもうまくいかない…が、結局、a = 1 b = 0 c = 0 すなわち 1 * np.exp(X) で、だいたいイイ感じになる。1, 5, 7 のときそれぞれ 2.7, 148, 1097 ぐらいになる。

レベル 3 未満は通常の死などに含まれるとして災害としては考えない感じでいいと思う。

で、水害に対して対策をするとどうなるか。防御レベルは最大 1.5 上がり、訓練レベルも最大 1.0 上がる感じか。防御が成功…堤防が防ぎきったときは、被害は堤防のみで、レベル 7 の水害で、堤防が レベル 1.5 下がるぐらい。堤防が防ぎ切れなかったときは、元々あ死者が 1000 人なのが、訓練レベルが最大の 3 で、500人ぐらいまで減る感じ…1 レベル下の災害になる感じか。