где
Т – абсолютная температура;
В каталоге продукции BetaTHERM рекомендуется более сложная модель (Модель 2):
где
Т – абсолютная температура;
В измерительных приборах и системах, разрабатываемых и выпускаемых НТЦ «ГеоПАЛС», применяется еще более сложная математическая модель термистора (Модель 3):
где
Т – абсолютная температура;
Для сравнения моделей 1, 2 и 3 были использованы таблицы [6], представляющие собой градуировочные характеристики термисторов из различных материалов, выпускаемых фирмой BetaTHERM. Ниже приводятся результаты анализа для термисторов из материала № 7 (из этого материала изготовлены миниатюрные microchip, используемые в приборах НТЦ «ГеоПАЛС»).
В диапазоне температур (0...100) ºС были методом наименьших квадратов определены постоянные коэффициенты для каждой из моделей 1, 2, 3, определены погрешности аппроксимации во всех точках диапазона температур (через 5 ºС) и максимальная погрешность аппроксимации.
Зависимости погрешности от температуры для моделей 1, 2, 3 представлены на рис. 4.1. Анализ этих графиков показывает, что модель 1 и модель 2 имеют существенные систематические погрешности. В то же время отклонения данных от модели 3 имеют случайный характер, а следовательно дальнейшее усложнение модели нецелесообразно. Максимальная погрешность для модели 1 составляет 0,6 ºС, для модели 2 – 0,25 ºС, для модели 3 – 0,012 ºС, причем погрешность для модели 3 имеет случайный характер и определяется главным образом погрешностью получения исходных экспериментальных данных.
Таким образом, наилучшей математической моделью термистора является модель 3. Аналогичные выводы получены при анализе характеристик термисторов фирмы BetaTHERM, изготавливаемых из других материалов.
Графики зависимости погрешности от температуры для моделей 1, 2, 3
Рис. 4.1
Модель 3 содержит четыре постоянных коэффициента
В зависимости от назначения, состава и структуры измерительной системы для термического каротажа вычисление значения измеряемой температуры по формуле (4.4) с использованием индивидуальных значений параметров
4.2 Математическая модель канала измерения содержания воды в нефти
Как уже отмечалось ранее функция преобразования датчика влажности имеет вид:
Ее можно записать в следующем виде:
где А1 – некоторый постоянный параметр:
СХ находится по формуле:
Или
где А2 – постоянный коэффициент:
Подставляя выражение (4.8) в (4.7) и учитывая выражение (4.5), получим:
Диэлектрическую проницаемость среды ε можно представить следующим образом:
где
α – относительное содержание воды.
Подставляя (4.10) в (4.9) получим:
Обозначим
Удобнее использовать не период повторения импульсов, а их частоту.
Окончательно математическая модель будет выглядеть следующим образом:
Константы а1 и а2определяются путем калибровки в воде и в нефти по следующему алгоритму.
При α=0, т.е. для чистой нефти выражение (4.13) принимает следующий вид:
При α=1, т.е. для воды (4.13) принимает следующий вид:
Для нахождения коэффициентов необходимо решить систему уравнений:
Вычитаем одно уравнение из другого:
Из выражения (4.17) находим а1:
Подставляя значение а1 в любое из уравнений системы, находим а2.
Найдем ориентировочные значения коэффициентов. При чистой воде частота
Тогда:
5 Анализ погрешностей
5.1 Основная погрешность канала измерения температуры
Основная погрешность проявляется при нормальных условиях (температура окружающей среды 20 ºС, относительная влажность воздуха от 45 до 75 %, давление 100 кПа). В ее состав входят погрешность первичного преобразователя и погрешность АЦП.
Погрешность первичного преобразователя включает в себя две составляющие:
1) Погрешность математической модели (из-за нелинейности модели преобразования). При использовании модели (4.4) погрешность в диапазоне температур от 0 до 120 ºС не превышает Δ1=±0,015 ºС.