3_Подсчёт запасов газа. Лекция 3 Подсчёт запасов газа
 Скачать 104.5 Kb.
|
|
Лекция 3 Подсчёт запасов газа При подсчете запасов газа различают свободный газ, то есть из газовых залежей и газовых шапок нефтегазовых (газонефтяных) залежей, и газ, растворенный в нефти (попутный газ). Объемный метод подсчета запасов свободного газа основан на тех же принципах определения объема залежи, что и объемный метод подсчета запасов нефти: Подсчет запасов газа и конденсата производился объемным методом по общепринятой методике в соответствии со следующими формулами: Для газа: Qг = F* hэфф г* Kп* Кг* (Рн* н – Рк* к)* f*.  или Кр-барический к-тгде: F - площадь газоносности, тыс.м2 - определяется в пределах контура газоносности на основе подсчетного плана, совмещенного со структурной картой по кровле пласта. Hэфф г – средняя (средневзвешенная) газонасыщенная толщина коллекторов, м - определяется по данным интерпретации материалов геофизических исследований скважин как средневзвешенное значение эффективной газонасыщенной толщины определяется как частное от деления объёма газонасыщенных пород, рассчитанного по карте эффективных газонасыщенных толщин, на соответствующую площадь. Кп - коэффициент открытой пористости, доли ед - определяется в результате комплексной интерпретации материалов ГИС. Рассчитывается средневзвешенная величина коэффициента открытой пористости по эффективной газонасыщенной толщине для коллекторов пластов. Кг - коэффициент газонасыщенности, доли. ед определяется лабораторным путём с учетом содержания связанной воды и соотносится с данными ГИС; при наличии нефтяных оторочек с учётом остаточной нефтенасыщенности Кг=1 – Кво – Кно; состав газа:СН4 - метан, С2Н6 - этан, С3Н8 – пропан,С4Н10 – бутан. При построении карт газонасыщенности в трёхмерных построениях нужно учесть, что значение коэффициента газонасыщенности на контуре не ноль, а 0,3-0,4. И кружочки в зоне предельного насыщения, а в переходной зоне насыщенность зависит ещё и от положения над ГВК. Рн и Рк - начальное и конечное пластовое давление в залежи, физ. ат; Пластовое давление, давление в недрах складывается из давления породы и давления насыщающего её флюида. Давление породы (геостатическое давление) создаётся весом породы и определяется плотностью пород и их мощностью. Средний градиент геостатического давления составляет 0,023 МПа на 1 м мощности (при плотности пород 2,3 г/см3). Давление насыщающего породу флюида (гидростатическое давление) создаётся весом воды, заключённой в системе. При средней плотности пластовой воды 1,23 г/см3 градиент давлений составляет 0,01-0,0133 МПа/м. Коэффициент аномальности давлений – отношение пластового давления (Рпл) к нормальному гидростатическому (Рн.г) Ка=РплРн.г, Коэффициент аномальности до 1,25 (для гиганских резервуаров), до 1,4 (крупные). Непско-Ботуобинская область имеет низкие давления. Теоретическая формула Рпл гидрост=ρ*g*H/100 (в физ. атм), где ρ – плотность воды в г/см3, g-скорость свободного падения 9,8 м/ч Н – глубина определения давления, м Среднее начальное пластовое давление рассчитывается на середине объёма в водоплавающей залежи, где за середину принимается отметка расположенная на 1/3 высоты залежи над поверхностью ГНК. Для газовых объёмов, давление приводится на отметку, соответствующую середине высоты залежи. Значение пластового давления получают по данным замеров глубинным манометром или измеренное давление на устье приведённое к глубине кровли пласта с учётом веса столба газа. 1293*10-9*Н*γ Рпл=Ру* е Где Ру – давление на устье закрытой скважины , МПа; е – основание натуральных логарифмов, е=2,718 Н – глубина кровли пласта в скважине, м Конечное пластовое давление рассчитывалось на момент, когда давление на устье скважины после извлечения запасов газа составит 1 атм. или 0.1 МПа. Запасы газа рассчитывают для условия полного истощения пластовой энергии, при котором Рк = 0.1МПа, коэффициент к при этом практически равен 1, что принято при подсчете запасов. (Перевод параметра пластового давления из атм в МПа умножением на 0,101322) Кр=Р00 - Ростост, Рст где Рост – среднее остаточное абсолютное давление на устье скважины до 1 атмосферы 1293*10-9*Н*γ Рост = е , где Н – глубина скважины в см, γ – плотность газа по воздуху н и к - поправки на сверхсжимаемость углеводородных газов для начальных и конечных пластовых условий (поправки на отклонение углеводородных газов от закона Бойля— Мариотта соответственно для давлений p0 и. pст). При подсчете геологических запасов пластового и сухого газа, состав и свойства газоконденсатных систем в данных термодинамических условиях учитываются коэффициентами н и к (поправки на сверхсжимаемость углеводородных газов для начальных и конечных пластовых условий). Коэффициенты н и к представляют обратные значения коэффициентов сжимаемости zн и zк, рассчитанных, соответственно, для начальной и конечной величины пластового давления. н=1/ zн, Z – зависит от пластового давления, температуры и состава газа. Значение коэффициентов сжимаемости zн устанавливается по лабораторным исследованиям состава газа или расчётным путём. Если количество неуглеводородных и тяжёлых углеводородных компонент в газе менее 10 %, коэффициент сверхсжимаемости определяется по двум параметрам: pкр – пластового давления и Ткр - температуры, по которым рассчитывают приведённые параметры pпр и Тпр. Если компонентный состав газа известен, то рассчитываются псевдокритические параметры давления и температуры, затем находят суммы псевдокритических давлений и суммы псевдокритических температур, по ним определяют приведённые давления и температуры, первое равно отношению пластового давления (атм) к сумме псевдокритических давлений: Рпр=Рпл/∑Рпс кр и приведённая температуря равна отношению температуры пластовой к сумме псевдокритической температуре: Тпр=Тпл/∑Тпс кр. О  пределяется компонентный состав газа, содержание всех элементов С1, С2, С3,С4, С5, С6+в, N2. CO2, относительная плотность газа, пластовое давление и температура. Вычисляются приведённое давления Рпр и приведённая температура Тпр. Расчётным путём или по графику Брауна (зависимость коэффициента сверхсжимаемости для природного газа от приведённого давления и температуры) определяется коэффициент сверхсжимаемости, через параметры (Рпр, Тпр). m, e, p, c, s – соответственно процентное содержание метана, этана, пропана, бутана и сероводорода H2S. В Западной Сибири Для неокома коэффициент сверхсжимаемости меньше 1. f - поправка на температуру, для приведения объема газа к поверхностным условиям; Поправка за температуру учитывает изменение объема газа при приведении к поверхностным условиям (стандартной температуре 200С) и определяется по формуле: f = Тст/ Тпл. где Тст - стандартная температура (200С) в К; Тпл - пластовая температура в К. или f =(273 +20)/(273+tпл) - коэффициент газоотдачи. Коэффициент извлечения газа условно принят равным единице. Газ в залежи – это пластовый газ, содержащий конденсат, газ без конденсата – это «сухой» газ. Запасы «сухого» газа рассчитываются как произведение пластового газа на величину мольной доли «сухиго» газа. Пример Qг =330000*14* 0,15*0,70*(290*1,15 – 1*1)*0,83= 133875472,5 тыс. м3) или 133875,5 млн. м3 F =330000 тыс. м2; Hэфф г=14 м; Кп =0,15 доли ед. Кг=0,70 доли ед Рн=290 физ. атм Рк=1 физ. атм н =1,15 к=1 f =(273 +20)/(273+80)=0,83 Балансовые запасы газа, растворенного в нефти, рассчитывают по формуле Qг.бал =Qн.бал r0, где Qг.бал, Qн.бал —балансовые запасы газа, м3 и нефти, т; r0—содержание газа в нефти при начальном пластовом давлении, м3/т. Величина извлекаемых запасов газа, растворенного в нефть, зависит от режима работы нефтегазоносных пластов. При водонапорном режиме (при котором разрабатывается подавляющее большинство месторождений) газовый фактор в процессе эксплуатации залежи мало изменяется во времени и извлекаемые запасы газа, растворенного в нефти, подсчитывают по упрощенной формуле Qг.изв =Qн.изв r, где r—газовый фактор, м3/т, замеренный на поверхности при давлении 0,1 МПа; Qн.изв —извлекаемые запасы нефти, т; Qг.изв —извлекаемые запасы газа, растворенного в нефти, м3. Конденсат(С5+в) Конденсат представляет собой газообразные УВ, состава С5+высш с растворёнными в них и находящимися в подчиненном состоянии жидких и твёрдых. УВ. Соединение твёрдых УВ в 1 л конденсата характеризуется конденсатным. фактором КФ=10-500 (чем больше КФ, тем больше конденсат похож на нефть). В пластовых условиях (Тпл=100-150, Рпл=20-30 МПа) конденсат находится в состоянии парогазовой системы. Поднимаясь на поверхность (Т=10-20, Р=0,1мПа) парогазовая система конденсируется и на поверхности становится жидкостью. В отличии от нефти цвет конденсата может менятся от темно-окрашенных тонов к бесцветному. Конденсат в отличии от нефтей не содержит асфальтенов. Плотность 0,77-0,80. Вязкость нефтей и конденсатов является характеристикой, отражающей способность молекул, движущихся флюидов, преодолевать сопротивления движению. Вязкость определяет скорость перемещения флюидов по пустотному пространству пород, она бывает динамическая и кинематическая. Кинематическая это отношение динамической вязкости к плотности воды, динамическая это произведение времени истечения флюида через капиляр определённого сечения и коэф-т пропорционально зависящий от диаметра вязкозиметра Для конденсата: Подсчет геологических запасов стабильного конденсата (С5+в) выполнен исходя из запасов свободного газа и среднего потенциального содержания С5+в. Qк =(Qг/106*П) *. где: Qг – запасы газа в залежи, м3; П - среднее потенциальное содержание С5+ в(стабильного конденсата) в залежи, г/м3 ; Потенциальное конденсатосодержание  Где К – содержание С5 и выше в сыром конденсате, L – в отсепарированном конденсате из расчета на 1 м3 пластового газа. - коэффициент конденсатоизвлечения. Пример Qк =((140000*106)/106*290)*0,5= 20300 тыс. т) Qг = 140000 млн м3; П = 290, г/м3; =0,5, доли ед. Подсчёт запасов этана, пропана и бутана Qкомп =(Qг/106*Пкомп) Qг – запасы газа в залежи, м3; Пкомп – содержание этана, пропана или бутана в газе, в г/м3 Пример Qэтана =(140000*106/106*75,5) = 10570 тыс. т) Qг = 140000 млн. м3; Пком = 75,5, г/м3; Метод по падению давления Основан на связи количества извлекаемого газа с величиной падения давления в процессе разработки газовой залежи. Считается, что для газовой залежи работающей на газовом режиме эта зависимость постоянна во времени, то есть количество газа, добываемого при снижении пластового давления на 0,1 МПа, постоянно в процессе разработки. (Qn(i+1) — Qгi)/(ii — i+1 i+1)=const, гдеQn(i+1) и Qгi – добытое количество газа за вторую и первую даты (измерения проводят раз в квартал), м3, i и i+1 – соответствующие первой и второй датам пластовые давления в залежи, МПа; i и i+1 – поправки на сжимаемость газа соответственно при давлениях i и i+1, обратные коэффициенту сжимаемости Z, то есть =1/Z Если на первую дату подсчета в начале разработки залежи добыто Q1объемов газа, при этом давление в залежи составило p1, а на вторую более позднюю дату отобрано Q2 объемов газа и давление равнялось р2, то добыча газа за этот период (от первого до второго подсчета) на единицу падения давления составит Q==(Q2-Q1)/(p1-p2) Исходя из того, что и в дальнейшем при падении пластового давления в залежи до некоторой его конечной величины будут добываться одинаковые количества газа на единицу падения давления, получают следующую формулу для подсчета запасов газа: Qг=(Q2 — Q1)(22 — 11) / (11 — 22), где Qг —промышленные запасы газа на дату, когда уже было отобрано газа Q2, м3. Начальные запасы свободного газа будут определяться Qг=(Qг(i+1) — Qгi)/(ii — i+1 i+1)/ 00 , где 0 – начальное пластовое давление в залежи, МПа; 0– соответствующая ему поправка на сжимаемость газа соответственно при давлениях 0 Если замеры отбора ведутся с начала разработки, то формула имеет вид: Qг0=Qгi 00 /(00— i i) Для залежей с водонапорным режимом метод по падению давления неприменим, так как при подсчете запасов газа этим методом предполагается, что первоначальный объем пор пласта, занятый газом, не меняется в процессе эксплуатации. Метод падения давления применим на залежах работающих на газовом режиме. Так как он позволяет определять запасы газа дренируемого объёма, то непременным условием отнесение подсчитанных запасов к начальным балансовым является вовлечение в разработку всего объёма залежи. Если в залежи начинается упруговодонапорный режим, сопровождающийся внедрением пластовой воды в залежь, в этом случаи при подсчёте запасов необходимо использовать данные того периода, когда залежь работала на газовом режиме, обычно это соответствует промежутку времени, за который отбирается 5-10 % начальных запасов залежи при равномерном дренировании всего объёма. Если объём залежи дренируется лишь частично, а в залежь начинает внедрятся вода, то применение метода падения давления может привести к существенным погрешностям. При газо-водонапорном режиме в формулу необходимо вводить поправку на количество газа, вытесненного за определенный период времени напором воды (Q'). Тогда формула для подсчета запасов примет следующий вид: Qг=(Q2 — Q1 — Q') 22 / (11 — 22). Остаточное давление в этом случае учитывать нет необходимости. Если количество газа, вытесненного напором воды, определить невозможно, запасы газа следует подсчитывать объемным методом. Если месторождение газоконденсатное, то после определения запасов газа подсчитывают запасы газоконденсата: Qк =(Qг/106*П) *, где Qг – запасы газа в залежи, м3; П - среднее потенциальное содержание С5+ в(стабильного конденсата) в залежи, г/м3 ; - коэффициент конденсатоизвлечения. Объемный метод подсчета запасов газа можно применять на любой стадии разведанности залежи. Для использования метода по падению давления необходимо иметь данные эксплуатации скважин. Объемный метод применяется при любом режиме работы пласта. Метод по падению давления эффективен лишь при газовом режиме, при водонапорном (газо-водонапорном) режиме точность расчета этим методом резко снижается. Для проверки возможности применения метода подсчета по падению давления рассчитывают количества добытого из залежи газа на единицу падения давления в разные периоды разработки. Если результаты этих расчетов совпадают, можно применить метод по падению давления. Увеличение количества добытого газа на единицу снижения давления в более поздние периоды разработки указывает на наличие напора вод и вытеснение ими части объема газа.  Q– количество добытого газа за период снижения давления от РН до Р в м3, РН – начальное абсолютное пластовое давление в атм., Р – среднее пластовое давление, после отбора Q газа, a1 и a2 – поправки на сжимаемость соответственно для РН и Р. Режим газовой залежи устанавливается с помощью графика, построенного в координатах иQг. Если фактические точки образуют прямую то в залежи газовый режим. Продолжив прямую до пересечения с осью Qг, получим начальные запасы газа в залежи, приведённые к стандартным условиям. Внедрение воды в залежь отражается в постепенном уменьшении угла наклона прямой вследствие, замедления темпа падения давления.  Однако бывают случаи, когда работа залежей на упруговодонапорном режиме характеризуется прямолинейной зависимостью, что обусловлено темпом отбора газа. Поэтому необходим комплексное исследование по контролю за внедрением воды в залежь. Оно включает наблюдение за изменением пластового давленияв пьезометрических скважинах. Снижение начального пластового давления свидетельствует о распространении воронки депрессии за границы залежи. На крупных месторождениях бурят контрольные скважины, где проводят геофизические исследования. Активность краевых вод при разработке может явиться причиной утечки газа. При подсчёте запасов газа методом падения предусматривается необходимость обоснования и расчёта - начальных и текущих пластовых давлений и температур; - начального и текущего положения ГВК; - изменение во времени устьевых и пластовых давлений; -установление гидродинамических связей между залежами месторождения; - степень дренируемости отдельных участков; - режимы работы залежи и отдельных участков залежи; - динамика вторжения пластовых вод в залежь; - потери газа при фонтанировании и исследовании; -перетоки газа; - величина отбора газа, конденсата и воды. На величину запасов газа влияет детальность замеров. Закон Дарси Рпл*Vпл= Ратм*Vг атм Q = V/T Рпл* Q пл= Ратм* Q г атм Q = Кпр*S (P12 –P22) 2L *μ*Pатм V= Q = Кпр*(P1 –P2) для жидкости S μ* L Линейный закон фильтрации Дарси Q = C*H1-H2* ώ L Ώ –площадь поперечного сечения Q – объёмный расход H1-H2 – потеря напора по линии напоры в начальном и конечном сечении образца пористой среды L – длина образца С – коэффициент фильтрации, зависит от свойств пористой среды и от свойств жидкости Q =k *(P1 –P2)* ώ μ* L P1 и P2 – давления приведённые к плоскости отсчёта геометр. высот (перепад давлений) к – коэффициент проницаемости зависит от свойств пористой среды μ- вязкость k = С μ ρg 1 Дарси = 1,02*10-12 м2 К= Д, С= м/с Список литературы М.М. Иванова, И.П. Чоловский, Ю.И. Брагин Нефтегазопромысловая геология. М., Недра, 2000, 414 с. |