W związku ze zmniejszającymi się zasobami gazu w złożach zlokalizowanych w utworach piaskowcowych i skałach węglanowych coraz większym zainteresowaniem cieszą się obiekty strukturalne związane z utworami piaskowcowo-mułowcowo-łupkowymi ulokowane w utworach miocenu. Złoża te charakteryzują się dużą zmiennością facjalną tych utworów, w tym wyklinowaniem się lamin i warstewek. W związku z powyższym porowatość oznaczona dla tych skał nie może być liczona do bilansu porowatości efektywnej, służącej do oszacowania zasobów wydobywalnych. Obliczając miąższość efektywną, powinno się pominąć miąższość skał łupkowych. Jednym z zasadniczych problemów, będących przedmiotem tego opracowania, jest wydzielenie miąższości skał należących do poszczególnych typów litologicznych. W celu rozpoznania zmienności litologicznej tego typu skał użyto pomiarów profilowania gamma. Należy jednak podkreślić, że ze względu na bardzo małą miąższość warstw piaskowcowych i mułowcowych nie są one możliwe do jednoznacznej identyfikacji za pomocą profilowania gamma. Profilowanie gamma wskazuje raczej wartości uśrednione w obrębie pomiaru niż konkretne warstwy. W związku z tym dokonano oszacowania procentowego udziału miąższości każdego z typów litologicznych w całkowitej miąższości struktury przewierconej przez każdy odwiert. Podział taki przeprowadzono za pomocą analizy rozkładów empirycznych pomiarów gamma. Założono, że niejednorodności litologiczne uwidocznią się na wykresach histogramów w formie rozkładów wielomodalnych lub wyraźnych zmian monotoniczności. Dane te posłużyły do wygenerowania przestrzennych rozkładów parametrów na siatce modelu w punktach nieobjętych pomiarem. Użyto do tego celu metody krigingu zwyczajnego oraz krigingu z trendem. Przy zastosowaniu tych samych metod interpolacyjnych wygenerowano mapy stropu i spągu modelu badanego złoża. Korzystając z wyników pomiarów porowatości, dokonano przyporządkowania wyniku tych pomiarów do wydzielonych wcześniej typów litologicznych. Dla każdego typu skały i przy użyciu obydwu metod krigingu wyznaczono rozkład porowatości w punktach siatki modelu złoża. Do oszacowania niepewności uzyskanych wyników posłużono się metodą kroswalidacji. Na zakończenie obliczono objętość porów efektywnych modelu złoża, zakładając, że tylko skały piaskowcowe i mułowcowe są źródłem porów efektywnych. Oszacowano niepewność uzyskanego wyniku. Na rzecz opracowania posłużono się rzeczywistymi danymi z pomiarów geofizycznych dla złoża w utworach miocenu.
EN
Due to decreasing gas reserves in deposits found in sandstone f and carbonate formations, geological structures composed of mudstone and shale layers found in Miocene formations are growing in popularity. These deposits are characterized by high facial variability. Shale rocks have very low permeability. Therefore, the porosity determined for these rocks cannot be taken into account in the balance of the effective porosity used to estimate the natural gas reserves. Shale rock thickness should be excluded when calculating effective thickness of the gas bearing formations. The main problem of this study is the proper separation of the rocks belonging to individual lithological facies. Gamma ray logs were used to identify lithological variability of this type of rock. It should be emphasized, however, that due to the very low thickness of sandstones and mudstones, it was not possible to identify them clearly using the archival gamma ray logs. The archival gamma ray logs indicate the average values of the layer rather than that of specific laminas, because of its measurement resolution. Therefore, an estimation was made in order to calculate the percentage share of the thickness of each lithological type in the total layer thickness in each well. This division was made using the analysis of empirical distributions of gamma ray logs. It was assumed that lithological heterogeneities would be visible on experimental histograms in the form of multimodal distributions or clear changes in monotonicity. These data were used to generate a spatial distribution of parameters on the model grid at points not covered by the measurement. For this purpose, the ordinary kriging method and the kriging with trend method were used. Using the same interpolation methods, the structure maps of the top and bottom of the model were generated. The results of the porosity measurements were assigned to the previously separated lithological types of rocks. The distribution of porosity at the grid points of the deposit model was determined for each type of rock and using both kriging methods. The cross-validation method was used to assess the uncertainty of the results. Finally, the effective pore volume of the deposit model was calculated, assuming that only sandstone and mudstone rocks are sources of the effective pores. The uncertainty of the analysis was estimated. Real data from geophysical measurements for the Miocene gas field were used for the study.
JavaScript jest wyłączony w Twojej przeglądarce internetowej. Włącz go, a następnie odśwież stronę, aby móc w pełni z niej korzystać.