Способ интерпретации сейсмических фаций с использованием текстурного анализа и нейронных сетей

1 Область техники Настоящее изобретение относится к геофизическим методам разведки, а более точно - к способу описания и картирования сейсмических фаций в сейсмических данных. Предшествующий уровень техники Анализ сейсмических фаций является важным этапом для интерпретации сейсмических данных при описании пласта. Интерпретации сейсмических фаций имеют существенное значение при начальной разведке коллектора,оценке результатов разведки, описании коллектора и, в конечном счете, при разработке месторождения. Сейсмические фации представляют собой стратиграфическую единицу или область,которая имеет характерную диаграмму отражения, отличающуюся от диаграмм отражения других областей. Области различающихся сейсмических фаций обычно описывают, используя описательные термины, которые отражают крупномасштабные сейсмические характеристики, такие как амплитуда отраженной волны, непрерывность и внутренняя конфигурация отражающих поверхностей, ограниченных стратиграфическими горизонтами. Применение и анализ сейсмических фаций существенно изменяются от случая применения ко всему бассейну до случая подробного описания пласта. В масштабе всего бассейна разведочный анализ сейсмических фаций применяют при исследовании углеводородных систем с целью широкой идентификации участков источника, пласта-коллектора и изолированных наклонных областей. Эти области обычно идентифицируют на основе геометрии отраженных волн, а также амплитуд сейсмических волн и непрерывности. Локально полунепрерывные отражающие поверхности с высокими амплитудами сейсмических волн часто используют для идентификации потенциальных коллекторов углеводородов, таких как глубоководные каналы, тогда как непрерывные или полунепрерывные области с низкими амплитудами сейсмических волн можно использовать для идентификации изолированных наклонных единиц. Анализ сейсмических фаций можно также применять в пределах одного пластаколлектора, чтобы получать подробное описание физических свойств. В применениях локального масштаба характеристики непрерывности и амплитуды обычно строго не определены и основаны на измерении свойств породы или на интерпретации отложений. Если можно выявить зависимость между особенностями сейсмической волны и физическими свойствами, объемы сейсмических фаций можно использовать для предсказания распределения свойств породы и геологических моделей. Обычная методика, используемая для анализа сейсмических фаций и картирования, представляет собой ручной процесс, где интерпретатор данных сейсмической разведки на основе 2 визуального наблюдения принимает решение о характере данных отражения сейсмических волн в пределах представляющего интерес интервала и наносит данные на карту. Затем сейсмические фации используют для различных целей, но, прежде всего, для интерпретации распределения литофаций или свойств породы. Восприятие, интуиция и опыт квалифицированного интерпретатора оказывают существенное влияние на успех анализа сейсмических фаций. Однако они могут стать причиной того, что проведение анализа сейсмических фаций станет субъективной, требующей много времени и трудной задачей. В нефтедобывающей промышленности используют несколько способов автоматизации и улучшения интерпретации сейсмических фаций из сейсмических данных.Leading Edge, Geophys. Explor., vol. 5,9, pp. 51-55, 1986, заложили основу для формирования алгоритма и, следовательно, для автоматизации процесса интерпретации сейсмических данных. Однако эти авторы не раскрыли действующую модель или описание каких-либо особенностей возможных атрибутов или алгоритмов классификации.seismic texture classification", Society of Petroleum Engineers 35482, 1996, рассмотрели текстурный анализ сейсмических данных, а также классификацию текстурных атрибутов с использованием версии анализа главных компонентов и распределений вероятности. Хотя в этих публикациях способы текстурного анализа применяются для сейсмических данных, но не реализованы преимущества вероятностных нейронных сетей и не используются вероятностные значения для оптимизации классификации. В этих способах также не используется схема интерактивного обучения, а текстурный анализ не управляется углом наклона. Процесс направления вычисления при стратиграфическом расслоении, определяемый углом наклона сейсмических отражающих поверхностей, называют управлением углом наклона.seismic textures: Implication for quantitative seismic interpretation and hydrocarbon exploration",1999, описал применение обычного текстурного анализа для получения сейсмических текстурных атрибутов, которые дают количественную оценку интенсивности отражения, непрерывности и геометрии. Однако в этой публикации не описаны способы классификации текстурных атрибутов. В частности, Gao, 1999, не использу 3 ет ни вероятностную нейронную сеть, ни интерактивный интерпретатор, обучающий нейронную сеть. Кроме того, текстурный анализ не является управляемым углом наклона.of Seismic Attributes разработал способ, в котором различные сейсмические атрибуты применяют для интерактивного обучения нейронной сети. Однако текстурные атрибуты не используют, а разработанная сеть является полносвязанной нейронной сетью с обратным распространением, а не вероятностной нейронной сетью. Р. Meldahl, R. Heggland, P.F.M. de GrootInterpretation", а также в публикации международной заявки WO 00/16125, "Method of seismicsignal processing", используются сейсмические атрибуты для интерактивного обучения нейронной сети и получения объема фаций. Однако при обучении и получении куба данных столба одновременно рассматривается и классифицируется только один класс группы данных, а не большое количество классов. Поэтому в архитектуре нейронной сети используются только два конечных выходных узла. Вычисляется вероятностный куб данных, а затем в качестве стадии постобработки формируются релейные пороговые значения, предназначенные для принятия решения о том, относится или нет объект к классу, представляющему интерес. Для управления сейсмическими атрибутами используется сложная схема преобразования Винера-Радона. Атрибуты для отдельных классов выбираются вручную.recognition", a также в патентах FR 2738920 19970321 и ЕР 808467 19971126, раскрыт способ распознавания сейсмических образов, основанный на использовании сейсмических трасс. Каждую сейсмическую трассу в пределах заданного пользователем интервала разлагают на заданное пользователем число эмпирических ортогональных функций. Затем полученные функции классифицируют, используя алгоритм классификации, реализуемый нейронной сетью, а не текстурный анализ с помощью обученного интерпретатора. Краткое изложение существа изобретения Задачей настоящего изобретения является разработка способа, который обеспечивает возможность быстрой объективной классификации сейсмических данных и который можно использовать в процессе картирования сейсмических фаций. В этом способе также должен имитироваться процесс, который используется в сейс 004551 4 мическом интерпретаторе, и должны имитироваться результаты, получаемые вручную с помощью сейсмического интерпретатора. Согласно настоящему изобретению предложен способ идентификации сейсмических фаций в объеме сейсмических данных. Сначала вычисляют большое количество начальных текстурных атрибутов, характерных для объема сейсмических данных. Затем конструируют вероятностную нейронную сеть на основании вычисленных начальных текстурных атрибутов. Далее вычисляют конечные текстурные атрибуты для всего объема сейсмических данных. В заключение, вычисленные текстурные атрибуты классифицируют, используя сконструированную вероятностную нейронную сеть. Краткое описание чертежей Настоящее изобретение и его преимущества поясняются нижеследующим подробным описанием со ссылками на сопровождающие чертежи, где фиг. 1 изображает схему алгоритма способа согласно изобретению; фиг. 2 - пример сейсмического поперечного разреза, полученного при исследовании участка, согласно изобретению; фиг. 3 А - сейсмический поперечный разрез, на котором показаны многоугольники 3 В,3 С, 3D, 3 Е, 3F, 3G, выбранные для текстурного анализа, согласно изобретению; фиг. 3 В-3G - классификации фаций, соответствующих многоугольникам на фиг. 3 А, согласно изобретению; фиг. 4 - разрез классификаций фаций, вычисленный с помощью заявленного способа,согласно изобретению; фиг. 5 - разрез уклона, используемый при управлении углом наклона, согласно изобретению; фиг. 6 - достоверный разрез, соответствующий разрезу классификаций фаций, показанному на фиг. 4,фиг. 7 - объем сейсмических фаций согласно изобретению. Описание предпочтительных вариантов осуществления изобретения Способ распознавания и картирования сейсмических фаций в сейсмических данных,преимущественно в сейсмических амплитудных данных, осуществляется согласно алгоритму,приведенному на фиг. 1. На шаге 101 выбирают трехмерный объем сейсмических данных. Хотя рассматривается трехмерный объем данных, способ работает одинаково хорошо и в случае множеств двумерных данных. Этот объем сейсмических данных будет использован для вычисления объема сейсмических фаций и соответствующих доверительных объемов. Предпочтительно, чтобы сейсмические данные были сейсмическим атрибутом или амплитудными данными, включаю 5 щими данные ближнего, дальнего и полного суммирования. На шаге 102 из объема сейсмических данных, выбранных на шаге 101, выбирают, по меньшей мере, один поперечный разрез. На шаге 103 на поперечных разрезах строят большое количество многоугольников. Нет необходимости в том, чтобы многоугольники были одного и того же размера, ориентации или из одного и того же стратиграфического интервала, они могут быть распределены на всем протяжении поперечных разрезов и объема при любой подходящей ориентации или геометрии. Предпочтительно, чтобы многоугольники строились путем представления их в цифровом виде на индикаторе выбранных поперечных разрезов. На шаге 103 многоугольники строят для того, чтобы на основании объема сейсмических данных получить изображение, характеризующее тип фаций. Предпочтительно предусматривать достаточное количество примеров каждого представляющего интерес типа фаций, чтобы описать вариацию, имеющуюся во входном объеме сейсмических данных, выбранных на шаге 101. Типы фаций описываются с помощью сейсмической текстуры. Сейсмическая текстура является характеристикой, которая дает количественную оценку многим аспектам описания сейсмических фаций, выполняемого сейсмическим интерпретатором. Сейсмическая текстура представляет собой количественную меру амплитуды отражения, непрерывности и внутренней конфигурации отражателей. Сейсмические текстурыможно охарактеризовать как гладкие или шероховатые, мелкомасштабные или крупномасштабные, описать как текстурные атрибуты и оценить количественно с помощью стандартных статистических методов. Сейсмическая текстура представляет собой по своей природе многотрассовый сейсмический атрибут и, тем самым, существенно отличается от многих вычисленных традиционным путем сейсмических атрибутов, основанных на одиночных трассах. Поэтому анализ сейсмической текстуры напоминает процесс основанного на визуальном представлении анализа, выполняемого сейсмическим интерпретатором, в том смысле, что традиционный анализ атрибута не делается. Интерпретатор не исследует одну или две трассы одновременно. Точнее, чтобы выполнить классификацию, интерпретатор исследует совокупность трасс в виде изображения. Это иной метод анализа, он предоставляет возможность сбора данных о геометрии отражения в пределах всей исследуемой области. В текстурном анализе пространственная организация значений пикселов внутри заданной области представляется указанными выше многоугольниками на поперечных разрезах. В общем случае эту область, окно текстурного анализа, называют "текселом". Для количественной оценки текстуры изображения в текселе 6 использовано преобразование изображения,которое приводит к получению яркостных матриц совпадения событий. Яркостные матрицы совпадения событий описывают пространственные зависимости между пикселами небольшой области в пределах большего изображения в текселе. На практике яркостные матрицы совпадения событий вычисляют в перекрывающихся текселах, так что можно полностью наблюдать любой переход между текстурными классами в пределах всего изображения. Перекрывающиеся текселы сканируют по горизонтали и вниз в пределах изображения до тех пор, пока не будет обработано все изображение. Яркостные матрицы совпадения событий представляют собой матрицы размером NxN,где N - число уровней яркости, используемых для количественной оценки изображения. Например, восьмиразрядные данные имеют 28=256 уровней яркости, а яркостная матрица совпадения событий, образованная из этого изображения, будет матрицей, которая имеет 256 строк и 256 столбцов. Вычисление и анализ яркостных матриц совпадения событий представляет собой задачу, требующую больших вычислительных ресурсов, при этом вычислительные потребности пропорциональны N2. Каждый элемент в яркостной матрице совпадения событий выражает относительную частоту появления двух точек с соответствующими значениями i и j пикселов на расстоянии D(d, ) друг от друга в пределах тексела. Например, если пиксел А имеет значение i и находится на расстоянии D от пиксела В со значением j, то позиции i, j яркостной матрицы совпадения событий необходимо увеличить на единицу. Этот процесс осуществляют для каждого существующего множества пикселов в пределах тексела. В самом общем случае применения яркостной матрицы совпадения событий вычисления могут отражать как переход значений пикселов, так и направление или "зернистость" в изображении. Текстурный анализ путем конструирования яркостной матрицы совпадения событий на основании тексела изображения является эффективным двумерным (или трехмерным) расширением одномерного последовательного анализа Маркова. Структуру яркостных матриц совпадения событий, полученных из сейсмических данных,можно понять эвристически. В однородных областях, в которых однородность или непрерывность точно определена в заданном направлении, различия между значениями пикселов будут меньше, и поэтому элементы, близкие к диагоналям яркостных матриц совпадения событий, будут иметь более высокие значения. В случае менее однородных областей различия между значениями соседних пикселов будут более высокими и поэтому получающиеся яркостные матрицы совпадения событий будут иметь более высокие значения на дополнитель 7 ном расстоянии от диагонали. Среднее значение пиксела также представлено в яркостной матрице совпадения событий. Области низкой амплитуды имеют яркостные матрицы совпадения событий со значениями, сгруппированными вблизи центра. С другой стороны, области с более высокой амплитудой имеют более распределенные значения яркостной матрицы совпадения событий либо вдоль диагонали в случае непрерывных текстур, либо на всем протяжении яркостной матрицы совпадения событий при более разрывных текстурах. На шаге 104 осуществляют классификацию начальных фаций для многоугольников,выбранных на шаге 103. Примеры классификаций фаций, обычно используемых в настоящем изобретении, включают непрерывную по высоким амплитудам, полунепрерывную по высоким амплитудам, непрерывную по средним амплитудам, полунепрерывную по средним амплитудам, непрерывную по низким амплитудам, полунепрерывную по низким амплитудам, хаотическую и по прозрачности. Примеры первых шести классификаций показаны соответственно на фиг. 3 В, 3 С, 3G, 3D, 3 Е и 3F. На шаге 105 яркостные матрицы совпадения событий конструируют на основании изображений в образованных на шаге 103 многоугольниках для фаций из шага 104. В этом случае статистические преобразования этих матриц описывают пространственные зависимости между пикселами небольшой области. На шаге 106 начальные текстурные атрибуты вычисляют на основании сконструированных на шаге 105 яркостных матриц совпадения событий, используя скользящее окно, задаваемое пользователем. Это вычисление, в целом, может быть названо формированием значений сейсмических текстур в сейсмических данных. Для полной имитации процесса, осуществляемого сейсмическим интерпретатором, предпочтительно вычислять двумерные текстурные атрибуты и затем фильтровать на временном интервале, чтобы полностью имитировать трехмерную операцию. В качестве варианта можно также вычислять трехмерные текстурные атрибуты и использовать для описания сейсмических фаций. Яркостные матрицы совпадения событий не являются эффективно интерпретируемыми непосредственно, а более эффективно описываются скалярными статистическими характеристиками, называемыми текстурными атрибутами. Текстурные атрибуты можно подразделить на дескрипторы первого и второго порядков. Статистические данные первого порядка дают количественное представление о глобальном распределении значений пикселов в пределах изображения и могут быть вычислены непосредственно на основании тексела с использованием стандартных статистических методов даже без промежуточного преобразования к яркостной матрице совпадения событий. Сред 004551 8 няя абсолютная амплитуда и среднеквадратическое отклонение значений амплитуды в пределах тексела являются примерами текстурных атрибутов первого порядка, полезными при описании амплитудных аномалий и интенсивности отражения. Производные атрибуты, например мгновенную амплитуду, фазу и частоту,также можно использовать для получения статистических данных первого порядка. Статистические данные первого порядка отражают начальный подход к детальной количественной оценке текстуры, и, хотя некоторые геофизические области можно грубо описать на основании различимых интервалов значений пикселов, в общем случае отдельные текселы не могут быть адекватно описаны на основе только их статистических данных первого порядка. Например, при использовании только средних значений амплитуд неупорядоченная область высоких амплитуд сейсмического изображения не может быть непременно отделена от непрерывной области высоких или даже умеренных амплитуд. Статистические данные второго порядка изображения дают количественную оценку пространственной взаимосвязи пикселов в пределах изображения и вычисляются путем промежуточного преобразования к яркостной матрице совпадения событий. Статистические данные второго порядка, статистические данные яркостной матрицы совпадения событий, включают в себя характеристики формы трассы, геометрию отражения и непрерывность отражения в дополнение к величине амплитуд. Статистические данные второго порядка текселов представляют собой многотрассовый атрибут изображения, который обеспечивает возможность получения геометрии отражения и непрерывности путем анализа управляемой углом наклона яркостной матрицы совпадения событий. Текстурные атрибуты, которые предпочтительно использовать в настоящем изобретении,включают в себя, но не ограничены ими, текстурную неоднородность, инерцию (также известную как момент различия элементов или контраст), энтропию и энергию (также известную как равномерность). Математические выражения для этих текстурных атрибутов имеют вид: где cij - i-ая и j-ая составляющая яркостной матрицы с совпадения событий, 9n - размер матрицы (квадрат целого числа уровней яркости в пределах изображения). Кроме того, для исключения элементов матрицы, которые больше 1, и, тем самым, статистического выравнивания составляющих любой единичной матрицы при вычислении конкретного атрибута входная яркостная матрица с совпадения событий нормируется так, что Первый текстурный атрибут - текстурная однородность - характеризует количество локальных подобий внутри тексела. Вследствие обратной пропорциональности (i-j)2 локальная текстурная однородность будет больше для яркостных матриц совпадения событий с элементами, сконцентрированными вблизи диагонали. Эти яркостные матрицы совпадения событий соответствуют текстурам с плохими характеристиками организации и контраста, но с небольшим числом уровней яркости на одинаковых расстояниях друг от друга и при одном и том же разнесении по азимуту. Меньшие значения текстурной однородности будут соответствовать большим значениям яркостной матрицы совпадения событий, дальше отстоящим от диагонали матрицы, то есть многочисленным различающимся уровням яркости, которые находятся на тех же расстояниях и азимутах. Эти характеристики делают текстурную однородность особенно полезной для количественной оценки непрерывности. Второй текстурный атрибут - текстурная инерция - является показателем контраста яркостной матрицы совпадения событий и характеристикой, противоположной текстурной однородности. В случае высококонтрастного изображения текстурная однородность низкая, текстурная инерция высокая. Третий текстурный атрибут - текстурная энтропия - характеризует потерю пространственной организации в пределах вычислительного окна. Текстурная энтропия является высокой,когда все элементы яркостной матрицы совпадения событий равны, что соответствует грубой текстуре, и является низкой, когда текстура более однородная или сглаженная. Четвертый текстурный атрибут - текстурная энергия - также характеризует пространственную организацию в пределах вычислительного окна. В противоположность текстурной энтропии, текстурная энергия является самой низкой, когда все элементы яркостной матрицы совпадения событий равны. В этом случае равновероятны все или большая часть яркостных уровней в пределах вычислительного окна. Это является признаком грубой текстуры. И, наоборот, самые высокие значения текстурной энергии показывают наличие высоких значений в яркостной матрице совпадения событий. В этом случае только некоторые из яркостных уровней являются доминирующими. Область внутри 10 вычислительного окна является более однородной или имеет до некоторой степени регулярный характер. На шаге 107 на основании начальных текстурных атрибутов и соответствующих классификаций фаций, на основании результатов шагов 105 и 106 соответственно, конструируют вероятностную нейронную сеть. Нейронная сеть представляет собой взаимосвязанный узел простых обрабатывающих элементов. Способность нейронной сети к обработке заключается в силе соединений или в весах, приобретаемых при осуществлении процесса адаптации к или обучения с помощью множества обучающих примеров. Одно из преимуществ нейронных сетей заключается в возможности обучения или модификации силы соединений внутри сети для получения желаемых результатов. В случае применения для классификации нейронная сеть может рассматриваться как особый случай классификации с контролируемой схемой, когда обучение нейронной сети является выполнением контролируемых упражнений. После достаточного обучения на ряде калибровочных изображений нейронную сеть можно применять для остальных изображений в объеме данных. Во время вычислительного процесса благодаря связанности узлов внутри обычной нейронной сети, то есть весам, осуществляются модификации входных векторов атрибутов, а модифицированные значения передаются на следующий слой сети. Посредством обучения веса сети модифицируются так, что на основе конкретного множества обучающих примеров модификация входных векторов атрибутов создает желаемый результат. Обучение сети и модификация весов соединений приводят к получению поверхности решений для сети. Поверхность решения является n-мерной поверхностью, которая придает сети возможность разделения входных обучающих данных на категории. Одно из преимуществ алгоритма нейронной сети по сравнению со схемами обычной классификации заключается в возможности получения нелинейных границ. При решении типичных задач классификации или предсказания обычно имеются только три слоя: первый входной слой; второй - "скрытый" слой; и третий - выходной слой. Вероятностные нейронные сети представляют собой параллельные реализации классификатора Байеса. Вероятностная нейронная сеть является трехслойной сетью, которая может эффективно выполнять классификацию образов. Математически вероятностные нейронные сети весьма сходны с кригингом, где близость к известным точкам используется для управления классификацией и предсказанием неизвестных точек. В обычном виде вероятностная нейронная сеть не обучается тем же способом, что и известная нейронная сеть, описанная выше. Точнее, обучающие векторы просто становятся 11 весовыми векторами в первом слое сети. Этот более простой способ придает вероятностной нейронной сети преимущество, заключающееся в том, что не требуется экстенсивное обучение. Например, при сейсмическом текстурном анализе текстурные атрибуты обучающих изображений представляют весовые векторы в первый слой сети. В результате, это дает преимущество в виде заметного ускорения обучения по сравнению с архитектурами нейронных сетей более традиционных типов, такими как полносвязанные архитектуры с обратным распространением. Кроме того, для вероятностной нейронной сети характерно обобщение, тогда как для известных сетей, даже при большом количестве обучающих данных, сходимость не гарантируется, а на стадии обучения обобщение данных не используется. Когда вероятностной нейронной сети представлен входной пример, в первом или во входном слое вычисляются отклонения входного вектора от обучающих входных векторов и формируется вектор, элементы которого указывают, насколько входные данные близки к обучающим данным. Во втором слое эти вклады суммируются для каждого класса входных данных для формирования в качестве выходного сигнала сети вектора вероятности. Это обеспечивает еще одно преимущество, вытекающее из использования вероятностной нейронной сети. Оно заключается в возможности извлечения вероятностной классификации непосредственно из второго или скрытого слоя в дополнение к классификации с максимальной вероятностью из третьего или выходного слоя. В настоящем изобретении входными обучающими точками для вероятностной нейронной сети, сконструированной на шаге 107, являются начальные текстурные атрибуты из шага 105 и соответствующие классификации начальных фаций из шага 106. Выходными данными вероятностной нейронной сети будут классификации фаций (и вероятностный объем, рассмотренный ниже). Затем вероятностную нейронную сеть можно использовать для классификации всего объема сейсмических данных. Однако в этом месте предпочтительно осуществлять контроль качества и при необходимости полностью повторно обучать вероятностную нейронную сеть. Поэтому на шаге 108 исходную вероятностную нейронную сеть используют для классификации фаций в части объема сейсмических данных из шага 101. Предпочтительно, чтобы эта часть была одним из поперечных разрезов,выбранных на шаге 102. На шаге 109 определяют, является ли удовлетворительной классификация фаций из части объема сейсмических данных. Если находят, что классификация фаций не является удовлетворительной, то процесс возвращают к шагу 103. Обучающее множество может быть модифицировано либо путем уда 004551 12 ления существующих многоугольников, либо путем добавления новых многоугольников. Затем воссоздают вероятностную нейронную сеть вместе с модифицированным обучающим множеством и снова осуществляют контроль. В настоящем изобретении эта возможность обучения и быстрого контроля вероятностной нейронной сети и последующей интерактивной модификации ранее существовавшего обучающего множества обеспечивает возможность повторного проведения классификации фаций, которая в случае интерпретатора выполняется вручную. Только после этого процесс будет продолжен до классификации всего объема сейсмических данных. Поэтому если на шаге 109 определяют, что частичная классификация фаций является удовлетворительной, то процесс продолжают. На шаге 110 вычисляют классификации конечных фаций во всем объеме сейсмических данных из шага 101, используя вероятностную нейронную сеть, сконструированную на шаге 107. Этим образуют классификационный объем сейсмических фаций, основанный на сейсмических текстурных атрибутах, полученных из первоначальных, заданных пользователем многоугольников. Качество объема сейсмических фаций зависит от качества входных данных. Снижение качества входных данных часто наблюдается при увеличении глубины нижнего горизонта. Использование окна единственного размера при вычислении яркостной матрицы совпадения событий для всего объема вносит вклад в этот отрицательный эффект. Результаты улучшаются при изменении размера окна на протяжении объема. Предпочтительно делать ширину окна большей по мере снижения частоты данных с увеличением глубины. Этот способ работает в сочетании с динамической регулировкой размера окна на основании задаваемого пользователем доверительного уровня. В другом альтернативном варианте осуществления изобретения,касающемся снижения качества сейсмических данных, до подачи на вход данные могут быть сначала подвергнуты фильтрации с помощью свертывающего или медианного фильтра для сглаживания данных. В заключение, на этапе 111 из выходных данных вероятностной нейронной сети также получают доверительный объем. В альтернативном варианте осуществления доверительный объем может динамически использоваться во время вычисления классификационного объема сейсмических фаций на шаге 110. Если степень достоверности падает ниже заданного уровня,процесс вычисления размера окна может автоматически корректироваться до тех пор, пока доверительный интервал не превысит приемлемый уровень, при этом фации вычисляют повторно и повторно классифицируют. Еще в одном альтернативном варианте осуществления процесс получения яркостных 13 матриц совпадения событий на шаге 105 может регулироваться в зависимости от угла падения. Стратиграфическое строение конкретного геологического отложения является важным аспектом, который всегда учитывается сейсмическим интерпретатором, хотя и неосознанно. Например, интерпретаторы сейсмических фаций не рассматривают непрерывность исключительно во временной плоскости. Скорее, они оценивают непрерывность, руководствуясь стратиграфическим расслоением, обозначенным углом наклона сейсмических отражающих поверхностей. Как описано выше, текстурный анализ и конструирование яркостной матрицы совпадения событий для тексела зависят от направления просмотра или азимута , на котором пикселы в пределах тексела являются связанными. Текстурный анализ в приложении к сейсмическим данным является чрезвычайно чувствительным к стратиграфическому строению тексела и, чтобы соответствующим образом имитировать процесс, осуществляемый человеком-интерпретатором, также должен подчиняться стратиграфическому углу наклона отражающих поверхностей. Слежение за стратиграфическим углом наклона при вычислении яркостной матрицы совпадения событий максимизирует непрерывность изображения, представленного в яркостной матрице совпадения событий. Процесс управления вычислением с помощью стратиграфического угла наклона называют управлением углом наклона. Чтобы должным образом управлять вычислением яркостной матрицы совпадения событий,для текстурного анализа необходима высокая степень разрешения стратиграфической геометрии. Для получения требуемого разрешения используют многотрассовый изобразительный характер тексела, а углы наклона в пределах изображения оценивают с помощью градиентного метода. В этом методе на первом шаге требуется вычисление горизонтального (dx) и вертикального (dy) градиента значений пикселов в пределах изображения. Затем вычисляют локальный уголнаклона отражающих поверхностей Отношение dy/dx имеет размерность времени, приходящегося на общую глубинную точку. Однако для удобства эту размерность можно проигнорировать, и угол наклона может быть выражен в фиктивных градусах, отнесенных к временному промежутку по горизонтали. В заключение, в качестве необязательного шага, для удаления шума может быть использован выбранный пользователем медианный фильтр. После того как угол наклона отражения станет известным повсюду в пределах тексела,при вычислении яркостной матрицы совпадения событий угол наклона используют для указания азимута просмотра для каждого пиксела для 14 сравнения пикселов. Участки резкого изменения угла наклона плохо отображаются при использовании не управляемого углом наклона способа вычисления яркостной матрицы совпадения событий. Отрицательные эффекты резкого изменения угла наклона можно минимизировать путем сглаживания или датирования объема вдоль стратиграфического слоя до осуществления анализа фаций. Для способа согласно настоящему изобретению не нужно использовать скважинные данные в качестве калибровочных данных. Это является преимуществом при разведке и на уже разработанных аренах, где имеется немного скважин, доступных для сейсмической калибровки на скважине. Калибровку всегда можно выполнить, в основном, ее осуществляют после вычисления. Однако она не нужна при использовании способа согласно настоящему изобретению. Для других способов с использованием сейсмических атрибутов и нейронных сетей необходима корреляция между сейсмическими и скважинными данными. Настоящее изобретение использовалось для образования объемов сейсмических фаций из обычных сейсмических амплитудных данных. Оно также использовалось применительно к объемным атрибутивным данным зависимости амплитуды отражения от удаления, например к объемам с пересечением склонов. Хотя для классификации фаций вычисляют и используют многочисленные текстурные атрибуты, в настоящем изобретении все необходимые атрибуты вычисляют по мере необходимости (на лету). Поэтому классифицируют только сейсмический объем, а фации и объемы вероятности запоминают в любое выделенное время. Текстурные атрибуты и другие объемы не создают. Это обеспечивает преимущество,заключающееся в том, что не требуется большое пространство памяти для данных. Настоящее изобретение пригодно для картирования сейсмических фаций на одной линии или на всем протяжении трехмерного объема. Возможность преобразования обычной амплитуды сейсмической волны или атрибутивных объемов в объемы сейсмических фаций приводит к существенному уменьшению затрат времени, к повышенной точности и к воспроизводимости в рамках процесса интерпретации сейсмических данных. Объемы сейсмических фаций используют для общего анализа геометрии и непрерывности пласта-коллектора, для выбора схемы расположения скважины и для улучшения геологических моделей, предназначенных для использования при планировании и организации работ по разработке пласта. Результаты исследования в заданном районе показали эффективность настоящего изобретения. На фиг. 2 показан сейсмический поперечный разрез, выбранный на шаге 102. На фиг. 3 А показаны многоугольники, выбранные 15 для текстурного анализа, используемые на шаге 103. Соответствующие классификации фаций показаны на фиг. 3 В-3G, например, используемые на шаге 106. На фиг. 4 показан полученный разрез классификаций фаций, например, вычисляемый на шаге 109. На фиг. 5 показан разрез уклона, используемого для управления углом падения. На фиг. 6 показан соответствующий доверительный разрез, вычисленный на шаге 111. Низкие значения достоверности можно наблюдать в зонах разлома и вблизи них, где стратиграфические и структурные взаимодействия усложняют интерпретацию фаций. Наконец, на фиг. 7 показан объем сейсмических фаций, вычисленный на этапе 110. ФОРМУЛА ИЗОБРЕТЕНИЯ 1. Способ идентификации сейсмических фаций в объеме сейсмических данных, заключающийся в том, что(a) вычисляют множество начальных текстурных атрибутов, характеризующих объем сейсмических данных,(b) конструируют вероятностную нейронную сеть на основании вычисленных начальных текстурных атрибутов,(c) вычисляют классификацию фаций в части объема сейсмических данных,(d) повторяют шаги с (а) по (с) до тех пор,пока вычисленные классификации фаций в части объема сейсмических данных не станут удовлетворительными, и(е) вычисляют классификацию фаций по всему объему сейсмических данных при использовании сконструированной вероятностной нейронной сети. 2. Способ по п.1, отличающийся тем, что при вычислении начальных текстурных атрибутов осуществляют выбор, по меньшей мере, одного поперечного разреза из объема сейсмических данных,конструирование множества многоугольников на выбранных поперечных разрезах и вычисление начальных текстурных атрибутов по изображениям в сконструированных многоугольниках. 3. Способ по п.2, отличающийся тем, что для вычисления текстурных атрибутов конструируют яркостные матрицы совпадения событий по изображениям в сконструированных многоугольниках и вычисляют начальные текстурные атрибуты на основании сконструированных яркостных матриц совпадения событий. 4. Способ по п.3, отличающийся тем, что для конструирования яркостных матриц совпадения событий дополнительно формируют объем углов наклона из объема сейсмических данных и применяют управление углом наклона при использовании значений углов наклона. 16 5. Способ по п.1, отличающийся тем, что для вычисления текстурных атрибутов осуществляют позиционирование скользящего окна на всем протяжении объема сейсмических данных и вычисление текстурных атрибутов в скользящем окне. 6. Способ по п.1, отличающийся тем, что дополнительно формируют объем доверительных значений, получаемых из сконструированной вероятностной нейронной сети. 7. Способ по п.6, отличающийся тем, что дополнительно осуществляют выбор доверительного уровня и корректируют размер скользящего окна для поддержания доверительных значений выше выбранного доверительного уровня. 8. Способ по п.1, отличающийся тем, что дополнительно отображают проклассифицированные текстурные атрибуты. 9. Способ по п.1, отличающийся тем, что используют сейсмические данные, которые содержат амплитуды сейсмических волн. 10. Способ по п.1, отличающийся тем, что используют сейсмические данные, которые содержат сейсмические атрибуты.