Blockage behavior analysis research method and system based on particle shape influence

Abstract

The invention relates to the technical field of particle blocking research, in particular to a blocking behavior analysis research method and system based on particle shape influence, and the method comprises the following steps: obtaining ellipsoidal particles with different particle shapes, and constructing a first blocking behavior coupling simulation model of the ellipsoidal particles; configuring and verifying the first blockage behavior coupling simulation model to obtain a second blockage behavior coupling simulation model; performing dynamic simulation on the ellipsoidal particles according to the second blockage behavior coupling simulation model to obtain a blockage behavior fitting result; and performing comprehensive analysis on the blockage behavior based on the blockage behavior fitting result, thereby realizing analysis research on the blockage behavior. According to the method, a full-analytic calculation fluid dynamics coupling discrete element method is adopted, and the particle blockage condition in the geometric shapes of the simplified pores of the ellipsoidal particles with different transverse-longitudinal ratios under the fluid flowing condition is simulated. The influence of different particle shapes on the pore blockage characteristics is evaluated through parameter analysis.

Description

Technical Field

[0001] The present invention relates to the technical field of particle plugging research, and in particular to a research method and system for analyzing plugging behavior based on the influence of particle shape. Background Art

[0002] In the prior art, for the particle plugging behavior in porous media under fluid driving conditions, numerical simulation methods are usually adopted to predict and analyze the plugging phenomenon. By constructing a mathematical model of the interaction between fluid and particles, the movement, contact, and plugging process of particles in pores are simulated. However, in order to simplify the calculation and improve the efficiency, the prior art generally uses spherical particle approximation to represent actual particles. This simplified method makes large-scale simulation possible by reducing the complexity of modeling, and reveals the basic laws of particle plugging to a certain extent.

[0003] However, using spherical particle approximation ignores the important influence of particle shape on plugging behavior. The particles in actual engineering often have irregular shapes, such as ellipsoids, polyhedrons, or sheet-like structures, and these shape characteristics play a key role in the interaction between particles and fluid, particle-particle contact, and the formation of plugging structures. For example, the rotation behavior, orientation distribution, and contact mode with pore walls of non-spherical particles are significantly different from those of spherical particles. Therefore, the numerical simulation method based on spherical particle approximation cannot accurately reflect the complex behavior of real particles, resulting in a significant reduction in the accuracy of simulation results and making it difficult to explain the diverse plugging phenomena observed in actual engineering.

[0004] In view of this, the present invention provides a research method and system for analyzing plugging behavior based on the influence of particle shape, which abandons the spherical particle hypothesis and directly simulates the plugging behavior of ellipsoidal particles with different aspect ratios under fluid flow conditions. Through the plugging behavior coupling simulation model, the influence of different particle shapes on plugging behavior is analyzed and studied, revealing the combined influence of particle shape, flow velocity, and pore geometry on plugging events. It not only improves the simulation accuracy of plugging behavior, but also provides a theoretical basis for understanding the plugging and unplugging of non-standard particles in pores, and has important application value for optimizing industrial filtration equipment, preventing formation particle migration plugging, and enhancing geotechnical engineering stability. Summary of the Invention

[0005] Aiming at the defects in the prior art, the present invention provides a research method and system for analyzing plugging behavior based on the influence of particle shape.

[0006] To achieve the above object, in a first aspect, the present invention provides a research method for analyzing plugging behavior affected by particle shape, and the method includes the following steps: obtaining ellipsoidal particles with different particle shapes, and constructing a first coupled simulation model of plugging behavior of the ellipsoidal particles; configuring and verifying the first coupled simulation model of plugging behavior to obtain a second coupled simulation model of plugging behavior; performing dynamic simulation on the ellipsoidal particles according to the second coupled simulation model of plugging behavior to obtain a fitting result of plugging behavior; comprehensively analyzing the plugging behavior based on the fitting result of plugging behavior to achieve the analytical research on the plugging behavior. Through a systematic simulation and analysis process, the present invention significantly improves the scientificity and practicality of the research on particle plugging behavior; the modeling based on ellipsoidal particles with different aspect ratios breaks through the limitation of the traditional spherical particle hypothesis, can accurately describe the differential influence of particle shape, and makes the simulation results closer to the engineering reality; the configuration and verification of the coupled model ensure the reliability of the model and effectively reduce the numerical error; the dynamic simulation of the plugging process based on the coupled simulation model of plugging behavior provides a reference basis for engineering problems such as the optimization of industrial filtration equipment and the prevention and control of formation particle migration.

[0007] Optionally, the step of obtaining ellipsoidal particles with different particle shapes and constructing a first coupled simulation model of plugging behavior of the ellipsoidal particles includes: determining different aspect ratios, and constructing the ellipsoidal particles with different particle shapes based on the aspect ratios; constructing a first control equation for the particle phase and solving a second control equation for the fluid phase; combining the first control equation and the second control equation to construct the first coupled simulation model of plugging behavior of the particle phase and the fluid phase. By introducing parametric design of particle shape, the present invention significantly improves the accuracy of plugging behavior simulation and the ability to reveal the mechanism; constructing ellipsoidal particles based on the aspect ratio breaks through the geometric limitation of the traditional spherical hypothesis, and can systematically analyze the differential influence of different-sized particles on the hydrodynamic response; establishing control equations for the particle phase and the fluid phase respectively realizes the full-scale analysis of fluid-solid interaction; the constructed coupled simulation model can accurately capture the plugging behavior process of ellipsoidal particles and provides an analysis basis for revealing complex plugging mechanisms.

[0008] Optionally, configuring and validating the first clogging behavior coupling simulation model to obtain a second clogging behavior coupling simulation model includes: obtaining the pore geometry parameters and fluid parameters corresponding to the ellipsoidal particles; configuring the parameters of the first clogging behavior coupling simulation model based on the pore geometry parameters and the fluid parameters to obtain model parameters; establishing a model validation test, and validating the first clogging behavior coupling simulation model based on the model validation test in combination with the ellipsoidal particles to obtain the second clogging behavior coupling simulation model. The present invention significantly improves the engineering applicability of the model through parametric configuration and multi-dimensional validation, calibrates the model based on the pore geometry parameters and the fluid parameters, so that the simulation scenario highly fits the actual working conditions; ensures the credibility of the simulation results through the model validation test; provides a benchmark framework for subsequent dynamic simulation and effectively improves the accuracy of dynamic fitting.

[0009] Optionally, comprehensively analyzing the clogging behavior based on the clogging behavior fitting result to realize the analytical research on the clogging behavior, including: analyzing the particle migration of the ellipsoidal particles and the clogging behavior based on the clogging behavior fitting result to obtain a first analysis result; analyzing the influence of the ellipsoidal particles on the clogging behavior in combination with the clogging behavior fitting result to obtain a second analysis result; obtaining the force condition and position distribution of the ellipsoidal particles based on the clogging behavior fitting result as a third analysis result; analyzing the fluid pressure drop in the clogging behavior based on the clogging behavior fitting result to obtain a fourth analysis result; analyzing the influence of the flow velocity and the pore cone angle on the clogging behavior based on the clogging behavior fitting result to obtain a fifth analysis result; realizing the analytical research on the clogging behavior by combining the first analysis result, the second analysis result, the third analysis result, the fourth analysis result and the fifth analysis result. The present invention realizes the analytical research on the clogging behavior through multi-dimensional comprehensive analysis. Based on the fitting result, the particle migration trajectory and the clogging evolution law are analyzed, and the clogging mechanism of the ellipsoidal particles is revealed; through the analysis of the influence of the particle shape, a theoretical basis is provided for the prevention and control of abnormal-shaped particles; combined with the analysis of the clogging structure characteristics, the constitutive conditions of the clogging behavior are clarified; the fluid pressure drop analysis establishes a dynamic correlation between the flow field distribution and the clogging event; through the flow velocity-pore cone angle parameter analysis, the influence on the clogging behavior is analyzed; and a comprehensive analytical research on the clogging behavior is carried out.

[0010] Optionally, analyzing the particle migration and the plugging behavior of the ellipsoidal particles based on the plugging behavior fitting result to obtain a first analysis result, including: in a single pore structure, obtaining the plugging process and the unclogging process of the ellipsoidal particles based on the plugging behavior fitting result; acquiring the flow field variation in the plugging process and the unclogging process. The present invention reveals the dynamic evolution mechanism of the plugging-unclogging of ellipsoidal particles through refined analysis of the single pore structure, accurately captures the whole process of the particle plugging behavior based on the fitting result, and clarifies the influence of the ellipsoidal particles on the flow field, providing a microscopic mechanical basis for analyzing the plugging behavior.

[0011] Optionally, analyzing the influence of the ellipsoidal particles on the plugging behavior by combining the plugging behavior fitting result to obtain a second analysis result, including: according to the plugging behavior fitting result, obtaining the evolution of the particle configuration of the ellipsoidal particles in the plugging behavior; in the plugging behavior, obtaining the variation of the particle erosion rate and the variation of the average particle velocity of the ellipsoidal particles; combining the variation of the particle erosion rate and the variation of the average particle velocity, and conducting a single-particle flushing test on the ellipsoidal particles to obtain a test result. The present invention reveals the dynamic regulation mechanism of the ellipsoidal particles on the plugging behavior through multi-dimensional quantitative analysis; by tracking the evolution of the particle configuration, clarifies the correlation between the particles and the pores during the plugging process; the analysis of the variation of the erosion rate and the average velocity more comprehensively analyzes the plugging behavior; the single-particle flushing test and the analysis of the evolution of the fluid velocity field provide a mechanical basis for optimizing the anti-plugging flushing process.

[0012] Optionally, obtaining the force condition and the position distribution of the ellipsoidal particles according to the plugging behavior fitting result as a third analysis result, including: obtaining the force chain distribution result of the ellipsoidal particles according to the plugging behavior fitting result, and analyzing the force chain distribution result to obtain the force condition; obtaining the plugging particle structure of the ellipsoidal particles according to the plugging behavior fitting result, and using the plugging particle structure as the position distribution. The present invention reveals the internal mechanical mechanism of the plugging structure of the ellipsoidal particles through the force chain analysis of the ellipsoidal particles, obtains the force chain distribution result based on the distribution of the inter-particle contact force extracted from the plugging behavior fitting result, quantifies the spatial configuration and stability of the plugging structure, provides microscopic mechanical evidence for understanding the interlocking effect caused by the orientation arrangement of the ellipsoidal particles with different particle shapes, and at the same time provides a structural weakening point for optimizing the anti-plugging strategy, which helps to deeply analyze the influence of the ellipsoidal particles with different particle shapes on the plugging behavior.

[0013] Optionally, analyzing the fluid pressure drop in the plugging behavior based on the plugging behavior fitting result to obtain a fourth analysis result includes: obtaining the fluid pressure drop of the ellipsoidal particles based on the plugging behavior fitting result; obtaining the evolution relationship between the average particle velocity of the ellipsoidal particles and the fluid pressure drop, and judging the plugging condition of the pores based on the evolution relationship. Through the comparison and verification of the measured and predicted fluid pressure drops, the present invention significantly improves the accuracy of plugging behavior analysis; as a direct quantitative index of the plugging degree, the dynamic change of the fluid pressure drop clearly depicts the characteristics of the plugging evolution stage; analyzing the relationship between the average particle velocity and the fluid pressure drop provides key data support for optimizing the anti-plugging design of porous media.

[0014] Optionally, analyzing the influence of the flow velocity and the pore cone angle on the plugging behavior based on the plugging behavior fitting result to obtain a fifth analysis result includes: determining different flow velocities and different pore cone angles, and obtaining a plugging simulation result based on the flow velocity and the pore cone angle by combining the control variable method; for the ellipsoidal particles, analyzing the influence of the flow velocity on the particle erosion rate under different pore structures to obtain an influence result. Through multi-parameter coupling analysis, the present invention reveals the cooperative regulation mechanism of the flow velocity and the pore cone angle on the plugging behavior; the control variable method can accurately quantify the non-linear influence of the flow velocity on the particle erosion rate; combined with the analysis of the change of the pore cone angle, it provides a quantitative basis for optimizing the anti-plugging operation; it lays a scientific foundation for the design and operation parameter regulation of industrial filtration equipment.

[0015] In a second aspect, the present invention provides a plugging behavior analysis and research system based on the influence of particle shape. The system executes the plugging behavior analysis and research method based on the influence of particle shape provided by the present invention. The system includes an input device, an output device, a processor, and a memory. The input device, the output device, the processor, and the memory are interconnected. Among them, the memory is used to store a computer program, and the computer program includes program instructions. The processor is configured to call the program instructions. Through the collaborative optimization of hardware, the present invention realizes the full-process research of plugging behavior. The high-performance processor and the large-capacity memory work together, significantly shortening the coupling simulation calculation cycle. The high-speed input and output devices ensure the instant interaction of experimental parameters and simulation results. The hardware parallel computing architecture effectively supports efficient response. This integrated design improves the efficiency of the mechanism research of plugging behavior and provides a high-timeliness decision support platform for industrial anti-plugging design. Description of the Drawings

[0016] Figure 1 It is a flowchart of a plugging behavior analysis and research method based on the influence of particle shape according to an embodiment of the present invention;

[0017] Figure 2 It is a schematic diagram of a coupling method according to an embodiment of the present invention;

[0018] Figure 3 Schematic diagram of the coupling simulation model for plugging behavior in the embodiments of the present invention;

[0019] Figure 4 Comparison diagram of the results between the coupling simulation model for plugging behavior and the particle sedimentation test in the embodiments of the present invention;

[0020] Figure 5 Comparison diagram of the results between the coupling simulation model for plugging behavior and the Ergun test in the embodiments of the present invention;

[0021] Figure 6 Configuration diagram of ellipsoidal particles during the simulation process in the embodiments of the present invention;

[0022] Figure 7 Variation trend diagram of particle erosion rate and average particle velocity in the embodiments of the present invention;

[0023] Figure 8 Evolution diagram of particle configuration with different particle shapes during the simulation process in the embodiments of the present invention;

[0024] Figure 9 Variation trend diagram of particle erosion rate and average particle velocity for three particle shapes in the embodiments of the present invention;

[0025] Figure 10 Schematic diagram of the results of the single-particle flushing test in the embodiments of the present invention;

[0026] Figure 11 Schematic diagram of the force chain of plugging particles in the embodiments of the present invention;

[0027] Figure 12 Schematic diagram of the plugging structure of spherical particles in the embodiments of the present invention;

[0028] Figure 13 Schematic diagram of the plugging structure of ellipsoidal particles at different flow velocities in the embodiments of the present invention;

[0029] Figure 14 Relationship diagram between fluid pressure drop and average particle velocity in the embodiments of the present invention;

[0030] Figure 15 Schematic diagram of the simulation results at different flow velocities and orifice cone angles in the embodiments of the present invention;

[0031] Figure 16 Variation trend diagram of particle erosion rate at different flow velocities in the embodiments of the present invention;

[0032] Figure 17 Framework diagram of an analytical research system for plugging behavior based on the influence of particle shape in the embodiments of the present invention. Detailed implementation manners

[0033] Specific embodiments of the present invention will be described in detail below. It should be noted that the embodiments described herein are for illustrative purposes only and are not intended to limit the present invention. In the following description, numerous specific details are set forth in order to provide a thorough understanding of the present invention. However, it will be apparent to those of ordinary skill in the art that the present invention may be practiced without these specific details. In other instances, well-known circuits, software, or methods have not been specifically described to avoid obscuring the present invention.

[0034] Throughout the specification, references to "one embodiment", "an embodiment", "one example", or "an example" mean that a particular feature, structure, or characteristic described in connection with the embodiment or example is included in at least one embodiment of the present invention. Thus, the phrases "in one embodiment", "in an embodiment", "one example", or "an example" appearing throughout the specification do not necessarily all refer to the same embodiment or example. Additionally, the specific features, structures, or characteristics may be combined in any suitable combination and / or sub-combination in one or more embodiments or examples. Moreover, those of ordinary skill in the art should understand that the diagrams provided herein are for illustrative purposes only and are not necessarily drawn to scale.

[0035] Please refer to Figure 1 , an embodiment of the present invention provides a research method for analyzing clogging behavior based on the influence of particle shape, and the method includes the following steps:

[0036] S1. Obtain ellipsoidal particles with different particle shapes and construct a first clogging behavior coupling simulation model for the ellipsoidal particles.

[0037] In this embodiment, different aspect ratios are determined, and ellipsoidal particles with different particle shapes are constructed based on the aspect ratios; a first control equation for the particle phase is constructed, and a second control equation for the fluid phase is solved; a first clogging behavior coupling simulation model for the particle phase and the fluid phase is constructed by combining the first control equation and the second control equation.

[0038] Specifically, different aspect ratios are determined, and the aspect ratios include 1.0, 1.5, and 2.0. It should be noted that when the aspect ratio is 1.0, the ellipsoidal particle is a perfect sphere.

[0039] Specifically, the discrete element method (DEM), which simulates the behavior of dispersed particles based on Newton's second law of motion and the contact model, is used to construct a first control equation for the particle phase, and the first control equation includes the translational motion equation and the rotational motion equation of the ellipsoidal particles.

[0040] The above translational motion equation satisfies the following relationship:

[0041]

[0042] Among them, is the mass of the particle, is the instantaneous velocity, is the total number of contacts of the particle, is the index variable, is the contact force exerted by the particle on the particle, is the interaction force between the particle and the fluid.

[0043] The above rotational motion equation satisfies the following relationship:

[0044]

[0045] Among them, is the inertia moment tensor of the particle, is the angular velocity, is the total number of contacts of the particle, is the index variable, is the torque exerted by the particle on the particle, is the torque between the particle and the fluid.

[0046] In this embodiment, the focus of the research is on single-phase incompressible fluid. The second control equation of the fluid phase is solved by using Computational Fluid Dynamics (CFD for short). The second control equation includes the Continuity Equation and the Navier-Stokes Equation.

[0047] The above continuity equation satisfies the following relationship:

[0048]

[0049] Among them, is the fluid density, is the divergence operator, is the average velocity of the fluid element.

[0050] The above Navier-Stokes equation satisfies the following relationship:

[0051]

[0052] Among them, is the fluid density, is the average velocity of the fluid element, is the divergence operator, represents the direct product of two vectors, is the gradient operator, is the fluid pressure, is the fluid viscosity, is the acceleration of gravity.

[0053] Furthermore, the influence of particle motion on the flow field is achieved by forcing the particle velocity to be equal to the average velocity of the fluid element.

[0054] The above-mentioned particle velocity satisfies the following relationship:

[0055]

[0056] where, is the particle velocity, is the instantaneous velocity, is the angular velocity, is the positive vector relative to the center point of the particle, is the particle domain, indicates inside the particle domain.

[0057] The average velocity of the above-mentioned fluid element satisfies the following relationship:

[0058]

[0059] where, is the average velocity of the fluid element, is the fluid velocity before correction, is the porosity, is the particle velocity, is the simulation domain, indicates inside the simulation domain.

[0060] Please refer to Figure 2 , which shows a schematic diagram of the coupling method. Figure 2 In (a) of Figure 2 is a schematic diagram of the un-resolved coupling method, and in (b) of

[0061] is a schematic diagram of the resolved coupling method; in this embodiment, the resolved coupling method (CFD-DEM method with full resolution) is adopted to simulate the three-dimensional (3D) interaction between fluid and particles. This method can accurately resolve the fluid flow around the particle surface and effectively capture the motion trajectory of the elliptical particle in the flow field.

[0062] Specifically, based on the particle-fluid interaction force and combining the first control equation and the second control equation, a first plugging behavior coupling simulation model (CFD-DEM) for the particle phase and the fluid phase is constructed, including:

[0063] The hydraulic force exerted by the fluid on the particle is obtained by integrating the fluid stress on the particle surface and satisfies the following relationship:

[0064]

[0065]

[0066] Among them, is the traction force vector of the fluid on the particle surface, is the fluid stress of the Newtonian fluid, is the outward normal vector of the particle surface, is the particle surface, indicates on the particle surface, is the fluid pressure, is the unit tensor, is the viscous stress tensor of the fluid, is the fluid viscosity, is the gradient operator, represents the direct product of two vectors, is the average velocity of the fluid element.

[0067] Furthermore, by using the divergence theorem and transferring all hydraulic forces to the particle center point, the hydraulic force received by the particle satisfies the following relationship:

[0068]

[0069]

[0070] Among them, is the resultant hydraulic force vector of the fluid on the particle, is the particle surface, is the traction force vector of the fluid on the particle surface, is the surface area of the particle, is the fluid stress of the Newtonian fluid, is the outward normal vector of the particle surface, is the particle domain, is the divergence operator, is the volume of the particle domain, is the gradient operator, is the fluid pressure, is the viscous stress tensor of the fluid, is the resultant hydraulic moment vector of the fluid on the particle, represents the position vector pointing from the particle center to a certain point on the particle surface.

[0071] In addition, the discretized calculations of the resultant hydraulic force vector and the resultant hydraulic moment vector satisfy the following relationship:

[0072]

[0073]

[0074]

[0075]

[0076] Among them, is the hydrodynamic resultant force vector of the fluid on the particle, is the gradient operator, is the fluid pressure, is the divergence operator, is the viscous stress tensor of the fluid, is the volume of the fluid element, is the porosity, is the particle-fluid interaction force density, is the hydrodynamic resultant torque vector of the fluid on the particle, represents the position vector pointing from the particle center to a point on the particle surface, is the fluid viscosity, is the average velocity of the fluid element.

[0077] S2. Configure and verify the first plugging behavior coupling simulation model to obtain the second plugging behavior coupling simulation model.

[0078] In this embodiment, obtain the pore geometry parameters and fluid parameters corresponding to the ellipsoidal particles; configure the parameters of the first plugging behavior coupling simulation model based on the pore geometry parameters and fluid parameters to obtain the model parameters; establish a model verification test, and verify the first plugging behavior coupling simulation model based on the model verification test in combination with the ellipsoidal particles to obtain the second plugging behavior coupling simulation model.

[0079] Please refer to Figure 3 , which shows a schematic diagram of the plugging behavior coupling simulation model; in this embodiment, as shown in (a) of Figure 3 , in order to study the influence of particle shape on pore plugging and unplugging, three ellipsoidal particles with different aspect ratios are adopted in the CFD-DEM simulation. It should be noted that when the aspect ratio of the ellipsoidal particle reaches 1.0, its geometric shape will become a perfect sphere. As a control case, the influence of particle shape is isolated and quantified by systematic comparison with the ellipsoidal particles.

[0080] In addition, the pore shape has an impact on plugging and unplugging, and the truncated conical pore is used to simplify the pore structure in the porous medium plugging process. For Figure 3 the truncated conical pore shown in (b) of c denoted by D cIt is expressed as (), in order to achieve a predetermined cone angle (expressed as α). Reducing the pore length will increase the cone angle, thereby reducing the pore outlet, squeezing out the particles and hindering the particle removal. On the contrary, increasing the pore contraction diameter will decrease the cone angle, increase the number of particles required for blockage, and thus reduce the possibility of blockage.

[0081] As Figure 3 shown in (c) of [], the pore structure is simplified to a tubular wall, and the outlet is contracted into a cone to induce particle blockage. At a distance of 2 times the particle diameter (expressed as D p from the pore inlet), 125 particles with a diameter of 100 mm are randomly generated in the tubular region of the pore structure. The fluid flows in the positive x-axis direction, and the position of the pore structure is fixed. To avoid the influence of particle volume, the volume of two ellipsoidal particles is equal to the volume of the spherical particle, and then 125 particles are generated in the specified area.

[0082] Specifically, the probability of pore blockage is proportional to the ratio of the pore contraction diameter to the particle diameter, and satisfies the following relationship:

[0083]

[0084] where is the ratio of the pore contraction diameter to the particle diameter, is the pore contraction diameter, is the particle diameter.

[0085] The preliminary simulation results show that when R0 > 3.0, pore blockage rarely occurs, and the present invention does not conduct subsequent simulations on R0 values exceeding this critical value. Therefore, the research focus is on analyzing the blockage behavior within the critical parameter range, where blockage phenomena will occur, and defining the ratio of the pore length to the particle diameter, which satisfies the following relationship:

[0086]

[0087] where is the ratio of the pore length to the particle diameter, is the pore length, is the particle diameter.

[0088] The pore geometric parameters and fluid parameters in the numerical simulation are shown in Table 1:

[0089] Table 1

[0090]

[0091] It should be noted that V f represents the fluid velocity, measured in D p / s. For example, 1D p / s means that the fluid velocity is one particle diameter per second.

[0092] Furthermore, fluid properties such as density and viscosity are assigned based on the properties of water at room temperature. Particle properties, including density, coefficient of friction, and Poisson's ratio, are defined according to the characteristics of quartz particles. In the numerical simulation of particle migration in soil, the Young's modulus of particles is usually set in the range of 10^6 Pa to 10^10 Pa. In this study, a value of 10^6 Pa was selected to balance accuracy and computational efficiency. In addition, the DEM time step was adaptively adjusted to ensure numerical stability and convergence. The model parameters of the first plugging behavior coupling simulation model (CFD-DEM) are shown in Table 2:

[0093] Table 2

[0094]

[0095] Please refer to Figure 4 , which shows the comparison diagram of the results of the plugging behavior coupling simulation model and the particle sedimentation experiment; Figure 4 It shows the comparison of the results of the particle sedimentation velocity in the CFD-DEM model and the experiment, Figure 4 where (a) in Figure 4 is the comparison diagram of the results of spherical particles,

[0096] Furthermore, to verify the ability of the CFD-DEM coupling model to capture the particle-fluid interaction of particle aggregates, the Ergun's Experiment (referred to as Ergun) was simulated. In this model, water flow flushes the particle bed upward from the bottom inlet, and the relationship between the pressure drop and the fluid velocity is as follows:

[0097]

[0098] where is the pressure drop, is the dynamic viscosity of the fluid, is the height of the particle bed, is the porosity of the particle bed, is the particle diameter, is the fluid velocity, is the fluid density.

[0099] Please refer to Figure 5 , which shows the comparison diagram of the results of the plugging behavior coupling simulation model and the Ergun's Experiment; In this study, a particle bed with a size of 100.0 mm × 30.0 mm × 90.0 mm was generated at the bottom of the fluid domain. As Figure 5As shown in (a) therein, the porosity of the particle assembly is set to 0.5, and spherical particles with a diameter of 10.0 mm are filled at the bottom. A constant-velocity fluid is applied at the bottom, and the pressure drop between the inlet and the outlet is monitored and recorded. The inlet flow velocity ranges from 0.02 m / s to 0.06 m / s to verify the accuracy of the model at different flow velocities. As Figure 5 As shown in (b) therein, the simulated pressure drops at different flow velocities and the pressure-drop - fluid velocity relationship show excellent consistency, indicating that the proposed CFD-DEM model can accurately capture the fluid flow behavior within the particle assembly.

[0100] S3. Performing dynamic simulation on the ellipsoidal particles according to the second plugging behavior coupling simulation model to obtain a plugging behavior fitting result.

[0101] In this embodiment, plugging behavior simulation is performed on ellipsoidal particles with different particle sizes according to the constructed second plugging behavior coupling simulation model to obtain a plugging behavior fitting result.

[0102] S4. Based on the plugging behavior fitting result, comprehensively analyzing the plugging behavior to achieve an analytical study of the plugging behavior.

[0103] Among them, S4 specifically includes the following steps:

[0104] S41. Analyzing the particle migration of the ellipsoidal particles and the plugging behavior based on the plugging behavior fitting result to obtain a first analysis result.

[0105] In this embodiment, the simulation results in the case of an aspect ratio of 2.0, R0 = 2.3, R c = 2.0, and V f = 2Dp / s are discussed, and the particle plugging and unclogging processes are illustrated by taking this as an example.

[0106] Please refer to Figure 6 , which shows the configuration diagram of the ellipsoidal particles during the simulation process; a series of snapshots capturing the movement of the ellipsoidal particles are presented, clearly showing the plugging and unclogging processes of the particles in the single-hole structure. As Figure 6 As shown in (a) therein, at t = 0 seconds, it is the initial stage of the simulation process, and the particles are aggregated in a closed area near the pore inlet. They are closely arranged but remain stationary due to zero flow velocity. As Figure 6 As shown in (b) therein, at t = 3.5 seconds, the fluid starts to flow, accompanied by the migration of the particles towards the pore. When the particles continue to migrate towards the outlet, the particles aggregate and interleave in the conical region of the pore, resulting in plugging. The formation of a stable oriented dome structure composed of multiple elongated ellipsoidal particles hinders the outward movement of the particles, and the system enters a quasi-static stage with a low particle velocity. As Figure 6As shown in (c), at t = 5.1 s, due to the continuous flow of the fluid, the particles are subjected to the relevant flow-induced forces, and a few particles near the outlet start to move again, indicating that the force chains between the particles are weakened or rearranged. As Figure 6 shown in (d), at t = 9.0 s, all the particles escape from the pores in sequence. The recovery of particle movement after blockage means that the blockage structure itself is not stable and is easily rearranged by fluid drive.

[0107] Please refer to Figure 7 , which shows the trend diagrams of the particle erosion rate and the average particle velocity; it shows the changes in the particle erosion rate and the average velocity of all particles. As Figure 7 shown in (a), the particle erosion rate remains zero within the first 2 s, indicating that the particles migrate with the fluid within the pore structure and have not reached the outlet yet. After that, within the time interval from 2.0 s to 3.5 s, the particle erosion rate increases in a stepwise manner, indicating that some particles escape from the outlet in sequence without serious blockage. At around 3.5 s, the particle erosion rate stops increasing and enters a relatively long blockage period of 1.6 s, during which no particles are washed out. Subsequently, the particle erosion rate increases rapidly between 5 s and 9 s, marking the start of the main dredging stage, during which most of the particles escape from the pores. During this stage, several intermittent blockage events are also observed. These temporary blockages are caused by the instable structures instantaneously formed between the particles near the pore outlet. As Figure 7 shown in (b), the average velocity of the particles initially rises and then drops significantly around 3.5 s, which corresponds to the pore blockage. After 5 s, as the flow channel reopens, the velocity rises again. The subsequent velocity fluctuations reflect the repeated cycles of blockage and dredging. In short, the blockage process and the particle structure are highly dynamic and are affected by instantaneous flow obstacles.

[0108] Furthermore, analyzing the changes in the flow field during the blockage process and the dredging process, at the beginning of the simulation, the fluid flows smoothly through the pore domain, and the streamlines show a symmetric and orderly pattern and converge towards the pores. During the blockage stage, the flow field is significantly disrupted. Dense particle clusters block the pores, resulting in a sharp drop in the fluid velocity within the blockage area. The streamlines in this area become sparse and chaotic, and at the same time, the upstream flow direction shifts laterally, forming a high-speed jet flow outside the blockage area. Although the blockage restricts the fluid flow, it generates a pressure drop, providing a pressure gradient force for pore dredging. As the dredging progresses, the streamlines gradually recover coherence, and the fluid begins to re-infiltrate the pore channels. Once all the particles escape, the flow field returns to the original pattern.

[0109] S42. Analyze the influence of the ellipsoidal particles on the blockage behavior in combination with the fitting results of the blockage behavior to obtain a second analysis result.

[0110] In this embodiment, to study the effect of particle shape on the clogging behavior, three sets of simulations were conducted using particles with aspect ratios of 1.0, 1.5, and 2.0, respectively. The pore geometry and the inlet fluid velocity were kept the same in all cases, R0 = 2.3, R c = 3.0, V f = 2D p / s.

[0111] Please refer to Figure 8 , which shows the evolution diagram of particle configuration for different particle shapes during the simulation process; it presents the evolution of particle configuration in three representative simulation results with different particle shapes. It can be seen that spherical particles tend to form a relatively dense and stable clogging structure near the outlet, and most particles remain stationary after the initial clogging. In contrast, ellipsoidal particles, especially those with a higher aspect ratio, exhibit more complex dynamic behaviors, including both clogging and unclogging. The elongated shape causes more rearrangement under the action of fluid forces, which not only delays the formation of a stable clogging structure but also introduces intermittent flow paths within the clogging region. It is worth noting that for ellipsoidal particles with an aspect ratio of 2.0, the particles will completely escape at the final stage, and intermittent clogging will occur around 4.5 seconds before that.

[0112] Please refer to Figure 9 , which shows the changing trends of particle erosion rate and average particle velocity for three particle shapes; it displays the changes in the particle erosion rate (denoted by P e ) and the average particle velocity for three particle shapes. As shown in (a) of Figure 9 , the particle erosion rates of spherical and ellipsoidal particles with an aspect ratio of 1.5 reach their peaks near 0.12 and 0.45, respectively, indicating the occurrence of clogging. In contrast, the particle erosion rate of ellipsoidal particles with an aspect ratio of 2.0 continues to increase and stops growing between 3.5 seconds and 5 seconds, indicating particle clogging. Due to the particle shape, unclogging occurs after 5 seconds, and the particle erosion rate continues to rise until it reaches 1.0. It is worth noting that both spherical and ellipsoidal particles with an aspect ratio of 1.5 will experience clogging, but the degree of particle loss varies. As shown in (b) of Figure 9 , which shows the average velocity of all particles, indicating an initial increase. The spherical particles reach their peak velocity at 2.0 seconds and decelerate to a complete stop at 3.0 seconds; the ellipsoidal particles with an aspect ratio of 1.5 also show a similar velocity decay pattern, but their velocity does not reach zero until 5 seconds because, even if the front part of the ellipsoidal particles is clogged, the particles at the rear can still rotate. For ellipsoidal particles with an aspect ratio of 2.0, obvious velocity fluctuations occur after 5 seconds because they cannot form a stable structure to resist the action of the fluid force. As the aspect ratio of the particles increases, the possibility of unclogging increases.

[0113] Please refer to Figure 10 , which is a schematic diagram showing the results of single-particle flushing tests; in this embodiment, as described above, compared with spherical particles, ellipsoidal particles have a lower tendency to block. To further explain this phenomenon, a series of single-particle flushing tests were conducted under the same conditions. As Figure 10 shown in (a) of

[0114] , an ellipsoidal particle was placed at the model inlet at different initial angles. Under the action of hydrodynamic forces, the particle gradually rotates and adjusts its direction during the leftward movement. Generally, when the particle flows out of the pore, its major axis is almost collinear with the flow direction. This behavior of adjusting the direction makes it easier for the ellipsoidal particle to escape from the pore, thus reducing the probability of blockage. Figure 10 Specifically, Figure 10 (a) of Figure 10 shows the particle rotation in different aspect ratio tests;

[0115] (b) of

[0116] shows the change in the angle between the major axis of the particle and the positive X direction when the initial angle is 30°;

[0117] Please refer to Figure 11 , which is a schematic diagram of the force chain of the blocked particles; it shows the force chain network in the blocked particle structure; as Figure 11 (a) of Figure 11 is the force chain distribution result of spherical particles,

[0118] (b) of

[0119] is the force chain distribution result of ellipsoidal particles with an aspect ratio of 1.5. Both spherical particles and ellipsoidal particles form stable force chains during the blocking process, and the front particles bear the main load. Figure 12, which is a schematic diagram of the clogging structure of spherical particles; based on the classification of the clogging structure, the present invention determines three types of spherical particle clogging structures in the pore space, including a single arch structure, a branched arch structure, and a dome structure. Figure 12 In (a) of Figure 12 is a schematic diagram of a single arch structure. Figure 12 In (c) of Figure 12 is a schematic diagram of a dome structure. In (d) of R c = 3.5.

[0120] Furthermore, compared with spherical particles, the clogging structure of ellipsoidal particles is more sensitive to the flow rate. The clogging structures of ellipsoidal particles are classified according to different flow rates.

[0121] Please refer to Figure 13 , which is a schematic diagram of the clogging structure of ellipsoidal particles at different flow rates; at low flow rates, Figure 13 in (a) of Figure 13 is a schematic diagram of the oriented dome structure of the first form of ellipsoidal particles. In (b) of Figure 13 is a schematic diagram of the oriented dome structure of the second form of ellipsoidal particles. Under low flow rate conditions, since the ellipsoidal particles have enough time to adjust their directions, the ellipsoidal particles flow out with their long axes perpendicular to the outlet, thus forming two different forms of oriented dome structures; however, due to the contraction diameter being too small, even at the optimal position, the particles cannot escape. At high flow rates, Figure 13 in (c) of Figure 13 is a schematic diagram of the interlocking bridging structure of the first form of ellipsoidal particles. In (d) of

[0122] S44. Analyze the fluid pressure drop in the clogging behavior based on the fitting result of the clogging behavior to obtain a fourth analysis result.

[0123] The fluid pressure drop is a sensitive indicator of pore clogging. Therefore, the present invention explores the relationship between the fluid pressure drop and particle clogging. The spatial distribution of particles greatly changes the flow field and pressure between the inlet and the outlet.

[0124] Please refer toFigure 14 , shown as a graph of fluid pressure drop versus average particle velocity; at R0 = 2.3, R c = 3.0 and V f = 2D p / s, the pressure drop and average particle velocity of three particle shapes were monitored. Figure 14 In (a) is a graph of the fluid pressure drop versus average particle velocity for spherical particles, Figure 14 in (b) is a graph of the fluid pressure drop versus average particle velocity for ellipsoidal particles (aspect ratio 1.5), Figure 14 in (c) is a graph of the fluid pressure drop versus average particle velocity for ellipsoidal particles (aspect ratio 2.0). During the particle clogging and unclogging processes, there is an obvious negative correlation between these two variables. When particles are trapped in the pores, the fluid flow obstruction increases, resulting in a corresponding increase in the fluid pressure drop. This trend can be observed in all particle shapes, indicating that the fluid pressure drop can reliably show the start and state of clogging and unclogging.

[0125] S45. Analyze the influence of flow velocity and pore cone angle on the clogging behavior based on the fitting results of the clogging behavior to obtain the fifth analysis result.

[0126] Clogging behavior may occur under different flow velocities and pore cone angles; in addition to particle shape, fluid velocity and pore cone angle also have a significant impact on the clogging process. Fluid velocity and pore shape affect the movement of particles, thus affecting potential clogging during the migration process. To systematically evaluate these effects, simulations were carried out at different flow velocities (i.e., V f = 2D p / s, 4D p / s, 6D p / s) and pore cone angles (i.e., α = 23.2°, 25.9°, 26.6°, 27.9°, 29.5°, 31.7°). In the simulations analyzing the influence of flow velocity, the pore cone angle remained constant. Conversely, in the simulations studying the influence of pore cone angle on clogging, the flow velocity remained constant.

[0127] Please refer to Figure 15 , shown as a schematic diagram of the simulation results at different flow velocities and pore cone angles; where the hollow dots indicate particle clogging and the solid dots indicate unclogging. Among the 54 groups of simulations involving three particle shapes, a total of 13 clogging events occurred. The influence of flow velocity on clogging is reflected in the shortening of particle migration time and the greater the flow force at higher flow velocities. Especially when the outlet particles have not formed a stable structure, the strong flow force will destroy the clogging structure. The pore cone angle also affects clogging. When particles move into the pores, the pore walls will hinder their movement, resulting in particle retention. As the pore cone angle increases, this effect becomes more obvious.

[0128] Such as Figure 15As shown in (a) of [reference], 8 clogging events were observed in 18 simulations involving spherical particles. The flow velocity and pore cone angle values were mainly concentrated in V f = 2D p / s, 4D p / s and α = 29.5°, 31.6°.

[0129] As Figure 15 shown in (b) of [reference], ellipsoidal particles with an aspect ratio of 1.5 experienced 5 clogging events in 18 simulations. The clogging events were mainly concentrated under low flow velocity and high pore cone angle conditions.

[0130] As Figure 15 shown in (c) of [reference], ellipsoidal particles with an aspect ratio of 2.0 did not experience any clogging in 18 simulations. All ellipsoidal particles successfully escaped from the pores due to their shape.

[0131] Please refer to Figure 16 , which shows the trend diagram of the particle erosion rate at different flow velocities; it shows the influence of flow velocity on the particle erosion rate of different pore structures, where R c = 3.0, R0 = 2.3, 2.6, 3.0.

[0132] Figure 16 (a) of [reference], Figure 16 (b) of [reference] and Figure 16 (c) of [reference] describe the influence of different flow velocities on the clogging and unclogging of spherical particles under different R0 conditions. In Figure 16 (b) of [reference], clogging occurred at V f = 2D p / s, while all spherical particles escaped from the pore space at V f = 4D p / s and V f = 6D p / s.

[0133] Figure 16 (d) of [reference], Figure 16 (e) of [reference] and Figure 16 (f) of [reference] show the behavior of ellipsoidal particles with an aspect ratio of 1.5 at different flow velocities under different R0 conditions. As the flow velocity increases, the escape rate of the particles also increases, especially in Figure 16 (e) of [reference] and Figure 16 (f) of [reference].

[0134] Figure 16 (g) of [reference], Figure 16 (h) of [reference] and Figure 16 (i) of [reference] show the movement of ellipsoidal particles with an aspect ratio of 2.0 at different flow velocities under different R0 conditions. In Figure 16In (g), blockage occurs at V f = 2D p / s. Due to the instability of the blockage structure, the blockage is cleared after 5 seconds. At higher flow velocities (V f = 4D p / s, 6D p / s), the particle escape process becomes stable. This is because at higher flow velocities, the particle escape speed increases significantly, reducing the possibility of particles remaining in the pore space and generating stronger fluid forces that break the intermittently formed blockage structure.

[0135] S46. Combine the first analysis result, the second analysis result, the third analysis result, the fourth analysis result, and the fifth analysis result to achieve an analytical study of the blockage behavior.

[0136] In this embodiment, to explore the influence of particle shape on pore blockage, a fully analytical CFD-DEM method is used to simulate the interaction between particles and fluids. In addition, a representative model is created by changing the aspect ratio of the particles. To comprehensively analyze the blockage behavior, a detailed comparison is made between spherical particles and ellipsoidal particles with different aspect ratios. In addition, key parameters are also studied, and the summary of the analytical research conclusions includes:

[0137] (1) The aspect ratio of the particles has a positive effect on the pore clearing behavior. Under the same flow conditions, ellipsoidal particles show a higher particle erosion rate. As the aspect ratio increases, the escape effect of these particles becomes more obvious. It should be noted that ellipsoidal particles with an aspect ratio of 2.0 can remain unblocked even under the most extreme simulation conditions.

[0138] (2) Ellipsoidal particles continuously adjust their directions during fluid flow and finally leave the pore with the long axis perpendicular to the outlet. This dynamic direction adjustment is beneficial for ellipsoidal particles to pass through the pore smoothly. However, ellipsoidal particles with an aspect ratio of 1.5 still cause blockage during the simulation. Blockage also occurs when the pore shrinks too narrow to accommodate multiple ellipsoidal particles simultaneously.

[0139] (3) The blockage structure of ellipsoidal particles is greatly affected by the flow velocity. Under low flow velocity conditions, ellipsoidal particles have enough time to adjust their directions so that their long axes are perpendicular to the outlet, thus forming an oriented dome structure. Under high water flow conditions, two structures will form at the outlet: one is an interlaced interlocking bridging structure, and the other is a single arch structure with ellipsoidal particles arranged in parallel. Under these conditions, ellipsoidal particles do not have enough time to optimize the arrangement.

[0140] (4) During the process of particle clogging and unclogging, there is an obvious negative correlation between fluid velocity and fluid pressure drop. When particles are retained in the pores, the obstruction of fluid flow increases, resulting in a corresponding increase in fluid pressure drop. This trend can be observed in all particle shapes, indicating that fluid pressure drop can reliably indicate the start and state of clogging and unclogging.

[0141] The present invention explores the influence of particle shape on pore clogging. The research results provide a theoretical basis for understanding the clogging and unclogging of non-spherical particles in pores. It reveals the pore clogging behavior of non-spherical particles and simulates the transport of particles in porous media under the influence of fluid. It has practical significance for preventing the loss of mechanical stability of soil caused by pore evacuation of non-spherical particles, controlling their migration in strata, and understanding the impact of particle clogging on reservoir permeability.

[0142] Please refer to Figure 17 , in an optional embodiment, the present invention provides a research system for analyzing clogging behavior based on the influence of particle shape. The system includes an input device, an output device, a processor, and a memory, which are interconnected with each other. Among them, the memory is used to store computer programs, and the computer programs include program instructions. The processor is configured to call the program instructions and execute the specific steps of the related embodiments of the method for analyzing clogging behavior based on the influence of particle shape provided by the present invention. The research system for analyzing clogging behavior based on the influence of particle shape provided by the present invention has a complete structure, is objective and stable, and improves the overall applicability and practical application ability of the present invention.

[0143] In summary, the method and system for analyzing clogging behavior based on the influence of particle shape provided by the present invention adopt the full-resolution computational fluid dynamics-discrete element method (CFD-DEM) to simulate the particle clogging situation of ellipsoidal particles with different aspect ratios in a simplified pore geometry under fluid flow conditions, and evaluate the influence of particle shape on pore clogging characteristics through parameter analysis; the method of the present invention is easy to understand, simple to calculate, has less workload, is convenient for engineering application, and provides a theoretical basis and technical support for the further development of the technical field of particle clogging research.

[0144] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, not to limit them; although the present invention has been described in detail with reference to the foregoing embodiments, those of ordinary skill in the art should understand that they can still modify the technical solutions recorded in the foregoing embodiments, or perform equivalent replacements on some or all of the technical features; and these modifications or replacements do not make the essence of the corresponding technical solutions deviate from the scope of the technical solutions of the embodiments of the present invention, and they should all be covered by the scope of the claims and the description of the present invention.

Claims

1. A research method for analyzing clogging behavior based on the influence of particle shape, characterized in that It includes the following steps: Obtain ellipsoidal particles with different particle shapes and construct a first coupled simulation model for the clogging behavior of the ellipsoidal particles; Configure and verify the first coupled simulation model for the clogging behavior to obtain a second coupled simulation model for the clogging behavior; Perform dynamic simulation on the ellipsoidal particles according to the second coupled simulation model for the clogging behavior to obtain a fitting result of the clogging behavior; Conduct a comprehensive analysis of the clogging behavior based on the fitting result of the clogging behavior to realize the analytical research on the clogging behavior.

2. The analytical research method for blocking behavior based on the influence of particle shape according to claim 1, wherein The obtaining of ellipsoidal particles with different particle shapes and the construction of the first coupled simulation model for the clogging behavior of the ellipsoidal particles include: Determine different aspect ratios and construct ellipsoidal particles with different particle shapes based on the aspect ratios; Construct a first control equation for the particle phase and solve a second control equation for the fluid phase; Combine the first control equation and the second control equation to construct the first coupled simulation model for the clogging behavior of the particle phase and the fluid phase.

3. The analytical research method for clogging behavior based on the influence of particle shape according to claim 1, characterized in that The configuration and verification of the first coupled simulation model for the clogging behavior to obtain a second coupled simulation model for the clogging behavior include: Obtain the pore geometry parameters and fluid parameters corresponding to the ellipsoidal particles; Perform parameter configuration on the first coupled simulation model for the clogging behavior based on the pore geometry parameters and the fluid parameters to obtain model parameters; Establish a model verification test, and verify the first coupled simulation model for the clogging behavior based on the model verification test in combination with the ellipsoidal particles to obtain the second coupled simulation model for the clogging behavior.

4. The method for analyzing plugging behavior based on the influence of particle shape according to claim 1, wherein The comprehensive analysis of the clogging behavior based on the fitting result of the clogging behavior to realize the analytical research on the clogging behavior includes: Analyze the particle migration of the ellipsoidal particles and the clogging behavior based on the fitting result of the clogging behavior to obtain a first analysis result; Analyze the influence of the ellipsoidal particles on the clogging behavior in combination with the fitting result of the clogging behavior to obtain a second analysis result; Obtain the force condition and position distribution of the ellipsoidal particles according to the fitting result of the clogging behavior as a third analysis result; Analyze the fluid pressure drop in the clogging behavior based on the fitting result of the clogging behavior to obtain a fourth analysis result; Analyze the influence of the flow velocity and the pore cone angle on the clogging behavior based on the fitting result of the clogging behavior to obtain a fifth analysis result; Realize the analytical research on the clogging behavior by combining the first analysis result, the second analysis result, the third analysis result, the fourth analysis result and the fifth analysis result.

5. The clogging behavior analysis research method based on the influence of particle shape according to claim 4, characterized in that, The analysis of the particle migration of the ellipsoidal particles and the clogging behavior based on the fitting result of the clogging behavior to obtain a first analysis result includes: In a single-hole structure, obtain the clogging process and the unclogging process of the ellipsoidal particles based on the fitting result of the clogging behavior; Obtain the flow field change conditions in the clogging process and the unclogging process.

6. The clogging behavior analysis research method based on the influence of particle shape according to claim 4, characterized in that, The analysis of the influence of the ellipsoidal particles on the clogging behavior in combination with the fitting result of the clogging behavior to obtain a second analysis result includes: According to the fitting result of the clogging behavior, obtain the evolution of the particle configuration of the ellipsoidal particles in the clogging behavior; During the described blocking behavior, the variation of the particle erosion rate and the variation of the average particle velocity of the ellipsoidal particles are obtained; Combining the variation of the particle erosion rate and the variation of the average particle velocity, a single-particle flushing test is performed on the ellipsoidal particles to obtain test results.

7. The clogging behavior analysis research method based on the influence of particle shape according to claim 4, characterized in that, The force condition and position distribution of the ellipsoidal particles obtained according to the fitting result of the blocking behavior as the third analysis result include: The force chain distribution result of the ellipsoidal particles is obtained according to the fitting result of the blocking behavior, and the force condition is obtained by analyzing the force chain distribution result; Based on the fitting result of the blocking behavior, the blocking particle structure of the ellipsoidal particles is obtained, and the blocking particle structure is used as the position distribution.

8. The method for analyzing plugging behavior based on the influence of particle shape according to claim 4, characterized in that The analysis of the fluid pressure drop in the blocking behavior based on the fitting result of the blocking behavior to obtain the fourth analysis result includes: Based on the fitting result of the blocking behavior, the fluid pressure drop of the ellipsoidal particles is obtained; The evolution relationship between the average particle velocity of the ellipsoidal particles and the fluid pressure drop is obtained, and the blocking condition of the pores is judged based on the evolution relationship.

9. The clogging behavior analysis and research method based on the influence of particle shape according to claim 4, wherein, The analysis of the influence of the flow velocity and the pore cone angle on the blocking behavior based on the fitting result of the blocking behavior to obtain the fifth analysis result includes: Different flow velocities and different pore cone angles are determined, and based on the flow velocity and the pore cone angle, a blocking simulation result is obtained by combining the control variable method; For the ellipsoidal particles, the influence of the flow velocity on the particle erosion rate is analyzed under different pore structures to obtain an influence result.

10. A research system for analyzing clogging behavior based on the influence of particle shape, characterized in that The system includes an input device, an output device, a processor, and a memory. The input device, output device, processor, and memory are interconnected. Among them, the memory is used to store a computer program, the computer program includes program instructions, and the processor is configured to call the program instructions to execute the method for analyzing the blocking behavior based on the influence of particle shape as described in any one of claims 1-9.

Images