6.3. Моделирование добычи газа из гидратов методами понижения давления, нагрева гидратосодержащих пород и комбинированным методом

Рассмотрим перспективность добычи газа из гидратов на примере задачи притока газа к вертикальной скважине, полностью вскрывшей гидратонасыщенный пласт. Тогда система уравнений, описывающих разложение гидрата в пористой среде, будет иметь вид:

а) закон сохранения массы для газа и воды:

(6.6)

где P — давление, T — температура, S — водонасыщенность, v — гидратонасыщенность, z — коэффициент сверхсжимаемости; r — радиальная координата; t — время; m — пористость, ρ_g, ρ_w, ρ_h — плотности газа, воды и гидрата соотвественно; k(v) — проницаемость пористой среды в присутствии гидратов; f_g(S), f_w(S) — функции относительных фазовых проницаемостей для газа и воды; μ_g, μ_w — вязкости газа и воды; — массовое содержание газа в гидрате;

б) уравнение сохранения энергии:

(6.7)

где: ρС_e — теплоемкость породы и вмещающих флюидов; c_g, c_w — теплоемкость газа и воды соответственно; H — теплота фазового перехода гидрата; η — дифференциальный адиабатический коэффициент; ψ — коэффициент дросселирования (коэффициент Джоуля-Томсона); λ_e — коэффициент теплопроводности породы и вмещающих флюидов.

В каждой точке пласта должно выполняться условие термодинамического равновесия:

Т = A ln P + B   (6.8) где A и B — эмпирические коэффициенты. Зависимость проницаемости породы от насыщенности гидратов принято представлять в виде степенной зависимости:

k (v) = k₀ (1 - v)N    (6.9)

где: k₀ — абсолютная проницаемость пористой среды при отсутствии гидратов; N — константа, характеризующая степень ухудшения проницаемости с ростом гидратонасыщенности.

В начальный момент времени однородный и единичной мощности пласт имеет давление Р₀, температуру Т₀ и насыщенность гидратами v₀. Метод понижения давления моделировался заданием на скважине постоянного дебита, а тепловой метод – тепловым источником постоянной мощности. Соответственно при комбинированном методе задавались постоянный расход газа и мощность теплового источника, необходимая для устойчивого разложения гидратов.

При моделировании добычи газа из гидратов рассматриваемыми методами учитывались следующие ограничения. При начальной пластовой температуре 10°С и давлении 5,74 МПа коэффициент Джоуля-Томсона составляет 3-4 градуса на 1 МПа депрессии. Таким образом, при депрессии 3-4 МПа забойная температура может достичь температуры замерзания воды. Как известно, замерзание воды в породе не только снижает проницаемость призабойной зоны, но и приводит к более катастрофическим последствиям — смятию обсадных колон, разрушению коллектора и т.д. Поэтому для метода понижения давления принималось, что за 100 суток работы скважины забойная температура не должна снизиться ниже 0°С. Для теплового метода ограничением является рост температуры на стенке скважины и самого нагревателя. Поэтому при расчетах принималось, что за 100 суток работы скважины забойная температура не должна превысить 110°С. При моделировании комбинированного метода учитывались оба ограничения.

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

Q* = Q — E/q    (6.10)

где Q — дебит газа на забое, м³/сут.; E — подводимая к забою тепловая энергия, Дж/сут.; q — теплота сгорания метана (33,28∙10⁶), Дж/м³.

Расчеты проводились при следующих параметрах: P₀ = 5,74 МПа; T₀ = 283 К; S = 0,20; m = 0,35; ρ_h = 910кг/м³, ρ_w = 1000 кг/м³; k₀ = 0,1 мкм²; N = 1 (коэффициент в формуле (4)); μ_g = 0,014 мПа∙с; μ_w = 1 мПа∙с; ε = 0,134; A = 7,28 К; B = 169,7 К; ρС_e = 1,48∙10⁶ Дж/(м³∙К); c_g = 2600 Дж/(кг∙К), c_w = 4200 Дж/(кг∙К); H = 0,5 МДж/кг; λ_e = 1,71 Вт/(м∙К). Результаты расчетов сведены в табл. 6.1.

Таблица 6.1

Дебиты газа и воды при добыче газа из гидратов различными методами

Насыщенность гидратами, %	20	40	60	80
	Метод понижения давления
Дебит газа, м3/сут	22800	15500	6700	0
Дебит воды, дм3/сут	0	0	0	0
	Тепловой метод
Дебит газа, м3/сут	- 0,22	6,9	7,8	7,95
Дебит воды, дм3/сут	0	0,2	1,38	5,2
	Комбинированный метод
Дебит газа, м3/сут	22950	15700	6980	9,45
Дебит воды, дм3/сут	0	2,1	6,5	14,0

Анализ этих результатов расчетов показывает, что метод понижения давления является пригодным для гидратных пластов, где насыщенность гидратами невелика, а газ или вода не потеряли свою подвижность. Естественно, что при увеличении гидратонасыщенности (а значит, сокращении проницаемости согласно уравнению (6.9)) эффективность этого метода резко падает. Так, при насыщенности пор гидратами более 80% получить приток из гидратов за счет снижения забойного давления практически невозможно. Таким образом, разработка гидратных залежей за счет понижения давления возможна только при закачке ингибиторов в призабойную зону, что значительно увеличит себестоимость добываемого газа.

Тепловой метод разработки газогидратных месторождений пригоден для пластов, имеющих высокое содержание гидратов в порах. Однако, как показывают результаты расчетов, тепловое воздействие через забой скважины малоэффективно. Это связано с тем, что процесс разложения гидратов сопровождается поглощением тепла с высокой удельной энтальпией 0,5 МДж/кг (для примера: теплота плавления льда составляет 0,34 МДж/кг). По мере удаления фронта разложения от забоя скважины все больше энергии тратится на прогрев вмещающих пород и кровли пласта, поэтому зона теплового воздействия на гидраты через забой скважины исчисляется первыми метрами.

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