Man Hyok Song, Jin Gyong Ho, Chol Kim, Yong O. Chol, Song Lyu
{"title":"A Method to Extract Image Features and Lineaments Based on a Multi-hillshade Continuous Wavelet Transform","authors":"Man Hyok Song, Jin Gyong Ho, Chol Kim, Yong O. Chol, Song Lyu","doi":"10.1007/s11004-024-10146-5","DOIUrl":"https://doi.org/10.1007/s11004-024-10146-5","url":null,"abstract":"<p>This paper presents a new method for extracting the image features and lineaments related to local extrema of an image or a digital elevation model (DEM) such as ridges and valleys based on the continuous wavelet transform (CWT) of a set of variously illuminated hillshades. The method originates from the principle that a hillshade can exactly reflect the lineaments nearly perpendicular to the illumination direction of the hillshade, but not other ones. The method consists of four steps: (1) preparation of a set of differently illuminated hillshades of the input data, (2) detection of directional edges nearly perpendicular to the illumination direction from each hillshade based on the CWT, (3) a combination of multidirectional edges into an omnidirectional feature image, and (4) identification of lineaments through linkage and linearization of image feature lines. CWT coefficients of each hillshade are used to calculate the gradient and its direction of the hillshade. For each hillshade, directional edge pixels where the gradient direction is parallel to the illumination direction are selectively detected to form accurate and solitary image feature lines related to local extrema of the input data. Directional edges of each hillshade are easily classified into positive and negative edges using the hillshade gradient. As they have similar directions, they are easily linked to form longer line raster objects, which are converted into vector objects, that is, directional lineaments. The multidirectional edges and lineaments given from all the hillshades are combined to form an omnidirectional feature image and a group of omnidirectional lineaments. Its application to real DEMs shows the demonstrated advantages of the proposed method over other methods and the similarity between detected lineaments and fault lines in the study area.</p>","PeriodicalId":51117,"journal":{"name":"Mathematical Geosciences","volume":"43 1","pages":""},"PeriodicalIF":2.6,"publicationDate":"2024-06-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141256075","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Oscar Ovanger, Jo Eidsvik, Jacob Skauvold, Ragnar Hauge, Ingrid Aarnes
{"title":"Addressing Configuration Uncertainty in Well Conditioning for a Rule-Based Model","authors":"Oscar Ovanger, Jo Eidsvik, Jacob Skauvold, Ragnar Hauge, Ingrid Aarnes","doi":"10.1007/s11004-024-10144-7","DOIUrl":"https://doi.org/10.1007/s11004-024-10144-7","url":null,"abstract":"<p>Rule-based reservoir models incorporate rules that mimic actual sediment deposition processes for accurate representation of geological patterns of sediment accumulation. Bayesian methods combine rule-based reservoir modelling and well data, with geometry and placement rules as part of the prior and well data accounted for by the likelihood. The focus here is on a shallow marine shoreface geometry of ordered sedimentary packages called bedsets. Shoreline advance and sediment build-up are described through progradation and aggradation parameters linked to individual bedset objects. Conditioning on data from non-vertical wells is studied. The emphasis is on the role of ‘configurations’—the order and arrangement of bedsets as observed within well intersections in establishing the coupling between well observations and modelled objects. A conditioning algorithm is presented that explicitly integrates uncertainty about configurations for observed intersections between the well and the bedset surfaces. As data volumes increase and model complexity grows, the proposed conditioning method eventually becomes computationally infeasible. It has significant potential, however, to support the development of more complex models and conditioning methods by serving as a reference for consistency in conditioning.</p>","PeriodicalId":51117,"journal":{"name":"Mathematical Geosciences","volume":"51 1","pages":""},"PeriodicalIF":2.6,"publicationDate":"2024-06-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141195979","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Generalized Solution for Double-Porosity Flow Through a Graded Excavation Damaged Zone","authors":"Kristopher L. Kuhlman","doi":"10.1007/s11004-024-10143-8","DOIUrl":"https://doi.org/10.1007/s11004-024-10143-8","url":null,"abstract":"<p>Prediction of flow to boreholes or excavations in fractured low-permeability rocks is important for resource extraction and disposal or sequestration activities. Analytical solutions for fluid pressure and flowrate, when available, are powerful, insightful, and efficient tools enabling parameter estimation and uncertainty quantification. A flexible porous media flow solution for arbitrary physical dimensions is derived and extended to double porosity for converging radial flow when permeability and porosity decrease radially as a power law away from a borehole or opening. This distribution can arise from damage accumulation due to stress relief associated with drilling or mining. The single-porosity graded conductivity solution was initially found for heat conduction, the arbitrary dimension flow solution comes from hydrology, and the solution with both arbitrary dimension and graded permeability distribution appeared in reservoir engineering. These existing solutions are combined and extended here to two implementations of the double-porosity conceptual model, for both a simpler thin-film mass transfer and more physically realistic diffusion between fracture and matrix. This work presents a new specified-flowrate solution with wellbore storage for the simpler double-porosity model, and a new, more physically realistic solution for any wellbore boundary condition. A new closed-form expression is derived for the matrix diffusion solution (applicable to both homogeneous and graded problems), improving on previous infinite series expressions.</p>","PeriodicalId":51117,"journal":{"name":"Mathematical Geosciences","volume":"12 1","pages":""},"PeriodicalIF":2.6,"publicationDate":"2024-05-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140927233","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Ivana Pavlů, Jitka Machalová, Raimon Tolosana-Delgado, Karel Hron, Kai Bachmann, Karl Gerald van den Boogaart
{"title":"Principal Component Analysis for Distributions Observed by Samples in Bayes Spaces","authors":"Ivana Pavlů, Jitka Machalová, Raimon Tolosana-Delgado, Karel Hron, Kai Bachmann, Karl Gerald van den Boogaart","doi":"10.1007/s11004-024-10142-9","DOIUrl":"https://doi.org/10.1007/s11004-024-10142-9","url":null,"abstract":"<p>Distributional data have recently become increasingly important for understanding processes in the geosciences, thanks to the establishment of cost-efficient analytical instruments capable of measuring properties over large numbers of particles, grains or crystals in a sample. Functional data analysis allows the direct application of multivariate methods, such as principal component analysis, to such distributions. However, these are often observed in the form of samples, and thus incur a sampling error. This additional sampling error changes the properties of the multivariate variance and thus the number of relevant principal components and their direction. The result of the principal component analysis becomes an artifact of the sampling error and can negatively affect the subsequent data analysis. This work presents a way of estimating this sampling error and how to confront it in the context of principal component analysis, where the principal components are obtained as a linear combination of elements of a newly constructed orthogonal spline basis. The effect of the sampling error and the effectiveness of the correction is demonstrated with a series of simulations. It is shown how the interpretability and reproducibility of the principal components improve and become independent of the selection of the basis. The proposed method is then applied on a dataset of grain size distributions in a geometallurgical dataset from Thaba mine in the Bushveld complex.</p>","PeriodicalId":51117,"journal":{"name":"Mathematical Geosciences","volume":"42 1","pages":""},"PeriodicalIF":2.6,"publicationDate":"2024-05-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140883589","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Geologically Constrained Convolutional Neural Network for Mineral Prospectivity Mapping","authors":"Fanfan Yang, Renguang Zuo","doi":"10.1007/s11004-024-10141-w","DOIUrl":"https://doi.org/10.1007/s11004-024-10141-w","url":null,"abstract":"<p>Various deep learning algorithms (DLAs) have been successfully employed for mineral prospectivity mapping (MPM) to support mineral exploration, due to their superior nonlinear extraction capabilities. DLAs algorithms are typically purely data-driven approaches that may ignore the geological domain knowledge. This renders the predictive results inconsistent with the mineralization mechanism and results in poor interpretation. In this study, a geologically constrained convolutional neural network (CNN) that involves soft and hard geological constraints was proposed for mapping gold polymetallic mineralization potential in western Henan Province of China. A penalty term based on the controlling equation of the spatial coupling relationship between the ore-controlling strata and gold deposits was constructed as a soft constraint to guide the CNN model training according to additional prior geological knowledge. In addition, domain knowledge related to mineralization processes and a geochemical indicator were simultaneously embedded as hard constraints in the feature extractor and classifier of the CNN, respectively, to control the model training based on the mineralization mechanism. The comparative experiments demonstrated that the geologically constrained CNN was superior to other models, thus indicating that the coupling of data and domain knowledge is effective for MPM and further improves the rationality and interpretability of the obtained results.</p>","PeriodicalId":51117,"journal":{"name":"Mathematical Geosciences","volume":"22 1","pages":""},"PeriodicalIF":2.6,"publicationDate":"2024-04-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140812500","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"And the 2024 Krumbein Medalist of the IAMG is…","authors":"Eric Grunsky","doi":"10.1007/s11004-024-10140-x","DOIUrl":"https://doi.org/10.1007/s11004-024-10140-x","url":null,"abstract":"","PeriodicalId":51117,"journal":{"name":"Mathematical Geosciences","volume":"6 1","pages":""},"PeriodicalIF":2.6,"publicationDate":"2024-04-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140588575","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Maxime Keutgen De Greef, Gert Jan Weltje, Irène Gijbels
{"title":"Estimating Rock Composition from Replicate Geochemical Analyses: Theory and Application to Magmatic Rocks of the GeoPT Database","authors":"Maxime Keutgen De Greef, Gert Jan Weltje, Irène Gijbels","doi":"10.1007/s11004-024-10138-5","DOIUrl":"https://doi.org/10.1007/s11004-024-10138-5","url":null,"abstract":"<p>Chemical analyses of powdered rocks by different laboratories often yield varying results, requiring estimation of the rock’s true composition and associated uncertainty. Challenges arise from the peculiar nature of geochemical data. Traditionally, major and trace elements have been measured using different methods, resulting in chemical analyses where the sum of the parts fluctuates around 1 rather than precisely totaling 1. Additionally, all chemical analyses contain an undisclosed mass fraction representing undetected chemical elements. Because of this undisclosed and unknown mass fraction, geochemical data represent a particular kind of compositional data in which closure to unity is not guaranteed. We argue that chemical analyses exist in the hypercube while being sampled from a true composition residing in the simplex. Therefore, we propose an algorithm that generates random chemical analyses by simulating the data acquisition protocol in geochemistry. Using the algorithm’s output, we measure the bias and mean squared error (MSE) of various estimators of the true mean composition. Additionally, we explore the impact of missing values on estimator performance. Our findings reveal that the optimized binary log-ratio mean, a new estimator, exhibits the lowest MSE and bias. It performs well even with up to 70% missing values, in contrast to other classical estimators such as the arithmetic mean or the geometric mean. Applying our approach to the GeoPT database, which contains replicate analyses of igneous rocks from numerous geochemical laboratories, we introduce an outlier detection technique based on the Mahalanobis distance between a laboratory’s logit coordinates and the optimized mean estimate. This enables a probabilistic ranking of laboratories based on the atypicality of their performance. Finally, we offer an accessible R implementation of our findings through the GitHub repository linked to this paper [subject classification numbers: 10 (compositions) 85 (statistics)].</p>","PeriodicalId":51117,"journal":{"name":"Mathematical Geosciences","volume":"56 1","pages":""},"PeriodicalIF":2.6,"publicationDate":"2024-04-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140588583","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"2.5D Hexahedral Meshing for Reservoir Simulations","authors":"","doi":"10.1007/s11004-023-10106-5","DOIUrl":"https://doi.org/10.1007/s11004-023-10106-5","url":null,"abstract":"<h3>Abstract</h3> <p>We present a new method for generating pure hexahedral meshes for reservoir simulations. The grid is obtained by extruding a quadrangular mesh, using ideas from the latest advances in computational geometry, specifically the generation of semi-structured quadrangular meshes based on global parameterization. Hexahedral elements are automatically constructed to smoothly honor the geometry of input features (domain boundaries, faults, and horizons), thus making it possible to be used for multiple types of physical simulations on the same mesh. The main contributions are as follows: the introduction of a new semi-structured hexahedral meshing workflow producing high-quality meshes for a wide range of fault systems, and the study and definition of weak verticality on triangulated surface meshes. This allows us to design better and more robust algorithms during the extrusion phase along non-vertical faults. We demonstrate (i) the simplicity of using such hexahedral meshes generated using the proposed method for coupled flow-geomechanics simulations with state-of-the-art simulators for reservoir studies, and (ii) the possibility of using such semi-structured hexahedral meshes in commercial structured flow simulators, offering an alternative gridding approach to handle a wider family of fault networks without recourse to the stair-step fault approximation. </p>","PeriodicalId":51117,"journal":{"name":"Mathematical Geosciences","volume":"56 1","pages":""},"PeriodicalIF":2.6,"publicationDate":"2024-04-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140588388","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Effects of Fluid Pressure Development on Hydrothermal Mineralization via Cellular Automaton Simulation","authors":"Yihui Xiong, Renguang Zuo, Oliver P. Kreuzer","doi":"10.1007/s11004-024-10139-4","DOIUrl":"https://doi.org/10.1007/s11004-024-10139-4","url":null,"abstract":"<p>The behavior and evolution trajectory of hydrofracture, which show a close relationship with the hydrothermal mineralization process, is greatly influenced by fluid flow and fluid pressure. However, further investigation is needed to achieve an in-depth understanding of the formation and evolution mechanisms behind the link between the rate of fluid pressure development and the occurrence of induced hydrofracture and mineralization process. We considered different fluid pressure development rates as the initial data for a cellular automaton model. With the increase in the fluid pressure increase rates, the corresponding hydrofracture became more focused, changing in scale from a large number of small-scale hydrofractures to a small number of large-scale hydrofractures. Episodes of fluid pressure fluctuation induced by either low or high fluid pressure increase rates were shown to trigger mineral precipitation and further contribute to the generation of strong spatially structured and enriched geochemical patterns. Moreover, the correlation length at the percolation threshold, which is of great significance to the degree and scale of mineralization, increased with the increasing fluid pressure increase rates. It was concluded that computational grids with high fluid pressure increase rates are much more prone to produce enriched geochemical patterns with strong spatial structures than grids with low fluid pressure increase rates owing to a larger correlation length at the percolation threshold. These results suggest that the way of fluid pressure development is a key factor for quantifying the behavior of hydrofracture and mineralization process.</p>","PeriodicalId":51117,"journal":{"name":"Mathematical Geosciences","volume":"61 1","pages":""},"PeriodicalIF":2.6,"publicationDate":"2024-04-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140588511","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Simulation Enhancement GAN for Efficient Reservoir Simulation at Fine Scales","authors":"","doi":"10.1007/s11004-024-10136-7","DOIUrl":"https://doi.org/10.1007/s11004-024-10136-7","url":null,"abstract":"<h3>Abstract</h3> <p>In this paper, an innovative approach for enhancing fluid transport modeling in porous media is presented, which finds application in various fields, including subsurface reservoir modeling. Fluid flow models are typically solved numerically by addressing a system of partial differential equations (PDEs) using methods such as finite difference and finite volume. However, these processes can be computationally demanding, particularly when aiming for high precision on a fine scale. Researchers have increasingly turned to machine learning to explore solutions for PDEs in order to improve simulation efficiency. The proposed method combines an adaptive multi-scale strategy with generative adversarial networks (GAN) to increase simulation efficiency on a fine scale. The devised model, called simulation enhancement GAN (SE-GAN), takes coarse-scale simulation results as input and generates fine-scale results in conjunction with the provided petrophysical properties. With this new approach, a deep learning model is trained to map coarse-scale results to fine-scale outcomes, rather than directly solving the fluid flow model. Case studies reveal that SE-GAN can achieve a significant improvement in accuracy while reducing computational time compared to the original fine-scale simulation solver. A comprehensive evaluation of numerical experiments is conducted to elucidate the benefits and limitations of this method. The potential of SE-GAN in accelerating the numerical solver for reservoir simulations is also demonstrated.</p>","PeriodicalId":51117,"journal":{"name":"Mathematical Geosciences","volume":"11 1","pages":""},"PeriodicalIF":2.6,"publicationDate":"2024-03-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140298068","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}