METHOD, DEVICE AND STORAGE MEDIUM FOR EVALUATING WIND ENERGY RESOURCES IN COMPLEX TERRAIN

Abstract

A method, a device and a storage medium for evaluating wind energy resources in complex terrain are provided, and the method includes: obtaining a climate field based on observation data of wind speed; obtaining an anomaly field; superimposing a climate field interpolation result and an outlier interpolation result with a consistent spatial resolution to obtain a wind speed interpolation result; performing a deviation correction on the wind speed interpolation result and the observation data of wind speed to obtain a final result; calculating an average effective wind power density; and estimating a wind power density based on a daily average wind speed. An accuracy of wind speed data is improved; the wind energy resources are evaluated in situations including complex terrain and lack of hourly wind speed data; and a high-precision data set of the wind energy resources is established to improve an evaluation accuracy of the wind energy resources.

Claims

1. A method for evaluating wind energy resources in complex terrain, comprising: step 1, obtaining a climate field based on observation data of wind speed; wherein the obtaining a climate field based on observation data of wind speed comprises: obtaining an average climate field based on the observation data of wind speed, and performing a spatial interpolation on the average climate field by using a thin-plate smoothing spline function of terrain covariates to obtain a climate field interpolation result, wherein an interpolation accuracy of the average climate field is consistent with an accuracy required for evaluating the wind energy resources; step 2, obtaining an anomaly field; wherein the obtaining an anomaly field comprises: obtaining a difference between each observation data of wind speed and the climate field interpolation result as an outlier, and performing a spatial interpolation on the outlier by using a thin- plate smoothing spline function of terrain covariates to obtain an outlier interpolation result, wherein an interpolation accuracy of the outlier is consistent with the accuracy required for evaluating the wind energy resources; step 3, superimposing the climate field interpolation result and the outlier interpolation result with a consistent spatial resolution to obtain a wind speed interpolation result; step 4, performing a deviation correction on the wind speed interpolation result and the observation data of wind speed to obtain a final result; wherein the deviation correction comprises: an equidistant cumulative distribution function method; and original observation data of wind speed is processed through steps 1-3 when target-precision observation data of wind speed is lacked; step 5, calculating an hourly average effective wind power density; wherein step 5 comprises: step 5.1, estimating an hourly wind power density based on an hourly wind speed, wherein a formula of the hourly wind power density is expressed as follows: $\mathrm{WP}=\frac{1}{2n}{?}_{\mathrm{air}}{v}^3;$ wherein WP represents the hourly wind power density, v represents the hourly wind speed, and ?.sub.air represents an air density; wherein a calculation formula of the air density is expressed as follows: ${?}_{\mathrm{air}}=\frac{P_{\mathrm{ave}}}{R{T}_{\mathrm{ave}}};$ wherein P.sub.ave represents annual average atmospheric pressure, R represents a gas constant, and T.sub.ave represents an annual average temperature; step 5.2, calculating an hourly average wind power density, wherein a formula of the hourly average wind power density is expressed as follows: $\stackrel{\_}{\mathrm{WP}}=\frac{1}{2n}\underset{i=1}{\overset{n}{.Math.}}{?}_{\mathrm{air}}{v}_i^3;$ wherein WP represents the hourly average wind power density, n represents a number of records in a set period, v.sub.i represents an hourly wind speed of an i-th record of the n records, and ?.sub.air represents the air density; step 5.3, calculating an hourly effective wind power density, wherein a formula of the hourly effective wind power density is expressed as follows: ${\mathrm{WP}}_E=\frac{1}{2}{?}_{\mathrm{air}}{?}_{v_{\mathrm{start}}}^{{?}_{\mathrm{stop}}}{\mathrm{Function}}_P\left(v\right).Math.\mathrm{vdv};$ wherein WP.sub.E represents the hourly effective wind power density, v.sub.start represents a start-up wind speed, v.sub.stop represents a shutdown wind speed, ?.sub.air represents the air density, and Function.sub.P(v) represents a probability density function of wind speed; and step 5.4, applying the formula of the hourly effective wind power density to calculate the hourly average wind power density to thereby obtain the hourly average effective wind power density, wherein the applying the formula of the hourly effective wind power density to calculate the hourly average wind power density to thereby obtain the hourly average effective wind power density comprises: assuming a number of records of an hourly effective wind speed within the n records in the set period as m, wherein the formula of the hourly average wind power density is expressed as follows: $\stackrel{\_}{\mathrm{WP}}=\frac{1}{2n}\underset{i=1}{\overset{n}{.Math.}}{?}_{\mathrm{air}}{v}_i^3=\frac{1}{2n}.Math.{?}_{\mathrm{air}}.Math.\left(v_1^3+v_2^3+v_3^3+.Math.+v_n^3\right)=\frac{1}{2n}.Math.{?}_{\mathrm{air}}.Math.\left(v_1^3+v_2^3+v_3^3+.Math.+v_m^3+v_{m+1}^3+.Math.+v_n^3\right);$ since a wind power density that does not belong to the hourly effective wind speed is zero in the calculation of the hourly effective wind power density, obtaining a formula as follows: $v_{m+1}^3+.Math.+v_n^3=0;$ wherein a formula of the hourly average effective wind power density is expressed as follows: $\stackrel{\_}{{\mathrm{WP}}_E}=\frac{1}{2n}.Math.{?}_{\mathrm{air}}.Math.\left(v_1^3+v_2^3+v_3^3+.Math.+v_m^3\right)=\frac{1}{2n}\underset{i=1}{\overset{m}{.Math.}}{?}_{\mathrm{air}}{v}_i^3;$ wherein WP.sub.E represents the hourly average effective wind power density, n represents the number of records in the set period, m represents the number of records of the hourly effective wind speed in the set period, and ?.sub.air represents the air density; and step 6, estimating a daily average effective wind power density based on a daily average wind speed; wherein step 6 comprises: assuming a wind speed of a i-th record being ?.sub.i times of the daily average wind speed v, wherein a formula of the daily average effective wind power density is expressed as follows: $\stackrel{\_}{{\mathrm{WP}}_E^{}}=\frac{1}{2n}.Math.{?}_{\mathrm{air}}.Math.\left[{\left({?}_1\stackrel{\_}{v}\right)}^3+{\left({?}_2\stackrel{\_}{v}\right)}^3+{\left({?}_3\stackrel{\_}{v}\right)}^3+.Math.+{\left({?}_m\stackrel{\_}{v}\right)}^3\right]=\frac{1}{2n}.Math.{?}_{\mathrm{air}}.Math.{\stackrel{\_}{v}}^3\left({?}_1^3+{?}_2^3+{?}_3^3+.Math.+{?}_m^3\right)=\frac{1}{2n}{?}_{\mathrm{air}}{\stackrel{\_}{v}}^3\underset{i=1}{\overset{m}{.Math.}}\left({?}_i^3\right);$ making $?=\underset{i=1}{\overset{m}{.Math.}}\left({?}_t^3\right),$ wherein the formula of the daily average effective wind power density is expressed as follows: $\stackrel{\_}{{\mathrm{WP}}_E^{}}=\frac{1}{2n}{?}_{\mathrm{air}}?{\stackrel{\_}{v}}^3;$ wherein WP.sub.E represents the daily average effective wind power density, n represents a number of records in a set period, v represents the daily average wind speed, ? represents a ratio set of the hourly wind speed and the daily average wind speed, the probability density function of wind speed Function.sub.P(v), a shape parameter k and a scaling parameter c are calculated by a Weibull distribution, and formulas of the probability density function of wind speed Function.sub.P(v), the shape parameter k and the scaling parameter c are expressed as follows: ${\mathrm{Function}}_p\left(v\right)=\frac{{\mathrm{Function}}_P\left(v\right)}{{\mathrm{Function}}_P\left(v_{\mathrm{start}}?v?v_{\mathrm{stop}}\right)}=\frac{\left(\frac{k}{c}\right){\left(\frac{v}{c}\right)}^{k-1}\exp \left[-{\left(\frac{v}{c}\right)}^k\right]}{\exp \left[-{\left(\frac{v_{\mathrm{start}}}{c}\right)}^k\right]-\exp \left[-{\left(\frac{v_{\mathrm{stop}}}{c}\right)}^k\right]};$ $k={\left[\frac{\mathrm{ave}\left(v\right)}{\mathrm{Std}.\mathrm{Deviation}\left(v\right)}\right]}^{1.086};$ $c=\frac{\mathrm{ave}\left(v\right)}{?\left(1+\frac{1}{k}\right)};$ wherein Function.sub.P(v) represents the probability density function of wind speed, k represents the shape parameter, c represents the scaling parameter, Std. Deviation(v) represents standard deviation of the wind speed, and ? represents a gamma function; and wherein the method for evaluating wind energy resources in complex terrain further comprises: estimating, based on the daily average effective wind power density, a total amount of exploitable wind energy resources of a target area, determining whether the target area is suitable for developing the wind energy resources, when the target area is determined to be suitable for developing the wind energy resources, developing wind energy resources of the target area by working personnel.

2. The method as claimed in claim 1, wherein the calculating an average climate field comprises: selecting data for calculating the climate filed from an average value of wind speed observation period for thirty years.

3. The method as claimed in claim 1, wherein the equidistant cumulative distribution function method comprises formulas expressed as follows: $?={\mathrm{Function}}_{\mathrm{obs}}^{-1}\left({\mathrm{Function}}_{\mathrm{output}2}\left(\mathrm{variable}\right)\right)-{\mathrm{Function}}_{\mathrm{output}1}^{-1}\left({\mathrm{Function}}_{\mathrm{output}2}\left(\mathrm{variable}\right)\right);$ ${\mathrm{variable}}_{\mathrm{correct}}=\mathrm{variable}+?;$ wherein variable represents input data of a climate element variable, variable.sub.correct represents a correction result of the climate element variable, Fuction represents an equidistant cumulative distribution function, Fuction.sup.?1 represents an inverse operation of the equidistant cumulative distribution function, obs represents observation data of wind speed during training, output1 represents an output result during training, and output2 represents an output result during correcting.

4-5. (canceled)

6. A device for evaluating wind energy resources in complex terrain, comprising a memory, a processor and a computer program stored in the memory and executed in the processor, wherein the computer program is configured to be executed by the processor to implement the steps of the method as claimed in claim 1.

7. (canceled)

Description

BRIEF DESCRIPTION OF DRAWING

[0041] FIGURE illustrates a functional block diagram of a method for evaluating wind energy resources in complex terrain according to an embodiment of the disclosure.

DETAILED DESCRIPTION OF EMBODIMENTS

[0042] Technique solutions of the disclosure will be further described in conjunction with drawings below.

[0043] An embodiment of the disclosure provides a method for evaluating wind energy resources in complex terrain, and the method includes the following steps 1-6. [0044] In step 1, a climate field is obtained based on observation data of wind speed. Firstly, an average climate field is obtained based on the observation data of wind speed, then a spatial interpolation is performed on the average climate filed by using a thin-plate smoothing spline function of terrain covariates (i.e., adding the value of the average climate filed into the thin-plate smoothing spline function for calculation) to obtain a climate field interpolation result, and an interpolation accuracy of the average climate filed is consistent with an accuracy required for evaluating the wind energy resources. A step for obtaining the average climate field includes that data for calculating the climate field is selected from an average value of wind speed observation period for thirty years. [0045] In step 2, an anomaly field is obtained. Firstly, a difference between each observation data of wind speed and the climate field interpolation result is obtained as an outlier, then a spatial interpolation is performed on the outlier by using a thin-plate smoothing spline function of terrain covariates (i.e., adding the outlier into the thin-plate smoothing spline function for calculation) to obtain an outlier interpolation result (i.e., anomaly field interpolation result), and an interpolation accuracy of the outlier is consistent with the accuracy required for evaluating wind energy resource. [0046] In step 3, the climate field interpolation result (i.e., climate field) and the outlier interpolation result (i.e., anomaly field) with a consistent spatial resolution are superimposed to obtain a wind speed interpolation result. [0047] In step 4, a deviation correction is performed on the wind speed interpolation result and the observation data of wind speed to obtain a final result; and the deviation correction includes an equidistant cumulative distribution function method; and original observation data of wind speed is processed through steps 1-3 when target-precision observation data of wind speed is lacked. Formulas of the equidistant cumulative distribution function method are expressed as follows:

[00012]$?={\mathrm{Function}}_{\mathrm{obs}}^{-1}\left({\mathrm{Function}}_{\mathrm{output}2}\left(\mathrm{variable}\right)\right)-{\mathrm{Function}}_{\mathrm{output}1}^{-1}\left({\mathrm{Function}}_{\mathrm{output}2}\left(\mathrm{variable}\right)\right);$${\mathrm{variable}}_{\mathrm{correct}}=\mathrm{variable}+?;$ [0048] where variable represents input data of a climate element variable, variable.sub.correct represents a correction result of the climate element variable, Fuction represents an equidistant cumulative distribution function, Fuction.sub.?1 represents an inverse operation of the equidistant cumulative distribution function, obs represents observation data of wind speed during training, output1 represents an output result during training, and output2 represents an output result during correcting. [0049] In step 5, an average effective wind power density (i.e., hourly average effective wind power density) is calculated, and the step 5 includes the following steps 5.1-5.4. [0050] In step 5.1, a wind power density is estimated based on an hourly wind speed, and a formula of the wind power density is expressed as follows:

[00013]$WP=\frac{1}{2}{?}_{\mathrm{air}}{v}^3;$ [0051] where WP represents the wind power density, v represents the hourly wind speed, and ?.sub.air represents an air density.

[0052] Furthermore, in step 5.1, a formula of the air density is expressed as follows:

[00014]${?}_{\mathrm{air}}=\frac{P_{\mathrm{ave}}}{{\mathrm{RT}}_{\mathrm{ave}}};$ [0053] where P.sub.ave represents annual average atmospheric pressure, R represents a gas constant, and T.sub.ave represents an annual average temperature;

[0054] In step 5.2, an average wind power density is calculated, and a formula of the average wind power density is expressed as follows:

[00015]$\stackrel{\_}{WP}=\frac{1}{2n}{.Math.}_{i=1}^n{?}_{\mathrm{air}}{v}_i^3;$ [0055] where WP represents the average wind power density, n represents a number of records in a set period, v.sub.i represents a wind speed of an i-th record of the n records, and ?.sub.air represents the air density.

[0056] In step 5.3, an effective wind power density is calculated, and a formula of the effective wind power density is expressed as follows:

[00016]$W{P}_E=\frac{1}{2}{?}_{\mathrm{air}}{?}_{v_{\mathrm{start}}}^{v_{\mathrm{stop}}}{\mathrm{Function}}_P\left(v\right).Math.\mathrm{vdv};$ [0057] where WP.sub.E represents the effective wind power density, v.sub.start represents a start-up wind speed, v.sub.stop represents a shutdown wind speed, ?.sub.air represents the air density, and Function.sub.P(v) represents a probability density function of wind speed.

[0058] In step 5.4, the formula of the effective wind power density is applied to calculate the average wind power density to thereby obtain the average effective wind power density, and a process for the applying includes the following steps.

[0059] A number of records of an effective wind speed within the n records in the set period is assumed as m, and the formula of the average wind power density is expressed as follows:

[00017]$\stackrel{\_}{WP}=\frac{1}{2n}\underset{i=1}{\overset{n}{.Math.}}{?}_{\mathrm{air}}{v}_i^3=\frac{1}{2n}.Math.{?}_{air}.Math.\left(v_1^3+v_2^3+v_3^3+.Math.+v_n^3\right)=\frac{1}{2n}.Math.{?}_{\mathrm{air}}.Math.\left(v_1^3+v_2^3+v_3^3+.Math.+v_m^3+v_{m+1}^3+.Math.+v_n^3\right);$

[0060] Since a wind power density that does not belong to the effective wind speed is zero in the calculation of the effective wind power density, a formula is obtained as follows:

[00018]$v_{m+1}^3+.Math.+v_n^3=0.$

[0061] Furthermore, in step 5.4, a formula of the average effective wind power density is expressed as follows:

[00019]$\stackrel{\_}{W{P}_E}=\frac{1}{2n}.Math.{?}_{\mathrm{air}}.Math.\left(v_1^3+v_2^3+v_3^3+.Math.+v_m^3\right)=\frac{1}{2n}{.Math.}_{i=1}^m{?}_{\mathrm{air}}{v}_i^3;$ [0062] where WP.sub.E represents the average effective wind power density, n represents the number of records in the set period, m represents the number of records of the effective wind speed in the set period, and ?.sub.air represents the air density.

[0063] In step 6, a wind power density (i.e., daily average effective wind power density) is estimated based on a daily average wind speed, and the step 6 includes the following steps.

[0064] A wind speed of a i-th record is assumed as ?.sub.i times of an average wind speed v (i.e., daily average wind speed), and a formula of the average effective wind power density (i.e., daily average effective wind power density) is expressed as follows:

[00020]$\stackrel{\_}{W{P}_E^{}}=\frac{1}{2n}.Math.{?}_{\mathrm{air}}.Math.\left[{\left({?}_1\stackrel{?}{v}\right)}^3+{\left({?}_2\stackrel{?}{v}\right)}^3+{\left({?}_3\stackrel{?}{v}\right)}^3+.Math.+{\left({?}_m\stackrel{?}{v}\right)}^3\right]=\frac{1}{2n}.Math.{?}_{\mathrm{air}}.Math.{\stackrel{?}{v}}^3\left({?}_1^3+{?}_2^3+{?}_3^3+.Math.+{?}_m^3\right)=\frac{1}{2n}{?}_{\mathrm{air}}{\stackrel{?}{v}}^3{.Math.}_{i=1}^m\left({?}_i^3\right);$

[0065] When a formula ?=?.sub.i=1.sup.m(?.sub.i.sup.3) is satisfied, the formula of the average effective wind power density is expressed as follows:

[00021]$\stackrel{\_}{{\mathrm{WP}}_E^{}}=\frac{1}{2n}{?}_{\mathrm{air}}?{\stackrel{?}{v}}^3;$ [0066] where WP.sub.E represents the daily average effective wind power density, n represents the number of records in a set period, v represents the average wind speed, i.e., the daily average wind speed and v=ave(v); ? represents a ratio set of the hourly wind speed and the average wind speed, a probability density function of wind speed Function.sub.P(v), a shape parameter k and a scaling parameter c are calculated by a Weibull distribution, and formulas of the probability density function of wind speed Function.sub.P(v), the shape parameter k and the scaling parameter c are expressed as follows:

[00022]${\mathrm{Function}}_P\left(v\right)=\frac{{\mathrm{Function}}_P\left(v\right)}{{\mathrm{Function}}_P\left(v_{\mathrm{start}}?v?v_{\mathrm{stop}}\right)}=\frac{\left(\frac{k}{c}\right){\left(\frac{v}{c}\right)}^{k-1}\exp \left[-{\left(\frac{v}{c}\right)}^k\right]}{\exp \left[-{\left(\frac{v_{\mathrm{start}}}{c}\right)}^k\right]-\exp \left[-{\left(\frac{v_{\mathrm{stop}}}{c}\right)}^k\right]};$$k={\left[\frac{\mathrm{ave}\left(v\right)}{\mathrm{Std}.\mathrm{Deviation}\left(v\right)}\right]}^{1.086};$$c=\frac{\mathrm{ave}\left(v\right)}{?\left(1+\frac{1}{k}\right)};$ [0067] where Function.sub.P(v) represents the probability density function of wind speed, k represents the shape parameter, c represents the scaling parameter, Std. Deviation(v) represents standard deviation of the wind speed, and ? represents a gamma function.

[0068] Specifically, a calculation result of the probability density function of wind speed Function.sub.P(v) is the corresponding ratio ?.sub.i of the hourly wind speed and the daily average wind speed. After the ratio ?.sub.i is obtained, ? is obtained according to ?=?.sub.i=1.sup.m(?.sub.i.sup.3), thus the daily average effective wind power density WP.sub.E is obtained according to

[00023]$\stackrel{\_}{{\mathrm{WP}}_E^{}}=\frac{1}{2n}{?}_{\mathrm{air}}?{\stackrel{?}{v}}^3.$

[0069] An embodiment of the disclosure further provides a device for evaluating wind energy resources in complex terrain, and the device includes a memory, a processor and a computer program stored in the memory and executed in the processor, and the computer program is configured to be executed by the processor to implement the steps of the method for evaluating wind energy resources in complex terrain.

[0070] An embodiment of the disclosure further provides a storage medium for evaluating wind energy resources in complex terrain, and the storage medium stores a computer program therein, and the computer program is configured to be executed to implement the steps of the method for evaluating wind energy resources in complex terrain. In some embodiments, the storage medium is a non-transitory storage medium.