Разработка системы прогнозирования катастрофических паводков в Восточно-Казахстанской области с применением данных ДЗЗ, ГИС-технологий и машинного обучения

Руководитель программы: Алия Нугуманова, доктор PhD по специальности «Информационные системы», директор научно-инновационного центра “Big Data and Blockchain Technologies” в Astana IT University. Выпускница Восточно-Казахстанского университета имени Сарсена Аманжолова. Защитила диссертацию по теме больших данных и NLP в Казахском национальном исследовательском университете имени К. Сатпаева в 2014 году.

Наименование конкурса: ПЦФ МНВО РК​

Годы реализации: 2024–2026

Объем финансирования: 639, 4 млн тенге​

Цель программы

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

Актуальность программы

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

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

Рисунок 1 – Пространственное распределение метеорологических станций по территории ВКО, сформированная в программной среде ArcGIS.

Рисунок 2 – Расположение гидрологических постов в пределах территории ВКО

Задачи программы

  1. Исследование, анализ и оценка паводковой ситуации на территории ВКО на основе многолетних данных об осадках, снегозапасах и процессах снеготаяния
  2. Дистанционное исследование территорий водосбора талых вод и зон затопления по материалам космической съемки и съемки с беспилотных летательных аппаратов (БПЛА)
  3. Создание цифровых моделей территорий водосбора талых вод и зон затопления ВКО
  4. Создание цифровой модели гидрографической сети ВКО
  5. Расчет аккумуляции талых вод, определение технологий и методов для их сбора
  6. Определение территорий, наиболее подверженных затоплениям в результате снеготаяния
  7. Создание методики расчёта и прогнозирования территорий затопления с использованием современных методов радарной космической съемки, съемки в оптическом (мульти- и гипер-спектральных) диапазонах, спутниковой альтиметрии и съемки с БПЛА
  8. Разработка согласованной базы геореляционной структуры и гидрографических данных на территорию Восточно-Казахстанской области на основе множественно-реляционной модели (ГИС-сервис)
  9. Рекомендации по использованию талых вод для нужд сельского хозяйства, орошения, энергетики и других видах водопользования

Ожидаемые результаты

  1. Публикации

По результатам программы будут опубликованы:  

1) не менее 9 (девяти) статей и (или) обзоров в рецензируемых научных изданиях по научному направлению программы, входящих в 1 (первый), 2 (второй) и (или) 3 (третий) квартиль по импакт-фактору в базе данных Web of Science и (или) имеющих процентиль по CiteScore в базе данных Scopus не менее 50 (пятидесяти). 

2) не менее 10 (десять) статей в журналах, рекомендованных КОКНВО. 

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

4) не менее 2 (двух) патентов в зарубежных патентных бюро (европейском, американском, японском) или не менее 2 (двух) зарубежных или международных патентов, включенных в базу данных Derwent Innovations Index (Web of Science, Clarivate Analytics) либо не менее 5 (пяти) объектов интеллектуальной собственности (патент; для заявок в области информационных технологий — авторское свидетельство), зарегистрированных в Национальном Институте интеллектуальной собственности Республики Казахстан. 

  1. Прямые результаты
  • По результатам программы будут получены следующие прямые результаты:
  • Рекомендации по использованию талых вод для нужд сельского хозяйства, орошения, энергетики и других видах водопользования.
  • Научно-обоснованная методика дистанционного исследования территорий затопления по материалам космической съемки и с беспилотных летательных аппаратов.
  • Прогнозная модель расчёта и моделирования территорий затопления с использованием современных методов радарной космической съемки, съемки в оптическом (мульти и гиперспектральных) диапазонах, спутниковой альтиметрии и съемки с БПЛА.
  • Цифровая модель территории ВКО, в том числе территории водосбора и зон затопления в ВКО;
  • Комплексная информационно-аналитическая система сбора, обработки, анализа данных по материалам космической съемки и съемки с БПЛА, разработка согласованной базы геореляционной структуры и гидрографических данных на территорию Восточно-Казахстанской области на основе множественно-реляционной модели (ГИС сервис).
  • Научно-методическая база и техническая документация по методике сбора и представления первичных данных;
  • Прогнозная модель территорий затопления в результате таяния снега и аккумуляции талых вод по ВКО.
  • Расчет возможных аккумуляционных запасов талых и паводковых вод на территории Восточно-Казахстанской области.

Рисунок 3– Положение сезонной снеговой границы в бассейне р. Уба в 2021 г. по данным спутников Sentinel 

Достигнутые результаты

  1. Публикации

Scopus:

  1. Hybrid feature-based neural network regression method for load profiles forecasting

В данном исследовании рассматривается острая необходимость в усовершенствованных моделях прогнозирования спроса, способных точно предсказывать потребление энергии, особенно с учётом различных географических и климатических условий. В работе предлагается новая модель прогнозирования спроса, которая объединяет методы кластеризации и инженерии признаков с регрессией на основе нейронных сетей, при этом особое внимание уделяется учёту корреляции с температурой воздуха. Эффективность модели оценивалась на основе эталонного набора данных из Тетуана (Марокко), где существующие методы прогнозирования демонстрировали значения RMSE в диапазоне от 6429 до 10 220 [МВт⋅ч]. В отличие от них, предлагаемая модель достигла существенно более низкого значения RMSE — 5168, что свидетельствует о её превосходстве. Последующее применение модели для прогнозирования спроса в Астане (Казахстан) в рамках тематического исследования дополнительно подтвердило её эффективность. Сравнительный анализ с базовой нейронной сетью показал заметное улучшение: предложенная модель достигла значения MAPE в 5,19%, тогда как у базовой модели этот показатель составил 17,36%. Эти результаты подчёркивают потенциал предлагаемого подхода в повышении точности прогнозирования спроса, особенно в разнообразных географических контекстах, за счёт использования климатических факторов. Методология также демонстрирует перспективность для более широких применений, таких как прогнозирование наводнений, оценка урожайности сельскохозяйственных культур и управление водными ресурсами.

 

Satan A. et al. Hybrid feature-based neural network regression method for load profiles forecasting //Energy Informatics. – 2025. – Т. 8. – №. 1. – С. 19. https://doi.org/10.1186/s42162-025-00481-0

 

КОКНВО: 

  1. Flood Risk Mapping in The Irtysh River Basin Using Satellite Data

Наводнения являются одними из самых частых и разрушительных природных катастроф, вызывая значительный экономический ущерб и человеческие жертвы по всему миру. Эффективное управление рисками наводнений зависит от точных методов моделирования, способных предсказывать уязвимые территории и оценивать потенциальные последствия. В данном исследовании моделируется динамика наводнений в бассейне реки Иртыш вблизи города Усть-Каменогорска (Восточный Казахстан), подверженного сезонным паводкам, с использованием спутниковых снимков высокого разрешения и цифровых моделей рельефа. Основной целью работы является визуальное моделирование рисков наводнений на основе характеристик рельефа. В исследовании используются спутниковые изображения, предоставляемые платформой Mapbox, которая объединяет данные MODIS, Landsat 7, Maxar и Google Earth Engine, обеспечивая доступ к изображениям Sentinel-2 с отражательной способностью поверхности при разрешении 10 метров. Для моделирования затоплений применяются данные о высотах из глобальной цифровой модели рельефа Copernicus с разрешением 30 метров.

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

Rakhymbek K., Zhomartkan N., Nurekenov D., Zhantassova Z. Flood Risk Mapping in The Irtysh River Basin Using Satellite Data //Scientific Journal of Astana IT University. – 2024. – Т. 19. – С. 140-149. https://doi.org/10.37943/19LRYW4856 

 

  1. High-Resolution Satellite Estimation of Snow Cover for Flood Analysis in East Kazakhstan Region

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

Для достижения этой цели используется мультиспектральная спутниковая съёмка высокого разрешения из набора данных Sentinel-2 Surface Reflectance. Анализ охватывает изображения, собранные в период с марта по октябрь за 2021–2024 годы. Обработка данных осуществляется в платформе Google Earth Engine с применением строгой фильтрации на основе пространственного пересечения с исследуемыми бассейнами и доли облачных пикселей, что обеспечивает высокое качество данных для анализа снежного покрова. В исследовании применяются несколько дистанционных индексов для оценки снежного покрова. Нормализованный дифференциальный снежный индекс (NDSI) рассчитывается с использованием зелёного и коротковолнового инфракрасного каналов для определения пикселей, покрытых снегом. Доля снежного покрова (fractional snow-covered area, fSCA) вычисляется на основе NDSI с применением эмпирического уравнения ‘FRA6T’, что позволяет получить более детальное представление о распределении снега по территории бассейнов. Дополнительно применяется пороговое значение отношения ближнего инфракрасного диапазона к коротковолновому инфракрасному, что минимизирует путаницу между снегом и водой, особенно вблизи водоёмов и в периоды активного таяния. Полученные карты снежного покрова и значения fSCA обеспечивают подробное представление о распределении снега и динамике его таяния, способствуя оценке роли снежного стока в развитии паводковых рисков. Полученные выводы могут использоваться для уточнения моделей прогнозирования наводнений, совершенствования систем раннего предупреждения и поддержки обоснованного управления водными ресурсами в уязвимых регионах.

Alzhanov A., Nugumanova A. High-Resolution Satellite Estimation of Snow Cover for Flood Analysis in East Kazakhstan Region //Scientific Journal of Astana IT University. – 2024. – Т. 19. – С. 118-127. https://doi.org/10.37943/19VUAO6399. 

 

  1. Interrelationships between snowpack dynamics and tree growth in the Tigiretsky Ridge (Altai): Implications for ecological responses to climate variability

 

В данном исследовании рассматривается сложное взаимодействие между динамикой снежного покрова и ростом деревьев в районе Тигирецкого хребта, с использованием дендрохронологических данных и данных о снежном покрове за период с 2013 по 2020 год. С применением метода температурного индекса таяния (Temperature-Based Melt-Index Method) была проведена точная оценка максимального водного эквивалента снежного покрова в зимние месяцы, что позволило выявить значительную пространственную изменчивость, обусловленную высотой, экспозицией склонов и близостью к водоразделам. Результаты исследования указывают на асимметричное распределение снежного покрова: на южных склонах на низких высотах наблюдаются более высокие запасы снега, тогда как на больших высотах ситуация меняется. Особенно примечательно, что вблизи водоразделов запасы снега на северных склонах могут превышать запасы на южных в 30 раз. Анализ также выявил положительную корреляцию между увеличением водного эквивалента снежного покрова и радиальным приростом Abies sibirica L. (пихты сибирской) в экотоне верхней границы леса, что указывает на значительный экологический отклик деревьев на изменение снежных условий. Полученные данные способствуют более глубокому пониманию воздействия климатической изменчивости на взаимодействие между снежным покровом и растительностью в горных экосистемах, формируя основу для дальнейших исследований, направленных на раскрытие механизмов этих взаимосвязей.

 

Bykov N. I., Birjukov R. J. Interrelationships between snowpack dynamics and tree growth in the Tigiretsky Ridge (Altai): Implications for ecological responses to climate variability // Acta Biologica Sibirica. – 2024. –V. 10. – P. 1319–1336. https://doi.org/10.5281/zenodo.14190443 (журнал индексируется в Scopus по направлению Environmental Science > Ecology; Global and Planetary Change, процентиль журнала — 29). 

 

  1. Система мониторинга и обеспеченность данными прогноза наводнений Восточно-Казахстанской области

Данная статья посвящена сбору доступной информации о водных объектах Восточно-Казахстанской области, полученной из открытых источников — данных, предоставленных Республиканской гидрометеорологической службой и Министерством по чрезвычайным ситуациям Республики Казахстан. Целью работы является создание пространственно привязанной базы атрибутивных данных в рамках ГИС-платформы, а также формирование основы для дальнейших исследований, направленных на разработку и внедрение моделей прогнозирования наводнений. Проведён обзор существующих исследований в области прогнозирования паводков как по регионам Казахстана, так и за рубежом, что позволило сформулировать критерии оценки данных и систем гидрологического мониторинга. В статье также дана характеристика гидрологического режима ключевых репрезентативных водных объектов Восточного Казахстана. На основе пространственного анализа сети гидрологических станций и известных зон затопления была проведена оценка пространственного охвата государственной системы наблюдений, по результатам которой разработаны рекомендации по её расширению. В заключение сделаны выводы о применимости собранных данных для построения прогностических моделей.

Павленко А.В., Мансурова А.К., Кызырканов А., Черных Д.В. Система мониторинга и обеспеченность данными прогноза наводнений Восточно-Казахстанской области //Вестник Карагандинского Университета. Серия Биология. Медицина. География. – 2024. – Т. 4. https://doi.org/10.31489/2024bmg4/183-196

 

  1. Feature selection methods for lstm-based river water level and discharge forecasting

Точное прогнозирование речного стока и уровней воды имеет решающее значение для эффективного управления водными ресурсами, смягчения последствий наводнений и обеспечения общественной безопасности. В данном исследовании проведено сравнение методов отбора признаков на основе корреляционного анализа и метода главных компонент (PCA) для моделей прогнозирования на базе LSTM в бассейне реки Уба, в пределах города Шемонаиха Восточно-Казахстанской области. Исходный набор данных охватывает период с 1995 по 2021 год, при этом данные за 1995–2019 годы использовались для обучения и валидации моделей, а данные за 2020–2021 годы — для тестирования. Оба метода отбора признаков сократили исходный набор предикторов до 13 переменных, при этом удалось в целом сохранить точность прогнозирования. Для повышения стабильности предсказаний и снижения дисперсии, связанной со случайной инициализацией, была обучена ансамблевая модель из 10 LSTM-сетей, использующих входные последовательности длиной 60 дней и осуществляющих прогноз на 10-дневный горизонт. Оценка производительности моделей проводилась с использованием коэффициента эффективности Нэша–Сатклиффа (Nash–Sutcliffe Efficiency, NSE). Результаты показали, что метод отбора признаков на основе корреляции дал сопоставимую точность с моделью, использующей полный набор признаков, при тестировании на данных 2020 года. Это говорит о том, что исключение высоко коррелированных признаков не снижает способность модели к краткосрочному прогнозированию. Модель с признаками, отобранными на основе PCA, продемонстрировала некоторое отставание на более длительных горизонтах в 2020 году, но показала преимущество на большинстве горизонтов в 2021 году. Однако общая точность прогнозов в 2021 году снизилась по сравнению с 2020 годом, что свидетельствует о большей изменчивости гидрологических условий и отклонении их от исторических данных обучения, указывая на необходимость периодического обновления модели с учётом новых данных. Оба метода отбора признаков эффективно снизили размерность данных, сохранив при этом предсказательную способность моделей. Однако ни один из методов не оказался универсально лучшим на всех временных горизонтах прогноза. Эти результаты подчеркивают важность систематического подхода к выбору признаков в гидрологических моделях, а также необходимость адаптации моделей к изменяющимся природным условиям.

Alzhanov A., Nugumanova A. FEATURE SELECTION METHODS FOR LSTM-BASED RIVER WATER LEVEL AND DISCHARGE FORECASTING //Scientific Journal of Astana IT University. – 2025. Т. 21. https://doi.org/10.37943/21EHLH9882

 

  1. Прямые результаты
  • Собран набор многолетних данных об осадках, снегозапасах, температурных режимах и солнечной радиации, водных объектах и набор спутниковых снимков водосборных территорий и зон затопления в ВКО.​
  • Разработаны цифровые модели рельефа водосборных территорий и зон затопления и концептуальная модель геореляционной базы данных водных ресурсов ВКО. ​
  • Разработана концептуальная модель геореляционной базы данных водных ресурсов ВКО на основе модели ODM (Observation Data Model). Для реализации использована СУБД PostgreSQL с расширением PostGIS, обеспечивающим поддержку геоданных и геопространственный анализ.​
  • Установлены 2 гидростанции и 2 метеостанции на территории ВКО.

 

Место реализации программы 

  1. Астана – Astana IT University.
  2. Восточно-Казахстанская область, Усть-Каменогорск – Восточно-Казахстанский университет имени Сарсена Аманжолова.

Рисунок 4 – Структура основных таблиц концептуальной модели для хранения метаданных 

Рисунок 5 – Структура основных таблиц концептуальной модели для хранения значений наблюдений

Команда программы

Ф.И.О., образование, степень, ученое звание 

Индекс Хирша, идентификаторы ResearcherID, ORCID, Scopus Author ID (при наличии) 

Роль в программе 

Нугуманова Алия Багдатовна, PhD по специальности «Информационные системы»

Индекс Хирша: Scopus: 6
ResearcherID: L-9616-2015
ORCID: 0000-0001-5522-4421
Scopus ID: 55864815200

Научный руководитель, Главный научный сотрудник

Муканова Балгайша Гафуровна, доктор физико-математических наук

Индекс Хирша: Scopus: 6, WoS: 5
ResearcherID: E-9945-2013
ORCID: 0000-0002-0823-6451
Scopus ID: 14828253100

Главный научный сотрудник

Черных Дмитрий Владимирович, АлтГУ (1994), доктор географических наук, доцент

Индекс Хирша: РИНЦ: 17, WoS: 3, Scopus: 6
ResearcherID: A-5129-2015
ORCID: 0000-0003-0151-2596
Scopus ID: 6666494700

Ведущий научный сотрудник

Гарцман Борис

Индекс Хирша: Scopus: 11
ORCID: 0000-0002-5876-7015
Scopus ID: 24438012800

Ведущий научный сотрудник

Жантасова Женискуль

Индекс Хирша: Scopus: 2
ORCID: 0000-0001-5550-7587
Scopus ID: 57197738233

Ведущий научный сотрудник

Морейдо Всеволод Михайлович, гидролог, кандидат географических наук

Индекс Хирша: WoS: 7
ResearcherID: D-7056-2014
ORCID: 0000-0003-1763-4096
Scopus ID: 56681224900

Ведущий научный сотрудник

Бондарович Андрей Александрович, АлтГУ (1993), кандидат географических наук, доцент

Индекс Хирша: РИНЦ: 7, WoS: 3, Scopus: 4
ResearcherID: AAH-5245-2020
ORCID: 0000-0001-9326-642X
Scopus ID: 57209656569

Ведущий научный сотрудник

Быков Николай Иванович, ТГУ (1984), кандидат географических наук, доцент

Индекс Хирша: РИНЦ: 15, WoS: 3, Scopus: 4
ResearcherID: HGT-8099-2022
ORCID: 0000-0003-3756-1859
Scopus ID: 35274333200

Ведущий научный сотрудник

Рахимжанова Анар Жанатовна, PhD

Индекс Хирша: Scopus: 1
ORCID: 0000-0003-4646-0603
Scopus ID: 57641574300

Старший научный сотрудник

Байбурин Ержан Мухаметкалиевич

Индекс Хирша: Scopus: 4
ResearcherID: GRT-5436-2022
ORCID: 0000-0002-1583-9912
Scopus ID: 56111999400

Старший научный сотрудник

Бирюков Роман

Индекс Хирша: Scopus: 6
ORCID: 0000-0002-5617-7206
Scopus ID: 57204692413

Старший научный сотрудник

Жакиев Нурхат

Индекс Хирша: Scopus: 5
ResearcherID: EBG-7696-2022
ORCID: 0000-0002-4904-2047
Scopus ID: 56043145000

Старший научный сотрудник

Қызырқанов Абзал Ермекбайұлы, магистр

Индекс Хирша: Scopus: 3
ResearcherID: AGL-3716-2022
ORCID: 0000-0002-8817-8617
Scopus ID: 57238488800

Научный сотрудник

Очередько Игорь

Индекс Хирша: Scopus: 4
ORCID: 0000-0001-8012-8589
Scopus ID: 57214718565

Научный сотрудник

Махамбетова Жансая Кайырбаевна

ORCID: 0000-0001-5024-0289

Научный сотрудник

Альжанов Алмас Миржанович, докторант

Индекс Хирша: Scopus: 1

ORCID: 0009-0007-8083-2366
Scopus ID: 58859587600

Научный сотрудник

Павленко Анатолий Владимирович, магистр естественных наук в области географии

ORCID: 0000-0001-8556-6633
Scopus ID: 59423326200

Научный сотрудник

Маулит Алмасбек, магистр технических наук

Индекс Хирша: Scopus: 3
ResearcherID: DFH-7960-2022
ORCID: 0000-0002-0519-3222
Scopus ID: 57220810948

Научный сотрудник

Нурекенов Даурен Махсутбекович, магистр технических наук

Индекс Хирша: Scopus: 1
ORCID: 0009-0003-5126-6319
Scopus ID: 57202929699

Научный сотрудник

Мансурова Айганым

Индекс Хирша: Scopus: 1

Scopus ID: 59233698800

ORCID 0009-0007-9076-0722

Научный сотрудник

Жомартқан Нұрасыл Қайратұлы, магистрант

ORCID: 0009-0006-3935-2013

Младший научный сотрудник

Рахымбек Камилла

ORCID: 0009-0008-7404-8433

Младший научный сотрудник

homescontents ataşehir escort ataşehir escort bostancı escort escort istanbul escort şişli escort istanbul eskort ataköy escort Marsbahis giriş Marsbahis ataşehir escort Marsbahis giriş Marsbahis küçükçekmece escort kadıköy escort marsbahis giris marsbahis casino marsbahis güncel adres marsbahis Şartsız deneme bonusu veren siteler Şartsız deneme bonusu veren siteler Deneme Bonusu Veren Siteler Yeni 2025 Deneme Bonusu Veren Siteler Deneme Bonusu Veren Siteler deneme bonusu veren siteler 2025 serifali eskort atasehir escort bayan bursa escort bursa eskort yenibosna escort umraniye escort teksert film izle film izle film izle film hd film sakarya escort sakarya escort
homescontents
https://www.fapjunk.com
gaziantep escort gaziantep escort
sakarya escort akyazı escort arifiye escort erenler escort eve gelen escort ferizli escort geyve escort hendek escort otele gelen escort sapanca escort söğütlü escort taraklı escort
sakarya escort akyazı escort arifiye escort erenler escort eve gelen escort ferizli escort geyve escort hendek escort karapürçek escort karasu escort kaynarca escort kocaali escort otele gelen escort pamukova escort sapanca escort söğütlü escort taraklı escort
Sakarya escort Sakarya escort Sakarya escort Sakarya escort Sakarya escort Sakarya escort Sapanca escort Sapanca escort Sapanca escort Sapanca escort Karasu escort
hd film izle