CN104268429A - Satellite-borne SAR (Synthetic Aperture Radar) based offshore wind energy resource remote sensing method and system - Google Patents

Satellite-borne SAR (Synthetic Aperture Radar) based offshore wind energy resource remote sensing method and system Download PDF

Info

Publication number
CN104268429A
CN104268429A CN201410542837.9A CN201410542837A CN104268429A CN 104268429 A CN104268429 A CN 104268429A CN 201410542837 A CN201410542837 A CN 201410542837A CN 104268429 A CN104268429 A CN 104268429A
Authority
CN
China
Prior art keywords
wind
wind energy
wind speed
sar
remote sensing
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.)
Pending
Application number
CN201410542837.9A
Other languages
Chinese (zh)
Inventor
范开国
于兴修
王瑶
傅斌
常俊芳
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hubei University
Original Assignee
Hubei University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Hubei University filed Critical Hubei University
Priority to CN201410542837.9A priority Critical patent/CN104268429A/en
Publication of CN104268429A publication Critical patent/CN104268429A/en
Pending legal-status Critical Current

Links

Landscapes

  • Wind Motors (AREA)

Abstract

The invention provides a satellite-borne SAR (Synthetic Aperture Radar) based offshore wind energy resource remote sensing method. The satellite-borne SAR based offshore wind energy resource remote sensing method comprises step 1, obtaining a long time sequence SAR sea surface wind speed value in every grid and a sea surface wind direction value in every grid after grid processing and generating a wind rose diagram; step 2, obtaining different altitudes of sea surface wind speed values, estimating a scale factor A and a shape factor k in every grid and calculating an average wind speed value of different altitudes of representative wind energy characteristics; step 3, obtaining a wind energy density value; step 4, generating a remote sensing result distribution diagram of the different altitudes of SAR offshore wind energy resources, obtaining the macroscopic and microscopic time and space characteristics of the offshore area wind energy resources according to the remote sensing result distribution diagram and giving out the microscopic time and space characteristics of a potential wind power station and potential wind power station wind energy resources; step 5, giving out fan space layout parameters of the potential wind power station and combining with a microcosmic research result of the win energy resources of the potential wind power station to give out the fan layout of the potential wind power station.

Description

Based on inshore offshore wind energy resource remote sensing technique and the system of satellite-borne SAR
Technical field
The present invention relates to wind energy resources remote sensing technology field, particularly a kind of inshore offshore wind energy resource remote sensing technique based on satellite-borne SAR and system.
Background technology
Inshore offshore wind energy resource because of its not land occupation resource, selection space, position large, be conducive to selecting place, few by environmental constraints, wind speed is higher, more horn of plenty, transport and lifting are with favourable conditions, Wind turbines single-machine capacity is larger for wind energy resources, the more high advantage of annual utilization hours, countries in the world are developing the offshore wind farm industry of this country one after another.
The exploitation of inshore offshore wind energy resource first to understand can utilize wind energy resources position, reserves and area coverage etc., the early stage of wind energy turbine set builds and then needs to know the spatial-temporal distribution characteristic of wind energy resources within the scope of construction and the blower fan spatial arrangement technical parameter of potential wind energy turbine set, and the laying for wind turbine provides Macrocosm and microcosm foundation.Therefore wind energy turbine set erection offshore wind energy resource research work in earlier stage directly has influence on the economy of project construction.
For the research of inshore offshore wind energy resource, survey according to land meteorological station the Wind Data that classic method that wind statistics result carries out extrapolating cannot obtain coverage density between inshore offshore high-altitude, and its precision of numerical simulation technology and resolution all lower.Along with the development of satellite remote sensing technology, scatterometer etc. for offshore wind energy resource research provide new technological means, but its resolution is too low, and because the impact being subject to land return cannot obtain the effective wind field information of inshore, this loses effective value and meaning to the research of inshore high resolving power offshore wind energy resource.
Synthetic-aperture radar (Synthetic Aperture Radar, SAR) operation wavelength is centimetre length, can not by the impact of the factor such as cloud layer, weather, there is round-the-clock, round-the-clock, high resolving power, high spatial coverage density advantage, and along with more and more day by day maturation with SAR high resolving power ocean surface wind retrieving technology that succeeds in sending up being loaded with SAR microwave remote sensor satellite, the high resolving power Wind Data of application SAR inverting carries out wind energy resources Remote Sensing Study, has become the new technical means of inshore offshore wind energy resource research.
When the existing wind energy resources remote sensing technology based on SAR image is only carried out macroscopical to inshore offshore wind energy resource, empty characteristic research, the Characteristics of Wind Field of wind turbine wake zone also only gives some qualitative descriptions, wind energy resources cannot be obtained and enrich the wind energy resources fine structure change in region and the quantitative variation characteristic data in space of wind turbine wake zone small scale wind field, for potential wind energy turbine set provides relevant blower fan spatial arrangement technical parameter, therefore the micro-analysis out of true of wind energy resources, is also unfavorable for the blower fan layout of wind energy turbine set.
Summary of the invention
The invention provides a kind of can the inshore offshore wind energy resource remote sensing technique based on satellite-borne SAR of the quantitative variation characteristic in space of quantitative test wind turbine wake zone small scale wind field and the change of wind energy resources fine structure and system.
Based on an inshore offshore wind energy resource remote sensing technique for satellite-borne SAR, it comprises the following steps:
S1, the process of long-term sequence SAR sea surface wind speed retrieval is carried out to the SAR remote sensing images obtained, to obtain after gridding process in each grid wind direction of ocean surface value in long-term sequence SAR ocean surface wind speed value and each grid; Wind rose map is generated according to wind direction of ocean surface value in each grid;
S2, obtain the ocean surface wind speed value of different altitude height by long-term sequence SAR ocean surface wind speed value in each grid, and estimate scale factor A in each grid in Weibull distributed model and form factor k by Maximum Likelihood Estimation Method; Pass through obtained scale factor A and form factor k and calculate the mean wind speed value that different altitude height represents wind power features;
The air density values of S3, calculating Different Altitude; Acquisition wind energy concentration value is calculated by the air density values of mean wind speed value and Different Altitude;
S4, generate the remote sensing distribution of results figure (mean wind speed figure, wind energy concentration figure, form factor k figure, scale factor A figure and wind rose map) of the SAR offshore offshore wind energy resource of different altitude height according to mean wind speed value, wind energy concentration value, form factor k, scale factor A and wind rose map, and according to remote sensing distribution of results figure obtain inshore Sea area wind energy resources Macrocosm and microcosm time, empty characteristic, when providing the microcosmic of potential wind energy turbine set and potential wind energy turbine set wind energy resources, empty characteristic;
S5, when analyzing the built single blower fan of wind park many, to the fine structure change of High Resolution SAR wind field Spatial Variation and potential wind energy turbine set wind energy resources, provide the blower fan spatial arrangement technical parameter of potential wind energy turbine set, in conjunction with the microexamination result of the wind energy resources of potential wind energy turbine set, provide the blower fan layout of potential wind energy turbine set.
Based on an inshore offshore wind energy resource remote sensing system for satellite-borne SAR, it comprises with lower module:
Inverting module, for carrying out the process of long-term sequence SAR sea surface wind speed retrieval to the SAR remote sensing images obtained, to obtain after gridding process in each grid wind direction of ocean surface value in long-term sequence SAR ocean surface wind speed value and each grid; And for generating wind rose map according to wind direction of ocean surface value in each grid;
Parameter determination module, for obtaining the ocean surface wind speed value of different altitude height by long-term sequence SAR ocean surface wind speed value in each grid, and for estimated different height above sea level in each grid by Maximum Likelihood Estimation Method Weibull distributed model in scale factor A and form factor k; Also calculate for the scale factor A by obtaining and form factor k the mean wind speed value that different altitude height represents wind power features;
Wind energy concentration determination module, for calculating the air density values of Different Altitude; And obtain wind energy concentration value for being calculated by the air density values of mean wind speed value and Different Altitude;
Analysis module, for generating the remote sensing distribution of results figure of the SAR inshore offshore wind energy resource of different height above sea level according to mean wind speed value, wind energy concentration value, form factor k, scale factor A and wind rose map, and for obtain according to remote sensing distribution of results figure coastal wind energy resources regional macro and microcosmic time, empty characteristic, provide potential wind energy turbine set; For according to the built single blower fan of wind energy turbine set many time, to the fine structure change of High Resolution SAR wind field Spatial Variation and potential wind energy turbine set wind energy resources, provide the blower fan spatial arrangement technical parameter of potential wind energy turbine set, in conjunction with the wind energy resources microexamination result of potential wind energy turbine set, provide the blower fan layout of potential wind energy turbine set.
Inshore offshore wind energy resource remote sensing technique based on satellite-borne SAR provided by the invention and system, by the remote sensing distribution of results figure by generating the SAR inshore offshore wind energy resource of different altitude height according to mean wind speed value, wind energy concentration value, form factor k, scale factor A and wind rose map, and can according to remote sensing distribution of results figure obtain coastal wind energy resources regional macro and microcosmic time, empty characteristic.Can the fine structure change of the Spatial Variation of quantitative test wind turbine wake zone small scale wind field and potential wind energy turbine set wind energy resources, provide the blower fan spatial arrangement technical parameter of potential wind energy turbine set, be conducive to the rational deployment in blower fan space.
Accompanying drawing explanation
Fig. 1 is the inshore offshore wind energy resource remote sensing technique process flow diagram based on satellite-borne SAR that embodiment of the present invention provides;
Fig. 2 is the sub-process figure of step S1 in Fig. 1;
Fig. 3 is the inshore offshore wind energy resource remote sensing system structural drawing based on satellite-borne SAR that embodiment of the present invention provides;
Fig. 4 is the minor structure block diagram of inverting module in Fig. 3.
Embodiment
As shown in Figure 1, embodiments provide a kind of inshore offshore wind energy resource remote sensing technique based on satellite-borne SAR, it comprises the following steps:
S1, the process of long-term sequence SAR sea surface wind speed retrieval is carried out to the SAR remote sensing images obtained, to obtain after gridding process in each grid wind direction of ocean surface value in long-term sequence SAR ocean surface wind speed value and each grid.Wind rose map is generated according to wind direction of ocean surface value in each grid.
Wind rose map (abbreviation wind rose) is also wind direction frequency rose diagram, and it is according to a certain area each wind direction of average statistics and the percent value of wind speed for many years, and draws by a certain percentage, and generally multiplex 8 or 16 compass azimuths represent.
Alternatively, as shown in Figure 2, described step S1 comprises following sub-step:
The SAR remote sensing images of S11, inverting long-term sequence are to obtain Wind Speed Inversion data.
SAR remote sensing images can be bought by the satellite remote sensing ground receiving station of Chinese Academy of Sciences's remote sensing and digital earth research institute and obtain.Utilize C-band SAR sea surface wind speed retrieval pattern conventional in the world and polarizability pattern inverting SAR ocean surface wind speed.
C-band SAR sea surface wind speed retrieval model conventional in the world mainly comprises CMOD4, CMOD5 and CMOD-IFR2, and combines for the different polarizability models of HH polarization, cross polarization conversion.
S12, the mode data of long-term sequence SAR ocean surface wind speed value and long-term sequence is carried out correlation analysis, remove typhoon data, retain relevant long-term sequence SAR Wind Speed Inversion data.
Because utilize the wind direction information of SAR image extracting directly will cause comparatively big error, (spatial resolution is 1 ° × 1 ° to the long-term sequence wind direction data (NOGAPS) that kriging analysis method can be adopted to provide USN's global atmosphere operational forecast system, temporal resolution 6 hours) or Environmental forecasting centre (NCEP) and American National Center for Atmospheric Research (NCAR) long-term sequence of combining release analyzes wind direction data (NCEP-NCAR) again, and (spatial resolution is 2.5 ° × 2.5 °, temporal resolution 6 hours) carry out space interpolation, obtain the wind direction of ocean surface information with SAR Pixel domain position phase registration.Mode data can pass through MM5 pattern acquiring, MM5 (Mesoscale Model5) has the ability of multinest ability, non-static(al) dynamic mode and four-dimensional assimilate, and can run on a computer platform, simulate or forecast the general circulation of mesoscale and regional scale.
Wherein, the long-term sequence SAR remote sensing images spatial coverage after reservation needs to reach more than 165 scapes or 165 scapes; The NOGAPS wind direction data of long-term sequence and NCEP-NCAR are analyzed wind direction data again and can be downloaded by Environmental forecasting centre.
S13, gridding process is carried out to the long-term sequence Wind Speed Inversion data retained, to obtain after gridding process in each grid wind direction of ocean surface value in long-term sequence SAR ocean surface wind speed value and each grid.And generate wind rose map according to wind direction of ocean surface value in each grid.
Sizing grid is set to the different resolution grade being not more than 1000 meters × 1000 meters, grid can be divided into 1000,900,800,700,600,500,400,300,200 and 100 meters etc., obtains wind direction of ocean surface value in long-term sequence SAR ocean surface wind speed value in each net point and each grid.The data fit WGS84 standard of gridding, shows with the projection pattern of longitude and latitude grid.
S2, obtain the ocean surface wind speed value of different altitude height by long-term sequence SAR ocean surface wind speed value in each grid, and estimate scale factor A in each grid in different height Weibull distributed model above sea level and form factor k by Maximum Likelihood Estimation Method.Pass through obtained scale factor A and form factor k and calculate the mean wind speed value representing different altitude height wind power features in each grid.
When the probability of occurrence of various grade wind speed differs widely, the wind energy calculated just has very big-difference, therefore needs to utilize wind velocity distributing paremeter model to calculate wind energy concentration.Research shows, the two parameter Weibull wind speed profile model containing scale factor A (m/s) and form parameter k (dimensionless) is a kind of more satisfactory pattern for wind energy calculates.Formula below gives the Weibull wind speed profile function of 0 ~ V integration, as long as and know parameter A and k, the distribution characteristics of wind speed just can be determined.
F ( V ) = P ( v ≤ V ) = 1 - exp [ - ( V A ) k ]
Alternatively, as follows by the formula of the ocean surface wind speed value of long-term sequence SAR ocean surface wind speed value acquisition different altitude height in each grid in described step S2:
wherein, V nand V 1be respectively Z nand Z 1the mean wind speed of At The Height, α is the shear index of mean wind speed along with height change, and the value of α is chosen as 0.09.
In Wind Power Generation, axial fan hub height generally build more than off sea 70m height or higher height in.And SAR inverting obtains is 10m At The Height Wind Data, obtain the wind energy resources distribution situation of marine differing heights more comprehensively, need the Wind Data of sea 10m At The Height SAR inverting obtained to be converted to different altitude height to calculate the ocean surface wind speed of differing heights, can above formula be passed through describe the exponential law Wind outline of VERTICAL SHEAR OF WIND feature, be similar to and realize wind speed vertical trimming change feature.
The computing formula of mean wind speed value is as follows:
wherein Γ is gamma function, and k is form factor (dimensionless), A is scale factor (m/s).
The air density values of S3, calculating Different Altitude.And calculate acquisition wind energy concentration value by the air density values of mean wind speed value and Different Altitude.
Alternatively, in described step S3, the computing formula of the air density values of Different Altitude is as follows:
ρ=(P 0/ (Rt)) e (-gz/ (Rt)), wherein P 0for standard flat atmospheric pressure (101.325kPa), g is acceleration of gravity (9.8m/s2), z is on-the-spot height (m), R is gas law constant J/ (KgK)), t is temperature value, alternatively, the relational expression dt/dz=-1/100 DEG C/m that temperature value t can be changed with the change of height by atmospheric temperature calculates.
On sea level general in international standard, air themperature and pressure are respectively 288.15k and 101.325kPa, so the atmospheric density of standard sea level is 1.225kg/m3, in each grid, the atmospheric density at different altitude height place also can obtain in approximate treatment.
The computing formula of wind energy concentration value is as follows:
W ‾ = 1 2 ρ ∫ 0 ∞ V 3 f ( V ) dV = 0.5 ρ A 3 Γ ( 3 k + 1 ) , Wherein ρ is air density values.
S4, mean wind speed value according in each grid, wind energy concentration value, form factor k, scale factor A and wind rose map generate the remote sensing result space distribution plan (mean wind speed figure, wind energy concentration figure, form factor k figure and scale factor A figure) of the SAR inshore offshore wind energy resource of whole inshore marine site wind rose map and different altitude height, and according to remote sensing distribution of results figure obtain coastal wind energy resources region time, empty characteristic.
Alternatively, in described step S4, the spatial resolution of the remote sensing result space distribution plan of inshore offshore wind energy resource is respectively 1000,900,800,700,600,500,400,300,200 and 100 meters; Coastal wind energy resources region time, empty characteristic comprise the horizontal-spatial distribution of offshore sea waters wind energy resources, the vertical space distribution of different above sea level height, season distribution, month distribution character.
(sizing grid is set to be not less than 500 meters × 500 meters the remote sensing distribution of results figure of SAR inshore offshore wind energy resource, comprise 1000,900,800,700,600 and 500 meters) result, can from macroscopically carrying out macrovisual study to inshore offshore wind energy resource distribution situation, mainly comprise on the horizontal-spatial distribution to inshore offshore wind energy resource, different elevation vertical space distribution highly, seasonal variations and Monthly changes and give research and analysis, find out inshore offshore wind energy resource and enrich region.
(sizing grid is set to be not more than 500 meters × 500 meters inshore offshore wind energy resource to be enriched to the remote sensing distribution of results figure of the SAR offshore wind energy resource in region, comprise 400,300,200 and 100 meters) result, from microcosmic, micro-analysis research is carried out to wind energy resources, mainly comprise to enrich on the horizontal space change in region, the vertical space change of different above sea level height, seasonal variations and Monthly changes wind energy resources and give research and analysis, provide the small scale fine structure variation characteristic that wind energy resources enriches region.
Based on satellite-borne SAR inshore offshore wind energy resource remote sensing technology and research flow process, the wind energy remote sensing distribution of results figure in seashore research marine site, coastal waters is obtained by the mean wind speed value of the long-term sequence in each net point, wind energy concentration value, form factor k and scale factor A, and eliminate grid burrs on edges by the medium filtering of 3 × 3 or 5 × 5, obtain the inshore SAR offshore wind energy resource remote sensing distribution results after smoothing processing.And then it is remote-sensing distributed based on high resolution SAR inshore offshore wind energy resource, the horizontal space of offshore sea waters wind energy resources is changed, the vertical space change of different above sea level height, seasonal variations and Monthly changes carry out macrovisual study, find out inshore offshore wind energy resource and enrich region, and can be the planning of possible candidate's wind energy turbine set reference frame is provided; Wind energy resources is enriched on the horizontal space change in region, different elevation vertical space change highly, seasonal variations and Monthly changes simultaneously and carry out microexamination, for the blower fan space rational deployment of potential wind energy turbine set lays the foundation.
S5, SAR Ocean Wind-field by long-term sequence, analyze the built single blower fan of wind energy turbine set many time, to the Small-scale Space variation characteristic of wind field, change with the fine structure of potential wind energy turbine set wind energy resources (mean wind speed, wind energy concentration, form factor k, scale factor A and wind direction), provide the blower fan spatial arrangement technical parameter of potential wind energy turbine set, in conjunction with the microexamination result of the wind energy resources of potential wind energy turbine set, provide the blower fan layout of potential wind energy turbine set.
Analyze the small scale wind field Spatial Variation of single blower fan (fan blade height is clear and definite) at different wind direction, different wind speed, alternatively, spacing laid by blower fan along wind direction direction is more than 75% distance returning to original wind speed or original wind speed, lay spacing along the blower fan perpendicular to wind direction direction and be similarly more than 75% distance returning to original wind speed or original wind speed, like this in conjunction with the wind energy resources microexamination result in potential wind energy turbine set marine site, rational blower fan spatial arrangement technical parameter just can be provided.
For the blower fan spatial arrangement technical parameter of the potential wind energy turbine set of China, can first for Donghai Bridge in Shanghai wind energy turbine set, carry out the space small scale variation characteristic of wind field in the different wind speed in wind electric field blower tail district, different wind direction situation, provide the blower fan spatial arrangement technical parameter of potential wind energy turbine set.Wherein, blower fan tail district wind speed mainly comprises wind speed along wind direction direction with perpendicular to the variation range in wind direction direction and change intensity with the wind speed variation characteristic of distance, provides the wind field Small-scale Space variation characteristic around wind electric field blower.The change distance of wind speed is generally to return to more than 75% of original wind speed or original wind speed for standard, the change intensity of wind speed is defined as ((originally wind speed-wind speed now)/(originally wind speed)), if wind speed is unchanged, be 0%, like this by the wind field space characteristics analysis of change intensity, give the small scale wind field space characteristics in single blower fan tail district.
For the microexamination of the wind energy resources of potential wind energy turbine set, by enriching the research and analysis of the horizontal space change in region, the vertical space change of different above sea level height, seasonal variations and Monthly changes to wind energy resources, provide the small scale fine structure variation characteristic that wind energy resources enriches region, mainly comprise the most centrostigma of potential wind energy turbine set wind energy resources, the fine structure variation characteristic of the wind direction in potential wind energy turbine set marine site, mean wind speed, wind energy concentration, form factor k and scale factor A.
Alternatively, in potential wind energy turbine set region, be basic point by reference to the longitude and latitude of the most centrostigma of wind energy resources and elevation information, by the fine structure variation characteristic of the wind direction in potential wind energy turbine set marine site, mean wind speed, wind energy concentration, form factor k and scale factor A, in conjunction with the small scale wind field space characteristics analysis to built wind field single blower fan tail district, provide potential wind field blower fan spatial arrangement technical parameter, be conducive to the rational deployment in blower fan space.
Inshore offshore wind energy resource remote sensing technique based on satellite-borne SAR provided by the invention, the remote sensing distribution of results figure of different SAR inshore offshore wind energy resource highly above sea level is generated by the mean wind speed value according to different spatial resolutions (spatial resolution comprise 1000,900,800,700,600,500,400,300,200 with 100 meters), wind energy concentration value, form factor k, scale factor A and wind direction, and the macroscopic property in coastal wind energy resources region can be obtained according to remote sensing distribution of results figure, the planning for wind energy turbine set provides possible microcosmic candidate site.Micro-analysis can be carried out according to remote sensing distribution of results figure to potential wind energy turbine set simultaneously, and quantitative test can build wind turbine wake zone small scale wind field Spatial Variation, provide the reference technique parameter of potential wind electric field blower spatial arrangement, be conducive to the rational deployment in wind electric field blower space.
As shown in Figure 3, the embodiment of the present invention also provides a kind of inshore offshore wind energy resource remote sensing system based on satellite-borne SAR, and it comprises with lower module:
Inverting module 10, for carrying out the process of SAR sea surface wind speed retrieval to the long-term sequence SAR remote sensing images obtained, to obtain wind direction of ocean surface value in the SAR ocean surface wind speed value of long-term sequence in each grid after gridding process and each grid.And for generating wind rose map according to wind direction of ocean surface value in each grid.
Preferably, as shown in Figure 4, described inverting module 10 comprises as lower unit:
Wind Speed Inversion data capture unit 11, for the SAR remote sensing images of inverting long-term sequence to obtain Wind Speed Inversion data.
Verification unit 12, for the mode data of long-term sequence SAR ocean surface wind speed value and long-term sequence is carried out correlation analysis, remove typhoon data, retain remaining long-term sequence SAR Wind Speed Inversion data, the long-term sequence SAR remote sensing images spatial coverage after reservation needs to reach more than 165 scapes or 165 scapes.
Gridding processing unit 13, for carrying out gridding process to the long-term sequence SAR Wind Speed Inversion data retained, to obtain after gridding process in each grid wind direction of ocean surface value in long-term sequence SAR ocean surface wind speed value and each grid.Gridding processing unit 13 is also for generating wind rose map according to wind direction of ocean surface value in each grid.
Parameter determination module 20, for obtaining the ocean surface wind speed value of different altitude height by long-term sequence SAR ocean surface wind speed value in each grid, and for by the scale factor A in Maximum Likelihood Estimation Method estimation Weibull distributed model and form factor k.Parameter determination module 20 also for by obtain the scale factor A of different altitude height and the mean wind speed value of form factor k computational representation wind power features value.
Preferably, as follows by the formula of the ocean surface wind speed value of long-term sequence SAR ocean surface wind speed value acquisition different altitude height in each grid in described parameter determination module 20:
wherein, V nand V 1be respectively Z nand Z 1the mean wind speed of At The Height, α is the shear index of mean wind speed along with height change, and the value of α is chosen as 0.09;
The computing formula of mean wind speed value is as follows:
wherein Γ is gamma function, and k is form factor, A is scale factor.
Wind energy concentration determination module 30, for calculating the air density values of Different Altitude.And obtain wind energy concentration value for being calculated by the air density values of mean wind speed value and Different Altitude.
Alternatively, in described wind energy concentration determination module 30, the computing formula of the air density values of Different Altitude is as follows:
ρ=(P 0/ (Rt)) e (-gz/ (Rt)), wherein P 0for standard flat atmospheric pressure, g is acceleration of gravity, and z is on-the-spot height, and R is gas law constant, and t is temperature value, and alternatively, the relational expression dt/dz=-1/100 DEG C/m that temperature value t can be changed with the change of height by atmospheric temperature calculates.
The computing formula of wind energy concentration value is as follows:
W ‾ = 1 2 ρ ∫ 0 ∞ V 3 f ( V ) dV = 0.5 ρ A 3 Γ ( 3 k + 1 ) , Wherein ρ is air density values.
Analysis module 40, the remote sensing distribution of results figure of SAR inshore offshore wind energy resource of marine site different altitude height, different spatial resolutions, different time resolution is covered for generating SAR image according to mean wind speed value, wind energy concentration value, form factor k, scale factor A and wind rose map, analysis module 40 also for obtain according to remote sensing distribution of results figure coastal wind energy resources regional macro, microcosmic time, empty characteristic, find out inshore offshore wind energy resource and enrich region, and can be the planning of possible candidate's wind energy turbine set reference frame is provided;
For by change the horizontal space of potential wind energy turbine set marine site offshore wind energy resource, the vertical space of different height above sea level changes, the research and analysis of seasonal variations and Monthly changes, provides the small scale fine structure variation characteristic that wind energy resources enriches region;
For building the small scale wind field Spatial Variation of the single blower fan of wind energy turbine set (fan blade height is clear and definite) at different wind direction, different wind speed by the analysis of long-term sequence SAR ocean surface wind speed, alternatively, spacing laid by blower fan along wind direction direction is more than 75% distance returning to original wind speed or original wind speed, lay spacing along the blower fan perpendicular to wind direction direction and be similarly more than 75% distance returning to original wind speed or original wind speed, to blowing machine spatial arrangement technical parameter;
For the microexamination result of the wind energy resources of the blower fan spatial arrangement technical parameter by obtaining according to built wind energy turbine set single blower fan analysis and potential wind energy turbine set, providing the blower fan spatial arrangement technical parameter of potential wind energy turbine set, carrying out the blower fan layout of wind energy turbine set.
Alternatively, in described analysis module 40, the characteristic in inshore offshore wind energy resource region comprises the macroscopic property of the horizontal space change of offshore sea waters wind energy resources, the vertical space change of different above sea level height, seasonal variations and Monthly changes, also comprises that horizontal space change, the Bu Tong vertical space of height above sea level that inshore offshore wind energy resource enriches the wind energy resources (wind energy resources most centrostigma, wind direction, mean wind speed, wind energy concentration, form factor k and scale factor A) in region changes, the microscopic characteristics of seasonal variations and Monthly changes small scale.
Alternatively, in described analysis module 40, blower fan spatial arrangement technical parameter comprises built single blower fan tail district wind speed and mainly comprises wind speed along wind direction direction and the scope changed with wind speed perpendicular to the distance in wind direction direction and wind speed change intensity with the wind speed variation characteristic of distance.
In this instructions, each embodiment adopts the mode of going forward one by one to describe, and what each embodiment stressed is the difference with other embodiments, between each embodiment identical similar portion mutually see.
Professional can also recognize further, in conjunction with unit and the algorithm steps of each example of embodiment disclosed herein description, can realize with electronic hardware, computer software or the combination of the two, in order to the interchangeability of hardware and software is clearly described, in the above description according to the functional composition and the step that generally describe each example.These functions perform with hardware or software mode actually, depend on application-specific and the design constraint of technical scheme.Professional and technical personnel can use distinct methods to realize described function to each specifically should being used for, but this realization should not exceed scope of the present invention.
The software module that the method described in conjunction with embodiment disclosed herein or the step of algorithm can directly use hardware, processor to perform, or the combination of the two is implemented.Software module can be placed in any other forms of storage medium known in random access memory, internal memory, ROM (read-only memory), electrically programmable ROM, electrically erasable ROM, register, hard disk, moveable magnetic disc, CD-ROM or technical field.
Be understandable that, for the person of ordinary skill of the art, other various corresponding change and distortion can be made by technical conceive according to the present invention, and all these change the protection domain that all should belong to the claims in the present invention with distortion.

Claims (11)

1., based on an inshore offshore wind energy resource remote sensing technique for satellite-borne SAR, it is characterized in that, it comprises the following steps:
S1, the process of long-term sequence SAR sea surface wind speed retrieval is carried out to the SAR remote sensing images obtained, to obtain after gridding process in each grid wind direction of ocean surface value in long-term sequence SAR ocean surface wind speed value and each grid; Wind rose map is generated according to wind direction of ocean surface value in each grid;
S2, obtain the ocean surface wind speed value of different altitude height by long-term sequence SAR ocean surface wind speed value in each grid, and estimate scale factor A in each grid in Weibull distributed model and form factor k by Maximum Likelihood Estimation Method; Pass through obtained scale factor A and form factor k and calculate the mean wind speed value that different altitude height represents wind power features;
The air density values of S3, calculating Different Altitude; Acquisition wind energy concentration value is calculated by the air density values of mean wind speed value and Different Altitude;
S4, according to mean wind speed value, wind energy concentration value, form factor k, scale factor A and wind rose map generate the remote sensing distribution of results figure of the SAR offshore offshore wind energy resource of different altitude height, remote sensing distribution of results figure comprises mean wind speed figure, wind energy concentration figure, form factor k schemes, scale factor A figure and wind rose map, and the horizontal-spatial distribution of inshore Sea area wind energy resources is obtained according to remote sensing distribution of results figure, the vertical space distribution of different height above sea level, season distribution, month the Macrocosm and microcosm such as distribution character time, empty characteristic, when providing the microcosmic of potential wind energy turbine set and potential wind energy turbine set wind energy resources, empty characteristic,
S5, when analyzing the built single blower fan of wind park many, to the fine structure change of High Resolution SAR wind field Spatial Variation and potential wind energy turbine set wind energy resources, provide the blower fan spatial arrangement technical parameter of potential wind energy turbine set, in conjunction with the microexamination result of the wind energy resources of potential wind energy turbine set, provide the blower fan layout of potential wind energy turbine set.
2., as claimed in claim 1 based on the inshore offshore wind energy resource remote sensing technique of satellite-borne SAR, it is characterized in that, described step S1 comprises following sub-step:
The SAR remote sensing images of S11, inverting long-term sequence are to obtain Wind Speed Inversion data;
S12, the mode data of long-term sequence SAR ocean surface wind speed value and long-term sequence is carried out correlation analysis, remove typhoon data, retain remaining Wind Speed Inversion data;
S13, gridding process is carried out to the Wind Speed Inversion data retained, to obtain after gridding process in each grid wind direction of ocean surface value in long-term sequence SAR ocean surface wind speed value and each grid; And generate wind rose map according to wind direction of ocean surface value in each grid.
3. as claimed in claim 2 based on the inshore offshore wind energy resource remote sensing technique of satellite-borne SAR, it is characterized in that, the formula being obtained the ocean surface wind speed value of different altitude height in described step S2 by long-term sequence SAR ocean surface wind speed value in each grid is as follows:
wherein, V nand V 1be respectively Z nand Z 1the mean wind speed of At The Height, α is the shear index of mean wind speed along with height change, and the value of α is chosen as 0.09;
The computing formula of mean wind speed value is as follows:
wherein Γ is gamma function, and k is form factor, A is scale factor.
4., as claimed in claim 3 based on the offshore offshore wind energy resource remote sensing technique of satellite-borne SAR, it is characterized in that,
In described step S3, the computing formula of the air density values of Different Altitude is as follows:
ρ=(P 0/ (Rt)) e (-gz/ (Rt))wherein P 0for standard flat atmospheric pressure, g is acceleration of gravity, and z is on-the-spot height, and R is gas law constant, and t is temperature value;
The computing formula of wind energy concentration value is as follows:
W ‾ = 1 2 ρ ∫ 0 ∞ V 3 f ( V ) dV = 0.5 ρ A 3 Γ ( 3 k + 1 ) , Wherein ρ is air density values.
5. as claimed in claim 4 based on the inshore offshore wind energy resource remote sensing technique of satellite-borne SAR, it is characterized in that, in described step S4, the characteristic in inshore offshore wind energy resource region comprises the macroscopic property of the horizontal space change of offshore sea waters wind energy resources, the vertical space change of different above sea level height, seasonal variations and Monthly changes, comprises the microscopic characteristics that wind energy resources enriches the horizontal space change of the wind energy resources in region, the vertical space change of different height above sea level, seasonal variations and Monthly changes simultaneously.
6. as claimed in claim 5 based on the inshore offshore wind energy resource remote sensing technique of satellite-borne SAR, it is characterized in that, in described step S5, the blower fan spatial arrangement technical parameter of potential wind energy turbine set comprises the wind speed variation characteristic of built wind energy turbine set single blower fan tail district wind speed with distance, mainly comprise wind speed along wind direction direction and perpendicular to the distance in wind direction direction with wind speed variation range and wind speed change intensity, wherein change intensity is chosen as below 25% or 25%, change distance then selects the blower fan along wind direction direction to lay apart from for wind speed being more than 75% distance returning to original wind speed or original wind speed, blower fan perpendicular to wind direction direction is laid and is similarly more than 75% space length returning to the original wind speed of original wind speed, calculate the different wind speed provided based on long-term sequence SAR ocean surface wind speed, blower fan spatial arrangement technical parameter in wind direction situation, in conjunction with the microexamination result of the wind energy resources of potential wind energy turbine set, microexamination result comprises the most centrostigma of wind energy resources, wind rose map and mean wind speed, wind energy concentration, the horizontal space change of form factor k and scale factor A, the vertical space change of different height above sea level, the microscopic characteristics of seasonal variations and Monthly changes small scale, carry out the blower fan layout of potential wind energy turbine set.
7., based on an inshore offshore wind energy resource remote sensing system for satellite-borne SAR, it is characterized in that, it comprises with lower module;
Inverting module, for carrying out the process of long-term sequence SAR sea surface wind speed retrieval to the SAR remote sensing images obtained, to obtain after gridding process in each grid wind direction of ocean surface value in long-term sequence SAR ocean surface wind speed value and each grid; And for generating wind rose map according to wind direction of ocean surface value in each grid;
Parameter determination module, the ocean surface wind speed value of different altitude height is obtained for the long-term sequence SAR ocean surface wind speed value by retaining in each grid, and for being estimated scale factor A in each grid in Weibull distributed model and form factor k by Maximum Likelihood Estimation Method; Also calculate for the scale factor A by obtaining and form factor k the mean wind speed value that different altitude height represents wind power features;
Wind energy concentration determination module, for calculating the air density values of Different Altitude; And obtain wind energy concentration value for being calculated by the air density values of mean wind speed value and Different Altitude;
Analysis module, for generating the remote sensing distribution of results figure of the SAR inshore offshore wind energy resource of different waters height according to the mean wind speed value in each grid, wind energy concentration value, form factor k, scale factor A and wind rose map, and for obtain according to remote sensing distribution of results figure coastal wind energy resources regional macro time, empty characteristic, provide the potential wind energy turbine set that wind energy resources is abundant; Enrich the high resolving power wind energy resources remote sensing distribution of results figure in region for generating wind energy resources, and for the microcosmic according to high-definition remote sensing distribution of results map analysis potential wind field region time, empty fine structure variation characteristic; For according to the built single blower fan of wind field many time, to High Resolution SAR wind field Spatial Variation, provide rational blower fan spatial arrangement technical parameter; For the blower fan spatial arrangement technical parameter provided according to the built single blower fan analysis of wind energy turbine set, in conjunction with the fine structure mutation analysis result of potential wind energy turbine set wind energy resources, provide the blower fan spatial arrangement technical parameter of potential wind energy turbine set, for the blower fan layout of carrying out potential wind energy turbine set provides spatial arrangement foundation.
8., as claimed in claim 7 based on the inshore offshore wind energy resource remote sensing system of satellite-borne SAR, it is characterized in that, described inverting module comprises as lower unit:
Wind Speed Inversion data capture unit, for many scapes SAR remote sensing images of inverting long-term sequence to obtain Wind Speed Inversion data;
Verification unit, for the mode data of long-term sequence SAR ocean surface wind speed value and long-term sequence is carried out correlation analysis, remove typhoon data, retain remaining SAR Wind Speed Inversion data, the long-term sequence SAR remote sensing images spatial coverage after reservation needs to reach more than 165 scapes or 165 scapes;
Gridding processing unit, for carrying out gridding process to the SAR Wind Speed Inversion data retained, to obtain after gridding process in each grid wind direction of ocean surface value in long-term sequence SAR ocean surface wind speed value and each grid; And for generating wind rose map according to wind direction of ocean surface value in each grid.
9. as claimed in claim 8 based on the inshore offshore wind energy resource remote sensing system of satellite-borne SAR, it is characterized in that, the formula being obtained the ocean surface wind speed value of different altitude height in described parameter determination module by long-term sequence SAR ocean surface wind speed value in each grid is as follows:
wherein, V nand V 1be respectively Z nand Z 1the mean wind speed of At The Height, α is the trimming index of mean wind speed along with height change, and the value of α is chosen as 0.09;
The computing formula of mean wind speed value is as follows:
wherein Γ is gamma function, and k is form factor, A is scale factor.
10., as claimed in claim 9 based on the inshore offshore wind energy resource remote sensing system of satellite-borne SAR, it is characterized in that,
In described wind energy concentration determination module, the computing formula of the air density values of Different Altitude is as follows:
ρ=(P 0/ (Rt)) e (-gz/ (Rt)), wherein P 0for standard flat atmospheric pressure, g is acceleration of gravity, and z is on-the-spot height, and R is gas law constant, and t is temperature value;
The computing formula of wind energy concentration value is as follows:
W ‾ = 1 2 ρ ∫ 0 ∞ V 3 f ( V ) dV = 0.5 ρ A 3 Γ ( 3 k + 1 ) , Wherein ρ is air density values.
11. as claimed in claim 10 based on the offshore offshore wind energy resource remote sensing system of satellite-borne SAR, it is characterized in that, in described analysis module the characteristic in coastal wind energy resources region comprise inshore offshore wind energy resource horizontal space change, the vertical space change of different above sea level height, seasonal variations and Monthly changes macroscopic property; Comprise the microscopic characteristics that wind energy resources enriches the horizontal space change of the high resolving power wind energy resources in region, the vertical space change of different above sea level height, seasonal variations and Monthly changes simultaneously; When comprising the single blower fan of built wind field many, to High Resolution SAR wind field Spatial Variation, provide rational blower fan spatial arrangement technical parameter simultaneously; Comprise the blower fan spatial arrangement technical parameter provided according to the built single blower fan analysis of wind energy turbine set simultaneously, in conjunction with the fine structure mutation analysis result of potential wind energy turbine set wind energy resources, provide the blower fan spatial arrangement technical parameter of potential wind energy turbine set, for the blower fan layout of carrying out potential wind energy turbine set provides spatial arrangement foundation.
CN201410542837.9A 2014-10-15 2014-10-15 Satellite-borne SAR (Synthetic Aperture Radar) based offshore wind energy resource remote sensing method and system Pending CN104268429A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410542837.9A CN104268429A (en) 2014-10-15 2014-10-15 Satellite-borne SAR (Synthetic Aperture Radar) based offshore wind energy resource remote sensing method and system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410542837.9A CN104268429A (en) 2014-10-15 2014-10-15 Satellite-borne SAR (Synthetic Aperture Radar) based offshore wind energy resource remote sensing method and system

Publications (1)

Publication Number Publication Date
CN104268429A true CN104268429A (en) 2015-01-07

Family

ID=52159950

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410542837.9A Pending CN104268429A (en) 2014-10-15 2014-10-15 Satellite-borne SAR (Synthetic Aperture Radar) based offshore wind energy resource remote sensing method and system

Country Status (1)

Country Link
CN (1) CN104268429A (en)

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105469326A (en) * 2015-12-24 2016-04-06 东北电力大学 Novel fan arrangement method for wind power plant
CN107767030A (en) * 2017-09-25 2018-03-06 浙江大学 A kind of offshore wind energy resource appraisal procedure based on multi-source remote sensing satellite wind speed correction
CN110322038A (en) * 2018-03-29 2019-10-11 北京金风科创风电设备有限公司 Method and equipment for automatically arranging fans based on mesoscale data
CN111868533A (en) * 2018-03-20 2020-10-30 三菱电机株式会社 Wind flow sensing system and method for determining a velocity field of a wind flow
CN112115406A (en) * 2020-09-28 2020-12-22 自然资源部第二海洋研究所 Ocean internal mesoscale vortex inversion method and system based on remote sensing sea surface data
CN112612916A (en) * 2020-12-29 2021-04-06 深圳航天宏图信息技术有限公司 Method and device for generating inspection error spatial distribution map of ocean satellite data
CN112946643A (en) * 2021-01-28 2021-06-11 华东师范大学 Offshore wind power extraction method and system based on time sequence radar remote sensing
CN115693666A (en) * 2022-12-30 2023-02-03 中国华能集团清洁能源技术研究院有限公司 Offshore wind farm generated energy determination method and system based on satellite inversion
CN116953703A (en) * 2023-07-21 2023-10-27 国家卫星海洋应用中心 Offshore wind energy assessment method, device and equipment
CN117111068A (en) * 2023-10-19 2023-11-24 南京信大卫星应用研究院有限公司 Sea surface wind field monitoring system based on satellite scatterometer data

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20110169686A1 (en) * 2008-05-29 2011-07-14 Gian Luca Eusebi Borzelli Target detection in a sar-imaged sea area
CN104036108A (en) * 2014-03-06 2014-09-10 浙江海洋学院 Joint inversion method of sea surface wind field and sea waves

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20110169686A1 (en) * 2008-05-29 2011-07-14 Gian Luca Eusebi Borzelli Target detection in a sar-imaged sea area
CN104036108A (en) * 2014-03-06 2014-09-10 浙江海洋学院 Joint inversion method of sea surface wind field and sea waves

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
常俊芳: "近岸海洋风能的星载SAR遥感技术及应用研究", 《中国博士学位论文全文数据库(基础科学辑)》 *

Cited By (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105469326A (en) * 2015-12-24 2016-04-06 东北电力大学 Novel fan arrangement method for wind power plant
CN107767030A (en) * 2017-09-25 2018-03-06 浙江大学 A kind of offshore wind energy resource appraisal procedure based on multi-source remote sensing satellite wind speed correction
CN107767030B (en) * 2017-09-25 2021-10-01 浙江大学 Offshore wind energy resource assessment method based on multi-source remote sensing satellite wind speed correction
CN111868533A (en) * 2018-03-20 2020-10-30 三菱电机株式会社 Wind flow sensing system and method for determining a velocity field of a wind flow
CN110322038A (en) * 2018-03-29 2019-10-11 北京金风科创风电设备有限公司 Method and equipment for automatically arranging fans based on mesoscale data
CN110322038B (en) * 2018-03-29 2022-07-15 北京金风科创风电设备有限公司 Method and equipment for automatically arranging fans based on mesoscale data
CN112115406B (en) * 2020-09-28 2024-01-12 自然资源部第二海洋研究所 Ocean internal mesoscale vortex inversion method and system based on remote sensing sea surface data
CN112115406A (en) * 2020-09-28 2020-12-22 自然资源部第二海洋研究所 Ocean internal mesoscale vortex inversion method and system based on remote sensing sea surface data
CN112612916A (en) * 2020-12-29 2021-04-06 深圳航天宏图信息技术有限公司 Method and device for generating inspection error spatial distribution map of ocean satellite data
CN112612916B (en) * 2020-12-29 2024-02-06 深圳航天宏图信息技术有限公司 Method and device for generating inspection error space distribution diagram of marine satellite data
CN112946643B (en) * 2021-01-28 2022-07-12 华东师范大学 Offshore wind power extraction method and system based on time sequence radar remote sensing
CN112946643A (en) * 2021-01-28 2021-06-11 华东师范大学 Offshore wind power extraction method and system based on time sequence radar remote sensing
CN115693666A (en) * 2022-12-30 2023-02-03 中国华能集团清洁能源技术研究院有限公司 Offshore wind farm generated energy determination method and system based on satellite inversion
CN116953703A (en) * 2023-07-21 2023-10-27 国家卫星海洋应用中心 Offshore wind energy assessment method, device and equipment
CN117111068A (en) * 2023-10-19 2023-11-24 南京信大卫星应用研究院有限公司 Sea surface wind field monitoring system based on satellite scatterometer data
CN117111068B (en) * 2023-10-19 2024-03-22 南京信大卫星应用研究院有限公司 Sea surface wind field monitoring system based on satellite scatterometer data

Similar Documents

Publication Publication Date Title
CN104268429A (en) Satellite-borne SAR (Synthetic Aperture Radar) based offshore wind energy resource remote sensing method and system
Bolch et al. Status and change of the cryosphere in the extended Hindu Kush Himalaya region
Mughal et al. Wind modelling, validation and sensitivity study using Weather Research and Forecasting model in complex terrain
Al-Yahyai et al. Review of the use of numerical weather prediction (NWP) models for wind energy assessment
Fierro et al. Application of a lightning data assimilation technique in the WRF-ARW model at cloud-resolving scales for the tornado outbreak of 24 May 2011
Koračin et al. Transition of stratus into fog along the California coast: Observations and modeling
Bright et al. The sensitivity of the numerical simulation of the southwest monsoon boundary layer to the choice of PBL turbulence parameterization in MM5
Gadad et al. Offshore wind power resource assessment using Oceansat-2 scatterometer data at a regional scale
Badger et al. Wind class sampling of satellite SAR imagery for offshore wind resource mapping
Garreaud et al. Orographic precipitation in coastal southern Chile: Mean distribution, temporal variability, and linear contribution
Yang et al. Evaluating the impact of the COSMIC RO bending angle data on predicting the heavy precipitation episode on 16 June 2008 during SoWMEX-IOP8
Favier et al. Atmospheric drying as the main driver of dramatic glacier wastage in the southern Indian Ocean
Helmis et al. A comparative study and evaluation of mixing-height estimation based on sodar-RASS, ceilometer data and numerical model simulations
Lagouvardos et al. The meteorological model BOLAM at the National Observatory of Athens: assessment of two-year operational use
Pérez-Landa et al. Mesoscale circulations over complex terrain in the Valencia coastal region, Spain–Part 2: Modeling CO 2 transport using idealized surface fluxes
Hegdahl et al. An event-based approach to explore selected present and future atmospheric river–induced floods in Western Norway
Bastin et al. Diurnal cycle of water vapor as documented by a dense GPS network in a coastal area during ESCOMPTE IOP2
Gascón et al. Analysis of a localized flash-flood event over the central Mediterranean
Yan et al. The benefit of GPS zenith delay assimilation to high‐resolution quantitative precipitation forecasts: A case‐study from COPS IOP 9
Teixeira et al. Surface to boundary layer coupling in the urban area of Lisbon comparing different urban canopy models in WRF
Squitieri et al. On the forecast sensitivity of MCS cold pools and related features to horizontal grid spacing in convection-allowing WRF simulations
Vel’tishchev et al. Short-range forecast of heavy precipitation and strong wind using the convection-allowing WRF models
Yushkov et al. Modeling an urban heat island during extreme frost in Moscow in January 2017
Chen et al. Comparative study on WRF model simulations from the viewpoint of optimum ship routing
Láska et al. High-resolution numerical simulation of summer wind field comparing WRF boundary-layer parametrizations over complex Arctic topography: case study from central Spitsbergen

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication

Application publication date: 20150107

RJ01 Rejection of invention patent application after publication