СИСТЕМА И СПОСОБ РЕКОНСТРУКЦИИ С ИСПОЛЬЗОВАНИЕМ "МЯГКОГО ПОЛЯ"
SYSTEM AND METHOD FOR RECONSTRUCTION WITH "SOFT FIELD"

Предпосылки создания изобретения
[0001] Настоящее изобретение относится в общем к системам и способам реконструкции данных и, в частности, к системам и способам оценки свойств исследуемых областей, в частности, при реконструкции объектов, содержащих множество материалов, с использованием «мягкого поля».
[0002] Томография с использованием «мягкого поля», например, электроимпедансная томография (Electrical Impedance Tomography, EIT), диффузная оптическая томография, эластография или связанные с ними технологии могут применяться для измерения внутренних свойств объекта, таких как электрические свойства материалов, включая внутренние структуры объекта. Например, в EIT-системах оценивают распределение электрической проводимости внутренних структур объекта. В упомянутых EIT-системах реконструируют проводимость и/или диэлектрическую проницаемость материалов в области или в объеме на основе прикладываемого возбуждения (например, тока) и измеренного отклика (например, напряжения), получаемого на поверхности упомянутой области или объема или вблизи этой поверхности. Впоследствии могут формироваться визуальные распределения этих оценок.
[0003] При томографии с использованием «мягкого поля» традиционные алгоритмы реконструкции могут получать распределение импеданса в объекте без использования априорной информации. Однако подобные процедуры реконструкции требуют больших вычислительных мощностей вследствие итераций, необходимых для сходимости решения. Таким образом, процедура реконструкции в традиционных алгоритмах может требовать много времени и высокоскоростной электроники или процессоров. В соответствии с этим, если необходимы быстрые измерения, например, для визуализации объекта, содержащего множество материалов, в реальном времени, к примеру, потока газа через трубу, традиционные алгоритмы не могут работать удовлетворительно. Кроме того, традиционные алгоритмы реконструкции не способны работать с высоким контрастом в реальном распределении импеданса объекта.
Сущность изобретения
[0004] В соответствии с одним из вариантов осуществления настоящего изобретения предлагается способ реконструкции для томографии с использованием «мягкого поля». Способ включает получение информации о приложенных входных сигналах и измеренных выходных сигналах для возбуждаемого объекта с использованием множества преобразователей и формирование матрицы полной проводимости на основе упомянутой информации о приложенных входных сигналах и измеренных выходных сигналах. Способ включает также определение множества моментов с использованием упомянутой матрицы полной проводимости и вычисление распределения свойств возбуждаемого объекта с использованием упомянутого множества моментов.
[0005] В соответствии с другим вариантом осуществления настоящего изобретения предлагается способ реконструкции для томографии с использованием «мягкого поля». Способ включает получение информации о приложенных входных сигналах и измеренных выходных сигналах для возбуждаемого объекта с использованием множества преобразователей и выполнение итеративной реконструкции путем преобразования симметричных компонентов (symmetrical component transform, SCT) с использованием упомянутой информации о приложенных входных сигналах и измеренных выходных сигналах. Способ включает также определение распределения свойств возбуждаемого объекта на основе итеративной реконструкции SCT.
[0006] В соответствии с еще одним вариантом осуществления настоящего изобретения предлагается система томографии с использованием «мягкого поля», которая включает множество преобразователей, сконфигурированных для их размещения вблизи поверхности объекта, и одно или более устройств подачи возбуждения, связанных с упомянутым множеством преобразователей и сконфигурированных для формирования сигналов возбуждения для упомянутого множества преобразователей. Упомянутая система томографии с использованием «мягкого поля» включает также один или более детекторов отклика, связанных с упомянутым множеством преобразователей и сконфигурированных для измерения в упомянутом множестве преобразователей отклика объекта на возбуждение, приложенное посредством упомянутого множества преобразователей, на основе упомянутых сигналов возбуждения. Система томографии с использованием «мягкого поля» включает также модуль реконструкции с использованием «мягкого поля», сконфигурированный для реконструкции распределения свойств объекта на основе упомянутых сигналов возбуждения и измеренного отклика с применением множества моментов, определенных из процедуры реконструкции матрицы полной проводимости.
Краткое описание чертежей
[0007] Настоящее изобретение станет более понятным из последующего описания вариантов его осуществления со ссылками на приложенные чертежи.
[0008] На фиг.1 представлена упрощенная структурная схема системы томографии с использованием «мягкого поля» в соответствии с различными вариантами осуществления настоящего изобретения.
[0009] На фиг.2 представлен вид в перспективе конфигурации преобразователя в соответствии с одним из вариантов осуществления настоящего изобретения.
[0010] На фиг.3 представлена упрощенная схема, иллюстрирующая реконструкцию упомянутого распределения свойств.
[0011] На фиг.4 представлена схема, иллюстрирующая информационный поток для томографии с использованием «мягкого поля» в соответствии с различными вариантами осуществления настоящего изобретения.
[0012] На фиг.5 представлена схема определения полной проводимости в соответствии с различными вариантами осуществления настоящего изобретения.
[0013] На фиг.6 представлена блок-схема способа реконструкции с использованием «мягкого поля» с применением симметричных компонентов в соответствии с различными вариантами осуществления настоящего изобретения.
[0014] На фиг.7 представлена таблица, в которой показаны значения моментов инерции, вычисленных в соответствии с различными вариантами осуществления настоящего изобретения.
[0015] На фиг.8 представлена блок-схема способа реконструкции с использованием «мягкого поля» и симметричных компонентов в соответствии с другими вариантами осуществления настоящего изобретения.
[0016] На фиг.9 представлена диаграмма координатной сетки, используемой в процедуре преобразования полярной координатной сетки, в соответствии с различными вариантами осуществления настоящего изобретения.
[0017] На фиг.10 переставлена диаграмма, иллюстрирующая процедуру преобразования полярной координатной сетки в соответствии с различными вариантами осуществления настоящего изобретения.
[0018] На фиг.11 представлена блок-схема способа выполнения итеративной реконструкции с использованием преобразования симметричных компонентов (SCT) в соответствии с различными вариантами осуществления настоящего изобретения.
Подробное описание изобретения
[0019] Более полное понимание предшествующего краткого описания, а также последующего подробного описания конкретных вариантов осуществления настоящего изобретения может быть обеспечено при их прочтении в комбинации с приложенными чертежами. В случаях, когда чертежи иллюстрируют схемы функциональных блоков различных вариантов осуществления изобретения, эти функциональные блоки не обязательно указывают на разделение между схемами аппаратного обеспечения. Соответственно, например, один или более функциональных блоков (например, процессоров, контроллеров, схем или блоков памяти) могут быть реализованы в одном элементе или в нескольких элементах аппаратного обеспечения. Следует понимать, что различные варианты осуществления настоящего изобретения не ограничены схемами, взаимосвязями компонентов/элементов или средствами, показанными на чертежах.
[0020] В настоящем описании модуль или шаг, упомянутый в единственном числе, не исключает множества упомянутых элементов или шагов, если только такое исключение явно не указано. Кроме того, выражение "один из вариантов осуществления изобретения" не исключает существование других вариантов осуществления настоящего изобретения, также включающих приведенные признаки. Кроме того, если иное явно не указано, варианты осуществления изобретения, "включающие" или "имеющие" модуль или множество модулей с определенным свойством, могут включать дополнительные подобные модули, не обладающие упомянутым свойством.
[0021] В различных вариантах осуществления настоящего изобретения предлагается система и способ томографии с использованием «мягкого поля», в частности объектов, включающих множество материалов, в которых применяют симметричные компоненты для оценки свойств объекта, например свойств текучих газов. В различных вариантах осуществления настоящего изобретения предлагается подход, в котором осуществляют итерацию по определяемым моментам (вместо итерации по измеренным токам). По меньшей мере одним техническим результатом различных вариантов осуществления настоящего изобретения является повышение точности и скорости визуализации свойств объектов, включающих множество материалов. Например, в промышленных приложениях путем практического использования по меньшей мере одного из вариантов осуществления настоящего изобретения может быть быстро получено реконструированное визуальное представление распределения потока газа в трубопроводе, например, в реальном времени при протекании газа.
[0022] Следует отметить, что в настоящем описании термин "томография с использованием «мягкого поля»" относится в общем к любому томографическому, или многомерному, расширению способа томографии, который не является томографией с использованием «жесткого поля».
[0023] Один из вариантов осуществления системы 20 томографии с использованием «мягкого поля» проиллюстрирован на фиг.1. Например, система 20 томографии с использованием «мягкого поля» может представлять собой систему электроимпедансной томографии (Е1Т), применяемую для определения электрических свойств материалов в объекте 22, в частности объекте, включающем множество материалов (как показано на фиг.2). В одном из вариантов осуществления настоящего изобретения внутри объекта 22 или в другой области или объеме может определяться пространственное распределение электрической проводимости (σ) и/или диэлектрической проницаемости (ε). Например, система 20 томографии с использованием «мягкого поля» обеспечивает EIT для многофазных измерений потока внутри объекта 22, таких как визуализация свойств или объемной скорости потока газа или нефти в потоке нефтепродуктов в трубе.
[0024] В проиллюстрированном варианте осуществления настоящего изобретения система 20 включает множество преобразователей 24 (например, электродов), расположенных в объекте, например, разнесенных друг от друга по внутренней окружности трубы 42 и находящихся в контакте с текучей средой, как показано на фиг.2. В одном из вариантов осуществления настоящего изобретения множество колец 40 преобразователей 24 расположены по внутренней длине, например, разнесены на расстояние D (например, один метр) друг от друга вдоль оси трубы 42. Например, преобразователи 24 (к примеру, электроды, тепловые источники, ультразвуковые преобразователи) могут быть расположены на поверхности внутренней окружности трубы 42, вблизи поверхности или проходить через поверхность снаружи трубы 42 внутрь нее (например, игольчатые электроды). Соответственно, преобразователи 24 могут иметь различную форму, например, представлять собой, помимо прочего, электроды поверхностного контакта, изолированные электроды, электроды с емкостной связью, а также проводящие витки, например, антенны. Разнесенные друг от друга кольца 40 могут формировать многофазный потоковый измеритель в соответствии с различными вариантами осуществления настоящего изобретения для определения, например, количества газа в трубе 42 (на основе визуализации газа и нефти в трубе 42) и скорости потока на основе разности измерений между кольцами 40 (в двух местах в трубе 42), например, путем выполнения взаимной корреляции. Таким образом, в соответствии с различными вариантами осуществления настоящего изобретения реконструкцию осуществляют в каждом из колец 40. Соответственно, могут быть обеспечены объемная визуализация распределения свойств газа или определение потока газа в трубе 42, например, для определения количества газа и нефти, протекающих через нее.
[0025] Следует отметить, что система 20 томографии с использованием «мягкого поля» может представлять собой систему другого типа. Например, система 20 томографии с использованием «мягкого поля» может представлять собой, помимо прочих, систему диффузной оптической томографии (Diffuse Optical Tomography, DOT), систему спектроскопии ближнего ИК-диапазона (Near InfraRed Spectroscopy, NIRS), систему термографии, систему эластографии или систему СВЧ-томографии.
[0026] Устройство 26 подачи возбуждения и детектор 28 отклика подключены к приемопередатчикам 24, каждый из которых соединен с модулем 30 реконструкции с использованием «мягкого поля». Модуль 30 реконструкции с использованием «мягкого поля» может представлять собой любой тип процессора или вычислительного устройства, выполняющего реконструкцию с использованием «мягкого поля» по меньшей мере частично на основе принятых откликов от преобразователей 24, как более подробно будет описано ниже. Например, модуль 30 реконструкции с использованием «мягкого поля» может представлять собой аппаратное обеспечение, программное обеспечение или их комбинацию. В одном из вариантов осуществления настоящего изобретения устройство 26 подачи возбуждения и детектор 28 отклика являются физически отдельными устройствами. В других вариантах осуществления настоящего изобретения устройство 26 подачи возбуждения и детектор 28 отклика физически интегрированы в одном элементе. Имеется также контроллер 33, который передает инструкции в устройство 26 подачи возбуждения, управляющее преобразователями 24 на основе этих инструкций. Следует отметить, что устройство 26 подачи возбуждения может быть соединено со всеми или некоторыми преобразователями 24.
[0027] Следует также отметить, что могут использоваться различные типы возбуждений для получения данных распределения свойств для их использования в процедуре реконструкции в соответствии с различными вариантами осуществления настоящего изобретения. Например, помимо прочих, в сочетании с различными вариантами осуществления настоящего изобретения могут быть использованы электрические, магнитные, оптические, тепловые или ультразвуковые возбуждения. В упомянутых различных вариантах осуществления настоящего изобретения преобразователи 24 могут быть связаны с объектом 22 различными способами и не обязательно должны находиться в непосредственном контакте или располагаться только на внутренней поверхности объекта 22 (например, они могут иметь электрическую, емкостную, гальваническую связь и т.п.).
[0028] В некоторых вариантах осуществления настоящего изобретения система 20 томографии с использованием «мягкого поля» может использоваться для формирования визуального представления распределения электрического импеданса в множестве различных приложений, например, для определения свойств материала в потоке смешанной текучей среды, включающей нефть и воду (или другие жидкости или газы, например, нефтепродукты), или для подземных областей для анализа почвы или проверки полотна дороги. Однако упомянутые варианты осуществления настоящего изобретения могут быть использованы и в других приложениях, например, когда объект 22 представляет собой часть тела человека, такую как голова, грудь или нога, в которых воздух и ткани имеют различные электрические проводимости.
[0029] В различных вариантах осуществления настоящего изобретения упомянутые преобразователи выполнены из любого подходящего материала. Например, тип применяемых преобразователей 24 может зависеть от конкретного применения, чтобы соответствующий тип приемопередатчика (например, электрод, катушка и т.п.) использовался для формирования возбуждений "мягкого поля" (например, электромагнитного поля) и для приема откликов объекта 22 на возбуждения для конкретного применения. В некоторых вариантах осуществления настоящего изобретения для формирования электрического тока может применяться проводящий материал. Например, преобразователи 24 могут быть выполнены из одного или более металлов, таких как медь, золото, платина, сталь, серебро и их сплавы. Другие примеры материалов для формирования преобразователей 24 включают неметаллические электропроводные материалы, например, материалы на основе кремния, использованные в сочетании с микросхемами. В одном из вариантов осуществления настоящего изобретения, в котором объект 22 представляет собой трубу 42 с потоком газа или текучей среды, протекающими через нее, преобразователи 24 выполнены влагоустойчивыми. Дополнительно, преобразователи 24 могут иметь различную форму и/или размеры, например, представлять собой структуры в форме стержня, плоской пластины или структуры игольчатой формы. Следует отметить, что в некоторых вариантах осуществления настоящего изобретения преобразователи 24 изолированы друг от друга. В других вариантах осуществления настоящего изобретения преобразователи 24 могут быть размещены в непосредственном омическим контакте с объектом 22 или могут иметь емкостную связь с объектом 22.
[0030] При функционировании системы преобразователи 24 или их подмножество могут быть использованы для передачи сигналов (например, доставки или модуляции сигналов), к примеру, для непрерывной доставки электрического тока или для доставки изменяющегося во времени сигнала таким образом, что возбуждение может быть приложено в заданном частотно-временном диапазоне (например, 1 кГц - 1 МГц) для формирования электромагнитного поля в объекте 22. В EIT-применениях измеряют результирующие поверхностные потенциалы, а именно напряжения на приемопередатчиках 24, для определения распределения электрической проводимости или диэлектрической проницаемости с использованием способов реконструкции согласно настоящему описанию. Например, может реконструироваться визуальное распределение на основе геометрии преобразователей 24, приложенных токов и измеренных напряжений.
[0031] Таким образом, в различных вариантах осуществления настоящего изобретения устройство 26 подачи возбуждения прикладывает возбуждение к каждому из преобразователей 24, а детектор 28 отклика измеряет отклик на каждом из преобразователей 24 (которые могут мультиплексироваться мультиплексором) в ответ на приложенное к приемопередатчикам 24 возбуждение. Следует отметить, что может использоваться возбуждение любого типа, например, помимо прочего, электрический ток, электрическое напряжение, магнитное поле, радиочастотная волна, тепловое поле, оптический сигнал, механическая деформация или ультразвуковой сигнал.
[0032] К примеру, в EIT-применениях, как показано на фиг.3, осуществляют реконструкцию для определения исследуемых областей 32 в объекте 22. Как показано на фиг.3, детектор 28 отклика (показанный на фиг.1) измеряет напряжение отклика (или ток отклика) на преобразователях 24 в ответ на ток или напряжение, прикладываемые устройством 26 подачи возбуждения (показанным на фиг.1) к преобразователям 24.
[0033] Следует отметить, что детектор 28 отклика может включать один или более элементов преобразования аналогового сигнала (не показаны на чертеже), которые усиливают и/или фильтруют измеренное напряжение или ток отклика. В других вариантах осуществления настоящего изобретения процессор системы 20 томографии с использованием «мягкого поля» включает элемент преобразования сигнала для усиления и/или фильтрации напряжения отклика или тока отклика, принятого от детектора 28 отклика.
[0034] Модуль 30 реконструкции с использованием «мягкого поля», соответственно, вычисляет отклик объекта 22 на приложенное возбуждение. В соответствии с этим, система 20 томографии с использованием «мягкого поля» может быть использована для реконструкции распределения свойств или для визуализации потока.
[0035] В системе 20 томографии с использованием «мягкого поля», в различных вариантах осуществления настоящего изобретения, используется подход с использованием симметричных компонентов (называемых также симметрическими). Например, информационный поток 48 EIT, в котором применяют подход с использованием симметричных компонентов и определяемых моментов для реконструкции, проиллюстрирован на фиг.4. В частности, карту 50 полной проводимости, сформированную из одной или более матриц (например, заранее вычисленных матриц), основанных на возбуждениях от вычислительного устройства 52, используют для предсказания напряжений (предсказанные данные), подаваемых в модуль 30 реконструкции с использованием «мягкого поля». Следует отметить, что в некоторых вариантах осуществления настоящего изобретения карту полной проводимости используют опционально. Возбуждение прикладывают к объекту 22 (показан на фиг.1-3) с помощью инструмента 54 томографии с использованием «мягкого поля», который может включать преобразователи 24 и другие компоненты для возбуждения и измерений, при этом измеренные напряжения (измеренные данные) также подают в модуль 30 реконструкции с использованием «мягкого поля». Модуль 30 реконструкции с использованием «мягкого поля» выполняет затем реконструкцию с использованием различных вариантов осуществления настоящего изобретения для формирования оценки распределения 56 свойств, например, распределения импеданса, для идентификации исследуемых областей 32 внутри объекта 22, например, содержимого различных материалов в потоке жидкости или газа. Следует отметить, что упомянутые различные компоненты могут быть физически отдельными компонентами или элементами или могут быть объединены. Например, модуль 30 реконструкции с использованием «мягкого поля» может составлять часть системы 20 томографии с использованием «мягкого поля» (как показано на фиг.1).
[0036] Посредством применения различных вариантов осуществления настоящего изобретения обеспечивают реконструкцию с использованием «мягкого поля», которая позволяет определить распределение импеданса без использования априорной информации (а также без использования модели предсказания). Например, в одном из вариантов осуществления настоящего изобретения может выполняться определение импеданса или полной проводимости в соответствии с алгоритмом 60 определения полной проводимости, показанным на фиг.5. Определение полной проводимости включает применение прикладываемых входных сигналов (например, возбуждений) на шаге 62 и измеряемых выходных сигналов (например, откликов) с помощью инструмента 54 томографии с использованием «мягкого поля» в качестве входных данных для построения матрицы 64 полной проводимости (или матрицы импеданса), которая задает карту 50 полной проводимости или карту импеданса, соответственно (показаны на фиг.4), что будет более подробно описано далее. Матрица 64 полной проводимости включает элементы полной проводимости, связанные с моментами (например, электромагнитными моментами) следующим образом:
<mrow><msub><mi>Y</mi><mi>s</mi></msub><mo>=</mo><mi>f</mi><mrow><mo>(</mo><mi>M</mi><mo>)</mo></mrow><mtext> </mtext><mi>В</mi><mi>ы</mi><mi>р</mi><mi>а</mi><mi>ж</mi><mi>е</mi><mi>н</mi><mi>и</mi><mi>е</mi><mtext> </mtext><mn>1</mn></mrow>
<mrow><msubsup><mi>Y</mi><mi>s</mi><mrow><mi>p</mi><mo>,</mo><mi>q</mi></mrow></msubsup><mo>=</mo><msub><mi>f</mi><mrow><mi>p</mi><mo>,</mo><mi>q</mi></mrow></msub><mrow><mo>(</mo><mrow><msub><mi>M</mi><mrow><mi>p</mi><mo>,</mo><mi>q</mi></mrow></msub></mrow><mo>)</mo></mrow><mo>=</mo><msub><mi>f</mi><mrow><mi>p</mi><mo>,</mo><mi>q</mi></mrow></msub><mstyle displaystyle="true"><mrow><mo>∬</mo><mrow><mi>σ</mi><mrow><mo>(</mo><mrow><mi>x</mi><mo>,</mo><mi>y</mi></mrow><mo>)</mo></mrow><msup><mrow><mrow><mo>(</mo><mrow><mi>x</mi><mo>+</mo><mi>i</mi><mi>y</mi></mrow><mo>)</mo></mrow></mrow><mrow><mi>p</mi><mo>−</mo><mn>2</mn></mrow></msup><msup><mrow><mrow><mo>(</mo><mrow><mi>x</mi><mo>−</mo><mi>i</mi><mi>y</mi></mrow><mo>)</mo></mrow></mrow><mrow><mi>q</mi><mo>−</mo><mn>2</mn></mrow></msup><mi>d</mi><mi>x</mi><mi>d</mi><mi>y</mi><mo>,</mo><mtext> </mtext><mi>В</mi><mi>ы</mi><mi>р</mi><mi>а</mi><mi>ж</mi><mi>е</mi><mi>н</mi><mi>и</mi><mi>е</mi><mtext> </mtext><mn>2</mn></mrow></mrow></mstyle></mrow>
s
*
где b - коэффициент, связывающий Y с упомянутыми моментами, который может быть вычислен заранее на основе моделирования и т.п. Следует отметить, что в различных вариантах осуществления настоящего изобретения σ=σ. Таким образом, в настоящем описании термин "проводимость" в различных вариантах осуществления настоящего изобретения включает следующие три электрических свойства: проводимость (σ), магнитную проницаемость (µ) и диэлектрическую проницаемость (ε). Соответственно, в различных вариантах осуществления настоящего изобретения различные выражения в настоящем описании включают σ, µ и ε.
[0037] Матрица моментов может быть задана следующим образом:
<mrow><mi>M</mi><mo>=</mo><mi>S</mi><mi>σ</mi><mo>⇒</mo><mi>σ</mi><mo>=</mo><msup><mrow><mrow><mo>(</mo><mrow><msup><mi>S</mi><mi>H</mi></msup><mi>S</mi></mrow><mo>)</mo></mrow></mrow><mrow><mo>−</mo><mn>1</mn></mrow></msup><msup><mi>S</mi><mi>H</mi></msup><mi>M</mi><mo>,</mo><mtext> </mtext><mi>В</mi><mi>ы</mi><mi>р</mi><mi>а</mi><mi>ж</mi><mi>е</mi><mi>н</mi><mi>и</mi><mi>е</mi><mtext> </mtext><mn>3</mn></mrow>
<mrow><mi>S</mi><mo>=</mo><mstyle displaystyle="true"><mrow><mo>∬</mo><mrow><msup><mrow><mrow><mo>(</mo><mrow><mi>x</mi><mo>+</mo><mi>i</mi><mi>y</mi></mrow><mo>)</mo></mrow></mrow><mrow><mi>p</mi><mo>−</mo><mn>2</mn></mrow></msup><msup><mrow><mrow><mo>(</mo><mrow><mi>x</mi><mo>−</mo><mi>i</mi><mi>y</mi></mrow><mo>)</mo></mrow></mrow><mrow><mi>q</mi><mo>−</mo><mn>2</mn></mrow></msup><mi>d</mi><mi>x</mi><mi>d</mi><mi>y</mi><mo>,</mo><mtext> </mtext></mrow></mrow></mstyle><mi>В</mi><mi>ы</mi><mi>р</mi><mi>а</mi><mi>ж</mi><mi>е</mi><mi>н</mi><mi>и</mi><mi>е</mi><mtext> </mtext><mn>4</mn></mrow>
H
-1
H
где S - структурная матрица, которая может быть вычислена заранее, например, исходя из геометрии инструмента 54 томографии с использованием «мягкого поля» (показан на фиг.4), σ - определяемая проводимость, (SS)S - псевдообращение структурной матрицы (например, многострочное обращение), М - матрица моментов. Таким образом, в одном из вариантов осуществления настоящего изобретения, экспериментально измеренные моменты связывают с использованием приведенных выше выражений для оценки распределения проводимостей.
[0038] Следует отметить, что в настоящем описании символ σ обозначает как проводимость в точке, так и вектор проводимости.
H
-1
H
-1
-1
s
transformed
s
[0039] Соответственно, значения полной проводимости (Y) могут быть использованы для определения информации о распределении, включая физические и геометрические свойства распределения, в соответствии с шагами 66 и 68, где в общем Y=I.(VV).V и Y=Y=P.Y.Q, соответственно. В одном из вариантов осуществления настоящего изобретения, Р=Q=матрица дискретного преобразования Фурье (discrete Fourier transform, DFT), С. Соответственно, в этом варианте осуществления изобретения Y=C.Y.C. Следует, однако, отметить, что может использоваться любой способ или алгоритм, основанный на матричном преобразовании, и С является лишь одним из примеров, использованным для описания одного из вариантов осуществления настоящего изобретения.
[0040] Матрицу 64 полной проводимости преобразуют на шаге 70 в преобразованную матрицу полной проводимости, как будет более подробно описано далее, с использованием подхода, заключающегося в преобразовании симметричных компонентов (SCT), в котором Y задает распределение в виде дискретной матрицы. Следует отметить, что помимо описанных в настоящем документе, могут использоваться и другие способы преобразования и обращения. Соответственно подход SCT используют для случая, где Р=Q=матрица DFT, С.
[0041] В одном из вариантов осуществления настоящего изобретения предлагается способ 80 реконструкции с использованием «мягкого поля» и симметричных компонентов, показанный на фиг.6. Вначале будет дано общее описание способа 80, затем конкретное описание различных шагов, включая различные реализации. Способ 80 включает приложение сигналов возбуждения на шаге 82 и измерение откликов на шаге 84. Например, к объекту могут прикладываться электрические токи с использованием множества преобразователей в соответствии с настоящим описанием, при этом отклики измеряют в каждом из упомянутых преобразователей. Следует отметить, что сигналы возбуждения, прикладываемые на шаге 82, могут представлять собой, например, любые ортонормированные сигналы.
[0042] После этого на шаге 86 строят матрицу полной проводимости с использованием приложенных сигналов возбуждения и измеренных откликов. Например, упомянутая матрица полной проводимости может задавать значения для приложенных и измеренных сигналов каждого из упомянутого множества преобразователей. Матрицу полной проводимости затем преобразуют на шаге 88 с использованием процедур умножения слева и умножения справа, например, обращают с использованием матрицы преобразования Фурье.
[0043] Следует отметить, что упомянутое распределение внутри объекта предполагается равномерным, при этом формируют сетку, например, из квадратных, прямоугольных или округлых элементов. Соответственно, упомянутая координатная сетка может быть задана аналитически. Следует, однако, отметить, что вместо подхода на основе координатной сетки может использоваться подход на основе формы, как будет более подробно описано далее.
[0044] Затем вычисленные моменты, заданные преобразованной матрицей, используют на шаге 90 в итеративной процедуре вычисления для оценки, например, электрической проводимости (σ) внутри объекта для реконструкции распределения свойств внутри объекта, например, внутри объекта, включающего множество материалов. Например, может определяться разность вычисленных моментов и обновляться оцененный параметр σ до тех пор, пока не будет достигнута сходимость (в пределах заранее заданного уровня).
[0045] Используемая сетка может быть уточнена на шаге 92. Например, могут быть выбраны области для уточнения сетки, при этом уточненную сетку вычисляют аналитически на основе ранее определенных ограниченных аномалий или различных материалов. Затем на шаге 94 выполняется итеративная процедура вычисления на уточненной сетке.
s
[0046] В отношении конкретного построения матрицы полной проводимости на шаге 86, в одном из вариантов осуществления настоящего изобретения каждый элемент Y этой матрицы соответствует полярному моменту проводимости, задаваемому следующим набором выражений:
<mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow><msub><mi>Y</mi><mi>s</mi></msub><mrow><mo>(</mo><mrow><mi>p</mi><mo>,</mo><mi>q</mi></mrow><mo>)</mo></mrow><mo>=</mo><mstyle displaystyle="true"><mrow><munder><mo>∬</mo><mi>Ω</mi></munder><mrow><mi>σ</mi><mrow><mo>(</mo><mrow><mi>x</mi><mo>,</mo><mi>y</mi></mrow><mo>)</mo></mrow><msup><mrow><mrow><mo>(</mo><mrow><mi>x</mi><mo>+</mo><mi>i</mi><mi>y</mi></mrow><mo>)</mo></mrow></mrow><mrow><mi>p</mi><mo>−</mo><mn>2</mn></mrow></msup><msup><mrow><mrow><mo>(</mo><mrow><mi>x</mi><mo>−</mo><mi>i</mi><mi>y</mi></mrow><mo>)</mo></mrow></mrow><mrow><mi>q</mi><mo>−</mo><mn>2</mn></mrow></msup><mi>d</mi><mi>x</mi><mi>d</mi><mi>y</mi></mrow></mrow></mstyle></mrow></mtd></mtr><mtr columnalign="left"><mtd columnalign="left"><mrow><msub><mi>Y</mi><mi>s</mi></msub><mrow><mo>(</mo><mrow><mi>p</mi><mo>,</mo><mi>q</mi></mrow><mo>)</mo></mrow><mo>=</mo><msub><mi>σ</mi><mn>0</mn></msub><mo>.</mo><mtext> </mtext><mstyle displaystyle="true"><mrow><munder><mo>∬</mo><mrow><msub><mi>Ω</mi><mrow><mi>hom</mi><mi>o</mi><mi>g</mi><mi>e</mi><mi>n</mi><mi>e</mi><mi>o</mi><mi>u</mi><mi>s</mi></mrow></msub></mrow></munder><mrow><msup><mrow><mrow><mo>(</mo><mrow><mi>x</mi><mo>+</mo><mi>i</mi><mi>y</mi></mrow><mo>)</mo></mrow></mrow><mrow><mi>p</mi><mo>−</mo><mn>2</mn></mrow></msup><msup><mrow><mrow><mo>(</mo><mrow><mi>x</mi><mo>−</mo><mi>i</mi><mi>y</mi></mrow><mo>)</mo></mrow></mrow><mrow><mi>q</mi><mo>−</mo><mn>2</mn></mrow></msup><mi>d</mi><mi>x</mi><mi>d</mi><mi>y</mi></mrow></mrow></mstyle><mtext> </mtext><mo>∵</mo><msub><mi>σ</mi><mn>0</mn></msub><msub><mi>σ</mi><mrow><mi>a</mi><mi>n</mi><mi>o</mi><mi>m</mi><mi>a</mi><mi>l</mi><mi>y</mi></mrow></msub></mrow></mtd></mtr><mtr columnalign="left"><mtd columnalign="left"><mrow><mtable columnalign="right"><mtr columnalign="right"><mtd columnalign="right"><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mtable columnalign="left"><mtr><mtd><msub><mi>Y</mi><mi>s</mi></msub><mrow><mo>(</mo><mrow><mi>p</mi><mo>,</mo><mi>q</mi></mrow><mo>)</mo></mrow><mo>=</mo><msub><mi>σ</mi><mn>0</mn></msub><mo>.</mo><mstyle displaystyle="true"><mrow><munder><mo>∬</mo><mi>Ω</mi></munder><mrow><msup><mrow><mrow><mo>(</mo><mrow><mi>x</mi><mo>+</mo><mi>i</mi><mi>y</mi></mrow><mo>)</mo></mrow></mrow><mrow><mi>p</mi><mo>−</mo><mn>2</mn></mrow></msup><msup><mrow><mrow><mo>(</mo><mrow><mi>x</mi><mo>−</mo><mi>i</mi><mi>y</mi></mrow><mo>)</mo></mrow></mrow><mrow><mi>q</mi><mo>−</mo><mn>2</mn></mrow></msup><mi>d</mi><mi>x</mi><mi>d</mi><mi>y</mi><mo>−</mo><msub><mi>σ</mi><mn>0</mn></msub></mrow></mrow></mstyle><mo>.</mo><mtext> </mtext><mstyle displaystyle="true"><mrow><munder><mo>∬</mo><mrow><msub><mi>Ω</mi><mrow><mi>a</mi><mi>n</mi><mi>o</mi><mi>m</mi><mi>a</mi><mi>l</mi><mi>y</mi></mrow></msub></mrow></munder><mrow><msup><mrow><mrow><mo>(</mo><mrow><mi>x</mi><mo>−</mo><mi>i</mi><mi>y</mi></mrow><mo>)</mo></mrow></mrow><mrow><mi>p</mi><mo>−</mo><mn>2</mn></mrow></msup><msup><mrow><mrow><mo>(</mo><mrow><mi>x</mi><mo>−</mo><mi>i</mi><mi>y</mi></mrow><mo>)</mo></mrow></mrow><mrow><mi>q</mi><mo>−</mo><mn>2</mn></mrow></msup><mi>d</mi><mi>x</mi><mi>d</mi><mi>y</mi></mrow></mrow></mstyle></mtd></mtr><mtr><mtd><mtable columnalign="right"><mtr columnalign="right"><mtd columnalign="right"><mrow><mi>В</mi><mi>ы</mi><mi>р</mi><mi>а</mi><mi>ж</mi><mi>е</mi><mi>н</mi><mi>и</mi><mi>е</mi><mtext> </mtext><mn>5,</mn><mtext> 6</mtext><mtext> </mtext><mi>и</mi><mtext> </mtext><mn>7</mn></mrow></mtd></mtr></mtable></mtd></mtr></mtable></mtd></mtr></mtable></mrow></mtd></mtr></mtable></mrow></mtd></mtr></mtable></mrow>
s
[0047] В приведенных выше выражениях p и q определяют строку и столбец матрицы, а х и у задают координаты пикселей. Соответственно, элементы Y задают моменты, которые используются в итеративной процедуре нахождения решения, описанной здесь. Также следует отметить, что предполагается равномерное распределение.
[0048] Следует также отметить, что следующая часть выражения 7 вычисляется заранее:
<mrow><msub><mi>σ</mi><mn>0</mn></msub><mo>.</mo><mstyle displaystyle="true"><mrow><munder><mo>∬</mo><mi>Ω</mi></munder><mrow><msup><mrow><mrow><mo>(</mo><mrow><mi>x</mi><mo>+</mo><mi>i</mi><mi>y</mi></mrow><mo>)</mo></mrow></mrow><mrow><mi>p</mi><mo>−</mo><mn>2</mn></mrow></msup><msup><mrow><mrow><mo>(</mo><mrow><mi>x</mi><mo>−</mo><mi>i</mi><mi>y</mi></mrow><mo>)</mo></mrow></mrow><mrow><mi>q</mi><mo>−</mo><mn>2</mn></mrow></msup><mi>d</mi><mi>x</mi><mi>d</mi><mi>y</mi></mrow></mrow></mstyle></mrow>
s
[0049] Элементы Y могут быть также выражены в полярных координатах следующим образом:
<mtable columnalign="left"><mtr><mtd><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow><mo>−</mo><msub><mi>σ</mi><mn>0</mn></msub><mo>.</mo><mstyle displaystyle="true"><mrow><munder><mo>∬</mo><mrow><msub><mi>Ω</mi><mrow><mi>a</mi><mi>n</mi><mi>o</mi><mi>m</mi><mi>a</mi><mi>l</mi><mi>y</mi></mrow></msub></mrow></munder><mrow><msup><mrow><mrow><mo>(</mo><mrow><mi>x</mi><mo>+</mo><mi>i</mi><mi>y</mi></mrow><mo>)</mo></mrow></mrow><mrow><mi>p</mi><mo>−</mo><mn>2</mn></mrow></msup><msup><mrow><mrow><mo>(</mo><mrow><mi>x</mi><mo>−</mo><mi>i</mi><mi>y</mi></mrow><mo>)</mo></mrow></mrow><mrow><mi>q</mi><mo>−</mo><mn>2</mn></mrow></msup><mi>d</mi><mi>x</mi><mi>d</mi><mi>y</mi></mrow></mrow></mstyle></mrow></mtd></mtr><mtr columnalign="left"><mtd columnalign="left"><mrow><mo>−</mo><msub><mi>σ</mi><mn>0</mn></msub><mo>.</mo><mstyle displaystyle="true"><mrow><munder><mo>∬</mo><mrow><msub><mi>Ω</mi><mrow><mi>a</mi><mi>n</mi><mi>o</mi><mi>m</mi><mi>a</mi><mi>l</mi><mi>y</mi></mrow></msub></mrow></munder><mrow><msup><mi>r</mi><mrow><mi>p</mi><mo>+</mo><mi>q</mi><mo>−</mo><mn>4</mn></mrow></msup><msup><mi>e</mi><mrow><mrow><mo>(</mo><mrow><mi>p</mi><mo>−</mo><mi>q</mi></mrow><mo>)</mo></mrow><mi>i</mi><mi>θ</mi></mrow></msup><mi>d</mi><mi>x</mi><mi>d</mi><mi>y</mi></mrow></mrow></mstyle></mrow></mtd></mtr></mtable><mo>=</mo></mtd></mtr><mtr><mtd><mi>В</mi><mi>ы</mi><mi>р</mi><mi>а</mi><mi>ж</mi><mi>е</mi><mi>н</mi><mi>и</mi><mi>е</mi><mtext> </mtext><mn>8</mn></mtd></mtr></mtable>
s
[0050] В различных вариантах осуществления настоящего изобретения свойства элемента Y могут быть организованы в таблицу 100, как показано на фиг.7, в которой представлены значения вычисленных моментов инерции. В таблице 100 строки представляют значения р, а столбцы - значения q для выражений 5-8, приведенных выше. Следует отметить, что значения в части 102 полезны, когда импеданс комплексный, то есть одновременно резистивный и реактивный. Дополнительно, значения в части 102 являются по существу избыточными вследствие периодичности коэффициентов, используемых в указанных выражениях.
H
-1
H
-1
-1
s
s
[0051] В различных вариантах осуществления настоящего изобретения предлагается способ 110, показанный на фиг.8 и определяющий подход для томографии EIT с использованием симметричных компонентов. В частности, способ 110 включает вычисление на шаге 112 значения Y полной проводимости на основе приложенных возбуждений и измеренных откликов, например, на основе приложенных токов (I) и измеренных напряжений (V) с использованием множества преобразователей 24 (показаны на фиг.1). Например, Y может быть вычислено как Y=I.(VV).V. После этого, на шаге 114 может быть вычислено Y на основе Y с использованием выражения Y=C.Y.C (задающего процедуру умножения слева и справа), где С и Спредставляют матрицу преобразования и обратную ей матрицу, соответственно. Матрица преобразования может представлять собой любую подходящую матрицу, например, ортогональную матрицу. В одном из вариантов осуществления настоящего изобретения матрица преобразования представляет собой матрицу преобразования Фурье, заданную следующим образом:
<mtable columnalign="left"><mtr><mtd><mtable><mtr><mtd><mrow><mi>C</mi><mo>=</mo><mrow><mo>[</mo><mrow><mtable><mtr><mtd><mn>1</mn></mtd><mtd><mn>1</mn></mtd><mtd><mn>1</mn></mtd><mtd><mo>…</mo></mtd><mtd><mn>1</mn></mtd></mtr><mtr><mtd><mn>1</mn></mtd><mtd><mi>α</mi></mtd><mtd><mrow><msup><mi>α</mi><mn>2</mn></msup></mrow></mtd><mtd><mo>…</mo></mtd><mtd><mrow><msup><mi>α</mi><mrow><mi>n</mi><mo>−</mo><mn>1</mn></mrow></msup></mrow></mtd></mtr><mtr><mtd><mn>1</mn></mtd><mtd><mrow><msup><mi>α</mi><mn>2</mn></msup></mrow></mtd><mtd><mrow><msup><mi>α</mi><mn>4</mn></msup></mrow></mtd><mtd><mo>…</mo></mtd><mtd><mrow><msup><mi>α</mi><mrow><mn>2</mn><mrow><mo>(</mo><mrow><mi>n</mi><mo>−</mo><mn>1</mn></mrow><mo>)</mo></mrow></mrow></msup></mrow></mtd></mtr><mtr><mtd><mo>⋮</mo></mtd><mtd><mo>⋮</mo></mtd><mtd><mo>⋮</mo></mtd><mtd><mrow /></mtd><mtd><mo>⋮</mo></mtd></mtr><mtr><mtd><mn>1</mn></mtd><mtd><mrow><msup><mi>α</mi><mrow><mi>n</mi><mo>−</mo><mn>1</mn></mrow></msup></mrow></mtd><mtd><mrow><msup><mi>α</mi><mrow><mn>2</mn><mrow><mo>(</mo><mrow><mi>n</mi><mo>−</mo><mn>1</mn></mrow><mo>)</mo></mrow></mrow></msup></mrow></mtd><mtd><mo>…</mo></mtd><mtd><mrow><msup><mi>α</mi><mrow><mrow><mo>(</mo><mrow><mi>n</mi><mo>−</mo><mn>1</mn></mrow><mo>)</mo></mrow></mrow></msup><msup><mi>a</mi><mrow><mrow><mo>(</mo><mrow><mi>n</mi><mo>−</mo><mn>1</mn></mrow><mo>)</mo></mrow></mrow></msup></mrow></mtd></mtr></mtable></mrow><mo>]</mo></mrow><mtext> </mtext><mi>г</mi><mi>д</mi><mi>е</mi><mtext> </mtext><mi>α</mi><mo>=</mo><msup><mi>e</mi><mrow><mi>i</mi><mn>.2</mn><mfrac bevelled="true"><mi>π</mi><mi>n</mi></mfrac></mrow></msup></mrow></mtd></mtr><mtr><mtd><mrow><msup><mi>α</mi><mrow><mi>n</mi><mtext> </mtext><mi>g</mi><mo>+</mo><mi>h</mi></mrow></msup><mo>=</mo><msup><mi>α</mi><mi>h</mi></msup><mtext> </mtext><mo>∵</mo><msup><mi>α</mi><mi>n</mi></msup><mo>=</mo><mn>1</mn></mrow></mtd></mtr></mtable></mtd></mtr><mtr><mtd><mi>В</mi><mi>ы</mi><mi>р</mi><mi>а</mi><mi>ж</mi><mi>е</mi><mi>н</mi><mi>и</mi><mi>е</mi><mtext> </mtext><mn>9</mn></mtd></mtr></mtable>
s
[0052] После этого, на шаге 116 на основе Y вычисляют моменты, в соответствии с предшествующим описанием. Вычисленные моменты используют в итеративной процедуре реконструкции SCT, которая может представлять собой реконструкцию на основе формы на шаге 118 или реконструкцию на основе координатной сетки на шаге 124. Следует отметить, что упомянутые процедуры реконструкции могут выполняться с использованием любого подходящего способа реконструкции распределения.
[0053] В отношении реконструкции на основе формы на шаге 118, в одном из вариантов осуществления настоящего изобретения на шаге 120 используют эквивалентные окружности и на шаге 122 выполняют идентификацию вершин. Следует, однако, отметить, что для использования при реконструкции вместо окружности может быть выбрана любая форма, в частности любая элементарная геометрическая фигура. Соответственно, упомянутая аномалия может иметь любую геометрическую форму.
c
c
[0054] Упомянутая реконструкция на основе формы, выполняемая на шаге 118, включает в одном из вариантов осуществления настоящего изобретения использование окружностей, которые полностью описываются тремя координатами (x,y,R). В данном варианте осуществления настоящего изобретения количество независимых переменных, деленное на 3, дает максимальное количество эквивалентных окружностей. Соответственно, можно допустить, что исследуемые области 32 (как показано на фиг.3) имеют форму окружностей, которые могут представлять собой, например, пузыри в потоке нефтепродуктов. В соответствии с этим может использоваться начальное допущение или аппроксимация множества окружностей для определения размера и местоположения исследуемых областей 32, которые могут представлять собой воздушные аномалии в потоке жидкости. В одном из вариантов осуществления настоящего изобретения могут предполагаться две σ-области, одна из которых представляет исследуемую область 32 (например, аномалию), а другая - фон. Однако могут быть использованы дополнительные окружности для предположения о дополнительных исследуемых областях 32. Использование предположения об округлых исследуемых областях 32 обеспечивает простую визуализацию, например, визуализацию измерения потока. Затем в данном способе реконструкции применяют вычисленные моменты для сходимости решения, как будет более подробно описано ниже.
[0055] В одном из вариантов осуществления настоящего изобретения реконструкция на основе координатной сетки, выполняемая на шаге 124, включает преобразование в полярную координатную сетку на шаге 126, за которой следует опциональное уточнение координатной сетки на шаге 128, которое может выполняться адаптивно. В частности, преобразование в полярную координатную сетку может осуществляться с использованием следующих выражений:
<mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow><msub><mi>M</mi><mrow><mi>m</mi><mi>n</mi></mrow></msub><mo>=</mo><mstyle displaystyle="true"><mrow><munder><mo>∬</mo><mrow><msub><mi>Ω</mi><mrow><mi>a</mi><mi>n</mi><mi>o</mi><mi>m</mi><mi>a</mi><mi>l</mi><mi>y</mi></mrow></msub></mrow></munder><mrow><msup><mi>r</mi><mi>m</mi></msup><msup><mi>e</mi><mrow><mi>i</mi><mi>θ</mi><mo>.</mo><mi>n</mi></mrow></msup><mi>d</mi><mi>A</mi></mrow></mrow></mstyle></mrow></mtd></mtr><mtr columnalign="left"><mtd columnalign="left"><mrow><msub><mi>M</mi><mrow><mi>m</mi><mi>n</mi></mrow></msub><mrow><mo>(</mo><mrow><mi>R</mi><mo>,</mo><mi>t</mi><mo>,</mo><msub><mi>θ</mi><mn>1</mn></msub><mo>,</mo><mi>Δ</mi><mi>θ</mi></mrow><mo>)</mo></mrow><mo>=</mo><mstyle displaystyle="true"><mrow><munderover><mo>∫</mo><mi>R</mi><mrow><mi>R</mi><mo>+</mo><mi>t</mi></mrow></munderover><mrow><mstyle displaystyle="true"><mrow><munderover><mo>∫</mo><mrow><msub><mi>θ</mi><mn>1</mn></msub></mrow><mrow><msub><mi>θ</mi><mn>1</mn></msub><mo>+</mo><mi>Δ</mi><mi>θ</mi></mrow></munderover><mrow><msup><mi>r</mi><mrow><mi>m</mi><mo>+</mo><mn>1</mn></mrow></msup><msup><mi>e</mi><mrow><mi>i</mi><mi>θ</mi><mo>.</mo><mi>n</mi></mrow></msup><mi>d</mi><mi>r</mi><mi>d</mi><mi>θ</mi></mrow></mrow></mstyle></mrow></mrow></mstyle></mrow></mtd></mtr><mtr columnalign="left"><mtd columnalign="left"><mrow><mi>В</mi><mi>ы</mi><mi>р</mi><mi>а</mi><mi>ж</mi><mi>е</mi><mi>н</mi><mi>и</mi><mi>е</mi><mtext> </mtext><mn>10</mn><mtext> </mtext><mi>и</mi><mtext> </mtext><mn>11</mn></mrow></mtd></mtr></mtable></mrow>
<mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow><mo>=</mo><mfrac><mn>1</mn><mrow><mi>m</mi><mo>+</mo><mn>2</mn></mrow></mfrac><mrow><mo>{</mo><mrow><msup><mrow><mrow><mo>(</mo><mrow><mi>R</mi><mo>+</mo><mn>2</mn></mrow><mo>)</mo></mrow></mrow><mrow><mi>m</mi><mo>+</mo><mn>2</mn></mrow></msup><mo>−</mo><msup><mi>R</mi><mrow><mi>m</mi><mo>+</mo><mn>2</mn></mrow></msup></mrow><mo>}</mo></mrow><mi>Δ</mi><mi>θ</mi><mtext> </mtext><mi>е</mi><mi>с</mi><mi>л</mi><mi>и</mi><mtext> </mtext><mi>n</mi><mo>=</mo><mn>0</mn></mrow></mtd></mtr><mtr columnalign="left"><mtd columnalign="left"><mrow><mo>=</mo><mfrac><mn>1</mn><mrow><mrow><mo>(</mo><mrow><mi>m</mi><mo>+</mo><mn>2</mn></mrow><mo>)</mo></mrow><mrow><mo>(</mo><mrow><mi>n</mi><mi>i</mi></mrow><mo>)</mo></mrow></mrow></mfrac><mrow><mo>{</mo><mrow><msup><mrow><mrow><mo>(</mo><mrow><mi>R</mi><mo>+</mo><mn>2</mn></mrow><mo>)</mo></mrow></mrow><mrow><mi>m</mi><mo>+</mo><mn>2</mn></mrow></msup><mo>−</mo><msup><mi>R</mi><mrow><mi>m</mi><mo>−</mo><mn>2</mn></mrow></msup></mrow><mo>}</mo></mrow><mrow><mo>{</mo><mrow><msup><mi>e</mi><mrow><mi>i</mi><mi>n</mi><mrow><mo>(</mo><mrow><msub><mi>θ</mi><mn>1</mn></msub><mo>+</mo><mi>Δ</mi><mi>θ</mi></mrow><mo>)</mo></mrow></mrow></msup><mo>−</mo><msup><mi>e</mi><mrow><mi>i</mi><mi>n</mi><mo>.</mo><msub><mi>θ</mi><mn>1</mn></msub></mrow></msup></mrow><mo>}</mo></mrow><mtext> </mtext><mi>е</mi><mi>с</mi><mi>л</mi><mi>и</mi><mtext> </mtext><mi>n</mi><mo>></mo><mn>0</mn></mrow></mtd></mtr><mtr columnalign="left"><mtd columnalign="left"><mrow><mtext>Выражение 12 и 13</mtext></mrow></mtd></mtr></mtable></mrow>
[0056] Следует отметить, что выражение 13 может быть вычислено заранее. Может обеспечиваться координатная сетка 140 для преобразования в полярную координатную сетку, показанная на фиг.9 и представляющая в данном варианте осуществления настоящего изобретения поперечное сечение трубы 42 (показанной на фиг.2) или другой структуры в виде трубы. В координатной сетке 140 каждый сегмент 142 представляет пиксель в реконструируемом представлении. В данном варианте осуществления настоящего изобретения имеется 8 октантов, каждый из которых задается 4 сегментами в радиально (R) различных местоположениях. Следует отметить, что после вычисления момента становится известным местоположение каждого из пикселей, представленных сегментами 142.
s
[0057] Таким образом, в одном из вариантов осуществления настоящего изобретения может выполняться процедура 150 преобразования в полярную координатную сетку, как показано на фиг.10. Как можно видеть, изначально предполагается, что два сегмента 142а и 142b имеют значительно меньшую проводимость (а именно, 1 вместо 0). Вклад каждого из сегментов 142 в каждый момент также может быть вычислен заранее. В частности, может быть заранее вычислен коэффициент, связывающий Y с упомянутыми моментами, например, на основе моделирования или конфигурации аномалий. Следует отметить, что предварительные вычисления могут быть проверены с помощью статических экспериментов.
[0058] Для получения первого приближения решают описанную систему линейных уравнений для заранее вычисленных моментов 142 и измеренных моментов. Соответственно, определяют значения σ и назначают их для каждого из сегментов 142. Например, координатная сетка 140 может быть модифицирована для включения единой области 144 фона и двух отдельных сегментов 142а и 142b.
[0059] Для того чтобы лучше задать границы, например, определить границы второй фазы, сегменты 142а и 142b дополнительно делят на секции 146 для уточнения измерений. Следует отметить, что все другие сегменты 142 (кроме сегментов 142а и 142b) больше не используют в вычислениях. Соответственно, для второй фазы определяют моменты для дополнительно разделенных сегментов 142а и 142b.
[0060] В соответствии с этим, в различных вариантах осуществления настоящего изобретения для сетки и ее преобразования предполагают распределение, аппроксимированное с использованием наложения геометрических форм. Например, в одном из вариантов осуществления настоящего изобретения, в качестве аппроксимации вторичной фазы (например, пузырей в нефти или газе) используют окружности.
[0061] В еще одном варианте осуществления настоящего изобретения реконструкция на основе координатной сетки, выполняемая на шаге 124, включает преобразование в декартову координатную сетку на шаге 130, за которой следует опциональное уточнение координатной сетки на шаге 132. Это преобразование выполняют аналогично процедуре преобразования в полярную координатную сетку, описанной выше.
0
1
2
[0062] Таким образом, в различных вариантах осуществления настоящего изобретения предлагается подход с матричным преобразованием, например, EIT на основе SCT, в котором применяют пространственно-частотную реконструкцию с использованием множества значений σ (-σ, σ, σ …). При использовании предварительных вычислений в процедуре реконструкции в различных вариантах осуществления настоящего изобретения не решают дифференциальные уравнения, а решают только алгебраические или полиномиальные уравнения, которые вычисляются быстрее. Соответственно, EIT на основе SCT может применяться при реконструкции для высококонтрастных материалов, например, материалов с объемным потоком.
exp
s
-1
[0063] В одном из вариантов осуществления настоящего изобретения может выполняться способ 160 итеративной реконструкции SCT, показанный на фиг.11. Способ 160 может обеспечивать грубое, однако характерное представление общего распределения аномалий. Способ 160 включает вычисление экспериментальных моментов М на шаге 162, например, с использованием инструмента 54 томографии с использованием «мягкого поля», где M=bY. После этого на шаге 164 оценивают проводимость, которая в одном из вариантов осуществления настоящего изобретения задается как:
<mrow><mi>σ</mi><mo>=</mo><msup><mrow><mrow><mo>(</mo><mrow><msup><mi>S</mi><mi>H</mi></msup><mi>S</mi></mrow><mo>)</mo></mrow></mrow><mrow><mo>−</mo><mn>1</mn></mrow></msup><msup><mi>S</mi><mi>H</mi></msup><msub><mi>M</mi><mrow><mi>exp</mi></mrow></msub><mtext> </mtext><mi>В</mi><mi>ы</mi><mi>р</mi><mi>а</mi><mi>ж</mi><mi>е</mi><mi>н</mi><mi>и</mi><mi>е</mi><mtext> </mtext><mn>14</mn></mrow>
[0064] Следует отметить, что для упомянутой оценки могут использоваться различные формы и координатные сетки, как будет более подробно описано далее. Например, может быть использована координатная сетка любого типа, такая как полярная координатная сетка, треугольная координатная сетка или прямоугольная координатная сетка, среди прочего.
comp
[0065] После этого на шаге 166 вычисляют моменты, которые могут быть определены как M=Sσ. Следует отметить, что S может быть изменена на основе решения точного преобразования, которое уточняют после идентификации исследуемых областей (например, аномалий).
[0066] Затем на шаге 168 определяют разность моментов, вычисленных для текущей итерации, и экспериментально вычисленных моментов. Упомянутая разность задает элемент обновления, используемый для обновления упомянутой оценки на шаге 164 в данной итеративной процедуре. Элемент обновления, вычисляемый на основе упомянутой разности моментов, может задаваться следующим образом:
<mrow><msup><mi>σ</mi><mrow><mi>i</mi><mo>+</mo><mn>1</mn></mrow></msup><mo>=</mo><msup><mi>σ</mi><mi>i</mi></msup><mo>+</mo><mi>α</mi><msup><mrow><mrow><mo>(</mo><mrow><msup><mi>S</mi><mi>H</mi></msup><mi>S</mi></mrow><mo>)</mo></mrow></mrow><mrow><mo>−</mo><mn>1</mn></mrow></msup><msup><mi>S</mi><mi>H</mi></msup><mrow><mo>(</mo><mrow><msub><mi>M</mi><mrow><mi>exp</mi></mrow></msub><mo>−</mo><msub><mi>M</mi><mrow><mi>c</mi><mi>o</mi><mi>m</mi><mi>p</mi></mrow></msub></mrow><mo>)</mo></mrow><mtext> </mtext><mi>В</mi><mi>ы</mi><mi>р</mi><mi>а</mi><mi>ж</mi><mi>е</mi><mi>н</mi><mi>и</mi><mi>е</mi><mtext> </mtext><mn>15</mn></mrow>
comp
exp
[0067] Упомянутый обновляемый элемент (например, ошибочный элемент) затем подают на вход шага 164 оценки, таким образом обеспечивается итеративно обновляемое распределение проводимости. Данную итеративную процедуру выполняют до тех пор, пока не будет достигнута сходимость решения, например, M≈M. Соответственно, сходимость может быть достигнута, когда разность моментов нулевая или равна заранее заданному значению (например, допустимому отклонению).
[0068] Таким образом, в различных вариантах осуществления настоящего изобретения при итеративной реконструкции SCT используют процедуру матричного умножения в соответствии с настоящим описанием. В одном из вариантов осуществления настоящего изобретения используют прямоугольную сетку для увеличения скорости реконструкции, например, приблизительно до 2500 кадров в секунду. Использование сетки другой формы может давать другие скорости реконструкции. Например, при использовании треугольной сетки в некоторых вариантах осуществления настоящего изобретения может обеспечиваться скорость реконструкции, равная около 80 кадров в секунду.
s
s
[0069] Соответственно, с использованием различных вариантов осуществления настоящего изобретения, определяют Y на основе приложенных возбуждений и измеренных откликов. Затем на основе Y может определяться Y. После этого моменты, которые используют при итеративной реконструкции, могут быть определены как Y/b.
[0070] Различные варианты осуществления настоящего изобретения и/или его компоненты, например, модули, элемент или компоненты и контроллеры, могут быть также реализованы как часть одного или более компьютеров или процессоров. Компьютер или процессор может включать вычислительное устройство, устройство ввода, блок отображения и интерфейс, например, для доступа к Интернету. Упомянутый компьютер или процессор может включать микропроцессор. Микропроцессор может быть соединен с шиной связи. Упомянутый компьютер или процессор может также включать память. Память может включать память с произвольным доступом (Random Access Memory, RAM) и память "только для чтения" (Read Only Memory, ROM). Компьютер или процессор может также включать устройство хранения, например привод жестких дисков или привод съемных носителей, такой как привод оптических дисков, привод твердотельного носителя (например, флэш-RAM) и т.п. Устройство хранения может представлять собой также другие аналогичные средства для загрузки компьютерных программ или других инструкций в компьютер или процессор.
[0071] Термин "компьютер" или "модуль" в настоящем описании может включать любую систему на базе процессора или микропроцессора, включая системы с применением микроконтроллеров, компьютеры с ограниченным набором команд (reduced instruction set computer, RISC), специализированные интегральные схемы (application specific integrated circuit, ASIC), программируемые вентильные матрицы (field-programmable gate array, FPGA), графические процессоры (graphical processing unit, GPU), логические схемы и любые другие схемы или процессоры, способные выполнять описанные функции. Указанные выше примеры приведены исключительно для иллюстрации и не ограничивают определение и/или значение термина "компьютер".
[0072] Упомянутый компьютер или процессор выполняет набор инструкций, которые хранятся в одном или более элементах памяти, для обработки входных данных. Упомянутые элементы памяти могут также хранить требуемые данные или другую информацию. Элемент памяти может выполнять роль источника информации или элемента физической памяти в машине обработки данных.
[0073] Упомянутый набор инструкций может включать различные команды, которые инструктируют компьютер или процессор как машину обработки данных выполнять определенные операции, например, способов и процедур различных вариантов осуществления настоящего изобретения. Упомянутый набор инструкций может быть выполнен в виде программного обеспечения, которое может формировать часть материального машиночитаемого носителя или машиночитаемой среды. Упомянутое программное обеспечение может быть выполнено в различных формах, например, в виде системного программного обеспечения или прикладного программного обеспечения. Упомянутое программное обеспечение может быть также выполнено в виде набора отдельных программ или модулей, программного модуля в составе более крупной программы или части программного модуля. Упомянутое программное обеспечение может включать использование модульного программирования в форме объектно-ориентированного программирования. Обработка входных данных упомянутой машиной обработки данных может выполняться в ответ на команды оператора или в ответ на результаты предшествующей обработки, или в ответ на запрос, осуществляемый другой машиной обработки данных.
[0074] Термины "программное обеспечение", "встроенное программное обеспечение" и "алгоритм" в настоящем описании взаимозаменяемы и включают любую компьютерную программу, хранимую в памяти для исполнения компьютером, включая RAM-память, ROM-память, EPROM-память, EEPROM-память и энергонезависимую RAM-память (non-volatile RAM, NVRAM). Указанные выше типы памяти являются лишь примерами и, соответственно, не ограничивают типы памяти, которые могут быть использованы для хранения компьютерных программ.
[0075] Следует понимать, что приведенное выше описание предназначено для иллюстрации, а не для ограничения изобретения. Например, описанные выше варианты осуществления настоящего изобретения (и/или их аспекты) могут быть использованы в сочетании друг с другом. Кроме того, может быть выполнено множество модификаций для приспособления конкретной ситуации или конкретных материалов к принципам различных вариантов осуществления настоящего изобретения в пределах сущности изобретения. Несмотря на то, что описанные размеры и типы материалов предназначены для задания параметров различных вариантов осуществления настоящего изобретения, эти варианты осуществления изобретения не ограничивают изобретение и являются лишь примерами его осуществления. Из данного описания специалистам очевидны множество других вариантов осуществления настоящего изобретения. Объем настоящего изобретения определяется приложенной формулой изобретения и всеми эквивалентами, которые она охватывает. В приложенной формуле изобретения термины "включает" и "в котором" используют как более простые эквиваленты термина "охватывает". Кроме того, в формуле изобретения термины "первый", "второй", "третий" и т.п. использованы только для обозначения и не накладывают числовых ограничений на объекты. При этом ограничительная часть формулы изобретения не представлена в форме "средство плюс функция" и не должна интерпретироваться согласно шестому абзацу § 112, 35 U.S.С, если только в ограничительной части формулы изобретения не использовано в явном виде выражение "средства", за которым следует описание функции без упоминания других структур.
[0076] Выше были приведены примеры для описания различных вариантов осуществления настоящего изобретения, включая предпочтительный вариант, а также для обеспечения возможности специалисту реализовать на практике различные варианты выполнения изобретения, включая создание и использование любых устройств или систем и способов. Объем настоящего изобретения определяется формулой изобретения и может включать другие примеры, которые могут быть очевидны специалистам. Упомянутые другие примеры находятся в рамках формулы изобретения, если они имеют структурные элементы, не отличающиеся от указанных в пунктах формулы изобретения, или включают эквивалентные структурные элементы с несущественными отличиями от указанных в пунктах формулы изобретения.