СПОСОБ ОПРЕДЕЛЕНИЯ ПРОФИЛЯ ПРИТОКА И ПАРАМЕТРОВ ОКОЛОСКВАЖИННОГО ПРОСТРАНСТВА В МНОГОПЛАСТОВОЙ СКВАЖИНЕ Российский патент 2014 года по МПК E21B47/10 E21B47/06 

Описание патента на изобретение RU2505672C1

Область техники

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

Предшествующий уровень техники

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

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

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

Так, если величина скин-фактора близка к нулю (практически с учетом погрешности определения: - 1<S<1), то приствольная зона пласта считается неизмененной, а скважина совершенной. Большая положительная величина скин-фактора S>1 свидетельствует о загрязнении приствольной зоны пласта и несовершенстве скважины, что требует проведения дополнительных мероприятий по интенсификации притока (дополнительная перфорация, гидроразрыв пласта и др.). Значительная отрицательная величина скин-фактора S<-1 наблюдается в случае повышенной проницаемости приствольной зоны пласта (трещины, каверны и т.д.).

В настоящее время известны способы гидродинамического исследования скважины (Бузинов С.Н., Умрихин И.Д. Исследование нефтяных и газовых скважин и пластов. М., Недра, 1984), позволяющие определить в том числе и величину скин-фактора. Однако эти способы, как правило, определяют величину среднего скин-фактора сразу для нескольких работающих продуктивных пластов и не позволяют установить факторы, определяющие величину этого скин-фактора, т.е. оценить параметры зоны перфорации и параметры загрязнений призабойной зоны.

С другой стороны, из уровня техники известно (см. Чекалюк Е.В. Термодинамика нефтяного пласта. М., Недра, 1965, стр.238), что температура втекающего в скважину флюида, даже из изначально изотермического пласта, изменяется со временем (в технической литературе этот эффект именуется нестационарным эффектом Джоуля-Томсона). Информация о скорости изменения температуры притока со временем может быть использована для определения параметров загрязненной околоскважинной области с пониженной проницаемостью (см. Yu.A.Popov, V.Р.Pimenov, V.V.Tertychnyi, Developments of Geothermal Investigations of Oil and Gas Fields, Oilfield review, spring 2001, pp.4-11).

Однако для всех продуктивных пластов, кроме самого нижнего, эту информацию практически невозможно получить из данных по температуре, получаемой при каротаже в стволе скважины. Кроме того, теория, приведенная в (см. Чекалюк Е.В. Термодинамика нефтяного пласта М., Недра, 1965, стр.238), справедлива только для цилиндрически симметричного потока и напрямую не может быть использована для сложной пространственной структуры потока флюида в зоне перфорации.

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

Раскрытие изобретения

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

Заявленный технический результат достигается тем, что изменяют дебит скважины и осуществляют измерение во времени температуры втекающего в скважину флюида для каждого пласта, определяют величину изменения температуры ΔTp начального этапа, и величину установившегося значения А логарифмической производной температуры от времени для каждого пласта, величину удельного дебита q каждого пласта определяют по формуле:

,

где ε0 - коэффициент Джоуля-Томсона, η - вязкость флюида, π - математическая константа, равная 3,14159, ke- проницаемость пласта, долю длины перфорационного канала, через которую осуществляется приток флюида каждого пласта, рассчитывают по формуле:

,

где a - безразмерный параметр, определяемый по формуле:

,

где ΔP - разность между давлением в пласте и давлением в скважине, , rp - радиус перфорационного канала, nр - количество перфорационных каналов на 1 метр длины скважины, Lp - длина перфорационного канала, re - радиус внешней границы пласта, rc - радиус внутренней границы потока цилиндрически симметричного относительно оси скважины, рассчитываемый из времени ts, спустя которое устанавливается постоянное значение логарифмической производной температуры по формуле:

,

где , ϕ - пористость, ρmсm и ρfcf - объемные теплоемкости материала скелета горной породы (в отсутствие поровых флюидов) и флюида соответственно,

причем дебит Q каждого пласта скважины определяют по формуле Q=q·h, где h - толщина пласта, а профиль притока скважины определяют как совокупность дебитов Q всех пластов.

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

Устройство для измерения параметров флюида притока (Фиг.1) состоит из корпуса 1, который содержит заднюю стенку 4, крышку 6 и боковые стенки 5, которые выполнены сплошными и герметично соединенными между собой. Задняя стена 4 выполнена из слабодеформируемого, теплоизолирующего материала; крышка 6 и боковые стенки 5 выполнены из упругодеформируемого теплоизолирующего материала, например, упругого пластика или из резины. Корпус 1 может быть дополнительно снабжен днищем 7 с как минимум одним отверстием 8. В этом случае, крышка 7 также герметично соединяется с задней 4 и боковыми стенками 5 корпуса 1.

Внутри корпуса 1 расположена объемная сетчатая структура измерительных ячеек 2, выполненная также из упругодеформируемого теплоизолирующего материала. Сетчатая структура имеет равную высоту измерительных ячеек. Каждая измерительная ячейка содержит датчик температуры 3 и является независимым измерительным элементом. Устройство снабжено прижимным элементом 9.

Устройство работает следующим образом.

Устройство опускают в скважину 10 и располагают на уровне продуктивного пласта в зоне поступления флюида притока в скважину 10 и осуществляют плотное прижатие устройства к стенке 11 обсадной колонны 10 с помощью прижимного элемента 9. В качестве прижимного элемента 9 может быть использована, например, пневмопружина или иное устройство с пневмо-, гидро- или электроприводом, способное осуществить плотное прижатие устройства к внутренней поверхности обсадной колонны 10. Радиус кривизны устройства в прижатом к стенке 11 обсадной колонны 10 состоянии совпадает с внутренним радиусом кривизны колонны 10.

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

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

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

Выполнение стенок корпуса из теплоизолирующего материала позволяет исключить влияние процессов передачи тепла от основного потока в скважине на точность измерения температуры флюида притока.

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

Сетчатая структура представляет из себя, по сути, набор равных по высоте полых ячеек (в случае, например, ее прямоугольной формы, когда структура может быть выполнена в виде решетки, образованной перпендикулярно-пересекающимися перегородками, в пространственном, объемном выражении это - параллелепипед). Форма измерительных ячеек и, следовательно, форма сетчатой структуры может быть достаточно многообразной. Измерительные ячейки могут иметь форму n-угольника, где n-количество его углов (квадрат, пентагон, гексагон и др.). В любом случае их размеры должны быть рассчитаны исходя из соразмерности с перфорационными отверстиями, а количество ячеек (размеры сетчатой структуры) - с учетом параметров зоны перфорации.

Так, постоянная времени Δt устройства определяется объемом измерительной ячейки Vcell средним дебитом отдельного перфорационного канала qp:

где ΔР разность между давлением в пласте и давлением в скважине, π· - математическая константа, равная 3,14159, η - вязкость флюида, nр - число перфорационных каналов на 1 метр длины скважины в зоне перфорации, kе - проницаемость пласта, rw и re - радиус скважины и радиус внешней границы пласта.

Для реализации заявленного способа определения параметров скважины изменяют дебит скважины и для каждого продуктивного пласта измеряют изменение температуры Tin втекающего в скважину флюида во времени Tin=Tin(t).

Отличительной особенностью зависимости Tm(t) при нестационарном режиме работы проперфарированной скважины является быстрое изменение температуры на начальном этапе (первые 20-60 минут) и медленное изменение (в десятки раз медленнее, чем на начальном этапе) температуры после 3-5 часов (Фиг.4).

Из измеренной зависимости Tin(t) определяют величину изменения температуры ΔTр начального этапа (Фиг.4).

Осуществляют построение зависимости Тin(t) от логарифма времени t, прошедшего после начала добычи или изменения дебита скважины (Фиг.5). Находят время ts, спустя которое угол наклона полученной кривой (т.е. логарифмическая производная ) становится постоянным:

и находят установившееся значение логарифмической производной А. Величина этой постоянной связана с удельным дебитом скважины q[м3/(м·с)] (приток на единицу длины скважины) и соответствующей величиной скин-фактора s соотношениями:

После этого, по величине ts из соотношения (5) можно оценить расстояние rs от оси скважины, где (r>rs) течение флюида можно считать цилиндрически симметричным.

Величина rs, как показывают численные расчеты, пропорциональна длине Lр перфорационного канала.

Далее по величине изменения температуры ΔTp начального этапа находят безразмерный параметр:

который используют, в свою очередь, для расчета параметра δ, который дает оценку длины перфорационного канала Lp_act, через которую осуществляется приток флюида в скважину (δ=Lp_act/Lp):

где , rp и Lp - радиус перфорационного канала.

Значение δ=1 означает, что никакой загрязненной зоны вокруг скважины нет и перфорационный канал работает по всей длине. При известной длине перфорационных каналов значения параметра δ в интервале от 0 до 1 представляют количественную оценку доли длины перфорационного канала, через которую осуществляется приток флюида в скважину. В случае неизвестной длины Lp перфорационных каналов из той же формулы (7) определяют эффективную длину перфорационного канала Lpδ, которая характеризует одновременно размеры загрязненной зоны в околоскважинной области и длину каналов перфорации.

Далее, дебит Q каждого пласта скважины определяют по:

где h - толщина (мощность) пласта. Профиль притока скважины определяют как совокупность дебитов Q всех пластов.

Краткое описание чертежей

На фиг.1 изображено устройство для измерения температуры втекающего в скважину флюида, горизонтальный разрез и вид устройства со стороны обсадной колонны; фиг.2 - схематичное изображение расчетной области для моделирования течения флюида в пористой среде, питающей один перфорационный канал; на фиг.3 - расчетная область (сетка) для реализации 3D моделирования течения в пористой среде, питающей 1/4 перфорационного канала; на фиг.4 - зависимость увеличения температуры втекающего в скважину флюида от времени, прошедшего с момента начала добычи; на фиг.5 - зависимость температуры втекающего в скважину флюида от натурального логарифма времени (в секундах).

Варианты осуществления изобретения

Реализация заявленного способа определения профиля притока и параметров зоны перфорации была осуществлена путем 3D численного моделирования процессов тепло- и массопереноса в пористой среде в зоне перфорации с использованием коммерческого программного обеспечения для гидродинамического и теплового моделирования STAR-CD, разработанного компанией CD-ADAPCO. Моделировалось течение флюида в области пористой среды, питающей 1/4 перфорационного канала (Фиг.2, Фиг.3).

При моделировании использовались следующие расчетные параметры. Параметры зоны перфорации: h=0.4 м - расстояние по высоте и φ=60° - угол между перфорационными каналами, np=15 - число перфорационных каналов, приходящихся на 1 метр длины скважины, rр=0.01 м Lp=0.5 м - радиус и длина перфорационного канала. Размеры расчетной области - высота 0.2 м и угол 30° (Фиг.3), rw=0.1 м - радиус скважины, rе=10 м - внешний радиус расчетной области, Pw=50 бар - давление в скважине, Ре=100 бар - давление на внешней границе расчетной области, ke=10-13 м2 проницаемость пласта, 0.3·ke - проницаемость пласта в слое толщиной 0.2 м вокруг скважины, φ=0.2 - пористость пласта, cf=1800 Дж/(кг·К) и cm=750 Дж/(кг·К) - удельные теплоемкости пластового флюида и породы, ρf=900 кг/м3 и ρm=2700 кг/м3 - плотности пластового флюида и породы, ε0=-5·10-7 К/Па - коэффициент Джоуля-Томсона, η=0.01 Па·с - вязкость флюида.

В результате численного моделирования была получена величина массового притока в скважину (1/4 притока перфорационного канала) gc=1.21·10-3 кг/с. Объемный дебит qc на 1 м длины скважины рассчитывался по формуле:

Для проверки предлагаемого способа определения локальной величины притока и параметров зоны перфорации эта величина была найдена по приведенным на Фиг.4, Фиг.5 зависимостям температуры притока от времени.

Из Фиг.5 видно, что спустя время ts≈e10 с≈6 часов логарифмическая производная температуры становится равной .

Подставив эту величину в формулу (8), определили:


Эта величина qT практически совпадает с приведенным выше значением qc.

Далее, по формуле (9) и величине ts находим внутренний радиус rs области, в которой течение в пласте можно считать цилиндрически симметричным: rs=0.63 м. Эта величина хорошо согласуется с размерами скважины и длиной перфорационного канала rw+rs=0.6 м.

Из Фиг.4 видно, что в течение первых 10-20 минут температура притока увеличивается приблизительно на ΔТр=0.5 К ("быстрая" стадия роста температуры на начальном этапе), а далее температура меняется относительно медленно. Подставив эту величину в формулу (10), находим а=5 и далее, из формулы (11) находим δ≈0.8, т.е. приток осуществляется через 80% от длины перфорационного канала, что соответствует существованию (в расчетной модели) слоя с пониженной проницаемостью вокруг скважины.

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

Для того чтобы в зону измерения попало не менее 3 перфорационных каналов, площадь сетки измерительных ячеек должна превышать 3 Sp. Если принять, что сетка измерительных ячеек перекрывает 30% периметра скважины (длина дуги lf=0.2 м), ее высота ht должна составлять:

Промышленная применимость

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

Похожие патенты RU2505672C1

название год авторы номер документа
СПОСОБ ОПРЕДЕЛЕНИЯ ПРОФИЛЯ ПРИТОКА ФЛЮИДОВ И ПАРАМЕТРОВ ОКОЛОСКВАЖИННОГО ПРОСТРАНСТВА 2010
  • Шако Валерий Васильевич
  • Пименов Вячеслав Павлович
  • Кучук Фикри Джон
RU2460878C2
УСТРОЙСТВО ДЛЯ ИЗМЕРЕНИЯ ПАРАМЕТРОВ ФЛЮИДА ПРИТОКА СКВАЖИНЫ 2009
  • Шако Валерий Васильевич
  • Пименов Вячеслав Павлович
RU2418947C1
СПОСОБ ОПРЕДЕЛЕНИЯ ПРОФИЛЯ ПРИТОКА ФЛЮИДОВ И ПАРАМЕТРОВ ОКОЛОСКВАЖИННОГО ПРОСТРАНСТВА 2010
  • Шако Валерий Васильевич
  • Пименов Вячеслав Павлович
  • Кучук Фикри Джон
RU2455482C2
СПОСОБ ОПРЕДЕЛЕНИЯ ПРОФИЛЯ ПРИТОКА ФЛЮИДА В МНОГОПЛАСТОВОЙ СКВАЖИНЕ 2016
  • Шако Валерий Васильевич
  • Пименов Вячеслав Павлович
  • Котляр Лев Андреевич
RU2645692C1
Способ назначения нагнетательных и добывающих скважин и изменения их интервалов перфораций 2017
  • Родионов Сергей Павлович
  • Пичугин Олег Николаевич
  • Мусакаев Эмиль Наильевич
  • Косяков Виталий Петрович
RU2658422C1
СПОСОБ ОПРЕДЕЛЕНИЯ ПРОФИЛЯ ПРИТОКА ФЛЮИДОВ МНОГОПЛАСТОВЫХ ЗАЛЕЖЕЙ В СКВАЖИНЕ 2013
  • Шако Валерий Васильевич
  • Пименов Вячеслав Павлович
  • Тевени Бертран
  • Сидорова Мария Викторовна
RU2531499C1
Система долговременного распределенного мониторинга профиля притока в горизонтальной скважине, оборудованной ЭЦН 2019
  • Яковлев Андрей Александрович
  • Сулейманов Айяр Гусейн Оглы
  • Файзуллин Ильдар Гаязович
  • Ипатов Андрей Иванович
  • Кременецкий Михаил Израилевич
  • Шурунов Андрей Владимирович
  • Сарапулов Николай Павлович
  • Симаков Сергей Михайлович
RU2703055C1
Способ и система моделирования трещин гидроразрыва пласта бесконечно-конечной проводимости и поперечно-продольного расположения относительно горизонтального ствола скважины 2020
  • Коваленко Игорь Викторович
RU2745142C1
СПОСОБ ОПРЕДЕЛЕНИЯ ПРОФИЛЯ ПРИТОКА ФЛЮИДОВ МНОГОПЛАСТОВЫХ ЗАЛЕЖЕЙ 2011
  • Шако Валерий Васильевич
  • Пименов Вячеслав Павлович
  • Тевени Бертран
RU2474687C1
СПОСОБ ОБРАБОТКИ ПРИЗАБОЙНОЙ ЗОНЫ СКВАЖИНЫ 2010
  • Хисамов Раис Салихович
  • Шафигуллин Ринат Ильдусович
  • Чупикова Изида Зангировна
  • Афлятунов Ринат Ракипович
  • Козихин Роман Анатольевич
  • Мухамадеева Лилия Наилевна
RU2425960C1

Иллюстрации к изобретению RU 2 505 672 C1

Реферат патента 2014 года СПОСОБ ОПРЕДЕЛЕНИЯ ПРОФИЛЯ ПРИТОКА И ПАРАМЕТРОВ ОКОЛОСКВАЖИННОГО ПРОСТРАНСТВА В МНОГОПЛАСТОВОЙ СКВАЖИНЕ

Изобретение относится к области геофизических исследований нефтяных и газовых скважин и может быть использовано, в частности, при определении профиля притока скважины и параметров околоскважинного пространства. Согласно способу изменяют дебит скважины и осуществляют измерение во времени температуры втекающего в скважину флюида для каждого пласта, определяют величину изменения температуры ΔТ р начального этапа, и величину установившегося значения А логарифмической производной температуры от времени для каждого пласта. Величину удельного дебита q каждого пласта определяют по приведенному математическому выражению. Определяют дебит Q каждого пласта скважины, а профиль притока скважины определяют как совокупность дебитов Q всех пластов. Техническим результатом является повышение точности определения параметров скважины. 5 ил.

Формула изобретения RU 2 505 672 C1

Способ определения профиля притока и параметров околоскважинного пространства в многопластовой скважине, при котором изменяют дебит скважины и осуществляют измерение во времени температуры втекающего в скважину флюида для каждого пласта, определяют величину изменения температуры ΔТ р начального этапа, и величину установившегося значения А логарифмической производной температуры от времени для каждого пласта, величину удельного дебита q каждого пласта определяют по формуле:
,
где ε 0 - коэффициент Джоуля-Томсона, η - вязкость флюида, π - математическая константа, равная 3,14159, k e - проницаемость пласта, долю длины перфорационного канала, через которую осуществляется приток флюида каждого пласта, рассчитывают по формуле:
,
где a - безразмерный параметр, определяемый по формуле:
,
где ΔP - разность между давлением в пласте и давлением в скважине, , r p - радиус перфорационного канала, n p - количество перфорационных каналов на 1 метр длины скважины, L p - длина перфорационного канала, r e - радиус внешней границы пласта, r s - радиус внутренней границы потока цилиндрически симметричного относительно оси скважины, рассчитываемый из времени t s, спустя которое устанавливается постоянное значение логарифмической производной температуры по формуле:
,
где , ϕ - пористость, ρ m с m и ρ f c f - объемные теплоемкости материала скелета горной породы и флюида соответственно, причем дебит Q каждого пласта скважины определяют по формуле Q=q·h, где h - толщина пласта, а профиль притока скважины определяют как совокупность дебитов Q всех пластов.

Документы, цитированные в отчете о поиске Патент 2014 года RU2505672C1

УСТРОЙСТВО ДЛЯ ИЗМЕРЕНИЯ ВНУТРИСКВАЖИННЫХ ПАРАМЕТРОВ 2003
  • Осадчий В.М.
  • Леонов В.А.
  • Перегинец В.А.
  • Полыгалов В.Ф.
  • Шарифов Махир Зафар Оглы
  • Мусаверов Р.Х.
  • Гарипов О.М.
  • Синёва Ю.Н.
  • Мокрый М.В.
RU2249108C1
СПОСОБ ИССЛЕДОВАНИЯ НАГНЕТАТЕЛЬНЫХ СКВАЖИН 1997
  • Назаров В.Ф.
  • Валиуллин Р.А.
  • Адиев Я.Р.
  • Азизов Ф.Ф.
RU2121572C1
Способ определения профиля притока флюида в действующей газовой скважине и устройство для его осуществления 1986
  • Дивеев Исмаил Исхакович
  • Кавтанюк Владимир Захарович
  • Макушев Федор Иванович
  • Глазов Георгий Константинович
  • Кондратьев Дмитрий Венидиктович
  • Голубев Игорь Александрович
  • Ермаков Геннадий Владимирович
SU1421858A1
US 5551287 A, 03.09.1996.

RU 2 505 672 C1

Авторы

Шако Валерий Васильевич

Пименов Вячеслав Павлович

Даты

2014-01-27Публикация

2009-12-31Подача