Студопедия — Методика расчета фильтрационных сопротивлений при неустановившемся осесимметричном притоке жидкости (газа) к несовершенной скважине в неограниченном пласте.
Студопедия Главная Случайная страница Обратная связь

Разделы: Автомобили Астрономия Биология География Дом и сад Другие языки Другое Информатика История Культура Литература Логика Математика Медицина Металлургия Механика Образование Охрана труда Педагогика Политика Право Психология Религия Риторика Социология Спорт Строительство Технология Туризм Физика Философия Финансы Химия Черчение Экология Экономика Электроника

Методика расчета фильтрационных сопротивлений при неустановившемся осесимметричном притоке жидкости (газа) к несовершенной скважине в неограниченном пласте.






В Гл. 8 приведено решение (8.2.27) для неустановившегося притока сжимаемой жидкости к несовершенной скважине в однородно-анизотропном неограниченном по протяженности пласте по линейному закону фильтрации. Это решение может быть использовано для определения понижения давления на забое скважины после ее пуска, повышения давления на забое после остановки скважины, а также с целью анализа распределения потенциала (давления) в пласте во время работы скважина. Трудность использования указанного уравнения для практических целей заключается в его сложности. В приведенном виде решение не может быть попользовано для инженерных задач по двум причинам. Во-первых, функция (8.2.27) очень сложная и требует своего табулирования. Во-вторых, вид такой функции исключает возможность выделить время в качестве слагаемого в решении (8.2.27) и свести его к уравнению прямой, что необходимо для интерпретации кривых понижения или восстановления давления в скважинах традиционными методами.

Чтобы избежать упомянутых трудностей, можно поступить следующим образом. Как известно, в теории и практике гидродинамических исследований широко используется интегральная показательная функция. Несовершенство скважины при этом учитывается путем введения добавочных фильтрационных сопротивлений С1, взятых из решения задач для установившегося притока. Насколько такое допущение правильно, пока еще ни теоретически, ни экспериментально не доказано. Сведем функцию (8.2.27) к интегральной показательной и произведем расчет функции сопротивления (аналог добавочного сопротивления С1) при неустановившемся притоке, зависящей не только от геометрии пласта, но и от параметра Фурье (f 0).

 

9.3.1. Обоснование функции фильтрационного сопротивления. За основу возьмем функцию (8.2.28). Следовательно, для условий на забое интеграл (8.2.27) примет вид:

. (9.3.1)

С другой стороны, запишем следующее выражение:

. (9.3.2)

Решая совместно (9.3.1) и (9.3.2), находим

(9.3.3)

или

. (9.3.4)

Анализ показывает следующее. Равенство (9.3.4) содержит в себе интеграл вероятности Ф (х) и ряд

. (9.3.5)

Для расчета на ЭВМ функции сопротивления (9.3.4) требуется рассмотреть поведение (9.3.5), определить правило вычисления интеграла вероятности и выяснить условия его сходимости при заданных значениях f 0, , , t *.

Интеграл вероятности, определяемый как

(9.3.6)

для всех значений <1, вычислялся с помощью ряда

. (9.3.7)

Число используемых членов ряда определялось из соотношения

. (9.3.8)

Для значений <3,767 использовалось аппроксимирующее выражение:

, (9.3.9)

где а 0=1; а 1=0,07052308; а 2=0,04228201; а 3=0,009270527; а 4=0,0001520143; а 5=0,0002765672; а 6=0,000430638; для всех >3,767 принималось Ф (Х)=1. Область определения параметров ряда (9.3.5) лежит в пределах: ; . Так как значения параметров фиксированы, то ряд (9.3.5) зависит только от п. Члены рассматриваемого ряда представляют собой разность интегралов вероятности.

Интеграл вероятности – это ограниченная, монотонно возрастающая функция на интервале от 0 до ¥. Тогда при некотором п будет выполнено условие [27]

. (9.3.10)

Это будет допустимо, когда

. (9.3.11)

Разрешая неравенство (9.3.11) относительно п, находим число членов ряда (9.3.5). Интеграл (9.3.1) при =1 вычислялся методом Гаусса. Контроль точности вычислений проводился по значениям интегральной функции.

Расчет безразмерной депрессии и функции сопротивления производился соответственно по формулам (9.3.1) и (9.3.4) на БСМ-Минск-32, ОС “Фортран” в СибНИИНП объединения “Главтюменнефтегаз” В.И. Леоновым. Для расчета был разработан алгоритм. Параметры менялись в широком диапазоне: ; . Функция (9.3.1) затабулирована (см. Прил. 1 [28]); расчетные значения приведены в Прил. 2.

 

9.3.2. Сопоставление значений добавочных фильтрационных сопротивлений для установившегося и квазинеустановившегося притоков. Анализы показали, что при 0,10 для любого функция сопротивления не зависит от параметра Фурье f 0, а зависит только от геометрии пласта (рис. 9.9). При <0,01 с некоторого значения параметра Фурье функция безразмерного сопротивления также не зависит от , что свидетельствует о квазиустановившемся притоке. Как уже упоминалось ранее, в нефтегазопромысловой практике расчетов, связанных с неустановившемся притоком и учетом несовершенства скважин, принимаются значения фильтрационных сопротивлений, относящихся к установившемуся притоку в ограниченном пласте. Приемлемость этого допущения пока еще не обоснована.

Нами произведены сопоставления значений функции сопротивления R () для квазиустановившегося притока при f 0=10 созначениями добавочных фильтрационных сопротивлений при для установившегося притока в ограниченном пласте. Графические зависимости и (рис. 9.10), а также (рис. 9.11) показывают, что всюду значения функции сопротивления превосходят значения фильтрационных сопротивлений при установившемся притоке (R>C 1). Например, при =0,2 имеем: R / С 1=2,14 для =1000 и R / С 1=3,47 для =10; при =0,9 имеем: R / С 1=2,7 для R =1000 и R / С 1=27,4 для =10. С другой стороны, следует, отметить, что качественный характер поведения функции и аналогичный (см. рис. 9.9-9.11).

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

Подводя итоги изложенному, можно сделать следующие выводы:

– Численные значения функции для неустановившегося притока при любом вскрытии пласта всегда меньше численных значений фильтрационных сопротивлений при квазиустановившемся движении;

– Численные значения функции фильтрационного сопротивления для неустановившегося и квазиустановившегося притоков превосходят в несколько раз значения добавочных фильтрационных сопротивлений для установившегося притока в однородном ограниченном пласте; значений фильтрационных сопротивлений при квазиустановившемся движении;

– Численные значения функции фильтрационного сопротивления для неустановившегося и квазиустановившегося притоков превосходят в несколько раз значения добавочных фильтрационных сопротивлений для установившегося притока в однородном ограниченном пласте;

– Полученное аналитическое решение для неустановившегося притока сжимаемой жидкости (газа) к несовершенной скважине в бесконечном по протяженности пласте преобразовано в прямолинейную анаморфозу, которая позволяет эффективно производить интерпретацию кривых восстановления забойного давления;

– Предложенный здесь прием преобразования уравнения притока в интегральной форме (9.3.1), произведенные расчеты и табулирование функций и , позволят успешно реализовать для практических целей любое из решений, приведенных в Гл. 9.

 

Рис.9.9. Характеристика изменения функции фильтрационного сопротивления

при параметре

 

 

Рис. 9.10. Сопоставление значений функции сопротивления для квазиустановившегося притока при f 0=10-9 со значениями добавочных фильтрационных сопротивлений для установившегося притока

 

Рис. 9.11. Зависимость отношения при параметрах:

 







Дата добавления: 2015-08-12; просмотров: 705. Нарушение авторских прав; Мы поможем в написании вашей работы!



Вычисление основной дактилоскопической формулы Вычислением основной дактоформулы обычно занимается следователь. Для этого все десять пальцев разбиваются на пять пар...

Расчетные и графические задания Равновесный объем - это объем, определяемый равенством спроса и предложения...

Кардиналистский и ординалистский подходы Кардиналистский (количественный подход) к анализу полезности основан на представлении о возможности измерения различных благ в условных единицах полезности...

Обзор компонентов Multisim Компоненты – это основа любой схемы, это все элементы, из которых она состоит. Multisim оперирует с двумя категориями...

Влияние первой русской революции 1905-1907 гг. на Казахстан. Революция в России (1905-1907 гг.), дала первый толчок политическому пробуждению трудящихся Казахстана, развитию национально-освободительного рабочего движения против гнета. В Казахстане, находившемся далеко от политических центров Российской империи...

Виды сухожильных швов После выделения культи сухожилия и эвакуации гематомы приступают к восстановлению целостности сухожилия...

КОНСТРУКЦИЯ КОЛЕСНОЙ ПАРЫ ВАГОНА Тип колёсной пары определяется типом оси и диаметром колес. Согласно ГОСТ 4835-2006* устанавливаются типы колесных пар для грузовых вагонов с осями РУ1Ш и РВ2Ш и колесами диаметром по кругу катания 957 мм. Номинальный диаметр колеса – 950 мм...

ТЕРМОДИНАМИКА БИОЛОГИЧЕСКИХ СИСТЕМ. 1. Особенности термодинамического метода изучения биологических систем. Основные понятия термодинамики. Термодинамикой называется раздел физики...

Травматическая окклюзия и ее клинические признаки При пародонтите и парадонтозе резистентность тканей пародонта падает...

Подкожное введение сывороток по методу Безредки. С целью предупреждения развития анафилактического шока и других аллергических реак­ций при введении иммунных сывороток используют метод Безредки для определения реакции больного на введение сыворотки...

Studopedia.info - Студопедия - 2014-2024 год . (0.01 сек.) русская версия | украинская версия