Способ и система моделирования углеводородсодержащего пласта

1 Область изобретения Настоящее изобретение в общем виде относится к моделированию углеводородсодержащего пласта и, в частности, к способу и системе моделирования углеводородсодержащего пласта при условиях, когда текучую среду закачивают в пласт для вытеснения пластовых углеводородов. Способ по данному изобретению особенно полезен при моделировании воздействия вязких языков и протоков в пласте, когда закачиваемая текучая среда протекает через углеводородсодержащий пласт. Предпосылки создания изобретения При первичном извлечении нефти из подземных нефтеносных пластов или коллекторов обычно можно извлечь только ограниченную часть первичной нефти, присутствующей в коллекторе. По этой причине используют разнообразные дополнительные способы извлечения для улучшения вытеснения нефти из породы коллектора. Данные способы обычно можно классифицировать как термические способы извлечения (такие как операции нагнетания пара), способы заводнения и способы на основе вытеснения газом, которые могут осуществляться либо при условиях закачивания смешивающейся, либо несмешивающейся текучей среды. При операциях закачивания смешивающегося агента нагнетаемую текучую среду или растворитель закачивают в коллектор, чтобы получить однофазный раствор с нефтью по месту так, чтобы затем можно было бы извлечь нефть в виде более высокоподвижной фазы из коллектора. Растворитель типично представляет собой легкий углеводород, такой как сжиженный нефтяной газ (СНГ), углеводородный газ,содержащий относительно высокие концентрации алифатических углеводородов в диапазоне от С 2 до С 6, азот или диоксид углерода. Операции добычи путем закачивания смешивающегося агента обычно выполняют процедурой вытеснения, в которой растворитель закачивают в коллектор через нагнетательную скважину для вытеснения нефти из коллектора в направлении эксплуатационной скважины, из которой добывают нефть. Это дает эффективное вытеснение нефти в областях, через которые протекает растворитель. К сожалению, растворитель часто течет неравномерно через коллектор. Поскольку нагнетаемый в коллектор растворитель типично является значительно менее вязким по сравнению с пластовой нефтью, растворитель часто образует языки и протоки через пласт, оставляя часть пласта не охваченной. Кроме того, образование языков является присущей тенденцией высокоподвижного растворителя течь предпочтительно через более проницаемые части породы или не принимать во внимание силу тяжести в коллекторе. Смешиваемость растворителя с пластовой нефтью также влияет на эффективность ее вы 003418 2 теснения внутри коллектора. Некоторые растворители, такие как СНГ, смешиваются непосредственно с пластовой нефтью во всех пропорциях и полученные в результате смеси остаются одной фазой. Говорят, что такой растворитель является смешивающимся во время первого контакта или "смешивающийся при первом контакте". Другие растворители, используемые для закачивания в пласт смешивающихся с нефтью текучих сред, такие как диоксид углерода или углеводородный газ, образуют две фазы при непосредственном смешивании с пластовой нефтью, поэтому они не являются смешивающимися при первом контакте. Однако при достаточно высоком давлении массоперенос компонентов на месте между пластовой нефтью и растворителем образует вытесняющую фазу с переходной зоной составов текучих сред, которые находятся в диапазоне от состава нефти до состава растворителя, и все составы внутри переходной зоны данной фазы являются близко смешивающимися. Смешиваемость, достигаемую массопереносом компонентов на месте,получаемую в результате повторяющегося контакта нефти и растворителя в ходе течения, называют смешиваемостью при "многократном контакте" или динамической. Давление, требуемое для достижения смешиваемости при многократном контакте, называют "минимальным давлением смешиваемости". Растворители ниже минимального давления смешиваемости, называемые "почти смешивающимися" растворителями, могут вытеснять нефть почти так же, как и смешивающиеся растворители. Прогноз эксплуатационных характеристик при закачивании в коллектор смешивающихся с нефтью агентов требует реалистической модели представления пласта. Численное моделирование моделей коллекторов широко используют в нефтяной промышленности в качестве способа использования компьютера для предсказания эффектов явления вытеснения смешивающимися агентами. В большинстве случаев существует желание промоделироватьпроцессы переноса,происходящие в коллекторе. Типично переносу подвергается масса, энергия, момент или некоторая их комбинация. Используя численное моделирование, можно воспроизвести и наблюдать физическое явление и определить расчетные параметры без действительных лабораторных экспериментов и промысловых испытаний. Моделирование коллектора дает представление о состоянии реального углеводородсодержащего коллектора из состояния численной модели данного коллектора. Цель заключается в достаточно хорошем понимании сложных химических, физических процессов и процессов течения текучей среды, происходящих в коллекторе, чтобы прогнозировать будущее состояние коллектора для максимального увеличения добычи нефти. Моделирование коллектора часто относится к гидродинамике потока внутри кол 3 лектора, но в более широком смысле моделирование коллектора может также относиться к общей нефтяной системе, которая включает коллектор, нагнетательные скважины, эксплуатационные скважины, трубопроводы на поверхности и наземное перерабатывающее оборудование. Принципы численного моделирования заключаются в численном решении уравнений,описывающих физическое явление с помощью компьютера. Такие уравнения, как правило,представляют собой обычные дифференциальные уравнения и дифференциальные уравнения в частных производных. Данные уравнения обычно решают, используя численные способы,такие как метод конечного элемента, метод конечной разности, метод конечного объема и аналогичные. В каждом из данных методов физическую систему, которую необходимо моделировать, разделяют на более мелкие элементы сетки или блоки (ряд которых называют сеткой или меш), и переменные состояния, непрерывно изменяющиеся в каждом элементе сетки, представлены набором значений для каждого элемента сетки. В методе конечной разности исходное дифференциальное уравнение заменяют рядом алгебраических уравнений для выражения фундаментальных принципов сохранения массы, энергии и/или момента внутри каждого элемента сетки и переноса массы, энергии и/или момента между элементами сетки. Данные уравнения могут исчисляться миллионами. Такая замена непрерывно изменяющихся значений ограниченным числом значений для каждого элемента сетки называют "дискретизацией". Чтобы проанализировать изменяющееся во времени явление, необходимо вычислить физические количества при дискретных интервалах времени, называемых временным шагом, независимо от непрерывно изменяющихся условий как функции времени. Моделирование зависящих от времени процессов переноса происходит в последовательности временных шагов. При типичном моделировании коллектора решение первичных неизвестных величин,обычно давления, фазовых насыщенностей и составов ищут в специфических точках в интересующей области. Такие точки называют "узлами сетки" или, в более общем случае, "узлами". Элементы сетки создаются вокруг таких узлов и сетку определяют как группу таких элементов сетки. Принимается, что свойства, такие как пористость и проницаемость, являются постоянными внутри элемента сетки. Другие переменные, такие как давление и фазовые насыщенности, определяют в узлах. Звено между двумя узлами называют "связью". Поток текучей среды между двумя узлами типично моделируют в виде потока вдоль связи между ними. Моделирование состава углеводородсодержащего коллектора необходимо для прогнозирования процессов, таких как нагнетание в 4 пласт газа, смешивающегося с нефтью при первом контакте, смешивающегося при многократном контакте и почти смешивающегося. Фазы нефти и газа представлены многокомпонентными смесями. При таком моделировании гетерогенность коллектора и образование вязких языков и протоков вызывают изменения в фазовых насыщенностях и составах в масштабе настолько малом, как несколько сантиметров или менее. Мелкомасштабная модель может представить детали данного поведения вытеснения с неблагоприятной подвижностью. Однако использование мелкомасштабных моделей для моделирования данных изменений, как правило,не является практичным, поскольку мелкий масштаб деталей накладывает запрещающие требования на компьютерные ресурсы. Поэтому, для моделирования нефтеносных пластов обычно разрабатывают грубые модели, имеющие значительно меньшее количество элементов сетки. Значительные усилия были нацелены на разработку моделей, подходящих для использования при прогнозировании поведения пласта при закачивании смешивающегося с нефтью агента. Разработка модели с крупной сеткой, которая эффективно моделирует процессы вытеснения газом, является особенно интересной. Для композиционного моделирования высококачественная модель с крупной сеткой должна эффективно характеризовать изменения в фазовом поведении и изменения в составах нефти и газа при протекании вытеснения нефти. Предлагались многие различные способы. Большинство из данных предложений используют эмпирические модели для представления образования вязких языков при вытеснении агентом, смешивающимся при первом контакте. См., например:Validation of an Empirical Model for Viscous Fingering with Gravity Effects," SPE Reservoir Engineering, стр. 542-550, May 1988. Из данных моделей модель смешения Тодда-Лонгстаффа ("Т-Л") является наиболее популярной, и она широко используется при моде 5 лировании коллекторов. При использовании надлежащим образом модель смешения Т-Л предоставляет достаточно точные средние характеристики вытеснения при неблагоприятной подвижности, когда нагнетаемый растворитель и нефть являются смешивающимися при первом контакте. Однако модель смешения Т-Л является менее точной в условиях смешения при многократном контакте. Были предложены модели, которые используют Т-Л модель для объяснения образования вязких языков при ситуациях вытеснения агентом, смешивающимся при многократном контакте (см., например, Todd, M.R. and Chase,С.A., "A Numerical Simulator for PredictingEngineers, Houston, Texas, 1979, иногда относимую к методике "Тодда-Чаза"). При моделировании вытеснения агентом, смешивающимся при многократном контакте, в добавление к образованию вязких языков, принимаемому во внимание в модели смешения Т-Л, также необходимо рассматривать обмен компонентов растворителя и нефти между фазами согласно соотношениям фазового поведения. Важность взаимодействия между фазовым поведением и образованием языков при вытеснении агентами,смешивающимися при многократном контакте,описывается в работе Gardner, J.W., and Ypma,J.G.J., "An Investigation of Phase Behavior/Macroscopic Bypassing Interaction in СО 2Flooding," Society of Petroleum Engineering Journal, стр. 508-520, Октябрь 1984. Однако, данные предложения не могут эффективно объединить модель смешения и модель фазового поведения. Другая предложенная модель для принятия во внимание образования языков и протоков при вытеснении агентом, смешивающимся при многократном контакте, предлагает делать дисперсность компонентов растворителя и нефти зависимой от градиента вязкости, таким образом обращаясь к макроскопическим эффектам образования вязких языков (см. Young, L.C.,"The Use of Dispersion Relationships to ModelSPE/DOE Enhanced Oil Recovery Symposium,Tulsa, Апрель 20-23). В другой модели было предложено расширить Т-Л модель до многофазного многокомпонентного потока с упрощенным прогнозом фазового поведения (см.Symposium, Tulsa, Апрель 17-20). В еще одной модели было предложено использовать составы текучей среды, вытекающей из элемента сетки большого размера, для компенсации неодно 003418 6 родности распределения текучей среды в элементе сетки (см. Barker, J.W., and Fayers, F.J.,"Transport Coefficients for Compositional Simulation with Coarse Grids in Heterogeneous Media",SPE 22591, представленную на SPE 66th AnnualTech. Conf., Dallas, TX, Октябрь 6-9, 1991). В еще одной модели было предположено, что неполное смешение между растворителем и нефтью можно представить, предполагая, что термодинамическое равновесие превалирует только на границе раздела между двумя фазами, и процесс диффузии ведет составы нефти и растворителя по направлению к данным равновесным значениям (см. работу Nghiem, L.X., and Sammon, P.H., "A Non-Equilibrium Equation-of-StateCompositional Simulator," SPE 37980, представленную на 1997 SPE Reservoir Simulation Symposium, Dallas, TX, Июнь 8-17, 1997). Элементы сетки в данных моделях не были подразделены. Были сделаны предложения представить образование языков и протоков при вытеснении агентом, смешивающимся при многократном контакте, с использованием двухзонных моделей. См., например:Fayers, F.J., Barker, J.W., and Newley,T.M.J., "Effects of Heterogeneities on Phase Behavior in Enhanced Oil Recovery", in The Mathematics of Oil Recovery", P.R. King, editor, стр. 115-150, Clarendon Press, Oxford, 1992. Данные модели разделяют элементы сетки моделирования на зону, где имеет место полное смешение между закаченным растворителем и частью пластовой нефти, и зону, где пластовую нефть обходят и она не контактирует с растворителем. Хотя оказывается, что концептуальная структура данных моделей предлагает лучшее представление неполного смешения в процессах вытеснения агентом, смешивающимся при многократном контакте, чем модели одной зоны,физическая основа уравнений, используемых для описания просачивания и смешения, не является ясной. В частности, данные модели используют эмпирические корреляции для представления подвижностей нефти/растворителя в каждой зоне, используют эмпирические корреляции для представления переноса компонентов между зонами и создает ограничивающие допущения по составу зон и направлению переноса компонентов между зонами. Предполагалось,что эмпирическую подвижность и функции массопереноса в данных моделях можно определить, подгоняя их к результатам моделирования с мелкой сеткой. В результате на практике калибровка данных моделей будет отнимать много времени и являться дорогим процессом. Более того, маловероятно, что данные модели точ 7 но прогнозируют состояние вне диапазонов параметров, исследуемых в указанном моделировании с мелкой сеткой. В то время, как предлагаемые в прошлом двухзонные подходы имеют определенные преимущества, продолжает существовать необходимость в улучшенных моделях, которые предложили бы лучшее физическое представление просачивания и смешения при вытеснении с неблагоприятной подвижностью и, таким образом, способных предоставить более точный и эффективный прогноз состояния при закачке вытесняющего агента. Сущность изобретения Предлагается способ и система моделирования одной или нескольких характеристик многокомпонентного, углеводород-содержащего пласта, в который для вытеснения пластовых углеводородов нагнетают вытесняющую текучую среду, имеющую, по меньшей мере, один компонент. Первой стадией способа является представление, по меньшей мере, части пласта в виде множества элементов сетки. Каждый элемент сетки затем делят на две зоны, причем первая зона представляет часть каждого элемента сетки, охваченную вытесняющей текучей средой, а вторая зона представляет часть каждого элемента сетки, по существу, не охваченную вытесняющей текучей средой. Предполагается,что распределение компонентов в каждой зоне является, по существу, однородным. Создана модель, которая представляет свойства текучей среды внутри каждой зоны, потока текучей среды между элементами сетки, используя принципы теории фильтрации, и перенос компонентов между областями. Далее модель используется в моделирующей программе для моделирования одной или нескольких характеристик пласта. Краткое описание чертежей Настоящее изобретение и его преимущества будут лучше поняты из нижеследующего детального описания прилагаемых чертежей, на которых изображено следующее. Фиг. 1 иллюстрирует двухмерную схему течения растворителя через нефтяной коллектор для вытеснения из него нефти и показывает пример образования языков растворителя в коллекторе; фиг. 2 - пример двухмерной мелкомасштабной сетки, которая может представлять область коллектора из фиг. 1; фиг. 3 - двухмерный элемент сетки, охватывающий ту же проблемную область, изображенную на фиг. 1, с разделенным на две зоны элементом сетки, причем одна зона представляет зону области, охваченную закаченной текучей средой, вторая зона представляет зону области, не охваченную закаченной текучей средой. фиг. 4 иллюстрирует элемент сетки, изображенный на фиг. 3, схематически показывая фазовые доли в двух зонах элемента сетки; 8 фиг. 5 А - влияние координационного числа, z, на общую добычу нефти при нагнетании агента, смешивающегося при многократном контакте, моделируемое с использованием способа по настоящему изобретению; фиг. 5 В - влияние координационного числа, z, на прорыв фронта растворителя при нагнетании агента, смешивающегося при многократном контакте, моделируемое с использованием способа по настоящему изобретению; фиг. 6 А-D - влияние чисел Damkhler нефти на кривые добычи тяжелой нефти и легкой нефти при нагнетании агента, смешивающегося при многократном контакте, моделируемое с использованием способа по настоящему изобретению; фиг. 7 графически сравнивает опубликованные данные по добыче нагнетанием агента,смешивающегося при первом контакте, и наилучшие приближения, полученные с использованием способа по настоящему изобретению; фиг. 8 иллюстрирует координационные числа, полученные приведением в соответствие модели, используемой в способе по настоящему изобретению, и опубликованных данных как функцию отношения вязкости нефти/растворителя; фиг. 9 - опубликованные экспериментальные данные по добыче CO2/Soltrol и СO2/нефтьWasson, полученные на кернах, и модельные оценки, использующие опубликованную модель одной зоны; фиг. 10 - опубликованные экспериментальные данные по извлечению CO2/Soltrol иCO2/нефть Wasson, полученные на кернах, и модельные оценки, использующие способ по настоящему изобретению. Фигуры иллюстрируют конкретные варианты осуществления на практике способа по настоящему изобретению. Фигуры не предназначены для исключения из рамок изобретения других вариантов осуществления изобретения,которые являются результатом обычных и ожидаемых модификаций конкретных вариантов осуществления изобретения на практике. Подробное описание изобретения Для более полного понимания настоящего изобретения предоставляются следующие вводные комментарии. Для увеличения добычи углеводородов из подземных пластов были разработаны разнообразные способы увеличения добычи углеводородов, при котором в подземный пласт закачивают текучую среду в одну или несколько нагнетательных скважин месторождения, и углеводороды (а также закачиваемую текучую среду) добывают из пласта с помощью одной или несколько эксплуатационных скважин месторождения. Нагнетательные скважины обычно располагаются отдельно от эксплуатационных скважин, но одна или несколько нагнетательных скважин может позднее использоваться в каче 9 стве эксплуатационных скважин. Закачиваемая текучая среда может являться, например, любым нагревающим агентом, используемым в способе термической добычи (таким как пар),любой, по существу, не смешивающейся текучей средой, используемой в процессе закачивания несмешивающегося агента (такой как природный газ, вода или соляной раствор), или любой смешивающейся текучей средой, используемой в процессе закачивания смешивающегося агента (например, текучая среда, смешивающаяся при первом контакте, такая как сжиженный нефтяной газ, или текучая среда, смешивающаяся при многократном контакте или почти смешивающаяся текучая среда, такая как углеводороды с более низкими молекулярными массами, диоксид углерода или азот). Фиг. 1 схематично иллюстрирует двухмерную область 5 коллектора, которая является частью более крупного нефтесодержащего геологического пласта (на чертеже не показан),который необходимо анализировать с использованием способа по данному изобретению. На фиг. 1 закачиваемая текучая среда 11, которая,как принимают, является газообразной в данном описании, вытесняет многокомпонентную пластовую нефть 12 в области 5 коллектора. Необходимо понимать, что настоящее изобретение не ограничивается газообразной закачиваемой текучей средой, закачиваемая текучая среда также может быть жидкостью или многофазной смесью. Закачиваемая текучая среда 11 течет слева направо на фигуре. Фиг. 1 изображает образование вязких языков, которое происходит, когда закачиваемая текучая среда 11 вытесняет пластовую нефть 12. Закачиваемая текучая среда 11 стремится образовать языки через нефть 12 по направлению к эксплуатационной скважине (не показана на фигуре), приводя к преждевременному прорыву закачиваемой текучей среды 11 и обходу некоторого количества пластовой нефти 12. Образование вязких языков в основном вызывается большими различиями в вязкостях нефти 12 и закачиваемой текучей среды 11, приводящими к коэффициенту подвижности закачиваемой текучей среды и нефти, который имеет неблагоприятное действие на эффективность охвата области или на эффективность вытеснения закачиваемой текучей средой. Посредством прогрессивного способа определения характеристики коллектора область 5 коллектора можно представить элементами сетки в масштабе от сантиметров до нескольких метров, иногда называемой мелкомасштабной сеткой. Каждый элемент сетки может быть занят свойствами коллектора, включая, например,тип породы, пористость, проницаемость, начальную внутрипоровую насыщенность текучей средой, функции относительной проницаемости и капиллярного давления. 10 Фиг. 2 показывает пример двухмерной мелкомасштабной сетки 10, которая может представлять область 5 коллектора из фиг. 1. Область коллектора 5 из фиг. 1 представлена на фиг. 2 84 элементами сетки. Элементы 11' сетки представляют геологические зоны, которые охвачены закаченной текучей средой 11, а элементы 12' сетки представляют геологические зоны,которые содержат, по существу, пластовую нефть 12, не вытесненную закаченной текучей средой. Однако моделирование коллектора обычно не проводят с помощью мелкомасштабных сеток. Прямое использование мелкомасштабных моделей для полномасштабного моделирования коллектора обычно не осуществимо,поскольку мелкий масштаб деталей накладывает запрещающие ограничения на компьютерные ресурсы. Поэтому при моделировании обычно используют сетки с крупным масштабом, в то же время сохраняя, насколько это возможно,характеристики потока текучей среды и фазовое состояние мелкомасштабной модели. Сетка с грубым масштабом может представлять, например, все 84 элемента сетки на фиг. 2 одним элементом сетки. Следовательно, необходим способ моделирования составов текучей среды и поведения потоков текучей среды, принимающий во внимание образование языков и протоков. Способ по данному изобретению предоставляет такую возможность. Способ по настоящему изобретению начинает с представления области коллектора, которую необходимо анализировать, в виде соответствующей сетчатой системы. Область коллектора, которую необходимо анализировать, представлена множеством элементов сетки, расположенных по соседству друг с другом так, чтобы иметь границу между каждой парой соседних элементов сетки. Данную пространственную дискретизацию области коллектора можно провести, используя способы конечной разности, конечного объема, конечного элемента или аналогичные хорошо известные способы, которые основаны на разделении физической системы, которую необходимо моделировать, на меньшие единицы. Настоящее изобретение прежде всего описывается касательно использования способа конечной разности. Специалисты в данной области легко поймут, что настоящее изобретение также можно применять в связи со способами конечного элемента и способами конечного объема. При использовании способов конечной разности и конечного объема, единицы меньшего размера типично называют элементами сетки, а при использовании способа конечного элемента единицы меньшего размера типично называют элементами. Данные элементы сетки или элементы могут насчитывать от менее чем сотни до миллионов. В данной заявке для простоты представления используют термин элемент сетки, но необходимо понимать, что если при моделировании используется способ 11 конечного элемента, термин элемент заменит термин элемент сетки, когда он используется в данном описании. При реализации на практике настоящего изобретения, элементы сетки могут быть любой геометрической формы, такой как параллелепипеды (или кубы) или шестигранники (имеющие четыре вертикальные угловые грани, которые могут различаться по длине) или четырехгранник, ромбоид, трапецоид или треугольник. Сетка может включать прямоугольные элементы сетки, организованные в правильную, структурированную форму (как иллюстрируется на фиг.2), или она может включать элементы сетки, имеющие различные формы, расположенные в неправильной, неструктурированной форме,или она может включать многообразие обоих типов структурированных и неструктурированных форм. Полностью неструктурированные сетки могут быть составлены так, чтобы принять почти любую форму. Все элементы сетки предпочтительно расположены по границе, посредством чего избегая контакта любой стороны элемента сетки со сторонами двух других элементов сетки. Один из типов гибкой сетки, которую можно использовать в модели по изобретению,представляет сетку Voronoi. Элемент сетки Voronoi характеризуют как область пространства,которая ближе к своему узлу, чем к любому другому узлу, и сетку Voronoi создают из таких элементов сетки. Каждый элемент сетки связан с узлом и серией соседних элементов сетки. Сетка Voronoi является локально ортогональной в геометрическом смысле, т.е. границы элемента сетки являются перпендикулярными по отношению к линиям, соединяющим узлы с двух сторон каждой границы. По этой причине сеткиVoronoi также называют сеткой с перпендикулярным делением пополам (ПДП). Прямоугольный сетчатый блок (сетка Cartesian) является особым случаем сетки Voronoi. Сетка ПДП имеет гибкость, чтобы представить изменяющуюся в широких пределах геометрию коллектора,поскольку положение узлов можно выбирать свободно. Сетки ПДП получают, устанавливая положение узлов в данной области и затем задавая границы элемента сетки таким образом, что каждый элемент сетки содержит все точки, которые ближе к положению его узла, чем к положению любого другого узла. Поскольку межузловые соединения в сетке ПДП являются перпендикулярно разделенными пополам границами элемента сетки, это значительно упрощает решение уравнений потоков. Для более детального описания создания сетки ПДП см. работу Palagi, C.L. and Aziz, К.: "Use of Voronoi 12 Следующим шагом в способе по настоящему изобретению является разделение каждого элемента сетки, который охвачен закаченной текучей средой, на две зоны, первую зону, которая представляет часть элемента сетки, охваченную закаченной текучей средой 11, и вторую зону, которая представляет часть элемента сетки, не охваченную закаченной текучей средой 11. Принимают, что распределение компонентов в каждой зоне является однородным. Далее принимают, что текучие среды внутри каждой зоны находятся в термодинамическом равновесии. Однако, две зоны элемента сетки не находятся в равновесии друг с другом, и, как результат, составы и доли фазовых объемов внутри каждой зоны типично будут различаться. Фиг. 3 иллюстрирует двухмерное схематическое решение одного элемента 15 сетки, который представляет ту же область коллектора,представленную 84 элементами сетки сети 10(фиг. 2). В то время, как это не показано на фигурах, необходимо понимать, что элемент 15 сетки делит границы с соседними элементами сетки. Следующее ниже описание, касающееся элемента 15 сетки, также приложимо к другим элементам сетки, в которой элемент сетки является только одним из множества элементов сетки. Как показано на фиг. 3, элемент 15 сетки разделяют на две зоны 16, 17. Зона 16 представляет часть элемента сетки, охваченную закаченной текучей средой 11, а зона 17 представляет часть элемента сетки, которая не была вытеснена закаченной текучей средой 11. Зоны 16, 17 разделены поверхностью раздела или разделом 18, который, как принимают, имеет бесконечно малую толщину. Принимают, что многокомпонентные текучие среды внутри каждой зоны находятся в термодинамическом равновесии,которое означает, что составы текучих сред и фазовые объемы зон 16, 17 могут быть различными и обычно являются различными. Составы текучих сред могут различаться от элемента сетки к элементу сетки внутри сетки, и составы текучих сред внутри каждой зоны могут изменяться со временем. Поэтому раздел 18 может двигаться как функция времени, когда закачиваемая текучая среда 11 контактирует с большей зоной, представленной элементом 15 сетки. Движение раздела 18 зависит прежде всего от обмена текучих сред между элементом 15 сетки и соседними с ним элементами сетки, массопереноса через раздел 18 и закачивания или извлечения текучих сред через нагнетательную или эксплуатационную скважины, которые могут проходить через геологическую область, представленную элементом сетки. Фиг.4 иллюстрирует пример фазовых долей текучих сред в зонах 16, 17. Доля паровой фазы, которая состоит из закаченной текучей среды плюс испарившаяся нефть, показана цифрой 11 а в зоне 16 и цифрой 11b в зоне 17. Доля 13 жидкой фазы, которая состоит из пластовой нефти плюс растворенная закаченная текучая среда, показана цифрой 12 а в зоне 16 и цифрой 12b в зоне 17. Доля воды показана цифрой 13 а в зоне 16 и цифрой 13b в зоне 17. В примере, показанном на фиг. 4, зона 16 содержит прежде всего высокоподвижную закаченную текучую среду 11, а зона 17 содержит прежде всего малоподвижную пластовую нефть 12. Стрелка 20 представляет поток текучей среды, втекающий в зону 16 из охваченных зон элементов сетки,соседних с элементом 15 сетки. Стрелка 21 представляет поток текучей среды, втекающий в зону 17 из резидентных зон элементов сетки,соседних с элементом 15 сетки. Стрелка 22 представляет поток текучей среды, вытекающий из зоны 16 в охваченные зоны элементов сетки,соседних с элементом 15 сетки. Стрелка 23 представляет поток текучей среды, вытекающий из зоны 17 в резидентные зоны элементов сетки,соседних с элементом 15 сетки. Хотя стрелки показывают текучую среду, перемещающуюся слева направо, текучая среда могла бы втекать и вытекать из элемента 15 сетки в других направлениях. Стрелки 24 представляют массоперенос между зонами 16, 17. Компонентам дают возможность переноситься в любом направлении через раздел 18. Хотя стрелки 24 показывают перенос через между фазами одного типа (из паровой в паровую, из жидкого углеводорода в жидкий углеводород и из воды в воду), компоненты могут переноситься из любой фазы в зоне, являющеейся источником, в любую фазу в другой зоне. Зона 16 имеет нулевой объем, пока закачиваемая текучая среда не втечет в элемент 15 сетки. Можно моделировать закачиваемую текучую среду 11 в виде закачиваемой либо в зону 16, охваченную процессом вытеснения,либо в резидентную зону 17, или можно моделировать закачиваемую текучую среду 11 в виде закачиваемой в обе зоны 16 и 17. Текучие среды могут извлекаться как из зоны 16, охваченной процессом вытеснения, так и из резидентной зоны 17. Также можно моделировать элемент 15 сетки в таком виде, что закачиваемая текучая среда 11 течет от одной или нескольких нагнетательных скважин непосредственно в элемент 15 сетки, и его можно моделировать в таком виде, что закачиваемая текучая среда 11 непосредственно вытекает из элемента 15 сетки в одну или несколько эксплуатационных скважин. Хотя это не показано на чертежах, если через область коллектора, представленную элементом 15 сетки, проходит нагнетательная скважина,вытесняющую текучую среду 11, закачиваемую в элемент 15 сетки, можно моделировать в виде текучей среды, закачиваемой только в зону 16,охваченную процессом вытеснения, и если через область коллектора, представленную элементом 15 сетки, проходит эксплуатационная скважина, элемент 15 сетки можно моделировать, имеющим текучие среды, получаемые как 14 из зоны, охваченной процессом вытеснения, так и резидентной зоны 17. Хотя чертежи не показывают узлы элементов сетки, специалисты в данной области поймут, что каждый элемент сетки будет иметь узел. Принимается, что в операциях моделирования поток текучей среды между элементами сетки происходит между узлами элементов сетки или, говоря другими словами, через межузловые соединения. При реализации на практике настоящего изобретения, охваченную процессом вытеснения зону данного элемента сетки(зона 16 на фиг. 3 и 4) соединяют с охваченными процессом вытеснения зонами элементов сетки, соседними с данным элементом сетки, а резидентные зоны данного элемента сетки (зона 17 на фиг. 2) соединяют с резидентными зонами элементов сетки, соседних с данным элементом сетки. Между зоной 16, охваченной процессом вытеснения, и резидентной зоной 17 нет межузловых соединений. Поэтому заявители иногда относят способ настоящего изобретения к модели с разделенными узлами (МРУ). Следующим шагом в способе по настоящему изобретению является создание прогнозирующей модели, которая представляет свойства текучих сред внутри каждой зоны каждого элемента сетки, поток текучей среды между каждым элементом сетки и соседних с ним элементов сетки, и перенос компонентов между зонами 16, 17 для каждого элемента сетки. В предпочтительном варианте осуществления на практике модель включает ряд конечно-разностных уравнений для каждого элемента сетки, имеющих функции, представляющие подвижность каждой фазы текучей среды в зонах 16, 17, функции,представляющие фазовое состояние внутри зон 16, 17, и функции, представляющие массоперенос каждого компонента между зонами 16, 17. Необязательно модель может далее содержать функции, представляющие перенос энергии между зонами 16, 17. Функции переноса энергии могут быть желательны, например, для моделирования тепловых эффектов в результате операций закачки пара. Функции подвижности используют для описания потока через связи, и функция подвижности создается для каждой фазы в каждой зоне. Подвижности потоков 22, 23, уходящих из элемента 15 сетки, зависят от многих факторов,включая состав текучих сред в зоне 16, охваченной процессом вытеснения, и резидентной зоне 17, относительный размер (или объемную долю) зоны 16, охваченной процессом вытеснения, и резидентной зоны 17, гетерогенность элемента сетки и коэффициент подвижности нефти/закаченной текучей среды. Конкретные функциональные зависимости определяют посредством использования теории фильтрации. Основные принципы теории фильтрации описываются в работе S. Kirkpatrick, "Percolationand Conduction," Rev. Modern. Physics, Vol. 45, 15 стр. 574-588, 1973, которая включена здесь в качестве ссылки. При предпочтительном варианте осуществления на практике модель подвижности эффективной среды представляет элемент сетки в виде пористой сетки так, чтобы характеризовать влияние образования языков и протоков, которые могут иметь место в элементе сетки в зависимости от условий, преобладающих в элементе сетки в течение временного интервала. Эффективная подвижность каждой жидкой фазы в каждой зоне элемента сетки может быть вычислена специалистом в данной области с помощью данного описания. Примеры уравнений фазовой подвижности, полученные из модели эффективной среды, предоставляются ниже в виде уравнений (18)-(20). Способ по настоящему изобретению предполагает, что равновесие существует внутри зоны 16, охваченной процессом вытеснения, и внутри резидентной зоны 17. В качестве части модели выполняется определение свойств фаз,которые сосуществуют внутри зон 16 и 17. Предпочтительно, используют соответствующее уравнение состояния для расчета фазового поведения зоны 16, зоны 17. В примерах, предоставленных ниже, одномерная модель использует упрощенную модель состояния псевдотроичной фазы, которая характеризует смеси растворителя и нефти с точки зрения трех псевдокомпонент, растворителя (CO2), легкого компонента нефти и тяжелого компонента нефти. Упрощенная модель фазового состояния способна моделировать характерные черты процесса вытеснения, включая различные степени смешиваемости, находящиеся в диапазоне от смешивающегося при первом контакте, через смешивающийся при многократном контакте и почти смешивающегося до несмешивающегося агентов. Свойства фазового состояния могут быть определены специалистом в данной области. Способ по настоящему изобретению не предполагает равновесие между зоной 16, охваченной процессом вытеснения, и резидентной зоной 17 элемента сетки. Функции массопереноса используют для описания скорости движения компонентов через поверхность раздела или раздел 18 между зонами 16, 17. Этот массоперенос изображен на фиг.4 стрелками 24. Механизм массопереноса включает, но не ограничивается ими, молекулярную диффузию, конвективную дисперсию и капиллярную дисперсию. Способ по изобретению предполагает, что скорость каждого компонента массопереноса пропорциональна времени сопротивления движущей силе. Примеры движущих сил включают,но не ограничиваются ими, различие в составе и разница в капиллярном давлении между двумя зонами. Как только создается функция массопереноса для каждого компонента текучей среды,скорости массопереноса зависят от факторов,включая, но не ограничиваясь ими, идентичность компонентов, степень смешиваемости 16 между газом и нефтью, размер каждой зоны,геометрию элемента сетки, отношение подвижности газа/нефти, скорость, гетерогенность и водонасыщенность. Данные функциональной зависимости могут быть построены в модель массопереноса специалистами в данной области. Примеры функций массопереноса предоставляются в виде уравнений (10) и (14)-(16) ниже. Одним из первых шагов при разработке модели является выбор числа размерности пространства, желательной для представления геометрии коллектора. Необходимо рассматривать как внешнюю, так и внутреннюю геометрии. Внешние геометрии включают пределы коллектора или водоносного пласта (или элемент симметрии) и верхнюю и нижнюю часть коллектора или водоносного пласта (включая сбросы). Внутренняя геометрия включает степень индивидуальной проницаемости по площади и по вертикали единиц и непродуктивных зон, которые являются важными для решения проблемы и определения геометрии скважины (например,диаметра скважины, расстояния между скважинами и присутствия гидравлических разрывов,происходящих из скважины). Модель по настоящему изобретению не ограничивается конкретным числом размеров. Может быть создана прогнозирующая модель для одномерного (1-D), двухмерного (2-D) и трехмерного (3-D) моделирования коллектора. Модель 1-D будет редко использоваться для широкомасштабного изучения коллектора, поскольку она не может моделировать охват по площади и по вертикали. Модель 1-D нагнетания газа для прогнозирования эффективности вытеснения не может эффективно представлять гравитационные эффекты, перпендикулярные направлению потока. Однако, модели 1-D нагнетания газа можно использовать для исследования чувствительности характеристик коллектора к изменениям в параметрах процесса и для интерпретации лабораторных тестов по вытеснению. Ареальные модели 2-D нагнетания текучей среды можно использовать, когда структура потока по площади имеет преобладающее влияние на характеристики коллектора. Например,ареальные модели нормально используются для сравнения возможных сеток размещения скважин или для оценки влияния гетерогенности области на состояние коллектора, поперечные или радиальные модели 2-D закачивания газа можно использовать, когда структура потоков в вертикальном разрезе имеет преобладающее влияние на состояние коллектора. Например,поперечные или радиальные модели обычно используют для моделирования процессов, на которые преобладающее влияние оказывает сила тяжести, такие как закачивание газа в сводную часть пласта или закачивание газа в коллекторы, имеющие высокую вертикальную про 17 ницаемость, и для оценки влияния вертикальной гетерогенности на состояние коллектора. Модели 3-D могут быть желательны для эффективного представления сложной геометрии коллектора или сложной механики текучих сред в коллекторе. Например, модель может представлять собой модель 3-D, включающую слои сеток ПДП, которая иногда называется в нефтедобывающей промышленности как 2,5-Д. Слоистые сетки ПДП являются неструктурированными по поверхности и структурированными (слоистыми) по вертикали. Конструкция слоистых сеток 3-D описывается в Heinemann,Z.E., et al., "Modeling Reservoir Geometry WithScheme for Flexible Grids in Reservoir Simulation," SPE 37999, SPE Reservoir Simulation Symposium, Dallas, TX, Июль, 1997. Настоящее изобретение не ограничивается делением элемента сетки только на две зоны. Способ по настоящему изобретению можно использовать с элементами сетки, имеющими множество частей, таким образом деля элементы сетки на три или более зоны. Например,трехзонный элемент сетки может иметь одну зону, представляющую зону коллектора, охваченную закаченной текучей средой, вторую зону, представляющую зону коллектора, не охваченную закаченной текучей средой, и третью зону, представляющую смешанную зону пластовой текучей среды и закаченной текучей среды коллектора. В другом примере при операции нагнетания пара одна зона может представлять зону коллектора, охваченную закаченным паром, вторая зона, может представлять зону коллектора, занятую газом, отличным от пара, и третья зона может представлять зону коллектора, не занятую закаченным паром или другим газом. Газ, отличный от пара, может представлять собой природный газ, растворенный в нефти коллектора, который выделился из пластовой нефти, когда давление коллектора упало ниже давления насыщения нефти, или второй закаченный газ, такой как обогащенный газ, легкий углеводородный газ или СO2. Способ по настоящему изобретению можно использовать для моделирования добычи нефти из коллекторов с вязкой нефтью, при которой тепловая энергия вводится в коллектор для нагрева нефти, посредством чего уменьшается ее вязкость до точки, когда нефть можно заставить течь. Тепловая энергия может находиться в различных формах, включая нагнетание горячей воды и закачивание пара. Закачивание можно проводить в одну или несколько нагнетательных скважин, а добычу нефти можно осуществлять из одной или нескольких находящихся на расстоянии друг от друга эксплуатационных скважин. Одну скважину также можно использовать как для закачивания текучей среды, так и для добычи нефти. 18 Например, в процессе "huff and puff" пар вводят через скважину (которая может являться вертикальной или горизонтальной скважиной) в месторождение вязких углеводородов в течение некоторого периода времени, скважину останавливают для нагрева паром углеводородов,затем скважину переводят на добычу. Как только прогнозирующая модель создана, ее можно использовать в моделирующей программе для моделирования одной или нескольких характеристик пласта как функции времени. Модель основных потоков состоит из уравнений, которые управляют неустановившимся потоком текучих сред в ячеистой сетке коллектора, скважинах и наземных объектах. Для решения уравнений основных потоков специалистами в данной области могут быть выбраны соответствующие численные алгоритмы. Примеры численных алгоритмов, которые можно использовать, описываются в ReservoirSimulation, Henry L. Doherty Series Monograph,Vol.13, Mattax, C.C. and Dalton, R.L., editors, Society of Petroleum Engineers, Richardson, TX,1990. Моделирующая программа представляет собой систему компьютерных программ, которые реализуют численный алгоритм на компьютере. Специалистам в данной области ясно, что при осуществлении на практике настоящего изобретения интенсивно используется компьютер. Соответственно, использование компьютера, предпочтительно цифрового компьютера,для осуществления на практике настоящего изобретения является фактической необходимостью. Компьютерные программы для различных частей способа моделирования имеются в продаже (см., например, программные продукты,продаваемые для разработки элементов сетки,показа результатов, вычисления свойств потока текучей среды и решения линейного ряда уравнений, которые используются в моделирующей программе). Компьютерные программы для других частей изобретения могут быть разработаны специалистами в данной области на основе сформулированного здесь описания. Практика настоящего изобретения может применяться к части или ко всем элементам сетки моделируемой сетчатой системы. Для экономии компьютерного времени дополнительные вычисления, связанные с разделением элементов сетки на две или более зоны, предпочтительно применяют только к тем моделируемым элементам сетки, которые охвачены вытесняющей жидкостью. Способ по настоящему изобретению имеет преимущество по сравнению с двухмерными моделями вытеснения, которые использовались ранее. Данное преимущество может быть приписано следующим ключевым отличиям. Первое, для характеристики воздействия образования языков и протоков на эффективные подвижности жидкости используется теория фильт 19 рации. Второе, скорость переноса компонента между областями пропорциональна времени сопротивления движущей силе. Третье, функции массопереноса объясняют действительные процессы смешения, такие как молекулярная диффузия, конвективная дисперсия и капиллярная дисперсия. Данные улучшения приводят к более точному и эффективному прогнозированию вытеснения при неблагоприятной подвижности. Примеры одномерного моделирования Создана одномерная модель по настоящему изобретению и данную модель испытали,используя собственную моделирующую программу. Имеющиеся коммерческие моделирующие программы могут быть легко модифицированы специалистами в данной области, используя описание данного изобретения и представленные здесь допущения, чтобы получить результаты, по существу, аналогичные тем, что представлены ниже. В модели распределение компонентов между резидентной и охваченной процессом вытеснения зонами определяют уравнениями переноса, которые объясняют конвекцию вытесняющей и пластовой текучих сред и скорость переноса каждого компонента между зонами. В моделирующей программе используют описание четырехкомпонентной текучей среды. Четырьмя компонентами являются растворитель (CO2), легкая фракция сырой нефти,тяжелая фракция сырой нефти и вода. Предполагается, что текучие среды являются несжимаемыми и имеет место идеальное смешивание,что позволяет разъединить уравнения давления и уравнения переноса компонентов и заменить объемные доли на мольные доли в качестве переменных состава. Специалисты в данной области хорошо знакомы со способами расчета сжимаемости текучих сред и неидеального смешения. Также предполагается, что растворитель не переносится в резидентную зону и что водонасыщенность является одной и той же в обеих зонах. Следующее ниже описание примеров моделирования адресует к уравнениям, имеющим значительное количество математических символов, многие из которых определяются, когда они встречаются в тексте. Кроме того, для цели полноты описания, представлена таблица, следующая за детальным описанием и содержащая обозначения используемых здесь символов. Моделирующая программа создана в терминах стандартных уравнений переноса для общего количества каждого компонента, к которым добавлены уравнения переноса для количества каждого компонента в резидентной зоне. Количество каждого компонента в зоне, охваченной процессом вытеснения, затем получают вычислением разности. При данных предположениях безразмерные уравнения переноса для общего растворителя, тяжелого компонента Общую объемную долю легкого компонента w3 получают изw3 = 1- w1 - w2 - Sw В уравнении (4) компонент 1 является растворителем, компонент 2 является тяжелой фракцией нефти, а компонент 3 является легкой фракцией нефти. В уравнениях (1)-(4)x/L,ut/L,k/uL, tive + ile + roe + w, L представляет собой длину керна, k является проницаемостью, является пористостью, Рc является капиллярным давлением между нефтью и водой, yj представляет собой объемную долю компонента j в паровой части зоны, охваченной процессом вытеснения, xj представляет собой объемную долю компонента j в жидкой части зоны, охваченной процессом вытеснения, и xrj представляет собой объемную долю компонента j в неводной части резидентной зоны. wjwrj + wij представляют общие объемные доли компонента j, где wij(SgYj + Slxj) представляет объемную долю компонента j в зоне, охваченной процессом вытеснения, и wrj(1-)(1 - Sw)xrj представляет объемную долю компонента j в резидентной зоне.представляет объемную долю зоны, охваченной процессом вытеснения, определяемую какSg и Sl, соответственно, представляют собой насыщенности паром и жидкостью в зоне,охваченной процессом вытеснения. rое представляет подвижность пластовой текучей среды,ive представляет подвижность паровой фазы в зоне, охваченной процессом вытеснения, ilе представляет подвижность жидкой фазы в зоне,охваченной процессом вытеснения, и w, представляет подвижность воды, причем все параметры вычислены с использованием теории эффективной среды, как описывается ниже. Предполагается, что общая скорость закачивания является постоянной. Безразмерные уравнения транспорта для пластового растворителя, тяжелой нефти и легкой нефти, соответственно, представляют собой(объем/время) компонента j из резидентной зоны в зону, охваченную процессом вытеснения. Первый член с правой стороны данных уравнений объясняет конвекцию каждого компонента внутри резидентной зоны, а второй член объясняет перенос каждого компонента из резидентной зоны в зону, охваченную процессом вытеснения. Уравнение для давления представляет собой В одномерной моделирующей программе уравнения (1)-(3) и (6)-(8) подвергают дискретизации, чтобы получить шесть рядов конечноразностных уравнений в , которые решают разумно по времени способом прогноза и коррекцииHamming, интегрируя ряд обычных дифференциальных уравнений первого порядка (методHamming хорошо известен специалисту в данной области). Предполагается, что до закачивания растворителя не присутствует никакой зоны, охваченной процессом вытеснения, и, что,следовательно,вначале равно нулю в модели. Образование зоны, охваченной процессом вытеснения, инициируется, предполагая, что растворитель уходит исключительно в зону, охваченную процессом вытеснения, у поверхности закачивания в керн. После того, как wi, wri и Swij вычисляют вышеуказанным интегрированием,корректируют с помощью уравнения (5) и интегрирование продолжают далее к следующему временному шагу. Затем определяют распределение давления на каждом временном шаге, интегрируя уравнение (9) относительно . Функция массопереноса Предполагается, что в первом приближении скорость переноса между зонами пропорциональна разнице между объемными долями компонентов в резидентной и охваченной процессом вытеснения зонах где kj представляет собой коэффициент массопереноса для компонента j [единицы:время-1 ], иxrj и xij(Sgyj + Slxj)(1 - Sw) представляют собой объемные доли компонента j в резидентной зоне и зоне, охваченной процессом вытеснения, соответственно. В уравнении (10) разница объемных долей является движущей силой массопереноса, и коэффициент массопереноса характеризует сопротивление массопереносу. С данным допущением уравнения (6)-(8) принимают вид где DjkjL/u, известный как число Damkhler, представляет собой безразмерный коэффициент массопереноса. Величина числа Damkhler представляет скорость смешения компонентов между зоной, охваченной процессом вытеснения, и резидентной зоной относительно времени пребывания жидкости в керне. ЧислоDamkhler, равное нулю для всех компонентов,подразумевает отсутствие смешения, а высокие числа Damkhler подразумевают быстрое смешение. Данная модель согласуется с предположением, что смешение вызывает перенос компонента из зон более высоких концентраций в зоны более низких концентраций, таким образом,имея тенденцию выравнивать концентрации между двумя зонами. Коэффициенты массопереноса могут являться функциями локальной степени смешиваемости, геометрии элемента сетки, охваченной доли , отношения подвижности (m), скорости (u), гетерогенности и водонасыщенностиkj = kj (степень смешиваемости, геометрия элемента сетки, , m, u, гетерогенность, Sw) (14) Конкретные функциональные зависимости зависят от способов, посредством которых смешиваются вытесняющая и вытесняемая жидкости. Gardner, J.W., and Ypma, J.G.J., "An Investigation of Phase-Behavior/Macroscopic BypassingEngineering Journal, стр. 508-520, Октябрь 1984,описывает эффекты макроскопического просачивания при смешивании в процессах вытеснения агентом, смешивающимся при многократном контакте. Заявители обнаружили, что представленные Gardner и Ypma данные подразумевают, что коэффициенты массопереноса должны быть обратно пропорциональны времени,требующемуся для исключения подсеточных языков посредством поперечной дисперсии где d представляет поперечную ширину элемента сетки, DTj представляет коэффициент поперечной дисперсии компонента j, F является параметром, объясняющим эффекты зоны, охваченной процессом вытеснения, и гетерогенности и Clj является константой, которая может зависеть от компонента j. В первом приближении, коэффициент поперечной дисперсии включает вклады от моле 23 кулярной диффузии, конвективной дисперсии и капиллярной дисперсии. Коэффициент массопереноса модели включает данные вклады и может быть записан в безразмерной форме как где Doj представляет коэффициент молекулярной диффузии для компонента j, T(d) представляет поперечную дисперсность, уmax представляет максимальное межфазное натяжение газ/нефть для несмешивающегося вытеснения,DMj представляет число Damkhler для вытеснения агентом, смешивающимся при первом контакте, С 2 и С представляют корректирующие константы. Термины в первых скобках представляют безразмерные скорости массопереноса вследствие молекулярной диффузии и конвективной дисперсии, соответственно. Молекулярная диффузия преобладает при низкой скорости и незначительной ширине системы, и конвективная дисперсия преобладает при высокой скорости и большой ширине системы (T(d) является возрастающей функцией от d). Термины во вторых скобках объясняют капиллярную дисперсию (необходимо обратить внимание, что когда С равно нулю, т.е. текучие среды являются смешивающимися, Dj и DMj являются синонимическими). Для начальных тестовых целей принимается допущение, что отношение подвижностей и водонасыщенность не оказывают действия на коэффициенты массопереноса. При вытеснении агентом, смешивающимся при многократном контакте, и почти смешивающимся агентом межфазное натяжение зависит от расположения состава элемента сетки внутри двухфазной зоны фазовой диаграммы,чем ближе состав находится к критической точке, тем ниже оказалось бы межфазное натяжение. В контексте настоящей модели, где межфазное натяжение является мерой степени смешиваемости между растворителем и нефтью,межфазное натяжение в уравнении (16) представляет собой натяжение, которое существовало бы между паром и жидкостью, если все содержание элемента сетки находилось бы при равновесии. Следующее уравнение парахора использовали для расчета межфазного натяжения: где Pj представляет параметр парахора для компонента j, xj и yj являются мольными долями компонента j в фазах вытесняемой жидкости и вытесняемого пара, соответственно, 1 и v являются молярными плотностями жидкости и пара и n является показателем степени в диапазоне от 3,67 до 4. Ключевой чертой механистической модели массопереноса, используемой в данном приме 003418 24 ре, является то, что степень смешиваемости между растворителем и нефтью имеет значительное влияние на скорость смешения между зоной, охваченной процессом вытеснения, и резидентной зоной. В предшествующем уровне техники было предположено, что коэффициенты дисперсии несмешивающихся текучих сред в пористой среде могут быть по порядку величины больше, чем смешивающиеся коэффициенты дисперсии при равных экспериментальных условиях. Поэтому смешение должно быть более быстрым при несмешивающихся условиях, чем при смешивающихся условиях. В используемой в примере модели данное наблюдение было вставлено включением зависимости межфазного натяжения при вычислении коэффициента поперечной дисперсии. Поскольку межфазное натяжение зависит от фазового поведения посредством уравнения парахора, уравнение (17), необходимым параметром в контексте модели является константа межфазного натяжения, С. Модель массопереноса вводит ряд параметров (например, коэффициенты диффузии,дисперсность, межфазовое натяжение) в прогнозирующую модель по настоящему изобретению, которые не имеют аналогов в модели смешения Тодда-Лонгстаффа. В то время, как данные дополнительные параметры увеличивают сложность вычислений, все параметры настоящей модели по изобретению имеют физическое значение, которое может быть либо измерено,либо оценено относительно четким образом. Функция подвижности эффективной среды Теория фильтрации и аппроксимация эффективной среды являются известными способами описания критических явлений, проводимости, диффузии и потока в неупорядоченных гетерогенных системах (см., например, Kirkpatrick, S., "Classical Transport in Disordered Media:J. Chem. Phys., 1983, 78, 6849-6864). В контексте проблем потока в гетерогенных системах аппроксимация эффективной среды представляет транспорт в случайной гетерогенной среде в виде транспорта в эквивалентной (эффективной) гомогенной среде. Заявители наблюдали, что соответствие между аппроксимацией эффективной среды и теоретическими результатами является достаточно точным при нахождении далеко от порога фильтрации. Модель подвижности эффективной среды была создана для оценки подвижностей текучих сред в гетерогенной среде. Это сделано при предположении, что распределение растворителя и нефти внутри зоны элемента сетки можно представить случайной перемешанной сеткой 25 двух текучих сред. Нижеследующие аналитические выражения для подвижностей неводной фазы получают, предполагая, что сетка является изотропной и некоррелированной:"разветвленности" перемешанных сеток текучей среды. Увеличение z ведет к большей сегрегации нефти и растворителя так, что прорыв растворителя ускоряется, а добыча нефти задерживается. Относительные проницаемости оценивают, используя насыщенность текучей среды внутри ее зоны. Модель подвижности эффективной среды предлагает приблизительные аналитические выражения для фазовых подвижностей, которые принимают во внимание необходимые свойства (долю охвата, гетерогенность,отношение подвижности) с физической точки зрения. Представленные ниже результаты показывают, что модель подвижности эффективной среды точно схватывает профили добычи при вытеснении смешивающимся агентом. Функция фазового поведения В примерах по настоящему изобретению используют упрощенную псевдотроичную модель фазового поведения для одномерной моделирующей программы. В данной модели составы смесей растворителя и нефти характеризуются с точки зрения трех псевдокомпонентов: СО 2, легкого компонента нефти и тяжелого компонента нефти. Двухфазную оболочку в данной фазовой модели описывают квадратичным уравнением, постоянные которого определяются составами для точки полного смешения частично смешивающихся жидкостей и двумя границами оболочки на граничных условиях. В то время, как она только приблизительно представляет реальную систему, данная фазовая модель успешно моделирует фазовые поведения, 003418 26 соответствующие различным степеням смешиваемости, таким как смешивающийся при первом контакте (СПК), смешивающийся при многократном контакте (СМК) и почти смешивающийся агент (НС). Параметры, определяющие используемую в примерах 1-3 двухфазную оболочку, суммируются в табл. 1. Параметры в табл. 1 для случая СМК определяют описание псевдотроичной фазы системы СO2-нефть Means при 2000 фунт/кв. дюйм абс. (13790 кПа) и 100F(37,78C). Параметры в табл. 1 для случаев СПК и НС определяют описание псевдотроичной фазы, которая может быть получена при 100F(37,78C) и давлениях выше и ниже чем 2000 фунт/кв. дюйм абс. (13790 кПа), соответственно. В составе пластовой нефти преобладают тяжелые компоненты, соответствующие доле тяжелой нефти 0,8434 и доле легкой нефти 0,1566. Параметр В табл. 1 подстрочные индексы 1, 2, 3 обозначают растворитель, тяжелую нефть и легкую нефть, соответственно. V1G и V1L представляют границы двухфазной оболочки. V1G и V1L представляют объемные доли растворителя в газовой и жидкой фазах, соответственно, для смеси растворитель-тяжелая фракция. V1P и V3P представляют объемные доли растворителя и легкой фракции в точке полного смешения частично смешивающихся жидкостей. Параметры, определяющие двухфазную оболочку, используемую в примере 4 (более детально обсуждаемую ниже), суммируются в табл. 2. Используемые в примере 4 параметры определяют описанием псевдотроичной фазы системы СО 2-нефть Wasson при 2000 фунт/кв. дюйм абс. (13,790 кПа) и 100F (37,78 С). Данные были получены из работы Gardner, J.W.,Orr, F.M., and Patel, P.D., "Effect of Phase Behavior on CO2 Flood Displacement Efficiency,"Journal of Petroleum Technology, Ноябрь 1981,стр. 2067-2081. Состав сырой нефти соответствует объемной доле тяжелой нефти, равной 0,72 и объемной доле легкой нефти, равной 0,28. Результаты моделирования Входные данные, использованные при моделировании четырех примеров, принимают данные относительной проницаемости нефтьсоляной раствор и капиллярное давление, являющиеся типичными для карбонатной горной породы Сан Андерс. Свойства керна представляли собой: длина = 1 фут (0,3048 м), пористость = 0,19% и проницаемость = 160 мд (0,1579 мкм 2). Пример 1. Координационное число z в аппроксимации эффективной среды для теории фильтрации означает "разветвленность" или связанность сетки. В контексте настоящего изобретения z представляет структуру языка в элементе сетки и вводит влияние свойств, таких как коэффициент подвижности нефти/растворителя, гетерогенность коллектора и тип горной породы. В общем виде z может быть аналогичным параметру смешенияв модели смешения ТоддаЛонгстаффа. Фиг. 5 А показывает, что увеличение z приводит к уменьшенной добыче нефти, а фиг. 5 В показывает, что увеличение z приводит к более раннему прорыву растворителя. Обе кривые и добычи нефти и прорыва растворителя являются чувствительными к значению z. В частности, изменение z между двумя и пятью снижает добычу нефти при полученных 1,5 объемах порового пространства от 93 до 52% и уменьшает точку, при которой добываемая жидкость достигает концентрации 50% растворителя, от 0,55 до 0,24 полученных объемов порового пространства. Описание фазового поведения СМК в табл. 1 используют в данном примере и принимают, что числа Damkhler равны D1 = 0, D2 = 0,1 и D3 = 0,1. Моделирование данного примера начинают при водонасыщенности остаточной нефти, равной 0,35, и используют 25 элементов сетки в одномерной модели. Увеличение значения z в модели эффективной среды дает эффект, аналогичный уменьшению в значении параметра смешенияв модели смешения Тодда-Лонгстаффа, оба параметра приводят к увеличенному просачиванию нефти (более низкая добыча) и более раннему прорыву растворителя. Координационному числу z может быть установлено значение более чем два или равное двум при осуществлении на практике способа по настоящему изобретению. Z = 2 представляет последовательный поток нефти и растворителя и характеризует 28 вытеснение типа поршневого без образования языков или протоков. Zпредставляет параллельный поток нефти и растворителя и характеризует вытеснение с интенсивным образованием языков или протоков. Основываясь на данных результатах, можно ожидать, что z является важным параметром сопоставления прорыва растворителя и истории добычи нефти. Пример 2. Числа Damkhler представляют скорость смешения компонентов между зоной, охваченной процессом вытеснения, и резидентной зоной. Показанные на фиг. 6 А-D результаты демонстрируют, что данное изобретение успешно воспроизводит точные ограничивающие характеристики. В данном примере используют описание фазового поведения СМК в табл. 1 и предполагают, что числа Damkhler равныD1=0 для компонента растворителя и D2=D3 для компонентов нефти. Моделирование данного примера начинают при водонасыщенности остаточной нефти, равной 0,35, и используют 25 элементов сетки в одномерной модели. Фиг. 6 А показывает, что там, где нет смешения (числа Damkhler нефти = 0), модель точно предсказывает, что имеет место чистое вытеснение нефти без обмена компонентов между зонами. На фиг. 6 А кривая 30 представляет долю добытого легкого компонента нефти, а кривая 31 (которая имеет точно такую же форму, как кривая 30) является долей добытого тяжелого компонента нефти. Кривые добычи легкого и тяжелого компонентов 30 и 31 являются идентичными, что показывает, что состав нефти не изменяется. Когда имеет место быстрое смешение(числа Damkhler нефти больше, чем примерно 5), две зоны быстро достигают примерно одинакового состава. Поэтому результаты моделирования, показанные на фиг. 6D, фактически идентичны результатам традиционной модели одной зоны. На фиг. 6D кривая 60 представляет долю добытого легкого компонента нефти, а кривая 61 является долей добытого тяжелого компонента нефти. Показанные на фиг. 6D результаты также демонстрируют, что при увеличении числаDamkhler в процессе добычи СМК происходит увеличение фракционирования легкого компонента нефти в газовую фазу. Следовательно,легкий компонент предпочтительно извлекают в виде охваченного вытесняющим (высокоподвижным) растворителем и остается остаточная нефть, обогащенная тяжелым компонентом. Фиг. 6B и 6 С показывают результаты для промежуточных скоростей смешения. На фиг. 6 В кривая 40 представляет долю извлекаемого легкого компонента нефти, а кривая 41 является долей извлекаемого тяжелого компонента нефти. На фиг. 6 С кривая 50 представляет долю извлекаемого легкого компонента нефти, а кри 29 вая 51 является долей извлекаемого тяжелого компонента нефти. Данные фигуры показывают,что количество и состав добытой нефти сильно зависят от чисел Damkhler. Таким образом,выбор времени извлечения каждого компонента может быть согласован корректировкой чиселDamkhler. Небольшие изменения в добыче нефти и согласование составов добытой нефти и газа можно осуществить, изменяя числа Damkhler. Пример 3. Фиг. 7 показывает экспериментальные данные, приведенные в статье Blackwell, R.J.,Rayne, J.R., and Terry, W.M., "Factors Influencingthe Efficiency of Miscible Displacement," Petroleum Transactions, AIME (1959) 216, 1-8 (на которую далее идет ссылка как "Blackwell et al.") для вытеснения агентом, смешивающимся при первом контакте, при различных значениях начального отношения вязкостей нефти/растворителя. Экспериментальные данные, которые показаны в виде точек на фиг. 7, получают с использованием гомогенной песчаной упаковки и текучих сред равной плотности (чтобы минимизировать сегрегацию из-за силы тяжести). Эксперименты проводят при отношении вязкостей, равном 5, 86, 150 и 375. Вода не присутствовала в данных экспериментах. На фиг. 7 также показаны линии, которые соответствуют добыче нефти, получаемой при моделировании с использованием способа по настоящему изобретению, в котором начальное отношение вязкостей нефти/растворителя устанавливают при экспериментальном значении, и координационное число корректируют, чтобы получить наиболее возможное соответствие с экспериментальными данными. Оценивают, что число Damkhler равно порядка 10-4 (основываясь на DT = 0,0045 фут 2/день (4,2 см 2/день),= 0,4, L = 6 футов (1,83 м), d = 2 фута (0,61 м) и u= 40 фут/день (12,2 м/день и поэтому принимают, что оно фактически равно нулю. Таким образом, имеется только один корректируемый параметр, использованный при моделировании координационное число z. Используют двадцать пять элементов сетки в одномерной модели. Фиг. 7 показывает точное соответствие между экспериментальными данными Blackwellet al и результатами, полученными по способу настоящего изобретения. В частности, способ настоящего изобретения успешно прогнозирует стабилизацию добычи нефти после начального прорыва. Более того, соответствие между точками данных для вытеснения с неблагоприятным соотношением вязкости оказывается исключительно точным. Поскольку система, использованная Blackwell et al., была смешивающейся при первом контакте и дисперсия была незначительной, ни фазовое состояние, ни массоперенос не играет роли в изменении при моделируемой добыче. Соответствие с экспериментом в данном примере, следовательно, явля 003418 30 ется только подтверждением модели эффективной среды по настоящему изобретению. В то время, как принятая выше процедура может быть приравнена с промысловыми данными с сопоставимой историей для способа по настоящему изобретению, чтобы иметь прогнозирующую способность, было бы необходимо иметь возможность прогнозировать значение z изначально. На выбор z влияет отношение подвижностей, гетерогенность коллектора и тип породы. Фиг. 8 показывает график значений z,которые были использованы для получения соответствия с экспериментальными данными на фиг. 7 в виде функции отношения вязкостей нефти/растворителя. Как показано на фиг. 8, z показывает монотонное изменение с отношением вязкости. Представленные в примерах 1 и 3 результаты показывают, что координационное число,z, является ключевым параметром при реализации на практике настоящего изобретения, поскольку его можно использовать при сопоставлении прорыва растворителя и характеристики добычи нефти. Пример 2 показывает, что тонкую настройку добычи нефти, а также согласование составов добываемых нефти и газа, можно выполнить с помощью модели массопереноса. Используя координационное число z и числа Damkhler в качестве корректирующих параметров и соответствующую фазовую модель для изучаемой системы, можно использовать прогнозирующую модель по настоящему изобретению для сопоставления существенных характеристик (включая добычу нефти, прорыв закачиваемой текучей среды и составы добываемой текучей среды) любого процесса закачивания газа. Пример 3 показывает, что модель подвижности эффективной среды, используемая в способе по настоящему изобретению, можно использовать для описания образования языков и протоков, которые широко распространены в процессах смешивающегося вытеснения. Пример 4. Пример 4 представлен для демонстрации выгоды моделей фазового поведения и массопереноса. Экспериментальные данные, представленные в работах Gardner, J.W., Orr, F.M., and(на которую в дальнейшем идет ссылка "GardnerPetroleum Engineers Journal, стр. 508-520, Октябрь 1984, описывают взаимосвязь между фазовым состоянием и эффективностью вытеснения (добычи нефти) для процессов закачивания смешивающегося газа. Данные статьи представляют результаты экспериментов с кернами наCO2 в системе, смешивающейся при первом контакте (СПК), и вытеснение нефти WassonCO2 в системе, смешивающейся при многократном контакте (МСМ). Soltrol представляет собой продукт, производимый Phillips Petroleum Company, а нефть Wasson получена из месторождения Wasson в западном Техасе. Отношение вязкости нефти/растворителя было 16 для системыCO2/Soltrol и 21 для системы СО 2/нефть Wasson- достаточно близкое, чтобы сделать фазовое поведение единственным основным отличием между двумя системами. Поэтому для всех практических целей единственная причина для различия в добыче для двух систем могла быть приписана изменению в фазовом поведении и макроскопическому просачиванию (как результату изменившегося фазового поведения). Фиг. 9 показывает экспериментальные кривые добычи, полученные для систем(кривая 71). Различные ряды символов обозначают данные, полученные в дублирующих экспериментах с кернами при аналогичных условиях. Все испытания проводят на одном и том же керне Веrеа. Окончательная эффективность добычи нефти оказывается ниже для системы СО 2/нефть Wasson, как и скорость добычи. Образование вязких языков почти полностью отвечает за форму кривой 70 СПК добычиCO2/Soltrol в то время, как образование языков в результате разности вязкостей, так и фазовое поведение отвечают за форму кривой 71 СМК добычи СO2/нефть Wasson. Чтобы протестировать влияние образования языков на добычу,сначала проводят одномерное моделирование,использующее обычную модель одной зоны. Для моделирования данного примера устанавливают такие параметры моделирования, чтобы близко соответствовать экспериментальным системам CO2/Soltrol и СО 2/нефть Wasson. Вязкость CO2 устанавливают при 0,063 сП(0,000063 Па/с) в соответствии с данными, предоставленными Gardner et al. Soltrol имеет нормальный диапазон кипения, эквивалентный диапазону C11-C14, который соответствует вязкости примерно 1,2 сП (0,0012 Па/с). Однако для точного соответствия экспериментальному отношению вязкостей нефти/растворителя, равному 16, принимают, что вязкость Soltrol равна 1,01 сП (0,00101 Па/с). Фазовые вязкости вычисляют по правилу смешения в степени одна четвертая, которое хорошо известно специалистам в данной области. При установлении взаимосвязи относительной проницаемости-насыщенности при моделировании используют экспериментальные отношения относительной проницаемости газ/нефть. Моделирование проводят с 30 элементами сетки. Количество элементов сетки выбирают так, чтобы аппроксимировать уровень продольной дисперсии в эксперименталь 003418 32 ных системах. В случае моделирования СО 2/нефть Wasson выбирают фазовую модель,соответствующую экспериментальной, показанную в табл. 2. Фиг. 9 демонстрирует кривые добычи 72 и 73, полученные моделированием модели одной зоны, наряду с экспериментальными данными (кривые 70 и 71). Кривая 72 иллюстрирует результаты моделирования системыCO2/Soltrol, а кривая 73 иллюстрирует результаты моделирования системы СО 2/нефть Wasson. Из фиг. 9 ясно, что образование языков в результате разницы вязкостей подавляет скорость добычи нефти. Также очевидно, что модель одной зоны предлагает неадекватное описание(как качественно, так и количественно) добычи нефти в системах CO2/Soltrol и СO2/нефть Wasson. Однако модели одной зоны находятся в хорошем согласии с экспериментами с трубами малого диаметра (Gardner et al.), в которых подавлены эффекты просачивания. Чтобы оценить способность способа по настоящему изобретению моделировать экспериментальные данные, полученные на керне, способ по настоящему изобретению сначала применяют к системе СПК CO2/Soltrol. Параметрыz, Dsolvent, DMheavy и DMlight корректируют так,чтобы получить наилучшее соответствие с экспериментальными данными. Принимают, что для простоты DMheavy равен DMlighn. Получают лучшее соответствие для выбора z=4,5,Dsolvent=0, DMheavy, light = 0,5. Используя такие же параметры и принимая С=10, проводят моделирование, используя способ по настоящему изобретению для системы СO2/нефть Wasson. Все параметры моделирования (фазовое состояние, взаимосвязь относительная проницаемость/насыщенность и уровень дисперсии) устанавливают так, чтобы соответствовать экспериментально определенным значениям (данные получены из Gardner et al.). Вязкость нефти при моделировании изменяют для имитации отношения вязкостей нефть Wasson и нефть/растворитель, равного 21. Данные результаты нанесены на фиг. 10. На фиг. 10 кривые 70 и 71, показанные на фиг. 9, снова показаны для сравнения с результатами моделирования, кривой 74 системыCO2/Soltrol, с использованием двухзонной модели по настоящему изобретению, кривой 75 системы CO2/нефть Wasson с использованием двухзонной модели, используемой в способе по настоящему изобретению. Способ по настоящему изобретению дает превосходную работу по сопоставлению СМКCO2/Wasson с использованием тех же параметров, которые применялись для СПК системыCO2/Soltrol. Основная причина для удерживанияz фиксированным из моделирования CO2/Soltrol заключается в том, что поскольку эксперименты с Soltrol и нефтью Wasson выполняют на тех же кернах (такая же степень гетерогенности и тип 33 породы), и при фактически том же отношении вязкости нефти/растворителя (такое же отношение подвижности), значение z должно оставаться, по существу, неизменным. Коэффициенты массопереноса увеличиваются от значений, используемых для наилучшего соответствия для системы CO2/Soltrol. Физически это претворяется в увеличение скоростей массопереноса с уменьшением смешиваемости (от СПК до СМК)- когда смешиваемость уменьшается, капиллярная дисперсия увеличивается, приводя к более высоким скоростям массопереноса. В моделях, представленных в предшествующих примерах, принимают, что резидентная зона остается однофазной жидкостью. Однако,состав резидентной области может входить в многофазную оболочку, если компонентам растворителя обеспечен перенос в данную зону,что может быть исполнено специалистами в данной области. Это потребовало бы дополнительных быстрых вычислений для резидентной зоны и необходимость определения проницаемости как паровой, так и жидкой фазы для данной области. Модель разделенных узлов, используемая в способе по настоящему изобретению, особенно привлекательна для использования при моделировании коллекторов, в которые закачивают растворитель, поскольку все параметры, используемые в модели, имеют физическое значение, которое может быть либо измерено, либо оценено специалистами в данной области. Координационное число z в модели эффективной среды можно скорректировать для соответствия расчету времени добычи закаченной текучей среды. Наблюдается, что z увеличивается с увеличением начального соотношения подвижности нефть/растворитель. Константы Clj в функции массопереноса можно скорректировать для согласования с характеристикой добычи индивидуальных компонентов. Коэффициенты молекулярной диффузии Doj, можно оценить стандартными корреляциями, известными специалистам в данной области. Дисперсностьи константа диффузии С 2 будут зависеть от свойств породы и будут определять масштабирование от лабораторного до масштаба месторождения. В большинстве приложений параметр межфазного натяжения С должен быть постоянным для хорошего приближения. Влияние силы тяжести на относительные подвижности, к которому не обращались в предшествующих примерах, также может быть принято во внимание специалистами в данной области. Например, можно ожидать, что внутри элемента сетки фаза с низкой плотностью будет стремиться к сегрегации вверху элемента сетки и будет иметь более высокую эффективную подвижность в направлении вверх. Также в примерах моделирования не принималась во 34 внимание анизотропия в проницаемости. В моделировании 3-D отсутствие такой анизотропии может иметь тенденцию к переоценке потока в вертикальном направлении. Анизотропная формулировка модели эффективной среды может быть включена в модель специалистами в данной области, но это значительно усложнило бы вычисления. Еще одним фактором, который не рассматривался в настоящих примерах, является присутствие воды в элементах сетки. При моделировании нагнетания вода-сменяющий газ(ВСГ), газ нагнетается только в зону, охваченную процессом вытеснения, а вода нагнетается только в резидентную зону. Таким образом, образование зоны, охваченной процессом вытеснения, можно инициировать только нагнетанием высокоподвижного газа, а не нагнетанием воды. Водонасыщенность также могла бы воздействовать на коэффициенты массопереноса нефть/газ - которые типично включаются в модель. Функция переноса для воды может быть разработана специалистами в данной области так, что вода также может быть разделена между охваченной процессом вытеснения и резидентной зонами. Описаны принцип изобретения и лучший режим, рассматриваемый для применения данного принципа. Для специалистов в данной области будет очевидно, что могут быть сделаны различные изменения в вариантах осуществления на практике, описанных выше, без отклонения от идеи и рамок данного изобретения, как оно сформулировано в следующей ниже формуле изобретения. Следовательно, необходимо понимать, что данное изобретение не ограничивается конкретными деталями, показанными и описанными. ФОРМУЛА ИЗОБРЕТЕНИЯ 1. Способ моделирования одной или нескольких характеристик многокомпонентного,углеводородсодержащего пласта, при котором текучую среду, содержащую, по меньшей мере,один компонент, закачивают в пласт, по меньшей мере, через одну скважину, для вытеснения содержащихся в пласте углеводородов, содержащий следующие стадии:(а) представление пласта, по меньшей мере, в одном направлении в виде множества элементов сетки;(б) деление, по меньшей мере, некоторых элементов сетки на две или более зоны, причем первая зона представляет часть каждого элемента сетки, охваченную вытесняющей текучей средой, а вторая зона представляет часть каждого элемента сетки, по существу, не охваченную закачиваемой текучей средой, причем распределение компонентов в каждой зоне является, по существу, однородным;(в) создание модели, представляющей свойства текучей среды внутри каждой зоны,поток текучей среды между элементами сетки,используя принципы теории фильтрации и скорость переноса компонентов между зонами;(г) использование модели для моделирования одной или нескольких характеристик пласта. 2. Способ по п.1, в котором стадия (г) прогнозирует свойства пласта и текучих сред, которые он содержит, как функцию времени. 3. Способ по п.1, в котором вытесняющую текучую среду смешивают с углеводородами в пласте. 4. Способ по п.1, в котором вытесняющую текучую среду смешивают с углеводородами,присутствующими в пласте, при многократном контакте. 5. Способ по п.1, в котором вытесняющая текучая среда представляет собой диоксид углерода. 6. Способ по п.1, в котором вытесняющая текучая среда содержит углеводородный газ. 7. Способ по п.1, в котором модель, созданная на стадии (в), дополнительно представляет перенос энергии между зонами элемента сетки. 8. Способ по п.1, в котором вытесняющая текучая среда представляет собой пар, и модель из стадии (в) дополнительно представляет перенос энергии между зонами элемента сетки. 9. Способ по п.1, в котором элементы сетки включают неструктурированные элементы сетки. 10. Способ по п.1, в котором элементы сетки являются трехмерными. 11. Способ по п.1, в котором элементы сетки являются двухмерными. 12. Способ по п.1, в котором модель дополнительно учитывает компонент диффузии,дисперсность и межфазное натяжение внутри каждой зоны. 13. Способ по п.1, в котором скорость переноса компонента между зонами пропорциональна временам сопротивления движущей силе. 14. Способ моделирования одной или нескольких характеристик многокомпонентного,углеводородсодержащего пласта, при котором вытесняющую текучую среду закачивают для вытеснения присутствующих в пласте углеводородов, содержащий следующие стадии:(а) представление, по меньшей мере, части пласта в виде множества элементов сетки;(б) деление каждого элемента сетки на две зоны, причем первая зона представляет охваченную растворителем часть каждого элемента сетки, а вторая зона представляет часть каждого элемента сетки, по существу, не охваченную растворителем, причем состав текучей среды в каждой зоне является, по существу, однородным;(в) создание модели, включающей функции, представляющие подвижность каждой фазы в каждой зоне, используя принципы теории фильтрации, функции, представляющие фазовое поведение внутри каждой зоны и функции,представляющие массоперенос каждого компонента между зонами;(г) использование модели в моделирующей программе для моделирования добычи из пласта и для определения одной или нескольких его характеристик. 15. Способ по п.14, в котором стадии от (а) до (г) повторяют в течение множества интервалов времени и используют результаты для прогнозирования свойств пласта и текучих сред,которые он содержит, как функцию времени. 16. Система для определения одной или нескольких характеристик многокомпонентного углеводородсодержащего пласта, в который закачивают вытесняющую текучую среду,имеющую, по меньшей мере, один компонент,для вытеснения пластовых углеводородов, содержащая модель имеющую множество элементов сетки, представляющих пласт, каждый из которых разделен на две зоны, причем первая зона представляет часть каждого элемента сетки, охваченную вытесняющей текучей средой, а вторая зона представляет часть каждого элемента сетки, по существу, не охваченную вытесняющей текучей средой, при этом распределение компонентов в каждой зоне является, по существу, однородным, и подвижность текучих сред в каждой зоне определена на основании принципов теории фильтрации, и моделирующую программу, предназначенную для использования модели для моделирования пласта с целью определения его характеристик. 17. Система по п.16, в которой модель представляет свойства текучей среды внутри каждой зоны, поток текучей среды между элементами сетки и перенос компонентов между зонами. 18. Способ моделирования изменения, по меньшей мере, одного компонента многокомпонентной системы текучей среды в углеводородсодержащем пласте,характерные черты которого описываются рядом уравнений, посредством моделирующей программы на компьютере, содержащий следующие стадии:(а) обеспечение модели, имеющей каждый элемент сетки разделенный на две зоны, причем первая зона представляет часть каждого элемента сетки, охваченную вытесняющей текучей средой, а вторая зона представляет часть каждого элемента сетки, по существу, не охваченную вытесняющей текучей средой, причем распределение компонентов в каждой зоне является,по существу, однородным и подвижность текучих сред в каждой зоне определяют на основании принципов теории фильтрации;