Load determination method, device and equipment for high slope wind turbine generator and medium

Abstract

The invention discloses a load determination method, device and equipment for a high-slope wind turbine generator and a medium, and relates to the technical field of wind turbine generators, and the method comprises the steps: determining the wind measurement data of a target wind power plant under a high-slope terrain, constructing a terrain physical model based on the terrain information of a high slope, and obtaining a target terrain grid based on the terrain physical model; in each sector of the target wind power plant, determining a first sector influenced by the high slope terrain and a second sector not influenced by the high slope terrain; respectively simulating the first sector and the second sector based on the target terrain grid, and determining wind regime parameters in the target wind power plant by using an obtained target simulation result; based on the wind regime parameters, a corresponding first turbulent wind field model and a corresponding second turbulent wind field model are constructed for the first sector and the second sector respectively, and the target load of the wind turbine generator is determined through the first turbulent wind field model and the second turbulent wind field model. Therefore, the calculation efficiency and accuracy of the wind turbine generator load under the complex terrain can be improved.

Description

Technical Field

[0001] The present invention relates to the technical field of wind turbine generator systems, and in particular to a method, device, equipment and medium for determining the load of a high-slope wind turbine generator system. Background Art

[0002] When constructing wind turbines on high-slope terrain, accurate load calculation is crucial for ensuring safe operation and optimizing design. However, current methods for determining target loads have numerous flaws and are unable to meet practical engineering requirements.

[0003] At present, there are many deficiencies in the load calculation of high-slope wind turbines. First, although computational fluid dynamics simulation technology is often used to obtain wind parameters, the lack of actual wind measurement data from wind farms to correct the target simulation results results in a large deviation between the simulated wind conditions and the actual wind conditions. This deviation will lead to inaccurate load calculations. Second, traditional wind resource assessment methods only calculate loads based on wind parameters at hub height, ignoring the three-dimensional non-uniform wind field characteristics of the entire wind rotor surface under high-slope terrain. The high-slope terrain is complex and will cause phenomena such as airflow acceleration and circumference. The wind field is unevenly distributed in space. Relying only on a single height data is difficult to reflect the actual wind conditions, resulting in low reliability of the load calculation results. Third, existing methods cannot accurately reflect the impact of high-slope terrain on wind turbine load calculations.

[0004] Therefore, how to improve the calculation efficiency and accuracy of wind turbine loads under complex terrain is a technical problem that needs to be solved urgently. Summary of the Invention

[0005] In view of this, the purpose of the present invention is to provide a method, device, equipment and medium for determining the load of high-slope wind turbines, which can improve the efficiency and accuracy of calculating the load of wind turbines in complex terrain. The specific solution is as follows:

[0006] In a first aspect, the present application provides a method for determining the load of a high-slope wind turbine generator system, comprising:

[0007] Obtaining the current wind speed and direction time series of a target wind farm located on a high slope terrain to determine the corresponding wind measurement data, constructing a terrain physical model based on the high slope terrain information, and performing unstructured grid division on the terrain physical model to obtain a target terrain grid;

[0008] In each sector of the target wind farm, determining a first sector affected by the high side slope terrain and a second sector not affected by the high side slope terrain based on a preset sector step size and the terrain information;

[0009] Simulating the first sector and the second sector respectively based on the target terrain grid to obtain corresponding target simulation results, and determining wind condition parameters of the wind turbines and wind towers in the target wind farm using the target simulation results and the wind measurement data;

[0010] A first turbulent wind field model and a second turbulent wind field model are constructed for the first sector and the second sector respectively based on the wind condition parameters, and a target load of the wind turbine is determined using the first turbulent wind field model and the second turbulent wind field model.

[0011] Optionally, constructing a terrain physical model based on terrain information of a high slope includes:

[0012] Acquiring a satellite image and a measured topographic map of the high slope, correcting the satellite image based on a construction coordinate system of the high slope, and spatially aligning the corrected image with the measured topographic map;

[0013] After the spatial alignment, a spatial interpolation algorithm is used to complete the elevation data of the missing areas in the measured topographic map to generate a first topographic map, and the topographic feature points and the coordinates of the wind turbine generator set are marked in the first topographic map to obtain a second topographic map;

[0014] Adjusting the second topographic map using the design parameters of the high slope construction platform and preset three-dimensional modeling software to obtain a third topographic map;

[0015] The preset horizontal sampling interval and the preset vertical sampling interval are determined as modeling parameters, and a terrain physical model is constructed based on the third topographic map and the modeling parameters through the preset three-dimensional modeling software; the sampling interval is the interval at which the preset three-dimensional modeling software performs terrain sampling based on the third topographic map.

[0016] Optionally, performing unstructured grid division on the terrain physical model to obtain a target terrain grid includes:

[0017] Performing unstructured grid division on the terrain physical model using preset computational fluid dynamics software to obtain a first terrain grid;

[0018] Densifying the terrain that meets the preset complexity condition in the first terrain grid based on the global grid size control parameter, and generating a preset number of prismatic layer grids in the surface boundary layer of the first terrain grid to obtain a second terrain grid;

[0019] In the second terrain grid, a target encryption area is established with the wind turbine generator set as the center, and the target encryption area is encrypted based on a preset grid encryption strategy to obtain a target terrain grid.

[0020] Optionally, simulating the first sector and the second sector respectively based on the target terrain grid to obtain corresponding target simulation results includes:

[0021] Determining the computational domain of the preset computational fluid dynamics software;

[0022] When simulating the first sector, the inlet boundary of the computational domain is set to a fluctuating wind speed profile including atmospheric boundary layer characteristics, the outlet boundary of the computational domain is set to convection, and a preset wall function is set in the ground boundary of the computational domain to complete the first boundary setting operation;

[0023] After completing the first boundary setting operation, simulating the first sector using the large eddy simulation method and the target terrain grid using the preset computational fluid dynamics software to obtain a three-dimensional wind velocity field and turbulence field with a preset time resolution within the swept surface of the wind turbine rotor, so as to determine a corresponding first simulation result;

[0024] When simulating the second sector, determining a current wind speed profile based on current wind measurement data, and setting an inlet boundary of the computational domain to the current wind speed profile of uniform inflow to complete a second boundary setting operation;

[0025] After completing the second boundary setting operation, simulating the second sector using the preset Reynolds time-averaged turbulence model and the target terrain grid through the preset computational fluid dynamics software to obtain average wind speed gradient data within the wind turbine rotor plane and a reference wind speed and turbulence intensity distribution at the wind turbine hub, so as to determine a corresponding second simulation result;

[0026] A target simulation result is determined based on the first simulation result and the second simulation result.

[0027] Optionally, the determining wind condition parameters of the wind turbines and wind towers in the target wind farm using the target simulation results and the wind measurement data includes:

[0028] determining a preset wind direction angle range based on a first central wind direction angle of the third sector and a second central wind direction angle of the fourth sector, and randomly selecting a wind direction angle from the preset wind direction angle range as the current inflow wind direction angle; wherein the third sector and the fourth sector are adjacent to each other in the target wind farm;

[0029] Determine, based on the target simulation result, a first simulated wind speed of the wind measuring tower at the first central wind direction angle and a second simulated wind speed of the wind measuring tower at the second central wind direction angle;

[0030] Determining a first target wind speed of the wind measurement tower using the first simulated wind speed, the second simulated wind speed, and the current inflow wind direction angle;

[0031] Determine, based on the target simulation result, a first simulated wind direction of the wind measurement tower at the first central wind direction angle and a second simulated wind direction of the wind measurement tower at the second central wind direction angle;

[0032] Determining a first target wind direction of the wind measurement tower using the first simulated wind direction, the second simulated wind direction, and the current inflow wind direction angle;

[0033] Determine, based on the target simulation result, a third simulated wind speed of the wind turbine generator set at the first central wind direction angle and a fourth simulated wind speed of the wind turbine generator set at the second central wind direction angle;

[0034] Determining a second target wind speed of the wind turbine generator system using the third simulated wind speed, the fourth simulated wind speed, and the current inflow wind direction angle;

[0035] Determining, based on the target simulation result, a third simulated wind direction of the wind turbine generator set at the first central wind direction angle and a fourth simulated wind direction of the wind turbine generator set at the second central wind direction angle;

[0036] Determining a second target wind direction of the wind turbine generator system using the third simulated wind direction, the fourth simulated wind direction, and the current inflow wind direction angle;

[0037] determining a wind acceleration factor based on the first target wind speed and the second target wind speed, and correcting the second target wind speed using the wind acceleration factor and the current wind speed to obtain a corrected wind speed;

[0038] Corresponding wind condition parameters are determined based on the first target wind speed, the first target wind direction, the corrected wind speed, and the second target wind direction.

[0039] Optionally, constructing a corresponding first turbulent wind field model and a second turbulent wind field model for the first sector and the second sector respectively based on the wind condition parameters, and determining a target load of the wind turbine generator set using the first turbulent wind field model and the second turbulent wind field model includes:

[0040] constructing a rotor surface wind speed distribution matrix and a turbulence intensity distribution matrix for the wind turbine in the first sector based on the wind condition parameters, and determining a rotor surface power spectrum density matrix using the rotor surface wind speed distribution matrix and the turbulence intensity distribution matrix;

[0041] Capturing wind frequency domain data from the wind rotor surface power spectrum density matrix using a preset time step, and converting the wind frequency domain data into wind time domain data using a Fourier transform method;

[0042] Constructing a first turbulent wind field model of the first sector based on the wind condition time domain data;

[0043] Constructing a second turbulent wind field model for the wind turbine in the second sector based on the wind condition parameters and preset dynamic simulation software;

[0044] The target load of the wind turbine generator set is determined by using the first turbulent wind field model and the second turbulent wind field model through a preset dynamic simulation platform.

[0045] Optionally, after determining the target load of the wind turbine generator set by using the first turbulent wind field model and the second turbulent wind field model, the method further includes:

[0046] Obtaining a maximum load under each wind speed condition from the target load, determining a probability distribution function of the maximum load using a Gumbel distribution method, and determining a limit load of the first sector based on the probability distribution function;

[0047] Determining a preset sector step size, determining a ratio of a wind direction range to the preset sector step size as a sector number, and dividing the wind measurement data based on the sector number to determine a wind frequency matrix;

[0048] Randomly select a load from the target load as the current load;

[0049] Obtaining the number of wind speed intervals corresponding to the sectors in the wind frequency matrix, determining the load cycle frequency of the current load in the wind speed interval and the duration of the current wind speed in the wind speed interval, and determining the maximum number of cycles of the current load based on a preset load spectrum;

[0050] A fatigue load is determined based on the number of sectors, the number of intervals, the load cycle frequency, the duration, and the maximum number of cycles.

[0051] In a second aspect, the present application provides a load determination device for a high-slope wind turbine generator system, comprising:

[0052] A gridding module is used to obtain the current wind speed and wind direction time series of the target wind farm located on the high slope terrain to determine the corresponding wind measurement data, build a terrain physical model based on the terrain information of the high slope, and perform unstructured grid division on the terrain physical model to obtain the target terrain grid;

[0053] a sector determination module, configured to determine, in each sector of the target wind farm, a first sector affected by the high side slope terrain and a second sector not affected by the high side slope terrain based on a preset sector step size and the terrain information;

[0054] a wind condition determination module, configured to simulate the first sector and the second sector respectively based on the target terrain grid to obtain corresponding target simulation results, and determine wind condition parameters of the wind turbines and wind towers in the target wind farm using the target simulation results and the wind measurement data;

[0055] A load determination module is used to construct a corresponding first turbulent wind field model and a second turbulent wind field model for the first sector and the second sector respectively based on the wind condition parameters, and use the first turbulent wind field model and the second turbulent wind field model to determine the target load of the wind turbine group.

[0056] In a third aspect, the present application provides an electronic device, comprising:

[0057] Memory, used to store computer programs;

[0058] The processor is used to execute the computer program to implement the above-mentioned load determination method for the high-slope wind turbine group.

[0059] In a fourth aspect, the present application provides a computer-readable storage medium for storing a computer program; wherein, when the computer program is executed by a processor, the aforementioned method for determining the load of a high-slope wind turbine is implemented.

[0060] In the present application, a time series of current wind speed and wind direction of a target wind farm located under a high-slope terrain is obtained to determine corresponding wind measurement data, a terrain physical model is constructed based on the terrain information of the high slope, and the terrain physical model is subjected to unstructured grid division to obtain a target terrain grid; in each sector of the target wind farm, a first sector affected by the high-slope terrain and a second sector not affected by the high-slope terrain are determined based on a preset sector step size and the terrain information; the first sector and the second sector are respectively simulated based on the target terrain grid to obtain corresponding target simulation results, and the wind condition parameters of the wind turbines and the wind measurement tower in the target wind farm are determined using the target simulation results and the wind measurement data; based on the wind condition parameters, corresponding first turbulent wind field models and second turbulent wind field models are respectively constructed for the first sector and the second sector, and the target load of the wind turbine is determined using the first turbulent wind field model and the second turbulent wind field model. As can be seen from the above, in this application, the current wind speed and wind direction time series of the wind farm are first obtained to determine the wind measurement data, and at the same time, a terrain physical model is constructed based on the high slope terrain information, and the terrain physical model is divided into an unstructured grid to obtain the target terrain grid; then, in each sector of the target wind farm, according to the preset sector step size and terrain information, the first sector affected by the high slope terrain and the second sector not affected by it are determined; then, based on the target terrain grid, the first sector and the second sector are simulated respectively to obtain the corresponding target simulation results, and then these simulation results and wind measurement data are used to determine the wind condition parameters of the wind turbine and the wind tower in the target wind farm; finally, according to the wind condition parameters, the first turbulent wind field model and the second turbulent wind field model are constructed for the first sector and the second sector respectively, and the target load of the wind turbine is determined with the help of these two models. In this way, the present application can improve the calculation efficiency and accuracy of the wind turbine load under complex terrain, which is beneficial to the safety assessment of high slope wind turbines. BRIEF DESCRIPTION OF THE DRAWINGS

[0061] In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the following briefly introduces the drawings required for use in the embodiments or the description of the prior art. Obviously, the drawings described below are merely embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on the provided drawings without paying any creative work.

[0062] Figure 1 This is a flow chart of a method for determining the load of a high-slope wind turbine generator system disclosed in this application;

[0063] Figure 2 This is a flow chart of a specific method for determining the load of a high-slope wind turbine generator system disclosed in this application;

[0064] Figure 3 This is a schematic structural diagram of a load determination device for a high-slope wind turbine generator system disclosed in this application;

[0065] Figure 4 This is a structural diagram of an electronic device disclosed in this application. DETAILED DESCRIPTION

[0066] The following will clearly and completely describe the technical solutions in the embodiments of the present invention in conjunction with the accompanying drawings. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without making creative efforts are within the scope of protection of the present invention.

[0067] Currently, there are many deficiencies in the load calculation of high-slope wind turbines. First, although computational fluid dynamics simulation technology is often used to obtain wind parameters, the lack of actual wind measurement data from wind farms to correct the target simulation results results in a large deviation between the simulated wind conditions and the actual wind conditions. This deviation will lead to inaccurate load calculations. Second, traditional wind resource assessment methods only calculate loads based on wind parameters at hub height, ignoring the three-dimensional non-uniform wind field characteristics of the entire wind wheel surface under high-slope terrain. The high-slope terrain is complex, which will cause phenomena such as airflow acceleration and vortex flow. The wind field is unevenly distributed in space. It is difficult to reflect the actual wind conditions by relying solely on a single height data, resulting in low reliability of the load calculation results. Third, existing methods cannot accurately reflect the impact of high-slope terrain on wind turbine load calculations. To this end, the present application provides a load determination method, device, equipment and medium for high-slope wind turbines, which can improve the calculation efficiency and accuracy of wind turbine loads under complex terrain.

[0068] See also Figure 1 As shown, an embodiment of the present invention discloses a method for determining the load of a high-slope wind turbine generator set, comprising:

[0069] Step S11: Obtain the current wind speed and wind direction time series of the target wind farm located under the high slope terrain to determine the corresponding wind measurement data, construct a terrain physical model based on the terrain information of the high slope, and perform unstructured grid division on the terrain physical model to obtain the target terrain grid.

[0070] In this embodiment, wind measurement data for a target wind farm located on a high-slope terrain must first be obtained. Specifically, wind towers and other equipment deployed within the target wind farm are used to collect real-time wind speed and direction data, as well as time series data such as temperature and air pressure, under the high-slope terrain. Furthermore, the wind measurement height covers the height range of the wind turbine hub. The acquisition time can be 10 minutes. After acquisition, the raw data is preprocessed, including operations such as removing outliers and filling in missing data, ultimately generating wind measurement data that can be used for subsequent analysis.

[0071] After obtaining wind measurement data, a terrain physical model must be constructed based on the terrain information of the high slope. Specifically, there are various ways to obtain high-slope terrain information, such as obtaining high-resolution satellite imagery through satellite remote sensing technology, combined with measured terrain maps obtained through field measurements. The measured terrain maps can be obtained through field measurements using measuring equipment such as total stations and GPS (Global Positioning System). The measurement range must cover the target wind farm and surrounding areas that have an impact on the wind farm. After obtaining the satellite imagery and measured terrain maps, the satellite imagery is first geometrically corrected based on the high slope's construction coordinate system to eliminate image distortion caused by factors such as the satellite shooting angle and the curvature of the earth. Once the correction is complete, the corrected satellite imagery is spatially aligned with the measured terrain map to align their coordinate systems and ensure consistency in terrain information.

[0072] After the spatial alignment is completed, the spatial interpolation algorithm is used to complete the elevation data for the missing areas that may exist in the measured topographic map. The spatial interpolation algorithm can estimate the elevation value of the missing area based on the data of the surrounding known elevation points, thereby generating a complete first topographic map. Among them, the horizontal accuracy of the first topographic map can be , the elevation accuracy can be Based on the first topographic map, we further annotate terrain features, such as mountain tops, valleys, ridgelines, and slope changes, and also mark the specific coordinates of the wind turbines within the terrain to form a second topographic map. As can be seen, this annotated information more clearly reflects the relationship between terrain features and wind turbine locations.

[0073] Next, the second topographic map was adjusted using the design parameters of the high-slope construction platform and pre-set 3D modeling software. These design parameters included, but were not limited to, platform elevation, slope ratio, and retaining wall location. By importing the second topographic map into the pre-set 3D modeling software and modifying the terrain based on the construction platform design parameters, a third topographic map was generated that better matched the actual project conditions. It should be noted that in the process of generating the third topographic map, the transition area between the slope excavation surface and the natural terrain was carefully considered to ensure that the modified model met the design specifications while retaining the actual terrain characteristics.

[0074] When constructing a terrain physical model, appropriate modeling parameters must be determined, including preset horizontal and vertical sampling intervals. These sampling interval parameters determine the density of terrain sampling performed by the preset 3D modeling software on the third topographic map. The horizontal sampling interval controls the discretization of the terrain in the horizontal direction, while the vertical sampling interval controls the sampling accuracy in the elevation direction. The horizontal sampling interval can be 2 meters, and the vertical sampling interval can be 0.5 meters. Based on the third topographic map and the determined modeling parameters, the preset 3D modeling software discretizes the continuous terrain surface into a 3D terrain physical model. This model accurately reflects the geometric characteristics of high-slope terrain.

[0075] After constructing the terrain physical model, it needs to be meshed using an unstructured grid. This is done using pre-defined computational fluid dynamics software, such as CFD (Computational Fluid Dynamics), using unstructured grid generation techniques. Unstructured grids offer excellent adaptability and can better conform to complex, high-slope terrain surfaces. After meshing the terrain physical model using an unstructured grid, the first terrain grid is generated, which serves as the initial unstructured grid.

[0076] To improve calculation accuracy, especially in areas with complex terrain and near wind turbines, the grid needs to be encrypted. Based on the global grid size control parameters, the local grid is encrypted in the terrain areas of the first terrain grid that meet the preset complex conditions to more precisely capture the terrain details in these areas. At the same time, a preset number of prismatic layer grids are generated in the surface boundary layer of the first terrain grid to accurately simulate the boundary layer flow characteristics near the surface, which is crucial for wind farm calculations. After this processing, the second terrain grid is obtained.

[0077] Finally, in the second terrain grid, a target encryption area is established with the location of the wind turbine as the center. The scope of the target encryption area can be determined according to the size of the wind turbine and the range of influence of the flow field. Based on the preset grid encryption strategy, the target encryption area is further meshed. That is, a multi-level local grid encryption strategy is implemented for the area where the wind turbine machine site is located. A cylindrical encryption domain is established with the machine site as the center, and a refined grid discretization is used within the swept range of the wind rotor, and a transition zone is set to achieve a smooth change in the grid size. In particular, boundary layer encryption is performed in key areas such as near the hub height to ensure that the three-dimensional flow characteristics of the wind rotor plane can be accurately captured, while controlling the overall grid size within a reasonable range. Through this series of meshing and encryption operations, a full target terrain grid is finally obtained, laying the foundation for subsequent wind field simulation and load calculation.

[0078] Step S12: Determine, in each sector of the target wind farm, a first sector affected by the high side slope terrain and a second sector not affected by the high side slope terrain based on a preset sector step size and the terrain information.

[0079] In this embodiment, the target wind farm needs to be divided into multiple sectors first. The sector division is performed based on the preset sector step size. Usually, the range is divided into several sectors at a certain angle interval to form a sector set covering the entire target wind farm. For example, the preset sector step size is set to , the calculation formula for the number of sectors is , that is, the preset sector step size is At this time, the number of sectors is 16.

[0080] After completing the sector division, the impact of high-slope terrain is analyzed for each sector. Specifically, the elevation difference between the highest point of the slope and the construction platform within each sector is calculated. If the calculated elevation difference is greater than a preset difference threshold, the sector is determined to be affected by the high-slope terrain, i.e., the first sector. Conversely, if the elevation difference is less than or equal to the preset difference threshold, the sector is determined to be unaffected by the high-slope terrain, i.e., the second sector. The preset difference threshold can be 5 meters.

[0081] Step S13: Simulate the first sector and the second sector respectively based on the target terrain grid to obtain corresponding target simulation results, and determine the wind condition parameters of the wind turbines and wind towers in the target wind farm using the target simulation results and the wind measurement data.

[0082] In this embodiment, after the target terrain grid is constructed and sectorized, different simulation strategies need to be used to simulate the flow field according to the characteristics of different sectors. First, the calculation domain of the preset computational fluid dynamics software is determined.

[0083] When simulating the first sector, specific boundary conditions need to be set to accurately simulate the airflow characteristics under complex terrain. The inlet boundary of the computational domain is set to a pulsating wind speed profile that includes the characteristics of the atmospheric boundary layer. This profile takes into account the temporal and spatial correlation of atmospheric turbulence and can more realistically reflect the wind speed changes in the actual atmospheric boundary layer. The outlet boundary of the computational domain is set to a convection boundary condition, allowing the fluid to flow smoothly out of the computational domain to avoid artificial reflection effects. At the same time, a preset wall function is set in the ground boundary of the computational domain to simulate the interaction between the ground and the airflow, accurately capturing the velocity gradient and turbulence characteristics within the boundary layer, thereby completing the first boundary setting operation.

[0084] After completing the first boundary setting operation, the first sector is simulated using the large eddy simulation (LES) method and the target terrain grid through the preset computational fluid dynamics software. The large eddy simulation method can directly analyze large-scale turbulent structures and has high simulation accuracy for phenomena such as flow separation, vortex generation and evolution under complex terrain. Combined with the encryption processing of the target terrain grid, it can more accurately capture the details of the flow field, especially near the wind turbine and in areas with complex terrain. Based on this simulation process, the three-dimensional wind velocity field and turbulence field with a preset time resolution within the swept area of ​​the wind turbine rotor are obtained through transient solution, focusing on capturing the terrain-induced flow separation and vortex structure, thereby determining the corresponding first simulation results.

[0085] When simulating the second sector, the boundary condition settings are relatively simple. The current wind speed profile is determined based on the current wind measurement data. This profile is typically a uniform inflow or a simple gradient distribution. The inlet boundary of the computational domain is set to the current wind speed profile for the uniform inflow. The outlet boundary is also set to a convection boundary condition. The ground boundary is set to a wall function to complete the second boundary setting.

[0086] After completing the second boundary setup, the second sector is simulated using the preset computational fluid dynamics software using the Reynolds-averaged turbulence model (RANS, a model based on numerical simulation methods for turbulence) and the target terrain grid. For sectors not affected by high slopes, the flow field structure is relatively simple, and the Reynolds-averaged turbulence model can quickly obtain simulation results while maintaining a certain level of accuracy. Based on this simulation process, a steady-state solution is used to obtain the average wind speed gradient data within the wind turbine rotor plane, as well as the baseline wind speed and turbulence intensity distribution at the wind turbine hub, thereby determining the corresponding second simulation results.

[0087] Furthermore, based on the first and second simulation results, the characteristics of different sectors are comprehensively considered to determine the target simulation result. The target simulation result includes the flow field characteristics of the entire target wind farm under different wind direction conditions, providing a basis for the subsequent determination of wind condition parameters.

[0088] When determining wind parameters for wind turbines and wind towers, a wind direction interpolation method is used to improve parameter accuracy. Specifically, a preset wind direction angle range is determined based on the first central wind direction angle of the third sector and the second central wind direction angle of the fourth sector, which are adjacent to each other in the target wind farm. A wind direction angle is randomly selected from the preset wind direction angle range as the current inflow wind direction angle, representing a specific wind direction in the actual wind farm.

[0089] Based on the target simulation results, the first simulated wind speed of the wind tower at the first central wind direction angle and the second simulated wind speed of the wind tower at the second central wind direction angle are determined respectively. The first target wind speed of the wind tower is determined using these two simulated wind speed values ​​and the current inflow wind direction angle. The calculation formula can be expressed as: Where, is the current wind direction angle, is the first center angle, is the first simulated wind speed of the wind tower at the first central wind direction angle, is the second center angle, is the second simulated wind speed of the wind tower under the second central wind direction angle, is the weighting coefficient, and The calculation formula is .

[0090] Similarly, the first simulated wind direction of the wind tower at the first central wind direction angle and the second simulated wind direction of the wind tower at the second central wind direction angle are determined based on the target simulation results. The first target wind direction of the wind tower is determined using these two simulated wind direction values ​​and the current inflow wind direction angle. The calculation formula can be expressed as: Where, is the first simulated wind direction of the wind tower under the first central wind direction angle, is the second simulated wind direction of the wind tower at the second central wind direction angle.

[0091] For wind turbines, the same method is used. Based on the target simulation results, the third simulated wind speed of the wind turbine at the first central wind direction angle and the fourth simulated wind speed of the wind turbine at the second central wind direction angle are determined. Using these two simulated wind speed values ​​and the current inflow wind direction angle, the second target wind speed of the wind turbine is determined. The calculation formula can be expressed as: Where, is the third simulated wind speed of the wind turbine under the first central wind direction angle, is the fourth simulated wind speed of the wind turbine generator system under the second central wind direction angle.

[0092] Similarly, based on the target simulation results, the third simulated wind direction of the wind turbine at the first central wind direction angle and the fourth simulated wind direction of the wind turbine at the second central wind direction angle are determined. Using these two simulated wind direction values ​​and the current inflow wind direction angle, the second target wind direction of the wind turbine is determined. The calculation formula for the second target wind direction can be expressed as: Where, is the second target wind direction, is the third simulated wind direction of the wind turbine under the first central wind direction angle, It is the fourth simulated wind direction of the wind turbine generator system under the second central wind direction angle.

[0093] In order to further improve the accuracy of the wind speed of the wind turbine, the wind acceleration factor is determined based on the first target wind speed and the second target wind speed. The calculation formula of the wind acceleration factor A is: The wind acceleration factor reflects the difference in wind speed between the wind turbine location and the wind tower location, which is mainly caused by factors such as terrain and obstacles. The wind acceleration factor and the current wind speed are used to correct the second target wind speed to obtain the corrected wind speed. Specifically, it can be expressed as , where is the current measured wind speed. In other words, the wind acceleration factor is multiplied by the current wind speed, and then added to the second target wind speed to obtain the final corrected wind speed.

[0094] Based on the first target wind speed, first target wind direction, corrected wind speed, and second target wind direction determined above, a comprehensive analysis is performed to determine the complete wind parameters for the wind turbines and wind towers in the target wind farm. These parameters, including wind speed, turbulence intensity, wind direction, and wind shear in the first sector, accurately reflect the actual wind characteristics of the target wind farm on high-slope terrain.

[0095] Step S14: constructing a corresponding first turbulent wind field model and a second turbulent wind field model for the first sector and the second sector respectively based on the wind condition parameters, and determining the target load of the wind turbine generator set using the first turbulent wind field model and the second turbulent wind field model.

[0096] In this embodiment, after obtaining the wind parameters of the wind turbine and the wind tower, it is necessary to construct turbulent wind field models for different sector characteristics, and then determine the target load of the wind turbine. For the first sector, the construction process is as follows: Based on the wind parameters, the wind wheel surface wind speed distribution matrix and the turbulence intensity distribution matrix of the wind turbine in the first sector are constructed. The wind wheel surface wind speed distribution matrix reflects the wind speed at different positions within the wind wheel swept surface, and the turbulence intensity distribution matrix reflects the distribution of turbulence intensity on the wind wheel surface. Then, the wind wheel surface wind speed distribution matrix and the turbulence intensity distribution matrix are used to determine the wind wheel surface power spectrum density matrix. The power spectrum density matrix describes the energy distribution of wind speed fluctuations at different frequencies, and is an important parameter for characterizing turbulent characteristics. Among them, the calculation formula for the power spectrum density value of any point in the wind wheel surface power spectrum density matrix is ​​as follows:

[0097] ;

[0098] Where, is the power spectrum density value of any point in the power spectrum density matrix of the wind wheel surface, is the value of any point in the turbulence intensity distribution matrix, is the value of any point in the wind speed distribution matrix, is the integral scale constant, is the frequency, and the value of this frequency is generally 0.01-1Hz.

[0099] In order to convert the frequency domain data into time domain data, the preset time step is used to capture the wind frequency domain data from the power spectrum density matrix of the wind wheel surface. Among them, the preset time step can be 0.05 seconds. The wind frequency domain data is converted into wind time domain data through the Fourier transform method, that is, the energy distribution of wind speed fluctuations at different frequencies is converted into a wind speed sequence that changes with time. This time domain data can more intuitively reflect the real-time effect of the turbulent wind field on the wind turbine, and generate a high-slope turbulent wind field with anisotropy of wind speed and turbulence intensity. Among them, any grid point in the turbulent wind field ,time The wind speed calculation formula is as follows:

[0100] ;

[0101] Where, is a frequency variable.

[0102] Furthermore, a first turbulent wind field model for the first sector was constructed based on the time-domain wind data. This model takes a time-domain wind speed series as input and simulates a turbulent wind field with specific statistical characteristics. This model takes into account the spatial and temporal correlation of wind speed, as well as parameters such as turbulence intensity and turbulence integral scale, and accurately reproduces the complex turbulent wind field characteristics under high-slope terrain.

[0103] For the second sector, a second turbulent wind field model for the wind turbine in the second sector is constructed based on wind condition parameters and preset dynamic simulation software. This process utilizes the wind speed, turbulence intensity, wind direction, and wind shear exponent at the turbine hub height, as defined in the wind condition parameters. Using the wind generation module of the whole-machine simulation dynamics software, a turbulent wind field with isotropic turbulence intensity can be generated in accordance with IEC 61400-1 (an international standard for wind turbine design) to construct the corresponding second turbulent wind field model.

[0104] After constructing the turbulent wind field models for both sectors, the target loads for each sector of the wind turbine were determined using the first and second turbulent wind field models via a pre-set dynamic simulation platform. Furthermore, due to the complexity of the terrain and the nonlinearity of the spatial variation of wind speed, precise interpolation calculations were required to calculate the wind speed at any characteristic point on the blade segment under conditions of blade rotation and large deformation.

[0105] During the simulation, the first and second turbulent wind field models were applied to corresponding sectors, simulating the forces acting on the wind turbine under different wind directions. By solving the wind turbine's equations of motion and structural mechanics, response parameters such as stress, strain, and displacement of each component were obtained. Statistical analysis of these response parameters ultimately determined the target load for the wind turbine. This target load is an important indicator for assessing the safety and reliability of the wind turbine.

[0106] As can be seen from the above, in this application, the current wind speed and wind direction time series of the wind farm are first obtained to determine the wind measurement data, and at the same time, a terrain physical model is constructed based on the high slope terrain information, and the terrain physical model is divided into an unstructured grid to obtain the target terrain grid; then, in each sector of the target wind farm, according to the preset sector step size and terrain information, the first sector affected by the high slope terrain and the second sector not affected by it are determined; then, based on the target terrain grid, the first sector and the second sector are simulated respectively to obtain the corresponding target simulation results, and then these simulation results and wind measurement data are used to determine the wind condition parameters of the wind turbine and the wind tower in the target wind farm; finally, according to the wind condition parameters, the first turbulent wind field model and the second turbulent wind field model are constructed for the first sector and the second sector respectively, and the target load of the wind turbine is determined with the help of these two models. In this way, the present application can improve the calculation efficiency and accuracy of the wind turbine load under complex terrain, which is beneficial to the safety assessment of high slope wind turbines.

[0107] See also Figure 2 As shown, in order to more accurately evaluate the safety and reliability of high-slope wind turbines under different operating conditions after determining the target load, the present application further provides a specific method for determining the load of high-slope wind turbines, including:

[0108] Step S21: Obtain the maximum load under each wind speed condition from the target load, determine the probability distribution function of the maximum load using the Gumbel distribution method, and determine the ultimate load of the first sector based on the probability distribution function.

[0109] In this embodiment, after obtaining the target load, the maximum load under each wind speed condition must first be extracted from the target load. Wind speed conditions are typically divided based on the actual operation and design requirements of the wind farm. For example, the wind speed range can be divided into multiple intervals, each corresponding to a specific wind speed condition. For each wind speed condition, the maximum load under that condition is filtered from the target load data.

[0110] In order to accurately describe the probability distribution characteristics of these maximum loads, the Gumbel distribution method (i.e., the Gumbel distribution method, a common method for calculating "once-in-a-year" ocean hydrological and meteorological elements using a theoretically based frequency curve) is used to determine the probability distribution function of the maximum load. The calculation formula of the probability distribution function can be Where, is the maximum load, is the fitting coefficient for the rate of change of the function, is the fitting coefficient for the function distribution position.

[0111] Based on the resulting probability distribution function, the ultimate load of the first sector is determined. Ultimate loads are crucial for the structural design and safety assessment of wind turbines. This ultimate load value provides key parameters for wind turbine structural design, ensuring maximum safety under extreme operating conditions.

[0112] Step S22: Obtain the number of wind speed intervals corresponding to the sectors in the wind frequency matrix, determine the load cycle frequency of the current load in the wind speed interval and the duration of the current wind speed in the wind speed interval, and determine the maximum number of cycles of the current load based on the preset load spectrum.

[0113] In this embodiment, a preset sector step size is determined, and the ratio of the wind direction range to the preset sector step size is determined as the sector number. The wind measurement data is divided based on the sector number, and the frequency of occurrence of different wind speed intervals in each sector is counted to form a wind frequency matrix.

[0114] Next, a load is randomly selected from the target loads as the current load, and fatigue load analysis begins. In the wind frequency matrix, the interval number for the wind speed interval corresponding to the current sector is obtained. For each wind speed interval, the load cycle frequency of the current load is determined. The load cycle frequency refers to the number of load fluctuations per unit time and can be determined through statistical analysis of the load time series.

[0115] The duration of the current wind speed within this wind speed range is also determined. The maximum number of cycles for the current load within this wind speed range is then determined based on a preset load spectrum, combined with the duration of the wind speed range and the load cycle frequency. The preset load spectrum specifies the fatigue strength of different load types at different cycle times. By querying the load spectrum and performing interpolation calculations, the maximum number of cycles corresponding to the current load can be determined.

[0116] Step S23: determining a fatigue load based on the number of sectors, the number of intervals, the load cycle frequency, the duration, and the maximum number of cycles.

[0117] In this embodiment, fatigue loads are determined by combining wind farm sector and wind speed range characteristics. First, all sectors of the target wind farm and the wind speed ranges within each sector are traversed. For each sector and wind speed range, corresponding parameters are extracted: the duration of the wind speed range is obtained from wind condition statistics, the load cycle frequency is determined by analyzing the load time series, and the maximum number of cycles under this load is determined based on a preset load spectrum.

[0118] Next, the damage value for each sector-wind speed interval combination is calculated by multiplying the duration of the wind speed interval by the load cycle frequency and dividing the result by the maximum number of cycles. The damage values ​​for all wind speed intervals within the same sector are summed to obtain the cumulative damage for that sector. Next, taking into account the frequency of sector occurrence, the cumulative damage for that sector is multiplied by the probability of the wind direction occurring in that sector to determine the contribution of that sector to the total fatigue damage.

[0119] Finally, the contributions of all sectors are summed to obtain the total fatigue load. This process comprehensively considers the wind direction, wind speed distribution, and load cycle characteristics of the wind field. By calculating the fatigue damage accumulation under various working conditions, it provides a key basis for the fatigue life assessment and structural optimization design of wind turbines, ensuring the long-term operation reliability of the units in complex environments with high slopes. The fatigue load calculation formula is as follows:

[0120] ;

[0121] Where, is the number of sectors, is the number of wind speed intervals in the sector, is the duration of time the current wind speed is in this wind speed range, is the load cycle frequency, is the maximum number of cycles.

[0122] As can be seen from the above, in this application, after determining the target load, the maximum load of each wind speed condition is first extracted from the target load, the probability distribution function is fitted using the Gumbel distribution, and the limit load of the first sector is extrapolated to lay a solid foundation for the structural safety design of the unit under extreme working conditions; then the sectors are divided to construct a wind frequency matrix, the current load is selected, and parameters such as the number of wind speed intervals, load cycle frequency, and duration are obtained. The maximum number of cycles is determined in combination with the preset load spectrum; finally, the sectors and wind speed intervals are traversed, and the parameters are extracted to calculate the damage values ​​of each combination, and the total fatigue load is accumulated. The entire process, from extreme working conditions to long-term fatigue, comprehensively considers the characteristics of high-slope wind fields, providing both extreme load boundaries for unit structural design and accurate quantification of long-term fatigue damage, helping to optimize unit structure and evaluate lifespan. It can, to a certain extent, ensure the long-term reliable operation of wind turbines in high-slope environments and improve the overall safety of wind power projects.

[0123] Accordingly, see Figure 3As shown, an embodiment of the present application provides a load determination device for a high-slope wind turbine generator system, comprising:

[0124] A gridding module 11 is used to obtain the current wind speed and wind direction time series of the target wind farm located on the high slope terrain to determine the corresponding wind measurement data, construct a terrain physical model based on the terrain information of the high slope, and perform unstructured grid division on the terrain physical model to obtain the target terrain grid;

[0125] a sector determination module 12, configured to determine, in each sector of the target wind farm, a first sector affected by the high side slope terrain and a second sector not affected by the high side slope terrain based on a preset sector step size and the terrain information;

[0126] a wind condition determination module 13, configured to simulate the first sector and the second sector respectively based on the target terrain grid to obtain corresponding target simulation results, and determine wind condition parameters of the wind turbines and wind towers in the target wind farm using the target simulation results and the wind measurement data;

[0127] The load determination module 14 is used to construct a corresponding first turbulent wind field model and a second turbulent wind field model for the first sector and the second sector respectively based on the wind condition parameters, and use the first turbulent wind field model and the second turbulent wind field model to determine the target load of the wind turbine group.

[0128] As can be seen from the above, in this application, the current wind speed and wind direction time series of the wind farm are first obtained to determine the wind measurement data, and at the same time, a terrain physical model is constructed based on the high slope terrain information, and the terrain physical model is divided into an unstructured grid to obtain the target terrain grid; then, in each sector of the target wind farm, according to the preset sector step size and terrain information, the first sector affected by the high slope terrain and the second sector not affected by it are determined; then, based on the target terrain grid, the first sector and the second sector are simulated respectively to obtain the corresponding target simulation results, and then these simulation results and wind measurement data are used to determine the wind condition parameters of the wind turbine and the wind tower in the target wind farm; finally, according to the wind condition parameters, the first turbulent wind field model and the second turbulent wind field model are constructed for the first sector and the second sector respectively, and the target load of the wind turbine is determined with the help of these two models. In this way, the present application can improve the calculation efficiency and accuracy of the wind turbine load under complex terrain, which is beneficial to the safety assessment of high slope wind turbines.

[0129] In some specific implementations, the grid division module 11 specifically includes:

[0130] a spatial alignment unit, configured to obtain a satellite image and a measured topographic map of the high slope, correct the satellite image based on a construction coordinate system of the high slope, and spatially align the corrected image with the measured topographic map;

[0131] an information labeling unit, configured to, after spatial alignment, use a spatial interpolation algorithm to complete the elevation data of the missing areas in the measured topographic map to generate a first topographic map, and label the topographic feature points and the coordinates of the wind turbine generator set in the first topographic map to obtain a second topographic map;

[0132] a topographic map adjustment unit, configured to adjust the second topographic map using the design parameters of the high slope construction platform and preset three-dimensional modeling software to obtain a third topographic map;

[0133] The first model construction unit is used to determine a preset horizontal sampling interval and a preset vertical sampling interval as modeling parameters, and to construct a terrain physical model based on the third topographic map and the modeling parameters using the preset three-dimensional modeling software; the sampling interval is the interval at which the preset three-dimensional modeling software performs terrain sampling based on the third topographic map.

[0134] In some specific implementations, the grid division module 11 specifically includes:

[0135] A first grid generation unit is configured to perform unstructured grid division on the terrain physical model using preset computational fluid dynamics software to obtain a first terrain grid;

[0136] a second grid generation unit, configured to densify the terrain that meets a preset complexity condition in the first terrain grid based on a global grid size control parameter, and generate a preset number of prismatic layer grids in a surface boundary layer in the first terrain grid to obtain a second terrain grid;

[0137] The third grid generation unit is configured to establish a target encrypted area in the second terrain grid with the wind turbine as the center, and encrypt the target encrypted area based on a preset grid encryption strategy to obtain a target terrain grid.

[0138] In some specific implementations, the wind condition determination module 13 specifically includes:

[0139] A calculation domain determining unit, configured to determine the calculation domain of the preset computational fluid dynamics software;

[0140] A first boundary setting unit is configured to, when simulating the first sector, set the inlet boundary of the computational domain to a fluctuating wind speed profile including atmospheric boundary layer characteristics, set the outlet boundary of the computational domain to convection, and set a preset wall function in the ground boundary of the computational domain to complete a first boundary setting operation;

[0141] a first sector simulation unit configured to simulate the first sector using a large eddy simulation method and the target terrain grid using the preset computational fluid dynamics software after completing the first boundary setting operation, to obtain a three-dimensional wind velocity field and turbulence field with a preset time resolution within a swept surface of a wind turbine rotor, so as to determine a corresponding first simulation result;

[0142] a second boundary setting unit, configured to determine a current wind speed profile based on current wind measurement data when simulating the second sector, and set an inlet boundary of the calculation domain to the current wind speed profile of uniform inflow to complete a second boundary setting operation;

[0143] a first sector simulation unit, configured to, after completing the second boundary setting operation, simulate the second sector using a preset Reynolds time-averaged turbulence model and the target terrain grid through the preset computational fluid dynamics software to obtain average wind speed gradient data within the wind turbine rotor plane and a reference wind speed and turbulence intensity distribution at the wind turbine hub, so as to determine a corresponding second simulation result;

[0144] A simulation result determining unit is configured to determine a target simulation result based on the first simulation result and the second simulation result.

[0145] In some specific implementations, the wind condition determination module 13 specifically includes:

[0146] an angle determination unit, configured to determine a preset wind direction angle range based on a first central wind direction angle of a third sector and a second central wind direction angle of a fourth sector, and randomly select a wind direction angle from the preset wind direction angle range as a current inflow wind direction angle; wherein the third sector and the fourth sector are adjacent to each other in the target wind farm;

[0147] a first simulated wind speed determining unit, configured to determine, based on the target simulation result, a first simulated wind speed of the wind measuring tower at the first central wind direction angle and a second simulated wind speed of the wind measuring tower at the second central wind direction angle;

[0148] a first target wind speed determining unit, configured to determine a first target wind speed of the wind measurement tower by using the first simulated wind speed, the second simulated wind speed, and the current inflow wind direction angle;

[0149] a first simulated wind direction determining unit, configured to determine, based on the target simulation result, a first simulated wind direction of the wind measurement tower at the first central wind direction angle and a second simulated wind direction of the wind measurement tower at the second central wind direction angle;

[0150] a first target wind direction determining unit, configured to determine a first target wind direction of the wind measurement tower by using the first simulated wind direction, the second simulated wind direction, and the current inflow wind direction angle;

[0151] a second simulated wind speed determining unit, configured to determine, based on the target simulation result, a third simulated wind speed of the wind turbine generator set at the first central wind direction angle and a fourth simulated wind speed of the wind turbine generator set at the second central wind direction angle;

[0152] a second target wind speed determining unit, configured to determine a second target wind speed of the wind turbine generator set by using the third simulated wind speed, the fourth simulated wind speed, and the current inflow wind direction angle;

[0153] a second simulated wind direction determining unit, configured to determine, based on the target simulation result, a third simulated wind direction of the wind turbine generator set at the first central wind direction angle and a fourth simulated wind direction of the wind turbine generator set at the second central wind direction angle;

[0154] a second target wind direction determining unit, configured to determine a second target wind direction of the wind turbine generator set by using the third simulated wind direction, the fourth simulated wind direction, and the current inflow wind direction angle;

[0155] a wind speed correction unit, configured to determine a wind acceleration factor based on the first target wind speed and the second target wind speed, and correct the second target wind speed using the wind acceleration factor and the current wind speed to obtain a corrected wind speed;

[0156] A wind condition parameter determination unit is used to determine corresponding wind condition parameters based on the first target wind speed, the first target wind direction, the corrected wind speed and the second target wind direction.

[0157] In some specific implementations, the load determination module 14 specifically includes:

[0158] a first matrix determination unit, configured to construct a rotor surface wind speed distribution matrix and a turbulence intensity distribution matrix for the wind turbine in the first sector based on the wind condition parameters, and determine a rotor surface power spectrum density matrix using the rotor surface wind speed distribution matrix and the turbulence intensity distribution matrix;

[0159] A data conversion unit is used to capture wind frequency domain data from the wind rotor surface power spectrum density matrix using a preset time step, and convert the wind frequency domain data into wind time domain data using a Fourier transform method;

[0160] A second model building unit is configured to build a first turbulent wind field model of the first sector based on the wind condition time domain data;

[0161] a third model building unit, configured to build a second turbulent wind field model for the wind turbine in the second sector based on the wind condition parameters and preset dynamics simulation software;

[0162] The target load determination unit is configured to determine the target load of the wind turbine generator set by using the first turbulent wind field model and the second turbulent wind field model through a preset dynamics simulation platform.

[0163] In some specific implementations, the load determination module 14 further includes:

[0164] a limit load determination unit, configured to obtain a maximum load under each wind speed condition from the target load, determine a probability distribution function of the maximum load using a Gumbel distribution method, and determine a limit load of the first sector based on the probability distribution function;

[0165] a second matrix determination unit, configured to determine a preset sector step size, determine a ratio of a wind direction range to the preset sector step size as a sector number, and divide the wind measurement data based on the sector number to determine a wind frequency matrix;

[0166] A load determination unit, configured to randomly select a load from the target loads as a current load;

[0167] a number determination unit, configured to obtain the number of wind speed intervals corresponding to sectors in the wind frequency matrix, determine the load cycle frequency of the current load in the wind speed interval and the duration of the current wind speed in the wind speed interval, and determine the maximum number of cycles of the current load based on a preset load spectrum;

[0168] A fatigue load determination unit is configured to determine a fatigue load based on the number of sectors, the number of intervals, the load cycle frequency, the duration, and the maximum number of cycles.

[0169] Furthermore, the embodiment of the present application also discloses an electronic device, Figure 4This is a structural diagram of an electronic device 20 according to an exemplary embodiment. The content in the diagram cannot be considered as any limitation on the scope of use of this application. The electronic device 20 may specifically include: at least one processor 21, at least one memory 22, a power supply 23, a communication interface 24, an input / output interface 25, and a communication bus 26. The memory 22 is used to store a computer program, which is loaded and executed by the processor 21 to implement the relevant steps in the load determination method for a high-slope wind turbine disclosed in any of the aforementioned embodiments. In addition, the electronic device 20 in this embodiment may specifically be an electronic computer.

[0170] In this embodiment, the power supply 23 is used to provide operating voltage for each hardware device on the electronic device 20; the communication interface 24 can create a data transmission channel between the electronic device 20 and the external device. The communication protocol it follows is any communication protocol that can be applied to the technical solution of this application and is not specifically limited here; the input and output interface 25 is used to obtain external input data or output data to the outside world. Its specific interface type can be selected according to specific application needs and is not specifically limited here.

[0171] In addition, the memory 22, as a carrier for resource storage, can be a read-only memory, random access memory, disk or CD, etc. The resources stored thereon can include an operating system 221, a computer program 222, etc., and the storage method can be temporary storage or permanent storage.

[0172] The operating system 221 is used to manage and control the hardware devices on the electronic device 20 and the computer program 222. The operating system 221 can be Windows Server, NetWare, Unix, Linux, etc. In addition to including a computer program capable of implementing the load determination method for a high-slope wind turbine generator system executed by the electronic device 20 as disclosed in any of the aforementioned embodiments, the computer program 222 can further include computer programs capable of implementing other specific tasks.

[0173] Furthermore, this application discloses a computer-readable storage medium for storing a computer program. When executed by a processor, the computer program implements the aforementioned method for determining the load of a high-slope wind turbine. The specific steps of this method can be found in the corresponding contents disclosed in the aforementioned embodiments and will not be further described here.

[0174] The various embodiments in this specification are described in a progressive manner, with each embodiment focusing on its differences from the other embodiments. Reference can be made to the descriptions of the identical or similar parts between the various embodiments. For the devices disclosed in the embodiments, since they correspond to the methods disclosed in the embodiments, the descriptions are relatively simple, and the relevant parts can be referred to the descriptions of the methods.

[0175] Professionals may further appreciate that the units and algorithm steps of each example described in conjunction with the embodiments disclosed herein can be implemented in electronic hardware, computer software, or a combination of the two. In order to clearly illustrate the interchangeability of hardware and software, the above description has generally described the components and steps of each example according to their functions. Whether these functions are performed in hardware or software depends on the specific application and design constraints of the technical solution. Professionals and technicians may use different methods to implement the described functions for each specific application, but such implementation should not be considered beyond the scope of this application.

[0176] The steps of the methods or algorithms described in conjunction with the embodiments disclosed herein may be implemented directly using hardware, a software module executed by a processor, or a combination of the two. The software module may be placed in random access memory (RAM), internal memory, read-only memory (ROM), electrically programmable ROM, electrically erasable programmable ROM, registers, a hard disk, a removable disk, a CD-ROM, or any other form of storage medium known in the art.

[0177] Finally, it should be noted that, in this document, relational terms such as first and second, etc., are used only 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 terms "comprises," "comprising," or any other variations thereof are intended to cover non-exclusive inclusion, such that a process, method, article, or device comprising a series of elements includes not only those elements, but also other elements not explicitly listed, or elements inherent to such process, method, article, or device. In the absence of further limitations, an element defined by the phrase "comprising a ..." does not exclude the presence of additional identical elements in the process, method, article, or device comprising the element.

[0178] The above is a detailed introduction to the technical solution provided by the present application. Specific examples are used herein to illustrate the principles and implementation methods of the present application. The description of the above embodiments is only used to help understand the method of the present application and its core idea. At the same time, for those skilled in the art, according to the ideas of the present application, there may be changes in the specific implementation methods and application scope. In summary, the content of this specification should not be understood as a limitation on the present application.

Claims

1. A method for determining the load of a high-slope wind turbine generator system, characterized in that: include: Obtaining the current wind speed and direction time series of a target wind farm located on a high slope terrain to determine the corresponding wind measurement data, constructing a terrain physical model based on the high slope terrain information, and performing unstructured grid division on the terrain physical model to obtain a target terrain grid; In each sector of the target wind farm, determining a first sector affected by the high side slope terrain and a second sector not affected by the high side slope terrain based on a preset sector step size and the terrain information; Simulating the first sector and the second sector respectively based on the target terrain grid to obtain corresponding target simulation results, and determining wind condition parameters of the wind turbines and wind towers in the target wind farm using the target simulation results and the wind measurement data; A first turbulent wind field model and a second turbulent wind field model are constructed for the first sector and the second sector respectively based on the wind condition parameters, and a target load of the wind turbine is determined using the first turbulent wind field model and the second turbulent wind field model.

2. The method for determining the load of a high-slope wind turbine generator system according to claim 1, characterized in that: The construction of a terrain physical model based on the terrain information of the high slope includes: Acquiring a satellite image and a measured topographic map of the high slope, correcting the satellite image based on a construction coordinate system of the high slope, and spatially aligning the corrected image with the measured topographic map; After the spatial alignment, a spatial interpolation algorithm is used to complete the elevation data of the missing areas in the measured topographic map to generate a first topographic map, and the topographic feature points and the coordinates of the wind turbine generator set are marked in the first topographic map to obtain a second topographic map; Adjusting the second topographic map using the design parameters of the high slope construction platform and preset three-dimensional modeling software to obtain a third topographic map; The preset horizontal sampling interval and the preset vertical sampling interval are determined as modeling parameters, and a terrain physical model is constructed based on the third topographic map and the modeling parameters through the preset three-dimensional modeling software; the sampling interval is the interval at which the preset three-dimensional modeling software performs terrain sampling based on the third topographic map.

3. The method for determining the load of a high-slope wind turbine generator system according to claim 1, characterized in that: The performing unstructured grid division on the terrain physical model to obtain a target terrain grid includes: Performing unstructured grid division on the terrain physical model using preset computational fluid dynamics software to obtain a first terrain grid; Densifying the terrain that meets the preset complexity condition in the first terrain grid based on the global grid size control parameter, and generating a preset number of prismatic layer grids in the surface boundary layer of the first terrain grid to obtain a second terrain grid; In the second terrain grid, a target encryption area is established with the wind turbine generator set as the center, and the target encryption area is encrypted based on a preset grid encryption strategy to obtain a target terrain grid.

4. The method for determining the load of a high-slope wind turbine generator system according to claim 1, characterized in that: The step of constructing a corresponding first turbulent wind field model and a second turbulent wind field model for the first sector and the second sector based on the wind condition parameters, and determining a target load of the wind turbine generator set by using the first turbulent wind field model and the second turbulent wind field model, includes: constructing a rotor surface wind speed distribution matrix and a turbulence intensity distribution matrix for the wind turbine in the first sector based on the wind condition parameters, and determining a rotor surface power spectrum density matrix using the rotor surface wind speed distribution matrix and the turbulence intensity distribution matrix; Capturing wind frequency domain data from the wind rotor surface power spectrum density matrix using a preset time step, and converting the wind frequency domain data into wind time domain data using a Fourier transform method; Constructing a first turbulent wind field model of the first sector based on the wind condition time domain data; Constructing a second turbulent wind field model for the wind turbine in the second sector based on the wind condition parameters and preset dynamic simulation software; The target load of the wind turbine generator set is determined by using the first turbulent wind field model and the second turbulent wind field model through a preset dynamic simulation platform.

5. The method for determining the load of a high-slope wind turbine generator system according to any one of claims 1 to 4, characterized in that: After determining the target load of the wind turbine generator set by using the first turbulent wind field model and the second turbulent wind field model, the method further includes: Obtaining a maximum load under each wind speed condition from the target load, determining a probability distribution function of the maximum load using a Gumbel distribution method, and determining a limit load of the first sector based on the probability distribution function; Determining a preset sector step size, determining a ratio of a wind direction range to the preset sector step size as a sector number, and dividing the wind measurement data based on the sector number to determine a wind frequency matrix; Randomly select a load from the target load as the current load; Obtaining the number of wind speed intervals corresponding to the sectors in the wind frequency matrix, determining the load cycle frequency of the current load in the wind speed interval and the duration of the current wind speed in the wind speed interval, and determining the maximum number of cycles of the current load based on a preset load spectrum; A fatigue load is determined based on the number of sectors, the number of intervals, the load cycle frequency, the duration, and the maximum number of cycles.

6. A load determination device for a high-slope wind turbine generator system, characterized in that: include: A gridding module is used to obtain the current wind speed and wind direction time series of the target wind farm located on the high slope terrain to determine the corresponding wind measurement data, build a terrain physical model based on the terrain information of the high slope, and perform unstructured grid division on the terrain physical model to obtain the target terrain grid; a sector determination module, configured to determine, in each sector of the target wind farm, a first sector affected by the high side slope terrain and a second sector not affected by the high side slope terrain based on a preset sector step size and the terrain information; a wind condition determination module, configured to simulate the first sector and the second sector respectively based on the target terrain grid to obtain corresponding target simulation results, and determine wind condition parameters of the wind turbines and wind towers in the target wind farm using the target simulation results and the wind measurement data; A load determination module is used to construct a corresponding first turbulent wind field model and a second turbulent wind field model for the first sector and the second sector respectively based on the wind condition parameters, and use the first turbulent wind field model and the second turbulent wind field model to determine the target load of the wind turbine group.

7. The load determination device for a high-slope wind turbine generator system according to claim 6, characterized in that: The grid division module includes: a spatial alignment unit, configured to obtain a satellite image and a measured topographic map of the high slope, correct the satellite image based on a construction coordinate system of the high slope, and spatially align the corrected image with the measured topographic map; an information labeling unit, configured to, after spatial alignment, use a spatial interpolation algorithm to complete the elevation data of the missing areas in the measured topographic map to generate a first topographic map, and label the topographic feature points and the coordinates of the wind turbine generator set in the first topographic map to obtain a second topographic map; a topographic map adjustment unit, configured to adjust the second topographic map using the design parameters of the high slope construction platform and preset three-dimensional modeling software to obtain a third topographic map; The first model construction unit is used to determine a preset horizontal sampling interval and a preset vertical sampling interval as modeling parameters, and to construct a terrain physical model based on the third topographic map and the modeling parameters using the preset three-dimensional modeling software; the sampling interval is the interval at which the preset three-dimensional modeling software performs terrain sampling based on the third topographic map.

8. The load determination device for a high-slope wind turbine generator system according to claim 6, characterized in that: The grid division module includes: A first grid generation unit is configured to perform unstructured grid division on the terrain physical model using preset computational fluid dynamics software to obtain a first terrain grid; a second grid generation unit, configured to densify the terrain that meets a preset complexity condition in the first terrain grid based on a global grid size control parameter, and generate a preset number of prismatic layer grids in a surface boundary layer in the first terrain grid to obtain a second terrain grid; The third grid generation unit is configured to establish a target encrypted area in the second terrain grid with the wind turbine as the center, and encrypt the target encrypted area based on a preset grid encryption strategy to obtain a target terrain grid.

9. An electronic device, characterized in that: include: Memory, used to store computer programs; A processor is configured to execute the computer program to implement the load determination method for a high-slope wind turbine generator system according to any one of claims 1 to 5.

10. A computer-readable storage medium, characterized in that Used to store a computer program; wherein, when the computer program is executed by a processor, the load determination method for a high-slope wind turbine generator set according to any one of claims 1 to 5 is implemented.

Images