Aplikacja 'PA2PPS' powstała w czasie prac nad doskonaleniem systemu GeoWin w ramach projektu celowego:
Nowe aplikacje w polskim systemie GeoWin do przetwarzania i interpretacji profilowań geofizyki otworowej
finansowanego przez Ministerstwo Nauki i Szkolnictwa Wyższego oraz Spółkę Geofizyka Kraków w okresie 02.01.2006-28.02.2007.
Aplikacja ta służy do poprawienia czasu interwałowego na podstawie wyników uzyskanych z profilowania prędkości średnich, PPS
lub tzw. pomiaru Check shots. W wyniku pracy aplikacji uzyskuje się poprawione profilowanie czasu interwałowego. W aplikacji
tej można także utworzyć i wydrukować raport prezentujący wyniki. Raport jest także zapisywany w postaci pliku ASCII zawierającego
nagłówek z ogólnymi informacjami o otworze oraz tabelę z nastepującymi kolumnami: L.P. - liczba porządkowa, H - głębokość,
DT - czas interwałowy z profilowania akustycznego, DTcorr - czas interwałowy poprawiony, Hodogr. - hodograf z PPS, oraz TTI' -
sumaryczny czas interwałowy z dodaną wartością czasu To, potrzebnego na przejśćie fali akustycznej na odcinku od powierzchni Ziemi
do początkowej głębokości, na której rozpoczęto rejestrację DT.
Algorytm korekty profilowania akustycznego został opracownany na podstawie książki:
Seismic Surveing and Well Logging, S. Boyer and J-L. Mari, 1997, Editions Technip.
W aplikacji PA2PPS wykorzystuje się następujące dane wejściowe:
Czas interwałowy z profilowania akustycznego,
Wynik Profilowania Prędkości Średnich lub Check Shot'ów,
Profilowanie krzywizny otworu,
Wartość prędkości VSMP zadawana przez interpretatora, wartość ta jest równa lub zbliżona do
prędkości w strefie małych prędkości badanej w PPS,
Wartość DTmin, zadawana przez interpretatora, dla DT>DTmin oblicza się poprawkę wg algorytmu 2. W aplikacji
PA2PPS wartość ta jest obliczana automatycznie podczas wyboru odcinków o tym samycm trendzie dryftu na wykresie Dryft.
Przetwarzanie profilowania akustycznego obejmuje następujące operacje:
1. Na podstawie zbioru próbek DT, zarejestrowanych dla zadanego kroku głębokościowego oblicza się czas sumaryczny,
TTI (Transit Time Integrated). Jest to zbiór próbek określających czas potrzebny na dojście fali akustycznej do zadanej głębokości
(zakładając, że fala biegnie po ścianie otworu).
TTI(ni*d)=SUMA(d*DT(Hi) (1).
Gdzie: d – krok pomiarowy przy PA, [m],
DT - wartość czasu interwałowego na danej głębokości, Hi, [mikrosekundy/m],
ni – bieżący numer próbki, Hi=ni*d.
W wyniku sumowania uzyskuje się krzywą zmian czasu sumarycznego w funkcji głębokości. Krzywa ta zaczyna się na głębokości, na której występuje pierwsza poprawna próbka DT, HPA.
2. Sprawdzenie, na jakiej głębokości, HPPS występuje pierwsza poprawnie zarejestrowana wartość czasu z PPS,
zwykle głębokość HPPS jest mniejsza od głębokości HPA,
3. obliczenie różnicy: H0 = HPA - HPPS (2)
4. obliczenie czasu T0 = H0/VSMP (3)
gdzie: VSMP – jest to prędkość w interwale między HPPS i HPA, którą zadaje interpretator .
5. obliczenie czasu TTI’: TTI' = TTI + T0 (4)
Jeśli otwór ma wyraźne nachylenie, należy podczas obliczeń TTI sprowadzić odcinek d do pionu, czyli obliczyć wartość d' = d*cosg ,
gdzie g jest kątem odchylenia osi otworu od pionu i do obliczeń we wzorze (1) przyjąć d’ zamiast d.
Porównanie czasu sumarycznego z PA - TTI’ i czasu TVC z PPS
6. obliczenie dryftu: D(ni*d) = TVC - TTI' (5)
7. odjęcie od dryftu D wartości T0: D' = D - T0 (6)
Dryft D’ przyjmuje wartości dodatnie i ujemne.
8. obliczenie korekty wartości DT na podstawie algorytmu 1, DTcorr1,
9. obliczenie korekty wartości DT na podstawie algorytmu 2, DTcorr2,
Algorytm1 obliczenia skorygowanej krzywej DTcorr1 , preferowany dla dryftu dodatniego,
11. analiza trendu zmiany dryftu D’: zaczynamy od pierwszej głębokości, HPA i sprawdzamy,
jakie jest nachylenie krzywej dryftu w kolejnych punktach, wybieramy interwały o stałym (w przybliżeniu) nachyleniu o długości ok. 100 - 200 m,
klikając głębokości na ekranie?
12. tworzymy dodatkowe profilowanie ZD(ni*d), które informuje interpretatora, jaki jest trend na danej głębokości;
jeśli nachylenie krzywej dryftu w wyznaczonym odcinku jest dodatnie, D'((ni+1)*d) jest większe od D'((ni)*d) , profilowanie ZD przyjmuje wartość dodatnią, (+), jeśli nachylenie krzywej
dryftu w wyznaczonym odcinku jest ujemne, D'((ni+1)*d) jest mniejsze od D'((ni)*d) , profilowanie ZD przyjmuje wartość ujemną, (-),
13. odczytujemy wartości dryftu D’ w górnym punkcie, D’1 (H), zaznaczonym na wykresie) i w dolnym punkcie, D’2 (H2) i obliczamy różnicę: dDT21 = D2' - D1' (8)
14. obliczamy poprawkę do wartości DT: (9)
15. obliczamy poprawioną wartość DTcorr1 w interwale H21:
(10)
Algorytm2 obliczenia skorygowanej krzywej DTcorr1 , preferowany dla dryftu ujemnego,
16.obliczamy poprawioną wartość DTcorr2=G*(DT-DTmin)+DTmin (11)
gdzie: DTmin jest to wybrana przez interpretatora wartość DT, poniżej której nie stosuje się poprawki,
G=dDT21/T
T jest wartością całki w granicach H1 i H2 z wyrażenia dDT(h), H1 i H2 są głębokościami wyznaczonymi na wykresie dryftu,
dDT(h) = DT - DTmin, h = d*n
17. obliczamy trzecią wartość skorygowanego czasu interwałowego, DTcorr3
DTCORR3 = DTCORR1 dla wartości ZD dodatnich,
DTCORR3 = DTCORR2 dla wartości ZD ujemnych.