CN111382488A

CN111382488A - Method and device for eliminating phase deviation of time sequence curve

Info

Publication number: CN111382488A
Application number: CN201811621142.4A
Authority: CN
Inventors: 邓建军; 何东亚; 李岩
Original assignee: Beijing Goldwind Science and Creation Windpower Equipment Co Ltd
Current assignee: Beijing Goldwind Science and Creation Windpower Equipment Co Ltd
Priority date: 2018-12-28
Filing date: 2018-12-28
Publication date: 2020-07-07

Abstract

The invention provides a method and a device for eliminating phase deviation of a time sequence curve, wherein the eliminating method comprises the following steps: determining a plurality of first extreme points on a first timing curve; determining a plurality of second extreme points on a second timing curve; determining a plurality of control point pairs based on the plurality of first extreme points and the plurality of second extreme points; based on each control point pair, the second timing curve is shifted by taking the first timing curve as a reference so as to enable the phases of the first timing curve and the second timing curve to be consistent. By adopting the method and the device for eliminating the phase deviation of the time sequence curve, the aim of eliminating the global phase deviation can be fulfilled.

Description

Method and device for eliminating phase deviation of time sequence curve

Technical Field

The present invention relates generally to the field of signal processing technology, and more particularly, to a method and apparatus for eliminating phase deviation of a timing curve.

Background

At present, in the wind turbine adaptability analysis and tower customization business, different wind speeds, air densities, turbulence degrees, inflow angles, wind shears, initial phase angles, different wind conditions (e.g., NTM, ECD, EOG, EWM, ETM, EWS, etc.) and different working states of the wind turbine (e.g., the conditions of normal power generation, and mutual coupling of faults) need to be considered for the load calculation of the whole wind turbine, generally, the number of the calculated working conditions is up to 874, even thousands, these working conditions can be generally organized into a working condition table, automatically analyzed by a program, then submitted to a cloud platform and automatically calculated by Bladed software, and finally, the ultimate load and the fatigue load under all the working conditions are solved, so that the efficient and safe operation of the wind turbine is ensured.

It can be seen that the load computing method has a huge computing amount, and under the background that the current cloud computing resources are expensive, the service needs a large expenditure, and at this time, it is very urgent to reduce the cost.

At present, various methods such as theoretical calculation, machine learning, interpolation calculation and the like can be adopted for load prediction of key components of the wind turbine generator, different methods need different calculated amounts, and different calculation accuracies can be obtained. The interpolation method has high calculation efficiency and relatively reliable precision, so that the interpolation method becomes a novel method which is innovative and challenging in the aspect of wind turbine load prediction.

The interpolation method needs to perform interpolation processing based on two load time sequence curves under different values of a single wind parameter, but the load of the variable-speed variable-pitch wind turbine generator in turbulent wind has obvious non-stationarity and inevitably generates phase deviation (as shown in fig. 1 and 2) under the single wind parameter with different values, wherein a curve S1 and a curve S2 are the time sequence load of the tower bottom of the wind turbine generator obtained by simulating through Bladed software, and it can be seen from fig. 2 that the two curves have large phase deviation, which will seriously reduce the calculation accuracy of an interpolation result. The current methods for eliminating the phase deviation all lose the information of the time sequence load, and seriously affect the forecasting precision of the limit load and the fatigue load.

Disclosure of Invention

An exemplary embodiment of the present invention is directed to a method and an apparatus for eliminating a phase deviation of a timing curve, so as to overcome at least one of the above-mentioned disadvantages.

In one general aspect, there is provided a method of canceling a phase deviation of a timing curve, the method comprising: determining a plurality of first extreme points on a first timing curve; determining a plurality of second extreme points on a second timing curve; determining a plurality of control point pairs based on the plurality of first extreme points and the plurality of second extreme points; based on each control point pair, the second timing curve is shifted by taking the first timing curve as a reference so as to enable the phases of the first timing curve and the second timing curve to be consistent.

Alternatively, each control point pair may be composed of a first control point on the first timing curve and a second control point on the second timing curve, the first control point being one of the plurality of first extreme points, the second control point being one of the plurality of second extreme points that forms a control point pair with the one first extreme point, wherein the step of shifting the second timing curve based on each control point pair with the first timing curve as a reference so that the phases of the first and second timing curves are identical may include: the method includes dividing a first timing curve into a plurality of first curve segments based on all first control points on the first timing curve, dividing a second timing curve into a plurality of second curve segments based on all second control points on the second timing curve, and translating, for each first curve segment and a second curve segment corresponding in order to each first curve segment among the plurality of second curve segments, the second curve segment corresponding in order to the first curve segment with reference to the first curve segment so that phases of the first curve segment and the second curve segment corresponding in order to the first curve segment coincide.

Optionally, a second curve segment corresponding in order to any first curve segment may be translated with reference to the any first curve segment in such a way as to bring the any first curve segment into phase with the second curve segment corresponding in order to the any first curve segment: and elastically translating the second curve segment corresponding to any one first curve segment in sequence into a time period corresponding to any one first curve segment, so that the sampling time of the first sampling point of the second curve segment corresponding to any one first curve segment in sequence is consistent with the sampling time of the first sampling point of any one first curve segment, and the sampling time of the last sampling point of the second curve segment corresponding to any one first curve segment in sequence is consistent with the sampling time of the last sampling point of any one first curve segment.

Optionally, a second curve segment corresponding in order to any first curve segment may be translated with reference to the any first curve segment in such a way as to bring the any first curve segment into phase with the second curve segment corresponding in order to the any first curve segment: dividing a time period corresponding to any one first curve segment into N +1 equal parts, wherein N can be the number of intermediate sampling points contained in a second curve segment corresponding to any one first curve segment in sequence; and taking the sampling value of each intermediate sampling point on the second time sequence curve corresponding to any one first curve segment in sequence before translation as the sampling value of the halving point corresponding to each intermediate sampling point on the second time sequence curve corresponding to any one first curve segment in sequence after translation.

Optionally, the elimination method may further include: and obtaining the sampling value of the middle sampling point contained in any first curve segment by the second time sequence curve after translation through interpolation processing based on the sampling value of the first sampling point, the sampling value of the last sampling point and the sampling values of all the equant points on the second time sequence curve corresponding to any first curve segment in sequence after translation.

Alternatively, the plurality of first extreme points may include a plurality of first peak points and a plurality of first valley points, the plurality of second extreme points may include a plurality of second peak points and a plurality of second valley points, and the plurality of control point pairs may include a plurality of peak point pairs and a plurality of valley point pairs, wherein determining a plurality of control point pairs based on the plurality of first extreme points and the plurality of second extreme points may include: generating the plurality of peak point pairs based on the plurality of first peak points and the plurality of second peak points, and generating the plurality of valley point pairs based on the plurality of first valley points and the plurality of second valley points.

Optionally, the step of determining a plurality of first extreme points on the first timing curve may include: the method may further include selecting a plurality of first original peak points as the plurality of first peak points at predetermined intervals from all first original peak points on the first timing curve, selecting a plurality of first original valley points as the plurality of first valley points at predetermined intervals from all first original valley points on the first timing curve, and/or determining a plurality of second peak points on the second timing curve may include: and selecting a plurality of second original peak points as the plurality of second peak points at the preset interval from all the second original peak points on the second time sequence curve, and selecting a plurality of second original valley points as the plurality of second valley points at the preset interval from all the second original valley points on the second time sequence curve.

Alternatively, each first original peak point and each first original valley point may be determined by: determining all sampling points on the first time sequence curve; for any sampling point, if the sampling value of the sampling point is greater than the sampling value of the previous sampling point and greater than the sampling value of the next sampling point, determining the sampling point as a first original peak point, and if the sampling value of the sampling point is less than the sampling value of the previous sampling point and less than the sampling value of the next sampling point, determining the sampling point as a first original valley point, and/or determining each second original peak point and each second original valley point by the following method: determining all sampling points on the second time sequence curve; and for any sampling point, if the sampling value of the sampling point is greater than the sampling value of the last sampling point and greater than the sampling value of the next sampling point, determining the sampling point as a second original peak point, and if the sampling value of the sampling point is less than the sampling value of the last sampling point and less than the sampling value of the next sampling point, determining the sampling point as a second original valley point.

Alternatively, each pair of peak points may be determined by: for each first peak point on the first timing curve, searching whether a second peak point with a time interval with the sampling time of the first peak point being less than a preset time interval exists on the second timing curve, if the second peak point with the time interval with the sampling time of the first peak point being less than the preset time interval exists, determining the first peak point and the searched second peak point as a peak point pair, and/or determining each valley point pair by the following method: and for each first valley point on the first timing curve, searching whether a second valley point with the time interval of the sampling time of the first valley point smaller than a preset time interval exists on the second timing curve, and if the second valley point with the time interval of the sampling time of the first valley point smaller than the preset time interval exists, determining the first valley point and the searched second valley point as a valley point pair.

Optionally, the elimination method may further include: and storing the related information of the plurality of control point pairs into a data linked list, wherein the related information comprises the sampling time and the sampling value of each first control point in the plurality of control point pairs and the sampling time and the sampling value of each second control point.

Optionally, the elimination method may further include: if the last sampling point on the first timing curve is not included in all the first control points of the plurality of control point pairs, the sampling time and the sampling value of the last sampling point on the first timing curve can be added into the data linked list, and if the last sampling point on the second timing curve is not included in all the second control points of the plurality of control point pairs, the sampling time and the sampling value of the last sampling point on the second timing curve can be added into the data linked list, wherein the relevant information of the plurality of control point pairs is obtained from the data linked list, and the second timing curve is shifted by taking the first timing curve as a reference based on each control point pair, so that the phases of the first timing curve and the second timing curve are consistent.

Optionally, the first time sequence curve may be a load time sequence curve corresponding to an upper boundary point adjacent to a numerical value of any preset wind parameter under the working condition to be solved, and the second time sequence curve may be a load time sequence curve corresponding to a lower boundary point adjacent to the numerical value of any preset wind parameter under the working condition to be solved, where the elimination method may further include: and obtaining a load time sequence curve corresponding to the numerical value of any preset wind parameter under the working condition to be solved through interpolation processing based on the first time sequence curve and the second time sequence curve.

In another general aspect, there is provided an apparatus for canceling a phase deviation of a timing curve, the apparatus comprising: the first extreme value determining module is used for determining a plurality of first extreme values on the first time sequence curve; the second extreme value determining module is used for determining a plurality of second extreme value points on the second time sequence curve; a control point pair determination module that determines a plurality of control point pairs based on the plurality of first extreme points and the plurality of second extreme points; and the phase deviation elimination module is used for shifting the second time sequence curve by taking the first time sequence curve as a reference on the basis of each control point pair so as to enable the phases of the first time sequence curve and the second time sequence curve to be consistent.

Alternatively, each control point pair may be composed of a first control point on a first timing curve and a second control point on a second timing curve, the first control point being one of the plurality of first extreme points, the second control point being one of the plurality of second extreme points that forms a control point pair with the one first extreme point, wherein the phase deviation elimination module may include: the first time sequence curve is divided into a plurality of first curve segments based on all first control points on the first time sequence curve, the second time sequence curve is divided into a plurality of second curve segments based on all second control points on the second time sequence curve, and the second curve segments corresponding to the first curve segments in sequence are translated by taking the first curve segments as a reference so that the phases of the first curve segments and the second curve segments corresponding to the first curve segments in sequence are consistent.

Optionally, the translation submodule may translate a second curve segment corresponding in order to any first curve segment with reference to the any first curve segment in such a way as to bring the phase of the any first curve segment and the second curve segment corresponding in order to the any first curve segment into agreement: and elastically translating the second curve segment corresponding to any one first curve segment in sequence into a time period corresponding to any one first curve segment, so that the sampling time of the first sampling point of the second curve segment corresponding to any one first curve segment in sequence is consistent with the sampling time of the first sampling point of any one first curve segment, and the sampling time of the last sampling point of the second curve segment corresponding to any one first curve segment in sequence is consistent with the sampling time of the last sampling point of any one first curve segment.

Optionally, the translation submodule may translate a second curve segment corresponding in order to any first curve segment with reference to the any first curve segment in such a way as to bring the phase of the any first curve segment and the second curve segment corresponding in order to the any first curve segment into agreement: dividing a time period corresponding to any one first curve segment into N +1 equal parts, wherein N is the number of middle sampling points contained in a second curve segment corresponding to any one first curve segment in sequence; and taking the sampling value of each intermediate sampling point on the second time sequence curve corresponding to any one first curve segment in sequence before translation as the sampling value of the halving point corresponding to each intermediate sampling point on the second time sequence curve corresponding to any one first curve segment in sequence after translation.

Optionally, the eliminating device may further include: and the interpolation module is used for obtaining the sampling value of the middle sampling point contained in any first curve segment by the translated second time sequence curve through interpolation processing based on the sampling value of the first sampling point, the sampling value of the last sampling point and the sampling values of all the equant points on the translated second time sequence curve corresponding to any first curve segment in sequence.

Alternatively, the plurality of first extreme points may include a plurality of first peak points and a plurality of first valley points, the plurality of second extreme points may include a plurality of second peak points and a plurality of second valley points, and the plurality of control point pairs may include a plurality of peak point pairs and a plurality of valley point pairs, wherein the control point pair determination module may generate the plurality of peak point pairs based on the plurality of first peak points and the plurality of second peak points, and generate the plurality of valley point pairs based on the plurality of first valley points and the plurality of second valley points.

Optionally, the first extreme value determining module may select a plurality of first original peak points as the plurality of first peak points at predetermined intervals from all first original peak points on the first timing curve, and select a plurality of first original valley points as the plurality of first valley points at the predetermined intervals from all first original valley points on the first timing curve, and/or the second extreme value determining module may select a plurality of second original peak points as the plurality of second peak points at the predetermined intervals from all second original peak points on the second timing curve, and select a plurality of second original valley points as the plurality of second valley points at the predetermined intervals from all second original valley points on the second timing curve.

Alternatively, the first extreme value determining module may determine each first original peak point and each first original valley point by: determining all sampling points on the first time sequence curve; for any sampling point, if the sampling value of the sampling point is greater than the sampling value of the previous sampling point and greater than the sampling value of the next sampling point, determining the sampling point as a first original peak point, and if the sampling value of the sampling point is less than the sampling value of the previous sampling point and less than the sampling value of the next sampling point, determining the sampling point as a first original valley point, and/or determining each second original peak point and each second original valley point by a second polarity determination module in the following way: determining all sampling points on the second time sequence curve; and for any sampling point, if the sampling value of the sampling point is greater than the sampling value of the last sampling point and greater than the sampling value of the next sampling point, determining the sampling point as a second original peak point, and if the sampling value of the sampling point is less than the sampling value of the last sampling point and less than the sampling value of the next sampling point, determining the sampling point as a second original valley point.

Alternatively, the control-point pair determination module may determine each pair of peak points by: for each first peak point on the first timing curve, searching whether a second peak point with a time interval of the sampling time of the first peak point being less than a preset time interval exists on the second timing curve, if the second peak point with the time interval of the sampling time of the first peak point being less than the preset time interval exists, determining the first peak point and the searched second peak point as a peak point pair, and/or determining each valley point pair by the control point pair determining module in the following way: and for each first valley point on the first timing curve, searching whether a second valley point with the time interval of the sampling time of the first valley point smaller than a preset time interval exists on the second timing curve, and if the second valley point with the time interval of the sampling time of the first valley point smaller than the preset time interval exists, determining the first valley point and the searched second valley point as a valley point pair.

Optionally, the eliminating device may further include: and the storage module stores the related information of the plurality of control point pairs into a data link list, wherein the related information comprises the sampling time and the sampling value of each first control point in the plurality of control point pairs and the sampling time and the sampling value of each second control point.

Alternatively, if the last sample point on the first timing curve is not included in all the first control points of the plurality of control point pairs, the data link list establishing module may add the sample time and the sample value of the last sample point on the first timing curve to the data link list, if the last sample point on the second timing curve is not included in all second control points of the plurality of control point pairs, the data link list creation module may add the sample time and the sample value of the last sample point on the second timing curve to the data link list, wherein the phase offset cancellation module is capable of obtaining the related information of the plurality of control point pairs from the data link list, and based on each control point pair, and shifting the second time sequence curve by taking the first time sequence curve as a reference so as to enable the phases of the first time sequence curve and the second time sequence curve to be consistent.

Optionally, the first time sequence curve may be a load time sequence curve corresponding to an upper boundary point adjacent to a numerical value of any preset wind parameter under the working condition to be solved, and the second time sequence curve may be a load time sequence curve corresponding to a lower boundary point adjacent to a numerical value of any preset wind parameter under the working condition to be solved, where the eliminating device may further include: and the time sequence load determining module is used for obtaining a load time sequence curve corresponding to the numerical value of any preset wind parameter under the working condition to be solved through interpolation processing based on the first time sequence curve and the second time sequence curve.

In another general aspect, there is provided a computer-readable storage medium storing a computer program which, when executed by a processor, implements the method of canceling a phase deviation of a timing curve described above.

In another general aspect, there is provided a computing device, comprising: a processor; a memory storing a computer program that, when executed by the processor, implements the method for eliminating the phase deviation of the timing curve described above.

By adopting the method and the device for measuring the tower clearance of the wind generating set, the tower clearance of the wind generating set can be determined in real time, so that the condition that the blades sweep the tower is effectively avoided.

Drawings

The above and other objects and features of the present invention will become more apparent from the following description taken in conjunction with the accompanying drawings, in which:

FIG. 1 shows a global schematic of two load timing curves at different values of a single wind parameter;

FIG. 2 shows a partial schematic view of two load timing curves shown in FIG. 1;

fig. 3 illustrates a flowchart of a method of eliminating a phase deviation of a timing curve according to an exemplary embodiment of the present invention;

fig. 4 shows a general flow diagram of a method of cancellation of phase deviations of a timing curve according to an exemplary embodiment of the invention;

FIG. 5 illustrates a schematic distribution of the first extreme point and the second extreme point after screening according to an exemplary embodiment of the present invention;

FIG. 6 illustrates a distribution diagram of pairs of control points of a first timing curve and a second timing curve according to an exemplary embodiment of the present invention;

FIG. 7 shows a flowchart of the step of translating a second curve segment in accordance with an exemplary embodiment of the present invention;

FIG. 8 illustrates a schematic diagram of translating a second curve segment based on control point pairs in accordance with an exemplary embodiment of the present invention;

FIG. 9 illustrates a density of 0.9kg/m at air according to an exemplary embodiment of the present invention³And 1.0kg/m³Under the working condition of (1), a global schematic diagram of a tower bottom Mx load time sequence;

FIG. 10 illustrates a partial schematic diagram of the header shift effect of the tower bottom Mx loading timing sequence illustrated in FIG. 8, according to an exemplary embodiment of the present invention;

FIG. 11 illustrates a partial schematic of the mid-bias effect of the tower bottom Mx load timing sequence illustrated in FIG. 8 in accordance with an exemplary embodiment of the present invention;

FIG. 12 illustrates a partial schematic of the tail shift effect of the tower bottom Mx load timing illustrated in FIG. 8 in accordance with an exemplary embodiment of the present invention;

fig. 13 is a block diagram illustrating a phase deviation removal apparatus of a timing curve according to an exemplary embodiment of the present invention;

fig. 14 illustrates a block diagram of a phase offset cancellation module according to an exemplary embodiment of the present invention.

Detailed Description

Various example embodiments will now be described more fully with reference to the accompanying drawings, in which some example embodiments are shown.

Fig. 3 is a flowchart of a method of removing a phase deviation of a timing curve according to an exemplary embodiment of the present invention. Fig. 4 shows a general flowchart of a method of eliminating a phase deviation of a timing curve according to an exemplary embodiment of the present invention. The process of eliminating the phase offset of the two timing curves is described below in conjunction with the general flow chart of fig. 4.

Referring to fig. 3, in step S10, a plurality of first extreme points on the first timing curve are determined.

For example, all the first original extreme points on the first timing curve may be found, and the plurality of first extreme points are obtained by screening all the first original extreme points. This is because the number of extreme points found is very large because of the large fluctuation of the timing curve, and therefore the extreme points found need to be filtered. As an example, all the first original extreme points may be filtered by setting the distance between two adjacent first extreme points.

Here, the first primitive extremum points may include a plurality of first primitive peak points and a plurality of first primitive valley points, and the plurality of first extremum points may include a plurality of first peak points and a plurality of first valley points.

For example, the plurality of first peak points may be obtained by filtering the plurality of first original peak points, and the plurality of first valley points may be obtained by filtering the plurality of first original valley points.

The process of determining each first raw peak point and each first raw valley point is described below.

Determining all sampling points on the first time sequence curve, and for any sampling point on the first time sequence curve, if the sampling value of the sampling point is greater than the sampling value of the last sampling point and greater than the sampling value of the next sampling point, determining the sampling point as a first original peak point, and if the sampling value of the sampling point is less than the sampling value of the last sampling point and less than the sampling value of the next sampling point, determining the sampling point as a first original valley point.

Here, if the arbitrary sampling point is the first intermediate sampling point (i.e., one of the first intermediate sampling points that is most forward in time series) of all the intermediate sampling points, the last sampling point of the first intermediate sampling point is the first sampling point on the first timing curve, and the next sampling point of the first intermediate sampling point is the second intermediate sampling point in time series. If the any one sampling point is the last intermediate sampling point (i.e., the one that is the most backward in time sequence) of all the intermediate sampling points, the last sampling point of the last intermediate sampling point is the second last intermediate sampling point in time sequence, and the next sampling point of the last intermediate sampling point is the last sampling point on the first time sequence curve.

The following describes a process of screening a plurality of first original peak points and a plurality of first original valley points.

For example, after obtaining the plurality of first original peak points, the plurality of first peak points may be obtained by performing a screening in the following manner: and selecting a plurality of first original peak points as a plurality of first peak points at preset intervals from all the first original peak points on the first timing curve.

For example, after obtaining the plurality of first original valley points, the plurality of first valley points may be obtained by performing a screening in the following manner: selecting a plurality of first original valley points as a plurality of first valley points at predetermined intervals from all the first original valley points on the first timing curve.

In a preferred embodiment, the value of the predetermined interval may include the number of sampling points, for example, if the value of the predetermined interval is 3 sampling points, the above-mentioned screening process may be: and selecting a first original peak point as a first peak point every 3 sampling points, and selecting a first original valley point as a first valley point every 3 sampling points. As an example, the 1 st, 5 th, 9 … … th first original peak points may be selected as the first peak points, and the 1 st, 5 th, 9 … … th first original valley points may be selected as the first valley points.

It should be understood that the above-mentioned manner of taking the number of sampling points as the value of the predetermined interval is only an example, and other parameters may also be selected as the value of the predetermined interval, such as time length, distance, and the like. In addition, a randomly selected mode can be adopted to screen the plurality of first original peak points and the plurality of first original valley points.

In step S20, a plurality of second extreme points on the second timing curve are determined.

For example, all the second original extreme points on the second timing curve may be found, and the plurality of second extreme points may be obtained by screening all the second original extreme points. As an example, all the second original extreme points may be screened by setting the distance between adjacent second extreme points.

Here, the second original extreme point may include a plurality of second original peak points and a plurality of second original valley points, and the plurality of second extreme points may include a plurality of second peak points and a plurality of second valley points.

For example, the plurality of second peak points may be obtained by filtering the plurality of second original peak points, and the plurality of second valley points may be obtained by filtering the plurality of second original valley points.

The process of determining each second original peak point and each second original valley point is described below.

And determining all sampling points on the second timing curve, and for any sampling point on the second timing curve, if the sampling value of the sampling point is greater than the sampling value of the last sampling point and greater than the sampling value of the next sampling point, determining the sampling point as a second original peak point, and if the sampling value of the sampling point is less than the sampling value of the last sampling point and less than the sampling value of the next sampling point, determining the sampling point as a second original valley point.

The process of screening the plurality of second original peak points and the plurality of second original valley points is described below.

For example, after obtaining the plurality of second original peak points, the plurality of second peak points may be obtained by performing a screening in the following manner: and selecting a plurality of second original peak points as a plurality of second peak points at preset intervals from all the second original peak points on the second time sequence curve.

For example, after obtaining the plurality of second original valley points, the plurality of second valley points may be obtained by performing a screening in the following manner: and selecting a plurality of second original valley points as a plurality of second valley points at preset intervals from all the second original valley points on the second time sequence curve. Here, in the above-described embodiment, the first peak point, the first valley point, the second peak point, and the second valley point are selected at the same predetermined interval, but the present invention is not limited thereto, and the predetermined interval for selecting each extreme point may be different, or it is only necessary to ensure that the predetermined interval for selecting the first peak point and the second peak point is the same and the predetermined interval for selecting the first valley point and the second valley point is the same, so as to determine the control point pair subsequently.

Fig. 5 shows a distribution diagram of the first extreme point and the second extreme point after the screening according to the exemplary embodiment of the present invention.

As can be seen from fig. 5, after the first original peak points and the first original valley points are screened, the number of the first peak points and the second valley points is significantly reduced.

That is, the purpose of the above screening process is to reduce the number of the first extreme points and the second extreme points, thereby reducing the amount of calculation in the phase deviation elimination process, but to reserve a certain number of the first extreme points and the second extreme points, so as to ensure that the accuracy of the timing curve after the phase deviation elimination can meet the requirement.

The skilled person can adjust the values of the above predetermined intervals according to actual requirements, so that the number of the screened extreme points can meet the actual engineering requirements. If it is desired to improve the accuracy of eliminating the phase deviation, the values of the predetermined intervals may be decreased to obtain a greater number of first extreme points and second extreme points, which of course increases the amount of calculation as the number of first extreme points and second extreme points increases.

Returning to fig. 3, in step S30, a plurality of control point pairs are determined based on the plurality of first extreme points on the first timing curve and the plurality of second extreme points on the second timing curve. As an example, the plurality of control point pairs may include a plurality of peak point pairs and a plurality of valley point pairs.

For example, in step S30, a plurality of peak point pairs may be generated based on a plurality of first peak points on the first timing curve and a plurality of second peak points on the second timing curve, and a plurality of valley point pairs may be generated based on a plurality of first valley points on the first timing curve and a plurality of second valley points on the second timing curve.

Preferably, each peak point pair may be determined in the following manner.

And for each first peak point on the first timing curve, searching whether a second peak point with a time interval with the sampling time of the first peak point smaller than a preset time interval exists on the second timing curve, and if the second peak point with the time interval with the sampling time of the first peak point smaller than the preset time interval exists, determining the first peak point and the searched second peak point as a peak point pair.

For example, for any first peak point on the first timing curve, the difference between the sampling time of each second peak point on the second timing curve and the sampling time of the any first peak point is calculated, and if there is a second peak point whose difference is smaller than a predetermined time interval, the found second peak point and the any first peak point are determined as a peak point pair. And sequentially executing the searching process aiming at each first peak point to obtain a plurality of peak point pairs.

Here, by traversing all the first peak points on the first timing curve, the number of peak point pairs of the two timing curves can be increased, and the effect of translation can be improved.

Preferably, each valley point pair may be determined in the following manner.

And for each first valley point on the first timing curve, searching whether a second valley point with the time interval of the sampling time of the first valley point smaller than a preset time interval exists on the second timing curve, and if the second valley point with the time interval of the sampling time of the first valley point smaller than the preset time interval exists, determining the first valley point and the searched second valley point as a valley point pair.

For example, for any first valley point on the first timing curve, a difference between a sampling time of each second valley point on the second timing curve and a sampling time of the any first valley point is calculated, and if there is a second valley point whose difference is smaller than a predetermined time interval, the found second valley point and the any first valley point are determined as a valley point pair. And executing the searching process for each first valley point in sequence to obtain a plurality of valley point pairs.

Here, by traversing all the first valley points on the first timing curve, the number of valley point pairs of the two timing curves can be increased, and the effect of translation can be improved.

It should be understood that if the value of the predetermined time interval is too large, the correlation between the found control point pairs corresponding to the two timing curves cannot be guaranteed, and if the value is too small, the number of the found control point pairs corresponding to the two timing curves is small, which is not favorable for the subsequent elimination process of the phase deviation. Preferably, the value of the predetermined time interval may be 0.2 seconds (fig. 6 shows a distribution diagram of a plurality of control point pairs determined when the value is 0.2 seconds), but the present invention is not limited thereto, and a person skilled in the art may adjust the value of the predetermined time interval according to actual engineering requirements.

Preferably, the method of removing the phase deviation of the timing curve according to an exemplary embodiment of the present invention may further include: and storing the related information of the control point pairs into a data linked list. Here, each control point pair may be composed of a first control point on the first timing curve, which is one of the plurality of first extreme points on the first timing curve, and a second control point on the second timing curve, which is a second extreme point forming a control point pair with the one first extreme point among the plurality of second extreme points on the second timing curve. For example, the first control point may be a first peak point forming a pair of peak points and the second control point may be a second peak point forming the pair of peak points, or the first control point may be a first valley point forming a pair of valley points and the second control point may be a second valley point forming the pair of valley points.

In this case, the related information of the plurality of control point pairs may include a sampling time and a sampling value of each first control point of the plurality of control point pairs, and a sampling time and a sampling value of each second control point.

That is, after the plurality of control point pairs are determined, the sampling time and the sampling value of each first control point and the sampling time and the sampling value of each second control point in the plurality of control point pairs are saved into the data link list with respect to the control point pairs.

It should be understood that the total time length of the first and second timing curves before the shift is the same (e.g., timing curves within 0-600 seconds of each other), and the time interval between each sampling point of the first and second timing curves is also the same, i.e., the sampling period of the two timing curves is the same.

The above process of extracting the extreme points and the process of determining the control point pairs are prepared for the subsequent phase deviation elimination, and after the extreme points are extracted and the control point pairs are determined, the total time length of the first timing curve and the second timing curve should be ensured to be consistent before the translation, that is, the value ranges of the abscissas of the first timing curve and the second timing curve are ensured to be the same.

After determining the plurality of control point pairs, if the total time length of the first timing curve and/or the second timing curve changes, i.e. the first sampling point or the last sampling point of the timing curve is not included in the plurality of control point pairs, the total time length of the timing curve is reduced. At this time, the total time length of the timing curves needs to be complemented so that the total time length of the two timing curves is still complete after the extreme point extraction and control point pair determination process. Particularly, when the time sequence curve is the load time sequence curve of the wind turbine generator, the time sequence load information can be ensured not to be lost, and the time sequence information of the load can be completely reserved, so that the prediction accuracy of the limit load and the fatigue load can be improved.

For example, when any one of the sampling points is the first sampling point on the first timing curve, the sampling value of the first sampling point only needs to be compared with the sampling value of the next sampling point to determine whether the sampling point is the first original valley point or the first original peak point. And (4) judging the first sampling point on the second time sequence curve in the same way, so that the first sampling points of the two time sequence curves are either a valley point or a peak point. Because the trends of the two timing curves at the initial time are consistent (the characteristics of the timing load curves), the types of the first extreme points of the two timing curves are consistent (namely, both the first extreme points and both the second extreme points are peak points or both the first extreme points and the sampling time of the first sampling point is 0, so that the first sampling points on the two timing curves are necessarily a control point pair (possibly a peak point pair or a valley point pair). That is, the data chain table must include the sampling time and the sampling value of the first sampling point of the two timing curves. But whether the last sample point of the two timing curves is an extreme point and can become a control point pair is relatively random.

In this case, to ensure the coincidence of the total time lengths of the two timing curves, the method of eliminating the phase deviation of the timing curves according to the exemplary embodiment of the present invention may further include: and if the first sampling point and/or the last sampling point on the first timing curve are not included in all the first control points of the plurality of control point pairs, adding the sampling time and the sampling value of the first sampling point and/or the sampling time and the sampling value of the last sampling point on the first timing curve into the data link list.

And if the first sampling point and/or the last sampling point on the second timing curve are not included in all the second control points of the plurality of control point pairs, adding the sampling time and the sampling value of the first sampling point and/or the sampling time and the sampling value of the last sampling point on the second timing curve into the data link list.

Based on the above, the related information of the plurality of control point pairs can be obtained from the data linked list, and the second timing curve is shifted based on the first timing curve based on each control point pair, so that the phases of the first timing curve and the second timing curve are consistent.

In step S40, the second timing curve is shifted based on the first timing curve so that the phases of the first timing curve and the second timing curve match with each other on a control point-by-control-point basis.

It should be understood that the phase offset cancellation method shown in fig. 3 is a method for canceling a phase offset between two timing curves, and the first timing curve may refer to one of the two timing curves, and the second timing curve may refer to the other of the two timing curves.

That is, one of the first and second timing curves may be shifted with reference to the other of the first and second timing curves in step S40 so that the phases of the first and second timing curves coincide. That is, in step S40, the first time series curve may be shifted with respect to the second time series curve so that the phases of the first time series curve and the second time series curve match.

In a preferred embodiment, the phase deviation between the first timing curve and the second timing curve can be eliminated by dividing the first timing curve and the second timing curve respectively based on each control point pair and then shifting the divided curve segments respectively.

The process of performing the shift for each divided curve segment to eliminate the phase deviation of the two timing curves is described below with reference to fig. 7.

FIG. 7 shows a flowchart of the step of translating the second curve segment according to an exemplary embodiment of the present invention.

Referring to fig. 7, in step S401, the first timing curve is divided into a plurality of first curve segments based on all the first control points on the first timing curve.

That is, the first timing curve is divided into a plurality of first curve segments with all the first control points as dividing points. For example, when a control point pairs are determined, there are a first control points, and at this time, the first timing curve is divided into a +1 first curve segments with the a first control points as dividing points.

In step S402, the second timing curve is divided into a plurality of second curve segments based on all of the second control points on the second timing curve.

That is, the second timing curve is divided into a plurality of second curve segments with all the second control points as dividing points. For example, when a control point pairs are determined, there are a second control points, and at this time, the second timing curve is divided into a +1 second curve segments with the a second control points as dividing points.

As a result, the number of the second curve segments obtained by the division is equal to the number of the first curve segments, and the plurality of second curve segments and the plurality of first curve segments are sequentially associated with each other.

In step S403, for each first curve segment and a second curve segment corresponding in order to each first curve segment among a plurality of second curve segments, the second curve segment corresponding in order to the first curve segment is shifted with reference to the first curve segment so that phases of the first curve segment and the second curve segment corresponding in order to the first curve segment are consistent.

Preferably, the second curve segment corresponding in order to any one of the first curve segments may be translated with reference to the any one of the first curve segments in such a manner that the phases of the any one of the first curve segments and the second curve segment corresponding in order to the any one of the first curve segments coincide: and elastically translating the second curve segment corresponding to any one first curve segment in sequence into a time period corresponding to any one first curve segment, so that the sampling time of the first sampling point of the second curve segment corresponding to any one first curve segment in sequence is consistent with the sampling time of the first sampling point of any one first curve segment, and the sampling time of the last sampling point of the second curve segment corresponding to any one first curve segment in sequence is consistent with the sampling time of the last sampling point of any one first curve segment.

FIG. 8 illustrates a schematic diagram of translating a second curve segment based on a control point pair according to an exemplary embodiment of the present invention.

In this example, the trends of the timing curve 1 and the timing curve 2 are the same, and there is a relatively significant phase deviation, and fig. 8 is a partial schematic diagram of the timing curve 1 and the timing curve 2. Assume that the first control point on timing curve 1 in the first control point pair is at t1, the second control point on timing curve 2 in the first control point pair is at t2, and the first sample points of timing curve 1 and timing curve 2 coincide, both at t 0.

In this case, the phase deviation is eliminated by translating the timing curve 2 with reference to the timing curve 1, and the curve segment of the timing curve 2 between t0 and t2 is elastically translated to between t0 and t1, that is, the second curve segment is elastically translated to the time period corresponding to the first curve segment, with reference to the spring characteristic during the translation). Accordingly, the curve segment of the timing curve 2 between t2 and t4 may also be elastically translated between t1 and t3 with reference to the spring characteristic. Alternatively, the curve segment of the timing curve 2 between t1 (after translation t2 coincides with t 1) to t4 may be elastically translated between t1 to t3 based on the result of the translation of the curve segment of the timing curve 2 between t0 to t 2.

It should be understood that the translation shown in fig. 8 is only an example, and those skilled in the art can select a specific translation and a translation sequence according to requirements.

In a preferred embodiment, the second curve segment corresponding in order to any one of the first curve segments may be translated with reference to the any one of the first curve segments in such a manner that the any one of the first curve segments coincides in phase with the second curve segment corresponding in order to the any one of the first curve segments: and dividing the time period corresponding to any one first curve segment into N +1 equal parts, wherein N is the number of intermediate sampling points contained in a second curve segment corresponding to any one first curve segment in sequence, and taking the sampling value of each intermediate sampling point on a second time sequence curve corresponding to any one first curve segment in sequence before translation as the sampling value of an equal division point corresponding to each intermediate sampling point on a second time sequence curve corresponding to any one first curve segment in sequence after translation.

For example, taking the timing curve 1 and the timing curve 2 shown in fig. 8 as an example, assuming that 4 intermediate sample points are included in the curve segment of t0 to t2 of the timing curve 2 (i.e., 4 intermediate sample points are included in the curve segment of t0 to t2 except for the first sample point t0 and the last sample point t 2), the time interval between t0 to t1 of the timing curve 1 is divided equally by 4+1, and at this time, there are 4 equally divided points on the curve segment of t0 to t1 of the timing curve 1, and the sample values at the 4 equally divided points on the curve segment of t0 to t2 of the timing curve 2 after the translation are the sample values at the 4 intermediate sample points on the curve segment of t0 to t2 of the timing curve 2 before the translation.

Here, the second timing curve can be made to coincide in phase with the first timing curve by the above-described shifting process, that is, two timing curves having a small phase deviation can be obtained. For the case that the time sequence curve is a load time sequence curve of the wind turbine generator, if two time sequence curves with the same phase are required to be applied to interpolation processing, interpolation processing is generally performed on a sampling value on a first time sequence curve and a sampling value on a second time sequence curve at the same sampling time, but after the translation, it cannot be guaranteed that the translated second time sequence curve has a sampling value at the sampling time of each sampling point (the translated second time sequence curve has a sampling value at each equal division point), and at this time, the sampling value at the sampling time of each sampling point of the first time sequence curve needs to be obtained based on the existing sampling value on the translated second time sequence curve.

For example, the sample values of the translated second timing curve at the sampling times of the intermediate sample points included in the any first curve segment may be obtained by interpolation processing based on the sample value of the first sample point, the sample value of the last sample point, and the sample values at the respective divided points on the translated second timing curve corresponding in order to the any first curve segment.

Taking fig. 8 as an example, the sampling time of the sampling point between t0 and t1 of the timing curve 1 and the time corresponding to the bisector between t0 and t2 of the translated timing curve 2 may not coincide, and at this time, the sampling value of the curve segment between t0 and t2 of the translated timing curve 2 at the sampling time of each sampling point may be obtained by interpolation, so that both the timing curve 1 and the translated timing curve 2 have the sampling value at the same sampling time, so as to facilitate the subsequent interpolation processing.

It should be understood that the above-described manner of obtaining the sample value at the sampling time of each sample point by interpolation processing is merely an example, and the present invention is not limited thereto, and the sample value at the sampling time of each sample point may be obtained in other manners, for example, by data fitting.

In a preferred embodiment, the first and second timing curves may both be load timing curves of the wind turbine.

For example, a load sample library may be established in advance, in which wind parameter value combinations and time series loads corresponding to the wind parameter value combinations are stored, the wind parameter value combinations in the load sample library being obtained by arranging and combining a plurality of value points corresponding to each preset wind parameter.

The interpolation process is as follows: based on a pre-established load sample library, the combination of different wind parameter values is used as input, time sequence loads (namely load time sequence curves) corresponding to upper and lower boundary points of the numerical value of a single preset wind parameter are found out in the load sample library, and the time sequence loads under the numerical values of the preset wind parameters under the working condition to be solved are finally obtained through continuous elimination and linear interpolation.

Preferably, the plurality of preset wind parameters may be determined by: according to the linear correlation between the wind parameters and the load, selecting a preset number of wind parameters with the linear correlation ranked in the top from all the wind parameters related to the load as a plurality of preset wind parameters. Namely, the wind parameter with high linear correlation is selected as the preset wind parameter to be used in the subsequent time sequence load analysis processing. Here, the wind parameter may refer to a parameter for describing wind conditions, and may include, as examples, wind speed, turbulence, inflow angle, air density, wind shear. In a preferred embodiment, the preset wind parameters may include inflow angle, air density, wind shear.

As an example, the first time sequence curve may be a first load time sequence curve corresponding to an upper boundary point adjacent to a numerical value of any one of the preset wind parameters under the to-be-solved condition, and the second time sequence curve may be a second load time sequence curve corresponding to a lower boundary point adjacent to a numerical value of any one of the preset wind parameters under the to-be-solved condition.

For example, a plurality of value points corresponding to any one preset wind parameter may be searched from the load sample library, and two value points closest to the value of any one preset wind parameter are selected from the searched value points as two boundary value points adjacent to the value of the preset wind parameter. Here, the two boundary value points of any one of the preset wind parameters may include an upper boundary point that is greater than the value of any one of the preset wind parameters and a lower boundary point that is less than the value of any one of the preset wind parameters. Then, a first load time sequence curve corresponding to the working condition of the upper boundary point of the numerical value of any preset wind parameter and a second load time sequence curve corresponding to the working condition of the lower boundary point of the numerical value of any preset wind parameter are searched from the load sample library. And obtaining a load time sequence curve corresponding to the numerical value of any preset wind parameter under the working condition to be solved through interpolation processing based on the first load time sequence curve and the second load time sequence curve.

Taking a certain wind turbine as an example, a specific embodiment of eliminating the phase deviation of the timing curve will be described with reference to fig. 9 to 12.

In an exemplary embodiment of the present invention, a 2.5MW wind turbine is taken as a research object, and FIGS. 9 to 12 show that the wind turbine has a wind speed of 4.5m/s, a turbulence level of 15%, and an air density of 0.9kg/m³And 1.0kg/m³Under the working condition of (2), the offset effect of the Mx load time sequence at the bottom of the tower.

Here, the time-series curve s1 shows that the air density is 0.9kg/m³The time series curve s2 shows that the air density is 1.0kg/m³By way of example, the load time curve at an air density of 0.9kg/m can be obtained using a Bladed software simulation³And 1.0kg/m³The column bottom Mx load timing under the operating conditions of (1). The invention is not limited thereto and the load timing may be obtained in other ways, for example, other dynamic simulation software may be used, such as Fast, Hawc2, Simpack, etc. s2_ new is the load timing curve translated to the timing curve s 2. By adopting the method for eliminating the phase deviation of the time sequence curve, more control point pairs can be found out, the correspondence of extreme points of the two load time sequence curves can be ensured, and the correspondence degree of the curves between the extreme points can be greatly improved. Verified, the elimination method is used for key components of the wind turbine generator under other wind parametersThe load has better offset effect.

Fig. 13 is a block diagram illustrating a phase deviation removing apparatus of a timing curve according to an exemplary embodiment of the present invention.

As shown in fig. 13, the apparatus for removing a phase deviation of a timing curve according to an exemplary embodiment of the present invention includes: a first extreme value determining module 10, a second extreme value determining module 20, a control point pair determining module 30 and a phase deviation eliminating module 40.

Specifically, the first extreme value determining module 10 determines a plurality of first extreme points on the first timing curve.

For example, the first extreme value determining module 10 may find all the first original extreme values on the first timing curve, and obtain a plurality of first extreme values by filtering all the first original extreme values.

For example, the first extreme value determining module 10 may obtain a plurality of first peak points by filtering a plurality of first original peak points, and the first extreme value determining module 10 may obtain a plurality of first valley points by filtering a plurality of first original valley points.

The first extreme value determining module 10 may select a plurality of first original peak points as the plurality of first peak points at predetermined intervals from all first original peak points on the first timing curve, and select a plurality of first original valley points as the plurality of first valley points at predetermined intervals from all first original valley points on the first timing curve.

The first extreme value determining module 10 may determine each first original peak point and each first original valley point by: determining all sampling points on the first time sequence curve; for any sampling point, if the sampling value of the sampling point is larger than the sampling value of the last sampling point and larger than the sampling value of the next sampling point, determining the sampling point as a first original peak point, and if the sampling value of the sampling point is smaller than the sampling value of the last sampling point and smaller than the sampling value of the next sampling point, determining the sampling point as a first original valley point.

The second extreme determining module 20 determines a plurality of second extreme points on the second timing curve.

For example, the second extreme determining module 20 may find all the second original extreme points on the second timing curve, and obtain a plurality of second extreme points by filtering all the second original extreme points.

For example, the second polarity determination module 20 may obtain a plurality of second peak points by filtering a plurality of second original peak points, and the second polarity determination module 20 may obtain a plurality of second valley points by filtering a plurality of second original valley points.

The second polarity determining module 20 may select a plurality of second original peak points as the plurality of second peak points at predetermined intervals from all second original peak points on the second timing curve, and select a plurality of second original valley points as the plurality of second valley points at predetermined intervals from all second original valley points on the second timing curve.

The second extreme determination module 20 may determine each second original peak point and each second original valley point by: determining all sampling points on the second timing curve; and for any sampling point, if the sampling value of the sampling point is greater than the sampling value of the last sampling point and greater than the sampling value of the next sampling point, determining the sampling point as a second original peak point, and if the sampling value of the sampling point is less than the sampling value of the last sampling point and less than the sampling value of the next sampling point, determining the sampling point as a second original valley point.

Control point pair determination module 30 determines a plurality of control point pairs based on the plurality of first extreme points and the plurality of second extreme points. As an example, the plurality of control point pairs may include a plurality of peak point pairs and a plurality of valley point pairs.

For example, control point pair determination module 30 may generate a plurality of peak point pairs based on a plurality of first peak points on a first timing curve and a plurality of second peak points on a second timing curve, and generate a plurality of valley point pairs based on a plurality of first valley points on the first timing curve and a plurality of second valley points on the second timing curve.

For example, the control point pair determination module 30 may determine each peak point pair by: and for each first peak point on the first timing curve, searching whether a second peak point with a time interval with the sampling time of the first peak point smaller than a preset time interval exists on the second timing curve, and if the second peak point with the time interval with the sampling time of the first peak point smaller than the preset time interval exists, determining the first peak point and the searched second peak point as a peak point pair.

For example, the control point pair determination module 30 may determine each valley point pair by: and for each first valley point on the first timing curve, searching whether a second valley point with the time interval of the sampling time of the first valley point smaller than a preset time interval exists on the second timing curve, and if the second valley point with the time interval of the sampling time of the first valley point smaller than the preset time interval exists, determining the first valley point and the searched second valley point as a valley point pair.

Preferably, the apparatus for removing a phase deviation of a timing curve according to an exemplary embodiment of the present invention may further include: and the storage module stores the related information of the control point pairs into a data linked list. Here, the related information includes a sampling time and a sampling value of each first control point of the plurality of control point pairs, and a sampling time and a sampling value of each second control point.

And if the first sampling point and/or the last sampling point on the first timing curve are not included in all the first control points of the plurality of control point pairs, the data link table establishing module adds the sampling time and the sampling value of the first sampling point and/or the sampling time and the sampling value of the last sampling point on the first timing curve in the data link table. And if the first sampling point and/or the last sampling point on the second timing curve are not included in all the second control points of the plurality of control point pairs, the data link list establishing module adds the sampling time and the sampling value of the first sampling point and/or the sampling time and the sampling value of the last sampling point on the second timing curve in the data link list.

In this case, the phase deviation elimination module 40 may obtain the related information of each control point pair from the data link table, and shift the second timing curve based on each control point pair by using the first timing curve as a reference, so that the phases of the first timing curve and the second timing curve are consistent.

The phase deviation elimination module 40 shifts the second timing curve based on the first timing curve based on each control point pair, so that the phases of the first timing curve and the second timing curve are consistent. Here, each control point pair may be composed of a first control point on the first timing curve and a second control point on the second timing curve.

In a preferred embodiment, the phase deviation elimination module 40 may divide the first timing curve and the second timing curve based on each control point, and then shift the divided curve segments to eliminate the phase deviation between the first timing curve and the second timing curve.

The process of performing the shift for each divided curve segment to eliminate the phase deviation of the two timing curves will be described with reference to fig. 14.

Fig. 14 illustrates a block diagram of the phase deviation elimination module 40 according to an exemplary embodiment of the present invention.

As shown in fig. 14, the phase deviation elimination module 40 according to an exemplary embodiment of the present invention may include:

a first segmentation submodule 401, a second segmentation submodule 402 and a translation submodule 403.

Specifically, the first division submodule 401 divides the first timing curve into a plurality of first curve segments based on all the first control points on the first timing curve.

The second partitioning submodule 402 partitions the second timing curve into a plurality of second curve segments based on all of the second control points on the second timing curve.

The translation submodule 403 translates, for each first curve segment and a second curve segment of the plurality of second curve segments corresponding in order to each first curve segment, the second curve segment corresponding in order to the first curve segment with reference to the first curve segment so that phases of the first curve segment and the second curve segment corresponding in order to the first curve segment are consistent.

As an example, the translation submodule 403 may translate the second curve segment corresponding in order to any of the first curve segments with reference to the any of the first curve segments in such a way as to make the phases of the any of the first curve segments coincide with the second curve segment corresponding in order to the any of the first curve segments: and elastically translating the second curve segment corresponding to any one first curve segment in sequence into a time period corresponding to any one first curve segment, so that the sampling time of the first sampling point of the second curve segment corresponding to any one first curve segment in sequence is consistent with the sampling time of the first sampling point of any one first curve segment, and the sampling time of the last sampling point of the second curve segment corresponding to any one first curve segment in sequence is consistent with the sampling time of the last sampling point of any one first curve segment.

In a preferred embodiment, the translation submodule 403 may translate the second curve segment corresponding in order to any of the first curve segments with reference to any of the first curve segments in such a way that the phases of the any of the first curve segments and the second curve segment corresponding in order to the any of the first curve segments coincide: dividing a time period corresponding to any one first curve segment into N +1 equal parts, wherein N is the number of middle sampling points contained in a second curve segment corresponding to any one first curve segment in sequence; and taking the sampling value of each intermediate sampling point on the second time sequence curve corresponding to any one first curve segment in sequence before translation as the sampling value of the halving point corresponding to each intermediate sampling point on the second time sequence curve corresponding to any one first curve segment in sequence after translation.

Preferably, the apparatus for removing a phase deviation of a timing curve according to an exemplary embodiment of the present invention may further include: and an interpolation module (not shown in the figure) for obtaining the sampling value of the middle sampling point contained in any first curve segment on the translated second time sequence curve through interpolation processing based on the sampling value of the first sampling point, the sampling value of the last sampling point and the sampling value of each equal dividing point on the translated second time sequence curve corresponding to any first curve segment in sequence.

As an example, the first time-series curve may be a load time-series curve (i.e. a first load time-series curve) corresponding to an upper boundary point adjacent to a value of any one of the preset wind parameters under the condition to be solved, and the second time-series curve may be a load time-series curve (i.e. a second load time-series curve) corresponding to a lower boundary point adjacent to a value of the any one of the preset wind parameters under the condition to be solved.

Preferably, the apparatus for removing a phase deviation of a timing curve according to an exemplary embodiment of the present invention may further include: and a time sequence load determining module (not shown in the figure) which obtains a load time sequence curve corresponding to the numerical value of any preset wind parameter under the working condition to be solved through interpolation processing based on the first time sequence curve and the second time sequence curve.

For example, the time-series load determination module may pre-establish a load sample library, in which the wind parameter value combinations and the time-series loads corresponding to the wind parameter value combinations are stored, where the wind parameter value combinations in the load sample library are obtained by arranging and combining a plurality of value points corresponding to each preset wind parameter.

Based on a pre-established load sample library, the combination of different wind parameter values is used as input, time sequence loads (namely load time sequence curves) corresponding to upper and lower boundary points of the numerical value of a single preset wind parameter are found out in the load sample library, and the time sequence loads under the numerical values of the preset wind parameters under the working condition to be solved are finally obtained through continuous elimination and linear interpolation.

For example, the time-series load determining module may search a plurality of numerical points corresponding to any one of the preset wind parameters from the load sample library, and select two numerical points closest to the numerical value of any one of the preset wind parameters from the searched plurality of numerical points as two boundary numerical points adjacent to the numerical value of the preset wind parameter. Here, the two boundary value points of any one of the preset wind parameters may include an upper boundary point that is greater than the value of any one of the preset wind parameters and a lower boundary point that is less than the value of any one of the preset wind parameters. Then, the time sequence load determining module searches a first load time sequence curve corresponding to the working condition of the upper boundary point of the numerical value of any preset wind parameter and a second load time sequence curve corresponding to the working condition of the lower boundary point of the numerical value of any preset wind parameter from the load sample library. And obtaining a load time sequence curve corresponding to the numerical value of any preset wind parameter under the working condition to be solved through interpolation processing based on the first load time sequence curve and the second load time sequence curve.

There is also provided, in accordance with an exemplary embodiment of the present invention, a computing device. The computing device includes a processor and a memory. The memory is for storing a computer program. The computer program is executed by a processor, and causes the processor to execute the method for eliminating the phase deviation of the timing curve described above.

There is also provided, in accordance with an exemplary embodiment of the present invention, a computer-readable storage medium storing a computer program. The computer-readable storage medium stores a computer program that, when executed by a processor, causes the processor to execute the method of eliminating the phase deviation of the timing curve described above. The computer readable recording medium is any data storage device that can store data read by a computer system. Examples of the computer-readable recording medium include: read-only memory, random access memory, read-only optical disks, magnetic tapes, floppy disks, optical data storage devices, and carrier waves (such as data transmission through the internet via wired or wireless transmission paths).

By adopting the method and the device for eliminating the phase deviation of the time sequence curves, the phase deviation between the two time sequence curves can be effectively eliminated, and the accuracy of linear interpolation based on the two time sequence curves is improved.

In addition, with the method and the device for eliminating the phase deviation of the time sequence curve according to the exemplary embodiment of the present invention, the method is suitable for eliminating the phase errors of the time sequence curve of each component and each variable of the wind turbine generator, and in addition, the method can also be used for eliminating the phase deviation in other fields.

In addition, by adopting the method and the device for eliminating the phase deviation of the time sequence curve, disclosed by the exemplary embodiment of the invention, the linear correlation of the two time sequence curves for interpolation processing can be obviously improved, the accuracy of linear interpolation is obviously improved, and the method and the device have important significance for data mining of load information of the wind turbine generator.

In addition, the method and the device for eliminating the phase deviation of the time sequence curve, provided by the exemplary embodiment of the invention, provide an effective method for searching and filtering extreme points, and provide an effective method for translating the time sequence curve by utilizing the characteristics of the spring.

While the invention has been shown and described with reference to the preferred embodiments, it will be understood by those skilled in the art that various changes and modifications may be made therein without departing from the spirit and scope of the invention as defined by the appended claims.

Claims

1. A method for eliminating phase deviation of a timing curve, the method comprising:

determining a plurality of first extreme points on a first timing curve;

determining a plurality of second extreme points on a second timing curve;

determining a plurality of control point pairs based on the plurality of first extreme points and the plurality of second extreme points;

based on each control point pair, the second timing curve is shifted by taking the first timing curve as a reference so as to enable the phases of the first timing curve and the second timing curve to be consistent.

2. The cancellation method of claim 1, wherein each control point pair is composed of a first control point on a first timing curve and a second control point on a second timing curve, the first control point being one of the plurality of first extreme points, the second control point being a second extreme point of the plurality of second extreme points that forms a control point pair with the one first extreme point;

wherein shifting the second timing curve based on the first timing curve based on each control point pair to make the phases of the first timing curve and the second timing curve consistent comprises:

dividing the first timing curve into a plurality of first curve segments based on all first control points on the first timing curve,

dividing the second timing curve into a plurality of second curve segments based on all second control points on the second timing curve,

for each first curve segment and a second curve segment of the plurality of second curve segments corresponding in order to each first curve segment, translating the second curve segment corresponding in order to the first curve segment with reference to the first curve segment to make the phases of the first curve segment and the second curve segment corresponding in order to the first curve segment coincide.

3. The cancellation method of claim 2, wherein a second curve segment corresponding in order to any one of the first curve segments is translated with reference to the any one of the first curve segments in such a manner that phases of the any one of the first curve segments and the second curve segment corresponding in order to the any one of the first curve segments coincide:

and elastically translating the second curve segment corresponding to any one first curve segment in sequence into a time period corresponding to any one first curve segment, so that the sampling time of the first sampling point of the second curve segment corresponding to any one first curve segment in sequence is consistent with the sampling time of the first sampling point of any one first curve segment, and the sampling time of the last sampling point of the second curve segment corresponding to any one first curve segment in sequence is consistent with the sampling time of the last sampling point of any one first curve segment.

4. The cancellation method of claim 3, wherein a second curve segment corresponding in order to any one of the first curve segments is translated with reference to the any one of the first curve segments in such a manner that phases of the any one of the first curve segments and the second curve segment corresponding in order to the any one of the first curve segments coincide:

dividing a time period corresponding to any one first curve segment into N +1 equal parts, wherein N is the number of intermediate sampling points contained in a second curve segment corresponding to any one first curve segment in sequence;

and taking the sampling value of each intermediate sampling point on the second time sequence curve corresponding to any one first curve segment in sequence before translation as the sampling value of the halving point corresponding to each intermediate sampling point on the second time sequence curve corresponding to any one first curve segment in sequence after translation.

5. The cancellation method of claim 4, wherein the cancellation method further comprises:

and obtaining the sampling value of the middle sampling point contained in any first curve segment by the second time sequence curve after translation through interpolation processing based on the sampling value of the first sampling point, the sampling value of the last sampling point and the sampling values of all the equant points on the second time sequence curve corresponding to any first curve segment in sequence after translation.

6. The cancellation method according to any one of claims 1 through 5, wherein the plurality of first extreme points include a plurality of first peak points and a plurality of first valley points, the plurality of second extreme points include a plurality of second peak points and a plurality of second valley points, the plurality of control point pairs include a plurality of peak point pairs and a plurality of valley point pairs,

wherein determining a plurality of control point pairs based on the plurality of first extreme points and the plurality of second extreme points comprises:

generating the plurality of pairs of peak points based on the plurality of first peak points and the plurality of second peak points,

generating the plurality of valley point pairs based on the plurality of first valley points and the plurality of second valley points.

7. The cancellation method of claim 6, wherein the step of determining a plurality of first extreme points on the first timing curve comprises:

selecting a plurality of first original peak points as the plurality of first peak points at predetermined intervals from all the first original peak points on the first timing curve,

selecting a plurality of first original valley points as the plurality of first valley points at the predetermined interval from all the first original valley points on the first timing curve,

and/or the step of determining a plurality of second extreme points on the second timing curve comprises:

selecting a plurality of second original peak points as the plurality of second peak points at the predetermined interval from all the second original peak points on the second timing curve,

and selecting a plurality of second original valley points as the plurality of second valley points at the preset interval from all the second original valley points on the second time sequence curve.

8. The cancellation method of claim 7, wherein each first original peak point and each first original valley point is determined by:

determining all sampling points on the first time sequence curve;

for any sampling point, if the sampling value of the sampling point is greater than the sampling value of the last sampling point and greater than the sampling value of the next sampling point, determining the sampling point as a first original peak point, and if the sampling value of the sampling point is less than the sampling value of the last sampling point and less than the sampling value of the next sampling point, determining the sampling point as a first original valley point,

and/or each second original peak point and each second original valley point is determined by:

determining all sampling points on the second timing curve;

and for any sampling point, if the sampling value of the sampling point is greater than the sampling value of the last sampling point and greater than the sampling value of the next sampling point, determining the sampling point as a second original peak point, and if the sampling value of the sampling point is less than the sampling value of the last sampling point and less than the sampling value of the next sampling point, determining the sampling point as a second original valley point.

9. The cancellation method of claim 6, wherein each pair of peak points is determined by:

searching whether a second peak point with a time interval of the sampling time of the first peak point being less than a preset time interval exists on a second timing curve aiming at each first peak point on the first timing curve, if the second peak point with the time interval of the sampling time of the first peak point being less than the preset time interval exists, determining the first peak point and the searched second peak point as a peak point pair,

and/or, each valley point pair is determined by:

10. The cancellation method of claim 2, wherein the cancellation method further comprises:

and storing the related information of the plurality of control point pairs into a data linked list, wherein the related information comprises the sampling time and the sampling value of each first control point in the plurality of control point pairs and the sampling time and the sampling value of each second control point.

11. The cancellation method of claim 10, wherein the cancellation method further comprises:

adding the sampling time and the sampling value of the last sampling point on the first timing curve to the data link list if the last sampling point on the first timing curve is not included in all the first control points of the plurality of control point pairs,

if the first sampling point and/or the last sampling point on the second timing curve are not included in all the second control points of the plurality of control point pairs, adding the sampling time and the sampling value of the last sampling point on the second timing curve in the data link list,

and translating the second time sequence curve by taking the first time sequence curve as a reference based on each control point pair so as to enable the phases of the first time sequence curve and the second time sequence curve to be consistent.

12. The elimination method of claim 1, wherein the first time series curve is a load time series curve corresponding to an upper boundary point adjacent to a value of any one of the predetermined wind parameters under the condition to be solved, the second time series curve is a load time series curve corresponding to a lower boundary point adjacent to a value of said any one of the predetermined wind parameters under the condition to be solved,

wherein the elimination method further comprises:

and obtaining a load time sequence curve corresponding to the numerical value of any preset wind parameter under the working condition to be solved through interpolation processing based on the first time sequence curve and the second time sequence curve.

13. An apparatus for canceling a phase deviation of a timing curve, the apparatus comprising:

the first extreme value determining module is used for determining a plurality of first extreme values on the first time sequence curve;

the second extreme value determining module is used for determining a plurality of second extreme value points on the second time sequence curve;

a control point pair determination module that determines a plurality of control point pairs based on the plurality of first extreme points and the plurality of second extreme points;

and the phase deviation elimination module is used for shifting the second time sequence curve by taking the first time sequence curve as a reference on the basis of each control point pair so as to enable the phases of the first time sequence curve and the second time sequence curve to be consistent.

14. The removing apparatus according to claim 13, wherein each control point pair is composed of a first control point on a first timing curve and a second control point on a second timing curve, the first control point being one of the plurality of first extreme points, the second control point being a second extreme point of the plurality of second extreme points forming a control point pair with the one first extreme point;

wherein, phase deviation eliminates the module and includes:

a first division submodule which divides the first timing curve into a plurality of first curve segments based on all the first control points on the first timing curve,

a second segmentation sub-module that segments the second timing curve into a plurality of second curve segments based on all second control points on the second timing curve,

a translation submodule that translates, for each first curve segment and a second curve segment of the plurality of second curve segments corresponding in order to each first curve segment, the second curve segment corresponding in order to the first curve segment with reference to the first curve segment so that phases of the first curve segment and the second curve segment corresponding in order to the first curve segment coincide.

15. The cancellation apparatus of claim 14, wherein the translation sub-module translates a second curve segment corresponding in order to any one of the first curve segments with reference to the any one of the first curve segments to bring the any one of the first curve segments into phase agreement with the second curve segment corresponding in order to the any one of the first curve segments by:

16. The cancellation apparatus of claim 15, wherein the translation sub-module translates a second curve segment corresponding in order to any one of the first curve segments with reference to the any one of the first curve segments to bring the any one of the first curve segments into phase agreement with the second curve segment corresponding in order to the any one of the first curve segments by:

dividing a time period corresponding to any one first curve segment into N +1 equal parts, wherein N is the number of middle sampling points contained in a second curve segment corresponding to any one first curve segment in sequence;

17. The cancellation apparatus of claim 16, wherein the cancellation apparatus further comprises: and the interpolation module is used for obtaining the sampling value of the middle sampling point contained in any first curve segment by the translated second time sequence curve through interpolation processing based on the sampling value of the first sampling point, the sampling value of the last sampling point and the sampling values of all the equant points on the translated second time sequence curve corresponding to any first curve segment in sequence.

18. The cancellation apparatus of any one of claims 13 through 17, wherein the plurality of first extreme points include a plurality of first peak points and a plurality of first valley points, the plurality of second extreme points include a plurality of second peak points and a plurality of second valley points, the plurality of control point pairs include a plurality of peak point pairs and a plurality of valley point pairs,

wherein the control-point pair determination module generates the plurality of peak-point pairs based on the plurality of first peak points and the plurality of second peak points, and generates the plurality of valley-point pairs based on the plurality of first valley points and the plurality of second valley points.

19. The removing apparatus of claim 18, wherein the first extreme value determining module selects a plurality of first original peak points as the plurality of first peak points at predetermined intervals from all first original peak points on the first timing curve, selects a plurality of first original valley points as the plurality of first valley points at the predetermined intervals from all first original valley points on the first timing curve,

and/or the second polarity determining module selects a plurality of second original peak points as the plurality of second peak points at the preset interval from all the second original peak points on the second time sequence curve, and selects a plurality of second original valley points as the plurality of second valley points at the preset interval from all the second original valley points on the second time sequence curve.

20. The cancellation apparatus of claim 19, wherein the first extremum determining module determines each of the first raw peak points and each of the first raw valley points by:

determining all sampling points on the first time sequence curve;

and/or the second extreme determining module determines each second original peak point and each second original valley point by:

determining all sampling points on the second time sequence curve;

21. The elimination apparatus of claim 18, wherein the control point pair determination module determines each pair of peak points by:

and/or, the control point pair determination module determines each valley point pair by:

22. The cancellation apparatus of claim 14, wherein the cancellation apparatus further comprises:

and the storage module stores the related information of the plurality of control point pairs into a data link list, wherein the related information comprises the sampling time and the sampling value of each first control point in the plurality of control point pairs and the sampling time and the sampling value of each second control point.

23. The elimination apparatus of claim 22, wherein the data link list establishment module adds the sample time and the sample value of the last sample point on the first timing curve to the data link list if the last sample point on the first timing curve is not included in all of the first control points of the plurality of control point pairs,

the data link table establishing module adds a sampling time and a sampling value of a last sampling point on the second timing curve to the data link table if the last sampling point on the second timing curve is not included in all the second control points of the plurality of control point pairs,

the phase deviation elimination module acquires the related information of the plurality of control point pairs from the data linked list, and shifts the second time sequence curve based on each control point pair by taking the first time sequence curve as a reference so as to enable the phases of the first time sequence curve and the second time sequence curve to be consistent.

24. The elimination apparatus of claim 13, wherein the first timing curve is a load timing curve corresponding to an upper boundary point adjacent to a value of any one of the predetermined wind parameters under the to-be-solved condition, the second timing curve is a load timing curve corresponding to a lower boundary point adjacent to a value of the any one of the predetermined wind parameters under the to-be-solved condition,

wherein the eliminating device further comprises: and the time sequence load determining module is used for obtaining a load time sequence curve corresponding to the numerical value of any preset wind parameter under the working condition to be solved through interpolation processing based on the first time sequence curve and the second time sequence curve.

25. A computer-readable storage medium, in which a computer program is stored which, when being executed by a processor, carries out a method of canceling a phase deviation of a timing curve according to any one of claims 1 to 12.

26. A computing device, the computing device comprising:

a processor;

a memory storing a computer program which, when executed by the processor, implements the method of canceling a phase deviation of a timing curve according to any one of claims 1 to 12.