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 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.

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.

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.

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.

Modelling analysis embodies drastic transition among global potential natural vegetations in face of changing climate

Potential natural vegetation(PNV)is a valuable reference for ecosystem renovation and has garnered increasing attention worldwide.However,there is limited knowledge on the spatio-temporal distributions,transitional processes,and underlying mechanisms of global natural vegetation,particularly in the case of ongoing climate warming.In this study,we visualize the spatio-temporal pattern and inter-transition procedure of global PNV,analyse the shifting distances and directions of global PNV under the influence of climatic disturbance,and explore the mechanisms of global PNV in response to temperature and precipitation fluctuations.To achieve this,we utilize meteorological data,mainly temperature and precipitation,from six phases:the Last Inter-Glacial(LIG),the Last Glacial Maximum(LGM),the Mid Holocene(MH),the Present Day(PD),2030(20212040)and 2090(2081–2100),and employ a widely-accepted comprehensive and sequential classification sy–stem(CSCS)for global PNV classification.We find that the spatial patterns of five PNV groups(forest,shrubland,savanna,grassland and tundra)generally align with their respective ecotopes,although their distributions have shifted due to fluctuating temperature and precipitation.Notably,we observe an unexpected transition between tundra and savanna despite their geographical distance.The shifts in distance and direction of five PNV groups are mainly driven by temperature and precipitation,although there is heterogeneity among these shifts for each group.Indeed,the heterogeneity observed among different global PNV groups suggests that they may possess varying capacities to adjust to and withstand the impacts of changing climate.The spatio-temporal distributions,mutual transitions and shift tendencies of global PNV and its underlying mechanism in face of changing climate,as revealed in this study,can significantly contribute to the development of strategies for mitigating warming and promoting re-vegetation in degraded regions worldwide.

Prospect Theory Based Individual Irrationality Modelling and Behavior Inducement in Pandemic Control

Understanding and modeling individuals’behaviors during epidemics is crucial for effective epidemic control.However,existing research ignores the impact of users’irrationality on decision-making in the epidemic.Meanwhile,existing disease control methods often assume users’full compliance with measures like mandatory isolation,which does not align with the actual situation.To address these issues,this paper proposes a prospect theorybased framework to model users’decision-making process in epidemics and analyzes how irrationality affects individuals’behaviors and epidemic dynamics.According to the analysis results,irrationality tends to prompt conservative behaviors when the infection risk is low but encourages risk-seeking behaviors when the risk is high.Then,this paper proposes a behavior inducement algorithm to guide individuals’behaviors and control the spread of disease.Simulations and real user tests validate our analysis,and simulation results show that the proposed behavior inducement algorithm can effectively guide individuals’behavior.

Modelling the viscoplastic behaviour of Callovo-Oxfordian claystone with consideration of damage effect

In order to evaluate the performance of deep geological disposal of radioactive waste,an underground research laboratory(URL)was constructed by Andra in the Callovo-Oxfordian(COx)claystone formation at the Meuse/Haute-Marne(MHM).The construction of URL induced the excavation damage of host formations,and the ventilation in the galleries desaturated the host formation close to the gallery wall.Moreover,it is expected that the mechanical behaviour of COx claystone is time-dependent.This study presents a constitutive model developed to describe the viscoplastic behaviour of unsaturated and damaged COx claystone.In this model,the unsaturation effect is considered by adopting the Bishop effective stress and the van Genuchten(VG)water retention model.In terms of the viscoplastic behaviour,the nonstationary flow surface(NSFS)theory for unsaturated soils is used with consideration of the coupled effects of strain rate and suction on the yield stress.A progressive hardening law is adopted.Meanwhile,a non-associated flow rule is used,which is similar to that in Barcelona basic model(BBM).In addition,to describe the damage effect induced by suction change and viscoplastic loading,a damage function is defined based on the crack volume proportion.This damage function contains two variables:unsaturated effective stress and viscoplastic volumetric strain,with the related parameters determined based on the mercury intrusion porosimetry(MIP)tests.For the model validation,different tests on COx claystone under different loading paths are simulated.Comparisons between experimental and simulated results indicated that the present model is able to well describe the viscoplastic behaviour of damaged COx claystone,including swelling/shrinkage,triaxial extension and compression,and triaxial creep.

Discrete Element Modelling of Damage Evolution of Concrete Considering Meso-Structure of ITZ

The mechanical properties of interfacial transition zones(ITZs)have traditionally been simplified by reducing the stiffness of cement in previous simulation methods.A novel approach based on the discrete element method(DEM)has been developed for modeling concrete.This new approach efficiently simulates the meso-structure of ITZs,accurately capturing their heterogeneous properties.Validation against established uniaxial compression experiments confirms the precision of thismodel.The proposedmodel canmodel the process of damage evolution containing cracks initiation,propagation and penetration.Under increasing loads,cracks within ITZs progressively accumulate,culminating in macroscopic fractures that traverse themortarmatrix,forming the complex,serpentine path of cracks.This study reveals four distinct displacement patterns:tensile compliant,tensile opposite,mixed tensile-shear,and shear opposite patterns,each indicative of different stages in concrete’s damage evolution.The widening angle of these patterns delineates the progression of cracks,with the tensile compliant pattern signaling the initial crack appearance and the shear opposite pattern indicating the concrete model’s ultimate failure.

Enhanced UAV Pursuit-Evasion Using Boids Modelling:A Synergistic Integration of Bird Swarm Intelligence and DRL

TheUAV pursuit-evasion problem focuses on the efficient tracking and capture of evading targets using unmanned aerial vehicles(UAVs),which is pivotal in public safety applications,particularly in scenarios involving intrusion monitoring and interception.To address the challenges of data acquisition,real-world deployment,and the limited intelligence of existing algorithms in UAV pursuit-evasion tasks,we propose an innovative swarm intelligencebased UAV pursuit-evasion control framework,namely“Boids Model-based DRL Approach for Pursuit and Escape”(Boids-PE),which synergizes the strengths of swarm intelligence from bio-inspired algorithms and deep reinforcement learning(DRL).The Boids model,which simulates collective behavior through three fundamental rules,separation,alignment,and cohesion,is adopted in our work.By integrating Boids model with the Apollonian Circles algorithm,significant improvements are achieved in capturing UAVs against simple evasion strategies.To further enhance decision-making precision,we incorporate a DRL algorithm to facilitate more accurate strategic planning.We also leverage self-play training to continuously optimize the performance of pursuit UAVs.During experimental evaluation,we meticulously designed both one-on-one and multi-to-one pursuit-evasion scenarios,customizing the state space,action space,and reward function models for each scenario.Extensive simulations,supported by the PyBullet physics engine,validate the effectiveness of our proposed method.The overall results demonstrate that Boids-PE significantly enhance the efficiency and reliability of UAV pursuit-evasion tasks,providing a practical and robust solution for the real-world application of UAV pursuit-evasion missions.