Introduction

The clay minerals are hydrated aluminosilicates categorized as phyllosilicates, characterized by two-dimensional infinite sheets. They denote substances primarily comprised of fine particles with a size smaller than 2µ1. The presence of clay minerals in hydrocarbon reservoirs poses a significant challenge to the extraction of hydrocarbon. These kinds of minerals that are mostly find in the clastic formations, playing a crucial role in the governing reservoir characteristics, including porosity, permeability, capillary pressure, density, natural radioactivity, wettability, and water saturation. Nevertheless, the capacity to generate energy from these reserves remains considerably below their potential capacity due to the role of clay minerals. A key reason behind this problem lies in the prevalence of clay minerals within the reservoir rock. These minerals, prevalent in the fine-grained sediments, gas hydrate-bearing sediments, sandstone, mudstone, and shale, notably impact the reservoir quality by introducing numerous micropores into the existing intra-granular pore space2. Comprehensive understanding of clay minerals is vital in the evaluation process of hydrocarbon reservoirs. Even in small quantities, clay minerals play a significant influence on petroleum reservoir properties and can impact the recordings from various well logging tools3,4,5,6,7. Hence, an extensive grasp of clay mineralogy becomes crucial for both reservoir modelling and logging data interpretation2,8,9.

Clay minerals, include kaolinite, chlorite, illite, and illite–smectite mixed layer (I/S), are remarkably complex, often intergrown and encompassing a wide array of solid solutions formed through diverse processes. They arise across varying pressure, temperature, rock, and fluid compositions10,11. Identifying clay minerals precisely poses challenges because of the natural presence of interlayered phases, where different clay species e.g., illite and smectite (I/S) coexist in the layers. The issue arises from overlapping peaks of these constituent minerals solely and the similarity in their diffraction patterns when not treated with ethylene glycol. For instance, illite and I/S clays significantly impair reservoirs by diminishing pore-throat size and connectivity10. This similarity complicates identification, particularly in cases of low clay contents12,13.

Illite resulting from the transformation of potassium feldspars or muscovite, possesses a 2:1 structure1. Illite can have distinct shapes based on the levels of potassium saturation during its formation. When there is a high level of potassium saturation, hexagonal illite tends to precipitate, while lower saturation leads to the formation of lath-shaped illite. The presence of fibrous illite particularly hinders fluid flow in sandstones. Consequently, the primary morphology of illite significantly impacts reservoir permeability, underscoring the importance of finding its distribution and morphology for predicting sandstone reservoir properties14. Chlorite, a three-layer non-swelling clay mineral, plays a crucial role in regulating reservoir quality in various siliciclastic rocks. In deeply buried reservoirs, chlorite coating on sand grains inhibits the formation of quartz cements, resulting in unexpectedly favorable reservoir quality. Conversely, an abundance of chlorite, as seen in clay-rich siltstone and sandstone, causes pore throats to become obstructed, leading to extremely low permeability15. Taking into account the influence of clay mineralogy on a clastic reservoir, illite, chlorite, and kaolinite have the potential to significantly impact reservoir characteristics, such as porosity7.

Montmorillonite remains a significant and extensively utilized clay mineral in various industries. It falls under the category of di-octahedral smectites and originates from the chemical modification of volcanic ashes due to weathering. Smectites represent phyllosilicates of the 2:1 type, encompassing minerals, such as montmorillonite, bentonite, nontronite, saponite among others. They exhibit diverse compositions of Na, Al, Fe, and Mg. Kaolinite, on the other hand, arises from the transformation of aluminum-rich silicates like feldspars and muscovite1. Reservoir quality primarily depends on lithological characteristics rather than chlorite coatings. Kaolinite in reservoir is mostly formed from the chemical alteration of feldspar, chlorite emerges as clay coatings transformed from clay precursors due to the flocculation of suspended detrital clays or crystallization of pore fluids, while illite forms from the illitization of smectite group, K-feldspars, and authigenic kaolinite10.

The unique features of clay minerals, such as abundance, generally low specific gravity, low electrical resistivity, and swelling problems make them highly important for hydrocarbon resource exploration. Research on these properties is conducted globally and applied in various fields, including identifying clay cap rocks or cap rocks in geothermal systems, improving understanding of geothermal reservoirs, and tracking their changes over time. In certain geothermal systems, they can reduce permeability, necessitating thorough investigation16.

Alansari et al.17 explored the significance of considering the microporosity of clay minerals in the formation evaluation stage. They advocated for an integrated approach involving core analysis, well logs, and petrography by SEM and thin section studies, rather than relying solely on well logs. Griffiths et al.18, Sun et al.10, and Wang et al.2 highlighted the characteristics, content, types, distribution, and impacts of clay minerals on reservoir properties. The layered arrangement of clay minerals can be essential for comprehending the formation of methane hydrate in marine environments which was a topic of interest for Li et al.19 in a research about \({\text{CH}}_{4}\) hydrate formation in kaolinite particles with different surface contacts. They stated that molecular dynamics simulations revealed that the nucleation of methane hydrate is significantly influenced by the surface interactions of kaolinite particles which play a crucial role in the formation of aggregates within natural sediments. Jiu et al.6 examined the influence of clay minerals and cementation on pore throat in tight sandstone gas reservoirs. Zhou et al.20 investigated two critical parameters in the assessment of reservoir quality, including porosity and permeability changes following water injection in low permeability sandstone reservoirs. Clays can pose serious challenges for well logging operations. Jozanikohan and Abarghooei7 documented studies on clay mineralogy in clastic reservoirs using Fourier transform infrared spectroscopy (FTIR). Hou et al.21 explored the effects of clay mineral types and content on the effectiveness of lacustrine organic-rich shale reservoirs. Belghazdis and Hachem1 investigated the impact of diverse production climates and clay sources on clay structures and types. Deon et al.13 employed XRD, SEM analyses, and reflectance spectroscopy to detect interlayered illite/smectite clay minerals. The presence of clay minerals in terms of quantity, distribution, and structure significantly influences the most of physical characteristics of reservoirs, such as porosity, permeability, density, natural radioactivity, electrical conductivity, water content in petroleum fields. Studies indicate that the porosity of sandstone reservoirs is remarkably impacted by the clay mineral13,22. Aigbadon et al.23 highlighted the dominance of calcite in formations from the southern Benue Trough and the Dahomey Basin in Nigeria, suggesting their potential as limestone reservoirs with quartz and chlorite indicating siliciclastic input and diagenetic alterations. The study provided valuable insights into the depositional environments, diagenetic processes, and the role of elements in the evolution of carbonate deposits. The primary mineral composition of shale lithofacies associations can be accurately predicted using neural network analysis, whole-rock X-ray diffraction data, and wire-line log data, offering a foundation for lithofacies identification in shale intervals at well locations without core and test data24. The surface groups of clay minerals play a crucial role in CO2 hydrate nucleation by altering the distribution of water molecules and CO2. This, in turn, affects CO2 storage in deposits and has significant impacts on geochemical processes, including the interfacial interactions between clay surfaces and water molecules25. Conventionally, estimating two critical reservoir parameters, including volume of shale and porosity involves a very high error rate which led to the use of the multilayer perceptron (MLP) artificial neural network (ANN) to minimize this error. Implementing MLP proved that the error percentages for volume of shale and porosity significantly reduced from 89.46% and 58.3% in the traditional petrophysical approaches to1.28701 E−9 and 2.78731 E−4, respectively. The validation of the research findings from the MLP method considering the core data indicates that the estimation of these parameters in the studied field has been evaluated to be highly accurate26. Porosity and permeability tests indicated that a 5% increase in clay content (i.e., kaolinite, smectite, and illite) leads to a reduction in both porosity and permeability of rocks. Moreover, the findings suggested that the alterations in clay content have a more pronounced effect on permeability compared to porosity. A rise in the clay content significantly enhances the formation of chemical products in the reactions between cement and clay minerals. Challenges arise in accurately determining clay volume parameter which is crucial for reservoir characterization, impacting estimations of water saturation and reservoir productivity. An elevation in clay volume estimation may diminish rock porosity calculations, crucial for net pay value assessment27,28,29,30.

Various analytical methods, including X-ray fluorescence (XRF), Fourier transform infrared spectroscopy (FTIR), and X-ray diffraction are crucial for comprehensive clay mineral characterization1. Among these, the XRD stands out as a popular and non-invasive method31,32, which has been employed in geological studies and also for the purpose of identifying and measuring minerals in earth science33,34,35,36,37,38,39,40. It aligns well logs with mineralogical data effectively, especially for fine-grained clay minerals41. However, traditional XRD methods require careful sample preparation, risking damages or alterations to clays7,32,42,43. Scanning electron microscope (SEM) images can aid in precise microstructure examination44.

Previously, Ardebili et al.26 conducted a study on estimating the shale volume of the Kashfroud Formation in the Khangiran Gas Field using an artificial neural network. In the present study, however, the authors address a crucial gap in the comprehensive understanding of clay mineralogy in the formation which is essential for optimizing hydrocarbon extraction from clastic reservoirs. While prior research has focused on various clay mineral properties and their effects on reservoir characteristics globally, detailed studies specific to the formation, particularly with a comparative approach across distinct laboratory techniques, e.g., powder X-ray diffraction (PXRD), polarizing light microscopy, and scanning electron microscopy (SEM/EDX) remain limited. The present study’s innovation lies in its integrative approach, leveraging these three analytical methods to offer a precise understanding of the distribution, type, and abundance of clay minerals within the Formation. By combining these methodologies, the investigation aims to provide more accurate insights into clay mineral impacts on reservoir properties, a knowledge that is pivotal for enhancing reservoir evaluation and extraction techniques in the region.

Geological setting

This study has been conducted on the Kopet-Dagh sedimentary basin extending from northeastern Iran into southern Turkmenistan and east of the South Caspian Basin45,46. Positioned within the collision zone of the Iran Microplate and the Laurasian plates, the Kopet-Dagh Range plays a crucial role47. The Kopet-Dagh sedimentary basin stands as Iran’s second most significant hydrocarbon region following the Zagros and provides two main reservoir formations: Shurijeh (contains sweet gas), and Mozduran (contains sour gas)45 has been located in this basin and the Kashafrud Formation is the last drilled formation in the wells of Khangiran Gas Field.

The Khangiran Gas Field, has been located in this basin and Kashafrud Formation45,46. According to the National Iranian Oil Company (NIOC) report, drilled formations in the Khangiran Gas Field from surface to depth, including Khangiran, Chehel Kaman, Pestehleigh, Kalat, Neyzar, Abtalkh, Abdaraz, Aitamir, Sanganeh, Sarcheshmeh, Tiregan, Shurijeh, Mozduran, and Kashafrud48. The stratigraphic columns of the understudied Well, which the authors called it well #3, in the Khangiran Gas Field is depicted via version 5 of Strater software (Fig. 1). These formations date back from around 35 million years ago in the late Eocene period to approximately 165 million years ago in the upper Middle Jurassic period45.

Fig. 1
figure 1

The stratigraphic columns of the well #3 in the Khangiran Gas Field, plotted via Strater software.

The Khangiran Gas Field, situated in northeastern Iran within the Kopet-Dagh basin, stands out as a significant field and holds a considerable importance from the hydrocarbon perspective. Its structure spans northwest to southeast within the northeastern region of Mashhad province of Iran, and bordering Turkmenistan. This Field is located approximately 180 km northeast of Mashhad city and 25 km northwest of Sarakhs49,50.

The research studies by Afsharharb51 and Poursoltani and Gibling52 demonstrates that the Kashafrud Formation, dating (aging Middle Jurassic), serves as a hydrocarbon source for gas fields in the Kopet-Dagh basin. Furthermore, sedimentological and geochemical studies conducted by Ershadinia et al.53 also confirm the Kashafrud Formation’s reservoir potential. This formation, primarily composed of shale, sandstone, and conglomerate, extending over a significant area in the Kopet-Dagh basin in the northeastern area of Iran52. In the Formation, Jurassic siliciclastic unit (Lower–Middle Jurassic) is considered the main source of hydrocarbon within the Kopet-Dagh basin54. From a diagenetic perspective, the Kashafrud sandstones exhibit signs of various diagenetic activities, such as compression, solidification, dissolution, and creation of minerals formed within the rock47.

Materials and methods

Core data

The present investigation was carried out on core samples extracted from well #3, originating from the Kashafrud Formation situated within the Khangiran Gas Field. The drilled well aimed to investigate the hydrocarbon potential of the Kashafrud Formation, and assess the production of hydrocarbon from the Formation. The well #3 was specifically drilled down to the Kashafrud Formation, reaching a drilled thickness of 433 m. According to the published report made by NIOC48, during the drilling phase within this gas formation, a total of nine core boxes with 0.9 m length were recovered, which belonged to the depths between 3080.5 and 4397.5 m. Subsequently, 10 representative core samples were selected and extracted. Furthermore, a database consists of the depth and code of samples belonging to each core sample was collected.

Core washing is a highly critical process to ensure the removal of any surface contaminants, achieved through a series of rigorous cleaning and crushing steps. Recent advancements in core washing techniques emphasize the importance of reducing particle size and using solvents, e.g. propanol to maintain sample integrity during analysis13,39. In order to avoid from any possible contaminations, the inner parts of each core sample were crushed and then analyzed. The crushing stage provided a size of less than millimeter, and then core fully mixed. After crushing stage, each sample was ground using propanol solution to avoid any possible phase changes. To avoid non-clay minerals contaminating the fine grain fraction, specific clay minerals tests were conducted on a coarser fraction (< mm)31. In contrast, the general XRD analysis was carried out on the milled sample (< 75 µm fraction).

Wireline logging data

From wireline logging stage conducted in the Kashafrud Gas Formation, a dataset of 2263 logging data with a depth distance of 0.061 m were available. Due to inconsistencies in the depths of given core samples, a depth matching process was conducted. This involved a straightforward approach of averaging between upper and lower depths for depths where well logging dataset was not available. These samples underwent a comprehensive analysis to distinguish the type of sand contents in the clay minerals, aiming to quantify their meticulous and precise proportions accurately.

Method

Core analyses

To analyze 10 core samples, the X-ray diffraction (XRD) device was used, located in the X-ray laboratory of the School of Mining Engineering, University of Tehran, to determine how much clay content existed in the Kashafrud Gas Formation. To accomplish this, in the first step of this study, the general X-ray diffraction analysis was performed from 4 to 70° 2θ on 10 core samples. An area of 4.91 cm2 and a depth of 1 mm were analyzed on a circular sample holder made of plastic for full-rock XRD. An advanced Bruker D8 diffractometer, manufactured by Bruker AXS company, a joint product of Germany (Siemens) and the United States (Bruker company), which works with copper radiation Kα = 1.5406 Å, with the θ–2θ goniometry was employed to run the measurements. The samples were exposed to Nickle-filtered Cu Kα radiation at 40 kV and 30 mA with an angular speed of 1.2 \(^\circ\)/min. In the next step, the clay fraction was extracted by the proper procedure described in the subsequent sections, and then the PXRD analysis was performed from 4 to 40° 2θ on the clay fraction to fully identify the type of clay minerals.

From the 10 core samples taken from well #3 of the Khangiran Gas Field, three selected samples included codes #1, #3, and #9 for the preparation of thin section and polarizing microscope studies were chosen. The thin sections were prepared from uncrushed core samples at the School of Mining Engineering. Following this, studies on these samples were carried out by the Zeiss Axioplan 2 polarizing light microscope device.

In order to detect the type of constituent minerals, texture, crystals’ size, and microscopic conflicts, four samples coded #5, #7 (two samples) and #10, respectively underwent preparation using the Scanning Electron Microscope (SEM/EDAX) Camscan MV2300 model equipped with elemental analysis model XFlash 6l10, made by Bruker, Germany. These samples were studied in the laboratory of the School of Metallurgical Engineering at University of Tehran. To enhance conductivity, thin coatings were applied to each sample with several rounds of gold–palladium coating. For each sample, imaging with the secondary electron (SE) method involving low energy and penetration depth was used to determine surface characteristics and show variations in the samples. Besides, the backscattered electron method in the compositional mode was used to perform point analysis. The voltage used by the device during imaging was 20 kV, with a depth of field ranging from 15.9 to 16.7 mm.

To identify clay minerals, laboratory studies, which will be discussed in detail in the subsequent section, included the preparation of core samples, imaging, general XRD analysis with interpretation, SEM/EDAX analysis on the selected samples, and preparation as well as thin sections analyze. Then, for the specific identification of clay minerals types by the PXRD test, special preparations were made (incl. the removal or minimization of non-clay minerals using both physical and chemical methods, saturated processing with ethylene glycol steam, and heating to a temperature of 550 °C).

Laboratory studies

The preparation of samples

The preparation of core samples from well #3 in the Khangiran Gas Field (Fig. 2A) in the Mineral Processing Laboratory of the School of Mining Engineering was accomplished by crushing the core samples into several centimeter-sized pieces with a geological hammer (Fig. 2B). A geological hammer was used with the aim of reducing the dimensions of the given core samples from a minimum of 10 cm to a maximum of 2 cm. Then, for the next tests, a powder with dimensions of 75 µm was ground into a powder device (polarizer) for 30 min (Fig. 2C). The powdered sample is shown in the watch glass in Fig. 2D. Also, a piece of lump was kept for conducting thin section studies from the selected core samples.

Fig. 2
figure 2

(A) the core sample before crushing, (B) the core sample after reducing dimensions with a geological hammer, (C) the powdered sample inside the watch glass, and (D) the powdered sample inside nylon.

PXRD analysis

The X-ray method is an approach for general analysis and identification of minerals and specialized identification of clay minerals which provides useful information about their composition, type, and content31. Using the advanced diffractometer 2θ–θ, D8-Advance model, in the X-ray laboratory of the School of Mining Engineering, University of Tehran, the minerals that made up 10 core samples were identified. It is worth to mention that the purpose of this analysis is to identify the type of clay minerals and their weight percentage. In this test, the weight of the samples was at least 3 g in the form of powder with a size of 75 µm (under sieve 200 mesh). Constituent minerals in the studied samples, based on the PXRD analysis report, in order of abundance, were: clay minerals, alkali feldspars, plagioclase, ankerite, and pyrite. In Fig. 3, the average weight percentage of minerals of all samples is drawn separately.

Fig. 3
figure 3

The average weight percentage of minerals of the samples, obtained from the PXRD analysis.

According to the pie chart, quartz mineral had the highest percentage of constituent weight among the eight minerals identified in the samples. Also, illite and alkali feldspar minerals stood at the second and third ranks, respectively. According to Table 1, samples #10 and #4 had the highest and lowest amounts of clay with a total of 33.5% and 6.20%, respectively. Among these minerals, samples with high clay content were chosen for the SEM studies, while those with low to medium clay content were assigned for the thin section studies.

Table 1 The weight percentages and total weight percentages of the most common clay minerals in the samples.

Specialized identification of clay minerals by the PXRD test

The preparation of samples for quantitative PXRD test involves eliminating potential interferences, especially when clay minerals are not the primary components. Different treatments were used to reduce non-clay mineral interference, enabling precise identification of clays for analysis. These methods involved separating minerals, removing impurities such as iron oxides, and employing settling techniques. The decantation method played a crucial role in isolating clay fractions for a successful PXRD identification.

The specialized identifications revealed the clay minerals in the sample in order of abundance were illite, chlorite, and kaolinite. Figures 4 and 5 present both the general XRD and comparative graphs for these tests, respectively, for two of the samples. The percentages of the most common clay minerals in all samples along with the average and total weight percentages values are presented in Table 1. The qualitative identification curves of the clay minerals from these two samples are drawn in Figs. 6 and 7. These samples had the highest and lowest amounts of clay compared to the other samples.

Fig. 4
figure 4

The general XRD spectrum of the sample with the highest clay content. The identified phases are as follow: Chlorite (Chl), Illite (Ill), Quartz (Qtz), Feldspar (Fsp), Kaolinite (Kln), Ankerite (Ank), Pyrite (Py), and Plagioclase (PI).

Fig. 5
figure 5

The general XRD spectrum of the samples with the highest clay content. The identified phases are as follow: Chlorite (Chl), Illite (Ill), Quartz (Qtz), Feldspar (Fsp), Kaolinite (Kln), Ankerite (Ank), Pyrite (Py), Plagioclase (PI), and Alkali Feldspar (K-Fsp).

Fig. 6
figure 6

Qualitative identification curves of clay minerals from one of the core samples in the studied formation with the highest clay content. (A) Diagram dried in air after removing or minimizing non-clay minerals, (B) diagram heated to 550 °C, (C) graph of sample saturated with ethylene glycol, and (D) graph after boiling in hydrochloric acid.

Fig. 7
figure 7

Qualitative identification curves of clay minerals from one of the core samples in the studied formation with the lowest clay content. (A) Diagram dried in air after removing or minimizing non-clay minerals, (B) diagram heated up to 550 °C, (C) graph of sample saturated with ethylene glycol, and (D) graph after boiling in hydrochloric acid.

The results of the PXRD test are generally based on weight percentage. In order to have a comprehensive understanding from the results of the laboratory examinations, the weight percentage of each samples converted into volume percentage (Table 2).

Table 2 The density of samples, volume percentage, and total volume percentage of clay minerals.

Based on the results of the experiments, it can be concluded that sample #10 with a total of 12.8% at a depth of 4390 m had the highest content of clay minerals. Conversely, sample #2 with a total of 6.7% at a depth of 4397 m accounted for the lowest content of clay minerals. Figure 8 shows the volume percentages of these common minerals in all samples.

Fig. 8
figure 8

The volume percentages of clay minerals in the selected samples.

Polarizing light microscopy studies

According to 10 core samples taken from well #3 of the Khangiran Gas Field, samples with codes #1, #3, and #9 were selected for the preparation of thin section and light microscope studies. This selection was based on their clay content which means that these three samples had low to medium clay content compared to other remained samples. Thin sections were prepared from crushed core samples. Subsequently, Zeiss Axioplan 2 polarizing light microscope was used for detailed examination in the Optical Mineralogy and Mineralography Laboratory at the School of Mining Engineering. The sandstones in the samples were classified following the Folk (1980) classification system. The samples then were analyzed based on petrographic observations considering lithological, textural, structural criteria, sedimentary facies, diagenetic changes, and sedimentary environment of sandstones.

The sample, labelled #1 (Fig. 9A and B) corresponds to the classification of arkose according to the Folk’s classification. The studied sample falls within the clastic sedimentary rock group, specifically in the subgroup of sandstones. The main grains comprising the stone’s framework consist of quartz and feldspar, while minerals with low abundance (3–5% by volume) include chlorite, biotite, and muscovite. The size of grains varies between 80 and 300 µm, with an average exceeding 125 µm. Therefore, the sample is categorized as fine to medium-grained sandstones with moderate sorting. The grains were semi-rounded to semi-angular shape. Some quartz fragments show stretch and extension, indicating a locally semi-oriented orientation and angular—stylolite bordersas a result of diagenesis processes and bearing compressive stress. Also, the effects of recrystallization can be seen in some quartz crystals in the form of polycrystalline fragments. Feldspar grains were mostly displayed in the form of argillite and sericite, their complete pieces were rarely found. Autogenous minerals resulting from feldspars alteration were mainly in situ and had a little movement. Small pieces of rock fragments were of the type of chert with interlayers of mica. The matrix was less than 5% by volume and it was considered as semi-mature to mature. The matrix includes clay minerals, phyllosilicate minerals e.g., chlorite, sericite, muscovite, and microcrystalline silica. The distribution of matrix in the rock was heterogeneous and was locally more concentrated in some parts. Furthermore, the cementation was small (less than 1%) and included anhydrite, iron oxide-hydroxides, clay minerals, sericite, muscovite and chlorite.

Fig. 9
figure 9

The general views of sample #1: (A) Quartz in the form of polycrystalline pieces, and (B) Quartz including clay plus phyllosilicate minerals.

Scanning electron microscope analysis (SEM/EDAX)

Among 10 available samples, according to the percentage of clay minerals, three samples were selected for SEM/EDAX analysis. These three samples coded #5, #7 (1 and 2), and #10 were studied after preparation with Camscan MV2300 scanning electron microscope equipped with XFlash 6l10 elemental analysis at the laboratory of School of Metallurgical Engineering, Faculty of Engineering, University of Tehran. Their images were recorded and are present in Fig. 10.

Fig. 10
figure 10

The SEM images of selected samples, including: (A) Sample #5, (B) Sample #7—the first image, (C) Sample #7—the second image, and (D) Sample #10.

The interpretation of SEM images of clay lam sample extracted from sample #5

The magnification of the studied image (Fig. 10A) is 2500 times and the depth of field is 20.4 mm. In the studied field, clay minerals of kaolinite, illite, and chlorite had stuck together in a coated form due to their adhesive properties and had created concretion structures (circle-like, spherical structures). It can be assumed that all kinds of clay minerals were placed inside these structures with a coated form. In addition, plate-shaped clay minerals were also seen in the background. Kaolinite clay minerals were visible in locations J-7, I-6, K-7, N-8, and O-8. Additionally, kaolinite clay mineral in the background was spread in the plates form in the area from M-3 to M-7. The maximum size of kaolinite clay pieces in the studied sample is 40 µm, with an average size of 4.7 µm and a minimum size of 2 µm. Regarding the sizes of illite and chlorite fragments, they cannot be measured because of being coating and inter-layered shape, forming spherical shapes. The distribution pattern of the clays in this sample was pore-filling type, and in a few cases, it can also be classified as covering. The pore-filling pattern can be clearly seen in different parts of the studied field.

The interpretation of SEM images of clay lam sample extracted from the first image of sample #7-1

The magnification of the studied image (Fig. 10B) is 5000 and the depth of field is 20 mm. In this sample, clays were better separated from each other and clay plates were completely visible. However, similar to the previous sample, there was a very different dimensional distribution. That is, from very small pieces of clay, e.g., kaolinite clay in C-7 and D-5 locations, to larger types in D-6 and C-5 locations. Clays were strongly intertwined and stuck together. The advantage of the image being studied was that the plates can be seen independently of each other with a high resolution. The distribution pattern of the particles was also similar to the previous example, i.e., as a pore-filling and in some cases as a covering pattern. The minimum, average, and maximum particle sizes are 1.4 μm, 2.7 μm, and 13.7 μm, respectively. Moreover, a highly heterogeneous distribution of clay minerals types e.g., illite, chlorite, and kaolinite with various sizes was available in the studied field. In locations G-4 and B-2, one or more illites can be seen scattered in a small number. In the H-2 location, illite can be observed in the form of fibrous. In the range from F-7 to H-4, there were areas where a mixture of chlorite and kaolinite which can be seen in an intertwined form.

The interpretation of SEM images of clay lam sample extracted from the second image of sample #7-2

The magnification of the studied image (Fig. 10C) is 5000 and the depth of field is 20.1 mm. In another field of the same sample (code #7), the coated clay minerals were significantly crushed, forming a concretion. It should be noted that concretions were a mixture of clay minerals e.g., chlorite, kaolinite, and illite.

Some parts of the image that appeared very bright in color, (e.g., position N-11), seemed to contain iron-bearing clay fragments, possibly of the ferric chlorite type. This brightness could be attributed to the presence of elements with high atomic numbers (e.g., iron), found in chlorite. In this image, there was a mixture of clays e.g., kaolinite, chlorite, and illite coated and inter-layered. The only kaolinite piece that existed in an open form is the O-15 piece, and measuring its dimensions was not possible due to the pieces being chipped.

The interpretation of SEM images of clay lam sample extracted from sample #10

The magnification of the studied image (Fig. 10D) is 5000 and the depth of field is 1.20 mm. In sample #10, which was related to a shallower depth than the other samples, there was a lower degree of scaling and kaolinite pieces were mainly in the form of rounded plates. Some of these places are: B-10, C-10, D-11, E-11, E-12. The clay minerals in this sample appeared separately, although they were actually in a location where the background was a mixture of different types of clay minerals, e.g., kaolinite in the C-13 location, which was on a matrix. Kaolinite was the dominant mineral observed in the study field, but illite and chlorite can also be identified in the matrix in an interlock form. The maximum, average, and minimum sizes of kaolinite particle are 12.3 μm, 5.5 μm, and 2 μm, respectively. An important issue was the rotation of the kaolinite plates in this sample, potentially attributed to the shallow depth of this sample or other diagenetic reasons, e.g., alteration because of changes in temperature and pressure. The dominant pattern of clay distribution was of the covering type and the sub-pattern in this field was the pore-filling pattern. Moreover, a very long linear illite fragment with kaolinite fragments was observed from locations F-11 to D-12. In general, the SEM studies confirm the identification of clays types from the PXRD test. Furthermore, the distribution pattern of clay minerals and their changes in particle size were also investigated.

The field seen in the SEM diagram of Fig. 11 was subjected to general elemental analysis and the results are as follows. According to the diagram, the main peaks for three elements namely silicon, aluminum, and oxygen are clearly noticeable (the gold peak appears due to the coating of the sample with a gold–palladium alloy), which are the constituents of aluminosilicates or clay minerals in the sample. The origin of silicon and oxygen elements can also be related to the presence of quartz in the sample. Because the SEM analysis was conducted by performing the test on the slide of the clay section extracted from the sample, the amount of quartz was minimized and not completely removed due to the resistant nature of this mineral. Therefore, its peaks can be seen in the form of Si and O elements.

Fig. 11
figure 11

The elemental analysis (EDAX) of sample #10.

The fourth abundant element is sodium, which is related to plagioclase or some alkaline feldspars. Therefore, given the simultaneous presence of plagioclase and alkali feldspar in the PXRD test, sodium might be associated with these two minerals either separately or together. Calcium ranks as the fifth most abundant element in Fig. 11. Considering the presence of carbonate phases, e.g., ankerite and calcite, as indicated by the PXRD results, hence calcium is associated with these carbonate phases. Thus, calcium can be attributed to ankerite minerals or various types of classical plagioclase.

The sixth most abundant element is sulfur. The presence of this element can be justified due to the appearance of the cubic pyrite mineral in the PXRD test. The next long peak is related to the Fe element, which is probably due to the presence of ancrete, pyrite, and chlorite-type clay minerals. The eighth abundant element is Mg, undoubtedly present in chlorite. In total, the main phases of the sample are: oxygen, silicon, to some extent sodium, aluminum and calcium. The general elemental analysis prepared from the sample in the form of a map (Fig. 12) indicates that the dominant constituent phases of the sample were made of oxygen, silicon and aluminum. This means that the main producers were clay and quartz minerals. The results of the elemental analysis of sample #10 are listed in Table 3.

Fig. 12
figure 12

The elemental analysis diagram in the form of a map belonged to sample #10.

Table 3 The elemental analysis results (EDAX) of sample 16.

It can be seen that the majority of samples in the studied field were composed of silicon, oxygen, and aluminum elements, which are the main producers of clay minerals (aluminosilicates). Figure 13 shows the general view of the sample, presenting the elements separately.

Fig. 13
figure 13

(A) The studied field analyzed for these elements: (B) Oxygen (O), (C) Potassium (k), (D) Sodium (Na), (E) Magnesium (Mg), (F) Calcium (Ca), (G) Iron (Fe), (H) Sulfur (S), (I) Aluminum (Al), (J) Silicon (Si) and (K) Phosphorus (P).

Results and discussion

The evaluation of the present investigation indicates that the burial depth exceeded 5 km, the burial basin temperature of the study area was above 150 °C and it is classified as the metagenesis stage, which only gas generation occurs during this stage55. Sandstones, as the most significant reservoir rocks for both oil and gas, represent porous siliciclastic properties and manifest a wide range of textures and compositions due to their deposition in diverse environments. These minerals, especially the fine-grained ones, often contribute to a substantial reduction in permeability of a reservoir by migrating and consequently obstructing interconnecting pore throats. Found in both sandstones and certain carbonate formations, these clay minerals are mobile, particularly when encountering fresh or low-salinity water due to the absence of confining pressure. The type, ratio, and distribution mode (be it as cement or pore filling) of clay minerals within sandstones exert a significant influence on petrophysical properties and overall reservoir qualities56,57,58,59,60,61,62.

When clay minerals are trapped in tiny pore throats, it can lead to a dramatic increase in fluid pressure needed for sustained flow, resulting in an overall reduction in permeability across porous media. Most hydrocarbon formations contain these minerals, and their interaction with water, especially low-salinity fluids, prompts fines migration and clay swelling. The porous structure of sandstone represents a network of pores connected through narrower channels. The dispersed presence of clay minerals within these pore spaces and permeability channels leads to provide a decrease in pore throat dimensions, subsequently, overall permeability within the reservoir rock. Among the prevalent clay minerals in sandstone reservoirs, authigenic kaolinite often forms flakes that are considered pore-filling type57,59,61,62. The accurate estimations of both the water saturation and porosity can be obtained by finding out how much clay is present within a reservoir. Bedsides, the calculation of the effective porosity can be done more accurately by having a good knowledge of the reservoir clay content63,64.

Based on the results of PXRD studies which were also confirmed by other laboratory methods, the minerals of Kashafrud Formation include aluminosilicates, quartz, plagioclase, alkali feldspar, ancrete, calcite, pyrite, and chlorite. The determination of clay content in the reservoir formations is of top priority, because of the importance of the subject. To identify the constituent clays in the studied samples accurately, specialized tests were performed to identify the clay minerals, and the results of these tests showed that clay minerals in the formation in terms of abundance were illite, chlorite, and kaolinite.

Based on the results of polarizing light microscopy, the sample was considered a sandstone classified by Folk’s system, featuring primarily quartz and feldspar grains. It had also fine to medium-grained structure, shows evidence of diagenesis, with semi-rounded to semi-angular grains, and a matrix comprising clay minerals and phyllosilicates, unevenly distributed within the rock. The SEM studies were performed to determine the distribution pattern of the clay minerals. Based on these studies, the distribution pattern of clay minerals was mainly of pore filling type and in some limited cases, pore coating type was identified. In addition, the size of clay minerals was a minimum of from 1.4 μm to a maximum of 13.7 μm. The performed EDAX analysis indicated that silicon, oxygen, and aluminum elements were witnessed in the majority of samples taken from the Kashafrud Formation which the minerals that were identified through the PXRD method were further validated through subsequent PXRD tests.

Conclusions

Prior to this study, comprehensive research on the identification and distribution of clay minerals within the Kashfroud Formation was lacking. This study addressed this gap and provided original insights valuable to petroleum exploration and mining engineering. Utilizing polarizing light microscopy, scanning electron microscopy equipped with energy dispersive X-ray spectroscopy (SEM/EDX), and X-ray diffraction (XRD), the clay mineralogy of the Kashafrud Gas Formation was analyzed. 10 core samples from well #3, positioned at depths of 4320–4753 m within the Khangiran Gas Field, northeast Kopet-Dagh sedimentary basin, Iran, were examined.

Qualitative analysis revealed that the Kashafrud Formation primarily consists of illite, chlorite, kaolinite, quartz, and alkali feldspar. The distribution of clay minerals, determined through specialized tests, followed the order: illite (1.11 wt%), kaolinite (3.4 wt%), and chlorite (4.7 wt%). Polarizing light microscopy indicated that the studied sample was a fine- to medium-grained sandstone composed mainly of quartz and feldspar grains, exhibiting diagenetic features, such as stretching and extension in quartz fragments. The SEM/EDX analysis further detailed the presence of kaolinite, illite, and chlorite forming concretion structures, with kaolinite averaging 4.7 µm in size. Moreover, variations in clay distribution patterns, including pore-filling and covering, were observed, highlighting heterogeneity in particle sizes and orientations.

The proposed approach not only provides a quantitative assessment of clay mineral distribution owing to advanced analytical techniques but also accounts for the inherent heterogeneity and diagenetic influences affecting mineral composition. The findings emphasize the role of diagenetic processes and environmental factors in mineral alterations within the formation.