A Transient-Pressure-Based Numerical Approach for Interlayer Identification in Sand Reservoirs

Almost all sandstone reservoirs contain interlayers. The identification and characterization of these interlayers iscritical for minimizing the uncertainty associated with oilfield development and improving oil and gas recovery.Identifying interlayers outside wells using identification methods based on logging data and machine learning isdifficult and seismic-based identification techniques are expensive. Herein, a numerical model based on seepageand well-testing theories is introduced to identify interlayers using transient pressure data. The proposed modelrelies on the open-source MATLAB Reservoir Simulation Toolbox. The effects of the interlayer thickness, position,and width on the pressure response are thoroughly investigated. A procedure for inverting interlayer parametersin the reservoir using the bottom-hole pressure is also proposed. This method uses only transient pressuredata during well testing and can effectively identify the interlayer distribution near the wellbore at an extremelylow cost. The reliability of the model is verified using effective oilfield examples.

Improving the Short-Range Precipitation Forecast of Numerical Weather Prediction through a Deep Learning-Based Mask Approach

Due to various technical issues,existing numerical weather prediction(NWP)models often perform poorly at forecasting rainfall in the first several hours.To correct the bias of an NWP model and improve the accuracy of short-range precipitation forecasting,we propose a deep learning-based approach called UNet Mask,which combines NWP forecasts with the output of a convolutional neural network called UNet.The UNet Mask involves training the UNet on historical data from the NWP model and gridded rainfall observations for 6-hour precipitation forecasting.The overlap of the UNet output and the NWP forecasts at the same rainfall threshold yields a mask.The UNet Mask blends the UNet output and the NWP forecasts by taking the maximum between them and passing through the mask,which provides the corrected 6-hour rainfall forecasts.We evaluated UNet Mask on a test set and in real-time verification.The results showed that UNet Mask outperforms the NWP model in 6-hour precipitation prediction by reducing the FAR and improving CSI scores.Sensitivity tests also showed that different small rainfall thresholds applied to the UNet and the NWP model have different effects on UNet Mask's forecast performance.This study shows that UNet Mask is a promising approach for improving rainfall forecasting of NWP models.

Induced pluripotent stem cell-related approaches to generate dopaminergic neurons for Parkinson's disease

The progressive loss of dopaminergic neurons in affected patient brains is one of the pathological features of Parkinson's disease,the second most common human neurodegenerative disease.Although the detailed pathogenesis accounting for dopaminergic neuron degeneration in Parkinson's disease is still unclear,the advancement of stem cell approaches has shown promise for Parkinson's disease research and therapy.The induced pluripotent stem cells have been commonly used to generate dopaminergic neurons,which has provided valuable insights to improve our understanding of Parkinson's disease pathogenesis and contributed to anti-Parkinson's disease therapies.The current review discusses the practical approaches and potential applications of induced pluripotent stem cell techniques for generating and differentiating dopaminergic neurons from induced pluripotent stem cells.The benefits of induced pluripotent stem cell-based research are highlighted.Various dopaminergic neuron differentiation protocols from induced pluripotent stem cells are compared.The emerging three-dimension-based brain organoid models compared with conventional two-dimensional cell culture are evaluated.Finally,limitations,challenges,and future directions of induced pluripotent stem cell–based approaches are analyzed and proposed,which will be significant to the future application of induced pluripotent stem cell-related techniques for Parkinson's disease.

Enhancing perianal disease management with integrated physical and psychological approaches

This article provides a comprehensive analysis of the study by Hou et al,focusing on the complex interplay between psychological and physical factors in the postoperative recovery(POR)of patients with perianal diseases.The study sheds light on how illness perception,anxiety,and depression significantly influence recovery outcomes.Hou et al developed a predictive model that demonstrated high accuracy in identifying patients at risk of poor recovery.The article explores the critical role of pre-operative psychological assessment,highlighting the need for mental health support and personalized recovery plans in enhancing POR quality.A multidisciplinary approach,integrating mental health professionals with surgeons,anesthesiologists,and other specialists,is emphasized to ensure comprehensive care for patients.The study’s findings serve as a call to integrate psychological care into surgical practice to optimize outcomes for patients with perianal diseases.

Chyle leak following root of mesentery dissection in pancreaticoduodenectomy with inferior infracolic superior mesenteric artery first approach

BACKGROUND The root of mesentery dissection is one of the critical maneuvers,especially in borderline resectable pancreatic head cancer.Intra-abdominal chyle leak(CL)including chylous ascites may ensue in up to 10%of patients after pancreatic resections.Globally recognized superior mesenteric artery(SMA)first approaches are invariably performed.The mesenteric dissection through the inferior infracolic approach has been discussed in this study emphasizing its post-operative impact on CL which is the cornerstone of this study.AIM To assess incidence,risk factors,clinical impact of CL following root of mesentery dissection,and the different treatment modalities.METHODS This is a retrospective study incorporating the patients who underwent dissection of the root of mesentery with inferior infracolic SMA first approach pancreat-oduodenectomy for the ventral body and uncinate mass of pancreas in the Department of Gastrointestinal and General Surgery of Kathmandu Medical College and Teaching Hospital from January 1,2021 to February 28,2024.Intraop-erative findings and postoperative outcomes were analyzed.RESULTS In three years,ten patients underwent root of mesentery dissection with inferior infracolic SMA first approach pancreatoduodenectomy.The mean age was 67.6 years with a male-to-female ratio of 4:5.CL was seen in four patients.With virtue of CL,Clavien-Dindo grade Ⅱ or higher morbidity was observed in four patients.Two patients had a hospital stay of more than 20 days with the former having a delayed gastric emptying and the latter with long-term total parenteral nutrition requirement.The mean operative time was 330 minutes.Curative resection was achieved in 100%of the patients.The mean duration of the intensive care unit and hospital stay were 2.55±1.45 days and 15.7±5.32 days,respectively.CONCLUSION Root of mesentery dissection with lymphadenectomy and vascular resection correlated with occurrence of CL.After complete curative resection,these were managed with total parenteral nutrition without adversely impacting outcome.

An evaluation of numerical approaches for S-wave component simulation in rock blasting

The shear wave(S-wave) component of the total blast vibration always plays an important role in damage to rock or adjacent structures.Numerical approach has been considered as an economical and effective tool in predicting blast vibration.However,S-wave has not yet attracted enough attention in previous numerical simulations.In this paper,three typical numerical models,i.e.the continuum-based elastic model,the continuum-based damage model,and the coupled smooth particle hydrodynamics(SPH)-finite element method(FEM) model,were first introduced and developed to simulate the blasting of a single cylindrical charge.Then,the numerical results from different models were evaluated based on a review on the generation mechanisms of S-wave during blasting.Finally,some suggestions on the selection of numerical approaches for simulating generation of the blast-induced S-wave were put forward.Results indicate that different numerical models produce different results of S-wave.The coupled numerical model was the best,for its outstanding capacity in producing S-wave component.It is suggested that the model that can describe the cracking,sliding or heaving of rock mass,and the movement of fragments near the borehole should be selected preferentially,and priority should be given to the material constitutive law that could record the nonlinear mechanical behavior of rock mass near the borehole.

A Numerical Study on Effects of Land-Surface Heterogeneity from' Combined Approach' on Atmospheric ProcessPart II: Coupling-Model Simulations

Two land surface schemes, one the standard Biosphere / Atmosphere Transfer Scheme Version ie (BOZ) and the other B1Z based on B0Z and heterogeneously-treated by' combined approach' , were co 'pled to the meso-scale model MM4, respectively. Through the calculations of equations from the companion paper, parameters representing land surface heterogeneity and suitable for the coupling models were found out. Three cases were simulated for heavy rainfalls during 36 hours, and the sensitivity of short-term weather modeling to the land surface heterogeneity was tested. Through the analysis of the simulations of the three heavy rainfalls, it was demonstrated that BIZ, compared with BOZ, could more realistically reflect the features of the land surface heterogeneity, therefore could more realistically reproduce the circulation and precipitation amount in the heavy rainfall processes of the three cases. This shows that even short-term weather is sensitive to the land surface heterogeneity, which is more obvious with time passing, and whose influence is more pronounced in the lower layer and gradually extends to the middle and upper layer. Through the analysis of these simulations with BlZ, it is suggested that the bulk effect of smaller-scale fluxes (i.e., the momentum, water vapor and sensible heat fluxes) near the s ig nificantly-heterogeneous land surface is to change the larger-scale (i.e., meso-scale) circulation, and then to influence the development of the low-level jets and precipitation. And also, the complexity of the land-atmosphere interaction was shown in these simulations.

A Numerical Study on Effects of Land-Surface Heterogeneity from "Combined Approach" on Atmospheric Process Part I: Principle and Method

A method based on Giorgi (1997a, 1997b) and referred to as ’ combined approach’, which is a combi-nation of mosaic approach and analytical-statistical-dynamical approach, is proposed. Compared with those of other approaches, the main advantage of the combined approach is that it not only can represent both interpatch and intrapatch variability, but also cost less computational time when the land surface heterogeneity is considered. Because the independent variable of probability density function (PDF) is ex-tended to the single valued function of basic meteorological characteristic quantities, which is much more universal, the analytical expressions of the characteristic quantities (e.g., drag coefficient, snow coverage, leaf surface aerodynamical resistance) affected by roughness length are derived , when the roughness length(and / or the zero plane displacement) heterogeneity has been mainly taken into account with the approach. On the basis of the rule which the PDF parameters should follow, we choose a function y of the roughness length z 0 as the PDF independent variable, and set different values of the two parameters width ratio αn and height ratio γ of PDF (here a linear, symmetric PDF is applied) for sensitivity experiments, from which some conclusions can be drawn, e.g., relevant characteristic terms show different sensitivities to the heterogeneous characteristic (i.e., roughness length), which suggests that we should consider the heterogeneities of the more sensitive terms in our model instead of the heterogeneities of the rest, and which also implies that when the land surface scheme is coupled into the global or regional atmospheric model, sensitivity tests against the distribution of the heterogeneous characteristic are very necessary; when the parameter αn is close to zero, little heterogeneity is represented, and αn differs with cases, which have an upper limit of about 0.6; in the reasonable range of αn, a peak-like distribution of roughness length can be depicted by a small value of γ, etc.. Key words Representation of land surface heterogeneity - “ Combined approach” - Numerical experiment This work was supported by the National Sciences Foundation of China, Grant No.49875005 and the State Key Project (973) G19990434 (03).

Assimilated Tidal Results of Tide Gauge and TOPEX/POSEIDON Data over the China Seas Using a Variational Adjoint Approach with a Nonlinear Numerical Model

In order to obtain an accurate tide description in the China Seas, the 2-dimensional nonlinear numerical Princeton Ocean Model （POM） is employed to incorporate in situ tidal measurements both from tide gauges and TOPEX/POSEIDON （T/P） derived datasets by means of the variational adjoint approach in such a way that unknown internal model parameters, bottom topography, friction coefficients and open boundary conditions, for example, are adjusted during the process. The numerical model is used as a forward model. After the along-track T/P data are processed, two classical methods, i.e. harmonic and response analysis, are implemented to estimate the tide from such datasets with a domain covering the model area extending from 0° to 41°N in latitude and from 99°E to 142°E in longitude. And the results of these two methods are compared and interpreted. The numerical simulation is performed for 16 major constituents. In the data assimilation experiments, three types of unknown parameters （water depth, bottom friction and tidal open boundary conditions in the model equations） are chosen as control variables. Among the various types of data assimilation experiments, the calibration of water depth brings the most promising results. By comparing the results with selected tide gauge data, the average absolute errors are decreased from 7.9 cm to 6.8 cm for amplitude and from 13.0° to 9.0° for phase with respect to the semidiurnal tide M2 constituent, which is the largest tidal constituent in the model area. After the data assimilation experiment is performed, the comparison between model results and tide gauge observation for water levels shows that the RMS errors decrease by 9 cm for a total of 14 stations, mostly selected along the coast of China's Mainland, when a one-month period is considered, and the correlation coefficients improve for most tidal stations among these stations.

Bending of Euler-Bernoulli nanobeams based on the strain-driven and stress-driven nonlocal integral models: a numerical approach

Eringen's nonlocal elasticity theory is extensively employed for the analysis of nanostructures because it is able to capture nanoscale effects.Previous studies have revealed that using the differential form of the strain-driven version of this theory leads to paradoxical results in some cases,such as bending analysis of cantilevers,and recourse must be made to the integral version.In this article,a novel numerical approach is developed for the bending analysis of Euler-Bernoulli nanobeams in the context of strain-and stress-driven integral nonlocal models.This numerical approach is proposed for the direct solution to bypass the difficulties related to converting the integral governing equation into a differential equation.First,the governing equation is derived based on both strain-driven and stress-driven nonlocal models by means of the minimum total potential energy.Also,in each case,the governing equation is obtained in both strong and weak forms.To solve numerically the derived equations,matrix differential and integral operators are constructed based upon the finite difference technique and trapezoidal integration rule.It is shown that the proposed numerical approach can be efficiently applied to the strain-driven nonlocal model with the aim of resolving the mentioned paradoxes.Also,it is able to solve the problem based on the strain-driven model without inconsistencies of the application of this model that are reported in the literature.

Air Gap Prediction for Floating Bodies Using a 3D Numerical Wave Tank Approach

Computations for air gap response of a semisubmersible platform based on a 3D numerical wave tank approach are presented.The developed method is in time domain and can consider nonlinearities associated with incident wave and hydrostatic forces exactly in determining the body response, but the interaction hydrodynamics of radiation and diffraction are based on simplified linearization assumptions. The incident wave can be defined by any suitable wave theory and here defined by a fully nonlinear numerical wave model. After verifying the present computations results in its degenerated linearized version against the usual linear 3D Green function–based frequency-domain results for air gap predictions, systematic comparative studies are undertaken between linear and the approximate nonlinear solutions. It is found that nonlinear computations can yield considerably conservative predictions as compared to fully linear calculations, amounting to a difference of up to 30%–40% in the minimum air gap in steep ambient incident waves at high and moderate frequencies.

An approach to estimating and extrapolating model error based on inverse problem methods:towards accurate numerical weather prediction

Model error is one of the key factors restricting the accuracy of numerical weather prediction （NWP）. Considering the continuous evolution of the atmosphere, the observed data （ignoring the measurement error） can be viewed as a series of solutions of an accurate model governing the actual atmosphere. Model error is represented as an unknown term in the accurate model, thus NWP can be considered as an inverse problem to uncover the unknown error term. The inverse problem models can absorb long periods of observed data to generate model error correction procedures. They thus resolve the deficiency and faultiness of the NWP schemes employing only the initial-time data. In this study we construct two inverse problem models to estimate and extrapolate the time-varying and spatial-varying model errors in both the historical and forecast periods by using recent observations and analogue phenomena of the atmosphere. Numerical experiment on Burgers＇ equation has illustrated the substantial forecast improvement using inverse problem algorithms. The proposed inverse problem methods of suppressing NWP errors will be useful in future high accuracy applications of NWP.

Science Letters:On numerical calculation in symplectic approach for elasticity problems

The symplectic approach proposed and developed by Zhong et al. in 1990s for elasticity problems is a rational analytical method, in which ample experience is not needed as in the conventional semi-inverse method. In the symplectic space, elasticity problems can be solved using the method of separation of variables along with the eigenfunction expansion technique, as in traditional Fourier analysis. The eigensolutions include those corresponding to zero and nonzero eigenvalues. The latter group can be further divided into α-and β-sets. This paper reformulates the form of β-set eigensolutions to achieve the stability of numerical calculation, which is very important to obtain accurate results within the symplectic frame. An example is finally given and numerical results are compared and discussed.

A Simplified Numerical Approach for the Prediction of Rainfall-Induced Retrogressive Landslides

Retrogressive landslides are common geological phenomena in mountainous areas and on onshore and offshore slopes. The impact of retrogressive landslides is different from that of other landslide types due to the phenomenon of retrogression. The hazards caused by retrogressive landslides may be increased because retrogressive landslides usually affect housing, facilities, and infrastructure located far from the original slopes. Additionally, substantial geomorphic evidence shows that the abundant supply of loose sediment in the source area of a debris flow is usually provided by retrogressive landslides that are triggered by the undercutting of water. Moreover, according to historic case studies, some large landslides are the evolution result of retrogressive landslides. Hence the ability to understand and predict the evolution of retrogressive landslides is crucial for the purpose of hazard mitigation. This paper discusses the phenomenon of a retrogressive landslide by using a model experiment and suggests a reasonably simplified numerical approach for the prediction of rainfall-induced retrogressive landslides. The simplified numerical approach, which combines the finite element method for seepage analysis, the shear strength reduction finite element method, and the analysis criterion for the retrogression and accumulation effect, is presented and used to predict the characteristics of a retrogressive landslide. The results show that this numerical approach is capable of reasonably predicting the characteristics of retrogressive landslides under rainfall infiltration, particularly the magnitude of each landslide, the position of the slip surface, and the development processes of the retrogressive landslide. Therefore, this approach is expected to be a practical method for the mitigation of damage caused by rainfall-induced retrogressive landslides.

A numerical modelling approach to assess the behaviour of underground cavern subjected to blast loads

The paper gives an insight into the behaviour of large underground caverns which are subjected to blast loads. Caverns are generally constructed in hard rock formation which compels us to use blasting methods for the excavation works. Comparative study was done between models with intact rock mass and discontinuities to assess the stability of cavern as a result of blast loads. Numerical modelling was performed with 3 dimensional distinct element code(3 DEC) to analyse the performance of cavern walls in terms of displacement and to compute peak particle velocities(PPV) both around the cavern periphery and at surface of models. Results showed that the velocity wave with higher frequency exhibited large displacements around the periphery of cavern. Computation of PPV showed that model with horizontal joint sets showed lower PPV in comparison to model with intact rock mass. PPV values were also analysed on the surface for model consisting vertical joints spaced at 4 m intervals. Comparative study of PPV on surface vertically above the blast location between models with horizontal joints spaced at 4 m and vertical joints at 4 m intervals were conducted. Results depicted higher magnitudes of PPV for model with vertical joints in comparison to model with horizontal joints.

Overhanging rock slope by design:An integrated approach using rock mass strength characterisation,large-scale numerical modelling and limit equilibrium methods

Overhanging rock slopes(steeper than 90°) are typically avoided in rock engineering design, particularly where the scale of the slope exceeds the scale of fracturing present in the rock mass. This paper highlights an integrated approach of designing overhanging rock slopes where the relative dimensions of the slope exceed the scale of fracturing and the rock mass failure needs to be considered rather than kinematic release of individual blocks. The key to the method is a simplified limit equilibrium(LE) tool that was used for the support design and analysis of a multi-faceted overhanging rock slope. The overhanging slopes required complex geometries with constantly changing orientations. The overhanging rock varied in height from 30 m to 66 m. Geomechanical modelling combined with discrete fracture network(DFN)representation of the rock mass was used to validate the rock mass strength assumptions and the failure mechanism assumed in the LE model. The advantage of the simplified LE method is that buttress and support design iterations(along with sensitivity analysis of design parameters) can be completed for various cross-sections along the proposed overhanging rock sections in an efficient manner, compared to the more time-intensive, sophisticated methods that were used for the initial validation. The method described presents the development of this design tool and assumptions made for a specific overhanging rock slope design. Other locations will have different geological conditions that can control the potential behaviour of rock slopes, however, the approach presented can be applied as a general guiding design principle for overhanging rock cut slope.

A Numerical Approach for Multicomponent Vapor Solid Equilibrium Calculations in Gas Hydrate Formation

A new numerical approach has been developed for vapor solid equilibrium calculations and for predicting vapor solid equilibrium constant and composition of vapor and solid phases in gas hydrate formation. Equation of state methods generally do a good job of determining vapor phase properties, but for solid phase it is much more difficult and inaccurate. This proposed new model calculates vapor solid equilibrium constant and vapor and solid phase composition as a function of temperature and partial pressure. The results of this proposed numerical approach, for vapor solid equilibrium, have a good agreement with the available reported data. This new numerical model also has an advantage to tune coefficients, to cover different sets of experimental data accurately.

Depth factors for undrained bearing capacity of circular footing by numerical approach

The undrained vertical bearing capacity of embedded foundation has been extensively studied using analytical and numerical methods.Through comparing the results of a circular embedded foundation in the literature,a significant difference between the bearing capacity factors and depth factors is observed.Based on the previous research findings,numerical computations using FLAC code are carried out in this study to evaluate the undrained bearing capacity of circular foundations with embedment ratios up to five for different base and side foundation roughness conditions.Unlike the foundation base,the roughness of the foundation side has a significant effect on the bearing capacity.The comparison of the present results with numerical studies available in the literature shows that the discrepancy is related to the procedures used to simulate the foundation side interface conditions and to the estimation of the bearing capacity.

A new approach to numerical simulation of source development of earthquake

By considering the heterogeneity of geomechanical materials, the source development of earthquake under compression boundary conditions is studied with a newly developed numerical method, Rock Failure Process Analysis code (RFPA2D). The process of fault forming and associated micro seismicities in a rectangle area with a inclusion but without any clear structural features of original fault is modeled. The modeling demonstrates the whole process of source development of earthquake from deformation, micro failure to collapse and the behavior of temporal spatial distribution of micro seismicities. The stress, strain and the temporal spatial distribution of micro seismicities vividly depict the phenomena of localization, temporal transitions, dilatation or rise, elastic rebound and conjugate (X type) deformation zone.

Modelling of blast-induced damage in tunnels using a hybrid finite-discrete numerical approach

This paper presents the application of a hybrid finite-discrete element method to study blast-induceddamage in circular tunnels. An extensive database of field tests of underground explosions above tunnelsis used for calibrating and validating the proposed numerical method; the numerical results areshown to be in good agreement with published data for large-scale physical experiments. The method isthen used to investigate the influence of rock strength properties on tunnel durability to withstand blastloads. The presented analysis considers blast damage in tunnels excavated through relatively weak（sandstone） and strong （granite） rock materials. It was found that higher rock strength will increase thetunnel resistance to the load on one hand, but decrease attenuation on the other hand. Thus, undercertain conditions, results for weak and strong rock masses are similar. 2014 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Production and hosting byElsevier B.V. All rights reserved.