Способ построения геологических моделей подземных осадочных объёмов

010285 Родственные заявки Заявка на данное изобретение представляет собой международную заявку, регистрационный номерPCT/US2005/29753 на национальной фазе, поданную 19 августа 2005 г, по которой испрашивается приоритет предварительной заявки 60/609057 на патент США, поданной 10 сентября 2004 г. Область техники, к которой относится изобретение В общем, настоящее изобретение относится к области геологического моделирования. В частности,изобретение представляет собой способ построения модели геологических свойств на протяжении подземного объема. Уровень техники Геологическая модель является компьютерным представлением подземного объема геологического разреза, такого как нефтяной коллектор или бассейн осадконакопления. В нефтегазовой отрасли геологические модели обычно используют при таких работах, как определение мест закладки скважин, оценивание запасов углеводородов или планирование стратегий разработки коллектора. Геологические модели обычно являются входными данными в моделях движения флюидов, которые используются для проверки схем освоения и добычи с целью оптимизации возврата капиталовложений. Ключевым параметром при моделировании потока является пространственное распределение проницаемости, которое совместно со свойствами углеводородов и других флюидов, находящихся в подземном коллекторе, определяет продуктивность коллектора. В процессе геологического моделирования можно использовать данные многих различных видов,включая, но без ограничения, данные о свойстве породы, получаемые на основании скважинных и сейсмических данных, а также данные о структурных и стратиграфических поверхностях, которыми ограничиваются отдельные области в модели. Цель процесса заключается в создании представления о геологическом разрезе, которое является реалистическим, точным и согласующимся с доступными типовыми данными. Геологические модели могут иметь многочисленные различные формы. Чаще всего построение описательных или статических геологических моделей, предназначенных для применения в нефтяной отрасли, осуществляют в виде трехмерного набора модельных блоков (ячеек), которым приписывают геологические и/или геофизические свойства, такие как литология, пористость, акустический импеданс,проницаемость или водонасыщенность (в настоящем описании такие свойства в совокупности будут называться свойствами породы). Набор ячеек образует геологическую модель и является эквивалентом представляющего интерес подземного объема геологической среды. Для того чтобы не создавать избыточного числа ячеек, размеры ячеек обычно выбирают так, чтобы свойства породы в пределах ячейки были относительно одинаковыми. Назначение процесса геологического моделирования заключается в приписывании свойств породы каждой ячейке в геологической модели так, чтобы результирующая модель была точным отражением представляющего интерес подземного объема геологической среды. Для основанных на свойствах геоклеточных способов и объектно-ориентированного моделирования существуют два основных пути заполнения дискретизированного геологического объема. При геоклеточном подходе обычно используют геостатистические методы нахождения оценок (которые могут быть детерминированными или вероятностными). В этих способах учитывают расстояние, направление и пространственную непрерывность моделируемого свойства породы. Обычно в детерминированных методах нахождения оценок для каждого блока вычисляют оценку свойства породы, имеющую минимальную дисперсию. Вероятностными методами нахождения оценок получают распределения значений свойств породы и образуют набор геологических моделей для моделируемого свойства породы, при этом теоретически каждая модель является равновероятной. Пространственная непрерывность свойства породы может быть зарегистрирована на вариограмме. Вариограмма является хорошо известным средством для количественного выражения изменчивости свойства породы в зависимости от разноса и направления. Геостатистические модели имеют несколько важных преимуществ, заключающихся в том, что в них можно использовать широкое семейство существующих статистических алгоритмов, с ними легко согласуются контрольные пункты данных, такие как скважины, и они пригодны для моделирования и оптимизации свойств породы с использованием геофизических ограничивающих параметров, таких как трехмерные сейсмические данные. В патентах США 5838634, 6381543 и 6480790 раскрываются геоклеточные способы моделирования, реализуемые при обработке потоков, которые включают в себя повторяющиеся этапы оптимизации, предназначенные для приведения геологической модели в соответствие с геологическими и геофизическими данными, такими как каротажные данные, данные сейсмических исследований и данные по добыче флюида из геологического разреза и данные о давлении в нем. Большая часть пакетов программ для имеющего промышленное значение моделирования, включая PETREL,GOCAD и STRATAMODEL, содержит широкий спектр геостатистических средств, рассчитанных на удовлетворение требований промыслового геолога или инженера. Основной недостаток геоклеточных моделей заключается в том, что в целом ими не воспроизводятся точно структуры, наблюдаемые в осадочных системах, таких как реки, дельты, глубоководные каньоны и конусы выноса. Эти выводы являются важными, поскольку внутренняя структура осадочной системы может оказывать значительное влияние на качество и непрерывность коллектора.-1 010285 В объектно-ориентированных моделях объемы подземных коллекторов рассматривают совокупность геологических объектов, таких как каналы и осадочные языки. В патенте США 6044328 раскрыт один алгоритм объектно-ориентированного моделирования, посредством которого геологу или промысловому инженеру предоставляется возможность выбора геологических объектов из библиотеки аналогов для наилучшего согласования с моделируемым коллектором. Соответствие аналога оценивает оператор процесса на основании своих знаний в области геологии. В большей части пакетов программ для имеющего промышленное значение моделирования, включая PETREL, IRAP-RMS и GOCAD, объекты реализуются в виде объемных элементов, которыми имитируются каналы и языки, при этом используются упрощенные элементы, основанные на деформируемых пользователем конфигурациях, таких как полукольца и эллипсы. Объектно-ориентированное моделирование является наиболее полезным, когда отсутствует трехмерная пространственная информация, такая, которую получают из объема трехмерных сейсмических данных, или она имеется с низким разрешением. В случае простых геометрий, предусматриваемых согласно существующим способам, невозможно без затруднений собрать совокупную пространственную информацию, встречающуюся при современных трехмерных сейсмических исследованиях. Кроме того, в большей части современных способов обращается особое внимание на каналы и канальные комплексы как на основные элементы коллектора. Несмотря на то что канализированные системы относят к важным типичным представителям углеводородного коллектора, идентификация или моделирование каналов будет неполным без точного оценивания характеристик потенциально пористых и проницаемых песков,которыми они заполнены. Сами пески упакованы в виде округлых тел, подобных песчаным барьерам,наблюдаемым в существующих на сегодняшний день реках или дельтах. Они являются телами, которые образуют большую часть кремнекислых кластических коллекторов. Геологические модели, в которых учитывают и используют эти, естественным путем образующиеся, основные элементы должны обеспечивать получение более точных моделей геологического разреза. В основанных на концепции процессов способах геологического моделирования делались попытки воспроизводить стратиграфию геологического разреза путем моделирования или аппроксимирования физических процессов переноса и отложения осадков, наращивания осадочных отложений в хронологическом порядке. Результаты моделирования могут быть сверены с данными геологического разреза, а моделирование осуществляться повторно с использованием новых управляющих переменных при итерационном процессе до тех пор, пока не будет достигнуто приближенное соответствие с данными геологического разреза. В патентах США 5844799, 6205402 и 6246963 описаны три таких способа, в которых используют диффузионные или основанные на системе правил модели процессов для построения моделей в масштабе бассейна с ограниченной пространственной детализацией. В основанных на концепции процессов моделях обычно можно формировать реалистично выглядящие моделируемые отложения, но их обычно не используют для промышленного геологического моделирования, поскольку трудно скорректировать входные данные модели таким образом, чтобы для моделируемого отложения соблюдались ограничивающие условия для данных геологической среды. В отличие от геоклеточных и объектно-ориентированных моделей основанные на концепции процессов модели обычно невозможно эффективно оптимизировать с помощью реализуемых компьютером алгоритмов в тех случаях, когда становятся доступными новые статические или динамические данные. Наконец,при современных трехмерных сейсмических исследованиях для геологических моделей получают пространственную структуру, которую трудно воспроизводить применительно к основанным на концепции процессов моделям. Из вышеизложенного можно видеть, что существует необходимость в способе, в котором учитываются конфигурации и распределения свойств естественным образом образующихся осадочных отложений, однако же, который может быть легко привязан к имеющимся сейсмическим и скважинным данным. В описанном выше способе могут быть использованы основанные на концепции процессов объектно-ориентированные и геостатистические подходы. Предпочтительно, чтобы способом мог обеспечиваться автоматизированный процесс оптимизации, выполняемый компьютером, и получалась более точная модель подземного объема геологической среды с меньшими затратами дополнительного времени и усилий. Настоящим изобретением эта необходимость удовлетворяется. Сущность изобретения Согласно изобретению предусмотрен способ построения геологической модели, посредством которой на основании сейсмических данных точно определяют свойства породы в подземном объеме. Согласно этому варианту осуществления (а) на основании сейсмических данных идентифицируют контурные формы по меньшей мере одного геологического тела в подземном объеме; (b) на основании контурной формы геологического тела оценивают распределение зерен по размерам по меньшей мере в одной точке внутри геологического тела; (с) определяют свойства породы, включающие в себя пористость и проницаемость, используя найденное распределение зерен по размерам вместе с информацией о минералогии и истории захоронения. Согласно второму варианту осуществления также раскрыт способ построения геологической модели, посредством которой на основании сейсмических данных точно определяют свойства породы в под-2 010285 земном объеме. Согласно этому варианту осуществления:(а) выделяют из геологических данных контурные формы возможных геологических тел;(b) исследуют контурные формы, выбирая только те контурные формы, которые достоверно соответствуют преимущественно песчанистым геологическим телам;(с) корректируют размер контурных форм, чтобы сделать их, по существу, согласованными с дополнительными ограничениями на геологическую модель, такими как геофизические свойства интервала;(d) определяют изменяющееся в пространстве распределение зерен по размерам внутри контурных форм;(е) определяют свойства породы внутри геологических тел, используя определенное распределение зерен по размерам. Кроме того, свойства породы могут быть определены путем использования информации об истории захоронения вместе с распределением зерен по размерам. Согласно третьему варианту осуществления также раскрыт способ построения геологической модели, посредством которой на основании сейсмических данных точно определяют свойства породы в подземном объеме. Способ содержит этапы, при выполнении которых:(а) идентифицируют контурную форму по меньшей мере одного геологического тела на основании геологических данных;(b) оценивают поле осаждающего потока, из которого образовано по меньшей мере одно геологическое тело, на основании контурной формы;(с) оценивают распределение зерен по размерам внутри по меньшей мере части геологического тела, используя зависимость между оцененным полем потока, из которого образовано геологическое тело,и распределением зерен по размерам;(d) оценивают свойства породы осадочного бассейна в пределах по меньшей мере части осадочного бассейна, основываясь на распределении зерен по размерам. Краткое описание чертежей На чертежах: фиг. 1 - иллюстрация начальной стадии отложения осадочной системы потоком точечного источника; фиг. 2 - иллюстрация более поздней стадии осаждения такой осадочной системы, как показана на фиг. 1; фиг. 3 - иллюстрация еще более поздней стадии осаждения такой системы, как показана на фиг. 2; фиг. 4 - блок-схема последовательности операций согласно варианту осуществления изобретения; фиг. 5 - блок-схема последовательности операций согласно второму варианту осуществления изобретения; фиг. 6 - трехмерный вид куба сейсмических данных с обнаруженным зародышем; фиг. 7 - картографическое изображение трехмерного объема, показанного на фиг. 6; фиг. 8A - разрез и картографические изображения интервала возможного коллектора; фиг. 8B - картографическое изображение интервала возможного резервуара; фиг. 9 - картографическое изображение одного песчаного тела возможного коллектора, которое заполнено зернами с оцененными распределениями по размерам; фиг. 10 - схематическое представление захоронения и диагенетического изменения пробного блока коллекторного осадка с ходом времени захоронения, повышением давления или температуры; фиг. 11 - схематическое трехмерное изображение интервала коллектора, при этом в геологическую модель добавлены неколлекторный и низкокачественный коллекторный материал и фиг. 12 - блок-схема последовательности операций согласно третьему варианту осуществления изобретения. Подробное описание В нижеследующем подробном описании изобретение будет описано применительно к его предпочтительному варианту осуществления. Однако в той части, в какой нижеследующее описание является характерным для конкретного варианта осуществления или конкретного использования изобретения, оно подразумевается только иллюстративным. Поэтому изобретение не ограничено конкретным вариантом осуществления, описанным ниже, а точнее, изобретение включает в себя все варианты, модификации и эквиваленты, попадающие в рамки истинного объема прилагаемой формулы изобретения. Отложение обломочных осадочных тел обычно начинается с вытекания насыщенной наносами воды из узкого канала, такого как, например, устье реки, в открытую область, где поток замедляется. Место, где насыщенный наносами поток может расширяться и замедляться, обычно известно как впускной канал. Первоначально такие потоки свободно расширяются, а осадок откладывается по мере замедления потока. Затем по мере того, как откладываемый осадок нарастает по высоте, он начинает изменять поле потока. В конечном счете отложение становится достаточно большим, так что поток направляется в обход отложения. Следствием этого являются новые путь и впускной канал для поля потока в открытую область на расстоянии от или вблизи старого отложения. После этого процесс осадконакопления повторяется, а в системе образуется второе тело. К тому же в системе со временем могут быть эффективно-3 010285 сформированы несколько таких тел. В итоге, в ходе процесса образуется отложение, состоящее из стэков осадочных тел, которые являются удобной моделью структуры углеводородных коллекторов. Специалистам в данной области техники должно быть понятно, что характеристики потока, при которых образуются такие тела, могут быть исследованы на основании уравнения Навье-Стокса и принципов сохранения, примененных к насыщенному наносами потоку под неподвижной чистой водой, которая может быть увлечена им, такому как мутьевые потоки в глубоком океане. Такие исследования могут быть распространены на потоки без вертикального увлечения, например, на случай тонкого водного слоя под воздухом, и это распространение также находится в рамках объема настоящего изобретения. Конфигурацией и свойствами осаждающего потока определяются конфигурация и внутренние свойства осаждаемого геологического тела. Кроме того, поле потока связано как с контурной формой отложения, так и с распределением зерен по размерам во всех точках внутри отложения. Заявитель обнаружил, что (а) контурные формы геологических тел могут быть получены по геологическим данным, например по объемам сейсмических данных; (b) контурные формы геологических тел связаны с полем осаждающего потока, из которого образованы тела; (с) поле осаждающего потока в свою очередь определяет изменяющееся в пространстве распределение зерен внутри тела по размерам; и(d) распределением зерен по размерам, индивидуально или в сочетании с минералогической информацией и информацией об истории захоронения, определяются свойства породы, включающие в себя пористость и проницаемость, во всех точках внутри геологического тела. Объединение результатов этих научных наблюдений приводит к способу построения геологической модели, в которой пористость и проницаемость на всем протяжении подземного объема точно определены на основании сейсмического изображения подземного объема. На фиг. 1-3 показан процесс образования откладываемых водой обломочных осадочных систем при наличии точечного источника. На фиг. 1 представлена самая ранняя стадия осаждения, когда обогащенный осадками поток 10 выходит из впускного канала 12 в открытую область 14. По мере того как поток расширяется и замедляется, он теряет способность переносить осадок, и под потоком образуется отложение 16. На фиг. 2 показана более поздняя стадия отложения системы, ранее изображенной на фиг. 1. На фиг. 2 отложение 16 из фиг. 1 стало более крупным, и оно обозначено позицией 16'. Отложение 16' является столь крупным, что оно отклоняет поток 20. Отложение 16' больше не нарастает, но теперь отклоненный поток создает два новых отложения 22 и 24. На фиг. 3 показана более поздняя стадия отложения системы, изображенной на фиг. 2. Теперь поток 30 заключен в каналы 32 и 34, по которым он переносится мимо законченных отложений 16', 22' и 24', и при этом наращиваются отложения 36 и 38. Теперь будет описан первый вариант осуществления. Согласно фиг. 4 в этом варианте осуществления изобретения предусмотрены процедуры для определения свойств породы в подземной области путем использования сейсмического изображения области. Как пояснено фиг. 4, идентифицируют (этап 401) контурную форму по меньшей мере одного геологического тела в подземной области. Определяют (этап 402) распределение зерен по размерам по меньшей мере в одной точке внутри геологических тел, идентифицированных на этапе 401. Наконец, определяют (этап 403) по меньшей мере одно свойство породы,используя распределение зерен по размерам, найденное на этапе 402. Свойства породы включают в себя пористость и проницаемость. Отдельные этапы будут описаны более подробно в нижеследующих абзацах. Этап 401 заключается в идентификации контурных форм геологических тел и может быть осуществлен по геологическим данным с помощью ряда известных из уровня техники способов. Для примера,геологические данные могут включать в себя, но без ограничения ими, сейсмические данные, данные или результаты исследования обнажения пород, данные анализа образцов керна, данные о свойствах коллектора и/или каротажные данные. Сейсмические данные являются предпочтительным видом данных для определения контурной формы геологического тела. Геологические данные могут быть получены путем использования многих различных способов сбора данных из известного уровня техники, которые должны быть известны специалистам в данной области техники. На коммерческой основе сейсмическими или геофизическими компаниями могут быть собраны и предоставлены эти данные, а также ранее собранные данные, при этом данные могут быть предоставлены в аналоговом или цифровом виде. Все такие способы получения геологических данных подразумеваются входящими в рамки объема формулы изобретения. Существуют три основных способа нахождения контурной формы геологического тела по сейсмическим данным. Первый способ заключается в ручной идентификации в данных стратиграфических поверхностей, которыми задаются границы стратиграфических интервалов и образуются верхняя и нижняя границы геологических тел. Обычно геологи вручную идентифицируют стратиграфические поверхности путем распознавания и идентификации сейсмических отражений как горизонтов, которые коррелируют с изменениями свойств породы, в том числе с изменениями стратиграфических поверхностей. Подробная информация относительно этого процесса хорошо известна специалистам в данной области техники. Второй способ заключается в автоматизированном продвижении стратиграфических поверхностей,которыми задаются границы стратиграфических интервалов. Обычно этот способ осуществляют по-4 010285 сейсмическим данным путем компьютерной экстраполяции известных или предполагаемых стратиграфических поверхностей. Примеры включают в себя программы расчета когерентности сейсмических волн и границ разделов. Подробная информация относительно этого процесса хорошо известна специалистам в данной области техники и содержится в коммерческих прикладных программах, таких какVOXELGEO, GEOFRAME, SEISWORKS3D и PETREL. Третий способ заключается в автоматической идентификации стратиграфически реалистических геологических тел как областей, имеющих подобные значения сейсмических атрибутов. Предпочтительное воплощение этапа 401 заключается в автоматической идентификации геологических тел на основании способов, описанных в патентах США 6674689, 6823266 и 7024021. Посредством этих способов обеспечиваются обнаружение и каталогизация геологических тел с более высокой разрешающей способностью, чем это возможно способами интерпретации наземных наблюдений. Кроме того, этими способами для исследователей обеспечивается возможность извлечения представительной стратиграфической информации из объемов в широком диапазоне пространственной разрешающей способности путем использования иерархического подхода, при котором эти данные исследуются в нескольких масштабах. В патенте США 7024021 раскрыт способ осуществления обнаружения на стратиграфической основе зародыша в сейсмических данных, при этом в способе удовлетворяются определенные характеристические критерии, являющиеся стратиграфически достоверными. Критерии отражают многослойный характер геологической среды, так что получаемые в результате сейсмические объекты являются стратиграфически корректными. Кроме того, этим способом можно точно определять размер и конфигурацию конкретного сейсмического объекта путем наращивания объекта вокруг выбранной зародышевой точки в объеме сейсмических данных. Объект наращивают путем продвижения объекта в соседние трассы, если соседняя трасса удовлетворяет выбранным критериям. В патенте США 6823266 раскрыт способ выполнения объектно-ориентированного анализа связности в объеме трехмерных данных. Этот способ обеспечивает возможность прогнозирования связности сейсмического объекта, идентифицированного путем обнаружения зародыша или аналогичным путем. В этом способе при условии сравнения данных двух или более проходов каждый идентифицируемый сейсмический объект может быть отслежен на предмет изменений размера, пространственного положения или связей с другими сейсмическими объектами в зависимости от выбранного порогового значения или другого задаваемого пользователем ограничивающего параметра. Это позволяет оператору выбирать пороговое значение связности и отслеживать осадочное тело по этому пороговому значению. В патенте США 6674689 раскрыт способ классификации морфологии сейсмических объектов. В этом способе вычисляют геометрические статистические параметры и используют их для классификации морфологии сейсмического объекта. Сейсмически определенные тела, которые имеют конфигурации и внутренние атрибуты, указывающие на сильный потенциал углеводородного коллектора (или другого элемента, представляющего интерес), могут быть оценены и выделены на фоне второстепенных сейсмических объектов. Согласно предпочтительному варианту осуществления осадочное тело получают путем наращивания тела от зародышевой точки, используя обнаружение на стратиграфической основе зародыша в объеме сейсмических данных. Наращивание тела удовлетворяет определенным характеристическим критериям, рассчитанным на то, чтобы сделать тело стратиграфически корректным. Размер и конфигурация конкретного сейсмического объекта могут быть точно определены. Затем осуществляют объектноориентированный анализ связности по отношению к полученному телу. Идентифицированный сейсмический объект может быть отслежен на предмет изменений размера, пространственного положения или связей с другими сейсмическими объектами в зависимости от выбранного порогового значения. В заключение осуществляют морфологическую классификацию полученных объектов. На этапе 401 могут использоваться объемы двумерных и трехмерных данных многих видов, включая, но без ограничения ими, амплитуды сейсмических волн (в ближнем поле, дальнем поле или полные суммы), сейсмические атрибуты (длина замкнутого контура из участков профилей, мгновенные значения амплитуды и фазы), сейсмический импеданс (ограниченный по полосе и полный), а также производные параметры, такие как данные о сейсмических фациях и границах раздела. Этап 402 включает в себя определение распределения зерен по размерам в любой точке внутри геологических тел. Этот этап может быть осуществлен с помощью ряда технологий. Для некоторых технологий необходимы дополнительные данные, и предпочтительное воплощение этого этапа зависит от доступности этих вспомогательных данных, а также от характера геологических тел. Одна технология представляет собой способ, согласно которому осуществляют связь трехмерного контура геологических тел с распределением находящихся внутри зерен по размерам. Этот способ описан в патенте США 7062383. Он является предпочтительным воплощением этапа, если по сейсмическим данным могут быть идентифицированы полностью сохранившиеся (неэродированные) тела, и должен применяться, если трехмерная форма тел может быть точно определена. В этом способе внутренние свойства отложенного водой осадочного тела могут быть определены по трехмерной конфигурации отложения. Согласно одному воплощению предусмотрены решение уравнений, связывающих идентифи-5 010285 цированную трехмерную конфигурацию осадочного тела с полем потока, посредством которого образовано тело, и использование поля потока для нахождения свойств осадочного тела по меньшей мере в одной точке внутри тела. Вторая технология представляет собой способ, согласно которому осуществляют связь двумерного контура или очертания в плане геологических тел с распределением находящихся внутри зерен по размерам. Этот способ описан в патенте США 6885941. Способ должен применяться только в случае,если двумерная конфигурация геологических тел может быть точно определена. Путем использования этого способа свойства отложенного водой осадочного отложения могут быть спрогнозированы в любом месте по изопахите отложения. Согласно одному воплощению способ включает в себя:(а) определение изопахиты осаждения для исследуемого отложения;(b) определение свойств потока флюида во впускном канале исследуемого отложения и(с) определение свойства отложения в любой точке внутри отложения на основании моделирования потока флюида. Третья технология представляет собой способ, согласно которому осуществляют связь мощности геологических тел с распределением находящихся внутри зерен по размерам. Этот способ описан в публикации международной заявки WO 2004/083896. Этот способ должен применяться только в случае, если останцы тел имеются в наличии или если тела продолжаются за пределы области сейсмических данных,и при условии, что имеются результаты измерения распределения зерен по размерам или распределение может быть оценено в одной точке внутри каждого тела. Этим способом прогнозируют распределение зерен по размерам в любом заранее намеченном месте внутри отложенного водой осадочного отложения. Во-первых, определяют вертикальную мощность осадочного отложения в заранее назначенном местоположении, а также вертикальную мощность и распределение зерен по размерам во втором местоположении, отличающемся от заранее назначенного местоположения. Во-вторых, определяют величину расстояния между двумя местоположениями. В заключение величину расстояния, вертикальную мощность в обоих местоположениях и распределение зерен по размерам во втором местоположении используют для вычисления распределения зерен по размерам в заранее назначенном местоположении. Четвертая технология включает в себя определение свойств геологических тел, через которые проходит скважина, и заполнение более крупной геологической области синтезированными геологическими телами, имеющими распределение свойств, аналогичное тому, которое имеют тела, через которые проходит скважина. Этот способ должен применяться в случае, когда самые малоразмерные геологические тела являются слишком малыми, чтобы быть видимыми в сейсмических данных, при условии, что могут быть получены данные анализа образца керна или другой показатель, относящийся к мощности пласта и распределению зерен по размерам. Способ определения свойств геологических тел, через которые проходит скважина, описан в публикации международной заявки WO 2004/093521. Этим способом можно определять свойства отложенного водой осадочного тела по результатам измерений распределения зерен по размерам в одном месте в теле и мощности отложения. Согласно одному воплощению способ включает в себя:(а) определение свойств потока в месте измерения;(b) экстраполяцию свойств потока обратно к впускному каналу, по которому был выпущен осаждающий поток;(с) определение по меньшей мере одного свойства отложенного водой осадка на всем протяжении осадочного тела путем моделирования свойств потока с использованием экстраполированных свойств потока во впускном канале из этапа (b) в качестве граничного условия. После того как свойства осадочного тела в окрестности скважины определены, выявленные сейсмическим способом геологические тела по всей системе могут быть заполнены более мелкими телами,имеющими соответствующее статистическое распределение свойств, определяемое на основании распределения свойств, получаемых при исследовании в скважине. Этот способ описан в публикации международной заявки WO 2005/104033. Пятая технология включает в себя моделирование потоков флюида, из которых образованы тела и корректировку граничных условий потока до тех пор, контур моделируемого отложения не станет очень близко походить на контур отложения, обнаруженный сейсмическим способом. Подходящая модель для этого способа раскрыта в находящейся в процессе одновременного рассмотрения предварительной заявке на патент США, поданной одновременно с заявкой на данное изобретение, под названием Способ оценивания свойств осадочных отложений путем числового моделирования процессов осаждения. В описании раскрыт способ оценивания свойств породы путем числового моделирования процессов осаждения. Согласно этой технологии решают уравнения динамики флюида по схеме контрольных объемов для имитирования физики потока и отложения. Согласно одному воплощению эта технология числового моделирования включает в себя:(а) решение системы двумерных уравнений картографического изображения с зависящими от времени коэффициентами, по меньшей мере, для количества движения потока, высоты потока, концентрации взвешенного осадка и уноса вышележащей воды;(b) вычисление чистого осадочного отложения для каждого места картографического изображения-6 010285 путем с использованием свойств потока и(с) регистрацию изменчивости во времени чистого осадочного отложения. Поскольку вычислительные возможности повышаются, эти модели могут обладать большой прогнозирующей способностью при моделировании свойств породы и потенциальных регуляторов движения флюидов в моделях углеводородных коллекторов. Предпочтительный способ воплощения этой пятой технологии сочетается с основанным на гидродинамике гридинге (гидросетках), предназначенным для геологического моделирования свойств породы,в том числе распределения зерен по размерам. Этот способ раскрыт в публикации международной заявкиWO 2006/007466. В этом способе верхние части ячеек геологической модели соответствуют поверхностям постоянного времени, определяемым в процессе числового моделирования. В таком случае свойства породы бассейна представляют в этой, основанной на гидродинамике сетчатой структуре. Этап 403 включает в себя прогнозирование свойств породы внутри геологических тел на основании распределения зерен по размерам. Это может быть осуществлено рядом способов, но без ограничения ими, включающих в себя использование прямого вычисления свойств породы на основании распределения зерен по размерам, минералогии и истории захоронения. В отсутствие данных о размерах зерен в современных способах геологического моделирования обычно используют сводные двумерные графики данных о пористости и проницаемости, полученные на основании вырезанных для исследования образцов керна, подобранных по фациям осадконакопления, для оценивания проницаемости на протяжении всей модели. Значения и пространственное распределение пористости обычно оценивают геостатистическими способами в сочетании с анализом фаций осадконакопления по двумерным или трехмерным сейсмическим данным и данным бурения. Специалистам в области геологического моделирования хорошо известны эти технологии. Предпочтительный способ описан в публикации международной заявки WO 2005/104002. Согласно этому способу моделируют отложение, уплотнение и цементацию песчаника. Осаждение, уплотнение и цементацию зерен моделируют на основании определенных распределения зерен по размеру, минералогии и истории захоронения осажденных зерен. По результатам моделирования могут быть вычислены свойства песчаника, такие как пористость и проницаемость. Способ обеспечивает возможность моделирования нескольких минералогии в течение истории захоронения. Теперь будет описан второй вариант осуществления. Обратимся к фиг. 5, на которой этот вариант осуществления включает в себя процедуры, предназначенные для определения свойств породы в подземной области путем использования сейсмического изображения области. Как показано на фиг. 5, идентифицируют (этап 501) контурные формы геологических тел в пределах интервала коллектора. Каждую контурную форму оценивают и отбрасывают (этап 502) при наличии малой вероятности принадлежности формы к преимущественно песчанистому геологическому телу. Контурные формы корректируют (этап 503) по размерам, чтобы сделать их соответствующими дополнительным ограничениям, накладываемым на коллектор. Определяют (этап 504) распределение зерен в любой точке внутри форм геологических тел. Используя распределение зерен по размерам, определяют (этап 505) свойства породы. Объем в пределах интервала коллектора, который не относится к форме геологического тела, заполняют (этап 506) материалом с низкой проницаемостью. Этап 501 включает в себя идентификацию контурных форм геологических тел по геологическим данным. Способы, описанные относительно этапа 401 из фиг. 4, являются теми же самыми и осуществляются аналогичным образом на этапе 501 из фиг. 5. Предпочтительный способ заключается в автоматическом выделении контурных форм геологических тел на основании объема трехмерных сейсмических данных. На этапе 502 оценивают каждый полученный контур с целью определения, является ли он преимущественно песчанистым геологическим телом. Это осуществляют способами с использованием сейсмических атрибутов и с помощью геометрического анализа морфологии выделенных сейсмических объектов. Те выделенные контуры, которые не считают преимущественно песчанистыми геологическими телами, отбрасывают. Этап 503 заключается в коррекции размера преимущественно песчанистых геологических тел для привязки к геофизическим свойствам интервала, определенным по сейсмическим или скважинным данным. Обычно этот этап включает в себя изменение контуров тел путем расширения или сокращения так,чтобы коэффициент песчанистости интервала коллектора, оцениваемый по сейсмическим или скважинным данным, был, по существу, равен части вертикального расстояния на протяжении интервала в границах контуров тел. Процесс расширения или сокращения может включать в себя добавление постоянной мощности к контуру геологического тела или изъятие из него. В качестве альтернативы постоянный процент мощности может быть добавлен в каждом месте или изъят из него. Более прогрессивные технологии включают в себя добавление или удаление ячеек до достижения близкого согласования с конфигурациями песчаных тел, о которых получено представление из имитационных моделей, из результатов исследований обнажения пород, результатов исследований существующих на сегодняшний день условий осадконакопления и из объемов трехмерных сейсмических данных с высоким разрешением. Специалисты в данной области техники должны признавать наличие других способов коррекции контуров тел с-7 010285 целью согласования с геофизическими свойствами интервала. Этап 504 включает в себя нахождение распределения зерен по размерам по меньшей мере в одной точке внутри геологического тела. Способы, описанные относительно этапа 402 из фиг. 4, являются теми же самыми и осуществляются аналогичным образом на этапе 504 из фиг. 5. Этап 505 включает в себя прогнозирование свойств породы внутри геологического тела путем использования информации о распределении зерен по размерам из этапа 504. Предпочтительное воплощение заключается в использовании распределения зерен по размерам и минералогии вместе с историей захоронения для моделирования осадкообразования, уплотнения и цементации зерен. Такое моделирование может обеспечивать получение оценок пористости, проницаемости и других свойств породы в точках, в которых известно распределение зерен по размерам. Способы, описанные относительно этапа 403 из фиг. 4, являются теми же самыми и осуществляются аналогичным образом на этапе 505 из фиг. 5. Этап 506 заключается в заполнении пространства в модели, не включенного в формы геологических тел, материалом с низкой проницаемостью. Такое заполнение может быть осуществлено путем использования обычных геостатистических технологий, хорошо известных специалистам в данной области техники. Эти зоны низкой проницаемости могут содержать мелкий песок или ил, но обычно они содержат значительные количества обломочной глины (в виде очень мелкозернистых частиц, осажденных с песком или алевритовыми зернами). Эта глина часто заполняет поры между песчаными зернами, уменьшая пористость и проницаемость песка. Глина может быть сконцентрирована в слоях (слоистая глина) или довольно равномерно распределена по образцу (рассеянная глина). Согласно одному воплощению для распределения глины в моделируемой породе заявителем использованы алгоритмы кластеризации из способа Flo-Рас (описанного в международной заявке WO 2005/104002). Пользователь определяет количество глины, подлежащей введению в породу, размер и плотность глинистых кластеров и конфигурацию распределения (слоистую или рассеянную). Каждый фактор из конфигурации распределения, распространенности и плотности кластеров влияет на поровую структуру и поэтому на характеристики потока осадка. Каждый фактор из внедрения или осаждения глины, роста цемента и физического уплотнения оказывает влияние на снижение пористости и проницаемости с увеличением погребения. Теперь будет рассмотрен гипотетический пример. В этом примере выполняют обнаружение зародыша в кубе сейсмических данных. На фиг. 6 представлено трехмерное изображение куба 61 сейсмических данных с обнаруженным зародышем, иллюстрирующее отдельные предполагаемые песчаные тела 63, а также верхний 65 и нижний 67 ограничивающие горизонты интервала 69 коллектора, заданные интерпретатором. После того как тела идентифицируют путем обнаружения зародыша, их классифицируют на основании их морфологии и сейсмических атрибутов. На фиг. 7 показано картографическое изображение 75 трехмерного объема, показанного на фиг. 6. Песчаные тела 63 из фиг. 6 классифицируют на основании их морфологии и их сейсмических атрибутов на категории, которые отражают относительный потенциал коллектора. На основании морфологии и сейсмических атрибутов полагают, что оба песчаных тела, показанных на фиг. 7, достоверно являются пористыми, и полагают, что с малой вероятностью песчаные тела 73 будут жизнеспособными телами коллектора. Согласно этому варианту осуществления после классификации морфологии тел выполняют объектно-ориентированный анализ связности относительно тел в объеме трехмерных данных. На фиг. 8A и 8B показаны соответственно разрез и картографическое изображение возможного интервала 81 коллектора,ограниченного на предшествующих этапах этого примера. На фиг. 8A представлен разрез, подтверждающий возможную картину чередования слоев в углеводородном коллекторе 83, содержащем несколько расположенных последовательно по вертикали песчаных интервалов 85 между верхним горизонтом 87 коллектора и нижним горизонтом 89 коллектора. Как показано на фиг. 8B, пределы коэффициентов песчанистости могут быть определены для каждого столбца ячеек в разрезе и для каждой точки карты. На этом чертеже показаны три диапазона песчанистости. Коэффициент 82 песчанистости в пределах от 30 до 70%, коэффициент 84 песчанистости от 5 до 30% и коэффициент 86 песчанистости от 0 до 5%. Картой предоставляется одна возможность уточнения песчанистости путем использования многоугольников для разделения интервала коллектора на области относительных потенциалов песчанистости. Такая карта может быть столь же простой, как и карта на фиг. 8 В, при этом крупными многоугольниками обозначаются зоны относительного потенциала коллектора, или на ней могут быть показаны намного более сложные распределения. В более сложных моделях песчанистости могут иметься изменяющиеся в пространстве картины, получаемые при измерениях атрибутов интервала от трассы к трассе, осуществляемых на сейсмических объемах или петрофизических объемах, определяемых сейсмическим способом. В случаях, когда скопление обнаруженных тел коллектора не является полным, чтобы соответствовать значению песчанистости интервала коллектора с карты, скопление может быть скорректировано путем итеративного изменения протяженности ареала и вертикальной мощности отдельных тел коллектора. Затем прогнозируют распределения зерен по размерам и связанные с ними свойства пород тел. На фиг. 9 показано картографическое изображение одного песчаного тела 91 возможного коллектора с показом впускного канала 93 потока и направления 95 древнего потока, из которого осаждено песчаное тело.-8 010285 Песчаное тело заполнено зернами с распределениями по размерам, основанными на ожидаемом размере зерен на входе и конфигурации тела. Распределение по размерам зерен, переносимых потоком во впускном канале, может быть оценено по скважинным данным или данным об обнажении пород в районе исследования. Оцененные распределения зерен по размерам получают, используя один из алгоритмов для решения обратной задачи с получением модели размера зерен по конфигурации, описанных выше для этапов 402 и 504 из фиг. 4 и 5 соответственно. На рассматриваемом чертеже изолиниями размеров зерен охвачены крупнозернистые 90, среднезернистые 92, мелкозернистые 94 и очень мелкозернистые 96 пески, но фактическим результатом способов решения обратной задачи с получением модели размера зерен по конфигурации являются полные распределения зерен по размерам во всех точках в теле. После решения обратной задачи получившуюся модель размеров зерен привязывают к глубине,давлению и температурному режиму, соответствующим исследуемому интервалу. Анализ истории захоронения вместе с информацией о минералогиях исходных зерен и распределениях зерен по размерам может быть использован для прогнозирования изменений осадков из захоронения, обусловленных ходом времени, воздействием давления и температуры. Новые способы моделирования, такие как Flo-Рас (описанный в публикации международной заявки WO 2005/104002), обеспечивают возможность детализированного моделирования осадкообразования, уплотнения, цементации и заполнения глиной, подобранной к тонким сечениям осадков. На фиг. 10 показано схематическое представление захоронения и диагенетического изменения пробного блока 101 коллекторных осадков 103 в зависимости от возрастающего захоронения 105 с течением времени при сопутствующих повышениях давления и температуры 107. Исходный блок 101 с упаковкой зерен имеет высокую пористость и проницаемость при наличии глин 109 низкой плотности и отсутствии цемента. Конечный блок 102 коллекторной породы имеет пониженные значения пористости и проницаемости вследствие механического уплотнения песчаных зерен 103, а также химического растворения и выпадения в осадок заполняющих поры цементов 104 и уплотненных глин 109. Эта история захоронения обеспечивает возможность получения детальных оценок свойств породы на протяжении всей истории осадочного бассейна. На фиг. 11 показано трехмерное изображение, соответствующее последнему этапу в этой примерной последовательности операций, на котором неколлекторный и низкокачественный коллекторный материал 111 введены в конечную геологическую модель путем использования геостатистических способов, увязанных со скважинными данными и геофизическими данными интервала. Кроме того, показаны песчаные тела 113, а также кровля 115 коллектора и подошва 117 коллектора. Эти способами обеспечивается получение детализированной модели интервала коллектора, и они могут использоваться для оценивания свойств породы и характеристик продуктивности интервала коллектора. Кроме того, способы,описанные выше, могут быть использованы для прогнозирования свойств осадочного бассейна, а не только интервала коллектора, как описывалось в этом примере. В дополнение к вариантам осуществления, описанным выше, на фиг. 12 дополнительно представлен третий вариант осуществления. Что касается фиг. 12, то этот вариант осуществления включает в себя процедуры, предназначенные для определения свойств породы в подземной области путем использования сейсмического изображения области. Как показано на фиг. 12, на этапе 1201 идентифицируют контурную форму по меньшей мере одного геологического тела в подземной области. На этапе 1202 на основании контурной формы тела определяют поле осаждающего потока, из которого образовано геологическое тело. На этапе 1203, используя зависимость между оцененным полем потока и распределением осажденных зерен по размерам, определяют распределение зерен по размерам внутри геологического тела. На этапе 1204 на основании оцененного распределения зерен по размерам определяют свойства породы внутри геологического тела. Этап 1201 заключается в идентификации контурных форм геологических тел и может быть выполнен в соответствии с рядом способов, описанных относительно этапа 401. Этап 1202 включает в себя определение поля потока, из которого осаждено геологическое тело,осуществляемое по контурной форме тела. Согласно одному варианту осуществления этот этап выполняют так, как описано в патенте США 7062383, который включен в настоящее описание посредством ссылки, при этом получают уравнение, связывающее трехмерную конфигурацию осадочного тела с полем потока, из которого осаждено тело. Согласно второму варианту осуществления этот этап выполняют способом, описанным в патенте США 6885941, который включен в настоящее описание посредством ссылки. В соответствии с этим патентом свойства потока во впускном канале определяют по двумерному картографическому изображению контура геологического тела, а поле потока, из которого осаждено тело, определяют на основании свойств потока во впускном канале. Согласно третьему варианту осуществления этот этап выполняют путем моделирования потоков флюида, из которых может образовываться тело, и уточнения граничных условий потока до тех пор, пока моделируемый контур отложения не будет в наибольшей степени походить на контур, полученный на этапе 1201. Подходящий способ моделирования таких потоков описан в публикации международной заявки WO 2006/007466, которая включена в настоящее описание посредством ссылки. Согласно этим воплощениям свойства потока, задаваемые в виде части поля потока, обычно включают в себя, по меньшей мере, скорость потока, концентрацию-9 010285 взвешенного осадка в потоке и распределение зерен взвешенного осадка по размерам. Этап 1203 включает в себя определение распределения зерен по размерам в отложении на основании свойств потока. Это определение осуществляют, используя эмпирические зависимости между распределением осажденных зерен по размерам и свойствами потока. Соответствующие эмпирические зависимости рассмотрены в публикации международной заявки WO 2006/007466 и в патенте США 6885941, которые упомянуты выше. Этап 1204 включает в себя прогнозирование свойств породы внутри геологического тела на основании распределения зерен по размерам. Способы осуществления такого прогнозирования рассмотрены применительно к этапу 403. ФОРМУЛА ИЗОБРЕТЕНИЯ 1. Способ построения геологической модели подземного объема, содержащий этапы, при выполнении которых:(a) идентифицируют по меньшей мере одну контурную форму по меньшей мере одного геологического тела в подземном объеме на основании геологических данных, относящихся к подземному объему;(b) оценивают распределение зерен по размерам внутри по меньшей мере части по меньшей мере одного геологического тела, используя по меньшей мере одну контурную форму по меньшей мере одного геологического тела;(c) оценивают свойства породы подземного объема по меньшей мере в части осадочного бассейна на основании распределения зерен по размерам. 2. Способ по п.1, в котором геологические данные выбирают из сейсмических данных, данных об обнажении пород, данных анализа образцов керна, каротажных данных и любого их сочетания. 3. Способ по п.1, в котором по меньшей мере одну контурную форму по меньшей мере одного геологического тела идентифицируют способом, выбираемым из ручной идентификации стратиграфических поверхностей, автоматизированного продвижения стратиграфических поверхностей, автоматической идентификации стратиграфически реалистических геологических тел в качестве областей, имеющих подобные значения сейсмических атрибутов, и любого их сочетания. 4. Способ по п.3, в котором по меньшей мере одну контурную форму по меньшей мере одного геологического тела идентифицируют с помощью автоматической идентификации стратиграфически реалистических геологических тел конкретным способом, выбираемым из осуществления на стратиграфической основе обнаружения зародыша в сейсмических данных, осуществления объектно-ориентированного анализа связности в сейсмических данных, классификации морфологий сейсмических объектов и выделения приемлемых морфологий и любого их сочетания. 5. Способ по п.1, в котором на этапе (а) дополнительно получают сейсмическое тело на основании сейсмических данных путем наращивания по меньшей мере одного сейсмического тела от зародышевой точки, используя обнаружение на стратиграфической основе зародыша в массиве сейсмических данных,осуществляют объектно-ориентированный анализ связности по отношению по меньшей мере к одному полученному телу и классифицируют по меньшей мере одно полученное тело с помощью морфологии и выбирают сейсмическое тело, если морфология является приемлемой. 6. Способ по п.1, в котором устанавливают связь по меньшей мере одной контурной формы по меньшей мере одного геологического тела с распределением зерен по размерам, при этом установление связи выбирают из группы, состоящей из установления связи трехмерного контура по меньшей мере одного геологического тела с распределением находящихся внутри зерен по размерам, установления связи двумерного контура на виде в плане по меньшей мере одного геологического тела с распределением находящихся внутри зерен по размерам, установления связи мощности по меньшей мере одного геологического тела с распределением находящихся внутри зерен, определения свойств отложенного водой осадочного тела по результатам измерения распределения зерен по размерам и мощности отложения в одном месте по меньшей мере в одном геологическом теле, заполнения сейсмически определенных геологических тел меньшими телами, имеющими распределение свойств, обнаруживаемое в скважине, моделирования потоков флюида, из которых образовано по меньшей мере одно геологическое тело, с уточнением граничных условий потоков флюида до тех пор, пока контур моделируемого отложения не будет походить на обнаруженный сейсмическим способом контур отложения, и любого сочетания из них. 7. Способ по п.1, в котором при оценивании свойств породы подземного объема используют, по меньшей мере, распределение зерен по размерам по меньшей мере в одной точке внутри подземного объема. 8. Способ по п.1, в котором при выполнении этапа (с) дополнительно оценивают распределение зерен по размерам и минералогический состав зерен в песчанике, моделируют осаждение зерен на основании распределения зерен по размерам и минералогического состава зерен, моделируют уплотнение зерен и моделируют цементацию зерен.- 10010285 9. Способ по п.1, в котором на этапе (b) дополнительно образуют верхние и нижние границы ячеек сетки, соответствующие поверхностям периода отложения осадка, определенным путем моделирования отложения осадка, образуют боковые границы ячеек сетки и моделируют свойства породы подземного объема, используя ячейки сетки, связанные с верхними границами ячеек сетки, нижними границами ячеек сетки и боковыми границами ячеек сетки. 10. Способ построения геологической модели подземного объема, содержащий этапы, при выполнении которых:(a) идентифицируют по меньшей мере одну контурную форму возможного геологического тела на основании геологических данных, относящихся к подземному объему;(b) выбирают из идентифицированной по меньшей мере одной контурной формы только по меньшей мере одну контурную форму, достоверно соответствующую преимущественно песчанистым геологическим телам;(c) корректируют по меньшей мере одну контурную форму возможного геологического тела, чтобы сделать возможное геологическое тело, по существу, согласованным с дополнительными ограничениями на свойства тела;(d) оценивают распределение зерен по размерам внутри по меньшей мере части возможного геологического тела, руководствуясь по меньшей мере одной контурной формой возможного геологического тела;(e) оценивают свойства породы внутри по меньшей мере части возможного геологического тела на основании распределения зерен по размерам. 11. Способ по п.10, в котором дополнительно заполняют по меньшей мере часть объема в осадочном отложении, которое не соответствует преимущественно песчанистым геологическим телам, материалом с низкой проницаемостью. 12. Способ по п.10, в котором геологические данные выбирают из группы, состоящей из сейсмических данных, данных об обнажении пород, данных анализа образцов керна, каротажных данных и любого сочетания из них. 13. Способ по п.10, в котором по меньшей мере одну контурную форму возможного геологического тела идентифицируют способом, выбираемым из группы, состоящей из ручной идентификации стратиграфических поверхностей, автоматизированного продвижения стратиграфических поверхностей, автоматической идентификации стратиграфически реалистических геологических тел в качестве областей,имеющих подобные значения сейсмических атрибутов, и любого сочетания из них. 14. Способ по п.10, в котором по меньшей мере одну контурную форму возможного геологического тела идентифицируют с помощью автоматической идентификации стратиграфически реалистических геологических тел конкретным способом, выбираемым из группы, состоящей из осуществления на стратиграфической основе обнаружения зародыша в сейсмических данных, осуществления объектноориентированного анализа связности в сейсмических данных, классификации морфологий сейсмических объектов и выделения приемлемых морфологий и любого сочетания из них. 15. Способ по п.10, в котором этап (с) выполняют путем изменения по меньшей мере одной контурной формы так, чтобы коэффициент песчанистости по меньшей мере одного интервала коллектора, оцениваемый на основании сейсмических или скважинных данных, был, по существу, равен доле вертикального расстояния на протяжении по меньшей мере одного интервала коллектора, заключенного в контурах тел. 16. Способ по п.10, в котором этап (d) выбирают из группы, состоящей из установления связи трехмерного контура по меньшей мере одного геологического тела с распределением находящихся внутри зерен по размерам, установления связи двумерного контура на виде в плане по меньшей мере одного геологического тела с распределением находящихся внутри зерен по размерам, установления связи мощности по меньшей мере одного геологического тела с распределением находящихся внутри зерен, определения свойств отложенного водой осадочного тела по результатам измерения распределения зерен по размерам и мощности отложения в одном местоположении внутри по меньшей мере одного геологического тела, заполнения сейсмически определенных геологических тел меньшими телами, имеющими распределение свойств, обнаруживаемое в скважине, моделирования потоков флюида, из которых образовано по меньшей мере одно геологическое тело, с уточнением граничных условий потоков флюида до тех пор, пока контур моделируемого отложения не будет походить на обнаруженный сейсмическим способом контур отложения, и любого сочетания из них. 17. Способ по п.10, в котором при выполнении этапа (е) дополнительно оценивают распределение зерен по размерам и минералогический состав зерен в песчанике, моделируют осаждение зерен на основании распределения зерен по размерам и минералогического состава зерен, моделируют уплотнение зерен и моделируют цементацию зерен. 18. Способ по п.10, в котором при выполнении этапа (d) дополнительно образуют верхние и нижние границы ячеек сетки, соответствующие поверхностям периода отложения осадка, определенным путем моделирования отложения осадка, образуют боковые границы ячеек сетки и моделируют свойства породы бассейна, используя ячейки сетки, связанные с верхними границами ячеек сетки, нижними границами- 11010285 ячеек сетки и боковыми границами ячеек сетки. 19. Способ построения модели подземного объема, содержащий этапы, при выполнении которых:(a) идентифицируют контурную форму по меньшей мере одного геологического тела на основании геологических данных, относящихся к подземному объему;(b) оценивают поле осаждающего потока, из которого образовано по меньшей мере одно геологическое тело, на основании контурной формы;(c) оценивают распределение зерен по размерам внутри по меньшей мере части по меньшей мере одного геологического тела, используя зависимость между оцененным полем потока, из которого образовано по меньшей мере одно геологическое тело, и распределением зерен по размерам;(d) оценивают свойства породы осадочного бассейна в пределах по меньшей мере части подземного объема, основываясь на распределении зерен по размерам. 20. Способ по п.19, в котором дополнительно используют информацию о минералогии и истории захоронения в сочетании с распределением зерен по размерам для оценивания свойств породы. 21. Способ по п.19, в котором оцененное поле осаждающего потока имеет пространственную изменчивость по меньшей мере одного из скорости потока, высоты потока, концентрации взвешенного осадка, распределения зерен по размерам во взвешенном осадке и любого сочетания из них.