способ автоматического выделения горизонтов из трехмерного массива сигналов сейсмических данных, способ автоматического выбора горизонтов подповерхностных образований из трехмерного массива сигналов сейсмических данных и компьютеризованный способ автоматического выделения горизонтов пластов из трехмерного массива сигналов сейсмических данных
Классы МПК: | G01V1/28 обработка сейсмических данных, например их анализ для интерпретации, коррекции G01V1/32 преобразование одной записи в другую |
Автор(ы): | Гарольд А.Хилдбранд[US], Афтаб Алам[US], Питер Карагоуниз[GB] |
Патентообладатель(и): | Лэндмарк Графикс Корпорейшн (US) |
Приоритеты: |
подача заявки:
1992-07-06 публикация патента:
27.03.1998 |
Использование: для быстрого составления карт горизонтов из трехмерного массива сейсмических сигналов, согласно которым в ходе пакетной обработки цифровых данных, представляющих сигналы импульсов для каждого сейсмического сигнала указанного массива, определяются "единицы" и "нули". Сущность изобретения: создается битовый массив, в котором "единица" на глубине сигнала означает, что на этой глубине, а также на других глубинах с "единицей" в пределах заранее выбранной глубины на смежных сигналах, существует горизонт. Геологоразведчик в режиме реального масштаба времени выбирает исходную точку, идентифицируя тем самым "единицу" на конкретном сигнале на конкретной глубине. Автоматизированная процедура сканирования находит все другие "единицы" на глубинах смежных сигналов. Каждая найденная "единица" на смежном сигнале рассматривается в качестве исходной точки, и процесс продолжается до тех пор, пока уже не будут обнаруживаться новые "единицы". На основе глубин каждой найденной "единицы" создается карта горизонта. Предлагается также способ, согласно которому битовый массив заменяется сжатым массивом сигналов, в котором каждый бит битового массива заменен числовой величиной, соответствующей точному положению на глубине экстремума атрибута, представленного этим битом. 5 табл., 10 ил.
Рисунок 1, Рисунок 2, Рисунок 3, Рисунок 4, Рисунок 5, Рисунок 6, Рисунок 7, Рисунок 8, Рисунок 9, Рисунок 10, Рисунок 11
Формула изобретения
1. Способ автоматического выделения горизонтов из трехмерного массива сигналов сейсмических данных, при котором горизонт является поверхностью раздела пластов земной коры и представляется в сигналах сейсмических данных, проходящих через поверхность раздела пластов в виде общего атрибута сейсмического импульса, и при котором каждый из сигналов сейсмических данных является последовательностью записей данных, причем каждая запись включает глубину и амплитуду импульса и отделена по глубине от соседней с ней записи на равную величину, предусматривающий создание массива уплотненных сигналов и включающий следующие операции: идентификацию для каждого сигнала сейсмических данных глубин, на которых находятся горизонты, сохранение для каждого сигнала сейсмических данных последовательности числовых величин, отражающих глубины, на которых находятся горизонты, для создания уплотненных сигналов массива уплотненных сигналов, и запись для каждой глубины данных, которая идентифицируется как горизонт, интерполяцию глубин выше и ниже идентифицированной глубины для определения и сохранения в памяти точной интерполированной глубины атрибута импульса, за счет чего интерполированная глубина отражает более точную идентификацию глубины, чем глубина, соответствующая записи глубины, отличающийся тем, что уплотненный сигнал определяют как последовательность чисел, характеризующих глубины сигнала сейсмических данных, на которых находятся горизонты, а массив уплотненных сигналов определяют как трехмерный массив уплотненных сигналов, где каждый уплотненный сигнал соответствует одному из сигналов сейсмических данных указанного трехмерного массива сейсмических сигналов. 2. Способ по п. 1, отличающийся тем, что этап интерполяции включает в себя этап подбора параболы по трем точкам, определенным указанным атрибутом амплитуды импульса на идентифицированной глубине и амплитудами импульса на глубинах выше и ниже указанной идентифицированной глубины. 3. Способ по п.2, отличающийся тем, что атрибут амплитуды импульса является максимумом указанной амплитуды. 4. Способ по п.1, отличающийся тем, что определяют числовую величину атрибута указанного импульса на каждой из интерполированных глубин сжатого сигнала и для каждого из сжатых сигналов массива сжатых сигналов. 5. Способ по п.4, отличающийся тем, что сохраняют числовую величину атрибутов указанных импульсов совместно с соответствующими идентифицированными глубинами атрибутов сжатых сигналов. 6. Способ по п.1, отличающийся тем, что осуществляют выбор какой-либо исходной глубины на каком-либо сжатом сигнале, приходящимся на какой-либо конкретный горизонт, и автоматизированное сканирование конкретного горизонта по сжатому массиву сигналов, начиная с исходной глубины на сжатом сигнале, и выделения глубин на смежных сжатых сигналах указанного сжатого массива сигналов, которые лежат на конкретном горизонте. 7. Способ по п.6, отличающийся тем, что автоматизированное сканирование включает в себя: 1) идентификацию, начиная с показателя глубины на сжатом сигнале этого же показателя как указанной исходной глубины на смежных сжатых сигналах, 2) поиск по смежным сигналам для выявления глубин, находящихся в пределах заранее установленного интервала для показателей, в качестве глубины указанной исходной глубины, 3) обозначение любой найденной глубины на смежном сжатом сигнале, как находящейся на том же самом горизонте, что и исходная глубина, либо, в качестве альтернативы, повторение шага (2) с удвоенным интервалом для показателей, 4) повторение этапов 1), 2), 3) за исключением поиска по сжатым сигналам, на которых поиск уже был произведен, до тех пор, пока на смежных сжатых сигналах уже не будут больше обнаруживаться новые глубины. 8. Способ по п.1, отличающийся тем, что осуществляют сравнение выделенных точек глубины Z центрального сигнала каждого контрольного массива Y с точками глубины Z центрального сигнала соседних массивов Y так, чтобы каждая выделенная точка глубины Z центрального сигнала соседних контрольных массивов Y находилась в пределах заданной глубины точки глубины Z контрольного массива Y и чтобы глубина точки Z соседнего контрольного массива Y была аналогичной соответствующему побочному сигналу сейсмических данных точки глубины Z проверочного массива Y. 9. Способ по п.1, отличающийся тем, что определенными сигналами сейсмических данных являются все сигналы сейсмических данных, которые могут служить центральным сигналом утроенного трехмерного массива сигналов данных указанного трехмерного массива сигналов сейсмических данных. 10. Способ по п.1, отличающийся тем, что характеристика импульса контрольного сигнала сейсмических данных представляет собой максимальную амплитуду этого импульса. 11. Способ по п.1, отличающийся тем, что характеристика импульса контрольного сигнала сейсмических данных представляет собой минимальную амплитуду этого импульса. 12. Способ по п.1, отличающийся тем, что операция автоматического сканирования конкретного горизонта через битовый массив, начиная с исходного бита и выделения битов "1" в соседних битовых сигналах указанного битового массива, включает операции а) начиная с исходного бита в битовом сигнале, идентификации такой же глубины, как и у исходного бита, в соседних битовых сигналах, и b) поиска в соседних битовых сигналах для обнаружения битов, находящихся в пределах заданной глубины, как глубины этого исходного бита, с) обозначение любого обнаруженного бита соседнего битового сигнала, как находящегося на том же горизонте, что и исходный бит, и обозначение этого обнаруженного бита, как исходного бита, и d) повторение операций a), b) и c), за исключением поиска битовых сигналов, по которым уже проводился поиск, до тех пор, пока уже невозможно будет обнаружить битов или соседних битовых сигналов. 13. Способ по п.12, отличающийся тем, что такую форму бита, как функция глубины, наносят на систему координат X - Y с отображением изменения параметров глубины. 14. Способ по п.12, отличающийся тем, что осуществляют регистрацию каждого исходного бита как главного бита и регистрацию каждого обнаруженного бита соседнего битового сигнала как производного бита по отношению к этому главному биту. 15. Способ по п.1, отличающийся тем, что определяют числовую величину атрибута сигналов сейсмических данных на каждой глубине, на которой зарегистрирована единица ("1"), и осуществляют регистрацию этих числовых величин как функции глубины, на которой зарегистрированы эти единицы, а для каждого обнаруженного бита соседнего битового сигнала регистрацию этого обнаруженного бита как производного бита этого главного бита. 16. Способ по п.1, отличающийся тем, что определяют числовую величину атрибута сигналов сейсмических данных на каждой глубине, на которой зарегистрирована единица, и осуществляют регистрацию этих числовых величин как функции глубин, на которых зарегистрированы эти единицы. 17. Способ автоматического выбора горизонтов подповерхностных образований из трехмерного массива сигналов сейсмических данных с регистрацией связей предшественник - последователь между исходными и обнаруженными точками, отличающийся тем, что осуществляют: 1) выбор исходной точки, характеризующей сейсмический атрибут на сигнале сейсмических данных, который приходится на горизонт подповерхностных образований, при котором горизонт является поверхностью раздела пластов земной коры и представляется в сигналах сейсмических данных, проходящих через поверхность раздела пластов, в виде общего атрибута сейсмического импульса и при котором общий атрибут может изменяться по глубине, как функция координат X - Y сигналов сейсмических данных, 2) идентификацию той же глубины, начиная с исходной точки на сигнале сейсмических данных, что и глубина исходной точки на соседних сигналах сейсмических данных, 3) поиск соседних сигналов для обнаружения общего сейсмического атрибута в пределах заданной глубины как глубины исходной точки, 4) обозначение любого общего атрибута соседнего сигнала, как обнаруженной точки на том же горизонте, на котором находится исходная точка, и обозначение такой обнаруженной точки на соседнем сигнале как новой исходной точки, 5) регистрацию исходной точки как точки-предшественника и регистрацию каждой обнаруженной точки на соседнем сигнале как последователя указанной исходной точки-предшественника, 6) повторение операций 2), 3), 4) и 5), указанных выше, за исключением поиска сигналов, по которым уже проводился поиск, до тех пор, пока уже невозможно будет найти обнаруженные точки на соседних сигналах. 18. Способ по п.17, отличающийся тем, что осуществляют визуализацию на дисплее карты горизонта, отражающей обнаруженные точки, и визуализацию на дисплее траекторий от точек-предшественников к точкам-последователям из любой заданной точки на карте горизонта назад к исходной точке. 19. Компьютеризованный способ автоматического выделения горизонтов пластов из трехмерного массива сигналов сейсмических данных, при котором горизонт является поверхностью раздела пластов земной коры и представляется в сигналах сейсмических данных, проходящих через поверхность раздела пластов, общим атрибутом сейсмического импульса, и при котором общий атрибут может изменяться по глубине как функция координат X - Y - сигналов сейсмических данных, и при котором выделение горизонта начинают с исходной точки на глубине признака по меньшей мере одного сейсмического импульса и выбирают сигналы данных, которые имеют абрибут сейсмического импульса исходной точки, и продолжают эту операцию по всему трехмерному объему сигналов сейсмических данных за счет использования выделенных точек как новых исходных точек для выбора новых выделенных точек, отличающийся тем, что способ регистрации связи между исходными точками и выделенными точками предусматривает следующие операции: регистрацию каждой из исходных точек как точки-предшественника, регистрацию каждой выделенной точки, возникающей при выделении от исходной точки, как точки-последователя указанной исходной точки-предшественника и визуализацию на отображении выделенного горизонта траектории от выделенной точки-последователя назад к точке-предшественнику, которая привела к ее выбору в качестве выделенной точки. 20. Способ по п.19, отличающийся тем, что указанная траектория включает участки траектории, характеризующие траектории от каждой точки-предшественника к каждой точке-последователю.Описание изобретения к патенту
Настоящее изобретение относится к расшифровке сейсмических данных. Конкретно настоящее изобретение относится к машинному процессу выбора трехмерных сейсмических данных для формирования у специалистов в области нефтеразведки более четкого представления о геологическом строении и геометрии разреза. Более конкретно изобретение представляет собой наиболее точный и быстрый автоматизированный способ "выделения", или "отслеживания" отдельных сейсмических событий, или горизонтов, посредством трехмерных массивов данных. На фиг. 1 - 4 проиллюстрированы особенности используемые в способах выделения, соответствующих предшествующему уровню техники. На фиг. 1 приведено изображение части гипотического массива трехмерных сейсмических данных, которое дано для того, чтобы пояснить трехмерные зависимости, встречающиеся в тексте описания и сопуствующих чертежах настоящей заявки; на фиг. 2 - изометрическая проекция пяти сейсмических сигналов, которые иллюстрируют соотношение между некоторой "исходной точкой" и четырьмя соседними сигналами; на фиг. 3 - способ автоматизированного выделения по предшествующему уровню техники; на фиг. 4 - "итеративный" способ автоматизированного отслеживания по предшествующему уровню техники. На фиг. 1 приведено изображение части гипотического массива трехмерных сейсмических данных. Маленькие окружности в верхней части этого массива изображает наземное положение отдельных сигналов. Вертикальные линии изображают сейсмические сигналы, которые замерены по полному времени пробега вдоль оси этого массива. Такой пробег соответствует расстоянию или глубине, на которой генерируется импульс. Каждый отдельный сигнал является отображением акустических отражений от напластований в недрах земли в координатах "амплитуда-время". Фиг. 1 является наглядной иллюстрацией способа фактического представления каждого сигнала сейсмограммы. Каждый сигнал хранится в виде последовательности дискретных чисел, представляющих амплитудное значение сигнала по отношению к нулевой величине. Для адекватного представления любого числа, соответствующего амплитуде, необходимо использовать несколько "разрядов" (разряд - это двоичное число, принимающее два значения: 0 и 1). Как правило используется восемь, шестнадцать или тридцать два разряда. Такое представление числа разрядами повторяется для каждой точки времени, например, через 2- или4-миллисекундный интервал в течении шести секунд. Горизонтальное сечение, или временной срез представляет собой горизонатальный срез, или сечение, трехмерного массива данных. Оно показывает различные напластования в один и тот же момент времени. С другой стороны, карту горизонта или просто горизонт получают построение графика какого-либо атрибута конкретного импульса (обычно длительности импульса, иногда - максимальной или минимальной амплитуды) в координатах X-Y. Горизонт напоминает топографическую карту поверхности, однако является графическим изображение подповерхностных напластований. Атрибуты горизонта могут выделятся цветом, видом и т. п. Менее чем за десять лет компьютеризированная геологоразведка революционизировала сейсморазведку и разработку месторождений. Однако до недавнего времени один аспект расшировки сейсмических данных - выделение подповерхностных горизонтов (или просто выделение) оставался практически на уровне работы с карандашом и бумагой без использования автоматизированных компьютерных методов выделения. Традиционно выделение считают трудоемким процессом, поскольку производилось вручную, путем вычерчивания цветными карандашами одного за другим сейсмических размеров или профилей. В начале 80-х годов автоматизированные рабочие места САЕХ (аббревиатура названия "Разведка с помощью компьютера") предоставили специалистам в области сейсморазведки возможность выделять трехмерные данные с большей скоростью и эффективностью. Хотя расшифровка сейсмических профилей (т.е. двухмерных вертикальных сечений или "вертикальных сейсмических разрезов") все еще осуществлялась визуально - выделением одного профиля за другим, это стало возможным делать с использованием "мыши" на экране компьютера, устанавливая курсор на некоторые выбранные точки на каком-либо горизонте с последующим выделением компьютером всех остальных точек на этом профиле. Это явилось первым примером автоматизированного выделения, которое по сравнению с ручным позволяло увеличить не только производительность но и точность. В одной из предшествующих автоматизированных систем для выделения плоскости напластования (или формаций, или просто "горизонта") в существенно горизонтальном направлении через трехмерный массив данных пользователь выбирал (или "вводил") по крайне мере одну "исходную точку", которая затем "распространялась" по всем четырем направлениям в массиве трехмерных данных так, как это показано на фиг. 2, до тех пор пока не достигала границ определенной пользователем зоны. Пользователи могли выбирать между выделением сейсмических данных в одномодовом в двухмодовом режиме. "Исходная точка" определяется своими координатами X и Y, а также своим временем или глубиной (т.е. осью Z на фиг. 1). Она также определяется характеристикой или параметром отражения в этой точке. Такой характеристикой обычно является максимальная амплитуда отражения в этом месте массиве данных. Могут быть использованы и другие параметры или характеристики отражения в точке X, Y, Z, например, минимальная амплитуда, фаза, частота и т.п. Как показано на фиг. 3, при неинтеративном выделение проводился поиск сейсмических сигналов со сходными значениями амплитуды, прилегающих к исходным точкам, выбирался наилучший из них, а затем анализировался следующий имеющийся сигнал без перепроверки точности выборки. Итеративный режим выделения верефицировал соседний сигнал в качестве выделенного с использованием предыдущего сигнала. Верифицированный соседний сигнал рассматривался в качестве исходной точки, и выделение соседних сигналов продолжалось с него. Такой итеративный процесс, известный из предшествующего уровня техники, показан на фиг. 4. Верификация при проверке означает, что если величина амплитуды выделенного сигнала находится в пределах заданного пользователем допуска, то это выделение принимается. Пользователи могут устанавливать степень схожести (по шкале 1 - 10). Если выделенный сигнал не выдерживал такой проверки, он объявлялся "потерянным" до появления по крайней мере одного непосредственно соседствующего с ним сигнала достаточно совпадающего с ним. Более конкретно, как только на сигнале выбиралась исходная точка, сигнал сканировался вверх и вниз по оси Z, или оси времени, с целью нахождения локальных максимальных амплитуд, или просто "экстремумов". Локальный экстремум величины Xi, где i - показатель отцифровывания, определяется как:Xi-1<X>Xi+1 или
Xi-1>XXi<X+1
В случае максимума или минимума такое сканирование ограничивается точками перехода амплитуды сигнала через ноль. Подобный экстемум обычно немного варьирует с течением времени. Например, если t0 представляет исходную точку, то t1 обычно представляет время экстремума. Затем время t0 начинается на целевом сигнале. На целевом сигнале время варьирует вверх и вниз в диапазоне между точками перехода амплитуды сигнала через ноль до тех пор, пока не будет найден ближайший экстремум t2. Наконец, время t2 используют на сигнале, на котором существует исходная точка, и на таком"исходном" сигнале опять осуществляют сканирование вдоль оси Z для поиска ближайшего экстремума E3. Если t3 равно t1, то цель интеративного отслеживания считается достигнутой, и поиск продолжают дальше. Критерий отбора для интеративного отслеживания по предшествующему уровню развития техники определялся функцией

SCORE = (S

Эта сумма сравнивается с контрольной величиной, выбранной интерпретатором или конечным потребителем данных, и лежащей в диапазоне между 1 и 10. Если сумма превышает контрольную величину, целевой сигнал не выделяется. Для того, чтобы подготовить или выделить карту горизонта описанным выше способом по предшествующему уровню техники, требуется обработка весьма больших массивов данных. Необходимо не только выполнить процедуры выделения, но и сами эти процедуры требуют оперирования цифровыми данными, содержащими много разрядов информации,представляющей аналоговые сейсмические сигналы. В результате, даже при использовании на автоматизированных рабочих местах очень мощных компьютеров геологии и геофизики пользующиеся автоматизированными рабочими местами, снабженными программой выделения горизонта из трехмерных массивов, должны ждать, пока программа пройдет по всем данным и выполнит описанные выше процедуры. Такое ожидание может серьезно помешать рабочему процессу в тех случаях, когда пользователь желает просмотреть несколько горизонтов за короткое время. Основной целью настоящего изобретения является создание способа хранения и выделения трехмерных сейсмических данных, позволяющего по команде пользователя получать карту горизонта за значительно более короткое время, в 10 - 50 раз быстрее, чем способы по предшествующему уровню техники, использующие сопоставимую вычислительную мощность. Другая цель настоящего изобретения заключается в уменьшении не менее чем в 8 раз объема памяти вычислительной машины, необходимой для выделения горизонта. Другая цель настоящего изобретения состоит в создании альтернативного способа хранения и выделения трехмерных сейсмических данных, отличающегося высокой точностью генерации карт горизонтов и их атрибутов. Еще одной целью настоящего изобретения является создание такого способа выявления и интерактивной коррекции ошибок трассировки, который был бы намного быстрее и удобнее для пользователя по сравнению с известным. Указанные выше цели, а также другие преимущества и особенности настоящего изобретения обеспечиваются способом и устройством по настоящему изобретению, согласно которым основная масса расчетов, необходимых для процедур выделения их трехмерного массива сейсмических данных, выполняется заранее, на этапе "пакетной" (не интерактивной) обработки данных. Результаты пакетной обработки данных хранятся в виде " битового массива горизонта", в котором каждая сейсмическая величина исходного трехмерного куба оцифрованных данных заменена одним разрядом эталонной информации. Разрядный (битовый) массив горизонта получают, выбрав и обработав все сетки сейсмических профилей размером три на три и приписав разряд, равный 1, каждой точке на глубине, соответствующей центральному сигналу на сетке, для которой какой-либо показатель горизонта (или "характеристика") может быть выделен согласно методу выделения локального горизонта. такой характеристикой может быть максимальная амплитуда импульса, минимальная амплитуда импульса, либо другие характерные признаки (атрибуты), например переход через ноль, частота и т.д. Далее, каждый локальный горизонт, представленный разрядом "I" на центральном сигнале для конкретной точки на глубине, координируется с соседними сигналами так, чтобы образовалась непрерывная последовательность локального горизонта - от одного локального горизонта на глубине для центрального сигнала к соседнему локальному горизонту для центрального сигнала. В результате получают окончательный битовый массив, который при обращении к нему в интерактивном режиме выдает горизонты, начиная с любого бита (разряда) в массиве, выбранного в качестве исходной точки. Такие горизонты, получаемые в результате сканирования окончательного битового массива, являются эталонными эквивалентами горизонтов, выделенных на основе исходного восьмиразрядного массива сейсмических данных. Также предлагается способ и устройство хранения атрибутной информации, соответствующей каждой единице или "I" разряда в окончательном битовом массиве. Во время интерактивной фазы пользователь обращается только к окончательному битовому массиву, а не к данным исходного сейсмического массива. Этот массив данных может храниться в оперативной памяти компьютера, а не на диске, как того требовал предшествующий уровень техники. Сканированием окончательного битового массива из исходной точки получают карту горизонта. Процесс идет настолько быстро, что пользователь может воспроизвести карту горизонта как бы мгновенно. Ведется запись пути через битовый массив, так что в случае обнаружения ошибки в части результирующей карты, пользователь может выявить и удалить эту карту горизонта. Вышеописанный окончательный битовый массив обычно содержит 10 - 20 % единичных битов, а остальные биты являются нулями. По другому способу и устройству по настоящему изобретению вместо окончательного битового массива создают "сжатый массив сигналов". В том случае вместо записи единицы по месту эталонного ближайшего характерного признака импульса (например, пика, провала, либо какой-либо специфической фазы) хранится запись точного времени (либо его эквивалента - глубины), при котором этот характерный признак имел место. Эти значения времени или положения на глубине последовательно записываются в индексный файл. Этот альтернативный способ для определения точного времени или положения на глубине характерной величины для каждого локального горизонта предпочтительно интерполируют запись сейсмограмм. Характерные величины, определенные интерполяцией, записываются непрерывно, однако их положение на глубине определяется на основе величин, хранящихся в сжатом массиве сигналов. В интерактивном режиме расширение за исходную точку в сжатом массиве сигналов производится также, как и при сканировании окончательного битового массива, однако при этом используется другая процедура поиска. На фиг. 1 - 4 показаны автоматизированные способы выделения, соответствующие предшествующие уровню техники, и показывают трехмерные массивы сейсмических данных; на фиг. 5 - массив сигналов размером три на три, определенный вокруг центрального сигнала, используемого для получения возможных разрядов; на фиг. 6 а, б - один из способов определения того, находится или нет какая-либо характеристика импульса с конкретной глубиной на локальном горизонте (для того чтобы решить, должен ли быть включен бит в возможный битовый массив); на фиг. 7 а, б - процедуры проверки в ходе создания окончательного битового массива (для того чтобы гарантировать, что при последующем сканировании с любого бита в битовом массиве будет получен какой-либо горизонт); на фиг. 8 - процедуры проверки в ходе создания окончательного битового массива, при которых наборы битов для точек глубин центральных линий тестируются по отношению к соседним битам - как по отношению к восточным и западным, так и по отношению к северным и южным; на фиг. 9 - как после пакетной обработки пользователь выбирает исходную точку, причем исходная точка соответствует единичному биту в окончательном битовом массиве, а также иллюстрирует работу сканера по получению карты горизонта при одновременном сохранении исходной информации, пригодной для выявления и исправления ошибок на карте; на фиг. 10 - способ создания массива атрибутов из характерных величин, который соответствует единичным битам в окончательном битовом массиве. Лучший вариант осуществления изобретения. Способ по настоящему изобретению подразделяется на две фазы. Предложено программное обеспечение для вычислительной машины, предназначенное для работы на аппаратных средствах в процессе реализации данного способа. Первая фаза является пакетной фазой (не интерактивной) обработки данных, на которой трехмерные сейсмические данные сигналов первоначально обрабатываются на небольших пробных участках для того, чтобы выявить "возможные биты", репрезентативные для локальных горизонтов. Такие "возможные биты" затем обрабатывают по процедуре "сканирования" для того, чтобы гарантировать, что соседние локальные горизонты имеют совпадающие общие границы. Конечным результатом такой обработки является окончательный битовый массив, в котором каждый сигнал конвертирован в серию нулей и единиц как функций глубины и который представляет все горизонты. Единичные биты помещаются для глубин каждой плоскости напластования или "горизонта" по весовому массиву сейсмических данных. Горизонт определяется - как горизонтально непрерывная характеристика геологической формации, и обычно представляет собой максимум либо минимум импульса на дорожках сейсмограмм. Сейсмический импульс часто представляет собой отражения от глубинных геологических формаций. Помимо максимумов и минимумов можно использовать и другие атрибуты, или характеристики, однако в дальнейшем настоящая заявка будет ограничена только ими. Вторая фаза представляет собой интерактивную фазу, в которой окончательный битовый массив загружается в оперативную память автоматизированного рабочего места геологоразведчика. Пользователь сначала выбирает двухмерный сейсмический профиль, который выводится на мотор. Затем пользователь перемещает курсор на желаемые формации (т.е. максимум импульса) и нажимает клавишу. Такое нажатие определяет координаты X, Y, Z одной точки этих формаций. Эти координаты соотносятся с окончательным битовым массивом, хранимым в оперативной памяти, и бит, записанный под этими координатами X, Y, Z, определяется в качестве исходного для картирования горизонта. Эта исходная точка используется автоматическим способом выделения (названным "сканер"; сходную процедуру используют для получения окончательных битов из возможных битов) для распостранения этой исходной точки по всему трехмерному битовому массиву. Горизонт вырабатывается и выводится на дисплей с исключительной точностью и скоростью. Ниже следует детальное описание каждой фазы. Пакетная фаза. 1. Создание возможного битового массива. Фиг. 5 иллюстрирует преобразование цифровых данных, представляющих амплитуды сейсмических импульсов в виде функции глубины, в соответствующие разряды (биты), т.е. в 0 и 1, являющиеся функцией глубины, в которой единицы присвоены тем значениям глубин, на которых выделена заранее определенная характеристика импульса. Например, трехмерный сейсмический массив (в левой части фиг. 5) изображает массив сигналов размеров три на три, построенный вокруг центрального сигнала, являющимся функцией глубин Z и исходящим из точки с координатами X=1, Y=1. На глубинах Z1 и Z2 показаны два импульса 40, 42. Фактический сигнал может иметь много импульсов. Массив сейсмических сигналов размерностью три на три, примером которого является массив с центром в X=1, Y=1, обрабатывается как группа процессом PICK или программой 110, в результате чего получают центральный битовый сигнал в виде функции от глубины в точке X=1, Y=1. Предполагается, что выделение программой 110 подтвердило, что максимальные амплитуды, представленные на глубине Z1 импульсом 40 и на глубине Z2 импульсом 42, соответствуют отдельным локальным горизонтам сетке сейсмических профилей размером три на три. Возможный битовый массив в правой части фиг. 5 показывает, что центральный сигнал в качестве функции глубины имеет единичные биты, помещенные в Z1 и Z2, а все остальные глубины на этом сигнале представлены нулями. Такая процедура может быть повторена для всех центральных сигналов трехмерного массива сейсмических данных до тех пор, пока весь возможный битовый массив не будет полностью заполнен нулевыми и единичными битами (в качестве альтернативы, для того чтобы выявить окончательные биты по ходу процесса, могут быть найдены и тестированы по отношению к соседним локальным горизонтам возможные биты центрального сигнала). Программа PICK 110 (фиг. 5) применяется к каждой максимальной амплитуде каждого импульса каждой центрального сигнала каждого сейсмического битового массива размером три на три, как это описано ниже. 2. Определение возможных разрядов для каждого центрального сигнала возможного битового массива. Предпочтительный вариант способа PICK (фиг. 5) представлен на фиг. 6а, б. По этому способу вначале определяют максимальные амплитуды импульса для каждой глубины на центральном сигнале "массива" сигналов размером три на три. Затем выбирают набор пяти тестирующих импульсов к "югу", "западу", "северу" и "востоку" (обозначены А, B, С, D,) от центрального сигнала (фиг. 6а, б). Южный, западный, северный и восточный сигналы играют роль тестовых сигналов вокруг центрального сигнала О, который имеет глубину Z, соответствующую точке глубины максимальной (либо каким-то другим характерами. например минимальной) амплитуды импульса, как это показано на фиг. 6а. Предпочтительный вариант способа PICK 110 выделяют и итеративно верифицируют все четыре боковых сигнала по отношению к центральному сигналу, прежде чем точка максимума О на глубине Z1 соответствующего центрального сигнала начнет использоваться как единичный бит на глубине Z1 соответствующего сигнала возможного битового массива. Если какое-либо сравнение центрального сигнала и соседнего бокового сигнала не удовлетворит установленному пользователем критерию отбора, то максимум импульса на глубине Z1 не считается принадлежащим локальному горизонту, как не прошедший тестирование. Предпочтительный способ выделения по настоящему изобретению продемонстрирует поэтапный множественный процесс верификации сигналов с использованием набора по пять сигналов, с выделением четырех соседних (боковых) сигналов (фиг. 6б) и верификации каждой из четырех точек А, B, C и D точно также, как и итеративный режим по предшествующему уровню техники, показанный на фиг. 4. Каждая итеративная верификация использует так же самые определяемые показателем критерии, как это описано выше. Единичный бит записывается для точки глубины О в возможном битовом массиве только в том случае, если все четыре боковых сигнала из упомянутого набора пройдут тестирование. Если все четыре боковых сигнала пройдут такие тесты, то любой из боковых сигналов A, B, C и D может быть снова соединен с центральной точкой О. Такое правило "все или ничего" гарантирует, что точка глубины на центральном сигнале и выбранные глубины на его боковых сигналах действительно лежат на одном локальном сесмическом горизонте. В определенных условиях правило "все или ничего" не является оптимальным или необходимым и может оказаться достаточным использованием менее жестокого критерия. Выше описанный процесс повторяется для каждого импульса на центральном сигнале, расположенном в X=I, Y=I. Обычный сигнал имеет много импульсов с максимальной амплитудой. В общем случае имеется N импульсов, расположенных на глубинах Z1, Z2,...Zi...Zn, где i означает i-тый импульс. Для каждого последовательного теста, как это было описано выше, глубина соседнего сигнала записывается и хранится как точка на возможном локальном горизонте. Другими словами, единичный биты записываются и хранятся для "успешных" глубин Zi, и одновременно записывается глубина соседнего сигнала в A, B, C и D. Нули ("нулевые" бита) записываются во все эталонные точки для тех глубин (или интервалов времени), для которых не были записаны единичные биты. Подобная запись в памяти может быть представлена так, что это сделано в табл. 5. Независимо от этого, окончательным утверждением вышеописанный процесс выделения повторяют для каждого массива сигнала размером три или три вместо того, чтобы выполнять его для каждого сигнала. Например, массив данных размером три на три, центрированный в X=2, Y=2 обрабатывается в соответствии с вышеописанными способом выделения. Единичные биты записывают для глубины Z (назначены Z", (табл. 2). Как и в табл. 1, соответствующие глубины на боковых сигналах A, B, C и D также записываются для хранения. 3. Определение окончательного битового массива. По способу сканирования, используемого в "интерактивной фазе" (описан ниже), глубина любого бита может быть использована в качестве исходной точки для поиска бита на соседнем сигнале (т. е. сигнале, лежащем непосредственно к северу, югу, востоку или западу). Если такой поиск завершается успешно, то говорят, что исходный бит и вновь найденный бит "принадлежат" к одному и тому же горизонту. В предыдущем разделе был описан выбор возможных битов, который касался того, как локальный горизонт, заданный сейсмическими данными, используют для того, чтобы определить, какой разряд следует записать в память. Процесс, описанный в этом разделе, определяет, сможет ли каждый из бита генерировать локальный горизонт. Если да, то этот бит сохраняется, если нет, бит сбрасывается. Таким путем можно гарантировать, что "окончательный" битовый массив будет верно воспроизводить горизонты, определенные сейсмическими данными. С учетом указанного требования, и продолжая описанный выше пример (табл. 1 и 2), модно контактировать, что для гарантии нахождения единичных битов, предписанных глубин Z", на горизонте определенном способом выделения для X= 1, Y=1, массив данных 3х3, необходимо тестировать биты на различных глубинах Z" для X=2, Y=2, массив данных 3х3. Другими словами, для того чтобы определить, находится ли бит в Z" I в пределах заданного расстояния от бита в Zi (X= 1, Y=1), и наоборот, выполняют тестирование. Такое тестирование проводят в соответствии с табл. 3. Таблица 3. Тестирование выбора Z" (например, для массива данных с X=2, Y=1, 3х3 по отношению к массиву данных с X=1, Y=1, 3х3)
1. поочередно оценивают каждый бит центрального сигнала в







Класс G01V1/28 обработка сейсмических данных, например их анализ для интерпретации, коррекции
Класс G01V1/32 преобразование одной записи в другую
способы и системы сжатия данных - патент 2371741 (27.10.2009) | ![]() |