Ошибка сингулярной матрицы означает, что решатель не может определить единственное решение стационарного состояния на основе предоставленных уравнений. В модели материального баланса пилотной установки это обычно происходит потому, что уравнения баланса линейно зависимы, отсутствует важная технологическая спецификация или число независимых уравнений не соответствует числу неизвестных переменных потоков. Практическое решение заключается в определении числа степеней свободы модели, проверке ранга матрицы коэффициентов, удалении избыточных уравнений и добавлении только действительно независимых ограничений.

Решателю необходимо достаточно независимой информации, чтобы определить каждую неизвестную величину. Сингулярная матрица указывает на то, что по крайней мере одно уравнение дублирует уже имеющуюся информацию либо в модели отсутствует необходимая спецификация, например разделение потока, состав, конверсия или стехиометрическая зависимость.

Почему матрица становится сингулярной

Избыточные общие балансы

Общий материальный баланс может быть прямой суммой отдельных компонентных балансов. Например, если для воды, азота и кислорода уже записаны отдельные балансы, добавление общего материального баланса может не дать нового независимого уравнения.

В результирующей матрице содержатся уравнения, описывающие одно и то же ограничение в разных формах. Поэтому её ранг оказывается ниже числа неизвестных.

Полные компонентные балансы и ограничения состава

Состав потоков является ещё одним распространённым источником зависимости. Мольные доли удовлетворяют условию:

[ \sum_i x_i = 1 ]

Если все доли компонентов рассматриваются как неизвестные, одна доля обычно определяется через остальные и условие суммирования. Включение каждого компонентного баланса, каждой доли как независимой переменной и нескольких общих балансов может привести к переопределённой или зависимой системе уравнений.

Отсутствующие технологические спецификации

Модель также может стать сингулярной из-за недостаточной определённости. Одни только материальные балансы могут не определять, как поток разделяется между выходами или какой состав выходит из технологического блока.

К типичным отсутствующим спецификациям относятся:

Расход потока или подачи
Коэффициент разделения потока
Состав выходного потока
Спецификация по степени извлечения или влажности
Конверсия реакции
Стехиометрическая зависимость
Условие фазового равновесия или насыщения

Без такой информации балансам могут соответствовать несколько стационарных решений.

Повторяющиеся балансы блока или установки в целом

Если несколько технологических блоков моделируются совместно, общий баланс установки может просто представлять собой сумму балансов, уже записанных для отдельных блоков. Он может быть полезен в качестве проверочного, но его не следует включать в основное решение, если он не добавляет независимой информации.

Как диагностировать модель

Подсчёт неизвестных и уравнений

Явно перечислите все неизвестные, включая:

Общие расходы потоков
Расходы компонентов
Мольные или массовые доли
Степени протекания реакций
Доли разделения потоков
Неизвестные рабочие параметры

Затем подсчитайте количество независимых уравнений, доступных для их определения. В корректно заданной квадратной системе число независимых уравнений равно числу неизвестных, что даёт нулевое число степеней свободы.

Основная зависимость имеет вид:

[ \text{Число степеней свободы}

\text{число неизвестных}

\text{число независимых уравнений} ]

Положительное значение указывает на недостаточно заданную модель. Отрицательное значение означает, что модель может содержать избыточные уравнения или противоречивые спецификации.

Проверка ранга матрицы

Для системы, записанной в виде:

[ A\mathbf{x}=\mathbf{b} ]

матрица коэффициентов (A) должна содержать достаточно независимых столбцов для определения вектора неизвестных (\mathbf{x}). Для квадратной системы единственное решение требует:

[ \operatorname{rank}(A)=\text{число неизвестных} ]

В MATLAB проверьте матрицу с помощью:

rank(A)
size(A)

В Python используйте:

import numpy as np

np.linalg.matrix_rank(A)
A.shape

Ранг, меньший числа неизвестных, подтверждает, что модель не содержит достаточно независимой информации для получения единственного решения.

Анализ приведённых уравнений

Один только ранг выявляет проблему, но не объясняет её. Проанализируйте уравнения после исключения переменных или приведения строк, чтобы найти зависимости, такие как:

[ \text{общий баланс}

\text{баланс воды} + \text{баланс азота} + \text{баланс кислорода} ]

Также проверьте, не является ли уравнение состава простой суммой определений компонентных расходов и не задают ли две технологические спецификации одну и ту же зависимость для разделения потока или конверсии.

Практический подход к решению балансов

Согласованное преобразование измеренных расходов

Данные пилотной установки часто включают объёмные, массовые и молярные расходы. Прежде чем формировать матрицу, приведите их к единой основе.

Для жидкого потока объёмный расход можно преобразовать в молярный расход с использованием плотности и молекулярной массы:

[ \dot n

\frac{\dot V \rho}{M} ]

Несогласованные единицы не обязательно вызывают математическую сингулярность, но могут привести к неправильным коэффициентам и скрыть реальную структуру модели.

Сначала решайте компоненты, определяемые напрямую

Если компонент поступает через один поток и выходит через один поток без реакции или накопления, его баланс может непосредственно определить неизвестный расход. Вода в не реагирующей операции испарения или увлажнения является распространённым примером.

Этот аналитический этап помогает выявить, какие переменные действительно определяются, и уменьшить размер численной системы.

Осторожное применение общего баланса

После определения компонентных расходов используйте общий массовый или молярный баланс для определения оставшихся величин потоков. Считайте общий баланс независимым уравнением только в том случае, если он уже не следует из компонентных балансов.

Для стационарного нерегирующего блока:

[ \text{общий вход}=\text{общий выход} ]

Баланс должен включать каждый входящий и выходящий поток на одной и той же массовой или молярной основе.

Использование соотношений состава без избыточного определения переменных

Если мольные доли выходного потока неизвестны, стратегически используйте компонентные балансы и соотношение суммирования долей. Для (N) компонентов часто предпочтительно решать для (N-1) долей, а последнюю долю вычислять по формуле:

[ x_N=1-\sum_{i=1}^{N-1}x_i ]

Это позволяет не рассматривать зависимую долю как дополнительную независимую переменную.

Понимание компромиссов

Удаление уравнений может скрыть недостающую физику

Удаление избыточного уравнения уместно, если оно дублирует уже имеющуюся информацию. Однако это не следует делать, если уравнение представляет отсутствующее физическое ограничение, которое случайно было записано в зависимой форме.

Например, замена общего баланса спецификацией коэффициента разделения может устранить дефицит ранга и одновременно предоставить информацию, которая действительно необходима модели.

Добавление уравнений может привести к переопределённости модели

Добавление всех доступных измеренных зависимостей может создать переопределённую систему. Если измерения несовместимы из-за ошибок датчиков, уравнения могут не иметь точного решения, даже если они независимы.

Для согласования данных используйте подходящий метод наименьших квадратов или ограниченной оптимизации вместо того, чтобы заставлять обычный решатель квадратных систем удовлетворять противоречивым измерениям.

Дефицит ранга отличается от плохой обусловленности

Действительно сингулярная матрица содержит недостаточно независимой информации. Почти сингулярная или плохо обусловленная матрица может формально иметь полный ранг, но при этом давать крайне чувствительные или ненадёжные результаты, поскольку уравнения почти дублируют друг друга или переменные имеют существенно разный масштаб.

Используйте подходящее масштабирование и проверяйте число обусловленности матрицы, если решатель сообщает о нестабильности, а не о точной сингулярности.

Отрицательные результаты часто указывают на отдельную ошибку моделирования

Отрицательные расходы потоков или составы обычно нельзя исправить заменой линейного решателя. Они часто указывают на перепутанные знаки входа/выхода, пропущенные потоки, неправильные единицы или недопустимую технологическую спецификацию.

После устранения проблем с рангом повторно проверьте сохранение вещества, знаки и физические ограничения, например неотрицательность расходов и значения составов от нуля до единицы.

Как сделать правильный выбор в соответствии с вашей целью

Перед повторным запуском решателя стационарного состояния выполните следующие проверки:

Если ваша основная цель — найти единственное математическое решение: подсчитайте неизвестные, определите ранг матрицы коэффициентов и убедитесь, что модель имеет нулевое число степеней свободы и полный столбцовый ранг.
Если ваша основная цель — удалить избыточные балансы: сохраните независимые компонентные балансы и технологические спецификации, а общие балансы используйте как проверочные уравнения, если они являются линейными комбинациями.
Если ваша основная цель — моделировать разделение потока или состав выходного потока: добавьте явную независимую спецификацию коэффициента разделения, степени извлечения, состава или фазового равновесия.
Если ваша основная цель — согласовать измеренные данные пилотной установки: отделите точные физические балансы от зашумлённых измерений и используйте метод ограниченных наименьших квадратов, если система переопределена или несовместима.
Если ваша основная цель — исправить нефизичные результаты: после проверки ранга матрицы проверьте знаки потоков, пропущенные технологические члены, единицы измерения, стехиометрию и ограничения состава.

Сингулярная матрица — это диагностический сигнал о том, что информационную структуру модели необходимо скорректировать до того, как можно будет доверять её численному решению.

Сводная таблица:

Причина	Диагностическая проверка	Рекомендуемое решение
Избыточные балансы	Сравнить зависимости между общими и компонентными балансами	Удалить дублирующие уравнения или оставить общие балансы для проверки
Отсутствующая спецификация	Проверить число степеней свободы и неизвестные переменные потоков	Добавить независимое условие разделения потока, степени извлечения, состава, конверсии или равновесия
Зависимые переменные состава	Проверить ограничения суммирования долей	Решать для N-1 долей и вычислять последнюю долю по соотношению суммы
Несогласованные единицы	Проверить массовые, молярные и объёмные основы расходов	Привести все данные к единой основе перед построением матрицы
Плохая обусловленность	Проверить масштабирование и число обусловленности	Изменить масштаб переменных и коэффициентов; отличать почти сингулярность от истинного дефицита ранга
Противоречивые измерения	Проверить, имеет ли система измеренных уравнений точное решение	Использовать согласование данных методом ограниченных наименьших квадратов вместо базового решателя квадратной системы
Нефизичные результаты	Проверить знаки, пропущенные потоки, стехиометрию и ограничения	Исправить технологическую модель и проверить неотрицательность расходов и допустимость составов

Планируете строительство или модернизацию пилотной установки для обучения, разработки процессов или промышленных исследований? Свяжитесь с LABPARK, чтобы обсудить учебные и профессиональные пилотные установки для единичных операций, адаптированные к вашим потребностям. LABPARK предоставляет модульные системы для химической инженерии, биопроцессов и биотехнологий, а также экологических технологий и очистки воды. Наша команда поможет университетам, исследовательским институтам и предприятиям выбрать подходящие балансы, контрольно-измерительные приборы, технологические спецификации и экспериментальную конфигурацию для надёжной эксплуатации, более эффективного обучения и воспроизводимых результатов исследований.