СПОСОБ ПРОГНОЗИРОВАНИЯ ИЗМЕНЕНИЯ СВОЙСТВ ПРИЗАБОЙНОЙ ЗОНЫ ПЛАСТА ПОД ВОЗДЕЙСТВИЕМ БУРОВОГО РАСТВОРА
PREDICTION OF BOTTOMHOLE FORMATION ZONE UNDER ACTION OF DRILLING MUD

Изобретение относится к способам прогнозирования изменения характеристик призабойной зоны нефте/газосодержащего пласта под воздействием проникающего бурового раствора (или промывочных жидкостей) в ходе бурения, заканчивания и очистки скважин.
Буровые растворы представляют собой сложные смеси глины, твердых частиц (размером от менее одного микрона до нескольких миллиметров) и органических добавок (полимеров, крахмала, ПАВ и т.п.), содержащихся в "несущей" жидкости - "основе" бурового раствора, в качестве которой может выступать вода, нефть, дизель или какая-либо синтетическая жидкость.
При бурении и заканчивании скважин на репрессии компоненты бурового раствора проникают в пласт и вызывают изменение его фильтрационно-емкостных свойств, формируется сложная структура призабойной зоны пласта. Обычно выделяется внешняя фильтрационная корка (образующаяся на стенке скважины и состоящая из отфильтрованных твердых частиц и глины бурового раствора), внутренняя фильтрационная корка (формируемая частицами бурового раствора, проникающими в пластовые породы) и зона проникновения фильтрата бурового раствора.
В ходе технологической процедуры очистки скважины (путем постепенного ее вывода на добычу) компоненты бурового раствора частично вымываются из околоскважинной зоны. Тем не менее часть компонентов остается удержанной в поровом пространстве породы коллектора, что вызывает ухудшения характеристик призабойной зоны, прежде всего, снижение проницаемости. Как результат, снижается продуктивность скважины. Для характеризации этого явления обычно используется термин "повреждение призабойной зоны пласта" или кратко "повреждение пласта".
В большей части математических моделей, используемых для описания процесса проникновения бурового раствора, учитывается только образование/удаление внешней фильтрационной корки и проникновение фильтрата бурового раствора (вода/нефть/дизель) в пластовые породы - считается, что частицы бурового раствора в пластовые породы не проникают.
Исторически численное моделирование динамики проникновения фильтрата бурового раствора впервые стало использоваться в методиках оценки свойств пласта путем сравнения рассчитанного и реального радиального распределения электросопротивления по данным многозондового индукционного каротажа, а позднее - по результатам одновременной инверсии такого каротажа в сочетании с динамикой давления, зарегистрированной в ходе отбора глубинных проб из скважины. Возможное изменение проницаемости и иных свойств призабойной зоны пласта под воздействием проникших компонент бурового раствора в таких методиках учитывается лишь как дополнительное гидравлическое сопротивление (связанное с образованием внутренней и остатками внешней фильтрационной корки).
Более детальный учет механизмов изменения свойств призабойной зоны пласта под воздействием проникших компонент бурового раствора содержится, например, в патенте US 7099811. Описанный в данном патенте способ базируется на лабораторных фильтрационных экспериментах с длинными образцами керна (40 см), включающих закачку бурового раствора в керн и последующую закачку пластовой жидкости с противоположного конца керна с целью моделирования процесса очистки. Профили поврежденной и восстановленной проницаемости вдоль образца керна, полученные в ходе лабораторных фильтрационных экспериментов, используются в качестве входных параметров для гидродинамического симулятора, который учитывает распределение проницаемости в призабойной зоне, используя цилиндрическую сетку с очень мелкими ячейками (порядка нескольких миллиметров) вокруг скважины. Динамика изменения проницаемости моделируется как функция суммарного расхода жидкостей через образец керна. Недостатком данного подхода является необходимость проведения отдельного, достаточно длительного, лабораторного эксперимента для анализа даже относительно небольшой вариации концентрации компонент бурового раствора, фильтрационных свойств образца или перепада давления при бурении/очистки скважины. Тем самым для выбора оптимального состава бурового раствора, а также выбора и прогноза режима бурения/очистки скважины требуется выполнить большое число лабораторных экспериментов, особенно если речь идет о бурении скважины в сложном многопластовом месторождении.
Предлагаемый способ обеспечивает возможность прогнозирования изменения свойств призабойной зоны пласта, таких как проницаемость и насыщенность, в ходе бурения/очистки скважины, возникающего из-за проникновения/удаления компонент бурового раствора. Указанный технический результат обеспечивается за счет комбинирования математического моделирования и лабораторных фильтрационных экспериментов и использования рентгеновской компьютерной микротомографии для получения кривых концентрации проникших твердых частиц бурового раствора, что накладывает дополнительное ограничение при решении обратной задачи определения параметров модели.
Предлагаемый способ осуществляется следующим образом.
Определяют реологические свойства бурового раствора, фильтрата бурового раствора и пластового флюида. Измеряют свойства внешней фильтрационной корки, а также пористость и проницаемость образца керна. Создают математическую модель внешней фильтрационной корки для описания роста внешней фильтрационной корки и просачивания через нее фильтрата бурового раствора.
Прокачивают буровой раствор через образец керна и регистрируют динамику перепада давления на образце и расхода истекающей из образца жидкости. Дополнительно может быть зарегистрирована динамика профиля насыщенности фильтрата бурового раствора посредством, например, измерений электросопротивления по длине образца в процессе прокачки бурового раствора.
Осуществляют рентгеновскую компьютерную микротомографию образца керна и на основе анализа полученных данных микротомографии определяют профиль концентрации проникших в образец твердых частиц бурового раствора. Создают математическую модель внутренней фильтрационной корки для описания динамики изменения концентрации частиц бурового раствора в поровом пространстве образца керна и сопутствующего изменения проницаемости образца керна.
Создают сцепленную математическую модель внешней и внутренней фильтрационных корок, включающую модель внешней фильтрационной корки и модель внутренней фильтрационной корки. Используя созданную сцепленную математическую модель внешней и внутренней фильтрационных корок, а также измеренные ранее свойства внешней фильтрационной корки, определяют параметры математической модели внутренней фильтрационной корки для описания динамики изменения концентрации частиц бурового раствора в поровом пространстве образца керна, при которых одновременно воспроизводятся данные эксперимента по прокачке бурового раствора через образец керна (динамика перепада давления на образце, расход истекающей из образца жидкости и, дополнительно, профиль насыщенности фильтрата бурового раствора) и профиль концентрации проникших частиц бурового раствора. Используют полученные параметры математической модели внешней фильтрационной корки и математической модели внутренней фильтрационной корки для прогнозирования изменения свойств призабойной зоны пласта под воздействием проникших компонент бурового раствора.
Определяемыми параметрами математической модели внутренней фильтрационной корки для описания динамики изменения концентрации частиц бурового раствора в поровом пространстве образца керна являются показатель, определяющий интенсивность снижения проницаемости породы с ростом доли удержанных в поровом пространстве частиц бурового раствора, коэффициент захвата частиц i-го типа, коэффициент мобилизации частиц i-го типа, критическая скорость мобилизации частиц i-го типа, эмпирический коэффициент, зависящий от структуры порового пространства и механизма захвата частиц в порах.
Прокачка бурового раствора через образец керна может осуществляться с циркуляцией раствора по касательной к входному торцу образца керна.
Дополнительно для моделирования процесса удаления ("отмыва") компонент бурового раствора через образец керна может быть осуществлена прокачка пластовой жидкости, при этом закачка пластовой жидкости осуществляется с торца, противоположного торцу, с которого осуществлялась закачка бурового раствора ("обратная закачка").
Предпочтительно, перед прокачкой бурового раствора через образец керна осуществляют предварительную рентгеновскою компьютерную микротомографию этого образца.
Если компоненты бурового раствора, влияющие на фильтрационные свойства керна, имеют слабый рентгеновский контраст (например, глина), в раствор добавляют контрастный агент.
Если размеры образца керна не позволяют провести рентгеновскую компьютерную микротомографию с требуемым разрешением, из образца высверливают образец малого размера, например цилиндр, и осуществляют микротомографию этого малого образца.
Используя описанную выше последовательность действий, могут быть получены параметры математических моделей внешней и внутренней фильтрационных корок для различных типов буровых растворов и образцов кернов различных литотипов, что позволяет построить зависимости параметров моделей от свойств образцов кернов и буровых растворов. Полученные зависимости используют для моделирования динамики проникновения или удаления компонентов различных буровых растворов в призабойной зоне и прогнозирования сопутствующего изменения свойств призабойной зоны в предполагаемых условиях бурения/очистки. Свойствами кернов являются проницаемость, характерный размер пор, смачиваемость.
Изобретение поясняется чертежами, где на фиг.1 показаны внешняя и внутренняя фильтрационные корки, на фиг.2 показана качественная схема типичного эксперимента по исследованию свойств внешней фильтрационной корки путем продавливания бурового раствора через фильтровальную бумагу или керамический диск, на фиг.3 приведен пример изменения проницаемости в ходе лабораторного фильтрационного эксперимента с глинистым буровым раствором и образцом керна, на графике обозначены типичные стадии процесса, на фиг.4 показаны примеры высверленного мини-цилиндра малого размера для проведения рентгеновской компьютерной микротомографии, на фиг.5 изображен поперечный разрез керна на рентгеновской микрокомпьютерной томографии, содержащий глинистые частицы бурового раствора, модифицированные добавлением специального контрастного элемента, на фиг.6 показан пример реконструированного профиля концентрации проникших частиц бурового раствора вдоль образца керна, на фиг.7 показаны регионы пространственной сетки в околоскважинной зоне пласта для моделирования проникновения в пласт компонентов бурового раствора, на фиг.8 показана рассчитанная динамика потерь бурового раствора в ходе бурения, а также рассчитанные профили распределения фильтрата бурового раствора и поврежденной проницаемости (в данном случае проведена нормировка на начальную проницамость) в призабойной зоне пласта после 4-х часов бурения.
Изобретение осуществляется следующим образом. В соответствии с заданной рецептурой подготавливают буровой раствор и берут образец керна. Рецептура любого бурового раствора включает в себя концентрацию и тип твердой фазы (глина, частицы), используемой для изготовления бурового раствора. Как показано на фиг.1, компоненты бурового раствора проникают в направлении от контакта со скважиной 3 в пласт 4 и вызывают изменение его фильтрационно-емкостных свойств, формируется сложная структура призабойной зоны пласта. Обычно выделяется внешняя фильтрационная корка 1 (образующаяся на стенке скважины 3 и состоящая из отфильтрованных твердых частиц и глины бурового раствора), внутренняя фильтрационная корка 2 (формируемая частицами бурового раствора, проникающими в пластовые породы) и зона проникновения фильтрата бурового раствора.
Определяют реологические свойства бурового раствора, фильтрата бурового раствора и пластового флюида согласно, например, РД 39-2-645-81 "Методика контроля параметров буровых растворов" и ГОСТ 33-2000 "Нефтепродукты. Прозрачные и непрозрачные жидкости. Определение кинематической вязкости и расчет динамической вязкости".
Проводят эксперимент по определению показателя фильтрации бурового раствора согласно РД 39-2-645-81 "Методика контроля параметров буровых растворов" или с помощью иных методик с фильтр-прессом по продавливанию бурового раствора через фильтровальную бумагу или керамический диск (см., например, Dewan J.T., Chenevert М.Е. A model for filtration of Water-base Mud During Drilling: Determination of Mudcake Parameters, Petrophysics, 2001, V. 42, №3, pp.237-250). В ходе данного эксперимента твердая фаза бурого раствора отфильтровывается на поверхности фильтровальной бумаги (или керамического диска в зависимости от типа эксперимента), образуя малопроницаемую фильтрационную корку - аналог внешней фильтрационной корки, формирующейся на стенке скважины в ходе бурения.
mc
На фиг.2 приведена качественная схема типичного эксперимента по исследованию свойств внешней фильтрационной корки путем продавливания бурового раствора 5 с помощью пресса 6 через фильтровальную бумагу (или керамический диск) 7, в ходе чего формируется слой отфильтровавшейся твердой фазы (внешняя фильтрационная корка) 8 толщины h, а фильтрат бурового раствора 9 собирается в специальном мерном стакане.
Создают математическую модель внешней фильтрационной корки для описания роста внешней фильтрационной корки и просачивания через нее фильтрата бурового раствора. Для математического моделирования используется система уравнений, основанных на балансе массы частиц и фильтрата:
<mrow><mfrac><mrow><mi>d</mi><msub><mi>h</mi><mrow><mi>m</mi><mi>c</mi></mrow></msub></mrow><mrow><mi>d</mi><mi>t</mi></mrow></mfrac><mo>=</mo><mo>−</mo><mfrac><mrow><msub><mi>C</mi><mi>Σ</mi></msub><msub><mi>u</mi><mrow><mi>m</mi><mi>f</mi></mrow></msub></mrow><mrow><mrow><mo>(</mo><mrow><mn>1</mn><mo>−</mo><msub><mi>C</mi><mi>Σ</mi></msub></mrow><mo>)</mo></mrow><mrow><mo>[</mo><mrow><mn>1</mn><mo>−</mo><msub><mi>ϕ</mi><mrow><mi>m</mi><mi>c</mi></mrow></msub><mrow><mo>(</mo><mrow><mi>Δ</mi><msub><mi>p</mi><mrow><mi>m</mi><mi>c</mi></mrow></msub></mrow><mo>)</mo></mrow></mrow><mo>]</mo></mrow></mrow></mfrac><mo>−</mo><msub><mi>f</mi><mi>e</mi></msub><mo>,</mo><mtext>  </mtext><mtable><mtr><mtd><mrow /></mtd><mtd><mrow /></mtd><mtd><mrow /></mtd><mtd><mrow /></mtd></mtr></mtable><mtext> (1)</mtext></mrow>
<mrow><mi>d</mi><mi>i</mi><mi>v</mi><mtext> </mtext><msub><mi>u</mi><mrow><mi>m</mi><mi>f</mi></mrow></msub><mo>=</mo><mn>0,</mn><mtext> </mtext><mtable><mtr><mtd><mrow /></mtd><mtd><mrow /></mtd><mtd><mrow /></mtd><mtd><mrow /></mtd></mtr></mtable><mtext>  (2)</mtext></mrow>
а также законе Дарси для скорости фильтрации через фильтрационную корку:
<mrow><msub><mi>u</mi><mrow><mi>m</mi><mi>f</mi></mrow></msub><mo>=</mo><mo>−</mo><mfrac><mrow><msub><mi>k</mi><mrow><mi>m</mi><mi>c</mi></mrow></msub><mrow><mo>(</mo><mrow><mi>Δ</mi><msub><mi>p</mi><mrow><mi>m</mi><mi>c</mi></mrow></msub></mrow><mo>)</mo></mrow></mrow><mrow><msub><mi>μ</mi><mrow><mi>m</mi><mi>f</mi></mrow></msub></mrow></mfrac><mfrac><mrow><mi>Δ</mi><msub><mi>p</mi><mrow><mi>m</mi><mi>c</mi></mrow></msub></mrow><mrow><msub><mi>h</mi><mrow><mi>m</mi><mi>c</mi></mrow></msub></mrow></mfrac><mo>,</mo><mtext> </mtext><mtable><mtr><mtd><mrow /></mtd><mtd><mrow /></mtd><mtd><mrow /></mtd><mtd><mrow /></mtd></mtr></mtable><mtext>  (3)</mtext></mrow>
mc
e
mf
mf
mc
mc
mc
Σ
здесь фильтрат бурового раствора предполагается несжимаемым; t - время; параметры с нижним индексом 'mc' соответствуют внешней фильтрационной корке; h - толщина внешней фильтрационной корки; f - скорость эрозии ("смыва") внешней фильтрационной корки из-за циркуляции бурового раствора в затрубном пространстве; u - скорость фильтрации фильтрата бурового раствора через фильтрационную корку; µ - вязкость фильтрата бурового раствора; Δp - перепад давления на внешней фильтрационной корке; k и φ - проницаемость и пористость внешней фильтрационной корки соответсвенно; C - объемная концентрация твердой фазы в буровом растворе (твердая фаза в общем случае может состоять из нескольких различных компонент).
Уравнение динамики роста внешней фильтрационной корки (1) включает член, описывающий ее эрозию ("смыв") под воздействием циркуляции бурового раствора в затрубном пространстве.
mc
mc
mc
mc
Система уравнений (1) - (3) включает два замыкающих соотношения k(Δp) и φ(Δp), описывающие изменение проницаемости и пористости внешней фильтрационной корки при изменении перепада давления на ней (т.е. сжимаемость внешней фильтрационной корки). Следующая степенная зависимость подтверждена экспериментально (см., например, Dewan J.T., Chenevert М.Е. A model for filtration of Water-base Mud During Drilling: Determination of Mudcake Parameters, Petrophysics, 2001, V.42, №3, pp.237-250):
<mrow><msub><mi>k</mi><mrow><mi>m</mi><mi>c</mi></mrow></msub><mo>=</mo><mfrac><mrow><msub><mi>k</mi><mrow><mi>m</mi><mi>c</mi><mn>0</mn></mrow></msub></mrow><mrow><msup><mrow><mrow><mo>(</mo><mrow><mi>Δ</mi><msub><mi>p</mi><mrow><mi>m</mi><mi>c</mi></mrow></msub></mrow><mo>)</mo></mrow></mrow><mi>ν</mi></msup></mrow></mfrac><mo>,</mo><mtext>  </mtext><msub><mi>ϕ</mi><mrow><mtext>mc</mtext></mrow></msub><mo>=</mo><mfrac><mrow><msub><mi>ϕ</mi><mrow><mi>m</mi><mi>c</mi><mn>0</mn></mrow></msub></mrow><mrow><msup><mrow><mrow><mo>(</mo><mrow><mi>Δ</mi><msub><mi>p</mi><mrow><mi>m</mi><mi>c</mi></mrow></msub></mrow><mo>)</mo></mrow></mrow><mi>δ</mi></msup></mrow></mfrac><mo>,</mo><mtext>   </mtext><mtable><mtr><mtd><mrow /></mtd><mtd><mrow /></mtd><mtd><mrow /></mtd><mtd><mrow /></mtd></mtr></mtable><mtext>(4)</mtext></mrow>
mc0
mc0
mc
где k, φ - опорные значения проницаемости и пористости, определяемые экспериментально при Δp=1 psi; ν, δ - показатели сжимаемости, также получаемые экспериментальным путем.
mc0
mc0
Следующие свойства внешней фильтрационной корки неизвестны: проницаемость k, пористость φ, показатели сжимаемости ν и δ.
Эти свойства определяют (см., например, Dewan J.T., Chenevert М.Е. А model for filtration of Water-base Mud During Drilling: Determination of Mudcake Parameters, Petrophysics, 2001, V. 42, №3, pp.237-250) исходя из данных проведенного ранее эксперимента по продавливанию бурового раствора через фильтровальную бумагу или керамический диск.
Далее измеряют открытую пористость и абсолютную проницаемость исходного образца керна, например, согласно ГОСТ 26450.1-85 и ГОСТ 26450.2-85.
Проводят лабораторный фильтрационный эксперимент по прокачиванию бурового раствора через образец керна. В состав лабораторного оборудования для фильтрационных экспериментов с прокачиванием бурового раствора входит ячейка фильтрации (кернодержатель), а также может входить блок измерения электросопротивления в различных частях образца керна вдоль оси фильтрации. Закачка бурового раствора в образец керна может быть реализована в двух режимах: с циркуляцией раствора по касательной к входному торцу образца керна и без таковой. В целях воспроизведения динамического процесса проникновения бурового раствора проводят лабораторные эксперименты в репрезентативных скважинных условиях (температура, репрессия на пласт и скорость сдвига применительно к буровому раствору). После прокачивания бурового раствора может быть осуществлена прокачка пластовой жидкости с противоположного конца керна ("обратная закачка") с целью моделирования процесса удаления ("отмыва") компонент бурового раствора.
В зависимости от конкретной конструкции кернодержателя, в нем могут быть исследованы образцы керна диаметром от 0.8 до 5 см и длиной до 30-40 см. Например, стандартный кернодержатель предназначен для исследования образцов керна диаметром 3 см и длиной 3-6 см, однако существуют также различные модификации более длинного кернодержателя с несколькими портами (обычно, 4-5) для измерения перепада давления в различных частях керна в ходе фильтрационного эксперимента (см. например, Longeron D., Alfenore J., Salehi N., Saintpére S., 2000, Experimental approach to characterize drilling mud invasion, formation damage and cleanup efficiency in horizontal wells with openhole completions. SPE 58737).
В ходе фильтрационного эксперимента с буровым раствором регистрируют динамику перепада давления Δp(t) на образце керна и объемного расхода истекающей из керна жидкости q(t). Могут быть произведены измерения электросопротивления по длине образца, что позволяет зарегистрировать динамику профиля насыщенности фильтрата бурового раствора. На практике, как правило, используется динамика снижения / восстановления проницаемости образца керна в зависимости от объема закачанной жидкости, которая рассчитывается исходя из зарегистрированной динамики перепада давления и расхода на основе закона Дарси:
<mrow><mi>q</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mo>−</mo><mfrac><mrow><msub><mi>k</mi><mrow><mi>c</mi><mi>o</mi><mi>r</mi><mi>e</mi></mrow></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow><mrow><msub><mi>μ</mi><mo>*</mo></msub></mrow></mfrac><mi>A</mi><mfrac><mrow><mi>Δ</mi><mi>p</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow><mi>L</mi></mfrac><mtext>   (5)</mtext></mrow>
core
*
где t - время; k(t) - проницаемость керна; Δp(t) - перепад давления на образце керна; q(t) - объемный расход истекающей из керна жидкости; L - длина керна; А - площадь поперечного сечения керна; µ - вязкость истекающай из керна жидкости (на практике при прокачивании бурового раствора эта величина предполагается равной вязкости фильтрата бурового раствора, а при обратной прокачке - вязкости закачиваемой пластовой жидкости).
Пример изменения проницаемости в ходе лабораторного фильтрационного эксперимента по прокачиванию глинистого бурового раствора через образец керна и последующей обратной прокачке пластовой жидкости с противоположного конца керна приведен на фиг.3. Буровой раствор прокачивали при фиксированном перепаде давления на режиме циркуляции раствора по касательной к входному торцу образца керна.
В ходе прокачивания бурового раствора (см. фиг.3) выделяются три стадии изменения проницаемости: быстрое падение на начальном этапе 10, обусловленное формированием внутренней фильтрационной корки, последующее более плавное снижение проницаемости 11, связанное с ростом внешней фильтрационной корки, и конечная стадия 12, когда проницаемость керна остается практически постоянной в силу стабилизации толщины внешней фильтрационной корки из-за ее эрозии под воздействием циркулирующего бурового раствора и слабого изменения параметров внутренней фильтрационной корки. В процессе обратной прокачки происходит полный или частичный срыв внешней фильтрационной корки (при превышении перепадом давления некоторой пороговой величины) и последующий отмыв части внутренней фильтрационной корки 13, благодаря чему проницаемость керна частично восстанавливается.
После лабораторного эксперимента с буровым раствором осуществляют рентгеновскую компьютерную микротомографию образца керна. Может быть проведена предварительная рентгеновская компьютерная микротомография исходного образца керна, до воздействия на него бурового раствора - это позволяет улучшить процедуру анализа данных рентгеновской микрокомпьютерной томографии. Если компоненты бурового раствора, влияющие на фильтрационные свойства керна, имеют слабый рентгеновский контраст (например, глина), в раствор добавляют контрастный агент (например, согласно методике, описанной в патенте РФ №2467316).
Если размеры образца керна не позволяют провести рентгеновскую компьютерную микротомографию с требуемым разрешением, из образца высверливают образец малого размера, например цилиндр, и осуществляют микротомографию этого малого образца (см. фиг.4).
Пример поперечного разреза образца керна на рентгеновской микрокомпьютерной томографии, содержащего глинистые частицы бурового раствора, модифицированные добавлением специального контрастного элемента, показан на фиг.5.
Специальный анализ данных рентгеновской микрокомпьютерной томографии (см. например, US №20130010919) позволяет получить профиль концентрации проникших в образец твердых частиц бурового раствора, пример которого приведен на фиг.6.
Создают математическую модель внутренней фильтрационной корки для описания динамики изменения концентрации частиц бурового раствора в поровом пространстве образца керна (иными словами, модель внутренней фильтрационной корки) и сопутствующего изменения проницаемости образца керна.
Математическая модель переноса и изменения концентрации (накопления/отмыва) твердых частиц бурового раствора в пористой среде, в отличие от традиционной модели (см. например Herzig J.P., Leclerc D.M., Le Goff P., 1970, Flow of Suspensions through Porous Media - Application to Deep Filtration, Industrial and Engineering Chemistry. Vol.62. № 5), предполагающей частицы одинаковыми и описывающей процесс их захвата/мобилизации в пористой среде с помощью единственного кинетического уравнения, обобщается для описания переноса двух разных видов твердых частиц. Различие физических свойств твердых частиц (например, размера) отражено путем введения двух кинетических уравнений с существенно разными коэффициентами захвата (и мобилизации). При необходимости, добавляется третий и последующие типы частиц.
Если проницаемость пласта на границе со скважиной снижается ниже определенной критической величины, считается, что все частицы, содержащиеся в буровом растворе, осаждаются (отфильтровываются) на стенке скважины и формируют растущую со временем внешнюю фильтрационную корку.
В настоящем изобретении мы рассматриваем проблему проникновения бурового раствора при следующих допущениях: два флюида (фильтрат бурового раствора и пластовый флюид (нефть, газ или вода)); линейный закон Дарси; буровой раствор включает два типа несжимаемых частиц; скорости частиц совпадают со скоростью несущей жидкости (отсутствие проскальзывания); диффузией частиц пренебрегается; частицы могут быть либо мобильными, либо захваченными поровыми ловушками.
Прежде чем сформулировать систему уравнений, введем основные определения.
bulk
Рассмотрим некоторый представительный объем V пористой среды (образец породы).
<mrow><msubsup><mi>V</mi><mrow><mi>p</mi><mi>o</mi><mi>r</mi><mi>e</mi></mrow><mn>0</mn></msubsup></mrow>
Начальный поровый объем (до проникновения компонент бурового раствора)
:
<mrow><msubsup><mi>V</mi><mrow><mi>p</mi><mi>o</mi><mi>r</mi><mi>e</mi></mrow><mn>0</mn></msubsup><mo>=</mo><msub><mi>V</mi><mrow><mi>b</mi><mi>u</mi><mi>l</mi><mi>k</mi></mrow></msub><mo>−</mo><msub><mi>V</mi><mrow><mi>m</mi><mi>a</mi><mi>t</mi><mi>r</mi><mi>i</mi><mi>x</mi></mrow></msub><mtext>   (6)</mtext></mrow>
bulk
matrix
где V - общий объем образца породы; V - объем твердой матрицы ("скелета") в образце породы.
0
Начальная пористость образца породы φ:
<mrow><msub><mi>ϕ</mi><mn>0</mn></msub><mo>=</mo><mfrac><mrow><msubsup><mi>V</mi><mrow><mi>p</mi><mi>o</mi><mi>r</mi><mi>e</mi></mrow><mn>0</mn></msubsup></mrow><mrow><msub><mi>V</mi><mrow><mi>b</mi><mi>u</mi><mi>l</mi><mi>k</mi></mrow></msub></mrow></mfrac><mtext>    (7)</mtext></mrow>
<mrow><msubsup><mi>V</mi><mrow><mi>p</mi><mi>o</mi><mi>r</mi><mi>e</mi></mrow><mrow><mi>r</mi><mi>e</mi><mi>d</mi></mrow></msubsup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow>
В ходе фильтрации бурового раствора происходит захват частиц поровыми ловушками, что приводит к снижению порового объема
:
<mrow><msubsup><mi>V</mi><mrow><mi>p</mi><mi>o</mi><mi>r</mi><mi>e</mi></mrow><mrow><mi>r</mi><mi>e</mi><mi>d</mi></mrow></msubsup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>V</mi><mrow><mi>b</mi><mi>u</mi><mi>l</mi><mi>k</mi></mrow></msub><mo>−</mo><msub><mi>V</mi><mrow><mi>m</mi><mi>a</mi><mi>t</mi><mi>r</mi><mi>i</mi><mi>x</mi></mrow></msub><mo>−</mo><mstyle displaystyle="true"><munderover><mo>∑</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mn>2</mn></munderover><mrow><msubsup><mi>V</mi><mrow><msub><mi>p</mi><mi>i</mi></msub></mrow><mi>t</mi></msubsup></mrow></mstyle><mtext>   (8)</mtext></mrow>
<mrow><msubsup><mi>V</mi><mrow><msub><mi>p</mi><mi>i</mi></msub></mrow><mi>t</mi></msubsup></mrow>
здесь
- суммарный объем захваченных частиц i-го типа (i=1, 2) в образце породы.
Отметим, что поровый объем становится зависящим от времени согласно динамике захвата частиц.
Уменьшенная ("поврежденная") пористость φ(t):
<mrow><mi>ϕ</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msubsup><mi>V</mi><mrow><mi>p</mi><mi>o</mi><mi>r</mi><mi>e</mi></mrow><mrow><mi>r</mi><mi>e</mi><mi>d</mi></mrow></msubsup></mrow><mrow><msub><mi>V</mi><mrow><mi>t</mi><mi>o</mi><mi>t</mi><mi>a</mi><mi>l</mi></mrow></msub></mrow></mfrac><mtext>   (9)</mtext></mrow>
Λ
Насыщенность S порового пространства фазой Λ
<mrow><msub><mi>s</mi><mi>Λ</mi></msub><mo>=</mo><mfrac><mrow><msub><mi>V</mi><mi>Λ</mi></msub></mrow><mrow><msub><mi>V</mi><mrow><mi>p</mi><mi>o</mi><mi>r</mi><mi>e</mi></mrow></msub></mrow></mfrac><mtext>   (10)</mtext></mrow>
Λ
где V - объем фазы Λ в образце породы.
i
Объемная концентрация мобильных частиц i-го типа в буровом растворе C:
<mrow><msub><mi>C</mi><mi>i</mi></msub><mo>=</mo><mfrac><mrow><msubsup><mi>V</mi><mrow><mi>p</mi><mi>i</mi></mrow><mi>m</mi></msubsup></mrow><mrow><mstyle displaystyle="true"><munder><mo>∑</mo><mi>i</mi></munder><mrow><msubsup><mi>V</mi><mrow><mi>p</mi><mi>i</mi></mrow><mi>m</mi></msubsup><mo>+</mo><msub><mi>V</mi><mrow><mi>m</mi><mi>f</mi></mrow></msub></mrow></mstyle></mrow></mfrac><mo>≡</mo><mfrac><mrow><msubsup><mi>V</mi><mrow><mi>p</mi><mi>i</mi></mrow><mi>m</mi></msubsup></mrow><mrow><msub><mi>V</mi><mrow><mi>m</mi><mi>u</mi><mi>d</mi></mrow></msub></mrow></mfrac><mo>=</mo><mfrac><mrow><msub><mi>s</mi><mrow><mi>p</mi><mi>i</mi></mrow></msub></mrow><mrow><msub><mi>s</mi><mrow><mi>m</mi><mi>u</mi><mi>d</mi></mrow></msub></mrow></mfrac><mtext>   (11)</mtext></mrow>
<mrow><msub><mi>V</mi><mrow><mi>m</mi><mi>u</mi><mi>d</mi></mrow></msub><mo>=</mo><mstyle displaystyle="true"><munder><mo>∑</mo><mi>i</mi></munder><mrow><msubsup><mi>V</mi><mrow><mi>p</mi><mi>i</mi></mrow><mi>m</mi></msubsup><mo>+</mo><msub><mi>V</mi><mrow><mi>m</mi><mi>f</mi></mrow></msub></mrow></mstyle></mrow>
<mrow><msubsup><mi>V</mi><mrow><mi>p</mi><mi>i</mi></mrow><mi>m</mi></msubsup></mrow>
mf
pi
mud
где
- общий объем бурового раствора в образце породы; V - объем фильтрата бурового раствора в образце породы;
- суммарный объем мобильных частиц i-го типа в образце породы; s - насыщенность порового пространства мобильными частицами i-го типа; S - насыщенность порового пространства буровым раствором.
i
Объемная фракция захваченных частиц i-го типа в образце σ
<mrow><msub><mi>σ</mi><mi>i</mi></msub><mo>=</mo><mfrac><mrow><msubsup><mi>V</mi><mrow><mi>p</mi><mi>i</mi></mrow><mi>t</mi></msubsup></mrow><mrow><msub><mi>V</mi><mrow><mi>b</mi><mi>u</mi><mi>l</mi><mi>k</mi></mrow></msub></mrow></mfrac><mtext>   (12)</mtext></mrow>
<mrow><msubsup><mi>V</mi><mrow><mi>p</mi><mi>i</mi></mrow><mi>t</mi></msubsup></mrow>
bulk
где
суммарный объем захваченных частиц i-го типа в объеме V.
В соответствии с определениями начальной и уменьшенной пористости (7), (9) и объемной фракцией захваченных частиц (12) получаем следующее соотношение для уменьшенной ("поврежденной") пористости:
<mrow><mi>ϕ</mi><mrow><mo>(</mo><mrow><msub><mi>σ</mi><mi>Σ</mi></msub></mrow><mo>)</mo></mrow><mo>=</mo><msub><mi>ϕ</mi><mn>0</mn></msub><mo>−</mo><msub><mi>σ</mi><mi>Σ</mi></msub><mtext>   (13)</mtext></mrow>
Σ
1
2
где σ=σ+σ - суммарная фракция захваченных частиц 1-го и 2-го типа.
Система уравнений, описывающих проникновение компонент бурового раствора в пласт, содержит уравнения массового баланса для пластовой жидкости, фильтрата бурового раствора и частиц двух типов:
<mrow><mfrac><mo>∂</mo><mrow><mo>∂</mo><mi>t</mi></mrow></mfrac><mrow><mo>[</mo><mrow><mi>ϕ</mi><mrow><mo>(</mo><mrow><msub><mi>σ</mi><mi>Σ</mi></msub></mrow><mo>)</mo></mrow><msub><mi>ρ</mi><mrow><mi>r</mi><mi>f</mi></mrow></msub><msub><mi>s</mi><mrow><mi>r</mi><mi>f</mi></mrow></msub></mrow><mo>]</mo></mrow><mo>+</mo><mi>d</mi><mi>i</mi><mi>ν</mi><mrow><mo>[</mo><mrow><msub><mi>ρ</mi><mrow><mi>r</mi><mi>f</mi></mrow></msub><msub><mi>u</mi><mrow><mi>r</mi><mi>f</mi></mrow></msub></mrow><mo>]</mo></mrow><mo>=</mo><mn>0</mn><mtext>   (14)</mtext></mrow>
<mrow><mfrac><mo>∂</mo><mrow><mo>∂</mo><mi>t</mi></mrow></mfrac><mrow><mo>[</mo><mrow><mi>ϕ</mi><mrow><mo>(</mo><mrow><msub><mi>σ</mi><mi>Σ</mi></msub></mrow><mo>)</mo></mrow><msub><mi>ρ</mi><mrow><mi>m</mi><mi>f</mi></mrow></msub><msub><mi>s</mi><mrow><mi>m</mi><mi>f</mi></mrow></msub></mrow><mo>]</mo></mrow><mo>+</mo><mi>d</mi><mi>i</mi><mi>ν</mi><mrow><mo>[</mo><mrow><msub><mi>ρ</mi><mrow><mi>m</mi><mi>f</mi></mrow></msub><msub><mi>u</mi><mrow><mi>m</mi><mi>f</mi></mrow></msub></mrow><mo>]</mo></mrow><mo>=</mo><mn>0,</mn><mtext>   (15)</mtext></mrow>
<mrow><mfrac><mo>∂</mo><mrow><mo>∂</mo><mi>t</mi></mrow></mfrac><mrow><mo>[</mo><mrow><mi>ϕ</mi><mrow><mo>(</mo><mrow><msub><mi>σ</mi><mi>Σ</mi></msub></mrow><mo>)</mo></mrow><msub><mi>s</mi><mrow><mi>p</mi><mi>i</mi></mrow></msub></mrow><mo>]</mo></mrow><mo>+</mo><mi>d</mi><mi>i</mi><mi>ν</mi><mrow><mo>[</mo><mrow><msub><mi>C</mi><mi>i</mi></msub><msub><mi>u</mi><mrow><mi>m</mi><mi>u</mi><mi>d</mi></mrow></msub></mrow><mo>]</mo></mrow><mo>=</mo><mo>−</mo><msubsup><mi>q</mi><mrow><msub><mi>p</mi><mrow><mi>t</mi><mi>r</mi><mi>a</mi><mi>p</mi><mi>p</mi><mi>i</mi><mi>n</mi><mi>g</mi></mrow></msub></mrow><mi>i</mi></msubsup><mo>+</mo><msubsup><mi>q</mi><mrow><msub><mi>p</mi><mrow><msub><mrow /><mrow><mi>m</mi><mi>o</mi><mi>b</mi><mi>i</mi><mi>l</mi><mi>i</mi><mi>z</mi><mi>a</mi><mi>t</mi><mi>i</mi><mtext> </mtext><mi>o</mi><mi>n</mi></mrow></msub></mrow></msub></mrow><mi>i</mi></msubsup><mo>,</mo><mtext>  i</mtext><mo>=</mo><mtext>1</mtext><mtext>,2   (16)</mtext></mrow>
<mrow><mfrac><mo>∂</mo><mrow><mo>∂</mo><mi>t</mi></mrow></mfrac><msub><mi>σ</mi><mi>i</mi></msub><mo>=</mo><msubsup><mi>q</mi><mrow><msub><mi>p</mi><mrow><msub><mrow /><mrow><mi>t</mi><mi>r</mi><mi>a</mi><mi>p</mi><mi>p</mi><mi>i</mi><mi>n</mi><mi>g</mi></mrow></msub></mrow></msub></mrow><mi>i</mi></msubsup><mo>−</mo><msubsup><mi>q</mi><mrow><msub><mi>p</mi><mrow><msub><mrow /><mrow><mi>m</mi><mi>o</mi><mi>b</mi><mi>i</mi><mi>l</mi><mi>i</mi><mi>z</mi><mi>a</mi><mi>t</mi><mi>i</mi><mi>o</mi><mi>n</mi></mrow></msub></mrow></msub></mrow><mi>i</mi></msubsup><mo>,</mo><mtext>   i</mtext><mo>=</mo><mtext>1</mtext><mtext>,2   (17)</mtext></mrow>
rf
mf
rf
mf
pi
rf
mf
mud
<mrow><msubsup><mi>q</mi><mrow><msub><mi>p</mi><mrow><mi>t</mi><mi>r</mi><mi>a</mi><mi>p</mi><mi>p</mi><mi>i</mi><mi>n</mi><mi>g</mi></mrow></msub></mrow><mi>i</mi></msubsup></mrow>
<mrow><msubsup><mi>q</mi><mrow><msub><mi>p</mi><mrow><mi>m</mi><mi>o</mi><mi>b</mi><mi>i</mi><mi>l</mi><mi>i</mi><mi>z</mi><mi>a</mi><mi>t</mi><mi>i</mi><mi>o</mi><mi>n</mi></mrow></msub></mrow><mi>i</mi></msubsup></mrow>
где t - время; ρ - плотность пластовой жидкости; ρ - плотность фильтрата бурового раствора; s - насыщенность порового пространства пластовой жидкостью; s - насыщенность порового пространства фильтратом бурового раствора; s - насыщенность порового пространства частицами i-го типа; i - тип частицы; u - скорость фильтрации пластовой жидкости; u - скорость фильтрации фильтрата бурового раствора; u - скорость фильтрации бурового раствора;
- интенсивность захвата частиц i-го типа в породе пласта;
- интенсивность мобилизации ("отмыва") частиц i-го типа.
Скорости фильтрации бурового раствора и пластовой жидкости определяются согласно обобщенному закону Дарси
<mrow><msub><mi>u</mi><mrow><mi>m</mi><mi>u</mi><mi>d</mi></mrow></msub><mo>=</mo><mo>−</mo><mfrac><mrow><mi>k</mi><mrow><mo>(</mo><mrow><msub><mi>σ</mi><mi>Σ</mi></msub></mrow><mo>)</mo></mrow><msub><mi>k</mi><mrow><mi>r</mi><mo>,</mo><mi>m</mi><mi>u</mi><mi>d</mi></mrow></msub><mrow><mo>(</mo><mrow><msub><mi>s</mi><mrow><mi>m</mi><mi>u</mi><mi>d</mi></mrow></msub><mo>,</mo><msub><mi>σ</mi><mi>Σ</mi></msub></mrow><mo>)</mo></mrow></mrow><mrow><msub><mi>μ</mi><mrow><mi>m</mi><mi>u</mi><mi>d</mi></mrow></msub></mrow></mfrac><mo>∇</mo><msub><mi>p</mi><mrow><mi>m</mi><mi>u</mi><mi>d</mi></mrow></msub><mo>;</mo><msub><mrow><mtext>   u</mtext></mrow><mrow><mtext>rf</mtext></mrow></msub><mo>=</mo><mo>−</mo><mfrac><mrow><mi>k</mi><mrow><mo>(</mo><mrow><msub><mi>σ</mi><mi>Σ</mi></msub></mrow><mo>)</mo></mrow><msub><mi>k</mi><mrow><mi>r</mi><mo>,</mo><mi>r</mi><mi>f</mi></mrow></msub><mrow><mo>(</mo><mrow><msub><mi>s</mi><mrow><mi>r</mi><mi>f</mi></mrow></msub><mo>,</mo><msub><mi>σ</mi><mi>Σ</mi></msub></mrow><mo>)</mo></mrow></mrow><mrow><msub><mi>μ</mi><mrow><mi>r</mi><mi>f</mi></mrow></msub></mrow></mfrac><mo>∇</mo><msub><mi>p</mi><mrow><mi>r</mi><mi>f</mi></mrow></msub><mo>,</mo><mtext>  (18)</mtext></mrow>
r,mud
r,rf
mud
rf
mud
rf
Σ
1
2
<mrow><msub><mi>s</mi><mrow><mi>m</mi><mi>u</mi><mi>d</mi></mrow></msub><mo>=</mo><msub><mi>s</mi><mrow><mi>m</mi><mi>f</mi></mrow></msub><mo>+</mo><mstyle displaystyle="true"><munder><mo>∑</mo><mi>i</mi></munder><mrow><msub><mi>s</mi><mrow><mi>p</mi><mi>i</mi></mrow></msub></mrow></mstyle></mrow>
где k - проницаемость пласта; k - относительная фазовая проницаемость бурового раствора; k - относительная фазовая проницаемость пластовой жидкости; µ - вязкость бурового раствора; µ - вязкость пластовой жидкости; p - давление бурового раствора; p - давление пластовой жидкости;
- насыщенность порового пространства буровым раствором; C=C+C - суммарная объемная концентрация мобильных частиц 1-го и 2-го типа в буровом растворе.
с
В соответствии с определением капиллярного давления р:
<mrow><msub><mi>p</mi><mi>c</mi></msub><mo>=</mo><msub><mi>p</mi><mrow><mi>n</mi><mi>w</mi></mrow></msub><mo>−</mo><msub><mi>p</mi><mi>w</mi></msub><mo>,</mo><mtext>   (19)</mtext></mrow>
где нижний индекс "nw" соответствует несмачивающей фазе; нижний индекс "w" соответствует смачивающей фазе.
Скорость фильтрации фильтрата бурового раствора определяется следующим выражением:
<mrow><msub><mi>u</mi><mrow><mi>m</mi><mi>f</mi></mrow></msub><mo>=</mo><mrow><mo>(</mo><mrow><mn>1</mn><mo>−</mo><msub><mi>C</mi><mi>Σ</mi></msub></mrow><mo>)</mo></mrow><msub><mi>u</mi><mrow><mi>m</mi><mi>u</mi><mi>d</mi></mrow></msub><mo>.</mo><mtext>   (20)</mtext></mrow>
Для замыкания системы воспользуемся известными соотношениями, определяющими интенсивность захвата и мобилизации ("высвобождения") частиц в поровом пространстве, а также динамику изменения проницаемости (см. например Ives K.J., Pienvichitr V., 1965, Kinetics of filtration of dilute suspensions, Chemical Engineering Science. Vol.20. № 11; Herzig J.P., Leclerc D.M., Le Goff P., 1970, Flow of Suspensions through Porous Media - Application to Deep Filtration, Industrial and Engineering Chemistry. Vol.62. No.5):
<mrow><msubsup><mi>q</mi><mrow><msub><mi>p</mi><mrow><mi>t</mi><mi>r</mi><mi>a</mi><mi>p</mi><mi>p</mi><mi>i</mi><mi>n</mi><mi>g</mi></mrow></msub></mrow><mi>i</mi></msubsup><mo>=</mo><msubsup><mi>λ</mi><mrow><mi>c</mi><mi>i</mi></mrow><mn>0</mn></msubsup><msub><mi>C</mi><mrow><mi>p</mi><mi>i</mi></mrow></msub><msub><mi>u</mi><mi>m</mi></msub><mo>,</mo><mtext>    (21)</mtext></mrow>
<mrow><msubsup><mi>q</mi><mrow><msub><mi>p</mi><mrow><mi>m</mi><mi>o</mi><mi>b</mi><mi>i</mi><mi>l</mi><mi>i</mi><mi>z</mi><mi>a</mi><mi>t</mi><mi>i</mi><mi>o</mi><mi>n</mi></mrow></msub></mrow><mi>i</mi></msubsup><mo>=</mo><msubsup><mi>λ</mi><mrow><mi>m</mi><mi>i</mi></mrow><mn>0</mn></msubsup><msub><mi>σ</mi><mi>Σ</mi></msub><mrow><mo>(</mo><mrow><mi>U</mi><mo>−</mo><msub><mi>U</mi><mrow><mi>c</mi><mi>i</mi></mrow></msub></mrow><mo>)</mo></mrow><mo>,</mo><mtext>   (22)</mtext></mrow>
<mrow><mi>k</mi><mrow><mo>(</mo><mrow><msub><mi>σ</mi><mi>Σ</mi></msub></mrow><mo>)</mo></mrow><mo>=</mo><msub><mi>k</mi><mn>0</mn></msub><msup><mrow><mrow><mo>(</mo><mrow><mn>1</mn><mo>−</mo><mfrac><mrow><msub><mi>σ</mi><mi>Σ</mi></msub></mrow><mrow><mi>β</mi><msub><mi>ϕ</mi><mn>0</mn></msub></mrow></mfrac></mrow><mo>)</mo></mrow></mrow><mi>M</mi></msup><mo>;</mo><mtext>   (23)</mtext></mrow>
<mrow><msubsup><mi>λ</mi><mrow><mi>c</mi><mi>i</mi></mrow><mn>0</mn></msubsup></mrow>
<mrow><msubsup><mi>λ</mi><mrow><mi>m</mi><mi>i</mi></mrow><mn>0</mn></msubsup></mrow>
ci
здесь
- коэффициент захвата частиц i-го типа;
- коэффициент мобилизации частиц i-го типа; U - критическая скорость мобилизации частиц i-го типа; M - показатель, определяющий интенсивность снижения проницаемости породы с ростом доли удержанных в поровом пространстве частиц бурового раствора; β - эмпирический коэффициент, зависящий от структуры порового пространства и механизма захвата частиц в порах.
<mrow><msubsup><mi>λ</mi><mrow><mi>c</mi><mi>i</mi></mrow><mn>0</mn></msubsup></mrow>
<mrow><msubsup><mi>λ</mi><mrow><mi>m</mi><mi>i</mi></mrow><mn>0</mn></msubsup></mrow>
ci
Эмпирические параметры M,
,
, U, β математической модели внутренней фильтрационной корки для описания динамики изменения концентрации частиц бурового раствора в поровом пространстве образца керна неизвестны и могут быть определены путем настройки расчетов на данные лабораторных фильтрационных экспериментов и профиль концентрации проникших частиц бурового раствора (по данным рентгеновской компьютерной микротомографии образца керна).
Поскольку данные, регистрируемые в лабораторных фильтрационных экспериментах по закачке бурового раствора в образец керна, учитывают совместное влияния как внутренней, так и внешней фильтрационных корок (см., например, фиг.3), создается сцепленная математическая модель внешней и внутренней фильтрационных корок.
Связь ("сцепка") моделей внешней и внутренней фильтрационных корок реализована с помощью следующих предположений (см. фиг.1):
1. На начальном этапе фильтрации бурового раствора все содержащиеся в нем частицы проникают в образец керна (или пласт) и частично удерживаются в поровом пространстве, что приводит к формированию внутренней фильтрационной корки 2 на фиг.1 и, как следствие, снижению проницаемости породы.
crit
2. Если проницаемость керна в области входного торца (или пласта на границе со скважиной) 3 на фиг.1 снижается до некоторой критической величины k, то частицы больше не могут проникнуть в породу и отфильтровываются на поверхности керна (или стенке скважины), образуя внешнюю фильтрационную корку 1 на фиг.1.
При дискретизации задачи для численных расчетов данный критерий зависит от размера первой ячейки, начиная от торца керна. Этот размер, определяется крутизной падения профиля концентрации проникших частиц бурового раствора, полученных с помощью рентгеновской компьютерной микротомографии. Рекомендуется выбирать размер ячейки таким малым, чтобы на ее длине концентрация проникших частиц бурового раствора снижалась не более чем на 10%. Иным критерием является сходимость решения с заданной точностью при уменьшении размера сетки в области контакта пласт-скважина или вблизи торца керна.
crit
3. Критическое значение k предполагается равным проницаемости внешней фильтрационной корки при текущем перепаде давления.
Сформулированные предположения связывают модели внутренней и внешней корок.
mc0
mc0
Свойства внешней фильтрационной корки (проницаемость k, пористость φ, показатели сжимаемости ν и δ) уже определены исходя из данных проведенного эксперимента по продавливанию бурового раствора через фильтровальную бумагу или керамический диск (см. фиг.2).
<mrow><msubsup><mi>λ</mi><mrow><mi>c</mi><mi>i</mi></mrow><mn>0</mn></msubsup></mrow>
<mrow><msubsup><mi>λ</mi><mrow><mi>m</mi><mi>i</mi></mrow><mn>0</mn></msubsup></mrow>
ci
Фиксируя свойства внешней фильтрационной корки в сцепленной модели внешней и внутренней фильтрационных корок, находят параметры внутренней фильтрационной корки (M,
,
, U, β), при которых одновременно воспроизводятся данные лабораторного фильтрационного эксперимента (динамика перепада давления на образце, расход истекающей из образца жидкости и, дополнительно, профиль насыщенности фильтрата бурового раствора) и профиль концентрации проникших частиц бурового раствора.
Эмпирические параметры модели считаются определенными, когда расчеты соответствуют и данным лабораторного фильтрационного эксперимента и профилю концентрации проникших частиц бурового раствора и, дополнительно, профилю насыщенности фильтрата бурового раствора.
Полученные эмпирические параметры модели при данных концентрациях частиц в закачиваемом буровом растворе используют для моделирования динамики проникновения/удаления компонентов бурового раствора в призабойной зоне и сопутствующего изменения ее свойств в предполагаемых условиях бурения/очистки (репрессия/депрессия на пласт, количество проходок/очисток и т.п.) и временного графика.
Для этого используют ту же сцепленную модель внешней и внутренней фильтрационных корок, но в цилиндрической системе координат.
Задают начальные условия (до начала бурения): начальное давление, насыщение, проницаемость и пористость пласта.
w
i
w
В качестве граничного условия на границе скважина-пласт задается динамика давления p(r=R,t) и концентрация частиц в буровом растворе C(r=R,t) в соответствии с предполагаемой или реализованной программой бурения/очистки и типом буровых растворов.
Задается временной график бурения/очистки, планируемый или реализованный.
Для уменьшения времени расчета, не искажая физические особенности процесса проникновения компонент бурового раствора в пласт, используется трехзональная пространственная сетка (фиг.7):
w
d
- Равномерная, дополнительно измельченная сетка 14 - в зоне накопления захваченных частиц бурового раствора (внутренняя фильтрационная корка) от радиуса скважины R до R
d
inv
- Равномерная, измельченная сетка 15 - в зоне проникновения фильтрата бурового раствора от границы внутренней фильтрационной корки R до R
inv
- Неравномерная экспоненциальная сетка 16 - в зоне, незатронутой воздействием компонент бурового раствора от границы зоны проникновения фильтрата бурового раствора R до внешней границы сетки.
w
Внешняя фильтрационная корка, формирующаяся на стенке скважины, учитывается в модели как дополнительное, изменяющееся во времени, гидравлическое сопротивление в граничных условиях на границе пласта и скважины (r=R).
В качестве примера на фиг.8 приведена рассчитанная динамика потерь бурового раствора 17, а также рассчитанные профили распределения фильтрата бурового раствора 18 и поврежденной проницаемости 19 (в данном случае проведена нормировка на начальную проницаемость) в призабойной зоне пласта после 4-х часов бурения.
Описанный выше способ может быть реализован для множества различных буровых растворов, воздействующих на различные образцы кернов, соответствующих различным литотипам.
Согласно описанной выше последовательности путем настройки на данные лабораторных экспериментов по прокачке бурового раствора через образец керна (динамика перепада давления на образце, расход истекающей из образца жидкости и, дополнительно, профиль насыщенности фильтрата бурового раствора) и рентгеновской компьютерной микротомографии, могут быть определены параметры математических моделей внешней и внутренней фильтрационных корок для различных комбинаций разных типов воздействующих буровых растворов и образцов кернов, соответствующих различным литотипам.
mc0
mc0
ci
ci
<mrow><msubsup><mi>λ</mi><mrow><mi>c</mi><mi>i</mi></mrow><mn>0</mn></msubsup><mo>=</mo><msubsup><mi>λ</mi><mrow><mi>c</mi><mi>i</mi></mrow><mn>0</mn></msubsup><mrow><mo>(</mo><mrow><mi>k</mi><mo>,</mo><mi>ξ</mi><mo>,</mo><mi>r</mi></mrow><mo>)</mo></mrow></mrow>
<mrow><msubsup><mi>λ</mi><mrow><mi>m</mi><mi>i</mi></mrow><mn>0</mn></msubsup><mo>=</mo><msubsup><mi>λ</mi><mrow><mi>m</mi><mi>i</mi></mrow><mn>0</mn></msubsup><mrow><mo>(</mo><mrow><mi>k</mi><mo>,</mo><mi>ξ</mi><mo>,</mo><mi>r</mi></mrow><mo>)</mo></mrow></mrow>
На основе полученного набора параметров моделей для воздействия различных типов буровых растворов на образцы кернов различных литотипов могут быть построены зависимости ("корреляции") параметров моделей от свойств керна (в первую очередь, проницаемость k, и/или характерный размер пор ξ, и/или смачиваемость) и параметров компонентов бурового раствора (в первую очередь, размера кольматирующих частиц r): k=k(k,ξ,r),
,
, U=U(k,ξ,r), β=β(k,ξ,r).
Далее, при известной литологической стратификации разбуриваемых пластов, параметры модели рассчитывают согласно построенным корреляциям отдельно для каждого литологического пропластка и проводят моделирование динамики проникновения/удаления компонентов бурового раствора в призабойной зоне и сопутствующего изменения ее свойств в предполагаемых условиях бурения/очистки (репрессия/депрессия на пласт, количество проходок/очисток и т.п.) и временного графика.