Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Geologic Modeling”

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,009 records · Page 56
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Impact of cement composition, brine concentration, diffusion rate, reaction rate and boundary condition on self-sealing predictions for cement-CO 2 systems

Geological CO 2 storage (GCS) plays an important role in curbing CO 2 emissions by reducing the carbon footprint of difficult to decarbonize operations and achieve negative CO 2 emissions through activities like Bioenergy with Carbon Capture and Storage (BECCS) and Direct Air Carbon Capture and Storage (DACCS). Leakage of CO 2 through wells is an important concern when it comes to deployment of large-scale GCS. Existing wells in sites that are otherwise suitable for GCS can act as conduits for stored CO 2 to escape the reservoir. There is broad consensus that the main risk of leakage through wellbores is via fractures/damaged pathways. The results from several studies evaluating the permeability evolution of cement fractures in wells upon leakage of CO 2 agree that smaller fracture apertures, slower brine velocities and higher brine residence times promote self-sealing of fractures by mineral precipitation. Quantitatively, however, the differences in sealing conditions are significant and are typically attributed to differences in experimental conditions or model assumptions. Here we examine the sensitivity of our model, describing CO 2 leakage through wellbores, to cement composition, brine concentration, diffusion rates, and reaction rates. We also evaluate the impact of the boundary condition to allow comparisons between observations from experiments performed at constant flow rate and model predictions made at constant pressure conditions. Our results show that diffusion and reactions rates have the most impact on the self-sealing criteria for cement-CO 2 systems. In addition, conditions associated with self-sealing of fractures at constant flow rate require longer fractures, smaller fracture apertures, and slower velocities than under constant pressure.

58 GEOSCIENCES↗
DOE OSTI2022

Insights into Controls of Mineralogy and Pore Structure on the Density of Methane Adsorption Phase in Shales under Supercritical Conditions

With gradually increasing enhancing population and industrialization, exploitation of shale gas is enhancing in the world to satisfy the growing demand for energy worldwide. To exploit shale gas efficiently, accurate assessment of shale gas in place (GIP) is necessary for determining production strategies. Unlike free and dissolved gas, adsorbed gas contributes to the shale GIP up to 85% due to the well-developed surface area and micropores. Measured excess adsorption of methane (up to over 95% of subsurface shale gas) is often corrected to absolute adsorption to obtain the actual amount of adsorption under supercritical (or geological) conditions. During the correction, adsorption phase density (APD) is critical. However, the APD of subsurface shale gas and the effects of shale properties (e.g., mineralogy and pore structure) on APD remain poorly understood. A series of high-pressure methane adsorption isotherms (HPMAI) on Caney Shales were collected and analyzed in conjunction with other United States and Chinese shales from the literature at 35–125 °C and up to 15 MPa. Here, a three-layer Ono–Kondo (OK3) model is utilized to derive the temperature- and pressure-dependent APD coupling low-pressure nitrogen adsorption isotherms (LPNAI) and HPMAI. X-ray diffraction and organic geochemistry are combined to reveal the mineralogy. Brunauer–Emmett–Teller, Barrett–Joyner–Halenda, and Horvath–Kawazoe analyses via LPNAI are used to investigate the pore structures. Results show that APD increases with organic matter (OM) proxied by total organic carbon and decreases with clay minerals and the sum of quartz and feldspar. OM dramatically contributes to the APD as multiple-layer adsorption exists, and the APD for OM could be 1.4–8.5 times that for clay minerals. Other inorganic minerals contribute less to APD. The properties and constitution of the surface area instead of the volume fraction contribute to the APD in shales. APD does not show an obvious correlation with micropore volume, likely related to the ratio of micropore volume to the total pore volume. Here, we provide a significant and comprehensive study of petrological factors that impact the APD of subsurface shale gas, which will improve the estimation of supercritical adsorption and shale GIP under reservoir conditions. Also, the findings in this work can provide applications for subsurface carbon dioxide adsorption and storage.

58 GEOSCIENCES↗
DOE OSTI2020

The impact of mineral reactive surface area variation on simulated mineral reactions and reaction rates

Reactive transport modeling is an essential tool to simulate complex geochemical reactions in porous media that can impact formation properties including porosity and permeability. However, simulating these reactions is challenging due to uncertainties in model parameters, particularly mineral surface areas. Imaging has emerged as a powerful means of estimating model parameters including porosity and mineral abundance, accessibility and accessible surface area. However, these parameters, particularly mineral accessible surface area, vary with image resolution. This work aims to enhance understanding of the impact of image resolution and other means of estimating mineral reactive surface area on simulated mineral reactions and reaction rates. Mineral surface areas calculated from images with resolutions of 0.34 μm and 5.71 μm were used to simulate mineral reactions in the context of geologic CO 2 sequestration in the Paluxy formation at the continuum scale. Additional simulations were carried out using BET surface areas collected from the literature and geometric surface areas. Simulations were run for 7300 days and mineral volume fractions and effluent ion concentrations tracked and compared. Variations in mineral surface areas measured from images are within 1 order of magnitude and yield similar simulation results, indicating the impact of image resolution on simulated reactions and reaction rates is minimum for the resolutions and sample considered. In comparison, surface areas obtained from BET and geometric approaches are 1–4 orders of magnitude higher than image-obtained surface areas and result in greater simulated reaction rates and extents. Minerals with high reaction rates (calcite and siderite) are most impacted by surface area values at short times where simulated mineral volume fractions at longer times agree relatively well, even for simulations with several orders of magnitude variation in surface area. Phases with lower reaction rates, such as K-feldspar and muscovite, are predominantly impacted over longer times where variations in surface areas impact reaction extents and porosity evolution. Overall, variations in surface areas due to image resolution are small and result in little variation in simulated reactions and reaction rates while there are significant variations in simulation results when other surface area estimates are used. These variations, however, depend on the reactivity of the mineral phase where surface areas of fast-reacting phases largely impact simulations at short (hours to days) timescales and surface areas of slower-reactive phases impact longer term simulations (years).

58 GEOSCIENCES↗
DOE OSTI2021

Significance and complications of drilling fluid rheology in geothermal drilling: A review

The harsh downhole conditions of high pressure and high temperature (HPHT) encountered in geothermal wells make the drilling operation challenging. Drilling in such environments requires a special drilling mud formulation with high thermal stability and good rheological properties to fulfill the drilling fluid functions. Therefore, great efforts should be put into selecting the suitable drilling fluid, optimize and monitor the drilling fluid properties throughout drilling operations, and predicting its performance under downhole conditions. Rheological properties significantly impact many drilling parameters such as hole cleaning, fluid and wellbore stability, wellbore hydraulics, torque and drag, and other drilling issues. This paper discusses water-based drilling fluids' flow behavior under HPHT conditions and highlights the significance of fluid rheology in geothermal drilling. The common challenges and complications related to fluid rheology encountered in geothermal drilling are addressed in this paper, such as hole cleaning, wellbore hydraulics, and drilling fluid stability. This article also reviews the recent advances in drilling mud systems, rheology enhancement, and rheological properties measurements at surface and subsurface conditions. Furthermore, the rheology models of drilling fluid at elevated temperatures are reviewed to fully understand their flow behavior and establish a method for drilling engineers to optimize fluid formulations for geothermal drilling.

58 GEOSCIENCES↗
DOE OSTI2023

Predicting Short-Term Deformation in the Central Valley Using Machine Learning

Land subsidence caused by excessive groundwater pumping in Central Valley, California, is a major issue that has several negative impacts such as reduced aquifer storage and damaged infrastructures which, in turn, produce an economic loss due to the high reliance on crop production. This is why it is of utmost importance to routinely monitor and assess the surface deformation occurring. Two main goals that this paper attempts to accomplish are deformation characterization and deformation prediction. The first goal is realized through the use of Principal Component Analysis (PCA) applied to a series of Interferomtric Synthetic Aperture Radar (InSAR) images that produces eigenimages displaying the key characteristics of the subsidence. Water storage changes are also directly analyzed by the use of data from the Gravity Recovery and Climate Experiment (GRACE) twin satellites and the Global Land Data Assimilation System (GLDAS). The second goal is accomplished by building a Long Short-Term Memory (LSTM) model to predict short-term deformation after developing an InSAR time series using LiCSBAS, an open-source InSAR time series package. The model is applied to the city of Madera and produces better results than a baseline averaging model and a one dimensional convolutional neural network (CNN) based on a mean squared error metric showing the effectiveness of machine learning in deformation prediction as well as the potential for incorporation in hazard mitigation models. The model results can directly aid policy makers in determining the appropriate rate of groundwater withdrawal while maintaining the safety and well-being of the population as well as the aquifers’ integrity.

58 GEOSCIENCES↗
DOE OSTI2021

Simulating measurable ecosystem carbon and nitrogen dynamics with the mechanistically defined MEMS 2.0 model

Abstract. For decades, predominant soil biogeochemical models have used conceptual soil organic matter (SOM) pools and only simulated them to a shallow depth in soil. Efforts to overcome these limitations have prompted the development of the new generation SOM models, including MEMS 1.0, which represents measurable biophysical SOM fractions, over the entire root zone, and embodies recent understanding of the processes that govern SOM dynamics. Here we present the result of continued development of the MEMS model, version 2.0. MEMS 2.0 is a full ecosystem model with modules simulating plant growth with above- and belowground inputs, soil water and temperature by layer, decomposition of plant inputs and SOM, and mineralization and immobilization of nitrogen (N). The model simulates two commonly measured SOM pools – particulate and mineral-associated organic matter (POM and MAOM, respectively). We present results of calibration and validation of the model with several grassland sites in the US. MEMS 2.0 generally captured the soil carbon (C) stocks (R2 of 0.89 and 0.6 for calibration and validation, respectively) and their distributions between POM and MAOM throughout the entire soil profile. The simulated soil N matches measurements but with lower accuracy (R2 of 0.73 and 0.31 for calibration and validation of total N in SOM, respectively) than for soil C. Simulated soil water and temperature were compared with measurements, and the accuracy is comparable to the other commonly used models. The seasonal variation in gross primary production (GPP; R2 = 0.83), ecosystem respiration (ER; R2 = 0.89), net ecosystem exchange (NEE; R2 = 0.67), and evapotranspiration (ET; R2 = 0.71) was well captured by the model. We will further develop the model to represent forest and agricultural systems and improve it to incorporate new understanding of SOM decomposition.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Estimation of hydraulic conductivity in a watershed using sparse multi-source data via Gaussian process regression and Bayesian experimental design

Enhanced water management systems depend on accurate estimation of subsurface hydraulic properties. However, geologic formations can vary significantly, so information from a single source (e.g., widely spaced boreholes) is insufficient in characterizing subsurface aquifer properties. Therefore, multiple sources of information are needed to complement the hydrogeology understanding of a region. Here, this study presents a numerical framework in which information from different measurement sources is combined to characterize the 3D random field in a multi-fidelity prediction model. Coupled with the model, a Bayesian experimental design was used to determine the best future sampling locations. The Upper Sangamon watershed in east-central Illinois was selected as the case study site, where the multi-fidelity Gaussian process model was used to estimate the hydraulic conductivity in the region of interest. Multi-source observation data were obtained from electrical resistivity and borehole pumping tests. The accuracy of the model prediction is dependent on the locations and the distribution of both high- and low-fidelity data. Furthermore, the multi-fidelity model was compared with the single-fidelity model. The uncertainties and confidence in the measurements and parameter estimates were quantified and used to design future cycles of data collection to further improve the confidence intervals.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Deep learning multiphysics network for imaging CO 2 saturation and estimating uncertainty in geological carbon storage

Multiphysics inversion exploits different types of geophysical data that often complement each other and aims to improve overall imaging resolution and reduce uncertainties in geophysical interpretation. Despite the advantages, traditional multiphysics inversion is challenging because it requires a large amount of computational time and intensive human interactions for preprocessing data and finding trade-off parameters. These issues make it nearly impossible for traditional multiphysics inversion to be applied as a real-time monitoring tool for geological carbon storage. In this paper, we present a deep learning (DL) multiphysics network for imaging CO 2 saturation in real time. The multiphysics network consists of three encoders for analysing seismic, electromagnetic and gravity data and shares one decoder for combining imaging capabilities of the different geophysical data for better predicting CO 2 saturation. The network is trained on pairs of CO 2 label models and multiphysics data so that it can directly image CO 2 saturation. Here we use the bootstrap aggregating method to enhance the imaging accuracy and estimate uncertainties associated with CO 2 saturation images. Using realistic CO 2 label models and multiphysics data derived from the Kimberlina CO 2 storage model, we evaluate the performance of the deep learning multiphysics network and compare its imaging results to those from the deep learning single-physics networks. Our modelling experiments show that the deep learning multiphysics network for seismic, electromagnetic, and gravity data not only improves the imaging accuracy but also reduces uncertainties associated with CO 2 saturation images. Our results also suggest that the deep learning multiphysics network for the non-seismic data (i.e., electromagnetic and gravity) can be used as an effective low-cost monitoring tool in between regular seismic monitoring.

58 GEOSCIENCES↗
DOE OSTI2021

The effects of well damage and completion designs on geoelectrical responses in mature wellbore environments

Well integrity is one of the major concerns in long-term geologic storage sites due to the potential risk of well leakage and groundwater contamination. Evaluating changes in electrical responses due to energized steel-cased wells has the potential to quantify and predict possible wellbore failures because any kind of breakage or corrosion along highly conductive well casings will have an impact on the distribution of the subsurface electrical potential. However, realistic wellbore-geoelectrical models that can fully capture fine-scale details of well completion design and the state of well damage at the field scale require extensive computational effort, or they can even be intractable to simulate. To overcome this computational burden while still keeping the model realistic, we have used the hierarchical finite-element method that represents electrical conductivity at each dimensional component (1D edges, 2D planes, and 3D cells) of a tetrahedral mesh. This allows well completion designs with real-life geometric scales and well systems with realistic, detailed, progressive corrosion and damage in our models. We have developed a comparison of possible discretization approaches of a multicasing completion design in the finite-element model. The effects of the surface casing and the coupling between concentric well casings as well as the effects of the degree and the location of well damage on the electrical responses are also examined. As a result, we analyze real surface electric field data to detect wellbore integrity failure associated with damage.

58 GEOSCIENCES↗
DOE OSTI2022

Underground sources of drinking water chemistry changes in response to potential CO 2 leakage

The purpose of this study was to quantify changes to underground sources of drinking water (USDW) quality in response to potential CO 2 leakage from geologic CO 2 sequestration (GCS) reservoirs. Here, the developed a framework of combined laboratory experiments and reactive transport simulations and used this framework to evaluate the Ogallala aquifer overlying the Farnsworth Unit (FWU), an active GCS site, as a case study. Using chemical reaction parameters obtained from laboratory experiments and numerical simulations, site-specific mechanisms of CO 2 -water-sediment interactions at the USDW aquifer were interpreted. Long-term risks of potential CO 2 leakage were then evaluated with field-scale numerical models using the regional hydrogeological characteristics and reaction parameters obtained from our experiments and simulations. Results suggest that carbonate mineral impurity and cation exchange are key mechanisms for interactions between CO 2 and the aquifer sediment. Additionally, for a large leakage rate of 0.1 % injection from one leaky well, the leakage plume might impact an area of 300 m in diameter and significantly affect the local water quality by changing pH and cation concentrations (e.g., Zn, Ba and Sr). After leakage ceases, the zone of impacted fluids would not migrate significantly in subsequent decades due to a low regional groundwater flowrate (for this case study). The relatively small area of impact might not be detected in a monitoring well given the broader spacing in a typical field scenario. Effective early leakage detection may require additional tools, e.g., borehole CO 2 movement, four-dimensional seismicity, CO 2 soil flux, samples from deeper aquifers, etc., to ensure effective leakage detection and long-term safety of GCS projects.

58 GEOSCIENCES↗
DOE OSTI2025

A flow-through gamma detection system for sub-surface radioxenon tracer measurements in LYNM PE1

A continuous flow-through, real-time radioxenon sensor array has been designed to monitor either 127 Xe or 133 Xe tracer gas as part of a series of underground chemical explosions in the Low-Yield Nuclear Monitoring (LYNM) Physics Experiment 1 (PE1). The sensors are pairs of shielded thallium-doped sodium iodide gamma-ray detectors that surround a gas measurement chamber. Each measurement chamber is connected to a sampling location in the geology surrounding the chemical explosion or in the access tunnel. The first experiment in the series, PE1 A, occurred in 2023 and the system quantified the 127 Xe tracer gas in the gas sampling lines for four weeks following the experiment execution. The array generated a set of radioxenon concentration time series with a 15-minute resolution that will be used to refine subsurface gas transport models.

Gamma-ray coincidence detection↗
DOE OSTI2022

Downwind dispersion of CO 2 from a major subsea blowout in shallow offshore waters

Growing interest in offshore geologic carbon sequestration (GCS) motivates risk assessment of large-scale subsea CO 2 well blowouts or pipeline ruptures. For major leaks of CO 2 from wells or pipelines, significant fluxes of CO 2 may occur from the sea surface depending on water depth. In the context of risk assessment of human health and safety, we have used previously simulated coupled well-reservoir and water column model results as a source term for dense gas dispersion of CO 2 above the sea surface. The models are linked together by one-way coupling, that is, output of one model is used as input to the next model. These first-of-their-kind coupled flow results are applicable to assessing the hazard of CO 2 to people at and downwind of the sea surface location of emission. Hazard is quantified by plotting the downwind dispersion length (DDL), which we define in the study as the distances from the emission source to the point at which the emitted CO 2 has been diluted to 5% and 1.5% in air by volume. Here results suggest that large-scale blowouts in shallow water (10 m) may cause hazardous CO 2 plumes extending on the order of several hundred meters downwind. Details of the modeling show DDL has a maximum for windspeed (at an elevation of 10 m) of approximately 5 m/s, with smaller DDL for both weaker and stronger winds. This is explained by the fact that wind favors transport but also causes dispersion; therefore there is a certain wind speed that maximizes DDL.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

The Effect of DPC Fillers on FEPs Relevant to Disposal of SNF

The US Department of Energy (DOE) is investigating the use of different materials that could be used to fill the void space inside a dual-purpose canister (DPC) loaded with spent nuclear fuel (SNF) just before it is emplaced in a deep geologic repository. The purpose of adding filler material is to maintain subcritical conditions in the repository during the postclosure period, which can span up to 1,000,000 years. Several types of materials have been proposed, including metals, cements, particulates, and glass. Part of this investigation addresses how the presence of filler material inside a DPC will affect the performance of the repository with respect to the repository features; the consequences of events that may occur; and the multiple thermal, hydrologic, chemical, and mechanical processes that may occur in a deep geologic repository over long timescales. This report describes some of the filler materials that have been proposed and studied; identifies 11 features, 6 events, and 25 processes that may be affected by the presence of filler materials; and discusses the effects that may require consideration for each feature, event, or process. The results of this study can be used to direct appropriate research and to develop suitable models if the DOE decides to use fillers to maintain subcritical conditions in DPCs used to dispose of SNF.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

PYSIMFRAC: A Python library for synthetic fracture generation and analysis

In this paper, we introduce PYSIMFRAC, an open-source python library for generating 3-D synthetic fracture realizations, integrating with fluid simulators, and performing analysis. PYSIMFRAC allows the user to specify one of three fracture generation techniques (Box, Gaussian, or Spectral) and perform statistical analysis including the autocorrelation, moments, and probability density functions of the fracture surfaces and aperture. This analysis and accessibility of a python library allows the user to create realistic fracture realizations and vary properties of interest. In addition, PYSIMFRAC includes integration examples to two different pore-scale simulators and the discrete fracture network simulator, dfnWorks. The capabilities developed in this work provides opportunity for quick and smooth adoption and implementation by the wider scientific community for accurate characterization of fluid transport in geologic media. We present PYSIMFRAC along with integration examples and discuss the ability to extend PYSIMFRAC from a single complex fracture to complex fracture networks.

58 GEOSCIENCES↗
DOE OSTI2024

In Situ Velocity‐Strain Sensitivity Near the San Jacinto Fault Zone Analyzed Through Train Tremors

Abstract We utilize train tremors as P‐wave seismic sources to investigate velocity‐strain sensitivity near the San Jacinto Fault Zone. A dense nodal array deployed at the Piñon Flat Observatory is used to detect and identify repeating train energy emitted from a railway in the Coachella valley. We construct P‐wave correlation functions across the fault zone and estimate the spatially averaged dt / t versus strain sensitivity to be 6.25 × 10 4 . Through numerical simulations, we explore how the sensitivity decays exponentially with depth. The optimal solution reveals a subsurface sensitivity of 1.2 × 10 5 and a depth decay rate of 0.05 km −1 . This sensitivity aligns with previous findings but is toward the higher end, likely due to the fractured fault‐zone rocks. The depth decay rate, previously unreported, is notably smaller than assumed in empirical models. This raises the necessity of further investigations of this parameter, which is crucial to study stress and velocity variations at seismogenic depth.

Geology↗
DOE OSTI2022

Using a Multi‐Institutional Ensemble of Watershed Models to Assess Agricultural Conservation Effectiveness in a Future Climate

This study investigates the combined impacts of climate change and agricultural conservation on the magnitude and uncertainty of nutrient loadings in the Maumee River Watershed, the second-largest watershed of the Laurentian Great Lakes. Two scenarios — baseline agricultural management and increased agricultural conservation — were assessed using an ensemble of five Soil and Water Assessment Tools driven by six climate models. The increased conservation scenario included raising conservation adoption rates from a baseline of existing conservation practices to feasible rates in the near future based on farmer surveys. This increased adoption of winter cover crops on 6%–10% to 60% of cultivated cropland; subsurface placement of phosphorus fertilizers on 35%–60% to 68% of cultivated cropland; and buffer strips intercepting runoff from 29%–34% to 50% of cultivated cropland. Increased conservation resulted in statistically significant (p ≤ 0.05) reductions in annual loads of total phosphorus (41%), dissolved reactive phosphorus (18%), and total nitrogen (14%) under the highest emission climate scenario (RCP 8.5). While nutrient loads decreased with increased conservation relative to baseline management for all watershed models, different conclusions on the true effectiveness of conservation under climate change may be drawn if only one watershed model was used.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Surface Deformation and Seismicity Linked to Fluid Injection in the Raton Basin

Abstract It is suggested that in addition to seismicity deep fluid injection may cause surface uplift and subsidence in oil and gas‐producing regions. This study uses the Raton Basin as an example to investigate the hydromechanical processes of surface uplift and subsidence during wastewater injection. The Raton Basin, in southern central Colorado and northern central New Mexico, has experienced wastewater injection related to coalbed methane and gas production starting in 1994. In this study, we estimate the extent and magnitude of total vertical deformation in the Raton Basin from 1994 to 2020 and incremental deformation between the years 2017 to 2020. Results indicate a modeled uplift as much as 15 cm occurring between 1994 and 2020. Between 2017 and 2020, up to 3 cm of uplift occurred, largely near the northwestern injection wells. Most modeled uplift between 1994 and 2020 occurred near the southern wells, where the greatest cumulative volume of wastewater was injected. However, modeled subsidence occurred around the southern and eastern wells between 2017 and 2020, after the rate of injection decreased. Modeling indicates that while the magnitude of modeled uplift corresponds to cumulative injection volume and maximum rate in the long‐term, short‐term incremental deformation (uplift or subsidence) is controlled by changes in the rate of injection. The number of yearly earthquake events follows periods of rapid modeled uplifting throughout the Basin, suggesting that measurable surface deformation may be caused by the same injection‐induced pore pressure perturbations that initiate seismicity.

Geology↗
DOE OSTI2021

Risk assessment of engineering diseases of embankment–bridge transition section for railway in permafrost regions

Abstract The embankment–bridge transition section (EBTS) is one of the zones where railway diseases occur frequently in permafrost regions. Disease risk assessment of EBTSs can provide guidance for maintenance. In this study, considering the engineering geological conditions, climate characteristics, and embankment structure types along the Qinghai–Tibet Railway (QTR) as well as based on the disease inventory of the QTR from 2010 to 2019, the logistic regression (LR), support vector machine (SVM), and combination‐weight‐based gay relation analysis (GRA) were used for disease risk assessment of the EBTSs along the QTR in permafrost regions. The results indicate that the LR and SVM models have a better capability for EBTS disease prediction than the GRA model, and the SVM model can select more disease samples in relatively larger regions than the LR model. Based on the SVM and LR models, the risk level of EBTSs is divided into four classes: low‐ (29.9%), moderate‐ (39.6%), high‐ (22.1%), and very high (8.4%) risk. Finally, we selected 272 EBTSs in high‐ and very‐high‐risk classes for key observation during the maintenance of the QTR in permafrost regions. This study provides a reference for the risk assessment of railways built in permafrost regions using data‐driven methods.

Zhang, Saize↗