Bhangmetry i podwójny błysk: optyczne wykrywanie detonacji jądrowych

Streszczenie

Najbardziej znany obraz wybuchu jądrowego to "oślepiający błysk", ale dla aparatury pomiarowej ważniejsza jest nie sama jasność, tylko przebieg jasności w czasie. Atmosferyczna kula ognista może dać charakterystyczny podwójny impuls: bardzo szybki pierwszy pik, chwilowe osłabienie i wolniejszy drugi wzrost związany z przezroczystością fali uderzeniowej.^1,2

Bhangmeter jest szybkim, nieobrazującym radiometrem zaprojektowanym właśnie do rozpoznawania takiego przebiegu optycznego. Ten artykuł wyjaśnia fizykę podwójnego błysku, sens toru fotometrycznego, różnicę między amplitudą i kształtem sygnału oraz ograniczenia takiej detekcji na przykładzie programu Vela i incydentu z 22 września 1979 roku.³

Rozszerzenie tematu

Błysk jako przebieg, nie fotografia

Oko widzi wybuch jako nagłe rozjaśnienie. Kamera widzi kolejne klatki. Bhangmeter widzi coś jeszcze prostszego i zarazem bardziej użytecznego: jedną liczbę jasności jako funkcję czasu.

To jest zasadnicza różnica między obrazowaniem a fotometrią. Bhangmeter nie musi tworzyć zdjęcia kuli ognistej ani rozpoznawać kształtu chmury. Ma mierzyć szybkie zmiany natężenia światła w wybranym polu widzenia i sprawdzić, czy przebieg przypomina znaną sygnaturę atmosferycznej detonacji jądrowej.

W praktyce oznacza to, że najważniejszy jest czasowy profil sygnału. Dwa zdarzenia o podobnej maksymalnej jasności mogą być różne: błysk pioruna, odbicie słoneczne, meteor i kula ognista mogą dawać różne czasy narastania, zaniku, przerwy i drugiego maksimum.

Skąd bierze się podwójny błysk

Podwójny błysk jest skutkiem sprzężenia transportu promieniowania i hydrodynamiki fali uderzeniowej. W bardzo wczesnej fazie energia wybuchu jest przenoszona głównie przez promieniowanie rentgenowskie, które ogrzewa otaczające powietrze i tworzy sferę izotermiczną. Na tym etapie temperatura jest ogromna, ale powietrze i produkty jonizacji nie są prostym, przezroczystym oknem dla obserwatora.^1,2

Gdy kula ognista rozszerza się i chłodzi, transport radiacyjny przestaje dominować nad hydrodynamiką. Powstaje fala uderzeniowa, której rozgrzana, zjonizowana powłoka świeci, ale zarazem zasłania gorętsze wnętrze. Obserwowana jasność spada, mimo że wewnątrz nadal znajduje się bardzo gorąca materia.¹

Później front fali chłodzi się i staje się bardziej przezroczysty. Wtedy gorętsza, większa sfera wewnętrzna znów staje się widoczna, tworząc drugi impuls optyczny. Właśnie ten układ: szybki pierwszy błysk, minimum i wolniejszy drugi wzrost, był podstawą optycznego rozpoznawania atmosferycznych detonacji jądrowych.^1,2

Breakaway

W źródłach opisujących skutki wybuchów jądrowych kluczowy jest moment breakaway. To chwila, gdy fala uderzeniowa po wcześniejszym zasłanianiu wnętrza kuli przestaje być optycznie dominującą, nieprzezroczystą powłoką. Wtedy obserwator zaczyna widzieć znowu gorące wnętrze kuli ognistej.¹

W artykule źródłowym NW FAQ dla przykładowego wybuchu 20 kt podano skalę: minimum jasności pojawia się po rzędu 10 ms, a drugi szczyt po rzędu 150 ms. Dla 1 Mt podane czasy są dłuższe, rzędu 100 ms dla minimum i 900 ms dla drugiego maksimum.¹ Nie są to tu instrukcje estymacji mocy ładunku, tylko liczby porządkujące zjawisko fizyczne.

Najważniejszy wniosek dydaktyczny jest taki: pierwszy pik jest ekstremalnie szybki i jasny, ale niesie małą część całkowitej energii cieplnej. Drugi impuls jest wolniejszy i odpowiada za większość promieniowania cieplnego, bo trwa znacznie dłużej.¹

Co mierzy bhangmeter

Bhangmeter można rozumieć jako szybki fotometr o dużym zakresie dynamicznym. W najprostszym schemacie ma:

• element światłoczuły, na przykład fotodiodę,
• optykę lub pole widzenia ograniczające kierunki zbierania światła,
• tor zamieniający prąd fotoelektryczny na napięcie,
• szybką rejestrację przebiegu,
• algorytm lub procedurę porównującą kształt z oczekiwanym wzorcem.

Ważne jest słowo "nieobrazujący". Taki przyrząd nie rozpoznaje obiektów jak kamera. Może nie wiedzieć, czy źródło było okrągłe, wydłużone, częściowo zasłonięte albo poza osią. Dostaje tylko krzywą jasności w czasie.

To jest zaleta i wada. Zaleta, bo prosty radiometr może być szybki, odporny i łatwiejszy do kalibracji niż kamera. Wada, bo utrata obrazu oznacza większą zależność od filtracji fałszywych alarmów i od korelacji z innymi sensorami.

Amplituda nie wystarcza

Gdyby liczyła się tylko amplituda, każdy bardzo jasny błysk byłby podejrzany. To nie działa. Atmosfera, powierzchnie chmur, ocean, pioruny, meteory i odbicia słoneczne potrafią wygenerować bardzo silne sygnały optyczne.

Dlatego bhangmeter pracuje na kształcie sygnału:

• jak szybko pojawia się pierwszy pik,
• czy występuje minimum,
• jak narasta drugi impuls,
• jak długo trwa ogon,
• czy proporcje obu części są fizycznie spójne,
• czy drugi kanał lub drugi przyrząd widzi podobny przebieg.

To jest klasyczne zadanie detekcji sygnału w tle. Sama liczba maksymalna jest zbyt uboga. Potrzebny jest wektor cech czasowych.

Dlaczego pierwszy pik jest trudny

Pierwszy pik jest krótki. Wymaga dużej szybkości toru, małej bezwładności elementu światłoczułego i rejestracji, która nie zgubi milisekundowych szczegółów. W elektronice oznacza to szerokie pasmo, małe pojemności pasożytnicze, dobrą liniowość albo kontrolowane nasycenie oraz stabilne znaczniki czasu.

W artykule o oscyloskopie w pomiarach impulsów nanosekundowych ten sam problem pojawia się w mniejszej skali: jeśli tor jest za wolny, impuls staje się sztucznie szerszy i niższy. Bhangmeter ma identyczną logikę, tylko zamiast impulsu elektrycznego w kablu mierzy impuls optyczny przetworzony przez fotodiodę.

Nie znaczy to, że trzeba publikować szczegółową architekturę wojskowego toru detekcji. Dla dydaktyki wystarczy rozumieć wymaganie: tor musi zachować relację czasową między pierwszym maksimum, minimum i drugim maksimum.

Zakres dynamiczny

Kula ognista może być ekstremalnie jasna, a tło przed zdarzeniem może być wielokrotnie słabsze. Przyrząd musi więc radzić sobie z dużym zakresem dynamicznym. Jeśli tor nasyci się na pierwszym piku i długo wraca do pracy, może przegapić minimum oraz drugi impuls.

To jest bardzo praktyczna lekcja dla elektroniki pomiarowej. W klasycznym liczniku impulsów problemem bywa pile-up, czas martwy albo saturacja wzmacniacza. W bhangmetrze analogiczny problem dotyczy fotometrii: zbyt jasny pierwszy impuls może "oślepić" tor, jeśli nie przewidziano odpowiedniego zapasu.

Dlatego w takim przyrządzie ważna jest nie tylko czułość, ale też odzyskiwanie po dużym sygnale, stabilność tła i odporność na jednorazowe zakłócenia.

Dwa kanały i redundancja

Incydent Vela dobrze pokazuje sens redundancji. Według lokalnie wykorzystywanego opisu Gdenarza, sygnał z Vela 6911 był podobny do małego wybuchu jądrowego, ale proporcje sygnałów z dwóch kanałów odbiegały od wcześniejszych, jednoznacznych przypadków.³

To nie unieważnia automatycznie sygnału, ale obniża pewność. Jeśli dwa niezależne lub półniezależne kanały widzą ten sam przebieg w zgodnych proporcjach, przypadek jest mocniejszy. Jeśli widzą coś podobnego, lecz nie w pełni zgodnego, trzeba rozważyć artefakt optyczny, geometrię, uszkodzenie, odbicie, cząstkę uderzającą w sensor albo inne rzadkie zjawisko.

Właśnie dlatego pojedynczy bhangmeter nie jest systemem weryfikacji traktatu. Jest jednym kanałem, który musi zostać zestawiony z sejsmiką, infradźwiękami, hydroakustyką i radionuklidami.

Fałszywe alarmy

Fałszywy alarm nie musi oznaczać błędu przyrządu. Może oznaczać, że przyrząd zobaczył prawdziwy, silny błysk, który nie był wybuchem jądrowym.

Do kandydatów należą:

• pioruny i bardzo silne wyładowania atmosferyczne,
• meteory i bolidy,
• odbicia światła słonecznego od chmur, oceanu albo elementów satelity,
• zdarzenia cząstkowe w elektronice,
• uszkodzenie lub chwilowe przesłonięcie sensora,
• nietypowe zjawiska optyczne związane z geometrią obserwacji.

Z punktu widzenia metrologii ważne jest, żeby fałszywe alarmy nie były traktowane jak wstydliwy dodatek. Są częścią projektu detektora. Każdy próg i każdy wzorzec klasyfikacji ma kompromis: zbyt czuły system zgłasza za dużo kandydatów, zbyt surowy może przeoczyć słabe zdarzenia.

Program Vela

Satelity Vela powstały w kontekście monitorowania zakazu prób w atmosferze i przestrzeni kosmicznej po traktacie częściowego zakazu prób z 1963 roku. Ich zadaniem było wykrywanie sygnałów związanych z detonacjami jądrowymi: promieniowania X, gamma, neutronów, impulsów elektromagnetycznych i sygnałów optycznych.³

W nowszych satelitach programu pojawiły się bhangmetry, czyli szybkie czujniki optyczne szukające podwójnej krzywej jasności. Ten program ma też znany efekt uboczny naukowy: satelity Vela przyczyniły się do odkrycia kosmicznych rozbłysków gamma, bo zarejestrowały sygnały niepasujące do ziemskich detonacji jądrowych.

Ta historia jest dobrym przykładem, jak aparatura tworzona do kontroli zbrojeń może stać się narzędziem nauki podstawowej, o ile dane są analizowane ostrożnie i nie wymusza się jednej interpretacji.

Incydent Vela

22 września 1979 satelita Vela 6911 zarejestrował nad południowym Atlantykiem sygnał optyczny podobny do podwójnego błysku atmosferycznej detonacji jądrowej. Zdarzenie znane jako Incydent Vela albo Alert 747 do dziś pozostaje przedmiotem sporu.³

Najważniejsze nie jest tu rozstrzygnięcie polityczne, lecz metrologia. Sygnał optyczny był bardzo interesujący, ale nie został publicznie domknięty przez równie silne potwierdzenie radionuklidowe, sejsmiczne lub inne. Dodatkowo sporne były szczegóły zgodności kanałów optycznych.³

Dlatego Vela jest świetnym studium przypadku. Pokazuje, że nawet sygnał uważany za "charakterystyczny" nie wystarcza, jeśli brakuje niezależnych danych lub jeśli tor pomiarowy daje nietypowe proporcje.

Relacja do IMS/CTBT

Dzisiejszy IMS CTBTO nie jest po prostu siecią bhangmetrów. Jego podstawowe technologie to sejsmika, hydroakustyka, infradźwięki i radionuklidy.⁴ Detekcja optyczna jest raczej historycznym i satelitarnym kanałem rozpoznania zdarzeń atmosferycznych, który dobrze uzupełnia obraz, ale sam nie daje pełnego potwierdzenia.

Różnica jest logiczna. Bhangmeter widzi światło. IMS próbuje zbudować spójny przypadek: fale lokalizują i klasyfikują zdarzenie, radionuklidy mogą potwierdzić naturę jądrową, a modelowanie atmosferyczne łączy detekcje materiałowe z możliwym obszarem źródła.

W praktyce optyczna krzywa podwójnego błysku jest jednym elementem większego świata detekcji. Może być szybka i bardzo sugestywna, ale w systemie kontroli zbrojeń musi wejść do fuzji danych.

Co odróżnia optykę od radionuklidów

Optyka jest szybka. Jeśli satelita widzi zdarzenie, krzywa jasności powstaje natychmiast. Problem polega na tym, że optyka widzi przede wszystkim skutki promieniowania cieplnego i hydrodynamiki w atmosferze. Nie widzi bezpośrednio składu jądrowego.

Radionuklidy są wolniejsze. Muszą powstać, wydostać się, przemieścić i zostać zebrane przez stację. Ale gdy zostaną wykryte odpowiednie produkty rozszczepienia, ich wartość dowodowa jest inna: mówią o materiale, a nie tylko o świetle.⁴

Dlatego te kanały odpowiadają na inne pytania. Bhangmeter pyta: "czy widzę optyczny profil zgodny z atmosferyczną detonacją jądrową?" Radionuklidy pytają: "czy w powietrzu jest materiał zgodny z reakcją jądrową?"

Tor pomiarowy jako lekcja elektroniki

Bhangmeter jest interesujący dla działu Elektronika, bo łączy kilka klasycznych problemów:

• fotodiodę jako przetwornik światło-prąd,
• szybki wzmacniacz transimpedancyjny,
• szerokie pasmo i stabilność,
• zakres dynamiczny,
• próbkowanie czasowe,
• filtrację tła,
• automatyczne wykrywanie cech przebiegu,
• synchronizację czasu z innymi sensorami.

To są te same problemy, które pojawiają się w aparaturze impulsowej: oscyloskopach, torach scyntylacyjnych, szybkich licznikach i akwizycji FPGA. Różni się tylko czujnik wejściowy i skala zdarzenia.

W bezpiecznej dydaktyce można badać te zagadnienia na diodach LED, lampach błyskowych, migawkach, sztucznych przebiegach i symulowanych danych. Nie ma potrzeby ani sensu budować aparatury do realnej detekcji detonacji.

Jak klasyfikować syntetyczny przebieg

Do ćwiczeń wystarczy prosty model:

S(t) = A1 * exp(-t / tau1) + A2 * (1 - exp(-(t - td) / tau_r)) * exp(-(t - td) / tau_d),

dla t > td w drugim składniku. Pierwszy składnik imituje szybki błysk i zanik. Drugi imituje opóźniony, wolniejszy impuls. Do tego dodaje się tło i szum.

Z takiego syntetycznego przebiegu można wyznaczyć:

• czas pierwszego maksimum,
• czas minimum,
• czas drugiego maksimum,
• stosunek amplitud,
• szerokości impulsów,
• całkę pod krzywą,
• odporność klasyfikacji na szum i nasycenie.

To ćwiczenie uczy detekcji sygnałów bez publikowania praktycznego algorytmu rozpoznawania realnych testów. Dane są sztuczne, a cel jest aparaturowy.

Ograniczenia optycznej detekcji

Optyczny kanał działa najlepiej dla zdarzeń atmosferycznych lub widocznych z przestrzeni. Nie pomoże przy dobrze ukrytym zdarzeniu podziemnym. Może być zasłonięty przez chmury, geometrię Ziemi, położenie satelity albo oślepiające tło.

Jest też wrażliwy na interpretację. Krzywa "podobna" do podwójnego błysku nie jest automatycznie wyrokiem. Trzeba sprawdzić:

• zgodność kanałów,
• geometrię obserwacji,
• stan sensora,
• obecność zjawisk atmosferycznych,
• dane z innych satelitów,
• sygnały falowe,
• ewentualne radionuklidy.

Ta lista jest ważniejsza niż sama nazwa detektora. Bhangmeter jest świetnym przykładem przyrządu, którego wynik ma sens dopiero w systemie.

Podsumowanie

Podwójny błysk nie jest magicznym kodem, lecz konsekwencją fizyki kuli ognistej: najpierw bardzo szybka emisja, potem zasłonięcie przez gorącą falę uderzeniową, potem breakaway i wtórny wzrost jasności. Bhangmeter jest przyrządem, który zamiast obrazu zapisuje tę historię jako krzywą jasności.

Najważniejsza lekcja metrologiczna brzmi: w szybkim torze optycznym liczy się kształt, czas, dynamika i korelacja z innymi kanałami. Dlatego Incydent Vela pozostaje tak dobrym studium przypadku - nie dlatego, że daje prostą odpowiedź, lecz dlatego, że pokazuje, jak ostrożnie trzeba traktować nawet bardzo charakterystyczny sygnał.

Historia programu Vela i ewolucja detekcji optycznej

Program Vela (od hiszpańskiego „czuwać") uruchomiony przez USA w latach 60. XX wieku miał na celu monitorowanie przestrzegania Układu o Częściowym Zakazie Prób Jądrowych (Partial Test Ban Treaty, PTBT, 1963), który zabraniał prób w atmosferze, pod wodą i w przestrzeni kosmicznej. Pierwsze satelity Vela Hotel (1963–1965) wyposażone były głównie w detektory promieniowania jonizującego (gamma, rentgen, neutrony) i czujniki impulsu elektromagnetycznego (EMP), a nie w optyczne bhangmetry.

Seria Advanced Vela, zaczynając od Vela 5A i 5B (1969) oraz 6A i 6B (1970), dodała bhangmetry jako dodatkowy kanał detekcji. Satelity zostały umieszczone na orbitach silnie eliptycznych lub kołowych o wysokości około 120 000 km (prawie jedna trzecia odległości Księżyca), co dawało szerokie pole widzenia i dużą odporność na zakłócenia atmosferyczne.

Wcześniejsze satelity rozpoznania jądrowego skupiały się na przestrzeni kosmicznej i wysokich warstwach atmosfery. Bhangmeter rozszerzył detekcję na atmosferyczne detonacje blisko powierzchni Ziemi, gdzie optyczna sygnatura podwójnego błysku jest możliwa do obserwacji z wysokiej orbity.

Następcami Vela w roli detekcji zdarzeń jądrowych były systemy DSP (Defense Support Program, 1970–2007) z detektorami podczerwieni do wczesnego ostrzegania rakietowego, a także SBIRS (Space-Based Infrared System, od 2011) i GPS Block II, które nosiły dodatkowe czujniki wykrywania zdarzeń nuklearnych (NDS — Nuclear Detonation Detection System). Satelity GPS zawierają do dziś systemy NUDET (Nuclear Detonation Detection), obejmujące zarówno bhangmetry optyczne, jak i detektory promieniowania jonizującego.

Fizyka kuli ognistej: skale czasu i energii

Aby rozumieć sygnaturę bhangmetryczną, konieczne jest bardziej szczegółowe omówienie faz ewolucji kuli ognistej. W ciągu pierwszych mikrosekund po detonacji temperatura wnętrza osiąga dziesiątki milionów kelwinów. Na tym etapie emisja promieniowania obejmuje głównie rentgen i UV — zakres niewidoczny, ale decydujący o transporcie energii.

Sfera izotermiczna (isothermal sphere) to obszar, w którym promieniowanie rentgenowskie opuszcza wnętrze i podgrzewa chłodniejsze warstwy. Jej granica rozszerza się z prędkością przewyższającą dźwięk, ale poniżej prędkości światła. Temperatura na granicy izotermicznej spada do zakresu 300 000–30 000 K, gdzie powietrze zaczyna emitować widzialne światło.

Właśnie gdy temperatura na zewnętrznej granicy spada poniżej ~6 000 K (mniej niż temperatura powierzchni Słońca), jasność optyczna tej zewnętrznej powłoki gwałtownie maleje, mimo że wnętrze jest nadal dużo gorętsze. To jest mechanizm pierwszego minimum w krzywej jasności: wnętrze emituje intensywnie, ale jest zasłonięte przez relatywnie chłodną i nieprzezroczystą powłokę fali uderzeniowej.

Faza breakaway następuje, gdy fala uderzeniowa wystarczająco się rozszerzy i ochłodzi, by stać się przezroczystą dla widzialnego światła. Wtedy obserwator nagle „widzi" gorące wnętrze — stąd gwałtowny wzrost jasności drugiego impulsu. Czas breakaway jest funkcją energii wybuchu: silniejszy wybuch daje dłuższy czas do breakaway (głębsze i szersze wnętrze musi się oziębiać).

Dydaktyczne przybliżenia (nie operacyjne parametry):

• Dla wybuchu o energii rzędu kilkadziesiąt kiloton: breakaway następuje typowo po kilkudziesięciu milisekundach od detonacji.
• Dla wybuchu o energii rzędu megaton: breakaway trwa od kilkuset milisekund do kilku sekund.
• Krzywa jasności jest silnie asymetryczna: pierwszy impuls trwa milisekundy, drugi jest szerszy i niesie więcej całkowitej energii cieplnej.

Te liczby służą wyłącznie do ilościowego zrozumienia skali, nie do charakteryzowania żadnego konkretnego zdarzenia ani urządzenia.

Typy fotodetektorów stosowanych w szybkiej fotometrii

W dydaktyce bhangmetru ważna jest znajomość właściwości elementów światłoczułych, które mogą być stosowane w szybkich torach optycznych. Bezpiecznym kontekstem jest porównanie ogólnych detektorów stosowanych w szybkich pomiarach fizycznych.

Fotodiody PIN (p-intrinsic-n): Szerokie pasmo (do 1–10 GHz przy małej powierzchni), niski poziom szumów prądowych, dobra liniowość w szerokim zakresie mocy optycznej. Ograniczenie: mała powierzchnia aktywna (frakcja mm²), co wymaga dobrej optyki koncentrującej.

Fotodetektory lawinowe (APD, Avalanche Photodiode): Wewnętrzne wzmocnienie przez lawinę nośników (typowo M = 10–100), wyższe SNR przy słabych sygnałach, ale bardziej złożona charakterystyka szumów i czułość na napięcie zasilania i temperaturę. Pasmo zwykle węższe niż PIN.

Fotopowielacze (PMT): Bardzo wysokie wzmocnienie (10⁶–10⁸), dobra odpowiedź na słabe sygnały, ale wolniejszy czas narastania (nanosekudy), czułość na pola magnetyczne i ograniczona odporność na duże sygnały (nasycenie).

Fototranzystory i sensory CMOS/CCD: Wolniejsze i bardziej podatne na szum, ale zintegrowane z elektroniką, tanie i nadające się do celów dydaktycznych.

W historycznym bhangmetrze satelitarnym używano prawdopodobnie specjalizowanych fotodetektorów z szerokim polem widzenia i odpornością na środowisko kosmiczne, ale szczegóły techniczne takich systemów pozostają chronione. Dla celów dydaktycznych wystarczy rozumieć fundamentalne ograniczenia każdej klasy detektorów.

Przetwarzanie sygnału: od prądu do wyniku klasyfikacji

Droga od prądu fotodetektora do decyzji klasyfikacyjnej obejmuje kilka etapów, każdy z własnymi możliwościami i ograniczeniami:

1. Transimpedancja: Prąd fotodetektora (rzędu mikroamperów lub miliamperów dla dużych sygnałów) jest przetwarzany na napięcie przez wzmacniacz transimpedancyjny (TIA — Transimpedance Amplifier). Pasmo TIA określa, jak dobrze zachowany jest kształt sygnału. Kompromis: wyższe pasmo → większy szum odniesienia (Johnson + szum wejściowy).

2. Próbkowanie: Przetwornik ADC próbkuje napięciowy przebieg z określoną częstotliwością. Twierdzenie Nyquista mówi, że częstotliwość próbkowania musi być co najmniej dwukrotnie wyższa od najwyższej częstotliwości w sygnale (f_s ≥ 2f_max). Dla przebiegu z czasem trwania pierwszego impulsu rzędu milisekund wystarczy częstotliwość próbkowania ~10 kHz. Dla rozdzielczości kilkudziesięciu mikrosekund potrzebne jest ~100 kHz. Dla nanosekundowych szczegółów — GHz.

3. Ekstrakcja cech: Następuje identyfikacja lokalnych ekstremów krzywej — czas t₁ pierwszego maksimum, czas t₂ minimum, czas t₃ drugiego maksimum, amplitudy A₁ i A₂, stosunek A₁/A₂, szerokości połówkowe impulsów.

4. Klasyfikacja: Porównanie wektora cech z oczekiwanymi wzorcami dla różnych klas zdarzeń. W najprostszym modelu jest to progowanie kilku warunków jednocześnie (klasyfikator regułowy). W bardziej zaawansowanych systemach może to być uczenie maszynowe, ale podstawową zasadą zawsze pozostaje: kształt i czas, nie sama amplituda.

5. Korelacja z innymi kanałami: Wynik jednego sensora wchodzi do systemu wielokanałowego. Pewność finalna zależy od zgodności wszystkich dostępnych kanałów.

Tło kosmiczne i fałszywe alarmy satelitarne

Satelity na wysokich orbitach są stale bombardowane przez cząstki promieniowania kosmicznego (protony, elektrony z pasa Van Allena, ciężkie jony z promieniowania galaktycznego). Cząstka taka, trafiając w detektor, może wywołać impuls elektryczny trudny do odróżnienia od sygnału optycznego. Zjawisko to znane jest jako single-event upset (SEU) lub single-event transient (SET) w elektronice kosmicznej.

Dla bhangmetrów kluczowym zagadnieniem jest odporność na te fałszywe impulsy. Mechanizmy filtracji obejmują:

• Redundancję: dwa identyczne detektory na tym samym satelicie — jednoczesny fałszywy impuls w obu jest mało prawdopodobny przy kosmicznym tle.
• Kształt impulsu: impuls z cząstki kosmicznej ma zazwyczaj inny kształt niż przebieg podwójnego błysku (brakuje struktury czasowej z dwoma ekstremami).
• Amplitudę i czas trwania: skrajnie krótkie impulsy (nanosekudy) mogą być filtrowane jako niespójne z oczekiwanym sygnałem optycznym.
• Korelację ze zdarzeniami cząstkowymi identyfikowanymi przez inne instrumenty na pokładzie.

Właśnie przez ten kontekst satelitarny fragment historii Vela jest metrologicznie tak bogaty: jeden z argumentów kwestionujących interpretację sygnału z 1979 roku wskazywał na możliwość, że jeden z czujników optycznych mógł zarejestrować zdarzenie cząstkowe zamiast sygnału optycznego.

Porównanie bhangmetru z innymi kanałami IMS CTBTO

Dla pełnego obrazu systemu detekcji detonacji jądrowych warto zestawić bhangmetr z czterema oficjalnymi kanałami Międzynarodowego Systemu Monitoringu (IMS), który funkcjonuje w ramach Traktatu o Całkowitym Zakazie Prób Jądrowych (CTBT):

Kanał	Technologia	Czas do wyniku	Zasięg geograficzny	Dowód na charakter jądrowy	Ograniczenia
Bhangmeter (satelity)	Fotometria optyczna	Sekundy	Globalna (atmosfera)	Pośredni (kształt błysku)	Chmury, fałszywe alarmy optyczne, kosmiczne
Sejsmika	Sieć 170 stacji IMS	Minuty	Globalny	Pośredni (sygnatura sejsmiczna)	Nie odróżnia jądrowego od konwencjonalnego bez danych uzupełniających
Hydroakustyka	11 stacji IMS	Minuty–godziny	Ocean globalny	Pośredni (sygnatura ciśnieniowa)	Tylko zdarzenia bliskie oceanów lub wody
Infradźwięki	60 stacji IMS	Godziny	Globalna (atmosfera)	Pośredni (długofalowa sygnatura)	Zmienna czułość zależna od pogody
Radionuklidy	80 stacji pomiarowych	Dni–tygodnie	Globalna (po dyspersji atmosferycznej)	Bezpośredni (produkty rozszczepienia/aktywacji)	Czas transportu, rozcieńczenie, maskowanie

Widać wyraźnie, że bhangmeter jest najszybszy, ale ma najmniej bezpośredni dowód na charakter jądrowy zdarzenia. Radionuklidy są najwolniejsze, ale dają najsilniejszą bezpośrednią informację o naturze zdarzenia. Dlatego właśnie wielokanałowa fuzja danych jest fundamentem systemu weryfikacyjnego.

Trzy przykłady obliczeniowe

Przykład 1: Wyznaczanie czasu narastania pierwszego impulsu

Syntyczny przebieg bhangmetryczny ma pierwszy impuls opisany jako:
S₁(t) = A₁ · exp(-t / τ₁) dla t > 0.

Przyjmij A₁ = 10 000 jednostek umownych, τ₁ = 5 ms. Oblicz czas, po którym sygnał spada do 50%, 10% i 1% wartości szczytowej.

Rozwiązanie: t = −τ₁ · ln(S/A₁).

• t₅₀% = −5 ms · ln(0,5) = 5 ms · 0,693 = 3,47 ms
• t₁₀% = −5 ms · ln(0,1) = 5 ms · 2,303 = 11,51 ms
• t₁% = −5 ms · ln(0,01) = 5 ms · 4,605 = 23,03 ms

Wniosek: przy τ₁ = 5 ms sygnał od wartości szczytowej do 1% amplitudy spada w ciągu ~23 ms. Tor elektroniczny musi poprawnie rejestrować sygnał przez co najmniej 50–100 ms, aby nie obciąć ogona pierwszego impulsu i uchwycić minimum. Jeżeli czas narastania detektora jest zbyt długi (np. 10 ms), kształt pierwszego piku zostanie zniekształcony.

Przykład 2: Detekcja minimum w obecności szumu

Przebieg ma minimum w chwili t_min = 50 ms z wartością S_min = 500 j.u. Standardowe odchylenie szumu wynosi σ = 150 j.u. Oblicz stosunek sygnał-szum (SNR) w minimum i oceń, ile pomiarów trzeba uśrednić, by SNR wzrósł do 10.

Rozwiązanie:
SNR = S_min / σ = 500 / 150 = 3,33.

SNR rośnie jak √N po uśrednieniu N pomiarów. Wymagane N = (10 / 3,33)² = (3,003)² ≈ 9.

Interpretacja: przy uśrednieniu 9 próbek (z okna np. ±4 próbki wokół t_min) SNR w minimum osiąga 10, co jest bezpiecznym progiem do identyfikacji struktury przebiegu. Dla okna 9 próbek przy próbkowaniu 1 kHz odpowiada to integracji przez 9 ms — akceptowalne, jeśli kształt przebiegu zmienia się na skali >20 ms.

Przykład 3: Klasyfikacja na podstawie stosunku amplitud

Baza danych sygnatur zawiera 50 historycznych zdarzeń optycznych skalsyfikowanych jako: detonacje jądrowe (N=8), pioruny kuliste (N=15), bolidy (N=12), odbicia słoneczne (N=15). Wyznaczono stosunek amplitud A₁/A₂ dla każdego zdarzenia:

Klasa	Średnia A₁/A₂	Odchylenie std.
Detonacja	0,35	0,08
Piorun kulisy	0,85	0,20
Bolid	1,20	0,40
Odbicie słoneczne	brak minimum	—

Nowy sygnał ma A₁/A₂ = 0,40. Oblicz prawdopodobieństwo przynależności do klasy „detonacja" i „piorun" przy założeniu rozkładu normalnego i równych a priori (P_prior = 0,25 dla każdej klasy).

Gęstość prawdopodobieństwa dla A₁/A₂ = 0,40:

• p(x|detonacja) = N(0,35; 0,08) → z = (0,40 − 0,35)/0,08 = 0,625 → φ(z)/σ ≈ 0,323/0,08 = 4,04
• p(x|piorun) = N(0,85; 0,20) → z = (0,40 − 0,85)/0,20 = −2,25 → φ(z)/σ ≈ 0,032/0,20 = 0,16
• p(x|bolid) ≈ N(1,20; 0,40) → z = (0,40 − 1,20)/0,40 = −2,0 → φ(z)/σ ≈ 0,054/0,40 = 0,135
• p(x|odbicie) = 0 (brak minimum → A₁/A₂ niezdefiniowane)

Prawdopodobieństwa a posteriori (Bayes):

• mianownik: 4,04 + 0,16 + 0,135 + 0 = 4,335
• P(detonacja|x) = (4,04 × 0,25) / (4,335 × 0,25 + ...) = 4,04 / 4,335 ≈ 93,2%
• P(piorun|x) ≈ 3,7%

Wniosek dydaktyczny: stosunek amplitud A₁/A₂ = 0,40 silnie wskazuje na detonację, ale nawet przy 93% prawdopodobieństwie bayesowskim jeden z kanałów diagnostycznych nie wystarczy do ostatecznej klasyfikacji — potrzebna jest korelacja z innymi kanałami IMS.

Pytania otwarte

1. W jaki sposób zachmurzenie wpływa na sygnaturę bhangmetryczną? Czy chmury mogą jednocześnie tłumić amplitudę i modyfikować kształt przebiegu? Jak możliwa jest korekcja takiego efektu?

2. Jak zmienia się sygnatura optyczna dla wybuchu wysokościowego (above-ground, altitude > 10 km) w porównaniu do wybuchu przy powierzchni? Czy pierwszy impuls i breakaway są podobne, czy fundamentalnie różne?

3. W jakich warunkach geometrycznych bhangmeter może obserwować tę samą detonację pod różnymi kątami? Jak kąt obserwacji wpływa na zmierzone amplitudy i czasy?

4. Jakie są fizyczne ograniczenia zasięgu bhangmetru? Czy możliwe jest wykrycie detonacji za horyzontem geometrycznym satelity? Jaką rolę gra refrakcja atmosferyczna?

5. Jak system klasyfikacyjny powinien traktować zdarzenia, w których jedno kryterium klasyfikacyjne wskazuje na detonację, a inne — na zdarzenie naturalne? Jak formalnie modelować taką niepewność?

6. W jaki sposób promieniowanie kosmiczne i cząstki energetyczne w środowisku Ziemi mogą uszkodzić detektor optyczny satelity po wieloletniej pracy? Jak monitoruje się degradację toru?

7. Czy bhangmetr mógłby wykryć podziemne zdarzenie, które częściowo wydostało się przez komin (vent) lub wywołało widoczny efekt na powierzchni? Jakie byłyby cechy takiego sygnału?

8. Jak interpretować sygnał bhangmetru dla zdarzenia o nieznanej lokalizacji, gdy brak jest korelacji z innymi kanałami IMS? Jaki jest minimalny zestaw niezależnych danych do wydania oceny?

Podsumowanie dydaktyczne

1. Kształt, nie amplituda — bhangmeter rozpoznaje zdarzenie przez czasowy profil jasności, nie przez samą wartość maksymalną. Podwójny błysk (pierwszy pik → minimum → drugi pik) jest sygnaturą wynikającą z fizyki kuli ognistej w atmosferze: sfera izotermiczna, zasłona nieprzezroczystej fali uderzeniowej i breakaway.

2. Fizyka podyktowuje skalę — czasy charakterystyczne podwójnego błysku (milisekundy do sekund) zależą od energii wybuchu, gęstości atmosfery i geometrii obserwacji. Rozumienie tych zależności pozwala poprawnie projektować wymagania dla toru pomiarowego (pasmo, czas martwy, zakres dynamiczny).

3. Zakres dynamiczny i nasycenie — tor bhangmetru musi radzić sobie z ogromną różnicą jasności między pierwszym pikiem a tłem. Nasycenie toru na pierwszym piku prowadzi do utraty informacji o minimum i drugim piku — klasyczny błąd aparaturowy niebezpieczny dla klasyfikacji.

4. Redundancja kanałów — jeden sensor optyczny nie wystarczy. Redundancja (dwa niezależne detektory na tym samym satelicie) filtruje zdarzenia cząstkowe i poprawia wiarygodność. Brak zgodności kanałów optycznych w przypadku Vela 1979 był jednym z głównych powodów wątpliwości co do interpretacji.

5. Fuzja wielokanałowa — bhangmeter jest szybki, ale daje tylko pośredni dowód na charakter jądrowy zdarzenia. Ostateczna interpretacja wymaga korelacji z sejsmiką, hydroakustyką, infradźwiękami i — najważniejsze — radionuklidami, które jako jedyne dają bezpośredni dowód na reakcję jądrową.

6. Incydent Vela jako studium przypadku — zdarzenie z 22 września 1979 r. nad południowym Atlantykiem pokazuje, jak niejednoznaczny może być nawet sygnał o pozornie charakterystycznym kształcie. Brak zgodności między kanałami i brak publicznego potwierdzenia radionuklidowego pozostawiły sprawę otwartą przez dekady.

7. Środowisko satelitarne — kosmiczne promieniowanie jonizujące jest stałym źródłem fałszywych impulsów w detektorach satelitarnych. Projektowanie bhangmetrów musi uwzględniać odporność na SEU/SET, co przekłada się na wybór redundancji, kształtu filtrowania i kryteriów przynależności do klasy sygnału.

8. Bezpieczna dydaktyka — wszystkie opisane tu zasady (kształt przebiegu, parametryzacja modelu, klasyfikacja bayesowska, redundancja) można przećwiczyć na syntetycznych danych z lampami LED, oscyloskopem i mikrokontrolerem. Ćwiczenia dydaktyczne nie potrzebują żadnej wiedzy o geometrii ani parametrach prawdziwych urządzeń jądrowych — liczy się rozumienie metrologii szybkich sygnałów optycznych.

Naukowy efekt uboczny programu Vela: odkrycie rozbłysków gamma

Historia programu Vela ma zaskakujące zakończenie naukowe. W 1973 roku astrofizycy Ray Klebesadel, Ian Strong i Roy Olson z Los Alamos National Laboratory opublikowali artykuł opisujący krótkie, intensywne rozbłyski promieniowania gamma, których źródło znajdowało się poza Układem Słonecznym. Dane pochodziły właśnie z satelitów Vela — instrumenty rejestrujące promieniowanie gamma przeznaczone do wykrywania prób jądrowych zarejestrowały zdarzenia zupełnie innego rodzaju.

Kosmiczne rozbłyski gamma (GRB — Gamma-Ray Bursts) okazały się jednym z najbardziej energetycznych zjawisk we Wszechświecie: krótkie (ułamki sekund do kilku minut) błyski promieniowania gamma, których źródłami są katastroficzne zdarzenia astrofizyczne (zderzenia gwiazd neutronowych, hiperowe supernowe). Do dziś GRB są intensywnie badane przez satelity takie jak Swift i Fermi-GBM.

Odkrycie GRB przez satelity Vela jest klasycznym przykładem nauki „z przypadku" — aparatura stworzona do celów wojskowych i kontroli zbrojeń przyniosła fundamentalne odkrycie w astrofizyce. Podobny wzorzec można znaleźć w historii detektorów podczerwieni (zbudowanych pierwotnie do wczesnego ostrzegania rakietowego, a stosowanych do monitorowania aktywności wulkanicznej i erupcji pożarów) czy w historii systemu GPS (militarny system nawigacji, który stał się globalną infrastrukturą cywilną).

Dydaktyczna lekcja z tego odkrycia jest prosta: przyrząd zaprojektowany na określony sygnał może napotkać zupełnie inne zjawisko. Dobra aparatura powinna rejestrować dane szeroko, a interpretacja powinna być otwarta na nieoczekiwane. Program Vela, próbując wykrywać jądrowe próby, niejako „przy okazji" odkrył nową klasę astronomicznych zjawisk i zrewolucjonizował astrofizykę wysokich energii.

Warto dodać, że GRB rejestrowane przez Vela miały profile czasowe i energetyczne zupełnie inne niż oczekiwana sygnatura detonacji jądrowej: brak charakterystycznego podwójnego impulsu optycznego, intensywna emisja w zakresie gamma (nie widzialnym), brak korelacji z jakimkolwiek wydarzeniem naziemnym. To właśnie te różnice — zarejestrowane precyzyjnie przez kalibrowaną aparaturę — pozwoliły na wykluczenie ziemskiego pochodzenia i sfomułowanie hipotezy kosmicznej. Metrologia i dobra dokumentacja danych okazały się kluczowe dla przełomowego odkrycia naukowego. Naukowe otwarcie archiwów Vela po dekadach utajnienia jest świadectwem wartości długoterminowego przechowywania surowych danych pomiarowych — nawet wówczas, gdy pierwotne pytanie badawcze zostało zastąpione przez zupełnie inne.

Dodatkowe materiały multimedialne

Ćwiczenia praktyczne

Pierwsze ćwiczenie: na podstawie syntetycznej tabeli t, S(t) znaleźć pierwszy pik, minimum i drugi pik. Następnie obliczyć stosunek amplitud oraz czas między pikami.

Drugie ćwiczenie: dodać do tego samego przebiegu szum gaussowski i sprawdzić, przy jakim stosunku sygnału do szumu minimum przestaje być jednoznaczne. Celem jest zrozumienie fałszywych alarmów, nie rozpoznawanie realnych detonacji.

Trzecie ćwiczenie: porównać trzy przebiegi syntetyczne: pojedynczy błysk LED, błysk z podwójnym impulsem oraz przebieg nasycony, w którym pierwszy pik "przykleja się" do maksimum toru. Wskazać, który przypadek najłatwiej błędnie zinterpretować.

Czwarte ćwiczenie: zbudować bezpieczny fotometr dydaktyczny z fotodiodą lub fototranzystorem i rejestrować błyski diody LED sterowanej mikrokontrolerem. Zmierzyć, jak zmienia się przebieg po dodaniu filtru RC o różnych stałych czasowych.

Piąte ćwiczenie: dla przypadku Vela wypisać, jakie niezależne dane byłyby potrzebne do zwiększenia pewności interpretacji: drugi satelita, zgodność kanałów, sejsmika, infradźwięki, hydroakustyka, radionuklidy i model atmosferyczny.

Szóste ćwiczenie: porównać bhangmeter z czterema technologiami IMS. Dla każdego kanału wskazać: co mierzy, jak szybko daje wynik, jakie ma tło i czy może samodzielnie potwierdzić jądrową naturę zdarzenia.

Przejdź do ćwiczenia interaktywnego

Powiązane artykuły

• Vela 1979 jako problem metrologii, a nie tylko historii
• Międzynarodowy System Monitoringu CTBT: sejsmika, infradźwięki, hydroakustyka i radionuklidy
• Kula ognista i Wilson Cloud
• Impuls elektromagnetyczny (EMP)
• Oscyloskop w pomiarach impulsów nanosekundowych: pasmo, czas narastania i sondy