CN106997407A - Wind-resources scene reduction method based on trend fitting - Google Patents
Wind-resources scene reduction method based on trend fitting Download PDFInfo
- Publication number
- CN106997407A CN106997407A CN201611248741.7A CN201611248741A CN106997407A CN 106997407 A CN106997407 A CN 106997407A CN 201611248741 A CN201611248741 A CN 201611248741A CN 106997407 A CN106997407 A CN 106997407A
- Authority
- CN
- China
- Prior art keywords
- wind
- scene
- resources
- trend
- control point
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 47
- 238000000926 separation method Methods 0.000 claims description 9
- 230000015572 biosynthetic process Effects 0.000 claims description 5
- 238000003786 synthesis reaction Methods 0.000 claims description 5
- 230000011218 segmentation Effects 0.000 claims description 3
- 230000000694 effects Effects 0.000 abstract description 5
- 238000000605 extraction Methods 0.000 abstract description 4
- 238000010586 diagram Methods 0.000 description 6
- 230000018109 developmental process Effects 0.000 description 4
- 238000006467 substitution reaction Methods 0.000 description 2
- 230000001360 synchronised effect Effects 0.000 description 2
- 230000007812 deficiency Effects 0.000 description 1
- 238000011946 reduction process Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Wind Motors (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a kind of wind-resources scene reduction method based on trend fitting, it is characterised in that comprises the following steps：Parsing wind-resources scene graph simultaneously carries out trend fitting, afterwards trend fitting is carried out for each wind-resources scene graph in wind-resources scene process cluster, obtain corresponding global trend sequence, scene reduction parameter assignment is carried out afterwards, then calculates the distance between each global trend sequence in the preceding wind-resources scene process cluster of epicycle scene reduction；Finally reduction scene and more new scene probability are carried out until being met scene set and its probability of scene reduction number requirement.Refine the trend of wind-resources scene graph to judge the similitude between wind-resources scene by automation, reduce similar scene, and then reach that scene reduces purpose, the randomness of wind-resources is avoided to disturb, be conducive to the judgement of wind-resources scene similitude, strengthen wind-resources typical scene extraction effect, utilize significant for wind energy development.
Description
Technical field
The invention belongs to wind-resources specificity analysis technical field, more particularly to a kind of wind-resources based on trend fitting
Scape reduction method.
Background technology
The research of wind-resources scene reduction is mainly used in the extraction of wind-resources typical scene, reduces the work of wind-resources analysis
Measure, improve analysis efficiency, and then technical support is provided for wind energy development utilization.Improving the effect of wind-resources scene reduction has
Help strengthen the recognition capability of wind-resources typical scene, lift wind-resources analytical technology level.Current wind-resources both domestic and external
Scape reduction technology is directly handled wind-resources time series, is similar by time series scene Recognition closer to the distance
Scene, and delete similar scene to reach that scene reduces purpose.Similar scene means the similar scene of trend.Wind-resources are same
When there is tendency and randomness.The existing technology in the field not yet considers the randomness of wind-resources to wind-resources time series
The interference that distance is calculated, have impact on the recognition effect of wind-resources similar scene.
The content of the invention
The technical problem to be solved in the present invention is to overcome the deficiencies in the prior art there is provided a kind of based on trend fitting
Wind-resources scene reduction method, the trend for wind-resources carries out scene reduction, it is to avoid the randomness interference of wind-resources, improves
Wind-resources typical scene extraction effect, preferably utilizes for wind energy development and provides technical support.
In order to solve the above technical problems, the present invention is adopted the following technical scheme that：
A kind of wind-resources scene reduction method based on trend fitting, it is characterised in that comprise the following steps：
Step 1：According to wind speed and the relation of time in rectangular coordinate system, wind-resources scene graph, setting solution are parsed
Wind-resources scene graph is set after analysis, and by N number of control point, connection is constituted, and N >=2, NN wind-resources scene graph constitutes wind
Resource scene process cluster, NN >=3；
Step 2：If N>3, wind-resources scene hydrograph separation is carried out, afterwards to place wind-resources scene process after segmentation
The specific subinterval of line carries out trend fitting respectively, obtains the respective corresponding trend sequence as local trend；Afterwards again
The global trend of wind-resources scene graph where local trend is fitted to；If N≤3, without carrying out wind-resources scene process
Line is split, and directly carries out trend fitting to wind-resources scene graph, the global trend of wind-resources scene graph where obtaining
Sequence；
Step 3：Enter for each wind-resources scene graph in wind-resources scene process cluster in the way of step 2
Row trend fitting, obtains the global trend sequence of corresponding all wind-resources scenes, shared NN global trend sequence of setting；
Step 4：Scene reduces parameter assignment step, sets the wind-resources scene number for intending reduction as MM, i.e., scene is reduced
The number of wind-resources scene set Scene afterwards is NN-MM；
Step 5, the wind money of each wind-resources scene graph in wind-resources scene process cluster before the reduction of epicycle scene is calculated
The distance between source overall situation trend sequence；
Step 6, reduction scene and more new scene probability：
The wind-resources scene of minimum probability is determined according to the probability right distance between the global trend sequence of wind-resources, will
The wind-resources scene of the minimum probability is deleted from the set of probability right distance and updates the general of not deleted wind-resources scene
Rate, obtains the scene set after epicycle reduction scene and its corresponding probability；Scene reduction is performed successively, is such as performed scene and is replaced
The number of times mm changed<MM, then go to step 5 and start new round reduction, otherwise end step 5 is met to the circulation of step 6
The scene set of scene reduction number requirement and its probability.
Further, step 1 parsing wind-resources scene graph mode specific as follows is carried out：
The wind-resources scene graph be in rectangular coordinate system, according to the process of wind speed, using time t as abscissa,
Obtained using wind speed w as ordinate；Resolving includes resolving to wind-resources scene graph by some control point connection groups
Into if having N number of control point, numbering is followed successively by 1,2 from left to right ..., N, and i-th of control point coordinates is designated as (ti,wi), i=
1,2,...,N；Wind-resources scene graph is designated as { t, w }；If having NN wind-resources scene graph, wind-resources are constituted
Scene process cluster, is designated as { TT, WW }, and wherein wind-resources scene graph numbering is designated as ii=1,2 ..., NN.
Further, step 2 is distinguished by following three step and carried out：
Step 2.1, wind-resources scene hydrograph separation：
If N≤3, without carrying out wind-resources scene hydrograph separation, step 2.3 is directly entered；
If N>3, by (sg-1) × m+1 control point to (sg+1) × m+1 of wind-resources scene graph { t, w }
Individual control point is designated as the sg subinterval, and wherein m is integer, and during N=4, m takes 1, N>Span is when 4Sg=1,2 ..., SG, int (*) represent to round " * "；As m=1, SG=N-3；Work as m>When 1, ifThenOtherwiseExtremely
This, can obtain SG subinterval；The SG × m+1 control point to the n-th control point of { t, w } is designated as the SG+1 sub-district
Between；So far, can obtain between two subintervals adjacent in SG+1 subinterval, and gained subinterval has m+1 point overlapping,
I.e. above the rear m+1 point in subinterval and the preceding m+1 point in subinterval below are overlapping；There is 2m+1 control in preceding SG subinterval
Processed, it controls point coordinates to be designated as：The SG+1 subinterval contain N-SG ×
M control point, it controls point coordinates to be designated as
Step 2.2, wind-resources local trend is fitted：
If N>3, SG+1 subinterval obtained by step 2.1 is subjected to trend fitting respectively, each subinterval is obtained relative
The trend sequence answered；The control point coordinates of the trend sequence in preceding SG subinterval is designated as：Sg=1,2 ..., SG；The trend sequence in the SG+1 subinterval
Control point coordinates be designated as：For ease of distinguishing,
Title works as N>Trend sequence is obtained for local trend in this step 2.2 when 3；
Step 2.3, the global trend synthesis of wind-resources：
For ease of distinguishing, the trend that this step 2.3 is obtained is called global trend；
If N≤3, trend fitting directly is carried out to wind-resources scene graph { t, w }, its global trend sequence { t is obtained
,glowf, subsequently into step 3；
If N>When 3, global trend must be synthesized according to local trend；Specifically building-up process is：Take the 1st subinterval part
The preceding m control points of trendIt is used as preceding m of the global trend sequence of wind-resources
Control point, is designated asTake the rear N- (SG+ of the SG+1 subinterval local trend
1) × m-1 control pointIt is global as wind-resources
(SG+1) × m+2 control point to the n-th control point of trend sequence, i.e.,The SG+1 sub-district obtained by step 2.1
Between there is SG lap, according to being weighted place to the influence degree of its lap fitted trend between adjacent subarea
Reason, the local trend of the part is weighted as the following formula, the weighting local trend sequence of lap is obtained：
Wherein：The sg lap refers to the rear m+1 control point and the sg+1 subinterval in the sg subinterval
Preceding m+1 control point is overlapping；For j-th of control point in the weighting local trend sequence of the sg lap
Ordinate, j=1,2 ..., m+1, λ1、λ2For weight coefficient, For the sg subinterval
+ m control point ordinates of jth of local trend sequence,For j-th of control in the fitting sequence in the sg+1 subinterval
System point ordinate；The control point coordinates of the weighting local trend sequence of the sg lap is designated as：SG can be obtained after processing is weighted to SG lap
Individual weighting local trend sequence, and last control point coordinates and the latter weighting of previous weighting local trend sequence
First control point coordinates of local trend sequence is identical；Carry out deduplication operation；The specific of duplicate removal is processed as previous weighting
Last control point coordinates of local trend sequence retains, and removes first control of its latter weighting local trend sequence
Point coordinates processed；After such duplicate removal, the 1st weighting local trend sequence includes m+1 control point, i.e.,Remaining weighting local trend sequence includes m control point, i.e.,SG weighting local trend sequence after duplicate removal is joined end to end, it is taken as wind
The m+1 control point to (SG+1) × m+1 control point of resource overall situation trend sequence, i.e.,
So far, the global trend sequence of the wind-resources comprising N number of control point can be obtained, be also designated as t,glowf}。
Further, step 5 is carried out as follows：
Wind-resources scene set before epicycle is reduced is designated as { TT, WW }bef, the number of wind-resources scene before epicycle reduction
KK is designated as, the sequence number of wind-resources scene is designated as kk, i.e. kk=1, and 2 ..., KK, the probability of correspondence wind-resources scene are designated as pkk；When
When performing this step the 1st time, KK=NN, it is equal to be regarded as the probability that NN wind-resources scene occur, and is designated as pkk=1/NN；According to
{TT,WW}befIn the corresponding global trend sequence of each wind-resources scene graph calculate the global trend sequence of wind-resources two-by-two
Between probability right distance, be designated as PDkk, formula is as follows：
WhereinKk ≠ jj, kk=1,2 ..., KK, jj=1,2 ..., KK；
For k scene of kth, the set containing KK-1 distance can be obtained, { D is designated askk-jj, by { Dkk-jjIn it is minimum
It is worth corresponding scene sequence number jj labeled as JLkk, can be obtained containing KK JL for KK scenekkSet, be designated as { JLkk}；It is right
Wind-resources scene set { TT, WW } before epicycle reductionbef, the set containing KK probability right distance can be obtained, is designated as
{PDkk, by { PDkkIn the corresponding scene sequence number kk of minimum value be labeled as PDmin.
Further, step 6 is carried out as follows：
By serial number PDmin scene from { TT, WW }befIt is middle to delete；In { JLkkIn find out PDmin element
JLPDmin, then by the Probability p of serial number PDmin scenePDminAdd to serial number JLPDminScene probabilityI.e.So far, it can obtain the scene set { TT, WW } after epicycle reduction sceneaftAnd its it is corresponding general
Rate；The number of times for performing this step is designated as mm, if mm<MM, then go to step 5 and start new round reduction, otherwise end step 5
To the circulation of step 6, scene set and its probability of scene reduction number requirement are met.
Wind-resources scene reduction technology scheme provided by the present invention based on trend fitting, wind is refined by automating
The trend of resource scene graph is to judge the similitude between wind-resources scene there is provided new determination methods, as a result simply
Understand, be easy to implement easy.In analysis of Wind Energy Resource application, the scene number that wind-resources scene process and needs are reduced
It is used as input, you can the similarity degree between automatic decision scene, reduces similar scene, and then reach that scene reduces purpose.It is right
Than prior art, propose to carry out the similar identification of scene by distinguishing rule of the trend of wind-resources first, so as to avoid wind-resources
Randomness is disturbed, and is the important innovations of the art, is conducive to the judgement of wind-resources scene similitude, strengthens wind-resources allusion quotation
Type scene extraction effect, for wind energy development using significant, with important popularizing value.
Brief description of the drawings
Fig. 1 is that schematic diagram is split in the subinterval for the wind-resources scene graph implemented according to the present invention.Fig. 2 is the present invention
The scene synthesis schematic diagram of embodiment.
Fig. 3 be the embodiment of the present invention reduction before initial scene set schematic diagram.
Fig. 4 is to use existing scene reduction technology --- synchronous back substitution "flop-out" method carries out the process of scene reduction (in figure
Scene number is the scene number after reduction).
Fig. 5 be the use technical solution of the present invention of the embodiment of the present invention scene reduction process (figure Scene number for contracting
Scene number after subtracting).
Embodiment
In order that the purpose of the embodiment of the present invention, technical scheme, advantage become apparent from, below in conjunction with present invention implementation
Example and accompanying drawing introduce technical scheme.
The present invention provides a kind of wind-resources scene reduction method based on trend fitting, comprises the following steps：
Step 1, wind-resources scene graph is parsed
The wind-resources scene graph be in rectangular coordinate system, according to the process of wind speed, using time t as abscissa,
Obtained using wind speed w as ordinate；Resolving includes resolving to wind-resources scene graph by some control point connection groups
Into if having N number of control point, numbering is followed successively by 1,2 from left to right ..., N, and i-th of control point coordinates is designated as (ti,wi), i=
1,2,...,N；Wind-resources scene graph is designated as { t, w }.If having NN wind-resources scene graph, wind-resources are constituted
Scene process cluster, is designated as { TT, WW }, and wherein wind-resources scene graph numbering is designated as ii=1,2 ..., NN.
Step 2.1, wind-resources scene hydrograph separation
If N≤3, without carrying out wind-resources scene hydrograph separation, step 2.3 is directly entered.
If N>3, by (sg-1) × m+1 control point to (sg+1) × m+1 of wind-resources scene graph { t, w }
Individual control point is designated as the sg subinterval, and wherein m is integer, and its span isSg=1,2 ..., SG,
Int (*) represents to round " * ".As m=1, SG=N-3；Work as m>When 1, ifThenOtherwiseSo far, it can obtain SG subinterval.By { t, w }
The SG × m+1 control point to n-th control point be designated as the SG+1 subinterval.So far, SG+1 subinterval is can obtain,
And have m+1 point overlapping between adjacent two subintervals in gained subinterval, i.e., above the rear m+1 point in subinterval with after
The interval preceding m+1 point of face is overlapping.There is 2m+1 control point in preceding SG subinterval, and it controls point coordinates to be designated as：N-SG × m control point is contained in the SG+1 subinterval, and its control point is sat
It is labeled asSubinterval segmentation schematic diagram can be found in Fig. 1, son
Interval 1 includes the 1st control point to the 2m+1 control point, and coordinate is designated asSub-district
Between 2 include the m+1 control point to the 3m+1 control point, coordinate is designated as
Subinterval sg includes (sg-1) × m+1 control point to (sg+1) × m+1 control point, and coordinate is designated asSubinterval SG includes (SG-1) × m+1 control point to (SG+
1) × m+1 control point, coordinate is designated asSubinterval SG+1 includes SG × m
+ 1 control point to n-th control point, coordinate is designated asThis step
M specific value is set by those skilled in the art in rapid.
Step 2.2, wind-resources local trend is fitted
If N>3, SG+1 subinterval obtained by step 2.1 is subjected to trend fitting respectively, its is obtained and corresponding becomes
Gesture sequence；The control point coordinates of the trend sequence in preceding SG subinterval is designated as：Sg=1,2 ..., SG；The trend sequence in the SG+1 subinterval
Control point coordinates be designated as：For ease of distinguishing,
Title works as N>Trend sequence is obtained for local trend in this step when 3.The trend fitting of this step takes polynomial fitting method,
Fitting exponent number is set by those skilled in the art.
Step 2.3, the global trend synthesis of wind-resources
For ease of distinguishing, the trend that this step is obtained is called global trend.If N≤3, directly to wind-resources scene process
Line { t, w } carries out trend fitting, obtain its global trend sequence t,glowf, subsequently into step 3.If N>, must basis when 3
The global trend of local trend synthesis.Specifically building-up process is：Take the preceding m control points of the 1st subinterval local trendAs the preceding m control point of the global trend sequence of wind-resources, it is designated asTake rear N- (SG+1) × m-1 of the SG+1 subinterval local trend
Individual control pointIt is used as wind
(SG+1) × m+2 control point to the n-th control point of resource overall situation trend sequence, i.e.,The SG+1 sub-district obtained by step 2.1
Between there is SG lap, the local trend of the part is weighted as the following formula, the weighting of lap is obtained
Local trend sequence.
Wherein：The sg lap refers to the rear m+1 control point and the sg+1 subinterval in the sg subinterval
Preceding m+1 control point is overlapping；For j-th of control point in the weighting local trend sequence of the sg lap
Ordinate, j=1,2 ..., m+1, λ1、λ2For weight coefficient, For the sg subinterval
+ m control point ordinates of jth of local trend sequence,For j-th of control in the fitting sequence in the sg+1 subinterval
System point ordinate.This formula is act as according to being added between adjacent subarea to the influence degree of its lap fitted trend
Power processing.
The control point coordinates of the weighting local trend sequence of the sg lap is designated as：SG can be obtained after processing is weighted to SG lap
Individual weighting local trend sequence, and last control point coordinates and the latter weighting of previous weighting local trend sequence
First control point coordinates of local trend sequence is identical.Carry out deduplication operation.The specific of duplicate removal is processed as previous weighting
Last control point coordinates of local trend sequence retains, and removes first control of its latter weighting local trend sequence
Point coordinates processed.After such duplicate removal, the 1st weighting local trend sequence includes m+1 control point, i.e.,Remaining weighting local trend sequence includes m control point, i.e.,SG weighting local trend sequence after duplicate removal is joined end to end, it is taken as wind
The m+1 control point to (SG+1) × m+1 control point of resource overall situation trend sequence, i.e.,
So far, the global trend sequence of the wind-resources comprising N number of control point can be obtained, be also designated as t,glowf}.Step 3, institute
There is the global trend sequence of wind-resources scene
According to step 2.1 to step 2.3, for each wind-resources scene graph in wind-resources scene process cluster
The global trend sequence of corresponding wind-resources is obtained, be designated as t,glowf}ii, common NN.
Step 4, scene reduction parameter assignment
The wind-resources scene number for intending reduction is set as MM, i.e., scene reduce after wind-resources scene set Scenes
Number is NN-MM.
Step 5, the probability right distance between the global trend sequence of wind-resources before the reduction of epicycle scene is calculated
Wind-resources scene set before epicycle is reduced is designated as { TT, WW }bef, the number of wind-resources scene before epicycle reduction
KK is designated as, the sequence number of wind-resources scene is designated as kk, i.e. kk=1, and 2 ..., KK, the probability of correspondence wind-resources scene are designated as pkk；When
When performing this step the 1st time, KK=NN, it is equal to be regarded as the probability that NN wind-resources scene occur, and is designated as pkk=1/NN；According to
{TT,WW}befIn the corresponding global trend sequence of each wind-resources scene graph calculate the global trend sequence of wind-resources two-by-two
Between probability right distance, be designated as PDkk, formula is as follows：
WhereinKk ≠ jj, kk=1,2 ..., KK, jj=1,2 ..., KK；
For k scene of kth, the set containing KK-1 distance can be obtained, { D is designated askk-jj, by { Dkk-jjIn it is minimum
It is worth corresponding scene sequence number jj labeled as JLkk, can be obtained containing KK JL for KK scenekkSet, be designated as { JLkk}；It is right
Wind-resources scene set { TT, WW } before epicycle reductionbef, the set containing KK probability right distance can be obtained, is designated as
{PDkk, by { PDkkIn the corresponding scene sequence number kk of minimum value be labeled as PDmin.
Step 6：Reduce scene and more new scene probability
By serial number PDmin scene from { TT, WW }befIt is middle to delete；In { JLkkIn find out PDmin element
JLPDmin, then by the Probability p of serial number PDmin scenePDminAdd to serial number JLPDminScene probabilityI.e.So far, it can obtain the scene set { TT, WW } after epicycle reduction sceneaftAnd its it is corresponding general
Rate；The number of times for performing this step is designated as mm, if mm<MM, then go to step 5 and start new round reduction, otherwise end step 5
To the circulation of step 6, scene set and its probability of scene reduction number requirement are met.
Fig. 2 is the generation schematic diagram of scene process cluster, and obedience standard normal is superimposed respectively by sinusoidal trend and ladder trend
The random process line of distribution synthesizes scene graph 1 and scene graph 2 with correspondence.By foregoing two kinds of trend respectively with random
Generation and the 200 graphs superposition for obeying standardized normal distribution, can obtain the initial fields containing 400 scene graphs
Scape set (initial scene set before reduction), as shown in Figure 3.
Fig. 4 is to use existing scene reduction technology --- synchronous back substitution "flop-out" method carries out the process of scene reduction.Fig. 5
To carry out the process of scene reduction using technical solution of the present invention.Comparison diagram 4 and Fig. 5 understand that existing scene reduction technology exists
Sine and two kinds of different trend of ladder are not identified in scene reduction, and the present invention puies forward technical scheme and then have identified
Sine and two kinds of different trend of ladder, demonstrate the validity that the present invention puies forward technical scheme.
Automatic running can be realized using computer software technology during present invention specific implementation.
By embodiment achievement, the present invention puies forward technical scheme and have identified different trend, illustrates the present invention
Validity.Understand that the present invention can automatically and efficiently extract the trend of wind-resources scene and based on this progress scene reduction, be
Wind energy development utilizes and provides decision support.
Present invention is mainly applied to the reduction of wind-resources scene, in analysis of Wind Energy Resource application, by wind-resources scene process
Input is used as with the scene number for needing to reduce, you can the similarity degree between automatic decision scene, reduce similar scene, enter
And reach scene and reduce purpose.Compared with existing correlation technique, innovation of the invention be by trend with judge scene it
Between similitude.In consideration of it, be applied to of the invention from existing technology by the random of different trend and same distribution simultaneously
It is superimposed in the reduction of the scene of obtained scene process cluster, can be used to verify the reasonability of technical solution of the present invention.
It is emphasized that embodiment of the present invention is illustrative, rather than it is limited, therefore the present invention
The embodiment described in embodiment is not limited to, it is every by those skilled in the art's technique according to the invention scheme
The other embodiment drawn, also belongs to the scope of protection of the invention.
Claims (5)
1. a kind of wind-resources scene reduction method based on trend fitting, it is characterised in that comprise the following steps：
Step 1：According to wind speed and the relation of time in rectangular coordinate system, wind-resources scene graph is parsed, after setting parsing
By N number of control point, connection is constituted setting wind-resources scene graph, N >=2, NN wind-resources scene graph composition wind-resources
Scape process cluster, NN >=3；
Step 2：If N>3, wind-resources scene hydrograph separation is carried out, afterwards to the spy of place wind-resources scene graph after segmentation
Stator interval carries out trend fitting respectively, obtains the respective corresponding trend sequence as local trend；Afterwards again will be local
Global trend of the trend fitting into place wind-resources scene graph；If N≤3, without carrying out wind-resources scene hydrograph separation,
Trend fitting directly is carried out to wind-resources scene graph, the global trend sequence of wind-resources scene graph where obtaining；
Step 3：Trend is carried out in the way of step 2 for each wind-resources scene graph in wind-resources scene process cluster
Fitting, obtains the global trend sequence of corresponding all wind-resources scenes, shared NN global trend sequence of setting；
Step 4：Scene reduces parameter assignment step, sets the wind-resources scene number for intending reduction as MM, i.e., after scene reduction
The number of wind-resources scene set Scene is NN-MM；
Step 5, the wind-resources for calculating each wind-resources scene graph in wind-resources scene process cluster before the reduction of epicycle scene are complete
The distance between office's trend sequence；
Step 6, reduction scene and more new scene probability：
The wind-resources scene of minimum probability is determined according to the probability right distance between the global trend sequence of wind-resources, by the minimum
The wind-resources scene of probability is deleted from the set of probability right distance and updates the probability of not deleted wind-resources scene, is obtained
Scene set and its corresponding probability after epicycle reduction scene；Scene reduction is performed successively, such as performs the number of times that scene is replaced
mm<MM, then go to step 5 and start new round reduction, and otherwise end step 5 is met scene reduction to the circulation of step 6
Scene set and its probability that number is required.
2. the wind-resources scene reduction method according to claim 1 based on trend fitting, it is characterised in that step 1 is parsed
Wind-resources scene graph step specific as follows is carried out：
The wind-resources scene graph is in rectangular coordinate system, according to the process of wind speed, using time t as abscissa, with wind
Fast w obtains for ordinate；Resolving includes resolving to wind-resources scene graph to be connected by some control points constituting, if altogether
There is N number of control point, numbering is followed successively by 1,2 from left to right ..., N, i-th of control point coordinates is designated as (ti,wi), i=1,2 ...,
N；Wind-resources scene graph is designated as { t, w }；If having NN wind-resources scene graph, wind-resources scene process is constituted
Cluster, is designated as { TT, WW }, and wherein wind-resources scene graph numbering is designated as ii=1,2 ..., NN.
3. the wind-resources scene reduction method according to claim 2 based on trend fitting, it is characterised in that step 2 is by such as
Lower three steps difference is carried out：
Step 2.1, wind-resources scene hydrograph separation：
If N≤3, without carrying out wind-resources scene hydrograph separation, step 2.3 is directly entered；
If N>3, by (sg-1) × m+1 control point to (sg+1) of wind-resources scene graph { t, w } × m+1 control
Point is designated as the sg subinterval, and wherein m is integer, and during N=4, m takes 1, N>Span is when 4Sg=1,
2 ..., SG, int (*) represent to round " * "；As m=1, S G=N-3；Work as m>When 1, ifThenOtherwiseSo far,
It can obtain SG subinterval；The SG × m+1 control point to the n-th control point of { t, w } is designated as the SG+1 subinterval；Extremely
This, can obtain between two subintervals adjacent in SG+1 subinterval, and gained subinterval has m+1 point overlapping, i.e., before
The rear m+1 point in subinterval is overlapping with the preceding m+1 point in subinterval below；There is 2m+1 control point in preceding SG subinterval, its
Control point coordinates is designated as：Contain N-SG × m control in the SG+1 subinterval
Point, it controls point coordinates to be designated as
Step 2.2, wind-resources local trend is fitted：
If N>3, SG+1 subinterval obtained by step 2.1 is subjected to trend fitting respectively, each subinterval is obtained corresponding
Trend sequence；The control point coordinates of the trend sequence in preceding SG subinterval is designated as：Sg=1,2 ..., SG；The trend sequence in the SG+1 subinterval
Control point coordinates be designated as：For ease of distinguishing,
Title works as N>Trend sequence is obtained for local trend in this step 2.2 when 3；
Step 2.3, the global trend synthesis of wind-resources：
For ease of distinguishing, the trend that this step 2.3 is obtained is called global trend；
If N≤3, trend fitting directly is carried out to wind-resources scene graph { t, w }, obtain its global trend sequence t,glowf,
Subsequently into step 3；
If N>When 3, global trend must be synthesized according to local trend；Specifically building-up process is：Take the 1st subinterval local trend
Preceding m control pointsAs the preceding m control point of the global trend sequence of wind-resources,
It is designated asTake rear N- (SG+1) × m-1 of the SG+1 subinterval local trend
Control pointIt is used as the global trend sequence of wind-resources
(SG+1) × m+2 control point to n-th control point, i.e.,The SG+1 subinterval obtained by step 2.1
There is SG lap, according to being weighted processing to the influence degree of its lap fitted trend between adjacent subarea,
The local trend of the part is weighted as the following formula, the weighting local trend sequence of lap is obtained：
Wherein：The sg lap refers to the rear m+1 control point in the sg subinterval and the preceding m+1 in the sg+1 subinterval
Individual control point is overlapping；For the vertical seat at j-th of control point in the weighting local trend sequence of the sg lap
Mark, j=1,2 ..., m+1, λ1、λ2For weight coefficient, For the office in the sg subinterval
+ m control point ordinates of jth of portion's trend sequence,For j-th of control in the fitting sequence in the sg+1 subinterval
Point ordinate；The control point coordinates of the weighting local trend sequence of the sg lap is designated as：SG can be obtained after processing is weighted to SG lap
Individual weighting local trend sequence, and last control point coordinates of previous weighting local trend sequence and latter weighting office
First control point coordinates of portion's trend sequence is identical；Carry out deduplication operation；The specific of duplicate removal is processed as previous weighting office
Last control point coordinates of portion's trend sequence retains, and removes first control point of its latter weighting local trend sequence
Coordinate；After such duplicate removal, the 1st weighting local trend sequence includes m+1 control point, i.e.,Remaining weighting local trend sequence includes m control point, i.e.,SG weighting local trend sequence after duplicate removal is joined end to end, takes it to be provided as wind
The m+1 control point to (SG+1) × m+1 control point of source overall situation trend sequence, i.e.,
So far, the global trend sequence of the wind-resources comprising N number of control point can be obtained, be also designated as t,glowf}。
4. the wind-resources scene reduction method according to claim 3 based on trend fitting, it is characterised in that step 5 is by such as
Lower step is carried out：
Wind-resources scene set before epicycle is reduced is designated as { TT, WW }bef, the number of wind-resources scene is designated as before epicycle reduction
KK, the sequence number of wind-resources scene is designated as kk, i.e. kk=1, and 2 ..., KK, the probability of correspondence wind-resources scene are designated as pkk；When the 1st time
When performing this step, KK=NN, the probability for being regarded as NN wind-resources scene appearance is equal, is designated as pkk=1/NN；According to { TT, WW
}befIn the corresponding global trend sequence of each wind-resources scene graph calculate general between the global trend sequence of wind-resources two-by-two
Rate weight distance, is designated as PDkk, formula is as follows：
WhereinKk ≠ jj, kk=1,2 ..., KK, jj=1,2 ..., KK；
For k scene of kth, the set containing KK-1 distance can be obtained, { D is designated askk-jj, by { Dkk-jjIn minimum value correspondence
Scene sequence number jj be labeled as JLkk, can be obtained containing KK JL for KK scenekkSet, be designated as { JLkk}；For epicycle
Wind-resources scene set { TT, WW } before reductionbef, the set containing KK probability right distance can be obtained, { PD is designated askk, will
{PDkkIn the corresponding scene sequence number kk of minimum value be labeled as PDmin.
5. the wind-resources scene reduction method according to claim 4 based on trend fitting, it is characterised in that step 6 is by such as
Lower step is carried out：
By serial number PDmin scene from { TT, WW }befIt is middle to delete；In { JLkkIn find out PDmin element JLPDmin, so
Afterwards by the Probability p of serial number PDmin scenePDminAdd to serial number JLPDminScene probabilityI.e.So far, it can obtain the scene set { TT, WW } after epicycle reduction sceneaftAnd its corresponding probability；
The number of times for performing this step is designated as mm, if mm<MM, then go to step 5 and start new round reduction, otherwise end step 5 is to step
Rapid 6 circulation, is met scene set and its probability of scene reduction number requirement.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611248741.7A CN106997407B (en) | 2016-12-29 | 2016-12-29 | Wind resource scene reduction method based on trend fitting |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611248741.7A CN106997407B (en) | 2016-12-29 | 2016-12-29 | Wind resource scene reduction method based on trend fitting |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106997407A true CN106997407A (en) | 2017-08-01 |
| CN106997407B CN106997407B (en) | 2022-08-16 |
Family
ID=59431381
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201611248741.7A Expired - Fee Related CN106997407B (en) | 2016-12-29 | 2016-12-29 | Wind resource scene reduction method based on trend fitting |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106997407B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110321799A (en) * | 2019-06-04 | 2019-10-11 | 武汉大学 | A kind of scene number selection method based on SBR and average between class distance |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20140177915A1 (en) * | 2012-12-26 | 2014-06-26 | Shengyin FAN | Method and apparatus for detecting object |
| CN104156543A (en) * | 2014-08-28 | 2014-11-19 | 国家电网公司 | Major power grid failure successive on-off disturbance setting method for whole-process simulation |
| CN104268436A (en) * | 2014-10-20 | 2015-01-07 | 武汉大学 | Trend fitting based wind power intraday fluctuation continuous period recognition method and system |
| CN105631520A (en) * | 2015-07-16 | 2016-06-01 | 长沙理工大学 | Novel modeling method for extracting random and fuzzy uncertainty characteristics of wind speed |
| CN106230028A (en) * | 2016-09-08 | 2016-12-14 | 安徽电气工程职业技术学院 | A kind of Multipurpose Optimal Method of wind-powered electricity generation water-storage association system |
-
2016
- 2016-12-29 CN CN201611248741.7A patent/CN106997407B/en not_active Expired - Fee Related
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20140177915A1 (en) * | 2012-12-26 | 2014-06-26 | Shengyin FAN | Method and apparatus for detecting object |
| CN104156543A (en) * | 2014-08-28 | 2014-11-19 | 国家电网公司 | Major power grid failure successive on-off disturbance setting method for whole-process simulation |
| CN104268436A (en) * | 2014-10-20 | 2015-01-07 | 武汉大学 | Trend fitting based wind power intraday fluctuation continuous period recognition method and system |
| CN105631520A (en) * | 2015-07-16 | 2016-06-01 | 长沙理工大学 | Novel modeling method for extracting random and fuzzy uncertainty characteristics of wind speed |
| CN106230028A (en) * | 2016-09-08 | 2016-12-14 | 安徽电气工程职业技术学院 | A kind of Multipurpose Optimal Method of wind-powered electricity generation water-storage association system |
Non-Patent Citations (2)
| Title |
|---|
| XIANXUN WANG,ET AL.: "Influence of Sub-Daily Variation on Multi-Fractal Detrended Fluctuation Analysis of Wind Speed Time Series", 《PLOS ONE》 * |
| 王现勋: "基于量维和形维的风水联合出力波动量化评估", 《华中科技大学学报》 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110321799A (en) * | 2019-06-04 | 2019-10-11 | 武汉大学 | A kind of scene number selection method based on SBR and average between class distance |
| CN110321799B (en) * | 2019-06-04 | 2022-11-18 | 武汉大学 | Scene number selection method based on SBR and average inter-class distance |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106997407B (en) | 2022-08-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107591800B (en) | Method for predicting running state of power distribution network with distributed power supply based on scene analysis | |
| CN112257366B (en) | CNF generation method and system for equivalence verification | |
| CN109145003B (en) | Method and device for constructing knowledge graph | |
| CN103761337B (en) | Method and system for processing unstructured data | |
| CN108335010A (en) | A kind of wind power output time series modeling method and system | |
| CN105005029A (en) | Multi-mode radar signal sorting method based on data field hierarchical clustering | |
| CN103268362A (en) | Auxiliary design method of virtual terminals on basis of general template and key character matching | |
| CN101055606A (en) | Multiple clock system integrative circuit plane layout method | |
| CN109976806A (en) | Detection method is cloned based on the matched Java statement block of byte code sequence | |
| CN102521713B (en) | Data processing equipment and data processing method | |
| CN103995873A (en) | Data mining method and data mining system | |
| CN106997407A (en) | Wind-resources scene reduction method based on trend fitting | |
| Jia et al. | KG-FLIP: Knowledge-guided Fashion-domain Language-Image Pre-training for E-commerce | |
| CN106326005A (en) | Automatic parameter tuning method for iterative MapReduce operation | |
| CN102298618A (en) | Method for obtaining matching degree to execute corresponding operations and device and equipment | |
| CN110362620A (en) | A kind of list data structure method based on machine learning | |
| CN107621562B (en) | A kind of single-phase meter shunting accurate judgement method for analyzing stealing | |
| CN106569734B (en) | The restorative procedure and device that memory overflows when data are shuffled | |
| CN104778253B (en) | A kind of method and apparatus that data are provided | |
| CN106874243A (en) | Formula processing method and device based on character string | |
| CN106919762A (en) | A kind of finite element grid array modeling methods | |
| CN110321799A (en) | A kind of scene number selection method based on SBR and average between class distance | |
| CN105654106A (en) | Decision tree generation method and system thereof | |
| CN111310121A (en) | New energy output probability prediction method and system | |
| CN105207233A (en) | Reactive power optimization method based on combination of Metropolis-Hastings and PSO (Particle Swarm Optimization) |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20220816 |