Monthly meteorological prediction data hourly downscaling method based on distribution adjustment time mapping

Abstract

The invention belongs to the field of building performance simulation, and provides a monthly meteorological prediction data hourly downscaling method based on distribution adjustment time mapping. Comprising the steps of representative city determination, typical meteorological year data collection, global climate mode and monthly average resolution data set determination, main variable data set acquisition, candidate statistical distribution determination, maximum likelihood estimation, optimal distribution screening, empirical cumulative distribution function construction, distribution correction, quantile calculation, quantile mapping, distribution adjustment and model prediction. According to the method, hourly downscaling of monthly prediction data of different weather variables in the atmospheric circulation model is realized through a specific distribution adjustment and quantile mapping method, so that the physical characteristics of the data and the statistical characteristics of historical data are kept, and the long-term trend and extreme weather conditions of GCM prediction can be reflected; and reliable future meteorological data is provided for building performance simulation.

Description

Technical Field

[0001] The invention relates to the field of building performance simulation, and in particular to an hourly downscaling method for monthly meteorological forecast data based on distribution-adjusted time mapping. Background Art

[0002] Buildings account for approximately 28% of global greenhouse gas emissions and 30% of global energy consumption during their operational phase. As climate conditions change, buildings are likely to experience increased cooling demand, altered heating demand, and increased stress on building materials and systems, necessitating adaptive design and retrofit strategies. The necessity and importance of downscaling future weather data for building performance simulations cannot be underestimated when assessing the impact of climate change on building performance. Building performance simulations, abbreviated as BPS, are an essential tool for predicting the thermal behavior, energy consumption, and indoor comfort of buildings under future climate conditions. Accurate future weather data are essential for these simulations as they provide the fine-grained information required for reliable BPS results. Downscaling global climate model outputs to high-resolution local weather data enables architects and engineers to design and retrofit buildings that are resilient, energy-efficient, and comfortable under future climate conditions, thereby promoting sustainable development and climate adaptation. Global climate models are abbreviated as GCMs; downscaling is required because the spatial resolution of GCM outputs is coarse, typically between 100 and 250 kilometers, which is far from sufficient to capture local climate changes and specific weather patterns that affect individual buildings or urban areas. Downscaling methods are used to bridge this gap and transform coarse-grid GCM outputs into high-resolution local weather information required for BPS under future climate conditions. Existing downscaling methods can be divided into two major categories: statistical downscaling and dynamical downscaling. Dynamical downscaling methods use high-resolution regional climate models to simulate local climate processes based on GCM boundary conditions. They provide a detailed physical representation of the climate system and are able to capture the complex interactions between atmospheric, oceanic, and land surface processes. However, it is computationally intensive and requires a lot of resources to run long-term climate simulations. Statistical downscaling methods involve developing empirical relationships between large-scale climate variables and local-scale climate characteristics. They are computationally efficient and can be applied to a variety of GCM outputs, including regression analysis, weather types, and simulated day methods. Most of these techniques assume that the historical relationship between large-scale and local-scale variables will remain unchanged in the future. Due to the computational efficiency and flexibility of statistical downscaling, most studies related to building performance use statistical downscaling instead of dynamical downscaling. Compared with dynamical downscaling methods, statistical downscaling methods require significantly less resources and time, which allows statistical techniques to be applied to multiple GCM outputs and various scenarios, enhancing the robustness of simulations. In addition, statistical methods can use historical weather data such as typical meteorological year data to adjust according to specific local conditions, thus providing more convenient and specific climate predictions. These advantages make statistical downscaling a more practical choice for extensive and detailed building performance simulations. In recent years, many studies have developed and adopted various statistical downscaling methods to improve the accuracy of future local hourly weather data. These methods focus on converting coarse-grid GCM outputs into high-resolution local weather information.However, most statistical downscaling methods are not suitable for application in BPS due to their complexity and computational requirements. The morphing method proposed by the researchers aims to adjust historical weather data to account for predicted future climate change by applying shifts, stretches, and interpolations. It creates future weather files by modifying temperature, humidity, and solar radiation based on GCM projections and typical meteorological year data. The method is simple, efficient, easy to implement, and has low computational requirements, making it suitable for BPS-related applications. Its flexibility allows it to adapt to a variety of climates and scenarios, and integrates well with typical meteorological year data using reliable historical baselines. In addition, the practicality of the method is obvious because it can work without the need for long-term historical data.

[0003] However, the morphing method still has some limitations. It relies on linear assumptions, cannot fully capture local climate patterns, and cannot accurately predict extreme meteorological conditions. Summary of the invention

[0004] In order to overcome the shortcomings of the prior art, the purpose of the present invention is to provide a method for hourly downscaling of monthly meteorological forecast data based on distribution adjustment time mapping, which realizes the hourly downscaling of monthly forecast data of different weather variables in the atmospheric circulation model through the adjustment of specific distribution and quantile mapping method, while maintaining the physical properties of the data and the statistical characteristics of historical data, reflecting the long-term trends and extreme weather conditions predicted by GCM, and providing reliable future meteorological data for building performance simulation.

[0005] To achieve the above object, the present invention provides the following solutions:

[0006] A method for hourly downscaling of monthly meteorological forecast data based on distribution-adjusted time mapping, comprising:

[0007] Identify representative cities for pre-defined climate zones;

[0008] Compile historical weather data based on the representative cities to obtain typical meteorological year data;

[0009] Determine global climate patterns and monthly mean resolution datasets;

[0010] Extract hourly typical meteorological year data of preset weather variables according to the typical meteorological year data to obtain a main variable data set; the main variable data set includes: hourly temperature data, hourly relative humidity data, hourly wind speed data and hourly solar radiation data;

[0011] Setting a candidate statistical distribution of each of the preset weather variables according to the primary variable data set;

[0012] Verifying the candidate statistical distribution of each of the preset weather variables using maximum likelihood estimation to obtain a verification data set;

[0013] Screening the candidate statistical distributions of the preset weather variables according to the validation data set to obtain an optimal distribution of each of the preset weather variables;

[0014] Constructing an empirical cumulative distribution function based on the monthly average resolution data set and the main variable data set;

[0015] Using a preset correction function to correct the optimal distribution to obtain a corrected distribution;

[0016] Determine the quantile corresponding to each value in the prediction result of the global climate model using the empirical cumulative distribution function;

[0017] Mapping the quantile to the calculation result of the empirical cumulative distribution function of the primary variable data set to obtain an original prediction model;

[0018] The prediction model is adjusted using a preset distribution adjustment formula to obtain a target prediction model;

[0019] The target prediction model is used to predict future hourly weather data to obtain a target prediction result.

[0020] Preferably, the screening criteria for the representative cities include: weather patterns, temperature ranges and humidity levels.

[0021] Preferably, the climate scenarios of the global climate model include: SSP126, SSP245, SSP370, SSP434 and SSP585.

[0022] Preferably, historical weather data is compiled based on the representative cities to obtain typical meteorological year data, including:

[0023] Collect raw meteorological data through public weather databases;

[0024] Using a preset program to perform integrity detection and data cleaning on the original meteorological data to obtain preprocessed data;

[0025] Formatting the preprocessed data according to a preset data structure to obtain formatted data;

[0026] The missing values ​​obtained by the integrity detection of the preset program are supplemented by using the difference technology, and the supplementary data of the missing values ​​are matched to the formatted data to obtain the typical meteorological year data.

[0027] Preferably, the candidate statistical distributions include: normal distribution, lognormal distribution and skewed normal distribution.

[0028] Preferably, the judgment data of the maximum likelihood estimation is the residual sum of squares; the calculation formula of the residual sum of squares is: Wherein, SSE is the calculation result of the residual sum of squares; is the i-th observed data point; x i is the i-th fitted distribution value with parameter θ; n is the number of data points.

[0029] Preferably, the empirical cumulative distribution function is: Among them, F X (x is the empirical cumulative distribution function; II() is the indicator function; N is the number of observations; X i The i-th set of meteorological parameters; x is the meteorological parameter value.

[0030] Preferably, the correction function is: Y adjusted is the parameter correction result; is the parameter correction value; μ fut is the average value of the prediction results of the global climate model; μ historical is the mean value of the main variable data set.

[0031] Preferably, the distribution adjustment formula includes:

[0032]

[0033]

[0034] Among them, loc and scale are distribution parameters; a and b are shape parameters; σ obs is the standard deviation of the main variable data set; Γ(·) is the gamma function.

[0035] The present invention discloses the following technical effects:

[0036] The present invention provides an hourly downscaling method for monthly meteorological forecast data based on distribution adjustment time mapping. Through the adjustment of specific distribution and quantile mapping methods, the problem that conventional technologies cannot reflect the long-term trends and extreme weather conditions predicted by GCM while maintaining the physical properties of the data and the statistical characteristics of historical data is solved, and the hourly downscaling of monthly forecast data of different weather variables in the atmospheric circulation model and the reliable simulation of future meteorological data are achieved. BRIEF DESCRIPTION OF THE DRAWINGS

[0037] In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings required for use in the embodiments will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative labor.

[0038] Figure 1 A schematic diagram of a meteorological forecast data downscaling process based on distribution adjustment time mapping provided in an embodiment of the present invention. DETAILED DESCRIPTION

[0039] The following will be combined with the drawings in the embodiments of the present invention to clearly and completely describe the technical solutions in the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0040] The purpose of the present invention is to provide a meteorological forecast data downscaling method based on distribution adjustment time mapping, which realizes the hourly downscaling of monthly forecast data of different weather variables in the atmospheric circulation model through the adjustment of specific distribution and quantile mapping method, while maintaining the physical properties of the data and the statistical characteristics of historical data, reflecting the long-term trends and extreme weather conditions predicted by GCM, and providing reliable future meteorological data for building performance simulation.

[0041] In order to make the above-mentioned objects, features and advantages of the present invention more obvious and easy to understand, the present invention is further described in detail below with reference to the accompanying drawings and specific embodiments.

[0042] Figure 1 A schematic diagram of a downscaling process of meteorological forecast data based on distribution adjustment time mapping provided by an embodiment of the present invention is as follows: Figure 1 As shown, the present invention provides a method for hourly downscaling of monthly meteorological forecast data based on distribution adjustment time mapping, comprising:

[0043] Identify representative cities for pre-defined climate zones;

[0044] Compile historical weather data based on the representative cities to obtain typical meteorological year data;

[0045] Determine global climate patterns and monthly mean resolution datasets;

[0046] Extract hourly typical meteorological year data of preset weather variables according to the typical meteorological year data to obtain a main variable data set; the main variable data set includes: hourly temperature data, hourly relative humidity data, hourly wind speed data and hourly solar radiation data;

[0047] Setting a candidate statistical distribution of each of the preset weather variables according to the primary variable data set;

[0048] Verifying the candidate statistical distribution of each of the preset weather variables using maximum likelihood estimation to obtain a verification data set;

[0049] Screening the candidate statistical distributions of the preset weather variables according to the validation data set to obtain an optimal distribution of each of the preset weather variables;

[0050] Constructing an empirical cumulative distribution function based on the monthly average resolution data set and the main variable data set;

[0051] Using a preset correction function to correct the optimal distribution to obtain a corrected distribution;

[0052] Determine the quantile corresponding to each value in the prediction result of the global climate model using the empirical cumulative distribution function;

[0053] Mapping the quantile to the calculation result of the empirical cumulative distribution function of the primary variable data set to obtain an original prediction model;

[0054] The prediction model is adjusted using a preset distribution adjustment formula to obtain a target prediction model;

[0055] The target prediction model is used to predict future hourly weather data to obtain a target prediction result.

[0056] Optionally, the screening criteria for the representative cities include: weather patterns, temperature ranges and humidity levels.

[0057] Preferably, the climate scenarios of the global climate model include: SSP126, SSP245, SSP370, SSP434 and SSP585.

[0058] Specifically, historical weather data of the representative cities are compiled to obtain typical meteorological year data, including:

[0059] Collect raw meteorological data through public weather databases;

[0060] Using a preset program to perform integrity detection and data cleaning on the original meteorological data to obtain preprocessed data;

[0061] Formatting the preprocessed data according to a preset data structure to obtain formatted data;

[0062] The missing values ​​obtained by the integrity detection of the preset program are supplemented by using the difference technology, and the supplementary data of the missing values ​​are matched to the formatted data to obtain the typical meteorological year data.

[0063] Preferably, the candidate statistical distributions include: normal distribution, lognormal distribution and skewed normal distribution.

[0064] Optionally, the judgment data of the maximum likelihood estimation is the residual sum of squares; the calculation formula of the residual sum of squares is: Wherein, SSE is the calculation result of the residual sum of squares; is the i-th observed data point; x i is the i-th fitted distribution value with parameter θ; n is the number of data points.

[0065] Preferably, the empirical cumulative distribution function is: F X (x); where F X (x) is the empirical cumulative distribution function; II() is the index function; N is the number of observations; X i The i-th set of meteorological parameters; x is the meteorological parameter value.

[0066] Furthermore, the correction function is: adjusted =y i corr +(μ fut -μ historical ), where Y adjusted is the parameter correction result; is the parameter correction value; μ fut is the average value of the prediction results of the global climate model; μ historical is the mean value of the main variable data set.

[0067] Specifically, the distribution adjustment formula includes:

[0068]

[0069] Among them, loc and scale are distribution parameters; a and b are shape parameters; σ obs is the standard deviation of the main variable data set; Γ(·) is the gamma function.

[0070] Specifically, this embodiment collects typical meteorological year (TMY) data for representative cities in different ASHRAE climate zones. TMY data provides a reliable historical baseline of weather conditions, which is crucial for developing accurate downscaling methods. Each city selected in this embodiment represents a specific climate zone, characterized by unique weather patterns, temperature ranges, and humidity levels. This diverse selection ensures that the downscaling method can be tested and verified under a wide range of climatic conditions, thereby enhancing the robustness and applicability of the distribution adjusted time mapping (DATM) method. The TMY data for these cities come from a reputable platform, providing high-quality and consistent historical weather data for developing and verifying downscaling methods. Table 1 lists the selected representative cities and their climate zones, respectively.

[0071] Table 1

[0072]

[0073]

[0074] The list of cities is referenced from a typical building model database to ensure comprehensive representation of different ASHRAE climate zones. By selecting these representative cities to cover a wide range of climate conditions, it is intended to provide a solid basis for the application and validation of the proposed DATM method under different climate conditions.

[0075] Furthermore, the global climate model (GCM) used in this embodiment is MRI-ESM2-0 developed by a meteorological research institute. This model is part of the Coupled Model Intercomparison Project Phase 6 (CMIP6) project and is widely recognized for its comprehensive characterization of the climate system. MRI-ESM2-0 integrates various atmospheric, oceanic and land surface processes to simulate the Earth's climate and predict future climate scenarios based on different greenhouse gas emission trajectories. This embodiment uses the latest GCM output under the Shared Socioeconomic Pathway (SSP). SSP is a set of standardized scenarios that describe different trajectories of future social development, including greenhouse gas emissions, land use and economic growth. The SSPs selected in this embodiment include:

[0076] SSP126: This scenario represents a low greenhouse gas emissions pathway that emphasizes sustainable development and limiting global warming to below 2°C;

[0077] SSP245: Medium emissions scenario, where global development continues along historical trends with some mitigation efforts, leading to medium levels of greenhouse gas concentrations;

[0078] SSP370: a high-emissions scenario characterized by dispersed development and minimal mitigation efforts, leading to large increases in greenhouse gas concentrations and significant global warming;

[0079] SSP434: This scenario represents a pathway with moderate levels of mitigation and adaptation measures, resulting in moderate greenhouse gas emissions and associated climate impacts;

[0080] SSP585: an extreme scenario where emissions are very high, resulting in severe climate impacts, due to rapid economic growth and continued reliance on fossil fuels;

[0081] Monthly outputs of MRI-ESM2-0 under these SSPs are obtained based on selected representative cities. These outputs include key meteorological variables such as temperature, relative humidity, wind speed, and global horizontal irradiance (GHI). In order to combine these GCM projections with local historical weather data, the monthly GCM outputs are downscaled using the proposed Distribution Adjusted Time Mapping (DATM) method. The use of MRI-ESM2-0 and SSP scenarios ensures that the generated future weather data encompass a wide range of potential future climate scenarios, providing a basis for comparing and evaluating the performance of downscaling methods.

[0082] Preferably, this embodiment compiles historical weather data from 2015 to 2023 based on selected representative cities. This period was selected to reflect the real situation of climate conditions affected by ongoing climate change. Historical data can be used as a benchmark for evaluating the accuracy of downscaling methods in capturing real-world climate change and their ability to accurately predict future weather scenarios. Historical weather data comes from publicly accessible weather databases such as the National Oceanic and Atmospheric Administration (NOAA) and other local meteorological agencies, ensuring high quality and consistent records. The collected data include hourly observations of key meteorological variables such as temperature, relative humidity, and wind speed. However, global horizontal irradiance (GHI) data is missing in most sites and is therefore not considered in this embodiment. Data processing is performed using Python 3.9 and involves several key steps:

[0083] 1) Check the raw data for missing values, outliers, and inconsistencies; automated scripts are used to identify and correct errors to ensure the integrity of the dataset;

[0084] 2) formatting the cleaned data to be consistent with the structure of TMY, thus achieving seamless integration with the DATM method;

[0085] 3) For the period of missing data, statistical methods and interpolation techniques are used to accurately estimate the missing values ​​and maintain the continuity and reliability of the data set.

[0086] Preferably, the distribution adjusted time mapping (DATM) method is a method proposed in this embodiment, which provides an alternative method for future weather data prediction by combining monthly GCM output with historical TMY data. The method includes fitting the statistical distribution of historical weather data and adjusting these distributions according to GCM estimates to generate high-resolution hourly future weather data. DATM can solve the downscaling problem of multiple variables. The three weather variables that are mainly reduced and analyzed in this embodiment include temperature, wind speed and relative humidity.

[0087] Further, the DATM method first identifies and parameterizes the statistical distribution of weather data from the TMY dataset. This process includes the following key steps to produce appropriate and correct distributions of weather variables. First, hourly TMY data for each variable, temperature, relative humidity, wind speed, and solar radiation are obtained. The data is then analyzed to find the distribution that best fits the data. Commonly used distributions are normal distribution for temperature, lognormal distribution for wind speed and relative humidity, and skewed or lognormal distribution for solar radiation. In order to determine the most appropriate distribution, the method requires comparing several potential distributions and estimating them on historical data through maximum likelihood estimation (MLE). The goodness of fit of each distribution is evaluated based on statistical measures, such as the sum of squared errors (SSE) between the actual data and the fitted distribution. SSE is calculated using the following formula:

[0088]

[0089] Where x i is the observed data point; For x i The parameter is the fitted distribution value of θ; n is the number of data points.

[0090] Further, the distribution with the minimum SSE is considered the best fit. After selecting the optimal distribution, the parameters of the distribution are estimated so that the mean of the fitted distribution is as close as possible to the mean of the historical data. Some modifications are made to the parameters, especially for the Weibull, Rayleigh, skew-normal and log-normal distributions, to remove any significant mean shift. This is essential for adjusting the historical data to the future climate conditions estimated by the GCM outputs. The DATM method uses certain probability distributions for each climate variable, which are chosen to fit the climate variable in question. Table 2 summarizes the assumed distributions for each variable.

[0091] Table 2

[0092]

[0093] These distributions were chosen for their ability to represent the physical characteristics and statistical behavior of each variable. For temperature, the normal distribution and the distribution were used to accommodate the nature of the temperature data. Relative humidity is the ratio of the amount of water vapor in the air to the maximum amount of water vapor the air can hold at a given temperature, expressed as a percentage, and its range from 0 to 100% is best described by a lognormal distribution or a beta distribution. Wind speed is non-negative and often right-skewed, therefore, it is modeled using a lognormal distribution, a Weibull distribution, or a Rayleigh distribution. For solar radiation, since only positive values ​​are considered, skewed normal distribution, normal distribution, and Rayleigh distribution are used to model the variability of this variable as affected by atmospheric conditions. This helps to maintain the statistical properties of each climate variable during the downscaling process to the greatest extent possible.

[0094] Preferably, the DATM method uses the so-called statistical mapping to match the frequency distribution of historical TMY data with the projected future climate conditions obtained from GCMs. This process makes it possible to obtain downscaled weather data that correctly describe the changes in climate variables that are expected to occur. The method consists of several steps:

[0095] 1) In order to make historical observations consistent with future predictions, the first step is to construct empirical cumulative distribution functions (CDFs). For the hourly resolution TMY data X, construct empirical CDFs F X (x)(i.e. F X (x or CDF), CDF represents the probability distribution of historical TMY data:

[0096]

[0097] 2) A set of candidate distributions is used to determine the best fitting distribution for each climate variable. Each candidate distribution is fitted to the observed data using maximum likelihood estimation (MLE). The best fitting distribution and its parameters are then selected based on the minimum sum of squared errors (SSE). To ensure that the scaled-down future data reflect the mean changes projected by the climate models, the revised values ​​are adjusted to match the projected mean values ​​of the future models. For a given variable Y (i.e., a meteorological variable such as temperature, wind speed, etc.), the revised value is and the average output μ of future GCM fut , the adjustment can be described as:

[0098]

[0099] Y adjusted represents the final adjusted variable value. This is the final result after the correction is completed, reflecting the predicted changes in future climate conditions; represents the corrected value before mean adjustment. This step ensures that the long-term trend of the model forecast is preserved while correcting the distribution. After the best fitting distribution is determined and corrected, the quantiles of the future model forecasts need to be mapped to the historical observations. For the future model M with monthly mean resolution future Each value m in the estimate i , CDFF of typical meteorological year data in TMY m (m i ) to determine the corresponding quantile:

[0100] q i =F m (m i )

[0101] F m Represents the future model prediction value M future The cumulative distribution function (CDF) of , which has a monthly mean resolution. This is similar to the CDF representing the historical TMY (Typical Meteorological Year) data, but F m Specifically refers to the distribution of future predicted values ​​of the model. m (m i ), it describes the use of the CDF to find the future model forecast for a specific value m i The corresponding probability (or quantile) process. This probability q i This is then used to map the Find the corresponding value in the historical distribution. The specific breakdown is as follows:

[0102] F m is the cumulative distribution function of future model data;

[0103] m i It is a specific value from future predictions;

[0104] F m (m i ) gives the probability / quantile q of this value appearing in the future model distribution i ;

[0105] This quantile is then used to map back to the historical distribution.

[0106] This process is a key part of the statistical downscaling method, which aims to adjust future climate projections while maintaining the statistical properties of historical observations.

[0107] Using the quantile, we can map this quantile to the observed history (F X (x) to obtain the reduced hourly value (The function value corresponding to the quantile in the probability distribution adjusted by the mean):

[0108]

[0109] This mapping ensures that the future model data distribution is as close as possible to the historical data distribution. For the DATM method, 95% of the mapped values ​​are contained within the calculated margin. While retaining the historical data characteristics and the predicted future mean, the 5% mapped quantiles are clipped to prevent the weather data from being shrunk beyond realistic limits due to distribution mapping.

[0110] Furthermore, the DATM method employs distribution-specific adjustment techniques in order to map historical data into future predictions based on monthly averaged GCM outputs. In this context, the mathematical expressions for modifying the Rayleigh, Weibull, and Beta distributions are discussed separately. For the commonly used Rayleigh distribution for wind speed, the adjustment function minimizes the mean and standard deviation errors:

[0111]

[0112] where μ fut is the predicted future mean, σ obs is the observed standard deviation, and loc and scale are the distribution parameters.

[0113] The Weibull distribution is also suitable for wind speed, and the scale parameter is optimized to adjust while maintaining the shape:

[0114]

[0115] Where Γ represents the gamma function.

[0116] For the Beta distribution used for relative humidity, the shape parameters a, b were optimized:

[0117]

[0118] These optimization problems are solved numerically using Sequential Least Squares Programming (SLSQP) for Rayleigh distribution, Brent’s method for Weibull distribution, and Limited Memory Broyden-Fletcher-Goldfarb-Shanno (L-BFGS-B) algorithm for Beta distribution. This approach ensures that the adjusted distribution is as close as possible to the expected mean and the variability of historical data. Depending on the distribution, different algorithms are used to numerically solve these optimization problems. For Rayleigh distribution, Sequential Least Squares Programming (SLSQP) is used. SLSQP is an iterative method for solving constrained nonlinear optimization problems. It works by solving a series of quadratic programming subproblems. At each iteration k, the objective function f(u) and the constraint c(u) are approximated using a quadratic model and a linear model, respectively:

[0119]

[0120] Where p is the search direction, H k is an approximation of the Hessian matrix, E and I are the sets of equality constraints and inequality constraints respectively; u k represents the solution vector of the kth iteration; when i belongs to E, it means it is under an equality constraint, and when i belongs to I, it means it is under an inequality constraint. Then the algorithm updates the solution:

[0121] u k+1 =u k +α k ·p k

[0122] where α k Determined by the row search routine.

[0123] For the Weibull distribution, Brunt's method is used. This is an algorithm for finding roots that combines bisection, the secant method, and inverse quadratic interpolation. Brunt's method starts with an interval [a,b] known to contain the root. At each iteration, it chooses between three methods:

[0124] Inverse quadratic interpolation:

[0125]

[0126] Secant method

[0127]

[0128] dichotomy

[0129]

[0130] This method selects the most efficient method among the three methods in each iteration, ensuring the reliability of the bisection method and the fast convergence of the inverse quadratic interpolation as much as possible.

[0131] For the Beta distribution, the L-BFGS-B algorithm is used. It is a quasi-Newton method for solving large-scale optimization problems with simple bounds on the variables. The algorithm approximates the inverse Hessian matrix using limited computer memory. The core idea is to use the following update formula for the inverse Hessian approximation:

[0132]

[0133] but:

[0134]

[0135] s k =u k+1 -u k

[0136]

[0137] The algorithm only stores and uses the last memory size parameter m (a preset value used to approximate the previous iteration information of the Hessian matrix) for (s k ,y k ) to approximate the Hessian, where m is typically between 3 and 20. This allows it to handle large-scale problems efficiently. The L-BFGS-B variant also handles box constraints (i.e., simple bounds on the variables) by using a projected gradient method. At each iteration, it identifies the active set of bounds and performs the optimization in a subspace of the free variables. The algorithm uses this approximation to determine the search direction and performs a line search, similar to other quasi-Newton methods, but with the added ability to handle bounds constraints efficiently. This approach is particularly well suited for optimizing the shape parameter of the Beta distribution while being mindful of the presence of bounds, allowing the distribution to be precisely adjusted to match the predicted future mean and preserving the statistical properties of the historical data. m effectively represents the length of memory or history that the algorithm uses to build its approximation of the inverse Hessian matrix. Larger values ​​of m generally provide better approximations but require more memory and computational effort per iteration, while smaller values ​​of m use less memory but may require more iterations to converge.

[0138] The beneficial effects of the present invention are as follows:

[0139] The present invention realizes hourly downscaling of monthly forecast data of different weather variables in the atmospheric circulation model through the adjustment of specific distribution and quantile mapping method, which not only maintains the physical properties of the data and the statistical characteristics of historical data, but also reflects the long-term trends and extreme weather conditions predicted by GCM, providing reliable future meteorological data for building performance simulation.

[0140] The various embodiments in this specification are described in a progressive manner, and each embodiment focuses on the differences from other embodiments. The same or similar parts between the various embodiments can be referenced to each other.

[0141] The principles and implementation methods of the present invention are described in this article using specific examples. The description of the above embodiments is only used to help understand the method and core idea of ​​the present invention. At the same time, for those skilled in the art, according to the idea of ​​the present invention, there will be changes in the specific implementation methods and application scope. In summary, the content of this specification should not be understood as limiting the present invention.

Claims

1. A method for hourly downscaling of monthly meteorological forecast data based on distribution-adjusted time mapping, characterized in that: include: Identify representative cities for pre-defined climate zones; Compile historical weather data based on the representative cities to obtain typical meteorological year data; Determine global climate patterns and monthly mean resolution datasets; Extract hourly typical meteorological year data of preset weather variables according to the typical meteorological year data to obtain a main variable data set; the main variable data set includes: hourly temperature data, hourly relative humidity data, hourly wind speed data and hourly solar radiation data; Setting a candidate statistical distribution of each of the preset weather variables according to the primary variable data set; Verifying the candidate statistical distribution of each of the preset weather variables using maximum likelihood estimation to obtain a verification data set; Screening the candidate statistical distributions of the preset weather variables according to the validation data set to obtain an optimal distribution of each of the preset weather variables; Constructing an empirical cumulative distribution function based on the monthly average resolution data set and the main variable data set; Using a preset correction function to correct the optimal distribution to obtain a corrected distribution; Determine the quantile corresponding to each value in the prediction result of the global climate model using the empirical cumulative distribution function; Mapping the quantile to the calculation result of the empirical cumulative distribution function of the primary variable data set to obtain an original prediction model; The prediction model is adjusted using a preset distribution adjustment formula to obtain a target prediction model; The target prediction model is used to predict future hourly weather data to obtain a target prediction result.

2. The method for downscaling monthly meteorological forecast data hourly based on distribution-adjusted time mapping according to claim 1, characterized in that: The representative cities were selected based on criteria including weather patterns, temperature ranges and humidity levels.

3. The method for downscaling monthly meteorological forecast data hourly based on distribution-adjusted time mapping according to claim 1, characterized in that: The climate scenarios of the global climate model include: SSP126, SSP245, SSP370, SSP434 and SSP585.

4. The method for downscaling monthly meteorological forecast data hourly based on distribution-adjusted time mapping according to claim 1, characterized in that: Based on the historical weather data of the representative cities, typical meteorological year data are obtained, including: Collect raw meteorological data through public weather databases; Using a preset program to perform integrity detection and data cleaning on the original meteorological data to obtain preprocessed data; Formatting the preprocessed data according to a preset data structure to obtain formatted data; The missing values ​​obtained by the integrity detection of the preset program are supplemented by using the difference technology, and the supplementary data of the missing values ​​are matched to the formatted data to obtain the typical meteorological year data.

5. The method for downscaling monthly meteorological forecast data hourly based on distribution-adjusted time mapping according to claim 1, characterized in that: The candidate statistical distributions include: normal distribution, lognormal distribution and skewed normal distribution.

6. The method for downscaling monthly meteorological forecast data hourly based on distribution-adjusted time mapping according to claim 1, characterized in that: The judgment data of the maximum likelihood estimation is the residual sum of squares; the calculation formula of the residual sum of squares is: Wherein, SSE is the calculation result of the residual sum of squares; is the i-th observed data point; x i is the i-th fitted distribution value with parameter θ; n is the number of data points.

7. The method for downscaling monthly meteorological forecast data hourly based on distribution-adjusted time mapping according to claim 1, characterized in that: The empirical cumulative distribution function is: Among them, F X( x ) is the empirical cumulative distribution function; II() is the indicator function; N is the number of observations; X i The i-th set of meteorological parameters; x is the meteorological parameter value.

8. The method for downscaling monthly meteorological forecast data hourly based on distribution-adjusted time mapping according to claim 1, characterized in that: The correction function is: Among them, Y adjusted is the parameter correction result; is the parameter correction value; μ fut is the average value of the prediction results of the global climate model; μ historical is the mean value of the main variable data set.

9. The method for downscaling monthly meteorological forecast data hourly based on distribution-adjusted time mapping according to claim 1, characterized in that: The distribution adjustment formula includes: Among them, loc and scale are distribution parameters; a and b are shape parameters; σ obs is the standard deviation of the main variable data set; ( · ) is the gamma function.

Images