Method and device for acquiring settlement data of high-voltage line tower in coal mining subsidence area

Abstract

The invention discloses a method and a device for acquiring settlement data of a high-voltage line tower in a coal mining subsidence area. The method comprises the following steps: acquiring observation data of the settlement of each monitoring point; determining a missing sink value and missing time according to the observation data; when the missing time is larger than or equal to the first time length, pile repairing monitoring is conducted on the missing monitoring points, and the sinking speed value after pile repairing is obtained; the sinking speed values before monitoring point missing and after pile repairing are calculated in a segmented mode; performing overall fitting on the sinking speed distribution curve according to the subsection calculated sinking speed values before the monitoring points are missing and after pile repairing to obtain a sinking speed time domain distribution curve of the missing monitoring points; according to the missing time of the missing monitoring point and the sinking speed time domain distribution curve, the sinking speed value of the missing monitoring point at the missing time is obtained; and calculating a missing subsidence value in the missing time according to the subsidence speed value of the missing monitoring point in the missing time.

Description

TECHNICAL FIELD

[0001] The present application belongs to the technical field of subsidence monitoring in coal mining subsidence area, and particularly relates to a method and device for acquiring subsidence data of high-voltage line tower in coal mining subsidence area. BACKGROUND

[0002] In the coal mining subsidence area, the stability of the high-voltage line tower is crucial for power transmission, and the traditional monitoring methods (total station measurement, satellite positioning technology, etc.) have many limitations. The total station measurement relies on the visibility condition and is difficult to monitor comprehensively in the complex terrain of the coal mining subsidence area. The GNSS positioning accuracy and stability are affected in the areas with serious signal obstruction. In the field of hydrostatic leveling measurement, although there are technologies for subsidence monitoring, the existing devices and methods based on this technology still have deficiencies in monitoring the high-voltage line tower in the coal mining subsidence area. On the one hand, in the face of electromagnetic interference, vibration and complex environment in the coal mining subsidence area, the adaptability of the existing devices is poor, and the accuracy and reliability of the measurement data are difficult to guarantee. On the other hand, the data processing and analysis means are simple, and it is difficult to accurately evaluate the safety state of the high-voltage line tower and timely alarm. SUMMARY

[0003] The purpose of the embodiments of the present application is to provide a method and device for acquiring subsidence data of high-voltage line tower in coal mining subsidence area, to solve the problem of difficult accurate monitoring of subsidence data of high-voltage line tower in coal mining subsidence area.

[0004] In a first aspect, the embodiments of the present application provide a method for acquiring subsidence data of high-voltage line tower in coal mining subsidence area, comprising:

[0005] Obtaining observation data of subsidence amount of each monitoring point;

[0006] Determining the missing subsidence value and the missing time according to the observation data;

[0007] When the missing time is greater than or equal to the first time length, performing pile supplement monitoring on the missing monitoring point to obtain the subsidence velocity value after pile supplement;

[0008] Segmented calculating the subsidence velocity value of the monitoring point before and after pile supplement;

[0009] Overall fitting the subsidence velocity distribution curve according to the segmented calculated subsidence velocity value of the monitoring point before and after pile supplement, to obtain the subsidence velocity time domain distribution curve of the missing monitoring point;

[0010] Obtaining the subsidence velocity value of the missing monitoring point at the missing time according to the missing time and the subsidence velocity time domain distribution curve of the missing monitoring point;

[0011] Calculating the missing subsidence value in the missing time according to the subsidence velocity value of the missing monitoring point at the missing time.

[0012] Optionally, the step of fitting the subsidence velocity distribution curve as a whole according to the subsidence velocity values of the monitoring points before missing and after filling in the segments comprises:

[0013] The subsidence velocity distribution curve is fitted as a whole according to the subsidence velocity values of the monitoring points before missing and after filling in the segments by using a normal density function.

[0014] Optionally, the acquisition method further comprises:

[0015] When the missing time is less than the first time length, the missing subsidence value is acquired according to the observation data of the monitoring points by using a linear interpolation method.

[0016] Optionally, the acquisition method further comprises:

[0017] The mean value of the observation data is acquired according to the observation data of the subsidence of the monitoring points;

[0018] The subsidence value whose deviation from the mean value in the observation data of the subsidence of the monitoring points is greater than or equal to 3 times the standard deviation is determined as a subsidence abnormal value; and / or

[0019] The observation data of the subsidence of the monitoring points is denoised by using an adaptive Kalman filtering method.

[0020] Optionally, the acquisition method further comprises:

[0021] The angle θ of the tilt angle between any two monitoring points is calculated, and the relationship between the angle θ of the tilt angle and the horizontal distance between the two monitoring points is:

[0022]

[0023] Δh is the difference between the subsidence values of any two monitoring points, and L is the horizontal value between any two monitoring points.

[0024] Optionally, the acquisition method further comprises:

[0025] The subsidence value of the monitoring point is predicted by using a subsidence value prediction model, and the calculation formula of the prediction model is:

[0026]

[0027] wherein, is the predicted value of the subsidence value at time t+1; y t is the actual observation value at time t; is the predicted value of the subsidence value at time t.

[0028] Optionally, the acquisition method further comprises:

[0029] According to the observation data of the subsidence amount of the monitoring point and the missing subsidence value, the subsidence value change speed of the monitoring point in a unit time is obtained;

[0030] The reminding information is sent according to the subsidence value change speed of the monitoring point.

[0031] In a second aspect, an embodiment of the present application provides a coal mining subsidence area high-voltage line tower subsidence data acquisition device, comprising:

[0032] A first acquisition module is configured to acquire observation data of the subsidence amount of each monitoring point;

[0033] A determination module is configured to determine the missing subsidence value and the missing time according to the observation data;

[0034] A second acquisition module is configured to acquire the subsidence speed value after the pile supplement of the missing monitoring point when the missing time is greater than or equal to the first time length;

[0035] A first calculation module is configured to calculate the subsidence speed value before and after the pile supplement of the missing monitoring point in sections;

[0036] A fitting module is configured to perform overall fitting on the subsidence speed distribution curve according to the subsidence speed value before and after the pile supplement of the missing monitoring point calculated in sections, to obtain the subsidence speed time domain distribution curve of the missing monitoring point;

[0037] A third acquisition module is configured to obtain the subsidence speed value of the missing monitoring point at the missing time according to the missing time of the missing monitoring point and the subsidence speed time domain distribution curve;

[0038] A second calculation module is configured to calculate the missing subsidence value in the missing time according to the subsidence speed value of the missing monitoring point at the missing time.

[0039] In a third aspect, an embodiment of the present application provides an electronic device, comprising a processor and a memory, wherein the memory stores programs or instructions executable on the processor, and the programs or instructions are executed by the processor to implement the steps of the acquisition method described in the above embodiments.

[0040] In a fourth aspect, an embodiment of the present application provides a readable storage medium, wherein the readable storage medium stores programs or instructions, and the programs or instructions are executed by a processor to implement the steps of the acquisition method described in the above embodiments.

[0041] In the coal mining subsidence area high-voltage line tower settlement data acquisition method, the settlement speed values before and after the missing monitoring point is supplemented are calculated in sections, the settlement speed distribution curve is overall fitted according to the settlement speed values before and after the missing monitoring point is supplemented calculated in sections, the settlement speed time domain distribution curve of the missing monitoring point is obtained, the settlement speed value of the missing monitoring point in the missing time period is obtained according to the missing time of the missing monitoring point and the settlement speed time domain distribution curve, and the missing settlement value in the missing time is calculated according to the settlement speed value of the missing monitoring point in the missing time. Through the method in the application, the missing settlement value of the missing monitoring point can be calculated according to the missing monitoring point data, the missing data is processed in time and space dimensions, compared with the traditional method of estimating the missing data only in the spatial dimension, the application has higher precision, can obtain more comprehensive and accurate settlement information, can accurately evaluate the safety state of the high-voltage line tower, and timely discovers and warns potential safety hazards, and effectively guarantees the safe and stable operation of the power system. BRIEF DESCRIPTION OF DRAWINGS

[0042] Figure 1 It is a connection schematic diagram of the coal mining subsidence area high-voltage line tower settlement data acquisition device.

[0043] Reference signs

[0044] The first acquisition module 11, the second acquisition module 12 and the third acquisition module 13

[0045] The determination module 20

[0046] The first calculation module 31 and the second calculation module 32

[0047] The fitting module 40 DETAILED DESCRIPTION

[0048] The technical solutions in the embodiments of the application will be clearly and completely described below with reference to the drawings in the embodiments of the application. Obviously, the described embodiments are some of the embodiments of the application, rather than all the embodiments of the application. Based on the embodiments in the application, all other embodiments obtained by those skilled in the art without creative labor fall within the protection scope of the application.

[0049] The terms "first", "second" and the like in the specification and claims of the application are used to distinguish similar objects, and are not used to describe a particular order or sequence. It should be understood that the data used in this way can be interchanged under appropriate circumstances, so that the embodiments of the application can be implemented in an order other than those illustrated or described herein. In addition, "and / or" in the specification and claims indicates at least one of the connected objects, and the character " / ", generally indicates that the front and rear associated objects are in an "or" relationship.

[0050] The coal mining subsidence area high-voltage line tower settlement data acquisition method of the embodiment of the application comprises:

[0051] Obtaining observation data of subsidence amounts of each monitoring point;

[0052] Determining missing subsidence values and missing times according to the observation data;

[0053] When the missing time is greater than or equal to a first time length, performing pile supplement monitoring on the missing monitoring point to obtain a subsidence speed value after pile supplement; the missing monitoring point can be supplemented with a pile at the original position, and the pile supplement monitoring is performed on the missing monitoring point to obtain the subsidence speed value after pile supplement;

[0054] Segmentally calculating the subsidence speed values of the missing monitoring point before and after pile supplement;

[0055] Overall fitting a subsidence speed distribution curve according to the segmentally calculated subsidence speed values of the missing monitoring point before and after pile supplement to obtain a subsidence speed time domain distribution curve of the missing monitoring point;

[0056] Obtaining a subsidence speed value of the missing monitoring point at the missing time according to the missing time of the missing monitoring point and the subsidence speed time domain distribution curve;

[0057] Calculating a missing subsidence value in the missing time according to the subsidence speed value of the missing monitoring point at the missing time.

[0058] The pile supplement is to reset a monitoring instrument at the original position of the missing monitoring point or adjust the monitoring instrument to continue to obtain the observation data of the subsidence amount of the missing monitoring point. The monitoring instrument can include a hydrostatic level instrument to continue to obtain the observation data of the subsidence amount of the missing monitoring point.

[0059] In the coal mining subsidence area high-voltage line tower settlement data acquisition method of the application, the subsidence speed time domain distribution curve of the missing monitoring point is obtained by overall fitting the subsidence speed distribution curve according to the segmentally calculated subsidence speed values of the missing monitoring point before and after pile supplement, the subsidence speed value of the missing monitoring point in the missing time period is obtained according to the curve, and the missing subsidence value in the missing time is calculated. Through the method in the application, the subsidence value of the missing monitoring point can be calculated according to the missing data of the monitoring point, and the subsidence information can be obtained more comprehensively and accurately from the time and space dimensions. Compared with the traditional method of estimating the missing data only from the spatial dimension, the accuracy of the application is higher, the settlement information can be obtained more comprehensively and accurately, the safety state of the high-voltage line tower can be accurately evaluated, potential safety hazards can be discovered and warned in time, and the safe and stable operation of the power system is effectively ensured.

[0060] In some embodiments, the step of overall fitting the subsidence speed distribution curve according to the segmentally calculated subsidence speed values of the missing monitoring point before and after pile supplement can comprise:

[0061] Based on the segmented calculation of the sinking velocity values ​​before and after the monitoring point was missing and after the pile was added, the sinking velocity distribution curve was fitted with a normal density function. This allows for a more accurate fitting of the sinking velocity distribution, allowing the velocity distribution during the missing time to be accurately fitted, and the amount of sinking to be calculated based on the sinking velocity.

[0062] In an embodiment of the present invention, the acquisition method may further include:

[0063] When the missing time is less than the first time length, the missing subsidence value is obtained by linear interpolation based on the observed subsidence data of the monitoring point. The first time length can be 0.5 hours, 1 hour, 2 hours or 3 hours, preferably 1 hour. The specific time can be selected according to actual needs.

[0064] In some embodiments of the present invention, the obtaining method may further include:

[0065] Obtain the mean of the observation data based on the observation data of the subsidence at the monitoring point;

[0066] A subsidence outlier is defined as a subsidence value whose deviation from the mean is greater than or equal to three standard deviations. A median filter can be used to process subsidence outliers. This filter takes the data before and after the outlier and calculates the median value, which is then used to replace the outlier.

[0067] In some embodiments, after the step of obtaining the observation data of the subsidence amount of each monitoring point, the method may further include: performing denoising processing on the obtained observation data of the subsidence amount of the monitoring point using an adaptive Kalman filter method.

[0068] After the missing monitoring points are filled in their original locations, the observed data of the subsidence of the monitoring points after the filling can be obtained. The observed data of the subsidence of the monitoring points after the filling can be denoised using the adaptive Kalman filter method. Then, the subsidence velocity values ​​before and after the missing monitoring points and the filling can be calculated in sections based on the denoised observed data. Denoising the observed data of the subsidence of the monitoring points obtained using the adaptive Kalman filter method can make data fluctuations smoother, reduce noise interference, and make the processed data more accurate.

[0069] In an embodiment of the present invention, the acquisition method may further include:

[0070] Calculate the inclination angle θ between any two monitoring points. The relationship between the inclination angle θ and the horizontal distance between the two monitoring points can be:

[0071]

[0072] Ah is the difference between the subsidence values of any two monitoring points, and L is the horizontal value between any two monitoring points.

[0073] By calculating the angle of the inclination angle between any two monitoring points, the subsidence and inclination of the high-voltage line tower in the coal mining subsidence area can be evaluated, so as to take measures to protect the safety of the high-voltage line tower.

[0074] In the embodiment of the present application, the acquisition method can further comprise:

[0075] The subsidence value of the monitoring point is predicted by using the subsidence value prediction model, and the calculation formula of the prediction model is:

[0076]

[0077] wherein, is the predicted value of the subsidence value at time t+1; y t is the actual observation value at time t; is the predicted value of the subsidence value at time t.

[0078] Optionally, the acquisition method can further comprise:

[0079] According to the observation data of the subsidence amount of the monitoring point and the missing subsidence value, the subsidence value change speed of the monitoring point per unit time is obtained;

[0080] The warning information is sent according to the subsidence value change speed of the monitoring point.

[0081] If the subsidence value of the subsidence amount of the plurality of monitoring points continuously rises in a period of time, and the subsidence value rising speed gradually increases, it indicates that the high-voltage line tower is in an accelerated subsidence state; if the subsidence value rising rate gradually decreases, it indicates that the high-voltage line tower is in a decelerated subsidence state; if the subsidence value of the subsidence amount remains basically stable, it indicates that the high-voltage line tower is in a relatively stable state. By comparing the subsidence amount difference of the monitoring points at different positions and combining the tower body parameters, the inclination angle can be estimated. When the inclination angle exceeds the safety threshold (for example, generally ±0.5°), it indicates that the high-voltage line tower may be inclined and there is a risk of collapse. When the monitoring data exceeds the set safety threshold, the system can automatically send a warning information, and the warning information can include at least one of the high-voltage line tower position, number, subsidence amount, subsidence rate and safety state level.

[0082] In the application process, multiple hydrostatic leveling instruments can be reasonably distributed around the foundation periphery of the high-voltage line tower, and at the same time, reference datum points are set away from the influence range of coal mining subsidence and in a geologically stable place, and the same hydrostatic leveling instruments are installed, and the multiple hydrostatic leveling instruments are connected to form a monitoring network through high-pressure-resistant and corrosion-resistant connecting pipes. Each hydrostatic leveling instrument can be equipped with a customized sealed enclosure with a protection level of IP68 to prevent the intrusion of dust, mud and other impurities and ensure the normal operation of the instrument in harsh environments. Each hydrostatic leveling instrument can be connected to a data collector, which has a high-performance microprocessor that can quickly and accurately collect liquid level change data output by the hydrostatic leveling instrument. The data collector can integrate a 4G or LoRa wireless communication module to transmit collected data in real time to a data processing center. To ensure data transmission security, an encrypted transmission protocol can be used to prevent data from being stolen or tampered with. To deal with electromagnetic interference in the coal mining subsidence area, multiple layers of electromagnetic shielding layers can be provided in the device to shield the hydrostatic leveling instrument and the data collector in all directions. A vibration isolation device composed of a rubber vibration isolation pad and a spring vibration isolator can be installed between the high-voltage line tower foundation and the installation location of the hydrostatic leveling instrument, which is optimized in design according to the vibration frequency of the high-voltage line tower and the sensitivity of the hydrostatic leveling instrument, effectively reducing the impact of vibration on measurement results.

[0083] After the data processing center receives the data, it uses the Kalman filter algorithm for noise reduction. The data can be detected and repaired for missing values and outliers. Missing value detection can be achieved by comparing adjacent time periods and adjacent measurement point data. If a measurement point has no data for multiple consecutive sampling periods, it is determined that the measurement point data is missing. Outlier detection can use statistical analysis methods such as the 3σ criterion. If the deviation of a data point from the mean exceeds three standard deviations, the data point is determined to be an outlier. For missing values, if the missing time is short (e.g., less than 1 hour), linear interpolation can be used for repair; if the missing time is long (e.g., more than 1 hour), the subsequent missing data processing method is used for processing. For outliers, median filtering can be used for repair.

[0084] Taking a high-voltage line tower in a mine area in Shanxi as an example, the mine area is located in a complex geological area in Shanxi, and frequent coal mining activities have had a significant impact on the stability of the surrounding high-voltage line towers.

[0085] (1) Device installation

[0086] 1. Hydrostatic leveling instrument installation: Dig an installation pit at the designed location around the foundation periphery of the high-voltage line tower, place the hydrostatic leveling instrument smoothly, ensure that the level error is within the allowable range, connect each hydrostatic leveling instrument with a connecting pipe in turn, and seal the connection to prevent liquid leakage. At the same time, install the same hydrostatic leveling instrument at the reference datum point.

[0087] 2. Data acquisition and transmission module installation: Correctly connect the data acquisition device to the hydrostatic level, set the wireless communication module parameters to match the data processing center, and perform communication tests to ensure normal data transmission. Install the electromagnetic shielding layer and vibration isolation device, and adjust and secure them as required.

[0088] 3. Device Debugging: Inject an appropriate amount of measurement liquid into the connecting tube to ensure there are no bubbles and the liquid can flow freely. Initialize the device settings, including parameters such as the data logger's sampling frequency and data transmission interval. Test the device's measurement accuracy and reliability through simulated sedimentation tests, and adjust parameters to ensure normal operation.

[0089] (2) Data acquisition and analysis process

[0090] 1. Data Acquisition: The data collector regularly collects the liquid level change data from the hydrostatic level at a set sampling frequency (e.g., every 10 minutes). After collection, the data is packaged, encrypted, and sent to the data processing center via the wireless communication module. The data collected from the monitoring points is shown in Table 1.

[0091] Table 1 Data of monitoring points

[0092]

[0093] 2. Data preprocessing: After receiving the data, the data processing center uses the Kalman filter algorithm to remove noise. For monitoring point 1, the Kalman filter parameters are set to the noise covariance matrix Q = 0.01 and the measurement noise covariance matrix R = 0.1. Taking 8:00 as an example, the initial state estimate is set to 98.43 (i.e., the measured value), and the state covariance is set to 1. The Kalman filter iterative formula is used to predict the state at the next moment and then update it based on the measured value. For example, at 8:10, the predicted value is combined with the measured value of 98.58, and the filtered value is calculated to be approximately 98.56. Subsequent data from monitoring point 1 is filtered in this way to smooth data fluctuations and reduce noise interference.

[0094] There is no data at monitoring point 2 at 8:40 and 8:50. By comparing the adjacent time periods (the liquid level at 8:30 is 101.17 and the liquid level at 9:00 is 101.28) and the data at adjacent measuring points, it is determined that the data is missing. Since the missing time is less than 1 hour, linear interpolation is used to repair it. 8:40 Liquid Level Calculation

[0095] for: The liquid level at 8:50 is calculated as:

[0096]

[0097] For monitoring point 4, first calculate the mean and standard deviation. Calculate the mean of 13 data points from 8:10 to 10:00 (excluding the outlier value of 122.00 at 8:00) Calculate standard deviation 8:00 data 122.00 mm deviates from the mean by 122.00-100.75=21.25 mm, which exceeds 3 times the standard deviation 3x0.22=0.66 mm, and is determined to be an abnormal value.

[0098] Using the median filtering method, select the median of 100.52 mm for the two data points before and after 8:00, and repair 8:00 to 100.52 mm.

[0099] 3. Missing data processing:

[0100] For monitoring point 2, the subsidence velocity is calculated:

[0101] Before missing, calculate the subsidence velocity from 8:00 to 8:30, the liquid level at 8:00 is 101.08 mm, the liquid level at 8:30 is 101.17 mm, ΔS=101.17-101.08=0.09 mm, Δt=0.5 h, the subsidence velocity

[0102] After filling, calculate the subsidence velocity from 9:00 to 9:30, the liquid level at 9:00 is 101.28 mm, the liquid level at 9:30 is 101.44 mm. ΔS=101.44-101.28=0.16 mm, Δt=0.5 h, the subsidence velocity Using the subsidence velocity values v1=0.18 mm / h and v2=0.32 mm / h, the normal density function is fitted. The parameters can be solved by least squares method, etc., to obtain the mean μ=0.25 mm / h and the standard deviation σ=0.05 mm / h.

[0103] 4. Settlement calculation: The liquid level height H0=100.00 mm at the reference point of the hydrostatic level gauge, and the preprocessed data is shown in Table 2.

[0104] Table 2. Processed data

[0105]

[0106] Relative settlement calculation: Taking monitoring point 1 at 10:00 as an example, the liquid level height H1=99.14 mm, then the settlement S1=H0-H1=100-99.14=0.86 mm, the liquid level height of monitoring point 3 at 10:00 is 100.40 mm, and the settlement S3=H0-H3=100.00-100.40=-0.4 mm. The settlement of each monitoring point at different times is calculated according to this formula.

[0107] ​Tilt calculation: Select monitoring points A (monitoring point 1) and B (monitoring point 3) on the high-voltage tower foundation, which are L = 3m apart. At a certain moment, the settlement of monitoring point 1 is 0.86mm, and the settlement of monitoring point 3 is -0.4mm. The difference in settlement between the two points is ΔS = S1-S3 = 0.86-(-0.4) = 1.26mm. Radians, which are converted to degrees

[0108] 5. Settlement trend analysis: Continuously collect settlement data at each monitoring point at different times to ensure accurate analysis of settlement trends. Use simple exponential smoothing to predict settlement trends. The model formula for the simple exponential smoothing method is: F t+1 =αY t +(1-α)F t , where F t+1 is the predicted value at time t+1, Y t is the actual observation value at time t, F t is the predicted value at time t, and the smoothing coefficient is α.

[0109] Taking monitoring point 1 as an example, by adjusting the smoothing coefficient to 0.7 and performing iterative calculations, the predicted liquid level height values ​​at monitoring point 1 at different times can be obtained. When t is 10:00, the actual value is 99.13mm and the predicted value is 99.10mm; when t is 10:10, F 10:10 =0.7Y 10：00 +0.3F 10:00 =0.7×99.13+0.3×99.10=99.12mm. This prediction result can help determine the subsequent state change trend of the high-voltage tower.

[0110] 6. Safety Assessment and Early Warning: The status of high-voltage towers is determined based on settlement volume, settlement rate, settlement trend, and missing data processing results. When data exceeds the safety threshold, the system automatically issues an early warning to notify relevant departments for action.

[0111] In the curve fitting process, in addition to the normal density function, other bilaterally symmetrical distribution morphology functions, such as Gaussian mixture models and Bessel functions, can be used. According to the actual data characteristics and monitoring requirements, the most suitable function is selected to fit the sinking velocity distribution curve to improve the accuracy and adaptability of the curve fitting. In addition to 4G and LoRa, the wireless communication technology of the data acquisition and transmission module can select other wireless communication technologies such as 5G and NB-IoT according to the network coverage and communication needs of the coal mining subsidence area to ensure the stability and efficiency of data transmission. In terms of anti-interference measures, in addition to the use of electromagnetic shielding layers, software algorithms can be combined for anti-interference processing, such as digital filtering algorithms and adaptive anti-interference algorithms, to further improve the device's anti-interference ability in complex electromagnetic environments.

[0112] During the data acquisition and analysis process, the data collector collects the liquid level change data of the hydrostatic level gauge at a set sampling frequency (e.g., once every 10 minutes) and sends it to the data processing center through the wireless communication module after packaging and encryption. The data collected at the monitoring points in a certain period of time is shown in Table 3.

[0113] Table 3 Data collected at the monitoring points

[0114]

[0115] (1) Abnormal value repair verification: measurement point 1 (8:10) original data 122.00 mm (abnormal value)

[0116] Repair method:

[0117] Calculate the mean μ = 100.74 mm and the standard deviation σ = 0.18 mm of other time points (8:00 to 10:00), 3σ = 100.74 ± 0.54 mm, and replace the original value outside the range with the median.

[0118] Median = median(100.45, 100.59) = 100.52 mm, which meets the expectation.

[0119] (2) Missing value interpolation:

[0120] Calculate the liquid level at 8:40:

[0121] Calculate the liquid level at 8:50:

[0122] (3) Adaptive Kalman filter (AKF):

[0123] Initialization parameters:

[0124] Initial state:

[0125] Initial covariance matrix: P0 = diag(0.1, 0.01)

[0126] Process noise covariance: Q0 = diag(0.01, 0.001)

[0127] Observation noise covariance: R0 = 0.1

[0128] Forgetting factor: b = 0.95

[0129] Filter iteration:

[0130] Prediction:

[0131] Update:

[0132] Adaptive adjustment:

[0133] Filtering result example (part):

[0134]

[0135] (4) Settling amount calculation: reference the liquid level height H0 of the hydrostatic level gauge at the reference point = 100.00 mm, after pretreatment, the relative settling amount calculation: take monitoring point 1 at 10:00 as an example, the liquid level height H1 = 101.02 mm, then the settling amount S1 = H0-H1 = 100-101.02 = -1.02 mm, the liquid level height of monitoring point 3 at 10:00 is 100.40 mm, the settling amount S3 = H0-H3 = 100.00-100.40 = -0.4 mm, the settling amount of each monitoring point at different time is calculated according to the formula.

[0136] Inclination calculation: select monitoring points A (monitoring point 1) and B (monitoring point 3) with a distance of L = 3 m on the high-voltage line tower foundation. At 10:00, the settling amount of monitoring point 1 is -1.02 mm, the settling amount of monitoring point 3 is -0.4 mm,

[0137] (5) Settling trend analysis: through accurate processing of liquid level height data at different times, the settling amount of each monitoring point is calculated by using the formula, and the inclination angle is estimated according to the difference value of the settling amount of multiple monitoring points and the tower body parameters. Then, the settling trend analysis is carried out, the settling amount data of each monitoring point at different times is sorted, and the settling amount-time curve is drawn with time as the horizontal coordinate and the settling amount as the vertical coordinate. The settling trend is judged by carefully observing the curve trend. If the settling amount continuously rises and the rising rate gradually increases in a period of time, the high-voltage line tower is in an accelerated settling state; if the rising rate gradually decreases, it is in a decelerated settling state; if the curve is relatively flat, it means that the high-voltage line tower is in a relatively stable state.

[0138] As Figure 1 shown, the coal mining subsidence area high-voltage line tower settlement data acquisition device of the embodiment of the application comprises:

[0139] The first acquisition module 11 is used for acquiring the observation data of the subsidence amount of each monitoring point;

[0140] The determination module 20 is used for determining the missing subsidence value and the missing time according to the observation data;

[0141] The second acquisition module 12 is used for acquiring the subsidence speed value after pile supplementing of the missing monitoring point when the missing time is greater than or equal to the first time length;

[0142] The first calculation module 31 is used for calculating the subsidence velocity values before and after the missing monitoring point is supplemented in segments;

[0143] The fitting module 40 is used for fitting the subsidence velocity distribution curve as a whole according to the subsidence velocity values before and after the missing monitoring point is supplemented in segments, so as to obtain the time domain distribution curve of the subsidence velocity of the missing monitoring point.

[0144] The third acquisition module 13 is used for obtaining the subsidence velocity value of the missing monitoring point at the missing time according to the missing time of the missing monitoring point and the time domain distribution curve of the subsidence velocity.

[0145] The second calculation module 32 is used for calculating the missing subsidence value in the missing time according to the subsidence velocity value of the missing monitoring point at the missing time.

[0146] The electronic device of the embodiment of the present application comprises a processor and a memory, the memory stores a program or instructions which can be run on the processor, and the program or instructions are executed by the processor to realize the steps of the acquisition method described in the above embodiment.

[0147] The readable storage medium of the embodiment of the present application stores a program or instructions, and the program or instructions are executed by the processor to realize the steps of the acquisition method described in the above embodiment.

[0148] The embodiments of the present application are described above in combination with the drawings, but the present application is not limited to the above specific embodiments, and the above specific embodiments are only illustrative but not restrictive, and those skilled in the art can make many forms under the inspiration of the present application without departing from the purpose of the present application and the scope protected by the claims, which all belong to the protection of the present application.

Claims

1. A method for acquiring settlement data of high-voltage line towers in coal mining subsidence areas, characterized in that: include: Obtain observation data of subsidence at each monitoring point; Determine the missing sinking values ​​and missing times based on the observed data; When the missing time is greater than or equal to the first time length, performing pile-filling monitoring on the missing monitoring point to obtain the sinking speed value after the pile-filling; Calculate the sinking velocity values ​​before the monitoring point is lost and after the pile is added in sections; The sinking velocity distribution curve is fitted as a whole based on the sinking velocity values ​​calculated in sections before and after the monitoring point is missing, and the sinking velocity time domain distribution curve of the missing monitoring point is obtained; According to the missing time and sinking speed time domain distribution curve of the missing monitoring point, the sinking speed value of the missing monitoring point at the missing time is obtained; The missing subsidence value during the missing time is calculated based on the subsidence velocity value of the missing monitoring point during the missing time.

2. The acquisition method according to claim 1, characterized in that The steps of overall fitting the sinking velocity distribution curve based on the segmented calculated sinking velocity values ​​before the monitoring point is lost and after the pile is added include: The normal density function is used to fit the sinking velocity distribution curve as a whole based on the segmented calculated sinking velocity values ​​before the monitoring point is lost and after the pile is added.

3. The acquisition method according to claim 1, characterized in that Also includes: When the missing time is less than the first time length, the missing subsidence value is obtained by linear interpolation based on the observation data of the subsidence amount of the monitoring point.

4. The acquisition method according to claim 1, characterized in that Also includes: Obtain the mean of the observation data based on the observation data of the subsidence at the monitoring point; Determine the subsidence values ​​in the observation data of the monitoring point whose deviation from the mean is greater than or equal to 3 times the standard deviation as subsidence outliers; and / or The adaptive Kalman filter method is used to denoise the observation data of the subsidence amount obtained from the monitoring points.

5. The acquisition method according to claim 1, characterized in that Also includes: Calculate the tilt angle θ between any two monitoring points; The relationship between the tilt angle θ and the horizontal distance between two monitoring points is: Δh is the difference in subsidence between any two monitoring points, and L is the level between any two monitoring points.

6. The acquisition method according to claim 1, characterized in that Also includes: The subsidence value prediction model is used to predict the subsidence value of the monitoring point. The calculation formula of the prediction model is: in, is the predicted value of the sinking value at time t+1; t is the actual observation value at time t; is the predicted value of the sinking value at time t, where the value range of α is (0, 1).

7. The acquisition method according to claim 1, characterized in that: Also includes: According to the observation data of the subsidence of the monitoring point and the missing subsidence value, the subsidence value change rate of the monitoring point per unit time is obtained; A reminder message is issued based on the changing speed of the subsidence value of the monitoring point.

8. A device for acquiring settlement data of high-voltage line towers in coal mining subsidence areas, characterized in that: include: The first acquisition module is used to obtain the observation data of the subsidence amount of each monitoring point; A determination module is used to determine the missing sinking value and missing time based on the observation data; The second acquisition module is used to obtain the sinking speed value of the missing monitoring point after the pile is added when the missing time is greater than or equal to the first time length; The first calculation module is used to calculate the sinking velocity value before the monitoring point is lost and after the pile is added in sections; A fitting module is used to perform an overall fitting of the sinking velocity distribution curve based on the segmented calculated sinking velocity values ​​before and after the monitoring point is missing, to obtain the sinking velocity time domain distribution curve of the missing monitoring point; The third acquisition module is used to obtain the sinking velocity value of the missing monitoring point at the missing time according to the missing time and sinking velocity time domain distribution curve of the missing monitoring point; The second calculation module is used to calculate the missing sinking value during the missing time according to the sinking speed value of the missing monitoring point during the missing time.

9. An electronic device, characterized in that: The method comprises a processor and a memory, wherein the memory stores a program or instruction that can be run on the processor, and when the program or instruction is executed by the processor, the steps of the acquisition method according to any one of claims 1 to 7 are implemented.

10. A readable storage medium, characterized in that: The readable storage medium stores a program or instruction, and when the program or instruction is executed by the processor, the steps of the acquisition method according to any one of claims 1 to 7 are implemented.

Images