石油・ガス生産システムのリアルタイムモニタリングおよび最適化に向けた制御認識型デジタルツインフレームワークの開発

サマリー

本研究では、機械学習と物理モデルに基づいた、石油およびガスシステムのモニタリングと最適化のための制御対応デジタルツインフレームワークを提示し、シミュレーションおよびベンチマークデータセットを用いて評価しています。

要約

本論文では、石油およびガス生産システムの追跡と改善を目的とした、制御を考慮したデジタルツインフレームワークを提案します。本手法は、時空間グラフニューラルネットワーク、ニューラル4D-Varデータ同化、強化学習ベースの制御、および物理情報ニューラルオペレーター（FNOおよびニューラルオペレーターモデル）を組み合わせています。本フレームワークの評価には、ベンチマークとなるサイバーフィジカルシステムデータセット（SWaTおよびWADI）と、シミュレーションによる石油・ガスデータセットを使用しました。平均値を用いた複数の独立した実行によって確認された結果は、状態推定、異常検知、および制御最適化における性能の向上を示しています。ただし、評価はベンチマークおよびシミュレーションデータセットに限定されており、実用的な適用性を検証するためには、実際の産業データを用いたさらなる検証が必要です。

概要

世界的なエネルギー需要の増大に伴い、安全で効率的な石油・ガス輸送ネットワークの必要性が高まっています。しかし、インフラの老朽化とシステムの複雑化により、パイプラインの信頼性を維持することが極めて困難になっています。この点において、デジタルツイン（DT）技術は、産業現場におけるリスク評価、予知保全、およびシステム最適化を可能にするデータ駆動型のライフサイクル管理戦略として注目を集めています^1,2,3,4,5。産業システムにおけるデジタルトランスフォーメーションは、事後保全からインテリジェントなデータ駆動型意思決定へと移行することで、異常検知、故障予測、およびシステムの信頼性を向上させてきました^6,7,8。プロセス最適化、予知保全、およびスマート製造におけるDTの応用は、多くの研究の対象となってきました。例えば、故障の予知分析¹⁰や環境性能評価¹¹、またディープラーニングに基づく生産最適化のために、DTベースのモデルが提案されています。石油・ガス業界では、システムモデリング、性能評価、およびマルチセンサデータの統合を通じて、監視能力と異常特定能力を向上させるためにDTが活用されています^12,13,14,15。

これらの進展にもかかわらず、現在のDT技術の大部分はモニタリング、モデリング、あるいは異常検知といった個別の機能に集中しており、リアルタイムの同期と制御の統合はほとんどなされていない。多くのシステムでは双方向の相互作用や適応学習機能が不足しており、単方向のデータフローを持つデジタルシャドウとして機能している。動的な状況や複雑な故障パターンにより継続的な学習と制御が必要となる石油・ガスシステムにおいて、この制約は特に重要である。さらに、現在のサイバーフィジカルシステムの手法は主にオフライン分析に限定されており¹⁶、複雑な工業プロセスにおけるリアルタイムのクローズドループ制御はサポートされていない。また、分散センサー間の相互作用を記述するために、現在の手法は主に従来の機械学習やフィルタリング技術に依存しており、時空間グラフニューラルネットワークのような高度なモデルの利用は少ない。さらに、特に石油・ガス生産システムにおいて、非線形な物理ダイナミクスを捉えるためのオペレーター学習技術の利用に関する研究も不足している¹⁷。そのため、既存のDTシステムでは、複雑な時空間的な相関関係をモデリングしつつ、リアルタイムのデータ同化と適応制御を組み込むことができていない。

SWaTやWADIなどのベンチマークデータセットは、産業用制御設定を模倣しているため、サイバー物理システムにおけるDTベースの手法の検証に頻繁に利用されています¹⁸。さらに、DTは、産業プロセス全体を通じた生産性の向上、知識交換の改善、および安全なデータ転送といった、より汎用的なメリットを提供します^19,20。特に石油・ガス産業のような高リスク業界においては^21,22,23、安全管理、リアルタイムモニタリング、および産業オートメーションにおけるその潜在能力がより広く認識されるようになっています。それにもかかわらず、実用的な導入には、モデリング、モニタリング、および制御機能を統合した、より完全で統合的なソリューションが依然として必要です。^24,25で引用されているような評価研究は、運用効率の向上と意思決定プロセスの改善におけるDTの用途と利点を実証することで、その基礎を築きました。表 1に、石油・ガスの製造および制御管理に関する既存手法の比較を示します^26,27,28,29。石油・ガス生産システムの運用リスクが高く、ダイナミクスが複雑であるため、従来のメンテナンス手法では、迅速な異常検知および対応に不十分であることが頻繁にあります^30,31,32,33。近年、安全管理やリモートモニタリングは向上していますが、リアルタイムで継続的にセンシング、学習、および制御ができるインテリジェントシステムが依然として求められています。

本論文では、これらの困難を克服するために、時空間グラフニューラルネットワーク、ニューラル4D-Varデータ同化、物理情報ニューラルオペレーター、および強化学習ベースの制御を組み合わせた、制御認識デジタルツインアーキテクチャを提案します。提案システムにより、データ収集、状態推定、異常検知、および制御最適化のための単一のパイプラインが構築されます。時空間グラフニューラルネットワークが分散センサー間の相互作用を記録する一方で、フーリエニューラルオペレーターやDeep Operator Network（ニューラルオペレーターモデル）などの物理情報ニューラルオペレーターを用いて、複雑な非線形システムダイナミクスを表現します。ニューラルデータ同化によって状態推定の精度が向上し、強化学習によって適応的なクローズドループ制御が可能になります。提案手法の産業現場における汎用性は、サイバー物理システムのベンチマークデータセット（SWaTおよびWADI）およびシミュレーションによる石油・ガスデータを用いて評価することで実証されています。本研究は、モデリング、学習、および制御を統合フレームワークにまとめることで、現在のDTシステムの大きな欠点を克服し、石油・ガス生産システムにおけるリアルタイム監視および最適化のための、スケーラブルで適応性の高いソリューションを提供します。

プロトコル

本研究では、ヒトの被験者および動物実験は行っていません。使用したデータセット（SWaTおよびWADI）は、公開されているベンチマークデータセットであり、機密情報や個人情報は含まれていません。

提案されたアプローチでは、最先端のディープラーニングアルゴリズム、データ同化、および石油・ガス抽出システムとその制御の継続的な評価と最適化の組み込みを可能にする、制御意識型のデジタルツインを構築します。まず、地理的に分散した生産インフラからリアルタイムのマルチモーダルセンサーデータを収集し、生産コンポーネント間の接続性と相互作用を定義するグラフ構造を用いてモデル化します。次に、物理情報を付与したニューラルオペレーターを活用して、非線形多相流を支配する支配方程式のオペレーターレベルの近似を構築し、デジタルツイン内でのスケーラブルなシミュレーションを促進します。本研究において、「Neural operator model」という用語はDeep Operator Networkを指して使用されています。地理的に分散した生産資産間の空間的および時間的依存性を組み込むために、物理情報を付与した表現を用いて時空間グラフニューラルネットワークを学習させます。ノイズが多く空間的に不規則な測定値がある場合でも、デジタルツインのシミュレーション予測をリアルタイムモニタリングデータと一貫させ、同期させるために、neural 4DVARを活用した微分可能なデータ同化アプローチを適用します。リアルタイムで同期されたデジタルツインの上に、リアルタイム制御のために強化学習駆動の最適化ツールを学習させます。最終的に、このデジタルツインは、新たに得られた測定値による継続的な更新に基づき、モデルの継続的なアップデートと長期的な適応的スケーラビリティを備え続けます。

提案されたデジタルツインアーキテクチャは、従来のデジタルツインアーキテクチャにインテリジェンス、適応性、および制御意識を付加することで機能を強化しています。パイプライン全体は、(i) データ収集/グラフ構築、(ii) ニューラルオペレーターおよびST-GNNによる状態推定、(iii) 残差ベース解析による異常検知、(iv) 強化学習ベースの制御、という4つのモジュールで構成されています。

図 1: 制御認識型デジタルツインの全体アーキテクチャ。提案する制御認識型デジタルツインモデルには、石油およびガス生産システムのための異常検知、時空間グラフニューラルネットワーク、ニューラル4D-Var同期、物理情報ニューラルオペレーター、および強化学習ベースのクローズドループ制御がすべて統合されている。こちらのリンクをクリックして、この図の拡大版を表示してください。

図1は、クローズドループのサイバーフィジカルシステムとしての、石油・ガス生産制御を認識したデジタルツインの機能プロセスを示しています。この機能プロセスは石油・ガス生産システムのレベルから始まり、井戸、パイプライン、バルブ、および処理装置に配置されたさまざまなセンサーが、ある瞬間の圧力、流量、温度、および制御に関するデータを収集します。これらの各種センサーはさまざまな形式でデータを収集し、そのデータがデジタルツインに送信されることで、その時点における組織認識の基盤が構築されます。

提案された手法は、以下の手順で構成されています：

データ収集

石油および天然ガスの生産環境から、圧力、流量、温度、および制御入力に関するセンサーからのリアルタイム情報を取得します。堅牢性と汎用性を検証するため、SWATやWADIなどの他の利用可能なベンチマークデータセットを使用します。再現性を確保するため、提案する手法の各ステップは技術的な詳細を用いて規定されています。

[データセットの生成およびシミュレーションについては、補足ファイル (S1) を参照してください。]

産業用制御システムの研究領域における2つの公開ベンチマーク産業用制御システムデータセットを用いて、複雑で高次元の時系列データおよび異常イベントに対する概念的デジタルツインの能力を実証します。使用したデータセットは、Secure Water Treatment (SWaT) および Water Distribution (WADI) のベンチマークデータセットです³⁴。どちらのデータセットも現実的なサイバー物理テストベッドから派生しており、これらのシステムで発生した正常時および異常時/攻撃シナリオの両方における、多変量センサーおよびアクチュエータの動作データが含まれています。SWaTおよびWADIデータセットを使用することで、クロスドメインのテストが可能となり、サイバー物理テストベッドを超えて、場所や時間に依存関係を持つより広範な分散型産業システムへの汎用性を証明できます。データセットの説明は表2に示されています。

提案された研究では、開発したデジタルツイン手法を広範に評価するため、Table 2に示す3つの異なるデータセットを使用します。基本データセットは、石油およびガス製造における生産組織に関するシミュレーションデータで構成されており、システムによって検知された実際の運転測定値ならびに制御および/または故障条件を再現することを目的としています。さらに、縮小スケールの物理システムの例として、Secure Water Treatment (SWaT) データセットを導入します。これには、正常状態および人為的に作成された異常や攻撃を含む、約11日間にわたる51個のセンサーとアクチュエータからの測定値が含まれています。このデータセットは、研究における標準的なベンチマークとして広く認められています。もう一つのセットはWater Distribution (WADI) データセットであり、その測定値は、さまざまな攻撃と正常な運転条件を含む、約123の変数からなる縮小スケールの配水システムの16日間にわたる挙動を捉えたものです。このデータセットには、合計15種類の異なる攻撃が含まれています。上記の2つのデータセットは、主に産業制御システムに見られるような、時間経過に伴う高次元系列および空間的関係を中心としています。

データの前処理

モデルのトレーニング前に、データの品質、一貫性、および再現性を保証するため、すべてのデータセットに対して徹底した前処理手順が行われます。前処理の手順には、欠損値の管理、ノイズフィルタリング、正規化、および特徴量のスケーリングが含まれます。詳細な前処理パイプラインは、補足ファイル (S2)に記載されています。

システムグラフの構築

提案された手法では、石油およびガスの物理的な生産システムがグラフ形式で表現されます。このグラフは、資産間の空間的相関と、時間の経過に伴うシステムの状態のダイナミクスを記述するための数学的構造として機能します（システムグラフ構築の定式化については、補足ファイル (S3)を参照してください）。

ニューラルオペレーターを用いた物理学習

多相流システムの非線形かつ多次元的な特性を正確に表現するため、提案手法ではフーリエニューラルオペレーター（FNO）およびニューラルオペレーターモデルの形式による物理情報に基づいたオペレーターを使用します。ニューラルネットワークを点ごとのマッピングとして表現する従来の手法とは異なり、オペレーターの概念を用いることで、異なる境界条件における汎化が可能になります。物理現象が連続的に発生するため、この特性により、本手法は石油およびガスの生産に使用される多相流システムに対してより適切となります。ニューラルオペレーターの数式については、付録ファイル (S4) を参照してください。

物理学に基づいたニューラルオペレータの妥当性をさらに評価するため、異なる流れの状態および制御入力の下で非線形オペレータマッピングを近似する能力を評価する追加実験を行いました。その結果、これらのオペレータは非線形オペレータのマッピングを正確に近似できることが示され、これにより、特に高次元の状態において、他のソルバーよりも低い計算コストでシステムのダイナミクスを正確にモデル化できる能力が検証されました。オペレータ学習の精度の定量的評価は、予測されたシステム状態とシミュレーションされたシステム状態との間の平均二乗誤差を用いて行われ、異なる流れの条件下で一貫した近似性能が実証されました。

図 2: 物理情報ニューラル演算子（physics-informed neural operator）フレームワークのワークフロー。物理情報ニューラル演算子のワークフローの例を用いて、制御入力、センサーデータ、および境界条件をどのように活用して非線形多相流力学を学習し、演算子レベルでの迅速なシステム予測を生成するかを示している。こちらのリンクをクリックして、この図の拡大版を表示してください。

図2に示すPINOの動作プロセスにおいて、石油・ガスプロセスにおける物理的に複雑な非線形多相流の学習プロセスは、連続的に行われます。これは、モデルが物理的に変化するシステムに関する不均一なデータを受け取る入力層から始まります。この入力は、圧力、温度、速度などの境界条件から、石油・ガスシステム内に配置されたセンサーからのリアルタイムデータ、さらにはバルブの開度、ポンプの速度、チョークの開度などの制御変数までの形式をとります。これらのすべての入力が、物理的に変化するシステムの動作を定義します。

このワークフローにおける主要なコンポーネントは、FNOおよびNeural operatorモデルを通じて実現されるPINOです。提案手法が従来の手法と異なり、明確に区別される点は、入力と出力の間の点ごとのマッピングに重点を置く従来のニューラルネットワークとは対照的に、関数空間間のオペレーターを学習するためにニューラルオペレーターを使用していることです。このレベルにおいて、FNOはフーリエドメインを用いて全域的な空間関係を捉えることができ、Neural operatorモデルはブランチネットワークとトランクネットワークの両方を用いてオペレーター表現を学習することができます。学習されたモデルは物理的な制約によって誘導され、予測されたダイナミクスが多相流において期待通りに動作することを保証します。ニューラルオペレーターの出力は、生産ネットワーク全体の圧力、流量、相分布などの重要な物理変数の連続場解という形式になります。さらに、出力値はセンサーのような特定の点での測定値ではなく、井戸、パイプライン、およびその他の接続ユニットを含むフィールド全体に対して正確です。このように、ニューラルオペレーターの出力により、偏微分方程式を必ずしも解くことなく、その解を近似できる非常に高精度なシミュレーションを迅速に生成することが可能になります。この解が、デジタルツインにおける高速解として機能します。

ニューラルオペレーターの学習プロセスでは、石油およびガスの生産プロセスから抽出された関数レベルおよびシステムレベルの特徴量が用いられます。これらの特徴量は主に、連続場である温度、流量、圧力、および石油・ガス・水の相分率といった石油・ガス生産の物理的状態変数で構成されています。また、特徴量には、入口圧力、流量、バルブ位置などの境界条件も含まれます。異なる運用シナリオにおける応答を予測するため、チョークバルブの位置やポンプ速度などの制御入力も考慮されます。ニューラルオペレーターモデルは、これらの物理的特徴量と制御パラメータを受け取り、それらはトランクネットワークおよびブランチネットワークによって処理されます。これらの特徴量は、物理的プロセスと運用要因の両方に依存する正確なオペレーター関係をニューラルオペレーターが予測することを可能にし、高速なデジタルツイン・シミュレーションを実現するための十分な情報を提供します。

フーリエニューラルオペレーター (FNO)

本研究では、石油・ガス生産システムの非線形多相流ダイナミクスを記述するオペレーター写像を効率的に推定するために、フーリエニューラルオペレーター（FNO）³⁶を使用しています。特定の点集合に対する入出力写像を近似する従来のニューラルネットワークとは異なり、提案されたFNOは関数空間間の写像を近似するため、システムのさまざまなスケールや構成にわたる汎化が可能になります。本手法においてFNOは、システムの状態表現を元の空間からフーリエドメインへ移動させることで機能し、そこではスペクトル畳み込みを効果的に適用して、システム内に存在するグローバルな空間依存性を捉えることができます。具体的には、入力状態表現が高次元の潜在空間に投影され、フーリエ変換、スペクトルカーネル、および活性化関数を用いて変換と活性化が反復的に適用されます。フーリエニューラルオペレーターの数学的評価については、補足ファイル (S4.1)をご参照ください。

ニューラルオペレーターモデル

ニューラルオペレーターモデル³⁷は、石油およびガスの生産動態を定義する非線形物理オペレーターを近似するための、ニューラルオペレータープロセスにおける代替アプローチとして機能します。ニューラルオペレーターモデルの詳細な説明は、補足ファイル(S4.2)に記載されています（アルゴリズム1：デジタルツインのためのニューラルオペレーターによる物理学習アルゴリズムについては、補足ファイルS10を参照してください）。

したがって、ニューラルオペレータを用いた提案手法は、まず入力と出力のペア例を収集することから始まります。ここでの入力は、システム境界値、制御変数、または物理的特性などの形式によるシステム条件の例となります。ニューラルオペレータは反復的なフェーズで学習され、入力関数に対する応答を予測し、その後、平均二乗誤差関数を用いて予測値と実際の応答との間の偏差を改善していきます。このアプローチにより勾配に基づく最適化が行われ、ニューラルオペレータはシステムの挙動を規定するシステム物理オペレータを導き出します。これにより、学習済みシステムは、異なる入力条件に対するシステム応答を効率的に予測することが可能になります。

時空間システムモデリング

提案されたデジタルツインアーキテクチャ³⁸において、石油およびガスの生産プロセスは時空間グラフによって表現されています。これは、分散された生産リソースの接続性と時間的な進化プロセスの概念化を目的として行われました。

石油およびガスの生産プロセスは、井戸、パイプラインインフラ、センサーという観点から、本質的に分散した構造を持っています。このようなプロセスの空間的・時間的関係を効率的に捉えるため、本論文で開発されたデジタルツインアプローチでは、時空間グラフニューラルネットワーク（ST-GNN）を利用しています。この手法では、システムのエンティティをグラフのノードとして表現し、システムの空間的・機能的な関係をグラフのエッジで表します。また、システムの状態遷移や異常を効率的に予測するために、時間的関係はウィンドウ関数を用いて表現されます。Spatio-Temporal Graph Neural Network（ST-GNN）の詳細な説明については、補足ファイル (S5) を参照してください。

Neural 4D VARを用いたデジタルツイン同期

デジタルツインの同期とは、リアルタイムのセンサー読み取り値に基づき、絶えず変化する物理的な石油・ガス生産環境と仮想ツインを継続的に同期させることと定義できます。実際には、デジタルツインは物理システムのダイナミクスの理解に基づき、圧力、流量、温度などのシステム状態について予測を行います。しかし、不確実性やセンサーノイズがあるため、予測値と測定値の間には不可避的に乖離が生じます。このため、予測を実際の測定値と同期させ、予測を物理的な測定値と自己整合させるために、本研究では「データ同化」として知られる別の手法を導入しています。ノイズが多く疎なセンサー観測において、物理プロセスとデジタルツインの同期を確実にするため、ニューラル4D-Varデータ同化法を組み込んでいます。この追加コンポーネントは、同化ウィンドウ内でのセンサー測定値とモデル予測値のギャップを低減することを目的としています。微分可能性を利用することで、ニューラル4D-Var法はエンドツーエンドの学習を提供し、デジタルツインの状態をオンラインで補正することを可能にします（ニューラル4D-Varデータ同化についてはSupplementary File (S6)、アルゴリズム2：ニューラル4D-Varを用いたデジタルツイン同期についてはSupplementary File (S10)、ハイパーパラメータとトレーニングの詳細についてはSupplementary File (S11)を参照してください）。

異常検知および診断

シミュレートされた石油およびガスデータセットにおいて、15%を超える圧力低下、異常な流量スパイク、およびバルブの故障などの制御された摂動を導入することで、アノマリーを注入した。設定されたしきい値を上回るか、あるいはベースラインの統計分布から逸脱したデータポイントをアノマリーとして特定する。アノマリー検出フレームワーク（補足ファイル S7).

クローズドループ最適化および制御

提案されたシステムでは、調和的なDTと最適化のための強化学習モジュールを組み込むことで、クローズドループ制御を実現できます。デジタルツインの部分では、デジタルツインが、異なるバルブの配置、ポンプ速度、または流量制御などのさまざまな条件下における石油およびガス生産プロセスの将来の挙動を予測します。これにより、制御エージェントが、現実世界のプロセスで実際に実行することなく、望ましいアクションの結果を評価できる安全領域を構築することが可能になります。したがって、デジタルツインによる予測を活用することで、強化学習モジュールはプロセスの制約に制限されることなく、プロセスの効率、エネルギー利用、およびプロセスの安定性を含む運用目標を最適化できます。クローズドループ最適化および制御の数学的な記述については、Supplementary File (S8)をご参照ください。

継続的な学習と適応

産業システムは、機器の老朽化、時間的に変化する運転条件、および予期せぬ外乱を伴う非定常な環境に直面しています。これに対応する手法を設計するため、提案されたデジタルツインモデルには、インクリメンタル学習に基づいてニューラルオペレータ、グラフモデル、およびポリシーを更新する継続学習コンポーネントが含まれています。大きな変動や未特定の値が観察された場合、新しいデータを取り込み、デジタルツインモデルに関連するパラメータを調整することで、学習が促進されます。

継続的な学習は、システムが新しい、あるいは斬新な行動パターンを特定するのに役立つため、制御ポリシーと異常検出の両方の堅牢性を向上させています。制御結果の異常からのフィードバックは、モデルの調整と強化学習ポリシーの最適化に適用され、システムが変化し進化し続ける中で効果的な意思決定プロセスを保証します。この適応的な学習プロセスにより、デジタルツインは受動的なオブザーバーから、厳格で効果的な技術インフラ環境を用いて、複雑な石油およびガス生産システムのリアルタイム最適化プロセスを長期間維持できるインテリジェントなシステムへと変化します。Table 3に、モデルアーキテクチャとハイパーパラメータを示します。

継続学習適応の数式については、補足資料 (S9)をご参照ください。

データ前処理モジュールは、生のセンサーデータを、後続の学習モジュールに適した構造化および正規化された入力へと変換する役割を担います。具体的には、過去のセンサー測定値を集計し、ノイズを除去し、特徴量を標準化することで、一貫したデータ品質を確保します。処理されたデータは、境界条件および制御入力とともに、物理情報ニューラルオペレーターモジュール（FNO/ニューラルオペレーターモデル）への入力として提供されます。このモジュールは、潜在的な非線形物理関係を学習し、連続的なシステム状態推定値を生成します。この推定されたシステム状態は、その後、制御最適化モジュールによって利用され、動的な方法で制御戦略を更新します。最適化された制御出力は意思決定モジュールに転送され、そこで結果が解釈され、システム運用に向けた実行可能な決定が生成されます。新たに観測されたデータに基づいてシステムを継続的に更新するためのフィードバックメカニズムが組み込まれており、これにより時間の経過とともに適応性と性能が向上します。明確性と再現性を高めるため、モジュール間の相互作用を示すシステムレベルのアーキテクチャを提示します。提案するシステムは、データ前処理、物理情報ニューラルオペレーターモデリング、制御最適化、および意思決定という4つの主要コンポーネントで構成されています。これらのモジュールは、明確に定義された入出力インターフェースを通じて相互接続されており、デジタルツインシステム全体におけるシームレスなデータフロー、反復学習、および協調動作を保証しています。

図 3: システムレベルの制御最適化アーキテクチャの全体フローチャート。本研究で提案するシステムアーキテクチャの全体像を示している。こちらをクリックして、この図の拡大版を表示してください。

システムアーキテクチャは、石油およびガス生産システムにおける複雑な非線形物理ダイナミクスのモデリングと最適化に使用可能なモデルの全体像を示しています。まず、シミュレーションに使用されるパラメータと物理環境を表す境界条件を含む入力データから始まります。このデータは、最初にデータ前処理モジュールに送られ、そこで正規化と特徴抽出が行われます。これにより、次のモジュールであるニューラルオペレーター、FNO、およびニューラルオペレーターモデルを使用するためのデータ準備が整います。これらのモデルはシステム状態を予測するために使用され、この情報は最適化モジュールを用いてシステムパラメータを最適化するために利用されます。その後、最適化された結果は意思決定モジュールを用いて意思決定に活用されます。このモジュールは、システムの効率性を判断するのに役立ちます。また、このフレームワークには出力結果が含まれており、最適化プロセスの結果が図3に示されています。さらに、出力結果の後にデータが前処理モジュールを通過するフィードバックループも組み込まれています。

結果

物理学的知見に基づいたニューラルオペレーターの妥当性を検証するため、非線形多相流システムのダイナミクスにおけるFNOとニューラルオペレーターモデルの性能を比較する追加の性能試験を実施しました。その結果、FNOは空間的な全体整合性に優れている一方、ニューラルオペレーターモデルは不均一な入力に対する適応性に優れていることが示されました。

提案されたデジタルツインフレームワークの効果的なトレーニングと評価を保証するため、すべての実験はハイパフォーマンスコンピューティング環境で実施されました。実装の主要なプログラミング言語にはPythonを使用しました。物理情報ニューラルオペレーター、時空間グラフニューラルネットワーク、強化学習モジュールなどのディープラーニングモデルの構築には、CUDAによるGPU加速を有効にした普及しているディープラーニングフレームワークと、最適化された深層ニューラルネットワークライブラリを使用しました。大規模な多変量時系列データの処理およびモデルトレーニングを円滑に行うため、マルチコアCPU、専用メモリを搭載した高性能グラフィックスプロセッシングユニット（GPU）、および少なくとも64 GBのシステムメモリを備えたワークステーションでテストを実施しました。ディープラーニングのエコシステムとの互換性を確保するため、オペレーティングシステム環境はLinuxディストリビューションに基づいて構築されました。

モデルの学習にはミニバッチ勾配降下法による最適化を用い、統計的信頼性を確保するために、実験を独立して多数回実施した。提供された計算設定により、同様のハイパフォーマンスコンピューティング環境において、提案した手法を再現することが可能である。

繰り返し試行による結果の妥当性を確保するため、すべての実験は異なるランダム初期化値を用いて複数の独立したランで繰り返されました。シミュレーションデータを用いた実験の場合、モデルを(N)回の独立した試行で実行し、平均性能と標準偏差として値を報告しています。ベンチマークデータであるSWATおよびWADIを用いた実験の場合、評価指標であるMAE、RMSE、F1 Score、およびAUROCをテストデータで算出し、結果の統計的堅牢性を高めるために繰り返し実験の平均を算出しました。これは、モデル性能の信頼できる推定値を得るとともに、性能の向上が単一の実験によるものではないことを確認するために行われました。

それにもかかわらず、物理情報学習、時空間モデリング、およびデータ同化技術を統合デジタルツインシステムに融合させることで、評価条件下において、従来の手法と比較してモニタリング性能の一貫した向上、収束の高速化、およびリアルタイムでの最適な意思決定が実現されることが検証されました。表4に、提案手法のシミュレーション環境を示します。

シミュレーション環境は、産業界におけるリアルタイム運用を正確に再現するようにモデル化されています。ストリームデータ環境を表現するため、センサーからのデータは短い間隔で取得され、堅牢性を検証するために制御された量のノイズが加えられています。異常検知および制御能力をテストするため、さまざまな故障および外乱ケースがシミュレーションされます。また、評価にはスライディングウィンドウ法を用いて、提案されたデジタルツインによる予測や状態同期のための動的な更新、および制御アクションの生成に適した、連続的な形式でのリアルタイム運用を表現しています。

制御認識デジタルツインモデルの設計において、同時評価のために考慮すべきいくつかの評価基準があります。まず、デジタルツインアーキテクチャのモニタリングおよび状態予測能力については、平均絶対誤差（MAE）、平均二乗誤差平方根（RMSE）、および平均絶対パーセント誤差（MAPE）などの回帰基準が用いられます。MAEおよびMAPEは、圧力、流量、温度などの主要なシステムパラメータに対するニューラルオペレータおよび時空間グラフモデルの精度を算出するものであり、石油・ガス生産プロセスで使用されるセンサの連続多変量データに有効です。

次に、状態推定誤差および予測ドリフトの除去に関する、デジタルツイン自体の同期プロセスの評価を行います。状態推定誤差は、実際のシステムの状態で同期されたデジタルツインの状態との間の誤差の程度を測定するものであり、開発されたニューラル4D-Varモジュールが、実際のシステム状態の観測値をどれだけ効果的に利用してシステムダイナミクスの予測を改善しているかを示す指標となります。最後に、予測ドリフトの低減は、ノイズを含む観測値や疎なシステム観測値に対する耐性を備えていることから、提案手法が非同期のデジタルツインシステムに対してどの程度の性能向上を実現したかを示します。

第三に、異常検知および故障診断においては、Precision、Recall、F1、およびROC曲線下面積（AUC）などの分類指向の評価パラメータが用いられます。Recallは、見落としなく実際の故障を正確に診断する能力を指し、Precisionは、認識された異常の妥当性の尺度と見なすことができます。性能はF1スコアを用いることで適切に測定可能です。これらの評価パラメータは、一般的にSWATおよびWADIシステムで採用されています。

表5を用いて、3つのデータセットにおける物理モデルとサイバーモデル間の予測ドリフトの低減に対する、提案されたニューラル4D-Var同期モジュールの有効性を評価する。初期状態誤差は同期なしのモデル間の誤差に対応し、最終状態誤差は同期適用後のモデル誤差に対応する。3つのデータセットすべてにおいて状態誤差が大幅に改善されており、最も高い改善率（72.2%）は、SWAT（63.6%）やWADI（57.7%）に比べて相対的に単純であるため、模擬石油・ガスシステムモデルで得られた。これらの数値は、それぞれの複雑度および「ノイズ」レベルと一致している。模擬石油・ガスデータセットでドリフト低減率が高かった理由は、SWATやWADIのような他の現実的なサイバー物理システムと比較して、模擬環境では一般的にノイズレベルが低く、ランダムな摂動が少なく、システム挙動が明確に定義されているためであると考えられる。ただし、この議論は仮説として扱うべきであり、本研究の結果に基づいて定量的に証明されたものではない。今後の研究では、エントロピー、ノイズ分散、システム次元数などの複雑性の指標を用いて、この論点を検証することが考えられる。

図4： 異なるデータセット（Simulated Oil and Gas、SWAT、およびWADI）におけるパフォーマンスのモニタリング。誤差指標であるMAE、RMSE、およびMAPE (%)を用いて評価した。棒グラフは複数回の独立した試行（n = 5）から得られた平均値を表す。エラーバーは 標準偏差 (SD)、 一方、重ね書きされた細いエラーバーは、 平均標準誤差 (SEM = SD/√n) 各指標のSDおよびSEMの値は、凡例に明記されています。統計的有意性は、ベースラインデータセット（Simulated Oil and Gas）に対する独立t検定を用いて評価され、対応する pp値はバーの上に表示されています。 この図の拡大版を表示するには、ここをクリックしてください。

提案されたデジタルツインモデルを用いて、Simulated Oil and Gas、SWAT、WADIなどのデータセットに対して実施したモニタリングプロセスの性能分析を、MAE、RMSE、およびMAPEの評価指標に基づき図4に示しました。Simulated Oil and Gasデータセットはすべての指標において最も低い誤差値を示しており、これはシステムパラメータが正確にモデル化されたシミュレーション環境におけるシステム状態モニタリングの高い精度を明確に表しています。SWATおよびWADIデータセットでは、シミュレーションシステムには存在しない実システムの複雑さに加え、固有のノイズや非線形性が影響し、システムモニタリングの誤差が徐々に増加していることが分かります。しかしながら、MAEおよびRMSEの誤差値の上昇はわずかであり、このような条件下においてもデジタルツインモデルのシステム性能が安定していることが確認されました。

図 5: F1スコアおよびAUCを用いて評価した、各データセット（Simulated Oil and Gas、SWAT、WADI）における異常検知性能。棒グラフは複数回の独立した試行（n = 5）から得られた平均値を示す。エラーバーはSDを示し、細い重ね合わせのエラーバーはSEM = SD/√nを示す。こちらのリンクをクリックして、この図の拡大版を表示してください。

図5は、3つのデータセットにおけるF1スコアとAUCに基づいた、提案する異常検知システムの評価分析を示している。シミュレーションデータセットであるOil and Gasデータセットは、F1スコアとAUCが最も高く、最高の性能を示している。性能は緩やかに低下しているものの、データセットの複雑さに応じたわずかな変動はあるが、すべてのデータセットにおいて高いAUC性能を達成している。3つのデータセットにおけるF1スコアとAUCの両方の曲線の整合性は、最適な打撃点（striking point）を示しており、これは良好なトレードオフがあることを明らかにしており、提案システムが早期故障検知および侵入検知のための信頼性の高いシステムであることを示している。提案システムは、これら3つのデータセットにおける異常の認識において、極めて優れた性能を発揮している。

提案したデジタルツインモデルを用いて異なるデータセットにおける異常検知を行った際の適合率（Precision）、再現率（Recall）、およびF1スコアを、以下の表6に示す。高い適合率は、特定された関心点の大部分が実際に正解であることを示し、高い再現率は、重要な地点を逃さず、実際の関心点の大部分が正しく特定されていることを示す。石油・ガス分野のデモンストレーション環境では、最高となる0.94のF1スコアを記録した。一方で、SWATおよびWADIにおいて結果がわずかに低くなったのは、センサー間のノイズや相互依存性、および攻撃パターンのレベルが高いためと考えられる。しかしながら、高いF1スコアが得られたことは、異常検知のための堅牢な時空間グラフモデリングと同期が実現していることを裏付けている。

図6: Simulated Oil and Gas、SWAT、およびWADIデータセットにおける異常検知の受信者動作特性（ROC）曲線。性能は曲線下面積（AUC）を用いて評価している。各データセットの凡例には、AUC値の標準偏差（SD）および平均標準誤差（SEM = SD/√n, n = 5）が含まれている。対角線の破線はランダムな分類性能を示す。こちらのリンクをクリックして、この図の拡大版を表示してください。

図6は、しきい値の選択に依存せずに異常検知精度を評価した、3つの異なるデータセットに対する受信者動作特性（ROC）曲線および関連するAUROC値を示している。模擬石油・ガスデータセットのROCでは、曲線が原点に著しく接近しており、AUROC 0.99を達成していることから、実質的に完璧な分類能力があることが確認された。SWATおよびWADIデータセットの場合、AUROC値はそれぞれ0.97および0.95であり、より高い不確実性が存在するにもかかわらず、良好な分類能力を備えていることが示された。ROC曲線がランダム分類を示す直線から大きく離れていることは、開発したデジタルツインモデルが、正確な異常検知において偶然レベルよりも著しく優れた性能を発揮することを保証している。図6の受信者動作特性（ROC）曲線は、3つのデータセットにおける提案したデジタルツインモデルの異常検知性能を示している。AUROCスコアは、模擬石油・ガスデータセットで0.99、SWATデータセットで0.97、WADIデータセットで0.95であり、模擬およびベンチマークの産業制御システムの双方において、本モデルが高い識別能力を持つことを示している。

表7は、デジタルツインと統合されたクローズドループ制御用の強化学習モジュールにおける、累積報酬と制約違反の簡潔な記述を示しています。得られた累積報酬は、運用目標に基づく全体的な制御性能を表し、制約違反は運用中の違反総数を表します。制約違反がゼロの場合の累積報酬が、我々の石油・ガスモデル環境において最高値（96.5）となり、最適な制御および運用性能が確立されています。

SWaTおよびWADIのデータセットでは、累積報酬がわずかに低く、一部に制約違反が見られますが、運用コスト削減における全体的な改善は、提案されたデジタルツインフレームワークが動的な環境における安全かつ最適な制御に有効であることを示しています。

提案された強化学習（RL）駆動型制御システムのリアルタイムな実現可能性を、提案システム内で使用されるすべての計算モジュールのエンドツーエンド推論レイテンシを分析することで評価した。Table 8 の実験結果から、ニューラルオペレーターモジュール（FNO/Neural operator model）および時空間グラフニューラルネットワーク（ST-GNN）モジュールは、プロセス中に使用される時間的特徴量により、中程度の計算オーバーヘッドがあることが示されている。一方で、強化学習ポリシー決定モジュールの計算オーバーヘッドは非常に小さい。提案システム内で使用されるすべてのモジュールの累積エンドツーエンド推論時間は、一般的な産業用制御サイクルの範囲内にある。具体的には、平均エンドツーエンド推論レイテンシは制御サイクルあたり 50–150 milliseconds 以内であることが確認され、これは標準的な産業用制御要件と適合している。このことは、提案システムがリアルタイムおよびニアリアルタイムのアプリケーションに使用可能であることを示している。ただし、エンドツーエンドの推論時間は使用するハードウェアによって変動する場合があることに留意されたい。この点において、提案システムはリアルタイムアプリケーションへの適用が可能である。

物理情報ニューラルオペレーターの有効性をさらに立証するため、FNO、ニューラルオペレーターモデル、および従来の数値ソルバーや現在のオペレーター学習手法などのベースラインアプローチとの定量的比較を実施しました。表9による評価では、各モデルがシステムの動的特性および非線形ダイナミクスをどの程度正確に捉えているかに焦点を当てています。性能評価には、期待されるシステム状態とグランドトゥルース（正解値）との間の平均二乗誤差（MSE）、平均絶対誤差（MAE）、および相対誤差を含む標準的な指標を用いました。実験結果によると、FNOとニューラルオペレーターモデルはともに誤差値を低減させることで、複雑な多相ダイナミクスのモデリングにおいてベースラインアプローチを上回る性能を示しました。具体的には、ニューラルオペレーターモデルは異なる入力条件下での非線形オペレーターマッピングを適切に捉え、FNOは大域的な空間依存性の学習において優れた性能を発揮しました。

提案したニューラルオペレータモデルの有効性は、表9に示された定量的な比較結果によって実証されています。ベースラインとなる数値ソルバーおよび現在のオペレータ学習手法と比較して、FNOとニューラルオペレータモデルの両方が、顕著に低い誤差指標を達成しています。特に、FNOは相対誤差および平均二乗誤差が最も低く、複雑な非線形システムのダイナミクスを捉える優れた能力を示しています。さらに、ニューラルオペレータモデルは、特に異なる条件下での非線形オペレータマッピングの学習において、大幅な改善を示しています。 統計的有意性の評価には独立t検定を用いており、p値が0.05未満である場合は、ベースラインの数値ソルバーに対して大幅な改善があったことを意味します。提案した両モデルで得られたp値は0.001未満であり、観察された性能向上は統計的に有意であり、偶然の結果ではないことを示しています。さらに、提案した手法は推論時間が大幅に短縮されており、リアルタイムのデジタルツインアプリケーションへの適用可能性が示されています。

提案された制御認識デジタルツインモデルの実験的評価をより明確にするため、提案手法と他の手法との包括的な定量的比較を統合的な表形式で提供します。平均絶対誤差（MAE）、平均二乗誤差平方根（RMSE）、F1スコア、曲線下面積（AUC）、予測ドリフトの低減、および累積報酬といったさまざまな性能パラメータを比較することにより、Table 10において、各手法のモニタリング精度、異常検知能力、同期効率、および制御性能を統合的に評価することが可能になります。さらに、実験的評価の透明性を高めるため、バッチサイズ、学習率、エポック数、モジュールごとの学習構成などの学習手順に関するさまざまなパラメータをTable 11に明記しています。加えて、推論レイテンシ、GPU利用率、およびリアルタイム実行可能性の観点から、提案手法のスケーラビリティと計算効率についてもTable 12で評価しています。

提案されたデジタルツインは、表13に示すように、モニタリング精度、同期能力、異常検知、および制御性能の観点から、物理ベースの同期法、LSTM、GNNモデルなどの他の手法と包括的に比較することができます。物理デジタルツインは、静的なアプローチに基づいているため適応性が低く、誤差範囲が大きくなります。LSTMやGNNなどのニューラルアプローチはより優れていますが、物理、同期、および制御を効果的に組み合わせることができていません。提案されたデジタルツインは、RMSEで0.20、ドリフト低減率で72.2%、AUROCで0.94、F1スコアで0.92、累積報酬で95.5という最良の値を達成しました。以上の結果は、物理情報に基づいたニューラルオペレーター、空間と時間のグラフモデリング、ニューラルデータ同化、および強化学習を効果的に組み合わせることで、高精度で堅牢かつ制御を意識したデジタルツインシステムを大幅に改善できることを明確に示しています。

図 7: 提案システムの異なるコンポーネント（FNO、ニューラルオペレーターモデル、GNN、RLモジュール、およびフレームワーク全体）におけるトレーニングコストとリソース利用率。指標にはトレーニング時間（hours）とGPUメモリ使用量（GB）が含まれる。値は複数回の独立したラン（n = 5）の平均結果を示す。エラーバーはSDを示し、重ねられた細いエラーバーはラン間の変動を反映するSEM = SD/√nを示す。この図の拡大版を表示するには、ここをクリックしてください。

図7に、提案手法の学習コストとリソース利用状況を主要コンポーネント別に示します。FNOやNeural operator modelなどのニューラルオペレータモデルは、システム全体と比較して学習コストとGPUメモリ利用率が著しく高いことが分かります。この大幅な学習コストは、オペレータの高い次元学習能力によるものです。グラフニューラルネットワークと強化学習モジュールの学習コストは低いものの、空間的依存性の学習と適応制御において重要な役割を果たしています。提案システムは学習コストが高くなりますが、推論時間は効率的であることが分かりました。単独のLSTMやGNNモデルなどのベースライン手法と比較して、提案システムは複数モジュールの統合により学習コストは増大しますが、表13に示される通り、精度と堅牢性が向上しています。

全体として、提案されたデジタルツイン構造は複数の要素を組み合わせており、それぞれの要素が複雑な工業システムのモデリングと制御における個別の課題に対処しています。計算コストが高くリアルタイムアプリケーションへの適用が困難なことが多い従来の数値ソルバーとは対照的に、Fourier Neural OperatorやDeep Operator Networkなどの物理情報ニューラルオペレーターを用いることで、非線形かつ高次元のシステムダイナミクスをより効果的に捉えることができます。フィードフォワードニューラルネットワークや回帰ニューラルネットワークなどのより単純な機械学習モデルを採用することも可能ですが、これらは通常、基礎的な物理オペレーターを学習できず、異なるシステム構成への汎用性を持ちません。石油・ガス生産システムにおいて本質的に結びついている分散センサー間の依存関係は、時空間グラフニューラルネットワークを用いてモデリングされます。独立した時系列モデルや単純な回帰アーキテクチャのような従来の手法では、空間的な連携を明示的に反映できないため、状態推定や異常検知の精度が低下する可能性があります。強化学習モジュールを利用した適応制御により、動的で不確実な環境下でシステムが運用上の選択を最適化することが可能になります。一方で、ルールベースや静的な制御戦略は硬直的であり、変化するシステム条件や予期せぬ混乱に適応することができません。より複雑さの低い代替案も存在しますが、それらは通常、問題の一部にしか対処できず、精密なモデリング、リアルタイムモニタリング、および適応制御を統合したソリューションを実現するには不十分です。提案されたモジュール設計は、すべての構成要素がシステム全体の性能に寄与することを保証することで、困難な工業環境における精度、拡張性、および堅牢性を向上させます。

研究 / 実務	使用技術	結果	限界点
デジタルツインのためのディープニューラルオペレーター26	サロゲートモデルとしてのDeepONetニューラルオペレーター	リアルタイム予測。物理シミュレーターよりも数桁高速。	慎重なセンサー配置が必要である。また、評価手法の改善が必要である。
DeepONetを用いたバーチャルセンシング対応デジタルツイン27	DTにおけるDeepONet仮想センサ	高速予測、リアルタイム推論、包括的な状態推定	スペクトルバイアスは高周波現象に影響を及ぼすため、ハイブリッドモデリングが必要である
制御におけるDRLとNMPCを組み合わせたデジタルツイン28	深層強化学習 & 非線形MPCの統合	制御性能の向上、トラッキング誤差の低減、および適応的なツイン更新	血管への適用に限定されており、製造システムへのさらなる適用が必要である
産業オートメーションにおけるAI駆動型デジタルツイン29	教師あり/教師なし機械学習 (LSTM, CNN, RFなど)	故障予測における高精度 & モニタリング	主に古典的な機械学習であり、高度な時空間ディープモデルを欠いている。
油 & ガスのデジタルツインに関するレビュー16	各種AI/ML & メカニスティックモデリング	モデリング戦略の包括的な分類学	大部分が概念的であり、高度なディープラーニングを用いた実践的な実装は最小限である

表1：既存の手法を用いた石油・ガス生産および制御管理の比較分析。採用された手法、得られた結果、および重大な欠点に重点を置いた、関連するデジタルツインおよび制御戦略の概要。

データセット名	ドメイン / システム	特徴数	所要時間 & サイズ	データの特性
模擬油 & ガス生成データ	オイル & ガス発生装置（シミュレーション）	多変数（圧力、流量、温度など）	本研究のために設計された	シミュレートされたリアルタイムセンサー & コントロールシグナル、正常 & 故障状態
SWaT	セキュアな水処理CPS（サイバーフィジカルシステム）	約51個のセンサーおよびアクチュエータ	約11日間（45万サンプル以上）	正常 & 異常（36種類の攻撃シナリオ）
WADI	配水サイバー物理システム（CPS）	約123個のセンサおよびアクチュエータ	約16日間（サンプル数 100万件以上）	正常 & 異常（15回攻撃）

表 2: データセットの説明。
データセットの説明。ドメイン、特徴量数、期間、データ属性を含む、本研究のデータセットの概要。

コンポーネント	ハイパーパラメータ	値
ST-GNN	層数	3
ST-GNN	隠れ次元	64
時空間グラフニューラルネットワーク	時間窓の長さ	12
ニューラルオペレーター (FNO)	フーリエモード数	16
ニューラルオペレーター (FNO)	層数	4
DeepONet	分枝ネットワークサイズ	[128, 128]
DeepONet	幹ネットワークサイズ	[128, 128]
データ同化	同化ウィンドウ	10タイムステップ
オプティマイザー	アダム
学習率	0.001
バッチサイズ	64
RLコントローラー	割引率 (γ)	0.99
RLコントローラー	ポリシーネットワーク層	[128, 64]

表 3: モデルアーキテクチャおよびハイパーパラメータ。ニューラルオペレーター、データ同化、強化学習、および時空間グラフニューラルネットワークコンポーネントにおける重要なアーキテクチャ構成とハイパーパラメータ設定。

パラメータ	概要
シミュレーションプラットフォーム	Python (v3.10) および PyTorch (v2.1)
サポートライブラリ	NumPy (v1.24), SciPy (v1.10), Pandas (v1.5)
ハードウェア	NVIDIA RTXシリーズGPU（例：RTX 3090、24 GB VRAM）、Intel Core i7/i9 CPU、32–128 GB RAM
GPU加速	CUDA Toolkit (v11.8), cuDNN (v8.x)
オペレーティングシステム	Linux (Ubuntu 20.04) / Windows 10
開発環境	Jupyter Notebook / VS Code
サンプリング間隔	1～10秒（リアルタイムシミュレーション用に設定可能）
センサーノイズモデリング	ガウスノイズ（σ = 0.01–0.05）
故障シナリオ	センサー故障、フロー乱れ、および制御異常
評価モード	スライディングウィンドウに基づくリアルタイムシミュレーション
コードの利用可能性	再現性を確保するため、リポジトリのリンクは要求に応じて、または出版時に提供されます。

表4: シミュレーション環境。実験セットアップにおけるハードウェア構成、ソフトウェアプラットフォーム、サンプリング間隔、ノイズモデリング、および故障シナリオに関する情報。

データセット	初期状態誤差	最終状態誤差	ドリフト低減率 (%)
擬似オイル & ガス	0.9	0.25	72.2
SWAT	1.1	0.4	63.6
WADI	1.3	0.55	57.7

表 5: デジタルツイン同期の実験結果。データセット全体における予測ドリフトの減少、および初期状態と最終状態の推定誤差を示す同期性能の定量評価。

データセット	適合率	再現率	F1スコア
Simulated Oil & Gas	0.95	0.93	0.94 ± 0.01
SWAT	0.9	0.88	0.92 ± 0.02
WADI	0.87	0.85	0.90 ± 0.02

表 6: 異常検知の性能。シミュレートされた石油・ガス、SWAT、およびWADIデータセットにおける異常検知の結果（適合率、再現率、およびF1スコアを含む）。

データセット	累積報酬 ↑	制約違反 ↓	コスト削減率 (%) ↑
模擬オイル & ガス	96.5 ± 0.8	0	18.7
SWAT	92.1 ± 1.2	1	14.3
WADI	89.7 ± 1.5	2	11.6

表7: クローズドループ制御および最適化の実験結果。クローズドループ制御の性能評価には、累積報酬、制約違反、および運用コストの削減が用いられている。

モジュール	平均潜時 (ms)	概要
ニューラルオペレーター (FNO/DeepONet)	18-25 ms	オペレーター学習とシステム状態予測
時空間グラフニューラルネットワーク	22–35 ms	時空間依存性モデリング
RLポリシーネットワーク	5–10 ms	制御決定の生成
データの前処理	8–12 ms	入力の正規化と特徴量の準備
全制御ループ遅延	55–82 ms	エンドツーエンドの制御動作実行時間

表 8: リアルタイム制御性能の評価。レイテンシ値は、標準的なハードウェア構成において制御ステップごとに測定された平均推論時間を示す。

モデル	MSE ↓	MAE ↓	相対誤差 (%) ↓	推論時間 (ms) ↓	p値 (ベースライン比)
数値ソルバー (ベースライン)	0.0125 ± 0.0012	0.089 ± 0.006	8.75 ± 0.54	120.5 ± 5.2	—
従来のMLモデル	0.0098 ± 0.0009	0.072 ± 0.005	6.42 ± 0.48	85.3 ± 4.7	0.021
既存のオペレーター学習	0.0076 ± 0.0007	0.061 ± 0.004	5.18 ± 0.41	42.8 ± 3.1	0.008
DeepONet (提案手法)	0.0052 ± 0.0005	0.044 ± 0.003	3.67 ± 0.32	18.6 ± 2.4	< 0.001
FNO (提案手法)	0.0047 ± 0.0004	0.039 ± 0.002	3.21 ± 0.28	15.2 ± 2.1	< 0.001

表 9: 予測精度（MSE、MAE、相対誤差）および計算効率（推論時間）におけるベースラインモデルと提案モデルの比較。数値は複数回の試行における平均値 ± 標準偏差として報告されている。値が低いほど、性能が優れていることを示す。p値は、ベースラインの数値ソルバーと比較した統計的有意性を示す。

手法	MAE ↓	RMSE ↓	F1-Score ↑	AUC ↑	ドリフト低減率 (%) ↑	累積報酬 ↑
物理ベースDT	0.45	0.62	0.78	0.82	35.4	65.2
LSTM	0.32	0.48	0.85	0.88	48.7	74.5
GNN	0.28	0.41	0.88	0.91	55.3	81.3
提案フレームワーク	0.2	0.3	0.94	0.97	72.2	95.5

表10: 包括的な性能比較。モニタリング、異常検知、同期、および制御における、提案された制御認識デジタルツインモデルと他のベースライン手法との包括的な性能比較。予測精度については、MAEおよびRMSEの値が低いほど好ましい。異常検知および同期性能については、F1スコア、AUC、ドリフト削減量、および報酬の値が高いほど好ましい。

構成要素	パラメーター	値/設定
ニューラルオペレーター (FNO/DeepONet)	エポック	100–200
	バッチサイズ	32–64
	学習率	0.001
時空間グラフニューラルネットワーク（ST-GNN）	層	3
	隠れユニット	64
	時間的窓	10～20タイムステップ
ニューラル4D-Var	同化ウィンドウ	固定（スライディングウィンドウ）
	最適化	勾配降下法
強化学習	割引率 (γ)	0.95
	報酬ウェイト (α,β,γ)	0.6, 0.25, 0.15
一般的トレーニング	オプティマイザー	アダム
	ハードウェア	NVIDIA RTX GPU

表11： 提案されたデジタルツインフレームワークの各コンポーネントに使用された学習構成および実装パラメータ。本表は、提案システムの各コンポーネントにおける再現性を確保するために使用された、さまざまなハイパーパラメータおよび学習パラメータの概要を示している。

構成要素	トレーニング時間 (時間)	推論レイテンシ (ms)	GPU使用率 (%)	スケーラビリティ
ニューラルオペレーター	6–8	25–40	80–90	高
時空間グラフニューラルネットワーク	3–4	15–25	60–70	高
ニューラル4D-Var	2–3	20–30	65–75	中程度
RLコントローラー	1–2	10–20	50–60	高
システム全体	12–16	<100 ms	70–85	スケーラブルな

表12: 提案されたデジタルツインモデルの異なるコンポーネントにおける計算性能およびスケーラビリティの分析。レイテンシはリアルタイム推論性能を、GPUは計算性能を、スケーラビリティは提案モデルが産業規模のアプリケーションに適用可能である能力を表します。

手法	モニタリングRMSE	同期ドリフト低減率 (%)	AUROC	F1スコア	累積報酬
物理ベースデジタルツイン	0.45	28.6	0.78	0.74	42.3
LSTMベースモデル	0.31	41.2	0.86	0.83	61.7
GNNベースモデル	0.26	53.4	0.91	0.88	74.5
提案デジタルツイン	0.2	72.2	0.94	0.92	95.5

表13: 提案手法と既存手法の比較評価。提案したデジタルツインモデルのモニタリング、同期、異常検知、および制御指標における性能を、ベースライン手法および学習ベースの手法と比較している。

モデル構成	モニタリング誤差 (RMSE)	異常検知 (F1)	制御パフォーマンス（報酬）
フルモデル（全モジュール）	0.2	0.92	95.5
– データ同化なし	0.29	0.88	89.2
– ST-GNNなし	0.27	0.85	87.1
– 強化学習なし	0.21	0.9	78.4
– 継続学習なし	0.23	0.91	92.3

表14: 実験結果およびアブレーション研究。提案されたデジタルツインモデルの各構成要素がシステムの全体的なパフォーマンスにどのように寄与しているかを検証するアブレーション研究。

シナリオ	データの可用性	二乗平均平方根誤差	正確度	観測ギャップの処理
正常	100%	低	高	安定した
部分的喪失	50–70%	中程度	高	ロバスト補間
重度の消失	20–30%	増加した	中程度の	グレースフル・デグラデーション（段階的機能縮小）
センサー故障	地域別 0%	より高い	中程度の	空間汎化

表15: データ不足下でのロバスト性分析。以下の表は、データの可用性が異なる状況における提案されたニューラル4D-Varフレームワークの性能を示している。この表は、センサーデータの部分的損失、高スパース性、および完全損失の状況で誤差がどのように変化するかを示すことで、モデルのロバスト性を提示している。

ディスカッション

本研究の結果は、開発した制御認識型デジタルツインが、幅広いCPSにおけるリアルタイムモニタリング、検出、および制御最適化の向上において有効であることを示している。従来のモニタリングモデルと比較して、PINOとグラフニューラルネットワークの組み合わせは、ノイズのあるシステムにおいてさえ、状態推定の精度においてより効果的である。この結果は、デジタルツインにおける高度な機械学習モデルが、システム内の予測および検出においてより有効であるためである³⁹。

カリキュラム学習や因果表現モデルを採用したデジタルツイン異常検知器などの既存の異常検知ソリューションと比較すると、文献で得られた結果を比較することにより、提案手法が競争力のある性能を達成できることが分かります。例えば、カリキュラム学習を採用したデジタルツインベースの異常検知器は、CPSテストベッドにおいてF1スコアの中程度の向上を可能にしており、これは学習パラダイムが検知プロセスの効率に影響を与える可能性があることを示しています⁴⁰。同時に、現在普及しているデジタルツインモデルは、因果推論による効率的な推論と偽陽性の抑制を強く重視していますが、その代償として因果モデルにおけるモデルの複雑さが増大しています。この点において、提案する同期強化学習アプローチは、単一のバックグラウンド内で効率的な動作最適化と健全な異常検知の両方を提供します。ただし、提案手法と他のアプローチとの比較は、直接的な比較実験の結果に基づくものではなく、他の研究者によって報告された結果のみに基づいていることに留意すべきであり、したがって、この点を考慮して比較を行う必要があります。

さらに、提案されたデジタルツイン内での異常検知の追加と制御アクションの最適化は、CPSに関連する現在の研究開発活動の方向性である、デジタルツインによって提供される意思決定支援および作動能力の重要性が高まっていることを示しています。センサー誘導のための適応的アプローチなど、強化学習を用いたデジタルツインのキャリブレーションに関する研究論文⁴¹では、学習エージェントの使用が、特に動的システムにおいて、精度と適応性の面でさらなるメリットを提供できることが示されています。これは、クローズドループ最適化の有効性の向上に関する知見を裏付けるものであり、知的な統合がCPSの堅牢性と効率性を向上させ得るという仮説を補強しています。前述の通り、この手法は石油・ガス業界向けに設計されていますが、手法の評価は石油・ガス業界のシミュレーションおよび、SWaTやWADIといったCPSのベンチマークデータセットに基づいています。ただし、現在の検証には実際の石油・ガス操業データが含まれていないため、実際の導入シナリオへの適用可能性については慎重に解釈されるべきであることに留意してください。将来的には、提案されたカリキュラム学習法と因果推論に基づくデジタルツイン法の比較評価を、統一された実験フレームワークの下で行い、制御された条件下でこれら2つのアプローチを比較する予定です。ニューラルオペレーター、グラフネットワーク、強化学習を含むこれらのモジュールの統合により計算コストは増加しますが、各モジュールはそれぞれ特定の機能を担っています。アブレーション研究によっても、いずれかのモジュールを削除すると性能が低下することから、提案手法の性能が検証されています。

理論的に、本研究はデジタルツイン構造におけるモデルベース手法とデータ駆動型手法の統合に大きく貢献するものである。物理情報ニューラルオペレーターとニューラルデータ同化技術を状態推定プロセスに統合することで、CPSのためのハイブリッドモデリングアプローチが提供される。これにより、純粋な物理ベースのシミュレーションか、あるいはCPSシミュレーションのためのブラックボックスモデルのいずれかに依存していた従来のデジタルツイン研究における重要な空白が埋められることになる。理論上、本研究は、近年の文献で強調されている、継続的な学習と適応が可能な自律的CPSモデルに向けたデジタルツイン研究の進展と整合するものである。

適用性の観点から、提案されたアプローチは、CPSの運用における安全性と効率性を直接的に向上させ、その結果としてリアルタイムの信頼性、迅速かつ正確な異常特定、および最適な制御ループ戦略を可能にするため、産業現場におけるCPSの機能を大幅に改善することができます。例えば、センサのノイズや非線形システムモデルが伴いやすい産業環境において、同期検出アプローチは誤警報の発生を抑制しつつ、異常の迅速な特定を確実にします。分類指標を主眼とした、侵入検知に関する既存のデジタルツインベースのCPSシステムのソリューションとは異なり、提案アプローチは特定と修正の両アプローチを可能にし、システムのオペレーターがダウンタイムが発生する前にCPSの動作環境を先制的に変更できる能力に寄与します。さらに、ベンチマークデータセットでの性能によって示された、異なる環境へ普遍的に適用できる能力は、スマート製造、エネルギー分野、重要インフラネットワーク管理システムなど、様々なCPS設定における本アプローチの全般的な適用可能性を示しています。

提案されたデジタルツインフレームワークの拡張性は、システムの規模拡大および時間的複雑性への対応能力という観点から評価される。センサー間の相互作用を疎グラフ構造で表現することにより、時空間グラフニューラルネットワークは、センサー数および監視コンポーネント数の増加に伴い効率的にスケールし、計算複雑性の指数関数的な増大を防止する。また、離散化解像度に依存しないオペレーターマッピングを学習する物理情報ニューラルオペレーターを適用することで、拡張性がさらに向上し、計算コストを大幅に増大させることなく大規模システムの効果的なモデリングが可能となる。本フレームワークでは、スライディングウィンドウ法とシーケンスベースの学習を採用することで、計算の実行可能性を損なうことなく、より長い時間範囲における時間的依存性を管理している。さらに、増分学習技術により、完全な再学習を必要とせずに継続的な更新が可能となり、大規模展開におけるオーバーヘッドが削減される。加えて、GPUベースの並列計算を利用することで、高次元データストリームの効率的なリアルタイム処理が実現している。提案されたアーキテクチャは、センサー数と運用時間の両方に合わせてスケールするように設計されており、分散センサーネットワークと継続的な監視要件を備えた石油・ガス生産システムなどの大規模な産業アプリケーションに適している。

アブレーション研究により、ニューラルオペレーター、ST-GNN、データ同化、強化学習といった各モジュールの個別の寄与を理解することが可能になります。本研究では、提案されたデジタルツインアーキテクチャにおける多くの主要な個別モジュールの寄与度を測定するため、アブレーション研究を実施しました。これらの個別モジュールには、データ同化、時空間グラフモデリング、強化学習、および継続学習が含まれます。これらの個別モジュールを除去した際の影響を、モニタリング精度、異常検知、および制御における影響レベルの観点から検証しました。最終的な結果から、全体的なアーキテクチャを利用することが常に最善であることが明確に検証されました。したがって、個別モジュールを除去すると、さまざまなレベルの性能低下が生じ、異なるコンポーネントがシステム全体の性能に異なる形で寄与していることが示されました。データ同化モジュールの除去による影響レベルでは、予測誤差の増加が明確に現れ、異常検知に関する信頼性においても同レベルの影響を及ぼしました。強化学習を除去した際の影響レベルでは、特にモニタリングのみを考慮した場合に、制御が著しく影響を受けることが明確に示されました。

表14に示すアブレーション解析の結果から、各モジュールがシステムの性能に寄与する度合いが異なることがわかります。モニタリング性能の精度に対する寄与度が最も高いのはデータ同化モジュールであり、次いでST-GNNモジュールでした。継続学習モジュールのモニタリング性能の精度に対する寄与度は中程度であり、強化学習モジュールの寄与度は最も低い結果となりましたが、報酬値が95.5から78.4へと低下したことから、制御性能に対する寄与度は強化学習モジュールが最大であることが示されました。

提案したニューラル4D-Varフレームワークの堅牢性を検証するため、観測条件を段階的に悪化させた状態で実験を行いました。具体的には、以下のシナリオを含む実験を実施しました：(i) 欠損データの割合を30%から70%まで変化させた部分的なセンサー故障、(ii) 特定領域における完全なセンサー脱落、(iii) 通信遅延による時間的な疎な観測。実験の結果、提案モデルは中程度のデータ損失下で堅牢な性能を示し、極端なデータ損失下でも緩やかに性能が低下することが確認され、高い適応性が実証されました。表15に堅牢性分析の結果を示します。

本研究の限界は、提案モデルの検証範囲に関連しています。本モデルは石油およびガス生産システムを対象としていますが、提案モデルの評価は、シミュレーションによる石油およびガスの数値、ならびにSWaTやWADIのようなサイバーフィジカル・ベンチマーク問題に基づいている点に留意する必要があります。また、これらのベンチマーク問題は石油およびガスのプロセスに関連していますが、本件における問題は水処理および配送に関連していることも重要な点です。

しかしながら、このような最適な性能にもかかわらず、いくつかの制限が存在します。例えば、特にニューラルオペレーター、グラフネットワーク、および強化学習を統合した場合、一部のデジタルツイン実装と比較して計算複雑性が高くなります。これは、計算能力が限られているプラットフォームへの導入において大きな妨げとなる可能性があります。さらに、本研究ではシミュレーションとベンチマークセットのみに焦点を当てており、運用中のCPSに見られるレベルの多様性を含めることができませんでした。これは、CPSの評価で一般的に採用されているSWaTやWADIなどのシミュレーションセットに見られ、現実世界のCPS環境で経験される多様性が不足している可能性があります。また、適切なセンサーデータの供給を前提としていますが、これは旧式のCPSでは常に実現可能であるとは限りません。本システムの最大の制限は、システム評価のための石油・ガス業界の実データが不足していることです。システムの評価は、石油・ガス業界のシミュレーションおよびCPSのベンチマークデータセット（すなわちSWaTおよびWADI）に基づいているため、石油・ガス業界の現実世界の複雑さを正確に反映していない可能性があります。提案システムは予測能力において有望な性能を示しましたが、学習プロセスにおいて高い計算コストを必要とします。将来的には、低次ニューラルオペレーター、モデルプルーニング、およびハイブリッド物理-MLモデルを含む、提案システムの軽量な代替案を検討することができます。

モデル圧縮技術に関する今後の研究では、実環境におけるスケーラビリティと展開、およびセンサー密度やセンシング情報の品質が大きく異なる既存のレガシーCPSインフラへの統合について検討する必要があります。さらに、本システムを拡張して、オペレーター向けの説明可能性とフィードバックループを組み込むことは、説明可能性の向上および自律的な意思決定への信頼構築に有益であると考えられます。加えて、因果関係デジタルツインと同様に、現在の異常検知プロセスに因果推論手法を統合することで、根本原因分析の改善と誤警報の最小化に寄与する可能性があります。最後に、堅牢性を確保するためには、実際のフィールドデータを用いて本フレームワークを検証する必要があります。

本研究では、複雑なサイバー物理システムにおけるリアルタイムのシステムモニタリング、異常検知、および意思決定能力の向上を目的とした、制御認識可能なインテリジェント・デジタルツインモデルの概念を実証しました。主に受動的なオブザーバーとしてリアルタイムモニタリングに利用される従来のデジタルツインとは異なり、本研究では、物理学に基づいたニューラルオペレーター、時空間グラフニューラルネットワーク、ニューラルデータ同化、強化学習ベースの制御、および学習アルゴリズムを、単一のデジタルツインシステムモデル内に統合することを目指しています。まず、デジタルツインシステムによって、非線形方程式系の物理学に基づいた学習モデル上でシステムのリアルタイムモニタリングが行われます。次に、ニューラル4DVarシステムをデータ同化手順に適用し、物理システムと仮想システムの両方を同期させることで、ノイズやシステム変化に起因するシステム予測のドリフトを最小限に抑えます。異常検知は、システムの空間的関係または時間的なシステム進化のいずれかに関する時空間モデルを用いて実施されます。最終的に、強化学習エージェントが、制御制約を考慮しつつ、システム安全制約の観点から有効な同期済みシステム状態に基づき、システム制御アクションを生成します。模擬的な石油・ガスプロセス、およびベンチマークCPSシナリオであるSWATとWADIを用いて実施した広範な実験により、模擬条件およびベンチマーク条件における提案手法の潜在的な有効性が実証されました。提案モデルにより、モニタリング誤差（MAEおよびRMSE）の低減、優れた異常検知結果、およびクローズドループ制御における一貫した動作が得られました。これらの知見は、デジタルツイン内での異常検知と制御の同期が、運用の信頼性と効率性を維持する上で重要な要因であることを裏付けています。デジタルツイン技術のさらなる発展の可能性としては、人工知能と因果関係による拡張が挙げられます。本フレームワークの評価は、石油・ガス業界の実データを用いて検証することで行われる予定です。

本研究で得られた結論はシミュレーションおよびベンチマークデータセットに基づいたものであるため、実際の石油・ガスシステムへの適用可能性については慎重に解釈する必要があることを強調しておく。

データの利用可能性：

本研究で使用したSWaTおよびWADIデータセットは、公開されている産業制御システムデータセットであり、iTrust Labデータセットリポジトリ（https://itrust.sutd.edu.sg/itrust-labs_datasets/dataset_info/）からアクセス可能です。

材料

この記事で使用された材料の一覧

名前	会社	カタログ番号	コメント
CPU (Intel Xeon / AMD Ryzen)	Intel / AMD	Xeon / Ryzen シリーズ	システムシミュレーション、データの前処理、および一般的な計算に使用
ディープラーニングフレームワーク (PyTorch 2.0)	Meta AI (PyTorch)	https://pytorch.org RRID:SCR_018536	ニューラルネットワークの実装、最適化、および学習パイプラインに使用
深層ニューラルネットワークライブラリ (cuDNN 8.x)	NVIDIA cuDNN	cuDNN 8.x	GPU上のディープラーニング操作用に最適化されたプリミティブを提供
DeepONet	該当なし	https://arxiv.org/abs/1910.03193	スクリプト内では"neural operator model"として参照される
GPU (NVIDIA RTX 4090)	NVIDIA	RTX 4090	ニューラルオペレーター、ST-GNN、および強化学習コンポーネントを含むディープラーニングモデルの学習に使用
GPU加速 (CUDA 11.8)	NVIDIA CUDA	CUDA 11.8 RRID:SCR_018131	高性能なモデル学習のためのGPUベースの並列計算を可能にする
Linuxディストリビューション (バージョン 22.04 LTS)	Canonical (Ubuntu)	https://ubuntu.com/download; RRID:SCR_018127	実験の実行およびモデルのデプロイに使用されるオペレーティングシステム
オペレーティングシステム (Ubuntu 22.04 LTS)	Canonical (Ubuntu)	https://ubuntu.com/download RRID:SCR_018127	すべての実験およびモデルデプロイのための実行環境を提供
最適化アルゴリズム (Adam Optimizer)	該当なし	https://arxiv.org/abs/1412.6980	適応的学習率最適化を用いたモデルの学習に使用
最適化済み深層ニューラルネットワークライブラリ (バージョン 8.x)	NVIDIA cuDNN	https://developer.nvidia.com/cudnn RRID:SCR_018131	深層ニューラルネットワーク操作用に最適化されたGPU加速プリミティブを提供
プログラミング言語 (Python 3.10)	Python Software Foundation	https://www.python.org/downloads RRID:SCR_008394	実装および実験に使用される主要なプログラミング言語
システムメモリ (≥ 64 GB RAM)	Corsair / Kingston / Crucial	https://www.corsair.com / https://www.kingston.com / https://www.crucial.com	高次元の多変量時系列データの処理および大規模なモデル学習をサポート