A method and device for determining the boundary value of aerosol extinction coefficient

Abstract

The present invention discloses a method and device for determining the boundary value of an aerosol extinction coefficient. The method comprises: constructing a calculation function for the boundary value of an aerosol extinction coefficient; iteratively calculating the calculation function according to a pre-constructed iterative formula to determine the upper boundary value of the aerosol extinction coefficient in the height interval to be calculated, wherein the iterative formula includes an adjustment factor; and calculating the boundary value of the aerosol extinction coefficient at each layer of the height interval to be calculated layer by layer according to the upper boundary value of the aerosol extinction coefficient in the height interval to be calculated and the Fernald inversion algorithm.

Description

Technical Field

[0001] The present invention relates to the field of atmosphere monitoring technology, and more specifically, to a method and device for determining a boundary value of an aerosol extinction coefficient. Background Art

[0002] Micro-pulse lidar can provide continuous and automatic observations of the optical properties and spatial distribution of clouds and aerosols, making it gradually become the main remote sensing tool for atmospheric aerosol research. In the process of inverting the aerosol extinction coefficient, the commonly used atmospheric extinction coefficient inversion algorithms are the Collis slope method, the Klett method, and the Fernald method. When the atmosphere is assumed to be uniformly distributed, the Collis slope method can be used to invert the extinction coefficient. But in fact, uniformly distributed atmospheric conditions do not occur, so most of the Klett method or Fernald method is used to solve the extinction coefficient. For these two algorithms, the boundary value of the extinction coefficient at different times and under different weather conditions is an uncertain value. Therefore, the key link in accurately inverting the aerosol extinction coefficient is the selection of the boundary value of the extinction coefficient. How to select the boundary value of the extinction coefficient has become a technical problem that needs to be solved urgently. Summary of the invention

[0003] In view of the deficiencies in the prior art, the present invention provides a method and device for determining a boundary value of an aerosol extinction coefficient.

[0004] According to one aspect of the present invention, a method for determining a boundary value of an aerosol extinction coefficient is provided, comprising:

[0005] Construct the calculation function of the boundary value of aerosol extinction coefficient;

[0006] Iteratively calculate the calculation function according to a pre-constructed iterative formula to determine the upper boundary value of the aerosol extinction coefficient in the height range to be calculated, wherein the iterative formula includes an adjustment factor;

[0007] According to the upper boundary value of the aerosol extinction coefficient in the height interval to be calculated and the Fernald inversion algorithm, the boundary value of the aerosol extinction coefficient at each height in the height interval to be calculated is calculated layer by layer.

[0008] Optionally, a calculation function for the boundary value of the aerosol extinction coefficient is constructed, including:

[0009] Select the boundary layer bottom reference point The height range where the boundary value of the aerosol extinction coefficient does not change dramatically ;

[0010] Calculate the height interval The mean aerosol extinction coefficient of

[0011] A calculation function for the boundary value of the aerosol extinction coefficient is constructed according to the mean value of the aerosol extinction coefficient.

[0012] Optionally, the height interval The formula for determining is:

[0013] ,

[0014] In the formula, The value is -15; is the distance resolution.

[0015] Optionally, the calculation function of the aerosol coefficient boundary value is:

[0016] ,

[0017] in,

[0018] ,

[0019] In the formula, x is the boundary value of aerosol extinction coefficient, For height range The mean aerosol extinction coefficient of D For height range Distance resolution Divided into equal intervals, ,in Indicates rounding down; height interval Height inside , d is the interval height, d =1,2,3,..., D， For height The aerosol extinction coefficient at is given by the boundary value of the aerosol extinction coefficient x and height Decide.

[0020] Optionally, the iteration formula is:

[0021] ,

[0022] In the formula, is the regulating factor, For the k +1 aerosol extinction coefficient boundary value; For the k The aerosol extinction coefficient boundary value of the iteration; For the k The height interval of the iteration The mean aerosol extinction coefficient of For the k The result of the aerosol coefficient boundary value calculation function at the iteration.

[0023] Optionally, the adjustment factor The value selection strategy is as follows:

[0024] The regulatory factor Perform initial assignment ;

[0025] In each iteration, the adjustment factor is adjusted according to the update effect of the current aerosol extinction coefficient boundary value size;

[0026] When the iteration end condition is met, the iteration ends and the boundary value of the aerosol extinction coefficient is determined.

[0027] Optionally, the adjustment factor is adjusted according to the update effect of the current aerosol extinction coefficient boundary value Size, including:

[0028] If the iteration produces a better solution for the aerosol extinction coefficient boundary value, then Reduce the adjustment factor ,in k is the current iteration number, the first adjustment coefficient ;

[0029] If the iteration result does not provide a better solution for the boundary value of the aerosol extinction coefficient, then Increase the regulation factor , where the second adjustment coefficient .

[0030] Optionally, the iteration end condition is:

[0031] or or ,

[0032] In the formula, is the pre-set iteration accuracy; K is the maximum number of iterations set in advance; k is the current iteration number.

[0033] Optionally, the aerosol extinction coefficient boundary value of each layer in the height interval to be calculated is calculated layer by layer according to the upper boundary value of the aerosol extinction coefficient in the height interval to be calculated and the Fernald inversion algorithm, including:

[0034] According to the upper boundary value of the aerosol extinction coefficient in the height interval to be calculated and the echo signal intensity at the upper boundary of the height interval to be calculated, the Fernald inversion algorithm is used to perform integral calculation layer by layer from top to bottom to calculate the boundary value of the aerosol extinction coefficient at each height in the height interval to be calculated.

[0035] According to another aspect of the present invention, there is provided a device for determining a boundary value of an aerosol extinction coefficient, comprising:

[0036] A construction module is used to construct a calculation function for the boundary value of the aerosol extinction coefficient;

[0037] An iteration module, used for iteratively calculating the calculation function according to a pre-constructed iteration formula to determine the upper boundary value of the aerosol extinction coefficient in the height interval to be calculated, wherein the iteration formula includes an adjustment factor;

[0038] The calculation module is used to calculate the boundary value of the aerosol extinction coefficient at each height in the height interval to be calculated layer by layer according to the upper boundary value of the aerosol extinction coefficient in the height interval to be calculated and the Fernald inversion algorithm.

[0039] According to another aspect of the present invention, a computer-readable storage medium is provided, wherein the storage medium stores a computer program, and the computer program is used to execute the method described in any one of the above aspects of the present invention.

[0040] According to another aspect of the present invention, an electronic device is provided, comprising: a processor; a memory for storing instructions executable by the processor; the processor is configured to read the executable instructions from the memory and execute the instructions to implement the method described in any one of the above aspects of the present invention.

[0041] Therefore, the present application finds a balance between convergence speed and stability by introducing a dynamic adjustment factor ω in the iterative calculation of the boundary value of the aerosol correlation coefficient to dynamically adjust the iterative result, thereby improving the efficiency and accuracy of the solved aerosol extinction coefficient boundary value, and then solving the aerosol extinction coefficient of the entire layer. BRIEF DESCRIPTION OF THE DRAWINGS

[0042] A more complete understanding of exemplary embodiments of the present invention may be obtained by referring to the following drawings:

[0043] Figure 1 is a flow chart of a method for determining a boundary value of an aerosol extinction coefficient provided by an exemplary embodiment of the present invention;

[0044] Figure 2 It is a structural schematic diagram of a device for determining a boundary value of an aerosol extinction coefficient provided by an exemplary embodiment of the present invention;

[0045] Figure 3 This is a structure of an electronic device provided by an exemplary embodiment of the present invention. DETAILED DESCRIPTION

[0046] Below, the exemplary embodiments according to the present invention will be described in detail with reference to the accompanying drawings. Obviously, the described embodiments are only part of the embodiments of the present invention, rather than all the embodiments of the present invention, and it should be understood that the present invention is not limited to the exemplary embodiments described here.

[0047] It should be noted that the relative arrangement of components and steps, the numerical expressions and numerical values ​​set forth in these embodiments do not limit the scope of the present invention unless specifically stated otherwise.

[0048] Those skilled in the art can understand that the terms "first" and "second" in the embodiments of the present invention are only used to distinguish different steps, devices or modules, etc., and neither represent any specific technical meaning nor indicate the necessary logical order between them.

[0049] It should also be understood that, in the embodiments of the present invention, “plurality” may refer to two or more than two, and “at least one” may refer to one, two or more than two.

[0050] It should also be understood that any component, data or structure mentioned in the embodiments of the present invention can generally be understood as one or more, unless explicitly limited or otherwise indicated in the context.

[0051] In addition, the term "and / or" in the present invention is only a description of the association relationship of the associated objects, indicating that there can be three relationships. For example, A and / or B can represent: A exists alone, A and B exist at the same time, and B exists alone. In addition, the character " / " in the present invention generally indicates that the associated objects before and after are in an "or" relationship.

[0052] It should also be understood that the description of the various embodiments of the present invention focuses on the differences between the various embodiments, and the same or similar aspects thereof can be referenced to each other, and for the sake of brevity, they will not be described one by one.

[0053] At the same time, it should be understood that for the convenience of description, the sizes of the various parts shown in the drawings are not drawn according to the actual proportional relationship.

[0054] The following description of at least one exemplary embodiment is merely illustrative in nature and is in no way intended to limit the invention, its application, or uses.

[0055] Technologies, methods, and equipment known to ordinary technicians in the relevant art may not be discussed in detail, but where appropriate, the technologies, methods, and equipment should be considered part of the specification.

[0056] It should be noted that like reference numerals and letters refer to similar items in the following figures, and therefore, once an item is defined in one figure, it need not be further discussed in subsequent figures.

[0057] Embodiments of the present invention can be applied to electronic devices such as terminal devices, computer systems, servers, etc., which can operate with many other general or special computing system environments or configurations. Examples of well-known terminal devices, computing systems, environments and / or configurations suitable for use with electronic devices such as terminal devices, computer systems, servers, etc. include, but are not limited to: personal computer systems, server computer systems, thin clients, thick clients, handheld or laptop devices, microprocessor-based systems, set-top boxes, programmable consumer electronics, network personal computers, small computer systems, large computer systems, and distributed cloud computing technology environments including any of the above systems, etc.

[0058] Electronic devices such as terminal devices, computer systems, servers, etc. can be described in the general context of computer system executable instructions (such as program modules) executed by computer systems. Generally, program modules can include routines, programs, object programs, components, logic, data structures, etc., which perform specific tasks or implement specific abstract data types. Computer systems / servers can be implemented in a distributed cloud computing environment, where tasks are performed by remote processing devices linked through a communication network. In a distributed cloud computing environment, program modules can be located on local or remote computing system storage media including storage devices.

[0059] Exemplary Methods

[0060] Figure 1 FIG. 1 is a flow chart of a method for determining a boundary value of an aerosol extinction coefficient provided by an exemplary embodiment of the present invention. This embodiment can be applied to electronic devices, such as Figure 1 As shown, the method 100 for determining the boundary value of the aerosol extinction coefficient includes the following steps:

[0061] Step 101, constructing a calculation function for the boundary value of the aerosol extinction coefficient;

[0062] Step 102, iteratively calculating the calculation function according to a pre-constructed iterative formula to determine the upper boundary value of the aerosol extinction coefficient in the height interval to be calculated, wherein the iterative formula includes an adjustment factor;

[0063] Step 103, calculating the boundary value of the aerosol extinction coefficient at each height in the height interval to be calculated layer by layer according to the upper boundary value of the aerosol extinction coefficient in the height interval to be calculated and the Fernald inversion algorithm.

[0064] Specifically, in view of the technical problems existing in the background technology, the present application proposes a method for determining the boundary value of the aerosol extinction coefficient to accurately obtain the boundary value of the aerosol extinction coefficient. The specific implementation steps are as follows:

[0065] (1) Construct the calculation function of the boundary value of the aerosol extinction coefficient: Assume that the reference point at the bottom of the boundary layer is The boundary value of the aerosol extinction coefficient at is x, Select a height interval near the aerosol extinction coefficient where there is no drastic change ;

[0066] ,

[0067] Value -15, height interval The mean aerosol extinction coefficient within the range is:

[0068] ,

[0069] Where D is According to the distance resolution within the interval Divided into equal intervals, ,(in indicates rounding down); The constructed equation is:

[0070] .

[0071] (2) By introducing the adjustment factor ω, the result can be avoided from fluctuating. The iterative formula is:

[0072] ,

[0073] Among them, the value strategy of ω is:

[0074] 1) Initial value ω=0.01;

[0075] 2) In each iteration, adjust the size of ω according to the current update effect;

[0076] 3) In each iteration, if the aerosol extinction coefficient boundary value of the current iteration is better than the boundary value of the previous iteration, the adjustment factor is updated as:

[0077] ,

[0078] where γ1∈(0,0.1] (e.g., γ1=0.05) is used to reduce the value of the adjustment factor in order to more finely tune the model parameters.

[0079] If the aerosol extinction coefficient boundary value of the current iteration is not better than the previous boundary value, the updated adjustment factor is:

[0080] ,

[0081] Among them, γ2∈(1,1.1] (for example, γ2=1.05) is used to increase the value of the adjustment factor to avoid the instability of the method and make the update process closer to the stability of the gradient descent method.

[0082] (3) Generate an iterative sequence according to (2) { }, the condition for determining the end of iteration is:

[0083] or or ,

[0084] Among them, ε is the pre-set iteration accuracy; K is the pre-set maximum number of iterations; each time an iteration is performed, the iteration end condition is judged and the ω balance speed and stability are dynamically adjusted. If the conditions are met, the iteration is stopped. After the iteration is completed, the boundary value of the aerosol extinction coefficient can be obtained;

[0085] (4) According to the above steps, we can get the upper boundary Extinction coefficient at .initialization: , ) and echo signal strength ); Integral calculation: from the upper boundary Starting from arrive interval, and calculate step by step downward; step-by-step inversion: use Fernald formula to calculate each Extinction coefficient at , iterate step by step, each time The value at the bottom reference value.

[0086] According to the Fernald inversion algorithm, the aerosol extinction coefficient at the boundary height Zc is determined as:

[0087] ,

[0088] In the formula, , , and are the extinction coefficients of aerosol and molecule at height z, respectively, and are the backscattering coefficients of aerosol and molecules at height z, is the atmospheric aerosol extinction backscattering ratio parameter that varies with altitude. In the prior art, it is generally assumed to be a quantity that does not vary with altitude. ; is the extinction backscattering ratio of atmospheric air molecules that varies with altitude. In the prior art, it is generally assumed to be a quantity that does not vary with altitude. ; is the distance square correction signal, , is the echo signal received by the laser radar at distance z, which can be expressed as: , where E is the output energy; C is the calibration constant of the instrument; is the atmospheric transmittance of the molecule, ; is the atmospheric transmittance considering only aerosol attenuation, . is the boundary value of the aerosol extinction coefficient at the reference height Zc.

[0089] Therefore, the present application finds a balance between convergence speed and stability by introducing a dynamic adjustment factor ω in the iterative calculation of the boundary value of the aerosol correlation coefficient to dynamically adjust the iterative result, thereby improving the efficiency and accuracy of the solved aerosol extinction coefficient boundary value, and then solving the aerosol extinction coefficient of the entire layer.

[0090] Exemplary Devices

[0091] Figure 2 FIG. 1 is a schematic diagram of a device for determining an aerosol extinction coefficient boundary value provided by an exemplary embodiment of the present invention. Figure 2 As shown, the device 200 includes:

[0092] A construction module 210 is used to construct a calculation function of an aerosol extinction coefficient boundary value;

[0093] An iteration module 220 is used to iteratively calculate the calculation function according to a pre-constructed iteration formula to determine the upper boundary value of the aerosol extinction coefficient in the height interval to be calculated, wherein the iteration formula includes an adjustment factor;

[0094] The calculation module 230 is used to calculate the boundary value of the aerosol extinction coefficient at each height in the height interval to be calculated layer by layer according to the upper boundary value of the aerosol extinction coefficient in the height interval to be calculated and the Fernald inversion algorithm.

[0095] Optionally, the construction module 210 includes:

[0096] The selection submodule is used to select the height interval where the boundary value of the aerosol extinction coefficient at the reference point Zc at the bottom of the boundary layer does not change dramatically ;

[0097] The first calculation submodule is used to calculate the height interval The mean aerosol extinction coefficient of

[0098] A submodule is constructed to construct a calculation function for the boundary value of the aerosol extinction coefficient based on the mean value of the aerosol extinction coefficient.

[0099] Optionally, the height interval The formula for determining is:

[0100]

[0101] In the formula, The value is -15; is the distance resolution.

[0102] Optionally, the calculation function of the aerosol coefficient boundary value is:

[0103] ,

[0104] in,

[0105] ,

[0106] Where x is the boundary value of aerosol extinction coefficient, For height range The mean value of aerosol extinction coefficient; D is the height range Distance resolution Divided into equal intervals, ,in Indicates rounding down; height interval Height inside , d is the interval height, d=1,2,3,...,D , For height The aerosol extinction coefficient at is given by the boundary value of the aerosol extinction coefficient x and height Decide.

[0107] Optionally, the iteration formula is:

[0108] ,

[0109] In the formula, is the regulating factor, is the boundary value of the aerosol extinction coefficient at the k+1th time; is the boundary value of aerosol extinction coefficient of the kth iteration; is the height interval at the kth iteration The mean aerosol extinction coefficient of The result of the aerosol coefficient boundary value calculation function at the kth iteration.

[0110] Optionally, the adjustment factor The value selection strategy is as follows:

[0111] The regulatory factor Perform initial assignment ;

[0112] In each iteration, the adjustment factor is adjusted according to the update effect of the current aerosol extinction coefficient boundary value size;

[0113] When the iteration end condition is met, the iteration ends and the boundary value of the aerosol extinction coefficient is determined.

[0114] Optionally, adjust the adjustment factor according to the update effect of the current aerosol extinction coefficient boundary value Size, including:

[0115] If the iteration produces a better solution for the aerosol extinction coefficient boundary value, then Reduce the adjustment factor ω, where k is the current iteration number, the first adjustment coefficient ;

[0116] If the iteration result does not provide a better solution for the boundary value of the aerosol extinction coefficient, then Increase the adjustment factor ω, where the second adjustment coefficient .

[0117] Optionally, the iteration end condition is:

[0118] or or ,

[0119] Where ε is the preset iteration accuracy; K is the preset maximum number of iterations; and k is the current number of iterations.

[0120] Optionally, the calculation module 230 includes:

[0121] The second calculation submodule is used to calculate the boundary value of the aerosol extinction coefficient at each height in the height interval to be calculated by using the Fernald inversion algorithm to perform integral calculation from top to bottom layer by layer according to the upper boundary value of the aerosol extinction coefficient in the height interval to be calculated and the echo signal strength at the upper boundary of the height interval to be calculated.

[0122] Exemplary Electronic Devices

[0123] Figure 3 This is a structure of an electronic device provided by an exemplary embodiment of the present invention. Figure 3As shown, the electronic device 30 includes one or more processors 31 and a memory 32 .

[0124] The processor 31 may be a central processing unit (CPU) or other forms of processing units having data processing capabilities and / or instruction execution capabilities, and may control other components in the electronic device to perform desired functions.

[0125] The memory 32 may include one or more computer program products, which may include various forms of computer-readable storage media, such as volatile memory and / or non-volatile memory. The volatile memory may include, for example, a random access memory (RAM) and / or a cache memory (cache), etc. The non-volatile memory may include, for example, a read-only memory (ROM), a hard disk, a flash memory, etc. One or more computer program instructions may be stored on the computer-readable storage medium, and the processor 31 may run the program instructions to implement the methods of the software programs of the various embodiments of the present invention described above and / or other desired functions. In one example, the electronic device may also include: an input device 33 and an output device 34, which are interconnected via a bus system and / or other forms of connection mechanisms (not shown).

[0126] In addition, the input device 33 may also include, for example, a keyboard, a mouse, and the like.

[0127] The output device 34 can output various information to the outside. The output device 34 can include, for example, a display, a speaker, a printer, a communication network and a remote output device connected thereto.

[0128] Of course, to simplify, Figure 3 Only some of the components related to the present invention in the electronic device are shown, and components such as a bus, an input / output interface, etc. are omitted. In addition, the electronic device may further include any other appropriate components according to specific application conditions.

[0129] Exemplary computer program products and computer-readable storage media

[0130] In addition to the above-mentioned methods and devices, an embodiment of the present invention may also be a computer program product, which includes computer program instructions, which, when executed by a processor, enable the processor to execute the steps of the method according to various embodiments of the present invention described in the above-mentioned "Exemplary Method" section of this specification.

[0131] The computer program product may be written in any combination of one or more programming languages ​​to write program code for performing the operations of the embodiments of the present invention, including object-oriented programming languages ​​such as Java, C++, etc., and conventional procedural programming languages ​​such as "C" or similar programming languages. The program code may be executed entirely on the user computing device, partially on the user device, as a separate software package, partially on the user computing device and partially on a remote computing device, or entirely on a remote computing device or server.

[0132] In addition, an embodiment of the present invention may also be a computer-readable storage medium having computer program instructions stored thereon, which, when executed by a processor, enable the processor to execute the steps of the method according to various embodiments of the present invention described in the above “Exemplary Method” section of this specification.

[0133] The computer readable storage medium can adopt any combination of one or more readable media. The readable medium can be a readable signal medium or a readable storage medium. The readable storage medium can include, for example, but is not limited to, a system, system or device of electricity, magnetism, light, electromagnetic, infrared, or semiconductor, or any combination of the above. More specific examples (non-exhaustive list) of readable storage media include: an electrical connection with one or more wires, a portable disk, a hard disk, a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), an optical fiber, a portable compact disk read-only memory (CD-ROM), an optical storage device, a magnetic storage device, or any suitable combination of the above.

[0134] The basic principle of the present invention is described above in conjunction with specific embodiments. However, it should be pointed out that the advantages, strengths, effects, etc. mentioned in the present invention are only examples and not limitations, and it cannot be considered that these advantages, strengths, effects, etc. must be possessed by each embodiment of the present invention. In addition, the specific details disclosed above are only for the purpose of illustration and facilitation of understanding, rather than limitation, and the above details do not limit the present invention to being implemented by adopting the above specific details.

[0135] Each embodiment in this specification is described in a progressive manner, and each embodiment focuses on the differences from other embodiments. The same or similar parts between the embodiments can be referred to each other. For the system embodiment, since it basically corresponds to the method embodiment, the description is relatively simple, and the relevant parts can be referred to the partial description of the method embodiment.

[0136] The block diagrams of the devices, systems, equipment, and systems involved in the present invention are only illustrative examples and are not intended to require or imply that they must be connected, arranged, and configured in the manner shown in the block diagram. As will be appreciated by those skilled in the art, these devices, systems, equipment, and systems can be connected, arranged, and configured in any manner. Words such as "including", "comprising", "having", etc. are open words, referring to "including but not limited to", and can be used interchangeably with them. The words "or" and "and" used here refer to the words "and / or" and can be used interchangeably with them, unless the context clearly indicates otherwise. The word "such as" used here refers to the phrase "such as but not limited to", and can be used interchangeably with it.

[0137] The method and system of the present invention may be implemented in many ways. For example, the method and system of the present invention may be implemented by software, hardware, firmware or any combination of software, hardware, firmware. The above order of steps for the method is only for illustration, and the steps of the method of the present invention are not limited to the order specifically described above, unless otherwise specifically stated. In addition, in some embodiments, the present invention may also be implemented as a program recorded in a recording medium, which includes machine-readable instructions for implementing the method according to the present invention. Thus, the present invention also covers a recording medium storing a program for executing the method according to the present invention.

[0138] It should also be noted that in the system, device and method of the present invention, each component or each step can be decomposed and / or recombined. These decompositions and / or recombinations should be regarded as equivalent schemes of the present invention. The above description of the disclosed aspects is provided to enable any technician in the field to make or use the present invention. Various modifications to these aspects are very obvious to those skilled in the art, and the general principles defined here can be applied to other aspects without departing from the scope of the present invention. Therefore, the present invention is not intended to be limited to the aspects shown here, but in accordance with the widest range consistent with the principles and novel features disclosed here.

[0139] The above description has been given for the purpose of illustration and description. In addition, this description is not intended to limit the embodiments of the present invention to the forms disclosed herein. Although a number of example aspects and embodiments have been discussed above, those skilled in the art will recognize certain variations, modifications, changes, additions and sub-combinations thereof.

Claims

1. A method for determining the boundary value of an aerosol extinction coefficient, characterized in that: include: Construct the calculation function of the boundary value of aerosol extinction coefficient; Iteratively calculating the calculation function according to a pre-constructed iterative formula to determine an upper boundary value of the aerosol extinction coefficient in the height interval to be calculated, wherein the iterative formula includes an adjustment factor; Calculate the boundary value of the aerosol extinction coefficient at each height in the height interval to be calculated layer by layer according to the upper boundary value of the aerosol extinction coefficient in the height interval to be calculated and the Fernald inversion algorithm; The iteration formula is: In the formula, ω is the adjustment factor, x k+1 is the boundary value of the aerosol extinction coefficient at the k+1th time; x k is the boundary value of aerosol extinction coefficient of the kth iteration; is the height interval (z c′ ,z c ) of the aerosol extinction coefficient; Φ(x k ) is the result of the aerosol coefficient boundary value calculation function at the kth iteration; The value strategy of the adjustment factor ω is as follows: The adjustment factor ω is initially assigned a value of ω=0.01; In each iteration, the magnitude of the adjustment factor ω is adjusted according to the update effect of the current aerosol extinction coefficient boundary value; When the iteration end condition is met, the iteration is ended, and the boundary value of the aerosol extinction coefficient is determined; Adjusting the size of the adjustment factor ω according to the update effect of the current aerosol extinction coefficient boundary value includes: If the iteration produces a better solution for the aerosol extinction coefficient boundary value, then Reduce the adjustment factor ω, where k is the current iteration number and the first adjustment coefficient γ1∈(0,0.1]; If the iteration result does not provide a better solution for the boundary value of the aerosol extinction coefficient, then Increase the adjustment factor ω, wherein the second adjustment coefficient γ2∈(1,1.1].

2. The method according to claim 1, characterized in that Construct the calculation function of the boundary value of the aerosol extinction coefficient, including: Select the reference point z at the bottom of the boundary layer c The height range where the boundary value of the aerosol extinction coefficient does not change dramatically (z c′ ,z c ); Calculate the height interval (z c′ ,z c )’s mean aerosol extinction coefficient; Constructing a calculation function of the boundary value of the aerosol extinction coefficient according to the mean value of the aerosol extinction coefficient; The height interval (z c′ ,z c ) is determined by: With c′ =from c +θ×Δz In the formula, θ is -15; Δz is the distance resolution.

3. The method according to claim 2, characterized in that The calculation function of the aerosol extinction coefficient boundary value is: in, Where x is the boundary value of aerosol extinction coefficient, is the height interval (z c′ ,z c ) is the mean value of the aerosol extinction coefficient; D is the height interval (z c′ ,z c ) is divided into equal intervals according to the distance resolution Δz, in Indicates rounding down; height interval (z c′ ,z c ) d =d×Δz+z c′ , d is the interval height, d=1,2,3,...,D,σ a (z d ,x) is the aerosol extinction coefficient at height zd, which is determined by the boundary value x of the aerosol extinction coefficient and the height z d Decide.

4. The method according to claim 1, characterized in that: The iteration end condition is: |x k+1 -x k |< ε or φ(x k ) < ε or k < K Where ε is the preset iteration accuracy; K is the preset maximum number of iterations; and k is the current number of iterations.

5. The method according to claim 1, characterized in that Calculating the boundary value of the aerosol extinction coefficient at each layer in the height interval to be calculated layer by layer according to the upper boundary value of the aerosol extinction coefficient in the height interval to be calculated and the Fernald inversion algorithm, including: According to the upper boundary value of the aerosol extinction coefficient in the height interval to be calculated and the echo signal intensity at the upper boundary of the height interval to be calculated, the Fernald inversion algorithm is used to perform integral calculation layer by layer from top to bottom to calculate the boundary value of the aerosol extinction coefficient at each height in the height interval to be calculated.

6. A device for determining a boundary value of an aerosol extinction coefficient, characterized in that: include: A construction module is used to construct a calculation function for the boundary value of the aerosol extinction coefficient; An iteration module, used for iteratively calculating the calculation function according to a pre-constructed iteration formula to determine the upper boundary value of the aerosol extinction coefficient in the height interval to be calculated, wherein the iteration formula includes an adjustment factor; A calculation module, used to calculate the boundary value of the aerosol extinction coefficient at each height in the height interval to be calculated layer by layer according to the upper boundary value of the aerosol extinction coefficient in the height interval to be calculated and the Fernald inversion algorithm; The iteration formula is: In the formula, ω is the adjustment factor, x k+1 is the boundary value of the aerosol extinction coefficient at the k+1th time; x k is the boundary value of aerosol extinction coefficient of the kth iteration; is the height interval (z c′ ,z c ) of the aerosol extinction coefficient; Φ(x k ) is the result of the aerosol coefficient boundary value calculation function at the kth iteration; The value strategy of the adjustment factor ω is as follows: The adjustment factor ω is initially assigned a value of ω=0.01; In each iteration, the magnitude of the adjustment factor ω is adjusted according to the update effect of the current aerosol extinction coefficient boundary value; When the iteration end condition is met, the iteration is ended, and the boundary value of the aerosol extinction coefficient is determined; Adjusting the size of the adjustment factor ω according to the update effect of the current aerosol extinction coefficient boundary value includes: If the iteration produces a better solution for the aerosol extinction coefficient boundary value, then Reduce the adjustment factor ω, where k is the current iteration number and the first adjustment coefficient γ1∈(0,0.1]; If the iteration result does not provide a better solution for the boundary value of the aerosol extinction coefficient, then Increase the adjustment factor ω, wherein the second adjustment coefficient γ2∈(1,1.1].

Images