Regional Weather Modification Sounding Data Processing Method

Abstract

The embodiment of the present invention relates to a method for processing regional sounding data for artificial weather modification, comprising: acquiring sounding data of multiple temperature layers; the sounding data comprising temperature, air pressure and dew point temperature; determining multiple air pressure standard layer information and temperature layer information according to the temperature and air pressure; each air pressure standard layer information corresponds to multiple temperatures and dew point temperatures, and each temperature layer information corresponds to multiple air pressure and dew point temperatures; calculating a first physical quantity set according to the air pressure standard layer information and the temperature layer information; calculating a K index, a modified K index, a total index, a modified total index and an A index; determining estimated weather information according to the K index, the modified K index, the total index, the modified total index and the A index; and determining target weather information according to the estimated weather information and the first physical quantity set.

Description

Technical Field

[0001] The invention relates to the field of meteorological technology, and in particular to a method for processing regional artificial weather influencing sounding data. Background Art

[0002] Meteorological sounding data in atmospheric engineering is an important basis for studying and evaluating the atmospheric environment, and the processing of meteorological sounding data directly affects the accuracy and reliability of the analysis results. In the existing technology, the data of sounding data often exist in isolation and cannot be well integrated, which leads to a limited application scope of sounding data. Summary of the invention

[0003] The purpose of the present invention is to provide a method for processing regional artificial weather modification sounding data in view of the defects of the prior art, so as to solve the problems existing in the prior art.

[0004] To achieve the above object, an embodiment of the present invention provides a method for processing regional artificial weather modification sounding data, the method comprising:

[0005] Acquiring sounding data of multiple temperature layers; the sounding data includes temperature, air pressure and dew point temperature;

[0006] According to the temperature and air pressure, multiple air pressure standard layer information and temperature layer information are determined; each air pressure standard layer information corresponds to multiple temperatures and dew point temperatures, and each temperature layer information corresponds to multiple air pressures and dew point temperatures;

[0007] Calculating a first physical quantity set according to the pressure standard layer information and the temperature layer information;

[0008] According to K = t 850 -t 500 +t d850 -t 700 +t d700 Calculate the K index; where K is the correction index, t 860 is the temperature at 850hPa, t 500 is the temperature at 500hPa, t d850 is the dew point temperature at 850hPa, t 700 is the temperature at 700hPa, t d700 The dew point temperature is 700hPa;

[0009] According to MK=t 地面 +t d850 -t d700 -t 500 Calculate the modified K index; where MK is the modified K index, t 地面 is the ground temperature;

[0010] According to TI=(t850 -t 500 )+(t d850 -t d500 ) calculates the total index; where TI is the total index;

[0011] According to tt=0.5*(t 地面 +t d地面 +t 850 +t d850 )-2t 500 Calculate the modified total index; where TT is the modified total index, t 地面 is the ground dew point temperature;

[0012] According to A=(t 850 -t 500 )-[(tt d ) 850 +(tt d ) 700 +(tt d ) 500 ]Calculate the A index; where A is the A index, tt d is the difference between the temperature and the dew point temperature;

[0013] Determining estimated weather information based on the K index, the modified K index, the total index, the modified total index and the A index;

[0014] Target weather information is determined according to the estimated weather information and the first set of physical quantities.

[0015] In a possible implementation, the first set of physical quantities includes saturated water vapor pressure, saturated specific humidity, mixing ratio and saturated mixing ratio, relative humidity, condensation temperature, specific humidity, humid air gas constant, humid air constant pressure specific volume and humid air constant volume specific volume, humid air potential temperature, virtual potential temperature, virtual temperature, vertical pressure gradient, condensation latent heat, wet static temperature, wet static temperature, wet bulb temperature, wet bulb potential temperature; the second set of physical quantities includes a first index, a first correction index, a second index, a second correction index and a third index.

[0016] In a possible implementation manner, determining target weather information according to the estimated weather information and the first set of physical quantities specifically includes:

[0017] According to the estimated weather information, a first condition factor and a second condition factor are selected from the first physical quantity set; the correlation degree between the first condition factor and the meteorological element is greater than the correlation degree between the second condition factor and the meteorological element;

[0018] Target weather information is determined according to the first condition factor and the second condition factor.

[0019] In a possible implementation manner, determining the meteorological condition corresponding to the target meteorological element according to the first condition factor and the second condition factor specifically includes:

[0020] determining a first identifier according to a first difference between the first condition factor and a first threshold;

[0021] determining a second identifier according to a second difference between the second condition factor and a second threshold;

[0022] Target weather information is determined according to a minimum value of the first identifier and the second identifier.

[0023] In a possible implementation manner, the calculating the first physical quantity set according to the pressure standard layer information and the temperature layer information specifically includes:

[0024] Using dew point temperature, calculate water vapor pressure and saturated water vapor pressure at different levels;

[0025] Calculating saturated specific humidity, mixing ratio and saturated mixing ratio at different levels by using the saturated water vapor pressure and the air pressure;

[0026] Calculating relative humidity at different levels using the mixing ratio and the saturated mixing ratio;

[0027] Calculating condensation temperatures at different levels using the temperature and the relative humidity;

[0028] Calculating specific humidity at different levels using the water vapor pressure and air pressure;

[0029] Using specific humidity, calculate the gas constant of moist air at different levels, the specific volume of moist air at constant pressure, and the specific volume of moist air at constant volume;

[0030] Calculate the moist air potential temperature using temperature, moist air gas constant, and moist air constant pressure specific volume;

[0031] Calculate virtual potential temperatures at different levels using specific humidity and the potential temperature of the humid air;

[0032] Using the temperature and the specific humidity, calculating virtual temperatures at different levels;

[0033] Calculating vertical pressure gradients at different levels using air pressure and the virtual temperature;

[0034] Using the temperatures, calculate the latent heat of condensation at different levels:

[0035] Calculate the wet static temperature at different levels using absolute temperature, potential height, specific volume of moist air at constant pressure, latent heat of condensation, and specific humidity;

[0036] Calculate the saturated wet static temperature at different levels using absolute temperature, potential height, constant pressure specific volume of moist air, latent heat of condensation, and saturated specific humidity;

[0037] Calculate the wet bulb temperature at different levels using temperature, specific humidity, saturated specific humidity, and specific volume of moist air at constant pressure;

[0038] Calculate the wet bulb potential temperature at different levels by using the wet bulb temperature, air pressure, humid air gas constant, and humid air constant pressure specific volume;

[0039] Use temperature, air pressure and water vapor pressure to calculate the dry air potential temperature at different levels.

[0040] In a possible implementation, the calculation of water vapor pressure and saturated water vapor pressure at different levels by using dew point temperature specifically includes:

[0041] According to the formula Calculate the water vapor pressure at different levels; e is the water vapor pressure, td is the dew point temperature; when td>-15℃, e is the water vapor pressure on the water surface, when td<-15℃, e is the water vapor pressure on the ice surface; a is the coefficient;

[0042] According to the formula Calculate the saturated water vapor pressure at different levels; where e s is the saturated water vapor pressure, when td>-15℃ a is the water vapor pressure on the water surface, when td<-15℃, e s is the water vapor pressure on the ice surface.

[0043] In a possible implementation, the calculating of the saturated specific humidity, the mixing ratio and the saturated mixing ratio at different levels by using the saturated water vapor pressure and the air pressure specifically includes:

[0044] According to the formula Calculate the saturated specific humidity, where q s is the saturated specific humidity; e s is the saturated water vapor pressure; P is the air pressure;

[0045] According to the formula Calculate the mixing ratio; where r is the mixing ratio and e is the water vapor pressure;

[0046] According to the formula Calculate the saturated mixing ratio, where r s is the saturated mixing ratio.

[0047] In a possible implementation, the using the mixing ratio and the saturated mixing ratio to calculate the relative humidity at different levels, using the temperature and the relative humidity to calculate the condensation temperature at different levels; and using the water vapor pressure and the air pressure to calculate the specific humidity at different levels specifically include:

[0048] According to the formula Calculate the relative humidity at different levels; where r is the mixing ratio, r s is the saturated mixing ratio, U is the relative humidity;

[0049] According to the formula Calculate the condensation temperature; where T L is the condensation temperature, T k is the absolute temperature corresponding to the temperature;

[0050] According to the formula Calculate the specific humidity; q is the specific humidity and e is the water vapor pressure.

[0051] In a possible implementation, the use of specific humidity to calculate the gas constant of humid air, the specific volume of humid air at constant pressure, and the specific volume of humid air at constant volume at different levels specifically includes:

[0052] According to the formula RR = R d (1+0.608*q) calculates the humid air gas constant; where RR is the humid air gas constant, q is the specific humidity, and R d is a constant;

[0053] According to formula C p =C pd *(1+0.86*q) calculates the constant pressure specific volume of moist air, where C p is the specific volume of moist air at constant pressure, C pd is a constant;

[0054] According to formula C v =C vd (1+0.96q) to calculate the specific volume of moist air at constant volume; where C v is the constant volume of moist air.

[0055] In a possible implementation, the calculating of virtual potential temperatures at different levels by using the specific humidity and the potential temperature of the humid air specifically includes:

[0056] According to θ v =(1+0.608*q)θ to calculate the virtual temperature; where θ v is the virtual potential temperature, θ is the moist air potential temperature, and q is the specific humidity.

[0057] By applying the regional artificial weather modification sounding data processing method provided by the embodiment of the present invention, more parameters can be obtained through the sounding data, and through further fusion of various parameters, the application range of the parameters is further expanded, the utilization rate of the sounding data is improved, and meteorological analysis can be better performed. BRIEF DESCRIPTION OF THE DRAWINGS

[0058] Figure 1A flow chart of a method for processing regional weather modification sounding data provided by an embodiment of the present invention;

[0059] Figure 2 for Figure 1 Flow chart of step 120 in FIG. DETAILED DESCRIPTION

[0060] In order to make the purpose, technical scheme and advantages of the present invention clearer, the present invention will be further described in detail below with reference to the accompanying drawings. Obviously, the described embodiments are only part of the embodiments of the present invention, rather than all 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.

[0061] The technical solution of the present invention is further described in detail below through the accompanying drawings and embodiments.

[0062] Figure 1 The flowchart of the method for processing regional artificial weather modification sounding data provided by the embodiment of the present invention is as follows. Figure 1 , the technical solution of the present invention is described with specific embodiments.

[0063] Step 100, sounding data of multiple temperature layers; the sounding data includes temperature, air pressure and dew point temperature;

[0064] Specifically, the temperature layers include 0°C, -5°C, -10°C, -15°C, -20°C, and -30°C. The sounding station includes multiple standard layers including ground, 850hPa, 700hPa, 500hPa, 400hPa, 300hPa, 250hPa, 200hPa, 150hPa, and 100hPa.

[0065] If the sounding data is in MICAPS Category 5 format, it is necessary to interpolate based on the standard layer and characteristic layer data to automatically calculate the air pressure, geopotential height, and dew point temperature of the 0℃, -5℃, -10℃, -15℃, -20℃, and -30℃ layers. Alternatively, draw a temperature logarithmic pressure diagram (T-LogP diagram) of a sounding station on the interface, and automatically determine the air pressure and geopotential height of the 0℃, -5℃, -10℃, -15℃, -20℃, and -30℃ layers based on the point coordinates, and the determination results are automatically saved to the system. Jump to step 4 to continue.

[0066] If the sounding data is L-band sounding second data, the air pressure and potential height of the 0℃, -5℃, -10℃, -15℃, -20℃ and -30℃ layers can be read directly.

[0067] Step 110, determining multiple air pressure standard layer information and temperature layer information according to the temperature and air pressure; each air pressure standard layer information corresponds to multiple temperatures and dew point temperatures, and each temperature layer information corresponds to multiple air pressures and dew point temperatures;

[0068] Step 120, calculating and obtaining a first physical quantity set according to the pressure standard layer information and the temperature layer information;

[0069] Among them, the first set of physical quantities includes saturated water vapor pressure, saturated specific humidity, mixing ratio and saturated mixing ratio, relative humidity, condensation temperature, specific humidity, humid air gas constant, humid air constant pressure specific volume and humid air constant volume specific volume, humid air potential temperature, virtual potential temperature, virtual temperature, vertical pressure gradient, condensation latent heat, wet static temperature, wet static temperature, wet bulb temperature, wet bulb potential temperature. Figure 2 As shown, according to the pressure standard layer information and the temperature layer information, the first physical quantity set is calculated to specifically include:

[0070] Step 201, using the dew point temperature, calculate the water vapor pressure and saturated water vapor pressure at different levels;

[0071] Step 202, using the saturated water vapor pressure and the air pressure, calculating the saturated specific humidity, the mixing ratio and the saturated mixing ratio at different levels;

[0072] Step 203, using the mixing ratio and the saturated mixing ratio, calculate the relative humidity at different levels;

[0073] Step 204, using the temperature and relative humidity, calculating the condensation temperature of different layers;

[0074] Step 205, using the water vapor pressure and the air pressure, calculate the specific humidity at different levels;

[0075] Step 206, using the specific humidity, calculating the gas constant of humid air, the specific volume of humid air at constant pressure and the specific volume of humid air at constant volume at different levels;

[0076] Step 207, calculating the humid air potential temperature using the temperature, the humid air gas constant, and the humid air constant pressure specific volume;

[0077] Step 208, using the specific humidity and the potential temperature of the humid air, calculate the virtual potential temperature at different levels;

[0078] Step 209, using temperature and specific humidity, calculating virtual temperatures at different levels;

[0079] Step 210, using air pressure and virtual temperature, calculate the vertical air pressure gradients at different levels;

[0080] Step 211, using the temperature, calculate the condensation latent heat at different levels:

[0081] Step 212, calculating the wet static temperature at different levels using absolute temperature, potential height, specific volume of moist air at constant pressure, latent heat of condensation, and specific humidity;

[0082] Step 213, calculating the saturated wet static temperature at different levels using absolute temperature, potential height, constant pressure specific volume of moist air, latent heat of condensation, and saturated specific humidity;

[0083] Step 214, using temperature, specific humidity, saturated specific humidity, and constant pressure specific volume of humid air to calculate wet bulb temperatures at different levels;

[0084] Step 215, using the wet bulb temperature, air pressure, humid air gas constant, and humid air constant pressure specific volume, calculate the wet bulb potential temperature at different levels;

[0085] Step 216, using temperature, air pressure, and water vapor pressure, calculate the dry air potential temperature at different levels.

[0086] The following formula is used to specifically describe steps 201 to 215:

[0087] 1) Using dew point temperature t d Calculate the water vapor pressure e at different levels according to the formula:

[0088]

[0089] Where e is the water vapor pressure, t d is the dew point temperature; when t d When td>-15℃, e is the water vapor pressure on the water surface, and when td<-15℃, e is the water vapor pressure on the ice surface; a is the coefficient, and b is the coefficient;

[0090] When t d >-15℃ e is the water vapor pressure on the water surface, that is, e=e 水 (a=17.2693882, b=35.86);

[0091] When t d When <-15℃, e is the water vapor pressure on the ice surface, that is, e=e 冰 (a=21.8745584, b=7.66);

[0092] 2) Calculate the specific humidity q at different levels using the water vapor pressure e and air pressure P according to the formula:

[0093]

[0094] 3) Use the temperature t to calculate the saturated water vapor pressure e at different levels according to the formula s :

[0095]

[0096] When td >-15℃ s is the water vapor pressure on the water surface, that is, e s =e s水 (a=17.2693882b=35.86)

[0097] When t d <-15℃ s is the water vapor pressure on the ice surface, that is, e s =e s冰 (a=21.8745584b=7.66)

[0098] 4) Using saturated water vapor pressure e s , air pressure P is used to calculate the saturated relative humidity qs at different levels according to the formula:

[0099]

[0100] 5) Use the water vapor pressure e and air pressure P to calculate the mixing ratio r of different layers according to the formula:

[0101]

[0102] 6) Using saturated water vapor pressure e s , air pressure P calculates the saturated mixing ratio r at different levels according to the formula s :

[0103]

[0104] 7) Using the mixing ratio r and the saturation mixing ratio r s Calculate the relative humidity uw at different levels:

[0105]

[0106] 8) Use the relative humidity q to calculate the gas constant RR of humid air at different levels according to the formula:

[0107] RR=R d (1+0.608*q)

[0108] Among them, R d =207.05J / kgK

[0109] 9) Using the specific humidity q, calculate the constant pressure specific volume C of humid air at different levels according to the formula p :

[0110] C p =C pd *(1+0.86*q)

[0111] C pd =1004J / kg K

[0112] 10) Using the specific humidity q, calculate the constant volume specific volume C of humid air at different levels according to the formula v :

[0113] C v =C vd (1+0.96q)

[0114] 11) Using temperature t and relative humidity q, calculate virtual temperature T at different levels according to the formula v :

[0115] T v =(1+0.608*q)(t+273.16)-273.16

[0116] 12) Calculate the latent heat of condensation at different levels using the temperature t according to the formula:

[0117] L = (2500-2.4*t) × 10 3

[0118] 13) Using temperature t, humid air gas constant RR, humid air constant pressure specific volume C p , calculate the moist air potential temperature θ (in Kelvin) according to the formula:

[0119]

[0120] 14) Calculate the dry air potential temperature θ at different levels using the formula: temperature t, air pressure P, and water vapor pressure e d (Units are in Kelvin):

[0121]

[0122] 15) Using the relative humidity q and the potential temperature of wet air θ, calculate the virtual potential temperature θ at different levels according to the formula V (Units are in Kelvin):

[0123] θ V =(1+0.608*q)θ

[0124] 16) Using temperature t, specific humidity q, saturated specific humidity q s 、Specific volume of moist air at constant pressure C p , calculate the wet bulb temperature Tw (in Kelvin) at different levels according to the formula:

[0125]

[0126] 17) Using wet bulb temperature T w , air pressure P, humid air gas constant RR, humid air constant pressure specific volume C p , calculate the wet bulb potential temperature θ at different levels according to the formula ω(Units are in Kelvin):

[0127]

[0128] 18) Using absolute temperature T (T = t + 273.15 (K)), potential height Z, constant pressure specific volume of moist air Cp, latent heat L, and specific humidity q, the wet static temperature T at different levels is calculated according to the formula σ :

[0129]

[0130] Where g is the acceleration due to gravity, g = 9.80665 m / s2

[0131] 19) Using absolute temperature T (T = t + 273.15 (K)), potential height Z, constant pressure specific volume of moist air Cp, latent heat L, saturated specific humidity q s , calculate the saturated wet static temperature at different levels according to the formula

[0132]

[0133] 20) Using the air pressure P and virtual temperature Tv, calculate the vertical pressure gradient G at different levels according to the formula z :

[0134]

[0135] Among them, gravitational acceleration g = 9.80665m / s2, R d =287.05J / kg K.

[0136] 21) Use temperature t (converted to absolute temperature) and relative humidity U to calculate the condensation temperature TL (i.e. the temperature at the elevated condensation height) at different levels:

[0137]

[0138] Among them, T k , U are the absolute temperatures of the air block in the initial state (T k =t+273.15(K)), relative humidity (i.e. the ratio of water vapor pressure to saturated water vapor pressure: U=e / e s *100), T L The unit is Kelvin (K).

[0139] Step 130, according to K=t 850 -t 500 +td 850 -t 700 +t d700 Calculate the K index; where K is the correction index, t850 is the temperature at 850hPa, t 500 is the temperature at 500hPa, t d850 is the dew point temperature at 850hPa, t 700 is the temperature at 700hPa, t d50 0 is the dew point temperature of 700hPa;

[0140] Specifically, the temperature t at 850 hPa and 700 hPa, the dew point temperature t d , and the temperature t at 500hPa, the K index is calculated according to the formula. The K index can reflect the stratification stability of the atmosphere.

[0141] Step 140, according to MK=t 地面 +t d850 -t d700 -t 500 Calculate the modified K index; where MK is the modified K index, t 地面 is the ground temperature;

[0142] Specifically, the modified K index MK is calculated according to the formula using the ground temperature t at 500 hPa and the dew point temperature td at 850 hPa and 700 hPa.

[0143] Step 150, according to TI=(t 850 -t 500 )+(t d850 -t d500 ) calculates the total index; where TI is the total index;

[0144] Specifically, the total index TI is calculated according to the formula using the temperature t at 850 hPa and 500 hPa and the dew point temperature td.

[0145] Step 160, according to TT = 0.5*(t 地面 +t d地面 +t 850 +t d850 )-2t 500 Calculate the modified total index; where TT is the modified total index, t 地面 is the ground dew point temperature;

[0146] Specifically, the corrected total index TT is calculated according to the formula using the temperatures t of the ground, 850 hPa, and 500 hPa, and the dew point temperatures td of the ground and 850 hPa.

[0147] Step 170, according to A=(t 850 -t 500 )-[(tt d ) 850 +(ttd ) 700 +(tt d ) 500 ]Calculate the A index; where A is the A index, tt d is the difference between the temperature and the dew point temperature;

[0148] Specifically, the temperature t of 850 hPa and 500 hPa, and the dew point temperature difference (tt) of 850 hPa, 700 hPa, and 500 hPa are used. d ), calculate the A index A according to the formula. Index A represents the weather conditions.

[0149] Step 180, determining estimated weather information according to the K index, the modified K index, the total index, the modified total index and the A index;

[0150] Specifically, weather information can be estimated based on the five indexes. For example, any number of the five indexes can be combined to obtain estimated weather information, and the estimated weather information can be displayed in the form of a dynamic graph for easy viewing.

[0151] Step 190, determining target weather information based on the estimated weather information and the first physical quantity set.

[0152] Specifically, according to the estimated weather information, a first condition factor and a second condition factor are selected from the first physical quantity set; the correlation degree between the first condition factor and the rain meteorological element is greater than the correlation degree between the second condition factor and the meteorological element; and the target weather information is determined according to the first condition factor and the second condition factor.

[0153] Determining the meteorological conditions corresponding to the target meteorological element according to the first condition factor and the second condition factor specifically includes: determining the first identifier according to the first difference between the first condition factor and the first threshold; determining the second identifier according to the second difference between the second condition factor and the second threshold; and determining the target weather information according to the minimum value of the first identifier and the second identifier.

[0154] Among them, the first condition factor is at least one of saturated water vapor pressure, saturated specific humidity, mixing ratio and saturated mixing ratio, relative humidity, condensation temperature, specific humidity, wet air gas constant, wet air constant pressure specific volume and wet air constant volume specific volume, wet air potential temperature, virtual potential temperature, virtual temperature, vertical pressure gradient, condensation latent heat, wet static temperature, wet static temperature, wet bulb temperature, and wet bulb potential temperature; the second condition factor is at least one of saturated water vapor pressure, saturated specific humidity, mixing ratio and saturated mixing ratio, relative humidity, condensation temperature, specific humidity, wet air gas constant, wet air constant pressure specific volume and wet air constant volume specific volume, wet air potential temperature, virtual potential temperature, virtual temperature, vertical pressure gradient, condensation latent heat, wet static temperature, wet static temperature, wet bulb temperature, and wet bulb potential temperature, and the second condition factor is different from the first condition factor.

[0155] In one example, the first condition factor may be a main factor, and the step of determining the main factor in the present application includes:

[0156] (1) Read the saturated wet static temperature at 400 hPa and the wet static temperature at 850 hPa, calculate the difference, and determine whether the difference is less than or equal to the threshold. If so, record the identifier X′1 as 1, otherwise record the identifier X′1 as 0;

[0157] (2) Read the vertical pressure gradient of 400 hPa and the ground, calculate the difference, and determine whether the difference is greater than or equal to the threshold. If so, record the identifier X′2 as 1, otherwise record the identifier X′2 as 0.

[0158] 2) The second condition factor can be judged as follows:

[0159] (1) Read the water vapor pressure at 700 hPa and the saturated water vapor pressure, calculate the difference, and determine whether the difference is less than or equal to the threshold. If so, record the identifier X″1 as 1; otherwise, record the identifier X″1 as 0;

[0160] (2) Read the dew point temperature and saturated wet static temperature at 700 hPa, calculate the difference, and determine whether the difference is less than or equal to a threshold value. If so, record the identifier X″2 as 1; otherwise, record the identifier X″2 as 0.

[0161] 3) Comprehensive judgment based on the first condition factor and the second condition factor:

[0162] (1) Find the minimum value of identifiers X′1 and X′2;

[0163] (2) Find the maximum value of identifiers X″1 and X″2;

[0164] (3) Find the minimum value of the output results of (1) and (2). If the result is 1, it can be determined that there is a hail warning in the area covered by station xx, and continue to step 7. If the result is 0, record the log of no hail warning for this sounding, and end this calculation.

[0165] By applying the regional artificial weather modification sounding data processing method provided by the embodiment of the present invention, more parameters can be obtained through the sounding data, and through further fusion of various parameters, the application range of the parameters is further expanded, the utilization rate of the sounding data is improved, and meteorological analysis can be better performed.

[0166] The professionals should further realize that the units and algorithm steps of each example described in conjunction with the embodiments disclosed herein can be implemented by electronic hardware, computer software, or a combination of the two. In order to clearly illustrate the interchangeability of hardware and software, the composition and steps of each example have been generally described in the above description according to the function. Whether these functions are performed in hardware or software depends on the specific application and design constraints of the technical solution. Professionals and technicians can use different methods to implement the described functions for each specific application, but such implementation should not be considered to be beyond the scope of the present invention.

[0167] The steps of the method or algorithm described in conjunction with the embodiments disclosed herein may be implemented using hardware, a software module executed by a processor, or a combination of the two. The software module may be placed in a random access memory (RAM), a memory, a read-only memory (ROM), an electrically programmable ROM, an electrically erasable programmable ROM, a register, a hard disk, a removable disk, a CD-ROM, or any other form of storage medium known in the art.

[0168] The specific implementation methods described above further illustrate the objectives, technical solutions and beneficial effects of the present invention in detail. It should be understood that the above description is only a specific implementation method of the present invention and is not intended to limit the scope of protection of the present invention. Any modifications, equivalent substitutions, improvements, etc. made within the spirit and principles of the present invention should be included in the scope of protection of the present invention.

Claims

1. A method for processing regional weather modification sounding data, characterized in that: The method comprises: Acquiring sounding data of multiple temperature layers; the sounding data includes temperature, air pressure and dew point temperature; According to the temperature and air pressure, multiple air pressure standard layer information and temperature layer information are determined; each air pressure standard layer information corresponds to multiple temperatures and dew point temperatures, and each temperature layer information corresponds to multiple air pressures and dew point temperatures; Calculating a first physical quantity set according to the pressure standard layer information and the temperature layer information; according to Calculate the K index; where K is the correction index, The temperature is 850hPa, is a temperature of 500hPa, The dew point temperature is 850hPa. The temperature is 700hPa, The dew point temperature is 700hPa; according to Calculate the modified K index; where, To correct the K index, is the ground temperature; according to Calculate the total index; where, is the total index; according to Calculate the modified total index; where, To correct the overall index, is the ground dew point temperature; according to Calculate the A index; where A is the A index, tt d is the difference between the temperature and the dew point temperature; Determining estimated weather information based on the K index, the modified K index, the total index, the modified total index and the A index; Determining target weather information according to the estimated weather information and the first set of physical quantities; According to the estimated weather information, a first condition factor and a second condition factor are selected from the first physical quantity set; the correlation degree between the first condition factor and the meteorological element is greater than the correlation degree between the second condition factor and the meteorological element; determining target weather information according to the first condition factor and the second condition factor; The determining of the target weather information according to the first condition factor and the second condition factor specifically includes: The first conditional factor is the main factor. The steps to identify the main factor include: Read the saturated wet static temperature of 400hPa and the wet static temperature of 850hPa, calculate the difference, and determine whether the difference is less than or equal to the threshold. If so, record the identifier 1, otherwise record the identifier is 0; Read the vertical pressure gradient of 400hPa and the ground, calculate the difference, and determine whether the difference is greater than or equal to the threshold. If so, record the identifier 1, otherwise record the identifier is 0; The second condition factor is judged as follows: Read the water vapor pressure at 700hPa and the saturated water vapor pressure, calculate the difference, and determine whether the difference is less than or equal to the threshold. If so, record the identifier 1, otherwise record the identifier is 0; Read the dew point temperature and saturated wet static temperature at 700hPa, calculate the difference, and determine whether the difference is less than or equal to the threshold. If so, record the identifier. 1, otherwise record the identifier is 0; Comprehensive judgment based on the first condition factor and the second condition factor includes: Pair Identifier , Find the minimum value; Pair Identifier , Find the maximum value; Pair Identifier , Minimum value and identifier , Continue to find the minimum value of the maximum value. If the result is 1, it is judged that there is a hail warning.

2. The method according to claim 1, characterized in that The first set of physical quantities includes saturated water vapor pressure, saturated specific humidity, mixing ratio and saturated mixing ratio, relative humidity, condensation temperature, specific humidity, humid air gas constant, humid air constant pressure specific volume and humid air constant volume specific volume, humid air potential temperature, virtual potential temperature, virtual temperature, vertical pressure gradient, condensation latent heat, wet static temperature, wet static temperature, wet bulb temperature, and wet bulb potential temperature.

3. The method according to claim 1, characterized in that The calculating the first physical quantity set according to the pressure standard layer information and the temperature layer information specifically includes: Using dew point temperature, calculate water vapor pressure and saturated water vapor pressure at different levels; Calculating saturated specific humidity, mixing ratio and saturated mixing ratio at different levels by using the saturated water vapor pressure and the air pressure; Calculating relative humidity at different levels using the mixing ratio and the saturated mixing ratio; Calculating condensation temperatures at different levels using the temperature and the relative humidity; Calculating specific humidity at different levels using the water vapor pressure and air pressure; Using specific humidity, calculate the gas constant of moist air at different levels, the specific volume of moist air at constant pressure, and the specific volume of moist air at constant volume; Calculate the moist air potential temperature using temperature, moist air gas constant, and moist air constant pressure specific volume; Calculate virtual potential temperatures at different levels using specific humidity and the potential temperature of the humid air; Using the temperature and the specific humidity, calculating virtual temperatures at different levels; Calculating vertical pressure gradients at different levels using air pressure and the virtual temperature; Using the temperatures, calculate the latent heat of condensation at different levels: Calculate the wet static temperature at different levels using absolute temperature, potential height, specific volume of moist air at constant pressure, latent heat of condensation, and specific humidity; Calculate the saturated wet static temperature at different levels using absolute temperature, potential height, constant pressure specific volume of moist air, latent heat of condensation, and saturated specific humidity; Calculate the wet bulb temperature at different levels using temperature, specific humidity, saturated specific humidity, and specific volume of moist air at constant pressure; Calculate the wet bulb potential temperature at different levels by using the wet bulb temperature, air pressure, humid air gas constant, and humid air constant pressure specific volume; Use temperature, air pressure and water vapor pressure to calculate the dry air potential temperature at different levels.

4. The method according to claim 3, characterized in that The calculation of the saturated specific humidity, mixing ratio and saturated mixing ratio at different levels by using the saturated water vapor pressure and the air pressure specifically includes: According to the formula Calculate the saturated specific humidity, where q s is the saturated specific humidity; e s is the saturated water vapor pressure; P is the air pressure; According to the formula Calculate the mixing ratio; where r is the mixing ratio and e is the water vapor pressure; According to the formula Calculate the saturated mixing ratio, where is the saturated mixing ratio.

5. The method according to claim 3, characterized in that: The method of calculating the relative humidity at different levels by using the mixing ratio and the saturated mixing ratio, calculating the condensation temperature at different levels by using the temperature and the relative humidity, and calculating the specific humidity at different levels by using the water vapor pressure and the air pressure specifically includes: According to the formula Calculate the relative humidity at different levels; where r is the mixing ratio, r s is the saturated mixing ratio, U is the relative humidity; According to the formula Calculate the condensation temperature; where T L is the condensation temperature, T k is the absolute temperature corresponding to the temperature; According to the formula Calculate the specific humidity; q is the specific humidity, e is the water vapor pressure, and p is the air pressure.

6. The method according to claim 3, characterized in that The method of calculating the gas constant of humid air, the specific volume of humid air at constant pressure and the specific volume of humid air at constant volume at different levels by using specific humidity specifically includes: According to the formula Calculate the humid air gas constant; where RR is the humid air gas constant, q is the specific humidity, is a constant; According to the formula Calculate the specific volume of moist air at constant pressure, where is the specific volume of moist air at constant pressure, is a constant.

7. The method according to claim 3, characterized in that The method of calculating virtual potential temperatures at different levels by using specific humidity and the potential temperature of the humid air specifically includes: according to Calculate the virtual temperature; where, For virtual temperature, is the potential temperature of moist air and q is the specific humidity.

Images