Кислотное воздействие на горные породы (стр. 35 )

Однако рассмотрение сложных неоднородных сред, в том числе трещиноватых, с применением аналитических базисных функций как, например, в указанной выше работе [85], представляется весьма затруднительным ввиду практической неосуществимости их построения для произвольной системы трещин. Поэтому базисные функции определяются численно из решения задачи на исходной подробной сетке, неперекрывающееся объединение ячеек которой образует ячейки грубой сетки.

Ключевым моментом [75] при определении базисных функций является выбор граничных условий, специальный вид которых может существенно повысить качество решения задачи в целом. Действительно [86], стандартные для методов конечных элементов линейные базисные функции не отражают особенности решения в случае, когда разрывы коэффициентов задачи выходят на границу ячейки грубой сетки. В этом случае можно использовать граничные условия, полученные из решения одномерных задач на ребрах [87]. Следует отметить, что данный подход очевидным образом обобщается с трехмерным случаем [51]. Как показано в [79, 88, 89], использование базисных функций, являющихся решениями локальных задач, может привести к неточным результатам при наличии в среде протяженных высокопроводящих каналов, которыми, в частности, являются рассматриваемые здесь трещины. В работе [79] предлагается метод построения граничных условий из решения полной задачи на подробной сетке в начальный момент времени. Таким образом, используется глобальная информация, что улучшает [79, 88, 89] качество получаемого с помощью построенного базиса решения.

Другой особенностью применения базисных функций, содержащих информацию о мелкомасштабных неоднородностях, является восстановление корректного распределения давления на подробной сетке, которое может быть использовано для построения поля скоростей и решения задач переноса флюидов на этой сетке. Такой подход получил название демасштабирование, по аналогии с прямым процессом ремасштабирования. Его применение дает, например, возможность напрямую учесть эффекты временной разномасштабности потоков флюидов в среде с наличием трещин, тем самым [90] значительно повысить качество получаемого результата по сравнению, например, с использованием стандартного метода осреднения.

В [91] отмечается, что движение жидкостей и газов, и в особенности многофазных ве­ществ, в трещиновато-пористом пласте обладает целым рядом особенностей. В трещиновато-пористом пласте емкостью и проводимостью обла­дают как блоки породы, так и сами трещины. Если блоки породы непроницаемые, то систему трещин можно считать своеобразной фильтрующей средой. Уравнение неустано­вившегося движения однородной жидкости в такой среде (в случае, конечно, слабой сжимаемости среды) будет вполне аналогичным уравнению движения однородной жидкости в обычной пористой среде, т. е. уравнению теплопроводности. Если же трещины в трещиновато-пористом пласте каким-то образом сделать непроницаемыми в продольном направлении, но проницаемыми в поперечном направлении, то, учитывая, что объем трещин обычно невелик по сравнению с поровым объемом блоков, трещиновато-пористый пласт превратится практически в обычную пористую среду. Таким образом, трещиновато-пористая среда может в пределе «превращаться» как в среду с чисто трещинной пористостью, так и в обычную пористую среду. В общем же случае эта среда содержит признаки как пористой, так и чисто трещинной среды. Поэтому при математическом описании движения однородной жидкости в трещиновато-пористой среде естественно представить эту среду в виде «вложенных» друг в друга пористой и трещинной сред. При установившемся движении жидкости в трещиновато - пористой среде эта среда будет вести себя как среда, проводимость которой равна сумме проводимостей пористой и трещинной сред. Если же движение жидкости в трещиновато-пористой среде неустановившееся, вступит в действие явление обмена жидкостью между системой блоков и системой трещин. Для математического описания движения жидкости в трещиновато-пористой среде вводят два понятия скорости фильтрации - скорость фильтрации в системе трещин и скорость фильтрации в системе блоков , два давления - давление в системе трещин и давление в блоках , обмен жидкостью между блоками и трещинами происходит квазистационарно, т. е. явно не зависит от времени t, получают следующую формулу для скорости этого обмена в элементарном объеме пласта:

                                       (6.1.1.1.1),

где ? - коэффициент, характеризующий интенсивность обмена жидкостью между системой блоков и системой трещин, зависящий от проницаемости блоков и их геометрической характеристики.

Уравнения неразрывности течения однородной жидкости в трещиновато-пористой среде в соответствии со сказанным выше имеют следующий вид

,        (6.1.1.1.2),

где и - пористости системы трещин и блоков; ? - плотность жидкости.

Принимая зависимости скоростей фильтрации в трещинах и блоках от соответствующих градиентов давлений в форме закона Дарси получают замкнутую систему уравнений, описывающих движение однородной жидкости в трещиновато-пористой среде.

Закон обмена жидкостью между блоками и трещинами может быть представлен, конечно, и в форме, явно учитывающей нестационарность, т. е. зависимость этого процесса от времени. Такой закон обмена жидкостью между блоками и трещинами был введен  [91] и [93]. В работе [94] предлагается формула обмена между блоками и трещинами без использования понятия давления в блоках

                               (6.1.1.1.3),

где ро - начальное давление в трещинах; ?2 - упругоемкость блоков.

Систему уравнений (1.1)-(1.3) можно также считать системой уравнений движения однородной жидкости в среде с двойной пористостью и использовать эту систему для описания соответствующего течения жидкости в пористом пласте с сильно развитой литологической неоднородностью.

При использовании уравнений (1.1)-(1.3) применительно к трещиновато-пористой среде в ряде случаев можно сделать дальнейшие упрощения. Так, например, при большей проводимости трещин по сравнению с проводимостью блоков распределение давле­ния в них можно считать квазистационарным, а также считать, что блоки являются своего рода «источниками», питающими систему трещин. Тогда уравнения (6.1)-(6.3) упростятся и примут следующий вид:

, .  (6.1.1.1.4),

При исключении из уравнений (1.4) давления получаем одно уравнение

.  (6.1.1.1.5),

При быстром изменении давления на границе трещиновато-пористого пласта давление в трещинах в непосредственной близости от границы пласта принимает значение, близкое к давлению на границе пласта. Давление же в блоках р2 вблизи границы пласта может, как это следует непосредственно из (6.5), существенно отличаться от давления в трещинах . Разность давлений в блоках и трещинах вблизи границы пласта - «скачок давления» [95] уменьшается со временем по экспоненциальному закону.

Механизм движения однородных и неоднородных жидкостей, а также газа и газированных жидкостей в трещиновато-пористых пластах в настоящее время изучен в соответствии с приведенными выше представлениями об этих пластах.

Большое значение для практики имеет процесс вытеснения нефти водой из трещиновато-пористых пластов. В тех случаях, когда пористый коллектор нефти является гидрофильным, при контакте воды с этим коллектором происходит капиллярная пропитка. Если пористые блоки хорошо смачиваются водой, то при закачке в трещиновато-пористый пласт воды она вытесняет нефть из трещин, а из блоков породы в трещины поступает нефть за счет капиллярной пропитки. Если закачивать в пласт воду, то она будет поступать в блоки породы за счет капиллярной пропитки. Экспериментальные данные показывают, что скорость капиллярного впитывания воды (расход воды, впитывающейся в породу в единицу времени, равный расходу нефти, выходящей из породы) можно представить следующим образом [96]:

,   (6.1.1.1.6),

где А - коэффициент, определяемый экспериментальным путем; , - проницаемость и пористость блоков породы.

При закачке воды в трещиновато-пористый пласт с расходом q (?) процессом капиллярной пропитки охватываются не одновременно все блоки пласта. Поэтому обозначим символом ?(?) время начала капиллярной пропитки блока породы, находящегося на расстоянии (- размер блока породы) от начала координат. Если капиллярная пропитка блока породы началась в момент времени ?(?), то скорость капиллярной пропитки этого блока в момент времени ? равна ?[? - ?(?)]. Считая, что в каждый момент времени расход воды идет на капиллярную пропитку блоков, получаем