Ile wynosi ryzyko DCS? Jak algorytm przewiduje chorobę dekompresyjną z profilu nurkowania

Zrozumienie ryzyka to nie to samo co strach. Przeciwnie — kiedy wiesz, skąd bierze się liczba i co nią rusza, zamiast lęku masz spokojną kontrolę. Ten artykuł jest dla osób, które chcą zajrzeć pod maskę: jak to możliwe, że z kilku liczb opisujących profil — głębokość, czas, gazy, tempo wynurzania — da się oszacować prawdopodobieństwo choroby dekompresyjnej (DCS). I dlaczego dwie pozornie podobne decyzje na wynurzeniu dają zupełnie różne ryzyko.

Dwa pytania, dwa rodzaje modeli

Klasyczny komputer nurkowy i tabele odpowiadają na pytanie zero-jedynkowe: „czy z tej głębokości wolno mi już wynurzyć się wyżej?". To model planujący. W rekreacji i technice króluje tu Bühlmann ZH-L16C z Gradient Factors: dzieli ciało na 16 „przedziałów tkankowych" o różnym tempie nasycania azotem i pilnuje, by w żadnym z nich napięcie gazu nie przekroczyło dopuszczalnej granicy (M-value). Wynik to sufit i harmonogram przystanków — ale wciąż odpowiedź „tak/nie".

Model predykcyjny odpowiada na inne pytanie: „jak bardzo prawdopodobne jest, że TEN konkretny profil skończy się DCS?". Bo DCS jest zdarzeniem losowym — ten sam profil u tej samej osoby raz kończy się dobrze, raz nie. Zamiast twardej granicy dostajemy prawdopodobieństwo: P(DCS). I właśnie taki model — probabilistyczny — jest bohaterem tego tekstu.

Jak algorytm liczy prawdopodobieństwo

Pomysł, który spina całą tę rodzinę modeli, sformułowali Weathersby, Homer i Flynn w pracy „On the likelihood of decompression sickness" (Journal of Applied Physiology, 1984). Brzmi tak: w każdej chwili nurkowania istnieje pewne chwilowe natężenie ryzyka (hazard), które rośnie, gdy tkanka jest przesycona gazem obojętnym względem ciśnienia otoczenia. Całkowite prawdopodobieństwo liczy się jak w analizie przeżycia — z sumy (całki) tego hazardu po całym nurkowaniu:

Tu G to skumulowany hazard, g_i to wagi poszczególnych przedziałów ryzyka, a r_i(t) to chwilowe ryzyko w przedziale i. To ostatnie rośnie z przesyceniem — z nadwyżką napięcia gazu w tkance ponad ciśnienie otoczenia, odniesioną do tego ciśnienia:

Ryzyko nie kumuluje się na dnie. Ono narasta przede wszystkim na wynurzeniu i jeszcze przez wiele godzin po wyjściu z wody — gdy ciśnienie spada, a tkanki są wciąż nasycone.

To kluczowa intuicja, którą widać wprost we wzorze: gdy P_i ≤ P_amb (jesteś głęboko, tkanka dopiero się nasyca), ryzyko jest zerowe. Rośnie dopiero, gdy ciśnienie otoczenia spada poniżej napięcia w tkance — czyli na wynurzaniu i na powierzchni. Dlatego całkę prowadzi się przez standardowe okno 24 godzin po nurkowaniu (większość objawów DCS ujawnia się w tym czasie; to okno jest zarazem parametrem kalibracji modelu). Konkretny algorytm, którego tu użyłem, to LE1nt — wariant modelu liniowo-wykładniczego (linear-exponential) z linii Thalmanna i amerykańskiej marynarki. Trzy „przedziały ryzyka" (osobne od tkanek Bühlmanna, z własnymi stałymi czasowymi) sumują wkłady do G, a wynik mapuje się na prawdopodobieństwo w wersji trójmianowej: brak DCS / postać łagodna / postać poważna. Pokazany wyżej wzór to historyczny, binomialny fundament (Weathersby, 1984); konkretne liczby w tym artykule pochodzą z jego trójmianowego uściślenia (model mnoży skumulowany hazard przez czynnik a+1, gdzie a≈0,116), co podnosi je o kilkanaście procent względem samego 1−e^(−G).

Jest tu jeszcze jeden niuans, który warto znać: w fazie liniowej modelu ryzyko przedziału pośredniego liczy się nie z samego napięcia gazu, lecz z „obciążenia", które po przekroczeniu pewnego progu wymywa się liniowo — wolniej niż wykładniczo. To właśnie dlatego tkanki pośrednie pozostają przesycone i groźne jeszcze długo po wyjściu z wody.

Skąd wiadomo, że liczby mają sens? Współczynniki dopasowano metodą największej wiarygodności do zbioru BIG292 — 3322 rzeczywistych nurkowań eksperymentalnych marynarek USA, Wielkiej Brytanii i Kanady, w których odnotowano 190 przypadków DCS. To realne dane o ludziach, nie teoria. Współczynniki LE1nt pochodzą z otwartej pracy Howle i in. (PLoS ONE, 2017), a dokładna kinetyka — z Howle, Weber i Vann (Computers in Biology and Medicine, 2009).

Ta sama rodzina opisuje też rolę tlenu (Parker, Survanshi, Weathersby, 1998) i nowsze, dokładniejsze modele sprzężonej perfuzji i dyfuzji (dwie części „Probabilistic pharmacokinetic models of DCS", 2017 i 2018). Wszystkie odpowiadają na to samo pytanie — z różną wiernością fizjologii:

• Weathersby 1984 — fundament: DCS jako zdarzenie losowe, ryzyko z całki hazardu, dopasowanie największej wiarygodności.
• Linia LE / LEM (Thalmann, NMRI98) — ryzyko z napięcia gazów obojętnych; LEM dokłada osobny wpływ tlenu.
• LE1nt (Howle 2017) — wersja trójmianowa: rozróżnia DCS łagodne i poważne; to jej używa symulator poniżej.
• Modele PK (2017/2018) — sprzężona perfuzja-dyfuzja, najbliżej fizyki pęcherzyków; bardziej złożone, ta sama logika.

Plan to nie ryzyko — i tu wchodzą Gradient Factors

Zanim policzymy ryzyko, potrzebujemy planu — harmonogramu przystanków. Tu sterują Gradient Factors (GF): dwa pokrętła konserwatyzmu nałożone na Bühlmanna. GF Low (np. 30%) mówi, jak blisko granicy M-value wolno Ci podejść na pierwszym, najgłębszym przystanku — niższa wartość wymusza wcześniejsze, głębsze przystanki. GF High (np. 70%) mówi, jak nasycony możesz być w chwili wyjścia na powierzchnię — niższa wartość wydłuża płytkie przystanki. „GF 30/70" to właśnie ta para. GF 100/100 to czysty Bühlmann bez marginesu; im niżej, tym ostrożniej i tym dłuższa dekompresja. GF zmienia plan, a model probabilistyczny przelicza każdy plan na jedną liczbę — P(DCS).

Przykład: 40 m przez 40 minut, powietrze + EAN50 + tlen

Weźmy konkretne, wymagające nurkowanie: 40 metrów, 40 minut na dnie na powietrzu, do dekompresji butla EAN50 (przepinana ok. 21 m) i czysty tlen na płytkich przystankach. Plan liczymy na bazie ZHL-16C + GF, ryzyko modelem LE1nt. Oto, jak ten sam profil wygląda przy różnych ustawieniach GF (czas dekompresji obejmuje przystanki i przepięcia gazu):

• GF 100/100 (bez marginesu): ~25 min dekompresji, P(DCS) ≈ 0,99%.
• GF 70/85: ~35 min, P(DCS) ≈ 0,78%.
• GF 50/85: ~37 min, P(DCS) ≈ 0,75%.
• GF 40/85 (typowe „techniczne"): ~38 min, P(DCS) ≈ 0,74%.
• GF 30/70 (ostrożne): ~47 min, P(DCS) ≈ 0,58%.
• GF 20/60 (bardzo ostrożne): ~56 min, P(DCS) ≈ 0,47%.

Widać prawidłowość: każda dodatkowa minuta dekompresji obniża ryzyko, ale z malejącym zyskiem. Przejście z GF 100/100 na 40/85 wydłuża wynurzanie o połowę (×1,5) i ścina ryzyko o jedną czwartą. Dalsze dokręcanie do 20/60 dokłada kolejne kilkanaście minut (i sporo ekspozycji tlenowej — CNS), a urywa już tylko ułamki procenta. Dla kontekstu: spokojne nurkowanie rekreacyjne (18 m przez 50 min na powietrzu, GF 85/85) ma w tym modelu P(DCS) rzędu 0,2% — kilkukrotnie mniej niż nasza techniczna „czterdziestka".

Najlepiej poczuć to samemu. Poniższy symulator liczy plan i ryzyko na żywo — tym samym silnikiem co liczby w tym artykule (Bühlmann ZH-L16C + GF oraz model LE1nt), zweryfikowanym wprost względem pełnego symulatora nurkowego. Pokręć głębokością, czasem i suwakami GF, włączaj i wyłączaj gazy deko, i patrz, co dzieje się z prawdopodobieństwem:

Co naprawdę dzieje się, gdy skracasz dekompresję

Tu zaczyna się najciekawsza część. Wróćmy do planu GF 30/70 (47 min dekompresji, P(DCS) ≈ 0,58%) i zróbmy trzy „grzechy" wynurzania, każdy policzony tym samym modelem. Najpierw scenariusz skrajny: wynurzenie wprost na powierzchnię, bez żadnych przystanków. Potem dwa subtelniejsze: usunięcie dokładnie ⅓ czasu dekompresji (ok. 16 minut) — raz z przystanków głębokich, raz z płytkich.

• Pełna dekompresja (47 min): P(DCS) ≈ 0,58% — punkt odniesienia.
• Pominięte ⅓ czasu z przystanków GŁĘBOKICH: P(DCS) ≈ 0,62% — ryzyko prawie się nie zmienia (×1,06).
• Pominięte ⅓ czasu z przystanków PŁYTKICH: P(DCS) ≈ 1,20% — ryzyko się podwaja (×2,1).
• Wynurzenie wprost, bez przystanków: P(DCS) ≈ 1,99% — ponad trzy razy więcej (×3,4).

Wykres pokazuje to samo od strony czasu. Przy pełnej dekompresji (linia spokojna) ryzyko pozostaje znikome nawet przez całe wynurzanie — przystanki trzymają tkanki poniżej progu przesycenia — i kończy łagodnie przy ~0,58%. Wynurzenie wprost (linia czerwona) startuje wcześniej i pnie się do ~1,99%, a gros tego przyrostu dzieje się dopiero PO wyjściu z wody. To nie dno jest groźne, lecz powrót.

Wróćmy do liczb: usunęliśmy dokładnie tyle samo minut, a efekt jest diametralnie różny w zależności od tego, KTÓRE minuty wycięliśmy. Dlaczego? Bo to płytkie przystanki domykają odsycanie przedziałów wiodących w momencie wyjścia. Ostatni skok ciśnienia — z 3–6 metrów do zera — jest największy względnie i najgroźniejszy; jeśli wykonasz go z wciąż nasyconymi wolnymi tkankami, hazard strzela w górę. Przystanki głębokie pilnują głównie szybkich tkanek, które i tak nie są limitujące przy wyjściu — ich skrócenie niemal nic nie zmienia.

To nie jest tylko artefakt modelu. Eksperyment amerykańskiej marynarki (NEDU, 2011) porównał wprost dwa harmonogramy o tym samym czasie: jeden z naciskiem na przystanki głębokie, drugi na płytkie. Wersja z głębokimi przystankami dała WIĘCEJ przypadków DCS — bo trzymała wolne tkanki nasycone dłużej. Stąd współczesny zwrot ku wyższym GF Low. Nasz symulator mówi dokładnie to samo: minuta spędzona płytko jest dla bezpieczeństwa warta więcej niż minuta spędzona głęboko.

Jeśli kiedykolwiek miałbyś skrócić dekompresję — a planem jest jej NIE skracać — to ostatnie metry są tymi, których nie wolno oddać.

Czego ten model nie wie

Uczciwość jest częścią bezpieczeństwa, więc jasno: to model dydaktyczny i ma granice. Liczby ilustrują mechanizm, nie są planem nurkowania ani poradą medyczną.

• Domena kalibracji to powietrze i nitroks na poziomie morza. Hel (trimiks), nurkowanie wysokogórskie i obiegi zamknięte (CCR) to już ekstrapolacja poza dane.
• Model nie zna Ciebie: nawodnienia, snu, wysiłku, temperatury, przetrwałego otworu owalnego (PFO), formy dnia. Te czynniki potrafią przesunąć ryzyko bardziej niż kilka minut przystanku.
• Bezwzględna wartość procentu jest niepewna co do czynnika rzędu 2–3× i zależy od zbioru kalibracyjnego — akurat ten wariant (LE1nt na zbiorze BIG292) odczytuje ryzyko NIŻEJ niż modele kalibrowane wprost na danych marynarki (NMRI98), więc nie traktuj tych procentów jako pesymistycznej, „bezpiecznej" górnej granicy. Wiarygodne i pouczające są PORÓWNANIA: o ile ryzyko rośnie, gdy zmieniasz profil.
• Żaden model nie zastępuje szkolenia, własnego komputera i tabel oraz zdrowego rozsądku w wodzie.

Co z tego wynika dla Ciebie

Nie po to liczymy ryzyko, żeby się bać, tylko żeby podejmować lepsze decyzje na spokojnie. Trzy rzeczy zostają z tej lekcji: dekompresja działa i opłaca się robić ją w całości; konserwatyzm (niższe GF) kupuje bezpieczeństwo, ale z malejącym zyskiem i kosztem tlenowym; a jeśli cokolwiek miałoby pójść nie tak — to płytkie, „nudne" przystanki są tymi, których pilnujesz najbardziej.

Tej spokojnej, świadomej relacji z liczbami uczę na kursach — od pierwszego oddechu po planowanie nurkowań, na których takie wykresy przestają być abstrakcją. Bo chodzi o to, żeby rozumieć, a nie tylko „wciskać OK na komputerze" — reszta to już Twoja decyzja, podjęta z głową.