CN113609639B - Evaporation waveguide correction method suitable for stable condition - Google Patents

Evaporation waveguide correction method suitable for stable condition Download PDF

Info

Publication number
CN113609639B
CN113609639B CN202110710507.6A CN202110710507A CN113609639B CN 113609639 B CN113609639 B CN 113609639B CN 202110710507 A CN202110710507 A CN 202110710507A CN 113609639 B CN113609639 B CN 113609639B
Authority
CN
China
Prior art keywords
temperature
wind speed
evaporation waveguide
characteristic parameter
specific humidity
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202110710507.6A
Other languages
Chinese (zh)
Other versions
CN113609639A (en
Inventor
王娜
金振中
张炳文
岳明桥
周玢
孙秉志
朱麟
王红光
张利军
慈旋
杨继坤
刘成
章晗
寇雨筱
隋景辉
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
General Research Institute Of Experimental Training Of 92493 Army Of Chinese Pla
Original Assignee
General Research Institute Of Experimental Training Of 92493 Army Of Chinese Pla
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by General Research Institute Of Experimental Training Of 92493 Army Of Chinese Pla filed Critical General Research Institute Of Experimental Training Of 92493 Army Of Chinese Pla
Priority to CN202110710507.6A priority Critical patent/CN113609639B/en
Publication of CN113609639A publication Critical patent/CN113609639A/en
Application granted granted Critical
Publication of CN113609639B publication Critical patent/CN113609639B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/15—Correlation function computation including computation of convolution operations
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/10—Numerical modelling
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/08—Thermal analysis or thermal optimisation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Mathematical Physics (AREA)
  • Algebra (AREA)
  • Software Systems (AREA)
  • Databases & Information Systems (AREA)
  • Computing Systems (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • Waveguides (AREA)

Abstract

The invention provides an evaporation waveguide correction model suitable for a stable condition, which improves the overall performance of the evaporation waveguide model by improving the prediction precision of the evaporation waveguide under the offshore stable condition. Introducing a new pervasive function into an evaporation waveguide model, correcting a wind speed pervasive function in a momentum flux relation under a stable condition in the existing evaporation waveguide model, and correcting a temperature and specific humidity pervasive function in a sensible heat and latent heat flux relation under a stable condition in the existing evaporation waveguide model; under the unstable condition, the wind speed, temperature and specific humidity pervasive function in the flux relation verified by the data of the offshore test is used.

Description

Evaporation waveguide correction method suitable for stable condition
Technical Field
The invention relates to evaporation waveguide in the field of radio meteorology, in particular to a model for predicting evaporation waveguide parameters based on sea surface hydrometeorology parameters, which is mainly used for evaluating the over-the-horizon performance of an offshore microwave band radio system.
Background
The atmospheric waveguide is an abnormal refraction structure which is generated in troposphere atmosphere, can change the normal propagation characteristic of electromagnetic waves, greatly influences the performance of radar, communication, electronic countermeasure and other electronic equipment which work in the marine environment, and particularly can enable the radar to realize over-the-horizon detection and the like.
1) The atmospheric waveguide is divided into an evaporation waveguide, a surface waveguide and a suspension waveguide. The evaporation waveguide is formed by the sharp decrease of humidity along with the height caused by the evaporation of sea water in the sea surface, and is permanently existed in almost any sea area, so that the evaporation waveguide is greatly concerned. With the development of micrometeorology, based on the similar theory of Monenin, different scholars at home and abroad successively put forward different evaporation waveguide models, such as Jesses model, Paulus-Jesses model, MGB model, BYC model, NPS model, pseudo refractive index model and the like. The method is characterized in that a stable condition is set when the temperature difference between the sea and the gas sea is greater than 0, an unstable condition is set when the temperature difference between the gas sea is less than 0, and a neutral condition is set when the temperature difference between the gas sea and the gas sea is equal to 0. In the existing models, evaporation waveguide models such as Jesses models, Paulus-Jesses models, MGB models, pseudo refractive index models and the like are obtained by directly constructing a universal function established by using land test data, and the large model error exists under stable and unstable conditions. The BYC model and the NPS model are obtained based on tropical sea area sea air coupling response test data, and the sea air coupling response test is located in the tropical sea area and mostly under an unstable condition, so that the NPS model can be ensured to have good performance under the unstable condition, but the situation of larger prediction error is probably caused by adopting a universal function established by land test data when the model is established under the stable condition. With the continuous accumulation of the offshore test data, the related scholars establish a universal function with wider application range under stable conditions, and the application range of the flux profile relation in the similar theory of the Morinin is expanded. As in 2007, a new pervasive function under stable conditions was developed based on the SHEBA (Surface Heat Budget of the arc Ocean Surface equilibrium test) test. Therefore, a modified evaporation waveguide model can be constructed by using the new universal functions.
The related patents are as follows:
the invention patent (a regional evaporation waveguide forecasting method, patent number: ZL201218007283.4) provides a method for realizing numerical forecasting of an offshore evaporation waveguide by using a numerical forecasting and evaporation waveguide mode coupling method, so as to expand the functions of a numerical forecasting mode and the evaporation waveguide and solve the regional forecasting problem of the evaporation waveguide.
The invention patent (method for monitoring all types of atmospheric waveguides by ground-based GNSS masker, ZL201518010851.X) provides a method for monitoring atmospheric waveguides by ground-based GNSS masker, which expands received masker signals from GPS to GNSS signals comprising Beidou, GPS and GLONASS on one hand; on the other hand, the evaporation waveguide monitoring method based on ground meteorological data is integrated, and passive monitoring of the full-type atmospheric waveguide is achieved. Overcomes the defects in the prior art.
Disclosure of Invention
The invention provides an evaporation waveguide correction model, and aims to improve the prediction precision of an evaporation waveguide under the marine stable condition, so that the overall higher prediction performance of the evaporation waveguide model is realized.
In order to achieve the purpose, the invention provides a pervasive function established based on measured data under the marine stable condition in an evaporation waveguide model, corrects a wind speed pervasive function in a momentum flux relation under the stable condition in the existing evaporation waveguide model, and corrects a temperature and specific humidity pervasive function in a sensible heat and latent heat flux relation under the stable condition in the existing evaporation waveguide model; and the wind speed, temperature and specific humidity pervasive functions in the flux relation verified by marine test data are used under the unstable condition. The modified evaporation waveguide model comprises the following specific steps:
1) input parameters of evaporation waveguide model
Firstly, inputting parameters of atmospheric temperature Ta, seawater surface temperature Ts, relative humidity RH, pressure P, wind speed u, temperature sensor height zt, humidity sensor height zq and wind speed sensor height zu;
2) parameter initialization
Secondly, initializing the temperature roughness z based on the sea-air coupling relation0tSpecific humidity roughness z0qCharacteristic parameter z of wind speed roughness0Temperature characteristic parameter theta*Specific humidity characteristic parameter q*Characteristic parameter u of wind speed*；
3) Calculation of wind speed pervasive function under unstable condition
Calculating the dimensionless quantity xi ═ zu/L again, namely the stability parameter ═ the erection height of the wind speed sensor/similar length of the tannin; if the value is greater than 0, the wind speed pervasive function takes the form:
Figure GDA0003503269230000031
4) calculation of unsteady condition temperature and specific humidity pervasive functions
The temperature and specific humidity universal functions take the form:
Figure GDA0003503269230000032
wherein, am＝5，bm＝am/6.5，ah＝bh＝5，ch＝3，x＝(1+ξ)13
Figure GDA0003503269230000033
5) Calculation of steady-condition pervasive functions
The pervasive function under the stable condition takes the same form as that in the NPS model, namely the wind speed pervasive function takes the following form:
Figure GDA0003503269230000041
Figure GDA0003503269230000042
Figure GDA0003503269230000043
the temperature and specific humidity universal functions take the form:
Figure GDA0003503269230000044
ψtk＝2ln[(1+zpt)/2],zpt＝(1-15ξ)0.5
Figure GDA0003503269230000045
6) recalculating feature parameters
Recalculating the temperature characteristic parameter theta*Specific dampness characteristic of ginsengNumber q of*Characteristic parameter u of wind speed*Temperature roughness z0tSpecific humidity roughness z0qCharacteristic parameter z of wind speed roughness0。
7) Calculating a feature parameter convergence value
The calculation of the bit temperature characteristic parameter theta is repeated for 3 times again*Specific humidity characteristic parameter q*Characteristic parameter u of wind speed*Temperature roughness z0tSpecific humidity roughness z0qCharacteristic parameter z of wind speed roughness0(ii) a And obtaining the convergence values of the characteristic parameters.
8) Calculating evaporation waveguide modified refractive index profile and height
Then establishing a temperature profile and a specific humidity profile by using the characteristic parameters, and further obtaining a corrected refractive index profile by combining a pressure profile; and finally, selecting the position corresponding to the minimum value of the corrected refractive index as the height of the evaporation waveguide.
The invention has the beneficial effects that: compared with the existing evaporation waveguide model, the method can build the corrected evaporation waveguide model by more universal functions established by marine actual measurement data, achieves higher prediction effect under stable and unstable conditions, and finally can enable the marine evaporation waveguide model result to be more reasonable and credible.
Drawings
FIG. 1 is a flow chart of evaporative waveguide model modification construction.
FIG. 2 shows the calculation results of the height of the evaporation waveguide at 25 ℃ on the sea surface using the NPS and the correction model, respectively, (a) the NPS model, and (b) the correction model.
FIG. 3 shows the calculation results of the height of the evaporation waveguide at a sea surface temperature of 15 ℃ using the NPS and the correction model, respectively, (a) the NPS model, and (b) the correction model.
Detailed Description
The following further describes a specific embodiment of the present invention with reference to the drawings and technical solutions.
Example 1:
the surface temperature of input seawater is 25 ℃, the relative humidity is 80%, the temperature difference between air and sea is 0-3 ℃, the wind speed is 1-16 m/s, and the results of an NPS model and a corrected evaporation waveguide model are shown in an attached figure 2.
As can be seen from the results in FIG. 2, the NPS results show that the height of the evaporation waveguide is abnormally large under the stable condition, such as the conditions of small wind speeds of 1m/s and 4 m/s. However, the height of the evaporation waveguide is generally considered to be within 0-40 m, and under the same input condition, although the result of the model after correction exceeds 40m, the model is still improved and is relatively more reasonable compared with the NPS model.
Example 2:
the surface temperature of the simulated input seawater is respectively 15 ℃, the relative humidity is 80%, the temperature difference of the air and the sea is 0-3 ℃, the wind speed is 1-16 m/s, and the results of an NPS model and a corrected evaporation waveguide model are shown in an attached figure 3.
As can be seen from the results in fig. 3, under a stable condition, the NPS result may have a condition that the height of the evaporation waveguide is abnormally large, such as a condition that the wind speed is 1m/s and a condition that the wind speed is 4m/s (the temperature difference between the air and the sea is about 2 ℃), and under the same input condition, the result of the model after the correction is within 0-40 m, which is relatively more reasonable.

Claims (1)

1. An evaporation waveguide correction method suitable for a stable condition is characterized in that a pervasive function established based on measured data under a marine stable condition is used for correcting a wind speed pervasive function in a momentum flux relation under the stable condition in an existing evaporation waveguide model, and correcting a temperature and specific humidity pervasive function in a sensible heat and latent heat flux relation under the stable condition in an existing evaporation waveguide model; the general function of wind speed, temperature and specific humidity in the flux relation verified by marine test data is used under the unstable condition; the method is characterized by comprising the following specific steps:
1) input parameters of evaporation waveguide model
Inputting an atmospheric temperature Ta, a seawater surface temperature Ts, relative humidity RH, pressure P, a wind speed u, a temperature sensor height zt, a humidity sensor height zq and a wind speed sensor height zu into an evaporation waveguide model;
2) parameter initialization
Initializing temperature roughness z based on sea-air coupling relation0tSpecific humidity roughness z0qWind speed roughness characteristicsParameter z0Temperature characteristic parameter theta*Specific humidity characteristic parameter q*And wind speed characteristic parameter u*；
3) Calculation of wind speed pervasive function under unstable condition
And calculating a dimensionless quantity xi, zu/L, and if the value is larger than 0, adopting the following form as a wind speed universal function:
Figure FDA0003503269220000011
4) calculation of unsteady condition temperature and specific humidity pervasive functions
The temperature and specific humidity universal functions take the form:
Figure FDA0003503269220000012
wherein, am＝5，bm＝am/6.5，ah＝bh＝5，ch＝3，x＝(1+ξ)1/3
Figure FDA0003503269220000013
5) Calculation of steady-condition pervasive functions
The pervasive function under the stable condition takes the same form as that in the NPS model, namely the wind speed pervasive function takes the following form:
Figure FDA0003503269220000021
Figure FDA0003503269220000022
zpu＝(1-15ξ)0.25
Figure FDA0003503269220000023
zpgu＝(1-10.15ξ)0.333，
the temperature and specific humidity universal functions take the form:
Figure FDA0003503269220000024
ψtk＝2ln[(1+zpt)/2],zpt＝(1-15ξ)0.5
Figure FDA0003503269220000025
6) recalculating feature parameters
Recalculating the temperature characteristic parameter theta*Specific humidity characteristic parameter q*Characteristic parameter u of wind speed*Temperature roughness z0tSpecific humidity roughness z0qCharacteristic parameter z of wind speed roughness0；
7) Calculating a feature parameter convergence value
The calculation of the bit temperature characteristic parameter theta is repeated for 3 times again*Specific humidity characteristic parameter q*Characteristic parameter u of wind speed*Temperature roughness z0tSpecific humidity roughness z0qCharacteristic parameter z of wind speed roughness0Obtaining the convergence values of the characteristic parameters;
8) calculating evaporation waveguide modified refractive index profile and height
Then establishing a temperature profile and a specific humidity profile by using the characteristic parameters, and further obtaining a corrected refractive index profile by combining a pressure profile; and finally, selecting the position corresponding to the minimum value of the corrected refractive index as the height of the evaporation waveguide.
CN202110710507.6A 2021-06-25 2021-06-25 Evaporation waveguide correction method suitable for stable condition Active CN113609639B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110710507.6A CN113609639B (en) 2021-06-25 2021-06-25 Evaporation waveguide correction method suitable for stable condition

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110710507.6A CN113609639B (en) 2021-06-25 2021-06-25 Evaporation waveguide correction method suitable for stable condition

Publications (2)

Publication Number Publication Date
CN113609639A CN113609639A (en) 2021-11-05
CN113609639B true CN113609639B (en) 2022-03-25

Family

ID=78336791

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110710507.6A Active CN113609639B (en) 2021-06-25 2021-06-25 Evaporation waveguide correction method suitable for stable condition

Country Status (1)

Country Link
CN (1) CN113609639B (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116610902A (en) * 2023-05-06 2023-08-18 中国人民解放军国防科技大学 Sea area evaporation waveguide daily change influence factor calculation method and system
CN117421601B (en) * 2023-12-19 2024-03-01 山东省科学院海洋仪器仪表研究所 Sea surface evaporation waveguide near-future rapid forecasting method

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107977479A (en) * 2016-10-24 2018-05-01 中国人民解放军信息工程大学 The computational methods and device of the spatial field of electromagnetic wave
CN112711899A (en) * 2020-11-27 2021-04-27 山东省科学院海洋仪器仪表研究所 Fusion prediction method for height of evaporation waveguide

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8886498B2 (en) * 2010-12-03 2014-11-11 University Of Delaware Simulator and method for simulating an acoustic field of an acoustic waveguide
CN108663727B (en) * 2017-03-31 2020-07-03 西北工业大学 Method for estimating height of evaporation waveguide in world sea area by using evaporation rate

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107977479A (en) * 2016-10-24 2018-05-01 中国人民解放军信息工程大学 The computational methods and device of the spatial field of electromagnetic wave
CN112711899A (en) * 2020-11-27 2021-04-27 山东省科学院海洋仪器仪表研究所 Fusion prediction method for height of evaporation waveguide

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
基于NPS蒸发波导预测模型的适应性研究;杨少波等;《电子测量与仪器学报》;20161215(第12期);全文 *
环渤海大气波导监测和试验方法;王娜等;《太赫兹科学与电子信息学报》;20190630;全文 *
蒸发波导RSHMU模型预测性能分析*;张利军等;《微波学报》;20210607;第37卷(第4期);全文 *
蒸发波导模型特征及其适用性研究;郭相明等;《海洋预报》;20131015(第05期);全文 *
蒸发波导预测模型研究;宋伟等;《华中科技大学学报(自然科学版)》;20130523(第05期);全文 *

Also Published As

Publication number Publication date
CN113609639A (en) 2021-11-05

Similar Documents

Publication Publication Date Title
CN113609639B (en) Evaporation waveguide correction method suitable for stable condition
Demuth et al. Evaluation of Advanced Microwave Sounding Unit tropical-cyclone intensity and size estimation algorithms
CN112711899B (en) Fusion prediction method for height of evaporation waveguide
CN103869383B (en) Based on air data computer and its implementation of Kalman filtering algorithm
CN109543356A (en) Consider the ocean interior temperature-salinity structure remote sensing inversion method of Space atmosphere
Sjöblom et al. The turbulent kinetic energy budget in the marine atmospheric surface layer
KR20130044551A (en) A system and method for sea-surface wind -detection using satellite observation
JP6521777B2 (en) Tsunami monitoring system
Shi et al. A new evaporation duct climatology over the South China Sea
CN109145344A (en) A kind of experience ZTD model refinement method based on sounding data
CN116794621A (en) Fusion diagnosis method for electromagnetic wave propagation loss
CN113534194B (en) Troposphere temperature and humidity profile inversion method combining GNSS and wind lidar
CN108663727B (en) Method for estimating height of evaporation waveguide in world sea area by using evaporation rate
CN116559977A (en) NPS-GRYANIK20 evaporation waveguide height prediction method based on BP neural network
CN113947035A (en) Data heaven-earth correlation method for transition of hypersonic velocity boundary layer
CN116068560A (en) Marine evaporation waveguide floating type detection system and method based on radar sea clutter
CN110793649A (en) Method for correcting brightness and temperature of rough sea surface by using backscattering cross section
Barthelmie et al. Best practice for measuring wind speeds and turbulence offshore through in-situ and remote sensing technologies
CN111597692A (en) Surface net radiation estimation method, system, electronic equipment and storage medium
Zhu et al. Roles of TAO/TRITON and Argo in Tropical Pacific Observing Systems: An OSSE Study for Multiple Time Scale Variability
CN110767322B (en) Ocean floating platform hot spot stress calculation method based on response surface model
CN114252875A (en) High-precision gridding method for imaging altimeter data
Weiss-Wrana Turbulence statistics in littoral area
Liu et al. Analysis on the performances of the GNSS tropospheric delay correction models
CN113670951B (en) Microwave radiometer self-adaptive inversion algorithm based on shipborne mobile platform

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant