Lithologic Logs Research Articles

Attention mechanism‐assisted recurrent neural network for well log lithology classification

AbstractLithology classification is a fundamental aspect of reservoir classification. Due to the limited availability of core samples, computational modelling methods for lithology classification based on indirect measurements are required. The main challenge for standard clustering methods is the complex vertical dependency of sedimentological sequences as well as the spatial coupling of well logs. Machine learning methods, such as recurrent neural networks, long short‐term memory and bidirectional long short‐term memory, can account for the spatial correlation of the measured data and the predicted model. Based on these developments, we propose a novel approach using two distinct models: a self‐attention‐assisted bidirectional long short‐term memory model and a multi‐head attention‐based bidirectional long short‐term memory model. These models consider spatial continuity and adaptively adjust the weight in each step to improve the classification using the attention mechanism. The proposed method is tested on a set of real well logs with limited training data obtained from core samples. The prediction results from the proposed models and the benchmark one are compared in terms of the accuracy of lithology classification. Additionally, the weight matrices from both attention mechanisms are visualized to elucidate the correlations between depth steps and to help analyse how these mechanisms contribute to improved prediction accuracy. The study shows that the proposed multi‐head attention‐based bidirectional long short‐term memory model improves classification, especially for thin layers.

Pattern recognition approach (PRA) for identifying oil reservoir lithology of Camaal oil field, Yemen

The accurate determination of reservoir lithology remains a challenge in petroleum engineering. There are some conventional techniques available to determine the lithology. However, the application of those techniques has been long and complex. So, the main goal of this study is to simplify the identification of reservoir lithology. This paper presents a Pattern Recognition Approach (PRA) to identify the reservoir lithology simply and accurately. It is type of artificial neural network. Four wells from the Camaal Field were chosen to develop this approach. Around 32400 data points from the previous wells were digitized. The PRA approach used depth, gamma ray, lithology, sonic, neutron, and density logs as inputs. The model classifies lithology into permeable and impermeable rocks, further categorizing them into clastic and carbonate rocks, and subsequently into specific types into sand, sandstone, dolomite and limestone. The results show that the proposed approach provides a suitable prediction of lithology with higher accuracy compared with actual lithology. The model demonstrates high accuracy rates in identifying various lithologies, with overall accuracies of 76.2% for permeable/impermeable rocks, 94.9 for clastic/carbonate rocks, 86.2% for sand/sandstone, and 92.8% for dolomite/limestone.

Evaluation of specific retention, specific yield, and storage-dependent drainability efficiency in a coastal milieu via geo-electrical technology

ABSTRACT The aim of this research was the application of geo-electrical technology in the determination of specific retention, specific yield, and storage-dependent drainability efficiency within a major coastal milieu in Nigeria. This improves on past work where expensive and time-intensive pumping tests were employed to determine specific retention and specific yield. In addition to the determination of these key aquifer geo-kinetic properties, other important aquifer characteristics such as aquifer potentiality, protectivity, and vulnerability to contamination within this major coastal system were determined. Geo-electrical technology employed Wenner and Schlumberger arrays to undertake Vertical Electrical Soundings (VES) and Electrical Resistivity Tomography (ERT) surveys and results obtained were constrained by ground truth from lithological logs. Results from the electrical resistivity surveys indicated that the lithological strata comprised motley topsoil, coarse sand, fine sand, and sandy clay. Measures of specific yield and specific retention were generated as secondary geo-electrical indices. Specific yield intersected with specific retention at 73% storage-dependent drainability efficiency (SDE), indicating that this percentage must be exceeded for optimal groundwater extraction from the rock matrix's pore spaces during pumping. Aquifer transmissivity measures indicated high aquifer potentiality; longitudinal conductance measures indicated poor aquifer protectivity implying increased vulnerability of the aquifer to contamination.

Estimation of Groundwater Potential and Aquifer Protective Capacity Within Old Ikenga Hotel, Nsukka, Enugu State, Nigeria

A geophysical survey using electrical resistivity method was conducted around Old Ikenga Hotel, Nsukka in Enugu State, Nigeria to investigate the groundwater potential and aquifer protective capacity of the area. The project area lies within latitudes 6°50'4.0''N – 6°57'52.0"N and longitudes 7°21'6.3"E – 7°28'12.0"E, and covers an area of about 89.6 km2. Vertical Electrical Soundings (VES) were carried out with a digital read out resistivity meter (ABEM SAS 1000). The VES points were marked at 25 m and 75 m along a 100 m line. A total of eight soundings were carried out in the area. The VES data collected were interpreted using INTERPEX software and the results presented in terms of resistivity, thickness, depth and lithology. The lithology was inferred by correlating the result to the lithology log of the borehole drilled in the hotel and the geology of the study area. The VES result shows lithologic layers varying from 4 to 5. Aquiferous sand and Aquiferous sandy shale constitute the aquifer units in the area at depth of 30.26 m to 188.20 m, with their thickness ranging from 30 m to 74 m as shown by their isopach map. The aquifer protective capacity was determined by calculating for longitudinal conductance and matching the values to known standards. The calculated longitudinal conductance varies from 0.0409 to 3.1235 mhos. The interpreted VES results reveal poor, moderate, good and very good aquifer protective capacities of the overburden layers.

PALYNOLOGICAL, PALEOENVIRONMENTAL AND PALEOCLIMATE ANALYSES OF EMI-5-WELL, OFFSHORE NIGER DELTA, NIGERIA

Palynological, Paleoenvironmental and Paleoclimatic investigations of Emi-5 well Offshore, Niger Delta, Nigeria was conducted to established palynozones, determine age, reconstruct paleoenvironmental conditions and paleoclimatic conditions to permit stratigraphic correlation, provide insights into the depositional settings and changes overtime within the basin. Fifty one (51) shaly ditch cutting samples were subjected to standard palynological method involving sample maceration with 10 % HCl and 40 % HF acids to remove carbonates and siliceous components, neutralization with distilled water and sieving through 5-micron sieve, followed by separation of palynomorphs from residue prior to mounting of the grains on glass slides for further analysis. The palynological composition is marked by dorminance of pollen and spores over dinoflagellate cysts. The intercalation of the sandy and shaly intervals across the section evidenced from lithologic, textural, and wire line log data suggested that the entire studied interval belong to the Agbada Formation. Three interval range zones were delineated using the international stratigraphic guide for the establishment of biozones, they are Stereisporites sp – Verrutricolporites rotundiporus; Verrutricolporites rotundiporus – Gemmemonoporites sp. and Retistephanocolpites gracilis – Pachydermites diederixi Zones and based on the aforementioned zones the interval is assigned middle Miocene to late Miocene. Coastal-deltaic environment of deposition is inferred for the studied interval on the basis of occurrence and dominance of diagnostic forms such as Zonocostites ramonae,Rhizophora, Psilatricoloporites crasssus and Acrostichum aureum while good representation of these forms and the identified palyno-ecological groups placed the interval on a wet climatic zone.

Correlation of the Coniacian and Santonian stages of the Upper Cretaceous in the Anglo-Paris Basin

The correlation of the Coniacian and Santonian chalks of the Anglo-Paris Basin is described on the basis of detailed lithological logs and extensive records of macrofossils and microcrinoids. In the almost complete absence of ammonites, inoceramid bivalves afford the highest resolution correlation of these stages in chalks, but their value here is limited by the absence of key genera and species, most notably in the Upper Coniacian and middle and Upper Santonian. Echinoids and other macrofossils (brachiopods, stalked crinoids, belemnites) have proved useful, but many are long-ranging or uncommon. Some marker beds, including flints and marl seams, provide useful correlations across the basin, but are locally absent. For the Upper Santonian, the stemless benthonic crinoids Uintacrinus and Marsupites provide high-resolution correlation, both within the basin and to other regions. The successions on the basin margins, in the far north of France (Nord, Pas de Calais) and the southwest (Touraine) are condensed and yield ammonites in association with important inoceramid species. The controls on sedimentation caused by sea-level changes are evaluated on a basinal and global scale, most especially for the Upper Santonian.

Assessment to locate potential deep aquifer systems using lithological logs, pumping tests and electromagnetic surveys in hard rock terrain of Dodoma urban area, Tanzania

The Dodoma region located in Central Tanzania is a semi-arid area characterised by scarce surface water resources. Because of climate change and dependence on groundwater resources in the Dodoma urban area, there is currently an increased demand for water due to low groundwater recharge and high groundwater withdrawals due to the growing population. There is therefore a need to explore potential deep aquifer systems around the Dodoma urban area to meet the increasing water demand. This study uses lithological logs, pumping tests, and electromagnetic survey data to delineate the subsurface lithologies to locate potential deep aquifer system in the area. Results show that the shallow aquifer system is unconfined in some areas with resistivity values ranging from 11 to 28 Ωm and semi-confined in other areas with resistivity values ranging from 19 to 27 Ωm. The unconfined aquifer extends up to 60 m, while the semi-confined aquifer extends between 50–55 m and 65–120 m. The study found a potential deep aquifer system at a depth of between 200 and 290 m, with resistivity values ranging from 11 to 20 Ωm. Currently, this deep aquifer system has not been exploited as most boreholes in the area are only up to 150 m deep. The main lithology in the study area is weathered and fractured granite, with different degrees of weathering and fracturing, indicating a hard-rock aquifer system. This study adds valuable knowledge on the location of potential deep aquifer systems in the area for proper groundwater utilisation and management.

Development and influencing factors of middle‐deep tight sandstone reservoir ‘sweet spots’ in fault basin—Taking the Eocene Shahejie Formation in Suning area of Raoyang Sag in Bohai Bay basin as an example

This study takes the tight sandstone of the lower sub‐member of the first member of the Eocene Shahejie Formation (Es1x) in the Suning area, Raoyang Sag, Bohai Bay Basin. Furthermore, it determines the rock type of the middle‐deep tight sandstone through lithologic profiling, core observation and log data analysis. The genetic control factors of tight sandstone in the study area are analysed using conventional thin section, pore‐casted thin section and scanning electron microscope, and the characteristics and genetic control factors of tight sandstone reservoirs are summarized. The research results show that the tight sandstones in this area are primarily feldspathic sandstone and feldspathic lithic sandstone, with high quartz contents, around 61.50%. The average porosity of the reservoir in the study area is 4.69%, and the average permeability is 2.92 mD, which are characteristics of low‐porosity and low‐permeability reservoirs. Through research, it has been found that the distribution of sedimentary facies and diagenesis in the study area control the development and distribution of ‘sweet spots’. Both compaction and cementation reduce the physical properties of the reservoir, while dissolution can increase porosity, thereby increasing the physical properties of the reservoir; dissolution resulted in an increase of approximately 3.57%–3.69% in the porosity of the tight reservoir in the study area, with an average porosity increase of 3.63%, and is the main controlling factor for the development of ‘sweet spots’ in the reservoir. The results provide an important basis for guiding the exploration and development of tight sandstone oil and gas in the lower sub‐member of the first member of the Eocene Shahejie Formation (Es1x)in the Suning area.

Timescales of mafic magmatic fractionation documented by paleosecular variation in basaltic drill core, Snake River Plain volcanic province, Idaho, USA

The timescales over which fractional crystallization and recharge work in mafic volcano-plutonic provinces is subject to great uncertainty. Currently modeled processes are subject to the scale of measurement: monogenetic basaltic fields accumulate over hundreds of thousands of years, consistent with U-Th-Ra isotopic variations that imply 50% crystallization of basic magmas on timescales of 100,000 years or more, whereas crystal diffusion modeling implies phenocryst residence times of ∼1−1000 years. Monogenetic basalts of the Snake River Plain in southern Idaho, USA, are up to 2 km thick and postdate passage over the Yellowstone−Snake River Plain hotspot. Detailed lithologic and geophysical logging of core from deep drill holes, along with chemical stratigraphy and high-resolution paleomagnetic inclination measurements, document individual eruptive units, compound lava flows, and basaltic flow groups that accumulated over 1−6 m.y. Hiatuses are commonly marked by loess or fluvial interbeds that vary from ∼0.1 m thick to 20 m thick. Radiometric (40Ar-39Ar, detrital zircon U-Pb) and paleomagnetic timescale ages show that the deepest hole (Kimama drill hole, 1912 m total depth) accumulated over ∼6 m.y. Cycles of fractional crystallization and recharge are recognized in the chemical stratigraphy as up-section shifts in major and trace elements; these fractionation cycles commonly represent 40%−50% fractionation. Individual fractionation cycles may comprise 20−40 eruptive units (8−17 lava flows) with little to no change in paleomagnetic inclination (0°−1°), whereas adjacent cycles may differ by several degrees from one another or reflect changes in polarity. Rates of paleosecular variation in Holocene lavas and sediments dated using 14C document significant shifts in magnetic inclination over short timescales, ranging from ∼0.05° to 2°/decade, with an average of ∼0.5°/decade and a minimum rate of 0.05°/decade. This implies that fractionation cycles with ≤1° variation in magnetic inclination formed on timescales of a few decades up to a few centuries (20−200 years). Thus, the lavas collectively represent only a few thousand years of eruptive activity, with major flow groups separated in time by tens to hundreds of thousands of years. We suggest that the rates defined by paleosecular variation capture the timescales of magmatic chamber evolution (fractionation/recharge) in the seismically imaged mid-crustal sill complex; in contrast, we suggest that crystal diffusion modeling captures the residence times in shallow subvolcanic magmatic chambers that underlie individual monogenetic volcanoes.

Geophysical characterization of groundwater aquifers in the Western Debrecen area, Hungary: insights from gravity, magnetotelluric, and electrical resistivity tomography

The recent study followed a multi-methodological approach integrating gravity, magnetotelluric (MT), and electrical resistivity tomography (ERT) to investigate the geometry and hydrological characteristics of the main hydrostratigraphical units in the Western Debrecen area, Eastern Hungary. The integration of these methods aims to delineate potential zones for groundwater development and guide effective extraction strategies. In the gravity investigation, the Bouguer anomaly map undergoes spectral analysis for the separation of shallow and deep features, offering a preliminary indication of basement rock depth. Subsequently, gravity data inversion is employed to map variations in basement rock topography, revealing a basin structure, with sediment thicknesses extending up to 2 km. On the other hand, the MT data are modeled using the 1D Occam inversion algorithm to validate the results of the gravity data analysis. This inversion, constrained with lithological logs is further utilized to delineate the main hydrostratigraphical units in the study area. Accordingly, four units are identified, including the Nagyalföld Aquifer, Algyő and Endrődi Aquitards, the Badenian Aquifer, and the Pre-Neogene Aquitard. Consequently, Dar Zarrouk parameters based transmissivity and the hydraulic conductivity of the aquifer units are measured. The Nagyalföld aquifer showed a hydraulic conductivity that ranged between 7.9 and 11.9 m/day, while the Badenian aquifer showed an average hydraulic conductivity of 13.1 m/day. The ERT data are employed to map the spatial distribution of the depth to the water table. The shallow water table is observed in regions characterized by an elevated thickness of sedimentary rocks, attributed to their high specific capacity. Integrating these hydrogeophysical methods provided a comprehensive understanding of the subsurface hydrology and enabled better-informed decision-making for groundwater development.

Joint interpretation and modeling of potential field data for mapping groundwater potential zones around Debrecen area, Eastern Hungary

The Debrecen area, as part of the Great Hungarian Plain (GHP), is associated with a multi-aquifer system that is overly exploited to fulfill the development plans. This research aims to jointly interpret and model gravity and magnetic data to map the subsurface geology and structures that govern groundwater occurrence. Various potential field techniques, including spectral analysis, anomaly derivatives, analytical signal, and Euler deconvolution were employed to map the distribution and depth of the buried geological structures. The combination of the potential field techniques enabled the construction of a detailed lineament map, providing valuable insights into the distribution of the subsurface structural features. It was indicated that the main structural trend is NW–SE and NE–SW, that coincides with the main structural trends in Hungary. Subsequently, a lineament density map is derived, indicating that the eastern, central, and northwestern parts of the area form the most promising zones for groundwater prospection. The joint inversion of gravity and magnetic data has further enhanced the understanding of subsurface geology. The depth to the basement rock varied between 1.18 and 2.2 km. The highest depth to the basement meets with thick sedimentary sequences bounded by normal faults forming graben and horst structures. Moreover, the distribution of these sediments is investigated using lithological logs indicating the thickness of the main hydrostratigraphic units in the Debrecen area. These units include Nagyalföld Aquifer, Algyő Aquitard, Endrőd Aquitard, and Miocene Badenian Aquifer units, which mainly consist of sand, silt, marl, and gravel. The recent study demonstrated the effectiveness of the joint interpretation in enhancing the knowledge of lithology and geological structures. However, a detailed geophysical survey is recommended to characterize the hydrostratigraphic units in the Debrecen area.

ERT data assimilation to characterize aquifer hydraulic conductivity heterogeneity through a heat‐tracing experiment

AbstractGeothermal energy systems, such as heat pumps relying on aquifers, use renewable sources of energy that are accessible in urban areas. It is necessary to characterize the subsurface hydraulic properties prior to the installation of such systems. In this context, a heat‐tracing experiment is a typical field test that can help with the characterization of the subsurface. During a heat‐tracing experiment, monitoring with downhole temperature sensors, water‐level pressure transducers and electrical resistivity tomography (ERT) can be used to help characterize the hydrogeological properties. Previous monitoring tools have shortcomings, such as low‐resolution data and over‐smoothing; thus, they fail to reproduce the heterogeneity of hydrogeological properties. Ensemble Kalman filter (EnKF) is a promising tool that can overcome the over‐smoothing problem to replicate the hydrogeological property heterogeneity. In this work, we proposed a new procedure to assimilate time‐lapse cross‐borehole ERT data into a numerical model of groundwater flow and heat transfer, where the groundwater is extracted and heated water is reinjected into an unconfined sandy‐gravel aquifer. The finite element model (FEFLOW 7.3) of groundwater flow and heat transfer is integrated with petrophysical relationship and electrical forward modelling (ResIPy) to estimate cross‐borehole ERT measurements. Then, the estimated apparent resistivity is assimilated to update the hydraulic conductivity model using EnKF. The results of the application of the proposed approach to an experimental site located in Quebec City (Canada) demonstrate that the heterogeneity of K is correctly reproduced as the updated K model is reasonably consistent with the lithological log. In addition, the proposed approach was able to replicate the cross‐borehole ERT field and temperature measurements. The comparison between prior and posterior distribution of K with slug test results shows that the EnKF made the final (assimilated) distribution of K move towards K values inferred with slug tests.

Geophysical, hydrogeological and laboratory assessment of groundwater quality and its health implication in some urban communities in cross river state

The electrical resistivity technique which involved the Schlumberger depth sounding method, electrical resistivity tomography and geochemical analyses of water samples collected from boreholes were used to investigate the quality of groundwater in some urban communities in Cross River State, also the health implication of the said groundwater was also investigated. Forty nine electrical resistivity data were collected, modeled, and interpreted after calibration with lithologic logs. Ten electrical resistivity tomography data was were also collected, modeled and interpreted. Forty nine borehole water samples were collected and analyzed to determine anion, cation concentrations and some physical and chemical parameters, such as water colour, temperature, total dissolved solids, and electrical conductivity. The results show that the lithostratigraphy of the study area is composed of sands, sandstones (fractured, consolidated and loosed), siltstones, shales (compacted and fractured) of the Asu River Group, Eze-Aku Formation which comprises the aquifer units, and the Nkporo Shale Formation. The aquifer conduits are known to be rich in silicate minerals, and the groundwater samples in some locations show a significant amount of Ca2+, Mg2+, and Na+. These cations balanced the consumption of H+ during the hydrolytic alteration of silicate minerals. The geochemical analysis of groundwater samples revealed dominant calcium–magnesium–carbonate– bicarbonate water faces. The groundwater quality was observed to be influenced by the interaction of some geologic processes but was classified to be good to excellent, indicating its suitability for domestic and irrigation purposes.

Application of the electrical resistivity method and the estimation of limestone volume: a case study

The present work used the electrical resistivity approach to conduct a three-dimensional modeling and initial volume estimation of the limestone layer in the Mintom region located in southern Cameroon. In order to achieve the objectives of the study, a total of 21 electrical soundings spaced 250 m were first collected in the field using the Schlumberger array. These soundings were conducted along three profiles oriented in an east–west direction, spaced 500 m. Additionally, a geological survey was conducted to identify and emphasize the presence of limestone formations within the designated study region. The interpretation of the sounding data was conducted based on the analysis of the sounding curves. The interpretation outcomes, specifically resistivity and thickness, were compared with the geological field data, resulting in the development of lithostratigraphic logs for each sounding. The geological sections were constructed using the logs of the designated profile. The lithological logs were utilized to establish a lithological interface model and calculate the volume of the limestone layer at 260 ± 13 × 106 m3, utilizing the inverse distance method built into RockWorks software. A resistivity value is assigned to each geological layer in a sounding curve, allowing for the development of a resistivity variation model specific to the limestone layer. The proposed model facilitates the categorization of limestone layers based on their resistivity variations, thus serving as a fundamental reference for prospective exploratory activities within the designated study region. Our integrated approach provides a replicable model for a better understanding of the limestone reserve and effective management of this valuable resource.

An ensemble-based machine learning solution for imbalanced multiclass dataset during lithology log generation

The lithology log, an integral component of the master log, graphically portrays the encountered lithological sequence during drilling operations. In addition to offering real-time cross-sectional insights, lithology logs greatly aid in correlating and evaluating multiple sections efficiently. This paper introduces a novel workflow reliant on an enhanced weighted average ensemble approach for producing high-resolution lithology logs. The research contends with a challenging multiclass imbalanced lithofacies distribution emerging from substantial heterogeneities within subsurface geological structures. Typically, methods to handle imbalanced data, e.g., cost-sensitive learning (CSL), are tailored for issues encountered in binary classification. Error correcting output code (ECOC) originates from decomposition strategies, effectively breaking down multiclass problems into numerous binary subproblems. The database comprises conventional well logs and lithology logs obtained from five proximate wells within a Middle Eastern oilfield. Utilizing well-known machine learning (ML) algorithms, such as support vector machine (SVM), random forest (RF), decision tree (DT), logistic regression (LR), and extreme gradient boosting (XGBoost), as baseline classifiers, this study aims to enhance the accurate prediction of underground lithofacies. Upon recognizing a blind well, the data from the remaining four wells are utilized to train the ML algorithms. After integrating ECOC and CSL techniques with the baseline classifiers, they undergo evaluation. In the initial assessment, both RF and SVM demonstrated superior performance, prompting the development of an enhanced weighted average ensemble based on them. The comprehensive numerical and visual analysis corroborates the outstanding performance of the developed ensemble. The average Kappa statistic of 84.50%, signifying almost-perfect agreement, and mean F-measures of 91.04% emphasize the robustness of the designed ensemble-based workflow during the evaluation of blind well data.

Integrated seismic inversion for clastic reservoir characterization: Case of the upper Silurian reservoir, Tunisian Ghadames Basin

Deltaic petroleum systems are often complex and heterogeneous. Delineation of sandstone bodies in such complex depositional settings constitutes a crucial challenge in reservoir exploration. However, reservoir characterization based on conventional well-logs and stratigraphic interpretation in complex settings often leads to uncertainties. In this work, we present an integrated quantitative interpretation to reveal heterogeneity and the reservoir properties in a complex deltaic environment. The study benefits from integrating borehole lithology logs, electrical well logging, and a 3D post-stack seismic cube. To demonstrate the practical application and effectiveness of our proposed approach, seismic inversion has been applied to the Acacus Formation in the Ghadames Basin of southern Tunisia. We estimate the acoustic impedance volume from 3D seismic data using a stochastic inversion algorithm. This algorithm enables the generation of multiple realizations, accounting for uncertainties and providing a more comprehensive understanding of subsurface impedance. Then, the porosity model is computed by establishing a functional relationship between acoustic impedance and porosity. The analysis of the relationship between acoustic impedance and different rock properties helped to differentiate lithological units of sandstone and claystone. This lithology classification is used to estimate the lithofacies model. Finally, the obtained porosity and lithofacies models are validated with well data. The crucial correlation between facies distribution and porosity allows for accurate mapping of individual thin sandstone bodies and their property distribution beyond well control. The lithofacies model shows significant sandstone bodies in the Acacus A and C units that constitute good reservoirs. The porosity model confirms this, as these sandstone facies exhibit high porosity (20%). The proposed approach has proven to be powerful in reservoir delineation, characterization, and exploration, and can be applied in similar geological setting basins, and frameworks.

Re-evaluation of the Mesozoic–Cenozoic biostratigraphy of the Laurentian Subbasin of the Scotian Basin, offshore eastern Canada

We use new and existing nannofossil, palynological, and microfossil biostratigraphic data in conjunction with lithologic and geophysical logs from four wells to establish a series of sequence-stratigraphic events in the Mesozoic–Cenozoic of the Laurentian Subbasin of offshore Newfoundland, eastern Canada. Well biostratigraphic events are integrated with reflection seismic in the area to correlate regional seismic stratigraphic surfaces. The four wells are: Bandol-1, Emerillon C-56, East Wolverine G-37, and Heron H-73. We extend the event stratigraphic scheme previously developed for the Scotian Margin, offshore Nova Scotia, into new areas to the east along the southern Grand Banks, where we recognize four new well-log sequence stratigraphic events, and we modify the definition of a previously recognized regional surface. The new and modified regional surfaces are the Early Albian Unconformity, the Late Bathonian Maximum Flooding Surface (MFS), the Late Bajocian MFS (renamed from Bathonian/Bajocian MFS), the ?Bajocian/Toarcian Unconformity, and the Late Pliensbachian MFS. We recognize the "Avalon Unconformity" and "Base-Tertiary Unconformity" of previous studies as amalgamations of multiple smaller-scale unconformities and refine their age in the studied wells. A major improvement over our earlier Scotian Margin event schemes is the extension of the event stratigraphy into the Early Jurassic using a suite of marine biostratigraphic markers. We compare the Early Jurassic event scheme to Deep Sea Drilling Project Site 547B on the conjugate Moroccan Margin to better constrain potential source rock intervals and the early history of the central Atlantic Ocean.

Shale Gas Reservoir Pore Pressure Prediction: A Case Study of the Wufeng–Longmaxi Formations in Sichuan Basin, Southwest China

Pore pressure prediction is critical for shale gas reservoir characterization and simulation. The Wufeng–Longmaxi shale, in the southeastern margin of the Sichuan Basin, is identified as a complex reservoir affected by overpressure generation mechanisms and variability in lithification. Thus, standard methods need to be adapted to consistently evaluate pore pressure in this basin. Based on wireline logs, formation pressure tests, and geological data, this study applied the Eaton–Yale approach, which extends the theoretical basis of Eaton and Bowers methods to reservoir geological conditions and basin history. The method was developed by integrating petrophysical properties, rock physics interpretations, and geology information. The essential steps include (1) a multi-mineral analysis to determine mineral and fluid volumes; (2) a determination of the normal pressure trend line and extending it to overpressured sections; (3) predicting pore pressure using the basic Eaton approach and identifying overpressured zones; (4) correcting compressional velocity using lithology logs and a rock physics model; (5) determining the Biot Alpha coefficient and vertical-effective stress and estimating the new pore pressure values using the Eaton–Yale method. Overpressure zones were corrected, and reservoir pore pressure varied between 30.354 and 34.959 MPa in the wells. These research results can provide a basis for building reservoir simulation models, identifying reservoir boundaries, and predicting relative permeability.

Multi-scale data joint inversion of minerals and porosity in altered igneous reservoirs—A case study in the South China Sea

There are abundant igneous gas reservoirs in the South China Sea with significant value of research, and lithology classification, mineral analysis and porosity inversion are important links in reservoir evaluation. However, affected by the diverse lithology, complicated mineral and widespread alteration, conventional logging lithology classification and mineral inversion become considerably difficult. At the same time, owing to the limitation of the wireline log response equation, the quantity and accuracy of minerals can hardly meet the exploration requirements of igneous formations. To overcome those issues, this study takes the South China Sea as an example, and combines multi-scale data such as micro rock slices, petrophysical experiments, wireline log and element cutting log to establish a set of joint inversion methods for minerals and porosity of altered igneous rocks. Specifically, we define the lithology and mineral characteristics through core slices and mineral data, and establish an igneous multi-mineral volumetric model. Then we determine element cutting log correction method based on core element data, and combine wireline log and corrected element cutting log to perform the lithology classification and joint inversion of minerals and porosity. However, it is always difficult to determine the elemental eigenvalues of different minerals in inversion. This paper uses multiple linear regression methods to solve this problem. Finally, an integrated inversion technique for altered igneous formations was developed. The results show that the corrected element cutting log are in good agreement with the core element data, and the mineral and porosity results obtained from the joint inversion based on the wireline log and corrected element cutting log are also in good agreement with the core data from X-ray diffraction. The results demonstrate that the inversion technique is applicable and this study provides a new direction for the mineral inversion research of altered igneous formations.

Application of audio-frequency magnetotelluric data to cover characterisation – validation against borehole petrophysics in the East Tennant region, Northern Australia

The characterisation of the thickness and geology of cover sequences significantly improves targeting for mineral exploration in buried terrains. Audio-frequency Magnetotelluric (AMT) data is applicable to characterise cover sequences, where their conductivity (inverse resistivity) can be differentiated. We present a regional study from the under-cover East Tennant region in the Northern Territory (Australia) where we have applied deterministic and probabilistic inversion methods to derive 2D and 1D resistivity models. We integrated these models with information of co-located basement penetrating boreholes (lithological and geophysical logs) to ground-truth and validate the models and to improve geophysical interpretations. In the East Tennant region, borehole lithology and wireline logging demonstrate that the modelled AMT response is largely controlled by the mineralogy of the cover and basement rocks. The bulk conductivity is due primarily to bulk mineralogy and the success of using the AMT models to predict cover thickness is shown to be dependent on whether the bulk mineralogy of cover and basement rocks are sufficiently different to provide a detectable conductivity contrast. Our investigation of a range of geological scenarios that differ in thickness, complexity and geology of the cover and basement rocks suggests that in areas where there is sufficient difference in bulk mineralogy and where the stratigraphy is relatively simple, AMT models predict the cover thickness with high certainty. In more complex scenarios interpretation of AMT models may be more ambiguous and requires integration with other data (e.g. drilling, wireline logging, potential field modelling). Overall, we conclude that the application of the method has been validated and the results compare favourably with borehole stratigraphy logs once geological (i.e. bulk mineralogical) complexity is understood. This demonstrates that the method is capable of identifying major litho-stratigraphic units with resistivity contrasts. Our results have assisted with the planning of regional drilling programs and have helped to reduce the uncertainty and risk associated with intersecting targeted stratigraphic units in covered terrains.