Evaporation waveguide height determination method based on domestic numerical weather forecast mode

Abstract

The invention relates to the technical field of aerospace meteorology, and particularly discloses a domestic numerical weather forecast mode-based evaporation waveguide height determination method, which comprises the following steps of: obtaining meteorological coupling data output by a domestic numerical weather forecast mode CMA-MESO; determining parameters according to the meteorological coupling data and a turbulent flux algorithm COARE; constructing a temperature vertical profile and a specific humidity vertical profile according to the parameters, the Monin-Obchhoff similarity theory and an NPS model; constructing a modified refractive index profile according to the temperature vertical profile and the specific humidity vertical profile; and determining the height of the evaporation waveguide according to the height corresponding to the lowest point of the corrected refractive index profile. The method is suitable for determining the evaporation waveguide height in a domestic business meteorological mode.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of space weather, and particularly relates to a method for determining evaporation duct height based on a domestic numerical weather prediction model. BACKGROUND

[0002] As a typical low-altitude atmospheric duct structure, the evaporation duct (EDH) has a significant impact on the propagation characteristics of electromagnetic waves in the marine environment. Its existence can significantly enhance the effectiveness of radar, communication, navigation and other systems in the open ocean and sea environment, and has important military and civilian value. However, the evaporation duct height is highly sensitive to changes in meteorological elements, especially dependent on high-resolution, real-time near-sea meteorological data.

[0003] In the open ocean and sea area, conventional observation data is scarce, and international reanalysis data has problems such as update lag, coverage limitations, and poor physical consistency, making it difficult to meet the needs of China's sea area for high-precision and controllable EDH prediction.

[0004] Therefore, it is urgent to develop an EDH calculation method based on the domestic operational numerical prediction system to realize the real combat support and meteorological support for key sea areas. SUMMARY

[0005] To solve the above problems, the purpose of the present application is to provide a method for determining the evaporation duct height based on a domestic numerical weather prediction model, which relies on the CMA-MESO model independently developed by the China Meteorological Administration to realize the determination of the evaporation duct height suitable for the domestic operational meteorological model.

[0006] The present application provides a method for determining the evaporation duct height based on a domestic numerical weather prediction model, comprising: obtaining meteorological coupling data output by a domestic numerical weather prediction model CMA-MESO; determining parameters according to the meteorological coupling data and the COARE turbulent flux algorithm; constructing temperature vertical profiles and specific humidity vertical profiles according to the parameters, the Monin-Obukhov similarity theory and the NPS model; constructing a modified refractive index profile according to the temperature vertical profiles and the specific humidity vertical profiles; determining the evaporation duct height according to the height corresponding to the lowest point of the modified refractive index profile.

[0007] In one possible implementation, the parameters include: Monin-Obukhov length, roughness length, temperature scale parameter and humidity scale parameter.

[0008] In a possible implementation, the determining the parameters according to the meteorological coupling data and the COARE turbulent flux algorithm includes: The Monin-Obukhov length is expressed according to the following formula : The Monin-Obukhov stability parameter is expressed according to the following formula : In the formula, is the von Karman constant, is the acceleration of gravity, is the measurement height of the wind speed, is the absolute air temperature observation value, is the temperature scale parameter, is the humidity scale parameter, is the friction velocity.

[0009] In a possible implementation, the momentum roughness length is expressed according to the following formula : The roughness Reynolds number is expressed according to the following formula : The humidity roughness length is expressed according to the following formula In the formula, is the momentum roughness length, is the wind speed-dependent Charnok parameter, is the acceleration of gravity, is the friction velocity, is the kinematic viscosity of air, is the humidity roughness length, is the heat roughness length.

[0010] In a possible implementation, the friction velocity is expressed according to the following formula : The humidity scale parameter is expressed according to the following formula : The temperature scale parameter is expressed according to the following formula : In the formula, is the friction velocity, The effective wind speed at 10m For momentum exchange coefficient, The difference in specific humidity between the sea surface and the air, Humidity exchange coefficient, The potential temperature difference between the sea surface and the air, The heat exchange coefficient is denoted as . In one possible implementation, the momentum exchange coefficient is calculated according to the following formula. : The humidity exchange coefficient can be calculated using the following formula. : Calculate the heat exchange coefficient using the following formula. : In the formula, Wind speed measurement height The height for measuring relative humidity For measuring atmospheric temperature, For momentum roughness length, For thermal roughness length, For humidity roughness length, To calculate the stability correction function for wind speed, To calculate the temperature stability correction function, is the von Kármán constant.

[0011] In one possible implementation, constructing the temperature vertical profile and the specific humidity vertical profile based on the parameters and the Moning-Obukhov similarity theory includes: The expression for the vertical temperature profile is determined using the following formula. : In the formula, It is the atmospheric temperature at the sea surface. It is the dry adiabatic lapse rate, It is the von Kármán constant, It is the SHEBA stability correction function. MOS scale parameters for potential temperature For the selected height, For temperature roughness length, The parameters are the Moning-Obukhov stability parameters.

[0012] In one possible implementation, the expression for the specific humidity vertical profile is determined according to the following formula. : In the formula, It is the specific humidity of the sea surface. It is the dry adiabatic lapse rate, It is the von Kármán constant, It is the SHEBA stability correction function. For the MOS scale parameters of specific humidity, For the selected height, Humidity roughness length, The parameters are the Moning-Obukhov stability parameters.

[0013] In one possible implementation, the SHEBA stability correction function is expressed by the following formula. : In the formula, , , , For constants The parameters are the Moning-Obukhov stability parameters.

[0014] In one possible implementation, constructing the modified refractive index profile includes: The atmospheric corrected refractive index profile is constructed according to the following formula. : In the formula, For the selected height, For air pressure, For the vertical temperature profile, It is the water vapor pressure; Water vapor pressure is expressed by the following formula. : In the formula, The ratio of the molar mass of water vapor to the average molar mass of dry air. This is the specific humidity vertical profile.

[0015] This invention provides a method for determining evaporation waveguide height based on a domestically developed numerical weather prediction model, achieving independent and domestically developed data-driven and algorithmic mechanisms. It fully utilizes a large amount of regional observation data (including ship, buoy, and remote sensing data) from the CMA-MESO model, significantly improving the simulation accuracy of boundary layer elements. The improved COARE3.0 algorithm and SHEBA stability correction function are introduced to enhance adaptability and accuracy under strongly stable stratification conditions in tropical seas. It supports seamless integration into the CMA-MESO forecast process, enabling hourly, high-resolution dynamic estimation and distribution mapping of EDH. This provides technical support for shipborne communication systems, electromagnetic environment prediction, and marine combat capability assessment, and has broad application prospects. Attached Figure Description

[0016] Figure 1 A flowchart illustrating the method for determining the height of an evaporating waveguide provided in an embodiment of the present invention. Detailed Implementation

[0017] The embodiments of the present invention will be further described in detail below with reference to the accompanying drawings and examples. The following detailed description of the embodiments and the accompanying drawings are used to illustrate the principles of the present invention by way of example, but should not be used to limit the scope of the present invention. That is, the present invention is not limited to the described preferred embodiments, and the scope of the present invention is defined by the claims.

[0018] In the description of this invention, it should be noted that, unless otherwise stated, "a plurality of" means two or more; the terms "first," "second," etc., are used for descriptive purposes only and should not be construed as indicating or implying relative importance; those skilled in the art can understand the specific meaning of the above terms in this invention as appropriate.

[0019] Figure 1 A flowchart illustrating the method for determining the height of an evaporating waveguide provided in an embodiment of the present invention is shown below. Figure 1 As shown, this invention provides a method for determining the evaporation waveguide height based on a domestic numerical weather prediction model, including: Step S1: Obtain meteorological coupling data output by the domestic numerical weather prediction model CMA-MESO; Step S2: Determine the parameters based on meteorological coupling data and the COARE turbulent flux algorithm; In one possible implementation, the parameters include: Moning-Obukhoff length, roughness length, temperature scale parameter, and humidity scale parameter.

[0021] In iterative calculations, the Moning-Obukhov length And its reciprocal dimensionless stability parameter, in order to characterize the core quantity of atmospheric stability, the atmospheric stability is updated at each step with the latest frictional wind speed and flux scale, ensuring that it eventually converges to a self-consistent solution.

[0022] The length of the Moning-Obukhoff curve is expressed by the following formula. : The Moning-Obukhov stability parameter is expressed by the following formula. : In the formula, For von Kármán constant, For gravitational acceleration, For measuring wind speed at height, For absolute air temperature observation values, For temperature scale parameters, For humidity scale parameters, The friction speed is denoted as .

[0023] The appropriate humidity roughness length is determined by the Reynolds number under different flow conditions.

[0024] The momentum roughness length is expressed by the following formula. : The roughness Reynolds number is expressed by the following formula. : The temperature and humidity roughness length is expressed by the following formula: In the formula, For momentum roughness length, For wind speed-dependent Chanok parameters, For gravitational acceleration, Friction speed, For the kinematic viscosity of air, For humidity roughness length, This represents the thermal roughness length.

[0025] We construct stable modified flux coefficients (momentum exchange coefficient, humidity exchange coefficient, and heat exchange coefficient) to correct the systematic influence of non-neutral stratification (stable / unstable) on turbulent fluxes. We then perform physical process corrections under different conditions of unstable or stable stratification.

[0026] Calculate the momentum transfer coefficient using the following formula. : The humidity exchange coefficient can be calculated using the following formula. : Calculate the heat exchange coefficient using the following formula. : In the formula, Wind speed measurement height The height for measuring relative humidity For measuring atmospheric temperature, For momentum roughness length, For thermal roughness length, For humidity roughness length, To calculate the stability correction function for wind speed, To calculate the temperature stability correction function, is the von Kármán constant.

[0027] The stability correction function is expressed by the following formula. : Iterative updates are performed using a stability correction function to achieve multi-process coupling, reduce errors, and iteratively calculate the updated turbulence scaling parameters, as follows: The friction speed is expressed by the following formula. : The humidity scale parameter is expressed by the following formula. : The temperature scale parameter is expressed by the following formula. : In the formula, Friction speed, The effective wind speed at 10m For momentum exchange coefficient, The difference in specific humidity between the sea surface and the air, Humidity exchange coefficient, The potential temperature difference between the sea surface and the air, The heat exchange coefficient is denoted as . In one possible implementation, when calculating the Moning-Obukhov length... At that time, the initial value of the friction velocity is calculated according to the following formula. The initial value of the friction speed Input the NPS model to obtain the friction velocity. .

[0028] This applies when wind speeds are very low, but factors such as sea surface temperature differences can still drive buoyancy flux. An empirical gust speed is introduced into the wind speed calculation. To address the numerical instability in air-sea flux calculations under low wind speed conditions, a corrected effective wind speed is obtained. The velocity of ocean currents at the sea surface is expressed by the following formula. : In the formula, The wind speed is the value observed at a height of 10m. For gust speed, The effective wind speed is at 10m.

[0029] Establishing a 10-m wind speed conversion is a fundamental step in accurately estimating air-sea flux, ensuring the physical consistency of the parameterization scheme and the universality of the calculation results. The corrected wind speed extrapolated to 10m is expressed by the following formula. : In the formula, The altitude at which wind speed is measured.

[0030] To provide a physically reasonable starting point for the entire parameterization scheme and avoid iterative divergence caused by random initial values, an initial frictional wind speed is established. The frictional speed is expressed by the following formula. : In one possible implementation, the kinematic viscosity of the air needs to be calculated using the following formula when calculating the humidity roughness length. Gravitational acceleration .

[0031] The surface roughness length includes the Charnock term plus a viscosity term. The Charnock term reflects the impact of wave-breaking turbulent disturbances on the air-sea interface, while the viscosity term characterizes molecular viscosity effects, describing the laminar subsurface characteristics under smooth sea surface conditions. It integrates two effects: dynamic roughness driven by wind-generated waves (the first term) and smooth flow dominated by air viscosity (the second term). It provides a unified characterization of wave-induced turbulence and molecular viscosity effects, offering a continuous description from smooth to rough sea states. The kinematic viscosity of air is expressed by the following formula. : The acceleration due to gravity is expressed by the following formula. : In the formula, The selected geographic latitude, in degrees. The gravitational acceleration at the equator. , , , is a coefficient.

[0032] . In one possible implementation, when calculating the humidity roughness length, the initial value of the roughness length is calculated according to the following formula. The initial value of the roughness length Input the NPS model to obtain the roughness length.

[0033] The aerodynamic roughness length is expressed by the following formula. : In the formula, These are atmospheric temperature observations. This is an initial estimate of the length of the sea surface aerodynamic roughness. It is the acceleration due to gravity. The kinetic viscosity of air.

[0034] In one possible implementation, when calculating the stable corrected flux coefficient, the overall exchange coefficient is calculated according to the following formula, and the overall exchange coefficient is... Input the NPS model to obtain the stable corrected flux coefficient.

[0035] The dimensionless drag coefficient / drag coefficient essentially quantifies the momentum exchange efficiency at the air-sea interface. Its main function is to connect wind speed observations with turbulent flux, integrating roughness and stability effects to drive ocean-atmospheric boundary layer processes. The neutral drag coefficient at a 10-meter reference height is expressed by the following formula. : In the formula, For von Kármán's constant, =0.4.

[0036] Estimating the dimensionless scalar exchange coefficient correlation coefficient for heat flux conversion mainly involves converting the dimensionless coefficient into an initial value for the thermal roughness length. This is expressed by the following formula: : The initial value of the thermal roughness length is expressed by the following formula. : In the formula, Neutral heat exchange coefficient at a reference height of 10 meters. .

[0037] The dimensionless flux coefficient for altitude measurement is used to correct for the vertical variation of physical quantities with altitude, eliminating systematic biases caused by differences in observed altitude. It also incorporates physical constraints to ensure that flux calculations satisfy boundary layer dynamics and maintain physical consistency between sensible and latent heat fluxes. For atmospheric turbulence, it quantifies the impact of atmospheric stability on flux transport, distinguishing the inhibitory or enhancing effects of unstable / neutral / stable stratification on turbulent exchange.

[0038] The wind speed is expressed as the measurement height according to the following formula. The momentum exchange coefficient on ; The following formula represents the temperature measurement height. scalar heat exchange coefficient : The comprehensive exchange coefficient is expressed by the following formula. : In the formula, , Used for subsequent calculations of initial stability parameters.

[0039] In one possible implementation, when calculating the Moning-Obukhov length... At this time, the initial value of the Moning-Obukhov length needs to be calculated according to the following formula. The initial value of the Moning-Obukhov length. Input the NPS model to obtain the Moning-Obukhov length. .

[0040] The Moning-Obukhov dimensionless height is used to determine the physical rationality of using a correction function to dynamically adjust flux estimates by quantitatively characterizing the stable state of the atmospheric boundary layer.

[0041] Richardson Number The relative importance of buoyancy and wind shear was quantified. The Richardson number is expressed by the following formula. : Using Richardson's number To estimate the dimensionless Monin-Obukhov stability parameter. The Monin-Obukhov stability parameter is expressed by the following formula. : The Moning-Obukhov length is a core physical scale describing stability. The initial value of the Moning-Obukhov length is expressed by the following formula. : In one possible implementation, the humidity scale parameter is calculated. At that time, the specific humidity difference between the sea surface and the air is calculated using the following formula. .

[0042] The following formula can automatically select between water and ice surfaces based on whether the temperature is below the freezing point. The following formula represents the saturated vapor pressure. : The freezing point of seawater is expressed by the following formula. : The saturated specific humidity is expressed by the following formula. : In the formula, Salinity correction factor, This refers to sea surface salinity.

[0043] The actual specific humidity is expressed by the following formula. : In the formula, For specific humidity observation values, This represents the potential temperature difference between the sea surface and the air. The specific humidity difference between the sea surface and the air is represented by T, the atmospheric temperature is the observed value, and the sea level pressure is the observed value.

[0044] In one possible implementation, the temperature-scale parameters are calculated. At that time, the potential temperature difference between the sea surface and the air is calculated using the following formula. .

[0045] By introducing the dry adiabatic lapse rate to correct the air temperature at the measured altitude, "false" temperature differences caused by changes in altitude and air pressure can be eliminated.

[0046] The lapse rate of dry adiabatic heat is expressed by the following formula. : In the formula, The specific heat constant at constant pressure of dry air. =1004.67 J / (kg·K); The potential temperature difference between the sea surface and the air can be expressed by the following formula. : In the formula, These are seawater temperature observation values. For atmospheric temperature, For the overall Richardson number, To measure the height of wind speed The Moning-Obukhov stability parameters at the location For the length of Mourning-Obukhov, This represents the observed absolute air temperature.

[0047] In one possible implementation, the initial value of the temperature-scale parameter is calculated according to the following formula: Temperature-scale parameter The initial values ​​of the initial values ​​of the humidity scale parameters. ,Will and Input the NPS model to obtain temperature-scale parameters and humidity-scale parameters. .

[0048] The initial scale parameter provides a physically reasonable initial value for turbulent flux iteration, preventing mathematical singularities under low wind speed conditions and suppressing computational oscillations under strong stable stratification.

[0049] The initial value of the temperature scale parameter is expressed by the following formula: The initial value of the humidity scale parameter is expressed by the following formula: Used to correct deviations in logarithmic profiles under non-neutral atmospheric conditions, the integral stability function of heat and water vapor is expressed by the following formula. : In one possible implementation, the momentum roughness length is calculated. At that time, calculate the wind speed-dependent Chanok parameters according to the following formula. .

[0050] The gust factor is expressed by the following formula. : The neutral wind speed at a reference height of 10 meters is expressed by the following formula. : The Channock parameter, which depends on wind speed, is expressed by the following formula. : In the formula, For empirical slope coefficients, These are empirical intercept coefficients.

[0051] Step S3: Construct the temperature vertical profile and the specific humidity vertical profile based on the parameters, the Moning-Obukhov similarity theory, and the NPS model. The expression for the vertical temperature profile is determined using the following formula. : In the formula, It is the atmospheric temperature at the sea surface. It is the dry adiabatic lapse rate, It is the von Kármán constant, It is the SHEBA stability correction function. MOS scale parameters for potential temperature For the selected height, For temperature roughness length, The parameters are the Moning-Obukhov stability parameters.

[0052] The expression for the specific humidity vertical profile is determined using the following formula. : In the formula, It is the specific humidity of the sea surface. It is the dry adiabatic lapse rate, It is the von Kármán constant, It is the SHEBA stability correction function. For the MOS scale parameters of specific humidity, For the selected height, Temperature roughness length, The parameters are the Moning-Obukhov stability parameters.

[0053] The SHEBA stability correction function is a correction function derived from the surface thermal budget of the Arctic marine experiment. It shows applicability under both weak and stable conditions and is widely used to calculate turbulent fluxes.

[0054] The SHEBA stability correction function can be expressed by the following formula. : In the formula, , , , For constants The parameters are the Moning-Obukhov stability parameters.

[0055] In the formula, = =5、 =3、 .

[0056] Step S4: Construct a modified refractive index profile based on the temperature vertical profile and the specific humidity vertical profile; The atmospheric corrected refractive index profile is constructed according to the following formula. : In the formula, For the selected height, For air pressure, For the vertical temperature profile, It is the water vapor pressure; Based on Dalton's law of partial pressures and the ideal gas law, moist air is considered an ideal gaseous mixture of dry air and water vapor. Specific humidity... It is the ratio of the mass of water vapor to the total mass of moist air, while water vapor pressure... It is the partial pressure produced by water vapor as a component of the gas mixture. The following formula serves as a bridge connecting these two different methods of measurement.

[0057] Water vapor pressure is expressed by the following formula. : In the formula, The ratio of the molar mass of water vapor to the average molar mass of dry air. This is the specific humidity vertical profile.

[0058] Among them, air pressure The barometric altimeter formula is calculated according to the following steps. It is the integral form of the hydrostatic equilibrium equation assuming a constant air layer temperature. It quantitatively describes the physical law governing the change in air pressure with geopotential height for an ideal gas in a gravitational field under hydrostatic equilibrium. To extend the hydrostatic equation, which applies only to dry air, to real moist air, an imaginary temperature is introduced. This allows the density change caused by water vapor to be equivalently incorporated into the temperature term, thus maintaining the consistent form of the barometric altimeter formula while using the gas constant for dry air.

[0059] The mean virtual temperature is expressed by the following formula. : (32) The target height is expressed by the following formula. air pressure : in, For reference height, =0、 For target height, The gas constant for dry air =287J / (kg·K) The sea level pressure observation value at an altitude of 0m. For sea surface specific humidity observation value, The temperature reading is taken at a height of 2m.

[0060] Step S5: Determine the height of the evaporation waveguide based on the height corresponding to the lowest point of the modified refractive index profile.

[0061] In one example, the vertical gradient of the corrected refractive index profile The capture height is 14.7m, and the evaporation waveguide height EDH is determined to be 14.7m.

[0062] The key protection points of this invention are as follows: 1) Combine the CMA-MESO forecast field with the COARE3.0 + SHEBA corrected physics algorithm; 2) Utilize the M-profile calculation strategy optimized by physical consistency; 3) A high-precision EDH estimation process that can be embedded into domestic business systems has been implemented.

[0063] Existing methods for determining EDH have the following drawbacks: 1) It heavily relies on foreign reanalysis data and algorithm models, lacking autonomy and controllability; 2) It is difficult to integrate into the domestic operational numerical weather prediction process, and it is impossible to achieve near real-time EDH forecasting; 3) The physical parameters of the model are poorly adapted to the characteristics of the boundary layer in the Chinese sea area, resulting in a large accuracy deviation. In particular, under strong stable stratification or weak disturbance conditions, EDH overestimation is likely to occur. 4) It is difficult to support the environmental adaptability design of domestic shipborne systems and radar electromagnetic propagation performance.

[0064] Existing methods for determining EDH have the following drawbacks: 1) It heavily relies on foreign reanalysis data and algorithm models, lacking autonomy and controllability; 2) It is difficult to integrate into the domestic operational numerical weather prediction process, and it is impossible to achieve near real-time EDH forecasting; 3) The physical parameters of the model are poorly adapted to the characteristics of the boundary layer in the Chinese sea area, resulting in a large accuracy deviation. In particular, under strong stable stratification or weak disturbance conditions, EDH overestimation is likely to occur. 4) It is difficult to support the environmental adaptability design of domestic shipborne systems and radar electromagnetic propagation performance.

[0065] Compared with the prior art, the present invention has the following advantages: 1) Significantly improves the accuracy of EDH calculations: especially in strongly stable marine stratification environments, it overcomes the overestimation problem of traditional models; 2) Supports seamless integration with domestic models: does not rely on foreign data, enhancing the autonomy and controllability of business systems; 3) High-resolution computing power: Supports EDH spatiotemporal distribution mapping with hourly and kilometer-level accuracy; 4) High scalability: It is easy to extend to multiple application scenarios such as radar wave propagation simulation and shipborne system optimization; 5) Excellent environmental adaptability: It is particularly suitable for waveguide calculations in complex tropical marine environments.

[0066] The above description is merely a specific embodiment of the present invention, but the scope of protection of the present invention is not limited thereto. Any variations or substitutions that can be easily conceived by those skilled in the art within the technical scope disclosed in the present invention should be included within the scope of protection of the present invention. Therefore, the scope of protection of the present invention should be determined by the scope of the claims.

Claims

1. A method for determining the height of an evaporation waveguide based on a domestic numerical weather prediction model, characterized in that, include: Acquire meteorological coupling data output by the domestic numerical weather prediction model CMA-MESO; The parameters are determined based on the meteorological coupling data and the COARE turbulent flux algorithm. Based on the parameters, the Moning-Obukhoff similarity theory, and the NPS model, construct the temperature vertical profile and the specific humidity vertical profile; A modified refractive index profile is constructed based on the temperature vertical profile and the specific humidity vertical profile. The height of the evaporation waveguide is determined based on the height corresponding to the lowest point of the modified refractive index profile.

2. The method for determining the height of an evaporating waveguide according to claim 1, characterized in that, The parameters include: Moning-Obukhoff length, roughness length, temperature scale parameter, and humidity scale parameter.

3. The method for determining the height of an evaporating waveguide according to claim 2, characterized in that, The parameters determined based on the meteorological coupling data and the COARE turbulent flux algorithm include: The length of the Moning-Obukhoff curve is expressed by the following formula. : ； The Moning-Obukhov stability parameter is expressed by the following formula. : ； In the formula, For von Kármán constant, For gravitational acceleration, For measuring wind speed at height, For absolute air temperature observation values, For temperature scale parameters, For humidity scale parameters, The friction speed is denoted as .

4. The method for determining the height of an evaporating waveguide according to claim 2, characterized in that, Also includes: The momentum roughness length is expressed by the following formula. : ； The roughness Reynolds number is expressed by the following formula. : ； The humidity roughness length is expressed by the following formula: ； In the formula, For momentum roughness length, For wind speed-dependent Chanok parameters, For gravitational acceleration, Friction speed, For the kinematic viscosity of air, For humidity roughness length, This represents the thermal roughness length.

5. The method for determining the height of an evaporating waveguide according to claim 2, characterized in that, Also includes: The friction speed is expressed by the following formula. : ； The humidity scale parameter is expressed by the following formula. : ； The temperature scale parameter is expressed by the following formula. : ； In the formula, Friction speed, The effective wind speed at 10m For momentum exchange coefficient, The difference in specific humidity between the sea surface and the air, Humidity exchange coefficient, The potential temperature difference between the sea surface and the air, The heat exchange coefficient is denoted as .

6. The method for determining the height of an evaporating waveguide according to claim 6, characterized in that, Also includes: Calculate the momentum transfer coefficient using the following formula. : ； The humidity exchange coefficient can be calculated using the following formula. : ； Calculate the heat exchange coefficient using the following formula. : ； In the formula, Wind speed measurement height The height for measuring relative humidity For measuring atmospheric temperature, For momentum roughness length, For thermal roughness length, For humidity roughness length, To calculate the stability correction function for wind speed, To calculate the temperature stability correction function, is the von Kármán constant.

7. The method for determining the height of an evaporating waveguide according to claim 1, characterized in that, The construction of the temperature vertical profile and the specific humidity vertical profile based on the parameters, the Moning-Obukhov similarity theory, and the NPS model includes: The expression for the vertical temperature profile is determined using the following formula. : ； In the formula, It is the atmospheric temperature at the sea surface. It is the dry adiabatic lapse rate, It is the von Kármán constant, It is the SHEBA stability correction function. MOS scale parameters for potential temperature For the selected height, For temperature roughness length, The parameters are the Moning-Obukhov stability parameters.

8. The method for determining the height of an evaporating waveguide according to claim 7, characterized in that, Also includes: The expression for the specific humidity vertical profile is determined using the following formula. : ； In the formula, It is the specific humidity of the sea surface. It is the dry adiabatic lapse rate, It is the von Kármán constant, It is the SHEBA stability correction function. For the MOS scale parameters of specific humidity, For the selected height, Roughness length due to humidity The parameters are the Moning-Obukhov stability parameters.

9. The method for determining the height of an evaporating waveguide according to claim 8, characterized in that, Also includes: The SHEBA stability correction function can be expressed by the following formula. : ； In the formula, , , , For constants The parameters are the Moning-Obukhov stability parameters.

10. The method for determining the height of an evaporating waveguide according to claim 1, characterized in that, The construction of the modified refractive index profile includes: The atmospheric corrected refractive index profile is constructed according to the following formula. : ； In the formula, For the selected height, For air pressure, For the vertical temperature profile, It is the water vapor pressure; Water vapor pressure is expressed by the following formula. : ； In the formula, The ratio of the molar mass of water vapor to the average molar mass of dry air. This is the specific humidity vertical profile.

Images