CN111368414B - Offshore oil spill track simulation method based on thickness information - Google Patents

Offshore oil spill track simulation method based on thickness information Download PDF

Info

Publication number
CN111368414B
CN111368414B CN202010127439.6A CN202010127439A CN111368414B CN 111368414 B CN111368414 B CN 111368414B CN 202010127439 A CN202010127439 A CN 202010127439A CN 111368414 B CN111368414 B CN 111368414B
Authority
CN
China
Prior art keywords
oil
volume
spill
water
water surface
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202010127439.6A
Other languages
Chinese (zh)
Other versions
CN111368414A (en
Inventor
陈曦
金文�
韩洪洪
王万财
李彦知
郭海楠
韩雅桐
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tianjin University
Original Assignee
Tianjin University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Tianjin University filed Critical Tianjin University
Priority to CN202010127439.6A priority Critical patent/CN111368414B/en
Publication of CN111368414A publication Critical patent/CN111368414A/en
Application granted granted Critical
Publication of CN111368414B publication Critical patent/CN111368414B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention discloses a marine oil spill track simulation method based on thickness information. The method and the device have no special requirements on the oil film thickness acquisition device and method, have wide application range, consider the realizable problem in computer operation, and have strong practical value and wide market prospect.

Description

Offshore oil spill track simulation method based on thickness information
Technical Field
The invention relates to simulation of an oil spill drift track after an offshore oil spill accident, in particular to a method for simulating an oil spill track by comprehensively utilizing oil film thickness information and an oil particle method acquired by an offshore oil film thickness measuring device.
Background
Offshore oil spill accidents are the most serious marine ecological pollution and ecological disaster, and the economic loss and pollution consequences caused by the marine oil spill accidents are huge. The timely and efficient emergency treatment of the oil spill pollution accident is an effective means for reducing the influence of the oil spill disaster. The oil spilling track simulation aims at simulating the irregular shape, the drifting track and the oil spilling amount of an oil film, and plays an essential important role in decision and implementation of an emergency treatment scheme of an oil spilling accident.
At present, an oil particle method is mainly adopted for simulating an oil spilling track, the method tracks the process of the oil particle micro-cluster moving along with the advection and turbulence diffusion of the surrounding water body, and the influence of various ocean power factors on the oil spilling diffusion can be accurately described. The oil particle method adopts iterative computation, the accuracy of subsequent iterative computation is directly determined by setting the initial value of the computation, and the closer the initial value is to the true value, the more accurate the simulation is, and the slower the error divergence rate is. The oil particle method takes the oil spill area and the oil spill thickness as initial values for calculation. In acquiring the oil spill area, a remote sensing image is generally used, and since the remote sensing image is one of the accepted methods for acquiring the best oil spill area, the oil spill area, which is one of the initial values of calculation, has high accuracy. However, the existing oil particle method estimates the oil spill thickness by using an empirical value, and has a large deviation with the actual oil spill thickness, and the thickness empirical value is used as a calculation initial value to carry out iterative operation of marine oil spill trajectory simulation, so that the simulation precision is reduced, and the oil spill trajectory prediction time meeting the precision requirement is shortened.
Disclosure of Invention
The invention aims to provide a method for simulating an offshore oil spill track based on thickness information aiming at the technical defects in the prior art, and the method is used for simulating the oil spill track by using oil film thickness information and an oil particle method acquired by an offshore oil film thickness measuring device.
The technical scheme adopted for realizing the purpose of the invention is as follows:
a marine oil spill track simulation method based on thickness information comprises the following steps:
s101, obtaining initial values of oil spilling areas of thin oil film and thick oil film by using remote sensing images
Figure BDA0002394830050000021
And
Figure BDA0002394830050000022
obtaining initial values of the thicknesses of the thin oil film and the thick oil film spilled oil by using a marine oil film thickness measuring device
Figure BDA0002394830050000023
And
Figure BDA0002394830050000024
calculating initial values for thin oil film, thick oil film, and total spill volume
Figure BDA0002394830050000025
Figure BDA0002394830050000026
And VS(t0) (ii) a And calculating the initial volume value v (n) of the oil particlesk,t0)，
Wherein,
Figure BDA0002394830050000027
Figure BDA0002394830050000028
wherein N is the total number of oil particles, NSThe number of oil particles on the water surface;
s102, updating the oil spill volume by using a weathering model, and calculating the oil spill volume on the water surface at the time t;
according to a calculation formula of the volume change rate of the spilled oil obtained by a differential equation, and by combining various initial values obtained in the step S101, the volume of the thin oil film on the water surface at the time t is obtained through calculation
Figure BDA0002394830050000029
Volume of thick oil film on water surface
Figure BDA00023948300500000210
And total volume V of oil spilled on the waterS(t)；
Wherein, the calculation formula of the volume change rate of the spilled oil is as follows;
Figure BDA00023948300500000211
in the formula,
Figure BDA00023948300500000212
is a differential symbol, TN represents a thin oil film, TK represents a thick oil film,
Figure BDA00023948300500000213
is the volume of the thin oil film on the water surface,
Figure BDA00023948300500000214
is the volume of the thick oil film on the water surface, VSIs the total volume of oil spills on the water surface;
s103, calculating the volume of the updated oil particles by using the initial value of the oil spill volume and the updated oil spill volume;
the coupling equation (13) uses the initial volume value of the oil particles obtained in step S101 and the volume of the thin oil film on the water surface at time t calculated in step S102
Figure BDA00023948300500000215
Volume of thick oil film on water surface
Figure BDA00023948300500000216
Calculating to obtain the volume v (n) of oil particles at the time tkT), calculating the volume v (n) of the oil particles at time tkAnd t) calculation:
Figure BDA00023948300500000217
in the formula,
Figure BDA0002394830050000031
the percentage of non-evaporable oil particles in the whole oil particles, NE is the number of oil particles affected by evaporation,
Figure BDA0002394830050000032
the vertical line (E) represents the evaporation;
s104, updating coordinate information of the oil particles by using the hydrodynamic model;
calculating to obtain coordinate information of the oil particles at the time t by calculating position increment of the oil particles on the water surface and dissolved in the water, and performing integrated operation on the coordinate information to realize simulation of an oil spill track;
and S105, calculating the oil spilling concentration by using a discretized numerical calculation method, and multiplying the calculated oil spilling concentration by the calculated oil spilling volume at the time t to obtain the oil spilling amount.
Wherein, among others,
Figure BDA0002394830050000033
and
Figure BDA0002394830050000034
the calculation method of (c) is as follows:
Figure BDA0002394830050000035
Figure BDA0002394830050000036
Figure BDA0002394830050000037
indicating the rate of change in volume resulting from evaporation,
Figure BDA0002394830050000038
indicating the rate of change of the volume of spill oil dissolved in water,
Figure BDA0002394830050000039
indicating the rate of change of the volume of spill oil on the water surface.
The method for calculating the position increment of the oil particles on the water surface and dissolved in the water comprises the following steps:
wherein the position increment of the oil particles on the water surface is calculated as follows:
dxkS(t)＝[UC(xk,yk,0,t)+UW(xk,yk,t)+US(xk,yk,t)dt+dx＇k(t)，
dxkS(t) represents the increase in the position of the oil particles on the water surface, UCIs the euler velocity component;
UW＝(UW,VW) Is the effect of ocean currents, related to the strength of sea surface wind and the angle between wind and ocean currents, Uw,VwIs the force of the ocean current along two perpendicular components;
Uw＝α(Wxcosβ+Wysinβ)，
Vw＝α(-Wxsinβ+Wycosβ)
Wxand WyRespectively, the intensity of the wind at the predetermined location in the latitude and longitude directions, alpha is the drift factor, beta is the drift angle,
US＝(US,VS) Is the wave action, the calculation method is shown as the following formula:
Figure BDA0002394830050000041
US＝DScosθ，Vs＝DSsinθ，
wherein ω is angular frequency, S is wave spectrum, k is wave number;
dx′k(t) is the random portion of the oil particle motion,
Figure BDA0002394830050000042
wherein r is a random number between 0 and 1, khIs a constant;
the position increment of oil particles dissolved in water is calculated as follows:
dxkD(t)＝UC(xk,yk,zk,t)dt+dx′k(t)，
dxkD(t) represents the increase in the position of oil particles dissolved in water.
The calculation method for calculating the oil spill concentration by using the discretization numerical calculation method comprises the following steps:
Figure BDA0002394830050000043
Figure BDA0002394830050000044
CS(xT,yTt) represents the concentration of spilled oil on the water surface, CD(xT,yTT) represents the concentration of spilled oil on the surface of the water dissolved in the concentration of spilled oil in the water, nSRepresenting the total number of oil particles on the water surface, nDRepresents the total number of oil particles dissolved in water, ρ represents the oil spill density, δ xTδyTThe grid length in the longitude and latitude directions when grid discretization is performed in grid calculation on a computer is shown.
According to the method, the thickness information of the offshore spilled oil film is combined with the existing oil particle method to simulate the offshore spilled oil track, the thickness of the oil film obtained by the offshore oil film thickness measuring device is used for replacing the thickness experience estimation value used by the original oil particle method, the defects of the existing oil particle method can be effectively overcome, the precision of the spilled oil track simulation is improved, the error diffusion rate is reduced, and the prediction time is prolonged.
The method and the device have no special requirements on the oil film thickness acquisition device and method, have wide application range, consider the realizable problem in computer operation, and have strong practical value and wide market prospect.
Drawings
FIG. 1 is a flow chart of a marine oil spill trajectory simulation method based on thickness information;
FIG. 2 is a schematic view of a processing apparatus according to an embodiment;
FIG. 3 is a schematic diagram showing the marine oil spill trajectory 8 hours after the simulated oil spill of the embodiment occurs;
FIG. 4 is a schematic diagram showing the marine oil spill trajectory 16 hours after the simulated oil spill of the embodiment occurs;
FIG. 5 is a schematic diagram showing the marine oil spill trajectory 24 hours after the simulated oil spill of the example.
Detailed Description
The invention is described in further detail below with reference to the figures and specific examples. It should be understood that the specific embodiments described herein are merely illustrative of the invention and are not intended to limit the invention.
As shown in FIG. 1, the method for simulating the marine oil spill track based on the thickness information comprises the following steps:
the method comprises the steps of firstly obtaining a calculation initial value of the oil spilling volume by using a remote sensing image and the marine oil spilling thickness, then updating the oil spilling volume by using a weathering model, then updating the oil particle volume by using the oil spilling volume and an updating formula, then updating coordinate information of oil particles by using a hydrodynamic model so as to obtain an oil spilling track, and finally calculating the oil spilling concentration and the oil spilling amount by using a discretization numerical calculation method.
As shown in fig. 2, this fig. 2 shows a simulation apparatus used for the simulation method of the present invention, including a satellite 1 for acquiring a remote sensing image of the area of a spill oil 3 on the water surface, an offshore spill oil thickness measurement buoy 2 for acquiring the thickness of the spill oil on the water surface; the lower part of the oil spill 3 on the water surface is provided with oil spill 4 dissolved in the water.
The following detailed calculation procedures for the simulation method of the present invention are described below:
and S101, acquiring a calculation initial value of the oil spill volume by using the remote sensing image and the marine oil spill thickness.
The method comprises the following steps: acquiring the oil spill area by using a satellite remote sensing image, acquiring the oil spill thickness by using a marine oil spill thickness measuring buoy, and calculating the initial value of the volume of the thin oil film on the water surface by using the formulas (10), (11) and (12)
Figure BDA0002394830050000051
Volume of thick oil film on water surface
Figure BDA0002394830050000052
And total volume V of oil spilled on the waterS(t0)。
In order to realize the invention, the following processes are firstly carried out:
11) and establishing a continuous time dynamic equation of the oil spill concentration in the marine environment.
Because the spilled oil in the marine environment can be acted by ocean currents to generate movement, the spilled oil can generate physical and chemical changes in the movement process. Thus, a continuous time dynamic equation for spill oil concentration in a marine environment is established as follows:
Figure BDA0002394830050000061
where C (x, y, z, t) is the oil spill concentration, x spatial coordinates are (x, y, z), t is time, U is the sea current field, consisting of elements (U, V, W), and K is the diffusion tensor, which is used to describe the turbulent effect. r isj(x, C (x, t), t) is used for describing the concentration change rate caused by the physical and chemical changes of the spilled oil, j represents the jth factor influencing the concentration change, U is a sea flow field, three-dimensional fields are generally represented by components under three orthogonal coordinates such as (U, V, W), the formula (1) is decomposed into a formula (2) and a formula (3) to enable the left side of the formula (1) to be subjected to decomposition
Figure BDA0002394830050000066
Moving to the right of the equal sign to obtain the product.
Figure BDA0002394830050000062
Figure BDA0002394830050000063
In the formula, C1Is the concentration of spilled oil under weathering, the factors that can cause concentration changes are many, and the weathering is mainly considered here, so C is used1Indicating the concentration of oil spill under weathering.
The invention does not consider the adsorption effect of the shoreline, so the formula (2) can decompose C into the concentration C of the spilled oil on the water surfaceSAnd the concentration of spilled oil C dissolved in waterD。CSIs calculated as shown in equation (4):
Figure BDA0002394830050000064
where ρ is the oil spill density, VSIs the volume of oil spilled on the water surface, ASIs the oil spill area on the water surface.
Concentration of spilled oil C dissolved in WaterDThe formula (5) is shown below:
Figure BDA0002394830050000065
where ρ is the oil spill density, VDIs the volume of oil spilled, A, dissolved in waterDIs the oil spill area dissolved in water.
Since the efflorescence changes the oil spill concentration with time, equations (4) and (5) are derived from time to give equations (6) and (7), respectively, for later use in grid computing in a computer.
Figure BDA0002394830050000071
Figure BDA0002394830050000072
To calculate equation (3), the volume V of oil spilled on the water surface is calculatedsDecomposed into N volumes of v (N)kT) oil particles of which each has a three-dimensional coordinate vector xk(nkT), a total of N oil particles, where N is assumedkThe oil particles have the same volume v (n) at time tk,t)。
xk(nk,t)＝(xk(nk,t),yk(nk,t),zk(nk,t)) (8)
22) Numerical solution of the continuous time equation of dynamics of the concentration of spilled oil in marine environments.
To solve the above continuous-time dynamical equations on a computer, discretization and numerical solution were performed. Firstly, the volume V of the oil spill on the water surfaceSIs decomposed into
Figure BDA0002394830050000073
And
Figure BDA0002394830050000074
as shown in formula (9).
Figure BDA0002394830050000075
In the formula, TN represents a thin oil film, TK represents a thick oil film,
Figure BDA0002394830050000076
is the volume of the thin oil film on the water surface,
Figure BDA0002394830050000077
is the volume of the thick oil film on the water surface.
According to the invention, the calculation of the initial value is to obtain the initial values of the oil spilling areas of the thin oil film and the thick oil film by using the remote sensing image
Figure BDA0002394830050000078
And
Figure BDA0002394830050000079
obtaining initial values of the thicknesses of the thin oil film and the thick oil film spilled oil by using a marine oil film thickness measuring device
Figure BDA00023948300500000710
And
Figure BDA00023948300500000711
then, initial values of the thin oil film, the thick oil film, and the total oil spill volume are calculated according to the following equations (10), (11), (12)
Figure BDA00023948300500000712
And VS(t0)
Figure BDA00023948300500000713
Figure BDA00023948300500000714
Figure BDA00023948300500000715
Initial value v (n) of oil particle volumek,t0) Can be calculated by the following formula (13)
Figure BDA00023948300500000716
Wherein N is the total number of oil particles, NSThe number of oil particles on the water surface.
And S102, updating the oil spill volume by using the weathering model.
The specific method comprises the following steps: obtaining the volume of the thin oil film on the water surface at the time t by numerical calculation according to the change rate of the volume of the spilled oil and the initial value obtained in the step S101 by using differential equations of the following formulas (15) and (16)
Figure BDA0002394830050000081
Volume of thick oil film on water surface
Figure BDA0002394830050000082
And total volume V of oil spilled on the waterS(t)。
Wherein the oil spill volume change rate can be obtained by differentiating the oil spill volume change rate by the formula (9)
Figure BDA0002394830050000083
Figure BDA0002394830050000084
Is a differential sign, due to
Figure BDA0002394830050000085
Not good for calculation, the oil film is divided into a thick oil film and a thin oil film, and the calculation is performed separately from the following equations (15) and (16)
Figure BDA0002394830050000086
And
Figure BDA0002394830050000087
Figure BDA0002394830050000088
Figure BDA0002394830050000089
Figure BDA00023948300500000810
indicating the rate of change in volume resulting from evaporation,
Figure BDA00023948300500000811
indicating the rate of change of the volume of spill oil dissolved in water,
Figure BDA00023948300500000812
indicating the rate of change of the volume of spill oil on the water surface.
And S103, updating the volume of the oil particles by using the oil spilling volume and an updating formula.
The method comprises the following steps: the initial volume value of the oil particles obtained by combining the following equations (17) and (18) with the equation (13), and the volume of the thin oil film on the water surface at time t calculated in S102
Figure BDA00023948300500000813
Volume of thick oil film on water surface
Figure BDA00023948300500000814
Is calculated toVolume of oil particles v (n) to time tk,t)。
Wherein the volume v (n) of oil particles at time tkT) is calculated as follows:
Figure BDA00023948300500000815
in the formula,
Figure BDA00023948300500000816
percentage of non-evaporable oil particles in the whole oil particles, NE representing the number of oil particles subjected to evaporation, f(E)(t) calculation according to equation (18)
Figure BDA00023948300500000817
The vertical line (E) shows the evaporation.
And S104, updating the coordinate information of the oil particles by using the hydrodynamic model.
The specific method comprises the following steps: the coordinate information of the oil particles at time t is calculated by equations (20), (21), (22) and (23), and the oil spill trajectory can be simulated by performing a collective operation on the coordinate information.
For oil particles on the water surface, the position increment dx thereofkS(t) is calculated as in formula (19)
dxkS(t)＝[UC(xk,yk,0,t)+UW(xk,yk,t)+US(xk,yk,t)dt+dx′k(t) (19)
UCIs the Euler velocity component, UW＝(UW,VW) Is the ocean current effect, which is related to the strength of sea surface wind and the included angle between wind and ocean current, and is calculated as shown in formula (20), UwAnd VwIs the force of the ocean current along two perpendicular components.
UW＝α(Wxcosβ+Wysinβ)
VW＝α(-Wxsinβ+Wycosβ) (20)
WxAnd WyThe intensity of the wind in the latitude and longitude directions at 10 meters, respectively, is the drift factor, and β is the drift angle.
US＝(US,VS) Is the action of sea waves, and the calculation method is shown as a formula (21)
Figure BDA0002394830050000091
US＝DScosθ
VS＝DSsinθ (21)
Where ω is the angular frequency, S is the wave spectrum, and k is the number of waves.
In formula (19), dx'k(t) is the random part of the oil particle motion, and the calculation method is shown in formula (22):
Figure BDA0002394830050000092
wherein r is a random number between 0 and 1, khIs a constant.
For oil particles dissolved in water, the position increment dx thereofkD(t) is calculated as shown in equation (23):
dxkD(t)＝UC(xk,yk,zk,t)dt+dx′k(t) (23)
and S105, calculating the oil spilling concentration and the oil spilling amount by using a discretized numerical calculation method.
The specific method comprises the following steps: the oil spill concentration can be calculated by the equations (26), (27), (28), and (29), and the oil spill amount can be obtained by multiplying the oil spill concentration by the oil spill volume at time t calculated in step S102.
Wherein, for calculation on a computer, the formulas (6) and (7) are subjected to grid discretization, and the length of each grid side in the longitude direction and the latitude direction is delta xTAnd δ yTLet AS＝AD＝δxTδyTThen the formulas (6) and (7) can representComprises the following steps:
Figure BDA0002394830050000101
Figure BDA0002394830050000102
the above calculation formulas are all continuous time equations, and a computer cannot directly calculate the continuous time equations and can calculate the continuous time equations only after discretization.
The oil particles on the water surface satisfying the formula (26) are counted to obtain the total number nS
Figure BDA0002394830050000103
The total number n was obtained by counting the oil particles dissolved in water satisfying the formula (27)D
Figure BDA0002394830050000104
According to formulae (24) and (25) and nSAnd nDThe concentration of oil spill in the available grids is shown as equations (28) and (29):
Figure BDA0002394830050000105
Figure BDA0002394830050000106
examples simulation of the marine oil spill trajectories after the oil spill occurs for 8 hours, 16 hours and 24 hours as shown in fig. 3, 4 and 5, in which the thick line regions are the oil spill trajectories and the arrow directions are the directions of the ocean current fields.
The foregoing is only a preferred embodiment of the present invention, and it should be noted that, for those skilled in the art, various modifications and decorations can be made without departing from the principle of the present invention, and these modifications and decorations should also be regarded as the protection scope of the present invention.

Claims (3)

1. A marine oil spill track simulation method based on thickness information is characterized by comprising the following steps:
s101, obtaining initial values of oil spilling areas of thin oil film and thick oil film by using remote sensing images
Figure FDA0003458250770000011
And
Figure FDA0003458250770000012
obtaining initial values of the thicknesses of the thin oil film and the thick oil film spilled oil by using a marine oil film thickness measuring device
Figure FDA0003458250770000013
And
Figure FDA0003458250770000014
calculating initial values for thin oil film, thick oil film, and total spill volume
Figure FDA0003458250770000015
Figure FDA0003458250770000016
And VS(t0) (ii) a And calculating the initial volume value v (n) of the oil particlesk，t0)，
Wherein,
Figure FDA0003458250770000017
Figure FDA0003458250770000018
wherein N is the total number of oil particles, NSThe number of oil particles on the water surface;
s102, updating the oil spill volume by using a weathering model, and calculating the oil spill volume on the water surface at the time t;
the weathering model is as follows:
Figure FDA0003458250770000019
Figure FDA00034582507700000110
wherein C (x, y, z, t) is the concentration of oil spill, x space coordinates are (x, y, z), t is time, U is the sea current field and is composed of elements (U, V, W), and K is the diffusion tensor and is used for describing the turbulent flow effect; c1Is the concentration of oil spill under weathering,
rj (x, C (x, t), t) is used to describe the concentration change rate caused by the physical and chemical changes of the oil spill;
decomposing C (x, y, z, t) into the concentration C of oil spilled on the water surfaceSAnd the concentration of spilled oil C dissolved in waterD；
Concentration of spilled oil C on the waterSThe calculation of (d) is as follows:
Figure FDA00034582507700000111
where ρ is the oil spill density, VSIs the volume of oil spilled on the water surface, ASIs the oil spill area on the water surface;
concentration of spilled oil C dissolved in WaterDThe calculation formula of (a) is as follows:
Figure FDA00034582507700000112
in the formula, VDIs the volume of oil spilled, A, dissolved in waterDIs the oil spill area dissolved in water;
since the concentration of the oil spill changes with time due to weathering, the time is derived to obtain the following formula for later grid calculation;
Figure FDA0003458250770000021
Figure FDA0003458250770000022
the volume V of the oil spill on the water surfacesDecomposed into N volumes of v (N)kT) oil particles of which each has a three-dimensional coordinate vector xk(nkT), a total of N oil particles, where N is assumedkThe oil particles have the same volume v (n) at time tk，t)；
xk(nk，t)＝(xk(nk，t)，yk(nk，t)，zk(nk，t))
The volume V of the oil spill on the water surfaceSIs decomposed into
Figure FDA0003458250770000023
And
Figure FDA0003458250770000024
the following were used:
Figure FDA0003458250770000025
in the formula, TN represents a thin oil film, TK represents a thick oil film,
Figure FDA0003458250770000026
is the volume of the thin oil film on the water surface,
Figure FDA0003458250770000027
is the water surfaceThe volume of the thick oil film above;
according to a calculation formula of the volume change rate of the spilled oil obtained by a differential equation, and by combining various initial values obtained in the step S101, the volume of the thin oil film on the water surface at the time t is obtained through calculation
Figure FDA0003458250770000028
Volume of thick oil film on water surface
Figure FDA0003458250770000029
And total volume V of oil spilled on the waterS(t)；
Wherein, the calculation formula of the volume change rate of the spilled oil is as follows;
Figure FDA00034582507700000210
in the formula,
Figure FDA00034582507700000211
is a differential sign;
s103, calculating the volume of the updated oil particles by using the initial value of the oil spill volume and the updated oil spill volume;
using the initial volume value of the oil particles obtained in step S101 and the volume of the thin oil film on the water surface at time t calculated in step S102
Figure FDA00034582507700000212
Volume of thick oil film on water surface
Figure FDA00034582507700000213
Calculating to obtain the volume v (n) of oil particles at the time tkT), calculating the volume v (n) of the oil particles at time tkAnd t) calculation:
Figure FDA0003458250770000031
in the formula,
Figure FDA0003458250770000032
the percentage of non-evaporable oil particles in the whole oil particles, NE is the number of oil particles affected by evaporation,
Figure FDA0003458250770000033
the vertical line (E) represents the evaporation;
s104, updating coordinate information of the oil particles by using the hydrodynamic model;
calculating to obtain coordinate information of the oil particles at the time t by calculating position increment of the oil particles on the water surface and dissolved in the water, and performing integrated operation on the coordinate information to realize simulation of an oil spill track;
s105, calculating the oil spilling concentration by using a discretized numerical calculation method, and multiplying the calculated oil spilling concentration by the calculated oil spilling volume at the time t to obtain the oil spilling amount;
the calculation method for calculating the oil spill concentration by using the discretization numerical calculation method comprises the following steps:
Figure FDA0003458250770000034
Figure FDA0003458250770000035
CS(xT，yTt) represents the concentration of spilled oil on the water surface, CD(xT，yTT) represents the concentration of spilled oil on the surface of the water dissolved in the concentration of spilled oil in the water, nSRepresenting the total number of oil particles on the water surface, nDRepresents the total number of oil particles dissolved in water, ρ represents the oil spill density, δ xT，δyTRepresenting the length of grid sides in the longitude and latitude directions, x, when grid discretization is performed during grid calculation on a computerT，yTRepresenting longitude and latitude coordinates of the grid.
2. The marine oil spill trajectory simulation method based on thickness information of claim 1, wherein,
Figure FDA0003458250770000036
and
Figure FDA0003458250770000037
the calculation method of (c) is as follows:
Figure FDA0003458250770000038
Figure FDA0003458250770000039
Figure FDA00034582507700000310
indicating the rate of change in volume resulting from evaporation,
Figure FDA00034582507700000311
indicating the rate of change of the volume of spill oil dissolved in water,
Figure FDA0003458250770000041
indicating the rate of change of the volume of spill oil on the water surface.
3. The marine oil spill trajectory simulation method based on thickness information of claim 1, wherein the position increment of the oil particles on the water surface and dissolved in the water is calculated as follows:
wherein the position increment of the oil particles on the water surface is calculated as follows:
dxkS(t)＝[UC(xk，yk，0，t)+UW(xk，yk，t)+US(xk，yk，t)]dt+dx′k(t)，
dxkS(t) represents the increase in the position of the oil particles on the water surface, UCIs the euler velocity component;
UW＝(UW，VW) Is the effect of ocean currents, related to the strength of sea surface wind and the angle between wind and ocean currents, UW，VWIs the force of the ocean current along two perpendicular components;
UW＝α(Wxcosβ+Wysinβ)，
VW＝α(-Wxsinβ+Wycosβ)
Wxand WyRespectively, the intensity of the wind at the predetermined location in the latitude and longitude directions, alpha is the drift factor, beta is the drift angle,
US＝(US，VS) Is the wave action, the calculation method is shown as the following formula:
Figure FDA0003458250770000042
wherein ω is angular frequency, S is wave spectrum, k is wave number;
dx′k(t) is the random portion of the oil particle motion,
Figure FDA0003458250770000043
wherein r is a random number between 0 and 1, khIs a constant; dt represents the sign of the differential; dt represents a time increment;
the position increment of oil particles dissolved in water is calculated as follows:
dxkD(t)＝UC(xk，yk，zk，t)dt+dx′k(t)，
dxkD(t) represents the increase in the position of oil particles dissolved in water.
CN202010127439.6A 2020-02-28 2020-02-28 Offshore oil spill track simulation method based on thickness information Active CN111368414B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010127439.6A CN111368414B (en) 2020-02-28 2020-02-28 Offshore oil spill track simulation method based on thickness information

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010127439.6A CN111368414B (en) 2020-02-28 2020-02-28 Offshore oil spill track simulation method based on thickness information

Publications (2)

Publication Number Publication Date
CN111368414A CN111368414A (en) 2020-07-03
CN111368414B true CN111368414B (en) 2022-04-08

Family

ID=71208239

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010127439.6A Active CN111368414B (en) 2020-02-28 2020-02-28 Offshore oil spill track simulation method based on thickness information

Country Status (1)

Country Link
CN (1) CN111368414B (en)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102156817A (en) * 2011-04-11 2011-08-17 国家海洋信息中心 Rapid prediction method for drifting path and diffusion of marine spilled oil
CN107944608A (en) * 2017-11-16 2018-04-20 福建四创软件有限公司 Sea drift thing and oil drift and diffusion forecasting procedure based on satellite remote sensing
CN108268751A (en) * 2018-01-26 2018-07-10 国家海洋局第海洋研究所 Underwater oil spilling motion value analogy method and underwater oil spilling Behavior modeling system
CN109284865A (en) * 2018-09-05 2019-01-29 蓝海天智（舟山）海洋科技有限公司 A kind of marine oil overflow motion profile prediction technique based on the oil fence

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102156817A (en) * 2011-04-11 2011-08-17 国家海洋信息中心 Rapid prediction method for drifting path and diffusion of marine spilled oil
CN107944608A (en) * 2017-11-16 2018-04-20 福建四创软件有限公司 Sea drift thing and oil drift and diffusion forecasting procedure based on satellite remote sensing
CN108268751A (en) * 2018-01-26 2018-07-10 国家海洋局第海洋研究所 Underwater oil spilling motion value analogy method and underwater oil spilling Behavior modeling system
CN109284865A (en) * 2018-09-05 2019-01-29 蓝海天智（舟山）海洋科技有限公司 A kind of marine oil overflow motion profile prediction technique based on the oil fence

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
Oil spill impact modeling: development and validation;Deborah.P等;《Environmental Toxicology and Chemistry》;20041031;第23卷(第10期);2441-2456 *
三维溢油数值模式研究及其在近海的应用;郭为军;《中国博士学位论文全文数据库工程科技Ⅰ辑》;20120615;B027-92 *
渤海海洋溢油数学模型基本理论及应用研究;吴丹;《中国优秀硕士学位论文全文数据库工程科技Ⅰ辑》;20120215;B027-334 *

Also Published As

Publication number Publication date
CN111368414A (en) 2020-07-03

Similar Documents

Publication Publication Date Title
CN107944608B (en) Sea surface drift and oil spill drift diffusion forecasting method based on satellite remote sensing
Lauritzen et al. Atmospheric transport schemes: desirable properties and a semi-Lagrangian view on finite-volume discretizations
CN111650672B (en) Method for realizing rapid prediction of atmospheric diffusion of air pollutants by adopting time stack
CN106777917B (en) Hydraulic structure calculates maritime affairs traffic control radar shaded areas and influences evaluation method
US20130304430A1 (en) Three-Dimensional Tracer Dispersion Model
CN109783846B (en) GNSS ocean buoy-based sea level measurement uncertainty evaluation method
Mahpeykar et al. Numerical modelling the effect of wind on Water Level and Evaporation Rate in the Persian Gulf
Hu et al. Numerical model of tides in Pearl River estuary with moving boundary
CN111368414B (en) Offshore oil spill track simulation method based on thickness information
Zhang et al. Regional ocean current field construction based on an empirical bayesian kriging algorithm using multiple underwater gliders
Allwine et al. MELSAR: a mesoscale air quality model for complex terrain. Volume 1. Overview, technical description and user's guide
Cai et al. Flux footprints in the convective boundary layer: large-eddy simulation and lagrangian stochastic modelling
CN108182320B (en) Inland waterway-oriented free drift object stranding probability distribution calculation method
Morita et al. Tsunami-induced drift motion and numerical simulation using various models
Oporto et al. Analysis of contaminant transport under wind conditions on the surface of a shallow lake
Menshutkin et al. Ecosystem modeling of freshwater lakes (Review): 1. Hydrodynamics of lakes
Deleersnijder et al. Num erical Sim ulation of the Three-Dim ensional Tidal Circulation in an Island’s wake
CN117556719B (en) Hydrodynamic model verification method based on natural tracer ions
Shi et al. Spatial interpolation of tidal data using a multiple-order harmonic equation for unstructured grids
Mu et al. The optimization of maritime search and rescue simulation system based on CPS
Zounemat-Kermani et al. Conjunction of 2D and 3D modified flow solvers for simulating spatio-temporal wind induced hydrodynamics in the Caspian Sea
KR102602962B1 (en) System and method for simulation of marine pollution dispersion using numerical tracer technique
Wulandari et al. Two dimensional Hydrodynamic Simulation in Labuan Bajo Waters, Indonesia.
Gao et al. A Method Using the Huber Function for Inversion of Tidal Open Boundary Conditions of the M 2 Constituent in the Bohai and Yellow Seas
US20240363201A1 (en) System and method for simulation of marine pollution dispersion using numerical tracer technique

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant