Расчет контура растеканиясточных вод
Для осуществления гидродинамических прогнозных расчетов необходимо прежде всего определиться с граничными условиями поглощающего горизонта. В практике подземного захоронения вод поглощающий горизонт принимается однородным неограниченным (бесконечным) в плане, поскольку для целей ПЗС выбираются выдержанные по литологии, мощности и простиранию пласты-коллекторы.
В гидродинамическом отношении поглощающие горизонты всегда залегают в зоне затрудненного водообмена, для которой характерны низкие скорости движения подземных вод, которыми при расчете можно пренебречь. Исключение составляют поглощающие горизонты, гидродинамически взаимосвязанные с активно обводняющимися разрабатываемыми залежами углеводородов. В таких поглощающих горизонтах скорость движения пластовых вод должна учитываться.
Таким образом, для расчета радиусов растекания захороняемых вод наобъектах газовой промышленности применяются два варианта, в зависимости от гидродинамических условий поглощающего горизонта: в поглощающем горизонт е отсутствует движение пластовых вод; движение пластовых вод в поглощающем горизонте имеется. Первый вариант характерен для поглощающих горизонтов, несвязанных гидродинамически с разрабатываемыми углеводородными залежами. Это преимущественно поглощающие горизонты над продуктивных отложений, в которые захороняются СВ на глубокозалегающих месторождениях (Астраханском, Карачаганакском, Тенгизском и др.). Второй вариант относится к поглощающим горизонтам, подстилающим разрабатываемые залежи и гидродинамически с ними взаимосвязанным. Такие условия имеют место при захоронении вод под газовые залежи на месторождениях Тюменского Севера, Оренбургском НГКМ, Вуктыльском ГКМи др.
Вариант 1. Движение сточных вод воднородном неограниченном пласте без учета скорости движения пластовых вод(закачка в условиях бассейна) представляет собой равномерное растекание по поглощающему горизонту во все стороны от нагнетательной скважины. Фронт распространения закачиваемых вод имеет в плане форму окружности. Расчет радиуса этой окружности, т.е. контура продвижения сточных вод в поглощающем горизонте, производится по формуле
, (14)
r= (15)
где г - радиус контура распространения закачиваемых сточных вод, м;
Т - время движения сточных вод от нагнетательных скважин, сут;
Q -объем закачки сточных вод, м3/сут;
m -эффективная мощность поглощающего горизонта, м;
n -коэффициент пористости коллекторов поглощающего горизонта, доли единицы.
Для создания «запаса прочности» прогнозных расчетов в знаменатель формулы (15)вводится коэффициент вытеснения Кв пластовой воды, численное значение которого из практики соответствует около 0,7. Тогда окончательная расчетная формула радиуса распространения сточных вод в поглощающем горизонте имеет вид
R= , (16)
В целях повышения надежности гидродинамических прогнозов рекомендуется принимать жесткие условия расчетных показателей посредством понижающих коэффициентов к значениям пористости и эффективной мощности пород. Значения этих коэффициентов могут колебаться от 0,9до 0,7 для средней пористости и от 0,9 до 0,5 для эффективной мощности.
Вариант 2. В поглощающем горизонте, гидродинамически взаимосвязанном с разрабатываемой газовой залежью, образовалась депрессионная воронка. Движение пластовых вод в пределах депрессионной воронки направлено к ее центру, где отмечается минимальное пластовое давление. В этом же направлении будут двигаться захороняемые сточные воды. Скорость движения пластовых вод в поглощающем горизонте определяются поформуле Дарси:
V=Kф J, (17)
где V - скорость фильтрации подземных вод, м/сут;
Кф – коэффициент фильтрации пород, м/сут;
J -гидравлический уклон потока.
При закачке СВ в поглощающий горизонт в случае наличия движения подземных вод (т.е. в условиях потока)образуется структура фильтрационного течения, резко отличающаяся от структуры движения жидкости, закачиваемой в условиях бассейна (когда V = 0). В условиях потока плоскорадиальное движение нагнетаемой жидкости, которое имеет место в бассейне, вследствие наложения на него плоскопараллельного движения деформируется и становится плоско криволинейным (рис. 6).
По В.А.Боревской, И.Т.Гаврилову, В.М. Гольдбергу и др. (1978), поле фильтрационного течения делится на две области: внутреннюю и внешнюю. Во внутренней области, примыкающей к скважине, все линии тока заканчиваются на скважине. Во внешней области линии тока обтекают внутреннюю область и минуют скважину. С удалением от внутренней области (вверх и вниз по оси Y) искривления линий тока уменьшаются, и они приближаются к прямым линиям естественного плоскопараллельного течения подземных вод. Внутренняя область отделена от внешней нейтральной линией тока MAN. Закачиваемая жидкость распространяется только во внутренней области, называемой областью растекания. Таким образом, при наличии естественного движения подземных вод сточные воды распространяются не равномерно во все стороны, как в условиях бассейна, а в пределах области растекания. Следовательно, для того чтобы оценить масштабы распространения сточных вод в пласте, следует определить размеры области растекания. Вверх по потоку границей этой области является лежащая на оси Х такназываемая водораздельная точка А, в которой противоположно направленная скорость естественного течения и скорость, обусловленная действием скважины, взаимно уравновешиваются, а скорость фильтрации в этой точке равна нулю. Вниз по потоку область питания распространяется неограниченно.
Размеры области растекания стоков определяются по формулам
XA = (18)
Y0= ; (19)
Y¥ = ; (20)
где ХA - расстояние от скважины до водораздельной точки, м;
Y0- ширина области растекания по линии скважины (X = 0), м;
Y¥- ширина области растекания в удалении от скважины вниз по потоку (Х=¥), м;
m -эффективная мощность поглощающего горизонта, м;
V – скорость фильтрации подземных вод, м/сут.
Расчет расстояний продвижения закачиваемых вод производится методом подбора, исходя из формул:
для главной линии тока
T0= (21)
Линии тока в неограниченном пласте в условиях бассейна
Линии тока в неограниченном потоке:
Г0- область квазирадиального течения (Rг0<0,2XА)
Перемещение фронта закачиваемой жидкости в неограниченном потоке
Рис. 6. Динамика движения жидкости в условиях бассейна и потока в однородном неограниченном пласте
для других линий тока
T= ,(22)
X= ,(23)
Y= ,(24)
Т0 – время перемещения закачиваемых вод от скважины до какой-либо фиксированной точки по главной линии тока, сут;
Т - то же, по другим линиям тока, сут;
n -коэффициент пористости поглощающего горизонта, доли ед.;
ХA- расстояние от скважины до водораздельной точки, м;
V - скорость фильтрации подземных вод потока, м/сут,
X, Y - координаты точки, до которой рассчитывается время перемещения закачиваемых вод, м.
Продвижение СВ определяется по главной линии тока и по ширине потока на линии скважин, с помощью чего графически устанавливается весь контур распространения вод. По главной линии тока прогноз перемещения стоков производится вверх и вниз по формуле (21). Продвижение СВ по ширине потока на линии скважин определяется исходя из формулы (22), которая после соответствующих преобразований имеет вид
, (25)
где Т - время перемещения сточных вод по ширине потока на линии скважины, м. Остальные обозначения те же, что и в предыдущих формулах.