A rapid oil spill simulation method for submarine pipeline oil leakage accidents

Abstract

The present invention discloses a method for rapid simulation of oil spills applied to submarine pipeline oil leakage accidents. First, the present invention adds a module for reading three-dimensional flow field data on the basis of a conventional two-dimensional oil spill model on the ocean surface; secondly, designs a stress-driven algorithm for the sea surface wind field on oil particles at different water depths below the sea surface; thirdly, designs an algorithm for the floating velocity of oil particles below the sea surface; fourthly, designs an algorithm for the vertical mixing layer height and vertical diffusion coefficient of the sea surface; finally, combines flow field forecast data and wind field forecast data to drive the calculation of a three-dimensional oil spill drift and diffusion model. On the basis of a conventional two-dimensional oil spill model, the present invention adds the vertical physical change process of oil particles, and by loading three-dimensional flow field data, realizes the simulation calculation of related physical processes such as floating, drifting and diffusion of oil particles in the three-dimensional water space from the seabed to the sea surface, which provides reference value for rapid simulation and emergency deduction decision-making of oil pollutant leakage in marine pollution.

Description

Technical Field

[0001] The invention belongs to the technical field of environmental monitoring, and in particular relates to a rapid oil spill simulation method applied to submarine pipeline oil leakage accidents. Background Art

[0002] Oil spills often occur over vast expanses of water with complex hydrological and meteorological conditions. Monitoring alone makes it difficult to quickly assess the impact of an oil spill and predict the development of oil film pollution. Currently, leading oil spill models, both domestically and internationally, employ a gridless Lagrangian particle tracing algorithm, treating oil spill contaminants as a collection of a fixed number of oil particles. These models directly incorporate flow and wind field forecasts, capturing real-time flow and wind speed data around oil film particles. This data is supplemented by a fixed shoreline to determine the landing status of the oil film particles. This simulation algorithm is computationally lightweight and fast, requiring only minutes to compute.

[0003] Existing oil spill simulation technology mainly targets oil spills on the ocean surface and is a numerical model of a two-dimensional plane field. It does not take into account the related physical processes such as the rising, drifting and diffusion of the oil film from the seabed to the sea surface in sudden submarine oil spill pollution accidents such as leakage from submarine oil pipelines and leakage from shipwreck oil tanks. Summary of the Invention

[0004] In view of the deficiencies in the prior art, the present invention provides a rapid oil spill simulation method for submarine pipeline oil leakage accidents.

[0005] A first aspect of the present invention provides a method for rapid simulation of oil spills in submarine pipeline oil leak accidents, the method comprising the following steps:

[0006] Step 1: Read 3D flow and wind field data:

[0007] Based on the two-dimensional oil spill model on the ocean surface, a module for reading three-dimensional flow field and wind field data is added;

[0008] Read 3D flow field data, confirm the flow field file structure, and extract velocity components, water level and water depth information;

[0009] Read wind field data, confirm the wind field file structure, and extract sea surface wind speed information;

[0010] Step 2: Calculate the force and motion parameters of the oil particles:

[0011] Based on the read three-dimensional flow field data, the drift velocity component of the oil particles driven by the ocean current in the three-dimensional water space is calculated;

[0012] Based on the read wind field data, the velocity components and deflection angles of oil particles driven by wind stress at different water depths are calculated;

[0013] Calculate the rising velocity of oil particles below the sea surface by combining the physical properties of oil particles (such as density and particle size) with the properties of water (such as water density and viscosity);

[0014] Calculate the vertical diffusion coefficient at different water depths below the sea surface based on water depth, surface wave wavelength, significant wave period, and breaking wave height;

[0015] Step 3: Simulate the three-dimensional movement of oil particles:

[0016] Combining 3D flow and wind field data with calculated oil particle forces and motion parameters to drive a 3D oil spill drift and diffusion model;

[0017] Based on the calculation time step, the actual displacement of each oil particle at each time step is calculated to simulate the floating, drifting and diffusion process of oil particles from the seabed to the sea surface.

[0018] In a second aspect of the present invention, a computer device is provided, comprising a memory, a processor, and computer executable instructions stored in the memory and executable on the processor. When the processor executes the instructions, the method for rapid oil spill simulation applied to submarine pipeline oil leakage accidents is implemented.

[0019] A third aspect of the present invention provides a computer storage medium storing computer executable instructions, which, when executed, implement the method for rapid oil spill simulation applied to submarine pipeline oil leakage accidents.

[0020] The present invention has the following beneficial effects: Existing mainstream oil spill simulation technologies are primarily focused on numerical simulations of surface oil spills, and are unable to accurately depict the three-dimensional (3D) release of oil from sudden submarine oil spills, such as leaks from submarine pipelines and shipwreck tanks. This invention addresses these shortcomings in existing oil spill prediction technology by simulating the physical processes involved in the rise, drift, and diffusion of oil particles in the three-dimensional water space between the seabed and the surface. BRIEF DESCRIPTION OF THE DRAWINGS

[0021] Figure 1 is a flow chart of a method according to an embodiment of the present application;

[0022] Figure 2 This is a visualization example diagram of the simulation of the vertical drift and diffusion process of petroleum pollutants in an embodiment of the present application;

[0023] Figure 3 This is a visualization example diagram of the simulation of the horizontal drift and diffusion process of petroleum pollutants in an embodiment of the present application. DETAILED DESCRIPTION

[0024] The following will clearly and completely describe the technical solutions in the embodiments of the present invention in conjunction with the accompanying drawings. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without making creative efforts are within the scope of protection of the present invention.

[0025] This application is based on the internationally widely used Lagrangian oil particle algorithm. On the basis of the conventional two-dimensional oil spill model, the vertical physical change process of oil particles is added. By loading three-dimensional flow field data, the simulation calculation of related physical processes such as the floating, drifting and diffusion of oil particles in the three-dimensional water space between the seabed and the sea surface is realized, thereby designing a rapid oil spill simulation method that can be applied to submarine pipeline oil leakage accidents.

[0026] like Figure 1 As shown, the embodiment of the present application provides a method for rapid simulation of oil spills applied to submarine pipeline oil leakage accidents, comprising the following steps:

[0027] Step S1: Based on the conventional two-dimensional ocean surface oil spill model, a module for reading three-dimensional flow field data is added.

[0028] In a preferred example, the steps include:

[0029] (1) Based on the data format protocol of three-dimensional flow field data in NetCDF format, read the header file of flow field data and confirm the flow field file structure;

[0030] (2) Identify the parameters representing the velocity components (three dimensions), water level, and water depth in the flow field file;

[0031] (3) Based on the position of each particle in the three-dimensional space of the water body and the three-dimensional coordinates of each grid point of the flow field data, the three-dimensional drift velocity of the oil particles driven by the ocean current is calculated using distance inverse weighted interpolation.

[0032] Step S2: Design a stress-driven algorithm for the sea surface wind field on oil particles at different water depths below the sea surface.

[0033] In a preferred example, the steps include:

[0034] (1) The classical formula for the effect of wind stress on water stress at different water depths is used to calculate the degree to which oil particles at different water depths are affected by wind stress:

[0035] (1)

[0036] Where, represents the wind stress surface driving factor of oil particles, to It is the partial coefficient of the regression curve of the empirical formula, which needs to be determined through field experiments or reference literature. Indicates the sea surface wind speed value, Indicates the water depth where the oil particles are located.

[0037] (2)

[0038] Where, represents the sub-deep driving factor of wind stress on oil particles, to It is the partial coefficient of the regression curve of the empirical formula. Its value changes with the change of the water depth range of the oil particles. It needs to be determined through field experiments or references. Indicates the sea surface wind speed value.

[0039] (3)

[0040] Where, represents the comprehensive impact factor of wind stress on oil particles, represents the wind stress surface driving factor of oil particles, Represents the sub-deep driving factor of wind stress on oil particles.

[0041] (2) The classic formula for the deflection angle of water flow at different water depths due to wind stress is used to calculate the deflection angle of oil particles at different water depths driven by wind stress:

[0042] (4)

[0043] Where, Indicates the influencing factor of the surface deflection angle of oil particle drift driven by wind stress, to It is the partial coefficient of the regression curve of the empirical formula, which needs to be determined through field experiments or reference literature. Indicates the sea surface wind speed value, Indicates the water depth where the oil particles are located.

[0044] (5)

[0045] Where, Influencing factors of the deep deflection angle of oil particles drifting driven by wind stress, Indicates the sea surface wind speed value.

[0046] (6)

[0047] Where, Influencing factors of the sub-deep deflection angle of oil particles drifting driven by wind stress, to It is the partial coefficient of the regression curve of the empirical formula. Its value changes with the change of the water depth range of the oil particles. It needs to be determined through field experiments or reference literature. Indicates the sea surface wind speed value.

[0048] (7)

[0049] Where, represents the deflection angle of oil particle drift driven by wind stress relative to the wind vector at different water depths, Indicates the influencing factor of the surface deflection angle of oil particle drift driven by wind stress, Influencing factor of the sub-deep deflection angle of oil particles drifting driven by wind stress.

[0050] (3) Based on the above calculation results and the water surface wind field data, the wind stress-driven oil particle drift velocity components at different water depths are calculated.

[0051] (8)

[0052] Where, represents the oil particle drift rate driven by wind stress, It is a comprehensive influencing factor of wind stress, which is used to correct the influence of wind stress on the drift rate of oil particles at different water depths. Indicates wind speed.

[0053] (9)

[0054] Where, represents the angle of the oil particle drift velocity vector driven by wind stress, is the wind vector angle, It represents the deflection angle of oil particle drift driven by wind stress relative to the wind vector at different water depths, and is used to correct the influence of wind stress on the angle of oil particle drift velocity vector at different water depths.

[0055] Step S3: Design an algorithm for the rising velocity of oil particles below the sea surface.

[0056] In a preferred example, the steps include:

[0057] (1) Combining the classical experimental parameter data of water viscosity coefficient, interpolation calculation is performed to obtain the water viscosity coefficient at different water temperatures;

[0058] (2) Based on the median diameter of oil particles, gravitational acceleration, oil particle density, water density and water viscosity, calculate the rising velocity of oil particles below the sea surface.

[0059] (10)

[0060] Where, represents the rising speed of oil particles below the sea surface, represents the acceleration due to gravity, represents the median diameter of oil particles, represents the density of oil particles, represents the water density, Represents the viscosity coefficient of water.

[0061] Step S4: Design an algorithm for calculating the vertical mixed layer height and vertical diffusion coefficient of the sea surface.

[0062] The height of the vertical mixed layer at the sea surface is mainly used to determine whether the vertical position of the oil particles belongs to the sea surface boundary mixing layer, which affects the calculation of whether the oil particles will carry out sea surface plane expansion, evaporation, emulsification and entrainment processes.

[0063] The vertical diffusion coefficient is mainly used to calculate the vertical diffusion velocity component of oil particles below the sea surface.

[0064] In a preferred example, the steps include:

[0065] (1) Calculate the vertical mixing layer height at the sea surface based on the sea surface wave height and sea surface wind speed;

[0066] (11)

[0067] Where, represents the vertical mixed layer height at the sea surface, represents the sea surface wind speed, represents the calculation time step, Indicates the height of breaking waves on the sea surface.

[0068] (2) Calculate the vertical diffusion coefficient at different water depths below the sea surface based on water depth, surface wave wavelength, effective wave period, and breaking wave height.

[0069] (12)

[0070] Where, represents the vertical diffusion coefficient at different water depths below the sea surface, Indicates the breaking wave height, represents the effective wave period, Indicates water depth. Represents the wavelength of sea surface waves.

[0071] Step S5: Combine the flow field forecast data and the wind field forecast data to drive the calculation of the three-dimensional oil spill drift and diffusion model.

[0072] In a preferred example, the steps include:

[0073] (1) Combining the flow field forecast data and the wind field forecast data, the actual drift velocity vector of the oil particles is calculated based on the three-dimensional drift velocity component of the oil particles driven by the ocean current, the velocity component driven by wind stress, the rising velocity component driven by seawater buoyancy, and the diffusion velocity component of the oil particles on and below the sea surface obtained in the first four steps;

[0074] (13)

[0075] Where, represents the actual drift velocity vector of the oil particles, represents the three-dimensional drift velocity component of oil particles driven by ocean currents, The velocity component of oil particles driven by wind stress, represents the rising velocity component of oil particles driven by the buoyancy of seawater, It represents the diffusion velocity component of oil particles on and below the sea surface.

[0076] (14)

[0077] Where, represents the diffusion velocity component of oil particles on and below the sea surface, is a random number between -1 and 1, indicating a random diffusion process. represents the diffusion coefficient vector of oil particles in horizontal and vertical directions, Indicates the calculation time step.

[0078] (2) The three-dimensional oil spill drift and diffusion model calculates the actual displacement of each oil particle at each time step based on the setting of the calculation time step and the actual drift speed of each oil particle, thereby realizing the simulation calculation of related physical processes such as the floating, drift and diffusion of oil particles in the three-dimensional water space from the seabed to the sea surface. Figure 2 The image shows the rise of oil particles from an oil leak starting at a point 15 meters below the seabed. The red dashed line indicates the depth of the mixed layer at the sea surface. Figure 3 What is shown is Figure 2 The horizontal plane in the example shows the results, with the blue parts corresponding to the oil particles.

[0079] Based on the same concept as the above method, an embodiment of the present application also provides a computer device, including a memory, a processor, and computer-executable instructions stored in the memory and executable on the processor. When the processor executes the instructions, it implements the above-mentioned method for rapid oil spill simulation applied to submarine pipeline oil leakage accidents.

[0080] Based on the same concept as the above method, an embodiment of the present application also provides a computer storage medium, which stores computer-executable instructions. When the computer-executable instructions are executed, the above-mentioned oil spill rapid simulation method applied to submarine pipeline oil leakage accidents is implemented.

[0081] Those skilled in the art will appreciate that all or part of the processes in the above-described method embodiments can be implemented by instructing the relevant hardware using a computer program. The computer program can be stored in a non-volatile computer-readable storage medium. When executed, the computer program can include the processes of the above-described method embodiments. Any reference to memory, storage, database, or other media used in the various embodiments provided herein may include non-volatile and / or volatile memory. Non-volatile memory may include read-only memory (ROM), programmable ROM (PROM), electrically programmable ROM (EPROM), electrically erasable programmable ROM (EEPROM), or flash memory. Volatile memory may include random access memory (RAM) or external cache memory. By way of illustration and not limitation, RAM is available in various forms, such as static RAM (SRAM), dynamic RAM (DRAM), synchronous DRAM (SDRAM), double data rate SDRAM (DDRSDRAM), enhanced SDRAM (ESDRAM), synchronous link DRAM (SLDRAM), RAMbus direct RAM (RDRAM), direct RAMbus dynamic RAM (DRDRAM), and RAMbus dynamic RAM (RDRAM).

[0082] In summary, this application is based on the particle tracer numerical calculation method. On the basis of the conventional two-dimensional oil spill model, the vertical physical change process of oil particles is added. By loading three-dimensional flow field data, the simulation calculation of related physical processes such as the floating, drifting and diffusion of oil particles in the three-dimensional water space from the seabed to the sea surface is realized. In this way, a rapid oil spill simulation method that can be applied to submarine pipeline oil leakage accidents is designed, which provides reference value for the rapid simulation and emergency deduction decision-making of petroleum pollutant leakage in marine pollution.

[0083] The technical solution of the present invention is not limited to the above-mentioned specific embodiments. Any technical variations made according to the technical solution of the present invention fall within the protection scope of the present invention.

Claims

1. A rapid oil spill simulation method for submarine pipeline oil leakage accidents, characterized in that The method comprises the following steps: Step 1: Read 3D flow and wind field data: Based on the two-dimensional oil spill model on the ocean surface, a module for reading three-dimensional flow field and wind field data is added; Read 3D flow field data, confirm the flow field file structure, and extract velocity components, water level and water depth information; Read wind field data, confirm the wind field file structure, and extract sea surface wind speed information; Step 2: Calculate the force and motion parameters of the oil particles: Based on the read three-dimensional flow field data, the drift velocity component of the oil particles driven by the ocean current in the three-dimensional water space is calculated; Based on the read wind field data, the velocity components and deflection angles of oil particles driven by wind stress at different water depths are calculated; Combining the physical properties of oil particles and water characteristics, the rising speed of oil particles below the sea surface is calculated; Calculate the vertical diffusion coefficient at different water depths below the sea surface based on water depth, surface wave wavelength, significant wave period, and breaking wave height; Step 3: Simulate the three-dimensional movement of oil particles: Combining 3D flow and wind field data with calculated oil particle forces and motion parameters to drive a 3D oil spill drift and diffusion model; Based on the calculation time step, the actual displacement of each oil particle at each time step is calculated to simulate the floating, drifting and diffusion process of oil particles from the seabed to the sea surface.

2. The method for rapid simulation of oil spills in submarine pipeline oil leaks according to claim 1, characterized in that: The three-dimensional flow field data and wind field data are read based on the NetCDF format, including flow velocity components, water level and water depth information, and sea surface wind speed information.

3. The rapid oil spill simulation method for submarine pipeline oil leakage accident according to claim 1 or 2, characterized in that: The velocity component of the oil particles driven by wind stress is calculated as follows: The surface driving factor of wind stress and the sub-deep driving factor of wind stress are determined by using the empirical formula regression curve; Combining the sea surface wind speed and the water depth where the oil particles are located, the velocity components of the oil particles driven by wind stress at different water depths are calculated.

4. The rapid oil spill simulation method for submarine pipeline oil leakage accidents according to claim 3 is characterized in that: It also includes determining the factors affecting the deflection angle driven by wind stress according to the water depth where the oil particles are located; The velocity component driven by wind stress and the deflection angle influencing factor are combined to calculate the drift velocity vector of the oil particles driven by wind stress at different water depths.

5. The rapid oil spill simulation method for submarine pipeline oil leakage accidents according to claim 1 is characterized in that: It also includes calculating the height of the vertical mixing layer on the sea surface, which is used to determine whether the vertical position of the oil particles belongs to the sea surface boundary mixing layer, and is obtained based on the wave height of the breaking waves and the sea surface wind speed.

6. The rapid oil spill simulation method for submarine pipeline oil leakage accidents according to claim 1 or 5, characterized in that: The three-dimensional oil spill drift and diffusion model calculates the actual displacement of each oil particle at each time step based on the calculation time step and the actual drift velocity vector of the oil particle.

7. The rapid oil spill simulation method for submarine pipeline oil leakage accidents according to claim 6 is characterized in that: The actual drift velocity vector of the oil particles includes a three-dimensional drift velocity component driven by ocean currents, a velocity component driven by wind stress, an upward velocity component driven by seawater buoyancy, and a diffusion velocity component of the oil particles on and below the sea surface.

8. The rapid oil spill simulation method for submarine pipeline oil leakage accidents according to claim 1 is characterized in that: The simulation also includes a visual display of the movement trajectory and distribution of oil particles in the three-dimensional water space.

9. A computer device comprising a memory, a processor, and computer executable instructions stored in the memory and executable on the processor, characterized in that: When the processor executes the instructions, the method for rapid oil spill simulation applied to submarine pipeline oil leakage accidents according to any one of claims 1 to 8 is implemented.

10. A computer storage medium storing computer executable instructions, characterized in that: When the computer executable instructions are executed, the oil spill rapid simulation method for submarine pipeline oil leakage accidents according to any one of claims 1 to 8 is implemented.

Images