Monitoring point distribution method and equipment based on industrial park gas monitoring and medium

Abstract

The invention discloses a monitoring point distribution method and device based on industrial park gas monitoring and a medium, and relates to the technical field of carbon monitoring, triangular mesh division is carried out based on DT triangular mesh division, and according to the difference of monitored areas, the DT triangular mesh division method is used for setting point distribution of the monitored areas, so that the monitoring point distribution accuracy is improved. At least one monitoring point is selected in each grid for carbon dioxide concentration monitoring, point distribution setting can be carried out in different area ranges of prevailing wind directions in different states, the optimal monitoring result can be obtained through the least carbon dioxide gas monitoring points, the accuracy and credibility of the final data result are improved, and the method is suitable for popularization and application. The problems that the current industrial park monitoring coverage range is not wide, the device cost is high, and improvement of the carbon monitoring device distribution position affects the accuracy of the carbon ground monitoring technology are solved.

Description

Technical Field

[0001] The present invention relates to the technical field of carbon monitoring, and in particular to a monitoring point arrangement method, equipment and medium based on industrial park gas monitoring. Background Art

[0002] As an important carrier of my country's industrial development, industrial parks contribute more than 50% of the country's industrial output; but at the same time, industrial parks gather a large number of enterprises in a relatively small space, with large resource and energy consumption, concentrated pollutant emissions, and carbon emissions accounting for about 31% of national emissions. Accurate monitoring of greenhouse gas emissions in industrial parks is the basic premise for industrial parks to formulate carbon reduction paths and an important support for the implementation of the dual carbon strategy in the industrial field. Whether the layout of monitoring points is representative will greatly affect the reliability and accuracy of the monitoring results.

[0003] In the existing industrial park monitoring system, most of the research on carbon dioxide monitoring focuses on the monitoring equipment level, and focuses on improving the hardware perspective of carbon dioxide sensor monitoring accuracy. A few studies consider the deployment perspective, and at the same time, there is no consideration of measurement accuracy and result confidence.

[0004] my country's industrial parks are numerous, diverse, and at different stages of development. While continuing to promote green development of industrial parks and green manufacturing of products, it is also very important to systematically and scientifically monitor the main carbon dioxide emission points in industrial parks to provide data support for carbon emission reduction.

[0005] Due to the limitations of manpower, equipment, funds and other costs required by the existing monitoring methods, it is impossible to conduct high-density monitoring of carbon dioxide in the entire industrial park. Therefore, in-depth research is urgently needed to make the layout of monitoring points representative, scientific, reasonable and feasible. Summary of the invention

[0006] The technical problem to be solved by the present invention is that the existing monitoring method for carbon dioxide gas monitoring does not cover all aspects, has low monitoring accuracy, and cannot effectively monitor gas. The purpose is to provide a monitoring point distribution method, equipment and medium based on industrial park gas monitoring, which performs triangular mesh division based on DT triangular mesh division, and sets the points for the monitoring area according to the differences in the monitored areas through the DT triangular mesh division method. By selecting at least one monitoring point in each grid to monitor the carbon dioxide concentration, the points can be set in different areas with different dominant wind directions. The best monitoring results can be obtained with the least carbon dioxide gas monitoring points, and the accuracy and credibility of the final data results can be improved, so as to solve the problems of the current industrial park monitoring coverage being limited, the device cost being high, and the location of the carbon monitoring device affecting the accuracy of the carbon ground monitoring technology.

[0007] The present invention is achieved through the following technical solutions: The first aspect of the present invention provides a monitoring point arrangement method based on industrial park gas monitoring, comprising the following specific steps: Obtain park characteristic data of the target monitoring area, and construct a characteristic data point set based on the characteristic data; Based on the feature data point set, the DT triangulation method is used to triangulate the target monitoring area and generate a triangular mesh; Obtain the measurement values ​​of the sensors at the three nodes of each triangular mesh and calculate the coefficient of variation of each triangular mesh; In the ith triangle mesh: If the coefficient of variation is within the set threshold, the center point of the i-th triangular mesh is selected for point distribution; If the coefficient of variation is not within the set threshold, it is determined whether the circumscribed circle of the i-th triangular mesh is smaller than the sensor monitoring radius; If it is not less than, the Delaumnay triangulation algorithm is used to encrypt the i-th triangular mesh, and the execution is jumped to obtain the measurement values ​​of the sensors on the three nodes of each triangular mesh; If it is less than, the center point of the i-th triangle mesh is selected for point distribution.

[0008] The present invention performs triangular mesh division based on DT triangular mesh division, and sets the points for the monitoring area according to the differences in the monitored areas through the DT triangular mesh division method. By selecting at least one monitoring point in each grid to monitor the carbon dioxide concentration, the points can be set in different areas with dominant wind directions in different states. The best monitoring results can be obtained with the least number of carbon dioxide gas monitoring points, and the accuracy and credibility of the final data results can be improved. The problems of the current industrial park monitoring coverage being limited, the device cost being high, and the location of the carbon monitoring device affecting the accuracy of the carbon ground monitoring technology can be solved.

[0009] Furthermore, the park characteristic data of the target monitoring area includes: the characteristics of the park itself in the target monitoring area, the distribution of emission data of the target monitoring gas in the area, and the meteorological data in the area.

[0010] Furthermore, the park characteristics of the target monitoring area include: the park's geographical location, regional boundaries, area, traffic conditions, connection with surrounding infrastructure, park scale, industry type and production process; The distribution of emission data of target monitored gases in the area includes: the location and size of workshops and other major buildings, organized and unorganized emission outlets and their emission parameters, types and emission rates of emission sources, the closedness of the target monitoring area of ​​the park, and the distribution of buildings and terrain in the park that affect airflow movement; The meteorological data in the area include: wind speed, wind direction, temperature, humidity, precipitation, surface radiation, temperature profile and tropospheric altitude.

[0011] Furthermore, when monitoring points are arranged, the meteorological data in the area are measured and determined, including: Collect wind speed and direction, obtain an instant wind direction and wind speed value every t time, perform continuous measurement within T period, obtain n wind direction values ​​and n wind speed values, and measure and judge meteorological data based on n wind direction values ​​and n wind speed values; During the collection of wind speed and direction, if the deviation of n wind direction values ​​is greater than the set threshold, the monitoring point is moved and resampled.

[0012] Furthermore, the DT triangulation method is used to triangulate the target monitoring area, specifically including: Obtain the boundary of the target monitoring area, and select the boundary point and the carbon emission source point in the center of the target monitoring area according to the boundary of the target monitoring area; According to the carbon emission source points at the boundary points and the center of the target monitoring area, the target monitoring area is divided into multiple grids using the Delaunay function.

[0013] Furthermore, the calculation of the coefficient of variation of each triangular mesh specifically includes: , where σ is the standard deviation and μ is the mean.

[0014] Furthermore, before arranging points on the i-th triangular grid, the spatial variability structure of carbon dioxide characteristics is fitted using a semivariance function to obtain a function model. The specific steps include: ; Where: N(h) is the number of sample logarithms of the interval; γ(h) is the sample variance function of the regionalized variable Z(X); h is the spatial separation distance between two sample points; Z(X i ) and Z(X i +h) are the localized variables Z(X) at the spatial position X i and X i +h's observed value.

[0015] Furthermore, after the points are obtained based on the triangular mesh, the points are optimized according to the actual situation of the park, including: Obtain the monitored gas emission data of the target monitoring area, and determine the first dominant wind direction and the second dominant wind direction of the carbon emission data based on the monitored gas emission source data and meteorological data; Based on the first dominant wind direction and the second dominant wind direction of the emission data, the gas concentration distribution of each grid is obtained; The points are optimized according to the gas concentration distribution of each grid.

[0016] A second aspect of the present invention provides an electronic device, including a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein when the processor executes the program, a monitoring point deployment method based on industrial park gas monitoring is implemented.

[0017] A third aspect of the present invention provides a computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements a monitoring point deployment method based on industrial park gas monitoring.

[0018] Compared with the prior art, the present invention has the following advantages and beneficial effects: The present invention performs triangular mesh division based on DT triangular mesh division, and sets the points for the monitoring area according to the differences in the monitored areas through the DT triangular mesh division method. By selecting at least one monitoring point in each grid to monitor the carbon dioxide concentration, the points can be set in different areas with different dominant wind directions. The best monitoring results can be obtained with the least number of carbon dioxide gas monitoring points, and the accuracy and credibility of the final data results can be improved, solving the problems of the current industrial park monitoring coverage being limited, the device cost being high, and the location of the carbon monitoring device affecting the accuracy of the carbon ground monitoring technology. The present invention comprehensively considers the carbon emissions of industrial parks, meteorological conditions, functional zoning, and actual conditions. Compared with the prior art, it is more comprehensive, clear, and accurate. BRIEF DESCRIPTION OF THE DRAWINGS

[0019] In order to more clearly illustrate the technical solutions of the exemplary embodiments of the present invention, the following briefly introduces the drawings required for use in the embodiments. It should be understood that the following drawings only illustrate certain embodiments of the present invention and should not be regarded as limiting the scope. For those of ordinary skill in the art, other relevant drawings can be obtained based on these drawings without creative work. In the drawings: Figure 1 This is a flow chart of the carbon monitoring device deployment solution in an embodiment of the present invention; Figure 2 The DT decomposition algorithm flow in the embodiment of the present invention; Figure 3 The distribution of industrial park areas and selected points in the embodiment of the present invention; Figure 4 The triangulation and numbering in the embodiment of the present invention; Figure 5 The encrypted triangulation and numbering in the embodiment of the present invention; Figure 6This is a distribution diagram of carbon dioxide monitoring devices in an industrial park in an embodiment of the present invention. DETAILED DESCRIPTION

[0020] In order to make the objectives, technical solutions and advantages of the present invention more clearly understood, the present invention is further described in detail below in conjunction with embodiments and drawings. The exemplary embodiments of the present invention and their description are only used to explain the present invention and are not intended to limit the present invention.

[0021] As a possible implementation, Figure 1 As shown, this embodiment provides a method for monitoring points of carbon dioxide in an industrial park based on industrial park gas monitoring, which specifically includes the following steps: 1. Analyze the characteristics of the industrial park and study the information of the industrial park, specifically: Obtain the scale, industry type, production process, and emission source possibilities of the industrial park. Collect technical parameters and data of relevant industrial processes and equipment, including: The geographical location, regional boundaries, area, traffic conditions, and connection with surrounding infrastructure of the park; the construction, operation and management of the park; the scale, industry type, production process, and emission source possibility of the park. Collect technical parameters and data of relevant industrial processes and equipment. Distribution of emission sources in the park. The park floor plan, marking the basic orientation; the location, name and size of workshops and other major buildings; organized and unorganized emission outlets and their main parameters; types and emission rates of emission sources; the height and nature of the perimeter wall of the park (enclosed or ventilated); the distribution of buildings and terrain in the park that affect airflow movement, and marked on the park floor plan.

[0022] 2. Collect meteorological data of the park and determine the monitoring scenarios, specifically: Meteorological data of the park in the past 10 years were collected and analyzed. Meteorological elements mainly include: wind speed, wind direction, temperature, humidity, precipitation, surface radiation, temperature profile and tropospheric height, etc.; the meteorological conditions on the site of the industrial park were simply measured and determined, and appropriate meteorological conditions should be selected for monitoring points, because the diffusion and dilution of carbon dioxide depend on the movement of the atmosphere due to the influence of meteorological environmental factors. Under different meteorological conditions, the carbon dioxide concentration of the same emission source at the same location may differ by dozens of times. It is very important to select appropriate on-site meteorological conditions. The specific suitability is shown in Table 1. According to the Technical Guidelines for Monitoring Unorganized Emissions of Atmospheric Pollutants (HJ / T55-2000), the more suitable wind direction and wind speed conditions are <29° and 1.0~3.0m / s, respectively, and the atmospheric stability is between DF levels.

[0023] Table 1 Classification of suitability of monitoring various meteorological factors ; Note: A, B, C, D, E, and F represent strong instability, unstable, weak instability, neutral, relatively stable, and stable, respectively; a, b, c, and d represent decreasing suitability for monitoring.

[0024] Usually, a 10,000-cup portable wind direction and anemometer can be used, or other portable wind direction and anemometers with the same function can be used. Prepare the instruments, equipment and reagents specified in the sampling part of the fugitive emission monitoring in the corresponding standard analysis method of the substance to be tested.

[0025] Place the portable wind direction and anemometer in an open area of ​​the unit to be measured. If there is no suitable open area on site, the portable wind direction and anemometer can be placed at a high place (but generally not more than 15m) to measure the wind direction and speed.

[0026] According to the instrument manual, turn on the brake switch of the portable wind direction and wind speed meter and start reading. Read the instant wind direction and wind speed value every 1 minute. Continue measuring for 10 minutes and get a total of 10 wind direction values ​​(accurate to 50) and 10 wind speed values ​​(accurate to 0.1m / s). If the wind direction and wind speed change greatly at that time and you feel that 10 minutes of measurement is still insufficient, you can appropriately extend the measurement time.

[0027] The 10-minute average wind speed is calculated from 10 wind speed readings; the average wind direction and the standard deviation of wind direction change (±S degrees) are calculated from 10 wind direction readings and filled in Table 2.

[0028] Table 2 Simple measurement results of wind direction and wind speed in open areas ; In addition to being measured before sampling, wind direction and wind speed should be repeated 1-2 times during the sampling process. If a significant change in wind direction is found, the monitoring point should be moved and re-sampled.

[0029] According to the characteristics of carbon dioxide emissions in industrial parks, the analysis is conducted, taking into account the fixed emission sources within the park. The monitoring time is selected as the peak production period, that is, when carbon dioxide emissions show a clear increasing trend, the impact of production activities in the park on carbon dioxide concentration can be fully captured. Determine the monitoring time to choose cloudy or cloudy weather, when the atmospheric stability is high and the diffusion and dilution capacity of carbon dioxide is relatively weak. Choose monitoring under conditions that are not conducive to the diffusion and dilution of carbon dioxide, which can more accurately reflect the distribution of carbon dioxide concentration within the park. In order to obtain accurate monitoring data, continuous monitoring is carried out for 3 days under normal operating conditions to ensure that the monitoring data is relatively accurate.

[0030] 3. Compare the environmental monitoring point layout methods and choose the grid method for monitoring point layout. Specifically: On the basis of following the above basic principles of point layout, it is necessary to flexibly select the basic method of atmospheric environment monitoring point layout to improve the effectiveness of atmospheric environment monitoring. Specifically, it mainly includes grid point layout method, functional area point layout method, concentric circle point layout method and fan point layout method.

[0031] (1) Functional area layout method The functional area distribution method divides an area into industrial areas, residential areas, densely trafficked areas, prosperous commercial areas, cultural areas, clean areas, control areas, etc. according to their functions. The functional area distribution method is a method used to plan and arrange functional areas, aiming to reasonably allocate the location and scope of different functional areas to meet specific needs and requirements.

[0032] (2) Geometric point distribution method The geometric figure distribution method includes grid distribution method, concentric circle distribution method and fan distribution method.

[0033] 1) Grid point method The grid point distribution method is to divide the detection area into grids, and place the sampling points at the center of the grid or the intersection of the straight lines according to the detection point requirements. Secondly, when setting the fixed points, the various conditions of the area must be fully considered. The detection personnel need to combine the actual carbon source emissions, population mobility, funds and other aspects. The sampling points are finally determined through a comprehensive analysis of carbon source emissions. The point grid distribution method is more suitable for areas with severe atmospheric environmental pollution or areas with uniform regional pollution distribution, so as to ensure the effectiveness of the grid distribution method, conform to the principle of consistency between the surrounding concentration and the concentration of the distribution point location, and maximize the accuracy of monitoring.

[0034] 2) Concentric circle point distribution method The concentric circle distribution method is a commonly used environmental monitoring sampling method. The main steps are as follows: ① Find the center of the pollution source and draw concentric circles (5-7) with it as the center.

[0035] ② List several radial lines from the center of the circle. The intersection of the radial lines and the circle is the sampling point position (at least five rays, usually 8).

[0036] ③The radii of the concentric circles are 4, 10, 20, and 40 km respectively, and 4 or 8 sampling points are set on each circle (which can be flexibly set depending on the upwind or downwind direction).

[0037] This method of point distribution is mainly used for pollution clusters composed of multiple pollution sources, or areas with concentrated pollution. Since there are many monitoring points arranged according to this method, the monitoring points should be selected according to certain principles: the points should be denser in the expected high-concentration areas and the junction areas of high-concentration and low-concentration, and sparser in other areas; the points should be denser in the downwind direction, and sparser in other directions. This method can more representatively reflect the concentration of pollutants and provide reliable samples for environmental monitoring.

[0038] 3) Fan-shaped point distribution method The fan-shaped point distribution method is a commonly used environmental monitoring sampling method. The main steps are as follows: ①Determine the source of pollution and take it as the apex.

[0039] ② With the dominant wind direction as the axis, draw a fan-shaped area on the ground in the downwind direction as the point distribution range 1. The fan-shaped angle is generally 45° to 90°.

[0040] ③ The sampling points are set on several arcs at different distances from the point source, and the angle between the line connecting two adjacent points and the vertex is generally 10°~20°.

[0041] This method of point distribution is mainly used for isolated elevated point sources and areas with obvious dominant wind direction. This method can more representatively reflect the concentration of pollutants and provide reliable samples for environmental monitoring.

[0042] Due to the large number of industrial parks in my country, their wide variety and different development stages, the current monitoring device deployment method is mostly used for uniform deployment, without considering the actual situation of each industrial park. Carbon dioxide concentration monitoring is carried out based on the coverage area of ​​the industrial park and the strong emission sources, and deployment points are set around.

[0043] Therefore, this embodiment will adopt the grid point distribution method to solve the problems existing in the prior art, divide the industrial park into grids, and select at least one monitoring point in each grid to monitor the carbon dioxide concentration, clarify the functions of each area in a park and the characteristics of the industrial chain structure, and distribute points based on the functional area grid point distribution method.

[0044] 4. Based on the spatial analysis theory, the DT triangulation method is used to determine the monitoring points. Specifically: The most commonly used triangulation in practice is DT, which is a special triangulation. DT means that the triangles of the triangulation contain only Delaunay edges (Delaunay). Delaunay edges must satisfy: there is a circle passing through the two endpoints of this edge, and the circle (note that it is the circle, and at most three points on the circle are co-circular) does not contain any other points in the point set. DT has two characteristics: ① Empty circle characteristic. The Delaunay triangulation is unique (any four points cannot be co-circular), and there will be no other points within the range of the circumscribed circle of any triangle in the Delaunay triangulation. ② Maximize the minimum angle characteristic. Among the triangulations that may be formed by the scattered point set, the minimum angle of the triangle formed by the DT triangulation is the largest. In this sense, the Delaunay triangulation is the "closest to regularization" triangulation. DT is unique, and when adding, deleting, or moving a vertex in an area, it will only affect the adjacent triangles. These excellent characteristics make DT well applied in many fields. The atmospheric environment is an uneven and changing continuum. For each divided triangular grid, selecting the center point of the corner grid as the layout point of the sensor can well reflect the concentration of a certain gas in this area.

[0045] At the same time, the idea of ​​spatial analysis is used to guide the optimal layout design of the monitoring network, so as to obtain representative data of the largest range of space with the least number of monitoring points, accurately and objectively reflect the status and trend of carbon dioxide emissions in industrial parks, and reflect more information. Classical statistical methods are used to analyze the spatial variability of carbon dioxide distribution. The spatial variability of carbon dioxide distribution is described by calculating the maximum value, minimum value, mean, standard deviation, variance, skewness, kurtosis, coefficient of variation, etc. of sample data. In general, C v When the spatial variability of carbon dioxide distribution is less than 10%, it is a weak variability. v Between 10% and 100%, the variation is moderate. v When ≥100%, it is a strong variation degree, and the specific calculation formula is: ; Where: C v is the coefficient of variation, σ is the standard deviation, and μ is the mean.

[0046] An important aspect in the process of collecting and calculating sample data of random variables is the spatial correlation of samples. This spatial continuity is an important characteristic of many natural phenomena. Geostatistics provides a theory to describe this spatial continuity and combines classical regression methods with spatial continuity to develop better methods.

[0047] Let x1,…,x n is a series of observation points on the region, z1,…,z nis the corresponding observation value. And any set of any vector, if z(x i ),…,z(x n ) and the mean of z(x1),…,z(x n +h) have the same mean, and the covariance Cov[Z(x),Z(y)] depends only on the distance between the sample points x and y, that is, |xy|, then we say that Z is second-order stationary, and the result is that the mean and variance of Z are constants that do not change in space.

[0048] If we further relax the requirements for random variables, we say that random functions have inherent assumptions under the following conditions: The existence of the mean does not depend on x, that is, E[Z(x)]=m. For any x, E[Z(x+h)-Z(x)]=0. For any distance h, the variable E[Z(x+h)-Z(x)]=m has a finite variance that does not depend on x.

[0049] For any x and h, the following holds: ; The following formula holds: where γ(h) is called the (semi)variogram (or variogram). The variogram is used to characterize the spatial variation structure or spatial continuity of random variables. It is the basis of geostatistics. It can describe both the spatial structural changes of regionalized variables and their random changes.

[0050] Listed below are several commonly used variogram models: 1) Spherical model: ; Where C0 is the nugget constant, C0+C1 is the base value, and a is the range.

[0051] 2) Exponential model: ; When h=3a, y(h)≈C0+C1, ​​so its range is 3a.

[0052] 3) Gaussian model: ;when When , y(h)≈C0+C1, ​​so its range is .

[0053] 4) Power function model (0<λ<2): ; When λ=1, the power function model is the commonly used linear model: y(h)=C0+C1h; In summary, the spatial variability structure of carbon dioxide characteristics is fitted using the semivariogram function to obtain a functional theoretical model. The calculation formula is as follows: ; Where: γ(h) is the sample variance function of the regionalized variable Z(X); h is the spatial separation distance between two sample points; Z(X i ) and Z(X i +h) are the localized variables Z(X) at the spatial position X i and X i +h is the observed value, N(h) is the number of sample pairs with spacing.

[0054] The variance function has three important parameters. Among them, the ratio of the nugget value to the sill, namely the nugget-sill ratio, is usually used to describe the spatial correlation of the variables. When the ratio is less than 25%, it means that there is a strong spatial correlation between the variables. When the ratio is between 25% and 75%, it means that there is a moderate degree of spatial correlation between the variables. When the value is greater than 75%, it means that there is a weak spatial correlation between the variables.

[0055] like Figure 2 As shown in the figure, it is the flow chart of the algorithm in this paper. The algorithm can realize the sensor layout based on Delaunay triangulation.

[0056] STEP 1: Select a suitable point set based on the given industrial park area and area. The distribution of the point set is generally required to cover the entire area, show the outline of the entire area, and the carbon emission source points in the industrial park, so as to ensure that certain areas are not missed during the triangulation process.

[0057] STEP2: For the selected point set, use the DT triangulation method to triangulate the entire area and generate a triangular mesh.

[0058] STEP 3: For each triangular mesh generated, use a handheld carbon monitoring device to detect the measurement values ​​of sensors at three nodes.

[0059] STEP 4: Based on the measurement value of each triangular mesh node, calculate the corresponding coefficient of variation (reflecting the severity of the fluctuation of the measurement value in this area).

[0060] STEP5: Compare the coefficient of variation with the set threshold. If the coefficient of variation is greater than the threshold, go to STEP 6; otherwise, go to STEP 8.

[0061] STEP6: Determine the size of the circumscribed circle of the triangular mesh and the sensor radius, consider the sensor's measurement range, and prevent the mesh from being too dense, which will lead to unnecessary waste and meaningless measurement. If the circumscribed circle of the triangular mesh is larger than the sensor radius, go to STEP7; otherwise, go to STEP8.

[0062] STEP7: Select appropriate points in the triangular mesh, use the DT triangular mesh division method to encrypt the mesh, and then go to STEP3.

[0063] STEP8: For the triangular mesh that meets the requirements, select its center as the sensor point. Continue in this way until the coefficient of variation of each triangular mesh meets the threshold requirement, and finally find the sensor point corresponding to each triangular mesh.

[0064] For different types of industrial parks, the point distribution process and method are the same, but the threshold values ​​are set differently to achieve different density of points. Carbon dioxide is a gas with a higher density than air. Specifically, the molecular weight of carbon dioxide is about 44, while the average molecular weight of air is 29. Theoretical analysis shows that in the absence of other power, gases heavier than air are more likely to deposit on the ground. Therefore, the height of the monitoring point is set between 1 and 4 meters in order to achieve accurate monitoring of the carbon dioxide concentration in the near-source area on the ground. Since DT encryption will only affect nearby triangles, it is only necessary to calculate the coefficient of variation of several affected triangles in the surrounding area, which can increase efficiency and save time.

[0065] Set the selected industrial park area as Figure 3 As shown, select the points on the contour and the carbon emission source point in the center, and the final effect is as follows Figure 3 The industrial park area and selected point distribution are shown.

[0066] Use matlab to call the Delaunay function to complete the segmentation. Figure 4 As shown in Table 3, it is divided into 13 grids in total and numbered. According to the position of each triangular grid, the carbon dioxide value at each node is measured, and then the coefficient of variation of each triangular grid is calculated, as shown in Table 3.

[0067] Table 3 Coefficient of variation of triangular mesh node measurements ; Assuming that the coefficient of variation is 30% and the nugget base ratio is 25% as the threshold, it can be seen that grids 1, 4, and 12 do not meet the requirements and need to be further encrypted. The encrypted effect is shown in the figure below. Figure 5 After numbering each grid, the coefficient of variation of each grid is calculated according to the value of the node. If the threshold requirement is met, the encryption is stopped.

[0068] like Figure 6 The layout scheme shown is formulated according to the functional area layout method. This scheme is mainly planned according to the actual needs of each functional area in the industrial park and the carbon dioxide emission data, so as to achieve comprehensive and accurate monitoring of the carbon dioxide emissions in the park and provide support for carbon emission monitoring data.

[0069] A second aspect of this embodiment provides an electronic device, including a memory, a processor, and a computer program stored in the memory and executable on the processor. When the processor executes the program, a monitoring point deployment method based on industrial park gas monitoring is implemented.

[0070] A third aspect of this embodiment provides a computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements a monitoring point deployment method based on gas monitoring in an industrial park.

[0071] The specific implementation methods described above further illustrate the objectives, technical solutions and beneficial effects of the present invention in detail. It should be understood that the above description is only a specific implementation method of the present invention and is not intended to limit the scope of protection of the present invention. Any modifications, equivalent substitutions, improvements, etc. made within the spirit and principles of the present invention should be included in the scope of protection of the present invention.

Claims

1. A monitoring point layout method based on industrial park gas monitoring, characterized in that: The specific steps include: Obtain park characteristic data of the target monitoring area, and construct a characteristic data point set based on the characteristic data; Based on the feature data point set, the DT triangulation method is used to triangulate the target monitoring area and generate a triangular mesh; Obtain the measurement values ​​of the sensors at the three nodes of each triangular mesh and calculate the coefficient of variation of each triangular mesh; In the ith triangle mesh: If the coefficient of variation is within the set threshold, the center point of the i-th triangular mesh is selected for point distribution; If the coefficient of variation is not within the set threshold, it is determined whether the circumscribed circle of the i-th triangular mesh is smaller than the sensor monitoring radius; If it is not less than, the Delaumnay triangulation algorithm is used to encrypt the i-th triangular mesh, and the execution is jumped to obtain the measurement values ​​of the sensors on the three nodes of each triangular mesh; If it is less than, the center point of the i-th triangle mesh is selected for point distribution.

2. The monitoring point layout method based on industrial park gas monitoring according to claim 1 is characterized in that: The park characteristic data of the target monitoring area includes: the park characteristics of the target monitoring area, the emission data distribution of the target monitoring gas in the area, and the meteorological data in the area.

3. The monitoring point layout method based on industrial park gas monitoring according to claim 2 is characterized in that: The park characteristics of the target monitoring area include: the park's geographical location, regional boundaries, area, transportation conditions, connection with surrounding infrastructure, park scale, industry type and production process; The distribution of emission data of target monitored gases in the area includes: the location and size of workshops and other major buildings, organized and unorganized emission outlets and their emission parameters, types and emission rates of emission sources, the closedness of the target monitoring area of ​​the park, and the distribution of buildings and terrain in the park that affect airflow movement; The meteorological data in the area include: wind speed, wind direction, temperature, humidity, precipitation, surface radiation, temperature profile and tropospheric altitude.

4. The monitoring point layout method based on industrial park gas monitoring according to claim 3 is characterized in that: When conducting monitoring deployment, the meteorological data in the area shall be measured and determined, including: Collect wind speed and direction, obtain an instant wind direction and wind speed value every t time, perform continuous measurement within T period, obtain n wind direction values ​​and n wind speed values, and measure and judge meteorological data based on n wind direction values ​​and n wind speed values; During the collection of wind speed and direction, if the deviation of n wind direction values ​​is greater than the set threshold, the monitoring point is moved and resampled.

5. The monitoring point layout method based on industrial park gas monitoring according to claim 1 is characterized in that: The DT triangulation method is used to triangulate the target monitoring area, specifically including: Obtain the boundary of the target monitoring area, and select the boundary point and the carbon emission source point in the center of the target monitoring area according to the boundary of the target monitoring area; According to the carbon emission source points at the boundary points and the center of the target monitoring area, the target monitoring area is divided into multiple grids using the Delaunay function.

6. The monitoring point layout method based on industrial park gas monitoring according to claim 1 is characterized in that: The calculation of the coefficient of variation of each triangular mesh specifically includes: , where σ is the standard deviation and μ is the mean.

7. The monitoring point layout method based on industrial park gas monitoring according to claim 1 is characterized in that: Before arranging points on the i-th triangular grid, the spatial variability structure of carbon dioxide characteristics is fitted using a semivariance function to obtain a function model. The specific steps include: ； Where: N(h) is the number of sample logarithms of the interval; γ(h) is the sample variance function of the regionalized variable Z(X); h is the spatial separation distance between two sample points; Z(X i ) and Z(X i +h) are the localized variables Z(X) at the spatial position X i and X i The observed value of +h is the number of sample pairs with spacing.

8. The monitoring point layout method based on industrial park gas monitoring according to claim 1 is characterized in that: After the points are obtained based on the triangular mesh, the points are also optimized according to the actual situation of the park, including: Obtain the monitored gas emission data of the target monitoring area, and determine the first dominant wind direction and the second dominant wind direction of the carbon emission data based on the monitored gas emission source data and meteorological data; Based on the first dominant wind direction and the second dominant wind direction of the emission data, the gas concentration distribution of each grid is obtained; The points are optimized according to the gas concentration distribution of each grid.

9. An electronic device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that: When the processor executes the program, the monitoring point deployment method based on industrial park gas monitoring as described in any one of claims 1 to 8 is implemented.

10. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the program is executed by a processor, the monitoring point layout method based on industrial park gas monitoring as described in any one of claims 1 to 8 is implemented.

Images