Multi-sites fisheries modeling with Itô stochastic differential equations and stochastic optimization problem

技術要約：伊藤型確率微分方程式と確率的最適化を用いたマルチサイト漁業モデリング

問題提起
本論文は、集魚装置（FADs）を利用するシステムに焦点を当てた、マルチサイト漁業のモデリングと最適化について取り組んでいる。決定論的なモデルが歴史的に個体群動態を支配してきたが、それらは気候変動や生態学的攪乱といった自然の変動性や不確実性を考慮できていない。著者らは、複数のサイト（$L$）および沖合のストック（$s$）における魚類バイオマスと漁獲努力量の進化をモデル化するために、確率過程を用いたフレームワークを提案している。核心となる課題は、標準的な確率微分方程式（SDE）の定式化において計算コストが膨大になりがちな行列の平方根の計算を回避しつつ、出生、死亡、移動、および漁獲活動のランダムな変動を捉えるシステムを構築することである。さらに、本研究は、これらの動的で不確実な制約条件下で漁獲利益を最大化するための確率的最適化問題を解くことを目的としている。

手法
本研究は、確率的モデリングと確率的最適化という二部構成の手法を採用している。

1. 確率的モデリング（SDEの定式化）:

   • 離散から連続へ: 著者らはまず、$L$ 個のサイトを持つシステムの離散時間確率過程を定義する。微小時間 $\Delta t$ におけるすべての可能な状態変化（出生、死亡、サイト間および沖合への移動、および船舶の移動）を列挙する。
   • ドリフトと拡散: これらの離散的な遷移から、期待値ベクトル（ドリフト項、$f$）と共分散行列（$V$）を導出する。
   • 等価なシステム: 標準的な SDE 定式化 $dX_t = f(t, X)dt + B(t, X)dW(t)$ は、$B = V^{1/2}$ を必要とする。著者らは、大きな $(2L+1) \times (2L+1)$ 行列の平方根を計算することは計算負荷が高く、大きな $L$ に対して明示的な公式が存在しないことを指摘している。
   • 行列分解: 直接的な行列平方根の計算を避けるため、著者らは等価な定式化 $dX_t = f(t, X)dt + H(t, X)dW^{**}(t)$ （ここで $V = HH^T$）を利用する。行列 $H$ は、出生、死亡、移動、および船舶の移動を表すブロック行列を用いて明示的に構築される。このアプローチは、基礎となる遷移確率から導出される $H$ の明示的な構造を利用することで、直接 $V^{1/2}$ を計算する必要をなくし、計算コストを回避するものである。$L \ge 3$ の場合、行列 $H$ の次元は $(2L+1) \times (L+1)^2$ となり、状態の次元よりも列が多くなるが、この明示的な構築によって全共分散行列 $V$ の分解に伴う数値的コストを回避できる。
   • 一般化: モデルはまず $L=2$ のサイトに対して詳細に記述され、その後、ロジスティック成長および特定の捕獲率を組み込んだ $L \ge 3$ のサイトへと一般化される。

2. 確率的最適化:

   • 目的関数: 時間ホライゾン $[0, T]$ における期待総利益の最大化を目的とする。これは、収益（漁獲価格 $\times$ 捕獲係数 $\times$ バイオマス $\times$ 努力量）とコスト（単位努力量あたりのコスト）の差として定義される。
   • 制約条件: 最適化は、導出された SDE の動力学に従う。停止時刻 $\tau$ は、ストックがゼロになるか、または時刻 $T$ に達した時のいずれか早い方の時刻として定義される。
   • 動的計画法: この問題は動的計画法の原理を用いて解かれ、ハミルトン・ヤコビ・ベルマン（HJB）方程式へと導かれる。価値関数 $U(x, t)$ は、将来の最大期待利益を表す。
   • 制御のための簡略化: 最適化フェーズでは、制御変数（漁獲努力量）に関連する生物学的および捕獲のダイナミクスに焦点を当てるため、船舶の移動コストと売上に関連する項を除外して、拡散行列 $H$ を $H_1$ に簡略化する。

主な貢献

• 等価な SDE システム: 本論文は、標準的な SDE 定式化（$V^{1/2}$ を必要とするもの）と、 $V=HH^T$ となる代替的な定式化（行列 $H$ を用いるもの）の間の等価性を提案し、確立している。この代替案は、ブロック行列から $H$ を明示的に構築することにより、高次元システムにおける行列平方根の計算負担を回避できるため、理論的に健全であり、かつ計算効率が高い。
• マルチサイト確率フレームワーク: 著者らは、サイト間および沖合のストック間の魚類の移動、ならびに漁船のサイト間の移動を明示的に考慮した、マルチサイト漁業のための包括的な確率モデルを開発した。
• 確率的最適化の解: 本研究は、マルチサイト漁業のための確率的最適制御問題を定式化し、解決している。SDE の動力学を HJB 方程式と結合させることで、不確実性下で利益を最大化する最適な漁獲努力量戦略を導き出している。

結果

• 数値的検証: 著者らは、$L=2$ および $L=10$ のサイトの場合について、オイラー・丸山法を用いて 500 個のサンプルパスによる数値シミュレーションを行った。
• 等価性の確認: 標準的な定式化（行列 $B=V^{1/2}$ を使用）と提案された定式化（行列 $H$ を使用）の比較により、バイオマスおよび漁獲努力量の期待値がほぼ同一であることが示された。両手法間の絶対誤差は低く、理論的な等価性と、計算の複雑さを軽減するための $H$ 定式化の有用性が検証された。
• 最適化の結果: HJB 方程式のシミュレーションにより、価値関数はストックおよびバイオマスのレベルとともに増加し、漁獲努力量とともに減少することが明らかになった。モデルから導出された最適な漁獲努力量の方策は、確率的なショックが存在する場合でも、経済的リターンと生態学的持続可能性のバランスを取ることが可能であることを示している。

意義と主張
本論文は、理論的モデリングと数値的効率性の間の溝を埋める厳密な確率的フレームワークを提供することで、漁業管理の分野に貢献すると主張している。

• 理論的貢献: 高次元のマルチサイト・システムにおける標準的な SDE 定式化に対し、数学的に等価でありながら計算負荷の低い代替案を提示しており、複雑で高次元な漁業のモデリングをより実行可能なものにしている。
• 実用的応用: 結果は、固有の不確実性に直面している現代の漁業において、堅牢な制御のために確率的最適化手法が必要であることを示唆している。提案されたフレームワークは、適応的な収穫戦略および持続可能なマルチサイト政策の開発を支援するものである。
• 限界と今後の課題: 著者らは、モデルが一定のパラメータ（例：移動率）を仮定しており、年齢構造や複雑な種間相互作用（捕食者-被食者関係）を含んでいないことを謙虚に認めている。彼らは、将来の研究として、エコシステムベースの管理への拡張、栄養段階の相互作用の組み込み、および漁業モデルと気候モデルの結合を挙げている。また、本研究で使用されたオイラー・丸山法を超えた、より高度な数値スキーム（ミルシュタイン法、ルンゲ・クッタ法など）の探索も提案している。