Method, system, and electronic device for determining a leak diffusion path for a mountainous oil pipeline

Abstract

The present application relates to oil pipeline technology field, especially to a kind of method, system and electronic equipment for determining mountainous area oil pipeline leakage diffusion path, method includes: obtaining the mountain model of mountainous area where oil pipeline is located, and calculating the resultant force of flow resistance and gravity component that each fluid particle is subjected to on the mountain surface of mountainous area where oil pipeline is located, and the acceleration of each fluid particle under the action of corresponding resultant force, determine the grid corresponding to each fluid particle in mountain model after each unit step length time, record and according to the grid corresponding to each fluid particle in mountain model after each unit step length time, obtain the diffusion path of fluid leaked from oil pipeline.When the oil pipeline in mountainous area leaks, the initial velocity of leakage is relatively large, the present application considers the influence of initial velocity of leakage on leakage diffusion path, so that the determined leakage diffusion path is more accurate and efficient.

Description

Technical Field

[0001] This invention relates to the field of oil pipeline technology, and in particular to a method, system, and electronic device for determining the leakage diffusion path of an oil pipeline in a mountainous area. Background Technology

[0002] In the process of transporting crude oil or refined oil over long distances, long-distance oil pipelines have advantages such as large transport capacity, low freight costs, stable transportation, and low energy consumption. However, long-distance oil pipelines span long distances and pass through areas with high mountains and deep valleys, crisscrossing rivers, dense forests, and frequent geological disasters. Once a natural disaster or human-caused damage causes a pipeline leak, a large amount of high-pressure, high-speed leakage will damage the surrounding ecological environment and even threaten the lives and property of the people.

[0003] The national standard "Code for Integrity Management of Oil and Gas Pipelines" (GB 32167-2015) classifies high-consequence zones of oil pipelines into three levels and specifies corresponding distance ranges. The standard also states that "when the terrain near the oil pipeline is highly undulating, the possible flow direction of the leaked oil can be determined based on the terrain, underground culverts, etc., and the distance can be adjusted accordingly." However, it does not specify how to adjust the distance. Currently, the determination of distances to high-consequence zones mainly relies on manual on-site identification, recording, and judgment, which requires a certain level of experience. Therefore, to improve the identification of high-consequence zones in pipeline integrity and to carry out emergency response to pipeline safety accidents, it is necessary to study the dynamic diffusion path and diffusion range of leaked oil under complex surface conditions in mountainous areas.

[0004] Current research on oil pipeline leakage and diffusion under complex terrain conditions mainly relies on CFD simulation technology. CFD simulation primarily uses fluid dynamics equations to control the flow process of oil. CFD simulation results of oil leakage and diffusion on the mountain surface can provide the spatial distribution of oil and the distribution of physical fields during the leakage and diffusion process, providing data support for the study of oil diffusion mechanisms. However, it has drawbacks such as complex modeling and excessively long simulation calculation time, and cannot quickly provide the diffusion path and diffusion range of leaked oil in large-scale spatial conditions.

[0005] Studies on surface leakage and diffusion in mountainous areas have shown that the topographic features of the mountains and the leakage rate of crude oil significantly influence the leakage and diffusion behavior of pipelines in mountainous regions. When the leakage rate is low, crude oil tends to accumulate and flow towards gullies and valleys. Hydrological research, by analyzing the natural law of surface water flow from high to low altitudes, with higher water volume in valleys and lower water volume on ridges, has gradually developed topographic feature extraction algorithms based on physical simulation of surface water flow, such as the SFD8 algorithm. However, the SFD8 algorithm is mainly applied to the extraction of mountain water systems and only considers the main direction of fluid movement under gravity, leading to problems such as inaccurate leakage and diffusion path determination. Summary of the Invention

[0006] The technical problem to be solved by the present invention is to address the shortcomings of the prior art by providing a method, system and electronic equipment for determining the leakage diffusion path of oil pipelines in mountainous areas.

[0007] The technical solution of the method for determining the leakage diffusion path of an oil pipeline in a mountainous area according to the present invention is as follows:

[0008] S1. Obtain the geographic elevation data of the mountainous area where the oil pipeline is located. Divide the image of the mountainous area where the oil pipeline is located into grids based on the geographic elevation data to obtain the mountain model of the mountainous area where the oil pipeline is located. Calculate the resultant force of the flow resistance and gravity component on the surface of the mountain where the oil pipeline is located for each fluid particle, as well as the acceleration of each fluid particle under the corresponding resultant force. Here, the fluid particle represents a preset virtual point position in the fluid leaking from the oil pipeline, and all fluid particles are used to characterize the spatial position of the fluid leaking from the oil pipeline.

[0009] S2. Based on the initial leakage velocity and acceleration of any fluid particle, calculate the position of the fluid particle after a unit step time, and determine the grid corresponding to the position of the fluid particle in the mountain model, until the grid corresponding to each fluid particle in the mountain model after a unit step time is determined.

[0010] S3. Repeat S2, record and obtain the diffusion path of the leaked fluid from the oil pipeline based on the grid corresponding to each fluid particle in the mountain model after each unit step time. The technical solution of the system for determining the diffusion path of a leaked oil pipeline in a mountainous area according to the present invention is as follows:

[0011] The technical solution of the system for determining the diffusion path of leaks in oil pipelines in mountainous areas according to the present invention is as follows:

[0012] It includes a mesh generation calculation module, a calculation determination module, and a repetitive call module;

[0013] The grid division calculation module is used to: acquire the geographic elevation data of the mountainous area where the oil pipeline is located, divide the image of the mountainous area where the oil pipeline is located into grids according to the geographic elevation data, obtain the mountain model of the mountainous area where the oil pipeline is located, and calculate the resultant force of the flow resistance and gravity component on the surface of the mountainous area where the oil pipeline is located for each fluid particle, as well as the acceleration of each fluid particle under the action of the corresponding resultant force. Here, the fluid particle represents the preset virtual point position in the fluid leaking from the oil pipeline, and all fluid particles are used to characterize the spatial position of the fluid leaking from the oil pipeline.

[0014] The calculation and determination module is used to: calculate the position of any fluid particle after a unit step time based on the initial leakage velocity and acceleration of any fluid particle, and determine the grid corresponding to the position of the fluid particle in the mountain model, until the grid corresponding to each fluid particle in the mountain model after a unit step time is determined.

[0015] The repeated call module is used to: repeatedly call the calculation and determination module, record and obtain the diffusion path of the fluid leaking from the oil pipeline based on the grid corresponding to each fluid particle in the mountain model after each unit step time.

[0016] The present invention provides a storage medium storing instructions, which, when read by a computer, cause the computer to execute any of the above-described methods for determining the diffusion path of a leak in an oil pipeline in a mountainous area.

[0017] An electronic device according to the present invention includes a processor and the above-described storage medium, wherein the processor executes instructions in the storage medium.

[0018] The beneficial effects of the technical solution of this invention are as follows:

[0019] When an oil pipeline in a mountainous area leaks, the initial velocity of the leak is relatively high. This invention considers the influence of the initial velocity on the leakage diffusion path, making the determined leakage diffusion path more accurate and efficient. Attached Figure Description

[0020] Figure 1 This is a flowchart illustrating a method for determining the leakage diffusion path of an oil pipeline in a mountainous area, according to an embodiment of the present invention.

[0021] Figure 2 This is a schematic diagram of a mountain model;

[0022] Figure 3 A schematic diagram of the force analysis of a fluid particle on a complex curved surface;

[0023] Figure 4 A schematic diagram of the leakage path at different times;

[0024] Figure 5 Schematic diagram of leakage paths under different surface conditions

[0025] Figure 6 This is a schematic diagram of the structure of a system for determining the leakage diffusion path of an oil pipeline in a mountainous area, according to an embodiment of the present invention. Detailed Implementation

[0026] like Figure 1 As shown in the figure, a method for determining the leakage diffusion path of an oil pipeline in a mountainous area according to an embodiment of the present invention includes the following steps:

[0027] S1. Obtain the geographic elevation data of the mountainous area where the oil pipeline is located. Based on the geographic elevation data, divide the image of the mountainous area into grids to obtain the mountain model. Calculate the resultant force of the flow resistance and gravity component on the surface of the mountain where the oil pipeline is located for each fluid particle, as well as the acceleration of each fluid particle under the corresponding resultant force. Here, the fluid particle represents a preset virtual point position in the fluid leaking from the oil pipeline, and all fluid particles are used to characterize the spatial location of the fluid leaking from the oil pipeline.

[0028] The specific process of obtaining the geographic elevation data of the mountainous area where the oil pipeline is located is as follows:

[0029] Geographic elevation data of the mountainous area where the oil pipeline is located can be obtained from the SRTM open-source database, ASTER open-source database, or other open-source databases. The geographic elevation data of the mountainous area where the oil pipeline is located can realize the digital simulation of the mountainous area where the oil pipeline is located, and the ground elevation of the mountainous area where the oil pipeline is located can be represented by a set of ordered numerical arrays.

[0030] Specifically, the image of the mountainous area where the oil pipeline is located is divided into grids based on geographic elevation data to obtain a mountain model of the area. The specific process is as follows:

[0031] Each elevation data point in the geographic elevation data is used as the elevation data point for each node in each grid, generating an elevation grid model where each grid area is a square, i.e., a mountain model of the mountainous area where the oil pipeline is located. The size of each grid can be set according to the actual situation.

[0032] The specific process for calculating the resultant force of flow resistance and gravitational component experienced by each fluid particle on the mountain surface in the mountainous area where the oil pipeline is located is as follows:

[0033] Based on the slope of the mountain where any fluid particle is located, determine the gravity component of the fluid particle, calculate the flow resistance experienced by the fluid particle according to the existing Chezy formula, and obtain the resultant force of the flow resistance and gravity component of the fluid particle, until the resultant force of the flow resistance and gravity component of each fluid particle is obtained, and then calculate the acceleration of each fluid particle under the action of the corresponding resultant force.

[0034] S2. Based on the initial leakage velocity and acceleration of any fluid particle, calculate the position of the fluid particle after a unit step time, and determine the grid corresponding to the position of the fluid particle in the mountain model, until the grid corresponding to each fluid particle in the mountain model after a unit step time is determined.

[0035] Specifically, based on the initial leakage velocity and acceleration of any fluid particle, the position of the fluid particle after a unit step time is calculated.

[0036] When initializing the positions of fluid particles, the initial distribution of the fluid particles in the mountain model is determined based on the location of the pipeline rupture point. Simultaneously, based on the operating pressure of the pipeline, the initial velocity of the fluid particles, i.e., the initial leakage velocity, is calculated using Bernoulli's equation. The position of the fluid particles after a unit step time is then calculated using kinematics. S3: Repeat S2, recording and determining the diffusion path of the leaked fluid from the pipeline based on the mesh corresponding to each fluid particle in the mountain model after each unit step time.

[0037] Optionally, the above technical solution also includes:

[0038] S4. Display the diffusion path of the leaked fluid from the oil pipeline in the mountain model as a dynamic diagram.

[0039] Optionally, the above technical solution also includes: interpolating the mesh on the mountain model to increase the resolution of the mesh on the mountain model, so as to more realistically display the diffusion path of the fluid leaking from the oil pipeline.

[0040] Optionally, the above technical solution also includes:

[0041] S5. Make advance preparations for the area covered by the diffusion path of the leaked fluid from the oil pipeline. In areas where the leaked oil is likely to flow, build earthen dams or dig pits to block the continued flow of the leaked oil along the ground surface. Also, prepare oil booms, oil booms, skimmers and other devices in advance in the area to reduce the losses from the leak accident.

[0042] This patent aims to provide a method for rapidly calculating the diffusion path of oil leaks in mountainous oil pipelines. It incorporates factors affecting leak flow, such as leak flow rate, leak velocity, and surface roughness, enabling a more scientific assessment of the impact range of leaked oil. The current method can quickly calculate the main oil diffusion path given DEM data of the mountain and the leak velocity. This method mainly includes the following aspects:

[0043] 1) Real terrain modeling (boundary conditions): This paper models real mountains using a digital elevation model. A digital elevation model can digitally simulate the ground terrain using limited terrain elevation data. It represents the ground elevation in the form of an ordered numerical array.

[0044] 2) Fluid Particle Force Analysis (Governing Equations): To simplify calculations, this paper simplifies the forces acting on fluid particles, considering only gravity and flow resistance during motion. The component of gravity along the mountain surface is the primary force driving the fluid particle motion. The flow resistance is calculated primarily using the Chezy formula (widely used for calculating head loss in open channel and surface flow). First, the head loss during the flow process is calculated, and then the average flow resistance is derived from the head loss.

[0045] 3) Provide relevant information about the leak (initial conditions): The algorithm requires initial conditions including the location of the leak point, the magnitude and direction of the leak velocity.

[0046] 4) Calculate the main path of oil spill diffusion (computational solution): The oil spill diffusion process can be regarded as the oil undergoing variable acceleration motion on the mountain surface, while the acceleration in each grid cell is constant. The position and velocity of the oil particle after time t can be obtained by using the constant acceleration formula and a specified time step. Repeating this calculation process will yield the main path of oil spill diffusion.

[0047] 5) Visualization of the main diffusion path: By taking a small time step and connecting the oil location points at different times, a relatively smooth oil spill diffusion path can be obtained. Displaying this path on a realistic terrain model can vividly demonstrate the dynamic diffusion process of oil after the leak occurs.

[0048] The following embodiment illustrates a method for determining the leakage diffusion path of an oil pipeline in a mountainous area according to the present invention:

[0049] S10. This invention models real mountainous terrain using a Digital Elevation Model (DEM). A DEM can digitally simulate ground topography using limited terrain elevation data, representing ground elevation through an ordered numerical array. The elevation data in this invention are sourced from open-source databases such as SRTM and ASTER. By interpolating the studied mountain elevation data to generate a surface, a realistic and complex ground model of the oil pipeline can be obtained. Figure 2 As shown.

[0050] S11. The dynamic model of the migration of leaked oil on the mountain surface is relatively complex. To simplify the calculation, the research method of this invention simplifies the forces acting on fluid particles, considering only the gravity and flow resistance experienced by the fluid particles during the movement. For example... Figure 3As shown, a fluid particle moves on a curved surface with velocity V, and the normal vector of the surface is n. G and f are the gravity and flow resistance acting on the fluid particle, respectively. G' is the component of gravity along the direction of the fastest descent along the surface, and also the component of gravity G on the tangent plane of the surface. F is the resultant force of the gravity component G' and the flow resistance f. The components of G' along the rows and columns of the numerical matrix are G... x and G y The calculation is as follows

[0051]

[0052] In the formula, G x,(i,j) With G y(i,j) These represent the components of gravity along the direction of fastest descent along the surface in the row and column directions of the numerical array, respectively, with units of N; DEM (i,j) represents the elevation of the location represented by each value in the numerical array of the digital elevation model, in meters (m); m represents the mass of a single fluid particle, in kilograms (kg); g represents the acceleration due to gravity, in meters per second (m / s²). 2 ; l represents the accuracy of the digital elevation model, m.

[0053] The second factor affecting the motion of fluid particles is flow resistance. During the process of oil moving along the mountain after a leak, the flow of oil is similar to open channel overflow or surface overflow. The frictional resistance experienced by the oil during the flow can be calculated using Formula 2 (also known as the Cheze formula).

[0054]

[0055] In the formula, h f ρ is the head loss of the fluid along the flow path, m; s is the distance the fluid travels, m; V is the average velocity of the fluid during the flow process, m / s; C is the Chezy coefficient, m. 0.5 / s -1 R is the hydraulic radius, in meters; the Chezy coefficient is calculated using the Manning formula.

[0056]

[0057] In the formula, n is the roughness coefficient, representing the roughness between the flowing liquid and the ground surface. For smooth concrete, sparsely grassed soil surfaces, heavily grassed soil surfaces, and light shrubs, n can be taken as 0.012, 0.07, 0.1, and 0.35, respectively. Formula 4 for calculating flow resistance can be derived from formulas 2 and 3. Flow resistance is mainly related to the fluid flow velocity, the hydraulic radius of the fluid, and the roughness coefficient between the fluid and the mountain surface.

[0058]

[0059] S12. Give the location of the leak (x0, y0) and the initial leakage velocity V. 0,x V 0,y .

[0060] S13. Calculate the resultant force of the fluid particles using formulas 1 and 4, and determine the resultant acceleration along the x and y directions as a. x and a y The velocity and distance traveled after time Δt are calculated using Equations 5 and 6, respectively.

[0061] V t,x =V 0,x +a x ·Δt,V t,y =V 0,y +a y ·Δt (5)

[0062]

[0063] Then the location of the fluid particle after time Δt is:

[0064] x = x0 + s t,x ,y=y0+s t,y (7)

[0065] S14. Repeat step S13 to obtain the positions of fluid particles at different times. Connecting the fluid position points yields a dynamic leakage path diagram of the leaking oil, such as... Figure 4 As shown. Furthermore, the roughness coefficient in Formula 3 can be modified to simulate the impact of different seasons on the mountain surface on the leakage path, such as... Figure 5 As shown.

[0066] In the above embodiments, although the steps are numbered S1, S2, etc., they are only specific embodiments given in this application. Those skilled in the art can adjust the execution order of S1, S2, etc. according to the actual situation, which is also within the protection scope of this invention. It can be understood that in some embodiments, some or all of the above embodiments may be included.

[0067] like Figure 6 As shown, a system 200 for determining the leakage diffusion path of an oil pipeline in a mountainous area according to an embodiment of the present invention includes a grid division calculation module 210, a calculation and determination module 220, and a repeated call module 230.

[0068] The grid division calculation module 210 is used to: acquire the geographic elevation data of the mountainous area where the oil pipeline is located, divide the image of the mountainous area where the oil pipeline is located into grids based on the geographic elevation data, obtain the mountain model of the mountainous area where the oil pipeline is located, and calculate the resultant force of the flow resistance and gravity component on the surface of the mountainous area where the oil pipeline is located for each fluid particle, as well as the acceleration of each fluid particle under the action of the corresponding resultant force. Here, the fluid particle represents the preset virtual point position in the fluid leaking from the oil pipeline, and all fluid particles are used to characterize the spatial position of the fluid leaking from the oil pipeline.

[0069] The calculation and determination module 220 is used to: calculate the position of any fluid particle after a unit step time based on the initial leakage velocity and acceleration of any fluid particle, and determine the grid corresponding to the position of the fluid particle in the mountain model, until the grid corresponding to each fluid particle in the mountain model after a unit step time is determined.

[0070] The repeated call module 230 is used to repeatedly call the calculation and determination module 220 to record and obtain the diffusion path of the fluid leaking from the oil pipeline based on the grid corresponding to each fluid particle in the mountain model after each unit step time.

[0071] Optionally, the above technical solution also includes a display module, which is used to display the diffusion path of the fluid leaking from the oil pipeline in a dynamic diagram on the mountain model.

[0072] Optionally, in the above technical solution, the mesh generation calculation module 210 is also used to: perform interpolation processing on the mesh on the mountain model to increase the resolution of the mesh on the mountain model.

[0073] Optionally, in the above technical solution, the mesh generation calculation module 210 is specifically used for:

[0074] Obtain geographical elevation data of the mountainous area where the oil pipeline is located from the SRTM or ASTER open-source database.

[0075] The parameters and steps of each unit module in the system 200 for determining the leakage diffusion path of an oil pipeline in a mountainous area, as described above, can be referred to the parameters and steps in the embodiments of the method for determining the leakage diffusion path of an oil pipeline in a mountainous area described above, and will not be repeated here.

[0076] An embodiment of the present invention provides a storage medium storing instructions, which, when read by a computer, cause the computer to execute any of the above-described methods for determining the diffusion path of a leak in an oil pipeline in a mountainous area.

[0077] An electronic device according to an embodiment of the present invention includes a processor and the aforementioned storage medium. The processor executes instructions in the storage medium. The electronic device may be a computer, a mobile phone, or the like.

[0078] Those skilled in the art will know that this invention can be implemented as a system, method, or computer program product.

[0079] Therefore, this disclosure can be implemented in the following forms: it can be entirely hardware, entirely software (including firmware, resident software, microcode, etc.), or a combination of hardware and software, generally referred to herein as a "circuit," "module," or "system." Furthermore, in some embodiments, the invention can also be implemented as a computer program product in one or more computer-readable media, the computer-readable medium containing computer-readable program code.

[0080] Any combination of one or more computer-readable media may be used. A computer-readable medium can be a computer-readable signal medium or a computer-readable storage medium. A computer-readable storage medium can be, for example, but not limited to, an electrical, magnetic, optical, electromagnetic, infrared, or semiconductor system, apparatus, or device, or any combination thereof. More specific examples (a non-exhaustive list) of computer-readable storage media include: an electrical connection having one or more wires, a portable computer disk, a hard disk, random access memory (RAM), read-only memory (ROM), erasable programmable read-only memory (EPROM or flash memory), optical fiber, portable compact disk read-only memory (CD-ROM), optical storage device, magnetic storage device, or any suitable combination thereof. In this document, a computer-readable storage medium can be any tangible medium that contains or stores a program that can be used by or in connection with an instruction execution system, apparatus, or device.

[0081] Although embodiments of the present invention have been shown and described above, it is understood that the above embodiments are exemplary and should not be construed as limiting the present invention. Those skilled in the art can make changes, modifications, substitutions and variations to the above embodiments within the scope of the present invention.

Claims

1. A method for determining the diffusion path of a leak in an oil pipeline in a mountainous area, characterized in that, include: S1. Obtain the geographic elevation data of the mountainous area where the oil pipeline is located. Divide the image of the mountainous area where the oil pipeline is located into grids based on the geographic elevation data to obtain the mountain model of the mountainous area where the oil pipeline is located. Calculate the resultant force of the flow resistance and gravity component on the surface of the mountain where the oil pipeline is located for each fluid particle, as well as the acceleration of each fluid particle under the corresponding resultant force. Here, the fluid particle represents a preset virtual point position in the fluid leaking from the oil pipeline, and all fluid particles are used to characterize the spatial position of the fluid leaking from the oil pipeline. S2. Based on the initial leakage velocity and acceleration of any fluid particle, calculate the position of the fluid particle after a unit step time, and determine the grid corresponding to the position of the fluid particle in the mountain model, until the grid corresponding to each fluid particle in the mountain model after a unit step time is determined. S3. Repeat S2, record and obtain the diffusion path of the fluid leaking from the oil pipeline based on the grid corresponding to each fluid particle in the mountain model after each unit step time. The mountain model of the mountainous area where the oil pipeline is located is constructed by digitally simulating the real mountainous terrain using a digital elevation model. When acquiring the geographic elevation data of the mountainous area where the oil pipeline is located, the geographic elevation data is obtained from the SRTM or ASTER open-source database. This geographic elevation data represents the ground elevation of the mountainous area in the form of an ordered numerical array. Based on the geographic elevation data, the image of the mountainous area where the oil pipeline is located is divided into grids, with each elevation data point used as the elevation data for each node of each grid, generating a square elevation grid model for each grid region. The elevation grid model is a mountain model of the mountainous area where the oil pipeline is located. Interpolation processing is performed on the grid in the mountain model to increase its resolution. The mountain model provides the position information and elevation data of each fluid particle moving on the mountain surface. The fluid particles are preset virtual points in the fluid leaking from the oil pipeline, and all fluid particles represent the spatial location of the leaking fluid. Based on the elevation data corresponding to each grid in the mountain model, the slope of the mountain at the location of each fluid particle is determined, and then the gravitational component experienced by each fluid particle on the mountain surface is calculated. The mountain model is also used to determine the grid corresponding to the fluid particle in the mountain model based on the position of the fluid particle after each unit step time during the calculation of the motion of each fluid particle, and to record the grid corresponding to each fluid particle in the mountain model after each unit step time, so as to obtain the diffusion path of the fluid leaking from the oil pipeline.

2. The method for determining the leakage diffusion path of an oil pipeline in a mountainous area according to claim 1, characterized in that, Also includes: The diffusion path of the fluid leaking from the oil pipeline is displayed in the mountain model as a dynamic diagram.

3. The method for determining the leakage diffusion path of an oil pipeline in a mountainous area according to claim 1, characterized in that, Also includes: Interpolation processing is performed on the mesh on the mountain model to increase the resolution of the mesh on the mountain model.

4. The method for determining the leakage diffusion path of an oil pipeline in a mountainous area according to claim 1, characterized in that, The acquisition of geographic elevation data for the mountainous area where the oil pipeline is located includes: Obtain the geographical elevation data of the mountainous area where the oil pipeline is located from the SRTM open-source database or the ASTER open-source database.

5. A system for determining the diffusion path of a leak in an oil pipeline in a mountainous area, characterized in that, It includes a mesh generation calculation module, a calculation determination module, and a repetitive call module; The grid division calculation module is used to: acquire the geographic elevation data of the mountainous area where the oil pipeline is located, divide the image of the mountainous area where the oil pipeline is located into grids according to the geographic elevation data, obtain the mountain model of the mountainous area where the oil pipeline is located, and calculate the resultant force of the flow resistance and gravity component on the surface of the mountainous area where the oil pipeline is located for each fluid particle, as well as the acceleration of each fluid particle under the action of the corresponding resultant force. Here, the fluid particle represents the preset virtual point position in the fluid leaking from the oil pipeline, and all fluid particles are used to characterize the spatial position of the fluid leaking from the oil pipeline. The calculation and determination module is used to: calculate the position of any fluid particle after a unit step time based on the initial leakage velocity and acceleration of any fluid particle, and determine the grid corresponding to the position of the fluid particle in the mountain model, until the grid corresponding to each fluid particle in the mountain model after a unit step time is determined. The repeated call module is used to: repeatedly call the calculation and determination module, record and obtain the diffusion path of the fluid leaking from the oil pipeline based on the grid corresponding to each fluid particle in the mountain model after each unit step time; The mountain model of the mountainous area where the oil pipeline is located is constructed by digitally simulating the real mountainous terrain using a digital elevation model (DEM). When acquiring the geographic elevation data of the mountainous area, the DEM is obtained from the SRTM or ASTER open-source database. This geographic elevation data represents the ground elevation of the mountainous area as an ordered numerical array. Based on the geographic elevation data, the image of the mountainous area is divided into grids. Each elevation data point is used as the elevation data for each node of each grid, generating a square elevation grid model. This elevation grid model constitutes the mountain model of the mountainous area where the oil pipeline is located. Interpolation processing is performed on the grids in the mountain model to increase the density of the grids. The resolution; the mountain model is used to provide the position information and elevation data of each fluid particle moving on the mountain surface. The fluid particles are preset virtual points in the fluid leaking from the oil pipeline. All fluid particles are used to characterize the spatial position of the fluid leaking from the oil pipeline. Based on the elevation data corresponding to each grid in the mountain model, the slope of the mountain at the position of each fluid particle is determined, and then the gravitational component of each fluid particle on the mountain surface is calculated. The mountain model is also used to determine the grid corresponding to the fluid particle in the mountain model based on the position of the fluid particle after each unit step time during the calculation of the movement of each fluid particle, and to record the grid corresponding to each fluid particle in the mountain model after each unit step time, so as to obtain the diffusion path of the fluid leaking from the oil pipeline.

6. The system for determining the leakage diffusion path of an oil pipeline in a mountainous area according to claim 5, characterized in that, It also includes a display module, which is used to display the diffusion path of the fluid leaking from the oil pipeline in the mountain model in the form of a dynamic diagram.

7. The system for determining the leakage diffusion path of an oil pipeline in a mountainous area according to claim 5, characterized in that, The grid division calculation module is also used to: perform interpolation processing on the grid on the mountain model to increase the resolution of the grid on the mountain model.

8. A system for determining the diffusion path of a leak in an oil pipeline in a mountainous area according to claim 5, characterized in that, The mesh generation calculation module is specifically used for: Obtain the geographical elevation data of the mountainous area where the oil pipeline is located from the SRTM open-source database or the ASTER open-source database.

9. A storage medium, characterized in that, The storage medium stores instructions that, when read by a computer, cause the computer to execute a method for determining the leakage diffusion path of an oil pipeline in a mountainous area as described in any one of claims 1 to 4.

10. An electronic device, characterized in that, It includes a processor and the storage medium of claim 9, wherein the processor executes instructions in the storage medium.

Images