Plate exchanger blockage cause prediction method, device, equipment, medium and product

Abstract

The application discloses a plate exchanger blockage cause prediction method, device, equipment, medium and product, relates to the coastal engineering technical field, and the method comprises the following steps: a three-dimensional baroclinic ocean numerical model in a preset space range where a target engineering sea area is located is constructed; the target engineering sea area comprises a plurality of plate exchangers; a calculation domain is determined according to boundary conditions of the three-dimensional baroclinic ocean numerical model; the boundary conditions comprise surface boundary conditions, bottom boundary conditions and fixed boundary conditions; the three-dimensional baroclinic ocean numerical model is configured based on the calculation domain; water body structure in the calculation domain is simulated by running the configured three-dimensional baroclinic ocean numerical model, and observation data in the target engineering sea area in a target time period is obtained; the observation data comprises tide level, flow velocity and flow direction of each observation point; and the cause of plate exchanger blockage is determined according to the observation data. The application can provide reliable data support for blockage cause analysis.

Description

Technical Field

[0001] This application relates to the field of coastal engineering technology, and in particular to a method, apparatus, equipment, medium and product for predicting the causes of blockage in plate heat exchangers. Background Technology

[0002] In coastal engineering projects, plate heat exchangers (Sea Water Cooling Heat Exchangers, SECs) often experience pressure fluctuations due to biological blockage, affecting equipment operating efficiency and stability. However, the causes of blockage in plate heat exchangers are lacking analysis, making it impossible to implement targeted prevention measures. Summary of the Invention

[0003] The purpose of this application is to provide a method, apparatus, equipment, medium, and product for predicting the causes of blockage in plate heat exchangers, which can provide reliable data support for the analysis of blockage causes.

[0004] To achieve the above objectives, this application provides the following solution:

[0005] Firstly, this application provides a method for predicting the causes of blockage in a plate switch, the method comprising:

[0006] A three-dimensional baroclinic marine numerical model is constructed within a predetermined spatial range of the target engineering sea area; the target engineering sea area includes multiple plate exchangers;

[0007] The computational domain is determined based on the boundary conditions of the three-dimensional baroclinic ocean numerical model; the boundary conditions include surface boundary conditions, bottom boundary conditions, and solid boundary conditions.

[0008] The three-dimensional baroclinic ocean numerical model is configured based on the computational domain;

[0009] The configured three-dimensional baroclinic ocean numerical model is used to simulate the water structure within the computational domain, and the observation data of the target engineering sea area within the target time period is obtained; the observation data includes the tide level, current velocity and current direction of each observation point.

[0010] The cause of the plate exchanger blockage was determined based on the observed data.

[0011] Optionally, configuring the three-dimensional baroclinic ocean numerical model based on the computational domain specifically includes:

[0012] Based on the computational domain, initial conditions and open boundary conditions are input into the three-dimensional baroclinic ocean numerical model;

[0013] Input atmospheric boundary conditions based on the initial conditions;

[0014] The bottom friction and horizontal eddy viscosity coefficients of the three-dimensional baroclinic ocean numerical model are determined based on the atmospheric boundary conditions.

[0015] The time steps of the inner and outer models of the three-dimensional baroclinic ocean numerical model are determined based on the bottom friction and the horizontal eddy viscosity coefficient.

[0016] Optionally, the initial condition is a temperature-salinity field; the temperature-salinity field is derived from HYCOM reanalysis data;

[0017] The open boundary conditions include applying tidal level and temperature-salinity driving data; the temperature-salinity driving data is derived from HYCOM reanalysis data, and the frequency of applying tidal level and temperature-salinity driving data is once every 3 hours.

[0018] The atmospheric boundary conditions are wind field data at a height of 10 meters above the sea surface per hour, which are derived from the ECMWF 5th generation reanalysis dataset.

[0019] Optionally, the surface boundary conditions include wind pressure components and net heat flux; the bottom boundary conditions include bottom pressure components and bottom groundwater flux; and the solid boundary conditions include the boundary normal velocity components, precipitation rate, and evaporation rate.

[0020] Optionally, the governing equations of the three-dimensional baroclinic ocean numerical model include the continuity equation, momentum equation, density equation, temperature equation, and salinity equation.

[0021] Optionally, the cause of plate heat exchanger blockage can be determined based on the observed data, specifically including:

[0022] When the flow velocity at the front end of the plate exchanger's intake is higher than the set value at the rear end of the intake, it is determined that the plate exchanger is blocked because sediment is easily drawn into the intake culvert.

[0023] Secondly, this application provides a plate heat exchanger blockage cause prediction device, characterized in that the plate heat exchanger blockage cause prediction device applies any one of the plate heat exchanger blockage cause prediction methods, and the plate heat exchanger blockage cause prediction device includes:

[0024] A three-dimensional baroclinic ocean numerical model construction module is used to construct a three-dimensional baroclinic ocean numerical model within a preset spatial range where the target engineering sea area is located; the target engineering sea area includes multiple plate exchangers;

[0025] The computational domain determination module is used to determine the computational domain based on the boundary conditions of the three-dimensional baroclinic ocean numerical model; the boundary conditions include surface boundary conditions, bottom boundary conditions, and solid boundary conditions.

[0026] A configuration module is used to configure the three-dimensional baroclinic ocean numerical model based on the computational domain;

[0027] The observation data simulation module is used to run the configured three-dimensional baroclinic ocean numerical model to obtain observation data of the target engineering sea area within the target time period.

[0028] The cause determination module is used to determine the cause of the plate exchanger blockage based on the observation data.

[0029] Thirdly, this application provides a computer device, including: a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the computer program to implement the steps of the plate switch congestion cause prediction method described in any one of the above.

[0030] Fourthly, this application provides a computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the steps of the board switch congestion cause prediction method described above.

[0031] Fifthly, this application provides a computer program product, including a computer program that, when executed by a processor, implements the steps of the board switch blockage cause prediction method described above.

[0032] According to the specific embodiments provided in this application, the following technical effects are disclosed:

[0033] This application provides a method, apparatus, equipment, medium, and product for predicting the causes of plate exchanger blockage. When plate exchanger blockage occurs, a configured three-dimensional baroclinic ocean numerical model (FVCOM) is run to simulate the water structure within the computational domain, obtaining observational data of the target engineering sea area within a target time period. The observational data includes tidal level, current velocity, and current direction at each observation point. Based on the observational data, the cause of plate exchanger blockage is determined. Since the three-dimensional baroclinic ocean numerical model is a model based on an unstructured triangular mesh, it can accurately simulate the tidal current movement patterns in areas with complex coastlines and topography, providing reliable data support for blockage cause analysis. Attached Figure Description

[0034] To more clearly illustrate the technical solutions in the embodiments of this application or the prior art, the drawings used in the embodiments will be briefly introduced below. Obviously, the drawings described below are only some embodiments of this application. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0035] Figure 1A flowchart illustrating a method for predicting the causes of blockage in a plate switch, provided in an embodiment of this application;

[0036] Figure 2 This application provides a computational domain grid and water depth distribution map as an embodiment.

[0037] Figure 3 This is a schematic diagram of the engineering area grid and terrain provided in one embodiment of this application;

[0038] Figure 4 This is a schematic diagram of the tide level verification results provided in an embodiment of this application;

[0039] Figure 5 A schematic diagram comparing the calculated and measured values ​​of the surface current velocity at the first observation point of the spring tide, provided in an embodiment of this application;

[0040] Figure 6 A schematic diagram comparing the calculated and measured values ​​of the mid-tidal current velocity at a first observation point provided in an embodiment of this application;

[0041] Figure 7 A schematic diagram comparing the calculated and measured values ​​of the bottom current velocity and direction of the tidal bore at a first observation point, provided in an embodiment of this application;

[0042] Figure 8 A schematic diagram comparing the calculated and measured values ​​of the surface current velocity and direction of the high tide at a second observation point, provided in an embodiment of this application;

[0043] Figure 9 A schematic diagram comparing the calculated and measured values ​​of the mid-tidal current velocity and direction at a second observation point provided in an embodiment of this application;

[0044] Figure 10 A schematic diagram comparing the calculated and measured values ​​of the bottom current velocity and direction of the tidal bore at a second observation point, provided in an embodiment of this application;

[0045] Figure 11 A schematic diagram comparing the calculated and measured values ​​of the surface current velocity and direction of the high tide at a third observation point, provided in an embodiment of this application;

[0046] Figure 12 A schematic diagram comparing the calculated and measured values ​​of the mid-tidal current velocity and direction at a third observation point provided in an embodiment of this application;

[0047] Figure 13 A schematic diagram comparing the calculated and measured values ​​of the bottom current velocity and direction of the tidal bore at a third observation point provided in an embodiment of this application;

[0048] Figure 14 This is a schematic diagram comparing the calculated and measured values ​​of the neap tide surface velocity at a first observation point, provided in an embodiment of this application.

[0049] Figure 15 A schematic diagram comparing the calculated and measured values ​​of the mid-water velocity at the first observation point during neap tides, provided in an embodiment of this application;

[0050] Figure 16 A schematic diagram comparing the calculated and measured values ​​of the bottom current velocity and direction of neap tide at a first observation point, provided in an embodiment of this application;

[0051] Figure 17 A schematic diagram comparing the calculated and measured values ​​of the surface velocity and direction of the neap tide at a second observation point, provided in an embodiment of this application;

[0052] Figure 18 A schematic diagram comparing the calculated and measured values ​​of mid-level flow velocity and direction at a second observation point during neap tides, provided in an embodiment of this application.

[0053] Figure 19 This is a schematic diagram comparing the calculated and measured values ​​of the bottom current velocity and direction of the neap tide at a second observation point, provided in an embodiment of this application.

[0054] Figure 20 A schematic diagram comparing the calculated and measured values ​​of the surface velocity and direction of neap tide at a third observation point provided in an embodiment of this application;

[0055] Figure 21 A schematic diagram comparing the calculated and measured values ​​of the mid-level flow velocity and direction at the third observation point during neap tides, provided in an embodiment of this application;

[0056] Figure 22 A schematic diagram comparing the calculated and measured values ​​of the bottom current velocity and direction of neap tide at a third observation point provided in an embodiment of this application;

[0057] Figure 23 A large-area surface turbulence flow field diagram provided for an embodiment of this application;

[0058] Figure 24 A large-area bottom-layer rapid flow field diagram provided in one embodiment of this application;

[0059] Figure 25 A large-area surface jet stream field diagram provided for an embodiment of this application;

[0060] Figure 26 A large-area bottom jet flow field diagram provided in an embodiment of this application;

[0061] Figure 27 A vertical average flow field diagram of the sea area of ​​the intake port basin during a rapid rise, provided in one embodiment of this application;

[0062] Figure 28 This application provides a vertical average flow field diagram of the sea area near the intake during a rapid rise in water level, as shown in one embodiment of the present application.

[0063] Figure 29 A vertical average flow field diagram of the sea area of ​​the intake port basin at the moment of rapid descent, provided in an embodiment of this application;

[0064] Figure 30 A vertical average flow field diagram of the sea area near the water intake at the moment of rapid descent, provided in one embodiment of this application;

[0065] Figure 31 A 15-day average vertical average flow field diagram of the intake port basin sea area provided in an embodiment of this application;

[0066] Figure 32 A 15-day average vertical mean flow field diagram of the sea area near the water intake provided in an embodiment of this application;

[0067] Figure 33 This is a schematic diagram of the functional modules of a plate switch blockage cause prediction device provided in an embodiment of this application.

[0068] Figure 34 This is a schematic diagram of the structure of a computer device provided in an embodiment of this application. Detailed Implementation

[0069] The technical solutions of the embodiments of this application will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of this application, and not all embodiments. Based on the embodiments of this application, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of this application.

[0070] To make the above-mentioned objectives, features and advantages of this application more apparent and understandable, the application will be further described in detail below with reference to the accompanying drawings and specific embodiments.

[0071] In one exemplary embodiment, this application provides a method for predicting the causes of blockage in a plate switch, such as... Figure 1 As shown, the method for predicting the causes of blockage in a plate heat exchanger includes:

[0072] Step 101: Construct a three-dimensional baroclinic marine numerical model within a preset spatial range of the target engineering sea area; the target engineering sea area includes multiple plate exchangers.

[0073] Step 102: Determine the computational domain based on the boundary conditions of the three-dimensional baroclinic ocean numerical model; the boundary conditions include surface boundary conditions, bottom boundary conditions, and solid boundary conditions.

[0074] Step 103: Configure the three-dimensional baroclinic ocean numerical model based on the computational domain.

[0075] Step 104: Run the configured three-dimensional baroclinic ocean numerical model to simulate the water structure within the computational domain and obtain observation data of the target engineering sea area within the target time period; the observation data includes the tide level, current velocity and current direction of each observation point.

[0076] Step 105: Determine the cause of the plate exchanger blockage based on the observed data.

[0077] In one exemplary embodiment, the three-dimensional baroclinic ocean numerical model is a three-dimensional baroclinic ocean numerical model based on an unstructured triangular mesh.

[0078] The governing equations of the three-dimensional baroclinic ocean numerical model include the continuity equation, momentum equation, density equation, temperature equation, and salinity equation.

[0079] By transforming the vertical coordinates, the vertical coordinates are converted into coordinates of the relative positions of the free surface of the seawater and the seabed, in order to better fit the seabed topography.

[0080] Continuity equation: .

[0081] Momentum equation:

[0082] ;

[0083] .

[0084] In the formula: t For time; x , y , z The coordinates are due east, due north, and vertical. u , v , w These represent the computational fluid dynamics unit in... x , y , z The velocity component in the direction; f Coriolis parameter; F u , F v They represent the horizontal momentum at... x , y Components in direction; K m The viscosity coefficient is the coefficient of vertical rotation. P For pressure, For reference density, it is usually taken as the average density of seawater, about 1025 kg / m³.

[0085] The formula for calculating pressure is: .

[0086] In the formula:g This is the acceleration due to gravity.

[0087] Density equation: .

[0088] It is generally assumed that water is incompressible and that its density is only related to temperature. T and salinity S Simple linear correlation.

[0089] Temperature equation: .

[0090] Salinity equation: .

[0091] In the formula: S Salinity; T For temperature; K h This is the vertical coefficient term for rotational dissipation; F S This is the dissipation term for salinity; F T This is the dissipation term for temperature.

[0092] FVCOM via vertical Coordinate transformation converts vertical coordinates into coordinates of the relative positions of the free surface of the seawater and the seabed, which can better fit the seabed topography. The coordinate transformation formula is:

[0093] .

[0094] in, D For water depth; H Still water is deep; ζ The height of the free surface; The range of change is from -1 to 0, which can be obtained by substituting the continuity equation and momentum equation from the previous text.

[0095] Continuity equation:

[0096] .

[0097] Momentum equation:

[0098]

[0099] ;

[0100]

[0101] .

[0102] in, Let be the integration variable, representing the vertical coordinate that changes during the integration process.

[0103] Temperature equation:

[0104] .

[0105] Salinity equation:

[0106] .

[0107] The horizontal diffusion term is defined as follows:

[0108] ;

[0109] ;

[0110] ;

[0111] In the formula: A m The horizontal eddy diffusion coefficient; A h The horizontal thermal diffusivity; K m The vertical eddy viscosity coefficient; K h is the thermal vertical eddy diffusion coefficient; F u and F v It is horizontal momentum; F T and F S For thermal and salinity diffusion terms; The gradient of shortwave radiation; q , l These are turbulent kinetic energy and turbulent scale, respectively.

[0112] In an exemplary embodiment, the surface boundary conditions include wind pressure components and net heat flux; the bottom boundary conditions include bottom pressure components and bottom groundwater flux; and the solid boundary conditions include the boundary normal velocity components, precipitation rate, and evaporation rate.

[0113] The boundary conditions are given below, at the surface layer, =0:

[0114] ;

[0115] ;

[0116] ;

[0117] At the bottom boundary, =-1:

[0118] ;

[0119] ;

[0120] ;

[0121] Fixed boundary: .

[0122] in, These are the surface wind pressure components in the x and y directions, respectively; These are the base pressure components in the x and y directions, respectively; Q b Ω represents the bottom groundwater flux; Ω represents the area of ​​the groundwater source. V n The normal velocity component is the boundary. Net heat flux at the surface; For surface short-wave heat flux; c p α represents the specific heat capacity of seawater; α represents the seabed slope; n represents the direction of the outer normal to the boundary. Precipitation rate; Evaporation rate This represents the horizontal eddy diffusion coefficient.

[0123] In one exemplary embodiment, a computational domain is established based on the characteristics of the engineering area. A triangular mesh system is employed, with larger meshes used in areas farther from the project and smaller meshes used near the project. More specifically, areas farther from the target engineering sea area use a first-size mesh (larger mesh), while areas closer to the target engineering sea area use a second-size mesh (smaller mesh). The entire simulation area (computational domain) consists of 17,982 nodes and 34,143 triangular elements, with a minimum spatial step size of approximately 5 meters.

[0124] In this example, the computational domain of the three-dimensional baroclinic ocean numerical model and the grid distribution within the computational domain are as follows: Figure 2 and Figure 3 As shown. In order to clearly understand the tidal conditions in the sea area near the target project, the sea area near the intake and discharge outlets was densified.

[0125] Water depth and shoreline were obtained from nautical charts, while local water depths near the project area were obtained from measured topographic data.

[0126] In an exemplary embodiment, step 103 specifically includes steps 301-304.

[0127] Step 301: Based on the computational domain, the initial conditions and open boundary conditions are input into the three-dimensional baroclinic ocean numerical model using the hot-start calculation method.

[0128] The initial condition is a temperature-salinity field, which is derived from HYCOM (HYBRID COORDINATEOCEAN MODEL) reanalysis data. For example, the temperature-salinity field from HYCOM reanalysis data on May 23, 2024, is used as the initial field, and the three-dimensional baroclinic ocean numerical model is a baroclinic model, considering the influence of temperature and salinity effects on the water density field to better simulate the water structure of the engineering sea area.

[0129] There are three open boundaries within the computational domain (e.g., Figure 2 (The three sea areas in the central, northern, and southwestern regions), the open boundary conditions include the application of tidal level and temperature-salinity driving data; the temperature-salinity driving data is derived from HYCOM reanalysis data, and the frequency of applying tidal level and temperature-salinity driving data is once every 3 hours. The open boundary tidal level is provided by the large model.

[0130] Step 302: Input atmospheric boundary conditions based on the initial conditions.

[0131] The atmospheric boundary condition is the hourly wind field data at 10 meters above sea level, which is derived from the ECMWF (European Centre for Medium-Range Weather Forecasts) Generation 5 reanalysis dataset (ERA5).

[0132] Step 303: Determine the bottom friction and horizontal eddy viscosity coefficients of the three-dimensional baroclinic ocean numerical model based on the atmospheric boundary conditions.

[0133] Specifically, step 303 includes: based on the FVCOM's drag coefficient. C d The expression for the bottom friction is given. The horizontal eddy viscosity coefficient Am is calculated using the Smagorinsky parametric model.

[0134] FVCOM's undercarriage drag coefficient C d The expression is:

[0135] .

[0136] In the formula: k=0.4, which is the KAMAN constant; z0 is the seabed roughness; Z ab It is the distance between the σ layer closest to the seabed and the seabed, and 0.0025 is the default minimum value.

[0137] In FVCOM, the horizontal eddy viscosity coefficient A m Calculated using the Smagorinsky parametric model:

[0138] .

[0139] In the formula: C is a constant parameter, The area of ​​a single momentum control element.

[0140] Step 304: Determine the inner model time step and outer model time step of the three-dimensional baroclinic ocean numerical model based on the bottom friction and the horizontal eddy viscosity coefficient.

[0141] Specifically, step 304 includes: calculating the stability conditions based on CFL (Courant–Friedrichs–Lewy), setting the time step of the inner model to 2s, the time step of the outer model to 0.2s, and the ratio of the inner model to the outer model to 10:1.

[0142] In one exemplary embodiment, the hydrodynamic model is validated based on the internal model time step and the external model time step. For example, tidal level, current velocity, and current direction data obtained during the spring and neap tide observations in June 2024 are used to validate the constructed three-dimensional baroclinic ocean numerical model. The measured and calculated results are compared to ensure that the model (the Hongyanhe three-dimensional hydrodynamic model constructed based on the FVCOM model) can well reflect the tidal current movement patterns of the engineering sea area.

[0143] The tidal current data used in the three-dimensional baroclinic ocean numerical model simulation were obtained during the spring and neap tide observations in June 2024, including tidal level, current velocity, and current direction data from three consecutive fixed-point current velocity and current direction observations, as well as tidal level data from a temporary tide gauge station. The three consecutive fixed-point observations were the first observation point, the second observation point, and the third observation point.

[0144] A comparison of calculated (simulated) and measured tide levels during spring and neap tide observations at three observation points is shown below. Figure 4 As shown, the comparison between calculated and measured values ​​of flow velocity and flow direction is as follows: Figures 5-22 As shown. Figures 5-22 Part (a) shows the results of the flow velocity comparison, and part (b) shows the results of the flow direction comparison.

[0145] A comparison of the measured and calculated results shows that, except for a few points where the calculated values ​​deviate significantly from the measured values, the calculated changes in tidal level, flow velocity, and flow direction at each observation point are basically consistent with the measured values, and meet the requirements of relevant regulations and engineering needs. A comparison of the numerical simulation results indicates that the three-dimensional tidal current mathematical model established in this application can well reflect the tidal current movement patterns in the engineering sea area, and can further provide necessary hydrodynamic conditions for analyzing and predicting the impact of thermal drainage generated by the engineering construction.

[0146] Figures 23-26 These are the high-current fields (surface and bottom layers) and low-current fields (surface and bottom layers) in the large-scale area and the sea area near the project. The tidal currents in the outer sea area of ​​the Hongyan River are generally NE-SW oriented. During periods of high tide, the flood current flows from south to north along the western coast of Changxing Island into Fuzhou Bay, continuing northward through the Hongyan River area where the project is located. In the project area, the flood current direction is slightly NE, and the flow direction is generally consistent between the surface and bottom layers. In terms of velocity distribution, the surface current has a higher velocity, while the bottom current has a lower velocity. During periods of low tide, the ebb current flows from north to south along the southeast coast of Liaodong Bay out of the bay, flowing through the Hongyan River area where the project is located. In the project area, the ebb current direction is slightly SW, and the flow direction is generally consistent between the surface and bottom layers. The project sea area is dominated by significant reciprocating currents.

[0147] The maximum surface current velocity in the vicinity of the project area can reach 1.5 m / s during high and low tides. In addition, in the sea area southwest of the nuclear power plant, due to the topography, the current velocity is significantly higher than in other sea areas, with a maximum current velocity of up to 1.8 m / s.

[0148] This application also includes an analysis of the flow field characteristics of the engineering sea area based on the verification results of a three-dimensional baroclinic ocean numerical model. Here, the verification results refer to the comparison between the simulation output of the numerical model and the actual data observed in the engineering sea area. The analysis covers the large-scale flow field characteristics of the engineering sea area, including the flow field conditions during periods of rapid rise and fall, as well as the flow field characteristics within the harbor basin.

[0149] This application also includes simulations of blockage causes based on flow field characteristics. For example, it simulates the causes of biological and residual blockage in plate exchangers and analyzes the relationship between flow velocity near the intake and the transport of blockage material. The simulation revealed that the flow velocity on the intake side is significantly higher than that of other units, resulting in an increased sediment-carrying capacity of the water, making it easier for sediment to be entrained into the intake culvert.

[0150] Based on the simulation results, blockage prevention measures to eliminate eddies in the generating unit are proposed. The causes of eddy formation on both sides of the water intake are analyzed, and based on the site conditions, tidal flow numerical simulation comparative analysis is used to propose site-specific suggestions for eliminating eddies.

[0151] In an exemplary embodiment, determining the cause of plate exchanger blockage based on the observation data specifically includes: when the flow velocity at the front end of the plate exchanger's intake is higher than the set value of the flow velocity at the rear end of the intake, it is determined that the cause of the plate exchanger blockage is that sediment is entrained into the intake culvert.

[0152] Analysis of the flow field characteristics inside the harbor basin: For example, the study investigated the causes of biological and debris blockage in plate exchangers, and presented vertically averaged flow field diagrams of the harbor basin and the sea area near the intake during periods of rapid rise and fall. Figures 27-30During periods of rapid rise, the average vertical velocity near the intake reached 0.25 m / s, 66.7% higher than the 0.15 m / s of Unit 6. During periods of rapid fall, the velocities increased to 0.30 m / s and 0.20 m / s respectively, with the difference remaining at 50%. At these two extreme moments, the velocity near the unit's intake was higher than that of other units. Long-period average flow field data from 15 days of continuous simulation shows that... Figures 31-32 The current velocity near the water intake reached 0.28 m / s, while the current velocity in other areas was about 0.18 m / s.

[0153] Cause analysis: The flow velocity at the intake of Unit 5 was significantly higher than that of other units (66.7% during rapid rises and 50% during rapid falls). The other units refer to multiple plate heat exchangers. This increased flow velocity exponentially increases the water's sediment-carrying capacity, making it easier for debris such as shell fragments (particle size > 0.5 mm) to be transported to the intake. Actual measurement data also shows that the amount of shell and other sediments on the seabed at the intake was significantly higher than at other units, confirming the existence of this transport mechanism.

[0154] Based on the analysis of the causes of congestion, congestion prevention and control measures are proposed.

[0155] Observations revealed eddies forming on both sides of the intakes of Phase II (Units 5 and 6), while none were observed in Phase I. The reason for this is that both sides of the Phase II intakes have an inner bay, which facilitates the formation of circulation currents. The Phase I intakes, being planar, disperse ocean currents and do not form eddies. The inner bays are likely one of the main causes of eddy formation. Therefore, based on the site conditions and comparative analysis of tidal current models, tailored recommendations for eliminating eddies were proposed.

[0156] This application constructs a three-dimensional baroclinic marine numerical model, combined with high-precision flow field simulation and boundary condition settings, to accurately analyze the ocean current field structure and blockage transport trajectory under different tidal conditions within the harbor basin. This method is highly accurate, adaptable, and eco-friendly, and is also applicable to the ecological optimization and maintenance of coastal engineering projects.

[0157] By simulating the transport trajectory of blockages, this study reveals the hydrodynamic causes of blockages in plate heat exchangers and clarifies the impact of velocity differences and eddy formation on blockage input. This method is applicable to various coastal engineering scenarios, allowing for adjustments to model parameters and boundary conditions based on specific engineering needs, demonstrating broad applicability and flexibility.

[0158] By optimizing the intake design and flow field structure, the risk of plate heat exchanger blockage is significantly reduced, equipment operating efficiency and stability are improved, and maintenance costs are reduced.

[0159] This application analyzes the blockage of plate heat exchangers from an ecosystem perspective, based on factors such as water flow, organisms, and the distribution of attached organisms, and utilizes natural forces such as water flow to resolve the blockage.

[0160] Based on the same inventive concept, this application also provides a plate switch blockage cause prediction device for implementing the plate switch blockage cause prediction method described above. The solution provided by this device is similar to the solution described in the above method. Therefore, the specific limitations of one or more plate switch blockage cause prediction device embodiments provided below can be found in the limitations of the plate switch blockage cause prediction method above, and will not be repeated here.

[0161] In one exemplary embodiment, such as Figure 33 As shown, a plate heat exchanger blockage cause prediction device is provided. The plate heat exchanger blockage cause prediction device applies the aforementioned plate heat exchanger blockage cause prediction method. The plate heat exchanger blockage cause prediction device includes:

[0162] A three-dimensional baroclinic ocean numerical model construction module is used to construct a three-dimensional baroclinic ocean numerical model within a preset spatial range where the target engineering sea area is located; the target engineering sea area includes multiple plate exchangers;

[0163] The computational domain determination module is used to determine the computational domain based on the boundary conditions of the three-dimensional baroclinic ocean numerical model; the boundary conditions include surface boundary conditions, bottom boundary conditions, and solid boundary conditions.

[0164] A configuration module is used to configure the three-dimensional baroclinic ocean numerical model based on the computational domain;

[0165] The observation data simulation module is used to run the configured three-dimensional baroclinic ocean numerical model to obtain observation data of the target engineering sea area within the target time period.

[0166] The cause determination module is used to determine the cause of the plate exchanger blockage based on the observation data.

[0167] In one exemplary embodiment, a computer device is provided, which may be a server or a terminal, and its internal structure diagram may be as follows. Figure 34As shown, this computer device includes a processor, memory, input / output (I / O) interfaces, and a communication interface. The processor, memory, and I / O interfaces are connected via a system bus, and the communication interface is also connected to the system bus via the I / O interfaces. The processor provides computational and control capabilities. The memory includes non-volatile storage media and internal memory. The non-volatile storage media stores the operating system, computer programs, and a database. The internal memory provides the environment for the operating system and computer programs in the non-volatile storage media to run. The database stores data for predicting the causes of board switch congestion. The I / O interfaces are used for exchanging information between the processor and external devices. The communication interface is used for communicating with external terminals via a network connection. When the computer program is executed by the processor, it implements a method for predicting the causes of board switch congestion.

[0168] Those skilled in the art will understand that Figure 34 The structures shown are merely block diagrams of some structures related to the present application and do not constitute a limitation on the computer device to which the present application is applied. Specific computer devices may include more or fewer components than shown in the figures, or combine certain components, or have different component arrangements. In an exemplary embodiment, a computer device is provided, including a memory and a processor. The memory stores a computer program, and the processor executes the computer program to implement the steps in the above-described method embodiments.

[0169] In one exemplary embodiment, a computer-readable storage medium is provided storing a computer program that, when executed by a processor, implements the steps in the above-described method embodiments.

[0170] In one exemplary embodiment, a computer program product is provided, including a computer program that, when executed by a processor, implements the steps in the above-described method embodiments.

[0171] It should be noted that the user information (including but not limited to user device information, user personal information, etc.) and data (including but not limited to data used for analysis, data stored, data displayed, etc.) involved in this application are all information and data authorized by the user or fully authorized by all parties, and the collection, use and processing of the relevant data must comply with relevant regulations.

[0172] Those skilled in the art will understand that all or part of the processes in the above embodiments can be implemented by a computer program instructing related hardware. The computer program can be stored in a non-volatile computer-readable storage medium. When executed, the computer program can include the processes of the embodiments described above. Any references to memory, databases, or other media used in the embodiments provided in this application can include at least one of non-volatile and volatile memory. Non-volatile memory can include read-only memory (ROM), magnetic tape, floppy disk, flash memory, optical memory, high-density embedded non-volatile memory, resistive random access memory (ReRAM), magnetic random access memory (MRAM), ferroelectric random access memory (FRAM), phase change memory (PCM), graphene memory, etc. Volatile memory can include random access memory (RAM) or external cache memory, etc. By way of illustration and not limitation, RAM can take many forms, such as Static Random Access Memory (SRAM) or Dynamic Random Access Memory (DRAM).

[0173] The databases involved in the embodiments provided in this application may include at least one type of relational database and non-relational database. Non-relational databases may include, but are not limited to, blockchain-based distributed databases. The processors involved in the embodiments provided in this application may be general-purpose processors, central processing units, graphics processing units, digital signal processors, programmable logic units, data processing logic units, etc., and are not limited to these.

[0174] The technical features of the above embodiments can be combined in any way. For the sake of brevity, not all possible combinations of the technical features in the above embodiments are described. However, as long as there is no contradiction in the combination of these technical features, they should be considered to be within the scope of this specification.

[0175] This document uses specific examples to illustrate the principles and implementation methods of this application. The descriptions of the above embodiments are only for the purpose of helping to understand the methods and core ideas of this application. Furthermore, those skilled in the art will recognize that, based on the ideas of this application, there will be changes in the specific implementation methods and application scope. Therefore, the content of this specification should not be construed as a limitation of this application.

Claims

1. A method for predicting the causes of blockage in a plate heat exchanger, characterized in that, The method for predicting the causes of blockage in plate switches includes: A three-dimensional baroclinic marine numerical model is constructed within a predetermined spatial range of the target engineering sea area; the target engineering sea area includes multiple plate exchangers; The computational domain is determined based on the boundary conditions of the three-dimensional baroclinic ocean numerical model; the boundary conditions include surface boundary conditions, bottom boundary conditions, and solid boundary conditions. The three-dimensional baroclinic ocean numerical model is configured based on the computational domain; The configured three-dimensional baroclinic ocean numerical model is used to simulate the water structure within the computational domain, and the observation data of the target engineering sea area within the target time period is obtained; the observation data includes the tide level, current velocity and current direction of each observation point. The cause of the plate heat exchanger blockage was determined based on the observed data. Specifically, configuring the three-dimensional baroclinic ocean numerical model based on the computational domain includes: Based on the computational domain, initial conditions and open boundary conditions are input into the three-dimensional baroclinic ocean numerical model; Input atmospheric boundary conditions based on the initial conditions; The bottom friction and horizontal eddy viscosity coefficients of the three-dimensional baroclinic ocean numerical model are determined based on the atmospheric boundary conditions. The time steps of the inner and outer models of the three-dimensional baroclinic ocean numerical model are determined based on the bottom friction and the horizontal eddy viscosity coefficient.

2. The method for predicting the causes of blockage in a plate heat exchanger according to claim 1, characterized in that, The initial condition is a temperature-salinity field; the temperature-salinity field is derived from HYCOM reanalysis data; The open boundary conditions include applying tidal level and temperature-salinity driving data; the temperature-salinity driving data is derived from HYCOM reanalysis data, and the frequency of applying tidal level and temperature-salinity driving data is once every 3 hours. The atmospheric boundary conditions are wind field data at a height of 10 meters above the sea surface per hour, which are derived from the ECMWF 5th generation reanalysis dataset.

3. The method for predicting the causes of blockage in a plate heat exchanger according to claim 1, characterized in that, The surface boundary conditions include wind pressure components and net heat flux; the bottom boundary conditions include bottom pressure components and bottom groundwater flux; the solid boundary conditions include the boundary normal velocity components, precipitation rate, and evaporation rate.

4. The method for predicting the causes of blockage in a plate heat exchanger according to claim 1, characterized in that, The governing equations of the three-dimensional baroclinic ocean numerical model include the continuity equation, momentum equation, density equation, temperature equation, and salinity equation.

5. The method for predicting the causes of blockage in a plate heat exchanger according to claim 1, characterized in that, The causes of plate heat exchanger blockage were determined based on the observed data, specifically including: When the flow velocity at the front end of the plate exchanger's intake is higher than the set value at the rear end of the intake, it is determined that the plate exchanger is blocked because sediment is easily drawn into the intake culvert.

6. A device for predicting the cause of blockage in a plate heat exchanger, characterized in that, The plate heat exchanger blockage cause prediction device applies the plate heat exchanger blockage cause prediction method according to any one of claims 1-5, and the plate heat exchanger blockage cause prediction device comprises: A three-dimensional baroclinic ocean numerical model construction module is used to construct a three-dimensional baroclinic ocean numerical model within a preset spatial range where the target engineering sea area is located; the target engineering sea area includes multiple plate exchangers; The computational domain determination module is used to determine the computational domain based on the boundary conditions of the three-dimensional baroclinic ocean numerical model; the boundary conditions include surface boundary conditions, bottom boundary conditions, and solid boundary conditions. A configuration module is used to configure the three-dimensional baroclinic ocean numerical model based on the computational domain; The observation data simulation module is used to run the configured three-dimensional baroclinic ocean numerical model to obtain observation data of the target engineering sea area within the target time period. The cause determination module is used to determine the cause of the plate exchanger blockage based on the observation data.

7. A computer device, comprising: A memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that the processor executes the computer program to implement the plate switch blockage cause prediction method according to any one of claims 1-5.

8. A computer-readable storage medium having a computer program stored thereon, characterized in that, When executed by a processor, the computer program implements the plate switch blockage cause prediction method as described in any one of claims 1-5.

9. A computer program product, comprising a computer program, characterized in that, When executed by a processor, the computer program implements the plate switch blockage cause prediction method as described in any one of claims 1-5.

Images