Geostatistical seismic inversion and 3D modelling of metric flow units,porosity and permeability in Brazilian presalt reservoir

Flow units(FU)rock typing is a common technique for characterizing reservoir flow behavior,producing reliable porosity and permeability estimation even in complex geological settings.However,the lateral extrapolation of FU away from the well into the whole reservoir grid is commonly a difficult task and using the seismic data as constraints is rarely a subject of study.This paper proposes a workflow to generate numerous possible 3D volumes of flow units,porosity and permeability below the seismic resolution limit,respecting the available seismic data at larger scales.The methodology is used in the Mero Field,a Brazilian presalt carbonate reservoir located in the Santos Basin,who presents a complex and heterogenic geological setting with different sedimentological processes and diagenetic history.We generated metric flow units using the conventional core analysis and transposed to the well log data.Then,given a Markov chain Monte Carlo algorithm,the seismic data and the well log statistics,we simulated acoustic impedance,decametric flow units(DFU),metric flow units(MFU),porosity and permeability volumes in the metric scale.The aim is to estimate a minimum amount of MFU able to calculate realistic scenarios porosity and permeability scenarios,without losing the seismic lateral control.In other words,every porosity and permeability volume simulated produces a synthetic seismic that match the real seismic of the area,even in the metric scale.The achieved 3D results represent a high-resolution fluid flow reservoir modelling considering the lateral control of the seismic during the process and can be directly incorporated in the dynamic characterization workflow.

Challenges in predictive modelling of chronic kidney disease:A narrative review

The exponential rise in the burden of chronic kidney disease(CKD)worldwide has put enormous pressure on the economy.Predictive modeling of CKD can ease this burden by predicting the future disease occurrence ahead of its onset.There are various regression methods for predictive modeling based on the distribution of the outcome variable.However,the accuracy of the predictive model depends on how well the model is developed by taking into account the goodness of fit,choice of covariates,handling of covariates measured on a continuous scale,handling of categorical covariates,and number of outcome events per predictor parameter or sample size.Optimal performance of a predictive model on an independent cohort is desired.However,there are several challenges in the predictive modeling of CKD.Disease-specific methodological challenges hinder the development of a predictive model that is cost-effective and universally applicable to predict CKD onset.In this review,we discuss the advantages and challenges of various regression models available for predictive modeling and highlight those best for future CKD prediction.

A practical guideline for univariate meta-analytic structural equation modelling

Background:Meta-analysis is a quantitative approach that systematically integrates results from previous research to draw conclusions.Structural equation modelling is a statistical method that integrates factor analysis and path analysis.Meta-analytic structural equation modeling(MASEM)combines meta-analysis and structural equation modeling.It allows researchers to explain relationships among a group of variables across multiple studies.Methods:We used a simulated dataset to conduct a univariate MASEM analysis,using Comprehensive Meta Analysis 3.3,Analysis of Moment Structures 24.0 software.Results:Despite the lack of concise literature on the methodology,our study provided a practical step-by-step guide on univariate MASEM.Conclusion:Researchers can employ MASEM analysis in applicable fields based on the description,principles,and practices expressed in this study and our previous publications mentioned in this study.

Modelling of internal ballistics of gun systems:A review

A deep understanding of the internal ballistic process and the factors affecting it is of primary importance to efficiently design a gun system and ensure its safe management. One of the main goals of internal ballistics is to estimate the gas pressure into the combustion chamber and the projectile muzzle velocity in order to use the propellant to its higher efficiency while avoiding over-pressure phenomena. Dealing with the internal ballistic problem is a complex undertaking since it requires handling the interaction between different constituents during a transient time lapse with very steep rise of pressure and temperature. Several approaches have been proposed in the literature, based on different assumptions and techniques. Generally, depending on the used mathematical framework, they can be classified into two categories: computational fluid dynamics-based models and lumped-parameter ones. By focusing on gun systems, this paper offers a review of the main contributions in the field by mentioning their advantages and drawbacks. An insight into the limitations of the currently available modelling strategies is provided,as well as some considerations on the choice of one model over another. Lumped-parameter models, for example, are a good candidate for performing parametric analysis and optimisation processes of gun systems, given their minimum requirements of computer resources. Conversely, CFD-based models have a better capacity to address more sophisticated phenomena like pressure waves and turbulent flow effects. The performed review also reveals that too little attention has been given to small calibre guns since the majority of currently available models are conceived for medium and large calibre gun systems.Similarly, aspects like wear phenomena, bore deformations or projectile-barrel interactions still need to be adequately addressed and our suggestion is to dedicate more effort on it.

Investigation of Maxima Assumptions in Modelling Tropical Cyclone- Induced Hazards in the South China Sea

The present study aims to examine the suitability of two commonly used assumptions that simplify modelling metoceanconditions for designing offshore wind turbines in the South China Sea (SCS). The first assumption assumes thatjoint N-year extreme wind and wave events can be independently estimated and subsequently combined. The secondone assumes peak wind and waves can be modelled as occurring simultaneously during a tropical cyclone (TC) event.To better understand the potential TC activity, a set of 10000 years synthetic TC events are generated. The wind fieldmodel and the Mike 21 spectral wave model are employed to model the TC-induced hazards. Subsequently, theeffect of the assumptions is evaluated by analyzing the peak structural response of both monopile and semisubmersibleoffshore wind turbines during TC events. The results demonstrate that the examined assumptions are generally accurate.By assessing the implications of these assumptions, valuable insights are obtained, which can inform andimprove the modelling of TC-induced hazards in the SCS region.

Numerical analysis of high‑speed railway slab tracks using calibrated and validated 3D time‑domain modelling

Concrete slabs are widely used in modern railways to increase the inherent resilient quality of the tracks,provide safe and smooth rides,and reduce the maintenance frequency.In this paper,the elastic performance of a novel slab trackform for high-speed railways is investigated using three-dimensional finite element modelling in Abaqus.It is then compared to the performance of a ballasted track.First,slab and ballasted track models are developed to replicate the full-scale testing of track sections.Once the models are calibrated with the experimental results,the novel slab model is developed and compared against the calibrated slab track results.The slab and ballasted track models are then extended to create linear dynamic models,considering the track geodynamics,and simulating train passages at various speeds,for which the Ledsgard documented case was used to validate the models.Trains travelling at low and high speeds are analysed to investigate the track deflections and the wave propagation in the soil,considering the issues associated with critical speeds.Various train loading methods are discussed,and the most practical approach is retained and described.Moreover,correlations are made between the geotechnical parameters of modern high-speed rail and conventional standards.It is found that considering the same ground condition,the slab track deflections are considerably smaller than those of the ballasted track at high speeds,while they show similar behaviour at low speeds.

New model of linkage evolution for the transtensional fault systems in the Nanpu Sag of Bohai Bay Basin:Insight from seismic interpretation and analogue modelling

The evolution of faults within the same stress field is frequently influenced by numerous factors,involving the reactivation of pre-existing structures,stress transmission through ductile detachment layers,and the growth,interaction,as well as linkage of new fault segments.This study analyses a complex multi-phase oblique extension fault system in the Nanpu Sag(NPS)of the Bohai Bay Basin(BBB),China.High-resolution three-dimensional(3D)seismic data and analogue modelling indicate that the oblique extensional reactivation of pre-existing structures governs the sequential arrangement of fault segments in the caprock,and they dip synthetically to the reactivated fault at depth.During the NW-SE extension in the Eocene,the predominant movement of the pre-existing fault is strike-slip.Subsequently,during the N-S extension since the Oligocene,inclined at 20.to the pre-existing fault,forming splay fault segments and ultimately creating large en-echelon arcuate faults linked by relay ramps.Using fault throw-distance(T-D)and laser scanning,we reconstructed the fault evolution model of oblique extension reactivation in the presence of a ductile detachment basement.Our study illustrates that the arcuate faults can be categorized into linear master fault segments controlled by pre-existing structures,bending splay faults in the termination zone,and normal fault segments responding to the regional stress field.The interaction between faults occurs among normal faults and strike-slip faults,and the kinematic unification of the two fault systems is accomplished in the intersection zone.As the faults continue to evolve,the new fault segments tend to relinquish the control of pre-existing structures and concentrate more on the development of planar and continuous major faults.The ductile detachment layer significantly contributes to the uniform distribution of strain,resulting in narrow shear zones and discontinuous normal faults in its absence.

Assessment of sediment connectivity using modelling and field-based approaches in the Slavíč River catchment(MoravskoslezskéBeskydy Mts,Czech Republic)

Mountain streams act as conveyors of sediments within the river continuum,where the physical transport of sediments between river reaches through the catchment or between individual parts(e.g.,between hillslopes and channels)of the catchment is assumed.This study focused on sediment connectivity analysis in the SlavíčRiver catchment in the MoravskoslezskéBeskydy Mts in the eastern part of the Czech Republic.The connectivity index and connectivity index target modelling were combined with an analysis of anthropogenic interventions.Additionally,field mapping,grain size of bed sediments and stream power analysis were used to obtain information about connectivity in the catchment.Based on the analysis and obtained results,terrain topography is the current main driving factor affecting the connectivity of sediment movement in the SlavíčRiver catchment.However,the modelling provided valuable information about high sediment connectivity despite different recent land use conditions(highly forested area of the catchment)than those in historical times from the 16th to 19th centuries when the SlavíčRiver catchment was highly deforested and sediment connectivity was probably higher.The analysis of anthropogenic interventions,field mapping,grain size of bed sediments and stream power analysis revealed more deceleration of sediment movement through the catchment,decreased sediment connectivity with bed erosion,and gradual river channel process transformation in some reaches.Field mapping has identified various natural formations and human-induced changes impacting the longitudinal and lateral connectivity in the SlavíčRiver.For instance,embankments along 48%of the river's length,both on the right and left banks,significantly hinder lateral sediment supply to the channel.Stream power index analysis indicates increased energy levels in the flowing water in the river's upper reaches(up to 404.8 W m^(-2)).This high energy is also observed in certain downstream sections(up to 337.6 W m^(-2)),where it is influenced by human activities.These conditions lead to intensified erosion processes,playing a crucial role in sediment connectivity.Similar observations were described in recent studies that pointed out the long-term human interventions on many streams draining European mountains,where a decrease in sediment connectivity in these streams is linked with sediment deficits and the transformation of processes forming channels.

A data-adaptive network design for the regional gravity field modelling using spherical radial basis functions

A high-precision regional gravity field model is significant in various geodesy applications.In the field of modelling regional gravity fields,the spherical radial basis functions(SRBFs)approach has recently gained widespread attention,while the modelling precision is primarily influenced by the base function network.In this study,we propose a method for constructing a data-adaptive network of SRBFs using a modified Hierarchical Density-Based Spatial Clustering of Applications with Noise(HDBSCAN)algorithm,and the performance of the algorithm is verified by the observed gravity data in the Auvergne area.Furthermore,the turning point method is used to optimize the bandwidth of the basis function spectrum,which satisfies the demand for both high-precision gravity field and quasi-geoid modelling simultaneously.Numerical experimental results indicate that our algorithm has an accuracy of about 1.58 mGal in constructing the gravity field model and about 0.03 m in the regional quasi-geoid model.Compared to the existing methods,the number of SRBFs used for modelling has been reduced by 15.8%,and the time cost to determine the centre positions of SRBFs has been saved by 12.5%.Hence,the modified HDBSCAN algorithm presented here is a suitable design method for constructing the SRBF data adaptive network.

A framework for dynamic modelling of railway track switches considering the switch blades,actuators and control systems

The main contribution of this paper is the development and demonstration of a novel methodology that can be followed to develop a simulation twin of a railway track switch system to test the functionality in a digital environment.This is important because,globally,railway track switches are used to allow trains to change routes;they are a key part of all railway networks.However,because track switches are single points of failure and safety-critical,their inability to operate correctly can cause significant delays and concomitant costs.In order to better understand the dynamic behaviour of switches during operation,this paper has developed a full simulation twin of a complete track switch system.The approach fuses finite element for the rail bending and motion,with physics-based models of the electromechanical actuator system and the control system.Hence,it provides researchers and engineers the opportunity to explore and understand the design space around the dynamic operation of new switches and switch machines before they are built.This is useful for looking at the modification or monitoring of existing switches,and it becomes even more important when new switch concepts are being considered and evaluated.The simulation is capable of running in real time or faster meaning designs can be iterated and checked interactively.The paper describes the modelling approach,demonstrates the methodology by developing the system model for a novel“REPOINT”switch system,and evaluates the system level performance against the dynamic performance requirements for the switch.In the context of that case study,it is found that the proposed new actuation system as designed can meet(and exceed)the system performance requirements,and that the fault tolerance built into the actuation ensures continued operation after a single actuator failure.

Key Issues for Modelling, Operation, Management and Diagnosis of Lithium Batteries: Current States and Prospects

1 Introduction.With the continuous growth of the global population,the energy demand continues to increase.However,due to the dominance of fossil fuels in global energy and fossil fuels are non-renewable,it has led to the global energy crisis[1].Besides,the use of fossil fuels will generate a mass of air pollutants(e.g.,carbon dioxide,sulfur dioxide,etc.),which will cause serious environmental pollution,climate change[2],etc.To resolve the aforementioned issues,countries around the world have implemented a variety of measures hoping to fundamentally adjust the global energy structure and achieve sustainable development.Thereinto,“Paris Agreement”reached in 2015 under the framework of“United Nations Framework Convention on Climate Change”aims to control the increase in the average temperature of the globe to within 2°C below preindustrial levels,and thereafter to peak global greenhouse gas emissions as soon as possible,continuously decreasing thereafter[3].United Kingdom plans to reduce the average exhaust emissions of“new cars”to approximately 50–70 g/km by 20230,which is roughly half of what it is now[4].In addition,China proposed a plan at“United Nations General Assembly”in 2020 to peak carbon dioxide emissions by 2030 and strive to achieve carbon neutrality by 2060.It is a fact that the whole world is committed to changing the current energy structure,protecting the Earth’s ecology,and achieving global sustainable development[5].

Utility of fluid inclusion paleo-temperature in petroleum system modelling:A case study from western offshore,India

The paleo-temperature(Th)data from fluid inclusions are utilized for thermal history modelling using PetroMod software.Generally,bottom hole temperature(BHT)and vitrinite reflectance(Ro)measurements are widely used in petroleum system modelling(PSM)in the oil industry for calibration purposes.Th representing the minimum temperature of fluid entrapment estimated from fluid-inclusion study provides extra support to build the thermal models for PSM.Fluid inclusion parameters along with Rock-Eval pyrolysis analysis have been used to predict the maturity of oil in terms of API gravity as well as the maturity of source rocks respectively.Two exploratory wells RV-1(Mumbai Offshore Basin)and KK4C-A-1(Kerala-Konkan Offshore Basin),India were examined and the T_(h)from most of the fluid inclusions of wells RV-1 and KK4C-A-1 fell in the oil window range of 60-140℃suggesting thermal conditions favourable for oil generation in both of the wells.T_(h)of coeval aqueous inclusions along with the Hydrocarbon Fluid inclusions(HCFIs)was used to calibrate PSM.Vital parameters show that source rocks of well RV-1 are mature and that of well KK4C-A-1 are immature.Two sets of PSM are created in terms of generation and expulsion for the dry wells RV-1 and KK4C-A-1 and calibrated each well using fluid inclusion Th and BHT.From the fluid inclusion analysis method,it is evident that hydrocarbon generation happened in both wells and the paleo-temperature indicates that the formations of both wells were subjected to temperatures in the oil window range,even though it was designated as dry wells in the present scenario.The present study highlights the application of fluid inclusion paleo-temperature(Th)during calibration instead of commonly used methods.We could obtain desirable and accurate data output from PSM using T_(h) calibration.

Numerical Modelling of Coupled Heat and Mass Transfer in Porous Materials: Application to Cinder Block Bricks

In this work, we present numerical modelling of coupled heat and mass transfer within porous materials. Our study focuses on cinder block bricks generally used in building construction. The material is assumed to be placed in air. Moisture content and temperature have been chosen as the main transfer drivers and the equations governing these transfer drivers are based on the Luikov model. These equations are solved by an implicit finite difference scheme. A Fortran code associated with the Thomas algorithm was used to solve the equations. The results show that heat and mass transfer depend on the temperature of the air in contact with the material. As this air temperature rises, the temperature within the material increases, and more rapidly at the material surface. Also, thermal conductivity plays a very important role in the thermal conduction of building materials and influences heat and mass transfer in these materials. Materials with higher thermal conductivity diffuse more heat.

Data for critical infrastructure network modelling of natural hazard impacts:Needs and influence on model characteristics

Natural hazards impact interdependent infrastructure networks that keep modern society functional.While a va-riety of modelling approaches are available to represent critical infrastructure networks(CINs)on different scales and analyse the impacts of natural hazards,a recurring challenge for all modelling approaches is the availability and accessibility of sufficiently high-quality input and validation data.The resulting data gaps often require mod-ellers to assume specific technical parameters,functional relationships,and system behaviours.In other cases,expert knowledge from one sector is extrapolated to other sectoral structures or even cross-sectorally applied to fill data gaps.The uncertainties introduced by these assumptions and extrapolations and their influence on the quality of modelling outcomes are often poorly understood and difficult to capture,thereby eroding the reliability of these models to guide resilience enhancements.Additionally,ways of overcoming the data avail-ability challenges in CIN modelling,with respect to each modelling purpose,remain an open question.To address these challenges,a generic modelling workflow is derived from existing modelling approaches to examine model definition and validations,as well as the six CIN modelling stages,including mapping of infrastructure assets,quantification of dependencies,assessment of natural hazard impacts,response&recovery,quantification of CI services,and adaptation measures.The data requirements of each stage were systematically defined,and the literature on potential sources was reviewed to enhance data collection and raise awareness of potential pitfalls.The application of the derived workflow funnels into a framework to assess data availability challenges.This is shown through three case studies,taking into account their different modelling purposes:hazard hotspot assess-ments,hazard risk management,and sectoral adaptation.Based on the three model purpose types provided,a framework is suggested to explore the implications of data scarcity for certain data types,as well as their reasons and consequences for CIN model reliability.Finally,a discussion on overcoming the challenges of data scarcity is presented.

Contribution to the Full 3D Finite Element Modelling of a Hybrid Stepping Motor with and without Current in the Coils

The paper presents our contribution to the full 3D finite element modelling of a hybrid stepping motor using COMSOL Multiphysics software. This type of four-phase motor has a permanent magnet interposed between the two identical and coaxial half stators. The calculation of the field with or without current in the windings (respectively with or without permanent magnet) is done using a mixed formulation with strong coupling. In addition, the local high saturation of the ferromagnetic material and the radial and axial components of the magnetic flux are taken into account. The results obtained make it possible to clearly observe, as a function of the intensity of the bus current or the remanent induction, the saturation zones, the lines, the orientations and the magnetic flux densities. 3D finite element modelling provide more accurate numerical data on the magnetic field through multiphysics analysis. This analysis considers the actual operating conditions and leads to the design of an optimized machine structure, with or without current in the windings and/or permanent magnet.

Modelling of blast furnace ironmaking process based on digital twin

The blast furnace ironmaking process is a crucial step in the steel industry.Effectively modelling the blast furnaces is significant in ensuring smooth operation and accelerating the digitalisation transformation of blast furnaces.The authors focus on the mechanism modelling of the blast furnace operation process,using a digital twin model development platform to simulate the main reaction processes inside the blast furnaces.Under on-site production conditions,Si,Mn and Ti contents of molten iron and the corresponding indicators for model accuracy are calculated.For Si,Mn,and Ti content,the Root Mean Square Error values are 0.0678,0.0108 and 0.0093 for dataset 1,while 0.0933,0.0120 and 0.0111 for dataset 2,respectively,indicating that the model has a small simulation error and high accuracy.By comparing the simulation results with accurate laboratory results,the model is validated to have satisfactory simulation reliability and compensates for the existing shortcomings in the blast furnace mechanism modelling field.

Application of 3D Static Modelling in Reservoir Characterization:A Case Study from the Qishn Formation in Sharyoof Oil Field,Masila Basin,Yemen

Three-dimensional(3 D)static modelling techniques are applied to the characterization of the Qishn Formation(Fm.)in the Sharyoof oil field locating within the Masila basin,southeastern Yemen.The present study was initiated by the seismic structural interpretation,followed by building a 3 D structural framework,in addition to analysing well log data and from these,3 D facies and petrophysical models are constructed.In the Sharyoof oil field,the Qishn Fm.exhibits depth values within the range of 400-780 m below sea level,with a general increase towards the SSE.A set of high dip angle normal faults with a general ENE-WSW trend dissect the rocks.The strata are also folded as a main anticline with an axis that is parallel to the fault trend,formed as a result of basement uplift.According to the facies models,the Qishn Fm.comprises 43.83% limestone,21.53% shale,21.26% sandstone,13.21% siltstone and 0.17% dolomite.The Qishn Carbonates Member has low porosity values making it a potential seal for the underlying reservoirs whereas the Upper Qishn Clastics SI A and C have good reservoir quality and SIB has fair reservoir quality.The Upper Qishn Clastics S2 and S3 also have fair reservoir quality,while the Lower Qishn Clastics zone has good reservoir quality.The water saturation decreases towards the west and east and increases towards north and south.The total original oil in-place(OOIP)of the Upper Qishn clastics is 106 million STB within the SI A,SIC and S2 zones.Drilling of development wells is recommended in the eastern study area,where good trapping configuration is exhibited in addition to the presence of a potential seal(Upper Qishn Carbonates Member)and reservoir(Qishn Clastics Member)with high porosity and low water saturation.

Implication of community-level ecophysiological parameterization to modelling ecosystem productivity:a case study across nine contrasting forest sites in eastern China

Parameterization is a critical step in modelling ecosystem dynamics.However,assigning parameter values can be a technical challenge for structurally complex natural plant communities;uncertainties in model simulations often arise from inappropriate model parameterization.Here we compared five methods for defining community-level specific leaf area(SLA)and leaf C:N across nine contrasting forest sites along the North-South Transect of Eastern China,including biomass-weighted average for the entire plant community(AP_BW)and four simplified selective sampling(biomass-weighted average over five dominant tree species[5DT_BW],basal area weighted average over five dominant tree species[5DT_AW],biomass-weighted average over all tree species[AT_BW]and basal area weighted average over all tree species[AT_AW]).We found that the default values for SLA and leaf C:N embedded in the Biome-BGC v4.2 were higher than the five computational methods produced across the nine sites,with deviations ranging from 28.0 to 73.3%.In addition,there were only slight deviations(<10%)between the whole plant community sampling(AP_BW)predicted NPP and the four simplified selective sampling methods,and no significant difference between the predictions of AT_BW and AP_BW except the Shennongjia site.The findings in this study highlights the critical importance of computational strategies for community-level parameterization in ecosystem process modelling,and will support the choice of parameterization methods.

Effect of dry-wet cycles on dynamic properties and microstructures of sandstone:Experiments and modelling

Underground pumped storage power plant(UPSP)is an innovative concept for space recycling of abandoned mines.Its realization requires better understanding of the dynamic performance and durability of reservoir rock.This paper conducted ultrasonic detection,split Hopkinson pressure bar(SHPB)impact,mercury intrusion porosimetry(MIP),and backscatter electron observation(BSE)tests to investigate the dynamical behaviour and microstructure of sandstone with cyclical dry-wet damage.A coupling FEM-DEM model was constructed for reappearing mesoscopic structure damage.The results show that dry-wet cycles decrease the dynamic compressive strength(DCS)with a maximum reduction of 39.40%,the elastic limit strength is reduced from 41.75 to 25.62 MPa.The sieved fragments obtain the highest crack growth rate during the 23rd dry-wet cycle with a predictable life of 25 cycles for each rock particle.The pore fractal features of the macropores and micro-meso pores show great differences between the early and late cycles,which verifies the computational statistics analysis of particle deterioration.The numerical results show that the failure patterns are governed by the strain in pre-peak stage and the shear cracks are dominant.The dry-wet cycles reduce the energy transfer efficiency and lead to the discretization of force chain and crack fields.

Experiment and constitutive modelling of creep deformation in the frozen silt-concrete interface

To ensure the long-term safety and stability of bridge pile foundations in permafrost regions,it is necessary to investigate the rheological effects on the pile tip and pile side bearing capacities.The creep characteristics of the pile-frozen soil interface are critical for determining the long-term stability of permafrost pile foundations.This study utilized a self-developed large stress-controlled shear apparatus to investigate the shear creep characteristics of the frozen silt-concrete interface,and examined the influence of freezing temperatures(−1,−2,and−5°C),contact surface roughness(0,0.60,0.75,and 1.15 mm),normal stress(50,100,and 150 kPa),and shear stress on the creep characteristics of the contact surface.By incorporating the contact surface’s creep behavior and development trends,we established a creep constitutive model for the frozen silt-concrete interface based on the Nishihara model,introducing nonlinear elements and a damage factor.The results revealed significant creep effects on the frozen silt-concrete interface under constant load,with creep displacement at approximately 2-15 times the instantaneous displacement and a failure creep displacement ranging from 6 to 8 mm.Under different experimental conditions,the creep characteristics of the frozen silt-concrete interface varied.A larger roughness,lower freezing temperatures,and higher normal stresses resulted in a longer sample attenuation creep time,a lower steady-state creep rate,higher long-term creep strength,and stronger creep stability.Building upon the Nishihara model,we considered the influence of shear stress and time on the viscoelastic viscosity coefficient and introduced a damage factor to the viscoplasticity.The improved model effectively described the entire creep process of the frozen silt-concrete interface.The results provide theoretical support for the interaction between pile and soil in permafrost regions.