Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “quantification”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,171 records · Page 65
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Defect and satellite characteristics of additive manufacturing metal powders

Metal additive manufacturing (AM) requires high-quality metal powders to three-dimensionally (3D) print metallic components with complex and customizable geometries. The lack of quantification of AM metal powders creates quality control challenges for 3D printed components, increases the uncertainty of printing reliability and net cost of inspected and certified printed components, and reduces the recyclability of used powders. However, critical characteristics of AM metal powders that are decisive factors for the 3D printing process, such as internal porosity, contamination, and satellite feature, remain ambiguous. In this work, we developed a novel approach to 3D quantify key characteristics of AM metal powders down to individual particles by using high-resolution synchrotron x-ray computed tomography. Empowered by the penetrative capability of high-energy x-ray, internal porosity and contamination within as-atomized metal powders from high-entropy alloys to nickel-based superalloys were evaluated. Additionally, the newly-developed dispersion method enables the homogeneous separation of individual particles, and consequently, results in the implementation of 3D particle shape analysis. To resolve a major challenge of identification and quantification of satellite-feature particles in as-atomized AM metal powders, the satellite features were quantitated by modeling and analyzing the shape parameter of local thickness variance. Furthermore, the 3D analytical methods of particle assessment in this study can be applied to other materials systems like rock, food, and pharmaceutical particles, and provide insights for process optimization across powder metallurgy, concrete, food and pharmaceutical manufacturing, and AM industries.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Electrochemical Characterization of Photo-Driven Hole-Scavenging by Cadmium Sulfide Quantum Dot-Nitrogenase Biohybrid Complexes

Light-driven biohybrid systems that couple semiconductor nanocrystals with enzymes offer a promising strategy for solar-to-chemical energy-conserving reduction reactions, yet are often limited by inefficient hole scavenging. Hole scavenging is critical for maintaining charge separation in the light-absorbing electron donor molecule and, thus, for sustaining catalytic turnover, making quantification essential to improving system efficiency. Here, we introduce a photoelectrochemical approach for real-time monitoring of hole scavenging in cadmium sulfide (CdS) quantum dot (QD)-nitrogenase MoFe protein biohybrids that catalyze the reduction of H+ and N2. Using hydroquinone (HQ) as the sacrificial electron donor (SED), oxidation of HQ to benzoquinone (BQ) by photogenerated valence-band holes is coupled to the electrochemical reduction of BQ at an electrode, enabling quantification of the hole-scavenging reaction via chronoamperometry. This approach provides a continuous, real-time readout of charge-transfer dynamics under illumination. Systematic variation of light intensity and donor concentration reveals that hole scavenging scales with photon flux and exhibits a non-linear dependence on SED concentration, with an optimum at intermediate HQ concentrations. These results establish chronoamperometry as a powerful tool for probing photoinduced charge transfer, advancing the quantitative understanding of hole-scavenging in biohybrid systems, and offering a generalizable framework for optimizing solar-driven biocatalysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Dimensionality effect of conductive carbon fillers in LiNi 1/3 Mn 1/3 Co 1/3 O 2 cathode

Developing advanced electrode architectures through modifying active materials, conductive fillers, binders, and electrolytes as well as processing methods has drawn significant research interest. Due to the insufficient electrical conductivity of many active materials, adding conductive carbon fillers to composite electrodes provides the necessary electrical conductivity. The dimensionality effect among different conductive fillers has a significant impact on electrochemistry, which can be associated with morphological and chemical heterogeneities of electrodes. In this study, synchrotron X-ray mosaic nanotomography and X-ray spectroscopy nanoimaging provided direct three-dimensional (3D) visualization and quantification capabilities to investigate the dimensionality effects of Super P (SP) and single-walled carbon nanotube (SWCNT) fillers on the capacity retention of LiNi 1/3 Mn 1/3 Co 1/3 O 2 (NMC111). The results indicate that NMC/SWCNT electrodes, with a wrapping effect from the SWCNTs, exhibited more homogeneous particle size distributions, morphological changes, and chemical states than NMC/SP electrodes, without the wrapping effect. This work developed a framework of 3D quantification methods to study the capacity fading behavior associated with morphological and chemical heterogeneities and paved the way toward designing electrodes for high rate energy storage applications.

25 ENERGY STORAGE↗
DOE OSTI2021

A Transient-Response methodology based on experiments and modeling for Cu-Redox Half-Cycle kinetic analysis on a Cu-SSZ-13 SCR catalyst

In this report the Standard SCR reaction catalyzed by Cu-SSZ-13 is a redox process consisting of a reduction half cycle (RHC) and an oxidation half cycle (OHC) that cycle the active Cu sites between the Cu(II) and Cu(I) states. In the current work, a transient-response methodology consisting of experimental transient response Cu-redox (TRCR) measurements and kinetic modeling was developed for detailed study of individual SCR Cu-redox half cycles. The TRCR protocol allows quantification of the reducible Cu density, Cu(II)-Cu(I) partitioning, relative native RHC and OHC rates, and limiting half cycle during SCR. The half-cycle kinetics are studied over a wide (170–450 °C) temperature, and in differential segments along the catalyst length using spatially resolved capillary inlet mass-spectrometry (SpaciMS). The protocol alone provides two independent measures indicating that OHC increases faster than RHC with temperature, and that at all but the lowest temperature SCR is RHC limited. Introducing the transient kinetic models allows half-cycle reaction pathways, orders and activation energies to be studied and determined; e.g., RHC involves surface and gas-phase NH 3 routes, and OHC with O 2 involves dimer formation. Application of the full methodology provides further confirmation and quantification of the protocol conclusions, and specifically the first measurements of RHC and OHC activation energies; specifically for parallel RHC surface and gas-phase routes, and experimentally obtained OHC activation energy vs. first-principle calculations in the literature. Beyond insights related to the degreened catalyst studied here, the methodology provides a broadly available approach for quantifying how aging processes selectively impact RHC and OHC.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Damage development in neutron-irradiated concrete in a test reactor: Hygro-thermal and mechanical simulations

This paper reports the development of a 3D mesoscale hygro-thermal-mechanical simulation approach to predict damage in concrete irradiated in a test reactor. This framework, developed in MOOSE, considers the effects of elevated temperature, moisture content, and high neutron fluence (energy threshold, E > 0.1 MeV) on the mortar and aggregates separately. The first-stage simulation implements hygro-thermal analysis to determine the temperature and RH inside the specimen as a function of imposed radiation energy. These are used as inputs to the second stage, which considers radiation-induced volumetric expansion (RIVE) of aggregates, and creep, shrinkage, and stress-strain response of mortar to predict the expansion, stresses, and damage in specimens made using different coarse aggregates and subjected to different irradiation times. The irradiation time-dependent damage in the mortar is expressed using an isotropic damage parameter. This multi-physics model serves as a predictive tool for damage quantification in concrete due to neutron irradiation.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Portable chemical detection platform for on-site monitoring of odorant levels in natural gas

The adequate odorization of natural gas is critical to identify gas leaks and to reduce accidents. To ensure odorization, natural gas utility companies collect samples to be processed at core facilities or a trained human technician smells a diluted natural gas sample. In this work, we report a detection platform that addresses the lack of mobile solutions capable of providing quantitative analysis of mercaptans, a class of compounds used to odorize natural gas. Detailed description of the platform hardware and software components is provided. Designed to be portable, the platform hardware facilitates extraction of mercaptans from natural gas, separation of individual mercaptan species, and quantification of odorant concentration, with results reported at point-of-sampling. The software was developed to accommodate skilled users as well as minimally trained operators. Detection and quantification of six commonly used mercaptan compounds (ethyl mercaptan, dimethyl sulfide, n-propylmercaptan, isopropyl mercaptan, tert–butyl mercaptan, and tetrahydrothiophene) at typical odorizing concentrations of 0.1–5 ppm was performed using the device. We demonstrate the potential of this technology to ensure natural gas odorizing concentrations throughout distribution systems.

03 NATURAL GAS↗
DOE OSTI2022

Motif-centric phosphoproteomics to target kinase-mediated signaling pathways

Identifying cellular phosphorylation pathways based on kinase-substrate relationships is a critical step to understanding the regulation of physiological functions in cells. Mass spectrometry-based phosphoproteomics workflows have made it possible to comprehensively collect information on individual phosphorylation sites in a variety of samples. However, there is still no generic approach to uncover phosphorylation networks based on kinase-substrate relationships in rare cell populations. Here, we describe a motif-centric phosphoproteomics approach combined with multiplexed isobaric labeling, in which in vitro kinase reactions are used to generate targeted phosphopeptides, which are spiked into one of the isobaric channels to increase detectability. Proof-of-concept experiments demonstrate selective and comprehensive quantification of targeted phosphopeptides by using multiple kinases for motif-centric channels. More than 7,000 tyrosine phosphorylation sites were quantified from several tens of micrograms of starting materials. This approach enables the quantification of multiple phosphorylation pathways under physiological or pathological regulation in a motif-centric manner.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

How does uncertainty of soil organic carbon stock affect the calculation of carbon budgets and soil carbon credits for croplands in the U.S. Midwest?

Cropland carbon budget depicts the amount of carbon flowing in and out of agroecosystems and the changes in carbon stocks of soil and living biomass during the same period. Soil carbon credit is the additional change in soil carbon stock under certain farming practices compared with the business-as-usual practices. Accurately calculating cropland carbon budget and soil carbon credit is critical to assessing climate change mitigation potential in agroecosystems. The calculation of cropland carbon budget and soil carbon credit is sensitive to local soil and climatic conditions, especially initial soil organic carbon (SOC) stock, which is determined by both SOC concentration (SOC%) and bulk density (Bulk_Density). SOC stock data are either from soil sampling or gridded public survey data. In agroecosystem models, SOC stock data are a key model input for quantifying cropland carbon budget and soil carbon credit. However, various types and degrees of uncertainties exist in SOC stock datasets, which propagate to the quantification of SOC stock change. In particular, a large discrepancy is found in two widely used SOC stock datasets — Rapid Carbon Assessment dataset (RaCA) and Gridded Soil Survey Geographic Database (gSSURGO) — in the U.S. Midwest, with a relative difference (quantified using Normalized Root Mean Square Error, NRMSE) of 48.0% for 0–30 cm SOC stock between the two datasets. It remains largely unclear how uncertainty in SOC stocks affects the calculation of cropland carbon budget and soil carbon credit. To address this question, we used a well-validated process-based agroecosystem model, ecosys, to assess the impacts of SOC stock uncertainty on carbon budget and soil carbon credit calculation in the U.S. Midwestern corn-soybean rotation systems. Our results reveal the following findings: (1) A sizable discrepancy exists in simulated cropland carbon budget between using gSSURGO and using RaCA for their SOC% and Bulk_Density as model inputs, with a Pearson correlation coefficient (r) of only 0.4 for simulated change of SOC stock (ΔSOC) using these two different soil datasets. (2) Simulated cropland carbon budget components were more sensitive to initial SOC% than to Bulk_Density. For example, the upper and lower quartiles of multi-year averaged ΔSOC were –29.8 and 4.8 gC/m 2 /year for the selected counties respectively, with an uncertainty of 13.7 and 0.7 gC/m 2 /year induced by uncertainties in initial SOC% and Bulk_Density, respectively. (3) Both simulated ΔSOC and its uncertainty were negatively correlated with initial SOC%, whereas ΔSOC was negatively correlated with air temperature, and ΔSOC uncertainty was positively correlated with air temperature. (4) The uncertainty of calculated soil carbon credits was much smaller compared with the uncertainty of calculated absolute carbon budgets assuming the same SOC stock uncertainty level in the inputs. Specifically, in our assessment comparing planting cover crops vs no cover crop, the uncertainty of calculated soil carbon credits induced by initial SOC% uncertainty was less than 4% (relative to the quantified value of the soil carbon credits) for 90% of the cases. Our analysis highlights that high accuracy measurement of SOC% as inputs is needed for the calculation of cropland carbon budgets; however, soil carbon credit quantification is much less sensitive to the initial SOC% inputs, and the current publicly available soil datasets (e.g., gSSURGO) are largely suitable for the calculation of soil

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

CO 2 storage site characterization using ensemble-based approaches with deep generative models

Estimating spatially distributed properties such as permeability from available sparse measurements is a great challenge in efficient subsurface CO 2 storage operations. In this paper, a deep generative model that can accurately capture complex subsurface structure is tested with an ensemble-based inversion method for accurate and accelerated characterization of CO 2 storage sites. We chose Wasserstein Generative Adversarial Network with Gradient Penalty (WGAN-GP) for its realistic reservoir property representation and Ensemble Smoother with Multiple Data Assimilation (ES-MDA) for its robust data fitting and uncertainty quantification capability. WGAN-GP are trained to generate high-dimensional permeability fields from a low-dimensional latent space and ES-MDA then updates the latent variables by assimilating available measurements. Several subsurface site characterization examples including Gaussian, channelized, and fractured reservoirs are used to evaluate the accuracy and computational efficiency of the proposed method and the main features of the unknown permeability fields are characterized accurately with reliable uncertainty quantification. Furthermore, the estimation performance is compared with a widely-used variational, i.e., optimization-based, inversion approach, and the proposed approach outperforms the variational inversion method in several benchmark cases. We explain such superior performance by visualizing the objective function in the latent space: because of nonlinear and aggressive dimension reduction via generative modeling, the objective function surface becomes extremely complex while the ensemble approximation can smooth out the multi-modal surface during the minimization. This suggests that the ensemble-based approach works well over the variational approach when combined with deep generative models at the cost of forward model runs unless convergence-ensuring modifications are implemented in the variational inversion.

42 ENGINEERING↗
DOE OSTI2022

A combined experimental and modelling study of granite hydrothermal alteration

Geochemical reactions can induce significant changes of rock reservoir porosity and permeability via mineral dissolution and precipitation processes, affecting the long-term fluid behaviour within various geological systems. Here, the understanding and quantification of these reactions rely on field and experimental studies and on the predictions of reactive transport models. The present study was aimed at assessing the extent to which current geochemical models integrating available mineral dissolution/precipitation rate equations can reproduce the experimental data obtained from 4 to 17-day long hydrothermal alteration experiments of a muscovite-biotite granite and, thus, help provide an accurate description of the evolution of geothermal systems within granitic reservoirs. The experiments were conducted at a constant temperature of 180 °C and over an aqueous fluid pH range of 2 to 8.5, using both mixed-flow and static batch reactors. Modelled major element (K, Al, Si, Ca, and Mg) concentrations were generally in satisfactory agreement with the corresponding measured elemental fluxes – the differences between modelled and experimental values were generally within a factor of 5 – and the predicted identity and mass of formed secondary phases were consistent with the microscopic observations of the reacted solids. However, larger differences between measured and modelled element concentrations were observed when significant amounts of secondary phases formed, notably at pH 2 to 3, and for longer-term batch experiments. Much of this concentration difference stems from the underestimation of the amounts of Al-phases formed at acid to near-neutral pH. Although an idealized rock composition was considered, the observed mismatch between model calculations and experimental data can be attributed to inadequate mineral precipitation reaction rates and a poor description of reactive surface areas in existing geochemical modelling codes. More accurate quantification of precipitation kinetics, including nucleation and growth, and improved descriptions of the temporal change of mineral surface area would enhance the predictive capabilities of reactive transport models and benefit, particularly, the efforts aimed at increasing the sustainability of EGS reservoirs.

58 GEOSCIENCES↗
DOE OSTI2021

Crack nucleation from non-metallic inclusions in aluminum alloys described by peridynamics simulations

Conventional engineering methods oftentimes have challenges in the quantification of crack nucleation processes from manufacturing defects that are relevant for engineering component lifing. In this work, we present peridynamic simulation framework for the description of the crack nucleation process from cluster of non-metallic inclusions in aluminum alloys. Our non-local simulation framework characterizes crack nucleation process as multiple micro-crack nucleation events from individual inclusions, and eventually one micro-crack dominates. We define individual stages of the crack nucleation process, i.e., nucleation, micro-crack, technical, and crack initiation, that allow a quantification and meta model development of both the individual stages and the entire crack nucleation process.

36 MATERIALS SCIENCE↗
DOE OSTI2025

ARCH: Large-scale knowledge graph via aggregated narrative codified health records analysis

Objective: Electronic health record (EHR) systems contain a wealth of clinical data stored as both codified data and free-text narrative notes (NLP). The complexity of EHR presents challenges in feature representation, information extraction, and uncertainty quantification. Here, to address these challenges, we proposed an efficient Aggregated naRrative Codified Health (ARCH) records analysis to generate a large-scale knowledge graph (KG) for a comprehensive set of EHR codified and narrative features. Methods: Using data from 12.5 million Veterans Affairs patients, ARCH first derives embedding vectors and generates similarities along with associated p-values to measure the strength of relatedness between clinical features with statistical certainty quantification. Next, ARCH performs a sparse embedding regression to remove indirect linkage between features to build a sparse KG. Finally, ARCH was validated on various clinical tasks, including detecting known relationships between entity pairs, predicting drug side effects, disease phenotyping, as well as sub-typing Alzheimer’s disease patients. Results: ARCH produces high-quality clinical embeddings and KG for over 60,000 codified and narrative EHR concepts. The KG and embeddings are visualized in the R-shiny powered web-API.3 ARCH achieved high accuracy in detecting EHR concept relationships, with AUCs of 0.926 (codified) and 0.861 (NLP) for similar EHR concepts, and 0.810 (codified) and 0.843 (NLP) for related pairs. It detected drug side effects with a 0.723 AUC, which improved to 0.826 after fine-tuning. Using both codified and NLP features, the detection power increased significantly. Compared to other methods, ARCH has superior accuracy and enhances weakly supervised phenotyping algorithms’ performance. Notably, it successfully categorized Alzheimer’s patients into two subgroups with varying mortality rates. Conclusion: The proposed ARCH algorithm generates large-scale high-quality semantic representations and knowledge graph for both codified and NLP EHR features, useful for a wide range of predictive modeling tasks.

Electronic health records↗
DOE OSTI2021

B-PINNs: Bayesian Physics-informal Neural Networks for Forward and Inverse PDE Problems with Noisy Data

We propose a Bayesian physics-informed neural network (B-PINN) to solve both forward and inverse nonlinear problems described by partial differential equations (PDEs) and noisy data. In this Bayesian framework, the Bayesian neural network (BNN) combined with a PINN for PDEs serves as the prior while the Hamiltonian Monte Carlo (HMC) or the variational inference (VI) could serve as an estimator of the posterior. B-PINNs make use of both physical laws and scattered noisy measurements to provide predictions and quantify the aleatoric uncertainty arising from the noisy data in the Bayesian framework. Compared with PINNs, in addition to uncertainty quantification, B-PINNs obtain more accurate predictions in scenarios with large noise due to their capability of avoiding overfitting. We conduct a systematic comparison between the two different approaches for the B-PINNs posterior estimation (i.e., HMC or VI), along with dropout used for quantifying uncertainty in deep neural networks. Our experiments show that HMC is more suitable than VI with mean field Gaussian approximation for the B-PINNs posterior estimation, while dropout employed in PINNs can hardly provide accurate predictions with reasonable uncertainty. Finally, we replace the BNN in the prior with a truncated Karhunen-Loève (KL) expansion combined with HMC or a deep normalizing flow (DNF) model as posterior estimators. The KL is as accurate as BNN and much faster but this framework cannot be easily extended to high-dimensional problems unlike the BNN based framework.

Non linear PDEs, Noisy data, Bysian physics inform↗
DOE OSTI2025

Microstructurally informed synchrotron x-ray analysis revealing helium defect transitions in ultrafine grained tungsten

The formation of insoluble gaseous defects in materials due to nuclear transmutation or ion implantation involves the diffusion of impurity atoms to form atomic defect clusters that coalesce into bubbles or cavities and ultimately degrade the material properties. Transmission electron microscopy (TEM) is limited in its ability to resolve sub-nanometer gas clusters whereas X-ray diffraction (XRD) provides information pertaining to local atomic changes. Here, in this study, helium (He) implanted ultrafine grained tungsten is explored through a multimodal defect characterization campaign combining TEM-informed Small Angle X-ray Scattering (SAXS) analysis, XRD lattice parameter measurements, and nanoscale He cluster quantification from a region of reciprocal space accessible via Wide Angle X-ray Scattering (WAXS). Moderate elevated temperature implantations are shown to produce high concentrations of sub-nanoscale He clusters and small, homogeneously distributed cavities, which collectively are linked to lattice expansion and further substantiated through complementary atomistic simulations. Increased implantation temperatures encourage the diffusion of these defects to the grain boundaries (GBs), leading to lattice relaxation and the growth of large GB cavities manifesting as bimodal size distributions in the SAXS analysis. Overall, our results demonstrate the utility of multimodal synchrotron X-ray analysis in bridging the gap between microscale He cavity quantification and atomic-scale defect analysis.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Analysis of phospholipids and triacylglycerols in intravenous lipid emulsions

Intravenous lipid emulsions (ILEs) are used for parenteral nutrition, providing a vital source of essential fatty acids and concentrated energy for patients who are unable to absorb nutrients via the digestive track. They are commonly used to treat local and non-local anesthetic toxicity, and lipophilic drug overdose. ILE are composed of natural lipids, and the composition of these natural lipids can be varied based on their source. The lipids are susceptible to hydrolytic degradation with time, resulting various lipid degradation products such as Lysophosphatidylcholines (LPs), affecting the actual composition of nutrients in the formulation. As a result, the identification and quantification of lipid components, including degradation products, in ILEs are crucial in quality control. In this study, lipids from different batches of ILE Intralipid® 20%, were separated and identified using a UHPLC-ESI-QTOF system and SimLipid® high throughput lipid identification software. Out of 47 lipids identified, 34 were phospholipids (PLs) and the others were triacylglycerols (TAGs). Most of the phospholipids detected were phosphatidylcholines (PC) and Lysophosphatidylcholines (LPC). A total of 9 LPCs, 18 PCs, 6 phosphoethanolamines (PEs), and 1 sphingomyelin (SM) were identified. The LPCs concentration changed with the manufacturing date and storage time. Furthermore, this UHPLC method enabled the identification and quantification of lipids and their decomposition products in complex ILE emulsion mixtures on a single 20-minute chromatographic run.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

An infrared spectral database for gas-phase quantitation of volatile per- and polyfluoroalkyl substances (PFAS)

We report the construction of a database of vetted infrared spectra specifically targeting volatile fluorocarbon gases that may be emitted during thermal treatment of per- and polyfluoroalkyl substances (PFAS) to assist understanding of treatment processes and improve quantification. To populate this database, protocols derived from the Pacific Northwest National Laboratory (PNNL) infrared spectral database are used, curtailing the species selection for this data set. Each spectrum in the database is a weighted average derived from 10 or more individual measurements at different partial pressures (static method) or flow rates (dissemination method) to yield good fidelity of both strong and weak infrared signatures, with each composite spectrum ranging from = 6500 cm-1 to = 600 cm-1 with an apodized resolution of 0.112 cm-1. This resolution was chosen to fully resolve all spectral features, recognizing that atmospheric pressure broadening results in nearly all ro-vibrational lines having linewidths = 0.1 cm-1. As an example case, application of the database is demonstrated via identification and quantification of dominant 1H-perfluoroheptane and perfluorohept-1-ene fluorocarbon products resulting from thermal decomposition of perfluorooctanoate (PFOA) below 450 °C.

Infrared, Gas-phase spectra, FTIR, Spectral databa↗
DOE OSTI2023

Multiscale spatial analysis of fracture nodes in two dimensions

Spatial arrangement of fractures as a function of scale is an important component of fracture quantification for inferential and predictive modeling. Available methods that analyze fracture spatial arrangement are based on one-dimensional spacing data; therefore, they are limited to semi-parallel fractures. Such methods cannot be applied to fracture networks in higher dimensions, particularly when fractures have different orientations. Here to characterize fracture arrangements in two dimensions, we propose using Ripley’s K-function, as a method of point pattern analysis, to quantify spatial arrangement of fracture nodes. Fracture nodes, such as barycenters, intersection points, and tips, are point-based representations of fracture locations and connectivity within the fracture network. We introduce formulations for isotropic as well as directional analyses of spatial arrangement. In addition, we derive formulations for edge correction in circular and rectangular study domains. Finally, we demonstrate applications of Ripley’s K-function on two natural fracture datasets. Our proposed method supports quantification and characterization of fracture spatial arrangements that allow practitioners to build representative models of fractures in the subsurface.

02 PETROLEUM↗