ALBERT

Description: Beaches are characterized by complex spatial and temporal patterns of erosion and accretion subjected to significant wave and tide influence. The objective of this study is to estimate the evolution of hydromorphodynamic conditions on the shoreline of Grand Popo Beach observed from two adjacent video camera setups. We have analyzed the impact of the variability of hydrodynamic parameters on the beach evolution and evaluated the variabilities of the hydrodynamic and morphologic parameters from the two cameras. Despite the nonhomogeneity within the cameras’ intrinsic properties, the various results obtained from the two systems indicate that wave conditions (peak period and significant height) from the cameras have the same variations, whereas the shoreline variations of camera A are not the same as those of camera B. It is generally during the summer that the Grand Popo Beach is exposed to an agitated environment with strong observed values of significant wave heights and wave energy flux, undoubtedly resulting in significant sediment transport along the beach leading a shoreline retreat. The results indicate that in 3.5 years the shoreline of Grand Popo Beach has retreated by 10 m.

Description: Coastal erosion has become an issue globally. As the sea level continues to rise due to global warming, projections are that erosion would increase. To address the issue sustainably, relevant scientific information such as sediment transport and shoreline dynamics is required at local scales. Over the years, erosion at the eastern coast of Ghana has been mainly evaluated using low-resolution imagery due to the challenges in retrieving high-resolution data in the nearshore region. An unmanned aerial vehicle (UAV) approach was used to assess beach evolution on a seasonal and short-term basis between May 2018 and December 2019. Using the UAV approach, shoreline changes, sediment volume changes, and profile elevation parameters were extracted for the Dzita beach of Ghana. Shoreline changes indicated the dominance of erosion during the first and third phases (May 2018–December 2018 and June 2019–December 2019, respectively) at a rate of −7.23 ± 0.23 and −4.85 ± 0.23 m/yr, whereas the second phase (January 2019–June 2019) showed accretion of +8.44 ± 0.23 m/yr. Beach profiles from the first, second, and third phases had steep and gentle slopes, respectively. From these observations, it was recommended that soft engineering approaches such as beach nourishment should be implemented to protect the shoreline and strict prevention of nearshore sand mining and gravel mining. It is also possible that the beach could go through a cycle of changes. Further studies using this same approach should be done as well as probing into other parameters such as nearshore bathymetry to have a better understanding of beach dynamics as envisaged.

Description: The connectivity of complex carbonate reservoirs has an essential impact on the exploration and development of these reservoirs. From geologic genesis, the connectivity of complex carbonate reservoirs is mainly controlled by faults and dissolution. Therefore, accurate identification of faults and karst caves is the key to studying reservoir connectivity. The Ordovician carbonate reservoir in the Hudson Oilfield of the Tarim Basin is used for our reservoir connectivity analysis study. First, we calculate the coherence and curvature attributes, respectively, and then merge the two attributes using a neural network algorithm. Finally, we use the ant-tracking method to track the faults for the merged data. The results show that the approach substantially enhances deterministic faults that can be seen directly on the seismic data, and subtle faults can also be identified. For reservoir identification, we use the diffraction imaging method to describe the karst reservoir in this study area. The results show that diffraction imaging can identify small-scale caves that cannot be well recognized on the seismic reflection data. Furthermore, the caves connected on the diffraction seismic data are isolated from each other on the seismic reflection data, making the connection between caves clearer. Based on the results of the fault and cave identification, we analyze the reservoir connectivity of the study area using the oil pressure and daily production data, which indicates that the north–northwest and near-north–south faults probably play a role in the connection of the reservoirs, whereas the northeast–east faults tend to block the connection of the reservoirs.

Description: The choice of an initial model for seismic waveform inversion is important. In matured exploration areas with adequate well control, we can generate a suitable initial model using well information. However, in new areas where well control is sparse or unavailable, such an initial model is compromised and/or biased by the regions with more well controls. Even in matured exploration areas, if we use time-lapse seismic data to predict dynamic reservoir properties, an initial model that we obtain from the existing preproduction wells could be incorrect. In this work, we outline a new methodology and workflow for a nonlinear prestack isotropic elastic waveform inversion. We call this method a data-driven inversion, meaning that we derive the initial model entirely from the seismic data without using any well information. By assuming a locally horizontal stratification for every common midpoint and starting from the interval P-wave velocity, estimated entirely from seismic data, our method generates pseudowells by running a two-pass 1D isotropic elastic prestack-waveform inversion that uses the reflectivity method for forward modeling and the genetic algorithm for optimization. We then use the estimated pseudowells to build the initial model for seismic inversion. By applying this methodology to real seismic data from two different geologic settings, we determine the usefulness of our method. We believe that our new method is potentially applicable for subsurface characterization in areas where well information is sparse or unavailable. Additional research is, however, necessary to improve the computational efficiency of the methodology.

Description: Lateral changes in velocity about faults can give rise to fault shadow (FS) zones on time-migrated data volumes, which can result in structural interpretation artifacts in the fault trap reservoir. To address this issue, we have adopted a new reconstruction method of FS distortion structures based on a deep learning fully connected network (FCN). We use the 3D stratigraphic dip attributes to quantitatively delineate the extent of the FS zone. Then, we train a model to construct a nonlinear trend surface based on the structures of the stratigraphic reflectors that fall outside of the shadow zone. Finally, we use this nonlinear trend surface to compensate for the distorted structure within the FS zone. We calibrate our method using synthetic data and find that the method can accurately recover structural data within the FS distortion zone. We then test the effectiveness of our workflow by applying it to recover real FS distortion structures in the Pearl River Mouth Basin of the South China Sea. The results confirm that our method significantly reduces drilling depth errors in the FS zone. Compared to the traditional polynomial fitting method, the multilayer, multiparameter, and flexible nonlinear activation function of FCN is more capable of reconstructing nonlinear geologic structures in the FS zone. We find the FCN-based geologic reconstruction method to be efficient and effective for exploring potential structures in the FS zone and thereby in avoiding the risks of structural failure.

Description: Enhanced hydrocarbon recovery is essential for continued economic development of unconventional reservoirs. We have focused on dynamic characterization of the Niobrara and Codell Formations in Wattenberg Field through the development and analysis of a full integrated reservoir model. We determine the effectiveness of the hydraulic fracturing and production with two seismic monitor surveys, surface microseismic, completion data, and production data. The two monitor surveys were recorded after stimulation and again after two years of production. Identification of reservoir deformation due to hydraulic fracturing and production improves reservoir models by mapping nonstimulated and nonproducing zones. Monitoring these time-variant changes improves the prediction capability of reservoir models, which in turn leads to improved well and stage placement. We quantify dynamic reservoir changes with time-lapse P-wave seismic data using prestack inversion and velocity-independent layer stripping for velocity and attenuation changes within the Niobrara and Codell reservoirs. A 3D geomechanical model and production data are history matched, and a simulation is run for two years of production. Results are integrated with time-lapse seismic data to illustrate the effects of hydraulic fracturing and production. Our analyses illustrate that chalk facies have significantly higher hydraulic fracture efficiency and production performance than marl facies. In addition, structural and hydraulic complexity associated with faults generate spatial variability in a well’s total production.

Description: Compaction effects can obscure the impedance separation between hydrocarbon-bearing and fully brine-saturated sandstones. We have improved their discrimination by removing depth-related trends from inverted seismic impedance. Although the ratio of compressional- to shear-wave velocity versus seismic compressional-wave impedance crossplots shows differences among pay, brine sand, and shale trends, using absolute inverted impedances only imperfectly distinguishes hydrocarbon sands from brine sands due to outliers. In a given locality, statistical comparison of well log and seismic-derived impedances enables us to obtain a shale impedance model for a lithology baseline to detrend the impedance from the effects of burial and overburden. This has the effect of unmasking anomalies associated with hydrocarbon-bearing sands and serves as a reliable fluid discriminator. For an offshore Gulf of Mexico data set on the flank of a salt dome, with pay occurring over a wide range of depths, we identify hydrocarbon-bearing sands with a greater success rate after detrending the absolute seismic impedance.

Description: The best petrophysical models are based on direct measurements from the core. Unfortunately, core is not available in many cases, either for economic, logistical, or historical reasons. In this study, we needed to construct a detailed Field Development Plan (FDP) for the small, marginal B-9 field in the Western Offshore Basin, India, which did not merit core acquisition. The objective is to propose a workflow for building a petrophysical model with limited data sets instead of a typical FDP workflow. After analyzing the assumptions, limitations, and uncertainties involved in the petrophysical model, we used advanced petrophysical logs to reduce uncertainty and create a robust petrophysical model. We carried out a log-based petrophysical study to determine the volume of shale, porosity, saturation, and permeability. The advanced petrophysical logs including spectroscopy, nuclear magnetic resonance (NMR), formation pressures, and well testing data are used to calibrate the petrophysical model. Spectroscopy data are used to calibrate the mineralogical volumes and grain density, whereas the porosity is calibrated from NMR data. We calibrated log-derived permeability results with NMR permeability and mobility from well test data. We used heterogeneous rock analysis on petrophysical outputs to carry out petrophysical rock typing (PRT). This has helped in establishing the porosity-permeability relationship and saturation-height model for each PRT. In the absence of irreducible water saturation ([Formula: see text]) information from the core, NMR-derived [Formula: see text] is calculated and then used to calibrate the saturation model. Log-derived permeability and saturation are estimated, which agrees well with the available testing data. This provided a robust petrophysical model that served as a basis for geologic static and reservoir dynamic models. The gas-down-to and water-up-to methods are used to establish the contacts. The resulting saturation height model agreed well with the saturations derived from the log, which gave us confidence in our dynamic model.

Description: The complete characterization of a reservoir requires accurate determination of properties such as the porosity, gamma ray, and density, among others. A common workflow is to predict the spatial distribution of properties measured by well logs to those that can be computed from the seismic data. In general, a high degree of scatter of data points is seen on crossplots between P-impedance and porosity, or P-impedance and gamma ray, suggesting great uncertainty in the determined relationship. Although for many rocks there is a well-established petrophysical model correlating the P-impedance to porosity, there is not a comparable model correlating the P-impedance to gamma ray. To address this issue, interpreters can use crossplots to graphically correlate two seismically derived variables to well measurements plotted in color. When there are more than two seismically derived variables, the interpreter can use multilinear regression or artificial neural network analysis that uses a percentage of the upscaled well data for training to establish an empirical relation with the input seismic data and then uses the remaining well data to validate the relationship. Once validated at the wells, this relationship can then be used to predict the desired reservoir property volumetrically. We have described the application of deep neural network (DNN) analysis for the determination of porosity and gamma ray over the Volve field in the southern Norwegian North Sea. After using several quality-control steps in the DNN workflow and observing encouraging results, we validate the final prediction of the porosity and gamma-ray properties using blind well correlation. The application of this workflow promises significant improvement to the reservoir property determination for fields that have good well control and exhibit lateral variations in the sought properties.

Description: Erosion is a major global challenge facing coastal regions, and it is projected to increase on a regional to global scale as sea levels continue to rise. Deltas, which are important ecosystems, are particularly vulnerable due to their low-lying nature, subsidence, reduction in sediment supply, population increase, and exposure to an increasing frequency of extreme events. The need for sustainable management of these systems requires accurate estimates of shoreline dynamics at the local scale and higher spatial resolutions for engineering and decision making. We have assessed the shoreline dynamics of the Volta River delta in Ghana for a medium term of 12 years using high-resolution satellite imagery. The shoreline change rates are correlated with wave dynamics to explain the observed shoreline evolution within the delta. Our results confirm that erosion dominates the studied coasts with rates reaching as high as 31 m/yr close to the mouth of the delta where the water level shows a strong relationship with the shoreline change. These rates are evidenced by the destruction of fishing villages such as Fuveme located close to the mouth. Anthropogenic factors, such as the construction of sea defense projects, are also influencing erosion patterns across the study area. We recommend a softer approach for coastal management within the delta.

Description: Recent advancements in distributed acoustic sensing (DAS) technology open new ways for borehole-based seismic monitoring of CO2 geosequestration. Compared to 4D surface seismic monitoring, repeated vertical seismic profiling (VSP) surveys with DAS receivers considerably reduce the cost and invasiveness of time-lapse CO2 monitoring. However, standard borehole imaging techniques cannot provide the same level of reservoir illumination as 3D surface seismic. The performance of VSP imaging can be significantly improved with interferometric utilization of free-surface multiples. We have developed a feasibility study of interferometric imaging with a synthetic walkaway VSP data set, followed by its application to field walkaway VSP data recorded by conventional borehole geophones and two types of DAS (standard and engineered fibers). Both experiments (synthetic and field) demonstrate that interferometric imaging is a viable method to extend the subsurface image beyond the coverage of standard VSP imaging. Specifically, the interferometry approach provides a more detailed upper section of the subsurface, whereas standard migration of primary reflections provides a more detailed bottom part of the image. Comparison of the standard and engineered fibers indicates that both fibers are sensitive to free-surface multiples, but the engineered fiber provides a much higher signal-to-noise ratio; thus, it is preferable for interferometric imaging with multiples. The result obtained with the engineered DAS cable indicates that in the depth range suitable for both methods, the VSP interferometric image of the reflectors is comparable to the surface seismic image. The experiment on the field DAS data proves that DAS is sensitive enough to record the nonprimary wavefield for imaging and monitoring of the subsurface.

Description: Distributed acoustic sensing (DAS) technologies are now becoming widespread, particularly in vertical seismic profiling (VSP). Being a spatially densely sampled recording of the seismic wavefield, DAS data provide an extended measurement compared with point geophone VSP. We have developed a basic theory that enables an intuitive geophysical understanding of DAS data amplitudes using the concepts of kinetic and potential energy and their fluxes. We start by relating DAS and geophone measurements to potential energy and kinetic energy, respectively. We use this relationship and energy balancing along the well to construct a scheme for inverting DAS and geophone wavefields for density and velocity simultaneously. Then, recognizing that it may be impractical to have geophones and DAS, we adopt a second inversion scheme that eliminates the need for geophones and uses up- and downgoing DAS wavefields instead. There is no need for first-break picking or windowing the data, and the full-length DAS records can be used in both inversion schemes. We test these inversion schemes on 2D elastic synthetics.

Description: Seismic amplitude-variation-with-offset (AVO) inversion from prestack seismic data plays a significant role in estimating elastic parameters and characterizing reservoir properties. In general, sparse regularization is widely used to solve ill-posed inverse problems by reducing the solution space of subsurface parameters, which makes seismic AVO inversion more stable. However, the traditional sparse constraint inversion only focuses on the vector sparsity of reflectivity, instead of the structural sparse characteristics of the estimated parameters. Consequently, various elastic parameters demonstrate different formation structural features in the same location of stratum. In this study, we have developed a novel approach that combines the structural sparsity and the vector sparsity of the model reflectivity to establish the posterior probability density distribution and solve the objective function of the model parameters. Based on the relationship among multiple elastic parameters, we divide the model parameters to be inverted into several groups according to intrinsic structural sparse characteristics of elastic parameters. In this case, all of the model parameters at the same sampling point are classified into the identical group, which ensures that different estimated parameters indicate the same characteristic in terms of stratigraphic structure. From the perspective of Bayesian inference, we use the modified Cauchy probability density function (PDF) to characterize the group sparsity and describe the relationship among model parameters in the same group by Gaussian PDF. Furthermore, we estimate the optimum solution corresponding to the maximum a posteriori probability under Bayesian inference. Synthetic experiments on a Marmousi model prove that the estimated P-velocity, S-velocity, and density are consistent with those of the real models, and the application of field data confirms the availability and feasibility of group sparse inversion.

Description: Mississippian carbonate and silica-rich reservoirs of northern and central Oklahoma formed along a regionally extensive carbonate ramp to basin transect. The stratigraphy, lithology, and porosity characteristics of the Mississippian Meramec and Osage series significantly vary as older ramp carbonates prograde southward and transition into younger calcareous and quartz-rich silt deposits of the Anadarko Basin. Lithofacies identified within the northern carbonate-dominated portion of the system commonly include altered chert, skeletal grainstones, peloidal packstones-grainstones, bioturbated wackestones-packstones, bioturbated mudstones-wackestones, glauconitic sandstones, and siliceous shale. Lithofacies within the southern siliciclastic-dominated portion of the system include structureless to bioturbated sandstones, siltstones, and laminated mudstones, each with varying degrees of carbonate content. We have grouped these core-based lithofacies into dominant lithologies/rock types, which tie to well-log properties. Electrofacies classification methods including the artificial neural network (ANN) and k-means clustering predict lithologies in noncored wells. ANN yielded the highest overall prediction accuracy of 85% for lithologies. Core, well log, and lithology log data establish the regional stratigraphic framework. The Mississippian interval of interest subdivides into 16 stratigraphic zones. A depositional-dip-oriented cross section and the associated reservoir models illustrate the proximal to distal and stratigraphic variability of the lithology and porosity. Lithology trends moving from north to south, from older to younger strata, reveal a carbonate-dominated succession capped by diagenetically altered chert northward shifting into a siliciclastic-dominated interval, which increases in clay content southward. Northward, prospective conventional reservoirs developed near cycle tops within diagenetically replaced cherts and cherty limestones associated with subaerial exposure and sea-level fluctuations. Southward, higher total porosity associates with increased clay content linked to the suppression of calcite cement, forming prospective unconventional targets near the bases of depositional cycles.

Description: Tight-gas sandstone reservoirs of the Ordos Basin in China are characterized by high rock-fragment content, dissimilar pore types, and a random distribution of fluids, leading to strong local heterogeneity. We model the seismic properties of these sandstones with the double-double porosity theory, which considers water saturation, porosity, and the frame characteristics. A generalized seismic wavelet is used to fit the real wavelet, and the peak frequency-shift method is combined with the generalized S-transform to estimate attenuation. Then, we establish rock-physics templates (RPTs) based on P-wave attenuation and impedance. We use the log data and related seismic traces to calibrate the RPTs and generate a 3D volume of rock-physics attributes for the quantitative prediction of saturation and porosity. The predicted values are in good agreement with the actual gas production reports, indicating that the method can be effectively applied to heterogeneous tight-gas sandstone reservoirs.

Description: The second member of the Kongdian Formation (usually abbreviated as the E k2 shale) is one of the most significant exploring targets for shale oil at the Cangdong Sag of the central Bohai Bay Basin. It consists of siliceous shale, mixed shale, and calcareous shale. To better understand why organic matter accumulated in the E k2 shale, we have analyzed major and trace elemental compositions to reconstruct the provenance and sedimentary environment. Tectonic discriminatory diagrams suggest that the tectonic setting of the parental rocks for the E k2 shale belonged to the Continental Island Arc. The distribution patterns of trace elements and rare earth elements + yttrium (REEs + Y) are close to the intermediate igneous rock. The ratios of Al2O3/TiO2 ranging from 21.41 to 27.59 with a mean value of 23.93 also demonstrate a parental rock of the intermediate igneous rock. Siliceous and mixed shales indicate K2O/Al2O3 of 0.17–0.29, chemical index of weathering of 28.79–97.79, plagioclase index of alteration of 38.24–95.57, and chemical index of alteration of 40.29–80.23. These weathering proxies denote that the E k2 shale underwent a low weathering degree in an arid climate and a high weathering degree in a semiarid climate. The V/(V + Ni) ratios and pyrite framboids indicate an anoxic sedimentary condition. The δ18O values of carbonate minerals in the E k2 shale range from −9.8‰ to 0.7‰, and they are positively correlated to the δ13C values. The Sr/Ba ratios, δ18O, and chemical mineral associations indicate that siliceous and mixed shales were deposited in a fresh to brackish anoxic water column under a semiarid climate. Whereas calcareous shale was deposited in a saline to hypersaline anoxic water column under an arid climate.

Description: Hydraulic fracturing (HF) and horizontal drilling are essential to the development of shale gas and oil. Production depends on the stimulation success. During fracture initiation, propagation, and closure, cracks emit acoustic waves; these can be monitored in real time as microseismics in the field and as acoustic emissions (AEs) in the laboratory. AEs are the laboratory equivalent of field-scale microseismics and contain detailed information about HF fracture mechanics. The number of acoustic events correlates with the number of induced fractures and hence the stimulation volume. Three HF protocols under dry conditions were carried out on Tennessee sandstone: (1) a constant injection rate, (2) a precyclic injection, and (3) a variable-rate injection test. All three tests were performed under the same principal stress conditions: vertical stress of 10.3 MPa (1500 psi), minimum horizontal stress of 3.5 MPa (500 psi), and maximum horizontal stress of 20.7 MPa (3000 psi). In total, 16 piezoelectric transducers were mounted around a cylindrical sample to record the AEs. We have performed postsignal processing to extract AE event attributes, including the amplitudes, signal-to-noise ratio, arrival time, event location (with the velocity-anisotropy input), and frequency analyses. The AE events associated with the constant-rate injection test possessed the lowest frequencies (150–270 kHz). The variable-rate test AE events possessed higher frequencies (160–310 kHz), whereas the precyclic injection had events with the highest frequencies, peaking at 330 kHz. Acoustic events before failure had lower amplitudes, but higher frequency compared to those recorded postbreakdown, suggesting different failure modes. Precyclic injection induced the greatest number of locatable events before and after failure.

Description: A stage of mafic magmatic activity occurred in Early Permian in the Dabasong Uplift of Junggar Basin, part of the magma intruded into the normal sedimentary and shallow buried fine sandstone to form diabase, and part of the magma erupted to form basalt. The surrounding fine sandstone just entered in the early diagenetic stage A when the magma intruded. The compaction of the surrounding clastic rock and rupture of a small number of clastic grains were caused by the extrusion of the magma intrusion. The presence of chemically deposited alkaline minerals such as calcite, dolomite, shortite, natural alkali, and northupite indicates an alkali lake sedimentary environment for the Fengcheng Formation. Primary alkaline minerals dissolved from the surrounding rocks were subsequently transported and precipitated to form cements. The formation of the calcite cements and calcite metasomatism resulted in considerable densification of the surrounding rock during early diagenesis and destruction of the reservoir quality. The mafic magma had abundant Fe2+ and Mg2+ ions and was deficient in K+ ions, resulting in large amounts of chlorite and iron precipitation in the surrounding rock mainly composed of clay. We have analyzed the influence of an ultrashallow intrusion on the surrounding clastic rock during the early diagenetic period, which provided a typical reference for establishing a systematic mechanistic model of how magmatic intrusions affect the surrounding rock.

Description: Integrated airborne and ground geophysical studies were conducted in parts of Ilesha schist belt, southwestern Nigeria. The goal was to provide a useful guide for mineral prospecting, with the hope of considerably narrowing down the future search for mineral deposits within the study area. Aeromagnetic and aeroradiometric data were analyzed for the reconnaissance study. In addition, the reduction-to-equator transform, analytic signal, tilt derivative, and Euler deconvolution filters were applied to the aeromagnetic data to enhance shallow and deep geologic features. The aeroradiometric data were used to determine spatial variations in the concentrations of uranium (U), thorium (Th), and potassium (K) in near-surface rocks and to map spatial lithologic changes. The 2D-magnetic sections, radiometric profiles, inverted resistivity, and induced polarization (IP) sections were generated from the integrated geophysical data. The electrical resistivity tomography (ERT) results reveal the subsurface heterogeneity (to a depth of approximately 197 m) and varied geoelectric layers (topsoil, lateritic-clay, weathered rock, and basement rock). The IP sections show varying degrees of chargeability and features that suggest the presence of disseminated mineralized bodies concealed in some areas. The overburden thickness varies between 4 and 85 m as determined from the 2D-magnetic and electric resistivity sections. Anomalous peaks on profiles of elemental ratios (eTh/K, eTh/eU, and K/eU) correlate with the results of IP and ERT. Data sets are well correlated and highlight areas with relevant structural and lithologic signatures favorable for mineral deposition. The methodology adopted in our research is well adapted, and the interpretation techniques provided insight into regional and local lithostructural settings. These anomalous areas are suggested as targets for future exploration works.

Description: We have evaluated the results of a receiver decimation study in a deepwater context using separated wavefield imaging (SWIM) algorithms to provide extended illumination for imaging without ocean-bottom node (OBN) positioning constraints. We carried out subsurface imaging using the SWIM imaging technique with a reduced OBN layout, and we compared the results with those from conventional one-way wave-equation migration. We found from the results that the SWIM algorithm makes it possible to reduce the OBN layout while obtaining a similar subsurface image with the same shot geometry, which allows a reduced receiver acquisition effort, offers more geometry flexibility without affecting the image quality, with a potentially significant reduction of acquisition cost and 4D processing turnaround time.

Description: With the rapid development of the high-speed railway industry, train detection and identification play a vital role in capacity improvement and safe operation in railway systems. Conventional detection methods such as track circuit and axle counting tend to be interfered with by severe weather conditions and irrelevant conductive objects, leading to false detections. Fiber-optic distributed acoustic sensing (DAS) technology is a prevailing sensing method in geophysics research, petroleum exploration, and structure inspection. Compared to traditional detection techniques, DAS is suitable for long-distance detection and is resistant to severe weather conditions and electrical interference. We have developed a train detection and classification system using DAS technology and have explored an effective classification method for train identification. Specifically, we conduct a field experiment by the side of a railroad over viaducts and the data are collected with the DAS detection system. To eliminate the impact of background noise, DC noise, and motor vehicle signals from the original data, we adopt a wavelet denoising method and Chebyshev filter to extract the features of three types of train signals. The vibration signals of these different trains indicate remarkable cyclical variations related to the number of wheelsets in the time domain and have similar narrow-band discrete spectrums with different characteristic peak frequencies. Furthermore, based on the features of the train signals, we select a support vector machine classifier to identify three types of trains, with accuracy greater than 97%.

Description: Several P-wave azimuthal anisotropy studies have been conducted for the SEG Advanced Modeling Program Phase II Barrett model data. However, these analyses provide fracture property estimation that is inconsistent with the actual model properties. Therefore, we perform a feasibility study to understand the influence of the overburden and reservoir properties, and processing and inversion steps, which together determine the success of the fracture interpretation from seismic data. The 1D model properties (orthorhombic for overburden and reservoir) are first extracted from the actual Barrett model properties at two locations. Anisotropic prestack reflectivity modeling exposes the true orthorhombic response of the 1D medium in the form of common-offset and common-azimuth gathers. The true anisotropic response is obscured in the Barrett data (generated by finite-element modeling) due to the mild lateral velocity variations and orthorhombic anisotropy in the overburden. We then expose the reservoir anisotropic response using an isotropic overburden in the reflectivity modeling. This indicates that the P-wave velocity variation with azimuth (VVAZ) responses generated by the reservoir itself are weak, which leads to an unstable VVAZ inversion to estimate the interval normal moveout (NMO) velocity anisotropy. The reservoir thickness (125 m or 65 ms of two-way traveltime) or NMO velocity anisotropy (6%–7%) needs to be at least doubled to obtain a stable VVAZ inversion. Anisotropic geometric-spreading correction improves the amplitude variation with azimuth inversion results when reflectivity modeling models the orthorhombic overburden. The converted wave (C-wave) has a stronger VVAZ response compared to that of the P-wave. We suggest that the C-wave data could be useful to constrain fracture interpretation in the Barrett model. We conclude that the results of previous studies are due to the combination of the residual influence of overburden after processing and imaging and the weak anisotropy responses from the reservoir.

Description: Pockmarks, as depression morphology related to fluid escape on the seafloor, are revealed by 3D seismic data on the northwestern South China Sea (SCS) margin. The pockmarks can be classified into two groups based on their various shapes in plan view: the circular group and the elongating group. These pockmarks in the study area could be defined as mega-pockmarks because their maximum diameters can reach to 7.5 km. They commonly develop more than one crater, which are called the central crater and the secondary crater. The seismic data illuminated their complicated internal architectures in the subsurface, as well as their evolution periods, such as the initiation stage, mature stage, and abandonment stage. According to the buried structures and their genesis mechanism, mega-pockmarks could be classified into linear faults-associated pockmarks and volcano-associated pockmarks. The linear fault-associated pockmarks root on the top Middle Miocene, where the linear faults are distributed. The linear faults on the top of fluid reservoir in Middle Miocene act as conduits for fluid seepage. The fluid seepage is driven by the break of balance between the hydrostatic and pore pressure. When the fluid seepage initiates, they will migrate along the linear faults, making the linear feature of pockmarks on the seafloor. Thermogenic gas from deep intervals and biogenic gas from shallow intervals may be fluid sources for the genesis of pockmarks. However, the volcanic activities control the genesis and evolution of volcano-associated pockmarks. Volcano-associated pockmarks root on the craters of volcanoes. The volcanoes underneath the pockmarks provide volcanic hydrothermal solutions, such as phreatomagmatic eruptions through the volcanic craters. The confined fluid seepages make the pockmarks exhibit a more circular shape on the seafloor. Long-term, multiepisode fluid expulsions generate the complicated internal architecture that leads to multicratered mega-pockmarks on the northwestern margin of SCS.

Description: The Middle to Upper Eocene series are characterized by multiple hiatuses related to erosion, nondeposition, or condensed series in the Cap Bon and Gulf of Hammamet provinces. We have performed an integrated study taking advantage of surface and subsurface geology, faunal content, borehole logs, electrical well logs, vertical seismic profiles, and surface seismic sections. Calibrated seismic profiles together with borehole data analysis reveal unconformities with deep erosion, pinch outs, normal faulting, and basin inversion that are dated Campanian, intra-Lutetian, and Priabonian compressive phases; these events were also described at the regional scale in Tunisia. Tectonics, sea-level fluctuations, and climate changes closely controlled the depositional process during the Middle to Upper Eocene time. The depositional environment ranges from internal to outer platform separated by an inherited paleo high. We determine eight third-order sequences characterizing the interaction among tectonic pulsations, sea-level changes, and the developed accommodation space within the Middle to Upper Eocene interval. We correlate the obtained results of the Cap Bon-Gulf of Hammamet provinces with the published global charts of sea-level changes, and we find a good correspondence across third-order cycles. Model-based 3D inversion proved to be a solution to model the lateral and vertical lithologic distribution of the Middle to Upper Eocene series.

Description: Data transformation, regional-residual separation, trend analysis, and analytic signal (AS) depth estimation were applied to aeromagnetic data covering the Anambra Basin, which is a major depocenter in the Benue Trough, southeast Nigeria, with the primary objectives of accentuating the attributes of magnetic sources and determining if sufficient sediment thickness exists for hydrocarbon generation, maturation, and expulsion. The application of data transformation techniques (such as map projection, merging, and reduction-to-pole) and regional-residual ensured the computation of a crustal magnetic field that would be suitable for magnetic analyses. Results indicate that the magnetic basement in the basin forms an undulating surface overlain by sediments with an average thickness ranging between 4 and 7.5 km, whereas the maximum thickness reaches 8 km in some areas. This depth range suggests a promising prospect for source-facies maturation and expulsion. We expect that areas in the study area with these appreciable sediment thicknesses, good preservation of graben fill, and suitable areal closures or fault structures would be favorable for hydrocarbon prospectivity.

Description: We have carried out a geomechanical study of three wells, one each in the Niobrara A, Niobrara C, and Codell Sandstone to investigate how the state of stress and stress variations with depth affect vertical hydraulic fracture growth and shear stimulation of preexisting fractures. We determine that the higher magnitudes of measured least principal stress values in the Niobrara A and C shales are the result of viscoplastic stress relaxation. Using a density log and a vertical transverse isotropy velocity model developed to accurately locate the microseismic events, we theoretically calculated a continuous profile of the magnitude of the least principal stress with depth. This stress profile explains the apparent vertical hydraulic fracture growth as inferred from the well-constrained depths of the associated microseismic events. Finally, we determine that because of the upward propagation of hydraulic fractures from the Niobrara C to the Niobrara A, the latter formation experienced considerably more shear stimulation, which may contribute to the greater production of oil and gas from that formation.

Description: The prediction of natural fracture networks and their geomechanical properties remains a challenge for unconventional reservoir characterization. Because natural fractures are highly heterogeneous and of subseismic scale, integrating petrophysical data (i.e., cores and well logs) with seismic data is important for building a reliable natural fracture model. Therefore, I have developed an integrated and stochastic approach for discrete fracture network modeling with field data experimentation. In the method, I first perform a seismic attribute analysis to highlight the discontinuity in the seismic data. Then, I extrapolate the well-log data that include localized but high-confidence information. By using the fracture intensity model including seismic and well logs, I build the final natural fracture model that can be used as a background model for the subsequent geomechanical analysis such as simulation of hydraulic fractures propagation. As a result, our workflow combining multiscale data in a stochastic approach constructs a reliable natural fracture model. I validate the constructed fracture distribution by its good agreement with the well-log data.

Description: The extensive development of unconventional reservoirs using horizontal drilling and multistage hydraulic fracturing has generated large volumes of reservoir characterization and production data. The analysis of this abundant data using statistical methods and advanced machine-learning (ML) techniques can provide data-driven insights into well performance. Most predictive modeling studies have focused on the impact that different well completion and stimulation strategies have on well production but have not fully exploited the available in situ rock property data to determine its role in reservoir productivity. We have used machine-learning techniques to rank rock mechanical properties, microseismic attributes, and stimulation parameters in the order of their significance for predicting natural gas production from an unconventional reservoir. The data for this study came from a hydraulically fractured well in the Marcellus Shale in Monongalia County, West Virginia. The data classes included measurements aggregated by well completion stage that included (1) gas production, (2) well-log-derived measurements including bulk density, elastic moduli, shear impedance, compressional impedance, brittleness, and gamma measurements, (3) microseismic attributes, (4) long-period long-duration (LPLD) event counts, (5) fracture counts, and (6) stimulation parameters that included the fluid injection volume and average pumping pressure. To identify observable proxies for the drivers of gas production, we evaluated five commonly used ML approaches including multivariate adaptive regression spline, Gaussian mixture model, random forest, gradient boosting, and neural network. We selected five variables including LPLD event count, seismogenic b-value, hydraulic diffusivity, cumulative moment, and fluid volume as the features most likely to impact gas productivity at the stage level in the study area. The data-driven selection of these parameters for their importance in determining gas production can help reservoir engineers design more effective hydraulic-fracture treatments in the Marcellus Shale and other similar unconventional reservoirs. Plain language summary: We use machine-learning methods and data-driven selection of reservoir parameters to rank and better understand their importance in determining gas production, which can help reservoir engineers design more effective hydraulic-fracture treatments in the Marcellus Shale and other similar unconventional reservoirs.

Description: Many applications in seismology involve the modeling of seismic-wave traveltimes in anisotropic media. We have developed homotopy solutions of the acoustic eikonal equation for P-wave traveltimes in attenuating transversely isotropic media with a vertical symmetry axis. Instead of the commonly used perturbation theory, we use the homotopy analysis method to express the traveltimes by a Taylor series expansion over powers of an embedding parameter. For the derivation, we first perform homotopy analysis of the eikonal equation and derive the linearized ordinary differential equations for the coefficients of the Taylor series expansion. Then, we obtain the homotopy solutions for the traveltimes by solving the linearized ordinary differential equations. Results of our investigation with approximate formulas demonstrate that the analytical expressions are efficient methods for the computation of traveltimes from the eikonal equation. In addition, these formulas are also effective methods for benchmarking approximate numerical solutions in strongly attenuating anisotropic media.

Description: The Triassic-Jurassic deep sandstone reservoirs in onshore Denmark are known geothermal targets that can be exploited for sustainable and green energy for the next several decades. The economic development of such resources requires accurate characterization of the sandstone reservoir properties, namely, volume of clay, porosity, and permeability. The classic approach to achieving such objectives has been to integrate well-log and prestack seismic data with geologic information to obtain facies and reservoir property predictions in a Bayesian framework. Using this prestack inversion approach, we can obtain superior spatial and temporal variations within the target formation. We then examined whether unsupervised facies classification in the target units can provide additional information. We evaluated several machine learning techniques and found that generative topographic mapping further subdivided intervals mapped by the Bayesian framework into additional subunits.

Description: In a shale gas reservoir, pore characterization is an important factor used to determine gas storage capacity. However, the nanometer (nm)-scale pore system in shale is difficult to explore by traditional optical, scanning electron microscopy, or even nuclear magnetic resonance well logging. We have investigated the pore structure and storage capacity of the Marcellus Shale through integration of petrophysical analysis from laboratory and well-logging data and nitrogen adsorption. The isotherm of Marcellus Shale is a composite isotherm, which has features of type I, type II, and type IV isotherms with type H4 of the hysteresis loop, suggesting slit-like pores developed in the Marcellus Shale. Quantitative analysis of pore volumes from the nitrogen adsorption indicates that density porosity may be more properly used to approximate the shale porosity and estimate the shale gas volume. In addition, the specific surface area, micropore, and mesopore volumes have a positive relationship with the kerogen volume and total organic content (TOC). By using the Langmuir and Brunauer-Emmet-Teller models, the simulated result indicates that the higher adsorbed quantity of the Marcellus Shale could be the result of the increase of micropore volume contributed, by the increase of kerogen or TOC content. The proposed equations rapidly compute TOC, a key parameter to predict gas storage capacity in overmature shale such as the Marcellus Shale.

Description: Reservoir damage is considered one of the major challenges in the oil and gas industry. Many studies have been conducted to understand formation damage mechanisms in borehole wells, but few studies have been conducted to analyze the data to detect the source, causes, and mitigations for each well where damage has occurred. I have investigated and quantified the reasons and mitigation of reservoir damage problems in the middle Miocene reservoir within the El Morgan oil field at the southern central Gulf of Suez, Egypt. I used integrated production, reservoir, and geologic data sets and their history during different operations to assess the reservoir damage in El Morgan-XX well. The collected data include the reservoir rock type, fluid, production, core analysis, rock mineralogy, geology, water chemistry, drilling fluids, perforations depth intervals, workover operations, and stimulation history. The integration of different sets of data gave a robust analysis of reservoir damage causes and helps to suggest suitable remediation. Based on these results, I conclude the following: (1) Workover fluid has been confirmed as the primary damage source, (2) the reservoir damage mechanisms could be generated by multisources including solids and filtrate invasions, fluid/rock interaction (deflocculating of kaolinite clay), water blockage, salinity chock, and the high sulfate content of the invaded fluid, and (3) multidata integration leads to appropriate reservoir damage analysis and effective design of the stimulation treatment. Furthermore, minimizing fluid invasion into the reservoir section by managing the overbalance during drilling and workover operations could be very helpful. Fluid types and solids should be considered when designing the stimulation treatment and compatibility tests should be performed. Long periods of completion fluid in boreholes are not recommended, particularly if the completion fluid pressure and reservoir pressure are out of balance, as well as the presence of sensitive formation minerals.

Description: We have characterized the petrophysical and geomechanical properties of the Late Cretaceous Turonian and Cenomanian carbonate reservoirs from the southeast Constantine Basin, northern Algeria. In general, Turonian carbonates exhibit a wide range of porosities (2%–15%) and permeabilities (0.001–10 mD), whereas the Cenomanian reservoir appears to be very tight (

Description: In the Wattenberg Field, the Reservoir Characterization Project at the Colorado School of Mines and Occidental Petroleum Corporation (Oxy) (formerly the Anadarko Petroleum Corporation) collected time-lapse seismic data for characterization of changes in the reservoir caused by hydraulic fracturing and production in the Niobrara Formation and Codell Sandstone member of the Carlile Formation. We have acquired three multicomponent seismic surveys to understand the dynamic reservoir changes caused by hydraulic fracturing and production of 11 horizontal wells within a 1 mi2 section (the Wishbone Section). The time-lapse seismic survey acquisition occurred immediately after the wells were drilled, another survey after stimulation, and a third survey after two years of production. In addition, we integrate core, petrophysical properties, fault and fracture characteristics, as well as P-wave seismic data to illustrate reservoir properties prior to simulation and production. Core analysis indicates extensive amounts of bioturbation in zones of high total organic content (TOC). Petrophysical analysis of logs and core samples indicates that chalk intervals have high amounts of TOC (〉2%) and the lowest amount of clay in the reservoir interval. Core petrophysical characterization included X-ray diffraction analysis, mercury intrusion capillary pressure, N2 gas adsorption, and field emission scanning electron microscopy. Reservoir fractures follow four regional orientations, and chalk facies contain higher fracture density than marl facies. Integration of these data assist in enhanced well targeting and reservoir simulation.

Description: Four controlled-source audio-frequency magnetotelluric (CSAMT) profiles were conducted in early 2011 to investigate the distribution rules of oil shale in the Tongchuan area of the southern Ordos Basin. After terrain correction, we high-pass filtered the resulting images to prepare the 2D CSAMT apparent resistivity profiles for further analysis. Specifically, we have correlated the high-pass-filtered apparent resistivity anomalies to the distribution of oil shale seen in three available wells. The results reveal a close relationship between the oil shale distribution and the anomalously high-apparent-resistivity belt-like anomalies in the Tongchuan area. Our analysis indicates two prospective areas with different depths for oil shale within the study area.

Description: Gas-bearing prediction of tight sandstone reservoirs is significant but challenging due to the relationship between the gas-bearing property and its seismic response being nonlinear and complex. Although machine learning (ML) methods provide potential for solving the issue, the major challenge of ML applications to gas-bearing prediction is that of generating accurate and interpretable intelligent models with limited training sets. The k nearest neighbor ( kNN) method is a supervised ML method classifying an unlabeled sample according to its k neighboring labeled samples. We have introduced a kNN-based gas-bearing prediction method. The method can automatically extract a gas-sensitive attribute called the gas-indication local waveform similarity attribute (GLWSA) combining prestack seismic gathers with interpreted gas-bearing curves. GLWSA uses the local waveform similarity among the predicting samples and the gas-bearing training samples to indicate the existence of an exploitable gas reservoir. GLWSA has simple principles and an explicit geophysical meaning. We use a numerical model and field data to test the effectiveness of our method. The result demonstrates that GLWSA is good at characterizing the reservoir morphology and location qualitatively. When the method applies to the field data, we evaluate the performance with a blind well. The prediction result is consistent with the geologic law of the work area and indicates more details compared to the root-mean-square attribute.

Description: The Shaximiao Formation in the Zhongjiang Gas Field of the Sichuan Basin was initially a high-productivity gas field with the bright spot channel as the vital exploration target. With further development, gas wells were obtained in some nonbright spot areas, which caused interpreters to pay great attention to the channels with nonbright spot abnormal amplitudes. We have developed a method to delineate nonbright spot channels from the complicated sand-mudstone contact relationship. First, we classified sandstone into types I, IIa, IIb, and III, depending on the responses of the amplitude variation with offset from the drilled data, to produce a forward model. We the explain why the hidden channel cannot be identified using the full-angle stack seismic data based on this model. Afterward, we put forward a difference, between the synthetic seismogram responses of bright and nonbright channels, in creating seismic-to-well ties for nonbright channels. This difference from bright channels is that the synthetic data’s wave peak is not corresponding to the peak of the real seismic data. The wave trough has the same situation. Finally, we used far-angle stack seismic data to calculate coherent energy and instantaneous spectral attributes (the latter produced for red-green-blue blending) to identify the hidden channel. We observed that parts of the channel are more clearly visible in the far-angle stack than in the full-angle stack data. In the latter situation, we cannot describe the geometric shape of the channel elaborately. The Shaximiao Formation example is a relatively effective analog for nonbright spot plays compared with elsewhere.

Description: Energy exploration is becoming increasingly complex worldwide, and tight sandstone gas is an important field for the future development of the oil and gas industry. For the reservoir properties of the Shaximiao Gas Reservoir on the eastern slope of the Western Sichuan Depression in the Sichuan Basin, western China, it was found that the low-resistance characteristics of the reservoir complicate the gray characteristics among reservoir fluid property parameters. Some commonly used fluid property identification techniques, such as the flow zone index method, correlation analysis method of logging parameters, and traditional mathematical statistical methods, have poor fluid property evaluation results. Therefore, how to eliminate the influence of the gray features among the reservoir parameters on the identification of reservoir fluid properties and how to accurately identify the reservoir fluid properties are urgent problems that need to be solved. In this paper, we have developed a new method for identifying the fluid properties of tight sandstone reservoirs by combining gray system theory and multivariate statistical theory. This method can perform gray correlation weight analysis on parameters (combined parameters) closely related to fluid properties; furthermore, the logging identification method based on gray correlation weight analysis is used to identify reservoir fluid properties. The results indicate that the gray correlation weight analysis can accurately characterize the gray characteristics of reservoir fluid parameters and that the gray comprehensive correlation weight results are in good agreement with the production status of the studied gas reservoir. We used the method to identify the fluid properties of the target layer in 58 wells in the study area, and the discrimination rate of the model was 86.5%. In addition, the new model was used to predict the reservoir fluid properties of 12 newly drilled wells in the study area and the accuracy of the reservoir fluid property prediction was 91.67%.

Description: We analyzed a synthetic transfer zone and its associated fault planes and relay ramp in Penobscot, a potential offshore field in the Scotian Basin. Transfer zones are structural areas where one fault dies out and another fault begins, forming a relay ramp in the middle. They can be categorized as divergent, convergent, and synthetic transfer zones depending on the relative location and dipping directions of the faults. These zones not only play an important role in fluid migration but also help interpreters delineate secondary features such as fractures, splay shears, and Riedel faults. Commonly, those faults would branch into smaller splays and the relay ramp can get “breached” with connecting faults with the increase of slip. The study area in the Scotian Basin is characterized by two major listric normal faults dipping in the same direction giving rise to a synthetic transfer zone. These faults are clearly visible on seismic attributes, including curvature and coherence slices extracted along the top of the Cretaceous Petrel Formation. However, when analyzing the seismic attributes along the overlying Wyandot Formation’s top, we observe channel-like features, which run parallel as well as at an angle to these faults. However, when we performed further analysis using seismic amplitude’s vertical slices, interpreted horizons, and seismic attributes, we found that these features are not channels. We divided the features into two types: The first type is parallel to the main faults and can be associated with the grabens formed by synthetic and antithetic secondary faults (northeast–southwest); the second type is related to the polygonal faulting associated with differential compaction and gravitational loading of the Wyandot Chalk Formation. Apart from the two lineations, there are north-northeast–south-southwest-oriented lineations, which are an impression of basement faulting, and north-northwest–south-southeast-oriented lineations representing the acquisition footprint. Geological feature: Synthetic and antithetic faults in a synthetic transfer zone Seismic appearance: Channel-like feature on seismic attribute horizon slices Alternative interpretations: Channels Features with similar appearance: Steep slope channels Formation: Wyandot Age: Late Cretaceous Location: Penobscot, Nova Scotia, Canada Seismic data: Open-source SEG Analysis tools: Coherence and curvature attributes

Description: The phenomenon that frequency decreases and amplitude increases near the bottom of a gas layer on a seismic profile is called a low-frequency shadow, but this phenomenon may not occur in all gas reservoirs. When the tight gas reservoir is thick enough, spectral decomposition data after Fourier transformation will indicate characteristics similar to those of low-frequency shadows, which we call generalized low-frequency shadows. Compared to the dominant frequency of the nongas-bearing zone spectral, the dominant frequency of a gas zone moves toward the low end of the frequency range and the low-frequency amplitude increases accordingly. By analyzing known gas reservoirs such as the Sulige and Yanchang tight sandstones in the Ordos Basin and tight carbonate rocks in the Tarim Basin, we can see that, with the visual dominant seismic frequency close to 30 Hz, the peak frequency of the gas-bearing tight sandstones and tight dolomite reservoirs will move from approximately 30 to 10–15 Hz. There is a certain correlation among the drop of the dominant frequency of a tight gas reservoir, the attenuation energy difference, and the thickness and productivity of the gas layer. Several cases indicate that nearly all tight gas layers thicker than 15 m exhibit attenuation characteristics of generalized low-frequency shadows.

Description: Studies of the new growth and re-distribution of Cu-rich phases in chondrites of different petrologic subtypes can potentially provide insights into post-accretionary parent-body processes. We present a systematic study of the distribution of Cu-rich phases and metallic Cu in Ornans-like carbonaceous chondrites (CO3) that underwent little aqueous alteration or shock (most with shock stages of S1) but exhibit a range of thermal metamorphism (subtype 3.0–3.7). A comparison to ordinary chondrites (OCs), which have undergone a larger range of shock levels, allows us to constrain the relative roles of radiogenic and shock heating in the origin of Cu distribution in chondrites. We found that the Cu content of Ni-rich metal and calculated bulk Cu content of CO3 chondrites (based on mass-balance calculations) show an increase from CO3.0 to CO3.2 chondrites. We speculate that some unidentified phases in the matrix account for a significant portion (nearly ~100 ppm) of the Cu budget in bulk samples of CO3.0 chondrites, while Ni-rich metal is the main Cu-carrier for CO3.2–3.7 chondrites. Within CO3.2–3.7 chondrites, Cu and Ni contents of Ni-rich metal are positively correlated, showing a systematic decrease from lower to higher subtype (~0.41 wt% Cu and ~45.0 wt% Ni in CO3.2 Kainsaz; ~0.28 wt% Cu and ~38.8 wt% Ni in CO3.7 Isna). Metallic Cu grains were found in every sample of CO3.2–3.7 chondrites, but not in any CO3.0–3.1 chondrites. Metallic Cu is: (1) present at metallic-Fe-Ni-pyrrhotite interfaces; (2) associated with fine irregular pyrrhotite grains in Ni-rich-metal-pyrrhotite nodules; (3) associated with fizzed pyrrhotite (fine-grained mixtures of irregularly shaped metal grains surrounded by pyrrhotite); (4) present at the edges of metallic Fe-Ni grains; and (5) present as isolated grains. In some metallic-Cu-bearing mineral assemblages, pyrrhotite has higher Cu concentrations than adjacent Ni-rich metal and shows a drop in Cu concentration at the interface between metallic Cu and Cu-rich pyrrhotite. This implies that the precipitation of metallic Cu grains could be related to the local Cu enrichment of pyrrhotite. We consider that radiogenic heating is mainly responsible for the formation of opaque phases in CO chondrites based on the relatively slow metallographic cooling rate (~0.1–5 °C/Ma), the increasing uniformity of Ni contents in Ni-rich metal with increasing CO subtype (44.3 ± 17.3 wt% in CO3.00 to 38.8 ± 3.4 wt% in CO3.7 chondrite), and the relatively narrow range of pyrrhotite metal/sulfur ratios (~0.976–0.999). Metal/sulfur ratios of pyrrhotite grains in most CO3.2–3.7 chondrites (mean = ~0.986–0.997; except Lancé) are slightly higher than those in CO3.0–3.1 chondrites (mean = ~0.981–0.987; except Y-81020), possibly indicative of a release and re-mobilization of sulfur during progressive heating as previously reported for type-3 chondrites. In this regard, we suggest most metallic Cu grains in CO3 chondrites may have precipitated from Cu-rich pyrrhotite due to sulfidation of Fe-Ni metal during parent-body thermal metamorphism. Locally, a few metallic Cu grains associated with fizzed pyrrhotite could have formed during transient shock-heating. Both thermal and shock metamorphism could be responsible for the formation of metallic Cu. Although the systematic decrease in the Ni contents of Ni-rich metal from subtype-3.2 to subtype-3.8 also occurs in OCs, the average Cu contents of Ni-rich metal grains are indistinguishable among type-3 OCs of different subtypes. The paucity of metallic Cu in weakly shocked type-3 OCs could be related to: (1) the relatively low-bulk Cu contents of OCs, and/or (2) the relatively rapid metallographic cooling rates at

Description: Antigorite, a high-pressure polymorph of serpentine, is considered to be the most abundant hydrous mineral in subduction zones. Although antigorite dehydration is presumed as one of the origins of intermediate-depth earthquakes in the subduction zone, the amount of antigorite is uncertain because the amount of water infiltrated into the oceanic lithosphere is still debated. To investigate whether antigorite can be formed even with limited water availability, we conducted the axial deformation experiments of magnesium germanate at 1.2 GPa and T = 500–800 °C using a Griggs-type deformation apparatus. Magnesium germanate is an analog material of magnesium silicate, and the starting material was dried prior to experimentation. Nevertheless, the samples had initially high porosity, and hence a small amount of water (about 200 ppm wt H2O) was retained in the samples. In the samples deformed at 600 °C, stable slip occurred, and TEM analysis revealed that fine-grained platelets of germanate antigorite existed along the faults. A sharp absorption band assigned to the OH-stretching vibration of antigorite in Fourier transform infrared spectroscopic (FTIR) analysis also implies that antigorite was formed in the samples deformed at a temperature lower than 600 °C. Our results indicate that strain-induced hydration of germanate olivine results in antigorite formation even with only a small amount of water present. Thus, partial serpentinization in the oceanic lithosphere can occur under slight water infiltration due to the high strain accumulated by subduction.

Description: Gold (Au) deposits have formed in orogenic belts throughout Earth’s history. However, the upper temperature limits of orogenic Au vein formation are difficult to constrain because measurements made on fluid inclusions focus on intermediate to late-stage minerals (e.g., quartz and calcite) or are based on P-T estimates for the metamorphic mineral assemblages of the host rocks. We conducted a study of TiO2 polymorphs that are among the earliest minerals that grew in Au-bearing veins of the Dongyuan deposit, Jiangnan orogenic Au belt, South China. Based on Raman analyzes, we identified TiO2 polymorphs of anatase (with Raman peaks at 396, 515, and 638 cm−1), rutile (with Raman peaks at 235, 447, and 613 cm−1), and anatase–rutile intergrowths. Transmission electron microscope (TEM) confirmed the polymorphs identifying the [111] zone axis of anatase, [110] zone axis of rutile, and [111] and [111] zone axes of rutile–anatase intergrowths. The TiO2 polymorphs in the Dongyuan Au veins constrain a temperature range for early mineral precipitation in the veins of 450–550 °C. The results show that ore-forming fluids for this orogenic Au deposit emplaced in the shallow crust originated from deeper and hotter crustal levels (e.g., high-grade metamorphic rocks in the middle to lower crust).

Description: The transport of calcium carbonate (CaCO3) into the Earth’s interior through subduction is one of the key processes in the global cycling of carbon. To develop a better understanding of the CaCO3 structural stability during subduction processes, the phase transitions among CaCO3-I (calcite), CaCO3-II, -III/IIIb, and aragonite under pressure-temperature (P-T) conditions up to 2.5 GPa and 600 °C, in hydrous and anhydrous environments, were investigated using a hydrothermal diamond-anvil cell. One displacive and two reconstructive processes during the phase transitions among CaCO3 polymorphs were confirmed from the results obtained from in situ observations and Raman spectroscopic measurements. Meanwhile, the effect of Ca-substitutional metal cations (e.g., Mg2+) in CaCO3 and the presence of an aqueous fluid on the phase transition processes have been determined. Specifically, the CaCO3-I ↔ -II phase transition is a displacive process, occurring instantly at pressures varying from 1.6 GPa at room temperature to 1.5 GPa at 500 °C with the phase equilibrium boundary having a minimum P-T point at ~1.4 GPa at 300 °C, and is completely reversible upon cooling and decompression. The CaCO3-II → -III phase transition is a reconstructive process, observed at P-T conditions from 2.0 GPa at room temperature to 2.5 GPa at 150 °C, and is accomplished by solid recrystallization starting from CaCO3-II, transitioning through an intermediate CaCO3-IIIb, and ending at the CaCO3-III structure. The phase transition between CaCO3-I or -II and aragonite, which is also a reconstructive process, was found to occur by progressive solid recrystallization under high P-T hydrous and anhydrous conditions, or alternatively, via dissolution-precipitation under low-P-T hydrous conditions, depending on the presence of aqueous fluids and the heating rate of the system. The substitution for Ca2+ by other metal cations (e.g., Mg2+, Mn2+, Fe2+) in CaCO3 results in a significant increase in the pressures for the displacive and solid recrystallization reconstructive phase transitions, but has no detectable influence on the CaCO3-I/II ↔ aragonite transformation via a dissolution-precipitation process under hydrous conditions. Our results show that the presence of Ca-substitutional metal cations in CaCO3 is a key factor controlling the phase stability of CaCO3 under high P-T conditions, and suggest that aragonite should be the predominant phase in the upper mantle in subduction zones where the heating rate is very low and slab dehydration is prevalent.

Description: The Mogok metamorphic belt (MMB), Myanmar, is one of the most well-known gemological belts on Earth. Previously, 40Ar/39Ar, K-Ar, and U-Pb dating have yielded Jurassic-Miocene magmatic and metamorphic ages of the MMB and adjacent areas; however, no reported age data are closely related to the sapphire and moonstone deposits. Secondary ion mass spectrometry (SIMS) U-Pb dating of acicular rutile inclusions in sapphire and furnace step-heating 40Ar/39Ar dating of moonstone (antiperthite) in syenites from the MMB yield ages of 13.43 ± 0.92 and 13.55 ± 0.08 Ma, respectively, indicating both Myanmar sapphire and moonstone formed at the same time, and the ages are the youngest published in the region. The ages provide insight into the complex histories and processes of magmatism and metamorphism of the MMB, the formation of gemstone species in this belt, and the collision between India and Asia. In addition, our high field strength element data for the oriented rutile inclusions suggest an origin by co-precipitation rather than exsolution. In situ age determination of this nature is particularly significant since rutile inclusions in other gemstones, such as rubies, can be used to help constrain the geological history of their host rocks elsewhere.

Description: X-ray computed microtomography (CT) of impact rock varieties from the Kara astrobleme is used to test the method’s ability to identify the morphology and distribution of the rock components. Three types of suevitic breccias, clast-poor melt rock, and a melt clast from a suevite were studied with a spatial resolution of 24 µm to assess CT data values of 3D structure and components of the impactites. The purpose is first to reconstruct pore space, morphology, and distribution of all distinguishable crystallized melt, clastic components, and carbon products of impact metamorphism, including the impact glasses, after-coal diamonds, and other carbon phases. Second, the data are applied to analyze the morphology and distribution of aluminosilicate and sulfide components in the melt and suevitic breccias. The technical limitations of the CT measurements applied to the Kara impactites are discussed. Because of the similar chemical composition of the aluminosilicate matrix, glasses, and some lithic and crystal clasts, these components are hard to distinguish in tomograms. The carbonaceous matter has absorption characteristics close to air, so the pores and carbonaceous inclusions appear similar. However, X-ray microtomography could be used to prove the differences between the studied types of suevites from the Kara astrobleme using structural-textural features of the whole rock, porosity, and the distributions of carbonates and sulfides.

Description: Fibrous amphibole and clay mineral inclusions that form striking trapiche-like star patterns within quartz crystals from Inner Mongolia, China, present a challenge to uncover how these crystals grow and incorporate inclusions in a geological context. We propose that the patterns formed as a result of protogenic clay (ferrosaponite or nontronite) inclusions that were preferentially trapped on rough surfaces during quartz crystal growth. The rough surface texture of these crystals is the result of multiple growth centers during 2D nucleation and spread and split crystal formation. Observations via optical microscopy, cathodoluminescence, and three-dimensional micro-CT scanning highlight how the exterior surface textures on the termination of a complete quartz crystal mimic its interior inclusion patterns. Cathodoluminescence images, as well as varying aluminum concentrations along a core-to-exterior transect in a quartz crystal slice, suggest that the formation fluid underwent a heterogeneous chemical history. Measurements of Ti and observations of fluid inclusions suggest the quartz formed at a temperature of under 348 °C. This study presents the details surrounding split crystal growth in quartz in a natural geological setting, which has implications for inspiring new materials and may serve as an indicator for turbid and highly supersaturated formation fluid conditions in geological formations.

Description: The concentration of sulfur that can be dissolved in a silicate liquid is of fundamental importance because it is closely associated with several major Earth-related processes. Considerable effort has been made to understand the interplay between the effects of silicate melt composition and its capacity to retain sulfur, but the dependence on pressure and temperature is mostly based on experiments performed at pressures and temperatures below 6 GPa and 2073 K. Here we present a study of the effects of pressure and temperature on sulfur content at sulfide saturation of a peridotitic liquid. We performed 14 multi-anvil experiments using a peridotitic starting composition, and we produced 25 new measurements at conditions ranging from 7 to 23 GPa and 2173 to 2623 K. We analyzed the recovered samples using both electron microprobe and laser ablation ICP-MS. We compiled our data together with previously published data that were obtained at lower P-T conditions and with various silicate melt compositions. We present a new model based on this combined data set that encompasses the entire range of upper mantle pressure-temperature conditions, along with the effect of a wide range of silicate melt compositions. Our findings are consistent with earlier work based on extrapolation from lower-pressure and lower-temperature experiments and show a decrease of sulfur content at sulfide saturation (SCSS) with increasing pressure and an increase of SCSS with increasing temperature. We have extrapolated our results to pressure-temperature conditions of the Earth’s primitive magma ocean, and show that FeS will exsolve from the molten silicate and can effectively be extracted to the core by a process that has been termed the “Hadean Matte.” We also discuss briefly the implications of our results for the lunar magma ocean.

Description: A new mineral of the beryl group, johnkoivulaite, Cs(Be2B)Mg2Si6O18, was recovered from the gem gravels in the Pein Pyit area of the Mogok region in Myanmar. Thus far, only a single crystal has been identified. It has dimensions of about 5.8 × 5.7 × 5.5 mm. This specimen has an irregular shape but still has discernible crystal form with geometric growth patterns observed on the crystal faces. The crystal of johnkoivulaite is grayish-violet in color and strongly pleochroic, going from nearly colorless with E┴c to dark bluish-violet with E||c. Johnkoivulaite has a Mohs hardness of about 7½ and a measured density of 3.01(10) g/cm3. It is uniaxial (–) with ω = 1.607(1) and ε = 1.605(1) (white light). Electron microprobe analyses gave the empirical formula of (Cs0.85K0.10Na0.01)(Be1.88B1.12)(Mg1.66Fe0.27Mn0.01Al0.05) (Si5.98)O18 with Be calculated by stoichiometry and confirmed by LA-ICP-MS measurements. Johnkoivulaite is hexagonal, P6/mmc (no. 192) with a = 9.469(2), c = 9.033(2) Å, V = 701.5(3) Å3, and Z = 2. Johnkoivulaite is isostructural with beryl and exhibits partial substitution of B for Be at the distorted tetrahedral site, Mg for Al at the octahedral site, and Cs in the channel sites within the stacked Si6O18 rings. This substitution can be written as (CsMg2B)(☐Al2Be)–1. Johnkoivulaite, the seventh member of the beryl group, is named in honor of gemologist John Koivula in recognition of his contributions to mineralogy and gemology.

Description: In this paper, we report a first-principles Molecular Dynamics (FPMD) study of interfacial structures and acidity constants of goethite. The pKa values of the groups on (010), (110), and (021) surfaces (space group Pbnm) are derived with the FPMD based vertical energy gap technique. The results indicate that major reactive groups include ≡Fe2OH2 and ≡FeOH2 on (010), ≡FeOH2, ≡Fe3OLH, and ≡Fe3OUH on (110), and ≡FeOhH2 and ≡Fe2OH on (021). The interfacial structures were characterized in detail with a focus on the hydrogen bonding environment. With the calculated pKa values, the point of zero charges (PZCs) of the three surfaces are derived and the overall PZC range of goethite is found to be consistent with the experiment. We further discuss the potential applications of these results in future studies toward understanding the environmental processes of goethite.

Description: The mobility of Ti, a member of high field strength elements, in metamorphic fluids is crucial to understand the recycling of commonly perceived nominally soluble elements and for mass-flux calculations during crustal processes. In this study, we present evidence for large-scale Ti mobility from a suite of clinohumite±spinel-bearing dolomitic marbles from the Makrohar area in central India. The studied rocks mostly contain dolomite and calcite (in subequal proportions) and a subordinate amount of forsterite. It commonly develops 1–5 cm thick, laterally continuous, mostly parallel, sometimes anastomosing, brown-colored clinohumite rich bands with variable spinel. Clinohumite has moderate Ti and F (TiO2 = 0.55–2.88 wt%; F = 0.94–1.88 wt%; n = 32). Textural and phase equilibria modeling indicate that clinohumite grew at the expense of forsterite + dolomite under static conditions due to infiltration of F- and Ti-bearing extremely H2O-rich fluids (XCO2 〈 0.03), at ~5–6 kbar pressure and ~650–700 °C temperature. The Ti and F were most likely supplied by highly channelized aqueous fluids restricted within the centimeter-thick bands. The negative volume change of the reactions further facilitated fluid ingress. The lateral continuity of the bands over several meters across multiple out-crops indicates that Ti was mobile at the meter to kilometer scale. The results are in accordance with experimental studies that solubility of Ti increases in the presence of halides and imply that Ti may be much more mobile in metamorphic fluids during regional metamorphism than previously anticipated.

Description: δ-AlOOH has emerged as a promising candidate for water storage in the lower mantle and could have delivered water into the bottom of the mantle. To date, it still remains unclear how the presence of iron affects its elastic, rheological, vibrational, and transport properties, especially across the spin crossover. Here, we conducted high-pressure X-ray emission spectroscopy experiments on a δ-(Al0.85Fe0.15) OOH sample up to 53 GPa using silicone oil as the pressure transmitting medium in a diamond-anvil cell. We also carried out laser Raman measurements on δ-(Al0.85Fe0.15)OOH and δ-(Al0.52Fe0.48)OOH up to 57 and 62 GPa, respectively, using neon as the pressure-transmitting medium. Evolution of Raman spectra of δ-(Al0.85Fe0.15)OOH with pressure shows two new bands at 226 and 632 cm−1 at 6.0 GPa, in agreement with the transition from an ordered (P21nm) to a disordered hydrogen bonding structure (Pnnm) for δ-AlOOH. Similarly, the two new Raman bands at 155 and 539 cm−1 appear in δ-(Al0.52Fe0.48)OOH between 8.5 and 15.8 GPa, indicating that the incorporation of 48 mol% FeOOH could postpone the order-disorder transition upon compression. On the other hand, the satellite peak (Kβ′) intensity of δ-(Al0.85Fe0.15)OOH starts to decrease at ~30 GPa and it disappears completely at 42 GPa. That is, δ-(Al0.85Fe0.15)OOH undergoes a gradual electronic spin-pairing transition at 30–42 GPa. Furthermore, the pressure dependence of Raman shifts of δ-(Al0.85Fe0.15)OOH discontinuously decreases at 32–37 GPa, suggesting that the improved hydrostaticity by the use of neon pressure medium could lead to a relatively narrow spin crossover. Notably, the pressure dependence of Raman shifts and optical color of δ-(Al0.52Fe0.48)OOH dramatically change at 41–45 GPa, suggesting that it probably undergoes a relatively sharp spin transition in the neon pressure medium. Together with literature data on the solid solutions between δ-AlOOH and ε-FeOOH, we found that the onset pressure of the spin transition in δ-(Al,Fe)OOH increases with increasing FeOOH content. These results shed new insights into the effects of iron on the structural evolution and vibrational properties of δ-AlOOH. The presence of FeOOH in δ-AlOOH can substantially influence its high-pressure behavior and stability at the deep mantle conditions and play an important role in the deep-water cycle.

Description: While much progress has been made in electron-probe microanalysis (EPMA) to improve the accuracy of point analysis, the same level of attention has not always been applied to the quantification of wavelength-dispersive spectrometry (WDS) X-ray intensity maps at the individual pixel level. We demonstrate that the same level of rigor applied in traditional point analysis can also be applied to the quantification of pixels in X-ray intensity maps, along with additional acquisition and quantitative processing procedures to further improve accuracy, precision, and mapping throughput. Accordingly, X-ray map quantification should include pixel-level corrections for WDS detector deadtime, corrections for changes in beam current (beam drift), changes in standard intensities (standard drift), high-accuracy removal of background intensities, quantitative matrix corrections, quantitative correction of spectral interferences, and, if required, time-dependent corrections (for beam and/or contamination sensitive materials). The purpose of quantification at the pixel level is to eliminate misinterpretation of intensity artifacts, inherent in raw X-ray intensity signals, that distort the apparent abundance of an element. Major and minor element X-ray signals can contain significant artifacts due to absorption and fluorescence effects. Trace element X-ray signals can contain significant artifacts where phases with different average atomic numbers produce different X-ray continuum (bremsstrahlung) intensities, or where a spectral interference, even an apparently minor one, can produce a false-positive intensity signal. The methods we propose for rigorous pixel quantification require calibration of X-ray intensities on the instrument using standard reference materials, as we already do for point analysis that is then used to quantify multiple X-ray maps, and thus the relative time overhead associated with such pixel-by-pixel quantification is small. Moreover, the absolute time overhead associated with this method is usually less than that required for quantification using manual calibration curve methods while resulting in significantly better accuracy. Applications to geological, synthetic, or engineering materials are numerous as quantitative maps not only show compositional 2D variation of fine-grained or finely zoned structures but also provide very accurate quantitative analysis, with precision approaching that of a single point analysis, when multiple-pixel averaging in compositionally homogeneous domains is utilized.

Description: Carletonmooreite (IMA 2018-68), Ni3Si, is a new nickel silicide mineral that occurs in metal nodules from the Norton County aubrite meteorite. These nodules are dominated by low-Ni iron (kamacite), with accessory schreibersite, nickelphosphide, perryite, and minor daubréelite, tetrataenite, taenite, and graphite. The chemical composition of the holotype carletonmooreite determined by wavelength-dispersive electron-microprobe analysis is (wt%) Ni 82.8 ± 0.4, Fe 4.92 ± 0.09, and Si 13.08 ± 0.08 (n = 6, total = 100.81) giving an empirical formula of (Ni2.87Fe0.18)Σ3.05Si0.95, with an end-member formula of Ni3Si. Further grains discovered in the specimen after the new mineral submission extend the composition, i.e., (wt%) Ni 81.44 ± 0.82, Fe 5.92 ± 0.93, Cu 0.13 ± 0.02, and Si 13.01 ± 0.1 (n = 11, total = 100.51 ± 0.41), giving an empirical formula (Ni2.83Fe0.22Cu0.004)Σ3.05Si0.95. The backscat tered electron-diffraction patterns were indexed by the Pm3m auricupride (AuCu3)-type structure and give a best fit to synthetic Ni3Si, with a = 3.51(1) Å, V = 43.2(4) Å3, Z = 1, and calculated density of 7.89 g/cm3. Carletonmooreite is silver colored with an orange tinge, isotropic, with a metallic luster and occurs as euhedral to subhedral crystals 1 × 5 µm to 5 × 14 µm growing on tetrataenite into kamacite. The dominant silicide in the Norton County aubrite metal nodules is perryite (Ni,Fe)8(Si,P)3, with carletonmooreite restricted to localized growth on rare plessite fields. The isolated nature of small euhedral carletonmooreite single crystals suggests low-temperature growth via solid-state diffusion from the surrounding kamacite and epitaxial growth on the tetrataenite. This new mineral is named in honor of Carleton B. Moore, chemist and geologist, and founding director of the Center for Meteorite Studies at Arizona State University, for his many contributions to cosmochemistry and meteoritics.

Description: Spectral features of hydrogen defects in natural mantle minerals derive from physico-chemical conditions of the lithosphere. Although hydrogen defects in synthetic orthopyroxene have been well investigated, their complex spectral features in natural orthopyroxenes are still difficult to decipher. To clarify this issue, it is indispensable to reveal what happens to hydrogen defects during high-temperature processes, thereby fingerprinting the origins of hydrogen defects observed in natural orthopyroxene. Here, we carry out Fourier transform infrared spectroscopic studies on hydrogen defects of three natural orthopyroxenes at elevated temperatures to 1000 °C. Hydrogen defects display reversible disordering at temperatures above 700 °C, which is different from those at ambient conditions. Moreover, hydrogen diffusivities are significantly different between the orthopyroxene samples from different tectonic settings despite their similar iron contents. Even for the same crystal, different hydrogen defects display different diffusion behaviors. Hydrogen defects corresponding to the 3420 cm−1 band have the fastest diffusivity relative to the other hydrogen defects. Most importantly, hydrogen defects can redistribute in the crystal, with new hydrogen defects produced at the cost of the initial hydrogen defects rather than involving a reaction with an external hydrogen source. Combining these findings with previously reported hydrogen defects in natural olivine and clinopyroxene at high temperatures, we propose that: (1) to correctly relate hydrogen defects features to geological processes, it is imperative to understand their behavior and origin, and (2) hydrogen disordering should be taken into account when predicting and extrapolating data on physical properties of the mantle from room-temperature measurements.

Description: Sound velocities of iron and iron-based alloys at high pressure and high temperature are crucial for understanding the composition and structure of Earth’s and other telluric planetary cores. In this study, we performed ultrasonic interferometric measurements of both compressional (νP) and shear (νS) velocities on a polycrystalline body-centered-cubic (bcc)-Fe90Ni10 up to 8 GPa and 773 K. The elastic moduli and their pressure and temperature derivatives are derived from least-square fits to third-order finite strain equations, yielding KS0 = 154.2(8) GPa, G0 = 73.2(2) GPa, KS0′ = 4.6(2), G0′ = 1.5(1), ∂KS/∂T = –0.028(1) GPa/K, and ∂G/∂T = –0.023(1) GPa/K. A comparison with literature data on bcc-Fe suggests that nickel not only decreases both P and S wave velocities but also weakens the temperature effects on the elastic moduli of Fe-Ni alloys.