Gdy inżynierowie poszukują ropy, gazu lub budują tunele, polegają na falach sejsmicznych — niewielkich drganiach wysyłanych przez grunt — aby ujawnić, co znajduje się poniżej. Jednak fale te nie przemieszczają się bez zmian: zwalniają i tracą energię, gdy przechodzą przez różne skały. W badaniu tym analizowano, jak suchy węgiel, a zwłaszcza drobne ziarna tworzące jego szkielet, przekształcają te fale. Łącząc precyzyjne pomiary laboratoryjne z symulacjami komputerowymi, autorzy pokazują, jak zderzenia ziaren, tarcie i mieszanka rozmiarów ziaren w węglu kontrolują prędkość i wygaszanie fal sejsmicznych, dając wskazówki do lepszej obrazowania podziemnego i bezpieczniejszego wydobycia zasobów. Rysunek 1.
Próbki poddane drganiom, by zbadać duże zagadnienia
Naukowcy zaczęli od prawdziwych kawałków węgla z dwóch basenów węglowych w Chinach: węgla wysokozasiarczonego (o większym stopniu metamorfozy), bardziej zwartego i dojrzałego, oraz węgla niskorankowego, młodszego i luźniej zgranulowanego. Pocięli te węgle na małe cylindry i wykonali odpowiadające cylindry z dwóch materiałów do drukowania 3D: żywicy fotoczułej o gumowatej charakterystyce oraz sztywniejszego tworzywa PLA. Wszystkie próbki były starannie wysuszone, uszczelnione i wyposażone w tensometry, a następnie zamocowane w niestandardowym stanowisku do badań niskoczęstotliwościowych, które delikatnie ściskało je w zakresie częstotliwości od 1 do 250 Hz — w przybliżeniu odpowiadającym pasmu stosowanemu w badaniach sejsmicznych. Porównując, jak próbki się rozciągały i ściskały, zespół mógł obliczyć, jak szybko fale sprężyste (fale P) przemieszczają się przez każdą próbkę oraz jak silnie te fale są tłumione.
Jak węgiel wygląda pod mikroskopem
Obrazy mikrostruktury węgla wyjaśniają, dlaczego różne węgle inaczej oddziałują na fale. Węgiel o wyższym stopniu metamorfozy ma ziarna podobnej wielkości, ściśle i równomiernie upakowane, pozostawiając głównie drobne, izolowane pory. Struktura ta odzwierciedla intensywne zagęszczenie i przemiany chemiczne w czasie. Węgiel niskorankowy natomiast wykazuje szeroką mieszankę rozmiarów ziaren, luźniejsze upakowanie i wiele większych, dobrze połączonych porów. To nieuporządkowane ułożenie pozwala ziarnom poruszać się, zderzać i ślizgać się łatwiej, gdy przechodzi fala, tworząc więcej możliwości wydzielania energii z fali. Te wizualne różnice pomagają wyjaśnić, dlaczego węgiel niskorankowy wykazuje silniejszą zależność prędkości fali od częstotliwości oraz większe tłumienie niż węgiel wysokorankowy.
Symulowanie zderzeń ziaren jedno cząstka po drugiej Rysunek 2.
Aby zajrzeć do środka procesu, autorzy zbudowali model komputerowy traktujący węgiel nie jako gładki blok, lecz jako tysiące małych sferycznych cząstek połączonych ze sobą. W tym modelu elementów dyskretnych każde ziarno może naciskać, ciągnąć i ślizgać się względem sąsiadów, a specjalne człony tłumiące reprezentują utratę energii podczas zderzeń normalnych i ruchu stycznego (ślizgowego). Przeprowadzając wirtualne testy ściskania w szerokim zakresie częstotliwości, stwierdzili, że zwiększanie tych członów tłumiących oraz nierównomierność rozkładu rozmiarów cząstek zarówno zmniejsza prędkość fali P, jak i znacznie zwiększa tłumienie. Tłumienie styczne — związane z tarciem ślizgowym — okazało się szczególnie ważne, powodując mniej więcej trzy do czterech razy większe straty energii niż tłumienie normalne. Gdy wszystkie tłumienia ustawiono na zero, fale przemieszczały się najszybciej i prawie nie wykazywały dyspersji ani tłumienia.
Wydrukowane skały jako kontrolowane stanowiska badawcze
Modele drukowane 3D pełnią rolę uproszczonych, kontrolowanych wersji skał. Wydruk z żywicy zachowuje się jak wysoce lepkie, gumowate ciało stałe: ma gęstą strukturę, wysoką wartość współczynnika Poissona i silne tarcie wewnętrzne, co prowadzi do wyraźnej zależności prędkości fali od częstotliwości oraz wysokiego tłumienia. Wydruk z PLA, wykonany metodą wytłaczania, jest bardziej sztywny i zachowuje się bliżej klasycznego ciała sprężystego, z mniejszym tarciem wewnętrznym i słabszym tłumieniem. W rezultacie wykazuje mniejsze zmiany prędkości fali wraz z częstotliwością i niższe tłumienie. Porównanie tych materiałów syntetycznych z naturalnymi węglami potwierdziło, że zarówno tłumienie na poziomie ziaren, jak i równomierność rozmiarów ziaren odgrywają kluczowe role w kształtowaniu odpowiedzi sejsmicznej. Symulacje z użyciem modelu zespalanych cząstek odtworzyły ogólne trendy obserwowane w eksperymentach, choć drobne szczegóły dotyczące tłumienia pozostają trudne do dokładnego dopasowania.
Co to oznacza dla interpretacji sygnałów sejsmicznych
Dla laików najważniejszy wniosek jest taki, że w suchym węglu to grzechotanie i ślizganie się stałych ziaren — nie tylko płyny w porach — może silnie spowalniać i osłabiać fale sejsmiczne, szczególnie przy pewnych częstotliwościach. Węgiel niskorankowy, luźno upakowany z szeroką mieszanką rozmiarów ziaren, zachowuje się jak lepszy „amortyzator” niż ściśle upakowany węgiel wysokorankowy. Zrozumienie, jak tarcie styczne, zderzenia normalne i rozkład rozmiarów ziaren kontrolują zachowanie fal, pomaga geofizykom wybierać lepsze modele przy interpretacji danych sejsmicznych w obszarach bogatych w węgiel, poprawiając oszacowania własności skał i zmniejszając niepewność w badaniach podpowierzchniowych.
Cytowanie: Chen, H., Zou, G., Feng, X. et al. Experimental and numerical investigation of elastic wave dispersion and attenuation induced by coal particle damping.
Sci Rep16, 6033 (2026). https://doi.org/10.1038/s41598-026-36113-0
Słowa kluczowe: mikrostruktura węgla, tłumienie fal sejsmicznych, tłumienie przez cząstki, modelowanie metodą elementów dyskretnych, wydrukowane 3D próbki skał