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

СПОСОБ ПОСТРОЕНИЯ ИЗОБРАЖЕНИЯ АКУСТИЧЕСКИХ СВОЙСТВ ПУТЕМ НЕЛИНЕЙНОГО НИЗКОЧАСТОТНОГО УПРАВЛЕНИЯ СВОЙСТВАМИ ВЫСОКОЧАСТОТНОГО РАССЕЯНИЯ И РАСПРОСТРАНЕНИЯ (ВАРИАНТЫ) И АКУСТИЧЕСКОЕ УСТРОЙСТВО ДЛЯ ОСУЩЕСТВЛЕНИЯ ЭТОГО СПОСОБА Представлены новый способ построения акустических изображений, который обеспечивает получение изображений со сниженным шумом реверберации; получение изображений параметров нелинейного рассеяния и распространения от объекта, а также способы определения коррекций для аберраций волнового фронта, вызванных пространственными изменениями акустической скорости распространения. Эти способы нашли применение в различных областях, таких как акустические изображения геологических структур, получение изображений подводных объектов с помощью гидролокации, а также получение медицинских ультразвуковых изображений. Способы основаны на обработке принятых сигналов от излученных групп акустических импульсов в двухполосном частотном диапазоне с перекрывающимися высокочастотными и низкочастотными импульсами. Высокочастотный импульс используется для реконструкции изображения, тогда как низкочастотный импульс используется для управления свойствами высокочастотного рассеяния и/или распространения высокочастотного импульса. Первый способ использует рассеянный сигнал от одиночной группы импульсов в двухполосном частотном диапазоне для фильтрации в области "быстрого времени" (времени глубины) для построения сигнала с подавленным шумом реверберации, с чувствительностью 1-й гармоники и улучшенной пространственной разрешимостью. В других способах излучают две или более групп импульсов в двухполосном частотном диапазоне, в которых частота, и/или фаза, и/или амплитуда низкочастотного импульса изменяется с каждой излученной группой импульсов. Путем фильтрации координаты номера импульса и корректировки задержек нелинейного распространения, а желательно также и амплитуд извлекают линейный обратно рассеянный сигнал, а также численные параметры нелинейного распространения и рассеяния. Подавленные сигналы реверберации в дальнейшем могут быть использованы для оценки коррекции аберрации волнового фронта, поскольку в них содержатся приблизительные оценки коррекции аберрации. Сигналы также могут быть использованы с широкими излучаемыми пучками для многочисленных параллельных принимаемых пучков для увеличения частоты кадров в 2D- и 3D-изображениях. Нелинейный сигнал полезен при построении изображений разницы в свойствах объекта, микрокальцификаций, в росте волокнистых тканей или пенных клеток или микропузырьков газа, которые наблюдаются при декомпрессии или которые вводят в качестве ультразвукового контрастного вещества. Способы также могут быть использованы для построения изображений в проходящих волнах для реконструкции изображений для томографии и дифракционной томографии путем генерации измеренных данных. 014167 Область изобретения Изобретение относится к способам и системам для построения изображений от пространственных изменений акустических свойств объекта, в частности пузырьков газа и центров рассеяния высокой плотности в этом объекте. Способы находят применения в разных областях с самыми разнообразными объектами, например для построения ультразвуковых изображений биологических тканей и текучих сред, акустических изображений геологических структур и для гидролокационного определения объектов в воде. Описание уровня техники Получение акустических изображений используется в самых разных приложениях, таких как получение медицинских ультразвуковых изображений внутренних органов, гидролокационные изображения рыб, морских животных и других объектов в море; получение изображений геологических структур для разных целей, таких как изучение археологических раскопок и наблюдение за нефтяными скважинами. Для разных применений используют широкий диапазон частот излучаемого акустического импульса: от инфразвука для построения изображений некоторых геологических структур до 100 МГц для построения ультразвуковых изображений некоторых биологических и микроскопических структур. Несмотря на большой выбор приложений и разброс частот, способы построения изображения очень похожи для всех приложений. В заявке на патент авторы в целом назвали эти способы построения изображения, как построение акустических изображений, хотя в большой части приложений, особенно в медицинских приложениях, используют неслышимые ультразвуковые частотыв диапазоне от 20 кГц до 100 МГц. Поскольку эти частоты для построения изображений, таких как медицинских изображений, лежат в ультразвуковой области, также будем использовать термин "построение ультразвуковых изображений", хотя и не ограничивая способы только ультразвуковыми частотами и медицинскими приложениями. Несмотря на широкое использование акустических изображений, изображения, получаемые в настоящее время, имеют шумовой фон, требуют значительных навыков для их интерпретации, а также дают ограниченную количественную информацию об объектах. Это создает проблемы для дифференциации структур объектов и количественной оценки свойств этих объектов. При получении гидроакустических изображений в воде из-за отражения от морского дна часто, например, трудно различать такие объекты, расположенные близко к морскому дну или непосредственно на морском дне, как рыбы или другие морские животные или мины. Также при получении геологических изображений может быть трудным определить состав материалов геологических структур. В медицинских приложениях может быть трудным отличить структуры типа опухолей или атеросклеротические ткани от здоровых тканей. Главные причины этого изложены ниже. Пространственные изменения в линейных акустических свойствах объекта (плотность и коэффициент сжатия) являются физической основой построения акустических изображений. Однако из-за больших изменений акустических свойств в сложных структурах существуют следующие эффекты, которые ухудшают качество изображений:i) границы между материалами со значительно различающимися акустическими свойствами могут привести к настолько сильному отражению акустического импульса, что многократные отражения будут иметь большую амплитуду. Такие многократные отражения называются реверберацией импульса и вносят вклад в распространяющийся акустический импульс в виде "хвоста", который в акустическом изображении проявляется в виде шума;ii) изменения в акустической скорости в сложных структурах объектов приводят к аберрации прямого распространения акустического волнового фронта, приводя к разрушению фокусировки основного максимума пучка и увеличению боковых максимумов пучка. Уменьшенная фокусировка основного максимума пучка путем аберрации волнового фронта уменьшает пространственную разрешающую способность системы построения акустических изображений. И реверберации импульса, и увеличение боковых максимумов пучка путем аберрации волнового фронта вносят кумулятивный шум в изображение, что, в свою очередь, уменьшает отношение самого сильного центра рассеяния к самому слабому, которые могут быть распознаны рядом друг с другом, что и определяет так называемое контрастное разрешение изображения. Этот шум также называется акустическим шумом, поскольку он создан самим излучаемым акустическим импульсом. Увеличение мощности излучаемого импульса не приведет, таким образом, к увеличению отношения сигнал-шум этого типа в отличие от шума электронных приемников. Например, в эхокардиографии шум от реверберации импульса может затенить изображения апикальной области сердца, создавая сложности нахождения апикального тромба и уменьшенного сжатия апикального миокарда. Еще, например, при построении изображения сонной артерии шум от реверберации может затенить нахождение и установление границ тромбоцида. Как и в этих примерах, шум от реверберации импульса ограничивает способности нахождения слабых объектов и дифференциацию небольших различий в контрасте изображений во всех аспектах построения акустических изображений. Получение изображений на 2-й гармонике - способ, который направлен на снижение ухудшения качества изображения за счет реверберации импульса в структурах, расположенных вблизи источника акустических импульсов, как, например, стенка тела человека, поскольку содержание второй гармоники в импульсе нарастает с глубиной и, следовательно, очень низкое, когда импульс проходит сквозь структу-1 014167 ры, расположенные близко с источником, как стенка тела. Однако чувствительность изображений на 2-й гармонике меньше (примерно на 20 дБ), чем изображений на 1-й гармонике, что и ограничивает максимальную глубину изображений, особенно в плотных структурах, таких как геологические структуры, и в биологических объектах, таких как печень, почки, грудные железы, а также при построении изображений скорости потока крови. Для построения 3D-изображений в реальном времени желателен широкий излученный пучок, который покрывается многочисленными параллельными принимаемыми пучками для увеличения объемной скорости изображения. Такие широкие излучаемые пучки на 2-й гармонике трудно получить из-за сниженной амплитуды первой гармоники в широком пучке, что создает трудности для построения изображений на 2-й гармонике с многочисленными параллельными принимаемыми пучками,используемыми для построения 3D-изображений в реальном времени. Это особенно верно для разреженных акустических решеток, в которых ограничено количество элементов, создающих излучаемый пучок. В медицинских приложениях поражения тканей, такие как опухоли и атеросклероз артериальной стенки, влияют на акустические параметры ткани, такие как модуль сдвига, коэффициент объемного сжатия и акустическое поглощение. Изменения этих свойств, главным образом, вызываются ростом пенных клеток, жира или волокнистых молекул соединительных тканей, а также благодаря сегрегации кальция в ткани. Рост соединительных тканей увеличивает акустическое поглощение и модуль сдвига, причем последний вызывает увеличенную жесткость к пальпации, что можно наблюдать путем прикосновения к ткани. Большая работа была проделана по оценке модуля сдвига путем использования объемных ультразвуковых волн для определения смещения волн сдвига в ткани, способами, часто называемыми эластографией, а также удаленной ультразвуковой пальпацией. Однако на сегодняшний день эти способы нашли ограниченное клиническое применение, и все еще существует необходимость улучшенной дифференциации таких изменений ткани под воздействием ультразвука. На сегодняшний день в опухолях грудных желез сегрегированные микрокальцификации обнаруживают с помощью рентгеновской маммографии, которые являются индикаторами злокачественных опухолей. Эти микрокальцификации настолько малы, что отраженный от них ультразвуковой сигнал пропадает в сигнале от окружающей их ткани и они не могут быть обнаружены существующими способами построения ультразвуковых изображений. Таким образом, существует необходимость в улучшении способов построения ультразвуковых изображений, чтобы также можно было обнаруживать такие микрокальцификации. Микрокальцификации в атеросклеротичной бляшке также дают информацию о ее стабильности и здесь также необходимы улучшенные способы построения ультразвуковых изображений таких микрокальцификаций. Некоторые заболевания также влияют на кровоснабжение ткани, например, через ангиогенез или некроз микрососудистой системы в злокачественных опухолях, или уменьшенного потока крови из-за сужения сосудов или тромбоза как в коронарных артериях сердца, так и в периферических сосудах. Скорость потока крови в микрососудистой системе и небольших сосудах так мала, что она не может быть определена с помощью обычных, неинвазивных ультразвуковых методов, основанных на эффекте Допплера. Для улучшения качества ультразвуковых изображений микрососудистой системы, а также для оценки кровоснабжения ткани были разработаны ультразвуковые контрастные вещества в форме растворов небольших микропузырьков (диаметром примерно 3 мкм). Эти микропузырьки вводят в кровяной поток, что приводит к сильно увеличенному и нелинейному рассеянию ультразвука от крови. Они, следовательно, сильно увеличивают нелинейное рассеяние от ткани, содержащей такие микропузырьки,причем в особых случаях в плотных тканях можно наблюдать отдельные микропузырьки, которые создают перспективу построения молекулярного ультразвукового изображения со специальными контрастными пузырьками в ткани. Такие микропузырьки могут обеспечить полезное усиление изображения при их введении в другие жидкости тела, например в межклеточную жидкость для наблюдения за выделением лимфы к сигнальным лимфоузлам, или в мочевыделительной системе для нацеленного присоединения пузырьков к опухолевым тканям, или в других местах. Во время декомпрессии в космосе или при погружении микропузырьки газа часто образуются в ткани спонтанно, вызывая кессонную болезнь, поэтому существует необходимость для раннего обнаружения таких пузырьков газа для улучшения профиля декомпрессии, чтобы избежать кессонной болезни персонала при выполнении таких работ, а также для наблюдения за образованием таких пузырьков, как предупреждение на ранней стадии во время таких работ. Во время добычи нефти из нефтяных скважин часто бывает желательно наблюдать за изменениями в геологических структурах вокруг скважин, например за количеством нефти или газа в песчанике, за границами между нефтью, газом и водой, а также за структурными сдвигами вблизи скважины. Акустические свойства структур и, в частности, нелинейные компоненты акустических свойств зависят от количества газа, нефти или воды в пористой горной породе. Акустические изображения структур, окружающих нефтяную скважину, могут быть получены с помощью акустических преобразователей в нефтяной скважине. Использование способов построения акустических изображений, которые обеспечивают численные акустические данные от объекта, позволяет, следовательно, определять количество газа, нефти или воды в структурах, окружающих нефтяную скважину.-2 014167 При обнаружении рыбы, или морских животных, или других объектов, находящихся вблизи морского дна, часто бывает трудно различать между отраженными сигналами от морского дна и от объекта,особенно при боковых пучках. У промысловой рыбы или морских животных имеются наполненные газом плавательный пузырь или легкие, акустические свойства которых сильно отличаются от акустических свойств морского дна, и при этом нелинейны. Эти различия в акустических свойствах вместе со способами в соответствии с настоящим изобретением могут быть использованы для дифференциации накладывающихся отраженных сигналов от таких объектов и от морского дна. Эти способы могут быть также использованы для усиления небольших твердых структур, таких как мины, на более мягком морском дне или в почве, как они используются для обнаружения микрокальцификаций в опухоли. Таким образом, существует потребность в улучшенном способе построения акустического изображения, в котором уменьшен шум изображения и усилен контраст изображения, связанный с изменением свойств объекта Способы, выполненные в соответствии с настоящим изобретением, уменьшают шум изображения и усиливают контраст изображения, связанный с изменением свойств объекта, путем излучения на объект групп акустических импульсов в двухполосном частотном диапазоне с перекрывающимися высокочастотными и низкочастотными импульсами. Акустические импульсы в двухполосном частотном диапазоне использовались ранее в построении медицинских ультразвуковых изображений для различных целей, где для кардиовмешательства в М-режиме и в допплеровском (Британский журнал по болезням сердца, январь 1984 г.,51 (1), с. 61-69) одновременном излучении использовались импульсы частотой 3 и 1,5 МГц с фиксированными фазовыми соотношениями между импульсами для построения оптимального изображения сердца в М-режиме (импульс 3 МГц) и для допплеровского измерения скорости крови(импульс 1,5 МГц). Использовался преобразователь с концентрическим кольцевым устройством, в котором 3-МГц ультразвуковой импульс в М-режиме излучался и принимался центральным диском преобразователя, тогда как 1,5-МГц допплеровский ультразвуковой импульс излучался и принимался окружающим его кольцевым элементом. Использование излучаемых импульсов в двухполосном частотном диапазоне для улучшенного обнаружения ультразвукового контрастного вещества, микропузырьков, также описано в патенте США 5410516. В этом патенте описано одновременное излучение двух ультразвуковых импульсов с разными центральными частотами, при этом рассеянные от микропузырьков импульсы содержат как суммарную, так и разностную частоты от излученных частот, полученные при нелинейном рассеянии от микропузырьков, т.е. эти суммарные и разностные частоты используются для обнаружения микропузырьков. Подобное использование импульсов в двухполосном частотном диапазоне для обнаружения ультразвукового контрастного вещества описано в патенте США 6312383, где фаза между двумя частотными полосами изменялась между излучениями. Этот случай можно рассматривать как частный случай описанного в патенте США 5410516, в котором изменение фазы низкочастотного импульса можно рассматривать как биение между низкой частотой и частотой следования импульсов. Однако, хотя оба из указанных выше патентов используют нелинейное рассеяние импульсов в двухполосном частотном диапазоне для обнаружения контрастного вещества в ткани, оба они не рассматривают нелинейный эффект влияния импульса низкочастотного диапазона на скорость прямого распространения импульса высокочастотного диапазона, который в практической ситуации ограничивает подавление сигнала от ткани относительно сигнала от контрастного вещества. В патентах также не отмечено, как может быть использована информация от нелинейного рассеяния от обычной ткани или других объектов. Накапливаемые нелинейные эффекты прямого распространения создают подобные характеристики сигналов для сильного линейного рассеяния от ткани, как и локальное нелинейное рассеяние от микропузырьков и тканей. Этот эффект затеняет локальное нелинейное рассеяние от микропузырьков и других частей объекта и ограничивает соотношение мощности контрастное вещество-ткань (CTR). Наличие пузырьков газа и микропузырьков в некоторой области также значительно увеличивает прямой накапливаемый эффект нелинейного распространения и приводит к тому, что линейное рассеяние от объектов, находящихся вне этой области, сильно затеняет рассеяние от пузырьков газа и микропузырьков в объекте. Это явление, например, сильно влияет на построение изображений от контрастного вещества в миокарде с помощью импульсов, которые проходят желудочек сердца с контрастным веществом, прежде чем войдут в миокард, и могут поэтому ошибочно определить кровоток в ишемическом миокарде. Это также влияет на дифференциацию между газом и нефтью за областью с газом в геологических структурах. Настоящее изобретение отличается от предшествующего уровня техники тем, что в нем использованы вызванные низкочастотным импульсом нелинейные эффекты скорости распространения высокочастотного импульса, и тем, что понимание этих эффектов использовано для формирования сигналов для построения изображений, основанных на высокочастотных прошедших и рассеянных сигналах. Это позволяет разделить накапливаемый нелинейный эффект, влияющий на сигналы, и локальные нелинейные параметры объекта, позволяя тем самым определить локальные нелинейные параметры объекта, что было невозможно в предшествующем уровне техники. Это изобретение также предлагает способ разделения накапливаемого влияния акустического поглощения на сигналы, позволяя определить параметры-3 014167 локального акустического поглощения объекта. Сущность изобретения Способы находят применение в построении акустических изображений как с помощью обратно рассеянных сигналов, так и с помощью компьютерной реконструкции, основанной на измерении углового рассеяния и/или прямого пропускания. Группы двухполосных акустических импульсов с компонентами импульсов как в низком, так и в высоком частотных диапазонах, перекрывающихся во временном пространстве, излучаются в направлении области, где находится объект, изображение которого надо построить. В процессе формирования сигнала, используемого для построения изображения, использовано нелинейное управление посредством низкочастотного импульса рассеивающих свойств объекта и свойств распространения высокочастотного импульса. Высокочастотные компоненты принятых сигналов обрабатывают для построения параметров сигналов изображения, а низкочастотные компоненты принятых сигналов могут быть удалены путем фильтрации, например непосредственно в решетке приемного преобразователя. Обработка сигналов согласно изобретению является частью полной обработки, необходимой для построения изображений, причем дополнительная обработка, необходимая для формирования окончательного изображения и не описанная в настоящем изобретении, является широко известной. Для построения изображения с помощью обратного рассеяния используют способы для формирования радиальных линий изображений, тогда как 2D- и 3D-изображения получают путем поперечного сканирования пучком согласно известным способам. Радиальные линии изображений могут быть огибающей сигнала для структурных изображений, допплеровских измерений радиального смещения и скоростей смещения центра рассеяния, деформаций радиального сдвига или скоростей деформаций сдвига соответствующих перемещений центра рассеяния или "быстрого времени" (времени глубины) спектральных параметров для построения характеристики объекта. Параллельные излученные и/или полученные пучки могут быть использованы для построения многочисленных радиальных линий изображения параллельно для ускорения частоты кадров. С помощью реконструкции изображения посредством компьютерной томографии(КТ) способы обеспечивают улучшенные измерения для указанной реконструкции с уменьшенным шумом реверберации импульса и нелинейными параметрами изображения, которые дают дополняющую информацию. Изобретение предусматривает несколько способов для улучшенного построения изображений с увеличенным числом импульсов, необходимых для формирования импульса, с дополняющим снижением частоты кадров изображения. Изобретение предусматривает также устройство для работы двух или более способов и процедур для оптимального выбора способов, обеспечивающих наилучшую работу устройства при заданных ограничениях, таких как частота кадров, качество изображения, их комбинация и т.д. В первом способе в соответствии с изобретением высокочастотный импульс распространяется на фоне отрицательного пространственного градиента низкочастотного импульсного колебания, так что спад высокочастотного импульса имеет более высокую скорость распространения, нежели фронт импульса. Это приводит к кумулятивному пространственному сжатию высокочастотного импульса по мере его распространения в объект, увеличивая частоту и полосу пропускания (т.е. укорачивая длину) высокочастотного импульса, вдобавок к нелинейному самоискажению высокочастотного импульса, которое приводит к появлению гармонических компонентов в импульсе. Это увеличение частоты, определяемое уменьшением длины импульса, противодействует понижению центральной частоты импульса из-за частотно-зависимого поглощения в объекте, обеспечивая, следовательно, получение более высокой принимаемой центральной частоты, нежели когда этот способ не используется. Поскольку амплитуда низкочастотного импульса сильно уменьшена после первого отражения, многочисленные рассеянные импульсы не будут иметь то же сжатие (укорочение) длины, вызванное нелинейным влиянием низкочастотного импульса на скорость распространения высокочастотного импульса,а из-за поглощения их частота упадет до более низких частот, чем частоты рассеянных импульсов первого порядка с тем же запаздыванием распространения, и, следовательно, могут быть отфильтрованы, приводя к отчетливому подавлению шума, вызванного реверберацией импульса (многократное рассеяние),как и при построении изображений с помощью второй гармоники, но с чувствительностью как с 1-й гармоникой, обеспечивая получение изображений с большей глубины и использование более высоких акустических частот, чем при построении изображений с помощью второй гармоники, улучшая тем самым пространственное разрешение. Также проще получить более широкие излучаемые пучки с возможностью использования более параллельных принимаемых пучков по сравнению с получением изображений с помощью второй гармоники, обеспечивая более высокую частоту кадров для 2D- и особенно для 3D-изображений. Это особенно верно, когда для излучаемого пучка используют разреженные решетки,поскольку в них трудно обеспечить достаточно большие амплитуды для адекватного гармонического самоискажения импульса из-за ограниченного числа элементов решетки. Когда для построения изображений используется 2-я гармоническая полоса принимаемого сигнала, также интересно принять во внимание нелинейное сжатие импульса, поскольку сдвиг частоты вниз из-за поглощения нейтрализуется,обеспечивая более высокие частоты и более короткие импульсы при глубоких диапазонах с улучшенным разрешением. Изобретение также использует этот тип сжатия импульсов для наблюдаемых импульсов-4 014167 сдвига объекта от излучательной силы осуществляющих нажим импульсов, для разделения частот между наблюдаемыми и осуществляющими нажим импульсами. Изобретение также использует расположение высокочастотного импульса близко к максимуму низкочастотного импульса излучаемой группы импульсов, чтобы обеспечить более высокую излучаемую амплитуду высокочастотного импульса с ограничением в механическом индексе (МИ) объекта. Во втором способе в соответствии с изобретением излучают две или более групп акустических импульсов в двухполосном частотном диапазоне последовательно для каждой радиальной линии изображения, причем высокочастотный импульс располагают близко к максимуму или к минимуму низкочастотного импульса, при этом частоту, и/или фазу, и/или амплитуду низкочастотного импульса изменяют с каждой излученной группой для того, чтобы нелинейным образом управлять акустическим рассеянием и свойствами прямого распространения объекта для высокочастотных компонентов. Нелинейное управление скоростью прямого распространения вместе с этим способом используют в процессе формирования сигналов для построения изображения. С помощью этого способа можно, например, сформировать 1-й сигнал для построения изображения, уравнение (14), со значительным подавлением шума реверберации импульса и с чувствительностью 1-й гармоники, чтобы иметь те же самые преимущества, что и одиночного импульса, описанного выше. Изобретение также предлагает определение задержек нелинейного распространения, которые дают 1-й численный нелинейный параметр изображения, уравнение (27), являющийся нелинейным параметром прямого распространения, в виде комбинации дифференциальных и оценочных задержек нелинейного распространения, а также определение амплитуды низкочастотного импульса. Частота этого импульса может быть выбрана настолько низкой (обычно 1/5-1/20 от частоты высокочастотного импульса), что различиями в поглощении акустической мощности между различными объектами и отдельными предметами можно пренебречь, а амплитуду низкочастотного импульса при этом можно определить путем моделирования или измерения в водных или нефтяных смесях. Изобретение также предлагает способ определения локального коэффициента поглощения через комбинацию определенных нелинейных задержек распространения в соответствии со вторым способом (а также в соответствии с 3- и 4-м способами, описанными ниже), а также радиального градиента нелинейных задержек распространения, центральной частоты в высокочастотном принятом сигнале согласно первому способу и радиального градиента указанной центральной частоты. Уменьшенный шум реверберации в принимаемых сигналах способствует в соответствии с изобретением оценке корректировок для аберрации волнового фронта, как, например, описано в патентах США 6485423 и 6905465, а также в заявке на патент США 1094387, на которые дана ссылка в связи с настоящим изобретением. Изобретение также дает приблизительную оценку корректировок задержки для аберрации волнового фронта, полученную из нелинейных задержек распространения, как оценки для сигналов от каждого элемента или субапертурного сигнала, определенного в описании ниже. В соответствии с изобретением далее в процессе формирования сигналов для построения изображений корректируют принятые высокочастотные сигналы, используя оценки нелинейных задержек распространения. Затем можно сильно подавить линейно рассеянный в этом процессе от объекта сигнал и получить 2-й сигнал изображения, уравнения (19) и (28), который является нелинейно рассеянным сигналом, отображающим локальные нелинейные свойства объекта в масштабе, меньшем, чем высокочастотная длина волны, тогда как нелинейные параметры распространения показывают нелинейные свойства объекта в масштабе, большем пары высокочастотных длин волн. Затем нелинейный сигнал обеспечивает контраст изображения к быстрым изменениям структур объекта с улучшенной дифференциацией этих структур. Нелинейное рассеяние особенно сильно на границах между материалами с большой разницей в податливости, таких как границы между биологическими тканями и более жесткими тканями,такими как соединительная и мускулатурная ткани, или твердые материалы, такие как кальцификации,или между частями объекта с низкой и высокой податливостью, такими как жир и микропузырьки газа в биологических тканях, газ в пористых горных породах, плавательные пузыри рыб или легкие морских животных и окружающая вода и ткань, улучшая тем самым характеризацию структур объекта. Пузырьки газа часто естественным образом присутствуют в объекте, например в плавательном пузыре рыбы или легких морского животного, пузырьки газа в пористой горной породе, микропузырьки,спонтанно образованные в ткани во время декомпрессии, или микропузырьки, введенные в объект в качестве контрастного вещества. Компрессионная динамика пузырьков в акустических волнах давления описывается дифференциальным уравнением, показывающим резонансное акустическое рассеяние с частотно зависящим запаздыванием фазы между падающей и рассеянной волнами, в противоположность рассеянию от твердых тел или текучих сред, где изменением с частотой этой фазы можно практически пренебречь. Резонансная частота и, следовательно, это запаздывание фазы рассеянного сигнала от пузырька газа для высокочастотного импульса также управляются низкочастотным импульсом, вдобавок к амплитуде сигнала, что позволяет выделить большую часть рассеянной мощности из пузырьков газа посредством второго способа (не только гармонических компонентов) и значительно увеличивает CNR(отношение контраста к шуму) по сравнению с существующими способами. При хорошем определении корректировок в нелинейных задержках распространения, а также корректировок амплитуды способы в-5 014167 соответствии с изобретением сильно подавляют линейно рассеянный от объекта сигнал и значительно увеличивают COR (отношение контраста к объекту) по сравнению с существующими способами. В отличие от способов обнаружения контрастного вещества, существующих в уровне техники, таких как получение гармонических изображений, инверсия импульса или допплеровская мощность, в способах, выполненных в соответствии с настоящим изобретением, могут быть использованы более высокие акустические частоты по сравнению с резонансной частотой пузырька с улучшенным пространственным разрешением. Также можно использовать более низкие амплитуды импульса (меньший механический индекс (МИ, в результате чего не разрушаются пузырьки, являющиеся медицинским контрастным веществом. Это важно при построении изображения определенных заданных в ткани микропузырьков, где только ограниченная группа микропузырьков пристает к выбранным тканям, например опухолевым тканям, атеросклеротичным бляшкам, тромбам, и т.д., где важно получить изображение пузырьков без их разрушения. Облако микропузырьков в биологической ткани, текучих средах или пористых горных породах имеет сильное нелинейное влияние на скорость распространения проходящего через него импульса, причем в этих случаях чрезвычайно важно обеспечить корректировки для нелинейной задержки распространения, чтобы хорошо подавить линейно рассеянный сигнал от объекта вне облака микропузырьков. С помощью этих корректировок задержки изобретение предлагает разделение между накопленной нелинейной задержкой прямого распространения и локальным нелинейным рассеянием, в отличие от того что предлагается в других способах, таких как гармоническое получение изображений или получение изображений с помощью инверсии импульса, и обеспечивает огромное преимущество для подавления сигнала для построения изображения от объекта, когда получают изображение от пузырьков газа, находящихся вне облака, например в периферическом миокарде при построении медицинских изображений,обнаружении газа вне области с газом при построении геологических изображений или обнаружении рыбы или морских животных в области вне стаи рыб или морских животных. Если в этом случае не сделаны никакие корректировки на нелинейные эффекты распространения или эти корректировки ограничены, то линейно рассеянный сигнал от объектов в областях, где отсутствуют пузырьки газа, будет обладать теми же свойствами, что и рассеянный от пузырьков газа сигнал, маскируя тем самым обнаружение пузырьков газа в этих областях. Это может ошибочно указать на кровоток в ишемической области миокарда, на газ в геологических структурах и на рыбу или других морских животных в воде. Помимо многочисленных медицинских применений изображений микропузырьков в качестве ультразвукового контрастного вещества, получение изображений микропузырьков в соответствии с этим способом в случаях декомпрессии, встречающихся при работе в космосе и при погружениях, может быть использовано для наблюдения за образованием таких пузырьков, чтобы изучать и разрабатывать профили декомпрессии, а также чтобы иметь раннее предупреждение о появлении кессонной болезни. При дальнейшей обработке данных согласно второму способу высокочастотные сигналы с корректировкой задержки объединяют вдоль координаты номера импульса, чтобы получить 3-й сигнал изображения, уравнения (17) и (29), который является линейно рассеянным сигналом. Линейно рассеянный сигнал имеет то же самое ослабление из-за поглощения мощности, что и 2-й сигнал изображения, который является нелинейно рассеянным сигналом. Посредством комбинации нелинейно и линейно рассеянных сигналов и определения амплитуды низкочастотного импульса, как показано выше, в изобретении приводится 2-й численный нелинейный параметр изображения, уравнение (30), который является численным нелинейным параметром рассеяния. Этот 2-й численный нелинейный параметр представляет собой пространственные флуктуации нелинейных параметров объекта в масштабе, приблизительно меньшем, чем длина волны в высокочастотной полосе, тогда как 1-й численный нелинейный параметр,уравнение (27), представляет собой усредненные по пространству нелинейные параметры объекта в масштабе, приблизительно большем, чем длина волны в высокочастотной полосе. Таким образом, обратное рассеяние и прямое распространение представляют собой, следовательно, два разных численных параметра изображения, которые можно наблюдать для большей информации о характеристиках объекта. Численные нелинейные параметры улучшают тем самым дифференциацию структур объекта, а также способствуют характеризации объекта с помощью этого способа, который весьма полезен при диагностике опухолей и атеросклеротичных бляшек в медицинских приложениях, при получении информации о газе, нефти или воде в геологических структурах, а также при получении информации о количестве и размерах рыбы и морских животных в воде. Калибровка изменений этих численных параметров с температурой также способствует оценке локальной температуры с помощью ультразвука, например, для выбора гипертермичного или гипотермичного лечения опухолей в медицине. Она также дает новые способы количественного определения объема контрастного вещества в ткани, кровотока через ткань, а также относительные объемы газа и нефти в геологических структурах. При наличии неподвижных, временно стационарных объектов можно, например, излучить два импульса с разной частотой, и/или фазой, и/или амплитудой низкочастотных компонентов, затем объединить рассеянные или излученные сигналы от этих импульсов для подавления реверберации импульсов и определить нелинейные параметры объекта. Когда объект и акустический зонд перемещаются относительно друг друга, преимущественно можно излучить более двух импульсов для каждой радиальной ли-6 014167 нии изображения, чтобы адекватно подавить линейно рассеянный сигнал или подавить шум реверберации импульса, связанный с многочисленными импульсами. Например, можно излучить ряд K-импульсов,имеющих одинаковую фазу высокочастотных компонентов, но с разными частотами, и/или фазами,и/или амплитудами низкочастотных компонентов для каждого импульса. Обратно рассеянные сигналы от этих импульсов объединяют в фильтре верхних частот с поимпульсным накоплением, в котором подавляются реверберации импульса и проходят компоненты 1-го порядка рассеянного сигнала. С помощью оценок и корректировок для нелинейных задержек распространения, а также желательно и корректировок для амплитуды, выполненных перед фильтром верхних частот, можно выделить локальный нелинейно рассеянный сигнал от объекта или рассеянный сигнал от пузырьков газа, а также численные нелинейные параметры распространения и рассеяния объекта. В этом втором способе шум реверберации импульса (и в малой степени сами нелинейные компоненты сигнала) вносят ошибки в корректировки нелинейных задержек распространения. Эти ошибки ограничивают подавление линейно рассеянного сигнала при оценке нелинейно рассеянного сигнала. Для эффективного устранения влияния шума реверберации импульса на оценку корректировок нелинейной задержки можно использовать 2-ю гармоническую компоненту рассеянных сигналов со вторым способом в соответствии с изобретением или же можно в соответствии с 3-м способом изобретения излучить по меньшей мере три импульса с разными частотами, и/или фазами, и/или амплитудами низкочастотного импульса, как описано со ссылкой на уравнения (42)-(44). Этот третий способ также проявляет влияние нелинейного рассеяния на определение корректировок нелинейных задержек. В четвертом способе в соответствии с изобретением, описанным со ссылкой на уравнения (45)-(48),можно излучить по меньшей мере четыре импульса с разными частотами, и/или фазами, и/или амплитудами низкочастотного импульса, обеспечивая определение корректировок нелинейных задержек, линейно рассеянного сигнала и нелинейно рассеянного сигнала с минимальной помехой друг от друга и от шума реверберации. Также можно использовать 2-ю гармоническую полосу (или любую гармоническую полосу) принятого сигнала для формирования изображения согласно любому из второго, третьего и четвертого способов вместо использования первой гармонической полосы принятого сигнала, как описано выше. Преимуществом этого является дополнительное подавление шума реверберации импульса как в сигналах для построения изображения, так и в оценке нелинейной задержки распространения, как было описано выше. Недостатком является меньшая чувствительность при более глубоких диапазонах, тем самым уменьшается глубина изображения для той же частоты изображения. С помощью электронного управления направлением пучка обычно используют то же направление пучка и излучают фокус для всех излученных импульсов для каждой радиальной линии изображения и диапазона глубины, причем полученные сигналы объединяют для подавления линейно рассеянных от объекта сигналов для этой линии изображения. Типичные используемые схемы фильтров включают КИХ-фильтры или фильтры с переменным временем отклика импульса, как ортогональная декомпозиция с использованием, например, полиномов Лежандра с фильтрацией, осуществляемой вдоль номера координаты импульса для каждой глубины. С помощью механического управления направлением пучка как с кольцевой решеткой или получением 3D-изображений обычно излучают импульсы с изменением в частоте, и/или фазе, и/или амплитуде низкочастотного импульса с непрерывной разверткой по направлению пучка, подавая сигнал для каждой глубины на фильтр верхних частот вдоль номера координаты импульса. Выходы фильтров верхних частот затем запоминаются для каждой глубины и радиальной линии изображения для оценки сигналов и параметров изображения, которые будут использованы для реконструкции изображения вдоль радиальной линии изображения в этом диапазоне глубины. Изобретение также предлагает основные конструкции устройств для построения изображений, которые работают в соответствии со способами согласно изобретению. Поскольку число импульсов для каждой радиальной линии изображения вместе с качеством изображения и информацией увеличивается с увеличением порядка способа, частота кадров уменьшается с увеличением порядка способа. В наиболее продвинутой модели устройство может работать, используя более чем один способ с процедурами для оптимального выбора способов для построения изображения с наилучшей эффективностью, принимая во внимание заданные ограничения. Типичными ограничениями являются минимальная частота кадров,минимальные требования к качеству изображения и т.д. Наконец, изобретение предлагает процедуру конструирования решеток преобразователей, которая сводит к минимуму нелинейное влияние низкочастотного импульса на задержку распространения высокочастотного импульса. При малых амплитудах (50 кПа) компонентов низкочастотного импульса такие решетки преобразователей обеспечивают построение изображений медицинских ультразвуковых контрастных веществ или микропузырьков при декомпрессии с ограниченным, но все-таки интересным подавлением линейно рассеянных сигналов от объекта, без корректировки нелинейных задержек распространения высокочастотного импульса, привнесенных низкочастотным импульсом.-7 014167 Краткое описание чертежей Фиг. 1 изображает первый тип излучаемых импульсов в соответствии с изобретением, причем этот первый тип содержит как низкочастотный, так и высокочастотный импульс, при этом высокочастотный импульс расположен на пространственном градиенте низкочастотного импульса. Фиг. 2 показывает, как полосовая пропускающая фильтрация по глубине принятого сигнала из импульса, показанного на фиг. 1, может быть использована для значительного подавления шума реверберации импульса. Фиг. 3 изображает второй тип излучаемых импульсов в соответствии с изобретением, содержащих как низкочастотный импульс, так и высокочастотный импульс, в котором высокочастотный импульс,например, расположен в максимуме либо положительного, либо отрицательного периода низкочастотного импульса. Фиг. 4 изображает задержки прямого распространения высокочастотного импульса, которые созданы низкочастотным импульсом, изображенным на фиг. 3. Фиг. 5 изображает ряд принятых высокочастотных сигналов от последовательных излученных импульсов как функцию быстрого времени (глубины) и медленного времени (номера координаты импульса). Фиг. 6 изображает принятые линейные и нелинейные частотные линии вдоль частотной координаты медленного времени. Фиг. 7 а и 7b показывают, как реверберации импульса испытывают меньшее управление нелинейным распространением посредством низкочастотного импульса, чем рассеянные сигналы первого порядка. Фиг. 7 с изображает выигрыш от обработки сигнала, зависящий от глубины, полученный благодаря способу подавления реверберации импульса в сигнале. Фиг. 8 показывает узел решетки преобразователей для излучения низкочастотных и высокочастотных компонентов. Фиг. 9 изображает основную решетку преобразователей и принцип создания устройства для одновременного измерения пропускания и углового рассеяния в объекте. Фиг. 10 показывает блок-схему блока оценки для сигналов и параметров изображения, которые могут быть получены этим способом. Фиг. 11 показывает еще одну блок-схему блока оценки для сигналов и параметров изображения, которые могут быть получены этим способом. Фиг. 12 показывает блок-схему блока устройства для построения изображения от рассеянных сигналов в соответствии с изобретением. Фиг. 13 показывает блок-схему блока устройства для реконструкции томографического изображения из измерений пропускания и углового рассеяния в соответствии с изобретением. Фиг. 14 показывает расположение акустических преобразователей в нефтяной скважине для наблюдения за геологическими структурами вокруг нефтяной скважины. Описание вариантов выполнения изобретения Акустические объемные волны в однородных материалах в линейном режиме подчиняются линейному волновому уравнению, где скорость с 0 распространения объемной волны определяется плотностью 0 и коэффициентом объемного сжатия 0 однородной среды распространения. Коэффициент объемного сжатия в линейном приближении объемной упругости определяется через относительный объем сжатия материала, как где V - относительный объем сжатия небольшого объема V, подверженного воздействию давления p;- смещение частицы в материале, так что - - относительное сжатие объема. Рассеяние акустических объемных волн от объектов осуществляется пространственными флуктуациями коэффициента сжатия и плотности объекта. Обозначим пространственно меняющиеся плотность и коэффициент сжатия для амплитуд низкого давления соответственно как 0(r) и 0(r), где r является пространственной координатой. Тогда линейный коэффициент обратного рассеяния от некоторой локальной точки r будет равен где 0a(r) и 0a(r) - пространственно усредненные плотность и коэффициент объемного сжатия в масштабе приблизительно пары длин волнакустического импульса;es - единичный вектор в направлении наблюдения рассеянной волны, как показано на фиг. 9. Скалярное произведение eies=cosis, где is - угол между направлением падающей волны и направлением рассеянной волны. Для обратного рассеяния is= и eies=-1. Волновое число падающей волныk=/с=2/, где- угловая частота; с - акустическая скорость распространения. Линейный обратно рассеянный сигнал в точке r от волны давления с амплитудой 1(r, ) при угловой частотепропорциона-8 014167 лен k20(r)1(r, ). Изображения обычно получают с помощью излученного пучка с центральной частотой 1 и полосой В 1, где сигнал изображения - прошедший полосой фильтр вариант k20(r) в координатахr диапазона около 2k1=21/с и полосой 2 В 1/с. Поскольку более мягкий материал (больший коэффициент объемного сжатия) обычно имеет меньшую плотность, члены с коэффициентом объемного сжатия и плотностью в уравнении (2) обычно имеют противоположные знаки, так что для обратного рассеяния эти величины складываются арифметически, причем член с коэффициентом объемного сжатия превышает при обратном рассеянии член с плотностью примерно в 2,5 раза. Пространственные изменения в средней плотности и коэффициенте объемного сжатия, 0a(r) и 0a(r), приводят к изменениям в акустической скорости распространения, т.е. Это пространственное изменение скорости распространения ответственно за аберрации волнового фронта, в частности, встречающиеся в биологических изображениях в стенке тела, но также и во всех объектах, как грудная железа и лимфатические железы, в которых есть область жира, или соединительной ткани. Большие изменения скорости распространения в геологических изображениях часто встречаются во всем изображении. Типичные значения коэффициента объемного сжатия для мягких тканей и жидкостей 04001012 Па-1 с типичной амплитудой акустического импульса p106 Па, что дает V/V0,410-3. Сжатие объема приводит к увеличению плотности как /rho0=0p0,410-3. Горные породы имеют меньший коэффициент сжатия, тогда как вода и нефтепродукты имеют коэффициент сжатия, сходный с тканями. Аналогично, когда объект подвергнут сжатию, коэффициент объемного сжатия уменьшается, что вместе с учетом акустического поглощения в объекте видоизменяет уравнение (1) [1]: где n=(1+В/2 А)5 - нелинейный параметр, соответствующий традиционно используемым параметрам В и А для нелинейного модуля упругости [1]. Свертка по времени между давлением волнового фронта и h представляет собой зависящее от частоты поглощение акустической мощности в материале. Первый член описывает нелинейный коэффициент объемного сжатия, зависящий от давления, тогда как производная этого члена по давлению при стандартном давлении 0 дает что дает относительное изменение коэффициента сжатия с давлением, так как n/0=-2n0p0. Нелинейные изменения в коэффициенте объемного сжатия, следовательно, равны в 2n10 раз больше, чем нелинейные изменения в плотности, созданные давлением, где для p1 получаемn/0=-2n0p0410-3. Нелинейное изменение плотности и коэффициента сжатия приводит к нелинейному изменению как скорости рассеяния, так и скорости прямого распространения волны, а изобретение использует эти эффекты для уменьшения шума, связанного с реверберацией импульса, увеличения контрастности изображения для различных структур объекта, микрокальцификаций и пузырьков газа и для создания численных акустических параметров изображения объекта, микрокальцификаций, пузырьков газа и областей,наполненных газом. Далее описаны иллюстративные варианты выполнения изобретения со ссылкой на чертежи. В первом способе в соответствии с изобретением используют сжатие и расширение во времени высокочастотного импульса низкочастотным импульсом для управления центральной частотой и шириной полосы прямо распространяющегося в объекте импульса. Это сжатие во времени создается благодаря зависимости от давления скорости прямого распространения, которая приблизительно может быть выражена как где na и 0a - локальные пространственно усредненные значения по паре длин волн при нулевом давлении, как было определено выше. Для дальнейшего объяснения этого принципа обратимся к фиг. 1 а, на которой изображен излученный импульс, состоящий из низкочастотной компоненты 101 с амплитудой p0 и добавленной высокочастотной компонентой 102 с амплитудой р 1, причем высокочастотная компонента "едет" на отрицательном пространственном градиенте низкочастотного импульса, с центром высокочастотного импульса, расположенного, например, около нуля низкочастотного импульса c. Высокочастотный импульс используют для построения изображения, а низкочастотный импульс удаляют путем фильтрации в приемнике, например в самом преобразователе приемника, или так, как описано ниже.-9 014167 Зависящая от давления скорость распространения вызывает возрастающее с накоплением искажение прямого распространения импульса, определяемое фактическим давлением импульса, которое складывается из давлений низкочастотного импульса и высокочастотного импульса, которые на расстоянии r распространения создает искаженный низкочастотный импульс 103 с искаженным высокочастотным импульсом 104. Искажение высокочастотного импульса может быть разделено на сжатие длины импульса нулевых точек высокочастотного импульса, созданное локальным давлением низкочастотного импульса, и на самоискажение формы импульса, вызванное самим мгновенным значением высокочастотного давления. Номер позиции 105 на фиг. 1b показывает сжатый во времени искаженный импульс (пунктирная линия), тогда как для сравнения показан неискаженный высокочастотный импульс 106 (точечнопунктирная линия), при этом добавление самоискажения импульса дает полностью искаженный импульс 107 (сплошная линия). Сжатие импульса происходит из-за того, что более высокое давление низкочастотного импульса на срезе высокочастотного импульса приводит к более высокому значению скорости распространения среза импульса по сравнению со скоростью распространения с более низким значением низкочастотного давления в передней части высокочастотного импульса. Это сжатие импульса приводит к увеличению центральной частоты и ширины полосы высокочастотного импульса, тогда как искажение формы импульса вносит гармонические компоненты фундаментальной частотной полосы высокочастотного импульса. И то, и другое используются в первом способе в соответствии с изобретением. Искажение распространения низкочастотного импульса вносит гармонические компоненты излученной низкочастотной полосы. Это нелинейное искажение прямого распространения импульса является тем же самым эффектом,что создает гармонические компоненты в прямо распространяющемся импульсе, которое линейным образом обратно рассеивается от объекта и используется в гармоническом построении изображений объектов, описанном далее в отношении уравнений (10)-(14) и фиг. 9. Амплитуда гармонических компонентов в импульсе сначала увеличивается с дальностью распространения, далее затухает с дальностью распространения из-за поглощения акустической мощности высокочастотного импульса и расходимости пучка. Низкочастотная полоса может быть выбрана настолько узкой (1/5-1/20 от ширины полосы высокой частоты), что поглощением низкочастотного импульса практически можно пренебречь для фактических диапазонов изображения. Уменьшение амплитуды низкочастотного импульса тем самым, главным образом, определяется расходимостью пучка, которую можно ограничить путем выбора конструкции решетки; влияние нелинейного распространения низкочастотного импульса на высокочастотный импульс посредством конструкции решетки может быть тем самым сделано достаточно большим во всем диапазоне изображения, а также и в других случаях в соответствии с изобретением, как, например, показано на фиг. 3. Все это обеспечивает повышенную чувствительность при глубоких диапазонах с помощью способов, выполненных в соответствии с изобретением, в сравнении с получением изображений с помощью второй гармоники, явления, к которому вернемся в связи с уравнением (14) и фиг. 7 е. В некоторых случаях использование второй гармонической полосы принятых высокочастотных сигналов вместе со способами согласно изобретению может иметь определенное преимущество, поскольку используется для улучшенного подавления шума реверберации импульса в изображениях и для определения параметров изображения, как это описано ниже. Чтобы прояснить существенные моменты, сначала проанализируем случай, когда отсутствует расходимость пучка (т.е. плоские волны) и нет поглощения ультразвука. Разница давлений низкочастотного импульса на длине волны 1 высокочастотного импульса дается выражением р 0=1 р 0/r. Градиент создает разницу в скорости распространения на длине волны высокой частоты, равную ca=-na0ac0aр 0, которая является причиной сжатия волны с длиной r распространения. Надо отметить, что эта же ca приводит к увеличению низкочастотного градиента давления с r пропорционально уменьшению высокочастотной длины волны, так что падение давления вдоль сжатой высокочастотной длины волны поддерживается равным р 0, как показано на фиг. 1 а. По истечении времениt=r/c0a распространения получаем, таким образом, сжатие высокочастотной длины волны 1=cat=rca/c0a=-n0ap0r. При отсутствии поглощения и расходимости пучка получаем р 0=1 р 0/r=р 000010=2p0010/00, где р 00 - амплитуда, 00=2/00 - волновое число и 00 - длина волны излученного низкочастотного импульса на поверхности решетки до искажения распространения. По истечении времени распространения t=r/c0a получаем, таким образом, сжатие высокочастотной длины волны 1=-n0 аp0r=-n0ap00k0010r, что дает увеличение сжатия на высокой частоте, равное Сжатие/расширение сопровождается пропорциональным изменением ширины полосы импульса(примерно обратной длине импульса). В соответствии с этими формулами частота становится бесконечной для zsh=1/na0ap00k00 - явление, которое называют "акустической встряской", при котором отрицательный градиент низкочастотных колебаний становится бесконечным. В реальности поглощение акустической мощности предотвращает эту ситуацию.- 10014167 Акустическое поглощение также ответственно за скольжение вниз по частоте центральной частоты высокочастотного импульса при сохранении ширины полосы импульса. Увеличение частоты из-за сжатия наблюдается только для уходящего импульса, в котором амплитуда низкочастотного импульса достаточно велика, тогда как скольжение вниз по частоте из-за поглощения наблюдается как для уходящего,так и для рассеянного импульса, в котором для обратного рассеяния получается длина распространения,равная 2r. Для гауссовой огибающей импульса скольжение вниз по частоте дается выражением где В 1 - ширина полосы на 6 дБ высокочастотного акустического импульса;=/8,686-1. Следует отметить, что с помощью построения изображения методом трансмиссионной компьютерной томографии, который описан далее в связи с фиг. 9, получаем только одностороннее распространение, в котором скольжение вниз по частоте из-за поглощения равно половине величины, даваемой уравнением (8). Типичный диапазон построения изображения R для построения изображения с помощью обратного рассеяния ограничен акустическим поглощением, которое возрастает линейно с частотой. Расстояние R,таким образом, соотносится с высокочастотной акустической длиной волны 1=c0a/f1 как R2001-3001. Для р 0=1 МПа, n=5, 0 а=40010-12 Па-1, R=2501 и 10/00=10 получаем f1p=0,458f10, что для f10=10 даетf1p=4,58. "Акустическую встряску" находят для zsh/10=796, что в 2,65 раза больше, чем максимальный диапазон построения изображения. Для В 10=5 получаем скольжение вниз по частоте из-за поглощения высокочастотной центральной частоты, равное f1a=4, в соответствии с уравнением (8), причем полученное скольжение компенсирует увеличение частоты за счет сжатия, так что принимаемая частота для обратного рассеяния первого порядка в поглощающей среде приблизительно равна f10=10 для всех глубин,т.е. как и излученная. Поскольку сжатие импульса создается низкочастотным импульсом, где, как описано выше, для реальных диапазонов построения изображения поглощением мощности можно пренебречь, то сжатие длины импульса и соответствующее ему уширение полосы практически не зависят от поглощения для всех реальных диапазонов построения изображения. Скольжение вниз по частоте центральной частоты высокочастотного импульса из-за поглощения происходит, однако, из-за поглощения высокочастотного импульса. Скольжение вниз по частоте пропорционально квадрату абсолютной ширины полосы В 1 импульса, причем скольжение вниз по частоте из-за поглощения сохраняет неизменной ширину полосы. Следовательно, в примере, приведенном выше, совместным эффектом нелинейного сжатия импульса и преобразования с понижением частоты из-за поглощения является импульс с приблизительно постоянной центральной частотой, которая поддерживает поперечное разрешение (ширину пучка), но с шириной полосы(приблизительно обратной длине импульса), которая увеличивается с глубиной, что улучшает разрешение диапазона с глубиной. Заметим, что из уравнения (8) следует, что увеличение В 1 с глубиной приводит к возрастающему скольжению вниз по частоте с глубиной из-за поглощения. Вернемся к более подробному анализу этого случая при рассмотрении уравнений (35) и (36). Преобразование с повышением частоты высокочастотного импульса, связанное с прямым распространением, может быть использовано для улучшения разрешения изображения при более глубоких диапазонах. Это также может быть использовано для улучшения проникновения с лучшим разрешением при глубоких диапазонах, где, например, можно излучить достаточно низкую частоту, которая повышается до более высокой частоты при более глубоких диапазонах посредством низкочастотного импульса, снижая тем самым общее поглощение вдоль пути импульса для полученной частоты импульса в высокочастотной полосе при глубоких диапазонах. В этом аспекте также может быть использовано скольжение фазы низкочастотного и высокочастотного импульсов с расстоянием распространения, с особыми решениями профиля низкочастотного пучка относительно профиля высокочастотного пучка, как описано ниже в связи с фиг. 8. Это скольжение обеспечивает, например, возможность того, что высокочастотный импульс, расположенный на отрицательном пространственном градиенте низкочастотных колебаний в близких/средних диапазонах, будет скользить к нулю или даже к положительному пространственному градиенту низкочастотных колебаний при далеком диапазоне. Это уменьшает сжатие импульса практически до нуля или даже растягивает высокочастотный импульс путем воздействия на него низкочастотного импульса при глубоких диапазонах, что (вместе со скольжением из-за поглощения) уменьшает центральную частоту и ширину полосы импульса в высокочастотной полосе для глубоких диапазонов с повышенным проникновением. Амплитуда низкочастотного импульса сильно падает при первом отражении, а влиянием низкочастотного импульса на сжатия/расширение высокочастотного импульса можно практически пренебречь после первого отражения, как описано в отношении фиг. 7 далее, тогда как скольжение вниз по частоте высокочастотной центральной частоты из-за поглощения преобладает для всего расстояния распространения многократно рассеянного импульса. Следовательно, расстояние по частоте между рассеянным сигналом 1-го порядка (однократное рассеяние) и шумом реверберации импульса (многократное рассеяние) может в изображении увеличиваться с глубиной, как показано на фиг. 2. На этом чертеже номер- 11014167 позиции 201 показывает ожидаемое изменение центральной частоты первой гармонической компоненты обратно рассеянного сигнала 1-го порядка, при этом ширина полосы B1(r) сигнала увеличивается с глубиной, как показано граничными линиями 202, из-за сжатия высокочастотного импульса низкочастотным импульсом. Центральная частота шума реверберации импульса из-за поглощения уменьшается с глубиной, что изображено линией fre(r), отмеченной номером позиции 203 на этом чертеже. Самые мощные источники, создающие шум реверберации импульса, обычно находятся близко к акустическому источнику, создавая ширину полосы Bre(r) шума реверберации импульса, которая имеет практически ту же самую величину, что и ширина полосы излученного высокочастотного импульса, что проиллюстрировано предельными линиями 204. Вторая гармоническая полоса обратно рассеянного сигнала 1-го порядка благодаря самоискажению прямо распространяющегося импульса (104) показана с центральной частотой 2f1(r) номером позиции 205 и шириной полосы B2(r)B1(r) предельными линиями 206. Амплитуда второй гармонической компоненты сначала возрастает с глубиной, а затем падает благодаря поглощению и расходимости пучка первой гармонической высокочастотной компоненты. Первый сигнал изображения с сильным подавлением шума реверберации импульса в соответствии с изобретением может быть получен с помощью фильтра приемника, который подавляет более низкие частоты с частотой отсечки, которая скользит с глубиной диапазона. На фиг. 2 путем примера показан полосовой фильтр приемника со скользящей центральной частотой fre(r), изображенный линией 207, с шириной полосы Bre(r), которая может меняться с глубиной, как показано граничными линиями 208 на этом чертеже. В ближнем и среднем диапазонах разность частоты между шумом реверберации импульса и первой гармонической полосой рассеянного сигнала 1-го порядка настолько мала, что невозможно разделить две компоненты в частотной области. Однако 2-я гармоническая полоса рассеянного сигнала 1-го порядка быстро возрастает по амплитуде и имеет низкий шум реверберации, так что установка частоты фильтра приемника вблизи второй гармонической полосы для малых и средних глубин, как показано на чертеже, дает принимаемый сигнал с сильным подавлением шума реверберации импульса на этих глубинах. Первая глубина изображения обычно является расстоянием в объект, где обнаруживаются вторые гармонические компоненты излученного высокочастотного импульса. В противном случае фильтр приемника может пропустить первую гармоническую компоненту при совсем малых глубинах и перейти ко второй гармонической компоненте в ближнем и среднем диапазонах изображения. С увеличением z разделение по частоте между первой гармонической компонентой рассеянного сигнала 1-го порядка и шумом реверберации импульса увеличивается, при этом можно осуществлять скольжение фильтра приемника по частоте, а также увеличить его ширину полосы, как показано на чертеже, чтобы включить в частоту принимаемого сигнала компоненты первой гармонической полосы рассеянного сигнала 1-го порядка. Это обеспечивает большую амплитуду принимаемого сигнала, поскольку первые гармонические компоненты мощнее и менее ослаблены с глубиной, чем вторые гармонические компоненты, что также поддерживает чувствительность системы построения изображения для больших глубин. Значение Bre может увеличиваться с глубиной, поскольку ширина первой гармонической полосы импульса увеличивается изза описанного сжатия длины импульса, а чтобы включить как первые, так и вторые гармонические компоненты для увеличения мощности сигнала, для глубоких диапазонов можно предусмотреть уменьшение ширины полосы для снижения шума приемника. Также в решетках преобразователей, в которых высокочастотный импульс скользит от отрицательного к положительному пространственному градиенту низкочастотного импульса, как это было описано выше, т.е. от сжатия импульса к его расширению, ширина полосы импульса уменьшается, что может быть сопоставимо с уменьшенным значением Bre. Следует отметить, что шум реверберации импульса лежит в низкочастотном диапазоне, так что вместо полосового фильтра можно использовать фильтр верхних частот приемника, в котором частота отсечки скользит сz, чтобы включить в сигнал большую часть первой гармонической полосы, поддерживая сильное подавление шума реверберации импульса. С помощью этого способа можно, таким образом, сохранить чувствительность первой гармоники для глубоких диапазонов с одновременным подавлением шума реверберации импульса, как при построении изображений с помощью второй гармоники, обеспечивая более глубокие изображения плотных объектов, таких как печень, почки, грудные железы и т.п., с более высокими частотами и лучшим разрешением. Уменьшение длины высокочастотного импульса за счет сжатия также улучшает разрешение диапазона в изображении, которое в известной степени ранее не наблюдалось. При решении профиля пучка,описанном выше, в котором фаза между низкочастотным и высокочастотным импульсами скользит с глубиной так, что высокочастотный импульс расширяется на больших глубинах, все еще возможно с помощью подходящего решения получить в достаточной степени разделенные шум реверберации импульса и рассеяние 1-го порядка в частотном пространстве, для того чтобы шум реверберации импульса мог быть в значительной степени подавлен с помощью фильтра приемника. Эта модификация указанного способа обеспечивает тем самым более глубокое проникновение, поддерживая высокие частоты для лучшей разрешающей способности в средней и близкой зоне. Следует, однако, отметить, что для некоторых применений преимуществом является использование второй гармонической полосы принимаемого сигнала во всем диапазоне изображения, чтобы получить наилучшее подавление шума реверберации им- 12014167 пульса, а также наилучшее пространственное разрешение изображения. Преимуществом сжатия импульса является, таким образом, улучшенная разрешающая способность диапазона при глубоких диапазонах благодаря сжатию длины импульса, а также возрастающие частоты второй гармоники при глубоких диапазонах, приводя к более узким пучкам и лучшей поперечной разрешающей способности. Улучшенное разделение между второй гармонической полосой и полосой шума реверберации импульса при глубоких диапазонах также улучшает подавление шума реверберации импульса по сравнению со стандартным получением изображения с помощью второй гармоники. Принимаемый высокочастотный сигнал выделяют из общего рассеянного/излученного сигнала(также включая низкочастотные компоненты) путем фильтрации, например непосредственно в преобразователе приемника. Однако искажение нелинейного распространения компоненты низкочастотного импульса, как показано на фиг. 1 а, создает гармонические компоненты низкочастотной полосы и с малым разделением между низкочастотной и высокочастотной полосами, причем эти гармонические компоненты могут создавать энергию из низкочастотного импульса в высокочастотной полосе. Такие гармонические компоненты могут быть удалены путем излучения низкочастотного импульса без высокочастотного импульса и помещения принятого сигнала в высокочастотную полосу. Принятый сигнал затем вычитают из принятых сигналов с высокочастотным импульсом, присутствующих в излученной группе импульсов. Эта процедура может быть использована также и с другими способами в соответствии с этим изобретением, чтобы уменьшить полученную энергию от излученного низкочастотного импульса в высокочастотной полосе. Помещая высокочастотный импульс на положительный пик низкочастотного импульса излученной группы, можно уменьшить механический индекс (МИ) высокочастотного импульса (благодаря меньшей отрицательной амплитуде общей группы импульсов), что позволяет излучать высокочастотные сигналы с большими амплитудами. Это приводит к большему гармоническому искажению высокочастотного импульса, повышает чувствительность гармонического изображения с высокочастотным импульсом, поскольку фильтр приемника, описанный выше, выполнен для выбора гармонических полос высокочастотного импульса для подавления шума реверберации импульса. Должно быть очевидно, что одно направление из каждого направления пучка может излучать более одной группы импульсов, представленных на фиг. 1, с последующей обработкой принятого сигнала от каждого импульса, как на фиг. 2, где указанные обработанные принятые сигналы используются с общеизвестной дальнейшей обработкой для создания сигналов для построения изображения, как, например,для структурных анатомических изображений, допплеровских изображений скорости движущихся центров рассеяния и всех сигналов, которые могут быть из них получены, а также частотного анализа по глубине/времени для характеризации центров рассеяния, и т.п., известного любому специалисту. 2D- и 3D-изображения формируются поперечным сканированием пучком с возможными параллельно излученными и/или приятыми пучками. При этом не существует таких же ограничений на ширину излученного пучка, как с получением изображений с помощью второй гармоники, что позволяет использовать более широкие излученные пучки и более параллельные принятые пучки для более высокой частоты кадров 2D- и 3D-изображений, по сравнению с изображениями, полученными с помощью гармоник. Способ сжатия импульса для высокочастотного импульса также полезен в случаях, в которых излучательная сила ультразвуковых импульсов используется для осуществления локального нажима на объект, например, для измерения сдвиговых деформаций, связанных с модулем сдвига объекта, или для улучшения прикрепления целевых пузырьков - контрастных веществ к выбранным тканям. С помощью настоящего способа можно увеличить разделение по частоте между осуществляющими нажим импульсами и импульсами для осуществления наблюдения путем помещения последних на отрицательный пространственный градиент низкочастотного импульса для увеличения приемной частоты импульсов для осуществления наблюдения. Осуществляющие нажим импульсы могут быть излучены как более длинные импульсы с нулевым низкочастотным импульсом или как последовательность коротких высокочастотных осуществляющих нажим импульсов, помещенных вблизи положительного пространственного градиента низкочастотного импульса для растягивания и преобразования с понижением частоты высокочастотных осуществляющих нажим импульсов. В обоих случаях импульсы для осуществления наблюдения приобретают более высокие частоты, нежели осуществляющие нажим импульсы, и могут быть излучены вскоре после или во время последовательности осуществляющих нажим импульсов, при этом разделение эхо-сигналов между осуществляющими нажим импульсами и импульсами для осуществления наблюдения осуществляют путем фильтрации в глубоком (быстром) временном пространстве, как было отмечено выше. Во втором способе в соответствии с изобретением излучают две или более групп импульсов с частотными компонентами в низкочастотной и высокочастотной полосах, которые перекрываются во времени и в которых амплитуда, и/или фаза, и/или частота низкочастотных импульсов меняется от импульса к импульсу. Этот способ дает другой тип принимаемого сигнала с сильно подавленным шумом реверберации импульса, аналогично предыдущему способу в соответствии с изобретением, а также обеспечивает построение изображений нелинейных параметров рассеяния в объекте, в особенности получение изображений микрокальцификаций, пузырьков газа и областей, наполненных газом, а также дает численные- 13014167 нелинейные параметры рассеяния и распространения объекта. Начнем описание этого способа со ссылки на фиг. 3 а, на которой показан излученный импульс, состоящий из низкочастотной компоненты 301 и высокочастотной компоненты 302, причем высокочастотная компонента "едет" на положительном гребне низкочастотного импульса с амплитудой р 0. Высокочастотный импульс используют для построения изображения, а низкочастотный импульс удаляют путем фильтрации в приемнике. Поскольку эти два импульса распространяются через объект вместе, амплитуда низкочастотного импульса влияет на коэффициент рассеяния объекта для высокочастотных компонентов посредством нелинейного изменения в коэффициенте сжатия и плотности как Член с нелинейным коэффициентом сжатия теперь в 2n2,525 раз больше, чем член с нелинейной плотностью, причем два этих члена обычно имеют противоположные знаки. Построение изображений с помощью высокочастотных импульсов с центральной частотой 1 дает, таким образом, полосовой фильтр в координате диапазона этого параметра в частотном диапазоне около 2k1=21/c, как было описано выше. Далее, например, излучают 2-й импульс, как показано на фиг. 3b, где полярность низкочастотного импульса 303 обратная по сравнению с импульсом 301, тогда как высокочастотный импульс 304 занимает такое же временное положение в группе импульсов, что и импульс 302, изображенный на фиг. 3 а, так что низкочастотная амплитуда в месте расположения высокочастотного импульса теперь равна -р 0. Поскольку нелинейный параметр рассеяния в уравнении (9) линеен по добавленному давлению, р 0, нелинейный рассеянный сигнал от высокочастотного импульса 304 на фиг. 3b имеет противоположный знак относительно нелинейного рассеянного сигнала от высокочастотного импульса 302, изображенного на фиг. 3 а. В то же самое время на линейный компонент высокочастотного рассеяния от объекта, как показано уравнением (2), не влияет низкочастотный импульс. Из-за нелинейных изменений скорости распространения с давлением согласно уравнению (6) скорость распространения высокочастотного импульса будет изменяться между величинами, соответствующими импульсам, показанным на фиг. 3 а и 3b, как с/с 0 а=2na0ap0. Обратно рассеянный сигнал от высокочастотных импульсов оказывается, таким образом, сдвинут по времени для положительных и отрицательных полярностей низкочастотных импульсов. Скорость с 0 а распространения для мягких тканей,воды и нефтепродуктов имеет средние значения 1,5 мм/мкс. Горные породы имеют примерно в 2 раза более высокие скорости. Время запаздывания обратно рассеянного сигнала от центра рассеяния на расстоянии r определяется как где s - расстояние вдоль оси пучка;t0(r) - время запаздывания для р 0=0;(r) - добавленное нелинейное время запаздывания распространения из-за нелинейного управления скорости распространения высокочастотного импульса низкочастотным импульсом; р 0 - амплитуда низкочастотного импульса в месте расположения высокочастотного импульса как функция глубины. В дальнейшем будем называть (r) временем запаздывания нелинейного распространения или задержкой нелинейного распространения. У высокочастотного импульса также будет наблюдаться накапливающееся самоискажение, как было описано в отношении фиг. 1b, которое увеличивает гармонические полосы высокочастотного импульса для определенного расстояния, сменяющееся их уменьшением из-за поглощения высокочастотного импульса при более глубоких диапазонах. Благодаря низкой частоте импульса низкочастотной полосы(обычно 1/5-1/20 от высокой частоты) запаздывание нелинейного распространения, привнесенное низкочастотным импульсом, будет преобладающим для гораздо больших глубин. Множитель 2 в t0(r) появляется из-за суммирования времени запаздывания распространения исходящего, излучаемого, импульса и времени запаздывания обратно рассеянного импульса. Низкочастотная компонента будет иметь достаточно высокую амплитуду, чтобы влиять на скорость распространения исходящего импульса, и, следовательно, этот множитель 2 отсутствует в (r). Поскольку нелинейное управление временем запаздывания осуществляют на исходящем импульсе, это управление временем запаздывания одинаково для рассеянного сигнала во всех направлениях, а также и в прямом направлении, что является управлением скоростью прямого распространения, к которому вернемся в связи с описанием фиг. 9.- 14014167 Когда фазовые соотношения высокочастотных и низкочастотных компонентов практически постоянны вдоль пучка, это время сдвига будет монотонно меняться с локальной пространственной средней величиной naа, как показано на фиг. 4, на котором номер позиции 401 обозначает +(r), где высокочастотный импульс "едет" на положительном гребне низкочастотного импульса, а номер позиции 402 обозначает -(r), где высокочастотный импульс "едет" на отрицательной долине низкочастотного импульса. Величина (r) обозначает разностную задержку между этими двумя импульсами и показана номером позиции 403. Для амплитуды низкочастотного импульса р 01 МПа получаем c/c0a=na0ap0210-3, что для расстояния R=3001=300c0aT1, где Т 1=l/f1 - период высокочастотного импульса, т.е. получаем добавленное время запаздывания (R) из уравнения (10), равное 0,6T1, т.е. близко к периоду при центральной частоте высокочастотной полосы. Этот диапазон изображений типичен для построения медицинских ультразвуковых изображений, тогда как в воде можно иметь больший диапазон, а в геологических применениях максимальный диапазон для построения изображений может быть меньше этой величины. Следует отметить, что разница в задержке нелинейного распространения для положительного и отрицательного низкочастотного импульса на фиг. 3 в два раза больше этой величины. Следовательно, даже с более низкими амплитудами p0, например до p050 кПа, получаем значительную задержку нелинейного распространения в соответствии с уравнением (10), которая должна быть скомпенсирована для адекватного выделения нелинейно рассеянного сигнала, как описано далее. Для построения изображений с помощью компьютерной томографии в проходящих пучках имеем распространение в одну сторону без уменьшения сигнала в процессе рассеяния, что обеспечивает по меньшей мере двойной диапазон построения изображений с помощью компьютерной томографии в проходящих пучках и, следовательно,двойную конечную величину задержки нелинейного распространения, к которой вернемся в связи с описанием фиг. 9. Теперь опишем, как изобретение устанавливает сигналы изображения путем нелинейного управления низкочастотным импульсом параметрами рассеяния и распространения объекта для высокочастотного импульса в соответствии со вторым способом изобретения. Сначала исследуем случай получения изображений с помощью обратного рассеяния, причем xk(t) обозначим принятый обратно рассеянный сигнал от излученного импульса с номером k, см. примеры на фиг. 5. Время t отражает глубину центров рассеяния и обозначено как быстрое время, тогда как номер k координаты импульса измеряет более медленные изменения в объекте и обозначается как координата медленного времени. Скорость выборки вдоль координаты медленного времени является частотой следования импульса fprf=1/Tprf, где Tprf - временной интервал между излученными импульсами и обычно выбирается немного большим, чем время Tmax, необходимое для получения рассеянного сигнала с самого глубокого диапазона запаздывания изображения. Для построения медицинских изображений Tmax300 Т 1, тогда как для построения изображений в воде получаем большие значения, а для построения геологических изображений обычно получаются более низкие значения Tmax. На чертеже схематически проиллюстрированы принятые сигналы для пяти выборок медленного времени 501-505 как функция быстрого времени. Сигналы изменяются как функция координаты медленного времени за счет следующих эффектов: изменения в p0k низкочастотного импульса как функции от k. В некоторых случаях согласно изобретению амплитуда низкочастотного импульса p0k(-1)k. Это приводит к изменению принятого сигнала в координате медленного времени для фиксированного t, вызванному нелинейным рассеянием и распространением с частотой медленного времени fprf/2, как описано далее; перемещение между центрами рассеяния и решеткой преобразователей в направлении диапазона высокочастотного пучка. Это приводит к допплеровскому сдвигу принятого сигнала в координате медленного времени для фиксированного t; перемещение между центрами рассеяния и решеткой преобразователей в направлении, поперечном направлению пучка. Это явление может, к примеру, наблюдаться при механическом сканировании в поперечном направлении к направлению акустического пучка или при перемещении центров рассеяния, таких как сердечная стенка или протекающая нефть, что приводит к расширению частоты сигнала вдоль координаты медленного времени. Математическая модель обратно рассеянной первой гармонической полосы высокочастотного сигнала xk(t), где низкочастотный импульс переключает полярность между излученными импульсами, т.е.p0k(-1)k, а центры рассеяния при этом перемещаются, может быть записана как где d=-21r/c0 - средний допплеровский сдвиг для центров рассеяния, перемещающихся со средней радиальной скоростью r от преобразователя в каждой ячейке диапазона вдоль пучка; величина ulk(t) - комплексная огибающая линейного обратно рассеянного сигнала;unk(t) - комплексная огибающая нелинейного обратно рассеянного сигнала от высокочастотного импульса с номером k с положительной амплитудой р 0 низкочастотного импульса.- 15014167 Огибающие изменяются с номером координаты k импульса, поскольку центры рассеяния и пучок перемещаются относительно друг друга, а центры рассеяния могут перемещаться с разными скоростями в пределах ячейки диапазона, причем оба эффекта приводят к расширению частоты сигнала в координате медленного времени. Величина (t) - запаздывание нелинейного распространения как функция координаты быстрого диапазона-времени для положительной амплитуды низкочастотного импульса. Для простоты использовали аналитическую форму принятого сигнала, где физический радиочастотный аку. стический сигнал Аналитический сигнал может быть получен из физического сигнала как где Н обозначает преобразование Гильберта сигнала;xk%(t) представляет собой комплексную огибающую сигнала. Вторая гармоническая полоса может быть представлена формулой, похожей на уравнение (11), где угловая частота 21, допплеровская частота 2d, а нелинейный рассеянный сигнал очень мал и им можно пренебречь, кроме случая рассеяния от микропузырьков. Вторая гармоническая полоса подавила шум реверберации импульса, который может помочь в определении задержки нелинейного распространения,к чему вернемся в отношении уравнения (21). Переключение нелинейной задержки, (-1)k(t) между импульсами наиболее сильно проявляется в фазе, поскольку переключение фазы равно (-1)k1(t), по сравнению с переключением огибающихt-(-1)k(t), так как ширина полосы сигнала ограничена. Для наблюдения влияния переключения запаздывания на огибающие разделим огибающие на четную и нечетную функции около t, что позволяет выразить комплексную огибающую принятого сигнала как где верхний индекс е обозначает четные компоненты, а верхний индекс o обозначает нечетные компоненты пооколо t. Четные компоненты не изменяются при смене знака , тогда как нечетные компоненты изменяют знак. Следует отметить, что (-1)k=ехрik=expikTprfprf/2, где prf=2/Tprf представляет собой угловую частоту следования импульса и, следовательно, угловую частоту выборки в координатах медленного времени. Введя это выражение, уравнение (12) может быть далее выражено как Для фиксированного быстрого времени t сигнал в координате k медленного времени состоит из 4 компонентов (частотных линий), как показано на фиг. 6 а, на котором номер позиции 601 изображает частотную линию А, равную (uelk-uonk) cos1(t)eidTprfk, центрированную около среднего допплеровского сдвига d, номер позиции 602 изображает частотную линию В, равную -i(uelk-uonksin1(t)ei(d+prf/2)Tprfk,центрированную около d+prf/2, номер позиции 603 изображает частотную линию С, равную(uenk-uolk) cos1(t) ei(d+prf/2)Tprfk, центрированную около d+prf/2, и номер позиции 604 изображает частотную линию D, равную -i(uenk-uolk) sin1(t)eidTprfk, центрированную около d. Линии получают путем смешивания линейного и нелинейного рассеяния с переключением задержки нелинейного распространения, (-1)k(t) между импульсами. При отсутствии переключения нелинейной задержки, т.е. при =0, нечетные компоненты обращаются в нуль, а четные компоненты равны изначальным огибающим. Линии 602 и 604 исчезают, а вся линейно рассеянная мощность содержится в линии 601, центрированной около d, тогда как нелинейно рассеянная мощность содержится в линии 603,центрированной около d+prf/2.- 16014167 Эффект переключения задержки нелинейного распространения наиболее сильно сказывается на фазе сигнала из-за ограниченной ширины полосы сигнала, поскольку это приводит к смещению частоты со сдвигом prf/2. Он также приводит к смещению частоты со сдвигом prf/2 путем участия в огибающихt-(-1)k(t) через нечетные компоненты огибающих для линейного рассеяния, ulk(t), и нелинейного рассеяния, unk(t), тогда как переключение задержки огибающей не обнаруживает эффекта сдвига частоты на четные компоненты огибающих. Сдвиг задержки в фазе приводит к сдвигу части линейно рассеянной мощности от линии, центрированной около d, к линии 602, центрированной около d+prf/2, представленной четной компонентой ulke. Совместное переключение в фазе и огибающей сдвигает часть мощности от центрированной около d к центрированной около d+prf/2 и обратно к центрированной около d как часть линии 604, представленной нечетной компонентой ulko. Аналогично, переключение нелинейной задержки сдвигает часть нелинейно рассеянной мощности от линии 603 к линии 604 через переключение фазы, представленой четной компонентой uenk, тогда как совместное переключение в фазе и огибающей сдвигает часть мощности от центрированной около d+prf/2 к центрированной около d и обратно к центрированной около d+prf/2 как часть линии 602, представленной нечетной компонентой uonk. Заметим, что с увеличениемамплитуда линии 601 падает cos1(t), a амплитуда линии 602 возрастает sin1(t) и достигает максимума, когда 1(t)=/2, что также дает нуль для линии 601. Это означает, что значительная часть линейно рассеянной мощности путем переключения задержки в фазе уходит от центрирования около d к центрированию около d+prf/2, но некоторая часть линейно рассеянной мощности через нечетные компоненты ulk(t) примешивается обратно с центрированием около d как часть линии 604. Такой же эффект наблюдается для нелинейного рассеяния, в котором с увеличениемамплитуда линии 603 падает cos1(t), а амплитуда линии 604 возрастает sin1(t) и достигает максимума, когда часть нелинейно рассеянной мощности путем переключения задержки в фазе уходит от центрирования около d+prf/2 к центрированию около d, но некоторая часть нелинейно рассеянной мощности через нечетные компоненты unk(t) примешивается обратно с центрированием около d+prf/2 как часть линии 602. Если ячейка диапазона также покрывает движущуюся кровь или другие текучие среды, такие как нефть, то результатом линейного рассеяния от движущейся текучей среды будут дополнительные и обычно более широкие спектры, в которых мощность распределена между линией 605 с начальными допплеровскими сдвигами d=-21r/c, где r - распределение скорости текучей среды, и линией 606 с частотами d+prf/2. Смешивание сигнала от текучей среды с переключением задержки нелинейного распространения подчиняется тем же законам, что и для линейного рассеяния от объекта в линиях 601 и 602. Однако нелинейное рассеяние от текучей среды настолько слабо, что исчезает в шуме. Модель сигнала в уравнениях (11)-(13) включает только рассеянный сигнал первого порядка, в которой исходящий высокочастотный импульс следует за низкочастотным импульсом. С учетом многократного рассеяния, называемого также реверберациями, исходящего высокочастотного импульса, в этой модели получаем некоторые модификации многократно рассеянного сигнала. Амплитуда низкочастотного рассеянного импульса низкая и, следовательно, влиянием его нелинейной задержки на скорость распространения рассеянного высокочастотного импульса, уравнения (6), (10), можно пренебречь. Это особенно важно для реверберации исходящего импульса в стенке тела, где фиг. 7 а показывает пример структуры решетки 701 преобразователей и отражателей стенки тела. Номер позиции 702 показывает сильный отражатель спереди решетки. Излученный импульс следует по пути, отмеченному номером позиции 703, на котором при первом попадании в отражатель 702 импульс частично проходит дальше как 704, а частично отражается как 705. Отраженный импульс затем опять отражается от поверхности преобразователя или других сильных отражателей для создания отраженного импульса 706, который снова частично проходит дальше и частично отражается и т.д. Отраженный изначальный импульс от более глубокого отражателя 707 показан номером позиции 708 на фиг. 7b для положительного низкочастотного импульса и номером позиции 709 для отрицательного низкочастотного импульса, причем эти импульсы имеют разную задержку, определяемую изменениями в задержке нелинейного распространения согласно уравнению (10). Дважды отраженный от рефлектора 702 импульс показан номером позиции 710 для положительного низкочастотного импульса. После первого отражения от 702 амплитуды импульсов уменьшаются, при этом уменьшение амплитуды низкочастотного импульса снижает управление временем запаздывания для импульсов реверберации в соответствии с уравнением (10), по сравнению с прямо распространяющимся импульсом. Дважды отраженный от рефлектора 702 импульс для отрицательного низкочастотного импульса будет, следовательно,лишь незначительно отличаться по задержке от 710, и поэтому он показан на чертеже как 711, перекрывающимся с 710. Реверберации импульса исходящего импульса в стенке тела будут, следовательно, наблюдаться как частотная линия, центрированная около =0, показанная как точечно-пунктирная линия 607 на фиг. 6 а, где также показана линия смешанной задержки, центрированная около prf/2, имеющая очень малую амплитуду. Авторы изобретения теперь в состоянии описать и обсудить, как несколько новых параметров изо- 17014167 бражения могут быть выделены из прошедших и рассеянных акустических сигналов в соответствии со вторым способом изобретения, чтобы сформировать новые акустические изображения структур объекта и скоростей центров рассеяния улучшенного качества. Выделение этих параметров может быть продемонстрировано, основываясь на принятой последовательности xk(t) сигналов. Первый сигнал изображения, выделенный в соответствии со вторым способом изобретения, основан на рассеянном сигнале с коррекцией реверберации, полученный путем полосовой фильтрации принятой последовательности в области медленного времени около prf/2, как, например, показано полосовым фильтром 610 на фиг. 6 а. Этот фильтр сильно ослабляет линию 607 реверберации и выделяет линии 602 и 603, в которых компоненты линейного рассеяния являются доминирующими по сравнению с компонентами нелинейного рассеяния. Путем объединения принятых высокочастотных сигналов, например в фильтре, который ослабляет низкочастотные компоненты медленного времени, пропуская, однако,компоненты медленного времени в полосе, как в традиционном допплеровском методе построения изображений, можно получить группу линейно рассеянных сигналов с коррекцией реверберации Амплитуды этих сигналов монотонно возрастают с . Путем выбора адекватной частоты низкочастотного поля (например, 01/10) можно получить очень низкое ослабление поглощения низкочастотного поля в пределах диапазона изображения, которое незначительно меняется между разными объектами и отдельными предметами в заданных диапазонах глубины. Можно, таким образом, задать низкочастотное поле и амплитуду импульса так, чтобы получить монотонное возрастание (t), как на фиг. 4, чтобы, в свою очередь, получить монотонно возрастающую,изменяющуюся с глубиной функцию выигрыша при обработке сигналов для линейного рассеяния с коррекцией реверберации, показанную номером позиции 712 на фиг. 7 с. Задавая низкочастотное поле так, что 1(Tmax)=/2, где Tmax - максимальный диапазон-время, получаем в Tmax кривую усиления, близкую к максимуму. Эта изменяющаяся с глубиной функция выигрыша при обработке сигналов вместе с контролируемым пользователем изменяющимся с глубиной усилением принимает участие в устройствах построения акустического изображения, снижая необходимость вмешательства пользователя в управление глубиной усиления. Перемещения объекта, такого как миокард,может привести к допплеровским сдвигам, которые могут пройти сквозь фильтр вместе с компонентами в уравнении (14). Это может быть преимуществом, поскольку шум реверберации импульса неподвижной стенки тела сильно ослаблен, перемещение способствует сигналу 1-го порядка от структур объектов проходу сквозь фильтр, улучшая тем самым качество изображения, например, апикальной области сердца. Выходы с фильтра будут также содержать нелинейно рассеянные компоненты линий 602 и 603, а именно которые имеют максимум для =0 и убывают с глубиной при возрастании . Однако нелинейно рассеянные компоненты сигнала в уравнении (15) будут пренебрежимо малы по сравнению с линейно рассеянными компонентами сигнала в уравнении (14). Сигналы после полосового фильтра могут быть использованы для дальнейшей допплеровской обработки для создания допплеровских спектров и радиальных допплеровских линий изображения скоростей центров рассеяния в соответствии с известными способами, в которых полное 2D- или 3Dизображение создают поперечным качанием пучка. Этот способ особенно полезен для допплеровской оценки смещений миокарда и деформаций смещения (радиального градиента смещения), поскольку шум реверберации сильно мешает таким оценкам. Для измерения скорости текучей среды следует отметить,что сигнал помех объекта в уравнении (14) находится около prf/2, и для того, чтобы подавить помехи объекта для оценки сигнала текучей среды, следует использовать либо полосовой задерживающий фильтр в медленном времени около prf/2, либо смещать сигнал по частоте от prf/2 до =0 и использовать при этом стандартный фильтр подавления помех верхних частот перед определением допплеровских частот сигнала текучей среды. Можно использовать несколько различных и перекрывающихся направлений пучков для каждого k,например, как получено непрерывным качанием пучка с помощью кольцевой решетки. Можно также затем использовать БИХ-фильтр для фильтрации в медленном времени (полосовой, нижних частот,верхних частот), где фильтр нижних частот проиллюстрирован в уравнении (77) для аналогичной обработки.- 18014167 Первый сигнал в соответствии со вторым способом изобретения для использования в радиальной линии изображения для усиления линейного обратного рассеяния с подавлением ревербераций импульса может, таким образом, быть получен как огибающая одного из сигналов z где Env представляет собой оператор огибающей;r=ct/2 - диапазон глубины вдоль пучка,при этом формула для z1(t) - пример полосовой комбинации, как в уравнении (19). Полное 2D- или 3D-изображение затем получают поперечным сканированием пучков. Реверберации импульса уменьшаются благодаря накапливающемуся влиянию управления нелинейной скоростью распространения на давление, т.е. тем же эффектом, который формирует гармоническое искажение в прямо распространяющемся импульсе, который используется в получении гармонических изображений. Однако в настоящем способе нелинейное распространение создают низкочастотным импульсом с таким низким поглощением, что чувствительность этого способа аналогична чувствительности для 1-го гармонического построения изображения. Этот факт позволяет использовать более высокие частоты изображения, чем в случае со 2-м гармоническим построением изображения, с улучшенным разрешением, и в особенности позволяет получать лучшие изображения на глубоких диапазонах в плотных объектах, таких как печень, почки и грудные железы. Поскольку излученный пучок в этом способе представляет собой 1-й гармонический пучок, в нем легче получить более широкий излученный пучок, чем в случае со 2-м гармоническим излученным пучком. Это позволяет использовать параллельные принимаемые пучки для увеличения частоты кадров с построением 3D акустических изображений. Для разделения компонентов нелинейного рассеяния и компонентов линейного рассеяния нужно правильно организовать задержку (компенсацию временного сдвига) принятых сигналов, так что смешение частотного сдвига переключения задержки нелинейного распространения исчезает для линейного сигнала, и линия 602 исчезает, при этом остается только линия 603 нелинейного рассеяния в полосе около prf/2. Коррекция задержки зависит от амплитуд и/или от фаз низкочастотных импульсов относительно высокочастотных импульсов и, кроме того, также изменяется с глубиной согласно уравнению (10) и фиг. 4. В соответствии с одним аспектом изобретения можно оценить коррекции задержки k(t) путем максимизации мощности сигнала Обычно не существует реперного сигнала, так что для определения коррекций задержки нужно использовать один из сигналов в качестве реперного, при этом коррекция задержки для этого сигнала становится равной нулю. Следовательно, можно определить только k-1 независимых коррекций задержки в соответствии со способами, описанными ниже. Сумма представляет собой низкочастотный фильтр в области медленного времени с частотной передаточной функцией Можно также использовать другие варианты низкочастотных фильтров, например, обозначенных номером позиции 611 на фиг. 6b. Этого, например, можно достичь путем введения различного веса k в сумму уравнения (17) в качестве КИХ-фильтра, но также в соответствии с известными способами может быть использован БИХ-фильтр, особенно в случае непрерывного качания пучка, как и с кольцевой решеткой. Коррекции задержки k(t) вносят смешение частот, что перемещает всю линейно рассеянную мощность в линию 621 на фиг. 6b, центрированную около =0, поскольку коррекции задержки в уравнении(17) также удаляют допплеровское отклонение d, как описано ниже. Смешение частот коррекций задержки также перемещает нелинейно рассеянную мощность в линию 622, центрированную около=prf/2. Нелинейно рассеянный сигнал может, следовательно, быть получен после коррекций задержки в виде где суммирование представляет собой полосовой пропускающий фильтр около prf/2, который, например, при hk(t)=1 принимает вид тогда как другие значения для hk(t)=1 приводят к модификации этого фильтра, например фильтра, обо- 19014167 значенного номером позиции 610 на фиг. 6b. Линейная комбинация в уравнениях (17) и (19) дает один выходной сигнал в области медленного времени, при этом специалисту очевидно, как изменить уравнения как КИХ- или БИХ-фильтры, которые создают последовательность выходных сигналов в области медленного времени аналогично тому, как описано в связи с уравнениями (14) и (16). В этом случае такие последовательности сигналов в области медленного времени используют для дальнейшей допплеровской обработки сигналов в области медленного времени для формирования изображений допплеровской скорости или деформации смещения или скорости деформации центров рассеяния объектов и текучей среды в соответствии с известными способами. Для оценки коррекций задержки можно в соответствии с одним вариантом выполнения изобретения разделить принимаемый интервал Т времени/глубины на более короткие временные интервалыTi, i=1, , I, которые с вероятностью перекрываются, так что ТiTi, и определить оптимальные коррекции задержки для каждого интервала Ti по отдельности. Мощность сигнала z1(t) в интервале Ti тогда дается функционалом а коррекции задержки определяют путем максимизации J1i по отношению к ip. Примеры процедур максимизации даны ниже. Нелинейное самоискажение высокочастотного импульса, описанное на фиг. 1 в отношении уравнения (6), приводит к падению гармонического преобразования в первой гармонической высокочастотной полосе с расстоянием распространения. Это падение преобразования слегка различается для отрицательных р 0 и положительных p0, приводя к небольшой разнице в амплитуде линейно рассеянного от объекта 1-го гармонического высокочастотного сигнала для положительных и отрицательных p0. Кроме того,неточности в излучающих усилителях создают трудности в излучении точно заданных амплитуд, создавая неточности в амплитуде принятых сигналов. Для улучшенного подавления линейно рассеянного высокочастотного сигнала в уравнении (19) можно организовать возможность изменения hk(t) от k, так что мощность сигнала znc(t) сведена к минимуму в каждом интервале Ti. Такие величины hk могут, например,быть найдены путем минимизации мощности в znc(t) в каждом интервале по отношению к hk, после корректировки задержки сигналов. Это сводится к минимизации следующего функционала: Примеры процедур минимизации даны в уравнениях (55)-(58). Изменения в величинах hk также ответственны за изменения в амплитудах излучения высокочастотных компонентов, когда амплитуда и/или полярность, и/или фаза, и/или частота низкочастотного импульса изменяется от одного излученного импульса к другому. Задержки ip можно затем очень просто получить путем максимизации Jli, тогда как hk получают путем минимизации Jni. Максимизация Jli дает оценку средней коррекции задержки для каждого интервала Ti. Для наилучшей коррекции в соответствии с уравнением (17) следует придать значения этих оценок задержки точке внутри соответствующих интервалов Ti и получить интерполированные оценки коррекции задержки k в каждой точке выборки быстрого времени t между выбранными точками. Выбранными точками могут быть, например, середины интервалов или центр тяжести мощности в принятых сигналах в соответствующих интервалах или аналогичные точки. Могут быть использованы несколько методов интерполяции, такие как линейная интерполяция, сплайн-интерполяция любой степени и Фурье-интерполяция. Линейно рассеянный от объекта сигнал, показанный линиями 605 и 606 на фиг. 6 а, в процессе коррекции задержки перемещается в линию 623, показанную на фиг. 6b, таким же образом, как и линейно рассеянный сигнал от текучей среды, после чего указанный сигнал может быть обработан в соответствии с хорошо известными способами допплеровской обработки сигналов. Для сигнала реверберации в линии 607 смешение частоты в коррекциях задержки приводит к распространению энергии в линию 624, центрированную на -d, и в линию 625, центрированную на prf/2-d. Шум реверберации будет, таким образом, вносить ошибки в оценку коррекций задержки нелинейного распространения максимизацией Jli в уравнении (21), а смешение частотных сдвигов коррекций задержки вносит шум реверберации в той же частотной полосе медленного времени, что и нелинейное рассеяние после коррекций задержки (линия 625), и, следовательно, вносит шум в оценку сигнала нелинейного рассеяния, например, в соответствии с уравнением (19). На сигналы изображения, описанные в уравнениях (24)-(30) ниже, будет тем самым больше влиять шум реверберации импульса, чем сигнал изображения, основанный на уравнениях (14) и- 20014167 Со ссылкой на обсуждение в связи с фиг. 1 и 2 должно быть очевидно, что можно использовать не только первую гармоническую полосу принятого сигнала, но и другие полосы, в частности вторую гармоническую полосу, что видно из анализа уравнений (14)-(22). Вторая гармоническая полоса, главным образом, формируется прямым самоискажением импульса, показанным номером позиции 104 на фиг. 1,который локально рассеивается как линейно, так и нелинейно, как показано в уравнении (9). Поскольку самоискажение мало в многократно рассеянных импульсах, вторая гармоническая полоса принятых высокочастотных сигналов будет обладать значительным подавлением шума реверберации, как описано в отношении фиг. 2. При сильном шуме реверберации для определения задержек нелинейного распространения можно преимущественно использовать вторые гармонические компоненты обратно рассеянных сигналов в уравнении (21). Использование второй гармонической полосы для сигнала в уравнениях (14) и (16) приводит к дополнительному подавлению шума реверберации импульса, в особенности когда один из многочисленных центров рассеяния движется. Использование второй гармонической полосы в уравнении (19) дает вторую гармоническую полосу локального нелинейно рассеянного сигнала. Можно также использовать вторую гармоническую полосу для определения задержек нелинейного распространения, а затем применить коррекции задержки к первой гармонической полосе принятого сигнала для определения первой гармонической полосы нелинейно рассеянного сигнала, как и выше. Способы, описанные в патентах США 6485423 и 6905465, также могут быть полезны при использовании совместно с настоящим изобретением для уменьшения влияния ревербераций импульса на определение нелинейных задержек. Далее будет показан третий, уравнения (42)-(44), и четвертый, уравнения (45)-(48), способы в соответствии с изобретением, в которых эффект шума реверберации импульса подавляют непосредственно при определении задержек распространения. До того как описывать эти способы, следует показать, как можно выделить многочисленные сигналы параметров изображения с помощью второго-четвертого способов изобретения. Максимизация Jli no zlc, как описано уравнением (17), вносит коррекции задержки как для времени запаздывания нелинейного распространения, так и для среднего допплеровского времени запаздывания,вызванного смещением объекта/преобразователя во временном интервале Ti. Причиной этого является то, что суммирующий фильтр, описываемый уравнениями (17) и (18), имеет частотный отклик медленного времени с максимумом в =0, а максимум Jli достигается, когда вся линейная энергия сдвинута и центрирована около =0. При использовании фильтра нижних частот медленного времени с плоским откликом около =0 и сильным ослаблением около -prf/2 можно получить максимум Jli, при котором имеются только задержки нелинейного распространения, а линейно рассеянная энергия линий 602 и 603 сдвинута в пропускающую полосу фильтра, без поправок на допплеровскую задержку. Однако такой фильтр трудно выполнить с ограниченным числом импульсов медленного времени, используемых для каждой радиальной линии изображения, и, следовательно, можно получить наиболее грубую оценку, используя фильтр нижних частот с выраженным максимумом при =0. Оценка коррекции задержки тогда будет представлять собой общую задержку распространения как сумму задержки нелинейного распространения и средних допплеровских задержек смещения (т.е. допплеровской задержки) относительного перемещения центров рассеяния объекта и преобразователя. Как было отмечено в связи с уравнением (17), можно определить K-1 задержек нелинейного распространения, т.е. на одно меньше, чем число сигналов. С тремя сигналами или более чем с тремя сигналами точное определение коррекций задержки для сигнала, смоделированного уравнениями (11)-(13), и фильтра с максимальным частотным откликом при =0 дает коррекции задержки, свободные от ошибок,как полные задержки распространения, которые являются суммой задержек нелинейного распространения и допплеровских задержек где d=2r(t)Tprf/c представляет собой допплеровскую задержку смещения из-за радиального смещения r(t)Tprf центра рассеяния между излученными группами импульсов. Допплеровский фазовый сдвиг и допплеровская частота могут быть найдены из Эта допплеровская оценка интересна при определении радиального смещения (из фазы в уравнении(24 а и скорости (из угловой частоты в уравнении (24b объектов, например миокарда, а также радиальной деформации смещения и скорости деформации в качестве радиального градиента радиального смещения и скоростей центров рассеяния. Для определения допплеровских сдвигов центров рассеяния в шуме помех, как кровь или другие текучие среды, можно сначала использовать фильтр подавления помех верхних частот перед допплеровскими коррекциями, как описано в отношении фиг. 12 и известно из уровня техники. Определение деформации смещения центров рассеяния вдоль радиального направления пучка может быть получено из производной dk(t) по быстрому времени. Аналогично можно получить оценку радиальной скорости деформации центров рассеяния из производной dk(t) по быстрому времени.- 21014167 Задержку нелинейного распространения находят как Поскольку определение ik содержит ошибки, можно уменьшить ошибку определения путем усреднения уравнений (24) и (25) по нескольким соседним значениям k. При поперечном перемещении центров рассеяния или при быстром механическом качании акустического пучка может иметь место характерное изменение i,k от k благодаря обмену структур объекта в пучке от импульса к импульсу, причем для усреднения следует ограничить число (k) импульсов. Поскольку р 0(r) может быть определено из предварительных измерений благодаря низкому поглощению низкочастотного импульса, из задержек нелинейного распространения можно определить первый численный нелинейный параметр изображения, представляющий собой свойства нелинейного прямого распространения в материале. Инкремент в коррекциях задержки между соседними интервалами Ti представляет собой параметр нелинейного прямого распространения, который может быть записан в виде где nia и ia представляют собой пространственные средние по диапазонному интервалу, соответствующему Ti; ap0ik - средняя амплитуда низкочастотной компоненты в том же диапазонном интервале, соответствующая излученному импульсу с номером k. Первый численный нелинейный параметр/сигнал изображения (параметр изображения нелинейного распространения) затем получают из уравнения (26) как Второй сигнал изображения, который используют для построения изображений в соответствии с изобретением, представляет собой огибающую anc(r) нелинейно рассеянного сигнала znc(2r/c) в уравнении(19). Эта огибающая соотносится с параметрами нелинейного рассеяния материала как где n(r) - полоса пропускания, фильтрованная около 2kl, как обсуждалось в связи с уравнениями (2) и (9), и усредненная в поперечном направлении вместе с амплитудой р 0 низкочастотного импульса поперек профиля высокочастотного пучка для диапазона r. Экспоненциальный член описывает ослабление поглощения высокочастотного акустического импульса в объекте и компенсируется регулируемой пользователем компенсацией G(r) усиления времени/глубины в акустическом устройстве. Поглощающий множитель может быть найден путем сравнения anc(r) с огибающей 3-го сигнала изображения, линейно рассеянного сигнала, zlc(t) в уравнении (17) после коррекций задержки, которая соотносится с параметрами линейного рассеяния и акустического поглощения в объекте как где l(r) - полоса пропускания, фильтрованная около 2kl и усредненная в поперечном направлении поперек профиля высокочастотного пучка для диапазона r, как обсуждалось выше. Когда р 0(r) известно, например, путем вычислений или измерений, как было описано выше, можно объединить сигналы в уравнениях (28) и (29), чтобы получить 2-й численный нелинейный параметр/сигнал изображения, численный параметр/сигнал нелинейного рассеяния объекта как Откорректированный временной сдвиг содержит тогда информацию распространения о величинеniaia, усредненной по локальному интервалу Ti, тогда как сигнал znc(t) нелинейного рассеяния содержит информацию о локальных, пространственных флуктуациях величины n в интервале Ti. Можно также использовать интерполированные значения ik вдоль быстрого времени до k(t), с учетом того, что уравнения (26) и (27) представляют собой производную вдоль выборок быстрого времени для получения более гладкого варианта np(t) нелинейного параметра npi изображения. Аналогично, значения npi можно приписать точкам внутри интервалов Ti и интерполировать значения между этими точками для построения изображения аналогично, как описано выше для интерполяции задержки распространения. Способы снижения реверберации импульса, как описано в связи с фиг. 2 и уравнением (14), полезны совместно со способами определения коррекций для аберрации волнового фронта, например, как описано в патентах США 6485423, 6905465 и заявки на патент США 10/894387, для снижения негативного влияния шума реверберации на определение коррекции аберрации. Для коррекции аберрации можно использовать решетку акустических преобразователей с двумерным распределением элементов,при этом коррекции применяют к сигналу каждого элемента перед окончательным суммированием в формирователе пучка либо во многих случаях можно объединить сигналы от соседних элементов в сиг- 22014167 налы субапертуры, в которых коррекции аберрации применены к сигналам субапертуры перед окончательным суммированием пучка, а не непосредственно к сигналам отдельных элементов. Наилучшие результаты в определении коррекций аберрации получают, когда подавление шума реверберации импульса применяют ко всем сигналам элемента (или субапертуры) перед определением коррекций аберрации,причем это подавление может быть также применено к суммированному сигналу пучка, поскольку влияния ревербераций импульса снижены при корреляции между сигналами элементов и суммированным сигналом пучка, как описано в указанных выше патентах. Представляя пространственное изменение акустической скорости распространения в виде с 0 а(r)=с 0+с 0 а(r), где с 0 - постоянная скорость распространения, принятая 1,54 мммкс при расчете задержек формирователя пучка согласно предположению однородного материала, можно аппроксимировать задержку аберрации волнового фронта как где r - вектор положения элемента на поверхности решетки, фактического элемента или субапертуры;rf - вектор положения фокуса пучка и(r, rf) - путь луча от центра r элемента до фокуса rf. Для корректировки аберраций волнового фронта для элементов, расположенных в r, следует организовать задержку излученных импульсов для излученного пучка и принятых сигналов для принятого пучка с коррекцией задержки cor(r)=-ab(r), как описано в патенте США 6485423. Существует сильная корреляция между изменениями плотности и коэффициента сжатия в разных материалах, что согласно уравнению (3) также говорит о корреляции между изменениями в коэффициенте сжатия и в скорости распространения как где с - коэффициент пропорциональности, определяемый экспериментально, причем следует заметить, что naа=np - численный параметр нелинейного распространения в уравнении (27), который теперь определен для каждого сигнала элемента или субапертуры, в зависимости от положения r на поверхности решетки. Параметры n00 и npag - пространственно усредненные параметры для всех элементов и интервалов задержки. Подставляя это выражение в уравнение (31), получаем приблизительную оценку для коррекций аберрации в виде где последняя сумма независима от r. Эта оценка также может быть использована как начальная оценка для процедур, описанных в приведенных выше патентах и заявках на патент. Изменение npi(r) с r главным образом возникает из-за распространения через стенку тела, и, поскольку именно изменения параметра ab(r) с положением элемента или субапертуры в r создают аберрации, можно получить хорошие результаты путем суммирования для интервалов i в уравнении (33) лишь слегка за стенку тела. В этой области близкого поля можно создать такое низкочастотное поле, что p0(s) приблизительно постоянно в стенке тела для каждого элемента, так что этот член можно вынести за пределы интеграла для (r) в уравнении (10). Таким образом, видно, что можно сопоставить задержки коррекции аберрации непосредственно с задержками нелинейного распространения за стенкой тела и вычесть пространственно усредненное значение (r), т.е. ag, по всем элементам, а именно где pag(r) - пространственно усредненное значение низкочастотного поля p0(s) вдоль пути распространения (r, Rb), начиная от положения r элемента через стенку тела толщиной Rb. Как отмечено выше, задержки нелинейного распространения с помощью второго способа могут быть определены исходя из второй гармонической полосы принятых сигналов для уменьшения ошибок,вызванных шумом реверберации импульса. В третьем и четвертом способах, представленных ниже, в соответствии с изобретением шума реверберации импульса избегают непосредственно при обработке сигналов, тогда как в практических случаях использование второй гармонической полосы улучшает точность с помощью этих способов. Из уравнения (8) видно, что из производной центральной частоты принятого высокочастотного сигнала можно оценить локальный коэффициент поглощения (r) как функцию диапазона r. Однако при отсутствии сжатия импульса скольжение частоты вниз может привести к тому, что преобразователь приемника будет ограничивать ширину полосы принимаемого сигнала. Следовательно, путем использования- 23014167 сжатия импульса в уравнении (7) так, что оно компенсирует скольжение частоты вниз из-за поглощения,можно поместить принимаемый сигнал в середину полосы преобразователя приемника, чтобы избежать отсечки ширины полосы сигнала снизу из-за ограничений полосы преобразователя приемника, обеспечивая возможность определения локального поглощения мощности из измерений принятой центральной частоты высокочастотного сигнала и задержек (r) нелинейного распространения, как определено в уравнении (10). Параметр (r) оценивают, например, согласно второму способу, описанному выше, или согласно третьему или четвертому способам, описанным ниже, в дополнение к излучению группы импульсов, как на фиг. 1, где высокочастотный импульс помещают вблизи максимального пространственного градиента низкочастотного импульса. Со ссылкой на обсуждение в связи с уравнениями (7) и (8) в этом последнем случае получаем сжатие длины высокочастотного импульса, равное d1=-na0ap00k0010p0(r)/p00dr на расстоянии распространения dr, где р 0(r) - локальная амплитуда низкочастотного импульса, включая изменения, связанные со сходимостью/расходимостью пучка и небольшим поглощением; р 00 и k00 - соответственно амплитуда и волновое число низкочастотного импульса на поверхности решетки; а 10 - высокочастотная длина волны на поверхности решетки. Принятая центральная частота высокочастотного сигнала соотносится с длиной волны как f1(r)=c0a/1(r), что даетdf1=-f12/c0ad1-0,36B12(r)(r)dr=na0ak0010p0(r)f12dr/с 0 а-0,36B12(r)(r)dr. Ширина полосы импульса соотносится к длине R1 импульса как В 1=c0a/R1. Согласно рассуждениям,приведенным выше, длина импульса сжимается на расстоянии dr распространения какdR1=-na0ap00k00R10p0(r)/p00dr,где R10 - длина импульса на поверхности решетки. Преобразование вниз по частоте из-за поглощения сохраняет ширину полосы сигнала и даетdB1=-B12/c0adR1=na0ap00k00R10p0(r)B12dr/c0a. В результате этого анализа получаем следующие дифференциальные уравнения для f1 и B1 Уравнение (35) может быть проинтегрировано непосредственно, что позволяет определить локальный коэффициент поглощения из измеренной задержки (r) нелинейного распространения, его градиент и градиент центральной высокой частоты f1(r) какd/dr=-na(r)0 а(r)р 0(r)/с 0 а(r) и R10B10=с 0 а(0) и k00c0a(0)=0,центральной угловой частоты низкочастотного импульса. Путем регулировки амплитуды низкочастотного импульса можно избежать модуляции принятого высокочастотного сигнала краями высокочастотной приемной полосы так, что уравнение (35b) выполняется, что позволяет его интегрирование в уравнение (36b), что, в свою очередь, позволяет определить(r) из уравнений (35 а) и (36 а). Нелинейный параметр n становится очень мал для твердых материалов, как и коэффициент сжатия. Таким образом, в частности, на границах между мягкими и более твердыми материалами, например тканями с высокой плотностью соединительных волокнистых молекул, кальцификациями или другими материалами с высокой плотностью, нелинейное рассеяние становится очень сильным. Аналогично получают сильное нелинейное рассеяние на границах в более мягкие материалы, такие как жир, пенистые клетки, и особенно микропузырьки газа в ткани, где нелинейное рассеяние еще больше усиливается, как описано ниже. Это же верно и для границы между твердой горной породой и текучей средой или газом в геологических структурах, плавательным пузырем у рыб или легких у морских животных в воде, мин в земле или на мягком дне моря и т.д. Нелинейное изображение усиливает, следовательно, визуализацию таких структур. Изобретение, таким образом, полезно для визуализации микрокальцификаций в мягкой ткани, например для построения изображений опухолей в грудной железе и других тканях или атеросклеротической бляшки, которые трудно наблюдать с помощью существующих акустических способов построения изображения. Кроме того, при менее резких изменениях в коэффициенте сжатия материала, а именно когда податливость уменьшается с ростом соединительной ткани или когда податливость увеличивается с ростом жировых или пенистых клеток, нелинейные параметры, определенные этими способами, дают увеличенный контраст изображения для изменений в ткани по сравнению с существующими способами построения изображений. Параметры изображения, представленные в уравнениях (27) и (30),позволяют получить численные оценки изменений в ткани.- 24014167 Для пузырьков газа, которые либо естественным образом присутствуют в объекте в виде углеводородов, или декомпрессионные пузырьки при погружении, либо вводятся в объект в виде ультразвукового контрастного вещества, динамика их рассеяния описывается дифференциальным уравнением, которое дает резонансное рассеяние с зависящим от частоты запаздыванием по фазе между падающей и рассеянной волной, совершенно противоположно рассеянию от текучей среды или твердых объектов, где изменение запаздывания по фазе с частотой очень мало. Низкочастотный импульс управляет диаметром микропузырьков (малый диаметр при положительном р 0 и большой диаметр при отрицательном р 0), а следовательно, и резонансной частотой микропузырька. Это, в свою очередь, управляет запаздыванием по фазе рассеянного сигнала для высокочастотного импульса в дополнение к амплитуде рассеянного сигнала. Управление особенно сильно для высокочастотных импульсов вблизи резонансной частоты микропузырьков, как описано в заявке на патент США 10/851820, поданной 21 мая 2004 г. Скорректированный нелинейный сигнал znc(t), например, в соответствии с уравнением (19) будет содержать почти всю высокочастотную рассеянную мощность от пузырьков - контрастного вещества(как линейные, так и нелинейные компоненты). Для рассеивающего объекта, который содержит микропузырьки газа, настоящее изобретение, таким образом, значительно увеличивает CNR (отношение контраста к шуму) по сравнению с существующими способами построения изображений таких пузырьков путем выделения почти всей рассеянной мощности высокочастотного сигнала из микропузырьков (в особенности сильные линейные компоненты, а не только нелинейные компоненты). Коррекции для переключения низкочастотным импульсом задержек нелинейного распространения обеспечивают подавление линейно рассеянной мощности от объекта, приводя к большому COR (отношение контраста к объекту). Способ, таким образом, разделяет эффекты нелинейного прямого распространения и локального нелинейного рассеяния и использует локальное управление изменением частотой фазы рассеянного от пузырьков газа сигнала для получения сильного локального сигнала от пузырьков газа с сильным подавлением локального сигнала от объекта. Это отличается от существующих способов построения изображений при помощи контрастных веществ, в которых нелинейное распространение приводит к накапливаемому влиянию на прямо распространяющийся импульс, которое, к тому же, усиливает линейное рассеяние от ткани в процессе детектирования так, что этот линейный сигнал от ткани маскирует сигнал от микропузырьков (а также и нелинейно рассеянный сигнал от ткани). Аналогичные эффекты имеют место при детектировании пузырьков газа за облаком пузырьков газа в геологических структурах и объектов за стаей рыб с плавательными пузырями или морских животных с легкими. Когда импульс проходит сквозь облако пузырьков газа (также сквозь плавательные пузыри), в них создается увеличенное, накапливаемое запаздывание нелинейного распространения для прямо распространяющегося импульса, а также нелинейное изменение в амплитуде импульса, явление, которое подчеркивает необходимость корректировки задержек нелинейного распространения и изменений в амплитуде импульса для получения хорошего подавления линейно рассеянного сигнала от объекта вне облака. Настоящее изобретение, таким образом, имеет явное преимущество перед известными способами построения изображения с помощью микропузырьков. Например, в гармоническом способе построения изображений увеличенное, накапливаемое гармоническое искажение для импульса, что проходит через облако микропузырьков, проявляет себя как сильные гармонические компоненты в линейном рассеянии от ткани вне облака. Это, например, может привести к сильному гармоническому рассеянию от миокарда для импульса, который прошел через облако контрастного вещества в желудочке сердца, маскируя рассеянный сигнал от контрастного вещества - микропузырьков в миокарде. Этот эффект может ошибочно указать на кровообращение в области миокарда с очень низкой перфузией или с ее отсутствием, а также указывать на газ в геологической структуре без газа. С помощью настоящего способа для рассеяния от миокарда за облаком исключают влияние облака микропузырьков в желудочке сердца на прямо распространяющийся импульс путем коррекции задержек нелинейного распространения. Изобретение разделяет локальное нелинейное рассеяние и накапливаемый эффект нелинейного прямого распространения и,следовательно, препятствует тому, что измерение локального нелинейного рассеяния приведет к таким ошибочным указаниям на отсутствие микропузырьков в миокарде. Относительно нелинейных гармонических способов построения изображений настоящее изобретение может использовать более широкополосный излучаемый импульс, и, следовательно, можно достичь более высокой разрешающей способности диапазона изображения. Кроме того, могут быть использованы более высокие частоты для построения изображений, в результате чего наблюдается значительное увеличение разрешения по сравнению с другими способами построения изображений пузырьков газа как в направлении диапазона, так и в поперечном направлении. Режим работы настоящего изобретения менее чувствителен к амплитуде импульсов, использующихся для построения изображений, по сравнению с гармоническими способами построения изображений. Вместе с указанным подавлением принятого линейно рассеянного сигнала с результирующим увеличением COR это облегчает неразрушающее детектирование и построение изображений с высоким разрешением отдельных пузырьков - контрастного вещества с низким механическим индексом (МИ).- 25014167 В медицинских применениях построение изображений повышенной чувствительности и высокого разрешения ультразвукового контрастного вещества имеет большую перспективу в получении изображений изменений в микрососудистой системе, например при неоангиогенезе или некрозе в опухолях или при сниженном кровообращении в миокарде, где были разработаны некоторые стандартные способы использования времени заполнения контрастным агентом. Численные параметры в уравнениях (27) и(30) дают количественную информацию о плотности контрастного вещества в ткани и, следовательно,обеспечивают улучшенные данные об относительном объеме микрососудистой системы. Способы также имеют применения в оценке относительного объема газа в геологических структурах, плотности и размерах рыбы с плавательным пузырем и морских животных с легкими. Путем разрушения контрастного вещества - пузырьков в некоторой области и измерения времени заполнения можно согласно известным правилам получить численные значения кровообращения в ткани. Построение изображений контрастного вещества высокой чувствительности и высокого разрешения также полезно, чтобы проследить за выделением лимфы для нахождения сигнальных лимфоузлов при хирургическом удалении опухолей. Поскольку уравнения (27) и (30) дают параметры изображения, которые не зависят от поглощения в объекте, можно использовать эти параметры объекта для характеризации объекта (например, содержание жира или соединительной ткани в мягкой ткани, содержание нефти, газа или воды в пористой горной породе или количество и размер рыбы или морских животных). Кроме того, можно определить локальные изменения скорости распространения с температурой, основываясь на экспериментах, для наблюдения за изменениями температуры ткани при термическом лечении больной ткани, например, посредством фокусированного ультразвука высокой мощности (ФУВМ) (HIFU), радиочастотной абляции или криохирургии. Температура может быть отслежена по изменениям в численных параметрах, а также по изменениям в скорости распространения, которые являются причиной времени запаздывания между обратно рассеянными сигналами от изображения к изображению по мере изменения температуры. Радиальные градиенты этого времени запаздывания и определяют локальную температуру. Поскольку низкочастотный импульс не влияет на локальное линейное рассеяние высокочастотного импульса, подразумевается, что изменения амплитуды, и/или фазы, и/или частоты низкочастотных компонентов в излучаемых импульсах, но не тех, что показаны на фиг. 3, могут привести к аналогичному результату в подавлении ревербераций импульса от ближнего источника для получения линейно рассеянного сигнала со сниженным шумом реверберации и в подавлении линейного обратного рассеяния для получения нелинейно рассеянного сигнала, как отмечено выше. Например, низкочастотная составляющая 2-го импульса на фиг. 3b может отсутствовать или высокочастотные импульсы не обязательно должны "ехать" точно на положительном гребне или отрицательной впадине низкочастотных импульсов. Эта универсальность важна, так как фаза между двумя частотными импульсами может возникать из-за дифракции или рассогласованности излучательных поверхностей решеток преобразователей, поскольку низкочастотные и высокочастотные компоненты меняются с расстоянием распространения вдоль пучка. Низкочастотные и высокочастотные полосы излученных импульсов часто так широко отстоят друг от друга, что иногда предпочитают использовать разделенные решетки преобразователей для излучения двух полос импульса. Такие решетки могут быть выполнены в виде концентрических колец с разными резонансными частотами, в которых пучки от решеток автоматически перекрываются, или же решетки могут быть установлены бок о бок с пересечением пучков под острым углом. Когда используют две разные решетки для низкочастотных и высокочастотных компонентов с пучками, которые перекрываются под острым углом, фазовые соотношения между низкочастотными и высокочастотными импульсами могут иметь сильную пространственную зависимость, в зависимости от геометрического расположения и размеров поверхностей двух излучающих решеток. Иллюстративное расположение раздельных низкочастотных и высокочастотных решеток преобразователей, излучающих вдоль z-оси, в соответствии с изобретением показано на фиг. 8 а, на котором номер позиции 801 обозначает высокочастотную решетку, а номер позиции 802 обозначает низкочастотную решетку, состоящую из двух частей на каждой стороне высокочастотной решетки. Чертеж может, например, показывать поперечное сечение линейного или криволинейного расположения решетки, где у-ось - вертикальная ось,перпендикулярная азимутальной плоскости сканирования, т.е. x-z-плоскости, перпендикулярнойy-z-плоскости. Чертеж может также показывать поперечное сечение расположения кольцевой решетки сz-осью в качестве излучающей оси, где номер позиции 801 обозначает поперечное сечение высокочастотной кольцевой решетки, а номер позиции 802 обозначает поперечное сечение низкочастотного кольцевого элемента. Как линейное, так и кольцевое расположения проявляют разные задержки распространения для низкочастотных и высокочастотных решеток, которые следует должным образом принять во внимание при разработке решетки и обработки сигнала, чтобы использовать все преимущество основных физических эффектов, лежащих в основе изобретения. Границы иллюстративного низкочастотного пучка показаны номером позиции 803, тогда как границы высокочастотного пучка показаны номером позиции 804. Следует отметить, что имеется ближняя зона 805 спереди решетки, где имеется ограниченное перекрывание между низкочастотным и высокочастотным пучками, приводя тем самым к особенно сильному подавлению стенкой тела ревербераций исходящего импульса, с помощью способов, описанных в связи с уравнениями (14) и (16) и фиг. 7.- 26014167 Также отметим, что расстояние распространения от низкочастотной решетки до центра рассеяния,находящегося вблизи z-оси, гораздо больше для низкочастотного импульса, чем для высокочастотного импульса, в зависимости от геометрических размеров низкочастотной и высокочастотной решеток. Аксиальное низкочастотное поле в точке 806 на оси z и с угловой частотой 0=ck0 будет равно где Plt - излучаемое давление на поверхности решетки;Rlo(z) - расстояние 807 от наружного края низкочастотной решетки до точки 806 на z-оси; aRli(z) - расстояние 808 от внутреннего края низкочастотной решетки до точки 806. Фазовый член представляет собой среднее запаздывание распространения где i(z) представляет собой запаздывание распространения от низкочастотной решетки до точкиh(z) представляет собой запаздывание распространения от высокочастотной решетки до точки 806;Rho(z) - расстояние 809 от наружного края высокочастотной решетки до точки 806. Кроме запаздывания фазы распространения, также будет наблюдаться изменение знака члена с синусом в уравнении (37), как часть фазы члена Hl(r), которая меняется скачками впри прохождении нуля в этом члене с синусом. Номером позиции 810 на фиг. 8b показана разница в запаздывании фазы между низкочастотным полем и запаздыванием распространения высокочастотного поля, которая может быть выражена как как функция от z для типичной геометрии с высокочастотной апертурой Dh=7 мм и с внутренними и наружными частями низкочастотной апертуры соответственно Dli=10 мм и Dlo=15 мм. Резкое падение запаздывания фазы нав точке 811 (z 6,6 мм) происходит, когда разностная фаза распространения l(z)=0(Rli(z)+Rlo(z/2c0 между наружным и внутренним краями низкочастотной решетки проходит через нуль, что приводит к изменению знака члена с синусом в уравнении (37) с -1 на+1. Амплитуда низкочастотного поля Hl(z) показана не в масштабе в логарифмических единицах номером позиции 812, при этом надо заметить, что нуль амплитуды совпадает со сдвигом нав фазе. Нули поля могут быть найдены, когда разница в запаздывании фазы распространения от наружных и внутренних краев решетки к точке поля есть нечетное число , что дает скачок нав выражении для l(z). Для сравнения номером позиции 813 показана разность в фазе распространения между низкочастотными и высокочастотными решетками, равная Следует заметить, что l(z) повторяет l(z) до скачка в - в точке 811 и повторяет это запаздывание с разницейпосле этого. Из-за большого и зависящего от z запаздывания фазы распространения между низкочастотными и высокочастотными решетками получаем зависящее от z относительное положение между высокочастотными и низкочастотными импульсами. Так, например, высокочастотный импульс, который обычно начинается на верхнем гребне низкочастотного импульса, в ситуации, показанной номером позиции 814, скользит к нижней впадине низкочастотного импульса, в ситуации, показанной номером позиции 815, в которой запаздывание фазы l(z) изменилось на , и т.д. Подставляя давление в уравнение (10), откуда вместе с n=5, 0 а=40010-12-1, с 0 а=1540 м/с и Plt=50 кПа, получаем запаздывание нелинейного распространения, показанное номером позиции 816 на фиг. 8 с. Так как запаздывание фазы l(z) (810) изменяется на нескольковдоль оси z, как показано на фиг. 8b, получаем сильно осциллирующую задержку нелинейного распространения, показанную номером позиции 816, с максимальной задержкой, равной max=2 нс. Если не скорректировать правильно задержку для этой задержки нелинейного распространения, получаем максимальное подавление линейного обратно рассеянного сигнала, равное 2sin(1max)-18 дБ для f1=10 МГц. Для построения изображения контрастного вещества обычно используют более высокое низкочастотное давление р 0, например р 0200 кПа, которое даст max=8 нс в примере выше и уменьшит подавление линейно рассеянного высокочастотного сигнала на 2sin(1max)-6 дБ для f1=10 МГц даже при большом скольжении фазового соотношения между высокочастотными и низкочастотными компонентами,как получается в этом примере. Следовательно, амплитуда 50 кПа для низкочастотного давления слишком мала, но может в особых случаях обеспечить интересные изображения контрастного вещества и других микропузырьков газа. Таким образом, можно в некоторой степени усилить сигнал от контрастного вещества без поправок на задержки нелинейного распространения при условии, что решетки выполнены- 27014167 так, что присутствует осциллирующее изменение задержек нелинейного распространения, как и в случае 816. Надо отметить, что, когда высокочастотный импульс находится вблизи нуля в низкочастотном колебании, эффекты нелинейного рассеяния и прямого распространения малы. Этого эффекта можно избежать путем использования разных фаз между низкочастотными и высокочастотными компонентами в последовательных излучаемых импульсах, что сдвигает пространственное положение, в котором высокочастотный импульс находится вблизи нуля низкочастотных колебаний. Нули в амплитуде Hl(s) низкочастотного поля могут быть сдвинуты в пространстве путем изменения низкочастотного центра между последовательными излучаемыми импульсами. Кроме того, благодаря ширине высокочастотного пучка сигнал получают и снаружи оси, где Hl(s)0, а импульс состоит из частотной полосы, которая усредняет нулевые точки для многих частот. Нулей амплитуды можно избежать путем уменьшения ширины низкочастотной решетки, что, однако, также снизит отношение амплитуды давления к управляющему напряжению решетки. Однако для того чтобы выделить нелинейное рассеяние от других объектов, таких как микрокальцификации в мягкой ткани и газ в пористой горной породе, и получить численные параметры объектов,необходимо использовать более высокие низкочастотные амплитуды, при которых корректировка задержек нелинейного распространения становится обязательной для подавления линейного рассеяния и выделения нелинейного рассеяния. Это также является случаем для использования способа подавления шума реверберации импульса, как описано в отношении фиг. 7. Кроме того, с более слабым осциллирующим изменением задержек нелинейного распространения имеется сильное преимущество корректировки задержек для обнаружения микропузырьков газа даже с малой амплитудой низкочастотного импульса. Часто выбирают такую конструкцию низкочастотных и высокочастотных излучающих апертур,чтобы получить минимальное скольжение задержки между низкочастотными и высокочастотными импульсами исходя из максимизации нелинейного управления параметрами рассеяния и распространения для высокочастотного импульса вдоль всего диапазона изображения. Это необходимо для получения монотонного увеличения задержки нелинейного распространения, как показано на фиг. 4. Чтобы избежать -сдвиг фазы в 810 (и соответствующего нуля амплитуды), следует использовать меньшую ширину низкочастотных элементов (т.е. меньшую разницу между Rlo и Rli), а чтобы снизить скольжение фазы между низкочастотными и высокочастотными импульсами, следует подтянуть низкочастотную излучающую поверхность как можно ближе и даже может быть с перекрыванием высокочастотной излучающей поверхности. Однако для многих применений можно смириться с осциллирующим поведением задержек нелинейного распространения и даже использовать этот факт для получения низких задержек нелинейного распространения, как на фиг. 7 с, или сменить сжатие длины импульса для ближнего и среднего полей на растяжение длины импульса для глубоких диапазонов, чтобы снизить частоту и улучшить чувствительность, как описано в отношении фиг. 2. Чтобы снизить влияние шума реверберации импульса на оценку коррекций для задержек нелинейного распространения, а также в нелинейно рассеянном сигнале, можно в третьем способе в соответствии с изобретением излучить больше двух импульсов более чем с двумя разными амплитудами, и/или фазами, и/или частотами низкочастотного импульса. В качестве первого примера согласно изобретению,в котором отсутствует перемещение между центрами рассеяния и решеткой преобразователей (т.е. нулевое допплеровское смещение), излучают низкочастотные импульсы с амплитудами р 0, 0, -р 0, в которых высокочастотный импульс следует близко к гребню или впадине низкочастотного импульса, как, например, показано на фиг. 3 а и 3b. Принятые высокочастотные сигналы от трех импульсов в интервале Ti могут быть записаны какxn(t) - нелинейно рассеянный сигнал с амплитудой р 0 низкочастотного импульса;r(t) - сигнал реверберации для высокочастотного импульса и-i, - задержка нелинейного распространения для высокочастотного импульса, когда тот передается на положительном гребне с давлением р 0 низкочастотного импульса. Согласно одному варианту выполнения изобретения формируем комбинацию из этих трех сигналов в виде где xL(t, i)=xl(t)-xl(t+i) представляет собой комбинацию линейно рассеянных сигналов.- 28014167 Эти комбинации исключили из z1(t) и z2(t) шум реверберации импульса, и, поскольку амплитуда нелинейно рассеянного сигнала гораздо ниже амплитуды линейно рассеянного сигнала, имеем приблизительно z2(t)=z1(t-i), и, следовательно, для интервала Т[i можно определить коррекцию задержки i путем максимизации следующего функционала: относительно i. Когда i определено, можно вычислить компоненты сигналов линейного и нелинейного рассеяния из уравнения (43), например, путем Фурье-преобразования уравнения (43) в быстром времени, которое приводит к системе линейных уравнений по XL и Xn. Сущность проведенной выше операции заключается в том, что три или более импульсов с тремя или более уровнями управления задержками нелинейного распространения позволяют скомбинировать измерения для получения по меньшей мере двух сигналов, в которых шум реверберации импульса в значительной степени подавлен, причем эти новые сигналы могут быть использованы для определения коррекций задержки. В описанной операции разница в амплитуде соседних низкочастотных импульсов должна быть постоянной, обеспечивая задержку нелинейного распространения для каждого сигнала, которая является целым кратным одной задержки. При постоянном во времени радиальном перемещении между решеткой преобразователей и центрами рассеяния, при котором имеется постоянная допплеровская задержка между принятыми сигналами для каждого импульса, также имеется четвертая неизвестная величина, которую нужно определить, что, в свою очередь, требует излучения по меньшей мере четырех импульсов с разными амплитудами низкочастотного импульса, или же могут быть использованы 5 импульсов с 3 амплитудами низкочастотного импульса, как в уравнении (42), в котором путем максимизации мощности сигнала, таким же образом, как и в уравнениях (21) и (44), определяют комбинированную задержку нелинейного распространения и допплеровскую задержку. Задержка нелинейного распространения и допплеровская задержка могут быть разделены посредством операций, как в уравнениях (24) и(25), в которых детали этих и подобных операций могут быть отработаны в соответствии с изобретением любым специалистом в этой области техники. Тем не менее, сама величина i должна быть определена из комбинаций, наподобие уравнения (44),или из их эквивалентов в форме преобразования Фурье, и, поскольку z1 и z2 содержат как линейно рассеянный, так и нелинейно рассеянный сигнал, последний вносит, хотя и очень незначительную, ошибку в определение коррекции задержки, что, в свою очередь, вносит ошибку в оценку нелинейно рассеянного сигнала тем же самым образом, что и оценки, даваемые уравнениями (17)-(22). Основная причина этой ошибки заключается в том, что для неподвижных объектов имеется четыре неизвестных: линейное рассеяние xl(t); нелинейное рассеяние xn(t); шум реверберации импульса r(t); задержка (t) нелинейного прямого распространения (т.е. неподвижные объекты при отсутствии допплеровской задержки). Так как (t) медленно изменяется с t, можно считать эту величину приблизительно постоянной во временном интервале Ti, как было показано выше. Для более точного определения всех четырех неизвестных надо выполнить по меньшей мере четыре измерения с четырьмя уровнями давления +р 0, +р 0/3,-р 0/3 и -р 0 низкочастотного импульса, чтобы получить Затем можно исключить r(t), как и в уравнении (43), чтобы получить три сигнала z1(t), z2(t) и z3(t) и определить i для синхронизации этих сигналов с максимизацией функционала, аналогичного функционалу в уравнении (44). Влияние xn(t) на оценку i будет меньше, чем в случае трех измерений, но ошибка,которую вносят измерения, все равно не будет полностью исключена. Лучший подход в соответствии с четвертым способом изобретения, который устраняет влияниеxn(t) на определение i, получают при преобразовании Фурье уравнения (45) на интервале Ti, где (t) может быть аппроксимирована постоянной i, которая обеспечивает переопределенную систему линейных уравнений в виде где YT=Y1, Y2, Y3,Y4; ХT=Xl, Хп, R, а А(i) - матрица размером 43, полученная из преобразования Фурье уравнения (43) в соответствии