Petrochemical park gas leakage monitoring and positioning method and system, electronic device and storage medium

Abstract

The application discloses a petrochemical park gas leakage monitoring and positioning method, comprising the following steps: determining the downwind boundary of the target park dominant wind direction as a monitoring boundary, and obtaining the light path integral concentration data and meteorological data on the monitoring boundary; reconstructing the light path integral concentration data, and obtaining a self-checking result according to the reconstruction result; and determining the position of a leakage source based on an optimization algorithm inversion according to the self-checking result. The application also discloses a petrochemical park gas leakage monitoring and positioning system, an electronic device and a storage medium. The application cooperates the monitoring boundary design and the active optical monitoring with the meteorological data, performs self-checking based on the light path integral concentration reconstruction, and then determines the leakage source through the optimization algorithm inversion. The application has strong implementability, can quickly and accurately position the leakage source under the premise of not affecting the normal production of the petrochemical park, and provides reliable guidance for subsequent emergency disposal.

Description

Technical Field

[0001] This invention relates to the field of atmospheric pollution source monitoring and location technology, and in particular to a method, system, electronic equipment and storage medium for monitoring and locating gas leaks in petrochemical industrial parks. Background Technology

[0002] Petrochemical enterprises are a pillar of the national economy, but they are also one of the most significant sources of air pollutant leaks. With the development of the national economy, the construction of petrochemical industrial parks is accelerating, and production activities are becoming increasingly frequent. Simultaneously, the frequency of air pollutant leaks in petrochemical industrial parks is also greatly increasing. Therefore, timely detection of gas leaks and their sources in the early stages of an accident is crucial for emergency response to pollution incidents and for minimizing losses.

[0003] In the area of ​​locating atmospheric pollutant leakage sources, current methods mainly include active olfaction based on mobile robots equipped with gas sensors and static gas source location methods relying on wireless sensor networks and atmospheric diffusion models. For example, patent document CN 109540141 A discloses a pollution source location mobile robot and method, which uses a pollution source location strategy combining a zigzag search algorithm and an evolutionary concentration gradient search algorithm through a pollution source location mobile robot equipped with a pollutant sensor array; patent document CN 104597212 A discloses an atmospheric pollution source location method that determines the location of pollution sources by deploying sensors at multiple monitoring points in a city and combining them with a Gaussian smoke cloud diffusion model.

[0004] However, due to the numerous equipment installations, complex terrain, and large area covered by petrochemical parks, the implementation of active olfactory methods for locating air pollutant leaks within these parks remains challenging due to limitations in the mobile robot's speed, obstacle avoidance, endurance, and explosion-proof capabilities. Furthermore, static gas source location methods relying on wireless sensor networks, due to the unique characteristics of petrochemical parks, cannot achieve full coverage of the entire park; they can only monitor certain key equipment online, and the accuracy of the location is closely related to the distribution of sensor points and the number of sensors used.

[0005] The information disclosed in this background section is intended only to enhance the understanding of the overall background of the invention and should not be construed as an admission or in any way implying that the information constitutes prior art known to those skilled in the art. Summary of the Invention

[0006] One of the objectives of this invention is to provide a method, system, electronic device, and storage medium for monitoring and locating gas leaks in petrochemical industrial parks, thereby improving the problem of accurately determining the location of gas leak sources in petrochemical industrial parks.

[0007] Another objective of this invention is to provide a method, system, electronic device, and storage medium for monitoring and locating gas leaks in petrochemical industrial parks, thereby improving the problem of high operational difficulty in existing monitoring technologies.

[0008] Another objective of this invention is to provide a method, system, electronic device, and storage medium for monitoring and locating gas leaks in petrochemical industrial parks, thereby improving the problem that existing monitoring technologies cannot comprehensively monitor target industrial parks.

[0009] To achieve the above objectives, according to a first aspect of the present invention, the present invention provides a method for monitoring and locating gas leaks in a petrochemical industrial park, comprising: determining the downwind boundary of the prevailing wind direction of the target park as the monitoring boundary; acquiring optical path integral concentration data and meteorological data on the monitoring boundary; reconstructing the optical path integral concentration data and obtaining a self-test result based on the reconstruction result; and determining the location of the leak source based on the self-test result and an optimization algorithm.

[0010] Furthermore, the above technical solution for the gas leak monitoring and location method in petrochemical parks also includes a step to verify the accuracy of the location: drawing a circle with the determined location of the leak source as the center and radius r, where r is 1% to 3% of the vertical distance between the location of the leak source and the monitoring boundary; drawing a straight line along the prevailing wind direction with the peak point of the reconstructed optical path integral concentration as the starting point; and determining whether the location of the leak source is accurate based on whether the straight line intersects the circle. If the straight line intersects the circle, the location is accurate; otherwise, the location is incorrect, and the process returns to the step of obtaining the optical path integral concentration data and meteorological data on the monitoring boundary.

[0011] Furthermore, in the above technical solution, the meteorological data includes wind direction (Wd), wind speed, and cloud cover.

[0012] Furthermore, in the above technical solution, an active optical monitoring device and an i-faceted pyramidal reflector are used to obtain optical path integral concentration data. The active optical monitoring device is deployed at the endpoint of the monitoring boundary, and the i-faceted pyramidal reflector is evenly distributed on the monitoring boundary and the other endpoint.

[0013] Furthermore, in the above technical solution, the step of reconstructing the optical path integral concentration data includes:

[0014] Construct the objective function S,

[0015]

[0016] In the formula, PIC i Let be the optical path integral density between the active optical monitoring device and the i-th pyramidal mirror, xi be the distance between the i-th pyramidal mirror and the active optical monitoring device, and B, m, and s be the fitting parameters; and

[0017] Minimize the objective function S to obtain the values ​​of B, m, and s. Then, the reconstructed optical path integral concentration is:

[0018]

[0019] Furthermore, in the above technical solution, the step of obtaining the self-test result based on the reconstruction result is as follows:

[0020] The self-test index FC is obtained based on the reconstructed optical path integral concentration.

[0021]

[0022] In the formula σ PIC PIC (Optical Path Integration Concentration) i The calculated standard deviation, σ c_pic The reconstructed optical path integral concentration c_pic i The calculated standard deviation, where R is the Pearson correlation coefficient. PIC (Optical Path Integration Concentration) i The calculated average value, The reconstructed optical path integral concentration c_pic i The calculated average value;

[0023] If FC > 0.8, the location of the leakage source is determined by inversion based on the optimization algorithm;

[0024] If FC≤0.8, then return to the steps of obtaining the optical path integral concentration data and meteorological data on the monitoring boundary.

[0025] Furthermore, in the above technical solution, the location of the leakage source is determined by inversion based on the optimization algorithm using an elevated continuous point source atmospheric diffusion Gaussian model.

[0026] Furthermore, in the above technical solution, determining the location of the leakage source based on the optimization algorithm includes coordinate transformation, where the coordinates are transformed into:

[0027] With the location of the active optical monitoring device as the origin P, and east as the x-axis and north as the y-axis, a rectangular coordinate system NPE is constructed. The coordinates of the i-th pyramidal reflector are (E... i N i ), in (E i-1 N i-1 ) and (E i N i Insert p-1 points evenly between the nth and nth points, where the coordinates of the nth point are... denoted as (E) k N k );

[0028] In the NPE coordinate system, the coordinates of the leak source are (E...L N L If the source of the leak is taken as the origin and the prevailing wind direction is taken as the X-axis, a rectangular coordinate system XOY is established according to the right-hand rule. Then, for any point (E... k N k The coordinates of ) in the XOY coordinate system are (X k Y k ),in

[0029] X k =(E k -E L )cosθ+(N k -N L sinθ

[0030] Y k =(N k -N L )cosθ-(E k -E L sinθ

[0031] θ = 90 - Wd.

[0032] Furthermore, in the above technical solution, in the XOY coordinate system, according to the elevated continuous point source atmospheric diffusion Gaussian model, the coordinates (X... k Y k The concentration at () is

[0033]

[0034] In the formula, Q represents the emission rate from the leakage source, u represents the average wind speed, and σ represents the average wind speed. y σ represents the horizontal diffusion parameter of the leakage source. z The vertical diffusion parameter of the leakage source is represented, z is the height of the active optical monitoring device, and H is the relative height between the leakage source and the active optical monitoring device.

[0035] Construct the objective function SSE.

[0036]

[0037] In the formula E L N L Q and H are the fitting parameters;

[0038] Minimize the objective function SSE to obtain the fitting parameters E. L N L The values ​​of Q and H determine the coordinates of the leak source as (E). L N L ,z+H).

[0039] According to a second aspect of the present invention, a gas leak monitoring and location system for a petrochemical industrial park is provided, comprising: a data acquisition unit for acquiring optical path integral concentration data and meteorological data in real time at a monitoring boundary, wherein the monitoring boundary is the downwind boundary of the prevailing wind direction of the target industrial park; a self-testing unit for reconstructing the optical path integral concentration data and obtaining a self-testing result based on the reconstructed optical path integral concentration; an inversion unit for determining the location of the leak source based on an optimization algorithm; and a ground control unit communicatively connected to the data acquisition unit, the self-testing unit, and the inversion unit.

[0040] Furthermore, in the above technical solution, the petrochemical park gas leak monitoring and positioning system also includes: a positioning accuracy verification unit, which is communicatively connected to the ground control module. The positioning accuracy verification unit draws a circle with the determined location of the leak source as the center and r as the radius, where r is 1% to 3% of the vertical distance between the location of the leak source and the monitoring boundary; it draws a straight line along the prevailing wind direction with the peak point of the optical path integral concentration reconstructed by the self-test unit as the starting point; and it determines whether the location of the leak source is accurately positioned based on whether the straight line intersects with the circle.

[0041] Furthermore, in the above technical solution, the data acquisition unit includes: an active optical monitoring device, which is deployed at the endpoint of the monitoring boundary; an i-faceted pyramidal reflector, which is evenly distributed on the monitoring boundary and its other endpoint; and a meteorological data acquisition module, which is set in the middle of the monitoring boundary and is unobstructed around it.

[0042] Furthermore, in the above technical solution, the active optical monitoring device can be adjusted horizontally by 180° and tilted by 45°.

[0043] Furthermore, in the above technical solution, the active optical monitoring equipment uses a pan-tilt unit to adjust the horizontal and vertical angles.

[0044] Furthermore, in the above technical solution, the ground control unit controls the operation mode and dwell time of the gimbal; the ground control unit controls the data acquisition cycle of the data acquisition unit.

[0045] According to a third aspect of the present invention, an electronic device is provided, comprising: at least one processor; and a memory communicatively connected to the at least one processor; wherein the memory stores instructions executable by the at least one processor, the instructions being executed by the at least one processor to cause the at least one processor to perform a petrochemical park gas leak monitoring and location method as described in any of the above technical solutions.

[0046] According to a fourth aspect of the present invention, the present invention provides a non-transitory computer-readable storage medium storing computer-executable instructions for causing a computer to perform a gas leak monitoring and location method in a petrochemical park as described in any of the above technical solutions.

[0047] Compared with the prior art, the present invention has one or more of the following beneficial effects:

[0048] 1. This invention combines monitoring boundary design and active optical monitoring with meteorological data, performs self-checking based on optical path integral concentration reconstruction, and then uses an optimization algorithm to determine the leak source. It is highly feasible and can quickly and accurately locate the leak source without affecting normal production in the petrochemical park, providing reliable guidance for subsequent emergency response.

[0049] 2. Compared with the active olfactory method based on mobile robots, this invention avoids the problems of difficult robot path planning and high practical operation difficulty caused by the dense buildings and complex terrain in the petrochemical park. Compared with the static gas source positioning method based on wireless sensor network, the active optical remote sensing monitoring equipment can monitor the boundary and achieve full coverage of the target park. Moreover, it is easy to install on site and gets rid of the geographical limitations of wireless sensor network.

[0050] 3. This invention can verify the accuracy of leak source location, ensuring the accuracy of leak source location.

[0051] 4. This invention employs an active optical monitoring device, which offers fast monitoring speed and high sensitivity, and can simultaneously monitor multiple components. It can achieve rapid and accurate qualitative and quantitative analysis of the target gas within seconds, enabling timely and early detection of gas leaks and further improving the speed of locating gas leak sources. In addition, the concentration data collected by this invention is optical path integrated concentration data, covering the concentration information along the entire monitoring optical path. This avoids the "outlier" problem caused by factors such as wind field and terrain in point monitoring, resulting in higher stability and accuracy.

[0052] The above description is merely an overview of the technical solution of the present invention. In order to better understand the technical means of the present invention and to implement it according to the contents of the specification, and to make the above and other objects, technical features and advantages of the present invention easier to understand, one or more preferred embodiments are listed below and described in detail with reference to the accompanying drawings. Attached Figure Description

[0053] Figure 1 This is a flowchart of a gas leak monitoring and location method in a petrochemical industrial park according to an embodiment of the present invention.

[0054] Figure 2This is a schematic diagram of a gas leak monitoring and location system in a petrochemical industrial park according to an embodiment of the present invention.

[0055] Figure 3 This is a schematic diagram of the data acquisition unit according to Embodiment 3 of the present invention.

[0056] Figure 4 This is a schematic diagram of the data acquisition unit according to Embodiment 4 of the present invention.

[0057] Figure 5 This is a schematic diagram of the data acquisition unit according to Embodiment 5 of the present invention.

[0058] Figure 6 This is a schematic diagram of the hardware structure of an electronic device for implementing a gas leak monitoring and location method in a petrochemical industrial park according to an embodiment of the present invention. Detailed Implementation

[0059] The specific embodiments of the present invention will now be described in detail with reference to the accompanying drawings, but it should be understood that the scope of protection of the present invention is not limited to the specific embodiments.

[0060] Unless otherwise expressly stated, throughout the specification and claims, the term "comprising" or its variations such as "including" or "comprises" shall be understood to include the stated elements or components without excluding other elements or other components.

[0061] In this document, for ease of description, spatial relative terms such as “below,” “under,” “down,” “above,” “above,” “up,” etc., are used to describe the relationship of one element or feature to another element or feature in the accompanying drawings. It should be understood that spatial relative terms are intended to encompass different orientations of an object in use or operation, in addition to those depicted in the figures. For example, if an object in the figure is flipped, an element described as “below” or “under” another element or feature would be oriented “above” that element or feature. Thus, the exemplary term “below” can encompass both the downward and upward orientations. An object may also have other orientations (rotated 90 degrees or other orientations), and the spatial relative terms used herein should be interpreted accordingly.

[0062] In this document, the terms "first," "second," etc., are used to distinguish two different elements or parts, and are not used to define specific positions or relative relationships. In other words, in some embodiments, the terms "first," "second," etc., can also be used interchangeably.

[0063] The method, system, electronic equipment and storage medium for monitoring and locating gas leaks in petrochemical industrial parks of the present invention are described in more detail below by way of specific embodiments. It should be understood that the embodiments are merely exemplary and the present invention is not limited thereto.

[0064] Example 1

[0065] Combination Figure 1 As shown, the flow of the gas leak monitoring and location method in the petrochemical park according to this embodiment is as follows:

[0066] S110 defines the monitoring boundary.

[0067] First, based on the target area's layout plan and the prevailing wind direction, conduct a patrol of the downwind boundary to ensure sufficient space for the deployment of the boundary monitoring system. The downwind boundary of the target area's prevailing wind direction is then designated as the monitoring boundary. If there is no prevailing wind direction, the entire boundary of the target area must be designated as the monitoring boundary.

[0068] S120 acquires optical path integral concentration data and meteorological data at the monitoring boundary.

[0069] Active optical monitoring equipment and i-faced pyramidal reflectors are used to acquire integrated optical path concentration data. The active optical monitoring equipment is deployed at one end of the monitoring boundary, and the i-faced pyramidal reflectors are evenly distributed at the monitoring boundary and the other end. Meteorological data may include wind direction (Wd), wind speed, and cloud cover.

[0070] The S130 reconstructs the optical path integral concentration data and obtains the self-test results based on the reconstruction results.

[0071] S131 Construct the objective function S,

[0072]

[0073] In the formula, PIC i Let be the optical path integral density between the active optical monitoring device and the i-th pyramidal mirror, xi be the distance between the i-th pyramidal mirror and the active optical monitoring device, and B, m, and s be the fitting parameters; and

[0074] Minimize the objective function S to obtain the values ​​of B, m, and s. Then, the reconstructed optical path integral concentration is:

[0075]

[0076] S132 obtains the self-test index FC based on the reconstructed optical path integral concentration.

[0077]

[0078] In the formula σ PIC PIC (Optical Path Integration Concentration) i The calculated standard deviation, σ c_pic The reconstructed optical path integral concentration c_pic iThe calculated standard deviation, where R is the Pearson correlation coefficient. PIC (Optical Path Integration Concentration) i The calculated average value, The reconstructed optical path integral concentration c_pic i The calculated average value.

[0079] S133 determines whether to continue to the next step based on the self-test index FC.

[0080] If FC > 0.8, the location of the leakage source is determined by inversion based on the optimization algorithm;

[0081] If FC≤0.8, then return to step S120.

[0082] Based on the self-test results, S140 determines the location of the leak source using an optimization algorithm.

[0083] S141 coordinate transformation:

[0084] With the location of the active optical monitoring device as the origin P, and east as the x-axis and north as the y-axis, a rectangular coordinate system NPE is constructed. The coordinates of the i-th pyramidal reflector are (E... i N i ), in (E i-1 N i-1 ) and (E i N i Insert p-1 points evenly between the nth and nth points, where the coordinates of the nth point are... denoted as (E) k N k );

[0085] In the NPE coordinate system, the coordinates of the leak source are (E... L N L If the source of the leak is taken as the origin and the prevailing wind direction is taken as the X-axis, a rectangular coordinate system XOY is established according to the right-hand rule. Then, for any point (E... k N k The coordinates of ) in the XOY coordinate system are (X k Y k ),in

[0086] X k =(E k -E L )cosθ+(N k -N L sinθ

[0087] Y k =(N k -N L )cosθ-(E k-E L sinθ

[0088] θ = 90 - Wd.

[0089] S142 used an elevated continuous point source atmospheric diffusion Gaussian model to determine the location of the leak source.

[0090] Based on the Gaussian model of atmospheric diffusion from an elevated continuous point source, the coordinates (X...) k Y k The concentration at () is

[0091]

[0092] In the formula, Q represents the emission rate of the leakage source (an unknown quantity to be fitted), u represents the average wind speed (the average wind speed within one cycle of cyclic monitoring by the active optical monitoring equipment using the i-face pyramidal reflector), and σ represents the average wind speed. y σ represents the horizontal diffusion parameter of the leakage source. z This represents the vertical diffusion parameter of the leak source, which is related to cloud cover. z is the height of the active optical monitoring device, and H is the relative height between the leak source and the active optical monitoring device.

[0093] Construct the objective function SSE.

[0094]

[0095] In the formula E L N L Q and H are the fitting parameters;

[0096] Minimize the objective function SSE to obtain the fitting parameters E. L N L The values ​​of Q and H determine the coordinates of the leak source as (E). L N L ,z+H).

[0097] S150 checks the positioning accuracy.

[0098] Draw a circle with the determined location of the leak source as the center and radius r, where r is 1% to 3% of the vertical distance between the location of the leak source and the monitoring boundary; draw a straight line along the prevailing wind direction starting from the point where the peak of the reconstructed optical path integral concentration is located; and determine whether the location of the leak source is accurate based on whether the straight line intersects the circle. If the straight line intersects the circle, the location is accurate, and the location of the leak source is output; otherwise, the location is incorrect, and return to step S120.

[0099] Specifically, the reconstructed optical path integral concentration c_pic i The concentration at any point on the monitoring boundary is

[0100] In the formula, x is the distance between any point on the monitoring boundary and the active optical monitoring device, then x = m is the point where the peak value of the reconstructed optical path integral concentration is located. In the NPE coordinate system, a straight line is drawn along the prevailing wind direction, starting from the point where the peak value of the reconstructed optical path integral concentration is located; with the leakage source (E L N L Draw a circle with center r and radius r. If the straight line passes through the circular area centered at the leak source point, the leak source is accurately located; otherwise, the location is incorrect.

[0101] Example 2

[0102] Combination Figure 2 As shown, the petrochemical park gas leak monitoring and location system of this embodiment includes: a data acquisition unit 10, which is used to acquire optical path integral concentration data and meteorological data of the monitoring boundary in real time, the monitoring boundary being the downwind boundary of the target park's prevailing wind direction; a self-testing unit 20, which reconstructs the optical path integral concentration data and obtains a self-test result based on the reconstructed optical path integral concentration; an inversion unit 30, which determines the location of the leak source based on an optimization algorithm; and a location accuracy verification unit 40, which is communicatively connected to the ground control module. The location accuracy verification unit draws a circle with the determined leak source location as the center and a radius of r, where r is 1% to 3% of the vertical distance between the leak source location and the monitoring boundary; draws a straight line along the prevailing wind direction starting from the peak point of the optical path integral concentration reconstructed by the self-testing unit; and determines whether the location of the leak source is accurate based on whether the straight line intersects the circle. The petrochemical park gas leak monitoring and location system of this embodiment is controlled by a ground control unit 50, which is communicatively connected to the data acquisition unit 10, the self-testing unit 20, the inversion unit 30, and the location accuracy verification unit 40. The data acquisition unit 10 may include: an active optical monitoring device, which is deployed at the endpoint of the monitoring boundary; an i-faced pyramidal reflector, which is evenly distributed on the monitoring boundary and its other endpoint; and a meteorological data acquisition module, which is set in the middle of the monitoring boundary and is unobstructed.

[0103] The ground control unit 50 can perform the following functions via a wireless transmission system: setting the operation mode of the pan-tilt unit, such as memorizing the positions of the multi-faceted pyramidal reflectors to enable sequential cyclic monitoring of the multi-faceted pyramidal reflectors; setting the monitoring order of the multi-faceted pyramidal reflectors; setting the dwell time on each pyramidal reflector; and setting the data acquisition cycle of the data acquisition unit 10. Generally, the dwell time can be set to be consistent with the data acquisition cycle. In addition, the ground control unit 50 can also acquire distance data between the pyramidal reflectors and the active optical monitoring equipment, as well as the height data of the active optical monitoring equipment above the ground.

[0104] For specific monitoring and positioning methods, please refer to Example 1, which will not be repeated here.

[0105] Example 3

[0106] like Figure 3 As shown, this embodiment applies the present invention to the gas leak monitoring and location of a target petrochemical industrial park 200. The target petrochemical industrial park 200 experiences significant wind direction fluctuations throughout the year, with no clear dominant wind direction; therefore, the monitoring boundary is the entire boundary of the petrochemical industrial park 200. The data acquisition unit of this embodiment is equipped with two active optical monitoring devices 310 and a multifaceted pyramidal reflector 320, thereby achieving full coverage monitoring of the entire target petrochemical industrial park 200.

[0107] The two active optical monitoring devices 310 of the data acquisition unit are equipped with pan-tilt units (not shown in the figure) and are positioned at the endpoints of the approximate diagonal of the boundary of the target petrochemical park 200. The target petrochemical park 200 is square with a side length of approximately 2.8 km. The leaked gas is mainly volatile organic compounds. Therefore, the active optical monitoring device 310 uses an active Fourier transform infrared spectrometer (OP-FTIR). A pyramidal reflector 320 is placed at 0.7 km intervals along each boundary, with 4 pyramidal reflectors 320 on each boundary (the pyramidal reflectors at the endpoints of the two boundaries are shared), for a total of 14 pyramidal reflectors 320 (only 6 pyramidal reflectors are shown as an example in the figure). The pyramidal reflectors 320 are at the same height as the OP-FTIR. To avoid the industrial park boundary wall from obstructing the leaked and diffused gas, the height z of the OP-FTIR is 2 m. The four pyramidal mirrors 320 on the same boundary should be approximately aligned with the OP-FTIR, but to avoid mutual obstruction affecting data acquisition, the positions of the four pyramidal mirrors can be adjusted appropriately. The adjustment should be based on the principle that the pan / tilt unit can control the OP-FTIR main unit and the pyramidal mirrors to complete focusing, while minimizing the adjustment offset of the four pyramidal mirrors. The pan / tilt unit can control the OP-FTIR main unit to perform 180° horizontal adjustment and 45° pitch adjustment, thereby sequentially focusing the four pyramidal mirrors from near to far, acquiring the integrated optical density data (PIC) for each of the four monitoring optical paths. 1-4 After the integrated concentration data of the four monitoring optical paths are fully collected, the pan-tilt unit can automatically adjust its angle to focus on the four-sided pyramidal reflector on the other boundary, and similarly acquire the integrated concentration data of the four monitoring optical paths. Ultimately, each active optical monitoring device 310 achieves periodic cyclic monitoring of the eight-sided pyramidal reflector. The meteorological data acquisition module 330 is located in the middle of the monitoring boundary, at a height of approximately 10m, with no obvious obstacles nearby. The meteorological data acquisition module 330 collects data such as wind direction, wind speed, and cloud cover.

[0108] In this embodiment, the dwell time and data acquisition cycle of the active optical monitoring device 310 are both set to 1 minute. The specific monitoring and positioning method can be found in Embodiment 1, and will not be repeated here.

[0109] Example 4

[0110] Combination Figure 4 As shown, this embodiment applies the present invention to gas leak monitoring and location of a target petrochemical industrial park 400. The prevailing wind direction 401 of the target petrochemical industrial park 400 is obvious; therefore, the monitoring boundary is the downwind boundary 402 of the prevailing wind direction 401. The data acquisition unit of this embodiment includes an active optical monitoring device 510 and a bifacial pyramidal reflector 520, thereby achieving full coverage monitoring of the entire target petrochemical industrial park 400. The meteorological data acquisition module 530 is located in the middle of the monitoring boundary, at a height of approximately 10m, with no obvious obstacles nearby. The meteorological data acquisition module 530 collects data such as wind direction, wind speed, and cloud cover. Specific monitoring and location methods can be found in Embodiment 1, and will not be repeated here.

[0111] Example 5

[0112] Combination Figure 5 As shown, this embodiment applies the present invention to gas leak monitoring and location at a target petrochemical park 700. Although the prevailing wind direction 701 at the target petrochemical park 700 is stable year-round, and the prevailing wind direction 701 is northerly, the downwind boundary 702 (i.e., the monitoring distance) is too long, approximately 4.8 km. Due to the limitation of the optical path length that the active optical monitoring equipment can detect, it is difficult to fully cover the downwind boundary in one direction. Therefore, the active optical monitoring equipment 810 is fixed at the middle position of the monitoring boundary, and eight pyramidal reflectors 820 are evenly distributed on the east and west sides at a spacing of 0.6 km. The meteorological data acquisition module 830 is set near the active optical monitoring equipment 810 at a height of approximately 10 m, with no obvious obstacles nearby. The meteorological data acquisition module 830 collects data such as wind direction, wind speed, and cloud cover. The specific monitoring and location method can be referred to in Embodiment 1, and will not be repeated here.

[0113] Example 6

[0114] This embodiment provides a non-transitory (non-volatile) computer storage medium that stores computer-executable instructions that can execute the methods in any of the above method embodiments and achieve the same technical effect.

[0115] Example 7

[0116] This embodiment provides a computer program product, which includes a computer program stored on a non-transitory computer-readable storage medium. The computer program includes program instructions, which, when executed by a computer, cause the computer to perform the methods described above and achieve the same technical effects.

[0117] Example 8

[0118] Figure 6 This is a schematic diagram of the hardware structure of the electronic device for performing the lightning strike pipeline simulation evaluation method according to this embodiment. The device includes one or more processors 610 and a memory 620. Taking one processor 610 as an example, the device may also include an input device 630 and an output device 640.

[0119] The processor 610, memory 620, input device 630, and output device 640 can be connected via a bus or other means. Figure 3 Taking the example of a connection between China and Israel via a bus.

[0120] The memory 620, as a non-transitory computer-readable storage medium, can be used to store non-transitory software programs, non-transitory computer-executable programs, and modules. The processor 610 executes various functional applications and data processing of the electronic device by running the non-transitory software programs, instructions, and modules stored in the memory 620, thereby implementing the processing method of the above-described method embodiments.

[0121] The memory 620 may include a program storage area and a data storage area, wherein the program storage area may store the operating system and applications required for at least one function; the data storage area may store data, etc. Furthermore, the memory 620 may include high-speed random access memory, and may also include non-transitory memory, such as at least one disk storage device, flash memory device, or other non-transitory solid-state storage device. In some embodiments, the memory 620 may optionally include memory remotely located relative to the processor 610, and these remote memories may be connected to the processing device via a network. Examples of such networks include, but are not limited to, the Internet, corporate intranets, local area networks, mobile communication networks, and combinations thereof.

[0122] Input device 630 can receive input digital or character information and generate signal input. Output device 640 may include display devices such as a display screen.

[0123] One or more modules are stored in memory 620 and, when executed by one or more processors 610, execute:

[0124] The downwind boundary of the target park's prevailing wind direction is defined as the monitoring boundary, and optical path integral concentration data and meteorological data on the monitoring boundary are obtained;

[0125] Reconstruct the optical path integral concentration data and obtain self-test results based on the reconstruction results; and

[0126] Based on the self-inspection results, the location of the leak source is determined by inversion using an optimization algorithm.

[0127] The above-described product can execute the methods provided in the embodiments of the present invention, and has the corresponding functional modules and beneficial effects for executing the methods. Technical details not described in detail in this embodiment can be found in the methods provided in other embodiments of the present invention.

[0128] The device embodiments described above are merely illustrative. The units described as separate components may or may not be physically separate, and the components shown as units may or may not be physical units; that is, they may be located in one place or distributed across multiple network units. Some or all of the modules can be selected to achieve the purpose of this embodiment according to actual needs.

[0129] Through the above description of the embodiments, those skilled in the art can clearly understand that each embodiment can be implemented by means of software plus a general-purpose hardware platform, and of course, it can also be implemented by hardware. Based on this understanding, the above technical solutions, in essence or the parts that contribute to the related technology, can be embodied in the form of a software product. This computer software product can be stored in a computer-readable storage medium, such as ROM / RAM, magnetic disk, optical disk, etc., including several instructions to cause a computer device (which may be a personal computer, server, or network device, etc.) to execute the methods of various embodiments or some parts of embodiments.

[0130] The foregoing description of specific exemplary embodiments of the present invention is for illustrative and explanatory purposes. These descriptions are not intended to limit the invention to the precise forms disclosed, and it will be apparent that many changes and variations can be made in accordance with the foregoing teachings. The exemplary embodiments were chosen and described in order to explain the specific principles of the invention and its practical application, thereby enabling those skilled in the art to implement and utilize various different exemplary embodiments of the invention, as well as various different choices and variations. Any simple modifications, equivalent changes, and alterations made to the foregoing exemplary embodiments should fall within the scope of protection of the present invention.

Claims

1. A method for monitoring and locating gas leaks in a petrochemical industrial park, characterized in that, include: The downwind boundary of the prevailing wind direction of the target park is defined as the monitoring boundary, and optical path integral concentration data and meteorological data on the monitoring boundary are acquired; the meteorological data includes wind direction Wd, wind speed and cloud cover. The optical path integral concentration data is reconstructed, and a self-test result is obtained based on the reconstruction result; The steps for reconstructing the optical path integral concentration data include: Construct the objective function S, ， In the formula, PIC i Let be the optical path integral concentration between the active optical monitoring device and the i-th pyramidal mirror, xi be the distance between the i-th pyramidal mirror and the active optical monitoring device, and B, m, and s be the fitting parameters; Minimize the objective function S to obtain the values ​​of B, m, and s; The reconstructed optical path integral concentration is: ； The steps to obtain the self-test result based on the reconstruction result are as follows: obtain the self-test index FC based on the reconstructed optical path integral concentration; if FC > 0.8, determine the location of the leakage source based on the optimization algorithm; if FC ≤ 0.8, return to the steps of obtaining the optical path integral concentration data and meteorological data on the monitoring boundary. Based on the self-test results, the location of the leak source is determined by inversion using an optimization algorithm; determining the location of the leak source by inversion using an optimization algorithm includes coordinate transformation to obtain the coordinates of the specific leak point.

2. The method for monitoring and locating gas leaks in petrochemical industrial parks according to claim 1, characterized in that, It also includes steps to verify the accuracy of the positioning: Draw a circle with the location of the determined leak source as the center and radius r, where r is 1% to 3% of the vertical distance between the location of the leak source and the monitoring boundary; Draw a straight line along the prevailing wind direction, starting from the point where the peak value of the reconstructed optical path integral concentration is located. as well as The location of the leak source can be determined accurately based on whether the straight line intersects the circle. If the straight line intersects the circle, the positioning is accurate; Otherwise, if the location is incorrect, return to the step of obtaining the optical path integral concentration data and meteorological data on the monitoring boundary.

3. The method for monitoring and locating gas leaks in petrochemical industrial parks according to claim 2, characterized in that, The optical path integral concentration data is obtained using an active optical monitoring device and an i-faceted pyramidal reflector. The active optical monitoring device is deployed at one end of the monitoring boundary, and the i-faceted pyramidal reflector is evenly distributed on the monitoring boundary and the other end.

4. The method for monitoring and locating gas leaks in petrochemical industrial parks according to claim 1, characterized in that, The self-test index ， In the formula PIC (Optical Path Integration Concentration) i The calculated standard deviation The integrated concentration of the reconstructed optical path The calculated standard deviation, where R is the Pearson correlation coefficient. PIC (Optical Path Integration Concentration) i The calculated average value, The integrated concentration of the reconstructed optical path The calculated average value.

5. The method for monitoring and locating gas leaks in petrochemical industrial parks according to claim 4, characterized in that, The location of the leakage source was determined by inversion using an optimization algorithm, employing an elevated continuous point source atmospheric diffusion Gaussian model.

6. The method for monitoring and locating gas leaks in petrochemical industrial parks according to claim 5, characterized in that, The coordinates are converted to: With the location of the active optical monitoring device as the origin P, and east as the x-axis and north as the y-axis, a rectangular coordinate system NPE is constructed. The coordinates of the i-th facet pyramidal reflector are (E... i N i ), in (E i-1 N i-1 ) and (E i N i Insert p-1 points evenly between the nth and nth points, where the coordinates of the nth point are... (E i-1+ (E i -E i-1 ), N i-1 + (Ni-N i-1 )), record (E k , N k ); In the NPE coordinate system, the coordinates of the leak source are (E... L N L ), then with the leak source as the origin and the prevailing wind direction as the X-axis, a rectangular coordinate system XOY is established according to the right-hand rule. Then any point (E) k N k The coordinates of in the XOY coordinate system are (X k Y k ),in ， ， 。 7. The method for monitoring and locating gas leaks in petrochemical industrial parks according to claim 6, characterized in that, In the XOY coordinate system, according to the Gaussian model of atmospheric diffusion from an elevated continuous point source, the coordinates (X...) k Y k The concentration at () is ， In the formula, Q represents the emission rate from the leakage source, and u represents the average wind speed. This represents the horizontal diffusion parameters of the leak source. The vertical diffusion parameter of the leakage source is represented, z is the height of the active optical monitoring device, and H is the relative height between the leakage source and the active optical monitoring device. Construct the objective function SSE. ， In the formula E L N L Q and H are the fitting parameters; Minimize the objective function SSE to obtain the fitting parameters E. L N L The values ​​of Q and H determine the coordinates of the leak source as (E). L N L ,z+H).

8. A gas leak monitoring and location system for petrochemical industrial parks, characterized in that, The method described in any one of claims 1 to 7 includes: The data acquisition unit is used to acquire optical path integral concentration data and meteorological data of the monitoring boundary in real time. The monitoring boundary is the downwind boundary of the prevailing wind direction of the target park. The self-test unit reconstructs the optical path integral concentration data and derives the self-test result based on the reconstructed optical path integral concentration. The inversion unit determines the location of the leakage source based on an optimization algorithm; and The ground control unit is communicatively connected to the data acquisition unit, the self-test unit, and the inversion unit.

9. The petrochemical industrial park gas leak monitoring and location system according to claim 8, characterized in that, Also includes: The positioning accuracy verification unit is communicatively connected to the ground control unit. The positioning accuracy verification unit draws a circle with the determined location of the leak source as the center and a radius of r, where r is 1% to 3% of the vertical distance between the location of the leak source and the monitoring boundary; it draws a straight line along the prevailing wind direction with the peak point of the optical path integral concentration reconstructed by the self-test unit as the starting point; and it determines whether the location of the leak source is accurate based on whether the straight line intersects with the circle.

10. The petrochemical industrial park gas leak monitoring and location system according to claim 8, characterized in that, The data acquisition unit includes: An active optical monitoring device is deployed at the endpoints of the monitoring boundary; i-faceted pyramidal reflectors, uniformly distributed along the monitoring boundary and its other endpoint; and The meteorological data acquisition module is located in the middle of the monitoring boundary and is unobstructed in the surrounding area.

11. The petrochemical industrial park gas leak monitoring and location system according to claim 10, characterized in that, The active optical monitoring device can be adjusted horizontally by 180° and tilted by 45°.

12. The petrochemical industrial park gas leak monitoring and location system according to claim 11, characterized in that, The active optical monitoring device uses a pan-tilt unit to adjust the horizontal and vertical angles.

13. The petrochemical industrial park gas leak monitoring and location system according to claim 12, characterized in that, The ground control unit controls the operation mode and dwell time of the gimbal; the ground control unit controls the data acquisition cycle of the data acquisition unit.

14. An electronic device, characterized in that, include: At least one processor; as well as A memory communicatively connected to the at least one processor; wherein, The memory stores instructions that can be executed by the at least one processor, which, when executed by the at least one processor, causes the at least one processor to perform the petrochemical park gas leak monitoring and location method as described in any one of claims 1 to 7.

15. A non-transitory computer-readable storage medium, characterized in that, The non-transitory computer-readable storage medium stores computer-executable instructions, which are used to cause the computer to perform the petrochemical park gas leak monitoring and location method as described in any one of claims 1 to 7.

Images