Water turbine parameter determination method, device and equipment based on sea area resource characteristics and storage medium

Abstract

The invention discloses a water turbine parameter determination method, device and equipment based on sea area resource characteristics and a storage medium. The method comprises the steps that historical year-round water depth data, tide level data and sea area shoreline data of a to-be-developed sea area are acquired; obtaining a sea area water depth initial value according to the water depth data; obtaining the diameter of the water turbine according to the sea water depth initial value, the sea surface safety distance and the seabed safety operation distance; according to a preset tidal current simulation model, obtaining an annual average flow velocity at a preset water depth and a corresponding hour number; target annual average flow rates are extracted, and the minimum annual average flow rate in the target annual average flow rates serves as the rated flow rate of the water turbine; calculating the initial rated power of the water turbine according to the rated flow velocity; and finally determining the chord length and the torsion angle of the turbine blade according to the turbine diameter, the rated flow velocity and the initial rated power. By implementing the method, the rated power of the water turbine can be determined according to local sea area conditions, and specific water turbine parameters are obtained.

Description

Technical Field

[0001] The present invention relates to the technical field of water turbine parameter determination, and in particular to a water turbine parameter determination method, device, equipment and storage medium based on sea area resource characteristics. Background Art

[0002] Tidal energy is a clean, renewable energy source that utilizes the flow energy of seawater. Tidal power generation is an emerging ocean power generation technology, and its core equipment is the tidal turbine. The design of tidal turbines needs to consider multiple factors, including flow velocity, water depth, etc., to ensure that they can have high power generation efficiency in the sea area where they are deployed. Therefore, optimization design technology plays a key role in the design of tidal turbines, so as to design the optimal turbine parameters according to the actual sea conditions and maximize the use of tidal energy resources.

[0003] The existing tidal energy turbine design scheme does not take into account the resource conditions of the local sea area. Usually, a total power is preset first, and then the specific turbine parameters are designed based on this goal. The advantage of this design scheme is that a set of turbine parameters can be quickly obtained, but its disadvantage is that the designed rated power may be too large or too small, and it cannot adapt to the resource conditions of a specific sea area, thus failing to maximize the development of the tidal energy resources in the sea area.

[0004] Therefore, it is an urgent problem to determine the rated power of the turbine according to the local sea conditions and then obtain the specific turbine parameters. Summary of the invention

[0005] The present invention provides a method, device, equipment and storage medium for determining turbine parameters based on sea resource characteristics, which can determine the rated power of the turbine according to local sea conditions and then obtain specific turbine parameters.

[0006] An embodiment of the present invention provides a method for determining turbine parameters based on marine resource characteristics, comprising:

[0007] Obtain the historical full-year water depth data, historical full-year tide level data, and historical full-year sea area coastline data of the sea area to be developed;

[0008] The theoretical lowest tide level is obtained based on the above water depth data, and the above theoretical lowest tide level is used as the initial value of the sea area water depth;

[0009] Obtaining the sea surface safety distance and the seabed safety operating distance, and calculating the turbine diameter according to the initial value of the sea area water depth, the sea surface safety distance and the seabed safety operating distance;

[0010] Input the water depth data, the sea area coastline data and the tide level data into a preset tidal current simulation model to obtain the annual average flow velocity at the preset water depth of the sea area to be developed in each simulated year and the number of hours in a year when the annual average flow velocity is reached;

[0011] Compare the number of hours corresponding to each simulated year with the preset number of hours, and extract the annual average flow rate of the simulated year with hours greater than the preset number of hours as the target annual average flow rate;

[0012] The smallest annual average flow rate among the above target annual average flow rates is used as the rated flow rate of the above turbine;

[0013] The initial rated power of the turbine is calculated based on the above rated flow rate;

[0014] The chord length and the twist angle of the turbine blades are determined according to the turbine diameter, the rated flow rate and the initial rated power.

[0015] Furthermore, the construction of the above preset power flow simulation model includes:

[0016] Acquire a number of data samples and a number of flow velocity data samples corresponding to the data samples; wherein the data samples include: shoreline data samples, water depth data samples and tide level data samples;

[0017] Inputting the above-mentioned data samples into the tidal current simulation model to be trained, so that the above-mentioned tidal current model to be trained simulates the flow velocity of the simulated year according to the above-mentioned tidal level data samples to obtain the simulated flow velocity corresponding to the simulated year;

[0018] According to the above simulated flow rate and the above several flow rate data samples, a relative error is calculated;

[0019] The relative error is compared with a preset relative error value. If the relative error is smaller than the preset relative error, the parameters of the tidal current simulation model to be trained are adjusted and the tidal current simulation model is continued to be trained; otherwise, a trained tidal current simulation model is obtained.

[0020] Furthermore, the initial rated power of the turbine is calculated based on the rated flow rate, including:

[0021] Obtaining several cross sections of the turbine blades;

[0022] Obtaining the initial chord length of the first section, the initial torsion angle of the first section, and the preset parameter strategy network;

[0023] Repeat the first parameter calculation operation until the initial chord lengths and initial torsion angles of all cross sections are obtained;

[0024] The initial rated power is calculated based on the initial chord length, initial twist angle and the rated flow velocity of each section;

[0025] The above parameter calculation operation includes:

[0026] Obtaining the initial chord length and initial torsion angle of the current section; wherein, when the current section is the first section, the initial chord length and initial torsion angle of the current section are the initial chord length and initial torsion angle of the first section;

[0027] According to the initial chord length of the current section, the initial torsion angle of the current section and the above-mentioned preset parameter strategy network, the initial chord length and initial torsion angle of the next section are calculated.

[0028] Furthermore, the initial rated power is calculated based on the initial chord length, initial twist angle and rated flow velocity of each section, including:

[0029] Obtain the distance from each section to the radius of the turbine blade, the turbine speed, and the turbine airfoil;

[0030] According to the distance, the rotation speed of the turbine and the rated flow rate, the initial inflow angle of each section and the initial relative velocity of each section are calculated;

[0031] According to the initial twist angle of each section and the corresponding initial incoming flow angle, the initial angle of attack of each section is calculated;

[0032] According to the airfoil and the initial angle of attack of each section, an initial lift coefficient of each section and an initial drag coefficient of each section are obtained;

[0033] The initial lift of each section and the initial drag of each section are calculated according to the initial relative velocity of each section, the initial chord length of each section, the initial lift coefficient of each section and the initial drag coefficient of each section;

[0034] The initial moment of each section is calculated according to the initial incoming flow angle of each section, each of the distances, the initial lift of each section and the initial drag of each section;

[0035] According to the number of turbine blades and the initial moment of each section, the initial total moment of the turbine is calculated;

[0036] The initial rated power of the turbine is calculated based on the initial total torque and the rotational speed.

[0037] Furthermore, the chord length and twist angle of the turbine blades are determined according to the turbine diameter, the rated flow rate and the initial rated power, including:

[0038] Repeating the parameter updating operation according to the turbine diameter, the rated flow rate and the initial rated power until the chord length and the twist angle of the turbine blade are determined;

[0039] The above parameter update operation includes:

[0040] Obtain the current rated power; wherein the initial rated power is the above initial rated power;

[0041] Calculate the current energy efficiency based on the turbine diameter, rated flow rate and current rated power;

[0042] Calculate the current instant reward based on the current chord length, current torsion angle and current energy efficiency of all cross sections;

[0043] According to the current instant reward and the chord lengths and torsion angles of all current sections, the updated chord length and torsion angle of the first section are calculated;

[0044] Determine whether the updated chord length and torsion angle of the first section converge, if so, recalculate the updated chord lengths and torsion angles of all other sections according to the updated chord length and torsion angle of the first section, and determine the chord length and torsion angle of the turbine blade;

[0045] Otherwise, based on the updated chord length and torsion angle of the first section, the updated chord lengths and torsion angles of all other sections are recalculated to obtain the updated chord lengths and torsion angles of all sections, and the updated rated power is calculated based on the updated chord lengths and torsion angles of all sections.

[0046] Based on the above method embodiment, the present invention provides a corresponding device embodiment;

[0047] The present invention provides a device for determining turbine parameters based on sea area resource characteristics, comprising:

[0048] Data acquisition module, water depth calculation module, diameter calculation module, hours determination module, target annual average flow rate determination module, rated flow rate determination module, initial rated power calculation module, and chord length and torsion angle calculation module;

[0049] The data acquisition module is used to obtain the historical water depth data, historical tide level data and historical coastline data of the sea area to be developed throughout the year;

[0050] The water depth calculation module is used to obtain the theoretical lowest tide level according to the water depth data, and use the theoretical lowest tide level as the initial value of the water depth of the sea area;

[0051] The diameter calculation module is used to obtain the sea surface safety distance and the seabed safety operating distance, and calculate the turbine diameter according to the initial value of the sea area water depth, the sea surface safety distance and the seabed safety operating distance;

[0052] The hour determination module is used to input the water depth data, the sea area coastline data and the tide level data into a preset tidal current simulation model to obtain the annual average flow velocity at the preset water depth of the sea area to be developed in each simulated year and the number of hours in a year to reach the annual average flow velocity;

[0053] The above-mentioned target annual average flow rate determination module is used to compare the number of hours corresponding to each simulated year with the preset number of hours, and extract the annual average flow rate of the simulated year with the number of hours greater than the preset number of hours as the target annual average flow rate;

[0054] The rated flow rate determination module is used to use the smallest annual average flow rate among the target annual average flow rates as the rated flow rate of the turbine;

[0055] The initial rated power calculation module is used to calculate the initial rated power of the turbine according to the rated flow rate;

[0056] The chord length and torsion angle calculation module is used to determine the chord length and torsion angle of the turbine blades according to the turbine diameter, the rated flow rate and the initial rated power.

[0057] Furthermore, the above-mentioned hour determination module includes:

[0058] Data sample acquisition unit, model training unit, relative error calculation unit and model determination unit;

[0059] The sample acquisition unit is used to acquire a number of data samples and a number of flow velocity data samples corresponding to the data samples; wherein the data samples include: shoreline data samples, water depth data samples and tide level data samples;

[0060] The model training unit is used to input the plurality of data samples into the tidal current simulation model to be trained, so that the tidal current model to be trained simulates the flow velocity of the simulation year according to the plurality of tidal level data samples to obtain the simulated flow velocity corresponding to the simulation year;

[0061] The relative error calculation unit is used to calculate the relative error based on the simulated flow rate and the plurality of flow rate data samples;

[0062] The above-mentioned model determination unit is used to compare the above-mentioned relative error with the preset relative error value. If the above-mentioned relative error is smaller than the above-mentioned preset relative error, the parameters of the above-mentioned tidal current simulation model to be trained are adjusted and the above-mentioned tidal current simulation model is continued to be trained; otherwise, a trained tidal current simulation model is obtained.

[0063] Furthermore, the above-mentioned initial rated power calculation module includes:

[0064] Cross-section acquisition unit, parameter acquisition unit, initial chord length and torsion angle calculation unit, initial rated power calculation unit;

[0065] The cross-section acquisition unit is used to acquire a plurality of cross-sections of the turbine blades;

[0066] The parameter acquisition unit is used to acquire the initial chord length of the first section, the initial torsion angle of the first section, and the preset parameter strategy network;

[0067] The initial chord length and torsion angle calculation unit is used to repeatedly perform the first parameter calculation operation until the initial chord lengths and initial torsion angles of all cross sections are obtained;

[0068] The initial rated power calculation unit is used to calculate the initial rated power according to the initial chord length, initial twist angle and the rated flow velocity of each section.

[0069] Based on the above method embodiment, the present invention provides a terminal device embodiment;

[0070] The present invention provides a terminal device, including a processor, a memory, and a computer program stored in the memory and configured to be executed by the processor. When the processor executes the computer program, it implements a method for determining turbine parameters based on sea area resource characteristics as described in any embodiment of the present invention.

[0071] Based on the above method embodiment, the present invention provides a storage medium embodiment;

[0072] The present invention provides a storage medium, including a processor, a memory, and a computer program stored in the memory and configured to be executed by the processor. When the processor executes the computer program, a method for determining turbine parameters based on sea area resource characteristics according to any embodiment of the present invention is implemented.

[0073] The embodiments of the present invention have the following beneficial effects:

[0074] The present invention provides a method, device, terminal equipment and storage medium for optimizing turbine parameters based on the characteristics of marine resources. The above method first obtains the historical water depth data, historical tide level data and historical sea area coastline data of the undeveloped sea area for the whole year; then obtains the theoretical lowest tide level based on the above water depth data, and uses the above theoretical lowest tide level as the initial value of the sea area water depth; obtains the sea surface safety distance and the seabed safe operating distance, and calculates the turbine diameter based on the difference between the initial value of the sea area water depth, the above sea surface safety distance and the above seabed safe operating distance; then inputs the above water depth data, the above sea area coastline data and the above tide level data into a preset tidal current simulation model, and obtains a number of preset water depths of the undeveloped sea area in the simulated years. The annual average flow rate and the number of hours corresponding to each annual average flow rate; then the preset number of hours is obtained, the above-mentioned number of hours is compared with the preset number of hours, and the target annual average flow rates corresponding to the above-mentioned number of hours is extracted when the above-mentioned number of hours is not less than the preset number of hours, and the smallest annual average flow rate among the above-mentioned target annual average flow rates is used as the rated flow rate of the above-mentioned turbine; then the initial rated power of the turbine is calculated according to the above-mentioned rated flow rate; finally, the turbine parameters are optimized according to the above-mentioned turbine diameter, the above-mentioned rated flow rate and the above-mentioned initial rated power to obtain the optimized turbine parameters; wherein the above-mentioned turbine parameters are the chord length and the twist angle of the turbine blades. Therefore, the present invention obtains historical water depth data, historical tide data and historical sea area coastline data of the sea area to be developed, and then uses the water depth data to determine the turbine diameter, and uses the water depth data and sea area coastline data to determine the rated flow rate of the turbine, and then obtains the initial rated power of the turbine based on the obtained turbine rated flow rate, and finally optimizes the turbine parameters according to the turbine diameter, rated flow rate and initial rated power, thereby solving the problem of determining the turbine rated power according to the actual conditions of the sea area to be planned. BRIEF DESCRIPTION OF THE DRAWINGS

[0075] Figure 1 It is a flow chart of a method for determining turbine parameters based on sea area resource characteristics provided by one embodiment of the present invention.

[0076] Figure 2 It is a schematic diagram of a tidal range cumulative frequency curve provided by an embodiment of the present invention.

[0077] Figure 3 It is a schematic diagram of the structure of a device for determining turbine parameters based on sea area resource characteristics provided by one embodiment of the present invention. DETAILED DESCRIPTION

[0078] The technical solution of the present invention will be described clearly and completely below in conjunction with the accompanying drawings of the present invention. 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 creative work are within the scope of protection of the present invention.

[0079] like Figure 1 As shown, an embodiment of the present invention provides a method for determining turbine parameters based on sea area resource characteristics, comprising:

[0080] Step S101: Acquire historical water depth data, historical tide level data, and historical coastline data of the sea area to be developed throughout the year;

[0081] Preferably, the above-mentioned full year history is at least 1 year.

[0082] In this step, the historical water depth data of the undeveloped sea area throughout the year, the historical tide level data of the whole year, and the historical coastline data of the whole year were obtained.

[0083] Step S102: obtaining a theoretical lowest tide level according to the water depth data, and using the theoretical lowest tide level as an initial value of the sea depth;

[0084] Specifically, after obtaining the above water depth data, the harmonic analysis method is used to obtain the theoretical lowest tide level of the sea area to be developed.

[0085] Specifically, after obtaining the above-mentioned water depth data, the obtained water depth data is first normalized, and then the normalized water depth data is Fourier transformed to obtain multiple simple harmonic tides. The simple harmonic tides include major tides and minor tides. The major tides include: M2 (semi-diurnal tide), S2 (semi-diurnal tide), N2 (semi-diurnal tide), N2 (semi-diurnal tide), K1 (diurnal tide), O1 (diurnal tide) and P1 (diurnal tide); the minor tides include: M4 (quarter diurnal tide), MS4 (quarter diurnal tide) and M6 (one-sixth diurnal tide).

[0086] Specifically, for each of the above simple harmonic tides, the least square method is used to fit it to obtain the amplitude and phase angle of each simple harmonic tide. The amplitude and phase angle obtained here are the harmonic constants. The amplitude represents the size of the simple harmonic tide, and the phase angle represents the time delay of the simple harmonic tide. The long-term tide level is reconstructed using the obtained harmonic constants, and the minimum value of the reconstructed tide level is identified. This minimum value is the lowest water level in the tide theory, that is, the theoretical lowest tide level mentioned above.

[0087] In this step, the above water depth data are analyzed by harmonic analysis method to obtain the theoretical lowest tide level of the sea area to be developed, and then the initial value of the sea depth of the sea area to be developed is obtained.

[0088] Step S103: obtaining the sea surface safety distance and the seabed safety operating distance, and calculating the turbine diameter according to the initial value of the sea area water depth, the sea surface safety distance and the seabed safety operating distance;

[0089] Specifically, the turbine diameter is obtained by subtracting the sea surface safety distance from the initial value of the sea area water depth and then subtracting the seabed safety operating distance.

[0090] In this step, the turbine diameter is calculated based on the initial value of the water depth in the sea area, the sea surface safety distance and the seabed safety operating distance.

[0091] Step S104: inputting the water depth data, the sea area coastline data and the tide level data into a preset tidal current simulation model to obtain the annual average flow velocity at the preset water depth of the sea area to be developed in each simulated year and the number of hours in a year when the annual average flow velocity is reached;

[0092] Preferably, the simulation years simulated by the tidal simulation model are 10 years, the preset water depth is 0.4 times the turbine diameter, and the annual average flow velocity at the preset water depth can be analyzed using the six-point method.

[0093] Specifically, the six-point method is as follows: As the tide level changes at different times with the rise and fall of the tide, the water depth data at each moment is set to 1.0H (bottom water depth), then 0.0H represents the surface water depth, 0.2H represents the 0.2H layer water depth (the depth from the surface to here is 0.2H), 0.4H represents the 0.4H layer water depth (the depth from the surface to here is 0.4H), 0.6H represents the 0.6H layer water depth (the depth from the surface to here is 0.6H), 0.8H represents the 0.8H layer water depth (the depth from the surface to here is 0.8H), so far, the continuously changing flow velocity sequence of each layer is obtained. After obtaining the preset water depth of 0.4 times the turbine diameter, substitute it into the specific water depth to obtain the continuously changing flow velocity sequence value of the water depth, and then the annual average flow velocity value at this depth can be calculated.

[0094] In a preferred embodiment, the construction of the preset power flow simulation model includes:

[0095] Acquire a number of data samples and a number of flow velocity data samples corresponding to the data samples; wherein the data samples include: shoreline data samples, water depth data samples and tide level data samples;

[0096] Inputting the above-mentioned data samples into the tidal current simulation model to be trained, so that the above-mentioned tidal current model to be trained simulates the flow velocity of the simulated year according to the above-mentioned tidal level data samples to obtain the simulated flow velocity corresponding to the simulated year;

[0097] According to the above simulated flow rate and the above several flow rate data samples, a relative error is calculated;

[0098] The relative error is compared with a preset relative error value. If the relative error is smaller than the preset relative error, the parameters of the tidal current simulation model to be trained are adjusted and the tidal current simulation model is continued to be trained; otherwise, a trained tidal current simulation model is obtained.

[0099] Specifically, the above tide data samples are harmonized and analyzed to obtain the harmonic constants of each tide, and then the tidal height and tidal time for the next year are obtained. Then, ADCP (Acoustic Doppler Current Profiler) is used to monitor the velocity and full tide in the sea area to be developed. A complete full tide includes spring tide, moderate tide and neap tide. The tide level method is used to determine the full tide observation date. The tidal height and tidal time are obtained according to the full tide observation date to obtain the daily average tidal range. Then, the tidal range cumulative frequency curve is drawn according to the daily average tidal range. Figure 2 shown.

[0100] Specifically, the forecast tidal range within the planned observation period should meet the following requirements: the cumulative frequency of tidal range on low tide observation dates is greater than 90%, the cumulative frequency of tidal range on mid-tide observation dates is between 40% and 50%, and the cumulative frequency of high tide range is less than 10%.

[0101] Specifically, the tidal current simulation model to be trained is first gridded, the above data samples are input into the tidal current simulation model, and the years to be simulated are set to obtain the flow velocity, flow direction, probability of occurrence of each level of flow velocity and annual average power density of the continuous time series corresponding to each simulated year. The simulated time period is overlapped with the time period corresponding to the flow velocity data sample, the simulated flow velocity curve is drawn, and the relative error between the simulated flow velocity data and the flow velocity data sample is calculated. If the relative error is greater than 85%, it is considered that the model simulation data is accurate and the model training is completed. Otherwise, it is necessary to adjust the grid division method, turbulence simulation form, driving data and other parameters to train the model, and finally obtain an accurate tidal current simulation model.

[0102] In this preferred embodiment, based on the tidal simulation model, the annual average flow velocity of the sea area to be developed at a water depth of 0.4 times the diameter of the turbine in each simulation year and the number of hours corresponding to each annual average flow velocity are obtained.

[0103] Step S105: Compare the number of hours corresponding to each simulated year with the preset number of hours, and extract the annual average flow rate of the simulated year with the number of hours greater than the preset number of hours as the target annual average flow rate;

[0104] Preferably, the preset number of hours is 3000 hours.

[0105] In this step, taking the above content as an example, the target annual average flow rate is determined by extracting the annual average flow rate corresponding to hours greater than 3,000 hours from several annual average flow rates.

[0106] Step S106: taking the smallest annual average flow rate among the above target annual average flow rates as the rated flow rate of the above turbine;

[0107] In this step, the rated flow rate of the turbine is determined.

[0108] Step S107: Calculate the initial rated power of the turbine according to the rated flow rate;

[0109] In a preferred embodiment, the initial rated power of the turbine is calculated based on the rated flow rate, including:

[0110] Obtaining several cross sections of the turbine blades;

[0111] Obtaining the initial chord length of the first section, the initial torsion angle of the first section, and the preset parameter strategy network;

[0112] Repeat the first parameter calculation operation until the initial chord lengths and initial torsion angles of all cross sections are obtained;

[0113] The initial rated power is calculated based on the initial chord length, initial twist angle and the rated flow velocity of each section;

[0114] The above parameter calculation operation includes:

[0115] Obtaining the initial chord length and initial torsion angle of the current section; wherein, when the current section is the first section, the initial chord length and initial torsion angle of the current section are the initial chord length and initial torsion angle of the first section;

[0116] According to the initial chord length of the current section, the initial torsion angle of the current section and the above-mentioned preset parameter strategy network, the initial chord length and initial torsion angle of the next section are calculated;

[0117] Specifically, the above preset parameter strategy network is:

[0118]

[0119] In the formula, θ' represents the initial chord length or initial torsion angle of the next section, θ represents the initial chord length or initial torsion angle of the current section, and lr represents the learning rate, which is a hyperparameter used to control the step size of each update. Represents the gradient of the chord length or torsion angle of the current section.

[0120] In another preferred embodiment, the initial rated power is calculated based on the initial chord length, the initial twist angle and the rated flow velocity of each cross section, including:

[0121] Obtain the distance from each section to the radius of the turbine blade, the turbine speed, and the turbine airfoil;

[0122] Preferably, the above distance is the hub radius of the turbine, which is preset to be 0.06 m.

[0123] Specifically, the turbine speed can be obtained by the preset tip speed ratio and the following formula:

[0124] λ=Ω×r / V;

[0125] Where λ represents the tip speed ratio, Ω represents the turbine speed, r represents the radius of each section, and V represents the rated flow velocity.

[0126] According to the distance, the rotation speed of the turbine and the rated flow rate, the initial inflow angle of each section and the initial relative velocity of each section are calculated;

[0127] Specifically, the initial flow angle of each section is calculated according to the following formula:

[0128]

[0129] In the formula, Φ i represents the incoming flow angle of the i-th section, r i Represents the distance from the ith section to the radius of the turbine blade.

[0130] Specifically, the initial relative velocity of each section is calculated according to the following formula:

[0131]

[0132] In the formula, represents the relative velocity of the i-th section.

[0133] According to the initial twist angle of each section and the corresponding initial incoming flow angle, the initial angle of attack of each section is calculated;

[0134] Specifically, the initial angle of attack of each section is calculated according to the following formula:

[0135] α i =Φ i -β i ;

[0136] In the formula, α i represents the angle of attack of the i-th section, β i represents the torsion angle of the i-th section.

[0137] According to the airfoil and the initial angle of attack of each section, an initial lift coefficient of each section and an initial drag coefficient of each section are obtained;

[0138] Specifically, according to the airfoil and the angle of attack, the lift coefficient and the drag coefficient corresponding to each section can be obtained by looking up the table.

[0139] The initial lift of each section and the initial drag of each section are calculated according to the initial relative velocity of each section, the initial chord length of each section, the initial lift coefficient of each section and the initial drag coefficient of each section;

[0140] Specifically, the initial lift of each section is calculated according to the following formula:

[0141]

[0142] Where, L i represents the lift of the ith section, ρ represents the density of water, in kg / m 3 , c i represents the chord length of the i-th section, C Li represents the lift coefficient of the ith section.

[0143] Specifically, the initial resistance of each section is calculated by the following formula:

[0144]

[0145] Where D i represents the resistance of the ith section, C Di represents the drag coefficient of the ith section.

[0146] The initial moment of each section is calculated according to the initial incoming flow angle of each section, each of the distances, the initial lift of each section and the initial drag of each section;

[0147] Specifically, the initial moment of each section is calculated according to the following formula:

[0148] M i =L i r i cos(Φ i )-D i r i sin(Φ i );

[0149] Where M i represents the moment of the i-th section.

[0150] According to the number of turbine blades and the initial moment of each section, the initial total moment of the turbine is calculated;

[0151] Specifically, the initial total torque is calculated according to the following formula:

[0152]

[0153] Where M total represents the total moment, N represents the number of blades, the number of blades is a preset value, and I represents the number of sections.

[0154] The initial rated power of the turbine is calculated based on the initial total torque and the rotational speed.

[0155] Specifically, the initial rated power is calculated according to the following formula:

[0156] P=ΩM total ;

[0157] Where P represents the rated power in W.

[0158] Step S108: determining the chord length and twist angle of the turbine blades according to the turbine diameter, the rated flow rate and the initial rated power.

[0159] In a preferred embodiment, the chord length and the twist angle of the turbine blades are determined according to the turbine diameter, the rated flow rate and the initial rated power, including:

[0160] Repeating the parameter updating operation according to the turbine diameter, the rated flow rate and the initial rated power until the chord length and the twist angle of the turbine blade are determined;

[0161] The above parameter update operation includes:

[0162] Obtain the current rated power; wherein the initial rated power is the above initial rated power;

[0163] Calculate the current energy efficiency based on the turbine diameter, rated flow rate and current rated power;

[0164] Specifically, the current energy acquisition efficiency is calculated according to the following formula:

[0165]

[0166] In the formula, C p represents the energy acquisition efficiency, A represents the swept area, and the unit is m 2 , H represents the turbine diameter.

[0167] The current instant reward is calculated based on the current chord length, current torsion angle and current energy efficiency of all cross sections. Specifically, the current instant reward is calculated based on the following formula:

[0168] r t =r CP +r β +r c ;

[0169] r Cp_reward =C p *10;

[0170] r Cp_penalty = -10*(C′ p -C p );

[0171] r β_monotony_reward =10;

[0172] r β_range_reward =10;

[0173] r β_monotony_penalty = -10 × T;

[0174] r β_range_penalty = -10 × N;

[0175] r c_range_reward =10;

[0176] r c_smoothness_reward =10;

[0177] r c_range_penalty = -10 × M;

[0178] r c_smoothness_penalty = -10 × k;

[0179]

[0180] In the formula, r t represents the immediate reward, r Cp represents the immediate reward of energy acquisition efficiency, r β represents the instant reward of the twist angle, r c represents the immediate bonus of chord length, C p represents the energy acquisition efficiency, C′ p represents the preset energy efficiency, T represents the number of times the monotonically decreasing condition is not met, N represents the number of torsion angles exceeding the range of [-5°, 45°], M represents the number of chord lengths exceeding the range of [0.05m, 0.5m], and k represents the chord length that does not meet |c i+1 -c i |≤0.1×c i The number of times this condition is i+1represents the chord length of the next section of the i-th section, β i represents the torsion angle of the i-th section, r Cp_reward Reward function representing energy acquisition efficiency, r Cp_penalty represents the penalty function of energy acquisition efficiency, r β_monotony_reward represents the monotonically decreasing reward function of the torsion angle, r β_range_reward represents the range-limited reward function for the torsion angle, r β_monotony_penalty represents the monotonically decreasing penalty function of the torsion angle, r β_range_penalty represents the range limit penalty function of the torsion angle, r c_range_reward represents the range-limited reward function for the chord length, r c_smoothness_reward Represents the reward function for the smoothness of the change in chord length, r c_range_penalty represents the range-limited penalty function for the chord length, r c_smoothness_penalty Represents the penalty function for the smoothness of the change in chord length.

[0181] Preferably, since the change of the chord length at each cross section should be smooth, a smoothing coefficient of the chord length is set, and the coefficient value is 1.1.

[0182] Preferably, the preset energy acquisition efficiency is 0.35.

[0183] According to the current instant reward and the chord lengths and torsion angles of all current sections, the updated chord length and torsion angle of the first section are calculated;

[0184] Specifically, the current instant reward is compared with the previous instant reward. If the current instant reward is greater than the previous instant reward, the updated chord length and torsion angle of the first section are calculated according to the current instant reward and the chord lengths and torsion angles of all current sections, in a preset direction and a preset step size; otherwise, the updated chord length and torsion angle of the first section are calculated according to the current instant reward and the chord lengths and torsion angles of all current sections, in the opposite direction of the preset direction and a preset step size.

[0185] In schematic form, if the current instant reward is 10 and the previous instant reward is 5, the preset direction is the direction of increasing chord length and torsion angle, and the preset step size is 2, then when calculating the updated chord length and torsion angle of the first section, the chord length and torsion angle of all current sections are increased by 2 to obtain the updated chord length and torsion angle of the first section; if the current instant reward is 10 and the previous instant reward is 15, then when calculating the updated chord length and torsion angle of the first section, the chord length and torsion angle of all current sections are reduced by 2 to obtain the updated chord length and torsion angle of the first section.

[0186] Determine whether the updated chord length and torsion angle of the first section converge, if so, recalculate the updated chord lengths and torsion angles of all other sections according to the updated chord length and torsion angle of the first section, and determine the chord length and torsion angle of the turbine blade;

[0187] Otherwise, based on the updated chord length and torsion angle of the first section, the updated chord lengths and torsion angles of all other sections are recalculated to obtain the updated chord lengths and torsion angles of all sections, and the updated rated power is calculated based on the updated chord lengths and torsion angles of all sections.

[0188] Specifically, in step S107, based on the updated chord length and torsion angle of the first section, the updated chord lengths and torsion angles of all other sections are recalculated to obtain the updated chord lengths and torsion angles of all sections, and the updated rated power is calculated based on the updated chord lengths and torsion angles of all sections.

[0189] In this preferred embodiment, the chord length and twist angle of the turbine blades are obtained according to the turbine diameter, rated flow rate and initial rated power.

[0190] Based on the above method embodiment, the present invention provides a corresponding device embodiment.

[0191] like Figure 3 As shown, an embodiment of the present invention provides a device for determining turbine parameters based on sea area resource characteristics, comprising:

[0192] Data acquisition module, water depth calculation module, diameter calculation module, hours determination module, target annual average flow rate determination module, rated flow rate determination module, initial rated power calculation module, and chord length and torsion angle calculation module;

[0193] The data acquisition module is used to obtain the historical water depth data, historical tide level data and historical coastline data of the sea area to be developed throughout the year;

[0194] The water depth calculation module is used to obtain the theoretical lowest tide level according to the water depth data, and use the theoretical lowest tide level as the initial value of the water depth of the sea area;

[0195] The diameter calculation module is used to obtain the sea surface safety distance and the seabed safety operating distance, and calculate the turbine diameter according to the initial value of the sea area water depth, the sea surface safety distance and the seabed safety operating distance;

[0196] The hour determination module is used to input the water depth data, the sea area coastline data and the tide level data into a preset tidal current simulation model to obtain the annual average flow velocity at the preset water depth of the sea area to be developed in each simulated year and the number of hours in a year to reach the annual average flow velocity;

[0197] The above-mentioned target annual average flow rate determination module is used to compare the number of hours corresponding to each simulated year with the preset number of hours, and extract the annual average flow rate of the simulated year with the number of hours greater than the preset number of hours as the target annual average flow rate;

[0198] The rated flow rate determination module is used to use the smallest annual average flow rate among the target annual average flow rates as the rated flow rate of the turbine;

[0199] The initial rated power calculation module is used to calculate the initial rated power of the turbine according to the rated flow rate;

[0200] The chord length and torsion angle calculation module is used to determine the chord length and torsion angle of the turbine blades according to the turbine diameter, the rated flow rate and the initial rated power.

[0201] In a preferred embodiment, the above-mentioned hour determination module includes:

[0202] Data sample acquisition unit, model training unit, relative error calculation unit and model determination unit;

[0203] The sample acquisition unit is used to acquire a number of data samples and a number of flow velocity data samples corresponding to the data samples; wherein the data samples include: shoreline data samples, water depth data samples and tide level data samples;

[0204] The model training unit is used to input the plurality of data samples into the tidal current simulation model to be trained, so that the tidal current model to be trained simulates the flow velocity of the simulation year according to the plurality of tidal level data samples to obtain the simulated flow velocity corresponding to the simulation year;

[0205] The relative error calculation unit is used to calculate the relative error based on the simulated flow rate and the plurality of flow rate data samples;

[0206] The above-mentioned model determination unit is used to compare the above-mentioned relative error with the preset relative error value. If the above-mentioned relative error is smaller than the above-mentioned preset relative error, the parameters of the above-mentioned tidal current simulation model to be trained are adjusted and the above-mentioned tidal current simulation model is continued to be trained; otherwise, a trained tidal current simulation model is obtained.

[0207] In another preferred embodiment, the initial rated power calculation module includes:

[0208] Cross-section acquisition unit, parameter acquisition unit, initial chord length and torsion angle calculation unit, initial rated power calculation unit;

[0209] The cross-section acquisition unit is used to acquire a plurality of cross-sections of the turbine blades;

[0210] The parameter acquisition unit is used to acquire the initial chord length of the first section, the initial torsion angle of the first section, and the preset parameter strategy network;

[0211] The initial chord length and torsion angle calculation unit is used to repeatedly perform the first parameter calculation operation until the initial chord lengths and initial torsion angles of all cross sections are obtained;

[0212] The initial rated power calculation unit is used to calculate the initial rated power according to the initial chord length, initial twist angle and the rated flow velocity of each section.

[0213] It should be noted that the device embodiments described above are merely schematic, wherein the modules described above as separate components may or may not be physically separated, and the components displayed as modules may or may not be physical modules, that is, they may be located in one place, or they may be distributed on multiple network units. Some or all of the modules may be selected according to actual needs to achieve the purpose of the scheme of this embodiment. In addition, in the accompanying drawings of the device embodiments provided by the present invention, the connection relationship between the modules indicates that there is a communication connection between them, which can be specifically implemented as one or more communication buses or signal lines. Ordinary technicians in this field can understand and implement it without paying creative labor. The above schematic diagram is only an example of a device for determining turbine parameters based on marine resource characteristics, and does not constitute a limitation on a device for determining turbine parameters based on marine resource characteristics, and may include more or fewer components than shown in the figure, or combine certain components, or different components.

[0214] Based on the above method item embodiments, the present invention provides corresponding terminal device item embodiments.

[0215] Another embodiment of the present invention provides a terminal device, including a processor, a memory, and a computer program stored in the memory and configured to be executed by the processor. When the processor executes the computer program, it implements a method for determining turbine parameters based on sea area resource characteristics as described in any embodiment of the present invention.

[0216] Exemplarily, in this embodiment, the computer program may be divided into one or more modules, which are stored in the memory and executed by the processor to implement the present invention. The one or more modules may be a series of computer program instruction segments capable of implementing specific functions, which are used to describe the execution process of the computer program in the device;

[0217] The terminal device may be a computing device such as a desktop computer, a notebook, a PDA, a cloud server, etc. The device may include, but is not limited to, a processor, a memory;

[0218] The processor referred to herein may be a central processing unit (CPU), other general-purpose processors, digital signal processors (DSP), application-specific integrated circuits (ASIC), field-programmable gate arrays (FPGA) or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc. The general-purpose processor may be a microprocessor or any conventional processor, etc. The above-mentioned processor is the control center of the above-mentioned device, and uses various interfaces and lines to connect various parts of the entire device;

[0219] The above-mentioned memory can be used to store the above-mentioned computer programs and / or modules. The above-mentioned processor realizes various functions of the above-mentioned device by running or executing the computer programs and / or modules stored in the above-mentioned memory, and calling the data stored in the memory. The above-mentioned memory can mainly include a program storage area and a data storage area, wherein the program storage area can store an operating system, at least one application required for a function, etc.; in addition, the memory can include a high-speed random access memory, and can also include a non-volatile memory, such as a hard disk, a memory, a plug-in hard disk, a smart memory card (Smart Media Card, SMC), a secure digital (Secure Digital, SD) card, a flash card (Flash Card), at least one disk storage device, a flash memory device, or other volatile solid-state storage devices.

[0220] Based on the above method item embodiments, the present invention provides a corresponding storage medium item embodiment.

[0221] Another embodiment of the present invention provides a storage medium, which includes a stored computer program, wherein when the computer program is running, the device where the storage medium is located is controlled to execute a method for determining turbine parameters based on sea area resource characteristics as described in any embodiment of the present invention.

[0222] In this embodiment, the storage medium is a computer-readable storage medium, the computer program includes computer program code, and the computer program code may be in source code form, object code form, executable file or some intermediate form, etc. The computer-readable medium may include: any entity or device capable of carrying the computer program code, recording medium, USB flash drive, mobile hard disk, magnetic disk, optical disk, computer memory, read-only memory (ROM), random access memory (RAM), electric carrier signal, telecommunication signal and software distribution medium, etc.

[0223] Compared with the prior art, by implementing the above-mentioned embodiments of the present invention, the rated power of the turbine can be determined according to the local sea conditions, and then the specific turbine parameters can be obtained.

[0224] The above is a preferred embodiment of the present invention. It should be pointed out that for ordinary technicians in this technical field, several improvements and modifications can be made without departing from the principle of the present invention. These improvements and modifications are also considered to be within the scope of protection of the present invention.

Claims

1. A method for determining turbine parameters based on marine resource characteristics, characterized in that: include: Obtain the historical full-year water depth data, historical full-year tide level data, and historical full-year sea area coastline data of the sea area to be developed; Obtaining a theoretical minimum tide level according to the water depth data, and using the theoretical minimum tide level as an initial value of the sea area water depth; Obtaining a sea surface safety distance and a seabed safety operating distance, and calculating a turbine diameter according to the initial value of the sea area water depth, the sea surface safety distance and the seabed safety operating distance; Input the water depth data, the sea area coastline data and the tide level data into a preset tidal current simulation model to obtain the annual average flow velocity at the preset water depth of the sea area to be developed in each simulated year and the number of hours in a year when the annual average flow velocity is reached; The number of hours corresponding to each simulated year is compared with the preset number of hours, and the annual average flow rate of the simulated year with the number of hours greater than the preset number of hours is extracted as the target annual average flow rate; Taking the smallest annual average flow rate among the target annual average flow rates as the rated flow rate of the turbine; Calculating the initial rated power of the turbine according to the rated flow rate; The chord length and the twist angle of the turbine blades are determined according to the turbine diameter, the rated flow rate and the initial rated power.

2. A method for determining turbine parameters based on marine resource characteristics according to claim 1, characterized in that: The construction of the preset power flow simulation model includes: Acquire a number of data samples and a number of flow velocity data samples corresponding to the data samples; wherein the data samples include: shoreline data samples, water depth data samples and tide level data samples; Inputting the plurality of data samples into a tidal current simulation model to be trained, so that the tidal current model to be trained simulates the flow velocity of the simulated year according to the plurality of tidal level data samples to obtain a simulated flow velocity corresponding to the simulated year; Calculating a relative error based on the simulated flow rate and the plurality of flow rate data samples; The relative error is compared with a preset relative error value. If the relative error is smaller than the preset relative error, the parameters of the power flow simulation model to be trained are adjusted and the power flow simulation model is continued to be trained; otherwise, a trained power flow simulation model is obtained.

3. A method for determining turbine parameters based on marine resource characteristics according to claim 2, characterized in that: The step of calculating the initial rated power of the turbine according to the rated flow rate comprises: Acquire several cross sections of the turbine blade; Obtaining the initial chord length of the first section, the initial torsion angle of the first section, and the preset parameter strategy network; Repeat the first parameter calculation operation until the initial chord lengths and initial torsion angles of all cross sections are obtained; The initial rated power is calculated according to the initial chord length, the initial twist angle and the rated flow velocity of each section; The parameter calculation operation includes: Obtaining an initial chord length and an initial torsion angle of the current section; wherein, when the current section is the first section, the initial chord length and the initial torsion angle of the current section are the initial chord length and the initial torsion angle of the first section; According to the initial chord length of the current section, the initial torsion angle of the current section and the preset parameter strategy network, the initial chord length and the initial torsion angle of the next section are calculated.

4. A method for determining turbine parameters based on marine resource characteristics according to claim 3, characterized in that: The initial rated power is calculated according to the initial chord length, the initial torsion angle and the rated flow velocity of each section, including: Obtaining the distance from each section to the radius of the turbine blade, the turbine speed, and the turbine airfoil; According to the distance, the rotation speed of the turbine and the rated flow rate, an initial inflow angle of each section and an initial relative velocity of each section are calculated; Calculate the initial angle of attack of each section according to the initial twist angle of each section and the corresponding initial incoming flow angle; According to the airfoil and the initial angle of attack of each section, an initial lift coefficient of each section and an initial drag coefficient of each section are obtained; Calculating the initial lift of each section and the initial drag of each section according to the initial relative velocity of each section, the initial chord length of each section, the initial lift coefficient of each section and the initial drag coefficient of each section; Calculating the initial moment of each cross section according to the initial incoming flow angle of each cross section, each of the distances, the initial lift of each cross section, and the initial drag of each cross section; Calculating the initial total moment of the turbine according to the number of turbine blades and the initial moment of each section; The initial rated power of the turbine is calculated according to the initial total torque and the rotational speed.

5. A method for determining turbine parameters based on marine resource characteristics according to claim 4, characterized in that: Determining the chord length and the twist angle of the turbine blades according to the turbine diameter, the rated flow rate and the initial rated power includes: Repeating the parameter updating operation according to the turbine diameter, the rated flow rate and the initial rated power until the chord length and the twist angle of the turbine blade are determined; The parameter updating operation includes: Obtaining the current rated power; wherein the initial rated power is the initial rated power; Calculate the current energy efficiency based on the turbine diameter, rated flow rate and current rated power; Calculate the current instant reward based on the current chord length, current torsion angle and current energy efficiency of all cross sections; According to the current instant reward and the chord lengths and torsion angles of all current sections, the updated chord length and torsion angle of the first section are calculated; Determine whether the updated chord length and torsion angle of the first section have converged. If so, recalculate the updated chord lengths and torsion angles of all other sections based on the updated chord length and torsion angle of the first section to determine the chord length and torsion angle of the turbine blade. Otherwise, based on the updated chord length and torsion angle of the first section, the updated chord lengths and torsion angles of all other sections are recalculated to obtain the updated chord lengths and torsion angles of all sections, and the updated rated power is calculated based on the updated chord lengths and torsion angles of all sections.

6. A device for determining turbine parameters based on marine resource characteristics, characterized in that: include: Data acquisition module, water depth calculation module, diameter calculation module, hours determination module, target annual average flow rate determination module, rated flow rate determination module, initial rated power calculation module, and chord length and torsion angle calculation module; The data acquisition module is used to obtain the historical water depth data of the sea area to be developed throughout the year, the historical tide level data of the year, and the historical sea area coastline data of the year; The water depth calculation module is used to obtain a theoretical minimum tide level according to the water depth data, and use the theoretical minimum tide level as an initial value of the sea area water depth; The diameter calculation module is used to obtain the sea surface safety distance and the seabed safety operating distance, and calculate the turbine diameter according to the initial value of the sea area water depth, the sea surface safety distance and the seabed safety operating distance; The hour determination module is used to input the water depth data, the sea area coastline data and the tide level data into a preset tidal current simulation model to obtain the annual average flow velocity at the preset water depth of the sea area to be developed in each simulated year and the number of hours in a year to reach the annual average flow velocity; The target annual average flow rate determination module is used to compare the number of hours corresponding to each simulated year with the preset number of hours, and extract the annual average flow rate of the simulated year with the number of hours greater than the preset number of hours as the target annual average flow rate; The rated flow rate determination module is used to take the smallest annual average flow rate among the target annual average flow rates as the rated flow rate of the turbine; The initial rated power calculation module is used to calculate the initial rated power of the turbine according to the rated flow rate; The chord length and torsion angle calculation module is used to determine the chord length and torsion angle of the turbine blades according to the turbine diameter, the rated flow rate and the initial rated power.

7. A device for determining turbine parameters based on marine resource characteristics according to claim 6, characterized in that: The hours determination module comprises: Data sample acquisition unit, model training unit, relative error calculation unit and model determination unit; The sample acquisition unit is used to acquire a number of data samples and a number of flow rate data samples corresponding to the data samples; wherein the data samples include: shoreline data samples, water depth data samples and tide level data samples; The model training unit is used to input the plurality of data samples into the tidal current simulation model to be trained, so that the tidal current model to be trained simulates the flow velocity of the simulation year according to the plurality of tidal data samples to obtain the simulated flow velocity corresponding to the simulation year; The relative error calculation unit is used to calculate the relative error according to the simulated flow rate and the plurality of flow rate data samples; The model determination unit is used to compare the relative error with a preset relative error value. If the relative error is smaller than the preset relative error, the parameters of the tidal current simulation model to be trained are adjusted and the tidal current simulation model is continued to be trained; otherwise, a trained tidal current simulation model is obtained.

8. The device for determining turbine parameters based on sea area resource characteristics according to claim 7, characterized in that: The initial rated power calculation module comprises: Cross-section acquisition unit, parameter acquisition unit, initial chord length and torsion angle calculation unit, initial rated power calculation unit; The cross-section acquisition unit is used to acquire a plurality of cross-sections of the turbine blade; The parameter acquisition unit is used to acquire an initial chord length of the first section, an initial torsion angle of the first section, and a preset parameter strategy network; The initial chord length and torsion angle calculation unit is used to repeatedly perform the first parameter calculation operation until the initial chord lengths and initial torsion angles of all cross sections are obtained; The initial rated power calculation unit is used to calculate the initial rated power according to the initial chord length, the initial torsion angle and the rated flow velocity of each cross section.

9. A terminal device, characterized in that: It comprises a processor, a memory, and a computer program stored in the memory and configured to be executed by the processor, and when the processor executes the computer program, it implements a method for determining turbine parameters based on sea area resource characteristics as described in any one of claims 1 to 5.

10. A storage medium, characterized in that: The storage medium includes a stored computer program, wherein when the computer program is running, the device where the storage medium is located is controlled to execute a method for determining turbine parameters based on sea area resource characteristics as described in any one of claims 1 to 5.

Images