Отпечатано с готовых диапозитивов
В. А. Евдокимова
И. Н. Кочина
Сборник задач
По подземной
Гидравлике
Учебное пособие для вузов
Издание второе, стереотипное
Перепечатка с издания 1979 г.
Первое издание допущено Министерством высшего
и среднего специального образования СССР в
качестве учебного пособия для студентов высших
учебных заведений, обучающихся по специальности
«Технология и комплексная механизация разработки нефтяных и газовых месторождений»
АльянС
Москва 2007
УДК 622.276 (0.75.8)
ББК
Е15
Евдокимова В.А., Кочина И.Н.
Е15 Сборник задач по подземной гидравлике:
Учебное пособие для вузов. — 2-е изд., стереотипное. Перепечатка с издания 1979 г. - М.: ООО ИД «Альянс», 2007. - 168 с.
Приведенные в учебном пособии ладами, можно использовать при проектировании разработки нефтяных и газовых месторождений. Первое, издание выпущено в 1979 году издательством «Недра». Учебное пособие рассчитано на студентов нефтяных вузов и факультетов.
ISBN 978-5-903034-13-0
© ООО ИД «Альянс», 2007
ЕВДОКИМОВА Вера Алексеевна
КОЧИНА Ираида Николаевна
Сборник задач по подземной гидравлике
Учебное пособие для вузов
Издание второе, стереотипное
Перепечатка с издания 1979 г.
Подписано в печать 27.10.06. Формат 60x90/16.
Гарнитура Литературная. Печати офсетная Исч. л. 10,5.
Тираж 1000 экз. Заказ 61061
ООО Издательский дом «Альянс»
105120, Москва, ул. Сергия Радонежского, д.9, стр. 5
Тел./факс (499) 973-06-80, 973-09-41
973-17-82, 973-17-96, 973-18-56
Отпечатано с готовых диапозитивов
Г. Саратов, ул. Волжская, д.28
ЗАО "Типография "Полиграфист"
ПРЕДИСЛОВИЕ
В сборник включены задачи, которые можно использовать при проектировании нефтяных и газовых месторождений, решении различных проблем гидротехники, инженерной геологии, гидрогеологии, ирригации и горного дела. Решение многих задач подземной гидравлики полезно также при расчете искусственных фильтров различных конструкций, пористых катализаторов и т. д.
При составлении сборника задач авторы использовали многолетний опыт преподавания курса «Подземная гидравлика» в Московском институте нефтехимической и газовой промышленности им. акад. И. М. Губкина. В сборник, в основном, вошли задачи, которые предлагались студентам на практических занятиях.
Настоящее пособие предназначено также для студентов специальностей «Геология и разведка нефтяных и газовых месторождений» и «Экономика и организация нефтяной и газовой промышленности».
Сборник задач состоит из 15 глав. К каждой главе дается краткая теория. Ко всем задачам имеются ответы. Типовые и наиболее сложные задачи приведены с решениями. В решениях некоторых задач даются выводы формул, отсутствующие в учебной литературе.
В сборник входят задачи на определение фильтрационных характеристик пластов, расчет производительности нефтяных игазовых эксплуатационных и нагнетательных скважин в однородных и неоднородных по проницаемости пористых пластах, а также в деформируемых трещиноватых пластах, учет интерференции скважин (совершенных и несовершенных), расчет продвижения водонефтяного контакта, определение высоты подъема конуса подошвенной воды при эксплуатации нефтяных или газовых месторождений с подошвенной водой, определение дебитов и распределения давления при движении газированной жидкости в пористой среде, изменение дебитов и давлений при нестационарном движении упругой жидкости и газа в деформируемой пористой среде, вытеснение нефти водой по теории Баклея — Леверетта и др.
ОСНОВНЫЕ ПОНЯТИЯ ТЕОРИИ ФИЛЬТРАЦИИ
Фильтрация
Фильтрацией называется движение жидкостей, газов и их смесей в пористых и трещиноватых средах, т. е. в твердых телах, пронизанных системой сообщающихся между собой пор и микротрещин.
Фильтрация жидкостей и газов по сравнению с движением в трубах и каналах обладает некоторыми специфическими особенностями.
Фильтрация происходит по чрезвычайно малым в поперечных размерах норовым каналам при очень малых скоростях движения жидкостей. Силы трения при движении жидкости в пористой среде очень велики, так как площади соприкосновения жидкостилощади соприкосновения жидкое с твердыми частицами огромны.
Пористая среда характеризуется коэффициентами пористости и проснетпости.
Коэффициент пористости тесть отношение объема пор (τпор) ко всему объему пористой среды (τ)
(I.1)
Под пористостью в теории фильтрации понимается активная пористость, которая учитывает только те поры и микротрещины, которые соединены между собой и через которые может фильтроваться жидкость.
Коэффициентом просветности п называется отношение площади просветов (ωпросв) в данном сечении пористой среды ко всей площади этого сечения (ω)
(I.2)
Можно показать, что среднее по длине пласта значение просветности равно пористости, т. е.
(I.3)
поэтому среднее значение площади просветов
Упрощенной моделью пористой среды является модель фиктивного грунта. Фиктивный грунт состоит из шариков одного диаметра, уложенных определенным образом. Основным элементом (основной ячейкой) фиктивного грунта является ромбоэдр, который получится, если принять центры восьми соприкасающихся частиц за вершины углов (рис. 1), В зависимости от острого угла θ боковой грани ромбоэдра укладка шаров более или менее плотная.
Угол θ изменяется в пределах от 60° до 90°. Углу θ = 60° соответствует наиболее плотная укладка шаров, углу θ = 90° — наиболее свободная.
Пористость фиктивного грунта определяется по формуле Ч. Слихтера
(I.4)
из которой следует, что пористость зависит не от диаметра частиц, а лишь от их взаимного расположения, которое определяется углом θ.
Чтобы формулы для фиктивного грунта можно было применять для естественного грунта, нужно заменить реальный грунт эквивалентным ему фиктивным, который должен иметь такое же гидравлическое сопротивление, как у естественного грунта.
Диаметр частиц такого фиктивного грунта называется эффективным диаметром (d3). Эффективный диаметр определяется в результате механического анализа грунта. Его просеивают через набросит с различной площадью отверстий и, таким образом, разделяют на фракции. За средний диаметр каждой фракции принимают среднее арифметическое крайних диаметров, т.е.
Затем строят кривую механического (фракционного) состава грунта, откладывая по оси абсцисс средние диаметры фракций di а но оси ординат — сумму масс фракций Δg1+ Δg2 + + ... + Δgi в % от общей массы.
Последняя точка кривой имеет абсциссу, равную dn, и ординату Δg1+ Δg2+... + Δgn =100% (рис. 2).
Существует много способов определения эффективного диаметра. По способу А. Газена dэ определяется по кривой механического состава. За эффективный принимается такой диаметр шарообразной частицы, который соответствует сумме масс всех фракций, начиная от нуля и кончая этим диаметром, равной 10%. Надо найти, кроме того, диаметр d0, который соответствует сумме масс фракций,, равной 60%. Коэффициент однородности do/dэ должен быть не более 5 (do/dэ ≤ 5) и dэ должен лежать в пределах от 0,1 до 3 мм.
По способу Крюгера — Цункера используют данные механического анализа грунта и определяют dэ по формуле
(I.5)
Скоростью фильтрации w называется отношение объемного расхода жидкости к площади поперечного сечения пласта, нормального к направлению движения жидкости
(I.6)
Скорость фильтрации представляет собой фиктивную скорость, с которой двигалась бы жидкость, если бы пористая среда отсутствовала (m=1).
Средняя скорость движения жидкости v равна отношению объемного расхода к площади просветов ωпросв (живому сечению потока)
(I.7)
Скорость фильтрации и средняя скорость движения связаны соотношением
Критерий Рейнольдса
Подобно тому, как в трубной гидравлике критерием режима движения служит число Рейнольдса
Re = υdρ/μ (II.1)
в теории фильтрации вводится безразмерный параметр
Re = (II. 2)
где и — некоторая характерная скорость; — линейный параметр, характеризующий среднее сечение поровых каналов; ρ — плотность жидкости; μ — динамический коэффициент вязкости.
Скорость фильтрации, при которой нарушается закон Дарси, называется критической скоростью фильтрации (wкр).
Однако нарушение линейного закона фильтрации еще не означает перехода от ламинарного движения к турбулентному. Закон Дарси нарушается вследствие того, что силы инерции, возникающие в жидкости за счет извилистости каналов и изменения площади их поперечных сечений, становятся при w > wкр соизмеримыми с силами трения.
В трубной гидравлике значение Re, при котором происходит смена режимов, равно Reкр= 2320, в теории фильтрации закон Дарси имеет место при значении безразмерного параметра Re, меньшего критического (Reкр), которое устанавливается из опыта.
Впервые число Рейнольдса для фильтрации жидкости было введено Н. Н. Павловским в виде
(II.3)
т. е. за характерную скорость была взята скорость фильтрации w, а линейный параметр представлен выражением
(II.4)
Критические значения Re по Павловскому заключены в интервале
Reкр = 7,5 9.
В. Н. Щелкачев предложил взять за линейный параметр выражение, пропорциональное корню квадратному из коэффициента проницаемости,
(II.5)
Число Рейнольдса по В. Н. Щелкачеву имеет вид
(II.6)
a критические значения лежат в интервале
l ≤ Reкр ≤ 12.
По М. Д. Миллионщикову за характерную скорость взята редкая скорость движения жидкости
υ = w/m,
аза линейный параметр — выражение т. е.
(II.7)
Если вычисленное по одной из формул (II.З), (II.6), (II.7) значение числа Re оказывается меньше нижнего критического значения Reкр, то закон Дарси справедлив, если Re больше верхнего значения Reкр, то закон Дарси заведомо нарушен.
Широкий диапазон изменения Reкр объясняется тем, что в формулы для числа Re входят параметры k и т, которые не полностью характеризуют микроструктуру породы. Как следует из опытов, для каждой горной породы можно указать более узкий диапазон значений Reкр [16].
Определение режима фильтрации жидкостей и газов имеет большое практическое значение, ибо без знания закона фильтрации в пласте нельзя правильно рассчитать дебиты скважин, распределение давления в пласте, а также невозможно определение параметров пласта (k, h, m и др.) по данным исследования нефтяных и газовых скважин.
Формула Дюпюи
При плоскорадиальном движении векторы скорости фильтрации направлены по радиусам к оси скважины, поэтому давление и скорость фильтрации зависят только от одной координаты r. При этом во всех горизонтальных плоскостях поле скоростей и давлений будет одинаковым.
Примером плоскорадиального фильтрационного потока является приток к гидродинамически совершенной скважине, покрывшей горизонтальный пласт бесконечной протяженности на всю мощность h и сообщающейся с пластом через полностью открытую боковую поверхность цилиндра, отделяющую cтвол скважины от продуктивного пласта.
Поток будет также плоскорадиальным при притоке к совершенной скважине радиуса rс (или оттоке от скважины), расположенной в центре ограниченного горизонтального цилиндрического пласта мощностью h и радиусом RK (рис. 7).
Если на внешней границе пласта, совпадающей с контуром литания, поддерживается постоянное давление рк, а на забое скважины постоянное давление pс, пласт однороден по пористости и проницаемости, фильтрация происходит по закону Дарси, то объемный дебит скважины определится по формуле Дюпюи:
(III.4)
где μ — динамический коэффициент вязкости.
Закон распределения давления определяется по одной из формул:
(III.5)
либо
(III.6)
либо
(III.7)
Линия р=р(r) называется депрессионной кривой давления. Характерно, что при приближении к скважине градиенты давления и скорости фильтрации резко возрастают. При построении карты изобар следует учитывать, что радиусы изобар изменяются в геометрической прогрессии, в то время, как давление на изобарах изменяется арифметической прогрессии.
Индикаторная линия — зависимость дебита скважины от депрессии Δр = рк—рс, при притоке к скважине в условиях справедливости закона Дарси представляет собой прямую линию, определяемую уравнением Q=KΔp.
Коэффициент продуктивности
(III.8)
численно равен дебиту при депрессии, равной единице.
Закон движения частиц вдоль линии тока, если при t = 0 частица находилась в точке с координатой r = r0, описывается уравнением
(III.9)
или
(III.9a)
Средневзвешенное по объему порового пространства Ω пластовое давление
где
Подставляя выражение для p (III.5), выполняя интегрирование и пренебрегая всеми членами, содержащими rc2, получим
(III.11)
Закон распределения давления и формула дебита при нарушении закона Дарси при притоке к совершенной скважине получаются из двучленной формулы
(III.12)
Подставляя выражение для скорости фильтрации
w = Q/2πrh
в(III.12) и разделяя переменные, получим
(III.13)
Интегрируя по р в пределах от рс до рк и по r в пределах от rс до Rk будем иметь
(III.14)
Решая полученное квадратное уравнение, находим дебит скважины Q. Интегрируя (III.13) по р в пределах от р до рк и по r в пределах от r до Rк, найдем закон распределения давления
(III.15)
Как видно из (III.14), индикаторная линия при нарушении закона Дарси является параболой.
Если фильтрация происходит по закону Краснопольского, то дебит определяется по формуле
(III. 16)
Жидкости по закону Дарси
Фильтрационный поток называется радиально-сферическим, «если векторы скорости фильтрации направлены в пространстве по прямым, радиально сходящимся к одной точке (или расходящимся от нее).
Благодаря центральной симметрии давление и скорость фильтрации зависят и в этом случае только от одной координаты r, отсчитываемой от центра (рис. 8). Примером потока, весьма близкого радиально-сферическому, является приток жидкости к гидродинамически несовершенной скважине малого диаметра, едва вскрывшей непроницаемую горизонтальную кровлю однородного пласта большой мощности (теоретически бесконечной). Если на забое скважины, представленной в виде полусферы радиуса rс, поддерживается постоянное приведенное давление, , а на достаточно большом расстоянии от скважины, на полусферической поверхности радиуса Rк сохраняется постоянное давление и фильтрация в однородном пласте происходит по закону Дарси, то объемный дебит скважины определяется по формуле
(III.17)
Приведенное давление в любой точке пласта определяется по формуле
(III.18)
азакон движения частиц вдоль линии тока от точки с координатой r0 до точки с координатой r описывается уравнением
(III.19)
Задача 20
Определить дебит дренажной галереи шириной B = 100 м, если мощность пласта h=10 м, расстояние до контура питания l = 10 км, коэффициент проницаемости пласта k=1 Д, динамический коэффициент вязкости жидкости μ = l сП, давление на контуре питания pк = 9,8 МПа (100 кгс/см2) и давление в галерее pг = 7,35 МПа (75 кгс/см2). Движение жидкости напорное, подчиняется закону Дарси.
Ответ: Q = 21,6 м3/сут.
Задача 21
Определить коэффициент проницаемости пласта (в различных системах единиц), если известно, что в пласте происходит одномерное, прямолинейно-параллельное установившееся движение однородной жидкости по закону Дарси. Гидравлический уклон i = 0,03, ширина галереи В = 500 м, мощность пласта h=6 м, плотность жидкости ρ = 850 кг/м3, динамический коэффициент вязкости μ = 5 сП и дебит галереи Q = 30 м3/сут.
Ответ: k=2,27 Д=32∙10-8 см2=2,32∙10-12 м2.
Задача 22
Показать графически распределение давления и найти градиент давления при прямолинейно-параллельном движении в пласте несжимаемой жидкости по линейному закону фильтрации, используя следующие данные: длина пласта lк = 5 км, мощность пласта h=10 м, ширина галереи B = 300 м, коэффициент проницаемости пласта k = 0,8 Д, давление в галерее рг = 2,94 МПа (30 кгс/см2), динамический коэффициент вязкости жидкости μ = 4 сП, дебит галереи Q = 30 м3/сут.
Ответ: p = 5,78 - 0,0568∙10-2х (х в м, р в МПа), - dp/dx = 0,0568∙10-2 МПа/м.
Задача 23
Определить дебит нефтяной скважины (в т/сут) в случае установившейся плоскорадиальной фильтрации жидкости по закону Дарси, если известно, что давление на контуре питания рк =9,8 МПа (100 кгс/см2), давление на забое скважины рс =7,35 МПа (75 кгс/см2), коэффициент проницаемости пласта k = 0,5 Д, мощность пласта h = 15 м, диаметр скважины Dc=24,8 см, радиус контура питания Rк=10 км, динамический коэффициент вязкости жидкости μ = 6 мПа∙с и плотность жидкости р = 850 кг/м3.
Ответ: Qm= 127 т/сут.
Задача 24
Определить давление на расстоянии 10 и 100 м от оси скважины при плоскорадиальном установившемся движении несжимаемой жидкости по линейному закону фильтрации, считая, что коэффициент проницаемости пласта k = 0,5 Д, мощность пласта h = 10 м, давление на забое скважины рс = = 7,84 МПа (80 кгс/см2), радиус скважины rс== 12,4см, динамический коэффициент вязкости нефти μ = 4∙10-3 кг/м∙с, плотность нефти ρ = 870 кг/м3 и массовый дебит скважины Qm = 200 т/сут.
Ответ: p1=9,28 МПа; p2 = 10,06 МПа.
Задача 25
Построить индикаторную линию (зависимость дебита Q от перепада давления Δр = рк —рс), имеющуюся при установившейся плоскорадиальной фильтрации жидкости по линейному закону, если известно, что давление на контуре питания pк = 8,82 МПа (90 кгс/см2), коэффициент проницаемости пласта k = 600 мД, мощность пласта h=10 м, диаметр скважины Dc = = 24,8 см, расстояние от оси скважины до контура питания Rk=10 км и динамический коэффициент вязкости нефти μ = 5 мПа∙с.
Ответ: индикаторная линия — прямая, описываемая уравнением Q = 5,77 Δр (Q в м3/сут, Δр в кгс/см2).
Задача 26
Определить коэффициент гидропроводности пласта kh/μ по данным о коэффициенте продуктивности скважины. Известно, что фильтрация происходит по закону Дарси, коэффициент продуктивности K=18 т/сут (кгс/см2), среднее расстояние между скважинами 2σ = 1400 м, плотность ρ=925 кг/м3, радиус скважины rс= 0,1 м.
Ответ:kh/μ = 3,18∙10-9 м4∙с/кг (318 Д∙см/сП).
Задача 27
Определить средневзвешенное по объему пластовое давление, если известно, что давление на контуре питания рк = 9,8 МПа (100 кгс/см2), давление на забое возмущающей скважины pс = 7,84 МПа (80 кгс/см2), расстояние до контура питания Rк = 25 км, радиус скважины rс = 10 см. В пласте имеет место установившееся плоскорадиальное движение несжимаемой жидкости по закону Дарси.
Ответ: р = 9,72 МПа (99,19 кгс/см2).
Задача 28
Определить относительное понижение sp/s= (Hк—Н)/(Нк—Hс) пьезометрического уровня в реагирующих скважинах, расположенных от возмущающей скважины на расстояниях 1 м, 100 м, 1 км, 10 км. Движение жидкости установившееся плоскорадиальное по закону Дарси. Радиус скважины rс = 0,1 м, расстояние до контура питания Rк=100 км.
Ответ: sp/s равно соответственно 0,83; 0,50; 0,33; 0,167.
Задача 29
Определить время отбора нефти из призабойной зоны скважины радиусом r0=100 м, если мощность пласта h=10 м, коэффициент пористости пласта m = 20%, массовый дебит нефти Qm =40 т/сут, плотность ее ρ= 920 кг/м3, rс = 0,1 м.
Ответ: Т = 1440 сут.
Задача 30
Определить время t, за которое частица жидкости подойдет к стенке скважины с расстояния r0 = 200 м, если коэффициент проницаемости пласта k=1 Д, динамический коэффициент вязкости нефти μ = 5 сП, депрессия во всем пласте радиусом Rк = 1 км составляет рк—рс= 10 кгс/см2; мощность пласта h=10м, коэффициент пористости пласта m = 15%, радиус скважины rc = 10 см.
Ответ: t = 1600 сут.
Задача 31
Как изменится дебит скважины Q при увеличении радиуса скважины вдвое?
1. Движение происходит по линейному закону фильтрации.
2. Фильтрация происходит по закону Краснопольского.
Начальный радиус скважины rс = 0,1 м. Расстояние до контуpa питания Rк = 5 км.
Ответ: 1) Q’: Q=l,07; 2) Q’: Q= 1,41, т. е. при движении жидкости по линейному закону фильтрации влияние изменения радиуса скважины менее интенсивно, чем при движении по закону Краснопольского.
Задача 32
Найти изменение перепада давления Δр при увеличении радиуса скважины вдвое, при котором дебит остается прежним. Рассмотреть два случая, как в предыдущей задаче. Начальный радиус скважины rс = 0,1 м, расстояние до контура питания Rк= 1 км.
Ответ: 1) Δр'/Δр = 0,925, 2) Δр'/Δр = 0,5.
Задача 33
Во сколько раз необходимо увеличить радиус скважины, чтобы дебит ее при прочих равных условиях удвоился?
1) Движение жидкости происходит по закону Дарси.
2) Жидкость фильтруется по закону Краснопольского. Начальный радиус скважины rс = 0,1 м. Расстояние до контура питания Rк= 1 км.
Ответ: 1) n=100, r'c=10 м; 2) n = 4, r'с = 0,4 м.
Задача 34
Скважина радиусом rс=10 см расположена в центре кругового пласта радиусом Rк = 350 м. Коэффициент проницаемости пласта k = 0,8 Д, мощность h=12м, динамический коэффициент вязкости нефти μ=5 сП. Определить дебит скважины, считая, что залежь по контуру радиуса Rк частично непроницаема (рис. 9). Контур питания представляет собой в плане дугу окружности радиусом Rк с центральным углом α = 120°. Давление на контуре питания рк = 27,9 МПа (285 кгс/см2), давление на забое скважины рс = 7,84 МПа (80 кгс/см2).
Решение. Задачу можно свести к плоскорадиальной, если в формуле Дюпюи за контурное давление принять средневзвешенное по длине окружности давление рк
МПа,
= 2,22-10~3 м3/с - 192 м3/сут.
Задача 35
Сколько жидкости следует закачивать в пласт в единицу времени через нагнетательную скважину, если необходимо, чтобы давление в скважине поддерживалось в процессе закачки на Δр=1,47 МПа (15 кгс/см2) выше давления, установившегося в пласте на расстоянии r = 2 км от скважины? Имеет место закон Дарси. Динамический коэффициент вязкости μ =1 сП, коэффициент проницаемости пласта k = l50 мД, мощность пласта h =10 м, радиус скважины rс = 10 см.
Ответ: Q=123 м3/сут.
Задача 36
Определить приведенное давление в точках, отстоящих на r = 20 м, 10 м, 5 м, 1,5 м, 1 м от центра забоя скважины, вскрывшей пласт бесконечной мощности на величину b = 0,5 м. На расстоянии Rк=1000м приведенное давление рк*= 9,8 МПа (100 кгс/см2), на забое скважины рс*=7,35 МПа (75кгс/см2), рллиус скважины r'c = 12,4см. Фильтрация к скважине происходит по закону Дарси.
Указание.Представляя забой скважины в виде полусферы, Равновеликой по площади забою действительной скважины, определить радиус полусферы rc(2πr'сb=2πr2c).
Ответ: соответственно р* = 9,77; 9,74; 9,68; 9,39; 9,19 МПа.
Задача 37
Скважина вскрывает пласт бесконечно большой мощности на небольшую глубину. Считая движение радиальносферическим, определить время перемещения частиц жидкости вдоль линий тока от точки с координатой r0=100 м до точки с координатой r = 5 м. Скважина эксплуатируется с постоянным дебитом Q = 120м3/сут, коэффициент пористости пласта m = 15%.
Ответ: t = 7,15 лет.
СВЯЗЬ ПЛОСКОЙ ЗАДАЧИ ТЕОРИИ
КОМПЛЕКСНОГО ПЕРЕМЕННОГО
В самом общем случае давление и скорость фильтрации зависят от трех координат точки в пласте. Если давление и скорость фильтрации зависят только от двух координат, то в каждой плоскости, перпендикулярной к третьей оси, поле скоростей и давлений будет одинаковым. В этом случае фильтрационный поток называется плоским. Плоский фильтрационный поток имеет место при работе одной или нескольких гидродинамически совершенных (эксплуатационных и нагнетательных) скважин в однородном горизонтальном пласте постоянной мощности. Именно такие потоки будут рассмотрены в настоящем разделе.
Принцип суперпозиции
Назовем точечным стоком па плоскости точку, поглощающую жидкость. Сток можно рассматривать как гидродинамически совершенную эксплуатационную скважину бесконечно малого радиуса в пласте единичной мощности. Точечный источник — это точка, выделяющая жидкость (аналог нагнетательной скважины). Заменяя источники и стоки скважинами конечного диаметра, мы практически не допускаем никакой ошибки, поэтому будем в дальнейшем отождествлять скважины с источниками и стоками.
При работе в бесконечном пласте одной скважины-стока фильтрация будет плоскорадиальной и давление в точке на расстоянии rот центра скважины определяется по формуле
(IV.1)
где q=Q/h — дебит скважины-стока, приходящийся на единицу мощности пласта; С — постоянная интегрирования.
Назовем потенциалом скорости фильтрации Ф выражение Ф = kp/μ. Переходя от давления к потенциалу, получим значение потенциала в точке на расстоянии rот центра скважины
(IV.2)
Дебиту источника (нагнетательной скважины) приписывается знак минус.
При совместной работе в пласте нескольких скважин результирующий потенциал в любой точке пласта М равен алгебраической сумме потенциалов Ф1, Ф2, ... , обусловленных работой каждой отдельной скважины.
(IV.3)
Скорости фильтрации при этом складываются геометрически (рис. 10, а, б). Это называется принципом суперпозиции или сложения течений.
Используя принцип суперпозиции, можно приближенно рассчитывать дебиты или забойные потенциалы (а следовательно, и забойные давления) для группы скважин, работающих в пласте с весьма удаленным контуром питания. Потенциал Фк на контуре питания считается известным, а расстояние от контура питания до всех скважин — одно и то же и приблизительно равно Rк.
Помещая мысленно точку М последовательно на забой каждой скважины, где Фм = Фсi, получим из общего уравнения (IV.3) систему п уравнений (п — число скважин). Постоянная интегрирования находится из условия на контуре питания. Окончательно система уравнений для определения дебитов или забойных потенциалов примет вид
(IV.4)
…………………………………………………………
здесь rij — расстояние между центрами i-той и j-той скважин.
Принцип суперпозиции можно использовать, если скважины работают в пласте, ограниченном контуром питания той или мной формы, или непроницаемыми границами (линии выклинивания, сбросы), но для выполнения тех или иных условий на границах приходится вводить фиктивные скважины за пределами пласта, которые создают в совокупности с реальными скважинами необходимые условия на границах.
При этом задача сводится к рассмотрению одновременной работы реальных и фиктивных скважин в неограниченном пласте. Этот метод называется методом отображения источников и стоков. Он широко применяется не только в подземной гидравлике и гидродинамике, но и при решении задач теории электричества, магнетизма и электропроводности.
Так, если эксплуатационная скважина находится в пласте с прямолинейным контуром питания на расстоянии а от контура, то ее надо зеркально отобразить относительно контура, т. е. поместить фиктивную скважину с другой стороны от контура на расстоянии а (рис. 11) и считать ее дебит отрицательным (скважина — источник). При этом потенциал в любой точке М равен
на контуре питания r1 = r2 и Ф = С=Фк, а дебит скважины определяется по формуле
(IV.5)
Метод отображения источников и стоков используется также в задачах 52, 53 для нахождения дебита скважины, работающей в пласте, ограниченном пересекающимися прямолинейными непроницаемыми границами. При помощи этого метода можно определить дебит скважины, эксцентрично расположенной в круговом пласте
(IV.6)
где δ — расстояние от центра скважины до центра кругового пласта (эксцентриситет).
Интерференция скважин
Дебит каждой скважины бесконечной цепочки, расположенной на расстоянии L от прямолинейного контура питания (рис. 12), выражается формулой
(IV.7)
где σ - половина расстояния между скважинами. Если L ≥ σ, то приближенно можно принять, что
и тогда
(IV.8)
Дебит одной скважины кольцевой батареи, состоящей из п скважин, в круговом пласте радиуса Rк (рис. 13) имеет вид
где R1 — радиус батареи; rс — радиус скважин.
Если число скважин батареи велико (больше пяти или шести), то (R1/Rк)2n ≤ 1 и этим выражением можно пренебречь по сравнению с единицей; если, кроме того, заменить
то получим приближенную формулу
(IV. 10)
Формулы (IV.7) и (IV.9) можно вывести, используя метод отображения.
Если в пласте эллиптической формы работает п равноотстоящих друг от друга скважин (рис. 14), то дебит одной скважины определяется по формуле, предложенной В. Т. Мироненко [11] (IV.11)
где β находится из уравнения
(IV. 12).
х — координата центра скважины; L — малая полуось эллипса; σ — половина расстояния между скважинами.
§ 3. Метод эквивалентных фильтрационных сопротивлений
Одним из методов расчета дебитов многорядных батарей или цепочек скважин является метод эквивалентных фильтрационных сопротивлений Ю. П. Борисова.
Суммарный дебит цепочки из п скважин равен
(IV.13)
Используя электрогидродинамическую аналогию и учитывая, что аналогом объемного расхода является сила тока, а аналогом разности давлений — разность электрических потенциалов, выражение, стоящее в знаменателе, можно назвать фильтрационным сопротивлением. Оно складывается из внешнего фильтрационного сопротивления
(IV.14)
которое представляет собой сопротивление потоку от контура питания до галереи длиной В = 2σп, расположенной на расстоянии L от контура питания, и из внутреннего фильтрационного сопротивления
(IV. 15)
которое выражает собой сопротивление, возникающее при подходе жидкости к скважинам в зоне радиуса σ/π, где фильтрация практически плоскорадиальная.
Формула (IV. 13) примет вид
(IV.16)
Электрическая схема, соответствующая последней формуле, представляет собой два последовательно соединенных проводника с сопротивлениями ρ и ρ', с разностью потенциалов pк и pc и силой тока Q' (рис. 15).
Если в пласте имеется три цепочки с