$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Для проверки достоверности нейронных операторов с учетом физических закономерностей были проведены дополнительные тесты производительности для сравнения эффективности FNO и модели нейронного оператора при моделировании динамики системы нелинейного многофазного потока. Результаты показали, что FNO демонстрирует лучшие показатели глобальной пространственной согласованности, в то время как модель нейронного оператора лучше адаптируется к гетерогенным входным данным.
Для обеспечения эффективного обучения и оценки предлагаемой структуры цифрового двойника все эксперименты проводились в высокопроизводительной вычислительной среде. В качестве основного языка программирования при реализации использовался Python. Для создания моделей глубокого обучения, таких как физически-информированные нейронные операторы, пространственно-временные графовые нейронные сети и модули обучения с подкреплением, использовался популярный фреймворк глубокого обучения с ускорением на GPU через CUDA и оптимизированная библиотека глубоких нейронных сетей. Чтобы облегчить обработку многомерных временных рядов большого объема и обучение моделей, испытания проводились на рабочей станции с многоядерным ЦП, высокопроизводительным графическим процессором (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). В то время как полнота характеризует способность точно диагностировать реальные ошибки, не пропуская ни одной из них, точность можно рассматривать как показатель достоверности распознанной аномалии. Эффективность работы может быть адекватно оценена с помощью метрики F1. Данные параметры оценки обычно применяются в системах SWAT и WADI.
Таблица 5 используется для оценки эффективности предлагаемого нейронного модуля синхронизации 4D-Var с целью уменьшения дрейфа прогнозирования между физической и кибермоделью для трех наборов данных. Начальная ошибка состояния соответствует ошибке между моделями без синхронизации, а конечная ошибка состояния — ошибке модели после применения синхронизации. Видно, что для всех трех наборов данных наблюдается значительное улучшение ошибки состояния, при этом наибольшее улучшение (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 для каждой метрики эксплицитно указаны в легенде. Статистическая значимость оценивалась с помощью независимых t-критериев по отношению к базовому набору данных (Simulated Oil and Gas); соответствующие значения p отображены над столбцами. Пожалуйста, нажмите здесь, чтобы просмотреть увеличенную версию этого рисунка.
Анализ эффективности процесса мониторинга, осуществляемого предложенной моделью цифрового двойника на таких наборах данных, как Simulated Oil and Gas, SWAT и WADI, с учетом критериев оценки MAE, RMSE и MAPE, представлен на рисунке 4. Набор данных Simulated Oil and Gas имеет самые низкие показатели ошибки по всем критериям, что наглядно подтверждает точность мониторинга состояния системы в симулированной среде, где параметры системы смоделированы точно. Наборы данных SWAT и WADI демонстрируют ошибки мониторинга системы, которые постепенно увеличиваются из-за сложности реальных систем, а также присущего им шума и нелинейностей, которые отсутствуют в симулированных системах. Тем не менее, незначительное повышение показателей ошибки MAE и RMSE подтверждает стабильность работы модели цифрового двойника в таких условиях.

Рисунок 5: Эффективность обнаружения аномалий в различных наборах данных (Simulated Oil and Gas, SWAT и WADI), оцененная с помощью F1-score и AUC. Столбцы представляют средние значения, полученные в результате нескольких независимых запусков (n = 5). Планки погрешностей указывают SD, а более тонкие наложенные планки погрешностей представляют SEM = SD/√n. Пожалуйста, нажмите здесь, чтобы просмотреть увеличенную версию этого рисунка.
Рисунок 5 демонстрирует результаты анализа эффективности предлагаемой системы обнаружения аномалий на основе показателей F1-score и AUC для трех наборов данных. Набор данных «Oil and Gas» (симулированный набор данных) показывает наилучшие результаты, имея самые высокие значения F1-score и AUC. Несмотря на постепенное снижение производительности, система демонстрирует высокие показатели AUC для всех наборов данных с незначительными колебаниями в зависимости от сложности данных. Совпадение кривых F1-score и AUC для трех наборов данных указывает на оптимальную точку пересечения, что свидетельствует о хорошем компромиссе и делает предлагаемую систему надежным инструментом для раннего обнаружения неисправностей и вторжений. Предложенная система демонстрирует выдающиеся результаты в распознавании аномалий для всех трех наборов данных.
Точность (Precision), полнота (Recall) и F1-меры предложенной модели цифрового двойника при обнаружении аномалий в различных наборах данных представлены в Таблице 6 ниже. Высокие значения точности свидетельствуют о том, что большинство выявленных интересующих точек действительно являются таковыми, а высокие значения полноты показывают, что большинство фактических интересующих точек определены правильно, без пропуска важных объектов. В демонстрационной среде для нефтегазовой отрасли зафиксирована самая высокая F1-мера, равная 0.94, в то время как несколько более низкие результаты на наборах данных SWAT и WADI могут быть обусловлены более высоким уровнем шума и взаимозависимостями датчиков и паттернов атак. Тем не менее, высокие значения F1-меры подтверждают достижение надежного пространственно-временного графового моделирования и синхронизации для обнаружения аномалий.

Рисунок 6: Кривые рабочих характеристик приемника (ROC-кривые) для обнаружения аномалий в наборах данных Simulated Oil and Gas, SWAT и WADI. Эффективность оценивается по площади под кривой (AUC). Стандартное отклонение (SD) и стандартная ошибка среднего (SEM = SD/√n, n = 5) для значений AUC указаны в легенде для каждого набора данных. Диагональная пунктирная линия соответствует эффективности случайной классификации. Нажмите здесь, чтобы просмотреть увеличенную версию этого рисунка.
Рисунок 6 демонстрирует кривые рабочих характеристик приемника (ROC) вместе с соответствующими показателями AUROC для трех различных наборов данных, что позволяет оценить точность обнаружения аномалий независимо от выбора порога. На ROC-кривой для смоделированного набора данных Oil and Gas область проходит очень близко к началу координат, что обеспечивает значение AUROC 0,99 и подтверждает практически безупречную способность к классификации. В случае с наборами данных SWAT и WADI значения AUROC составляют 0,97 и 0,95 соответственно, что свидетельствует о хорошей способности к классификации, несмотря на наличие более высоких неопределенностей. Тот факт, что ROC-кривые значительно удалены от линий, соответствующих случайной классификации, гарантирует, что разработанная модель цифрового двойника работает значительно лучше случайного уровня при точном обнаружении аномалий. Кривые рабочих характеристик приемника (ROC) на Рисунке 6 показывают эффективность предлагаемой модели цифрового двойника при обнаружении аномалий в трех наборах данных. Значения AUROC составляют 0,99 для смоделированного набора данных Oil and Gas, 0,97 для набора данных SWAT и 0,95 для набора данных WADI, что указывает на дискриминационную способность модели как в смоделированных, так и в эталонных промышленных системах управления.
Таблица 7 содержит краткое описание совокупных вознаграждений и нарушений ограничений для модуля обучения с подкреплением для управления в замкнутом контуре, интегрированного с цифровым двойником. Полученное совокупное вознаграждение отражает общую эффективность управления на основе эксплуатационных целей, в то время как нарушения ограничений представляют собой общее количество нарушений в процессе работы. Совокупные вознаграждения при нулевых нарушениях ограничений являются максимальными (96,5) для нашей модели среды нефтегазового сектора, что подтверждает оптимальное управление и эксплуатационные показатели.
Несмотря на то, что наборы данных SWaT и WADI демонстрируют несколько более низкие кумулятивные награды и некоторые нарушения ограничений, общее улучшение экономии эксплуатационных расходов доказывает, что предложенная платформа цифрового двойника остается эффективной для безопасного и оптимального управления в динамических средах.
Возможность реализации предложенной системы управления на основе обучения с подкреплением (RL) в режиме реального времени была оценена путем анализа сквозной задержки вывода всех вычислительных модулей, используемых в данной системе. Экспериментальные данные в Таблице 8 показывают, что модуль нейронного оператора (модель FNO/Neural operator) и модуль пространственно-временной графовой нейронной сети (ST-GNN) имеют умеренные вычислительные затраты из-за использования временных признаков в процессе обработки. С другой стороны, модуль принятия решений на основе политики обучения с подкреплением имеет очень низкие вычислительные затраты. Совокупное время сквозного вывода всех модулей, используемых в предложенной системе, находится в пределах типичных циклов промышленного управления. В частности, средняя сквозная задержка вывода составляет от 50–150 миллисекунд за один цикл управления, что соответствует стандартным требованиям к промышленным системам управления. Это свидетельствует о том, что предложенная система может быть использована в приложениях реального и почти реального времени. Однако следует отметить, что время сквозного вывода может варьироваться в зависимости от используемого оборудования. В этом отношении предложенная система пригодна для применения в режиме реального времени.
Для дальнейшего подтверждения эффективности физически-информированных нейронных операторов проводится количественное сравнение FNO, модели нейронного оператора и базовых подходов, таких как традиционные численные решатели и современные методы обучения операторов. Оценка, представленная в Таблице 9, сосредоточена на том, насколько хорошо каждая модель описывает динамику системы и ее нелинейные характеристики. Для оценки производительности используются стандартные метрики, включая среднюю квадратичную ошибку (MSE), среднюю абсолютную ошибку (MAE) и относительную ошибку между ожидаемыми и истинными состояниями системы. Согласно результатам экспериментов, и FNO, и модель нейронного оператора превосходят базовые подходы при моделировании сложной многофазной динамики, обеспечивая более низкие значения ошибки. В частности, модель нейронного оператора успешно воспроизводит нелинейные отображения операторов при различных входных условиях, в то время как FNO эффективно выявляет глобальные пространственные зависимости.
Эффективность предложенных моделей нейронных операторов подтверждается количественными сравнительными результатами, приведенными в Таблице 9. По сравнению с базовыми численными решателями и современными методами обучения операторов, как FNO, так и модель нейронного оператора достигают заметно более низких показателей ошибки. В частности, FNO имеет наименьшую относительную ошибку и среднеквадратичную ошибку, что демонстрирует ее превосходную способность улавливать динамику сложных нелинейных систем. Кроме того, модель нейронного оператора демонстрирует значительный прирост эффективности, особенно при обучении нелинейным отображениям операторов в различных условиях. Для оценки статистической значимости использовались независимые t-критерии; значения p меньше 0,05 указывают на существенное улучшение по сравнению с базовым численным решателем. Значения p, полученные для обеих предложенных моделей, составляют менее 0,001, что свидетельствует о том, что наблюдаемое повышение производительности является статистически значимым и не является результатом случайности. Более того, предложенные методы демонстрируют значительно меньшее время вывода, что подтверждает их применимость в приложениях цифровых двойников в режиме реального времени.
Для дальнейшего повышения ясности экспериментальной оценки предлагаемой модели цифрового двойника с учетом управления в интегрированном табличном формате представлено всестороннее количественное сравнение предлагаемого подхода с другими методами. Это сравнение различных параметров эффективности, таких как средняя абсолютная ошибка (MAE), среднеквадратичная ошибка (RMSE), F1-мера, площадь под кривой (AUC), снижение дрейфа прогнозирования и кумулятивное вознаграждение, позволит провести комплексную оценку точности мониторинга, способности обнаружения аномалий, эффективности синхронизации и качества управления методами в Таблице 10. Кроме того, для повышения прозрачности экспериментальной оценки в Таблице 11 подробно указаны различные параметры процедуры обучения, такие как размер пакета, скорость обучения, количество эпох и конфигурация обучения по модулям. Далее в Таблице 12 также оцениваются масштабируемость и вычислительная эффективность предлагаемого подхода с точки зрения задержки вывода, загрузки GPU и применимости в режиме реального времени.
Предложенный цифровой двойник может быть всесторонне compared с другими методами, такими как метод синхронизации на основе физических принципов, модели LSTM и GNN, с точки зрения точности мониторинга, способности к синхронизации, обнаружения аномалий и эффективности управления, как показано в Таблице 13. Физические цифровые двойники менее адаптивны и имеют более высокую погрешность из-за использования статических подходов. Нейросетевые подходы, такие как LSTM и GNN, работают лучше, но не обеспечивают эффективного сочетания физики, синхронизации и управления. Предложенный цифровой двойник достиг наилучших показателей: 0.20 для RMSE, 72.2% в снижении дрифта, 0.94 для AUROC, 0.92 для F1-Score и 95.5 для кумулятивного вознаграждения. Приведенные результаты наглядно демонстрируют, что эффективное сочетание нейронных операторов с учетом физики, графового моделирования пространства и времени, нейронного усвоения данных и обучения с подкреплением позволяет значительно улучшить точность, надежность и управляемость системы цифрового двойника.

Рисунок 7: Стоимость обучения и использование ресурсов различными компонентами предлагаемой системы, включая FNO, модель нейронного оператора, GNN, модуль RL и всю структуру в целом. Метрики включают время обучения (часы) и использование памяти GPU (GB). Значения представляют собой средние результаты по нескольким независимым запускам (n = 5). Планки погрешностей указывают SD, в то время как более тонкие наложенные планки погрешностей представляют SEM = SD/√n, отражая вариабельность между запусками. Пожалуйста, нажмите здесь, чтобы просмотреть увеличенную версию этого рисунка.
На Рисунке 7 представлены основные компоненты затрат на обучение и использования ресурсов предлагаемого подхода. Можно заметить, что модели нейронных операторов, такие как FNO и нейронная модель оператора, имеют значительные затраты на обучение и потребление памяти GPU по сравнению с остальными частями системы. Высокие затраты на обучение обусловлены способностью операторов к многомерному обучению. Было установлено, что модули графовых нейронных сетей и обучения с подкреплением характеризуются более низкими затратами на обучение, однако они имеют важное значение для изучения пространственных зависимостей и адаптивного управления. Несмотря на то, что предлагаемая система требует больших затрат на обучение, время вывода оказывается эффективным. По сравнению с базовыми методами, такими как отдельные модели LSTM или GNN, предлагаемая система влечет за собой более высокие затраты на обучение из-за интеграции нескольких модулей, но обеспечивает повышенную точность и устойчивость, что продемонстрировано в Таблице 13.
В целом, предлагаемая структура цифрового двойника объединяет несколько элементов, каждый из которых решает определенную задачу в моделировании и управлении сложными промышленными системами. В отличие от традиционных численных решателей, которые зачастую требуют больших вычислительных затрат и менее пригодны для приложений реального времени, физически-информированные нейронные операторы, такие как Fourier Neural Operator и Deep Operator Network, используются для более эффективного захвата нелинейной и многомерной динамики системы. Могут быть применены более простые модели машинного обучения, такие как нейронные сети прямого или рекуррентного распространения, однако они обычно не осваивают фундаментальные физические операторы и не способны к обобщению для различных конфигураций системы. Зависимости между распределенными датчиками, которые неразрывно связаны в системах добычи нефти и газа, моделируются с помощью пространственно-временной графовой нейронной сети. Традиционные методы, такие как независимые модели временных рядов или простые рекуррентные архитектуры, могут привести к снижению точности оценки состояния и обнаружения аномалий, поскольку они не отражают пространственные связи явным образом. Адаптивное управление, использующее модуль обучения с подкреплением, позволяет системе оптимизировать эксплуатационные решения в динамических и неопределенных средах. С другой стороны, стратегии управления на основе правил или статические стратегии являются ригидными и не способны адаптироваться к изменяющимся условиям системы или непредвиденным сбоям. Хотя существуют менее сложные альтернативы, они обычно решают лишь отдельные части проблемы и недостаточны для создания интегрированного решения, сочетающего точное моделирование, мониторинг в реальном времени и адаптивное управление. Предложенная модульная архитектура повышает точность, масштабируемость и устойчивость в сложных промышленных условиях, обеспечивая вклад каждого компонента в общую эффективность системы.
| Обучение / Работа | Используемые методы | Результаты | Ограничения |
| Глубокий нейросетевой оператор для цифрового двойника26 | Нейросетевой оператор DeepONet в качестве суррогатной модели | Прогнозирование в режиме реального времени; на несколько порядков быстрее физических симуляторов | Требуется тщательное размещение датчиков; методы оценки нуждаются в совершенствовании |
| Цифровой двойник с поддержкой виртуального сенсоринга на базе DeepONet27 | Виртуальные датчики DeepONet в цифровых двойниках | Быстрое прогнозирование; вывод в реальном времени; комплексная оценка состояния | Спектральный сдвиг влияет на высокочастотные явления; необходимо гибридное моделирование |
| Цифровой двойник с применением DRL и NMPC в системе управления28 | Глубокое обучение с подкреплением & Интеграция нелинейного MPC | Улучшенные рабочие характеристики управления; снижение ошибки слежения и адаптивное обновление цифрового двойника | Применение ограничено сосудами; необходима дальнейшая апробация на производственных системах |
| Цифровой двойник на базе ИИ в промышленной автоматизации29 | Обучаемое/необучаемое машинное обучение (LSTM, CNN, RF и т. д.) | Высокая точность прогнозирования неисправностей & мониторинг | Преимущественно классическое машинное обучение; отсутствуют продвинутые пространственно-временные глубокие модели |
| Масло & Обзор цифровых двойников газовых систем16 | Различные методы ИИ/МО & механизм моделирования | Комплексная таксономия стратегий моделирования | Преимущественно концептуальный характер; минимальное количество практических реализаций с использованием продвинутого глубокого обучения (DL) |
Таблица 1: Сравнительный анализ добычи нефти и газа и управления регулированием с использованием существующих методов. Обзор соответствующих стратегий цифровых двойников и управления с акцентом на применяемые методы, достигнутые результаты и существенные недостатки.
| Название набора данных | Область / Система | Кол-во признаков | Продолжительность & Размер | Характеристики данных |
| Моделирующее масло & Данные о выходе газа | Масло & Система генерации газа (симуляция) | Многофакторный анализ (давление, поток, температура и т. д.) | Разработано специально для данного исследования | Имитация датчика реального времени & контрольные сигналы, норма & состояния неисправности |
| SWaT | Защищенная киберфизическая система (КФС) очистки воды | ~51 датчиков и исполнительных механизмов | ~11 дней (~450 тыс. свыше образцов) | Нормальный & аномалия (36 сценария атаки) |
| WADI | Киберфизическая система распределения воды | ~123 датчиков и актуаторов | ~16 дней (~1 млн+ образцов) | Нормальный & аномалия (15 атак) |
Таблица 2: Описание набора данных.
Описание набора данных. Сводка наборов данных исследования, включая область, количество признаков, продолжительность и атрибуты данных.
| Компонент | Гиперпараметр | Значение |
| ST-GNN (Spatio-Temporal Graph Neural Network — пространственно-временная графовая нейронная сеть) | Количество слоев | 3 |
| ST-GNN | Скрытое измерение | 64 |
| ST-GNN (пространственно-временная графовая нейронная сеть) | Длительность временного окна | 12 |
| Нейросетевой оператор (FNO) | Количество мод Фурье | 16 |
| Нейронный оператор (FNO) | Количество слоев | 4 |
| DeepONet (Глубокая операторная сеть) | Размер разветвленной сети | [128, 128] |
| DeepONet | Размер магистральной сети | [128, 128] |
| Ассимиляция данных | Окно ассимиляции | 10 временных шагов |
| Оптимизатор | Адам | |
| Скорость обучения | 0.001 | |
| Размер партии | 64 | |
| RL-контроллер | Коэффициент дисконтирования (γ) | 0.99 |
| Контроллер с обучением с подкреплением | Слои сети политик | [128, 64] |
Таблица 3: Архитектура модели и гиперпараметры. Важные конфигурации архитектуры и настройки гиперпараметров для нейронных операторов, ассимиляции данных, обучения с подкреплением и компонентов пространственно-временной графовой нейронной сети.
| Параметр | Описание |
| Симуляционная платформа | Python (v3.10) с PyTorch (v2.1) |
| Вспомогательные библиотеки | NumPy (v1.24), SciPy (v1.10), Pandas (v1.5) |
| Оборудование | ГП NVIDIA серии RTX (например, RTX 3090, 24 ГБ VRAM), ЦП Intel Core i7/i9, 32–128 ГБ ОЗУ |
| 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, включая точность (precision), полноту (recall) и F1-меру.
| Набор данных | Кумулятивное вознаграждение ↑ | Нарушения ограничений ↓ | Снижение затрат (%) ↑ |
| Моделированное масло & Газ | 96.5 ± 0.8 | 0 | 18.7 |
| SWAT (сокращение от Surface Water Assessment Tool) — инструмент оценки поверхностных вод | 92.1 ± 1.2 | 1 | 14.3 |
| WADI | 89.7 ± 1.5 | 2 | 11.6 |
Таблица 7: Экспериментальные результаты управления с замкнутым контуром и оптимизации. Для оценки эффективности управления с замкнутым контуром используются совокупная награда, количество нарушений ограничений и снижение эксплуатационных расходов.
| Модуль | Средняя латентность (мс) | Описание |
| Нейросетевой оператор (FNO/DeepONet) | 18–25 мс | Обучение операторов и прогнозирование состояния системы |
| ST-GNN | 22–35 мс | Моделирование пространственно-временных зависимостей |
| Сеть стратегий обучения с подкреплением | 5–10 мс | Генерация управляющих решений |
| Предобработка данных | 8–12 мс | Нормализация входных данных и подготовка признаков |
| Общая задержка контура управления | 55–82 мс | Время выполнения управляющего воздействия от начала до конца |
Таблица 8: Оценка производительности управления в режиме реального времени. Значения задержки представляют собой среднее время вывода, измеренное за один шаг управления при стандартной конфигурации оборудования.
| Модель | MSE ↓ | MAE ↓ | Относительная погрешность (%) ↓ | Время вывода (мс) ↓ | p-значение (относительно исходного уровня) |
| Численный решатель (базовый вариант) | 0.0125 ± 0.0012 | 0.089 ± 0.006 | 8.75 ± 0.54 | 120.5 ± 5.2 | — |
| Традиционная модель машинного обучения | 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-мера ↑ | AUC ↑ | Снижение дрейфа (%) ↑ | Кумулятивное вознаграждение ↑ |
| Физически обоснованный цифровой двойник | 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 |
| Графовые нейронные сети | 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-score, 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 |
| Общий тренинг | Оптимизатор | Адам |
| Оборудование | GPU NVIDIA RTX |
Таблица 11: Конфигурация обучения и параметры реализации, использованные для различных компонентов предлагаемой структуры цифрового двойника. В таблице представлен обзор различных гиперпараметров и параметров обучения, использованных для обеспечения воспроизводимости предлагаемой системы для разных компонентов.
| Компонент | Время обучения (ч) | Задержка вывода (мс) | Загрузка GPU (%) | Масштабируемость |
| Нейронные операторы | 6–8 | 25–40 | 80–90 | Высокая |
| ST-GNN | 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 (площадь под ROC-кривой) | F-мера (F1-score) | Кумулятивное вознаграждение |
| Цифровой двойник на основе физических моделей | 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: Экспериментальные результаты и анализ влияния компонентов (ablation study). Анализ влияния компонентов, исследующий вклад каждой части предлагаемой модели цифрового двойника в общую эффективность системы.
| Сценарий | Доступность данных | Среднеквадратическая ошибка (RMSE) | Точность | Обработка пропусков в наблюдениях |
| Нормальный | 100% | Низкий | Высокий | Стабильный |
| Частичная потеря | 50–70% | Умеренный | Высокий | Робастная интерполяция |
| Тяжелая утрата | 20–30% | Увеличенный | Умеренный | Постепенная деградация |
| Отказ датчика | 0% по регионам | Выше | Умеренный | Пространственная генерализация |
Таблица 15: Анализ робастности в условиях дефицита данных. В следующей таблице представлены показатели эффективности предложенной нейросетевой структуры 4D-Var при различных сценариях доступности данных. Таблица демонстрирует робастность модели, показывая, как изменяются ошибки в ситуациях частичной потери, высокой разреженности и полной потери данных с датчиков.