К оглавлению журнала

 

УДК 622.276

© Ю.П. Коротаев, 1992

НОВАЯ МЕТОДИКА ОПРЕДЕЛЕНИЯ ФИЛЬТРАЦИОННЫХ ПАРАМЕТРОВ ПРИ СТАЦИОНАРНЫХ РЕЖИМАХ ФИЛЬТРАЦИИ

Ю. П. КОРОТАЕВ (ГАНГ)

Сегодня на практике при решении таких задач по фильтрации газа и жидкости, как обработка результатов гидродинамических исследований кернов, скважин и проектирование разработки месторождений природного газа применяется двучленный закон, в котором зависимость между градиентом давления dp/dx и скоростью фильтрации v представляется в виде

где m и r соответственно динамическая вязкость и плотность флюида; k и l – соответственно проницаемость и макрошероховатость пористой среды.

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

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

С целью проверки указанного противоречия были поставлены специальные прецизионные исследования многочисленных естественных и искусственных кернов с применением акустико-гидродинамического метода (АГДМ) при одновременном измерении гидродинамических и акустических характеристик. Исследования АГДМ проводились в широком диапазоне изменения скоростей (дебитов), снималось до 20–40 точек.

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

Анализ исследований показал, что каждой пористой среде соответствует свое критическое значение скорости фильтрации Vкp, которое соответствует верхней границе применимости закона Дарси. В диапазоне изменения V<=Vкр, т. е. при малых скоростях фильтрации справедлив закон Дарси:

При этом практически кроме акустического фона при фильтрации отсутствует аэродинамический шум.

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

В диапазоне изменения скоростей фильтрации VxVкр зависимость между dp/dx и v, как показали результаты проведенных исследований, имеет вид

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

P1–P2=Aq–Bqкрq+Bq2, (4)

где A=mL/kF; B=rL/lF2; L и F – соответственно длина и площадь поперечного сечения образца пористой среды; P1 и P2 соответственно давление на входе и выходе из образца пористой среды; q – дебит жидкости; qкркритический дебит жидкости, соответствующий Vкр. При режимах, когда q<=qкр:

P1–P2=Aq. (5)

Решение (3) для прямолинейной одномерной фильтрации газа дает аналогичную формулу

P12–P22=aQ–bQкр+bQ2, (6)

где a=2mLpат/kF; b=2LPатrат/lF2; Q – текущий дебит газа при атмосферном давлении Pат; Qкркритический дебит газа при Рат.

Критический дебит соответствует максимальному дебиту, когда еще справедлив закон Дарси. При режимах, когда Q<=Qкp:

P12–P22=aQ. (7)

Таким образом, согласно уравнениям (4)–(7) изменяется методика проведения и обработки результатов исследований кернов, состоящая в нахождении не только коэффициентов а и b, но и критической скорости (критического дебита Qкр ). Методика же проведения исследований кернов состоит в том, что исследования проводятся во всем диапазоне изменения режимов от закона Дарси (7) до трехчленного закона (6).

Рассмотрим пример исследования и обработки кернов по предлагаемой методике. Разделив обе части уравнений (6) и (7) на Q получим для режимов, когда Q>Qкр:

где аф–bQкр, отрезок отсекаемой на оси ординат при Q=0. Для режимов, когда Q>=Qкр получим

Обработка результатов исследований кернов приводит к удельной индикаторной кривой, состоящей из двух прямых (рис. 1). Начальный горизонтальный участок (7) соответствует закону Дарси, и второй участок (2) соответствует трехчленному закону (8).

Для повышения точности целесообразно проводить исследования кернов АГДМ с установкой микрофона на соответствующем расстоянии на выходе газа из керна.

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

При плоскорадиальной фильтрации зависимость между dp/dr и v с учетом сказанного будет

Когда V<=Vкр для всего интервала от радиуса скважины Rс до радиуса контура питания Rк получим формулу при фильтрации газа согласно закону Дарси

где Pк и Pз соответственно давление на контуре питания и забое скважины.

Когда q<=qкр при фильтрации жидкости на режимах:

Pк–Pз=Aq, (13)

где A=(m/2pkh) Ln(Rk/Rc) (14)

Когда V>Vkp, то в интервале пласта от Rc до Ro справедливо уравнение (10) и от Ro до Rk закон Дарси. Радиус укрупненной скважины Ro соответствует критической скорости фильтрации Vkp.

Тогда решение уравнения (10) приводит при фильтрации газа к формуле

Pк2-Pз2= (a-bQкр)Q + bQ (Q - Qкр Ln(Q/Qкр)) (15)

где а соответствует (12), b=rатРат/(2p2h2lRc)

Поэтому при фильтрации жидкости имеем

Pк-Рз=(A-Bqкр)q+Bq(q-qкрLn(q/qкр)), (16)

где А соответствует (14),В=rатРат/(2p2h2lRc) Обработку результатов исследования скважин с целью определения коэффициентов фильтрационного сопротивления а, b и Qкр, проводим следующим образом. Поделив левую и правую части уравнения (11) на Q, получим для газовых скважин

DP2/Q=а. (17)

Обрабатывая согласно (17) результаты исследований скважин в координатах DP2/Q и Q в пределах изменения дебитов при Q<=Qкр, получим горизонтальную прямую, параллельную оси Q, которая от секает на оси ординат отрезок, равный коэффициенту а (прямая 1 на рис.2). При дебитах Q>Qкр из уравнения (15) получим

где

Обрабатывая результаты исследовании в координатах DP2/Q и Q получим прямую 3 (см. рис. 2), тангенс угла наклона которой равен b и которая на оси ординат отсекает отрезок, равный аbQкр.

Таким образом, применяя для обработки двучленную формулу, допускается ошибка в определении коэффициента а на величину bQкр, т.е. часто значительно завышается величина проницаемости по сравнению с фактическим ее значением и занижается b. Величину критического дебита определяем по точке пересечения прямых (см. рис. 2). Аналогичным образом проводится обработка исследований нефтяных скважин, исходя из (13) и (16).

Из практики обработки результатов исследований газовых скважин известно, что более 50 % скважин не удается обработать по двучленной формуле из-за того, что в координатах DP2/Q и Q вместо прямой получается кривая, иногда сходится с гиперболой. Для обработки таких аномального вида индикаторных кривых автором в свое время была предложена формула

DP2=aQ+bQ2+c, (19)

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

Сегодня на практике формула (19) используется повсеместно без соответствующей проверки наличия или отсутствия жидкости на забое, когда она теряет физический смысл. Предложенная методика обработки результатов исследований согласно формулам (12) и (18) позволяет интерпретировать индикаторные кривые наличием двух законов фильтрации, считавшиеся ранее аномального вида и без инструментального обоснования обрабатываемые по формуле (19).

Для скважин, гидродинамически несовершенных по степени вскрытия, а также несовершенных по характеру вскрытия, получены формулы, подобные (15) и (16) с соответственно иной структурой коэффициентов а и b. Обработка же результатов подобна совершенным скважинам.

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

Нефтяные и водяные скважины работают и исследуются обычно в диапазоне при Q<=Qкр.

Методика была подтверждена сотрудниками и студентами кафедры разработки и эксплуатации газовых и газоконденсатных месторождений МИНГ им. И.М. Губкина по многим скважинам Уренгойского, Карачаганакского, Астраханского, Шебелинского и других месторождений. Например, величина Qкр для скважин Еняхинской площади Уренгойского месторождения находится в пределах 300– 700 тыс. м3/сут.

Критический дебит соответствует максимальному энергосберегающему дебиту надежной эксплуатации скважин, так как при Qкр потери давления DP2 прямопропорциональны Q, а при Q>Qкр потери давления растут более интенсивно за счет влияния члена с QQ, что приводит к интенсивному росту градиента на стенке забоя скважин, сопровождается возникновением акустических колебаний, способствующих разрушению пласта и выносу песка из скважины.

На основе Qкp предлагается новый технологический режим эксплуатации скважин, состоящий в поддержании постоянного энергосберегающего дебита Qкp=const в течение всего или основного периода разработки месторождения. Разработка месторождения с применением Qкp=const обеспечивает эксплуатацию скважин с минимальными потерями пластовой энергии и без выноса песка из скважины.

Считаем целесообразным при проектировании и разработке газовых месторождений рассмотрение с технико-экономической оценкой варианта с режимом Qкр=const.

Несмотря на необходимость опережающего бурения значительного числа скважин при их эксплуатации на режиме Qкр=const могут быть получены благоприятные технико-экономические показатели за счет резкого снижения мощности ДКС, переноса срока их ввода практически на конец периода постоянной добычи газа, роста коэффициента газоотдачи особенно в неоднородных коллекторах и снижения затрат на капитальный ремонт скважин за счет уменьшения их осложнений и аварий.

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

ABSTRACT

Stationary investigations of gas filtration at filtration rates lower or higher than critical have been conducted experimentally on natural and artificial cores. A new well' survey procedure, as well as data interpretation methods are recommended.

РИС. 1. РЕЗУЛЬТАТЫ ИССЛЕДОВАНИЙ КЕРНА № 297:

1 – по формуле (9) и 2 – (8)

РИС. 2. ОБРАБОТКА РЕЗУЛЬТАТОВ ИССЛЕДОВАНИЯ СКВАЖИН ПРИ ПЛОСКОРАДИАЛЬНОМ ПРИТОКЕ ГАЗА:

1 – по формуле (17), 2 – (18) в координатах DP2/Q и Q, 3 – (18) в координатах DP2/Q и Q–Qкр Ln (Q/Qкр)

Сайт создан в системе uCoz