Method and device for calculating received irradiance of offshore floating photovoltaic platform

Abstract

The invention relates to the technical field of offshore photovoltaic power generation, in particular to a calculation method and device for receiving irradiance of an offshore floating photovoltaic platform. According to the method, the geometric parameters and the material attributes of the photovoltaic platform are obtained, the geometric model of the photovoltaic platform is determined according to the geometric parameters and the material attributes, the potential flow-boundary element geometric model of the photovoltaic platform is determined according to the geometric model of the photovoltaic platform, and the potential flow-boundary element geometric model of the photovoltaic platform is used for analogue simulation. The method comprises the steps of obtaining a wave field, obtaining the speed potential of the wave field, determining incident wave force, reflected wave force and radiated wave force borne by a photovoltaic platform according to the speed potential of the wave field, and substituting the incident wave force, the reflected wave force and the radiated wave force into a time-domain inclination angle dynamic change equation set to obtain a time-domain inclination angle dynamic change equation set. According to the method, the inclination angle between the photovoltaic panel of the photovoltaic platform and the horizontal plane along with the time change is determined, then the platform receiving irradiance is determined according to the inclination angle, and through the configuration mode, the offshore photovoltaic platform receiving irradiance can be accurately calculated.

Description

Technical Field

[0001] The present invention relates to the technical field of offshore photovoltaic power generation, and particularly relates to a calculation method and device for the received irradiance of an offshore floating photovoltaic platform. Background Art

[0002] With the continuous growth of the global demand for clean energy, offshore photovoltaic platforms, as an important part of the development of clean energy, are gradually becoming a hot topic in research and application. In China, a certain number of offshore photovoltaic platforms have successfully achieved project implementation, showing good development prospects. However, as an emerging offshore structure, offshore photovoltaic platforms have unique structural characteristics and working environments. Their windward area is relatively large, and compared with traditional offshore engineering structures, the influence of wind loads is more significant, which needs to be considered particularly in the design and construction process.

[0003] In addition, due to the offshore moving state of the offshore photovoltaic platform, the existing technology cannot accurately calculate the received irradiance of it. And accurately obtaining irradiance data is crucial for evaluating the power generation efficiency of the offshore photovoltaic platform, optimizing the layout of photovoltaic modules, and improving energy utilization efficiency.

[0004] Based on this, the present invention proposes a calculation method and device for the received irradiance of an offshore floating photovoltaic platform to solve the problem of how to accurately calculate the received irradiance of an offshore photovoltaic platform. Summary of the Invention

[0005] To solve the problem of how to accurately calculate the received irradiance of an offshore photovoltaic platform, an embodiment of the present invention provides a calculation method and device for the received irradiance of an offshore floating photovoltaic platform.

[0006] In a first aspect, an embodiment of the present invention provides a calculation method for the received irradiance of an offshore floating photovoltaic platform, the method comprising:

[0007] Obtain the geometric parameters and material properties of the photovoltaic platform;

[0008] Based on the geometric parameters and the material properties, determine the geometric model of the photovoltaic platform;

[0009] Based on the geometric model of the photovoltaic platform, determine the potential flow-boundary element geometric model of the photovoltaic platform;

[0010] Use the potential flow-boundary element geometric model of the photovoltaic platform to perform simulation to obtain the velocity potential of the wave field;

[0011] Based on the velocity potential of the wave field, determine the incident wave force, reflected wave force, and radiation wave force received by the photovoltaic platform;

[0012] Substitute the incident wave force, the reflected wave force, and the radiation wave force into the time-domain dynamic change equation set of the tilt angle to determine the tilt angle of the photovoltaic panel of the photovoltaic platform with respect to the horizontal plane over time; wherein, the time-domain dynamic change equation set of the tilt angle is constructed by using the lift coefficient and the drag coefficient of the photovoltaic platform under the action of wind load, and the lift coefficient and the drag coefficient change dynamically according to the change of the pitch angle of the photovoltaic platform;

[0013] Based on the tilt angle, determine the irradiance received by the platform.

[0014] In a first aspect, an embodiment of the present invention provides a calculation device for the irradiance received by an offshore floating photovoltaic platform, including:

[0015] An acquisition module, configured to acquire the geometric parameters and material properties of the photovoltaic platform;

[0016] A first data processing module, configured to determine the geometric model of the photovoltaic platform based on the geometric parameters and the material properties;

[0017] A second data processing module, configured to determine the potential flow-boundary element geometric model of the photovoltaic platform based on the geometric model of the photovoltaic platform;

[0018] A third data processing module, configured to perform simulation using the potential flow-boundary element geometric model of the photovoltaic platform to obtain the velocity potential of the wave field;

[0019] A fourth data processing module, configured to determine the incident wave force, the reflected wave force, and the radiation wave force received by the photovoltaic platform based on the velocity potential of the wave field;

[0020] A fifth data processing module, configured to substitute the incident wave force, the reflected wave force, and the radiation wave force into the time-domain dynamic change equation set of the tilt angle to determine the tilt angle of the photovoltaic panel of the photovoltaic platform with respect to the horizontal plane over time; wherein, the time-domain dynamic change equation set of the tilt angle is constructed by using the lift coefficient and the drag coefficient of the photovoltaic platform under the action of wind load, and the lift coefficient and the drag coefficient change dynamically according to the change of the pitch angle of the photovoltaic platform;

[0021] A sixth data processing unit, configured to determine the irradiance received by the platform based on the tilt angle.

[0022] In a third aspect, an embodiment of the present invention further provides an electronic device, including a memory and a processor. A computer program is stored in the memory, and when the processor executes the computer program, the method described in any embodiment of the present invention is implemented.

[0023] Fourthly, an embodiment of the present invention further provides a computer-readable storage medium, on which a computer program is stored. When the computer program is executed on a computer, the computer is made to execute the method described in any embodiment of the present invention.

[0024] An embodiment of the present invention provides a method and device for calculating the irradiance received by an offshore floating photovoltaic platform. In the method for calculating the irradiance received by an offshore floating photovoltaic platform, the primary step is to collect the geometric parameters of the photovoltaic platform, including the size specifications, shape characteristics of the platform, and material properties, such as data on the density, elastic modulus, etc. of the material. Based on these data, a geometric model of the photovoltaic platform that conforms to the actual situation is constructed. This model visually presents the physical form of the photovoltaic platform. Subsequently, relying on the geometric model, through specific mathematical transformations and boundary condition settings, a potential flow-boundary element geometric model of the photovoltaic platform is determined. This model can effectively consider the interaction between the fluid and the structure, providing strong support for simulating the impact of waves on the photovoltaic platform. With the help of simulation software, a simulation is carried out using the potential flow-boundary element geometric model of the photovoltaic platform to solve the velocity potential of the wave field. This velocity potential reflects the energy distribution and flow state of the waves in space and time. Based on the obtained velocity potential of the wave field, the incident wave force, reflected wave force, and radiation wave force received by the photovoltaic platform are determined. These three wave forces describe the mechanical effects of waves on the photovoltaic platform from different angles and are important parameters for analyzing the motion state of the platform. Substitute the incident wave force, reflected wave force, and radiation wave force into the time-domain dynamic change equation set of the tilt angle. The construction of this equation set is based on the lift coefficient and drag coefficient of the photovoltaic platform under the action of wind load, and the lift coefficient and drag coefficient will be dynamically adjusted according to the real-time change of the pitch angle of the photovoltaic platform. This dynamic change characteristic enables the time-domain dynamic change equation set of the tilt angle to more truly and accurately reflect the motion state of the photovoltaic platform in a complex marine environment, and then more accurately solve the tilt angle of the photovoltaic panel of the photovoltaic platform with respect to the horizontal plane over time. By accurately grasping the change of the tilt angle, the accuracy of determining the irradiance received by the platform is significantly improved, providing more reliable data support for the evaluation of the power generation efficiency of the offshore photovoltaic platform and the optimization of the layout of photovoltaic modules. Description of the Drawings

[0025] In order to more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the following will briefly introduce the drawings required for use in the description of the embodiments or the prior art. Obviously, the drawings in the following description are some embodiments of the present invention. For those of ordinary skill in the art, without creative efforts, other drawings can also be obtained based on these drawings.

[0026] Figure 1 Shows a flowchart of a method for calculating the irradiance received by an offshore floating photovoltaic platform according to an embodiment;

[0027] Figure 2 It is a hardware architecture diagram of an electronic device provided by an embodiment of the present invention;

[0028] Figure 3 It shows a structural diagram of a calculation device for receiving irradiance of an offshore floating photovoltaic platform according to an embodiment. Specific embodiments

[0029] To make the objectives, technical solutions, and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present invention. Apparently, the described embodiments are some, but not all, of the embodiments of the present invention. All other embodiments obtained by those of ordinary skill in the art based on the embodiments of the present invention without creative efforts shall fall within the protection scope of the present invention.

[0030] Please refer to Figure 1 , an embodiment of the present invention provides a method for calculating the received irradiance of an offshore floating photovoltaic platform, and the method includes:

[0031] Step 100: Obtain the geometric parameters and material properties of the photovoltaic platform;

[0032] Step 102: Determine the geometric model of the photovoltaic platform based on the geometric parameters and material properties;

[0033] Step 104: Determine the potential flow-boundary element geometric model of the photovoltaic platform based on the geometric model of the photovoltaic platform;

[0034] Step 106: Use the potential flow-boundary element geometric model of the photovoltaic platform for simulation to obtain the velocity potential of the wave field;

[0035] Step 108: Determine the incident wave force, reflected wave force, and radiation wave force received by the photovoltaic platform based on the velocity potential of the wave field;

[0036] Step 110: Substitute the incident wave force, reflected wave force, and radiation wave force into the time-domain tilt angle dynamic change equation set to determine the tilt angle of the photovoltaic panel of the photovoltaic platform with respect to the horizontal plane over time; wherein, the time-domain tilt angle dynamic change equation set is constructed based on the lift coefficient and drag coefficient of the photovoltaic platform under wind load, and the lift coefficient and drag coefficient change dynamically according to the change of the pitch angle of the photovoltaic platform;

[0037] Step 112: Determine the irradiance received by the platform based on the tilt angle.

[0038] In this embodiment, in the method for calculating the irradiance received by an offshore floating PV platform, the primary step is to collect the geometric parameters of the PV platform, including the size specifications, shape characteristics of the platform, and material properties such as the density and elastic modulus of the material. Based on these data, a geometric model of the PV platform that conforms to the actual situation is constructed. This model visually presents the physical form of the PV platform. Subsequently, relying on the geometric model, through specific mathematical transformations and boundary condition settings, a potential flow-boundary element geometric model of the PV platform is determined. This model can effectively consider the interaction between the fluid and the structure, providing strong support for simulating the impact of waves on the PV platform. With the help of simulation software, a simulation is carried out using the potential flow-boundary element geometric model of the PV platform to solve the velocity potential of the wave field. This velocity potential reflects the energy distribution and flow state of the waves in space and time. Based on the obtained velocity potential of the wave field, the incident wave force, reflected wave force, and radiation wave force acting on the PV platform are determined. These three wave forces describe the mechanical effects of waves on the PV platform from different angles and are important parameters for analyzing the motion state of the platform. Substitute the incident wave force, reflected wave force, and radiation wave force into the time-domain dynamic change equation set of the tilt angle. The construction of this equation set is based on the lift coefficient and drag coefficient of the PV platform under the action of wind load, and the lift coefficient and drag coefficient will be dynamically adjusted according to the real-time change of the pitch angle of the PV platform. This dynamic change characteristic enables the time-domain dynamic change equation set of the tilt angle to more realistically and accurately reflect the motion state of the PV platform in a complex marine environment, and further more accurately solve the tilt angle of the PV panels of the PV platform with respect to the horizontal plane over time. By precisely grasping the change of the tilt angle, the accuracy of determining the irradiance received by the platform is significantly improved, providing more reliable data support for the evaluation of the power generation efficiency of the offshore PV platform and the optimization of the layout of PV modules.

[0039] In one embodiment of the present invention, the velocity potential of the wave field is determined by the following formula:

[0040]

[0041] a w =h / 2

[0042] ω=2π / T

[0043] k=2π / L

[0044] In the formula, a w is the wave amplitude, i is the imaginary number, ω is the wave circular frequency, is the velocity potential of the wave field, is the position vector of the PV platform, g is the acceleration due to gravity, t is the time, Z is the coordinate of the PV platform along the water depth direction, d is the water depth, k is the wave number, R is the radial radius of the PV platform, H is the wave height, T is the wave period, and L is the wave length.

[0045] In this embodiment, those skilled in the art can customize the amplitude, complex number, circular frequency of the wave, position vector of the photovoltaic platform, gravitational acceleration, time, coordinate of the photovoltaic platform along the water depth direction, water depth, wave number, radial radius of the photovoltaic platform, wave height, wave period, and wavelength according to the actual usage situation. Through the above equations, the velocity potential of the wave field of the photovoltaic platform at different times can be accurately determined.

[0046] In an embodiment of the present invention, the incident wave force, reflected wave force, and radiation wave force are determined by the following formulas:

[0047]

[0048] In the formula, F I is the sum of the incident wave force and the hydrostatic pressure, F d is the reflected wave force, F r is the radiation wave force, is the incident velocity potential, is the reflected velocity potential, is the radiation velocity potential, ρ is the density of water, z 0 is the depth of the photovoltaic platform immersed in water, and S is the surface area of the photovoltaic platform.

[0049] In this embodiment, using the boundary element method, the control equations of each unit after grid division are discretely solved. The incident wave force, reflected wave force, and radiation wave force acting on the structure under the action of waves can be obtained.

[0050] In an embodiment of the present invention, the time-domain tilt angle dynamic change equation set is constructed by the following formulas:

[0051]

[0052] F(t) = F I (t) + F d (t) + F e (t) + F t (t)

[0053]

[0054] In the formula, β is the tilt angle, X is the displacement of the photovoltaic platform, M s is the mass matrix of the photovoltaic platform, C is the hydrodynamic damping matrix, A ∞ is the added mass matrix of the fluid, K hys is the hydrostatic stiffness matrix, h(t - τ) is the convolution integral of the wave radiation force, F(t) is the total external force acting on the photovoltaic platform, and F t (t) is the sum of the connection force between the mooring force and the photovoltaic platform, Fe (t) is the wind load force, ρ is the air density, A is the area subjected to wind force, V wind is the wind velocity, C is the wind coefficient which is determined by the lift coefficient and the drag coefficient, θ 0 is the initial tilt angle of the photovoltaic panel to the horizontal plane, is the rotational component of the photovoltaic platform displacement in the Z-axis direction, V s is the component of the platform velocity along the X-axis, ω p is the rotational component of the photovoltaic platform velocity along the Y axis, H+h COG It is the height of the center of gravity of the photovoltaic panel from the overall center of gravity of the platform.

[0055] In one embodiment of the present invention, the irradiance received by the platform is determined by the following formula:

[0056] G η =B η +D η +R η

[0057] B η =DNI×r b

[0058] r b =max(0,cosθ / cosθ z )

[0059] cosθ=cosθ z cosβ+sinθ z sinβcos(γ s -γ)

[0060]

[0061] a=max(0,cosθ)

[0062] b=max(cos85°,cosθ z )

[0063]

[0064] In the formula, G η is the irradiance received by the platform, B η For direct irradiation, D η is diffuse irradiance, R η is reflected irradiance, DNI is direct irradiance, DHI is diffuse irradiance, GHI is radiant irradiance, AM is air mass, E a is the extraterrestrial radiation, η is the azimuth angle of the photovoltaic panel, and γ is the direction of the photovoltaic panel toward θ z The solar zenith angle, γ sis the solar azimuth angle, f 11 , f 12 , f 13 , f 21 , f 22 and f 23 are the first correlation coefficient, the second correlation coefficient, the third correlation coefficient, the fourth correlation coefficient, the fifth correlation coefficient and the sixth correlation coefficient in sequence.

[0065] In this embodiment, those skilled in the art can customize the first correlation coefficient, the second correlation coefficient, the third correlation coefficient, the fourth correlation coefficient, the fifth correlation coefficient and the sixth correlation coefficient according to the actual usage situation, and accumulate the irradiance of each time step within a day to obtain the total sunshine amount received by the floating photovoltaic platform at sea throughout the day.

[0066] In an embodiment of the present invention, the lift coefficient and the drag coefficient are determined through the following steps:

[0067] Construct a geometric model of the meshed photovoltaic panel array and the wind field calculation domain;

[0068] Based on the geometric model of the meshed photovoltaic panel array and the wind field calculation domain and the photovoltaic panel-wind field coupling control equation, simulate the action of the wind field on the photovoltaic panel array to determine the velocity distribution and pressure distribution of the wind field;

[0069] Substitute the velocity distribution and pressure distribution into the preset Morison equation to determine the lift coefficient and the drag coefficient.

[0070] In this embodiment, relying on the geometric model of the meshed photovoltaic panel array and the wind field calculation domain, combined with the photovoltaic panel-wind field coupling control equation, a numerical simulation of the interaction between the photovoltaic panel array and the wind field is carried out. During the simulation process, complex factors such as the flow characteristics of the wind, the blocking and interference effects of the photovoltaic panels are fully considered. By solving the coupling control equation, the velocity distribution and pressure distribution of the wind field around the photovoltaic panel array are accurately determined. These distribution data intuitively reflect the energy transfer and mechanical action mechanism between the wind field and the photovoltaic panels. Substitute the simulated wind field velocity distribution and pressure distribution data into the preset Morison equation. The Morison equation is a classic formula widely used in the field of fluid mechanics, which can effectively describe the relationship between the force exerted on an object in a fluid and the fluid characteristics. Through the calculation of this equation, the lift coefficient and the drag coefficient of the photovoltaic panel under the action of the wind field are accurately determined.

[0071] In an embodiment of the present invention, the photovoltaic panel-wind field coupling control equation is as follows:

[0072]

[0073] In the formula, is the velocity vector, x i , x j is the spatial tensor component, t is the time, is the pressure, ρ is the fluid density, v t is the kinematic viscosity, is the product of the velocity pulsation components, σ k is a constant related to the turbulence model, C 1， C 2 are all empirical constants, ε is the turbulent kinetic energy dissipation rate, K is the turbulent kinetic energy, P ij is the tensor component related to the pressure, δ ij is the Kronecker symbol.

[0074] As Figure 2 , Figure 3 shown, the embodiment of the present invention provides a calculation device for the received irradiance of an offshore floating photovoltaic platform. The device embodiment can be implemented by software, or by hardware or a combination of software and hardware. In terms of the hardware level, as Figure 2 shown, it is a hardware architecture diagram of an electronic device where the calculation device for the received irradiance of an offshore floating photovoltaic platform provided by the embodiment of the present invention is located. In addition to Figure 2 shown processor, memory, network interface, and non-volatile memory, the electronic device where the device is located in the embodiment usually may also include other hardware, such as a forwarding chip responsible for processing packets, and so on. Taking software implementation as an example, as Figure 3 shown, as a device in a logical sense, it is formed by the CPU of its corresponding electronic device reading the computer program in the non-volatile memory into the memory and running it.

[0075] As Figure 3 shown, the calculation device for the received irradiance of an offshore floating photovoltaic platform provided in this embodiment, the device includes:

[0076] An acquisition module 300, configured to acquire the geometric parameters and material properties of the photovoltaic platform;

[0077] A first data processing module 302, configured to determine the geometric model of the photovoltaic platform based on the geometric parameters and the material properties;

[0078] A second data processing module 304, configured to determine the potential flow-boundary element geometric model of the photovoltaic platform based on the geometric model of the photovoltaic platform;

[0079] A third data processing module 306, configured to perform simulation using the potential flow-boundary element geometric model of the photovoltaic platform to obtain the velocity potential of the wave field;

[0080] The fourth data processing module 308 is configured to determine the incident wave force, the reflected wave force, and the radiation wave force received by the photovoltaic platform based on the velocity potential of the wave field;

[0081] The fifth data processing module 310 is configured to substitute the incident wave force, the reflected wave force, and the radiation wave force into the time-domain dynamic change equation set of the tilt angle to determine the tilt angle of the photovoltaic panel of the photovoltaic platform with respect to the horizontal plane over time; wherein, the time-domain dynamic change equation set of the tilt angle is constructed by the lift coefficient and the drag coefficient of the photovoltaic platform under the action of wind load, and the lift coefficient and the drag coefficient change dynamically according to the change of the pitch angle of the photovoltaic platform;

[0082] The sixth data processing unit 312 is configured to determine the irradiance received by the platform based on the tilt angle.

[0083] In an embodiment of the present invention, the velocity potential of the wave field is determined by the following formula:

[0084]

[0085] a w = H / 2

[0086] ω = 2π / T

[0087] k = 2π / L

[0088] In the formula, a w is the wave amplitude, i is the complex number, ω is the wave circular frequency, is the velocity potential of the wave field, is the position vector of the photovoltaic platform, g is the acceleration due to gravity, t is the time, Z is the coordinate of the photovoltaic platform along the water depth direction, d is the water depth, k is the wave number, R is the radial radius of the photovoltaic platform, H is the wave height, T is the wave period, and L is the wave length.

[0089] In an embodiment of the present invention, the incident wave force, the reflected wave force, and the radiation wave force are determined by the following formula:

[0090]

[0091] In the formula, F I is the sum of the incident wave force and the hydrostatic pressure, F d is the reflected wave force, F r is the radiation wave force, is the incident velocity potential, is the reflected velocity potential, is the radiation velocity potential, ρ is the density of water, z 0 is the depth of the photovoltaic platform immersed in water, and S is the surface area of the photovoltaic platform.

[0092] In one embodiment of the present invention, the time domain tilt angle dynamic change equation group is constructed by the following formula:

[0093]

[0094] F(t)=F I (t)+F d (t)+F e (t)+F t (t)

[0095]

[0096] Where β is the tilt angle, X is the displacement of the photovoltaic platform, and M s is the mass matrix of the photovoltaic platform, C is the fluid dynamic damping matrix, A ∞ is the fluid additional mass matrix, K hys is the hydrostatic stiffness matrix, h(tt) is the convolution integral of the wave radiation force, F(t) is the sum of the external forces acting on the photovoltaic platform, and F t (t) is the sum of the mooring force and the connection force between the photovoltaic platform, F e (t) is the wind load force, ρ is the air density, A is the area subjected to wind force, V wind is the wind velocity, C is the wind coefficient which is determined by the lift coefficient and the drag coefficient, θ 0 is the initial tilt angle of the photovoltaic panel to the horizontal plane, is the rotational component of the photovoltaic platform displacement in the Z-axis direction, V s is the component of the platform velocity along the X-axis, ω p is the rotational component of the photovoltaic platform velocity along the Y axis, H+h COG It is the height of the center of gravity of the photovoltaic panel from the overall center of gravity of the platform.

[0097] In one embodiment of the present invention, the platform received irradiance is determined by the following formula:

[0098] G η =B η +D η +R η

[0099] B η =DNI×r b

[0100] r b =max(0,cosθ / cosθ z )

[0101] cosθ=cosθ z cosβ+sinθz sinβcos(γ s -γ)

[0102]

[0103] a = max(0, cosθ)

[0104] b = max(cos85°, cosθ z )

[0105]

[0106] Wherein, G η is the irradiance received by the platform, B η is the direct irradiation, D η is the diffuse irradiation, R η is the reflected irradiation, DNI is the direct normal irradiance, DHI is the diffuse horizontal irradiance, GHI is the global horizontal irradiance, AM is the air mass, E a is the extraterrestrial radiation, η is the azimuth angle of the photovoltaic panel, γ is the orientation of the photovoltaic panel towards the solar zenith angle θ z of, γ s is the solar azimuth angle, f 11 , f 12 , f 13 , f 21 , f 22 and f 23 are the first correlation coefficient, the second correlation coefficient, the third correlation coefficient, the fourth correlation coefficient, the fifth correlation coefficient and the sixth correlation coefficient in sequence.

[0107] In an embodiment of the present invention, the lift coefficient and the drag coefficient are determined by the following steps:

[0108] Construct a geometric model of the segmented photovoltaic panel array and the wind field calculation domain;

[0109] Based on the geometric model of the segmented photovoltaic panel array and the wind field calculation domain and the photovoltaic panel - wind field coupling control equation, simulate the action of the wind field on the photovoltaic panel array to determine the velocity distribution and pressure distribution of the wind field;

[0110] Substitute the velocity distribution and the pressure distribution into the preset Morrison equation to determine the lift coefficient and the drag coefficient.

[0111] In an embodiment of the present invention, the photovoltaic panel - wind field coupling control equation is as follows:

[0112]

[0113] Wherein, is the velocity vector, x i, x j is a spatial tensor component, t is time, is pressure, ρ is fluid density, v t is kinematic viscosity, is the product of velocity fluctuation components, σ k is a constant related to the turbulence model, C 1， C 2 are all empirical constants, ε is the dissipation rate of turbulent kinetic energy, K is the turbulent kinetic energy, P ij is a tensor component related to pressure, δ ij is the Kronecker symbol.

[0114] It can be understood that the structure illustrated in the embodiments of the present invention does not constitute a specific limitation on a calculation and measurement device for the received irradiance of an offshore floating photovoltaic platform. In other embodiments of the present invention, a calculation and measurement device for the received irradiance of an offshore floating photovoltaic platform may include more or fewer components than those illustrated, or combine certain components, or split certain components, or have different component arrangements. The illustrated components can be implemented in hardware, software, or a combination of software and hardware.

[0115] Regarding the information interaction, execution process, etc. between the various modules in the above device, since they are based on the same concept as the method embodiments of the present invention, the specific content can be referred to the description in the method embodiments of the present invention and will not be elaborated here.

[0116] The embodiments of the present invention also provide an electronic device, including a memory and a processor. A computer program is stored in the memory. When the processor executes the computer program, it implements a calculation method for the received irradiance of an offshore floating photovoltaic platform in any embodiment of the present invention.

[0117] The embodiments of the present invention also provide a computer-readable storage medium. A computer program is stored on the computer-readable storage medium. When the computer program is executed by a processor, it causes the processor to execute a calculation method for the received irradiance of an offshore floating photovoltaic platform in any embodiment of the present invention.

[0118] Specifically, a system or device equipped with a storage medium can be provided. Software program code for implementing the functions in any of the above embodiments is stored on the storage medium, and the computer (or CPU or MPU) of the system or device reads and executes the program code stored on the storage medium.

[0119] In this case, the program code read from the storage medium itself can implement the functions of any one of the above embodiments. Therefore, the program code and the storage medium storing the program code constitute a part of the present invention.

[0120] Examples of storage media for providing program code include floppy disks, hard disks, magneto-optical disks, optical disks (such as CD-ROM, CD-R, CD-RW, DVD-ROM, DVD-RAM, DVD-RW, DVD+RW), magnetic tapes, non-volatile memory cards, and ROMs. Optionally, the program code can be downloaded from a server computer via a communication network.

[0121] In addition, it should be clear that not only can the functions of any of the above embodiments be realized by executing the program code read by a computer, but also by causing an operating system or the like operating on the computer based on the instructions of the program code to complete part or all of the actual operations.

[0122] Furthermore, it can be understood that the program code read from the storage medium is written into the memory provided in an expansion board inserted into the computer or into the memory provided in an expansion module connected to the computer, and then based on the instructions of the program code, a CPU or the like installed on the expansion board or the expansion module is caused to execute part or all of the actual operations, thereby realizing the functions of any of the above embodiments.

[0123] It should be noted that in this article, relational terms such as first and second are only used to distinguish one entity or operation from another entity or operation, and do not necessarily require or imply any actual relationship or order between these entities or operations. Moreover, the term "comprising", "including" or any other variant thereof is intended to cover non-exclusive inclusion, so that a process, method, article or device comprising a series of elements includes not only those elements but also other elements not expressly listed, or elements inherent to such process, method, article or device.

[0124] Those of ordinary skill in the art can understand that all or part of the steps of implementing the above method embodiments can be completed by hardware related to program instructions. The foregoing program can be stored in a computer-readable storage medium. When the program is executed, it performs the steps including the above method embodiments; and the foregoing storage medium includes various media such as ROM, RAM, magnetic disks, or optical disks that can store program code.

[0125] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention and are not intended to limit them; although the present invention has been described in detail with reference to the foregoing embodiments, those of ordinary skill in the art should understand that they can still modify the technical solutions described in the foregoing embodiments, or perform equivalent replacements for some of the technical features; and these modifications or replacements do not make the essence of the corresponding technical solutions deviate from the spirit and scope of the technical solutions of the embodiments of the present invention.

Claims

1. A method for calculating the irradiance received by an offshore floating photovoltaic platform, characterized in that: The method comprises: Obtain the geometric parameters and material properties of the photovoltaic platform; Determining a geometric model of a photovoltaic platform based on the geometric parameters and the material properties; Based on the geometric model of the photovoltaic platform, determining a potential flow-boundary element geometric model of the photovoltaic platform; The photovoltaic platform potential flow-boundary element geometric model is used to perform simulation to obtain the velocity potential of the wave field; Based on the velocity potential of the wave field, determining the incident wave force, reflected wave force and radiated wave force on the photovoltaic platform; Substituting the incident wave force, the reflected wave force and the radiation wave force into the time domain tilt angle dynamic change equation group, determining the tilt angle of the photovoltaic panel of the photovoltaic platform with the horizontal plane over time; wherein the time domain tilt angle dynamic change equation group is constructed by the lift coefficient and the drag coefficient of the photovoltaic platform under the action of wind load, and the lift coefficient and the drag coefficient change dynamically according to the change of the pitch angle of the photovoltaic platform; Based on the tilt angle, the irradiance received by the platform is determined.

2. The method according to claim 1, characterized in that The velocity potential of the wave field is determined by the following formula: a w =H / 2 ω=2π / T k=2π / L In the formula, a w is the amplitude, i is a complex number, ω is the wave circular frequency, is the velocity potential of the wave field, is the position vector of the photovoltaic platform, g is the gravitational acceleration, t is the time, Z is the coordinate of the photovoltaic platform along the water depth direction, d is the water depth, k is the wave number, R is the radial radius of the photovoltaic platform, H is the wave height, T is the wave period, and L is the wavelength.

3. The method according to claim 2, characterized in that The incident wave force, the reflected wave force and the radiated wave force are determined by the following formula: In the formula, F I is the sum of the incident wave force and the hydrostatic pressure, F d is the reflected wave force, F r is the radiation wave force, is the incident velocity potential, is the reflection velocity potential, is the radiation velocity potential, ρ is the density of water, z0 is the depth of the photovoltaic platform immersed in water, and S is the surface area of ​​the photovoltaic platform.

4. The method according to claim 3, characterized in that: The time domain tilt angle dynamic change equation group is constructed by the following formula: F(t)=F I (t)+F d (t)+F e (t)+F t (t) Where β is the tilt angle, X is the displacement of the photovoltaic platform, and M s is the mass matrix of the photovoltaic platform, C is the fluid dynamic damping matrix, A ∞ is the fluid additional mass matrix, K hys is the hydrostatic stiffness matrix, h(t-τ) is the convolution integral of the wave radiation force, F(t) is the sum of the external forces acting on the photovoltaic platform, and F t (t) is the sum of the mooring force and the connection force between the photovoltaic platform, F e (t) is the wind load force, ρ is the air density, A is the area subjected to wind force, V wind is the wind velocity, C is the wind coefficient which is determined by the lift coefficient and the drag coefficient, θ0 is the initial tilt angle between the photovoltaic panel and the horizontal plane, is the rotational component of the photovoltaic platform displacement in the Z-axis direction, V s is the component of the platform velocity along the X-axis, ω p is the rotational component of the photovoltaic platform velocity along the Y axis, H+h COG It is the height of the center of gravity of the photovoltaic panel from the overall center of gravity of the platform.

5. The method according to claim 4, characterized in that The platform received irradiance is determined by the following formula: G η =B η +D η +R η B η =DNI×r b r b =max(0,cosθ / cosθ z ) cosθ=cosθ z cosβ+sinθ z sinβcos(γ s -c) a=max(0,cosθ) b=max(cos85°,cosθ z ) In the formula, G η is the irradiance received by the platform, B η For direct irradiation, D η is diffuse irradiance, R η is reflected irradiance, DNI is direct irradiance, DHI is diffuse irradiance, GHI is radiant irradiance, AM is air mass, E a is the extraterrestrial radiation, η is the azimuth angle of the photovoltaic panel, and γ is the direction of the photovoltaic panel toward θ z The solar zenith angle, γ s is the solar azimuth, f 11 , f 12 , f 13 , f 21 , f 22 and F 23 They are the first correlation coefficient, the second correlation coefficient, the third correlation coefficient, the fourth correlation coefficient, the fifth correlation coefficient and the sixth correlation coefficient.

6. The method according to claim 1, characterized in that The lift coefficient and the drag coefficient are determined by the following steps: Construct the geometric model of the photovoltaic panel array and wind farm calculation domain after segmentation; Based on the geometric model of the photovoltaic panel array and wind field calculation domain after the segmentation and the photovoltaic panel-wind field coupling control equation, the photovoltaic panel array is simulated to be affected by the wind field, and the velocity distribution and pressure distribution of the wind field are determined; The velocity distribution and the pressure distribution are substituted into a preset Morrison equation to determine the lift coefficient and the drag coefficient.

7. The method according to claim 6, characterized in that The photovoltaic panel-wind field coupling control equation is as follows: In the formula, is the velocity vector, x i , x j is the mean spatial tensor component, t is the time, is the pressure, ρ is the fluid density, v t is the kinematic viscosity, is the product of the velocity pulsation components, σ k is a constant related to the turbulence model, C 1， C2 are empirical constants, ε is the turbulent kinetic energy dissipation rate, K is the turbulent kinetic energy, P ij is the tensor component related to pressure, δ ij is the Kronecker symbol.

8. A device for calculating the irradiance received by an offshore floating photovoltaic platform, characterized in that: include: An acquisition module, used to obtain geometric parameters and material properties of the photovoltaic platform; A first data processing module, configured to determine a geometric model of a photovoltaic platform based on the geometric parameters and the material properties; A second data processing module is used to determine a potential flow-boundary element geometric model of the photovoltaic platform based on the geometric model of the photovoltaic platform; A third data processing module is used to perform simulation using the photovoltaic platform potential flow-boundary element geometric model to obtain the velocity potential of the wave field; a fourth data processing module, for determining incident wave force, reflected wave force and radiated wave force exerted on the photovoltaic platform based on the velocity potential of the wave field; A fifth data processing module is used to substitute the incident wave force, the reflected wave force and the radiation wave force into a time domain tilt angle dynamic change equation group to determine the tilt angle of the photovoltaic panel of the photovoltaic platform with the horizontal plane over time; wherein the time domain tilt angle dynamic change equation group is constructed by the lift coefficient and the drag coefficient of the photovoltaic platform under the action of wind load, and the lift coefficient and the drag coefficient change dynamically according to the change of the pitch angle of the photovoltaic platform; The sixth data processing unit is used to determine the irradiance received by the platform based on the tilt angle.

9. An electronic device, characterized in that: The method comprises a memory and a processor, wherein the memory stores a computer program, and when the processor executes the computer program, the method according to any one of claims 1 to 7 is implemented.

10. A computer-readable storage medium, characterized in that: A computer program is stored thereon, and when the computer program is executed in a computer, the computer is caused to execute the method according to any one of claims 1 to 7.

Images