Обучение по природни науки и върхови технологии

От изследователските лаборатории

A STUDY ON HOW PERMEABILITY HETEROGENEITY INFLUENCES INITIAL WATER CUT IN UNDERSATURATED OIL FORMATIONS

https://doi.org/10.53656/nat2024-5-6.03

Резюме. The exploration of undersaturated oil-bearing, layer-by-layer heterogeneous oil reservoirs has become increasingly important due to the overall depletion of oil reserves. It is associated with non-standard changes in well production and increased values of water cut when these hard-to-recover oil reserves are developed. In its classical understanding, an undersaturated deposit is a conventionally large transitional water-oil zone of the reservoir. The initial water cut is generally expected to correspond to the high initial water saturation of the formation. However, in some instances, this may not be the case. Predicting the trajectory of water cut during the initial stages of oil field development presents a considerable challenge, ascertaining whether it will increase, stabilize, or decrease. Hence, the primary objective of this study is to identify the key parameter and analyze the range of its variation using statistical methods. This analysis will enable us to make highly accurate predictions regarding the anticipated variations in water cut during the feasibility assessment of new well drilling projects. In this study, the graphical dependence of the deviation of water cut during the early stages of well operation on the layer-by-layer heterogeneity of the formation was found for the conditions of the selected object, the Vatyogan oil field. The study's findings strongly support the development of innovative technologies for selectively isolating reservoirs. These technologies would enable the reversible restriction of permeability in the low-permeability, low-saturated part of the reservoir. This is a significant advancement considering that current technologies primarily target selective isolation of highly permeable and water-saturated interlayers.

Ключови думи: permeability heterogeneity; water cut; undersaturated oil formation; water oil zone, OWC

1. Introduction

The initial geological and physical properties of reservoirs and the fluids saturating them largely determine the future performance of oil-producing wells and the oil production facility (Petrakov, Kupavykh & Kupavykh 2020; Yu et al. 2022; Agzamov, et al. 2023; Al-Obaidi & Chang 2024). However, some processes in the bottom hole formation zone of wells that have opened young, unformed oil deposits characterized by low initial oil saturation are not so obvious.

The prevailing assumption is that the initial fluid present in the rock of the prospective oil deposit was water. Microscopic accumulations of hydrocarbons in much smaller volumes after their formation were evenly distributed throughout the entire thickness of sedimentary rocks, the size of which was many times greater than the thickness of the prospective (future) oil reservoir (Zou et al. 2013; He et al. 2023; Zhang et al. 2022; Al-Obaidi, Patkin & Guliaeva 2003). During the prolonged process of oil deposit formation, oil accumulates in geological traps, which are impermeable rock layers shaped like an anticline. The oil rises upwards due to the difference in density between oil and water. This movement happens rapidly in large pores and cracks due to the dominance of gravitational forces over capillary forces. The percentage of oil in the trap increases, and water saturation decreases accordingly. The rate of oil saturation growth over time is uneven, as it is uneven for reservoirs of different permeability and depths. Regardless of the deposit's formation time, hydrocarbons cannot completely replace water (capillary trapping of water in small dead-end pores) (Pentland et al. 2012; Lawrence 2019; Mccarter et al. 2020; Khalaf & Al-Obaidi 2019).

The oil field under study is Vatyogan, a large oil field located in Russia, specifically in the Khanty-Mansi Autonomous Okrug (Vaguet 2013). The proven reserves of the field as of the end of 2004 amount to almost 1.5 billion barrels of oil. Cumulative oil production has exceeded 135 million tons. The deposit is of the layer-arch type and is multilayered. The reservoir is of the terrigenous type. Sandy-clayey deposits of the Quaternary, Neogene, Paleogene, Cretaceous, and Jurassic ages represent the lithological and stratigraphic section of the Vatyogan field (Sharma et al. 2000; Chang, AlObaidi & Smirnov 2023).

For example, the oil deposits of the Vatyogan field are characterized by a high initial oil saturation of 98% (formed oil deposit). Low oil saturation can be found only in the transition zones of the formation, in areas where water and oil zones of deposits come into contact (OWC). In this case, the transition zones are relatively small in thickness (conditionally 0.5 – 1 m). In wells that have penetrated such a formation, a prolonged period without water is a frequent event (the transition zone is typically left unperforated) (Korolev et al. 2023; Wang et al. 2023; Melzer, Koperna & Kuuskraa 2006; Al-Obaidi, Smirnov & Alwan 2021).

The well can produce more than 70-80% of recoverable oil reserves with no more than a 5% water cut. In the drainage area of the well, stagnant, residual waters do not affect the extraction of primary oil reserves.

The primary cause of water flooding is typically water from oil-bearing formations, injection wells, formation water, and breakthroughs due to behind-the-casing flows (Haines et al. 2024; Gambolati & Teatini 2015; Muggeridge et al. 2013; Al-Obaidi 2016). In the first two cases, the water cut of well production increases to 60-70% (breakthrough of remote waters) and stabilizes for production from 5 to 15% of recoverable oil reserves (Rajbongshi & Gogoi 2024; Xue, Liu & Zhang 2023; Jia et al. 2023; Al-Obaidi, Chang & Khalaf 2024). Water is produced from the highly permeable part of the reservoir, but oil continues to flow from the oil-bearing part. In the end, all layers are flooded with water. This process involves a series of water-isolation operations along the well, which gradually increases the coverage of the formation by allowing drainage in the still oil-saturated part of the formation over time (Zhao & Gates 2015; Lee et al. 2023; Katende & Sagala, 2019; Chang, Al-Obaidi & Patkin 2021).

The study examined how permeability heterogeneity within oil reservoirs affects the initial water cut and its subsequent evolution. These factors are crucial in evaluating the viability of new well-drilling projects.

2. Methodology and materials

For oil production objects (facilities) , it is scarce to observe a prolonged initial water-free period in layers characterized by low initial oil saturation ranging from 35 to 70%.These are the so-called young oil production objects with unformed deposits, in which mass exchange processes are still intensive and are far from the final slowdown. The entire oil deposit is one large transition zone (Male 2019; Long 2023). Man-made human intervention (well drilling, creation of differential pressure in the formation, formation of a flow field in the well drainage area) at this stage of the deposit development is accompanied by processes that are not typical for formed deposits with an initial oil saturation of more than 90% and a limited transition zone in the OWC.

Figure 1 depicts a section of the initial oil-saturated formations map for the Vatyogan field, specifically the Aptian stage, AB1-3 formation (southwest), showing accumulated oil and liquid extractions and water injection. This object (AB1-3 formation) is characterized by low (from 33 to 69%) initial oil saturation of productive formations (Fu, Wang & Xiao 2024; Fanchi, Christiansen & Heymans 2002; Guo et al. 2021; Miel, Hameed & Hussein 2022).

Figure 1. Map of initial oil-saturated formations with accumulated oil and liquid extractions and water injection. Vatyogan field, Aptian stage, AB1-3 reservoir (southwest)

Wells with low initial oil saturation are often associated with high initial water cut due to the phase permeability of fluids for a particular reservoir type (see Fig. 2). However, it is not possible to immediately determine whether the water cut will increase, remain stable, or decrease over time. However, this issue is crucial in the feasibility study for efficiently developing such objects during the drilling stage, where it is necessary to know the starting and typical characteristics of future wells (Wang et al. 2018; Rudyk et al. 2024; Zhang et al. 2019; Smirnov & Al-Obaidi 2008).

Figure 2. Distribution of initial water cut of wells (first month) and the average initial oil saturation of formations. Vatyogan field, Aptian stage, formation AB1-3 (southwest)

Given the substantial impact of layer-by-layer heterogeneity on the formation of strata, it is reasonable to expect that this factor will also influence the sequence of reservoir saturation during deposit formation (see Fig. 3). As the layers of different permeability in an oil reservoir gradually reach almost equal levels of oil saturation, approaching maximum levels of 97-98%, the permeability plays a crucial role in determining how the oil is displaced. This, in turn, influences the average water cut of the well production during the initial period.

Figure 3. Dependence of the initial oil saturation of interlayers on their permeability within the section of one well. Vatyogan field, Aptian stage, AB1-3 formation (southwest).

In most cases, formation water and injected water have greater mobility than oil and tend to break through highly permeable interlayers, while low-permeability oil-saturated interlayers are gradually deactivated. However, this phenomenon is primarily observed during the late development stage and in the formed deposits.

In the case of the object under study, the deposit's formation process led to advanced saturation (maturation), especially in the highly permeable section of the reservoir (refer to Figure 3). This highlights the heterogeneity of the reservoir's permeability across the section. Next, a study was conducted on the influence of permeability heterogeneity on the initial water cut and its subsequent change, which are the indicators of our interest.

3. Results and discussions

Figures 4 – 9 show the performance indicators of six production wells of the Vatyogan field, operating in the deposit section shown in Figure 1. The initial indicators and the dynamics of water cut in the production of these wells are of interest to our study.The wells were chosen to showcase the wide range of variations in the nature of the initial water flooding. The first three cases (refer to Fig. 4 – 6) are uncommon, while the following three cases are much more typical (refer to Fig. 7 – 9). From the graphs (see Fig. 4-9) it is clear that after the wells are put into operation, after about six months their water cut changes, but not always in a positive direction.

Figure 4. Indicators of well number 3522 of the Vatyogan field

Figure 5. Indicators of well number 1437 of the Vatyogan field

Figure 6. Indicators of well number 1438Y of the Vatyogan field

Figure 7. Indicators of well number 1435 of the Vatyogan field

Figure 8. Indicators of well number 1436 of the Vatyogan field

Figure 9. Indicators of well number 1377 of the Vatyogan field

Table 1 displays the main geological and physical characteristics of the productive formations of the wells in the area under consideration. It also illustrates the changes in their initial water cut under low initial oil saturation conditions of the AV1-3 reservoir of the Vatyogan field.

Table 1. Comparison of geological and physical characteristics of heterogeneous productive formations of wells and changes in their initial water cut under conditions of low initial oil saturation of the reservoir AV1-3 of the Vatyogan field

Well NoOil formationthickness, mPorosity,%Oil saturation, %Permeability, μm2Layer-by-layerheterogeneity, rel. unitInitial water cut, %Water cut in the rst sixmonths of operation, %Deviation of water cutfrom the initial level inthe initial period, %518315,721,758,80,08524,20622,15,31‒16,8137620,422,056,00,15897,65820,95,9‒15,0137518,322,958,10,19303,9089,92,9‒7,0144013,221,356,70,10082,8755,62,3‒3,3131615,120,247,30,04854,52114,46,4‒8,0143510,420,857,20,04945,89844,624,7‒19,9518219,323,059,40,55248,8515633‒23,0137718,322,064,60,10323,585030,3‒19,7143916,822,360,80,13504,6482,11,3‒0,814368,420,655,20,06651,52339,929,4‒10,5137810,823,961,10,47562,169,77,18‒2,5131516,922,163,30,47258,2816045,2‒14,8518018,623,754,50,42154,738,56,6‒1,9209012,220,556,10,03673,10575,284,18,935227,421,145,80,08661,72656,468,812,4131719,119,954,20,03652,29135,146,911,851786,321,651,50,08471,6859,116,67,51438Y1821,960,30,10262,3462,963,1143714,921,850,00,06901,6575,929,623,7517916,222,151,60,07780,7862,31310,7518118,923,262,20,28342,61811312,0

The table compares Oil-saturated formations' thicknesses, porosity, oil saturation, permeability, layer-by-layer heterogeneity compared with the initial water content, the water cut in the first six months of operation, and the deviation of the water cut from the initial level in the initial period.

Upon pairwise examination of the parameters, it was found that the initial water cut's graphical dependence on the layer-by-layer heterogeneity of the formation is the most significant (refer to Fig. 10). The deviation of the initial water cut is the deviation from the initial water cut level during the initial period (Yang 2024; Wang et al. 2013).

Figure 10. Estimation of the dependence of the deviation of the initial water cut in the initial period on the layer-by-layer heterogeneity for undersaturated oil reservoirs

Three distinct zones (I, II, and III) have been identified within the given construction. These zones reflect varying degrees of layer-by-layer heterogeneity in the formation. Zone I corresponds to the range of 0 to 2, Zone II spans from 2 to 4, and Zone III includes values of 4 or more relative units (see Fig. 10).These zones were selected according to the delimitation principle of areas with stable (in most cases) facts of positive and negative deviations of water cut. The third zone represents both positive and negative values between these zones. The absolute values of the zone boundaries (layer-by-layer heterogeneity) have been rounded to whole numbers. Further investigation involving a larger number of wells will help us to determine the exact value of these boundaries.

4. Conclusions

The research presented in this work was initiated due to unusual initial water cut indicators in several wells of the Vatyogan field of the AV1-3 formation.The graphical dependence of the deviation of the initial water cut on the layer-by-layer heterogeneity of the formation for undersaturated oil reservoirs is the result of the pairwise comparative analysis.

Upon delving into the data from our research, a compelling trend was observed. It turns out that the greater the layer-by-layer heterogeneity of the reservoir formation in terms of permeability (specifically, when the square of the variation coefficient exceeds 4 relative units), the more advantageous it becomes for drilling in that particular section. This insight could have significant implications for drilling operations.

When drilling relatively homogeneous formations (heterogeneity less than 2 relative units), rapid water cut up to 75-85% is predicted with high probability. This is typical for the conditions of the AV1-3 formation of the Vatyogan field. Moreover, it is difficult to make an unambiguous prediction of the initial water cut within the heterogeneity range of 2 to 4 relative units.

Selective reservoir isolation technologies aimed at isolating the waterflooded, highly permeable part of the reservoir are ineffective in such conditions. New solutions are needed that allow reversible limitation of the permeability of the lowpermeability, low-saturated part of the reservoir.

REFERENCES

AGZAMOV, A. K. et al., 2023. Analysis of the Degree of Influence of Geological and Technological Factors on The Initial Flow Rates of Wells Commissioned After Their Downtime, E3S Web Conf., vol. 411. 01023. https://doi.org/10.1051/e3sconf/202341101023.

AL-OBAIDI, S., 2016. High Oil Recovery Using Traditional Waterflooding Under Compliance of The Planned Development Mode. Journal of Petroleum Engineering & Technology, vol. 6, no. 2, pp. 48 – 53.

AL-OBAIDI, S. & CHANG, W. J., 2024. Techniques to Boost Oil Production in the Development of Multi-Reservoir Fields. Journal of Earth & Environmental Waste Management, vol. 2, no. 3, pp. 1 – 7.

AL-OBAIDI, S., CHANG, W. J., & KHALAF, F. H., 2024. Offshore Heavy Oil Displacement Using Water Flooding: Flow Characteristics and Efficiency. Natural Sciences and Advanced Technology Education, vol. 33, no. 1 – 2, pp. 91 – 110. https://doi.org/10.53656/nat2024-1-2.01.

AL-OBAIDI, S., PATKIN, A. A. & GULIAEVA, N. I., 2003. Advance Use for The NMR Relaxometry To Investigate Reservoir Rocks. JoPET. vol. 2, no. 3, pp. 45 – 48.

AL-OBAIDI, S, SMIRNOV, V. & ALWAN, H. H, 2021. Experimental Study About Water Saturation Influence on Changes in Reservoirs Petrophysical Properties. Walailak J Sci Tech, vol. 18, no. 13, pp. 1 – 10.

CATHLES, L. 2019. On the Processes that Produce Hydrocarbon and Mineral Resources in Sedimentary Basins. Geosciences, vol. 9, no. 12, p. 520. https://doi.org/10.3390/geosciences9120520.

CHANG, W. J, AL-OBAIDI, S. & PATKIN, A. A., 2021. The Use of Oil soluble Polymers to Enhance Oil Recovery in Hard to Recover Hydrocarbons Reserves. Int Res J Mod Eng Technol Sci., vol. 3, no. 1, pp. 982 – 987.

CHANG, W. J., AL-OBAIDI, S. & SMIRNOV, V. I., 2023. Application of Multi-Stage Hydraulic Fracturing (MSFH) in the Development of Deep-Lying Deposits with Low Permeability of Hydrocarbon Fields. J Geol Geophys. vol. 12, p. 1154.

FANCHI, J. R., CHRISTIANSEN, R. L., & HEYMANS, M. J., 2002. Estimating Oil Reserves of Fields with Oil/Water Transition Zones. SPE Res Eval & Eng, vol. 5, pp. 311 – 316. https://doi.org/10.2118/79210-PA.

FU, J., WANG, G. & XIAO, Z., 2024. Petrophysical Analysis of Low-Contrast Oil Reservoirs with Low Porosity and Low Formation-Water Salinity in Sapu Intercalation of Longxi Area, Daqing Oilfield, China. Applied Sciences, vol. 14, no. 8, p. 3210. https://doi.org/10.3390/app14083210.

GAMBOLATI, G., & TEATINI, P., 2015. Geomechanics Of Subsurface Water Withdrawal and Injection. Water Resources Research, vol. 51, no. 6, pp. 3922 – 3955. https://doi.org/10.1002/2014WR016841.

GUO, M., et al., 2021. Rapid Evaluation of Probable Recoverable Reserves in Tight Reservoirs: A Case Study of Chang 8 Reservoir (Eighth Member of Yanchang Formation) in Honghe Oil Field. Petroleum Geology & Experiment, vol. 1, pp. 154 – 160.

HAINES, S. S., VARELA, B. A., TENNYSON, M. E. et al., 2024. How Quickly Do Oil And Gas Wells “Water Out”? Quantifying And Contrasting Water Production Trends. Nat Resour Res, vol. 33, pp. 591 – 608. https://doi.org/10.1007/s11053-024-10308-6.

HE, P., et al., 2023. The Relationship Between Fluid Evolution and Hydrocarbon Accumulation and Metallization in The NanpanjiangYoujiang Basin: Evidence from Calcite Petrography and Fluid Inclusions, Energies., vol 16, no. 1, p. 328. https://doi.org/10.3390/ en16010328.

JIA, D., et al., 2023. Collaboration Between Oil Development and Water/Power Consumption in High-Water-Cut Oilfields. Sustainability, vol. 15, no. 14, p. 11405. https://doi.org/10.3390/su151411405.

KOROLEV, E. et al., 2023. Oil Deposits of The Visean Stage of The Northeastern Part of The Tatarstan Republic. BIO Web Conf., vol. 71, p. 02037. https://doi.org/10.1051/bioconf/20237102037.

KATENDE, A. & SAGALA, F., 2019. A Critical Review of Low Salinity Water Flooding: Mechanism, Laboratory and Field Application, Journal of Molecular Liquids, vol. 278, pp. 627 – 649. https://doi.org/10.1016/j. molliq.2019.01.037.

KHALAF, F. H. & AL-OBAIDI, S., 2019. Development of Traditional Water Flooding to Increase Oil Recovery. Int J Sci Tech Res., vol. 8, pp. 177 – 181.

LEE, P. et al., 2023. Laboratory Evaluation on The Water Flooding Characteristics in Bottom Water Reservoir Containing Interbeds. ACS Omega, vol. 8, no. 45, pp. 42409 – 42416. https://doi.org/10.1021/ acsomega.3c04873.

LONG, X., 2023. Adjustment Of Infill in Transition Zone to Improve Oilfield Productivity, E3S Web of Conferences, vol. 375, 01024. https://doi. org/10.1051/e3sconf/202337501024.

MALE, F., 2019. Using A Segregated Flow Model to Forecast Production of Oil, Gas, And Water in Shale Oil Plays. Journal of Petroleum Science and Engineering, vol. 180, pp. 48 – 61. https://doi.org/10.1016/j. petrol.2019.05.010.

MCCARTER, C. P. R., et al., 2020. Pore-Scale Controls on Hydrological and Geochemical Processes in Peat: Implications on Interacting Processes. Earth-Science Reviews, vol. 207, p. 103227. https://doi.org/10.1016/j. earscirev.2020.103227.

MELZER, L. S., KOPERNA, G. J., & KUUSKRAA, V. A., 2006. The Origin and Resource Potential of Residual Oil Zones. In: SPE Annual Technical Conference and Exhibition. https://doi. org/10.2118/102964-MS.

MIEL, H., HAMEED, A. O. & HUSSEIN, K. F., 2022. Modeling And Monitoring the Development of An Oil Field Under Conditions of Mass Hydraulic Fracturing. Trend Sci., vol. 19, no. 8, p. 3436.

MUGGERIDGE, A., et al., 2013. Recovery Rates, Enhanced Oil Recovery and Technological Limits. Philos. Trans. A Math. Phys. Eng. Sci., vol. 372, no. 2006. https://doi.org/10.1098/rsta.2012.0320.

PENTLAND, C. H. et al., 2012. The Influence of Pore Space Geometry on The Entrapment of Carbon Dioxide by Capillary Forces. In: SPE Asia Pacific Oil and Gas Conference and Exhibition. https://doi. org/10.2118/158516-MS.

PETRAKOV, D., KUPAVYKH, K. & KUPAVYKH, A., 2020. The Effect of Fluid Saturation on The Elastic-Plastic Properties of Oil Reservoir Rocks. Curved and Layered Structures. vol. 7, pp. 29 – 34. https://doi. org/10.1515/cls-2020-0003.

RAJBONGSHI, A. & GOGOI, S. R., 2024. A Review on Oilfield Produced Water and Its Treatment Technologies. Petroleum Research. https://doi. org/10.1016/j.ptlrs.2024.06.003.

RUDYK, S. et al., 2024. Matching Of Water Breakthroughs in A LowResistivity Oil Reservoir Using Permeability Anisotropy. Applied Sciences, vol. 14, no. 11, p. 4618. https://doi.org/10.3390/app14114618.

SHARMA, R. et al., 2000. Integrated Study of Vatyogan Oil Field, WestSiberian Region of Russia. In: SPE Annual Technical Conference and Exhibition. https://doi.org/10.2118/62900-MS.

SMIRNOV, V. I. & AL-OBAIDI, S., 2008. Innovative Methods of Enhanced Oil Recovery. Oil Gas Res., vol. 1, no. 1.

VAGUET, Y., 2013. Oil and Gas Towns in Western Siberia: Past, Present and Future Challenges. Proceedings from the First International Conference on Urbanisation in the Arctic, vol. 6, no. 218, pp. 125 – 132.

WANG, D.-Y. et al., 2023. Petroleum Geochemistry and Origin of Shallow-Buried Saline Lacustrine Oils in The Slope Zone of The Mahu Sag, Junggar Basin, NW China. Petroleum Science, vol. 20, no. 6, pp. 3363 – 3378. https://doi.org/10.1016/j.petsci.2023.08.024.

WANG, J. et al., 2018. Characteristics and Influential Factors of Water Cut Change in Low Permeability Water-Drive Oil Reservoirs in Shengli Oilfield. Oil and Gas Geology, vol. 39, pp. 833 – 838.

WANG, R. R. et al., 2013. A New Water Displacement Curve for The High Water-Cut Stage. Petroleum Science and Technology, vol. 31, no. 13, pp. 1327 – 1334. https://doi.org/10.1080/10916466.2012.756015.

XUE, L., LIU, P. & ZHANG, Y., 2023. Status and Prospect of Improved Oil Recovery Technology of High Water Cut Reservoirs. Water, vol. 15, no. 7, p. 1342. https://doi.org/10.3390/w15071342.

YANG, R., 2024. Further Study on Oil/Water Relative Permeability Ratio Model and Waterflooding Performance Prediction Model for High Water Cut Oilfields Sustainable Development. J Petrol Explor Prod Technol, vol. 14, pp. 1297 – 1309. https://doi.org/10.1007/s13202-024-01753-3.

YU, Z., et al., 2022. Analysis of Factors of Productivity of Tight Conglomerate Reservoirs Based on Random Forest Algorithm. ACS Omega, vol. 7, no. 23, pp. 20390 – 20404. https://doi.org/10.1021/acsomega.2c02546.

ZHANG, Y., et al., 2022. Hydrocarbon Accumulation Conditions and Favorable Exploration Zones in The Lunpola Basin, Tibet. Frontiers in Earth Science, vol. 10. https://doi.org/10.3389/feart.2022.1034069.

ZHANG, D. et al., 2019. The Residual Potential Of Bottom Water Reservoir Based Upon Genetic Algorithm For The Relative Permeability Inversion. Journal of Geoscience and Environment Protection, vol. 7, pp. 192 – 201. https://doi.org/10.4236/gep.2019.74012.

ZHAO, D. W. & GATES, I. D., 2015. On Hot Water Flooding Strategies for Thin Heavy Oil Reservoirs. Fuel, vol. 153, pp. 559 – 568. https://doi. org/10.1016/j.fuel.2015.03.024.

ZOU, C. H. et al., 2013. Continuous Hydrocarbon Accumulation Over a Large Area as a Distinguishing Characteristic of Unconventional Petroleum: The Ordos Basin, North-Central China, Earth Science Reviews, vol. 126, pp. 358 – 369. https://doi.org/10.1016/j.earscirev.2013.08.006.

Година XXXIII, 2024/5-6 Архив

стр. 347 - 361 Изтегли PDF