Для проверки достоверности нейронных операторов с учетом физических закономерностей были проведены дополнительные тесты производительности с целью сравнения эффективности FNO и модели нейронного оператора при моделировании динамики системы нелинейного многофазного потока. Результаты показали, что FNO обеспечивает лучшую глобальную пространственную согласованность, в то время как модель нейронного оператора демонстрирует более высокую адаптивность к гетерогенным входным данным.
Для обеспечения эффективного обучения и оценки предлагаемой структуры цифрового двойника все эксперименты проводились в высокопроизводительной вычислительной среде. В качестве основного языка программирования для реализации использовался Python. Для создания моделей глубокого обучения, таких как физически-информированные нейронные операторы, пространственно-временные графовые нейронные сети и модули обучения с подкреплением, применялся популярный фреймворк глубокого обучения с поддержкой GPU-ускорения через CUDA и оптимизированная библиотека глубоких нейронных сетей. Для облегчения обработки многомерных временных рядов большого объема и обучения моделей испытания проводились на рабочей станции с многоядерным 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). В то время как полнота характеризует способность точно диагностировать реальные ошибки, не пропуская ни одной из них, точность можно рассматривать как меру достоверности распознанной аномалии. Эффективность может быть адекватно оценена с помощью метрики 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-критериев по отношению к базовому набору данных (моделируемая нефть и газ) и соответствующие pЗначения 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-меры и AUC. Столбцы представляют средние значения, полученные в результате нескольких независимых запусков (n = 5). Планки погрешностей указывают Стандартное отклонение, в то время как более тонкие наложенные планки погрешностей обозначают SEM = SD/√n. Пожалуйста, нажмите здесь, чтобы просмотреть увеличенную версию этого рисунка.
Рисунок 5 демонстрирует результаты анализа эффективности предлагаемой системы обнаружения аномалий на основе показателей F1-score и AUC для трех наборов данных. Набор данных «Нефть и газ» (Oil and Gas), являющийся симулированным, показал наилучшие результаты, имея самые высокие значения F1-score и AUC. Несмотря на постепенное снижение эффективности, система демонстрирует высокие показатели AUC для всех наборов данных с незначительными отклонениями в зависимости от сложности данных. Совпадение кривых F1-score и AUC для трех наборов данных указывает на оптимальную точку пересечения, что свидетельствует о хорошем компромиссе и делает предлагаемую систему надежным инструментом для раннего обнаружения неисправностей и вторжений. Предложенная система демонстрирует выдающиеся результаты в распознавании аномалий для всех трех наборов данных.
Показатели точности (precision), полноты (recall) и F1-меры предложенной модели цифрового двойника при обнаружении аномалий на различных наборах данных представлены в Таблица 6 ниже. Хотя высокие значения точности (precision) показывают, что большинство выявленных интересующих точек действительно являются таковыми, высокие значения полноты (recall) свидетельствуют о том, что большая часть фактических интересующих точек определена верно, без пропуска важных объектов. В демонстрационной среде для нефтегазового сектора зафиксирован самый высокий показатель 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.9 и подтверждает практически безупречную способность классификации. В случае наборов данных SWAT и WADI значения AUROC составляют 0.97 и 0.95 соответственно, что свидетельствует о хорошей способности классификации, несмотря на наличие более высоких неопределенностей. Тот факт, что ROC-кривые находятся далеко от линий, представляющих случайную классификацию, гарантирует, что разработанная модель цифрового двойника работает значительно лучше случайного уровня при точном детектировании аномалий. Кривые рабочих характеристик приемника (ROC) на рисунке 6 показывают эффективность предложенной модели цифрового двойника при детектировании аномалий в трех наборах данных. Показатели AUROC составляют 0.9 для симулированного набора данных «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,01, что свидетельствует о том, что наблюдаемое повышение производительности статистически значимо и не является результатом случайности. Кроме того, предложенные методы характеризуются значительно более коротким временем вывода, что подтверждает их применимость в приложениях цифровых двойников в реальном времени.
Для дальнейшего повышения ясности экспериментальной оценки предлагаемой модели цифрового двойника с учетом управления в интегрированном табличном формате представлено исчерпывающее количественное сравнение предлагаемого подхода с другими методами. Такое сравнение различных параметров эффективности, таких как средняя абсолютная ошибка (MAE), среднеквадратичная ошибка (RMSE), F1-мера, площадь под кривой (AUC), снижение дрейфа прогноза и совокупная награда, позволит провести комплексную оценку точности мониторинга, способности к обнаружению аномалий, эффективности синхронизации и качества управления методами в Таблице 10. Кроме того, для повышения прозрачности экспериментальной оценки в Таблице 11 подробно указаны различные параметры процедуры обучения, такие как размер пакета, скорость обучения, количество эпох и конфигурация обучения по модулям. Кроме того, в Таблице 12 оцениваются масштабируемость и вычислительная эффективность предлагаемого подхода с точки зрения задержки вывода, использования GPU и возможности реализации в режиме реального времени.
Предлагаемый цифровой двойник может быть всесторонне сравнен с другими методами, такими как метод синхронизации на основе физических принципов, модели LSTM и GNN, с точки зрения точности мониторинга, способности к синхронизации, обнаружения аномалий и эффективности управления, как показано в Таблице 13. Физические цифровые двойники менее адаптивны и имеют более высокую погрешность из-за использования статических подходов. Нейросетевые подходы, такие как LSTM и GNN, работают лучше, но не обеспечивают эффективного объединения физики, синхронизации и управления. Предлагаемый цифровой двойник достиг наилучших значений: 0,20 для RMSE, 72,2% при снижении дрейфа, 0,94 для AUROC, 0,92 для F1-Score и 95,5 для совокупной награды. Приведенные результаты четко указывают на то, что эффективное объединение физически-информированных нейронных операторов, графового моделирования пространства и времени, нейронной ассимиляции данных и обучения с подкреплением позволяет значительно улучшить точность, надежность и управляемость системы цифрового двойника.

Рисунок 7: Затраты на обучение и использование ресурсов различными компонентами предлагаемой системы, включая FNO, модель нейронного оператора, GNN, модуль RL и общую архитектуру. Показатели включают время обучения (часы) и объем используемой памяти GPU (ГБ). Значения представляют собой средние результаты по нескольким независимым запускам (n = 5). Планки погрешностей указывают на СО (стандартное отклонение), в то время как более тонкие наложенные планки погрешностей представляют собой SEM = SD/√n, что отражает вариабельность между сериями экспериментов. Пожалуйста, нажмите здесь, чтобы просмотреть увеличенную версию этого рисунка.
На рисунке 7 представлены основные компоненты стоимости обучения и использования ресурсов предлагаемого подхода. Можно заметить, что модели нейронных операторов, такие как FNO и нейронная операторная модель, характеризуются значительными затратами на обучение и высоким использованием памяти GPU по сравнению с остальными частями системы. Столь высокая стоимость обучения обусловлена способностью операторов к многомерному обучению. Установлено, что модули графовых нейронных сетей и обучения с подкреплением имеют более низкую стоимость обучения, однако они играют важную роль в изучении пространственных зависимостей и адаптивном управлении. Несмотря на более высокие затраты на обучение предлагаемой системы, время вывода оказалось эффективным. По сравнению с базовыми методами, такими как отдельные модели LSTM или GNN, предлагаемая система влечет за собой более высокую стоимость обучения из-за интеграции нескольких модулей, но обеспечивает повышенную точность и устойчивость, что продемонстрировано в таблице 13.
В целом, предлагаемая структура цифрового двойника объединяет несколько элементов, каждый из которых решает определенную задачу моделирования и управления сложными промышленными системами. В отличие от традиционных численных методов решения, которые зачастую требуют значительных вычислительных затрат и менее пригодны для приложений реального времени, для более эффективного описания нелинейной и многомерной динамики системы используются физически информированные нейронные операторы, такие как Fourier Neural Operator и Deep Operator Network. Могут применяться более простые модели машинного обучения, например, сети прямого распространения или рекуррентные нейронные сети, однако они обычно не позволяют выучить фундаментальные физические операторы и не способны к обобщению для различных конфигураций системы. Зависимости между распределенными датчиками, которые неразрывно связаны в системах добычи нефти и газа, моделируются с помощью пространственно-временной графовой нейронной сети. Традиционные методы, такие как независимые модели временных рядов или простые рекуррентные архитектуры, могут привести к снижению точности оценки состояния и обнаружения аномалий, поскольку в них явным образом не отражены пространственные связи. Адаптивное управление, реализуемое с помощью модуля обучения с подкреплением, позволяет системе оптимизировать эксплуатационные решения в динамических и неопределенных средах. С другой стороны, стратегии управления на основе правил или статические стратегии являются жесткими и неспособны адаптироваться к меняющимся условиям системы или непредвиденным сбоям. Несмотря на наличие менее сложных альтернатив, они обычно затрагивают лишь отдельные аспекты проблемы и недостаточны для создания интегрированного решения, сочетающего в себе точное моделирование, мониторинг в реальном времени и адаптивное управление. Предложенный модульный дизайн повышает точность, масштабируемость и надежность в сложных промышленных условиях, обеспечивая вклад каждого компонента в общую эффективность системы.
| Обучение / Работа | Используемые методы | Результаты | Ограничения |
| Глубокий нейронный оператор для цифрового двойника26 | Нейронный оператор DeepONet в качестве суррогатной модели | Прогнозирование в режиме реального времени; на несколько порядков быстрее физических симуляторов | Требуется тщательное размещение датчиков; методы оценки нуждаются в совершенствовании |
| Цифровой двойник с поддержкой виртуального сенсинга на базе DeepONet27 | Виртуальные датчики DeepONet в цифровых двойниках (DT) | Быстрое прогнозирование; вывод в реальном времени; комплексная оценка состояния | Спектральное смещение влияет на высокочастотные явления; требуется гибридное моделирование |
| Цифровой двойник с использованием глубокого обучения с подкреплением (DRL) и нелинейного прогностического управления (NMPC) в системе регулирования28 | Глубокое обучение с подкреплением & Интеграция нелинейного MPC | Улучшенные характеристики управления; меньшая ошибка слежения и адаптивное обновление цифрового двойника | Применение ограничено сосудами; необходима дальнейшая адаптация для производственных систем |
| Цифровой двойник на базе ИИ в промышленной автоматизации29 | Обучение с учителем/без учителя в машинном обучении (LSTM, CNN, RF и др.) | Высокая точность прогнозирования отказов & мониторинг | Преимущественно классическое машинное обучение; отсутствуют продвинутые пространственно-временные глубокие модели |
| Масло & Обзор цифровых двойников газовых систем16 | Различные методы ИИ/МО & механистическое моделирование | Комплексная таксономия стратегий моделирования | Преимущественно концептуальный характер; минимальная практическая реализация с применением продвинутого глубокого обучения |
Таблица 1: Сравнительный анализ добычи нефти и газа и систем управления по сравнению с существующими методами. Обзор соответствующих стратегий управления и цифровых двойников с акцентом на используемые методы, достигнутые результаты и основные недостатки.
| Название набора данных | Область / Система | Кол-во признаков | Длительность и размер | Характеристики данных |
| Смоделированные данные по добыче нефти и газа | Система добычи нефти и газа (моделирование) | Многомерные (давление, поток, температура и т. д.) | Разработано для данного исследования | Смоделированные сигналы датчиков и управления в реальном времени, нормальные состояния и состояния сбоя |
| SWaT | Киберфизическая система (КФС) безопасной очистки воды | ~51 датчик + исполнительные механизмы | ~1 дней (~450k+ выборок) | Норма и аномалии (36 сценариев атак) |
| WADI | КФС распределения воды | ~123 датчик + исполнительные механизмы | ~16 дней (~1M+ выборок) | Норма и аномалии (15 атак) |
Таблица 2: Описание набора данных.
Описание набора данных. Сводка наборов данных исследования, включая область применения, количество признаков, продолжительность и атрибуты данных.
| Компонент | Гиперпараметр | Значение |
| ST-GNN | Количество слоев | 3 |
| ST-GNN (Spatio-Temporal Graph Neural Network — пространственно-временная графовая нейронная сеть) | Скрытое измерение | 64 |
| ST-GNN (Пространственно-временная графовая нейронная сеть) | Продолжительность временного окна | 12 |
| Нейронный оператор (FNO) | Количество мод Фурье | 16 |
| Нейросетевой оператор (FNO) | Количество слоев | 4 |
| DeepONet | Размер сети ветвления | [128, 128] |
| DeepONet | Размер магистральной сети | [128, 128] |
| Ассимиляция данных | Окно ассимиляции | 10 временных шагов |
| Оптимизатор | Адам | |
| Скорость обучения | 0.001 | |
| Размер серии | 64 | |
| Контроллер на основе обучения с подкреплением | Коэффициент дисконтирования (γ) | 0.99 |
| RL-контроллер | Слои сети стратегии | [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 (v1.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: Результаты эксперимента по синхронизации цифрового двойника. Количественная оценка эффективности синхронизации, демонстрирующая снижение дрейфа прогноза по наборам данных, а также ошибки оценки начального и конечного состояний.
| Набор данных | Точность | Полнота (Recall) | F-мера (F1-score) |
| Имитация нефти & Газ | 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: Экспериментальные результаты управления с замкнутым контуром и оптимизации. Для оценки эффективности управления с замкнутым контуром используются совокупная награда, количество нарушений ограничений и снижение эксплуатационных расходов.
| Модуль | Средняя латентность (мс) | Описание |
| Нейронный оператор (FNO/DeepONet) | 18–25 мс | Обучение операторов и прогнозирование состояния системы |
| ST-GNN | 2–35 мс | Моделирование пространственно-временных зависимостей |
| Сеть политик RL (обучения с подкреплением) | 5–10 мс | Генерация управляющих решений |
| Предобработка данных | 8–12 мс | Нормализация входных данных и подготовка признаков |
| Общая задержка в контуре управления | 5–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, относительная ошибка) и вычислительной эффективности (время вывода). Значения представлены как среднее ± стандартное отклонение по результатам нескольких запусков. Более низкие значения указывают на лучшую производительность. pp-значения указывают на статистическую значимость по сравнению с базовым численным решателем.
| Метод | 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 (Spatio-Temporal Graph Neural Network) — пространственно-временная графовая нейронная сеть | Слои | 3 |
| Скрытые слои | 64 |
| Временное окно | 10–20 временных шагов |
| Нейросетевой 4D-Var | Окно ассимиляции | Фиксированное (скользящее окно) |
| Оптимизация | Градиентный спуск |
| Обучение с подкреплением | Коэффициент дисконтирования (γ) | 0.95 |
| Веса вознаграждения (α,β,γ) | 0.6, 0.25, 0.15 |
| Общее обучение | Оптимизатор | Адам |
| Оборудование | Графический процессор 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 | <10 мс | 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: Экспериментальные результаты и анализ вклада компонентов (абляционный анализ). Анализ вклада каждого элемента предлагаемой модели цифрового двойника в общую эффективность системы.
| Сценарий | Доступность данных | Среднеквадратическая ошибка (RMSE) | Точность | Обработка пропусков в наблюдениях |
| Нормальный | 100% | Низкий | Высокий | Стабильный |
| Частичная потеря | 50–70% | Умеренный | Высокий | Робастная интерполяция |
| Тяжелая степень потери | 20–30% | Увеличенный | Умеренный | Постепенная деградация |
| Отказ датчика | 0% по регионам | Более высокий | Умеренный | Пространственная генерализация |
Таблица 15: Анализ робастности в условиях дефицита данных. В следующей таблице представлены показатели эффективности предлагаемой архитектуры нейронного 4D-Var в различных ситуациях доступности данных. Таблица демонстрирует робастность модели, показывая, как изменяются ошибки при частичной потере, высокой разреженности и полной потере данных с датчиков.