Method and system for predicting oil film migration path after oil contamination leakage in reservoir construction and related equipment

Abstract

The invention discloses a method for predicting an oil film migration path after oil contamination leakage in reservoir construction, and the method comprises the steps: employing a three-dimensional reservoir hydrodynamic model to calculate a reservoir flow field according to the specific water flow condition of a reservoir, and enabling the reservoir flow field to serve as a flow field foundation for driving an oil film to move; on the basis, oil film evaporation calculation is introduced, and meanwhile various condition factors such as crude oil types, air temperature, oil spilling area, wind speed, solar radiation and oil film thickness are introduced, so that the accuracy of the total amount in the oil film migration process is effectively improved.

Description

Technical Field

[0001] The invention relates to the field of water conservancy projects, and in particular to a prediction technology for an oil film migration path after oil pollution leakage during reservoir construction. Background Art

[0002] With the rapid development of reservoir construction, it plays a vital role in water resource allocation, flood prevention and disaster reduction, and improvement of the ecological environment. However, during the construction and operation of reservoirs, the potential risk of oil leakage cannot be ignored. Oil leakage will not only cause direct pollution to the water quality of the reservoir, but may also affect the downstream ecological environment and human health through water flow. Therefore, accurately predicting the migration path of the oil film after oil leakage during reservoir construction is of great significance for taking emergency measures in time and controlling the scope of pollution spread.

[0003] At present, there is a lack of a specific prediction method for the migration path of oil film after oil leakage during reservoir construction. The existing prediction method for oil migration in water bodies is mainly based on the flow field calculated by the convection-diffusion equation. This oil migration prediction method is mostly applicable to open rivers, lakes or marine environments, and fails to fully consider the unique water flow conditions and environmental factors of reservoirs.

[0004] If the existing scheme based on the convection-diffusion equation to calculate the flow field to predict the migration path of oil pollution is used to calculate and predict the migration path of oil film after oil pollution leakage during reservoir construction, the following problems will exist:

[0005] (1) Incomplete consideration of factors and insufficient prediction accuracy. The traditional method of predicting oil migration based on the flow field calculated by the convection-diffusion equation often uses two-dimensional water flow simulation to calculate the flow field, and uses the vertical average flow to replace the surface flow velocity. Although it can simulate the diffusion trend of oil pollution to a certain extent, when applied to the reservoir environment, it often ignores the unique water flow conditions of the reservoir. The water flow velocity of the reservoir may show complex changes due to various factors such as topography, storage capacity, and water release. In the deep reservoir area, the water flow velocity may be relatively slow, which is conducive to the sedimentation of suspended matter in the water body; while near the water outlet or spillway, the water flow velocity will increase significantly, forming a turbulent water flow. Therefore, a three-dimensional refined flow field calculation method is required when calculating the reservoir flow field.

[0006] (2) Existing solutions often ignore the evaporation factor of the oil film when calculating the oil film migration path in the river, resulting in large errors in the calculation of the oil film area.

[0007] The above factors have a significant impact on the prediction results of the oil film migration path, which will greatly affect the accuracy and reliability of the prediction results of the oil film migration path after oil leakage during reservoir construction.

[0008] Due to the failure to fully consider the above-mentioned key factors, existing technologies often find it difficult to accurately simulate the diffusion range, speed and direction of the oil film when predicting the migration path of the oil film after an oil leak in a reservoir, resulting in a large deviation between the predicted results and the actual situation, and failing to meet the needs of reservoir environmental risk assessment and emergency management. Summary of the invention

[0009] In view of the problem that the prediction results of the prior art for predicting the migration path of the oil film after an oil leak in a reservoir have large errors in accuracy and reliability, the purpose of the present invention is to provide a prediction scheme for the migration path of the oil film after an oil leak during reservoir construction. The prediction scheme innovatively introduces factors such as reservoir-specific water flow conditions and evaporation of the oil film, effectively improving the accuracy and reliability of the prediction results.

[0010] In order to achieve the above object, the present invention provides a method for predicting the migration path of oil film after oil leakage in reservoir construction, the prediction method comprising:

[0011] (1) Construct a three-dimensional reservoir hydrodynamic model;

[0012] On the basis of the three-dimensional hydrodynamic model, the three-dimensional topographic data of the reservoir and the basic water conservancy data are introduced to form a three-dimensional reservoir hydrodynamic model, and the real-time monitoring data of the reservoir is integrated into the three-dimensional reservoir hydrodynamic model to perform real-time correction on the three-dimensional reservoir hydrodynamic model;

[0013] (2) calculating the flow field of the reservoir area based on the three-dimensional reservoir hydrodynamic model constructed in step (1);

[0014] (3) coupled oil spill module;

[0015] Based on the historical data of oil spill accidents and the corresponding physical and chemical properties of oil spilled substances, an oil spill module is constructed, and the constructed oil spill module is coupled with the three-dimensional reservoir hydrodynamic model. By exchanging the reservoir flow field data calculated by the three-dimensional reservoir hydrodynamic model in real time with the corresponding oil spill data information in the oil spill module, the oil spill module and the three-dimensional reservoir hydrodynamic model work together to jointly simulate the movement trajectory and fate of the oil spill in the water.

[0016] (4) Integrate the oil film evaporation calculation module;

[0017] Based on the oil spill accident data, the physical and chemical properties of the corresponding oil film and the corresponding meteorological data, an oil film evaporation calculation module is constructed, and the constructed oil film evaporation calculation module is integrated with the oil spill module and the three-dimensional reservoir hydrodynamic model. By exchanging the corresponding oil spill data information in the oil spill module with the corresponding calculated oil film evaporation data information in the oil film evaporation calculation module in real time, the collaborative work of the two modules is realized, thereby simulating the evaporation process of the oil film.

[0018] In some embodiments of the present invention, the prediction method is configured to use a particle tracking method to simulate the convection, expansion and turbulent diffusion process of the oil spill in the oil spill module constructed in step (3), and to use a random walk method to simulate the turbulent diffusion process of the oil spill in the horizontal direction.

[0019] In some embodiments of the present invention, the oil film evaporation calculation module constructed in step (4) of the prediction method is configured to simulate the evaporation process of the oil film based on the mass conservation and energy conservation modes.

[0020] In order to achieve the above object, the present invention provides a prediction system for the migration path of oil film after oil leakage during reservoir construction, the prediction system comprising:

[0021] A three-dimensional reservoir hydrodynamic model, wherein the three-dimensional reservoir hydrodynamic model is formed by introducing three-dimensional reservoir terrain data and water conservancy basic data on the basis of the three-dimensional hydrodynamic model, and the three-dimensional reservoir hydrodynamic model is configured to calculate the flow field of the reservoir area;

[0022] An oil spill module, which is formed based on historical data of oil spill accidents and corresponding physical and chemical property data of oil spill substances, and is configured to be coupled with a three-dimensional reservoir hydrodynamic model, so that the oil spill module and the three-dimensional reservoir hydrodynamic model can work together to jointly simulate the movement trajectory and fate of the oil spill in the water by exchanging reservoir area flow field data calculated by the three-dimensional reservoir hydrodynamic model and corresponding oil spill data information in the oil spill module in real time;

[0023] The oil film evaporation calculation module is composed of oil spill accident data, physical and chemical property data of the corresponding oil film and corresponding meteorological data, and is configured to be integrated with the oil spill module and the three-dimensional reservoir hydrodynamic model. The corresponding oil spill data information in the oil spill module and the corresponding oil film evaporation data information calculated and determined in the oil film evaporation calculation module are exchanged in real time to achieve the collaborative work of the two modules, thereby simulating the evaporation process of the oil film.

[0024] In some embodiments of the present invention, the oil spill module is configured to use a particle tracking method to simulate the convection, expansion and turbulent diffusion process of the oil spill, and a random walk method to simulate the turbulent diffusion process of the oil spill in the horizontal direction.

[0025] In some embodiments of the present invention, the oil film evaporation calculation module is configured to simulate the evaporation process of the oil film based on mass conservation and energy conservation modes.

[0026] In order to achieve the above object, the present invention also provides a computer-readable storage medium having a program stored thereon, which, when executed by a processor, implements the steps of the method for predicting the migration path of the oil film after oil leakage during reservoir construction.

[0027] In order to achieve the above object, the present invention also provides a processor, which is used to run a program, and when the program is running, the steps of the method for predicting the migration path of oil film after oil leakage in reservoir construction are executed.

[0028] In order to achieve the above-mentioned purpose, the present invention also provides a terminal device, which includes a processor, a memory, and a program stored in the memory and executable on the processor, wherein the program code is loaded and executed by the processor to implement the steps of the method for predicting the migration path of the oil film after oil leakage during the construction of the above-mentioned reservoir.

[0029] In order to achieve the above object, the present invention also provides a computer program product, which, when executed on a data processing device, is suitable for executing the steps of the method for predicting the migration path of oil film after oil leakage in reservoir construction.

[0030] The prediction scheme of the oil film migration path after oil leakage in reservoir construction provided by the present invention has the following technical features compared with the prior art:

[0031] (1) Highly targeted: The solution of the present invention is specifically aimed at predicting the migration path of oil film after oil leakage during reservoir construction. It fully considers and introduces the unique water flow conditions, topography and environmental factors of the reservoir, providing more accurate predictions for reservoir environmental risk assessment and emergency management. On this basis, the three-dimensional refined flow field calculation method is adopted to more accurately simulate the water flow velocity and direction at different locations and depths in the reservoir, thereby more realistically reflecting the migration process of the oil film in the reservoir and improving the accuracy of the prediction.

[0032] 2. Comprehensive consideration: In the prediction process, the scheme of the present invention not only introduces the water flow conditions of the reservoir, but also fully considers the evaporation factor of the oil film and introduces the oil film evaporation calculation, which helps to more accurately calculate the mass change and area change of the oil film in the diffusion process and reduce the prediction error.

[0033] 3. High prediction accuracy: The scheme of the present invention comprehensively considers multiple key factors and adopts innovative calculation methods, so that the prediction results of the present invention have high accuracy and reliability, and can better meet the needs of reservoir environmental risk assessment and emergency management. BRIEF DESCRIPTION OF THE DRAWINGS

[0034] The present invention is further described below in conjunction with the accompanying drawings and specific embodiments.

[0035] Figure 1 It is a reservoir basin grid and an enlarged view of the reservoir area in the example of the present invention;

[0036] Figure 2 This is an example diagram of a model after interpolating terrain in an example of the present invention;

[0037] Figure 3 An example diagram of the surface hydrodynamic flow field calculated in an example of the present invention;

[0038] Figure 4 This is an example diagram of oil film distribution 3 hours after the reservoir oil spill accident in the example of the present invention;

[0039] Figure 5 It is an example diagram of the oil film thickness process line at the water intake of the water plant in the example of the present invention;

[0040] Figure 6 TABLE 1 is a table showing the maximum thickness (μm) of the oil film over time in the examples of the present invention. DETAILED DESCRIPTION

[0041] In order to make the technical means, creative features, objectives and effects achieved by the present invention easy to understand, the present invention is further described below with reference to specific diagrams.

[0042] In view of the unique water flow conditions of reservoirs, the scheme of the present invention provides a method for predicting the migration path of the oil film after oil leakage during reservoir construction. This prediction method can be effectively used to predict the migration path of the oil film after oil leakage during reservoir construction by comprehensively introducing factors such as reservoir water flow characteristics and evaporation of the oil film, and ensure the accuracy and reliability of the prediction results.

[0043] The scheme of the present invention provides a method for predicting the migration path of an oil film after an oil spill during reservoir construction. Aiming at the unique water flow conditions of the reservoir, a three-dimensional reservoir hydrodynamic model is used to calculate the reservoir flow field, which is used as the flow field basis for driving the movement of the oil film. On this basis, an oil film evaporation calculation is introduced, and at the same time, multiple conditional factors such as crude oil type, temperature, oil spill area, wind speed, solar radiation and oil film thickness are introduced, thereby effectively improving the accuracy of the total amount during the oil film migration process.

[0044] Based on this, the method for predicting the migration path of oil film after oil leakage in reservoir construction provided by the present invention specifically includes the following steps:

[0045] (1) Construct a three-dimensional reservoir hydrodynamic model;

[0046] On the basis of the three-dimensional hydrodynamic model, the three-dimensional topographic data of the reservoir and the basic water conservancy data are introduced to form a three-dimensional reservoir hydrodynamic model, and the real-time monitoring data of the reservoir is integrated into the three-dimensional reservoir hydrodynamic model to perform real-time correction on the three-dimensional reservoir hydrodynamic model;

[0047] (2) calculating the flow field of the reservoir area based on the three-dimensional reservoir hydrodynamic model constructed in step (1);

[0048] (3) coupled oil spill module;

[0049] Based on the historical data of oil spill accidents and the corresponding physical and chemical properties of oil spilled substances, an oil spill module is constructed, and the constructed oil spill module is coupled with the three-dimensional reservoir hydrodynamic model. By exchanging the reservoir flow field data calculated by the three-dimensional reservoir hydrodynamic model in real time with the corresponding oil spill data information in the oil spill module, the oil spill module and the three-dimensional reservoir hydrodynamic model work together to jointly simulate the movement trajectory and fate of the oil spill in the water.

[0050] (4) Integrate the oil film evaporation calculation module;

[0051] Based on the oil spill accident data, the physical and chemical properties of the corresponding oil film and the corresponding meteorological data, an oil film evaporation calculation module is constructed, and the constructed oil film evaporation calculation module is integrated with the oil spill module and the three-dimensional reservoir hydrodynamic model. By exchanging the corresponding oil spill data information in the oil spill module with the corresponding calculated oil film evaporation data information in the oil film evaporation calculation module in real time, the collaborative work of the two modules is realized, thereby simulating the evaporation process of the oil film.

[0052] As a further explanation, in step (1), the prediction method first constructs an accurate three-dimensional reservoir hydrodynamic model based on high-precision topographic data, hydrological and meteorological data, and real-time monitoring data.

[0053] Then, based on the three-dimensional reservoir hydrodynamic model, the water flow movement of the reservoir is simulated, including flow velocity, flow direction, water level changes, etc., to provide an accurate flow field background for subsequent oil spill simulation.

[0054] When constructing a three-dimensional reservoir hydrodynamic model, the basic data of the reservoir is first obtained, specifically the three-dimensional terrain data of the reservoir.

[0055] Specifically, high-precision reservoir geospatial data can be obtained through UAV oblique photography and satellite remote sensing images, including oblique photography models, DEM (digital elevation model) and DOM (digital orthophoto), etc. The three-dimensional terrain data of the reservoir obtained in this way will be able to provide accurate topographic information of the reservoir and can serve as the basis for building a three-dimensional reservoir model.

[0056] Next, further obtain the basic water conservancy data corresponding to the reservoir.

[0057] Specifically, here we obtain the basic water conservancy data corresponding to the reservoir by integrating a local water conservancy map, including the reservoir's geometric shape, storage capacity curve, water level-reservoir capacity relationship, etc.; it also includes the reservoir's inflow, outflow, leakage and other hydrological data, as well as meteorological, geological and other information around the reservoir.

[0058] On this basis, a three-dimensional hydrodynamic model is used, and the collected three-dimensional topographic data of the reservoir and the basic water conservancy data are input into the three-dimensional hydrodynamic model. The three-dimensional hydrodynamic model is calculated to form a three-dimensional reservoir hydrodynamic model, which can simulate the water flow movement in the reservoir, including flow velocity, flow direction, water level changes, etc. The three-dimensional hydrodynamic model used here can be a DELFT3D model or a similar three-dimensional hydrodynamic model.

[0059] Finally, the reservoir data is monitored in real time to obtain real-time monitoring data of the reservoir, such as water level, flow, rainfall, etc.; and the real-time monitoring data of the reservoir is input into the three-dimensional reservoir hydrodynamic model to integrate the real-time monitoring data of the reservoir into the three-dimensional reservoir hydrodynamic model, thereby performing real-time correction and verification of the three-dimensional reservoir hydrodynamic model.

[0060] The three-dimensional reservoir hydrodynamic model thus formed can realize accurate simulation of reservoir water flow movement by combining the real-scene three-dimensional model with the hydrodynamic model.

[0061] As a further explanation, in step (2), when calculating the reservoir flow field based on the three-dimensional reservoir hydrodynamic model, the prediction method calculates the reservoir flow field by solving the control equation of water flow movement, thereby obtaining the water flow movement data in the reservoir, including flow velocity, flow direction, water level changes, etc.

[0062] As an example, the governing equation of the water flow motion here may be the Navier-Stokes equation, etc. In the calculation process, an appropriate numerical method (such as finite element method, finite volume method, etc.) is used to discretize the governing equation, and a numerical solution of the flow field is obtained through iterative solution.

[0063] As a further explanation, step (3) of the prediction method constructs a corresponding oil spill module based on the reservoir flow field calculated in step (2), thereby accurately simulating the movement trajectory and destination of the oil spill in the water.

[0064] In this step, when constructing the oil spill module, the historical data of oil spill accidents and the physical and chemical properties data of the corresponding oil spill substances in the oil spill accidents are first obtained and integrated.

[0065] The historical data of oil spill accidents here include the type, quantity, location, time, etc. of oil spill;

[0066] The physical and chemical properties of oil spill materials include density, viscosity, volatility, solubility, etc.

[0067] Next, a corresponding oil spill module is constructed based on these data. The constructed oil spill module can effectively simulate the diffusion, volatilization, emulsification and dissolution processes of oil spills.

[0068] Specifically, when constructing the oil spill module, the solution of the present invention is configured to adopt a particle tracking method in combination with a random walk method.

[0069] Among them, the particle tracking method is used to simulate the movement trajectory of the oil spill in the water, including convection, expansion and turbulent diffusion. The random walk method is used to simulate the turbulent diffusion of the oil spill in the horizontal direction, making the simulation results closer to the actual situation;

[0070] On this basis, the Lagrange equation is further introduced to solve the vertical diffusion, that is, the Lagrange equation is used to solve the diffusion process of the oil spill in the vertical direction at the same time.

[0071] The oil spill module thus formed comprehensively introduces the influence of factors such as water flow velocity, diffusion coefficient and gravity on the vertical diffusion of oil spill, making the simulation results more accurate.

[0072] The solution of the present invention is based on the constructed oil spill module, which is further coupled with the three-dimensional reservoir hydrodynamic model constructed in step (1), thereby accurately simulating the movement trajectory and destination of the oil spill in the water.

[0073] The solution of the present invention couples the oil spill module in the three-dimensional reservoir hydrodynamic model and configures the reservoir water flow movement data information (i.e., reservoir flow field) calculated by the three-dimensional reservoir hydrodynamic model to exchange the oil spill data information in the oil spill module in real time, so that the two modules work together to jointly simulate the movement trajectory and destination of the oil spill in the water. This coupling method not only improves the accuracy of the model, but also makes the simulation results more valuable for practical applications.

[0074] The scheme of the present invention couples the oil spill module in the three-dimensional reservoir hydrodynamic model, and introduces a particle tracking method to simulate the convection, expansion and turbulent diffusion process of the oil spill on the basis of the three-dimensional reservoir hydrodynamic model. The random walk method is simultaneously used to simulate the turbulent diffusion in the horizontal direction and the Langeven equation is solved in the vertical direction, thereby accurately simulating the movement trajectory and destination of the oil spill in the water.

[0075] As a further explanation, in step (4), this prediction method further introduces oil film evaporation calculation on the basis of step (3), thereby further improving the accuracy of the simulation prediction results of the entire oil film migration path.

[0076] This step first collects and integrates relevant data on oil spill accidents, physical and chemical properties of the corresponding oil film in the oil spill accident, and corresponding meteorological data, thereby constructing a corresponding oil film evaporation calculation module.

[0077] The relevant data of the integrated oil spill accidents here include basic information such as oil spill type, quantity, location, time, etc.

[0078] Collecting data on the physical and chemical properties of the oil film, such as oil density, viscosity, volatility, solubility, etc., is crucial for simulating the evaporation process of the oil film.

[0079] Obtain meteorological data, including wind speed, wind direction, temperature, humidity, etc., which will affect the evaporation rate of the oil film.

[0080] The oil film evaporation calculation module constructed here is configured as an evaporation model based on mass conservation and energy conservation modes, thereby simulating the evaporation process of the oil film.

[0081] Then, the constructed oil film evaporation calculation module was introduced into the three-dimensional reservoir hydrodynamic model that had been coupled with the oil spill module.

[0082] On this basis, the corresponding oil spill data information in the oil spill module is exchanged in real time with the corresponding oil film evaporation data information calculated and determined in the oil film evaporation calculation module to achieve the collaborative work of the two modules, thereby simulating the evaporation process of the oil film.

[0083] In the scheme of the present invention, the scheme based on steps (1)-(3) realizes the simulation prediction of the oil film migration path after oil leakage during reservoir construction. On this basis, the oil film evaporation calculation is further introduced through step (4) to realize the introduction of oil film evaporation conditions in the simulation prediction process of the oil film migration path, thereby further improving the accuracy of the simulation prediction of the oil film migration path after oil leakage during reservoir construction.

[0084] Here, the implementation scheme of this scheme, which further incorporates the oil film evaporation factor on the basis of constructing the oil spill module to improve the simulation accuracy, is further explained.

[0085] When constructing the oil spill module, this solution adds the collection of parameters related to oil film evaporation, such as oil temperature, ambient humidity, wind speed, solar radiation, etc., on the basis of data collection. These parameters have a direct impact on the evaporation rate of the oil film.

[0086] Preprocess the collected data to ensure its accuracy and completeness, and provide a reliable basis for subsequent simulations.

[0087] On this basis, based on the convection-diffusion equation and incorporating the physical mechanism of oil film evaporation, a comprehensive model is formed that can simultaneously simulate oil diffusion, drift and evaporation. Furthermore, according to parameters such as oil type, environmental conditions and leakage, the model is carefully optimized, especially for parameters related to oil film evaporation rate, to ensure the accuracy of the simulation results.

[0088] Specifically, the scheme of the present invention constructs a corresponding oil film evaporation submodule (i.e., the aforementioned oil film evaporation calculation module) in the oil spill module to simulate the evaporation process of the oil film in the reservoir. At the same time, through the organic coordination between the oil film evaporation submodule and the oil spill module, the diffusion, drift and evaporation of oil pollution can be simulated simultaneously.

[0089] The oil film evaporation submodule here introduces a variety of factors that affect oil film evaporation, such as oil temperature, ambient humidity, wind speed, solar radiation, etc., and calculates the evaporation rate of the oil film based on these factors.

[0090] Based on the organic coordination between the oil film evaporation submodule and the oil spill module, the oil film evaporation rate can be combined with the oil diffusion, drift and other processes to achieve a more accurate prediction of the oil migration path.

[0091] Furthermore, the solution of the present invention organically couples the oil spill module of the integrated oil film evaporation submodule with the three-dimensional reservoir hydrodynamic model to ensure real-time exchange and sharing of data between the modules.

[0092] Finally, the scheme of the present invention further conducts a comprehensive verification of the integrated coupled model, and evaluates the accuracy and reliability of the model by comparing it with actual oil spill incidents or known data.

[0093] The method for predicting the migration path of the oil film after oil leakage in the construction of the reservoir provided by the solution of the present invention can constitute a corresponding software program in specific application, forming a corresponding prediction system for the migration path of the oil film after oil leakage in the construction of the reservoir. When the software program is running, it will execute the above-mentioned smooth curve construction method and store it in the corresponding storage medium for the processor to call and execute.

[0094] The resulting prediction system for the migration path of oil film after oil leakage during reservoir construction mainly includes three functional modules: a three-dimensional reservoir hydrodynamic model, an oil spill module, and an oil film evaporation calculation module.

[0095] Specifically, the three-dimensional reservoir hydrodynamic model in the prediction system is formed by introducing three-dimensional reservoir terrain data and water conservancy basic data on the basis of the three-dimensional hydrodynamic model, and is configured to be able to calculate the flow field of the reservoir area.

[0096] The three-dimensional reservoir hydrodynamic model thus formed can be used to calculate the flow field in the reservoir area, which can more accurately simulate the water flow velocity and direction at different locations and depths in the reservoir, thereby improving the accuracy of predicting the oil film migration path.

[0097] The oil spill module in this prediction system is formed based on the historical data of oil spill accidents and the corresponding physical and chemical properties data of the oil spilled substances, and is configured to be coupled with the three-dimensional reservoir hydrodynamic model. By real-time exchange of the reservoir flow field data calculated by the three-dimensional reservoir hydrodynamic model and the corresponding oil spill data information in the oil spill module, the oil spill module and the three-dimensional reservoir hydrodynamic model work together to jointly simulate the movement trajectory and destination of the oil spill in the water.

[0098] The oil film evaporation calculation module in this prediction system is based on oil spill accident data, the physical and chemical properties of the corresponding oil film, and the corresponding meteorological data, and is configured to be integrated with the oil spill module and the three-dimensional reservoir hydrodynamic model. By exchanging the corresponding oil spill data information in the oil spill module with the corresponding calculated oil film evaporation data information in the oil film evaporation calculation module in real time, the two modules work together to simulate the evaporation process of the oil film.

[0099] Therefore, the oil overflow module and the oil film evaporation calculation module are introduced into the prediction system, so that the evaporation factor of the oil film is effectively introduced into the prediction system for prediction calculation, which can more comprehensively consider the mass change of the oil film during the diffusion process, thereby reducing the prediction error.

[0100] As a further explanation, when constructing the three-dimensional reservoir hydrodynamic model in this prediction system, the basic data of the reservoir, specifically the three-dimensional terrain data of the reservoir, is first obtained.

[0101] Specifically, high-precision reservoir geospatial data can be obtained through UAV oblique photography and satellite remote sensing images, including oblique photography models, DEM (digital elevation model) and DOM (digital orthophoto), etc. The three-dimensional terrain data of the reservoir obtained in this way will be able to provide accurate topographic information of the reservoir and can serve as the basis for building a three-dimensional reservoir model.

[0102] Next, further obtain the basic water conservancy data corresponding to the reservoir.

[0103] Specifically, here we obtain the basic water conservancy data corresponding to the reservoir by integrating a local water conservancy map, including the reservoir's geometric shape, storage capacity curve, water level-reservoir capacity relationship, etc.; it also includes the reservoir's inflow, outflow, leakage and other hydrological data, as well as meteorological, geological and other information around the reservoir.

[0104] On this basis, a three-dimensional hydrodynamic model is used, and the collected three-dimensional topographic data of the reservoir and the basic water conservancy data are input into the three-dimensional hydrodynamic model. The three-dimensional hydrodynamic model is calculated to form a three-dimensional reservoir hydrodynamic model, which can simulate the water flow movement in the reservoir, including flow velocity, flow direction, water level changes, etc. The three-dimensional hydrodynamic model used here can be a DELFT3D model or a similar three-dimensional hydrodynamic model.

[0105] Finally, the reservoir data is monitored in real time to obtain real-time monitoring data of the reservoir, such as water level, flow, rainfall, etc.; and the real-time monitoring data of the reservoir is input into the three-dimensional reservoir hydrodynamic model to integrate the real-time monitoring data of the reservoir into the three-dimensional reservoir hydrodynamic model, thereby performing real-time correction and verification of the three-dimensional reservoir hydrodynamic model.

[0106] The three-dimensional reservoir hydrodynamic model thus formed can realize accurate simulation of reservoir water flow movement by combining the real-scene three-dimensional model with the hydrodynamic model.

[0107] Furthermore, when calculating the flow field of the reservoir area, the three-dimensional reservoir hydrodynamic model calculates the flow field of the reservoir area by solving the control equation of water flow movement, thereby obtaining the water flow movement data in the reservoir, including flow velocity, flow direction, water level changes, etc.

[0108] As a further explanation, when constructing the oil spill module in this prediction system, the historical data of oil spill accidents and the physical and chemical properties data of the corresponding oil spill substances in the oil spill accidents are first obtained and integrated.

[0109] The historical data of oil spill accidents here include the type, quantity, location, time, etc. of oil spill;

[0110] The physical and chemical properties of oil spill materials include density, viscosity, volatility, solubility, etc.

[0111] Next, a corresponding oil spill module is constructed based on these data. The constructed oil spill module can effectively simulate the diffusion, volatilization, emulsification and dissolution processes of oil spills.

[0112] Specifically, when constructing the oil spill module, it is configured to adopt a particle tracking method in combination with a random walk method.

[0113] Among them, the particle tracking method is used to simulate the movement trajectory of the oil spill in the water, including convection, expansion and turbulent diffusion. The random walk method is used to simulate the turbulent diffusion of the oil spill in the horizontal direction, making the simulation results closer to the actual situation;

[0114] On this basis, the Lagrange equation is further introduced to solve the vertical diffusion, that is, the Lagrange equation is used to solve the diffusion process of the oil spill in the vertical direction at the same time.

[0115] The oil spill module thus formed comprehensively introduces the influence of factors such as water flow velocity, diffusion coefficient and gravity on the vertical diffusion of oil spill, making the simulation results more accurate.

[0116] The solution of the present invention is based on the constructed oil spill module, which is further coupled with a three-dimensional reservoir hydrodynamic model, thereby accurately simulating the movement trajectory and destination of the oil spill in the water.

[0117] In this prediction system, the oil spill module is coupled to the three-dimensional reservoir hydrodynamic model, and configured to exchange the reservoir water flow movement data information (i.e., reservoir flow field) calculated by the three-dimensional reservoir hydrodynamic model with the oil spill data information in the oil spill module in real time, so that the two modules can work together to simulate the movement trajectory and destination of the oil spill in the water. This coupling method not only improves the accuracy of the model, but also makes the simulation results more valuable for practical applications.

[0118] In this prediction system, the oil spill module is coupled to the three-dimensional reservoir hydrodynamic model. On the basis of the three-dimensional reservoir hydrodynamic model, the particle tracking method is introduced to simulate the convection, expansion and turbulent diffusion process of the oil spill. The random walk method is used to simulate the turbulent diffusion in the horizontal direction and the Langeven equation is solved in the vertical direction, thereby accurately simulating the movement trajectory and destination of the oil spill in the water.

[0119] As a further explanation, when constructing the oil film evaporation calculation module in this prediction system, first collect and integrate the relevant data of the oil spill accident, the physical and chemical properties data of the corresponding oil film in the oil spill accident, and the corresponding meteorological data, so as to construct the corresponding oil film evaporation calculation module.

[0120] The relevant data of the integrated oil spill accidents here include basic information such as oil spill type, quantity, location, time, etc.

[0121] Collecting data on the physical and chemical properties of the oil film, such as oil density, viscosity, volatility, solubility, etc., is crucial for simulating the evaporation process of the oil film.

[0122] Obtain meteorological data, including wind speed, wind direction, temperature, humidity, etc., which will affect the evaporation rate of the oil film.

[0123] The oil film evaporation calculation module constructed here is configured as an evaporation model based on mass conservation and energy conservation modes, thereby simulating the evaporation process of the oil film.

[0124] Then, the constructed oil film evaporation calculation module was introduced into the three-dimensional reservoir hydrodynamic model that had been coupled with the oil spill module.

[0125] On this basis, the corresponding oil spill data information in the oil spill module is exchanged in real time with the corresponding oil film evaporation data information calculated and determined in the oil film evaporation calculation module to achieve the collaborative work of the two modules, thereby simulating the evaporation process of the oil film.

[0126] Compared with the prior art in predicting the migration path of oil film after oil leakage during reservoir construction, the prediction scheme of the migration path of oil film after oil leakage during reservoir construction provided by the present invention has the following characteristics and innovations:

[0127] (1) For specific environments:

[0128] Existing technologies: Mostly applicable to open rivers, lakes or ocean environments, and fail to fully consider the unique flow conditions and environmental factors of reservoirs. Mostly use two-dimensional water flow simulation calculations, using vertical average flow instead of surface flow velocity, which cannot accurately reflect the flow changes in the three-dimensional space of the reservoir.

[0129] The solution of the present invention fully considers the unique water flow conditions and environmental factors of the reservoir, and adopts a three-dimensional refined flow field calculation method, which can more accurately simulate the water flow velocity and direction at different positions and depths in the reservoir, thereby improving the accuracy of oil film migration path prediction.

[0130] (2) Introducing oil film evaporation factors:

[0131] Prior art: In the calculation of the migration path of river oil film, the evaporation factor of the oil film is often ignored, resulting in large errors in the calculation of the oil film area.

[0132] The solution of the present invention is to incorporate the evaporation factor of the oil film into the prediction model, and more comprehensively consider the mass change of the oil film during the diffusion process, thereby reducing the prediction error.

[0133] (3) Improved prediction accuracy and reliability:

[0134] Existing technology: Due to the failure to fully consider the above key factors, the prediction results often deviate greatly from the actual situation and cannot meet the needs of reservoir environmental risk assessment and emergency management.

[0135] The solution of the present invention: by comprehensively considering the reservoir-specific water flow conditions, environmental factors and the evaporation effect of the oil film, a more accurate simulation of the oil film migration path is achieved, the accuracy and reliability of the prediction results are improved, and the actual needs of reservoir environmental risk assessment and emergency management are better met.

[0136] It can be seen that compared with the prior art, the scheme of the present invention further constructs three-dimensional refined calculations by specifically introducing the unique water flow conditions, topography and environmental factors of the reservoir, and introduces the oil film evaporation factor, which can significantly improve the prediction accuracy and reliability, thereby more effectively dealing with the problem of oil leakage during the reservoir construction process.

[0137] In order to further illustrate the solution of the present invention, the implementation process of the solution of the present invention and the corresponding technical features are further illustrated below through specific application examples.

[0138] In specific applications, a prediction system for the migration path of oil film after oil leakage during reservoir construction is constructed based on the above-mentioned scheme, and is preferably presented in the form of a prediction model.

[0139] The prediction model here mainly includes three parts: (1) calculation of reservoir flow field based on three-dimensional reservoir hydrodynamic model; (2) coupling oil spill module; (3) calculation of oil film evaporation.

[0140] Specifically, (1) the flow field in the reservoir area is calculated based on the three-dimensional reservoir hydrodynamic model.

[0141] The three-dimensional reservoir hydrodynamic model here is specifically implemented by a three-dimensional hydrodynamic equation group. Under the Boussinesq and static assumptions, the incompressible viscous fluid is transformed from the Cartesian coordinate system to the horizontal non-orthogonal curve and vertical coordinate system, and its hydrodynamic control equation group is:

[0142]

[0143]

[0144] Among them, ξ=ξ(x,y), η=η(x,y), x, y and z are the east, north and vertical directions in the Cartesian coordinate system. ξ and η are the two axes of the horizontal non-orthogonal curvilinear coordinate system; σ is the axis of the vertical σ coordinate, which varies from -1 at the bottom of the river (z = -h) to 0 at the surface of the river (z = ζ).

[0145] In the above equations, (1) and (2) are momentum equations in the ξ and η directions respectively, and (3) is the continuity equation. Where ζ is the height of the river surface fluctuation, and D is the total water depth (D = h + ζ).

[0146] u1, v1 are the velocity components in the ξ and η directions respectively, and their expressions are:

[0147]

[0148] where u and v are the velocity components in the x and y directions respectively, and J is the Jacobian function (J = x ξ y η -x η y ξ , ), the subscript symbols ξ and η represent differential operations.

[0149] ω is the vertical flow velocity under the σ coordinate, and its expression is:

[0150]

[0151] In addition, some variables in the formula have the following meanings:

[0152] h3=y ξ y η +x ξ x η .

[0153] g is the acceleration due to gravity, f is the Coriolis parameter, K h is the vertical eddy diffusion coefficient, K m is the vertical eddy viscosity coefficient, ρ and ρ o are the perturbation and reference densities respectively. ξ and F η represent the horizontal diffusion terms of momentum in the ξ and η directions respectively, and can be expressed by equations (7) to (8) respectively:

[0154]

[0155] Among them A m is the horizontal viscosity coefficient, A h A is the horizontal diffusion coefficient, and generally the two values ​​are equal. m It can be calculated by the formula of Smagorinsky parameterization method:

[0156]

[0157] Where A is a constant, which is 0.02.

[0158] K m is the vertical viscosity coefficient, K h is the vertical diffusion coefficient, and both can be calculated using the 2.5-order turbulence closure model improved by Mellor and Yamada. The turbulent kinetic energy transport equation can be expressed as:

[0159]

[0160] The right side of the above equation is the turbulent kinetic energy shear generation term, buoyancy generation term, dissipation term, horizontal diffusion term and vertical diffusion term. 2 =u′ 2 +v′ 2 +ω′ 2 is one times of the turbulent kinetic energy, and l is the turbulent mixing length.

[0161] The turbulent mixing length transport equation is as follows:

[0162]

[0163] is the wall proximity function, which is defined as:

[0164]

[0165] Where κ is the Karman constant, which is taken as 0.4; d s and d bare the distances to the river surface and to the river bottom respectively; E1, E2 and B1 are all empirical constants, which are 1.8, 1.33 and 16.6 respectively.

[0166] Further, K m =qlS M , K h =qlS H , K q =qlS q .

[0167] Where S H , S M and S q are stable functions, and they satisfy the following relationship:

[0168]

[0169] S q =0.20 (15)

[0170] in Empirical constants: A1=0.92, A2=0.74, B2=10.1, C1=0.08.

[0171] q 2 and q 2 The boundary conditions on the river surface and river bottom are:

[0172]

[0173] Among them, u τs and u τb are the friction velocities on the river surface and river bottom, respectively.

[0174] This example solution uses the above three-dimensional hydrodynamic equations to build a three-dimensional reservoir hydrodynamic model. Compared with the existing technology, it has the following technical features:

[0175] (1) Application of three-dimensional non-orthogonal curvilinear coordinate system:

[0176] This example scheme uses a three-dimensional non-orthogonal curvilinear coordinate system to describe the water flow movement in the reservoir, which can more accurately reflect the complex topographic features of the reservoir and the changes in water flow in different directions.

[0177] By converting the Cartesian coordinate system into horizontal non-orthogonal curves and vertical coordinate systems, the water flow movement problem in the reservoir can be handled more flexibly, and the precision and accuracy of the simulation can be improved.

[0178] (2) Completeness of the hydrodynamic control equations:

[0179] This example solution provides a complete set of hydrodynamic control equations, including momentum equations and continuity equations, which can fully describe the laws of water flow in the reservoir.

[0180] The equations take into account a variety of physical factors, such as gravity, Coriolis force, horizontal diffusion and vertical diffusion, which have an important impact on water movement.

[0181] (3) Introduction of turbulent closure model:

[0182] This example scheme uses the improved 2.5-order turbulence closure model of Mellor and Yamada to calculate the vertical viscosity coefficient and vertical diffusion coefficient. This model can more accurately describe the turbulent motion in the reservoir.

[0183] By introducing the turbulent kinetic energy transport equation and the turbulent mixing length transport equation, we can have a deeper understanding of the impact of turbulent motion on water flow motion and further improve the simulation accuracy.

[0184] (3) Detailed processing of boundary conditions:

[0185] This example solution has made detailed treatment of the boundary conditions of the river surface and river bottom, and considered the influence of parameters such as friction velocity on the movement of water flow.

[0186] Through such fine processing, the actual situation of water flow movement in the reservoir can be reflected more accurately, and the reliability of the simulation results can be improved.

[0187] Therefore, this example solution will have the following technical effects:

[0188] (1) Can improve simulation accuracy:

[0189] By adopting a three-dimensional non-orthogonal curvilinear coordinate system, a complete set of hydrodynamic control equations, and sophisticated boundary condition processing, this scheme can more accurately simulate the water flow movement in the reservoir.

[0190] This is of great significance for predicting the migration path of the oil film after an oil spill, and can improve the accuracy and reliability of the prediction results.

[0191] (2) Ability to enhance adaptability:

[0192] This example solution can be applied to different types of reservoirs and different water flow conditions, and has strong adaptability and flexibility.

[0193] This makes this example solution have a wider application prospect and higher practical value in actual applications.

[0194] (3) Promoting scientific research:

[0195] This scheme provides new ideas and methods for the research in the field of reservoir hydrodynamics, which will help promote scientific research and technological progress in this field.

[0196] In summary, constructing a three-dimensional reservoir hydrodynamic model based on the above calculation formula not only improves the precision and accuracy of the simulation, but also enhances the adaptability and flexibility of the scheme, providing strong support for the emergency response and disposal of oil spills in the reservoir. It has significant technical characteristics and innovations, and embodies a high degree of technical effect.

[0197] (2) coupled oil spill module;

[0198] On the basis of calculating the flow field, the oil is divided into a finite number of particles, and the position of each particle is calculated according to the following formula:

[0199] X=X0+(U+αW 10 cosA+rcosB)Δt (18)

[0200] Y=Y0+(V+αW 10 sinA+rsinB)Δt (19)

[0201] Where: X0, Y0 are the coordinates of a particle at the previous moment; X1, Y1 are the coordinates of a particle at the next moment; U, V are the velocity components in the X and Y directions respectively; W 10 is the wind speed at 10m above the reservoir surface; A is the wind direction; α is the wind drift coefficient. Under the influence of wind, the increase in the drift speed of the oil film is about 2%-3% of the wind speed, and the drift direction forms an angle of 0 to 40° with the wind direction. Since the drift speed and direction of the center of the oil film are the sum of the velocity vectors caused by the surface flow and the wind, the wind drift coefficient α is used in the simulation as the angle correction coefficient between the wind and the flow, which is generally between 0.02-0.03, and 0.03 is taken in this paper; r is the random diffusion term, r=RE, R is a random number between 0 and 1, E is the diffusion coefficient, A α is the horizontal diffusion coefficient, provided by the hydrodynamic model; B is the random diffusion direction, B = 2πR.

[0202] In this example, the coupled oil spill module uses the particle tracking method to simulate the migration of the oil film. This method regards the oil film as consisting of a finite number of particles, and the motion trajectory of each particle is calculated by a specific formula. The particle tracking method can accurately simulate the diffusion and drift process of the oil film in a complex flow field, while taking into account the influence of multiple factors such as wind speed, wind direction, and wind drift coefficient.

[0203] In this example, the coupled oil spill module further introduces the influence of wind speed and wind direction on oil film drift. By introducing the wind drift coefficient α, the module can simulate the drift speed and direction of the oil film under the action of wind.

[0204] The value range of the wind drift coefficient α is determined based on experience and is reasonably adjusted in the simulation to ensure the accuracy of the simulation results.

[0205] In this example, the coupled oil spill module further introduces a random diffusion term. By introducing the random diffusion term r, the diffusion process of the oil film can be simulated more realistically. This random term is calculated by the random number R and the diffusion coefficient E, which can simulate the uncertainty of the oil film in the diffusion process.

[0206] The introduction of random diffusion terms makes the simulation results more consistent with the actual situation and improves the accuracy and reliability of the simulation.

[0207] In this example, the coupled oil spill module further introduces an angle correction coefficient. In the simulation, the angle correction coefficient α of wind and flow is used to correct the angle between the drift direction of the oil film and the wind direction. This correction coefficient can more accurately simulate the drift direction of the oil film under the combined action of wind and flow.

[0208] From the above, it can be seen that the coupled oil spill module in this example scheme comprehensively considers the influence of multiple factors such as flow velocity, wind speed, wind direction, wind drift coefficient, random diffusion term, etc. on the migration of oil film, so that the module can more accurately simulate the migration process of oil film in a complex environment.

[0209] In this example solution, the oil spill module is coupled, and the method of combining particle tracking with random diffusion can not only simulate the deterministic motion trajectory of the oil film, but also consider the uncertainty of the oil film in the diffusion process, making the simulation results more comprehensive and accurate, and more in line with the actual situation, providing strong support for emergency response and disposal of oil spills.

[0210] In this example, the coupled oil spill module has been reasonably adjusted in the simulation to ensure the accuracy of the simulation results. This method takes into account both the empirical value and the actual situation, making the simulation results more realistic.

[0211] The coupled oil spill module in this example scheme can be applied to different types of reservoirs and different oil leakage situations. By adjusting relevant parameters and coefficients, the module can flexibly simulate the oil film migration process under different conditions.

[0212] (3) Oil film evaporation calculation model is used to calculate the oil film evaporation.

[0213] The evaporation of oil film is affected by factors such as crude oil type, temperature, oil spill area, wind speed, solar radiation and oil film thickness. This model considers the evaporation of oil spill, and the evaporation term formula used in the program is:

[0214] F v =ln[1+B(TG / T)θexp(A-BT0 / T)][T / BT G ] (20)

[0215] Where TG is the slope of the distillation curve, T is the initial distillation point, T0 is the ambient temperature (air temperature), and A and B are dimensionless constants.

[0216] θ=KAt / V (21)

[0217] K is the air end mass transfer coefficient, A is the oil spill area, t is the time, and V is the volume of the oil spill.

[0218]

[0219] UW is the wind speed at 10 meters above the river surface, and R ranges from (0.1 to 10)×10 -4 m / s level (the K value in this paper is 10 -4 ), t is the time after the oil spill occurs, in hours.

[0220] The oil film evaporation calculation model in this example scheme comprehensively considers multiple influencing factors such as crude oil type, temperature, oil spill area, wind speed, solar radiation and oil film thickness. Multiple influencing factors are coupled together to form a comprehensive oil film evaporation calculation model, which enables the model to more accurately simulate the actual situation of oil film evaporation and improves the accuracy and reliability of the simulation.

[0221] The oil film evaporation calculation model in this example solution introduces two parameters: the slope of the distillation curve and the initial distillation point. These two parameters can reflect the evaporation characteristics of crude oil. With the introduction of these two parameters, the model can more accurately simulate the oil film evaporation process of different types of crude oil. This innovation makes the model have a wider application prospect and can be applied to different types of crude oil spills.

[0222] The oil film evaporation calculation model in this example calculates the air-end mass transfer coefficient and describes the mass transfer process between the oil film and the air by using the air-end mass transfer coefficient K. The calculation of this coefficient takes into account multiple factors such as wind speed, oil spill area and time, which can more accurately reflect the rate of oil film evaporation and improve the accuracy of the simulation.

[0223] The oil film evaporation calculation model in this example solution refines the wind speed to the wind speed of 10 meters on the river surface and gives a specific value range. This refinement allows the model to more accurately simulate the impact of wind speed on oil film evaporation.

[0224] This example uses the prediction model formed by the above scheme to predict the migration path of the oil film after oil leakage during reservoir construction. The implementation process is as follows:

[0225] As an example, a reservoir in a certain river basin is selected as the research scope. Taking the oil leakage that occurred during the construction of the reservoir as an example, the migration changes of the oil film in the future construction period of the water area are simulated and predicted based on the constructed prediction model.

[0226] First, according to the scope of the study area, the unstructured grid is drawn for the construction sea area, and the model calculation range is mainly the river section upstream and downstream of the reservoir. The model uses a curved non-orthogonal grid in the horizontal direction, with the number of grids in the horizontal direction being 978×53 and three layers in the vertical direction; among them, the horizontal grid size varies from 30 to 300m, and the grid in the reservoir construction area is locally encrypted, such as Figure 1 shown.

[0227] Next, based on the collected topographic data of the reservoir basin, the grid established in step 2 is interpolated to obtain the topographic boundary conditions of the area, such as Figure 2 shown.

[0228] Next, we collected the measured hydrodynamic data, wind conditions and solar radiation data of the reservoir basin, set up a prediction model based on the measured data, and assigned boundary values. The boundary conditions include upstream flow boundaries, wind field boundaries, and light conditions.

[0229] 5. Calibrate and verify the hydrodynamic results calculated by the prediction model based on measured data, adjust the bottom friction coefficient and other parameters of the prediction model to ensure the accuracy of the prediction model.

[0230] 6. After completing the calibration verification, perform three-dimensional hydrodynamic calculations to obtain the hydrodynamic field distribution of the reservoir area and the entire upper and lower basins, such as Figure 3 As shown, it is an example diagram of the calculated surface hydrodynamic flow field.

[0231] 7. Based on the calculation results of the flow field, according to the type of oil pollution and the amount of leakage, calculate the migration and diffusion results of the oil film in different time periods, and analyze its impact on sensitive points upstream and downstream of the reservoir area, such as the prediction results of the change of oil film thickness at the water plant intake, the position, area and maximum thickness of the oil film at different times, etc.

[0232] According to the type of oil pollution and the amount of leakage, set the relevant parameters of the oil film diffusion model, such as the initial position, area, thickness, etc. of the oil film. Import the flow field calculation results as the driving force for the oil film diffusion.

[0233] The oil film diffusion model is used in combination with the flow field calculation results to calculate the migration and diffusion of the oil film in different time periods. The position, area, and maximum thickness of the oil film in each time period are recorded.

[0234] According to the calculation results, the impact of the oil film on the sensitive points upstream and downstream of the reservoir area is analyzed, such as the prediction results of the change of the oil film thickness at the water plant water intake. Draw the oil film distribution field map, such as Figure 4As shown; show the distribution of oil film 3 hours after the reservoir oil spill accident, draw the process line of oil film thickness at the water plant water intake, such as Figure 5 As shown in Table 1, the maximum thickness of the oil film at different time periods is recorded.

[0235] 8. Based on the calculation results, formulate oil pollution interception and treatment measures.

[0236] It can be seen from the above examples that compared with the prior art, the scheme of the present invention further constructs three-dimensional refined calculations by specifically introducing the unique water flow conditions, topography and environmental factors of the reservoir, and introduces the oil film evaporation factor, thereby significantly improving the prediction accuracy and reliability, thereby more effectively dealing with the problem of oil leakage during the reservoir construction process.

[0237] Based on the above-mentioned prediction scheme for the oil film migration path after oil leakage during reservoir construction, an embodiment of the present invention also provides a computer-readable storage medium on which a program is stored. When the program is executed by a processor, the steps of the method for predicting the oil film migration path after oil leakage during reservoir construction are implemented.

[0238] An embodiment of the present invention further provides a processor, which is used to run a program, wherein when the program is run, the steps of the method for predicting the migration path of the oil film after oil leakage in reservoir construction are executed.

[0239] An embodiment of the present invention also provides a terminal device, which includes a processor, a memory, and a program stored in the memory and executable on the processor, wherein the program code is loaded and executed by the processor to implement the steps of the method for predicting the migration path of the oil film after oil leakage during reservoir construction.

[0240] The present invention also provides a computer program product, which, when executed on a data processing device, is suitable for executing the steps of the above-mentioned method for predicting the migration path of the oil film after oil leakage in reservoir construction.

[0241] In the above embodiments, the description of each embodiment has its own emphasis. For parts that are not described in detail in a certain embodiment, reference can be made to the relevant descriptions of other embodiments.

[0242] Those skilled in the art can clearly understand that, for the convenience and brevity of description, the specific working processes of the systems, devices and modules described above can refer to the corresponding processes in the aforementioned method embodiments and will not be repeated here.

[0243] Those skilled in the art will appreciate that embodiments of the present invention may be provided as methods, systems, or computer program products. Therefore, the present invention may take the form of a complete hardware embodiment, a complete software embodiment, or an embodiment combining software and hardware. Moreover, the present invention may take the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code.

[0244] The present invention is described with reference to the flowcharts and / or block diagrams of the methods, devices (systems), and computer program products of the embodiments of the present invention. It should be understood that each process and / or block in the flowchart and / or block diagram, as well as the combination of the processes and / or blocks in the flowchart and / or block diagram, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing device to generate a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate instructions for implementing the processes in the flowchart and / or block diagram. Figure 1 A process or multiple processes and / or boxes Figure 1 A device that provides the functions specified in a block or multiple blocks.

[0245] These computer program instructions may also be stored in a computer-readable memory capable of directing a computer or other programmable data processing device to operate in a specific manner, so that the instructions stored in the computer-readable memory produce an article of manufacture comprising an instruction device, which implements the process Figure 1 A process or multiple processes and / or boxes Figure 1 A function specified in one or more boxes.

[0246] These computer program instructions can also be loaded onto a computer or other programmable data processing device so that a series of operating steps are executed on the computer or other programmable device to produce a computer-implemented process, thereby providing instructions for implementing the process. Figure 1 A process or multiple processes and / or boxes Figure 1 The steps for the functions specified in one or more boxes.

[0247] In a typical configuration, a computing device includes one or more processors (CPU), input / output interfaces, network interfaces, and memory.

[0248] The memory may include non-permanent memory in a computer-readable medium, random access memory (RAM) and / or non-volatile memory in the form of read-only memory (ROM) or flash RAM. The memory is an example of a computer-readable medium.

[0249] Computer readable media include permanent and non-permanent, removable and non-removable media that can be implemented by any method or technology to store information. Information can be computer readable instructions, data structures, program modules or other data. Examples of computer storage media include, but are not limited to, phase change memory (PRAM), static random access memory (SRAM), dynamic random access memory (DRAM), other types of random access memory (RAM), read-only memory (ROM), electrically erasable programmable read-only memory (EEPROM), flash memory or other memory technology, compact disk read-only memory (CD-ROM), digital versatile disk (DVD) or other optical storage, magnetic cassettes, magnetic tape disk storage or other magnetic storage devices or any other non-transmission media that can be used to store information that can be accessed by a computing device.

[0250] It should also be noted that the terms "include", "comprises" or any other variations thereof are intended to cover non-exclusive inclusion, so that a process, method, commodity or device including a series of elements includes not only those elements, but also other elements not explicitly listed, or also includes elements inherent to such process, method, commodity or device. In the absence of more restrictions, the elements defined by the sentence "comprises a ..." do not exclude the existence of other identical elements in the process, method, commodity or device including the elements.

[0251] It will be appreciated by those skilled in the art that embodiments of the present invention may be provided as methods, systems or computer program products. Therefore, the present invention may take the form of a complete hardware embodiment, a complete software embodiment or an embodiment combining software and hardware. Moreover, the present invention may take the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) containing computer-usable program codes.

[0252] The above-mentioned method of the present invention, or a specific system unit, or a part of the same, is a pure software architecture, which can be arranged on a physical medium, such as a hard disk, an optical disk, or any electronic device (such as a smart phone, a computer-readable storage medium) through program code. When the machine loads the program code and executes it (such as a smart phone loading and executing it), the machine becomes a device for implementing the present invention. The above-mentioned method and device of the present invention can also be transmitted in the form of program code through some transmission media, such as cables, optical fibers, or any transmission mode. When the program code is received, loaded and executed by a machine (such as a smart phone), the machine becomes a device for implementing the present invention.

[0253] The above shows and describes the basic principles, main features and advantages of the present invention. It should be understood by those skilled in the art that the present invention is not limited to the above embodiments, and the above embodiments and descriptions are only for explaining the principles of the present invention. Without departing from the spirit and scope of the present invention, the present invention may have various changes and improvements, which fall within the scope of the present invention to be protected. The scope of protection of the present invention is defined by the attached claims and their equivalents.

Claims

1. A method for predicting the migration path of oil film after oil leakage in reservoir construction, characterized in that: The prediction method comprises: (1) Construct a three-dimensional reservoir hydrodynamic model; On the basis of the three-dimensional hydrodynamic model, the three-dimensional topographic data of the reservoir and the basic water conservancy data are introduced to form a three-dimensional reservoir hydrodynamic model, and the real-time monitoring data of the reservoir is integrated into the three-dimensional reservoir hydrodynamic model to perform real-time correction on the three-dimensional reservoir hydrodynamic model; (2) calculating the flow field of the reservoir area based on the three-dimensional reservoir hydrodynamic model constructed in step (1); (3) coupled oil spill module; Based on the historical data of oil spill accidents and the corresponding physical and chemical properties of oil spilled substances, an oil spill module is constructed, and the constructed oil spill module is coupled with the three-dimensional reservoir hydrodynamic model. By exchanging the reservoir flow field data calculated by the three-dimensional reservoir hydrodynamic model in real time with the corresponding oil spill data information in the oil spill module, the oil spill module and the three-dimensional reservoir hydrodynamic model work together to jointly simulate the movement trajectory and fate of the oil spill in the water. (4) Integrate the oil film evaporation calculation module; Based on the oil spill accident data, the physical and chemical properties of the corresponding oil film and the corresponding meteorological data, an oil film evaporation calculation module is constructed, and the constructed oil film evaporation calculation module is integrated with the oil spill module and the three-dimensional reservoir hydrodynamic model. By exchanging the corresponding oil spill data information in the oil spill module with the corresponding calculated oil film evaporation data information in the oil film evaporation calculation module in real time, the collaborative work of the two modules is realized, thereby simulating the evaporation process of the oil film.

2. The method for predicting the migration path of oil film after oil leakage in reservoir construction according to claim 1 is characterized in that: The prediction method is configured to use a particle tracking method to simulate the convection, expansion and turbulent diffusion process of the oil spill in the oil spill module constructed in step (3), and a random walk method to simulate the turbulent diffusion process of the oil spill in the horizontal direction.

3. The method for predicting the migration path of oil film after oil leakage in reservoir construction according to claim 1 is characterized in that: The oil film evaporation calculation module constructed in step (4) of the prediction method is configured to simulate the evaporation process of the oil film based on the mass conservation and energy conservation modes.

4. A prediction system for the migration path of oil film after oil leakage during reservoir construction, characterized in that: The prediction system comprises: A three-dimensional reservoir hydrodynamic model, wherein the three-dimensional reservoir hydrodynamic model is formed by introducing three-dimensional reservoir terrain data and water conservancy basic data on the basis of the three-dimensional hydrodynamic model, and the three-dimensional reservoir hydrodynamic model is configured to calculate the flow field of the reservoir area; An oil spill module, which is formed based on historical data of oil spill accidents and corresponding physical and chemical property data of oil spill substances, and is configured to be coupled with a three-dimensional reservoir hydrodynamic model, so that the oil spill module and the three-dimensional reservoir hydrodynamic model can work together to jointly simulate the movement trajectory and fate of the oil spill in the water by exchanging reservoir area flow field data calculated by the three-dimensional reservoir hydrodynamic model and corresponding oil spill data information in the oil spill module in real time; The oil film evaporation calculation module is composed of oil spill accident data, physical and chemical property data of the corresponding oil film and corresponding meteorological data, and is configured to be integrated with the oil spill module and the three-dimensional reservoir hydrodynamic model. The corresponding oil spill data information in the oil spill module and the corresponding oil film evaporation data information calculated and determined in the oil film evaporation calculation module are exchanged in real time to achieve the collaborative work of the two modules, thereby simulating the evaporation process of the oil film.

5. The system for predicting the migration path of oil film after oil leakage in reservoir construction according to claim 4 is characterized in that: The oil spill module is configured to use a particle tracking method to simulate the convection, expansion and turbulent diffusion process of the oil spill, and a random walk method to simulate the turbulent diffusion process of the oil spill in the horizontal direction.

6. The system for predicting the migration path of oil film after oil leakage in reservoir construction according to claim 4, characterized in that: The oil film evaporation calculation module is configured to simulate the evaporation process of the oil film based on mass conservation and energy conservation modes.

7. A computer-readable storage medium having a program stored thereon, characterized in that: When the program is executed by a processor, the steps of the method for predicting the migration path of an oil film after oil leakage during reservoir construction described in any one of claims 1 to 3 are implemented.

8. A processor, the processor being used to run a program, characterized in that: When the program is running, the steps of the method for predicting the migration path of the oil film after oil leakage in reservoir construction described in any one of claims 1 to 3 are executed.

9. A terminal device, comprising a processor, a memory, and a program stored in the memory and executable on the processor, characterized in that: The program code is loaded and executed by the processor to implement the steps of the method for predicting the migration path of oil film after oil leakage in reservoir construction as described in any one of claims 1-3.

10. A computer program product, characterized in that When executed on a data processing device, it is suitable for executing the steps of the method for predicting the migration path of an oil film after oil leakage in reservoir construction as described in any one of claims 1 to 3.

Images