Jaiyq
Ғылыми-техникалық әдістеме · Jaiyq экологиялық мониторинг платформасы

JAIYQ-MRI

Тасқын-импульсті жұмыртқа банкіне негізделген маса тәуекел моделі
Атырау қаласы және Жайық өзенінің атырауы, Солтүстік Каспий

Толық аты: Jaiyq Flood-pulse Mosquito Risk Intelligence
Ядросы: FPEB — Flood-Pulse Egg-Bank engine
Түрі: механистік компартмент моделі + қашықтан зондтау
Валидация күйі: валидацияланбаған
Аумағы: Атырау облысы, 90 есептеу нүктесі
Жаңару жиілігі: сағат сайын
Құжаттың мәртебесі. Бұл — жұмыс істеп тұрған жүйенің техникалық сипаттамасы, жарнамалық материал емес. Модельдің әлі тексерілмегені, параметрлерінің калибрленбегені және қандай тұжырым жасауға жарамайтыны ашық жазылған. Әр сан қай дереккөзден келетіні көрсетілген.

1. Модель нені өлшейді және нені өлшемейді

Өлшейді: маса көбеюіне жағдайдың қаншалық қолайлы екенін — уақыт пен кеңістік бойынша салыстыруға жарайтын 0–100 индекс.

Өлшемейді: маса САНЫН. Индекс 70 деген жерде бір шаршы метрде қанша маса бар екенін бұл модель айта алмайды. Ол үшін далалық тұзақ деректері қажет, олар Атырау бойынша жинақталмаған.

Сондықтан индекс реттік (ordinal): «А ауданы Б ауданынан қауіптірек» деген тұжырым жасауға жарайды, «А ауданында 500 маса бар» деген тұжырымға жарамайды.

Неге бұл модель керек болды

Дайын модельдердің басым бөлігі ыдыс-су (container-breeding) масаларына — Aedes aegypti, Ae. albopictus — арналған. Олардың ажыратқышы: адам жинаған су (шина, ыдыс, науа). Ал Жайық атырауында басты драйвер мүлдем басқа: көктемгі қар суы тасқыны. Су басқанда құрғақ жатқан диапаузалы жұмыртқа банкі жарылады. Ыдыс-су модельдері бұл механизмді мүлдем қамтымайды.

2. Экологиялық негіз

2.1 Түрлер

ТүрҮлесіЭкологиясыМаңызы
Culex modestus56,2 %Тұрақты су: қамыс көлшіктері, суару каналдары, қала дренажыБатыс Нил вирусының (WNV) басты тасымалдаушысы
Aedes caspius (және vexans, flavescens, cinereus)басым топТасқын-су, тұздылау батпақ. Десикацияға төзімді диапаузалы жұмыртқаМассалық мазалау — тұрғындарға тікелей әсері
Anopheles (maculipennis, hyrcanus)1,8 %Тұрақты су, жылы климатАз, бірақ жылынумен таралуда

Батыс Қазақстан / Жайық-Ақжайық аулау жинақтары бойынша. Атырау қаласының өз тұзақ базасы жоқ.

Модельге тікелей әсер еткен түзету. Әдетте «Атырауда Aedes caspius басым» деп болжанады. Ал аулау деректерінде Culex modestus доминант (56 %). Сондықтан модель екі экологияны БӨЛЕК ұстайды: тасқын-су Aedes (импульсті шыңдар) және тұрақты-су Culex (WNV, жаз ортасы), айлық салмақпен ауысады.

2.2 Басты драйверлер (дәлелденген реті бойынша)

  1. Көктемгі тасқынның көлемі мен уақыты — жұмыртқа су басумен жарылады
  2. Гидропериод — көлшік дернәсілдің дамуын аяқтауға жеткілікті ұзақ тұруы қажет
  3. Қамыс алқаптары (Phragmites) — Culex modestus үшін ең күшті мекен предикторы
  4. Температура — буын санын және даму жылдамдығын белгілейді
  5. Топырақ ылғалы, жауын — қосалқы су көздері
  6. Қалалық амплификация — дренаж, подвал, суару арықтары

3. Модельдің архитектурасы

Жеті қабат. Әрқайсысының іске асу күйі ашық көрсетілген — жоба құжатында сипатталғанның бәрі кодта бар емес.

ҚабатМазмұныКүйіІске асуы
L1Гидро-мекен: су басуды анықтауістейдіSentinel-1 SAR + GloFAS → тасқын импульсі
L1Каспий деңгейі модуляторыжоқТексерілген тегін дереккөз табылмады
L2FPEB ядросы: жұмыртқа банкі динамикасыістейдіdE/dt, dL/dt, dA/dt тәуліктік интеграция
L2ГидропериодістейдіS1 өтулері бойынша «кемінде N күн»
L3Температура гейтіістейдіMordecai 2019 термиялық шектері
L4Culex / WNV тармағыістейдіSentinel-2 NDVI → қамыс мекені
L5ML болжам (LSTM/GRU + RF/SHAP)жоқВалидациясыз қосу жалған дәлдік берер еді
L6Bayesian белгісіздік интервалыжоқСол себеп
L7Digital twin + ассимиляцияішінараКүй қайта интегралданады; ассимиляциялайтын бақылау жоқ

4. Формулалар

4.1 Жиынтық индекс

MRI = 100 · Φ_T(T) · (0,15 + 0,85 · S) · (1 + 0,4 · U + 0,5 · W)

  Φ_T — температура гейті (0…1)
  S   — түрлер қоспасы: S = clamp(w_A(ай) · A + w_C(ай) · C)
  A   — тасқын-су Aedes ересек индексі (FPEB интеграциясынан)
  C   — тұрақты-су Culex индексі
  U   — қалалық амплификация (елді мекенге жақындық, 0…1)
  W   — тасқын мүшесі = бейімділік × өлшенген импульс

Тұрақты мүше 0,15 — температура рұқсат еткен кездегі фондық белсенділік: су режимі нашар болса да маса мүлдем нөл болмайды.

4.2 Температура гейті Φ_T

Φ_T(T) = (T − T₀)(Tm − T) / ((Tm − T₀)/2)²,   T₀ < T < Tm
Φ_T(T) = 0                                    әйтпесе

  T₀ = 16,8 °C     Tm = 34,9 °C     шыңы = 25,9 °C

Шектер Mordecai т.б. (2019) Батыс Нил вирусы жүйесі үшін жариялаған мәндерден. Нормалау шыңда дәл 1 береді. Мақаладағы Brière фитінің дәл көшірмесі емес — симметриялы жуықтау, бірақ шыңы (25,9 °C) мақалада хабарланған ~25 °C оптимумына сәйкес келеді.

4.3 FPEB ядросы — тасқын-импульсті жұмыртқа банкі

Модельдің жүрегі. Тәуліктік қадаммен интегралданатын үш теңдеу:

dE/dt = f · (φ_egg(ай) − E)  −  h(W) · E
dL/dt = h(W) · E  −  L/τ(T)  −  μ_L · L
dA/dt = (L/τ(T)) · s(hydro)  −  μ_A · A

  E — жұмыртқа банкі (диапаузалы, құрғақта тірі қалады)
  L — дернәсіл
  A — ересек
  h(W) = k_h · W            жарылу жылдамдығы (су басу қосады)
  τ(T) = DD / (T − T_base)  даму ұзақтығы, тәулік
  s(hydro)                  гидропериодтан шыққан тірі қалу
  φ_egg(ай)                 маусымдық банк сыйымдылығы
Осы теңдеулердің мәні — КІДІРІС. Су басқан күні маса пайда болмайды: алдымен жұмыртқа жарылады, содан кейін дернәсіл τ(T) күн дамиды. 25 °C-та бұл ≈ 10 тәулік. Кідіріссіз модель шыңды дұрыс емес күнге қояды әрі «массалық шығу қашан күтіледі» деген ең пайдалы сұраққа жауап бере алмайды.

4.4 Параметрлер

ПараметрМәніМағынасыНегізі
DD150 °C·тәулікЖұмыртқа→ересек градус-күнКулициндік әдебиет
T_base10 °CДаму табалдырығыКулициндік әдебиет
μ_L0,10 /тәулікДернәсіл өліміӘдебиеттегі типтік шама
μ_A0,12 /тәулікЕресек өліміӘдебиеттегі типтік шама
k_h0,6Жарылу коэффициентіМодель баптауы
f0,15 /тәулікБанктің толығуыМодель баптауы
T₀ / Tm16,8 / 34,9 °CТермиялық шектерMordecai 2019 ✓
Жеті параметрдің тек біреуі (термиялық шектер) нақты жарияланған зерттеуден алынған. Қалғандары — кулициндік әдебиеттегі жалпы шамалар мен модель баптаулары. Aedes caspius үшін жергілікті калибрлеу ЖОҚ. Бұл — модельдің басты әлсіздігі.

4.5 Тасқын мүшесі W

W = БЕЙІМДІЛІК × ИМПУЛЬС

  БЕЙІМДІЛІК — Жайық арнасына/атырауға жақындық (география, тұрақты)
  ИМПУЛЬС    — 0,65 · SAR + 0,35 · GloFAS  (өлшенген, күн сайын өзгереді)
  SAR        = min(1, артық су ауданы % / 5%)
  GloFAS     = ағынның өз терезесіндегі қатынасы

Бұл бөліну маңызды: бейімділік «қай жер су басуға бейім» дегенді, импульс «бүгін су басты ма» дегенді көрсетеді. Бұрын екеуі араласып, сәуірде де, қаңтарда да бірдей мән беретін.

4.6 Culex тармағы

C = clamp(0,55 · R + 0,30 · W + 0,15 · U) · H(rh)

  R     — қамыс мекені (Sentinel-2 NDVI > 0,4 үлесі / 0,3)
  H(rh) — ылғалдылық көбейткіші: (rh − 40) / 45

5. Кіріс деректері

ДерекДереккөзТүріКідіріс
Су басқан аумақ, гидропериодCopernicus Sentinel-1 GRD (VV, GAMMA0)🛰 өлшем1–6 тәулік
Қамыс мекеніCopernicus Sentinel-2 L2A (NDVI)🛰 өлшем1–5 тәулік
Өзен ағыныGloFAS (Open-Meteo Flood API)📊 модель1 тәулік
Температура, ылғал, жауынECMWF (Open-Meteo)📊 модель1 сағат
Топырақ ылғалыECMWF топырақ моделі📊 модель1 сағат
Қала/елді мекен тізіліміЖобаның өз тізілімі (8 елді мекен)тізілім
Есептеу нүктелеріАтыраудың 65 нүктесі + облыстық 5×5 тортізілім

Есептеу тізбегі

  1. 90 нүкте бойынша ECMWF метеорологиясы алынады (өткен 30 + алдағы 14 тәулік)
  2. GloFAS күндік ағыны → күндік тасқын импульсі W(t)
  3. Sentinel-1 бақылау терезелерінде артық су ауданы мен гидропериод өлшенеді
  4. Sentinel-2 NDVI бойынша қамыс мекені бағаланады
  5. Әр нүкте үшін FPEB теңдеулері 30 тәуліктік «жүгіріспен» интегралданады
  6. Ересек индексі A(t) идеал жағдайдағы теориялық шыңға нормаланады
  7. Culex тармағы қосылып, айлық салмақпен араластырылады
  8. Температура гейті мен амплификация қолданылып, 0–100 шкаласына келтіріледі
Күй дерекқорда сақталмайды. Әр сұраныста нөлден қайта интегралданады. Артықшылығы: нәтиже детерминистік, қайталанады, күй жылжып кетпейді. Жобадағы өрт қаупі индексі (FWI) де дәл осылай есептеледі — бұл салада қалыптасқан тәсіл.

6. Әлемдік модельдермен салыстыру

МодельТәсіліКүшті жағыАтырауға жарамауының себебі
MoLSМеханистік өмір цикліТолық биологиялық тізбекПараметрлері ыдыс-су Ae. aegypti үшін. Тасқын-импульс механизмі жоқ
Aedes-AIНейрожелі (LSTM/GRU)Жылдам, ықтималдық шығысMoLS-тің aegypti биологиясын мұрагерлейді, АҚШ-қа калибрленген
MAMOTH / EYWAСпутник + метео + тұзақОперациялық, Еуропада қолданыста. Ең жақын аналогОқыту үшін тұзақ деректері МІНДЕТТІ — Атырауда олар жоқ
VECTRIКлиматтық-динамикалықКонтиненттік ауқым, су айдынын қоса есептейдіБезгек (Anopheles) үшін, Африкаға бағдарланған
Mordecai термиялық модельдеріФизиологиялық-термиялықЖергілікті дерексіз ауысады, рецензияланғанТЕК температура. Су, мекен, гидрология жоқ
Камарг Ae. caspius моделіТасқынға негізделген популяция динамикасыДәл біздің түрге арналған, тасқын-импульс механизмі барСпутникпен байланыспаған, операциялық жүйе емес
RF / XGBoost + SHAPСтатистикалық оқытуТүсіндірмелі, бейсызықОқыту үшін бақылау керек — жоқ. Экстраполяциясы нашар

Салыстырмалы қасиеттер

ҚасиетJAIYQ-MRIMAMOTHMoLSAedes-AIКамаргMordecai
Тасқын-импульс ажыратқышыиәжоқжоқжоқиәжоқ
Диапаузалы жұмыртқа банкіиәжоқішінаражоқиәжоқ
Спутниктен өлшенген суиәиәжоқжоқжоқжоқ
Гидропериод (ұзақтығы)иәжоқжоқжоқжанамажоқ
Қамыс мекені (NDVI)иәиәжоқжоқжоқжоқ
Қос түр, маусымдық салмақиәішінаражоқжоқжоқжоқ
Тұзақ деректерінсіз жұмысиәжоқиәжоқжоқиә
ВАЛИДАЦИЯЛАНҒАНЖОҚиәиәиәиәиә
Белгісіздік интервалыжоқішінаражоқиәжоқиә

7. Уникальділік — шыншыл баға

Не жаңа ЕМЕС

Модельдің бірде-бір компоненті жаңа емес. Барлығы жарияланған, белгілі әдістер:

Не жаңа

Жаңалық — компоненттерде емес, олардың ЖИНАҚТАЛУЫНДА және ОПЕРАЦИЯЛЫҚ ІСКЕ АСУЫНДА:

  1. Спутниктен ӨЛШЕНГЕН тасқын импульсі мен гидропериод механистік жұмыртқа банкі теңдеулерінің тікелей драйвері ретінде. Камарг моделі тасқынды есептейді, бірақ спутникпен байланыспаған. MAMOTH спутникті қолданады, бірақ жұмыртқа банкі динамикасы жоқ әрі тұзақпен оқытылады. Екеуін біріктірген жүйе әдебиетте табылмады.
  2. Қос түрдің фазаға тәуелді салмағы — көктемде тасқын-су Aedes, жаз ортасында WNV-тасымалдаушы Culex — бір архитектурада.
  3. Тұзақ деректерінсіз операциялық жұмыс. Ең жақын аналог (MAMOTH) тұзақсыз іске қосыла алмайды.
  4. Солтүстік Каспий атырауына арнайы бейімделу — бұл аймақ үшін жарияланған операциялық маса моделі табылмады.
Бірақ: «үздік» немесе «дәлірек» деген тұжырым ЖАСАЛМАЙДЫ. Модель валидацияланбаған. Салыстыру кестесіндегі «иә» белгілері мүмкіндіктің бар-жоғын көрсетеді, дәлдікті емес. Валидацияланған модельмен (MAMOTH, MoLS) дәлдік бойынша салыстыру үшін алдымен өз тұзақ базамыз болуы керек. Ғылыми жарияланымда бұл модельді «жаңа архитектура ұсынылды» деп сипаттауға болады, «нақтырақ болжайды» деп сипаттауға болмайды.

8. Не жетіспейді және оны шешу жолы

8.1 Валидация — ең үлкен олқылық

Мәселе: Модель ешқашан нақты маса санымен салыстырылмаған. «MRI 70 нүктесінде MRI 40 нүктесіне қарағанда шынымен көп маса бар ма?» деген сұраққа жауап жоқ.

Шешу жолы:

  • 6–10 тұзақ (CDC light trap немесе BG-Sentinel) индексі әртүрлі нүктелерге: 2 — Жайық жағасы, 2 — қамыс алқабы, 2 — қала ішінде, 2 — құрғақ бақылау нүктесі
  • Кезеңі: мамыр–қыркүйек, аптасына 2 рет санау
  • Ұзақтығы: кемінде 1 маусым (шешуші), 2 маусым (сенімді)
  • Нәтижесі: MRI мен тұзақ саны арасындағы Спирман корреляциясы. ρ > 0,6 болса модель «валидацияланған» деп белгіленеді
  • Кім істей алады: облыстық СЭС, ветеринария қызметі немесе университет биология кафедрасы — жабдық қымбат емес

8.2 Түрге тән параметрлер

Мәселе: Жеті параметрдің алтауы Aedes caspius үшін калибрленбеген.

Шешу жолы:

  • Тұзақ деректері жиналған соң параметрлерді кері фиттеу (least-squares немесе Bayesian)
  • Немесе зертханада дернәсілдің даму ұзақтығын әртүрлі температурада өлшеу — бір маусымдық жұмыс

8.3 Каспий деңгейінің әсері

Мәселе: Каспий деңгейі 1996 жылдан бері түсуде. Ол тұздылау батпақтардың ауданын қайта пішіндейді — бірақ модельде ескерілмеген.

Шешу жолы:

  • Спутниктік альтиметрия деректері (Hydroweb, G-REALM, Copernicus Marine)
  • Ай сайынғы деңгей қатары модельге баяу модулятор ретінде қосылады
  • Кедергі: тексерілген тегін API әлі табылмады — қолмен жүктелетін файлдарды автоматтандыру керек

8.4 Гидропериодтың дәлдігі

Мәселе: Sentinel-1 қайталауы ~6 тәулік, сондықтан «су кемінде N күн тұрды» деген төменгі шек қана шығады. Екі өту арасында су кеуіп, қайта басуы мүмкін.

Шешу жолы:

  • Sentinel-2 су индексін (MNDWI) қосу — бұлтсыз күндері қосымша өтулер береді, уақыттық тығыздық шамамен екі есе артады
  • Инфрақұрылым дайын: S2 сұранысы қамыс мекені үшін жазылған

8.5 Қамыс классификациясының тексерілуі

Мәселе: NDVI > 0,4 — тығыз өсімдік, бірақ суармалы егіс пен ағаш екпелері де сол шекке кіреді.

Шешу жолы:

  • Жоғары ажыратымдылықты сурет бойынша 50–100 нүктені қолмен тексеру
  • Sentinel-2 көп арналы классификациясы (қамыстың су үстіндегі спектрі ерекше)
  • Далалық маршрутпен растау

8.6 Белгісіздік интервалы

Мәселе: Модель нүктелік сан береді, ал шын мәнінде белгісіздігі үлкен.

Шешу жолы:

  • Параметрлердің ықтимал диапазонымен ансамбль жүргізу (50–100 жүгіріс) → пайыздық интервал
  • Бұл ЖАҢА дерек талап етпейді, тек есептеу
  • Бірақ интервалдың ені калибрлеусіз шындыққа сай болмайды — валидациядан КЕЙІН істеген дұрыс

8.7 Басымдық реті

ЖұмысКүшТәуелділігі
1Тұзақ желісін орнату (валидация)орташаДалалық жұмыс, серіктес керек
2Гидропериодқа S2 қосуазТәуелділігі жоқ — бірден істеуге болады
3Қамыс классификациясын тексеруазТәуелділігі жоқ
4Каспий деңгейі модуляторыорташаДереккөзді автоматтандыру
5Параметрлерді фиттеуорташа№1 аяқталуы керек
6Ансамбль + белгісіздікорташа№5 аяқталуы керек
7ML қабаты (L5)үлкен№1 — оқыту үшін жапсырма керек

9. Шектеулер — қысқаша тізім

10. Дереккөздер

  1. Mordecai E.A. et al. (2019) — Thermal biology of mosquito-borne disease, PLoS NTD
  2. Mordecai E.A. et al. (2017) — Detecting the impact of temperature on transmission, PLoS NTD
  3. Aedes caspius flooding population-dynamics model (Камарг), Bulletin of Entomological Research
  4. Батыс Қазақстан WNV / Culex modestus, Frontiers in Public Health (2020)
  5. WNV Қазақстан жылқыларында, Атырауды қоса (8,7 %)
  6. MAMOTH — Earth Observation data-driven model for mosquito abundance; EYWA жүйесі
  7. MoLS — Mosquito Life Cycle Simulation
  8. Aedes-AI нейрожелілері, PLoS Comput Biol
  9. Van Wagner C.E. (1987) — Canadian FWI жүйесі (spin-up тәсілінің үлгісі)
  10. Ақжайық Рамсар алаңы
Jaiyq · ecojaiyq.com · Қазақстан мен Каспий жағалауының экологиялық AI мониторинг платформасы
Бұл құжат жүйенің ағымдағы техникалық күйін сипаттайды. Модель валидацияланбаған: шығысы тексеру тағайындауға негіз болады, бірақ сот немесе әкімшілік іс үшін дәлел ретінде қолданылмайды.