Syntetyczne modelowanie ciśnienia kapilarnego przy użyciu fal sonicznych dipolowych do poprawy charakteryzacji złóż

Dlaczego pory w skale mają znaczenie dla energii i wody

Głęboko pod ziemią ropa, gaz i woda przemieszczają się przez drobne przestrzenie między ziarnami skały. Rozmiar i kształt tych porów kontrolują, jak łatwo płyny mogą być przechowywane i wydobywane; są one odwzorowane w „krzywych ciśnienia kapilarnego”, na które inżynierowie polegają przy planowaniu odwiertów i zarządzaniu złożami. Mierzenie tych krzywych zwykle wymaga kosztownych, czasochłonnych badań laboratoryjnych na rzadkich próbkach rdzeni. W artykule przedstawiono sposób syntetycznego generowania takich krzywych kapilarnych, wykorzystujący jedynie fale dźwiękowe rejestrowane w odwiertach oraz matematycznie realistyczny obraz sieci porów skały.

Słuchając skał zamiast je krając

Tradycyjne badania złóż w dużym stopniu opierają się na rdzeniach: cylindrycznych fragmentach skały pobieranych z podłoża i badanych metodami takimi jak iniekcja rtęci czy wirowanie. Testy te ujawniają, jak płyny wchodzą i wychodzą z porów przy różnych ciśnieniach, lecz wymagają nienaruszonych próbek, specjalistycznego sprzętu oraz znacznego czasu i kosztów. W wielu złożach próbki rdzeniowe są nieliczne albo ich brak, co pozostawia kluczowe pytania o porowatość, przepuszczalność i rozmieszczenie płynów bez odpowiedzi. Z drugiej strony, logi wiertnicze — szczególnie logi soniczne dipolowe — są rutynowo rejestrowane na całej długości otworu wiertniczego. Narzędzia te wysyłają fale ściskające (P) i ścinające (S) w formację i rejestrują ich prędkość, dostarczając bogatego, choć pośredniego obrazu sztywności i wewnętrznej struktury skały.

Budowanie realistycznego obrazu sieci porów

Autorzy łączą pomiary soniczne z geometrią porów za pośrednictwem modelu komputerowego traktującego pory jako drobne wnęki osadzone w elastycznym ciele stałym. Kształty porów opisują mieszaniną prostych form „euklidesowych” (takich jak gładkie okrągłe lub gwiaździste otwory) oraz bardziej złożonych struktur „fraktalnych”, które naśladują chropowate, rozgałęzione obrysy widoczne na obrazach mikroskopu elektronowego piaskowców. Z użyciem metody elementów brzegowych symulują, jak te mieszane populacje porów odkształcają się przy zmianach ciśnienia, a z tej odpowiedzi wyprowadzają właściwości sprężyste masy i porowatość. Kluczowym postępem jest narzucenie fraktalnych praw skalowania dla populacji małych porów: zamiast dopuszczać wiele różnych rozkładów porów, które wszystkie pasują do tych samych danych sprężystych, model zmusza liczbę porów o danym rozmiarze do podążania za zmierzoną, powtarzalną zależnością. To w dużym stopniu ogranicza długoletni problem, że bardzo różne struktury porów mogą wyglądać identycznie dla fal sonicznych.

Od kształtów porów do zachowania kapilarnego

Gdy uzyska się wiarygodną sieć porów, ramy metody przekształcają ją w krzywą ciśnienia kapilarnego. Najpierw łączą skumulowaną objętość porów z nasyceniem płynu, traktując całkowitą porowatość jako 100% wypełnienia. Następnie obwód każdego poru przeliczają na „równoważny promień gardzieli” — największe koło, które można wcisnąć w otwór porowy — korzystając z metod numerycznych śledzących wpisane okręgi wzdłuż długiej osi poru. Zastosowanie standardowego związku fizycznego między promieniem gardzieli a ciśnieniem kapilarnym przekształca cały rozkład wielkości porów w przewidywaną krzywą kapilarną. Ponieważ reguły fraktalne wiążą obfitość małych, odkształcalnych porów z obserwowanym zachowaniem sprężystym, powstała krzywa nie jest jedynie dopasowaniem do danych; jest bezpośrednią konsekwencją geometrii porów wymaganej do zgodności z logami sonicznymi.

Testowanie metody na rzeczywistych polach naftowych

Naukowcy zastosowali swoje podejście do dwóch piaskowcowych złóż o bardzo różnych dostępnych danych. W polu naftowym Tapi–TTT w Ekwadorze brakowało niektórych logów, więc rekonstrukowali porowatość i gęstość, używając ustalonych zależności między czasem przelotu fali sonicznej a właściwościami skały. W formacji Las Piedras w Wenezueli dostępny był bardziej kompletny i wyższej jakości zestaw logów. W obu przypadkach skalibrowali model sprężysty za pomocą danych sonicznych dipolowych, zbudowali mieszane populacje porów euklidesowo–fraktalnych i wygenerowali syntetyczne krzywe kapilarne dla rtęci. Krzywe te ściśle odzwierciedlały pomiary laboratoryjne, chwytając nie tylko ogólną porowatość, lecz także kluczowe cechy, takie jak udział dużych „megaporów”, nieusuwalne nasycenie wodą oraz przejście od płynów związanych do swobodnych. Błędy mieściły się w kilku procentach dla fraktalnych wymiarów typowych dla piaskowców, a odchylenia rosły tylko wtedy, gdy parametry fraktalne przesunięto poza zakres wspierany przez niezależne badania obrazowe.

Co to oznacza dla przyszłych badań podpowierzchniowych

Dla czytelnika niebędącego specjalistą praktyczny wniosek jest taki, że inżynierowie coraz częściej mogą „usłyszeć” to, co skała powiedziałaby w laboratorium, bez konieczności wynoszenia jej na powierzchnię. Łącząc realistyczne fraktalne opisy ścian porów z dobrze ugruntowaną fizyką akustyczną, metoda ta przekształca rutynowe logi soniczne w syntetyczne krzywe ciśnienia kapilarnego zgodne z benchmarkami laboratoryjnymi. Przy odpowiednim pokryciu logami oferuje szybszą, tańszą i szerzej stosowalną alternatywę dla tradycyjnych badań opartych na rdzeniach, szczególnie w złożach, gdzie próbki są rzadkie lub trudne w obróbce. Chociaż obecne badanie koncentruje się na czystych, dobrze zagęszczonych piaskowcach, autorzy wskazują, jak tę samą strategię można by rozszerzyć na niektóre wapienie i łupki po zmierzeniu ich właściwości fraktalnych, co obiecuje szersze zastosowanie tego „wirtualnego laboratorium” w przyszłych projektach związanych z energią i zasobami wodnymi.

Cytowanie: Galarza-Alava, J., Mendoza-Sanz, J. Synthetic capillary pressure modeling with dipole sonic waves for enhanced reservoir characterization. Sci Rep 16, 11697 (2026). https://doi.org/10.1038/s41598-026-47717-x

Słowa kluczowe: ciśnienie kapilarne, logi soniczne dipolowe, fraktalna struktura porów, charakteryzacja złóż, porowatość piaskowca