CN-122023616-A - Offshore oil spill diffusion simulation method and device, storage medium and computer equipment

CN122023616ACN 122023616 ACN122023616 ACN 122023616ACN-122023616-A

Abstract

The invention relates to the technical field of marine environment science, and discloses a marine oil spill diffusion simulation method, a device, a storage medium and computer equipment, wherein the method comprises the steps of creating a two-dimensional image of a simulated sea area and drawing an initial oil film density image; the method comprises the steps of constructing an integrated speed field, executing image remapping operation to obtain a advection oil film density image, constructing an anisotropic Gaussian diffusion kernel according to wind field data, performing convolution operation with the advection oil film density image to obtain a diffusion oil film density image, executing morphological expansion operation to obtain an expansion oil film density image, performing mask operation based on a binary coastline mask image, separating a stranded oil film region and a floating oil film region, updating a stranded oil amount, repeatedly executing the steps until simulation is finished, and outputting a space-time evolution result of oil film diffusion. The method converts the computational fluid dynamics problem into the image processing problem, improves the computational efficiency and the system flexibility and reduces the economic cost on the premise of ensuring the simulation rationality of the physical process.

Inventors

LIU ZIJING
LIU SHENGQIANG
ZHANG LIGUO
CHENG JINXIANG
YAO HAIYUAN
ZHENG CHAOHUI
LV JINPENG
ZHANG GAOLING

Assignees

交通运输部规划研究院

Dates

Publication Date: 20260512
Application Date: 20251223

Claims (10)

1. An offshore oil spill diffusion simulation method, comprising: creating a two-dimensional image of a simulated sea area, and drawing an initial oil film density image at the current moment corresponding to an initial oil spilling area on the two-dimensional image; Constructing a comprehensive speed field based on flow field data and wind field data, and performing image remapping operation on the initial oil film density image based on the comprehensive speed field to obtain a advection oil film density image; constructing an anisotropic Gaussian diffusion kernel according to the wind field direction and the wind field intensity in the wind field data, and performing convolution operation on the anisotropic Gaussian diffusion kernel and the advection oil film density image to obtain a diffusion oil film density image; Performing morphological expansion operation on the diffusion oil film density image to obtain an expansion oil film density image, performing mask operation on the expansion oil film density image based on a preset binary coastline mask image, separating a stranded oil film region and a floating oil film region, and updating the stranded oil amount corresponding to the stranded oil film region in the coastline mask image; And taking the floating oil film area as an initial oil film density image at the next simulation moment, repeatedly executing the steps until the simulation is finished, and outputting a space-time evolution result of oil film diffusion.
2. The method of claim 1, wherein creating a two-dimensional image of the simulated sea area and drawing an initial oil film density image of the current time corresponding to the initial oil spill area on the two-dimensional image comprises: creating a two-dimensional image corresponding to the simulated sea area, and setting all pixel values in the two-dimensional image as background values; drawing a solid circular area on the two-dimensional image based on the position coordinates and the radius of the influence range of the oil spill accident, and setting the pixel value in the solid circular area as an oil film value; Calculating the sum of pixel values in the solid circular area, obtaining the total oil leakage amount of the oil overflow accident, and calculating the ratio of the total oil leakage amount to the sum of pixel values to obtain a normalization factor; and multiplying each pixel value in the two-dimensional image by the normalization factor to obtain an initial oil film density image at the current moment.
3. The method of claim 1, wherein constructing a composite velocity field based on flow field data and wind field data and performing an image remapping operation on the initial oil film density image based on the composite velocity field to obtain a advection oil film density image comprises: Linearly superposing the flow field data and the wind field data based on a preset weight factor to construct a comprehensive speed field, wherein the value range of the weight factor of the wind field data is 0.01 to 0.05; Calculating a source coordinate corresponding to each pixel in the initial oil film density image at the last moment by adopting a reverse particle tracking strategy based on the comprehensive speed field and a preset time step; Performing image remapping sampling operation on the initial oil film density image according to the source coordinates to obtain a first intermediate oil film density image after advection; Performing point multiplication operation on the first intermediate oil film density image and a preset binary coastline mask image to obtain a second intermediate oil film density image after advection, and calculating to obtain stranded oil quantity according to an operation result of the point multiplication operation; and carrying out normalized correction on the total mass of the second intermediate oil film density image and the stranded oil mass based on the total mass of the initial oil film density image, and taking the corrected second intermediate oil film density image as a final advection oil film density image.
4. The method of claim 1, wherein constructing an anisotropic gaussian diffusion kernel from the wind field direction and the wind field intensity in the wind field data, and convolving the anisotropic gaussian diffusion kernel with the advection oil film density image, to obtain a diffusion oil film density image, comprises: Determining a diffusion parameter matrix for representing diffusion anisotropy according to the wind field direction and the wind field intensity in the wind field data and a preset basic diffusion coefficient, wherein the diffusion intensity in the downwind direction is larger than that in the crosswind direction; Generating an anisotropic Gaussian convolution kernel based on the diffusion parameter matrix, and carrying out normalization processing on the anisotropic Gaussian convolution kernel, wherein the direction of the anisotropic Gaussian convolution kernel is aligned with the direction of the wind field; Performing two-dimensional convolution operation on the advection oil film density image and the anisotropic Gaussian convolution kernel after normalization processing to obtain an initial diffusion result image; Performing dot multiplication operation on the initial diffusion result image and a preset binary coastline mask image to obtain a first intermediate diffusion image; And carrying out normalization correction on the first intermediate diffusion image based on the total mass of the advection oil film density image to obtain a final diffusion oil film density image.
5. The method of claim 1, wherein performing a morphological dilation operation on the diffuse oil film density image results in an extended oil film density image, comprising: Generating a spreading core for executing morphological expansion operation according to a wind field direction, a wind field intensity and a preset basic diffusion coefficient in the wind field data, wherein when the wind field intensity exceeds a preset threshold value, the spreading core is elliptical, the major axis direction of the spreading core is aligned with the wind field direction, the major axis radius of the ellipse increases along with the increase of the wind field intensity, and when the wind field intensity is lower than the preset threshold value, the spreading core is circular; Performing morphological expansion operation on the diffusion oil film density image and the spreading core to obtain an initial expansion result image; and calculating the ratio of the total mass of the diffusion oil film density image to the total mass of the initial expansion result image to obtain a correction factor, and multiplying the initial expansion result image by the correction factor to obtain an expansion oil film density image.
6. The method of claim 1, wherein masking the extended oil film density image based on a predetermined binary coastline mask image to separate a stranded oil film region and a floating oil film region, and updating a stranded oil amount corresponding to the stranded oil film region in the coastline mask image comprises: performing dot multiplication operation on the extended oil film density image and a preset binary coastline mask image to obtain a potential stranded oil mass distribution image; Multiplying the potential stranded oil mass distribution image by a preset stranded probability parameter to obtain a newly increased actual stranded oil mass image at the current moment, determining a stranded oil film area based on the newly increased actual stranded oil mass image, and multiplying the potential stranded oil mass distribution image by a non-stranded probability parameter to obtain an oil mass image with the possibility of backflow, wherein the sum of the stranded probability parameter and the non-stranded probability parameter is 1; Subtracting the potential stranded oil mass distribution image from the extended oil film density image to obtain a first intermediate image, performing point multiplication operation on an oil mass image which is likely to flow back and the complement of the binary coastline mask image to obtain a backflow oil mass image which is redistributed to an adjacent water area, adding the first intermediate image and the backflow oil mass image to obtain an updated floating oil film image, and determining a floating oil film area based on the floating oil film image; and accumulating the newly increased actual stranded oil mass image at the current moment in the accumulated stranded oil mass image for recording the total stranded oil mass, and updating the stranded oil mass corresponding to the stranded oil film region.
7. The method of claim 1, wherein after the separating out of stranded oil film region and floating oil film region, updating a stranded oil amount corresponding to the stranded oil film region in the coastline mask image, the method further comprises: Determining the evaporation rate and the emulsification rate of the current time step based on environmental parameters and oil properties, wherein the environmental parameters comprise temperature, wind speed and wave height, and the oil properties comprise oil types; according to the evaporation rate, the emulsification rate and a preset time step, calculating to obtain an evaporation loss coefficient and an emulsification loss coefficient; Multiplying the floating oil film image corresponding to the floating oil film region by the evaporation loss coefficient and the emulsification loss coefficient simultaneously to obtain an oil film density image after quality loss; And calculating the mass loss of the floating oil film image caused by evaporation and emulsification, and accumulating the mass loss into the total evaporation amount and the total emulsification amount respectively.
8. An offshore oil spill diffusion simulation device, comprising: The image initialization module is used for creating a two-dimensional image of the simulated sea area and drawing an initial oil film density image at the current moment corresponding to the initial oil spilling area on the two-dimensional image; The advection migration module is used for constructing a comprehensive speed field based on flow field data and wind field data, and executing image remapping operation on the initial oil film density image based on the comprehensive speed field to obtain an advection oil film density image; The turbulence diffusion module is used for constructing an anisotropic Gaussian diffusion core according to the wind field direction and the wind field intensity in the wind field data, and carrying out convolution operation on the anisotropic Gaussian diffusion core and the advection oil film density image to obtain a diffusion oil film density image; The shoreline interaction module is used for performing morphological expansion operation on the diffusion oil film density image to obtain an expansion oil film density image, performing mask operation on the expansion oil film density image based on a preset binary shoreline mask image, separating a stranded oil film region and a floating oil film region, and updating the stranded oil amount corresponding to the stranded oil film region in the shoreline mask image; and the iteration output module is used for taking the floating oil film area as an initial oil film density image at the next simulation moment, repeatedly executing the steps until the simulation is finished, and outputting a space-time evolution result of oil film diffusion.
9. A storage medium having stored thereon a computer program, which when executed by a processor, implements the steps of the method of any of claims 1 to 7.
10. A computer device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, characterized in that the computer program when executed by the processor implements the steps of the method according to any one of claims 1 to 7.

Description

Offshore oil spill diffusion simulation method and device, storage medium and computer equipment Technical Field The invention relates to the technical field of marine environment science, in particular to a marine oil spill diffusion simulation method, a device, a storage medium and computer equipment. Background The offshore oil spill accident has the characteristics of strong burst, high diffusion speed, profound ecological influence and the like, and the rapid and accurate prediction of the migration path and coverage range of the oil film is a precondition for carrying out emergency treatment, resource scheduling and ecological damage evaluation; currently, the mainstream oil spill diffusion simulation technology is based on Computational Fluid Dynamics (CFD) framework, discretizes the ocean area into a structured or unstructured grid, simulates the space-time evolution of an oil film concentration field by solving a convection-diffusion equation set containing advection, turbulence diffusion, wind drag, evaporation, emulsification and other physical and chemical processes, and the method is usually integrated in a professional numerical simulation platform such as MIKE, GNOME, OSISAF, relies on high-precision ocean and meteorological forced field input, and performs multi-step iterative computation under physical constraints such as conservation of quality. However, the conventional technology faces serious challenges in practical application, firstly, the calculation complexity is high, single simulation often needs 6 to 12 hours or more, real-time requirements on prediction results in a golden window period of oil spill emergency response are difficult to meet even in a high-performance computing environment, secondly, the system highly depends on commercial software authorization and special hardware resources, the purchase and operation costs are high, a core algorithm is closed, first-line departments such as local maritime affairs and environmental protection are difficult to bear or autonomously master, and more importantly, the existing model architecture lacks flexibility, is difficult to efficiently support complex tasks such as multi-source leakage, batch scene deduction or man-machine interaction parameter adjustment, and is weak in regional risk assessment and dynamic decision support scenes. Therefore, an offshore oil spill diffusion simulation method with strong timeliness, low cost, strong controllability and strong flexibility is urgently needed to be researched. Disclosure of Invention In view of the above, the application provides a method, a device, a storage medium and a computer device for simulating oil spill diffusion at sea, which mainly aims to solve the technical problems of poor timeliness, high cost, poor controllability and poor flexibility of the oil spill diffusion simulation technology in the prior art. According to a first aspect of the present invention, there is provided an offshore oil spill diffusion simulation method comprising: creating a two-dimensional image of a simulated sea area, and drawing an initial oil film density image at the current moment corresponding to an initial oil spilling area on the two-dimensional image; Constructing a comprehensive speed field based on flow field data and wind field data, and performing image remapping operation on the initial oil film density image based on the comprehensive speed field to obtain a advection oil film density image; constructing an anisotropic Gaussian diffusion kernel according to the wind field direction and the wind field intensity in the wind field data, and performing convolution operation on the anisotropic Gaussian diffusion kernel and the advection oil film density image to obtain a diffusion oil film density image; Performing morphological expansion operation on the diffusion oil film density image to obtain an expansion oil film density image, performing mask operation on the expansion oil film density image based on a preset binary coastline mask image, separating a stranded oil film region and a floating oil film region, and updating the stranded oil amount corresponding to the stranded oil film region in the coastline mask image; And taking the floating oil film area as an initial oil film density image at the next simulation moment, repeatedly executing the steps until the simulation is finished, and outputting a space-time evolution result of oil film diffusion. Optionally, the method comprises the steps of creating a two-dimensional image of a simulated sea area, drawing an initial oil film density image of the current moment corresponding to an initial oil spilling area on the two-dimensional image, setting all pixel values in the two-dimensional image as background values, drawing a solid circular area on the two-dimensional image based on position coordinates and an influence range radius of an oil spilling accident, setting the pixel values in the solid circular area as oil film values,