CN110766291A - Horizontal plane daily total radiation data acquisition method based on solar radiation subareas - Google Patents

Horizontal plane daily total radiation data acquisition method based on solar radiation subareas Download PDF

Info

Publication number
CN110766291A
CN110766291A CN201910930612.3A CN201910930612A CN110766291A CN 110766291 A CN110766291 A CN 110766291A CN 201910930612 A CN201910930612 A CN 201910930612A CN 110766291 A CN110766291 A CN 110766291A
Authority
CN
China
Prior art keywords
solar radiation
data
radiation
observation
region
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
Application number
CN201910930612.3A
Other languages
Chinese (zh)
Other versions
CN110766291B (en
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.)
Xian University of Architecture and Technology
Original Assignee
Xian University of Architecture and Technology
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 Xian University of Architecture and Technology filed Critical Xian University of Architecture and Technology
Priority to CN201910930612.3A priority Critical patent/CN110766291B/en
Publication of CN110766291A publication Critical patent/CN110766291A/en
Application granted granted Critical
Publication of CN110766291B publication Critical patent/CN110766291B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063—Operations research, analysis or management
    • G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
    • G06Q10/06315—Needs-based resource requirements planning or analysis
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J1/00—Photometry, e.g. photographic exposure meter
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/08—Construction

Landscapes

  • Engineering & Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Human Resources & Organizations (AREA)
  • Theoretical Computer Science (AREA)
  • Strategic Management (AREA)
  • Economics (AREA)
  • Mathematical Physics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Entrepreneurship & Innovation (AREA)
  • General Business, Economics & Management (AREA)
  • Tourism & Hospitality (AREA)
  • Marketing (AREA)
  • Computational Mathematics (AREA)
  • Operations Research (AREA)
  • Mathematical Optimization (AREA)
  • Data Mining & Analysis (AREA)
  • Mathematical Analysis (AREA)
  • Software Systems (AREA)
  • Development Economics (AREA)
  • Game Theory and Decision Science (AREA)
  • Algebra (AREA)
  • General Engineering & Computer Science (AREA)
  • Quality & Reliability (AREA)
  • Educational Administration (AREA)
  • Databases & Information Systems (AREA)
  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Photometry And Measurement Of Optical Pulse Characteristics (AREA)

Abstract

The invention discloses a method for acquiring horizontal plane daily total radiation data based on solar radiation subareas, which is characterized in that the solar radiation subareas are acquired by using a clustering analysis subarea method based on meteorological data, so that the data quantity participating in the subareas is increased, and the accuracy of the subarea result is improved, thereby improving the accuracy of the solar radiation data acquisition; by adopting a partitioning method combining a clustering algorithm and geographical distribution, the result of partitioning has both objectivity and practicability, so that the accuracy of obtaining solar radiation data is improved; a daily total radiation calculation model for a radiation area is provided, a determination method of a constant coefficient of the area model is provided, a new acquisition method is provided for generating daily total radiation in a non-radiation observation data area in China, and the accuracy of acquiring solar radiation data is improved.

Description

Horizontal plane daily total radiation data acquisition method based on solar radiation subareas
Technical Field
The invention relates to a method for acquiring total solar radiation data, in particular to a method for acquiring total horizontal plane daily radiation data based on solar radiation partitions.
Background
Solar radiation is used as a key meteorological element influencing the indoor thermal environment of a building, the thermal comfort of a human body and the energy consumption of the building, and data of the solar radiation is an important basis for formulating a design strategy of the thermal environment of the building. Along with the continuous deep urbanization process and increasingly highlighted fine and precise construction requirements of China, higher requirements are put forward on radiation data.
Generally, solar radiation data are obtained by direct observation of solar radiation observation stations, but are limited to problems of capital, maintenance and the like, only 98 solar radiation ground observation stations are available in China, the number of the solar radiation ground observation stations is not matched with the development requirement of urbanization in China, the solar radiation data of regions without the solar radiation observation stations cannot be obtained by measurement, and the depth of building energy-saving design is greatly restricted by the shortage of the solar radiation data.
Disclosure of Invention
The invention aims to provide a method for acquiring horizontal plane daily total radiation data based on solar radiation subareas, which is used for solving the problem that solar radiation data cannot be acquired in certain areas without solar radiation observation sites in the prior art.
In order to realize the task, the invention adopts the following technical scheme:
a horizontal plane total solar radiation data acquisition method based on solar radiation subareas is used for acquiring total solar radiation of a target area according to meteorological data acquired by a meteorological observation station, wherein the target area has no radiation observation station, and the method is implemented according to the following steps:
step 1, acquiring meteorological data of a plurality of meteorological observation sites in a geographic range, and acquiring an observation data matrix according to the meteorological data;
clustering the observation data matrix to obtain a plurality of labels;
dividing a geographical range according to the labels to obtain a plurality of solar radiation subareas, wherein each solar radiation subarea comprises at least one region, and one solar radiation subarea comprises a target region;
step 2, acquiring the average clear sky index of the solar radiation subarea of the target area in the year-round;
step 3, selecting a total daily radiant quantity calculation model according to the annual average clear sky index of the solar radiation subarea where the target area is located, and specifically comprising the following steps:
if the average clear sky index of the sun in the successive years is less than or equal to 0.4, the calculation model of the total solar radiation amount isWherein G is total solar radiation in MJ/m2，G0The total daily astronomical radiance is MJ/m2，TmaxThe daily maximum temperature obtained for the meteorological observation site in degrees C.TminThe daily minimum air temperature obtained by the meteorological observation station is shown in the unit of DEG C, S is the sunshine duration obtained by the meteorological observation station and is shown in the unit of h, S0The number of illuminable hours obtained by the meteorological observation station is h; a. b and c are constant coefficients;
if the average clear sky index of the sun in the successive years is more than 0.4 and less than or equal to 0.55, the calculation model of the total solar radiation amount is
Figure BDA0002220185560000022
Wherein E is the average air pressure obtained by the meteorological observation station, the unit is hPa, and d is a constant coefficient;
if the average clear sky index of the sun in the successive years is more than 0.55, the calculation model of the total solar radiation amount is
Figure BDA0002220185560000023
And 4, obtaining the total solar radiation of the target area according to the total daily radiation calculation model.
Further, obtaining a constant coefficient of the daily total radiation calculation model in the step 2 specifically includes:
selecting at least one sample area from all areas included in a solar radiation subarea where a target area is located, wherein the sample area is an area with a meteorological observation station and a radiation observation station;
obtaining meteorological observation data and radiation observation data in each sample area to obtain sample data, wherein the meteorological observation data comprises an illuminable time S0And the sunshine hours S obtained by the meteorological observation station, and the daily maximum air temperature T obtained by the meteorological observation stationmaxDaily minimum air temperature T obtained at the weather observation siteminAnd/or average air pressure E obtained at meteorological observation site, said radiometric observationThe measured data comprises the total solar radiation G of the specimen region;
and performing regression on the selected daily total radiation calculation model by using the sample data to obtain a constant coefficient of the daily total radiation calculation model.
Further, the step 1 specifically includes:
step 1.1, acquiring observation data of N regions with ground meteorological observation sites in a geographic range from t1 to t2, and acquiring cumulative monthly mean value data of sunshine duration and cumulative monthly mean value data of daily average temperature of each region with the meteorological observation sites, wherein N is the number of the ground meteorological observation sites and is a positive integer;
wherein the mean data of the year-round months of the sunshine hours of any region with a meteorological observation site comprises
Figure BDA0002220185560000031
Wherein
Figure BDA0002220185560000032
Represents the annual monthly mean of the number of sunshine hours in the mth month; the data of the annual average temperature and the monthly average value of the daily average temperature in any region with a meteorological observation site comprises
Figure BDA0002220185560000041
Wherein
Figure BDA0002220185560000042
Representing the mean annual average temperature of the mth month;
step 1.2, filling the chronological monthly average value data of the sunshine hours of each region with the meteorological observation site and the chronological monthly average value data of the average air temperature into an Nx 25 matrix to obtain an observation data matrix;
each row in the observation data matrix represents an observation station, the first column represents the number of the observation station, the second column to the thirteenth column represent the annual monthly mean value of sunshine hours from the 1 st month to the 12 th month, and the fourteenth column to the twenty-fifth column represent the annual monthly mean value of the average air temperature from the 1 st month to the 12 th month;
step 1.3, normalizing the numerical values of the second column to the twenty-fifth column in the observation data matrix to obtain a normalized observation data matrix;
step 1.4, carrying out hierarchical clustering on the normalized observation data matrix to obtain the number of label classes;
according to the label number, distributing labels for each region with the weather observation station;
step 1.5, dividing the map corresponding to the geographic range according to the label, and specifically comprising the following steps:
and dividing the region corresponding to the same label and at least comprising one meteorological observation station into a solar radiation subarea to obtain a plurality of solar radiation subareas.
Further, the step 1.4 of performing hierarchical clustering on the normalized observation data matrix to obtain a tag class number specifically includes:
performing hierarchical clustering on the normalized observation data matrix to obtain inflection points of a dispersion square sum curve in the hierarchical clustering;
and taking the value of the deviation square and the quantity of the curve inflection points as a label class number.
Further, at least one region with a meteorological observation site corresponding to the same tag is divided into a solar radiation subarea, a plurality of solar radiation subareas are obtained, and the method specifically comprises the following steps:
1.5.1, dividing regions with the same labels and meteorological observation sites into the same type, and taking the same type of region as a region to obtain a plurality of regions, wherein each region comprises at least one region with meteorological observation sites;
step 1.5.2, mapping the plurality of areas to a map corresponding to a geographical range to obtain a partition map;
step 1.5.3, searching two adjacent areas in the partition map, obtaining a boundary of the two adjacent areas, and specifically performing the following steps:
step 1.5.3.1, searching two adjacent areas in the partition map, namely area A and area B, and respectively obtaining area position point sets of weather observation stations at the edge of area A and area B;
the region position point set of the meteorological observation station at the region edge A comprises I region edge points A, the region position point set of the meteorological observation station at the region edge B comprises J region edge points B, and I and J are positive integers;
step 1.5.3.2, calculating the geodetic distance between the ith A area edge point and the J B area edge points, and taking the midpoint of the geodetic distance between the jth B area edge point which is closest to the geodetic distance and the ith A area edge point as a boundary point, wherein I is less than or equal to I, and J is less than or equal to J;
step 1.5.3.3, repeating step 1.5.3.2 until all boundary points are obtained, connecting all boundary points into a line, and obtaining a boundary line of two adjacent areas;
and 1.5.4, repeating the step 1.5.3 until all boundaries are obtained, and dividing the geographical range according to the boundaries to obtain a plurality of solar radiation partitions.
Compared with the prior art, the invention has the following technical effects:
1. the method for acquiring the total daily radiation data of the horizontal plane based on the solar radiation subareas, provided by the invention, considers that the number of ground meteorological observation stations in China is far more than that of the radiation observation stations, and the method for acquiring the radiation subareas by using the meteorological data increases the data quantity participating in the subareas and improves the accuracy of the subarea results, thereby improving the accuracy of acquiring the solar radiation data;
2. the method for acquiring the horizontal plane daily total radiation data based on the solar radiation subareas adopts a subarea method combining a clustering algorithm and geographical distribution, and the result of division takes both objectivity and practicability into consideration, so that the accuracy of acquiring the solar radiation data is improved;
3. the method for acquiring the horizontal plane total daily radiation data based on the solar radiation subareas, which is provided by the invention, provides a total daily radiation calculation model for a radiation area and a determination method for a constant coefficient of the area model, provides a new acquisition method for generating the total daily radiation in a non-radiation observation data area in China, and improves the accuracy of acquiring the solar radiation data.
Drawings
FIG. 1 is a dispersion sum of squares curve provided in an embodiment of the present invention;
FIG. 2 is a diagram illustrating a hierarchical clustering result provided in an embodiment of the present invention;
FIG. 3 is a schematic view of a demarcation line acquisition provided in an embodiment of the present invention;
fig. 4 is a schematic view of solar radiation zoning provided in an embodiment of the invention.
Detailed Description
The present invention will be described in detail below with reference to the accompanying drawings and examples. So that those skilled in the art can better understand the present invention. It is to be expressly noted that in the following description, a detailed description of known functions and designs will be omitted when it may obscure the subject matter of the present invention.
The following definitions or conceptual connotations relating to the present invention are provided for illustration:
total solar radiation: total solar radiation during a day.
Clear sky index: the ratio of the total solar radiation incident on the horizontal plane to the astronomical radiation.
Sunshine hours: the time of direct sunlight irradiating the ground in one day is obtained by observation of a meteorological observation station.
Number of illuminable hours (astronomical illuminable hours): in the absence of any shading, the time that the sun's center passes from the east to the west horizon, the number of hours it takes for the light to strike the ground, is a fixed number calculated from geographic and astronomical parameters for each area.
Total daily astronomical radiation: the total daily amount of astronomical radiation for each region is a fixed value calculated from geographic and astronomical parameters, without considering atmospheric effects on the earth's surface and determined only by the solar radiation of the solar-terrestrial astronomical relations.
The embodiment discloses a method for acquiring total daily radiation data of a horizontal plane based on solar radiation partitions, which is used for calculating the total daily solar radiation of the horizontal plane of a target area according to meteorological data acquired by a meteorological observation station, wherein the target area is free of the meteorological observation station.
The method comprises the following steps:
step 1, acquiring meteorological data of a plurality of meteorological observation sites in a geographic range, and acquiring an observation data matrix according to the meteorological data;
clustering the observation data matrix to obtain a plurality of labels;
dividing a geographical range according to the labels to obtain a plurality of solar radiation subareas, wherein each solar radiation subarea comprises at least one region, and one solar radiation subarea comprises a target region;
in the invention, the geographical range can be worldwide, Asian or nationwide, and can be determined according to the location of the target area, for example, the target area is a Linzhi area, and the geographical range can be China or Asian; or the target area is the hokkaido, the geographic scope may be all japan or all asia.
In this embodiment, the whole area of China (except Hongkong and Australia) is taken as an example.
In this step, the meteorological data includes sunshine hours, daily average air temperature, average air pressure, average wind speed, average relative humidity, daily deterioration, and the like, but in order to improve the accuracy of the method for finally obtaining the sunshine radiation data, the meteorological data is first screened.
In the embodiment, the offset correlation and the complex correlation coefficient of each meteorological parameter and the total daily radiation are calculated one by using the data of 91 observation stations which have meteorological observation and solar radiation observation in China in 2000-2013 (sunshine hours, daily average temperature, average air pressure, average wind speed, average relative humidity and poor daily) and the total daily radiation observation data, the correlation test result shows that the correlation between the sunshine hours and the daily average temperature and the total daily radiation is higher than that of other meteorological parameters, and the sunshine hours and the daily average temperature are selected as the solar radiation climate partition indexes.
Optionally, the step 1 specifically includes:
step 1.1, acquiring observation data of N regions with ground meteorological observation sites in a geographic range from t1 to t2, and acquiring cumulative monthly mean value data of sunshine duration and cumulative monthly mean value data of daily average temperature of each region with the meteorological observation sites, wherein N is the number of the ground meteorological observation sites and is a positive integer;
wherein the mean data of the year-round months of the sunshine hours of any region with a meteorological observation site comprises
Figure BDA0002220185560000081
Wherein
Figure BDA0002220185560000082
Represents the annual monthly mean of the number of sunshine hours in the mth month; the monthly mean value data of the mean temperature of any region with meteorological observation sites comprises
Figure BDA0002220185560000091
Wherein
Figure BDA0002220185560000092
Representing the mean annual average temperature of the mth month;
in the present embodiment, the cluster data is from the daily data set of the Chinese ground climate data. 641 observation sites with national weather observation data continuously recorded for more than 20 years between 1 month and 2013 and 12 months and 31 days are selected.
And respectively calculating the sunshine hours of 641 stations from 1984 to 2013 and the annual and monthly average value of daily average air temperature.
Figure BDA0002220185560000093
Figure BDA0002220185560000094
The average of the sun hours of the mth month, i.e. m is 1,2 … …, 12. Unit: hours (h).
Sy(m): year y month mSun hours monthly mean value, y 1,2 … …, n, 20<n is less than or equal to 30. Unit: hours (h).
Figure BDA0002220185560000095
Figure BDA0002220185560000096
The average temperature in the mth month and the average temperature in the mth month, m is 1,2 … …, 12. Unit: DEG C.
Ty(m): mean monthly mean temperature of mth month and mth month in year y, y is 1,2 … …, n, 20<n is less than or equal to 30. Unit: DEG C.
Step 1.2, filling the chronological monthly average value data of the sunshine hours of each region with the meteorological observation site and the chronological monthly average value data of the average air temperature into an Nx 25 matrix to obtain an observation data matrix;
each row in the observation data matrix represents an observation station, the first column represents the number of the observation station, the second column to the thirteenth column represent the annual monthly mean value of sunshine hours from the 1 st month to the 12 th month, and the fourteenth column to the twenty-fifth column represent the annual monthly mean value of the average air temperature from the 1 st month to the 12 th month;
in this embodiment, a 641 × 25 observation matrix is established, the number of rows of the matrix is 641, each row represents an observation station, the number of columns of the matrix is 25, the 1 st column is a station number, the 2 nd to 13 th columns are the annual average values of the daily hours of 1 month to 12 months, and the 14 th to 25 th columns are the annual average values of the average air temperatures of 1 month to 12 months.
Step 1.3, normalizing the numerical values of the second column to the twenty-fifth column in the observation data matrix to obtain a normalized observation data matrix;
in the present embodiment, the data in the 2 nd to 25 th columns of the observation data matrix are normalized to eliminate the dimensional difference between the sunshine hours and the daily average air temperature.
Step 1.4, carrying out hierarchical clustering on the normalized observation data matrix to obtain the number of label classes;
according to the label number, distributing labels for each region with the weather observation station;
optionally, the performing hierarchical clustering on the normalized observation data matrix in step 1.4 to obtain the number of the tag classes specifically includes:
performing hierarchical clustering on the normalized observation data matrix to obtain the quantity of dispersion square and curve inflection points in the hierarchical clustering;
and taking the value of the deviation square and the quantity of the curve inflection points as a label class number.
In this embodiment, the solar radiation partitioning is completed by a bottom-up condensation classification method in hierarchical clustering, each sample is a class at the beginning of clustering, the Ward algorithm is selected for class-to-class merging, and the distance between samples is the squared Euclidean (Euclidean) distance. When the classes are sequentially merged, a new dispersion square sum curve is generated, the dispersion square sum curve is shown in figure 1, the best classification number of solar radiation in China is judged to be 8 according to the inflection point of the curve, and the clustering result is shown in figure 2.
Step 1.5, according to the labels, dividing at least one region with a meteorological observation site corresponding to the same label into a solar radiation subarea to obtain a plurality of solar radiation subareas.
In this embodiment, as shown in fig. 2, the site areas included in the same type of tag are mostly concentrated in one geographical area in the spatial distribution, but some site areas are far from the concentrated area of the tag of this type and are located in the concentrated area of other types of tags, which are called outliers, for example, the geographical position of the class 3 outlier is located in the concentrated area of the class 2 tag, and the geographical area where the same type of tag is concentrated is divided into one radiation area without considering the outliers when dividing the solar radiation area.
Optionally, in step 1.5, according to the tag, at least one region with a weather observation station corresponding to the same tag is divided into a solar radiation partition, so as to obtain a plurality of solar radiation partitions, which is specifically executed according to the following steps:
1.5.1, dividing regions with the same labels and meteorological observation sites into the same type, and taking the same type of region as a region to obtain a plurality of regions, wherein each region comprises at least one region with meteorological observation sites;
step 1.5.2, mapping the plurality of areas to a map corresponding to a geographical range to obtain a partition map;
step 1.5.3, searching two adjacent areas in the partition map, obtaining a boundary of the two adjacent areas, and specifically performing the following steps:
step 1.5.3.1, searching two adjacent areas in the partition map, namely area A and area B, and respectively obtaining area position point sets of weather observation stations at the edge of area A and area B;
the region position point set of the meteorological observation station at the region edge A comprises I region edge points A, the region position point set of the meteorological observation station at the region edge B comprises J region edge points B, and I and J are positive integers;
step 1.5.3.2, calculating the geodetic distance between the ith A area edge point and the J B area edge points, and taking the midpoint of the geodetic distance between the jth B area edge point which is closest to the geodetic distance and the ith A area edge point as a boundary point, wherein I is less than or equal to I, and J is less than or equal to J;
step 1.5.3.3, repeating step 1.5.3.2 until all boundary points are obtained, connecting all boundary points into a line, and obtaining a boundary line of two adjacent areas;
and 1.5.4, repeating the step 1.5.3 until all boundaries are obtained, and dividing the geographical range according to the boundaries to obtain a plurality of solar radiation partitions.
In this embodiment, a plurality of points in the 5-class and 6-class edge regions are randomly selected, as shown in fig. 3, i, i +1, i +2, … are edge points of class 5, j, j +1, j +2, j +3, … are edge points of class 6, and the geodetic distance between two points is calculated by using the inverse gaussian mean norm formula according to the longitude and latitude of the edge points.
Taking point i as an example, the distance from point j is Si,jIt is shown that similar distances from j +1, j +2, j +3, … are denoted S, respectivelyi,j+1，Si,j+2，Si,j+3… are provided. Taking Min (S)i,j,Si,j+1,Si,j+2,Si,j+3…), assuming the result is Si,jThen, the midpoint of the distance between the I point and the j point is calculated and used as a boundary point, the calculation of the (I + 1) th, I + 2) th and (…) th points is sequentially completed according to the above steps to obtain a plurality of boundary points, and the points are connected by a smooth curve to be used as the boundary lines of the 5 th class and the 6 th class.
In this example, 8 solar radiation sections are finally obtained, as shown in fig. 4.
Step 2, acquiring the annual average clear sky index of the solar radiation subarea of the target area;
in the present embodiment, it is preferred that,
Figure BDA0002220185560000131
Kt: sun clear sky index (dimensionless), G is total daily radiation (MJ/m)2)，G0Is the total daily astronomical emission (MJ/m)2)。
In the formula ISCFor solar constant, 4.921MJ/m is taken2；E0Is an eccentricity correction factor for the earth orbit.
Figure BDA0002220185560000134
Figure BDA0002220185560000135
The average clear sky index of the day of the year,
Figure BDA0002220185560000136
and (4) the average clear sky index of the day of the y year, and n is the recording age limit of the radiation observation data.
Taking the target region, namely the tender river, as an example, the target region belongs to a solar radiation 4 region, and the target region comprises 7 sample regions with solar radiation observation data, namely a desert river, Aihui, Hiragel, Fuyu, Jia Mus, Harbin and Yangji. The average clear sky index of 7 sample areas in the year-round days is calculated by using the formula and is shown in table 1, and the average value is taken as the average clear sky index of the area in the year-round.
TABLE 1 average clear sky index of sample area over the years
Figure BDA0002220185560000137
Step 3, selecting a daily solar total radiation quantity calculation model according to the annual average clear sky index of the solar radiation subarea where the target area is located, and specifically comprising the following steps:
if the average clear sky index of the sun in the successive years is less than or equal to 0.4, the calculation model of the total solar radiation amount is
Figure BDA0002220185560000141
Wherein G is total solar radiation in MJ/m2，G0The total daily astronomical radiance is MJ/m2，TmaxThe daily maximum temperature obtained for the meteorological observation site in degrees C.TminThe daily minimum air temperature obtained by the meteorological observation station is shown in the unit of DEG C, S is the sunshine duration obtained by the meteorological observation station and is shown in the unit of h, S0The number of illuminable hours obtained by the meteorological observation station is h; a. b and c are constant coefficients;
if the average clear sky index of the sun in the successive years is more than 0.4 and less than or equal to 0.55, the calculation model of the total solar radiation amount is
Figure BDA0002220185560000142
Wherein E is the average air pressure obtained by the meteorological observation station, the unit is hPa, and d is a constant coefficient;
if the average clear sky index of the sun in the successive years is more than 0.55, the calculation model of the total solar radiation amount is
Figure BDA0002220185560000143
Optionally, obtaining a constant coefficient of the daily total radiation calculation model in step 2 specifically includes:
selecting at least one sample area from all areas included in a solar radiation subarea where a target area is located, wherein the sample area is an area with a meteorological observation station and a radiation observation station;
obtaining meteorological observation data and radiation observation data in each sample area to obtain sample data, wherein the meteorological observation data comprises an illuminable time S0And the sunshine hours S obtained by the meteorological observation station, and the daily maximum air temperature T obtained by the meteorological observation stationmaxDaily minimum air temperature T obtained at the weather observation siteminAnd/or average atmospheric pressure E obtained by a meteorological observation station, wherein the radiation observation data comprise the total solar radiation G of the sample area;
and performing regression on the selected daily total radiation calculation model by using the sample data to obtain a constant coefficient of the daily total radiation calculation model.
In the embodiment, taking each solar radiation subarea as an example, the constant coefficients of the daily total radiation calculation model are obtained by regression of data of 2000-2013 of a station having meteorological and solar radiation observation data in an area, and are shown in table 2.
TABLE 2 solar radiation area model
Figure BDA0002220185560000151
Total daily radiance G (MJ/m) in the table2) (ii) a Total daily astronomical radiance G0(MJ/m2) (ii) a Sunshine duration S (h); illuminable time S0(h) (ii) a Average air pressure E (0.1 hPa); daily minimum and maximum temperature Tmin，Tmax(0.1℃)；δ: the declination angle (deg) of the weft yarn,
Figure BDA0002220185560000162
latitude (deg).
And 4, acquiring the daily total solar radiation of the target area according to the daily total radiation calculation model.
The total solar radiation quantity of a target area day by day is obtained by using the calculation model provided by the invention, and the sunshine hours S (h) of a meteorological station in the area are used; average air pressure E (0.1 hPa); and observing data at the lowest daily temperature, the highest daily temperature Tmin and Tmax (0.1 ℃) to generate the total daily radiant quantity of the non-radiative observation site region in 2000-2013.
In the present embodiment, the lijiang area is taken as an example, and the sunshine hours S of the area in 5, 18 and 2017 are 4.5h and the sunshine hours S are013.4h, total daily astronomical radiance G0＝40.04MJ/m2. The region belongs to the radiation VII region, so that a is 0.208 and b is 0.517, then
Figure BDA0002220185560000163
Namely the total daily radiant emittance of 15.57MJ/m in the Lijiang area in 5-18 months in 20172。
Example two
In order to prove that the method provided by the patent can be used for calculating the total daily radiation of the areas without the radiation observation stations, the estimation error and the national total daily radiation distribution condition are verified respectively.
And selecting a station for verification in each radiation area, repeating the selecting and verifying stations in each area for three times, wherein the selected verifying stations are distributed in different directions of the edge of the area each time, the verifying stations have radiation observation data and do not participate in regression of model coefficients, 8 radiation areas are formed, and the number of the verifying stations is 24 in total. And calculating the average absolute error percentage MAE% and the root mean square error percentage RMSE% of the area model and the verification station self-built model to the daily total radiation of the verification station respectively. The mean value of the area model MAE% was 11.9%, the RMSE% was 15.9%, and the corresponding results for the station model were 11.4% and 15.3%, respectively. The estimation error of the area model is slightly higher than the station model but differs by less than 1%.
The method provided by the patent is used for calculating the daily total radiant emittance of 819 stations for radiationless observation data in China from 2000 to 2013, and the result shows thatThe highest solar radiation center appears in the river of lion-quan and the river of Yalu-Tibetan-Bujiang in the southwest Tibet, and the average annual total radiation amount is 7500MJ/m2The above; the second highest region appears in the Qinghai Chauda basin and extends to the northeast, reaches the junction areas of three provinces including northwest part of Gansu, west part of inner Mongolia and east part of Xinjiang, and the average annual total radiation amount can reach 6500MJ/m2The above; then the area of the Tarim basin and the Turpan basin, the annual total radiation quantity is 6000MJ/m on average2Left and right; the Tianshan mountain in Xinjiang is a low-value region of solar radiation in the west, and the average annual total radiation amount is 5000MJ/m2～5500MJ/m2. The total solar radiation is highest in the east region, the solar radiation is relatively low in the southeast and northeast regions, the Sichuan basin is a low-value solar radiation region in China, and the average total annual radiation is 4000MJ/m2The following. Comparing the above results with the distribution diagram of the solar total irradiance in the middle year of the building climate zone standard (GB50178-93), it is found that the distribution rule and the variation trend of the total irradiance are completely consistent with the standard, and the numerical value is slightly different from the standard because the used data years are different, to sum up, the method for acquiring the horizontal plane daily total irradiance data based on the solar irradiance zone provided by the present patent is completely applicable in China.

Claims (5)

1. A horizontal plane total solar radiation data acquisition method based on solar radiation subareas is used for acquiring total solar radiation of a target area according to meteorological data acquired by a meteorological observation station, wherein the target area is not provided with a radiation observation station, and the method is characterized by comprising the following steps of:
step 1, acquiring meteorological data of a plurality of meteorological observation sites in a geographic range, and acquiring an observation data matrix according to the meteorological data;
clustering the observation data matrix to obtain a plurality of labels;
dividing a geographical range according to the labels to obtain a plurality of solar radiation subareas, wherein each solar radiation subarea comprises at least one region, and one solar radiation subarea comprises a target region;
step 2, acquiring the average clear sky index of the solar radiation subarea of the target area in the year-round;
step 3, selecting a total daily radiant quantity calculation model according to the annual average clear sky index of the solar radiation subarea where the target area is located, and specifically comprising the following steps:
if the average clear sky index of the sun in the successive years is less than or equal to 0.4, the calculation model of the total solar radiation amount is
Figure FDA0002220185550000011
Wherein G is total solar radiation in MJ/m2，G0The total daily astronomical radiance is MJ/m2，TmaxThe daily maximum temperature obtained for the meteorological observation site in degrees C.TminThe daily minimum air temperature obtained by the meteorological observation station is shown in the unit of DEG C, S is the sunshine duration obtained by the meteorological observation station and is shown in the unit of h, S0Is an illuminable number and has the unit of h; a. b and c are constant coefficients;
if the average clear sky index of the sun in the successive years is more than 0.4 and less than or equal to 0.55, the calculation model of the total solar radiation amount is
Figure FDA0002220185550000012
Wherein E is the average air pressure obtained by the meteorological observation station, the unit is hPa, and d is a constant coefficient;
if the average clear sky index of the sun in the successive years is more than 0.55, the calculation model of the total solar radiation amount is
And 4, obtaining the total solar radiation of the target area according to the total daily radiation calculation model.
2. The method for acquiring horizontal plane daily total radiation data based on solar radiation subareas according to claim 1, wherein obtaining constant coefficients of the daily total radiation calculation model in step 2 specifically comprises:
selecting at least one sample area from all areas included in a solar radiation subarea where a target area is located, wherein the sample area is an area with a meteorological observation station and a radiation observation station;
obtaining meteorological observation data and radiation observation data in each sample area to obtain sample data, wherein the meteorological observation data comprises an illuminable time S0And the sunshine hours S obtained by the meteorological observation station, and the daily maximum air temperature T obtained by the meteorological observation stationmaxDaily minimum air temperature T obtained at the weather observation siteminAnd/or average atmospheric pressure E obtained by a meteorological observation station, wherein the radiation observation data comprise the total solar radiation G of the sample area;
and performing regression on the selected daily total radiation calculation model by using the sample data to obtain a constant coefficient of the daily total radiation calculation model.
3. The method for acquiring horizontal plane daily total radiation data based on solar radiation subareas according to claim 1, wherein the step 1 specifically comprises:
step 1.1, acquiring observation data of N regions with ground meteorological observation sites in a geographic range from t1 to t2, and acquiring cumulative monthly mean value data of sunshine duration and cumulative monthly mean value data of daily average temperature of each region with the meteorological observation sites, wherein N is the number of the ground meteorological observation sites and is a positive integer;
wherein the mean data of the year-round months of the sunshine hours of any region with a meteorological observation site comprises
Figure FDA0002220185550000031
Wherein
Figure FDA0002220185550000032
Represents the annual monthly mean of the number of sunshine hours in the mth month; the data of the annual average temperature and the monthly average value of the daily average temperature in any region with a meteorological observation site comprises
Figure FDA0002220185550000033
Wherein
Figure FDA0002220185550000034
Representing the mean annual average temperature of the mth month;
step 1.2, filling the chronological monthly average value data of the sunshine hours of each region with the meteorological observation site and the chronological monthly average value data of the average air temperature into an Nx 25 matrix to obtain an observation data matrix;
each row in the observation data matrix represents an observation station, the first column represents the number of the observation station, the second column to the thirteenth column represent the annual monthly mean value of sunshine hours from the 1 st month to the 12 th month, and the fourteenth column to the twenty-fifth column represent the annual monthly mean value of the average air temperature from the 1 st month to the 12 th month;
step 1.3, normalizing the numerical values of the second column to the twenty-fifth column in the observation data matrix to obtain a normalized observation data matrix;
step 1.4, carrying out hierarchical clustering on the normalized observation data matrix to obtain the number of label classes;
according to the label number, distributing labels for each region with the weather observation station;
step 1.5, dividing the map corresponding to the geographic range according to the label, and specifically comprising the following steps:
and dividing the region corresponding to the same label and at least comprising one meteorological observation station into a solar radiation subarea to obtain a plurality of solar radiation subareas.
4. The method for acquiring horizontal plane daily total radiation data based on solar radiation subareas as claimed in claim 3, wherein the step 1.4 is to perform hierarchical clustering on the normalized observation data matrix to obtain the number of label classes, and specifically comprises:
performing hierarchical clustering on the normalized observation data matrix to obtain a dispersion square and a curve inflection point in the hierarchical clustering;
and taking the value of the deviation square and the quantity of the curve inflection points as a label class number.
5. The method for acquiring total daily radiation data in a horizontal plane based on solar radiation partitions as claimed in claim 3, wherein at least one region with a meteorological observation station corresponding to the same tag is divided into a solar radiation partition, and a plurality of solar radiation partitions are acquired, and the method is specifically executed according to the following steps:
1.5.1, dividing regions with the same labels and meteorological observation sites into the same type, and taking the same type of region as a region to obtain a plurality of regions, wherein each region comprises at least one region with meteorological observation sites;
step 1.5.2, mapping the plurality of areas to a map corresponding to a geographical range to obtain a partition map;
step 1.5.3, searching two adjacent areas in the partition map, obtaining a boundary of the two adjacent areas, and specifically performing the following steps:
step 1.5.3.1, searching two adjacent areas in the partition map, namely area A and area B, and respectively obtaining area position point sets of weather observation stations at the edge of area A and area B;
the region position point set of the meteorological observation station at the region edge A comprises I region edge points A, the region position point set of the meteorological observation station at the region edge B comprises J region edge points B, and I and J are positive integers;
step 1.5.3.2, calculating the geodetic distance between the ith A area edge point and the J B area edge points, and taking the midpoint of the geodetic distance between the jth B area edge point which is closest to the geodetic distance and the ith A area edge point as a boundary point, wherein I is less than or equal to I, and J is less than or equal to J;
step 1.5.3.3, repeating step 1.5.3.2 until all boundary points are obtained, connecting all boundary points into a line, and obtaining a boundary line of two adjacent areas;
and 1.5.4, repeating the step 1.5.3 until all boundaries are obtained, and dividing the geographical range according to the boundaries to obtain a plurality of solar radiation partitions.
CN201910930612.3A 2019-09-29 2019-09-29 Horizontal plane daily total radiation data acquisition method based on solar radiation subareas Expired - Fee Related CN110766291B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910930612.3A CN110766291B (en) 2019-09-29 2019-09-29 Horizontal plane daily total radiation data acquisition method based on solar radiation subareas

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910930612.3A CN110766291B (en) 2019-09-29 2019-09-29 Horizontal plane daily total radiation data acquisition method based on solar radiation subareas

Publications (2)

Publication Number Publication Date
CN110766291A true CN110766291A (en) 2020-02-07
CN110766291B CN110766291B (en) 2022-07-26

Family

ID=69330756

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910930612.3A Expired - Fee Related CN110766291B (en) 2019-09-29 2019-09-29 Horizontal plane daily total radiation data acquisition method based on solar radiation subareas

Country Status (1)

Country Link
CN (1) CN110766291B (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111967655A (en) * 2020-07-28 2020-11-20 中国南方电网有限责任公司 Short-term load prediction method and system
CN115097492A (en) * 2022-08-10 2022-09-23 湖南北云科技有限公司 Ionospheric error elimination method and related equipment

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109034626A (en) * 2018-07-30 2018-12-18 西安建筑科技大学 A kind of evaluation method that west area heating in solar energy building utilizes
JP6552077B1 (en) * 2019-01-18 2019-07-31 株式会社ヒデ・ハウジング Solar radiation normalization statistical analysis system, solar radiation normalization statistical analysis method and solar radiation normalization statistical analysis program

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109034626A (en) * 2018-07-30 2018-12-18 西安建筑科技大学 A kind of evaluation method that west area heating in solar energy building utilizes
JP6552077B1 (en) * 2019-01-18 2019-07-31 株式会社ヒデ・ハウジング Solar radiation normalization statistical analysis system, solar radiation normalization statistical analysis method and solar radiation normalization statistical analysis program

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
于瑛等: "太阳辐射两级区化方法及其应用", 《土木建筑与环境工程》 *
刘大龙等: "以晴空指数为主要依据的太阳辐射分区", 《建筑科学》 *
艾士博等: "不同气候区建筑获得太阳辐射的动态分布特征", 《建筑节能》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111967655A (en) * 2020-07-28 2020-11-20 中国南方电网有限责任公司 Short-term load prediction method and system
CN115097492A (en) * 2022-08-10 2022-09-23 湖南北云科技有限公司 Ionospheric error elimination method and related equipment

Also Published As

Publication number Publication date
CN110766291B (en) 2022-07-26

Similar Documents

Publication Publication Date Title
CN106909722B (en) A kind of accurate inversion method of large area of temperature near the ground
CN106908415B (en) A kind of big region crops time of infertility Soil Moisture Monitoring method based on amendment NDVI time series
Zhang et al. Spatial and temporal variability in the net primary production of alpine grassland on the Tibetan Plateau since 1982
Zhang et al. Land cover classification of the North China Plain using MODIS_EVI time series
CN105184423B (en) A kind of wind power plant cluster wind speed forecasting method
CN111598045A (en) Remote sensing farmland change detection method based on object spectrum and mixed spectrum
CN110766291B (en) Horizontal plane daily total radiation data acquisition method based on solar radiation subareas
CN104657935A (en) Near-surface air temperature inversion method
CN115204691B (en) Urban artificial heat emission estimation method based on machine learning and remote sensing technology
CN112857583A (en) Vegetation temperature and bare ground temperature estimation method based on satellite remote sensing data
CN116127871A (en) Dynamics method for reconstructing tropical ocean near-surface flow field by utilizing multi-source satellite data
Tian et al. Spatiotemporal differentiation and attribution of land surface temperature in China in 2001–2020
Zhang et al. Variation in vegetation greenness in spring across eastern China during 1982–2006
Beyazit et al. Modeling of the hourly horizontal solar diffuse radiation in Sanliurfa, Turkey
He et al. Retrieval of rugged mountainous areas Land Surface Temperature from high-spatial-resolution thermal infrared remote sensing data
Liao et al. Fast and accurate estimation of solar irradiation on building rooftops in Hong Kong: A machine learning-based parameterization approach
Wang et al. Influence of complex topography on global solar radiation in the Yangtze River Basin
CN110705010B (en) Method for simulating ground surface heat island at night on the basis of remote sensing
CN109885808B (en) Near-surface meteorological element calculation method
Sun et al. Construction of a spatially gridded heat flux map based on airborne flux Measurements using remote sensing and machine learning methods
Malekpour et al. Modeling of relationship between land use/cover and land surface temperature using ASTER datasets
Yang et al. Terrain revised model for air temperature in mountainous area based on DEMs: A case study in Yaoxian county
Zhou et al. Land surface albedo estimation with Chinese GF-1 WFV data in Northwest China
CN106777668A (en) The method that a kind of raising region humidity simulates effect
CN113190633A (en) Dual-distance-based climate space partitioning method

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

Granted publication date: 20220726

CF01 Termination of patent right due to non-payment of annual fee