Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Телков, А. П. Подземная гидрогазодинамика

.pdf
Скачиваний:
34
Добавлен:
20.10.2023
Размер:
7.79 Mб
Скачать

воряла последнему граничному условию в IX (2). В конечном сче­ те получается решение для распределения потенциала в пласте, свободное от указанных ограничений:

Ф(о . р', 1, h) = Ф0

- 4

f (р,р',Е,А);

 

IX (8 )

 

 

 

П0

 

 

 

 

 

 

 

 

 

у

z

i b

 

1Х(9)

а П0

;/ -

л г

 

Ао

й =

тh0-

где F — некоторая функция,

выражаемая

в рядах

и интегралах

через функции

Бесселя

[16].

 

 

имеем Ф ==

Фс. Тогда

При г — гс

(второе

граничное условие)

из IX (8 ) следует формула

IX (5),

где

 

 

 

 

Е

=

■F(o, h)

 

 

IX (1 0 )

3.Точное решение задачи о потенциале точечного стока

(источника) в однородно-анизотропном круговом

пласте. Приток

к несовершенной

скважине

 

 

Как известно, в области, содержащей источники или стоки,

потенциал удовлетворяет уравнению Пуассона [16]

 

Д<р = ф (*, у,

г):

 

IХ (1 1 )

Здесь А— оператор Лапласа;

<р— потенциал;

<1> (х, у,

г)—

плотность стоков как функция координат.

 

г =

Рассмотрим приток к точечному стоку с координатами

=0 , 2 = т), расположенному в круговом осесимметричном одно­

родно-анизотропном

пласте

(рис. 43).

 

 

 

 

В этом случае уравнение IX(1) будет иметь вид:

 

 

<Э3ф ,

1

д ср

1 д-ю

Ь (г) „ /

.

2

kr

 

+ -L £» +

й

-

 

*’ =

f

,Х(12>

дг2

г

дг

 

Р и с. 43. Приток к точечному стоку и линии стоков

91

где

q —-мощность

точечного

стока;

б (х) — функция Дирака.

 

Будем считать

кровлю ( z=

0) и подошву (z =

ft0)

непроницае­

мыми, т. е.

 

 

 

 

 

 

 

 

 

—■— 0

при

z = 0

и

z — h.Q

I Х(13)

На

внешнем контуре

(г =

R0)

примем

 

 

 

 

 

<р =

сро

при

г — R0

 

1Х(14)

 

Математическое

решение задачи

IX (12)—IX

(14)

выполнено

Сткляниным Ю. И. с помощью метода интегральных преобразова­ ний [16].

Точное решение для потенциала точечного стока имеет вид:

~ ch [т * (I_ А)1 ch т - ^ 7о 0 ?о. к)

ср (р, 6, Л)

= ср0 ------q^-

-

2

L Р_____________ J

Р_____________________

IX (15)

 

 

14 sh |Xj/p /? (н-i)

 

£ < ft

 

 

 

*Ro <=1

 

 

 

cp(p, £, ft)

=cp0

Д 7 1

 

chtjT

(1~ £)] ch7 ' /t/o(^

^

IX(16)

l > h

 

 

 

- R 0 i=i

 

l^i sh|Aj/p/? (|лг)

 

 

Здесь

[1/ — положительный

корень

уравнения

J 0 (p-i) =

0;

J0 (x) — функция

Бесселя

первого рода нулевого порядка;

J x (х) — функция

Бесселя

первого рода первого порядка;

безразмерные

координаты и

параметры:

 

 

 

 

Р

 

%h0

 

 

 

 

 

 

 

 

IХ( 17)

 

 

 

 

 

 

 

 

 

 

 

Найдем потенциал для линии стоков (рис. 43). В формулах IX (15)

и IX (16)

примем

 

 

 

 

 

 

 

 

q = q (ft) = const

 

при

0

<

ft <

1;

ft = -£-

IX (18)

 

 

 

 

 

 

 

 

 

 

 

«0

 

Очевидно,

что

потенциал для линии стоков будет

 

 

 

 

 

 

ф = ftol ? (Р,l,h)dh,

 

1Х(19)

 

 

 

 

 

 

 

о

 

 

 

 

 

где ср

определяется по формулам

IX (15) и IX (16). Интегрируя

IX (19),

после некоторых преобразований получим окончательно:

 

 

 

 

 

In R0 +

«о

c h ^ e s h !^ -(i_ /l) /0( ^

|)

Ф — Ф0

=

±

 

2 2

— i — s-J!-----5-----------

 

I < h

 

 

2nhho

 

 

i= l

 

И? sh (Aj/p /? (fif)

IX (20)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

92

Рассматривая приток к элементарному кольцу, интегрируя полученное уравнение вдоль кольца — стоков и усредняя потен­ циал вдоль стенки цилиндра — стоков, т. е. вдоль вскрытой части пласта, получим формулу для дебита:

^

2тс А0 (Ф0 — Фс)

_ 2к h0(Фб — Фс)

IX (2 2 )

 

 

 

 

где

 

 

 

IX (23)

 

h 1=1

^ 3 sh tL /?(N)

 

 

93

 

 

о

Ro

 

 

_ |

 

 

q

r

1

К*, г

^

1—<

1

 

 

*

Кg

1

=с:

С......

7

Р и с . 45. Схема притока к несовершенной скважине в двухслойном однородноанизотропном пласте

Функция IX (23) рассчитана на ЭВМ, графическое изображе' ние которой представлено на рис. 44.

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

Переходя от потенциала к давлению в формуле IX(22) и при­ нимая вместо давления р функцию Лейбензона Р, а вместо объем­ ного расхода Q— весовой расход G = уamQ. а затем делая обрат­ ный переход, получим

Q =

. р* ~ р*

IX (22')

ЦРат

Во

 

4. Потенциал несовершенной скважины в двухслойном однородно-анизотропном круговом пласте

Решение данной задачи (рис. 45) выполнено Ю. И. Сткляниным. Решения для распределения потенциала (Фх— потенциал для верхнего пласта, Ф2— потенциал для нижнего пласта) выра­ жаются в сложной интегральной форме. Для скважины радиуса гс потенциал на скважине определится интегралом

ь

ф с = I q (tj) <р (гс, г, т]) d г},

IX (24)

О

 

94

где т), q и ®— ордината, мощность

и потенциал точечного стока

соответственно.

можно добиться постоянства

Соответствующим подбором q

потенциала на забое вскрытой части

пласта, т. е. Фс = constant

при 0 < Z < Ь.

 

Удобнее за потенциал на стенке скважины принять среднее значение его, вычисленное по формуле:

IX (25)

о

где Oj (g, р, К) определяется сложным выражением [16], которое из-за громоздкости здесь не приводится.

5. Потенциал точечного стока (источника), горизонтальной дрены и несовершенной галереи в полосообразном однородно­

анизотропном пласте

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

 

Рассмотрим

приток к точечному

стоку

с координатами (tj,

/ х ) ,

расположенному несимметрично

в

полосообразном горизон­

тальном

однородно-анизотропном

пласте с

подошвенной водой

(рис. 46).

исходное

принимается уравнение

 

 

За

 

 

 

 

 

 

 

 

 

IX (26)

где

б (ztj)

и

б (х— 1г) — функции

Дирака.

 

Введем потенциал

 

 

 

 

 

 

 

 

 

 

 

 

IX (27)

и характеристику

анизотропности

пласта

 

 

 

 

 

 

 

 

 

IX (28)

 

Тогда

уравнение IX (26) примет

вид

 

 

 

 

 

 

 

 

 

 

IX (29)

Кровля и подошва считаются непроницаемыми, т. е.

— = 0

при

z = 0,

z = h0

IX (30)

дг

 

 

 

 

95

Р и с. 46. Схема притока к точечному стоку, горизонтальной дрене, несовер­ шенной щели в полосообразном пласте

На контуре питания для простоты принимается

9 = 0

при

х = 0,

х = /

IX(31)

Уравнение IX (29) с граничными условиями IX (30) и IX (31) можно решить методом интегральных преобразований, применяя последовательно косинус - и синус - преобразование Фурье с конечными пределами и формулы обращения. Такое решение дано в работе [16]:

«, ch m л 2L (й0 — г) ch /я то) —

?! =

Щ-Ъ ----------------------—т----------- —

IX (32)

г>т]

11 m~i т sh т ъ Л„ ~ I sin т я fi sin т к

 

 

0 / '

I

I

 

ооch т п ~ (/г0— T | ) c h m i z i

92 = - f f X

-----------Ц ------------

p - i

-----р

IX (33)

г<7:

m=1 т sh т тс n0 — /sin т и is i n

т я —

 

Формулы IX (32) и IX (33) дают распределение потенциала, вызванного точечными стоком и источником, в элементе анизо­ тропного полосообразного пласта. Их можно использовать при экспериментировании на щелевом лотке и для горизонтальной дрены. Можно получить потенциал и для несовершенной щели в том же пласте. При этом щель будем рассматривать как линию стоков, расположенных вдоль прямой х = от г = 0 до г = Ь, при постоянной мощности q (рис. 46). Тогда потенциал линии сто­ ков определится интегралом

Ф = [ 9 (х >г,п) d у}

IX (34)

о

 

96

Подставляя IX (32) и IX (33) в IX (34), интегрируя последнее уравнение и делая некоторые преобразования, получим

 

 

Ф ,=

I ср^х,

z , rt) d

7] =

 

 

 

 

zyb

о

 

 

 

 

 

 

 

 

* тих ,,

 

. ,

т их ,

 

 

= _

ch -у— (h0— г) sh —j~

b

 

 

_ , 0 , m г:у. , ,

mu e*

. m и x

 

 

 

1m2 s h -----

A0/s

n -----1 s i n ---------

 

 

 

 

 

/

0

 

;

l

 

 

 

, ffln

 

, min

,,

,. .

m - l. .

mu. x

 

 

ch------ z sh

-----

(h,— b) ■sin----- i- s i n -------

Ф» =

-^ -2

 

l

 

 

______ l

/ _

 

rrfi sh

—/г„

 

 

z < i

n

m = l

 

 

 

IX (35)

+

( —

2l? " 7 F “ x

ПР Н

 

0 <

X <

l 1

+

I

(l — x)

при

/, <

x <

IX(36)

2qj-

/

Е с л и b = h0,

то

из формулы

IX (36)

получается выражение

для потенциала

совершенней

щели в полосообразней залежи:

(

 

 

ПРИ 0 < х <

Ф ,=

 

 

 

 

 

 

1X (37)

V— 2q *2 (/ — х)

 

при < Х<(/

Если за горизонтальною дрену принять горизонтальную сква­ жину диаметром d, произвольно расположенную в пласте, то потенциал такой дрены определится интегралом

ь

Ф = ) С£г (х, г, т,) d У\

1Х(38),

е—d

 

Ряды IX (32) и IX (33) сходятся медленно, однако после улуч­ шения сходимости [16] для них можно получить приближенные выражения с достаточной степенью точности, более удобные для вычислений.

6. Анализ распределения потенциала в однородно-анизотропном: круговом пласте

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

Впервые детальный анализ распределения потенциала на поверхности забоя несовершенной по степени вскрытия скважины

4 Заказ 612

97

 

в однородном пласте с осевой симметрией, а также в пласте по всей мощности относительно несовершенной скважины, дан М. Мас-

кетом [3].

Представляет теоретический и практический интерес (напри­ мер, при решении некоторых статических и динамических задач теории конусообразования) показать более наглядно характер распределения потенциала в однородно-анизотропном ограничен­ ном пласте, вскрытом несовершенной скважиной (рис. 43). Исполь­ зуем решение IX (21), которое в безразмерных параметрах прини­ мает следующий вид:

 

 

 

=

^

( 1 ,

Ро,

Л ,

 

 

 

IX (39)

 

 

 

 

ch 1--- - (ij sh J L iHh(\4 -FA

 

F ( £, Ро, h,

р -

Е

 

Po__________ P

'

R p J

IX (40)

 

 

A sh X i- i U ih)

 

 

 

 

i—1

 

 

 

Расчет функции IX (40) производился

на ЭВМ для параметра

^ = 10 и относительного вскрытия

h =

0,3

при

фиксированных

значениях

r/R0■ Результаты

представлены

на рис. 47. Степень

отклонения

от

плоскорадиального

потока

для

каждого значе­

ния r/Rо можно

показать из

соотношения

 

 

 

 

 

5% =

- ^ ° ’5

/

Ь1

• 10096

 

I X(41)

 

 

 

Г£=0,°

 

 

 

 

 

 

Р и с . 47. Распределение потенциала скоростей фильтрации в круговом пласте (к определению зоны пространственного движения)

9 8

Такое отклонение составляет: 17,2; 11,2; 1,18 и 0,25% для отношения радиусов 0,01; 0,05; 0,10 и 0,20 соответственно.

Таким образом,

для зоны r0 ~ rJR 0 > 0,10 практически

можно считать поток

плоскорадиальным, поскольку на расстоянии

г0 0,10 отклонение составляет величину порядка одного процента.

На рис. 47 показаны линии равных потенциалсв, где линия г0 — = 0,10 представляет ссбой почти прямую, параллельную оси абс­ цисс.

Чтобы нагляднее представить себе размеры зон пространст­

венного и

плоскорадиального течения,

рассмотрим

один

пример.

Пусть

R 0 =

200 м, h0 =

20 м и х ^

1

(пласт сднорсднс-изотроп-

ный),

что соответствует

параметру с0 =

10. Для

 

вскрытия

h =

; ^ 0,3,

принимая га

0,10 как

границу

между

указанными

зо­

нами,

находим т--

0,10

R 0

0,10 •

200

20

м.

Как

видим,

зона пространственного притока оказалась равной одной мощности пласта. С увеличением относительного вскрытия, как это явствует из физических соображений, эта зона будет уменьшаться, прибли­ жаясь к скважине, и в предельном случае вскрытия пласта (сква­ жина совершенная) во всей области дренажа будет иметь место стро­ го плоскорадиальный приток. Что касается анизотропных пластов, то, очевидно, зона пространственного движения будет меньше по

сравнению с однородно-изотропными пластами и

в общем случае

г0 = / (р, Щ

 

Проведенный анализ подтверждает выводы

М. Маскета и

И. А. Чарного о размерах зоны пространственного движения при притоке к несовершенной скважине по степени вскрытия пласта. Маскет показывает, что в однородно-изотропном пласте ПрИ от'

носительном вскрытии

h = 0,5

плоскорадиальность

потока на­

рушается лишь при г <

2h0. И.

А. Чарный предлагает принимать

зону пространственного потока

радиуса г = (1 -ь 1,5)

h0.

Произведенная количественная оценка позволяет принять за критерий, характеризующий приток к несовершенной скважине, параметр р. При 10 и ft > 0,3 (что выполняется в большинст­ ве практических случаев) зона пространственного притока с высо­

кой

степенью точности может быть принята равной r0 = h0. При:

р >

Ю этот радиус будет несколько больше и для практических:

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

ный на схеме условного разделения притока на две зоны.

4*

99

7. Приток неньютоновсксй жидкости к несоверненной скважине

Нефти некоторых месторождений СССР (Башкирии, Азер' байджана, Узбекистана) содержат большое количество пара" фино-смолоасфальтеновых компонентов. Наличие таких компонен" тов может обусловливать проявление структурно-механических свойств нефти (динамическое напряжение сдвига и структурная вязкость) при ее фильтрации в пористой среде. Исследования по­ казали, что чем больше нефть содержит смолоасфальтеновых ком­ понентов и чем ниже температура, тем большее значение приобре­ тают ее структурно-механические свойства.

Известно, что для движения обычных (ньютоновских) жидко­ стей в пористой среде, когда абсолютная вязкость жидкости ос­ тается постоянной, справедлив закон Дарси. Движение вязко­ пластичных (неньютоноЕСких) жидкостей характеризуется дина­ мическим напряжением сдвига и структурной вязкостью. В этом случае фильтрация жидкости описывается так называемым обобщен­ ным законом Дарси или законом фильтрации с преденльым гра­ диентом давления сдвига, достаточно полно исследован м в ра­ ботах А. X. Мирзаджанзаде и др., Султанова Б. И., В. М. Ентова. Обобщенный закон Дарси можно записать по аналогии, вводя коэффициент структурной вязкости Дх'0) как функцию динамиче­ ского напряжения сдвига (т0) [17]

и

grad

— р + р grad

k

I *

Л Со)

IX (42)

 

 

Ti Со)

 

!

При некоторых условиях,

когда величина второго слагаемого

в уравнении

IX (42)

составляет незначительную долю* от всей

суммы и ей можно пренебречь без существенного ущерба для

практических расчетов, обобщенное уравнение

Дарси IX (42)

может быть

записано приближенно

в виде

 

 

и х — grad —— р =

— grad Ф

IX (43)

 

 

лСо)

 

 

Отсюда следует, что

если известно решение для

распределения

потенциала

скорости

фильтрации

Ф ньютоновской жидкости

при притоке к совершенной или несовершенной скважине, то оно может быть использовано для приближенных расчетов и для слу­

чая фильтрации вязко-пластичной

жидкости

путем введения

т) (т0) вместо коэффициента абсолютной

вязкости

р. ньютоновской

* Оценка величины grad

—- — .может

быть

произведена по формуле с [17]

 

Wo)

 

 

 

grad

Г[1+ к/и (г)]

 

х = А т

 

т,(хо)

 

ц

 

100

Соседние файлы в папке книги из ГПНТБ