A method, device, equipment, medium and product for obtaining a vertical distribution profile of atmospheric heating rate by aerosol

Abstract

The application discloses a method and device for obtaining an aerosol-atmosphere heating rate vertical distribution profile, equipment, medium and product, and relates to the technical field of atmospheric environment monitoring. The method comprises the following steps: for each height corresponding aerosol parameters and meteorological parameters, the following processing is performed: fitting aerosol absorption coefficients of several wavelengths to obtain continuous wavelength aerosol absorption coefficients; calculating an aerosol absorption rate based on actinic flux and the continuous wavelength aerosol absorption coefficients; calculating air density and constant-pressure specific heat capacity based on the meteorological parameters; calculating the aerosol-atmosphere heating rate according to the aerosol absorption rate, the air density and the constant-pressure specific heat capacity; and determining the vertical distribution profile according to the aerosol-atmosphere heating rate calculated at each height. The application realizes high-precision and high-resolution vertical heating rate vertical distribution profile acquisition by simultaneously measuring the actinic flux, the aerosol absorption coefficient and the meteorological parameters at each height.

Description

Technical Field

[0001] This application relates to the field of atmospheric environmental monitoring technology, and in particular to a method, apparatus, equipment, medium and product for obtaining the vertical distribution profile of the atmospheric heating rate of aerosols. Background Technology

[0002] Aerosols are a crucial component of Earth's atmosphere, profoundly impacting climate, air quality, ecosystems, and human health. Aerosols influence the atmospheric energy balance by absorbing and scattering solar radiation. Absorbent aerosols (such as black carbon, brown carbon, and dust) heat the atmosphere after absorbing solar radiation, and their heating rate is a key parameter for quantifying the aerosol radiation effect. Traditional methods for obtaining the vertical distribution profile of aerosol heating rates include ground-based observations and satellite remote sensing. However, traditional ground-based observations use atmospheric radiative transfer models to obtain photochemical flux, and many parameters are based on assumptions rather than actual measurements, resulting in significant uncertainty in the obtained vertical distribution profile of the heating rate. Satellite remote sensing inversion methods have low accuracy in the vertical direction, making it difficult to meet the needs for refined acquisition of heating rate profiles. Summary of the Invention

[0003] The purpose of this application is to provide a method, apparatus, equipment, medium, and product for obtaining the vertical distribution profile of the heating rate of aerosols to the atmosphere, which can achieve high-precision and high-resolution acquisition of the vertical distribution profile of the heating rate.

[0004] To achieve the above objectives, this application provides the following solution:

[0005] In a first aspect, this application provides a method for obtaining the vertical distribution profile of the atmospheric heating rate of aerosols, comprising:

[0006] Acquire a number of aerosol parameters and meteorological parameters corresponding to several altitudes; the aerosol parameters include: photochemical flux and aerosol absorption coefficients at several wavelengths;

[0007] For each altitude, the corresponding aerosol and meteorological parameters are processed as follows:

[0008] By fitting the aerosol absorption coefficients at several wavelengths, the aerosol absorption coefficients at continuous wavelengths are obtained.

[0009] The aerosol absorption rate is calculated based on the photochemical flux and the aerosol absorption coefficient of the continuous wavelength.

[0010] Based on the meteorological parameters, calculate the air density and specific heat capacity at constant pressure;

[0011] The atmospheric heating rate of the aerosol is calculated based on the aerosol's absorption rate, the air density, and the specific heat capacity at constant pressure.

[0012] Based on the calculated aerosol heating rate to the atmosphere at each altitude, the vertical distribution profile of the aerosol heating rate to the atmosphere as a function of altitude is determined.

[0013] Optionally, the fitting of aerosol absorption coefficients at several wavelengths to obtain aerosol absorption coefficients at continuous wavelengths specifically includes:

[0014] By fitting the aerosol absorption coefficient to the wavelength at several wavelengths, a functional relationship between the absorption coefficient and the wavelength is obtained; the functional relationship between the absorption coefficient and the wavelength is as follows: Where C is a constant value obtained after fitting, and AAE is the wavelength absorption index obtained after fitting. b abs ( λ () represents the aerosol absorption coefficient at wavelength λ;

[0015] Based on the functional relationship between the absorption coefficient and wavelength, the aerosol absorption coefficient for continuous wavelengths is calculated.

[0016] Optionally, based on the photochemical flux and the aerosol absorption coefficient of the continuous wavelength, the formula for calculating the aerosol absorption rate is as follows:

[0017] ;

[0018] Among them: A areosol Indicates the aerosol absorption rate. This represents the aerosol absorption coefficient at wavelength λ, where λ represents the wavelength. This represents the photochemical flux at wavelength λ.

[0019] Optionally, the formula for calculating the air density is:

[0020] ;

[0021] in, ρ represents the density of moist air; P represents the atmospheric pressure of moist air; Pv represents the water vapor pressure; Rd represents the dry air specific gas constant; T represents the temperature; Rv represents the moist air specific gas constant; RH represents the relative humidity; Pvs represents the saturated water vapor partial pressure corresponding to the measured temperature.

[0022] The formula for calculating the isobaric specific heat capacity is:

[0023] ;

[0024] in, ; q represents specific humidity; Pd represents dry air pressure; C P C represents the specific heat capacity of moist air at constant pressure. Pd C represents the specific heat capacity of dry air at constant pressure. PVThis represents the specific heat capacity of water vapor at constant pressure.

[0025] Optionally, based on the aerosol's absorption rate, the air density, and the isobaric specific heat capacity, the formula for calculating the aerosol's heating rate to the atmosphere is as follows:

[0026] ;

[0027] Wherein: H areosol A represents the rate at which aerosols heat the atmosphere. areosol The aerosol absorption rate is represented by ρ, and the air density is represented by C. P This indicates the specific heat capacity at constant pressure.

[0028] Optionally, calculating the aerosol absorption rate based on the photochemical flux and the aerosol absorption coefficient of the continuous wavelength specifically includes:

[0029] Based on the solar elevation angle, solar azimuth angle and UAV attitude three-dimensional tilt angle when the photochemical flux is obtained, the cosine correction coefficient is determined;

[0030] Based on the cosine correction coefficient, the photochemical flux is cosine corrected to obtain the corrected photochemical flux, calculated using the following formula: ;in, This indicates the correction for photochemical flux. Represents photochemical flux. This represents the cosine correction factor;

[0031] The aerosol absorption rate is calculated based on the corrected photochemical flux and the aerosol absorption coefficient of the continuous wavelength.

[0032] Secondly, this application provides an apparatus for obtaining the vertical distribution profile of the atmospheric heating rate of aerosols, including: an unmanned aerial vehicle platform, a photochemical flux measuring instrument, a multi-wavelength aerosol absorption coefficient measuring instrument, a meteorological sensor, and a processor;

[0033] The drone platform supports rotary-wing or fixed-wing aircraft.

[0034] The photochemical flux measuring instrument is equipped with a hemispherical omnidirectional probe;

[0035] The photochemical flux measuring instrument can measure the following wavelength ranges: ultraviolet, visible light, and infrared.

[0036] The multi-wavelength aerosol absorption coefficient measuring instrument is used to measure at least five wavelengths in the ultraviolet, visible, and infrared bands.

[0037] The meteorological sensor is used to measure relative humidity, temperature, and atmospheric pressure.

[0038] The processor is used for:

[0039] Acquire a number of aerosol parameters and meteorological parameters corresponding to several altitudes; the aerosol parameters include: photochemical flux and aerosol absorption coefficients at several wavelengths;

[0040] For each altitude, the corresponding aerosol and meteorological parameters are processed as follows:

[0041] By fitting the aerosol absorption coefficients at several wavelengths, the aerosol absorption coefficients at continuous wavelengths are obtained.

[0042] The aerosol absorption rate is calculated based on the photochemical flux and the aerosol absorption coefficient of the continuous wavelength.

[0043] Based on the meteorological parameters, calculate the air density and specific heat capacity at constant pressure;

[0044] The atmospheric heating rate of the aerosol is calculated based on the aerosol's absorption rate, the air density, and the specific heat capacity at constant pressure.

[0045] Based on the calculated aerosol heating rate to the atmosphere at each altitude, the vertical distribution profile of the aerosol heating rate to the atmosphere as a function of altitude is determined.

[0046] Thirdly, this application provides a computer device, including: a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the computer program to implement the method described above for obtaining the vertical distribution profile of the aerosol heating rate to the atmosphere.

[0047] Fourthly, this application provides a computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the method described above for obtaining the vertical distribution profile of the aerosol heating rate to the atmosphere.

[0048] Fifthly, this application provides a computer program product, including a computer program that, when executed by a processor, implements the method described above for obtaining the vertical distribution profile of the aerosol heating rate to the atmosphere.

[0049] According to the specific embodiments provided in this application, the following technical effects are disclosed:

[0050] This application provides a method, apparatus, equipment, medium, and product for obtaining the vertical distribution profile of the aerosol heating rate of the atmosphere. The method includes: acquiring aerosol parameters and meteorological parameters corresponding to several altitudes; the aerosol parameters include: photochemical flux and aerosol absorption coefficients at several wavelengths; for the aerosol parameters and meteorological parameters corresponding to each altitude, the following processing is performed: fitting the aerosol absorption coefficients at several wavelengths to obtain aerosol absorption coefficients at continuous wavelengths; calculating the aerosol absorption rate based on the photochemical flux and the aerosol absorption coefficients at continuous wavelengths; calculating the air density and isobaric specific heat capacity based on the meteorological parameters; calculating the aerosol heating rate of the atmosphere based on the aerosol absorption rate, the air density, and the isobaric specific heat capacity; and determining the vertical distribution profile of the aerosol heating rate of the atmosphere as a function of altitude based on the aerosol heating rate calculated at each altitude. This application achieves high-precision, high-resolution acquisition of the vertical distribution profile of aerosol heating rate by simultaneously measuring the photochemical flux, aerosol absorption coefficient, and meteorological parameters at each altitude, and by fitting continuous wavelengths. Attached Figure Description

[0051] To more clearly illustrate the technical solutions in the embodiments of this application or the prior art, the drawings used in the embodiments will be briefly introduced below. Obviously, the drawings described below are only some embodiments of this application. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0052] Figure 1 This is an application environment diagram of a method for obtaining the vertical distribution profile of the atmospheric heating rate of aerosols according to an embodiment of this application.

[0053] Figure 2 This is a flowchart illustrating a method for obtaining the vertical distribution profile of the atmospheric heating rate of aerosols, as provided in an embodiment of this application.

[0054] Figure 3 This is a schematic diagram of the structure of a computer device provided in an embodiment of this application. Detailed Implementation

[0055] The technical solutions of the embodiments of this application will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of this application, and not all embodiments. Based on the embodiments of this application, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of this application.

[0056] This application proposes a method for obtaining the vertical distribution profile of aerosol heating rate based on measured photochemical flux, absorption coefficient, and meteorological elements from a UAV. The method includes: acquiring the latitude and longitude coordinates of the location to be tested, the solar altitude angle, solar azimuth angle, photochemical flux, and three-dimensional attitude data of the UAV corresponding to the observation time; acquiring aerosol absorption coefficients at multiple wavelengths; determining the aerosol absorption coefficients for continuous wavelengths within the measurement range; determining the aerosol absorption rate within the measurement band; acquiring meteorological parameters within the detection altitude range, including temperature, relative humidity, and atmospheric pressure; and calculating the aerosol heating rate corresponding to the measurement altitude based on the above data, thereby determining the profile of the aerosol heating rate as a function of altitude.

[0057] This application relates to the field of atmospheric environment monitoring technology, and in particular to a method for calculating the vertical profile of aerosol heating rate based on measured photochemical flux, absorption coefficient and meteorological elements (atmospheric temperature, relative humidity and atmospheric pressure) data from an unmanned aerial vehicle platform. This method is applicable to atmospheric pollution research, climate effect assessment and environmental governance.

[0058] To make the above-mentioned objectives, features and advantages of this application more apparent and understandable, the application will be further described in detail below with reference to the accompanying drawings and specific embodiments.

[0059] The method for obtaining the vertical distribution profile of the aerosol heating rate to the atmosphere provided in this application embodiment can be applied to, for example... Figure 1 In the application environment shown, terminal 102 communicates with server 104 via a network. A data storage system can store the data that server 104 needs to process. The data storage system can be set up independently, integrated into server 104, or placed in the cloud or on another server.

[0060] The terminal 102 can be, but is not limited to, various desktop computers, laptops, smartphones, tablets, IoT devices, and portable wearable devices. IoT devices can include smart speakers, smart TVs, smart air conditioners, and smart in-vehicle devices. Portable wearable devices can include smartwatches, smart bracelets, and head-mounted devices. The server 104 can be implemented using a standalone server or a server cluster composed of multiple servers, or it can be a cloud server.

[0061] In one exemplary embodiment, such as Figure 2 As shown, a method for obtaining the vertical distribution profile of the atmospheric heating rate of aerosols is provided. This method is executed by a computer device, specifically by a terminal or server alone, or by both a terminal and a server. In this embodiment, the method is applied to... Figure 1 Taking server 104 as an example, the explanation includes the following steps S1 to S6. Wherein:

[0062] S1. Obtain several aerosol parameters and meteorological parameters corresponding to several altitudes; the aerosol parameters include: photochemical flux and aerosol absorption coefficients at several wavelengths.

[0063] In this embodiment, the photochemical flux corresponding to the detection location and detection time is obtained. w / (m 2 •nm, watts per (m²·nanometer). Obtaining aerosol absorption coefficients (Mm) at multiple wavelengths. -1 (megameters).

[0064] The wavelength measurement range for photochemical flux and aerosol absorption coefficient must cover the ultraviolet, visible, and infrared bands, with photochemical flux observations being continuous wavelengths and absorption coefficient observations including at least five wavelengths. The number of wavelengths depends on the number of wavelengths that the aerosol absorption coefficient measuring instrument onboard the UAV can detect.

[0065] The meteorological parameters include temperature (T, K), relative humidity (RH, %), and atmospheric pressure (P, Pa). The accuracy of the meteorological parameters depends on the meteorological sensors they are equipped with.

[0066] The photochemical flux observed by the UAV needs to be calculated by looking up the corresponding solar altitude angle and solar azimuth angle using the latitude and longitude obtained from the GPS on the UAV, and then combined with the UAV's IMU body attitude data (the angles of the body in three dimensions) to perform photochemical flux cosine correction. The specific correction includes:

[0067] Based on the solar altitude angle, solar azimuth angle, and UAV three-dimensional attitude tilt angle when the photochemical flux is obtained, the cosine correction coefficient is determined.

[0068] Based on the cosine correction coefficient, the photochemical flux is cosine corrected to obtain the corrected photochemical flux, calculated using the following formula: ;in, This represents the corrected photochemical flux (the actual photochemical flux after correction, in mW⋅m). −3 (milliwatts per cubic meter) Photochemical flux (raw value obtained from spectrometer measurement, unit: mW⋅m) −3 (milliwatts per cubic meter) This represents the cosine correction factor, which is determined by the solar altitude angle, solar azimuth angle, and the three-dimensional attitude of the UAV.

[0069] ;

[0070] in: Indicates the solar altitude angle. Indicates the azimuth angle of the sun. Indicates the pitch angle of the drone. Indicates the roll angle of the drone. This indicates the yaw angle of the drone.

[0071] The core purpose of this correction is to adjust the values ​​measured by the UAV in any attitude to the values ​​that should be measured under the assumption that the sensor plane is always horizontal (i.e. parallel to the ground).

[0072] The solar altitude angle and azimuth angle can be calculated from ephemeris tables. The roll, pitch, and heading angles of the UAV are provided by its onboard inertial measurement unit (IMU) or integrated navigation system. It is worth noting that high-quality spectral measurement systems will have a cosine corrector (usually a diffuser) installed at the probe end. Its main function is to optimize the sensor's response to light at different incident angles, making it conform as closely as possible to the cosine law, thus laying a good physical foundation for subsequent mathematical corrections.

[0073] For each altitude, the corresponding aerosol and meteorological parameters are processed as follows:

[0074] S2. Fit the aerosol absorption coefficients at several wavelengths to obtain the aerosol absorption coefficients for continuous wavelengths. In this embodiment, the aerosol absorption coefficients for continuous wavelengths within the measurement range are determined. The wavelength range depends on matching the wavelength range of the photochemical flux.

[0075] First, the aerosol absorption coefficient at several wavelengths is fitted to the wavelength to obtain the functional relationship between the absorption coefficient and the wavelength; the functional relationship between the absorption coefficient and the wavelength is as follows: Where C is a constant value obtained after fitting, and AAE is the wavelength absorption index obtained after fitting. b abs ( λ () represents the aerosol absorption coefficient at wavelength λ;

[0076] Then, based on the functional relationship between the absorption coefficient and the wavelength, the aerosol absorption coefficient corresponding to the continuous wavelength (resolution 1 nm) in the measurement band range is calculated.

[0077] S3. Based on the photochemical flux and the aerosol absorption coefficient of the continuous wavelength, calculate the aerosol absorption rate. In this embodiment, the aerosol absorption rate (A) is determined within the measurement wavelength range. aerosol ,mW⋅m −3 (milliwatts per cubic meter). The measurement band range covers ultraviolet, visible, and infrared bands, with typical values ​​of 350~2500 nm (but not limited to), and the specific value range can be determined by the measurement personnel according to actual needs.

[0078] The formula for calculating the absorption rate of aerosols is:

[0079] ;

[0080] Among them: A areosol The aerosol absorption rate (mW / m) is the integral over a specific wavelength range. 3 (milliwatts per cubic meter) The aerosol absorption coefficient (Mm) at wavelength λ is expressed as follows. -1 (megameter), λ represents wavelength (nm, nanometer). This represents the photochemical flux at wavelength λ.

[0081] S4. Based on the meteorological parameters, calculate the air density and specific heat capacity at constant pressure.

[0082] The formulas for calculating air density and specific heat capacity are as follows:

[0083] ;

[0084] ;

[0085] ;

[0086] ;

[0087] Where: ρ represents the density of moist air, kg / m³ 3 The values ​​are: kilograms per cubic meter; P represents the measured atmospheric pressure of moist air; Pv represents the water vapor pressure; Pd represents the dry air pressure; T represents the temperature; RH represents the relative humidity, in percent; Pvs represents the saturated water vapor partial pressure corresponding to the measured temperature; Rd represents the dry air specific gas constant; Rv represents the moist air specific gas constant, which is 287.05 under normal temperature conditions (-253.15 K~323.15 K). Rv=461.5 ; q indicates specific humidity; C P C represents the specific heat capacity of moist air at constant pressure. Pd C represents the specific heat capacity of dry air at constant pressure. PV C represents the isobaric specific heat capacity of water vapor at room temperature (-253.15 K~323.15 K). Pd =1005 C PV =1860 .

[0088] The parameters observed by the meteorological sensors are temperature (T, K), relative humidity (RH, %), and atmospheric pressure (P, Pa). Other parameters used in the calculations are constants or intermediate parameters obtained from literature references.

[0089] Moist air density ( The density of moist air can be calculated using the ideal gas law and the composition of moist air. Moist air consists of dry air and water vapor, and its density is the total mass divided by the volume.

[0090] Calculation steps:

[0091] 1. Calculate the saturated water vapor pressure ( ): Calculate the saturated vapor pressure of water at a given temperature using empirical formulas (such as the Tetens formula):

[0092] .

[0093] 2. Calculate the actual partial pressure of water vapor ( Relative humidity (Expressed as a percentage, 0-100) is defined as: ,therefore:

[0094] .

[0095] 3. Calculate the partial pressure of dry air ( ):

[0096] .

[0097] 4. Substitute into the formula for the density of moist air :

[0098] .

[0099] Where: RH represents relative humidity (%); It indicates temperature in degrees Celsius (°C); Represents absolute temperature (K), that is ,in Temperature in degrees Celsius (°C); Indicates total air pressure. ; This indicates the partial pressure of dry air (Pa); This represents the partial pressure of water vapor (Pa). The gas constant representing dry air. ; The gas constant representing water vapor. .

[0100] Specific heat capacity of moist air ( The specific heat capacity of moist air is a weighted average of the specific heat capacities of dry air and water vapor, calculated by their mass fractions. The specific heat capacity of moist air is typically specified as its specific pressure heat capacity. .

[0101] Calculation steps:

[0102] 1. Calculate the saturated water vapor pressure ( ): Calculate the saturated vapor pressure of water at a given temperature using empirical formulas (such as the Tetens formula):

[0103] .

[0104] 2. Calculate the actual partial pressure of water vapor ( Relative humidity (Expressed as a percentage, 0-100) is defined as: ,therefore:

[0105] .

[0106] 3. Calculate the partial pressure of dry air ( ):

[0107] .

[0108] 4. Calculate absolute humidity (water vapor mass fraction) :

[0109] 。

[0110] 5. Substitute into the specific heat capacity formula ( ):

[0111] 。

[0112] in: The value represents temperature in degrees Celsius (°C); RH represents relative humidity (%). Indicates total air pressure. ; This indicates the partial pressure of dry air (Pa); This represents the partial pressure of water vapor (Pa). This indicates absolute humidity, which is the mass fraction of water vapor (unitless). The gas constant representing dry air. ; The gas constant representing water vapor. ; This represents the specific heat capacity of water vapor at constant pressure. ; This represents the specific heat capacity of dry air at constant pressure. .

[0113] S5. Calculate the atmospheric heating rate of the aerosol based on the aerosol absorption rate, the air density, and the specific heat capacity at constant pressure.

[0114] The calculation formula is:

[0115] ;

[0116] Wherein: H areosol A represents the rate at which aerosols heat the atmosphere (K / d, K / day); areosol Indicates the aerosol absorption rate; ρ represents the air density; C P This indicates the specific heat capacity at constant pressure.

[0117] S6. Based on the aerosol heating rate to the atmosphere calculated at each altitude, determine the vertical distribution profile of the aerosol heating rate to the atmosphere as a function of altitude.

[0118] This embodiment calculates H at different heights, for example, at 1-meter intervals. areosol Ultimately, we can obtain H. areosol Vertical profiles that vary with height.

[0119] The aerosol heating rate corresponding to the measurement height is obtained, thereby determining the profile of the aerosol heating rate as a function of height. The accuracy of the vertical distribution profile of the heating rate depends on the photochemical flux measuring instrument, the aerosol absorption coefficient measuring instrument, the meteorological sensor, and the hovering accuracy of the UAV.

[0120] The resolution of the vertical distribution profile of aerosol heating rate depends on the ascent and descent rate of the UAV and the data refresh rate of the observation instrument. Under the current technical conditions, it is possible to obtain a vertical distribution profile of aerosol heating rate with a resolution of at least 1m.

[0121] In one exemplary embodiment, an apparatus is provided for obtaining a vertical distribution profile of the atmospheric heating rate of aerosols, comprising: an unmanned aerial vehicle (UAV) platform and measuring instruments; the UAV platform supports rotary-wing or fixed-wing aircraft.

[0122] The measuring instrument integrates a photochemical flux measuring instrument, a multi-wavelength aerosol absorption coefficient measuring instrument, a meteorological sensor, and a processor.

[0123] The photochemical flux measuring instrument can measure wavelengths covering ultraviolet, visible, and infrared bands, with typical values ​​ranging from 350 to 2500 nm (but not limited to). The specific range can be determined by the measuring personnel according to actual needs, and it is equipped with a hemispherical omnidirectional probe.

[0124] The multi-wavelength aerosol absorption coefficient measuring instrument can measure wavelengths covering the ultraviolet, visible, and infrared bands. The typical reference wavelength range is at least five wavelengths within 370~950nm (e.g., including but not limited to 370nm, 430nm, 470nm, 520nm, 565nm, 590nm, 660nm, 700nm, 880nm, and 950nm).

[0125] The meteorological sensor can measure temperature (T, K), relative humidity (RH, %), and atmospheric pressure (P, Pa) to calculate air density (ρ, kg / m³). 3 ) and specific heat capacity at constant pressure (Cp).

[0126] The processor is used for:

[0127] Acquire a number of aerosol parameters and meteorological parameters corresponding to several altitudes; the aerosol parameters include photochemical flux and aerosol absorption coefficients at several wavelengths.

[0128] For each altitude, the corresponding aerosol and meteorological parameters are processed as follows:

[0129] By fitting the aerosol absorption coefficients at several wavelengths, the aerosol absorption coefficients at continuous wavelengths can be obtained.

[0130] The aerosol absorption rate is calculated based on the photochemical flux and the aerosol absorption coefficient for the continuous wavelength.

[0131] Based on the meteorological parameters, calculate the air density and specific heat capacity at constant pressure.

[0132] The atmospheric heating rate of the aerosol is calculated based on the aerosol's absorption rate, the air density, and the specific heat capacity at constant pressure.

[0133] Based on the calculated aerosol heating rate to the atmosphere at each altitude, the vertical distribution profile of the aerosol heating rate to the atmosphere as a function of altitude is determined.

[0134] In one exemplary embodiment, a computer device is provided, which may be a server or a terminal, and its internal structure diagram may be as follows. Figure 3As shown, the computer device includes a processor, memory, input / output (I / O) interfaces, and a communication interface. The processor, memory, and I / O interfaces are connected via a system bus, and the communication interface is also connected to the system bus via the I / O interfaces. The processor provides computational and control capabilities. The memory includes non-volatile storage media and internal memory. The non-volatile storage media stores the operating system, computer programs, and a database. The internal memory provides the environment for the operation of the operating system and computer programs stored in the non-volatile storage media. The I / O interfaces are used for exchanging information between the processor and external devices. The communication interface is used for communicating with external terminals via a network connection. When executed by the processor, the computer program implements a method for obtaining the vertical distribution profile of the atmospheric heating rate of aerosols.

[0135] Those skilled in the art will understand that Figure 3 The structure shown is merely a block diagram of a portion of the structure related to the present application and does not constitute a limitation on the computer device to which the present application is applied. Specific computer devices may include more or fewer components than those shown in the figure, or combine certain components, or have different component arrangements.

[0136] In one exemplary embodiment, a computer device is also provided, including a memory and a processor, wherein the memory stores a computer program, and the processor executes the computer program to implement the steps in the above-described method embodiments.

[0137] In one exemplary embodiment, a computer-readable storage medium is provided storing a computer program that, when executed by a processor, implements the steps in the above-described method embodiments.

[0138] In one exemplary embodiment, a computer program product is provided, including a computer program that, when executed by a processor, implements the steps in the above-described method embodiments.

[0139] It should be noted that the user information (including but not limited to user device information, user personal information, etc.) and data (including but not limited to data used for analysis, data stored, data displayed, etc.) involved in this application are all information and data authorized by the user or fully authorized by all parties, and the collection, use and processing of the relevant data must comply with relevant regulations.

[0140] Those skilled in the art will understand that all or part of the processes in the above embodiments can be implemented by a computer program instructing related hardware. The computer program can be stored in a non-volatile computer-readable storage medium. When executed, the computer program can include the processes of the embodiments described above. Any references to memory, databases, or other media used in the embodiments provided in this application can include at least one of non-volatile and volatile memory. Non-volatile memory can include read-only memory (ROM), magnetic tape, floppy disk, flash memory, optical memory, high-density embedded non-volatile memory, resistive random access memory (ReRAM), magnetic random access memory (MRAM), ferroelectric random access memory (FRAM), phase change memory (PCM), graphene memory, etc. Volatile memory can include random access memory (RAM) or external cache memory, etc. By way of illustration and not limitation, RAM can take many forms, such as Static Random Access Memory (SRAM) or Dynamic Random Access Memory (DRAM).

[0141] The databases involved in the embodiments provided in this application may include at least one type of relational database and non-relational database. Non-relational databases may include, but are not limited to, blockchain-based distributed databases. The processors involved in the embodiments provided in this application may be general-purpose processors, central processing units, graphics processing units, digital signal processors, programmable logic devices, quantum computing-based data processing logic devices, etc., and are not limited to these.

[0142] The technical features of the above embodiments can be combined in any way. For the sake of brevity, not all possible combinations of the technical features in the above embodiments are described. However, as long as there is no contradiction in the combination of these technical features, they should be considered to be within the scope of this specification.

[0143] This document uses specific examples to illustrate the principles and implementation methods of this application. The descriptions of the above embodiments are only for the purpose of helping to understand the methods and core ideas of this application. Furthermore, those skilled in the art will recognize that, based on the ideas of this application, there will be changes in the specific implementation methods and application scope. Therefore, the content of this specification should not be construed as a limitation of this application.

Claims

1. A method for obtaining the vertical distribution profile of the atmospheric heating rate of aerosols, characterized in that, include: Acquire a number of aerosol parameters and meteorological parameters corresponding to several altitudes; the aerosol parameters include: photochemical flux and aerosol absorption coefficients at several wavelengths; For each altitude, the corresponding aerosol and meteorological parameters are processed as follows: By fitting the aerosol absorption coefficients at several wavelengths, the aerosol absorption coefficients at continuous wavelengths are obtained. The aerosol absorption rate is calculated based on the photochemical flux and the aerosol absorption coefficient of the continuous wavelength. Based on the meteorological parameters, calculate the air density and specific heat capacity at constant pressure; The atmospheric heating rate of the aerosol is calculated based on the aerosol's absorption rate, the air density, and the specific heat capacity at constant pressure. Based on the aerosol heating rate to the atmosphere calculated at each altitude, determine the vertical distribution profile of the aerosol heating rate to the atmosphere as a function of altitude. The process of fitting the aerosol absorption coefficients at several wavelengths to obtain the aerosol absorption coefficients at continuous wavelengths specifically includes: By fitting the aerosol absorption coefficient to the wavelength at several wavelengths, a functional relationship between the absorption coefficient and the wavelength is obtained; the functional relationship between the absorption coefficient and the wavelength is as follows: Where C is a constant value obtained after fitting, and AAE is the wavelength absorption index obtained after fitting. b abs ( λ () represents the aerosol absorption coefficient at wavelength λ; Based on the functional relationship between the absorption coefficient and the wavelength, the aerosol absorption coefficient for continuous wavelengths is calculated. The calculation of the aerosol absorption rate based on the photochemical flux and the aerosol absorption coefficient of the continuous wavelength specifically includes: Based on the solar elevation angle, solar azimuth angle and UAV attitude three-dimensional tilt angle when the photochemical flux is obtained, the cosine correction coefficient is determined; Based on the cosine correction coefficient, the photochemical flux is cosine corrected to obtain the corrected photochemical flux, calculated using the following formula: ;in, This indicates the correction for photochemical flux. Represents photochemical flux. This represents the cosine correction factor; The aerosol absorption rate is calculated based on the corrected photochemical flux and the aerosol absorption coefficient of the continuous wavelength.

2. The method for obtaining the vertical distribution profile of the atmospheric heating rate of aerosols according to claim 1, characterized in that, Based on the photochemical flux and the aerosol absorption coefficient for the continuous wavelength, the formula for calculating the aerosol absorption rate is as follows: ； Among them: A areosol Indicates the aerosol absorption rate. This represents the aerosol absorption coefficient at wavelength λ, where λ represents the wavelength. This represents the photochemical flux at wavelength λ.

3. The method for obtaining the vertical distribution profile of the atmospheric heating rate of aerosols according to claim 1, characterized in that, The formula for calculating air density is: ； in, ρ represents the density of moist air; P represents the atmospheric pressure of moist air; Pv represents the water vapor pressure; Rd represents the dry air specific gas constant; T represents the temperature; Rv represents the moist air specific gas constant; RH represents the relative humidity; Pvs represents the saturated water vapor partial pressure corresponding to the measured temperature. The formula for calculating the isobaric specific heat capacity is: ； in, ; q represents specific humidity; Pd represents dry air pressure; C P C represents the specific heat capacity of moist air at constant pressure. Pd C represents the specific heat capacity of dry air at constant pressure. PV This represents the specific heat capacity of water vapor at constant pressure.

4. The method for obtaining the vertical distribution profile of the atmospheric heating rate of aerosols according to claim 1, characterized in that, Based on the aerosol's absorption rate, the air density, and the isobaric specific heat capacity, the formula for calculating the aerosol's atmospheric heating rate is as follows: ； Wherein: H areosol A represents the rate at which aerosols heat the atmosphere. areosol The aerosol absorption rate is represented by ρ, and the air density is represented by C. P This indicates the specific heat capacity at constant pressure.

5. An apparatus for obtaining the vertical distribution profile of the atmospheric heating rate of aerosols, characterized in that, include: Unmanned aerial vehicle platform, photochemical flux measuring instrument, multi-wavelength aerosol absorption coefficient measuring instrument, meteorological sensor and processor; The drone platform supports rotary-wing or fixed-wing aircraft. The photochemical flux measuring instrument is equipped with a hemispherical omnidirectional probe; The photochemical flux measuring instrument can measure the following wavelength ranges: ultraviolet, visible light, and infrared. The multi-wavelength aerosol absorption coefficient measuring instrument is used to measure at least five wavelengths in the ultraviolet, visible, and infrared bands. The meteorological sensor is used to measure relative humidity, temperature, and atmospheric pressure. The processor is used for: Acquire a number of aerosol parameters and meteorological parameters corresponding to several altitudes; the aerosol parameters include: photochemical flux and aerosol absorption coefficients at several wavelengths; For each altitude, the corresponding aerosol and meteorological parameters are processed as follows: By fitting the aerosol absorption coefficients at several wavelengths, the aerosol absorption coefficients at continuous wavelengths are obtained. The aerosol absorption rate is calculated based on the photochemical flux and the aerosol absorption coefficient of the continuous wavelength. Based on the meteorological parameters, calculate the air density and specific heat capacity at constant pressure; The atmospheric heating rate of the aerosol is calculated based on the aerosol's absorption rate, the air density, and the specific heat capacity at constant pressure. Based on the aerosol heating rate to the atmosphere calculated at each altitude, determine the vertical distribution profile of the aerosol heating rate to the atmosphere as a function of altitude. The process of fitting the aerosol absorption coefficients at several wavelengths to obtain the aerosol absorption coefficients at continuous wavelengths specifically includes: By fitting the aerosol absorption coefficient to the wavelength at several wavelengths, a functional relationship between the absorption coefficient and the wavelength is obtained; the functional relationship between the absorption coefficient and the wavelength is as follows: Where C is a constant value obtained after fitting, and AAE is the wavelength absorption index obtained after fitting. b abs ( λ () represents the aerosol absorption coefficient at wavelength λ; Based on the functional relationship between the absorption coefficient and the wavelength, the aerosol absorption coefficient for continuous wavelengths is calculated. The calculation of the aerosol absorption rate based on the photochemical flux and the aerosol absorption coefficient of the continuous wavelength specifically includes: Based on the solar elevation angle, solar azimuth angle and UAV attitude three-dimensional tilt angle when the photochemical flux is obtained, the cosine correction coefficient is determined; Based on the cosine correction coefficient, the photochemical flux is cosine corrected to obtain the corrected photochemical flux, calculated using the following formula: ;in, This indicates the correction for photochemical flux. Represents photochemical flux. This represents the cosine correction factor; The aerosol absorption rate is calculated based on the corrected photochemical flux and the aerosol absorption coefficient of the continuous wavelength.

6. A computer device, comprising: A memory, a processor, and a computer program stored in the memory and capable of running on the processor, characterized in that the processor executes the computer program to implement the method for obtaining the vertical distribution profile of the aerosol heating rate to the atmosphere as claimed in any one of claims 1-4.

7. A computer-readable storage medium having a computer program stored thereon, characterized in that, When executed by a processor, the computer program implements the method for obtaining the vertical distribution profile of the aerosol heating rate to the atmosphere as described in any one of claims 1-4.

8. A computer program product, comprising a computer program, characterized in that, When executed by a processor, the computer program implements the method for obtaining the vertical distribution profile of the aerosol heating rate to the atmosphere as described in any one of claims 1-4.

Images