способ определения температурных поправок при гравиметрических измерениях

Формула изобретения

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

Описание изобретения к патенту

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

Известен принятый за прототип способ определения температурных поправок при измерениях нетермостатированным статическим гравиметром на основе учета инерционного запаздывания показаний гравиметра от изменения температуры (Гуляев П. Ю. Методика учета температурной динамики среды по результатам измерений. Цифровое картографирование, городской кадастр и ГИС: Научн.-техн. сб. по геодезии, аэрокосм. съемкам и картографии. М.: ЦНИИГАиК, 1996. - С. 78-83). Способ включает проведение эталонных испытаний с определением реакции гравиметра на заданные изменения температурного режима, последующий выбор (построение) прогнозной рекуррентной динамической модели температурной инерционности гравиметра, используемой при проведении гравиметрических измерений для вычисления температурных поправок по измеренным до рейса и во время рейса температурам прибора и определенным при эталонных испытаниях коэффициентам инерционности, и последующее внесение вычисленных температурных поправок в результаты измерений.

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

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

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

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

Модель представлена формулами

t_экв,i=(t_i-t₀)-F_i, (1)



q_i=Kt_экв,i, (3)

где t_экв,i - вычисленная эквивалентная температура;

F_i - переходная функция состояния;

- время измерения;

t_i - измеренная температура в момент i;

t₀ - измеренная температура в начальный момент;

- коэффициент инерции, определяемый по результатам эталонирования;

j - индекс (⁺ или ^-) при коэффициенте (_j = ⁺ при значении t_j-t₀0, _j = ^- при значении t_j-t₀<0).
q_i - поправка в измеренное значение силы тяжести за эквивалентную температуру;

К - температурный коэффициент для вычисления поправок за температуру, определяемый по результатам эталонирования;

i - порядковый номер момента, на который вычисляется температура.

При i=1 значение функции F_i-1=0.

Для построения модели используются совокупности G измеренных значений ускорений силы тяжести g()G и совокупности Т_вн значений внутренних температур гравиметра t_вн()T_вн. Значения принадлежат совокупности ряда моментов измерений.

Коэффициенты инерции и температурный коэффициент К определяются по результатам эталонирования гравиметров методом подбора путем подстановки в формулы (1)-(3) различных численных значений этих коэффициентов и выбора из них вероятнейших путем сопоставления временного ряда вычисленных значений поправок q()Q с рядом измеренных значений ускорений силы тяжести g()G. Условие выбора - минимальное расхождение этих рядов.

На фиг. 1-3 представлены графики, иллюстрирующие способ согласно изобретению.

Способ осуществляется следующим образом. До начала работ при проведении метрологических исследований приборов производятся температурные исследования в режиме суточных колебаний температуры, близком по амплитуде и периодичности к режиму предстоящих работ. Измерения начинают спустя 2-3 суток после начала заданного режима. В процессе исследования в течение 3-5 суток непрерывно с интервалами 0,5-1 час снимаются показания гравиметра (ускорение силы тяжести g и внутренняя температура t_вн) и фиксируется время снятия показаний. Полученные временные ряды результатов измерений используются для определения эталонных значений коэффициентов инерции ⁺ и ^- и температурного коэффициента К по приведенным выше формулам (1)-(3). Определение коэффициентов ⁺, ^- и К производится методом подбора. Численные значения коэффициентов подставляются в формулы (1)-(3), вычисляются ряды значений t_экв,i и поправок q_i. Затем временные ряды поправок q()Q, вычисленных при различных значениях ⁺, ^- и К, сравниваются с временным рядом измеренных значений g()G, оценивается точность построенных моделей и выбираются вероятнейшие значения коэффициентов ⁺, ^- и К, которые принимаются за эталонные. Условием выбора служит минимальное расхождение рядов Q и G. Все вычисления производятся на компьютере.

При проведении гравиметрических работ необходимо начинать измерение температуры гравиметра (с фиксацией моментов времени измерений) за несколько часов до начала рейса. Продолжительность предшествующего рейсу периода измерения температуры должна быть не меньше величины (в часах), определенной при эталонировании. Частота измерений до начала рейса и во время рейса должна обеспечивать линейность изменения температуры между моментами измерений. По окончании рейса вычисляются ускорения силы тяжести g и поправки q по формулам (1)-(3), причем вычисление эквивалентных температур начинается с момента t₀ до начала рейса, отстоящего от начала рейса на интервал времени, не меньший величины . После введения поправок q в измеренные значения силы тяжести g смещение нульпункта учитывается линейно по времени по стандартной методике.

Эффективность предлагаемого способа подтверждается приведенными на фиг. 1-3 графиками, полученными при реализации изобретения. Измерения выполнялись гравиметром типа ГНУ. В продолжение всего периода измерений гравиметр неподвижно находился на постаменте. Внутренняя температура гравиметра менялась от 7 до 35^oС. Отсчеты по гравиметру производились с интервалом 1 час. На фиг. 1 представлены графики изменений во времени внутренней температуры t_вн (кривая 1) и показаний g гравиметра (кривая 2). Ускорения силы тяжести g вычислялись по стандартной методике. Судя по графикам (кривые 1 и 2), запаздывание изменения показаний g относительно изменения температуры t_вн составляет более 10 часов. На фиг.2 представлены графики изменения во времени показаний гравиметра g (кривая 3) и вычисленных согласно изобретению значений поправок q (кривая 4). Графики g и q (кривые 3 и 4 на фиг.2) синфазны, то есть запаздывание изменения показаний практически отсутствует.

На фиг.3 представлены графики ошибок g, оставшихся в показаниях гравиметра после учета температурных поправок и поправок за линейное смещение нульпункта. Кривая 5 - график ошибок, оставшихся после исправления согласно изобретению показаний g температурными поправками q, вычисленными по формулам (1)-(3). Для сравнения на фиг.3 приведен график ошибок, оставшихся после учета смещения нульпункта в неисправленных показаниях g (кривая 6), а также график ошибок, оставшихся после учета температурных поправок, вычисленных по внутренней температуре t_вн без учета инерционного запаздывания показаний гравиметра от температуры (кривая 7). Из сравнения между собой кривых 5-7 на фиг. 3, видно, что после введения в измеренные значения g поправок q, вычисленных согласно изобретению, ошибки измерений (кривая 5) уменьшились почти в 10 раз.

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