A highly integrated ocean atmosphere corrected refractive index profile detection device and method

By designing a highly integrated marine atmospheric correction refractive index profile detection device, using iterative calculations to similar theory to Monin Obhof, the problem of inconvenience in obtaining the refractive index profile at low altitude atmospheric correction in the existing technology is solved, and efficient and convenient detection effects are achieved.

CN119620067BActive Publication Date: 2025-04-22OCEANOGRAPHIC INSTR RES INST SHANDONG ACAD OF SCI
View PDF 2 Cites 0 Cited by

Patent Information

Application Number
CN202510161653.6
Authority / Receiving Office
CN · China
Patent Type
Patents(China)
Current Assignee / Owner
Filing Date
2025-02-14
Publication Date
2025-04-22
Estimated Expiration
2045-02-14

AI Technical Summary

Technical Problem

When obtaining the refractive index profile of low-altitude atmospheric at sea, the prior art has disadvantages such as low integration, difficult operation and maintenance, high cost, and simple calculation methods, resulting in inconvenient detection.

Method used

A highly integrated marine atmospheric correction refractive index profile detection device is designed, including main support body, louver case shell, wind speed and wind direction sensor, air intake disc, atmospheric pressure sensor, GNSS module, sea skin temperature sensor and millimeter wave radar. By iteratively calculating the stability correction function, roughness and near-formation scale parameters, the atmospheric correction refractive index profile is calculated based on the Mooningobukhov similarity theory.

Benefits of technology

Convenient handheld atmospheric corrected refractive index profile detection is realized, which improves detection efficiency, reduces operational difficulty and cost, and enhances calculation accuracy and robustness.

✦ Generated by Eureka AI based on patent content.

Smart Images

  • Figure CN119620067B_ABST
    Figure CN119620067B_ABST
Patent Text Reader

Abstract

The present invention relates to the field of marine atmosphere detection, and discloses a highly integrated marine atmosphere corrected refractive index profile detection device and method, the device comprising a main support body, a temperature and humidity sensor installed on the main support body, four wind speed and direction sensors, an atmospheric pressure sensor, a GNSS module and a GNSS antenna, a sea surface temperature sensor, an LCD display and display buttons, and a millimeter wave radar. The method comprises data acquisition; iterative calculation of stability correction function, roughness and near-surface layer scale parameters; atmospheric temperature, specific humidity and water vapor pressure profile calculation; atmospheric corrected refractive index profile calculation. The device and method disclosed by the present invention can conveniently detect the atmospheric corrected refractive index profile of the area to be measured at the current moment, thereby improving the detection efficiency, and providing fast, convenient and effective data support for aspects such as marine electromagnetic wave propagation path analysis, communication link channel research, radar target detection and remote sensing information acquisition.
Need to check novelty before this filing date? Find Prior Art

Description

Technical Field

[0001] The present invention relates to the field of ocean atmosphere detection, and in particular to a highly integrated ocean atmosphere corrected refractive index profile detection device and method. Background Art

[0002] The atmospheric correction refractive index profile determines the propagation path of electromagnetic waves in low-altitude areas at sea. Affected by the evaporation process of seawater and the inversion process near the sea surface, the atmospheric refractive index near the sea surface shows an exponential downward trend with increasing altitude, which causes the low-altitude atmospheric stratification and abnormal refraction of electromagnetic waves at sea, causing the electromagnetic waves near the sea surface to be trapped in the "waveguide layer" and produce "over-the-horizon" propagation with low loss. Because of this, a "blind spot" for electromagnetic wave propagation is formed in some areas above the waveguide layer. Therefore, it is crucial to quickly perceive the meteorological and hydrological conditions at low altitudes at sea and obtain the atmospheric correction refractive index profile to accurately predict the detection distance and detection blind spot of marine targets and realize over-the-horizon target detection and identification.

[0003] At present, the atmospheric correction refractive index profile is generally obtained by meteorological gradient observation. The meteorological gradient observation method is to set up a meteorological gradient tower near the sea surface or install a meteorological gradient detection device on a ship, use at least three layers of sensors to obtain meteorological and hydrological gradient data near the sea surface, and use a simple logarithmic function to fit and calculate the atmospheric correction refractive index profile. This kind of meteorological gradient observation requires large equipment such as offshore towers or ships, and has the disadvantages of low integration, high difficulty in operation and maintenance, high cost, and simple calculation method, which brings inconvenience to the detection of atmospheric correction refractive index profile. Therefore, it is urgent to propose a highly integrated device and method that can conveniently and effectively detect the atmospheric correction refractive index profile to solve the above problems. Summary of the invention

[0004] In order to solve the above technical problems, the present invention provides a highly integrated ocean atmosphere corrected refractive index profile detection device and method, so that the atmospheric corrected refractive index profile of the test area at the current moment can be directly detected in a convenient handheld manner, thereby improving the detection efficiency of the atmospheric corrected refractive index profile.

[0005] To achieve the above object, the technical solution of the present invention is as follows:

[0006] A highly integrated ocean atmosphere corrected refractive index profile detection device, comprising a main support body, a shutter box shell connected to the top of the main support body through a support rod 1, a temperature and humidity sensor installed inside the shutter box shell, four wind speed and direction sensors installed on the top of the shutter box shell, the top of the shutter box shell connected to an air intake disc through a support rod 2, an air intake hole is provided on the side of the air intake disc, and an atmospheric pressure sensor is installed inside; a GNSS shell is installed on the top of the air intake disc, a GNSS module is installed inside the GNSS shell, and a GNSS antenna is installed on the top;

[0007] A controller is installed in the upper part of the main support body, and a sea surface temperature sensor is installed on the outside. The sea surface temperature sensor is installed tilted downward and forms an angle of 45° with the main support body. An LCD display and display buttons are installed in the middle part of the main support body, and a processor is installed in the LCD display. A power supply is installed in the lower part of the main support body, and a device handle is installed on the outside. The bottom of the main support body is connected to a radar casing, and a millimeter-wave radar pointing downward is installed in the radar casing.

[0008] A highly integrated ocean atmosphere corrected refractive index profile detection method, using the highly integrated ocean atmosphere corrected refractive index profile detection device as described above, comprises the following steps:

[0009] Step 1: Data acquisition:

[0010] The detection device is held on the offshore platform and placed in the measured area so that the detection device is located at a surface layer 0 to 40 m above the sea surface, and the atmospheric temperature, relative humidity, atmospheric pressure, sea surface temperature, wind speed and the height of the detection device from the sea surface in the measured area are obtained;

[0011] Step 2: Iterative calculation of stability correction function, roughness and near-surface layer scale parameters:

[0012] Calculate the stability factor, and use the data in step 1 to calculate the stability correction function, roughness and near-surface layer scale parameters according to an iterative algorithm;

[0013] Step 3: Calculation of atmospheric temperature profile, specific humidity profile and water vapor pressure profile:

[0014] Based on the current stability correction function, roughness and near-surface layer scale parameters, as well as atmospheric temperature, relative humidity and atmospheric pressure, the vertical profiles of atmospheric temperature, specific humidity and water vapor pressure in the sea surface layer within a vertical height of 0 to 40 m are calculated based on the Monin Obukhov similarity theory;

[0015] Step 4: Calculation of atmospheric correction refractive index profile:

[0016] Based on the calculated vertical profiles of atmospheric temperature, specific humidity and water vapor pressure, the atmospheric refractive index formula is used to calculate the atmospheric refractive index profile under the current environment, and it is corrected for the curvature of the earth to obtain the atmospheric corrected refractive index profile in the sea surface layer.

[0017] In the above scheme, the specific method of step 2 is as follows:

[0018] S1: Data input: input the atmospheric temperature AT, relative humidity RH, atmospheric pressure AP, sea surface temperature SST, wind speed WS and the height H of the detection device from the sea surface collected in step 1;

[0019] S2: Calculation of stability correction function: Calculate the stability factor and calculate the stability correction function based on the stability factor; the stability correction function includes the wind speed stability correction function , Temperature stability correction function and humidity stability correction function ;

[0020] S3: Roughness calculation: The roughness includes temperature roughness , humidity roughness and dynamic roughness ;

[0021] S4: Calculation of near-surface scale parameters: Near-surface scale parameters include friction velocity , characteristic potential temperature and characteristic specific humidity Three types;

[0022] S5: Iteration: Recalculate the stability factor , the stability of the current atmospheric environment is judged according to the stability factor. When the atmospheric environment is in a strong stable state, the calculation is terminated and the current stability correction function, roughness and near-surface scale parameter values ​​are obtained. Otherwise, S2, S3 and S4 are iteratively calculated to obtain the stability correction function, roughness and near-surface scale parameter values ​​in a non-strong stable state. When the number of iterations is When , the iteration ends and the current stability correction function, roughness and near-surface layer scale parameter values ​​are output.

[0023] In a further technical solution, the stability factor calculation method is as follows:

[0024] ;

[0025] in, It is any height between 0 and 40 m in the sea surface layer, in m. is the Monin Obukhov length, which is calculated as follows:

[0026] ;

[0027] in, is the friction speed, is the acceleration due to gravity, is the Karman constant, is the temperature scale parameter.

[0028] In a further technical solution, for the wind speed stability correction function , Temperature stability correction function and humidity stability correction function , select one of the BH91, Businger-Dyer, CB05, SHEBA07 and HYQ92 schemes under neutral and stable atmospheric conditions, and select one of the Grachev00, Fairall96, Fairall03, Edson04 and HYQ92 schemes under unstable atmospheric conditions.

[0029] In a further technical solution, the dynamic roughness Adopt one of the schemes YT96, TY01, O02, GW06 and PS07, the temperature roughness and humidity roughness The same calculation method is as follows:

[0030] ;

[0031] Among them, the Reynolds number , is the sea surface viscosity coefficient, is the friction speed.

[0032] In a further technical solution, the friction speed , characteristic potential temperature and characteristic specific humidity The calculation method is as follows:

[0033] ;

[0034] ;

[0035] ;

[0036] in, is the Karman constant, , and are the differences in wind speed, temperature and specific humidity at height z and at the roughness point, respectively.

[0037] In the above scheme, in step 3, the atmospheric temperature profile is calculated based on the Monin-Obukhov similarity theory Humidity Profile , the calculation method is as follows:

[0038] ;

[0039] ;

[0040] in, Any height between 0 and 40m, is the sea surface temperature, K / m, is the specific humidity at the sea surface; is the Karman constant, is the characteristic potential temperature, is the temperature roughness, is the temperature stability correction function, is the characteristic specific humidity, is the humidity roughness, is the humidity stability correction function;

[0041] The atmospheric pressure profile is obtained by integrating the fluid statics equation and the ideal gas law. , the calculation method is as follows:

[0042] ;

[0043] Where R is the gas constant, Any height z between 0 and 40m and the height of the device The average value of the virtual temperature; The height of the device The atmospheric pressure at

[0044] Through the atmospheric pressure profile Calculating the Vapor Pressure Profile , the calculation method is as follows:

[0045] ;

[0046] in, It is the ratio of the gas constant of dry air to the gas constant of water vapor.

[0047] In the above scheme, in step 4, the atmospheric temperature profile, specific humidity profile and water vapor pressure profile are used to calculate the atmospheric refractive index profile within the vertical height of 0 to 40m in the sea surface layer. , the calculation method is as follows:

[0048] ;

[0049] in, is the atmospheric temperature profile, is the specific humidity profile, is the water vapor pressure profile.

[0050] In the above scheme, in step 4, the atmospheric refractive index profile is The earth curvature correction is performed to obtain the atmospheric corrected refractive index profile within the vertical height of 0 to 40m in the sea surface layer. , the correction method is as follows:

[0051] ;

[0052] The radius of curvature of the earth km, Any height between 0 and 40m.

[0053] Through the above technical solution, the highly integrated ocean atmosphere corrected refractive index profile detection device and method provided by the present invention have the following beneficial effects:

[0054] The present invention reveals the logarithmic rate characteristics of the vertical change of near-surface wind speed based on the Monin-Obukhov similarity theory. This theory provides an important basis for the design of the calculation framework of the present invention. Through iterative calculation of the stability correction function, roughness and near-surface layer scale parameters, the present invention uses a highly integrated marine atmosphere corrected refractive index profile detection device to replace the complex and large offshore gradient system, greatly reducing the operating difficulty of the system, significantly improving the ease of use, flexibility and maintainability of the system, and providing a more concise and effective way to obtain the near-sea surface atmospheric corrected refractive index profile.

[0055] The present invention can effectively reduce the deviations or errors that may exist in a single solution through the flexible combination of stability correction functions and roughness, thereby improving the accuracy of the overall calculation, enhancing the robustness, and improving the generalization ability, so that the device can achieve accurate detection of the atmospheric corrected refractive index profile in multiple sea areas. BRIEF DESCRIPTION OF THE DRAWINGS

[0056] In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings required for use in the embodiments or the description of the prior art are briefly introduced below.

[0057] Figure 1 A schematic diagram of a highly integrated ocean atmosphere corrected refractive index profile detection device disclosed in an embodiment of the present invention;

[0058] Figure 2 A schematic flow chart of a highly integrated ocean atmosphere corrected refractive index profile detection method disclosed in an embodiment of the present invention;

[0059] Figure 3 Flowchart for iterative calculation of stability correction function, roughness and near-surface scale parameters;

[0060] Figure 4 This is an example diagram of the detection results of the atmospheric correction refractive index profile;

[0061] Figure 5 This is a schematic diagram of the LCD display interface.

[0062] In the figure, 1. main support body; 2. support rod 1; 3. louver box shell; 4. wind speed and direction sensor; 5. support rod 2; 6. air intake disc; 7. air intake hole; 8. GNSS shell; 9. GNSS antenna; 10. sea surface temperature sensor; 11. LCD display; 12. display button; 13. device handle; 14. radar shell. DETAILED DESCRIPTION

[0063] The technical solutions in the embodiments of the present invention will be described clearly and completely below in conjunction with the accompanying drawings in the embodiments of the present invention.

[0064] The present invention provides a highly integrated ocean atmosphere corrected refractive index profile detection device, such as Figure 1 As shown, it includes a main support body 1, a louvered box housing 3 made of plastic material connected to the main support body 1 through a support rod 2, and a temperature and humidity sensor is installed inside the louvered box housing 3 to detect the temperature AT (K) and relative humidity RH (%) in the atmospheric environment where the device is located when in use. The louvered box housing is used to prevent direct sunlight and rain from sticking to the temperature and humidity sensor probe, and the temperature and humidity sensor probe uses the HMP155 model.

[0065] Four wind speed and direction sensors 4 are installed on the top of the shutter box shell to detect the wind speed WS (m / s) at the location where the device is in use. The four wind speed and direction sensors use ultrasonic generators and are distributed at 90° on the top of the shutter box shell.

[0066] The top of the shutter box shell is connected to the air intake disc 6 through the support rod 2 5. The air intake disc 6 is made of plastic. An air intake hole 7 is opened on the side of the air intake disc 6. An atmospheric pressure sensor is installed inside. The atmospheric pressure sensor is used to detect the atmospheric pressure AP (hPa) at the location where the device is used.

[0067] A GNSS housing 8 is installed on the top of the air intake disc 6, a GNSS module is installed inside the GNSS housing, and a GNSS antenna 9 is installed on the top; in this embodiment, a Beidou satellite navigation system is adopted, including a Beidou antenna and a Beidou chip, and an ATK1218-BD Beidou chip is used based on the NMEA-0183 protocol to obtain the longitude, latitude and real-time time of the location of the device when in use.

[0068] A controller is installed on the upper part of the main support body 1, and a sea surface temperature sensor 10 is installed on the outside. The sea surface temperature sensor uses a BRW600-400A infrared temperature measuring probe. The sea surface temperature sensor is installed tilted downward at a 45° angle to the main support body 1 and is used to detect the sea surface temperature SST (K) below the position where the system is used.

[0069] An LCD display 11 and display buttons 12 are installed in the middle of the main support body 1, and a processor is installed in the LCD display.

[0070] A power supply is installed at the lower part of the main support body 1, with a voltage of 12V. The current varies with the load size of each unit when it is running. The main support body 1 is located outside the power supply as a device handle 13. The bottom of the main support body 1 is connected to a radar housing 14. A millimeter-wave radar pointing downward is installed in the radar housing 14 to detect the height H (m) from the sea surface when the system is in use. The sensors are connected to the processor using RS485 communication.

[0071] A highly integrated ocean atmosphere corrected refractive index profile detection method, using the detection device as described above, such as Figure 2 As shown, the following steps are included:

[0072] Step 1: Data acquisition:

[0073] The detection device is held on the offshore platform and placed in the measured area so that the detection device is located in the surface layer 0 to 40 m above the sea surface, and the atmospheric temperature AT (K), relative humidity RH (%), atmospheric pressure AP (hPa), sea surface temperature SST (K), wind speed WS (m / s) and the height H (m) of the detection device from the sea surface are obtained.

[0074] The device should be held above the head to prevent the body from blocking the wind speed measurement. The millimeter wave radar under the device should be pointed vertically toward the sea level, and the sea surface temperature sensor in the device should be pointed toward the sea surface.

[0075] Offshore platforms include but are not limited to offshore meteorological towers, various types of ships and vessels, offshore buoy platforms, and other offshore platforms that can meet the requirements of handheld detection.

[0076] Step 2: Iterative calculation of stability correction function, roughness and near-surface layer scale parameters:

[0077] Calculate the stability factor, and use the data in step 1 to calculate the stability correction function, roughness and near-surface layer scale parameters according to the iterative algorithm.

[0078] like Figure 3 As shown, the specific method is as follows:

[0079] S1: Data input: input the atmospheric temperature AT, relative humidity RH, atmospheric pressure AP, sea surface temperature SST, wind speed WS and the height H of the detection device from the sea surface collected in step 1;

[0080] S2: Stability correction function calculation: Calculate the stability factor and calculate the stability correction function based on the stability factor; the stability correction function includes the wind speed stability correction function , Temperature stability correction function and humidity stability correction function ; For the three stability correction functions of wind speed, temperature and humidity, one of the schemes including but not limited to BH91, Businger-Dyer, CB05, SHEBA07 and HYQ92 is selected under neutral and stable atmospheric conditions, and one of the schemes including but not limited to Grachev00, Fairall96, Fairall03, Edson04 and HYQ92 is selected under unstable atmospheric conditions. The embodiment of the present invention uses the BH91 scheme under neutral and stable atmospheric conditions, and uses Grachev00 under unstable atmospheric conditions for example:

[0081] Stability Factor The calculation method is as follows:

[0082] ;

[0083] in, It is any height between 0 and 40 m in the sea surface layer, in m. is the Monin Obukhov length, which is calculated as follows:

[0084] ;

[0085] in, is the friction speed, is the acceleration due to gravity, is the Karman constant, is the temperature scale parameter.

[0086] BH91 wind speed stability correction function The calculation method is as follows:

[0087] ;

[0088] The temperature stability correction function of the BH91 scheme is the same as the humidity stability correction function, and the calculation method is as follows:

[0089] .

[0090] Wind speed stability correction function of Grachev00 scheme The calculation method is as follows:

[0091] ;

[0092] in:

[0093] ;

[0094] ;

[0095] ;

[0096] .

[0097] The temperature stability correction function of the Grachev00 scheme is the same as the humidity stability correction function, and the calculation method is as follows:

[0098] ;

[0099] in:

[0100] ;

[0101] ;

[0102] ;

[0103] .

[0104] S3: Roughness calculation: Roughness includes temperature roughness , humidity roughness and dynamic roughness , where the dynamic roughness Including but not limited to one of the solutions YT96, TY01, O02, GW06 and PS07, temperature roughness and humidity roughness The same calculation method is as follows:

[0105] ;

[0106] Where: Reynolds number , is the sea surface viscosity coefficient, is the friction speed.

[0107] The dynamic roughness of the present invention Using the YT96 scheme as an example:

[0108] ;

[0109] in:

[0110] .

[0111] S4: Calculation of near-surface scale parameters: Near-surface scale parameters include friction velocity , characteristic potential temperature and characteristic specific humidity Three types; the calculation method is as follows:

[0112] ;

[0113] ;

[0114] ;

[0115] in, is the Karman constant, , and are the differences in wind speed, temperature and specific humidity at height z and at the roughness point, respectively.

[0116] S5: Iteration: Recalculate the stability factor , the stability of the current atmospheric environment is judged based on the stability factor. When the atmospheric environment is in a strong stable state ( ), end the calculation and obtain the current stability correction function, roughness and near-surface scale parameter values, otherwise iterate S2, S3 and S4 and obtain the stability correction function, roughness and near-surface scale parameter values ​​in the non-strong stable state. When the number of iterations is When , the iteration ends and the current stability correction function, roughness and near-surface layer scale parameter values ​​are output.

[0117] Step 3: Calculation of atmospheric temperature profile, specific humidity profile and water vapor pressure profile:

[0118] According to the current stability correction function, roughness and near-surface layer scale parameters as well as atmospheric temperature, relative humidity and atmospheric pressure, the vertical profiles of atmospheric temperature, specific humidity and water vapor pressure within the vertical height of 0 to 40 m in the sea surface layer are calculated based on the Monin Obukhov similarity theory.

[0119] Calculation of atmospheric temperature profile based on Monin-Obukhov similarity theory Humidity Profile , the calculation method is as follows:

[0120] ;

[0121] ;

[0122] in, Any height between 0 and 40m, is the sea surface temperature, K / m, is the specific humidity at the sea surface; is the Karman constant, is the characteristic potential temperature, is the temperature roughness, is the temperature stability correction function, is the characteristic specific humidity, is the humidity roughness, is the humidity stability correction function.

[0123] The atmospheric pressure profile is obtained by integrating the fluid statics equation and the ideal gas law. , the calculation method is as follows:

[0124] ;

[0125] Where R is the gas constant, Any height z between 0 and 40m and the height of the device The average value of the virtual temperature; The height of the device The atmospheric pressure at.

[0126] Through the atmospheric pressure profile Calculating the Vapor Pressure Profile , the calculation method is as follows:

[0127] ;

[0128] in, is the ratio of the gas constant of dry air to the gas constant of water vapor, .

[0129] Step 4: Calculation of atmospheric correction refractive index profile:

[0130] Based on the calculated vertical profiles of atmospheric temperature, specific humidity and water vapor pressure, the atmospheric refractive index formula is used to calculate the atmospheric refractive index profile under the current environment, and it is corrected for the curvature of the earth to obtain the atmospheric corrected refractive index profile in the sea surface layer.

[0131] The atmospheric refractivity profile within the vertical height of 0 to 40 m in the sea surface layer is calculated using the atmospheric temperature profile, specific humidity profile and water vapor pressure profile. , the calculation method is as follows:

[0132] ;

[0133] in, is the atmospheric temperature profile, is the specific humidity profile, is the water vapor pressure profile.

[0134] Atmospheric refractive index profile The earth curvature correction is performed to obtain the atmospheric corrected refractive index profile within the vertical height of 0 to 40m in the sea surface layer. , the correction method is as follows:

[0135] ;

[0136] The radius of curvature of the earth km, The atmospheric correction refractive index profile is obtained for any altitude between 0 and 40 m. Schematic diagram as Figure 4 shown.

[0137] The longitude, latitude, real time, data collected in step one, and atmospheric correction refractive index profile information in the vertical height of 0 to 40 m in the sea surface layer obtained by precise calculation in step four are completely saved to the SD card in the LCD display through the SDIO protocol.

[0138] like Figure 5 As shown, the longitude, latitude, real-time time, data collected in step one, and the atmospheric correction refractive index profile information in the vertical height of 0 to 40m in the sea surface layer obtained by precise calculation in step four are displayed on the LCD display in an intuitive and clear manner through the HDMI protocol to facilitate further analysis and viewing. The upper part of the LCD display is used to display the current ambient atmospheric temperature AT (K), relative humidity RH (%), atmospheric pressure AP (hPa), sea surface temperature SST (K), wind speed WS (m / s) and the height H (m) of the detection device from the sea surface. The middle part of the LCD display shows the state of the atmospheric correction refractive index profile, and the lower part of the LCD display shows the longitude, latitude and real-time time of the current position.

[0139] The above description of the disclosed embodiments enables one skilled in the art to implement or use the present invention. Various modifications to these embodiments will be apparent to one skilled in the art, and the general principles defined herein may be implemented in other embodiments without departing from the spirit or scope of the present invention. Therefore, the present invention will not be limited to the embodiments shown herein, but rather to the widest scope consistent with the principles and novel features disclosed herein.

Claims

1. A highly integrated ocean atmosphere corrected refractive index profile detection method, characterized in that: The steps include: Step 1: Data acquisition: A handheld detection device is placed on the offshore platform in the measured area so that the detection device is located at a surface layer 0 to 40 m above the sea surface, and the atmospheric temperature, relative humidity, atmospheric pressure, sea surface temperature, wind speed and the height of the detection device from the sea surface in the measured area are obtained; Step 2: Iterative calculation of stability correction function, roughness and near-surface layer scale parameters: Calculate the stability factor, and use the data in step 1 to calculate the stability correction function, roughness and near-surface layer scale parameters according to an iterative algorithm; Step 3: Calculation of atmospheric temperature profile, specific humidity profile and water vapor pressure profile: Based on the current stability correction function, roughness and near-surface layer scale parameters, as well as atmospheric temperature, relative humidity and atmospheric pressure, the vertical profiles of atmospheric temperature, specific humidity and water vapor pressure in the sea surface layer within a vertical height of 0 to 40 m are calculated based on the Monin Obukhov similarity theory; Step 4: Calculation of atmospheric correction refractive index profile: Based on the calculated vertical profiles of atmospheric temperature, specific humidity and water vapor pressure, the atmospheric refractive index formula is used to calculate the atmospheric refractive index profile under the current environment, and the atmospheric refractive index profile is corrected for the curvature of the earth to obtain the atmospheric corrected refractive index profile in the sea surface layer; In step 3, the atmospheric temperature profile is calculated based on the Monin-Obukhov similarity theory. Humidity Profile , the calculation method is as follows: ; ; in, Any height between 0 and 40m, is the sea surface temperature, K / m, is the specific humidity at the sea surface; is the Karman constant, is the characteristic potential temperature, is the temperature roughness, is the temperature stability correction function, is the characteristic specific humidity, is the humidity roughness, is the humidity stability correction function; The atmospheric pressure profile is obtained by integrating the fluid statics equation and the ideal gas law. , the calculation method is as follows: ; Where R is the gas constant, Any height z between 0 and 40m and the height of the device The average value of the virtual temperature; The height of the device The atmospheric pressure at is the acceleration due to gravity; Through the atmospheric pressure profile Calculating the Vapor Pressure Profile , the calculation method is as follows: ; in, It is the ratio of the gas constant of dry air to the gas constant of water vapor.

2. A highly integrated ocean atmosphere corrected refractive index profile detection method according to claim 1, characterized in that: The specific method of step 2 is as follows: S1: Data input: input the atmospheric temperature AT, relative humidity RH, atmospheric pressure AP, sea surface temperature SST, wind speed WS and the height H of the detection device from the sea surface collected in step 1; S2: Calculation of stability correction function: Calculate the stability factor and calculate the stability correction function based on the stability factor; the stability correction function includes the wind speed stability correction function , Temperature stability correction function and humidity stability correction function ; S3: Roughness calculation: The roughness includes temperature roughness , humidity roughness and dynamic roughness ; S4: Calculation of near-surface scale parameters: Near-surface scale parameters include friction velocity , characteristic potential temperature and characteristic specific humidity Three types; S5: Iteration: Recalculate the stability factor , the stability of the current atmospheric environment is judged according to the stability factor. When the atmospheric environment is in a strong stable state, the calculation is terminated and the current stability correction function, roughness and near-surface scale parameter values ​​are obtained. Otherwise, S2, S3 and S4 are iteratively calculated to obtain the stability correction function, roughness and near-surface scale parameter values ​​in a non-strong stable state. When the number of iterations is When , the iteration ends and the current stability correction function, roughness and near-surface layer scale parameter values ​​are output.

3. A highly integrated ocean atmosphere corrected refractive index profile detection method according to claim 2, characterized in that: The stability factor calculation method is as follows: ; in, It is any height between 0 and 40 m in the sea surface layer, in m. is the Monin Obukhov length, which is calculated as follows: ; in, is the friction speed, is the acceleration due to gravity, is the Karman constant, is the temperature scale parameter.

4. The highly integrated ocean atmosphere corrected refractive index profile detection method according to claim 2, characterized in that: Correction function for wind speed stability , Temperature stability correction function and humidity stability correction function , select one of the BH91, Businger-Dyer, CB05, SHEBA07 and HYQ92 schemes under neutral and stable atmospheric conditions, and select one of the Grachev00, Fairall96, Fairall03, Edson04 and HYQ92 schemes under unstable atmospheric conditions.

5. The highly integrated ocean atmosphere corrected refractive index profile detection method according to claim 3, characterized in that: The dynamic roughness Adopt one of the schemes YT96, TY01, O02, GW06 and PS07, the temperature roughness and humidity roughness The same calculation method is as follows: ; Among them, the Reynolds number , is the sea surface viscosity coefficient, is the friction speed.

6. The highly integrated ocean atmosphere corrected refractive index profile detection method according to claim 2, characterized in that: The friction speed , characteristic potential temperature and characteristic specific humidity The calculation method is as follows: ; ; ; in, is the Karman constant, , and are the differences in wind speed, temperature and specific humidity at height z and at the roughness point, respectively.

7. The highly integrated ocean atmosphere corrected refractive index profile detection method according to claim 1, characterized in that: In step 4, the atmospheric temperature profile, specific humidity profile and water vapor pressure profile are used to calculate the atmospheric refractive index profile within the vertical height of 0 to 40m in the sea surface layer. , the calculation method is as follows: ; in, is the atmospheric temperature profile, is the specific humidity profile, is the water vapor pressure profile.

8. The highly integrated ocean atmosphere corrected refractive index profile detection method according to claim 1, characterized in that: In step 4, the atmospheric refractive index profile The earth curvature correction is performed to obtain the atmospheric corrected refractive index profile within the vertical height of 0 to 40m in the sea surface layer. , the correction method is as follows: ; The radius of curvature of the earth km, Any height between 0 and 40m.

9. A highly integrated ocean atmosphere corrected refractive index profile detection device using the method as claimed in claim 1, characterized in that: It includes a main support body, a shutter box shell connected to the main support body through a support rod 1, a temperature and humidity sensor is installed inside the shutter box shell, four wind speed and direction sensors are installed on the top of the shutter box shell, the top of the shutter box shell is connected to an air intake disc through a support rod 2, an air intake hole is provided on the side of the air intake disc, and an atmospheric pressure sensor is installed inside; a GNSS shell is installed on the top of the air intake disc, a GNSS module is installed inside the GNSS shell, and a GNSS antenna is installed on the top; A controller is installed in the upper part of the main support body, and a sea surface temperature sensor is installed on the outside. The sea surface temperature sensor is installed tilted downward and forms an angle of 45° with the main support body. An LCD display and display buttons are installed in the middle part of the main support body, and a processor is installed in the LCD display. A power supply is installed in the lower part of the main support body, and a device handle is installed on the outside. The bottom of the main support body is connected to a radar casing, and a millimeter-wave radar pointing downward is installed in the radar casing.

Citation Information

Patent Citations

  • Evaporation waveguide correction model suitable for stable conditions

    CN113609639A

  • Optimization method for evaporation waveguide prediction model parameters

    CN116822567A