В настоящее время в области информационно-измерительных систем, предназначенных для контроля пространственного положения скважины, намечаются тенденции к разработке гибридных инклинометрических систем (ИС) [1, 5, 6].
Гибридная ИС является результатом интеграции традиционных гироскопических и магнитометрических инклинометров в единый измерительный модуль. Структурная избыточность, организованная путем объединения первичных датчиков, основанных на различных физических принципах работы, позволяет организовать процедуру резервирования и расширить функциональные возможности ИС. Помимо основной задачи пространственной ориентации скважины гибридная ИС позволяет контролировать целостность стенок обсадных труб при восстановлении старых месторождений, осуществлять поиск тел с повышенной магнитной восприимчивостью и определять положение магнитных масс [1].
Очевидно, что для анализа возможности проведения измерений, полноценного функционирования ИС со структурной избыточностью, последующей оценки достоверности измерений на базе гироскопических или магнитометрических датчиков необходима реализация комплексной процедуры принятия решения. Однако в опубликованных на сегодняшний день научных работах, посвященных гибридным ИС, упоминания о подобных системах отсутствуют.
Авторами настоящей работы предлагается система поддержки принятия решения для гибридной ИС, алгоритм работы которой основан на процедуре регистрации недопустимых магнитных аномалий, и критерии выбора достоверного измерения.
Процедура регистрации недопустимых магнитных аномалий
Предлагаемая концепция информационной поддержки принятия решения с привлечением технологии картографического веб-сервиса основана на методике решения задач аналитического контроля, исследования, моделирования и визуализации параметров геомагнитного поля (ГМП) и его вариаций, и оценки их влияния на информационно-измерительную технику и оборудование.
Оценка текущего значения ГМВ реализуется посредством процедуры сравнения. Для этого текущие параметры ГМП сравниваются с эталонными значениями, например с рассчитанными значениями вектора ГМП внутриземных источников в заданных пространственно-временных координатах [2, 3]. Проведенный анализ позволяет судить о степени отклонения и характере распределения полученных данных от известных эталонных значений, наличии или отсутствии ГМВ в исследуемой области, а также оценить размеры и силовые характеристики выявленных аномалий.
Разработанное с этой целью веб-приложение отличается возможностью 2D/3D-визулизации параметров ГМП и его вариаций в режиме реального времени, а также включает в себя функцию автоматизированного анализа параметров ГМВ с целью решения задачи их идентификации по природе происхождения и автоматического формирования географических карт специального назначения для поддержки принятия решения при планировании и проведении инклинометрических исследований.
Созданное веб-приложение функционирует в строгом соответствии с требованиями, предъявляемыми к приложениям данного класса. Единственное, что необходимо пользователю для его эффективного использования, – это так называемый пользовательский агент, или, как его чаще называют, веб-браузер [3].
С помощью радиокнопок пользователь выбирает режим ввода геопространственных координат: вручную или из списка из базы данных. Так, к примеру, выбрав режим ввода значений широты и долготы, пользователь нажимает на кнопку «Поиск». В результате выполняется перерисовка карты применительно к введенным координатам (в информационной области) и расчет параметров геомагнитного поля (в аналитической области).
Гибкость представленного программного обеспечения обусловлена его кросс-платформенностью: приложение одинаково эффективно работает на уровне как персонального стационарного компьютера, так и для мобильных устройств.
Решение о возможности проведения инклинометрических работ принимается по результатам сравнения допустимых значений ГМВ с текущими. Оценку допустимых ГМВ предложено реализовывать следующим образом.
В гибридной ИС блок акселерометров измеряет проекции ускорения свободного падения gi на оси, связанные с корпусом скважинного модуля, блок гироскопов – проекции угловой скорости Земли ωi, блок магнитометров – проекции напряженности магнитного поля Земли hi.
Известно, что погрешность магнитометрических датчиков зависит от наличия намагниченных масс, что в свою очередь обуславливает невозможность работы магнитометрических инклинометров в обсаженных скважинах и в условиях магнитных аномалий [5, 6].
Таким образом, говорить о достоверности измерений на базе магнитометрического модуля можно только при условии отсутствия недопустимых ГМВ.
Выходной сигнал магнитометра при условии наличия ГМВ принимает вид
(1)
где Ti – «чистый» сигнал в условиях отсутствия помех; ΔT – помехи; kT – коэффициент преобразования магнитометра; Hi – проекция напряженности ГМП на ось чувствительности i-го магнитометра; ΔH – ГМВ, обусловленные геомагнитными аномалиями, наличием обсадных колонн или магнитонеустойчивыми зонами; T0 – нулевой сигнал.
Значение невозмущенного ГМП H рассчитывается для каждой точки Земли и для определенного периода времени [4]. Отклонения наблюдаемых значений магнитного поля Земли H′ от нормального поля H являются аномалиями магнитного поля или ГМВ:
ΔH = H’ – H. (2)
При скважинных исследованиях отклонения напряженности ГМП от прогнозируемого значения могут свидетельствовать о наличии геомагнитных аномалий, т.е. о расположении рудных пород вблизи скважины. Следовательно, помимо основной задачи определения пространственной ориентации нефтяных и газовых скважин можно использовать гибридные ИС для решения геологических задач, а именно осуществлять поиск тел с повышенной магнитной восприимчивостью и определять положение магнитных масс в околоскважинном пространстве.
В результате нормирования выражения (1) получим
(3)
где – выходной сигнал с магнитометра при условии отсутствия ГМВ в безразмерном виде, где εh – относительное отклонение магнитного поля Земли от нормы, выраженное в процентах.
Подставляя в формулу для расчета азимутального угла значения выходных сигналов магнитометров при условии отсутствия и наличия ГМВ, получим абсолютную погрешность измерения азимута:
(4)
где θ – зенитный угол; φ – визирный угол.
Согласно РД 153-39.0-072-01 «Техническая инструкция по проведению геофизических исследований и работ приборами на кабеле в нефтяных и газовых скважинах», основная погрешность измерения азимута должна удовлетворять условию град. Таким образом, анализ полученных значений Δα позволяет получить допустимое значение ГМВ, при котором гибридная ИС удовлетворяет точностным требованиям.
Превышение заданного значения свидетельствует о наличии недопустимых для магнитометрического блока ГМВ и проведении измерений только на базе гироскопов. Картина распределения допустимых ГМВ по всей поверхности земного шара представлена на рис. 1.
Алгоритм принятия решения о достоверности измерений
Если согласно результатам, представленным на специализированном веб-сервисе, уровень текущих или прогнозируемых ГМВ не превышает допустимый уровень, необходимо проводить измерения на базе всех трех типов датчиков с последующим анализом и оценкой степени достоверности результатов.
В качестве условий, позволяющих судить о достоверности измерений на базе гироскопических и магнитометрических датчиков, предложено использовать отклонения экспериментальных показаний гироскопов ΔD и магнитометров ΔT от опорных значений, абсолютную погрешность широты месторождения Δφш и угла магнитного наклонения Δv, а также условие пересечения интервалов:
где αω, αh – азимут, рассчитанный для комбинаций и
соответственно. С обоснованием выбранных условий и результатами расчетов согласно представленным критериям можно ознакомиться в работе [5].
Рис. 1. Изолинии допустимых ГМВ при инклинометрии скважин
Разработанный алгоритм основан на попарном сравнении величин ΔD и ΔT, Δφш и Δv, а также проверке наличия или отсутствия пересечения интервалов. Результатом принятия решения являются следующие варианты:
α = αh; α = αω; [2].
Алгоритм принятия решения реализуется следующим образом. В наземный блок гибридной ИС поступают сигналы, измеренные скважинным модулем (Ai, Di, Ti (i = 1, 2, 3) – выходные сигналы, поступающие с триады акселерометров, гироскопов и магнитометров). Далее происходит расчет пространственных углов (αω, αh, φ, θ), а также широты месторождения и угла магнитного наклонения v′. На следующем шаге рассчитываются отклонения ΔD и ΔT, абсолютные погрешности Δφш, Δv. Следующий этап – реализация процедуры выбора достоверного значения α из полученных αω, αh.
Результаты тестирования алгоритма подтвердили эффективность включения первичных датчиков различной физической природы в измерительный модуль ИС. Так, например, на широте Δφш = 56,05° точность восстановления азимута при ГМВ εh = 1,5% увеличилась в 1,66 раза, а при εh = 3% в 3,83 раза (рис. 2), на широте Δφш = 67,45° в 2,62 и 6,56 раза соответственно (рис. 3). В среднем при уровне ГМВ εh = 3% точность восстановления азимута на средних широтах увеличивается в 5,19 раза, а на высоких в 7,68 раза.
Рис. 2. Отношение (i) абсолютных погрешностей азимута до и после реализации алгоритма принятия решения в зависимости от α (φш = 56,05°)
Рис. 3. Отношение (i) абсолютных погрешностей азимута до и после реализации алгоритма принятия решения в зависимости от α (φш = 67,45°)
Заключение и выводы
Таким образом, методика расчета параметров ГМП и его вариаций, картографическая визуализация областей геомагнитных аномалий в двух- и трехмерном виде с ранжированием по степени возможного воздействия на показания информационно-измерительной техники, сравнение допустимых значений ГМВ с текущими позволяет в режиме реального времени предоставить информацию, необходимую для поддержки принятия решения при планировании и проведении инклинометрических исследований. В свою очередь, комплексный подход, основанный на привлечении информации с картографического веб-сервиса, и применение синтезированного алгоритма принятия решения, позволяют повысить достоверность измерений, надежность и точность ИС, сократить временные и стоимостные затраты.
Работа поддержана грантами РФФИ и грантом Президента Российской Федерации: № 14-07-00260-а, 14-07-31344-мол_а, 15-17-20002-Д_С, 15-07-02731-а, МК-5340.2015.9.
Рецензенты:
Веревкин А.П., д.т.н., профессор, зав. кафедрой автоматизации технологических процессов и производств, ФГБОУ ВПО «Уфимский государственный нефтяной технический университет», г. Уфа;
Даринцев О.В., д.т.н., ведущий научный сотрудник, зав. лабораторией «Робототехника и управление в технических системах», ФГБУН «Институт механики им. Р.Р. Мавлютова», г. Уфа.